南海西部珠江口盆地珠江组一段发育砂泥岩薄互层，这些地层属海相沉积，分布范围广，勘探潜力大^[1]。由于薄层厚度呈现厘米级甚至更薄，远低于常规测井仪器1 ~ 4 ft（0.3 ~ 1.2 m）的垂向分辨率，来源于薄层的综合响应信号将被围岩物理性质响应所钝化，大部分响应反映地层的围岩物理性质^[2]，测井资料表现为低电阻率、高自然伽马等低阻油层特征^[3]。单个砂层厚度越薄测井响应值偏离真实值的差异越大，所以这样的油气储层极可能因没有被正确识别而被忽略，而其测井评价长期以来也面临着巨大挑战。

在过去几十年，针对高分辨率测量研究出几种方法，例如垂向分辨率高达0.2 in（0.5 cm）的地层微电阻率成像（FMI）、先进的砂泥岩薄互层处理解释模型（LSA）、高分辨率分析和岩石学参数重构等。这些方法通常需要高分辨率测量的辅助，并且在地层建模或测井评价中需要花费较长时间。由于海上作业的特点，需要在裸眼测井完成之后不久，对目标层段进行快速合理的评价以便做出正确及时的决策，前述这些方法往往不能够满足现场快速决策的需求。

1 砂泥岩薄互层的岩石特性与核磁响应机理

核磁共振T2谱反映了地层的孔隙结构，在仪器测量体积范围内T2值的大小与岩石孔径大小相关。对于均质的碎屑岩，典型的T2谱是随岩石孔隙大小而连续分布的。对于砂泥薄互层，通常包含有两组显著的大小孔隙分布：黏土水、束缚流体等微孔隙主要存在于泥岩薄层中，而大孔隙仅存于砂岩薄层中，大孔隙填充的主要是自由流体。因此在核磁测井仪器探测体积内，来自于储层砂泥岩薄互层的孔隙信息表现值在T2分布谱上通常为比较明显的双峰分布。

Cao Minh等人在2006年展示了砂泥岩薄互层的核磁岩石物理特征并用地面试验模拟的方法得到了双峰谱。实验显示，核磁测量是仪器垂向分辨率范围内，在仪器测量空间内所有测量信号的线性累加。由于泥岩比砂岩通常具有更小的孔隙，T2谱的分布将呈现双峰分布的特征，核磁响应的重点是它与泥岩薄层的几何形态无关，只要总的体积部分不变化，2个或3个泥岩层分层或连接到一起，将产生相同的T2双峰分布^[4]。

图 1是中国南海西部的一个实例，显示了在砂泥岩薄互层中核磁T2谱的岩石物理特征。在最右边道的微电阻率成像给出了精细的砂泥岩薄互层直观图；从左边数第三道核磁T2谱也清楚地呈现出双峰现象。高的T2峰态反应的是砂岩纹层，低的反应的是泥质或粉砂纹层。应指出的是薄层厚约1 in（2.5 cm）或者更薄；体积密度和中子孔隙度测井曲线不能够分辨出这些薄层地层。

在文献中，也给出了计算砂岩含量和薄层砂岩内在渗透率（K_sand）的公式^[4]。在砂泥岩薄互层条件下，用于现场砂岩比（F_sand）计算方法简化公式如下：

(1)

(2)

式中：Φ_t —核磁测量的总孔隙度；

Φ_sand—划分的砂岩核磁孔隙度；

K_t—地层的真正有效渗透率；

泥岩被认为是不渗透的。在文献中，K_t用来自于CPMG回波总数方法计算的核磁渗透率来替代。

2 渗透率计算方法研究

尽管F_sand和K_sand的计算在文献中已给出，精确的核磁砂岩含量和渗透率仍不能在薄储层直观显现。目前常用的核磁仪器的垂向分辨率是18 in（45.72 cm）或更低。尽管过去研发的核磁测井的高分辨率方法以及对砂泥岩薄互储层的评价得到了改善，这种级别的垂向分辨率仍不足于直接表现薄层的特征^[5]。因此，对砂泥岩薄互层来说，需要建立一种新的渗透率模型评价地层物理特性。

