An approximation for the one-way wave operator takes the form of separated space and wave-number variables and makes it possible to use the FFT, which results in a great improvement in the computational efficiency. Fr...An approximation for the one-way wave operator takes the form of separated space and wave-number variables and makes it possible to use the FFT, which results in a great improvement in the computational efficiency. From the function approximation perspective, the OSA method shares the same separable approximation format to the one-way wave operator as other separable approximation methods but it is the only global function approximation among these methods. This leads to a difference in the phase error curve, impulse response, and migration result from other separable approximation methods. The difference is that the OSA method has higher accuracy, and the sensitivity to the velocity variation declines with increasing order.展开更多
When seismic waves propagate through the geological formation,there is a significant loss of energy and a decrease in imaging resolution,because of the viscoacoustic properties of subsurface medium.This profoundly imp...When seismic waves propagate through the geological formation,there is a significant loss of energy and a decrease in imaging resolution,because of the viscoacoustic properties of subsurface medium.This profoundly impacts seismic wavefield propagation,imaging and interpretation.To accurately image the true structure of subsurface medium,the consensus among geophysicists is to no longer treat subsurface medium as ideal homogeneous medium,but rather to incorporate the viscoacoustic properties of subsurface medium.Based on the generalized screen propagator using conventional acoustic wave equation(acoustic GSP),our developed method introduces viscoacoustic compensation strategy,and derives a one-way wave generalized screen propagator based on time-fractional viscoacoustic wave equation(viscoacoustic GSP).In numerical experiments,we conducted tests on two-dimensional multi-layer model and the Marmousi model.When comparing with the acoustic GSP using the acoustic data,we found that the imaging results of the viscoacoustic GSP using the viscoacoustic data showed a significant attenuation compensation effect,and achieved imaging results for both algorithms were essentially consistent.However,the imaging results of acoustic GSP using viscoacoustic data showed significant attenuation effects,especially for deep subsurface imaging.This indicates that we have proposed an effective method to compensate the attenuated seismic wavefield.Our application on a set of real seismic data demonstrated that the imaging performance of our proposed method in local areas surpassed that of the conventional acoustic GSP.This suggests that our proposed method holds practical value and can more accurately image real subsurface structures while enhancing imaging resolution compared with the conventional acoustic GSP.Finally,with respect to computational efficiency,we gathered statistics on running time to compare our proposed method with conventional Q-RTM,and it is evident that our method exhibits higher computational efficiency.In summary,our proposed viscoacoustic GSP method takes into account the true properties of the medium,still achieves migration results comparable to conventional acoustic GSP.展开更多
Various migration methods have been proposed to image high-angle geological structures and media with strong lateral velocity variations; however, the problems of low precision and high computational cost remain unres...Various migration methods have been proposed to image high-angle geological structures and media with strong lateral velocity variations; however, the problems of low precision and high computational cost remain unresolved. To describe the seismic wave propagation in media with lateral velocity variations and to image high-angle structures, we propose the generalized screen propagator based on particle swarm optimization （PSO-GSP）, for the precise fitting of the single-square-root operator. We use the 2D SEG/EAGE salt model to test the proposed PSO-GSP migration method to image the faults beneath the salt dome and compare the results to those of the conventional high-order generalized screen propagator （GSP） migration and split-step Fourier （SSF） migration. Moreover, we use 2D marine data from the South China Sea to show that the PSO-GSP migration can better image strong reflectors than conventional imaging methods.展开更多
The absorption effect of actual subsurface media can weaken wavefield energy, decrease the dominating frequency, and further lead to reduced resolution. In migration, some actions can be taken to compensate for the ab...The absorption effect of actual subsurface media can weaken wavefield energy, decrease the dominating frequency, and further lead to reduced resolution. In migration, some actions can be taken to compensate for the absorption effect and enhance the resolution. In this paper, we derive a one-way wave equation with an attenuation term based on the time- space domain high angle one-way wave equation. A complicated geological model is then designed and synthetic shot gathers are simulated with acoustic wave equations without and with an absorbing term. The derived one-way wave equation is applied to the migration of the synthetic gathers without and with attenuation compensation for the simulated shot gathers. Three migration profiles are obtained. The first and second profiles are from the shot gathers without and with attenuation using the migration method without compensation, the third one is from the shot gathers with attenuation using the migration method with compensation. The first and third profiles are almost the same, and the second profile is different from the others below the absorptive layers. The amplitudes of the interfaces below the absorptive layers are weak because of their absorption. This method is also applied to field data. It is concluded from the migration examples that the migration method discussed in this paper is feasible.展开更多
Walkaway VSP cannot obtain accurate velocity field,as it asymmetrically reflects ray path and provides uneven coverage to underground target,thereby presenting issues related to imaging quality.In this study,we propos...Walkaway VSP cannot obtain accurate velocity field,as it asymmetrically reflects ray path and provides uneven coverage to underground target,thereby presenting issues related to imaging quality.In this study,we propose combining traveltime tomography and prestack depth migration for VSP of an angle-domain walkaway,in a bid to establish accurate two-dimensional and three-dimensional(3 D)velocity models.First,residual curvature was defined to update velocity,and an accurate velocity field was established.To establish a high-precision velocity model,we deduced the relationship between the residual depth and traveltime of common imaging gathers(CIGs)in walkaway VSP.Solving renewal velocity using the least squares method,a four-parameter tomographic inversion equation was derived comprising formation dip angle,incidence angle,residual depth,and sensitivity matrix.In the angle domain,the reflected wave was divided into up-and down-transmitted waves and their traveltimes were calculated.The systematic cumulative method was employed in prestack depth migration of a complex surface.Through prestack depth migration,the offset-domain CIGs were obtained,and dip angle was established by defining the stack section horizon.Runge–Kutta ray tracing was employed to calculate the ray path from the reflection point to the detection point,to determine the incident angle,and to subsequently calculate the ray path from the reflection point to the irregular surface.The offset-domain residual depths were mapped to the angle domain,and a new tomographic equation was established and solved.Application in the double complex area of the Tarim Basin showed the four-parameter tomographic inversion equation derived in this paper to be both correct and practical and that the migration algorithm was able to adapt to the complex surface.展开更多
In this paper, we developed the theory and algorithm of an elastic one-way boundary element method（BEM） and a corresponding hybrid elastic thin-slab propagator for earth media with sharp boundaries between strong co...In this paper, we developed the theory and algorithm of an elastic one-way boundary element method（BEM） and a corresponding hybrid elastic thin-slab propagator for earth media with sharp boundaries between strong contrast media. This approach can takes the advantage of accurate boundary condition of BEM and completely overcomes the weak contrast limitation of the perturbationtheory based one-way operator approach. The one-way BEM is a smooth boundary approximation, which avoids huge matrix operations in exact full BEM. In addition, the one-way BEM can model the primary-only transmitted and reflected waves and therefore is a valuable tool in elastic imaging and inversion. Through numerical tests for some simple models,we proved the validity and efficiency of the proposed method.展开更多
Ru-Shan Wu has made seminal contributions in many research areas in geophysics,such as seismic-wave propagation,scattering,imaging,and inversion.We highlight some of his research in holography imaging,diffraction tomo...Ru-Shan Wu has made seminal contributions in many research areas in geophysics,such as seismic-wave propagation,scattering,imaging,and inversion.We highlight some of his research in holography imaging,diffraction tomography,seismic-wave scattering and its applications to studying Earth’s heterogeneity,oneway wave propagation and one-return wave modeling,beamlet and dreamlet applications,strong non-linear full-waveform inversion,and direct envelop inversion.展开更多
基金sponsored by the National Natural Science Foundation of China (Nos. 40774069 and 40974074)the State Key Program of National Natural Science of China (No. 40830424)the National 973program (No. 007209603)
文摘An approximation for the one-way wave operator takes the form of separated space and wave-number variables and makes it possible to use the FFT, which results in a great improvement in the computational efficiency. From the function approximation perspective, the OSA method shares the same separable approximation format to the one-way wave operator as other separable approximation methods but it is the only global function approximation among these methods. This leads to a difference in the phase error curve, impulse response, and migration result from other separable approximation methods. The difference is that the OSA method has higher accuracy, and the sensitivity to the velocity variation declines with increasing order.
基金financially supported by the National Natural Science Foundation of China (grant Nos.42004103,42374149)Sichuan Science and Technology Program (grant No.2023NSFSC0257)CNPC Innovation Found (2022DQ02-0306)。
文摘When seismic waves propagate through the geological formation,there is a significant loss of energy and a decrease in imaging resolution,because of the viscoacoustic properties of subsurface medium.This profoundly impacts seismic wavefield propagation,imaging and interpretation.To accurately image the true structure of subsurface medium,the consensus among geophysicists is to no longer treat subsurface medium as ideal homogeneous medium,but rather to incorporate the viscoacoustic properties of subsurface medium.Based on the generalized screen propagator using conventional acoustic wave equation(acoustic GSP),our developed method introduces viscoacoustic compensation strategy,and derives a one-way wave generalized screen propagator based on time-fractional viscoacoustic wave equation(viscoacoustic GSP).In numerical experiments,we conducted tests on two-dimensional multi-layer model and the Marmousi model.When comparing with the acoustic GSP using the acoustic data,we found that the imaging results of the viscoacoustic GSP using the viscoacoustic data showed a significant attenuation compensation effect,and achieved imaging results for both algorithms were essentially consistent.However,the imaging results of acoustic GSP using viscoacoustic data showed significant attenuation effects,especially for deep subsurface imaging.This indicates that we have proposed an effective method to compensate the attenuated seismic wavefield.Our application on a set of real seismic data demonstrated that the imaging performance of our proposed method in local areas surpassed that of the conventional acoustic GSP.This suggests that our proposed method holds practical value and can more accurately image real subsurface structures while enhancing imaging resolution compared with the conventional acoustic GSP.Finally,with respect to computational efficiency,we gathered statistics on running time to compare our proposed method with conventional Q-RTM,and it is evident that our method exhibits higher computational efficiency.In summary,our proposed viscoacoustic GSP method takes into account the true properties of the medium,still achieves migration results comparable to conventional acoustic GSP.
基金supported by the 863 Program of China(No.2013AA064201)National Science and Technology Major Project(No.2016ZX05003-003)
文摘Various migration methods have been proposed to image high-angle geological structures and media with strong lateral velocity variations; however, the problems of low precision and high computational cost remain unresolved. To describe the seismic wave propagation in media with lateral velocity variations and to image high-angle structures, we propose the generalized screen propagator based on particle swarm optimization （PSO-GSP）, for the precise fitting of the single-square-root operator. We use the 2D SEG/EAGE salt model to test the proposed PSO-GSP migration method to image the faults beneath the salt dome and compare the results to those of the conventional high-order generalized screen propagator （GSP） migration and split-step Fourier （SSF） migration. Moreover, we use 2D marine data from the South China Sea to show that the PSO-GSP migration can better image strong reflectors than conventional imaging methods.
基金supported in part by the National Natural Science Foundation of China(No.40974069,41174119)the Research of Novel Method and Technology of Geophysical Prospecting,CNPC(No.2011A-3602)the National Major Science and Technology Program(No.2011ZX05010,2011ZX05024)
文摘The absorption effect of actual subsurface media can weaken wavefield energy, decrease the dominating frequency, and further lead to reduced resolution. In migration, some actions can be taken to compensate for the absorption effect and enhance the resolution. In this paper, we derive a one-way wave equation with an attenuation term based on the time- space domain high angle one-way wave equation. A complicated geological model is then designed and synthetic shot gathers are simulated with acoustic wave equations without and with an absorbing term. The derived one-way wave equation is applied to the migration of the synthetic gathers without and with attenuation compensation for the simulated shot gathers. Three migration profiles are obtained. The first and second profiles are from the shot gathers without and with attenuation using the migration method without compensation, the third one is from the shot gathers with attenuation using the migration method with compensation. The first and third profiles are almost the same, and the second profile is different from the others below the absorptive layers. The amplitudes of the interfaces below the absorptive layers are weak because of their absorption. This method is also applied to field data. It is concluded from the migration examples that the migration method discussed in this paper is feasible.
基金supported by the national project "Geophysical Complex Technologies for Reservoirs and Unconventional Gas Reservoirs"(No.2017 ZX05018-004-003)
文摘Walkaway VSP cannot obtain accurate velocity field,as it asymmetrically reflects ray path and provides uneven coverage to underground target,thereby presenting issues related to imaging quality.In this study,we propose combining traveltime tomography and prestack depth migration for VSP of an angle-domain walkaway,in a bid to establish accurate two-dimensional and three-dimensional(3 D)velocity models.First,residual curvature was defined to update velocity,and an accurate velocity field was established.To establish a high-precision velocity model,we deduced the relationship between the residual depth and traveltime of common imaging gathers(CIGs)in walkaway VSP.Solving renewal velocity using the least squares method,a four-parameter tomographic inversion equation was derived comprising formation dip angle,incidence angle,residual depth,and sensitivity matrix.In the angle domain,the reflected wave was divided into up-and down-transmitted waves and their traveltimes were calculated.The systematic cumulative method was employed in prestack depth migration of a complex surface.Through prestack depth migration,the offset-domain CIGs were obtained,and dip angle was established by defining the stack section horizon.Runge–Kutta ray tracing was employed to calculate the ray path from the reflection point to the detection point,to determine the incident angle,and to subsequently calculate the ray path from the reflection point to the irregular surface.The offset-domain residual depths were mapped to the angle domain,and a new tomographic equation was established and solved.Application in the double complex area of the Tarim Basin showed the four-parameter tomographic inversion equation derived in this paper to be both correct and practical and that the migration algorithm was able to adapt to the complex surface.
基金supported by National Scientific Foundation of China with Grant No. 41774067
文摘In this paper, we developed the theory and algorithm of an elastic one-way boundary element method（BEM） and a corresponding hybrid elastic thin-slab propagator for earth media with sharp boundaries between strong contrast media. This approach can takes the advantage of accurate boundary condition of BEM and completely overcomes the weak contrast limitation of the perturbationtheory based one-way operator approach. The one-way BEM is a smooth boundary approximation, which avoids huge matrix operations in exact full BEM. In addition, the one-way BEM can model the primary-only transmitted and reflected waves and therefore is a valuable tool in elastic imaging and inversion. Through numerical tests for some simple models,we proved the validity and efficiency of the proposed method.
文摘Ru-Shan Wu has made seminal contributions in many research areas in geophysics,such as seismic-wave propagation,scattering,imaging,and inversion.We highlight some of his research in holography imaging,diffraction tomography,seismic-wave scattering and its applications to studying Earth’s heterogeneity,oneway wave propagation and one-return wave modeling,beamlet and dreamlet applications,strong non-linear full-waveform inversion,and direct envelop inversion.
