水侵量-高等渗流力学

1、精选优质文档-倾情为你奉上精选优质文档-倾情为你奉上专心-专注-专业专心-专注-专业精选优质文档-倾情为你奉上专心-专注-专业 MACROBUTTON MTEditEquationSection2 Equation Chapter 1 Section 1 SEQ MTEqn r h * MERGEFORMAT SEQ MTSec r 1 h * MERGEFORMAT SEQ MTChap r 1 h * MERGEFORMAT 油藏天然水侵量的计算1问题简述图1 油藏示意图如图1所示，圆形油藏周围有边水，在油藏生产过程中，油藏平均压力是变化的，求解在这种情况下的累积水侵量。2技术思路将油藏看

2、成一口井，但这口井的井底压力（即油藏的平均压力）是不断变化的。因此，该问题是变压力条件下求产量的定解问题。先求得定压力下的产量解，然后根据杜哈美原理求得变压力下的累积产量，即求得了油藏的累积水侵量。3理论公式推导第一步：油藏边界上压力为常数时的压降解设油藏边界上压力为常数时的解（初始压力为）为，其无因次数学模型为： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 1)其中，无因次半径： MACROBUTTON MTPlaceRef * MERG

3、EFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 2)无因次时间： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 3)对上述数学模型进行Laplace变换，并令： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFO

4、RMAT 4)先对渗流方程进行变换：得到： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 5)对内边界条件进行变换，得到： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 6)对无限大边水区域情况的外边界条件进行变换，得到： MACROBUTTON MTPlaceRef * MERGEFORM

5、AT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 7)对有限封闭边水区域情况的外边界条件进行变换，得到： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 8)对 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (5)，两边同乘，然后第一项分子分母同乘，第二项分子分母

6、同乘： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 9)为0阶虚宗量的Bessel方程，通解为： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 10)第一种外边界条件：无限大边水区域的情况有 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Ch

7、arformat ! * MERGEFORMAT (7)，且，由的渐近性知，当，得到，得到： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 11)将内边界条件 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (6)带入 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * M

8、ERGEFORMAT (11)得到： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 12)第二种外边界条件：有限封闭边水区域的情况将内边界条件 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (6)和外边界条件 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGE

9、FORMAT (8)带入 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (10)，得到：第二步：水侵量与压力差的关系由达西定律： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 13)物理意义为：单位压差、单位厚度条件下的水侵速度。则有，累积水侵量为： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn

10、 h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 14)根据 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (2) GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (3)，则有： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * A

11、rabic * MERGEFORMAT 15)其中， MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 16)为单位压差条件下的水侵量，若压差，则有。第三步：计算对 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (16)进行Laplace变换，得到： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h

12、 * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 17)第一种外边界：对于无限大边水区域的情况：联立 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (12)和 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (17)得到： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORM

13、AT ( SEQ MTEqn c * Arabic * MERGEFORMAT 18)第二种外边界：对于有限封闭边水区域的情况：联立 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (10)和 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (17)得到： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MT

14、Eqn c * Arabic * MERGEFORMAT 19)用Stehfest数值反演方法对 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (18)和 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (19)进行反演得到（matlab程序见附件）。第四步：计算变压差条件下的水侵量由杜哈美原理： MACROBUTTON MTPlaceRef * MERGEFORMAT SEQ MTEqn h

15、* MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 20)上面求得的累积产量是在压力降为常数的条件下的计算的累积水侵量。但是在一般情况下，随着油藏的开采其平均压力是变化的。若压力降变化可由解析表达式描述，则可由 GOTOBUTTON ZEqnNum * MERGEFORMAT REF ZEqnNum * Charformat ! * MERGEFORMAT (20)来计算，若压力差不能由表达式表示，则将油藏平均压力的变化用一系列稳定压力台阶表示，每一台阶内认为油藏压力是不变的，则水侵量可由下式计算： MACROBUTTON MTPlaceRef

16、 * MERGEFORMAT SEQ MTEqn h * MERGEFORMAT ( SEQ MTEqn c * Arabic * MERGEFORMAT 21)4计算步骤框图图2 水侵量计算步骤框图5计算结果分析图3 无因次累积水侵量曲线（有限水区和无限大水区）结果分析：在时间一定的情况下，越大，即边水区域越大，无因次累积水侵量越大，当时，无因次累积水侵量最大。当一定时，在前期，无因次时间越大，无因次累积水侵量越大，当到达一定时间后，无因次累积水侵量不再增加，此时油藏平均压力和水区平均压力达到稳定，水侵过程停止。图4 变压差下的累积水侵量曲线结果分析：在时间一定时，越大，即边水区域越大，累积

