Full-Waveform Modeling of Complex Media Seismic Waves for Irregular Topography and Its Application in Metal Ore Exploration

Wenchao Su,Shoudong Huo,Xuhui Zhou
DOI: https://doi.org/10.3390/min14070664
IF: 2.5
2024-06-28
Minerals
Abstract:Seismic exploration has caught widespread attention in metal ore exploration due to its higher resolution. However, the presence of topography and complex underground structures in metal ore exploration complicates seismic records. Therefore, it is essential to apply a numerical simulation method suitable for metal ore exploration to study the propagation law of seismic waves in shallow and ore-forming zones, providing reliable theoretical support for multi-component seismic techniques. In particular, the presence of topography generates strong-amplitude surface waves, scattered waves, and converted waves, which consistently distort seismic records and affect the imaging accuracy of the metallogenic belts. Additionally, the propagation of seismic waves is also affected by the anisotropy and viscoelasticity of the underground medium. This paper proposes an elastic wave finite-difference numerical simulation method suitable for irregularly topographical and complex medium conditions, named the comprehensive parameter correction method, which implements a free-surface boundary condition based on the concept of medium averaging. It is algorithmically simple and implies no additional computational costs. Meanwhile, the results obtained by this method are highly consistent with those of the spectral element method, demonstrating its accuracy. By presenting several numerical simulation cases and illustrating the impact of topography and medium conditions on seismic records, this paper demonstrates the necessity of considering irregularly topographical and complex medium conditions in metal ore exploration. In conclusion, the numerical simulation method we propose provides a solid theoretical foundation for the application of seismic exploration methods in metal ore exploration.
geochemistry & geophysics,mineralogy,mining & mineral processing
What problem does this paper attempt to address?
This paper attempts to solve the problem of the influence of complex terrain and underground medium conditions on seismic wave propagation in metal mineral exploration. Specifically, the paper focuses on: 1. **Distortion of seismic records by complex terrain**: In metal mineral exploration, the complexity of the terrain can lead to the generation of surface waves, scattered waves and converted waves, which can significantly distort seismic records and affect the imaging accuracy. 2. **Influence of complex underground medium conditions**: In addition to terrain factors, the anisotropy and viscoelasticity of the underground medium can also affect the propagation of seismic waves. Anisotropy makes the propagation speed of seismic waves different in different directions, while viscoelasticity will cause the attenuation of seismic waves, manifested as changes in amplitude and phase. To solve these problems, the paper proposes an elastic - wave finite - difference numerical simulation method suitable for irregular terrain and complex medium conditions - the **Comprehensive Parameter Correction Method (CPCM)**. By implementing the free - surface boundary condition, combined with the step - approximation method and the constitutive relation and density correction at the free surface, this method can effectively handle the numerical simulation problem under irregular terrain conditions. In addition, for cases involving anisotropic, viscoelastic or anisotropic - viscoelastic media, the paper also derives the corresponding mathematical expressions of the free - surface boundary condition, and verifies the feasibility and accuracy of CPCM through numerical experiments. Finally, the paper shows the application potential of this method in actual metal mineral exploration. In particular, the test on the Half Mile Lake deposit model shows that CPCM has good performance in dealing with practical problems. ### Key formulas 1. **Constitutive relation correction under free - surface boundary conditions**: \[ \begin{bmatrix} \sigma_{xx} \\ \sigma_{zz} \\ \sigma_{xz} \end{bmatrix} = \begin{bmatrix} \frac{2\mu(\lambda + \mu)}{\lambda + 2\mu} & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{bmatrix} \begin{bmatrix} \varepsilon_{xx} \\ \varepsilon_{zz} \\ \varepsilon_{xz} \end{bmatrix} \] where \(\sigma\) represents the stress tensor, \(\varepsilon\) represents the strain tensor, and \(\lambda\) and \(\mu\) are Lame parameters. 2. **Constitutive relation correction in anisotropic media** (taking VTI media as an example): \[ \begin{bmatrix} \sigma_{xx} \\ \sigma_{zz} \\ \sigma_{xz} \end{bmatrix} = \begin{bmatrix} \frac{c_{11}c_{33}-c_{13}^2}{2c_{33}} & 0 & 0 \\ 0 & 0 & 0 \\ 0 & 0 & 0 \end{bmatrix} \begin{bmatrix} \varepsilon_{xx} \\ \varepsilon_{zz} \\ \varepsilon_{zx} \end{bmatrix} \] 3. **Constitutive relation correction in viscoelastic media**: \[ \begin{bmatrix} \sigma_{xx} \\ \sigma_{zz} \\ \sigma_{xz} \end{bmatrix} = \begin{bmatrix} \end{bmatrix} \]