CN104819923B - 基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法 - Google Patents

Abstract

本发明公开了一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法，包括以下步骤：第一、挑选基础样品；第二、划分孔隙结构类型；第三、构建识别孔隙结构的测井图版；第四、分孔隙结构类型，建立孔喉分布反演公式；第五、验证公式的正确性；第六、利用测井图版识别井筒剖面的孔隙结构类型；第七、利用核磁共振测井T₂谱定量反演自由流体孔喉分布。本发明的有益效果如下：提供定量评价孔隙结构的一种可靠方法和技术，打破现行方法的研究思想。沿井眼剖面，通过核磁共振测井反演出连续的自由流体孔喉分布，为低渗透砂岩储层有效性评价提供了直接依据，也是测井信息用于定量反演储层微观孔隙结构的有益探索，推动了核磁测井技术的应用与发展。

Description

基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法

技术领域

本发明涉及岩心核磁共振、压汞实验及核磁共振测井领域，特别涉及一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法。

背景技术

核磁共振的基础是原子核的磁性及其外加磁场的相互作用。在外加磁场作用下使得地层流体中的质子极化，核自旋空间取向从无序变为有序。极化后，关闭外加磁场，测量自旋系统从不平衡状态向平衡状态过程的时间(弛豫时间)，即测量纵向弛豫时间T₁、横向弛豫时间T₂。

磁共振测井的原始数据是随时间衰减的自旋回波串，自旋回波串包含了储层物性、孔隙类型、孔径大小、流体类型及其分布等十分丰富的信息。自旋一回波串可以表示为多个衰减指数的迭加：

式中：Echo(t)，核磁测井测量得到的随时间t变化的回波幅度；

T_2i，第i个组分的横向弛豫时间，i＝l，…，N；

ф_i，对应于T_2i组分的核磁共振孔隙度，i＝l，…，N；

通过对核磁测井仪器测量得到的自旋回波串(即式(1))进行多指数拟合，可以求出每一测量点的T_2i-ф_i的分布曲线，即所谓的核磁共振测井T₂分布谱。该T₂分布谱的横坐标为T₂值，纵坐标为信号相对幅度，总的信号幅度与岩石孔隙度有关。

在岩石孔隙中，流体分子弛豫时间不仅与本身弛豫特点有关，还与岩石孔隙结构、成分关系密切，弛豫过程较复杂。目前普遍认为有三种机制控制孔隙流体的横向弛豫过程，即体积弛豫、表面弛豫和扩散弛豫，分别用T_2B、T_2S、T_2D表示相应的弛豫时间常数。进而横向弛豫时间表示为：

式中：T_2B为流体的体积(自由)弛豫时间，ms；

T_2S为横向表面弛豫时间，ms；

T_2D为横向扩散弛豫时间，ms；

D为扩散系数,μm2/ms；

G为磁场梯度，G/cm；

T_E为回波间隔，ms；

S为孔隙表面积，cm²；

V为孔隙体积，cm³；

ρ₂为岩石的横向表面弛豫强度，μm/ms；

T_2B的数值通常为2～3s，比T₂大得多。即T_2B>>T₂。因此式3中右边的第一项可忽略。当磁场均匀时(相对G很小)，且T_E足够短时，式2右边第三项也可忽略，于是T₂为：

得到T₂与孔径r_c的关系式为：

式中：F_S为形状几何因子。对于球形孔隙，F_S＝3,；对柱状管道，F_S＝2。由式5可知，孔隙内的流体的弛豫时间与孔隙空间的大小及形状有关，孔隙越小，比表面积越大，表面相互作用的影响就越强烈，T₂时间也越短。观测的弛豫时间T₂和平均孔径r_c是一一对应的，因此，可以利用T₂谱来评价岩石孔隙大小及其孔喉分布。

