真三向应力作用下深部储层砂岩渗透率各向异性实验研究

渗透率各向异性是沉积岩层理结构中一个非常典型的现象,一方面它由原生沉积结构决定,即原生各向异性,另一方面它受到总应力和孔隙压力影响,即诱发各向异性.为研究真三向应力条件下储层砂岩渗透率原生与诱发各向异性,以中国东北部S6储气库储层砂岩为研究对象,采用东北大学自主研制的硬岩真三轴应力-渗流耦合装置对储层砂岩进行渗流实验,通过稳态法完成同一砂岩3个相互垂直方向的渗透率测试.实验结果表明:储层砂岩在施加的应力和孔隙压力范围内,平行层理方向的渗透率k x为100.94～113.98 m D,k y为98.34～111.41 mD,垂直层理方向的渗透率k z为54.98～63.29 mD;储层砂岩3个正交方向渗透率均随主应力增加而减少,随孔隙压力加载而递增;与气体渗流方向垂直的应力对渗透率的影响大于与气体渗流方向平行的应力对渗透率的影响;当外部应力方向都垂直于气体渗流方向时,与层理垂直的应力对渗透率的影响大于与层理平行的应力对渗透率的影响;孔隙压力对储层砂岩渗透率的线弹性响应并不是各向同性的,孔隙压力对水平层理方向所产生的渗透率增量超过了垂直层理方向.研究结论为地下储气库储层砂岩渗透率准确预测提...     

波尔兹曼数字岩芯致密砂岩渗透率研究

致密砂岩渗透率在油气勘探开发、地应力测量及水库地质灾害等领域具有重要意义,但鉴于目前实验技术局限性,温压耦合渗透率测量尚无法通过实验手段实现. 在已有研究工作基础上,首次建立了基于D3Q27数字岩芯模型,并计算了高温压耦合低渗砂岩渗透率问题. 首先,以鄂尔多斯盆地某油田延长组致密砂岩为例,利用X射线CT断层成像技术岩芯获取10μm, 5μm, 2μm分辨率致密砂岩内部结构数据,应用基于量子力学第一性原理的D3Q27格子波尔兹曼数字岩芯模型建立数值模型. 进而,利用自编3DLBM程序分别计算了不同分辨率渗透率随围压(0～200MPa)、孔隙压(0～65MPa)和温度(25℃～180℃)变化规律,通过与Inc AUTOLAB2 000C岩石测试分析系统实验结果对比,验证了程序的可靠性,得到低渗砂岩断层最佳分辨率;最后,在并行CPU-GPU平台上计算了高温高压耦合(0℃～400℃, 0～1.4 GPa)下致密砂岩渗透率值及其各向异性随温压变化规律,并讨论了致密砂岩中水在达到超临界状态后对致密岩石内部结构的影响.     

页岩黏土孔隙含水饱和度分布及其对甲烷吸附的影响

考虑储层原始含水特征,甲烷在页岩的吸附特征属于气液固三相复杂作用结果,水分在很大程度上影响页岩吸附能力,将成为制约页岩气资源量评估可靠性的主要原因之一.鉴于页岩水分主要分布于黏土等无机矿物孔隙内部,分析了甲烷-水膜-页岩黏土三相作用特征,结果表明:甲烷在干燥黏土表面吸附满足气固界面Langmuir吸附特征,在黏土水膜表面吸附满足气液界面Gibbs吸附特征,在气液固三相作用下满足"气固"与"气液"界面混合吸附特征;同时研究还发现:不同尺度孔隙内含水饱和度分布特征存在差异,部分小孔隙可以被水分充满,而大孔隙仅吸附一定厚度水膜.因此,水分对甲烷吸附能力的影响主要表现为两个方面:小孔隙被水分阻塞而失去吸附能力;大孔隙表面水膜改变甲烷吸附特征(气固界面吸附转变为气液界面吸附),以黏土样品为例,两者综合效应可以致使甲烷吸附能力降低约90%.从微观角度揭示了水分对页岩吸附能力的影响机理,将为建立合理评价页岩吸附气含量的方法奠定理论基础.     

