Based on the analysis of the conjugate gradient algorithm, we implement a threedimensional （3D） conjugate gradient inversion algorithm with magnetotelluric impedance data. During the inversion process, the 3D conjug...Based on the analysis of the conjugate gradient algorithm, we implement a threedimensional （3D） conjugate gradient inversion algorithm with magnetotelluric impedance data. During the inversion process, the 3D conjugate gradient inversion algorithm doesn＇ t need to compute and store the Jacobian matrix but directly updates the model from the computation of the Jacobian matrix. Requiring only one forward and four pseudo-forward modeling applications per frequency to produce the model update at each iteration, this algorithm efficiently reduces the computation of the inversion. From a trial inversion with synthetic magnetotelluric data, the validity and stability of the 3D conjugate gradient inversion algorithm is verified.展开更多
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.展开更多
Computation of impedance tensor elements is one of the important steps in magnetotelluric data processing. Conventionally, the impedance tensor is defined as a 2 x 2 matrix with Zxx, Zxy, Zyx, and Zyy as elements. In ...Computation of impedance tensor elements is one of the important steps in magnetotelluric data processing. Conventionally, the impedance tensor is defined as a 2 x 2 matrix with Zxx, Zxy, Zyx, and Zyy as elements. In the present study, the six-element impedance tensor is computed with a 2 × 3 matrix using Zxx, Zxy, Zyx, Zyy, Zxz, and Zyz. The properties of the impedance tensor elements have been analyzed for these above two types. The methodology has been tested with five component magnetotelluric data from the Kutch sedimentary basin, Gujarat, India. From the computation of apparent resistivity computation and phase we observed that there is small difference between the four and six impedance elements of Zxy and Zyx for most of the frequency band. However for longer period data, more than 100 sec, an increase in the apparent resistivity and decrease in the phase is observed. We also note that the tipper magnitude is nearly zero for most of the periods, but gradually shows an increasing trend for longer periods （〉100 see）. The Kutch sedimentary basin geoeleetric section shows near horizontal layers at shallow depths and anomalous high conductivity heterogeneous layers at deeper depths may be responsible for the large Hz component at longer periods. This indicates that the vertical component of the magnetic field, Hz, does play an important role in the estimation of electric field parameters in the region with large 2D/3D structures.展开更多
基金sponsored by National Natural Science Foundation of China (Grant Nos. 40774029, 40674037, and 40374024)the National Hi-tech Research and Development Program of China (863 Program) (No. 2007AA09Z310)the Program for New Century Excellent Talents in University (NCET).
文摘Based on the analysis of the conjugate gradient algorithm, we implement a threedimensional （3D） conjugate gradient inversion algorithm with magnetotelluric impedance data. During the inversion process, the 3D conjugate gradient inversion algorithm doesn＇ t need to compute and store the Jacobian matrix but directly updates the model from the computation of the Jacobian matrix. Requiring only one forward and four pseudo-forward modeling applications per frequency to produce the model update at each iteration, this algorithm efficiently reduces the computation of the inversion. From a trial inversion with synthetic magnetotelluric data, the validity and stability of the 3D conjugate gradient inversion algorithm is verified.
基金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.
文摘Computation of impedance tensor elements is one of the important steps in magnetotelluric data processing. Conventionally, the impedance tensor is defined as a 2 x 2 matrix with Zxx, Zxy, Zyx, and Zyy as elements. In the present study, the six-element impedance tensor is computed with a 2 × 3 matrix using Zxx, Zxy, Zyx, Zyy, Zxz, and Zyz. The properties of the impedance tensor elements have been analyzed for these above two types. The methodology has been tested with five component magnetotelluric data from the Kutch sedimentary basin, Gujarat, India. From the computation of apparent resistivity computation and phase we observed that there is small difference between the four and six impedance elements of Zxy and Zyx for most of the frequency band. However for longer period data, more than 100 sec, an increase in the apparent resistivity and decrease in the phase is observed. We also note that the tipper magnitude is nearly zero for most of the periods, but gradually shows an increasing trend for longer periods （〉100 see）. The Kutch sedimentary basin geoeleetric section shows near horizontal layers at shallow depths and anomalous high conductivity heterogeneous layers at deeper depths may be responsible for the large Hz component at longer periods. This indicates that the vertical component of the magnetic field, Hz, does play an important role in the estimation of electric field parameters in the region with large 2D/3D structures.
