Next Article in Journal
Universal Adversarial Training Using Auxiliary Conditional Generative Model-Based Adversarial Attack Generation
Next Article in Special Issue
Special Issue on Technological Advances in Seismic Data Processing and Imaging
Previous Article in Journal
Design of a Miniaturized Electron Cyclotron Resonance Ion Source for High-Voltage Proton Accelerator
Previous Article in Special Issue
Three-Dimensional Processing of Reflections for Passive-Source Seismology Based on Geometric Design
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Application of Iterative Virtual Events Internal Multiple Suppression Technique: A Case of Southwest Depression Area of Tarim, China

1
Northwest Branch, Research Institute of Petroleum Exploration & Development of Petro China, Lanzhou 730020, China
2
School of Earth Sciences, Northeast Petroleum University, Daqing 163318, China
3
National Engineering Research Center of Offshore Oil and Gas Exploration, Beijing 100027, China
4
Heilongjiang Provincial Key Laboratory of Oil and Gas Geophysical Exploration, Daqing 163318, China
*
Author to whom correspondence should be addressed.
Appl. Sci. 2023, 13(15), 8832; https://doi.org/10.3390/app13158832
Submission received: 11 March 2023 / Revised: 22 May 2023 / Accepted: 31 May 2023 / Published: 31 July 2023
(This article belongs to the Special Issue Technological Advances in Seismic Data Processing and Imaging)

Abstract

:
The seismic records in the Cambrian southwest depression of the Tarim Basin exhibit discrepancies when compared to the actual geological setting, which is caused by the presence of multiples. Despite the application of the Radon transform, multiple interferences persist beneath the Cambrian salt in the pre-stack data, with significant variations in energy and frequency across the horizontal direction. In addition, other multiple suppression methods are also difficult to handle this problem. To address this issue, we have developed an iterative virtual event internal multiple suppression method for post-stack data. This novel algorithm extends the traditional virtual event internal multiple suppression approach, eliminating the need for data regularization and avoiding the problem of the traditional virtual events method requiring sequential extraction of primaries from relevant layers, which greatly improves computational efficiency and simplifies the implementation steps of the traditional method. Numerical experiments demonstrate the efficacy of our method in suppressing internal multiples in both synthetic and field data while preserving primary signals. When applied to real seismic data profiles, the iterative method yields structural characteristics that align more closely with sedimentary laws and reduces disparities in energy and frequency of multiples along the horizontal axis. Consequently, our method provides a robust foundation for subsequent hydrocarbon source rock prediction.

1. Introduction

