To extend our successful Magnetotelluric(MT) experiments in the area from Yadong to Bamocuo in 1995 (Chen et al., 1996), two new Magnetotelluric(MT)experiments have been carried out along two main profiles in the cent...To extend our successful Magnetotelluric(MT) experiments in the area from Yadong to Bamocuo in 1995 (Chen et al., 1996), two new Magnetotelluric(MT)experiments have been carried out along two main profiles in the central and Northern Tibetan Plateau in the summer of 1998 and 1999. While the 1995 MT work mainly focused on the study of the electrical structure of the Yarlung\|Zangbo River Suture , the 1998 and 1999 experiments have been designed with following purpose:(1) Study whether partially molten layer widely exists in the crust of the central and northern Tibet.(2) Study the electrical structures of crust and upper mantle in central and northern Tibet which may relate to the attenuated Shear Seismic waves.(3) Study the detail electrical structures of Bangong\|Nujiang Suture and Jinsha Suture.In 1998, the MT team (China University of Geosciences, University of Washington and Geological Survey of Canada) recorded MT data along the Tibet 500 Line which extends about three hundred and eighty kilometers from Deqing to Longweicuo. We used LIMS system to record the long period(20~20000s) MT data at twenty\|six sites and used MT24 to record broadband(320Hz,2000s) MT data at fifty\|eight sites.展开更多
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.展开更多
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 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.展开更多
This paper discusses use of approximations and the Integral Mean Value Theorem to show that 6 coefficients approximately describe the distortions of near surface inhomogeneities on the MT field of a horizontally layer...This paper discusses use of approximations and the Integral Mean Value Theorem to show that 6 coefficients approximately describe the distortions of near surface inhomogeneities on the MT field of a horizontally layered earth model. When these 6 coefficients are considered together with those of the magnetic field of a horizontally layered earth model,the analytic and approximate wave impedance equations can be derived for the MT response of a horizontally layered earth model with near-surface 2-D and 3-D inhomogeneities. These approximate wave impedance equations are used with inverted MT data for 2-D and 3-D forward modelling. Although these 6 coefficients cannot be determined before inversion,initial estimates can be used. The 6 coefficients and the asistivity and thickness of each layer of a horizontally layered earth can be obtained by using published inversion methods. The 6 coefficients give important informaion (depths and resistivities) on the near-surface inhomogenelties.The authors inverted 2-D and 3-D theoretical models for Fast Approximate Inversion of Magnetotelluric (FAIMT) data for a horizontally layered earth with near-surface inhomogeneities compares favorably with traditional invrsion methods, especially for inverting regional or basin structures. This method simplifies computation and gives a reasonable 1 -D geological model with fewer nonuniquenas problems.展开更多
The INDEPTH MT results show that there are no electrical features of deep fractures along the Yarlung Zangbo River, but a large high conductivity body exists in the area between Gyangze and Rinbung. It dips northwar...The INDEPTH MT results show that there are no electrical features of deep fractures along the Yarlung Zangbo River, but a large high conductivity body exists in the area between Gyangze and Rinbung. It dips northward, extends downward up to the depth of about 55 km and indicates the exposure of the possible real position of the Yarlung Zangbo suture. There are three sets of electrical gradient and distorted zones reflecting the structure of faults in the high conductivity region. These three fault belts, which dip northward and gradually converge downward to the main fault belt, and a series of south dipping faults in the north side form the exhaled structural characteristics of the Yarlung Zangbo suture. There is a close relationship between the large high conductivity body and the underground thermal state.展开更多
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.展开更多
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.展开更多
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.……展开更多
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.展开更多
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.展开更多
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.展开更多
In recent years, the data collecting and processing of magnetotelluric sounding (MT) have progressed greatly. How to improve the quality of field surveyed data and obtain the high quality parameters such as apparent...In recent years, the data collecting and processing of magnetotelluric sounding (MT) have progressed greatly. How to improve the quality of field surveyed data and obtain the high quality parameters such as apparent resistivity and phase is the most important link during the overall flow. It is closely related to the signal-to-noise ratio of magnetic sensors, the resolution of data collecting units and relevant processing methods as well. Produced by Germany Metronix Measurement Instruments and Electronics Ltd., the new type multi-channel geophysical measurement system GMS-06 was already adopted and used by domeslic units such as Institute of Geology and Geophysics of Chinese Academy of Sciences and Tongji University. Based on the brief introduction of GMS-06, this paper will give the comparisons and analyses of the synchronous surveyed MT data by GMS-06, MMS-03 MT system which was also produced by Metronix and V5-2000 MT system produced by Canada Phoenix Geophysics Ltd.. Then the preliminary processing and interpretation results of the long period MT data observed by GMS-06 in the area of Shanghai and Zhejiang are also discussed in the paper.展开更多
文摘To extend our successful Magnetotelluric(MT) experiments in the area from Yadong to Bamocuo in 1995 (Chen et al., 1996), two new Magnetotelluric(MT)experiments have been carried out along two main profiles in the central and Northern Tibetan Plateau in the summer of 1998 and 1999. While the 1995 MT work mainly focused on the study of the electrical structure of the Yarlung\|Zangbo River Suture , the 1998 and 1999 experiments have been designed with following purpose:(1) Study whether partially molten layer widely exists in the crust of the central and northern Tibet.(2) Study the electrical structures of crust and upper mantle in central and northern Tibet which may relate to the attenuated Shear Seismic waves.(3) Study the detail electrical structures of Bangong\|Nujiang Suture and Jinsha Suture.In 1998, the MT team (China University of Geosciences, University of Washington and Geological Survey of Canada) recorded MT data along the Tibet 500 Line which extends about three hundred and eighty kilometers from Deqing to Longweicuo. We used LIMS system to record the long period(20~20000s) MT data at twenty\|six sites and used MT24 to record broadband(320Hz,2000s) MT data at fifty\|eight sites.
文摘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.
基金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 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.
文摘This paper discusses use of approximations and the Integral Mean Value Theorem to show that 6 coefficients approximately describe the distortions of near surface inhomogeneities on the MT field of a horizontally layered earth model. When these 6 coefficients are considered together with those of the magnetic field of a horizontally layered earth model,the analytic and approximate wave impedance equations can be derived for the MT response of a horizontally layered earth model with near-surface 2-D and 3-D inhomogeneities. These approximate wave impedance equations are used with inverted MT data for 2-D and 3-D forward modelling. Although these 6 coefficients cannot be determined before inversion,initial estimates can be used. The 6 coefficients and the asistivity and thickness of each layer of a horizontally layered earth can be obtained by using published inversion methods. The 6 coefficients give important informaion (depths and resistivities) on the near-surface inhomogenelties.The authors inverted 2-D and 3-D theoretical models for Fast Approximate Inversion of Magnetotelluric (FAIMT) data for a horizontally layered earth with near-surface inhomogeneities compares favorably with traditional invrsion methods, especially for inverting regional or basin structures. This method simplifies computation and gives a reasonable 1 -D geological model with fewer nonuniquenas problems.
基金The study is supported by National Natural Science Foundation of China
文摘The INDEPTH MT results show that there are no electrical features of deep fractures along the Yarlung Zangbo River, but a large high conductivity body exists in the area between Gyangze and Rinbung. It dips northward, extends downward up to the depth of about 55 km and indicates the exposure of the possible real position of the Yarlung Zangbo suture. There are three sets of electrical gradient and distorted zones reflecting the structure of faults in the high conductivity region. These three fault belts, which dip northward and gradually converge downward to the main fault belt, and a series of south dipping faults in the north side form the exhaled structural characteristics of the Yarlung Zangbo suture. There is a close relationship between the large high conductivity body and the underground thermal state.
基金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.
基金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.
基金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.……
基金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.
基金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.
基金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.
基金National Natural Science Foundation of China (40104005) and National High Technical Research and Development Project (2004AA615010).
文摘In recent years, the data collecting and processing of magnetotelluric sounding (MT) have progressed greatly. How to improve the quality of field surveyed data and obtain the high quality parameters such as apparent resistivity and phase is the most important link during the overall flow. It is closely related to the signal-to-noise ratio of magnetic sensors, the resolution of data collecting units and relevant processing methods as well. Produced by Germany Metronix Measurement Instruments and Electronics Ltd., the new type multi-channel geophysical measurement system GMS-06 was already adopted and used by domeslic units such as Institute of Geology and Geophysics of Chinese Academy of Sciences and Tongji University. Based on the brief introduction of GMS-06, this paper will give the comparisons and analyses of the synchronous surveyed MT data by GMS-06, MMS-03 MT system which was also produced by Metronix and V5-2000 MT system produced by Canada Phoenix Geophysics Ltd.. Then the preliminary processing and interpretation results of the long period MT data observed by GMS-06 in the area of Shanghai and Zhejiang are also discussed in the paper.