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

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.

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.

FAST APPROXIMATE INVERSION OF MAGNETOTELLURIC(FAIMT) DATA FOR A HORIZONTALLY LAYERED EARTH WITH NEAR-SURFACE INHOMOGENEITIES

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.

基于VM-MT3D的污染物质运移数值模拟预测分析

加强地下水溶质运移的数值模拟研究,对于地下水资源的保护、利用具有重要的现实意义。利用VM-MT3D数值模型,对研究区地下水NH 3-N溶质渗漏污染进行多年模拟、预测,分析渗漏处理站和处理站下游厂界区域的污染浓度历时变化过程,发现NH 3-N污染羽最远影响距离为496.2 m,污染羽范围为36800 m 2;另外,处理站和处理站下游厂界分别在7.4、14.8年后才能达到地下水Ⅲ类标准限值。最后,对地下水溶质运移模型就行验证识别,发现模型模拟的地下水位和实际水位差均在0~1.2 m之间,水位误差较小,满足模拟预测所需精度要求。

Modeling the 3D Terrain Effect on MT by the Boundary Element Method

A numerical method is put forward in this paper, using the boundary element method （BEM） to model 3D terrain effects on magnetotelluric （MT） surveys, Using vector integral theory and electromagnetic field boundary conditions, the boundary problem of two electromagnetic fields in the upper half space （air） and lower half space （earth medium） was transformed into two vector integral equations just related to the topography ： one magnetic equation for computing the magnetic field and the other electrical equation for computing the electrical field. The topography integral is decomposed into a series of integrals in a triangle element. For the integral in a triangle element, we suppose that the electromagnetic field in it is the stack of the electromagnetic field in the homogeneous earth and the topography response which is a constant; so the computation becomes simple, convenient and highly accurate. By decomposition and computation, each vector integral equation can be calculated by solving three linear equations that are related to the three Cartesian directions. The matrix of these linear equations is diagonally dominant and can be solved using the Symmetric Successive Over-Relaxation （SSOR） method. The apparent resistivity curve of MT on two 3D terrains calculated by BEM is shown in this paper.

Three-dimensional magnetotelluric forward modeling using structural-point EFGM and FEM coupling method

In this paper,the authors propose a method of three-dimensional(3D)magnetotelluric(MT)forward modeling algorithm based on the meshfree and finite element coupling method.The model is discretized by regular nodes in the central area,and the radial point interpolation method(RPIM)based on the global weakness is utilized to construct the meshfree shape function.The Governing equations in each background gird are solved by Gaussian integration.In the extended area where the points are sparsely distributed,to avoid the instability of the meshfree method,finite element method(FEM)with regular grids is used to solve the governing equation.Finally,the meshfree and finite element governing equations are coupled by the continuity of the field at the interfaces,and the direct solution technique is used to realize the 3D MT forward modeling.Numerical experiments of several typical electrical models are used to verify the effectiveness of the method.

Thre-dimensional topographic responses in MT using finite diference method

In this paper, the magnetotelluric response to 3-D earth topography was studied using the finite difference method for rectangular cells. Three numerical modeling calculations were carried out; the results indicate that the finite difference method is applicable to the magnetotelluric problem in estimating the 3-D topographic responses.

关于MT二维反演中数据旋转方向的选择问题初探

文中通过二维及相应的三维理论模型响应在不同方向上反演结果的对比分析,探讨了大地电磁一、二维反演中数据旋转方向的选择问题。发现在纯二维情况下,宜采用主轴方向的数据进行反演。此时如采用旋转到测线方向的数据,则要求测线与构造倾向之间的夹角不能太大(对于本文的模型不能大于30°)。而在三维情况下,采用旋转到最佳主轴方向的数据进行二维反演未必能得到比旋转到测线方向更优的反演结果。由于实测资料一般都可视为三维响应,故本文认为在实际资料处理中,对于那些构造走向明显,而测线是尽量垂直于构造的情况下,可直接将数据旋转到测线方向后进行二维反演,从而避免最佳主轴方向上极化模式选取的困难。

庐枞矿集区矾山—将军庙地区AMT三维反演及地质结构解释