The increasing complexity of geological structures necessitates advancements in oil and gas exploration and development technologies, making multiple suppression an essential and challenging problem to address. The presence of multiples interferes with the identification of primary reflections, reducing the signal-to-noise ratio of seismic data, compromising the accuracy and reliability of seismic imaging, and potentially leading to erroneous interpretations and source rock investigations. Consequently, effective suppression of multiples is crucial for enhancing seismic data quality and facilitating oil and gas exploration and development. In northwest China, there are abundant resources of oil and natural gas. However, the complexity of its geological features and strong reflective interfaces in underground media, such as rock mounds, coal seams, basalt, etc., can generate strong internal multiples. The understanding of internal geological structures is constrained. Therefore, the demand for effective suppression of seismic multiples is very urgent.
Seismic multiples are generally classified into two categories based on the location of reflective layers: surface-related multiples and internal multiples. Surface-related multiples exhibit relatively stronger energy, particularly in offshore seismic data. Consequently, numerous researchers have focused on predicting and suppressing surface-related multiples, resulting in the maturation of suppression techniques for this category. In 1992, Verschuur et al. employed the concept of series expansion to decompose the seismic wavefield into primaries and multiples of all orders, effectively suppressing surface-related multiples [1]. This method obviates the need for subsurface geological model information, is a data-driven method, delivers superior results in field data processing, and exhibits high computational efficiency. Many geophysicists have expressed interest in and conducted research on this approach, leading to its widespread adoption in the exploration and development field. In 2009, Van Groenestijn and Verschuur proposed using sparse inversion to estimate primaries, eliminating the requirement for adaptive matched subtraction and near offset extrapolation [2,3]. This advancement significantly optimized the traditional surface-related multiple elimination (SRME) methods. In 2015, Ma et al. developed a suppression method for 3D surface-related multiples [4]. This technique assumes that the multiple contribution trace set corresponding to each seismic record in 3D seismic data is hyperbolic, enabling multiple predictions and suppression in 3D seismic data through sparse inversion. In 2022, Zhu et al. introduced a least squares datum continuation method for suppressing surface-related multiples within the least squares inversion theoretical framework [5]. This approach targets seafloor-separated up- and down-wave records, and, without multiple separation or fractional extraction, iterative inversion yields pre-stack seismic data from virtual observations on the seafloor, enhancing the resolution of effective signals and improving seismic imaging quality.
In comparison to surface-related multiples, internal multiples exhibit smaller differences in frequency, stacking velocity, and normal moveout (NMO) correction, making their suppression more challenging. In 1998, Jakubowicz introduced a data-driven method for internal multiple suppression [6]. This approach directly employs surface observed seismic data to construct internal multiples, enhancing computational efficiency. However, the method necessitates the accurate selection of primary reflection events from seismic data, which is difficult to implement in complex field data. In 2005, Berkhout and Verschuur proposed an internal multiple prediction method based on common focus technology, which is grounded in wavefield extrapolation [7,8]. This method is better suited to complex geological conditions but relies on a velocity model, limiting its application in field data. In 2006, Weglein et al. presented the inverse scattering series method based on the point scattering model [9]. As a data-driven, wave equation-based method, it does not require prior information such as velocity models. However, the method involves substantial computation and is only effective for near offsets, making it insufficient for current data processing demands. In 2013, Ypma and Verschuur introduced an internal multiple suppression method based on sparse inversion [10]. This approach minimizes the damage caused by adaptive matched subtraction to primary signals and effectively protects effective waves. Nevertheless, the method still entails considerable computation and yields unstable wavelets, complicating large-scale field data processing. In 2015, inspired by Marchenko’s imaging method, Meles et al. proposed a self-focusing-based Marchenko internal multiple suppression method [11]. This technique does not require accurate velocity models, only an inaccurate macro velocity model to forward direct waves. Subsequently, direct waves and original data are used to construct the up- and down-going Green functions for relevant virtual source points, generating the associated internal multiples. However, this method demands significant computation due to the need for virtual source points at various depths, making it difficult to apply widely to large scale, complex structured seismic data. In 2018, Zhang and Staring proposed a one-step internal multiple suppression method based on Marchenko’s self-focusing approach, also known as the data domain internal multiple suppression method [12]. This technique does not require macro velocity models for direct wave estimation and directly outputs primaries without internal multiples. However, this method imposes high requirements on the observation system, and its large-scale application in field data remains limited. In 2020, Zhang et al. applied the Marchenko internal multiple suppression method to field seismic data processing, achieving favorable results [13]. Nonetheless, in low signal-to-noise ratio scenarios, energy may not be sufficiently focused, rendering the Marchenko internal multiple suppression method ineffective. In 2022, Zhang et al. compared the self-focusing-based Marchenko method for internal multiple suppression with the data domain method [14]. In the same year, Zhang et al. introduced a multiple suppression method based on self-attention convolutional auto encoders [15]. This method can reduce artificial cost, reduce dependence on unknown prior information, and improve data processing efficiency.
In 2006 and 2009, Ikelle proposed a virtual seismic event construction method to predict internal multiples [16,17]. This approach does not require prior information, such as velocity models or subsurface structures, enabling accurate internal multiple predictions. However, this method requires sequential extraction of primaries to construct internal multiples generated by relevant layers, which increase the computational complexity. In 2013, Wu et al. utilized a multi-channel L1 norm adaptive matching algorithm to suppress multiples, considering the differences in amplitude and phase between virtual event method-predicted internal multiples and actual internal multiples, achieving favorable results in model data [18]. In 2018, Liu et al. introduced an adaptive virtual event method to address the excessive dependency on the matching algorithm in the traditional virtual event method, enhancing internal multiple suppression capabilities [19]. Although researchers have improved and modified the virtual event method to address various issues, the traditional virtual event method still requires sequential internal multiple constructions on subsurface interfaces, incurring high computational costs for seismic data derived from complex structures and making primary extraction difficult. To overcome this challenge, we propose an iterative virtual event method. Compared to the traditional virtual event method, the iterative virtual event technique incorporates an iterative calculation process, mitigates the influence of false frequencies generated by data space convolution on predicted multiple models, and significantly improves computational efficiency. When applied to actual onshore post-stack seismic data from the southwest depression of the Tarim Basin, our method effectively suppresses internal multiples that are challenging to eliminate in pre-stack traces, emphasizes the energy of deep effective waves, and further enhances the quality of seismic profiles.
The paper is organized as follows. After the introduction, we provide an overview of the study region. Subsequently, we present the theory of virtual events for internal multiple eliminations, followed by synthetic and field data experiments that validate the effectiveness of our approach. Finally, we draw conclusions based on detailed analyses and discussions.

