Spatial structural characteristics of the Deda ancient landslide in the eastern Tibetan Plateau:Insights from Audio-frequency Magnetotellurics and the Microtremor Survey Method

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.

Imaging the Architecture of Mineral Systems and the Pathways of Ore-forming Fluids across Mongolia with Magnetotellurics

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).

3D joint inversion of controlled-source audio-frequency magnetotelluric and magnetotelluric data

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.

Improved Responses with Multitaper Spectral Analysis for Magnetotelluric Time Series Data Processing:Examples from Field Data

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.

Tectonic Evolution and Lithospheric Structure of the Beishan Orogen:Insights from Magnetotelluric Studies

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.

Three-dimensional magnetotellurics modeling using edgebased finite-element unstructured meshes

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.

Proterozoic tectonothermal processes imaged with magnetotellurics and seismic reflection in southern Australia

Over the last two decades,co-located seismic and magnetotelluric(MT) profiles provided fundamental geophysical data sets to image the Australian crust.Despite their complimentary nature,the data are processed and often interpreted separately without common processes in mind.We here qualitatively compare 2 D resistivity inversion models derived from MT and seismic reflection profiles across a region of Archean-Proterozoic Australia to address the causes of variations in seismic response and anomalous conductivity in the crust.We find that there exists a spatial association between regions of low reflectivity in seismic sections and low resistivity in co-located2 D MT modelled sections.These relationships elucidate possible signatures of past magmatic and fluid-related events.Depending on their diffuse or discrete character,we hypothesize these signatures signify fossil melting of the crust due to mafic underplating,magma movement or hydrothermal fluid flow through the crust.The approach discussed herein is a process-oriented approach to interpretation of geophysical images and a significant extension to traditional geophysical methods which are primarily sensitive to a singular bulk rock property or state.

3D Step-by-step inversion strategy for audio magnetotellurics data based on unstructured mesh

A three-dimensional(3D)step-by-step inversion strategy for audio magnetotellurics(AMT)is investigated in this study.The objective function is minimized by iteratively solving the Gauss-Newton normal equation,and the inversion region is discretized with unstructured tetrahedral elements.The inversion proceeds step-by-step from a coarse mesh to a fine mesh.In the inversion iteration process,a mesh is adaptively optimized according to the spatial gradient information about the model resistivity to fine delineate the boundaries of abnormal bodies.In the early stage of inversion execution,a coarse mesh is used for inversion,and the inversion stability is improved by reducing the number of inversion elements.In addition,mesh refinement is performed in the iterative inversion process.The inversion results obtained from the previous mesh are used as the reference and initial models for the next mesh iterative inversion.The step-by-step inversion strategy can ensure that the inversion is performed in the correct direction,improving the inversion stability and results gradually.Synthetic results show that the step-by-step inversion strategy with a Gauss-Newton method for 3D AMT inversion is stable and reliable,which lays a foundation for further practical 3D AMT data inversion.

Intra-continental Orogeny:Insights from Magnetotelluric Data into the Mesozoic Uplift History of the Eastern Jiangnan Orogen in South China

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.

3D model of Svecofennian Accretionary Orogen and Karelia Craton based on geology,reflection seismics,magnetotellurics and density modelling:Geodynamic speculations

A 3D model of deep crustal structure of the Archaean Karelia Craton and late Palaeoproterozoic Svecofennian Accretionary Orogen including the boundary zone is presented.The model is based on the combination of data from geological mapping and reflection seismic studies,along profiles 1-EU,4 B,FIRE-1-2 a-2 and FIRE-3-3 a,and uses results of magnetotelluric soundings in southern Finland and northern Karelia.A seismogeological model of the crust and crust-mantle boundary is compared with a model of subhorizontal velocity-density layering of the crust.The TTG-type crust of the Palaeoarchaean and Mesoarchaean microcontinents within the Karelia Craton and the Belomorian Province are separated by gently dipping greenstone belts,at least some of which are palaeosutures.The structure of the crust was determined mainly by Palaeoproterozoic tectonism in the intracontinental settings modified by a strong collisional compression at the end of the Palaeoproterozoic.New insights into structure,origin and evolution of the Svecofennian Orogen are provided.The accretionary complex is characterized by inclined tectonic layering:the tectonic sheets,~15 km thick,are composed of volcanic-sedimentary rocks,including electro-conductive graphite-bearing sedimentary rocks,and electro-resistive granitoids,which plunge monotonously and consecutively eastward.Upon reaching the level of the lower crust,the tectonic sheets of the accretionary complex lose their distinct outlines.In the seismic reflection pattern they are replaced by a uniform acoustically translucent medium,where separate sheets can only be traced fragmentarily.The crust-mantle boundary bears a diffuse character:the transition from crust to mantle is recorded by the disappearance of the vaguely drawn boundaries of the tectonic sheets and in the gradual transition of acoustically homogeneous and translucent lower crust into transparent mantle.Under the effect of endogenic heat flow,the accretionary complex underwent high-temperature metamorphism and partial melting.Blurring of the rock contacts,which in the initial state created contrasts of acoustic impedance,was caused by partial melting and mixing of melts.The 3 D model is used as a starting point for the evolutionary model of the Svecofennian Accretionary Orogen and for determination of its place in the history of the Palaeoproterozoic Lauro-Russian intracontinental orogeny,which encompassed a predominant part of the territory of Lauroscandia,a palaeocontinent combining North American and East European cratons.The model includes three stages in the evolution of the Lauro-Russian Orogen(-2.5,2.2-2.1 and 1.95-1.87 Ga).The main feature of the Palaeoproterozoic evolution of the accretionary Svecofennian Orogen and Lauroscandia as a whole lay in the causal link with evolution of a superplume,which initiated plate-tectonic events.The Svecofennian-Pre-Labradorian palaeo-ocean originated in the superplume axial zone;the accretionary orogens were formed along both continental margins due to closure of the palaeo-ocean.

