Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor...Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor CSAMT forward modeling using the vector finite-element method. To verify the feasibility of the algorithm, we calculate the electric field, magnetic field, and tensor impedance of the 3D CSAMT far-zone field in layered media and compare them with theoretical solutions. In addition, a three-dimensional anomaly in half-space is also simulated, and the response characteristics of the impedance tensor and the apparent resistivity and impedance phase are analyzed. The results suggest that the vector finite-element method produces high-precision electromagnetic field and impedance tensor data, satisfies the electric field discontinuity, and does not require divergence correction using the vector finite-element method.展开更多
How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem ...How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem was obtained by the regularization methods in which some desired structures were imposed to stabilize the inverse problem. By the smoothness-constrained model and approximate sensitivity method, the stable subsurface resistivity structures were reconstructed. The synthetic examples show that the smoothness-constrained regularized inversion method is effective and can be reasonable to reconstruct three-dimensional subsurface resistivity structures.展开更多
Geophysical inversion under different stabilizers has different descriptions of the target body boundary,especially in complex geological structures.In this paper,we present an extremum boundary inversion algorithm ba...Geophysical inversion under different stabilizers has different descriptions of the target body boundary,especially in complex geological structures.In this paper,we present an extremum boundary inversion algorithm based on different stabilizers for electrical interface recognition.Firstly,we use the smoothest and minimum-support stabilizing functional to study the applicability of adaptive regularization inversion algorithm.Then,an electrical interface recognition method based on different stabilizers is developed by introducing extremum boundary inversion algorithm.The testing shows that the adaptive regularization inversion method does work for different stabilizers and has a low dependence on the initial models.The ratio of the smooth and focusing upper and lower boundaries obtained using the extremum boundary inversion algorithm can clearly demarcate electrical interfaces.We apply the inversion algorithm to the magnetotelluric(MT)data collected from a preselected area of a high-level-waste clay-rock repository site in the Tamusu area.We recognized regional structures with smooth inversion and the local details with focusing inversion and determined the thickness of the target layer combined with the geological and drilling information,which meets the requirement for the site of the high-level waste clay-rock repository.展开更多
基金supported by the National Natural Science Foundation of China(No.41104068)the Deep Exploration in China,Sino Probe-03-05
文摘Scalar CSAMT is only suitable for measurements in one and two dimensions perpendicular to geological structures. For complex 3D geoelectric structure, tensor CSAMT is more suitable. In this paper, we discuss 3D tensor CSAMT forward modeling using the vector finite-element method. To verify the feasibility of the algorithm, we calculate the electric field, magnetic field, and tensor impedance of the 3D CSAMT far-zone field in layered media and compare them with theoretical solutions. In addition, a three-dimensional anomaly in half-space is also simulated, and the response characteristics of the impedance tensor and the apparent resistivity and impedance phase are analyzed. The results suggest that the vector finite-element method produces high-precision electromagnetic field and impedance tensor data, satisfies the electric field discontinuity, and does not require divergence correction using the vector finite-element method.
基金Project(20110162120064)supported by Higher School Doctor Subject Special Scientific Research Foundation of ChinaProject(10JJ6059)supported by the Natural Science Foundation of Hunan Province,China
文摘How to get the rapid and stable inversion results and reconstruct the clear subsurface resistivity structures is a focus problem in current magnetotelluric inversion. A stable solution of an ill-posed inverse problem was obtained by the regularization methods in which some desired structures were imposed to stabilize the inverse problem. By the smoothness-constrained model and approximate sensitivity method, the stable subsurface resistivity structures were reconstructed. The synthetic examples show that the smoothness-constrained regularized inversion method is effective and can be reasonable to reconstruct three-dimensional subsurface resistivity structures.
基金supported by the National Natural Science Foundation of China(Nos.41604104,41674077 and 41404057)PRC High-level Radioactive Waste Geological Disposal Project([2014] No.1578)+2 种基金Open Fund of State Key Laboratory of Marine Geology(Tongji University)(MGK1704)Jiangxi Province Youth Science Fund(No.20171BAB213031)Scientific Research Starting Foundation for Doctors of East China University of Technology(DHBK201403)
文摘Geophysical inversion under different stabilizers has different descriptions of the target body boundary,especially in complex geological structures.In this paper,we present an extremum boundary inversion algorithm based on different stabilizers for electrical interface recognition.Firstly,we use the smoothest and minimum-support stabilizing functional to study the applicability of adaptive regularization inversion algorithm.Then,an electrical interface recognition method based on different stabilizers is developed by introducing extremum boundary inversion algorithm.The testing shows that the adaptive regularization inversion method does work for different stabilizers and has a low dependence on the initial models.The ratio of the smooth and focusing upper and lower boundaries obtained using the extremum boundary inversion algorithm can clearly demarcate electrical interfaces.We apply the inversion algorithm to the magnetotelluric(MT)data collected from a preselected area of a high-level-waste clay-rock repository site in the Tamusu area.We recognized regional structures with smooth inversion and the local details with focusing inversion and determined the thickness of the target layer combined with the geological and drilling information,which meets the requirement for the site of the high-level waste clay-rock repository.
