利用地磁台站数据和站间传递函数估算长周期大地电磁测深的本地阻抗

大地电磁测深的有效信号为天然电磁场,容易受到人文电磁噪声干扰,如何提高大地电磁资料的信噪比,已经成为制约大地电磁应用和发展的一个重要因素.本文利用地磁台站数据,首先将地磁台站数据转化为长周期大地电磁测深所采集的磁场数据格式,利用远参考方法估算长周期大地电磁测点与地磁台站之间的站间传递函数,包括长周期大地电磁测点电场与地磁台站水平磁场之间的拟阻抗,以及长周期大地电磁测点水平磁场与地磁台站水平磁场之间的水平磁场传递函数,将拟阻抗与水平磁场传递函数相乘便可以得到长周期大地电磁测点本地阻抗的无偏估计.实测数据处理结果表明,该方法可以提高本地阻抗估算的数据质量,减小阻抗估算误差.     

瞬变电磁测深早期数据的修正

通过时间-频率转换关系,TEM数据可以转换成平面波场测深数据,从而可以对TEM资料进行拟平面波场处理解释.在对瞬变电磁视电阻率数据向平面波场测深视电阻率数据转换时,发现由于瞬变电磁使用晚期计算公式及装置问题,使测深曲线早期数据发生畸变.文中建立了视电阻率曲线进入晚期所满足的关系式,从理论上给出不同情况下瞬变电磁测深视电阻率曲线进入晚期的临界点.以瞬变电磁、大地电磁、CSAMT为例,对大量的模型进行正演计算,对计算结果进行对比分析,建立不同地表电性结构、不同时间延迟情况下,瞬变电磁早期数据误差的校正量板.     

不同参数条件下远参考大地电磁处理效果对比分析

研究了在不同噪声强度、不同参考道、噪声仅存在基站或参考站时,远参考处理对测深曲线的影响.通过对标准时间序列进行仿真模拟,发现基站含不同信噪比的非相关噪声,远参考大地电磁对测深曲线提升效果不一:高信噪比(SNR=10.58 d B)条件下,远参考大地电磁能够较完美的恢复测深曲线,低信噪比(SNR=0.2 d B)条件下,远参考大地电磁对测深曲线质量提升效果有限;相同噪声仅存在于基站或者参考站对资料处理效果有区别:噪声仅存在于参考站的处理效果比仅存在于基站的效果要好;参考站无噪声情况下,电道或者磁道作为参考道对处理结果影响不大;基站与参考站之间的磁场相关性、基站磁场噪声与参考站磁场的非相关性与远参考处理效果关系密切.通过对实际资料的计算也有上述结论.     

从瞬变电磁测深数据到平面电磁波场数据的等效转换

中心回线源装置的瞬变电磁场是一种涡流场,在地下主要以扩散形式传播.平面电磁波在地下介质中传播时具有反射、折射等电磁学特性,为了借鉴平面电磁波场测深成熟的解释方法,文中开展了从瞬变电磁涡流场测深数据到平面波场数据的快速等效转换的研究.通过大量的理论模型的正演计算、曲线对比,误差统计,对瞬变电磁场、平面波场在地下传播机制特性的分析,建立了一种由瞬变电磁测深视电阻率数据向平面波场测深视电阻率数据转换的时间-频率等效对应关系.利用这一关系式,把扩散场数据快速转换成等效平面波场测深数据. 为进一步进行瞬变电磁拟平面波成像解释打下理论基础.     

庐枞矿集区大地电磁法的远参考效果研究

庐枞矿集区大地电磁测深(MT)数据受强烈的工业、通讯、矿山、民用等电磁干扰,致使资料数据离差大、曲线形态畸变、可信度差,本文从大地电磁测深张量阻抗计算表达式出发,分析了噪声对张量阻抗估计产生的影响.选用本区质量较好的MT数据作为"远参考点",应用远参考方法处理庐枞MT数据,研究其对卡尼亚电阻率-相位曲线的改善效果,以及对数据相干度和信噪比的改善效果,并提出磁场信号的相关性是远参考处理的关键因素.     

