一种适用于强电磁干扰地区的深部采空区探测方法与流程

本发明涉及一种适用于强电磁干扰地区的深部采空区探测方法。

背景技术：

近年来，由于不明采空区所引起的塌陷、突水、有毒气体泄漏等事故频发，已经严重影响到城市建设、资源开发和环境保护等领域，甚至还威胁到人员和财产安全，带来了巨大的损失。因此，如何能够准确探测煤层采空区的分布范围，具有重要的现实意义。

对此国内外学者将物探方法应用于采空区预测，目前常用的有浅层高分辨地震勘探、高密度电法勘探、探地雷达精细探测、瞬变电磁和eh4等方法，均已取得良好的预测结果。

可控源音频大地电磁发(controlled-sourceaudiomagnetotellurics，简称csamt)具有勘探深度范围大、抗干扰能力强、高阻屏蔽作用小、分辨率高、工作效率高等优点，被广泛应用于地热资源探测、桥梁隧道工程、古墓探测、煤矿采空区预测、水文地质、溶洞预测、有色金属矿床、石油勘探开发等领域。近年来，csamt逐渐被应用于采空区探测，并且在采空区探测时也取得了良好的效果，因此，csamt为采空区的预测提供了新的方法。

然而，大部分采空区位于煤矿以及村庄等附近，采空区分布的地方往往有较多的动力电线和高压线，因此电磁干扰异常强烈、各类随机噪声干扰，从而导致有效信号被隐藏。综上，如何确保在强电磁干扰地区能够获得高精度的csamt数据，去除浅部异常信息，同时增强深部采空区的电性特征，以便能够实现采空区的准确预测，具有重要的现实意义。

技术实现要素：

本发明的目的在于提出一种适用于强电磁干扰地区的深部采空区探测方法，以便对强电磁干扰地区因多干扰因素影响产生的电性异常进行滤除，从而提高测区采空区探测精度。

本发明为了实现上述目的，采用如下技术方案：

一种适用于强电磁干扰地区的深部采空区探测方法，包括如下步骤：

i.基于深度神经网络学习与预测方法对原始采集的csamt数据进行异常值剔除；

具体过程如下：

i.1.设计深度神经网络结构

设计两个输入层，输入参数分别为频率、电阻率；设计两个输出层，输出为0或者1，若输出为0，则代表为异常值，进行剔除；若输出为1，代表正常值，进行保留；

深度神经网络的隐藏层数设为k层，每个隐藏层的神经元的数量由公式(1)得到：

l＝2+α1(1)

式中，α1为正整数，且α1的取值范围为[1,10]，k的取值范围是[5,10]；

其中，深度神经网络的激励函数采用relu函数，深度神经网络优化算法采用adagrad算法，深度神经网络采用dropout正则化方法以防止过拟合现象；

i.2.基于样本数据曲线拟合进行标签制作，具体过程如下：

从原始采集的csamt数据中随机选取几组点号的数据作为样本数据；

已知某一点号获得的一系列csamt数据为首先将电阻率转换为对数值，则：

其中，i1代表某一点号的第i1个数据点；

代表某一点号的第i1个卡尼亚电阻率，表示第i1个卡尼亚电阻率对应的频率；

随后，对数据点通过最小二乘法进行曲线拟合，获得函数解析式为：

最后，依据曲线拟合结果进行标签制作，则标签函数：

式中，ε为误差；若输出f＝0，则第i1个卡尼亚电阻率对应的点代表为异常值，进行剔除，若输出f＝1，则第i1个卡尼亚电阻率对应的点代表正常值，进行保留；

i.3.将步骤i.2中的样本数据进行划分，80％的作为学习样本，20％的作为验证样本；

首先将学习样本输入步骤i.1设计的深度神经网络结构中，进行深度神经网络模型训练；

假设深度神经网络某一神经元的输出为：

式中，j为大于0的整数；

为深度神经网络中第k个隐藏层的第j个神经元的输出，k∈k；

