The electrical structure of the Slave craton

The Slave craton in northwestern Canada, a relatively small Archean craton (600×400 km), is ideal as a natural laboratory for investigating the formation and evolution of Mesoarchean and Neoarchean sub-continental lithospheric mantle (SCLM). Excellent outcrop and the discovery of economic diamondiferous kimberlite pipes in the centre of the craton during the early 1990s have led to an unparalleled amount of geoscientific information becoming available.

Over the last 5 years deep-probing electromagnetic surveys were conducted on the Slave, using the natural-source magnetotelluric (MT) technique, as part of a variety of programs to study the craton and determine its regional-scale electrical structure. Two of the four types of surveys involved novel MT data acquisition; one through frozen lakes along ice roads during winter, and the second using ocean-bottom MT instrumentation deployed from float planes.

The primary initial objective of the MT surveys was to determine the geometry of the topography of the lithosphere–asthenosphere boundary (LAB) across the Slave craton. However, the MT responses revealed, completely serendipitously, a remarkable anomaly in electrical conductivity in the SCLM of the central Slave craton. This Central Slave Mantle Conductor (CSMC) anomaly is modelled as a localized region of low resistivity (10–15 Ω m) beginning at depths of 80–120 km and striking NE–SW. Where precisely located, it is spatially coincident with the Eocene-aged kimberlite field in the central part of the craton (the so-called “Corridor of Hope”), and also with a geochemically defined ultra-depleted harzburgitic layer interpreted as oceanic or arc-related lithosphere emplaced during early tectonism. The CSMC lies wholly within the NE–SW striking central zone defined by Grütter et al. [Grütter, H.S., Apter, D.B., Kong, J., 1999. Crust–mantle coupling; evidence from mantle-derived xenocrystic garnets. Contributed paper at: The 7th International Kimberlite Conference Proceeding, J.B. Dawson Volume, 1, 307–313] on the basis of garnet geochemistry (G10 vs. G9) populations.

Deep-probing MT data from the lake bottom instruments infer that the conductor has a total depth-integrated conductivity (conductance) of the order of 2000 Siemens, which, given an internal resistivity of 10–15 Ω m, implies a thickness of 20–30 km. Below the CSMC the electrical resistivity of the lithosphere increases by a factor of 3–5 to values of around 50 Ω m. This change occurs at depths consistent with the graphite–diamond transition, which is taken as consistent with a carbon interpretation for the CSMC.

Preliminary three-dimensional MT modelling supports the NE–SW striking geometry for the conductor, and also suggests a NW dip. This geometry is taken as implying that the tectonic processes that emplaced this geophysical–geochemical body are likely related to the subduction of a craton of unknown provenance from the SE (present-day coordinates) during 2630–2620 Ma. It suggests that the lithospheric stacking model of Helmstaedt and Schulze [Helmstaedt, H.H., Schulze, D.J., 1989. Southern African kimberlites and their mantle sample: implications for Archean tectonics and lithosphere evolution. In Ross, J. (Ed.), Kimberlites and Related Rocks, Vol. 1: Their Composition, Occurrence, Origin, and Emplacement. Geological Society of Australia Special Publication, vol. 14, 358–368] is likely correct for the formation of the Slave's current SCLM.     

苏北地区大地电磁测深采集方法技术研究

根据苏北地区大地电磁勘探环境特点进行了电极距长度、磁场位置和采集时长的试验工作，确定了大地电磁测深勘探技术参数设置下限，电极距长度不能小于50m，采集时间要大于8h，磁距应控制在2km范围内。     

MT时间序列的FIR数字滤波分析

大地电磁测深法(MT)野外观测得到的原始资料是包含多种频谱的时间序列,这种时间序列中常包含不需要的信息和50 Hz及其谐波的干扰。结合FIR数字滤波技术与MT时间序列的特征,提出了一种基于FIR数字滤波技术的MT时间序列去噪方法,讨论了FIR数字滤波器的原理和实现,分析了理论数据通过FIR数字滤波的带通和陷波处理,探讨了对实测时间序列的带通和陷波所取得的效果。     

利用有限单元法模拟二维MT倾子响应

在二维地电模型中,考虑到大地电磁测深(MT)倾子响应依赖于横电(TE)模式的数值模拟问题,笔者采用矩形网格单元和双二次插值对MT的倾子响应进行求解。首先,给出了二维地电条件下的边值问题和变分问题,并通过有限单元法对模型进行单元剖分、插值、积分和整体合成;其次,通过求解复系数方程组得到了每个节点的电场值;最后,采用差分方法求得电场沿横向和垂向的偏导数值,从而计算出二维MT的倾子响应。通过对2个二维地电模型的倾子响应进行数值模拟,获得了倾子响应的实部、虚部和振幅,结果很好地反映了不均匀体的横向分布情况。     

可控源音频大地电磁法深部找矿实验效果

为了考查电磁法在探测深部金属矿及地质构造方面的实际效果,用可控源音频大地电磁法分别在埋深为500 m和1 100 m的两个已知隐伏铁、铜矿体上进行了实验工作。实验结果表明地电异常显示了矿体、断层以及岩层接触带的分布。在矿体埋藏深度小于500 m、矿层较厚时,地电异常十分明显;若矿体埋藏深度大于1 000 m,地电异常主要反应了控矿、容矿构造。     

