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 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.展开更多
On December 18, 2023, the M_(S)6.2 Jishishan earthquake occurred in the northeastern region of the QinghaiXizang Plateau, causing heavy casualties and property damage in Gansu and Qinghai Provinces. In this study,we i...On December 18, 2023, the M_(S)6.2 Jishishan earthquake occurred in the northeastern region of the QinghaiXizang Plateau, causing heavy casualties and property damage in Gansu and Qinghai Provinces. In this study,we integrate space imaging geodesy, finite fault inversion, and back-projection methods to decipher its rupture property, including fault geometry, coseismic slip distribution, rupture direction, and propagation speed. The results reveal that the seismogenic fault dips to the southwest at an angle of 29°. The major slip asperity is dominated by reverse slip and is concentrated within a depth range of 7–16 km, which explains the significant uplift near the epicenter observed by both the Sentinel-1 ascending and descending In SAR data. Moreover, the teleseismic array waveforms indicate a northwest propagating rupture with an overall slow rupture velocity of~1.91 km/s(AK array) or 1.01 km/s(AU array).展开更多
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).展开更多
Based on the repeated gravity observation data from 1996 to 2007 from the Longmenshan gravity network, which has been dealt with by adjustment processing, the benchmark interference removal and impact of elevation cha...Based on the repeated gravity observation data from 1996 to 2007 from the Longmenshan gravity network, which has been dealt with by adjustment processing, the benchmark interference removal and impact of elevation changes removal, and by using the 3-D inversion method to reflect underground density, we analyze the characteristics of Longmenshan regional dynamic crustal density at depths of 25km, 20km and 15kin. The results show that in the Wenchuan earthquake preparation process, the regional density field showed marked characteristics both in time and space distribution. From the point of time process, the density change trend in the ten years before the earthquake presents a periodic change pattern： steady phase, dramatic stage, slow reducing phase and slow increase phase. The degree of density changes is from large to small, which means that earthquake gestation has reached the final stage. From the point of space distribution, density change distribution has a tendency of ＂dispersion--relative concentration＂, this shows that before the earthquake, the entropy of the underground density field was decreased. In addition, dramatic density changes often occur in the Longmenshan fault zone and western Sichuan plateau. Also, with the increase of depth, the trend of density change is more and more obvious. Through comparative analysis, the influence of density change on gravity is much bigger than that from height change.展开更多
基金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.
基金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.
基金supported by the Open Fund of Hubei Luojia Laboratory(230100015)the Strategic Priority Research Program of the Chinese Academy of Sciences(XDB41000000)the Knowledge Innovation Program of Wuhan-Shuguang Project(2023010201020281).
文摘On December 18, 2023, the M_(S)6.2 Jishishan earthquake occurred in the northeastern region of the QinghaiXizang Plateau, causing heavy casualties and property damage in Gansu and Qinghai Provinces. In this study,we integrate space imaging geodesy, finite fault inversion, and back-projection methods to decipher its rupture property, including fault geometry, coseismic slip distribution, rupture direction, and propagation speed. The results reveal that the seismogenic fault dips to the southwest at an angle of 29°. The major slip asperity is dominated by reverse slip and is concentrated within a depth range of 7–16 km, which explains the significant uplift near the epicenter observed by both the Sentinel-1 ascending and descending In SAR data. Moreover, the teleseismic array waveforms indicate a northwest propagating rupture with an overall slow rupture velocity of~1.91 km/s(AK array) or 1.01 km/s(AU array).
基金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).
基金funded by the National Natural Science Foundation of China(41330314)Projects of Science for Earthquake Resilience(XH15049Y)+1 种基金National Science and Technology Support Program of China(2012BAK19B02,2012BAK19B03)Special Research Foundation for Seismology(201108009)
文摘Based on the repeated gravity observation data from 1996 to 2007 from the Longmenshan gravity network, which has been dealt with by adjustment processing, the benchmark interference removal and impact of elevation changes removal, and by using the 3-D inversion method to reflect underground density, we analyze the characteristics of Longmenshan regional dynamic crustal density at depths of 25km, 20km and 15kin. The results show that in the Wenchuan earthquake preparation process, the regional density field showed marked characteristics both in time and space distribution. From the point of time process, the density change trend in the ten years before the earthquake presents a periodic change pattern： steady phase, dramatic stage, slow reducing phase and slow increase phase. The degree of density changes is from large to small, which means that earthquake gestation has reached the final stage. From the point of space distribution, density change distribution has a tendency of ＂dispersion--relative concentration＂, this shows that before the earthquake, the entropy of the underground density field was decreased. In addition, dramatic density changes often occur in the Longmenshan fault zone and western Sichuan plateau. Also, with the increase of depth, the trend of density change is more and more obvious. Through comparative analysis, the influence of density change on gravity is much bigger than that from height change.
