技术领域
[0001] 本
发明属于油气田开发过程中的页岩气藏动态储量计算领域，具体为一种考虑多种因素的页岩气藏储量求解方法，在储量计算过程中提供有效指导。
背景技术
[0002] 页岩气主要以自由气、
吸附气和溶解气3种形式储存在高
碳泥页岩或暗色泥岩中。自由气主要分布在基质孔隙和裂缝孔隙中；吸附气主要分布在基质系统微孔隙的内表面；
溶解气相对较少，且主要分布在
沥青质和少量页岩油中。通常页岩气储量计算方法有3种：
类比法、容积法和物质平衡法。在页岩气勘探开发过程中，储量评价十分重要，物质平衡法是目前评价页岩气藏储量的常用方法。
[0003] 类比法主要用于勘探初期粗略评价页岩气藏地质储量，张金华、李宏勋等分别提出了类比法的适用条件。容积法适用于计算页岩气藏静态地质储量，徐海霞、陈元千等提出了页岩气藏地质储量的计算方法，但未考虑吸附相占据的孔隙度。Ambrose等建立了考虑吸附相孔隙度页岩气藏储量计算模型，但该模型未考虑多组分吸附及溶解气。曹延宽、曲占庆等虽研究了考虑多组分吸附的页岩气藏储量，但未考虑溶解气的储量。姜瑞忠等人从分子
角度推导了页岩气储量计算公式，但也只考虑了单一甲烷分子构成的吸附相储量。物质平衡法常用来计算页岩气藏动用地质储量，计算时需要大量生产数据，众多学者在采用物质平衡方程计算页岩气藏储量时，未将溶解在
干酪根中的溶解气考虑在内，导致计算结果不准确。要准确地计算页岩气的储量，应该减少假设条件，建立一个无限接近真实气藏的完
整理论模型。本发明建立同时考虑吸附相视孔隙度干酪根中溶解气储量、吸附相
密度、基质和裂缝孔隙体积随压
力变化的物质平衡方程。
发明内容
[0004] 本发明主要是克服
现有技术中的不足之处，提出一种计算更准确的考虑多种因素的页岩气藏储量求解方法。
[0005] 本发明解决上述技术问题所提供的技术方案是：一种考虑多种因素的页岩气藏储量求解方法，包括以下步骤：
[0006] 步骤S100、测试、收集页岩气藏的基本参数和气藏历史生产数据；
[0007] 步骤S200、考虑基质和裂缝自由气、吸附相视孔隙度、吸附相密度、吸附气临界
解吸压力及吸附气解吸对基质固相产生的
变形情况，修正
岩石压缩系数，并同时考虑了吸附及溶解气，建立页岩气藏物质平衡方程：
[0008]
[0009] 式中：Vp为累计产气量，108m3；Vm为页岩基质中自由气体积，108m3；Vf为裂缝中自由气储量，108m3；Bg为当前地
层压力下气体的体积系数，m3/m3；Bgi页岩气的体积系数，m3/m3；Cx为修正的岩石压缩系数，MPa-1；Cw为
地层水压缩系数，MPa-1；Cf为裂缝压缩系数，MPa-1；Smwi为基质中束缚水
饱和度，f；Sfwi为裂缝中束缚水饱和度，f；Pi为气藏原始压力，MPa；PL为兰氏压力，MPa；Pd为临界解吸压力，MPa；φt为基质孔隙度，f；φa为吸附相视孔隙度，f；T为气层
热力学温度，K；b为吸附量等于饱和量一半时的气体压力(兰氏压力常量)，MPa-1；p为储层压力，MPa；R为气体常数，MPa·m3/(K·kmol)；V0为气体摩尔体积,10-3m3/mol；E为
杨氏模量,MPa；ρr为干酪根的相对密度，g/m3；φads为吸附相的孔隙度，f；φorg为有机质中自由气所占据的孔隙度，f；CTOC为总有机碳含量，％；
[0010] 步骤S300、根据基本参数、气藏生产数据和页岩气藏物质平衡方程得到Y/M和F/M的一次线性关系，其关系曲线的斜率为页岩基质中自由气储量Vm，截距为裂缝岩石中自由气储量Vf；
[0011] 其中Y＝VpBg，M＝Bg-Bgi+CcfBgiΔp，
[0012]
[0013] Δp＝pi-p，
[0014] 式中：Ccm为基质孔隙系统有效压缩系数，MPa-1；Ccf为裂缝系统有效压缩系数，MPa-1；
[0015] 步骤S400、最后得到页岩气藏总储量：
[0016]
[0017] 式中：c(pi)为原始压力下，干酪根中甲烷的
溶解度，m3/m3；
[0018] 进一步的技术方案是，所述气藏生产数据包括多组数据，每组数据包括地层压力p、地层压力为p时的气体体积系数Bg、累计产气量Vp。
[0019] 进一步的技术方案是，所述步骤S200页岩气藏物质平衡方程的具体建立过程为：
[0020] 步骤S201、计算当地层压力下降到p时的岩石基质及束缚水膨胀体积ΔVm：
[0021]
[0022] 式中：Smwi为基质中束缚水饱和度，f；ΔVm基质孔隙体积减小量,108m3；
[0023] 步骤S202、计算当地层压力下降到p时的裂缝岩石及束缚水膨胀体积ΔVf：
[0024]
[0025] 式中：Sfwi裂缝中束缚水饱和度，f；ΔVf裂缝孔隙体积减小量,108m3；
[0026] 步骤S203、当地层压力下降到p时，在地层条件下，溶解气的扩散量可由下式计算：
[0027]
[0028] 步骤S204、根据页岩气藏体积守恒定律可得到：累计产气量＝基质内自由气膨胀体积+基质颗粒变形和束缚水发生弹性膨胀体积+吸附气解析出的气体体积+扩散的溶解气体积+吸附气解吸导致基质收缩增加的孔隙体积+裂缝内自由气的膨胀体积+裂缝内孔隙体积减小和束缚水膨胀体积，则页岩气藏物质平衡方程为：
[0029]
[0030] 本发明的有益效果是：本发明综合考虑了页岩气吸附及溶解气、基质和裂缝自由气、吸附相视孔隙度、吸附气临界解吸压力及吸附气解吸对基质产生的收缩效应，修正了岩石压缩系数，建立了一种新的页岩气藏物质平衡方程，该方程对于合理计算页岩气藏动态储量具有重要指导意义。
附图说明
[0031] 图1为新、旧岩石物理模型对比图；
[0032] 图2为有机质中干酪根的空间分布图；
[0033] 图3为物质平衡方程储量回归图；
具体实施方式
[0035] 下面结合实例和附图对本发明做更进一步的说明。
[0036] 本发明是一种考虑多种因素的页岩气藏储量求解方法，包括以下步骤：
[0037] 步骤S100、测试、收集页岩气藏的基本参数和气藏历史生产数据；所述气藏生产数据包括多组数据，每组数据包括地层压力p、地层压力为p时的气体体积系数Bg、累计产气量Vp。
[0038] 步骤S200、考虑基质和裂缝自由气、吸附相视孔隙度、吸附相密度、吸附气临界解吸压力及吸附气解吸对基质固相产生的变形情况，修正岩石压缩系数，并同时考虑了吸附及溶解气，建立页岩气藏物质平衡方程：
[0039]
[0040] 式中：Vp为累计产气量，108m3；Vm为页岩基质中自由气体积，108m3；Vf为裂缝中自由8 3 3 3 3 3
气储量，10m ；Bg为当前地层压力下气体的体积系数，m/m；Bgi页岩气的体积系数，m/m ；Cx为修正的岩石压缩系数，MPa-1；Cw为地层水压缩系数，MPa-1；Cf为裂缝压缩系数，MPa-1；Smwi为基质中束缚水饱和度，f；Sfwi为裂缝中束缚水饱和度，f；Pi为气藏原始压力，MPa；PL为兰氏压力，MPa；Pd为临界解吸压力，MPa；φt为基质孔隙度，f；φa为吸附相视孔隙度，f；T为气层热力学温度，K；b为吸附量等于饱和量一半时的气体压力(兰氏压力常量)，MPa-1；p为储层压力，MPa；R为气体常数，MPa·m3/(K·kmol)；V0为气体摩尔体积,10-3m3/mol；E为杨氏模量,MPa；ρr为干酪根的相对密度，g/m3；φads为吸附相的孔隙度，f；φorg为有机质中自由气所占据的孔隙度，f；CTOC为总有机碳含量，％；
[0041] 步骤S300、根据基本参数、气藏历史生产数据和页岩气藏物质平衡方程得到Y/M和F/M的一次线性关系，其关系曲线的斜率为页岩基质中自由气储量Vm，截距为裂缝岩石中自由气储量Vf；
[0042] 其中Y＝VpBg，M＝Bg-Bgi+CcfBgiΔp，
[0043]
[0044]Δp＝pi-p，
[0045] 式中：Ccm为基质孔隙系统有效压缩系数，MPa-1；Ccf为裂缝系统有效压缩系数，-1MPa ；
[0046] (1)式通过上述化简，可得：
[0047]
[0048] 根据基本参数、气藏生产数据和页岩气藏物质平衡方程得到Y/M和F/M的一次线性关系，其关系曲线的斜率为页岩基质中自由气储量Vm，截距为裂缝岩石中自由气储量Vf；
[0049] 步骤S400、最后得到页岩气藏总储量：
[0050] 式中：c(pi)为原始压力下，干酪根中甲烷的溶解度，m3/m3；CTOC为总有机碳含量，％；φads为吸附相的孔隙度，f；φorg为有机质中自由气所占据的孔隙度，f；b为兰氏压力常量，MPa-1；
[0051] 由于页岩气与
煤层气吸附机理相近，因此可借鉴研究
煤层中吸附气的方法对页岩气藏中吸附气的机理进行分析。对于微孔类吸附，常可用Langmuir等温模型处理，吸附形式：
[0052] 式中：VL为兰氏体积常量，m3/t；b为吸附量等于饱和量一半时的气体压力(兰氏压-1力常量)，MPa ；p为储层压力。
[0053] 微孔内壁，由于存在分子吸附层，且因这种吸附层在一般情况下可视为不流动吸附层，因此在气井打开储层的早期，自由气孔隙度比实际孔隙度小。定义被吸附气占据的吸附相视孔隙度计算式为：
[0054]
[0055] 式中：φa为吸附相视孔隙度，f；M为表观
天然气相对分子
质量，g/mol；ρs吸附相密度g/cm3；a为单位换算常数值，取1.318×10-6；ρb为页岩岩石密度，g/cm3；
[0056] 页岩气藏压力下降到p时的吸附气解吸量为：
[0057] 式中：pd为临界解吸压力，MPa；φt为基质孔隙度，f；
[0058] 根据Bangham固体变形理论和Langmuir等温吸附模型，得到吸附气解析导致岩石形变量与地层压力之间的关系：
[0059] 式中：R为气体常数,MPa·m3/(K·kmol)；T为气层热力学温度,K；V0为气体摩尔体积10-3m3/mol；Δε为吸附气解吸引起的收缩基质变形程度；E为杨氏模量，MPa；pd为临界解吸压力，MPa；
[0060] 随着地层压力降低，吸附气发生解吸，基质开始收缩；同时
流体产出，有效
应力增加，岩石的骨架体积被收缩，因此页岩基质总变形量为基质收缩形变量和岩石弹性膨胀形变量之和。即
[0061] 式中：Δεz为有效应力下岩石弹性膨胀和基质收缩形变量之和；Cm为岩石弹性压缩系数，MPa-1；pi为原始地层压力，MPa；
[0062] 由Seidle提出的孔隙度与基质形变量的关系，可得基质孔隙度与压力的表达式为:
[0063]
[0064] 式中：φ为变形后基质孔隙度，f；φt为基质孔隙度，f；
[0065] 由压缩系数的定义可得:
[0066]
[0067] 联立式(9)和式(10)可得
[0068] 随着地层压力的下降，岩石颗粒变形和束缚水发生弹性膨胀，则岩石基质中孔隙体积的变化值为:
[0069]
[0070] 地层压力下降到p时，裂缝岩石及束缚水膨胀体积为:
[0071]
[0072] 式中：Cx修正后岩石的压缩系数，MPa-1；Cw地层水压缩系数，MPa-1；Cf裂缝压缩系数，MPa-1；ΔVm基质中孔隙体积的变化量，108m3；ΔVf裂缝孔隙体积的变化量，108m3；
[0073] 吸附气解吸导致的收缩基质变形程度和储层压力的关系：
[0074] 式中Δεd吸附气解吸导致的收缩基质变形程度；
[0075] 吸附气解吸导致的基质孔隙体积的变化量为：
[0076]
[0077] 沥青质中甲烷溶解度的经验公式可由下式来计算，为了避免在压力较低时溶解度出现负值，在温度为25-100℃、压力为0时，溶解度也为0。
[0078]
[0079] 式中：b1＝-0.018931；b2＝-0.850480；b3＝827.26；b4＝-635.26
[0080] 甲烷在干酪根中的溶解度和在沥青质中一样，则溶解气体积的计算公式为：
[0081] Vdsk(p)＝c(p)Vsk (17)
[0082] 式中：Vsk＝VrockVdiff；Vsk为干酪根总体积，m3；
[0083] 又因为总干酪根占基质的体积分数Vtker可以表示成：
[0084]
[0085] 式中：Vdiff为固体干酪根占基质总体积的体积分数，f；
[0086] 干酪根的总体积为：
[0087]
[0088] 式中
[0089] 溶解气的储量为：
[0090] 式中：φads为吸附相的孔隙度，f；φorg为有机质中自由气所占据的孔隙度，f；
[0091] 当地层压力下降到p时，在地层条件下，溶解气的扩散量可由下式计算：
[0092]
[0093] 根据体积守恒定律可得到：累计产气量＝基质内自由气膨胀体积+基质颗粒变形和束缚水发生弹性膨胀体积+吸附气解析出的气体体积+扩散的溶解气体积+吸附气解吸导致基质收缩增加的孔隙体积+裂缝内自由气的膨胀体积+裂缝内孔隙体积减小和束缚水膨胀体积
[0094]
[0095] 上式为一个非线性的方程，为使其简单实用，必须进行线性化处理。
[0096] 令：Δp＝pi-p，
[0097] 化简得
[0098]
[0099] 令Y＝VpBg，M＝Bg-Bgi+CcfBgiΔp，
[0100]
[0101] (23)式通过化简，可得：
[0102]
[0103] 根据基本参数、气藏历史生产数据和页岩气藏物质平衡方程得到Y/M和F/M的线性关系，其关系曲线的斜率为页岩基质中自由气储量Vm，截距为裂缝岩石中自由气储量Vf；
[0104] S400、最后得到页岩气藏总储量：
[0105] 上述式中：Vp为累计产气量，108m3；Vm为页岩基质中自由气体积，108m3；Vf为裂缝中自由气储量，108m3；Bg为当前地层压力下气体的体积系数，m3/m3；Bgi页岩气的体积系数，m3/m3；Cx为修正的岩石压缩系数，MPa-1；Cw为地层水压缩系数，MPa-1；Cf为裂缝压缩系数，MPa-1；Smwi为基质中束缚水饱和度，f；Sfwi为裂缝中束缚水饱和度，f；Pi为气藏原始压力，MPa；PL为兰氏压力，MPa；Pd为临界解吸压力，MPa；φt为基质孔隙度，f；φa为吸附相视孔隙度，f；T为气层热力学温度，K；b为吸附量等于饱和量一半时的气体压力(兰氏压力常量)，MPa-1；p为储层压力，MPa；R为气体常数，MPa·m3/(K·kmol)；V0为气体摩尔体积,10-3m3/mol；E为杨氏模
3
量,MPa；ρr为干酪根的相对密度，g/m ；φads为吸附相的孔隙度，f；φorg为有机质中自由气所占据的孔隙度，f；CTOC为总有机碳含量，％；
[0106] 实例分析
[0107] 某页岩气藏的基本参数如下：Pi＝24.13MPa，Bgi＝4.82×10-3m3/m3，Smwi＝0.25，Sfwi＝0，Cm＝4.35×10-4MPa-1，Cf＝2.0×10-2MPa-1，Cw＝4.35×10-2MPa-1，φt＝0.112，VL＝11.32m3/t，b＝0.088，PL＝2.41MPa，ρsc＝7.7×10-4g/cm3，ρs＝0.34g/cm3，ρko＝1.325g/cm3，T＝366.49K，E＝26800MPa，ρb＝2.65g/cm3，M＝20g/mol，R＝8.134×10-3MPa·m3/(K·kmol),V0＝24.5m3/kmol，φads＝0.0023，φorg＝0.012，CTOC＝5％，假设临界解吸压力Pd＝
23.13MPa，生产数据如表1所示。
[0108] 表1某页岩气藏生产数据
[0109] 利用以上数据，采用本文推导的平衡方程进行气藏储量回归，结果如图3和表2。
[0110] 表2储量计算结果 108m3
[0111] 以上所述，并非对本发明作任何形式上的限制，虽然本发明已通过上述
实施例揭示，然而并非用以限定本发明，任何熟悉本专业的技术人员，在不脱离本发明技术方案范围内，当可利用上述揭示的技术内容作出些变动或修饰为等同变化的等效实施例，但凡是未脱离本发明技术方案的内容，依据本发明的技术实质对以上实施例所作的任何简单
修改、等同变化与修饰，均仍属于本发明技术方案的范围内。
