基于深度分辨率的隧道电阻率超前探测优化方法及系统与流程

[0001]
本公开属于隧道电阻率超前探测观测技术领域，具体涉及一种基于深度分辨率的隧道电阻率超前探测优化方法及系统。


背景技术：

[0002]
本部分的陈述仅仅是提供了与本公开相关的背景技术信息，不必然构成在先技术。
[0003]
随着隧道建设重点逐步向地形地质复杂的西部山区和水域阻隔的东部海峡地区转移，艰险山区深埋隧道和跨江越海隧道大量涌现，对于各类具有洞线长、埋深大等特点的深长隧道，对其进行地质调查更加困难。对于深埋长大隧道，由于各类不良地质赋存更隐蔽、致灾性更高，要求分辨率更高、目标更难探，对于隧道超前探测精度提出了更高的要求和挑战。目前传统的地表勘查条件和技术水平难以满足工程对探测深度和精细化探查的需求，而当前主流的隧道电阻率超前探测方法对含水体成像边界识别模糊、探不到、探不准，且对于一些亚米级甚至分米级的含导水构造难以有效识别，因此，需要一种新的隧道电阻率超前探测装置，能够对亚米级含水构造进行精细成像，即“孔中电极供电、掌子面阵列式测量”的隧道钻孔电阻率精细超前探测观测装置。
[0004]
与隧道常规观测模式相比，隧道钻孔电阻率超前探测观测装置在精细化成像方面更具优势：1)具有更远的测深；2)可以获取钻孔周边异常体更丰富更直接的信息。因此，隧道钻孔电阻率超前探测在隧道电阻率超前探测中具有广阔的应用前景。
[0005]
然而，据发明人了解，目前隧道钻孔电阻率超前探测基本观测模式仍存在以下关键问题尚未解决：
[0006]
为了保证反演成像的效果，需要通过布设大量的电极来获得大量的观测数据，观测工作需要耗费大量的时间和人力，基本观测模式得到的大量数据对反演计算造成了极大的挑战，降低了观测效率。
[0007]
同时，当前常规优化方法是提升整体的模型分辨率，是均一的提升，不能解决深部模型分辨率过分损失的问题。


技术实现要素：