2.1 薄层渗透率测量误差分析

有效渗透率与测量尺度相关，不同测井工具和不同方法由于自身的垂向分辨率和探测深度不同，对同一段薄互层地层，会得出不同的测量结果和渗透率计算结果。如图 2所示，建立一个砂岩比F_sand=1 : 2，砂岩自身渗透率K_sand=1 μm²的砂泥岩薄互层正演模型^[6]。假定砂泥岩都是均质的，根据达西定律，在考察（测量）空间内的地层平均渗透率（K_b）和实际有效渗透率都是0.5 μm²。假定砂岩和泥岩薄层都具有等同的总核磁孔隙度（Φ_t），所有的束缚流体仅存于泥岩中，用隐含参数和Timur-Coates公式进行核磁渗透率计算^[7]，则得到的核磁视渗透率（KTIM）为0.08 μm²，计算得到的相对误差为（500-80）/80=525%。如果砂岩比例更低，F_sand=1 : 3，则误差变得更大，相对误差达到（333-20）/20=1 565%。在砂泥薄互层这种情况下，利用常用公式计算的核磁视渗透率比实际地层的真实渗透率就大大降低了，从而低估了砂泥岩薄互储层的质量和产能，导致有经济价值的层段被忽略。

2.2 薄互层鉴别

利用核磁测量识别薄互层是基于这样的逻辑进行计算的：有生产能力的砂泥岩薄互层通常是由分开的大孔隙（砂岩纹层）和小孔隙（泥岩纹层）组成，这就导致T2谱分布的双峰形态。

在研究中首先建立了一个岩石类型指数SQI，以用来进一步区分地层是薄互层或常规层。由核磁T2谱分布形态计算产生的SQI定义为高于某个截止值的大孔隙与总孔隙体积比的函数。薄互层区分指示线SQI_disc，从核磁视渗透率KTIM计算得到，计算方法如下；式中的视渗透率KTIM由Timur-Coates渗透率计算公式^[7]所得。

(3)

式中：A、B分别为单独的偏差量和系数，它们可以从SQI-KTIM交会图推断所得；对于多数现场研究例子，简单估算的隐含值为0和0.33；图 3中在常规砂岩和泥质砂岩中指示出砂泥岩薄互层，高的SQI和低的KTIM反映了高概率的砂泥岩互层。如果SQI大于SQI_disc，岩石类型认为是互层的，相反则认为是常规砂泥岩储层。

2.3 渗透率计算

如前所述，由于采用的是均质砂岩模型，在砂泥岩薄互层中测井计算的视渗透率受泥质纹层影响很大；同时由于受测井仪器分辨率的局限，如何确定砂岩纹层的真实渗透率进而准确计算砂泥岩薄互层有效渗透率是一个难题。相关文献^[4-8]中给出的计算方法都无法真正排除泥质纹层的影响。如果认为岩石类型是砂泥互层，比较准确的砂层纹层内在渗透率可以通过修正的SDR公式计算所得。

(4)

式中：T _{2LM, sand}为与核磁砂岩划分相关的T2对数平均值，Φ_{t, sand}是与薄层砂岩相关的100%砂岩总孔隙度，对于SDR渗透率转化的隐含常数是：a=4，m=2，n=4。需要指出的是，T _{2LM, sand}并不是通常意义的整个核磁T2谱的对数平均值，而是被认为是砂岩响应部分的T2谱峰的对数平均值。这可以通过进一步分析核磁测井T2谱得到，但现场测井通常无法直接提供。

本文中阐述的方法开发了一种相对简单易行的用以评价砂岩纹层内在渗透率（K_sand）的替代方法。从油田实际生产角度来说，在砂泥岩薄互层中真正感兴趣的仅仅是高渗透的有经济性的砂岩纹层。在研究的目标油田高渗透砂岩中，大孔隙和总孔隙体积比Ratio_L2T与SQI呈典型的线性关系^[9]，如图 4所示。因此，通过SQI和图 4确定的参数可以计算Ratio_L2T；利用Ratio_L2T对核磁测量的薄互层T2谱中的砂岩纹层孔隙度进行校正，得到砂岩纹层的总孔隙Φ_sand（式5 ~式7）；结合式1：F_sand=Φ_sand/Φ_t，可以进一步计算F_sand。

