(完整word)油藏数值模拟隐式差分MATLAB源程序.doc

％隐式求解方法

％t 为投产后某一时刻,单位：天

％d:迭代时间；

%Pwf1：W1 井底流压；

％Q2：W2 井产油量；

function ［P，d,Pwf1,Q2]=implict(t)

%油藏参数

dy=200；

dt=24;

n=t*24/dt；％迭代时间步数

re=0。208＊dx；

rw=0。1;

％渗透率

K=[0 259 222 200 190 180 185 0 0 0 0;259 259 222 200 190 180 185 185 0 0 0；310 310 240 235

228 210 195 195 0 0 0;330 330 290 270 250 230 205 197.5 180 185 0;350 350 300 280 259 222 200

190 180 185 185;340 340 320 290 310 240 235 228 210 195 195；355 355 335 315 310 290 270 250

230 205 205;0 0 0 0 325 300 280 240 210 215 215；0 0 0 0 340 320 290 260 235 225 225;0 0 0 0

P=Pini＊ones（10,11）;

％P（1，:)=0；

％P（2:10，1）=0；

％P（2：10，11）=0;

％P(2：4,8:10)=0；

％P(8：10,2：4）=0;

%系数矩阵

%初始化

a=zeros（10,11）;

b=zeros(10,11);

a（i，j)=3600e-9*dt＊2＊H(i,j—1)＊K（i,j-1)＊K(i，j）/u/C/Fai(i，j）/dx/dx/（H（i，

j—1）＊K（i，j-1）+H(i,j）＊K（i,j)）；

b(i，j）=3600e—9*dt＊2＊H(i，j+1）*K（i,j+1）＊K（i，j)/u/C/Fai（i，j)/dx/dx/（H(i,j+1)*K