油气层渗流力学

油气层渗流力学编辑ppt第五章：单相液体的稳定渗流理论5.1单相液体稳定渗流模型的典型解及其应用单相稳定渗流的数学模型为：本章主要研究三种典型的几何模型条带模型圆形模型球形模型5.1编辑ppt将5.1式化成极坐标形式为：

n=0:平面平行流n=1：平面径向流n=2：球面径向流

编辑ppt一单相液体平面平行渗流1假设条件：水平均质等厚（长度L，宽度W，厚度h2微分方程及其定解条件（n=1,r=x)5.3Pe pB编辑ppt3方程的求解--采用分离变量法求解LPePB方程的解压力分布编辑ppt3压力梯度分布4渗流速度分布编辑ppt5单相稳定流的产量公式编辑ppt6单相流的水动力学场图1)几个概水动力学场图----流线与等压线组成的图形等压面-----在渗流场中压力相等的空间点组成的面等压线----在同一渗流平面内,压力相等的点组成的线流线-----与等压线垂直的线编辑ppt2)等压线,流线的绘制规定---在作等压线时,必须使任意两条相邻等压线间的压力相等;作流线时,必须使任意两条相邻流线间的流量相等.3)绘制流线,等压线的意义----形象地描绘流体的流向及能量损耗规律,流速分布规律,比计算公式更直观和更生动具体编辑ppt4)单相单度流的水动力学场图等压线方程为:流线的方程为:水动力学场图----等间隔的水平线与垂直线形成的网格X=常数Y=常数编辑ppt7)应用确定产量确定生产压差编辑ppt二单相不可压缩流体平面径向流1假设条件1)圆形----2)达西流动,---编辑ppt2渗流微分方程3微分方程的求解编辑ppt4压力分布编辑ppt5压力梯度分布rrw编辑ppt6渗流速度分布编辑ppt7井产量计算编辑ppt8水动力学场图等压线方程-----r=常数流线方程------Q=常数水动力学场图编辑ppt9应用研究影响产量的因素研究提高产量的方法-----编辑ppt三、平均地层压力及质点运动规律1、问题的提出2、计算平均地层压力的方法面积加权法3、径向流平均地层压力的计算方法编辑ppt4、近似计算5、实例计算详见书p43编辑ppt四、质点移动规律1、问题的提出2、单向流质点移动规律问题的描述计算公式应用实例编辑ppt3、平面径向流质点移动规律问题的描述计算应用编辑ppt五、单相不可压缩的球面径向流1、假设条件（1)供给区是以re为半径的球面，。。。。（2)不可压缩流体按达西定律稳定渗流。2、渗流微分方程3、方程的求解方法----分离变量方法编辑ppt4、压力分布编辑ppt5、球形流压力梯度分布6、渗流速度分布编辑ppt7、产量公式8、应用条件厚油层、钻开程度不完善等编辑ppt5.2、气体稳定渗流微分方程的典型解气体渗流微分方程编辑ppt一、服从达西定律的气体单向稳定渗流1、假设条件单向流动等温的达西渗流气体保持稳定渗流状态2、方程及定解条件编辑ppt3、压力分布特点：压力平方与距离成线性关系4、压力梯度分布编辑ppt5、速度分布6、水动力学场图编辑ppt产量公式编辑ppt二、服从达西渗流定律的气体平面径向稳定渗流1、假设条件2、渗流方程编辑ppt3、压力分布（设uz=常数）编辑ppt5、压力梯度分布6、渗流速度分布编辑ppt7、水动力学场图8、产量公式编辑ppt第一节弹性不稳定渗流的物理过程第七章弹性液体在弹性多孔介质中的渗流第二节弹性液体不稳定渗流的微分方程第三节弹性不稳定渗流数学模型的典型解第四节单位制与法定计算单位第五节弹性不稳定渗流时油井间的干扰第六节弹性不稳定渗流压力波的传播规律编辑ppt第七章弹性液体在弹性多孔介质中的渗流第一节弹性不稳定渗流的物理过程一、力在刚性及弹性体中的传播编辑ppt1、关闭油井时，井底压力不是一下子恢复到静止地层压力Pe，而是逐步上升至Pe的。tPPePw2、甲井投产引起地层压力下降，乙井不会立刻感受到，而是要经过一个时阶段t之后，乙井井底压力及产量才受到影响而发生变化。3、开井时，井底压力不是一下子降至某个值就稳定下来，而是逐渐下降到某个值才稳定下来。编辑ppt二、弹性驱动时，流动的基本特征井筒QPe原始地层压力开井生产井筒压力下降井壁附近地层压力下降地层压力下降范围向外扩大压力波的传播非瞬时完成编辑ppt1、由于流体及岩层具有弹性，因而井底压力变化在地层中的传播是非瞬时完成的。2、地层中任意一点上的压力、流速将随时间而变。 P=f1(t,r); V=f2(t,r)3、地层压力的下降，将引起流体和岩石颗粒膨胀，孔隙度下降，从而迫使一部分流体在弹性力的作用下，从孔隙中排除流入井底。 =f1(P); =f2(P)弹性驱动的特点：编辑ppt三、井以定产量投产时地层压力变化规律ABEHR(t1)WUGR(t2)W1U1H1R(t3)（一）压力传播的过程编辑pptR(t1)R(t2)R(t3)地层中任意一点上的流速的变化：XX点处压降漏斗的斜率越来越大，表明该点的流速在逐渐增加，通过此点向井流动的液体越来越多。编辑ppt（二）外边界对压降漏斗的影响1、地层为无限大时压降漏斗不断扩大与加深的过程将是无限制地持续下去。编辑ppt2、地层边界有供给时ReGst=tB设t=tB时刻压力波到达供给边界。油井产量由两部分液体组成：Q1——通过供给边界流入地层的液量；Q2——供给半径Re以内，地层及流体弹性膨胀排除的液量。编辑pptReGst=tBt=编辑ppt3、地层边界无供给时Gst=tB设t=tB时刻压力波到达供给边界。在供给边缘上压力梯度为零：编辑ppt（三）常用的名词激动区域（影响区）——井底压力变化的波及区域。条件影响边缘半径（影响半径）——压力降波及区域的半径，即激动区半径。压降漏斗前缘——激动区前缘。弹性驱动第一相（不稳定流动初期）——压降漏斗前缘到达地层边界以前称为弹性驱动第一相。编辑ppt四、井以定井底压力投产时，地层压力及产量变化规律1、地层边界有液源供给的情况（1）地层压力变化分析GsRet=tBPw编辑ppt（2）井产量变化规律GsRePw假设在t=0时，井底压力瞬间从Pe降到Pw，地层中不发生压力降。Q实际生产：编辑ppt2、地层边界无供给的情况（1）地层压力变化分析GsRet=tBPwt=编辑ppt（2）井产量变化规律Qt编辑ppt第二节弹性液体不稳定渗流的微分方程可压缩液体在弹性多孔介质中按达西定律渗流的微分方程：编辑ppt在极坐标系中：式中：n=0 单向流 一维n=1 径向流 二维n=2 球形流 三维编辑ppt第三节弹性不稳定渗流数学模型的典型解一、无限大地层，井以定产量投产的典型解渗流微分方程：边界条件：初始条件：编辑ppt查数学手册知，幂积分函数所以无限大地层以定产量Q投产后，地层中任一点M上压降计算公式。编辑ppt幂积分函数可展开成级数：当时可以只保留级数的前两项。编辑ppt井底压降：考虑井壁污染或异常产生的附加阻力编辑ppt二、圆形封闭地层中心，一定产量井投产的典型解编辑pptn是方程J1(R)Y1()-Y1(R)J1()=0的根。Jo