2. The Overview of the Study Region

The Southwest Depression of the Tarim Basin, situated in the southwestern part of the basin, is bordered by the South Tianshan orogenic belt to the north and the West Kunlun orogenic belt to the south. Encompassing an area of approximately 100 × 100 km2, it is the largest sub-basin within the Tarim Basin and a favorable area for oil and gas exploration in Cambrian Ordovician carbonate rocks. The traps formed during the Himalayan movement predominantly consist of effective hydrocarbon source rocks distributed in the Bachu fault uplift and Maigaiti slope area (Qu et al., 2000) [20]. After extensive research, a consensus has emerged concerning the Cambrian subsalt of the Maigaiti slope: the Maigaiti slope contains Ordovician Carboniferous oil and gas sources, while drilling through the Cambrian south of Bachu revealed no Cambrian source rock, indicating the development of Cambrian subsalt source rock in the slope area. The Middle-Lower Cambrian serves as an important reservoir and regional cap rock for the area’s deep carbonate rocks and is a significant horizon for hydrocarbon source rock development (Cui et al., 2017; Zhang et al., 2015) [21,22]. Its seismic sequence features three strong reflective interfaces within the Middle Cambrian. The energy at the bottom boundary of the Lower Cambrian is not strong, whereas the bottom of the upper salt section displays a continuous and stable strong wave crest. A noticeable phase change occurs at the bottom of the lower salt section, and the sedimentary pattern from the Early Cambrian to Middle Cambrian exhibits significant changes.
The three strong reflective interfaces in the Middle Cambrian provide favorable conditions for internal multiple formations. However, internal multiples can distort seismic waveforms, decrease the signal-to-noise ratio, and affect effective wave identification, thereby complicating seismic processing, diminishing seismic imaging authenticity and reliability, influencing hydrocarbon source rock studies, and increasing the difficulty of oil and gas exploration and development in the northwest region. To better understand and identify the internal multiples in the study area, well-seismic calibration was performed using existing seismic and logging information, as illustrated in Figure 1. The well-seismic calibration results reveal that the bottom interface of the upper salt section at the Middle Cambrian is stable. The lithology of the lower salt section of the Middle Cambrian comprises pure gypsum salt rock with a strong wave crest at the bottom, while the lower salt section contains gypsum-bearing dolomite and exhibits a trough at the bottom boundary. The bottom boundary of the Lower Cambrian primarily appears as a trough, and the time thickness of the Lower Cambrian varies significantly across different drilling wells, as does the mismatch between well synthetic records and seismic records. Subsequently, velocity spectra were generated for multiple identification and confirmation, as depicted in Figure 2. Figure 2 displays the velocity spectra of two distinct CMP gathers, revealing numerous low-velocity energy clusters. Following NMO correction of CMP gathers using primary velocity, downward-bending events indicate the presence of a large number of multiples in the original seismic data, necessitating multiple suppression processing.
Considering the characteristics of multiples in the study area, we initially employed the Radon transform to suppress internal multiples in the pre-stack data. However, the Radon transform cannot entirely eliminate all internal multiples present in the data, and this method can only suppress some long-period multiples. Multiples interference still exists under the Cambrian salt, as shown by the red circle in the velocity spectrum in Figure 3. Additionally, there are significant differences in energy and frequency of multiples in the horizontal direction. It is challenging to suppress short-path internal multiples in the pre-stack data, necessitating the selection of an appropriate multiple suppression method from the post-stack data for this purpose. The iterative method for suppressing internal multiples based on virtual events is capable of accurately predicting internal multiples without requiring prior information, such as velocity models and underground structures, while maintaining high computational efficiency. Consequently, we opt for the iterative method to suppress internal multiples in the seismic data for the study area.

3. Internal Multiple Elimination Using the Theory of Virtual Events

