CN102053280A - 带有参考线圈的核磁共振地下水探测***及探测方法 - Google Patents

带有参考线圈的核磁共振地下水探测***及探测方法 Download PDF

Info

Publication number
CN102053280A
CN102053280A CN2010105374652A CN201010537465A CN102053280A CN 102053280 A CN102053280 A CN 102053280A CN 2010105374652 A CN2010105374652 A CN 2010105374652A CN 201010537465 A CN201010537465 A CN 201010537465A CN 102053280 A CN102053280 A CN 102053280A
Authority
CN
China
Prior art keywords
emission
signal
magnetic resonance
reference coil
nuclear magnetic
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
Application number
CN2010105374652A
Other languages
English (en)
Other versions
CN102053280B (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 CN2010105374652A priority Critical patent/CN102053280B/zh
Publication of CN102053280A publication Critical patent/CN102053280A/zh
Application granted granted Critical
Publication of CN102053280B publication Critical patent/CN102053280B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/30—Assessment of water resources

Landscapes

  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及一种带有参考线圈的核磁共振地下水探测***及探测方法。通过多路A/D采集单元同步采集发射/接收线圈中的核磁共振信号以及参考线圈中噪声信号的全波形数据，通过计算参考线圈采集的噪声信号与核磁共振信号的最大相关性，实现参考线圈最佳位置和数量的布设，在信号和噪声统计特性未知的情况下，采用变步长自适应算法，最大限度对消发射/接收线圈获得核磁共振信号中的噪声，实现多场源复杂强噪声干扰下核磁共振信号的提取，有效地解决了村庄附近以及城市周边地区核磁共振探测干扰多、多种干扰噪声数据难以分离的问题，提高了仪器抗干扰性能，为在村庄附近以及城市周边地区寻找地下水提供了一种可靠的探测装置和方法。

Description

带有参考线圈的核磁共振地下水探测***及探测方法
技术领域
本发明涉及一种用于寻找地下水的核磁共振探测***，尤其是带有参考线圈的核磁共振地下水探测***及探测方法。
背景技术
核磁共振(Nuclear Magnetic Resonance，简称NMR)地下水探测方法是一种直接勘查地下水的地球物理方法，与传统的地球物理勘查地下水的方法相比具有高分辨力、高效率、信息量丰富和解唯一性等优点，不打钻就可以确定出含水层的深度、厚度、含水率、含水层平均孔隙度，且完成一个核磁共振测深点的费用仅为一个水文地质勘探钻孔费用的1/10。
但是，在NMR地下水探测方法中，探测对象是地下原子核的自旋磁矩，其在地面产生的电磁场十分微小，导致仪器接收装置需要具有检测极微弱信号(nV级)的能力才能实现核磁共振探测信号的提取。然而，高灵敏度的接收装置极易受到电磁噪声的干扰。目前在我国，电网覆盖率已达到96.4％，除西藏等偏远地区由于地理条件限制未架电网以外，电网已遍及全国各地。所以，在村庄附近以及城市周边地区，由于电力线噪声(50Hz及其谐波)以及偶发的奇异噪声影响非常严重，NMR地下水探测方法的应用受到限制。特别是拉莫尔频率非常靠近工频噪声的谐波频点时，NMR探水仪很难甚至不能获得可靠的核磁共振信号。
US6177794公开了一种了利用宏观表现出的核磁共振现象来寻找存在核磁共振现象的地下液体矿产的新技术，运用一组相位可控的接收天线阵列以及井中的接收天线来比较精确地测定地下存在核磁共振现象的液体，该方法利用地面上和地下的线圈组合同时接收信号，用来对核磁共振信号进行良好接收，可以实现对地下石油，水等矿产资源的探测。CN01278229.7公开了一种核磁共振地下水层探测仪，包括信号检测器，其特征在于，该信号检测器的一输入端连接一第一开关，该第一开关的号端串接一变速电阻器，该变速电阻器的另端串接第二和第三两开关，第二开关的另端串接一限流电阻，并依序串接有整流器和发电机，该发电机的另端与信号检测器的另一端连接，其中该第二和第三两开关之间接有一第一电容器，该第一电容器的另端与信号检测器的另一端连接，该第一开关与变速电阻器之间接有一线圈，该线圈的另端与信号检测器的另一端连接，该信号检测器的两端之间有一电容器。该发明将核磁共振的方法应用于地下水层的探测中，设计了适用于寻找500米以上地下水的核磁共振探测仪器。CN1936621公开了一种“核磁共振与瞬变电磁联用仪及其方法”，是将核磁共振与瞬变电磁组合成一体的核磁共振与瞬变电磁联用仪及其方法。首先将联用仪选择在瞬变电磁工作模式下，在测线上铺设发射线圈和接收线圈，对测区内每一个测点进行测量，测量完毕后，对瞬变电磁数据进行初步的处理，找出低电阻率点，并对低电阻率异常测点标定；再将联用仪切换到核磁共振工作模式下，并以已标定的异常测点为中心铺设发射线圈进行核磁共振测量工作，以测得数据或图形与瞬变电磁所测的电阻率异常进行比较，用以判断瞬变电磁所测的电阻率异常目的层的真伪。用一套设备实现两种仪器的功能，减少了设备投资，发挥了两种仪器各自优点，提高了探测效率和精度。CN101285895公开了一种“线电源激发多道接收地面核磁共振方法和***”，包括以下步骤：1)利用发射机通过两个打入地下的发射电极向地下供入拉摩尔频率的交变电流，通过改变电极的距离和电流的大小产生不同的激发磁场；2)断开激发场，在被激发的氢质子旋进到正常平衡态过程中，用高灵敏探头组接收呈指数衰减规律的信号；3)将所述的信号通过信号传输线传送给接收机，再将所述接收机接收的信号发送给信号处理与成像***，进行信号的预处理、去噪、反演和成像等处理。该发明能够提供地下小体积的精细勘察；多道接收工作效率高，探测空间范围大，有利于三维反演成像等特点；线电源激发，施工灵活简单；多道接收工作效率高，通过改变发射电流实现探测目的。US7466128B2公开了一种多通道核磁共振采集装置和处理方法，用于实现在静态磁场的三维核磁共振成像，不受控于应用静态磁场梯度。推荐的应用是用地磁场对地下水的精确探测和定位。多线圈阵列均可以工作在发射和接收模式，相应的数据处理算法应用到产生3维核磁共振自旋数据。为处理这种多线圈数据来估计3维核磁共振自旋密度分布，采用至少2个发射与接收线圈采集数据。CN101251606公开了“一种弱信号检测仪器中有用信号频带内工频谐波干扰抑制电路”，是利用频率、相位、幅度全面跟踪有用信号频带内工频谐波造成干扰的补偿电路。是主微处理器控制信号传感器与阻抗匹配网络的接通与断开，阻抗匹配网络经前置放大器、宽带滤波器与程控放大器连接，输出信号包络检测器通过控制线和数据总线与从微处理器连接，程控放大器经锁相环连接DDS信号发生器。有效地抑制检测仪器中有用信号频带内的工频谐波干扰，实时地跟踪这种干扰的频率、相位和幅度，将其从接收信号中滤除，而不影响其它有用信号，不影响和这种干扰频率重叠的有用信号频率成分。上述发明的核磁共振地下水探测仪器及方法中的现有技术大部分是采用核磁共振技术实现对地下水的精确探测以及采用多通道方式用于提高仪器工作效率、探测精度及实现3维核磁共振反演成像，但对于复杂多场源强干扰下核磁共振地下水探测***及方法则没有涉及；而发明的工频谐波干扰抑制电路则只能解决有用信号频带内的某一个工频谐波的滤除，对于其他工频谐波和奇异噪声等电磁干扰则无能为力。
发明内容
本发明的目的就是针对上述现有技术的不足，提供一种带有参考线圈的核磁共振地下水探测***及探测方法，以解决村庄附近以及城市周边地区复杂多场源强干扰下核磁共振信号有效提取难题。
本发明的目的是通过以下技术方案实现的：
计算机1通过串口线与控制***2连接，计算机1通过控制线与发射电源3连接，控制***2通过信号线与发射控制驱动4连接，发射控制驱动4通过信号线经大功率H桥路7与发射/接收线圈9连接，控制***2通过控制线与接收控制单元5连接，接收控制单元5通过控制线与多路A/D采集单元16连接，多路A/D采集单元16通过数据线与计算机1连接，控制***2通过控制线与收发切换开关6、第一切换开关13、第二切换开关18……乃至第N切换开关22连接，发射电源3通过交互信号线与大功率H桥路7连接，H桥路7输出端与配谐电容8和发射/接收线圈9连接，发射/接收线圈9通过信号线与收发切换开关6、信号调理电路10、信号放大器11和多路A/D采集单元16连接，第一参考线圈12经第一切换开关13、第一调理电路14、第一放大器15与多路A/D采集单元16连接，第二参考线圈17经第二切换开关18、第二调理电路19、第二放大器20与多路A/D采集单元16连接，第N参考线圈21经第N切换开关22、第N调理电路23、第N放大器24与多路A/D采集单元16连接构成。
带有参考线圈的核磁共振地下水探测***的探测方法，包括以下顺序和步骤：
a、在测区内选择探测地点，以探测点为中心，铺设方形发射/接收线圈9，根据探测点附近电磁干扰源情况确定铺设至少一个参考线圈乃至N个参考线圈，参考线圈要尽可能靠近干扰源铺设，与发射/接收线圈9线框边界的距离至少要大于方形发射/接收线圈9的边长，以此确保参考线圈采集的是噪声而不是核磁共振信号；
b、开启带有参考线圈的核磁共振地下水探测***，发射控制驱动4控制大功率H桥路向发射/接收线圈9中供入频率为拉莫尔频率的交变电流，通过计算机1控制改变发射电源3电压输出的大小，来改变H桥路在发射/接收线圈9中输出电流的大小，即改变激发脉冲矩的大小，从而激发地下不同深度的水体；
c、控制***2控制断开激发场，控制收发切换开关6、第一切换开关13、第二切换开关18乃至第N切换开关22同步将发射/接收线圈9中核磁共振信号和第一参考线圈12、第二参考线圈17和第N参考线圈21中噪声信号的全波形数据通过多路A/D采集单元16采集；
d、将多路A/D采集单元16采集的核磁共振信号和噪声信号全波形数据上传至计算机1，计算参考线圈采集的噪声信号与发射/接收线圈9采集的核磁共振信号的最大相关性，最佳值是0.7以上，如果计算值小于此值，则重新布设参考线圈的位置和数量，重复步骤b和步骤c；
e、基于自适应噪声抵消原理，采用带有一个或N个参考线圈的自适应噪声抵消***，变步长自适应最小均方(LMS)算法进行噪声对消，实现核磁共振信号的有效提取；
f、将上述提取的核磁共振全波形信号进行特征参数提取，获得初始振幅、弛豫时间、相位、频率信息；
g、根据获得的初始振幅、弛豫时间、相位、频率信息进行正反演处理，估算出地下水储量和渗透率。
有益效果：采用参考线圈实现对探测点附近多种干扰噪声数据的采集，在信号和噪声统计特性未知的情况下，采用变步长自适应算法，自适应跟踪学习参考线圈中的噪声，基于相关抵消原理最大限度对消发射/接收线圈中的噪声，实现多场源复杂强噪声干扰下核磁共振信号的提取，提高了仪器抗干扰性能，有效地解决了村庄附近以及城市周边地区核磁共振探测干扰多、多种干扰噪声数据难以分离的问题，为在村庄附近以及城市周边地区寻找地下水提供了一种可靠的探测装置和方法。
附图说明
图1是带有参考线圈的核磁共振地下水探测***结构框图。
图2是带有参考线圈的核磁共振地下水探测***信号自适应噪声抵消***框图。
1计算机，2控制***，3发射电源，4发射控制驱动，5接收控制单元，6收发切换开关，7大功率H桥路，8配谐电容，9发射/接收线圈，10信号调理电路，11信号放大器，12第一参考线圈，13第一切换开关，14第一调理电路，15第一放大器，16多路A/D采集单元，17第二参考线圈，18第二切换开关，19第二调理电路，20第二放大器，21第N参考线圈，22第N切换开关，23第N调理电路，24第N放大器
具体实施方式
下面结合附图和实施例作进一步的详细说明：
带有参考线圈的核磁共振地下水探测***是由计算机1通过串口线与控制***2连接，计算机1通过控制线与发射电源3连接，控制***2通过信号线与发射控制驱动4连接，发射控制驱动4通过信号线经大功率H桥路7与发射/接收线圈9连接，控制***2通过控制线与接收控制单元5连接，接收控制单元5通过控制线与多路A/D采集单元16连接，多路A/D采集单元16通过数据线与计算机1连接，控制***2通过控制线与收发切换开关6、第一切换开关13、第二切换开关18……乃至第N切换开关22连接，发射电源3通过交互信号线与大功率H桥路7连接，H桥路7输出端与配谐电容8和发射/接收线圈9连接，发射/接收线圈9通过信号线与收发切换开关6、信号调理电路10、信号放大器11和多路A/D采集单元16连接，第一参考线圈12经第一切换开关13、第一调理电路14、第一放大器15与多路A/D采集单元16连接，第二参考线圈17经第二切换开关18、第二调理电路19、第二放大器20与多路A/D采集单元16连接，第N参考线圈21经第N切换开关22、第N调理电路23、第N放大器24与多路A/D采集单元16连接构成。
具体工作过程是：
在发射阶段：计算机1控制发射电源3产生某一数值电压，并通过控制***2，控制发射控制驱动4，使之驱动大功率H桥路7，大功率H桥路7的两个输出端与配谐电容8和发射/接收线圈9连接构成串联谐振回路，使线圈中产生具有某一拉莫尔频率的交变电流，激发地下某一深度和范围内的地下水中的氢质子。激发深度取决于激发脉冲矩的大小，即发射/接收线圈中的电流大小。激发范围取决于发射/接收线圈9的大小。收发切换开关6采用高压干簧继电器，在发射阶段使继电器断开，使发射/接收线圈9不与接收***相连，第一切换开关13、第二切换开关18乃至第N切换开关均置于断开状态。
在接收阶段：控制***2同步控制收发切换开关6与第一切换开关13、第二切换开关18乃至第N切换开关22，使多路A/D采集单元16同步采集发射/接收线圈9中的核磁共振信号以及第一参考线圈12、第二参考线圈17…第N参考线圈21中的噪声信号。鉴于被测核磁共振信号的频率范围1.3-2.6KHz以及频谱特征，加之信号及其微弱，信号调理电路10、信号放大器11被采用，第一参考线圈12、第二参考线圈17…第N参考线圈21采集的噪声数据亦采用了第一调理电路14、第二调理电路19…第N调理电路23、第一放大器15、第二放大器20…第N放大器24被送至多路A/D采集单元16。
带有参考线圈的核磁共振地下水探测方法，具体实施步骤包括：
a、在测区内选择探测地点，以探测点为中心，铺设方形发射/接收线圈9，根据探测点附近电磁干扰源情况确定铺设至少一个参考线圈乃至N个参考线圈，参考线圈要尽可能靠近干扰源铺设，与发射/接收线圈9线框边界的距离至少要大于方形发射/接收线圈9的边长，以此确保参考线圈采集的是噪声而不是核磁共振信号；
b、开启带有参考线圈的核磁共振地下水探测***，发射控制驱动4控制大功率H桥路向发射/接收线圈9中供入频率为拉莫尔频率的交变电流，通过计算机1控制改变发射电源3电压输出的大小，来改变H桥路在发射/接收线圈9中输出电流的大小，即改变激发脉冲矩的大小，从而激发地下不同深度的水体；
c、控制***2控制断开激发场，控制收发切换开关6、第一切换开关13、第二切换开关18乃至第N切换开关22同步将发射/接收线圈9中核磁共振信号和第一参考线圈12、第二参考线圈17和第N参考线圈21中噪声信号的全波形数据通过多路A/D采集单元16采集；
d、将多路A/D采集单元16采集的核磁共振信号和噪声信号全波形数据上传至计算机1，计算参考线圈采集的噪声信号与发射/接收线圈9采集的核磁共振信号的最大相关性，最佳值是0.7以上，如果计算值小于此值，则重新布设参考线圈的位置和数量，重复步骤b和步骤c；
e、基于自适应噪声抵消原理，采用带有一个或N个参考线圈的自适应噪声抵消***，变步长自适应LMS算法进行噪声对消，实现核磁共振信号的有效提取；
f、将上述提取的核磁共振全波形信号进行特征参数提取，获得初始振幅、弛豫时间、相位、频率信息；
g、根据获得的初始振幅、弛豫时间、相位、频率信息进行正反演处理，估算出地下水储量和渗透率。
结合附图2，本发明设计多参考线圈方式实现对采集到的核磁共振信号的有效提取。在核磁共振信号中常见的对信号影响较大的噪声主要是工频谐波干扰以及奇异噪声，其数学模型表达形式均可以直接或者间接表示成正余弦信号，基于这一特点，采用90°移相方式相当于设计了多个自适应陷波器来进行自适应噪声对消。90°移相方式是通过Hilbert变换实现的。
所采用的变步长自适应LMS算法的迭代公式：
zk(n)＝Wk1(n)vk(n)+Wk2(n)vk⊥(n)
e(n)＝s(n)+v(n)-[z1(n)+...+zN(n)]
Wk1(n+1)＝Wk1(n)+2μ(n)e(n)vk(n)
Wk2(n+1)＝Wk2(n)+2μ(n)e(n)vk⊥(n)
(k＝1...N)
其中，vk(n)和vk⊥(n)代表第k个参考线圈采集到的噪声数据以及经过90°移相后的噪声数据。μ(n)表示变步长，其取决于野外测点的信噪比，当被测区域信噪比较低时，μ(n)取值要大一些，以便尽快跟踪噪声，相反，当被测区域信噪比较高时，μ(n)取值要小一些。
实施例1
带有参考线圈的核磁共振地下水探测***是由计算机1通过串口线与控制***2连接，计算机1通过控制线与发射电源3连接，控制***2通过信号线与发射控制驱动4连接，发射控制驱动4通过信号线经大功率H桥路7与发射/接收线圈9连接，控制***2通过控制线与接收控制单元5连接，接收控制单元5通过控制线与多路A/D采集单元16连接，多路A/D采集单元16通过数据线与计算机1连接，控制***2通过控制线与收发切换开关6、第一切换开关13连接，发射电源3通过交互信号线与大功率H桥路7连接，H桥路7输出端与配谐电容8和发射/接收线圈9连接，发射/接收线圈9通过信号线与收发切换开关6、信号调理电路10、信号放大器11和多路A/D采集单元16连接，第一参考线圈12经第一切换开关13、第一调理电路14、第一放大器15与多路A/D采集单元16连接构成。
带有参考线圈的核磁共振地下水探测方法，具体实施步骤包括：
a、在测区内选择探测地点，以探测点为中心，铺设方形发射/接收线圈9，根据探测点附近电磁干扰源情况确定铺设一个参考线圈，参考线圈要尽可能靠近干扰源铺设，与发射/接收线圈9线框边界的距离至少要大于方形发射/接收线圈9的边长，以此确保参考线圈采集的是噪声而不是核磁共振信号；
b、开启带有参考线圈的核磁共振地下水探测***，发射控制驱动4控制大功率H桥路向发射/接收线圈9中供入频率为拉莫尔频率的交变电流，通过计算机1控制改变发射电源3电压输出的大小，来改变H桥路在发射/接收线圈9中输出电流的大小，即改变激发脉冲矩的大小，从而激发地下不同深度的水体；
c、控制***2控制断开激发场，控制收发切换开关6、第一切换开关13同步将发射/接收线圈9中核磁共振信号和第一参考线圈12中噪声信号的全波形数据通过多路A/D采集单元16采集；
d、将多路A/D采集单元16采集的核磁共振信号和噪声信号全波形数据上传至计算机1，计算参考线圈采集的噪声信号与发射/接收线圈9采集的核磁共振信号的最大相关性，最佳值是0.7以上，如果计算值小于此值，则重新布设参考线圈的位置，重复步骤b和步骤c；
e、基于自适应噪声抵消原理，采用带有一个参考线圈的自适应噪声抵消***，变步长自适应LMS算法进行噪声对消，实现核磁共振信号的有效提取；
f、将上述提取的核磁共振全波形信号进行特征参数提取，获得初始振幅、弛豫时间、相位、频率信息；
g、根据获得的初始振幅、弛豫时间、相位、频率信息进行正反演处理，估算出地下水储量和渗透率。
实施例2
计算机1通过串口线与控制***2连接，计算机1通过控制线与发射电源3连接，控制***2通过信号线与发射控制驱动4连接，发射控制驱动4通过信号线经大功率H桥路7与发射/接收线圈9连接，控制***2通过控制线与接收控制单元5连接，接收控制单元5通过控制线与多路A/D采集单元16连接，多路A/D采集单元16通过数据线与计算机1连接，控制***2通过控制线与收发切换开关6、第一切换开关13、第二切换开关18连接，发射电源3通过交互信号线与大功率H桥路7连接，H桥路7输出端与配谐电容8和发射/接收线圈9连接，发射/接收线圈9通过信号线与收发切换开关6、信号调理电路10、信号放大器11和多路A/D采集单元16连接，第一参考线圈12经第一切换开关13、第一调理电路14、第一放大器15与多路A/D采集单元16连接，第二参考线圈17经第二切换开关18、第二调理电路19、第二放大器20与多路A/D采集单元16连接构成。
带有参考线圈的核磁共振地下水探测方法，具体实施步骤包括：
a、在测区内选择探测地点，以探测点为中心，铺设方形发射/接收线圈9，根据探测点附近电磁干扰源情况确定铺设二个参考线圈，参考线圈要尽可能靠近干扰源铺设，与发射/接收线圈9线框边界的距离至少要大于方形发射/接收线圈9的边长，以此确保参考线圈采集的是噪声而不是核磁共振信号；
b、开启带有参考线圈的核磁共振地下水探测***，发射控制驱动4控制大功率H桥路向发射/接收线圈9中供入频率为拉莫尔频率的交变电流，通过计算机1控制改变发射电源3电压输出的大小，来改变H桥路在发射/接收线圈9中输出电流的大小，即改变激发脉冲矩的大小，从而激发地下不同深度的水体；
c、控制***2控制断开激发场，控制收发切换开关6、第一切换开关13、第二切换开关18同步将发射/接收线圈9中核磁共振信号和第一参考线圈12、第二参考线圈17中噪声信号的全波形数据通过多路A/D采集单元16采集；
d、将多路A/D采集单元16采集的核磁共振信号和噪声信号全波形数据上传至计算机1，计算参考线圈采集的噪声信号与发射/接收线圈9采集的核磁共振信号的最大相关性，最佳值是0.7以上，如果计算值小于此值，则重新布设参考线圈的位置，重复步骤b和步骤c；
e、基于自适应噪声抵消原理，采用带有二个参考线圈的自适应噪声抵消***，变步长自适应LMS算法进行噪声对消，实现核磁共振信号的有效提取；
f、将上述提取的核磁共振全波形信号进行特征参数提取，获得初始振幅、弛豫时间、相位、频率信息；
g、根据获得的初始振幅、弛豫时间、相位、频率信息进行正反演处理，估算出地下水储量和渗透率。
实施例3
计算机1通过串口线与控制***2连接，计算机1通过控制线与发射电源3连接，控制***2通过信号线与发射控制驱动4连接，发射控制驱动4通过信号线经大功率H桥路7与发射/接收线圈9连接，控制***2通过控制线与接收控制单元5连接，接收控制单元5通过控制线与多路A/D采集单元16连接，多路A/D采集单元16通过数据线与计算机1连接，控制***2通过控制线与收发切换开关6、第一切换开关13、第二切换开关18……乃至第八切换开关22连接，发射电源3通过交互信号线与大功率H桥路7连接，H桥路7输出端与配谐电容8和发射/接收线圈9连接，发射/接收线圈9通过信号线与收发切换开关6、信号调理电路10、信号放大器11和多路A/D采集单元16连接，第一参考线圈12经第一切换开关13、第一调理电路14、第一放大器15与多路A/D采集单元16连接，第二参考线圈17经第二切换开关18、第二调理电路19、第二放大器20与多路A/D采集单元16连接，第八参考线圈21经第八切换开关22、第八调理电路23、第八放大器24与多路A/D采集单元16连接构成。
带有参考线圈的核磁共振地下水探测方法，具体实施步骤包括：
a、在测区内选择探测地点，以探测点为中心，铺设方形发射/接收线圈9，根据探测点附近电磁干扰源情况确定铺设八个参考线圈，参考线圈要尽可能靠近干扰源铺设，与发射/接收线圈9线框边界的距离至少要大于方形发射/接收线圈9的边长，以此确保参考线圈采集的是噪声而不是核磁共振信号；
b、开启带有参考线圈的核磁共振地下水探测***，发射控制驱动4控制大功率H桥路向发射/接收线圈9中供入频率为拉莫尔频率的交变电流，通过计算机1控制改变发射电源3电压输出的大小，来改变H桥路在发射/接收线圈9中输出电流的大小，即改变激发脉冲矩的大小，从而激发地下不同深度的水体；
c、控制***2控制断开激发场，控制收发切换开关6、第一切换开关13、第二切换开关18乃至第八切换开关22同步将发射/接收线圈9中核磁共振信号和第一参考线圈12、第二参考线圈17和第八参考线圈21中噪声信号的全波形数据通过多路A/D采集单元16采集；
d、将多路A/D采集单元16采集的核磁共振信号和噪声信号全波形数据上传至计算机1，计算参考线圈采集的噪声信号与发射/接收线圈9采集的核磁共振信号的最大相关性，最佳值是0.7以上，如果计算值小于此值，则重新布设参考线圈的位置，重复步骤b和步骤c；
e、基于自适应噪声抵消原理，采用带有八个参考线圈的自适应噪声抵消***，变步长自适应LMS算法进行噪声对消，实现核磁共振信号的有效提取；
f、将上述提取的核磁共振全波形信号进行特征参数提取，获得初始振幅、弛豫时间、相位、频率信息；
g、根据获得的初始振幅、弛豫时间、相位、频率信息进行正反演处理，估算出地下水储量和渗透率。