页岩储层纳微米孔隙CO2/CH4吸附及驱替特性研究进展

利用CO₂开采页岩气不仅能够提高页岩气采收率, 还能够节省水资源并且对CO₂进行地质封存, 有助于实现页岩气开采过程的碳中和. 富有机质页岩储层纳微米孔隙中气体运移机制不同于常规储层, CO₂在储层中具有超临界特性, 致使开采机理复杂, 无法得到CO₂开采页岩气微观机理的准确认识, 所以研究CH₄, CO₂及其二元混合物在页岩储层纳微米孔隙中的吸附及驱替特性对准确评估和高效开采页岩气至关重要. 本文从实验、理论以及模拟方面对页岩储层纳微米孔隙中CH₄的吸附特性、CO₂/CH₄二元混合物竞争吸附特性以及驱替特性进行了综合分析, 对气体在纳微米孔隙中吸附及驱替特性的基础研究及关键问题进行讨论分析并提出了展望. 研究表明CH₄在页岩储层中表现为物理吸附, 有机质特征(丰度、成熟度、类型)、孔隙结构、无机矿物组成、温度和压力、含水率对页岩的CH₄吸附能力均有一定程度的影响. 在相同条件下, CO₂比CH₄更易被页岩储层吸附, 在页岩储层中注入CO₂可以促进CH₄的解吸, 并有利于CO₂的地质埋存. 开采方案的部署可采用井网形式的注采方式, 可以通过调整注入井的位置、数量以及CO₂注入速率对开采方案进行优化.      

页岩气和致密砂岩气藏微裂缝气体传输特性

页岩气和致密砂岩气藏发育微裂缝,其开度多在纳米级和微米级尺度且变化大,因此微裂缝气体传输机理异常复杂.本文基于滑脱流动和努森扩散模型,分别以分子之间碰撞频率和分子与壁面碰撞频率占总碰撞频率的比值作为滑脱流动和努森扩散的权重系数,耦合这两种传输机理,建立了微裂缝气体传输模型. 该模型考虑微裂缝形状和尺度对气体传输的影响. 模型可靠性用分子模拟数据验证.结果表明:(1)模型能够合理描述微裂缝中所有气体传输机理,包括连续流动,滑脱流动和过渡流动;(2)模型能够描述不同开发阶段,微裂缝中各气体传输机理对传输贡献的逐渐变化过程;(3)微裂缝形状和尺度影响气体传输,相同开度且宽度越大的微裂缝,气体传输能力越强,且在高压和微裂缝大开度的情况下表现更明显.      

页岩气井组分比例变化规律研究

页岩气开发过程中,生产井产出气的组分比例会随时间发生变化.本文基于组分模型数值模拟研究了生产井中甲烷组分比例变化的规律.研究表明,吸附气、渗透率与孔隙度影响页岩气组分比例的瞬态响应特征. 吸附气显著影响组分比例的变化规律,吸附量的大小决定组分比例的变化值及组分比例导数曲线的上下位置. 渗透率影响组分比例初期变化规律,但在后期,不同渗透率对瞬态组分比例规律的影响基本一致.孔隙度对组分比例变化及其导数曲线的影响与吸附气的影响类似,但在生产初期,孔隙度对组分比例的影响要小于吸附气对组分比例的影响. 本文的研究提供了一种进行页岩地层参数评价的新方法.      

考虑多重运移机制耦合页岩气藏压裂水平井数值模拟

页岩作为典型的微纳尺度多孔介质,游离气与吸附气共存,传统的达西定律已无法准确描述气体在页岩微纳尺度的运移规律.基于双重介质模型和离散裂缝模型构建页岩气藏分段压裂水平井模型,其中基岩中考虑气体的黏性流、Knudsen 扩散以及气体在基岩孔隙表面的吸附解吸,吸附采用Langmuir等温吸附方程;裂缝中考虑黏性流和Knudsen扩散,在此基础上建立基岩-裂缝双重介质压裂水平井数学模型并采用有限元方法对模型进行求解.结果表明,基岩固有渗透率越小,表面扩散和Knudsen扩散的影响越大,反之则越小;人工裂缝的性质包括条数、开度、半长以及间距,主要影响压裂水平井生产早期,随着人工裂缝参数值的增加,压裂水平井产能增加,累产气量也越大.其次,页岩气藏压裂诱导缝和天然裂缝的发育程度对页岩气藏的产能有很大的影响,水平井周围只有人工裂缝,周围天然裂缝不开启或不发育时,页岩气藏的水平井的产能较低.      

