Себестоимость добычи в модели газового месторождения: исследование и применение Текст научной статьи по специальности «Математика»

УДК 519.86

А. К. Скиба, А.Н. Соломатин, В. Р. Хачатуров

Федеральный исследовательский центр «Информатика и управление» РАН

Себестоимость добычи в модели газового месторождения: исследование и применение

Исследуется непрерывная динамическая модель газового месторождения с взаимо-влияюгцими скважинами. Основным критерием разработки газового месторождения является минимальная себестоимость добычи газа. Устанавливается зависимость минимальной себестоимости от горизонта планирования. Рассматриваются проекты разработок двух месторождений, из которых необходимо выбрать одно с наименьшей себестоимостью для последующего освоения, причем пересмотр выбора в дальнейшем невозможен. Решается задача о правильности сделанного выбора для различных значений горизонта планирования. Найдены условия, при которых выбор был сделан правильно и при которых он был ошибочным.

Ключевые слова: непрерывная динамическая модель, модель газового месторождения, себестоимость добычи газа, постоянные и переменные затраты, минимизация себестоимости, проект разработки.

А. К. Skiba, A.N. Solomatin, V.R. Khachaturov Federal Research Center «Computer Science and Control» of the Russian Academy of Sciences

Prime cost of production in the model of the gas field: research and application

The continuous dynamic model of a gas field with mutually influencing wells is investigated. The main criterion for development of the gas field is a minimum prime cost of gas production. The dependence of a minimum prime cost on the planning horizon is established. We consider projects of the development of two fields of which it is necessary to choose one with the smallest prime cost for the subsequent development, and revision of the choice is further impossible. The problem of the correctness of the made choice for various values of the planning horizon is solved. We find conditions such that the choice is made correctly and such that it is erroneous.

Key words: continuous dynamic model, gas field model, prime cost of gas production, constant and variable expenses, minimization of prime cost, a project of development.

1. Введение

В отделе Методов проектирования развивающихся систем ФИЦ ИУ РАН на протяжении многих лет ведутся работы по динамическому планированию разработки нефтяных, газовых и газоконденсатных месторождений как для одного месторождения, так и для группы в целом [1-3]. С учетом современных основ эксплуатации газовых месторождений [4] на основе имитационной модели группы газовых месторождений была разработана и внедрена Система планирования добычи газа, которая широко использовалась при расчетах для месторождений Западной и Восточной Сибири. Система позволяет на заранее заданном временном периоде рассчитывать объемы добычи газа и технико-экономические показатели добычи в динамике по месторождениям. Был разработан ряд оригинальных математических моделей газовых месторождений и решено большое количество задач оптимального

© Скиба А. К., Соломатин А. Н., Хачатуров В. Р., 2018

(с) Федеральное государственное автономное образовательное учреждение высшего образования «Московский физико-технический институт (государственный университет)», 2018

управления [5-6]. Одной из них является проблема оптимального экономического роста [7], решение которой основано на предложениях К. Эрроу [8].

Рассмотренные в [7] экономические показатели добычи газа в основном касались прибыли, дохода и других показателей, связанных с рыночной ценой, для чего имеются вполне обоснованные причины. Однако рыночная цена является слабо прогнозируемой экономической величиной даже на небольшом промежутке времени, что явно было продемонстрировано за последние несколько лет. В формулах, описывающих себестоимость, отсутствует рыночная цена,что дает себестоимости значительные преимущества перед другими вышеупомянутыми показателями. Поэтому можно отнести себестоимость к важнейшим экономическим показателям и уделить ей более серьезное внимание.

В последние годы был проведен ряд исследований по оптимизации себестоимости в смежных областях. Например, оптимальное планирование и развитие инфраструктуры для добычи сланцевого газа обсуждалось в [9], а в [10] рассматривалась общая оптимизационная модель для планирования производства и потребления энергоресурсов в различных отраслях промышленности. Повышенный интерес к проблематике себестоимости подтверждает актуальность темы, рассматриваемой в настоящем работе.

2. Исследование себестоимости газового месторождения

Рассмотрим динамическую модель функционирования газового месторождения с взаи-мовлияющими скважинами [2]. Введём следующие обозначения:

Т— горизонт планирования;

N — действующий фонд добывающих скважин;

