Traditional deep learning methods pursue complex and single network architectures without considering the petrophysical relationship between different elastic parameters.The mathematical and statistical significance o...Traditional deep learning methods pursue complex and single network architectures without considering the petrophysical relationship between different elastic parameters.The mathematical and statistical significance of the inversion results may lead to model overfitting,especially when there are a limited number of well logs in a working area.Multitask learning provides an eff ective approach to addressing this issue.Simultaneously,learning multiple related tasks can improve a model’s generalization ability to a certain extent,thereby enhancing the performance of related tasks with an equal amount of labeled data.In this study,we propose an end-to-end multitask deep learning model that integrates a fully convolutional network and bidirectional gated recurrent unit for intelligent prestack inversion of“seismic data to elastic parameters.”The use of a Bayesian homoscedastic uncertainty-based loss function enables adaptive learning of the weight coeffi cients for diff erent elastic parameter inversion tasks,thereby reducing uncertainty during the inversion process.The proposed method combines the local feature perception of convolutional neural networks with the long-term memory of bidirectional gated recurrent networks.It maintains the rock physics constraint relationships among diff erent elastic parameters during the inversion process,demonstrating a high level of prediction accuracy.Numerical simulations and processing results of real seismic data validate the eff ectiveness and practicality of the proposed method.展开更多
Prestack seismic inversion methods adopt approximations of the Zoeppritz equations to describe the relation between reflection coefficients and P-wave velocity, S-wave velocity, and density. However, the error in thes...Prestack seismic inversion methods adopt approximations of the Zoeppritz equations to describe the relation between reflection coefficients and P-wave velocity, S-wave velocity, and density. However, the error in these approximations increases with increasing angle of incidence and variation of the elastic parameters, which increases the number of inversion solutions and minimizes the inversion accuracy. In this study, we explore a method for solving the reflection coefficients by using the Zoeppritz equations. To increase the accuracy of prestack inversion, the simultaneous inversion of P-wave velocity, S-wave velocity, and density by using prestack large-angle seismic data is proposed based on generalized linear inversion theory. Moreover, we reduce the ill posedness and increase the convergence of prestack inversion by using the regularization constraint damping factor and the conjugate gradient algorithm. The proposed prestack inversion method uses prestack large-angle seismic data to obtain accurate seismic elastic parameters that conform to prestack seismic data and are consistent with logging data from wells.展开更多
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.展开更多
Bottom-simulating reflectors (BSRs) in seismic profile always indicate the bottom of gas hydrate stability zone, but is difficult to determine the distribution and features of gas hydrate sediments (GHS). In this stud...Bottom-simulating reflectors (BSRs) in seismic profile always indicate the bottom of gas hydrate stability zone, but is difficult to determine the distribution and features of gas hydrate sediments (GHS). In this study, based on AVA forward modeling and angle-domain common-image gathers we use prestack AVA parameters consistency inversion in predicting gas hydrate sediments in the Shenhu area at northern slope of South China Sea, and obtain the vertical and lateral features and saturation of GHS.展开更多
The three parameters of P-wave velocity, S-wave velocity, and density have remarkable differences in conventional prestack inversion accuracy, so study of the consistency inversion of the ＂three parameters＂ is very ...The three parameters of P-wave velocity, S-wave velocity, and density have remarkable differences in conventional prestack inversion accuracy, so study of the consistency inversion of the ＂three parameters＂ is very important. In this paper, we present a new inversion algorithm and approach based on the in-depth analysis of the causes in their accuracy differences. With this new method, the inversion accuracy of the three parameters is improved synchronously by reasonable approximations and mutual constraint among the parameters. Theoretical model calculations and actual data applications with this method indicate that the three elastic parameters all have high inversion accuracy and maintain consistency, which also coincides with the theoretical model and actual data. This method has good application prospects.展开更多
Economic shale gas production requires hydraulic fracture stimulation to increase the formation permeability. Hydraulic fracturing strongly depends on geomechanical parameters such as Young＇s modulus and Poisson＇s r...Economic shale gas production requires hydraulic fracture stimulation to increase the formation permeability. Hydraulic fracturing strongly depends on geomechanical parameters such as Young＇s modulus and Poisson＇s ratio. Fracture-prone sweet spots can be predicted by prestack inversion, which is an ill-posed problem; thus, regularization is needed to obtain unique and stable solutions. To characterize gas-bearing shale sedimentary bodies, elastic parameter variations are regarded as an anisotropic Markov random field. Bayesian statistics are adopted for transforming prestack inversion to the maximum posterior probability. Two energy functions for the lateral and vertical directions are used to describe the distribution, and the expectation-maximization algorithm is used to estimate the hyperparameters of the prior probability of elastic parameters. Finally, the inversion yields clear geological boundaries, high vertical resolution, and reasonable lateral continuity using the conjugate gradient method to minimize the objective function. Antinoise and imaging ability of the method were tested using synthetic and real data.展开更多
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.展开更多
In this paper, we use the method of pseudo-offset migration （POM） to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into co...In this paper, we use the method of pseudo-offset migration （POM） to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into common conversion scatter point （CCSP） gathers directly by POM, which simplifies the conventional processing procedure for converted waves. The POM gather fold and SNR are high, which is favorable for velocity analysis and especially suitable for seismic data with low SNR. We used equivalent anisotropic theory to compute anisotropic parameters. Based on the scattering wave traveltime equation in a VTI medium, the POM pseudo-offset migration in anisotropic media was deduced. By amplitude-preserving POM gather mapping, velocity analysis, stack processing, and so on, the anisotropic migration results were acquired. The forward modeling computation and actual data processing demonstrate the validity of converted wave pre-stack time migration with amplitude-preservation using the anisotropic POM method.展开更多
叠前反演是获取地下介质弹性参数的一种重要手段,马尔可夫蒙特卡洛(Markov Chain Monte Carlo,MCMC)算法是叠前反演求解的经典方法。相比于传统的数值优化算法和线性反演方法,MCMC反演算法具备更高的精度,但仍然存在依赖初始模型、计算...叠前反演是获取地下介质弹性参数的一种重要手段,马尔可夫蒙特卡洛(Markov Chain Monte Carlo,MCMC)算法是叠前反演求解的经典方法。相比于传统的数值优化算法和线性反演方法,MCMC反演算法具备更高的精度,但仍然存在依赖初始模型、计算耗时长和不确定性大等问题。为此,对常规MCMC反演算法进行改进,提出基于构造倾角约束的BLI-MCMC叠前随机反演方法。首先,将地质构造倾角加入先验约束信息中,提高反演的采样效率,降低反演结果的不确定性;然后,利用贝叶斯线性反演(Bayesian Linear Inversion,BLI)算法为MCMC反演提供良好的初始模型,并作为迭代起点,缩短马尔科夫链的燃烧时间,从初始模型角度提高反演的效率。模拟数据和实际资料应用结果均表明,改进后的方法能够显著提高反演精度和效率,保持了地下介质较高的横向连续性。该方法可为地下起伏介质的反演提供技术支撑。展开更多
Wavelets are critical to inversion methods. Incorrect phase estimation will affect the objective function and cause convergence to local minima, and thus produce biased or incorrect results. Based on two simple models...Wavelets are critical to inversion methods. Incorrect phase estimation will affect the objective function and cause convergence to local minima, and thus produce biased or incorrect results. Based on two simple models and ignoring all other factors, we studied the variation of the wavelet phase as a function of frequency and its effect on the prestack waveform inversion. Numerical experiments show that an incorrect phase may result in large deviations from the real solution, even if there is a high similarity between the model and real wavelets. The precision of the inversion slightly improves by using the constant-phase rotation; however, the effect of phase inaccuracy is not eliminated, which limits the precision of prestack inversion.展开更多
基金supported by National Key R&D Program of China(2018YFA0702501)National Natural Science Foundation of China (41974140)+1 种基金Science and Technology Management Department,China National Petroleum Corporation(2022DQ0604-01)China National Petroleum Corporation-China University of Petroleum (Beijing) Strategy。
文摘Traditional deep learning methods pursue complex and single network architectures without considering the petrophysical relationship between different elastic parameters.The mathematical and statistical significance of the inversion results may lead to model overfitting,especially when there are a limited number of well logs in a working area.Multitask learning provides an eff ective approach to addressing this issue.Simultaneously,learning multiple related tasks can improve a model’s generalization ability to a certain extent,thereby enhancing the performance of related tasks with an equal amount of labeled data.In this study,we propose an end-to-end multitask deep learning model that integrates a fully convolutional network and bidirectional gated recurrent unit for intelligent prestack inversion of“seismic data to elastic parameters.”The use of a Bayesian homoscedastic uncertainty-based loss function enables adaptive learning of the weight coeffi cients for diff erent elastic parameter inversion tasks,thereby reducing uncertainty during the inversion process.The proposed method combines the local feature perception of convolutional neural networks with the long-term memory of bidirectional gated recurrent networks.It maintains the rock physics constraint relationships among diff erent elastic parameters during the inversion process,demonstrating a high level of prediction accuracy.Numerical simulations and processing results of real seismic data validate the eff ectiveness and practicality of the proposed method.
基金supported by the 973 Program of China(No.2011CB201104 and 2011ZX05009)the National Science and the Technology Major Project(No.2011ZX05006-06)
文摘Prestack seismic inversion methods adopt approximations of the Zoeppritz equations to describe the relation between reflection coefficients and P-wave velocity, S-wave velocity, and density. However, the error in these approximations increases with increasing angle of incidence and variation of the elastic parameters, which increases the number of inversion solutions and minimizes the inversion accuracy. In this study, we explore a method for solving the reflection coefficients by using the Zoeppritz equations. To increase the accuracy of prestack inversion, the simultaneous inversion of P-wave velocity, S-wave velocity, and density by using prestack large-angle seismic data is proposed based on generalized linear inversion theory. Moreover, we reduce the ill posedness and increase the convergence of prestack inversion by using the regularization constraint damping factor and the conjugate gradient algorithm. The proposed prestack inversion method uses prestack large-angle seismic data to obtain accurate seismic elastic parameters that conform to prestack seismic data and are consistent with logging data from wells.
基金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.
文摘Bottom-simulating reflectors (BSRs) in seismic profile always indicate the bottom of gas hydrate stability zone, but is difficult to determine the distribution and features of gas hydrate sediments (GHS). In this study, based on AVA forward modeling and angle-domain common-image gathers we use prestack AVA parameters consistency inversion in predicting gas hydrate sediments in the Shenhu area at northern slope of South China Sea, and obtain the vertical and lateral features and saturation of GHS.
基金sponsored by the National Major Program (No. 2011ZX05006-006)the 973 Program of China (No. 2011CB201104)Technical Research of Elastic Flooding Boundary and Well Network Optimization at the Development Late Stage of Low Permeable Oil Field (No. 2011ZX05009)
文摘The three parameters of P-wave velocity, S-wave velocity, and density have remarkable differences in conventional prestack inversion accuracy, so study of the consistency inversion of the ＂three parameters＂ is very important. In this paper, we present a new inversion algorithm and approach based on the in-depth analysis of the causes in their accuracy differences. With this new method, the inversion accuracy of the three parameters is improved synchronously by reasonable approximations and mutual constraint among the parameters. Theoretical model calculations and actual data applications with this method indicate that the three elastic parameters all have high inversion accuracy and maintain consistency, which also coincides with the theoretical model and actual data. This method has good application prospects.
基金supported by CNPC fundamental research project(No.2014E-3204)
文摘Economic shale gas production requires hydraulic fracture stimulation to increase the formation permeability. Hydraulic fracturing strongly depends on geomechanical parameters such as Young＇s modulus and Poisson＇s ratio. Fracture-prone sweet spots can be predicted by prestack inversion, which is an ill-posed problem; thus, regularization is needed to obtain unique and stable solutions. To characterize gas-bearing shale sedimentary bodies, elastic parameter variations are regarded as an anisotropic Markov random field. Bayesian statistics are adopted for transforming prestack inversion to the maximum posterior probability. Two energy functions for the lateral and vertical directions are used to describe the distribution, and the expectation-maximization algorithm is used to estimate the hyperparameters of the prior probability of elastic parameters. Finally, the inversion yields clear geological boundaries, high vertical resolution, and reasonable lateral continuity using the conjugate gradient method to minimize the objective function. Antinoise and imaging ability of the method were tested using synthetic and real data.
基金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.
文摘In this paper, we use the method of pseudo-offset migration （POM） to complete converted wave pre-stack time migration with amplitude-preservation in an anisotropic medium. The method maps the original traces into common conversion scatter point （CCSP） gathers directly by POM, which simplifies the conventional processing procedure for converted waves. The POM gather fold and SNR are high, which is favorable for velocity analysis and especially suitable for seismic data with low SNR. We used equivalent anisotropic theory to compute anisotropic parameters. Based on the scattering wave traveltime equation in a VTI medium, the POM pseudo-offset migration in anisotropic media was deduced. By amplitude-preserving POM gather mapping, velocity analysis, stack processing, and so on, the anisotropic migration results were acquired. The forward modeling computation and actual data processing demonstrate the validity of converted wave pre-stack time migration with amplitude-preservation using the anisotropic POM method.
文摘叠前反演是获取地下介质弹性参数的一种重要手段,马尔可夫蒙特卡洛(Markov Chain Monte Carlo,MCMC)算法是叠前反演求解的经典方法。相比于传统的数值优化算法和线性反演方法,MCMC反演算法具备更高的精度,但仍然存在依赖初始模型、计算耗时长和不确定性大等问题。为此,对常规MCMC反演算法进行改进,提出基于构造倾角约束的BLI-MCMC叠前随机反演方法。首先,将地质构造倾角加入先验约束信息中,提高反演的采样效率,降低反演结果的不确定性;然后,利用贝叶斯线性反演(Bayesian Linear Inversion,BLI)算法为MCMC反演提供良好的初始模型,并作为迭代起点,缩短马尔科夫链的燃烧时间,从初始模型角度提高反演的效率。模拟数据和实际资料应用结果均表明,改进后的方法能够显著提高反演精度和效率,保持了地下介质较高的横向连续性。该方法可为地下起伏介质的反演提供技术支撑。
基金financially supported by the National Basic Research Program of China(973 Program)(No.2013CB228600)National Natural Science Foundation of China(No.41304108)+1 种基金Science Foundation of China University of Petroleum,Beijing(No.KYJJ2012-05-06)Science Foundation of China,University of Petroleum,Beijing(No.2462013YJRC007)
文摘Wavelets are critical to inversion methods. Incorrect phase estimation will affect the objective function and cause convergence to local minima, and thus produce biased or incorrect results. Based on two simple models and ignoring all other factors, we studied the variation of the wavelet phase as a function of frequency and its effect on the prestack waveform inversion. Numerical experiments show that an incorrect phase may result in large deviations from the real solution, even if there is a high similarity between the model and real wavelets. The precision of the inversion slightly improves by using the constant-phase rotation; however, the effect of phase inaccuracy is not eliminated, which limits the precision of prestack inversion.
