RU2709047C1 - Способ адаптации гидродинамической модели продуктивного пласта нефтегазоконденсатного месторождения с учетом неопределенности геологического строения - Google Patents

Способ адаптации гидродинамической модели продуктивного пласта нефтегазоконденсатного месторождения с учетом неопределенности геологического строения Download PDF

Info

Publication number
RU2709047C1
RU2709047C1 RU2019100279A RU2019100279A RU2709047C1 RU 2709047 C1 RU2709047 C1 RU 2709047C1 RU 2019100279 A RU2019100279 A RU 2019100279A RU 2019100279 A RU2019100279 A RU 2019100279A RU 2709047 C1 RU2709047 C1 RU 2709047C1
Authority
RU
Russia
Prior art keywords
development
geological
model
hydrodynamic
parameters
Prior art date
Application number
RU2019100279A
Other languages
English (en)
Inventor
Сергей Александрович Кирсанов
Александр Иосифович Пономарев
Анатолий Васильевич Меркулов
Тимур Владимирович Сопнев
Анатолий Кузьмич Арабский
Руслан Леонидович Кожухарь
Original Assignee
Общество с ограниченной ответственностью "Газпром добыча Ямбург"
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Общество с ограниченной ответственностью "Газпром добыча Ямбург" filed Critical Общество с ограниченной ответственностью "Газпром добыча Ямбург"
Priority to RU2019100279A priority Critical patent/RU2709047C1/ru
Application granted granted Critical
Publication of RU2709047C1 publication Critical patent/RU2709047C1/ru

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B43/00Methods or apparatus for obtaining oil, gas, water, soluble or meltable materials or a slurry of minerals from wells
    • EFIXED CONSTRUCTIONS
    • E21EARTH OR ROCK DRILLING; MINING
    • E21BEARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

Изобретение относится к способу адаптации гидродинамической модели с учетом неопределенности геологического строения. Техническим результатом является минимизация погрешности расчета технологических показателей разработки месторождения с применением гидродинамических моделей. Способ включает проведение геофизических, газодинамических, гидродинамических исследований скважин, отбор керна, проведение петрофизических исследований, обобщение материалов по изучению геологического строения, построение геологической модели месторождения, определение распределения фаций по площади месторождения по данным петрофизических и геофизических исследований, определение минимального, максимального и наиболее вероятного значений коэффициентов песчанистости, пористости, нефтегазонасыщенности, проницаемости, эффективных нефтегазонасыщенных толщин для каждой ячейки (блока) трехмерной модели месторождения, оценку геологических показателей в межскважинном пространстве, расчет показателей разработки на гидродинамической модели. В базу данных автоматической системы управления технологическими процессами (АСУ ТП) и/или информационно-управляющей системы (ИУС) вводят рассчитанные на гидродинамической модели показатели разработки и допустимые отклонения. Осуществляют контроль фактических показателей разработки по приборам, установленным на скважинах, и записывают результаты этих измерений в свою базу данных. Проводят проверку отклонения расчетных показателей от фактически измеренных. Формируют блок информации о фактических текущих и исторических показателях разработки, которые передаются по каналам связи на дополнительную обработку. Строят трехмерное распределение погрешности расчета адаптируемого показателя разработки по гидродинамической модели. Определяют аналитическую зависимость между геологическими параметрами, участвующими в адаптации, и адаптируемым показателем разработки. Определяют значения геологических параметров для каждой ячейки (блока) трехмерной модели. Для каждой ячейки проверяют соответствие значений заданных параметров. Если величина геологического параметра выходит за границы вероятностных отклонений, ему присваивают соответственно максимальное или минимальное вероятное значение. Продолжают процесс адаптации для других ячеек (блоков) модели до достижения заданной погрешности. Передают уточненные расчетные данные в базу данных АСУ ТП и/или ИУС для дальнейшего контроля за разработкой месторождения. В случае невозможности обеспечения заданной точности расчета показателей определяют геологические параметры с наибольшей степенью неопределенности и высоким влиянием на показатели разработки. Определяют зоны месторождения с максимальной погрешностью расчета показателей разработки и в выявленных зонах месторождения проводят дополнительные геофизические, петрофизические, гидродинамические исследования для локального уточнения геологических параметров, и по результатам дополнительных исследований проводят повторную процедуру адаптации гидродинамической модели, после чего соответствующие параметры загружают в базу данных для дальнейшего контроля за разработкой месторождения. 1 з.п. ф-лы, 2 ил., 1 табл.

Description

