allgosts.ru35. ИНФОРМАЦИОННЫЕ ТЕХНОЛОГИИ. МАШИНЫ КОНТОРСКИЕ35.080. Программное обеспечение

ГОСТ Р 56449-2015 Месторождения газовые, газоконденсатные, нефтегазовые и нефтегазоконденсатные. Программное обеспечение для гидродинамического моделирования месторождений. Основные функциональные и технические требования

Обозначение:
ГОСТ Р 56449-2015
Наименование:
Месторождения газовые, газоконденсатные, нефтегазовые и нефтегазоконденсатные. Программное обеспечение для гидродинамического моделирования месторождений. Основные функциональные и технические требования
Статус:
Действует
Дата введения:
05/01/2016
Дата отмены:
-
Заменен на:
-
Код ОКС:
35.080

Текст ГОСТ Р 56449-2015 Месторождения газовые, газоконденсатные, нефтегазовые и нефтегазоконденсатные. Программное обеспечение для гидродинамического моделирования месторождений. Основные функциональные и технические требования

ФЕДЕРАЛЬНОЕ АГЕНТСТВО ПО ТЕХНИЧЕСКОМУ РЕГУЛИРОВАНИЮ И МЕТРОЛОГИИ

ГОСТР

56449—

2015



НАЦИОНАЛЬНЫЙ

СТАНДАРТ

РОССИЙСКОЙ

ФЕДЕРАЦИИ

МЕСТОРОЖДЕНИЯ ГАЗОВЫЕ, ГАЗОКОНДЕНСАТНЫЕ, НЕФТЕГАЗОВЫЕ И НЕФТЕГАЗОКОНДЕНСАТНЫЕ

Программное обеспечение для гидродинамического моделирования месторождений Основные функциональные и технические

требования

Издание официальное

СпиАфтшфоф*

201*

Предисловие

1 РАЗРАБОТАН Открытым акционерным обществом «Газпром» (ОАО «Газпром») и Обществом с ограниченной ответственностью «Газпром георесурс» (ООО «Газпром георесурс»)

2 ВНЕСЕН Техническим комитетом по стандартизации ТК 23 «Нефтяная и газовая промышленность»

3 УТВЕРЖДЕН И 8ВЕДЕН В ДЕЙСТВИЕ Приказом Федерального агентства по техническому регулированию и метрологии от 15 июня 2015 г. No 670-ст

4 ВВЕДЕН ВПЕРВЫЕ

Правила применения настоящего стандарта установлены в ГОСТ Р 1.0—2012 (раздел 8). Информация об изменениях к настоящему стандарту публикуется в ежегодном (ло состоянию на 1 января текущего года) информационном указателе «Национальные стандарты». а официальный текст изменений и поправок — в ежемесячном информационном указателе «Национальные стандарты». В случае пересмотра (замены) или отмены настоящего стандарта соответствующее уведомление будет опубликовано в ближайшем выпуске ежемесячного информационного указателя «Национальные стандарты». Соответствующая информация, уведомление и тексты размещаются также е информационной системе общего пользования — на официальном сайте Федерального агентства по техническому регулированию и метрологии в сети Интернет ()

© Стамдартинформ. 2016

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

II

Содержание

7 Функциональные требования к программному обеспечению гидродинамического

8 Требования к документированию программного обеспечения гидродинамического

9 Технические требования к программному обеспечению гидродинамического моделирования .... 10

10 Подтверждение соответствия программного обеспечения для гидродинамического

Приложение Б (справочное) Задание условий на границах расчетной области в рамках

Приложение Г (справочное) Матрица применения технологий программного обеспечения

НАЦИОНАЛЬНЫЙ СТАНДАРТ РОССИЙСКОЙ ФЕДЕРАЦИИ

МЕСТОРОЖДЕНИЯ ГАЗОВЫЕ. ГАЗОКОНДЕНСАТНЫЕ. НЕФТЕГАЗОВЫЕ И НЕФТЕГАЗОКОНДЕНСАТНЫЕ

Программное обеспечение для гидродинамического моделирования месторождений Основные функциональные и технические требования