可控源音频大地电磁测深法勘查深度研究

可控源音频大地电磁测深法的勘查深度问题在理论研究和实际工作中都得到了广泛的关注,本文分析了水平层状大地电磁正反演方法的基本原理,说明了可控源音频大地电磁测深法在深部勘查方面具有可靠性.通过建立两个三层地电模型,并对吉林某地区的CSAMT实测数据进行了带源反演成像,与已知钻井资料对比分析表明,可控源音频大地电磁测深法对地下1000~2000 m深度的反演较为准确,能够反映地下电性的垂向变化.     

试论大地电磁测深法在地震前兆探索中的进一步应用

从大地电磁测深理论分析和实际资料数字试验出发,我们探讨了有关提高大地电磁测深法监测深部电性变化前兆效能的某些问题.研究结果表明,在所定义的几种大地电磁视电阻率中,尚未加以利用的阻抗实部视电阻率具有某些独特的优点,可能成为监测深部电性变化的一个较好的参数指标;应根据各记录道误差的分配,计算过程误差的传播和累加对不同大地电磁响应函数产生不同的影响,选择相对精度较高的参数作为主要监测参量;而识别深部电性变化的多参数综合判别方法应是基本的分析方法.实际观测到的反映地震前兆的某些大地电磁资料证实了上述论证.      

大地电磁场源效应特征分析及其校正研究

大地电磁法一直以来基本都遵循场源为平面电磁波的假设,这种假设普遍适用于中高频的大地电磁研究,但在具有显著非平面波场特性的区域或开展长周期大地电磁研究时,则对场源平面波的理论基础带来了挑战,若不满足理论假设,将会导致深部探测的可靠性降低,勘探目标定位风险增高.本文通过选择线状、片状及非规则场源模型来模拟极光电集流、赤道环电流和横向波动的波数场源,开展了大地电磁场源效应的模拟研究,对比分析了非平面波模式和平面波模式的响应差别.模拟表明大地电磁场源效应强度与模型平均电阻率正相关,与频率和测点与源水平距离负相关,且与场源的高度、宽度和横向波长等参数均有关系,非平面波模式和平面波模式正演响应在一些情况下具有超过10％的差异.基于场源效应特征,归纳并提出了五种场源效应校正方法,一是计算极限模型来截断场源效应影响频段的频域截断法,二是对多个同步观测的邻近测点进行平均的测点平均法,三是延长测点观测时间的时间域延长法,四是剔除较强垂直磁场分量所对应场源效应较强段数据的时间域剔除法,五是直接考虑场源的反演,五种方法均能一定程度上降低大地电磁场源效应的影响.对于受场源效应影响的大地电磁数据,若直接采用平面波方式进行反演则可能在深部产生较大畸变,若在获知场源相关参数的基础上,可开展考虑场源的大地电磁反演来消除场源效应影响.     

南黄海海底大地电磁测深试验研究

针对南黄海盆地深层地震地质条件差,地震波受屏蔽,能量严重衰减,在海洋反射地震资料中很难识别反映中-古生界的反射地震波组,不能满足南黄海盆地前第三系油气资源前景评价的需要等问题,于2006年5月,开展了南黄海海域海底大地电磁测深试验,在盆地内五莲斜坡和胶州凹陷区布设了3个海底大地电磁测深点,取得了可靠的海底观测数据.本文详细阐述了海底大地电磁测深与常规大地电磁测深在数据处理方法上的差别,论述了海洋电磁噪声对海底大地电磁场观测的影响,提出了海底大地电磁场实测资料“水平姿态”校正、“方位”校正和“海洋电磁噪声”校正的方法.通过对试验观测资料的处理和反演,并结合地质、物性资料进行综合分析,结果表明在南黄海盆地海底以下深部还存留有古生代地层,这对于南黄海前第三系含油气前景的评价具有重要意义.试验结果也证明了海底大地电磁测深技术在研究古生代残留盆地方面的问题可以发挥相当好的作用.     