Three-dimensional magnetotelluric modeling in the spherical and Cartesian coordinate systems:A comparative study

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.

基于地球物理数据融合的大兴安岭中段东缘地壳结构研究

位于中亚造山带东段的兴安地块和松辽锡林浩特地块被晚古生代贺根山-黑河缝合带分开。随后,大兴安岭与松辽盆地被中新生代嫩江-八里罕断裂分开。在大兴安岭中段,这两构造带的位置及板块拼合时间尚不明确。为研究大兴安岭中段东缘重要构造断裂、构造单元及结晶基底属性特征,在扎兰屯-阿荣旗地区采集处理一条长60 km的重力、磁法和大地电磁数据,重新处理了1:20万扎兰屯幅、阿荣旗幅区域布格重力异常数据。通过多种地球物理数据融合,解译了7条重要断裂,其中庞家街断裂和红星断裂是控制贺根山-黑河缝合带的深大断裂。北西向阿伦河断裂为深大断裂,至少切割到中地壳。发现了两个中地壳低阻带,其中C2异常解释为贺根山-黑河缝合带。通过1:25万阿荣旗幅地质填图资料,证实了地球物理数据融合的推断结果,为大兴安岭中段东缘构造带演化研究提供了重要支撑。

地磁感应电流多学科交叉研究进展

地磁感应电流(geomagnetically induced currents,GIC)是指地球变化磁场在地球表面导体中产生的感应电流,可对电网、油气管道、高铁等人造基础设施造成严重破坏,包括电力系统失效、通信干扰和电子设备损坏等,影响其稳定运行和功能。研究GIC的意义在于深入理解其产生机理及影响因素,为制定相应的预防和应对策略提供科学依据。论文在深入总结GIC研究进展的基础上,详细阐述了GIC的产生机理及影响因素,包括宇宙空间、地球环境和电力系统对GIC的影响;概述了GIC的研究现状,分别从空间物理、地球物理和电气工程领域对GIC进行了综合评述;进而详细介绍了GIC估算的步骤,包括地磁场数据与建模、地电场和导体中GIC的计算;并重点分析了不同地球物理参数对GIC估算结果的影响以及互相之间的差异;最后,对GIC研究面临的挑战、地球物理在相关问题解决中的潜在贡献以及研究未来进行展望。

准噶尔盆地乌伦古凹陷电性结构特征

弄清准噶尔盆地乌伦古凹陷地下深部电性结构,对于准噶尔盆地乌伦古凹陷的深部油气勘探具有重要意义。综合钻井、地震、测井、地质露头观测及分析化验等资料,总结了研究区的地球物理特征,包括地层岩石密度特征、磁性特征、电性特征。在此基础上,对乌伦古凹陷的一条大地电磁测深原始时间序列的电场数据进行视电阻率曲线平滑处理和静态偏移校正处理,并采用Occam法进行反演,得到乌伦古凹陷地下10km深度范围内的电阻率剖面。依据该电阻率剖面,结合地震及钻井资料,初步确定了准噶尔盆地乌伦古凹陷地层的构造形态,发现乌伦古凹陷西北端和东南端地下存在高阻异常,且基底中存在低阻地层。研究结论对准噶尔盆地乌伦古凹陷的油气勘探具有重要的价值。

综合物探在矿山深部资源勘查中的应用

矿产资源是社会和经济发展必不可少的物质基础,目前国内多数老矿山面临已探明储量枯竭的严峻形势,现有的勘查知识和技术难以解决矿山找矿勘查所遇到的特殊难题。针对矿山找矿的特殊性,本文结合安徽省霍邱某矿区的地质及物探资料,在采矿区选取典型剖面进行重磁、微动及可控源音频大地电磁法(Controllable Source Audio Magnetotelluric Method,CSAMT)。以微动和音频大地电磁数据分析推测浅部矿体采空及塌陷区;在深部将微动勘探与重力数据结合,取长补短,通过剪切模量将微动视横波速度与地下介质密度直接联系起来,经已知矿体位置剪切模量的分布规律并结合重磁反演在300 m以深推测了两条近平行层状分布的矿体,矿体长200~250 m,最厚处约50 m;在标高100~200 m的类似S形褶皱拐点处可能为成矿有利位置;同时推测已知主矿体在底板以深还有延伸。通过本次工作,证实了在当前老矿区内,将新技术微动勘探与传统物探方法相结合,各取其长,能提高矿山深边部及外围隐伏矿床(体)定位预测的有效性,挖掘新的找矿潜力。