The iterative virtual event internal multiple suppression approach extends the traditional virtual event internal multiple suppression method and avoids constructing internal multiple from every layer, which greatly improves computational efficiency [23,24]. Seismic virtual events are not directly recorded in standard seismic data acquisition, but their existence can help us construct internal multiples with scattering points at the sea surface [16,17]. Therefore, the total seismic wave field data P ( x R , x S , ω ) received from the surface can be used to successfully construct internal multiples. The process of predicting internal multiples mainly includes two parts. For the internal multiples related to the first underground interface, it is necessary to first extract the primary event P i ( x R , x S , ω ) of ith interface and the primary P i ( x R , x S , ω ) generated by the interface below this interface. Then, correlating the primary P i ( x R , x S , ω ) generated by other interfaces below this interface with the time reversal P i H ( x R , x S , ω ) of the primary event P i ( x R , x S , ω ) to build the virtual event P i v ( x R , x S , ω ) of this interface (Figure 4a). Finally, the convolution operation is carried out by using the constructed virtual event P i v ( x R , x S , ω ) and the primary P i ( x R , x S , ω ) generated by the interface below this interface to construct the internal multiples about this interface (Figure 4b). Namely,
P i v ( x R , x S , ω ) = x P i ( x , x S , ω ) P i H ( x R , x , ω )
M i ( x R , x S , ω ) = x P i v ( x R , x , ω ) P i ( x , x S , ω )
where x represents any point on the surface, R is the receiver point, S is the source, H represents complex conjugation, and ω represents angular frequency. The subscripts represent the number of layers, namely, the lower subscript i in Formulas (1) and (2) represents the ith layer. P i v ( x R , x S , ω ) represents the constructed post-stack virtual events. M i ( x R , x S , ω ) is the internal multiple of the constructed post-stack data.
Figure 4 shows the construction process of virtual events and related internal multiples; constructing post-stack virtual events P v ( x R , x S , ω ) for the correlation of delay wave P ( x , x S , ω ) and lead wave P H ( x R , x , ω ) . Then, the constructed virtual events and primary convolution is used to obtain the internal multiples of the post-stack-related layers M ( x R , x S , ω ) , where the red star represents the convolution operation.
Assuming that the median of seismic data does not include surface-related multiples but only internal multiples, the effective wave after internal multiples suppression is as follows:
P 0 ( x R , x S , ω ) = P ( x R , x S , ω ) M ( x R , x S , ω ) ,
where P 0 ( x R , x S , ω ) is the data of internal multiple suppression and P ( x R , x S , ω ) is the original seismic data. According to the estimation of multiple scattering by iterative inversion [25,26], the multiple can be written as follows:
M ( x R , x S , ω ) = P 0 ( x R , x S , ω ) } 0 A P ( x S , x R , ω )
where A is the surface factor. Substituting Formula (4) into Formula (3):
P 0 ( x R , x S , ω ) = P ( x R , x S , ω ) P 0 ( x R , x S , ω ) } 0 A P 0 ( x S , x R , ω )
In terms of the Neumann series, Formula (5) can be written as follows:
{ P 0 ( x R , x S , ω ) } ( n ) = P ( x R , x S , ω ) { P 0 ( x R , x S , ω ) } 0 ( n 1 ) A ( n ) P ( x R , x S , ω ) ,
where n represents the number of iterations and the expression of surface factor A is as follows:
A ( n ) = ( { P 0 ( x R , x S , ω ) } 0 ( n 1 ) ) 1 ( P ( x R , x S , ω )     { P 0 ( x R , x S , ω ) } n ) ( P ( x R , x S , ω ) ) 1
In order to obtain more abundant internal multiple information, the multiple iteration theory can be used for multiple models. The multiples after iteration n + 1th can be expressed as follows:
{ M ( x R , x S , ω ) } ( n + 1 ) = { P 0 ( x R , x S , ω ) } 0 ( n ) A ( n ) P ( x R , x S , ω )
The multiples after n iterations are expressed as follows:
{ M ( x R , x S , ω ) } ( n ) = { P 0 ( x R , x S , ω ) } 0 ( n 1 ) [ { P 0 ( x R , x S , ω ) } 0 ( n 2 ) ] 1 × ( P ( x R , x S , ω ) { P 0 ( x R , x S , ω ) } 0 ( n 1 ) )
Bringing Formula (7) into Formula (9), the internal multiple model after n iterations is as follows:
{ M ( x R , x S , ω ) } ( n ) = { P 0 ( x R , x S , ω ) } 0 ( n 1 ) ( { P 0 ( x R , x S , ω ) } 0 ( n 1 ) ) T × [ { P 0 ( x R , x S , ω ) } 0 ( n 2 ) ( { P 0 ( x R , x S , ω ) } 0 ( n 2 ) ) T ] 1 × ( P ( x R , x S , ω ) { P 0 ( x R , x S , ω ) } 0 ( n 1 ) )
where T represents matrix transpose.
The primary focus of this article is on the suppression of internal multiples in post-stack profiles. In the post-stack profile, firstly, we can employ the tracing of horizons method to obtain the primary events. Secondly, by shifting the window up and down, we can acquire the primary data volume and other primary data volume below this primary used to construct the virtual events. Subsequently, we input the seismic data obtained, the iterative virtual event internal multiple suppression technique initially utilizes the virtual events internal multiple methods, such as Formulas (1) and (2), to construct the internal multiple models for the corresponding layer. Finally, through iterative updates, such as Formula (10), a more comprehensive multiple model is obtained, ultimately leading to the suppression of the multiple results. This approach not only enhances the computational efficiency, accuracy, and applicability but also reduces the requirements for the original input data, thereby improving imaging accuracy, which can make it widely applicable in the processing of field data.

