Rongjiang Tang

and 3 more

The 1D inversion method is widely used in airborne transient electromagnetic method (ATEM) but often encounter inaccurate imaging and overlook important geological information. The three-dimensional inversion is relatively robust and provide more detailed information on subsurface resistivity distributions but limited by significant computational burdens. In this study, we develop a rapid 3D imaging scheme based on deep learning for subsurface resistivity structures with lateral heterogeneity, which performs better than traditional 1D inversion but slightly inferior to 3D inversion, but its efficiency makes it more practical than traditional 3D inversion. To create a reasonably large training set of different resistivity models that generalizes well. Here we propose three strategies that focus on training datasets, to improve the performance of DL models, including divide-and-conquer strategy, random models generating and depth constraints. Through a large of reasonable structural models, appropriate network setups, a generalized result can be obtained through our proposed UNet framework, which has been demonstrated to be effective on both synthetic and field data. The results of gradient-based 3D and 1D inversion method are compared and analyzed, respectively, which demonstrates the reliability of the proposed method. To flexibly apply the deep learning model to field data, we propose a model scaling scheme to make pre-trained UNet model compatible with different inversion specifications and survey configurations, without retraining a new DL model. The inverted structures from field data clearly delineate the geometries of the lakes, faults and surrounding mountains. The inversion operator can support instantaneous 3D subsurface resistivity imaging for ATEM observations.

Rongjiang Tang

and 1 more

The 2008 Wenchuan Mw 7.9 mainshock has caused catastrophic destruction to cities along the northwestern margin of the Sichuan Basin. The Wenchuan-Maoxian Fault (WMF) on the hinterland side, along with a conjugate buried Lixian fault (LXF) was not activated by this earthquake but is likely to experience large earthquakes in the future. We perform 3D dynamic earthquake rupture simulations on the WMF and LXF to access the possibility of the earthquake occurrence and further explore the possible size of earthquakes and the distribution of high seismic risk in the future. We firstly invert focal mechanism solutions to get a heterogeneous tectonic stress field as the initial stress of simulation. Then we develop a new method to refine fault geometry through inverting long-term slip rates. Several fault nucleation points, friction coefficients, and initial stress states are tested, and the general rupture patterns for these earthquake scenarios are evaluated and could fall into three groups. Depending on initial conditions, the dynamic rupture may start in the LXF, leading to magnitude-7.0 earthquakes, or start in the WMF, then cascades through the LXF, leading to magnitude-7.5 earthquakes, or both start and arrest in the WMF, leading to around magnitude-6.5 or 7.0 earthquakes. We find that the rupture starting on the reverse oblique-slip tends to jump to the strike-slip fault, but the reverse process is suppressed. The rupture propagating eastward causes larger coseismic displacements than the westward propagation, and relatively high peak ground velocities are distributed near the northeastern end of WMF.