一种单个电阻率异常体的频率域传导电法的工作参数优化方法与流程

本发明涉及一种勘查地球物理领域的传导电法勘探工作参数的优化新方法。

背景技术：

目前，频率域传导电法勘探方法的工作参数一般是勘探工作人员基于各自的经验和自己对勘探区情况的大体认识进行大致估计选择。由于每个勘探工作人员的经验、对勘探区的大体情况预估等方面存在差异，从而勘探中的工作参数选择随意性较强，缺乏一个统一的衡量标尺，容易造成设计的勘探工作量过多或过少。勘探工作量过多则容易造成工作量浪费、勘探成本高、效率低等问题；勘探工作量过少则容易造成勘探效果难以满足勘探目标要求、丢失勘探目标体等问题。如何比较有效地解决具体勘探工作的工作参数设计问题，并作为实际勘探工作的具体指导，变得十分必要。

技术实现要素：
：

本发明是基于目前常规的频率域传导电法勘探方法的工作参数选择中所存在的问题，提出通过分析勘探区的地形等级、频率域传导电法的装置类型、发送机同时供电的通道数量、设计的勘探深度、欲勘探电阻率异常体的尺寸要求(最小水平长度、最小水平宽度、最小铅直厚度)、欲勘探电阻率异常体与围岩的电阻率比值、勘探时间要求系数等方面的信息，通过特定的公式获取频率域传导电法优化后的工作参数，如勘探区范围(勘探区的水平长度和水平宽度)、测点水平点距、测线水平线距、供电点水平间距、供电点数量、接收机的最小总道数等，从而实现以更合理的工作参数和勘探工作量完成勘探区的勘探工作。该方法能获得频率域传导电法更优的工作参数，从而以更合理的时间和经济成本获得可靠的勘探效果，即实现勘探成本与勘探效果的最优平衡。

一种单个电阻率异常体的频率域传导电法的工作参数优化方法，其具体步骤如下：

a)收集勘探区内的地形地貌地物情况，确定勘探区的地形等级(terrain_level)，根据公式(1)确定地形影响参数k，k＝terrain-level(1)；确定勘探区的地形等级可以采用相应的勘探标准或规范进行确定，如优选地根据中国地质调查局于2009年10月颁布的《地质调查项目预算标准》中的关于物化探的地形等级划分标准进行确定。

b)收集勘探区设计的勘探深度(prospecting_depth，单位为米，prospecting_depth>0)；即确定需要了解勘探区内多少米勘探深度范围内的电阻率异常体的分布情况。优选地，勘探区设计的勘探深度(prospecting_depth)小于等于1000米。

c)收集勘探区欲勘探电阻率异常体的最小水平长度、最小水平宽度、最小铅直厚度的尺寸要求(欲勘探电阻率异常体的最小水平长度为body_length，单位为米，body_length>0；欲勘探电阻率异常体的最小水平宽度为body_width，单位为米，body_width>0；欲勘探电阻率异常体的最小铅直厚度为body_thick，单位为米，body_thick>0；)；欲勘探电阻率异常体是指想在勘探区内寻找最小多大规模的电阻率异常体；由于勘探区内实际存在的电阻率异常体的形态一般为不规则体，且属于未知信息，故假定需要寻找的电阻率异常体为规则的长方体，并对其长宽厚的最小尺寸进行限定。优选地，欲勘探电阻率异常体的最小水平长度、最小水平宽度、最小铅直厚度尺寸均大于等于1米。

d)收集勘探区欲勘探电阻率异常体的电阻率与围岩的电阻率，并计算两者电阻率比值(resistivity_ratio，无量纲)，resistivity_ratio等于欲勘探电阻率异常体的电阻率与围岩的电阻率之比，resistivity_ratio不等于1，resistivity_ratio>0；优选地，欲勘探电阻率异常体与围岩的电阻率参数通过测试相应标本的电阻率确定；若欲勘探电阻率异常体或围岩的电阻率存在多个测试结果，则求取欲勘探电阻率异常体的电阻率平均值或围岩的电阻率平均值；即若欲勘探电阻率异常体的电阻率存在多个测试结果，则求取欲勘探电阻率异常体的电阻率平均值；若围岩的电阻率存在多个测试结果，则求取围岩的电阻率平均值。若勘探区内暂时无法获取到相应电阻率异常体和围岩的标本，则结合勘探区前期的资料，并综合其他相邻地区的相应岩性标本的电阻率参数，对本区内的电阻率异常体和围岩的电阻率参数进行估算，并把估算值作为相应的电阻率数值。优选地，欲勘探目标体与围岩的电阻率比值满足条件：resistivity_ratio大于等于2或resistivity_ratio小于等于0.5。即欲勘探电阻率异常体与围岩之间存在比较明显的电阻率差异，以便通过电法勘探工作能比较明显地区分电阻率异常体和围岩。