庐江—枞阳矿集区是长江中下游成矿带内的一个重要的多矿种成矿区.为揭示矾山—将军庙地区的地质结构以指导找矿,首先,在庐枞矿集区北部的矾山—将军庙地区开展了三维音频大地电磁数据采集工作,采用Rhoplus方法进行AMT数据死频带校正保证数据质量,利用相位张量分析地下介质的维性特征;其次,利用模块化三维反演并行代码ModEM开展了实测数据的三维反演,获得了一个典型的具有双层结构特征的火山岩盆地三维地下电性模型;最后通过电性模型及其他地质资料的综合解释,依据反演模型中近地表的低阻层和深部的高阻隆起区,分别勾画出了火山岩地层的厚度和深部侵入岩的分布范围,并在小岭地区发现地表浅部下方存在一个含矿次级火山机构的高导体.

Three-dimensional electrical structure model of the Yangbajain geothermal field in Tibet:Evidence obtained from magnetotelluric data

The Yangbajain Geothermal Field in Tibet is located in the fault subsidence basin of the central Yadong-Gulu Rift Valley.The spatial distribution of the field is controlled by mountain-front fault zones on the northwestern and southeastern sides of the basin.Geothermal power has been generated in Yangbajain for more than 40 years.However,owing to the lack of threedimensional(3D) geophysical exploration data,key geological issues related to the partial melt body of the Yangbajain Geothermal Field,such as its location,burial depth,and geometric form,as well as the ascending channel of the geothermal fluid,have for a long time been controversial.In this study,3D inversion was performed using measured geo-electromagnetic total impedance tensor data from 47 survey points.The extracted horizontal sections at different depths and profiles,and at different lines,reflect the 3D electrical structure model of the geothermal field in the study area.Subsequently,three findings were obtained.First,the partial melt body,located below the China-Nepal Highway extending along the northeast direction,is the heat source of the Yangbajain Geothermal Field.The burial depth range of the molten body was determined to range between approximately 6.2 and 14 km.Moreover,the geothermal fluid ascended a horn-shaped circulation channel with an up-facing opening,located in the northern section of the sulfur ditch area.The study results revealed that deep rock fissures(>2 km) were not well developed and had poor permeability.In addition,no layered heat reservoirs with high water richness were observed in the northern part of the study area.However,the application of enhanced geothermal system(EGS) technology in the northern region would be essential to improving the power generation capacity of the Yangbajain Geothermal Field.In addition,the study found no deep high-temperature heat storage areas in the southern region of the study area.

Magnetotelluric data reveals subduction polarity and reactivation of the Mudanjiang suture zone,Northeast China

The Jiamusi and Songnen blocks converged in the easternmost segment of the Central Asian Orogenic Belt as a result of the subduction and subsequent closure of the Mudanjiang oceanic plate during the Permian-Jurassic.The Mudanjiang suture zone was later directly affected by subductions of the Paleo-Pacific plate and Pacific plate and is therefore an ideal place to study the subduction polarity and later transformation of a paleo-suture zone.Using three-dimensional inversion of magnetotelluric data collected along a 160-km-long profile across the Mudanjiang suture zone,we established a resistivity model of the suture zone and adjacent area.Our results reveal the subduction polarity and subduction trace of the Mudanjiang oceanic plate and provide geoelectrical evidence for reactivation of the Mudanjiang suture zone induced by the(Paleo-)Pacific plate subduction.The suture zone shows a complex conductive structure.The west-dipping crustal-scale conductor beneath the Songnen-Jiamusi collision zone represents the fossil subduction zone and indicates the westward subduction polarity of the Mudanjiang oceanic plate.Furthermore,the Mudanjiang fault identified by surface geology does not fully represent the deep structure of the Mudanjiang suture zone.The definition of the suture zone should be extended to the whole conductive region with a lateral extent of~70 km.Solid conductive minerals beneath the arc in front of the subduction zone were exhumated up from deep to the upper crust.The“chimney”-shaped conductor connected with the mantle represents the intrusive pathways of mantle-derived materials,suggesting that the Mudanjiang suture zone was reactivated by subductions of the Paleo-Pacific plate and Pacific plate,leading to remelting of the cooled and crystallized materials in the pathways.Therefore,subduction of the(Paleo-)Pacific plate destroyed the lithospheric structure of the paleo collision zone in the eastern segment of the Central Asian orogenic belt,and the large-scale crustal conductor beneath the suture zone reflects reactivation of the paleo-suture zone.

NEC MT860＋投影机直面普教用户

