吉林大学学报(地球科学版)

• 地球探测与信息技术 • 上一篇    下一篇

基于不等式约束的磁共振信号T2谱多指数分解法及算法改进

林君,慧芳,孙淑琴,蒋川东,林婷婷   

  1. 地球信息探测仪器教育部重点实验室/吉林大学仪器科学与电气工程学院,长春130026
  • 收稿日期:2012-12-29 出版日期:2013-11-26 发布日期:2013-11-26
  • 通讯作者: 孙淑琴(1969-),女,副教授,硕士生导师,主要从事数字信号处理研究,Tel:0431-88502987 E-mail:sunsq@jlu.edu.cn
  • 作者简介:林君(1954-),男,教授,博士生导师,主要从事地球物理探测技术及仪器研究,Tel:0431-88502473,E-mail:lin_jun@jlu.edu.cn
  • 基金资助:

    科技部科技创新工作方法项目(3B810Z950537);国家重大科学仪器设备开发和应用专项(2011YQ030133)

MRS Signal T2 Inversion by Multi-Exponential Decomposition with Inequality Constraint

Lin Jun, Hui Fang, Sun Shuqin, Jiang Chuandong, Lin Tingting   

  1. Lab. of Geo-Exploration and Instrumentation Ministry of Education/College of Instrumentation and Electrical Engineering, Jilin University, Changchun130026,China
  • Received:2012-12-29 Online:2013-11-26 Published:2013-11-26

摘要:

磁共振地下水探测信号的弛豫时间T2和振幅与地下岩石孔隙大小和地下水的含水量直接相关。目前,磁共振地下水探测方法中,信号的衰减特性由平均弛豫时间T*2来描述。然而,对于复杂的多孔介质组成的岩石结构,仅由T*2无法准确描述地下孔隙情况。详细研究了基于不等式约束的多指数分解法,进行了信号T2谱的提取,并重新计算多指数形式下的信号初始振幅。结果表明:对比信号在不同采集时间下的T2谱提取结果,延长信号的采集时间有助于T2谱的准确提取;通过磁共振仿真信号计算并分析不同信噪比下的T2谱提取结果,信噪比大于30 dB时相对误差平均值小于5%;进行算法改进,得到信噪比大于10 dB时相对误差平均值小于5%的结果,且幅值和T2的相对误差平均值降低了约1/3。由此可知:采用多指数分解法能够获得磁共振信号的T2谱,正则化技术和加权修正能够提高T2谱的准确度和抗干扰能力;与单指数方式相比,采用多指数拟合能够更准确地恢复信号初始振幅。

关键词: 磁共振测深, T2谱, 多指数分解, 不等式约束, 加权修正

Abstract:

The amplitude and relaxation time T2 of magnetic resonance sounding (MRS) signal are directly related with water content and pore size of aquifers. Average T2, denoted by T*2, is often used to describe the decay characteristic of MRS signal. However, in the case of porous media rock, pore size should be described by T2 spectrum of MRS signal. To obtain the spectrum and estimate the initial amplitude more accurately, we studied multi-exponential decomposition method based on inequality constraint in this paper. Numerical results show that extending signal length helps the accuracy of T2 spectrum estimation. With traditional method, average relative error is less than 5% when SNR is greater than 30 dB. With our modified algorithm, a SNR about 10 dB could lead to the same error. meanwhile, the error of the estimated amplitude is reduced by at least 1/3. It can be concluded that multi-exponential decomposition method could obtain T2 spectrum of MRS signal, and adding regularization and weighted correction to the method could improve the accuracy of T2 spectrum inversion. Comparing to mono-exponential fitting, multi-exponential inversion can estimate the initial amplitude more accurately.

Key words: magnetic resonance sounding, T2 spectrum, multi-exponential decomposition, inequality constraint, weighted correction

中图分类号: 

  • TN911.7
[1] 林君, 张扬, 张思远, 舒旭, 杜文元, 林婷婷. 强电磁干扰下磁共振地下水探测噪声压制方法研究进展[J]. 吉林大学学报(地球科学版), 2016, 46(4): 1221-1230.
[2] 林君, 蒋川东, 段清明, 王应吉, 秦胜伍, 林婷婷. 复杂条件下地下水磁共振探测与灾害水源探查研究进展[J]. J4, 2012, 42(5): 1560-1570.
[3] 李鹏举, 姜大鹏, 李鹏飞, 宋延杰. 核磁共振T2谱M-P广义逆反演研究[J]. J4, 2011, 41(2): 579-585.
[4] 林峰, 王祝文, 刘菁华, 丁阳, 李长春. 核磁共振T2谱奇异值反演改进算法[J]. J4, 2009, 39(6): 1150-1155.
Viewed
Full text


Abstract

Cited

  Shared   
  Discussed   
No Suggested Reading articles found!