[发明专利]T2谱反演方法、装置和存储介质

说明书

本发明公开了一种T2谱反演方法、装置和存储介质，该方法包括：接收用户输入的布点数据；根据所述布点数据及相关采集参数，确定反演系数矩阵；采集多组回波数据，并根据所述多组回波数据分别对所述反演系数矩阵进行扩充处理，得到扩充后的反演系数矩阵；根据所述扩充后的反演系数矩阵，使用正负扫描算子，进行非负约束确定横向弛豫分量的初始幅度P。本发明实现了T2谱的反演，并提高了T2谱反演的速度和准确率。

技术领域

本发明涉及地质勘探的测井技术领域，尤其涉及一种T2谱反演方法、装置和存储介质。

背景技术

采用核磁共振(Nuclear magnetic resonance，NMR)方式测井时，实际仪器测得的每个回波信号都是多种弛豫组分的总体效应。这可以用一个多指数函数来表示：式中，A(t)为t时刻测到的回波幅度；T_2i为第i种弛豫分量的横向弛豫时间，其用于表征孔隙尺寸大小；P_i为第i种横向弛豫分量的初始幅度，其用于表征某种特征弛豫尺寸孔隙的孔隙度大小。其中，T_2i在反演中，是先期假定的一系列的值，经过确定各特征弛豫和特征弛豫组分i后，结合回波串，j＝1，…，k(k为回波数)可以构成超定方程组，通过该方程组求出P_i的过程叫做T2谱的反演，它是核磁共振测井资料处理的关键。T2谱的反演是一个典型的病态问题，或称非适定性问题。

近年来，NMR的T2谱反演方法研究取得了很大进展，国内外出现了许多不同的反演算法，目前，基于奇异值分解法的MAP-II仍是目前国际上广泛采用的T2谱反演方法。该方法通过构建T2谱反演线性方程组，利用奇异值的分解法针对实际问题求解该线性方程组的解，最终获得连续性较好的反演T2谱。

但是，在现有技术的具体实施过程中，由于实测回波串比较多，内存开销大，反演速度较慢。为了提高信噪比与反演速度，通常需要对回波串进行滤波，然后抽取一部分数据进行反演，同时压缩T2布点个数，最终对反演结果进行插值与平滑，其中间环节影响大，效率也低。

发明内容

本发明提供一种T2谱反演算法方法、装置和存储介质，以实现简洁高效的T2谱反演。

第一方面，本发明实施例提供一种T2谱反演方法，包括：

接收用户输入的布点数据；

根据所述布点数据，确定反演系数矩阵；

采集多组回波数据，并根据所述多组回波数据分别对所述反演系数矩阵进行扩充处理，得到扩充后的反演系数矩阵；

根据所述扩充后的反演系数矩阵，确定横向弛豫分量的初始幅度P。

在本方案中，通过从接收用户输入的布点数据，并根据所述布点数据，确定反演系数矩阵，避免了回波参数中根据每个处理点来回建立反演矩阵，然后，采集多组回波数据，并根据所述多组回波数据分别对所述反演系数矩阵进行扩充处理，得到扩充后的反演系数矩阵，最后根据所述扩充后的反演系数矩阵，确定横向弛豫分量的初始幅度P，实现了T2谱的反演，并通过建立反演系数矩阵提高了T2谱反演的效率。

可选的，所述根据所述布点数据，确定反演系数矩阵，包括：

接收用户输入的对数最小值T_2min和对数最大值T_2max；

根据所述布点数据、所述对数最小值T_2min和所述对数最大值T_2max，确定第i种弛豫分量的横向弛豫时间T_2i；

根据所述T_2i，确定所述反演系数矩阵。

可选的，所述根据所述布点数据，确定反演系数矩阵，包括：

接收用户输入的指数最小值T_2emin和指数步长增量ΔT_2e；