解剖黑油模型源程序汇报坤

1、油藏数值模拟之黑油模型解剖答辩人：孟凡坤2014.05子函数UINIT1（4）子函数CGIN及MF1（10、11）子函数MATBAL（12）汇报内容油藏数值模拟之黑油模型解剖油藏数值模拟之黑油模型解剖（UINIT1)UINIT1作用：给模型赋初值开始KREST=1读入P(I,J,K), SO(I,J,K),SW(I,J,K)SG(I,J,K),PBOT(I,J,K)读入重启运算起始时间ETIR,再启动起始步数NRESTT,再启动初始平均压PAVGO读入第I口井在第K层的产油量CUMO，产水量CUMW，产气量CUMG结束读入NVQN1活动井数YES读入OOIP, OWIP, ODGIP, OFG

2、IPKPI=0 读入PI, PGOC, WOC, GOCEL(I,J,K)GOCWOC计算PI, BBW,RSW等计算PGOC,BBG等计算BPT, BBO, RSO等计算PN(I,J,K)KSI=0 读入SO,SW,SG=1-SO-SWSO=SOI,SW=SWI,SG=SGI结束NO读入PN321NOYES初始启动(KREST=0)：主要以不同方式给各个网格点赋初始压力和饱和度。KPI=0，原始压力根据油水界面压力、高度差和流体密度折算得到:KPI=1，各网格点原始压力由数据文件直接读入。KSI=0,则各层油气水均为定值； KSI=1,则各网格点原始油气水饱和度不同，由数据文件直接读入。再启

3、动方式(KREST=1)：饱和度和压力值由再启动文件读入，还应读入再启动起始点的储量、累计产量，起始时刻及步数等。1：ELGOC2：GOCELWOC油藏数值模拟之黑油模型解剖（CGIN、MF1)CGIN：预处理共轭梯度法(ORTHOMIN)求解压力其中，M为矩阵A进行不完全LU分解后得到的近似矩阵，M为预处理矩阵，根据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矩阵形式：求解公式油藏数值模拟之黑油模型解剖（MATBAL)对每一步末计算结果物

5、质平衡检验，确定计算误差。这一计算分别对油、气、水进行。开始变量声明i=1,j=1,k=1计算PAVGO,PAVG和OP计算时步内WP，WI，GP和GIi=i+1,j=j+1,k=k+1NOiII,jJJ,kKK YES计算油藏COP, CWP, CGP, CWI, CGI计算油藏OPR, WPR, GPR, WIR, GIR计算物质平衡误差MBEO, MBEW, MBEG计算时步内GOR,WOR, CGOR, CWOR结束PAVG0=0 时步初油藏平均压力；PAVG=0 时步末油藏平均压力OP=0 时步内累积产油量；WP=0 时步内累积产水量GP=0 时步内累积产气量；WI=0 时步内累积注

6、水量GI=0 时步内累积注气量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): 各网格节点孔隙度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:时步初油相体积；STBO:时步末油相体积MBEW=(STBW/(STBWI-WP-WI)- 1.0)*100.0STBWI:时步初油相体积；STBW:时步末