致密滩坝砂储集层孔隙分形特征、预测及应用——以东营凹陷为例

侯庆杰，刘显太，韩宏伟，刘浩杰，魏国华，陈雨茂，于文政，王奇韵

1.中国石油化工股份有限公司胜利油田分公司物探研究院，山东东营 257022

2.胜利石油管理局博士后科研工作站，山东东营 257000

3.中国石油化工股份有限公司胜利油田分公司，山东东营 257000

0 引言

致密油是继页岩气之后全球非常规油气勘探开发的新热点[1-3]，被石油工业界誉为“黑金”[4]。美国Bakken 致密油已实现规模化开发，扭转了北美石油产量递减趋势，改变了全球能源格局，成为全球重要油气接替能源[5]。中国致密油勘探近两年在渤海湾、鄂尔多斯、准噶尔和塔里木等盆地相继取得了一定突破，以胜利油田为例，致密油控制储量高达2.7×108t，是常规油气藏接替能源[6]，其中，致密滩坝砂储量为1.7×108t，占致密油藏储量的63%，为最主要的致密油储层类型。滩坝砂体是发育于陆相湖盆滨浅湖地区非常重要的一种沉积砂体，多分布在湖泊边缘、湖湾、湖中局部隆起的缓坡一侧，是波浪和湖流将各处带来的砂质组分再改造后在有利场所堆积而成[7]。沉积物粒度细，一般为细砂岩、粉砂岩、泥质粉砂岩等[8]，总体具有薄互层特征。由于非均质性强，储量动用难度较大，动用率较低，如何实现在致密滩坝砂储层中“差中选优”是提高致密油储量动用率的关键问题之一。从其他地区致密油藏开发经验来看[9-10]，孔隙结构不仅决定了储集层储集和渗流能力，也影响着最终可动油气总量[11]。

中高渗透储层通常物性较好，开发难度低，储层评价参数通常以孔隙度、渗透率及丰度为主，评价参数简单。与常规砂岩储层相比，致密砂岩储层具有孔喉尺度小、非均质性强和孔隙结构复杂的特点[12-14]，评价时需要加入孔喉半径、分选系数和退汞效率等孔隙结构参数。然而，根据致密储层的特点，仅仅使用有限的样品实验参数表示整个评价层段和区域，忽略致密储层非均质性，使评价结果与实际情况有很大的偏差；如果使用测井资料预测各个评价参数，参数种类繁多导致工作量巨大，均无法综合全面地评价储层优劣。孔隙分形概念的出现在一定程度上解决了这些困难。前人研究表明，砂岩的微观孔隙结构具有很好的统计自相似性[15-16]。与经典的欧式几何相比，利用分形几何能够对孔隙结构的非均质性和复杂性进行更真实地表征，而分形维数即为定量表征孔隙分形程度的参数。孔隙分形与储层物性、孔隙结构的关系已有相关讨论，但是孔隙分形如何影响开发效果，如何利用孔隙分形评价储层优劣还少见讨论。此外，目前对于孔隙分形维数的计算多停留在单点实验测试分析的层面上，无法形成单井或者平面上连续性的预测。因此，针对以上问题，本文以东营凹陷沙四上纯下次亚段的致密滩坝砂储集层为例，利用薄片、物性及压汞等相关测试数据，计算致密滩坝砂孔喉分形维数，并探讨分形维数与储层物性、孔隙结构参数之间的关系，然后优选测井数据，建立了一种致密滩坝砂储层孔隙分形维数的测井预测模型，对东营凹陷西部地区致密滩坝砂储集层分形维数的平面分布进行了预测。基于分形维数储层评价标准，对典型区的致密滩坝砂储层进行了分类与评价，优选出有利储层分布区域，为确定开发区域动用次序提供依据。

1 研究区地质概况

研究区东营凹陷为渤海湾盆地中受陈南断裂控制的一个二级构造单元，是总体呈“北断南超”的中新生代箕状断陷盆地，南接广饶凸起，北部与陈家庄、林樊家和滨县凸起相接，东部为青坨子凸起，西部跨过青城凸起与惠民洼陷相接，面积约为5 850 km2（图1）[17-20]。古近系沙四段，主要分为沙四上亚段和沙四下亚段，其中，沙四上亚段又可进一步分为纯上和纯下两个次亚段。根据前人研究成果，致密滩坝砂储层主要分布于纯下次亚段，其中，坝砂单层厚度相对较大，一般大于2 m，滩砂单层厚度通常小于2 m，坝砂整体物性优于滩砂，分布区域如图1所示。