——第一类零阶贝塞尔函数；Yo——第二类零阶贝塞尔函数；J1

——第一类一阶贝塞尔函数；Y1——第二类一阶贝塞尔函数；编辑ppt圆形封闭地层不稳定渗流典型解P(r,t)无因次化进行拉普拉斯变化常微分方程编辑ppt弹性驱动不稳定渗流第二相初期的实用公式。若进入弹性驱动第二相晚期，可简化为：编辑ppt视稳定流期地层任意一点上的压力计算公式：视稳定流时期，地层中各点压力降落速度相等，且等于一个常数。编辑ppt三、弹性不稳定渗流有界定压边界的典型解编辑pptn是方程J0(R)Y1()-Y0(R)J1()=0的根。Jo

——第一类零阶贝塞尔函数；Yo——第二类零阶贝塞尔函数；J1——第一类一阶贝塞尔函数；Y1——第二类一阶贝塞尔函数；编辑ppt编辑ppt四、无限大地层井以变产量投产时，不稳定渗流的典型解其中：——时间变量；Q()——随生产时间而变的产量；t——生产时间。编辑ppt假设=f(x,y,z,,t)代表D域内点M(x,y,z)处当时间为t的势，并满足热传导方程：初始条件：t=0=0边界条件：=g(x,y,z,t)编辑ppt无限大地层中，当井以变产量投产后，时刻t在离井r处所形成的压力降。井以定产量生产时，Q()=Q=常数；引入新变量：编辑ppt当=0时，当=t时,编辑ppt编辑ppt第四节单位制与法定计量单位1987年原石油工业部颁发了《油气藏工程常用参数代号及计算单位》，SY6255-87编辑ppt第五节弹性不稳定渗流时油井间的干扰设地层中有n口井同时在弹性驱动下投产，地层中任意一点M上的压力降，应等于每口井单独投产时，在该点形成的压力降的迭加。Pi-Pm——n口井同时投产后，时刻t在点M形成的压力降；Qj——j井的产量；rj——点M至j井的距离；tj——到时刻t为止，j井的生产时间。编辑ppt任意一口井的井底压力降：其中：K=1，2，3，，nr11=rw1，r22=rw2，，rkk=rwkrjk——j井至k井的距离；rwk——k井井半径编辑ppt第六节弹性不稳定渗流压力波的传播规律ReGst=tBt=为使研究方便，把流体渗流的不稳定渗流过程视为若干稳定状态的依序替换。这种研究方法称为稳定状态依次替换法。编辑ppt按照稳定状态依次替换法，任意任意时刻t的流动视为稳定的，压力分布可表示为：井产量Q可表示为：编辑ppt用物质平衡法求解：假设到时刻t为止，原油的总采出量为G，地层条件下流体的平均重率为ordr编辑pptrdr投产t时刻后，因弹性膨胀从小单元体中排除的液体重量为：式中：t=0时，=()e对应压力Pe

