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

说明书

本发明提供一种快速二维核磁共振测井信号处理及T1T2谱反演方法，包括：S1，从每个组分的回波信号中按照设定规则抽取部分回波信号，得到抽取的回波信号；S2，根据核磁共振采集模式并利用抽取的回波信号构建系数矩阵A；S3，对系数矩阵A进行T1T2谱求解；S4，输出T1T2谱。本发明方法针对核磁共振回波信号特点，通过信号抽取，实现了大矩阵的压缩，缩短了反演时间。

技术领域

本发明属于地球物理测井领域，涉及一种快速二维核磁共振测井信号处理及T1T2谱反演方法。

背景技术

一维核磁共振测井在油气识别和孔隙结构评价中已经得到广泛应用，但在实际应用中也存在一系列不足，比如多相流体存在时流体识别率偏低等问题凸显。由此发展了二维核磁共振测井方法，通过二维核磁信号处理得到二维谱来实现对于多相流体的识别。

采用不同等待时间进行回波信号测量，对回波信号构建方程式，通过反演可以得到T1T2谱，对应理论图谱或岩心实验图谱可以完成流体的准确识别。

目前二维核磁测井资料的处理技术并不成熟，较一维核磁回波信号相比，二维核磁回波信号数据量成倍增长，在处理效率上很低，无法投入正常生产应用。

发明内容

本发明的主要目的为实现一种快速二维核磁共振测井信号处理及T1T2谱反演方法，解决二维核磁反演时间长、效率低问题。

本发明是通过以下技术方案来实现：

一种快速二维核磁共振测井信号处理及T1T2谱反演方法，包括：

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

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

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

S4，输出T1T2谱。

优选的，S1之前，对回波信号进行滤波处理，S1中从每个组分的滤波处理后的回波信号中按照设定规则抽取部分回波信号，得到抽取的回波信号。

进一步的，所述滤波为小波阀值滤波。

进一步的，小波阀值滤波采用小波软阀值滤波，并采用heursure阀值选择规则。

进一步的，滤波处理具体为：

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

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

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

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

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

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

进一步的，N为50，N’为150。

优选的，S2中，系数矩阵A是n×m矩阵。

进一步的，S3具体为：

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

AX＝Y

A＝UVS^T