图1 东营凹陷沙四上纯下次亚段滩坝沉积体系分布图（据唐东[17]）Fig.1 Map of EsU4 beach-bar deposition system (after Tang[17])

1.1 岩石学特征

本次研究选取了12块东营凹陷沙四上纯下次亚段致密滩坝砂储集层样品进行岩石铸体薄片观察、物性和高压压汞测试，并基于压汞实验测试数据进行孔喉分形计算。根据薄片统计显示（表1），滩砂与坝砂的矿物碎屑组成近似，无明显差别。12 个样品岩性主要为岩屑质长石砂岩。石英含量为37%～51%，平均值为46.25%；长石含量为30%～38%，平均值为34%；岩屑类型主要为变质岩、沉积岩和岩浆岩为主，含量为16%～25%，平均值为19.67%。填隙物以泥质杂基为主，含量为2%～8%，平均值为5.5%。胶结物含量为4%～21%（平均值为9.42%），主要为方解石和白云石（平均含量分别为5.08%和3.42%），其次为硅质胶结（平均含量为0.92%）。砂岩中碎屑颗粒粒径为0.06～0.32 mm，平均值为0.15 mm，主要为粉砂和细砂，成分成熟度中等，分选中等，碎屑磨圆度主要为次棱角状。

表1 致密滩坝砂样品岩石学特征Table 1 Petrological properties of tight beach-bar sand samples

1.2 物性特征

滩坝砂储集层孔隙度和渗透率的统计结果表明（图2）：坝砂储集层孔隙度为1.0%～30.1%（主体分布在5%～18%，平均值为11.8%）；渗透率主体分布在0.1×10-3µm2～10×10-3µm2（平均值为0.92×10-3µm2）；滩砂储集层孔隙度为0.8%～26%（主体分布在3%～16%，平均值为8.9%），渗透率主体分布在0.01×10-3µm2～1×10-3µm2（平均值为0.67×10-3µm2），其孔隙度和渗透率均小于坝砂储集层。从整体上看，73%的样品孔隙度小于15%，渗透率小于10×10-3µm2，孔隙度平均值小于12%，渗透率平均值小于1×10-3µm2，为典型的致密砂岩储集层。

图2 滩坝砂储集层孔隙度和渗透率分布直方图Fig.2 Porosity and permeability histograms for beach-bar sand reservoir

1.3 压汞特征

12个样品的压汞曲线特征如图3所示，其中，样品编号1、2、4、9、10和11为坝砂样品，样品编号3、5、6、7、8和12为滩砂样品。压汞参数如表2所示，坝砂样品整体体现为粗歪度、分选较好的特点，排驱压力0.131～0.881 MPa（平均值为0.375 MPa），最大孔喉半径0.834～5.611µm（平均值为2.954µm），饱和度中值压力0.83～4.59 MPa（平均值为1.87 MPa），饱和度中值喉道半径0.16～0.886µm（平均值为0.541µm）；滩砂样品整体体现为细歪度、分选较差的特点，排驱压力0.334～3.562 MPa（平均值为1.753 MPa），最大孔喉半径0.206～2.201µm（平均值为0.752µm），饱和度中值压力1.95～21.09 MPa（平均值为11.82 MPa），饱和度中值喉道半径0.035～0.377 µm（平均值为0.117 µm）。滩砂样品孔喉半径整体小于坝砂样品。

表2 样品压汞参数表Table 2 Mercury intrusion parameters for 12 samples

图3 研究样品压汞曲线特征Fig.3 Mercury intrusion curves for 12 samples

根据压汞曲线可以获得孔喉半径频率分布图（图4），如图所示，孔喉半径多呈单峰分布，主要分布区间为0.1～1µm，但坝砂与滩砂的孔喉峰值半径（分布频率最高对应的孔喉半径）有一定的区别，其中，坝砂孔喉峰值半径为0.183～1.529 µm，平均值为0.8µm；滩砂孔喉峰值半径明显小于坝砂，为0.059～0.189 µm，平均值为0.119 µm。随孔喉峰值半径增加，孔隙度相应增大，当孔喉峰值半径小于0.4 µm时，孔隙度均小于10%。