Claims (2)

1.一种带有参考线圈的核磁共振地下水探测***，其特征在于，计算机1通过串口线与控制***2连接，计算机1通过控制线与发射电源3连接，控制***2通过信号线与发射控制驱动4连接，发射控制驱动4通过信号线经大功率H桥路7与发射/接收线圈9连接，控制***2通过控制线与接收控制单元5连接，接收控制单元5通过控制线与多路A/D采集单元16连接，多路A/D采集单元16通过数据线与计算机1连接，控制***2通过控制线与收发切换开关6、第一切换开关13、第二切换开关18……乃至第N切换开关22连接，发射电源3通过交互信号线与大功率H桥路7连接，H桥路7输出端与配谐电容8和发射/接收线圈9连接，发射/接收线圈9通过信号线与收发切换开关6、信号调理电路10、信号放大器11和多路A/D采集单元16连接，第一参考线圈12经第一切换开关13、第一调理电路14、第一放大器15与多路A/D采集单元16连接，第二参考线圈17经第二切换开关18、第二调理电路19、第二放大器20与多路A/D采集单元16连接，第N参考线圈21经第N切换开关22、第N调理电路23、第N放大器24与多路A/D采集单元16连接构成。
2.按照权利要求1所述的带有参考线圈的核磁共振地下水探测***，其特征在于，按以下顺序和步骤工作：
a、在测区内选择探测地点，以探测点为中心，铺设方形发射/接收线圈9，根据探测点附近电磁干扰源情况确定铺设至少一个参考线圈乃至N个参考线圈，参考线圈要尽可能靠近干扰源铺设，与发射/接收线圈9线框边界的距离至少要大于方形发射/接收线圈9的边长，以此确保参考线圈采集的是噪声而不是核磁共振信号；
b、开启带有参考线圈的核磁共振地下水探测***，发射控制驱动4控制大功率H桥路向发射/接收线圈9中供入频率为拉莫尔频率的交变电流，通过计算机1控制改变发射电源3电压输出的大小，来改变H桥路在发射/接收线圈9中输出电流的大小，即改变激发脉冲矩的大小，从而激发地下不同深度的水体；
c、控制***2控制断开激发场，控制收发切换开关6、第一切换开关13、第二切换开关18乃至第N切换开关22同步将发射/接收线圈9中核磁共振信号和第一参考线圈12、第二参考线圈17和第N参考线圈21中噪声信号的全波形数据通过多路A/D采集单元16采集；
d、将多路A/D采集单元16采集的核磁共振信号和噪声信号全波形数据上传至计算机1，计算参考线圈采集的噪声信号与发射/接收线圈9采集的核磁共振信号的最大相关性，最佳值是0.7以上，如果计算值小于此值，则重新布设参考线圈的位置和数量，重复步骤b和步骤c；
e、基于自适应噪声抵消原理，采用带有一个或N个参考线圈的自适应噪声抵消***，变步长自适应最小均方算法进行噪声对消，实现核磁共振信号的有效提取；
f、将上述提取的核磁共振全波形信号进行特征参数提取，获得初始振幅、弛豫时间、相位、频率信息；
g、根据获得的初始振幅、弛豫时间、相位、频率信息进行正反演处理，估算出地下水储量和渗透率。
CN2010105374652A 2010-11-10 2010-11-10 带有参考线圈的核磁共振地下水探测***及探测方法 Expired - Fee Related CN102053280B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2010105374652A CN102053280B (zh) 2010-11-10 2010-11-10 带有参考线圈的核磁共振地下水探测***及探测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2010105374652A CN102053280B (zh) 2010-11-10 2010-11-10 带有参考线圈的核磁共振地下水探测***及探测方法