为深度神经网络中第k个隐藏层的第j个神经元的输入；

f为深度神经网络的激活函数；

为第k-1层的第i个神经元到第k层的第j个神经元的连接权重；

为深度神经网络第k层第j个神经元的偏置项；

在误差反向传播中，需要对深度神经网络的预测输出和期望输出进行计算，得到深度神经网络的预测输出和期望输出之间的误差，误差计算公式为：

式中，dic为期望输出，yic为深度神经网络的预测输出，q为训练样本数目；

由公式(5)和(6)可知，误差函数ek是由各个神经元决定的；不断对深度神经网络的连接权重和偏置项进行修改，使深度神经网络的误差函数ek满足设定的误差函数阈值；

连接权重和偏置的更新过程如下：

bj^k(h+1)＝bj^k-1(h)+ηδi(8)

式中，h表示第h个训练步数，h为大于0的整数；

为第h个训练步数得到的权系数；η为学习率，η的取值范围为(0,0.05]；δi＝(di-yi)y′i，其中，y′i为yi的导数；bj^k(h)为第h个训练步数得到的偏置；

通过上述步骤，完成深度神经网络模型的训练；

随后，用检验样本对深度神经网络模型的异常值剔除效果以及稳定性进行测试，若测试结果的准确率达到精度要求，则完成模型的训练与测试；

若不能达到精度要求，重复上述训练过程，直至深度神经网络模型达到精度要求为止；

i.4.对野外的不同剖面不同点号处的数据输入训练好的模型，自动进行异常点的剔除；

ii.基于多方向时变二维中值空间滤波法对异常值剔除后的csamt数据进行静态校正；

具体过程如下：

ii.1.确定滤波窗口

选取c作为中值滤波的基准窗口，对异常值剔除后的csamt数据进行中值滤波；

位置(m,n)处的中值输出表示为：

其中，表示以c作为基准窗口(m,n)位置处的中值；

ia＝m-(c-1)/2,…,m+(c-1)/2；ja＝n-(c-1)/2,…,n+(c-1)/2；

随后引入下式，计算阈值t：

ρix,iy表示原始电阻率值，ix＝1,…,nx，iy＝1,…,ny，nx为电阻率剖面横向的点数，ny为电阻率剖面纵向的点数；再根据阈值t，分别设定时变窗口的大小：

式中，cia,ja为滤波时的应用窗口大小，α,β,γ,δ分别表示常数，且α＞β,δ＞γ；

ii.2.确定滤波方向

对于cia,ja＝2n+1为滤波窗口，当中心点位于(m,n)时，滤波函数被定义为：

其中，n为大于0的常数，w1[m,n]、w2[m,n]、w3[m,n]、w4[m,n]分别表示经过中心点(m,n)的四个不同方向上的各个电阻率值；

ρ(m+ka,n)为点(m+ka,n)处的电阻率值；

ρ(m,n+ka)为点(m,n+ka)处的电阻率值；

ρ(m+ka,n+ka)为点(m+ka,n+ka)处的电阻率值；

ρ(m+ka,n-ka)为点(m+ka,n-ka)处的电阻率值；

中值滤波的输出中值ρmlm(m,n)被定义为：

ρmlm(m,n)＝median[ymax(m,n),ymin(m,n),ρ(m,n)](13)

其中：

ymax(m,n)＝max1≤ib≤4[zib(m,n)](14)

ymin＝min1≤ib≤4[zib(m,n)](15)

zib(m,n)＝median[wib[(m,n)]],ib＝1,2,3,4(16)

随后计算cia,ja＝2n+1滤波窗口内的平均值ρave(m,n)，则：

然后将中心点(m,n)处的平均值ρave(m,n)和中心点(m,n)处的中值ρmlm(m,n)进行比值，即得到静态校正系数k(m,n)，k(m,n)数学表达式如下：

