CN112415602B - 基于深度分辨率的隧道电阻率超前探测优化方法及*** - Google Patents

基于深度分辨率的隧道电阻率超前探测优化方法及*** Download PDF

Info

Publication number
CN112415602B
CN112415602B CN202011105343.6A CN202011105343A CN112415602B CN 112415602 B CN112415602 B CN 112415602B CN 202011105343 A CN202011105343 A CN 202011105343A CN 112415602 B CN112415602 B CN 112415602B
Authority
CN
China
Prior art keywords
resolution
initial set
model
resistivity
depth
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202011105343.6A
Other languages
English (en)
Other versions
CN112415602A (zh
Inventor
刘征宇
白鹏
张永恒
高雪池
许孝滨
庞永昊
马川义
聂利超
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Shandong University
Shandong High Speed Group Co Ltd
Original Assignee
Shandong University
Shandong High Speed Group Co Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Shandong University, Shandong High Speed Group Co Ltd filed Critical Shandong University
Priority to CN202011105343.6A priority Critical patent/CN112415602B/zh
Publication of CN112415602A publication Critical patent/CN112415602A/zh
Priority to PCT/CN2021/124212 priority patent/WO2022078516A1/zh
Priority to US18/031,715 priority patent/US20230384473A1/en
Application granted granted Critical
Publication of CN112415602B publication Critical patent/CN112415602B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/38—Processing data, e.g. for analysis, for interpretation, for correction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01R—MEASURING ELECTRIC VARIABLES; MEASURING MAGNETIC VARIABLES
    • G01R27/00—Arrangements for measuring resistance, reactance, impedance, or electric characteristics derived therefrom
    • G01R27/02—Measuring real or complex resistance, reactance, impedance, or other two-pole characteristics derived therefrom, e.g. time constant
    • G01R27/20—Measuring earth resistance; Measuring contact resistance, e.g. of earth connections, e.g. plates
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V3/00—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation
    • G01V3/02—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with propagation of electric current

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)
  • Investigating Or Analyzing Materials By The Use Of Electric Means (AREA)

Abstract

本公开提供了一种基于深度分辨率的隧道电阻率超前探测优化方法及***，将采集的所有电极测量数据整合为综合集，计算综合集的模型分辨率矩阵，从综合集中选择若干电极测量数据形成初始集；将不在初始集中电极的测量数据添加到初始集中，形成多个临时子集；根据模型分辨率矩阵，计算每个临时子集的深度分辨率平衡优度函数，选择深度分辨率平衡优度函数值最优的临时子集作为新的初始集，判断此时的初始集的平均相对模型分辨率是否满足优化要求；若不满足要求，则不断更新临时子集，否则输出此时的初始集；根据更新后的初始集确定掌子面测量电极的数目和位置，得到优化后的钻孔电阻率超前探测的有效观测模式。

Description

