Научный инжиниринг как основа процессов моделирования при разработке месторождений

Хасанов М.М.; Ситников А.Н.; Пустовских А.А.; Рощектаев А.П.; Исмагилов Н.С.; Падерин Г.В.; Шель Е.В.

МАТЕРИАЛЫ КОНФЕРЕНЦИИ

DOI: https://doi.Org/10.18599/grs.2018.3.142-148

УДК 550.8:622.276

Научный инжиниринг как основа процессов моделирования

при разработке месторождений

М.М. Хасанов, А.Н. Ситников, А.А. Пустовских*, А.П. Рощектаев, Н.С. Исмагилов, Г.В. Падерин, Е.В. Шель

ООО «Газпромнефть НТЦ», Санкт-Петербург, Россия

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

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

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

Третий пример демонстрирует применение теории размерностей к задаче моделирования гидроразрыва пласта (ГРП) для определения значимости тех или иных параметров на дизайн ГРП. С помощью обезразмеривания зависимости длины трещины от объёма закачки, получается эмпирическая формула для длины трещины ГРП от его параметров, по которой определяется уровень их значимости.

Ключевые слова: геологическое моделирование; геостатистика; геомеханика; гидроразрыв пласта; гидродинамическое моделирование; метод источников

Для цитирования: Хасанов М.М., Ситников А.Н., Пустовских А.А., Рощектаев А.П., Исмагилов Н.С., Падерин Г.В., Шель Е.В. (2018). Научный инжиниринг как основа процессов моделирования при разработке месторождений. Георесурсы, 20(3), Ч.1, с. 142-148. DOI: https://doi.org/10.18599/grs.2018.3.142-148

Задачи контроля и управления технологическими процессами разработки месторождений часто приводят к необходимости моделирования геофизических полей, процессов, связанных с фильтрацией в пласте и движением по трубам многофазных смесей, геомеханических задач в процессе бурения и при проведении гидроразрыва пласта, и многих других. Традиционно, описание процессов нефтегазодобычи проводится на основе дифференциальных уравнений движения жидкостей и газов в пористых средах и трубах. Однако такой подход не позволяет описать многие существенные свойства пласта. Как всякие большие системы, объекты нефтегазодобычи требуют использования целой иерархии моделей - от дифференциальных до интегральных, от детерминированных до адаптивных, способных описать не только различные уровни организации систем, но и взаимодействие между этими уровнями (Дьячук, 2003). Решением всех этих задач по иерархичному моделированию процессов нефтегазовой отрасли занимается научный инжиниринг (англ. SriencEngmeermg) - область науки, находящаяся на стыке системного инжиниринга, фундаментальных наук и кибернетики (Рис. 1).

В данной работе будут рассмотрены примеры подходов научного инжиниринга к решению задач, которые

* Ответственный автор: Алексей Анатольевич Пустовских E-mail: [email protected]

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

Геологическое моделирование: применение спектрального подхода для моделирования геофизических полей

Разработка спектрального подхода к задачам геологического моделирования была впервые представлена в работе (Байков и др., 2010), в которой авторами предложен

Рис. 1. Научный инжиниринг (SE) личных походов

на стыке областей раз-

SCIEN71FICAND TECHNICAL JOURNAL

GEDRESGURCES

способ, основанный на разложении каротажа на коэффициенты по ортонормированному базису в пространстве функций, интегрируемых с квадратом L2. Дальнейшее развитие метод получил в работах (Байков и др., 2012; Хасанов и др., 2015), в которых осуществлена более детальная проработка теоретических основ спектрального метода, а также представлены некоторые результаты расчётов, основанных на спектральном методе геологического моделирования. Стоит также упомянуть о том, что спектральный анализ каротажа получил применение для целей литолого-фаци-ального анализа (Хасанов и др., 2014).

Спектральный метод геологического моделирования -метод моделирования трёхмерных кубов геофизических свойств на основе скважинных данных. Принятая в спектральном методе геологического моделирования математическая модель представлена моделируемой областью й с Л3 стохастическим полем С(х,у,И), определенным на этой области и заданным на некотором вероятностном пространстве. Для такой модели скважинные данные для вертикальных скважин с фиксированной координатой по латерали (х*, у*) представляют собой случайный процесс, параметризованный переменной, характеризующей глубину С(х*,у*,И) = /(И). Если в моделировании участвуют N скважин с координатами (х,у ), на которых задано моделируемое свойство функциями /(И), то эти функции можно рассматривать как известные значения некоторой реализации стохастического поля С(х,у,И).

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

Известно, что в гильбертовом пространстве интегрируемых с квадратом функций (Ь.) существуют орто-нормированные базисы, которые позволяют разложить любую функцию, заданную в этом пространстве, в ряд по коэффициентам разложения, определяемым однозначно. Функции, определяющие моделируемое свойство, как функции с конечной энергией принадлежат этому пространству. Выберем в пространстве Ь2 базис многочленов Лежандра Р\И) который используем для разложения функций /(И):

Ш = (1)

где коэффициенты с', определяются как скалярное произведение в пространстве Ь. Набор коэффициентов разложения с', для каждого уровня разложения]определён в точках (х,у ), которые являются координатами скважин, и является известными значениями некоторой реализации стохастического поля с (х,у).

Для моделирования коэффициентов с' для всех (х,у) используется спектральный метод моделирования (Пригарин, Михайлов, 2005), основанный на теореме об интегральном представлении случайных процессов и полей.

Построение обусловленной реализации смоделированного поля с (ху) известными значениями с' осуществляется общепринятым для двухшаговых методов моделирования - строится разность кригингов по известным данным и данным, полученным в результате (Дюбрюль, 2002).

В результате осуществления вышеприведённых шагов получается набор реализаций стохастических полей коэффициентов разложения c (x,,y), которые обусловлены сква-жинными значениями cj. Для восстановления реализации всего поля G(x,y,h) в каждой точке моделируемой области D осуществляется суммирование базисных функций по смоделированным коэффициентам:

G{x,y,h)=YjCj{x,y)pjQi). (2)

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

Спектральный метод моделирования обладает существенными преимуществами перед классическими методами моделирования, такими как, например, последовательное гауссово моделирование. В частности, этот метод лишён главного ограничения классических методов - гипотезы о стационарности моделируемого свойства. Известно, что большинство геологических процессов, происходящих в природе, нестационарны. Существенным преимуществом метода является непараметрический периодограммный статистический анализ, вместо классически используемого параметрического вариограммного анализа, который не позволяет в полной мере оценить весь спектр изменчивости моделируемого свойства. Спектральный метод в силу алгоритма своего построения является принципиально параллелизуемым, что позволяют легко масштабировать вычисления. Кроме того, метод реализует построение модели без привязки к сетке (т.н. grid-free simulation), что позволяет осуществлять моделирование на сетках любой конфигурации и сложности, производить измельчение части сетки и уточнение модели на этой области и, наконец, параллельно генерировать реализации моделируемого свойства в различных областях модели. Прирост производительности спектрального моделирования на типовой рабочей станции геолога-модельера составляет от 2 до 4 раз по сравнению с классическими методами в зависимости от размера генерируемой модели.

Применение спектрального метода на практике продемонстрировало его способность хорошо воспроизводить геофизические поля в межскважинном пространстве и в неразбуренных областях, в особенности при сравнении с традиционными методами моделирования (Рис. 2) (Хасанов и др., 2014; Пригарин, Михайлов, 2005).

Гидродинамическое моделирование: метод источников в оптимизационных технологических расчетах

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

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

НАУЧНО-ТЕХНИЧЕСКИЙ ЖУРНАЛ

www.geors.ru ГЕОРЕСУРСЫ

ИЗ

Рис. 2. Сравнение реальной и синтетической кривых гамма каротажа. DGK, красная кривая, - реальный каротаж, SP_ DGK, синяя кривая, - синтетически каротаж, смоделированный спектральным методом, SGS_DGK, чёрная кривая, - синтетический каротаж, смоделированный методом SGS.

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

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

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

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

В представленной работе рассмотрен пример использования метода источников/стоков для решения задачи выбора оптимального времени отработки нагнетательных скважин (Ситников и др., 2015). Время отработки - это период работы скважины в режиме добычи перед переводом

ее в нагнетание.

Распределение давления от вертикальной скважины, для которой задана динамика жидкости, записывается в виде интеграла Дюамеля:

Wr,*)-/»+-!=-f-

'А0<

W! - г)

-dt

t- т

(3)

где ^

- коэффициент пьезопроводности, P - на-

чальное пластовое давление, ц - вязкость нефти, К - проницаемость, h - толщина пласта, С - общая сжимаемость, т - пористость, г - расстояние от источника, t - время расчета, т - переменная интегрирования.

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

Для кусочно-постоянного дебита жидкости выражение (3) трансформируется в уравнение (4):

где С - кусочно-постоянный дебит на каждом шаге, Рп -забойное давление, г - радиус скважины, п - количество

временных шагов, ЕЦ-х) = ( —еЬс - интегральная по-

х

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

Используя принцип суперпозиции (что обеспечивается линейностью уравнения пьезопроводности), можно рассчитать динамику нестационарных дебитов для системы из N скважин:

j = 1..К (5)

Для трещин ГРП конечной проводимости система уравнений (5) должна быть дополнена системой уравнений по расчету забойных давлений на источниках, имитирующих трещины.

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

пю^МгЛ. (6)

Под дисконтированной добычей с элемента разработки

SCIENT1FK AND TECHNICAL JOURNAL

GEGRESGURCES

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

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

Достоверность расчетов предложенным методом проверялась сравнительным анализом динамики добычи жидкости с результатами расчетов в коммерческом программном продукте и аналитическими зависимостями, полученными в статье (Хасанов и др., 2013). Апробация выполнялась на основе данных о строении одного из месторождений компании ПАО «Газпромнефть». Полученная зависимость безразмерной накопленной дисконтированной добычи от времени отработки нагнетательной скважины показана на рис. 3. Видно, что наибольшая дисконтированная добыча может быть получена при условии использования нагнетательных скважин в добыче на протяжении примерно девяти месяцев. Более ранний или поздний перевод скважин окажется экономически менее выгодным. Данный результат не только подтвердился сериями расчетов на полномасштабной гидродинамической модели, но и находится в согласии с характерным временем перевода скважин, которое было получено опытным путем при эксплуатации месторождения (Ситников и др., 2015).

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

[I ЕО IПО I ..Р ¿пп 2КО ]№ 4<й1

Т.сут

Рис. 3. Зависимость безразмерной дисконтированной добычи (¿^тахф^) от времени перевода скважины в нагнетание Т (Ситников и др., 2015)

Геомеханика: сравнение параметров ГрП в безразмерных переменных по дизайну и гидродинамическим исследованиям

Ввиду частого применения технологии гидроразрыва пласта (ГрП) на нефтяных и газовых месторождениях, имеется большой объем статистической информации по

проведенным операциям. По результатам обработки данной информации можно сделать вывод об эффективности проведённых ГрП, что позволит внести корректировки в дизайн и выработать дальнейшие рекомендации. однако данная задача осложнена тем, что при проведении гидродинамических исследований (ГДИС) на скважинах с ГрП, обнаруживается, что значения параметров трещины, в частности полудлина, существенно отличаются от тех, что запланированы по дизайну.

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

Из графика, приведенного на рис. 4 видно, что явной зависимости между длиной трещины и объемом закачки в общем случае не наблюдается. Причины, которые к этому привели, достаточно очевидны, и заключаются в том, что помимо объёма закачки сшитого геля (и связанной с ним массы проппанта) на длину, очевидно, могут влиять такие параметры пласта, как модули Юнга по разрезу, мощности пластов, сжимающие напряжения, перпендикулярные трещине, коэффициенты трещиностойкости, а также технологические параметры дизайна ГрП - реология жидкости, скорость закачки и концентрация проппанта.

Для анализа полученных данных требуется найти выборку скважин, которая была бы полностью идентична по всем параметрам ГрП кроме какого-то одного, после чего данную процедуру следовало бы повторить для всех параметров ГРП, чтобы установить степень влияния каждого из них. определив данные степени влияния, можно было бы установить наиболее значимые параметры для геометрии трещины ГРП. Это бы позволило предположить параметры, ошибка в которых могла привести к значительному систематическому расхождению между длиной трещины, запланированной по дизайну, и длиной трещины, полученной по гидродинамическим исследованиям скважины. На практике, однако, подобрать такие выборки трещин ГРП, которые не отличались бы по всем параметрам кроме какого-то одного, не представляется возможным.

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

0 1 00 200 300 .400 600 600

V, м3

Рис. 4. Зависимость длины трещин по дизайну и ГДИС от объема закаченной жидкости

НЮЧНО-ТЕХНИЧЕСЩЙ ЖУРНАЛ

www.geors.ru ГЕОРЕСУРСЫ

Для метода, предложенного в этой статье, это является необходимым условием.

Скорости закачки жидкости ГРП при дизайне также изменяются слабо, так что на большой разброс и образование «облака» на рис. 4 этот фактор не мог оказать решающее воздействие.

Таким образом, остаются в основном разнообразные геомеханические параметры, которые могут меняться вдоль пласта достаточно значительно (особенно мощность пласта). Для успешного статистического анализа требуется метод, который понизит размерность задачи по этим параметрам, и приведёт трещины ГРП в пластах с разными геомеханическими свойствами к одному «знаменателю». Таким методом является введение безразмерных параметров задачи о развитии трещины ГРП.

Развитие трещины по модели Planar3D в случае закачки неньютоновской жидкости описывается тремя законами (Хасанов и др., 2017):

- законом Гука;

- законом вязкого трения;

- законом сохранения массы.

Обезразмеривание уравнений проведем аналогично тому, как описано в работе (José I. Adachi et al., 2010) для случая трещины без проппанта; также будем считать, что утечки в пласт отсутствуют, модуль Юнга - однородный, а литология - трехслойная и симметричная. В результате получаем следующий список безразмерных параметров:

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

У =

k'E'2n+1Q" H3nàa2n+2 ,2п+2

V = V

кЕ

Qn

н3п+3д(72п+3 '271+1

L = L-

кЕ

Qn

С =

нЗп+1Аа2п+2

3n 1 1 СгН 2 +2АаП+2

Ï тт~

k'4'nQ2+2

К = I**

у H Лег

(7)

(8)

(9)

(10)

(11)

где Е =

1-v2

- модуль плоской деформации, E - модуль

Юнга, v - коэффициент Пуассона, n - показать поведения жидкости, k - коэффициент густоты потока, Q - расход жидкости, H- мощность пласта, Дст - контраст напряжений, L - полудлина трещины, V- объем закаченной жидкости, C- коэффициент утечек по Картеру, K - коэффициент трещиностойкости.

сравнив полученные безразмерные переменные с приведенными в статье (José I. Adachi et al., 2010) для модели Pseudo3D, можно увидеть, что они идентичны с точностью до констант. Единственным важным отличием является параметр g. Из работы видно, что если поделить мощность пласта на масштабный фактор длины из этой работы, то с точностью до константы получается безразмерный параметр g. Его физический смысл - отношение мощности пласта к достигаемым длинам трещины.

В итоге, обезразмеривание понижает размерность задачи на 5, так как модуль Юнга, коэффициент Пуассона,

контраст напряжений, мощность пласта и вязкость жидкости заменяются одним безразмерным параметром у (7).

Получение безразмерных параметров позволяет провести анализ уже проведенных работ ГРП, а именно, выявить некоторые закономерности при сравнении полудлин, полученных по дизайну и по ГДИС, в сопоставлении с общим объёмом закачки. Таким образом, необходимыми данными для анализа помимо значений полудлин трещин являются: мощность пласта, модуль Юнга и коэффициент Пуассона (входят в модуль плоской деформации), контраст напряжений в пласте, реология закачиваемой жидкости, расход жидкости и общий объём закачки жидкости ГРП. Прочностные свойства горных пород в данном анализе не учитываются. На рис. 5-7 представлены зависимости безразмерной длины от безразмерного объема для трех типов закачиваемой жидкости.

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

Приведенные выше результаты доказывают, что зависимость безразмерной длины трещины от безразмерного объема - степенная:

1 = А-Уа (12)

тогда подставив выражения (7) и (8) в (12), можно получить формулу, отражающую зависимость размерной длины трещины от объема закачки для произвольной степени а:

уЗп+1-(Зп+3)ад0.2п+2-(2п+3)а

Ь = А к,1-а(гП(1_а:)Е,2П+1-(2П+2)а ^ (13)

Из графиков (Рис. 5-7) можно заключить, что степень зависимости безразмерной длины трещины от объема закачки а находится в диапазоне 0.6-0.7. Приняв степень зависимости равной 0.6, а показатель поведения жидкости п равным 0.5, можно получить эмпирическую формулу для расчета полудлины трещины:

з

7 . Д<Т5 1

L = A ï 2 ï

HSk'SQSE'5

(14)

Данная формула позволяет оценить влияние ошибки (геомеханические параметры) или изменения входных параметров (технологические параметры) на изменение значения полудлины трещины (Табл. 1).

Параметр ^параметра Геомеханические параметры Технологические параметры

H E' Äff k' Q

+10% - 2% - 2% 6% - 4% - 2%

+20% - 4% - 4% 12% - 8% - 4%

+30% - 6% - 6% 18% - 12% - 6%

Табл. 1. Влияние изменения параметров на изменение длины трещины

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

SCIENT1FK AND TECHNICAL JOURNAL

GEDRESGURCES

Рис. 5. Зависимость безразмерной длины трещины от безразмерного объема: жидкость 1, м/р А; а = 0.69

Рис. 6. Зависимость безразмерной длины трещины от безразмерного объема: жидкость 2, м/р А; а = 0.63

Рис. 7. Зависимость безразмерной длины трещины от безразмерного объема: жидкость 3, м/р А; а = 0.66

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

литература

Байков ВА., Бакиров H.K., Яковлев A.A. (2010). Швые подходы в теории геостатистического моделирования. Вестник УГЛТУ, 37(2), с. 209-215.

Байков ВА., Бакиров H.K., Яковлев A.A. (2012). Maтемaтичеcкaя геология. Ъэм 1. Введение в геостатистику. M.-Ижевск: ИЕИ, Библиотека нефтяного инжиниринга, 228 с.

Дьячук A.H (2003). Hекоторые аспекты реконструкции систем сбора продукции скважин на поздних стадиях разработки нефтяных месторождений. Проблемы сбора, подготовки и транспорта нефти и нефтепродуктов, 62, с. 87-97.

Дюбрюль О. (2002). Использование геостатистики для включения в геологическую модель данных. EAGE: SEG. 295 c.

Пригарин CM., Mихaйлов ГА. (2005). Mетоды численного моделирования случайных процессов и нолей. Швосибирск: ИВMИMГ СО PAH, 259 c.

Ситников A.H., Пустовских A.A., Рощектаев A.H, Aнджyкaев Ц.В. (2015). Mетод определения оптимального времени отработки нагнетательных скважин в системе разработки месторождения. Нефтяное хозяйство, 3, c. 84-87.

Хасанов M.M., Белозеров Б.В., Бочков A.C, Ушмаев О.С., Фукс О^. (2014). Применение спектральной теории для анализа и моделирования фильтрационно-емкостных свойств нласта. Нефтяное хозяйство, 12, c. 60-64.

Хасанов M.M., Белозеров Б.В., Бочков A.C, Фукс О^., Tенгелиди Д.И. (2015). Aвтомaтизaция литолого-фациального анализа на основе спектральной теории. Нефтяное хозяйство, 12, c. 48-51.

Хасанов M.M., Mельчaевa О.Ю., Ситников A.H., Рощектаев A.H (2013). Динамика добычи из скважин с гидроразрывом пласта в экономически оптимальных системах разработки. Нефтяное хозяйство, 12, c. 36-39.

Хасанов M.M., Падерин Г.В., Шель Е.В., Яковлев A.A., Пустовских A.A. (2017). Подходы к моделированию гидроразрыва нласта и направления их развития. Нефтяное хозяйство, 12, c. 37-41.

José I. Adachi, Emmanuel Detournay, Anthony P. Peirce. (2010). Analysis of the classical pseudo-3D model for hydraulic fracture with equilibrium height growth across stress barriers. International Journal of Rock Mechanics and Mining Sciences, 47(4), pp. 625-639.

Сведения об авторах

М.М. Хасанов - доктор тех. наук, директор дирекции по технологиям «Газпромнефти», генеральный директор ООО «Газпромнефть ЭТЦ»

Россия, 190000, Санкт-Петербург, наб. реки Mойки, д. 75-79, литер Д

A.Н. Ситников - заместитель генерального директора по научному инжинирингу, ООО «Газпромнефть Hm,» Россия, 190000, Санкт-Петербург, наб. реки Mойки, д. 75-79, литер Д

A.A. Пустовских - канд. физ.-мат. наук., начальник департамента научно-методического сопровождения разработки и профессионального развития ООО «Газпромнефть ЭТЦ»

Россия, 190000, Санкт-Петербург, наб. реки Mойки, д. 75-79, литер Д

E-mail: [email protected]

A.П. Рощектаев - канд. физ.-мат. наук, ведущий эксперт департамента научно-методического сопровождения разработки и профессионального развития ООО «Газпромнефть ЭТЦ»

Россия, 190000, Санкт-Петербург, наб. реки Mойки, д. 75-79, литер Д

Н.С. Исмагилов - канд. физ.-мат. наук, руководитель направления Управления комплексного проектирования разработки, ООО «Газпромнефть ЭТЦ»

Россия, 190000, Санкт-Петербург, наб. реки Mойки, д. 75-79, литер Д

НАУЧНО-ТЕХНИЧЕСКИЙ ЖУРНАЛ

www.geors.ru ГЕОРЕСУРСЫ

Г.В. Падерин - главный специалист Управления комплексного проектирования разработки ООО «Газпромнефть НТЦ»

Россия, 190000, Санкт-Петербург, наб. реки Мойки, д. 75-79, литер Д

Е.В. Шель - ведущий специалист Управления комплексного проектирования разработки ооо «Газпромнефть НТЦ»

Россия, 190000, Санкт-Петербург, наб. реки Мойки, д. 75-79, литер Д

Статья поступила в редакцию 01.08.2018; Принята к публикации 16.08.2018; Опубликована 30.08.2018

Scientific engineering as the basis of modeling processes in field development

M.M. Khasanov, A.N. Sitnikov, A.A. Pustovskikh*, A.P. Roshchektayev, N.S. Ismagilov, G.V. Paderin, E.V. Shel

Gazprom Neft Science and Technology Center, Saint-Petersburg, Russian Federation Corresponding author: Alexey A. Pustovskikh, e-mail: [email protected]

Abstract. Three characteristic examples of the use of scientific engineering approaches for managing the technological processes of reservoir modeling at different hierarchical levels are presented in the article.

The first example demonstrates the application of the spectral approach for modeling geophysical fields - the available log data is decomposed from the spectrum of Legendre polynomials, after which a stochastic field of the expansion coefficients is constructed. The results obtained by this method of realizing geophysical fields correspond to real data in a wider area of modeling than in classical methods. The speed of building models also increases due to the convenience of parallelization.

The second example demonstrates the use of the source method to optimize the transfer of wells into injection. A flow rate of each well is found by simulating the wells in the development system as linear sources or sinks and recording the resulting system of equations for the flows at each time. According to the optimum of discounted extraction, there is an economically efficient time for well development.

The third example demonstrates the application of the theory of dimensions to the problem of hydraulic fracturing modeling to determine the significance of certain parameters for the design of fracturing. By measuring the dependence of the fracture length on the injection volume, we obtain an empirical formula for the fracture length from its parameters, which determines the level of their significance.

Keywords: geological modeling; geostatistics; geomechanics; hydraulic fracturing; hydrodynamic modeling; source method

Recommended citation: Khasanov M.M., Sitnikov A.N., Pustovskikh A.A., Roshchektayev A.P., Ismagilov N.S., Paderin G.V., Shel E.V. (2018). Scientific engineering as the basis of modeling processes in field development. Georesursy = Georesources, 20(3), Part 1, pp. 142-148. DOI: https://doi.org/10.18599/ grs.2018.3.142-148

References

Baykov V. A., Bakirov N.K., Yakovlev A.A. (2010). New approaches to the theory of geostatistical modeling. Vestnik UGATU, 37(2), pp. 209-215. (In Russ.)

Baykov V.A., Bakirov N.K., Yakovlev A.A. (2012). Mathematical Geology. V. 1. Introduction in geostatistics. Moscow-Izhevsk: IKI, Biblioteka nefftyanogo inzhiniringa, 228 p. (In Russ.)

Dyubryul O. (2002). Use of geostatistics to include data in the geological model. EAGE: SEG. 295 p. (In Russ.)

José I. Adachi, Emmanuel Detournay, Anthony P. Peirce. (2010). Analysis of the classical pseudo-3D model for hydraulic fracture with equilibrium height growth across stress barriers. International Journal of Rock Mechanics and Mining Sciences, 47(4), pp. 625-639.

Khasanov M.M., Belozerov B.V., Bochkov A.S., Fuks O.M., Tengelidi D.I. (2015). Automated lithologic facies analysis based on the spectral theory. Neftyanoe khozyaistvo = Oil Industry, 12, pp. 48-51. (In Russ.)

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Khasanov M.M., Belozerov B.V., Bochkov A.S., Ushmaev O.S., Fuks O.M. (2014). Application of the spectral theory to the analysis and modelling

of the rock properties of the reservoir. Neftyanoe khozyaistvo = Oil Industry, 12, pp. 60-64. (In Russ.)

Khasanov M.M., Melchaeva O.Yu., Sitnikov A.N., Roshchektaev A.P. (2013). Dynamics of hydraulically fractured wells production for economically optimal development systems. Neftyanoe khozyaistvo = Oil Industry, 12, pp. 36-39. (In Russ.)

Khasanov M.M., Paderin G.V., Shel E.V., Yakovlev A.A., Pustovskikh A.A. (2017). Approaches to modeling hydraulic fracturing and their development. Neftyanoe khozyaistvo = Oil Industry, 12, pp. 37-41. (In Russ.)

Prigarin S.M., Mikhailov G.A. (2005). Metods of numerical modelling of random processes and fields. Novosibirsk: IVMIMG SO RAN, 259 p. (In Russ.)

Sitnikov A.N., Pustovskikh A.A., Roshchektaev A.P., Andzhukaev T.V. (2015). A method to determine optimal switching time to injection mode for field development system. Neftyanoe khozyaistvo = Oil Industry, 3, pp. 84-87. (In Russian)

About the Authors

M.M. Khasanov - DSc (Engineering), Director General, Gazprom Neft Science and Technology Center