N — общий пробуренный и готовый к эксплуатации фонд добывающих скважин, являющийся естественным ограничением на действующий фонд добывающих скважин;

д — средний дебит добывающих скважин;

V — извлекаемый запас природного газа;

с — стоимость природного газа.

Естественно, что все введенные обозначения являются положительными величинами. Между описанными выше переменными устанавливается взаимосвязь, которую запишем в виде системы дифференциальных уравнений

V = -Ид, д = -а Ид, а = (1)

при начальных условиях V0 > 0 д° > 0. На фонд добывающих скважин накладывается следующее ограничение:

0 < N < N. (2)

Предполагается, что N > 0. Капитальные затраты описываются в виде линейной функции г + к]Ч, где г и к - постоянные коэффициенты. В настоящей работе рассматриваются две основные характеристики разработки месторождений: себестоимость добытого газа и прибыль, полученная от его реализации. Себестоимость газа определяется здесь как отношение капитальных затрат в период обустройства месторождения к добытому за плановый период объему природного газа:

8 = ^. (3)

¡дММ

0

Прибыль есть разница между стоимостью всей реализованной продукции и всеми имеющимися издержками:

т

Р = с [ дЫМ - г -Ш. (4)

В издержки включены только капитальные затраты, произведенные в предплановый период. Совокупные текущие затраты малы по сравнению с капитальными затратами, поэтому в дальнейшем их учитывать не будем.

Следует обратить внимание на то, что если вместо цены в (4) подставить себестоимость, то прибыль обратится в нуль. Таким образом, введенное понятие себестоимости оправдывает свое название «себестоимость». Однако если в определении прибыли будем учитывать коэффициент дисконтирования, то получаемая прибыль не будет равна нулю.

Проведём следующую замену переменных:

— — X с

И' = аИ, ТУ' = аТУ, V' = аУ = д, г' = ат, с' = -. (5)

к к

При написании указанных переменных штрихи при них в дальнейшем опускаем. Двойное неравенство (2) будем понимать в новых переменных. Перепишем формулы (1), (3) и (4) в более удобной для дальнейшего анализа форме:

д = -Кд; (6)

5 = к ; (7)

¡дММ о

Р = ^ (с У" дЫМ - г - N | . (8)

Ставятся две задачи оптимального управления: минимизация себестоимости добычи природного газа и максимизации прибыли от его реализации.

Задача 1. Требуется минимизировать функционал (7) при дифференциальной связи (6), начальном, условии д° > 0 и ограничении на управление (2).

Задача 2. Требуется максимизировать функционал (8) при дифференциальной связи (6), начальном, условии д° > 0 и ограничении на управление (2).

Легко показать, что минимум функционала (7) и максимум функционала (8) достигаются при N(¿) = ТУ на всем периоде планирования (£ £ [0,Т]). Преобразовав формулы (7) и (8), получим две функции, зависящие от ТУ и Т:

2 + N

Р(У, Т) = - (сд0(1 - е-Ят) - г - N) .

а \ )

Далее найдем минимум функции (9) при фиксированном Т. Для этого найдем частную производную функцию (9) по Ту и результат приравняем к нулю. После преобразований получаем следующее уравнение:

емт - 1 - РУТ - гТ = 0. (10)

Обозначим левую часть уравнения (10) через Ф(Т, ТУ), где Т и ТУ - положительные переменные. При этом параметр г считаем известной положительной фиксированной величиной. функция Ф(Т, ТУ) определена и непрерывна при Т > 0 и при ТУ > 0, также данная функция имеет непрерывные частные производные дФ/дТ и дФ/дТУ при положительных значениях Т и ТУ, причем дФ/дТХ = Т(емт - 1) > 0. Легко показать, что для любого То > 0 существует единственное ТУо > 0, при котором Ф(Т0,ТУо) = 0. Следовательно, по теореме о неявной функции при Т > 0 существует единственная положительно определенная функция ТУ = N*(Т), удовлетворяющая уравнению (10).

Подставим функцию ТУ*(Т) в уравнение (10). В результате приходим к тождеству

ем*(т)т - ! *(Т)Т - хТ = о. (11)

Воспользовавшись последним тождеством, упростим формулу (9):

5 (М *(Т),Т) = ^ е* *(т)т. (12)

Интересно исследовать поведение функций IV*(Т)Т, IV*(Т) и 5(Ж*(Т),Т) в зависимости от параметра Т € (0, те). Исследование поведения этих функций будет полезно при решении ниже поставленной задачи.

Докажем последовательно три утверждения.

Утверждение 1. Функция N*(Т)Т определена, непрерывна и строго возрастает от, 0 до те при Т > 0.

Разложим экспоненту в тождестве (11) в ряд Маклорена, и после преобразований получаем следующее двойное неравенство:

0 < Ы*(Т)Т < /2гТ. (13)

Устремляя в двойном неравенстве (13) параметр Т к нулю, получаем Ишт^о(Й*(Т)Т) = 0. Продифференцируем обе части тождества (11) по Т, и после преобразований получим

(М*(Т)Т)' = *(т)т _ 1 > (14)

Это означает, что функция IV*(Т)Т строго возрастает. Из тождества (11) вытекает следующий результат: *(Т)Т) = те. Утверждение 1 доказано.