Gas, gas condensate, oil. gas and condensate deposits. Software for flow simulation o( hydrocarbon fields.

Main functional and technical requirements

Дата введения — 2016—05—01

1 Область применения

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

1.2 Положениями настоящего стандарта руководствуются субъекты хозяйственной деятельности:

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

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

2 Нормативные ссылки

8 настоящем стандарте использованы нормативные ссылки на следующие стандарты:

ГОСТ 8.417—2002 Государственная система обеспечения единства измерений. Единицы величин

ГОСТ 19.105—78 Единая система программной документации. Общие требования к программным документам

ГОСТ 19.601—78 Единая система программной документации. Общие правила дублирования, учета и хранения

ГОСТ 19.603—78 Единая система программной документации. Общие правила внесения изменений

ГОСТ 2939—63 Газы. Условия для определения объема

ГОСТ 28195—89 Оценка качества программных средств. Общие положения

ГОСТ Р 8.645—2008 Государственная система обеспечения единства измерений. Метрологическое обеспечение работ по геологическому изучению, использованию и охране недр в Российской Федерации. Основные положения

РОСТР 53712—2009 Месторождения нефтяные и газонефтяные. Программные средства для проектирования и оптимизации процесса разработки месторождений. Основные требования

ГОСТ Р ИСО/МЭК 9126—93 Информационная технология. Оценка программной продукции. Характеристики качества и руководства по их применению

ГОСТ Р ИСО/МЭК 12119—2000 Информационная технология. Пакеты программ. Требования к качеству и тестирование

ГОСТ Р ИСО/МЭК 15910—2002 Информационная технология. Процесс создания документации пользователя программного средства

Издание официальное

Примечание — При пользовании настоящим стандартом целесообразно проверить действие ссылочных стандартов в информационной системе общего пользования — на официальном сайте Федерагъного агентства по техническому регулированию и метрологии в сети Интернет или по ежегодному информационному указателю «Национальные стандарты», который опубликован по состоянию на 1 января текущего года, и по выпускам ежемесячного информационного указателя «Национальные стандарты» за текущий год. Если заменен ссылочный стандарт, на который дана недатированная ссылка, то рекомендуется использовать действующую версию этого стандарта с учетом всех внесенных е данную версию изменений. Если заменен ссылочный стандарт, на который дана датированная ссылка, то рекомендуется использовать версию этого стандарта с указанным выше годом утверждения (принятия). Если после утверждения настоящего стандарта в ссылочный стандарт, на который дана датированная ссылка, внесено изменение, затрагивающее положение, на которое дана ссылка, то это положение рекомендуется применять без учета данного изменения. Если ссылочный стандарт отменен без замены, то положение, в котором дана ссылка на него, рекомендуется применять в части, не затрагивающей эту ссылку.

3 Термины и определения

В настоящем стандарте применены следующие термины с соответствующими определениями:

3.1 адаптация параметра: Процесс минимизации отклонения расчетного значения параметра от его фактического значения.

3.2

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

[ГОСТ Р 53554—2009. статья 2.69)

3.3 водоносный горизонт: Осадочная горная порода, представленная одним или несколькими переслаивающимися подземными слоями горных пород с различной степенью водопроницаемости, пустоты {поры, трещины) которых заполнены гравитационными водами.

3.4

газовый фактор: Отношение полученного количества нефтяного газа, приведенного к атмосферному давлению и температуре 20 *С. к количеству полученной за то же время нефти при том же давлении и той же температуре.

[ГОСТ Р 53554—2009. статья 2.94)

3.5

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

[ГОСТ Р 53712—2009. пункт 3.12]

Примечание — В данном стандарте для разделения программного обеспечения (ПО) сборных сетей и ПО гидродинамического моделирования фильтрационная (гидродинамическая) модель рассмотрена до забоя скважин.

3.6

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

[ГОСТ Р 53554—2009. статья 2.103]

3.7 давление насыщения: Давление, при достижении которого в процессе снижения давления начинается выделение газа из нефти.

3.8

дебит эксплуатационной скважины: Количество нефти и (или) газа, полученное из эксплуатационной скважины е единицу времени.

(ГОСТ Р 53554—2009. статья 2.113]

3.9 дифференциальная конденсация: Процесс образования жидкой фазы в газоконденсатной (многокомпонентной) смеси при ступенчатом выпуске из сосуда (бомбы PVT) паровой фазы (ступенчатом изменении давления) и неизменной температуре.

3.10

забойное давление: Давление флюида на забое буровой скважины при добыче нефти и газа.

(ГОСТ Р 53554—2009. статья 2.42]

Примечание — Давление на забое простаивающей (или временно остановленной) скважины называется пластовым давлением.

3.11 капиллярное давление: Разность значений давления на поверхности раздела фаз.

3.12 композиционная модель: Модель, в которой углеводородная пластовая система представлена в виде заданного количества компонент, каждый из которых характеризуется своими параметрами.

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

3.14 локальные измельчения: Разукрупнение геологической сетки с целью детализации выбранной области.

3.15 массив (проницаемости/песчанистости/пористости/связной водонасыщенности/оста-точной нефтенасыщенкости/остаточной газонасыщенкости) данных о характере насыщения:

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

3.16

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

(ГОСТ Р 53554—2009. статья 2.38]

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

3.18 модель нелетучей нефти: Модель, в которой углеводородная система аппроксимирована двумя компонентами: нелетучим (нефтью) и летучим (газом).

3.19 неоднородность объекта: Изменение фильтрационно-емкостных свойств по площади и разрезу объекта.

3.20 непротиворечивость гидродинамической модели: Соответствие гидродинамической модели физическим представлениям о геологическом объекте и протекающих в нем фильтрационных процессах.

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

3.22

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

[ГОСТ Р 53554—2009. статья 2.41]

3.23

программное обеспечение: Совокупность программ системы обработки информации и про* граммных документов, необходимых для эксплуатации этих программ.

[ГОСТ 19781—90. статья 5.2]

3.24 сланцевый газ: Природный газ. добываемый из горючих сланцев и состоящий преимуще* ственно из метана.

3.25

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

[ГОСТ Р 53712—2009, пункт 3.11]

3.26

уравнение состояния природного газа: Уравнение, которое связывает фактор сжимаемости, температуру, плотность и молярные доли компонентов.

[ГОСТ 30319.0—96. пункт 3.5]

3.27

устьевое давление: Давление, измеренное в верхней точке буровой скважины.

[ГОСТ Р 53554—2009. статья 2.45}

3.28 фазовые превращения: Переход вещества из одной термодинамической фазы в другую при изменении внешних условий.

3.29 фазовое равновесие: Одновременное существование термодинамически равновесных фаз в многофазной системе.

4 Сокращения

В настоящем стандарте применены следующие сокращения:

внк

водонефтяной контакт:

Г8К

газоводяной контакт:

там

гидродинамическое моделирование;

таи

гидродинамические исследования:

ГНК

газонефтяной контакт:

ГРП

гидравлический разрыв пласта:

ГС

горизонтальный ствол;

ГТМ

геолого-технические мероприятия:

дкс

дожимная компрессорная станция;

КГФ

конденсатно*гаэовый фактор;

КРС

капитальный ремонт скважин;

нкт

насосно*компрессорные трубы;

ОФП

относительные фазовые проницаемости;

ПАВ

поверхностно-активные вещества;

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

пд

пзп

по

пхг

У8

PVT


pressure volume temperature (давление объема температуры).

5 Общие положения

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

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

6 Исходные данные для программного обеспечения гидродинамического

моделирования

6.1 Исходными данными для ПО ГДМ являются:

6.1.1 Статическая основа:

• геометрия гидродинамической сетки:

• массивы проницаемости (диагональные компоненты тензора проницаемости);

• массив песчанистости;

• массив пористости;

• массив связанной водонасыщенности;

• массив остаточной насыщенности углеводородными фазами (нефть, газ);

• массив начальных насыщенностей;

• массив гидродинамически изолированных сегментов модели (геометрические координаты нарушений пласта):

• массив порового объема;

• массив проводимостей;

• массив регионов подсчета запасов;

• массив регионов сжимаемости породы;

• массив регионов классов коллекторов;

• массив PVT-свойств флюида:

• отметки ГНК. ВНК или Г8К по месторождению в цепом или по гидродинамически изолированным сегментам модели;

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

• газовый и/или конденсатно-газовый фактор (или начальный компонентный состав) по месторождению в целом, или по гидродинамически изолированным сегментам модели, или в зависимости от абсолютной отметки;

• результаты оценки начальных запасов и подсчетных параметров.

6.1.2. Сжимаемость породы (в том числе в виде таблиц, описывающих нелинейные зависимости пористости и проницаемости от пластового давления).

6.1.3 Функции относительных фазовых проницаемостей.

6.1.4 Зависимости давления в системе нефть—газ и вода—нефть, газ—углеводородный конденсат. углеводородный конденсат—вода (газ—вода при отсутствии жидкой углеводородной фазы) от насыщенности или модель J-функции.

Примечание — ./-функция является по своей сути нормированной величиной капиллярного давления. В качестве критериев нормирования выступают силы поверхностного натяжения между флюидами и фильтрационно-емкостные свойства матрицы:


(D

где P,(S„) — капиллярное давление (функция S0), Па;

Se — еодонасыщенмость. д.ед.; о — натяжение на границе раздела фаз. дин/см; о — угол смачивания (краевой угол смачивания), градус; к — проницаемость, определенная по образцу керна, мД; т — пористость, определенная по образцу керна, д.ед.

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

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

Примечание — При больших потоках флюида зависимость градиента давления от скорости течения переходит из линейной, описанной законом Дарси, е квадратичную модель Форхгеймера.

(2)

где ц — вязкость. Па-с;

К — проницаемость. мД; р — плотность, кг/м5; v — средняя скорость потока, м/с; в уравнении — обратная длина (1/м).

Модель Форхгеймера необходима в зонах с большими скоростями потоков, например около скважин.

6.1.7 PVT-свойства нефти и газа в виде:

- компонентного состава и свойств компонентов (критические параметры, коэффициенты бинарного взаимодействия и т. д.);

* таблиц, описывающих изменение плотности, вязкости, и объемных коэффициентов давления для разных значений газового фактора или КГФ. а также зависимость давления насыщения или начала конденсации от значения газового фактора или КГФ.

6.1.8 Свойства пластовой и закачиваемой воды.

6.1.9 Данные по скважинам:

- траектории скважин;

- интервалы перфорации скважин в измеренных или абсолютных глубинах;

- даты проведения перфорации, изоляции, дополнительной перфорации.

6.1.10 Данные о конструкциях скважин: диаметр скважины, диаметр и глубина спуска НКТ. устройства контроля притока и т. д.

6.1.11 Данные о проведенных ГРП. КРС и других ГТМ.

6.1.12 История разработки:

- дебиты газа, нефти, конденсата и воды по добывающим скважинам и/или группам добывающих скважин;

- расходы газа и воды по нагнетательным скважинам и/или группам нагнетательных скважин;

• забойные и устьевые давления по скважинам.

6.2 Единицы измерения исходных данных

Исходные данные для ПО ГДМ измеряют е системе СИ в соответствии с ГОСТ 8.417 и ГОСТ Р 8.645.

7 Функциональные требования к программному обеспечению

гидродинамического моделирования

7.1 Импорт, подготовка и проверка исходных данных, перечисленных в разделе 6.

7.2 Создание численной модели пластового флюида

7.2.1 Моделирование следующих лабораторных экспериментов:

• однократная сепарация;

• ступенчатая сепарация:

. дифференциальное разгазирование:

• контактная конденсация (исследование при постоянной массе);

• дифференциальная конденсация;

• определение давления насыщения;

• однократный расчет фазового равновесия.

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

7.2.3 Настройка параметров уравнения состояния по результатам лабораторных PVT-экспериментов. в том числе с использованием процедуры многомерной регрессии.

7.2.4 Возможность использования различных уравнений состояния.

Прим еча нив — Как правило, используют уравнения состояния Пенга-Робинсона. Соаее-Редлиха-Каонга. Зюдкевича-Йоффе. описание приведено в приложении А.

7.3 Задание начального состояния модели заключается:

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

• прямом задании неравновесного начального состояния модели.

7.4 Задание условий на границах расчетной области, в том числе в рамках моделей водоносных горизонтов Фетковича и Картера-Трейси.

Прим еча н и в — Методы определения водоносного пласта Фетковича и Картера-Трейси приведены в приложении Б.

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

7.6 Оценка начальных и текущих (остаточных) запасов УВ.

7.7 Контроль непротиворечивости модели.

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

7.9 Выполнение прогнозных расчетов для различных вариантов эксплуатации месторождения, учитывающих:

• время, в течение которого осуществляют прогнозирование;

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

• предполагаемый график перевода скважин на другие объекты и зарезки боковых стволов;

• планируемую модернизацию систем сбора, подготовки и нагнетания:

• технологические ограничения на добычу нефти, воды и газа;

• технологические ограничения на нагнетание воды и газа;

- ограничения на промысловый газовый фактор (для отдельных скважин, групп скважин и месторождения);

• ограничения на обводненность продукции (для отдельных скважин, групп скважин и месторождения);

- минимальный экономический дебит нефти, конденсата или газа;

- максимальный дебит нефти или газа по скважинам (если применим);

• ограничения на суммарный уровень добычи нефти газа или воды;

• ограничения по минимальному и максимальному давлению на устье скважин.

7.10 Параллельные вычисления — это:

• ускорение гидродинамических расчетов за счет использования параллельных процессоров;

• возможности запуска с использованием системы очередей (создание и управление очередями) удаленного запуска на расчет и мониторинга расчета на удаленном кластере.

7.11 Учет притока к газовым скважинам подразумевает:

• учет реальных свойств газа при расчете притока (нелинейность зависимости дебита от депрессии);

• наличие модели эффекта блокировки газоконденсатной скважины конденсатом при снижении забойного давления ниже давления начала конденсации:

• возможность учета отклонения от линейного закона при фильтрации газа в пласте;

• возможность вести упрощенный учет отклонения от линейного закона фильтрации в призабойной зоне скважины (скик-фактор. зависящий от дебита. — так называемый О-фактор).

7.12 Возможность детализации сетки в околоскважинной зоне с помощью:

• прямоугольных локальных измельчений сетки;

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

7.13 Моделирование закачки полимеров с возможностью учета:

- изменения вязкости закачиваемой воды с увеличением концентрации полимера:

• адсорбции полимера со снижением проницаемости;

• десорбции полимера.

7.14 Моделирование закачки С02 представляет собой:

• учет растворимости закачиваемого С02 в воде и нефти;

• моделирование закачки углекислого газа в нефтяные пласты для увеличения нефтеотдачи.

7.15 Моделирование водогазового воздействия — моделирование поочередной закачки воды и газа, а также совместной одновременной закачки.

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

7.17 Учет течения флюида в стволах ГС дает возможность:

• учета трения в стволе скважины при моделировании горизонтальных скважин;

• использования модели многосегментной скважины для точного описания поведения многофазного потока в стволе скважины, учета геометрии скважины в процессе расчета;

- моделирования многозабойных (раэветвленно-гориэонтальных) скважин в рамках модели многосегментной скважины.

7.18 Моделирование ГРП представляет собой моделирование:

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

• многостадийного ГРП. оценка результатов.

7.19 Композиционное моделирование — это возможность:

• детально учитывать состав пластовых флюидов и фазовые превращения на основе композиционной модели;

• моделирования процессов закачки газа, смешивающегося вытеснения;

• задания стандартных условий для расчета объемов УВ в соответствии с российскими стандартами (0.101325 МПа. 293.15 К).

7.20 Псевдокомпозиционное моделирование — это возможность упрощенно учитывать композиционные эффекты при закачке газа, отличного от газа сепарации в рамках модели нелетучей нефти.

7.21 Модели двойной пористости/двойной проницаемости.

7.22 Моделирование добычи метана из угольных пластов подразумевает:

• учет эффектов адсорбции/десорбции в зависимости от давления;

- возможность моделирования закачки C02/N2;

• возможность моделирования изменения пористости и проницаемости с изменением давления.

7.23 Моделирование добычи сланцевых газов — учет:

• эффектов адсорбции/десорбции в зависимости от давления;

• образования сети проницаемых трещин при проведении ГРП.

7.24 Визуализация и анализ результатов моделирования включает в себя:

• возможность работы с разными сетками (например, блочно-центрированная, геометрия угловой точки);

- возможность работы с радиальными и декартовыми локальными измельчениями;

• возможность визуализации свойств модели на любой временной шаг;

• возможность манипуляции различными секциями модели;

- возможность работы с геологическими разрезами модели;

• возможность установления диапазона различных свойств модели при визуализации;

• возможность построения графиков технологических показателей разработки, полученных по результатам расчетов по скважинам, группам скважин, объектам, месторождениям:

• возможность цветокодирования. при которой используют ту или иную цветовую палитру.

• построение 20-карт по результатам моделирования и возможность экспорта в современные графические редакторы;

• построение кроссплотов:

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

7.25 Возможность назначения свойств закачиваемой воды, отличных от свойств воды в резервуаре [2).

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

7.27 Возможность модификации исходных массивов, перечисленных в 6.1.1.

7.28 Возможность создания новых скважин.

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

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

7.31 Онлайн-интерфейс для совместного использования с системами:

• моделирования течения в скважинах и поверхностных сетях сбора и транспорта (в том числе для моделирования газосборной сети и ДКС на каждой расчетной итерации, при этом наземное оборудование должно являться ограничителем добычи);

• геомеханического моделирования.

7.32 ПО ГДМ должно обеспечивать возможность его использования под управлением систем автоматизированной адаптации и расчета неопределенности показателей разработки месторождения, а также систем интерпретации ГДИ скважин.

8 Требования к документированию программного обеспечения

гидродинамического моделирования

8.1 Для ПО ГДМ разрабатывают ПД в соответствии с требованиями ГОСТ Р ИСО/МЭК 15910 и ГОСТ Р 53712. Общие требования к ПД соответствуют ГОСТ 19.105.

8.2 Дублирование, учет и хранение ПД проводят в соответствии с требованиями ГОСТ 19.601. Изменения в ПД вносят в соответствии с ГОСТ 19.603.

8.3 ПД на программное средство и (или) его компоненты должна предусматривать наличие подробного руководства пользователя на русском языке и содержать следующие сведения:

• область применения;

• данные о структуре (перечень модулей, базовый комплект, список дополнительных модулей);

• описание пользовательского интерфейса:

• описание опций;

- описание ключевых слов;

• описание порядка экспорта и импорта данных;

• описание применяемых алгоритмов и физических моделей.

9 Технические требования к программному обеспечению

гидродинамического моделирования

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

- функционирование на современных операционных системах и аппаратных платформах;

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

• обеспечение возможности конвертации моделей, созданных в более ранних версиях ПО ГДМ. в формат обновленной версии ПО.

10 Подтверждение соответствия программного обеспечения

для гидродинамического моделирования

Подтверждение соответствия ПО ГДМ оценивают в соответствии с ГОСТ 28195. ГОСТ Р ИСО/ МЭК 9126 и ГОСТ Р ИСО/МЭК 12119 следующими способами:

• экспертный анализ руководства пользователя с проверкой наличия описания функциональности. реализующей положения настоящего стандарта;

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

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

Примечание — Обычно компании — разработчики ПО ГДМ используют испытательные тесты международной некоммерческой профессиональной организации технических специалистов нефтегазовой промышленности Society of Petroleum Engineers (тесты SPE).

Уравнения состояния

Как правило, используют уравнения состояния Пвнга-Робинсона. Соаее-Редлиха-Кеонгэ. Зюдкевича-Иоффе. В общей форме кубическое уравнение состояния может быть представлено таким образом:

(А.1)


ЯГ л

где р — абсолютное давление. Па:

R— универсальная газовая постоянная. Дж/(моль-К):

Г — абсолютная температура. К: v — молярный объем, мкмоль, а и б — параметры уравнения оосгояния: с, и с2 — константы, принимающие следующие величины 8 уравнениях состояния:

• Редлиха-Квонга 1 и 0.

• Соаве-Редлиха-Квонга 1 и 0.

• Пента-Робинсона 1 + \^2 и 1 - ^2 соответственно.

Параметры а и б уравнения состояния для компонентов выражены следующим образом:

(A.2)

(A.3)

(A-4)

(A.5)

оде функция 0(£а)едина для уравнений состояния; и» — ацентрический фактор: си г — критические и приведенные условия, a. Q,, и Ц, зависят от выбора уравнения состояния: Редлиха-Квонга: где А* = 0.42748;

Qj = 0.08664:

(А.6)

Соаве-Рвдлиха-Квонга:

(А.7)

(А.8)


где

* = 0.48+ 1.57ш-0,17о>2;

= 0.42748:

Oj= 0.08664.

Пенга-Роби неона: где Gj = 0.45724; flj= 0,07780;

(А.9)

к = 0.37464+ 1.5422о. - 0,26992<.>*. (АЛО)

Для более тяжелых углеводородных компонентов применяют другое значение ik вуравненин Пенга-Робинсона:

к = 0.379642 + 1.4803.U - 0.164423ч>2 + 0.016666«А (А.11)

Уравнение состояния Зодкевича-Иоффе — это да ух параметрическое уравнение состояния, которое было получено этими учеными с целью улучшения уравнения состояния Редлиха-Квонга:

- константы Qa, Оь зависят от температуры, т. е. Qe (7), Q6(T):

- fla (7). Оъ(Т) получают для каждого компонента, основываясь на плотности насыщенной жидкости и уравнения для летучести:

- Qe, Оъ берутся как не зависящие от температуры значения при температуре выше критической.

Задание условий на границах расчетной области в рамках моделей Фетковича

и Картера-Трейси

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

(Б.1)

«-


среднее давление в водоносном пласте в момент времени Па; коэффициент продуктивности для водоносного пласта: площадь ячейки /, м?

После завершения шага по времени рассчитывают полный приток воды из водоносного пласта следующим уравнением;

(Б.2)


то

Затем, используя уравнение материального баланса, пересчитывают новое значение для среднего давления

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

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

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

Модель Картера-Трейси — подход аналитического определения водоносного пласга, основанный на применении «функций влияния» для расчета полного притока воды из водоносного пласта и расчета изменения давления с течением времени между водоносным пластом и средним давлением на границе водоносного пласта и залежи. Функции влияния зависят от геометрии и свойств водоносного пласта. Два основных параметра определяют поведение водоносного пласта:

Константа времени (имеющая размерность времени)

(Б.З)

где — вязкость воды в водоносном пласте. Па-с:

Ф— пористость водоносного пласта. Д.ВД.:

С( — полная (порода + вода) сжимаемость. 1/Па;

Гф — внешний радиус пласта (или внутренний радиус водоносного пласта), м; кв — проницаемость водоносного пласта. мД:

с, = 0.008527 (константа для перевода величин из одной системы единиц в другую, значение приведено для Международной системы единиц СИ).

Консганта притока водоносного пласта (имеющая размерность полного притока на единицу падения давления)

ЬяС^МРС,#. (Б.4)

где h—толщина водоносного пласта, м;

0 — угол, под которым граница водоносного пласга видна из центра залежи, в градусах, деленный на 360: С2 = 6.283 (константа для перевода величин из одной системы единиц 8 другую, значение приведено для Международной системы единиц СИ).

Константу времени используют для преобразования времени в безразмерную форму с помощью соотношения


(Б.5)

Модель Картера-Трейси выражает падение давления на границе водоносного пласта через функцию влияния безразмерного давления следующим образом:

где Р*-Р —

оа-


начальное давление воды в водоносном пласте. Па:

среднее давление воды на границе водоносного ппасга/залежи, Па:

скорость притока водоносного пласта, м3/с.

Тесты корректности реализации специальных уравнений притока к скважине, учитывающих изменение PVT-свойств флюидов в призабойной зоне, включая возможные фазовые превращения

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

• уравнение Russell Goodrich позволяет учесть эффект призабойной зоны для идеального газа постоянной вязкости;

• псевдодавлемие газа позволяет учесть эффект призабойной зоны для реального газа (без учета фазовых превращений);

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

Однако в системе тестов SPE отсутствуют тесты на проверку корректности реализации этих уравнений.

В.1 Тест 1 Случай идеального газа постоянной вязкости

Полный перечень параметров модели приведен в таблице В.1.

В.1.1 Описание модели пласта

Рассматривают горизонтальный однородный изотропный пласт круглой формы (в плане). Размерность сетки фильтрационной модели — 21x21x1. Радиус круга (по центрам ячеек) — 1000 м. Размер ячейки в плане — 100 м. Мощность пласта — 10 м. Пористость — 20 %. проницаемость — 0,01 мкм2.

В.1.2 Свойства флюидов

Однофазная модель содержит только газ. Идеальный газ с постоянной вязкостью — 0.01 мПа-с.

В.1.3 Начальные и граничные условия

В граничных ячейках либо за счет увеличения норового объема, либо иным способом задают условие постоянного давления (контур питания). Начальное пластовое давление — 15 МПа. Вертикальная добывающая скважина расположена 8 центре модели е ячейке (11,11,1) и вскрывает перфорацией всю толщину пласта. Забойное давление — 10 МПа. Расчет выполнен на один год. Дебит газа на конец расчета считают установившимся.

Рисунок В.1 — Начальные и граничные условия

Таблица В.1 — Параметры модели теста 1

Параметр

Значение

Размерность модели

21x21x1

Размер ячейки сетки DX*DY*DZ, м

100x100x10

Мощность пласта h, м

10

Пористость. д.ед.

0.2

Проницаемость к. мкм2

0.01

Температура. "С

20

Вязкость газа р, мПа-с

0.01

Объемный коэффициент газа. д.ед.

1/р. где р — давление е 0.1 МПа

Радиус контура питания Rc, м

1000

Радиус скважины м

0.1

Давление на контуре питания рс, МПа

15

Параметр

Значение

Окончание таблицы В. 1

Параметр

Значение

Забойное давление pw МПа

10

Сжимаемость пор породы. 1/МПа

з-ю-4

Время расчета для получения дебита газа на установившемся режиме

В. 1.4 Эталонные результаты моделирования

Объемный дебит газа в стандартных условиях при установившемся плосхорадиагъном притоке (температура е пласте и на поверхности совпадает) вычисляют по формуле

<В.1)

где —давление в стандартных условиях.

Учитывая формулу объемного коэффициента газа, р^ = 0.1 МПа (не 0.101325 МПа. как обычно).

Тогда эталонное значение объемного дебита газа рвано: О = 368.382 тыс. м3/сут

В.1.5 Критерии оценки

Оценивают отклонение расчетного значения дебита газа через один год после начала моделирования от эталонного значения 368.382 тыс. м3/сут.

Допустимая относительная погрешность — 5 %.

В.2 Тест 2 Случай реального газа

Полный перечень параметров модели приведен в таблице В.2.

8.2.1 Описание модели пласта

Рассматривают горизонтальный однородный изотропный пласт круглой формы (в плане). Размерность сетки фильтрационной модели — 21x21x1. Радиус круга (по центрам ячеек) — 1000 м. Размер ячейки а плане — 100 м. Мощность пласта — 10 м. Пористость — 20 %. проницаемость — 0.01 мкм2.

В.2.2 Свойства флюидов

Модель однофазная и содержит только газ.

Сверхсжимаемосгь Z задана зависимостью: Z(p) = {р/рЛ) 0.5.

Вязкость газа. мПа-с, задана зависимостью: ц(р) = 0,01 p/pw-

В.2.3 Начальные и граничные условия

В граничных ячейках либо за счет увеличения пороеого объема, либо иным способом задают условие постоянного давления (контур питания). Начальное пластовое давление — 15 МПа. Вертикальная добывающая скважина расположена а центре модели в ячейке (11,11,1) и вскрывает перфорацией всю толщину пласта. Забойное давление — 10 МПа. Расчет выполнен на один год. Дебит газа на конец расчета считают установившимся.

Рисунок В.2 — Начальные и граничные условия Таблица В.2 — Параметры модели теста 2

Параметр

Значение

Размерность модели

21x21x1

Размер ячейки сетки DX*DY*DZ. м

100x100x10

Мощность пласта Л. м

10

Пористость, д.ед.

0,2

Проницаемость к. мкм2

0,01

Температура, *С

20

Вязкость газа р, мПа-с

H(p)=0.01-p/pw

Коэффициент сверхсжимаемосги газа. д.ед.

Z{p) = {ptpw?*

Радиус контура питания R^ м

1000

Радиус скважины м

0.1

Давление на контуре питания рс, МПа

15

Забойное давление МПа

10

Сжимаемость пор породы. 1/МПа

3-10-4

Время расчета для получения дебита газа на установившемся режиме

В.2.4 Эталонные результаты моделирования Функцией Лейбом зона называется следующая функция:

*ФШ

ш


ф


(В.2)


где А(р} — проницаемость; ц(р) — вязкость флюида; р(р) — плотность флюида.

Функция Лейбенэона определена с точностью до константы (как потенциал поля скоростей потока).

Для плоскорадиального потока (с нулевым скин-фактором) справедлива следующая формула массового дебита (жидкости или газа):

(В.З)


**Ъ

аде Qm h

rw

Pc

Риг


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

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

«=-у* } (В.4)

•w

Таким образом, значение дебита газа, рассчитанное аналитически, равно: Q = 264.934 тыс. м3/сут.

В.2.5 Критерии оценки

Оценивают отклонение расчетного значения дебита газа через один год после начала моделирования от эталонного значения 264.934 тыс. м3/сут.

Допустимая относительная погрешность — 5 %.

В.З Тест 3 Нефтяная залежь, разрабатываемая на забойных давлениях ниже давления насыщения

Полный перечень параметров модели приведен в таблице В.З.

В.3.1 Описание модели пласта

Рассматривают горизонтальный однородный изотропный пласт круглой формы (в плане). Размерность сетки фильтрационной модели — 21x21x1. Радиус круга (по центрам ячеек) — 1000 м. Размер ячейки а плане — 100 м. Мощность пласта — 10 м. Пористость — 20 %, проницаемость — 0,1 мкм2.

В.3.2 Свойства флюидов

Модель двухфазная: нефть, газ.

Свойства флюидов и зависимости PVT-свойств от давления представлены в таблице В.З.

8.3.3 Относительные фазовые проницаемости

Остаточная нвфтвнзсыщенность (при вытеснении нефти газом) равна 10 %. ОФП заданы квадратичными зависимостями:

(1-^ т» sw

a m*>w.

где s — гаэокасыщекность.


m^*,h***w,

6*-4, Гр** >0,9,


(В.6)


ЯИРвтт


Рисунок В.З — О гноем тельные фазовые проницаемости В.3.4 Начальные и граничные условия

В граничных ячейках либо за счет увеличения норового объема, либо иным способом задают условие постоянного давления {конгур питания). Начальное пластовое давление —15 МПа. Вертикальная добывающая скважина расположена в центре модеш в ячейке (11.11.1) и вскрывает перфорацией всю толщину пласта. Забойное давление — 3 МПа. Расчет выполнен на один год. Дебит нефти на конец расчета считают установившимся.

Рисунок В.4 — Начальные и граничные условия

Таблица В.З — Парэмегры модели теста 3. 8 зависимостях свойств от давления, единица измерения давления бар (0.1 МПа)

Параметр

Значение

Размерность модели

21x21x1

Размер ячейки сетки DX*DY*DZ. м

100x100x10

Мощность пласта h. м

10

Пористость. д.ед.

0.2

Проницаемость к, мкм2

0.1

Радиус контура литания Яе, м

1000

Радиус скважины и

0.108

Давление на контуре питания МПа

15

Давленые насыщения нефти газом МПа

9.5

Забойное давление р„. МПа

3

Объемный коэффициент нефти предельно насыщенной газом SQ. д.ед

0.95-рО.ОЭЗ

Сжимаемость нефти. 1/МПа

8-10-4

Вязкость нефти предельно насыщенной газом мПа с

7000'рг*-5

Сжимаемость вязкости нефти. 1/МПа

0.01

Окончание таблицы В.З

Параметр

Значение

Гаэосодвржание Rs нефти, предельно насыщенной газом. м33

0.056-р1-5

Объемный коэффициент газа Вд. дед.

р-1Д2

Вязкость газа ц^. мПа-с

8.8-10-4- рО-1168

Плотность разгазироваиной нефти, кг/м3

863

Молярная масса нефти, г/моль

212

Молярная масса газа, г/моль

31.3


В.3.5 Эталонные результаты моделирования

Обобщенное псеедодааление позволяет свести задачу установившегося движения многокомпонентной многофазной смеси к задаче движения (однофазной) жидкости.

Основа этого метода разработана С.А. Христианоаичем в конце 40-х годов XX века. Позднее Whitson и Fevang предложили использовать обобщенное псевдодавление для моделирования эффекта образования в призабойной зоне газоконденсатной банки.

Рассматривают установившееся плоскорадиальное течение многокомпонентной многофазной смеси в однородном изотропном пласте.

Многофазное течение состоит из лр фаз. которые содержат пс компонентов. Не обязагегъно каждая фаза должна содержать асе компоненты, но в общем случае каждый компонент может в определенных пропорциях присутствовать во всех фазах. Используют индекс г для обозначения компонентов, а индексу — для обозначения фаз.

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

Через X0 обозначают молярную долю компонента / в фазе /.

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

Из обобщенного закона Дарси следует, что молярная скорость фильтрации пропорциональна градиенту давления

(В.7)

где Vf— молярная плотность потока комитента г в фазе/. моль/м2с;

S = (S,. S2.....Snp) — вектор насыщенности фаз.

Все его компоненты неотрицагегъные. а их сумма равна 1;

А4у (p.S) — обобщенная молярная подвижность компонента; в фазе /

где bfp) —

т-Vs) -

Pj(P) —


молярная плотность фазы у: абсолютная проницаемость: относительная фазовая проницаемость фазы/, вязкость фазы /

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

обобщенная молярная подвижность компонента F,

Аналогично:


через V.


Обозначают:


Г'

л.


суммарную молярную плотность потока компонента t.


— полную молярную плотность потока.


дfm 2jUf— полная обобщенная молярная подвижность.

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

.Д.Ш


ООП*.


(В.9)


Это означает, что вектор насыщенности фаз S не может быть произвольным, а имеет при каждом давлении строго определенное значение. Таким образом, если известна доля каждого компонента в потоке, то можно выразить насыщенность фаз S как функцию давления S(p). Часто это трудно сделать аналитически, но. по крайней мере, внутри симуляторов возможно численное решение этой задачи. В настоящем стандарте не будет рассмотрен вопрос существования и единственности для каждого давления р такого вектора S(p). что выполнено (В.9). Обычно ОФП и свойства флюидов заданы таким образом, что это условие выполнено.

Таким образом, обобщенные подвижности становятся известными функциями давления р. Тогда из (В.7)дпя каждого компонента г

-10)

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