t时，=()对应压力P单元体体积流体密度及岩石孔隙的变化投产t刻后，地层的总采出量G1为：编辑ppt由状态方程可得：编辑ppt编辑ppt编辑ppt原油的总采出重量G应等于地层产液量G1与井筒产液量G2

之和编辑ppt编辑ppt导压系数的值反映了压力波传播的速度。编辑ppt油井井底压力的变化规律：编辑ppt第七节气体不稳定渗流微分方程的典型解气体不稳定渗流方程：考虑气体不稳定渗流时的导压系数为常数。气体渗流微分方程：编辑ppt一、无限大地层中一口井以定产量投产编辑ppt用压力平方表示为：编辑ppt二、圆形封闭地层中心一口井以定产量投产编辑ppt液体气体编辑ppt如果气体渗流达到拟稳定期之后，生产时间t已经很大：编辑ppt三、圆形供给边缘地层中心一口井以定产量投产编辑ppt液体：气体：编辑ppt当投产时间趋于无穷大时：当生产时间t>re2/4时，求和项可以略去不计。编辑ppt第八章油水两相渗流理论问题的提出前几章的假设条件：均质流体不考虑油和水在粘度和重度上的差别不考虑毛管力的影响地层压力必须高于饱和压力不产生溶解气从油中分离的过程单相流体的渗流问题。编辑ppt第一节非活塞式水驱油分析一、油水两相渗流的基本概念1、两相渗流区的形成形成两相区的原因：毛细管力的影响重率差的影响粘度差的影响编辑ppt（A）毛细管力的影响由于界面张力和岩石的润湿性所产生的毛管力有时是流动的阻力，有时是动力。（a）若岩石表面是亲油的，毛管力是阻力。P1P2流动方向水油Pc式中：——表面张力——润湿接触角r——毛管半径编辑ppt（b）若岩石表面是亲油的，毛管力是动力。P1P2流动方向水油Pc当毛管两端没有建立压差时（P1-P2=0），水在毛管力作用下也能渗入毛管。小毛管中毛管压力大，水首先渗入小毛管形成非活塞式推进。若两端建立压差P1-P2〉0，这种差别仍有可能存在。只有当所建立的压差P1-P2〉0大大地超过毛管力时，水主要靠外来压差渗入毛管，毛管力的影响就不明显了。编辑ppt（B）重率差的影响水比油重，因此油水相遇时，水向下，油向上，形成上油下水的两相区。当油水重率差很大，油层很厚，液流速度不大时，这种上油下水的两相区很容易形成。编辑ppt（C）粘度差的影响通常，油水粘度差异是比较大的。W=1mpa•so=3~10mpa•s水的流动比油的流动要容易得多。在外压差的作用下，由于大毛管通道横截面积大，阻力小，因而水首先渗入大毛管；又由于o>>W，水渗入的毛管中，总阻力下降，因而水窜越来越快，形成严重的指进现象。编辑ppt2、两相区的存在增大了渗流阻力（与活塞式驱油相比）水油11活塞式水油1122非活塞式水油1122油+水<渗流阻力相渗透率与绝对渗透率的关系：

