Геолого-промысловая характеристика Уренгойского месторождения, страница 11

на период 2000-2005г.г.

Таблица 7.4.- Расчет показателей разработки УКПГ-15 на период2000-2005г.г..

Годы	Qст	Qдоб	P, МПа	Z	qв	QВ	H
2000	15	233,5	5,752286	0,897485	45	407	13,9
2001	15	248,5	5,173172	0,907859	51	458	15,6
2002	15	263,5	4,585421	0,91946	53	511	16,8
2003	15	278,5	3,986832	0,932384	52	563	18,02
2004	15	293,5	3,375006	0,946752	54	607	19,57
2005	15	308,5	2,984261	0,94893	54	661	21,1

Рис.   -Прогноз изменнения пластового давления на 2000-2005г.г..

Рис    .-Прогноз добычи газа на 2000-2005г.г.

Рис.   -Прогноз сренего дебита скважин на период с 2000=2005гг.

Приложение 1. Программа для расчета показателей разработки при упруговодонапорном режиме.

USES CRT;

Const

  Qzap=634552;

  m=0.3;

  a_kriv=0.75;

  a_ost=0.21;

  k=3.5e-13;

  x=0.9;

  Kv=0.21;

  Mv=8e-4;

  Tpl=301;

  h=50;

  Acp=0.00357;

  aOmegaN=5150.24;

  Bcp=5.25e-6;

  Tct=273;

  Pn=11.6;

  Scp=0.47;

  E=0.05e+6;

  a=4.26e-5;

  b=-0.0095;

  g=0.95;

  Pat=0.1013;

  God=31536000;

Label

1,2;

Var

                  Qv,Rz,qvMinus,PLast,Fo,PFo,R,qvPlus:Real;

                  P_kriv,Zn,Z:Real;

                  zxc,Discriminant:Real;{promezhutochnie peremennie}

                  Qdob:longint;

procedure qwe;

begin

                    zxc:=1/(aOmegaN-Qv) *(Pn/Zn*aOmegaN-Pat*Tpl/Tct*Qdob);{P/Z}

                    Discriminant:=SQRT(SQR(b-1/zxc)-4*a*g);

                    P_kriv:=(-Discriminant+1/zxc-b)/2/a;

                    Z:=a*P_kriv*P_kriv+b*P_kriv+g;

                    R:=SQRT(Rz*Rz-Qv/(3.14*m*h*(a_kriv-a_ost)));

                    Fo:=x/Rz/Rz*(God);

                    PFo:=(sqrt(12*Fo)*ln(sqrt(12*Fo)-1)+1)/(sqrt(12*Fo)-1);

end;

Begin

                  clrscr;

                  Rz:=SQRT((aOmegaN*1000000)/(a_kriv*3.14*m*h));

                  WRITELN('Rz=',Rz);

                  readln;

                  Qv:=0;

                  qvMinus:=0;

                  PLast:=0;

                  Zn:=a*Pn*Pn+b*Pn+g;

1:               write('VVEDITE Qdob     ');

                  readln(Qdob);

                  writeln;

                  {cycle}

          Fo:=x/Rz/Rz*(God);

                  PLast:=P_kriv;

                  PFo:=(sqrt(12*Fo)*ln(sqrt(12*Fo)-1)+1)/(sqrt(12*Fo)-1);

                    { Find P_kriv i Z}

                    zxc:=1/(aOmegaN-Qv) *(Pn/Zn*aOmegaN-Pat*Tpl/Tct*Qdob);{P/Z}

                    Discriminant:=SQRT(SQR(b-1/zxc)-4*a*g);

                    P_kriv:=(-Discriminant+1/zxc-b)/2/a;

                    Z:=a*P_kriv*P_kriv+b*P_kriv+g;

                    {Radius zalezhi}

                    R:=SQRT(Rz*Rz-Qv/(3.14*m*h*(a_kriv-a_ost)));

2:               qvPlus:=(2*3.14*k*h*(Pn-P_kriv)*1000000)/(Mv*(PFo+ln(Rz/R)/Kv));

                  Qv:=Qv+(qvMinus+qvPlus)/2;

                  qvMinus:=qvPlus;

                  QWE;

                  if  ABS(P_kriv-PLast)<=E then

                  begin

                  writeln;

                  writeln('qv=',qvPlus);

                  writeln('Qv=',Qv);

                  writeln('P=',P_kriv);

                  writeln('R=',R);

                  writeln('Z=',z);

                  writeln('Fo=',Fo);

                  writeln('PFo=',PFo);

                  writeln;

                  Goto 1

                  end

                  else Goto 2;

End.