热-流-固耦合作用下页岩气储层渗透率的演化机制.pdf

DOI:10.11858/gywlxb.20230615热-流-固耦合作用下页岩气储层渗透率的演化机制张宏学1，刘卫群2,3（1.安徽理工大学力学与光电物理学院,安徽淮南232001；2.中国矿业大学力学与土木工程学院,江苏徐州221116；3.中国矿业大学深部岩土力学与地下工程国家重点实验室,江苏徐州221116）摘要：为了研究热-流-固耦合作用下页岩渗透率的演化机制，考虑热解吸、有效应力和热膨胀对页岩渗透率的影响，提出了页岩的有效应力-渗透率模型，该模型能够分析吸附应变和热膨胀应变对页岩渗透率的影响机制。基于该模型和多孔介质弹性理论，建立了单轴应变条件下页岩气储层的热解吸渗透率模型，该模型能够探讨页岩渗透率随温度和孔隙压力的演化规律。利用室内实验观测的页岩岩样渗透率实验数据，验证了该模型的有效性和准确性。结果表明：(1)热解吸渗透率模型能较好地拟合恒压变温条件下的 Marcellus 页岩渗透率。(2)探讨了恒温条件下页岩渗透率随孔压的演化机制，发现恒温条件下渗透率的演化规律呈“U 形”，温度越高，渗透率随孔压下降的反弹现象越不明显。(3)分析了恒压条件下页岩渗透率随温度的演化机制，发现恒压条件下渗透率随温度的演化规律呈“倒 U 形”，孔隙压力越大，温度对渗透率的影响越小。(4)分别在恒温和恒压条件下对热解吸渗透率模型进行敏感性分析，发现泊松比越大，渗透率比值梯度越大，孔隙体积模量越大，渗透率比值梯度越小。恒压条件下，当线胀系数大于临界值或朗缪尔体应变小于临界值，渗透率的演化规律不呈现“倒 U 形”。恒温条件下，当朗缪尔体应变小于临界值时，渗透率的演化规律不呈现“U 形”。关键词：热-流-固耦合；页岩渗透率；热解吸；有效应力；吸附应变中图分类号：O354.9;TP028.8文献标识码：A我国页岩气可采资源量为 2.5081013m31，四川盆地具有丰富的页岩气资源，已成为我国页岩气开发的主战场。四川盆地的页岩气储层埋深较大，其中龙马溪组的地层埋深普遍在 40004500m 之间，埋深较浅处约为 30004000m2。由于地下的地温梯度和静水压力梯度分别为 0.03K/m 和 0.01MPa/m3，所以龙马溪组的地层温度范围为 378423K（地表温度按 288K 计算）。Gensterblum 等3测定了甲烷在 303、318、333 和 348K（对应的储层深度分别为 500、1000、1500 和 2000m）时的吸附等温线，发现甲烷的吸附能力与温度和孔隙压力有关。在不同温度（373873K）和不同压力（0.115.0MPa）条件下，Geng 等4对油页岩进行了热解反应实验发现，随着温度的升高，孔隙率逐渐增大，裂隙的最大张开度也随之增大；最大孔隙率为室温下的 13.52 倍，裂隙总数为 813 条，裂隙最大孔径为 0.383mm。对于页岩储层，裂隙是气体流动的主要通道，气体解吸引起的体应变对页岩气的开采有着显著的影响。因此，建立热-流-固耦合作用下页岩的渗透率模型对准确预测深部页岩气的渗流规律尤为重要。目前，国内外学者主要通过理论分析和室内实验相结合的方法研究温度对页岩渗透率的影响。利用岩石三轴渗流仪，赵瑜等5研究了孔压和温度对页岩渗流的影响机理。利用自制的煤岩稳态实验装置，Wang 等6测试了我国无烟煤在不同气体压力和温度下的渗透率。吴迪等7利用三轴渗流实验装*收稿日期：2023-02-14；修回日期：2023-03-03 基金项目：安徽省教育厅科研基金（KJ2020A0329，KJ2016A207）作者简介：张宏学（1982），男，博士，副教授，主要从事非常规天然气开采中的关键力学问题研究.E-mail：第37卷第3期高压物理学报Vol.37,No.32023年6月CHINESEJOURNALOFHIGHPRESSUREPHYSICSJun.,2023035303-1置，研究了孔隙压力、体积应力和温度对页岩中 CH4气体渗流规律的影响。张道川等8采用多场耦合作用下不同相态 CO2致裂驱替 CH4实验装置，研究了有效应力、孔隙压力、温度以及吸附膨胀等因素对裂隙页岩渗透特性的影响。基于典型 S&D 模型，李波波等9建立了考虑温度和恒定外应力条件下的渗透率模型。Wang 等10采用高温三轴渗透率实验设备研究了现场应力条件下油页岩的各向异性渗透率。Schwartz 等11认为在恒定孔隙压力和温度变化下，页岩渗透率的演化规律可以用热解吸应变和页岩的热膨胀应变来描述，由此提出了热解吸渗透率增强模型。Teng 等12提出了变温度下煤气相互作用的热敏感性渗透率模型，该模型能够解释热膨胀、热压裂、基质吸附和裂隙表面的热挥发对煤渗透率的影响。在围压 1025MPa、温度 298343K、注入压力 612MPa 作用下，Ju 等13测试了煤样的 CO2渗透率，建立了表征裂隙煤 CO2渗透率的力学模型。Sinha 等14利用稳态渗透率测试设备研究了应力、压力、温度、吸附作用和气体类型对完整页岩岩样渗透率的影响。综上所述，国内外学者主要通过实验研究温度、围压、孔压和吸附等因素对页岩渗透率的影响，能够解释上述因素（特别是温度）对页岩渗透率影响的解析模型凤毛麟角。为此本研究将深部页岩储层简化为多孔弹性介质模型，考虑温度和吸附作用对渗透率的影响，推导页岩气开采过程中页岩储层的有效应力-渗透率模型，阐述围压、孔隙压力、有效应力、热解吸和温度对页岩渗透率的影响机理。然后，结合多孔介质弹性理论（考虑温度和解/吸附对变形的影响）建立页岩储层在单轴应变条件下的渗透率解析模型，利用实验数据验证该模型的有效性。最后，通过该模型分析不同温度下页岩渗透率的演化规律。1 温度作用下的有效应力-渗透率模型页岩的孔隙率可以表示为=VpVb(1)式中：Vp为页岩孔隙的体积，m3；Vb为页岩的体积，m3。对式(1)微分，获得d=dp+db(2)式中：dp为孔隙的体应变增量，db为页岩的体应变增量。考虑有效应力、气体的解吸附和温度对页岩变形的影响，db和 dp分别为db=dVbVb=1K(ddp)ds3TdT(3)dp=dVpVp=1Kpd(1Kp1Ks)dpds(4)ds式中：dVb和 dVp分别为页岩的体积增量和孔隙的体积增量，m3；K、Kp和 Ks分别为页岩的体积模量、孔隙体积模量和固体颗粒的体积模量，Pa；d 为平均法向应力增量，Pa；dp 为孔隙压力增量，Pa；为吸附引起的体应变增量；为 Biot 系数，=1K/Ks；T为页岩的线胀系数，K1；dT 为温度增量，K。联立式(2)、式(3)和式(4)，得到d=(1K1Kp)(ddp)3TdT(5)对式(5)积分0=exp(1K1Kp)(0)(p p0)3T(T T0)(6)0式中：为页岩的初始孔隙率；和 0分别为平均法向应力和初始平均法向应力，Pa；p 和 p0分别为孔隙压力和初始孔隙压力，Pa；T 和 T0分别为页岩气储层的温度和初始温度，K。页岩气储层的渗透率和孔隙率的关系15第37卷张宏学等：热-流-固耦合作用下页岩气储层渗透率的演化机制第3期035303-2kk0=(0)3(7)式中：k 和 k0分别为页岩气储层的渗透率和初始渗透率，m2。将式(6)代入式(7)，考虑页岩的体积模量远大于其孔隙体积模量，得到考虑温度、吸附作用的有效应力-渗透率模型kk0=exp3Kp(ee0)9T(T T0)(8)式中：e和 e0分别为有效应力和初始有效应力，单位 Pa，e=p，e0=0p0。该模型第 1 项为有效应力变化引起的体应变增量，第 2 项为温度（储层深度）变化引起的体应变增量。页岩气开采过程中储层渗透率的演化本质是储层的应变演化影响孔隙体积，该渗透率模型综合考虑了有效应力、温度对储层渗透率的影响。当不考虑温度对页岩气储层渗透率的影响时，式(8)可简化为文献 1617 推导的煤的有效应力-渗透率模型kk0=exp3Cp(ee0)(9)式中：Cp为孔隙压缩系数，单位 Pa1，Cp=1/Kp。2 单轴应变条件下的热解吸渗透率模型 2.1 单轴应变条件下的热解吸渗透率模型将页岩气储层简化为线弹性多孔介质，随着温度的升高，孔隙率和裂隙张开度均随之增大4。因此，考虑温度作用，规定压应力和压应变为正，页岩气储层的本构关系ij=kkij+2Gij+pij+Ksij+3KTTij(10)式中：ij为应力分量，单位 Pa；为拉梅常数，单位 Pa，=E/(1+)(12)；kk为体应变；ij为 Kronecker符号；G 为剪切模量，单位 Pa，G=E/2(1+)；E 为岩体的弹性模量，单位 Pa；为泊松比；ij为应变分量；为 Biot 系数，=1K/Ks；K 为体积模量，单位 Pa，K=E/3(12)；Ks为固体颗粒体积模量，单位 Pa。当考虑温度作用时，气体吸附/解吸引起的体应变利用下式计算11s=Lb0(T/T0)1/2expQ/(RT)p1+b0(T/T0)1/2expQ/(RT)p(11)b0=1pL1expQ/(RT0)式中：L为 Langmuir 体应变；b0为在无限大压力下的亲和常数，单位 Pa1，；pL为Langmuir 压力，单位 Pa；Q 为吸附热，单位 J/mol；R 为普适气体常数，单位 J/(molK)。联立式(8)、式(10)和式(11)，得到页岩气储层在单轴应变条件下的热解吸渗透率模型kk0=exp3Kp1+3(1)(p p0)2E9(1)(Lb0(T/T0)1/2exp(Q/(RT)p1+b0(T/T0)1/2exp(Q/(RT)pLb0exp(Q/(RT0)p01+b0exp(Q/(RT0)p0)exp3Kp2E3(1)+3KpT(T T0)(12)Kp=0K0=2b0/a0式中，18，1920。2.2 热解吸渗透率模型的验证文献 11 测试了 Marcellus 裂隙页岩和完整页岩岩样的渗透率比随温度的变化，分别见表 1 和表 2。实验所用气体为 CH4，实验过程中保持围压（24MPa）和孔隙压力（6MPa）不变。将式(12)热解吸渗透率模型的理论解和不同温度下 Marcellus 裂隙页岩的实验数据进行对比，当模型中的参数 T0=295K11、L=1.710318、b0=1.83104MPa111、Q=1.7104J/mol11、p=p0=6MPa、第37卷张宏学等：热-流-固耦合作用下页岩气储层渗透率的演化机制第3期035303-3R=8.31J/(molK)11、=0.35、E=32.75GPa21、Kp=24.26MPa、T=1.8106K1时，该模型的理论解与Marcellus 裂隙页岩的实验数据对比如图 1 所示。图 1 中，当温度从 295K 升高到 323K 时，Marcellus 裂隙页岩的渗透率增大了 21.1%。由于黏土和有机质的解吸以及岩桥的热膨胀，裂隙页岩的渗透率比值随着温度的升高而逐渐增大。由此可见，该热解吸渗透率模型的理论解与 Marcellus 裂隙页岩的实验数据基本一致。将式(12)的热解吸渗透率模型理论解与不同温度下 Marcellus 完整页岩（含天然裂隙）的实验数据进行对比，当模型中 T0=295K11、L=1.710318、b0=1.83104MPa111、Q=1.7104J/mol11、p=p0=6MPa、R=8.31J/(molK)11、=0.35，E=32.75GPa21、Kp=43.53MPa、T=1.8106K1时，该模型的理论解与 Marcellus 完整页岩的实验数据对比见图 2。由图 2 可知，随着温度的升高，完整页岩的渗透率比值呈反弹趋势，且整体逐渐增大，这是由于热解吸引起的孔隙膨胀量大于固体颗粒膨胀引起的孔隙闭合量。该渗透率模型的解析解能够与不同温度下 Marcellus 完整页岩的实验数据较好地吻合。2953003053103153203251.001.051.101.151.201.25Temperature/KPermeability ratioMarcellus fractured shaleProposed model图1不同温度下热解吸模型解析解与裂隙页岩的实验数据对比Fig.1Comparisonofanalyticalresultsofthermal-sorptivemodelandexperimentaldataoffracturedshaleunderdifferenttemperatures3003103203303401.201.151.101.051.000.95Temperature/KMarcellus intact shaleProposed modelPermeability ratio图2不同温度下热解吸模型解析解与完整页岩的实验数据对比Fig.2Comparisonofanalyticalresultsofthermal-sorptivemodelandexperimentaldataofintactshaleunderdifferenttemperatures表 1 裂隙页岩的渗透率比Table 1 Permeability ratio of fractured shaleTemperature/Kk/k0offracturedshale2951.0003031.0443131.1133231.211表 2 完整页岩的渗透率比Table 2 Permeability ratio of intact shaleTemperature/Kk/k0ofintactshaleTemperature/Kk/k0ofintactshale307.11.000332.71.095308.00.978333.31.045313.81.051337.71.058321.11.037337.91.061321.21.008338.01.084325.91.085第37卷张宏学等：热-流-固耦合作用下页岩气储层渗透率的演化机制第3期035303-4 3 温度作用下渗透率的演化机理 3.1 渗透率随温度和孔压变化的演化机理由式(12)计算得到，页岩储层在单轴应变条件下的渗透率随温度和孔隙压力的演化规律如图 3 所示，储层初始温度为 300K，初始孔压为 6MPa。如图 3 所示，当储层温度介于 300420K 范围内且保持不变时，随着孔隙压力的增大，渗透率比值先下降然后反弹，这是因为温度不变时，影响渗透率的因素主要为有效应力和气体吸附/解吸，随着孔隙压力的增大，有效应力减小，孔隙体积变大，导致渗透率增大；另一方面，孔隙压力增大，气体吸附量增大，孔隙体积减小，导致渗透率减小。孔隙压力增大前期，有效应力引起的渗透率增大量小于气体吸附引起的渗透率减小量，导致渗透率下降；孔隙压力增大后期，有效应力引起的渗透率增大量大于气体吸附引起的渗透率减小量，从而导致渗透率反弹。温度越高，渗透率比值下降越平缓，反弹越急剧。当储层温度介于 420440K 范围内，且保持不变时，随着孔压的增大，渗透率逐渐增大。这是因为温度较高时，随着孔隙压力的增大，有效应力引起的渗透率增大量始终大于气体吸附引起的渗透率减小量，导致渗透率逐渐增大。从图 3 还可以看出，当孔隙压力保持不变时，随着温度升高，渗透率比值先增大然后逐渐减小；孔隙压力越大，渗透率比值开始下降时的温度越大。这是因为当孔隙压力不变时，影响渗透率的主要因素为气体热解吸/热吸附和岩体热膨胀，随着温度升高，气体热解吸量增大，孔隙体积增大，导致渗透率增大；另一方面，随着温度升高，页岩固体颗粒膨胀，孔隙体积缩小，导致渗透率减小。温度升高前期，热解吸引起的渗透率增大量大于热膨胀引起的渗透率减小量，从而导致渗透率增大；而温度升高后期，热解吸引起的渗透率增大量小于热膨胀引起的渗透率减小量，导致渗透率减小。3.2 不同孔压作用下渗透率随温度变化的演化机理不同孔隙压力作用下，页岩气储层渗透率随温度变化的演化规律如图 4 所示。当孔隙压力保持不变（5、15和 25MPa）时，随着温度升高，渗透率比值先增大，然后逐渐减小，孔隙压力越小，热解吸引起的渗透率增大的速率越大，导致渗透率比值开始减小时的温度越小；当渗透率比值开始减小以后，孔隙压力越小，热膨胀引起的渗透率减小的速率越大，从而导致渗透率比值减小的越急剧。当孔隙压力保持 35 和 45MPa 不变时，随着温度升高，渗透率比值轻微下降，然后反弹，最后逐渐下降。孔隙压力较小时，反弹现象较为明显。从图 4 可以看出，当温度在 300440K 范围内变化，孔隙压力为 5MPa 时，渗透率比值变化范围为 0.71.2。而孔隙压力为 5 和 40MPa 时，渗透率比值在 0.91.1 范围内波动。由此可见，孔隙压力越大，温度对储层渗透率的影响越小。051015203003504004500.51.01.52.02.5Pressure/MPaTemperature/KPermeability ratiop0=6 MPa T0=300 K图3渗透率随温度和孔压的演化规律Fig.3Evolutionofpermeabilitywithtemperatureandporepressure3003203403603804004204400.850.900.951.001.051.101.151.201.251.301.35Temperature/K Permeability ratiop=5 MPa p=15 MPap=25 MPap=35 MPap=45 MPa图4不同孔压下渗透率随温度的演化规律Fig.4Evolutionofpermeabilitywithtemperatureunderdifferentporepressures第37卷张宏学等：热-流-固耦合作用下页岩气储层渗透率的演化机制第3期035303-5 3.3 不同温度作用下渗透率随孔压的演化机理不同温度作用下，页岩气储层渗透率随孔压的演化规律如图 5 所示。当储层温度保持不变（318、348、378 和 408K，对应的储层深度分别为1000、2000、3000 和 4000m）时，随着孔隙压力的减小，渗透率比值先减小然后逐渐增大，其演化规律呈反弹现象。储层温度（储层深度）越大，渗透率的反弹现象越不明显，当储层深度达到 3000m时，渗透率基本不反弹。当储层温度均保持不变时，影响储层渗透率的主要因素为有效应力和孔压降低引起的气体解吸。当孔隙压力下降时，一方面吸附于有机质颗粒表面和孔隙表面的气体解吸，孔隙体积和裂隙张开度增大，导致储层渗透率增大；另一方面，孔隙压力下降导致储层有效应力增大，孔隙和裂隙被压缩，导致储层渗透率减小。当气体解吸引起的基质收缩量小于有效应力增大引起的孔隙和裂隙的压缩量时，渗透率下降；当气体解吸引起的基质收缩量大于有效应力增大引起的孔隙和裂隙压缩量时，渗透率增大，这也正是恒温作用下储层渗透率随着孔隙压力的下降呈现反弹现象的原因。4 讨论在式(12)所示的热解吸渗透率模型中，泊松比、孔隙体积模量、热膨胀系数和 Langmuir 体应变对储层渗透率的影响较为显著，下面将分别讨论恒压、恒温作用下上述参数对储层渗透率的影响。4.1 恒压下渗透率的敏感性分析当孔隙压力保持 10MPa 不变时，泊松比对储层渗透率的影响如图 6 所示。可以看出，泊松比不变时，随着温度升高，储层渗透率比值先增大然后减小，呈现“倒 U 形”演化规律。泊松比越大，渗透率比值的梯度越大。究其原因，储层孔隙压力不变，热解吸/吸附和热膨胀是影响页岩储层渗透率的主要因素。当温度升高时，页岩岩体膨胀导致储层孔隙和裂隙体积压缩，使得储层渗透率下降；另外，由于温度升高，页岩的吸附能力逐渐降低22，从而导致储层孔隙体积逐渐增大，裂隙张开度随之增大4，储层渗透率也随之增大。综上所述，热解吸和热膨胀对页岩储层渗透率的影响是相反的。温度升高初期，热解吸引起的渗透率变化量大于热膨胀引起的渗透率变化量；温度升高后期，热解吸引起的渗透率变化量小于热膨胀引起的渗透率变化量，从而导致恒压下渗透率随温度呈现“倒 U 形”变化。泊松比越大，热解吸引起的渗透率变化量越大，热膨胀引起的渗透率变化量越小，因此渗透率比值梯度随之增大。当孔隙压力保持 10MPa 不变时，孔隙体积模量对储层渗透率的影响见图 7。如图 7 所示，当孔隙压力恒定时，储层渗透率比值随温度升高仍然呈现“倒 U 形”的演化规律，而且随着孔隙体积模量的增大，储层渗透率比值梯度减小。这是因为孔隙体积模量越大，热膨胀引起的渗透率变化量越大，而热解05101520250.40.50.60.70.80.91.0Pressure/MPaPermeability ratioT=318 KT=348 KT=378 KT=408 K图5不同温度下渗透率随孔压的演化规律Fig.5Evolutionofpermeabilitywithporepressureunderdifferenttemperatures3003504004501.01.11.21.31.4Temperature/KPermeability ratio=0.30=0.35=0.40图6恒压下泊松比对渗透率的影响Fig.6EffectofPoissonsratioonpermeabilityunderconstantporepressure第37卷张宏学等：热-流-固耦合作用下页岩气储层渗透率的演化机制第3期035303-6吸引起的渗透率变化量越小，从而导致渗透率比值梯度随孔隙体积模量的增大而减小。当储层的孔隙压力不变时，页岩的线胀系数对储层渗透率的影响见图 8。如图 8 所示，当线胀系数小于临界值（渗透率随温度的演化规律不再呈现“倒 U 形”时的线胀系数）时，渗透率比值随温度的演化规律呈现“倒 U 形”现象，可通过式(12)计算出该临界值；当线胀系数大于该临界值时，渗透率比值随温度的升高而逐渐下降，不再呈现“倒U 形”。在温升初期，线胀系数越大，渗透率比值梯度越大；在温升后期，渗透率比值梯度随线胀系数的增大而减小。这是因为孔隙压力不变时，只有气体的热解吸和岩体的热膨胀影响储层渗透率，线胀系数越大，页岩岩体因温度升高而引起的体积膨胀量越大，挤压的孔隙体积和裂隙体积越大，渗透率下降的越快。当热膨胀引起的渗透率变化量总是大于热解吸引起的渗透率变化量时，渗透率比值随温度的升高呈现逐渐下降的趋势。恒压作用下朗缪尔体应变对储层渗透率的影响见图 9。如图 9 所示，当朗缪尔体应变大于临界值（渗透率随温度的演化规律不呈现“倒 U 形”时的朗缪尔体应变）时，储层渗透率随温度的变化曲线呈“倒 U 形”趋势，可以通过式(12)计算出此临界值；当朗缪尔体应变小于临界值时，渗透率随着温度的升高而逐渐降低。在温升初期，朗缪尔体应变越小，渗透率比值梯度越大；而在温升后期，朗缪尔体应变越小，渗透率比值梯度越小。这是因为朗缪尔体应变发生变化时，热解吸引起的体应变随之变化，热膨胀引起的体应变不发生变化，因此当升温过程中热解吸引起的渗透率变化量始终小于热膨胀引起的渗透率变化量时，渗透率比值随温度的升高逐渐下降。朗缪尔体应变越大，温升初期热解吸引起的渗透率变化量越大，此渗透率比值梯度越小，温升后期热解吸引起的渗透率变化量越小，渗透率梯度比值越大。4.2 恒温下渗透率的敏感性分析当储层温度不变时，泊松比对储层渗透率的影响见图 10。如图 10 所示，恒温下储层渗透率随孔隙压力的下降，先逐渐降低然后反弹逐渐增大。泊松比越大，有效应力变化引起的孔隙和裂隙的体积变化越大，从而导致渗透率变化量越大；泊松比越大，气体解吸引起的孔隙和裂隙的体积变化越小，进而导致渗透率变化量越小。因此，泊松比越大，渗透率比值梯度越大。当储层温度不变时，孔隙体积模量对储层渗透率的影响见图 11。如图 11 所示，当储层温度保持不3003504004501.01.21.41.61.8Temperature/K Kp=15 MPaKp=20 MPaKp=25 MPaPermeability ratio图7恒压下孔隙体积模量对渗透率的影响Fig.7EffectofporevolumemodulusonpermeabilityunderconstantporepressurePermeability ratio30035040045000.20.40.60.81.01.21.4Temperature/KT=5.8106 K1T=3.8106 K1T=1.8106 K1图8恒压下线胀系数对渗透率的影响Fig.8Evolutionofcoefficientoflinearexpansiononpermeabilityunderconstantporepressure30035040045000.51.01.52.0Temperature/KL=2.0103L=2.0104Permeability ratioL=1.1103图9恒压下 Langmuir 体应变对渗透率的影响Fig.9EvolutionofLangmuirvolumestrainonpermeabilityunderconstantporepressure第37卷张宏学等：热-流-固耦合作用下页岩气储层渗透率的演化机制第3期035303-7变时，渗透率比值随孔隙压力的减小均呈现先变小再反弹增大的趋势。孔隙体积模量越大，孔隙抵抗变形的能力越强，相同的有效应力变化量引起的孔隙体积变化量越小，所以有效应力变化引起的渗透率变化量越小。另外，孔隙体积模量越大，气体解吸引起的渗透率变化量越小。因此，孔隙体积模量越大，渗透率比值梯度越小。当储层温度不变时，朗缪尔体应变对渗透率的影响见图 12。如图 12 所示，朗缪尔体应变越大，渗透率随孔隙压力变化呈现的先变小再反弹现象越明显，随着朗缪尔体应变的减小，反弹现象逐渐消失。这是因为有效应力引起的渗透率变化量与朗缪尔体应变无关，而气体解吸引起的渗透率变化量取决于朗缪尔体应变。朗缪尔体应变越小，气体解吸引起的渗透率变化量越小，从而导致在孔隙压力的下降过程中有效应力变化引起的渗透率变化量（减小）始终大于气体解吸引起的渗透率变化量（增大）。5 结论(1)考虑温度、气体压力、有效应力和热解吸对储层渗透率的影响，推导了页岩气储层的有效应力-渗透率模型，该模型能解释有效应力不变时渗透率随孔隙压力的演化机制。(2)建立了储层在单轴应变条件下的热解吸渗透率模型，该模型能准确地阐述 Marcellus 裂隙页岩和完整页岩岩样的渗透率随温度的演化机制。(3)热解吸/吸附、热膨胀和孔隙压力是影响页岩气储层渗透率的主要因素。当孔隙压力恒定时，渗透率随温度呈“倒 U 形”演化规律。孔隙压力越大，温度对渗透率的影响越小；当温度恒定时，渗透率随孔压变化呈“U 形”演化规律，温度（对应储层深度）越高，渗透率随孔压下降的反弹现象越不明显。(4)储层孔压恒定时，泊松比越大，渗透率比值梯度越大；孔隙体积模量越大，渗透率比值梯度越小；当线胀系数大于临界值时，渗透率比值随温度的演化规律不呈现“倒 U 形”；当朗缪尔体应变小于临界值时，渗透率比值随温度的演化规律也不呈现“倒 U 形”。(5)储层温度恒定时，泊松比越大，渗透率比值梯度越大；孔隙体积模量越大，渗透率比值梯度越小。当朗缪尔体应变小于临界值时，渗透率比值随孔压的变化不呈现“U 形”演化规律。05101520250.40.50.60.70.80.91.0Pressure/MPa=0.30=0.35=0.40Permeability ratio图10恒温下泊松比对渗透率的影响Fig.10EffectofPoissonsratioonpermeabilityunderconstanttemperature05101520250.20.40.60.81.0Pressure/MPaKp=15 MPaKp=20 MPaKp=25 MPaPermeability ratio图11恒温下孔隙体积模量对渗透率的影响Fig.11Effectofporevolumemodulusonpermeabilityunderconstanttemperature51015202500.20.40.60.81.0Pressure/MPaL=2.0104L=2.0103L=1.1103Permeability ratio图12恒温下 Langmuir 体应变对渗透率的影响Fig.12EffectofLangmuirvolumestrainonpermeabilityunderconstanttemperature第37卷张宏学等：热-流-固耦合作用下页岩气储层渗透率的演化机制第3期035303-8参考文献：张金川,林腊梅,李玉喜,等.页岩气资源评价方法与技术:概率体积法J.地学前缘,2012,19(2):184191.ZHANGJC,LINLM,LIYX,etal.Themethodofshalegasassessment:probabilityvolumemethodJ.EarthScienceFrontiers,2012,19(2):184191.1李敏,庞雄奇,罗冰,等.生烃潜力法在深层页岩气资源评价中的应用以四川盆地五峰龙马溪组优质烃源岩为例J.中国矿业大学学报,2021,50(6):10961107.LIM,PANGXQ,LUOB,etal.Applicationofhydrocarbongenerationpotentialmethodtodeepshalegasresourceevaluation:a case study of high-quality source rocks of the Wufeng-Longmaxi formation in the Sichuan Basin J.Journal of ChinaUniversityofMining&Technology,2021,50(6):10961107.2GENSTERBLUM Y,MERKEL A,BUSCH A,et al.Gas saturation and CO2 enhancement potential of coalbed methanereservoirsasafunctionofdepthJ.AAPGBulletin,2014,98(2):395420.3GENGYD,LIANGWG,LIUJ,etal.EvolutionofporeandfracturestructureofoilshaleunderhightemperatureandhighpressureJ.Energy&Fuels,2017,31(10):1040410413.4赵瑜,王超林,曹汉,等.页岩渗流模型及孔压与温度影响机理研究J.煤炭学报,2018,43(6):17541760.ZHAOY,WANGCL,CAOH,etal.InfluencingmechanismandmodellingstudyofporepressureandtemperatureonshalepermeabilityJ.JournalofChinaCoalSociety,2018,43(6):17541760.5WANGK,DUF,WANGGD,etal.Investigationofgaspressureandtemperatureeffectsonthepermeabilityandsteady-statetimeofChineseanthracitecoal:anexperimentalstudyJ.JournalofNaturalGasScienceandEngineering,2017,40:179188.6吴迪,王挺,刘雪莹,等.页岩渗透特性受热力条件影响的实验研究J.实验力学,2020,35(3):539546.WU D,WANG T,LIU X Y,et al.Experimental study on the influence of thermal conditions on shale permeabilitycharacteristicsJ.JournalofExperimentalMechanics,2020,35(3):539546.7张道川,周军平,鲜学福,等.多场耦合作用下页岩渗透特性实验研究J.地下空间与工程学报,2018,14(3):613621.ZHANGDC,ZHOUJP,XIANXF,etal.Experimentstudyonthecouplingmulti-fieldeffectonthedynamicvariationofpermeabilityinshaleJ.ChineseJournalofUndergroundSpaceandEngineering,2018,14(3):613621.8李波波,高政,杨康,等.考虑温度、孔隙压力影响的煤岩渗透性演化机制分析J.煤炭学报,2020,45(2):626632.LIBB,GAOZ,YANGK,etal.AnalysisofcoalpermeabilityevolutionmechanismconsideringtheeffectoftemperatureandporepressureJ.JournalofChinaCoalSociety,2020,45(2):626632.9WANGGY,YANGD,ZHAOYS,etal.ExperimentalinvestigationonanisotropicpermeabilityanditsrelationshipwithanisotropicthermalcrackingofoilshaleunderhightemperatureandtriaxialstressJ.AppliedThermalEngineering,2019,146:718725.10SCHWARTZB,ELSWORTHD.InvertedU-shapedpermeabilityenhancementduetothermallyinduceddesorptiondeterminedfromstrain-basedanalysisofexperimentsonshaleatconstantporepressureJ.Fuel,2021,302:121178.11TENGT,WANGJG,GAOF,etal.Athermallysensitivepermeabilitymodelforcoal-gasinteractionsincludingthermalfracturingandvolatilizationJ.JournalofNaturalGasScienceandEngineering,2016,32:319333.12JU Y,WANG J G,WANG H J,et al.CO2 permeability of fractured coal subject to confining pressures and elevatedtemperature:experimentsandmodelingJ.ScienceChinaTechnologicalSciences,2016,59(12):19311942.13SINHAS,BRAUNEM,DETERMANMD,etal.Steady-statepermeabilitymeasurementsonintactshalesamplesatreservoirconditions-effectofstress,temperature,pressure,andtypeofgasC/SPEMiddleEastOilandGasShowandConference.Manama,Bahrain:SPE,2016.14张宏学,刘卫群.非平衡解吸状态下页岩气储层渗透率演化机制J.岩土力学,2021,42(10):26962704.ZHANGHX,LIUWQ.Permeabilityevolutionmechanismofshalegasreservoirinnon-equilibriumdesorptionstateJ.RockandSoilMechanics,2021,4