Fast 1-D Velocity Optimization Inversion to 3D Velocity Imaging: A Case Study of Sichuan Maerkang Earthquake Swarm in 2022

Abstract

To obtain an accurate one-dimensional velocity model, we developed the EA_VELEST method based on the evolutionary algorithm and the VELEST program. This method can quickly generate a suitable 1D velocity model and finally input it into the 3D velocity inversion process using the TomoDD method. We adopt TomoDD methods to inverse the high-resolution three-dimension velocity structure and relative earthquake hypocenters for this sequence. This system processing flow was applied to the Sichuan Maerkang earthquake swarm in 2022. By collecting the seismic phase data of the Maerkang area between 1 January 2009 and 15 June 2022, we relocated the historical earthquakes in the area and obtained accurate 3D velocity imaging results. The relocated hypocenters reveal a SE-trending secondary fault, which is located ~5 km NW of the Songgang fault. In the first ten-hour of the sequence, events clearly down-dip migrated toward the SE direction. The inverted velocity structure indicates that the majority of earthquakes during the sequence occurred along the boundaries of the high and low-velocity zones or high and low-VP/VS anomalies. Especially both the two largest earthquakes, MS 5.8 and MS 6.0, occurred at the discontinuities of high and low-velocity zones. The EA_VELEST method proposed in this paper is a novel method that has played a very good enlightenment role in the optimization of the one-dimensional velocity model in geophysics and has certain reference significance. The 3D velocity results obtained in this paper and the analysis of tectonic significance provide a reference for the seismogenic environment of this Maerkang earthquake and the deep 3D velocity of the Ganzi block.