3-D finite-volume forward modeling of wide-field EM using scattered potentials

Peng RongHua,Hu XiangYun,Li JianHui,Hu ZuZhi,Yang Hui
DOI: https://doi.org/10.6038/cjg2018L0363
2018-01-01
Abstract:Controlled-source audio-frequency magnetotelluric (CSAMT) method has developed into an established technique and been extensively used for a broad range of applications in geophysical exploration and near-surface investigations over the past decades. Conventionally, CSAMT measurements are carried out in the far-field zones to meet the plane wave approximation of CSAMT data, where the separation between transmitters and receivers should be sufficiently large. However, near- and intermediate-field CSAMT data are often collected due to improper employment of transmitter-receiver separation and limitations of logistics access in some prospecting areas such as mountainous regions. Under such circumstance, Cagniard apparent resistivity defined for magnetotelluric plane waves is no longer reasonable in these areas. In order to overcome the limitation of far-field measurements, wide field electromagnetic method (WFEM) introduces the wide field apparent resistivity which is defined from the analytical expressions of excited electromagnetic fields for homogeneous halfspace. The use of wide field apparent resistivity instead of Cagniard apparent resistivity greatly improves measurable regions and reduces distortion for data collected in near- and intermediate-field zones. Considering both the three dimensional (3D) nature of the excited electromagnetic fields and complex subsurface structures, in this paper, we present an efficient finite volume 3D forward modeling algorithm for frequency-domain controlled-source electromagnetic simulations. To circumvent ill-conditioning of the system of linear equations resulting from discretization, Maxwell's equations are reformulated in terms of coupled scalar and vector potentials under Coulomb gauge condition, and an explicit divergence is enforced to avoid non-physical solutions introduced during solution process. In addition, the well-known scattering field formulation is employed to mitigate the numerical singularity in the vicinity of the impressed current source. The accuracy of the algorithm is validated against quasi-static solutions for 1D layered model. Besides, comparisons between CSAMT and WFEM measurements for typical 3D conductive structure under the same survey configuration demonstrate that the wide field apparent resistivity is less distorted in the non-plane wave zones and reveals better resolution capability than Cagniard apparent resistivity.
What problem does this paper attempt to address?