К проблеме оптимального регулирования процесса извлечения нефти с учетом вероятностного закона распределения некоторого параметра слоистого пласта

У

правление техническими системами и технологическими процессами

УДК 05.13.18, 01.02.05

К ПРОБЛЕМЕ ОПТИМАЛЬНОГО РЕГУЛИРООАНИЯ ПРОЦЕССА ИЗВЛЕЧЕНИЯ НЕФТИ С УЧЕТОМ ВЕРОЯТНОСТНОГО ЗАКОНА РАСПРЕДЕЛЕНИЯ НЕКОТОРОГО ПАРАМЕТРА СЛОИСТОГО ПЛАСТА1

Р.А. Султанов, В.Д. Слабнов, И.А. Фукин, В.В. Скворцов

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

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

ВВЕДЕНИЕ

Необходимость повышения объемов нефти, извлекаемой из пластов сложного геологического строения, характеристики которых изменяются по толщине и простиранию пласта, делает актуальной проблему оптимального регулирования процесса её извлечения. Для её решения используются усредненные математические модели. Преимущество двухмерных моделей по сравнению с трехмерными связано с недостатком геологической информации о строении пласта и значительным сокращением времени численных расчетов благодаря уменьшению числа ячеек конечно-разностной сетки [1].

В данной работе используется модель Баклея — Леверетта процесса вытеснения нефти водой из слоистых пластов, не учитывающей гравитационные и капиллярные силы [2]. Принимаются во внимание четыре гипотезы о вытеснении нефти

1 Основное содержание статьи доложено на IV Международной конференции по проблемам управления, Москва, Институт проблем управления им. В.А. Трапезникова РАН, 26—30 января 2009 г.

водой из элементарного объема слоистого пласта. Предполагается, что абсолютная проницаемость слоистого пласта может изменяться по толщине и простиранию. Вероятностный закон ее распределения в слоистом пласте определяется на основе заданных параметров распределения: математического ожидания, дисперсии и некоторого дополнительного условия. Полученный вероятностный закон используется для построения «фиктивных» (модифицированных) относительных фазовых проницаемостей в слоистом пласте [2]. Усреднение трехмерной модели по толщине пласта позволяет уменьшить размерность исходной модели на единицу.

Численное решение задачи управления для построенной усредненной математической модели позволяет определить симплекс-методом оптимальные дебиты и забойные давления источников и стоков при задании для них двусторонних ограничений и некоторых вариантов граничных условий в элементе заводнения. Критерием оптимизации служит максимизация целевой функции — суммарного отбора нефти [3, 4, 6, 7].

Цель данной работы — построение математической модели двухфазной фильтрации в слоистом

пласте с учетом четырех гипотез о вытеснении нефти водой из элементарного объема; постановка задачи оптимального управления; разработка численного алгоритма ее решения; анализ результатов решения задачи оптимального регулирования извлечения нефти из слоистого пласта при варьировании параметров распределения для дискретных аналогов равномерного и логарифмически-нормального законов распределения абсолютной проницаемости; сравнение результатов численного решения задачи без учета и с учетом регулирования.

1. МАТЕМАТИЧЕСКАЯ МОДЕЛЬ ПРОЦЕССА 20

ДВУХФАЗНОЙ ФИЛЬТРАЦИИ В СЛОИСТЫХ ПЛАСТАХ, НЕОДНОРОДНЫХ ПО ТОЛЩИНЕ И ПРОСТИРАНИЮ

Рассмотрим нефтяной пласт, в котором происходит процесс вытеснения нефти водой. Кровлю и подошву пласта будем считать непроницаемыми. Введем прямоугольную декартову систему координат х, у, z. Выделим элементарный объем в виде столбца, ограниченного сверху кровлей, а снизу — подошвой. Его основанием будет служить элементарная площадка йх, йу, а высотой — толщина пласта Н.

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

др др

и -£ , где р — пластовое давление, постоянны

дх ду

по вертикали и, в частности, не зависят от у [2].

Распределение абсолютной проницаемости в слоистом пласте с учетом дискретного аналога вероятностного закона распределения случайной величины можно задавать с помощью гистограммы или ряда распределения, приведенного в табл. 1, где в . — вероятность того, что абсолютная прони-

цаемость примет значение к.