最后，将公式(18)所得到的静态校正系数k(m,n)与实际电阻率值ρ(m,n)相乘，即可得到校正后的视电阻率值ρ′(m,n)，数学表达式如下：

ρ′(m,n)＝k(m,n)×ρ(m,n)(19)

至此，完成基于多方向时变二维中值空间滤波法对csamt数据的静态校正；

iii.经过上述步骤ii得到静态校正后的csamt数据，首先对该静态校正后的csamt数据进行一维反演，然后进行二维反演以获得卡尼亚电阻率反演断面图；

进一步对反演后的卡尼亚电阻率反演结果进行阈值分割处理，去除浅部异常信息；

具体过程如下：

将电阻率剖面自上而下划分为g层，计算每层电阻的平均值

按照公式(20)计算均值阈值t1，即：

其中，g为电阻率剖面所划分的层数，为第g层电阻的平均值，g的取值为[1,g]；

设定ρ(x,y)为静态校正后的电阻率值，假设阈值分割后的电阻率值用表示，将高于阈值t1的电阻率值保持不变，其它电阻率用常数s表示；

则阈值分割公式如下式(21)所示：

利用公式(21)对反演后的卡尼亚电阻率反演结果进行阈值分割，去除浅部异常信息；

iv.对经过步骤iii阈值分割后获得的卡尼亚电阻率反演断面图，进行基于幂函数的加权滤波处理，增强深部采空区电性特征，从而刻画出清晰的深部采空区边界；

具体过程如下：

设定阈值为t2，t2为大于0的常数，将阈值分割后的csamt数据中，低于阈值t2的数值全部变为原数值的a次方，大于t2的数值保持不变，如公式(22)所示：

其中，a表示大于1的常数；利用公式(22)对阈值分割处理后的csamt数据进行处理，增强深部采空区电性特征，从而刻画出清晰的深部采空区边界。

本发明具有如下优点：

如上所述，本发明提供了一种适用于强电磁干扰地区的深部采空区探测方法，该方法通过依次对原始获得的csamt数据进行异常点的剔除、静态效应的消除、近地表低阻异常的去除、深部不合理低阻异常的压制等四个层级的处理后，使得测区深部采空区有效电性特征得到了明显增强，所刻画的采空区边界更加清晰，从而提高测区采空区的探测精度。

附图说明

图1为本发明实施例中适用于强电磁干扰地区的深部采空区探测方法的流程图。

图2为本发明实施例中csamt数据异常点剔除深度神经网络结构示意图。

图3为本发明实施例中多方向滤波窗口示意图。

图4为本发明实施例中电阻率剖面层数划分示意图。

图5为本发明实施例中卡尼亚电阻率曲线异常数据点分布图。

图6为本发明实施例中实测数据静态校正示意图。

图7为本发明实施例中某一测线静态校正前后对比示意图。

图8为本发明实施例中均值阈值分割前后对比示意图。

图9为本发明实施例中幂函数法处理后的电阻率断面图。

具体实施方式

下面结合附图以及具体实施方式对本发明作进一步详细说明：

结合图1所示，一种适用于强电磁干扰地区的深部采空区探测方法，包括如下步骤：

i.基于深度神经网络学习与预测方法对原始采集的csamt数据进行异常值剔除。

csamt在野外数据采集工作过程中，由于受到各种电磁干扰，采集到的数据存在噪音，导致个别数据出现“飞点”现象，因此必须对这些异常点进行编辑或去噪。

然而，传统的手工剔除方法具有效率低、精度低等缺点，面对海量的csamt数据，本发明实施例采用深度神经网络进行异常值剔除，以提高异常值剔除精度与效率。

深度神经网络(deepneuralnetworks,简称dnn)与传统的浅层神经网络(bp神经网络)相比具有更多的隐藏层，因此dnn网络具有强适应性、强泛化能力，具有更强大的表征能力。深度神经网络具有多种多样的网络结构，但没有统一的网络架构，经典的深度神经网络一般是由一个输入层、多个隐藏层、输出层组成，每层都包括多个神经元。因此，本发明针对csamt数据分布以及异常点特征，设计一套适用于剔除csamt异常数据点的深度神经网络网络结构，以实现异常数据点的智能剔除，提高生产效率，降低勘探成本。