Ko+Kw<K两相的渗流阻力大于纯水区（纯油区）的渗流阻力编辑ppt3、影响两相区渗流阻力的因素两相区渗流阻力的大小取决于流体的粘度和两相区的渗透率。两相区的渗透率用相渗透率表示。相渗透率是含油含水饱和度的函数。编辑ppt二、两相区中油水饱和度的变化规律——贝克莱-列维尔特驱油理论1、实验观测结果当原始油水界面垂直于流线，含油区内束缚水含量为常数时，两相区内油水饱和度分布如图所示，沿流程含水饱和度Sw逐渐变小，含油饱和度逐渐升高，在两相区前缘X=Xf处，含水饱和度曲线突然降落。编辑pptSX水区两相区油区SorSoSwSwrSwfSofSor：残余油饱和度So：可流动的含油饱和度Sw：含水饱和度Swr：束缚水饱和度Sof：油水前缘可流动的含油饱和度编辑ppt2、数理分析（A）任一过水断面上的含水率fw油水两相的运动方程：Ko、Kw——岩层对油相和对水相的渗透率。编辑ppt在岩层的过水断面上，水流量和油流量可表示为：任一过水断面上的总流量可表示为：编辑ppt任一过水断面上，水流量与总流量之比称为该过水断面上的含水率。编辑ppt任一过水断面上含油率可表示为：含水率、含油率是含水饱和度的函数。编辑ppt（B）两相区中含水饱和度的分布从两相区中取一微小单元体，先考查小单元体中含水饱和度的变化，再考虑两相区中Sw的分布。XAdx编辑pptAdxfw1fw2在dt时间内流入单元体中的水量Qw1与流出水量Qw2之差，应等于dt时间内，单元体中含水量的变化。式中：q(t)——时刻t，通过任一过水断面的液流量；fw1、fw2——单元体入口、出口端面上的含水率；Sw2-Sw1=Sw——dt时间内，单元体中含水饱和度的变化。编辑ppt上式表示某一固定含水饱和度的前移速度，称为贝克莱-列维尔特方程。进行积分：编辑ppt——从两相区开始形成（t=0）到时刻t为止，渗入油区的总水量。它实际上等于排液道（或井排）生产至t时刻的总产量。给定Swfw’(Sw)对应的X编辑ppt计算两相区中含水饱和度分布的步骤：1、求fw~Sw关系曲线由相对渗透率曲线求相渗透率。相对渗透率SwKroKrw求含水率。编辑ppt2、绘制fw’(Sw)~Sw的关系曲线Swfwfw’3、计算两相区中含水饱和度分布Swfw’XSw1fw1’Sw2fw2’Sw3fw3’编辑ppt（C）求前缘饱和度及前缘位置水油+水油XoX