目前现有技术提出横向驰豫时间和毛管压力之间的转换关系，如式5，该方法假设孔隙结构可以简化为球形和柱状管道。其比表面积与孔径成线性关系，Pc＝C/T₂线性关系进行转换获得毛细管压力曲线，其转换系数C由相似原理确定，该方法的一个关键假设条件是要求核磁T₂谱曲线与毛细管压力曲线的形态一致性高，而从大量的岩心核磁共振测量与压汞测量结果对比发现，T₂谱与压汞微分曲线形态存在一定差异。并且，未考虑T₂谱与压汞曲线反应的孔隙空间不同这一关键。故而，这种方法计算误差就会较大。

在上述的研究基础上，基于最大相似性原理的方法确定T₂与Pc之间的转换系数C。首先假设一个C，将与P_c～S_Hg两条曲线重合在一张图上(图1)。通过调节C值来移动曲线。明显存在唯一值C，是使与P_c～S_Hg的相关性最大。去其相关性最大时的C值为T₂与Pc之间的转换系数，进而将核磁T₂谱转换为压汞微分曲线。这一方法的计算精度不高。原因在于忽略了孔喉空间中薄膜束缚水对小孔核磁信号的影响。即核磁T₂谱测量过程中，大孔喉空间薄膜束缚水信号叠加到小孔核磁信号上，使得小孔信号增强。因此，核磁T₂谱与压汞分布曲线在小孔后部分不成线性关系。故而，在小孔吼部分转换系数并不适用。

现有技术还提出了2种核磁共振分布评价毛细管压力曲线的改进方法。

第1种方法：岩心核磁T₂分布与相应的压汞毛管半径分布进行的相关分析。发现在大孔、小孔处两者呈现明显不同两种幂函数关系。所以，在T₂分布为双峰时需要用不同幂函数分段构造毛管压力曲线对T₂分布呈单峰分布时，采用单一幂函数来构建伪毛管压力曲线。

第2种方法：何雨丹从压汞、NMR核磁所反映的岩石孔隙空间入手，提出薄膜束缚水对核磁T₂谱的影响，并提出与压汞信息对应的是核磁自由流体谱而非核磁饱和谱。故而，消除薄膜束缚水影响，利用自由流体分布构造毛细管压力曲线。

该技术中方法2中利用转换系数C将自由流体谱转换为伪毛管压力曲线，未考虑大小孔喉的转换关系不同这一特征。

方法1是对方法2的改进，利用不同幂函数构建伪毛管压力曲线，与压汞得到的毛细管压力曲线对某些岩样较以前吻合度有了一定的提高，但它还是通过对饱和水T₂谱进行积分所得，其最大饱和度为100％，没有考虑到最大进汞饱和度问题，存在一定的局限性。

为了解决核磁测井定量评价孔隙结构的难题，在现有的研究基础上，即回到核磁原理，直接研究核磁饱和T₂谱与压汞孔喉分布曲线的对应关系，及核磁自由流体T₂谱与压汞孔喉分布≥0.0248μm部分的曲线(<0.0248μm的压汞孔喉分布曲线主要包含未连通孔喉的信息)的对应关系，并在分孔隙结构、分孔喉尺度的情况下构建反演核磁孔喉分布、自由流体孔喉分布公式。从而定量反演核磁孔喉分布和自由流体孔喉分布。提供定量评价孔隙结构的一种可靠方法和技术，打破现行方法的研究思想。推动核磁测技术的应用与发展。

发明内容

本发明针对现有技术的缺陷，提供了一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法，能有效的解决上述现有技术存在的问题。

为了实现以上发明目的，本发明采取的技术方案如下：

一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法，以下步骤：

第一步、挑选基础样品，基础样品的选取依据：

a、样品需以配套核磁实验和压汞实验的数据为依据；

b、样品的常规测井信息能真实反映地层信息；

c、样品范围涵盖超低孔、特低孔、低孔和中孔，非渗、超低渗、特低渗和低渗孔隙结构类型；

第二步、以压汞数据为基础，划分低渗砂岩储层孔隙结构类型，将低渗砂岩储层孔隙结构类型分为五种类型，分为I₁型、I₂型、II₁型、II₂型和III型；