4. Numerical Examples

To verify the effectiveness of the iterative virtual events internal multiple suppression method, first of all, we use a simple horizontal layered model as a test case, with the velocity model displayed in Figure 5. The velocity model size is 4000 m × 1600 m, and forward modeling is carried out using the higher order finite difference method based on the acoustic equation. An absorption boundary is applied at the top of the model. At the acquisition surface, seismic data are modeled with 501 sources and 501 receivers on a fixed spread with a spacing of 20 m, and the grid size is 2.5 m × 2.5 m. The source emits a Ricker wavelet with a 25 Hz center frequency. The number of time sampling points is 1024, and the time sampling interval is 4 ms. Figure 6 shows the seismic record of one shot, including primaries and internal multiples in forward modeling, with the direct wave removed. Since the primary focus of this paper is the suppression of internal multiples in post-stack data, zero offset gathers are extracted from the forward shot records for internal multiple suppression. The extracted zero offset gathers are shown in Figure 7a.
The iterative virtual event multiple suppression method is applied to suppress the internal multiples, and the resulting zero offset gathers are displayed in Figure 7b. The figure reveals that all internal multiples are essentially suppressed without damaging the primaries, indicating that the iterative virtual event method can effectively suppress stacked multiples in model data. Figure 8a displays the predicted initial internal multiples model, while Figure 8b shows the internal multiples after iteration. The internal multiple information after iteration is more abundant and matches the actual multiples better. The suppression effect is sufficiently ideal, demonstrating that the iterative virtual event multiple suppression method can effectively suppress the internal multiples in synthetic data.

5. Field Example

Per the interpreter’s request, the primary focus is on suppressing internal multiples below 6 s to observe the original seismic profile more clearly. The seismic profile below 4 s is shown in Figure 9. From Figure 9a, it can be seen that the events below 6 s are similar in wave shape to the strong energy primaries, and the energy and frequency of these events are significantly different in the horizontal direction. Based on the relevant characteristics of multiples and multiples identification methods, these events are determined as internal multiples generated by the three sets of strong reflection interfaces of Middle Cambrian. Proof by facts, these internal multiples are challenging to suppress in pre-stack data using the filtering methods, such as the Radon transform method. Therefore, on top of the pre-stack suppression of multiples, we suppressed those internal multiples that are difficult to suppress in post-stack data using the iterative virtual events internal multiples suppression method.
Figure 9 shows the profile comparison before and after the internal multiple suppression. Figure 9a displays the partial stack profile before internal multiple suppression, while Figure 9b shows the corresponding stack profile after the internal multiple suppression. As seen in Figure 9b, the internal multiples below 6 s are mostly suppressed, and the horizontal difference in energy and frequency of the events below 6 s is smaller after the internal multiples are suppressed, particularly at the position indicated by the red arrow. To more clearly observe the effect of internal multiple suppression before and after, the part in the yellow box in Figure 9 is enlarged, as shown in Figure 10. In Figure 10, the event, which is similar to the upper strong reflection waveform with the opposite phase, periodic appearance, and significant lateral variation, is essentially suppressed. This is consistent with the analysis and identification results of multiples discussed earlier (Figure 2 and Figure 3). The results demonstrate that the iterative virtual event multiple suppression method can effectively suppress the internal multiples in field data.
To further analyze the suppression effect and accuracy of internal multiples, some sections before and after internal multiple suppression are selected for spectrum analysis, as shown in Figure 11. Figure 11a is the partial superimposed section before internal multiple suppression, Figure 11b is the partial superimposed section after internal multiple suppression, and Figure 11c shows the frequency distribution curve before and after internal multiple suppression. It can be seen from the figure that after suppressing internal multiples, the frequency spectrum is significantly broadened, and the spectral energy of the seismic profile is increased. This indicates that while using the iterative virtual event technology to suppress the internal multiples, the effective wave is protected, and the resolution of seismic data is improved. The method can provide reliable seismic data for subsequent studies of hydrocarbon source rocks.

6. Conclusions

