СИСТЕМЫ И СПОСОБЫ ДЛЯ ОЦЕНКИ МОМЕНТОВ ПРОРЫВА ФЛЮИДА В МЕСТОНАХОЖДЕНИЯХ ДОБЫВАЮЩИХ СКВАЖИН
SYSTEMS AND METHODS FOR ASSESSMENT OF MOMENTS OF PENETRATION OF FLUID IN LOCATIONS OF PRODUCTION WELLS

ПЕРЕКРЕСТНАЯ ССЫЛКА НА РОДСТВЕННЫЕ ЗАЯВКИ
Нет.
ЗАЯВЛЕНИЕ ОТНОСИТЕЛЬНО ИССЛЕДОВАНИЯ, ФИНАНСИРУЕМОГО ИЗ ФЕДЕРАЛЬНОГО БЮДЖЕТА
Неприменимо.
ОБЛАСТЬ ТЕХНИКИ, К КОТОРОЙ ОТНОСИТСЯ ИЗОБРЕТЕНИЕ
Настоящее изобретение в целом относится к оценке моментов прорыва флюида в местонахождении добывающей скважины. Более конкретно, данное изобретение относится к оценке моментов прорыва флюида в местонахождении добывающей скважины на основании моделирования распространения флюида.
УРОВЕНЬ ТЕХНИКИ ИЗОБРЕТЕНИЯ
Разные системы и способы известны для оценки момента прорыва флюида в местонахождении добывающей скважины, включая в себя адаптацию модели (HM). Адаптация Модели (HM) является систематической процедурой изменения имитационной модели пласта для воспроизведения динамического отклика месторождения. В приложениях HM и адаптации модели месторождения к данным по добыче основными задачами являются: a) интеграция данных по добыче в модели месторождения; b) гибкость, рентабельность и вычислительная эффективность и c) полная загрузка динамических данных.
В последние десятилетие технология HM имела резкое развитие и получила значительное признание и расширение от традиционного (то есть ручного, детерминированного) подхода, главным образом основанного на стратиграфических способах, к новым разработкам, подобным вероятностному моделированию, HM на основе метода трубок тока, основанным на чувствительности/градиенте и экспериментальной модели.
Схема работы HM в значительной мере рассматривает минимизацию несовпадения между измеренным и смоделированным динамическим откликом флюида (например, нефти или воды) в отдельной добывающей скважине как одну из инверсионных основных целей. В исследованиях по увеличению нефтеотдачи пласта (EOR) посредством закачки воды, к примеру, несоответствие отклика представляет собой дифференциальные или интегральные кривые обводненности с двумя основными характеристиками: 1) момент прорыва флюида и 2) направленность и форма отклика. В то время как обе характеристики представляют собой важные переменные в процессе минимизации несоответствия, момент прорыва флюида является тем, что раскрывает наибольшее влияние на экономику производительности. Более того, интервал (то есть временные рамки) прорыва флюида всегда обременены неуверенностью, что делает усилия по оценке с наиболее возможной достоверностью еще более значимыми. Фактически, целесообразно в HM динамических данных скважины рассматривать момент прорыва как эффект первого порядка, колебания в кривой направленность/форма как эффект второго порядка, поскольку они в основном отражаются на условиях эксплуатации.
Несмотря на прогресс в технологии HM, это по-прежнему несомненно наиболее времязатратный аспект построения модели/имитационного моделирования, и схема работы HM сталкивается со многими проблемами, которые включают в себя:
i) нелинейные результаты между откликом производства и характеристиками пласта;
ii) неоднозначные решения, которые требуют определение некоторого вида ″однозначности″;
iii) относительное влияние ключевых параметров может не быть очевидным;
iv) условия не ограничены, а неопределенности и в переменных редко известны; и
v) данные по добыче могут быть ошибочно и изначально искажены.
СУЩНОСТЬ ИЗОБРЕТЕНИЯ
Настоящее изобретение таким образом отвечает вышеприведенным потребностям и преодолевает один или больше недостатков в предшествующем уровне техники, предоставляя системы и способы для оценки моментов прорыва флюида в местонахождении добывающей скважины на основании моделирования распространения флюида.
В одном варианте осуществления настоящее изобретение включает в себя способ для оценки момента прорыва флюида в добывающей скважине, на основании данных моделирования распространения флюида, содержащий: i) идентификацию данных отслеживания линий тока; ii) вычисление среднего времени в пути в каждой ячейке сети на основании данных отслеживания линий тока; iii) идентификацию кратчайшей или быстрейшей линии тока для данной добывающей скважины, используя среднее время в пути линии тока в каждой ячейке сети; iv) вычисление среднего времени пролета для кратчайшей или быстрейшей линии тока через каждую пройденную ячейку сети, используя процессор вычислительной машины, и v) оценку момента прорыва флюида в добывающей скважине, используя данные моделирования распространения флюида и среднее время пролета для кратчайшей или быстрейшей линии тока.
В другом варианте осуществления, настоящее изобретение включает в себя энергонезависимое устройство хранения программ, материально несущее исполняемые компьютером команды для оценки момента прорыва флюида в добывающей скважине. Команды исполняются для выполнения: i) идентификации данных отслеживания линий тока; ii) вычисления среднего времени в пути в каждой ячейке сети на основании данных отслеживания линий тока; iii) идентификацию кратчайшей или быстрейшей линии тока для данной добывающей скважины, используя среднее время в пути линии тока в каждой ячейке сети; iv) вычисление среднего времени пролета для кратчайшей или быстрейшей линии тока через каждую пройденную ячейку сети, используя [процессор вычислительной машины]; и v) оценки момента прорыва флюида в добывающей скважине, используя данные моделирования распространения флюида и среднее время пролета для кратчайшей или быстрейшей линии тока.
Дополнительные аспекты, преимущества и варианты осуществления данного изобретения станут очевидны специалистам в данной области техники из нижеследующего описания различных вариантов осуществления и сопутствующих чертежей.
КРАТКОЕ ОПИСАНИЕ ЧЕРТЕЖЕЙ
Настоящее изобретение описано ниже со ссылкой на прилагаемые чертежи, на которых схожие элементы обозначены схожими номерами позиций и на которых:
Фиг. 1 является блок-схемой, иллюстрирующей некоторый вариант осуществления способа для реализации настоящего изобретения.
Фиг. 2A иллюстрирует скорость и направление распространения флюида через широкий песчаный карман.
Фиг. 2B - скорость и направление распространения флюида через узкий песчаный карман.
Фиг. 3 иллюстрирует пример распространения флюида через песчаную фракцию фациальной модели во время начального этапа моделирования.
1
4
Фиг. 4A иллюстрирует синтезированную 2D модель проницаемости с 2500 ячейками сети (50×50) и пятиточечной сеткой скважин (1 нагнетательная скважина (I) и 4 добывающих скважины (P-P)).
Фиг. 4B иллюстрирует моделирование распространения флюида через 2D модель проницаемости с Фиг. 4A от нагнетательной скважины (I) относительно числа итераций (2500), которым было запущено моделирование.
Фиг. 5 иллюстрирует возможное распределение линий тока в пятиточечной сетке скважин на Фиг. 4B.
Фиг. 6 иллюстрирует время в пути линии тока вдоль ее длин дуги в заданной ячейке сети (i, j, k) 2D модели проницаемости.
1
Фиг. 7A иллюстрирует наблюдаемую (измеренную) кривую обводненности для добывающей скважины P на FIG. 4A.
2
Фиг. 7B иллюстрирует наблюдаемую (измеренную) кривую обводненности для добывающей скважины P на FIG. 4A.
3
Фиг. 7C иллюстрирует наблюдаемую (измеренную) кривую обводненности для добывающей скважины P на FIG. 4A.
4
Фиг. 7D иллюстрирует наблюдаемую (измеренную) кривую обводненности для добывающей скважины P на FIG. 4A.
Фиг. 8 является функциональной схемой, иллюстрирующей некоторый вариант осуществления системы для реализации настоящего изобретения.
ПОДРОБНОЕ ОПИСАНИЕ ПРЕДПОЧТИТЕЛЬНЫХ ВАРИАНТОВ ОСУЩЕСТВЛЕНИЯ
Объект патентования настоящего изобретения описан с конкретностью, однако само описание не предназначено ограничивать объем изобретения. Объект патентования, таким образом, может также быть осуществлен другими путями, включая в себя другие этапы или сочетания этапов, сходных с теми, что описаны в данном документе, в сочетании с другими современными или перспективными технологиями. Кроме того, хотя термин ″этап″ может быть использован в данном документе для описания разных элементов примененных способов, данный термин не должен быть интерпретирован как подразумевающий любой конкретный порядок среди или между разными этапами, раскрываемыми в данном документе, за исключением случаев, когда иное явно ограничено данным описанием по отношению к тому или иному порядку. В то время как настоящее изобретение может быть применено в нефтяной и газовой индустрии, оно не ограничивается этим и может также быть применено в других индустриях для достижения схожих результатов.
Настоящее изобретение включает в себя системы и способы для оценки моментов прорыва флюида в местоположениях добывающих скважин на основании моделирования распространения флюида. Настоящее изобретение включает в себя моделирование распространения флюида, являющегося обычно статическим и представляющего момент(ы) вторжения для флюида, закачанного в нагнетательную скважину(ы) для достижения добывающей скважины(-). Моделирование дает иметь возможность полное рассмотрение фациального моделирования, которое сохраняет контроль над литологической целостностью геологических моделей посредством прямого ограничения моделирования при фациальном распространении. Моделирование также сохраняет стохастичность распространения фронта флюида. Несмотря на статичную природу данного моделирования стохастическая выборка движения фронта флюида выполняется посредством использования равномерного распределения.
Настоящее изобретение преобразует момент(ы) вторжения флюида (заданные моделированием в единицах повторений) в интервал физического времени (заданный, например, в днях, неделях, месяцах…), являющийся совместимым с историей эксплуатации скважины. Таким образом, настоящее изобретение предоставляет новые возможности для быстрой оценки ценных характеристик производительности скважины быстрым и рентабельным способом. К примеру, быстрая и точная оценка момент(ов) прорыва флюида, относящаяся к отдельной модели месторождения, может быть достигнута до начала полной инверсии. Подобная оценка предоставит ценную информацию операторам скважины относительно динамики клапанов скважины, в особенности в EOR проектах с нагнетанием воды/газа в которых управление производством нефти и воды/газа открывает существенный экономический эффект.
BT
С целью достижения быстрой оценки момента(ов) прорыва флюида (T), настоящее изобретение использует сочетание отслеживания линии тока и соответствующего времени пролета (″TOF″) по ходу моделирования. Таким образом, настоящее позволяет быструю аппроксимацию моментов прорыва флюида, выполняя запуск моделирования и одно повторение отслеживания линии тока в процессе автоматизированной адаптации модели (″AHM″), учитывающей линии тока, для модели месторождения.
ОПИСАНИЕ СПОСОБА
Со ссылкой на Фиг. 1, блок-схема иллюстрирует некоторый вариант осуществления способа 100 для реализации настоящего изобретения.
На этапе 102 осуществляется моделирование распространения флюида (″FPS″). Одна техника для выполнения FPS основывается на алгоритме в пакете программного обеспечения RGeoS, разработанного D. Renard. Данный алгоритм FPS моделирует распространение нескольких флюидов, известных в нагнетательных и/или добывающих скважинах, которые предусмотрены фациальной информацией, известной в узлах регулярной сети, и стремится позволить флюиду, встреченному в скважинах (например, в нагнетательной скважине) расти или расширятся пространственно. Скорость и направление роста зависит от размера песчаных карманов, которые могут быть заполнены. На Фиг. 2, к примеру, действительная скорость и направление распространения флюида через широкий песчаный карман (Фиг. 2A) и узкий песчаный карман (Фиг. 2B) показаны. Чем больше карманы 206, 208, там быстрее рост. Векторы скорости 202, 204 используются в алгоритме FPS. Данный алгоритм разработан для выполнения некоторого моделирования числовой переменной, используя технику моделирования Eden. Данная техника предоставляет более быстрое альтернативное решение для программы моделирования потока многофазного флюида. Данная техника объединяет пример двойной среды ″черное и белое″, где белое изображает песок, а черное изображает глину с одной или несколькими нагнетательными скважинами и одной или несколькими добывающими скважинами, как показано на Фиг. 3. В этом примере, показаны местоположения песчаных фаций 302, 304, 306 и двух нагнетательных скважин 307, 308.
1
4
Теперь, со ссылкой на Фиг. 4A, синтезированная 2D модель проницаемости показана с 2500 ячейками сети (50×50) и пятиточечной сеткой скважин (1 нагнетательная скважина (I) и 4 добывающих скважины (P-P)). Данный алгоритм FPS был использован в 2500 повторениях, поскольку одна ячейка данной модели заполняется за одно повторение. На Фиг. 4B моделирование распространения флюида через 2D модель распространения на Фиг. 4A от нагнетательной скважины (I) показано относительно числа повторений (2500), которым было запущено данное моделирование. На Фиг. 5 показано одно возможное распространение линий тока в пятиточечной сетке скважин на Фиг. 4B.
С целью реализации алгоритма FPS в качестве быстрой доверенной оценки момента прорыва флюида в инверсии AHM в кривых обводнения, преобразования момента(ов) вторжения флюида в интервал физического времени(-) должен быть рассмотрен со следующими основными допущениями:
i) TOF линии тока представляет собой ключевой фактор нормализации;
ii) отслеживание TOF от добывающей скважины(-) представляет собой дренируемый скважиной объем; и
iii) отслеживание флюида от нагнетательной скважины дает оценку рабочего объема.
Для оценки момента прорыва флюида в добывающей скважине предполагается, что следующие расчеты завершены для данной модели месторождения, используя любую технику, хорошо известную в данной области техники, для отслеживания линий тока по результатам прямого моделирования давления и скорости флюида: a) расчет момента вторжения флюида (то есть этап 102); и b) первое повторение отслеживания линии тока и расчета TOF (то есть этап 106). Эти расчеты дадут a) момент вторжения флюида от алгоритма FPS, заданного числом повторений моделирования (принимая 1 повторение на ячейку сети); и b) общее число линий тока, проходящих через любую ячейку сети модели месторождения с координатами (i, j, k).
m
На этапе 104 идентифицируют результаты данных FPS из этапа 102, которые включают в себя момент вторжения флюида, заданный числом повторений моделирования, необходимым флюиду для достижения любой добывающей скважины (P) от нагнетающей скважины через одну или несколько ячеек сети, представляющих модель свойств пласта.
SLN
SLN
p
<mrow><mo>∂</mo><mi>τ</mi></mrow>
<mrow><mi>ψ</mi><mtable><mtr><mtd><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></mtd></mtr></mtable></mrow>
<mrow />
На этапе 106 идентифицируют данные отслеживания линий тока, используя любую хорошо известную технику, которые включают в себя число сегментов линий тока, пересекающих каждую ячейку сети (N), время пробега (
) для каждого сегмента линий тока (
) в каждой ячейке сети, индекс ячейки сети и общее число ячеек сети, пересекаемых всеми линиями тока, соединяющими нагнетательную скважину с добывающей скважиной.
Со ссылкой на Фиг. 6, показано время пробега линии тока вдоль по ее длине дуги в данной ячейке сети 2D модели проницаемости. Индексы (n) и (m) пробегают все сегменты линий тока в каждой ячейке сети и во всех добывающих скважинах, соответственно (n=[1..N], а m=[l..N]). Время пробега для сегмента линий тока в каждой ячейке сети может быть вычислено посредством интегрирования ″медленности″ индикатора линии тока вдоль каждой траектории линии тока, используя следующее уравнение:
<mrow><mo>∂</mo><mi>τ</mi><mrow><mo>(</mo><mrow><mi>ψ</mi><mtable><mtr><mtd><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></mtd></mtr></mtable></mrow><mo>)</mo></mrow><mo>=</mo><mstyle displaystyle="true"><mrow><munder><mo>∫</mo><mi>ψ</mi></munder><mrow><mo>∂</mo><mi>s</mi><mrow><mo>(</mo><mi>x</mi><mo>)</mo></mrow></mrow></mrow></mstyle><mi>d</mi><mi>r</mi></mrow>
(1)
<mrow><mo>∂</mo><mi>s</mi><mo stretchy="false">(</mo><mi>x</mi><mo stretchy="false">)</mo></mrow>
<mrow><mi>d</mi><mi>r</mi></mrow>
<mrow><mi>ψ</mi><mtable><mtr><mtd><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></mtd></mtr></mtable></mrow>
где
соответствует ″медленности″ индикатора линии тока (определяемой как обратное от скорости индикатора),
соответствует длине дуги сегмента линии тока (
) между местами входа и выхода на ограничивающей поверхности ячейки сети с координатами (i, j, k).
<mrow><mo>∂</mo><mover><mi>τ</mi><mo>~</mo></mover></mrow>
На этапе 108 среднее время пробега линии тока в каждой ячейке сети (
) вычисляют, учитывая все сегменты линий тока, пересекающих каждую ячейку сети, что может быть вычислено, используя следующее уравнение:
<mrow><mo>∂</mo><mover accent="true"><mi>τ</mi><mo>˜</mo></mover><mo>=</mo><mfrac><mn>1</mn><mrow><msub><mi>N</mi><mrow><mi>S</mi><mi>L</mi><mi>N</mi></mrow></msub></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>n</mi><mo>=</mo><mn>1</mn></mrow><mrow><msub><mi>N</mi><mrow><mi>S</mi><mi>L</mi><mi>N</mi></mrow></msub></mrow></munderover><mrow><mo>∂</mo><mi>τ</mi><mrow><mo>(</mo><mrow><mi>ψ</mi><mtable><mtr><mtd><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></mtd></mtr></mtable></mrow><mo>)</mo></mrow></mrow></mstyle></mrow>
(2)
SLN
<mrow><mo>∂</mo><mi>τ</mi><mrow><mo>(</mo><mrow><mi>ψ</mi><mtable><mtr><mtd><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></mtd></mtr><mtr><mtd><mrow><mi>m</mi><mo>,</mo><mi>n</mi></mrow></mtd></mtr></mtable></mrow><mo>)</mo></mrow></mrow>
где (N) является числом сегментов линий тока, пересекающих каждую ячейку сети из этапа 106, а
является временем пробега для каждого сегмента линий тока в каждой ячейке сети из этапа 106.
m
m
<mrow><mo>∂</mo><msup><mover accent="true"><mi>τ</mi><mo>˜</mo></mover><mrow><mi>min</mi></mrow></msup></mrow>
На этапе 114 кратчайшую/быстрейшую линию тока идентифицируют для каждой добывающей скважины (P), используя среднее время пробега линии тока в каждой ячейке сети из этапа 108 и любой хорошо известный алгоритм поиска. Кратчайшая/быстрейшая линия тока является линией тока с наименьшей суммой средних времен пробега линий тока (
) в ячейках сети, пересекаемых линией тока между нагнетательной скважиной (I) и добывающей скважиной (P).
<mrow><msub><mover accent="true"><mi>N</mi><mo>^</mo></mover><mrow><mi>G</mi><mi>C</mi></mrow></msub><msup><mrow /><mrow><mi>min</mi></mrow></msup></mrow>
На этапе 116 общее число всех ячеек сети (
), пересекаемых кратчайшей/быстрейшей линией тока, идентифицированное на этапе 114, и их индексы из этапа 116 сохраняют.
min
<mrow><mo>∂</mo><msup><mover accent="true"><mi>τ</mi><mo>˜</mo></mover><mrow><mi>min</mi></mrow></msup></mrow>
<mrow><msub><mover accent="true"><mi>N</mi><mo>^</mo></mover><mrow><mi>G</mi><mi>C</mi></mrow></msub><msup><mrow /><mrow><mi>min</mi></mrow></msup></mrow>
На этапе 118 среднее TOF (<TOF>) для кратчайшей/быстрейшей линии тока, идентифицированной на этапе 114, вычисляют по каждой пересекающей ячейку сети, используя наименьшую сумму средних времен пробега линий тока (
) для кратчайшей/быстрейшей линии тока, идентифицированной на этапе 114, и общее число всех ячеек сети (
), сохраненное на этапе 116, которое может быть вычислено, используя следующее уравнение:
<mrow><msup><mrow><mrow><mo>〈</mo><mrow><mi>T</mi><mi>O</mi><mi>F</mi></mrow><mo>〉</mo></mrow></mrow><mrow><mi>min</mi></mrow></msup><mo>=</mo><mfrac><mn>1</mn><mrow><msubsup><mover accent="true"><mi>N</mi><mo>^</mo></mover><mrow><mi>G</mi><mi>C</mi></mrow><mrow><mi>min</mi></mrow></msubsup></mrow></mfrac><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>u</mi><mo>=</mo><mn>1</mn></mrow><mrow><msubsup><mover accent="true"><mi>N</mi><mo>^</mo></mover><mrow><mi>G</mi><mi>C</mi></mrow><mrow><mi>min</mi></mrow></msubsup></mrow></munderover><mrow><mo>∂</mo><msubsup><mover accent="true"><mi>τ</mi><mo>˜</mo></mover><mi>u</mi><mrow><mi>min</mi></mrow></msubsup></mrow></mstyle></mrow>
(3)
m
2
3
1
4
где индекс (u) представляет собой число пробегов по всем индексам ячеек сети, пересекаемых кратчайшей/быстрейшей линией тока. Различие между ″быстрейшей″ и ″медленнейшей″ линией тока от распределения линий тока, относящихся к каждой добывающей скважине (P), соответствует отличию между однородным и неоднородным пространственным распределением свойств пласта, таких как, к примеру, пустоты. Разница между распределением линий тока на Фиг. 5 показывает, что добывающие скважины P и P соединены с нагнетательной скважиной (I) через явно разные геологические формации, чем добывающие скважины P и P, что может быть сопоставленным с подстилающей структурой пустот.
На этапе 120 способ 100 определяет все ли ячейки сети, пересекаемые кратчайшей/быстрейшей линией тока, были рассмотрены. Если все пересекаемые ячейки сети не были рассмотрены, то тогда способ 100 возвращается к этапу 118. Если все пересекаемые ячейки сети были рассмотрены, то тогда способ 100 переходит к этапу 124. В качестве варианта, этапы со 118 по 120 могут быть выполнены единовременно для каждой пересекаемой ячейки сети.
На этапе 124 оценку момента прорыва флюида для каждой добывающей скважины (Pm) вычисляют, объединяя данные отслеживания линий тока из этапа 106 с данными FPS из этапа 104, что может быть вычислено, используя следующее уравнение:
<mrow><msup><mi>T</mi><mrow><mi>B</mi><mi>T</mi></mrow></msup><mo>=</mo><msup><mrow><mrow><mo>〈</mo><mrow><mi>T</mi><mi>O</mi><mi>F</mi></mrow><mo>〉</mo></mrow></mrow><mrow><mi>min</mi></mrow></msup><mo>×</mo><mfrac><mrow><msubsup><mi>t</mi><mrow><mi>I</mi><mi>N</mi><mi>V</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msubsup></mrow><mrow><msub><mi>N</mi><mi>p</mi></msub></mrow></mfrac><mo>×</mo><mfrac><mrow><msubsup><mi>N</mi><mrow><mi>S</mi><mi>L</mi><mi>N</mi></mrow><mi>m</mi></msubsup></mrow><mrow><msub><mi>N</mi><mrow><mi>x</mi><mi>y</mi><mi>z</mi></mrow></msub></mrow></mfrac></mrow>
(4)
xyz
p
m
min
<mrow><msubsup><mi>N</mi><mrow><mi>S</mi><mi>L</mi><mi>N</mi></mrow><mi>m</mi></msubsup></mrow>
<mrow><msubsup><mi>t</mi><mrow><mi>I</mi><mi>N</mi><mi>V</mi></mrow><mrow><mi>i</mi><mo>,</mo><mi>j</mi><mo>,</mo><mi>k</mi></mrow></msubsup></mrow>
где (N) и (N) представляют собой общий размер модели свойств пласта и общее число добывающих скважин, соответственно, (<TOF>) представляет собой среднее TOF для кратчайшей/быстрейшей линии тока, вычисленное на этапе 118, (
) представляют собой общее число ячеек сети, пересекаемых всеми линиями тока, соединяющими нагнетательную скважину (I) с добывающей скважиной (P), а (
)представляют собой время вторжения флюида из этапа 104.
m
m
m
На этапе 126 способ 100 определяет, все ли добывающие скважины были рассмотрены. Если не все добывающие скважины (P) были рассмотрены, то тогда способ 100 возвращается к этапу 104. Если все добывающие скважины (P) были рассмотрены, то тогда способ 100 завершается. В качестве варианта, этапы с 104 по 126 могут быть выполнены единовременно для каждой добывающей скважины (P).
ПРИМЕР
1
2
3
4
Со ссылкой на синтезированную 2D модель проницаемости на Фиг. 4A, наблюдаемые (измеренные) кривые обводнения для конфигурации модели на FIG 4A даны на Фиг. 7A, 7B, 7C и 7D для каждой из четырех добывающих скважин (P, P, P и P).
Точки данных дата/время по оси x на Фиг. 7A-7D, соответствующие физическим датам, относящимся к плану нагнетания воды (точки данных прорыва воды), представлены ниже в Таблице 1.
Таблица 1
Точка даных
Физическая дата (дд/мм/гггг)
1
17/9/2000
2
4/6/2001
3
19/2/2002
4
6/10/2002
5
24/7/2003
6
9/4/2004
7
25/12/2004
8
11/9/2005
B7
Наблюдаемые моменты прорыва воды, выведенные из Фиг. 4A, даны ниже в Таблице 2. Более того, Таблица 2 перечисляет моменты вторжения воды, вычисленные посредством алгоритма FPS, моменты прорыва воды (T), вычисленные с использованием способа, предложенного на Фиг. 1 и погрешность, соответствующую результату, полученному посредством способа, предложенного на Фиг. 1.
Таблица 2
Производитель
Наблюдаемое T<sup>BT</sup> (дней)
Время вторжения (повторений)
T<sup>BT</sup> предложенным способом (дней)
Погрешность (дней/%)
P1
263
2272
281,41
18,41/+7%
P2
121
2027
124,63
15,73/+3%
P3
263
2268
239,33
-23,67/-9%
P4
1043
2491
1105,58
62,58/+6%
Результаты показывают, что способ, предложенный на Фиг. 1, способен быстро предсказать момент прорыва флюида с погрешностью менее чем 10% для данной пятиточечной сетки скважин. Достигнутая погрешность может быть разной (больше/меньше), когда распространение флюида применено через область со значительно большей геологической сложностью, и динамическая модель объединяет значительно больше число добывающих скважин.
ОПИСАНИЕ СИСТЕМЫ
Настоящее изобретение может быть реализовано через исполняемые компьютером последовательности команд, такие как программные модули, обычно называемые программными приложениями или прикладными программами, исполняемыми компьютером. Программное обеспечение может включать в себя, для примера, подпрограммы, программы, объекты, компоненты, структуры данных и так далее, которые выполняют отдельные задачи или реализуют отдельные абстрактные виды данных. Интерактивная среда DecisionSpace, которая является коммерческим программным обеспечением, поставляемым Landmark Graphics Corporation, может быть использована как приложение интерфейса для реализации настоящего изобретения. Данное программное обеспечение может также взаимодействовать с другими сегментами кода для запуска различных задач в ответ на данные, принятые в сочетании с источником принятых данных. Программное обеспечение может быть сохранено и/или перенесено на любом виде запоминающих устройств, таких как CD-ROM, накопитель на магнитных дисках, запоминающее устройство на ЦМД и полупроводниковое запоминающее устройство (например, различные виды ПЗУ или ОЗУ). Кроме того, данное программное обеспечение и его результаты могут быть переданы по различным средам передачи, таким как оптоволокно, металлический провод, и/или через любую из множества сетей передачи данных, таких как Интернет.
Более того, специалисты в данной области техники оценят, что данное изобретение может быть применено со множеством конфигураций вычислительных устройств, включая ручные устройства, мультипроцессорные системы, микропроцессорная или программируемая бытовая электроника, мини-ЭВМ, основные ЭВМ, и тому подобное. Любое число вычислительных устройств и компьютерных сетей являются подходящими для использования с настоящим изобретением. Данное изобретение может быть применено в средах распределенных вычислений, где задачи выполняются дистанционно обрабатывающими устройствами, которые объединены через сеть связи. В среде распределенных вычислений программные модули могут быть расположены как на локальных, так и на удаленных запоминающих устройствах ЭВМ, включая накопительные запоминающие устройства. Настоящее изобретение может, таким образом, быть реализовано с использованием различного аппаратного обеспечения, программного обеспечения или их сочетания, в вычислительном устройстве или другой системе обработки.
Со ссылкой на Фиг. 8, функциональная схема иллюстрирует некоторый вариант осуществления системы для реализации настоящего изобретения на компьютере. Система включает в себя вычислительный блок, иногда названный вычислительной системой, который содержит запоминающее устройство, прикладные программы, клиентский интерфейс, видеоинтерфейс, и обрабатывающий модуль. Вычислительный блок является только некоторым примером подходящей вычислительной среды и не предназначен предлагать какие-либо ограничения относительно объема использования или функциональных возможностей данного изобретения.
Запоминающее устройство, в первую очередь, хранит прикладные программы, которые также могут быть описаны как программные модули, содержащие команды, исполняемые компьютером, выполняемые вычислительным блоком для реализации настоящего изобретения, описанного в настоящем документе и показанном на Фиг. 2. Запоминающее устройство таким образом, содержит модуль оценки момента прорыва флюида, который обеспечивает способ, показанный и описанный со ссылкой на Фиг. 1, и объединяет функциональные возможности из оставшихся прикладных программ, показанных на Фиг. 8. Модуль оценки момента прорыва флюида, к примеру, может быть использован для исполнения многих из функций, описанных со ссылкой на способ 100 на Фиг. 1. Интерактивная среда Decision Space может быть использована, к примеру, в качестве интерфейсного приложения для реализации модуля оценки момента прорыва флюида и использовать результаты способа 100 на Фиг. 1.
Хотя вычислительный блок показан как имеющий обобщенное запоминающее устройство, вычислительный блок обычно включает в себя различные машиночитаемые носители. В качестве примера, а не ограничения, машиночитаемый носитель может содержать запоминающее устройство ЭВМ. Запоминающее устройство ЭВМ может включать в себя запоминающее устройство ЭВМ в виде энергозависимого и/или энергонезависимого запоминающего устройства, такого как постоянное запоминающее устройство (ПЗУ) и оперативное запоминающее устройство (ОЗУ). Базовая система ввода/вывода (BIOS), содержащая базовые подпрограммы, которые помогают передавать информацию между элементами в вычислительном блоке, например, во время запуска, обычно содержится в ПЗУ. ОЗУ обычно содержит данные и/или программные модули, которые немедленно доступны и/или в настоящий момент обрабатываются обрабатывающим модулем. В качестве примера, а не ограничения, вычислительный блок включает в операционную систему прикладные программы, другие программные модули и данные программ.
Компоненты, показанные в запоминающем устройстве, могут также быть включены в другую съемную/несъемную энергозависимую/энергонезависимую среду компьютерного хранения или они могут быть реализованы в вычислительном блоке через прикладной программный интерфейс (″API″) или облачные вычисления, которые могут располагаться на отдельном вычислительном блоке, подключенном к вычислительному блоку через вычислительное устройство или сеть передачи данных. Только для примера, накопитель на жестком диске может читать с или записывать на несъемный энергонезависимый магнитный носитель, накопитель на магнитных дисках может читать с или записывать на съемный, энергонезависимый магнитный диск, а накопитель на оптических дисках может читать с или записывать на съемный, энергонезависимый оптический диск, такой как CD-ROM или другой оптический носитель. Другие съемные/несъемные, энергозависимые/энергонезависимые среды компьютерного хранения, которые могут быть использованы в примерной операционной среде, могут включать в себя, но не ограничиваться этим, кассеты с магнитной лентой, карты флэш-памяти, универсальные цифровые диски, цифровые видеоленты, твердотельные ОЗУ, твердотельные ПЗУ и так далее. Устройства и сопутствующие им среды компьютерного хранения, рассмотренные выше, обеспечивают хранение машиночитаемых команд, структур данных, программных модулей и других данных для вычислительного блока.
Пользователь может вводить команды и информацию в вычислительный блок через пользовательский интерфейс, который может быть устройствами ввода, такими как клавиатура и позиционирующее устройство, обычно называемый ″мышь″, шаровой указатель или сенсорная панель. Устройства ввода могут включать в себя микрофон, рычажный указатель, антенну спутниковой связи, сканер и тому подобное. Эти и другие устройства ввода часто подключают к обрабатывающему модулю через системную шину, но они могут быть подключены через другой интерфейс и шинные структуры, такие как порт параллельного ввода-вывода или универсальная последовательная шина (USB).
Монитор или другой вид устройства отображения может быть подключен к системной шине через интерфейс, такой как видеоинтерфейс. Графический интерфейс пользователя (″GUI″) может также быть использован с видеоинтерфейсом для приема команд от пользовательского интерфейса и передачи команд в обрабатывающий модуль. Дополнительно к монитору, компьютеры могут также включать в себя другие периферийные устройства вывода, такие как акустическая система и печатающее устройство, которые могут быть подключены через интерфейс выходных периферийных устройств.
Несмотря на то, что многие другие внутренние компоненты вычислительного блока не показаны, специалисты в данной области техники признают, что подобные компоненты и их схема соединений являются хорошо известными.
Хотя настоящее изобретение было описано применительно к настоящим предпочтительным вариантам осуществления, специалисты в данной области техники поймут, что это не предназначено для ограничения данного изобретения этими вариантами осуществления. Из этого следует, что предполагается, что различные альтернативные варианты осуществления и модификации могут быть сделаны к раскрытым вариантам осуществления, не покидая сущность и объем данного изобретения, определенные прилагаемой формулой изобретения и ее эквивалентами.