MT数据处理中静校正方法对比

在大地电磁测深数据处理中,静态效应会影响解释结果.这里主要采用有限单元法进行正演计算,采用阻抗相位校正法、空间滤波法以及视电阻率导数解释法,对存在静态效应的理论模型进行数据处理.对比结果表明,前二种方法对静态效应起到了一定的压制和改善作用,且阻抗相位校正法效果相对较好,而通过视电阻率导数解释法,可以准确地反映地下真实异常形态.     

波西米亚块体西部长周期电磁研究(英文)

在波西米亚块体西南部做了一条长周期电磁剖面 ,它靠近德国超深钻的位置。用一个二维合成模型对数据做了大体的拟合 ,模型走向东西 ,近地表盖层高度地不均匀畸变。畸变层显示强烈的各向异性 ,优势电导方向为NW -SE至NNW -SSE。讨论了深部区域结构的几个不同的模型 ,其中一个模型可以解释由MT和GDS数据显示的主方向的偏差 ,还讨论了沿特别剖面得到的结果与区域地电结构的关系     

The electrical conductivity of Posidonia black shales — from magnetotelluric exploration to rock samples

We carried out a magnetotelluric field campaign in the South–East Lower Saxony Basin, Germany, with the main goal of testing this method for imaging regional Posidonia black shale sediments. Two‐dimensional inversion results of the magnetotelluric data show a series of conductive structures correlating with brine‐saturated sediments but also with deeper, anthracitic Westphalian/Namurian coals. None of these structures can be directly related with the Posidonia black shale, which appears to be generally resistive and therefore difficult to resolve with the magnetotelluric method. This assumption is supported by measurements of electrical resistivity on a set of Posidonia shale samples from the Hils syncline in the Lower Saxony basin. These rock samples were collected in shallow boreholes and show immature (0.53% Ro), oil (0.88% Ro), and gas (1.45% Ro) window thermal maturities. None of the black shale samples showed low electrical resistivity, particularly those with oil window maturity show resistivity exceeding 10⁴ Ωm. Moreover, we could not observe a direct correlation between maturity and electrical resistivity; the Harderode samples showed the highest resistivity, whereas the Haddessen samples showed the lowest. A similar trend has been seen for coals in different states of thermal maturation. Saturation of the samples with distilled and saline water solutions led to decreasing electrical resistivity. Moreover, a positive correlation of electrical resistivity with porosity is observed for the Wickensen and Harderode samples, which suggests that the electrical resistivity of the Posidonia black shale is mainly controlled by porosity.     

TWO DIMENSIONAL FINITE ELEMENT-BOUNDARY ELEMENT COUPLING MODEL OF MAGNETOTELLURIC SOUNDING

This paper describes how to apply the boundary element method to solve 2-D magneto-telluric sounding problem and then couple the BEM with finite element method to deal with inhomogeneous geological structure and its relation to regional stress regime. The BEM is able to calculate precisely the derivative of field variables. Such a new approach can more flexibly than others set up the upper boundary condition of the air region. Two practical examples are given. One is the study of topographical effect on the MT fields using boundary element method. There is an inverse relationship between the change of apparent resistivity and the topographical variation. The other is a coupled model of a basin in Nei Monngol Autonomous Region. A lower resistivity zone must exist in the crust and should be thicker below the central portion of the basin. This seems to agree with the regional stress regime changed from extension in Mesozoic to compression in Cenozoic     

祁连山中段深部电性结构及潜在震源危险区的研究

位于青藏高原北缘的祁连山中段是地震十分活跃的地区之一。本课题采用ＭＴ重复测量方法，在该区建立的长约４００ｋｍ的剖面上进行监测。本文根据监测所获取的ＭＴ参数，从静态的角度探讨了该区深部电性结构及其电性横向变化特征，进而评估了该区的潜在震源区。研究发现，该区深部介质电性的物质特征有：１．部分地带呈现出剧烈的电性横向变化，２．具有极为发育的特殊增厚的壳内高导层地段，３具有高、低热流的过渡边缘地带，４．具有与活动大断裂相交汇的隐伏断裂构造的局部地段。近年来该区地震活动增强，数次中强地震都发生在具有上述特征的地带中，该地带将可能是祁连山中段的强震－大震潜在震源危险区。这四种电性特征将是潜在震源区的判别标志。     