The iterative virtual events method, when used to suppress internal multiples, involves first using the traditional virtual events method to build the initial multiple models, and then iterating the initial multiple models to complete the suppression of internal multiples. This approach not only simplifies the process of multiple suppression and improves computational efficiency and accuracy, but also obtains more abundant multiple information. From the model test, field data processing, and profile and spectrum analysis before and after processing, we can conclude the following:
(1)
The internal multiples in the post-stack data are suppressed using the iterative virtual event internal multiples suppression method, significantly reducing the number of events close to the primary shape.
(2)
After suppressing the internal multiples, the spectrum is significantly broadened, the overall energy is raised, and the resolution is improved.
(3)
The phenomenon of multiples’ energy and frequency varying greatly in the horizontal direction is weakened. The obtained seismic profile is more conforming to the actual geologic structure.
These results indicate that this method can effectively suppress the internal multiples in seismic data and can be widely used in the internal multiples suppression of field data.

Author Contributions

Conceptualization, M.X. and W.W.; methodology, J.X. and W.L.; software, J.S.; validation, W.L., J.S. and B.J.; formal analysis, Y.Z. and W.L.; investigation, T.Z. and Y.Z.; resources, Y.W. and B.J.; data curation, J.X., Y.W. and T.Z.; writing—original draft preparation, M.X. and Y.W.; writing—review and editing, W.W. and J.X.; visualization, B.J. and M.X.; supervision, W.W. and W.L.; project administration, Y.Z.; funding acquisition, T.Z. and W.W. All authors have read and agreed to the published version of the manuscript.

Funding