(В-11)

Выражая из(В.11) и подставляя в (В. 10). получают:

<В12>

Интегрируя левую и правую части данного уравнения, получают формулу для молярного дебита компонента /

*-* w Ш


(В.13)


'чф

Г¥

Функцию Щр) в российской литературе называют функцией Хрисгиановича.

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

При этом дебит каждого компонента можно найти по его известной доле в потоке

(В.15)


Q, = 2,'0.

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

<В.16>


а-


ш

•я


[wW-«W!=^C«W*


(В.14)


_

В формуле (В.16) подчеркнута возможность зависимости молярной плотности компонента от температуры. При выводе формулы (В.14) предполагалось, что температура в пласте постоянна, но она может не совладать с температурой в стандартных условиях. Если в стандартных условиях компонент входит более чем в одну фазу, то формула (В.16) теряет смысл, либо молярная плотность Ьсот(и должна быть рассчитана с учетом долей компонента во всех фазах.

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

■**


№+**)


(В.17)


Вектор насыщенностей S = (s0> sff) выражен как функция давления. Так как sg « 1 - s^. то достаточно это сделать только для газонасыщенносги Очевидно, что при давлении выше давления насыщения рь гаэокасы-щенность равна нулю. При давлении ниже давления насыщения рь в плате будут присутствовать как фаза нефти, содержащая нефть и растворенный газ, так и фаза газа, содержащего только компонент газа. Учитывая принятые зависимости (см. таблщу В.З) и ОФП, молярные подвижности при этом имеют вид

и„-тгтш-4?

(В.18)

(В.19)

*,-*«т*

(В.20)

(В.21)

гав

(В.22)

чищ.

^ шш

(В.23)

ц и>№)