Утверждение 2. Функция N*(Т) определена, непрерывна и строго убывает от те до 0 при Т > 0.

Разделив обе части тождества (11) на Т, продифференцируем его по параметру Т. После преобразований получаем

рN*(т)т _ 1 лг*(Т)ТрМ*(т)т

(^юу=-—ТЧе-^ -) ^—. (15)

Знаменатель дроби функции (15) положителен при Т > 0. Покажем, что числитель этой дроби отрицателен при тех же значениях Т. Числитель дроби равен нулю при Т = 0. Продифференцируем числитель дроби по параметру Т. В результате с учетом утверждения 1 при Т > 0 получаем

- N*(Т)Те^*(т)т(Й*(Т)Т)' < 0.

Следовательно, числитель не обращается в нуль на интервале Т € (0, те) при Т > 0 и он отрицателен. Значит, функция N*(Т) строго убывает. Разделим все части двойного неравенства (13) на Т:

0 < N*(Т)

Перейдя к пределу при Т ^ те, приходим к следующему результату: Ншт^^ N*(Т) = 0.

Разложим функцию е1^*(т)т по формуле Маклорена 2-го порядка с остаточным членом в форме Пеано, и разложение подставим в тождество (11). В результате получаем

*(т )т)2 + ^ *(Т)Т)2 - хТ = 0. ^

(/V* (т)т)2 _

Разделим обе части тождества (16) на (—(2) ) и, учитывая Ншт^о( N*(Т)Т) = 0, перейдем к пределу при Т ^ 0:

1 - Иш

т 0

(( N *(т))2т) 0.

В результате получаем limT^о N*(Т) = те. Утверждение 2 доказано. Утверждение 3. Функция S(NN*(Т),Т), описывающая себестоимость добычи газа, определена, непрерывна, и строго убывает, от те до при Т > 0.

Подставим функцию N*(Т) в функцию (9). Далее продифференцируем функцию S(N*(Т), Т) по параметру Т и, воспользовавшись соотношением (10), получим

(S(N-(Т).Т))' = -kN•(Т)е-"^ЯЪт)2 < 0.

Следовательно, функция S(N*(Т), Т) строго убывает при Т > 0. Воспользовавшись утверждениями 1 и 2, вычисляем пределы:

z + N *(Т)

lim S(N*(Т), Т) = к-—-- ' = те;

т^о v v ' до(i _ e-N*(т)т) '