NEC公司今年新推出的MT860+投影机采用LCD显示技术,整机界面友好,SVGA分辨率高达3000ANSI流明,仅重5.9公斤.这款投影机具备多项先进功能,例如支持有线/无线LAN、3D校正(实现斜角投影)以及背景色校正等新功能元素.

NEC MT系列投影机新一代多媒体教学的得力助手

NEC公司的MT系列投影机在教育行业一直有着广泛的应用和良好的口碑。目前,新产品MT1065+/1060+/860+已经问世。MT1065+/1060+/860+不但高亮度(3500/3000/3000ANSI流明)、高对比度800∶1、最高可兼容至UXGA(1600×1200)分辨率,还具备以下功能:

基于水量水质耦合模拟优化的渠井结合灌区多目标水资源优化配置模型与方法

为探究灌区考虑水质水量耦合模拟优化的水资源优化配置模型与方法,基于FloPy和Pymoo构建了水量水质耦合的多目标模拟优化模型。采用地下水数值模拟模型MODFLOW和溶质输移模型MT3D模拟了灌溉引起的地下水位变化和污染物的扩散过程,考虑灌溉收益、作物水分生产率、地下水位变幅和污染扩散4个目标,分别采用多目标遗传算法NSGA-Ⅱ和NSGA-Ⅲ对该问题进行求解,并对比分析了两种方法的实际应用效果。采用多准则分析评价方法TOPSIS(Technique for Order Preference by Similarity to Ideal Solution),推荐了灌溉方案。推荐灌溉方案在平水年和枯水年分别可以最大减少地下水位变幅56%和43%,最大减少污染扩散面积54%和40%。

Three-Dimensional Conjugate Gradient Inversion of Magnetotelluric Impedance Tensor Data

We developed a three-dimensional（3D） conjugate gradient inversion algorithm for in-verting magnetotelluric impedance tensor measurements.In order to show the importance of including diagonal components of magnetotelluric impedance tensor in 3D inversion,synthetic data were inverted using the 3D conjugate gradient inversion,and the inversion results were compared and analyzed.The results from the 3D inversion of synthetic data indicate that both the off-diagonal and the diagonal components are required in inversions to obtain better inversion results when there are no enough data sites to recover the target resistivity structure.These examples show that lots of information about 3D structure is also contained in the diagonal components;as a result,diagonal components should be in-cluded in 3D inversions.The inversion algorithm was also used to invert the impedance tensor data ac-quired in the Kayabe area in Japan.Inversions with the synthetic and real data demonstrated the va-lidity and practicability of the inversion algorithm.

Numerical modeling of 3-D terrain effect on MT field

Using the boundary element method, the numerical modeling problem of three-dimensional terrain effect on magnetotelluric (MT) field is solved. This modeling technique can be run on PC in the case of adopting special net division. The result of modeling test for 2-D terrain by this modeling technique is basically coincident with that by 2-D modeling technique, but there is a great difference between the results of 3-D and 2-D modeling for 3-D terrain.

Modeling of 3-D terrain effect on MT field using BEM

THE modeling of 2-D terrain effect on magnetotelluric(MT)has been solved using finite ele-ment method(FEM)and boundary element method(BEM).This note uses BEM tosolve the problem.The solution of H_x- and E_x-polarization wave can be obtained by the samemethod.Thus,this note only discusses the solution of H_x-polarization wave.

球坐标系下三维大地电磁正演研究

大地电磁正演理论研究热点一直以来主要集中在如何提高计算效率和精度,但在剖面足够长、探测深度足够大的情况下,传统的笛卡尔坐标系数值模拟方式难以准确拟合地球曲率形态.本文研究了基于球坐标系的三维大地电磁正演,推导了交错网格有限差分三维正演公式,与一维解析解和三维标准模型测试对比,验证了正演算法的正确性.通过理论模型计算,对比分析球坐标和笛卡尔坐标系正演结果表明:球坐标系模拟更合理,避免了传统笛卡尔坐标拉伸投影所引入的误差,可代替目前的笛卡尔坐标模拟方法.基于球坐标和笛卡尔坐标系的三维大地电磁正演响应值随着频率变低差异越明显.球坐标和笛卡尔坐标计算结果差异度与频率、模型结构和电阻率有关.本文模型计算结果在数万秒周期处已出现接近10%的差异,对于较大尺度的长周期大地电磁,地球曲率的影响不能忽略.