页岩气多场耦合渗透率计算模型赵谦平1王博涛2姜磊3付顺勋4.DOC

页岩气多场耦合渗透率计算模型 赵谦平1，王博涛2，姜磊3，付顺勋4，曹成1 （1.陕西延长石油（集团）有限责任公司，陕西 西安 710075； 2.中国石油长庆油田分公司，陕西 西安 710299； 3.陕西延长石油（集团）有限责任公司，陕西 延安 716000； 4.中国石油长庆油田分公司，陕西 西安 710018） 摘要：为表征扩散、应力敏感、基质收缩和吸附层对页岩气渗透率的影响，通过加权系数耦合过渡扩散和努森扩散，建立了扩散渗透率数学模型，并通过实验进行验证，同时分析了基质收缩应力敏感Schepers[3]和Zuber[4]等人页岩气主要方式为达西流Rushing[5]、Dahaghi[6]等人页岩气基质主要方式为扩散；Javadpour和Ozkan等人和扩散同时存在页岩气 页岩气渗透率包括自由气黏性流和扩散的贡献。自由气的扩散包括Fick扩散、过渡扩散以及努森扩散[21-22]，扩散类型可根据克努森数进行判断。通过质量流量和渗透率的关系，将Fick扩散、过渡扩散和努森扩散以渗透率的形式表达： （1） 式中：KF为Fick扩散渗透率，m2；KT为过渡扩散渗透率，m2；KKn为努森扩散渗透率，m2；φ为孔隙度；kB为Boltzmann常数，kB=1.38×10-23J/K；T为温度，K；π取3.14；rA为气体分子半径，m；μ为气体黏度，Pa·s；p为孔隙压力，pa；R为气体常数，8314Pa·m3·K-1·kmol-1；D为孔隙直径，m；M为气体分子量，kg/mol；Kn为克努森数。 图1为克努森数对扩散渗透率的影响。由图1a可知，努森扩散渗透率远大于Fick扩散渗透率和过渡扩散渗透率。由图1b可以看出，当Kn＜101时，Fick扩散渗透率和过渡扩散渗透率较为接近；当Kn＞102时，过渡扩散渗透率和努森扩散渗透率较为接近。由此可知，若仅使用努森扩散来描述自由气的扩散，将夸大扩散的作用。实际应用时，随着压力和孔径的变化，3种扩散均有可能出现，需要具体判断扩散类型，给计算带来不便。 a 通过压力变化改变克努森数 b 通过孔径变化改变克努森数 图1 克努森数对扩散渗透率的影响 过渡扩散渗透率和努森扩散渗透率处于2个边界，采用加权系数对过渡扩散渗透率和努森扩散渗透率进行加权求和： （2） 式中：f为加权系数；Km、、Km=7、a=4、b=16，得到具有唯一表达式的扩散渗透率数学模型： （3） 式中：为扩散渗透率，m2。 拟合效果见图2。由图2可知：当Kn＜10-1时，建立的数学模型与Fick扩散模型重合；当101＞Kn＞10-1时，建立的数学模型与过渡扩散模型重合；当Kn＞101时，建立的数学模型与克努森扩散模型重合。由此可见，建立的数学模型可以准确地描述3个扩散阶段。 图2 扩散渗透率拟合 2 渗透率影响因素分析 2.1 基质收缩和应力敏感的影响 当页岩气藏开发时，孔隙压力随气体的采出而逐渐降低。通常情况下，上覆岩层的压力是不变的，可形成上覆岩层压力和孔隙压力差，即应力敏感效应。当孔隙压力降至一定程度时，孔隙表面的吸附气开始解吸，引起孔隙壁面表面自由能的变化，该值的升高会使孔隙变大[23]，即基质收缩效应。应力敏感和基质收缩共同作用的孔隙度变化可表示为： （4） 式中：φst-shr为应力敏感和基质收缩影响后的孔隙度；φmi为初始基质孔隙度；pi为原始地层压力，Pa；pm为目前地层压力，Pa；VL为兰格缪尔体积，m3/kg；pL为兰格缪尔压力，Pa；ρs为页岩密度；kg/m3；E为杨氏模量，Pa；V0为气体分子体积，m3/mol；cmi为基质压缩系数，Pa-1。 应力敏感和基质收缩产生的孔径变化为： （5） 式中：Dst-shr为应力敏感和基质收缩影响后的基质孔径，m；Dmi为初始基质孔径，m。 2.2 吸附层的影响 根据Li[24]等人的研究，页岩气吸附层厚度随孔隙压力而变化： （6） 式中：h为吸附层厚度，m。 原始状态下，吸附气和自由气达到平衡，当页岩气藏开发后，孔隙压力下降，吸附层厚度改变，从而对页岩气渗透率产生影响[25]。 2.3 基质收缩、应力敏感、吸附层的耦合 假设页岩受基质收缩、应力敏感和气吸附层影响后的有效孔隙直径为Deff，吸附孔隙度为φads，有效孔隙度为φeff。考虑Deff=Dst-shr-2h和φeff=φst-shr-φads，可推导得到有效孔隙直径和有效孔隙度表达式： （7） （8） 式中：Deff为有