解剖黑油模型源程序匯報(bào)坤

1、油藏?cái)?shù)值模擬之黑油模型解剖答辯人：孟凡坤2014.05子函數(shù)UINIT1（4）子函數(shù)CGIN及MF1（10、11）子函數(shù)MATBAL（12）匯報(bào)內(nèi)容油藏?cái)?shù)值模擬之黑油模型解剖油藏?cái)?shù)值模擬之黑油模型解剖（UINIT1)UINIT1作用：給模型賦初值開始KREST=1讀入P(I,J,K), SO(I,J,K),SW(I,J,K)SG(I,J,K),PBOT(I,J,K)讀入重啟運(yùn)算起始時(shí)間ETIR,再啟動(dòng)起始步數(shù)NRESTT,再啟動(dòng)初始平均壓PAVGO讀入第I口井在第K層的產(chǎn)油量CUMO，產(chǎn)水量CUMW，產(chǎn)氣量CUMG結(jié)束讀入NVQN1活動(dòng)井?dāng)?shù)YES讀入OOIP, OWIP, ODGIP, OFG

2、IPKPI=0 讀入PI, PGOC, WOC, GOCEL(I,J,K)GOCWOC計(jì)算PI, BBW,RSW等計(jì)算PGOC,BBG等計(jì)算BPT, BBO, RSO等計(jì)算PN(I,J,K)KSI=0 讀入SO,SW,SG=1-SO-SWSO=SOI,SW=SWI,SG=SGI結(jié)束NO讀入PN321NOYES初始啟動(dòng)(KREST=0)：主要以不同方式給各個(gè)網(wǎng)格點(diǎn)賦初始?jí)毫惋柡投取PI=0，原始?jí)毫Ω鶕?jù)油水界面壓力、高度差和流體密度折算得到:KPI=1，各網(wǎng)格點(diǎn)原始?jí)毫τ蓴?shù)據(jù)文件直接讀入。KSI=0,則各層油氣水均為定值； KSI=1,則各網(wǎng)格點(diǎn)原始油氣水飽和度不同，由數(shù)據(jù)文件直接讀入。再啟

3、動(dòng)方式(KREST=1)：飽和度和壓力值由再啟動(dòng)文件讀入，還應(yīng)讀入再啟動(dòng)起始點(diǎn)的儲(chǔ)量、累計(jì)產(chǎn)量，起始時(shí)刻及步數(shù)等。1：ELGOC2：GOCELWOC油藏?cái)?shù)值模擬之黑油模型解剖（CGIN、MF1)CGIN：預(yù)處理共軛梯度法(ORTHOMIN)求解壓力其中，M為矩陣A進(jìn)行不完全LU分解后得到的近似矩陣，M為預(yù)處理矩陣，根據(jù)By=b，其中B= AM-1，求得y，之后再由y=Mx，求得x。開始DE(IC)=Diag(A)=E(I,J,K)QVEC(IC)=B(I,J,K)PV(IC)=P(I,J,K)N2=0,A M=LUN1=0R= QVEC-A* PVEC, Q=M-1RN1=N1+1N2=N2+

4、1SV1=SV2=0AQ=A*Q, SV3=QTAQMR= M-1R, SV1=R*MRAL= SV1/SV3, DP=0PVEC = PVEC +AL*Q, R=R-AL*AQDP=max(AL*Q)iDP200NON1=EFC & NKNK r=b-Ax;q=M-1r; do k=0,KNK-1 ak=(r,M-1r)/(q,Aq) x=x+aq; r=r-aAq; dx=ak*q; if(dxKNK) exit. b=(r,M-1r)/(r,M-1r) q=M-1r+bq N=N+1NZ=111NN1=501矩陣形式：求解公式油藏?cái)?shù)值模擬之黑油模型解剖（MATBAL)對(duì)每一步末計(jì)算結(jié)果物

5、質(zhì)平衡檢驗(yàn)，確定計(jì)算誤差。這一計(jì)算分別對(duì)油、氣、水進(jìn)行。開始變量聲明i=1,j=1,k=1計(jì)算PAVGO,PAVG和OP計(jì)算時(shí)步內(nèi)WP，WI，GP和GIi=i+1,j=j+1,k=k+1NOiII,jJJ,kKK YES計(jì)算油藏COP, CWP, CGP, CWI, CGI計(jì)算油藏OPR, WPR, GPR, WIR, GIR計(jì)算物質(zhì)平衡誤差MBEO, MBEW, MBEG計(jì)算時(shí)步內(nèi)GOR,WOR, CGOR, CWOR結(jié)束PAVG0=0 時(shí)步初油藏平均壓力；PAVG=0 時(shí)步末油藏平均壓力OP=0 時(shí)步內(nèi)累積產(chǎn)油量；WP=0 時(shí)步內(nèi)累積產(chǎn)水量GP=0 時(shí)步內(nèi)累積產(chǎn)氣量；WI=0 時(shí)步內(nèi)累積注

6、水量GI=0 時(shí)步內(nèi)累積注氣量PAVG0=PAVG0+PN(I,J,K)*VP(I,J,K)PAVG=PAVG+P(I,J,K)*VP(I,J,K)OP=OP+QO(I,J,K)*FACT；VP(I,J,K): 各網(wǎng)格節(jié)點(diǎn)孔隙度QW0,WP=WP+QW(I,J,K)*FACT；QW0,GP=GP+QG(I,J,K)*DELT0QG0,GI=GI+QG(I,J,K)*DELT0MBEO=(STBO/(STBWI-OP)-1.0)*100.0STBOI:時(shí)步初油相體積；STBO:時(shí)步末油相體積MBEW=(STBW/(STBWI-WP-WI)- 1.0)*100.0STBWI:時(shí)步初油相體積；STBW:時(shí)步末