Feed aggregator

Combined Boundary Element and Finite Volume Methods for Modeling Fluid-Induced Seismicity in Fault Networks within Low-Permeability Rocks

Geophysical Journal International - Mon, 10/06/2025 - 00:00

SummaryTo better understand the mechanics of injection-induced seismicity, we developed a two-dimensional numerical code to simulate both seismic and aseismic slip on non-planar faults and fault networks driven by fluid diffusion along permeable faults, in an impervious host rock. Our approach integrates a boundary element method to model fault slip governed by rate-and-state friction with a finite-volume method to simulate fluid diffusion along fault networks. We demonstrate the capabilities of the method with two illustrative examples: (1) fluid injection inducing slow slip on a primary rough, rate-strengthening fault, which subsequently triggers microseismicity on nearby secondary, smaller faults, and (2) fluid injection on a single fault in a network of intersecting faults, leading to fluid diffusion and reactivation of slip throughout the network. This work highlights the importance of distinguishing between mechanical and hydrological processes in the analysis of induced seismicity, providing a powerful tool for improving our understanding of fault behavior in response to fluid injection, in particular when a network of geometrically complex faults is involved.

3D joint inversion of induced polarization and self-potential data for ore body localization

Geophysical Journal International - Mon, 10/06/2025 - 00:00

SUMMARYIn mineral exploration, induced polarization and self-potential are two broadly used active and passive geophysical methods, respectively. In the case of ore bodies, both methods are associated with charge distributions associated with a secondary electrical field (induced polarization) and a source current density (self-potential). Both the chargeability and volumetric source current density distributions bring information regarding the shape of ore bodies. Therefore the joint inversion of these datasets is expected to better tomograms of ore bodies. A joint inversion approach is developed to combine both methods. The objective function to minimize includes two independent components plus a cross-gradient joint function. The use of the cross-gradient is justified from the underlying physics of the two geophysical problems at play. The structure of the cost function is tailored to overcome some problems like convergence and parameter determination in the inverse process. Two synthetic tests and a laboratory experiment are used to benchmark the proposed algorithm. We demonstrate that the joint inversion algorithm performs better than the localizations obtained from independent inversion approaches. To refine the interpretation of the shape of ores, we introduce an ore presence index using the chargeability and source current density resulting from the joint inversion algorithm. The K-Medoids clustering algorithm is used to automatically categorize the calculated ore presence index into different clusters. The cluster with larger values successfully identifies the ore bodies associated with strong chargeability and/or volumetric source current density.