气体渗流机理
- 格式:doc
- 大小:101.50 KB
- 文档页数:4
气体在页岩储层中的渗流机理气体在页岩储层中的流动主要经历 3 个过程:(1) 吸附在页岩储层基质表面的天然气解吸附后形成自由气存储在基质孔隙中;(2) 基质孔隙中的自由气(包括游离态 、溶解态气体和解吸附后形成的气体)向低压区(裂缝网络系统)扩散 ;(3) 天然裂缝和压裂诱导裂缝中的自由气以渗流的方式流向井底 。
1. 解吸附作用页岩气藏与常规天然气藏最主要的区别是部分天然气以吸附状态存储于页岩基质中 。
气体在页岩储层基质颗粒表面上的吸附主要受温度 、压力 、吸附物(气体类型和性质)、吸附体(储层类型 、比表面积 、固体吸附能力)等的影响。
对于给定的页岩气藏 ,吸附体和吸附物性质保持不变 ,气藏内温度变化范围较小 ,气体吸附量是压力的函数 。
在钻井 、完井和开采过程中 ,孔隙压力下降 ,吸附在基质颗粒表面的气体开始解吸附 。
在平衡状态和特定温度条件下描述岩石表面气体吸附量的函数形式主要有 3 种 :Henry 线性等温吸附定律 、Freundlich 指数等温吸附定律和Langmuir 等温吸附定律 。
Henry 等温吸附定律 :E H V V p =Freundlich 等温吸附定律 :n E FV V p =Langmuir 等温吸附定律 2 扩散页岩储层中的扩散作用是指在浓度差的作用下 ,游离相天然气从高浓度区域向低浓度区域运动 ,即天然气由基质向裂缝系统进行扩散 ,当区域浓度平衡时 ,扩散现象停止 。
依据扩散过程可以分为拟稳态扩散和非稳态扩散。
式(5)给出了拟稳态扩散方程(Fick 第一定律) ,即单位时间内通过垂直于扩散方向的单位截面积的扩散物质流量(扩散通量)与该面积处的浓度梯度成正比 。
式(6)给出了非稳态扩散方程(Fick 第二定律) ,即扩散过程中扩散物质的浓度随时间变化 。
拟稳非稳页岩气藏开发过程中 ,基质内的天然气浓度随时间变化 ,非稳态扩散方程能够更准确地描述页岩气的扩散过程 。
致密气渗流机理研究综述致密气渗流是指天然气在煤层气、页岩气等微孔介质中的流动过程。
由于储层孔隙度很低、孔隙连通性差,导致气体无法通过常规的气水平分离和垂直驱替作用实现开采。
研究致密气渗流机理对于高效开发致密气资源具有重要意义。
从气体渗流的角度来看,致密气的渗流过程可以分为弥散渗流和煤层气吸附解吸两种机制。
弥散渗流是指气体在孔隙中依靠浓度差异进行的扩散过程,其主要影响因素为气体浓度梯度和孔隙结构特征。
煤层气吸附解吸是指气体分子在煤体孔隙中与煤质表面发生吸附和解吸作用,其主要影响因素为煤质特征和温度压力条件。
煤质特征是影响致密气渗流机理的重要因素之一。
煤质的孔隙结构特征直接影响气体在孔隙中的扩散速率。
煤质孔隙主要包括无序毛细孔、裂隙、聚并孔和溶蚀孔等,其中无序毛细孔是致密气渗流的主要通道。
煤质还会影响气体在煤质表面的吸附量和解吸速率。
温度压力条件也对致密气渗流机理有重要影响。
随着温度和压力的升高,气体在煤体孔隙中的吸附量减小,解吸速率增加,渗流能力提高。
适当的温度和压力条件还有助于改善煤质的吸附解吸性能,提高煤层气产量。
研究致密气渗流机理还需要考虑气体相态变化对渗流行为的影响。
在低压条件下,气体主要以气态存在,流动行为以扩散为主。
随着压力的增加,气体逐渐进入连续相态,流动行为更加接近于流体流动。
气体的压缩系数和黏度等物性参数发生变化,也会影响渗流行为。
数值模拟方法在研究致密气渗流机理中发挥了重要作用。
通过建立储层物理模型和数学模型,可以传统实验难以观测到的气体渗流行为进行模拟和预测。
常用的数值模拟方法包括等效介质模型、多孔介质模型和离散元模型等,其能够定量分析孔隙结构对渗流行为的影响,提高气田开发效率。
致密气渗流机理研究已经取得了一定的进展,但仍存在许多问题需要进一步研究。
未来的研究重点将更加注重储层物性特征、温度压力条件和气相态变化对致密气渗流机理的影响,加强数值模拟方法的应用,以实现对致密气开采的更加准确预测和高效开发。
致密气渗流机理研究综述
致密气渗流是指气体在石层孔隙结构非常均匀且孔隙半径非常小(通常小于0.1微米)
的石层中的渗流现象。
由于石层孔隙结构的特殊性质和气体分子之间的相互作用力的特殊
性质,致密气渗流具有一些独特的流体力学特性,包括非达西行为、非线性渗流特性和压
力-渗透率非线性关系等。
致密气渗流的机理研究是理解和优化致密气开发的关键。
在过去的几十年里,研究人
员通过实验室模拟、数值模拟和实际生产数据分析等方法,对致密气渗流机理进行了广泛
的研究。
研究人员通过实验室模拟石层孔隙结构,研究致密气在孔隙中的渗流行为。
一些研究
使用纳米孔隙膜进行压力滞回实验,研究致密气在纳米孔隙中的渗流特性。
实验结果表明,致密气渗流呈现非线性行为,即气体渗流速率不随压力梯度线性增加。
研究人员利用数值模拟方法对致密气渗流进行了研究。
数值模拟方法允许研究人员在
不同条件下模拟致密气渗流的行为,例如不同孔隙结构、不同渗透性和不同温度等条件下
的致密气渗流。
研究结果表明,致密气的渗流行为受到多种因素的影响,如岩石渗透率、
渗透率非线性关系、气体分子之间的相互作用力等。
研究人员通过分析实际生产数据,研究了致密气渗流的机理。
实际生产数据可以提供
有关致密气渗流行为的珍贵信息,例如渗透率和压裂效果。
通过对生产数据的分析,研究
人员可以评估不同开发方案的有效性,并优化致密气开发。
致密气渗流机理研究综述致密气渗流是一种重要的储层类型,具有储层容易堵塞、连通性差等特点,因此对致密气渗流机理的研究具有极为重要的现实意义。
目前研究表明致密气渗流在气体吸附、孔隙流动、渗流压力、气体分子间作用力等方面存在多种机理,下面将分别介绍。
第一,气体吸附机理。
由于致密储层的孔隙度小、孔隙尺寸狭小,因此致密储层表面存在大量的吸附点,可以吸附气体分子。
这些吸附点在孔内形成一层分子屏障,限制了渗流的通透性和扩散性。
同时,吸附点也会减缓气体在储层内的分布和传递速度,导致气体渗流速度缓慢、流量小,对于气藏的开发和生产是不利的。
第二,孔隙流动机理。
致密储层的孔隙大小与气体分子尺寸相当,因此孔隙径向渗流失去了统治地位,更倾向于表面流动和弥散扩散。
这种孔隙流动机理使得气体在致密储层的渗流受到了严重的阻力,使得致密储层开发、生产成本极高。
第三,渗流压力机理。
致密储层与传统储层不同,渗流通道系统堆积成为一个互相隔离的结构,使得气体流动时,所有孔隙都需要参与到渗流过程中来。
而在这个过程中,孔隙内部的压力梯度即是渗流压力,其细微变化将使得气体渗流速率发生不小的波动。
而随着渗流压力逐渐增大,渗流阻力也不断增大,使得气体从致密储层中流出变得困难。
第四,气体分子间作用力机理。
致密储层中的气体分子之间存在一种粘合作用力,是渗流速率下降的重要因素之一。
因此,要掌握气体的分子间作用力机理,才能深入研究致密储层中的渗流特征,为高效开采致密气藏提供更具有实际意义的理论依据。
总之,致密气渗流机理涉及到气体吸附、孔隙流动、渗流压力、气体分子间作用力等多个方面的问题,需要综合考虑才能得出准确的结论。
这项研究工作对于深化致密气藏的认识、发掘、评价和利用都有着极其重要的作用。
高温高压气储层渗流机理研究在能源领域的发展中,天然气作为一种清洁、高效的化石能源,受到越来越多的关注。
而在天然气的生产与开采过程中,高温高压气储层的渗流机理研究显得尤为重要。
本文将围绕高温高压气储层渗流机理展开论述,重点关注渗流模型、流体性质、渗流机制及影响因素等方面的研究内容。
首先,高温高压气储层的渗流模型是研究的重点之一。
渗流模型的建立是分析储层内气体流动行为的基础,由于高温高压储层的特殊性,需要考虑气体的非线性行为、非稳定性以及非均匀性等因素。
目前常用的渗流模型主要包括达西定律、布尔斯定律以及Forchheimer扩展模型等,这些模型可以帮助我们预测储层气体的渗流规律。
其次,高温高压储层中的流体性质对渗流机理也有着重要的影响。
高压下,气体的多种物理性质会发生变化,如气体的压缩性、黏度以及热物性等,这些性质的变化将会对气体在储层中的流动行为产生显著影响。
因此,对高温高压储层气体的性质进行准确的测量和模型拟合是研究渗流机理的重要前提。
渗流机制是高温高压气储层研究的核心内容之一。
在高温高压条件下,渗流机制主要包括机械渗流、重力驱替、压差驱替、吸附解吸以及渗流剪切等。
这些机制的相互作用将决定气体在储层中的流动方式和产能特性。
研究这些机制对于深入理解高温高压气储层储集层的特性以及气藏的开采方式具有重要意义。
最后,影响高温高压气储层渗流机理的因素也是研究的重点之一。
温度、压力、含水饱和度、岩石裂缝度等因素都对气体在储层中的渗流行为产生着重要影响。
其中,温度和压力是最主要的因素,高温会导致气体的物性变化,增加流动阻力;高压则会增加气体的压缩度,改变渗流规律。
因此,研究不同因素对渗流机理的影响,可以为天然气勘探与开采提供科学依据。
综上所述,高温高压气储层渗流机理的研究对于理解储层内气体的流动行为、优化开采方案以及提高天然气产量具有重要意义。
通过建立准确的渗流模型、研究气体流体性质、深入探究渗流机制,同时考虑影响因素的作用,可以为高温高压气储层的开发与应用提供科学依据。
致密气渗流机理研究综述【摘要】摘要:本篇文章主要介绍了致密气渗流机理的研究现状和进展。
首先从基本概念与特点入手,介绍了致密气藏的特点和其渗流机理。
接着介绍了致密气渗流机理研究的方法与技术,以及研究的主要进展与成果。
进行了国内外致密气渗流机理研究的比较分析,指出了存在的问题与挑战。
探讨了致密气渗流机理研究的重要性和发展前景,并提出了未来研究方向和建议。
致密气渗流机理的研究不仅对油气开发具有重要意义,还有着广阔的应用前景。
【关键词】致密气渗流,机理研究,综述,基本概念,方法与技术,进展与成果,比较分析,问题与挑战,重要性,发展前景,未来研究方向1. 引言1.1 致密气渗流机理研究综述的背景介绍致密气渗流机理研究的背景可追溯至20世纪80年代末,当时我国首次在下扬子地区发现了具有较高产气潜力的致密气储层。
随后,国际上也陆续发现了类似的致密气资源,引起了人们对致密气储层的广泛关注。
致密气渗流机理研究在相关领域中具有重要意义,不仅可以帮助提高致密气藏的开发效率,还可以为相关工程应用提供理论支持。
在此背景下,对致密气渗流机理进行深入综述和研究具有重要的理论和实践意义。
通过对致密气渗流机理的系统分析,可以更好地理解致密气储层的特性和行为规律,为未来的致密气资源开发与利用提供技术支持和理论指导。
1.2 研究意义及目的研究致密气渗流的机理具有重要的理论和应用价值。
对致密气渗流机理的深入研究能够揭示气体在致密储层中的迁移规律,为气体的有效开采提供科学依据。
致密气渗流的研究对于了解油气运移规律、油气藏储层特性有着重要意义。
研究致密气渗流机理还可以为气体藏地质储层工程的优化设计和实际生产提供重要参考。
本文旨在系统总结致密气渗流机理研究的现状和进展,为深入探讨这一领域提供参考。
希望通过本文的研究能够对致密气渗流机理的研究方向和发展提出一些建设性的建议,推动该领域的进一步发展和应用。
2. 正文2.1 致密气渗流机理的基本概念与特点致密气渗流是指气体在致密储层中的渗流过程,是一种复杂的非常规气藏开采技术。
页岩气藏气体流动机理及数值模拟研究页岩气是一种以页岩为主要储层的天然气资源,由于其在储层中的特殊性质,其流动机理和数值模拟研究对于有效开发和利用页岩气具有重要意义。
在页岩气藏中,气体流动的机理主要包括渗流机理和吸附机理。
渗流机理是指气体在页岩储层中的渗流过程,主要受到渗透率、孔隙度和渗透率分布等因素的影响。
吸附机理是指气体在页岩储层中与页岩表面发生吸附作用,主要受到吸附等温线和吸附解吸速率等因素的影响。
为了研究页岩气藏中气体的流动机理,数值模拟成为一种重要的研究手段。
数值模拟可以通过建立数学模型和计算方法,模拟气体在页岩储层中的流动过程,对气体的渗流和吸附行为进行定量描述。
数值模拟可以通过改变渗透率、孔隙度和吸附等温线等参数,研究它们对气体流动的影响,从而为页岩气藏的开发和利用提供科学依据。
在数值模拟研究中,常用的方法包括有限差分法、有限元法和边界元法等。
这些方法可以通过离散化储层模型,将连续的流动方程转化为离散的代数方程,然后通过迭代求解,得到气体在储层中的流动状态。
数值模拟可以通过改变模型的边界条件和参数,模拟不同的开发方案和条件,评估其对气体产量和开发效果的影响。
然而,数值模拟研究也存在一些挑战和限制。
首先,页岩气藏储层复杂多变,储层参数的确定和模型的建立存在一定的不确定性。
其次,数值模拟需要大量的计算资源和时间,对计算机性能和算法效率提出了较高的要求。
此外,数值模拟结果的可靠性和准确性也需要通过与实际生产数据和实验结果进行验证。
尽管存在一些挑战,但数值模拟研究对于页岩气藏的开发和利用具有重要意义。
通过数值模拟,可以评估不同的开发方案和条件对气体产量和开发效果的影响,优化开发策略,降低开发成本。
此外,数值模拟还可以预测页岩气藏的产量潜力和剩余资源,为储量评价和资源管理提供科学依据。
页岩气藏气体流动机理和数值模拟研究对于有效开发和利用页岩气具有重要意义。
通过研究气体在页岩储层中的渗流和吸附行为,可以揭示气体流动的机理,为开发策略的制定和优化提供依据。
致密气渗流机理研究综述摘要:致密气储层是指气体渗流能力较差且渗透率极低的储层。
致密气渗流机理的研究对于油气开发和资源利用具有重要意义。
本文综述了致密气渗流机理的研究现状和进展,并对未来的研究方向进行了展望。
关键词:致密气;渗流机理;研究进展;研究方向1. 引言致密气储层是指渗透率极低且渗流能力较差的气体储层。
由于储层孔隙度低、孔隙连通性差以及岩石多孔介质的复杂性,致密气储层的渗流机理与常规气田有很大的区别。
致密气储层的渗流机理研究可以提供油气开发和资源利用的理论依据,对于合理开发致密气储层具有重要意义。
2. 渗流机理研究方法致密气渗流机理的研究可以通过实验室实验、数值模拟和现场观测等方法进行。
实验室实验可以模拟储层条件,通过测量渗流参数和分析气体渗流规律来研究致密气渗流机理。
数值模拟可以通过建立合适的数值模型,模拟渗流过程和预测渗流行为。
现场观测可以通过地震勘探、压力监测和产能测试等方法,直接观测和分析储层中的渗流现象和规律。
3. 渗流机理的主要影响因素致密气渗流机理受到多种因素的影响,包括储层岩石性质、裂缝和节理、孔隙连通性、孔隙流体性质和应力状态等。
这些因素之间相互作用,共同决定了储层的渗流能力和渗透率。
裂缝和节理对渗流能力的影响最为显著,裂缝和节理的存在可以提高储层的渗透率和渗流能力。
4. 渗流机理的研究进展致密气渗流机理研究已经取得了一些重要的进展。
在储层岩石性质方面,通过岩石力学性质的研究,揭示了岩石裂缝和节理对渗流的影响机理。
通过实验室实验和数值模拟,研究了孔隙连通性和孔隙流体性质对渗流的影响。
通过现场观测,研究了储层中裂缝和节理的分布和演化。
这些研究进展为致密气渗流机理的研究提供了理论基础和实验依据。
5. 研究方向展望致密气渗流机理的研究还存在一些尚未解决的问题。
致密气渗流机理的多尺度和非线性特性需要进一步研究和分析。
致密气渗流机理的化学反应和物理现象对渗流行为的影响也需要深入研究。
致密气藏水平井压裂气体示踪剂渗流机理石油钻采工艺致密气藏是一种储量巨大但渗透性很差的天然气储层,开采致密气藏通常需要采用水平井压裂技术。
而为了了解水平井压裂后气体在储层中的渗流机理,研究人员使用气体示踪剂进行实验和模拟。
本文将探讨致密气藏水平井压裂气体示踪剂渗流机理的石油钻采工艺。
在致密气藏开采过程中,水平井压裂技术被广泛应用。
压裂是指通过注入高压液体或气体来破裂储层岩石,从而增加储层的渗透性,提高气体的产量。
而在压裂过程中,气体示踪剂的使用可以帮助研究人员了解气体在储层中的流动路径、渗流速度以及储层的渗透性。
气体示踪剂是一种添加在压裂液中的可追踪的气体。
常用的气体示踪剂包括氦气、氮气和二氧化碳等。
在压裂过程中,这些气体示踪剂会随着压裂液一同注入储层中,并通过地层压力的驱动下,沿着裂缝或渗流通道向储层中扩散。
研究人员可以通过在不同位置设置传感器或监测装置,来实时监测气体示踪剂的运动轨迹和浓度变化,从而获得关于储层渗流机理的信息。
通过对气体示踪剂的监测,研究人员可以获得压裂液的扩散速度、扩散距离以及渗透性的变化情况。
例如,当气体示踪剂的浓度变化不明显或变化较慢时,可能意味着压裂液在储层中的扩散速度较慢,渗透性较低。
而当气体示踪剂的浓度变化快速而明显时,可能意味着储层存在高渗透性的通道或裂缝。
这些信息对于评估压裂效果、优化压裂设计以及预测气体产量具有重要意义。
研究人员还可以利用气体示踪剂的数据来构建数学模型,模拟储层中气体的渗流行为。
通过对气体示踪剂的渗流机理进行建模,可以更好地理解储层中气体的分布、运动路径以及渗流速度。
这些模拟结果可以用于指导油气田的开发和管理,提高开采效率。
在石油钻采工艺中,致密气藏的水平井压裂是一项复杂的技术。
通过使用气体示踪剂进行实验和模拟,可以深入了解储层中气体的渗流机理。
这些信息对于压裂设计、油气田的开发和管理以及气体产量的预测都具有重要意义。
因此,石油钻采工艺中的致密气藏水平井压裂气体示踪剂渗流机理的研究是一个重要的领域,值得进一步深入研究和应用。
页岩气渗流机理
页岩气是指那些聚集在暗色泥页岩或高碳泥页岩中,以吸附或游离状态为主要存在方式的天然气。
它与常规天然气的理化性质完全一样,只不过赋存于渗透率、孔隙度极低的泥页岩之中,气流的阻力比常规天然气大,很大程度上增加了页岩气的开采难度,因此被业界归为非常规油气资源。
页岩自身的有效孔隙度很低,页岩气藏主要是由于大范围发育的区域性裂缝,或热裂解生气阶段产生异常高压在沿应力集中面、岩性接触过渡面或脆性薄弱面产生的裂缝提供成藏所需的最低限度的储集孔隙度和渗透率。
通常孔隙度最高仅为4% ~5%,渗透率小于1×10-3μm 。
页岩气藏有特殊的产气机制。
与常规低渗气藏不同,天然气在页岩中的流动主要有4种机理,这4种机理覆盖了从分子尺度到宏观尺度的流动。
主要表现为游离气渗流、解吸附、扩散和自吸。
第一 ,由于气体滑脱效应的存在 ,游离气在有机质和无机质基岩中的流动属非达西渗流,但在天然或水力裂缝中的流动为达西渗流。
第二,有机质上的吸附气对渗透率有不利的影响,这是由于有机质的天然气吸附层对天然气分子的引力增大所致,但是,如果有机质不属于多孔介质,仅作为连接基质孔隙或为裂缝之用,那么,在生产时,远离孔隙和裂缝的吸附气只能沿有机质表面易扩散的方式进行运移。
如果有机质属于多孔介质,部分吸附气能够直接释放进入有机质孔隙,并且,这样会使扩散的重要性被减弱 。
第三,自吸作用是当压裂水在致密气藏流动时发生的一种现象,在页岩储层压裂时,由于自吸作用和重力分异作用,导致压裂水的返排率不足50% 。
因此,气水两相在裂缝中共同流动时,往往气在裂缝的上部流动,此时,在裂缝的下部留有大量的水。
在钻井液和增产措施作业水的冷却作用下,储层接触面附近会聚集更多的束缚水,因而也会恶化自吸现象的影响。
1 Langmuir 单分子层吸附状态方程
假定固体表面是均匀的,对气体分子只做单分子层吸附.设气体的压力为p,未被气体分子吸附的表面积百分数为
θ.气体分子吸附的速度与气体的压力成正比,
也与未被气体分子吸附的表面积成正比,则吸附速度
a R cp θ=
式中,c 为比例系数.
气体脱附的速度与吸附气体分子所覆盖的表面积的百分数成正比,也与被吸附的气体分子中那些具备脱离表面逸向空间所需能量的分子所占的比例成正比.设吸附气体分子所覆盖的表面积的百分数为θ,设εa 为脱离表面逸向空间所需的最低能量,即吸附热εa,被吸附在表面的总分子数为Na,其中能量超过εa 的分子数为N*a,则有
/
*/a k T
a
a N
N fe
ε=
式中,f 为比例系数;k 为玻尔兹曼常数.
则脱附速度
/
a k T
d R d
e ε
θ=
式中,d 为比例系数
达到吸附平衡时,吸附速度应等于脱附速度,即Ra= Rd,所以
/0a kT
cp d e
εθθ=
未被气体分子吸附的表面积百分数θ0与吸附气体分子所覆盖的表面积的百分数θ之和应等于1,即
01
θθ+=. 可得单分子层吸附方程
1bp
bp θ=
+ 式中,
/a kT
c b
d e
εθ=
如果以Q 表示单位固体表面上吸附的气体的量,a 表示单位固体表面上饱和吸附气体的量,则Langmuir 方程转化为常用的形式:
1abp
Q bp =
+
在压力很低时,上式分母中的bp 相对于1可以忽略不计,吸附气体量Q 与压力p 成正比;在压力很高时, 上式分母中的1相对于bp 可忽略不计,吸附气体量Q 达到饱和,即发生饱和吸附。
若给定吸附剂、吸附质和平衡温度后,则吸附量Q 就只是吸附质的平衡压力P 的函数。
2 页岩气渗流的连续方程
在地层中取微小六面体单元,单元体中心点M 的质量流速为v ρ,在各坐标方向的分量分别为x v ρ,y v ρ和z v ρ。
在d t 时间内六面体内流入流出的气体总量差为:
[
]y x z v v v dxdydzdt x
y
z
ρρρ∂∂∂-+
+
∂∂∂
经过六面体流入与流出的质量之所以不同,是因为在六面体内岩石和气体弹性能的作用下释放或储存一部分质量的结果。
六面体内的气体质量变化为:
游离气体质量变化:
()
()
gf gf gm gm dxdydzdt dxdydzdt t
t
ρφρφ∂∂+
∂∂
吸附气体质量变化:()
gsc Q dxdydzdt t
ρ∂∂
标记:下标gf ,gm ,gsc 分别表示裂缝气体,基岩气体,标准气体。
根据质量守恒定理,d t 时间内六面体总的质量变化应等于同一时间内流入
流出的质量差,即:
()
()
()
[]y gf gf gm gm gsc x z v Q v v dxdydzdt dxdydzdt dxdydzdt dxdydzdt
x
y
z
t
t
t
ρρφρφρρρ∂∂∂∂∂∂-++=
+
+
∂∂∂∂∂∂ ()
()()
[
]y gf gf gm gm
gsc x z v Q v v x
y
z
t
t t
ρρφρφρρρ∂∂∂∂∂∂-+
+
=
+
+
∂∂∂∂∂∂ ()()()()0gf gf gm gm gsc Q div v t t t
ρφρφρρ∂∂∂+++=∂∂∂
3 气体不稳定渗流的微分方程 对于页岩气不稳定渗流,若假定 (1) 渗流为等温过程 (2) 气体粘度为常数
(3) 地层均质即k ,φ为常数,且不随压力变化 (4) 压力梯度很小,服从线性渗流定律 气体满足以下方程式:
1. 运动方程:k v gradP u =
2. 状态方程:P R T z
ρ=
3. 气体等温压缩系数: 111()||T
T
dV dz C P V dP
P
z dP
=-
=
-
那么
111()|[()]11[]()T P P d
P
p
t t R Tz
R T
z t
dP z t
P
dz P P C P R Tz P z dP t t
ρ
ρ∂∂∂∂==
+=∂∂∂∂∂∂-=∂∂
将上式带入到连续性方程()()()
()0gf gf gm gm gsc Q div v t t t
ρφρφρρ∂∂∂+
++=∂∂∂ 中可得:
()()()()01gf gf gf gm gsc
gf gf gf
gm gm
gsc
gf
gsc k P P abP P C P C P u t
t
t bP ρφρφρρ∂∂∂
∇∇+++=∂∂∂+由于
ρ是压力的函数,因而上式是一个非线性的微分方程,解决这类问题的方法之一是先
将因变量P 换成压力函数:
令 P P P
dP C ρρ=
+∇=∇⎰
或,将 P 代入上式,则得到真实气体不稳定等温渗流的数学模型(综合微分方程)为: 2
2
()()
(1)g f
g f g m
g s c
gf gf
gm gsc
gf
p p k P ab P C P C P u t
t
bP t
φφρ∂∂∂∇+++
=
∂∂+∂ 2
2
()()
0(1)
gsc
gf gf gm gf gf
gm gf
p p k ab P P C P C P u t
t
bP t
φφ∂∂∂∇+++
=∂∂+∂。