[0008]
本公开为了解决上述问题，提出了一种基于深度分辨率的隧道电阻率超前探测优化方法及系统，本公开对观测系统进行优化，优选出对模型分辨率贡献较大的测量电极点，在保证模型分辨率的基础上最终得到优化后的掌子面电极排列，精简了电极数量，提高了探测效率。
[0009]
根据一些实施例，本公开采用如下技术方案：
[0010]
一种基于深度分辨率的隧道电阻率超前探测优化方法，包括以下步骤：
[0011]
(1)将采集的所有电极测量数据整合为综合集；
[0012]
(2)计算综合集的模型分辨率矩阵，从综合集中选择若干电极测量数据形成初始
集；
[0013]
(3)将不在初始集中电极的测量数据添加到初始集中，形成多个临时子集；
[0014]
(4)根据模型分辨率矩阵，计算每个临时子集的深度分辨率平衡优度函数，选择深度分辨率平衡优度函数值最优的临时子集作为新的初始集，判断此时的初始集的平均相对模型分辨率是否满足优化要求；若不满足要求，则返回步骤(3)，否则输出此时的初始集；
[0015]
(5)根据更新后的初始集确定掌子面测量电极的数目和位置，得到优化后的钻孔电阻率超前探测的有效观测模式。
[0016]
作为可选择的实施方式，所述步骤(1)中，在隧道掌子面上布置钻孔，钻孔内布置电极，确定钻孔长度和电极间距；依据探测精度要求，确定正反演网格大小和掌子面上的电极数量和排列。
[0017]
作为可选择的实施方式，所述步骤(1)中，使用全空间孔隧电阻率方法进行数据采集。
[0018]
作为可选择的实施方式，所述步骤(2)中，选取掌子面上下两个测量电极所能采集的全部电位数据作为初始集。
[0019]
作为可选择的实施方式，所述步骤(4)中，深度分辨率矩阵由分辨率平衡因子组成，且各分辨率平衡因子根据反演深度确定。
[0020]
作为可选择的实施方式，所述步骤(4)中，每个临时子集的深度分辨率平衡优度函数为：
[0021][0022]
作为可选择的实施方式，所述步骤(4)中，当前初始集的相对模型分辨率为此时的初始集模型分辨率矩阵的主对角线元素与综合集模型分辨率矩阵的主对角线元素的比值。
[0023]
作为可选择的实施方式，所述步骤(4)中，平均相对模型分辨率为对相对模型分辨率中的元素求平均得到的。
[0024]
作为可选择的实施方式，所述步骤(4)中，判断此时的初始集的平均相对模型分辨率是否满足优化要求的具体过程是，判断此时的初始集的平均相对模型分辨率的值是否大于设定值。
[0025]
一种基于深度分辨率的隧道电阻率超前探测优化系统，包括：
[0026]
用于将采集的所有电极测量数据整合为综合集的模块；
[0027]
用于计算综合集的模型分辨率矩阵，从综合集中选择若干电极测量数据形成初始集的模块；
[0028]
用于将不在初始集中电极的测量数据添加到初始集中，形成多个临时子集的模块；
[0029]
用于根据模型分辨率矩阵，计算每个临时子集的深度分辨率平衡优度函数，选择深度分辨率平衡优度函数值最优的临时子集作为新的初始集，判断此时的初始集的平均相对模型分辨率是否满足优化要求；若不满足要求，则返回重新添加另外电极的测量数据形成新的临时子集，否则输出此时的初始集的模块；
[0030]
用于根据更新后的初始集确定掌子面测量电极的数目和位置的模块。
[0031]
一种隧道电阻率超前探测观测系统，包括多个测量电极，所述测量电极的数目和位置根据所述一种基于深度分辨率的隧道电阻率超前探测优化方法确定。
[0032]
与现有技术相比，本公开的有益效果为：
[0033]
1、本公开一种隧道电阻率超前探测装置形式优化方法，本公开在孔中电极等间距均匀布置、掌子面阵列电极测量的基本型式的基础上，使用深度分辨率平衡优化的drb方法对观测装置进行优化，优选出对模型分辨率贡献较大的测量电极点，在保证模型分辨率的基础上最终得到优化后的掌子面电极排列，精简了电极数量，提高了探测效率。
[0034]
2、本公开在深度分辨率平衡优度函数中，引入了深度分辨率平衡矩阵h，可以结合先验信息，调整不同网格模型分辨率在drb计算过程中的分辨率平衡度因子，可以通过适当增大某一深度区域的分辨率平衡度，使得这一区域网格的模型分辨率在drb的计算中更为重要，最终使得算法能够优先提升这一深度区域的模型分辨率。
附图说明
[0035]
构成本公开的一部分的说明书附图用来提供对本公开的进一步理解，本公开的示意性实施例及其说明用于解释本公开，并不构成对本公开的不当限定。
[0036]
图1是基于隧道钻孔电阻率超前探测装置深度分辨率平衡优化方法流程图；
[0037]
图2是钻孔电阻率观测装置基本模式电极分布示意图；
[0038]
图3是平均相对模型分辨率随迭代次数的变化曲线；
[0039]
图4是优化后观测装置的掌子面测量电极分布示意图；
[0040]
图5是深度分辨率平衡因子分布示意图；
[0041]
图6是优化后的相对模型分辨率图；
[0042]
图7是进行数值模拟时使用的地电模型图；
[0043]
图8是基于优化后的隧道钻孔观测模型装置的反演结果图。
具体实施方式：
[0044]
下面结合附图与实施例对本公开作进一步说明。
[0045]
应该指出，以下详细说明都是例示性的，旨在对本公开提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本公开所属技术领域的普通技术人员通常理解的相同含义。
[0046]
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本公开的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
[0047]
在本公开中，术语如“上”、“下”、“左”、“右”、“前”、“后”、“竖直”、“水平”、“侧”、“底”等指示的方位或位置关系为基于附图所示的方位或位置关系，只是为了便于叙述本公开各部件或元件结构关系而确定的关系词，并非特指本公开中任一部件或元件，不能理解为对本公开的限制。
[0048]
本公开中，术语如“固接”、“相连”、“连接”等应做广义理解，表示可以是固定连接，也可以是一体地连接或可拆卸连接；可以是直接相连，也可以通过中间媒介间接相连。对于
本领域的相关科研或技术人员，可以根据具体情况确定上述术语在本公开中的具体含义，不能理解为对本公开的限制。
[0049]
一种隧道电阻率超前探测装置深度分辨率平衡优化方法，如图1所示，包括以下步骤：
[0050]
(1)在隧道掌子面上布置钻孔，钻孔内布置电极，确定钻孔长度和电极间距；依据探测精度要求，确定正反演网格大小和掌子面上的电极数量和排列；
[0051]
在具体实例中，当然，需要说明的是，本实施例仅为示例，并不仅代表其他实施例的钻孔深度、电极数量、间距和数据总数等需要和本实施例一致，上述数据都是可以根据具体观测要求和环境进行合理变更的，此为本领域技术人员容易想到的，理应属于本公开的保护范围。
[0052]
钻孔最大深度设定为60m，孔内供电电极以等间距的方式排布，电极间间距为2m，则钻孔内共有30个供电电极。掌子面测量电极则呈阵列式排布，如图2所示共设置有8排电极，其中第1与第8排各自仅有3个电极，第4与第5排各自有9个电极，其余各排各自拥有7个电极。
[0053]
(2)使用全空间孔隧电阻率方法进行数据采集；
[0054]
具体地，数据采集采取钻孔内供电、掌子面接收的二极法型式，共采集到1560个测量数据。
[0055]
(3)对掌子面测量电极进行优化和精简：
①
计算综合集的模型分辨率矩阵m
c
，选取某个子集作为优化的初始集s
i
；
②
将每个其他的观测装置与初始集结合形成众多临时子集s
t
，计算出相应的模型分辨率矩阵m
t
；
③
确定深度分辨率平衡矩阵h，计算深度分辨率平衡优度函数drb优选出排位靠前的观测装置,并将其添加到构成初始集的该子集中形成新的子集；
④
计算相对模型分辨率m
r
，求出平均相对模型分辨率，判断是否满足优化要求；
⑤
若不满足要求，将新的子集作为优化的初始集s
i
，返回步骤
②
继续重复以上步骤，若满足条件，则输出当前的子集作为优化的最终结果。
[0056]
具体地，优化所依赖的基础综合集s
c
即为82个电极全部参与供电与测量时所产生的1560个电位数据。本次优化选取掌子面上下两个测量电极所能采集的全部电位数据作为初始集s
i
,其中包含有60个电位数据。
[0057]
计算模型分辨率矩阵m，其线性计算方法与估算公式为：
[0058]
m＝(g
t
g+c)-1
g
t
g
[0059]
式中，m为所要计算的模型分辨率矩阵，g为雅各比矩阵，c为约束矩阵。实施例中约束项使用常规约束c＝λi，其中约束因子λ取值为2.5
×
10-6
。
[0060]
将其余每个测量电极作为一个观测装置分别添加到初始集中，形成80个临时子集s
t
，并依据模型分辨率矩阵计算公式计算出相应的模型分辨率矩阵。
[0061]
所述的深度分辨率矩阵具体为：
[0062][0063]
深度分辨率平衡矩阵h中包含了所有网格的分辨率平衡因子，通过加入这一矩阵调整不同网格模型分辨率在drb计算过程中的分辨率平衡度。在钻孔电阻率观测装置模型
分辨率分布的初步研究中发现，随着深度的不断增加，深部区域网格的模型分辨率会不断减小。为了保证深部的模型分辨率不致过分损失，加入深度分辨率平衡矩阵。在实施例中反演深度为60m，因此h取为30m，当深度大于30m时，令所有网格的分辨率平衡度因子h为1.2，而隧道前方30m范围内网格的分辨率平衡度因子取值为1.0，具体的示意图如图5所示。
[0064]
利用drb函数将除初始集内的所有观测装置进行排序，深度分辨率平衡优度函数具体为：
[0065][0066]
m
t
储存着临时子集模型分辨率的主对角线元素,表示临时子集的分辨率向量，m
b
储存着初始集模型分辨率向量的主对角线元素，表示初始集的分辨率向量,h为深度分辨率平衡矩阵向量。m表示上述向量中的元素个数，m
t
(j)表示临时子集分辨率向量中的第j个元素。drb函数的值表征了新加入观测装置对原有基本集模型分辨率的提升,在一次迭代中选取排序最靠前的一个观测装置加入初始集中形成新的初始集。即在实施例中每次优化迭代结束时将会有一个新的测量电极点加入新的初始集中。
[0067]
计算新的初始集的相对模型分辨率，其公式为：
[0068][0069]
其中m
b
与m
c
分别储存着此时的初始集与综合集模型分辨率矩阵的主对角线元素，分别表示为初始集分辨率向量和综合集分辨率向量。公式中的除法代表m
b
中的每个元素分别除以m
c
对应位置的元素。所得结果m
r
即衡量了此时的初始集模型分辨率与综合集模型分辨率的接近程度，其元素均是介于0-1之间的数值，且显然当这些值越接近于1时表明此时的初始集的模型分辨率与综合集的越接近，其模型分辨率越好。
[0070]
计算新的初始集的平均相对模型分辨率，判断新的初始集的平均相对模型分辨率是否满足优化要求，即平均模型分辨率大小是否不小于0.9。
[0071]
具体地，在本实施例中，一次优化不能满足优化要求，将其余的每个观测电极分别作为一个观测装置添加到该更新后的初始集中，构成79个临时子集，重复以上步骤。
[0072]
具体地，本实施例中当优化迭代次数为20次时，平均相对模型分辨率刚好满足了大于0.9的要求，平均相对模型分辨率随优化迭代次数变化的曲线如附图3所示。图6为优化20次后的相对模型分辨率分布图，深处各部分的相对模型分辨率数值接近于1，表明优选出的子集的模型分辨率与综合集非常接近，深部的模型分辨率得到很大的提升。
[0073]
(4)根据优化结果确定掌子面测量电极的数目和位置，得到优化后的钻孔电阻率超前探测的有效观测模式。
[0074]
得到优化后的掌子面电极位置和数目，如图4所示。该图表明，优化后的掌子面测量电极相较于优化前，去除了较多多余电极，中间竖直的一列测量电极全被选中，其余被选中电极则较为对称地分布在四周边缘处，掌子面上共包含有22个测量电极，这样的测量电极布置能够较为均匀的提升整体的模型分辨率，而且结合钻孔中的30个供电电极共能采集660个电位数据，小于综合集数据的一半，观测效率大大提升，提高了工程适用性。
[0075]
图7是进行反演计算的地电模型，反演区域为30m*30m*60m，背景电阻率设置为1000ω.m，其中设置了一个低阻异常体，尺寸为4m*5m*6m，其电阻率为10ω.m。得到反演结果图，如图8所示，该数值模拟结果表明，该隧道钻孔电阻率超前探测装置深度分辨率平衡优化方法在精简了电极数目，减少了数据量，提高了反演效率的基础上，依然保证了反演成像的效果。
[0076]
以上所述仅为本公开的优选实施例而已，并不用于限制本公开，对于本领域的技术人员来说，本公开可以有各种更改和变化。凡在本公开的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本公开的保护范围之内。
[0077]
上述虽然结合附图对本公开的具体实施方式进行了描述，但并非对本公开保护范围的限制，所属领域技术人员应该明白，在本公开的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本公开的保护范围以内。