Изобретение относится к области добычи природного газа, нефти и конденсата, а именно к процедуре адаптации гидродинамических моделей для корректного прогноза технологических показателей разработки месторождений.
При построении гидродинамических моделей с историей разработки возникают отклонения расчетных показателей разработки от фактических показателей, наблюдаемых на месторождении. Применение математических моделей, не соответствующих фактическим показателям разработки недопустимо, поскольку ведет к принятию неверных технических решений при проектировании разработки и обустройстве месторождения.
Известен способ адаптации фильтрационных потоков к фактическим показателям разработки [Красовский А.В., Свентский С.Ю., Лысов А.О., Атеполихин В.В. «Комплексное геолого-технологическое моделирование крупного газового месторождения на примере Заполярного НГКМ» - сборник докладов XVII науч. - практ. конф. ТюменНИИгипрогаз «Проблемы развития газовой промышленности», 2014, С. 101-103]
Способ включает введение коэффициентов множителей порового объема или коэффициента проницаемости в определенных областях месторождения для снижения погрешности расчета показателей разработки на гидродинамической модели.
Существенным недостатком данного способа является то, что при его реализации не учитывают геологические особенности месторождения и в процессе адаптации могут возникнуть зоны месторождения с аномальными значениями пористости и проницаемости, которые будут выходить за границы максимальных и минимальных значений, полученных при петрофизических и геофизических исследованиях скважин и пластов.
Известен способ получения трехмерного распределения проницаемости пласта, в котором результаты гидродинамических исследований скважин используют для построения геологических и гидродинамических моделей (Патент RU №2479714).
Способ включает определение по результатам гидродинамических исследований скважин, вскрывающих пласт, осредненных по разрезу значений текущей фазовой проницаемости, пересчет текущей фазовой проницаемости на первоначальную проницаемость по нефти в присутствии остаточной воды с учетом термобарических условий в пласте и информации об относительных фазовых проницаемостях для каждой исследованной скважины, расчет зависимости изменения первоначальной проницаемости по разрезу с учетом результатов геофизических исследований в открытом стволе и профилей притоков пластовых флюидов, их корреляцию и построение трехмерного распределения проницаемости.
Представленный способ позволяет строить трехмерное распределение коэффициента проницаемости по данным гидродинамических исследований скважин.
Существенным недостатком указанного способа является то, что он не позволяет уточнить большинство геологических параметров, влияющих на фильтрацию нефти и газа в пласте. Также представленный способ не может использоваться для определения коэффициента проницаемости в межскважинном пространстве.
Наиболее близким по технической сущности к заявляемому является способ построения геологической и гидродинамической модели (Патент RU №2475646).
Суть известного способа заключается в том, что для построения геологической и гидродинамической моделей производят определение условий формирования пород по вещественному составу, а также по текстурным и структурным диагностическим признакам (литолого-фациальный анализ (ЛФА)), минералого-петрографический анализ осадочных пород исследуемого объекта, интерпретацию материалов геофизических исследований скважин (ГИС), обработку данных методами многомерной математической статистики.
Существенным недостатком известного способа является то, что он не позволяет уточнять показатели разработки и выявлять геологические особенности строения месторождений, для месторождения с историей разработки, по которым показатели разработки, рассчитанные по гидродинамической модели, существенно отклоняются от фактических показателей разработки.
Задачей, на решение которой направлено настоящее изобретение, является минимизация погрешности расчета технологических показателей разработки месторождения с применением гидродинамических моделей за счет снижения неопределенности геологического строения месторождения и оперативного реагирования на отклонение расчетных показателей разработки от фактических.
Техническим результатом, достигаемым от реализации изобретения, является повышение точности определения сроков и объема ввода производственных мощностей, снижение затрат на топливный газ дожимных компрессорных станций и химические реагенты, используемые при добыче и промысловой обработке скважинной продукции за счет точного прогноза динамики пластового давления и состава добываемого флюида; повышение конечного коэффициента извлечения УВС за счет корректного распределения добычи углеводородного сырья по скважинам и площадям месторождения.
Указанная задача решается, а технический результат достигается тем, что проводят геофизические исследования скважин и петрофизические исследования кернового материала. Собирают информацию обо всех гидродинамических исследованиях скважин. Собирают данные обо всех показателях разработки и всех работах, проводимых на скважинах за всю историю разработки. Рассчитывают погрешности определения коэффициентов песчанистости, пористости, нефтегазонасыщенности, проницаемости, эффективных нефтегазонасыщенных толщин, распределения фаций по скважинам.
Производят построение геологической модели месторождения. Построение геологической модели включает в себя этапы создания:
1) структурной модели;
2) литологической модели;
3) фациальной модели;
4) модели флюидальных контактов;
5) трехмерного распределения коэффициента песчанистости;
6) трехмерного распределения коэффициента пористости;
7) трехмерного распределения коэффициента насыщения;
8) трехмерного распределения коэффициента проницаемости.
Выбирают алгоритм построения геологической модели, обеспечивающий возможность производить многовариантное геологическое моделирование и другие математические операции, позволяющие оценить максимальное, минимальное и наиболее вероятное значение геологических параметров для каждой области трехмерного пространства месторождения с учетом погрешности определения геологических параметров по скважинным данным и с учетом погрешностей математического моделирования в межскважинном пространстве. Соотношение между минимальными и максимальными значениями и наиболее вероятными значениями геологических параметров в области трехмерного пространства будет сильно различаться в зависимости от степени удаленности точки пространства от фактических данных по скважинам. Вблизи скважин с выполненным комплексом ГИС неопределенность геологического строения будет, в первую очередь, обусловлена погрешностями определения исходных данных (ГИС, лабораторные исследования и т.д.). На удалении от скважин диапазон изменения геологических параметров будет больше, чем вблизи скважин, поскольку на погрешности определения параметров пластов в скважинах по комплексу петрофизических, геофизических и гидродинамических исследований накладывается неопределенность интерполяции/экстраполяции значений геологических параметров в межскважинном пространстве.
В гидродинамический симулятор загружают наиболее вероятные пространственно-распределенные значения геологических параметров. Используют данные о свойствах горных пород и насыщающих их флюидов, а также данные об истории разработки месторождения, проводят гидродинамический расчет. Добавляют в гидродинамическую модель плановые объемы добычи и рассчитывают прогнозные показатели разработки.
Прогнозные расчеты показателей разработки загружают в базу данных АСУ ТП и/или ИУС. АСУ ТП и/или ИУС с заданным шагом дискретизации регистрирует информацию с датчиков, установленных на скважинах, с показаниями давлений, расходов и составов добываемых флюидов. Эти данные передаются в базу данных АСУ ТП и/или ИУС с помощью телеметрии. В базу данных АСУ ТП и/или ИУС также загружают предельно допустимые отклонения каждого показателя разработки для каждой скважины/группы скважин. Базовые значения погрешностей расчета показателей разработки приведены во «Временном регламенте оценки качества и приемки трехмерных цифровых геолого-гидродинамических моделей, предоставляемых пользователями недр в составе технических проектов разработки месторождений углеводородного сырья на рассмотрение ЦКР Роснедр по УВС», принятом расширенным заседанием ЦКР Роснедр по УВС от 19.04.2012. При необходимости и соответствующем обосновании, в зависимости от решаемых производственных задач, связанных с проектированием разработки и обустройства месторождения, допускается установка более строгих ограничений на погрешности расчета показателей разработки исходя из требований нормативной документации, технологических и геологических особенностей месторождения. При определении допустимых отклонений также учитывают погрешность расчета фактических технологических показателей средствами телеметрии.
Параллельно с передачей данных телеметрии в базу данных, АСУ ТП и/или ИУС анализирует поступающие данные и определяет расхождение между расчетными и фактическими показателями разработки. При отклонении фактических значений параметров разработки от расчетных свыше предельно допустимого значения АСУ ТП и/или ИУС выдает оператору сообщение об ошибке и выгружает текущие и исторические данные о показателях разработки, которые передает по каналу связи для дальнейшей специальной обработки в программное обеспечение геологического и гидродинамического моделирования.
Используя эти данные, определяют чувствительность ключевых показателей разработки к геологическим параметрам. Составляют стратегию адаптации гидродинамической модели, в которой описано, какие показатели разработки какими геологическими параметрами будут адаптированы. По результатам геологического и гидродинамического моделирования определяют показатели разработки, характеризующиеся максимальной погрешностью расчета, а также максимальным влиянием на планирование разработки. Исходя из выбранных показателей разработки, выбирают геологические параметры для адаптации. Геологические параметры выбирают на основании данных о чувствительности показателя разработки к конкретному геологическому параметру и по степени неопределенности его расчета.
После каждой итерации производят сравнение рассчитанного значения показателя разработки с фактическим показателем, определяют погрешность его расчета и делают заключение о достижении требуемой погрешности расчета. При достижении требуемой погрешности расчета адаптацию заканчивают. Если требуемая погрешность не достигнута, производят модификацию следующего геологического параметра. Алгоритм адаптации гидродинамической модели представлен на фиг. 1.
Модификацию геологических параметров производят с использованием следующих трехмерных распределений:
1) погрешность показателя разработки;
2) минимально возможное значение геологического параметра;
3) максимально возможное значение геологического параметра;
4) наиболее вероятное значение геологического параметра.
Трехмерное распределение погрешности показателя разработки получают интерполяцией невязки показателя разработки, с высокой степенью дискретизации, определенной по работающему интервалу продуктивного пласта в скважинах. Для этого замеренные фактические показатели разработки интерполируют: пересчитывают по физическим формулам на глубину работающего интервала каждой из скважин и определяют значение показателя по всему интервалу продуктивного пласта/объекта разработки с шагом дискретизации, превышающим размерность геологической сетки. Рассчитывают невязку показателя разработки по работающему интервалу продуктивного пласта/объекта разработки в скважинах с высокой степенью дискретизации по формуле:
ΔT(MD)=Тзамер(MD)-Трасчет(MD),
где ΔТ(MD) - невязка показателя разработки по работающему интервалу продуктивного пласта в скважине;
Тзамер(MD) - замеренное значение показателя разработки по работающему интервалу продуктивного пласта в скважине;
Трасчет(MD) - рассчитанное по гидродинамической модели значение показателя разработки по работающему интервалу продуктивного пласта в скважине;
MD - глубина вдоль ствола скважины.
Затем по скважинным данным строят трехмерное распределение невязки показателя разработки. Трехмерное распределение невязки, построенное таким образом, учитывает не только изменение невязки по площади залежи, но и по разрезу. Это важно для многопластовых залежей и залежей с большим этажом нефтегазоносности.
Трехмерное распределение максимального и минимального значения каждого из геологических параметров получают при анализе геологической модели. Данные параметры получают, например, при помощи многовариантного моделирования. Полученные в результате многовариантного моделирования ансамбли геологических параметров анализируют и определяют допустимые диапазоны варьирования геологических параметров для каждой области трехмерного пространства месторождения. Таким образом, максимальные и минимальные значения геологических параметров будут определены отдельно для каждой области залежи и будут учитывать особенности геологического строения и геологической изученности месторождения.
В зависимости от геологических особенностей месторождения все этапы построения геологической модели разделяют на «влияющие» на результат построения других этапов моделирования и «не влияющие».
Для «зависимых» геологических параметров невозможно заранее определить максимальное и минимальное значение геологического параметра. Расчет трехмерных распределений максимального и минимального значения геологических параметров будет зависеть от математической связи между «зависимым» и «влияющим» параметрами. Связь «зависимого» и «влияющего» параметров бывает двух видов:
1) зависимый геологический параметр пересчитывают напрямую из «влияющего» параметра (например, пересчет коэффициента проницаемости по петрофизической формуле (корреляции) из коэффициента пористости);
2) «влияющий» геологический параметр используют при построении трехмерного распределения «зависимого» геологического параметра (например, расчет коэффициента нефтегазонасыщенности с использованием коэффициента пористости).
В первом случае трехмерное распределение максимального и минимального значения «зависимого» геологического параметра определяют путем расчета из «влияющего» параметра, и поэтому минимальное и максимальное возможное значение параметра отсутствует. Такой «зависимый» параметр используют для адаптации при существовании нескольких формул зависимостей, в которых определяют минимальные и максимальные возможные значения «зависимого» параметра при одинаковых значениях «влияющего» параметра. Во втором случае трехмерное распределение максимального и минимального значения зависимого геологического параметра определяют после проведения адаптации «влияющего» геологического параметра.
Типовая схема построения геологической модели представлена на фиг. 2. В зависимости от степени изученности месторождения и геологических особенностей месторождения, представленная схема может варьироваться для различных месторождений, в том числе и на одном месторождении на различных этапах его разработки. В данной схеме приведены основные связи между геологическими параметрами, позволяющие раскрыть суть реализации способа.
На предварительном этапе адаптации гидродинамической модели осуществляют: построение геологической модели, определение минимальных, максимальных и наиболее вероятных значений каждого пространственно-распределенного геологического параметра, расчет гидродинамической модели с использованием наиболее вероятных значений геологических параметров, определение погрешности расчета адаптируемых показателей разработки.
На первом этапе адаптации гидродинамической модели производят корректировку карты эффективных нефтегазонасыщенных толщин. Диапазон изменения эффективных нефтегазонасыщенных толщин зависит от структурного каркаса месторождения, глубин контактов флюидов и литологической модели месторождения. В зависимости от вклада каждого параметра в погрешность построения карты эффективных нефтегазонасыщенных толщин при уточнении карты эффективных нефтегазонасыщенных толщин в геологическую модель вносят соответствующие правки в структурную модель, модель насыщения и литологическую модель. В зависимости от условий построения геологической модели, контур залежи либо фиксируют, либо рассчитывают по результатам адаптации. При использовании фиксированного контура залежи, определяют линию изопахит, в границах которой производят уточнение карты эффективных нефтегазонасыщенных толщин. При расчете уточненной карты эффективных нефтегазонасыщенных толщин, соблюдают допустимый диапазон значений для каждой области месторождения. Расчет гидродинамической модели с каждым уточненным вариантом карты эффективных нефтегазонасыщенных толщин требует перестроения всех геологических параметров, поскольку карта эффективных толщин влияет на все последующее геологическое моделирование. Для снижения влияния от изменения фациальной модели, коэффициентов песчанистости, пористости, проницаемости и нефтегазонасыщенности на первом этапе адаптации, сохраняют средние значения коэффициентов песчанистости, пористости, проницаемости, нефтегазонасыщенности по месторождению и по площади залежи, а также сохраняют наиболее вероятное распределение фаций по площади залежи.
На втором этапе проводят уточнение фациальной модели месторождения. При наличии на месторождении разделения коллекторов на фации определяют распределение фаций по объему коллектора. В соответствии с уточненной картой эффективных нефтегазонасыщенных толщин рассчитывают соотношение фаций по площади месторождения. После того, как рассчитана уточненная литологическая и структурная модели, оценивают граничные значения содержания каждой фации по площади месторождения. Вместе с новым распределением фаций перестраивают все трехмерные распределения коэффициентов пористости, проницаемости и нефтегазонасыщенности с наиболее вероятными значениями, поскольку эти параметры варьируются от одной фации к другой. Настройку распределения фаций по площади месторождения реализуют, например, с использованием двумерных трендов распределения фации. Такое построение модели распространения фаций по объему месторождения позволяет гибко настраивать соотношение фаций по площади месторождения.
На третьем этапе адаптации проводят уточнение коэффициента песчанистости. Значение коэффициента песчанистости рассчитывают для каждой области месторождения исходя из наиболее вероятного значения и в диапазоне между минимальным и максимальным возможными значениями в конкретной области месторождения. Если построение трехмерного распределения коэффициента пористости или любого другого коэффициента, участвующего в расчете гидродинамической модели, зависит от коэффициента песчанистости, то проводят перестроение всех зависимых трехмерных распределений. Если коэффициент пористости и остальные параметры не зависят от коэффициента песчанистости, то уточнение коэффициента песчанистости проводят без изменения остальных геологических параметров.
На четвертом этапе производят расчет коэффициента пористости. Его расчет аналогичен расчету коэффициента песчанистости. Трехмерное распределение коэффициента нефтегазонасыщенности строят с учетом коэффициента пористости, поэтому проводят перестроение коэффициента нефтегазонасыщенности при каждой итерации расчета коэффициента пористости. Коэффициент проницаемости напрямую пересчитывают по корреляции с коэффициентом пористости, поэтому на каждой итерации пересчитывают трехмерное распределение проницаемости с учетом изменений пористости по объему залежи.
На пятом этапе проводят расчет трехмерного распределения коэффициента проницаемости. При этом учитывают то, что использование петрофизической зависимости лишает коэффициент проницаемости «гибкости» по отношению к коэффициенту пористости, поскольку для одного значения коэффициента пористости существует только одно возможное значение коэффициента проницаемости. При этом определяют погрешности расчета коэффициента проницаемости с этапа определения петрофизической формулы (корреляции). Допустимые минимальные и максимальные значения коэффициента проницаемости для одинаковых значений коэффициента пористости определяют по коэффициенту корреляции между коэффициентом пористости и коэффициентом проницаемости, определенным по керновым измерениям.
На шестом этапе проводят расчет коэффициента нефтегазонасыщенности. Методика расчета аналогична расчету коэффициента пористости и песчанистости.
Расчет модифицированных трехмерных распределений геологических параметров с использованием результатов расчета минимального и максимального значения геологического параметра, невязки показателя разработки при определении геологических параметров гидродинамической модели с наиболее вероятными их значениями, производят по формуле:
Figure 00000001
Figure 00000002
Figure 00000003
Figure 00000004
Figure 00000005
Figure 00000006
Figure 00000007
Figure 00000008
Figure 00000009
Figure 00000010
Figure 00000011
Figure 00000012
Figure 00000013
Figure 00000014
где Р - матрица пространственного распределения геологического параметра;
Pmin - матрица пространственного распределения минимальных геологических параметров во всех реализациях геологической модели;
Pmax - матрица пространственного распределения максимальных геологических параметров во всех реализациях геологической модели;
K1 - матрица пространственного распределения коэффициента, учитывающего погрешность расчета пластового давления;
K2 - матрица пространственного распределения коэффициента, учитывающего варьирование геологического параметра;
F - матрица пространственного распределения качественного параметра, погрешности расчета адаптируемого показателя разработки;
Tfact - матрица пространственного распределения фактического значения показателя разработки;
Trasc - матрица пространственного распределения расчетного значения показателя разработки;
Tnev - матрица пространственного распределения невязки расчетного значения показателя разработки;
х - порядковый номер реализации геологической модели;
r - общее количество реализаций геологической модели;
v - порядковый номер геологической реализации модели с наиболее вероятным геологическим строением, определенным экспертным путем;
i - номер столбца трехмерной сетки геологической модели;
j - номер строки трехмерной сетки геологической модели;
k - номер слоя трехмерной сетки геологической модели;
m - общее количество столбцов трехмерной сетки геологической модели;
n - общее количество строк трехмерной сетки геологической модели;
р - общее количество слоев трехмерной сетки геологической модели;
p(x)ijk - значение геологического параметра в ячейке геологической модели с координатами i;j;k, в реализации геологической модели номер х;
k1ijk - значение коэффициента, учитывающего погрешность расчета пластового давления, в ячейке геологической модели с координатами i;j;k;
k2ijk - значение коэффициента, учитывающего варьирование геологического параметра, в ячейке геологической модели с координатами
fijk - значение качественного параметра, погрешности расчета адаптируемого показателя разработки, в ячейке геологической модели с координатами i;j;k;
Figure 00000015
- значение фактического показателя разработки в ячейке геологической модели с координатами i;j;k;
Figure 00000016
- значение расчетного показателя разработки в ячейке геологической модели с координатами i;j;k;
Figure 00000017
- элемент матрицы Tdelta;
Figure 00000018
- элемент матрицы Tdelta с максимальным значением;
а - переменный адаптационный коэффициент;
b - переменный адаптационный коэффициент;
с - переменный адаптационный коэффициент;
d - переменный адаптационный коэффициент.
Коэффициент K2 определяет величину отклонения геологического параметра от наиболее вероятного значения в абсолютном выражении. Коэффициент K1 описывает закон определения доли коэффициента K2, на которую будет изменяться наиболее вероятное значение геологического параметра. Параметр F равен отношению погрешности адаптируемого показателя разработки в области трехмерного пространства к максимальной погрешности показателя разработки по залежи. В зависимости от значения параметра F залежь делится на 3 области. В первой области коэффициент K1 равен нулю, поэтому изменение базового геологического параметра не происходит. В данную область попадают участки месторождения с низкими погрешностями расчета гидродинамической модели. Во второй области коэффициент 0<K1<1, в данной области происходит плавная модификация базового геологического параметра. В третьей области коэффициент K1 равен 1. В данную область попадают участки залежи с максимальной погрешностью и поэтому значение модифицированного геологического параметра принимают равным минимальному/максимальному значению геологического параметра. В формуле (12) зависимость между коэффициентом K1 и погрешностью показателя разработки линейная. При различных геологических условиях используют и другие виды зависимости, например, в формуле (14) представлена квадратичная зависимость:
Figure 00000019
В зависимости от условий конкретного месторождения, используют различные виды зависимости между коэффициентом K1 и погрешностью показателя разработки: линейная, квадратичная, кубическая, логарифмическая и т.д.
Адаптацию гидродинамической модели производят подбором адаптационных коэффициентов а, b, с, d, для каждого геологического параметра. Влияние коэффициентов адаптации на геологические параметры показано в таблице 1:
Figure 00000020
Figure 00000021
Далее проводят серию расчетов гидродинамических моделей с различными адаптационными коэффициентами. Определяют значения адаптационных коэффициентов, при которых значения погрешности адаптируемого показателя разработки минимальны.
При достижении требуемой погрешности расчета адаптируемого показателя разработки, на гидродинамической модели рассчитывают уточненные прогнозные показатели разработки при плановых объемах добычи, которые загружают в базу данных АСУ ТП для дальнейшего контроля разработки месторождения.
Если же после процедуры адаптации гидродинамической модели погрешность расчета выходит за допустимые рамки, проводят дополнительные геофизические, петрофизические и гидродинамические исследования скважин. Для проведения исследований выбирают зоны с наиболее высоким расхождением расчетного и фактического показателей разработки месторождения. Выбирают геологический параметр для дополнительного исследования с учетом вероятностного диапазона изменения геологического параметра в конкретной области месторождения, а также с учетом влияния геологического параметра на адаптируемый показатель разработки.
После проведения дополнительных исследований, производят повторную процедуру адаптации гидродинамической модели с учетом новых базовых значений уточненных геологических параметров, а также с учетом нового диапазона вероятного значения геологического параметра.
Далее, при выполнении требований к точности расчета исторических показателей разработки на гидродинамической модели производят расчет прогнозных показателей разработки, которые загружают в АСУ ТП в качестве уточненных значений уставок для дальнейшего мониторинга разработки месторождения.
Применение данного способа позволяет повысить качество прогноза показателей разработки, месторождений нефти и газа с использованием гидродинамических моделей. В процессе адаптации гидродинамической модели происходит уточнение геологической модели и снижение неопределенностей геологических параметров по всему объему месторождения. Также заявленный способ позволяет увеличить точность подсчета запасов и расчета распределения запасов по площади месторождения. На основании уточненных данных о запасах УВС и прогнозах показателей разработки производят уточнение определения оптимального срока ввода необходимых производственных мощностей, прогноз объемов добычи флюидов, требуемых ГТМ и т.д.