基于深度分辨率的隧道电阻率超前探测优化方法及***
技术领域
本公开属于隧道电阻率超前探测观测技术领域，具体涉及一种基于深度分辨率的隧道电阻率超前探测优化方法及***。
背景技术
本部分的陈述仅仅是提供了与本公开相关的背景技术信息，不必然构成在先技术。
随着隧道建设重点逐步向地形地质复杂的西部山区和水域阻隔的东部海峡地区转移，艰险山区深埋隧道和跨江越海隧道大量涌现，对于各类具有洞线长、埋深大等特点的深长隧道，对其进行地质调查更加困难。对于深埋长大隧道，由于各类不良地质赋存更隐蔽、致灾性更高，要求分辨率更高、目标更难探，对于隧道超前探测精度提出了更高的要求和挑战。目前传统的地表勘查条件和技术水平难以满足工程对探测深度和精细化探查的需求，而当前主流的隧道电阻率超前探测方法对含水体成像边界识别模糊、探不到、探不准，且对于一些亚米级甚至分米级的含导水构造难以有效识别，因此，需要一种新的隧道电阻率超前探测装置，能够对亚米级含水构造进行精细成像，即“孔中电极供电、掌子面阵列式测量”的隧道钻孔电阻率精细超前探测观测装置。
与隧道常规观测模式相比，隧道钻孔电阻率超前探测观测装置在精细化成像方面更具优势：1)具有更远的测深；2)可以获取钻孔周边异常体更丰富更直接的信息。因此，隧道钻孔电阻率超前探测在隧道电阻率超前探测中具有广阔的应用前景。
然而，据发明人了解，目前隧道钻孔电阻率超前探测基本观测模式仍存在以下关键问题尚未解决：
为了保证反演成像的效果，需要通过布设大量的电极来获得大量的观测数据，观测工作需要耗费大量的时间和人力，基本观测模式得到的大量数据对反演计算造成了极大的挑战，降低了观测效率。
同时，当前常规优化方法是提升整体的模型分辨率，是均一的提升，不能解决深部模型分辨率过分损失的问题。
发明内容
本公开为了解决上述问题，提出了一种基于深度分辨率的隧道电阻率超前探测优化方法及***，本公开对观测***进行优化，优选出对模型分辨率贡献较大的测量电极点，在保证模型分辨率的基础上最终得到优化后的掌子面电极排列，精简了电极数量，提高了探测效率。
根据一些实施例，本公开采用如下技术方案：
一种基于深度分辨率的隧道电阻率超前探测优化方法，包括以下步骤：
(1)将采集的所有电极测量数据整合为综合集；
(2)计算综合集的模型分辨率矩阵，从综合集中选择若干电极测量数据形成初始集；
(3)将不在初始集中电极的测量数据添加到初始集中，形成多个临时子集；
(4)根据模型分辨率矩阵，计算每个临时子集的深度分辨率平衡优度函数，选择深度分辨率平衡优度函数值最优的临时子集作为新的初始集，判断此时的初始集的平均相对模型分辨率是否满足优化要求；若不满足要求，则返回步骤(3)，否则输出此时的初始集；
(5)根据更新后的初始集确定掌子面测量电极的数目和位置，得到优化后的钻孔电阻率超前探测的有效观测模式。
作为可选择的实施方式，所述步骤(1)中，在隧道掌子面上布置钻孔，钻孔内布置电极，确定钻孔长度和电极间距；依据探测精度要求，确定正反演网格大小和掌子面上的电极数量和排列。
作为可选择的实施方式，所述步骤(1)中，使用全空间孔隧电阻率方法进行数据采集。
作为可选择的实施方式，所述步骤(2)中，选取掌子面上下两个测量电极所能采集的全部电位数据作为初始集。
作为可选择的实施方式，所述步骤(4)中，深度分辨率平衡矩阵由分辨率平衡因子组成，且各分辨率平衡因子根据反演深度确定。
作为可选择的实施方式，所述步骤(4)中，每个临时子集的深度分辨率平衡优度函数为：
Figure GDA0003191270210000031
作为可选择的实施方式，所述步骤(4)中，当前初始集的相对模型分辨率为此时的初始集模型分辨率矩阵的主对角线元素与综合集模型分辨率矩阵的主对角线元素的比值。
作为可选择的实施方式，所述步骤(4)中，平均相对模型分辨率为对相对模型分辨率中的元素求平均得到的。
作为可选择的实施方式，所述步骤(4)中，判断此时的初始集的平均相对模型分辨率是否满足优化要求的具体过程是，判断此时的初始集的平均相对模型分辨率的值是否大于设定值。
一种基于深度分辨率的隧道电阻率超前探测优化***，包括：
用于将采集的所有电极测量数据整合为综合集的模块；
用于计算综合集的模型分辨率矩阵，从综合集中选择若干电极测量数据形成初始集的模块；
用于将不在初始集中电极的测量数据添加到初始集中，形成多个临时子集的模块；
用于根据模型分辨率矩阵，计算每个临时子集的深度分辨率平衡优度函数，选择深度分辨率平衡优度函数值最优的临时子集作为新的初始集，判断此时的初始集的平均相对模型分辨率是否满足优化要求；若不满足要求，则返回重新添加另外电极的测量数据形成新的临时子集，否则输出此时的初始集的模块；
用于根据更新后的初始集确定掌子面测量电极的数目和位置的模块。
一种隧道电阻率超前探测观测***，包括多个测量电极，所述测量电极的数目和位置根据所述一种基于深度分辨率的隧道电阻率超前探测优化方法确定。
与现有技术相比，本公开的有益效果为：
1、本公开一种隧道电阻率超前探测装置形式优化方法，本公开在孔中电极等间距均匀布置、掌子面阵列电极测量的基本型式的基础上，使用深度分辨率平衡优化的DRB方法对观测装置进行优化，优选出对模型分辨率贡献较大的测量电极点，在保证模型分辨率的基础上最终得到优化后的掌子面电极排列，精简了电极数量，提高了探测效率。
2、本公开在深度分辨率平衡优度函数中，引入了深度分辨率平衡矩阵H，可以结合先验信息，调整不同网格模型分辨率在DRB计算过程中的分辨率平衡度因子，可以通过适当增大某一深度区域的分辨率平衡度，使得这一区域网格的模型分辨率在DRB的计算中更为重要，最终使得算法能够优先提升这一深度区域的模型分辨率。
附图说明
构成本公开的一部分的说明书附图用来提供对本公开的进一步理解，本公开的示意性实施例及其说明用于解释本公开，并不构成对本公开的不当限定。
图1是基于隧道钻孔电阻率超前探测装置深度分辨率平衡优化方法流程图；
图2是钻孔电阻率观测装置基本模式电极分布示意图；
图3是平均相对模型分辨率随迭代次数的变化曲线；
图4是优化后观测装置的掌子面测量电极分布示意图；
图5是深度分辨率平衡因子分布示意图；
图6是优化后的相对模型分辨率图；
图7是进行数值模拟时使用的地电模型图；
图8是基于优化后的隧道钻孔观测模型装置的反演结果图。
具体实施方式：
下面结合附图与实施例对本公开作进一步说明。
应该指出，以下详细说明都是例示性的，旨在对本公开提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本公开所属技术领域的普通技术人员通常理解的相同含义。
需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本公开的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
在本公开中，术语如“上”、“下”、“左”、“右”、“前”、“后”、“竖直”、“水平”、“侧”、“底”等指示的方位或位置关系为基于附图所示的方位或位置关系，只是为了便于叙述本公开各部件或元件结构关系而确定的关系词，并非特指本公开中任一部件或元件，不能理解为对本公开的限制。
本公开中，术语如“固接”、“相连”、“连接”等应做广义理解，表示可以是固定连接，也可以是一体地连接或可拆卸连接；可以是直接相连，也可以通过中间媒介间接相连。对于本领域的相关科研或技术人员，可以根据具体情况确定上述术语在本公开中的具体含义，不能理解为对本公开的限制。
一种隧道电阻率超前探测装置深度分辨率平衡优化方法，如图1所示，包括以下步骤：
(1)在隧道掌子面上布置钻孔，钻孔内布置电极，确定钻孔长度和电极间距；依据探测精度要求，确定正反演网格大小和掌子面上的电极数量和排列；
在具体实例中，当然，需要说明的是，本实施例仅为示例，并不仅代表其他实施例的钻孔深度、电极数量、间距和数据总数等需要和本实施例一致，上述数据都是可以根据具体观测要求和环境进行合理变更的，此为本领域技术人员容易想到的，理应属于本公开的保护范围。
钻孔最大深度设定为60m，孔内供电电极以等间距的方式排布，电极间间距为2m，则钻孔内共有30个供电电极。掌子面测量电极则呈阵列式排布，如图2所示共设置有8排电极，其中第1与第8排各自仅有3个电极，第4与第5排各自有9个电极，其余各排各自拥有7个电极。
(2)使用全空间孔隧电阻率方法进行数据采集；
具体地，数据采集采取钻孔内供电、掌子面接收的二极法型式，共采集到1560个测量数据。
(3)对掌子面测量电极进行优化和精简：①计算综合集的模型分辨率矩阵Mc，选取某个子集作为优化的初始集Si；②将每个其他的观测装置与初始集结合形成众多临时子集St，计算出相应的模型分辨率矩阵Mt；③确定深度分辨率平衡矩阵H，计算深度分辨率平衡优度函数DRB优选出排位靠前的观测装置，并将其添加到构成初始集的该子集中形成新的子集；④计算相对模型分辨率Mr，求出平均相对模型分辨率，判断是否满足优化要求；⑤若不满足要求，将新的子集作为优化的初始集Si，返回步骤②继续重复以上步骤，若满足条件，则输出当前的子集作为优化的最终结果。
具体地，优化所依赖的基础综合集Sc即为82个电极全部参与供电与测量时所产生的1560个电位数据。本次优化选取掌子面上下两个测量电极所能采集的全部电位数据作为初始集Si，其中包含有60个电位数据。
计算模型分辨率矩阵M，其线性计算方法与估算公式为：
M＝(GTG+C)-1GTG
式中，M为所要计算的模型分辨率矩阵，G为雅各比矩阵，C为约束矩阵。实施例中约束项使用常规约束C＝λI，其中约束因子λ取值为2.5×10-6。
将其余每个测量电极作为一个观测装置分别添加到初始集中，形成80个临时子集St，并依据模型分辨率矩阵计算公式计算出相应的模型分辨率矩阵。
所述的深度分辨率平衡矩阵具体为：
Figure GDA0003191270210000081
深度分辨率平衡矩阵H中包含了所有网格的分辨率平衡因子，通过加入这一矩阵调整不同网格模型分辨率在DRB计算过程中的分辨率平衡度。在钻孔电阻率观测装置模型分辨率分布的初步研究中发现，随着深度的不断增加，深部区域网格的模型分辨率会不断减小。为了保证深部的模型分辨率不致过分损失，加入深度分辨率平衡矩阵。在实施例中反演深度为60m，因此h取为30m，当深度大于30m时，令所有网格的分辨率平衡度因子H为1.2，而隧道前方30m范围内网格的分辨率平衡度因子取值为1.0，具体的示意图如图5所示。
利用DRB函数将除初始集内的所有观测装置进行排序，深度分辨率平衡优度函数具体为：
Figure GDA0003191270210000082
Mt储存着临时子集模型分辨率的主对角线元素，表示临时子集的分辨率向量，Mb储存着初始集模型分辨率向量的主对角线元素，表示初始集的分辨率向量，H为深度分辨率平衡矩阵向量。m表示上述向量中的元素个数，Mt(j)表示临时子集分辨率向量中的第j个元素。DRB函数的值表征了新加入观测装置对原有基本集模型分辨率的提升，在一次迭代中选取排序最靠前的一个观测装置加入初始集中形成新的初始集。即在实施例中每次优化迭代结束时将会有一个新的测量电极点加入新的初始集中。
计算新的初始集的相对模型分辨率，其公式为：
Figure GDA0003191270210000083
其中Mb与Mc分别储存着此时的初始集与综合集模型分辨率矩阵的主对角线元素，分别表示为初始集分辨率向量和综合集分辨率向量。公式中的除法代表Mb中的每个元素分别除以Mc对应位置的元素。所得结果Mr即衡量了此时的初始集模型分辨率与综合集模型分辨率的接近程度，其元素均是介于0-1之间的数值，且显然当这些值越接近于1时表明此时的初始集的模型分辨率与综合集的越接近，其模型分辨率越好。
计算新的初始集的平均相对模型分辨率，判断新的初始集的平均相对模型分辨率是否满足优化要求，即平均模型分辨率大小是否不小于0.9。
具体地，在本实施例中，一次优化不能满足优化要求，将其余的每个观测电极分别作为一个观测装置添加到该更新后的初始集中，构成79个临时子集，重复以上步骤。
具体地，本实施例中当优化迭代次数为20次时，平均相对模型分辨率刚好满足了大于0.9的要求，平均相对模型分辨率随优化迭代次数变化的曲线如附图3所示。图6为优化20次后的相对模型分辨率分布图，深处各部分的相对模型分辨率数值接近于1，表明优选出的子集的模型分辨率与综合集非常接近，深部的模型分辨率得到很大的提升。
(4)根据优化结果确定掌子面测量电极的数目和位置，得到优化后的钻孔电阻率超前探测的有效观测模式。
得到优化后的掌子面电极位置和数目，如图4所示。该图表明，优化后的掌子面测量电极相较于优化前，去除了较多多余电极，中间竖直的一列测量电极全被选中，其余被选中电极则较为对称地分布在四周边缘处，掌子面上共包含有22个测量电极，这样的测量电极布置能够较为均匀的提升整体的模型分辨率，而且结合钻孔中的30个供电电极共能采集660个电位数据，小于综合集数据的一半，观测效率大大提升，提高了工程适用性。
图7是进行反演计算的地电模型，反演区域为30m*30m*60m，背景电阻率设置为1000Ω.m，其中设置了一个低阻异常体，尺寸为4m*5m*6m，其电阻率为10Ω.m。得到反演结果图，如图8所示，该数值模拟结果表明，该隧道钻孔电阻率超前探测装置深度分辨率平衡优化方法在精简了电极数目，减少了数据量，提高了反演效率的基础上，依然保证了反演成像的效果。
以上所述仅为本公开的优选实施例而已，并不用于限制本公开，对于本领域的技术人员来说，本公开可以有各种更改和变化。凡在本公开的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本公开的保护范围之内。上述虽然结合附图对本公开的具体实施方式进行了描述，但并非对本公开保护范围的限制，所属领域技术人员应该明白，在本公开的技术方案的基础上，本领域技术人员不需要付出创造性劳动即可做出的各种修改或变形仍在本公开的保护范围以内。