MMSO2E大地电磁测深仪资料处理解释软件系统研制

研制了一套与MMSO2E大地电磁测深仪配套的资料处理解释软件系统,支持MMSO2E大地电磁测深仪在野外的条件下完成常规的现场资料处理解释。经对比表明,本软件系统的实际使用效果并不劣于国外同类软件,且由于研制过程中注意吸收了近年来大地电磁测深方法研究的成果,而使之具有自己一定的特色。     

Coast effect distortion of marine magnetotelluric data: Insights from a pilot study offshore northeastern Japan

We report on strong coast effect distortions observed for broadband marine magnetotelluric (MT) data collected on the forearc offshore northeastern Japan. Eight days of horizontal electric and magnetic fields recorded at eight seafloor stations and the horizontal magnetic fields from a land remote station were processed with a robust multiple-station algorithm, yielding good MT responses and inter-station transfer functions at periods of 7–10,000 s. Transverse electric (TE) mode responses have cusps in apparent resistivity and negative phases at periods around 1000 s, while the transverse magnetic (TM) mode responses are galvanically depressed below the TE responses. An analysis of inter-station transfer functions confirms that the apparent resistivity cusps are a magnetic field, rather than electric field, phenomenon, consisting of an amplitude minimum and rapid phase change around a characteristic frequency. Poynting vectors for a TE coast effect model study illustrate that the anomalous phases are associated with energy diffusing back up to the seafloor from below, after being turned around from its usual downward propagating trajectory by inductive coupling between the conductive ocean and the resistive seafloor along the continental margin. We show that the characteristic frequency and position of the TE mode apparent resistivity cusps are determined by a relatively simple combination of the electrical resistivity of the seafloor, the depth of the ocean, and the distance from the coastline. By including coastlines and bathymetry in 2D inversion, we recover the seafloor conductivity structure along the forearc, demonstrating that broadband data can constrain the thickness of conductive forearc sediments and the underlying high resistivity associated with the mantle wedge and subducting oceanic lithosphere.     

CSAMT单分量数据解释方法

可控源音频大地电磁法(CSAMT)一直沿用大地电磁法(MT)的办法,通过计算电场分量与磁场分量的比值,求取卡尼亚视电阻率.而CSAMT场源已知,电场分量和磁场分量都与地下电阻率存在一定的关系,可以单独采用CSAMT电场分量或者磁场分量提取地下介质的视电阻率.本文通过分析电场分量与磁场分量的数据特性,提出利用CSAMT电场单分量数据进行视电阻率的计算,用改进的广义逆矩阵反演方法,使初始模型中的地电层数等于频道个数,克服了以往反演计算中层数较少的问题;实现全场区电场分量视电阻率曲线的拟合反演.同时对单分量视相位计算方法进行分析,结合山西大同地区积水采空区探测及数据解释结果,论证本文提出的单分量解释方法的有效性.     

震源深度范围内视电阻率随时间变化的观测试验及其与地震活动的关系

本文讨论了青海祁连俄博测点和甘肃山丹霍城测点1979年7月——1986年10月间的大地电磁测深重复观测试验结果。测点所在地区的电磁干扰背景较小、地震活动性较强。文中分析了观测记录和资料处理中的误差。结果表明,以初始值为基准,俄博测点1984年两次复测的视电阻率曲线在部分频带范围内出现过15——30%左右的系统变化,与此期间测区附近发生的一些M5的地震活动有着一定的对应关系;1985年复测的视电阻率曲线又恢复至初始值,而此期间测区附近没有M5及其以上的地震活动,1986年霍城测点和俄博测点复测的视电阻率曲线在较宽的频带范围内发生了40——70%左右的系统变化,与此期间震中距为60km左右的青海门源M=6.7,5.7和5.0地震群密切相关。      

THE DISPERSION RELATION, OF EARTH''S ELECTROMAGNETIC RESPONSE FUNCTIONS AND THE JOINT INVERSION OF DATA