图4 不同孔隙度滩坝砂样品的孔喉半径分布图Fig.4 Pore throat radius graphs for beach-bar sand samples with different porosity

2 孔隙结构分形特征与预测

2.1 孔喉结构分形特征

分形通常被定义为“一个粗糙或零碎的几何形状，可以分成数个部分，且每一部分都（或近似）是整体缩小后的形状”，即具有“自相似或标度不变性”[21]。前人的研究表明，致密砂岩的孔隙结构存在分形特征，分形理论可用于表征致密砂岩储集层的结构特征及其非均质性[22-24]；多孔岩石的分形维数通常介于2～3[25-27]，分形维数越小，储层孔喉分布越均匀，均质性越强[28-29]。

根据前人的研究成果[30]，储层孔隙分布表示式为：

根据油层物理学中对于饱和度的定义，孔径小于r的累积体积分数s即为润湿相饱和度[31]，且对于致密储层，孔喉半径的非均质性强，因此，rmin＜＜rmax，则有下式：

式中：Sw为润湿相饱和度，%；rmax为最大孔喉半径，µm。

对式6两边取对数可得：

式中：SHg为累积进汞饱和度，%。

根据12 块样品的lg（1-SHg）-lgr关系图（图5），可以计算出样品的分形维数。从图中可以看出，12 块致密滩坝砂样品的lg（1-SHg）-lgr关系图具有明显两段性，将lg（1-SHg）-lgr关系图中分段点对应的孔喉半径命名为孔喉转折半径，大于孔喉转折半径对应的孔隙为大孔喉，用ϕ1表示，小于孔喉转折半径对应的孔隙为小孔喉，用ϕ2表示。根据各孔隙空间的孔隙度加权平均得到整个孔喉的总分形维数D[27,32]，为

图5 12 块样品的lg（1-SHg）-lgr 关系图Fig.5 Diagrams of log(1-SHg) vs. logr for 12 samples

根据计算结果，两段分形维数D1、D2以及总分形维数D均在2～3范围内（表3），对坝砂和滩砂的分形维数分别统计，坝砂分形维数为2.222 4～2.531 9，平均值为2.441 2，而滩砂分形维数为2.585 0～2.742 4，平均值为2.650 9，明显高于坝砂。

表3 物性参数与分形维数Table 3 Physical properties and fractal dimensions

2.2 分形维数、储层物性和孔隙结构参数间的关系

2.2.1 分形维数与储层物性关系

分别建立分形维数与储层孔隙度和渗透率的关系。根据样品物性统计数据（表3），总分形维数D与储层孔隙度ϕ和渗透率k之间呈负向关关系（图6），R2分别为0.819和0.902，相关程度高。说明分形维数可以有效表征储层孔渗性能，即随着分形维数增大，储层的孔隙度与渗透率均变小，储层的孔渗性变差。

图6 分形维数D-孔隙度ϕ 和分形维数D-渗透率k 交会图Fig.6 Plots of (a) fractal dimension D vs. porosity ϕ and (b) D vs. permeability k

2.2.2 分形维数与孔隙结构参数关系

孔隙结构包含孔喉大小、分选性和连通性等信息，排驱压力PT、最大孔喉半径rmax、饱和度中值压力Pc50及孔喉中值半径r50表示孔喉大小，分选系数Sp表示孔喉分选性，退汞效率WE可表示孔喉连通性。根据分形维数与孔隙结构参数数据统计发现（图7），分形维数D与PT、rmax、Pc50、r50、Pc50、Sp和WE均具有较强相关性，其中，D与PT和Pc50呈正相关关系，与rmax和r50呈负相关关系，R2分别为0.566、0.519、0.823 和0.785，表明随着分形维数升高，孔喉半径有减小的趋势；D与Sp呈正相关关系，R2为0.732，表明随着分形维数升高，孔喉的分选性变差；D与WE呈负相关关系，R2为0.656，表明随着分形维数升高，孔喉的连通性变差。因此，分形维数可以作为综合反映储层孔隙结构特征的指标，分形维数越大，储层的储集条件越差。