所述I₁型排驱压力为P_d<0.2，最大进汞饱和度为SH_gmax≥80；

所述I₂型排驱压力为0.2≤P_d<0.4，最大进汞饱和度为70≤SH_gmax<80；

所述II₁型排驱压力为0.4≤P_d<0.7，最大进汞饱和度为60≤SH_gmax<70；

所述II₂型排驱压力为0.7≤P_d<1.6，最大进汞饱和度为50≤SH_gmax<60；

所述III型排驱压力为P_d≥1.6，最大进汞饱和度为SH_gmax<50；

第三步、提取基础样品点对应的常规测井信息，常规测井中表现出AC、HRID值越高，而GR值越低的特征，以此特征为依据，构建AC-HRID-GR气泡图图版、雷达图图版识别孔隙结构；

第四步、挑选标准样品，标准样品的选取依据如下：

a、从孔隙结构分类结果中筛选出标准样品；

b、标准样品能代表该种类型孔隙结构样品的特征；

第五步、针对每个标准样品，提取自由流体T₂谱与压汞孔喉分布的对应数据点，在分孔喉尺度<0.1μm、0.1-1μm、1-4μm的基础上，建立自由流体孔喉分布、饱和T₂谱与压汞孔喉分布的对应关系，从而计算核磁孔喉分布、核磁自由流体孔喉分布；

第六步、依据第五步中所构建的对应关系，反演自由流体孔喉分布、核磁孔喉分布，并建立核磁反演孔喉分布直方图，将其与压汞孔喉分布直方图比较、验证；

第七步、在划分孔隙结构的基础上，建立核磁饱和谱与核磁自由流体谱的对应关系，将核磁共振测井T2谱转换为核磁共振自由流体T2谱，进而计算地层自由流体孔喉分布；

第八步、沿井眼剖面定量反演低渗砂岩储层的自由流体孔喉分布。

作为优选，所述第五步中饱和T₂谱与压汞孔喉分布的对应关系，如下式：

r_c＝a₁×T₂+b₁ 1≤r_c<4μm

r_c：压汞孔喉分布，μm；

T₂：与r_c对应的自由流体谱横向弛豫时间，ms。

作为优选，所述步骤七中建立核磁饱和谱与核磁自由流体谱的对应关系，如下式：

I₁、I₂、II₁型孔隙结构：

y＝c₁×x+d₁T_2饱和谱≥40ms

y＝c₂×x²+d₂×x+e₂ 4.5≤T_2饱和谱<40ms

y＝c₃×x²+d₃×x+e₃T_2饱和谱<4.5ms

II₂、III型孔隙结构：

y＝c₄×ln(x)+d₄T_2饱和谱≥2.25ms

y＝c₅×x²+d₅×x+e₄ 0.7≤T_2饱和谱<2.25ms

y＝c₆×x²+d₆×x+e₅T_2饱和谱<0.7ms

y：为核磁自由流体谱分量；

x：核磁饱和谱分量。

与现有技术相比本发明的优点在于：在前人研究基础上，采用新的研究思路与方法，即回到核磁原理，直接研究核磁饱和T₂谱与压汞孔喉分布曲线的对应关系，及核磁自由流体T₂谱与压汞孔喉分布≥0.0248μm部分的曲线(<0.0248μm的压汞孔喉分布曲线主要包含未连通孔喉的信息)的对应关系，并在分孔隙结构、分孔喉尺度的情况下构建反演核磁孔喉分布、自由流体孔喉分布公式。从而定量反演核磁孔喉分布和自由流体孔喉分布。提供定量评价孔隙结构的一种可靠方法和技术，打破现行方法的研究思想。自由流体T₂与压汞孔喉分布对应关系更好，以此建立了不同孔隙结构类型二者之间不同孔喉尺度对应的关系式可在井筒剖面上通过识别孔隙结构类型，进而利用NMR定量反演核磁孔喉分布、自由流体孔喉分布，省去了构建伪毛管曲线环节，为低渗透砂岩储层有效性评价提供了直接依据，也是测井信息用于定量反演储层微观孔隙结构的有益探索，推动核磁测技术的应用与发展。

