[发明专利]一种快速二维核磁共振测井信号处理及T1T2谱反演方法

权利要求书

1.一种快速二维核磁共振测井信号处理及T1T2谱反演方法，其特征在于，包括：

S1，从每个组分的回波信号中按照设定规则抽取部分回波信号，得到抽取的回波信号；

S2，根据核磁共振采集模式并利用抽取的回波信号构建系数矩阵A；

S3，对系数矩阵A进行T1T2谱求解；

S4，输出T1T2谱；

S1中，对于每个组分，设定规则为：

(1)若回波个数ECHONUMN，则保留全部回波信号；

(2)若回波个数ECHONUM＝N，保留前N个点的回波信号，对于NECHONUM＝N’部分的回波信号，保留回波数为奇数的回波信号，对于ECHONUMN’部分的回波信号，保留ECHONUM/4为整数的回波信号；

S2中，系数矩阵A是n×m矩阵；

S3具体为：

S31，对系数矩阵A奇异值分解，得到两个正交矩阵U,S和对角阵V；其中U是n×m列正交矩阵，S是m×m阶正交矩阵，是m阶对角矩阵，其中是奇异值；

AX＝Y

A＝UVS^T

式中，X为向量型式的方程解，Y是抽取的回波信号；

S32，计算回波信号信噪比SNR，X在处进行信号截断，保留小于的奇异值，初始解表示为

S33，设定迭代计算最大次数inter，设定迭代步进因子alpha，精度要求eprecison；令迭代计算次数k＝1；

S34，进行第k次迭代计算：对X进行非负判断，保留X(i)＞＝0对应的系数矩阵A中的列，得到新的系数矩阵A；同时令X中小于零的项改为0.001，如果||AX-Y||₂＜epricison，转至S36；如果A为空，转至S36；如果k＝inter，转至S36；

S35，计算矩阵M＝A*(A)^T，计算新的X＝(M+alpha*eye(mRow))\Y，其中eye(mRow)表示阶数为mRow的单位矩阵，令k＝k+1，转S34；

S36，输出X。

2.根据权利要求1所述的快速二维核磁共振测井信号处理及T1T2谱反演方法，其特征在于，S1之前，对回波信号进行滤波处理，S1中从每个组分的滤波处理后的回波信号中按照设定规则抽取部分回波信号，得到抽取的回波信号。

3.根据权利要求2所述的快速二维核磁共振测井信号处理及T1T2谱反演方法，其特征在于，所述滤波为小波阈值滤波。

4.根据权利要求3所述的快速二维核磁共振测井信号处理及T1T2谱反演方法，其特征在于，小波阈值滤波采用小波软阈值滤波，并采用heursure阈值选择规则。

5.根据权利要求4所述的快速二维核磁共振测井信号处理及T1T2谱反演方法，其特征在于，滤波处理具体为：

S01，对回波信号进行小波分解，得到不同分解层数的小波，进行小波滤波计算，得到不同分解层数下的小波滤波结果；

S02，计算不同分解层数下小波滤波结果的小波定量评价指标IDC；

S03，输出最大IDC值对应的小波滤波结果作为滤波处理后的回波信号。

6.根据权利要求1所述的快速二维核磁共振测井信号处理及T1T2谱反演方法，其特征在于，N为50，N’为150。

7.根据权利要求1所述的快速二维核磁共振测井信号处理及T1T2谱反演方法，其特征在于，设定迭代计算最大次数inter＝500，设定迭代步进因子alpha＝0.01，精度要求eprecison＝0.0001。