CN111077581B - 一种隧道突水三维核磁共振超前探测装置及成像方法 - Google Patents

一种隧道突水三维核磁共振超前探测装置及成像方法 Download PDF

Info

Publication number
CN111077581B
CN111077581B CN201910414235.8A CN201910414235A CN111077581B CN 111077581 B CN111077581 B CN 111077581B CN 201910414235 A CN201910414235 A CN 201910414235A CN 111077581 B CN111077581 B CN 111077581B
Authority
CN
China
Prior art keywords
tunnel
receiving
magnetic resonance
dimensional
coils
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
CN201910414235.8A
Other languages
English (en)
Other versions
CN111077581A (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.)
Jilin University
Original Assignee
Jilin University
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 Jilin University filed Critical Jilin University
Priority to CN201910414235.8A priority Critical patent/CN111077581B/zh
Publication of CN111077581A publication Critical patent/CN111077581A/zh
Application granted granted Critical
Publication of CN111077581B publication Critical patent/CN111077581B/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/14—Electric or magnetic prospecting or detecting; Measuring magnetic field characteristics of the earth, e.g. declination, deviation operating with electron or nuclear magnetic resonance

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Engineering & Computer Science (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明属于地球物理勘探领域，具体来讲为一种隧道突水三维核磁共振超前探测装置及成像方法。该装置包括：球形支架，所述球形支架上以球心为圆心建立直角坐标系，以隧道掘进方向为X轴正方向，垂直于X轴，指向隧道拱顶面的方向为Z轴正方向，与X轴和Z轴夹角均为90°，指向隧道边墙面的方向为Y轴正方向，每个轴向上设置有多个相同尺寸、相同匝数的接收线圈，相邻线圈之间夹角相同。采用一个固定方向的发射线圈，三轴多角度接收天线，准确、快速的得出隧道全空间内的含水构造的核磁共振响应，并利用成像方法得到含水构造分布及含水量大小，大大减小隧道突水等事故误报、漏报的风险，减小经济及生命财产损失。

Description

一种隧道突水三维核磁共振超前探测装置及成像方法
技术领域
本发明属于地球物理勘探领域，具体来讲为一种隧道突水三维核磁共振超前探测装置及成像方法。
背景技术
在隧道开凿过程中，需要提前探明掌子面前方赋存的断层和暗河等地质构造，以防止突水等事故发生。核磁共振作为一种直接探测地下水的技术，具有直接、定量和解释唯一等优点，三维核磁共振成像结果不仅能够对灾害水源进行定量评估，而且能实现准确的定位。因此，隧道突水三维磁共振探测方法具有重要意义。
CN102262247A公开了一种隧道突水超前预测装置及方法，该方法利用收发一体多匝矩形或方形线圈，平行于掌子面布设。通过计算机控制高压电源，主控芯片通过驱动电路驱动发射桥路，产生激发磁场。放大电路将核磁共振信号放大后送到采集电路，采集电路通过转换器将放大器输出的模拟信号转换成数字信号后送至计算机，对数据进行显示、存储、滤波处理和解释，以实现掌子面附近含水体的一维探测。
CN102062877A公开了一种对前方水体超前探测的核磁共振装置及探测方法。计算机通过串口算机通过串口总线分别与***控制器、大功率电源、信号采集单元相连，***控制器经桥路驱动器、大功率型发射桥路和配谐电容与发射线圈连结构成。该方法利用收发分离线圈，垂直布设模式，虽然可以在狭小的空间中展开勘探工作，但该方法仍然属于一维探测，将隧道空间简化为层状模型，与实际情况不符。
CN104765072A公开了一种用环形天线旋转进行磁共振超前探测的方法。该方法将环形天线固定在可旋转支架上，采用垂直和水平旋转两种测量方式，在隧道和矿井中对掌子面前方、顶板和左右侧壁的含水体进行探测，实现掌子面前方扇形区域内的含水体二维探测。但是该方法采用收发同一线圈多次旋转装置，耗费大量测量时间，且不适用于溶洞等无法简化为二维模型的地质构造。
发明内容
本发明所要解决的技术问题在于提供一种隧道灾害水源三维核磁共振超前探测及成像方法，旨在采用一个固定方向的发射线圈，三轴多角度接收天线，准确、快速的得出隧道全空间内的含水构造的核磁共振响应，并利用成像方法得到含水构造分布及含水量大小，大大减小隧道突水等事故误报、漏报的风险，减小经济及生命财产损失。
本发明是这样实现的，
一种隧道突水三维核磁共振超前探测装置，该装置包括：
球形支架，所述球形支架上以球心为圆心建立直角坐标系，以隧道掘进方向为X轴正方向，垂直于X轴，指向隧道拱顶面的方向为Z轴正方向，与X轴和Z轴夹角均为90°，指向隧道边墙面的方向为Y轴正方向，每个轴向上设置有多个相同尺寸、相同匝数的接收线圈，相邻线圈之间夹角相同。
进一步地，X轴方向接收线圈具体摆放为：首先放置1个接收线圈，使其法向方向指向X轴方向，然后依次放置另外的接收线圈，使多个接收线圈间在Z轴方向的夹角均相同，且中心位置重叠。
进一步地，Y轴方向接收线圈具体摆放为：首先放置1个接收线圈，使其法向方向指向Y轴方向，然后依次放置另外的接收线圈，使多个接收线圈间在X轴方向的夹角均相同，且中心位置重叠。
进一步地，Z轴方向接收线圈具体摆放为：首先放置1个接收线圈，使其法向方向指向Z轴方向，然后依次放置另外的接收线圈，使多个接收线圈间在Y轴方向的夹角均相同，且中心位置重叠。
进一步地，设置一个多匝发射线圈平行于掌子面放置，使得法向方向指向X轴方向，同时，与球形支架结构保持距离，且发射线圈中心与球形支架中心点间为平行于X轴的直线，形成收发分离探测结构。
进一步地，每个轴向上设置6个接收线圈，相邻接收线圈之间夹角均为30°。
一种隧道突水三维核磁共振超前探测成像方法，该方法包括：从小到大依次在发射线圈内通入Q个不同脉冲矩的交变电流，每次发射停止后，三个轴向上多个接收线圈同时接收三维隧道空间内含水体产生的磁共振信号，共计Q次，则接收得到电压数据为：
Vobs＝[V1，1，V1,2，…，V1,Q，V2，1，V2,2，…，V2,Q，…，VRx，1，VRx,2，…，VRx,Q]T，T为转置运算。其中V1，1＝[V1,V2,…,VNS]，为第一个接收线圈对应第一个脉冲距接收到的电压数据，NS＝1f，f为采样频率；
将接收线圈Q个脉冲矩对应的电压数据分别进行尖峰噪声和工频谐波噪声消除，再经过包络检测和门积分，得到观测磁共振数据复数矩阵
Figure GDA0002998562820000031
其中，R和I分别代表实部和虚部，N代表门积分个数，Rx为接收线圈个数。
Figure GDA0002998562820000032
Figure GDA0002998562820000033
将不同接收方向线圈的观测磁共振数据复数矩阵V同时参与三维复数域QT反演，得到隧道全空间含水量三维成像结果。
进一步地，
三维复数域QT反演方法实现如下：
三维复数域QT反演的目标函数为
Figure GDA0002998562820000034
其中D＝1/ε是观测电压数据噪声水平估计加权值，Vcalc是利用正演方法计算得到的接收线圈多次旋转对应的磁共振信号集合，m是三维模型网格每个单元的含水量集合，C是三维模型平滑度矩阵，λ是正则化参数，用来平衡数据拟合度和模型平滑度。
其中
Figure GDA0002998562820000041
是经相位偏移修正后观测值，
Figure GDA0002998562820000042
为相位偏移量。利用共轭梯度方法求解反演目标函数，经过多次迭代，得到隧道三维模型空间的含水量分布。
本发明与现有技术相比，有益效果在于：本发明实现三维磁共振成像，克服一维探测及成像结果只能得到含水构造一个维度的层状信息，二维旋转天线探测方法需多次旋转测量，不仅耗费大量时间，且其成像结果只适用于含水断层等可以简化为二维模型的地质构造，例如含水断层。采用三轴多角度接收并进行三维成像，可以准确得到任何典型含水构造(例如含水断层、陷落的水柱和溶洞等)的位置及含水规模，减少隧道开凿过程带来的经济和生命财产损失。
附图说明
图1为本发明的测量装置结构示意图；
图2为X(a)、Y(b)、Z(c)三方向线圈结构分解示意图；
图3为X方向6个不同角度线圈接收到的电压值，a为第一接收线圈、b为第二接收线圈、c为第三接收线圈、d为第四接收线圈、e为第五接收线圈、f为第六接收线圈；
图4为隧道全空间内含水量三维分布图结果。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合实施例，对本发明进行进一步详细说明。应当理解，此处所描述的具体实施例仅仅用以解释本发明，并不用于限定本发明。
考虑到开凿隧道的空间尺寸，本实施例中隧道突水三维核磁共振超前探测装置设置共18个多匝接收线圈的直径为6m，匝数为30匝。首先，将18个接收天线固定于球形支架1上，X、Y、Z方向各6个，见图1形结构线圈所示，掌子面前方为含水构造3见图1。令隧道掘进方向为X轴正方向，垂直于X轴，指向隧道拱顶面的方向为Z轴正方向，与X轴和Z轴夹角均为90°，指向隧道边墙面的方向为Y轴正方向。每个方向上天线间的夹角为30°。并将固定好的球形天线结构放置于掌子面正前方。
X轴方向6个接收线圈具体摆放方式如下：首先放置1个接收线圈，使其法向方向指向X轴方向，然后依次放置另外5个接收线圈，使6个接收线圈间在Z轴方向的夹角均为30°，且中心位置重叠，见附图2a。Y轴方向6个接收线圈具体摆放方式如下：首先放置1个线圈，使其法向方向指向Y轴方向，然后依次放置另外5个接收线圈，使6个接收线圈间在X轴方向的夹角均为30°，且中心位置重叠，见附图2b。Z轴方向6个接收线圈具体摆放方式如下：首先放置1个接收线圈，使其法向方向指向Z轴方向，然后依次放置另外5个接收线圈，使6个接收线圈间在Y轴方向的夹角均为30°，且中心位置重叠，见附图2c。
将1个多匝发射线圈2平行与掌子面放置，即法向方向指向X轴方向，本实施例中发射线圈边长为6m，匝数为8匝，见图1中结构2。同时，与接收球形线圈结构保持6m距离，且发射线圈中心与球形发射中心点间为平行与X轴的直线，形成收发分离探测结构。
一种隧道突水三维核磁共振超前探测成像方法包括：从小到大依次在发射线圈内通入脉冲矩为0.01As～4As等对数间隔分布的20个值对应的交变电流，即Q＝20。每次发射停止后，18个多匝接收线圈同时接收三维隧道空间内含水体产生的磁共振信号，共计20次，则接收得到电压数据为
Vobs＝[V1，1，V1,2，…，V1,20，V2，1，V2,2，…，V2,20，…，V18，1，V18,2，…，V18,20]T，T为转置运算。其中V1，1＝[V1,V2,…,VNS]，为第一个接收线圈对应第一个脉冲距接收到的电压数据，NS＝1/f，f为采样频率，本实施例中，采样频率为10kHz。实验地点的Larmor频率为2300Hz，地磁倾角和地磁偏角分别为150°和90°。
将18个线圈Q个脉冲矩对应的电压数据分别进行尖峰噪声，工频谐波噪声消除，经过包络检测和门积分方法，得到观测核共振数据复数矩阵
Figure GDA0002998562820000061
其中，R和I分别代表实部和虚部，门积分个数N为30，接收线圈个数Rx等于18。
Figure GDA0002998562820000062
Figure GDA0002998562820000063
其中X方向6个接收线圈对应的实部和虚部初始振幅值见图3。
将18个不同接收方向线圈的数据V同时参与三维复数域QT反演，得到隧道全空间含水量三维成像结果，见图4所示。由图4可以看出，利用基于球形结构天线的三维核磁共振超前探测装置和方法得到的隧道全空间含水量三维成像结果与模型含水构造基本一致，验证了方法的有效性。
三维复数域QT反演方法实现如下：
三维复数域复数QT反演的目标函数为
Figure GDA0002998562820000064
其中D＝1/ε是观测信号噪声水平估计加权值，Vcalc是根据正演方法计算得到的接收线圈多次旋转对应的NMR信号集合，m是三维模型网格每个单元的含水量集合，C是三维模型平滑度矩阵，λ是正则化参数，用来平衡数据拟合度和模型平滑度，本实施例中设为100。
其中
Figure GDA0002998562820000065
是经相位偏移修正后观测值，
Figure GDA0002998562820000066
为相位偏移量。利用共轭梯度方法求解反演目标函数，经过多次迭代，即可得到隧道三维模型空间的含水量分布。
以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等，均应包含在本发明的保护范围之内。

Claims (6)

1.一种隧道突水三维核磁共振超前探测装置，其特征在于，该装置包括：
球形支架，所述球形支架上以球心为圆心建立直角坐标系，以隧道掘进方向为X轴正方向，垂直于X轴，指向隧道拱顶面的方向为Z轴正方向，与X轴和Z轴夹角均为90°，指向隧道边墙面的方向为Y轴正方向，每个轴向上设置有多个相同尺寸、相同匝数的接收线圈，相邻线圈所形成的平面之间夹角相同；
X轴方向接收线圈具体摆放为：首先放置1个接收线圈，使其法向方向指向X轴方向，然后依次放置另外的接收线圈，使相邻两个接收线圈所形成的平面之间的夹角均相同，且中心位置重叠；
Y轴方向接收线圈具体摆放为：首先放置1个接收线圈，使其法向方向指向Y轴方向，然后依次放置另外的接收线圈，使相邻两个接收线圈所形成的平面之间的夹角均相同，且中心位置重叠；
Z轴方向接收线圈具体摆放为：首先放置1个接收线圈，使其法向方向指向Z轴方向，然后依次放置另外的接收线圈，使相邻两个接收线圈所形成的平面之间的夹角均相同，且中心位置重叠。
2.按照权利要求1所述的隧道突水三维核磁共振超前探测装置，其特征在于，设置一个多匝发射线圈平行于掌子面放置，使得法向方向指向X轴方向，同时，与球形支架结构保持一定距离，且发射线圈中心与球形支架中心点间为平行于X轴的直线，形成收发分离探测结构。
3.按照权利要求1所述的隧道突水三维核磁共振超前探测装置，其特征在于，每个轴向上设置6个接收线圈，相邻接收线圈所形成的平面之间夹角均为30°。
4.采用权利要求1-3任意一项所述的装置的隧道突水三维核磁共振超前探测成像方法，其特征在于，该方法包括：从小到大依次在发射线圈内通入Q个不同脉冲矩的交变电流，每次发射停止后，三个轴向上多个接收线圈同时接收三维隧道空间内含水体产生的磁共振信号，共计Q次，则接收得到电压数据；
将接收线圈Q个脉冲矩对应的电压数据分别进行尖峰噪声和工频谐波噪声消除，再经过包络检测和门积分，得到观测磁共振数据复数矩阵V；
将不同接收方向线圈的观测磁共振数据复数矩阵V同时参与三维复数域QT反演，得到隧道全空间含水量三维成像结果。
5.按照权利要求4所述的方法，其特征在于，
三维复数域QT反演方法实现如下：
三维复数域QT反演的目标函数为
Figure FDA0002998562810000021
其中D＝1/ε是观测电压数据噪声水平估计加权值，Vcalc是利用正演方法计算得到的接收线圈多次旋转对应的磁共振信号集合，m是三维模型网格每个单元的含水量集合，C是三维模型平滑度矩阵，λ是正则化参数，用来平衡数据拟合度和模型平滑度；
其中
Figure FDA0002998562810000022
是经相位偏移修正后观测值，
Figure FDA0002998562810000023
为相位偏移量；利用共轭梯度方法求解反演目标函数，经过多次迭代，得到隧道三维模型空间的含水量分布。
6.按照权利要求4所述的方法，其特征在于，观测磁共振数据复数矩阵
Figure FDA0002998562810000024
其中，R和I分别代表实部和虚部，N代表门积分个数，Rx为接收线圈个数，T为转置运算，其中：
Figure FDA0002998562810000025
Figure FDA0002998562810000026
CN201910414235.8A 2019-05-17 2019-05-17 一种隧道突水三维核磁共振超前探测装置及成像方法 Active CN111077581B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910414235.8A CN111077581B (zh) 2019-05-17 2019-05-17 一种隧道突水三维核磁共振超前探测装置及成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910414235.8A CN111077581B (zh) 2019-05-17 2019-05-17 一种隧道突水三维核磁共振超前探测装置及成像方法

Publications (2)

Publication Number Publication Date
CN111077581A CN111077581A (zh) 2020-04-28
CN111077581B true CN111077581B (zh) 2021-08-27

Family

ID=70310301

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910414235.8A Active CN111077581B (zh) 2019-05-17 2019-05-17 一种隧道突水三维核磁共振超前探测装置及成像方法

Country Status (1)

Country Link
CN (1) CN111077581B (zh)

Families Citing this family (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113309575B (zh) * 2021-06-25 2022-03-25 中国矿业大学 基于核磁传感器的煤矿井下含水量实时监测***及方法
CN117075212B (zh) * 2023-10-16 2024-01-26 吉林大学 一种隧道磁共振裂隙结构成像方法
CN117233849B (zh) * 2023-11-14 2024-01-30 吉林大学 一种基于稳态进动序列的磁共振隧道灾害水快速定位方法

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102819046B (zh) * 2012-09-06 2015-08-05 吉林大学 双d型线圈核磁共振仪随掘进机巷道前方突水探测方法
US9417068B2 (en) * 2013-05-01 2016-08-16 Massachusetts Institute Of Technology Stable three-axis nuclear spin gyroscope
CN104597510B (zh) * 2015-01-29 2017-03-08 中国科学院电子学研究所 一种获取张量旋转不变量的频率域电磁法***
CN104776865B (zh) * 2015-03-21 2017-07-07 复旦大学 基于最大磁感应强度矢量旋转角快速测定的电磁跟踪***及方法
US9986630B2 (en) * 2015-05-26 2018-05-29 Krunoslav Subotic Superconducting magnet winding structures for the generation of iron-free air core cyclotron magnetic field profiles
CN205565905U (zh) * 2016-03-29 2016-09-07 中惠创智无线供电技术有限公司 谐振线圈耦合的3d球体模型
CN106971806B (zh) * 2017-04-18 2020-05-08 上海交通大学 一种核磁共振陀螺仪的三维Braunbeck驱动线圈

Also Published As

Publication number Publication date
CN111077581A (zh) 2020-04-28

Similar Documents

Publication Publication Date Title
CN111077581B (zh) 一种隧道突水三维核磁共振超前探测装置及成像方法
CN110007350B (zh) 一种磁探测方法盲区的分析方法
CN108873086B (zh) 一种采用地磁总场梯度阵列对磁性目标定位的方法
CN110308490B (zh) 基于磁传感器阵列的旋转运动磁异常目标探测装置及方法
CN104765072B (zh) 一种用环形天线旋转进行磁共振超前探测的方法
US11487038B2 (en) Operating method of a metal detector capable of measuring target depth
CN109696710A (zh) 基于多场信息反馈的矿井瞬变电磁自动化观测方法
CN201397397Y (zh) 多参量地下管道轨迹探测仪
CN109633490B (zh) 一种全张量磁梯度测量组件的标定方法
CN109141408A (zh) 一种长输地下管道定位***的误差补偿方法及实施***
CN103353612B (zh) 一种地下目标物体的测量定位设备及测量定位方法
CN109725360A (zh) 基于磁梯度张量不变量的单点定位方法
CN106873044B (zh) 阵列式squid核磁共振地下水探测装置及成像方法
Liu et al. Magneto-inductive magnetic gradient tensor system for detection of ferromagnetic objects
CN109633540B (zh) 一种磁源的实时定位***及实时定位方法
CN108287367A (zh) 一种基于时域电磁法的地面-孔中联合探测***及应用方法
CN103499810A (zh) 一种用于电磁定位的装置和方法
CN106772635A (zh) 一种基于扫描交汇的地下电缆三维定位装置及定位方法
CN104199056B (zh) 定位检测仪
CN116449442B (zh) 一种钻孔瞬变电磁法三维展开处理方法及***
CN112859185A (zh) 一种基于地球场核磁共振的非侵入性远程物质探测装置
CN109633510B (zh) 基于地面备份单机的卫星磁强计磁环境分析方法
CN109633539B (zh) 一种磁源的静态定位装置及静态定位方法
US10663614B1 (en) Magnetic anomaly sensing system using two triaxial magnetometer sensors
CN106908058B (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