[发明专利]一种核磁共振测井弛豫反演方法

权利要求书

1.一种核磁共振测井弛豫反演方法，其特征在于，包括如下步骤：

S1，利用高斯概率分布建立大孔占优双峰模型，设置核磁共振T2谱模型及模拟采集回波串参数，正演得到衰减的回波串；

S2，采用吉洪诺夫正则化方法，结合广义交叉校验准则，利用最小二乘拟合使得病态方程求解最小化；

S3，设置核磁共振回波间隔、横向弛豫时间T2谱布点方式，加入不同白噪声，采用基于广义交叉校验的吉洪诺夫方法对回波串进行反演得到T2谱；

所述步骤S2中，采用吉洪诺夫正则化方法，结合广义交叉校验广义交叉校验准则，具体步骤如下：

将解的极小范数作为附加条件代入到病态方程中，此时相当于原来的问题中加入约束条件下的最小问题，原问题则可转化为：

正则化方法的实质是在解的范数最小和数据的拟合程度之间达到最优折衷解；关于正则参数的选择，为了兼顾算法的收敛与速度，以及收敛问题的真实解，采用基于数据误差后验取法中的广义交叉较验准则，是对算子方程Ax＝y的离散形式给出的，

其中A＝(a_ij)_m×n；y＝(y₁,y₂,Λ,y_m)^T；

2.如权利要求1所述的一种核磁共振测井弛豫反演方法，所述步骤S1中，设置核磁共振T2谱模型及模拟采集回波串参数，具体步骤如下：

通过高斯概率分布建立双孔模型，设置预先要采集的回波串个数，简化的回波串采集模式为：

式中，M(t_i)为第i个回波的幅度值(i＝1,K,n)；

t_i为n个同样回波间隔的递减时间；

ε_i为第n个回波的噪声；

m为事先的横向弛豫时间；

T₂为横坐标上的对数分布；

a(T_j)为病态方程解出来的弛豫时间T_j对应的T₂幅度。

3.根据权利要求1所述的一种核磁共振测井弛豫反演方法，其特征在于，所述步骤S3中，对回波串进行反演，具体步骤如下：

采用对数均匀布点，个数为32，分布0.5-2048ms，回波间隔：TE＝0.6ms，之前采集的回波个数为1024个，设置4种不同的信噪比噪声，分别为100、30、20、10，反演的方程离散化后为：