Математическое моделирование
УДК 517.97
Смешанное ограничение в прикладной задаче оптимального управления
А. К. Скиба
Вычислительный центр им. А.А. Дородницына РАН ул. Вавилова, 40, Москва, Россия, 119333
В статье формулируется модель функционирования группы газовых месторождений и решается задача оптимального управления со смешанным ограничением. В качестве оптимизационного критерия задачи выбран максимум накопленной добычи для группы газовых месторождений.
Ключевые слова: принцип максимума, предложение К. Эрроу, смешанное ограничение, модель газовых месторождений, максимум накопленной добычи, максимум длины «полки».
1. Введение
В 1968 году вышла статья лауреата Нобелевской премии К. Эрроу [1], сыгравшая значительную роль в дальнейшем продвижении использования оптимального управления в теории экономического роста. В работе подчёркивалась особая важность применения оптимального управления в экономической теории. На основе принципа максимума Понтрягина были без достаточно строгого доказательства сформулированы предложения, являющиеся основой для решения широкого класса задач оптимального управления.
К. Эрроу продемонстрировал использование своих предложений на примере решения некоторых задач теории экономического роста. Полученным результатам была дана экономическая интерпретация.
Для решения задач оптимального управления со смешанными ограничениями К. Эрроу сформулировал предложения, являющиеся модификацией принципа максимума Понтрягина. В предложения дополнительно включены элементы нелинейного программирования — лагранжиан, множители Лагранжа и условия дополняющей нежёсткости. Также в упомянутой статье уделялось большое внимание конечному и бесконечному интервалу времени, теоремам существования, необходимым и достаточным условиям оптимальности.
После выхода в печать статьи К. Эрроу и её широкого обсуждения в научных кругах обозначились, по крайней мере, два направления в дальнейших исследованиях. С одной стороны, несмотря на не вполне строгие доказательства предложений К. Эрроу, использовать их для дальнейших научных работ не только в теории экономического роста, но и в других прикладных дисциплинах. Например, основываясь на предложениях К. Эрроу с успехом было проведено исследование односекторной модели на бесконечном интервале времени с выпукло-вогнутой производственной функцией [2]. Были получены достаточно интересные результаты.
С другой стороны, акцентировать большее внимание на создание строгой математической базы предложений К. Эрроу. В 1983 году вышла работа Е. Бол-дера [3], где была доказана теорема существования на бесконечном интервале времени. Значительный прорыв в формировании теоретических основ был совершён в 2007 году С.М. Асеевым и А.В. Кряжимским в их совместной публикации монографии «Принцип максимума Понтрягина и задачи оптимального экономического роста» [4]. Данная монография посвящалась теории принципа максимума Понтрягина для определённого класса задач оптимального управления с бесконечным интервалом времени. Были получены многие уникальные результаты,
Статья поступила в редакцию 17 мая 2012 г.
среди которых особо хотелось выделить исследование поведения сопряжённых переменных и гамильтониана на бесконечности.
Прежде чем перейти к описанию модели, сформулируем предложение К. Эр-роу, которое в дальнейшем будем использовать в качестве основного аппарата оптимального управления при исследовании задачи максимизации накопленной добычи.
2. Формулировка предложения К. Эрроу
В предложении К. Эрроу используются понятия: вспомогательные переменные и выбор инструментов. Такие словосочетания соответствуют стандартным понятиям, принятыми в теории оптимального управления, а именно: сопряжённым переменным и управлениям.
Для лучшего восприятия предложения К. Эрроу будем пользоваться стандартными понятиями, принятыми в теории оптимального управления.
Предложение 1. Пусть у(Ь) — управления (0 ^ Ь ^ Т), максимизирующие функционал
т
I и [х(г),у(г),г]<и
0
при условии
(a) х = с[х(г),у(г),г], ограничениях на управления
(b) ^[х(г),у(г),г] > о,
которые, возможно, включают переменные состояния, начальные условия на переменные состояния и граничные условия
х(Т) > о.
Если условие регулярности выполнено, то существуют сопряжённые переменные р(1) такие, что для каждого момента t
(c) ь(1) максимизирует Н[х(1),ь,р(1),Ц относительно ограничений (Ь), где Н(х, у,р,1) = и(х,ь,Ь) + рО(х, V, I);
(д) р = —дЬ/дхг при х = х(£), у(£) = у(£), р = р(1), где
(е) Ь(х, у,р, д, £) = Н(х, у,р, £) + дР(х, V, £) и множители Лагранжа д такие, что
(/) дь/дУк = о при х = х(г), у'(г) = ь(г), р = р(г), д > о, дР[х(г),д(г),г] = о,
и
(д) р(Т) > о, р(Т)х(Т) = о.
В предложении вектор х(1) имеет размерность п, а вектор у(1) — размерность т.
Перейдём к описанию модели.
3. Описание модели и постановка задачи
В 2009 году в журнале Вестник РУДН была опубликована статья [5], в которой описывалась рассматриваемая в настоящей работе модель. Несмотря на это,
заново опишем её, преследуя при этом две цели. Во-первых, все же существуют некоторые отличия в описаниях моделей. Во-вторых, желательно сохранить целостность описания проведённых исследований.
Рассмотрим модель функционирования группы газовых месторождений с вза-имовлияющими скважинами [6-9]. Группа состоит из п месторождений. Введём следующие обозначения: Нг — фонд добывающих скважин на г-м месторождении, 1Чг — ограничение сверху на фонд добывающих скважин на г-м месторождении, дг — средний дебит добывающих скважин на г-м месторождении, Уг — извлекаемый запас на г-м месторождении.
Между описанными переменными устанавливается дифференциальная взаимосвязь, которую запишем в виде системы дифференциальных уравнений
И = -Мгдг, ъ = -д^)т(г), г = 1,2,... ,п,
*г
при начальных условиях > 0, д° > 0. На фонд добывающих скважин накладывается следующее ограничение 0 ^ Мг ^ . Здесь Л^ > 0. Предполагаем также
д° _
различными между собой величины —г° Л^.
уг
Имеется общее ограничение на общую «полку» месторождений
п
г=1
Накопленную добычу газа вычисляем по формуле:
т
п
г=1 0
Задача максимизации накопленной добычи является актуальной как с практической, так и с теоретической точки зрения, о чем неоднократно подчёркивалось во многих работах, в том числе в работе [10]. Перейдём к математическому описанию задачи оптимального управления со смешанными ограничениями.
Задача 1. О максимизации накопленной добычи для группы газовых месторождений
Требуется максимизировать функционал
т
п
^ (1)
г=1 0
при дифференциальных связях
п°
4г = -I = 1,2, . . . ,П, (2)
* г
начальных условиях
д° > 0, г = 1, 2,...,п (3)
и ограничениях на управления
0 < Нг < Йг, г = 1, 2,...,п. (4)
£ ФЩ^) < я. (5)
г= 1
Управления Нг(г = 1,2,...,п) принадлежат множеству измеримых функций. Правый конец фазовой траектории свободен. Между собой различны величины
а° -
V01 г.
Существование оптимального управления следует, например,из теоремы, приведённой в [11, § 4.2].
Перейдём к исследованию задачи 1.
4. Исследование задачи
Сначала введём обозначения = 770 Для всех i = 1,2,... ,п. Далее упоря-
„о _ Я
у?
дочим фазовые переменные в порядке возрастания величин а0Кг. Заметим, что такой порядок единственный, поскольку согласно условиям в постановке задачи различны величины а0Мг.
Согласно предложению 1 выпишем гамильтониан и лагранжиан:
п
Н(д, М, -ф) = ^[дгМг - -ф^Щ, (6)
г=1
L(q, N, -ф, -у, 5, ¡3) = ^[giNi - фга0гдгМг + Ъ(Пг - Мг)] + 5гЩ + ¡3
г=1
Q - Y^tiiNi)
i=1
Далее опишем множество С допустимых управлений, на котором осуществляется максимизация гамильтониана (6):
G =j N е Rn | 0 < N < N, ^ giNi < . (7)
Основная цель исследования заключается в поиске на основании предложения 1 вектор-функций ф(t) и управлений N(t), которые удовлетворяют ниже приведённой сопряжённой системе дифференциальных уравнений (10), условиям трансверсальности
ф(Т )q(T) = 0, (8)
и при каждом t е [0,Т] управления N(t) максимизируют гамильтониан
n
н(q(t),N(t),m) = maxHm,N,m) = max- «оштщ. (9)
iv G^ i g^
г= 1
Система сопряжённых уравнений представится в виде:
*k(t) = -dL/dgi = -Ni(t) + a0ipi(t)Ni (t) + P(t)Nt(t), i = 1, 2,...,n или в более удобной форме
ijji(t) = [a0Mt) + p (t) - 1]Ni(t), i = 1, 2,... ,n. (10)
Положим в системе сопряжённых уравнениях (10)
= а° щг), г = 1, 2,...,п, (11)
Ьг(г) = [т - 1]ЩЬ), г = 1,2,...,п, (12)
тогда система сопряжённых уравнений (10) перепишется в следующем виде:
г'Ш = (ц(€)ф1(€) + Ь(€), ъ = 1, 2,...,п. (13)
Согласно предложению 1 множители Лагранжа /3), 6(1) должны удовлетворять равенствам:
дЬ/дКг = дг(г) - а%фг®Ш - тш - + 5г(г) = 0, г = 1,2,.
, п
или
дг(Ф - - т] = Ъ® - Ш, г = 1,2,...,п;
р (г)
я - ^ тт)
г= 1
= 0, р(г) > 0;
ъ(Шг - Йг ф] = 0, Ъ(г) > 0, г = 1, 2,...,п; бг(г)Щг) = 0, 5г(Ь) > 0, г = 1,2,...,п.
(14)
(15)
(16) (17)
В нелинейном программировании выписанные выше соотношения (15), (16), (17) определяются как условия дополняющей нежёсткости.
Как видно из дифференциальных уравнений (2), фазовые переменные д(Т) положительны при любых конечных значениях Т. Поэтому условия трансверсальности (8) представятся в следующем виде:
фг(Т) = 0, г = 1,2,...,п.
(18)
Разрешив сопряжённые уравнения (13) с учётом условий трансверсальности (18), в результате получаем
Фг^) = - ехр
1
- I аг($)М
г( ) ехр
( >, г = 1, 2,... ,п,
где аг(1) и Ъг(1) определяются формулами (11) и (12), соответственно.
Пусть д(1) — оптимальная траектория. Далее определим два множества:
н =\ <<1, 0 < i < т
{п
=1
{п
чЦ
=1
;
м = М V тт) = Я, 0 < г < Т\ .
(19)
(20) (21)
Справедливы следующие утверждения.
Утверждение 1. Пусть на оптимальной траектории при некотором значении I = ¿* £ [0, Т] выполняется неравенство
Я - ^ ШЩ^ > 0,
(22)
т.е. множество (20) непусто, тогда значение вектора оптимальных управлений определяется равенством N(£*) = N.
Доказательство. Пусть для к-й компоненты вектора управлений N(^ выполняется строгое неравенство: N0(£*) < Л^. Из (15) и (22) при Ь = £* вытекает, что = 0. С учётом (14), (16) и (17) получаем 5к(г*) > 0,^к(г*) = 0 и 1 < акфк(П.
Перепишем к-е сопряжённое уравнение (10)
оРк'Фк(€) = акМк (Макфк (V - 1] + а0^ (*)№. (23)
Положим х(г) = а0кфк(г) - 1; ф) = а0кЙк(¿); с1(?) = а0кЙк(г)/3(£). В новых переменных сопряжённое уравнение (23) перепишется в следующем виде: х = с(1)х+с1(1). Разрешим это сопряжённое уравнение с начальным условием х(Ъ*) ^ 0. В результате получим
х(Ь) = < х(Ь*) + I ё(г)ехр
- [ с(§)ё§ ёО ехр
)
(24)
Из (24) с учётом неотрицательности коэффициентов ф) и получаем х(1) = ак'Фк(¿) — 1 ^ х(Ъ*) > 0. Отсюда следует, что условие трансверсальности (18) для к компоненты вектор-функции ф(1) не выполнено. □
Следствие 1. На множестве (20) значение вектора оптимальных управлений однозначно определяется равенством N(р) = IV.
Утверждение 2. Пусть Г* £ [0,Т] и Г* £ Н, тогда [Г*,Т] С Н.
Доказательство. Предположим, что при некотором значении ^ £ (Р*,Т] выполняется равенство
п
$ - ^ ШЩг) = 0, (25)
=1
т.е. ^ £ М, тогда существует хотя бы одна к-я компонента вектора управлений N(£) такая, что Йк(Ь*) < Йк. Из утверждения 1 вытекает Ко(£*) = Йк. Пришли к противоречию. □
Утверждение 3. Пусть при некотором значении Ь = £* £ [0, Т] выполняется неравенство
п
Я - ^ ШЙг^) > 0, (26)
=1
тогда N(£) = N при всех Ь £ [£*,Т] и £* £ Н, где Н замыкание множества Н.
Доказательство. В связи с убыванием фазовых переменных и с учётом неравенства (28) связующее ограничение (5) неактивно на любом допустимом управлении при £ > ¿*. В этом случае задача 1 разбивается на п независимых оптимизационных задач, решение которых тривиально и определяется значениями вектора оптимальных управлений N(Ь) = N. Принадлежность I* замыканию множества (20) очевидна. □
Утверждение 4. Пусть Г* £ [0, Т] и Г* £ М, тогда [0, Г*] С М.
Доказательство утверждения 4 вытекает из утверждения 2 при использовании, например, метода «от противного».
Утверждение 5. Существует такое значение Т, при котором множество (20), построенное на оптимальной траектории, непусто.
п _
Доказательство. В качестве Т достаточно взять Т' = V°)/Q. Легко
г= 1
показать,что существует такое t £ (0,Т'), при котором неравенство (22) будет выполняться. Значит множество (20) непусто. □
Утверждение 6. Пусть хотя бы для одной к-й компоненты оптимального вектора управлений N(t) в момент t* £ (0, Т) выполняется строгое неравенство Nk(t*) < Nk, тогда существует такое t' £ (t*,Т], что [0, t'] С М.
Доказательство. Доказательство того, что [0, t*] С М, легко получить из утверждений 1 и 2, используя, например, метод «от противного». Следуя этим же методом «от противного», докажем, что [t*, t'] С М. Пусть не существует такого t', тогда с учётом этих же утверждений 1 и 2 вытекает, что [t*, t'] С Н и N(t) = N. Значит множество (20) непусто и inf Н = t*. В силу убывания и непрерывности фазовых переменных q(t) получаем N(t*) = N. Пришли к противоречию. □
Утверждение 7. Пусть множества (20) и (21) не пусты, тогда существует такое Ттах £ М, при котором 1(Ттах) = N.
Доказательство. В качестве Ттах берём inf Н, причём Ттах = тахМ. В силу убывания и непрерывности фазовых переменных q(t) получаем И(Ттах) = IN. □
Следствие 2. Пусть множества (20) и (21) не пусты, тогда оптимальная траектория проходит через плоскость, описываемой уравнением
п
X) Чг! = Q, (27)
=1
где вектор q принадлежит положительному ортанту пространства Rп.
Теорема 1. Существует только одно и.з трёх возможных вариантов поведений оптимальной траектории:
п
1) если в начальный момент выполнено неравенство Y1 QiN < Q, то множе-
=1
ство (20) непусто, вектор управлений N(t) = N при всех t £ [0, Т] и сумма
п
Y^ &(t)N{, < Q при t £ (0,Т] и со временем данная сумма убывает;
=1
п
2) если в начальный момент выполнено неравенство Y1 с£1г > Q и существует
=1
допустимая траектория, на которой выполняется равенство
п
Q - ^ v(t)^(t) = 0, (28)
=1
при всех £ [0, Т], то данная траектория является оптимальной, и максимальное .значение функционала (1) равно QТ. Заметим, что к числу допустимых траекторий можно отнести траектории, удовлетворяющие тем же самым условиям текущего пункта теоремы, но при большем .значении Т. Эти траектории также будут оптимальными.
3) если в начальный момент выполнено неравенство ^ д®N° > С? и множе-
г=1
ство (20) непусто, то существует такое Тшах, что равенство (25) выпол-
п _ _
няется при всех £ е [0, Тшах], сумма ^ < при всех £ е (Ттах,Т] и со
г=1
временем данная сумма убывает.
Доказательство теоремы 1 следует из вышеприведённых утверждений.
Теорема 2. Пусть множество (20) непусто и пусть фазовые переменные дг упорядочены в порядке возрастания величин а°г Йг, тогда в каждый момент I: (а) оптимальная траектория вычисляется по формулам
Ш) = ехР
г = 1, 2,...,щ (29)
(b) если выполнено условие ^ < , то N(I) = N;
г=1
(c) в противном случае существует такое целое число к е {1,..., п}, при ко-
к—1__к тором выполняется двойное неравенство ^ < < ^ и
г=1 г=1
Щ€) = Й при г = 1, 2,..., к - 1, (30)
Йг(г) = 0 при % = к + 1,...,п, (31)
Йк(1)
к-1
а - £ ит
г=1
/Ш- (32)
(й) Для любого V е ^,Т] выполняется неравенство N(£) < N
Доказательство. Введём обозначение °(£) = 'фг(Ъ) для г = 1, 2,...,п. В новых обозначениях сопряжённые уравнения (10) перепишутся в следующем виде:
фi(t) = [¥i (*)+ т - 1]а° Йг(€), г = 1, 2,...,п. (33)
Пусть £ е Н. Из утверждения 1 оптимальные управления N(^ = N, а из условия (15) множитель Лагранжа = 0 на множестве (20). Подставим эти значения в формулы (11) и (12). Из решений (19) системы сопряжённых уравнений (13) с учётом формул (11) и (12) получаем
(рф) = фф) = 1 - ехр[-а° ЩТ - ¿)], г = 1, 2,...,п. (34)
Далее с учётом упорядочивания фазовых переменных в порядке возрастания значений N° сравниваем величины функций при фиксированном ¿. В ре-
зультате приходим к следующему выводу: с возрастанием номера компоненты строго возрастают при фиксированном £ значения функций фг(р).
Положим (р\(0) < <р2(0) < ... < <рп(0). При этом £ = 0 е М. В следующей части доказательства теоремы 2 покажем, что такое же соотношение порядка функций фг^р) сохраняется на всем множестве М, т.е. для любого £ е М выполнены неравенства (р\(1) < р2(¿) < ... < <рп(1).
В соответствии с условием £ еМ введём множество С', на котором осуществляется максимизация гамильтониана (9),
С = е К" | 0 , ^qiNi = д| . (35)
Преобразуем гамильтониан (9)
п п
Н(■) = тах У^шт - = Я+ тах ^[-^Ш^Щ. (36)
N ее *—' N ее *—'
i=1 i=1
Максимизация гамильтониана (36) сводится к минимизации функции
п
^ЫШЪЪ]. (37)
N €С
г= 1
Сделаем замену переменных: и = ^; и = (¡N1. В результате максимизация гамильтониана (36) свелась к следующей задаче линейного программирования:
^^ißiUi ^ min (38)
г= 1
при условиях
п
£u = Q, (39)
=1
0 < и < U. (40)
При этом коэффициенты pi упорядочены в порядке строгого возрастания.
С учётом (39) преобразуем линейную функцию (38). В результате получаем следующую задачу линейного программирования на максимум:
п— 1
у^ ( СРп - <Pi)ui ^ max (41)
=1
Заметим, что все коэффициенты в описании линейной функции (41) положительны и убывают с увеличением номера компоненты вектора и.
Следующая процедура позволяет найти такой оптимальный вектор и, который доставляет максимум линейной функции (41) и соответственно минимум линейной функции (38). Сначала на первом шаге выбираем U1. Если U1 ^ Q, то U1 = Q и U = 0, i = 2,... ,п. В этом случае мы определили все оптимальные значения вектора U. Если U1 < Q, то U1 = U1 и делаем замену Q = Q -U1. Далее переходим ко второй компоненте U2 и повторяем процедуру первого шага. Если мы остановились на к-м шаге, то
U = U, г = 1, 2,..., к — 1, Uk G (0,Uk ], U = 0, i = k + 1,...,п.
Если мы прошли все шаги от 1 до п — 1, то U = U, i = 1, 2,... ,п — 1, Un = Q.
Из вышеприведённого доказательства можно утверждать, что для любой к-й компоненты вектора управлений справедливо следующее:
1) если Uk = Uik, то Nk(t) = Nk;
2) если Uk = 0, то Nk(t) = 0;
3) если ик е (0,и), то Йк(г) е (0,14).
На рис. 1 схематично изображены поведения управлений N1^), N2^), N3^), функций ф1(1),ф2(1),ф3(1) и множителя Лагранжа /3(Ь) на отрезке [0,Т] для задачи 1 с размерностью п = 3.
г- т. Г„..... г
Рис. 1. Поведения управлений Nl(t),N2(t),N3(t), функций <pi(t),<p2(t),<p3(t) и множителя Лагранжа ¡3(t) на отрезке [0, Т] для задачи 1 c размерностью
п =3
Рассмотрим в динамике поведение оптимальной траектории на всем множестве М, т.е. от 0 до Tmax. Пусть в начальный момент q0Ni > Q, тогда Ni(0) Е (0,Ni), Ni(0) = 0 для г = 2,3,...,п. В связи c убыванием фазовой переменной (t) существует такой момент п, что:
1) на полуинтервале t Е [0, т\) выполнены соотношения (t) Е (0, N\) и Ni(t) = 0 для г = 2,3,... ,п;
2) Ni (п) = iVi и Ni(Ti) = 0 для г =2, 3,...,п.
Из системы сопряжённых уравнений (33) с учётом (14), (16), (17) при t Е [0, ti) получаем: 71(t) = ¿1 (t) = 0; @(t) = 1 — щ(t) = Pi = const; 'Pi(t) = const для i = 1, 2,... ,n. При этом соотношения порядка, установленного среди компонент вектор-функции (р(t) по их значениям в момент t = 0, не изменяются на отрезке [0,ti ].
Далее существует такой момент T2, что:
1) на интервале (т\,т2) выполнены соотношения N1 (t) = N\, N2(t) Е $,N2) и Ni(t) = 0 для i = 3,... ,п;
2) Щ (T2) = , Й2(Т2) = N2 и Й1(Т2) = 0 для i = 3,...,П.
Из системы сопряжённых уравнений (33) с учётом (14), (16), (17) при t Е (т\,т2) получаем: 72(t) = S2(t) = 0; @(t) = 1 — <p2(t) = @2 = const;'Pi(t) = const для
г = 2,33,...,п; функция p1(t) убывает на интервале (т\, т2). Последнее вытекает из неравенства P1 = 1 -р1(т1) > P2 = 1 -р2(т1). При этом соотношения порядка, установленные среди компонент вектор-функции р(t) по их значениям в момент t = 0, сохраняются на отрезке [0, т2].
В момент Tmax согласно утверждению 7 N(Tmax) = N. Анализируя динамику изменений множителя Лагранжа ft(t) и динамику вектор-функции p(t), приходим к следующим заключениям: множитель Лагранжа ft(t) является кусочно-постоянной убывающей функцией; установленные соотношения порядка в момент t = 0 среди значений компонент вектор-функции p(t) с учётом строгих неравенств остаются без изменений на отрезке [0,Tmax].
Важно, что установленные соотношения строгого порядка среди значений компонент вектор-функции p(t) в момент t = 0 совпадают в момент t = Tmax с соотношениями строгого порядка среди значений компонент (34) этой же вектор-функции , но полученной с учётом условий трансверсальности при t > Tmax. Поэтому можно осуществить синтез вектор-функций p(t) в момент Tmax. Если бы имелись отличия в соотношениях строгого порядка, то синтез вектор-функций р( ) осуществить в момент Tmax было бы нельзя.
Теорема 2 доказана. □
Следствие 3. Пусть множество (20) непусто, тогда оптимальные траектории, полученные при решении задачи 1 и нижеприведённой задачи 2, совпадают на отрезке [0, Tmax] и максимальное значение функционала задачи 2 равно Tmax.
Задача 2. О максимизации длины общей «полки» для группы газовых месторождений
Требуется максимизировать функционал T —> max при дифференциальных связях (2), начальных условиях (3) и ограничениях на управления (4),
п
£ b(t)Nt(t) = (V, (42)
г=1
Управления 1г(г = 1, 2,... ,п) принадлежат множеству измеримых функций.
Доказательство следствия теоремы 2 состоит в следующем. Имеются две задачи 1 с одинаковыми исходными данными. Единственное отличие состоит в значении T. В первой задаче T = Ti, во второй задаче T = T2 и Ti <T2. Для каждой задачи выполнено условие не пустоты множества Н. Сравним оптимальные траектории. Они совпадают на периоде от 0 до T1, также совпадают значения Tmax. При этом Tmax < Ti.
Далее устремим T к Tmax (T > Tmax). В результате задача максимизации накопленной добычи для группы газовых месторождений сводится к очень важной прикладной задаче перспективного планирования для группы газовых месторождений, а именно: к задаче максимума длины общей «полки» для группы месторождений, и эта величина равна Tmax.
5. Обсуждение результатов исследования
Обсудим полученные результаты исследования с точки зрения модели группы газовых месторождений. Все п газовых месторождений отсортированы в порядке
0 _
возрастания величин Ш>1г. Оптимальная политика разработок газовых место-
1
рождений с точки зрения максимума накопленной добычи состоит в следующем. Сначала в момент = 0 вводится в разработку первое месторождение с наи-
о _
меньшим значением величины ^oN . Если мощностей у первого месторождения недостаточно, чтобы осуществлять добычу газа в течениие некоторого времени на
уровне величины общей «полки» месторождений, т.е. д10!^! ^ (, то в разработку вводится второе месторождение. Если мощностей у первого месторождения достаточно, чтобы осуществлять добычу газа на уровне величины общей «полки» месторождений, т.е. д0]У1 > (, то разработка первого месторождения осуществляется с добычей на уровне ((. При этом с самого начала часть скважин отключена. Однако с извлечением запаса газа на первом месторождении средний дебит уменьшается, а значит падает добыча газа. Для компенсации его падения в разработку включаются новые незадействованные скважины. Наступает момент I = Т\, когда в добыче задействованы все имеющиеся у первого месторождения скважины, т.е. д 1(т\)!1 = ((. Дальнейшая компенсация падения добычи осуществляется за счёт ввода в разработку следующего по порядку месторождения. В момент = в разработку вступает второе месторождение, которое за счёт ввода новых скважин компенсирует падение добычи на первом месторождении. Средний дебит на втором месторождении по мере извлечения своего запаса газа также падает. Наступает момент £ = , когда в добыче задействованы все имеющиеся как
2
у первого, так и второго месторождения скважины,т.е. ^ дг(т2)Лг = (. Данная
г=1
процедура продолжается до тех пор, пока в разработку в момент t — Тт&х ни
п _ _
будут включены все скважины на всех месторождениях, т. е.^ (Цг(Ттах)!!г = (.
=1
п
Далее совокупная добыча на всех месторождениях падает и ^ ¿¡г^)Лг < ( при
ге (Тшах,Т]. г=
Объяснение полученных результатов состоит в следующем. Предположим, что имеются два одинаковых месторождения, отличающиеся только количеством действующих скважин. Легко показать, что чем меньше скважин на месторождении, тем больше добытого газа за один и тот же период приходится на одну скважину. Отсюда следует, что сначала необходимо разрабатывать месторождения с меньшим количеством скважин, а затем — с большим количеством скважин.
Автором статьи на примере 2, 3, 4 и 5 газовых месторождений был проведён численный эксперимент. Методом перебора решалась задача на максимум длины «полки» месторождения. Для каждого случая рассматривалось 1000 вариантов. При этом исходные данные для каждого варианта выбирались случайным образом. Полученные в результате численного эксперимента данные подвергались анализу и исследованию. Численный эксперимент подтвердил правильность сделанных в статье выводов.
Литература
1. Эрроу К. Применение теории управления к экономическому росту // Математическая экономика. — М.: Мир, 1974. — С. 7-45. [Ehrrou K. Primenenie teorii upravleniya k ehkonomicheskomu rostu // Matematicheskaya ehkonomika. — M.: Mir, 1974. — S. 7-45. ]
2. Skiba A. K. Optimal Growth with a Convex-Concave Production Function // Econometrica. — 1978. — Vol. 46. — Pp. 527-539.
3. Balder E. J. On Existence Result for Optimal Economic Growth Problems // A Journal of Mathematical Analysis and Applications. — 1983. — Vol. 95. — Pp. 195213.
4. Асеев С. М., Кряжимский А. В. Принцип максимума Понтрягина и задачи оптимального экономического роста // Труды Математического Института им. В.А. Стеклова. — М.: Наука, 2007. — Т. 257. — С. 1-272. [Aseev S. M., Kryazhimskiyj A. V. Princip maksimuma Pontryagina i zadachi optimaljnogo ehkonomicheskogo rosta // Trudih Matematicheskogo Instituta im. V.A. Steklova. — M.: Nauka, 2007. — T. 257. — S. 1-272. ]
5. Скиба А. К. Принцип максимума в задаче максимизации дохода для модели газового месторождения // Вестник РУДН. Серия «Математика. Информатика. Физика». — 2009. — Т. 1. — С. 14-22. [Skiba A. K. Princip maksimuma v zadache maksimizacii dokhoda dlya modeli gazovogo mestorozhdeniya // Vestnik RUDN. Seriya «Matematika. Informatika. Fizika». — 2009. — T. 1. — S. 14-22. ]
6. Федосеев А. В., Хачатуров В. Р. Постановка и исследование задач оптимального управления для анализа перспективных планов в нефтегазодобывающей промышленности // Имитационное моделирование и математические методы анализа перспективных планов развития нефтедобывающей промышленности. — М.: ВЦ АН СССР, 1984. — С. 66-112. [Fedoseev A. V., Khachaturov V. R. Postanovka i issledovanie zadach optimaljnogo upravleniya dlya analiza perspektivnihkh planov v neftegazodobihvayutheyj promihshlennosti // Imitacionnoe modelirovanie i matematicheskie metodih analiza perspektivnihkh planov razvitiya neftedobihvayutheyj promihshlennosti. — M.: VC AN SSSR, 1984. — S. 66-112. ]
7. Маргулов Р. Д., Хачатуров В. Р., Федосеев А. В. Системный анализ в перспективном планировании добычи газа. — М.: Недра, 1991. — 288 с. [Margulov R. D., Khachaturov V. R., Fedoseev A. V. Sistemnihyj analiz v perspektivnom planirovanii dobihchi gaza. — M.: Nedra, 1991. — 288 s. ]
8. Крюков В. А, Скиба А. К., Федосеев А. В. Задачи оптимального управления разработкой газоконденсатного месторождения. — М.: ВЦ АН СССР, 1990. — 40 с. [Kryukov V. A., Skiba A. K., Fedoseev A. V. Zadachi optimaljnogo upravleniya razrabotkoyj gazokondensatnogo mestorozhdeniya. — M.: VC AN SSSR, 1990. — 40 s. ]
9. Моделирование освоения газовых месторождений на заключительной стадии эксплуатации / А. К. Скиба, В. Р. Хачатуров, А. В. Злотов, А. Н. Со-ломатин. — М.: ВЦ РАН, 2006. — 54 с. [Modelirovanie osvoeniya gazovihkh mestorozhdeniyj na zaklyuchiteljnoyj stadii ehkspluatacii / A. K. Skiba, V. R. Khachaturov, A. V. Zlotov, A. N. Solomatin. — M.: VC RAN, 2006. — 54 s. ]
10. Соломатин А. Н. Некоторые оптимизационные задачи для группы газовых месторождений. — М.: вЦ РАН, 2009. — 44 с. [Solomatin A. N. Nekotorihe optimizacionnihe zadachi dlya gruppih gazovihkh mestorozhdeniyj. — M.: VC RAN, 2009. — 44 s. ]
11. Ли Э. Б., Маркус Л. Основы теории оптимального управления. — М.: Наука, 1972. — 576 с. [Li Eh. B., Markus L. Osnovih teorii optimaljnogo upravleniya. — M.: Nauka, 1972. — 576 s. ]
UDC 517.97
An Optimal Control Applied Problem with a Mixed Constraint
A. K. Skiba
Dorodnicyn Computing Centre of the Russian Academy of Sciences 40 Vavilov str., Moscow, Russia, 119333
We consider a gas deposits group functioning model. The optimum control problem with a mixed constraint is stated and solved over a finite horizon. A maximum of accumulated production for a group of gas deposits is taken as the optimization criterion.
Key words and phrases: maximum principle, K. Arrow proposition, mixed constraint, model of gas deposits, maximum accumulated production, maximum length of the "ceiling".