[发明专利]基于TwIST的二维核磁共振谱的反演方法和装置

说明书

本发明提供一种基于TwIST的二维核磁共振谱的反演方法和装置，该方法包括：通过N组核磁共振测量参数，生成核矩阵，建立用于反演二维核磁共振谱的目标函数；在每一次对正则化参数进行取值时，采用TwIST算法，求取目标函数的解；根据目标函数的每一个解，确定与每一个解对应的数据比值；根据每一个数据比值、以及正则化参数每一次的取值，建立曲线关系；计算曲线关系中与正则化参数每一次的取值对应的斜率，并根据斜率确定正则化参数的最优取值；确定与最优取值对应的目标函数的解，并根据与最优取值对应的目标函数的解，生成二维核磁共振谱。算法迭代过程中利用了前面两次迭代的结果，可以提高计算精度，并且加快了算法的收敛速度。

技术领域

本发明涉及油气藏勘探技术领域，尤其涉及一种基于TwIST的二维核磁共振谱的反演方法和装置。

背景技术

随着油气藏勘探开发技术的发展，可以采用核磁共振测井进行储层岩石信息的探测。例如，可以采用核磁共振测井直接探测储层岩石孔隙中流体的氢核信息。核磁共振测井可以获取到回波数据；然后需要对回波数据进行反演，才能得到反映地层流体信息的二维核磁共振谱，该谱也称为弛豫谱，回波数据也称为核磁共振数据。

现有技术中，在对回波数据进行反演，以得到二维核磁共振谱的时候，可以采用L1正则化方法。L1正则化方法可以约束核磁共振数据反演结果的稀疏性，能得到分辨率更高的谱，具体的，可以采用原始-对偶混合梯度法，去求解核磁共振数据反演中的L1正则化问题，并利用广义交叉验证的方法求取正则化参数，最后得到二维核磁共振谱。

然而现有技术中的方法在求解L1正则化问题时，上述这些L1正则化方法在求解L1正则化问题时，每次迭代的结果只和前一次迭代的结果有关，算法收敛的速度较低，进而得到二维核磁共振谱的速度较慢。

发明内容

本发明提供一种基于TwIST的二维核磁共振谱的反演方法和装置，用以解决算法收敛的速度较低，进而得到二维核磁共振谱的速度较慢的问题。

一方面，本发明提供一种基于TwIST的二维核磁共振谱的反演方法，包括：

根据进行N次核磁共振测量时得到的N组核磁共振测量参数，生成核矩阵，其中，N为正整数；

根据所述核矩阵，建立用于反演二维核磁共振谱的目标函数，其中，所述目标函数中包括正则化参数；

在每一次对所述正则化参数进行取值时，采用TwIST算法，求取所述目标函数的解，其中，所述目标函数的解与所述正则化参数每一次的取值是一一对应的；

根据所述目标函数的每一个解，确定与每一个解对应的数据比值，其中，所述数据比值为解所对应的残差与核磁共振测量得到的回波数据的噪声方差之间的比值；

根据每一个所述数据比值、以及所述正则化参数每一次的取值，建立曲线关系，所述曲线关系为所述数据比值的对数与所述正则化参数的对数的曲线关系；

计算所述曲线关系中与所述正则化参数每一次的取值对应的斜率，并根据所述斜率，确定正则化参数的最优取值；

确定与所述最优取值对应的所述目标函数的解，并根据与所述最优取值对应的所述目标函数的解，生成二维核磁共振谱。

另一方面，本发明提供一种基于TwIST的二维核磁共振谱的反演装置，包括：

第一生成模块，用于根据进行N次核磁共振测量时得到的N组核磁共振测量参数，生成核矩阵，其中，N为正整数；

第一建立模块，用于根据所述核矩阵，建立用于反演二维核磁共振谱的目标函数，其中，所述目标函数中包括正则化参数；

求取模块，用于在每一次对所述正则化参数进行取值时，采用TwIST算法，求取所述目标函数的解，其中，所述目标函数的解与所述正则化参数每一次的取值是一一对应的；