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.展开更多
This paper presents a novel view of the impact of electron collision off-axis positions on the dynamic properties and relativistic nonlinear Thomson inverse scattering of excited electrons within tightly focused, circ...This paper presents a novel view of the impact of electron collision off-axis positions on the dynamic properties and relativistic nonlinear Thomson inverse scattering of excited electrons within tightly focused, circularly polarized laser pulses of varying intensities. We examine the effects of the transverse ponderomotive force, specifically how the deviation angle and speed of electron motion are affected by the initial off-axis position of the electron and the peak amplitude of the laser pulse. When the laser pulse intensity is low, an increase in the electron's initial off-axis distance results in reduced spatial radiation power, improved collimation, super-continuum phenomena generation, red-shifting of the spectrum's harmonic peak, and significant symmetry in the radiation radial direction. However, in contradiction to conventional understandings,when the laser pulse intensity is relatively high, the properties of the relativistic nonlinear Thomson inverse scattering of the electron deviate from the central axis, changing direction in opposition to the aforementioned effects. After reaching a peak, these properties then shift again, aligning with the previous direction. The complex interplay of these effects suggests a greater nuance and intricacy in the relationship between laser pulse intensity, electron position, and scattering properties than previously thought.展开更多
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.展开更多
As an important indicator parameter of fluid identification,fluid factor has always been a concern for scholars.However,when predicting Russell fluid factor or effective pore-fluid bulk modulus,it is necessary to intr...As an important indicator parameter of fluid identification,fluid factor has always been a concern for scholars.However,when predicting Russell fluid factor or effective pore-fluid bulk modulus,it is necessary to introduce a new rock skeleton parameter which is the dry-rock VP/VS ratio squared(DVRS).In the process of fluid factor calculation or inversion,the existing methods take this parameter as a static constant,which has been estimated in advance,and then apply it to the fluid factor calculation and inversion.The fluid identification analysis based on a portion of the Marmousi 2 model and numerical forward modeling test show that,taking the DVRS as a static constant will limit the identification ability of fluid factor and reduce the inversion accuracy.To solve the above problems,we proposed a new method to regard the DVRS as a dynamic variable varying with depth and lithology for the first time,then apply it to fluid factor calculation and inversion.Firstly,the exact Zoeppritz equations are rewritten into a new form containing the fluid factor and DVRS of upper and lower layers.Next,the new equations are applied to the four parameters simultaneous inversion based on the generalized nonlinear inversion(GNI)method.The testing results on a portion of the Marmousi 2 model and field data show that dynamic DVRS can significantly improve the fluid factor identification ability,effectively suppress illusion.Both synthetic and filed data tests also demonstrate that the GNI method based on Bayesian deterministic inversion(BDI)theory can successfully solve the above four parameter simultaneous inversion problem,and taking the dynamic DVRS as a target inversion parameter can effectively improve the inversion accuracy of fluid factor.All these results completely verified the feasibility and effectiveness of the proposed method.展开更多
For the problem of sensor faults and actuator faults in aircraft attitude control,this paper proposes a fault tolerant control(FTC)scheme based on extended state observer(ESO)and nonlinear dynamic inversion(NDI).First...For the problem of sensor faults and actuator faults in aircraft attitude control,this paper proposes a fault tolerant control(FTC)scheme based on extended state observer(ESO)and nonlinear dynamic inversion(NDI).First,two ESOs are designed to estimate sensor faults and actuator faults respectively.Second,the angular rate signal is reconstructed according to the estimation of sensor faults.Third,in angular rate loop,NDI is designed based on reconstruction of angular rate signals and estimation of actuator faults.The FTC scheme proposed in this paper is testified through numerical simulations.The results show that it is feasible and has good fault tolerant ability.展开更多
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.展开更多
In this paper,by means of the classical Lagrange inversion formula,the authors establish a general nonlinear inverse relation as the solution to the problem proposed in the paper[J.Wang,Nonlinear inverse relations for...In this paper,by means of the classical Lagrange inversion formula,the authors establish a general nonlinear inverse relation as the solution to the problem proposed in the paper[J.Wang,Nonlinear inverse relations for the Bell polynomials via the Lagrange inversion formula,J.Integer Seq.,Vol.22(2019),Article 19.3.8].As applications of this inverse relation,the authors not only find a short proof of another nonlinear inverse relation due to Birmajer,et al.(2012),but also set up a few convolution identities concerning the Mina polynomials.展开更多
Strong-scattering inversion or the inverse problem for strong scattering has different physical-mathematical foundations from the weak-scattering case.Seismic inversion based on wave equation for strong scattering can...Strong-scattering inversion or the inverse problem for strong scattering has different physical-mathematical foundations from the weak-scattering case.Seismic inversion based on wave equation for strong scattering cannot be directly solved by Newton’s local optimization method which is based on weak-nonlinear assumption.Here I try to illustrate the connection between the Schr̈odinger inverse scattering(inverse problem for Schr̈odinger equation)by GLM(Gel’fand-Levitan-Marchenko)the-ory and the direct envelope inversion(DEI)using reflection data.The difference between wave equation and Schr̈odinger equation is that the latter has a potential independent of frequency while the former has a frequency-square dependency in the potential.I also point out that the traditional GLM equation for potential inversion can only recover the high-wavenumber components of impedance profile.I propose to use the Schr̈odinger impedance equation for direct impedance inversion and introduce a singular impedance function which also corresponds to a singular potential for the reconstruction of impedance profile,including discontinuities and long-wavelength velocity structure.I will review the GLM theory and its application to impedance inversion including some numerical examples.Then I analyze the recently developed multi-scale direct envelope inversion(MS-DEI)and its connection to the inverse Schr̈odinger scattering.It is conceivable that the combination of strong-scattering inversion(inverse Schr̈odinger scattering)and weak-scattering inversion(local optimization based inversion)may create some inversion methods working for a whole range of inversion problems in geophysical exploration.展开更多
A CFD-based Numerical Virtual Flight(NVF)simulator is presented,which integrates an unsteady flow solver on moving hybrid grids,a Rigid-Body Dynamics(RBD)solver and a module of the Flight Control System(FCS).A techni...A CFD-based Numerical Virtual Flight(NVF)simulator is presented,which integrates an unsteady flow solver on moving hybrid grids,a Rigid-Body Dynamics(RBD)solver and a module of the Flight Control System(FCS).A technique of dynamic hybrid grids is developed to control the active control surfaces with body morphing,with a technique of parallel unstructured dynamic overlapping grids generating proper moving grids over the deflecting control surfaces(e.g.the afterbody rudders of a missile).For the flow/kinematic coupled problems,the 6 Degree-Of-Freedom(DOF)equations are solved by an explicit or implicit method coupled with the URANS CFD solver.The module of the control law is explicitly coupled into the NVF simulator and then improved by the simulation of the pitching maneuver process of a maneuverable missile model.A nonlinear dynamic inversion method is then implemented to design the control law for the pitching process of the maneuverable missile model.Simulations and analysis of the pitching maneuver process are carried out by the NVF simulator to improve the flight control law.Higher control response performance is obtained by adjusting the gain factors and adding an integrator into the control loop.展开更多
A general stochastic algorithm for solving mixed linear and nonlinear problems was introduced in[11].We show in this paper how it can be used to solve the fault inverse problem,where a planar fault in elastic half-spa...A general stochastic algorithm for solving mixed linear and nonlinear problems was introduced in[11].We show in this paper how it can be used to solve the fault inverse problem,where a planar fault in elastic half-space and a slip on that fault have to be reconstructed from noisy surface displacement measurements.With the parameter giving the plane containing the fault denoted by m and the regularization parameter for the linear part of the inverse problem denoted by C,both modeled as random variables,we derive a formula for the posterior marginal of m.Modeling C as a random variable allows to sweep through a wide range of possible values which was shown to be superior to selecting a fixed value[11].We prove that this posterior marginal of m is convergent as the number of measurement points and the dimension of the space for discretizing slips increase.Simply put,our proof only assumes that the regularized discrete error functional for processing measurements relates to an order 1 quadrature rule and that the union of the finite-dimensional spaces for discretizing slips is dense.Our proof relies on trace class operator theory to show that an adequate sequence of determinants is uniformly bounded.We also explain how our proof can be extended to a whole class of inverse problems,as long as some basic requirements are met.Finally,we show numerical simulations that illustrate the numerical convergence of our algorithm.展开更多
文摘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.
基金Project supported by the National Natural Science Foundation of China (Grant Nos.10947170/A05 and 11104291)the Natural Science Fund for Colleges and Universities in Jiangsu Province (Grant No.10KJB140006)+2 种基金the Natural Sciences Foundation of Shanghai (Grant No.11ZR1441300)the Natural Science Foundation of Nanjing University of Posts and Telecommunications (Grant No.NY221098)the Jiangsu Qing Lan Project for their sponsorship。
文摘This paper presents a novel view of the impact of electron collision off-axis positions on the dynamic properties and relativistic nonlinear Thomson inverse scattering of excited electrons within tightly focused, circularly polarized laser pulses of varying intensities. We examine the effects of the transverse ponderomotive force, specifically how the deviation angle and speed of electron motion are affected by the initial off-axis position of the electron and the peak amplitude of the laser pulse. When the laser pulse intensity is low, an increase in the electron's initial off-axis distance results in reduced spatial radiation power, improved collimation, super-continuum phenomena generation, red-shifting of the spectrum's harmonic peak, and significant symmetry in the radiation radial direction. However, in contradiction to conventional understandings,when the laser pulse intensity is relatively high, the properties of the relativistic nonlinear Thomson inverse scattering of the electron deviate from the central axis, changing direction in opposition to the aforementioned effects. After reaching a peak, these properties then shift again, aligning with the previous direction. The complex interplay of these effects suggests a greater nuance and intricacy in the relationship between laser pulse intensity, electron position, and scattering properties than previously thought.
基金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 National Natural Science Foundation of China(41904116,41874156,42074167 and 42204135)the Natural Science Foundation of Hunan Province(2020JJ5168)the China Postdoctoral Science Foundation(2021M703629)for their funding of this research.
文摘As an important indicator parameter of fluid identification,fluid factor has always been a concern for scholars.However,when predicting Russell fluid factor or effective pore-fluid bulk modulus,it is necessary to introduce a new rock skeleton parameter which is the dry-rock VP/VS ratio squared(DVRS).In the process of fluid factor calculation or inversion,the existing methods take this parameter as a static constant,which has been estimated in advance,and then apply it to the fluid factor calculation and inversion.The fluid identification analysis based on a portion of the Marmousi 2 model and numerical forward modeling test show that,taking the DVRS as a static constant will limit the identification ability of fluid factor and reduce the inversion accuracy.To solve the above problems,we proposed a new method to regard the DVRS as a dynamic variable varying with depth and lithology for the first time,then apply it to fluid factor calculation and inversion.Firstly,the exact Zoeppritz equations are rewritten into a new form containing the fluid factor and DVRS of upper and lower layers.Next,the new equations are applied to the four parameters simultaneous inversion based on the generalized nonlinear inversion(GNI)method.The testing results on a portion of the Marmousi 2 model and field data show that dynamic DVRS can significantly improve the fluid factor identification ability,effectively suppress illusion.Both synthetic and filed data tests also demonstrate that the GNI method based on Bayesian deterministic inversion(BDI)theory can successfully solve the above four parameter simultaneous inversion problem,and taking the dynamic DVRS as a target inversion parameter can effectively improve the inversion accuracy of fluid factor.All these results completely verified the feasibility and effectiveness of the proposed method.
基金supported by the Chinese Aviation Science Fund(20160757001)the National Natural Science Foundation of China(10577012)。
文摘For the problem of sensor faults and actuator faults in aircraft attitude control,this paper proposes a fault tolerant control(FTC)scheme based on extended state observer(ESO)and nonlinear dynamic inversion(NDI).First,two ESOs are designed to estimate sensor faults and actuator faults respectively.Second,the angular rate signal is reconstructed according to the estimation of sensor faults.Third,in angular rate loop,NDI is designed based on reconstruction of angular rate signals and estimation of actuator faults.The FTC scheme proposed in this paper is testified through numerical simulations.The results show that it is feasible and has good fault tolerant ability.
文摘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.
基金supported by the National Natural Science Foundation of China under Grant Nos.11971341 and 12001492the Natural Science Foundation of Zhejiang Province under Grant No.LQ20A010004.
文摘In this paper,by means of the classical Lagrange inversion formula,the authors establish a general nonlinear inverse relation as the solution to the problem proposed in the paper[J.Wang,Nonlinear inverse relations for the Bell polynomials via the Lagrange inversion formula,J.Integer Seq.,Vol.22(2019),Article 19.3.8].As applications of this inverse relation,the authors not only find a short proof of another nonlinear inverse relation due to Birmajer,et al.(2012),but also set up a few convolution identities concerning the Mina polynomials.
文摘Strong-scattering inversion or the inverse problem for strong scattering has different physical-mathematical foundations from the weak-scattering case.Seismic inversion based on wave equation for strong scattering cannot be directly solved by Newton’s local optimization method which is based on weak-nonlinear assumption.Here I try to illustrate the connection between the Schr̈odinger inverse scattering(inverse problem for Schr̈odinger equation)by GLM(Gel’fand-Levitan-Marchenko)the-ory and the direct envelope inversion(DEI)using reflection data.The difference between wave equation and Schr̈odinger equation is that the latter has a potential independent of frequency while the former has a frequency-square dependency in the potential.I also point out that the traditional GLM equation for potential inversion can only recover the high-wavenumber components of impedance profile.I propose to use the Schr̈odinger impedance equation for direct impedance inversion and introduce a singular impedance function which also corresponds to a singular potential for the reconstruction of impedance profile,including discontinuities and long-wavelength velocity structure.I will review the GLM theory and its application to impedance inversion including some numerical examples.Then I analyze the recently developed multi-scale direct envelope inversion(MS-DEI)and its connection to the inverse Schr̈odinger scattering.It is conceivable that the combination of strong-scattering inversion(inverse Schr̈odinger scattering)and weak-scattering inversion(local optimization based inversion)may create some inversion methods working for a whole range of inversion problems in geophysical exploration.
基金supported partially by National Key Research and Development Program (No. 2016YFB0200701)National Natural Science Foundation of China (Nos. 11532016 and 11672324)
文摘A CFD-based Numerical Virtual Flight(NVF)simulator is presented,which integrates an unsteady flow solver on moving hybrid grids,a Rigid-Body Dynamics(RBD)solver and a module of the Flight Control System(FCS).A technique of dynamic hybrid grids is developed to control the active control surfaces with body morphing,with a technique of parallel unstructured dynamic overlapping grids generating proper moving grids over the deflecting control surfaces(e.g.the afterbody rudders of a missile).For the flow/kinematic coupled problems,the 6 Degree-Of-Freedom(DOF)equations are solved by an explicit or implicit method coupled with the URANS CFD solver.The module of the control law is explicitly coupled into the NVF simulator and then improved by the simulation of the pitching maneuver process of a maneuverable missile model.A nonlinear dynamic inversion method is then implemented to design the control law for the pitching process of the maneuverable missile model.Simulations and analysis of the pitching maneuver process are carried out by the NVF simulator to improve the flight control law.Higher control response performance is obtained by adjusting the gain factors and adding an integrator into the control loop.
基金This work was supported by Simons Foundation Collaboration Grant[351025]。
文摘A general stochastic algorithm for solving mixed linear and nonlinear problems was introduced in[11].We show in this paper how it can be used to solve the fault inverse problem,where a planar fault in elastic half-space and a slip on that fault have to be reconstructed from noisy surface displacement measurements.With the parameter giving the plane containing the fault denoted by m and the regularization parameter for the linear part of the inverse problem denoted by C,both modeled as random variables,we derive a formula for the posterior marginal of m.Modeling C as a random variable allows to sweep through a wide range of possible values which was shown to be superior to selecting a fixed value[11].We prove that this posterior marginal of m is convergent as the number of measurement points and the dimension of the space for discretizing slips increase.Simply put,our proof only assumes that the regularized discrete error functional for processing measurements relates to an order 1 quadrature rule and that the union of the finite-dimensional spaces for discretizing slips is dense.Our proof relies on trace class operator theory to show that an adequate sequence of determinants is uniformly bounded.We also explain how our proof can be extended to a whole class of inverse problems,as long as some basic requirements are met.Finally,we show numerical simulations that illustrate the numerical convergence of our algorithm.