The Beishan orogen,located in the central segment of the Tianshan–Solonker suture within the southern Central Asian Orogenic Belt(CAOB),is crucial for understanding the accretionary processes and continental growth i...The Beishan orogen,located in the central segment of the Tianshan–Solonker suture within the southern Central Asian Orogenic Belt(CAOB),is crucial for understanding the accretionary processes and continental growth in Central Asia.This orogen developed through the episodic amalgamation and accretion of continental margin arcs,island arcs,ophiolites,and accretionary wedges,undergoing a complex process of accretion and evolution.Since the Phanerozoic,the Beishan orogen has experienced multiple phases of magmatic and collision events.The intricate distribution of magmatic arc rocks has obscured the complete basement traces,and the spatial superposition of multiple magmatic arc phases has complicated the study of its evolutionary history.展开更多
In order to attain good quality transfer function estimates from magnetotelluric field data(i.e.,smooth behavior and small uncertainties across all frequencies),we compare time series data processing with and without ...In order to attain good quality transfer function estimates from magnetotelluric field data(i.e.,smooth behavior and small uncertainties across all frequencies),we compare time series data processing with and without a multitaper approach for spectral estimation.There are several common ways to increase the reliability of the Fourier spectral estimation from experimental(noisy)data;for example to subdivide the experimental time series into segments,taper these segments(using single taper),perform the Fourier transform of the individual segments,and average the resulting spectra.展开更多
In the framework of a mineral system approach,a combination of components is required to develop a mineral system.This includes the whole-lithosphere architecture,which controls the transport of ore-forming fluids,and...In the framework of a mineral system approach,a combination of components is required to develop a mineral system.This includes the whole-lithosphere architecture,which controls the transport of ore-forming fluids,and favorable tectonic and geodynamic processes,occurring at various spatial and temporal scales,that influence the genesis and evolution of ore-forming fluids(Huston et al.,2016;Groves et al.,2018;Davies et al.,2020).Knowledge of the deep structural framework can advance the understanding of the development of a mineral system and the emplacement of mineral deposits.Deep geophysical exploration carried out with this aim is increasingly important for targeting new ore deposits in unexplored and underexplored regions(Dentith et al.,2018;Dentith,2019).展开更多
It is of crucial importance to investigate the spatial structures of ancient landslides in the eastern Tibetan Plateau’s alpine canyons as they could provide valuable insights into the evolutionary history of the lan...It is of crucial importance to investigate the spatial structures of ancient landslides in the eastern Tibetan Plateau’s alpine canyons as they could provide valuable insights into the evolutionary history of the landslides and indicate the potential for future reactivation.This study examines the Deda ancient landslide,situated in the Chalong-ranbu fault zone,where creep deformation suggests a complex underground structure.By integrating remote sensing,field surveys,Audio-frequency Magnetotellurics(AMT),and Microtremor Survey Method(MSM)techniques,along with engineering geological drilling for validation,to uncover the landslide’s spatial feature s.The research indicates that a fault is developed in the upper part of the Deda ancient landslide,and the gully divides it into Deda landslide accumulation zoneⅠand Deda landslide accumulation zoneⅡin space.The distinctive geological characteristics detectable by MSM in the shallow subsurface and by AMT in deeper layers.The findings include the identification of two sliding zones in the Deda I landslide,the shallow sliding zone(DD-I-S1)depth is approximately 20 m,and the deep sliding zone(DD-I-S2)depth is 36.2-49.9 m.The sliding zone(DD-Ⅱ-S1)depth of the DedaⅡlandslide is 37.6-43.1 m.A novel MSM-based method for sliding zone identification is proposed,achieving less than 5%discrepancy in depth determination when compared with drilling data.These results provide a valuable reference for the spatial structural analysis of large-deepseated landslides in geologically complex regions like the eastern Tibetan Plateau.展开更多
Despite extensive efforts to understand the tectonic evolution of the Jiangnan Orogen in South China,the orogenic process and its mechanism remain a matter of dispute.Previous geodynamic studies have mostly focused on...Despite extensive efforts to understand the tectonic evolution of the Jiangnan Orogen in South China,the orogenic process and its mechanism remain a matter of dispute.Previous geodynamic studies have mostly focused on collisional orogeny,which is commonly invoked to explain the Jiangnan Orogen.However,it is difficult for such hypotheses to reconcile all the geological and geophysical data,especially the absence of ultrahigh-pressure metamorphic rocks.Based on the magnetotelluric data,we present a group of resistivity models produced through the combination of two-dimensional and three-dimensional inversions,revealing the geo-electrical structures of Jiangnan and a typical collisional orogen.In our models,the resistive crust is separated into three parts by a prominent conductive layer with opposite dipping directions on both sides.A special thrust-nappe system,which is different from that developed in a typical collisional process,is revealed in the Jiangnan Orogen.This structure suggests a process different from the simple collisional orogeny.To interpret our observations,an'intra-continental orogeny'is proposed to address the development of the Jiangnan Orogen in the Mesozoic.Furthermore,this'reworked'process may contain at least two stages caused by the decoupling of the lithosphere,which is revealed by an extra conductive layer beneath Jiangnan.展开更多
With the increase in the coverage area of magnetotelluric data,three-dimensional magnetotelluric modeling in spherical coordinates and its differences with respect to traditional Cartesian modeling have gradually attr...With the increase in the coverage area of magnetotelluric data,three-dimensional magnetotelluric modeling in spherical coordinates and its differences with respect to traditional Cartesian modeling have gradually attracted attention.To fully understand the influence of the Earth’s curvature and map projection deformations on Cartesian modeling,qualitative and quantitative analyses based on realistic three-dimensional models need to be examined.Combined with five representative map projections,a type of model conversion method that transforms the original spherical electrical conductivity model to Cartesian coordinates is described in this study.The apparent resistivity differences between the spherical western United States electrical conductivity model and the corresponding five Cartesian models are then compared.The results show that the cylindrical equal distance map projection has the smallest error.A meridian convergence correction resulting from the deformation of the map projection is introduced to rotate the Cartesian impedance tensor from grid north to geographic north,which reduces differences from the spherical results.On the basis of the magnetotelluric field data,the applicability of the Cartesian coordinate system to western and contiguous United States models is quantitatively evaluated.Precise interpretations of the contiguous United States model were found to require spherical coordinates.展开更多
We implement a parallel algorithm with the advantage of MPI （Message Passing Interface） to speed up the rapid relaxation inversion for 3D magnetotelluric data. We test the parallel rapid relaxation algorithm with sy...We implement a parallel algorithm with the advantage of MPI （Message Passing Interface） to speed up the rapid relaxation inversion for 3D magnetotelluric data. We test the parallel rapid relaxation algorithm with synthetic and real data. The execution efficiency of the algorithm for several different situations is also compared. The results indicate that the parallel rapid relaxation algorithm for 3D magnetotelluric inversion is effective. This parallel algorithm implemented on a common PC promotes the practical application of 3D magnetotelluric inversion and can be suitable for the other geophysical 3D modeling and inversion.展开更多
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.展开更多
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.展开更多
Different geophysical exploration methods have significant differences in terms of exploration depth,especially in frequency domain electromagnetic(EM)exploration.According to the definition of skin depth,this differe...Different geophysical exploration methods have significant differences in terms of exploration depth,especially in frequency domain electromagnetic(EM)exploration.According to the definition of skin depth,this difference will increase with the effective detection frequency of the method.As a result,when performing three-dimensional inversion on single type of EM data,it is not possible to effectively distinguish the subsurface geoelectric structure at the full scale.Therefore,it is necessary to perform joint inversion on different type of EM data.In this paper we combine the magnetotelluric method(MT)with the controlled-source audio-magnetotelluric method(CSAMT)to study the frequency-domain three-dimensional(3D)joint inversions,and we use the unstructured finite-element method to do the forward modeling for them,so that the numerical simulation accuracies of different electromagnetic methods can be satisfied.By combining the two sets of data,we can obtain the sensitivity of the electrical structure at different depths,and depict the full-scale subsurface geoelectric structures.In actual mineral exploration,the 3D joint inversion is more useful for identifying subsurface veins in the shallow part and blind mines in the deep part.It can delineate the morphological distribution of ore bodies more completely and provide reliable EM interpretations to guide the mining of minerals.展开更多
Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficien...Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.展开更多
The workload of the 3D magnetotelluric forward modeling algorithm is so large that the traditional serial algorithm costs an extremely large compute time. However, the 3D forward modeling algorithm can process the dat...The workload of the 3D magnetotelluric forward modeling algorithm is so large that the traditional serial algorithm costs an extremely large compute time. However, the 3D forward modeling algorithm can process the data in the frequency domain, which is very suitable for parallel computation. With the advantage of MPI and based on an analysis of the flow of the 3D magnetotelluric serial forward algorithm, we suggest the idea of parallel computation and apply it. Three theoretical models are tested and the execution efficiency is compared in different situations. The results indicate that the parallel 3D forward modeling computation is correct and the efficiency is greatly improved. This method is suitable for large size geophysical computations.展开更多
Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured...Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured-grid-based methods can model complex underground structures with high accuracy and overcome the defects of traditional methods, such as the high computational cost for improving model accuracy and the difficulty of inverting with topography. In this paper, we used the limited-memory quasi-Newton(L-BFGS) method with an unstructured finite-element grid to perform 3D MT inversions. This method avoids explicitly calculating Hessian matrices, which greatly reduces the memory requirements. After the first iteration, the approximate inverse Hessian matrix well approximates the true one, and the Newton step(set to 1) can meet the sufficient descent condition. Only one calculation of the objective function and its gradient are needed for each iteration, which greatly improves its computational efficiency. This approach is well-suited for large-scale 3D MT inversions. We have tested our algorithm on data with and without topography, and the results matched the real models well. We can recommend performing inversions based on an unstructured finite-element method and the L-BFGS method for situations with topography and complex underground structures.展开更多
For magnetotelluric sounding (MT), many processing methods based on power spectrum have put forward lots of hypotheses, such as MT signals are Gaussian, linear and minimum-phase. If practical signals do not satisfy th...For magnetotelluric sounding (MT), many processing methods based on power spectrum have put forward lots of hypotheses, such as MT signals are Gaussian, linear and minimum-phase. If practical signals do not satisfy these requirements, the results will have a few problems as follows. Firstly, when signals are non-linear and non-Gaussian, the information of the earth contained in the MT signals cannot be sufficiently extracted; Secondly, when signals are non-Gaussian and non-minimum phase, the processed results cannot reflect the minimum phase characteristics of the signals. Hence, it is necessary for us to do further research on characteristics of MT signals (YAO, SUN, 1999; LI, CHENG, 2002; Nikias, Petropulu, 1993; ZHANG, 1996). Otherwise, we cannot judge the reliability of the processed results based on power spectrum.……展开更多
Based on the analysis of impedance tensor data, tipper data, and the conjugate gradient algorithm, we develop a three-dimensional （3D） conjugate gradient algorithm for inverting magnetotelluric full information data...Based on the analysis of impedance tensor data, tipper data, and the conjugate gradient algorithm, we develop a three-dimensional （3D） conjugate gradient algorithm for inverting magnetotelluric full information data determined from five electric and magnetic field components and discuss the method to use the full information data for quantitative interpretation of 3D inversion results. Results from the 3D inversion of synthetic data indicate that the results from inverting full information data which combine the impedance tensor and tipper data are better than results from inverting only the impedance tensor data （or tipper data） in improving resolution and reliability. The synthetic examples also demonstrate the validity and stability of this 3D inversion algorithm.展开更多
To improve magnetotelluric(MT)nonlinear inversion accuracy and stability,this work introduces the deep belief network(DBN)algorithm.Firstly,a network frame is set up for training in different 2D MT models.The network ...To improve magnetotelluric(MT)nonlinear inversion accuracy and stability,this work introduces the deep belief network(DBN)algorithm.Firstly,a network frame is set up for training in different 2D MT models.The network inputs are the apparent resistivities of known models,and the outputs are the model parameters.The optimal network structure is achieved by determining the numbers of hidden layers and network nodes.Secondly,the learning process of the DBN is implemented to obtain the optimal solution of network connection weights for known geoelectric models.Finally,the trained DBN is verified through inversion tests,in which the network inputs are the apparent resistivities of unknown models,and the outputs are the corresponding model parameters.The experiment results show that the DBN can make full use of the global searching capability of the restricted Boltzmann machine(RBM)unsupervised learning and the local optimization of the back propagation(BP)neural network supervised learning.Comparing to the traditional neural network inversion,the calculation accuracy and stability of the DBN for MT data inversion are improved significantly.And the tests on synthetic data reveal that this method can be applied to MT data inversion and achieve good results compared with the least-square regularization inversion.展开更多
Broad-band and long period magnetotelluric measurements made at 63 locations along -500 km long Chikmagalur-Kavali profile,that cut across the Dharwar craton (DC) and Eastern Ghat Mobile Belt (EGMB) in south India,is ...Broad-band and long period magnetotelluric measurements made at 63 locations along -500 km long Chikmagalur-Kavali profile,that cut across the Dharwar craton (DC) and Eastern Ghat Mobile Belt (EGMB) in south India,is modelled to examine the lithosphere architecture of the cratonic domain and define tectonic boundaries.The 2-D resistivity model shows moderately conductive features that intersperse a highly resistive background of crystalline rocks and spatially connect to the exposed schist belts or granitic intrusions in the DC.These features are therefore interpreted as images of fossil pathways of the volcanic emplacements associated with the greenstone belt and granite suite formation exposed in the region.A near vertical conductive feature in the upper mantle under the Chitradurga Shear Zone represents the Archean suture between the western and eastern blocks of DC.Although thick (-200 km) cratonic (highly resistive) lithosphere is preserved,significant part of the cratonic lithosphere below the western DC is modified due to plume-continental lithosphere interactions during the CretaceouseTertiary period.A west-verging moderately conductive feature imaged beneath EGMB lithosphere is interpreted as the remnant of the Proterozoic collision process between the Indian land mass and East Antarctica.Thin (-120 km) lithosphere is seen below the EGMB,which form the exterior margin of the India shield subsequent to its separation from East Antarctica through rifting and opening of the Indian Ocean in the Cretaceous.展开更多
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.展开更多
A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forwar...A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.展开更多
To better retain useful weak low-frequency magnetotelluric(MT)signals with strong interference during MT data processing,we propose a SVM-CEEMDWT based MT data signal-noise separation method,which extracts the weak MT...To better retain useful weak low-frequency magnetotelluric(MT)signals with strong interference during MT data processing,we propose a SVM-CEEMDWT based MT data signal-noise separation method,which extracts the weak MT signal affected by strong interference.First,the approximate entropy,fuzzy entropy,sample entropy,and Lempel-Ziv(LZ)complexity are extracted from the magnetotelluric data.Then,four robust parameters are used as the inputs to the support vector machine(SVM)to train the sample library and build a model based on the different complexity of signals.Based on this model,we can only consider time series with strong interference when using the complementary ensemble empirical mode decomposition(CEEMD)and wavelet threshold(WT)for noise suppression.Simulation results suggest that the SVM based on the robust parameters can distinguish the time periods with strong interference well before noise suppression.Compared with the CEEMD WT,the proposed SVM-CEEMDWT method retains more low-frequency low-variability information,and the apparent resistivity curve is smoother and more continuous.Moreover,the results better reflect the deep electrical structure in the field.展开更多
基金supported by the China Geological Survey(Grant No.DD20230254)。
文摘The Beishan orogen,located in the central segment of the Tianshan–Solonker suture within the southern Central Asian Orogenic Belt(CAOB),is crucial for understanding the accretionary processes and continental growth in Central Asia.This orogen developed through the episodic amalgamation and accretion of continental margin arcs,island arcs,ophiolites,and accretionary wedges,undergoing a complex process of accretion and evolution.Since the Phanerozoic,the Beishan orogen has experienced multiple phases of magmatic and collision events.The intricate distribution of magmatic arc rocks has obscured the complete basement traces,and the spatial superposition of multiple magmatic arc phases has complicated the study of its evolutionary history.
文摘In order to attain good quality transfer function estimates from magnetotelluric field data(i.e.,smooth behavior and small uncertainties across all frequencies),we compare time series data processing with and without a multitaper approach for spectral estimation.There are several common ways to increase the reliability of the Fourier spectral estimation from experimental(noisy)data;for example to subdivide the experimental time series into segments,taper these segments(using single taper),perform the Fourier transform of the individual segments,and average the resulting spectra.
文摘In the framework of a mineral system approach,a combination of components is required to develop a mineral system.This includes the whole-lithosphere architecture,which controls the transport of ore-forming fluids,and favorable tectonic and geodynamic processes,occurring at various spatial and temporal scales,that influence the genesis and evolution of ore-forming fluids(Huston et al.,2016;Groves et al.,2018;Davies et al.,2020).Knowledge of the deep structural framework can advance the understanding of the development of a mineral system and the emplacement of mineral deposits.Deep geophysical exploration carried out with this aim is increasingly important for targeting new ore deposits in unexplored and underexplored regions(Dentith et al.,2018;Dentith,2019).
基金supported by the National Natural Science Foundation of China(42372339)the China Geological Survey Project(DD20221816,DD20190319)。
文摘It is of crucial importance to investigate the spatial structures of ancient landslides in the eastern Tibetan Plateau’s alpine canyons as they could provide valuable insights into the evolutionary history of the landslides and indicate the potential for future reactivation.This study examines the Deda ancient landslide,situated in the Chalong-ranbu fault zone,where creep deformation suggests a complex underground structure.By integrating remote sensing,field surveys,Audio-frequency Magnetotellurics(AMT),and Microtremor Survey Method(MSM)techniques,along with engineering geological drilling for validation,to uncover the landslide’s spatial feature s.The research indicates that a fault is developed in the upper part of the Deda ancient landslide,and the gully divides it into Deda landslide accumulation zoneⅠand Deda landslide accumulation zoneⅡin space.The distinctive geological characteristics detectable by MSM in the shallow subsurface and by AMT in deeper layers.The findings include the identification of two sliding zones in the Deda I landslide,the shallow sliding zone(DD-I-S1)depth is approximately 20 m,and the deep sliding zone(DD-I-S2)depth is 36.2-49.9 m.The sliding zone(DD-Ⅱ-S1)depth of the DedaⅡlandslide is 37.6-43.1 m.A novel MSM-based method for sliding zone identification is proposed,achieving less than 5%discrepancy in depth determination when compared with drilling data.These results provide a valuable reference for the spatial structural analysis of large-deepseated landslides in geologically complex regions like the eastern Tibetan Plateau.
基金funded by the Anhui Province Science Program(2018-g-1-4)the National Science Program(Nos.41630320,42174087)the National Key R&D Program of China(No.2016YFC0600201)。
文摘Despite extensive efforts to understand the tectonic evolution of the Jiangnan Orogen in South China,the orogenic process and its mechanism remain a matter of dispute.Previous geodynamic studies have mostly focused on collisional orogeny,which is commonly invoked to explain the Jiangnan Orogen.However,it is difficult for such hypotheses to reconcile all the geological and geophysical data,especially the absence of ultrahigh-pressure metamorphic rocks.Based on the magnetotelluric data,we present a group of resistivity models produced through the combination of two-dimensional and three-dimensional inversions,revealing the geo-electrical structures of Jiangnan and a typical collisional orogen.In our models,the resistive crust is separated into three parts by a prominent conductive layer with opposite dipping directions on both sides.A special thrust-nappe system,which is different from that developed in a typical collisional process,is revealed in the Jiangnan Orogen.This structure suggests a process different from the simple collisional orogeny.To interpret our observations,an'intra-continental orogeny'is proposed to address the development of the Jiangnan Orogen in the Mesozoic.Furthermore,this'reworked'process may contain at least two stages caused by the decoupling of the lithosphere,which is revealed by an extra conductive layer beneath Jiangnan.
基金the National Natural Science Foundation of China(Nos.42220104002,42104073,and 41630317).
文摘With the increase in the coverage area of magnetotelluric data,three-dimensional magnetotelluric modeling in spherical coordinates and its differences with respect to traditional Cartesian modeling have gradually attracted attention.To fully understand the influence of the Earth’s curvature and map projection deformations on Cartesian modeling,qualitative and quantitative analyses based on realistic three-dimensional models need to be examined.Combined with five representative map projections,a type of model conversion method that transforms the original spherical electrical conductivity model to Cartesian coordinates is described in this study.The apparent resistivity differences between the spherical western United States electrical conductivity model and the corresponding five Cartesian models are then compared.The results show that the cylindrical equal distance map projection has the smallest error.A meridian convergence correction resulting from the deformation of the map projection is introduced to rotate the Cartesian impedance tensor from grid north to geographic north,which reduces differences from the spherical results.On the basis of the magnetotelluric field data,the applicability of the Cartesian coordinate system to western and contiguous United States models is quantitatively evaluated.Precise interpretations of the contiguous United States model were found to require spherical coordinates.
基金sponsored by National Natural Science Foundation of China(Grant No.40774029,40374024)the National Hi-tech Rsearch and Development Program of China(863 Program)(No.2007AA09Z310,)the Program for New Century Excellent Talents in University(NCET)
文摘We implement a parallel algorithm with the advantage of MPI （Message Passing Interface） to speed up the rapid relaxation inversion for 3D magnetotelluric data. We test the parallel rapid relaxation algorithm with synthetic and real data. The execution efficiency of the algorithm for several different situations is also compared. The results indicate that the parallel rapid relaxation algorithm for 3D magnetotelluric inversion is effective. This parallel algorithm implemented on a common PC promotes the practical application of 3D magnetotelluric inversion and can be suitable for the other geophysical 3D modeling and inversion.
基金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.
基金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(No.42074120).
文摘Different geophysical exploration methods have significant differences in terms of exploration depth,especially in frequency domain electromagnetic(EM)exploration.According to the definition of skin depth,this difference will increase with the effective detection frequency of the method.As a result,when performing three-dimensional inversion on single type of EM data,it is not possible to effectively distinguish the subsurface geoelectric structure at the full scale.Therefore,it is necessary to perform joint inversion on different type of EM data.In this paper we combine the magnetotelluric method(MT)with the controlled-source audio-magnetotelluric method(CSAMT)to study the frequency-domain three-dimensional(3D)joint inversions,and we use the unstructured finite-element method to do the forward modeling for them,so that the numerical simulation accuracies of different electromagnetic methods can be satisfied.By combining the two sets of data,we can obtain the sensitivity of the electrical structure at different depths,and depict the full-scale subsurface geoelectric structures.In actual mineral exploration,the 3D joint inversion is more useful for identifying subsurface veins in the shallow part and blind mines in the deep part.It can delineate the morphological distribution of ore bodies more completely and provide reliable EM interpretations to guide the mining of minerals.
基金National High Technology Research and Development Program(863 Program)(No.2006AA06Z105,2007AA06Z134)
文摘Three-dimensional forward modeling magnetotellurics （MT） problems. We present a is a challenge for geometrically complex new edge-based finite-element algorithm using an unstructured mesh for accurately and efficiently simulating 3D MT responses. The electric field curl-curl equation in the frequency domain was used to deduce the H （curl） variation weak form of the MT forward problem, the Galerkin rule was used to derive a linear finite-element equation on the linear-edge tetrahedroid space, and, finally, a BI-CGSTAB solver was used to estimate the unknown electric fields. A local mesh refinement technique in the neighbor of the measuring MT stations was used to greatly improve the accuracies of the numerical solutions. Four synthetic models validated the powerful performance of our algorithms. We believe that our method will effectively contribute to processing more complex MT studies.
基金This research is sponsored by the National Natural Science Foundation of China （No. 40374024）.
文摘The workload of the 3D magnetotelluric forward modeling algorithm is so large that the traditional serial algorithm costs an extremely large compute time. However, the 3D forward modeling algorithm can process the data in the frequency domain, which is very suitable for parallel computation. With the advantage of MPI and based on an analysis of the flow of the 3D magnetotelluric serial forward algorithm, we suggest the idea of parallel computation and apply it. Three theoretical models are tested and the execution efficiency is compared in different situations. The results indicate that the parallel 3D forward modeling computation is correct and the efficiency is greatly improved. This method is suitable for large size geophysical computations.
基金financially supported by the National Natural Science Foundation of China(No.41774125)Key Program of National Natural Science Foundation of China(No.41530320)+1 种基金the Key National Research Project of China(Nos.2016YFC0303100 and 2017YFC0601900)the Strategic Priority Research Program of Chinese Academy of Sciences Pilot Special(No.XDA 14020102)
文摘Traditional 3D Magnetotelluric(MT) forward modeling and inversions are mostly based on structured meshes that have limited accuracy when modeling undulating surfaces and arbitrary structures. By contrast, unstructured-grid-based methods can model complex underground structures with high accuracy and overcome the defects of traditional methods, such as the high computational cost for improving model accuracy and the difficulty of inverting with topography. In this paper, we used the limited-memory quasi-Newton(L-BFGS) method with an unstructured finite-element grid to perform 3D MT inversions. This method avoids explicitly calculating Hessian matrices, which greatly reduces the memory requirements. After the first iteration, the approximate inverse Hessian matrix well approximates the true one, and the Newton step(set to 1) can meet the sufficient descent condition. Only one calculation of the objective function and its gradient are needed for each iteration, which greatly improves its computational efficiency. This approach is well-suited for large-scale 3D MT inversions. We have tested our algorithm on data with and without topography, and the results matched the real models well. We can recommend performing inversions based on an unstructured finite-element method and the L-BFGS method for situations with topography and complex underground structures.
基金National Natural Science Foundation of China (40274039).
文摘For magnetotelluric sounding (MT), many processing methods based on power spectrum have put forward lots of hypotheses, such as MT signals are Gaussian, linear and minimum-phase. If practical signals do not satisfy these requirements, the results will have a few problems as follows. Firstly, when signals are non-linear and non-Gaussian, the information of the earth contained in the MT signals cannot be sufficiently extracted; Secondly, when signals are non-Gaussian and non-minimum phase, the processed results cannot reflect the minimum phase characteristics of the signals. Hence, it is necessary for us to do further research on characteristics of MT signals (YAO, SUN, 1999; LI, CHENG, 2002; Nikias, Petropulu, 1993; ZHANG, 1996). Otherwise, we cannot judge the reliability of the processed results based on power spectrum.……
基金supported by the National Hi-tech Research and Development Program of China（863Program）（No.2007AA09Z310） National Natural Science Foundation of China（Grant No.40774029 40374024）+1 种基金 the Fundamental Research Funds for the Central Universities（Grant No.2010ZY53） the Program for New Century Excellent Talents in University（NCET）
文摘Based on the analysis of impedance tensor data, tipper data, and the conjugate gradient algorithm, we develop a three-dimensional （3D） conjugate gradient algorithm for inverting magnetotelluric full information data determined from five electric and magnetic field components and discuss the method to use the full information data for quantitative interpretation of 3D inversion results. Results from the 3D inversion of synthetic data indicate that the results from inverting full information data which combine the impedance tensor and tipper data are better than results from inverting only the impedance tensor data （or tipper data） in improving resolution and reliability. The synthetic examples also demonstrate the validity and stability of this 3D inversion algorithm.
基金Project(41304090)supported by the National Natural Science Foundation of ChinaProject(2016YFC0303104)supported by the National Key Research and Development Project of ChinaProject(DY135-S1-1-07)supported by Ocean 13th Five-Year International Marine Resources Survey and Development of China
文摘To improve magnetotelluric(MT)nonlinear inversion accuracy and stability,this work introduces the deep belief network(DBN)algorithm.Firstly,a network frame is set up for training in different 2D MT models.The network inputs are the apparent resistivities of known models,and the outputs are the model parameters.The optimal network structure is achieved by determining the numbers of hidden layers and network nodes.Secondly,the learning process of the DBN is implemented to obtain the optimal solution of network connection weights for known geoelectric models.Finally,the trained DBN is verified through inversion tests,in which the network inputs are the apparent resistivities of unknown models,and the outputs are the corresponding model parameters.The experiment results show that the DBN can make full use of the global searching capability of the restricted Boltzmann machine(RBM)unsupervised learning and the local optimization of the back propagation(BP)neural network supervised learning.Comparing to the traditional neural network inversion,the calculation accuracy and stability of the DBN for MT data inversion are improved significantly.And the tests on synthetic data reveal that this method can be applied to MT data inversion and achieve good results compared with the least-square regularization inversion.
基金under the Supra Institutional Project (SIP-0012)carried under INDEX (PSC0204) project,funded by Council of Scientific and Industrial Research (CSIR), New Delhi, India
文摘Broad-band and long period magnetotelluric measurements made at 63 locations along -500 km long Chikmagalur-Kavali profile,that cut across the Dharwar craton (DC) and Eastern Ghat Mobile Belt (EGMB) in south India,is modelled to examine the lithosphere architecture of the cratonic domain and define tectonic boundaries.The 2-D resistivity model shows moderately conductive features that intersperse a highly resistive background of crystalline rocks and spatially connect to the exposed schist belts or granitic intrusions in the DC.These features are therefore interpreted as images of fossil pathways of the volcanic emplacements associated with the greenstone belt and granite suite formation exposed in the region.A near vertical conductive feature in the upper mantle under the Chitradurga Shear Zone represents the Archean suture between the western and eastern blocks of DC.Although thick (-200 km) cratonic (highly resistive) lithosphere is preserved,significant part of the cratonic lithosphere below the western DC is modified due to plume-continental lithosphere interactions during the CretaceouseTertiary period.A west-verging moderately conductive feature imaged beneath EGMB lithosphere is interpreted as the remnant of the Proterozoic collision process between the Indian land mass and East Antarctica.Thin (-120 km) lithosphere is seen below the EGMB,which form the exterior margin of the India shield subsequent to its separation from East Antarctica through rifting and opening of the Indian Ocean in the Cretaceous.
基金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.
基金Project(60672042) supported by the National Natural Science Foundation of China
文摘A finite element algorithm combined with divergence condition was presented for computing three-dimensional(3D) magnetotelluric forward modeling. The finite element equation of three-dimensional magnetotelluric forward modeling was derived from Maxwell's equations using general variation principle. The divergence condition was added forcedly to the electric field boundary value problem, which made the solution correct. The system of equation of the finite element algorithm was a large sparse, banded, symmetric, ill-conditioned, non-Hermitian complex matrix equation, which can be solved using the Bi-CGSTAB method. In order to prove correctness of the three-dimensional magnetotelluric forward algorithm, the computed results and analytic results of one-dimensional geo-electrical model were compared. In addition, the three-dimensional magnetotelluric forward algorithm is given a further evaluation by computing COMMEMI model. The forward modeling results show that the algorithm is very efficient, and it has a lot of advantages, such as the high precision, the canonical process of solving problem, meeting the internal boundary condition automatically and adapting to all kinds of distribution of multi-substances.
基金funded by the National Key R&D Program of China(No.2018YFC0603202)the National Natural Science Foundation of China(No.41404111)+1 种基金Natural Science Foundation of Hunan Province(No.2018JJ2258)Hunan Provincial Science and Technology Project Foundation(No.2018TP1018)
文摘To better retain useful weak low-frequency magnetotelluric(MT)signals with strong interference during MT data processing,we propose a SVM-CEEMDWT based MT data signal-noise separation method,which extracts the weak MT signal affected by strong interference.First,the approximate entropy,fuzzy entropy,sample entropy,and Lempel-Ziv(LZ)complexity are extracted from the magnetotelluric data.Then,four robust parameters are used as the inputs to the support vector machine(SVM)to train the sample library and build a model based on the different complexity of signals.Based on this model,we can only consider time series with strong interference when using the complementary ensemble empirical mode decomposition(CEEMD)and wavelet threshold(WT)for noise suppression.Simulation results suggest that the SVM based on the robust parameters can distinguish the time periods with strong interference well before noise suppression.Compared with the CEEMD WT,the proposed SVM-CEEMDWT method retains more low-frequency low-variability information,and the apparent resistivity curve is smoother and more continuous.Moreover,the results better reflect the deep electrical structure in the field.
