基于ICA的核磁共振地下水探测信号噪声消除方法技术

本发明专利技术公开了一种核磁共振地下水探测信号噪声消除方法，包括以下步骤：录入三组核磁共振响应数据，分别对这三组数据进行傅里叶变换，确定每组数据核磁共振中心频率附近所含工频谐波，然后构造与工频谐波同频率，与核磁共振相应数据同长度的正弦函数、余弦函数，并与核磁共振响应数据组成观测信号，采用独立分量分析算法对每组观测信号进行分离得到解混信号，进行数据重构以消除工频谐波的干扰，将三组去除工频谐波的核磁共振数据作为观测信号，再利用ICA算法处理，削弱剩余随机噪声干扰。本发明专利技术在消除工频谐波噪声干扰的同时，不会破坏信号的任何细节，不需要铺设参考线圈，操作简单，压制随机噪声不需要大量的数据，减少处理时间。

【技术实现步骤摘要】
基于ICA的核磁共振地下水探测信号噪声消除方法
：本专利技术涉及一种核磁共振数据预处理方法，具体涉及一种基于独立分量分析(ICA)的核磁共振地下水探测信号噪声消除方法。
技术介绍
：寻找和评价地下淡水资源的方法有很多种，常规的物探找水方法都是通过勘查含水构造的层位间接找水，不能解决何处有水、有多少水等一些与地下水紧密相关的基本问题，核磁共振地下水探测技术作为唯一的直接探测地下水的地球物理方法，具有信息量丰富、解唯一性等优点。但是核磁共振信号极其微弱，而且核磁共振探测仪器在接收核磁共振响应信号的同时不可避免地接收到复杂的噪声干扰，诸如雷电等自然干扰源引起的奇异噪声，由电力线、发电机和变电器等人为干扰源引起的工频谐波噪声，以及由其他噪声源引起的平稳噪声和时变噪声。复杂的噪声干扰使得核磁共振地下水探测信号基本被淹没，严重阻碍磁共振信号的提取，影响后期反演解释结果的准确性，因此核磁共振地下水探测信号消噪方法的研究意义重大。目前已经应用到核磁共振地下水探测信号噪声消除的方法有很多种，其中利用非线性能量算子能够对尖峰噪声进行较好地消除，基于大量数据的统计叠加方法能够对随机噪声进行一定的压制，但是需要的数据量大，叠加次数多，工作时间长，效率较低，工频谐波噪声因为与核磁共振地下水探测信号的中心频率非常接近，难以消除。现有的工频谐波噪声消除的方法有陷波器、基于参考线圈和变步长自适应算法、多通道核磁响应信号的自适应抵消，但是陷波器在消除工频谐波噪声的同时也会破坏核磁共振信号中一些重要的细节信息，基于参考线圈和变步长自适应算法，以及多通道核磁响应信号的自适应抵消方法不仅需要铺设参考线圈，操作复杂，而且只能消除与参考通道相关的噪声，对参考通道的依赖性较大。CN203759264U公开“一种基于工频整周期触发的核磁共振双极性叠加消噪装置”，是由双极性叠加消噪装置是由输出可调的大功率电源与大功率发射桥路连接，双极性主控制单元经发射桥路驱动和大功率发射桥路与谐振电容连接，双极性主控制单元分别与工频整周期同步触发单元、高压切换开关、信号调理电路和A/D采集单元连接，高压切换开关经信号调理电路和放大器电路与A/D采集单元连接构成。是通过分别采集两组方向相反的核磁共振信号和方向一致的噪声信号，通过相减的方法消除噪声，但是此方法不能保证两次采集的噪声完全一样，会有噪声残余，影响探测精度。
技术实现思路
：本专利技术的目的就在于针对上述现有技术的不足，结合ICA和核磁共振信号与噪声特性，提供一种基于ICA的核磁共振地下水探测信号噪声消除方法。本专利技术的目的是通过以下技术方案实现的：一种基于ICA的核磁共振地下水探测信号噪声消除方法，包括以下步骤：A、录入三组核磁共振响应数据；B、利用ICA独立成分分析，依次去除每组数据核磁共振信号中心频率附近的工频谐波干扰，然后分别对每组处理后的数据进ICA逆变换后再进行数据重构，得到消除工频谐波干扰的数据；C、将去除工频谐波噪声的三组数据作为观测信号，然后对其进行ICA处理，以消除剩余的随机噪声，并利用ICA逆变换进行数据重构，得到最终的消噪数据。步骤A所述的三组，是指三个源信号。因为，ICA算法的使用前提是观测信号的个数大于等于源信号的个数，为了计算简单通常取相同个数，三组数据相对合理。核磁共振数据包含MRS信号和噪声，其中噪声又包括工频谐波干扰、随机噪声和尖峰干扰等。本专利技术仅针对工频谐波干扰和随机噪声，因此，源信号分为MRS信号、工频谐波干扰和随机噪声。步骤B包括以下步骤：a、对第一组核磁共振响应数据进行快速傅里叶变换，确定该核磁共振地下水探测信号中心频率附近的工频谐波干扰频率f1,f2,…,fn；b、结合所确定频率，分别构造与工频谐波同频率，与核磁共振响应数据同长度的正弦函数、余弦函数，例如工频谐波频率为f1，则正弦函数为：sin(2πf1k)，余弦函数为：cos(2πf1k)，并与核磁共振响应数据构成观测信号；c、对观测信号进行预处理；d、利用快速独立分量分析算法对预处理后的观测信号进行分离，得到幅度有很大衰减的解混信号和解混矩阵w；e、将步骤d分离出的工频谐波噪声屏蔽清零，利用ICA逆变换进行数据重构，恢复核磁共振信号的幅度，其中独立成分yj是分解出来的信号，只保留独立成分yj，其他独立成分置零，是每个通道内所含的信号成分，由于核磁共振响应数据通道所含的核磁共振信号明显，正弦函数及余弦函数所在的通道基本不含核磁共振信号，因此选择核磁共振信号强的作为去工频谐波后的核磁共振信号，记为f、针对第二组核磁共振响应数据，重复步骤a～e得到去除工频谐波的第二组核磁共振响应数据，记为g、针对第三组核磁共振响应数据，重复步骤a～e得到去除工频谐波的第三组核磁共振响应数据，记为步骤C包括以下步骤：Ⅰ、将三组去除工频谐波噪声的数据组成观测信号，并对其进行预处理；Ⅱ、利用快速独立分量分析算法对预处理后的观测信号进行分离，得到幅度有很大衰减的解混信号和解混矩阵A；Ⅲ、将解混信号中噪声成分少、核磁共振信号明显的成分保留，其他成分屏蔽，利用ICA逆变换进行数据重构，恢复核磁共振信号幅度。ICA算法包括以下步骤：第一步、观测信号去均值，首先按照公式求解观测信号每一行的均值，然后利用公式对观测信号进行去均值处理，使数据中心化，满足零均值的假设；第二步、对去均值后的数据进行白化处理，根据公式Cx＝E{xxT}求观测信号的协方差矩阵，并求协方差矩阵的特征值d＝diag(d0,d1,…,d2n)，特征向量e＝(e0,e1,…e2n)，最后由公式x＝vx＝d1/2eTx对零均值观测信号进行白化处理，使数据具有单位方差，以消除数据各分量之间的相关性；第三步、求解混矩阵w，首先令j＝1，初始化解混向量w1，利用根据负熵最大化独立判据和牛顿迭代优化算法推导出递推式：求w1，利用公式对w1进行正交化，在根据wj＝wj/||wj||标准化w1，当w1收敛时，第一个独立成分对应的解混向量w1求解完毕，判断w1是否收敛，如果w1不收敛重新利用公式求w1，对w1进行正交化和标准化，直到w1收敛，第一个独立成分对应的解混向量w1求解完毕，j＝j+1，判断j≤m是否成立，如果j≤m，按照上述的步骤求出w2，直到j＞m时，独立成分对应的所有的解混向量w求解完毕，输出解混矩阵w；第四步、利用公式y＝wTx求解独立成分y；第五步、数据重构，由于利用ICA算法求出的输出信号y具有幅度的不确定性，与源信号相比初始幅度明显减小，将第四步分离出的噪声屏蔽清零，利用ICA逆变换进行数据重构，恢复核磁共振信号的幅度，其中独立成分yj是分解出来的信号，只保留独立成分yj，其他独立成分置零，是每个通道内所含的信号成分。有益效果：本专利技术对源信号和传输通道的先验知识没有要求，在试验过程中不需要铺设参考线圈，使工频噪声的消除效果不受参考线圈噪声采集情况的影响，在随机噪声消除的过程中不需要大量的数据进行叠加，操作简单、方便、效率高，由于采用了ICA算法，在分离噪声的过程中对信号的细节几乎没有破坏。附图说明：图1为基于ICA的核磁共振地下水探测信号工频谐波噪声消除方法流程图图2为用ICA算法对去除工频噪声后随机噪声消除流程图图3为ICA算法流程图图4为对第一组数据利用ICA对工频谐波噪本文档来自技高网...