Publications (2)

Publication Number Publication Date
CN102053280A true CN102053280A (zh) 2011-05-11
CN102053280B CN102053280B (zh) 2013-05-01

Family

ID=43957805

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2010105374652A Expired - Fee Related CN102053280B (zh) 2010-11-10 2010-11-10 带有参考线圈的核磁共振地下水探测***及探测方法

Country Status (1)

Country Link
CN (1) CN102053280B (zh)

Cited By (37)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102221711A (zh) * 2011-05-23 2011-10-19 吉林大学 核磁共振差分探测坑道突水超前预测装置及探测方法
CN103033849A (zh) * 2012-12-15 2013-04-10 吉林大学 多通道核磁共振地下水探测仪及其野外工作方法
CN103344995A (zh) * 2013-06-24 2013-10-09 吉林大学 引入人工磁场的核磁共振定向探测装置及探测方法
CN103412341A (zh) * 2013-08-22 2013-11-27 吉林大学 冷线圈核磁共振地下水探测装置及探测方法
CN103823244A (zh) * 2014-03-11 2014-05-28 吉林大学 磁共振三分量消噪装置及消噪方法
CN103913778A (zh) * 2014-04-10 2014-07-09 吉林大学 多个近端参考线圈的核磁共振信号实时噪声抵消装置
CN103969612A (zh) * 2013-01-28 2014-08-06 北京大学 一种用于核磁共振成像***的多功能采集模块
CN103984032A (zh) * 2014-04-28 2014-08-13 吉林大学 实时信噪特征提取的多工作模式核磁信号探测装置及方法
CN104407392A (zh) * 2014-12-09 2015-03-11 吉林大学 一发三收式对充水采空区的探测装置及探测方法
CN104407391A (zh) * 2014-12-05 2015-03-11 吉林大学 磁性源非调制式发射机及控制方法
CN104502984A (zh) * 2014-12-15 2015-04-08 吉林大学 特定频率噪声对消地下核磁共振探测装置及探测方法
CN104777442A (zh) * 2015-04-07 2015-07-15 吉林大学 一种核磁共振测深fid信号噪声抑制方法
CN104796110A (zh) * 2015-05-05 2015-07-22 吉林大学 模拟工频梳状陷波器及其调整方法
CN105223622A (zh) * 2015-09-24 2016-01-06 吉林大学 抗饱和高增益自适应抵消型地面核磁共振信号放大装置
CN105652328A (zh) * 2016-03-09 2016-06-08 吉林大学 核磁共振探水仪的自动配谐装置
CN105824053A (zh) * 2016-05-23 2016-08-03 吉林大学 自适应滤波的磁共振信号抗饱和消噪装置及消噪方法
CN106291724A (zh) * 2016-10-08 2017-01-04 重庆大学 一种用于地下核磁共振探水的发射/接收线圈
CN104216021B (zh) * 2014-09-09 2017-03-22 吉林大学 一种基于分步式发射的地下核磁共振探测方法
CN106772161A (zh) * 2017-02-28 2017-05-31 吉林大学 无线传感网络的多通道阵列式接收线圈及探测方法
CN106814400A (zh) * 2017-03-23 2017-06-09 吉林大学 基于阵列逆变充电的核磁共振探水发射装置及工作方法
CN107102280A (zh) * 2017-06-13 2017-08-29 武汉普瑞通科技有限公司 一种核磁共振信号接收线圈、装置以及核磁共振超前探测方法
CN107942397A (zh) * 2017-12-29 2018-04-20 吉林大学 用预极化场增强信号幅度的磁共振多通道探测方法及装置
CN108919365A (zh) * 2018-05-29 2018-11-30 福州华虹智能科技股份有限公司 一种无线电坑道***
CN109031435A (zh) * 2018-08-13 2018-12-18 吉林大学 基于相控阵的核磁共振地下水定向探测***及探测方法
CN109358369A (zh) * 2018-12-19 2019-02-19 冯家森 一种立体线圈式抗电磁干扰脉冲金属探测器
CN110109399A (zh) * 2019-07-01 2019-08-09 南京理工大学 基于fpga与stm32的多通道磁信号采集存储装置
CN110456420A (zh) * 2019-08-13 2019-11-15 吉林大学 基于近端参考线圈核磁共振地下水探测信号噪声消除方法
CN111290037A (zh) * 2020-03-13 2020-06-16 吉林大学 磁共振无损原位探测地下有机污染物的装置及探测方法
CN111796331A (zh) * 2020-08-24 2020-10-20 吉林大学 面向浅层地下水和烃类物质的地面磁共振探测装置及方法
CN113376704A (zh) * 2021-06-07 2021-09-10 电子科技大学 一种基于电发射-磁接收的井间电磁探测***及方法
CN113866837A (zh) * 2021-09-18 2021-12-31 吉林大学 电性源核磁共振与激发极化联用装置以及探测方法
CN115097533A (zh) * 2022-05-05 2022-09-23 吉林大学 一种基于tls-esprit算法的磁共振测深信号提取方法
CN115808648A (zh) * 2022-11-18 2023-03-17 无锡鸣石峻致医疗科技有限公司 一种磁共振***振铃噪声的测量装置和方法
CN116449438A (zh) * 2023-06-16 2023-07-18 中国地质大学（武汉） 一种瞬变电磁法参考消噪装置及方法
CN116856920A (zh) * 2023-07-06 2023-10-10 中国科学院地质与地球物理研究所 一种随钻方位电磁波电阻率仪器使用方法及仪器
CN117347416A (zh) * 2023-10-12 2024-01-05 无锡鸣石峻致医疗科技有限公司 一种物质成分的非侵入式检测方法、装置及存储介质
CN118067770A (zh) * 2024-04-18 2024-05-24 吉林大学 一种基于低场磁共振技术的粮仓水分监测装置与方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5578920A (en) * 1994-07-14 1996-11-26 Spectrospin Ag Compensation of interfering fields in NMR measurements in the earth's magnetic field
CN200997000Y (zh) * 2006-09-29 2007-12-26 吉林大学 地面核磁共振找水仪发射装置
CN101251606A (zh) * 2008-04-09 2008-08-27 吉林大学 弱信号检测仪器中有用信号频带内工频谐波干扰抑制电路
CN101285895A (zh) * 2008-06-06 2008-10-15 李振宇 线电源激发多道接收地面核磁共振方法和***
CN201278026Y (zh) * 2008-10-30 2009-07-22 吉林大学 地面核磁共振信号源

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US5578920A (en) * 1994-07-14 1996-11-26 Spectrospin Ag Compensation of interfering fields in NMR measurements in the earth's magnetic field
CN200997000Y (zh) * 2006-09-29 2007-12-26 吉林大学 地面核磁共振找水仪发射装置
CN101251606A (zh) * 2008-04-09 2008-08-27 吉林大学 弱信号检测仪器中有用信号频带内工频谐波干扰抑制电路
CN101285895A (zh) * 2008-06-06 2008-10-15 李振宇 线电源激发多道接收地面核磁共振方法和***
CN201278026Y (zh) * 2008-10-30 2009-07-22 吉林大学 地面核磁共振信号源

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
田宝凤,段清明: "核磁共振信号工频谐波的自适应滤除方法", 《吉林大学学报（信息科学版）》 *

