- •По дисцеплине «инженерные методы расчета в нефтяной и газовой промышленности»
- •Содержание
- •2. Задание №2…………………………………………………………………10
- •3. Задание №3…………………………………………………………………16
- •4. Задание №4…………………………………………………………………22
- •1 Задание № 1
- •1.1 Условие задачи
- •1.2 Основные теоретические выкладки
- •1.3 Ручной счёт
- •1.4 Блок - схема алгоритма решения задачи
- •1.5 Текст программы
- •1.6 Результаты вычислений
- •1.7 График
- •2 Задание № 2
- •2.1 Условие задачи
- •2.2 Основные теоретические выкладки
- •2.3 Ручной счет
- •2.4 Блок - схема алгоритма решения задачи
- •2.5 Текст программы
- •3.2 Основные Теоретические выкладки
- •3.3 Ручной счет
- •3.4 Блок-схема алгоритма решения задачи
- •3.5 Текст программы
- •3.6 Результаты вычислений
- •4 Задание №4
- •4.1 Условие задачи
- •4.2 Основные теоретческие выкладки
- •4.3 Ручной счет
- •4.4 Блок-схема алгоритма решения задачи
- •4.5 Текст программы
- •4.6 Результаты вычисления
- •4.7 График
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 |