e)收集勘探区的勘探时间要求系数；勘探时间要求系数是指完成本项勘探工作的整个野外数据采集时间要求(time，单位为天，time>0)；优选地，勘探时间要求系数大于等于0.1天。由于有些勘探区(如城市)对于勘探效率和勘探时间有比较严格的要求，希望能尽快完成相关的勘探工作，以便尽快对电阻率异常体进行评估，并尽量减少勘探工作对勘探区内的人文活动带来的干扰。

f)收集选用的频率域传导电法仪器的发送机同时供电的通道数量(number_transmit，单位为个，number_transmit>0)；发送机同时供电的通道数量是指频率域传导电法仪器的发送机能同时供电的通道数量，且每个通道所供出信号的频率相互不同；相应的频率域传导电法仪器的接收机能同时接收发送机同时供出的所有频率信号；优选地，频率域传导电法仪器的发送机能同时发送1个或多个单频信号；频率域传导电法仪器的接收机能同时接收发送机同时发出的所有信号；发送机能发送的单频数量与接收机能接收的单频数量相同；发送机能发送的频率值与接收机能接收的频率值相同。即若发送机能同时发送n个单频信号，则接收机也能同时接收n个单频信号。即本类仪器能同时发送或接收1个或多个电信号，从而能提高勘探效率、降低勘探成本，也能满足特殊勘探区的时间等方面的要求。

g)确定本次勘探工作中选用的电法勘探装置类型并赋值(array_type)，其中二极装置的array_type＝3，三极装置的array_type＝2，四极装置的array_type＝1；根据勘探区的具体要求，选用电法勘探装置类型，并根据装置类型进行赋值。

h)采用公式(2)获取频率域传导电法优化后的工作参数，包括勘探区的水平长度和水平宽度(勘探区的水平长度为exploration_length，单位为米，exploration_length>0；勘探区的水平宽度为exploration_width，单位为米，exploration_width>0)、测点水平点距(station_space，单位为米，station_space>0)、测线水平线距(line_space，单位为米，line_space>0)、供电点水平间距(supply_space，单位为米，supply_space>0)、供电点数量(supply_number，单位为个，supply_number>0)、选用的频率域传导电法接收机的最小总道数(number_measure，单位为道，number_measure>0)；

其中max()为求取数值中最大值的函数；roundup()为向上取整函数；log10()表示以10为底的对数；||表示绝对值函数；

i)根据公式(2)获取频率域传导电法优化后的工作参数，并根据获得的工作参数开展勘探区的频率域传导电法勘探工作，获取频率域传导电法勘探数据，并进行数据处理和解译，获得勘探成果。采用优化后的工作参数，能满足勘探区对本次勘探工作的各项特殊要求，从而能实现以更低的时间、经济等成本获得可靠的勘探成果。

附图说明：

图1为本发明的一种单个电阻率异常体的频率域传导电法的工作参数优化方法流程图；

具体实施方式：

以下参照图1，结合具体实施方式对本发明做进一步说明。

a)收集勘探区内的地形地貌地物情况，并根据中国地质调查局于2009年10月颁布的《地质调查项目预算标准》中的关于物化探的地形等级划分标准进行确定，从而确定勘探区的地形等级(terrain_level)，根据公式(1)确定地形影响参数k，k＝terrain-level(1)；假定根据2009年10月颁布的《地质调查项目预算标准》中的关于物化探的地形等级划分标准确定某勘探区的地形等级为3级，即terrain_level＝3，则k＝3；

b)收集勘探区设计的勘探深度(prospecting_depth，单位为米)，假定步骤a)选定的某勘探区的勘探深度限定在20米，则prospecting_depth＝20米；

