УДК 519.688:552.578.2.061.4
МЕТОД АДАПТИВНОЙ ИДЕНТИФИКАЦИИ ГИДРОДИНАМИЧЕСКИХ ИССЛЕДОВАНИЙ СКВАЖИН С УЧЕТОМ АПРИОРНОЙ ИНФОРМАЦИИ
В.Л. Сергеев, В.С. Аниканов
Томский политехнический университет E-mail: SergeevVL@ignd.tpu.ru
Рассматриваются метод адаптивной идентификации параметров нефтяных пластов в процессе проведения гидродинамических исследований скважин с учетом априорной информации. Приводятся пример определения пластового давления и фильтрационных параметров однородно пористых нефтяных пластов по данным исследований скважины нефтяного месторождения на неустановившихся режимах фильтрации по кривой восстановления давления.
Ключевые слова:
Идентификация, гидродинамические исследования скважин, априорная информация, нефтяные пласты. Key words:
Identification, hydrodynamic analysis of oil wells, a priori information, oil pools.
Гидродинамические исследования скважин (ГДИС) обеспечивают получение важнейших параметров нефтяных пластов, на основании которых осуществляются процессы добычи нефти, составляются технологические проекты разработки месторождений, создаются геолого-технологические модели процессов нефтегазодобычи. Задача идентификации для ГДИС состоит в построении оптимальной, в смысле заданных показателей качества, модели гидродинамических параметров скважины (забойное давление, дебит, динамический уровень, температура и т. п.) и оценке неизвестных фильтрационных параметров, энергетического состояния и геометрических параметров нефтяных пластов.
Следует отметить, что многие классические методы решения задач идентификации для ГДИС (метод касательных, наилучшего совмещения, детерминированных моментов и т. д.) не гарантируют устойчивость оценок параметров пласта и допустимую точность решений [1, 2]. Неустойчивость и низкая точность решений часто проявляются в условиях малых выборок, когда по ряду технических причин, в том числе и в целях сокращения времени простоя скважины, уменьшения потери добычи нефти, требуется прервать исследования. Изменение режимов работы скважин, влияние границ пласта, неоднородностей приводит к появлению дефектных (аномальных) значений забойных давлений и дебитов исследуемых скважин, что является причиной неустойчивости решения. Другая особенность классического подхода связана с планированием исследований и проведением интерпретации результатов ГДИС после их завершения, что часто приводит к простоям скважин и значительным затратам.
Наиболее общим продуктивным решением задачи идентификации ГДИС с использованием системного подхода является предложенный в [3, 4] метод интегрированных моделей, позволяющий объединить модели гидродинамических параметров исследуемой скважины, модели дополнительных априорных сведений и экспертные оценки в единую систему моделей, что обеспечивает устойчивость и повышает точность оценок параме-
тров нефтяных пластов, позволяет получать оптимальные, согласованные решения.
Данная работа является развитием технологии интегрированных моделей при идентификации и интерпретации ГДИС, где для сокращения времени простоя скважин предлагаются модели и алгоритмы, позволяющие определять параметры нефтяных пластов в процессе проведения исследований, не планируя заранее время их завершения.
В основе предлагаемого метода адаптивной идентификации использована интегрированная система моделей гидродинамических параметров скважин с учетом априорной информации вида:
| У = f0(an, О + 4, n = 1,2,3,...
[z = F (an,y,x) +Ц, (1)
где y', x„* измеренные на скважине в моменты времени tn значения гидродинамических параметров у, x=(xbx2,...,xr); z=(z k, k=1,p), - вектор значений дополнительных априорных данных параметров нефтяных пластов и скважин;f0, F=(f1,/2,...,fp) - известные с точностью до вектора параметров a„=(o!,k„ k=l,m) функции; у=1,/>) - слу-
чайные величины, представляющие погрешности измерений гидродинамических параметров скважин, ошибки априорных данных, экспертных оценок, неточности моделейf0, Fa и т. п.
Метод адаптивной идентификации ГДИС с учетом дополнительной априорной информации заключается в последовательном, по мере получения информации у„*, x„*, вычислении вектора оценок параметров an(1) путем решения оптимизационной задачи an (ß) = arg min Ф( J0 (an), ßkJk (an), k = 1, p,),
aneRm
n = 1, nk, (2)
анализа точности полученных приближений и принятии решения о прекращении исследований. Здесь запись arg min fx) означает точку минимума x* функции f (x)(f (x*)=minf(x)); Ф - векторный показатель качества (1), со стоящий из частного показателя качества J0(an) модели гидродинами-
Управление, вычислительная техника и информатика
ческих параметров и взвешенных весами ßk частных показателей качества Jk(an) моделей дополнительных априорных сведений; nk - объем измерений забойного давления, необходимых для обеспечения требуемой точности оценок an'(ß).
Решение о прекращении исследований может быть принято на основе визуального анализа графика зависимости приближений (2) от времени (см. рис. 1-3) либо по критерию стабилизации оценок, где, например, за момент их завершения nk, принимается то значение n, при котором выполняется неравенство
|(аи*_,- -ап) /ап | < eps, i = 1,2,3,...
(eps - заданная точность). (3)
Следует отметить, что при значениях управляющих параметров, равных нулю, априорная информация не учитывается, и оценки (2) при использовании квадратичного показателя качества /0(а„) совпадают с классическими оценками метода наилучшего совмещения [1, 2]
а'п(0) = arg min(УД^) = £(я - /а, V))2). (4)
п m г=1
Рассмотрим решение задачи (2) на примере идентификации и интерпретации нестационарных ГДИС на неустановившихся режимах фильтрации по кривой восстановления давления (КВД) однородного пористого нефтяного пласта, где в качестве модели забойного давления скважины часто используется функция [1, 2]:
p (t)=/о (t, а, а2) = p (to)+а +а2 in(t). (5) В (5)
а, =
Чо
4ла
-In
v у
а2 =
Чо
4ла
где а и % - гидропроводность и пьезопроводность нефтяного пласта; д0 - дебит скважины в момент ее остановки /0; Р(/0) - начальное значение забойного давления; гс=0,1 м - радиус скважины. Для (5) интегрированная система моделей КВД по аналогии с (1) примет вид
|Уп = «1п + «2п 1п(П) + 4 > П = 1 Пк,
|«Ь = «1п = а-п (6)
где уп*=Р„*-Р(/0); п=1,пк - измеренные в моменты времени Iп значения забойного давления скважи-
ны; ai =
Чо
4ла
1n(% а 2 =
Чо - -
—= о; а, % - допол-4л о
нительные априорные данные о пьезопроводности и гидропроводности нефтяного пласта;
При выборе комбинированного показателя качества Ф(ап) для системы моделей (6) в виде свертки частных показателей
Ф(ап) = Jо а) + ZßtJt (ап) = = \\У' - Ра„
k=1
+ \\а-ап
оптимизационная задача (2) сводится к решению системы линейных алгебраических уравнений
(ATW0 A + W )а'п (ß) = (AT W0Y *+W\ ),
(8)
(7)
где запись ||Х||Ж2 означает квадратичную форму ХТШХТ; А=(1, 1п(п), п=1,пк) - матрица размерности (пкх2); У*=(у„*, п=1,пк) - вектор; сё=(а1,а2) - априорные данные о параметрах модели КВД (6); Ж=ё1а§ (w(n), п=1,пк) - диагональная матрица весовой функции определяющая веса забойных давлений Р(Т„) в текущий момент времени /„;
(Р,Р2) - диагональная матрица управляющих параметров Р=(Р,Р2), определяющих значимость (вес) априорных данных а. Для получения (8) достаточно взять частные производные по параметрам ап от комбинированного функционала (7) и приравнять их к нулю.
Следует отметить, что выбор весовой функции V (0 в частном функционале качества /0(а„) (7) зависит от принятой стратегии обработки КВД. Например, при весовой функции вида м(п)=1, Уп=1,пк все значения забойного давления «равноправны», что целесообразно в случае адекватной модели КВД (5). При м(п)=1, пе[пьп2], м(п)=0, пг[пьп2], п2>п1 в обработке участвуют данные, попавшие в выбранный для интерпретации участок КВД [п„+п1,п„+п2]. Для метода адаптивной идентификации целесообразно использовать стратегию «скользящего интервала» [п„+п1+/, п„+п2+/], 1</<пк-пн-п2. Здесь I - номер интервала, пн - количество измерений забойного давления в начальном участке КВД [/0,4], который, как правило, исключают в связи с проявлением там спектра трудно формализуемых эффектов влияния ствола, скин-фактор скважины и т. п. [1, 2].
На рис. 1-3 приведены оценки пластового давления, гидропроводности и пьезопроводности, полученные при интепретации результатов гидродинамических исследований скважины № 323 нефтяного месторождения Тюменской области на неустановившихся режимах фильтрации по КВД. Обработка результатов ГДИС проводилась с использованием метода адаптивной идентификации (АИ) (2) для модели КВД (5) и метода наилучшего совмещения (НС) (4), который, по аналогии с (8), сводится к решению системы линейных уравнений:
(ЛТЖ0 Л)ап (0) = ЛТЖ0У
Оценки пластового давления Рш(п), гидропро-водности ап* и пьезопроводности %п* рассчитаны с использованием стратегии «скользящего интервала» при п2=10, п1=0 и пн=6
Р* (п)=Р (О+«Пп (Р*)+а1 (Р) 1п(Т),
ап = Чо14™'^ (Р'\
%% = ГС ехР(аП1 (р)/ а\п(р)),
где а1п(Р*), а'2п(Р*) - оценки параметров модели КВД (5), полученные к моменту времени /„; Т=200 ч; Рп* - оценка управляющего параметра Р матрицы Wa=pI (где I - единичная матрица) из решения оптимизационной задачи
р'п = ащшт \У* - ¥ап (рп
методом золотого сечения. В качестве априорных сведений использовались данные о гидропроводно-сти <г=1 Дсм/сП и пьезопроводности _-=4.10-4 м2/с.
Таблица. Результаты обработки КВД скважины № 323 нефтяного месторождения Тюменской области
Методы Время, ч Пластовое давление, ат Пьезопро-водность, 10-4, м2/с Гидропро-водность, Дсм/сП
«Раг6у$1ет» 99 152,3 4,53 1,3
Адаптивной идентификации 40 153,4 5,00 1,4
60 153,7 4,97 1,4
80 153,6 4,97 1,4
99 153,4 4,99 1,4
Наилучшего совмещения 40 161,1 3,11 0,9
60 154,4 4,88 1,4
80 148,4 13,20 3,7
99 152,3 6,18 1,7
Рис. 1. КВД (кривая 1). Оценки пластового давления методами АИ и НС (кривые 2, 3)
В таблице приведены оценки пластового давления, пьезопроводности и гидропроводности, полученные с использованием зарубежного программного комплекса «Рап8у81еш» при обработке всей КВД за 99 ч исследований и методами адаптивной идентификации и наилучшего совмещения за разные периоды исследований.
Из рис. 1-3 и таблицы видно, что предложенный метод адаптивной идентификации с учетом априорной информации при интерпретации ГДИС дает устойчивые и более точные оценки параметров нефтяного пласта по сравнению с оценками, полученными методом наилучшего совмещения. Из таблицы видно, что для получения оценок параметров методом адаптивной идентификации
с заданной по критерию (3) точностью (ер8=0,01), достаточно 50 ч, что позволяет в два раза сократить время исследований по сравнению с технологией обработки данных, реализованной в программе «Рап8у81еш».
Время, ч
Рис. 2. Оценки гидропроводности. Метод: 1) АИ; 2) НС
Время, ч
Рис. 3. Оценки пьезопроводности. Метод: 1) АИ; 2) НС
Выводы
Предложен метод адаптивной идентификации нефтяных пластов в процессе гидродинамических исследований скважин с учетом априорной информации.
Метод заключается в последовательном, по мере получения исходных данных, вычислении вектора оценок параметров моделей нефтяных пластов, с анализом точности полученных приближений и последующем принятии решения о прекращении исследований.
Показано, что использование метода адаптивной идентификации при интерпретации нестационарных гидродинамических исследований скважин на неустановившихся режимах фильтрации по кривой восстановления давления дает возможность повысить точность оценок фильтрационных параметров и энергетического состояния однородно пористых нефтяных пластов, сокращает прости скважин и потери добычи нефти.
4
5 2
СПИСОК ЛИТЕРАТУРЫ
1. Кульпин Л.Г., Мясников Ю.А. Гидродинамические методы исследований нефтегазовых пластов. - М.: Недра, 1974. -200 с.
2. Шагиев Р.Г. Исследование скважин по КВД. - М.: Наука, 1998, 304 с.
3. Сергеев П.В., Сергеев В.Л. Идентификация гидродинамических параметров скважин на неустановившихся режимах
фильтрации с учетом априорной информации // Известия Томского политехнического университета. - 2006. - Т. 309. -№5. - С. 156-161. 4. Сергеев В.Л. Интегрированные системы идентификации. -Томск: Изд-во НТЛ, 2004. - 240 с.
Поступила 04.05.2010г.