Magnetotelluric investigations for imaging electrical structure of Garhwal Himalayan corridor,Uttarakhand, India

Magnetotelluric investigations have been carried out in the Garhwal Himalayan corridor to delineate the electrical structure of the crust along a profile extending from Indo-Gangetic Plain to Higher Himalayan region in Uttarakhand, India. The profile passing through major Himalayan thrusts: Himalayan Frontal Thrust (HFF), Main Boundary Thrust (MBT) and Main Central Thrust (MCT), is nearly perpendicular to the regional geological strike. Data processing and impedance analysis indicate that out of 44 stations MT data recorded, only 27 stations data show in general, the validity of 2D assumption. The average geoelectric strike, N70°W, was estimated for the profile using tensor decomposition. 2D smooth geoelectrical model has been presented, which provides the electrical image of the shallow and deeper crustal structure. The major features of the model are (i) a low resistivity (<50Ωm), shallow feature interpreted as sediments of Siwalik and Indo-Gangetic Plain, (ii) highly resistive (> 1000Ωm) zone below the sediments at a depth of 6 km, interpreted as the top surface of the Indian plate, (iii) a low resistivity (< 10Ωm) below the depth of 6 km near MCT zone coincides with the intense micro-seismic activity in the region. The zone is interpreted as the partial melting or fluid phase at mid crustal depth. Sensitivity test indicates that the major features of the geoelectrical model are relevant and desired by the MT data.     

三维地形大地电磁场的边界元模拟方法

提出了一种用边界元法计算大地电磁场三维地形影响的数值模拟方法.首先用矢量积分理论和电磁场边界条件, 将上半空间(空气)和下半空间(地下介质)两个区域电磁场边值问题变为仅对地形界面的两个矢量面积分方程, 其中一个计算磁场, 称磁场方程; 另一个计算电场, 称电场方程.然后将对地形界面的积分剖分为一系列的三角单元积分.在三角单元积分中, 假设单元中电磁场为水平均匀大地空间电磁场与地形影响的迭加, 并假设地形影响为常项, 这样既保证了计算精度又使得计算方法简便.通过分解和计算, 每一个矢量面积分方程分解为对应3个坐标方向的3个常量线性方程, 这些线性方程组成了对角占优的线性方程组, 可用SSOR方法求解.文中给出了2个三维地形上大地电磁视电阻率曲线的计算结果.      

地下电性非均匀结构对地电场信号的影响

在数值模拟的基础上探讨了地下电性非均匀性问题对地电场信号的影响程度。这些影响主要包括地下电性非均匀体对天然电磁感应电场的影响和对产生于震源区附近的异常地电场的影响。结果表明, 地下电性不均匀体的电性特征与周围介质的电性特征差别越大, 则其对地电场的影响也会越大; 地下电性不均匀体的空间范围越大, 则其对地电场的影响也越大。在孕震区作为可能的电荷源模型中, 地下电性不均匀体的存在对电荷源所形成的异常地电场信号在地表的表现形式也有很大的影响。在地表所观测到的地电场异常空间分布不均匀性与地下电性不均匀体的空间位置有着密切关系。     

吉林龙岗火山区形变特征与隆升机理分析

吉林龙岗火山区是第四纪活动火山,具有潜在喷发危险。通过GNSS和水准资料分析区域三维地壳运动,得出2011—2019年相对欧亚板块水平运动以东南向为主,速度＜10 mm/a。敦化-密山断裂以东受日本2011年“3·11地震”影响强烈,现以拉张运动为主;近年来水准资料揭示火山区垂直运动以隆升为主。2014—2019年的InSAR资料显示隆升集中在靖宇一带。结合MT剖面反演得到的深部电性结构,龙岗火山区西侧高阻体分布在深度18 km以上,金龙顶子火山下方最浅。中部为早期喷发形成的火山,下方高阻体分布在深度40 km以上,东侧抚松一带分布在深度约20 km以上,地壳范围内的高阻结构表明岩浆已固结。高阻结构层可分辨出断层两侧电阻的差异性。高阻体下方存在大规模低阻结构,推测为中下地壳岩浆系统。金龙顶子火山深度10 km以下的低阻结构可能为岩浆通道,并与中下地壳岩浆系统相连。东侧区域岩浆平均深度约30 km,相对较浅。龙岗火山区幔源物质的上涌及间断性的向上运移引起地壳隆升及地震活动。     

大地电磁偏移成像技术的研究进展及其解决复杂地质问题的能力

总结和分析了大地电磁场成像技术的研究现状 ,在Zhdanov等对偏移电磁场的研究和成像技术基础上 ,提出了改进的有限差分法大地电磁场偏移成像技术 ,使差分方程精度和成像分辨率有了明显的提高 ,将前人大地电磁场偏移成像的模型试验研究推进到适用于复杂地电模型和实际资料处理解释的阶段 ,并通过松辽盆地大地电磁测深资料的偏移成像处理和解释展示了该方法的有效性和实用性