CN112882111B - 一种基于循环相关的磁共振响应信号参数提取方法及系统 - Google Patents
一种基于循环相关的磁共振响应信号参数提取方法及系统 Download PDFInfo
- Publication number
- CN112882111B CN112882111B CN202110059395.2A CN202110059395A CN112882111B CN 112882111 B CN112882111 B CN 112882111B CN 202110059395 A CN202110059395 A CN 202110059395A CN 112882111 B CN112882111 B CN 112882111B
- Authority
- CN
- China
- Prior art keywords
- magnetic resonance
- resonance response
- cyclic correlation
- correlation function
- response signal
- 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
Links
- 125000004122 cyclic group Chemical group 0.000 title claims abstract description 99
- 230000004044 response Effects 0.000 title claims abstract description 94
- 238000000605 extraction Methods 0.000 title claims abstract description 18
- 238000005314 correlation function Methods 0.000 claims abstract description 64
- 238000000034 method Methods 0.000 claims description 60
- 239000011159 matrix material Substances 0.000 claims description 39
- 238000000354 decomposition reaction Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 8
- 238000010276 construction Methods 0.000 claims description 6
- 102000002274 Matrix Metalloproteinases Human genes 0.000 claims description 5
- 108010000684 Matrix Metalloproteinases Proteins 0.000 claims description 5
- 230000008569 process Effects 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 230000005358 geomagnetic field Effects 0.000 claims description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 abstract description 32
- 238000005481 NMR spectroscopy Methods 0.000 abstract description 17
- 238000005516 engineering process Methods 0.000 abstract description 7
- 238000003672 processing method Methods 0.000 abstract 1
- 238000001514 detection method Methods 0.000 description 8
- 238000001914 filtration Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 2
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000005553 drilling Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000005672 electromagnetic field Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 125000004435 hydrogen atom Chemical group [H]* 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000007781 pre-processing Methods 0.000 description 1
- 230000000717 retained effect Effects 0.000 description 1
- 230000007704 transition Effects 0.000 description 1
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
-
- 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
-
- 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
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (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)
- High Energy & Nuclear Physics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明涉及一种基于循环相关的磁共振响应信号的参数提取方法及系统,用于核磁共振找水仪的磁共振响应信号中参数的提取。通过循环相关处理方法,抑制了磁共振响应信号中的高斯白噪声成分和工频谐波噪声成分,从而提高了磁共振响应信号参数的提取精度。然后对得到的待解算循环相关函数进行参数估计,提高了对地下水的勘探精度,为核磁共振测深技术在地下水勘探的应用提供了技术支持。
Description
技术领域
本发明涉及核磁共振技术领域,特别是涉及一种基于循环相关的磁共振响应信号参数提取方法及系统。
背景技术
核磁共振测深技术是目前世界上唯一直接有效探测地下水的地球物理勘探方法。核磁共振探测技术的原理是地下水中的氢核通过人工感应电磁场形成核磁矩。与过去的地下水探测办法相比,利用核磁共振技术找水有如下几点优势:一是与间接探测方法相比,只要有水且水层深度在探测的范围内,基于核磁共振技术的探测方法都能直接探测出结果,找水的效率更高,速度明显更快;二是经过设计,利用核磁共振找水仪能获得更多的地下水的有关信息,比如含水层深度,含水量的多少,地下含水层的孔隙度等等地质信息参数;三是与传统探测技术相比,经济性更强,勘测找水的全过程只需要很短的时间,如果选择钻孔的勘测方案,不仅仅浪费十倍以上的时间,也需要花费近十倍的人力物力。
核磁共振找水仪的基本原理是通过发射线圈向地下水激发能量,使地下水中的氢质子核外电子发生能级跃迁,然后再用接收线圈接收核外电子由高能级向低能级回迁时释放的能量,从而获取核磁共振响应信号,在核磁共振响应信号中提取关键参数,应用反演软件即可反演得到地下水的相关信息。
然而,核磁共振找水仪采集得到的观测数据中除了有磁共振响应信号之外,还有种类多样的人为噪声和环境噪声,产生原因多且十分复杂。噪声干扰无处不在,和我们需要的磁共振响应信号共存,当噪声比较大的时候,甚至信号会淹没在噪声中,给信号的提取和信号特征的分析带来很大的困难。
基于此,本发明提供一种能够有效滤除各种噪声的高精度磁共振响应信号参数提取方法及系统。
发明内容
本发明的目的是提供一种能够有效滤除各种噪声的磁共振响应信号参数提取方法及系统,基于循环相关方法,能够极大地提升磁共振响应数据的信噪比,继而有效提升磁共振响应信号的参数提取精度。
为实现上述目的,本发明提供了如下方案:
一种基于循环相关的磁共振响应信号参数提取方法,所述方法包括:
接收磁共振响应数据;
利用循环相关方法对所述磁共振响应数据进行处理,得到待解算循环相关函数;
利用参数估计方法对所述待解算循环相关函数进行求解,得到磁共振响应信号的平均衰减时间,实现磁共振响应信号参数的提取。
可选的,所述利用循环相关方法对磁共振响应数据进行处理,得到待解算循环相关函数具体包括:
对所述磁共振响应信号进行循环相关处理,得到所述磁共振响应信号的第一循环相关函数;
选取预设循环频率带入所述第一循环相关函数中,得到待解算循环相关函数。
其中,x(n)为磁共振响应数据,n为采样时刻,N为采样点数,τ为移位个数,α为循环频率,n=1,2,…,N。
可选的,所述预设循环频率为仅与磁共振响应信号相关的循环频率。
可选的,所述利用参数估计方法对所述待解算循环相关函数进行求解时,具体采用旋转不变参数估计方法对所述待解算循环相关函数进行求解。
可选的,所述采用旋转不变参数估计方法对所述待解算循环相关函数进行求解具体包括:
根据所述子空间旋转算符Φ计算得到平均衰减时间T2。
可选的,所述根据所述子空间旋转算符Φ计算得到平均衰减时间T2具体包括:
本发明还提供了一种基于循环相关的磁共振响应信号参数提取系统,所述系统包括:
接收模块,用于接收磁共振响应信号;
循环相关函数计算模块,用于利用循环相关方法对所述磁共振响应信号进行处理,得到待解算循环相关函数;
参数估计模块,用于利用参数估计方法对所述待解算循环相关函数进行求解,得到所述磁共振响应信号的平均衰减时间,完成磁共振响应信号参数的提取。
可选的,所述参数估计模块包括旋转不变参数估计单元,所述旋转不变参数估计单元用于具体采用旋转不变参数估计方法对所述待解算循环相关函数进行求解。
可选的,所述旋转不变参数估计单元具体包括:
参数求解子单元,用于根据所述子空间旋转算符Φ计算得到平均衰减时间T2。
根据本发明提供的具体实施例,本发明公开了以下技术效果:
本发明提供了一种磁共振响应信号参数提取方法及系统,通过采用循环相关方法,设置合适的循环频率,有效滤除了不同频率干扰,最大程度的保留了信号信息,极大地提升了磁共振响应信号的信噪比,从而提高了磁共振响应信号参数提取的精确度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种基于循环相关的磁共振响应信号参数提取方法流程图;
图2为本发明实施例提供的一种基于循环相关的磁共振响应信号参数提取系统框图;
图3为本发明实施例提供的旋转不变参数估计单元结构示意图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的目的是提供一种磁共振响应信号参数提取方法及系统,能够有效滤除各种噪声,提高对磁共振响应信号参数提取的精度。
采用本发明的磁共振响应信号参数提取方法及系统能够解决现有技术中利用磁共振响应信号探寻地下水的关键难点问题,即磁共振响应信号极其微弱(通常只有几十纳伏),而工频谐波噪声以及白噪声强度较高,严重影响磁共振响应信号重要参数提取精度的问题。
为使本发明的上述目的、特征和优点能够更加明显易懂,下面结合附图和具体实施方式对本发明作进一步详细的说明。
实施例1
如图1所示,本实施例提供一种基于循环相关的磁共振响应信号参数提取方法,具体包括:
步骤101:接收磁共振响应数据;
步骤102:利用循环相关方法对所述磁共振响应数据进行处理,得到待解算循环相关函数;
步骤103:利用参数估计方法对所述待解算循环相关函数进行求解,得到磁共振响应信号的平均衰减时间,完成磁共振响应信号参数的提取。
本实施例中通过循环相关方法抑制了磁共振响应数据中的高斯白噪声和工频谐波噪声成分,最大程度的保留了磁共振响应信号信息,从而提高了核磁共振找水仪的测量精度。
为了更清楚的展示本实施例提供的参数提取方法的优越性,下面对所述方法进行进一步的解释说明。
首先,将核磁共振找水仪磁共振响应信号的复数信号模型表示如下:
磁共振响应信号中的关键参数与地下水的含量、深度、含水层孔隙度等密切相关。如,初始振幅的大小和地下含水量成正比,也包含了地下含水层深度、厚度、单位体积含水量等信息;平均衰减时间反应了地下含水层平均孔隙度的信息;初始相位反映了含水层的导电率,通常,导电率为一个常数。因此,通过对磁共振响应信号参数进行提取,并通过软件反演就能获取地下水的相关信息。
然而,在实际应用过程中,核磁共振找水仪获取的实测磁共振响应观测信号数据中,往往会附加有白噪声和工频谐波噪声。因此,接收到的核磁共振找水仪的磁共振响应数据x(n)可以表示为:
x(n)=s(n)+η(n)+ε(n) (2)
其中,η(n)为高斯白噪声,ε(n)为工频谐波噪声。工频谐波噪声主要来源于变压器、电线、工厂及居民用电。离散的工频谐波噪声表达形式为其中fc为工频谐波的基频,通常为50Hz或60Hz,Ab和φb分别为工频谐波分量的幅值和相位,B为工频谐波分量个数。
在接收到磁共振响应数据x(n)后,可以对x(n)进行预处理,滤除x(n)中明显易去除的噪声,提高精确度的同时为后续步骤的运算处理降低复杂度。
然后利用循环相关方法对磁共振响应数据x(n)进行处理,抑制x(n)中的高斯白噪声和工频谐波噪声,具体包括:
其中,α为循环相关函数的循环频率,τ为移位个数。
然后将包含不含噪声的磁共振响应信号s(n)、白噪声η(n)和工频谐波噪声ε(n)三部分的磁共振响应数据x(n)带入到(3)式中:
由(4)式可以看出,具有相关性的子式仅有s(n)s*(n+τ)、η(n)η*(n+τ)和ε(n)ε*(n+τ),根据循环相关函数性质,可知其余项的循环相关函数均为0。由此对第一循环相关函数可以变换为:
其中,s(n)s*(n+τ)子项对应的循环相关函数为:
由公式(6)可以看出,当循环频率α取不同值时,有:
取α=2fL为预设循环频率,将预设循环频率带入所述公式(7)中,就得到了仅与s(n)s*(n+τ)子项相关的循环函数:
其中,fL为地磁场的拉摩尔频率,在采集磁共振响应信号前由地磁仪测量得到。
若工频谐波频率和白噪声频率均与磁共振响应信号频率不同,则取循环频率为α=2fL时,公式(5)中的η(n)η*(n+τ)和ε(n)ε*(n+τ)子项所对应的循环相关函数均为0。因此,当循环频率为α=2fL时,公式(5)中的所述第一循环相关函数为:
公式(9)即为待解算循环相关函数。
本实施例中,具体选择预设循环频率为2fL。但本领域的技术人员应当理解,预设循环频率的选择是与磁共振响应信号s(n)的特性密切相关的,本实施例中将2fL作为预设循环频率仅仅是为了更加清楚的对本发明的技术方案进行说明,并不应理解为对本发明的具体限制,实际上,选择任何能够实现噪声滤除的预设循环频率都将落入本发明的保护范围内。
本实施例中通过结合磁共振响应信号的自身特性,对磁共振响应数据进行划分来选择预设循环频率,将循环相关函数中与循环频率无关的白噪声和工频谐波噪声部分进行剔除,从而仅保留了与循环频率相关的磁共振响应信号信息,为进一步的参数估计提供了准确、简洁的数据基础,从而提高了磁共振响应信号参数提取的精度,进而提高了核磁共振找水仪对地下水的检测精度。
在得到待解算循环相关函数后,利用参数估计方法对待解算循环相关函数进行求解。为了更清楚的对求解过程进行解释,本实施例以具体采用旋转不变参数估计方法对所述待解算循环相关函数进行求解为例进行详细说明。
ESPRIT方法借助旋转不变技术估计信号参数(Estimation of SignalParameters via Rotational Invariance Techniques),采用旋转不变参数估计方法对所述待解算循环相关函数进行求解具体包括:
当然,采用旋转不变信号参数估计方法对待解算循环相关函数进行求解来得到磁共振响应信号的平均衰减时间仅是本实施例所提供的一个具体实施方式,并不是对本发明保护范围的进一步限制。其他参数估计方法如将循环相关函数的结构特征作为先验信息,构建冗余字典来描述循环频率结构来对循环相关函数进行估计等也落入本发明的保护范围内,只要能够实现对待解算循环相关函数进行求解估计即可。
实施例2
本实施例提供了一种基于循环相关的磁共振响应信号参数提取系统,如图2所示,所述系统包括:
接收模块M1,用于接收磁共振响应数据;
循环相关函数计算模块M2,用于利用循环相关方法对所述磁共振响应数据进行处理,得到待解算循环相关函数;
参数估计模块M3,用于利用参数估计方法对所述待解算循环相关函数进行求解,得到所述磁共振响应信号的平均衰减时间,完成磁共振响应信号参数的提取。
为了更好的实现功能,所述参数估计模块M2包括旋转不变参数估计单元,所述旋转不变参数估计单元用于具体采用旋转不变参数估计方法对所述待解算循环相关函数进行求解。
如图3所示,所述旋转不变参数估计单元具体包括:
参数求解子单元M34,用于根据所述子空间旋转算符Φ计算得到平均衰减时间T2。
本说明书中每个实施例重点说明的都是与其他实施例的不同之处,各个实施例之间相同相似部分互相参见即可。对于实施例公开的系统而言,由于其与实施例公开的方法相对应,所以描述的比较简单,相关之处参见方法部分说明即可。
本文中应用了具体个例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处。综上所述,本说明书内容不应理解为对本发明的限制。
Claims (8)
1.一种基于循环相关的磁共振响应信号参数提取方法,其特征在于,所述方法包括:
接收磁共振响应数据;
利用循环相关方法对所述磁共振响应数据进行处理,得到待解算循环相关函数,具体包括:对所述磁共振响应数据进行循环相关处理,得到所述磁共振响应数据的第一循环相关函数,所述磁共振响应数据的第一循环相关函数为:其中,x(n)为磁共振响应数据,n为采样时刻,N为采样点数,τ为移位个数,α为循环频率,n=1,2,…,N;选取预设循环频率带入所述第一循环相关函数中,得到待解算循环相关函数,当循环频率α取不同值时,有:
取α=2fL为预设循环频率,将预设循环频率代入以上公式中,就得到了仅与s(n)s*(n+τ)子项相关的循环函数:
其中fL为地磁场的拉摩尔频率,在采集磁共振响应信号前由地磁仪测量得到;
利用参数估计方法对所述待解算循环相关函数进行求解,得到子空间旋转算符Φ,所述子空间旋转算符Φ计算得到磁共振响应信号的平均衰减时间,完成所述磁共振响应信号参数的提取。
2.根据权利要求1所述的基于循环相关的磁共振响应信号参数提取方法,其特征在于,所述预设循环频率为仅与磁共振响应信号相关的循环频率。
3.根据权利要求1所述的基于循环相关的磁共振响应信号参数提取方法,其特征在于,所述利用参数估计方法对所述待解算循环相关函数进行求解时,具体采用旋转不变参数估计方法对所述待解算循环相关函数进行求解。
6.一种基于循环相关的磁共振响应信号参数提取系统,其特征在于,所述系统包括:
接收模块,用于接收磁共振响应数据;
循环相关函数计算模块,用于利用循环相关方法对所述磁共振响应数据进行处理,得到待解算循环相关函数,具体包括:对所述磁共振响应数据进行循环相关处理,得到所述磁共振响应数据的第一循环相关函数,所述磁共振响应数据的第一循环相关函数为:其中,x(n)为磁共振响应数据,n为采样时刻,N为采样点数,τ为移位个数,α为循环频率,n=1,2,…,N;选取预设循环频率带入所述第一循环相关函数中,得到待解算循环相关函数,当循环频率α取不同值时,有:
取α=2fL为预设循环频率,将预设循环频率代入以上公式中,就得到了仅与s(n)s*(n+τ)子项相关的循环函数:
其中fL为地磁场的拉摩尔频率,在采集磁共振响应信号前由地磁仪测量得到;
参数估计模块,用于利用参数估计方法对所述待解算循环相关函数进行求解,得到子空间旋转算符Φ,所述子空间旋转算符Φ计算得到所述磁共振响应信号的平均衰减时间,完成磁共振响应信号参数的提取。
7.根据权利要求6所述的基于循环相关的磁共振响应信号参数提取系统,其特征在于,所述参数估计模块包括旋转不变参数估计单元,所述旋转不变参数估计单元用于具体采用旋转不变参数估计方法对所述待解算循环相关函数进行求解。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110059395.2A CN112882111B (zh) | 2021-01-18 | 2021-01-18 | 一种基于循环相关的磁共振响应信号参数提取方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110059395.2A CN112882111B (zh) | 2021-01-18 | 2021-01-18 | 一种基于循环相关的磁共振响应信号参数提取方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112882111A CN112882111A (zh) | 2021-06-01 |
CN112882111B true CN112882111B (zh) | 2022-05-03 |
Family
ID=76048638
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110059395.2A Active CN112882111B (zh) | 2021-01-18 | 2021-01-18 | 一种基于循环相关的磁共振响应信号参数提取方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112882111B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115097533B (zh) * | 2022-05-05 | 2023-06-30 | 吉林大学 | 一种基于tls-esprit算法的磁共振测深信号提取方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103944174A (zh) * | 2014-04-28 | 2014-07-23 | 四川大学 | 基于互相关函数滤噪算法的低频振荡在线辨识方法 |
CN105301354A (zh) * | 2015-07-12 | 2016-02-03 | 九江学院 | 一种乘性和加性噪声中谐波信号频率估计方法 |
CN106991647A (zh) * | 2017-03-29 | 2017-07-28 | 华中科技大学 | 一种低信噪比图像重构方法和系统 |
CN109613617A (zh) * | 2019-01-24 | 2019-04-12 | 吉林大学 | 基于磁共振响应信号参数提取的地下水探测方法及系统 |
CN111913034A (zh) * | 2020-06-18 | 2020-11-10 | 江苏方天电力技术有限公司 | 一种基于高阶累积量与esprit算法的功率振荡检测方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3966860B2 (ja) * | 2004-02-12 | 2007-08-29 | 株式会社国際電気通信基礎技術研究所 | 電波到来方向探知装置及び方法 |
JP2013152239A (ja) * | 2013-03-29 | 2013-08-08 | Honda Elesys Co Ltd | 電子走査型レーダ装置、その制御方法及びプログラム |
CN105044453B (zh) * | 2015-08-11 | 2018-05-01 | 杨世永 | 一种适用于复杂噪声背景的谐波信号频率估计方法 |
CN107607998B (zh) * | 2017-09-25 | 2019-02-26 | 吉林大学 | 一种核磁共振找水仪磁共振响应信号参数提取方法及系统 |
CN109239648A (zh) * | 2018-10-24 | 2019-01-18 | 同方电子科技有限公司 | 基于对称循环频率的谱相关子空间测向方法 |
CN111308438B (zh) * | 2020-03-06 | 2021-11-02 | 中国人民解放军海军航空大学 | 一种回波数据散射特征提取方法及系统 |
-
2021
- 2021-01-18 CN CN202110059395.2A patent/CN112882111B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103944174A (zh) * | 2014-04-28 | 2014-07-23 | 四川大学 | 基于互相关函数滤噪算法的低频振荡在线辨识方法 |
CN105301354A (zh) * | 2015-07-12 | 2016-02-03 | 九江学院 | 一种乘性和加性噪声中谐波信号频率估计方法 |
CN106991647A (zh) * | 2017-03-29 | 2017-07-28 | 华中科技大学 | 一种低信噪比图像重构方法和系统 |
CN109613617A (zh) * | 2019-01-24 | 2019-04-12 | 吉林大学 | 基于磁共振响应信号参数提取的地下水探测方法及系统 |
CN111913034A (zh) * | 2020-06-18 | 2020-11-10 | 江苏方天电力技术有限公司 | 一种基于高阶累积量与esprit算法的功率振荡检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112882111A (zh) | 2021-06-01 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN106772646B (zh) | 一种地面核磁共振信号提取方法 | |
US8581587B2 (en) | SNMR pulse sequence phase cycling | |
CN104330433B (zh) | 一种获得目的储层横向弛豫时间分布的方法及装置 | |
CN103823244A (zh) | 磁共振三分量消噪装置及消噪方法 | |
CN104537714A (zh) | 磁共振与瞬变电磁空间约束联合反演方法 | |
CN112882111B (zh) | 一种基于循环相关的磁共振响应信号参数提取方法及系统 | |
CN103513285A (zh) | 一种确定横向表面弛豫速率的方法及装置 | |
CN110955982B (zh) | 一种变质岩储层渗透率计算方法、装置及计算机存储介质 | |
CN107607998B (zh) | 一种核磁共振找水仪磁共振响应信号参数提取方法及系统 | |
Larsen et al. | Processing of surface-nuclear magnetic resonance data from sites with high noise levels | |
CN109613617B (zh) | 基于磁共振响应信号参数提取的地下水探测方法及系统 | |
CN110454153A (zh) | 一种核磁共振测井弛豫反演方法 | |
Grombacher et al. | Steady‐State Surface NMR for Mapping of Groundwater | |
CN106443785A (zh) | 绕射波场获取方法和装置 | |
CN111335882B (zh) | 一种融合自然γ探测器的铀裂变瞬发中子测井控制方法 | |
CN112649887B (zh) | 基于钻井资料定量划分层序地层的方法及装置 | |
CN116956119A (zh) | 一种基于沉积噪声模拟的泥页岩层序地层划分方法及装置 | |
CN103995293B (zh) | 磁共振测深信号检测方法 | |
CN115826068A (zh) | 一种基于自适应高斯滤波的mrs信号包络提取方法 | |
CN112946760B (zh) | 基于正则化方法的未爆弹三维立体成像方法、装置及系统 | |
CN108897046B (zh) | 远探测声波逆时偏移成像方法及装置 | |
CN112710688B (zh) | 核磁共振纵向弛豫获取方法及系统 | |
CN110609332A (zh) | 地层数据采集方法、装置以及系统 | |
Li et al. | Research on the detection method of MRS Signal initial amplitude based on chaotic detection system | |
Lai et al. | Monte Carlo CLEAN spectral analysis method for detecting the stratigraphic cycles based on high-resolution electrical image log data |
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 |