Amplitude dissipation and phase dispersion occur when seismic waves propagate in attenuated anisotropic media,affecting the quality of migration imaging.To compensate and correct for these effects,the fractional Lapla...Amplitude dissipation and phase dispersion occur when seismic waves propagate in attenuated anisotropic media,affecting the quality of migration imaging.To compensate and correct for these effects,the fractional Laplacian pure viscoacoustic wave equation capable of producing stable and noise-free wavefields has been proposed and implemented in the Q-compensated reverse time migration(RTM).In addition,the second-order Taylor series expansion is usually adopted in the hybrid finite-difference/pseudo-spectral(HFDPS)strategy to solve spatially variable fractional Laplacian.However,during forward modeling and Q-compensated RTM,this HFDPS strategy requires 11 and 17 fast Fourier transforms(FFTs)per time step,respectively,leading to computational inefficiency.To improve computational efficiency,we introduce two high-efficiency HFDPS numerical modeling strategies based on asymptotic approximation and algebraic methods.Through the two strategies,the number of FFTs decreased from 11 to 6 and 5 per time step during forward modeling,respectively.Numerical examples demonstrate that wavefields simulated using the new numerical modeling strategies are accurate and highly efficient.Finally,these strategies are employed for implementing high-efficiency and stable Q-compensated RTM techniques in tilted transversely isotropic media,reducing the number of FFTs from 17 to 9 and 8 per time step,respectively,significantly improving computational efficiency.Synthetic data examples illustrate the effectiveness of the proposed Q-compensated RTM scheme in compensating amplitude dissipation and correcting phase distortion.展开更多
Reasonable field acquisition geometry can not only guide seismic exploration to obtain sufficient geological information of target body,but also reduce acquisition cost to the maximum.In this study,building on convent...Reasonable field acquisition geometry can not only guide seismic exploration to obtain sufficient geological information of target body,but also reduce acquisition cost to the maximum.In this study,building on conventional ray-based geometry design methods,we incorporate imaging results as a constraint to optimize the geometry design and evaluate its effectiveness.Firstly,the geological model of the target layer is established based on the geological data of the study area and surface seismic data combined with exploration tasks.Then,the ray-tracing method is employed to simulate and assess the proposed geometry design,verifying whether its parameters meet the exploration requirements.Finally,the imaging effect of the designed geometry on the target layer is tested by the cross-well seismic reverse time migration method.This methodology was applied to design the cross-well seismic acquisition geometry for offshore deviated wells in the X Oilfield.The simulation results demonstrate that the imaging-driven geometry design approach effectively guides field operations,enhances the imaging quality of the target layer,and reduces acquisition costs.展开更多
The East China Sea Shelf Basin,especially the Xihu depression in its northeastern region,has long been regarded as an important target area for oil and gas exploration in China.Since the development of the late Cretac...The East China Sea Shelf Basin,especially the Xihu depression in its northeastern region,has long been regarded as an important target area for oil and gas exploration in China.Since the development of the late Cretaceous,this region has experienced complex and multistage tectonic movements,such as extensional faulting and compressional inversion,forming its current unique geological structure pattern.As one of the largest Mesozoic–Cenozoic sedimentary basins offshore China,the East China Sea Shelf Basin not only contains abundant oil and gas resources but also occupies a key strategic position.In addition,it is gradually becoming an important base and strategic core area for China's oil and gas resource development.Although remarkable progress has been achieved in oil and gas exploration,seismic exploration in the Xihu depression still faces many challenges.Especially in deep target layers,weak seismic reflection energy,low signal-to-noise ratio,and poor wave group continuity seriously affect the efficiency and accuracy of oil and gas exploration.In addition,given the considerable influence of folds and faults in the East China Sea Xihu Basin,transverse anisotropy is widespread,which causes difficulty for conventional imaging technologies based on isotropic assumptions to achieve accurate imaging.Therefore,to address these challenges and meet the need for accurate imaging of complex structures in the Xihu depression of the East China Sea Basin,it is essential to advance research on anisotropic imaging technologies.This study responds to that need by integrating the reverse time migration method under a tilted transversely isotropic(TTI)medium model and applying it to the imaging of the Xihu depression.The imaging profiles based on the TTI medium can substantially improve the imaging accuracy and signal-to-noise ratio,effectively improve the imaging quality of the fault systems and mid-deep strata in the area,and provide new technical support and progress momentum for oil and gas exploration in the Xihu depression of the East China Sea.展开更多
The numerical dispersion phenomenon in the finite-difference forward modeling simulations of the wave equation significantly affects the imaging accuracy in acoustic reflection logging.This issue is particularly prono...The numerical dispersion phenomenon in the finite-difference forward modeling simulations of the wave equation significantly affects the imaging accuracy in acoustic reflection logging.This issue is particularly pronounced in the reverse time migration(RTM)method used for shear-wave(S-wave)logging imaging.This not only affects imaging accuracy but also introduces ambiguities in the interpretation of logging results.To address this challenge,this study proposes the use of a least-squares difference coefficient optimization algorithm aiming to suppress the numerical dispersion phenomenon in the RTM of S-wave reflection imaging logging.By optimizing the difference coefficients,the high-precision finite-difference algorithm serves as an effective operator for both forward and backward RTM processes.This approach is instrumental in eliminating migration illusions,which are often caused by numerical dispersion.The effectiveness of this optimized algorithm is demonstrated through numerical results,which indicate that it can achieve more accurate forward imaging results across various conditions,including high-and low-velocity strata,and is effective in both large and small spatial grids.The results of processing real data demonstrate that numerical dispersion optimization effectively reduces migration artifacts and diminishes ambiguities in logging interpretations.This optimization offers crucial technical support to the RTM method,enhancing its capability for accurately modeling and imaging S-wave reflections.展开更多
Seismic forward-prospecting in tunnels is an important step to ensure excavation safety. Nowadays, most advanced imaging techniques in seismic exploration involve calculating the solution of elastic wave equation in a...Seismic forward-prospecting in tunnels is an important step to ensure excavation safety. Nowadays, most advanced imaging techniques in seismic exploration involve calculating the solution of elastic wave equation in a certain coordinate system. However, considering the cylindrical geometry of common tunnel body, Cartesian coordinate system seemingly has limited applicability in tunnel seismic forwardprospecting. To accurately simulate the seismic signal received in tunnels, previous imaging method using decoupled non-conversion elastic wave equation is extended from Cartesian coordinates to cylindrical coordinates. The proposed method preserves the general finite-difference time-domain(FDTD)scheme in Cartesian coordinates, except for a novel wavefield calculation strategy addressing the singularity issue inherited at the cylindrical axis. Moreover, the procedure of cylindrical elastic reverse time migration(CERTM) in tunnels is introduced based on the decoupled non-conversion elastic wavefield. Its imaging effect is further validated via numerical experiments on typical tunnel models. As indicated in the synthetic examples, both the PP-and SS-images could clearly show the geological structure in front of the tunnel face without obvious crosstalk artifacts. Migration imaging using PP-waves can present satisfactory results with higher resolution information supplemented by the SS-images. The potential of applying the proposed method in real-world cases is demonstrated in a water diversion tunnel. In the end, we share our insights regarding the singularity treatment and further improvement of the proposed method.展开更多
Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of t...Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.展开更多
Anisotropy correction is necessary during the processing of converted PS- wave seismic data to achieve accurate structural imaging, reservoir prediction, and fracture detection. To effectively eliminate the adverse ef...Anisotropy correction is necessary during the processing of converted PS- wave seismic data to achieve accurate structural imaging, reservoir prediction, and fracture detection. To effectively eliminate the adverse effects of S-wave splitting and to improve PS- wave imaging quality, we tested methods for pre-stack migration imaging and anisotropic correction of PS-wave data. We based this on the propagation rules of seismic waves in a horizontal transverse isotropy medium, which is a fractured medium model that reflects likely subsurface conditions in the field. We used the radial （R） and transverse （T） components of PS-wave data to separate the fast and slow S-wave components, after which their propagation moveout was effectively extracted. Meanwhile, corrections for the energies and propagation moveouts of the R and T components were implemented using mathematical rotation. The PS-wave imaging quality was distinctly improved, and we demonstrated the reliability of our methods through numerical simulations. Applying our methods to three-dimensional and three-component seismic field data from the Xinchang-Hexingchang region of the Western Sichuan Depression in China, we obtained high-quality seismic imaging with continuous reflection wave groups, distinct structural features, and specific stratigraphic contact relationships. This study provides an effective and reliable approach for data processing that will improve the exploration of complex, hidden lithologic gas reservoirs.展开更多
In this paper, we use the method of pseudo-offset migration （POM） to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into co...In this paper, we use the method of pseudo-offset migration （POM） to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into common conversion scatter point （CCSP） gathers directly by POM, which simplifies the conventional processing procedure for converted waves. The POM gather fold and SNR are high, which is favorable for velocity analysis and especially suitable for seismic data with low SNR. We used equivalent anisotropic theory to compute anisotropic parameters. Based on the scattering wave traveltime equation in a VTI medium, the POM pseudo-offset migration in anisotropic media was deduced. By amplitude-preserving POM gather mapping, velocity analysis, stack processing, and so on, the anisotropic migration results were acquired. The forward modeling computation and actual data processing demonstrate the validity of converted wave pre-stack time migration with amplitude-preservation using the anisotropic POM method.展开更多
Although conventional reverse time migration can be perfectly applied to structural imaging it lacks the capability of enabling detailed delineation of a lithological reservoir due to irregular illumination. To obtain...Although conventional reverse time migration can be perfectly applied to structural imaging it lacks the capability of enabling detailed delineation of a lithological reservoir due to irregular illumination. To obtain reliable reflectivity of the subsurface it is necessary to solve the imaging problem using inversion. The least-square reverse time migration （LSRTM） （also known as linearized refleetivity inversion） aims to obtain relatively high-resolution amplitude preserving imaging by including the inverse of the Hessian matrix. In practice, the conjugate gradient algorithm is proven to be an efficient iterative method for enabling use of LSRTM. The velocity gradient can be derived from a cross-correlation between observed data and simulated data, making LSRTM independent of wavelet signature and thus more robust in practice. Tests on synthetic and marine data show that LSRTM has good potential for use in reservoir description and four-dimensional （4D） seismic images compared to traditional RTM and Fourier finite difference （FFD） migration. This paper investigates the first order approximation of LSRTM, which is also known as the linear Born approximation. However, for more complex geological structures a higher order approximation should be considered to improve imaging quality.展开更多
Engineering seismic exploration aims at shallow imaging which is confused by statics if the surface is uneven. Direct pre-stack depth migration （DPDM） is based on accurate elevations of sources and receivers, by whi...Engineering seismic exploration aims at shallow imaging which is confused by statics if the surface is uneven. Direct pre-stack depth migration （DPDM） is based on accurate elevations of sources and receivers, by which static correction is completely abandoned before migration and surely the imaging quality is remarkably improved. To obtain some artificial shot gathers, high-order staggered-grid finite-difference （FD） method is adapted to model acoustic wave propagation. Since the shot gathers are always disturbed by regular interferences, the statics still must be applied to supporting the interference elimination by apparent velocity filtering method. Then all the shot gathers should be removed back to their original positions by reverse statics. Finally, they are migrated by pre-stack reverse-time depth migration and imaged. The numerical experiments show that the DPDM can ideally avoid the mistakes caused by statics and increase imaging precision.展开更多
To support amplitude variation with offset （AVO） analysis in complex structure areas, we introduce an amplitude-preserving plane-wave prestack time migration approach based on the double-square-root wave equation in...To support amplitude variation with offset （AVO） analysis in complex structure areas, we introduce an amplitude-preserving plane-wave prestack time migration approach based on the double-square-root wave equation in media with little lateral velocity variation. In its implementation, a data mapping algorithm is used to obtain offset-plane-wave data sets from the common-midpoint gathers followed by a non-recursive phase-shift solution with amplitude correction to generate common-image gathers in offset-ray-parameter domain and a structural image. Theoretical model tests and a real data example show that our prestack time migration approach is helpful for AVO analysis in complex geological environments.展开更多
Elastic migration has been widely paid attention by employing the vector processing of mul- ticomponent seismic data. Ray based elastic Kirchhoff migration has such properties as high flexibility and high efficiency. ...Elastic migration has been widely paid attention by employing the vector processing of mul- ticomponent seismic data. Ray based elastic Kirchhoff migration has such properties as high flexibility and high efficiency. However, it has failed to solve many problems caused by multipath. On the other hand, elastic reverse-time migration （RTM） based on the two-way wave equation is known to be capable of dealing with these problems, but it is extremely expensive when applied in 3D cases and velocity model building. Based on the elastic Kirchhoff-Helmholtz integral, we calculate deeoupled backward-continued wavefields by introducing elastic Green functions for P- and S-waves, which is expressed by a summation of elastodynamic Gaussian beams. The PP and polarity-corrected PS images are obtained by calculating the correlation between downward and deeoupled backward-continued vector wavefields, where polarity correction is performed by analyzing the relation between the polarization direction of converted PS waves and incident angle on the interface. To a large extent, our method combines the high efficiency of ray-based migration with the high accuracy of wave-equation based reverse-time migration. Application of this method to multicomponent synthetic datasets from the fault model and Marmousi 2 model demonstrates the validity, flexibility and accuracy of the new method.展开更多
With an increasing demand for high-resolution imaging of complex subsurface structures, thin layers and hidden reservoirs, borehole and cross-well seismic migration methods have become important. However, large differ...With an increasing demand for high-resolution imaging of complex subsurface structures, thin layers and hidden reservoirs, borehole and cross-well seismic migration methods have become important. However, large differences are observed in the frequency bandwidth between the surface, borehole, and cross-well surveys. Thus, variable-grid- based algorithms have been adapted to reverse-time migration. Further, we introduce Lanczos filtering to ensure the stability of wavefield calculations as well as to decrease the artificial reflections that are caused due to the variable grid size. Finally, we observe that the application of this method to surface survey, borehole, and cross-well seismic data suggests improvements in the delineation of minor fractures and steeply dipping faults.展开更多
Pseudo-offset migration （POM） is a new method for prestack time migration of converted waves that improves on equivalent-offset migration （EOM）. The mapping of POM is different than EOM but the purpose of the two ...Pseudo-offset migration （POM） is a new method for prestack time migration of converted waves that improves on equivalent-offset migration （EOM）. The mapping of POM is different than EOM but the purpose of the two methods is to map the input samples to the common conversion scatter point （CCSP） gathers. This paper introduces the principles of the two migration methods and model parameter sensitivity tests for both POM and EOM. At large offset-to-depth ratios the hyperbolic approximation, the three-term approximation and the double square root （DSR） equation are used to NMO-correct the mapped POM gathers to obtain more accurate migration velocities. These equations were derived and calculated by small pseudo offset. POM is then used to image complex structure and prestack time migration.展开更多
With the development of seismic exploration,passive-source seismic data has attracted increasing attention.Ambient noise passive seismic sources exists widely in nature and industrial production.Passive seismic data i...With the development of seismic exploration,passive-source seismic data has attracted increasing attention.Ambient noise passive seismic sources exists widely in nature and industrial production.Passive seismic data is important in logging while drilling(LWD),large-scale structural exploration,etc.In this paper,we proposed a passive multiple reverse time migration imaging(PMRTMI)method based on wavefield decomposition and normalized imaging conditions method.This method differs from seismic interferometry in that it can use raw passive seismic data directly in RTM imaging without reconstruction of virtual active gather,and we use the wavefield decomposition method to eliminate the low frequency noise in RTM.Further,the energy normalized imaging condition is used in full wavefield decomposition,which can not only enhance the image quality of both edge and deep information but also overcome the wrong energy problem caused by uneven distribution of passive sources;furthermore,this method exhibits high efficiency.Finally,numerical examples with the Marmousi model show the effectiveness of the method.展开更多
Intrinsic attenuation of the earth causes energy loss and phase distortion in seismic wave propagation.To obtain high-resolution imaging results,these negative effects must be considered during reverse time migration(...Intrinsic attenuation of the earth causes energy loss and phase distortion in seismic wave propagation.To obtain high-resolution imaging results,these negative effects must be considered during reverse time migration(RTM).We can easily implement attenuation-compensated RTM using the constant Q viscoacoustic wave equation with decoupled amplitude attenuation and phase dispersion terms.However,the nonphysical amplitude-compensation process will inevitably amplify the high-frequency noise in the wavefield in an exponential form,causing the numerical simulation to become unstable.This is due to the fact that the amplitude of the compensation grows exponentially with frequency.In order to achieve stable attenuation-compensated RTM,we modify the analytic expression of the attenuation compensation extrapolation operator and make it only compensate for amplitude loss within the effective frequency band.Based on this modified analytic formula,we then derive an explicit time-space domain attenuation compensation extrapolation operator.Finally,the implementation procedure of stable attenuation-compensated RTM is presented.In addition to being simple to implement,the newly proposed attenuation-compensated extrapolation operator is superior to the conventional low-pass filter in suppressing random noise,which will further improve the imaging resolution.We use two synthetic and one land seismic datasets to verify the stability and effectiveness of the proposed attenuationcompensated RTM in improving imaging resolution in viscous media.展开更多
In marine seismic exploration,ocean bottom cable technology can record multicomponent seismic data for multiparameter inversion and imaging.This study proposes an elastic multiparameter lease-squares reverse time migr...In marine seismic exploration,ocean bottom cable technology can record multicomponent seismic data for multiparameter inversion and imaging.This study proposes an elastic multiparameter lease-squares reverse time migration based on the ocean bottom cable technology.Herein,the wavefield continuation operators are mixed equations:the acoustic wave equations are used to calculate seismic wave propagation in the seawater medium,whereas in the solid media below the seabed,the wavefields are obtained by P-and S-wave separated vector elastic wave equations.At the seabed interface,acoustic–elastic coupling control equations are used to combine the two types of equations.P-and S-wave separated elastic migration operators,demigration operators,and gradient equations are derived to realize the elastic least-squares reverse time migration based on the P-and S-wave mode separation.The model tests verify that the proposed method can obtain high-quality images in both the P-and S-velocity components.In comparison with the traditional elastic least-squares reverse time migration method,the proposed method can readily suppress imaging crosstalk noise from multiparameter coupling.展开更多
In seasonally breeding birds, the annual cycle of photoperiod is a principal environmental cue for temporal arrange- ment of different life-history stages, such as migration and breeding. In the past, most research ha...In seasonally breeding birds, the annual cycle of photoperiod is a principal environmental cue for temporal arrange- ment of different life-history stages, such as migration and breeding. In the past, most research has focused on the mechanisms of photoperiodic control of breeding with less attention paid to migration. In Gambel＇s white-crowned sparrow Zonotrichia leuco- phrys gambelii （GWCS）, photoreceptors for induction of breeding are known to reside in the basal hypothalamus. However, it is unknown whether the sites of photoperiodic reception for vernal migration are the same as those for breeding. Therefore, we hy- pothesized that they may be controlled separately. In this study, we exposed photosensitive GWCSs to low-penetration green light （wavelength at 510 nm） under a regime of 1 lux during the day and 〈0.1 lux at night, and switched the photoperiodic conditions from short day （10 h daytime） to long day （18 h daytime）. The results showed that the experimental birds developed traits associ- ated with vernal migration including mass increase, fat deposition and migratory restlessness behavior when transferred from short day to long day green light cycles, while control birds maintained continuously on short day green light conditions did not ex- press any migration related characteristics. Neither experimental nor control groups showed gonadal recrudescence trader either green light cycles. In support of our hypothesis, we were able to apparently dissociate the photoperiodic responses regulating vernal migra- tion and breeding, which suggests separate mechanisms of photoperiodic time measurement. Such distinct photoperiodic mechanisms may drive the free-tuned temporal arrangement of the two life history stages.展开更多
With the increasing complexity of prospecting objectives,reverse time migration( RTM) has attracted more and more attention due to its outstanding imaging quality. RTM is based on two-way wave equation,so it can avoid...With the increasing complexity of prospecting objectives,reverse time migration( RTM) has attracted more and more attention due to its outstanding imaging quality. RTM is based on two-way wave equation,so it can avoid the limits of angle in traditional one-way wave equation migration,image reverse branch,prism waves and multi-reflected wave precisely and obtain accurate dynamic information. However,the huge demands for storage and computation as well as low frequency noises restrict its wide application. The normalized cross-correlation imaging conditions based on wave field decomposition are derived from traditional cross-correlation imaging condition,and it can eliminate the low-frequency noises effectively and improve the imaging resolution. The practical procedure includes separating source and receiver wave field into one-way components respectively,and conducting cross-correlation imaging condition to the post-separated wave field. In this way,the resolution and precision of the imaging result will be promoted greatly.展开更多
The high-frequency electromagnetic waves of ground-penetrating radar(GPR)attenuate severely when propagated in an underground attenuating medium owing to the influence of resistivity,which remarkably decreases the res...The high-frequency electromagnetic waves of ground-penetrating radar(GPR)attenuate severely when propagated in an underground attenuating medium owing to the influence of resistivity,which remarkably decreases the resolution of reverse time migration(RTM).As an effective high-resolution imaging method,attenuation-compensated RTM(ACRTM)can eff ectively compensate for the energy loss caused by the attenuation related to media absorption under the influence of resistivity.Therefore,constructing an accurate resistivity-media model to compensate for the attenuation of electromagnetic wave energy is crucial for realizing the ACRTM imaging of GPR data.This study proposes a resistivity-constrained ACRTM imaging method for the imaging of GPR data by adding high-density resistivity detection along the GPR survey line and combining it with its resistivity inversion profile.The proposed method uses the inversion result of apparent resistivity data as the GPR RTM-resistivity model for imposing resistivity constraints.Moreover,the hybrid method involving image minimum entropy and RTM is used to estimate the medium velocity at the diff raction position,and combined with the distribution characteristics of the reflection in the GPR profile,a highly accurate velocity model is built to improve the imaging resolution of the ACRTM.The accuracy and eff ectiveness of the proposed method are verified using the ACRTM test of the GPR simulated data of a typical attenuating media model.On this basis,the GPR and apparent resistivity data were observed on a field survey line,and use the GPR resistivity-constrained ACRTM method to image the observed data.A comparison of the proposed method with the conventional ACRTM method shows that the proposed method has better imaging depth,stronger energy,and higher resolution,and the obtained results are more conducive for subsequent data analysis and interpretation.展开更多
基金support this research during the 14th Fiveyear Plan period under contract number 2021QNLM020001the Major Scientific and Technological Projects of Shandong Energy Group under contract number SNKJ2022A06-R23+2 种基金the National Natural Science Foundation of China under contract number 42374164the Funds for Creative Research Groups of China under contract number 41821002the basic theoretical research of seismic wave imaging technology in complex oilfield of Changqing Oilfield Company under contract number 2023e10502.
文摘Amplitude dissipation and phase dispersion occur when seismic waves propagate in attenuated anisotropic media,affecting the quality of migration imaging.To compensate and correct for these effects,the fractional Laplacian pure viscoacoustic wave equation capable of producing stable and noise-free wavefields has been proposed and implemented in the Q-compensated reverse time migration(RTM).In addition,the second-order Taylor series expansion is usually adopted in the hybrid finite-difference/pseudo-spectral(HFDPS)strategy to solve spatially variable fractional Laplacian.However,during forward modeling and Q-compensated RTM,this HFDPS strategy requires 11 and 17 fast Fourier transforms(FFTs)per time step,respectively,leading to computational inefficiency.To improve computational efficiency,we introduce two high-efficiency HFDPS numerical modeling strategies based on asymptotic approximation and algebraic methods.Through the two strategies,the number of FFTs decreased from 11 to 6 and 5 per time step during forward modeling,respectively.Numerical examples demonstrate that wavefields simulated using the new numerical modeling strategies are accurate and highly efficient.Finally,these strategies are employed for implementing high-efficiency and stable Q-compensated RTM techniques in tilted transversely isotropic media,reducing the number of FFTs from 17 to 9 and 8 per time step,respectively,significantly improving computational efficiency.Synthetic data examples illustrate the effectiveness of the proposed Q-compensated RTM scheme in compensating amplitude dissipation and correcting phase distortion.
基金funded by the Young Scientists Fund of the National Natural Science Foundation of China(42304135)the scientific research project of Gansu Coal Geology Bureau(2023-07).
文摘Reasonable field acquisition geometry can not only guide seismic exploration to obtain sufficient geological information of target body,but also reduce acquisition cost to the maximum.In this study,building on conventional ray-based geometry design methods,we incorporate imaging results as a constraint to optimize the geometry design and evaluate its effectiveness.Firstly,the geological model of the target layer is established based on the geological data of the study area and surface seismic data combined with exploration tasks.Then,the ray-tracing method is employed to simulate and assess the proposed geometry design,verifying whether its parameters meet the exploration requirements.Finally,the imaging effect of the designed geometry on the target layer is tested by the cross-well seismic reverse time migration method.This methodology was applied to design the cross-well seismic acquisition geometry for offshore deviated wells in the X Oilfield.The simulation results demonstrate that the imaging-driven geometry design approach effectively guides field operations,enhances the imaging quality of the target layer,and reduces acquisition costs.
基金the Major Science and Technology Project of CNOOC Limited Project‘Exploration Direction and Key Technology Study of Large and Medium-Sized Gas Fields in Haixi Sag,East China Sea’(No.KJZX-2023-0101)the China National Offshore Oil Corporation(CNOOC)during the‘14th Five-Year Plan’(No.KJGG2022-0104)+1 种基金the National Natural Science Foundation of China(Nos.42074138,42206195)the Science and Technology Innovation Project of the LaoShan Laboratory(No.2021WHZZB0703)。
文摘The East China Sea Shelf Basin,especially the Xihu depression in its northeastern region,has long been regarded as an important target area for oil and gas exploration in China.Since the development of the late Cretaceous,this region has experienced complex and multistage tectonic movements,such as extensional faulting and compressional inversion,forming its current unique geological structure pattern.As one of the largest Mesozoic–Cenozoic sedimentary basins offshore China,the East China Sea Shelf Basin not only contains abundant oil and gas resources but also occupies a key strategic position.In addition,it is gradually becoming an important base and strategic core area for China's oil and gas resource development.Although remarkable progress has been achieved in oil and gas exploration,seismic exploration in the Xihu depression still faces many challenges.Especially in deep target layers,weak seismic reflection energy,low signal-to-noise ratio,and poor wave group continuity seriously affect the efficiency and accuracy of oil and gas exploration.In addition,given the considerable influence of folds and faults in the East China Sea Xihu Basin,transverse anisotropy is widespread,which causes difficulty for conventional imaging technologies based on isotropic assumptions to achieve accurate imaging.Therefore,to address these challenges and meet the need for accurate imaging of complex structures in the Xihu depression of the East China Sea Basin,it is essential to advance research on anisotropic imaging technologies.This study responds to that need by integrating the reverse time migration method under a tilted transversely isotropic(TTI)medium model and applying it to the imaging of the Xihu depression.The imaging profiles based on the TTI medium can substantially improve the imaging accuracy and signal-to-noise ratio,effectively improve the imaging quality of the fault systems and mid-deep strata in the area,and provide new technical support and progress momentum for oil and gas exploration in the Xihu depression of the East China Sea.
基金supported by Scientific Research and Technology Development Project of CNPC(2021DJ4002,2022DJ3908).
文摘The numerical dispersion phenomenon in the finite-difference forward modeling simulations of the wave equation significantly affects the imaging accuracy in acoustic reflection logging.This issue is particularly pronounced in the reverse time migration(RTM)method used for shear-wave(S-wave)logging imaging.This not only affects imaging accuracy but also introduces ambiguities in the interpretation of logging results.To address this challenge,this study proposes the use of a least-squares difference coefficient optimization algorithm aiming to suppress the numerical dispersion phenomenon in the RTM of S-wave reflection imaging logging.By optimizing the difference coefficients,the high-precision finite-difference algorithm serves as an effective operator for both forward and backward RTM processes.This approach is instrumental in eliminating migration illusions,which are often caused by numerical dispersion.The effectiveness of this optimized algorithm is demonstrated through numerical results,which indicate that it can achieve more accurate forward imaging results across various conditions,including high-and low-velocity strata,and is effective in both large and small spatial grids.The results of processing real data demonstrate that numerical dispersion optimization effectively reduces migration artifacts and diminishes ambiguities in logging interpretations.This optimization offers crucial technical support to the RTM method,enhancing its capability for accurately modeling and imaging S-wave reflections.
基金funded by the National Natural Science Foundation of China (Grant Nos. 52021005 and 51739007)the Key Research and Development Plan of Shandong Province (Grant No. 2020ZLYS01)。
文摘Seismic forward-prospecting in tunnels is an important step to ensure excavation safety. Nowadays, most advanced imaging techniques in seismic exploration involve calculating the solution of elastic wave equation in a certain coordinate system. However, considering the cylindrical geometry of common tunnel body, Cartesian coordinate system seemingly has limited applicability in tunnel seismic forwardprospecting. To accurately simulate the seismic signal received in tunnels, previous imaging method using decoupled non-conversion elastic wave equation is extended from Cartesian coordinates to cylindrical coordinates. The proposed method preserves the general finite-difference time-domain(FDTD)scheme in Cartesian coordinates, except for a novel wavefield calculation strategy addressing the singularity issue inherited at the cylindrical axis. Moreover, the procedure of cylindrical elastic reverse time migration(CERTM) in tunnels is introduced based on the decoupled non-conversion elastic wavefield. Its imaging effect is further validated via numerical experiments on typical tunnel models. As indicated in the synthetic examples, both the PP-and SS-images could clearly show the geological structure in front of the tunnel face without obvious crosstalk artifacts. Migration imaging using PP-waves can present satisfactory results with higher resolution information supplemented by the SS-images. The potential of applying the proposed method in real-world cases is demonstrated in a water diversion tunnel. In the end, we share our insights regarding the singularity treatment and further improvement of the proposed method.
基金This research was supported by the National Nature Science Foundation of China （No. 41074100） and the Program for NewCentury Excellent Talents in the University of the Ministry of Education of China （No. NCET- 10-0812）.
文摘Prestack reverse time migration （RTM） is an accurate imaging method ofsubsurface media. The viscoacoustic prestack RTM is of practical significance because itconsiders the viscosity of the subsurface media. One of the steps of RTM is solving thewave equation and extrapolating the wave field forward and backward; therefore, solvingaccurately and efficiently the wave equation affects the imaging results and the efficiencyof RTM. In this study, we use the optimal time-space domain dispersion high-order finite-difference （FD） method to solve the viscoacoustic wave equation. Dispersion analysis andnumerical simulations show that the optimal time-space domain FD method is more accurateand suppresses the numerical dispersion. We use hybrid absorbing boundary conditions tohandle the boundary reflection. We also use source-normalized cross-correlation imagingconditions for migration and apply Laplace filtering to remove the low-frequency noise.Numerical modeling suggests that the viscoacoustic wave equation RTM has higher imagingresolution than the acoustic wave equation RTM when the viscosity of the subsurface isconsidered. In addition, for the wave field extrapolation, we use the adaptive variable-lengthFD operator to calculate the spatial derivatives and improve the computational efficiencywithout compromising the accuracy of the numerical solution.
基金supported by the National Natural Science Foundation of China(Grant No.41574099)the National Key Science and Technology Special Projects(grant No.2016ZX05002004-005)
文摘Anisotropy correction is necessary during the processing of converted PS- wave seismic data to achieve accurate structural imaging, reservoir prediction, and fracture detection. To effectively eliminate the adverse effects of S-wave splitting and to improve PS- wave imaging quality, we tested methods for pre-stack migration imaging and anisotropic correction of PS-wave data. We based this on the propagation rules of seismic waves in a horizontal transverse isotropy medium, which is a fractured medium model that reflects likely subsurface conditions in the field. We used the radial （R） and transverse （T） components of PS-wave data to separate the fast and slow S-wave components, after which their propagation moveout was effectively extracted. Meanwhile, corrections for the energies and propagation moveouts of the R and T components were implemented using mathematical rotation. The PS-wave imaging quality was distinctly improved, and we demonstrated the reliability of our methods through numerical simulations. Applying our methods to three-dimensional and three-component seismic field data from the Xinchang-Hexingchang region of the Western Sichuan Depression in China, we obtained high-quality seismic imaging with continuous reflection wave groups, distinct structural features, and specific stratigraphic contact relationships. This study provides an effective and reliable approach for data processing that will improve the exploration of complex, hidden lithologic gas reservoirs.
文摘In this paper, we use the method of pseudo-offset migration （POM） to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into common conversion scatter point （CCSP） gathers directly by POM, which simplifies the conventional processing procedure for converted waves. The POM gather fold and SNR are high, which is favorable for velocity analysis and especially suitable for seismic data with low SNR. We used equivalent anisotropic theory to compute anisotropic parameters. Based on the scattering wave traveltime equation in a VTI medium, the POM pseudo-offset migration in anisotropic media was deduced. By amplitude-preserving POM gather mapping, velocity analysis, stack processing, and so on, the anisotropic migration results were acquired. The forward modeling computation and actual data processing demonstrate the validity of converted wave pre-stack time migration with amplitude-preservation using the anisotropic POM method.
基金sponsored by The National Natural Science Fund(No.41574098)Sinopec Geophysical Key Laboratory Open Fund(No.wtyjy-wx2016-04-2)
文摘Although conventional reverse time migration can be perfectly applied to structural imaging it lacks the capability of enabling detailed delineation of a lithological reservoir due to irregular illumination. To obtain reliable reflectivity of the subsurface it is necessary to solve the imaging problem using inversion. The least-square reverse time migration （LSRTM） （also known as linearized refleetivity inversion） aims to obtain relatively high-resolution amplitude preserving imaging by including the inverse of the Hessian matrix. In practice, the conjugate gradient algorithm is proven to be an efficient iterative method for enabling use of LSRTM. The velocity gradient can be derived from a cross-correlation between observed data and simulated data, making LSRTM independent of wavelet signature and thus more robust in practice. Tests on synthetic and marine data show that LSRTM has good potential for use in reservoir description and four-dimensional （4D） seismic images compared to traditional RTM and Fourier finite difference （FFD） migration. This paper investigates the first order approximation of LSRTM, which is also known as the linear Born approximation. However, for more complex geological structures a higher order approximation should be considered to improve imaging quality.
文摘Engineering seismic exploration aims at shallow imaging which is confused by statics if the surface is uneven. Direct pre-stack depth migration （DPDM） is based on accurate elevations of sources and receivers, by which static correction is completely abandoned before migration and surely the imaging quality is remarkably improved. To obtain some artificial shot gathers, high-order staggered-grid finite-difference （FD） method is adapted to model acoustic wave propagation. Since the shot gathers are always disturbed by regular interferences, the statics still must be applied to supporting the interference elimination by apparent velocity filtering method. Then all the shot gathers should be removed back to their original positions by reverse statics. Finally, they are migrated by pre-stack reverse-time depth migration and imaged. The numerical experiments show that the DPDM can ideally avoid the mistakes caused by statics and increase imaging precision.
文摘To support amplitude variation with offset （AVO） analysis in complex structure areas, we introduce an amplitude-preserving plane-wave prestack time migration approach based on the double-square-root wave equation in media with little lateral velocity variation. In its implementation, a data mapping algorithm is used to obtain offset-plane-wave data sets from the common-midpoint gathers followed by a non-recursive phase-shift solution with amplitude correction to generate common-image gathers in offset-ray-parameter domain and a structural image. Theoretical model tests and a real data example show that our prestack time migration approach is helpful for AVO analysis in complex geological environments.
基金financially co-supported by the National 973 Project of China(Nos.2014CB239006,2011CB202402)the National Natural Science Foundation of China(Nos.41104069,41274124)+1 种基金the Natural Science Foundation of Shandong Province(No.ZR2011DQ016)the Fundamental Research Funds for the Central Universities of China(No.R1401005A)
文摘Elastic migration has been widely paid attention by employing the vector processing of mul- ticomponent seismic data. Ray based elastic Kirchhoff migration has such properties as high flexibility and high efficiency. However, it has failed to solve many problems caused by multipath. On the other hand, elastic reverse-time migration （RTM） based on the two-way wave equation is known to be capable of dealing with these problems, but it is extremely expensive when applied in 3D cases and velocity model building. Based on the elastic Kirchhoff-Helmholtz integral, we calculate deeoupled backward-continued wavefields by introducing elastic Green functions for P- and S-waves, which is expressed by a summation of elastodynamic Gaussian beams. The PP and polarity-corrected PS images are obtained by calculating the correlation between downward and deeoupled backward-continued vector wavefields, where polarity correction is performed by analyzing the relation between the polarization direction of converted PS waves and incident angle on the interface. To a large extent, our method combines the high efficiency of ray-based migration with the high accuracy of wave-equation based reverse-time migration. Application of this method to multicomponent synthetic datasets from the fault model and Marmousi 2 model demonstrates the validity, flexibility and accuracy of the new method.
基金supported by the National Natural Science Foundation of China(No.41574098)National Key Research and Development Program(No.2016YFC06011050)National Oil and Gas Major Project(No.2016ZX05006-002)
文摘With an increasing demand for high-resolution imaging of complex subsurface structures, thin layers and hidden reservoirs, borehole and cross-well seismic migration methods have become important. However, large differences are observed in the frequency bandwidth between the surface, borehole, and cross-well surveys. Thus, variable-grid- based algorithms have been adapted to reverse-time migration. Further, we introduce Lanczos filtering to ensure the stability of wavefield calculations as well as to decrease the artificial reflections that are caused due to the variable grid size. Finally, we observe that the application of this method to surface survey, borehole, and cross-well seismic data suggests improvements in the delineation of minor fractures and steeply dipping faults.
基金This work is sponsored by the Nation's Key Scientific & Technologic Projects (2003 BA613-09).
文摘Pseudo-offset migration （POM） is a new method for prestack time migration of converted waves that improves on equivalent-offset migration （EOM）. The mapping of POM is different than EOM but the purpose of the two methods is to map the input samples to the common conversion scatter point （CCSP） gathers. This paper introduces the principles of the two migration methods and model parameter sensitivity tests for both POM and EOM. At large offset-to-depth ratios the hyperbolic approximation, the three-term approximation and the double square root （DSR） equation are used to NMO-correct the mapped POM gathers to obtain more accurate migration velocities. These equations were derived and calculated by small pseudo offset. POM is then used to image complex structure and prestack time migration.
基金sponsored by the Natural Science Foundation of China(No.41874139)the Natural Science Foundation of China(No.41674124)Jilin Province Foundation for Excellent Youths(No.20190103139JH)
文摘With the development of seismic exploration,passive-source seismic data has attracted increasing attention.Ambient noise passive seismic sources exists widely in nature and industrial production.Passive seismic data is important in logging while drilling(LWD),large-scale structural exploration,etc.In this paper,we proposed a passive multiple reverse time migration imaging(PMRTMI)method based on wavefield decomposition and normalized imaging conditions method.This method differs from seismic interferometry in that it can use raw passive seismic data directly in RTM imaging without reconstruction of virtual active gather,and we use the wavefield decomposition method to eliminate the low frequency noise in RTM.Further,the energy normalized imaging condition is used in full wavefield decomposition,which can not only enhance the image quality of both edge and deep information but also overcome the wrong energy problem caused by uneven distribution of passive sources;furthermore,this method exhibits high efficiency.Finally,numerical examples with the Marmousi model show the effectiveness of the method.
基金supported by the Marine S&T Fund of Shandong Province for Pilot National Laboratory for Marine Science and Technology(Qingdao)(No.2021QNLM020001)the Major Scientific and Technological Projects of Shandong Energy Group(No.SNKJ2022A06-R23)the Major Scientific and Technological Projects of CNPC(No.ZD2019-183-003).
文摘Intrinsic attenuation of the earth causes energy loss and phase distortion in seismic wave propagation.To obtain high-resolution imaging results,these negative effects must be considered during reverse time migration(RTM).We can easily implement attenuation-compensated RTM using the constant Q viscoacoustic wave equation with decoupled amplitude attenuation and phase dispersion terms.However,the nonphysical amplitude-compensation process will inevitably amplify the high-frequency noise in the wavefield in an exponential form,causing the numerical simulation to become unstable.This is due to the fact that the amplitude of the compensation grows exponentially with frequency.In order to achieve stable attenuation-compensated RTM,we modify the analytic expression of the attenuation compensation extrapolation operator and make it only compensate for amplitude loss within the effective frequency band.Based on this modified analytic formula,we then derive an explicit time-space domain attenuation compensation extrapolation operator.Finally,the implementation procedure of stable attenuation-compensated RTM is presented.In addition to being simple to implement,the newly proposed attenuation-compensated extrapolation operator is superior to the conventional low-pass filter in suppressing random noise,which will further improve the imaging resolution.We use two synthetic and one land seismic datasets to verify the stability and effectiveness of the proposed attenuationcompensated RTM in improving imaging resolution in viscous media.
基金supported by National Natural Science Foundation of China(Nos.41904101,41774133)Natural Science Foundation of Shandong Province(ZR2019QD004)+1 种基金Fundamental Research Funds for the Central Universities(No.19CX02010A)the Open Funds of SINOPEC Key Laboratory of Geophysics(Nos.wtyjy-wx2019-01-03,wtyjywx2018-01-06)
文摘In marine seismic exploration,ocean bottom cable technology can record multicomponent seismic data for multiparameter inversion and imaging.This study proposes an elastic multiparameter lease-squares reverse time migration based on the ocean bottom cable technology.Herein,the wavefield continuation operators are mixed equations:the acoustic wave equations are used to calculate seismic wave propagation in the seawater medium,whereas in the solid media below the seabed,the wavefields are obtained by P-and S-wave separated vector elastic wave equations.At the seabed interface,acoustic–elastic coupling control equations are used to combine the two types of equations.P-and S-wave separated elastic migration operators,demigration operators,and gradient equations are derived to realize the elastic least-squares reverse time migration based on the P-and S-wave mode separation.The model tests verify that the proposed method can obtain high-quality images in both the P-and S-velocity components.In comparison with the traditional elastic least-squares reverse time migration method,the proposed method can readily suppress imaging crosstalk noise from multiparameter coupling.
文摘In seasonally breeding birds, the annual cycle of photoperiod is a principal environmental cue for temporal arrange- ment of different life-history stages, such as migration and breeding. In the past, most research has focused on the mechanisms of photoperiodic control of breeding with less attention paid to migration. In Gambel＇s white-crowned sparrow Zonotrichia leuco- phrys gambelii （GWCS）, photoreceptors for induction of breeding are known to reside in the basal hypothalamus. However, it is unknown whether the sites of photoperiodic reception for vernal migration are the same as those for breeding. Therefore, we hy- pothesized that they may be controlled separately. In this study, we exposed photosensitive GWCSs to low-penetration green light （wavelength at 510 nm） under a regime of 1 lux during the day and 〈0.1 lux at night, and switched the photoperiodic conditions from short day （10 h daytime） to long day （18 h daytime）. The results showed that the experimental birds developed traits associ- ated with vernal migration including mass increase, fat deposition and migratory restlessness behavior when transferred from short day to long day green light cycles, while control birds maintained continuously on short day green light conditions did not ex- press any migration related characteristics. Neither experimental nor control groups showed gonadal recrudescence trader either green light cycles. In support of our hypothesis, we were able to apparently dissociate the photoperiodic responses regulating vernal migra- tion and breeding, which suggests separate mechanisms of photoperiodic time measurement. Such distinct photoperiodic mechanisms may drive the free-tuned temporal arrangement of the two life history stages.
文摘With the increasing complexity of prospecting objectives,reverse time migration( RTM) has attracted more and more attention due to its outstanding imaging quality. RTM is based on two-way wave equation,so it can avoid the limits of angle in traditional one-way wave equation migration,image reverse branch,prism waves and multi-reflected wave precisely and obtain accurate dynamic information. However,the huge demands for storage and computation as well as low frequency noises restrict its wide application. The normalized cross-correlation imaging conditions based on wave field decomposition are derived from traditional cross-correlation imaging condition,and it can eliminate the low-frequency noises effectively and improve the imaging resolution. The practical procedure includes separating source and receiver wave field into one-way components respectively,and conducting cross-correlation imaging condition to the post-separated wave field. In this way,the resolution and precision of the imaging result will be promoted greatly.
基金supported by the National Natural Science Foundation of China (No.41604102)the Guangxi Natural Science Foundation project (No.2020GXNSFAA159121).
文摘The high-frequency electromagnetic waves of ground-penetrating radar(GPR)attenuate severely when propagated in an underground attenuating medium owing to the influence of resistivity,which remarkably decreases the resolution of reverse time migration(RTM).As an effective high-resolution imaging method,attenuation-compensated RTM(ACRTM)can eff ectively compensate for the energy loss caused by the attenuation related to media absorption under the influence of resistivity.Therefore,constructing an accurate resistivity-media model to compensate for the attenuation of electromagnetic wave energy is crucial for realizing the ACRTM imaging of GPR data.This study proposes a resistivity-constrained ACRTM imaging method for the imaging of GPR data by adding high-density resistivity detection along the GPR survey line and combining it with its resistivity inversion profile.The proposed method uses the inversion result of apparent resistivity data as the GPR RTM-resistivity model for imposing resistivity constraints.Moreover,the hybrid method involving image minimum entropy and RTM is used to estimate the medium velocity at the diff raction position,and combined with the distribution characteristics of the reflection in the GPR profile,a highly accurate velocity model is built to improve the imaging resolution of the ACRTM.The accuracy and eff ectiveness of the proposed method are verified using the ACRTM test of the GPR simulated data of a typical attenuating media model.On this basis,the GPR and apparent resistivity data were observed on a field survey line,and use the GPR resistivity-constrained ACRTM method to image the observed data.A comparison of the proposed method with the conventional ACRTM method shows that the proposed method has better imaging depth,stronger energy,and higher resolution,and the obtained results are more conducive for subsequent data analysis and interpretation.