(Я '*Т ),Т ) = к ^ + ^ ,т) =«

Утверждение 3 доказано.

Теперь перейдем к основной части настоящей работы.

3. Сравнительный анализ себестоимости двух газовых месторождений

Предположим, что имеются два газовых месторождения, из которых необходимо выбрать одно для последующей разработки и добычи газа в течение периода планирования То. Основным критерием выбора для конкурсной комиссии является минимум себестоимости продукции, который определяется на основании предоставленных проектов разработки и обустройства месторождений. Решение конкурсной комиссии является окончательным и пересмотру не подлежит. Выбор конкретного месторождения для осуществления добычи может быть обусловлен и другими причинами. К ним, например, можно отнести возможность подведения газопровода только к одному из двух месторождений.

Очевидно, что в данной ситуации каждый из собственников месторождений желает победить в конкурсе. Поэтому на каждом месторождении с учетом периода планирования То минимизируется себестоимость продукции, что определяет тем самым величину фонда добывающих скважин N*(То). В этой связи ставится естественный вопрос о правильности сделанного выбора: что будет, если вдруг в будущем период планирования изменится с величины То на другое значение Т?

Тогда возникает следующая задача. На основании имеющихся данных проекты разработок пересчитываются. Если выбор был сделан правильно, то к конкурсной комиссии нет вопросов. Если выбор сделан неправильно, то возникает конфликтная ситуация вплоть до судебных разбирательств.

Таким образом, на основе нового значения периода планирования Т создается новый проект обустройства, в основе которого лежит задача минимизации себестоимости, и определяется другая величина фонда добывающих скважин N * (Т). Себестоимость в этом случае описывается функцией 5(ТУ*(Т),Т).

Задача 3. Даны, пять положительных чисел ]Х*(То), 3\(ТУ**(То),То), Щ(То), 32(Щ(То),То) и То. Первые два, числа, от,носятся, к первому месторождению, а последующие два числа, - ко второму месторождению. Последнее число является периодом планирования добычи газа. Предполагается выполнение неравенства

Бг(Щ(То), То) < Б2(К*(То),То). (17)

Необходимо найти т,акие значения Т, при которых неравенство

31(ТУ**(Т),Т) < Б2(ТУ*(Т),Т)

выполняется, и такие значения Т, при которых оно не выполняется.

Согласно (12), значения себестоимости рассчитываются по формулам

Зг(Мг*(Т),Т) = ^(т)т, 1 = 1,2-

Переменные с индексом г = 1 относятся к первому месторождению, а с индексом г = 2 - ко второму месторождению. На основании исходных данных определяем конкретные значения двух чисел:

Ъг = Зг(Мг*(То), То)Тое(т°)т°, 1 = 1,2. Тогда себестоимости вычисляются по формулам

8г(Й*(Т),Т) = (т)т, г = 1,2. (18)

Разделив себестоимости (18) друг на друга, введем обозначение, которое в дальнейшем будем называть отношением себестоимостей К*(Т):

Я*(Т) = Я2(Щ(Т ),Т) = рЩ(т )т-IV* (т )т

( ) Б^^Т),Т) Ьг .

Используя новое обозначение, из неравенства (17) получаем ограничение на значение функции К*(Т) в точке То:

1 < К*(То) = ЫЩТд П) = ъ_2ещ(т°)т°-м?(т°)т°. (19)

( о) Зг(Щ(То),То) Ьг к ;

ем*(т)т - 1 - х*(Т)Т - г.т = 0, ¿ = 1,2. (20)

Конкретные числовые значения параметров и 22 можно определить из исходных данных по формулам

еЩ(т°)т° _ 1_ ^*(то)То

гг = --1 +\(1о)+о, 1 = 1,2. (21)

^о

Справедливо следующее утверждение.

Утверждение 4. Для каждого Т > 0 большему значению г соответствует большее значение N*(Т).

Из тождества (20) при г = 2 вычтем тождество (20) при 1 = 1. После элементарных

2>

е щ(т)т - Щ(Т)Т - ещ(т)т + Щ(Т)Т = (22 - 2г)Т > 0.

Учитывая то, что функция ех - х возрастает в зависимости от своего положительного аргумента, получаем Й*(Т)Т > ]Ч*(Т)Т. Отсюда ]Ч%(Т) > ]Ч*(Т). Если же хг > х2, то IV* (Т) > Щ(Т). Утверждение 4 доказано.

Следующее утверждение вытекает из предыдущего.

Утверждение 5. Если ЩТо) < Я*(То), то Щ(Т) < Щ(Т) при Т > 0, и если Щ(То) > Й*(То),то Щ(Т) > Щ(Т) при Т> 0

Утверждение 6. Функция ¡(х) = ^Хг убывает на интервале (0, те) от, 1 до 0.

( х)