а(і')'

Таблица 1

Распределение абсолютной проницаемости

к*

в Рі Р2 Рз Рпк

Сделаем еще ряд предположений относительно механизма вытеснения нефти водой. Пусть каждая у-я среда с абсолютной проницаемостью ка. имеет свою остаточную нефтенасыщенность 5*н(/ и долю связанной воды 5*/ которые считаются известными. В этом случае максимальные значения объемов воды и нефти в каждом элементарном столбце имеют следующий вид

5в(/) 1 %(/')’

0') 1 5*в(/)’

где буквенный индекс «в» относится к водной фазе, «н» — к нефтяной фазе.

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

Пусть и — водо- и нефтенасыщенность у-й среды соответственно.

Определим активные нефте- и водонасыщен-ность у-й среды как

1 - І =

Яв

Отсюда нефте- и водонасыщенность у-й среды можно представить в виде

5н(/) = 8*н(У)(1 - / + 5**(/),

5в(/) = 5*в(/)9./ + 5в№ (1 - е/').

Тогда нефте- и водонасыщенность в элементарном объеме имеют следующий вид:

Ян = I %(/)(! - 0і)Рі + I 4(і) ЄД,

і= і

"к

У'1

"к

і= і

ГГ

Яв = I Ч

і= і і= і

Очевидно, что

"к "к

1 іі + 1 (1 - (ІД- = 1-і = і і= і

Обозначим через к^^) и к^^) относительные фазовые проницаемости воды и нефти в у-й среде соответственно. Как видно из формулы (1), они будут зависеть от величины 9...

Составим уравнение материального баланса для /-й жидкости:

д Нка кФ др + 3 'Нка кФ др

д х і 1 д у 1 1= У

"

"

к

к

где / = в, н, ка — осредненная абсолютная проницаемость, т — пористость пласта; цв и цн — вязкости воды и нефти.

В данной математической модели «фиктивные»

относительные фазовые проницаемости кф Ц), / = в, н, построены на основе известных максимальных значений относительных фазовых проницаемостей, активных насыщенностей сред и вероятностного закона распределения абсолютной проницаемости по толщине пласта.

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

Систему уравнений (2) после несложных преобразований и учета закона Дарси можно привести к системе уравнений относительно усредненной

скорости фильтрации V и нефтенасыщенности 5: Шу

где V = — г (я^гаёр — суммарная усредненная скорость фильтрации нефти и воды, 9 — усреднен-

f k£

ный дебит скважины, г (5) = ка

V М-н М-в

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

пласте, і (s) = ka

k^ ( s -

Цн

/ т (s) — усредненная доля

нефти в потоке жидкости;

k^ (sJ = f I РЛАлЮ,

k ^ а(/) н(/)' н'’

kа /= І

k^ (sb) = k1 I PA/WO.

kа / = І

Для решения системы уравнений (3) требуется знать вид «фиктивных» относительных фазовых проницаемостей (4) по некоторому вероятностному закону распределения случайной величины или вид ряда распределения абсолютной проницаемости в каждой точке усредненного слоистого пласта.

2. ПОСТАНОВКА ЗАДАЧИ ОПТИМАЛЬНОГО РЕГУЛИРОВАНИЯ ПРОЦЕССА ИЗВЛЕЧЕНИЯ НЕФТИ

Постановка задачи регулирования при двухфазной фильтрации с учетом ввода скважин в эксплуатацию на основе математической модели Бак-лея — Леверетта (3) с начальным и граничными условиями и алгоритм ее решения была дана в работе [5]. Поставим в соответствие размещению к-й скважины в области фильтрации Б ее номер. Номера добывающих скважин отнесем к множеству / *, нагнетательных — к множеству /\/*, где / — множество номеров всех скважин. Предполагается, что срок Т разработки залежи определяется временем отключения последней действующей добывающей скважины. Отключение к-й добывающей скважины происходит при Тк < ? < Т, где Тк — время достижения к-й скважиной заданного предельно-допустимого значения текущей обводненности ее продукции 1 — Д Тк) = ^в*к. Будем рассматривать регулирование процесса нефтедобычи из неоднородной залежи в условиях принятой модели по критерию максимизации текущей нефтеотдачи на заданные моменты времени = 0, 1, 2,...

Математическая задача формулируется следующим образом: найти максимальный суммарный отбор нефти

+ 1

£ 1 4к(^нк(?г)й?г ^ та^ К е К + 1^

к е /*( гг) *и

для заданного планового периода [?л, 1п + 1] и ограничений на режимы действующего фонда скважин

р* т рзк т р*, 9* т |9к| т q*, у к е /*ф,

с учетом момента времени tr начала задачи прогнозирования, который определяется из условия

*г = ^кф) 1 t(^в*к (tr)), Ук е /*(tr),

где номера добывающих и нагнетательных скважин соответствуют оптимальному подмножеству / *ф и /(д\/ *(tr).

Здесь 9к(^) — дебит к-й скважины на tr-ю дату регулирования, ^вк(^) — доля воды в потоке жидкости для к-й скважины на ^-ю дату регулирования, рзк(^) — забойное давление к-й скважины на ^-ю дату регулирования.

При численном решении задачи управления нефтеизвлечением на каждую контрольную дату tr восстанавливается оптимальное поле давления в области фильтрации Б в результате численного

n

n

решения задачи линеиного программирования (5)—(7) [3, 5, 7, 8]:

£ aCTk(s(tr))( p ^ (tr) - psk(tr))FHk(tr)dtr ^ max,

k є J*

tr Є [tn, tn + І^

при ограничениях

Т - upi - У + °/ + Vpi + V + GJ- - ^ - І + aJ' + ^ + І +

J - ^ J -

+ (ctJ' + + ^зк =0

p* m pзk m p*, q* m jaa (s(tr))(pj (tr) - Jpзk(tr»j m q*,

vk є J*

где i = 1, Mx, j = 1, My, Mx и M — число узлов по оси OX и 07 соответственно, рзк — забойное давление скважин с координатами (хк, ук), р^. — сеточное давление по области фильтрации Б, 8 ^ — сеточный аналог функции Дирака, рф и р* — нижнее и верхнее ограничения на забойные давления сква-

жин соответственно, q* и q *

нижнее и верхнее

ограничения на дебиты скважин соответственно.

Найденное поле давления вместе со значениями оптимальных забойных давлений и дебитов действующих скважин обеспечивают при решении задачи управления максимальный суммарный отбор нефти из неоднородного пласта на ^-ю контрольную дату. Здесь линейное программирование используется одновременно в качестве метода решения системы алгебраических уравнений, полученных из разностного аналога (6) для первого уравнения системы (3).

3. ВЫЧИСЛИТЕЛЬНЫЕ ЭКСПЕРИМЕНТЫ

Рассмотрена модель неоднородного пласта прямоугольной конфигурации Б = {0 < x < а, 0 < у < Ь}, вскрытого девятиточечным элементом заводнения (рис. 1), для a = Ь = 850 м, И = 5 м, m = 0,2, мн = 3 мПа-с, мв = 1 мПа-с, 5* = 0,26, 5 * = 0,8, рк = 22 МПа, р* = 12,5 МПа, р* = 22 МПа,

3 3

9Ф = 1 м /сут, 9* = 500 м /сут, йх = йу = 25 м. Здесь йх, — шаги разностной сетки по осям координат.

х у

Относительные фазовые проницаемости кн(5) и кв(5) — линейные полиномы [5]. Проведена серия вычислительных экспериментов для элемента заводнения в соответствии со следующими гипотезами о характере вытеснения нефти водой.

1. Вода вытесняет нефть из участков (зон) пласта в первую очередь с наилучшей абсолютной проницаемостью, т. е. каждый элементарный объем разбивается на две зоны: водную и нефтяную

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

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

3. Вода вытесняет нефть пропорционально найденному переменному коэффициенту.

4. Обобщение гипотез 1 и 3. Каждый элементарный объем разбивается пороговым значением на две зоны: водную и нефтяную. Причем в зоне с наихудшей абсолютной проницаемостью нефтенасыщенность максимальна, а в зоне с наилучшей абсолютной проницаемостью — происходит неполное вытеснение нефти водой пропорционально найденному переменному коэффициенту.

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

На рис. 1 показано размещение 9 скважин в элементе заводнения для следующего варианта задания граничных условий: на одной стороне задано граничное условие I рода, контур питания КП, на трех остальных сторонах — граничное условие II рода, непроницаемая граница НГ Математическое ожидание абсолютной проницаемости ка(х, у) в области двухфазной фильтрации Б элемента за-

НГ

НГ

Рис. 1. Элемент заводнения слоистого пласта

a

Рис. 2. Графики 1—4 коэффициента нефтеотдачи ц(т) для гипотез 1—4 соответственно. Дискретный аналог логарифмичес-

ки-нормального закона, ка = 0,1 мкм2; а2 = 0,001

воднения принимает следующие значения: &а = 0,1;

2 2 0,3 мкм , дисперсия — а = 0,001; 0,01.

Вычислительные эксперименты проводились для дискретных аналогов двух законов распределения:

— ряд распределения с одинаковыми вероятностями того, что случайная величина ка по толщине пласта примет данное значение (дискретный аналог равномерного закона распределения случайной величины, табл. 2);

— ряд распределения с различными вероятностями того, что случайная величина ка по толщине пласта примет данное значение (дискретный аналог логарифмически-нормального закона распределения случайной величины, табл. 3).

В табл. 4 приведены значения коэффициента прокачки в поровых объемах и соответствующие им значения коэффициента конечной нефтеотдачи для двух законов распределения случайной величины с учетом заданий математического ожидания, дисперсии и гипотез 1, 2 и 4 о вытеснении нефти водой из элементарного объема слоистого пласта. Анализируя данные табл. 4 и графики, представленные на рис. 2, можно сделать вывод, что гипотезы 1 и 2 соответствуют нижнему и верхнему графикам значений коэффициента нефтеотдачи п от коэффициента прокачки поровых объемов т.

Таблица 2

~а = 0,1 мкм2 а2 = 0,001

ка(і) 0,051 0,063 0,075 0,087 0,100 0,112 0,124 0,136 0,149

в 0,111 0,111 0,111 0,111 0,111 0,111 0,111 0,111 0,111

Таблица 3

Логарифмически-нормальный закон

~а = 0,1 мкм а2 = 0,001

0,058 0,078 0,098 0,118 0,137 0,157 0,177 0,196 0,216

Ру 0,050 0,400 0,300 0,100 0,050 0,040 0,030 0,020 0,010

Таблица 4

Значения коэффициентов прокачки и нефтеотдачи

Вид закона распределения Дискретный аналог равномерного закона Дискретный аналог логарифмически-нормального закона

Номер гипотезы 1 2 4 1 2 4

м к м 01 ,1 ,0 0 0, = = ^ 2 Т 3,744 3,624 3,697 3,737 3,635 3,675

П 0,791 0,847 0,831 0,807 0,853 0,837

2 м к м 01 ,3 0 0, 0, = = ^ 2 Т 3,719 3,702 3,719 3,725 3,72 3,719

П 0,819 0,836 0,829 0,823 0,839 0,829

2 м к м 1 ,3 01 0, 0, = = ^ 2 Т 3,747 3,617 3,719 3,738 3,625 3,677

П 0,788 0,848 0,829 0,806 0,854 0,836

Рис. 3. Графики 1, 4 и 1_орг, 4_орг коэффициента нефтеотдачи п(т) для гипотез 1 и 4 без учета и с учетом оптимального регулирования соответственно. Дискретный аналог равномерного за-

22

кона, ка = 0,1 мкм ; а = 0,001

Характер расположения графиков сохраняется

как для одинаковых ка, но разных а , так и для

2 ~ одинаковых а , но разных ка. Графики коэффициента нефтеотдачи п от коэффициента прокачки поровых объемов т для гипотез 3 и 4 расположены между графиками для гипотез 1 и 2. Показатели извлечения нефти из слоистого пласта для гипотезы 4 лучше, чем для гипотезы 3 в случае дискретного аналога равномерного и логарифмически-нормального законов.

На рис. 3 приведены кривые зависимости коэффициента нефтеотдачи от коэффициента прокачки поровых объемов без учета и с учетом оптимального регулирования для дискретного аналога логарифмически-нормального закона в случае принятия гипотез 1 и 4. Анализ графиков показывает значительное превышение коэффициента текущей нефтеотдачи с учетом оптимального регулирования над аналогичным показателем без его учета.

Таблица 5

Дискретный аналог равномерного закона с учетом оптимального регулирования

Номер гипотезы 1 2 3 4

~а = 0,1 мкм2 а2 = 0,001 Т 2,907 2,48 2,659 2,405

Ап 0,01 0,015 0,011 0,084

АВНФ 2,046 2,593 2,449 3,367

Рис. 4. Графики 1_орг — 4_орг коэффициента нефтеотдачи п(т) с учетом оптимального регулирования для гипотез 1—4 соответ-

2

ственно. Дискретный аналог равномерного закона, ка = 0,1 мкм ; а2 = 0,001

В табл. 5 для гипотез 1—4 и дискретного аналога равномерного закона приведены значения т коэффициента прокачки поровых объемов и соответствующие им значения абсолютных разностей для коэффициентов конечной нефтеотдачи Ап и водо-нефтяного фактора АВНФ = 2А(0в/2А10н (2А(0в, 2А(0н — суммарная добыча воды и нефти) с учетом оптимального регулирования и без его учета, соответственно.

На рис. 4 приведены графики зависимости коэффициента нефтеотдачи п от коэффициента прокачки поровых объемов т с учетом дискретного аналога равномерного закона и гипотез 1—4 для оптимального регулирования. Анализируя кривые, можно сделать вывод о том, что наилучшие показатели получены для гипотезы 4, наихудшие — для гипотезы 1. Для гипотез 2 и 3 данные показатели расположены между двумя графиками, соответствующими гипотезам 1 и 4. Это объясняется видом относительных фазовых проницаемостей для гипотез 1—4 и динамикой изменения поля скоростей нефти в области двухфазной фильтрации элемента заводнения слоистого пласта.

ЗАКЛЮЧЕНИЕ

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

Разработан вычислительный алгоритм для определения модифицированных относительных фазовых проницаемостей на основе этих гипотез для произвольного дискретного аналога вероятностного закона распределения абсолютной проницаемости.

Предложена постановка двумерной нестационарной задачи оптимального регулирования извлечения нефти из слоистого пласта. Приведен алгоритм ее решения. Критерий качества решения — максимальный суммарный отбор нефти. Использование оптимального регулирования улучшило показатели извлечения нефти и водонефтяного фактора по сравнению с показателями прогнозирования.

Описанный подход позволяет решать задачи оптимального управления в элементе заводнения с нерегулярной сеткой скважин и заданием граничных условий 1-го или 11-го рода и их комбинации.

ЛИТЕРАТУРА

1. Северов Я.А. Сравнение гидродинамических показателей двухмерного с трехмерным моделированием пласта // Науч. тр. Оренб. нефтян. акцион. компании «ОНАКО». — 2001. — № 3. — С. 165—169.

2. Булыгин Д.В., Булыгин В.Я. Геология и имитация разработки залежей нефти. — М.: Недра, 1996. — 382 с.

3. Мееров М.В. Исследование и оптимизация многосвязных систем управления. — М.: Наука, 1986. — 236 с.

4. Закиров Э.С. Трехмерные многофазные задачи прогнозирования, анализа и регулирования разработки месторождений нефти и газа. — М.: Грааль, 2000. — 303 с.

5. Слабнов В.Д., Волков Ю.А., Скворцов В.В. Влияние некоторых факторов регулирования на основные показатели не-фтеизвлечения из неоднородного пласта // Математическое моделирование. — 2002. — Т. 14, № 1. — С. 3—15.

6. Dolle N., Brouwer D.R., Jansen J.D. Dynamic Optimization of Water Flooding with Multiple Injectors and Producers using Optimal Control Theory // Proc. XIV Intern. Conf. on Computational Methods in Water Resources. — Delft, 2002. — P. 23—28.

7. Slabnov V.D., Volkov Yu.A. To a problem of optimum oil production process control from seams, development with application of horizontal wells // Вторая междунар. конф. по проблемам управления (Москва, 17 — 19 июня 2003 г. / ИПУ РАН) // Избранные труды в двух томах. — М., 2003. — Т. 2. — С. 57—67.

8. Азиз Х, Сеттари Э. Математическое моделирование пластовых систем. — М.: Недра, 1982. — 407 с.

Статья представлена к публикации членом редколлегии Л.П. Боровских.

Султанов Рустем Анварович — канд. физ.-мат. наук, доцент, Казанский государственный финансово-экономический институт, (843) 291-13-71, И sultanov@mail.ru,

Слабнов Виктор Дмитриевич — канд. физ.-мат. наук, доцент, Институт экономики, управления и права, г. Казань,

(843) 231-92-90, И slabnov@ieml.ru,

Фукин Игорь Анатольевич — канд. физ.-мат. наук, доцент, зав. кафедрой, Институт экономики, управления и права, г. Казань, (843) 231-92-90, И figor@ieml.ru,

Скворцов Владимир Викторович — д-р техн. наук, профессор, Казанский государственный технологический университет,

И vskvortsov@mail.ru.

“Новая

книга

Новиков Д.А. Теория управления образовательными системами. — М.: Народное образование, 2009. — 452 с.

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

Управление образовательными системами рассмотрено с общекибернетических позиций, что позволяет изложить компоненты и структуру теории систематически. Значительное внимание уделено программно-целевому (проектному) подходу к управлению образовательными системами. Детально рассмотрены функции, задачи и методы управления образовательными системами. Их изложение сопровождено анализом практических ситуаций, с которыми сталкиваются руководители образовательных систем и на которых иллюстрируется применение теории.

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

По вопросам приобретения обращайтесь в ООО «НИИ школьных технологий»

(109202, Москва, шоссе Фрезер, д. 17, 3-й этаж, офис 1.

Я/1 (499) 270-20-65; Н market@narodnoe.org).

34

CONTROL SCIENCES № б • 2009