渗流力学实验.doc

《渗流力学实验.doc》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、渗流力学上机实验报告实验三实验名称：理想系统不稳定渗流井壁压力计算实验目的：了解渗流力学中最有代表性的三种渗流方式：平面线性渗流、平面径向渗流、空间球形渗流的求解方法、井壁压力及其导数的变化规律。实验原理：1、渗流力学中常用的无量纲量（SI单位制）举例：2、按此无量纲定义，以径向渗流为例，其无量纲控制方程组为：初始条件：内边界条件：外边界条件：3、压力导数渗流力学常用的压力导数定义为：对于不稳定渗流过程，压力导数比压力更敏感，能够比较清晰地反映流体渗流的方式。实验内容：1、平面线性渗流压力分布井壁压力及其压力导数，2、平面径向渗流压力及压力导数为3、球形空间渗流附

2、图如下实验四：实验名称：一维单向饱和度场数学模拟实验目的：理解油水两相渗流的物理意义，以及在水驱油过程中相对渗透率变化规律、含水率的变化规律以及饱和度的分布规律。原理公式：附图如下附：实验三代码OptionExplicitPublicFunctionVBPwt2(tdAsDouble)AsDoubleVBPwt2=0.5*VBEix(0.25/td)EndFunctionPublicFunctionVBEix(xAsDouble)AsDoubleDimoAsDouble,uAsDoubleIf(x<=1)Theno=-0.57721566+0.99999193*x

3、-0.24991055*x^2+0.05519968*x^3-0.00976885*x^4+0.00127857*x^5o=o-Log(x)Elseo=0.2677737343+x*(8.6347608925+x*(18.059016973+x*(8.5733287401+x)))u=3.9584969228+x*(21.0996530827+x*(25.6329561486+x*(9.5733223454*x)))o=(o/u)*Exp(-x)/xEndIfVBEix=oEndFunctionPublicFunctiondppei(tdAsDouble)AsDo

4、ubledppei=0.5*Exp(-0.25/td)EndFunctionPublicFunctionpderf(xAsDouble)AsDoublepderf=1-Exp(-x^2)*((0.254829592/(1+0.327951*x))+(-0.284496736/(1+0.327951*x)^2)+(1.421413741/(1+0.327951*x)^3)+(-1.453152027/(1+0.327951*x)^4)+(1.061405429/(1+0.327951*x)^5))+1.5*10^(-7)EndFunctionPublicFuncti

5、onpdcir(tdAsDouble)AsDoublepdcir=1-pderf(1/2/Sqr(td))EndFunctionPublicFunctiondpcir(tdAsDouble)AsDoubledpcir=td^(-0.5)/Sqr(2*3.141592653)*Exp(-1/(4*td))EndFunctionPublicFunctionPd(tdAsDouble)AsDoublePd=2*Sqr(td/3.141592653)EndFunctionPublicFunctiondPd(tdAsDouble)AsDoubledPd=Sqr(td/3.1

6、41592653)EndFunction