Г(х) = еХ((1ех-_ХI- 1. (22)

Знаменатель дроби (22) при х > 0 положителен. Покажем, что числитель при х > 0 отрицателен. При х = 0 числитель дроби (22) равен нулю. Продифференцируем числитель дроби. В результате получаем -хех. Следовательно, на интервале (0, те) числитель

х > 0 ( х)

( х)

зано.

Утверждение 7. Функция К*(Т) монотонно изменяется на, интервале (0, те) от, ^ до . Она строго возрастает, при ]Х*(То) < ]Х*(То) и строго убывает при ]Х*(То) > Щ(То). Далее будем исследовать функцию

¡(Т) = е**(т)т-**(т)т, (23)

которая отличается от В*(Т) на постоянную величину. Продифференцируем функцию

¡(ту.

/'(Т) = еЩ(т)т-Щ(т)т[Т(Щ(Т)Т)' - Т(¡У*(Т)Т)']/Т. (24)

Воспользовавшись равенством (14) и тождествами (20), преобразуем выражение в квадратных скобках соотношения (24):

Т(ТУ* (Т)Т)' - Т(ТУ**(Т)Т)' = Х2Т ХгТ

ем*(т)т _ 1 ещ(т)т _ 1 ]Х**(Т )Т Щ(Т )Т

(25)

ем*(т)т _ 1 еЩ(т)т _ 1

Из утверждений 5 и 6 с учетом соотношений (25) и (24) получаем, что при Р1*(Т0) < Р12*(Т0) функция (23) строго возрастает, а при Р1*(Т0) > Р12*(Т0) - строго убывает.

Легко показать справедливость следующих пределов:

lim ¡(Т) = 1; lim ¡(Т) = —. (26)

Следует отметить, что при поиске значения второго предела в (26) использовались пределы

е "*(т )т

lim —-— = Zi i = 1,2, т^ж Т

которые можно получить из тождеств (20) и утверждений 1 и 2. Утверждение 7 доказано. Переформулируем задачу 3 в терминах функции К*(Т).

Т

функция R*(r) больше 1 и при которых она меньше 1.

Легко показать, что из исходного предположения (19), числовых значений параметров (21), утверждений 4, 5 и 7 вытекает теорема 1, которая является решением задачи 3. Теорема 1.

Случай 1. Пусть 1У**(Т0) < 1У*(Т0). Б этом случае ^ > 1 и > 1, тогда:

а) функция R*(r) строго возрастает, на, интервале (0, те) от ^ до

б) если ^ < 1, то существует, единственное положительное число Т*, меньшее Т0, такое, что R*(r) = 1, причем при Т > Т* выполняется неравенство R*(r) > 1, а, при 0 <Т <Т* оно не выполняется;

в) если > 1 , т,о неравенство R*(Т) > 1 выполняется, при всех положительных

Т

Случай, 2. Пусть ТУ**(Т0) > ТУ*(Т0). В этом случае ^ < 1 и ^ > 1, тогда: а) функция, R*(Т) строго убывает на, интервале (0, те) от ^ до

б) если ^ < 1 ; то существует, единственное положительное число Т*, большее Т0, такое, что выполняется неравенство В*(Т) = I, причем при 0 < Т < Т* выполняется неравенство В*(Т) > 1, а при Т > Т* оно не выполняется;

в) если > 1, т,о первенство К*(Т) > 1 выполняется при всех положительных

Т

Случаи 3. Пусть К*(Т0) = М*(Т0). В этом случае = 1 и ^ > 1, тогда:

а) функция R* (Т) на интервале (0, те) принимает, постоянное значение, равное

б) выполняется неравенство R*(Т) > 1 при всех положительных значениях Т.

4. Заключение

В настоящей работе мы рассмотрели непрерывную динамическую модель газового месторождения. В рамках модели были выписаны формулы для расчета прибыли и себестоимости добычи газа. Прибыль считается самым важным экономическим показателем, только прибыль может полностью отражать уровень эффективности деятельности любого предприятия. Однако, в отличие от себестоимости добычи, цена на газ (которая включена в расчеты прибыли) представляет собой слабо прогнозируемое динамическое значение. Этот факт значительно усложняет прогноз совокупной прибыли. Себестоимость определяется как отношение капитальных затрат при строительстве месторождения к объему природного газа, добываемого в плановый период. Капитальные затраты описываются как линейная функция одной переменной с постоянной и переменной частями, причем переменная часть определяется фондом добывающих скважин.

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

