Tensor controlled-source audio-frequency magnetotellurics （CSAMT） can yield information about electric and magnetic fields owing to its multi-transmitter configuration compared with the common scalar CSAMT. The most...Tensor controlled-source audio-frequency magnetotellurics （CSAMT） can yield information about electric and magnetic fields owing to its multi-transmitter configuration compared with the common scalar CSAMT. The most current theories, numerical simulations, and inversion of tensor CSAMT are based on far-field measurements and the assumption that underground media have isotropic resistivity. We adopt a three-dimensional （3D） staggered-grid finite difference numerical simulation method to analyze the resistivity in axial anisotropic and isotropic media. We further adopt the limited-memory Broyden- Fletcher-Goldfarb-Shanno （LBFGS） method to perform 3D tensor CSAMT axial anisotropic inversion. The inversion results suggest that when the underground structure is anisotropic, the isotropic inversion will introduce errors to the interpretation.展开更多
A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric （CSAMT） method is developed to invert data in the entire region （near, transition, and far） and deal with ...A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric （CSAMT） method is developed to invert data in the entire region （near, transition, and far） and deal with the effects of artificial sources. First, a regularization factor is introduced in the 2D magnetic inversion, and the magnetic susceptibility is updated in logarithmic form so that the inversion magnetic susceptibility is always positive. Second, the joint inversion of the CSAMT and magnetic methods is completed with the introduction of the cross gradient. By searching for the weight of the cross-gradient term in the objective function, the mutual influence between two different physical properties at different locations are avoided. Model tests show that the joint inversion based on cross-gradient theory offers better results than the single-method inversion. The 2D forward and inverse algorithm for CSAMT with source can effectively deal with artificial sources and ensures the reliability of the final joint inversion algorithm.展开更多
The practical application of 3D inversion of gravity data requires a lot of computation time and storage space.To solve this problem,we present an integrated optimization algorithm with the following components:(1)tar...The practical application of 3D inversion of gravity data requires a lot of computation time and storage space.To solve this problem,we present an integrated optimization algorithm with the following components:(1)targeting high accuracy in the space domain and fast computation in the wavenumber domain,we design a fast 3D forward algorithm with high precision;and(2)taking advantage of the symmetry of the inversion matrix,the main calculation in gravity conjugate gradient inversion is decomposed into two forward calculations,thus optimizing the computational efficiency of 3D gravity inversion.We verify the calculation accuracy and efficiency of the optimization algorithm by testing various grid-number models through numerical simulation experiments.展开更多
Based on an analytical solution for the current point source in an anisotropic half-space,we study the apparent resistivity and apparent chargeability of a transversely isotropic medium with vertical and horizontal ax...Based on an analytical solution for the current point source in an anisotropic half-space,we study the apparent resistivity and apparent chargeability of a transversely isotropic medium with vertical and horizontal axes symmetry,respectively.We then provide a simple derivation of the anisotropy paradoxes in direct current resistivity and time-domain induced polarization methods.Analogous to the mean resistivity,we propose a formulation for deriving the mean polarizability.We also present a three-dimensional finite element algorithm for modeling the direct current resistivity and time-domain induced polarization using an unstructured tetrahedral grid.Finally,we provide the apparent resistivity and apparent chargeability curves of a tilted,transversely isotropic medium with diff erent angles,respectively.The subsequent results illustrate the anisotropy paradoxes of direct current resistivity and time-domain induced polarization.展开更多
The double-phase-shift filtering method,which is based on the traditional purephase-shift filtering method,is a novel approach to harmonic elimination that can be applied to more complicated signals such as white nois...The double-phase-shift filtering method,which is based on the traditional purephase-shift filtering method,is a novel approach to harmonic elimination that can be applied to more complicated signals such as white noise and slip-sweep.Nonetheless,any type of phase-shift filtering method necessitates a relationship between the frequency of fundamental sweep and time,which may cost necessitate an enormous amount of human and physical resources to achieve inaccurate results with low efficiency.This paper combines deep learning with harmonic elimination to produce a double-phase-shift filtering method based on AR2UNet,a type of U-Net with attention gates structure and recurrent residual blocks for improving accuracy and function while simplifying computational complexity.The input of the AR2UNet structure in this paper is seismic data of slip-sweep signals in vibroseis,and the output is signal frequency variation with the time of the fundamental waves,which are required to eliminate the harmonic waves and adjacent signals using a double-phase-shift method to obtain the fundamental sweep.The training sets and test sets are formed by forward models,and a Log-Cosh loss function is used to monitor the process,during which the results of AR2U-Net and traditional U-Net are compared to demonstrate the eminent function of AR2UNet.Following that,the outcomes’Log-Cosh loss functions and accuracy are also compared to validate the conclusion.AR2U-Net,when applied to raw data and combined with the doublephase-shift method,tends to polish the filtering effects and is worth promoting.展开更多
基金sponsored by National Natural Science Foundation of China(No.41374078)
文摘Tensor controlled-source audio-frequency magnetotellurics （CSAMT） can yield information about electric and magnetic fields owing to its multi-transmitter configuration compared with the common scalar CSAMT. The most current theories, numerical simulations, and inversion of tensor CSAMT are based on far-field measurements and the assumption that underground media have isotropic resistivity. We adopt a three-dimensional （3D） staggered-grid finite difference numerical simulation method to analyze the resistivity in axial anisotropic and isotropic media. We further adopt the limited-memory Broyden- Fletcher-Goldfarb-Shanno （LBFGS） method to perform 3D tensor CSAMT axial anisotropic inversion. The inversion results suggest that when the underground structure is anisotropic, the isotropic inversion will introduce errors to the interpretation.
基金jointly sponsored by the Fundamental Research Funds for the Central Universitiesthe National Natural Science Foundation of China(No.41374078)
文摘A two-dimensional forward and backward algorithm for the controlled-source audio-frequency magnetotelluric （CSAMT） method is developed to invert data in the entire region （near, transition, and far） and deal with the effects of artificial sources. First, a regularization factor is introduced in the 2D magnetic inversion, and the magnetic susceptibility is updated in logarithmic form so that the inversion magnetic susceptibility is always positive. Second, the joint inversion of the CSAMT and magnetic methods is completed with the introduction of the cross gradient. By searching for the weight of the cross-gradient term in the objective function, the mutual influence between two different physical properties at different locations are avoided. Model tests show that the joint inversion based on cross-gradient theory offers better results than the single-method inversion. The 2D forward and inverse algorithm for CSAMT with source can effectively deal with artificial sources and ensures the reliability of the final joint inversion algorithm.
基金Financial support by the China Geological Survey Project(Nos.DD20190030,DD20190032)
文摘The practical application of 3D inversion of gravity data requires a lot of computation time and storage space.To solve this problem,we present an integrated optimization algorithm with the following components:(1)targeting high accuracy in the space domain and fast computation in the wavenumber domain,we design a fast 3D forward algorithm with high precision;and(2)taking advantage of the symmetry of the inversion matrix,the main calculation in gravity conjugate gradient inversion is decomposed into two forward calculations,thus optimizing the computational efficiency of 3D gravity inversion.We verify the calculation accuracy and efficiency of the optimization algorithm by testing various grid-number models through numerical simulation experiments.
基金the special funding of Guiyang science and technology bureau and Guiyang University[GYUKY-[2021]]the National Key Research and Development Program of China-Geophysical Comprehensive Exploration and Information Extraction of Deep Mineral Resources(2016YFC0600505)the National K&D Program(2018YFC1504901,2018YFC1504904).
文摘Based on an analytical solution for the current point source in an anisotropic half-space,we study the apparent resistivity and apparent chargeability of a transversely isotropic medium with vertical and horizontal axes symmetry,respectively.We then provide a simple derivation of the anisotropy paradoxes in direct current resistivity and time-domain induced polarization methods.Analogous to the mean resistivity,we propose a formulation for deriving the mean polarizability.We also present a three-dimensional finite element algorithm for modeling the direct current resistivity and time-domain induced polarization using an unstructured tetrahedral grid.Finally,we provide the apparent resistivity and apparent chargeability curves of a tilted,transversely isotropic medium with diff erent angles,respectively.The subsequent results illustrate the anisotropy paradoxes of direct current resistivity and time-domain induced polarization.
基金supported by the National Science and Technology Major Project of China(No.2016ZX05003-003).
文摘The double-phase-shift filtering method,which is based on the traditional purephase-shift filtering method,is a novel approach to harmonic elimination that can be applied to more complicated signals such as white noise and slip-sweep.Nonetheless,any type of phase-shift filtering method necessitates a relationship between the frequency of fundamental sweep and time,which may cost necessitate an enormous amount of human and physical resources to achieve inaccurate results with low efficiency.This paper combines deep learning with harmonic elimination to produce a double-phase-shift filtering method based on AR2UNet,a type of U-Net with attention gates structure and recurrent residual blocks for improving accuracy and function while simplifying computational complexity.The input of the AR2UNet structure in this paper is seismic data of slip-sweep signals in vibroseis,and the output is signal frequency variation with the time of the fundamental waves,which are required to eliminate the harmonic waves and adjacent signals using a double-phase-shift method to obtain the fundamental sweep.The training sets and test sets are formed by forward models,and a Log-Cosh loss function is used to monitor the process,during which the results of AR2U-Net and traditional U-Net are compared to demonstrate the eminent function of AR2UNet.Following that,the outcomes’Log-Cosh loss functions and accuracy are also compared to validate the conclusion.AR2U-Net,when applied to raw data and combined with the doublephase-shift method,tends to polish the filtering effects and is worth promoting.
