[发明专利]多维核磁共振测量方法

权利要求书

1.一种多维核磁共振测量方法，所述方法用于分析材料孔隙结构和内部填充流体特性，其特征在于，所述方法包括：

(1)多维核磁共振脉冲序列的设计及多维核磁共振数据采集；

(2)对多维核磁共振数据进行反演和解释；其中，

所述步骤(1)包括如下步骤：

1.1、在TRS通道向被测样品施加90°射频脉冲将宏观磁化强度矢量M₀扳转90°；

1.2、在GRD通道施加一个持续时间为δ的梯度脉冲；

1.3、随后在TRS通道施加180°射频脉冲后，在GRD通道施加第二个持续时间为δ的梯度脉冲，其中步骤1.3中的梯度脉冲起始时间与步骤1.2中的梯度脉冲起始相差Δ时间；

1.4、等待前两个射频脉冲间隔相等的时间后，向被测样品施加90°射频脉冲，将一定量的磁化矢量从垂直静磁场B₀相同的方向扳转90°；

1.5、施加小角度射频脉冲，在ACQ通道采集自由衰减信号；

1.6、等待时间t_a后，重复施加小角度射频脉冲，会在ACQ通道重复采集到自由衰减信号；

1.7、改变GRD通道中的两个梯度脉冲的强度值，重复所述步骤1.4-1.6以分别采集不同脉冲强度下产生的自由衰减信号。

2.根据权利要求1所述的一种多维核磁共振测量方法，其特征在于，所述步骤1.7得到的自由衰减信号磁化矢量矩阵为：

M(g，nt_a，mt_s)＝∫∫∫K₁K₂K₃F(D，T₁，Δχ)dD·dT₁·dΔχ

其中，g为静磁场梯度值，n为180°脉冲个数，t_a为等待时间，m为FID采集点数，t_s为FID采集点时间间隔，F(D，T₁，Δχ)为被测样品的三维D-T₁-Δχ特性矩阵，T₁为纵向弛豫时间，D为被测样品自扩散系数，Δχ为被测样品与内部填充流体的磁化系数差异，K₁，K₂，K₃为三个核函数，其具体形式为：

K₁＝exp(-Dγ²g²δ²(Δ-δ/3))

K₂＝exp(-nt_a/T₁)

K₃＝exp(-γ·Δχ ·B₀·mt_s)

其中，γ为质子的磁旋比，B₀为静磁场强度。

3.根据权利要求1所述的一种多维核磁共振测量方法，其特征在于，所述步骤(2)采用快速三维数据处理算法对多维核磁共振数据进行反演。

4.根据权利要求2所述的一种多维核磁共振测量方法，其特征在于，所述步骤(2)具体为：

2.1、采用数学张量乘积将前两个核矩阵函数K₁和K₂耦合为一个新的核函数矩阵K₁₂：

2.2、将测量得到的回波串信号磁化矢量矩阵重新表达为：

M＝K₁₂FK₃

2.3、对上述核函数矩阵K₁₂和K₃进行SVD分解及奇异值截取，进而对采集数据进行压缩处理，对核函数矩阵进行奇异值分解可得：

K₁₂＝U₁₂·S₁₂·V′₁₂

K₃＝U₃·S₃·V′₃

其中S₁₂和S₃对角线元素值从大到小排列，且为对角矩阵，大小分别为s₁₂×s₁₂和s₃×s₃，其中s₁₂为K₁₂非零奇异值个数，s₃为K₃非零奇异值个数；U₁₂、V₁₂和U₃、V₃为正交单位阵；对对角矩阵S₁₂和S₃进行截取，使得核函数矩阵的条件数满足设定值C，即：

假设C为1000；和分别对应K₁₂和K₃最大的奇异值，即对角矩阵S₁₂和S₃的第一个对角线元素，表示K₁₂的第i个奇异值，表示K₃的第j个奇异值；

2.4：用截取之后的奇异值分解的单位矩阵对回波串信号磁化矢量矩阵M进行压缩，降低数据内存，压缩后的磁化矢量为：

其中，为压缩后的磁化矢量矩阵，为两个核函数K₁₂和K₃经过SVD分解和奇异值截取后的残余矩阵，U′₁₂、V′₁₂分别为K₁₂经过SVD分解和奇异值截取后的矩阵U₁₂、V₁₂的转置，V′₃为K₃经过SVD分解和奇异值截取后的矩阵V₃的转置；

2.5：数据压缩完成之后，采用Tikhonov正则化方法对数据矩阵进行反演，其中，正则化项为：

其中，α是正则化因子，与采集数据的信噪比相关，K为反演过程中的核矩阵的一般形式，||·||项代表矩阵的Frobenius范数；

2.6：通过选取最优的正则化因子α，得到最终解f_r，求解公式如下：

其中，为压缩截取后残余矩阵的张量积矩阵，为的转置矩阵，α_opt为最优正则化因子，I是单元矩阵，为压缩后的磁化矢量矩阵的一维矩阵形式，其中，