Полученные результаты позволяют принимать правильные решения в следующей ситуации. В конкурсе на разработку участвуют два месторождения; для разработки выбирается одно их них с наименьшей себестоимостью добычи газа, причем пересмотр выбора невозможен. В дальнейшем период планирования изменяется, происходит перерасчет себестоимо-стей и других показателей для каждого месторождения. Возникает вопрос: правильным ли был сделанный выбор? В случае отрицательного ответа возникает конфликтная ситуация со всеми вытекающими последствиями.

Литература

1. Маргулов Р.Д., Хачатуров В.Р., Федосеев A.B. Системный анализ в перспективном планировании добычи газа. М.: Недра, 1992.

2. Хачатуров В.Р., Солом am, ин А.Н., Злот,ов A.B. \и др.]. Планирование и проектирование освоения нефтегазодобывающих регионов и месторождений: Математические модели, методы, применение / под ред. В.Р. Хачатурова. М.: УРСС: ЛЕНАНД, 2015.

3. Хачатуров В.Р., Солом,am,ин А.Н., Скиба А.К. Моделирование разработки газовых месторождений // Труды МФТИ. 2017. Т. 9, № 3. С. 162-171.

4. Вяхирев Р.И., Коротаев Ю.П., Кабанов Н.И. Теория и опыт добычи газа. М.: Недра, 1998.

5. Крюков В.А., Скиба А.К., Федосеев A.B. Задачи оптимального управления разработкой газоконденсатного месторождения. М.: ВЦ АН СССР, 1990.

6. Скиба, А.К. Принцип максимума в задаче максимизации дохода для модели газового месторождения // Вестник РУДН. Серия Математика. Информатика. Физика. 2009. № 1. С. 14-22.

7. Skiba А.К. Optimal growth with a convex-concave production function // Econometrica. 1978. V. 46, N 3. P. 527-539.

8. Эрроу К. Применение теории управления к экономическому росту / Матем. экономика. М.: Мир, 1974. С. 7-45.

9. Arredondo-Ramirez К., Ponce-Ortega J.M., El-Halwagi М.М. A Optimal planning and infrastructure development for shale gas production // Energy Conversion and Management. 2016. 119. P. 91-100.

10. Elsholkami M., Elkamel A. General optimization model for the energy planning of industries including renewable energy: A case study on oil sands // AIChE Journal. 2017. 2(63). P. 610-638.

References

1. Margulov R.D., Khachaturov V.R., Fedoseyev A.V. The system analysis in advance planning of gas production. M.: Nedra, 1992. (in Russian).

2. Khachaturov V.R., Solomatin A.N., Zlotov A.V., et al., Planning and design of development of oil and gas extraction regions and fields: Mathematical models, methods, application. Ed. by V.R. Khachaturov. M.: URSS: LENAND, 2015. (in Russian).

3. Khachaturov V.R, Solomatin A.N., Skiba A.K. Modeling of gas fields development. Proceedings of MIPT. 2017. V. 9, N 3. P. 162-171. (in Russian).

4. Vyakhirev R.I., Korotaev Y.P., Kabanov N.I. Theory and experience of gas production. M.: Nedra, 1998. (in Russian).

5. Kryukov V. A., Skiba A.K., Fedoseyev A.V. Problems of optimal control of the development of the gas-condensate field. M.: CC AS USSR, 1990. (in Russian).

6. Skiba A.K. Maximum principle in a problem of income maximization for the model of gas deposit. Bulletin of RUDN. Series Mathematics. Informatics. Physics. 2009. N 1. P. 14-22. (in Russian).

7. Skiba A.K. Optimal growth with a convex-concave production function. Econometrica. 1978. V. 46, N 3. P. 527-539.

8. Arrow K. The application of control theory to economic growth. Mathematical economics. M.: Mir, 1974. (in Russian).

9. Arredondo-Ramirez K., Ponce-Ortega J.M., El-Halwagi M.M. Optimal planning and infrastructure development for shale gas production. Energy Conversion and Management. 2016. 119. P. 91-100.

10. Elsholkami M., Elkamel A. General optimization model for the energy planning of industries including renewable energy: A case study on oil sands. AIChE Journal. 2017. 2(63). P. 610-638.

Поступим в редакцию 10.09.2018