附图说明

图1为东营凹陷南斜坡沙河街组沙四段低渗致密砂岩压汞数据孔隙结构分类图；

图2为基本样品压汞特征参数排驱压力-物性交会图；

图3为基本样品压汞特征参数SH_gmax-物性交会图；

图4为基本样品压汞特征参数平均孔喉半径-物性交会图；

图5为常规测井“三性”雷达图识别图版；

图6为常规测井AC-GR-HRID气泡图识别图版；

图7为FX井AC-GR-HRID气泡图；

图8为I₁型标准样品核磁T₂谱特征；

图9为I₁型标准样品压汞孔喉分布特征；

图10为I₂型标准样品核磁T₂谱特征；

图11为I₂型标准样品压汞孔喉分布特征；

图12为II₁型标准样品核磁T₂谱特征；

图13为II₁型标准样品压汞孔喉分布特征；

图14为II₂型标准样品核磁T₂谱特征；

图15为II₂型标准样品压汞孔喉分布特征；

图16为III型标准样品核磁T₂谱特征；

图17为III型标准样品压汞孔喉分布特征；

图18为I₁型标准样品核磁自由流体谱与压汞孔喉分布对应图；

图19为I₂型标准样品核磁自由流体谱与压汞孔喉分布对应图；

图20为II₁型标准样品核磁自由流体谱与压汞孔喉分布对应图；

图21为II₂型标准样品核磁自由流体谱与压汞孔喉分布对应图；

图22为III型标准样品核磁自由流体谱与压汞孔喉分布对应图；

图23为I₁型标准样品自由流体孔喉分布反演公式图；

图24为I₂型标准样品自由流体孔喉分布反演公式图；

图25为II₁型标准样品自由流体孔喉分布反演公式图；

图26为II₂型标准样品自由流体孔喉分布反演公式图；

图27为III型标准样品自由流体孔喉分布反演公式图；

图28为I₁型标准样品核磁孔喉分布反演公式图；

图29为I₂型标准样品核磁孔喉分布反演公式图；

图30为II₁型标准样品核磁孔喉分布反演公式图；

图31为II₂型标准样品核磁孔喉分布反演公式图；

图32为III型标准样品核磁孔喉分布反演公式图；

图33为I₁型标准样品岩石孔喉分布反演图；

图34为I₂型标准样品岩石孔喉分布反演图；

图35为II₁型标准样品岩石孔喉分布反演图；

图36为II₂型标准样品岩石孔喉分布反演图；

图37为III型标准样品岩石孔喉分布反演图；

图38为低渗透砂岩储层孔隙结构定量反演实例分析图。

具体实施方式

为使本发明的目的、技术方案及优点更加清楚明白，以下参照附图并举实施例，对本发明做进一步详细说明。

一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法，包括以下步骤：

第一步、挑选基础样品，基本样品的选取依据：

a、样品需用配套核磁实验和压汞实验的数据，

b、样品必须是常规测井信息，

c、样品包括不同孔隙结构类型；

第二步、以压汞数据为基础，划分低渗砂岩储层孔隙结构类型：

由于低渗砂岩孔隙结构复杂，不同孔隙结构样品的压汞曲线特征差异较大。经对比分析，将低渗砂岩储层分为三大类五小类，如图1，分别为I₁型(排驱压力低P_d<0.2，最大进汞饱和度高SH_gmax≥80)；I₂型(排驱压力较低，0.2≤P_d<0.4，最大进汞饱和度较高，70≤SH_gmax<80)；II₁型(排驱压力集中在0.4≤P_d<0.7，最大进汞饱和度分布在60≤SH_gmax<70)；II₂型(排驱压力偏高，0.7≤P_d<1.6，最大进汞饱和度偏低，50≤SH_gmax<60)；III型(排驱压力高，P_d≥1.6，最大进汞饱和度低，SH_gmax<50)。并结合物性资料，建立压汞特征参数-物性参数交会图，样品孔隙度、渗透率与排驱压力、最大进汞饱和度等参数存在良好的对应关系，如图2、3、4。当孔隙结构越好，即孔隙度、渗透率越大时，排驱压力越小，最大进汞饱和度越大。结合核磁T₂谱分析发现，孔隙结构越好的样品束缚水饱和度越低。以此为依据，建立不同孔隙结构物性参数划分范围，见表1。