具体过程如下：

i.1.设计深度神经网络结构

本发明针对csamt数据异常特征设计了如图2所示的深度神经网络结构。

其中，设计两个输入层，输入参数分别为频率、电阻率；设计两个输出层，输出为0或者1，若输出为0，则代表为异常值，进行剔除；若输出为1，代表正常值，进行保留。

深度神经网络的隐藏层数设为k层，其中，k为大于0的正整数。

由于深度神经网络的隐藏层数k一般大于5层，因此本发明将隐藏层数k设为5-10层。

隐藏层节点数的确定是一个相对复杂的过程，如果隐藏层的节点数过多会使深度网络进入漫长的学习过程，耗费大量的时间，从而导致模型的泛学习能力下降；

如果隐藏层的节点数过少会使神经网络的学习时间过短，导致模型的容错能力降低，两者均会使神经网络的预测结果不准确。

为了确定最佳的隐藏层节点数，需要对同一数据集设置不同数目的节点进行神经网络的训练，从而找到能够使神经网络的误差降到最低的隐藏层节点数目。

本发明实施例给出了每个隐藏层神经元的数量的计算公式，如公式(1)所示：

l＝2+α1(1)

式中，α1为正整数，取值范围为[1,10]，α1的值由深度神经网络的具体训练情况来设定。

其中，深度神经网络的激励函数采用relu函数，深度神经网络优化算法采用adagrad算法，深度神经网络采用dropout正则化方法以防止过拟合现象。

i.2.基于样本数据曲线拟合进行标签制作，具体过程如下：

从原始采集的csamt数据中随机选取几组点号的数据作为样本数据。

已知某一点号获得的一系列csamt数据为首先将电阻率转换为对数值，则：

其中，i1代表某一点号的第i1个数据点。

代表某一点号的第i1个卡尼亚电阻率，表示第i1个卡尼亚电阻率对应的频率。

随后，对数据点通过最小二乘法进行曲线拟合，获得函数解析式为：

最后，依据曲线拟合结果进行标签制作，则标签函数：

式中，ε为误差；若输出f＝0，则第i1个卡尼亚电阻率对应的点代表为异常值，进行剔除，若输出f＝1，则第i1个卡尼亚电阻率对应的点代表正常值，进行保留。

该步骤i.2是基于深度神经网络进行异常值剔除的关键一步，也是制作csamt数据样本集的重要环节，所制作样本的有效性，直接影响到异常值剔除的效果。

i.3.将步骤i.2中的样本数据进行划分，80％的作为学习样本，20％的作为验证样本。

首先将学习样本输入步骤i.1设计的深度神经网络结构中，进行深度神经网络模型训练。

假设深度神经网络某一神经元的输出为：

式中，j为大于0的整数；

为深度神经网络中第k个隐藏层的第j个神经元的输出，k∈k；

为深度神经网络中第k个隐藏层的第j个神经元的输入；

f为深度神经网络的激活函数；

为第k-1层的第i个神经元到第k层的第j个神经元的连接权重；

为深度神经网络第k层第j个神经元的偏置项。

在误差反向传播中，需要对深度神经网络的预测输出和期望输出进行计算，得到深度神经网络的预测输出和期望输出之间的误差，误差计算公式为：

式中，dic为期望输出，yic为深度神经网络的预测输出，q为训练样本数目。

由公式(5)和(6)可知，误差函数ek是由各个神经元决定的；不断对深度神经网络的连接权重和偏置项进行修改，使深度神经网络的误差函数ek满足设定的误差函数阈值。

连接权重和偏置的更新过程如下：

bj^k(h+1)＝bj^k-1(h)+ηδi(8)