We thank the prospective basic projects of CNPC under grant no. 2021DJ-0304 and the National Natural Science Foundation of China under grant nos. 41974116 and 4227040468 for supporting this work.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Verschuur, D.J.; Berkhout, A.J.; Wapenaar, C.P.A. Adaptive surface-related multiple elimination. Geophysics 1992, 57, 1166–1177. [Google Scholar] [CrossRef]
  2. Van Groenestijn, G.J.A.; Verschuur, D.J. Estimation primaries by sparse inversion and application to near-offset data recon-struction. Geophysics 2009, 74, A23–A28. [Google Scholar] [CrossRef] [Green Version]
  3. Van Groenestjin, G.J.A.; Verschuur, D.J. Estimation primaries by sparse inversion near-offset reconstruction: Marine data applications. Geophysics 2009, 74, R119–R128. [Google Scholar] [CrossRef]
  4. Ma, J.; Chen, X.; Xue, Y. 3D surface-related multiple elimination. Oil Geophys. Prospect. 2015, 50, 14+33–40. [Google Scholar] [CrossRef]
  5. Zhu, F.; Cheng, J.; Wang, T.; Geng, J. Elimination of free-surface multiples using least-squares redatuming II: Marine seismic data acquired by streamers. Chin. J. Geophys. Chin. 2022, 65, 3123–3138. [Google Scholar] [CrossRef]
  6. Jakubowicz, H. Wave equation prediction and suppression of interbed multiples. In SEG Technical Program Expanded Abstracts; Society of Exploration Geophysicists: New Orleans, LA, USA, 1998; pp. 1527–1530. [Google Scholar] [CrossRef]
  7. Berkhout, A.J.; Verschuur, D.J. Removal of internal multiples with the common-focus-point (CFP) approach: Part 1—Explanation of the theor. Geophysics 2005, 70, V45–V60. [Google Scholar] [CrossRef] [Green Version]
  8. Verschuur, D.J.; Berkhout, A.J. Removal of internal multiples with the common-focus-point (CFP) approach: Part 2—Application strategies and data examples. Geophysics 2005, 70, V61–V72. [Google Scholar] [CrossRef] [Green Version]
  9. Weglein, A.B.; Nita, B.; Innanen, K.; Otnes, E.; Shaw, S.A.; Liu, F.; Zhang, H.; Ramírez, A.C.; Zhang, J.; Pavlis, G.L.; et al. Using the inverse scattering series to predict the wave-field at depth and the transmitted wave-field without an assumption about the phase of the measured reflection data or back propagation in the overburden. Geophysics 2006, 71, 1125–1137. [Google Scholar] [CrossRef] [Green Version]
  10. Ypma, F.H.C.; Verschuur, D.J. Estimating primaries by sparse inversion, a generalized approach. Geophys. Prospect. 2013, 61, 94–108. [Google Scholar] [CrossRef]
  11. Meles, G.A.; Löer, K.; Ravasi, M.; Curtis, A.; Filho, C.A.D.C. Internal multiple prediction and removal using Marchenko autofocusing and seismic interferometry. Geophysics 2015, 80, A7–A11. [Google Scholar] [CrossRef] [Green Version]
  12. Zhang, L.; Staring, M. Marchenko scheme based internal multiple reflection elimination in acoustic wavefield. J. Appl. Geophys. 2018, 159, 429–433. [Google Scholar] [CrossRef]
  13. Zhang, L.; Slob, E. A field data example of Marchenko multiple elimination. Geophysics 2020, 82, S65–S70. [Google Scholar] [CrossRef] [Green Version]
  14. Zhang, L.; Shao, J.; Zheng, Y.; Wang, Y.; Slob, E. A comparison of imaging domain and data domain Marchenko multiple elimination schemes. Chin. J. Geophys. Chin. 2022, 65, 3123–3138. [Google Scholar] [CrossRef]
  15. Zhang, M. A multiple suppression method based on self attention convolutional auto encoder. Geophys. Prospect. Pet. 2022, 61, 454–462. [Google Scholar]
  16. Ikelle, L.T. A construct of internal multiples from surface data only: The concept of virtual seismic events. Geophys. J. Int. 2006, 164, 383–393. [Google Scholar] [CrossRef] [Green Version]
  17. Ikelle, L.T.; Erez, I.; Yang, X. Scattering diagrams in seismic imaging: More insights into the construction of virtual events and internal multiples. J. Appl. Geophys. 2009, 67, 150–170. [Google Scholar] [CrossRef]
  18. Wu, J.; Wu, Z.; Hu, T.; He, Y.; Wang, P.; Yan, G.; Li, L. Seismic internal multiple attenuation based on constructing virtual events. Chin. J. Geophys. Chin. 2013, 56, 985–994. [Google Scholar] [CrossRef]
  19. Liu, J.; Hu, T.; Peng, G. Suppressing seismic inter-bed multiples with the adaptive virtual events method. Chin. J. Geophys. Chin. 2018, 61, 1196–1210. [Google Scholar] [CrossRef]
  20. Qu, Q.; Tang, Y.; Wang, H.; Lang, W. The petroleum system and objective interval for exploration of cambrian in southwest de-pression of Tarim Basin. Xing Jiang Pet. Geol. 2000, 21, 101–104+168–169. [Google Scholar]
  21. Cui, H.; Tian, L.; Liu, J.; Zhang, N. Ordovician tectonic activities and oil and gas accu-mulation in the Southwest Depression of Tarim Basin. Oil Geophys. Prospect. 2017, 52, 834–840+850. [Google Scholar]
  22. Zhang, L.; Han, E.; Zhu, L.; Zeng, C.; Fan, Q.; Wu, K.; Cao, Y.; Jiao, P. Characteristics of Evaporites Sedimentary Cycles and Its Controlling Factors of Paleocene Aertashi Formation in the Southwestern Tarim Depression. Acta Geol. Sin. 2015, 89, 2161–2170. [Google Scholar]
  23. Bao, P.; Wang, X.; Xie, J.; Wang, W.; Shi, Y.; Xu, J.; Yang, Y. Internal multiple suppression method based on iterative inversion. Chin. J. Geophys. Chin. 2021, 64, 2061–2072. [Google Scholar] [CrossRef]
  24. Bao, P.; Shi, Y.; Wang, W.; Zhang, W.; Pan, Z. An improved method for internal multiple elimination using the theory of virtu-al events. Pet. Sci. 2022, 19, 2663–2674. [Google Scholar] [CrossRef]
  25. Berkhout, A.J.; Verschuur, D.J. Estimation of multiple scattering by iterative inversion, Part I: Theoretical considerations. Geophysics 1997, 62, 1586–1595. [Google Scholar] [CrossRef]
  26. Verschuur, D.J.; Berkhout, A.J. Estimation of multiple scattering by iterative inversion, Part II: Practical aspects and examples. Geophysics 1997, 62, 1596–1611. [Google Scholar] [CrossRef] [Green Version]