Parallel finite element algorithm with domain decomposition for three-dimensional magnetotelluric modelling

We present a new finite element (FE) method for magnetotelluric modelling of three-dimensional conductivity structures. Maxwell's equations are treated as a system of first-order partial differential equations for the secondary fields. Absorbing boundary conditions are introduced, minimizing undesired boundary effects and allowing the use of small computational domains. The numerical algorithm presented here is an iterative, domain decomposition procedure employing a nonconforming FE space. It does not use global matrices, therefore allowing the modellization of large and complicated structures. The algorithm is naturally parallellizable, and we show results obtained in the IBM SP2 parallel supercomputer at Purdue University. The accuracy of the numerical method is verified by checking the computed solutions with the results of COMMEMI, the international project on the comparison of modelling methods for electromagnetic induction.     

大地电磁一维连续介质反演的曲线对比法

根据电磁波的传播特性,把视电阻率随周期变化的曲线转化为电阻率随深度变化的曲线,并以此作为初始反演的地电模型．通过初始地电模型得到的视电阻率曲线与真实模型的视电阻率曲线的对比,对初始地电模型的电阻率值进行校正．校正后的反演模型的视电阻率曲线与真实模型的视电阻曲线的拟合程度有所提高．然后如此反复进行多次校正,获得与真实模型更为接近的反演模型,反演的拟合误差一般小于ｌ％．模型试验和实际例子表明,该方法的拟合程度优于Ｂｏｓｔｉｃｋ法．     

大陆岩石圈导电性的研究方法

随着地球科学的进展 ,大陆岩石圈导电性结构的研究越来越引起人们的重视 ,而超宽频带大地电磁测深则是目前用于探测岩石圈导电性最有效的地球物理先进技术。它把现代性能优良的宽频带大地电磁系统 (MT 2 4NS或V5 2 0 0 0 )与长周期智能化大地电磁系统 (LIMS)配套使用 ,采集地面上频率范围为n× 10 -4～n× 10 2 Hz的天然电磁场信号 ,并通过一系列数据处理和反演计算 ,获得深达下地壳和上地幔的地下导电性结构模型。此模型不仅可以提供有关岩石圈地质构造轮廓的信息 ,更重要的是可以间接反映现今地下深部的热结构特征和物质状态分布特点。多年来 ,应用超宽频带大地电磁测深对青藏高原岩石圈导电性结构的研究表明 ,高原的中、下地壳确实是良导电性的 ;它可能说明西藏地壳中普遍存在岩石的局部熔融 ,或地热流体。沿着应县—商河剖面的大地电磁测深研究结果从电性的角度证明了太行山山前断裂为一组向东倾斜的深断裂 ,华北地区岩石圈以其为界划分为东、西两区 ;东区为低阻区 ,与构造活动区的岩石圈导电性特征相符 ;西区为高阻区 ,表现出稳定大陆区岩石圈导电性结构的特点。     

用边界单元法模拟二维地形对大地电磁场的影响

本文用边界单元法模拟二维地形对大地电磁场的影响.首先用格林公式将二维大地电磁的边值问题转变成积分方程,然后用边界单元法解积分方程,得到地形上的大地电磁场和它的法向导数,由此可计算电阻率.与有限单元法相比,本方法剖分后的地形与实际地形的拟合程度高,向计算机输入原始数据的工作简单,可以在微机上计算地形引起的视电阻率.本方法适合于在野外生产现场进行大地电磁法的地形改正.计算表明,地形对H_x型波的影响比对E_x型波严重得多;随着周期的增长,地形对H_x型波的影响变得稳定.     

由实测资料看天然电磁场的平面波特性

在大地电磁测深资料处理中发现 ,同时观测的多个测点的磁道信号在一定距离范围内互换 ,电磁测深曲线形态变化微小。通过实际资料处理结果研究了天然电磁场平面波的某些量化特性 ,其结果显示 ,在一定范围内的天然磁场信号满足平面波的特性。该研究结果为大地电磁测深勘探中合理布设磁场探测器 ,提供了实际依据