Cited By (59)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US9638779B2 (en) 2011-05-23 2017-05-02 Jilin University Equipment and method forecasting tunnel water inrush using magnetic resonance differential
WO2012159546A1 (zh) * 2011-05-23 2012-11-29 吉林大学 利用核磁共振差分超前预测坑道突水的装置及方法
CN102221711B (zh) * 2011-05-23 2013-05-01 吉林大学 核磁共振差分探测坑道突水超前预测装置及探测方法
CN102221711A (zh) * 2011-05-23 2011-10-19 吉林大学 核磁共振差分探测坑道突水超前预测装置及探测方法
CN103033849A (zh) * 2012-12-15 2013-04-10 吉林大学 多通道核磁共振地下水探测仪及其野外工作方法
CN103033849B (zh) * 2012-12-15 2015-06-10 吉林大学 多通道核磁共振地下水探测仪及其野外工作方法
CN103969612A (zh) * 2013-01-28 2014-08-06 北京大学 一种用于核磁共振成像***的多功能采集模块
CN103344995A (zh) * 2013-06-24 2013-10-09 吉林大学 引入人工磁场的核磁共振定向探测装置及探测方法
CN103344995B (zh) * 2013-06-24 2016-03-02 吉林大学 引入人工磁场的核磁共振定向探测装置的探测方法
CN103412341A (zh) * 2013-08-22 2013-11-27 吉林大学 冷线圈核磁共振地下水探测装置及探测方法
CN103412341B (zh) * 2013-08-22 2015-10-14 吉林大学 冷线圈核磁共振地下水探测装置及探测方法
CN103823244A (zh) * 2014-03-11 2014-05-28 吉林大学 磁共振三分量消噪装置及消噪方法
CN103823244B (zh) * 2014-03-11 2015-04-22 吉林大学 磁共振三分量消噪装置及消噪方法
CN103913778A (zh) * 2014-04-10 2014-07-09 吉林大学 多个近端参考线圈的核磁共振信号实时噪声抵消装置
CN103984032A (zh) * 2014-04-28 2014-08-13 吉林大学 实时信噪特征提取的多工作模式核磁信号探测装置及方法
CN103984032B (zh) * 2014-04-28 2017-02-15 吉林大学 实时信噪特征提取的多工作模式核磁信号探测装置及方法
CN104216021B (zh) * 2014-09-09 2017-03-22 吉林大学 一种基于分步式发射的地下核磁共振探测方法
CN104407391A (zh) * 2014-12-05 2015-03-11 吉林大学 磁性源非调制式发射机及控制方法
CN104407391B (zh) * 2014-12-05 2017-02-22 吉林大学 磁性源非调制式发射机及控制方法
CN104407392A (zh) * 2014-12-09 2015-03-11 吉林大学 一发三收式对充水采空区的探测装置及探测方法
CN104502984A (zh) * 2014-12-15 2015-04-08 吉林大学 特定频率噪声对消地下核磁共振探测装置及探测方法
CN104777442A (zh) * 2015-04-07 2015-07-15 吉林大学 一种核磁共振测深fid信号噪声抑制方法
CN104777442B (zh) * 2015-04-07 2017-06-16 吉林大学 一种核磁共振测深fid信号噪声抑制方法
CN104796110A (zh) * 2015-05-05 2015-07-22 吉林大学 模拟工频梳状陷波器及其调整方法
CN104796110B (zh) * 2015-05-05 2017-11-28 吉林大学 模拟工频梳状陷波器及其调整方法
CN105223622A (zh) * 2015-09-24 2016-01-06 吉林大学 抗饱和高增益自适应抵消型地面核磁共振信号放大装置
CN105652328A (zh) * 2016-03-09 2016-06-08 吉林大学 核磁共振探水仪的自动配谐装置
CN105824053A (zh) * 2016-05-23 2016-08-03 吉林大学 自适应滤波的磁共振信号抗饱和消噪装置及消噪方法
CN105824053B (zh) * 2016-05-23 2018-09-28 吉林大学 自适应滤波的磁共振信号抗饱和消噪装置及消噪方法
CN106291724A (zh) * 2016-10-08 2017-01-04 重庆大学 一种用于地下核磁共振探水的发射/接收线圈
CN106772161A (zh) * 2017-02-28 2017-05-31 吉林大学 无线传感网络的多通道阵列式接收线圈及探测方法
CN106772161B (zh) * 2017-02-28 2019-04-30 吉林大学 无线传感网络的多通道阵列式接收线圈及探测方法
CN106814400A (zh) * 2017-03-23 2017-06-09 吉林大学 基于阵列逆变充电的核磁共振探水发射装置及工作方法
CN107102280A (zh) * 2017-06-13 2017-08-29 武汉普瑞通科技有限公司 一种核磁共振信号接收线圈、装置以及核磁共振超前探测方法
CN107942397A (zh) * 2017-12-29 2018-04-20 吉林大学 用预极化场增强信号幅度的磁共振多通道探测方法及装置
CN108919365A (zh) * 2018-05-29 2018-11-30 福州华虹智能科技股份有限公司 一种无线电坑道***
CN109031435A (zh) * 2018-08-13 2018-12-18 吉林大学 基于相控阵的核磁共振地下水定向探测***及探测方法
CN109031435B (zh) * 2018-08-13 2023-06-30 吉林大学 基于相控阵的核磁共振地下水定向探测***及探测方法
CN109358369A (zh) * 2018-12-19 2019-02-19 冯家森 一种立体线圈式抗电磁干扰脉冲金属探测器
CN109358369B (zh) * 2018-12-19 2024-04-26 冯家森 一种立体线圈式抗电磁干扰脉冲金属探测器
CN110109399A (zh) * 2019-07-01 2019-08-09 南京理工大学 基于fpga与stm32的多通道磁信号采集存储装置
CN110456420A (zh) * 2019-08-13 2019-11-15 吉林大学 基于近端参考线圈核磁共振地下水探测信号噪声消除方法
CN111290037B (zh) * 2020-03-13 2021-05-04 吉林大学 磁共振无损原位探测地下有机污染物的装置及探测方法
CN111290037A (zh) * 2020-03-13 2020-06-16 吉林大学 磁共振无损原位探测地下有机污染物的装置及探测方法
CN111796331B (zh) * 2020-08-24 2021-07-02 吉林大学 面向浅层地下水和烃类物质的地面磁共振探测装置及方法
CN111796331A (zh) * 2020-08-24 2020-10-20 吉林大学 面向浅层地下水和烃类物质的地面磁共振探测装置及方法
CN113376704A (zh) * 2021-06-07 2021-09-10 电子科技大学 一种基于电发射-磁接收的井间电磁探测***及方法
CN113866837A (zh) * 2021-09-18 2021-12-31 吉林大学 电性源核磁共振与激发极化联用装置以及探测方法
CN113866837B (zh) * 2021-09-18 2023-12-19 吉林大学 电性源核磁共振与激发极化联用装置以及探测方法
CN115097533B (zh) * 2022-05-05 2023-06-30 吉林大学 一种基于tls-esprit算法的磁共振测深信号提取方法
CN115097533A (zh) * 2022-05-05 2022-09-23 吉林大学 一种基于tls-esprit算法的磁共振测深信号提取方法
CN115808648A (zh) * 2022-11-18 2023-03-17 无锡鸣石峻致医疗科技有限公司 一种磁共振***振铃噪声的测量装置和方法
CN116449438A (zh) * 2023-06-16 2023-07-18 中国地质大学（武汉） 一种瞬变电磁法参考消噪装置及方法
CN116449438B (zh) * 2023-06-16 2023-09-12 中国地质大学（武汉） 一种瞬变电磁法参考消噪装置及方法
CN116856920A (zh) * 2023-07-06 2023-10-10 中国科学院地质与地球物理研究所 一种随钻方位电磁波电阻率仪器使用方法及仪器
CN116856920B (zh) * 2023-07-06 2024-04-02 中国科学院地质与地球物理研究所 一种随钻方位电磁波电阻率仪器使用方法及仪器
CN117347416A (zh) * 2023-10-12 2024-01-05 无锡鸣石峻致医疗科技有限公司 一种物质成分的非侵入式检测方法、装置及存储介质
CN117347416B (zh) * 2023-10-12 2024-05-31 无锡鸣石峻致医疗科技有限公司 一种物质成分的非侵入式检测方法、装置及存储介质
CN118067770A (zh) * 2024-04-18 2024-05-24 吉林大学 一种基于低场磁共振技术的粮仓水分监测装置与方法

