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.展开更多
Multi-component seismic exploration technology, combining reflected PP- and converted PS-waves, is an effective tool for solving complicated oil and gas exploration problems. The improvement of converted wave resoluti...Multi-component seismic exploration technology, combining reflected PP- and converted PS-waves, is an effective tool for solving complicated oil and gas exploration problems. The improvement of converted wave resolution is one of the key problems. The main factor affecting converted wave resolution is the absorption of seismic waves in overlying strata. In order to remove the effect of absorption on converted waves, inverse Q filtering is used to improve the resolution. In this paper, we present a method to estimate the S-wave Q values from prestack converted wave gathers. Furthermore, we extend a stable and effective poststack inverse Q filtering method to prestack data which uses wave field continuation along the ray path to compensate for attenuation in prestack common shot PP- and PS-waves. The results of theoretical modeling prove that the method of estimating the S-wave Q values has high precision. The results from synthetic and real data prove that the stable inverse Q filtering method can effectively improve the resolution of prestack PP- and PS-waves.展开更多
For the complicated reservoir description of the GD oilfield, P-wave and S-wave elastic impedance inversion was carried out using pre-stack seismic data to accurately identify the lithology of the reservoir. The joint...For the complicated reservoir description of the GD oilfield, P-wave and S-wave elastic impedance inversion was carried out using pre-stack seismic data to accurately identify the lithology of the reservoir. The joint inversion was performed using three or more partial stacks to overcome the singularity of post-stack seismic inversion that can not satisfy the requirements of complex reservoir description and to avoid the instability of the inversion result caused by low signal-noise ratio in the pre-stack gather. The basic theory of prestack elastic impedance inversion is briefly described in this paper and, using real data of the GD oilfield, the key steps of angle gather wavelet extraction, horizon calibration, S-wave velocity prediction, and elastic parameter extraction were analyzed and studied. The comprehensive interpretation of multiple elastic parameters determined from log analysis is a key to improving the effect ofprestack seismic inversion.展开更多
Mathematical geophone （MG） and equal-time stacking （ETS） principles are used to implement seismic prestack forward modeling with irregular surfaces using the oneway acoustic wave-equation. This method receives sei...Mathematical geophone （MG） and equal-time stacking （ETS） principles are used to implement seismic prestack forward modeling with irregular surfaces using the oneway acoustic wave-equation. This method receives seismic primary reflections from the subsurface using a set of virtual MGs. The receivers can be located anywhere on an irregular observing surface. Moreover, the ETS method utilizes the one-way acoustic wave equation to easily and quickly image and extrapolate seismic reflection data. The method is illustrated using high single-noise ratio common shot gathers computed by numerical forward modeling of two simple models, one with a flat surface and one with an irregular surface, and a complex normal fault model. A prestack depth migration method for irregular surface topography was used to reoroduce the normal fault model with high accuracy.展开更多
The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a m...The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a method to calculate the Q-factor based on the prestack Q-factor inversion and the generalized S-transform. The proposed method specifies a standard primary wavelet and calculates the cumulative Q-factors; then, it finds the interlaminar Q-factors using the relation between Q and offset(QVO) and the Dix formula. The proposed method is alternative to methods that calculate interlaminar Q-factors after horizon picking. Because the frequency spectrum of each horizon can be extracted continuously on a 2D time–frequency spectrum, the method is called the continuous spectral ratio slope(CSRS) method. Compared with the other Q-inversion methods, the method offers nearly effortless computations and stability, and has mathematical and physical significance. We use numerical modeling to verify the feasibility of the method and apply it to real data from an oilfield in Ahdeb, Iraq. The results suggest that the resolution and spatial stability of the Q-profile are optimal and contain abundant interlaminar information that is extremely helpful in making lithology and fluid predictions.展开更多
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.展开更多
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.展开更多
Variable-depth streamer seismic data are characterized by low and high frequencies and can be used to obtain high-quality and resolution images of complex subsurface structures. Taking advantage of the frequency range...Variable-depth streamer seismic data are characterized by low and high frequencies and can be used to obtain high-quality and resolution images of complex subsurface structures. Taking advantage of the frequency range in the variable-depth streamer data, we propose the simultaneous inversion of prestack data from variable-offset stack gathers to obtain the P-wave impedance, S-wave impedance, and density. Next, we validate the method by using model and actual variable-depth streamer data from the Huizhou block. The results suggest that the broadband data recorded by variable-depth streamers improve the signal-to-noise ratio and quality of the inversion results and outperform the constant-depth streamer data in delineating the underground stratigraphy.展开更多
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.展开更多
Seismic velocity is important to migration of seismic data, interpretation of lithology and lithofacies as well as prediction of reservoir. The information of shear wave velocity is required to reduce the uncertainty ...Seismic velocity is important to migration of seismic data, interpretation of lithology and lithofacies as well as prediction of reservoir. The information of shear wave velocity is required to reduce the uncertainty for discriminating lithology, identifying fluid type in porous material and calculating gas saturation in reservoir prediction. Based on Zoeppritz equations, a numeral and scanning method was proposed in this paper. Shear wave velocity can be calculated with prestack converted wave data. The effects were demonstrated by inversion of theoretical and real seismic data.展开更多
Organic reef reservoirs in the platform margin of Kaijiang-Liangping trough in Damaoping area, Sichuan Basin are thin in single layer, fast in lateral variation, and have small P-impedance difference from the surround...Organic reef reservoirs in the platform margin of Kaijiang-Liangping trough in Damaoping area, Sichuan Basin are thin in single layer, fast in lateral variation, and have small P-impedance difference from the surrounding rock, it is difficult to identify and predict the reservoirs and fluid properties by conventional post-stack inversion. Through correlation analysis of core test data and logging P-S wave velocity, this work proposed a formula to calculate the shear wave velocity in different porosity ranges, and solved the issue that some wells in the study area have no S-wave data. AVO forward analysis reveals that formation porosity is the main factor affecting the variation of AVO type, the change of water saturation cannot affect the AVO type, but it has an effect on the change range of AVO. Through cross-plotting analysis of elastic parameters, it is found that fluid factor is a parameter sensitive to gas-bearing property of organic reef reservoir in the study area. By comparing results of post-stack impedance inversion, post-stack high frequency attenuation property, pre-stack simultaneous inversion and AVO anomaly analysis of angle gathers, it is found that the gas-bearing prediction of organic reef reservoirs by using fluid factor derived from simultaneous pre-stack inversion had the highest coincidence rate with actual drilling data. At last, according to the characteristics of fluid factor distribution, the favorable gas-bearing area of the organic reef reservoir in Changxing Formation was predicted, and the organic reef trap at the top of Changxing Formation in Block A of Damaoping area was sorted out as the next exploration target.展开更多
基金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.
基金supported by the 863 Program(Grant No.2007AA06Z218)
文摘Multi-component seismic exploration technology, combining reflected PP- and converted PS-waves, is an effective tool for solving complicated oil and gas exploration problems. The improvement of converted wave resolution is one of the key problems. The main factor affecting converted wave resolution is the absorption of seismic waves in overlying strata. In order to remove the effect of absorption on converted waves, inverse Q filtering is used to improve the resolution. In this paper, we present a method to estimate the S-wave Q values from prestack converted wave gathers. Furthermore, we extend a stable and effective poststack inverse Q filtering method to prestack data which uses wave field continuation along the ray path to compensate for attenuation in prestack common shot PP- and PS-waves. The results of theoretical modeling prove that the method of estimating the S-wave Q values has high precision. The results from synthetic and real data prove that the stable inverse Q filtering method can effectively improve the resolution of prestack PP- and PS-waves.
文摘For the complicated reservoir description of the GD oilfield, P-wave and S-wave elastic impedance inversion was carried out using pre-stack seismic data to accurately identify the lithology of the reservoir. The joint inversion was performed using three or more partial stacks to overcome the singularity of post-stack seismic inversion that can not satisfy the requirements of complex reservoir description and to avoid the instability of the inversion result caused by low signal-noise ratio in the pre-stack gather. The basic theory of prestack elastic impedance inversion is briefly described in this paper and, using real data of the GD oilfield, the key steps of angle gather wavelet extraction, horizon calibration, S-wave velocity prediction, and elastic parameter extraction were analyzed and studied. The comprehensive interpretation of multiple elastic parameters determined from log analysis is a key to improving the effect ofprestack seismic inversion.
基金This work was funded by National Natural Science Foundation of China （No. 40474044）.
文摘Mathematical geophone （MG） and equal-time stacking （ETS） principles are used to implement seismic prestack forward modeling with irregular surfaces using the oneway acoustic wave-equation. This method receives seismic primary reflections from the subsurface using a set of virtual MGs. The receivers can be located anywhere on an irregular observing surface. Moreover, the ETS method utilizes the one-way acoustic wave equation to easily and quickly image and extrapolate seismic reflection data. The method is illustrated using high single-noise ratio common shot gathers computed by numerical forward modeling of two simple models, one with a flat surface and one with an irregular surface, and a complex normal fault model. A prestack depth migration method for irregular surface topography was used to reoroduce the normal fault model with high accuracy.
基金supported by The National Key Research and Development Program Plane(No.2017YFC0601505)National Natural Science Foundation(No.41672325)Science&Technology Department of Sichuan Province Technology Project(No.2017GZ0393)
文摘The attenuation factor or quality factor(Q-factor or Q) has been used to measure the energy attenuation of seismic waves propagating in underground media. Many methods are used to estimate the Q-factor. We propose a method to calculate the Q-factor based on the prestack Q-factor inversion and the generalized S-transform. The proposed method specifies a standard primary wavelet and calculates the cumulative Q-factors; then, it finds the interlaminar Q-factors using the relation between Q and offset(QVO) and the Dix formula. The proposed method is alternative to methods that calculate interlaminar Q-factors after horizon picking. Because the frequency spectrum of each horizon can be extracted continuously on a 2D time–frequency spectrum, the method is called the continuous spectral ratio slope(CSRS) method. Compared with the other Q-inversion methods, the method offers nearly effortless computations and stability, and has mathematical and physical significance. We use numerical modeling to verify the feasibility of the method and apply it to real data from an oilfield in Ahdeb, Iraq. The results suggest that the resolution and spatial stability of the Q-profile are optimal and contain abundant interlaminar information that is extremely helpful in making lithology and fluid predictions.
基金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 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.
基金sponsored by Comprehensive Research of CNOOC(China) Limited(No.YXKY–2013–SZ–02)
文摘Variable-depth streamer seismic data are characterized by low and high frequencies and can be used to obtain high-quality and resolution images of complex subsurface structures. Taking advantage of the frequency range in the variable-depth streamer data, we propose the simultaneous inversion of prestack data from variable-offset stack gathers to obtain the P-wave impedance, S-wave impedance, and density. Next, we validate the method by using model and actual variable-depth streamer data from the Huizhou block. The results suggest that the broadband data recorded by variable-depth streamers improve the signal-to-noise ratio and quality of the inversion results and outperform the constant-depth streamer data in delineating the underground stratigraphy.
基金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.
基金National 973 Key Basic Research Development Program (2005CB422104)SINOPEC's Scientific and Techno-logical Development Program(P05063)
文摘Seismic velocity is important to migration of seismic data, interpretation of lithology and lithofacies as well as prediction of reservoir. The information of shear wave velocity is required to reduce the uncertainty for discriminating lithology, identifying fluid type in porous material and calculating gas saturation in reservoir prediction. Based on Zoeppritz equations, a numeral and scanning method was proposed in this paper. Shear wave velocity can be calculated with prestack converted wave data. The effects were demonstrated by inversion of theoretical and real seismic data.
基金Supported by the National Natural Science Foundation of China(41430316)China National Science and Technology Major Project(2017ZX05008-004-008).
文摘Organic reef reservoirs in the platform margin of Kaijiang-Liangping trough in Damaoping area, Sichuan Basin are thin in single layer, fast in lateral variation, and have small P-impedance difference from the surrounding rock, it is difficult to identify and predict the reservoirs and fluid properties by conventional post-stack inversion. Through correlation analysis of core test data and logging P-S wave velocity, this work proposed a formula to calculate the shear wave velocity in different porosity ranges, and solved the issue that some wells in the study area have no S-wave data. AVO forward analysis reveals that formation porosity is the main factor affecting the variation of AVO type, the change of water saturation cannot affect the AVO type, but it has an effect on the change range of AVO. Through cross-plotting analysis of elastic parameters, it is found that fluid factor is a parameter sensitive to gas-bearing property of organic reef reservoir in the study area. By comparing results of post-stack impedance inversion, post-stack high frequency attenuation property, pre-stack simultaneous inversion and AVO anomaly analysis of angle gathers, it is found that the gas-bearing prediction of organic reef reservoirs by using fluid factor derived from simultaneous pre-stack inversion had the highest coincidence rate with actual drilling data. At last, according to the characteristics of fluid factor distribution, the favorable gas-bearing area of the organic reef reservoir in Changxing Formation was predicted, and the organic reef trap at the top of Changxing Formation in Block A of Damaoping area was sorted out as the next exploration target.
