资源描述：
《水侵量--高等渗流力学.docx》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、油藏天然水侵量的计算1问题简述图1油藏示意图如图1所示，圆形油藏周围有边水，在油藏生产过程中，油藏平均压力是变化的，求解在这种情况下的累积水侵量。2技术思路将油藏看成一口井，但这口井的井底压力（即油藏的平均压力）是不断变化的。因此，该问题是变压力条件下求产量的定解问题。先求得定压力下的产量解，然后根据杜哈美原理求得变压力下的累积产量，即求得了油藏的累积水侵量。3理论公式推导第一步：油藏边界上压力为常数时的压降解设油藏边界上压力为常数时的解（初始压力为）为，其无因次数学模型为：(1)其中，无因次半径：(2)无因次时间：(3)对上述数学模型
2、进行Laplace变换，并令：(4)先对渗流方程进行变换：得到：(5)对内边界条件进行变换，得到：(6)对无限大边水区域情况的外边界条件进行变换，得到：(7)对有限封闭边水区域情况的外边界条件进行变换，得到：(8)对，两边同乘，然后第一项分子分母同乘，第二项分子分母同乘：(9)为0阶虚宗量的Bessel方程，通解为：(10)第一种外边界条件：无限大边水区域的情况有，且，由的渐近性知，当，，得到，得到：(11)将内边界条件带入得到：(12)第二种外边界条件：有限封闭边水区域的情况将内边界条件和外边界条件带入，得到：第二步：水侵量与压力差的
3、关系由达西定律：(13)物理意义为：单位压差、单位厚度条件下的水侵速度。则有，累积水侵量为：(14)根据，则有：(15)其中，(16)为单位压差条件下的水侵量，若压差，则有。第三步：计算对进行Laplace变换，得到：(17)第一种外边界：对于无限大边水区域的情况：联立和得到：(18)第二种外边界：对于有限封闭边水区域的情况：联立和得到：(19)用Stehfest数值反演方法对和进行反演得到（matlab程序见附件）。第四步：计算变压差条件下的水侵量由杜哈美原理：(20)上面求得的累积产量是在压力降为常数的条件下的计算的累积水侵量。但是
4、在一般情况下，随着油藏的开采其平均压力是变化的。若压力降变化可由解析表达式描述，则可由来计算，若压力差不能由表达式表示，则将油藏平均压力的变化用一系列稳定压力台阶表示，每一台阶内认为油藏压力是不变的，则水侵量可由下式计算：(21)4计算步骤框图图2水侵量计算步骤框图5计算结果分析图3无因次累积水侵量曲线（有限水区和无限大水区）结果分析：在时间一定的情况下，越大，即边水区域越大，无因次累积水侵量越大，当时，无因次累积水侵量最大。当一定时，在前期，无因次时间越大，无因次累积水侵量越大，当到达一定时间后，无因次累积水侵量不再增加，此时油藏平均
5、压力和水区平均压力达到稳定，水侵过程停止。图4变压差下的累积水侵量曲线结果分析：在时间一定时，越大，即边水区域越大，累积水侵量越大。在一定时，时间越大，累积水侵量越大。6结论本文通过理论推导和数值求解得到了无因次累积水侵量曲线和变压差条件下的累积水侵量曲线并进行了分析。源程序%第三步：画出真实空间中的无因次水侵量随无因次时间变化的关系曲线%clearall%clc;functionwaterrDw=[1.5:0.5:56:10];fori=1:13tD=linspace(0.1,i^2.5,300);form=1:300QtD(m)=Q
6、tD(tD(m),rDw(i));endloglog(tD,QtD)if(i==1)text(i^2.5,QtD(m),'rDw=Rw/Re=1.5');elses=strcat('text(i^2.5,',num2str(QtD(m)),',''',num2str(rDw(i)),''')');eval(s);end%s=strcat('text(i^2.5,',num2str(QtD(m)),',''rDw=',num2str(rDw(i)),''')');%eval(s);holdonendtD=linspace(0.1,13^2
7、.5,300);form=1:300QtD1(m)=QtD1(tD(m));endloglog(tD,QtD1)text(13^2.5,QtD1(m),'rDw无穷大');holdon%这个地方非常重要！！！xlabel('无因次时间tD');ylabel('无因次累积水侵量QtD');title('无因次累积水侵量曲线QtD-tD（有限水区和无限大水区）');QT();end%第一步：无因次水侵量Laplace空间解（第一种外边界条件：有限封闭边水区域的情况）functionQs=Qs(s,rDw)C=s*(besselk(0,sqr
8、t(s))*besseli(1,sqrt(s)*rDw)+besseli(0,sqrt(s))*besselk(1,sqrt(s)*rDw));A=besselk(1,sqrt(s)*rDw)/C;B=be