(В.24)

(В.28)


где


Отсюда


(В.26)


(В-27)


(В-25)


Тогда


* %лю-ц?

* и


После несложных алгебраических преобразований получают


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

В.3.6 Критерии оценки

Оценивают отклонение расчетного значения дебита нефти через один год после начала моделирования от эталонного значения 62.455 м3/сут.

Допустимая относительная погрешность — 5 %.

В.4 Тест 4 Газоконденсатная залежь, разрабатываемая на забойных давлениях ниже давления

начала конденсации

Полный перечень параметров модегы приведен в таблице В.4.

В.4.1 Описание модели пласта

Рассматривают горизонтальный однородный изотропный пласт круглой формы (в плане). Размерность сетки фильтрационной модели — 21x21x1. Радиус круга (по центрам ячеек) — 1000 м. Размер ячейки е плане — 100 м. Мощность пласта — 10 м. Пористость — 20 %. проницаемость — 0,01 мкм2.

В.4.2 Свойства флюидов

Модель двухфазная: нефть, газ.

Свойства флюидов и зависимости PVT-свойсгв от давления представлены в таблице В.4.

В.4.3 Относительные фазовые проницаемости

Остаточная нефгенасыщенность {при вытеснении нефти газом) равна 10 %. ОФП заданы квадратичными зависимостями:

где 5— газонасыщенность.


[ 0,1*«


|Г, гчм**0А


(В.29)



(В.30)


FMPerni


Рисунок В.5 — Относительные фазовые проницаемости В.4.4 Начальные и граничные условия

В граничных ячейках либо за счет увеличения поровото объема, либо иным способом задают условие постоянного давления {контур питания). Начальное пластовое давление — 15 МПа. Вертикальная добывающая скважина расположена в центре модели в ячейке (11.11,1) и вскрывает перфорацией всю толщину пласта. Забойное давление — 3 МПа. Расчет выполнен на один год. Дебит нефти на юнец расчета считают установившимся.

^И*11 H4I ЮфОИ ШЬ®

Н■ мпивгуог 2ооо. si«o2 и oe«v*>

^-т—

«00 2tO«e »00«S ?»«S «00*S


Рисунок В.6 — Начальные и граничные условия

Таблица В.4 — Параметры модели теста 4

Параметр

Значение или зависимость, здесь р — давление. 0.1 МПа (бар)

Мощность пласта h, м

10

Пористость. д.ед.

0.2

Проницаемость К, мкм2

0.01

Радиус контура питания Рс. м

1000

Радиус скважины rw м

0.108

Давление на контуре питания р^ МПа

30

Давление точки росы р^, МПа

23

Забойное давление р^ МПа

15

Объемный коэффициент нефти, предельно насыщенной газом В0. д.ед

5.59Ю-*р1в7

Вязкость нефти, предельно насыщенной газом це. мПа-с

0.162-р-0'000476

Гаэосодержание Rx нефти, предельно насыщенной газом. м33

з.овю-з-р236

Объемный коэффициент газа, предельно на съеденного нефтью В7 д.ед.

0.392Р*0-846

Объемный коэффициент газа, кедонасыщенного нефтью с равным начагъному, д.ед.

0.392 р^-^/р

Вязкость газа, предельно насыщенного нефтью. р9< мПа-с

8.66 Ю^р’6

Вязкость газа, недонасыщечного нефтью, с Ry. равным начальному. мПа-с

e.ee-IO^Pt,16^! +0.001(р-ре))

