Landfill leaks pose a serious threat to environmental health,risking the contamination of both groundwater and soil resources.Accurate investigation of these sites is essential for implementing effective prevention an...Landfill leaks pose a serious threat to environmental health,risking the contamination of both groundwater and soil resources.Accurate investigation of these sites is essential for implementing effective prevention and control measures.The self-potential(SP)stands out for its sensitivity to contamination plumes,offering a solution for monitoring and detecting the movement and seepage of subsurface pollutants.However,traditional SP inversion techniques heavily rely on precise subsurface resistivity information.In this study,we propose the Attention U-Net deep learning network for rapid SP inversion.By incorporating an attention mechanism,this algorithm effectively learns the relationship between array-style SP data and the location and extent of subsurface contaminated sources.We designed a synthetic landfill model with a heterogeneous resistivity structure to assess the performance of Attention U-Net deep learning network.Additionally,we conducted further validation using a laboratory model to assess its practical applicability.The results demonstrate that the algorithm is not solely dependent on resistivity information,enabling effective locating of the source distribution,even in models with intricate subsurface structures.Our work provides a promising tool for SP data processing,enhancing the applicability of this method in the field of near-subsurface environmental monitoring.展开更多
Most earth-dam failures are mainly due to seepage,and an accurate assessment of the permeability coefficient provides an indication to avoid a disaster.Parametric uncertainties are encountered in the seepage analysis,...Most earth-dam failures are mainly due to seepage,and an accurate assessment of the permeability coefficient provides an indication to avoid a disaster.Parametric uncertainties are encountered in the seepage analysis,and may be reduced by an inverse procedure that calibrates the simulation results to observations on the real system being simulated.This work proposes an adaptive Bayesian inversion method solved using artificial neural network(ANN)based Markov Chain Monte Carlo simulation.The optimized surrogate model achieves a coefficient of determination at 0.98 by ANN with 247 samples,whereby the computational workload can be greatly reduced.It is also significant to balance the accuracy and efficiency of the ANN model by adaptively updating the sample database.The enrichment samples are obtained from the posterior distribution after iteration,which allows a more accurate and rapid manner to the target posterior.The method was then applied to the hydraulic analysis of an earth dam.After calibrating the global permeability coefficient of the earth dam with the pore water pressure at the downstream unsaturated location,it was validated by the pore water pressure monitoring values at the upstream saturated location.In addition,the uncertainty in the permeability coefficient was reduced,from 0.5 to 0.05.It is shown that the provision of adequate prior information is valuable for improving the efficiency of the Bayesian inversion.展开更多
We assembled approximately 328 seismic records. The data set was from 4 digitally recording long-period and broadband stations of CDSN. We carried out the inversion based on the partitioned waveform inversion (PWI). I...We assembled approximately 328 seismic records. The data set was from 4 digitally recording long-period and broadband stations of CDSN. We carried out the inversion based on the partitioned waveform inversion (PWI). It partitions the large-scale optimization problem into a number of independent small-scale problems. We adopted surface waveform inversion with an equal block (2((2() discretization in order to acquire the images of shear velocity structure at different depths (from surface to 430 km) in the crust and upper-mantle. The resolution of all these anomalies has been established with (check-board( resolution tests. These results show significant difference in velocity, lithosphere and asthenosphere structure between South China Sea and its adjacent regions.展开更多
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.展开更多
Prticle swarm optimization(PSO)is adopted to invert the self-potential anomalies of simple geometry.Taking the vertical semi-infinite cylinder model as an example,the model parameters are first inverted using standard...Prticle swarm optimization(PSO)is adopted to invert the self-potential anomalies of simple geometry.Taking the vertical semi-infinite cylinder model as an example,the model parameters are first inverted using standard particle swarm optimization(SPSO),and then the searching behavior of the particle swarm is discussed and the change of the particles’distribution during the iteration process is studied.The existence of different particle behaviors enables the particle swarm to explore the searching space more comprehensively,thus PSO achieves remarkable results in the inversion of SP anomalies.Finally,six improved PSOs aiming at improving the inversion accuracy and the convergence speed by changing the update of particle positions,inertia weights and learning factors are introduced for the inversion of the cylinder model,and the effectiveness of these algorithms is verified by numerical experiments.The inversion results show that these improved PSOs successfully give the model parameters which are very close to the theoretical value,and simultaneously provide guidance when determining which strategy is suitable for the inversion of the regular polarized bodies and similar geophysical problems.展开更多
Estimating the spatial distribution of coseismic slip is an ill-posed inverse problem, and solutions may be extremely oscillatory due to measurement errors without any constraints on the coseismic slip distribution. I...Estimating the spatial distribution of coseismic slip is an ill-posed inverse problem, and solutions may be extremely oscillatory due to measurement errors without any constraints on the coseismic slip distribution. In order to obtain stable solution for coseismic slip inversion, regularization method with smoothness-constrained was imposed. Trade-off parameter in regularized inversion, which balances the minimization of the data misfit and model roughness, should be a critical procedure to achieve both resolution and stability. Then, the active constraint balancing approach is adopted, in which the trade-off parameter is regarded as a spatial variable at each model parameter and automatically determined via the model resolution matrix and the spread function. Numerical experiments for a synthetical model indicate that regularized inversion using active constraint balancing approach can provides stable inversion results and have low sensitivity to the knowledge of the exact character of the Gaussian noise. Regularized inversion combined with active constraint balancing approach is conducted on the 2005 Nias earthquake. The released moment based on the estimated coseismic slip distribution is 9.91×1021 N·m, which is equivalent to a moment magnitude of 8.6 and almost identical to the value determined by USGS. The inversion results for synthetic coseismic uniform-slip model and the 2005 earthquake show that smoothness-constrained regularized inversion method combined with active constraint balancing approach is effective, and can be reasonable to reconstruct coseismic slip distribution on fault.展开更多
To monitor growth and predict the yield of rice over a large area, the chlorophyll contents in the rice canopy were estimated using the unmanned aerial vehicle(UAV) remote sensing technology. In this work, multi-spect...To monitor growth and predict the yield of rice over a large area, the chlorophyll contents in the rice canopy were estimated using the unmanned aerial vehicle(UAV) remote sensing technology. In this work, multi-spectral image information of the rice crop was obtained using a 6-channel multi-spectral camera mounted on a fixed wing UAV, which was flown 600 m above the ground, between 11: 00-14: 00 on a sunny day in summer. The measured chlorophyll values were collected as sample sets. The s-REP index was screened out to estimate chlorophyll contents through the analysis of six kinds of spectral indexes of chlorophyll estimated capacity. An inversion model of the chlorophyll contents was then built using the least square support vector regression(LS-SVR)algorithm, with calibration and prediction R-square values of 0.89 and 0.83, respectively. Finally, remote sensing mapping for a UAV image of the Fangzheng County Dexter Rice Planting Park was accomplished using the inversion model. The inversion and measured values were then compared using regression fitting. R-square and root-mean-square error of the fitting model were 0.79 and 2.39,respectively. The results demonstrated that accurate estimation of rice-canopy chlorophyll contents was feasible using the LS-SVR inversion model developed using the s-REP vegetation index.展开更多
To improve magnetotelluric(MT)nonlinear inversion accuracy and stability,this work introduces the deep belief network(DBN)algorithm.Firstly,a network frame is set up for training in different 2D MT models.The network ...To improve magnetotelluric(MT)nonlinear inversion accuracy and stability,this work introduces the deep belief network(DBN)algorithm.Firstly,a network frame is set up for training in different 2D MT models.The network inputs are the apparent resistivities of known models,and the outputs are the model parameters.The optimal network structure is achieved by determining the numbers of hidden layers and network nodes.Secondly,the learning process of the DBN is implemented to obtain the optimal solution of network connection weights for known geoelectric models.Finally,the trained DBN is verified through inversion tests,in which the network inputs are the apparent resistivities of unknown models,and the outputs are the corresponding model parameters.The experiment results show that the DBN can make full use of the global searching capability of the restricted Boltzmann machine(RBM)unsupervised learning and the local optimization of the back propagation(BP)neural network supervised learning.Comparing to the traditional neural network inversion,the calculation accuracy and stability of the DBN for MT data inversion are improved significantly.And the tests on synthetic data reveal that this method can be applied to MT data inversion and achieve good results compared with the least-square regularization inversion.展开更多
According to the effective sampling of particles and the particles impoverishment caused by re-sampling in particle filter,an improved particle filtering algorithm based on observation inversion optimal sampling was p...According to the effective sampling of particles and the particles impoverishment caused by re-sampling in particle filter,an improved particle filtering algorithm based on observation inversion optimal sampling was proposed. Firstly,virtual observations were generated from the latest observation,and two sampling strategies were presented. Then,the previous time particles were sampled by utilizing the function inversion relationship between observation and system state. Finally,the current time particles were generated on the basis of the previous time particles and the system one-step state transition model. By the above method,sampling particles can make full use of the latest observation information and the priori modeling information,so that they further approximate the true state. The theoretical analysis and experimental results show that the new algorithm filtering accuracy and real-time outperform obviously the standard particle filter,the extended Kalman particle filter and the unscented particle filter.展开更多
Nonlinear resistivity inversion requires efficient artificial neural network(ANN)model for better inversion results.An evolutionary BP neural network(BPNN)approach based on differential evolution(DE)algorithm was pres...Nonlinear resistivity inversion requires efficient artificial neural network(ANN)model for better inversion results.An evolutionary BP neural network(BPNN)approach based on differential evolution(DE)algorithm was presented,which was able to improve global search ability for resistivity tomography 2-D nonlinear inversion.In the proposed method,Tent equation was applied to obtain automatic parameter settings in DE and the restricted parameter Fcrit was used to enhance the ability of converging to global optimum.An implementation of proposed DE-BPNN was given,the network had one hidden layer with 52 nodes and it was trained on 36 datasets and tested on another 4 synthetic datasets.Two abnormity models were used to verify the feasibility and effectiveness of the proposed method,the results show that the proposed DE-BP algorithm has better performance than BP,conventional DE-BP and other chaotic DE-BP methods in stability and accuracy,and higher imaging quality than least square inversion.展开更多
Geological structures often exhibit smooth characteristics away from sharp discontinuities. One aim of geophysical inversion is to recover information about the smooth structures as well as about the sharp discontinui...Geological structures often exhibit smooth characteristics away from sharp discontinuities. One aim of geophysical inversion is to recover information about the smooth structures as well as about the sharp discontinuities. Because no specific operator can provide a perfect sparse representation of complicated geological models, hyper-parameter regularization inversion based on the iterative split Bregman method was used to recover the features of both smooth and sharp geological structures. A novel preconditioned matrix was proposed, which counteracted the natural decay of the sensitivity matrix and its inverse matrix was calculated easily. Application of the algorithm to synthetic data produces density models that are good representations of the designed models. The results show that the algorithm proposed is feasible and effective.展开更多
A new ultrasound contrast imaging technique was proposed for eliminating the harmonic components from the emission signal transmitted by the broadband ultrasonic system.Reversal phase-inversion pulse was used for the ...A new ultrasound contrast imaging technique was proposed for eliminating the harmonic components from the emission signal transmitted by the broadband ultrasonic system.Reversal phase-inversion pulse was used for the first time to separate the contrast harmonics from the harmonics in the emission signal to improve the detection of contrast micro-bubbles.Based on the nonlinear acoustic theory of finite-amplitude effects and the associated distortion of the propagating wave,the Bessel-Fubini series model was applied to describe the nonlinear propagation effects of the reversal phase-inversion pulse,and the Church's equation for zero-thickness encapsulation model was used to produce the scattering-pulse of the bubble.For harmonic imaging,the experiment was performed using a 64-element linear array,which was simulated by Field II.The results show that the harmonic components from the emission signal can be completely cancelled,and the harmonics generated by the nonlinear propagation of the wave through the tissue,can be reduced by 15-30 dB.Compared with the short pulse,the reversal phase-inversion pulse can improve the contrast and definition of the harmonic image significantly.展开更多
To improve the global search ability and imaging quality of electrical resistivity imaging(ERI) inversion, a two-stage learning ICPSO algorithm of radial basis function neural network(RBFNN) based on information crite...To improve the global search ability and imaging quality of electrical resistivity imaging(ERI) inversion, a two-stage learning ICPSO algorithm of radial basis function neural network(RBFNN) based on information criterion(IC) and particle swarm optimization(PSO) is presented. In the proposed method, IC is applied to obtain the hidden layer structure by calculating the optimal IC value automatically and PSO algorithm is used to optimize the centers and widths of the radial basis functions in the hidden layer. Meanwhile, impacts of different information criteria to the inversion results are compared, and an implementation of the proposed ICPSO algorithm is given. The optimized neural network has one hidden layer with 261 nodes selected by AKAIKE's information criterion(AIC) and it is trained on 32 data sets and tested on another 8 synthetic data sets. Two complex synthetic examples are used to verify the feasibility and effectiveness of the proposed method with two learning stages. The results show that the proposed method has better performance and higher imaging quality than three-layer and four-layer back propagation neural networks(BPNNs) and traditional least square(LS) inversion.展开更多
For the problem of sensor faults and actuator faults in aircraft attitude control,this paper proposes a fault tolerant control(FTC)scheme based on extended state observer(ESO)and nonlinear dynamic inversion(NDI).First...For the problem of sensor faults and actuator faults in aircraft attitude control,this paper proposes a fault tolerant control(FTC)scheme based on extended state observer(ESO)and nonlinear dynamic inversion(NDI).First,two ESOs are designed to estimate sensor faults and actuator faults respectively.Second,the angular rate signal is reconstructed according to the estimation of sensor faults.Third,in angular rate loop,NDI is designed based on reconstruction of angular rate signals and estimation of actuator faults.The FTC scheme proposed in this paper is testified through numerical simulations.The results show that it is feasible and has good fault tolerant ability.展开更多
The classical elastic impedance （EI） inversion method, however, is based on the L2-norm misfit function and considerably sensitive to outliers, assuming the noise of the seismic data to be the Guassian-distribution....The classical elastic impedance （EI） inversion method, however, is based on the L2-norm misfit function and considerably sensitive to outliers, assuming the noise of the seismic data to be the Guassian-distribution. So we have developed a more robust elastic impedance inversion based on the Ll-norm misfit function, and the noise is assumed to be non-Gaussian. Meanwhile, some regularization methods including the sparse constraint regularization and elastic impedance point constraint regularization are incorporated to improve the ill-posed characteristics of the seismic inversion problem. Firstly, we create the Ll-norm misfit objective function of pre-stack inversion problem based on the Bayesian scheme within the sparse constraint regularization and elastic impedance point constraint regularization. And then, we obtain more robust elastic impedances of different angles which are less sensitive to outliers in seismic data by using the IRLS strategy. Finally, we extract the P-wave and S-wave velocity and density by using the more stable parameter extraction method. Tests on synthetic data show that the P-wave and S-wave velocity and density parameters are still estimated reasonable with moderate noise. A test on the real data set shows that compared to the results of the classical elastic impedance inversion method, the estimated results using the proposed method can get better lateral continuity and more distinct show of the gas, verifying the feasibility and stability of the method.展开更多
In this paper, broadband multi-frequencies matched-field inversion method is used to determine the environmental parameters in shallow water. According to different conditions, several broadband objective functions ar...In this paper, broadband multi-frequencies matched-field inversion method is used to determine the environmental parameters in shallow water. According to different conditions, several broadband objective functions are presented. Using ASIAEX2001 experiment data and genetic algorithms, environmental parameters are obtained, especially in sediment.展开更多
Optimization efficiencies and mechanisms of simulated annealing, genetic algorithm, differential evolution and downhill simplex differential evolution are compared and analyzed. Simulated annealing and genetic algorit...Optimization efficiencies and mechanisms of simulated annealing, genetic algorithm, differential evolution and downhill simplex differential evolution are compared and analyzed. Simulated annealing and genetic algorithm use a directed random process to search the parameter space for an optimal solution. They include the ability to avoid local minima, but as no gradient information is used, searches may be relatively inefficient. Differential evolution uses information from a distance and azimuth between individuals of a population to search the parameter space, the initial search is effective, but the search speed decreases quickly because differential information between the individuals of population vanishes. Local downhill simplex and global differential evolution methods are developed separately, and combined to produce a hybrid downhill simplex differential evolution algorithm. The hybrid algorithm is sensitive to gradients of the object function and search of the parameter space is effective. These algorithms are applied to the matched field inversion with synthetic data. Optimal values of the parameters, the final values of object function and inversion time is presented and compared.展开更多
The study of induced polarization （IP） information extraction from magnetotelluric （MT） sounding data is of great and practical significance to the exploitation of deep mineral, oil and gas resources. The linear i...The study of induced polarization （IP） information extraction from magnetotelluric （MT） sounding data is of great and practical significance to the exploitation of deep mineral, oil and gas resources. The linear inversion method, which has been given priority in previous research on the IP information extraction method, has three main problems as follows： 1） dependency on the initial model, 2） easily falling into the local minimum, and 3） serious non-uniqueness of solutions. Taking the nonlinearity and nonconvexity of IP information extraction into consideration, a two-stage CO-PSO minimum structure inversion method using compute unified distributed architecture （CUDA） is proposed. On one hand, a novel Cauchy oscillation particle swarm optimization （CO-PSO） algorithm is applied to extract nonlinear IP information from MT sounding data, which is implemented as a parallel algorithm within CUDA computing architecture; on the other hand, the impact of the polarizability on the observation data is strengthened by introducing a second stage inversion process, and the regularization parameter is applied in the fitness function of PSO algorithm to solve the problem of multi-solution in inversion. The inversion simulation results of polarization layers in different strata of various geoelectric models show that the smooth models of resistivity and IP parameters can be obtained by the proposed algorithm, the results of which are relatively stable and accurate. The experiment results added with noise indicate that this method is robust to Gaussian white noise. Compared with the traditional PSO and GA algorithm, the proposed algorithm has more efficiency and better inversion results.展开更多
基金Projects(42174170,41874145,72088101)supported by the National Natural Science Foundation of ChinaProject(CX20200228)supported by the Hunan Provincial Innovation Foundation for Postgraduate,China。
文摘Landfill leaks pose a serious threat to environmental health,risking the contamination of both groundwater and soil resources.Accurate investigation of these sites is essential for implementing effective prevention and control measures.The self-potential(SP)stands out for its sensitivity to contamination plumes,offering a solution for monitoring and detecting the movement and seepage of subsurface pollutants.However,traditional SP inversion techniques heavily rely on precise subsurface resistivity information.In this study,we propose the Attention U-Net deep learning network for rapid SP inversion.By incorporating an attention mechanism,this algorithm effectively learns the relationship between array-style SP data and the location and extent of subsurface contaminated sources.We designed a synthetic landfill model with a heterogeneous resistivity structure to assess the performance of Attention U-Net deep learning network.Additionally,we conducted further validation using a laboratory model to assess its practical applicability.The results demonstrate that the algorithm is not solely dependent on resistivity information,enabling effective locating of the source distribution,even in models with intricate subsurface structures.Our work provides a promising tool for SP data processing,enhancing the applicability of this method in the field of near-subsurface environmental monitoring.
基金Project(202006430012)supported by the China Scholarship Council。
文摘Most earth-dam failures are mainly due to seepage,and an accurate assessment of the permeability coefficient provides an indication to avoid a disaster.Parametric uncertainties are encountered in the seepage analysis,and may be reduced by an inverse procedure that calibrates the simulation results to observations on the real system being simulated.This work proposes an adaptive Bayesian inversion method solved using artificial neural network(ANN)based Markov Chain Monte Carlo simulation.The optimized surrogate model achieves a coefficient of determination at 0.98 by ANN with 247 samples,whereby the computational workload can be greatly reduced.It is also significant to balance the accuracy and efficiency of the ANN model by adaptively updating the sample database.The enrichment samples are obtained from the posterior distribution after iteration,which allows a more accurate and rapid manner to the target posterior.The method was then applied to the hydraulic analysis of an earth dam.After calibrating the global permeability coefficient of the earth dam with the pore water pressure at the downstream unsaturated location,it was validated by the pore water pressure monitoring values at the upstream saturated location.In addition,the uncertainty in the permeability coefficient was reduced,from 0.5 to 0.05.It is shown that the provision of adequate prior information is valuable for improving the efficiency of the Bayesian inversion.
基金State Natural Scientific Foundation (49734150) and National High Performance Computation Foundation.
文摘We assembled approximately 328 seismic records. The data set was from 4 digitally recording long-period and broadband stations of CDSN. We carried out the inversion based on the partitioned waveform inversion (PWI). It partitions the large-scale optimization problem into a number of independent small-scale problems. We adopted surface waveform inversion with an equal block (2((2() discretization in order to acquire the images of shear velocity structure at different depths (from surface to 430 km) in the crust and upper-mantle. The resolution of all these anomalies has been established with (check-board( resolution tests. These results show significant difference in velocity, lithosphere and asthenosphere structure between South China Sea and its adjacent regions.
基金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.
基金Projects(41874145,72088101)supported by the National Natural Science Foundation of China。
文摘Prticle swarm optimization(PSO)is adopted to invert the self-potential anomalies of simple geometry.Taking the vertical semi-infinite cylinder model as an example,the model parameters are first inverted using standard particle swarm optimization(SPSO),and then the searching behavior of the particle swarm is discussed and the change of the particles’distribution during the iteration process is studied.The existence of different particle behaviors enables the particle swarm to explore the searching space more comprehensively,thus PSO achieves remarkable results in the inversion of SP anomalies.Finally,six improved PSOs aiming at improving the inversion accuracy and the convergence speed by changing the update of particle positions,inertia weights and learning factors are introduced for the inversion of the cylinder model,and the effectiveness of these algorithms is verified by numerical experiments.The inversion results show that these improved PSOs successfully give the model parameters which are very close to the theoretical value,and simultaneously provide guidance when determining which strategy is suitable for the inversion of the regular polarized bodies and similar geophysical problems.
基金Projects(41604111,41541036) supported by the National Natural Science Foundation of China
文摘Estimating the spatial distribution of coseismic slip is an ill-posed inverse problem, and solutions may be extremely oscillatory due to measurement errors without any constraints on the coseismic slip distribution. In order to obtain stable solution for coseismic slip inversion, regularization method with smoothness-constrained was imposed. Trade-off parameter in regularized inversion, which balances the minimization of the data misfit and model roughness, should be a critical procedure to achieve both resolution and stability. Then, the active constraint balancing approach is adopted, in which the trade-off parameter is regarded as a spatial variable at each model parameter and automatically determined via the model resolution matrix and the spread function. Numerical experiments for a synthetical model indicate that regularized inversion using active constraint balancing approach can provides stable inversion results and have low sensitivity to the knowledge of the exact character of the Gaussian noise. Regularized inversion combined with active constraint balancing approach is conducted on the 2005 Nias earthquake. The released moment based on the estimated coseismic slip distribution is 9.91×1021 N·m, which is equivalent to a moment magnitude of 8.6 and almost identical to the value determined by USGS. The inversion results for synthetic coseismic uniform-slip model and the 2005 earthquake show that smoothness-constrained regularized inversion method combined with active constraint balancing approach is effective, and can be reasonable to reconstruct coseismic slip distribution on fault.
基金Supported by the National Key R&D Program of China(2016YFD0300610)
文摘To monitor growth and predict the yield of rice over a large area, the chlorophyll contents in the rice canopy were estimated using the unmanned aerial vehicle(UAV) remote sensing technology. In this work, multi-spectral image information of the rice crop was obtained using a 6-channel multi-spectral camera mounted on a fixed wing UAV, which was flown 600 m above the ground, between 11: 00-14: 00 on a sunny day in summer. The measured chlorophyll values were collected as sample sets. The s-REP index was screened out to estimate chlorophyll contents through the analysis of six kinds of spectral indexes of chlorophyll estimated capacity. An inversion model of the chlorophyll contents was then built using the least square support vector regression(LS-SVR)algorithm, with calibration and prediction R-square values of 0.89 and 0.83, respectively. Finally, remote sensing mapping for a UAV image of the Fangzheng County Dexter Rice Planting Park was accomplished using the inversion model. The inversion and measured values were then compared using regression fitting. R-square and root-mean-square error of the fitting model were 0.79 and 2.39,respectively. The results demonstrated that accurate estimation of rice-canopy chlorophyll contents was feasible using the LS-SVR inversion model developed using the s-REP vegetation index.
基金Project(41304090)supported by the National Natural Science Foundation of ChinaProject(2016YFC0303104)supported by the National Key Research and Development Project of ChinaProject(DY135-S1-1-07)supported by Ocean 13th Five-Year International Marine Resources Survey and Development of China
文摘To improve magnetotelluric(MT)nonlinear inversion accuracy and stability,this work introduces the deep belief network(DBN)algorithm.Firstly,a network frame is set up for training in different 2D MT models.The network inputs are the apparent resistivities of known models,and the outputs are the model parameters.The optimal network structure is achieved by determining the numbers of hidden layers and network nodes.Secondly,the learning process of the DBN is implemented to obtain the optimal solution of network connection weights for known geoelectric models.Finally,the trained DBN is verified through inversion tests,in which the network inputs are the apparent resistivities of unknown models,and the outputs are the corresponding model parameters.The experiment results show that the DBN can make full use of the global searching capability of the restricted Boltzmann machine(RBM)unsupervised learning and the local optimization of the back propagation(BP)neural network supervised learning.Comparing to the traditional neural network inversion,the calculation accuracy and stability of the DBN for MT data inversion are improved significantly.And the tests on synthetic data reveal that this method can be applied to MT data inversion and achieve good results compared with the least-square regularization inversion.
基金Project(60634030) supported by the Key Project of the National Natural Science Foundation of ChinaProject(60702066) supported by the National Natural Science Foundation of China+1 种基金Project (2007ZC53037) supported by Aviation Science Foundation of ChinaProject(CASC0214) supported by the Space-Flight Innovation Foundation of China
文摘According to the effective sampling of particles and the particles impoverishment caused by re-sampling in particle filter,an improved particle filtering algorithm based on observation inversion optimal sampling was proposed. Firstly,virtual observations were generated from the latest observation,and two sampling strategies were presented. Then,the previous time particles were sampled by utilizing the function inversion relationship between observation and system state. Finally,the current time particles were generated on the basis of the previous time particles and the system one-step state transition model. By the above method,sampling particles can make full use of the latest observation information and the priori modeling information,so that they further approximate the true state. The theoretical analysis and experimental results show that the new algorithm filtering accuracy and real-time outperform obviously the standard particle filter,the extended Kalman particle filter and the unscented particle filter.
基金Project(20120162110015)supported by the Research Fund for the Doctoral Program of Higher Education,ChinaProject(41004053)supported by the National Natural Science Foundation of ChinaProject(12c0241)supported by Scientific Research Fund of Hunan Provincial Education Department,China
文摘Nonlinear resistivity inversion requires efficient artificial neural network(ANN)model for better inversion results.An evolutionary BP neural network(BPNN)approach based on differential evolution(DE)algorithm was presented,which was able to improve global search ability for resistivity tomography 2-D nonlinear inversion.In the proposed method,Tent equation was applied to obtain automatic parameter settings in DE and the restricted parameter Fcrit was used to enhance the ability of converging to global optimum.An implementation of proposed DE-BPNN was given,the network had one hidden layer with 52 nodes and it was trained on 36 datasets and tested on another 4 synthetic datasets.Two abnormity models were used to verify the feasibility and effectiveness of the proposed method,the results show that the proposed DE-BP algorithm has better performance than BP,conventional DE-BP and other chaotic DE-BP methods in stability and accuracy,and higher imaging quality than least square inversion.
基金Projects(41174061,41374120)supported by the National Natural Science Foundation of China
文摘Geological structures often exhibit smooth characteristics away from sharp discontinuities. One aim of geophysical inversion is to recover information about the smooth structures as well as about the sharp discontinuities. Because no specific operator can provide a perfect sparse representation of complicated geological models, hyper-parameter regularization inversion based on the iterative split Bregman method was used to recover the features of both smooth and sharp geological structures. A novel preconditioned matrix was proposed, which counteracted the natural decay of the sensitivity matrix and its inverse matrix was calculated easily. Application of the algorithm to synthetic data produces density models that are good representations of the designed models. The results show that the algorithm proposed is feasible and effective.
基金Project(20070533131) supported by the National Research Foundation for the Doctoral Program of Higher Education of ChinaProject(50275150) supported by the National Natural Science Foundation of China
文摘A new ultrasound contrast imaging technique was proposed for eliminating the harmonic components from the emission signal transmitted by the broadband ultrasonic system.Reversal phase-inversion pulse was used for the first time to separate the contrast harmonics from the harmonics in the emission signal to improve the detection of contrast micro-bubbles.Based on the nonlinear acoustic theory of finite-amplitude effects and the associated distortion of the propagating wave,the Bessel-Fubini series model was applied to describe the nonlinear propagation effects of the reversal phase-inversion pulse,and the Church's equation for zero-thickness encapsulation model was used to produce the scattering-pulse of the bubble.For harmonic imaging,the experiment was performed using a 64-element linear array,which was simulated by Field II.The results show that the harmonic components from the emission signal can be completely cancelled,and the harmonics generated by the nonlinear propagation of the wave through the tissue,can be reduced by 15-30 dB.Compared with the short pulse,the reversal phase-inversion pulse can improve the contrast and definition of the harmonic image significantly.
基金Project(41374118)supported by the National Natural Science Foundation,ChinaProject(20120162110015)supported by Research Fund for the Doctoral Program of Higher Education,China+3 种基金Project(2015M580700)supported by the China Postdoctoral Science Foundation,ChinaProject(2016JJ3086)supported by the Hunan Provincial Natural Science Foundation,ChinaProject(2015JC3067)supported by the Hunan Provincial Science and Technology Program,ChinaProject(15B138)supported by the Scientific Research Fund of Hunan Provincial Education Department,China
文摘To improve the global search ability and imaging quality of electrical resistivity imaging(ERI) inversion, a two-stage learning ICPSO algorithm of radial basis function neural network(RBFNN) based on information criterion(IC) and particle swarm optimization(PSO) is presented. In the proposed method, IC is applied to obtain the hidden layer structure by calculating the optimal IC value automatically and PSO algorithm is used to optimize the centers and widths of the radial basis functions in the hidden layer. Meanwhile, impacts of different information criteria to the inversion results are compared, and an implementation of the proposed ICPSO algorithm is given. The optimized neural network has one hidden layer with 261 nodes selected by AKAIKE's information criterion(AIC) and it is trained on 32 data sets and tested on another 8 synthetic data sets. Two complex synthetic examples are used to verify the feasibility and effectiveness of the proposed method with two learning stages. The results show that the proposed method has better performance and higher imaging quality than three-layer and four-layer back propagation neural networks(BPNNs) and traditional least square(LS) inversion.
基金supported by the Chinese Aviation Science Fund(20160757001)the National Natural Science Foundation of China(10577012)。
文摘For the problem of sensor faults and actuator faults in aircraft attitude control,this paper proposes a fault tolerant control(FTC)scheme based on extended state observer(ESO)and nonlinear dynamic inversion(NDI).First,two ESOs are designed to estimate sensor faults and actuator faults respectively.Second,the angular rate signal is reconstructed according to the estimation of sensor faults.Third,in angular rate loop,NDI is designed based on reconstruction of angular rate signals and estimation of actuator faults.The FTC scheme proposed in this paper is testified through numerical simulations.The results show that it is feasible and has good fault tolerant ability.
基金Projects(U1562215,41674130,41404088)supported by the National Natural Science Foundation of ChinaProjects(2013CB228604,2014CB239201)supported by the National Basic Research Program of China+1 种基金Projects(2016ZX05027004-001,2016ZX05002006-009)supported by the National Oil and Gas Major Projects of ChinaProject(15CX08002A)supported by the Fundamental Research Funds for the Central Universities,China
文摘The classical elastic impedance （EI） inversion method, however, is based on the L2-norm misfit function and considerably sensitive to outliers, assuming the noise of the seismic data to be the Guassian-distribution. So we have developed a more robust elastic impedance inversion based on the Ll-norm misfit function, and the noise is assumed to be non-Gaussian. Meanwhile, some regularization methods including the sparse constraint regularization and elastic impedance point constraint regularization are incorporated to improve the ill-posed characteristics of the seismic inversion problem. Firstly, we create the Ll-norm misfit objective function of pre-stack inversion problem based on the Bayesian scheme within the sparse constraint regularization and elastic impedance point constraint regularization. And then, we obtain more robust elastic impedances of different angles which are less sensitive to outliers in seismic data by using the IRLS strategy. Finally, we extract the P-wave and S-wave velocity and density by using the more stable parameter extraction method. Tests on synthetic data show that the P-wave and S-wave velocity and density parameters are still estimated reasonable with moderate noise. A test on the real data set shows that compared to the results of the classical elastic impedance inversion method, the estimated results using the proposed method can get better lateral continuity and more distinct show of the gas, verifying the feasibility and stability of the method.
文摘In this paper, broadband multi-frequencies matched-field inversion method is used to determine the environmental parameters in shallow water. According to different conditions, several broadband objective functions are presented. Using ASIAEX2001 experiment data and genetic algorithms, environmental parameters are obtained, especially in sediment.
文摘Optimization efficiencies and mechanisms of simulated annealing, genetic algorithm, differential evolution and downhill simplex differential evolution are compared and analyzed. Simulated annealing and genetic algorithm use a directed random process to search the parameter space for an optimal solution. They include the ability to avoid local minima, but as no gradient information is used, searches may be relatively inefficient. Differential evolution uses information from a distance and azimuth between individuals of a population to search the parameter space, the initial search is effective, but the search speed decreases quickly because differential information between the individuals of population vanishes. Local downhill simplex and global differential evolution methods are developed separately, and combined to produce a hybrid downhill simplex differential evolution algorithm. The hybrid algorithm is sensitive to gradients of the object function and search of the parameter space is effective. These algorithms are applied to the matched field inversion with synthetic data. Optimal values of the parameters, the final values of object function and inversion time is presented and compared.
基金Projects(41604117,41204054)supported by the National Natural Science Foundation of ChinaProjects(20110490149,2015M580700)supported by the Research Fund for the Doctoral Program of Higher Education,China+1 种基金Project(2015zzts064)supported by the Fundamental Research Funds for the Central Universities,ChinaProject(16B147)supported by the Scientific Research Fund of Hunan Provincial Education Department,China
文摘The study of induced polarization （IP） information extraction from magnetotelluric （MT） sounding data is of great and practical significance to the exploitation of deep mineral, oil and gas resources. The linear inversion method, which has been given priority in previous research on the IP information extraction method, has three main problems as follows： 1） dependency on the initial model, 2） easily falling into the local minimum, and 3） serious non-uniqueness of solutions. Taking the nonlinearity and nonconvexity of IP information extraction into consideration, a two-stage CO-PSO minimum structure inversion method using compute unified distributed architecture （CUDA） is proposed. On one hand, a novel Cauchy oscillation particle swarm optimization （CO-PSO） algorithm is applied to extract nonlinear IP information from MT sounding data, which is implemented as a parallel algorithm within CUDA computing architecture; on the other hand, the impact of the polarizability on the observation data is strengthened by introducing a second stage inversion process, and the regularization parameter is applied in the fitness function of PSO algorithm to solve the problem of multi-solution in inversion. The inversion simulation results of polarization layers in different strata of various geoelectric models show that the smooth models of resistivity and IP parameters can be obtained by the proposed algorithm, the results of which are relatively stable and accurate. The experiment results added with noise indicate that this method is robust to Gaussian white noise. Compared with the traditional PSO and GA algorithm, the proposed algorithm has more efficiency and better inversion results.