页岩气藏压裂水平井试井分析

页岩气藏资源丰富,开发潜力巨大,已成为目前研究的热点.与常规气藏相比,页岩气藏运移机制复杂,流动模式呈非线性,有必要考虑页岩气的吸附解吸,天然微裂缝的应力敏感性,人工裂缝内的非达西流等非线性因素对压裂水平井压力响应的影响. 基于双重介质和离散裂缝混合模型,分别采用Langmuir等温吸附方程描述吸附解吸,渗透率指数模型描述应力敏感,Forchheimer方程描述非达西效应,建立页岩气藏压裂水平井数值试井模型. 运用伽辽金有限元法对模型进行求解.根据试井特征曲线,划分流动阶段,着重分析非线性因素对压力响应的影响.结果表明:页岩气藏压裂水平井存在压裂裂缝线性流、压裂裂缝径向流、地层线性流、系统径向流及封闭边界影响5 种流动阶段.吸附解吸的影响发生窜流之后,Langmuir吸附体积增大,拟压力导数曲线凹槽更加明显,系统径向流出现时间与压力波传播到边界时间均延迟;天然裂缝系统的应力敏感性主要影响试井曲线的晚期段,拟压力和拟压力导数曲线均表现为上翘,应力敏感效应越强,上翘幅度越大;高速非达西效应对早期段影响较大,非达西效应越强,拟压力降幅度越大,试井曲线上翘.与解析解的对比以及矿场实例验证了模型的正确性与适用性.      

页岩气致密储层渗流模型及压裂直井产能分析

页岩气储层中存在大量的纳微米孔隙,其中气体的流动规律不同于常规气藏. 文章考虑克努森扩散及解吸作用的影响,建立页岩气稳态条件下产能公式,利用渗流阻力法得到了页岩气储层压裂直井产能方程.结合生产实例,进行了数值模拟计算. 结果表明：当裂缝导流能力达到0.12μm2·cm 后,产气量增幅减小,由此对裂缝导流能力进行优化. 游离气产量占总产气量的85%～90%,游离气对总产气量贡献较大. 该模型为页岩气产能预测及开发指标优化提供了理论依据.     

页岩基质解吸-扩散-渗流耦合实验及数学模型

采用川南地区龙马溪组页岩样品,设计了页岩基质解吸-扩散-渗流耦合物理模拟实验,揭示了页岩基质气体流动特征以及压力传播规律.推导了页岩气解吸-扩散-渗流耦合数学模型并且利用有限差分法对数学模型进行数值求解,与实验结果相比较表明该数学模型能够很好地描述气体在页岩基质中的流动规律.同时对页岩基质气体流动的影响因素进行了分析,认为页岩基质的渗透率、扩散系数、解吸附常数等因素均能影响页岩基质气体的流量和压力传播规律,在页岩气藏的开发过程中需要考虑这些参数的影响,该数学模型为页岩气井产能计算提供了更准确的计算方法.      

不同滑移边界下的页岩渗透率修正模型