Claims (8)

1.一种基于深度分辨率的隧道电阻率超前探测优化方法，其特征是：包括以下步骤：
(1)将采集的所有电极测量数据整合为综合集；
(2)计算综合集的模型分辨率矩阵，从综合集中选择若干电极测量数据形成初始集；
(3)将不在初始集中电极的测量数据添加到初始集中，形成多个临时子集；
(4)根据模型分辨率矩阵，计算每个临时子集的深度分辨率平衡优度函数，选择深度分辨率平衡优度函数值最优的临时子集作为新的初始集，判断此时的初始集的平均相对模型分辨率是否满足优化要求；若不满足要求，则返回步骤(3)，否则输出此时的初始集；
(5)根据更新后的初始集确定掌子面测量电极的数目和位置，得到优化后的钻孔电阻率超前探测的有效观测模式；
所述步骤(4)中，每个临时子集的深度分辨率平衡优度函数为：
Figure FDA0003191270200000011
Mt储存着临时子集模型分辨率的主对角线元素,表示临时子集的分辨率向量，Mb储存着初始集模型分辨率向量的主对角线元素，表示初始集的分辨率向量,H为深度分辨率平衡矩阵向量；m表示上述向量中的元素个数，Mt(j)表示临时子集分辨率向量中的第j个元素；
所述步骤(4)中，此时的初始集的相对模型分辨率为此时的初始集模型分辨率矩阵的主对角线元素与综合集模型分辨率矩阵的主对角线元素的比值；
所述步骤(4)中，平均相对模型分辨率计算为对相对模型分辨率中的元素求平均得到的。
2.如权利要求1所述的一种基于深度分辨率的隧道电阻率超前探测优化方法，其特征是：所述步骤(1)中，在隧道掌子面上布置钻孔，钻孔内布置电极，确定钻孔长度和电极间距；依据探测精度要求，确定正反演网格大小和掌子面上的电极数量和排列。
3.如权利要求1所述的一种基于深度分辨率的隧道电阻率超前探测优化方法，其特征是：所述步骤(1)中，使用全空间孔隧电阻率方法进行数据采集。
4.如权利要求1所述的一种基于深度分辨率的隧道电阻率超前探测优化方法，其特征是：所述步骤(2)中，选取掌子面上下两个测量电极所能采集的全部电位数据作为初始集。
5.如权利要求1所述的一种基于深度分辨率的隧道电阻率超前探测优化方法，其特征是：所述步骤(4)中，深度分辨率平衡矩阵由分辨率平衡因子组成，且各分辨率平衡因子根据反演深度确定。
6.如权利要求1所述的一种基于深度分辨率的隧道电阻率超前探测优化方法，其特征是：所述步骤(4)中，判断此时的初始集的平均相对模型分辨率是否满足优化要求的具体过程是，判断此时的初始集的平均相对模型分辨率的值是否大于设定值。
7.一种基于深度分辨率的隧道电阻率超前探测优化***，其特征是：包括：
用于将采集的所有电极测量数据整合为综合集的模块；
用于计算综合集的模型分辨率矩阵，从综合集中选择若干电极测量数据形成初始集的模块；
用于将不在初始集中电极的测量数据添加到初始集中，形成多个临时子集的模块；
用于根据模型分辨率矩阵，计算每个临时子集的深度分辨率平衡优度函数，选择深度分辨率平衡优度函数值最优的临时子集作为新的初始集，判断此时的初始集的平均相对模型分辨率是否满足优化要求；若不满足要求，则返回重新添加另外电极的测量数据形成新的临时子集，否则输出此时的初始集的模块；
每个临时子集的深度分辨率平衡优度函数为：
Figure FDA0003191270200000031
Mt储存着临时子集模型分辨率的主对角线元素,表示临时子集的分辨率向量，Mb储存着初始集模型分辨率向量的主对角线元素，表示初始集的分辨率向量,H为深度分辨率平衡矩阵向量；m表示上述向量中的元素个数，Mt(j)表示临时子集分辨率向量中的第j个元素；
此时的初始集的相对模型分辨率为此时的初始集模型分辨率矩阵的主对角线元素与综合集模型分辨率矩阵的主对角线元素的比值；
平均相对模型分辨率计算为对相对模型分辨率中的元素求平均得到的；
用于根据更新后的初始集确定掌子面测量电极的数目和位置的模块。
8.一种隧道电阻率超前探测观测***，其特征是：包括多个测量电极，所述测量电极的数目和位置根据权利要求1-6中任一项所述的一种基于深度分辨率的隧道电阻率超前探测优化方法确定。
CN202011105343.6A 2020-10-15 2020-10-15 基于深度分辨率的隧道电阻率超前探测优化方法及*** Active CN112415602B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN202011105343.6A CN112415602B (zh) 2020-10-15 2020-10-15 基于深度分辨率的隧道电阻率超前探测优化方法及***
PCT/CN2021/124212 WO2022078516A1 (zh) 2020-10-15 2021-10-15 基于深度分辨率的隧道电阻率超前探测优化方法及***
US18/031,715 US20230384473A1 (en) 2020-10-15 2021-10-15 Method and system for advanced detection and optimization of tunnel resistivity based on depth resolution

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011105343.6A CN112415602B (zh) 2020-10-15 2020-10-15 基于深度分辨率的隧道电阻率超前探测优化方法及***