Плотность разгззированной нефти, кг/м3

780

Молярная масса нефти, г/моль

132

Молярная масса газа, г/моль

25.83

В.4.5 Эталонные результаты моделирования

Анагытическое значение дебита газа находится тем же методом и по аналогичным формулам, что и 8 тесте 3 (подраздел В.З). С учетом принятых зависимостей свойств от давления получают следующее аналитическое значение дебита газа: 300,308 тыс. м3/сут.

В.4.6 Критерии оценки

Оценивают отклонение расчетного значения дебита газа через один год после начала моделирования от эталонного значения 300,308 тыс. м3/сут.

Допустимая относительная погрешность — 5 %.

• • ‘ s ?

I|§

ill8 8

§ 1 i s 3 |

* 8

5

-

S W

II

II 1

* § в

Ф э

® 3

к>

Моделирование пластового флюида

и>

Расчет равновесного начального состояния модели

*

Задание условий иа транииах расчетной области

VI

Оценка начальных запасов УВ

о

Контроль непротиворечивости модели

Если применимо

ы

Адаптация мелели к истории разработки

с»

Выполнение прогнозных расчетов

Ф

Параллельные вычисления

О

Учет притока х тазовым скважинам

-

Возможность детализации сетки а охолосхаажиииой зоне

i

