一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法

1.本发明涉及数据处理技术领域，尤其涉及一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法。


背景技术：

2.复电阻率探测技术通过布置大量的电极组向地下供电，利用测量阵列获取成百上千组测量数据集。采用层析反演技术计算并获取地下复电阻率分布剖面，进而解译分析污染区分布。在污染场地修复与管控工程中，常采用浅地表覆盖防渗膜的方法来阻隔降雨入渗，防止污染场地的二次污染扩散风险。地表覆盖防渗膜常在其施工或钻孔取样过程中遭到局部破损。对于这类污染场地的复电阻率探测中，浅地表覆盖的防渗膜区域会阻断电信号的传递，影响复电阻率测量数据质量，而防渗膜破损区域会产生电信号传递通道，可获取相对可靠的复电阻率测量数据。在浅地表覆膜污染场地如何利用复电阻率探测技术识别防渗膜破损区域并识别可靠的复电阻率测量数据，进而提高复电阻率探测技术对污染场地调查及其解译的准确性是急需解决的问题。


技术实现要素：

3.本发明克服了现有技术的不足，提供了一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法。
4.为达上述目的，本发明采用的技术方案为：本发明提供了一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，包括以下步骤：构建不同破损比的防渗膜覆盖条件下的污染区复电阻率二维模型；对所述电阻率二维模型进行正演与反演数值模拟，计算所述电阻率二维模型反演误差以及污染区反演误差；根据所述电阻率二维模型反演误差以及污染区反演误差建立不同防渗膜破损比的拟合曲线，并确定有效防渗膜破损比；提取所述电阻率二维模型的浅层电阻率数据，并计算防渗膜破损区电阻率阈值，进而确定有效防渗膜破损比阈值；提取并计算污染场地浅地表电阻率数据，根据所述防渗膜破损区电阻率阈值计算污染场地浅地表覆膜破损比；根据有效防渗膜破损比阈值筛选有效复电阻率数据区域。
5.进一步地，本发明的一个较佳实施例中，防渗膜设置深度为实际场地防渗膜的平均埋深，在防渗膜的延伸方向等间距设置相同尺寸的破损区，防渗膜破损比的计算公式为：；
其中n为单个防渗膜破损区的数量，其取值为0 ~ ne的正整数，ne为测线电极数的一半，为单个防渗膜破损区的长度，de为测线电极间距，n取1~4的正整数，l为复电阻率测线长度，为防渗膜破损比。
6.进一步地，本发明的一个较佳实施例中，所述的一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，还包括：所述电阻率二维模型的宽度为现场实际测线长度l，模型单元的宽度相同，取值范围，深度为场地测量数据的反演复电阻率剖面深度，模型单元纵向尺寸从上到下按相同的倍数逐层放大，模拟电极等间距布置，与现场实际布置一致；污染体分布区设置为矩形区域，其横向宽度设置为，其埋深范围根据现场污染区平均分布区域确定，防渗膜设置为一层单元，其单元厚度为现场铺设防渗膜的实际厚度；其中，de为测线电极间距。
7.进一步地，本发明的一个较佳实施例中，所述的一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，还包括：所述的复电阻率二维模型，其防渗膜的复电阻率参数设置包括：防渗膜电阻率值大于等于1
×
105ω*m，防渗膜极化率值小于等于1
×
10-2
mv/v；其模型中污染区的复电阻率参数设置包括：场地复电阻率测线反演的污染区的平均电阻率值和极化率值；其模型中背景复电阻率参数设置包括：场地复电阻率测线反演的背景电阻率和极化率平均值。
8.进一步地，本发明的一个较佳实施例中，在对所述电阻率二维模型进行正演与反演数值模拟，计算所述电阻率二维模型反演误差以及污染区反演误差的步骤中，具体包括：其中，数值模拟中所采用的数据与现场复电阻率测量采用的一致，模型正演模拟采用有限差分或有限元方法，反演模拟采用高斯牛顿最小二乘法迭代，其中模型演误差计算公式如下：；其中为模型电阻率反演误差，m为复电阻率模型的单元数，为无防渗膜覆盖的复电阻率模型的第k个单元反演电阻率，为反演的复电阻率模型的第k个单元的电阻率值；为模型极化率反演误差，为无防渗膜覆盖的复电阻率模型的第k个单元反演极化率值，为反演的复电阻率模型的第k个单元的极化率值；其模型的污染区反演误差计算如下：
；其中为模型污染区电阻率反演误差，c为复电阻率模型污染区的单元数，为无防渗膜覆盖的复电阻率模型污染区的第i个单元反演电阻率，为反演的复电阻率模型污染区的第i个单元的电阻率值；为模型污染区极化率反演误差，为无防渗膜覆盖的复电阻率模型污染区的第i个单元反演极化率值，为反演的复电阻率模型污染区的第i个单元的极化率值。
9.进一步地，本发明的一个较佳实施例中，在根据所述电阻率二维模型反演误差以及污染区反演误差建立不同防渗膜破损比的拟合曲线，并确定有效防渗膜破损比的步骤中，具体包括：分别以模型电阻率反演误差、模型极化率反演误差、模型污染区电阻率反演误差以及模型污染区极化率反演误差为纵轴，以防渗膜破损比为横轴，进行曲线拟合，获得模型电阻率反演误差曲线、模型极化率反演误差曲线、模型污染区电阻率反演误差曲线以及模型污染区极化率曲线；当所述模型电阻率反演误差曲线、模型极化率反演误差曲线、模型污染区电阻率反演误差曲线以及模型污染区极化率曲线在同一防渗膜破损比处的斜率变化均小于1%时，则满足该条件的防渗膜破损比的最大值为有效防渗膜破损比。
10.进一步地，本发明的一个较佳实施例中，在提取所述电阻率二维模型的浅层电阻率数据，并计算防渗膜破损区电阻率阈值，进而确定有效防渗膜破损比阈值的步骤中，具体包括：其防渗膜破损区电阻率阈值的计算公式为：；其中为沿测线方向计算水平坐标为x处，模型表面以下第j个单元的电阻率值，h为复电阻率模型浅层单元数；计算所有防渗膜破损区的模型浅层单元的平均电阻率值的均值，以此作为防渗膜破损区电阻率阈值；以平均电阻率值为纵轴，以单元水平位置坐标值x为横轴，构建浅层模型平均电阻率曲线，在曲线中，提取小于防渗膜破损区电阻率阈值的曲线部分，计算该曲线部分水平长度与模型宽度l的比值，即有效防渗膜破损比阈值。
11.进一步地，本发明的一个较佳实施例中，在提取并计算污染场地浅地表电阻率数据，根据模型防渗膜破损区电阻率阈值计算污染场地浅地表覆膜破损比的步骤中，具体包括：污染场地浅地表电阻率数据计算公式如下：
；其中为沿测线方向测线水平坐标为x处；测线区地表以下第j个参数单元的电阻率值；h为浅地表参数单元数；为沿测线方向计算其在测线同一水平位置的参数单元的平均电阻率值，亦代表污染场地浅地表电阻率数据。
12.进一步地，本发明的一个较佳实施例中，所述的一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，还包括：以浅地表平均电阻率值为纵轴，以单元水平位置坐标值x为横轴，构建场地浅地表平均电阻率曲线；在曲线中，提取小于防渗膜破损区电阻率阈值的曲线部分，定义为有效参数曲线；当相邻有效参数曲线间隔小于时，相邻有效参数曲线及其所夹曲线部分可合并为同一曲线区域，并对不同曲线区域沿测线方向依次标号，并计算不同曲线区域中有效参数曲线水平长度之和与该曲线区域水平长度的比值，即为曲线区域的污染场地浅地表覆膜破损比；其中为有效复电阻率模型的单个防渗膜破损区的长度。
13.进一步地，本发明的一个较佳实施例中，在根据有效防渗膜破损比阈值筛选有效复电阻率数据区域的步骤中，具体包括：当曲线区域的污染场地浅地表覆膜破损比大于等于有效防渗膜破损比阈值时，则定义该曲线区域水平范围的场地测量数据的反演复电阻率剖面区域为有效复电阻率数据区域。
14.本发明解决了背景技术中存在的缺陷，本发明具备以下有益效果：本发明提出系统地计算污染场地浅地表防渗膜破损区复电阻率特征参数的方法，有效地识别并量化了污染场地浅地表防渗膜破损区的分布。本发明提出浅地表覆膜污染场地可靠的复电阻率测量数据评价和筛选技术，有助于提升复电阻率探测技术在污染场地调查及其解译工作中的准确性。
附图说明
15.为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他实施例的附图。
16.图1为一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法的流程图；图2为本发明中不同破损比的防渗膜覆盖条件下的复电阻率模型及参考模型示意图；图3为本发明中不同破损比的防渗膜覆盖条件下的污染区复电阻率反演剖面；图4为本发明中模型和污染区反演误差与不同防渗膜破损比的拟合曲线；图5为本发明中有效复电阻率模型的浅层平均电阻率曲线及其计算的防渗膜破损区电阻率阈值和有效防渗膜破损比阈值；
图6为本发明中污染场地浅地表电阻率曲线及有效复电阻率数据区域。
具体实施方式
17.为了能够更清楚地理解本发明的上述目的、特征和优点，下面结合附图和具体实施方式对本发明进行进一步的详细描述。需要说明的是，在不冲突的情况下，本技术的实施例及实施例中的特征可以相互组合。
18.在下面的描述中阐述了很多具体细节以便于充分理解本发明，但是，本发明还可以采用其他不同于在此描述的其他方式来实施，因此，本发明的保护范围并不受下面公开的具体实施例的限制。
实施例1
19.为实现上述目的，本发明提供了一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，所述方法如图1所示，具体包括：s1：构建如图2所示的不同破损比（2.38%~50.79%）的防渗膜覆盖条件下的污染区复电阻率二维模型；s2：对复电阻率二维模型进行正演与反演数值模拟，计算复电阻率二维模型反演误差及污染区反演误差，如图3所示；s3：根据电阻率二维模型反演误差以及污染区反演误差建立不同防渗膜破损比的拟合曲线，如图4所示，确定有效防渗膜破损比；s4：如图5所示，提取所述电阻率二维模型的浅层电阻率数据，计算防渗膜破损区电阻率阈值，进而确定有效防渗膜破损比阈值；s5：提取并计算污染场地浅地表电阻率数据，根据所述防渗膜破损区电阻率阈值计算污染场地浅地表覆膜破损比，如图6(1)所示；s6：如图6(2)和(3)所示，根据有效防渗膜破损比阈值筛选有效复电阻率数据区域。
20.构建不同破损比的防渗膜覆盖条件下的污染区复电阻率二维模型，如图2所示，防渗膜设置深度为0.6 m，为实际场地防渗膜的平均埋深。在防渗膜的延伸方向等间距设置相同尺寸的破损区，防渗膜破损比（）的计算公式为：；其中n为单个防渗膜破损区的数量，在本实例中其取值为0 ~ 32的正整数，de为测线电极间距2 m，n取2，单个防渗膜破损区的长度值为1 m，实例中模型的破损比范围为2.38%~50.79%，见图2和图3所示。
21.复电阻率二维模型其宽度为现场实际测线长度126 m，模型单元的宽度相同为0.5 m，深度为场地测量数据的反演复电阻率剖面深度21.3 m，模型单元纵向尺寸从上到下按1.1倍逐层放大；模拟电极等间距与现场实际布置一致，均为2 m；污染体分布区设置为矩形区域，其横向宽度设置为42 m，其埋深范围为2 m ~12 m；防渗膜设置为一层单元，其单元厚度为0.5 mm。
22.复电阻率二维模型，其防渗膜的复电阻率参数设置：防渗膜电阻率值为1
×
10
5 ω*m，防渗膜极化率值为1
×
10-2
mv/v；其模型污染区的电阻率值为200 ω*m，极化率值为1 mv/v；模型背景电阻率值为10 ω*m，极化率值为40 mv/v；数值模拟的测量数组阵列与现场复电阻率测量采用的一致，均为梯度阵列；模型正演模拟采用有限差分法计算，反演模拟采用高斯牛顿最小二乘法迭代；模型反演误差计算公式如下：；其中为模型电阻率反演误差，m为复电阻率模型的单元数，为无防渗膜覆盖的复电阻率模型的第k个单元反演电阻率，为反演的复电阻率模型的第k个单元的电阻率值；为模型极化率反演误差，为无防渗膜覆盖的复电阻率模型的第k个单元反演极化率值，为反演的复电阻率模型的第k个单元的极化率值。
23.模型的污染区反演误差计算如下：；其中为模型污染区电阻率反演误差，c为复电阻率模型污染区的单元数，为无防渗膜覆盖的复电阻率模型污染区的第i个单元反演电阻率，为反演的复电阻率模型污染区的第i个单元的电阻率值；为模型污染区极化率反演误差，为无防渗膜覆盖的复电阻率模型污染区的第i个单元反演极化率值，为反演的复电阻率模型污染区的第i个单元的极化率值。在本实例中，反演结果的见图3所示。
24.分别以为纵轴，以为横轴，进行曲线拟合，获得四条拟合曲线分别进行曲线拟合，获得四条拟合曲线分别为得四条拟合曲线分别进行曲线拟合，获得四条拟合曲线分别为，如图4所示；当四条拟合曲线在同一防渗膜破损比处的斜率变化均小于1%时，则满足该条件的防渗膜破损比的最大值为有效防渗膜破损比，如图4所示，有效防渗膜破损比为0.381。
25.有效防渗膜破损比0.381所对应的复电阻率模型为图5(1)所示模型为有效复电阻率模型，其防渗膜破损比为0.254；提取浅层0.6 m以上的单元的电阻率数据，沿测线方向计算同一水平位置的单元的平均电阻率值，其计算公式如下：
；其中为沿测线方向计算水平坐标为x处，模型表面以下第j个单元的电阻率值，h为复电阻率模型浅层单元数。以平均电阻率值为纵轴，以单元水平位置坐标值x为横轴，构建浅层模型平均电阻率曲线，如图5(2)所示。计算所有防渗膜破损区的模型浅层单元的平均电阻率值的均值，以此作为防渗膜破损区电阻率阈值，其值为150 ohm
·
m。在曲线中，提取小于防渗膜破损区电阻率阈值的曲线部分，计算该曲线部分水平长度与模型宽度l的比值，即有效防渗膜破损比阈值，其值为51.59%。
26.计算场地测量数据的反演电阻率剖面中2 m以上参数单元的电阻率数据，并沿测线方向计算其在测线同一水平位置的参数单元的平均电阻率值。其计算公式如下：；其中为沿测线方向测线水平坐标为x处，测线区地表以下第j个参数单元的电阻率值，h为浅地表参数单元数。以浅地表平均电阻率值为纵轴，以单元水平位置坐标值x为横轴，构建场地浅地表平均电阻率曲线，如图6(1)所示。在曲线中，提取小于防渗膜破损区电阻率阈值150 ohm
·
m的曲线部分，定义为有效参数曲线；当相邻有效参数曲线间隔小于3 m时，相邻有效参数曲线及其所夹曲线部分可合并为同一曲线区域，并对不同曲线区域沿测线方向依次标号（1，2
……
z）；计算不同曲线区域中有效参数曲线水平长度之和与该曲线区域水平长度的比值，即为曲线区域的防渗膜破损比，如图6(1)所示，的值为62%，其大于有效防渗膜破损比阈值51.59%，则定义该曲线区域水平范围的场地测量数据的反演复电阻率剖面区域为有效复电阻率数据区域，如图6(2)和(3)中黑色曲线区域所示。
27.综上所述，本发明第一方面提出系统地计算污染场地浅地表防渗膜破损区复电阻率特征参数的方法，有效地识别并量化了污染场地浅地表防渗膜破损区的分布。本发明第二方面提出浅地表覆膜污染场地可靠的复电阻率测量数据评价和筛选技术，有助于提升复电阻率探测技术在污染场地调查及其解译工作中的准确性。
28.在本技术所提供的几个实施例中，应该理解到，所揭露的设备和方法，可以通过其它的方式实现。以上所描述的设备实施例仅仅是示意性的，例如，所述单元的划分，仅仅为一种逻辑功能划分，实际实现时可以有另外的划分方式，如：多个单元或组件可以结合，或可以集成到另一个系统，或一些特征可以忽略，或不执行。另外，所显示或讨论的各组成部分相互之间的耦合、或直接耦合、或通信连接可以是通过一些接口，设备或单元的间接耦合或通信连接，可以是电性的、机械的或其它形式的。
29.上述作为分离部件说明的单元可以是、或也可以不是物理上分开的，作为单元显示的部件可以是、或也可以不是物理单元；既可以位于一个地方，也可以分布到多个网络单元上；可以根据实际的需要选择其中的部分或全部单元来实现本实施例方案的目的。
30.另外，在本发明各实施例中的各功能单元可以全部集成在一个处理单元中，也可以是各单元分别单独作为一个单元，也可以两个或两个以上单元集成在一个单元中；上述
集成的单元既可以采用硬件的形式实现，也可以采用硬件加软件功能单元的形式实现。
31.本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：移动存储设备、只读存储器（rom，read-only memory）、随机存取存储器（ram，random access memory）、磁碟或者光盘等各种可以存储程序代码的介质。
32.或者，本发明上述集成的单元如果以软件功能模块的形式实现并作为独立的产品销售或使用时，也可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明实施例的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令用以使得一台计算机设备（可以是个人计算机、服务器、或者网络设备等）执行本发明各个实施例方法的全部或部分。而前述的存储介质包括：移动存储设备、rom、ram、磁碟或者光盘等各种可以存储程序代码的介质。
33.以上，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何熟悉本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应以权利要求的保护范围为准。

