CN107179556A - 一种核磁共振的超前探测方法、系统及电子设备 - Google Patents
一种核磁共振的超前探测方法、系统及电子设备 Download PDFInfo
- Publication number
- CN107179556A CN107179556A CN201710444249.5A CN201710444249A CN107179556A CN 107179556 A CN107179556 A CN 107179556A CN 201710444249 A CN201710444249 A CN 201710444249A CN 107179556 A CN107179556 A CN 107179556A
- Authority
- CN
- China
- Prior art keywords
- detectable signal
- magnetic resonance
- nuclear magnetic
- advance
- kernel function
- 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.)
- Granted
Links
- 238000005481 NMR spectroscopy Methods 0.000 title claims abstract description 101
- 238000000034 method Methods 0.000 title claims abstract description 48
- 239000000523 sample Substances 0.000 title claims abstract description 38
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 29
- 230000006870 function Effects 0.000 claims description 69
- 238000000354 decomposition reaction Methods 0.000 claims description 16
- 230000005284 excitation Effects 0.000 claims description 14
- 238000005267 amalgamation Methods 0.000 claims description 12
- 238000012545 processing Methods 0.000 claims description 9
- 238000004804 winding Methods 0.000 claims description 7
- 239000004744 fabric Substances 0.000 claims description 6
- 238000000205 computational method Methods 0.000 claims 2
- 230000005311 nuclear magnetism Effects 0.000 claims 1
- 238000001514 detection method Methods 0.000 abstract description 43
- 238000005516 engineering process Methods 0.000 description 7
- 210000004209 hair Anatomy 0.000 description 6
- 238000011835 investigation Methods 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000005259 measurement Methods 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 230000008859 change Effects 0.000 description 3
- 238000013316 zoning Methods 0.000 description 3
- 235000013399 edible fruits Nutrition 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000010295 mobile communication Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000004590 computer program Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
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
本发明公开了一种核磁共振的超前探测方法、系统及电子设备，其中，所述核磁共振的超前探测方法包括：获得核磁共振装置的当前探测信号；根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水层先验信息，获得所述掌子面前方区域的含水率分布。本发明解决了现有技术中存在探测的水体分布情况的精度不高的技术问题。
Description
技术领域
本发明涉及地球物理勘探技术领域，尤其涉及一种核磁共振的超前探测方法、系统及电子设备。
背景技术
核磁共振测深技术(Nuclear Magnetic Resonance NMR)由于其具有测量信息丰富、解释直接等特点，被广泛应用于勘探领域中。
现有技术中，利用NMR技术进行勘探时，例如隧道、巷道、矿井等大型的地下工程，为了保证安全性，掘进头的探测需要分段进行，且掘进头前进的距离小于探测深度，然后通过对掘进头每次探测的信息进行分析和计算，从而得到水体的分布情况。
本申请发明人在实现本发明的技术方案时，发现现有技术中至少存在如下问题：
现有方法单独对掘进头每次探测的信息进行分析和计算，从而得到水体的分布情况，而忽视了每次探测信息之间的关联，因此导致得出的水体分布情况的精度不高。
可见，现有技术的探测方法存在探测的水体分布情况的精度不高的技术问题。
发明内容
本发明实施例提供一种核磁共振的超前探测方法、系统及电子设备，用以解决现有技术中存在探测的水体分布情况的精度不高的技术问题。
本发明第一方面提供了一种核磁共振的超前探测方法，所述方法包括：
获得核磁共振装置的当前探测信号；
根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；
根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水层先验信息，获得所述掌子面前方区域的含水率分布。
本发明提供的核磁共振的超前探测方法，所述根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息，包括：
预先获取所述巷道的电阻率环境；
根据所述电阻率环境，采用核磁共振正演计算获得每个接收线圈的核函数分布，正演计算方法具体为：
其中，qi为激发脉冲矩，为发射线圈中电流持续时间t与电流幅值I0的乘积；为计算点位置；为位置处的网格体积；为激发线圈在某处产生的垂直场的正向分量；为接收线圈在某处产生的虚拟垂直场的反向分量；分别为激发和接收线圈在某处产生垂直磁场的相位； 分别为恒定磁场、激发场垂直分量、接收线圈虚拟磁场的垂直分量在某处的单位向量；
将所述每个接收线圈的核函数分布合并，获得合并结果，以所述合并结果作为所述核函数分布信息。
本发明提供的核磁共振的超前探测方法，在所述根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布之前，还包括：
获得预先探测得到的噪声信号；
将所述噪声信号作为预先探测信号的一部分。
本发明提供的核磁共振的超前探测方法，所述根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布，具体包括：
将所述核函数分布信息进行分解，得到分解量；
对所述分解量、所述预先的探测信号和预先获得的含水率分布进行反演，获得反演结果，以所述反演结果作为所述掌子面前方区域的含水率分布，其中反演公式具体为：
其中，mest为反演结果，K1、K2、K3、 K4为分解量，d1、d2为预先的探测信号，m1、m2为之前探测获得的含水率分布， d3为当前探测信号。
基于与第一方面相同的发明构思，本发明第二方面提供了一种核磁共振的超前探测系统，所述系统包括：
第一获得模块，用于获得核磁共振装置的当前探测信号；
计算模块，用于根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；
第二获得模块，用于根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布。
本发明提供的核磁共振的超前探测系统中，所述计算模块还用于：
预先获取所述巷道的电阻率环境；
根据所述电阻率环境，采用核磁共振正演计算获得每个接收线圈的核函数分布，正演计算方法具体为：
其中，qi为激发脉冲矩，为发射线圈中电流持续时间t与电流幅值I0的乘积；为计算点位置；为位置处的网格体积；为激发线圈在某处产生的垂直场的正向分量；为接收线圈在某处产生的虚拟垂直场的反向分量；分别为激发和接收线圈在某处产生垂直磁场的相位； 分别为恒定磁场、激发场垂直分量、接收线圈虚拟磁场的垂直分量在某处的单位向量；
将所述每个接收线圈的核函数分布合并，获得合并结果，以所述合并结果作为所述核函数分布信息。
本发明提供的核磁共振的超前探测系统中，所述系统还包括处理模块，所述处理模块用于：在所述根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布之前，
获得预先探测得到的噪声信号；
将所述噪声信号作为预先探测信号的一部分。
本发明提供的核磁共振的超前探测系统中，所述第二获得模块还用于：
将所述核函数分布信息进行分解，得到分解量；
对所述分解量、所述预先的探测信号和预先获得的含水率分布进行反演，获得反演结果，以所述反演结果作为所述掌子面前方区域的含水率分布，其中反演公式具体为：
其中，mest为反演结果，K1、K2、K3、 K4为分解量，d1、d2为预先的探测信号，m1、m2预先获得的含水率分布，d3为当前探测信号。
基于与第一方面相同的发明构思，本发明第三方面提供了一种电子设备，包括存储器，处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现以下步骤：
获得核磁共振装置的当前探测信号；
根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；
根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布。
本发明实施例中提供的一个或多个技术方案，至少具有如下技术效果或优点：
本发明公开的一种核磁共振的超前探测方法，通过获得核磁共振装置的当前探测信号；并根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；然后根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布。由于在进行核磁共振超前探测方法中，在计算巷道的含水率分布时，利用了预先获得的探测信号和预先获得的含水率分布，充分考虑了各次探测信号和探测结果之间的关联，并可以综合利用多次测量的重复信息，在不增加探测工作量的情况下，可以提升探测的精确性和有效性，从而得到精确的水体分布情况，解决了现有技术的探测方法存在探测的水体分布情况的精度不高的技术问题。
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的具体实施方式。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作一简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明实施例中一种核磁共振的超前探测方法的流程图；
图2为本发明实施例中一种核磁共振的超前探测系统的结构图；
图3为本发明实施例中核磁共振的超前探测方法的探测模式图；
图4为图3中探测模式中核函数的位置示意图；
图5为本发明提供的一种电子设备的结构图。
具体实施方式
本发明实施例提供一种核磁共振的超前探测方法、系统及电子设备，用以解决现有技术中存在探测的水体分布情况的精度不高的技术问题。
本申请实施例中的技术方案，总体思路如下：
一种核磁共振的超前探测方法，所述方法包括：首先获得核磁共振装置的当前探测信号；并根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；然后根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水层先验信息，获得所述掌子面前方区域的含水率分布。
上述方法中，由于在进行核磁共振超前探测方法中，在计算巷道的含水率分布时，利用了预先获得的探测信号和预先获得的含水率分布，充分考虑了各次探测信号和探测结果之间的关联，并可以综合利用多次测量的重复信息，在不增加探测工作量的情况下，可以提升探测的精确性和有效性，从而得到精确的水体分布情况，解决了现有技术的探测方法存在探测的水体分布情况的精度不高的技术问题。
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
实施例一
本实施例提供了一种核磁共振的超前探测方法，请参见图1，所述方法包括：
步骤S101：获得核磁共振装置的当前探测信号；
步骤S102：根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；
步骤S103：根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布。
上述方法中，由于在进行核磁共振超前探测方法中，在计算巷道的含水率分布时，利用了预先获得的探测信号和预先获得的含水率分布，充分考虑了各次探测信号和探测结果之间的关联，并可以综合利用多次测量的重复信息，在不增加探测工作量的情况下，可以提升探测的精确性和有效性，从而得到精确的水体分布情况，解决了现有技术的探测方法存在探测的水体分布情况的精度不高的技术问题。
需要说明的是上述步骤S101和步骤S102的执行顺序不分先后，可以是先执行步骤S101，也可以是先执行步骤S102。
下面，结合图1对本申请提供的核磁共振的超前探测方法进行详细介绍：
首先，执行步骤S101，获得核磁共振装置的当前探测信号。
在具体的实施过程中，可以采用核磁共振装置进行探测，核磁共振装置采用激发线圈和接收线圈相分离的方式，并且由于探测的空间多为狭窄的巷道或矿井，接收线圈采用较小的线圈。
然后执行步骤S102：根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息。
具体地，所述根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息，包括：
预先获取所述巷道的电阻率环境；
根据所述电阻率环境，采用核磁共振正演计算获得每个接收线圈的核函数分布，正演计算方法具体为：
其中，qi为激发脉冲矩，为发射线圈中电流持续时间t与电流幅值I0的乘积；为计算点位置；为位置处的网格体积；为激发线圈在某处产生的垂直场的正向分量；为接收线圈在某处产生的虚拟垂直场的反向分量；分别为激发和接收线圈在某处产生垂直磁场的相位； 分别为恒定磁场、激发场垂直分量、接收线圈虚拟磁场的垂直分量在某处的单位向量；
将所述每个接收线圈的核函数分布合并，获得合并结果，以所述合并结果作为所述核函数分布信息。
在具体的实现过程中，可以依据现场采样测量的方法或者直接利用前期的探测资料，确定巷道周围的电阻率，然后根据电阻率经核磁共振正演计算得到核函数分布；具体来说，为了得到不同位置的核函数分布，一般采用多个接收线圈，各个接收线圈分布于巷道的不同位置，上述表示的是当激发线圈的激发脉冲矩为qi时，位于地下位置的核函数。在实际计算时，可以将地下空间(即巷道)进行网格化，分布计算每个网格节点上的核函数，即代表着第j个网格节点的位置。等的具体计算方法为：
其中：
其中：上述以上参数均可由有限元算法求得。
最后执行步骤S103：根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布。
具体地，所述根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布，具体包括：
将所述核函数分布信息进行分解，得到分解量；
对所述分解量、所述预先的探测信号和预先获得的含水率分布进行反演，获得反演结果，以所述反演结果作为所述掌子面前方区域的含水率分布，其中反演公式具体为：
其中，mest为反演结果，K1、K2、K3、 K4为分解量，d1、d2为预先的探测信号，m1、m2为之前探测获得的含水率分布，d3为当前探测信号。
可选地，为了进一步提高探测结果的精度，在所述根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布之前，还包括：
获得预先探测得到的噪声信号；
将所述噪声信号作为预先探测信号的一部分。
在具体的实施过程中，可以采用垂直布置线圈的方式来获得噪声信号。
为了更详细地介绍本发明实施例提供的一种核磁共振的超前探测方法，下面通过具体的例子来进行说明。
请参见图3，掘进头的探测是分段进行的，掘进头前进的距离小于探测深度，而对掘进头获取的信息的计算时，计算区域的长度会大于探测深度，例如，本实施李忠掘进头的探测深度为75m，掘进头前进的距离为25m,计算区域的长度为100m,即掘进头每前进25m，实际计算的长度为100m，这样会导致每次掘进前，掘进头前方的待掘进区域将会被前期多次探测所覆盖到，而现有的方法中无法综合利用各次测量中被重复探测覆区域的信息，导致得到的水体分布情况的精度不高，并造成了勘探工作的冗余和浪费。在本实施例中，掌子面前方 0-25m的范围都会被3次测量所覆盖，25-50m的范围会被2次测量覆盖。
对于连续的三次观测，其核函数的计算区域都是以线圈前方100米的区域。如图4所示，该探测方案包括三个线圈：线圈1、线圈2和线圈3，以线圈3 为例，将其二维核函数在水平方向分解为K1至K4。对于线圈1至线圈3，由于当地电条件基本相似，则它们的核函数的大小及相对线圈的分布都相同，因此当线圈2激发时，K1至K4位于m2至m5位置；当线圈1激发时，K1至K4位于m1至m4位置，于是可以得到下述公式：
其中，d1为第一次的探测信号，m1-m4为线圈1得出的含水率分布，n1为第一次探测的噪声信号，d2为第二次的探测信号，m2-m5为线圈2得出的含水率分布，n2为第二次探测的噪声信号，d3为第三次的探测信号，即为当前次的探测信号，m3-m6为线圈3得出的含水率分布，n3为第三次探测的噪声信号。
当线圈3进行测量时，由于m1、m2区域已经开挖及勘测，前两次探测的结果为已知，所以上述公式(1)-(3)可整合为如下形式：
然后对核磁共振超前探进行反演，可以得到反演结果：
实施例二
基于与实施例一同样的发明构思，本发明实施例二提供了一种核磁共振的超前探测系统，请参见图2，所述系统包括：
第一获得模,201，用于获得核磁共振装置的当前探测信号；
计算模块202，用于根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；
第二获得模块203，用于根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布。
本发明实施例提供的核磁共振的超前探测系统中，所述计算模块202还用于：
预先获取所述巷道的电阻率环境；
根据所述电阻率环境，采用核磁共振正演计算获得每个接收线圈的核函数分布，正演计算方法具体为：
其中，qi为激发脉冲矩，为发射线圈中电流持续时间t与电流幅值I0的乘积；为计算点位置；为位置处的网格体积；为激发线圈在某处产生的垂直场的正向分量；为接收线圈在某处产生的虚拟垂直场的反向分量；分别为激发和接收线圈在某处产生垂直磁场的相位； 分别为恒定磁场、激发场垂直分量、接收线圈虚拟磁场的垂直分量在某处的单位向量；
将所述每个接收线圈的核函数分布合并，获得合并结果，以所述合并结果作为所述核函数分布信息。
本发明实施例提供的核磁共振的超前探测系统中，所述系统还包括处理模块，所述处理模块用于：在所述根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布之前，
获得预先探测得到的噪声信号；
将所述噪声信号作为预先探测信号的一部分。
本发明实施例提供的核磁共振的超前探测系统中，所述第二获得模块203 还用于：
将所述核函数分布信息进行分解，得到分解量；
对所述分解量、所述预先的探测信号和预先获得的含水率分布进行反演，获得反演结果，以所述反演结果作为所述掌子面前方区域的含水率分布，其中反演公式具体为：
其中，mest为反演结果，K1、K2、K3、 K4为分解量，d1、d2为预先的探测信号，m1、m2为之前探测获得的含水率分布， d3为当前探测信号。
实施例一中的核磁共振的超前探测方法的各种变化方式和具体实例同样适用于本实施例的核磁共振的超前探测系统，通过前述对核磁共振的超前探测方法的详细描述，本领域技术人员可以清楚的知道本实施例中中的核磁共振的超前探测系统，所以为了说明书的简洁，在此不再详述。
实施例三
基于与实施例一同样的发明构思，本发明实施例三提供了一种电子设备，请参见图5，包括存储器301，处理器302及存储在存储器301上并可在处理器302上运行的计算机程序，所述处理器302执行所述程序时实现以下步骤：
获得核磁共振装置的当前探测信号；
根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；
根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布。
为了便于说明，图5仅示出了与本发明实施例相关的部分，具体技术细节未揭示的，请参照本发明实施例方法部分。其中，存储器301可用于存储软件程序以及模块，处理器302通过运行执行存储在存储器301的软件程序以及模块，从而执行移动终端的各种功能应用以及数据处理。
存储器301可主要包括存储程序区和存储数据区，其中，存储程序区可存储操作系统、至少一个功能所需的应用程序等；存储数据区可存储根据电子设备的使用所创建的数据等。处理器302移动通信终端的控制中心，利用各种接口和线路连接整个移动通信终端的各个部分，通过运行或执行存储在存储器 301内的软件程序和/或模块，以及调用存储在存储器301内的数据，执行移动终端机的各种功能和处理数据，从而对移动终端机进行整体监控。可选的，处理器302可包括一个或多个处理单元。
实施例一中的核磁共振的超前探测方法的各种变化方式和具体实例同样适用于本实施例的电子设备，通过前述对核磁共振的超前探测方法的详细描述，本领域技术人员可以清楚的知道本实施例中的核磁共振的电子设备，所以为了说明书的简洁，在此不再详述。
本发明实施例中提供的一个或多个技术方案，至少具有如下技术效果或优点：
本发明公开的一种核磁共振的超前探测方法，通过获得核磁共振装置的当前探测信号；并根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；然后根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布。由于在进行核磁共振超前探测方法中，在计算巷道的含水率分布时，利用了预先获得的探测信号和预先获得的含水率分布，充分考虑了各次探测信号和探测结果之间的关联，并可以综合利用多次测量的重复信息，在不增加探测工作量的情况下，可以提升探测的精确性和有效性，从而得到精确的水体分布情况，解决了现有技术的探测方法存在探测的水体分布情况的精度不高的技术问题。
尽管已描述了本发明的优选实施例，但本领域内的技术人员一旦得知了基本创造性概念，则可对这些实施例作出另外的变更和修改。所以，所附权利要求意欲解释为包括优选实施例以及落入本发明范围的所有变更和修改。
显然，本领域的技术人员可以对本发明实施例进行各种改动和变型而不脱离本发明实施例的精神和范围。这样，倘若本发明实施例的这些修改和变型属于本发明权利要求及其等同技术的范围之内，则本发明也意图包含这些改动和变型在内。
Claims (9)
1.一种核磁共振的超前探测方法，其特征在于，所述方法包括：
获得核磁共振装置的当前探测信号；
根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；
根据所述当前探测信号、所述核函数分布信息、预先的探测信号和含水层先验信息，获得所述掌子面前方区域的含水率分布。
2.如权利要求1所述的方法，其特征在于，所述根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息，包括：
预先获取所述巷道的电阻率环境；
根据所述电阻率环境，采用核磁共振正演计算获得每个接收线圈的核函数分布，正演计算方法具体为：
其中，qi为激发脉冲矩，为发射线圈中电流持续时间t与电流幅值I0的乘积；为计算点位置；为位置处的网格体积；为激发线圈在某处产生的垂直场的正向分量；为接收线圈在某处产生的虚拟垂直场的反向分量；分别为激发和接收线圈在某处产生垂直磁场的相位； 分别为恒定磁场、激发场垂直分量、接收线圈虚拟磁场的垂直分量在某处的单位向量；
将所述每个接收线圈的核函数分布合并，获得合并结果，以所述合并结果作为所述核函数分布信息。
3.如权利要求1所述的方法，其特征在于，在所述根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水层先验信息，获得所述掌子面前方区域的含水率分布之前，还包括：
获得预先探测得到的噪声信号；
将所述噪声信号作为预先探测信号的一部分。
4.如权利要求1所述的方法，其特征在于，所述根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水层先验信息，获得所述掌子面前方区域的含水率分布，具体包括：
将所述核函数分布信息进行分解，得到分解量；
对所述分解量、所述预先的探测信号和预先获得的含水率分布进行反演，获得反演结果，以所述反演结果作为所述掌子面前方区域的含水率分布，其中反演公式具体为：
其中，mest为反演结果，K1、K2、K3、K4为分解量，d1、d2为预先的探测信号，m1、m2预先获得的含水率分布，d3为当前探测信号。
5.一种核磁共振的超前探测系统，其特征在于，所述系统包括：
第一获得模块，用于获得核磁共振装置的当前探测信号；
计算模块，用于根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；
第二获得模块，用于根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水率分布，获得所述掌子面前方区域的含水率分布。
6.如权利要求5所述的系统，其特征在于，所述计算模块还用于：
预先获取所述巷道的电阻率环境；
根据所述电阻率环境，采用核磁共振正演计算获得每个接收线圈的核函数分布，正演计算方法具体为：
其中，qi为激发脉冲矩，为发射线圈中电流持续时间t与电流幅值I0的乘积；为计算点位置；为位置处的网格体积；为激发线圈在某处产生的垂直场的正向分量；为接收线圈在某处产生的虚拟垂直场的反向分量；分别为激发和接收线圈在某处产生垂直磁场的相位； 分别为恒定磁场、激发场垂直分量、接收线圈虚拟磁场的垂直分量在某处的单位向量；
将所述每个接收线圈的核函数分布合并，获得合并结果，以所述合并结果作为所述核函数分布信息。
7.如权利要求5所述的系统，其特征在于，所述系统还包括处理模块，所述处理模块用于：在所述根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水层先验信息，获得所述掌子面前方区域的含水率分布之前，获得预先探测得到的噪声信号；
将所述噪声信号作为预先探测信号的一部分。
8.如权利要求5所述的系统，其特征在于，所述第二获得模块还用于：
将所述核函数分布信息进行分解，得到分解量；
对所述分解量、所述预先的探测信号和之前探测获得的含水率分布进行反演，获得反演结果，以所述反演结果作为所述掌子面前方区域的含水率分布，其中反演公式具体为：
其中，mest为反演结果，K1、K2、K3、K4为分解量，d1、d2为预先的探测信号，m1、m2为之前探测获得的含水率分布，d3为当前探测信号。
9.一种电子设备，包括存储器，处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现以下步骤：
获得核磁共振装置的当前探测信号；
根据核磁共振正演计算获得所述核磁共振装置的接收线圈的核函数分布信息；
根据所述当前探测信号、所述核函数分布信息、预先的探测信号和预先获得的含水层先验信息，获得所述掌子面前方区域的含水率分布。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710444249.5A CN107179556B (zh) | 2017-06-13 | 2017-06-13 | 一种核磁共振的超前探测方法、系统及电子设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710444249.5A CN107179556B (zh) | 2017-06-13 | 2017-06-13 | 一种核磁共振的超前探测方法、系统及电子设备 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN107179556A true CN107179556A (zh) | 2017-09-19 |
| CN107179556B CN107179556B (zh) | 2019-05-10 |
Family
ID=59836388
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710444249.5A Active CN107179556B (zh) | 2017-06-13 | 2017-06-13 | 一种核磁共振的超前探测方法、系统及电子设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN107179556B (zh) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110515131A (zh) * | 2019-08-13 | 2019-11-29 | 吉林大学 | 一种基于时间推移技术的地面磁共振反演方法 |
| CN114594123A (zh) * | 2022-03-07 | 2022-06-07 | 武汉普瑞通科技有限公司 | 一种地下水污染的核磁共振探测方法、装置以及处理设备 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120074934A1 (en) * | 2010-09-29 | 2012-03-29 | Derrick Green | Nmr measurements and methods of analyzing nmr data |
| WO2012103397A3 (en) * | 2011-01-28 | 2012-10-04 | Schlumberger Technology B.V. | Two dimensional nmr of diffusion and relaxation for material characterization |
| CN104199114A (zh) * | 2014-09-03 | 2014-12-10 | 长安大学 | 一种基于瞬变电磁与核磁共振联合解释的隧道超前预报技术 |
| CN106707351A (zh) * | 2017-01-11 | 2017-05-24 | 吉林大学 | 一种tbm施工隧道超前探测磁共振装置系统 |
| CN106772645A (zh) * | 2016-12-15 | 2017-05-31 | 中国石油大学(北京) | 基于一般先验信息约束的核磁共振数据反演方法和装置 |
-
2017
- 2017-06-13 CN CN201710444249.5A patent/CN107179556B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20120074934A1 (en) * | 2010-09-29 | 2012-03-29 | Derrick Green | Nmr measurements and methods of analyzing nmr data |
| WO2012103397A3 (en) * | 2011-01-28 | 2012-10-04 | Schlumberger Technology B.V. | Two dimensional nmr of diffusion and relaxation for material characterization |
| CN104199114A (zh) * | 2014-09-03 | 2014-12-10 | 长安大学 | 一种基于瞬变电磁与核磁共振联合解释的隧道超前预报技术 |
| CN106772645A (zh) * | 2016-12-15 | 2017-05-31 | 中国石油大学(北京) | 基于一般先验信息约束的核磁共振数据反演方法和装置 |
| CN106707351A (zh) * | 2017-01-11 | 2017-05-24 | 吉林大学 | 一种tbm施工隧道超前探测磁共振装置系统 |
Non-Patent Citations (1)
| Title |
|---|
| 张楠: "核磁共振探测隧道掌子面前方含水体方法及应用研究", 《中国优秀硕士学位论文全文数据库•工程科技Ⅱ辑》 * |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110515131A (zh) * | 2019-08-13 | 2019-11-29 | 吉林大学 | 一种基于时间推移技术的地面磁共振反演方法 |
| CN110515131B (zh) * | 2019-08-13 | 2021-04-20 | 吉林大学 | 一种基于时间推移技术的地面磁共振反演方法 |
| CN114594123A (zh) * | 2022-03-07 | 2022-06-07 | 武汉普瑞通科技有限公司 | 一种地下水污染的核磁共振探测方法、装置以及处理设备 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN107179556B (zh) | 2019-05-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CA2671434A1 (en) | Systems and methods for measuring sea-bed resistivity | |
| CN105259584B (zh) | 一种水域电法勘探系统 | |
| CN113341467A (zh) | 基于多插值方法的矿井瞬变电磁三维显示方法 | |
| CN103454678A (zh) | 一种地震切片等时性的确定方法及系统 | |
| CN109143378A (zh) | 一种在煤矿巷道内顺层超前探测含水构造的二次时差方法 | |
| Chen et al. | An overview of ahead geological detection technologies in tunnels | |
| CN105467460B (zh) | 电磁勘探方法和装置 | |
| CN110320576A (zh) | 超前地质预报方法、装置及存储介质 | |
| CN107179556B (zh) | 一种核磁共振的超前探测方法、系统及电子设备 | |
| CN106291698B (zh) | 地震相沉积相确定方法和装置 | |
| Li et al. | Calibration method of azimuth electromagnetic boundary detection logging while drilling based on air–seawater double-layer dielectric as reflection interface | |
| EP3943985A1 (en) | System and method for mapping and monitoring reservoirs by electromagnetic crosswell and optimizing production | |
| Xie et al. | Fracture identification and characteristics of carbonate underground gas storage: an example from the eastern area of Sulige gas field, ordos Basin, China | |
| CN108549109A (zh) | 一种基于核磁共振技术的地面-巷道超前探测方法 | |
| Zang et al. | Forward and inverse models of magnetically-susceptible grout in rock fracture grouting | |
| Li et al. | From Streamline to Pathline: Visualizing Particle Trajectories Under Changing Velocity Fields | |
| CN104533405B (zh) | 泥页岩储层有机碳含量的检测方法和装置 | |
| CN105301638A (zh) | 一种提取风化层底界面的方法和装置 | |
| Guoqiang et al. | Detecting shallow caverns in China using TEM | |
| CN115929292A (zh) | 基于成像测井和常规测井融合的测井方法、装置及设备 | |
| WO2024082793A1 (zh) | 预测原位构造裂缝延伸长度的方法及装置 | |
| CN109884714A (zh) | 一种可控源电磁测量方法、装置及其存储介质 | |
| Thiel et al. | Side fault mapping enabled by 2D transverse inversion on new deep directional resistivity measurements | |
| Li et al. | Subsurface Injection Monitoring in Complex Geologic Media Using Pathline, Source Cloud and Time Cloud | |
| CN107783198A (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 |