页岩气直井非稳态非线性渗流的数值计算及产能预测

结合我国南方海相页岩气藏的实际,针对页岩气藏的特殊孔隙结构及多尺度流动的特征,首先由Knudsen数判断出储层主体流态为滑移流,并指出页岩气在页岩多孔介质中的传输过程是一个等温过程.然后修正了Beskok-Karniadakis二阶近似模型,建立了致密页岩气直井平面径向非稳态非线性渗流问题的数学模型.进一步应用玻尔兹曼变换将问题简化,并给出了求解压力场的数值离散显式迭代格式.最后数值计算得到了内边界定压条件下直井的压力场随时间、空间分布曲线,预测了产量随时间变化的规律,并分析了解吸、滑移和扩散效应对产量的影响.     

纳米尺度下页岩基质中的页岩气渗流及渗吸特征

基于页岩基质的纳米尺度孔隙特征,分析了页岩气在纳米孔隙中渗流的扩散、滑脱和达西渗流等对页岩气流动流量的影响.压差作用下页岩气流量的特征研究表明:达西流动产生的流量与地层压力、压力梯度和渗透率成正比;扩散引发的流量与压力梯度成正比,与平均压力成反比,与渗透率无关;滑脱引发的流量与压力梯度、渗透率成正比,压力梯度和渗透率一定时,滑脱流量为定值,与地层平均压力无关.低压下,气体的扩散效应和滑脱效应明显,扩散产生的流量所占比例较大,高压下渗流以达西流为主;在整个从低压到高压的地层压力区间,在低压下由于扩散作用,总的页岩气有一定的流量,随着地层压力的增大扩散产生的流量减小,到一定值后随着达西流量的进一步增大,总流量随地层压力的增加而增大,说明页岩气生产中存在一个最低效生产的压力区域.水在页岩薄片上吸水过程近似分为三个阶段:(1)水在页岩薄片表面的快速吸附;(2)水在裂缝中较快速度的渗吸;(3)水沿着微裂缝的缓慢渗吸.为页岩气开发中的压裂水的反排特征研究提供了理论基础.     

页岩气复杂孔裂隙真实气体传输机理和数学模型

基于滑脱流动和努森扩散,分别以分子之间碰撞频率和分子与壁面碰撞频率占总碰撞频率的比值作为滑脱流动和努森扩散的权重系数,进行权重叠加,建立了页岩气复杂孔裂隙气体传输模型.该模型综合考虑了滑脱效应和真实气体效应,同时还分别考虑了截面类型(圆形和矩形)和形状对气体传输的影响.用公开发表的分子模拟数据验证模型.结果表明:(1)本文模型能够合理地描述页岩气复杂孔裂隙气体传输机理,包括连续流动、滑脱流动和过渡流动;(2)页岩气孔裂隙截面类型和形状影响气体传输能力,相同截面面积,圆形截面孔裂隙气体传输能力大于矩形截面孔裂隙气体传输能力,矩形截面孔裂隙气体传输能力随纵横比增大而减小;与截面类型相比,截面形状对气体传输能力的影响更大;(3)真实气体效应提高了气体传输能力,且这种影响随压力增大而增大,随孔裂隙尺度减小而增大;(4)与圆形截面相比,真实气体效应对矩形截面气体传导率影响更大,且随矩形截面纵横比增大而增大.本文模型能为页岩气准确数值模拟奠定一些理论基础.     

耦合基质-天然裂缝渗流机理的页岩气表观渗透率

页岩储层微纳米孔隙、天然裂缝发育,不同类型孔隙中气体赋存状态和传输机理各异.本文在对天然裂缝定量表征的基础上,兼顾页岩气渗流特征和实际工程应用需要,综合考虑页岩岩芯孔-缝发育特征,基于渗透率串-并联模型建立耦合基质微纳米孔隙气体滑脱效应、扩散效应和天然裂缝渗流特征的页岩气表观渗透率模型,利用四川盆地牛蹄塘组天然裂缝发育的页岩岩样对理论模型进行验证.研究结果表明,该渗透率模型即能高度拟合天然页岩岩芯渗透率实验测量结果,又准确描述了岩芯内部基质孔隙和天然裂缝共同渗流的特点,更加符合实际渗流情况.总之,本文建立的页岩天然裂缝定量表征方法和综合考虑微纳米孔隙、天然裂缝渗流特征的表观渗透率模型不仅为页岩气藏复杂裂缝网络建模提供了一种新的手段,而且进一步推动了复杂缝网页岩气藏数值模拟研究工程应用的实现.     

