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.展开更多
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.展开更多
An accurate and wide-angle one-way propagator for wavefield extrapolation is an important topic for research on wave-equation prestack depth migration in the presence of large and rapid velocity variations. Based on t...An accurate and wide-angle one-way propagator for wavefield extrapolation is an important topic for research on wave-equation prestack depth migration in the presence of large and rapid velocity variations. Based on the optimal separable approximation presented in this paper, the mixed domain algorithm with forward and inverse Fourier transforms is used to construct the 3D one-way wavefield extrapolation operator. This operator separates variables in the wavenumber and spatial domains. The phase shift operation is implemented in the wavenumber domain while the time delay for lateral velocity variation is corrected in the spatial domain. The impulse responses of the one-way wave operator show that the numeric computation is consistent with the theoretical value for each velocity, revealing that the operator constructed with the optimal separable approximation can be applied to lateral velocity variations for the case of small steps. Imaging results of the SEG/EAGE model and field data indicate that the new method can be used to image complex structure.展开更多
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.展开更多
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.展开更多
垂直地震剖面(Vertical Seismic Profiling,VSP)资料处理中波场分离是关键问题之一.随着属性提取技术的发展,新的属性参数(例如Q值)提取技术对波场分离的保真性要求越来越高.本文改进了传统奇异值分解(Singular Value Decomposition,SVD...垂直地震剖面(Vertical Seismic Profiling,VSP)资料处理中波场分离是关键问题之一.随着属性提取技术的发展,新的属性参数(例如Q值)提取技术对波场分离的保真性要求越来越高.本文改进了传统奇异值分解(Singular Value Decomposition,SVD)法,给出了一种对波场的动力学特征具有更好的保真性,可以作为Q值提取的预处理步骤的零偏VSP资料上下行波场分离方法.该方法通过两步奇异值分解变换实现:第一步,排齐下行波同相轴,利用SVD变换压制部分下行波能量;第二步,在剩余波场中排齐上行波同相轴,使用SVD变换提取上行波场.在该方法的实现过程中,压制部分下行波能量后的剩余波场中仍然存在较强的下行波干扰,使得上行波同相轴的排齐比较困难.本文给出了一种通过极大化多道数据线性相关程度(Maximize Coherence,MC)排齐同相轴的算法,在一定程度上解决了低信噪比下排齐同相轴的问题.将本文提出的方法用于合成数据和实际资料的处理,并与传统SVD法的处理结果进行对比,结果表明本文提出的波场分离方法具有良好的保真性,得到波场的质量明显优于传统SVD法.通过对本文方法和传统SVD法处理合成数据得到的下行波场提取Q值,然后进行对比可知,本文方法可以有效提高所提取Q值的准确性,适合作为Q值提取的预处理步骤.展开更多
基金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(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 China National Natural Science Foundation （N0. 40474047）.
文摘An accurate and wide-angle one-way propagator for wavefield extrapolation is an important topic for research on wave-equation prestack depth migration in the presence of large and rapid velocity variations. Based on the optimal separable approximation presented in this paper, the mixed domain algorithm with forward and inverse Fourier transforms is used to construct the 3D one-way wavefield extrapolation operator. This operator separates variables in the wavenumber and spatial domains. The phase shift operation is implemented in the wavenumber domain while the time delay for lateral velocity variation is corrected in the spatial domain. The impulse responses of the one-way wave operator show that the numeric computation is consistent with the theoretical value for each velocity, revealing that the operator constructed with the optimal separable approximation can be applied to lateral velocity variations for the case of small steps. Imaging results of the SEG/EAGE model and field data indicate that the new method can be used to image complex structure.
基金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.
基金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.
