УДК 517.97
Принцип максимума в задаче максимизации дохода для модели газового месторождения
А. К. Скиба
Вычислительный центр имени А. А. Дородницына Российской академии наук ул. Вавилова, 40, Москва, 119333, Россия
Рассматривается модель функционирования газового месторождения. Ставится и решается на бесконечном интервале задача оптимального управления. В описание критерия включена стоимость реализации газа, заданная вогнутой функцией.
Ключевые слова: принцип максимума; модель газового месторождения; бесконечный интервал; фазовые диаграммы.
1. Введение
Рассмотрим модель функционирования газового месторождения с взаимовли-яющими скважинами [1—4]. Введём следующие обозначения: N — фонд добывающих скважин, IV — ограничение сверху на фонд добывающих скважин, д — средний дебит добывающих скважин, ё — коэффициент дисконтирования, V — извлекаемый запас, с — стоимость природного газа.
Введём понятие закупочной функции и (Я), которая устанавливает зависимость между стоимостью газа и его объёмом. Закупочная функция обладает следующими свойствами. Это вогнутая строго возрастающая функция, определённая на положительной полуоси. В дальнейшем мы будем пользоваться функциями из известного класса вогнутых функций
и (Я) = сОР, 0 < ц< 1. (1)
Здесь и далее с> 0 и Я > 0.
Между описанными выше переменными устанавливается взаимосвязь, которую запишем в виде системы дифференциальных уравнений
V = -Мд, д = - Ю Мд (2)
Уо
при начальных условиях Уо > 0, д0 > 0. На фонд добывающих скважин накладывается следующее ограничение 0 ^ N ^ N. Здесь N > 0. Добычу газа и доход вычисляем по формулам:
сю
я = дМ, I и(Я)е-8г(И.
Проведём следующую замену переменных:
М' = ^ К, N' = ^ N, V' = |0 V = д'. (3)
Уо У о Уо
В новых переменных д' = V'. При написании указанных переменных штрихи при них в дальнейшем опускаем.
Для иллюстрации возможностей модели рассмотрим две стратегии разработки газового месторождения. При этом не учитывая ограничения на фонд добывающих скважин.
Статья поступила в редакцию 23 сентября 2008 г.
Стратегия 1. Пусть газовое месторождение разбуривается с постоянным темпом N = п на периоде [0, ¿1]. Далее разбуривание месторождения прекращается, и добыча газа осуществляется при постоянном фонде N. В этом случае величина добычи газа изменяется по следующему закону:
= Г (I + Ш)Яо ехр(-1о£ - Ш2), при £ < ^,
= \ (N0 + Ш 1)qо ехр(-МоЬ -Ш2 - (1о + Ш 1)^ - Ь)), при г> и.
Анализ функции (4) показывает, что на начальном этапе разбуривания месторождения добыча газа на месторождении возрастает. Это связано тем, что прирост добычи за счёт ввода новых скважин превышает темп падения добычи, связанный с извлечением запасов газа. Если на этом этапе происходит прекращение добычи газа, то в момент ¿1 добыча Я^) достигает максимума, и на графике функции Я^) наблюдается характерный излом. Далее добыча падает. Если продолжать разбуривание месторождения с постоянным темпом достаточно долго, то прирост добычи газа за счёт ввода новых скважин будет уменьшаться, и со временем он совпадёт с темпом падения добычи. В этот момент величина добычи ( ( ) достигает максимума. Далее добыча газа падает.
Такие характерные особенности поведения добычи газа наблюдались на практике при разбуривании газовых месторождений.
Стратегия 2. Необходимо в динамике определить объем действующего фонда скважин, обеспечивающий извлечение всего запаса газа при постоянном заранее заданном уровне его добычи ((. Интерес в постоянном уровне добычи газа проявляется как со стороны промысловиков, так и со стороны потребителей газа. Промысловикам постоянный уровень добычи необходим для закупки и настройки промыслового и транспортного оборудования под обеспечение данного уровня добычи газа. Потребителям газа постоянный уровень добычи необходим при решении вопроса о стабилизации объёмов закупаемого газа. Приведённые выше доводы делают интересными исследование данной стратегии разработки месторождения. Перейдём к формальному описанию и исследованию упомянутой стратегии. Объем действующего фонда скважин изменяется по следующему закону:
N (€) = —^. (5)
о - (п
Таким образом, функция (5) определена на полуинтервале [0,£*), где ¿* = %/((. В конечный момент £ = ¿* происходит полное извлечение запасов газа (^о = Уо). В то же время Иш N(р) = <Х). Значит, сетка скважин уплотняется
до бесконечности. Такую стратегию невозможно полностью реализовать на практике. Однако она реализована частично. Этапу постоянной добычи газа предшествует этап нарастающей добычи (для небольших газовых месторождений он может отсутствовать) и следует этап падающей добычи. Проблема увеличения периода постоянной добычи особо актуальна в настоящее время, когда многие газовые месторождения завершают период постоянной добычи.
Далее ставим задачу поиска максимума функционала на бесконечном интервале. В качестве функционала берётся подверженный дисконтированию накопленный доход.
Задача 1. Требуется максимизировать функционал
сю
с J дч^е(6)
о
при дифференциальной связи
4= -Я = -1N, (7)
начальном условии
о>0
(8)
и ограничении на управление
0 < N < N. (9)
Управление N принадлежит множеству измеримых функций. Правый конец свободен.
Прежде чем приступить к поиску решения задачи 1, необходимо ответить на вопрос о его существовании. Для этого воспользуемся теоремой 3.6 статьи [5]. В связи с труднодоступностью данной статьи приведём в следующем параграфе общую постановку задачи оптимального управления на бесконечном интервале, основные определения и формулировку теоремы существования.
2. Общая постановка задачи оптимального управления на бесконечном интервале и теорема
существования
Рассмотрим общую постановку задачи оптимального управления
inf I (х,и).
Здесь П обозначает множество допустимых пар (х,и), х G АСЦОс([0,œ)) и и : [0, <Х)) ^ Мг является С-измеримой функцией такой, что почти для каждого t
x(t) G A(t), u(t) G U(t, x(t)), x = f (t, x(t), u(t)).
Здесь A обозначает множество многозначных функций из [0, œ) в с С х Вт-измеримым графиком A; U обозначает множество многозначных функций из А в Rm. График M = {(t,x,u) : t G [0, œ),x G A(t), и G U (t, x)} является С х Вт+Г -измеримым. Также f : M ^ является С х Вт+Г измеримой функцией. Наконец, I задаётся
сю
I(х,и) = J f0(t,x(t),u(t))dt, (х,и) G П,
0
где fo : M ^ (-œ, œ) является С х Вт+Г-измеримой функцией. Определим для (t, x) G А многозначную функцию
Q(t, x) = {(z0, z) : z0 > f0(t, x,u), z = f (t, x,u), и G U (t, x)}.
Определение. Говорят, что многозначная функция Q из А в Rm+1 обладает свойством (К) в (t, x) G А, если
Q(t, х)= Q cl {UQ(t, x!) : x! G A(t), lx' - xl < .
¿>0
В литературе это свойство также называют свойством (U). Заметим, что свойство (К) в (t,x) G А означает замкнутость множества Q(t,x).
Наконец, определим множество Па = {(х,и) G П : I(х,и) ^ а} для какого-нибудь a G R.
Перейдём к формулировке теоремы существования.
Теорема. Пусть для каждого t > 0 f (t, ■, ■) непрерывна на M (t), f0(t, ■, ■) полунепрерывна снизу на M (t), A(t) замкнуто, M (t) замкнуто. Предположим также, что существует a G R такое, что {ж(0) : (х,и) G Па} ограничено, {f (■,&(■), Ц^))|[0, T] : (х,и G Па)} равномерно интегрируемо для каждого Т > 0, Q обладает свойством (К) в каждой (t, x) G A, Q(t, х) выпукло в каждой (t,x) G А, {/-(■, х(^),и(^)) : (х,и) G Па} является сильно равномерно интегрируемой, Па непусто.
Тогда существует допустимая пара (x,ù) G П такая, что
I(х,и) = inf I(х,и).
3. Исследование задачи
Приведённая в предыдущем параграфе теорема является теоремой существования решения для задачи оптимального управления на бесконечном интервале.
Рассмотрим множество
Ü(t, q) = {(z°, z) : > -( qN)*e-5t, z= -qN, N e [0,N]} .
Множество Q(t, q) = {(z°, z) : z° ^ — ( —z)ve-St} является выпуклым множеством для каждого (t, q). Для исследуемой задачи также справедливо выполнение всех остальных условий теоремы. Таким образом, существует такое управление N(t) e [0, N], которое при дифференциальной связи (7) и начальном условии (8) доставляет максимум функционала (6).
Для поиска решения задачи будем использовать принцип максимума Понтря-гина [6,7] со свободным правым концом. Управление N, удовлетворяющее необходимым условиям оптимальности, ищем на множестве измеримых функций. Выпишем гамильтониан, сопряжённое уравнение и условие трансверсальности в формах, используемых в работе К. Эрроу [8]:
Н = c(q N — ^qN, (10)
ip = ^5 — сщп-1№п + (11)
lim ip(t)q(t)e-st = 0. (12)
t—
Здесь параметр r] e (0,1).
Следует отметить, что гамильтониан, сопряжённая переменная и условие трансверсальности в классическом представлении отличаются в их написании от гамильтониана, сопряжённой переменной и условия трансверсальности, используемых в настоящей работе и в работах [8,9], на величину e-St. Использование в таком виде переменных облегчает решение задач.
Остановимся подробнее на условиях трансверсальности, связанных с ними проблемах и способами их решения. Для решения задач оптимального управления на бесконечном интервале со свободным правым концом без достаточно строгой обоснованности был принят общий вид условий трансверсальности в усиленной форме [8]. До сих пор неизвестны в общем случае достаточно строгие доказательства условий трансверсальности на бесконечном интервале. Однако в работе С.М. Асеева и др. [10] для одной задачи оптимального управления на бесконечном интервале они были получены. Кроме того, нам также неизвестны примеры, противоречащие условиям трансверсальности в усиленной форме.
Вместе с тем в последнее время к задачам оптимального управления на бесконечном интервале проявляется все больший интерес. Это связано не только с тем, что уменьшается количество вводимых параметров при постановке задачи, но также с тем, что в динамических моделях при долгосрочном планировании не всегда можно строго точно указать время окончания планирования. Кроме того, при решении задач на бесконечном интервале легче определить и проанализировать магистральные свойства модели. Для решения этой и других аналогичных задач можно предложить несколько способов выхода из создавшейся ситуации:
— несмотря на отсутствие в общем случае достаточно строгих доказательств условий трансверсальности на бесконечном интервале времени, применить его формально с целью поиска решения задачи оптимального управления. Как правило, научные работники, сталкиваясь со сходной проблемой, до сих пор использовали именно данный подход;
— полностью отказаться от применения условия трансверсальности на бесконечном интервале времени, и поиск оптимального решения осуществить на множестве всех траекторий, удовлетворяющих необходимым условиям оптимальности. Далее сравнить значения функционалов с целью поиска наибольшего значения. Иногда удобно предварительно произвести исследования поведения траекторий на фазовой диаграмме;
— отказаться от поиска решения задачи оптимального управления на бесконечном интервале со свободным правым концом и заменить её другой постановкой задачи, максимально приближенной к исходной, но уже со строго доказанными условиями трансверсальности. Это возможно сделать, заменив бесконечный интервал конечным и изменив целевую функцию.
Мы воспользовались первым способом. Перейдём к обсуждению основных сложностей, связанных с поиском оптимального решения. Принцип максимума даёт нам необходимые условия оптимальности. Чтобы убедиться, что данная траектория действительно является оптимальной, следует применить достаточные условия оптимальности. Если для рассматриваемой задачи достаточные условия не доказаны, то имеется иной подход, а именно: необходимо доказать существование оптимального решения и единственность траектории, удовлетворяющей всем необходимым условиям оптимальности. Как правило, в большинстве практических задач оптимальное решение существует. Единственность иногда вытекает из самой аналитической процедуры поиска траектории, удовлетворяющей всем необходимым условиям оптимальности. Часто бывает, что мы угадываем такую траекторию, но это не даёт гарантий отсутствия других траекторий, также удовлетворяющих всем этим условиям. Поэтому необходимо рассмотреть остальные понтрягинские траектории и убедиться, что среди них отсутствуют траектории, удовлетворяющие условиям трансверсальности. С помощью математических выкладок сделать это достаточно сложно. Однако процедура значительно облегчается, если её осуществить на фазовой диаграмме. Если получены несколько траекторий, удовлетворяющих необходимым условиям оптимальности, то остаётся только вариант сравнения значений функционалов на этих траекториях.
Продолжим поиск оптимального решения для задачи 1. Для этого формально применим принцип максимума к задаче 1 без учёта ограничения (9).
Продифференцируем гамильтониан (10) по N и полученную производную приравняем к нулю:
H'N = cr]qvNV-1 — фq = 0.
Отсюда получаем
ф = crj(q N У-1. (13)
С учётом соотношения (13) разрешим сопряжённое уравнение (11). В Результате получаем
ф = ф0 (14)
и
ф = ф° est. (15)
Из условия трансверсальности (12) с учётом (15) вытекает
lim ( ) = 0. (16)
i—
Используя дифференциальные уравнения (14) и (7), с учётом (13) приходим к следующему дифференциальному уравнению:
? = — Ч-г = (СП)-^ ф ^+1. (17)
dq qN v " r v 7
С учётом (15) и (16) проинтегрируем дифференциальное уравнение (17) от t до ж:
1 . „ 1
ф = (crj) -q. (18)
1
Из (18) с учётом (13) получаем
с
т = —. (19)
Продолжим решение задачи 1 с учётом ограничения (9). Пусть выполнено неравенство
с
-г^— < N. (20)
1 — г]
Покажем, что максимум гамильтониана (10) на траектории с управлением (19) при всех значениях £ достигается внутри отрезка [0, Л]. Для этого, воспользовавшись (18), (19) и дифференциальным уравнением (7), после несложных преобразований получаем
НШ),ф(*),Й(г)) = тах [НШ= тах [сдп!п — фд • I] =
N е[0,лт] N е[0,лт]
тах
N е[0,м]
■п - ( Ь \П 1 „ _ 4.
сд'^е г — — сг] -^ д'^е 1
V - т^ t
= сд0е тах
N е[0,лг]
/ <5 V-1
Отсюда и из (20) вытекает, что при всех £ £ (0, то) максимум гамильтониана достигается внутри отрезка [0, IV]. Легко показать, что условие трансверсальности (12) выполнено. Для того, чтобы траектория с управлением (19) была оптимальной, необходимо убедиться в отсутствии других траекторий, удовлетворяющих необходимым условиям оптимальности. Рассмотрим траекторию с начальной точкой ( сг}(]о>1'п-1, д0) и с управлением N = N. Воспользовавшись (11) и (7), получаем
фд = сщт11-1 [(5 + ф — IV) + IV] /(ё + г]1).
Отсюда с учётом (20) вытекает следующее неравенство:
фд < сщ^И*п-1.
Значит максимум гамильтониана (10) достигается при N = N. Таким образом, если при некотором т управление И(т) = N, то N(1) = N при £ ^ т. Легко проверить, что условие трансверсальности (12) в этом случае не выполняется.
Переходим ко второй части доказательства. Пусть
с
N < 1-. (21)
1—
Покажем, что управление N = N удовлетворяет необходимым условиям оптимальности на интервале (0, то). Умножим дифференциальное уравнение (11) на дехр(—61) и после несложных преобразований приходим к следующему дифференциальному уравнению:
d ( • -81
('фде= —спЫЮ* е-н. (22)
Принимая во внимание условие трансверсальности (12), дифференциальное уравнение (7) и полагая N = N, проинтегрируем обе части дифференциального уравнения (22) от £ до то:
сю _
Н = сп I (ф) • = У+Ш. (23)
Отсюда
ф = С^Vд0- ^^к (24)
5 + Г] N х '
Подставим функцию (18) в гамильтониан (10):
Н {тШ *) = " - ^я) = сдЦ " - ^ ■ е-"«
Вычислим первую производную гамильтониана в точке N = N
Н'м) = сдЦе4^^ (5 - (1 - ^) .
Отсюда и из предположения (21) вытекает положительность производной на правом конце отрезка [0, И]. Кроме того, гамильтониан является вогнутой возрастающей функцией от N. Значит, при управлении N = N достигается максимум гамильтониана.
Доказательство единственности решения задачи 1 при выполнении неравенства
> IV (25)
1-
проведём на плоскости (фд, д).
На рис. 1 изображена кривая С, которая описывается уравнением
фд = сг] д*!*-1. (26)
Фя Траектория I Кривая О
Траектория III (оптим. тр.) Траектория II
* у' / / / /' я
Рис. 1. Фазовая диаграмма
Данная кривая делит фазовую плоскость на две части. Траектории «подозрительные» на оптимальность описываются дифференциальным уравнением
^ = - ^ + о?д^-1. (27)
Ад д1
Из максимума гамильтониана (10) вытекает, что если траектория, описываемая дифференциальным уравнением (27) лежит выше кривой С, то управление N < N, в противном случае N = N.
Пусть траектория, описываемая дифференциальным (27), лежит ниже кривой С, тогда, разрешив дифференциальное уравнение (27), получим
фд = В д-ъ + сЙ^-1 д ч , (28)
r]N + 5 х '
где
В =
ъ* - ^
Чо . (29)
Пусть выполнено неравенство (25), тогда траектория, описываемая функцией (23), лежит ниже кривой С (см. рис. 1). Возьмём на кривой С при д > 0 произвольную точку (фдо, д0) и определим поведения траекторий в точке (фд0, д0) кривой С. Для этого определим угол наклона траекторий, лежащих выше и ниже кривой С, а также угол наклона самой кривой С. В результате получаем
А(фд) 5(фдГ^ + ффд = с7д0-1й,-1(- V + 1), (30)
d q ( er] q) — q N
= -5сг]q--1NV~2 + cqq--1N--1 = сщ--1 NV--1 N + . (31)
^ = стГяГ1^- (32)
С учётом неравенства (25) вытекает, что угол наклона траекторий, лежащих как выше, так и ниже кривой G, меньше чем угол наклона самой кривой G. Значит, в кривую G входят траектории, лежащие ниже кривой G, и из неё исходят траектории, лежащие выше кривой G. Получен достаточно важный факт — траектория, оказавшаяся в области выше кривой G, никогда в дальнейшем её не покидает. Таким образом, если в некоторый момент т управление N(т) < N, то N(t) < IN при всех t > т. Из этого факта вытекает, что если траектория, удовлетворяющая всем необходимым условиям оптимальности, существует в области выше кривой G, то N(t) = j^-. Данное управление является недопустимым, так
как нарушается ограничение (25). Таким образом, траектории (траектории I), достигающие кривую G, не удовлетворяют необходимым условиям оптимальности.
Рассмотрим траекторию, начало которой лежит ниже кривой G. В описании траекторий, лежащих ниже кривой G, постоянное интегрирования В в (28) принимает как отрицательные, так и положительные значения. Рассмотрим случай В = 0, тогда (28) перепишется в виде
н = cN--1q-^. (33)
Данная траектория (траектория III) не достигает кривой G и удовлетворяет условию трансверсальности. Легко убедиться, что при отрицательных значениях постоянной интегрирования В траектория (28) также не достигает кривой G и условию трансверсальности не удовлетворяет. На рис. 1 эта траектория изображена как траектория II. При положительных значениях постоянной интегрирования В траектория (28) достигает кривой G при некотором положительном q. Как было показано выше, она не удовлетворяет необходимым условиям оптимальности.
Из вышесказанного вытекает существование единственной сопряжённой переменной (в форме, предложенной К. Эрроу)
m = cN--1q---1 -е(1--Ж (34)
и единственного управления N = N, удовлетворяющих всем необходимым условиям оптимальности. Из последнего предложения и из существования оптимальной траектории для задачи 1 вытекает то, что траектория III является оптимальной.
Для задачи 1 справедлива следующая теорема.
Теорема. Если выполнено неравенство < N, то оптимальным решением задачи 1 при г] е (0,1) на всем интервале (0, то) является управление,
заданное формулой N(t) = -—j. В противном случае оптимальным решением
является управление N = N.
Литература
1. Федосеев А. В., Хачатуров В. Р. Постановка и исследование задач оптимального управления для анализа перспективных планов в нефтегазодобывающей промышленности // Имитационное моделирование и матем. методы анализа перспективных планов развития нефтедобывающей промышденности. — М.: ВЦ АН СССР, 1984. — С. 66-112.
2. Маргулов Р. Д., Хачатуров В. Р., Федосеев А. В. Системный анализ в перспективном планировании добычи газа. — М.: Недра, 1991.
3. Крюков В. А., Скиба А. К., Федосеев А. В. Задачи оптимального управления разработкой газоконденсатного месторождения. — М.: ВЦ АН СССР, 1990.
4. Моделирование освоения газовых месторождений на заключительной стадии эксплуатации / А. К. Скиба, В. Хачатуров, А. В. Злотов, А. Н. Соломатин. — ВЦ РАН, 2006.
5. Balder E. J. On Existence Result for Optimal Economic Growth Problems // Journal of Mathematical Analysis and Applications. — 1983. — Т. 95. — С. 195213.
6. Флеминг У., Ришел Р. Оптимальное управление детерминированными и стохастическими системами. — М.: Мир, 1978.
7. Моисеев Н. Н. Элементы теории оптимальных систем. — М.: Наука, 1975.
8. Эрроу К. Применение теории управления к экономическому росту. Матем. экономика. — М.: Мир, 1974.
9. Skiba A. K. Optimal Growth with a Convex-Concave Production Function // Econometrica. — 1978. — Т. 46, № 3. — С. 527-539.
10. Асеев С. М., Кряжимский А. В., Тарасьев А. М. Принцип максимума Понт-рягина и условия трансверсальности для одной задачи оптимального управления на бесконечном интервале // Тр. МИ РАН. — 2001. — Т. 233. — С. 71-88.
UDC 517.97
Maximum Principle in a Problem of Maximization of the Income for Model of a Gas Deposit
A. K. Skiba
Dorodnicyn Computing Centre of the Russian Academy of Sciences Vavilov st. 40, 119333 Moscow, Russia
The model of functioning of a gas deposit is considered. The optimum control problem is put and dares over an Infinite Horizon. Cost of realisation of the gas is included in the criterion description, set by concave function.