式中，h表示第h个训练步数，h为大于0的整数；

为第h个训练步数得到的权系数；η为学习率，η的取值范围为(0,0.05]；δi＝(di-yi)y′i，其中，y′i为yi的导数；bj^k(h)为第h个训练步数得到的偏置。

通过上述步骤，完成深度神经网络模型的训练。

随后，用检验样本对深度神经网络模型的异常值剔除效果以及稳定性进行测试，若测试结果的准确率达到精度要求，则完成模型的训练与测试。

若不能达到精度要求，即测试结果的准确率达不到精度要求，则重复上述模型训练过程，直至深度神经网络模型达到精度要求为止。

i.4.对野外的不同剖面不同点号处的数据输入训练好的模型，自动进行异常点的剔除。

ii.基于多方向时变二维中值空间滤波法对异常值剔除后的csamt数据进行静态校正。

中值滤波作为一项成熟的信号处理技术，在各个领域具有广泛的应用。但传统的一维中值滤波存在滤波方向、滤波窗口大小单一等局限性，滤波方向以及滤波窗口的选择不当，既影响随机噪声较强位置处对噪声的削弱，又损伤有效信号较强位置处的有效信号。

基于以上原因，本发明提出了一种多方向自适应时变二维中值空间滤波法。首先根据csamt数据分布特征计算自适应时变滤波窗口，最大限度地保护有效信号，使噪音衰减能力与保护有效信号细节能力进一步加强；随后进行多方向滤波改善对滤波噪音的压制效果。

利用多方向时变二维中值空间滤波法进行静态校正。具体过程如下：

ii.1.确定滤波窗口

选取c作为中值滤波的基准窗口，对异常值剔除后的csamt数据进行中值滤波。

位置(m,n)处的中值输出表示为：

其中，表示以c作为基准窗口(m,n)位置处的中值。

ia＝m-(c-1)/2,…,m+(c-1)/2；ja＝n-(c-1)/2,…,n+(c-1)/2。

随后引入下式，计算阈值t：

ρix,iy表示原始电阻率值，ix＝1,…,nx，iy＝1,…,ny。nx为电阻率剖面横向的点数，ny为电阻率剖面纵向的点数。再根据阈值t，分别设定时变窗口的大小：

式中，cia,ja为滤波时的应用窗口大小，α,β,γ,δ分别表示常数，且α＞β,δ＞γ。

ii.2.确定滤波方向

对于cia,ja＝2n+1为滤波窗口，当中心点位于(m,n)时，滤波函数被定义为：

其中，n为大于0的常数，w1[m,n]、w2[m,n]、w3[m,n]、w4[m,n]分别表示经过中心点(m,n)的四个不同方向的各个电阻率值，如图3所示。其中：

w1[m,n]例如可以表示经过中心点(m,n)的90度方向上的各个电阻率值。

w2[m,n]例如可以表示经过中心点(m,n)的0度方向上的各个电阻率值。

w3[m,n]例如可以表示经过中心点(m,n)的135度方向上的各个电阻率值。

w4[m,n]例如可以表示经过中心点(m,n)的45度方向上的各个电阻率值。

ρ(m+ka,n)为点(m+ka,n)处的电阻率值。

ρ(m,n+ka)为点(m,n+ka)处的电阻率值。

ρ(m+ka,n+ka)为点(m+ka,n+ka)处的电阻率值。

ρ(m+ka,n-ka)为点(m+ka,n-ka)处的电阻率值。

中值滤波的输出中值ρmlm(m,n)被定义为：

ρmlm(m,n)＝median[ymax(m,n),ymin(m,n),ρ(m,n)](13)

其中：

ymax(m,n)＝max1≤ib≤4[zib(m,n)](14)

ymin＝min1≤ib≤4[zib(m,n)](15)

zib(m,n)＝median[wib[(m,n)]],ib＝1,2,3,4(16)

