坑（井）—地电阻率成像技术数据处理方法研究.doc

坑（井）—地电阻率成像技术数据处理方法研究( 谢成良1) 魏文博1),2),3) 叶高峰1),2),3) 金胜1),2),3) 景建恩1),2),3) 张乐天1) 董浩1) 张帆1) 1) 中国北京 100083 中国地质大学（北京）地球物理与信息技术学院 2) 中国北京 100083 中国地质大学地下信息探测技术与仪器教育部重点实验室 3) 中国北京 100083 地质过程与矿产资源国家重点实验室 摘 要 开发坑（井）—地电阻率成像技术数据处理软件包，对分布式地面—坑道电磁信号接收仪记录的时间序列进行分析处理，最终得到相应的视电阻率和相位参数。实际资料处理结果表明，在远区情况下，需考虑信噪比水平，选择合适的计算方法。当信噪比很低或较高时，可选择根据卡尼亚电阻率定义式直接比值的方法，该方法计算简便、效率高；当信噪比一般时，采用最小二乘或Robust等估算方法，能在一定程度上改善数据处理质量。 关键词 坑（井）—地电阻率成像技术；卡尼亚电阻率；数据处理软件；功率谱矩阵；估算频点 0 引言 目前矿产勘查面临向深部拓展找矿空间的难题，常规电法勘探技术显出有效探测深度不足、分辨率低和抗干扰能力差等缺陷。针对常规电法勘探技术的缺陷，课题组提出结合矿山或勘探区现有的钻孔和坑道工程，研发大功率坑（井）—地电阻率成像技术和仪器。其中，研发一套坑（井）—地电阻率成像技术的数据处理软件包（下文简称为软件包）是课题研究的主要内容之一。分析认为，该软件包的实现可借鉴可控源音频大地电磁法（controlled source audio-frequency magnetotelluric, CSAMT）的数据处理过程，目前研究人员对CSAMT的数据处理方法进行了广泛研究（安志国等，2006；底青云等，2008）。在远区情况下，可根据电磁场的各分量分别求出各场强单分量的视电阻率公式（黄皓平等，1992；毛先进等，1996；底青云等，2008；佟铁钢等，2009）。基于单分量的视电阻率计算方法，优点在于能够得到较多参数信息，有助于综合解释能力的提高，但缺点在于单道观测的资料受干扰因素较多（底青云等，2008）。为此，可延用大地电磁数据的处理方式，通过相互正交的电、磁场分量的比值求得卡尼亚电阻率（Cagniard, 1953；底青云等，2008）。在求得卡尼亚电阻率后，可将受到场源影响的低频部分截去，得到未受场源影响的数据参与反演工作，也可根据观测参数，对近区和过渡区的数据进行近场校正（周茂军等，1993；邝广，1996；张巍等，2008）。此外，研究人员对全区视电阻率的计算做了大量的工作（汤井田等，1994a，1994b；田继枫，2011）。全区视电阻率的计算可有效运用过渡区、甚至是近区的数据，是CSAMT数据处理的重要研究方向之一。然而，由于其定义形式一般较为复杂，简明的求解算法尚待研究。综上，根据“坑（井）—地电阻率成像技术及系统研制”课题的研究内容要求，针对分布式地面—坑道电磁信号接收仪的数据记录格式，在总结目前数据处理研究现状基础上，研究开发一套针对远区数据，集数据分析、数据处理、结果显示等功能于一体的坑（井）—地电阻率成像技术数据处理软件包，从而能够根据工区干扰情况，选择合适的计算方法，得到具有较好数据质量，并能较真实反映地下介质垂向电性变化的视电阻率、相位信息。 1 方法原理 1.1 卡尼亚电阻率定义 采集两对相互正交的电磁场分量，即和，和。通过计算二者比值的方式求得卡尼亚电阻率，电场单位使用mv/km,磁场单位使用nT。以xy模式为例，（1）式为xy模式的视电阻率和相位表达式。 这里T为信号的周期，为视电阻率，为相位，其中，分别为与的相位。通常，对于人工源频率域电磁法而言，在信噪比较高的情况下，可以在频率域对不同频率的电磁场信号计算视电阻率和相位；但在工区干扰水平较高，信噪比一般的情况下，通过采集多次叠加信号，计算信号的功率谱，最后采用最小二乘、Robust等估算算法计算每个测点不同频率的阻抗信息，将一定程度上增加资料处理的可靠性。 1.2 视电阻率计算流程 远区卡尼亚电阻率的估算主要包括以下步骤：（1）场分量旋转。当不满足测量电偶极平行或垂直于发送电偶极的条件时，需要进行场分量旋转；（2）傅里叶系数的计算。对原始信号进行加窗处理，选择合适的窗口长度，进行傅里叶变换；（3）傅里叶系数的校正。设仪器的A/D格值为，电场增益为，电场极距为，磁传感器灵敏度为，磁场增益为，又设对于频率，标定系数为，磁场分量傅里叶系数为，电场分量傅里叶系数为，则对、的校正如（2）式所示（Fiona等, 2005），其中，，分别为校正后磁、电场分量的傅里叶系数。（4）功率谱矩阵的计算。得到经校正的傅里叶系数后，需确定估算频点