【技术保护点】
一种基于ICA的核磁共振地下水探测信号噪声消除方法，其特征在于：包括以下步骤：A、录入三组核磁共振响应数据；B、利用ICA独立成分分析，依次去除每组数据核磁共振信号中心频率附近的工频谐波干扰，然后分别对每组处理后的数据进ICA逆变换后再进行数据重构，得到消除工频谐波干扰的数据；C、将去除工频谐波噪声的三组数据作为观测信号，然后对其进行ICA处理，以消除剩余的随机噪声，并利用ICA逆变换进行数据重构，得到最终的消噪数据。

【技术特征摘要】
1.一种基于ICA的核磁共振地下水探测信号噪声消除方法，其特征在于：包括以下步骤：A、录入三组核磁共振响应数据；B、利用ICA独立成分分析，依次去除每组数据核磁共振信号中心频率附近的工频谐波干扰，然后分别对每组处理后的数据进ICA逆变换后再进行数据重构，得到消除工频谐波干扰的数据；C、将去除工频谐波噪声的三组数据作为观测信号，然后对其进行ICA处理，以消除剩余的随机噪声，并利用ICA逆变换进行数据重构，得到最终的消噪数据。2.按照权利要求1所述的基于ICA的核磁共振地下水探测信号噪声消除方法，其特征在于，步骤B包括以下步骤：a、对第一组核磁共振响应数据进行快速傅里叶变换，确定该核磁共振地下水探测信号中心频率附近的工频谐波干扰频率f1，f2，…，fn；b、结合所确定频率，分别构造与工频谐波同频率，与核磁共振响应数据同长度的正弦函数、余弦函数，工频谐波频率为fn，n＝1，2，…，n，则正弦函数为：sin(2πfnk)，余弦函数为：cos(2πfnk)，并与核磁共振响应数据构成观测信号；c、对观测信号进行预处理；d、利用快速独立分量分析算法对预处理后的观测信号进行分...

【专利技术属性】
技术研发人员：林婷婷，张思远，张扬，万玲，蒋川东，林君，
申请(专利权)人：吉林大学，
类型：发明
国别省市：吉林;22