随后计算cia,ja＝2n+1滤波窗口内的平均值ρave(m,n)，则：

然后将中心点(m,n)处的平均值ρave(m,n)和中心点(m,n)处的中值ρmlm(m,n)进行比值，即得到静态校正系数k(m,n)，k(m,n)数学表达式如下：

最后，将公式(18)所得到的静态校正系数k(m,n)与实际电阻率值ρ(m,n)相乘，即可得到校正后的视电阻率值ρ′(m,n)，数学表达式如下：

ρ′(m,n)＝k(m,n)×ρ(m,n)(19)

至此，完成基于多方向时变二维中值空间滤波法对csamt数据的静态校正。

iii.经过上述步骤ii得到静态校正后的csamt数据，首先对该静态校正后的csamt数据进行一维反演，然后进行二维反演以获得卡尼亚电阻率反演断面图。

为压制浅部低阻异常干扰，突出深部采空区电性特征，需要进一步对反演后的卡尼亚电阻率反演结果进行阈值分割处理，以去除浅部异常信息。

具体过程如下：

如图4所示，将电阻率剖面自上而下划分为g层，计算每层电阻的平均值

按照公式(20)计算均值阈值t1，即：

其中，g为电阻率剖面所划分的层数，为第g层电阻的平均值，g的取值为[1,g]。

设定ρ(x,y)为静态校正后的电阻率值，假设阈值分割后的电阻率值用表示，将高于阈值t1的电阻率值保持不变，其它电阻率用常数s表示。

则阈值分割公式如下式(21)所示：

首先根据电阻率反演结果和近地表低阻体分布状况大致估算出浅部异常区的深度范围。然后利用公式(21)对浅部异常区深度范围内的卡尼亚电阻率反演结果进行阈值分割：

若电阻率值大于t1，则电阻率保持不变；若电阻率值小于t1，则低于t1的电阻率值全部变为某一常数s，可以有效压制低阻异常，去除浅部异常信息。

iv.对经过步骤iii阈值分割后获得的卡尼亚电阻率反演断面图，进行基于幂函数的加权滤波处理，增强深部采空区电性特征，从而刻画出清晰的深部采空区边界。

具体过程如下：

设定阈值为t2，t2为大于0的常数，将阈值分割后的csamt数据中，低于阈值t2的数值全部变为原数值的a次方，大于t2的数值保持不变，如公式(22)所示：

其中，a表示大于1的常数；利用公式(22)对阈值分割处理后的csamt数据进行处理，一般来说，采空区的电性异常相对比于巷道等非采空区域而言一般为低值，若a取大于1的常数，经过a次方处理后，采空区的电性异常值相对比于非采空区区域将会明显降低，因此可以有效增强深部采空区电性特征，刻画出清晰的深部采空区边界。

相比于原始电性特征图，经过本发明方法处理后，csamt深部采空区有效电性特征得到了明显增强，所刻画的采空区边界更加清晰，从而提高了测区深部采空区的探测精度。

本发明方法中各个步骤的组合以及他们之间的特定的组合顺序共同构成了本发明对现有技术做出贡献的技术特征，缺一不可，具体分析如下：

步骤i首先利用深度神经网络进行异常点剔除，保证csamt数据的稳定性及反演效果，有助于反演过程的收敛，为下文中值空间滤波法进行静态校正等提供高信噪比数据，提高地质解释的准确性；步骤ii在经过基于多方向时变二维中值空间滤波法的静态校正后，基本消除垂向条带状虚假异常，电阻率更加接近真实电性特征，为一维及二维高精度反演提供保证。步骤iii在经过静态校正后，为了降低浅层低阻异常(如由地表河流、池塘、近地表地下水汇集区域等引起的可能强低阻异常)对反演精度的影响，采用阈值分割方法以去除浅部异常信息，经过均值阈值处理后的反演断面图，浅部区域低阻异常消失，即浅部异常低阻信息得到压制，深部信息得到增强。由于深部采空区的电阻率值较小，可能会产生采空区分布范围扩散现象，步骤iv对获得的电阻率反演断面图进行基于幂函数的加权滤波处理，经过幂函数法处理后的断面图可以更清晰刻画采空区范围展布，同时深部采空区电性异常得到增强。