Figure 1. Well A synthetic record calibration and well seismic profile.
Figure 1. Well A synthetic record calibration and well seismic profile.
Applsci 13 08832 g001
Figure 2. Velocity spectrum of the different original CMP trace sets; (a) velocity spectrum of CMP1; (b) velocity spectrum of CMP2.
Figure 2. Velocity spectrum of the different original CMP trace sets; (a) velocity spectrum of CMP1; (b) velocity spectrum of CMP2.
Applsci 13 08832 g002
Figure 3. Velocity spectrum of the different CMP gathers after multiple suppression by the Radon transform (corresponds to Figure 2); (a) velocity spectrum of CMP1; (b) velocity spectrum of CMP2.
Figure 3. Velocity spectrum of the different CMP gathers after multiple suppression by the Radon transform (corresponds to Figure 2); (a) velocity spectrum of CMP1; (b) velocity spectrum of CMP2.
Applsci 13 08832 g003
Figure 4. Construction of the virtual events and the internal multiples; (a) construction of the virtual events; (b) construction of the internal multiples. There pentagram represents convolution, the black line and black arrow represent the actual detectable seismic wave paths, red line and red arrow is used to assist in constructing virtual events and internal multiples.
Figure 4. Construction of the virtual events and the internal multiples; (a) construction of the virtual events; (b) construction of the internal multiples. There pentagram represents convolution, the black line and black arrow represent the actual detectable seismic wave paths, red line and red arrow is used to assist in constructing virtual events and internal multiples.
Applsci 13 08832 g004
Figure 5. Velocity model.
Figure 5. Velocity model.
Applsci 13 08832 g005
Figure 6. Original seismic data with the internal multiple.
Figure 6. Original seismic data with the internal multiple.
Applsci 13 08832 g006
Figure 7. Zero-offset data with the internal multiples suppression. (a) Original zero-offset data with the internal multiples; (b) result after the internal multiples are suppressed.
Figure 7. Zero-offset data with the internal multiples suppression. (a) Original zero-offset data with the internal multiples; (b) result after the internal multiples are suppressed.
Applsci 13 08832 g007
Figure 8. The predicted internal multiple. (a) The predicted internal multiples using the traditional virtual events method; (b) the predicted internal multiples after iterating.
Figure 8. The predicted internal multiple. (a) The predicted internal multiples using the traditional virtual events method; (b) the predicted internal multiples after iterating.
Applsci 13 08832 g008
Figure 9. Profile comparison before and after the internal multiple suppression. (a) Original stacked profile; (b) stacked section after the internal multiple suppression.
Figure 9. Profile comparison before and after the internal multiple suppression. (a) Original stacked profile; (b) stacked section after the internal multiple suppression.
Applsci 13 08832 g009
Figure 10. Part of profile comparison before and after the internal multiple suppression. (a) Original stacked profile; (b) stacked section after the internal multiple suppression.
Figure 10. Part of profile comparison before and after the internal multiple suppression. (a) Original stacked profile; (b) stacked section after the internal multiple suppression.
Applsci 13 08832 g010
Figure 11. Spectrum analysis before and after the internal multiple suppression. (a) Partial stacked sections before the internal multiple suppression; (b) partial stack section after the internal multiple suppression; (c) frequency distribution curve before and after the internal multiple suppression.
Figure 11. Spectrum analysis before and after the internal multiple suppression. (a) Partial stacked sections before the internal multiple suppression; (b) partial stack section after the internal multiple suppression; (c) frequency distribution curve before and after the internal multiple suppression.
Applsci 13 08832 g011
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Xiao, M.; Xie, J.; Wang, W.; Liu, W.; Sun, J.; Jin, B.; Zhang, T.; Zhao, Y.; Wang, Y. Application of Iterative Virtual Events Internal Multiple Suppression Technique: A Case of Southwest Depression Area of Tarim, China. Appl. Sci. 2023, 13, 8832. https://doi.org/10.3390/app13158832

AMA Style

Xiao M, Xie J, Wang W, Liu W, Sun J, Jin B, Zhang T, Zhao Y, Wang Y. Application of Iterative Virtual Events Internal Multiple Suppression Technique: A Case of Southwest Depression Area of Tarim, China. Applied Sciences. 2023; 13(15):8832. https://doi.org/10.3390/app13158832

Chicago/Turabian Style

Xiao, Mingtu, Junfa Xie, Weihong Wang, Wenqing Liu, Jiaqing Sun, Baozhong Jin, Tao Zhang, Yuhe Zhao, and Yihui Wang. 2023. "Application of Iterative Virtual Events Internal Multiple Suppression Technique: A Case of Southwest Depression Area of Tarim, China" Applied Sciences 13, no. 15: 8832. https://doi.org/10.3390/app13158832

APA Style

Xiao, M., Xie, J., Wang, W., Liu, W., Sun, J., Jin, B., Zhang, T., Zhao, Y., & Wang, Y. (2023). Application of Iterative Virtual Events Internal Multiple Suppression Technique: A Case of Southwest Depression Area of Tarim, China. Applied Sciences, 13(15), 8832. https://doi.org/10.3390/app13158832

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop