© Н.С. Акелян, С.Б. Бекетов, 2011
УДК 622.245+622.279.7
Н.С. Акелян, С.Б. Бекетов
ОСОБЕННОСТИ ПОСТРОЕНИЯ МОДЕЛЕЙ ПРОГНОЗА ЭФФЕКТИВНОСТИ ГЕОЛОГОТЕХНИЧЕСКИХ МЕРОПРИЯТИЙ ПРИ РАЗРАБОТКЕ МЕСТОРОЖДЕНИЙ НЕФТИ И ГАЗА
Приведено описание методического подхода к оценке прогнозных показателей эффективности проведения геолого-технических мероприятий с целью повышения производительности скважин нефтяных и газовых месторождений.
Ключевые слова: нефть, газ, скважина, прогноз, модель.
Одной из основных проблем разработки нефтяных и газовых месторождений является прогнозирование прироста добычи нефти при проведении различных мероприятий по увеличению нефтегазоотдачи. Особое значение такой прогноз приобретает на поздней стадии разработки месторождений, когда в продукции скважин присутствует большое количество воды, их дебиты снижены по сравнению с начальными, а в продуктивных пластах происходят снижающие фильтрацион-но-емкост-ные свойства техногенные процессы, связанные с отбором флюидов. Несмотря на значительные успехи нефтепромысловой науки, достигнутые в последние годы, информация о текущем состоянии разработки эксплуатируемых пластов обычно остается неполной, что затрудняет создание и постоянное обновление многомерных детерминированных моделей фильтрации продуктивных пластов. Поэтому при прогнозировании эффективности проведения мероприятий по увеличению нефте-газоотдачи на практике специалисты добывающих предприятий пользуются в основном специально разработанными
кривыми вытеснения, а также корреляционным анализом.
Таким образом, проблемы планирования геолого-технических мероприятий (ГТМ) по увеличению производительности скважин нефтегазодобывающих предприятий, являются достаточно сложными и многоаспектными, особенно на поздней стадии разработки месторождений.
Как показывает опыт, в практических задачах не встречаются случайные величины, распределения которых точно соответствовали бы теоретическим распределениям. Последние являются математическими моделями реальных распределений. Подбор таких моделей и анализ их адекватности моделируемым случайным величинам, является одной из основных задач математической статистики, которая, в свою очередь, сводится к проверке предположений (гипотез) о виде модели распределения и о его параметрах.
Из-за наличия неучтенных факторов и неоднозначного действия учтенных, функция зависимости прироста добычи нефти является функцией лишь в статистическом смысле:
Характер изменения уровней добычи нефти и жидкости на различных стадиях разработки месторождений
описываемая ею математическая зависимость проявляется только в общем и среднем в массе наблюдений. Соответственно и аппаратом исследования производственных функций должны служить вероятностно-статис-тические методы.
При прогнозировании эффективности ГТМ, целесообразно использовать многофакторные функции, позволяющие измерить характер и силу совместного, комбинированного влияния нескольких показателей-факто-ров на величину изучаемого результативного показателя. Общий вид многофакторной функции может быть записан в виде:
у= ^ХЬ..., X!, ...,хп), (1)
где у - результативный показатель, Х1,..., X!, .,Хп - показатели-факторы.
Как известно, различают статические и динамические виды функции. В статических функциях не учитывается время как фактор, изменяющий основные характеристики изучаемой зависимости. Динамические производственные функции включают фактор времени: время может фигурировать в них как самостоятельная переменная величина, влияющая на результативный показатель, параметры и показатели-факторы также могут, рассматриваться как функции времени.
В нефтегазодобывающей промышленности важнейшую роль играет особенность изменения во времени основных эксплуатационных показателей (снижение дебита скважин, увеличение обводненности продукции, падение пластового давления, уменьшение запасов нефти и газа в пласте, выпадение в пласте солей, разрушение скелета породы и т.д.), что определяет ф актическую нео бходи -мость использования функции времени, как одного из определяющих факторов этапа разработки месторождений.
Как уже отмечалось, процесс разработки нефтяных и газовых месторождений носит стадийный характер. Условно можно выделить четыре стадии, выделенные на примере одного из газовых месторождений на юге России (рисунок). На первой стадии (участок I), когда происходит активное разбуривание и обустройство месторождения, ввод скважин и промысловых сооружений в эксплуатацию, добыча углеводородов растет, что обусловлено в значительной степени скоростью строительства скважин и обустройства месторождения которая зависит от работы буровых и промыслово-строительных подразделений.
Вторая стадия (участок II) характеризуется максимальной добычей углеводородов. Этот период соответствует «расцвету» эксплуатации месторождения. Скважины работают еще без осложнений, наземные коммуникации не нуждаются в ремонте. Прибыль, получаемая на этом этапе, как правило, является максимальной.
Третья стадия (участок III) характеризуется падением добычи углеводо-
родов и значительным ростом обводненности продукции скважин.
На четвертой стадии (участок IV) наблюдаются сравнительно мед-лен-ное, постепенное падение добычи углеводородов, дальнейший рост обводненности продукции скважин, вывод скважин из эксплуатационного фонда по причине их обводненности. В итоге четвертая стадия завершается полным прекращением добычи углеводородов и ликвидацией всей промысловой инфраструктуры на месторождении. Четвертую стадию называют поздней или завершающей стадией разработки.
Следует особо отметить, что развитие широкого спектра методов повышения нефтегазоотдачи пластов, применение эффективного воздейст-вия на продуктивные отложения позволяет на III и IV стадиях разработки месторождений стабилизировать до-бычу углеводородов, или даже до-биться увеличения добычи, увеличить коэффициент извлечения углеводо-родов из пласта, продлить период эксплуатации месторождения. Ис-пользование того или иного метода воздействия на пласт должно предва-рительно апробироваться на моделях, в т.ч. и математических с целью про-гноза фактически получаемых пара-метров на скважинах.
Количественное прогнозирование основано на переносе прошлых данных на будущее, то есть на предположении, что условия, породившие полученные ранее данные, мало отличаются от будущих условий. Исключение составляют те факторы, которые точно распознаны моделью прогнозирования, что в нефтегазодобыче встречается относительно редко. Можно выделить следующие этапы построения модели прогнозирования [1]:
Остановимся подробнее на важнейшем этапе - построение и оценка прогнозной модели.
Вначале исследуется модель, которая включает все возможные, поддающиеся изучению факторов влияя-ющих на цель исследования. Затем, проводится упрощение начальной модели, которая по принципу "простоты" будет являться "наилучшей" (практичной) интерпретацией цели исследования. Простые модели легче тестировать во время кросспроверки и повторного использования. Кроме того, простые модели намного легче использовать в будущем с точки зрения стоимости получения отклика и значений предикторов.
Прогнозирование результатов проведения того или иного вида ГТМ разрабатываемых месторождений нефти и газа с высокой достоверностью возможно при решении ряда задач, основными из которых являются следующие.
1. Определение наиболее значимых факторов, влияющих на прирост уровня добычи нефти и газа при проведении рассматриваемого ГТМ на данном объекте разработки.
2. Определение величины изменения производительности скважины (определение эффективности) в результате проведения ГТМ.
3. Определение вида функции изменения производительности скважи-ны и накопленной добычи углеводородов как функции от времени (соответственно и определение времени действия положительного эффекта от проведения ГТМ).
Для решения первых двух из поставленных задач могут служить такие математические методы и модели как вероятностные и статистические методы, имитационное моделирование и др.
Таблица 1
Основные этапы построения модели прогнозирования.
№ этапа Описание Примечание
I Сбор исходных данных и их проверка. Определение набора возможных независимых переменных, содержащих новую информацию о зависимой переменной. Выявление переменных, которые могут оказаться неадекватными данному прогнозу Независимая переменная исключается из рассмотрения, если она не имеет существенного значения для решения поставленной задачи или характеризуется значительными ошибками измерения. Проводится проверка на мультиколлинеарность - дубляж других независимых переменных. Если между независимыми факторами (предикторами) существует высокая корреляция (линейный коэффициент корреляции свыше 0,7-0,8), то нельзя определить их изолированное влияние на результирующий (прогнозируемый) показатель. Следует учесть, что данный подход не учитывает возможность ложной корреляции.
II Редукция или уплотнение исходных данных. Этап необходим в связи с тем, что может быть собрано как избыточное так и недостаточное число исходных данных.
III Построение и оценка прогнозной модели. Данный этап заключается в подборе модели прогноза, соответствующей особенностям собранных исходных данных с точки зрения минимизации ошибок прогноза. При этом, в процессе выбора прогнозной модели необходимо использовать методы качественной оценки.
IV Экстраполяция модели. Модель проходит тестирование на уже имеющихся данных, то есть применяется апробация модели - прогнозирование на сравнительно недавно полученных данных, для которых результаты уже известны, что позволяет проанализировать наблюдаемые ошибки.
Для определения наиболее значимых факторов, влияющих на прирост уровня добычи, вначале отбираются факторы на основе данных теоретического исследования. Далее отобранные факторы подвергаются проверке существенности их влияния на изучаемый показатель, как правило, такая проверка включает анализ матрицы парных корреляций, частных корреляций, проверку существенности (значимости) коэффициентов регрессии на основе ^критерия, анализ остатков (отклонений) и т. д. Особенностью множественной регрессии и корреляции является необходимость различать случаи корреляционной множественной связи, когда переменные х1, х2 ... хп являются случайные величинами; регрессионной, если переменные х1, х2 ... хп - неслучайными величинами, а
также смешанный случай, когда некоторые из переменных - случайные величины, другие - неслучайные. В случае корреляционной зависимости следует вычислять и интерпретировать коэффициенты корреляции, при регрессионной зависимости это не имеет смысла, а при наличии как случайных, так и неслучайных переменных коэффициенты корреляции следует вычислять только между случайными переменными [2]. Поскольку при прогнозировании эффективности ГТМ в модель включаются как случайные (пористость продуктивных пластов, коэффициент глинистости, вязкость пластовой нефти и др.), так и неслучайные величины (расход реагентов на 1 м перфорации, давление закачки и др., определяемые в зависимости от проводимого ГТМ), необходимо проводить про-
верку значимости факторов с использованием различных инструментов, учитывающих «смешанный» характер набора факторов.
Задачу определения основных показателей эффективности ГТМ (увеличения производительности скважины, времени действия положительного эффекта, накопленной добычи углеводородов) целесообразно решать с помощью такого инструмента выявления и анализа связи между несколькими переменными как множественная регрессия. Множественная регрессия представляет собой регрессию результативного признака с двумя или большим числом независимых переменных вида у= Дх^,..., х1, .. .,хп). Вид уравнения множественной регрессии может быть как линейный, так и нелинейными, в зависимости от первоначальных условий. К основным видам нелинейной регрессии относятся:
а) Полиномиальная регрессия, содержащая «главные эффекты» (эквивалентно тому что, значения различных предикторов независимы и входят в модель в первой степени) и эффекты высшего порядка для непрерывных предикторов (факторов), но не содержит взаимодействия эффектов. Например, полиномиальная регрессия второй степени для трех непрерывных предикторов р, q и г будет включать главные эффекты для р, q и г и их квадратичные (то есть, второго порядка) эффекты, но эта модель не будет содержать 2-фактор-ные эффекты или 3-факторное взаимодействие р на q и на г:
у = Ьо + Ь1р + Ь2р2 + bзq +
+Ь^2 + Ь5г + Ь6г2
б) Факторная регрессия. Для упрощения факторную регрессию определяют как функции, в которых присутствуют все возможные произведения непрерывных предикторов. Например, полная факторная регрессия для двух
непрерывных предикторов р и q будет включать главные эффекты (то есть эффекты первого порядка) для р и q, а также 2-факторное взаимодействие р по q, которое представлено произведением значений р и q: у = Ьо + Ьф + b2q + bзp*q
в) Регрессии квадратичной поверхности отклика - это смешанный тип регрессии со свойствами функций полиномиальной регрессии и факторных регрессионных функций. Регрессии квадратичной поверхности отклика содержат все эффекты из функций полиномиальной регрессии вплоть до 2 степени, а также 2-факторные взаимодействия предикторов. Уравнение регрессии с 3 непрерывными предикторами р, q и г:
у = Ь0 + Ь1р + Ь2р2 + + Ь^2 +
+Ь5г + Ьбг2 + b7p*q + Ь8р*г + b9q*г
г) Регрессионные модели с точками разрыва (кусочно-линей-ная регрессия). Данный вид модели следует использовать в тех случаях, когда вид зависимости между предикторами и переменной различается в разных областях значений независимых переменных. Общий вид функции кусочно-линейной регрессии имеет вид:
у = (Ьо1 +ЬП *Х1 +...+Ьт1 *хт )*(у <Ьп) + +(Ьо2 +Ь12 *Х1 +...+Ьт2 *Хт )*(у > Ьп ).
Такие функции нелинейной регрессии как степенная, показательная, экспоненциальная, логарифмическая удобно применять в тех случаях, когда количество предикторов (факторов) мало, в противном случае, обработка с использованием данных функций окажется громоздкой [3].
После составления модели множественной регрессии, проводится оценка регрессионных параметров. В множественной регрессии оценивание параметров регрессии происходит подбором коэффициентов, минимизирующих дис-
персию остатков (сумму квадратов остатков). Любые отклонения наблюдаемых величин от предсказанных означают некоторые потери в точности модели, например, из-за случайного шума (ошибок). Поэтому цель метода наименьших квадратов (на котором в большинстве случаев основывается построение регрессионных моделей) заключается в минимизации функции потерь. При достижении этой функцией минимума, получаем оптимальные оценки для параметров (свободного члена, коэффициентов регрессии), т.е. оценки по методу наименьших квадратов [3].
Следует также обратить особое внимание на так называемые выбросы (т.е. экстремальные наблюдения). Выбросы могут вызвать серьезное смещение оценок, "сдвигая" линию регрессии в определенном направлении и тем самым, вызывая смещение регрессионных коэффициентов. Часто исключение всего одного экстремального наблюдения приводит к совершенно другому результату. Важно различать действительные ошибки (настоящие выбросы) и корректные, но резко отличающиеся наблюдения.
В нефтегазодобывающей промышленности приходится иметь дело с трудно прогнозируемыми параметрами и отличить действительные ошибки от корректных значений часто бывает очень сложно. Можно просто исключить то или иное наблюдение из выборки и анализировать оставшиеся более согласованные между собой данные. Но, к сожалению, часто трудно найти убедительное подтверждение тому, что выбросы не соответствуют изучаемому процессу. Например, это касается дополнительного прироста добычи углеводородов после проведения какого-либо ГТМ. Исключив какое-либо экс-
тремальное значение этого прироста, можно потерять очень важное для прогнозирования значение. В области нефтедобычи на вопросы учитывать или нет то или иное наблюдение, как правило, не бывает однозначных ответов. При этом необходимо учитывать достаточно большой процент неудачных (безуспешных) воздействий на продуктивные пласты, так называемых бросовых работ, среди которых, конечно, наблюдаются явные технологические, а иногда и геологические браки.
После оценивания регрессионных параметров необходимо провести оценку прогнозной модели в целом. Рассмотрим некоторые методы проверки пригодности модели:
Объясненная доля дисперсии - отношение полной дисперсии зависимой переменной к доле дисперсии, приходящейся на остатки. Данное значение служит индикатором степени подгонки модели к исходным данным.
Метод максимального правдоподобия (хи-квадрат). Вычисление основано на методе максимального правдоподобия, т.е. максимизации функции правдоподобия. На практике статистика максимального правдоподобия хи-квадрат очень близка по величине к обычной статистике Пирсона хи-квадрат.
Также полезным бывает использование диаграммы рассеяния наблюдаемых и предсказанных значений. Если модель хорошо соответствует данным, можно ожидать, что точки расположатся вдоль прямой линии.
Кроме того для сравнения двух выборок, в нашем случае фактических данных и данных, полученных с помощью модели, используют икритерий (t-критерий является наиболее часто используемым методом обнаружения различия между средними двух выборок). и критерий может применяться, если пе-
ременные нормально распределены (внутри выборки), а дисперсии наблюдений в выборках не слишком различны. Если условия применимости /-критерия не выполнены, следует использовать непараметрические альтернативы /критерия: критерий серий Вальда-
Вольфовица, и критерий Манна-Уит-ни и критерий Колмогорова-Смир-нова и другие [4].
Проблемы прогнозирования технологических результатов проводимых гео-
1. Галиуллин Ф.Р., Карпов В.Г., Котов Д.В. Вопросы совершенствования методов прогнозирования эффективности мероприятий по увеличению нефтеотдачи пластов // Проблемы экономики и управления. - 2004. - № 3. - С. 2428.
лого-технологических мероприятий на сегодняшний день для нефтегазодобывающих предприятий, разрабатывающих нефтяные месторождения на поздней стадии их эксплуатации, являются достаточно сложными и многоаспектными. Для решения данных проблем и совершенствования инструментов прогнозирования с целью последующей оптимизации технологических показателей рекомендуется использование корреляционно-регрессионных зависимостей.
--------------- СПИСОК ЛИТЕРАТУРЫ
2. Красс М.С., Чупрынов Б.П. Математика для экономистов. - СПб: Питер, 2005. С. 464.
3. STATISTICA (Версия 6.1), Электронное руководство.
4. Справочник по математическим методам в геологии.-М.: Недра, 1987. С. 335.ЕШ
— Коротко об авторах --------------------------------------------------------------------1
Бекетов С.Б. - доктор технических наук, профессор СевКавГТУ, КубГТУ (г. Краснодар), почетный работник науки и техники РФ. bsb.gt@rambler.ru Акелян Н.С. - инженер 1-й категории ООО «НК «Роснефть»-НТЦ».
А
----------------------------------- ДИССЕРТАЦИИ
ТЕКУЩАЯ ИНФОРМАЦИЯ О ЗАЩИТАХ ДИССЕРТАЦИЙ ПО ГОРНОМУ ДЕЛУ И СМЕЖНЫМ ВОПРОСАМ
Автор Название работы Специальность Ученая степень
КАРАГАНДИНСКИЙ ГОСУДАРСТВЕННЫЙ ТЕХНИЧЕСКИЙ УНИВЕРСИТЕТ
ИСАБЕКОВ Ерен Туякович Технологические схемы ведения очистных и подготовительных работ при разработке выбросоопасных угольных пластов 25.00.22 к.т.н.