Claims (1)

  1. Способ адаптации гидродинамической модели с учетом неопределенности геологического строения нефтегазоконденсатных месторождений, включающий: проведение геофизических, газодинамических, гидродинамических исследований скважин, отбор керна, проведение петрофизических исследований, обобщение материалов по изучению геологического строения, построение геологической модели месторождения, определение распределения фаций по площади месторождения по данным петрофизических и геофизических исследований, определение минимального, максимального и наиболее вероятного значений коэффициентов песчанистости, пористости, нефтегазонасыщенности, проницаемости, эффективных нефтегазонасыщенных толщин для каждой ячейки (блока) трехмерной модели месторождения, которые определяют на основании погрешностей геофизических и петрофизических исследований, а также с учетом вероятностного характера оценки геологических показателей в межскважинном пространстве, расчет показателей разработки на гидродинамической модели с использованием наиболее вероятных значений геологических параметров, оценку погрешности расчета показателей разработки на гидродинамической модели путем сравнения с фактическими данными, выбор показателей разработки месторождения, оказывающих наибольшее влияние на планирование добычи и вызывающих наибольшую погрешность ее расчета на гидродинамической модели, выбор геологических параметров, оказывающих наибольшее влияние на расчет показателей разработки месторождения и вызывающих наибольшую погрешность их расчета, регулярное сопоставление фактических показателей разработки с рассчитанными на гидродинамической модели, отличающийся тем, что с целью повышения точности расчета показателей разработки месторождения в базу данных автоматической системы управления технологическими процессами (АСУ ТП) и/или информационно-управляющей системы (ИУС) вводят рассчитанные на гидродинамической модели показатели разработки и допустимые отклонения рассчитанных показателей разработки месторождения по каждой скважине, после чего с заданной дискретизацией АСУ ТП и/или ИУС осуществляет контроль фактических показателей разработки по приборам, установленным на скважинах, и записывает результаты этих измерений в свою базу данных, одновременно АСУ ТП и/или ИУС проводит проверку отклонения расчетных показателей от фактически измеренных и при превышении предельно допустимых отклонений расчетных показателей разработки от фактических АСУ ТП и/или ИУС выдает сообщение оператору с указанием показателя разработки, для которого нарушено допустимое отклонение, и одновременно АСУ ТП и/или ИУС формирует блок информации о фактических текущих и исторических показателях разработки, который передает по каналам связи на дополнительную обработку в блок программ трехмерного геологического и гидродинамического моделирования для оперативного запуска процесса анализа сложившейся ситуации, в ходе которого строят трехмерное распределение погрешности расчета адаптируемого показателя разработки по гидродинамической модели, определяют аналитическую зависимость между геологическими параметрами, участвующими в адаптации, и адаптируемым показателем разработки, рассчитанным на гидродинамической модели, определяют значения геологических параметров для каждой ячейки (блока) трехмерной модели, позволяющие минимизировать погрешность расчета адаптируемого показателя разработки на гидродинамической модели, для каждой ячейки (блока) проверяют соответствие значений заданных параметров установленным для данной ячейки (блока) пределам вероятностной вариации значений соответствующих геологических параметров, величина которых обусловлена погрешностями петрофизических и геофизических измерений, погрешностями определения коэффициентов песчанистости, пористости, нефтегазонасыщенности, проницаемости, эффективной нефтегазонасыщенной толщины, распределением фаций в межскважинном пространстве, и, если величина геологического параметра выходит за границы вероятностных отклонений, ему присваивают соответственно максимальное или минимальное вероятное значение и продолжают процесс адаптации для других ячеек (блоков) модели до достижения заданной погрешности расчетных показателей истории разработки при соблюдении допустимых отклонений расчетных показателей разработки от фактических, определяют прогнозные показатели разработки и передают уточненные расчетные данные в базу данных АСУ ТП и/или ИУС для дальнейшего контроля за разработкой месторождения, а в случае невозможности обеспечения заданной точности расчета показателей разработки по результатам проведенной адаптации гидродинамической модели определяют геологические параметры с наибольшей степенью неопределенности и высоким влиянием на показатели разработки, определяют зоны месторождения с максимальной погрешностью расчета показателей разработки и в выявленных зонах месторождения проводят дополнительные геофизические, петрофизические, гидродинамические исследования для локального уточнения геологических параметров, и по результатам дополнительных исследований проводят повторную процедуру адаптации гидродинамической модели, после чего соответствующие параметры загружают в базу данных АСУ ТП и/или ИУС для дальнейшего контроля за разработкой месторождения.