图7 分形维数与孔隙结构参数交会图Fig.7 Intersection graphs of fractal dimension D and pore structure parameters

2.3 分形维数的测井预测方法

致密滩坝砂储层非均质性强，连续系统的钻井取心和实验分析费用高昂且难以操作，以有限的样品取其平均值代表整套砂体的做法不能反映细小差别，得到的结论很可能误导开发决策。测井信息可以间接地反映地层岩性、物性及其流体等特征，部分测井曲线与孔隙结构具有很强的相关性，因此，可尝试建立一种准确定量的、易于操作的孔隙分形维数测井预测模型，预测孔隙分形维数的单井纵向分布，解决实验资料纵向连续性差的问题，提高预测精度。

2.3.1 测井参数优选

通过对致密滩坝砂样品的分形维数与常用测井曲线进行相关性分析（表4），发现微电位与微梯度差值的绝对值（|RLML-RNML|）、声波时差（AC）、自然伽马（GR）和电阻率（RT）相关程度较高，相关系数分别为-0.907、-0.896、0.877和-0.867，随着分形维数D增大，|RLML-RNML|值减小、AC 值减小、GR 值增大、RT值减小。

表4 致密滩坝砂样品分形维数与测井类型相关系数Table 4 Correlation coefficient between fractal dimension and logging type for samples of dense beach-bar sand

（1）微电位与微梯度差值绝对值（|RLMLRNML|）和声波时差（AC）主要反映物性变化，D值增大，孔隙度和渗透率均呈减小的趋势（图6），导致|RLML-RNML|和AC值减小。

（2）自然伽马（GR）主要反映岩石的放射性。D为大值的样品多为滩砂样品，粒度较细，且泥质含量较D为小值的坝砂样品大，导致放射性增强，GR 值增大。

（3）电阻率（RT）主要反映含油性，随着D升高，孔喉半径减小，分选性变差（图7），油气运移至储层困难加大，导致含油量减小，RT值降低。

2.3.2 测井预测模型的建立与检验

建立东营凹陷致密滩坝砂测井预测模型的具体方法为：首先根据滩坝砂沉积相特点，剔除非滩坝砂层段；为了抵消各类测井参数数量级不同带来的影响，将优选出的测井曲线（|RLML-RNML|、AC、GR 和RT），进行标准归一化处理，标准归一化公式为：

式中：测井样本数据为x，xmax=max{x}，xmin=min{x}，测井样本归一化数据为X；然后，将归一化后的|RLMLRNML|、AC、GR和RT值对D进行回归拟合，得到D的测井预测公式：

式中：D为分形维数；|RLML-RNML|为微电极曲线差值绝对值，Ω·m；AC 为声波时差值，µs/m；GR为自然伽马值，API；RT为电阻率值，Ω·m。

将通过测井预测模型预测的分形维数与通过压汞数据计算得到的分形维数进行比对（图8），R2值达0.936，基本满足实际生产和科研的需要。

图8 分形维数D 测井预测值与压汞数据计算值相关图Fig.8 Correlation of fractal dimension D log predicted value and value calculated from mercury intrusion data

3 基于分形维数储层评价标准建立及有利储层预测

3.1 储层评价标准

分形维数D可以有效表征储层孔渗性能，综合反映储层孔隙结构特征，因此，推测分形维数D值大小与油井产能具有一定的相关性。为了揭示分形维数D对储层开发效果的影响，本次研究选取已动用区块无压裂措施下的日产油，累积产油与分形维数D之间的关系（图9），建立了基于分形维数D的东营凹陷致密滩坝砂储层评价标准（表5）。

图9 分形维数D 值与日产油、累积产油关系Fig.9 Relationship between fractal dimension D and (a) daily oil production and (b) cumulative oil production

