2. 中国石油华北油田分公司山西煤层气勘探开发分公司, 山西晋城 048026;
3. 油气藏地质及开发工程国家重点实验室(成都理工大学), 四川成都 610059
2. Shanxi CBM Methane Exploration and Development Branch of PetroChina Huabei Oilfield Company, Jincheng, Shanxi, 048026, China;
3. State Key Laboratory of Oil and Gas Reservoir Geology and Exploration(Chengdu University of Technology), Chengdu, Sichuan, 610059, China
致密碎屑岩地层通常经历了强压实及复杂破坏性成岩作用,其原始孔隙度在地质演化过程中发生较大程度降低,现今孔隙度一般小于10%[1],覆压渗透率小于0.1 mD,但在一定条件下仍可形成具有勘探开发价值的工业性油气藏[2]。该类油气藏(如页岩油气及致密砂岩油气)具有烃类分布面积大、聚集程度差异大的特征[3, 4],裂缝的存在能够显著提升其储层质量[5],是其勘探、开发的“甜点区”。
裂缝除对地层岩石物性及含气性具有重要影响外,对岩石弹性也具有重要影响,其充填物还含有裂缝形成时间、期次及流体充注等方面的信息[6]。研究裂缝对岩石弹性的影响,对于研究岩石变形、破裂机制、断层区生长演化及地震各向异性等方面均具有一定的意义[7, 8]。前人多从岩心尺度开展力学及声学测试,研究裂缝对岩石动静态弹性的影响,未考虑裂纹的影响[9]。岩石中裂纹的开度通常小于50 μm,孔隙纵横比较小,因而仅能在显微镜下观察到[10]。裂纹是致密储层中最为重要的一种裂缝类型,但目前其对岩石弹性影响的研究还较少[11],主要集中在利用CT等测试技术刻画裂纹分布及富集程度[12],但这些测试通常费用较高,裂纹参数定量化较难,且这些测试技术均未探讨裂纹与岩心尺度裂缝间的关系。
DEM理论模型可以对岩石中裂纹的分布特征进行定量研究,J.G.Berryman等人[13]对DEM理论模型中岩石模量参数进行了求解;Kuster-Toksöz建立了干岩石骨架弹性模量与岩石孔隙形态间的定量关系[14];J.D.Eshelby[15]探讨了椭球状充填组分对岩石弹性性质的影响;Hongbing Li等人[16]利用DEM理论模型获得了干岩石体积模量及剪切模量的近似解;Xu-White模型考虑了岩石孔隙纵横比α对岩石纵横波波速的影响[17];Y.B.Budiansky[18]提出的自洽模型考虑了α对岩石剪切模量和体积模量的影响。致密碎屑岩储层的裂缝类型较为复杂,利用常规裂缝识别方法的效果通常较差。为此,笔者采用DEM理论模型对碎屑岩地层裂缝进行预测。
1 DEM理论模型碎屑岩地层中含气显示较好的层段,通常裂缝较为发育。根据裂缝识别结果,可以将单井目的层划分为裂缝发育段和欠发育段,以分析不同地层段的测井信息。利用DEM理论模型[19, 20]获得岩石孔隙纵横比α及干岩石体积模量,来分析裂缝对致密碎屑岩地层弹性的影响和不同尺度裂缝间的关联性。
DEM理论模型主要基于Berryman早期给出的有关岩石体积模量和剪切模量的耦合常微分方程[13],分别如下:


对于φ<10%的致密岩石,岩石孔缝可以分为球形(α>0.9)、针形(0.1≤α≤0.9)和硬币形(α<0.1)3类等效介质模型(见图1),通过二阶近似可以得到具有球形、针形及硬币形孔缝的干岩石体积模量表达方程[17]:




![]() |
图1 含不同形态孔缝的等效介质模型 Fig.1 Equivalent medium model with different pore-fracture forms |
碎屑岩地层的主要基质矿物为石英,此时Km=38.0 GPa,μm=44.4 GPa[21]。采用线性化方法对式(6)进行一阶求导,可以确定a和b的值[22]。
利用DEM理论模型可以确定不同类型α取值条件下的Kd值,α>0.1时可表示孔隙,α≤0.1时可表示裂纹[13]。根据Gassman方程[23],利用干岩石体积模量确定岩石体积模量:

根据声波测井资料可以计算岩石体积模量:

将由式(7)确定的体积模量与式(8)计算的K′进行对比[24],即可确定最小误差条件下的Kd及α。
2 实例分析 2.1 研究区概况及资料获取沁水盆地南部樊庄区块煤层气井区以开发二叠系山西组3号煤为主,沉积地层单元划分见表1。
地层系统 | 地层厚度/m | 主要岩性 | |
系组 | |||
第四系 | 0~18 | 黄土、亚黏土、卵砾石、亚黏土夹钙质结核层 | |
三叠系 | 0~72 | 石英砂岩、砾岩、砂质泥岩、暗紫色泥岩 | |
石千峰组 | > 0~281 | 砂质泥岩、泥岩夹细砂岩 | |
上石盒子组 | 342~521 | 长石石英砂岩及砂质泥岩、泥岩 | |
三叠系 | 下石盒子组 | 68~101 | 长石质砂岩、粉砂岩及紫红色泥岩、铝质泥岩 |
山西组 | 39~79 | 细粒砂岩、粉砂岩、泥岩互层,含煤4—5层 | |
石炭系 | 太原组 | 80~110 | K2—K6灰岩、中—粗粒砂岩、粉砂岩、泥岩,含煤7—9层 |
本溪组 | 0~14 | 铝土岩、褐铁矿、石英砂岩、泥岩局部夹薄层灰岩和煤线 | |
奥陶系 | 峰峰组 | 90~150 | 上部白云质泥灰岩、泥灰岩、泥质灰岩夹石膏,下部厚层状灰岩 |
老井复查及气测录井显示,太原组、山西组及下石盒子组碎屑岩地层气测异常极为普遍,具有致密气勘探、开发潜力。虽然目前研究区晚古生代沉积碎屑岩地层埋深较浅,但该地区石炭系、二叠系碎屑岩地层在晚三叠世之前一直持续沉降,最大古埋深一般大于3 000 m,强烈压实是造成地层致密的重要原因[20]。
X1井位于沁南樊庄区块,为直井,山西组地层埋深491.00~514.00 m,具有高密度及低孔隙度的致密特征。该井具有齐全的常规、阵列及成像测井资料。
通过成像测井及岩心观察可以识别出X1井5个裂缝发育段,其中未充填缝占52%,半充填缝占46%,全充填高阻缝占2%,该井裂缝类型以有效缝(前2类)为主。提取X1井所研究井段的测井数据,总样本数为145个;其中,裂缝发育段样本数据75个,裂缝不发育段样本数据70个。该地区山西组碎屑岩地层中较好的含气性与裂缝密切相关,当致密储层中裂缝不发育时,烃类的充注压力较高,很难实现烃类的规模有效富集及成藏[20]。研究区致密碎屑岩储层应定义为裂缝性致密储层,储层孔隙度低于8.0%,渗透率低于0.3 mD。
2.2 裂缝对致密碎屑岩储层弹性的影响 2.2.1 裂纹形态对岩石弹性的影响利用DEM模型对X1井测井数据进行处理分析,根据地层微电阻率扫描成像(FMI)测井资料及岩心裂缝识别结果将地层划分为裂缝发育段和不发育段。一般来说,低伽马值砂岩段地层的裂缝发育程度要比高伽马值泥岩地层高一些;同时,裂缝的存在,可引起井径、电阻率、补偿中子及声波时差等常规测井参数的响应[25]。从常规测井参数交汇图可以看出(见图2),由于研究区海陆过渡相碎屑岩沉积地层岩性变化频繁,利用常规测井参数不能较好地识别出裂缝发育段地层。
处理后的X1井目的层Ks、μs与孔隙度之间的关系如图3所示。从图3可以看出,裂缝发育段岩石的Ks和μs要略小于裂缝不发育段地层,随着孔隙度增大,地层岩石的Ks和μs均呈逐渐降低的趋势,说明利用弹性力学参数能较好地识别出裂缝发育段地层。
以α=0.10作为岩石中主要发育孔隙或裂纹的界线:当α>0.1时,岩石中主要发育孔隙,如图4(a)所示,显微薄片显示该石英砂岩孔隙空间以粒内溶孔为主;当α≤0.1时,岩石中主要发育裂纹,如图4(b)所示,该岩屑石英砂岩中裂纹发育,主要包括粒缘裂纹和穿粒裂纹。这些裂纹的分布通常具有一定方向性,且主裂纹通常具有较大的开度和较长的延伸距离,表明岩石破裂前集聚的应变能较大;主裂纹在延伸末端往往出现分支,图4(b)中的主裂纹末端出现2条分支,一条穿过石英颗粒并发生尖灭,另一条则沿着石英颗粒边缘延伸并逐渐发生尖灭。
![]() |
图2 X1井目的层岩石常规测井参数电阻率与伽马的关系 Fig.2 Relationship of conventional logging parameters vs. GR of the target zones in Well X1 |
![]() |
图3 考虑裂缝发育情况的地层岩石 Fig.3 Relationship of porosity vs. rock Ks and μs based on fracture development situations |
![]() |
图4 山西组致密碎屑岩孔缝特征 Fig.4 Characteristic of pores and fractures in tight clastic rocks of Shanxi Formation |
图3中3条虚线分别代表利用DEM理论模型确定的含球形(α=1.0)、针形(α=0.1)及硬币形(α=0.01)孔缝介质的Ks、μs与孔隙度的理论取值。从图3可以看出,对于裂缝不发育段碎屑岩地层,主要发育有孔隙;而对于裂缝发育段碎屑岩地层,主要发育有裂纹;裂缝发育段的α绝大部分小于0.01。因此,利用DEM理论模型可以较好地识别出裂纹发育地层段(对应α≤0.1),这些层段也恰好对应着利用岩心尺度识别出的裂缝发育段(图3中的橙红色数据点)。这表明,岩心尺度裂缝与岩石内部裂纹间具有较好的关联性。
从图3还可以看出,随着孔隙纵横比α由1.0向0.01逐渐降低的过程中,地层岩石的Ks和μs也逐渐降低。这表明岩石内部在由孔隙逐渐向裂纹占主导转变的过程中,岩石体积模量和剪切模量会降低,裂纹相比孔隙而言更容易造成岩石的弹性性质发生改变。这主要是因为岩石中的裂纹在应力作用下易于扩展和张开[26]。DEM理论模型主要从力学机制出发研究岩石的弹性性质,根据前述分析结果,利用岩石弹性参数表征碎屑岩地层裂缝发育程度也同样有效。
2.2.2 岩石是否饱和地层流体对裂缝影响储层弹性的影响利用DEM理论模型可以确定地层岩石的Kd,如图5所示。
![]() |
图5 地层岩石Ks和Kd与孔隙度的关系 Fig.5 Relationship of porosity vs. rock Ks and Kd |
从图5可以看出,与Ks类似,随着孔隙度增大,Kd逐渐降低,岩石Kd<Ks。因此,岩石在由饱和地层流体到转变为干岩石的过程中,体积模量会逐渐降低。
为了定量表征这种降低程度,绘制地层岩石体积模量增加幅度(Ks-Kd)/Ks与孔隙度的关系图,见图6。
![]() |
图6 地层岩石体积模量增加幅度与孔隙度的关系 Fig.6 Relationship of porosity vs. the decreasing amplitude of rock bulk modulus |
从图6可以看出,对于裂缝不发育段地层,随着孔隙度从4.15%增加到4.85%,体积模量的增加幅度变化不大,为0.8%~8.5%,平均为3.1%。同时,该致密碎屑岩储层裂缝不发育。一般来说,裂缝能增加岩石中孔隙团形成的概率[27]。当裂缝不发育时,岩石内部孔隙间连通性相对较差,因而饱和地层流体与否对岩石体积模量的影响程度较小。
从图6还可以看出,对于裂缝发育段地层,随地层岩石孔隙度增大,体积模量增加幅度逐渐增大,但呈现出明显的分段(两段)性。孔隙度较低时,体积模量的增加幅度较小,为1.5%~10.0%,平均为7.0%;孔隙度较高时,体积模量的增加幅度显著提高,主要分布在10.0%~40.0%,平均为23.0%。
将两段数据点分别拟合(见图6),可得拐点处所对应的地层岩石孔隙度为4.7%。出现拐点的原因主要是由于裂缝具有较强的沟通能力,随着地层孔隙度增大,其沟通能力越来越强。岩石中裂纹的形成经历了初始稳定扩展及后期非稳定指数型扩展的过程[28, 29, 30, 31]。图6中拐点的出现表明,当孔隙度超过某一值时,岩石内部裂纹极为发育,因此,该孔隙度可作为岩石内部裂纹发育程度的评价指标。
同时,研究结果表明,岩石内部裂纹的发育程度与物性具有一定正相关性,大量裂纹的存在沟通了致密储层内部原本孤立的孔隙,从而使岩石有效孔隙度增大。
2.2.3 不同类型地层岩石Kd/μd与孔隙度的关系由DEM理论模型提取Kd,由于岩石中饱和流体与否对其剪切模量无影响,因而μd=μs[32, 33]。X1井目的层不同类型碎屑岩地层岩石Kd/μd与孔隙度之间的关系如图7所示。
![]() |
图7 地层岩石Kd/μd与孔隙度的关系 Fig.7 Relationship of rock Kd/μd vs. porosity |
对于裂缝不发育段地层,其α绝大多数大于0.1,大部分大于0.5,只有少量数据点小于0.1,岩石Kd/μd与孔隙度在整体上具有较好的正相关性;对于裂缝发育段地层,其α均小于0.1,大部分数据值点小于0.01,其Kd/μd与孔隙度整体上具有较好的负相关性(见图7)。文献[34]基于DEM理论给出了考虑不同α条件下碎屑岩Kd/μd与孔隙度的变化关系,如图8所示。本文裂缝不发育地层数据大致相当于图8中0.5<α≤0.9的情况,而裂缝发育地层数据大致相当于图8中α=0.01的情况。这2种情况下的岩石Kd/μd与孔隙度变化趋势与本文结果均较为一致,表明本文基于DEM理论所获得的各弹性参数结果可靠,方法可行。
![]() |
图8 不同α条件下碎屑岩Kd/μd与孔隙度的关系 Fig.8 Relationship of porosity vs. Kd/μd of clastic rocks under different conditions of α |
1) 利用DEM理论模型对致密岩石裂缝参数进行反演,结果表明所确定的裂纹发育段与岩心尺度裂缝发育段具有一定关联性。随着岩石孔隙纵横比α由1.0向0.01转变过程中,地层岩石体积模量和剪切模量也逐渐发生降低,表明裂纹相比孔隙而言更容易引起岩石弹性发生改变。
2) 岩石是否饱和地层流体对裂缝影响储层弹性具有一定影响。在岩石由干岩石饱和地层流体过程中,对于裂缝不发育地层段,其体积模量增加幅度不大,平均为3.12%;对于裂缝发育地层段,其体积模量随孔隙度增大而增大。孔隙度较低时增大幅度较小,平均为7.0%;孔隙度较高时,增大幅度显著提高,平均为23.0%。拐点处所对应岩石孔隙度可作为岩石内部裂纹发育程度的评价指标,裂纹的发育程度与物性具有一定正相关性。
3) 裂缝不发育地层段的岩石Kd/μd与孔隙度之间具有较好的正相关性;裂缝发育地层段的岩石Kd/μd与孔隙度之间具有较好的负相关性。
4) 研究表明,通过DEM理论模型所获得的地层岩石孔隙纵横比虽然为有效值,但可以对地层岩石孔缝特征进行定量表征;裂缝的形成源于力学机制,利用岩石弹性性质定量表征碎屑岩地层裂缝发育特征具有有效性。
[1] | FARRELL N J C,HEALY D,TAYLOR C W .Anisotropy of permeability in faulted porous sandstones[J].Journal of Structural Geology,2014,63:50-67. |
[2] | DANA E,SKOCZYLAS F.Gas relative permeability and pore structure of sandstones[J].International Journal of Rock Mechanics and Mining Sciences,1999,36(5):613-625. |
[3] | NELSON P H.Pore-throat sizes in sandstones,tight sandstones,and shales[J].AAPG Bulletin,2009,93(3):329-340. |
[4] |
张映红,路保平,陈作,等.中国陆相致密油开采技术发展策略思考[J].石油钻探技术,2015,43(1):1-6. ZHANG Yinghong,LU Baoping,CHEN Zuo,et al.Technical strategy thinking fordeveloping continental tight oil in China[J].Petroleum Drilling Techniques,2015,43(1):1-6. |
[5] | LASH G G,ENGELDER T.An analysis of horizontal microcracking during catagenesis:example from the catskill delta complex[J].AAPG Bulletin,2005,89(11):1433-1449. |
[6] |
张晋言,刘海河,刘伟,等.基于双“泥质指示因子”的砂泥岩薄互层饱和度模型[J].石油钻探技术,2015,43(2):59-62. ZHANG Jinyan,LIU Haihe,LIU Wei,et al.Thin sand and mudstone interbed saturation model based on a dual mudstone indicator[J].Petroleum Drilling Techniques,2015,43(2):59-62. |
[7] | ANDERS M H,LAUBACH S E,SCHOLZ C H.Microfractures:a review[J].Journal of Structural Geology,2014,69(B):377-394. |
[8] |
尹帅,单钰铭,周文,等.Hoek-Brown准则在致密砂岩弹性参数测井解释中的应用[J].石油钻探技术,2015,43(1):88-95. YIN Shuai,SHAN Yuming,ZHOU Wen,et al.Application of Hoek-Brown criterion for tight sandstone elastic parameters in log interpretation[J].Petroleum Drilling Techniques,2015,43(1):88-95. |
[9] | GALE J F W,LAUBACH S E,OLSON J E,et al.Natural fractures in shale:a review and new observations[J].AAPG Bulletin,2014,98(11):2165-2216. |
[10] | ZHANG Jiaxiang,WONG Tengfong,DAVIS D M.Micromechanics of pressure-induced grain crushing in porous rocks[J].Journal of Geophysical Research,1990,95(B1):341-352. |
[11] | LAMA R D,VUTUKURI V S.Handbook on mechanical properties of rocks:vol.Ⅱ[M].Clausthal:Transport Technology Publications,1978. |
[12] | MENNDEZ Beatriz,ZHU Wenlu,WONG Tengfong.Micromechanics of brittle faulting and cataclastic flow in Berea sandstone[J].Journal of Structural Geology,1996,18(1):1-16. |
[13] | BERRYMAN J G,PRIDE S R,WANG H F.A differential scheme for elastic properties of rocks with dry or saturated cracks[J].Geophysical Journal International,2002,151(2):597-611. |
[14] | BERRYMAN J G.Single-scattering approximations for coefficients in Biot’s equations of poroelasticity[J].Journal Acoust Society American,1992,91(2):551-570. |
[15] | ESHELBY J D.The determination of the elastic field of an ellipsoidal inclusion,and related problems[J].Proceedings of the Royal Society A,Mathematical,Physical and Engineering Sciences,1957,241:376-394. |
[16] | LI Hongbing,ZHANG Jiajia.Analytical approximations of bulk and shear moduli for dry rock based on the differential effective medium theory[J].Geophysical Prospecting,2012,60(2):281-292. |
[17] | KEYSR G,XU Shiyu.An approximation for the Xu-White velocity model[J].Geophysics,1999,67(5):1406-1414. |
[18] | BUDIANSK Y B.On the elastic moduli of some heterogeneous materials[J].Journal of the Mechanics and Physics of Solids,1965,13(4):223-227. |
[19] | BREEYMAN J G,Berge P A.Critique of two explicit schemes for estimating elastic properties of multiphase composites[J].Mechanics of Materials,1996,22(2):149-164. |
[20] | SU Xianbo,LIN Xiaoying,ZHAO Mengjun,et al.The upper Paleozoic coalbed methane system in the Qinshui Basin, China[J].AAPG Bulletin,2005,89 (1):81-100. |
[21] |
尹帅,丁文龙,孙雅雄,等.泥页岩单轴抗压破裂特征及UCS影响因素[J].地学前缘,2016,23(2):75-95. YIN Shuai,DING Wenlong,SUN Yaxiong,et al.Shale uniaxial compressive failure property and the affecting factors of UCS[J].Earth Science Frontiers,2016,23(2):75-95. |
[22] | LI Hongbing,ZHANG Jiajia.Elastic moduli of dry rocks containing spheroidal pores based on differential effective medium theory[J].Journal of Applied Geophysics,2011,75(4):671-678. |
[23] | MAVKO G,MUJERKI T,DVORKIN J.Rock physicshandbook:tools for seismic analysis in porous media[M].Cambridge:Cambridge University Press,2003. |
[24] | YASAR E.Failure and failure theories for anisotropic rocks//Proceedings of the 17th International Mining Congress and Exhibition of Turkey-IMCET 2001[M].Ankara:Chamber of Mining Engineers of Turkey,2001:417-424. |
[25] |
丁文龙,尹帅,王兴华,等.致密砂岩气储层裂缝评价方法与表征[J].地学前缘,2015,22(4):173-187. DING Wenlong,YIN Shuai,WANG Xinghua,et al.Assessment method and characterization of tight sandstone gas reservoir fractures[J].Earth Science Frontiers,2015,22(4):173-187. |
[26] | BIENIAWSKI Z T.Mechanism of brittle fracture of rock:partⅡ:experimental studies[J].International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts,1967,4(4):407-423. |
[27] |
李伟,要惠芳,刘鸿福,等.基于显微CT的不同煤体结构煤三维孔隙精细表征[J].煤炭学报,2014,39(6):1127-1131. LI Wei,YAO Huifang,LIU Hongfu,et al.Advanced characterization of three-dimensional pores in coals with different coal-body structure by micro-CT[J].Journal of China Coal Society,2014,39(6):1127-1131. |
[28] | LAJTAI E Z.Microscopic fracture processes in a granite[J].Rock Mechanics and Rock Engineering,1998,31(4):237-250. |
[29] | KONCAGUL E C,SANTI P M.Predicting the unconfined compressive strength of the Breathitt shale using slake durability,shore hardness and rock structural properties[J].International Journal of Rock Mechanics and Mining Science, 1999,36(2):139-153. |
[30] | STEIGER R P,LEUNG P K.Quantitative determination of the mechanical properties of shale[R].SPE 18024,1992. |
[31] |
尹帅,丁文龙,赵威,等.基于阵列声波测井的海陆过渡相碎屑岩地层裂缝识别方法[J].石油钻探技术,2015,43(5):75-82. YIN Shuai,DING Wenlong,ZHAO Wei,et al.Fracture identification method for marine-continental transitional clastic rocks based on the array acoustic logging[J].Petroleum Drilling Techniques,2015,43(5):75-82. |
[32] | BATZLE M,WANG Zhijing.Seismic properties of pore fluids[J].Geophysics,1992,57(11):1396-1408. |
[33] |
尹帅,丁文龙,黄昌杰,等.非饱和致密砂岩储层含气饱和度测井评价[J].天然气地球科学,2016,27(1):1-10. YIN Shuai, DING Wenlong, HUANG Changjie, et al.Gas saturation logging evaluation of partially saturated tight sandstone reservoir[J].Natural Gas Geoscience,2016,27(1):1-10. |
[34] | LI Hongbing,ZHANG Jiajia.Modulus ratio of dry rock based on differential effective medium theory[J].Geophysics,2010, 75(2):N43-N50. |