On the basis of the dispersion relations of MT field, the necessity and applied prospects of the joint inversions using a pair of MT response functions which are correlative with the dispersion relations, are infered. A filter coefficient algorithm is made, with which the corresponding impedance phase data can be estimated using a set of apparent resistivities. The tests for the observed MT data show that when comparing the impedance phase estimated using the dispersion relation with the ob served phase, it can be checked whether the dispersion relation between observed apparent resistivity and phase data is satisfied or not, and that the use of the phase data corrected using the dispersion relation in the joint inversion is advantageous to obtain more confident results. It is shown that joint inversions are more advantageous than single parameter inversions, and that in the most case the joint inversion using the apparent resistivities of impedance real and imaginary parts is more advantageous than the jointinversion using the normal apparent resistivity and impedance phase. The existence of the dipersion relations between the ratio apparent resistivity and corresponding impedance phase of the orthogonal electric and magnetic field horizontal Components in the frequency EM sounding with horizontal electric dipole(FEMS) are discussed, the better effect of the joint inversion using the pair of EM response functions is obtained. The problems on the one-dimensional joint inversion for the MT and FEMS apparent resistivities, for which the observed frequency bands partly overlape each other, are studied. It is shown that this joint inversion is applicable and effective:the joint inversions of the practical data for two kinds of EM methods at two sites give the results well corresponding to the drilling data. The simulated MT inversions for the data of two kinds of EM methods are made, and more confident results also are obtained.     

南北地震带南段川滇黔接壤区电性结构特征和孕震环境

本文利用大地电磁测深数据,对穿过兰坪—思茅地块和川滇菱形地块以及进入扬子地块的云南兰坪—贵州贵阳大地电磁测深剖面展开了深部电性结构研究.采用大地电磁数据处理分析以及反演技术,对观测资料进行了由定性到定量全面地分析,通过二维非线性共轭梯度反演得到了沿剖面的较为详细的地壳上地幔电性结构,结合其他地质和地球物理资料的分析,对该剖面的二维电性结构进行解释,确定了主要断裂带和边界带的位置以及深部延伸情况,同时确定了壳内低阻层的分布位置,最后进行了区域动力学和孕震构造环境的探讨.研究表明:剖面壳幔电性结构分块性特征与区域地质构造分布特征基本一致,不同地块的电性结构存在显著差异,其中川滇菱形地块的结构相对复杂,上地壳的电性结构为高低阻相间分布特征,电阻率的突变带与地表断裂具有很好的对应关系;兰坪—思茅地块存在中上地壳低阻层,川滇菱形地块中西部存在下地壳低阻层,川滇菱形地块东部和华南地块西部存在中上地壳的低阻层;川滇菱形地块中部攀枝花附近的低阻层埋深最深,而华南地块西部会泽附近的低阻层埋深则最浅;兰坪—思茅地块和川滇菱形地块的中下地壳的低阻层可能与青藏高原物质的东南逃逸有关;华南块体的宣威以东的下地壳不存在低阻层,华南块体下地壳和上地幔的电阻率较高;攀枝花附近的高阻体可能是峨眉山玄武岩喷发导致底侵作用及幔源物质上侵的结果.     

基于站间天然电磁场单位脉冲响应的大地电磁时间序列去噪方法

针对天然大地电磁场信号在人文活动密集地区易受噪声干扰的问题,本文提出利用两个同步测点天然电磁场时间序列之间的单位脉冲响应,合成本地点受干扰时段的数据,从而去除大地电磁噪声.首先,选择高信噪比时段的数据,采用最小二乘法,估算本地点与参考点之间的单位脉冲响应,再根据卷积定律,结合参考磁场合成本地点的磁场和电场.最后用合成数据替换含噪声时段数据,实现时间域去噪.实测高信噪比数据和含噪数据的处理结果表明,该方法可以高精度合成本地点磁场与电场信号,有效去除本地点电场和磁场噪声,包括相关噪声,提高大地电磁数据质量.     

A data variance technique for automated despiking of magnetotelluric data with a remote reference

