我国淡水资源极为匮乏,高效准确地探测到地下可利用的淡水资源是解决水资源缺乏的有效途径。磁共振探测技术能直接探测水中氢质子的拉莫尔进动,具有非侵入式、能直接探测地下水的优势,且能反演得到地下水含水量、深度,孔隙度和渗透率等水文地质参数,获得越来越广泛的应用。然而,磁共振探测接收信号强度是纳伏级的,极易受到环境中的电磁噪声干扰。如果噪声干扰不能得到有效地去除,会导致后期数据反演解释结果不准确。影响磁共振地下水探测的噪声主要包括尖峰噪声、工频谐波噪声和随机噪声。国内外已有的相关文献,大多是在信号和噪声特性已知的情况下对噪声进行去除的,但是某些噪声并没有统一的规律,且混合机理具有不确定性,这会导致算法的噪声消除效率低、效果不理想。针对以上问题,由于盲源分离方法能够在源信号与混合过程均未知的情况下,从观测得到的多源混合信号中分离并恢复出源信号的各独立分量,因此采取盲源分离的方式对磁共振中的信号和噪声进行分离。在众多的盲源分离方法中,由于快速独立分量分析(Fast Independent Component Analysis,FastICA)算法具有良好的稳定性,并且收敛速度快,可以批处理数据,所以选择FastICA作为磁共振信号与噪声分离方法。本文首先针对FastICA幅度不确定性问题,提出利用数据重构法还原独立成分的幅度,并将该方法与频谱校正法进行了对比,仿真结果显示,数据重构法还原幅度的误差小于频谱校正法,尤其是在随机噪声较大的情况下,例如信噪比低于-37.1 d B时,频谱校正法的误差较大,甚至大于10%,超出误差允许的范围,整个实验中数据重构法误差均小于5%,满足实际应用需求。其次针对工频谐波噪声进行处理,并完成了仿真实验。基于FastICA的磁共振工频谐波噪声去除分为两种情况,一种是当工频谐波噪声频率与磁共振信号拉莫尔频率有一定的偏差,即两者的相关性不强时,将此时的工频噪声记为非同频噪声;另外一种是当工频谐波噪声频率与磁共振信号拉莫尔频率重合,即两者相关性很强时,将此时的工频谐波噪声记为同频噪声。非同频噪声满足FastICA的假设条件,可以通过高斯插值法对工频谐波的频率进行估计,并构造与工频谐波噪声同频率的正弦信号与余弦信号作为参考通道,以实现磁共振信号与工频谐波噪声的盲源分离,将其结果与陷波器进行对比,FastICA去除工频谐波噪声后数据的信噪比提高了43.42 d B,陷波器消噪后数据的信噪比提高了36.75 d B。结果表明FastICA算法不仅消噪效果好,而且不会对磁共振信号产生破坏,能够很好地保留磁共振信号携带的信息。同频噪声存在的情况也是磁共振工频消除的一个国际难点,此时需要先利用FastICA算法对其他频率的工频谐波噪声进行去除,然后截取数据的后一段,这段数据中磁共振信号已经衰减为零,数据中只含有同频噪声和随机噪声,采用高斯插值法对该工频的频率进行估计,利用相关检测的方法对该工频谐波噪声的幅度和相位进行计算,通过估计出的频率、幅度、相位可构造出该工频谐波噪声波形,且误差较小,最后将估计出的工频谐波噪声从去噪后的数据中减掉,即可实现同频噪声分离。去除工频谐波噪声后的数据中仍然残留着一些随机噪声。利用FastICA对磁共振随机噪声进行压制,仿真实验中仅需要两组观测数据,就能实现对高斯随机噪声的分离,这不仅能够提高数据的信噪比,还能大大地缩短野外探测时间,提高仪器的探测效率。最后分别在北湖公园和烧锅镇这两个地方进行野外探测,并利用FastICA算法对探测数据中的工频谐波噪声和随机噪声进行处理,最终北湖公园实验信噪比提高了27.18 d B,烧锅镇实验信噪比提高了33.53 d B,验证了该算法的有效性和实用性。
【学位单位】:吉林大学
【学位级别】:硕士
【学位年份】:2017
【中图分类】:TH834
【部分图文】:
图 1.1 关中平原 2016 年地下水埋深等值面图使用地球物理方法进行探测可以提高找水定井准确率,探测方法可分为非直接式和直接式找水[5]。非直接式找水方法主要包括电阻率法、激发极化法、时间域电磁法和频率域电磁法等[6,7],他们通过比较判断含水裂隙和断层的深度、宽度

图 1.2 河西走廊平原 2016 年地下水埋深等值面图磁共振测量信号的幅度与静态磁场强度的平方成正比,由于地磁场强度一般位于 25 – 65 μT,磁共振探测接收信号强度为 nV 级别,极其微弱。在进行磁共振野外观测时,不仅由氢质子在激发磁场的作用下产生的弛豫场在接收线圈中产

原子磁矩在磁场中进动示意图
【参考文献】
相关期刊论文 前10条
1 吕睿;;浅谈我国水资源保护[J];黑河学刊;2017年01期
2 万玲;张扬;林君;蒋川东;林婷婷;;基于能量运算的磁共振信号尖峰噪声抑制方法[J];地球物理学报;2016年06期
3 田宝凤;周媛媛;王悦;李振宇;易晓峰;;基于独立成分分析的全波核磁共振信号噪声滤除方法研究[J];物理学报;2015年22期
4 刘佳骏;董锁成;李泽红;;中国水资源承载力综合评价研究[J];自然资源学报;2011年02期
5 苗硕;;中国淡水资源现状与保护措施探讨[J];现代商贸工业;2010年17期
6 林君;;核磁共振找水技术的研究现状与发展趋势[J];地球物理学进展;2010年02期
7 王中兴;荣亮亮;林君;;地面核磁共振找水信号中的奇异干扰抑制[J];吉林大学学报(工学版);2009年05期
8 王平;宋希庚;薛冬新;唐斌;石磊;丁春雨;胡朝阳;;基于独立分量分析和时频分析的柴油机燃烧噪声研究Ⅰ:燃烧噪声分析[J];内燃机学报;2009年04期
9 丁二峰;刘颖;;利用浅层地下水缓解淡水资源短缺的探讨[J];地下水;2009年01期
10 魏巍;刘学伟;;基于独立分量分析的工频干扰消除技术[J];计算机应用研究;2009年01期
相关博士学位论文 前1条
1 武智瑛;多组分气体光声光谱的独立成分分析方法研究[D];大连理工大学;2013年
相关硕士学位论文 前8条
1 周媛媛;基于ICA的全波MRS信号消噪方法研究[D];吉林大学;2016年
2 曲永星;阵列式宽带核磁共振全波采集系统的研制[D];吉林大学;2015年
3 訾彦勇;基于EMD的磁共振探测信号噪声抑制方法研究[D];吉林大学;2015年
4 王微;基于小波分析和统计理论的核磁共振信号消噪方法研究[D];吉林大学;2013年
5 王健鹏;核磁共振地下水探测仪模拟跟踪滤波器研制[D];吉林大学;2012年
6 蒋川东;核磁共振地下水探测系统数据处理软件的设计与应用[D];吉林大学;2009年
7 任鹏飞;基于对消法的核磁共振找水仪工频谐波干扰滤除模块的设计[D];吉林大学;2009年
8 吴逊;基于独立成分分析的特征提取方法研究[D];厦门大学;2007年
本文编号:
2869308
本文链接:https://www.wllwen.com/kejilunwen/yiqiyibiao/2869308.html