1Ч>

Моделирование закачки полимеров

i

W

Моделирование закачки С02

*

Моделирование еодогаэоеого воздействия — моделирование поочередной закачки воды и газа, а также совместной одновременной закачки

VI

Моделирование меизотермической фильтрации

Если применимо

о»

Учет течения флюида в стволах ГС

Если применимо

ы

Моделирование ГРП

.

Композиционное моделирование

.

ф

Псевдокомпоэиционное моделирование

Если применимо

ы

о

Модели двойной пористости/даойиой проницаемости

I

N>

Моделирование добычи метена из угольных пластов

1

м

Моделирование добычи сланцевых тазов

к>

и

Визуализация и анализ результатов моделирования

к

Совместное использование с системами моделирования течения в скважинах и поверхностных сетях сбора и транспорта

ю

VI

Адаптация и расчет неопределенности показателей разработки месторождения

0)

Ol

a

s

c

0)

G

2


g

Q>

a

S

s

s

?

s

x

®


3.


Матрица применения технологий программного обеспечения гидродинамического моделирования


w

Ss

li

gj

i

Моделирование нефтяных залежей

§

В

1

1

1

2

А

-

S * ж-'""" ♦ 0»

II $ lf I

i е / г

* о

*ч>

Моделирование пластового флюида

<•>

Расчет равновесного начального состояния модели

*

Задание условий на границах расчетной области

Оценка начальных запасов Ув

Контроль непротиворечивости модели

Если применимо

Адаптация модели к истории разработки

с»

Выполнение прогнозных расчетов

ю

Параллельные вычисления

1

г

Учет притока к газовым скважинам

-

Возможность детализации сетки в околосхеажииной зоне

1

t

с

Моделирование закачки полимеров

1

1

Моделирование закачки С02

1

1

*

Моделирование водогазового воздействия — моделирование поочередной закачки воды и таза, а также совместной одновременной закачки

1

1

и*

Моделирование неизотермнческой фильтрации

1

Если применимо

S

Учет течения флюида в стволах ГС

1

Моделирование ГРП

.

1

5

Композиционное моделирование

.

1

2

Псевдокомпозициоииое моделирование

Если применимо

1

к>

о

Модели двойной порисгост и'двойиой проницаемости

l

1

V

Моделирование добычи метана из угольных пластов

l

t

V

Моделирование добычи сланцевых газов

+

*>

<й>

Визуализация и анализ результатов моделирования

к

Совместное использование с системами моделирования течения в скважинах и поверхностных сетях сбора и транспорта

м

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


Продолжение таблицы Г1


ос


Водогазовое

воздействие

£

1

3

-

S V v*/

в Ъ -4-^

0

1 1 5

2 Е £

* * 1 ! ё / 1

*

+

N

Моделирование пластовою флюида

*

+

и

Расчет равновесною начального состояния модели

*

+

*

Задание условий иа границах расчетной области

4

Оценка начальных запасов УВ

4

Cl

Контроль непротиворечивости модели

Если применимо

Если применимо

«ч

Адаптация модели к истории разработки

+

4

Выполнение прогнозных расчетов

4

*

Параллельные вычисления

1

1

о

Учет притока к газовым скважинам

4

-

Возможность детализации сетки а окопоскважинной зоне

1

1

N

Моделирование эахачм полимеров

1

1

и

Моделирование закачки COj

1

Моделирование водогазовою воздействия — моделирование поочередной закачки воды и газа, а также совместной одновременной закачки

1

1

Моделирование иеиэотермической фильтрации

1

1

о

Учет течения флюида а стволах ГС

4

■ч

Моделирование ГРП

1

1

2

Композиционное моделирование

1

1

<0

Псевдокомпозициоиное моделирование

Если применимо

Если применимо

м

о

Модели двойной порисюсти/двоймой проницаемости

1

1

IS)

Моделирование добычи метана из угольных пластов

1

1

N

IS)

Моделирование добычи сланцевых газов

4

м

<й)

Визуализация и анализ результатов моделирования

4

4

%

Совместное использование с системами моделирования течения а скважинах и поверхностных сетях сбора и транспорта

4

4

N

О*

Адаптация и расчет и ««пределен нос? и показателей разработки месторождения


Продолжение таблицы Г. 1


ic


i = it I §

I?

Si

Is

3 s

1?

Закачка ПАВ. полимеров

£

й

§

й*

-

S *

II $ lf I

I е / г

* о