RU2019100279A 2019-01-09 2019-01-09 Способ адаптации гидродинамической модели продуктивного пласта нефтегазоконденсатного месторождения с учетом неопределенности геологического строения RU2709047C1 (ru)

Priority Applications (1)

Application Number Priority Date Filing Date Title
RU2019100279A RU2709047C1 (ru) 2019-01-09 2019-01-09 Способ адаптации гидродинамической модели продуктивного пласта нефтегазоконденсатного месторождения с учетом неопределенности геологического строения

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
RU2019100279A RU2709047C1 (ru) 2019-01-09 2019-01-09 Способ адаптации гидродинамической модели продуктивного пласта нефтегазоконденсатного месторождения с учетом неопределенности геологического строения

Publications (1)

Publication Number Publication Date
RU2709047C1 true RU2709047C1 (ru) 2019-12-13

Family

ID=69006525

Family Applications (1)

Application Number Title Priority Date Filing Date
RU2019100279A RU2709047C1 (ru) 2019-01-09 2019-01-09 Способ адаптации гидродинамической модели продуктивного пласта нефтегазоконденсатного месторождения с учетом неопределенности геологического строения

Country Status (1)

Country Link
RU (1) RU2709047C1 (ru)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112485174A (zh) * 2020-10-19 2021-03-12 中国地质大学(北京) 基于堆叠立方体模型计算含水合物储层渗透率的方法
RU2754741C1 (ru) * 2021-03-12 2021-09-07 Общество с ограниченной ответственностью «Газпромнефть Научно-Технический Центр» (ООО «Газпромнефть НТЦ») Способ адаптации геолого-гидродинамической модели пласта
WO2022108490A1 (ru) * 2020-11-18 2022-05-27 Общество С Ограниченной Ответственностью "Газпромнефть Научно-Технический Центр" (Ооо "Газпромнефть Нтц) Способ локализации остаточных запасов
CN115375867A (zh) * 2022-10-24 2022-11-22 山东省地质调查院(山东省自然资源厅矿产勘查技术指导中心) 利用网格模型计算地热资源量方法、系统、设备及介质
CN116029136A (zh) * 2023-01-10 2023-04-28 中国矿业大学 一种非均质含水层中渗透系数随机场与孔隙度随机场协同构建方法
RU2809029C1 (ru) * 2023-06-27 2023-12-06 Общество с ограниченной ответственностью "Газпром добыча Ямбург" Способ уточнения фильтрационно-емкостных свойств пласта в межскважинном пространстве путем численной адаптации гидродинамической модели на результаты масштабных газодинамических исследований
CN117575106A (zh) * 2024-01-15 2024-02-20 中国地质科学院地质力学研究所 一种煤层气井产气剖面预测方法、系统、电子设备及介质

Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2123582C1 (ru) * 1998-03-11 1998-12-20 Научно-производственное предприятие "Нефтегеотехнология" Способ разработки сложнопостроенной залежи нефти
US20020082815A1 (en) * 2000-12-22 2002-06-27 Isabelle Rey-Fabret Method for forming an optimized neural network module intended to simulate the flow mode of a multiphase fluid stream
RU2346148C1 (ru) * 2008-02-01 2009-02-10 Общество с Ограниченной Ответственностью "ТНГ-Групп" Способ разработки нефтяных или нефтегазоконденсатных месторождений на поздней стадии
US20090194274A1 (en) * 2008-02-01 2009-08-06 Schlumberger Technology Corporation Statistical determination of historical oilfield data
RU2475646C1 (ru) * 2011-08-17 2013-02-20 Министерство образования и науки РФ Государственное образовательное учреждение высшего профессионального образования "Уральский государственный горный университет" Способ построения геологической и гидродинамической моделей месторождений нефти и газа
RU2014146091A (ru) * 2014-11-17 2016-06-10 Государственное бюджетное образовательное учреждение высшего профессионального образования "Альметьевский государственный нефтяной институт" Способ оптимизации процесса разработки и выработки запасов нефтяных и газовых месторождений
RU2657917C1 (ru) * 2017-03-21 2018-06-18 Общество с ограниченной ответственностью "Газпром добыча Ямбург" Способ уточнения геолого-газодинамической модели газовой залежи по данным эксплуатации

Patent Citations (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2123582C1 (ru) * 1998-03-11 1998-12-20 Научно-производственное предприятие "Нефтегеотехнология" Способ разработки сложнопостроенной залежи нефти
US20020082815A1 (en) * 2000-12-22 2002-06-27 Isabelle Rey-Fabret Method for forming an optimized neural network module intended to simulate the flow mode of a multiphase fluid stream
RU2346148C1 (ru) * 2008-02-01 2009-02-10 Общество с Ограниченной Ответственностью "ТНГ-Групп" Способ разработки нефтяных или нефтегазоконденсатных месторождений на поздней стадии
US20090194274A1 (en) * 2008-02-01 2009-08-06 Schlumberger Technology Corporation Statistical determination of historical oilfield data
RU2475646C1 (ru) * 2011-08-17 2013-02-20 Министерство образования и науки РФ Государственное образовательное учреждение высшего профессионального образования "Уральский государственный горный университет" Способ построения геологической и гидродинамической моделей месторождений нефти и газа
RU2014146091A (ru) * 2014-11-17 2016-06-10 Государственное бюджетное образовательное учреждение высшего профессионального образования "Альметьевский государственный нефтяной институт" Способ оптимизации процесса разработки и выработки запасов нефтяных и газовых месторождений
RU2657917C1 (ru) * 2017-03-21 2018-06-18 Общество с ограниченной ответственностью "Газпром добыча Ямбург" Способ уточнения геолого-газодинамической модели газовой залежи по данным эксплуатации

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112485174A (zh) * 2020-10-19 2021-03-12 中国地质大学(北京) 基于堆叠立方体模型计算含水合物储层渗透率的方法
CN112485174B (zh) * 2020-10-19 2021-09-14 中国地质大学(北京) 基于堆叠立方体模型计算含水合物储层渗透率的方法
WO2022108490A1 (ru) * 2020-11-18 2022-05-27 Общество С Ограниченной Ответственностью "Газпромнефть Научно-Технический Центр" (Ооо "Газпромнефть Нтц) Способ локализации остаточных запасов
RU2754741C1 (ru) * 2021-03-12 2021-09-07 Общество с ограниченной ответственностью «Газпромнефть Научно-Технический Центр» (ООО «Газпромнефть НТЦ») Способ адаптации геолого-гидродинамической модели пласта
CN115375867A (zh) * 2022-10-24 2022-11-22 山东省地质调查院(山东省自然资源厅矿产勘查技术指导中心) 利用网格模型计算地热资源量方法、系统、设备及介质
CN115375867B (zh) * 2022-10-24 2023-01-31 山东省地质调查院(山东省自然资源厅矿产勘查技术指导中心) 利用网格模型计算地热资源量方法、系统、设备及介质
CN116029136A (zh) * 2023-01-10 2023-04-28 中国矿业大学 一种非均质含水层中渗透系数随机场与孔隙度随机场协同构建方法
RU2809029C1 (ru) * 2023-06-27 2023-12-06 Общество с ограниченной ответственностью "Газпром добыча Ямбург" Способ уточнения фильтрационно-емкостных свойств пласта в межскважинном пространстве путем численной адаптации гидродинамической модели на результаты масштабных газодинамических исследований
CN117575106A (zh) * 2024-01-15 2024-02-20 中国地质科学院地质力学研究所 一种煤层气井产气剖面预测方法、系统、电子设备及介质
CN117575106B (zh) * 2024-01-15 2024-04-09 中国地质科学院地质力学研究所 一种煤层气井产气剖面预测方法、系统、电子设备及介质

Similar Documents

Publication Publication Date Title
RU2709047C1 (ru) Способ адаптации гидродинамической модели продуктивного пласта нефтегазоконденсатного месторождения с учетом неопределенности геологического строения
CA2692425C (en) Method, program and computer system for scaling hydrocarbon reservoir model data
CA2640725C (en) Methods, systems, and computer-readable media for fast updating of oil and gas field production models with physical and proxy simulators
US20070055447A1 (en) Method for updating a geological reservoir model by means of dynamic data
US20030028325A1 (en) Method of constraining by dynamic production data a fine model representative of the distribution in the reservoir of a physical quantity characteristic of the subsoil structure
US8364447B2 (en) Method, program and computer system for conciliating hydrocarbon reservoir model data
US5142471A (en) Method for determining the pressure or stress of a geological formation from acoustic measurement
US20110308792A1 (en) Method for operating an oil pool based on a reservoir model gradually deformed by means of cosimulations
CN106126936B (zh) 一种致密低渗透储层裂缝有效性的综合评价方法
CN105093331B (zh) 获取岩石基质体积模量的方法
Korde et al. Probabilistic decline curve analysis in the permian basin using Bayesian and approximate Bayesian inference
CN111967677B (zh) 一种非常规资源甜点分布的预测方法及装置
CN110566196B (zh) 一种储层连通性分析方法
CN108765562B (zh) 基于三维地质模型的油气产能评价方法
NO20200978A1 (en) Optimized methodology for automatic history matching of a petroleum reservoir model with ensemble kalman filter
CN111488666A (zh) 一种气藏水平井产能预测模型建立、预测方法及装置
CN112765527A (zh) 页岩气资源量计算方法及系统
RU2754741C1 (ru) Способ адаптации геолого-гидродинамической модели пласта
Worthington et al. Optimizing the value of reservoir simulation through quality-assured initialization
CN108875109B (zh) 一种异常地层压力的预测方法及系统
CN108647417B (zh) 一种确定页岩气储层含气饱和度的简易方法
CN105893674A (zh) 采用全局协方差进行地质属性预测的方法
CN113189677B (zh) 一种三维油藏物性参数模型自动更新方法
CN109858196B (zh) 一种油气藏参数解释方法及系统
Bowie et al. Applying numerical RTA to public data: Enabling field-wide property calibration and improved public data EUR forecasts