页岩气储层微观渗流机理研究

页岩气储层孔隙结构复杂,纳米孔隙所占比例大,而纳米孔内部气体流动机理不同于宏观流体流动,因此认识页岩气在纳米孔隙中的流动机理对页岩气的高效开采具有重要的科学意义.页岩气开采过程中,纳米孔的吸附解吸、应力敏感效应及滑脱效应使渗透率发生显著变化.为此,基于毛细管模型耦合考虑吸附变形修正应力应变的渗透率模型,在此基础上,考虑...     

页岩气藏多重介质流体跨尺度传质表征方法研究

页岩气藏多尺度孔隙介质发育、非均质性强,不同尺度介质流体运移机制复杂,明确多重介质流体跨尺度传质表征方法的适用范围对开展页岩气藏数值模拟研究具有重要意义.本文基于体积压裂改造页岩气藏储层多尺度介质分布特征及多尺度介质气体运移特性,以储层不同区域微小单元体为研究对象,分别建立离散介质、拟稳态窜流双重介质、瞬态窜流双重介质表征单元体模型,对比分析不同储层条件和裂缝参数时离散介质模型与不同窜流模式双重介质表征单元体模型内流体产量变化规律,得到页岩气藏多重介质流体跨尺度传质表征方法适用范围.研究结果表明:生产早期不同窜流模式对窜流量影响较大,多重介质渗透率相差越大,生产早期不同窜流模式得到的窜流量相差越大,窜流量差异持续时间越长;采用瞬态窜流或拟稳态窜流双重介质模型可以较为准确地描述干酪根有机质-无机基岩之间流体跨尺度传质规律;无机基岩-次生缝网流体跨尺度运移规律需要采用离散介质模型进行描述.离散缝网耦合干酪根-无机基岩双重介质模型可以比较精确描述体积改造页岩气藏流体运移规律,该模型对开展页岩气藏数值模拟器研究具有重要的理论意义.     

页岩气非线性渗流的有限元模拟方法

基于考虑解吸、扩散和滑移作用的页岩气非线性渗流理论模型,构建了以压力为自变量的页岩气非线性渗流的三维有限元方程.通过一个忽略解吸效应的页岩气稳定渗流算例和一个考虑解吸附效应的页岩气非稳定渗流算例,求解井筒附近压力场,以及压力和气井产量随时间的变化,并与文献中的理论解进行对比,验证数值方法的合理性.结果表明,本文提供的页岩气非稳定渗流数值计算方法,其结果与理论解吻合,且收敛性好,可用于页岩气开发中气井周围压力变化和气井产量的预测,为页岩气井设计和页岩气开采提供重要的手段.     

页岩气藏数值模拟方法评价及选择

复杂的裂缝网络给页岩气藏数值模拟研究带来巨大挑战.本文利用离散裂缝模型(DFM),分别建立微尺度模型和矿场尺度模型,研究了不同基质渗透率和裂缝开度对应条件下,不同类型裂缝对页岩气总产量的贡献率.研究结果表明:在量级约为(10–8~10–4)×10–3μm2的页岩气藏基质渗透率范围内,且存在SRV改造区的条件下,页岩气产量以SRV内人工裂缝的贡献为主,SRV内以及SRV外天然裂缝的贡献量相对较小.因此,综合考虑计算复杂度和计算精度,可通过构建基质-天然裂缝表观渗透率和连续介质模型的方法表征天然裂缝,采用离散介质模型中计算网格少、效率高的离散裂缝网络(DFN)模型表征人工裂缝.     

页岩纳米孔吸附气表面扩散机理和数学模型