*ч>

Моделирование пластового флюида

<•>

Расчет равновесного начального состояния модели

*

Задание условий на границах расчетной области

Оценка начальных запасов Ув

Контроль непротиворечивости модели

Если применимо

Если применимо

Адаптация модели к истории разработки

с»

Выполнение прогнозных расчетов

ю

Параллельные вычисления

1

1

г

Учет притока к газовым скважинам

-

Возможность детализации сетки в околосхеажииной зоне

1

с

Моделирование закачки полимеров

1

Моделирование закачки С02

1

*

Моделирование водогазового воздействия — моделирование поочередной закачки воды и таза, а также совместной одновременной закачки

1

i*

Моделирование неизотермнческой фильтрации

1

I

S

Учет течения флюида в стволах ГС

Моделирование ГРП

.

5

Композиционное моделирование

.

.

2

Псевдокомпозициоииое моделирование

Если применимо

Если применимо

ь>

о

Модели двойной порисгост и'двойиой проницаемости

1

I

V

Моделирование добычи метана из угольных пластов

1

1

V

*>

Моделирование добычи сланцевых газов

*>

<й>

Визуализация и анализ результатов моделирования

к

Совместное использование с системами моделирования течения в скважинах и поверхностных сетях сбора и транспорта

А»

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


Продолжение таблицы Г1


гг


§1

*1

И

$ 1

**

ff

5‘

а

I

?

£

£

8

Первичные методы извлечения (моделирование истощения)

-

S V v*/

1 1 5

2 E £

8 s *

! £ / 1 • a

1

+

*

N

Моделирование пластовою флюида

1

+

*

и

Расчет равновесною начального состояния модели

1

+

*

*

Задание условий иа границах расчетной области

1

+

Оценка начальных запасов УВ

1

+

Cl

Контроль непротиворечивости модели

1

Если применимо

Если применимо

«ч

Адаптация модели к истории разработки

1

+

Выполнение прогнозных расчетов

1

+

*

Параллельные вычисления

1

1

о

Учет притока к газовым скважинам

1

+

-

Возможность детализации сетки а окопоскважинной зоне

1

1

1

N

Моделирование эахачм полимеров

1

1

1

и

Моделирование закачки COj

1

+

1

Моделирование водогазовою воздействия — моделирование поочередной закачки воды и газа, а также совместной одновременной закачки

1

1

1

Моделирование кеиэотермической фильтрации

1

1

1

о

Учет течения флюида а стволах ГС

1

+

■ч

Моделирование ГРП

1

1

2

Композиционное моделирование

1

+

<0

Псевдокомпозиционное моделирование

1

Если применимо

Если применимо

м

о

Модели двойной пористости‘двойной проницаемости

1

1

1

IS)

Моделирование добычи метана из угольных пластов

1

1

1

N

IS)

Моделирование добычи сланцевых газов

1

+

м

<й)

Визуализация и анализ результатов моделирования

1

%

Совместное использование с системами моделирования течения а скважинах и поверхностных сетях сбора и транспорта

1

+

4

N

О*

Адаптация и расчет и ««пределен нос? и показателей разработки месторождения


Продолжение таблицы Г. 1


се

Моделирование залежей тяхелых нефтей и природных битумов

Газовые залежи в сланцах

Добыча метана из угольных пластов

-

S *

4» 0»

II $ lf I

i е / i 11 / %

* О

1

+

Моделирование пластового флюида

1

+

<•>

Расчет равновесного начального состояния модели

1

+

*

Задание условий на границах расчетной области

1

*

4*

Оценка начальных запасов Ув

1

*

0t

Контроль непротиворечивости модели

1

Если применимо

Если применимо

Адаптация модели к истории разработки

1

+

с»

Выполнение прогнозных расчетов

1

+

ю

Параллельные вычисления

1

1

+

5

Учет притока к газовым скважинам

1

+

-

Возможность детализации сетки в околоскеажимной зоне

1

1

1

С

Моделирование мкачхи полимеров

1

Моделирование закачки С02

1

1

1

*

Моделирование водогазового воздействия — моделирование поочередной закачки воды и газа, а также совместной одновременной закачки

1

1

i*

Моделирование неизотермнческой фильтрации

1

|

1

S

Учет течения флюида а стволах ГС

1

+

ы

Моделирование ГРП

1

1

5

Композиционное моделирование

1

.

1

2

Псеадокомлозиционное моделирование

1

Если применимо

Если применимо

ы

о

Модели двойной лористости'двойной проницаемости

1

1

+

V

Моделирование добычи метана из угольных пластов

1

1

V

Моделирование добычи сланцевых газов

1

+

hi

4й>

Визуализация и анализ результатов моделирования

1

к

Совместное использование с системами моделирования течения в скважинах и поверхностных сетях сбора и транспорта

1

+

hi

4*

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


Продолжение таблицы Г1


£

§

S

э j « 8 0) а

CJ

S

X

0)

if

ж

1

01

3

тэ

01

-

S V v*/

0 *о

1 1 5

2 Е £

8 в *

! 6 / 1 • а

4

4-

N

Моделирование пластовою флюида

4

4

и

Расчет равиоеесиою начального состояния модели

4

4

*

Задание условий иа границах расчетной области

4

Оценка начальных аапасое УВ

4

Cl

Контроль непротиворечивости модели

Если применимо

Если применимо

«ч

Адаптация модели к истории разработки

+

4

Выполнение прогнозных расчетов

+

4

*

Параллельные вычисления

1

1

о

Учет притока к газовым скважинам

+

-

возможность детализации сетки а окопоскважинной зоне

1

1

N

Моделирование закачки полимеров

1

1

и

Моделирование закачки COj

+

4

Моделирование водогазовою воздействия — моделирование поочередной закачки воды и газа, а также совместной одновременной закачки

4

4

Моделирование кеиэотермической фильтрации

1

1

о

Учет течения флюида а стволах ГС

1

1

■ч

Моделирование ГРП

4

2

Композиционное моделирование

1

1

<0

Псевдокомпозициомное моделирование

Если применимо

Если применимо

м

о

Модели двойной лорисюсти/двойной проницаемости

1

1

IS)

Моделирование добычи метана из угольных пластов

1

1

N

IS)

Моделирование добычи сланцевых газов

+

м

<й)

Визуализация и анализ результатов моделирования

+

4

%

Совместное использование с системами моделирования течения а скважинах и поверхностных сетях сбора и транспорта

+

4

N

О*

Адаптация и расчет и ««пределен нос? и показателей разработки месторождения


Скончание таблицы Г. 1


[1] РД 153-39.0-047—00

[2] РД 153-39.0-109—01


Библиография

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

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

УДК 004.45:006.354 ОКС 35.080 ОКП 42 5400

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

Редактор С А. Кузьмин Технический редактор В.Н. Прусакова Корректор Ю.М. Прокофьева Компьютерная верстка ЕА. Кондрашовой

Сдано в набор 01.03.2016. Подписано а печать 1S.03.2016 Формат 60*64%. Гарнитура Ариел.

Уел. леч. п. 4.65 Уч.-изд. п. 4.20 Тираж 35 из. Зак. 740.

Издано и отпечатано ео . 12399S Москва. Гранатный пер.. 4.