Gaussian beam prestack depth migration is an accurate imaging method of subsurface media. Prestack depth migration of multicomponent seismic data improves the accuracy of imaging subsurface complex geological structur...Gaussian beam prestack depth migration is an accurate imaging method of subsurface media. Prestack depth migration of multicomponent seismic data improves the accuracy of imaging subsurface complex geological structures. Viscoelastic prestack depth migration is of practical significance because it considers the viscosity of the subsurface media. We use Gaussian beam migration to compensate for the attenuation in multicomponent seismic data. First, we use the Gaussian beam method to simulate the wave propagation in a viscoelastic medium and introduce the complex velocity Q-related and exact viscoelastic Zoeppritz equation. Second, we discuss PP- and PS-wave Gaussian beam prestack depth migration algorithms for common-shot gathers to derive expressions for the attenuation and compensation. The algorithms correct the amplitude attenuation and phase distortion caused by Q, and realize multicomponent Gaussian beam prestack depth migration based on the attenuation compensation and account for the effect of inaccurate Q on migration. Numerical modeling suggests that the imaging resolution of viscoelastic Gaussian beam prestack depth migration is high when the viscosity of the subsurface is considered.展开更多
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.展开更多
The most popular hardware used for parallel depth migration is the PC-Cluster but its application is limited due to large space occupation and high power consumption. In this paper, we introduce a new hardware archite...The most popular hardware used for parallel depth migration is the PC-Cluster but its application is limited due to large space occupation and high power consumption. In this paper, we introduce a new hardware architecture, based on which the finite difference （FD） wavefield-continuation depth migration can be conducted using the Graphics Processing Unit （GPU） as a CPU coprocessor. We demonstrate the program module and three key optimization steps for implementing FD depth migration： memory, thread structure, and instruction optimizations and consider evaluation methods for the amount of optimization. 2D and 3D models are used to test depth migration on the GPU. The tested results show that the depth migration computational efficiency greatly increased using the general-purpose GPU, increasing by at least 25 times compared to the AMD 2.5 GHz CPU.展开更多
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.展开更多
With increased computational power, reverse-time prestack depth migration（RT-PSDM） has become a preferred imaging tool in seismic exploration, yet its use has remained relatively limited in ground-penetrating radar...With increased computational power, reverse-time prestack depth migration（RT-PSDM） has become a preferred imaging tool in seismic exploration, yet its use has remained relatively limited in ground-penetrating radar（GPR） applications. Complex topography alters the wavefield kinematics making for a challenging imaging problem. Model simulations show that topographic variation can substantially distort reflection amplitudes due to irregular wavefield spreading, attenuation anomalies due to irregular path lengths, and focusing and defocusing effects at the surface. The effects are magnified when the topographic variations are on the same order as the depth of investigation––a situation that is often encountered in GPR investigations. Here, I use a full wave-equation RT-PSDM algorithm to image GPR data in the presence of large topographic variability relative to the depth of investigation. The source and receiver wavefields are propagated directly from the topographic surface and this approach inherently corrects for irregular kinematics, spreading and attenuation. The results show that when GPR data are acquired in areas of extreme topography, RT-PSDM can accurately reconstruct reflector geometry as well as reflection amplitude.展开更多
Beamlet sources have strong local and directional character and can easily accomplish local illumination and migration. Besides, they provide better migration results than conventional migration methods. We introduce ...Beamlet sources have strong local and directional character and can easily accomplish local illumination and migration. Besides, they provide better migration results than conventional migration methods. We introduce the basic principles of beamlet prestack depth migration that includes a windowed Fourier transform and frame theory. We explain the Gabor-Daubechies （G-D） frame based on a Gaussian function. Beamlet decomposition provides information on the local space and direction of wavefield. We synthesize the beamlet source and beamlet records in the wavelet domain using both rectangle and Gaussian windows and then extrapolate the synthesized data with a Fourier finite-difference operator. We test the method using the standard Marmousi model. By comparing and analyzing the migration results of single directional beamlet and beamlets with different windows and directions, we demonstrate the validity of the prestack depth migration with Gaussian beamlets method.展开更多
Least-squares migration （LSM） is applied to image subsurface structures and lithology by minimizing the objective function of the observed seismic and reverse-time migration residual data of various underground refl...Least-squares migration （LSM） is applied to image subsurface structures and lithology by minimizing the objective function of the observed seismic and reverse-time migration residual data of various underground reflectivity models. LSM reduces the migration artifacts, enhances the spatial resolution of the migrated images, and yields a more accurate subsurface reflectivity distribution than that of standard migration. The introduction of regularization constraints effectively improves the stability of the least-squares offset. The commonly used regularization terms are based on the L2-norm, which smooths the migration results, e.g., by smearing the reflectivities, while providing stability. However, in exploration geophysics, reflection structures based on velocity and density are generally observed to be discontinuous in depth, illustrating sparse reflectance. To obtain a sparse migration profile, we propose the super-resolution least-squares Kirchhoff prestack depth migration by solving the L0-norm-constrained optimization problem. Additionally, we introduce a two-stage iterative soft and hard thresholding algorithm to retrieve the super-resolution reflectivity distribution. Further, the proposed algorithm is applied to complex synthetic data. Furthermore, the sensitivity of the proposed algorithm to noise and the dominant frequency of the source wavelet was evaluated. Finally, we conclude that the proposed method improves the spatial resolution and achieves impulse-like reflectivity distribution and can be applied to structural interpretations and complex subsurface imaging.展开更多
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.展开更多
The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward c...The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward continues the source and the receiver wavefields simultaneously. The mixed domain algorithm with forward Fourier and inverse Fourier transform is used to construct the double square root equation wavefield extrapolation operator. This operator separates variables in the wave number domain and variables in the space domain. The phase operation is implemented in the wave number domain, whereas the time delay for lateral velocity variation is corrected in the space domain. The migration algorithm is efficient since the seismic data are not computed shot by shot. The data set test of the Marmousi model indicates that the offset-domain migration provides a satisfied seismic migration section on which complex geologic structures are imaged in media with large and rapid lateral velocity variations.展开更多
Prestack depth migration for seismic reflection data is commonly used tool for imaging complex geological structures such as salt domes, faults, thrust belts, and stratigraphic structures. Phase shift plus interpolati...Prestack depth migration for seismic reflection data is commonly used tool for imaging complex geological structures such as salt domes, faults, thrust belts, and stratigraphic structures. Phase shift plus interpolation (PSPI) algorithm is a useful tool to directly solve a wave equation and the results have natural properties of the wave equation. Amplitude and phase characteristics, in particular, are better preserved. The PSPI algorithm is widely used in hydrocarbon exploration because of its simplicity, efficiency, and reduced efforts for computation. However, meaningful depth image of 3D subsurface requires parallel computing to handle heavy computing time and great amount of input data. We implemented a parallelized version of 3D PSPI for prestack depth migration using Open-Multi-Processing (Open MP) library. We verified its performance through applications to 3D SEG/EAGE salt model with a small scale Linux cluster. Phase-shift was performed in the vertical and horizontal directions, respectively, and then interpolated at each node. This gave a single image gather according to shot gather. After summation of each single image gather, we got a 3D stacked image in the depth domain. The numerical model example shows good agree- ment with the original geological model.展开更多
The alternately directional implicit （ADI） scheme is usually used in 3D depth migration. It splits the 3D square-root operator along crossline and inline directions alternately. In this paper, based on the ideal of ...The alternately directional implicit （ADI） scheme is usually used in 3D depth migration. It splits the 3D square-root operator along crossline and inline directions alternately. In this paper, based on the ideal of data line, the four-way splitting schemes and their splitting errors for the finite-difference （FD） method and the hybrid method are investigated. The wavefield extrapolation of four-way splitting scheme is accomplished on a data line and is stable unconditionally. Numerical analysis of splitting errors show that the two-way FD migration have visible numerical anisotropic errors, and that four-way FD migration has much less splitting errors than two-way FD migration has. For the hybrid method, the differences of numerical anisotropic errors between two-way scheme and four-way scheme are small in the case of lower lateral velocity variations. The schemes presented in this paper can be used in 3D post-stack or prestack depth migration. Two numerical calculations of 3D depth migration are completed. One is the four-way FD and hybrid 3D post-stack depth migration for an impulse response, which shows that the anisotropic errors can be eliminated effectively in the cases of constant and variable velocity variations. The other is the 3D shot-profile prestack depth migration for SEG/EAEG benchmark model with two-way hybrid splitting scheme, which presents good imaging results. The Message Passing Interface （MPI） programme based on shot number is adopted.展开更多
Pre-stack depth migration velocity analysis is one of the keys to influencing the imaging quality of pre-stack migration.In this paper we cover a residual curvature velocity analysis method on angle-domain common imag...Pre-stack depth migration velocity analysis is one of the keys to influencing the imaging quality of pre-stack migration.In this paper we cover a residual curvature velocity analysis method on angle-domain common image gathers（ADCIGs） which can depict the relationship between incident angle and migration depth at imaging points and update the migration velocity.Differing from offset-domain common image gathers（ODCIGs）,ADCIGs are not disturbed by the multi-path problem which contributes to imaging artifacts,thus influencing the velocity analysis.On the basis of horizontal layers,we derive the residual depth equation and also propose a velocity analysis workflow for velocity scanning.The tests to synthetic and field data prove the velocity analysis methods adopted in this paper are robust and valid.展开更多
Traditional pre-stack depth migration can only provide subsurface structural information. However, simple structure information is insufficient for petroleum exploration which also needs amplitude information proporti...Traditional pre-stack depth migration can only provide subsurface structural information. However, simple structure information is insufficient for petroleum exploration which also needs amplitude information proportional to reflection coefficients. In recent years, pre-stack depth migration algorithms which preserve amplitudes and based on the one- way wave equation have been developed. Using the method in the shot domain requires a deconvolution imaging condition which produces some instability in areas with complicated structure and dramatic lateral variation in velocity. Depth migration with preserved amplitude based on the angle domain can overcome the instability of the one-way wave migration imaging condition with preserved amplitude. It can also offer provide velocity analysis in the angle domain of common imaging point gathers. In this paper, based on the foundation of the one-way wave continuation operator with preserved amplitude, we realized the preserved amplitude prestack depth migration in the angle domain. Models and real data validate the accuracy of the method.展开更多
An important research topic for prospecting seismology is to provide a fast accurate velocity model from pre-stack depth migration. Aiming at such a problem, we propose a quadratic precision generalized nonlinear glob...An important research topic for prospecting seismology is to provide a fast accurate velocity model from pre-stack depth migration. Aiming at such a problem, we propose a quadratic precision generalized nonlinear global optimization migration velocity inversion. First we discard the assumption that there is a linear relationship between residual depth and residual velocity and propose a velocity model correction equation with quadratic precision which enables the velocity model from each iteration to approach the real model as quickly as possible. Second, we use a generalized nonlinear inversion to get the global optimal velocity perturbation model to all traces. This method can expedite the convergence speed and also can decrease the probability of falling into a local minimum during inversion. The synthetic data and Mamlousi data examples show that our method has a higher precision and needs only a few iterations and consequently enhances the practicability and accuracy of migration velocity analysis （MVA） in complex areas.展开更多
Gaussian beam migration （GBM） is an effec- tive and robust depth seismic imaging method, which overcomes the disadvantage of Kirchhoff migration in imaging multiple arrivals and has no steep-dip limitation of one-wa...Gaussian beam migration （GBM） is an effec- tive and robust depth seismic imaging method, which overcomes the disadvantage of Kirchhoff migration in imaging multiple arrivals and has no steep-dip limitation of one-way wave equation migration. However, its imaging quality depends on the initial beam parameters, which can make the beam width increase and wave-front spread with the propagation of the central ray, resulting in poor migration accuracy at depth, especially for exploration areas with complex geological structures. To address this problem, we present an adaptive focused beam method for shot-domain prestack depth migration. Using the infor- mation of the input smooth velocity field, we first derive an adaptive focused parameter, which makes a seismic beam focused along the whole central ray to enhance the wave- field construction accuracy in both the shallow and deep regions. Then we introduce this parameter into the GBM, which not only improves imaging quality of deep reflectors but also makes the shallow small-scale geological struc- tures well-defined. As well, using the amplitude-preserved extrapolation operator and deconvolution imaging condi- tion, the concept of amplitude-preserved imaging has been included in our method. Typical numerical examples and the field data processing results demonstrate the validity and adaptability of our method.展开更多
Angle-domain common-image gathers（ADCIGs） are the basic data in migration velocity analysis（MVA） and amplitude variation with angle（AVA） analysis. We propose a common-angle gather-generating scheme using Kirchho...Angle-domain common-image gathers（ADCIGs） are the basic data in migration velocity analysis（MVA） and amplitude variation with angle（AVA） analysis. We propose a common-angle gather-generating scheme using Kirchhoff PSDM based on the traveltime gradient field. The scheme includes three major operations：（1） to calculate the traveltime field of the source and the receiver based on the dynamic programming approach;（2） to obtain the refl ection angle according to the traveltime gradient field in the image space; and（3） to generate the ADCIGs during the migration process. Because of the computation approach, the method for generating ADCIGs is superior to conventional ray-based methods. We use the proposed ADCIGs generation method in 3D large-scale seismic data. The key points of the method are the following.（1） We use common-shot datasets for migration,（2） we load traveltimes based on the shot aperture, and（3） we use the MPI and Open Mp memory sharing to decrease the amount of input and output（I/O）. Numerical examples using synthetic data suggest that the ADCIGs improve the quality of the velocity and the effectiveness of the 3D angle-gather generation scheme.展开更多
The Pre-Stack Depth Migration （PSDM） method based on wavefield continuation is the most reliable method for imaging complex structure in the subsurface, although there are large computational costs and poorly adapti...The Pre-Stack Depth Migration （PSDM） method based on wavefield continuation is the most reliable method for imaging complex structure in the subsurface, although there are large computational costs and poorly adaptive geometry. Plane wave shot migration is another method to perform exact wave equation prestack imaging with high computational efficiency and without the migration aperture problem. Moreover, wavefield energy can be compensated at the target zone by controlled illumination. In this paper, plane wave shot PSDM was implemented by the control of the plane down-going wavefield and selection of number and range of the raypaths in order to optimize the imaging effect. In addition, controlled illumination techniques are applied to enhance the imaging precision of interesting areas at different depths. Numerical calculation indicates that plane wave shot imaging is a rapid and efficient method with less computational cost and easy parallel computation compared to the single-square-root operator imaging for common shot gathers and double- square-root operator imaging for common midpoint gathers.展开更多
To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, t...To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.展开更多
Fault block reservoirs are one of the main types of hydrocarbon reservoirs found in offshore basins,and they are widely distributed within the Mesozoic and Cenozoic basins of the northern South China Sea.Conventional ...Fault block reservoirs are one of the main types of hydrocarbon reservoirs found in offshore basins,and they are widely distributed within the Mesozoic and Cenozoic basins of the northern South China Sea.Conventional seismic profiles of complex fault areas often contain obvious fragmentation and distortion of seismic events that is corresponding to geological structure under the fault.This phenomenon is known as a fault shadow;it occurs in relation to rapid changes in velocity near the fault that deviate the ray path of waves,and it seriously affects the ability to determine the geological structure and subsequently evaluate the reserves of fault reservoirs.In the current conventional tomography inversion method,the velocity model is over-smoothed,which results in distortion of the reflection layer under the fault.Based on the velocity tomography inversion of imaging gathers method and the concept of regularization,this paper first introduces the anisotropy Gauss regularization operator.A high-resolution tomography inversion method is then developed,and the fault-controlled geological guidance is constrained.This technology is then applied to a complex fault block reservoir basin in the South China Sea,and the results show that it can significantly solve the problem of fault shadow imaging and determine the geological structures in the target area.The newly developed method thus has very good application prospects.展开更多
Based on perturbation theory, the wave equation extrapolation operator with mixed domains has the ability to deal with lateral velocity variations. It is the image method that has undergone much research in seismology...Based on perturbation theory, the wave equation extrapolation operator with mixed domains has the ability to deal with lateral velocity variations. It is the image method that has undergone much research in seismology. All extrapolation operators face the problem of choosing the reference velocity due to continuation in depth. The wavefield extrapolation operator with a single reference velocity is suitable for media with weak lateral variation. The multi-reference velocity extrapolation operator can cope with severe lateral velocity variations and improve image accuracy. However, the calculation cost is large. We present a self-adaptive approach to automatically determine the number of selected reference velocities according to the complexity of structure and the given velocity threshold value. The approach can be used to construct the SSF, FFD, WXFD, and GSP multi-reference velocity wavefield extrapolation image algorithms. The result of a salt-dome model data test demonstrates that the self-adoptive multi-reference wavefield extrapolation algorithm has the ability to deal with severe lateral velocity variations and can also be used for structure edge detection. The method is flexible and computationally cost-effective.展开更多
基金financially supported by the National Natural Science Foundation of China(No.U1262207)the National Science and Technology Major Project of China(Nos.2011 ZX05023-005-005 and 2011 ZX05019-006)the PetroChina Innovation Foundation(No.2013D-5006-0303)
文摘Gaussian beam prestack depth migration is an accurate imaging method of subsurface media. Prestack depth migration of multicomponent seismic data improves the accuracy of imaging subsurface complex geological structures. Viscoelastic prestack depth migration is of practical significance because it considers the viscosity of the subsurface media. We use Gaussian beam migration to compensate for the attenuation in multicomponent seismic data. First, we use the Gaussian beam method to simulate the wave propagation in a viscoelastic medium and introduce the complex velocity Q-related and exact viscoelastic Zoeppritz equation. Second, we discuss PP- and PS-wave Gaussian beam prestack depth migration algorithms for common-shot gathers to derive expressions for the attenuation and compensation. The algorithms correct the amplitude attenuation and phase distortion caused by Q, and realize multicomponent Gaussian beam prestack depth migration based on the attenuation compensation and account for the effect of inaccurate Q on migration. Numerical modeling suggests that the imaging resolution of viscoelastic Gaussian beam prestack depth migration is high when the viscosity of the subsurface is considered.
基金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 Natural Science Foundation of China (Nos. 41104083 and 40804024) Fundamental Research Funds for the Central Universities (No, 2011YYL022)
文摘The most popular hardware used for parallel depth migration is the PC-Cluster but its application is limited due to large space occupation and high power consumption. In this paper, we introduce a new hardware architecture, based on which the finite difference （FD） wavefield-continuation depth migration can be conducted using the Graphics Processing Unit （GPU） as a CPU coprocessor. We demonstrate the program module and three key optimization steps for implementing FD depth migration： memory, thread structure, and instruction optimizations and consider evaluation methods for the amount of optimization. 2D and 3D models are used to test depth migration on the GPU. The tested results show that the depth migration computational efficiency greatly increased using the general-purpose GPU, increasing by at least 25 times compared to the AMD 2.5 GHz CPU.
文摘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.
基金The Herbette Fondation at the University of Lausanne, Switzerland
文摘With increased computational power, reverse-time prestack depth migration（RT-PSDM） has become a preferred imaging tool in seismic exploration, yet its use has remained relatively limited in ground-penetrating radar（GPR） applications. Complex topography alters the wavefield kinematics making for a challenging imaging problem. Model simulations show that topographic variation can substantially distort reflection amplitudes due to irregular wavefield spreading, attenuation anomalies due to irregular path lengths, and focusing and defocusing effects at the surface. The effects are magnified when the topographic variations are on the same order as the depth of investigation––a situation that is often encountered in GPR investigations. Here, I use a full wave-equation RT-PSDM algorithm to image GPR data in the presence of large topographic variability relative to the depth of investigation. The source and receiver wavefields are propagated directly from the topographic surface and this approach inherently corrects for irregular kinematics, spreading and attenuation. The results show that when GPR data are acquired in areas of extreme topography, RT-PSDM can accurately reconstruct reflector geometry as well as reflection amplitude.
基金This project is sponsored by the National Natural Science Foundation （40474041）, CNPC Young Innovation Fund （04E7040）, the Post-doctoral Research Station of Zhongyuan 0ilfield, Jiangsu 0ilfield, and CNPC Geophysical Laboratories at the China University of Petroleum （East China）.
文摘Beamlet sources have strong local and directional character and can easily accomplish local illumination and migration. Besides, they provide better migration results than conventional migration methods. We introduce the basic principles of beamlet prestack depth migration that includes a windowed Fourier transform and frame theory. We explain the Gabor-Daubechies （G-D） frame based on a Gaussian function. Beamlet decomposition provides information on the local space and direction of wavefield. We synthesize the beamlet source and beamlet records in the wavelet domain using both rectangle and Gaussian windows and then extrapolate the synthesized data with a Fourier finite-difference operator. We test the method using the standard Marmousi model. By comparing and analyzing the migration results of single directional beamlet and beamlets with different windows and directions, we demonstrate the validity of the prestack depth migration with Gaussian beamlets method.
基金supported by the National Natural Science Foundation of China(No.41422403)
文摘Least-squares migration （LSM） is applied to image subsurface structures and lithology by minimizing the objective function of the observed seismic and reverse-time migration residual data of various underground reflectivity models. LSM reduces the migration artifacts, enhances the spatial resolution of the migrated images, and yields a more accurate subsurface reflectivity distribution than that of standard migration. The introduction of regularization constraints effectively improves the stability of the least-squares offset. The commonly used regularization terms are based on the L2-norm, which smooths the migration results, e.g., by smearing the reflectivities, while providing stability. However, in exploration geophysics, reflection structures based on velocity and density are generally observed to be discontinuous in depth, illustrating sparse reflectance. To obtain a sparse migration profile, we propose the super-resolution least-squares Kirchhoff prestack depth migration by solving the L0-norm-constrained optimization problem. Additionally, we introduce a two-stage iterative soft and hard thresholding algorithm to retrieve the super-resolution reflectivity distribution. Further, the proposed algorithm is applied to complex synthetic data. Furthermore, the sensitivity of the proposed algorithm to noise and the dominant frequency of the source wavelet was evaluated. Finally, we conclude that the proposed method improves the spatial resolution and achieves impulse-like reflectivity distribution and can be applied to structural interpretations and complex subsurface imaging.
基金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.
基金This paper is supported by the National Natural Science Foundation of China (No. 40474047)State Key Laboratory of Geological Processes and Mineral Resources (No. GPMR200654)the Focused Subject Program of Beijing (No. XK104910598).
文摘The offset-domain prestack depth migration with optimal separable approximation, based on the double square root equation, is used to image complex media with large and rapid velocity variations. The method downward continues the source and the receiver wavefields simultaneously. The mixed domain algorithm with forward Fourier and inverse Fourier transform is used to construct the double square root equation wavefield extrapolation operator. This operator separates variables in the wave number domain and variables in the space domain. The phase operation is implemented in the wave number domain, whereas the time delay for lateral velocity variation is corrected in the space domain. The migration algorithm is efficient since the seismic data are not computed shot by shot. The data set test of the Marmousi model indicates that the offset-domain migration provides a satisfied seismic migration section on which complex geologic structures are imaged in media with large and rapid lateral velocity variations.
文摘Prestack depth migration for seismic reflection data is commonly used tool for imaging complex geological structures such as salt domes, faults, thrust belts, and stratigraphic structures. Phase shift plus interpolation (PSPI) algorithm is a useful tool to directly solve a wave equation and the results have natural properties of the wave equation. Amplitude and phase characteristics, in particular, are better preserved. The PSPI algorithm is widely used in hydrocarbon exploration because of its simplicity, efficiency, and reduced efforts for computation. However, meaningful depth image of 3D subsurface requires parallel computing to handle heavy computing time and great amount of input data. We implemented a parallelized version of 3D PSPI for prestack depth migration using Open-Multi-Processing (Open MP) library. We verified its performance through applications to 3D SEG/EAGE salt model with a small scale Linux cluster. Phase-shift was performed in the vertical and horizontal directions, respectively, and then interpolated at each node. This gave a single image gather according to shot gather. After summation of each single image gather, we got a 3D stacked image in the depth domain. The numerical model example shows good agree- ment with the original geological model.
基金This research is supported by the Major State Basic Research Program of Peoples＇s Republic of China （No.C1999032803）, the National Key Nature Science Foundation （No.40004003） and ICMSEC Institute Director Foundation.
文摘The alternately directional implicit （ADI） scheme is usually used in 3D depth migration. It splits the 3D square-root operator along crossline and inline directions alternately. In this paper, based on the ideal of data line, the four-way splitting schemes and their splitting errors for the finite-difference （FD） method and the hybrid method are investigated. The wavefield extrapolation of four-way splitting scheme is accomplished on a data line and is stable unconditionally. Numerical analysis of splitting errors show that the two-way FD migration have visible numerical anisotropic errors, and that four-way FD migration has much less splitting errors than two-way FD migration has. For the hybrid method, the differences of numerical anisotropic errors between two-way scheme and four-way scheme are small in the case of lower lateral velocity variations. The schemes presented in this paper can be used in 3D post-stack or prestack depth migration. Two numerical calculations of 3D depth migration are completed. One is the four-way FD and hybrid 3D post-stack depth migration for an impulse response, which shows that the anisotropic errors can be eliminated effectively in the cases of constant and variable velocity variations. The other is the 3D shot-profile prestack depth migration for SEG/EAEG benchmark model with two-way hybrid splitting scheme, which presents good imaging results. The Message Passing Interface （MPI） programme based on shot number is adopted.
基金supported by the National 863 Program (Grant No.2006AA06Z206,Sustained supported)the National Science and Technology Major Project (Grant No.2008ZX05006-004)SinoPec Group Marine Facies Research (Grant No.08370502000410)
文摘Pre-stack depth migration velocity analysis is one of the keys to influencing the imaging quality of pre-stack migration.In this paper we cover a residual curvature velocity analysis method on angle-domain common image gathers（ADCIGs） which can depict the relationship between incident angle and migration depth at imaging points and update the migration velocity.Differing from offset-domain common image gathers（ODCIGs）,ADCIGs are not disturbed by the multi-path problem which contributes to imaging artifacts,thus influencing the velocity analysis.On the basis of horizontal layers,we derive the residual depth equation and also propose a velocity analysis workflow for velocity scanning.The tests to synthetic and field data prove the velocity analysis methods adopted in this paper are robust and valid.
基金supported by the National 863 Program(Grant No.2006AA06Z206)the National 973 Program(Grant No.2007CB209605)CNPC geophysical laboratories and Ph.D innovative funding in China University of Petroleum(East China)
文摘Traditional pre-stack depth migration can only provide subsurface structural information. However, simple structure information is insufficient for petroleum exploration which also needs amplitude information proportional to reflection coefficients. In recent years, pre-stack depth migration algorithms which preserve amplitudes and based on the one- way wave equation have been developed. Using the method in the shot domain requires a deconvolution imaging condition which produces some instability in areas with complicated structure and dramatic lateral variation in velocity. Depth migration with preserved amplitude based on the angle domain can overcome the instability of the one-way wave migration imaging condition with preserved amplitude. It can also offer provide velocity analysis in the angle domain of common imaging point gathers. In this paper, based on the foundation of the one-way wave continuation operator with preserved amplitude, we realized the preserved amplitude prestack depth migration in the angle domain. Models and real data validate the accuracy of the method.
基金This work is supported by National Natural Science Foundation of China （Grant No.40839905）.
文摘An important research topic for prospecting seismology is to provide a fast accurate velocity model from pre-stack depth migration. Aiming at such a problem, we propose a quadratic precision generalized nonlinear global optimization migration velocity inversion. First we discard the assumption that there is a linear relationship between residual depth and residual velocity and propose a velocity model correction equation with quadratic precision which enables the velocity model from each iteration to approach the real model as quickly as possible. Second, we use a generalized nonlinear inversion to get the global optimal velocity perturbation model to all traces. This method can expedite the convergence speed and also can decrease the probability of falling into a local minimum during inversion. The synthetic data and Mamlousi data examples show that our method has a higher precision and needs only a few iterations and consequently enhances the practicability and accuracy of migration velocity analysis （MVA） in complex areas.
文摘Gaussian beam migration （GBM） is an effec- tive and robust depth seismic imaging method, which overcomes the disadvantage of Kirchhoff migration in imaging multiple arrivals and has no steep-dip limitation of one-way wave equation migration. However, its imaging quality depends on the initial beam parameters, which can make the beam width increase and wave-front spread with the propagation of the central ray, resulting in poor migration accuracy at depth, especially for exploration areas with complex geological structures. To address this problem, we present an adaptive focused beam method for shot-domain prestack depth migration. Using the infor- mation of the input smooth velocity field, we first derive an adaptive focused parameter, which makes a seismic beam focused along the whole central ray to enhance the wave- field construction accuracy in both the shallow and deep regions. Then we introduce this parameter into the GBM, which not only improves imaging quality of deep reflectors but also makes the shallow small-scale geological struc- tures well-defined. As well, using the amplitude-preserved extrapolation operator and deconvolution imaging condi- tion, the concept of amplitude-preserved imaging has been included in our method. Typical numerical examples and the field data processing results demonstrate the validity and adaptability of our method.
基金funded by the National Basic Research Program of China(973 Program)(No.2011 CB201002)the National Natural Science Foundation of China(No.41374117)the great and special projects(No.2011ZX05003-003,2011ZX05005-005-008 HZ,and 2011ZX05006-002)
文摘Angle-domain common-image gathers（ADCIGs） are the basic data in migration velocity analysis（MVA） and amplitude variation with angle（AVA） analysis. We propose a common-angle gather-generating scheme using Kirchhoff PSDM based on the traveltime gradient field. The scheme includes three major operations：（1） to calculate the traveltime field of the source and the receiver based on the dynamic programming approach;（2） to obtain the refl ection angle according to the traveltime gradient field in the image space; and（3） to generate the ADCIGs during the migration process. Because of the computation approach, the method for generating ADCIGs is superior to conventional ray-based methods. We use the proposed ADCIGs generation method in 3D large-scale seismic data. The key points of the method are the following.（1） We use common-shot datasets for migration,（2） we load traveltimes based on the shot aperture, and（3） we use the MPI and Open Mp memory sharing to decrease the amount of input and output（I/O）. Numerical examples using synthetic data suggest that the ADCIGs improve the quality of the velocity and the effectiveness of the 3D angle-gather generation scheme.
基金This project is sporspored by Fund item:the National Development and Innovation Committee Program (2005) 2372the National High-tech R&D Program (863 Program) 2006AA06Z241 of ChinaYouth Innovation Fund of CNPC (Program:Prestack Imaging Integral Study for Complex near Surface)
文摘The Pre-Stack Depth Migration （PSDM） method based on wavefield continuation is the most reliable method for imaging complex structure in the subsurface, although there are large computational costs and poorly adaptive geometry. Plane wave shot migration is another method to perform exact wave equation prestack imaging with high computational efficiency and without the migration aperture problem. Moreover, wavefield energy can be compensated at the target zone by controlled illumination. In this paper, plane wave shot PSDM was implemented by the control of the plane down-going wavefield and selection of number and range of the raypaths in order to optimize the imaging effect. In addition, controlled illumination techniques are applied to enhance the imaging precision of interesting areas at different depths. Numerical calculation indicates that plane wave shot imaging is a rapid and efficient method with less computational cost and easy parallel computation compared to the single-square-root operator imaging for common shot gathers and double- square-root operator imaging for common midpoint gathers.
基金National Natural Science Foundation of China (49894190-024) and Geophysical Prospecting Key Laboratory Foundation of China National Petroleum Corporation.
文摘To the most of velocity fields, the traveltimes of the first break that seismic waves propagate along rays can be computed on a 2-D or 3-D numerical grid by finite-difference extrapolation. Under ensuring accuracy, to improve calculating efficiency and adaptability, the calculation method of first-arrival traveltime of finite-difference is de- rived based on any rectangular grid and a local plane wavefront approximation. In addition, head waves and scat- tering waves are properly treated and shadow and caustic zones cannot be encountered, which appear in traditional ray-tracing. The testes of two simple models and the complex Marmousi model show that the method has higher accuracy and adaptability to complex structure with strong vertical and lateral velocity variation, and Kirchhoff prestack depth migration based on this method can basically achieve the position imaging effects of wave equation prestack depth migration in major structures and targets. Because of not taking account of the later arrivals energy, the effect of its amplitude preservation is worse than that by wave equation method, but its computing efficiency is higher than that by total Green′s function method and wave equation method.
基金the National Science and Technology Major Project of China(No.2016ZX05026-002)the National Natural Science Founda-tion of China(Nos.42106072,42074138)+4 种基金the Shandong Provincial Natural Science Foundation(No.ZR2020QD071)the Major Scientific and Technological Innovation Project of Shandong Province(No.2019JZZY010803)the Fundamental Research Funds for the Central Universities(No.201964016)the Shandong Province post-doctoral in-novation projects of special funds(No.201903079)the China Scholarship Council(No.201906335010).
文摘Fault block reservoirs are one of the main types of hydrocarbon reservoirs found in offshore basins,and they are widely distributed within the Mesozoic and Cenozoic basins of the northern South China Sea.Conventional seismic profiles of complex fault areas often contain obvious fragmentation and distortion of seismic events that is corresponding to geological structure under the fault.This phenomenon is known as a fault shadow;it occurs in relation to rapid changes in velocity near the fault that deviate the ray path of waves,and it seriously affects the ability to determine the geological structure and subsequently evaluate the reserves of fault reservoirs.In the current conventional tomography inversion method,the velocity model is over-smoothed,which results in distortion of the reflection layer under the fault.Based on the velocity tomography inversion of imaging gathers method and the concept of regularization,this paper first introduces the anisotropy Gauss regularization operator.A high-resolution tomography inversion method is then developed,and the fault-controlled geological guidance is constrained.This technology is then applied to a complex fault block reservoir basin in the South China Sea,and the results show that it can significantly solve the problem of fault shadow imaging and determine the geological structures in the target area.The newly developed method thus has very good application prospects.
基金The work is financially supported by the "973" project "Large-scale Scientific Calculation and its Use in the Petroleum Industry (G19990328)" by the "863" Project "High Accuracy Seismic Exploration Technologies in the Transition Area (2000AA602018)" also by the Geophysical Research Institute, Shengli Oilfield Ltd. SINOPEC.
文摘Based on perturbation theory, the wave equation extrapolation operator with mixed domains has the ability to deal with lateral velocity variations. It is the image method that has undergone much research in seismology. All extrapolation operators face the problem of choosing the reference velocity due to continuation in depth. The wavefield extrapolation operator with a single reference velocity is suitable for media with weak lateral variation. The multi-reference velocity extrapolation operator can cope with severe lateral velocity variations and improve image accuracy. However, the calculation cost is large. We present a self-adaptive approach to automatically determine the number of selected reference velocities according to the complexity of structure and the given velocity threshold value. The approach can be used to construct the SSF, FFD, WXFD, and GSP multi-reference velocity wavefield extrapolation image algorithms. The result of a salt-dome model data test demonstrates that the self-adoptive multi-reference wavefield extrapolation algorithm has the ability to deal with severe lateral velocity variations and can also be used for structure edge detection. The method is flexible and computationally cost-effective.