页岩中的孔隙直径通常为纳米量级,基于连续流的达西定律已不能描述纳米级孔隙内的气体流动规律,一般采用附加滑移边界条件的Navier-Stokes方程对其进行描述. 由此可导出与压力相关的渗透率公式(称为"视渗透率"),并用来修正达西定律.因而,渗透率修正方法研究成为页岩气流动研究的热点之一.首先,基于Hagen-Poiseuille 流推导出一般形式二阶滑移模型下的速度分布和流量公式,并推导出相应的渗透率修正公式.该渗透率修正公式基本能将现有的滑移速度模型统一表达为对渗透率的修正. 基于一般形式的渗透率修正公式,重点研究了Maxwell, Hsia, Beskok与Ng 滑移模型速度分布渗透率修正系数、及其对井底压力的影响;提出了基于Ng 滑移速度模型的渗透率修正公式. 基于页岩实际储层温压系统及孔隙分布,计算了Kn 范围及储层条件下页岩气的流动形态,表明页岩气流动存在滑移流、过渡流与分子自由流. 而Ng 模型能描述Kn<88 的滑移流、过渡流、自由分子流的流量规律,因此可以用于描述页岩实际储层中页岩气的流动特征. 计算表明,随着Kn 的增加,不同滑移模型下的渗透率修正系数差异增大.Maxwell与Hsia模型适用于滑移流与过渡流早期,Beskok与Ng 模型可描述自由分子流下的流动规律,但二者在虚拟的孔径均为10nm页岩中,井底压力的差别开始显现;在虚拟的孔径均为1nm页岩中,井底压力的差别开始明显.      

页岩气储层压裂数值模拟技术研究进展

页岩气储层水力压裂数值模拟既要考虑页岩储层岩石的特性,又要兼顾水平井分段压裂施工工艺,是一个非常棘手的力学难题.本文简述了页岩气储层岩石具有的地质力学特征和页岩气储层开发常用的水平井分段压裂技术;详述了扩展有限元、边界元、离散元在水力压裂裂缝模拟上的应用现状,指出了它们在处理裂缝问题的局限性和优越性,总结出边界元三维位移不连续法是模拟多裂缝扩展的有效方法.     

页岩储层体积压裂的地应力变化研究

页岩储层属于致密超低渗透储层,需改造形成复杂缝网才有经济产能.体积压裂是页岩储层增产改造的主要措施,而地应力场特别是水平主应力差值是体积压裂的关键控制因素. 理论研究表明:(1)当初始两向水平主应力差较小时,容易形成缝网,反之不易产生缝网;(2)人工裂缝的形成能够改变地层初始应力场. 因此应在前人研究的基础上优化设计压裂方式,以克服和翻转初始水平主应力差值,产生体积缝网.基于此,建立了页岩气藏水平井体积压裂数值模型,模型中采用多孔介质流固耦合单元模拟页岩基质的行为,采用带有孔压的"cohesive"单元描述水力裂缝的性质,模型对"Texas Two-Step" 压裂方法进行了数值模拟,模拟结果得到了压裂过程中地层应力场的分布及其变化,模拟结果和解析公式计算结果吻合良好.模拟结果表明:(1)裂缝的产生减弱了地层应力场的各向异性;(2 对于低水平应力差页岩储层,采用"Texas Two-Step"压裂方法可以产生缝网. 对于采用"Texas Two-Step"压裂方法无法产生缝网的高应力差页岩储层,提出了三次应力"共振" 和四次应力"共振" 压裂方法并进行了数值模拟,模拟结果得到了压裂过程中页岩储层应力场的分布及其变化,得到了缝网形成的区域,模拟结果表明:(1)对于高应力差页岩储层,采用"Texas Two-Step" 压裂方法无法产生缝网;(2)对于高应力差页岩储层,三次应力"共振" 和四次应力"共振"压裂方法是有效的体积压裂缝网形成的方法.      

致密砂岩气藏纳微观结构及渗流特征

随着常规油气资源的减少,致密砂岩气藏逐渐成为勘探开发的热点.然而致密砂岩气高效产出机理研究还相对滞后,已经成为目前制约致密砂岩气大规模经济有效开发的瓶颈.致密砂岩由于其特殊的地质特征, 基质微观结构复杂,天然裂缝一定程度发育, 投产一般需要水力压裂等增产措施,气体在其中渗流存在跨尺度效应.研究揭示出致密砂岩气在多孔介质中的流动分为解吸、扩散、渗流等几个方式,包括浓度场下的扩散、压力场中的渗流等.综述了致密砂岩存在的非常规地质特征及其研究描述方法,并依据克努森数大小, 对致密砂岩气的渗流状态进行流态区域划分,为更好的理解和模拟致密砂岩气的流动状态提供了理论依据,并指出建立一种高效、简洁的微观结构精细描述方法,以及解决致密砂岩气在多尺度条件下的传质优化等是今后研究的重点方向.