Prestack reverse time migration（PSTM） is a common imaging method; however low-frequency noises reduce the structural imaging precision. Thus, the suppression of migration noises must be considered. The generation me...Prestack reverse time migration（PSTM） is a common imaging method; however low-frequency noises reduce the structural imaging precision. Thus, the suppression of migration noises must be considered. The generation mechanism of low-frequency noises is analyzed and the up-, down-, left-, and right-going waves are separated using the Poynting vector of the acoustic wave equation. The computational complexity and memory capacitance of the proposed method are far smaller than that required when using the conventional separation algorithm of 2D Fourier transform. The normalized wavefield separation crosscorrelation imaging condition is used to suppress low-frequency noises in reverse time migration and improve the imaging precision. Numerical experiments using the Marmousi model are performed and the results show that the up-, down-, left-, and right-going waves are well separated in the continuation of the wavefield using the Poynting vector. We compared the imaging results with the conventional method, Laplacian filtering, and wavefield separation with the 2D Fourier transform. The comparison shows that the migration noises are well suppressed using the normalized wavefield separation cross-correlation imaging condition and higher precision imaging results are obtained.展开更多
Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumb...Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.展开更多
In multi-component seismic exploration, the horizontal and vertical components both contain P- and SV-waves. The P- and SV-wavefields in a seismic record can be separated by their horizontal and vertical displacements...In multi-component seismic exploration, the horizontal and vertical components both contain P- and SV-waves. The P- and SV-wavefields in a seismic record can be separated by their horizontal and vertical displacements when upgoing P- and SV-waves arrive at the sea floor. If the sea floor P wave velocity, S wave velocity, and density are known, the separation can be achieved in ther-p domain. The separated wavefields are then transformed to the time domain. A method of separating P- and SV-wavefields is presented in this paper and used to effectively separate P- and SV-wavefields in synthetic and real data. The application to real data shows that this method is feasible and effective. It also can be used for free surface data.展开更多
Internal multiple interference,affecting both seismic data processing and interpretation,has been observed for long time.Although great progress has been achieved in developing a variety of internal-multiple-eliminati...Internal multiple interference,affecting both seismic data processing and interpretation,has been observed for long time.Although great progress has been achieved in developing a variety of internal-multiple-elimination(IME)methods,how to increase accuracy and reduce cost of IME still poses a significant challenge.A new method is proposed to effectively and efficiently eliminate internal multi-ples,along with its application in internal-multiple-eliminated-migration(IMEM),addressing this issue.This method stems from two-way wave equation depth-extrapolation scheme and associated up/down wavefield separation,which can accomplish depth-extrapolation of both up-going and down-going wavefields simultaneously,and complete internal-multiple-elimination processing,adaptively and effi-ciently.The proposed method has several features:(1)input data is same as that for conventional migration:source signature(used for migration only),macro velocity model,and receiver data,without additional requirements for source/receiver sampling;(2)method is efficient,without need of iterative calculations(which are typically needed for most of IME algorithms);and(3)method is cost effective:IME is completed in the same depth-extrapolation scheme of IMEM,without need of a separate pro-cessing and additional cost.Several synthesized data models are used to test the proposed method:one-dimensional model,horizontal layered model,multi-layer model with one curved layer,and SEG/EAGE Salt model.Additionally,we perform a sensitivity analysis of velocity using smoothed models.This analysis reveals that although the accuracy of velocity measurements impacts our proposed method,it significantly reduces internal multiple false imaging compared to traditional RTM techniques.When applied to actual seismic data from a carbonate reservoir zone,our method demonstrates superior clarity in imaging results,even in the presence of high-velocity carbonate formations,outperforming conven-tional migration methods in deep strata.展开更多
The conventional fast converted-wave imaging method directly uses backward Pand converted S-wavefield to produce joint images. However, this image is accompanied by strong background noises, because the wavefi elds in...The conventional fast converted-wave imaging method directly uses backward Pand converted S-wavefield to produce joint images. However, this image is accompanied by strong background noises, because the wavefi elds in all propagation directions contribute to it. Given this issue, we improve the conventional imaging method in the two aspects. First, the amplitude-preserved P-and S-wavef ield are obtained by using an improved space-domain wavef ield separation scheme to decouple the original elastic wavef ield. Second, a convertedwave imaging condition is constructed based on the directional-wavefield separation and only the wavefields propagating in the same directions used for cross-correlation imaging, resulting in effectively eliminating the imaging artifacts of the wavefields with different directions;Complex-wavefi eld extrapolation is adopted to decompose the decoupled P-and S-wavefield into directional-wavefields during backward propagation, this improves the eff iciency of the directional-wavef ield separation. Experiments on synthetic data show that the improved method generates more accurate converted-wave images than the conventional one. Moreover, the improved method has application potential in micro-seismic and passive-source exploration due to its source-independent characteristic.展开更多
Fracture-cave reservoirs in carbonate rocks are characterized by a large difference in fracture and cavity size,and a sharp variation in lithology and velocity,thereby resulting in complex diffraction responses.Some s...Fracture-cave reservoirs in carbonate rocks are characterized by a large difference in fracture and cavity size,and a sharp variation in lithology and velocity,thereby resulting in complex diffraction responses.Some small-scale fractures and caves cause weak diffraction energy and would be obscured by the continuous reflection layer in the imaging section,thereby making them difficult to identify.This paper develops a diffraction wave imaging method in the dip domain,which can improve the resolution of small-scale diffractors in the imaging section.Common imaging gathers(CIGs)in the dip domain are extracted by Gaussian beam migration.In accordance with the geometric differences of the diffraction being quasilinear and the reflection being quasiparabolic in the dip-domain CIGs,we use slope analysis technique to filter waves and use Hanning window function to improve the diffraction wave separation level.The diffraction dip-domain CIGs are stacked horizontally to obtain diffraction imaging results.Wavefield separation analysis and numerical modeling results show that the slope analysis method,together with Hanning window filtering,can better suppress noise to obtain the diffraction dip-domain CIGs,thereby improving the clarity of the diffractors in the diffraction imaging section.展开更多
基金supported by the National Natural Science Foundation of China(No.41174087,41204089)the National Oil and Gas Major Project(No.2011ZX05005-005)
文摘Prestack reverse time migration（PSTM） is a common imaging method; however low-frequency noises reduce the structural imaging precision. Thus, the suppression of migration noises must be considered. The generation mechanism of low-frequency noises is analyzed and the up-, down-, left-, and right-going waves are separated using the Poynting vector of the acoustic wave equation. The computational complexity and memory capacitance of the proposed method are far smaller than that required when using the conventional separation algorithm of 2D Fourier transform. The normalized wavefield separation crosscorrelation imaging condition is used to suppress low-frequency noises in reverse time migration and improve the imaging precision. Numerical experiments using the Marmousi model are performed and the results show that the up-, down-, left-, and right-going waves are well separated in the continuation of the wavefield using the Poynting vector. We compared the imaging results with the conventional method, Laplacian filtering, and wavefield separation with the 2D Fourier transform. The comparison shows that the migration noises are well suppressed using the normalized wavefield separation cross-correlation imaging condition and higher precision imaging results are obtained.
基金supported by the National Science Foundation of China(No.41174100)the Large-scale Oil and Gas Field and Coalbed Methane Development Major Projects(No.2011ZX05019-008-08)the China National Petroleum Corporation(No.2014A-3609)
文摘Wavefield separation of multicomponent seismic data to image subsurface structures can be realized in either the space domain or the wavenumber domain. However, as the particle velocity components used in the wavenumber-domain wavefield separation are not defined at the same grid point with the staggered-grid finite-difference method for elastic wavefield simulation, we propose the wavenumber-domain interpolation method to estimate the required values at the common grid points prior to the wavenumber-domain true-amplitude wavefield separation. Moreover, numerical experiments show that the wavenumber-domain interpolation method has high interpolation accuracy and the trueamplitude wavefield separation method shows good amplitude preservation. The application of the proposed methodology to elastic reverse-time migration can obtain good amplitudepreserved images even in the case of some velocity error.
基金This research is sponsored by National Natural Science Foundation of China （No. 40272041） and Innovative Foundation of CNPC （N0. 04E702）.
文摘In multi-component seismic exploration, the horizontal and vertical components both contain P- and SV-waves. The P- and SV-wavefields in a seismic record can be separated by their horizontal and vertical displacements when upgoing P- and SV-waves arrive at the sea floor. If the sea floor P wave velocity, S wave velocity, and density are known, the separation can be achieved in ther-p domain. The separated wavefields are then transformed to the time domain. A method of separating P- and SV-wavefields is presented in this paper and used to effectively separate P- and SV-wavefields in synthetic and real data. The application to real data shows that this method is feasible and effective. It also can be used for free surface data.
基金supported by the National Natural Science Foundation of China(Grant No.42004103)Sichuan Science and Technology Program(2023NSFSC0257)the CNPC Innovation Found(2022DQ02-0306).
文摘Internal multiple interference,affecting both seismic data processing and interpretation,has been observed for long time.Although great progress has been achieved in developing a variety of internal-multiple-elimination(IME)methods,how to increase accuracy and reduce cost of IME still poses a significant challenge.A new method is proposed to effectively and efficiently eliminate internal multi-ples,along with its application in internal-multiple-eliminated-migration(IMEM),addressing this issue.This method stems from two-way wave equation depth-extrapolation scheme and associated up/down wavefield separation,which can accomplish depth-extrapolation of both up-going and down-going wavefields simultaneously,and complete internal-multiple-elimination processing,adaptively and effi-ciently.The proposed method has several features:(1)input data is same as that for conventional migration:source signature(used for migration only),macro velocity model,and receiver data,without additional requirements for source/receiver sampling;(2)method is efficient,without need of iterative calculations(which are typically needed for most of IME algorithms);and(3)method is cost effective:IME is completed in the same depth-extrapolation scheme of IMEM,without need of a separate pro-cessing and additional cost.Several synthesized data models are used to test the proposed method:one-dimensional model,horizontal layered model,multi-layer model with one curved layer,and SEG/EAGE Salt model.Additionally,we perform a sensitivity analysis of velocity using smoothed models.This analysis reveals that although the accuracy of velocity measurements impacts our proposed method,it significantly reduces internal multiple false imaging compared to traditional RTM techniques.When applied to actual seismic data from a carbonate reservoir zone,our method demonstrates superior clarity in imaging results,even in the presence of high-velocity carbonate formations,outperforming conven-tional migration methods in deep strata.
基金supported by the National Science and Technology Major Project of China(No.2017ZX05018-005)National Natural Science Foundation of China(No.41474110)
文摘The conventional fast converted-wave imaging method directly uses backward Pand converted S-wavefield to produce joint images. However, this image is accompanied by strong background noises, because the wavefi elds in all propagation directions contribute to it. Given this issue, we improve the conventional imaging method in the two aspects. First, the amplitude-preserved P-and S-wavef ield are obtained by using an improved space-domain wavef ield separation scheme to decouple the original elastic wavef ield. Second, a convertedwave imaging condition is constructed based on the directional-wavefield separation and only the wavefields propagating in the same directions used for cross-correlation imaging, resulting in effectively eliminating the imaging artifacts of the wavefields with different directions;Complex-wavefi eld extrapolation is adopted to decompose the decoupled P-and S-wavefield into directional-wavefields during backward propagation, this improves the eff iciency of the directional-wavef ield separation. Experiments on synthetic data show that the improved method generates more accurate converted-wave images than the conventional one. Moreover, the improved method has application potential in micro-seismic and passive-source exploration due to its source-independent characteristic.
基金funded jointly by the National Natural Science Foundation of China(No.41104069)Shandong Province Higher Educational Science and Technology Program(No.J17KA197)+1 种基金Open Foundation of Shandong Provincial Key Laboratory of Depositional Mineralization&Sedimentary Minerals of Shandong University of Science and Technology(No.DMSM2018018)Chunhui Research Foundation of Shengli College,China University of Petroleum(No.KY2017007)。
文摘Fracture-cave reservoirs in carbonate rocks are characterized by a large difference in fracture and cavity size,and a sharp variation in lithology and velocity,thereby resulting in complex diffraction responses.Some small-scale fractures and caves cause weak diffraction energy and would be obscured by the continuous reflection layer in the imaging section,thereby making them difficult to identify.This paper develops a diffraction wave imaging method in the dip domain,which can improve the resolution of small-scale diffractors in the imaging section.Common imaging gathers(CIGs)in the dip domain are extracted by Gaussian beam migration.In accordance with the geometric differences of the diffraction being quasilinear and the reflection being quasiparabolic in the dip-domain CIGs,we use slope analysis technique to filter waves and use Hanning window function to improve the diffraction wave separation level.The diffraction dip-domain CIGs are stacked horizontally to obtain diffraction imaging results.Wavefield separation analysis and numerical modeling results show that the slope analysis method,together with Hanning window filtering,can better suppress noise to obtain the diffraction dip-domain CIGs,thereby improving the clarity of the diffractors in the diffraction imaging section.
