УДК 519.688:622.276.5.001.42
АДАПТИВНАЯ ИНТЕРПРЕТАЦИЯ НЕСТАЦИОНАРНЫХ ГИДРОДИНАМИЧЕСКИХ ИССЛЕДОВАНИЙ В СИСТЕМЕ «ПЛАСТ-СКВАЖИНА» МЕТОДОМ ИНТЕГРИРОВАННЫХ МОДЕЛЕЙ
К.С. Гаврилов, В.Л. Сергеев
Томский политехнический университет E-mail: [email protected]
Рассматриваются проблемы повышения информативности методов интерпретации нестационарных гидродинамических исследований в системе «пласт-скважина», и предлагается подход к их решению, основанный на использовании метода интегрированных моделей и метода адаптивной интерпретации. Приводятся примеры, показывающие преимущества предлагаемого подхода.
Ключевые слова:
Интерпретация, идентификация, гидродинамические исследования скважин, априорная информация, нефтяные пласты.
Key words:
Interprétation, identification, hydrodynamic analysis of oii wells, a-priori information, oilpools.
Введение
Нестационарные гидродинамические исследования скважин (ГДИС) являются наиболее информативным методом определения коллекторских и геометрических параметров пластов нефтяных и газовых месторождений. Они обеспечивают получение важнейших параметров залежей, на основании которых осуществляются процессы добычи нефти, составляются технологические проекты разработки месторождений, создаются геолого-техно-логические модели процессов нефтегазодобычи.
Существенным недостатком традиционных методов идентификации и интерпретации ГДИС является их затратный характер, поскольку интерпретация результатов производится после завершения заранее спланированных исследований, что часто приводит к продолжительным простоям скважин, значительной потере добычи углеводородов [1].
Однако в настоящее время, в связи с возможностью получения информации в режиме реального времени стационарными измерительными системами, требуется иная технология, позволяющая определять фильтрационные параметры и энергетическое состояние залежей в процессе гидродинамических исследований, не планируя заранее время их завершения [2].
С другой стороны, для повышения информативности, обеспечения устойчивости и повышения точности методов интерпретации ГДИС, как отмечалось рядом авторов, требуется интеграция информации о параметрах пласта и скважины [2-4]. Одним из перспективных способов интеграции информации, зарекомендовавшим себя при решении ряда практических задач мониторинга и управления разработкой нефтяных месторождений, является технология интегрированных моделей, которая позволяет объединить исходные данные, дополнительные априорные сведения, накопленный опыт и знания в единую систему моделей, что обеспечивает устойчивость оценок и существенно повышает их точность.
В данной работе на основе метода интегрированных моделей [4] и метода адаптивной идентификации [3] предлагаются и исследуются модели и алгоритмы, позволяющие проводить интерпретацию нестационарных ГДИС в процессе поступления информации, наиболее полно учитывать дополнительные априорные данные и экспертные оценки о параметрах пласта и скважины
Технология адаптивной интерпретации гидродинамического исследования скважин в системе «пласт-скважина»
Основой алгоритмов адаптивной интерпретации нестационарных ГДИС в системе «пласт-скважина» является стохастическая интегрированная система моделей вида:
Р,п = (/п ,ап , О) + €п ,
■ Рп, = Р (Т,а„,О)+П„, Га, = Га„ +vnl,I = 1,р,
V = V (Т, Оап) + еп, п = 1~П„, (1)
где РР3(а„,дп) - векторы фактических значений забойных давлений в текущий момент времени 4 и вычисленные на основе модели КВД /„ заданной с точностью до параметров а„=(а1„,а2„,а„т); ¥„, Г(/„,Оа) - экспертная оценка накопленной жидкости в стволе скважины после ее остановки, вычисленная на основе модели ¥(Т,Оа); Рш, Га, - векторы дополнительных априорных данных о пластовом давлении и фильтрационных параметрах пласта (полученных из ранее проведенных ГДИС, геофизических исследований скважин, лабораторных исследований керна и т. п.); Г, О - диагональные индикаторные матрицы нулей либо единиц (символ «0» используется для указания отсутствующих элементов векторов параметров а, а); Т -время максимальной длительности исследований скважины, гарантированное для восстановления в скважине пластового давления (7»34); £„, п, V,, е - векторы случайных величин, представляющих погрешности измерений забойных давлений, дебита притока продукции в скважине, ошибки допол-
нительных данных и экспертных оценок, неточность моделей и т. п.
Технологию адаптивной идентификации и интерпретации ГДИС по кривой восстановления давления (КВД), основанную на методе интегрированных моделей, можно представить как процесс принятия решений, содержащий этапы:
1) формирование исходных данных, дополнительных априорных сведений, экспертных оценок, выбор модели КВД в рамках системы моделей (1)
{Р*п,Рж, а1,1 = 1,p,Vn,р (a„,qn),n = 1,nk}; (2)
2) выбор показателя качества системы моделей (1) и решение оптимизационных задач вида
аn (ßn) = arg min ф(«п , ßnX n =1 nk,
ßl = arg min J0(fl'n(ß„)),
(3)
(4)
Ф = P
,- P3 (tn ,an, q)||^ + e \pm - P3 (T,a, q
+ß2 \Га, - Га + ß3 V - V(T, Ga
(6)
оптимизационная задача (3) с использованием метода Гаусса-Ньютона сводится к последовательно-
му решению системы линейных алгебраических уравнений
а = <1 + Ь,п^а‘п1’ * = 1>2Д...,
[Аг-1 -Да;-1 = В-1,
где запись ||2|Уозначает квадратичную форму ;
Аі-1 = (ОзШБз + РРТШВШ + ^0^ + &/))
в‘-1 = (Втзжез + + в о: е + в Да)1,
(7)
D =
Г дРз (t, ,а, q) . — . -—^
3 ' ■, i = 1,n , j = 1,m
да
(
D =
D =
дРз (T,а) . —
3 , J = 1,m
да
dV (T, Ga)
da
, J = 1, m
где запись arg min/x) означает точку минимума x* функции /(x)(/xQ=min/(x)); Ф(а,Д)= =Ф(/0(ап),АЛ(а),к=1,р) - комбинированный показатель качества системы моделей КВД (1), представляющий заданную функцию (функционал) Ф от частного показателя качества /0(а) модели КВД и от взвешенных весами ßn частных показателей качества /к(а) моделей дополнительных априорных данных и экспертных оценок;
3) анализ точности полученных приближений и принятие решения о прекращении либо о продолжении исследований с возвратом к этапу (2), либо (1).
Решение о прекращении исследований может быть принято на основе визуального анализа графика зависимости оценок параметров пласта и скважины (3) от времени либо по критерию стабилизации полученных приближений, где, например, за момент их завершения пк принимается то значение n, при котором выполняется неравенство
|(«l(n-/) - «jn ) / «nj | ^ Sj , j = 1,m> i = 1 2 3, ...
(Sj - заданная точность) (5)
Предлагаемая технология (2)-(5) позволяет синтезировать достаточно широкий спектр алгоритмов адаптивной идентификации и интерпретации ГДИС для линейных, нелинейных и непараметрических моделей КВД и моделей объектов аналогов [4].
Отметим, что для нелинейной модели КВД при выборе комбинированного показателя качества в виде суммы частных квадратичных показателей качества
- матрица и векторы частных производных от моделей КВД, пластового давления и накопленной добычи жидкости;
е = (Р* - Р (а)), е„, = (Р - Р (Т,а)),
Да = (Га - Га), е* = (V* - V(Т, Са))
- векторы и переменные невязок; Р’м, Га, - ус-
редненные значения дополнительных априорных данных и экспертных оценок пластового давления, фильтрационных параметров пласта и накопленной добычи продукции; W=diag(w(n), п=1,пк) - диагональная матрица весовой функции м>(1), определяющая вес забойного давления Р'з п в текущий момент времени п Д, вг, в3 - управляющие параметры. Для получения (7) достаточно взять частные производные по параметрам ап от комбинированного функционала (6) и приравнять их к нулю.
Следует также отметить, что задача (4) по определению оптимальных значений вектора управляющих параметров Д* не имеет аналитического решения и решается методами последовательных приближений.
Результаты идентификации и интерпретации
Результаты идентификации и интерпретации ГДИС скважины № 335 однородно-пористого нефтяного пласта месторождения Тюменской области на неустановившихся режимах фильтрации по КВД приведены на рис. 1-3 и в таблице. В качестве модели КВД использована зависимость Рз(0, позволяющая учитывать приток продукции в скважину после ее остановки [1]:
Р3 ($) = Р&) +
1 :) • йт, (8)
+ — J(q - q(т)) exp(- ,
4п<г t0 t-т к (t-т)
где д0 - дебит скважины до ее остановки; д(т) - модель дебита притока продукции в скважину после ее остановки: Р3(/0) - забойное давление в момент остановки скважины /0; а, к1=4х/г2щ - гидропроводность нефтяного пласта и комплексный параметр;
2
V X - приведенный радиус скважины и пъезопро-водность пласта.
На рис. 1 приведены значения забойного давления р„, ат, и дебита притока продукции в скважине д,, м3/сут.
Время,ч
Рис. 1. Значения забойного давления (1) идебита притока (2)
Обработка результатов ГДИС проводилась с использованием программных комплексов PanSystem, Saphir, адаптивного метода интегрированных моделей (7), адаптивного метода наилучшего совмещения, который, по аналогии с (7), сводится к последовательному решению системы линейных уравнений вида
a i = 1,2,3,...,
\(DlWD3 )i-1 -ЛаГ1 = (DT3We3 )M, (9)
и адаптивного метода интегрированных моделей с учетом априорной информации о фильтрационных параметрах пласта (АМИМ_1)
< =<Х + КпЛап "', 1 = 1,2,3,...,
<( dTwd3 + pdTdv )i-1 Ла\-1 =
^ = ( DT3We3 +pDfev ) Л (10)
который следует из (7) при вг=вз=0-
Оценки гидропроводности trn'=a{(P), комплексного параметра kp'(n)=a2'(P) и пластового давления
P* (n) = P(t0) +
+j(q° q(r)) exp(- —1—-) • dz
4n°n t„ n-T kp(n-t)
получены с использованием стратегии «скользящего интервала», где в обработке участвуют данные из интервала [п„+Лп+г], 1</<nk-bn, nH=0, Лп=5, что достигается соответствующим выбором значений весовых функций w(n), n=1,nk в матрице W(7), (9), (10). Здесь i - номер текущего интервала обработки; nH, Лп - количество измерений забойного давления в начальном участке КВД [t0,tH] и, соответственно, в интервале обработки. Оценки управляющих параметров ft=(P\,Pi,P) получены путем решения оптимизационной задачи
РП = argmin ||p;_n - P3 (tn rítfnX
методом деформированного многогранника [5].
В качестве дополнительных априорных сведений использовались данные о гидропроводности -=1 Дсм/сП, комплексном параметре ^=1 10-4с-1, пластовом давлении Рш=245 ат, и накопленной жидкости в стволе скважины после ее остановки ¥=2 м3 с последующей их корректировкой по схеме приведенной в [6].
Таблица. Результаты обработки КВД скважины № 335
Метод Время, ч Пластовое давление, ат Гидропро- водность, Дсм/сП Комплексный параметр, c-1,10-4
«PanSystem» 98 243,6 0,252 0,47807
«Saphir» 98 242,8 0,246 0,48134
АМНС 25 260,1 0,223 0,40235
50 257,9 0,245 0,74412
75 256,6 0,247 1,08129
98 254,6 0,268 2,31790
АМИМ_1 25 258,6 0,231 0,48 573
50 258,0 0,234 0,48409
75 257,2 0,233 0,48340
98 256,1 0,232 0,48314
АМИМ 25 246,9 0,252 0,48716
50 245,4 0,251 0,48407
75 244,6 0,253 0,48334
98 243,6 0,250 0,48309
На рис. 1-3 и в таблице приведены оценки пластового давления, гидропроводности и комплексного параметра, полученные адаптивным методом интегрированных моделей в системе «пласт-скважина» (АМИМ), адаптивным методом интегрированных моделей с учетом информации о фильтрационных параметрах пласта (АМИМ_1) и адаптивным методом наилучшего совмещения (АМНС) за разные периоды исследований.
Время,ч
Рис 2. Оценки пластового давления сучетом притока; Метод: 1) АМИМ; 2) АМИМ_1; 3) АМНС
Время,ч
Рис. 3. Оценки гидропроводности пласта; Метод: 1) АМИМ; 2) АМИМ_1; 3) АМНС
Из рис. 2, 3 и таблицы видно, что оценки параметров нефтяного пласта, полученные адаптивным методом интегрированных моделей в системе «пласт-скважина», устойчивые и более точные по сравнению с оценками, полученными адаптивным методом наилучшего совмещения и адаптивным методом интегрированных моделей с учетом информации о фильтрационных параметрах пласта.
Выводы
1. Предложены модели и алгоритмы адаптивной идентификации и интерпретации нестационарных гидродинамических исследований скважин в системе «пласт-скважина», основанные на технологии метода интегрированных моделей и позволяющие наиболее полно учитывать дополнительную априорную информацию о фильтрационных параметрах, энергетическом состоянии нефтяных пластов и накопленной продукции в скважине после ее остановки.
СПИСОК ЛИТЕРАТУРЫ
1. Шагиев Р.Г. Исследование скважин по КВД. - М.: Наука, 1998. - 304 с.
2. Ипатов А.И., Кременецкий М.И., Гуляев Д.Н. Современные технологии гидродинамических исследований скважин и их возрастающая роль в разработке углеводородов // Нефтяное хозяйство. - 2о09. - № 1. - С. 52-55.
3. Сергеев В.Л., Аниканов А.С. Метод адаптивной идентификации гидродинамических исследований скважин с учетом априорной информации // Известия Томского политехнического университета. - 2010. - Т. 317. - № 5. - С. 50-52.
2. На примере гидродинамических исследований скважины нефтяного месторождения Тюменской области показано, что метод адаптивной интерпретации в системе «пласт-скважина» дает более устойчивые и точные оценки параметров однородно пористого нефтяного пласта, сокращает время простоя скважины по сравнению с традиционным методом наилучшего совмещения, технологиями, реализованными в программных комплексах «PanSy-stem» и Saphir, и технологией адаптивной идентификации, где учитывается априорная информация только о фильтрационных параметрах пласта.
3. Наиболее полный учет дополнительной априорной информации о фильтрационных параметрах и энергетическом состоянии пласта и скважины позволяет повысить качество алгоритмов интерпретации, сократить простои скважин и увеличить добычу нефти.
4. Сергеев В.Л. Интегрированные системы идентификации. -Томск: Изд-во Томского политехнического университета, 2011. -198 с.
5. Пантелеев А.В., Летова Т.А. Методы оптимизации в примерах и задачах. - М.: Высшая школа, 2002. - 544 с.
6. Кемерова П.А., Сергеев В.Л., Аниканов А.С. Адаптивная идентификация и интерпретация нестационарных гидродинамических исследований с учетом притока продукции в скважине // Известия Томского политехнического университета. - 2011 -Т. 319. -№ 5. - С. 43 -46.
Поступила 21.04.2012 г.