Publications (2)

Publication Number Publication Date
CN112415602A CN112415602A (zh) 2021-02-26
CN112415602B true CN112415602B (zh) 2021-11-19

Family

ID=74854980

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011105343.6A Active CN112415602B (zh) 2020-10-15 2020-10-15 基于深度分辨率的隧道电阻率超前探测优化方法及***

Country Status (3)

Country Link
US (1) US20230384473A1 (zh)
CN (1) CN112415602B (zh)
WO (1) WO2022078516A1 (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112415602B (zh) * 2020-10-15 2021-11-19 山东大学 基于深度分辨率的隧道电阻率超前探测优化方法及***
CN113309506B (zh) * 2021-05-18 2023-02-03 山东大学 基于孔中电偶极子发射的超前观测方法与装置
CN114778948A (zh) * 2022-06-17 2022-07-22 中铁大桥科学研究院有限公司 动水隧道岩体电阻率监测方法及相关设备

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6353801B1 (en) * 1999-04-09 2002-03-05 Agilent Technologies, Inc. Multi-resolution adaptive solution refinement technique for a method of moments-based electromagnetic simulator
CN102759751A (zh) * 2012-07-30 2012-10-31 山东大学 地下工程高分辨率三维电阻率ct成像超前预报***和方法
CN106570227A (zh) * 2016-10-20 2017-04-19 湖南师范大学 一种超高密度电法的电极排列优化方法及装置
CN106772621A (zh) * 2017-01-24 2017-05-31 山东大学 一种近全方位电阻率隧道超前地质预报方法
CN107742015A (zh) * 2017-09-30 2018-02-27 中南大学 基于任意偶极‑偶极装置的直流激电法三维数值模拟方法
CN108776355A (zh) * 2018-07-20 2018-11-09 山东大学 隧道聚焦测深型三维激发极化超前探测仪器***
CN111323830A (zh) * 2020-01-14 2020-06-23 东华理工大学 一种基于大地电磁和直流电阻率数据的联合反演方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2402217A1 (fr) * 1977-09-01 1979-03-30 Geophysique Cie Gle Prospection electrique du sous-sol, a electrodes lineiques
JP2001166061A (ja) * 1999-12-10 2001-06-22 Ohbayashi Corp トンネル切羽前方の探査方法
US9146330B2 (en) * 2011-03-29 2015-09-29 Westerngeco L.L.C. Selecting a survey setting for characterizing a target structure
CN102798896B (zh) * 2011-05-27 2015-11-18 中国石油天然气集团公司 一种阵列感应测井仪器的测井信号合成处理方法及其***
US9746577B2 (en) * 2013-06-03 2017-08-29 Exxonmobil Upstream Research Company Uncertainty estimation of subsurface resistivity solutions
CN103645514B (zh) * 2013-12-25 2015-01-07 山东大学 多同性源电极阵列电阻率的地下工程超前探测方法及***
US9715034B2 (en) * 2015-12-18 2017-07-25 Schlumberger Technology Corporation Method for multi-tubular evaluation using induction measurements
CN108169801B (zh) * 2018-01-16 2020-09-15 陕西铁道工程勘察有限公司 高分辨地电阻率快速成像方法
CN110908000B (zh) * 2019-11-07 2021-10-19 吉林大学 基于变维贝叶斯的隧道瞬变电磁数据解释方法
CN112415602B (zh) * 2020-10-15 2021-11-19 山东大学 基于深度分辨率的隧道电阻率超前探测优化方法及***

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6353801B1 (en) * 1999-04-09 2002-03-05 Agilent Technologies, Inc. Multi-resolution adaptive solution refinement technique for a method of moments-based electromagnetic simulator
CN102759751A (zh) * 2012-07-30 2012-10-31 山东大学 地下工程高分辨率三维电阻率ct成像超前预报***和方法
CN106570227A (zh) * 2016-10-20 2017-04-19 湖南师范大学 一种超高密度电法的电极排列优化方法及装置
CN106772621A (zh) * 2017-01-24 2017-05-31 山东大学 一种近全方位电阻率隧道超前地质预报方法
CN107742015A (zh) * 2017-09-30 2018-02-27 中南大学 基于任意偶极‑偶极装置的直流激电法三维数值模拟方法
CN108776355A (zh) * 2018-07-20 2018-11-09 山东大学 隧道聚焦测深型三维激发极化超前探测仪器***
CN111323830A (zh) * 2020-01-14 2020-06-23 东华理工大学 一种基于大地电磁和直流电阻率数据的联合反演方法

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
"Practical aspects of applied optimized survey design for electrical resistivity tomography";Wilkinson P B等;《Geophysical Journal International》;20121231;第189卷(第1期);第428-440页 *
"不同噪声水平高密度电法的分辨率和勘探深度研究";邢润林 等;《物探化探计算技术》;20190131;第41卷(第1期);第68-79页 *
"电极布置方式对高密度电法探测分辨率的影响";杨磊 等;《地球物理学进展》;20191231;第34卷(第1期);第0406-0411页 *

Also Published As

Publication number Publication date
CN112415602A (zh) 2021-02-26
WO2022078516A1 (zh) 2022-04-21
US20230384473A1 (en) 2023-11-30

Similar Documents

Publication Publication Date Title
CN112415602B (zh) 基于深度分辨率的隧道电阻率超前探测优化方法及***
CN110702587B (zh) 一种基于温纳联合反演的土石坝渗漏诊断方法
CN113703068B (zh) 确定铀矿分布的方法
CN112949134B (zh) 基于非结构有限元方法的地-井瞬变电磁反演方法
CN110187398A (zh) 一种寻找井间目标体的多电极系探测方法
CN110988999A (zh) 基于跨孔电阻率法ct反演成像分析桩基的探测方法及***
CN106526329A (zh) 一种测量地表到深层大地电阻率的方法及装置
CN104350984A (zh) 堤坝白蚁防治技术
CN105136866A (zh) 基于高密度电阻率法的污染土检测方法
CN103065051A (zh) 一种对岩体自动进行分级分区的方法
CN206378448U (zh) 基于综合物探方法的污染土检测测线布置结构
CN105321204A (zh) 一种三维地质栅格体模型构建方法
CN106802432A (zh) 土层结构的探查方法和装置
He et al. Analysis of the anisotropic spatial variability and three-dimensional computer simulation of agricultural soil bulk density in an alluvial plain of north China
CN106772620B (zh) 十字形三维电法勘探方法
CN109343131B (zh) 一种侧向供电的地井直流电法探测方法与装置
CN111721831A (zh) 基于电刺激的三维层析成像堤坝隐伏渗漏通道扫描方法
CN112485829A (zh) 一种电阻率法计算动态含水量方法
CN114236624B (zh) 基于电磁法估算压裂改造空间体积的方法和***
CN115982826A (zh) 一种深埋长大隧道纵剖面地应力参数值估计方法、装置、设备及存储介质
CN106646622B (zh) 井字形三维电法勘探方法
CN105572732A (zh) 一种属性变化率逐次晋级的裂缝发育带检测方法
CN112882124B (zh) 一种三维铀成矿有利部位的圈定方法
CN113516347B (zh) 基于地球物理与水文数据融合的地下水水质评估方法
CN205139068U (zh) 基于高密度电阻率法的污染土检测测线布置结构

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant