Воспроизводственные затраты оцениваются по удельным замыкающим нормативам, с учетом интегрального фактора времени. Затем определяются полные затраты.
В блоке 5 производятся вычисления доходной части критерия, которая формируется исходя из годовой и суммарной добычи газа и сопутствующих компонетов, оцениваемых по оптовым (стартовым) ценам, а в последующие годы учитывается их рост во времени в соответствии с реалистичным учетом фактора времени. Вторым слагаемым доходной части является оценка необходимой для уравнивания вариантов продукции по замыкающим затратам, исчисленная с интегральным учетом фактора времени.
В блоке 6 производится расчет максимума приведенного эффекта.
В блоке 7 организована печать результатов вычисления.
Описанная блок-схема, построенная на основе интегрального критерия, представляет собой экономическую модель, позволяющую учесть специфику разработки газовых и газоконден-сатных месторождений совместно с приходной и расходной частями перспективного ТЭБ страны.
РЕАЛИЗАЦИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ НА ЭВМ
Метод решения уравнений (28), (29) и (30) - итерационный. Метод модифицирован с целью улучшения сходимости итерационного процесса. Начальные приближения для каждого временного шага определяются по расчетной схеме В.К. Саульева. Для каждого узла сеточной области используется свой коэффициент релаксации, зависящий от характера сходимости в данном узле. Сходимость итерационного процесса проверяется по материальному балансу в каждом узле с точностью 0,05%.
Временной шаг Д* дробится на полушаги по А//2, на каждом из которых расчет ведется по разным уравнениям:
1) t * t + A//2; обход сетки, аппроксимирующей пласт, совершается в порядке возрастания индексов i = 1, imwc;
-*■"!, ■'max*
При этом обходе решаются уравнения: для газовой залежи
°>5си
['
и+1/2 г-in
// L Г Г J //
138
х
£ а.
/-1.3 '
/2 |
[ |
п+1/2*|2 Г ih-1/212 г/ J |_ ,„ J
/-2,4
а, х
(31)
для водоносного бассейна
п+1/2 п
,"+1/2 в/у
,/t+l/2 .л+1/2 л.. +л |
dl 1Рв/ "Рв// ' /-2,4 ' ( В/ BI/J
уравнение связи .. +л..
П/ В J/
, . » n+l/2
я+1/2 „п+1/2
(32)
(33)
Коэффициенты сг, Л^, af, rfL получаются при переходе от
уравнений (28), (29), (30) к их разностным аналогам (31), (32), (33).
Преобразуя уравнение (31), получаем
♦'•-»•
где
А = X а;
Bt. - 0,5c|y/ZU
X
(34) 139
и р |
Поскольку входящие в р.. величины р^ и р3
как
давления уже просчитанных узлов.^звестны, то уравнения решаются явно (допуская, что Z(p.. «■ Zip.));
В. . + V (В2 . - 4 А. .р. .)
U____________ I]___________ I Г 11
2А
(35)
Коэффициенты а.> Ь.., с, d.. имеют вид:
а
1 kh
[iz |
1dx2
/Д/
или
dr
kh \lz
*«/-
' a)muhbn.
kh
kh
kh
4/1
(amh)
LL
l dx
kh
I,U
Индекс J соответствует смежным с iif) узлам при этом / - 1 соответствует узел (Г1, />; / - 2 - (ц / V, I - о - и ,
7 ' 2)~ ^ + а//2 -►* + -♦<; обход сетки гпюизводитсяобратным путем - порядке убывания индексов: i - w, 1; / " ^max> !•
При этом решаются уравнения:
0,5 с,, |
/-2,4 |
для газовой залежи
• |
п+1 |
п+1 /2] |
||||
p/Z(p) |
p/Z(p) |
|||||
т |
и |
'/ ) |
||||
( |
1 |
2 |
А |
|||
п+1 Pi |
р' |
г+1 |
}■ |
f S -Of х |
140
n+l/2 Pi
n+1/2
- с
чЛ+1
(36)
для водоносного бассейна
Уважаемый посетитель!
Чтобы распечатать файл, скачайте его (в формате Word).
Ссылка на скачивание - внизу страницы.