3-单相液体渗流数学模型

1第三章单相液体渗流数学模型3单相液体渗流数学模型2“油气层渗流数学模型”就是用数学语言综合描述油气层渗流过程中全部力学现象和物理化学现象的内在联系和运动规律的方程式（或方程组）。一个完整的油气层渗流数学模型包括两部分：一部分是描述油气渗流的综合微分方程，另一部分是其相应的定解条件（包括初始条件和边界条件）。3单相液体渗流数学模型3第一节渗流数学模型的建立原则第二节渗流数学模型的微分方程第三节渗流数学模型的定解条件3单相液体渗流数学模型4第一节渗流数学模型的建立原则3单相液体渗流数学模型5一、建立数学模型的基础油气渗流力学的研究方法：是把一定地质条件下油气渗流的力学问题转换为数学问题，对数学问题利用合适的方法求解后，运用到油气田开发的实际生产中。将渗流过程中的各种力学、物理、化学现象和规律，用数学语言进行描述，实际上就是用微分方程或微分方程组对这些现象和规律加以表达。由于渗流的形态和类型不同，它们遵循的力学规律有差异，伴随渗流过程出现的物理化学现象也不同，因此油气渗流数学模型的类型很多。3单相液体渗流数学模型62.实验基础利用渗流物理基础实验认识力学现象和规律，是建立数学模型的关键。一、建立数学模型的基础1.地质基础孔隙结构：与数学模型相对应边界条件：几何形状、边界性质、参数分布、初始条件：原始状况3单相液体渗流数学模型7无穷小微元体上：分析力学现象，物理量之间内在联系，建立微分方程（数学模型）。数学模型建立后，用数学理论论证是否有解？连续？唯一？3.科学的数学方法•微分法3单相液体渗流数学模型8运动方程（必须）状态方程（弹性）质量守恒方程（连续性方程）（必须）能量守恒方程（非等温）附加方程（如：扩散方程）初始条件和边界条件（必须）二、渗流数学模型的结构渗流数学模型要综合反映渗流过程中，各种现象（力学、物理学、化学及相互作用）的内在联系。3单相液体渗流数学模型9数学模型解决的问题有四类：（1）压力分布，p=f(x,y,z,t)（2）渗流速度分布，v=f(x,y,z,t)（3）饱和度分布，S=f(x,y,z,t)（4）界面移动规律，分界面~t三、建立渗流数学模型的步骤1、确定建立模型的目的和要求明确数学模型解决什么问题。3单相液体渗流数学模型102、研究各物理量的条件和情况（1）过程状况：等温？非等温？（2）系统状况：油藏？气藏？单组分？多组分？（3）相态状况：单相？多相？混相？（4）流态状况：线性渗流？非线性渗流？非牛顿渗流？物理化学渗流？3单相液体渗流数学模型113、确定未知量和其它物理量之间关系（1）渗流速度与压力梯度：写出运动方程（2）物性参数与压力：写出状态方程（3）渗流速度或饱和度与时间：写出连续性方程（4）其它物理化学作用的函数关系：能量转换方程、扩散方程3单相液体渗流数学模型124、写出数学模型综合微分方程组连续性方程作为综合方程，其它方程代入连续性方程，得到描述渗流过程全部物理现象的统一微分方程。5、量纲分析量纲分析可以检验所建数学模型是否正确。检查所建数学模型量纲是否一致，是否是齐次的。3单相液体渗流数学模型136、确定数学模型的适定性数学模型建立后，用数学理论论证是否有解？连续？唯一？7、给出初始条件，边界条件3单相液体渗流数学模型14第二节渗流数学模型的微分方程3单相液体渗流数学模型15一、运动方程服从达西定律：zpKvypKvxpKvzyxgradpKpKv3单相液体渗流数学模型16121ndLdpCvn)(2vvKdLdp服从非达西定律：3单相液体渗流数学模型17二、状态方程描述由于弹性而引起力学性质随状态而变化的方程式称为状态方程。1、液体的状态方程dpdVVCLL1定义3单相液体渗流数学模型18Vp设：液体的质量为MLVMdMdVL2dpdC1dpdVVCLL13单相液体渗流数学模型19分离变量，C取为常数，积分：00ddpCppdpdC1)(ln00ppC0()0Cppe3单相液体渗流数学模型20将)(0ppCe展开成麦可劳林级数，只取前两项：00[1()]Cpp式中：p0—大气压力；ρ0—大气压力p0下流体的密度；ρ—任一压力p时流体的密度。)(00ppCe3单相液体渗流数学模型212、岩石的状态方程压力变化会引起孔隙大小发生变化孔隙大小变化引起渗透率的变化pip3单相液体渗流数学模型22∵dpdVVC1dVVVVdVdVbb)/()/(dpdC1定义3单相液体渗流数学模型2300ddpCpp弹性介质的状态方程描述了孔隙介质在符合弹性状态变化范围内，孔隙度的变化规律。当压力降低时，孔隙缩小，将孔隙中部分流体排挤出去，而成为驱动流体的弹性能量。分离变量积分：dpdC1)(00ppCe)](1[00ppC3单相液体渗流数学模型2414woMPa10~~nccccccwo式中：p0—大气压力；φ0—大气压力下岩石的孔隙度；φ—任一压力p时岩石的孔隙度。)](1[00ppC3单相液体渗流数学模型25质量守恒定律（连续性原理）：在地层中任取一个微小的单元体，在单元体内若没有源和汇存在，那么包含在单元体封闭表面内的流体质量变化应等于同一时间间隔内流体流入和流出质量之差。三、质量守恒(连续性)方程xyz21MMM3单相液体渗流数学模型26在渗流数学模型中，用它来描述渗流过程各种力学规律和物理化学规律之间的内在联系，通过把运动方程、状态方程和其它方程在质量守恒定律上联系起来，成为一个描述渗流过程全部力学过程的微分方程组。3单相液体渗流数学模型27连续性方程的表现形式：描述运动要素（速度、密度、饱和度、浓度）随时间和坐标的变化关系，在稳定渗流时则是描述这些要素和坐标之间的变化。3单相液体渗流数学模型28渗流过程常见的连续性方程单相流体渗流的连续性方程多相流体渗流的连续性方程带传质扩散过程的连续性方程建立连续性方程的方法微分法积分法3单相液体渗流数学模型291.用微分法建立连续性方程在地层中取微小六面体单元xyzM3单相液体渗流数学模型30xyzM在dt时间内，小单元体中流体质量的变化M等于在同一时间内流入小单元体中的流体质量M1与流出小单元体的流体质量M2之差。21MMM3单相液体渗流数学模型31•M点的流速：vx质量流速：vx•M点的质量流速：dydzdtdxxvvxx)2)((•dt时间经M点（ab）流入的质量流量：dxMMMxababdxdydzMMM2)(dxxvvxx3单相液体渗流数学模型32•M点的质量流速：dydzdtdxxvvxx)2)((dxMMMxabab•dt时间经M点（ab）流入的质量流量：2)(dxxvvxx3单相液体渗流数学模型33六面体在dt时间内x方向流入流出的质量流量差：dxdydzdtxvx)(dydzdtdxxvvxx)2)((dydzdtdxxvvxx)2)((3单相液体渗流数学模型34同理，在dt时间内y、z方向流入流出的质量流量差分别为：dydxdzdtyvy)(dzdxdydtzvz)(3单相液体渗流数学模型35在dt时间内六面体内流入与流出的总质量流量差为：dxdydzdtzvyvxvzyx)()()(3单相液体渗流数学模型36经过六面体流入与流出的质量不同，是因在六面体岩石和液体弹性能量的作用下，释放或储存一部分质量的结果（岩石的弹性表现为孔隙度的变化，液体的弹性表现为液体密度的变化）。六面体内的质量变化应为：六面体内的孔隙体积：dxdydz六面体内的流体质量：dxdydz单位时间内流体质量变化率：dxdydzt)(dt时间流体质量的总变化率：dxdydzdtt)(3单相液体渗流数学模型37dt时间内六面体总质量变化：等于在dt时间内流入与流出的质量差：tzvyvxvzyx)()()()(单相均质可压缩流体在弹性介质中渗流的质量守恒方程（连续性方程）dtdxdydztdxdydzdtzvyvxvzyx)()()()(流入流出质量差六面体内质量变化3单相液体渗流数学模型380)()(vdit散度()yxzvvvdivvxyztzvyvxvzyx)()()()(3单相液体渗流数学模型39如果不可压缩流体（ρ=常数）在刚性介质中流动（=常数），则连续性方程为：0)(vdi物理意义：六面体流出流入质量差为零，即流入六面体的质量与流出的质量相等，它仍然是一个质量守恒方程式——稳定渗流的连续性方程。0)()(vdit3单相液体渗流数学模型402.用微分法建立连续性方程（自学）3单相液体渗流数学模型41四、典型渗流数学模型的微分方程1、单相不可压缩液体稳定渗流微分方程假设单相液体在均质介质中的渗流为满足线性渗流规律的等温稳定渗流过程，不考虑多孔介质及流体的压缩性。运动方程gradpKv3单相液体渗流数学模型420)(vdiv连续性方程由上述两个方程得：0)(gradpKdi3单相液体渗流数学模型432222220pppxyz拉普拉斯方程0)()()(zpzypyxpx因K/μ=常数，则：02p单相不可压缩液体在均质地层中稳定渗流的数学模型3单相液体渗流数学模型4402222222zpypxpp2Laplace算子0)(p0v0)(vdiv等价于3单相液体渗流数学模型45使用条件：单相均质液体线性运动规律不考虑多孔介质及流体的压缩性稳定渗流渗流过程等温2222220pppxyz3单相液体渗流数学模型46坐标系三维问题一维问题直角坐标(x,y,z)2222222zpypxpp222xpp圆柱坐标(r,,z)2222221)(1zpprrprrrp)(12rprrrp球坐标(r,,)2222222211()sin1(sin)()sinpprrrrrppr)(1222rprrrp各种坐标系下的拉普拉斯算子的表达式3单相液体渗流数学模型471）运动方程2）状态方程)](1[00ppc)](1[00PPc2、弹性多孔介质中单相可压缩液体的不稳定渗流微分方程gradpKv3单相液体渗流数学模型483）质量守恒方程（连续性方程）0)(vdit0)(vt或左端第一项？左端第二项？3单相液体渗流数学模型49200000000)())((ppCCppCC由状态方程得：))((00000ppCC∵C、C很小，CC可忽略，则：左端第一项？3单相液体渗流数学模型50令Ct=C+C，综合压缩系数。物理意义：单位岩石孔隙体积在降低单位压力时，由于孔隙收缩和液体膨胀总共排挤出来的液体体积，可以看成常数。)(00000ppCt))((00000ppCC3单相液体渗流数学模型51将其代入质量守恒方程)(00000ppCttpCtt00)(质量守恒方程左端第一项0)(vdit3单相液体渗流数