The traditional genetic algorithm(GA)has unstable inversion results and is easy to fall into the local optimum when inverting fault parameters.Therefore,this article considers the combination of GA with other non-line...The traditional genetic algorithm(GA)has unstable inversion results and is easy to fall into the local optimum when inverting fault parameters.Therefore,this article considers the combination of GA with other non-linear algorithms in order to improve the inversion precision of GA.This paper proposes a genetic Nelder-Mead neural network algorithm(GNMNNA).This algorithm uses a neural network algorithm(NNA)to optimize the global search ability of GA.At the same time,the simplex algorithm is used to optimize the local search capability of the GA.Through numerical examples,the stability of the inversion algorithm under different strategies is explored.The experimental results show that the proposed GNMNNA has stronger inversion stability and higher precision compared with the existing algorithms.The effectiveness of GNMNNA is verified by the BodrumeKos earthquake and Monte Cristo Range earthquake.The experimental results show that GNMNNA is superior to GA and NNA in both inversion precision and computational stability.Therefore,GNMNNA has greater application potential in complex earthquake environment.展开更多
The great Tancheng earthquake of M81/2 occurred in 1668 was the largest seismic event ever recorded in history in eastern China. This study determines the fault geometry of this earthquake by inverting seismological ...The great Tancheng earthquake of M81/2 occurred in 1668 was the largest seismic event ever recorded in history in eastern China. This study determines the fault geometry of this earthquake by inverting seismological data of present-day moderate-small earthquakes in the focal area. We relocated those earthquakes with the double-difference method and found focal mechanism solutions using gird test method. The inversion results are as follows： the strike is 21.6°, the dip angle is 89.5°, the slip angle is 170°, the fault length is about 160 km, the lower-boundary depth is about 32 km and the buried depth of upper boundary is about 4 km. This shows that the seismic fault is a NNE-trending upright right-lateral strike-slip fault and has cut through the crust. Moreover, the surface seismic fault, intensity distribution of the earthquake, earthquake-depth distribution and seismic-wave velocity profile in the focal area all verified our study result.展开更多
The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault par...The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault parameters inversion has nonlinear characteristics, and the gradient-based optimizer(GBO) has the characteristics of fast convergence speed and falling into local optimum hardly. This paper applies GBO algorithm to simulated earthquakes and real LuShan earthquakes in the nonlinear inversion of the Okada model to obtain the source parameters. The simulated earthquake experiment results show that the algorithm is stable, and the seismic source parameters obtained by GBO are slightly closer to the true value than the multi peak particle swarm optimization(MPSO). In the 2013 LuShan earthquake experiment, the root mean square error between the deformation after forwarding of fault parameters obtained by the introduced GBO algorithm and the surface observation deformation was 3.703 mm, slightly better than 3.708 mm calculated by the MPSO. Moreover, the inversion result of GBO algorithm is better than MPSO algorithm in stability. The above results show that the introduced GBO algorithm has a certain practical application value in seismic fault source parameters inversion.展开更多
A Mw6.4 earthquake occurred in L＇Aquila, central Italy at 1：32：42 （UTC）, April 6, 2009. We quickly obtained the moment tensor solution of the earthquake by inverting the P waveforms of broadband recordings from t...A Mw6.4 earthquake occurred in L＇Aquila, central Italy at 1：32：42 （UTC）, April 6, 2009. We quickly obtained the moment tensor solution of the earthquake by inverting the P waveforms of broadband recordings from the global seismographic network （GSN） stations using the quick technique of moment tensor inversion, and further inferred that the nodal plane of strike 132°, dip 53° and rake -103° is the seismogenic fault.展开更多
We review three derivative-free methods developed for uncertainty estimation of non-linear error propagation, namely, MC(Monte Carlo), SUT(scaled unscented transformation), and SI(sterling interpolation). In order to ...We review three derivative-free methods developed for uncertainty estimation of non-linear error propagation, namely, MC(Monte Carlo), SUT(scaled unscented transformation), and SI(sterling interpolation). In order to avoid preset parameters like as these three methods need, we introduce a new method to uncertainty estimation for the first time, namely, SCR(spherical cubature rule), which is no need for setting parameters. By theoretical derivation, we prove that the precision of uncertainty obtained by SCR can reach second-order. We conduct four synthetic experiments, for the first two experiments, the results obtained by SCR are consistent with the other three methods with optimal setting parameters, but SCR is easier to operate than other three methods, which verifies the superiority of SCR in calculating the uncertainty. For the third experiment, real-time calculation is required, so the MC is hardly feasible. For the forth experiment, the SCR is applied to the inversion of seismic fault parameter which is a common problem in geophysics, and we study the sensitivity of surface displacements to fault parameters with errors. Our results show that the uncertainty of the surface displacements is the magnitude of ±10 mm when the fault length contains a variance of 0.01 km^(2).展开更多
基金This manuscript is supported by the National Natural Science Foundation of China(No.42174011,41874001 and 42174011).
文摘The traditional genetic algorithm(GA)has unstable inversion results and is easy to fall into the local optimum when inverting fault parameters.Therefore,this article considers the combination of GA with other non-linear algorithms in order to improve the inversion precision of GA.This paper proposes a genetic Nelder-Mead neural network algorithm(GNMNNA).This algorithm uses a neural network algorithm(NNA)to optimize the global search ability of GA.At the same time,the simplex algorithm is used to optimize the local search capability of the GA.Through numerical examples,the stability of the inversion algorithm under different strategies is explored.The experimental results show that the proposed GNMNNA has stronger inversion stability and higher precision compared with the existing algorithms.The effectiveness of GNMNNA is verified by the BodrumeKos earthquake and Monte Cristo Range earthquake.The experimental results show that GNMNNA is superior to GA and NNA in both inversion precision and computational stability.Therefore,GNMNNA has greater application potential in complex earthquake environment.
基金supported by the National Natural Science Foundation of China(No.90814002)the Natural Science Foundation of Shandong Province(No.Y2005E02)
文摘The great Tancheng earthquake of M81/2 occurred in 1668 was the largest seismic event ever recorded in history in eastern China. This study determines the fault geometry of this earthquake by inverting seismological data of present-day moderate-small earthquakes in the focal area. We relocated those earthquakes with the double-difference method and found focal mechanism solutions using gird test method. The inversion results are as follows： the strike is 21.6°, the dip angle is 89.5°, the slip angle is 170°, the fault length is about 160 km, the lower-boundary depth is about 32 km and the buried depth of upper boundary is about 4 km. This shows that the seismic fault is a NNE-trending upright right-lateral strike-slip fault and has cut through the crust. Moreover, the surface seismic fault, intensity distribution of the earthquake, earthquake-depth distribution and seismic-wave velocity profile in the focal area all verified our study result.
基金the National Natural Science Foundation of China(Nos.42174011and 41874001).
文摘The use of geodetic observation data for seismic fault parameters inversion is the research hotspot of geodetic inversion, and it is also the focus of studying the mechanism of earthquake occurrence. Seismic fault parameters inversion has nonlinear characteristics, and the gradient-based optimizer(GBO) has the characteristics of fast convergence speed and falling into local optimum hardly. This paper applies GBO algorithm to simulated earthquakes and real LuShan earthquakes in the nonlinear inversion of the Okada model to obtain the source parameters. The simulated earthquake experiment results show that the algorithm is stable, and the seismic source parameters obtained by GBO are slightly closer to the true value than the multi peak particle swarm optimization(MPSO). In the 2013 LuShan earthquake experiment, the root mean square error between the deformation after forwarding of fault parameters obtained by the introduced GBO algorithm and the surface observation deformation was 3.703 mm, slightly better than 3.708 mm calculated by the MPSO. Moreover, the inversion result of GBO algorithm is better than MPSO algorithm in stability. The above results show that the introduced GBO algorithm has a certain practical application value in seismic fault source parameters inversion.
基金No.09FE3007 of Institute of Geophysics,China Earthquake Administration
文摘A Mw6.4 earthquake occurred in L＇Aquila, central Italy at 1：32：42 （UTC）, April 6, 2009. We quickly obtained the moment tensor solution of the earthquake by inverting the P waveforms of broadband recordings from the global seismographic network （GSN） stations using the quick technique of moment tensor inversion, and further inferred that the nodal plane of strike 132°, dip 53° and rake -103° is the seismogenic fault.
基金supported by the National Natural Science Foundation of China (41721003, 41974022, 41774024, 41874001)Open Research Fund Program of the Key Laboratory of Geospace Environment and Geodesy, Ministry of Education, China(20-02-05)
文摘We review three derivative-free methods developed for uncertainty estimation of non-linear error propagation, namely, MC(Monte Carlo), SUT(scaled unscented transformation), and SI(sterling interpolation). In order to avoid preset parameters like as these three methods need, we introduce a new method to uncertainty estimation for the first time, namely, SCR(spherical cubature rule), which is no need for setting parameters. By theoretical derivation, we prove that the precision of uncertainty obtained by SCR can reach second-order. We conduct four synthetic experiments, for the first two experiments, the results obtained by SCR are consistent with the other three methods with optimal setting parameters, but SCR is easier to operate than other three methods, which verifies the superiority of SCR in calculating the uncertainty. For the third experiment, real-time calculation is required, so the MC is hardly feasible. For the forth experiment, the SCR is applied to the inversion of seismic fault parameter which is a common problem in geophysics, and we study the sensitivity of surface displacements to fault parameters with errors. Our results show that the uncertainty of the surface displacements is the magnitude of ±10 mm when the fault length contains a variance of 0.01 km^(2).