页岩富含纳米孔,且吸附气占总气量可高达85%,因此页岩气表面扩散对气体传输具有重要的作用.页岩气藏压力高,页岩表面能量非均质性强,吸附气非等温解吸附等,均加剧了吸附气表面扩散模拟的复杂性.基于低压条件下推导的Hwang模型,考虑高压条件下吸附气覆盖度的影响,建立了页岩吸附气表面扩散模型.同时,该模型还考虑了页岩表面能量非均质性、等量吸附热和非等温解吸附对表面扩散的影响.研究表明:1)表面扩散系数随压力的增大而增大,随温度的升高而增大,随表面活化能的减小而增大,随气体分子量的减小而增大;2)黏性流动、努森扩散和表面扩散对气体传输的贡献是此消彼长的,主要受孔隙尺度和压力的控制;3)表面扩散在微孔(半径2 nm)中,对气体传输贡献大,可高达92.95%;在宏孔(半径50 nm)中,贡献低于4.39%,可忽略;在介孔(2半径50 nm)中,表面扩散的贡献介于微孔和宏孔之间.     

考虑页岩气赋存及非线性流动机理的产能预测半解析方法

常规油气井可较容易求得解析产能,而页岩气的解吸、滑脱和扩散等非线性流动机理导致产能模型具有严重的非线性.本文运用三线性流模型,引入新的拟压力和拟时间处理页岩气的解吸和压缩性的非线性,运用逐次替换法处理滑脱和扩散的非线性.引入动用范围的概念,结合物质平衡方法和牛顿迭代法求得不同时间下的平均压力,以此不断更新模型参数,获得产能的半解析解.结果表明:(1)该方法能高效处理模型的非线性,准确快速地获得页岩气井的产能;(2)忽略气体的压缩性会严重低估气井产能;滑脱和扩散提高了气井产能,但加剧了产量递减;吸附解吸提高了中后期的产能,减缓了产量递减;(3)该方法能很好地拟合矿场生产数据,并预测气井的产能和可采储量.     

基于数字岩心和LBM的页岩气流固耦合数值模拟

页岩具有很强的压力敏感性,围压和孔压的变化会改变页岩孔隙的大小,从而对页岩气的流动规律产生影响,利用数字岩心结合格子Boltzmann方法(lattice Boltzmann method,LBM)来研究页岩气微观渗流规律得到越来越多学者的重视.本文建立了应力条件下的数字岩心应力应变模型和页岩气渗流LBM模型,研究了应力对页岩气渗流的影响规律.研究结果表明:有机质中的纳米孔隙对应力更加敏感,随应力变化的程度相比矿物骨架孔隙更大,从而影响页岩气在纳米孔隙中的解吸和扩散;孔压对各渗流机理的影响要比围压的影响大,是因为孔压的变化不仅影响了孔隙的尺寸还影响了气体的平均分子自由程;当平均孔压从17 MPa降低到5 MPa时,解吸的气体量和通过扩散流动的气体量占总气体流量的比例不断增加,分别增加了2%和1.9%,而通过滑脱流动的气体量占总气体流量的比例不断减少,减少了3.8%.利用应力条件下的数字岩心和格子Boltzmann方法可以更精确地模拟页岩气在储层中的流动规律,更好地理解页岩气的产出机理.     

致密砂岩气水相对渗透率模型

致密砂岩富含微米级孔隙,其气、水两相流动复杂,如何建立可靠的相对渗透率模型是目前亟待解决的问题.基于流体单管流动方程耦合二阶滑移模型,结合分形理论,建立了致密砂岩储层气、水相对渗透率模型,该模型考虑了气体滑脱效应、孔喉结构参数、含水饱和度分布,模型的可靠性通过多个已发表的气水相对渗透率实验数据验证.结果表明:1)建立的模型可用以描述致密砂岩储层微米尺度孔隙中的流体流动行为; 2)压力和含水饱和度增加都会减小滑脱效应对气相相对渗透率的提高作用,当压力小于1 MPa时,气相相对渗透率对压力极为敏感;温度对气相相对渗透的影响有限,可忽略; 3)孔径分布分形维数和孔隙迂曲度分形维数增加,有利于气相相对渗透率的提高,但会降低水相相对渗透率,且孔径分布分形维数的影响更大.     

