The amplitude versus offset/angle(AVO/AVA)inversion which recovers elastic properties of subsurface media is an essential tool in oil and gas exploration.In general,the exact Zoeppritz equation has a relatively high a...The amplitude versus offset/angle(AVO/AVA)inversion which recovers elastic properties of subsurface media is an essential tool in oil and gas exploration.In general,the exact Zoeppritz equation has a relatively high accuracy in modelling the reflection coefficients.However,amplitude inversion based on it is highly nonlinear,thus,requires nonlinear inversion techniques like the genetic algorithm(GA)which has been widely applied in seismology.The quantum genetic algorithm(QGA)is a variant of the GA that enjoys the advantages of quantum computing,such as qubits and superposition of states.It,however,suffers from limitations in the areas of convergence rate and escaping local minima.To address these shortcomings,in this study,we propose a hybrid quantum genetic algorithm(HQGA)that combines a self-adaptive rotating strategy,and operations of quantum mutation and catastrophe.While the selfadaptive rotating strategy improves the flexibility and efficiency of a quantum rotating gate,the operations of quantum mutation and catastrophe enhance the local and global search abilities,respectively.Using the exact Zoeppritz equation,the HQGA was applied to both synthetic and field seismic data inversion and the results were compared to those of the GA and QGA.A number of the synthetic tests show that the HQGA requires fewer searches to converge to the global solution and the inversion results have generally higher accuracy.The application to field data reveals a good agreement between the inverted parameters and real logs.展开更多
The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to eff...The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to effectively invert these spectral parameters has become the focused area of the complex resistivity inversion. An optimized BP neural network (BPNN) approach based on Quantum Particle Swarm Optimization (QPSO) algorithm was presented, which was able to improve global search ability for complex resistivity multi-parameter nonlinear inversion. In the proposed method, the nonlinear weight adjustment strategy and mutation operator were used to enhance the optimization ability of QPSO algorithm. Implementation of proposed QPSO-BPNN was given, the network had 56 hidden neurons in two hidden layers (the first hidden layer has 46 neurons and the second hidden layer has 10 neurons) and it was trained on 48 datasets and tested on another 5 synthetic datasets. The training and test results show that BP neural network optimized by the QPSO algorithm performs better than the BP neural network without initial optimization on the inversion training and test models, and the mean square error distribution is better. At the same time, a double polarized anomalous bodies model was also used to verify the feasibility and effectiveness of the proposed method, the inversion results show that the QPSO-BP algorithm inversion clearly characterizes the anomalous boundaries and is closer to the values of the parameters.展开更多
The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault par...The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault parameters inversion has nonlinear characteristics, and the gradient-based optimizer(GBO) has the characteristics of fast convergence speed and falling into local optimum hardly. This paper applies GBO algorithm to simulated earthquakes and real LuShan earthquakes in the nonlinear inversion of the Okada model to obtain the source parameters. The simulated earthquake experiment results show that the algorithm is stable, and the seismic source parameters obtained by GBO are slightly closer to the true value than the multi peak particle swarm optimization(MPSO). In the 2013 LuShan earthquake experiment, the root mean square error between the deformation after forwarding of fault parameters obtained by the introduced GBO algorithm and the surface observation deformation was 3.703 mm, slightly better than 3.708 mm calculated by the MPSO. Moreover, the inversion result of GBO algorithm is better than MPSO algorithm in stability. The above results show that the introduced GBO algorithm has a certain practical application value in seismic fault source parameters inversion.展开更多
Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain M...Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo （MCMC） method is proposed in this paper, combining conventional MCMC method based on global optimization with a preconditioned conjugate gradient （PCG） algorithm based on local optimization, so this method does not depend strongly on the initial model. It converges to the global optimum quickly and efficiently on the condition that effi- ciency and stability of inversion are both taken into consid- eration at the same time. The test data verify the feasibility and robustness of the method, and based on this method, we extract the effective pore-fluid bulk modulus, which is applied to reservoir fluid identification and detection, and consequently, a better result has been achieved.展开更多
In order to solve the problems of multi-parameter,multi-extreme and multi-solution in the nonlinear iterative optimization process of Rayleigh wave inversion,the artificial bee colony(ABC)algorithm is selected for glo...In order to solve the problems of multi-parameter,multi-extreme and multi-solution in the nonlinear iterative optimization process of Rayleigh wave inversion,the artificial bee colony(ABC)algorithm is selected for global nonlinear inversion.The global nonlinear inversion method does not rely on a strict initial model and does not need to calculate the derivative of the objective function.The ABC algorithm uses the local optimization behavior of each individual artificial bee to finally highlight the global optimal value in the colony,and the convergence speed is faster.While searching for the global optimal solution,an effective local search can also be performed to ensure the reliability of the inversion results.This paper uses the ABC algorithm to perform Rayleigh wave dispersion inversion on the actual seismic data to obtain a clear undergrounding of shear wave velocity profile and accurately identify the location of the high-velocity interlayer.It is verified that the ABC algorithm used in the inversion of the Rayleigh wave dispersion curve is stable and converges quickly.展开更多
This paper describes an innovative, genetic algorithm based inverse model of nonlinear transducer. In the inverse modeling, using a genetic algorithm, the unknown coefficients of the model are estimated accurately. T...This paper describes an innovative, genetic algorithm based inverse model of nonlinear transducer. In the inverse modeling, using a genetic algorithm, the unknown coefficients of the model are estimated accurately. The simulation results indicate that this technique provides greater flexibility and suitability than the existing methods. It is very easy to modify the nonlinear transducer on line. Thus the method improves the transducer's accuracy. With the help of genetic algorithm (GA), the model coefficients' training are less likely to be trapped in local minima than traditional gradient based search algorithms.展开更多
The conventional Markov chain Monte Carlo （MCMC） method is limited to the selected shape and size of proposal distribution and is not easy to start when the initial proposal distribution is far away from the target ...The conventional Markov chain Monte Carlo （MCMC） method is limited to the selected shape and size of proposal distribution and is not easy to start when the initial proposal distribution is far away from the target distribution. To overcome these drawbacks of the conventional MCMC method, two useful improvements in MCMC method, adaptive Metropolis （AM） algorithm and delayed rejection （DR） algorithm, are attempted to be combined. The AM algorithm aims at adapting the proposal distribution by using the generated estimators, and the DR algorithm aims at enhancing the efficiency of the improved MCMC method. Based on the improved MCMC method, a Bayesian amplitude versus offset （AVO） inversion method on the basis of the exact Zoeppritz equation has been developed, with which the P- and S-wave velocities and the density can be obtained directly, and the uncertainty of AVO inversion results has been estimated as well. The study based on the logging data and the seismic data demonstrates the feasibility and robustness of the method and shows that all three parameters are well retrieved. So the exact Zoeppritz-based nonlinear inversion method by using the improved MCMC is not only suitable for reservoirs with strong-contrast interfaces and longoffset ranges but also it is more stable, accurate and antinoise.展开更多
常规非线性反演方法虽然对初始模型的依赖大为减弱,但局部收敛现象和计算速度慢仍然是瓶颈.本文提出了一种新的反演方法——量子路径积分算法(Quantum Path Integral Algorithm,简称QPIA).该方法引入量子力学的横向场、传播子等概念,并...常规非线性反演方法虽然对初始模型的依赖大为减弱,但局部收敛现象和计算速度慢仍然是瓶颈.本文提出了一种新的反演方法——量子路径积分算法(Quantum Path Integral Algorithm,简称QPIA).该方法引入量子力学的横向场、传播子等概念,并充分利用量子隧穿效应,大大提高反演的效率,具体是通过对反演目标函数的构建,并以Feynman的传播子来构成模型的接收概率来实现.在对一维大地电磁模型和实际数据进行试验后,表明该方法比常规反演方法更能够精确、稳定和快速地逼近真实模型.展开更多
基金supported by the National Natural Science Foundation of China(U19B6003,42122029)the Strategic Cooperation Technology Projects of CNPC and CUPB(ZLZX 202003)partially supported by SEG/WesternGeco Scholarship,SEG Foundation/Chevron Scholarship,and SEG/Norman and Shirley Domenico Scholarship
文摘The amplitude versus offset/angle(AVO/AVA)inversion which recovers elastic properties of subsurface media is an essential tool in oil and gas exploration.In general,the exact Zoeppritz equation has a relatively high accuracy in modelling the reflection coefficients.However,amplitude inversion based on it is highly nonlinear,thus,requires nonlinear inversion techniques like the genetic algorithm(GA)which has been widely applied in seismology.The quantum genetic algorithm(QGA)is a variant of the GA that enjoys the advantages of quantum computing,such as qubits and superposition of states.It,however,suffers from limitations in the areas of convergence rate and escaping local minima.To address these shortcomings,in this study,we propose a hybrid quantum genetic algorithm(HQGA)that combines a self-adaptive rotating strategy,and operations of quantum mutation and catastrophe.While the selfadaptive rotating strategy improves the flexibility and efficiency of a quantum rotating gate,the operations of quantum mutation and catastrophe enhance the local and global search abilities,respectively.Using the exact Zoeppritz equation,the HQGA was applied to both synthetic and field seismic data inversion and the results were compared to those of the GA and QGA.A number of the synthetic tests show that the HQGA requires fewer searches to converge to the global solution and the inversion results have generally higher accuracy.The application to field data reveals a good agreement between the inverted parameters and real logs.
文摘The significant advantage of the complex resistivity method is to reflect the abnormal body through multi-parameters, but its inversion parameters are more than the resistivity tomography method. Therefore, how to effectively invert these spectral parameters has become the focused area of the complex resistivity inversion. An optimized BP neural network (BPNN) approach based on Quantum Particle Swarm Optimization (QPSO) algorithm was presented, which was able to improve global search ability for complex resistivity multi-parameter nonlinear inversion. In the proposed method, the nonlinear weight adjustment strategy and mutation operator were used to enhance the optimization ability of QPSO algorithm. Implementation of proposed QPSO-BPNN was given, the network had 56 hidden neurons in two hidden layers (the first hidden layer has 46 neurons and the second hidden layer has 10 neurons) and it was trained on 48 datasets and tested on another 5 synthetic datasets. The training and test results show that BP neural network optimized by the QPSO algorithm performs better than the BP neural network without initial optimization on the inversion training and test models, and the mean square error distribution is better. At the same time, a double polarized anomalous bodies model was also used to verify the feasibility and effectiveness of the proposed method, the inversion results show that the QPSO-BP algorithm inversion clearly characterizes the anomalous boundaries and is closer to the values of the parameters.
基金the National Natural Science Foundation of China(Nos.42174011and 41874001).
文摘The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault parameters inversion has nonlinear characteristics, and the gradient-based optimizer(GBO) has the characteristics of fast convergence speed and falling into local optimum hardly. This paper applies GBO algorithm to simulated earthquakes and real LuShan earthquakes in the nonlinear inversion of the Okada model to obtain the source parameters. The simulated earthquake experiment results show that the algorithm is stable, and the seismic source parameters obtained by GBO are slightly closer to the true value than the multi peak particle swarm optimization(MPSO). In the 2013 LuShan earthquake experiment, the root mean square error between the deformation after forwarding of fault parameters obtained by the introduced GBO algorithm and the surface observation deformation was 3.703 mm, slightly better than 3.708 mm calculated by the MPSO. Moreover, the inversion result of GBO algorithm is better than MPSO algorithm in stability. The above results show that the introduced GBO algorithm has a certain practical application value in seismic fault source parameters inversion.
基金the sponsorship of the National Basic Research Program of China (973 Program,2013CB228604,2014CB239201)the National Oil and Gas Major Projects of China (2011ZX05014-001-010HZ,2011ZX05014-001-006-XY570) for their funding of this research
文摘Elastic impedance inversion with high efficiency and high stability has become one of the main directions of seismic pre-stack inversion. The nonlinear elastic impedance inversion method based on a fast Markov chain Monte Carlo （MCMC） method is proposed in this paper, combining conventional MCMC method based on global optimization with a preconditioned conjugate gradient （PCG） algorithm based on local optimization, so this method does not depend strongly on the initial model. It converges to the global optimum quickly and efficiently on the condition that effi- ciency and stability of inversion are both taken into consid- eration at the same time. The test data verify the feasibility and robustness of the method, and based on this method, we extract the effective pore-fluid bulk modulus, which is applied to reservoir fluid identification and detection, and consequently, a better result has been achieved.
文摘In order to solve the problems of multi-parameter,multi-extreme and multi-solution in the nonlinear iterative optimization process of Rayleigh wave inversion,the artificial bee colony(ABC)algorithm is selected for global nonlinear inversion.The global nonlinear inversion method does not rely on a strict initial model and does not need to calculate the derivative of the objective function.The ABC algorithm uses the local optimization behavior of each individual artificial bee to finally highlight the global optimal value in the colony,and the convergence speed is faster.While searching for the global optimal solution,an effective local search can also be performed to ensure the reliability of the inversion results.This paper uses the ABC algorithm to perform Rayleigh wave dispersion inversion on the actual seismic data to obtain a clear undergrounding of shear wave velocity profile and accurately identify the location of the high-velocity interlayer.It is verified that the ABC algorithm used in the inversion of the Rayleigh wave dispersion curve is stable and converges quickly.
文摘This paper describes an innovative, genetic algorithm based inverse model of nonlinear transducer. In the inverse modeling, using a genetic algorithm, the unknown coefficients of the model are estimated accurately. The simulation results indicate that this technique provides greater flexibility and suitability than the existing methods. It is very easy to modify the nonlinear transducer on line. Thus the method improves the transducer's accuracy. With the help of genetic algorithm (GA), the model coefficients' training are less likely to be trapped in local minima than traditional gradient based search algorithms.
基金sponsorship of the National Natural Science Foundation of China (41674130, 41404088)the National Basic Research Program of China (973 Program, 2013CB228604, 2014CB239201)+1 种基金the National Oil and Gas Major Projects of China (2016ZX05027004-001, 2016ZX05002005-09HZ)the Fundamental Research Funds for the Central Universities (14CX02113A, 15CX08002A) for their funding in this research
文摘The conventional Markov chain Monte Carlo （MCMC） method is limited to the selected shape and size of proposal distribution and is not easy to start when the initial proposal distribution is far away from the target distribution. To overcome these drawbacks of the conventional MCMC method, two useful improvements in MCMC method, adaptive Metropolis （AM） algorithm and delayed rejection （DR） algorithm, are attempted to be combined. The AM algorithm aims at adapting the proposal distribution by using the generated estimators, and the DR algorithm aims at enhancing the efficiency of the improved MCMC method. Based on the improved MCMC method, a Bayesian amplitude versus offset （AVO） inversion method on the basis of the exact Zoeppritz equation has been developed, with which the P- and S-wave velocities and the density can be obtained directly, and the uncertainty of AVO inversion results has been estimated as well. The study based on the logging data and the seismic data demonstrates the feasibility and robustness of the method and shows that all three parameters are well retrieved. So the exact Zoeppritz-based nonlinear inversion method by using the improved MCMC is not only suitable for reservoirs with strong-contrast interfaces and longoffset ranges but also it is more stable, accurate and antinoise.
文摘常规非线性反演方法虽然对初始模型的依赖大为减弱,但局部收敛现象和计算速度慢仍然是瓶颈.本文提出了一种新的反演方法——量子路径积分算法(Quantum Path Integral Algorithm,简称QPIA).该方法引入量子力学的横向场、传播子等概念,并充分利用量子隧穿效应,大大提高反演的效率,具体是通过对反演目标函数的构建,并以Feynman的传播子来构成模型的接收概率来实现.在对一维大地电磁模型和实际数据进行试验后,表明该方法比常规反演方法更能够精确、稳定和快速地逼近真实模型.
