Оптимизация темпов разработки газовых месторождений с учетом ресурсосбережения: Производственно-практическое издание, страница 88

Воспроизводственные затраты оцениваются по удельным за­мыкающим нормативам, с учетом интегрального фактора време­ни. Затем определяются полные затраты.

В блоке 5 производятся вычисления доходной части крите­рия, которая формируется исходя из годовой и суммарной до­бычи газа и сопутствующих компонетов, оцениваемых по опто­вым (стартовым) ценам, а в последующие годы учитывается их рост во времени в соответствии с реалистичным учетом факто­ра времени. Вторым слагаемым доходной части является оценка необходимой для уравнивания вариантов продукции по замы­кающим затратам, исчисленная с интегральным учетом фактора времени.

В блоке 6 производится расчет максимума приведенного эф­фекта.

В блоке 7 организована печать результатов вычисления.

Описанная блок-схема, построенная на основе интегрально­го критерия, представляет собой экономическую модель, поз­воляющую учесть специфику разработки газовых и газоконден-сатных месторождений совместно с приходной и расходной частями перспективного ТЭБ страны.

РЕАЛИЗАЦИЯ МАТЕМАТИЧЕСКОЙ МОДЕЛИ НА ЭВМ

Метод решения уравнений (28), (29) и (30) - итерационный. Метод модифицирован с целью улучшения сходимости итераци­онного процесса. Начальные приближения для каждого времен­ного шага определяются по расчетной схеме В.К. Саульева. Для каждого узла сеточной области используется свой коэф­фициент релаксации, зависящий от характера сходимости в данном узле. Сходимость итерационного процесса проверяется по материальному балансу в каждом узле с точностью 0,05%.

Временной шаг Д* дробится на полушаги по А//2, на каж­дом из которых расчет ведется по разным уравнениям:

1) t * t + A//2; обход сетки, аппроксимирующей пласт, совершается  в  порядке  возрастания  индексов   =   1,   imwc;

-*■"!, ■'max*

При этом обходе решаются уравнения: для газовой залежи


°>и


['


и+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


(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)


для водоносного бассейна