滲流力學(xué)實驗

滲流力學(xué)上機實驗報告實驗三實驗名稱：理想系統(tǒng)不穩(wěn)定滲流井壁壓力計算實驗?zāi)康模毫私鉂B流力學(xué)中最有代表性的三種滲流方式：平面線性滲流、平面徑向滲流、空間球形滲流的求解方法、井壁壓力及其導(dǎo)數(shù)的變化規(guī)律。實驗原理：1、滲流力學(xué)中常用的無量綱量(SI單位制)舉例：kh(p-p) p-p 1.842x10-3qpB"D1.842x10-3qpB，「Dp-p&d kh(p-p)iw iwt=3知kt;C=尸板;r='D 9pcr2D9chr2Drtw tw w2、按此無量綱定義，以徑向滲流為例，其無量綱控制方程組為:d2p*1dpdp

dr2rdrdt初始條件:pD")=0內(nèi)邊界條件:dpI1

r—d-I =-1DdrIDrD外邊界條件:p(3,t)=03、壓力導(dǎo)數(shù)滲流力學(xué)常用的壓力導(dǎo)數(shù)定義為:

dp_tdp

dIntdt對于不穩(wěn)定滲流過程，壓力導(dǎo)數(shù)比壓力更敏感，能夠比較清晰地反映流體滲流的方式。實驗內(nèi)容:實驗內(nèi)容:1、平面線性滲流壓力分布p(x,t1、平面線性滲流壓力分布p(x,t)_xm_2\Sexp(-Xd)-xerfc(*」)DDDDd 氣..兀 4t D 2jtD D井壁壓力及其壓力導(dǎo)數(shù)tD兀dp't

w^_■-D

dInt ''兀D2、平面徑向滲流壓力及壓力導(dǎo)數(shù)為dpD—dlntDp(r,t)_-1Ei(-，)

DDD2 4tDr-Le4%”2 1-r2——t ( —)_e4tDDrD 4tD 2D3、球形空間滲流rp(rp(r,t)_1-erf(DDDD'wD_ =t-0.5e4tDdt。 2&d,附圖如下wooC.Ol0.i1IDJGO1000 100G0IDOOOD實驗四：實驗名稱：一維單向飽和度場數(shù)學(xué)模擬實驗?zāi)康模豪斫庥退畠上酀B流的物理意義，以及在水驅(qū)油過程中相對滲透率變化規(guī)律、含水率的變化規(guī)律以及飽和度的分布規(guī)律。原理公式：犬=%［心沖=0。為累計注入量E注=吾血"鈕⑺附圖如下附：實驗三代碼OptionExplicitPublicFunctionVBPwt2(tdAsDouble)AsDoubleVBPwt2=0.5*VBEix(0.25/td)EndFunctionPublicFunctionVBEix(xAsDouble)AsDoubleDimoAsDouble,uAsDoubleIf(x<=1)Theno=-0.57721566+0.99999193*x-0.24991055*xA2+0.05519968*xA3-0.00976885*xA4+0.00127857*xa5o=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)AsDoubledppei=0.5*Exp(-0.25/td)EndFunctionPublicFunctionpderf(xAsDouble)AsDoublepderf=1-Exp(-xa2)*((0.254829592/(1+0.327951*x))+(-0.284496736/(1+0.327951*x)a2)+(1.421413741/(1+0.327951*x)a3)+(-1.453152027/(1+0.327951*x)a4)+(1.061405429/(1+0.327951*x)a5))+1.5*10a(-7)EndFunctionPublicFunctionpdcir(tdAsDouble)AsDoublepdcir=1-pderf(1/2/Sqr(td))EndFunctionPublicFunctiondpcir(tdAsDouble)AsDoubledpcir=tdA(-0.5)/Sqr(2*3.141592653)*Exp(-1/(4*td))EndFunctionPublicFunctionPd(tdAsDouble)AsDoublePd=2*Sq