CC BY
157
УДК 519.71:622.3
ИНТЕГРИРОВАННЫЕ МОДЕЛИ И АЛГОРИТМЫ ИДЕНТИФИКАЦИИ ДОБЫЧИ НЕФТИ С УЧЕТОМ АПРИОРНОЙ ИНФОРМАЦИИ
Д.В. Севостьянов
Томский государственный университет систем управления и радиоэлектроники E-mail: [email protected]
Рассматривается задача прогноза добычи нефти и оценки извлекаемых запасов с учетом дополнительной априорной информации и экспертных оценок технологических параметров разработки месторождений. Приводятся результаты анализа точности оценок в зависимости от числа лет разработки.
Введение
В последние годы в нефтяной промышленности в связи с широким внедрением информационных систем, современных технологий исследования скважин, систем регистрации разнообразной информации о состоянии объектов разработки, скважин, нефтяных пластов и месторождений возрос интерес к задачам построения математических моделей процессов нефтегазодобычи [1-6].
Проблемы построения математических моделей процессов нефтегазодобычи связаны с решением задач идентификации, которые часто называют обратными задачами. Задача идентификации заключается в построении оптимальных, в смысле заданных критериев качества, математических моделей технологических параметров разработки (добычи нефти, жидкости, воды, пластовых давлений, обводненности продукции и т. д.) с использованием промысловых данных, результатов комплекс -ных исследований скважин и нефтяных пластов. Задачи идентификации технологических параметров подразделяют на два крупных направления, имеющие свои цели и задачи.
Первое направление представляет задачи уровня проектирования разработки нефтяных месторождений, которые решаются крупными коллективами в научных центрах нефтяных компаний, проектных институтах [6]. На стадии проектирования разработки месторождения создаются цифровые геолого-технологические модели нефтяных месторождений, позволяющие прогнозировать показатели разработки на достаточно длительный период (20-30 лет), создавать технологические схемы и проекты разработки, определять стратегию развития компании.
Второе направление представляют задачи идентификации уровня мониторинга и оперативного управления разработкой нефтяных месторождений. На стадии разработки месторождений актуальным является моделирование технологических параметров для оперативного решения задач прогнозирования добычи нефти, оценки эффективности геолого-технических мероприятий, определения оптимальных режимов работы скважин и т. п.
Здесь более мобильными, легко адаптируемыми (настраиваемыми) на основе промысловых данных и результатов исследований скважин по сравнению
с геолого-технологическими моделями уровня проектирования являются различные регрессионные статические и динамические модели технологических параметров разработки. Наиболее широкое применение получили промыслово-технологи-ческие модели, основанные на характеристиках вытеснения, уравнениях фильтрации флюидов, малопараметрических регрессионных моделях добычи нефти, жидкости и воды [7, 8].
Однако реальные условия нефтегазодобычи характеризуются неполнотой, неоднородностью, недостаточностью исходных данных о фильтрационных параметрах и энергетическом состоянии нефтяных пластов, различными ошибками регистрации технологических и геологических параметров разработки, неадекватностью моделей и т. п. [1-5]. В данных условиях использование классических методов идентификации технологических параметров разработки часто связано с проблемами устойчивости и низкой точностью получаемых решений [5].
Перспективным направлением преодоления реально существующей сложности процессов нефтегазодобычи, недостаточности данных, повышения точности и устойчивости решений, является использование интегрированных систем идентификации с учетом разного рода дополнительной априорной информации, накопленного опыта и знаний [9, 10].
Интегрированные системы идентификации предоставляют возможность объединять математические (физически содержательные) модели технологических параметров разработки месторождений и формализованные экспертные оценки лица, принимающего решения, в единую интегрированную систему моделей, отражающую целостные свойства реальных систем, что дает интегральный (си-нергетический) эффект при решении разнообразных задач идентификации и управления [11].
Данная работа является развитием технологии интегрированных моделей при идентификации одного из важных технологических параметров разработки месторождений - добычи нефти с использованием малопараметрических регрессионных моделей отбора нефти и непараметрического метода аналогов [10] в условиях априорной неопределенности о статистических характеристиках дополнительных априорных данных и экспертных оценок.
Интегрированные модели и алгоритмы идентификации добычи нефти
Рассмотрим интегрированную систему моделей добычи нефти на основе малопараметрической модели, отражающей эволюцию текущего отбора нефти во времени, с учетом дополнительной априорной информации об извлекаемых запасах, параметрах модели добычи нефти и прогнозных значениях добычи, которую представим в удобном для исследований матричном виде:
е*=е* а)+4,
S = s(T ,а) +П = f f (T,a)dr+n,
0
a, = a +vt, l = l...d, QH = Q(T,a) +s,
(1)
где 0„*=(0*(/г), г=1...и), й(а)=(Д/;,а), /=1...и) - векторы значений добычи нефти, полученные в процессе разработки месторождения и на основе модели /(/, а) за соответствующие промежутки времени Д1=1—1-_х, /=1...и (год, меся—); £=(-, ]=1...р), а=(а]Ь ]=1...ш), /=1...й, -=(&(+), ]=1...г) -векторы экспертных оценок извлекаемых запасов нефти за время разработки Т, дополнительных априорных данных о параметрах модели добычи а=(а1,а2,...,ат) и прогнозных значений добычи нефти, полученные из разных источников информации, например, из проектных технологических документов на разработку месторождений; 4,п,У1,/=1...с1,е - векторы случайных величин с нулевыми математическими ожиданиями, представляющие погрешности измерения добычи нефти, ошибки, вызванные неадекватностью модели добычи нефти, ошибки дополнительных априорных сведений и экспертных оценок.
Рассмотрим решение задачи параметрической идентификации модели (1) с использованием комбинированного критерия качества
Ф(а, р) = ||е* - дя (а)||2 +1- s(T,а )|£ +
+Х W2i(a, -a)2 + 1 \Qh - f (a)||
(2)
непараметрического метода аналогов [10] ,соста-вленного из частных квадратичных критериев качества, где ||Х||2 - норма X,
W = diag( К (
s(T,a0) -sj
ß
L), j = 1... p),
W2l = diag( К (■
a0ß -aj,
ß2
j = 1...m),
W3 = diag( К (■
f (K +Tj,a°) -Qi (tn +rt)
ß3
), j = 1..r) -
диагональные матрицы весовых функций K(u) с управляющими параметрами ß,ß2,ß3, учитываю-
щих значимость (вес) дополнительных априорных данных и экспертных оценок, а0=(а10,а2(',...ат0) -вектор начальных приближений параметров модели добычи нефти.
Функции (ядра) К(и) введены по аналогии с непараметрическими оценками плотности вероятности и регрессии, обладают следующими свойствами [12]:
К(х - и)/к) ^ 1, к К(х - и)/к) ^ 0, к ^ 0.
Задача параметрической идентификации (1) сводится к решению оптимизационной задачи по определению вектора параметров модели добычи нефти
а'(Р) = а^шт Ф(а, в) (3)
а
и оптимизационной задачи по определению вектора управляющих параметров
||2
ß' = argmm(|QH - QH(a\ß)|| ),
(4)
где arg min fix) обозначает точку минимума x* функции fix).x
Приведем решение оптимизационной задачи (3) с использованием метода Ньютона [13]:
a' = a'-1 + h'Aa'-1, i = 1,2,3,..., (5)
где оценки приращений вектора параметров Aa'-1 на каждой итерации определяются путем решения системы линейных уравнений
A'-1Aa'-1 = B'
(6)
где
A'-1 = (DTDH + DT WD + JWy + DWD),
j=1
B'-1 = (DTHeH + 'esWpT +fw2j aO, + ^-WpT )'-
Dh =
(
DT =
df (t' ,a) da..
i = 1...n, j = 1...m
df (tn + zt ,a)
da t
i = 1...r, j = 1...m
матрицы
D. =
ds(T ,a) da,
частных
\
производных;
j = \...тп
- вектор частных произ-
водных от функции извлекаемых запасов л(Г,а); е„=(е„*-е„), ё^Я-з^а0)), Да/=(а/-а?), ет=(б„-б(г,а0)) - векторы невязок.
Для получения системы линейных уравнений (6) необходимо:
1. Разложить функции¡(1,а), й(Т,а) в ряд Тейлора и ограничиться членами первого порядка малости относительно приращения вектора параметров Да":
i=1
i=1
/ м=/ (а+£ (Щ^па
]=1
да,
,'(Г , а) = ,'(Г , а0) + £ (^^"Аа0.
.=1
За..
(7)
2. Используя разложение (7) перейти от модели (1) к линейной интегрированной системе моделей добычи нефти относительно приращения вектора параметров А а0:
№ = В*Аа0 +§,
е, = В0Аа0 +ц,
—О 0 Аа I =Аа +v,, I = ,
е = ВТ Аа0 + е.
(8)
3. Для модели (8) по аналоги с (2) составить комбинированный функционал качества
2
Ф(Аа ) = е - Бк А а +
е, - Б, Аа
+£ W2I (Аа0 - Аа0 )2 + Це^О - В?Аа(
(9)
дФ(Аа )
5Аа0
о
= -(П0)г^(е0 -БОАа0) -
(11)
где 0*,8, -;,/=1..Д0„ - векторы имитируемых значений добычи нефти, дополнительных априорных дан-
ных и экспертных оценок извлекаемых запасов, параметров модели и прогнозных значений добычи, полученные с использованием функции регрессии [7, 10],
Щ (/, а) = а1 ехр(-а 2 /) • ^
(12)
широко используемой при моделировании промысловых данных в качестве представления зависимости годовой добычи нефти от времени разработки; £,п,^,/=1...С,е - векторы случайных величин, полученные с использованием датчика псевдослучайных чисел, распределенных по нормальному закону с нулевыми средними значениями и единичной дисперсией; с, у'=1...4 константы, представляющие относительный уровень ошибок.
Для определения точности оценок прогноза добычи нефти
вЖ +т) = Щ +Т, а\р'))
и извлекаемых запасов
Б * =| Щ (Г, а'
(13)
(14)
4. Вычислить частные производные от функционала (9) по параметрам Аа" и приравнять их к нулю
использованы их относительные ошибки
5. ©Ж + т)) = аЬ8(Щ +т -Щ С. + т (15)
(Б *) = аЪ8(
Щ (гп +т, а)
,(Г,а*(Р') - ,(Г,а)) ,'(Г, а )
)
(16)
-(D0)гW1(e, -Б0Аа))-^21 (А а , -Аа0)-
I =1
-(D:)гW3(e0т- в: А а0) = О. (10)
После проведения несложных преобразований из (10) следует система линейных уравнений (6).
Следует отметить, что для решения оптимизационной задачи (4) по определению вектора управляющих параметров целесообразно использовать методы оптимизации без вычисления производных [13].
Точность оценок прогноза добычи нефти
и извлекаемых запасов
Рассмотрим вопросы точности алгоритмов прогнозирования добычи нефти и оценки извлекаемых запасов на примере интегрированной стохастической системы моделей (1), которую представим в удобном для статистического моделирования матричном виде:
е* = а (а )(1 + Ф,
Б = ',(г, а )(1 + с2п), а I = а (1 + с3У 1, I = 1..Д = в(т, а )(1 + с4)е,
в зависимости от длительности разработки ¡„=п (годы) при значениях параметров добычи нефти а;=50000, а2=2,1, а3=0,26. Здесь а*(в) - оценки параметров модели добычи нефти (12), полученные по алгоритму (6). Для определения оптимальных значений управляющих параметров в' (4) использован метод деформированного многогранника [13].
В табл. 1 приведены результаты расчета относительных ошибок прогноза добычи нефти (15) (т=3) и извлекаемых запасов (16) в зависимости от вида интегрированной модели добычи нефти (11) и уровня ошибок дополнительных априорных сведений. Для указания вида модели добычи нефти введены следующие обозначения: {0*} - модель добычи нефти (апр-орная информация не учитывается); {0„*,£}, {0„*,£, -}, {0„*,£, -,Ш - интегрированные системы моделей с учетом экспертных оценок извлекаемых запасов, дополнительных данных о параметрах моделей и прогнозных значениях добычи нефти, для которых введены условные обозначения М1-М4.
Для оценок параметров а* модели добычи нефти М1 и начальных приближений параметров а0 в алгоритме (5) использован метод Гаусса-Ньютона
[14]
а' = а' + И, А а1
, = 1,2,3,...,
(17)
где на каждом шаге I решается система линейных уравнений вида
(ВтнВн)'-1 •Аа'-1 = (Вгнен)'-1, которая следует из (6) при в=0, /=1...3.
W
I=1
Уровни ошибок дополнительных априорных данных и экспертных оценок имитировались следующие: добыча нефти 5 % (с1=0,05), извлекаемые запасы 2, 5 и 10 % (с2=0,02; 0,05; 0,1), дополнительные априорные данные параметров модели добычи - 2 % (с3=0,02), а ошибка прогноза добычи нефти равна 5 %, что соответствует приемлемым на практике уровням точности оценок. Длительность разработки составляла три года (п=3). Количество дополнительных данных и экспертных оценок в (11) для каждого типа данных равнялось трем (р=й=г=3). Вес дополнительных априорных данных и экспертных оценок при формировании матриц Ж2;, Ж3 (6) рассчитывался с использованием функции К(х)=ехр(-х2).
В табл. 2, 3 приведены значения относительных ошибок прогноза добычи нефти (13) (т=3) и извлекаемых запасов (14) в зависимости от числа лет разработки и вида интегрированной модели добычи нефти. Уровень ошибок имитируемых значений дополнительных данных об извлекаемых запасах и прогнозных значениях добычи составлял 5 %. Уровень ошибок имитируемых значений дополнительных априорных сведений о параметрах модели добычи нефти составлял 2 %.
Таблица 1. Относительные ошибки оценок прогноза добычи нефти и извлекаемых запасов, %
Оценки
Вид интегрированной модели добычи нефти
2
М2
3
М3
4
М4
Уровень ошибок экспертных оценок извлекаемых запасов
2 5 10 2 5 10 2 5 10
Прогноз добычи 53 2,6 3,2 5,1 2,2 2,8 4,5 1,9 2,3 3,1
Извлекаемые запасы 49 2,3 3,1 4.6 1,8 2,6 3,5 1,4 2,2 2,8
Модели добычи нефти
М1
М2
М3
М4
Длительность разработки (номер года)
1 2 3 4 5 6 7 8 9 10
- - 53 42 23 16,1 11,3 8,7 6,9 5,3
51 39 23 15,6 12,4 9,7 7,1 5,6 4,8 4,2
45 34 19 13,6 9,7 7,4 5,2 4,7 4,1 3,6
26 14 8,4 6,9 5,2 3,9 3,1 2,6 2,3 2,1
Таблица 3. Относительная ошибка оценок извлекаемых запасов, %
Модели добычи нефти
М1
М2
М3
М4
Длительность разработки (номер года)
1 2 3 4 5 6 7 8 9 10
- - 49 37 21 14,6 9,3 7,2 5,7 4,9
46 34 20 15,6 11,6 9,5 7,2 5,9 4,6 3,8
40 20 17 12,3 10,2 8,1 5,9 5,2 4,6 3,3
25 13 8,8 6,5 5,1 3,7 3,2 2,8 2,3 1,9
прогноза добычи нефти и извлекаемых запасов на ранних стадиях разработки месторождения (3-6 лет) по сравнению с оценками (17) на основе модели добычи нефти М1, где априорная информация не учитывается. Например, относительная ошибка извлекаемых запасов, полученная к третьему году разработки с использованием метода Гаусса-Ньютона (17), составляет 49 % (табл. 3), а с использованием метода интегрированных моделей с учетом априорной информации (5) порядка 9 %.
Приведем результаты прогноза добычи нефти с использованием фактических промысловых данных годовой добычи нефти одного из месторождений Тюменской области за 15 лет его разработки с 1990 г., изображенных на рисунке (линия 1). Линия 2 представляет прогнозные значения добычи нефти (13) по данным трех лет разработки п=3 и т=1,2,...,12, полученные с использованием модели (12) и алгоритма (17). Линия 3 отражает прогноз добычи нефти, с использованием дополнительной априорной информации об извлекаемых запасах 6,3 млн т, рассчитанных до начала разработки месторождения.
Таблица 2. Относительная ошибка оценок прогноза добычи нефти, %
1 2 3 4 5 6 7 8 9 10 11 12 13 14 15
Номер года разработки
Рисунок. Фактические значения добычи нефти (линия 1), прогнозные значения добычи с использованием модели М1 и М2 (линии 2, 3)
В табл. 4 отражены значения относительных ошибок прогноза добычи нефти, начиная с четвертого года разработки месторождения, полученные с использованием модели М1, алгоритма (17) и модели добычи нефти с учетом априорной информации об извлекаемых запасах М2 и алгоритма (5).
Таблица 4. Относительная ошибка оценок прогноза добычи нефти, %
Модели добычи нефти
М1
М2
Длительность разработки (номер года)
3 4 5 6 7 8 9 10 11 12 13 14
14 11 9 8,5 7,4 4 3 2,9 2,5 1,8 1,7 1,6
6 4,4 3,8 4,1 3,6 2,5 1,8 2,1 1,9 1,1 0,9 0,8
Из табл. 1-3 видно, что учет дополнительных априорных данных и экспертных оценок приводит к существенному увеличению точности оценок
Из рисунка и табл. 4 видно, что использование дополнительных априорных данных об извлекаемых запасах позволяет существенно повысить точность прогноза добычи нефти на ранней (первой) стадии разработки месторождения.
Выводы
Предложенные интегрированные модели добычи нефти и алгоритмы идентификации позволяют:
1) учитывать дополнительные априорные данные и экспертные оценки технологических параметров разработки месторождений: извлекаемые запасы, прогнозные значения добычи нефти, параметры модели добычи нефти и т. п.;
2) получать оценки прогноза добычи нефти и извлекаемых запасов в условиях априорной нео-
пределенности о статистических характеристиках ошибок дополнительных априорных сведений и экспертных оценок;
3) значительно, в два раза и более, повысить точность оценок прогноза добычи нефти и извлекаемых запасов при малом объеме промысловых данных на ранней (первой) стадии разработки месторождения в течение первых пяти лет по сравнению с оценками метода Гаусса-Ньютона, где априорная информация не учитывается.
СПИСОК ЛИТЕРАТУРЫ
1. Мирзаджанзаде А.Х., Хасанов М.М., Бахтизин Р.Н. Этюды о моделировании сложных систем нефтегазодобычи. Нелинейность, неравновесность, неоднородность. - Уфа: Гилем, 1999. - 462 с.
2. Булыгин В.Я., Булыгин Д.В. Имитация разработки залежей нефти. - М.: Недра, 1990. - 240 с.
3. Костюченко С.В., Ямпольский В.З. Мониторинг и моделирование нефтяных месторождений. - Томск: Изд-во НТЛ, 2000. - 246 с.
4. Хургин Я.И. Проблемы неопределенности в задачах нефти и газа. - М.-Ижевск: Институт компьютерных технологий, 2004. - 320 с.
5. Севостьянов Д.В., Сергеев В.Л. Интегрированные системы идентификации показателей разработки нефтяных месторождений // Доклады ТУСУР. - 2004. - № 2 (10). - С. 87-93.
6. Регламент составления проектно технологических документов на разработку нефтяных и газонефтяных месторождений. РД 153-39-007-95. - М.: ВНИИОЭНГ, 1996. - 202 с.
7. Шахвердиев А.Х. Унифицированная методика расчета эффективности геолого-технических мероприятий // Нефтяное хозяйство. - 2001. - № 5. - С. 44-48.
8. Пьянков В.Н. Алгоритмы идентификации параметров модели Баклея-Леверетта в задачах прогноза добычи нефти // Нефтяное хозяйство. - 1997. - № 10. - С. 62-65.
9. Сергеев В.Л. Интегрированные системы идентификации. -Томск: Изд-во НТЛ, 2004. - 240 с.
10. Сергеев В.Л. Идентификация систем с учетом априорной информации. - Томск: Изд-во НТЛ, 1999. - 146 с.
11. Кориков А.М., Севостьянов Д.В., Сергеев В.Л., Сергеев П.В. Адаптивные интегрированные системы идентификации и управления: Вопросы проектирования и развития // Электронные средства и системы управления: Труды Междунар. на-учно-практ. конф. - Томск: Изд-во ИОА СО РАН, 2005. - Т. 2. - С. 58-61.
12. Тарасенко Ф.П. Непараметрическая статистика. - Томск: Изд-во Том. гос. ун-та, 1975. - 292 с.
13. Пантелеев А.В., Летова Т.А. Методы оптимизации в примерах и задачах. - М.: Высшая школа, 2002. - 544 с.
14. Демиденко Е.З. Линейная и нелинейная регрессии. - М.: Финансы и статистика, 1981. - 300 с.
УДК 622.692
АППРОКСИМАЦИЯ СИГНАЛОВ ПРИ ПОСТРОЕНИИ АВТОМАТИЧЕСКИХ СИСТЕМ ОБНАРУЖЕНИЯ УТЕЧЕК НА НЕФТЕПРОВОДАХ
В.А. Якупов, А.М. Кориков*
Томский политехнический университет *Томский государственный университет систем управления и радиоэлектроники E-mail: [email protected], [email protected]
Рассмотрены проблемы построения систем автоматического обнаружения утечек, предложены варианты решения задачи доставки обрабатываемых данных в центры их сбора и обработки. Показана целесообразность использования аппроксимации данных для сжатия и фильтрации текущих измерений.
Введение
Развитие транспортной сети нефте- и нефте-продуктопроводов (НПС) сопровождается непрерывным повышением уровня их технической оснащенности, широким внедрением современных средств компьютеризации и автоматизации процессов перекачки. Трубопроводная транспортировка неразрывно связана с выполнением комплекса мероприятий по охране окружающей среды на принципиально новых научно-технических ос-
новах проектирования, строительства и эксплуатации магистральных трубопроводов. Несмотря на очевидные преимущества трубопроводного транспорта при транспортировании нефти и нефтепродуктов по магистральным трубопроводам, имеет место загрязнение атмосферы, водоемов и почвы. В связи с этим серьезное внимание уделяется сокращению потерь на магистральных трубопроводах. Как правило, потери связаны с нарушением правил эксплуатации, повреждением трубопрово-