气体吸附对页岩裂缝表观渗透率和页岩气采收率的影响

页岩气以及无水压裂技术采用的二氧化碳都可使页岩发生吸附变形,影响其裂缝表观渗透率,同时表观渗透率还受到有效应力和流体流态的影响.页岩气渗流过程中多种物理场相互耦合使得吸附应变对表观渗透率的影响不易被分析,常被忽略.本文基于多孔弹性理论和页岩吸附变形的特点,建立了适用于多种边界条件的页岩裂缝表观渗透率模型,提出了分析吸附变形对表观渗透率影响的方案,并分析了一组富有机质页岩样品中吸附应变对其表观渗透率的影响,最终结合页岩气多尺度渗流模型分析了吸附变形对页岩气采收率的影响.研究结果表明:采用不同边界条件和气体测量页岩表观渗透率时,其主控机理并不相同;基于此特点,结合"固定围压"和"固定孔隙压力"条件下的非吸附性气体和吸附性气体测量所得的页岩表观渗透率数据,可以确定内部吸附变形对表观渗透率有较大影响,不可忽略;此外,相较于内部吸附变形,气体吸附引起的页岩整体变形对表观渗透率和页岩气采收率有更大影响.     

页岩气多尺度传质特征及过程协调机制研究

页岩气开发过程中的多尺度传质行为显著,综合考虑页岩多尺度孔隙结构、多种气体传输机理和多种微尺度效应,构建了页岩表观等效渗透率数学模型,分析了实际生产过程中的页岩渗透性变化特征,明确了页岩气多尺度传质行为主控因素(气-固系统性质、孔隙类型及连通性、温度压力条件)的影响,提出了页岩气多尺度传质行为评价指标,即孔隙直径介于50~200 nm孔隙所占比例,从介科学的角度形成了页岩气多尺度传质行为评价方法,并结合页岩气开发所面临的主要问题,凝练了页岩气多尺度传质过程协调机制,即降低水相圈闭损害、沟通基质纳米孔及纳米孔扩容、促进解吸-扩散.     

页岩气纳米孔真实气体传输模型

页岩富含纳米孔,纳米孔气体传输不同于宏观流体流动.基于滑脱流动和克努森扩散两种传输机理,分别以分子之间碰撞频率和分子与壁面碰撞频率占总碰撞频率的比值作为滑脱流动和克努森扩散的权重系数,耦合这两种机理,建立了理想气体传输模型.同时考虑高压条件下真实气体分子间相互作用力和气体分子自身体积对气体传输的影响,建立了页岩纳米孔真实气体传输模型.模型可靠性通过分子模拟结果验证.结果表明:纳米孔真实气体传输模型能够更合理地描述所有的气体传输机理,包括连续流动、滑脱流动和过渡流动;真实气体效应对气体传输的影响可高达23%,其受压力、温度、纳米孔尺度和气体类型的控制;在室内实验条件下模拟页岩纳米孔气体传输时,用氦气代替甲烷,低估了甲烷的传输能力65.09%;用氮气代替甲烷,高估了甲烷的传输能力106.27%.     

一种表征煤储层压力与流体饱和度关系的数学模型

目前,针对生产过程中煤储层内部气、水两相渗流的研究尚薄弱,已建立的煤层气产能预测模型或产能分析方法中,一部分假设煤储层内部为单相渗流即单相水或单相气,一部分忽略流体饱和度和压力变化梯度,另一部分通过物质平衡方程建立平均地层压力与平均含水饱和度间关系,进而研究储层饱和度分布对气井产能的影响,但平均含水饱和度无法表征储层内部的饱和度分布.故现有模型均未考虑储层流体饱和度分布特征对气井产能的影响,导致现有产能模型预测值与实际煤层气井产出动态始终存在较大误差.鉴于此,根据气、水两相渗流微分方程,考虑应力敏感、基质收缩与气体解吸效应,建立了一种能够表征煤储层压力与流体饱和度间关系的数学模型.结合气、水两相相渗曲线与各生产时刻的储层压力分布曲线,提出迭代算法求解该模型,得到各生产时刻的流体饱和度分布规律.该模型的可靠性通过与商业数值模拟软件结果对比得到验证.该模型的提出填补了目前煤层气开发理论的空白,有助于准确计算气、水两相拟压力,为煤层气藏产能的准确预测奠定理论基础.     