陆域大地电磁基础网规划与展望

阵列式大地电磁测深是区域地球物理调查的重点发展方向之一,对标国外先进水平和国家社会的基础性需求,中国阵列式大地电磁测深工作程度依然较低,缺乏整体规划,开展全国性观测计划亟需提上日程。系统总结了当前区域地球物理学的发展趋势及国内外大地电磁网建设现状,提出了合作共建陆域多层级大地电磁基础网和大数据中心的顶层设计和规划布局,并对中国大地电磁基础网未来的重点发展方向进行了展望和思考。相关工作的顺利开展将为中国的区域地球物理调查拓展一项新属性参数,推动构建物性地球、结构地球和地球物理场四维模型,进而变革区域地球物理研究范式和工作模式,支撑中国深地探测走向世界前列。

基于深度学习的大地电磁二维反演研究

为了解决传统卷积神经网络反演由于层数过多而影响准确度的问题,开展了基于残差神经网络的大地电磁二维反演研究。通过大地电磁二维正演建立大量数据集,以TE和TM模式下的视电阻率和相位数据作为四通道网络输入,以对应的地电模型作为标签及输出进行有监督的学习,利用残差神经网络实现二维大地电磁反演。不同噪声水平的地电模型反演结果表明残差网络不仅可以很好地消除层数过多带来的准确度下降问题,还具有很强的抗噪作用。对冀中坳陷实测电磁资料进行反演,获得深部碳酸盐岩电阻率分布,据此分析了工区热储构造特征。理论模型和实测数据反演结果表明该方法具有良好的学习能力和抗噪性能,反演效果稳定可靠。

基于GMT脚本的AMT数据曲线图批量绘制方法与应用研究

AMT法观测天然交变电磁场,通过数据处理和反演可获得地下介质的电性分布特征,近年来得到了广泛应用。AMT数据曲线图包括测点信息、视电阻率和阻抗相位曲线及数据,可直观反映测点数据质量和曲线类型,是重要的原始资料。加拿大Phoenix Geophysics公司的MTU-5A和V8大地电磁仪是国内最常见的AMT仪器,原始数据经预处理,用MT-Editor软件编辑后可以批量导出PLT格式文件,使用awk命令可从PLT文件中提取测点信息和视电阻率、阻抗相位数据,利用GMT绘图软件将测点信息和数据绘制成图,结合循环语句可达到批量成图的目的。AMT数据曲线图可应用于原始数据质量评价和AMT测深曲线类型判别,在工作区地电规律研究中具有重要意义。

综合物探方法在碎屑岩地区找水中的应用——以云南宣威大路边村为例

为探索地球物理方法在碎屑地区找水应用的有效性,文章以云南宣威大路边村为研究对象,选取高密度电阻率法、联合剖面法、音频大地电磁法对碎屑岩层进行找水研究。结果表明:三叠系上统飞仙关组(T1f)碎屑岩富水性中等,含水性较均匀;受地面场地狭窄限制,高密度电法及联合剖面法测线长度有限,探测深度较浅,受碎屑岩电阻率小的影响,其找水应用效果不佳,在可探测的深度范围内较难识别基岩裂隙等找水有效异常,需结合其他物探方法甄别异常;音频大地电磁法找水应用效果相对较好,具有较高的横向分辨率,可有效识别碎屑岩中的断裂破碎带及构造裂隙带,表现为相对低阻异常特征;视电阻率单支测深曲线能有效分辨不同深度的岩土层结构,具有较高的纵向分辨率,对裂隙破碎带有较好指示意义。在地面场地较为狭窄的碎屑岩地区开展找水工作,音频大地电磁法结合视电阻单支测深曲线能取得较好效果,以音频大地电磁法确定碎屑岩地区中的裂隙发育带位置,以视电阻单支测深曲线确定裂隙发育带的深度。

非均匀网格有限差分法大地电磁测深法正演模拟

为解决有限差分法大地电磁正演模拟中模拟精度与计算效率之间的矛盾,采用非均匀网格有限差分进行了大地电磁响应计算。基于二维大地电磁方程,推导出非均匀网格剖分下的TM极化模式的大地电磁方程的差分格式。通过均匀半空间模型和低阻异常体模型的正演模拟,验证了非均匀网格差分法的正确性,非均匀网格差分法可以准确反映出复杂模型的大地电磁响应特征,识别异常体分布范围。