17、水侵量越大。在一定时，时间越大，累积水侵量越大。6结论本文通过理论推导和数值求解得到了无因次累积水侵量曲线和变压差条件下的累积水侵量曲线并进行了分析。源程序%第三步：画出真实空间中的无因次水侵量随无因次时间变化的关系曲线%clear all%clc;function waterrDw=1.5:0.5:5 6:10;for i=1:13 tD=linspace(0.1,i2.5,300); for m=1:300 QtD(m)=QtD(tD(m),rDw(i); end loglog(tD,QtD) if(i=1) text(i2.5,QtD(m),rDw=Rw/Re=1.5); else s=

18、strcat(text(i2.5,num2str(QtD(m),num2str(rDw(i),); eval(s); end% s=strcat(text(i2.5,num2str(QtD(m),rDw=,num2str(rDw(i),);% eval(s); hold onendtD=linspace(0.1,132.5,300);for m=1:300 QtD1(m)=QtD1(tD(m);end loglog(tD,QtD1) text(132.5,QtD1(m),rDw无穷大);hold on %这个地方非常重要！xlabel(无因次时间tD);ylabel(无因次累积水侵量QtD);

19、title(无因次累积水侵量曲线QtD-tD（有限水区和无限大水区）);QT();end%第一步：无因次水侵量Laplace空间解（第一种外边界条件：有限封闭边水区域的情况）function Qs=Qs(s,rDw)C=s*(besselk(0,sqrt(s)*besseli(1,sqrt(s)*rDw)+besseli(0,sqrt(s)*besselk(1,sqrt(s)*rDw);A=besselk(1,sqrt(s)*rDw)/C;B=besseli(1,sqrt(s)*rDw)/C;Qs=(-A*besseli(1,sqrt(s)+B*besselk(1,sqrt(s)/sqrt(s

20、);end% %第一步：无因次水侵量Laplace空间解（第二种外边界条件：无限大边水区域的情况）function Qs1=Qs1(s) %只是s的函数，数值反演以后也只是tD的函数Qs1=besselk(1,sqrt(s)/(s1.5*besselk(0,sqrt(s);endfunction QtD1=QtD1(tD)sum1=0;N=8;for i=1:N sum2=0; for k=floor(i+1)/2):min(i,N/2) %sum2=sum2+k(N/2)*factorial(2*k+1)/(factorial(k+1)*factorial(k)*factorial(N/2-

21、k+1)*factorial(i-k+1)*factorial(2*k-i+1); sum2=sum2+k(N/2)*factorial(2*k)/(factorial(k)*factorial(k-1)*factorial(N/2-k)*factorial(i-k)*factorial(2*k-i); end Vi=(-1)(N/2+i)*sum2; s=i*(log(2)/tD; sum1=sum1+Vi*Qs1(s);endQtD1=(log(2)/tD*sum1;end%第二步：Stehfest数值反演（我的编程思路）function QtD=QtD(tD,rDw)sum1=0;N=8

22、;for i=1:N sum2=0; for k=floor(i+1)/2):min(i,N/2) %sum2=sum2+k(N/2)*factorial(2*k+1)/(factorial(k+1)*factorial(k)*factorial(N/2-k+1)*factorial(i-k+1)*factorial(2*k-i+1); sum2=sum2+k(N/2)*factorial(2*k)/(factorial(k)*factorial(k-1)*factorial(N/2-k)*factorial(i-k)*factorial(2*k-i); end Vi=(-1)(N/2+i)*

23、sum2; s=i*(log(2)/tD; sum1=sum1+Vi*Qs(s,rDw);endQtD=(log(2)/tD*sum1;end%第三步：画出真实空间中的无因次水侵量随无因次时间变化的关系曲线（用不到，把它放在最后不对就是不对？）% clear all;% clc;% rDw=1.5:0.5:5 6:10;% for i=1:13% t=linspace(0.1,i2.5,300);% for m=1:300% QtD(i,m)=QtD(t(m),rDw(i);% end% loglog(t,QtD(i,:)% hold on% end% hold off% xlabel(无因次时间tD);% ylabel(水侵量QtD);% title(无因次水侵量曲线QtD-tD);%第四