Also Published As

Publication number Publication date
CN102053280B (zh) 2013-05-01

Similar Documents

Publication Publication Date Title
CN102053280B (zh) 带有参考线圈的核磁共振地下水探测***及探测方法
CN102012525B (zh) 分布式多参数深部电磁断面成像***及测量方法
CN102096111B (zh) 收发天线分离式核磁共振找水装置及找水方法
CN103823244B (zh) 磁共振三分量消噪装置及消噪方法
CN103018781B (zh) 2d/3d核磁共振与瞬变电磁联用仪及野外工作方法
CN102565862B (zh) 一种瞬变电磁响应信号梯度测量方法及观测装置
CN103033849B (zh) 多通道核磁共振地下水探测仪及其野外工作方法
CN104297807B (zh) 地下灾害水源探测磁共振成像装置及探测和成像方法
CN103852794B (zh) 烃类污染浅层地下水磁共振检测装置及检测方法
CN203502602U (zh) 消除工频谐波干扰的核磁共振探测装置
CN103809206A (zh) 核磁共振与瞬变电磁联用地下水探测装置及探测方法
CN102096112A (zh) 基于阵列线圈的核磁共振地下水探测仪及野外探测方法
CN104280780B (zh) 核磁共振与瞬变电磁联用仪及工作方法
CN103955004A (zh) 四通道核磁共振信号全波采集***及采集方法
CN104375191A (zh) 一种大地电磁数据采集***及其采集方法
CN102062876A (zh) 一种全区电偶源频率域电测深方法
CN103344996A (zh) 串联谐振式核磁共振探测装置及探测方法
CN107942397A (zh) 用预极化场增强信号幅度的磁共振多通道探测方法及装置
CN104502984B (zh) 特定频率噪声对消地下核磁共振探测装置及探测方法
CN102520451B (zh) 高精度双通道大地音频电磁物探仪
CN103383466A (zh) 一种煤层水射电测量仪
CN114527512A (zh) 一种频率域无人机多频电磁探测水平梯度采集***
CN104749643B (zh) 一种用于电磁干扰地区的csamt测量方法和装置
CN111856601B (zh) 一种分布式磁共振地下水探测装置及探测方法
CN203759264U (zh) 基于工频整周期触发的核磁共振双极性叠加消噪装置

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20130501

Termination date: 20151110