c)收集勘探区欲勘探电阻率异常体的最小水平长度、最小水平宽度、最小铅直厚度的尺寸要求(欲勘探电阻率异常体的最小水平长度为body_length，单位为米，body_length>0；欲勘探电阻率异常体的最小水平宽度为body_width，单位为米，body_width>0；欲勘探电阻率异常体的最小铅直厚度为body_thick，单位为米，body_thick>0；)；假定步骤a)选定的某勘探区想了解勘探区内存在的水平长度至少20米、水平宽度至少10米，铅直厚度至少2米的电阻率异常体，则body_length＝20米，body_width＝10米，body_thick＝2米；

d)收集勘探区欲勘探电阻率异常体与围岩的电阻率比值(resistivity_ratio)，resistivity_ratio等于欲勘探电阻率异常体的电阻率与围岩的电阻率之比，resistivity_ratio不等于1，无量纲，resistivity_ratio>0；假定步骤a)选定的某勘探区欲勘探的电阻率异常体的电阻率为50ω·m，围岩的电阻率为500ω·m，两者电阻率比值为0.1，即本勘探区想寻找相对于围岩而言的低阻异常体，则resistivity_ratio＝50/500＝0.1；

e)收集勘探区的勘探时间要求系数；勘探时间要求系数是指完成本项勘探工作的整个野外数据采集时间要求(time，单位为天，time>0)；假定步骤a)选定的某勘探区的频率域传导电法勘探工作希望在1天内完成，则time＝1天；

f)收集选用的频率域传导电法仪器的发送机同时供电的通道数量(number_transmit，单位为个)；发送机同时供电的通道数量是指频率域传导电法仪器的发送机能同时供电的通道数量，且每个通道所供出信号的频率相互不同；相应的频率域传导电法仪器的接收机能同时接收发送机同时供出的所有频率信号；优选地，频率域传导电法仪器的发送机能同时发送1个或多个单频信号；频率域传导电法仪器的接收机能同时接收发送机同时发出的所有信号；发送机能发送的单频数量与接收机能接收的单频数量相同；发送机能发送的频率值与接收机能接收的频率值相同。优选地，选用的频率域传导电法仪器的发送机同时供电的通道数量大于等于1个，即(number_transmit≥1)；假定步骤a)选定的某勘探区的勘探工作选用的频率域传导电法仪器的发送机同时供电的通道数量为5个，则number_transmit＝5；频率域传导电法仪器的接收机也能同时接收5个与发送机相应的信号。

g)确定本次勘探工作中选用的电法勘探装置类型并赋值(array_type)，其中二极装置的array_type＝3，三极装置的array_type＝2，四极装置的array_type＝1；假定步骤a)选定的某勘探区的频率域传导电法勘探工作装置类型为三极装置，则array_type＝2；

h)根据以上步骤的相关数据，采用公式(2)获取频率域传导电法优化后的工作参数，包括勘探区的水平长度和水平宽度(勘探区的水平长度为exploration_length，单位为米，exploration_length>0；勘探区的水平宽度为exploration_width，单位为米，exploration_width>0)、测点水平点距(station_space，单位为米，station_space>0)、测线水平线距(line_space，单位为米，line_space>0)、供电点水平间距(supply_space，单位为米，supply_space>0)、供电点数量(supply_number，单位为个，supply_number>0)、选用的频率域传导电法接收机的最小总道数(number_measure，单位为道，number_measure>0)；

把以上步骤的相关数据代入公式(2)，从而获得如下的工作参数：

即本勘探区如果想对20米埋深范围内的水平长度不少于20米，水平宽度不少于10米，水平厚度不少于2米的低阻异常体(resistivity_ratio＝0.1)进行控制，则需要布置60米*60米的勘探范围，测点水平点距4米，测线水平线距7米，供电点水平间距6米，供电点数量86个，且发送机同时供电的通道数量为5个的频率域传导电法勘探工作，频率域传导电法接收机的最小总道数为67道，从而能实现在1天内完成本次勘探工作；

i)根据以上获取频率域传导电法优化后的工作参数，开展勘探区的频率域传导电法勘探工作，获取频率域传导电法勘探数据，并进行数据处理和解译，获得勘探成果。采用优化后的工作参数，能满足勘探区对本次勘探工作的各项特殊要求，从而能实现以更低的时间、经济等成本获得可靠的勘探成果。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明内。