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.展开更多
With the development of parallel computing technology,non-linear inversion calculation efficiency has been improving.However,for single-point search-based non-linear inversion methods,the implementation of parallel al...With the development of parallel computing technology,non-linear inversion calculation efficiency has been improving.However,for single-point search-based non-linear inversion methods,the implementation of parallel algorithms is a difficult issue.We introduce the idea of group search to the single-point search-based non-linear inversion algorithm, taking the quantum Monte Carlo method as an example for two-dimensional seismic wave velocity inversion and practical impedance inversion and test the calculation efficiency of using different node numbers.The results show the parallel algorithm in theoretical and practical data inversion is feasible and effective.The parallel algorithm has good versatility. The algorithm efficiency increases with increasing node numbers but the algorithm efficiency rate of increase gradually decreases as the node numbers increase.展开更多
Rayleigh waves have high amplitude, low frequency, and low velocity, which are treated as strong noise to be attenuated in reflected seismic surveys. This study addresses how to identify useful shear wave velocity pro...Rayleigh waves have high amplitude, low frequency, and low velocity, which are treated as strong noise to be attenuated in reflected seismic surveys. This study addresses how to identify useful shear wave velocity profile and stratigraphic information from Rayleigh waves. We choose the Firefly algorithm for inversion of surface waves. The Firefly algorithm, a new type of particle swarm optimization, has the advantages of being robust, highly effective, and allows global searching. This algorithm is feasible and has advantages for use in Rayleigh wave inversion with both synthetic models and field data. The results show that the Firefly algorithm, which is a robust and practical method, can achieve nonlinear inversion of surface waves with high resolution.展开更多
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.展开更多
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.展开更多
In order to obtain accurate probability integration method(PIM) parameters for surface movement of multi-panel mining, a genetic algorithm(GA) was used to optimize the parameters. As the measured surface movement is a...In order to obtain accurate probability integration method(PIM) parameters for surface movement of multi-panel mining, a genetic algorithm(GA) was used to optimize the parameters. As the measured surface movement is affected by more than one mining panel, traditional PIM parameter inversion model is difficult to ensure the reliability of the results due to the complexity of rock movement. With crossover,mutation and selection operators, GA can perform a global optimization search and has high computation efficiency. Compared with the pattern search algorithm, the fitness function can avoid falling into local minima traps. GA reduces the risk of local minima traps which improves the accuracy and reliability with the mutation mechanism. Application at Xuehu colliery shows that GA can be used to inverse the PIM parameters for multi-panel surface movement observation, and reliable results can be obtained. The research provides a new way for back-analysis of PIM parameters for mining subsidence under complex conditions.展开更多
A group of statistical algorithms are proposed for the inversion of the three major components of CaseII waters inthe coastal area of the Huanghai Sea and the East China Sea. The algorithms are based on the in situ da...A group of statistical algorithms are proposed for the inversion of the three major components of CaseII waters inthe coastal area of the Huanghai Sea and the East China Sea. The algorithms are based on the in situ data collected inthe spring of 2003 with strict quality assurance according to NASA ocean bio-optic protocols. These algorithms arethe first ones with quantitative confidence that can be applied for the area. The average relative error of the inversedand in situ measured components' concentrations are: Chl-a about 37%, total suspended matter (TSM) about 25%,respectively. This preliminary result is quite satisfactory for CaseII waters, although some aspects in the modelneed further study. The sensitivity of the input error of 5% to remote sensing reflectance (Rrs) is also analyzed andit shows the algorithms are quite stable. The algorithms show a large difference with Tassans local SeaWiFSalgorithms for different waters, except for the Chl-a algorithm.展开更多
The applications of intelligent techniques have increased exponentially in recent days to study most of the non-linear parameters. In particular, the behavior of earth resembles the non- linearity applications. An eff...The applications of intelligent techniques have increased exponentially in recent days to study most of the non-linear parameters. In particular, the behavior of earth resembles the non- linearity applications. An efficient tool is needed for the interpretation of geophysical parameters to study the subsurface of the earth. Artificial Neural Networks （ANN） perform certain tasks if the structure of the network is modified accordingly for the purpose it has been used. The three most robust networks were taken and comparatively analyzed for their performance to choose the appropriate network. The single- layer feed-forward neural network with the back propagation algorithm is chosen as one of the well- suited networks after comparing the results. Initially, certain synthetic data sets of all three-layer curves have been taken tk^r training the network, and the network is validated by the field datasets collected from Tuticorin Coastal Region （78°7＇30＂E and 8°48＇45＂N）, Tamil Nadu, India. The interpretation has been done successfully using the corresponding learning algorithm in the present study. With proper training of back propagation networks, it tends to give the resistivity and thickness of the subsurface layer model of the field resistivity data concerning the synthetic data trained earlier in the appropriate network. The network is trained with more Vertical Electrical Sounding （VES） data, and this trained network is demon- strated by the field data. Groundwater table depth also has been modeled.展开更多
Precisely estimating the state of health(SOH)of lithium-ion batteries is essential for battery management systems(BMS),as it plays a key role in ensuring the safe and reliable operation of battery systems.However,curr...Precisely estimating the state of health(SOH)of lithium-ion batteries is essential for battery management systems(BMS),as it plays a key role in ensuring the safe and reliable operation of battery systems.However,current SOH estimation methods often overlook the valuable temperature information that can effectively characterize battery aging during capacity degradation.Additionally,the Elman neural network,which is commonly employed for SOH estimation,exhibits several drawbacks,including slow training speed,a tendency to become trapped in local minima,and the initialization of weights and thresholds using pseudo-random numbers,leading to unstable model performance.To address these issues,this study addresses the challenge of precise and effective SOH detection by proposing a method for estimating the SOH of lithium-ion batteries based on differential thermal voltammetry(DTV)and an SSA-Elman neural network.Firstly,two health features(HFs)considering temperature factors and battery voltage are extracted fromthe differential thermal voltammetry curves and incremental capacity curves.Next,the Sparrow Search Algorithm(SSA)is employed to optimize the initial weights and thresholds of the Elman neural network,forming the SSA-Elman neural network model.To validate the performance,various neural networks,including the proposed SSA-Elman network,are tested using the Oxford battery aging dataset.The experimental results demonstrate that the method developed in this study achieves superior accuracy and robustness,with a mean absolute error(MAE)of less than 0.9%and a rootmean square error(RMSE)below 1.4%.展开更多
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.展开更多
Seismic anisotropy provides important constraints on deformation patterns of Earth＇s material. Rayleigh wave dispersion data with azimuthal anisotropy can be used to invert for depth-dependent shear wavespeed azimuth...Seismic anisotropy provides important constraints on deformation patterns of Earth＇s material. Rayleigh wave dispersion data with azimuthal anisotropy can be used to invert for depth-dependent shear wavespeed azimuthal anisotropy, therefore reflecting depth-varying deformation patterns in the crust and upper mantle. In this study, we propose a two-step method that uses the Neighborhood Algorithm （NA） for the point-wise inversion of depth-dependent shear wavespeeds and azimuthal anisotropy from Rayleigh wave azimuthally anisotropic dispersion data. The first step employs the NA to estimate depth- dependent Vsv （or the elastic parameter L） as well as their uncertainties from the isotropic part Rayleigh wave dispersion data. In the second step, we first adopt a difference scheme to compute approximate Rayleigh-wave phase velocity sensitivity kernels to azimuthally anisotropic parameters with respect to the velocity model obtained in the first step. Then we perform the NA to estimate the azimuthally anisotropic parameters Gc/L and Gs/L at depths separately from the corresponding cosine and sine terms of the azimuthally anisotropic dispersion data. Finally, we compute the depth-dependent magnitude and fast polariza- tion azimuth of shear wavespeed azimuthal anisotropy. The use of the global search NA and Bayesian analysis allows for more reliable estimates of depth-dependent shear wavespeeds and azimuthal anisotropy as well as their uncertainties.We illustrate the inversion method using the azimuthally anisotropic dispersion data in SE Tibet, where we find apparent changes of fast axes of shear wavespeed azimuthal anisotropy between the crust and uppermost mantle.展开更多
The fi rst arrival waveform inversion(FAWI)has a strong nonlinearity due to the objective function using L2 parametrization.When the initial velocity is not accurate,the inversion can easily fall into local minima.In ...The fi rst arrival waveform inversion(FAWI)has a strong nonlinearity due to the objective function using L2 parametrization.When the initial velocity is not accurate,the inversion can easily fall into local minima.In the full waveform inversion method,adding a cross-correlation function to the objective function can eff ectively reduce the nonlinearity of the inversion process.In this paper,the nonlinearity of this process is reduced by introducing the correlation objective function into the FAWI and by deriving the corresponding gradient formula.We then combine the first-arrival wave travel-time tomography with the FAWI to form a set of inversion processes.This paper uses the limited memory Broyden-Fletcher-Goldfarb-Shanno(L-BFGS)algorithm to improve the computational effi ciency of inversion and solve the problem of the low effi ciency of the FAWI method.The overthrust model and fi eld data test show that the method used in this paper can eff ectively reduce the nonlinearity of inversion and improve the inversion calculation effi ciency at the same time.展开更多
Smooth constraint is important in linear inversion, but it is difficult to apply directly to model parameters in genetic algorithms. If the model parameters are smoothed in iteration, the diversity of models will be g...Smooth constraint is important in linear inversion, but it is difficult to apply directly to model parameters in genetic algorithms. If the model parameters are smoothed in iteration, the diversity of models will be greatly suppressed and all the models in population will tend to equal in a few iterations, so the optimal solution meeting requirement can not be obtained. In this paper, an indirect smooth constraint technique is introduced to genetic inversion. In this method, the new models produced in iteration are smoothed, then used as theoretical models in calculation of misfit function, but in process of iteration only the original models are used in order to keep the diversity of models. The technique is effective in inversion of surface wave and receiver function. Using this technique, we invert the phase velocity of Raleigh wave in the Tibetan Plateau, revealing the horizontal variation of S wave velocity structure near the center of the Tibetan Plateau. The results show that the S wave velocity in the north is relatively lower than that in the south. For most paths there is a lower velocity zone with 12-25 km thick at the depth of 15-40 km. The lower velocity zone in upper mantle is located below the depth of 100 km, and the thickness is usually 40-80 km, but for a few paths reach to 100 km thick. Among the area of Ando, Maqi and Ushu stations, there is an obvious lower velocity zone with the lowest velocity of 4.2-4.3 km/s at the depth of 90-230 km. Based on the S wave velocity structures of different paths and former data, we infer that the subduction of the Indian Plate is delimited nearby the Yarlung Zangbo suture zone.展开更多
The Pamir Plateau,at the northwestern margin of the Tibetan Plateau,is a key region for investigating continental collision and plateau uplifting.To probe its deep structure,we collected seismic data from 263 stations...The Pamir Plateau,at the northwestern margin of the Tibetan Plateau,is a key region for investigating continental collision and plateau uplifting.To probe its deep structure,we collected seismic data from 263 stations across 11 research projects.We applied cross-correlation to noise data and extracted surface wave dispersion data from cross-correlation functions.The extracted dispersion data were subsequently inverted using a 3-D transdimensional Bayesian inversion method(rj-3 DMcMC).The inversion result reveals several crustal low-velocity zones(LVZs).Their formation is likely related to crustal thickening,the exposure of gneiss domes,and thicker sedimentary sequences compared to surrounding areas.In the lower crust and upper mantle,the LVZs in southern Pamir and southeastern Karakoram evolve into high-velocity zones,which expand northeastward with increasing depth.This suggests northward underthrusting of the Indian Plate.We also analyzed the Moho using both the standard deviation of S-wave velocity and the S-wave velocity structure.Results show that significant variations in velocity standard deviation reliably delineate the Moho interface.展开更多
Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this cha...Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.展开更多
To improve the efficiency and accuracy of path planning for fan inspection tasks in thermal power plants,this paper proposes an intelligent inspection robot path planning scheme based on an improved A^(*)algorithm.The...To improve the efficiency and accuracy of path planning for fan inspection tasks in thermal power plants,this paper proposes an intelligent inspection robot path planning scheme based on an improved A^(*)algorithm.The inspection robot utilizes multiple sensors to monitor key parameters of the fans,such as vibration,noise,and bearing temperature,and upload the data to the monitoring center.The robot’s inspection path employs the improved A^(*)algorithm,incorporating obstacle penalty terms,path reconstruction,and smoothing optimization techniques,thereby achieving optimal path planning for the inspection robot in complex environments.Simulation results demonstrate that the improved A^(*)algorithm significantly outperforms the traditional A^(*)algorithm in terms of total path distance,smoothness,and detour rate,effectively improving the execution efficiency of inspection tasks.展开更多
基金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.
基金supported by National Key S&T Special Projects of Marine Carbonate(No.2008ZX05000-004)CNPC Projects(No.2008E-0610-10)
文摘With the development of parallel computing technology,non-linear inversion calculation efficiency has been improving.However,for single-point search-based non-linear inversion methods,the implementation of parallel algorithms is a difficult issue.We introduce the idea of group search to the single-point search-based non-linear inversion algorithm, taking the quantum Monte Carlo method as an example for two-dimensional seismic wave velocity inversion and practical impedance inversion and test the calculation efficiency of using different node numbers.The results show the parallel algorithm in theoretical and practical data inversion is feasible and effective.The parallel algorithm has good versatility. The algorithm efficiency increases with increasing node numbers but the algorithm efficiency rate of increase gradually decreases as the node numbers increase.
基金supported by the National Basic Research Program of China(No.2013CB228602)the National Science and Technology Major Project of China(No.2011ZX05004-003)the National High Technology Research Program of China(No.2013AA064202)
文摘Rayleigh waves have high amplitude, low frequency, and low velocity, which are treated as strong noise to be attenuated in reflected seismic surveys. This study addresses how to identify useful shear wave velocity profile and stratigraphic information from Rayleigh waves. We choose the Firefly algorithm for inversion of surface waves. The Firefly algorithm, a new type of particle swarm optimization, has the advantages of being robust, highly effective, and allows global searching. This algorithm is feasible and has advantages for use in Rayleigh wave inversion with both synthetic models and field data. The results show that the Firefly algorithm, which is a robust and practical method, can achieve nonlinear inversion of surface waves with high resolution.
文摘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.
基金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.
基金provided by the National Natural Science Foundation of China(No.51404272)the Hunan Province Key Laboratory of Coal Resources Clean-Utilization and Mine Environment Protection(No.E21224)
文摘In order to obtain accurate probability integration method(PIM) parameters for surface movement of multi-panel mining, a genetic algorithm(GA) was used to optimize the parameters. As the measured surface movement is affected by more than one mining panel, traditional PIM parameter inversion model is difficult to ensure the reliability of the results due to the complexity of rock movement. With crossover,mutation and selection operators, GA can perform a global optimization search and has high computation efficiency. Compared with the pattern search algorithm, the fitness function can avoid falling into local minima traps. GA reduces the risk of local minima traps which improves the accuracy and reliability with the mutation mechanism. Application at Xuehu colliery shows that GA can be used to inverse the PIM parameters for multi-panel surface movement observation, and reliable results can be obtained. The research provides a new way for back-analysis of PIM parameters for mining subsidence under complex conditions.
基金The work was supported by the Subsystem of Calibration and ValidationHY-I Ground Application System+1 种基金National Satellite Ocean Application Service(NSOAS)China High-Tech“863"Project under contract Nos 2001AA636010 and 2001AA637010/7030.
文摘A group of statistical algorithms are proposed for the inversion of the three major components of CaseII waters inthe coastal area of the Huanghai Sea and the East China Sea. The algorithms are based on the in situ data collected inthe spring of 2003 with strict quality assurance according to NASA ocean bio-optic protocols. These algorithms arethe first ones with quantitative confidence that can be applied for the area. The average relative error of the inversedand in situ measured components' concentrations are: Chl-a about 37%, total suspended matter (TSM) about 25%,respectively. This preliminary result is quite satisfactory for CaseII waters, although some aspects in the modelneed further study. The sensitivity of the input error of 5% to remote sensing reflectance (Rrs) is also analyzed andit shows the algorithms are quite stable. The algorithms show a large difference with Tassans local SeaWiFSalgorithms for different waters, except for the Chl-a algorithm.
文摘The applications of intelligent techniques have increased exponentially in recent days to study most of the non-linear parameters. In particular, the behavior of earth resembles the non- linearity applications. An efficient tool is needed for the interpretation of geophysical parameters to study the subsurface of the earth. Artificial Neural Networks （ANN） perform certain tasks if the structure of the network is modified accordingly for the purpose it has been used. The three most robust networks were taken and comparatively analyzed for their performance to choose the appropriate network. The single- layer feed-forward neural network with the back propagation algorithm is chosen as one of the well- suited networks after comparing the results. Initially, certain synthetic data sets of all three-layer curves have been taken tk^r training the network, and the network is validated by the field datasets collected from Tuticorin Coastal Region （78°7＇30＂E and 8°48＇45＂N）, Tamil Nadu, India. The interpretation has been done successfully using the corresponding learning algorithm in the present study. With proper training of back propagation networks, it tends to give the resistivity and thickness of the subsurface layer model of the field resistivity data concerning the synthetic data trained earlier in the appropriate network. The network is trained with more Vertical Electrical Sounding （VES） data, and this trained network is demon- strated by the field data. Groundwater table depth also has been modeled.
基金supported by the National Natural Science Foundation of China(NSFC)under Grant(No.51677058).
文摘Precisely estimating the state of health(SOH)of lithium-ion batteries is essential for battery management systems(BMS),as it plays a key role in ensuring the safe and reliable operation of battery systems.However,current SOH estimation methods often overlook the valuable temperature information that can effectively characterize battery aging during capacity degradation.Additionally,the Elman neural network,which is commonly employed for SOH estimation,exhibits several drawbacks,including slow training speed,a tendency to become trapped in local minima,and the initialization of weights and thresholds using pseudo-random numbers,leading to unstable model performance.To address these issues,this study addresses the challenge of precise and effective SOH detection by proposing a method for estimating the SOH of lithium-ion batteries based on differential thermal voltammetry(DTV)and an SSA-Elman neural network.Firstly,two health features(HFs)considering temperature factors and battery voltage are extracted fromthe differential thermal voltammetry curves and incremental capacity curves.Next,the Sparrow Search Algorithm(SSA)is employed to optimize the initial weights and thresholds of the Elman neural network,forming the SSA-Elman neural network model.To validate the performance,various neural networks,including the proposed SSA-Elman network,are tested using the Oxford battery aging dataset.The experimental results demonstrate that the method developed in this study achieves superior accuracy and robustness,with a mean absolute error(MAE)of less than 0.9%and a rootmean square error(RMSE)below 1.4%.
基金The Project Advanced Buildings technology in a Dense Urban Environment of Hong Kong Polytechnic University the Application Lab of Digital Seismic Wave Data in Center for Analysis and Prediction+1 种基金 China Seismological Bureau and the Earthquake Prediction
基金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.
基金supported by National Natural Science Foundation of China (Grant number 41222028)China National Special Fund for Earthquake Scientific Research in Public Interest (Grant number 201508008)the Fundamental Research Funds for the Central Universities (WK2080000053)
文摘Seismic anisotropy provides important constraints on deformation patterns of Earth＇s material. Rayleigh wave dispersion data with azimuthal anisotropy can be used to invert for depth-dependent shear wavespeed azimuthal anisotropy, therefore reflecting depth-varying deformation patterns in the crust and upper mantle. In this study, we propose a two-step method that uses the Neighborhood Algorithm （NA） for the point-wise inversion of depth-dependent shear wavespeeds and azimuthal anisotropy from Rayleigh wave azimuthally anisotropic dispersion data. The first step employs the NA to estimate depth- dependent Vsv （or the elastic parameter L） as well as their uncertainties from the isotropic part Rayleigh wave dispersion data. In the second step, we first adopt a difference scheme to compute approximate Rayleigh-wave phase velocity sensitivity kernels to azimuthally anisotropic parameters with respect to the velocity model obtained in the first step. Then we perform the NA to estimate the azimuthally anisotropic parameters Gc/L and Gs/L at depths separately from the corresponding cosine and sine terms of the azimuthally anisotropic dispersion data. Finally, we compute the depth-dependent magnitude and fast polariza- tion azimuth of shear wavespeed azimuthal anisotropy. The use of the global search NA and Bayesian analysis allows for more reliable estimates of depth-dependent shear wavespeeds and azimuthal anisotropy as well as their uncertainties.We illustrate the inversion method using the azimuthally anisotropic dispersion data in SE Tibet, where we find apparent changes of fast axes of shear wavespeed azimuthal anisotropy between the crust and uppermost mantle.
基金supported by the Major Scientific and Technological Project of PetroChina (ZD2019-183-003)Project of National Natural Science Foundation of China (42074133)+1 种基金the Fundamental Research Funds for the Central Universities (19CX02056A)Project of State Key Laboratory of Shale Oil and Gas Enrichment Mechanisms and Effective Development (33550000-21-FW0399-0009)
文摘The fi rst arrival waveform inversion(FAWI)has a strong nonlinearity due to the objective function using L2 parametrization.When the initial velocity is not accurate,the inversion can easily fall into local minima.In the full waveform inversion method,adding a cross-correlation function to the objective function can eff ectively reduce the nonlinearity of the inversion process.In this paper,the nonlinearity of this process is reduced by introducing the correlation objective function into the FAWI and by deriving the corresponding gradient formula.We then combine the first-arrival wave travel-time tomography with the FAWI to form a set of inversion processes.This paper uses the limited memory Broyden-Fletcher-Goldfarb-Shanno(L-BFGS)algorithm to improve the computational effi ciency of inversion and solve the problem of the low effi ciency of the FAWI method.The overthrust model and fi eld data test show that the method used in this paper can eff ectively reduce the nonlinearity of inversion and improve the inversion calculation effi ciency at the same time.
基金State Natural Science Foundation (49874021).Contribution No. 01FE2002, Institute of Geophysics, China Seismological Bureau.
文摘Smooth constraint is important in linear inversion, but it is difficult to apply directly to model parameters in genetic algorithms. If the model parameters are smoothed in iteration, the diversity of models will be greatly suppressed and all the models in population will tend to equal in a few iterations, so the optimal solution meeting requirement can not be obtained. In this paper, an indirect smooth constraint technique is introduced to genetic inversion. In this method, the new models produced in iteration are smoothed, then used as theoretical models in calculation of misfit function, but in process of iteration only the original models are used in order to keep the diversity of models. The technique is effective in inversion of surface wave and receiver function. Using this technique, we invert the phase velocity of Raleigh wave in the Tibetan Plateau, revealing the horizontal variation of S wave velocity structure near the center of the Tibetan Plateau. The results show that the S wave velocity in the north is relatively lower than that in the south. For most paths there is a lower velocity zone with 12-25 km thick at the depth of 15-40 km. The lower velocity zone in upper mantle is located below the depth of 100 km, and the thickness is usually 40-80 km, but for a few paths reach to 100 km thick. Among the area of Ando, Maqi and Ushu stations, there is an obvious lower velocity zone with the lowest velocity of 4.2-4.3 km/s at the depth of 90-230 km. Based on the S wave velocity structures of different paths and former data, we infer that the subduction of the Indian Plate is delimited nearby the Yarlung Zangbo suture zone.
基金supported by the National Natural Science Foundation of China(Grant No.42174126)the Alliance of International Science Organizations(ANSO)Project(Grant No.ANSO-CR-PP2022-04)+1 种基金the Deep Earth Probe and Mineral Resources Exploration National Science and Technology Major Project(Grant Nos.2024ZD1002206,2024ZD1002201)Key R&D Program of Xinjiang Uyghur Autonomous Region(Grant No.2024B03013-2)。
文摘The Pamir Plateau,at the northwestern margin of the Tibetan Plateau,is a key region for investigating continental collision and plateau uplifting.To probe its deep structure,we collected seismic data from 263 stations across 11 research projects.We applied cross-correlation to noise data and extracted surface wave dispersion data from cross-correlation functions.The extracted dispersion data were subsequently inverted using a 3-D transdimensional Bayesian inversion method(rj-3 DMcMC).The inversion result reveals several crustal low-velocity zones(LVZs).Their formation is likely related to crustal thickening,the exposure of gneiss domes,and thicker sedimentary sequences compared to surrounding areas.In the lower crust and upper mantle,the LVZs in southern Pamir and southeastern Karakoram evolve into high-velocity zones,which expand northeastward with increasing depth.This suggests northward underthrusting of the Indian Plate.We also analyzed the Moho using both the standard deviation of S-wave velocity and the S-wave velocity structure.Results show that significant variations in velocity standard deviation reliably delineate the Moho interface.
基金funded by the National Key R&D Program of China(Grant No.2022YFC2903904)the National Natural Science Foundation of China(Grant Nos.51904057 and U1906208).
文摘Due to the heterogeneity of rock masses and the variability of in situ stress,the traditional linear inversion method is insufficiently accurate to achieve high accuracy of the in situ stress field.To address this challenge,nonlinear stress boundaries for a numerical model are determined through regression analysis of a series of nonlinear coefficient matrices,which are derived from the bubbling method.Considering the randomness and flexibility of the bubbling method,a parametric study is conducted to determine recommended ranges for these parameters,including the standard deviation(σb)of bubble radii,the non-uniform coefficient matrix number(λ)for nonlinear stress boundaries,and the number(m)and positions of in situ stress measurement points.A model case study provides a reference for the selection of these parameters.Additionally,when the nonlinear in situ stress inversion method is employed,stress distortion inevitably occurs near model boundaries,aligning with the Saint Venant's principle.Two strategies are proposed accordingly:employing a systematic reduction of nonlinear coefficients to achieve high inversion accuracy while minimizing significant stress distortion,and excluding regions with severe stress distortion near the model edges while utilizing the central part of the model for subsequent simulations.These two strategies have been successfully implemented in the nonlinear in situ stress inversion of the Xincheng Gold Mine and have achieved higher inversion accuracy than the linear method.Specifically,the linear and nonlinear inversion methods yield root mean square errors(RMSE)of 4.15 and 3.2,and inversion relative errors(δAve)of 22.08%and 17.55%,respectively.Therefore,the nonlinear inversion method outperforms the traditional multiple linear regression method,even in the presence of a systematic reduction in the nonlinear stress boundaries.
文摘To improve the efficiency and accuracy of path planning for fan inspection tasks in thermal power plants,this paper proposes an intelligent inspection robot path planning scheme based on an improved A^(*)algorithm.The inspection robot utilizes multiple sensors to monitor key parameters of the fans,such as vibration,noise,and bearing temperature,and upload the data to the monitoring center.The robot’s inspection path employs the improved A^(*)algorithm,incorporating obstacle penalty terms,path reconstruction,and smoothing optimization techniques,thereby achieving optimal path planning for the inspection robot in complex environments.Simulation results demonstrate that the improved A^(*)algorithm significantly outperforms the traditional A^(*)algorithm in terms of total path distance,smoothness,and detour rate,effectively improving the execution efficiency of inspection tasks.