技术特征：
1.一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，其特征在于，包括以下步骤：构建不同破损比的防渗膜覆盖条件下的污染区复电阻率二维模型；对所述电阻率二维模型进行正演与反演数值模拟，计算所述电阻率二维模型反演误差以及污染区反演误差；根据所述电阻率二维模型反演误差以及污染区反演误差建立不同防渗膜破损比的拟合曲线，并确定有效防渗膜破损比；提取所述电阻率二维模型的浅层电阻率数据，并计算防渗膜破损区电阻率阈值，进而确定有效防渗膜破损比阈值；提取并计算污染场地浅地表电阻率数据，根据所述防渗膜破损区电阻率阈值计算污染场地浅地表覆膜破损比；根据有效防渗膜破损比阈值筛选有效复电阻率数据区域。2.根据权利要求1所述的一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，其特征在于，还包括：防渗膜设置深度为实际场地防渗膜的平均埋深，在防渗膜的延伸方向等间距设置相同尺寸的破损区，防渗膜破损比的计算公式为：；其中n为单个防渗膜破损区的数量，其取值为0 ~ n
e
的正整数，n
e
为测线电极数的一半，为单个防渗膜破损区的长度，d
e
为测线电极间距，n取1~4的正整数，l为复电阻率测线长度，为防渗膜破损比。3.根据权利要求1所述的一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，其特征在于，还包括;所述电阻率二维模型的宽度为现场实际测线长度l，模型单元的宽度相同，取值范围，深度为场地测量数据的反演复电阻率剖面深度，模型单元纵向尺寸从上到下按相同的倍数逐层放大，模拟电极等间距布置，与现场实际布置一致；污染体分布区设置为矩形区域，其横向宽度设置为，其埋深范围根据现场污染区平均分布区域确定，防渗膜设置为一层单元，其单元厚度为现场铺设防渗膜的实际厚度；其中，d
e
为测线电极间距。4.根据权利要求1所述的一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，其特征在于，所述的复电阻率二维模型，其防渗膜的复电阻率参数设置包括：防渗膜电阻率值大于等于1
×
105ω*m，防渗膜极化率值小于等于1
×
10-2 mv/v；其模型中污染区的复电阻率参数设置包括：场地复电阻率测线反演的污染区的平均电阻率值和极化率值；其模型中背景复电阻率参数设置包括：场地复电阻率测线反演的背景电阻率和极化率平均值。5.根据权利要求1所述的一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，
其特征在于，在对所述电阻率二维模型进行正演与反演数值模拟，计算所述电阻率二维模型反演误差以及污染区反演误差的步骤中，具体包括：其中，数值模拟中所采用的数据与现场复电阻率测量采用的一致，模型正演模拟采用有限差分或有限元方法，反演模拟采用高斯牛顿最小二乘法迭代，其中模型演误差计算公式如下：；其中为模型电阻率反演误差，m为复电阻率模型的单元数，为无防渗膜覆盖的复电阻率模型的第k个单元反演电阻率，为反演的复电阻率模型的第k个单元的电阻率值；为模型极化率反演误差，为无防渗膜覆盖的复电阻率模型的第k个单元反演极化率值，为反演的复电阻率模型的第k个单元的极化率值；其模型的污染区反演误差计算如下：；其中为模型污染区电阻率反演误差，c为复电阻率模型污染区的单元数，为无防渗膜覆盖的复电阻率模型污染区的第i个单元反演电阻率，为反演的复电阻率模型污染区的第i个单元的电阻率值；为模型污染区极化率反演误差，为无防渗膜覆盖的复电阻率模型污染区的第i个单元反演极化率值，为反演的复电阻率模型污染区的第i个单元的极化率值。6.根据权利要求1所述的一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，其特征在于，在根据所述电阻率二维模型反演误差以及污染区反演误差建立不同防渗膜破损比的拟合曲线，并确定有效防渗膜破损比的步骤中，具体包括：分别以模型电阻率反演误差、模型极化率反演误差、模型污染区电阻率反演误差以及模型污染区极化率反演误差为纵轴，以防渗膜破损比为横轴，进行曲线拟合，获得模型电阻率反演误差曲线、模型极化率反演误差曲线、模型污染区电阻率反演误差曲线以及模型污染区极化率曲线；当所述模型电阻率反演误差曲线、模型极化率反演误差曲线、模型污染区电阻率反演误差曲线以及模型污染区极化率曲线在同一防渗膜破损比处的斜率变化均小于1%时，则满足该条件的防渗膜破损比的最大值为有效防渗膜破损比。7.根据权利要求1所述的一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，其特征在于，在提取所述电阻率二维模型的浅层电阻率数据，并计算防渗膜破损区电阻率
阈值，进而确定有效防渗膜破损比阈值的步骤中，具体包括：其防渗膜破损区电阻率阈值的计算公式为：；其中为沿测线方向计算水平坐标为x处，模型表面以下第j个单元的电阻率值，h为复电阻率模型浅层单元数；计算所有防渗膜破损区的模型浅层单元的平均电阻率值的均值，以此作为防渗膜破损区电阻率阈值；以平均电阻率值为纵轴，以单元水平位置坐标值x为横轴，构建浅层模型平均电阻率曲线，在曲线中，提取小于防渗膜破损区电阻率阈值的曲线部分，计算该曲线部分水平长度与模型宽度l的比值，即有效防渗膜破损比阈值。8.根据权利要求1所述的一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，其特征在于，在提取并计算污染场地浅地表电阻率数据，根据模型防渗膜破损区电阻率阈值计算污染场地浅地表覆膜破损比的步骤中，具体包括：污染场地浅地表电阻率数据计算公式如下：；其中为沿测线方向测线水平坐标为x处；测线区地表以下第j个参数单元的电阻率值；h为浅地表参数单元数；为沿测线方向计算其在测线同一水平位置的参数单元的平均电阻率值，亦代表污染场地浅地表电阻率数据。9.根据权利要求8所述的一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，其特征在于，还包括：以浅地表平均电阻率值为纵轴，以单元水平位置坐标值x为横轴，构建场地浅地表平均电阻率曲线；在曲线中，提取小于防渗膜破损区电阻率阈值的曲线部分，定义为有效参数曲线；当相邻有效参数曲线间隔小于时，相邻有效参数曲线及其所夹曲线部分可合并为同一曲线区域，并对不同曲线区域沿测线方向依次标号，并计算不同曲线区域中有效参数曲线水平长度之和与该曲线区域水平长度的比值，即为曲线区域的污染场地浅地表覆膜破损比；其中为有效复电阻率模型的单个防渗膜破损区的长度。10.根据权利要求1-9任一项所述的一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，在根据有效防渗膜破损比阈值筛选有效复电阻率数据区域的步骤中，具体包括：当曲线区域的污染场地浅地表覆膜破损比大于等于有效防渗膜破损比阈值时，则定义该曲线区域水平范围的场地测量数据的反演复电阻率剖面区域为有效复电阻率数据区域。

技术总结
本发明涉及一种浅地表覆膜污染场地有效复电阻率数据区域筛选方法，属于数据处理技术领域，本发明通过提取所述电阻率二维模型的浅层电阻率数据，并计算防渗膜破损区电阻率阈值，进而确定有效防渗膜破损比阈值，提取并计算污染场地浅地表电阻率数据，根据所述防渗膜破损区电阻率阈值计算污染场地浅地表覆膜破损比，最后根据有效防渗膜破损比阈值筛选有效复电阻率数据区域。本发明提出系统地计算污染场地浅地表防渗膜破损区复电阻率特征参数的方法，有效地识别并量化了污染场地浅地表防渗膜破损区的分布。本发明提出浅地表覆膜污染场地可靠的复电阻率测量数据评价和筛选技术，有助于提升复电阻率探测技术在污染场地调查及其解译工作中的准确性。其解译工作中的准确性。其解译工作中的准确性。


技术研发人员：王文峰 毛德强 马新民 李书鹏 张家铭 郭丽莉 邱景琮
受保护的技术使用者：山东大学
技术研发日：2023.07.19
技术公布日：2023/8/16