newdata.m

来自「Petroleum Reservoir Simulation」· M 代码 · 共 80 行

dpkp1=triblocksolve(JM,-fk,2); %get the dPk+1
%dpkp1=gaus(JM,-fk'); %get the dPk+1
fkim1=fk;   %fki-1


%re_fill the data

%%%%%%%%%%%%%%%%%%%%%%pwit=pwi;
%%%%%%%%%%%%%%%%%%%%%%swit=swi;
%%%%%%%%%%%%%%%%%%%%%%phit=phi;
   
for ii=1:nr
   dpk(ii)=dpkp1(2*ii-1);
   dsk(ii)=dpkp1(2*ii);
 
   pwi(ii)=pwi(ii)+dpk(ii);
   swi(ii)=swi(ii)+dsk(ii);
   phi(ii)=phir.*exp(cf.*(pwi(ii)-pr));
   
   
   
   
%change this if something is wrong
   %if ii==nr
      %pwi(ii)=pwf;
   %end
   
   
   
   
   
   
   
   
   %if ii==1
      %swi(ii)=1;
   %end
   
end


get_start;

dpkp1=dpkp1';
for ii=1:(2*nr)
   if abs(dpkp1(ii))>ero(ii)|abs(fki(ii)-fkim1(ii))>erof(ii)
      e=1;
      break;
   else
      e=0;
   end
   
end