表1 ES₄层段5种孔隙结构分类表

第三步、建立常规测井识别孔隙结构的方法：

在划分孔隙结构的情况下，提取基本样品点对应的常规测井信息，发现孔隙结构越好，常规测井中表现出AC、HRID值越高，而GR值越低的特征。以此特征为依据，构建AC-HRID-GR气泡图图版、雷达图图版识别孔隙结构，如图5、6，沿井筒剖面识别低渗砂岩孔隙结构类型。

以FX井为例，如图7，3241.25-3243.2m层段为II₁孔隙结构，3243.2-3244.8m层段为在II₂孔隙结构，3246.13-3248.7m层段为II₁孔隙结构，3249.63-3251.8m层段为I₂孔隙结构。评价准则见表2。

表2常规测井识别孔隙结构准则表

第四步、挑选标准样品：

(1)基本样品的选取依据：

a、每种标准样品可代表该类孔隙结构类型；

b、孔隙结构特征明显；

5种标准样品特征描述：

a、I₁型孔隙结构

连通性较好，排驱压力值低，一般<0.2MPa，最大进汞饱和度值高，往往>80％；核磁自由流体T₂谱主要集中在10-300ms之间(图8)；孔喉分布集中在0.73-1.495μm范围内(图9)。其CPOR＝17.1％，CPERM＝10.991mD，物性级别属于中孔低渗；此孔隙结构层段具有良好的渗透性。

b、I₂型孔隙结构

连通性较差，P_d较低，在0.2-0.4MPa，SH_gmax分布于70-80％；；核磁自由流体T₂谱主要集中在10-300ms之间，但其含量低于I₁型(图10)；孔喉分布集中在0.1374-0.3721μm范围内(图11)；CPOR＝17.12％，CPERM＝8.233mD，属于低孔特低渗结构，此孔隙结构层段的渗透性较好。

c、II₁型孔隙结构

P_d分布于0.4-0.7MPa，SH_gmax分布于60-70％；核磁自由流体T₂谱主要集中在5-300ms之间(见图12)；孔喉分布集中在0.1497-0.3801μm范围内，但含量低于I₂型(图13)。CPOR＝11.46％，CPERM＝0.3627mD，属于特低孔超低渗结构。

d、II₂型孔隙结构

P_d分布于0.7-1.6MPa，SH_gmax分布于50-60％；核磁自由流体T₂谱分布广，但自由孔隙含量低(见图14)；孔喉分布在0.0936-0.1477μm范围内有一单峰，但其微孔含量高(图15)。CPOR＝10％，CPERM＝0.168mD，属属于特低孔超低渗结构，此孔隙结构层段的渗透性较差，为差渗透层。

e、III型孔隙结构

P_d值高，往往>1.6MPa，SH_gmax值低，一般<50％；T₂谱可见束缚孔隙度值大，而有效孔隙度值小，S_BVI则>75％。(见图16)，孔喉分布以微孔(<0.0248μm)分布为主(见图17)；CPOR＝5.87％，CPERM＝0.0141mD，属于超低孔非渗结构，此孔隙结构层段的渗透性差，为典型非渗层。

第五步、分孔隙结构类型、孔喉尺度构建反演公式