下面给出具体的实例，以验证本发明提出的深部采空区探测方法的有效性。

1.剔除原始数据异常点

取原始的csamt数据中的某一侧线，其卡尼亚电阻率随频率变化曲线分布如图5所示，采用深度神经网络对异常值进行剔除，如图5中的离散点为异常值。

由图5可以看出，当频率f＝4000hz时，该频点出现较明显异常数据点(“跳点”)，卡尼亚电阻率明显低于其同一频率的其他数据。

因此，通过剔除图5中突出“跳点”降低随机干扰引起的误差，以保证反演效果。

此即完成了第一级的异常点处理。

2.消除静态效应

由于地表或近地表存在电性不均匀差异或者存在高低起伏的地形，并且当不均匀的地质体相对于电磁波波长小很多时，会造成不均匀体表面电荷积聚，进而导致电场畸变，畸变后的数值主要表现在视电阻率曲线在视电阻率轴发生向上或向下平行移动，即为静态效应。

图6为本发明实施例中实测数据静态校正示意图。

其中，图6(a)、图6(b)分别为对上述测线中1345点和2120点的高阻、低阻异常进行空间滤波进行静态校正后视电阻率-频率结果示意图。

由上述结果表明静态校正后的卡尼亚电阻率曲线更符合实际情况，消除了虚假异常。

图7为某一测线静态校正前后对比图。

从图7(a)中可以看出，在静态校正前，局部“挂面条”现象严重，且浅部有杂乱圈闭异常，电阻率值跳跃较大，部分电阻率值发生畸变。如图7(a)圆圈内所示，出现陡立的等值线，呈现虚假的低阻异常，掩盖实际的地质情况。如图7(b)所示，在经过静态校正后“挂面条”现象减轻，基本消除了垂向条带状虚假异常，电阻率更加接近真实电性特征。

3.去除浅部低阻异常

在经过静态校正后，本发明采用阈值分割方法以去除浅部异常信息。

图8(a)是均值阈值分割处理前的反演断面图，可以看出椭圆区域存在多个低阻异常区域，图8(b)则是经过均值阈值处理后的反演断面图，可以看出浅部椭圆区域低阻异常消失，即浅部异常低阻信息(椭圆)得到压制，深部信息(方形区域)得到增强。

4.增强深部采空区电性特征

图9为经过幂函数法处理后的断面图，图8(b)为未经过幂函数法处理的电阻率断面图。

由图9和图8(b)对比不难看出，经过幂函数法处理后的断面图可以更清晰刻画采空区范围的边界，同时深部采空区电性异常得到增强。

本发明通过对由csamt方法在潍日高速公路测区(下伏存在隐伏煤矿采空区)所采集数据进行处理的应用，得出如下认识：

(1)本发明方法通过对原始csamt数据进行异常点的剔除、静态效应的消除、近地表低阻异常的去除、深部采空区电性特征增强等四个层级的处理后，测区深部采空区有效电性特征得到了明显增强，所刻画的采空区边界更加清晰。

以上说明本发明方法具有压制浅部异常干扰，突出深部采空区电性特征的作用。

(2)利用处理后的反演断面图，结合测区地质资料，布设钻孔并实施钻探进行验证，结果表明，实际采空区的埋深与物探异常具有很高的吻合度。

依此，进一步查明测区采空区的位置、分布范围，实现测区采空区更高精度的探测。

当然，以上说明仅仅为本发明的较佳实施例，本发明并不限于列举上述实施例，应当说明的是，任何熟悉本领域的技术人员在本说明书的教导下，所做出的所有等同替代、明显变形形式，均落在本说明书的实质范围之内，理应受到本发明的保护。