(5)

式5中，c=0.9147，d=0.078 1，校正的Φ_sand计算如下：

(6)

式6中的cutoff1是一个核磁测井T2截止值参数，用以划分出砂泥岩薄互层核磁T2谱中砂岩/泥岩的双峰部分，结合图 3，对砂岩纹层的渗透率进行校正和估算：

(7)

式7中A和B即式3中的偏移量A和系数B，由SQI-KTIM交会确定。

3 现场实例应用

南海西部珠江口盆地珠江组一段富集砂泥岩薄互层沉积。多年来，已经采用多种手段，诸如核磁、岩石物理分析、地层测试（MDT）、钻杆测试（DST）等方法对砂泥岩薄互层进行了测量和渗透率解释，所得结果缺乏一致性。尤其是基于核磁渗透率的产能估计比地层测试（MDT）和钻杆测试（DST）结果要低得多。原因就是对于砂泥岩薄互层用传统核磁渗透率转换而来的渗透率计算产能，会带来较大的误差。现方法应用于被常规评价方法忽视的砂泥岩薄互层产层，这种方法输出的薄互层砂岩比F_sand、体积平均有效孔隙度CMFF、体积平均渗透率K_bulk是评价油藏和储量计算的主要岩石特性参数，最终输出成果实例如图 5。图 5中，道1为深度道，道2 ~道4为常规测井资料，道5为薄互层标志，道6为地层砂岩比（纯砂岩层相对于薄互层的总体厚度比率），道7为体积平均渗透率，道8为薄互层识别置信度，道9为FMI成像；从这个图中可以很好的识别砂泥岩薄互层。新方法计算的薄互层平均渗透率比Timur-Coates公式计算的渗透率明显要高要好。

通过现场研究对本方法和常规方法计算的核磁渗透率进行对比，从地层测试MDT解释的有效渗透率用来刻度点，很明显与地层测试MDT结论相比的新方法所得的渗透率比标准的核磁转换所得的渗透率要好得多。图 6左图显示的是用Timur-Coates公式计算的核磁渗透率与MDT直接测量的渗透率的相关对比，右图为新方法得到的核磁渗透率与MDT直接测量的渗透率相关对比。从现场实例可知，新方法得到的核磁渗透率的精确度有明显的改善（图 6）。

4 结论

论述了一种利用核磁测量的T2谱来探测和评价砂泥岩薄互层的岩石物理属性的新方法。它利用在砂泥岩薄互层中核磁响应的双峰分布特征的优点来克服常规测井仪器垂向分辨率测量的局限。

这种方法定义了一些有用的参数，例如岩石类型、砂岩纹层比例F_sand、体积平均渗透率K_bulk等。它是基于仪器探测范围内薄互层的整体响应来设定和计算这些参数，而不是为了寻求解决某个独立的、单一的砂岩纹层的属性。文中论述的方法和认识可以用于在类似地区识别和评价新井或老井中可能被忽视的砂泥岩薄互层油气藏。

本文描述的新方法通过核磁测井为现场决策提供了快速直观的结论。案例研究证实从新方法所得的渗透率与MDT渗透率关系比标准核磁转换所得渗透率要好得多，砂泥岩薄互层的产能评价得到了显著改善，为确定经济性产层提供了可靠依据。

术语含义：

T2：横向弛豫时间；T_2LM，_sand：T2谱砂岩划分对数平均值；F_sand：砂岩含量；F_clay：黏土含量；F _shale：泥岩含量；Φ _t：总孔隙度；Φ_t，_sand：100%砂岩薄层的总孔隙度；Φ_sand：砂岩孔隙度；Φ_clay：黏土孔隙度；Φ_shale：泥岩孔隙度；BFV：束缚流体体积；KTIM：Timur-Coates渗透率；KSDR：SDR渗透率。