针对每个标准样品，提取自由流体T₂谱与压汞孔喉分布的对应数据点，如图18、19、20、21、22，在分孔喉尺度(<0.1μm、0.1-1μm、1-4μm)的基础上，如图23、24、25、26、27，自由流体孔喉分布反演公式见表3，同时建立饱和T₂谱与压汞孔喉分布的对应关系，如图28、29、30、31、32，从而建立核磁孔喉分布反演公式见表4。

表3 5种孔隙结构自由流体谱与压汞孔喉分布的对应关系表

表4 5种孔隙结构核磁饱和谱与压汞孔喉分布的对应关系

第六步、模型验证

依据所构建的反演公式，反演自由流体孔喉分布、核磁孔喉分布，并建立核磁反演孔喉分布直方图，将其与压汞孔喉分布直方图比较、验证，如图33、34、35、36、37。

第七步、核磁自由流体谱反演

实际核磁测井不能直接提取核磁共振自由流体谱。故而，在分孔隙结构的基础上，建立核磁饱和谱与核磁自由流体谱的对应关系(如表5)。进而可利用核磁共振测井T2谱计算地层自由流体孔喉分布。

表5 5种孔隙结构核磁饱和谱与自由流体谱的对应关系

(注：y为自由流体谱分量；x为饱和谱分量)

第八步、沿井眼剖面定量反演低渗砂岩储层的自由流体孔喉分布。

以FX井为例，在常规测井信息沿井筒剖面识别孔隙结构的基础上，运用相应孔隙结构的反演公式将核磁共振测井T₂谱沿井筒剖面方向定量反演为岩石孔隙空间分布，如图38。分析FX井核磁共振测井T₂谱，发现3238.125-3241.2m、3244.8-3246.13m、3248.7-3249.625m层段T₂谱位于最左边，明显孔隙结构差，为III型。3241.25-3243.2m(87号)、3246.13-3248.7m(89号)T₂谱分布范围广，明显双峰特征为II1型，解释结论为差油层，且在核磁T₂谱中高值分布含量高。3243.2-3244.8m(88号)T₂谱单峰分布，为II₂型，解释结论为油干层。3249.63-3251.8m(90号)T₂谱分布范围广，并且分布在T₂高值范围，为I₂型，解释结论为油层。与图6AC-GR-HRID识别图结果一致，说明核磁T₂谱识别孔隙结构和有效储层的适用性良好。进而并沿井筒剖面，通过核磁共振T₂谱反演的核磁自由流体谱来计定量反演FX井核磁孔喉分布和自由流体孔喉分布。

需要注意的是，针对本技术方案，可以通过别的替代方案同样能完成发明目的，

如研究对象转换，可将低渗砂岩换为其他岩性，如砂岩、泥页岩、碳酸盐岩；

孔隙结构类型划分可利用薄片、电镜、CT扫面等其他资料建立；

井筒剖面孔隙结构分类可采用时频分析方法或直接用核磁测井进行分类。

本领域的普通技术人员将会意识到，这里所述的实施例是为了帮助读者理解本发明的实施方法，应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合，这些变形和组合仍然在本发明的保护范围内。

Claims (3)

1.一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法，其特征在于包括以下步骤：

第一步、挑选基础样品，基础样品的选取依据：

a、样品需以配套核磁实验和压汞实验的数据为依据；

b、样品的常规测井信息能真实反映地层信息；

c、样品范围涵盖超低孔、特低孔、低孔和中孔，非渗、超低渗、特低渗和低渗孔隙结构类型；

第二步、以压汞数据为基础，划分低渗砂岩储层孔隙结构类型，将低渗砂岩储层孔隙结构类型分为五种类型，分为I₁型、I₂型、II₁型、II₂型和III型；

所述I₁型排驱压力为P_d<0.2，最大进汞饱和度为SH_gmax≥80；

所述I₂型排驱压力为0.2≤P_d<0.4，最大进汞饱和度为70≤SH_gmax<80；

所述II₁型排驱压力为0.4≤P_d<0.7，最大进汞饱和度为60≤SH_gmax<70；

