In Rayleigh wave exploration,the inversion of dispersion curves is a crucial step for obtaining subsurface stratigraphic information,characterized by its multi-parameter and multi-extremum nature.Local optimization al...In Rayleigh wave exploration,the inversion of dispersion curves is a crucial step for obtaining subsurface stratigraphic information,characterized by its multi-parameter and multi-extremum nature.Local optimization algorithms used in dispersion curve inversion are highly dependent on the initial model and are prone to being trapped in local optima,while classical global optimization algorithms often suffer from slow convergence and low solution accuracy.To address these issues,this study introduces the Osprey Optimization Algorithm(OOA),known for its strong global search and local exploitation capabilities,into the inversion of dispersion curves to enhance inversion performance.In noiseless theoretical models,the OOA demonstrates excellent inversion accuracy and stability,accurately recovering model parameters.Even in noisy models,OOA maintains robust performance,achieving high inversion precision under high-noise conditions.In multimode dispersion curve tests,OOA effectively handles higher modes due to its efficient global and local search capabilities,and the inversion results show high consistency with theoretical values.Field data from the Wyoming region in the United States and a landfill site in Italy further verify the practical applicability of the OOA.Comprehensive test results indicate that the OOA outperforms the Particle Swarm Optimization(PSO)algorithm,providing a highly accurate and reliable inversion strategy for dispersion curve inversion.展开更多
Combining the adaptive shrinkage genetic algorithm in the feasible region with the imaging of apparent vertical conductance differential, we have inverted the TEM conductive thin layer. The result of the inversion dem...Combining the adaptive shrinkage genetic algorithm in the feasible region with the imaging of apparent vertical conductance differential, we have inverted the TEM conductive thin layer. The result of the inversion demonstrates that by adaptive shrinkage in the feasible region, the calculation speed accelerates and the calculation precision improves. To a certain extent, in this method we surmount the transient electromagnetic sounding equivalence and reduced equivalence scope. Comparison of the inverted result with the forward curve clearly shows that we can image the conductive thin layer.展开更多
For density inversion of gravity anomaly data, once the inversion method is determined, the main factors affecting the inversion result are the inversion parameters and subdivision scheme. A set of reasonable inversio...For density inversion of gravity anomaly data, once the inversion method is determined, the main factors affecting the inversion result are the inversion parameters and subdivision scheme. A set of reasonable inversion parameters and subdivision scheme can, not only improve the inversion process efficiency, but also ensure inversion result accuracy. The gravity inversion method based on correlation searching and the golden section algorithm is an effective potential field inversion method. It can be used to invert 2D and 3D physical properties with potential data observed on flat or rough surfaces. In this paper, we introduce in detail the density inversion principles based on correlation searching and the golden section algorithm. Considering that the gold section algorithm is not globally optimized. we present a heuristic method to ensure the inversion result is globally optimized. With a series of model tests, we systematically compare and analyze the inversion result efficiency and accuracy with different parameters. Based on the model test results, we conclude the selection principles for each inversion parameter with which the inversion accuracy can be obviously improved.展开更多
This paper presents a new algorithm based on the power inversion （PI） and the linearly constrained minimum variance （LCMV）. This algorithm is capable of adjusting the weights of the antenna array in real time to r...This paper presents a new algorithm based on the power inversion （PI） and the linearly constrained minimum variance （LCMV）. This algorithm is capable of adjusting the weights of the antenna array in real time to respond to and improve the global positioning system （GPS） received signals coming from the desired directions and at the same time to highly suppress the jammers coming from the other directions. The simulation is performed for fixed and moving jammers. It indicates that this structure can give deeper nulls, more than 115 dB depths for fixed jammers and more than 94 dB depths for moving jammers.展开更多
The self-potential method is widely used in environmental and engineering geophysics. Four intelligent optimization algorithms are adopted to design the inversion to interpret self-potential data more accurately and e...The self-potential method is widely used in environmental and engineering geophysics. Four intelligent optimization algorithms are adopted to design the inversion to interpret self-potential data more accurately and efficiently: simulated annealing, genetic, particle swarm optimization, and ant colony optimization. Using both noise-free and noise-added synthetic data, it is demonstrated that all four intelligent algorithms can perform self-potential data inversion effectively. During the numerical experiments, the model distribution in search space, the relative errors of model parameters, and the elapsed time are recorded to evaluate the performance of the inversion. The results indicate that all the intelligent algorithms have good precision and tolerance to noise. Particle swarm optimization has the fastest convergence during iteration because of its good balanced searching capability between global and local minimisation.展开更多
To reduce the computational complexity of matrix inversion, which is the majority of processing in many practical applications, two numerically efficient recursive algorithms （called algorithms I and II, respectively...To reduce the computational complexity of matrix inversion, which is the majority of processing in many practical applications, two numerically efficient recursive algorithms （called algorithms I and II, respectively） are presented. Algorithm I is used to calculate the inverse of such a matrix, whose leading principal minors are all nonzero. Algorithm II, whereby, the inverse of an arbitrary nonsingular matrix can be evaluated is derived via improving the algorithm I. The implementation, for algorithm II or I, involves matrix-vector multiplications and vector outer products. These operations are computationally fast and highly parallelizable. MATLAB simulations show that both recursive algorithms are valid.展开更多
At present, near-surface shear wave velocities are mainly calculated through Rayleigh wave dispersion-curve inversions in engineering surface investigations, but the required calculations pose a highly nonlinear globa...At present, near-surface shear wave velocities are mainly calculated through Rayleigh wave dispersion-curve inversions in engineering surface investigations, but the required calculations pose a highly nonlinear global optimization problem. In order to alleviate the risk of falling into a local optimal solution, this paper introduces a new global optimization method, the shuffle frog-leaping algorithm （SFLA）, into the Rayleigh wave dispersion-curve inversion process. SFLA is a swarm-intelligence-based algorithm that simulates a group of frogs searching for food. It uses a few parameters, achieves rapid convergence, and is capability of effective global searching. In order to test the reliability and calculation performance of SFLA, noise-free and noisy synthetic datasets were inverted. We conducted a comparative analysis with other established algorithms using the noise-free dataset, and then tested the ability of SFLA to cope with data noise. Finally, we inverted a real-world example to examine the applicability of SFLA. Results from both synthetic and field data demonstrated the effectiveness of SFLA in the interpretation of Rayleigh wave dispersion curves. We found that SFLA is superior to the established methods in terms of both reliability and computational efficiency, so it offers great potential to improve our ability to solve geophysical inversion problems.展开更多
This paper discusses the inversion of velocity structure and hypocenters location in the Beijing Tianjin Tangshan Zhangjiakou area by genetic algorithm. The hypocenters location of sele...This paper discusses the inversion of velocity structure and hypocenters location in the Beijing Tianjin Tangshan Zhangjiakou area by genetic algorithm. The hypocenters location of selected earthquakes and crustal structure of this area are obtained using the travel time data of local earthquakes acquired by the Telemetered Seismic Network of Northern China. The mean and standard residuals of hypocenter location acquired by this method are much less than those provided by the report of respective earthquakes. The crustal structure of the first and the second layers obtained interpret the outline of the plain and mountain area in the region successfully and the crustal structure of the third layer nearly coincides with the Moho discontinuity obtained by artificial seismic sounding. These show the genetic algorithm is effective to the inversion of hypocenter location and three dimensional velocity structure.展开更多
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.展开更多
Viscoelastic parameters are becoming more important and their inversion algorithms are studied by many researchers. Genetic algorithms are random, self-adaptive, robust, and heuristic with global search and convergenc...Viscoelastic parameters are becoming more important and their inversion algorithms are studied by many researchers. Genetic algorithms are random, self-adaptive, robust, and heuristic with global search and convergence abilities. Based on the direct VSP wave equation, a genetic algorithm （GA） is introduced to determine the viscoelastic parameters. First, the direct wave equation in frequency is expressed as a function of complex velocity and then the complex velocities estimated by GA inversion. Since the phase velocity and Q-factor both are functions of complex velocity, their values can be computed easily. However, there are so many complex velocities that it is difficult to invert them directly. They can be rewritten as a function of Co and C∞ to reduce the number of parameters during the inversion process. Finally, a theoretical model experiment proves that our algorithm is exact and effective.展开更多
With the continuous development of full tensor gradiometer （FTG） measurement techniques, three-dimensional （3D） inversion of FTG data is becoming increasingly used in oil and gas exploration. In the fast processin...With the continuous development of full tensor gradiometer （FTG） measurement techniques, three-dimensional （3D） inversion of FTG data is becoming increasingly used in oil and gas exploration. In the fast processing and interpretation of large-scale high-precision data, the use of the graphics processing unit process unit （GPU） and preconditioning methods are very important in the data inversion. In this paper, an improved preconditioned conjugate gradient algorithm is proposed by combining the symmetric successive over-relaxation （SSOR） technique and the incomplete Choleksy decomposition conjugate gradient algorithm （ICCG）. Since preparing the preconditioner requires extra time, a parallel implement based on GPU is proposed. The improved method is then applied in the inversion of noise- contaminated synthetic data to prove its adaptability in the inversion of 3D FTG data. Results show that the parallel SSOR-ICCG algorithm based on NVIDIA Tesla C2050 GPU achieves a speedup of approximately 25 times that of a serial program using a 2.0 GHz Central Processing Unit （CPU）. Real airbome gravity-gradiometry data from Vinton salt dome （south- west Louisiana, USA） are also considered. Good results are obtained, which verifies the efficiency and feasibility of the proposed parallel method in fast inversion of 3D FTG data.展开更多
We propose a bio-optical inversion model that retrieves the absorption contributions of phytoplankton and colored detrital matter(CDM),as well as the phytoplankton size classes(PSCs),from total minus water absorption ...We propose a bio-optical inversion model that retrieves the absorption contributions of phytoplankton and colored detrital matter(CDM),as well as the phytoplankton size classes(PSCs),from total minus water absorption spectra.The model is based on three-component separation of phytoplankton size structure and a genetic algorithm.The model performance was tested on two independent datasets(the NASA bio-Optical Marine Algorithm Dataset(NOMAD) and the northern South China Sea(NSCS) dataset).The relationships between the estimated and measured values were strongly linear,especially for aCDM(412),and the Root Mean Square Error(RMSE) of the CDM exponential slope(SCDM) was relatively low.Next,the inversion model was directly applied to in-situ total minus water absorption spectra determined by an underwater meter during a cruise in September 2008,to retrieve the phytoplankton size structure in the seawater.By comparing the measured and retrieved chlorophyll a concentrations,we demonstrated that total and size-specific chlorophyll a concentrations could be retrieved by the model with relatively high accuracy.Finally,we applied the bio-optical inversion model to investigate changes in phytoplankton size structure induced by an anti-cyclonic eddy in the NSCS.展开更多
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.展开更多
The practical application of 3D inversion of gravity data requires a lot of computation time and storage space.To solve this problem,we present an integrated optimization algorithm with the following components:(1)tar...The practical application of 3D inversion of gravity data requires a lot of computation time and storage space.To solve this problem,we present an integrated optimization algorithm with the following components:(1)targeting high accuracy in the space domain and fast computation in the wavenumber domain,we design a fast 3D forward algorithm with high precision;and(2)taking advantage of the symmetry of the inversion matrix,the main calculation in gravity conjugate gradient inversion is decomposed into two forward calculations,thus optimizing the computational efficiency of 3D gravity inversion.We verify the calculation accuracy and efficiency of the optimization algorithm by testing various grid-number models through numerical simulation experiments.展开更多
A genetic algorithm of body waveform inversion is presented for better understanding of crustal and upper mantle structures with deep seismic sounding （DSS） waveform data. General reflection and transmission synthet...A genetic algorithm of body waveform inversion is presented for better understanding of crustal and upper mantle structures with deep seismic sounding （DSS） waveform data. General reflection and transmission synthetic seismogram algorithm, which is capable of calculating the response of thin alternating high and low velocity layers, is applied as a solution for forward modeling, and the genetic algorithm is used to find the optimal solution of the inverse problem. Numerical tests suggest that the method has the capability of resolving low-velocity layers, thin alternating high and low velocity layers, and noise suppression. Waveform inversion using P-wave records from Zeku, Xiahe and Lintao shots in the seismic wide-angle reflection/refraction survey along northeastern Qinghai-Xizang （Tibeteau） Plateau has revealed fine structures of the bottom of the upper crust and alternating layers in the middle/lower crust and topmost upper mantle.展开更多
The genetic algorithm is useful for solving an inversion of complex nonlinear geophysical equations. The multi-point search of the genetic algorithm makes it easier to find a globally optimal solution and avoid fall...The genetic algorithm is useful for solving an inversion of complex nonlinear geophysical equations. The multi-point search of the genetic algorithm makes it easier to find a globally optimal solution and avoid falling into a local extremum. The search efficiency of the genetic algorithm is a key to producing successful solutions in a huge multi-parameter model space. The encoding mechanism of the genetic algorithm affects the searching processes in the evolution. Not all genetic operations perform perfectly in a search under either a binary or decimal encoding system. As such, a standard genetic algorithm (SGA) is sometimes unable to resolve an optimization problem such as a simple geophysical inversion. With the binary encoding system the operation of the crossover may produce more new individuals. The decimal encoding system, on the other hand, makes the mutation generate more new genes. This paper discusses approaches of exploiting the search potentials of genetic operations with different encoding systems and presents a hybrid-encoding mechanism for the genetic algorithm. This is referred to as the hybrid-encoding genetic algorithm (HEGA). The method is based on the routine in which the mutation operation is executed in decimal code and other operations in binary code. HEGA guarantees the birth of better genes by mutation processing with a high probability, so that it is beneficial for resolving the inversions of complicated problems. Synthetic and real-world examples demonstrate the advantages of using HEGA in the inversion of potential-field data.展开更多
The firework algorithm（FWA） is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model...The firework algorithm（FWA） is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model show that the FWA performs well in both solution quality and efficiency. We apply the FWA in this study to crustal velocity structure inversion using regional seismic waveform data of central Gansu on the northeastern margin of the Qinghai-Tibet plateau. Seismograms recorded from the moment magnitude（MW） 5.4 Minxian earthquake enable obtaining an average crustal velocity model for this region. We initially carried out a series of FWA robustness tests in regional waveform inversion at the same earthquake and station positions across the study region,inverting two velocity structure models, with and without a low-velocity crustal layer; the accuracy of our average inversion results and their standard deviations reveal the advantages of the FWA for the inversion of regional seismic waveforms. We applied the FWA across our study area using three component waveform data recorded by nine broadband permanent seismic stations with epicentral distances ranging between 146 and 437 km. These inversion results show that the average thickness of the crust in this region is 46.75 km, while thicknesses of the sedimentary layer, and the upper, middle, and lower crust are 3.15,15.69, 13.08, and 14.83 km, respectively. Results also show that the P-wave velocities of these layers and the upper mantle are 4.47, 6.07, 6.12, 6.87, and 8.18 km/s,respectively.展开更多
The Film Copy Deliverer Problem (FDP),much more difficult than TSP,is a new problem in the combination optimization. In this paper,a new algorithm is introduced. First,the FDP is converted to TSP. Then an evolutionary...The Film Copy Deliverer Problem (FDP),much more difficult than TSP,is a new problem in the combination optimization. In this paper,a new algorithm is introduced. First,the FDP is converted to TSP. Then an evolutionary algorithm based on inversion operator is adopted. Compared with other Genetic Algorithms,it is not only more simple and easier to realize but also faster and more accurate.展开更多
The velocity distribution of layers from surface wave dispersion curve is a severely nonlinear program. Base on the Metropolis rule,we improved the simulated annealing algorithm to simultaneously inverse the velocitie...The velocity distribution of layers from surface wave dispersion curve is a severely nonlinear program. Base on the Metropolis rule,we improved the simulated annealing algorithm to simultaneously inverse the velocities and thicknesses using the dispersion data and identified the Moho and the bottom of lithosphere. The application to the numerical examples with 5% noise shows the velocity RMS is 1. 56% between the non-linear results and the original models when the condition of selecting method for temperature parameters and initial temperature are satisfied. Using the pure dispersions of Rayleigh wave,the nonlinear inversion has been carried out for S-wave velocities and thicknesses of the vertical profile crossing the Indian Plate,the Qinghai-Tibetan Plateau,and the Tarim Basin. It indicated that the crustal thickness is about 70 km in the Qiangtang block,while in the hinterland of the Qinghai-Tibetan Plateau the lithosphere is relatively thin(~ 130 km)from the velocity values and their offsets.展开更多
The traditional genetic algorithm(GA)has unstable inversion results and is easy to fall into the local optimum when inverting fault parameters.Therefore,this article considers the combination of GA with other non-line...The traditional genetic algorithm(GA)has unstable inversion results and is easy to fall into the local optimum when inverting fault parameters.Therefore,this article considers the combination of GA with other non-linear algorithms in order to improve the inversion precision of GA.This paper proposes a genetic Nelder-Mead neural network algorithm(GNMNNA).This algorithm uses a neural network algorithm(NNA)to optimize the global search ability of GA.At the same time,the simplex algorithm is used to optimize the local search capability of the GA.Through numerical examples,the stability of the inversion algorithm under different strategies is explored.The experimental results show that the proposed GNMNNA has stronger inversion stability and higher precision compared with the existing algorithms.The effectiveness of GNMNNA is verified by the BodrumeKos earthquake and Monte Cristo Range earthquake.The experimental results show that GNMNNA is superior to GA and NNA in both inversion precision and computational stability.Therefore,GNMNNA has greater application potential in complex earthquake environment.展开更多
基金sponsored by China Geological Survey Project(DD20243193 and DD20230206508).
文摘In Rayleigh wave exploration,the inversion of dispersion curves is a crucial step for obtaining subsurface stratigraphic information,characterized by its multi-parameter and multi-extremum nature.Local optimization algorithms used in dispersion curve inversion are highly dependent on the initial model and are prone to being trapped in local optima,while classical global optimization algorithms often suffer from slow convergence and low solution accuracy.To address these issues,this study introduces the Osprey Optimization Algorithm(OOA),known for its strong global search and local exploitation capabilities,into the inversion of dispersion curves to enhance inversion performance.In noiseless theoretical models,the OOA demonstrates excellent inversion accuracy and stability,accurately recovering model parameters.Even in noisy models,OOA maintains robust performance,achieving high inversion precision under high-noise conditions.In multimode dispersion curve tests,OOA effectively handles higher modes due to its efficient global and local search capabilities,and the inversion results show high consistency with theoretical values.Field data from the Wyoming region in the United States and a landfill site in Italy further verify the practical applicability of the OOA.Comprehensive test results indicate that the OOA outperforms the Particle Swarm Optimization(PSO)algorithm,providing a highly accurate and reliable inversion strategy for dispersion curve inversion.
文摘Combining the adaptive shrinkage genetic algorithm in the feasible region with the imaging of apparent vertical conductance differential, we have inverted the TEM conductive thin layer. The result of the inversion demonstrates that by adaptive shrinkage in the feasible region, the calculation speed accelerates and the calculation precision improves. To a certain extent, in this method we surmount the transient electromagnetic sounding equivalence and reduced equivalence scope. Comparison of the inverted result with the forward curve clearly shows that we can image the conductive thin layer.
基金supported by Specialized Research Fund for the Doctoral Program of Higher Education of China(20110022120004)the Fundamental Research Funds for the Central Universities
文摘For density inversion of gravity anomaly data, once the inversion method is determined, the main factors affecting the inversion result are the inversion parameters and subdivision scheme. A set of reasonable inversion parameters and subdivision scheme can, not only improve the inversion process efficiency, but also ensure inversion result accuracy. The gravity inversion method based on correlation searching and the golden section algorithm is an effective potential field inversion method. It can be used to invert 2D and 3D physical properties with potential data observed on flat or rough surfaces. In this paper, we introduce in detail the density inversion principles based on correlation searching and the golden section algorithm. Considering that the gold section algorithm is not globally optimized. we present a heuristic method to ensure the inversion result is globally optimized. With a series of model tests, we systematically compare and analyze the inversion result efficiency and accuracy with different parameters. Based on the model test results, we conclude the selection principles for each inversion parameter with which the inversion accuracy can be obviously improved.
文摘This paper presents a new algorithm based on the power inversion （PI） and the linearly constrained minimum variance （LCMV）. This algorithm is capable of adjusting the weights of the antenna array in real time to respond to and improve the global positioning system （GPS） received signals coming from the desired directions and at the same time to highly suppress the jammers coming from the other directions. The simulation is performed for fixed and moving jammers. It indicates that this structure can give deeper nulls, more than 115 dB depths for fixed jammers and more than 94 dB depths for moving jammers.
基金Project(41574123)supported by the National Natural Science Foundation of ChinaProject(2015zzts250)supported by the Fundamental Research Funds for the Central Universities,ChinaProject(2013FY110800)supported by the National Basic Research Scientific Program of China
文摘The self-potential method is widely used in environmental and engineering geophysics. Four intelligent optimization algorithms are adopted to design the inversion to interpret self-potential data more accurately and efficiently: simulated annealing, genetic, particle swarm optimization, and ant colony optimization. Using both noise-free and noise-added synthetic data, it is demonstrated that all four intelligent algorithms can perform self-potential data inversion effectively. During the numerical experiments, the model distribution in search space, the relative errors of model parameters, and the elapsed time are recorded to evaluate the performance of the inversion. The results indicate that all the intelligent algorithms have good precision and tolerance to noise. Particle swarm optimization has the fastest convergence during iteration because of its good balanced searching capability between global and local minimisation.
文摘To reduce the computational complexity of matrix inversion, which is the majority of processing in many practical applications, two numerically efficient recursive algorithms （called algorithms I and II, respectively） are presented. Algorithm I is used to calculate the inverse of such a matrix, whose leading principal minors are all nonzero. Algorithm II, whereby, the inverse of an arbitrary nonsingular matrix can be evaluated is derived via improving the algorithm I. The implementation, for algorithm II or I, involves matrix-vector multiplications and vector outer products. These operations are computationally fast and highly parallelizable. MATLAB simulations show that both recursive algorithms are valid.
基金supported by the National Natural Science Foundation of China(No.41374123)
文摘At present, near-surface shear wave velocities are mainly calculated through Rayleigh wave dispersion-curve inversions in engineering surface investigations, but the required calculations pose a highly nonlinear global optimization problem. In order to alleviate the risk of falling into a local optimal solution, this paper introduces a new global optimization method, the shuffle frog-leaping algorithm （SFLA）, into the Rayleigh wave dispersion-curve inversion process. SFLA is a swarm-intelligence-based algorithm that simulates a group of frogs searching for food. It uses a few parameters, achieves rapid convergence, and is capability of effective global searching. In order to test the reliability and calculation performance of SFLA, noise-free and noisy synthetic datasets were inverted. We conducted a comparative analysis with other established algorithms using the noise-free dataset, and then tested the ability of SFLA to cope with data noise. Finally, we inverted a real-world example to examine the applicability of SFLA. Results from both synthetic and field data demonstrated the effectiveness of SFLA in the interpretation of Rayleigh wave dispersion curves. We found that SFLA is superior to the established methods in terms of both reliability and computational efficiency, so it offers great potential to improve our ability to solve geophysical inversion problems.
文摘This paper discusses the inversion of velocity structure and hypocenters location in the Beijing Tianjin Tangshan Zhangjiakou area by genetic algorithm. The hypocenters location of selected earthquakes and crustal structure of this area are obtained using the travel time data of local earthquakes acquired by the Telemetered Seismic Network of Northern China. The mean and standard residuals of hypocenter location acquired by this method are much less than those provided by the report of respective earthquakes. The crustal structure of the first and the second layers obtained interpret the outline of the plain and mountain area in the region successfully and the crustal structure of the third layer nearly coincides with the Moho discontinuity obtained by artificial seismic sounding. These show the genetic algorithm is effective to the inversion of hypocenter location and three dimensional velocity structure.
基金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.
文摘Viscoelastic parameters are becoming more important and their inversion algorithms are studied by many researchers. Genetic algorithms are random, self-adaptive, robust, and heuristic with global search and convergence abilities. Based on the direct VSP wave equation, a genetic algorithm （GA） is introduced to determine the viscoelastic parameters. First, the direct wave equation in frequency is expressed as a function of complex velocity and then the complex velocities estimated by GA inversion. Since the phase velocity and Q-factor both are functions of complex velocity, their values can be computed easily. However, there are so many complex velocities that it is difficult to invert them directly. They can be rewritten as a function of Co and C∞ to reduce the number of parameters during the inversion process. Finally, a theoretical model experiment proves that our algorithm is exact and effective.
基金the Sub-project of National Science and Technology Major Project of China(No.2016ZX05027-002-003)the National Natural Science Foundation of China(No.41404089)+1 种基金the State Key Program of National Natural Science of China(No.41430322)the National Basic Research Program of China(973 Program)(No.2015CB45300)
文摘With the continuous development of full tensor gradiometer （FTG） measurement techniques, three-dimensional （3D） inversion of FTG data is becoming increasingly used in oil and gas exploration. In the fast processing and interpretation of large-scale high-precision data, the use of the graphics processing unit process unit （GPU） and preconditioning methods are very important in the data inversion. In this paper, an improved preconditioned conjugate gradient algorithm is proposed by combining the symmetric successive over-relaxation （SSOR） technique and the incomplete Choleksy decomposition conjugate gradient algorithm （ICCG）. Since preparing the preconditioner requires extra time, a parallel implement based on GPU is proposed. The improved method is then applied in the inversion of noise- contaminated synthetic data to prove its adaptability in the inversion of 3D FTG data. Results show that the parallel SSOR-ICCG algorithm based on NVIDIA Tesla C2050 GPU achieves a speedup of approximately 25 times that of a serial program using a 2.0 GHz Central Processing Unit （CPU）. Real airbome gravity-gradiometry data from Vinton salt dome （south- west Louisiana, USA） are also considered. Good results are obtained, which verifies the efficiency and feasibility of the proposed parallel method in fast inversion of 3D FTG data.
基金Supported by the Key Projects of the National Natural Science Foundation of China(Nos.41076014,U0933005,41176035,40906022,41206029)
文摘We propose a bio-optical inversion model that retrieves the absorption contributions of phytoplankton and colored detrital matter(CDM),as well as the phytoplankton size classes(PSCs),from total minus water absorption spectra.The model is based on three-component separation of phytoplankton size structure and a genetic algorithm.The model performance was tested on two independent datasets(the NASA bio-Optical Marine Algorithm Dataset(NOMAD) and the northern South China Sea(NSCS) dataset).The relationships between the estimated and measured values were strongly linear,especially for aCDM(412),and the Root Mean Square Error(RMSE) of the CDM exponential slope(SCDM) was relatively low.Next,the inversion model was directly applied to in-situ total minus water absorption spectra determined by an underwater meter during a cruise in September 2008,to retrieve the phytoplankton size structure in the seawater.By comparing the measured and retrieved chlorophyll a concentrations,we demonstrated that total and size-specific chlorophyll a concentrations could be retrieved by the model with relatively high accuracy.Finally,we applied the bio-optical inversion model to investigate changes in phytoplankton size structure induced by an anti-cyclonic eddy in the NSCS.
基金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.
基金Financial support by the China Geological Survey Project(Nos.DD20190030,DD20190032)
文摘The practical application of 3D inversion of gravity data requires a lot of computation time and storage space.To solve this problem,we present an integrated optimization algorithm with the following components:(1)targeting high accuracy in the space domain and fast computation in the wavenumber domain,we design a fast 3D forward algorithm with high precision;and(2)taking advantage of the symmetry of the inversion matrix,the main calculation in gravity conjugate gradient inversion is decomposed into two forward calculations,thus optimizing the computational efficiency of 3D gravity inversion.We verify the calculation accuracy and efficiency of the optimization algorithm by testing various grid-number models through numerical simulation experiments.
基金National Nature Science Foundation of China (40334040) & Joint Seismological foundation of CEA (101026)
文摘A genetic algorithm of body waveform inversion is presented for better understanding of crustal and upper mantle structures with deep seismic sounding （DSS） waveform data. General reflection and transmission synthetic seismogram algorithm, which is capable of calculating the response of thin alternating high and low velocity layers, is applied as a solution for forward modeling, and the genetic algorithm is used to find the optimal solution of the inverse problem. Numerical tests suggest that the method has the capability of resolving low-velocity layers, thin alternating high and low velocity layers, and noise suppression. Waveform inversion using P-wave records from Zeku, Xiahe and Lintao shots in the seismic wide-angle reflection/refraction survey along northeastern Qinghai-Xizang （Tibeteau） Plateau has revealed fine structures of the bottom of the upper crust and alternating layers in the middle/lower crust and topmost upper mantle.
文摘The genetic algorithm is useful for solving an inversion of complex nonlinear geophysical equations. The multi-point search of the genetic algorithm makes it easier to find a globally optimal solution and avoid falling into a local extremum. The search efficiency of the genetic algorithm is a key to producing successful solutions in a huge multi-parameter model space. The encoding mechanism of the genetic algorithm affects the searching processes in the evolution. Not all genetic operations perform perfectly in a search under either a binary or decimal encoding system. As such, a standard genetic algorithm (SGA) is sometimes unable to resolve an optimization problem such as a simple geophysical inversion. With the binary encoding system the operation of the crossover may produce more new individuals. The decimal encoding system, on the other hand, makes the mutation generate more new genes. This paper discusses approaches of exploiting the search potentials of genetic operations with different encoding systems and presents a hybrid-encoding mechanism for the genetic algorithm. This is referred to as the hybrid-encoding genetic algorithm (HEGA). The method is based on the routine in which the mutation operation is executed in decimal code and other operations in binary code. HEGA guarantees the birth of better genes by mutation processing with a high probability, so that it is beneficial for resolving the inversions of complicated problems. Synthetic and real-world examples demonstrate the advantages of using HEGA in the inversion of potential-field data.
基金supported by the National Natural Science Foundation of China (No. 41174034)
文摘The firework algorithm（FWA） is a novel swarm intelligence-based method recently proposed for the optimization of multi-parameter, nonlinear functions. Numerical waveform inversion experiments using a synthetic model show that the FWA performs well in both solution quality and efficiency. We apply the FWA in this study to crustal velocity structure inversion using regional seismic waveform data of central Gansu on the northeastern margin of the Qinghai-Tibet plateau. Seismograms recorded from the moment magnitude（MW） 5.4 Minxian earthquake enable obtaining an average crustal velocity model for this region. We initially carried out a series of FWA robustness tests in regional waveform inversion at the same earthquake and station positions across the study region,inverting two velocity structure models, with and without a low-velocity crustal layer; the accuracy of our average inversion results and their standard deviations reveal the advantages of the FWA for the inversion of regional seismic waveforms. We applied the FWA across our study area using three component waveform data recorded by nine broadband permanent seismic stations with epicentral distances ranging between 146 and 437 km. These inversion results show that the average thickness of the crust in this region is 46.75 km, while thicknesses of the sedimentary layer, and the upper, middle, and lower crust are 3.15,15.69, 13.08, and 14.83 km, respectively. Results also show that the P-wave velocities of these layers and the upper mantle are 4.47, 6.07, 6.12, 6.87, and 8.18 km/s,respectively.
文摘The Film Copy Deliverer Problem (FDP),much more difficult than TSP,is a new problem in the combination optimization. In this paper,a new algorithm is introduced. First,the FDP is converted to TSP. Then an evolutionary algorithm based on inversion operator is adopted. Compared with other Genetic Algorithms,it is not only more simple and easier to realize but also faster and more accurate.
基金sponsored by the National Natural Science Foundation of China(41774069,41504047 and 41604054)
文摘The velocity distribution of layers from surface wave dispersion curve is a severely nonlinear program. Base on the Metropolis rule,we improved the simulated annealing algorithm to simultaneously inverse the velocities and thicknesses using the dispersion data and identified the Moho and the bottom of lithosphere. The application to the numerical examples with 5% noise shows the velocity RMS is 1. 56% between the non-linear results and the original models when the condition of selecting method for temperature parameters and initial temperature are satisfied. Using the pure dispersions of Rayleigh wave,the nonlinear inversion has been carried out for S-wave velocities and thicknesses of the vertical profile crossing the Indian Plate,the Qinghai-Tibetan Plateau,and the Tarim Basin. It indicated that the crustal thickness is about 70 km in the Qiangtang block,while in the hinterland of the Qinghai-Tibetan Plateau the lithosphere is relatively thin(~ 130 km)from the velocity values and their offsets.
基金This manuscript is supported by the National Natural Science Foundation of China(No.42174011,41874001 and 42174011).
文摘The traditional genetic algorithm(GA)has unstable inversion results and is easy to fall into the local optimum when inverting fault parameters.Therefore,this article considers the combination of GA with other non-linear algorithms in order to improve the inversion precision of GA.This paper proposes a genetic Nelder-Mead neural network algorithm(GNMNNA).This algorithm uses a neural network algorithm(NNA)to optimize the global search ability of GA.At the same time,the simplex algorithm is used to optimize the local search capability of the GA.Through numerical examples,the stability of the inversion algorithm under different strategies is explored.The experimental results show that the proposed GNMNNA has stronger inversion stability and higher precision compared with the existing algorithms.The effectiveness of GNMNNA is verified by the BodrumeKos earthquake and Monte Cristo Range earthquake.The experimental results show that GNMNNA is superior to GA and NNA in both inversion precision and computational stability.Therefore,GNMNNA has greater application potential in complex earthquake environment.