The magnetotelluric method employs co‐located surface measurements of electric and magnetic fields to infer the local electrical structure of the earth. The frequency dependent ‘apparent resistivity’ curves can be inaccurate at long periods if input data are contaminated – even when robust remote reference techniques are employed. Data despiking prior to processing can result in significantly more reliable estimates of long period apparent resistivities. This paper outlines a two‐step method of automatic identification and replacement for spike‐like contamination of magnetotelluric data; based on the simultaneity of natural electric and magnetic field variations at distant sites. This simultaneity is exploited both to identify windows in time when the array data are compromised as well as to generate synthetic data that replace observed transient noise spikes. In the first step windows in data time series that contain spikes are identified according to an intersite comparison of channel ‘activity’– such as the variance of differenced data within each window. In the second step, plausible data for replacement of flagged windows are calculated by Wiener filtering coincident data in clean channels. The Wiener filters – which express the time‐domain relationship between various array channels – are computed using an uncontaminated segment of array training data. Examples are shown where the algorithm is applied to artificially contaminated data and to real field data. In both cases all spikes are successfully identified. In the case of implanted artificial noise, the synthetic replacement time series are very similar to the original recording. In all cases, apparent resistivity and phase curves obtained by processing the despiked data are much improved over curves obtained from raw data.     

Prospects of magnetotelluric sounding in some large structures in the Arctic ocean

In the work, the results of modeling of electric resistivity distribution at the Gakkel and Alpha-Mendeleev underwater ridges are given. Based on the seismostratigraphic models, the apparent resistivity curves at a few points of the profile at the crossing of the ridge were constructed. The advantages and disad-vantages of the MT sounding method are presented. The results obtained can be used for future experiments.     

The inversion of the data of magnetotelluric sounding and electromagnetic frequency sounding with horizontal electric dipole

On the basis of the dispersion relation of magnetotelluric response functions (MTRF), a filter coefficient algorithm has been made, with which the corresponding impedance phase data can be estimated using a set of apparent resistivity data. The tests of theoretical models and observed magnetotelluric (MT) data show that this algorithm is effective. Comparing the impedance phase estimated using dispersion relation with the observed phase, it can be checked whether the dispersion relation between the observed apparent resistivities and phase data was satisfied. The use of phase data corrected using the dispersion relation in the joint inversion for MT impedance is advantageous to obtain more reliable inversion results. The problems on the one-dimensional joint inversion for the (MT) apparent resistivity and the apparent resistivity of the frequency electromagnetic sounding (FEMS) with horizontal electric dipole, whose observed frequency bands are linked up each other, are studied. The observed data of two kinds of electromagnetic (EM) methods at two sites are used to inverse, the comparison with the drilling data show the results are more reliable. To supply the phase data of FEMS using the dispersion relation, for the apparent resistivity-phase data and impedance real part-imaginary part apparent resistivities of two kinds of EM methods the imitated MT joint inversions are made, and more similar results also are obtained. The Chinese version of this paper appeared in the Chinese edition ofActa Seismologica Sinica,15, 91–96, 1993. The projects sponsored by the Chinese Joint Seismological Science Foundation.     

地震监测人工源极低频电磁技术(CSELF)新试验

2009年利用大功率人工源极低频电磁技术(CSELF)在位于华北、华南、东北、西北和西南的12个地震台站和几个流动观测点进行了连续30天新的观测试验.结果表明,利用CSELF技术可以在1700 km之外测量到人工源电磁场信号,计算得到的电、磁场功率谱密度和视电阻率与天然源信号相比,抗干扰能力更强,观测信号更稳定,特别有利于识别和捕捉地震等诱发的电磁异常现象,在地震预测监测中具有很大的研究应用潜力.试验还发现,各地震台站和流动测点观测的CSELF信号的强度,与台站或者测量的电、磁场分量相对于发射源的距离、方位有明显的关系;发现在离开发射源的相对近区,场强随距离的衰减,比在相对远区更快;发现接收磁场信号比电场信号较易获得更高信噪比的数据.