所述II₂型排驱压力为0.7≤P_d<1.6，最大进汞饱和度为50≤SH_gmax<60；

所述III型排驱压力为P_d≥1.6，最大进汞饱和度为SH_gmax<50；

第三步、提取基础样品点对应的常规测井信息，常规测井中表现出AC、HRID值越高，而GR值越低的特征，以此特征为依据，构建AC-HRID-GR气泡图图版、雷达图图版识别孔隙结构；

第四步、挑选标准样品，标准样品的选取依据如下：

a、从孔隙结构分类结果中筛选出标准样品；

b、标准样品能代表该种类型孔隙结构样品的特征；

第五步、针对每个标准样品，提取自由流体T₂谱与压汞孔喉分布的对应数据点，在分孔喉尺度<0.1μm、0.1-1μm、1-4μm的基础上，建立自由流体孔喉分布、饱和T₂谱与压汞孔喉分布的对应关系，从而计算核磁孔喉分布、核磁自由流体孔喉分布；

第六步、依据第五步中所构建的对应关系，反演自由流体孔喉分布、核磁孔喉分布，并建立核磁反演孔喉分布直方图，将其与压汞孔喉分布直方图比较、验证；

第七步、在划分孔隙结构的基础上，建立核磁饱和谱与核磁自由流体谱的对应关系，将核磁共振测井T2谱转换为核磁共振自由流体T2谱，进而计算地层自由流体孔喉分布；

第八步、沿井眼剖面定量反演低渗砂岩储层的自由流体孔喉分布。

2.根据权利要求1的一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法，其特征在于：所述第五步中饱和T₂谱与压汞孔喉分布的对应关系，如下式：

r_c＝a₁×T₂+b₁ 1≤r_c<4μm

<mrow> <msub> <mi>r</mi> <mi>c</mi> </msub> <mo>=</mo> <msub> <mi>a</mi> <mn>2</mn> </msub> <mo>&amp;times;</mo> <msup> <msub> <mi>T</mi> <mn>2</mn> </msub> <msub> <mi>b</mi> <mn>2</mn> </msub> </msup> <mo>,</mo> <mn>0.1</mn> <mo>&amp;le;</mo> <msub> <mi>r</mi> <mi>c</mi> </msub> <mo>&lt;</mo> <mn>1</mn> <mi>&amp;mu;</mi> <mi>m</mi> </mrow>

<mrow> <msub> <mi>r</mi> <mi>c</mi> </msub> <mo>=</mo> <msub> <mi>a</mi> <mn>3</mn> </msub> <mo>&amp;times;</mo> <msup> <msub> <mi>T</mi> <mn>2</mn> </msub> <msub> <mi>b</mi> <mn>3</mn> </msub> </msup> <mo>,</mo> <msub> <mi>r</mi> <mi>c</mi> </msub> <mo>&lt;</mo> <mn>0.1</mn> <mi>&amp;mu;</mi> <mi>m</mi> </mrow>

r_c：压汞孔喉分布，μm；

T₂：与r_c对应的自由流体谱横向弛豫时间，ms。

3.根据权利要求2的一种基于核磁共振的低渗透砂岩储层孔隙结构定量反演方法，其特征在于：所述步骤七中建立核磁饱和谱与核磁自由流体谱的对应关系，如下式：

I₁、I₂、II₁型孔隙结构：

y＝c₁×x+d₁ T_2饱和谱≥40ms

y＝c₂×x²+d₂×x+e₂ 4.5≤T_2饱和谱<40ms

y＝c₃×x²+d₃×x+e₃ T_2饱和谱<4.5ms

II₂、III型孔隙结构：

y＝c₄×ln(x)+d₄ T_2饱和谱≥2.25ms

y＝c₅×x²+d₅×x+e₄ 0.7≤T_2饱和谱<2.25ms

y＝c₆×x²+d₆×x+e₅ T_2饱和谱<0.7ms

y：为核磁自由流体谱分量；

x：核磁饱和谱分量。