Добавил:
Upload Опубликованный материал нарушает ваши авторские права? Сообщите нам.
Вуз: Предмет: Файл:
отчет по имр.doc
Скачиваний:
23
Добавлен:
25.04.2019
Размер:
494.08 Кб
Скачать

3.4 Блок-схема алгоритма решения задачи

µ1 = µi

Рисунок 3.2 – Блок-схема к заданию №3

Аналогичный алгоритм для вычисления µ2.

3.5 Текст программы

program imr3;

const n=10;

type mas1=array[1..n] of real;

var a1, b1, c1, a2, b2, c2, p1, p2, NY1, NY2:real;

f1, f2, f3:text;

P, NY:mas1;

i:integer;

begin

assign(f1,'dano.txt');

assign(f2,'tab.txt');

assign(f3,'vivod.txt');

reset(f1);

reset(f2);

rewrite(f3);

read (f1, p1, p2);

for i:=1 to n do

begin

readln(f2, P[i], NY[i]);

if (P1=P[i]) then

begin

NY1:=NY[i];

writeln(f3,' ny1=',ny1:6:6);

end;

if (P2=P[i]) then

begin

NY2:=NY[i];

writeln(f3,' ny2=',ny2:6:6);

end;

end;

for i:=2 to n-1 do

begin

if (P1>P[i-1]) and (P1<P[i+1]) then

begin

a1:=((ny[i+1]-ny[i-1])/(p[i+1]-p[i-1])-(ny[i]-ny[i-1])/(p[i]-p[i-1]))/(p[i+1]-p[i]);

b1:=(ny[i]-ny[i-1])/(p[i]-p[i-1])-a1*(p[i]+p[i-1]);

c1:=ny[i-1]-a1*sqr(p[i-1])-b1*p[i-1];

NY1:=a1*sqr(p1)+b1*p1+c1;

writeln(f3,' a1=',a1:6:6,' b1=',b1:6:6,' c1=',c1:6:6,' ny1=',ny1:6:6);

end;

if (P2>P[i-1]) and (P2<P[i+1]) then

begin

a2:=((ny[i+1]-ny[i-1])/(p[i+1]-p[i-1])-(ny[i]-ny[i-1])/(p[i]-p[i-1]))/(p[i+1]-p[i]);

b2:=(ny[i]-ny[i-1])/(p[i]-p[i-1])-a2*(p[i]+p[i-1]);

c2:=ny[i-1]-a2*sqr(p[i-1])-b2*p[i-1];

NY2:=a2*sqr(p2)+b2*p2+c2;

writeln(f3,' a2=',a2:6:6,' b2=',b2:6:6,' c2=',c2:6:6,' ny2=',ny2:6:6);

end;

end;

close(f1);

close(f2);

close(f3);

end.

3.6 Результаты вычислений

Ответ:

μ1=0,9740 мПа·с

μ2=0,6495 мПа·с

4 Задание №4

4.1 Условие задачи

Найти распределение давления в остановленной газовой скважине от поверхности до глубины км. При решении использовать метод Рунге-Кутта.

Использованы следующие обозначения: - зависимость плотности газа от давления, h-шаг, км; -давление на поверхности, МПа; с - коэффициент пропорциональности в дифференциальном уравнении.

h

c

0,0009

0,74

-1,7198

1

60

0,1

0,08

4.2 Основные теоретческие выкладки

Распределение давление по стволу газовой скважины определяется из следующего отношения

,

Пусть , тогда

4.3 Ручной счет

Найдем давление, Р1 на глубине 0,1 км:

K1=h∙c∙ (a0∙p02+a1∙p0+a2)=0,1∙0,08∙ (0,0009∙602+0,74∙60-1,7198)= 0,3674

K2=h∙c∙ (a0∙ (p0+K1/2)2+a1∙ (p0+K1/2)+a2)=

=0,1∙0,08∙ (0,0009∙ (60+0,3674/2)2+0,74∙ (60+0,3674/2)-1,7198)=0,3686

K3:=h∙c∙ (a0∙ (p0+K2/2)2+a1∙ (p0+K2/2)+a2)=

=0,1∙0,08∙ (0,0009∙ (60+0,3686/2)2+0,74∙ (60+0,3686/2)-1,7198)= 0,3686

K4:=h∙c∙ (a0∙ (p0+K3)2+a1∙ (p0+K3)+a2)=

= 0,1∙0,03∙ (0,0009∙ (60+0,3686)2+0,6767∙ (60+0,3686)-11,131)= 0,3699

∆P1=(K1+2∙K2+2∙K3+K4)/6=(0,3674+2∙0,3686+2∙0,3686+0,3699)=

=0,369

p1 = p0+ ∆p1=60+0,369=60,740

Аналогично считается давление на остальных глубинах:

P0

60

∆P

K1

K2

K3

K4

P1

60,369

0,369

0,3674

0,3686

0,3686

0,3699

P2

60,740

0,371

0,3699

0,3711

0,3711

0,3724

P3

61,113

0,374

0,3724

0,3736

0,3737

0,3749

P4

61,490

0,376

0,3749

0,3762

0,3762

0,3775

P5

61,868

0,379

0,3775

0,3788

0,3788

0,3801

P6

62,250

0,381

0,3801

0,3814

0,3814

0,3827

P7

62,634

0,384

0,3827

0,3840

0,3840

0,3853

P8

63,020

0,387

0,3853

0,3866

0,3866

0,3879

P9

63,410

0,389

0,3879

0,3892

0,3892

0,3906

P10

63,801

0,392

0,3906

0,3919

0,3919

0,3933