疏松砂岩稠油油藏压裂井裂缝参数优新方法

世界上稠油油藏可采量超过了规原油.由于稠油油藏地层一般比较疏松生产易出砂,压裂防砂完井是疏松砂岩稠油油藏用造措施之一.而裂缝参数选取对于压裂防砂施工来说是至关重要.针对目前裂缝参数优化仅仅以产能为目标缺陷,本文引用裂缝流比重念来表征人工裂缝对油藏流体传输量贡献大小,裂缝流比重越大双线性流模式越标准,越有利于防砂生产.建立了压裂井数值模拟模,以计算不同宽比矩泄油区域中压裂井期生产情况.统一压裂设计方法基上,以增产和防砂为双重优化目标提出了一种优化疏松砂岩稠油油藏压裂井裂缝参数新方法：即对于定条件油藏和压裂用支撑剂数量,存一个特定无因次裂缝流能力使缝流比重数取最大值,此时裂缝参数即为最优裂缝参数.回归了正方泄油区域压裂防砂井参数优化公式,可以据油藏和支撑剂条件方便地进行压裂井裂缝参数优化.对于矩泄油区域,可以用等效支撑剂数消除宽比影响,将其等效为正方泄油区域之后进行裂缝参数优化.     

致密油气藏多段压裂水平井渗流模型研究

建立不同边界条件下的负表皮多段压裂水平井不稳定渗流模型,综合考虑压裂裂缝呈非等间距、非等半长、非等缝高、非等任意倾角的实际情况。利用源函数理论、三维特征值法和正交变换等现代数学分析方法对模型进行求解,获得三种外边界下的考虑井储和表皮影响的多段压裂水平井的井底压力动态响应解。最后,编程绘制标准的压力动态响应样版曲线,并对各参数进行敏感性分析。实例应用表明,该模型可以用作多段压裂水平井渗流特征研究及试井分析。     

深层油气藏水力压裂裂缝扩展数值模拟

深层油气藏水力压裂过程中的岩石变形,不仅要考虑弹性变形还要考虑塑性变形的问题,以往的裂缝扩展数值模拟都没有考虑塑性变形.文章基于Drucker-Prager屈服准则和相关联的流动法则建立了弹塑性水力裂缝扩展的数学模型.模型使用润滑理论描述裂缝内流体流动,耦合岩体变形,采用内聚力模型对裂缝扩展过程进行描述.基于有限元方法编制了数值求解器,通过与经典模型计算结果比较,验证了数学模型和算法的正确性.研究结果表明:相比于弹性水力裂缝,在流体注入量相同时,弹塑性水力裂缝长度短,裂缝宽度大,扩展时需要更高的流体压力.由于塑性应变的存在,缝内流体返排后,弹塑性水力裂缝不会完全闭合.随着岩石摩擦角的增大,弹塑性水力裂缝的宽度减小,所需流体压力增大.相对于理想弹塑性,硬化岩石介质中的水力裂缝宽度较小,所需流体压力较大.     

微纳通道内稠密气体流动和换热的Monte Carlo模拟

使用Enskog模拟Monte Carlo法(ESMC)对稠密气体在微纳通道内的流动和换热进行了模拟, 并与直接模拟Monte Carlo法(DSMC)和一致性Boltzmann算法(CBA)所得的结果进行了对比, 对微纳通道内稠密气体流动和换热特性进行了分析. 结果表明, 当气体密度较大时, 稠密气体效应对流动和换热的影响不可忽略, 这种效应使得通道内壁面阻力系数减小, 且壁面换热特性也与不考虑稠密气体效应时有所不同.