;;:=true; i:=N - 1 downto 2 doj:=M - 1 downto 2 doh [i, j] < ;> 0 then ^ [i, j]:=V ^ [i, j]-b1 ^ [i, j] * dp ^ [i +1, j]-f1 ^ [i, j] * dp ^ [ i, j +1]; abs (dp ^ [i, j])> 1e - 4 then flag:=false; ^ [i, j]:=pp ^ [i, j] + dp ^ [i, J];; i:=1 to N doj:=1 to M do begin b1 ^ [i, j]:=0; f1 ^ [i, j]:=0 end;; flag;:=0; Pcp :=0; Hcp:=0; i:=1 to N doj:=1 to M do
Продовження додатка А
gi1 [i, j]:=pp ^ [i, j];
ZapT:=ZapT+h[i,j]*dx*dy*kp*al1*p[i,j]/pat/z(p[i,j]/1e6)*z(pat/1e6);:=Qdob+qm[i,j]*dt;h[i,j]<>0 then begin Pcp:=Pcp + pp ^ [i, j] * h [i, j]; Hcp:=Hcp + h [i, j] end;;; (ZapN-ZapT-Qdob) / Qdob * 100> delta then delta:=(ZapN-ZapT-Qdob)/Qdob*100;:=Pcp/Hcp;:=Pcp;_1:=pp^[9,13];_1:=pp^[9,26];_1:=pp^[5,15];_1:=pp^[6,29];_1:=pp^[8,5];(kh,M*N*SizeOf(tipR));(b1,M*N*SizeOf(tipR));(pp,M*N*SizeOf(tipR));(dp,M*N*SizeOf(tipR));(V,M*N*SizeOf(tipR));(g1,M*N*SizeOf(tipR));(f1,M*N*SizeOf(tipR));;soprVIIa(gi1:m1;dt,t:real; var u5_1, u8_1, u10_1, u16_1, u75_1: real;
e5_1, e5_2, e8_1, e8_2, e10_1, e10_2, e16_1, e16_2, e75_1, e75_2: real);
const=16; {strok}=34; {stolbcov}
Ppl=7.8e6; {Pa}=1.877e - 12; {mkm2, IV block}=2.062e - 12; {mkm2, III block}=0.3;=0.78;=0.76;=0.011e - 3; {Pa * sec}
dx=100; {m}=100; {m}=100000;=86400;=real;=array [1 .. 10] of tipR;=array [1 .. N, 1 .. M] of tipR;
dmas=^ mas2;, dat1: text;
al:mas1;,j,count,k,cell:integer;,p,h,u:mas2;,c,a,b,f,d,R,sum,ZapN,ZapT,Qdob,Pcp,delta,Hcp,ggg,fff,bbb,v1:tipR;,g1,f1,kh,V,dp,pp:dmas;:boolean;alf(var al: mas1);
var, j: integer;, par1, par11, par2: tipR;
begin:=100; i:=1 to N - 1 doj:=1 to M - 1 doh [i, j] <> 0 then (h [i, j +1] <> 0 ) and (h [i +1, j] <> 0) then:=dx * dx * (kh ^ [i, j +1] + kh ^ [i, j]) / (kh ^ [i + 1, j] + kh ^ [i, j]) / dy / dy;
par1:=pi * pi / (2 * N * N * (1 + par));:=pi * pi / (2 * M * M * (1 +1 / par));
if par1
{5} [9,13]:=e5_1 {/ 86.4};
{75}
qm [8,5]:=e75_1 {/ 86.4};
{10}
Продовження додатка А
qm [5,15]:=e10_1 {/ 86.4};
{8}
qm [9,26]:=e8_1 {/ 86.4};
{16}
qm [6,29]:=e16_1 {/ 86.4};
end;
BEGIN
GetMem(kh,M*N*SizeOf(tipR));(b1,M*N*SizeOf(tipR));(pp,M*N*SizeOf(tipR));(dp,M*N*SizeOf(tipR));(V,M*N*SizeOf(tipR));(g1,M*N*SizeOf(tipR));(f1,M*N*SizeOf(tipR));:=0;i:=1 to N doj:=1 to M doh [i, j] <> 0 thenj <= 21 then kh^[i,j]:=k2*h[i,j]kh^[i,j]:=k1*h[i,j]kh^[i,j]:=0;[i,j]:=gi1[i,j];[i,j]:=0;^[i,j]:=u[i,j];[i,j]:=0;^[i,j]:=0;^[i,j]:=0;^[i,j]:=0;^[i,j]:=0;^[i,j]:=0;;;(al);:=0;Qdob:=0;delta:=0;:=0;
k:=0;;:=count +1; odd (count) then {obratnoe yporyadochivanie} i:=2 to N - 1 doj:=M - 1 downto 2 do:=0; c:=0; b:=0; f:=0; a:=0; d:=0; R:=0; h [i, j] <> 0 thenh [i, j - 1]=0 then g:= 0g:=(dx / dy) * p [i, j] * (kh ^ [i, j - 1] + kh ^ [i, j]) / mu / z (p [i, j] / 1e6) / 2; h [i - 1, j]=0 then c:=0c:=(dy / dx) * p [i, j] * (kh ^ [i - 1, j] + kh ^ [i, j] ) / mu / z (p [i, j] / 1e6) / 2; h [i +1, j]=0 then b:=0b:=(dy / dx) * p [i, j] * (kh ^ [i +1, j] + kh ^ [i, j]) / mu / z (p [i, j] / 1e6) / 2; h [i, j +1]=0 then f:=0f:=(dx/dy)*p[i,j]*(kh^[i,j+1]+kh^[i,j])/mu/z(p[i,j]/1e6)/2;:=-(g+c+b+f)-dx*dy/dt*al1*kp*h[i,j];:=pat*qm[i,j]/dx/dy-dx*dy/dt*al1*kp*h[i,j]*p[i,j]/z(p[i,j]/1e6)*u[i,j];:=d-...