根据日产油、累积产油与分形维数D的数据统计结果显示，以日产油量为主要分类指标，当2＜D＜2.35 时，可以认定为相对优质的储层，此时，对应层段的日产油量大于10 t/d，累积产油量大于90 t，孔隙度大于14%，渗透率大于0.9×10-3µm2，PT小于0.18 MPa，rmax大于4.1 µm，Pc50小于1.05 MPa，r50大于0.7µm，Sp小于0.19，WE大于33%；当2.35＜D＜2.55 时，认定为相对一般的储层，日产油量为5～10 t/d，累积产油量为50～90 t，孔隙度为10%～14%，渗透率为0.4×10-3µm2～0.9×10-3µm2，PT介于0.18～0.37 MPa，rmax为2.0～4.1µm，Pc50为1.05～2.3 MPa，r50为0.32～0.7µm，Sp为0.19～0.39，WE为26%～33%；当2.55＜D＜3 时，认定为相对较差的储层，日产油量小于5 t，累积产油小于50 t，孔隙度小于10%，渗透率小于0.4×10-3µm2，PT大于0.37 MPa，rmax小于2.0µm，Pc50大于2.3 MPa，r50小于0.32µm，Sp大于0.39，WE小于26%。

3.2 有利储层预测

本次研究以东营凹陷西部区块沙四上纯下次亚段的致密滩坝砂分布区域为例，利用分形维数D测井预测结果对致密滩坝砂体进行了储层评价。根据分形维数D平面分布等值线图（图10），结合东营凹陷致密滩坝砂储层综合评价标准（表5），优质储层的分形维数D为2～2.35，孔隙度大于14%，主要集中在两个区域，分别为B424井区和C107-C66-B3井区，全部分布于坝砂，该区域单井日产能高，日产油量10 t以上，为最有利的开发区域；一般储层的分形维数D为2.35～2.55，孔隙度为10%～14%，主要集中在四个区域，分别为B424 井区、F119 井区、C108-L108-C44井区和F145-F144-F134井区，主要为坝砂，少量为滩砂，该区域的单井日产油量为5～10 t，为较好的开发区域。利用分形维数D对致密滩坝砂储层进行综合评价，为下一步有利开发区域的选取及动用次序的确定提供了理论依据。

图10 东营凹陷西部区块分形维数D 值平面分布及储层分类Fig.10 Plane distribution of fractal dimension D and reservoir classification in western block, Dongying Sag

表5 东营凹陷致密滩坝砂储层综合评价标准Table 5 Comprehensive evaluation criterion for tight beach-bar sand reservoirs in Dongying Sag

表5 所列的储层评价标准适用于东营凹陷滩坝砂型致密储层，现场应用效果显著。本文所述基于分形维数的储层评价方法也在东营凹陷的砂砾岩和浊积岩型致密储层中进行了推广和应用，但是储层的优劣对应的分形维数范围与滩坝砂有一定的差异。因此，基于分形维数的储层评价方法可以应用于致密储层评价，但是标准的制定还需要根据致密储层的实际地质特点作相应的调整。

4 结论

（1）东营凹陷沙四上纯下次亚段分布典型的致密滩坝砂体沉积，利用压汞曲线数据计算孔隙分形维数，分析其与表征储层孔渗性能与孔隙结构参数的相关性，发现分形维数与总孔隙度ϕ、渗透率k、最大孔喉半径rmax、中值喉道半径r50、退汞效率WE呈负相关关系，与排驱压力PT、中值压力Pc50、分选系数Sp呈正相关关系，表明随着分形维数增大，储层孔渗性能、孔喉半径、分选性及连通性均变差。

（2）优选微电位与微梯度差值的绝对值（|RLML-RNML|）、声波时差（AC）、自然伽马（GR）和电阻率（RT）测井曲线参数，标准归一化拟合建立分形维数测井预测模型，将通过测井预测模型预测的分形维数与通过压汞数据计算得到的分形维数进行比对，R2值达0.936，基本满足实际生产和科研的需要。

（3）选取已动用区块无压裂措施下的日产油，累积产油与分形维数D之间的关系，建立基于分形维数D的东营凹陷致密滩坝砂储层评价标准，当2＜D＜2.35 时为优质储层，当2.35＜D＜2.55 时为一般储层，当2.55＜D＜3时为差储层。根据测井预测模型预测东营凹陷西部区块沙四上纯下次亚段的致密滩坝砂有利储层，认为优质储层主要集中在B424 井区和C107-C66-B3井区，为下一步储层动用次序的确定提供了理论依据。