The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy ...The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.展开更多
2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcani...2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcanic region. It is used in this study that seismic traveltime inversion for simultaneous determination of 2-D velocity and interface structure of the crust and upper mantle. The result shows that, under Changbaishan-Tianchi crater, there exists a low-velocity body in the shape of an inverted triangle, and the crustal reflecting boundaries and Moho all become lower by a varying margin of 2-6 km, forming a crustal root which is assumed to be the Changbaishan-Tianchi volcanic system. Finally, we make a comparison between our 2-D velocity model and the result from the studies by using trial-and-error forward modeling with SEIS83.展开更多
In this study,a deep learning algorithm was applied to two-dimensional magnetotelluric(MT)data inversion.Compared with the traditional linear iterative inversion methods,the MT inversion method based on convolutional ...In this study,a deep learning algorithm was applied to two-dimensional magnetotelluric(MT)data inversion.Compared with the traditional linear iterative inversion methods,the MT inversion method based on convolutional neural networks(CNN)does not rely on the selection of the initial model parameters and does not fall into the local optima.Although the CNN inversion models can provide a clear electrical interface division,their inversion results may remain prone to abrupt electrical interfaces as opposed to the actual underground electrical situation.To solve this issue,a neural network with a residual network architecture(ResNet-50)was constructed in this study.With the apparent resistivity and phase pseudo-section data as the inputs and with the resistivity parameters of the geoelectric model as the training labels,the modified ResNet-50 model was trained end-to-end for producing samples according to the corresponding production strategy of the study area.Through experiments,the training of the ResNet-50 with the dice loss function effectively solved the issue of over-segmentation of the electrical interface by the cross-entropy function,avoided its abrupt inversion,and overcame the computational inefficiency of the traditional iterative methods.The proposed algorithm was validated against MT data measured from a geothermal field prospect in Huanggang,Hubei Province,which showed that the deep learning method has opened up broad prospects in the field of MT data inversion.展开更多
In the paper, a 2D symmetrical anisotropic medium whose strike agrees with one of the horizontal principal axes is considered to develop a corresponding inversion technique. In the specified conditions, if we assume a...In the paper, a 2D symmetrical anisotropic medium whose strike agrees with one of the horizontal principal axes is considered to develop a corresponding inversion technique. In the specified conditions, if we assume an equivalent conductivity anisotropy in both the vertical and dipping directions, i.e., σzz=σyy, the differential equations obtained are formally the same as that for TE and TM modes in the 2D isotropic geoelectrical media. The same inversion technique as that in the 2D isotropic media can be employed to obtain the anisotropic conductivities. It means that the TE and TM inversion results in the isotropic media can be respectively thought as the resistivities in the two principal directions of the symmetrically anisotropic media, which has offered a new approach and a theoretical guidance for interpreting magnetotelluric data. And the inversion technique developed here is used to test the magnetotelluric data in the area of Tianzhu and Yongdeng in Gansu Province, so that the crust anisotropic geoelectrical structures in this region can be obtained.展开更多
Based on the fact that it is diffi cult to implement optimum inversion using 2D and 3D forward modeling with magnetic-source transient electromagnetics(TEM),this paper explores a novel approach to the implementation o...Based on the fact that it is diffi cult to implement optimum inversion using 2D and 3D forward modeling with magnetic-source transient electromagnetics(TEM),this paper explores a novel approach to the implementation of 2D magnetic-source TEM inversion.In particular,we converted magnetic-source TEM data into magnetotelluric(MT)data and then used a 2D MT inversion method to implement a 2D magnetic-source TEM inversion interpretation.First,we studied the similarity between magnetic-source TEM waves and MT waves and between magnetic-source TEM all-time apparent resistivity and MT Cagniard apparent resistivity.Then,we selected an optimal time-frequency transformation coeffi cient to implement rapid time-frequency transformation of all-time TEM apparent resistivity to MT Cagniard apparent resistivity.Afterward,we conducted 1D pseudo-MT inversions of magnetic-source 1D TEM theoretical models.The 1D inversion results demonstrated that the diff erence between the inversion parameters and model parameters was small,while the MT 1D inversion method could be used to conduct magnetic 1D TEM inversion within a certain margin of error.We further conducted 2D pseudo-MT inversions of 3D magnetic-source TEM theoretical models,and the 2D inversion results indicated that selecting a joint 2D pseudo-MT transverse-electric(TE)and transverse-magnetic(TM)inversion method based on measuring the line above a 3D anomalous body can help to accurately implement a 2D inversion interpretation of the 3D TEM response.展开更多
This work essentially consists in inverting in an exact, explicit, and original way the pentadiagonal Toeplitz matrix or tridiagonal block matrix resulting from the discretization of the two-dimensional Laplace operat...This work essentially consists in inverting in an exact, explicit, and original way the pentadiagonal Toeplitz matrix or tridiagonal block matrix resulting from the discretization of the two-dimensional Laplace operator. This method is an algorithm facilitating the resolution of a large number of problems governed by PDEs involving the Laplacian in two dimensions. It guarantees high precision and high efficiency in solving various differential equations.展开更多
基金Projects(41074085,41374118)supported by the National Natural Science Foundation of ChinaProject(20120162110015)supported by Doctoral Fund of Ministry of Education of ChinaProject(NCET-12-0551)supported by Program for New Century Excellent Talents in University,China
文摘The 2D data processing adopted by the high-density resistivity method regards the geological structures as two degrees, which makes the results of the 2D data inversion only an approximate interpretation;the accuracy and effect can not meet the precise requirement of the inversion. Two typical models of the geological bodies were designed, and forward calculation was carried out using finite element method. The forward-modeled profiles were obtained. 1% Gaussian random error was added in the forward models and then 2D and 3D inversions using a high-density resistivity method were undertaken to realistically simulate field data and analyze the sensitivity of the 2D and 3D inversion algorithms to noise. Contrast between the 2D and 3D inversion results of least squares inversion shows that two inversion results of high-density resistivity method all can basically reflect the spatial position of an anomalous body. However, the 3D inversion can more effectively eliminate the influence of interference from Gaussian random error and better reflect the distribution of resistivity in the anomalous bodies. Overall, the 3D inversion was better than 2D inversion in terms of embodying anomalous body positions, morphology and resistivity properties.
基金Key Project (95-11-02-01) from China Seismological Bureau.Contribution No. RCEG200129, Research Center of Exploration Geophysi
文摘2-D velocity structure and tectonics of the crust and upper mantle is revealed by inversion of seismic refraction and wide-angle reflection traveltimes acquired along the profile L1 in the Changbaishan-Tianchi volcanic region. It is used in this study that seismic traveltime inversion for simultaneous determination of 2-D velocity and interface structure of the crust and upper mantle. The result shows that, under Changbaishan-Tianchi crater, there exists a low-velocity body in the shape of an inverted triangle, and the crustal reflecting boundaries and Moho all become lower by a varying margin of 2-6 km, forming a crustal root which is assumed to be the Changbaishan-Tianchi volcanic system. Finally, we make a comparison between our 2-D velocity model and the result from the studies by using trial-and-error forward modeling with SEIS83.
基金co-funded by the National Natural Science Foundation of China(No.42220104002,42174095,U1812402,and 41630317)the Opening Fund of Key Laboratory of Geological Survey and Evaluation of Ministry of Education(No.GLAB2022ZR10)the Fundamental Research Funds for the Central Universities,China University of Geosciences(Wuhan).
文摘In this study,a deep learning algorithm was applied to two-dimensional magnetotelluric(MT)data inversion.Compared with the traditional linear iterative inversion methods,the MT inversion method based on convolutional neural networks(CNN)does not rely on the selection of the initial model parameters and does not fall into the local optima.Although the CNN inversion models can provide a clear electrical interface division,their inversion results may remain prone to abrupt electrical interfaces as opposed to the actual underground electrical situation.To solve this issue,a neural network with a residual network architecture(ResNet-50)was constructed in this study.With the apparent resistivity and phase pseudo-section data as the inputs and with the resistivity parameters of the geoelectric model as the training labels,the modified ResNet-50 model was trained end-to-end for producing samples according to the corresponding production strategy of the study area.Through experiments,the training of the ResNet-50 with the dice loss function effectively solved the issue of over-segmentation of the electrical interface by the cross-entropy function,avoided its abrupt inversion,and overcame the computational inefficiency of the traditional iterative methods.The proposed algorithm was validated against MT data measured from a geothermal field prospect in Huanggang,Hubei Province,which showed that the deep learning method has opened up broad prospects in the field of MT data inversion.
基金National Natural Science Foundation of China (40074010).
文摘In the paper, a 2D symmetrical anisotropic medium whose strike agrees with one of the horizontal principal axes is considered to develop a corresponding inversion technique. In the specified conditions, if we assume an equivalent conductivity anisotropy in both the vertical and dipping directions, i.e., σzz=σyy, the differential equations obtained are formally the same as that for TE and TM modes in the 2D isotropic geoelectrical media. The same inversion technique as that in the 2D isotropic media can be employed to obtain the anisotropic conductivities. It means that the TE and TM inversion results in the isotropic media can be respectively thought as the resistivities in the two principal directions of the symmetrically anisotropic media, which has offered a new approach and a theoretical guidance for interpreting magnetotelluric data. And the inversion technique developed here is used to test the magnetotelluric data in the area of Tianzhu and Yongdeng in Gansu Province, so that the crust anisotropic geoelectrical structures in this region can be obtained.
基金this research project is funded by a major science and technology project of Gansu province,“research on the complete set technology for highway construction in collapsible loess region of Gansu province”(No.1302GKDA009).
文摘Based on the fact that it is diffi cult to implement optimum inversion using 2D and 3D forward modeling with magnetic-source transient electromagnetics(TEM),this paper explores a novel approach to the implementation of 2D magnetic-source TEM inversion.In particular,we converted magnetic-source TEM data into magnetotelluric(MT)data and then used a 2D MT inversion method to implement a 2D magnetic-source TEM inversion interpretation.First,we studied the similarity between magnetic-source TEM waves and MT waves and between magnetic-source TEM all-time apparent resistivity and MT Cagniard apparent resistivity.Then,we selected an optimal time-frequency transformation coeffi cient to implement rapid time-frequency transformation of all-time TEM apparent resistivity to MT Cagniard apparent resistivity.Afterward,we conducted 1D pseudo-MT inversions of magnetic-source 1D TEM theoretical models.The 1D inversion results demonstrated that the diff erence between the inversion parameters and model parameters was small,while the MT 1D inversion method could be used to conduct magnetic 1D TEM inversion within a certain margin of error.We further conducted 2D pseudo-MT inversions of 3D magnetic-source TEM theoretical models,and the 2D inversion results indicated that selecting a joint 2D pseudo-MT transverse-electric(TE)and transverse-magnetic(TM)inversion method based on measuring the line above a 3D anomalous body can help to accurately implement a 2D inversion interpretation of the 3D TEM response.
文摘This work essentially consists in inverting in an exact, explicit, and original way the pentadiagonal Toeplitz matrix or tridiagonal block matrix resulting from the discretization of the two-dimensional Laplace operator. This method is an algorithm facilitating the resolution of a large number of problems governed by PDEs involving the Laplacian in two dimensions. It guarantees high precision and high efficiency in solving various differential equations.
