СПОСОБ МОДЕЛИРОВАНИЯ ПОДЗЕМНОГО ОБЪЕМА

1. Способ мониторинга изменений насыщенности подземного объема, содержащий этапы, на которых:
получают данные наблюдений поведения насыщенности из упомянутого подземного объема на протяжении периода времени;
используют одну или более моделей, получающих смоделированные данные поведения насыщенности из упомянутого подземного объема на протяжении периода времени;
преобразуют каждые из упомянутых данных наблюдений и упомянутых смоделированных данных таким образом, что в каждом случае данные описывают показателями:
параметра расположения фронта, в котором ячейка подземного объема связывается со значением в соответствии с ее расположением относительно фронта флюида, для которого выполняется мониторинг, и
параметра интенсивности вытеснения, в котором ячейка подземного объема связывается со значением в соответствии с наблюдаемой вариацией насыщенности на протяжении периода времени, или оценочной вариацией насыщенности на протяжении периода времени; и
вычисляют несоответствие между поведением насыщенности в упомянутых преобразованных данных наблюдений и поведением насыщенности в упомянутых преобразованных смоделированных данных на протяжении периода времени.
2. Способ по п. 1, в котором упомянутый параметр расположения фронта содержит бинарный индикаторный параметр, причем значение, привязанное к ячейке, зависит от того, находится ли она позади упомянутого фронта или впереди упомянутого фронта.
3. Способ по п. 2, в котором ячейка определяется как находящаяся позади упомянутого фронта, если она имеет определенный начальный уровень насыщенности выше порога насыщенности или определенную вариацию в уровне насыщенности выше порога вариации насыщенности; иначе она определяется, как находящаяся впереди упомянутого фронта.
4. Способ по п. 3, содержащий, для каждой обрабатываемой фазы, этап, на котором выбирают упомянутый порог вариации насыщенности, который ясно разделяет кластеры ячеек, не заполненных рассматриваемым флюидом, от ячеек, заполненных таким флюидом.
5. Способ по п. 3 или 4, содержащий этап, на котором для каждой обрабатываемой фазы выбирают упомянутый порог насыщенности.
6. Способ по любому из пп. 3-5, в котором там, где ячейка имеет определенную вариацию в уровне насыщенности выше упомянутого порога вариации насыщенности, упомянутый параметр интенсивности вытеснения имеет наблюдаемую вариацию насыщенности для соответствующей ячейки, и где ячейка не имеет определенной вариации в уровне насыщенности выше упомянутого порога вариации насыщенности, упомянутый параметр интенсивности вытеснения привязан к оценочной вариации насыщенности, ожидаемой для рассматриваемого флюида, если ячейка была перед упомянутым фронтом для рассматриваемого флюида.
7. Способ по п. 6, в котором упомянутую оценочную вариацию насыщенности выводят из среднего из наблюдаемых или смоделированных вариаций насыщенности на протяжении некоторых или всех ячеек позади упомянутого фронта, которые не находятся внутри начального расположения фронта для рассматриваемого флюида.
8. Способ по п. 6 или 7, в котором упомянутую оценочную вариацию насыщенности оценивают с использованием аналитической модели, которая обеспечивает возможность оценивать насыщенность позади фронта насыщенности как функции свойств породы и флюида.
9. Способ по любому из пп. 1-8, в котором параметр интенсивности вытеснения является непрерывным переменным параметром.
10. Способ по любому из пп. 1-9, в котором значение упомянутой оценочной вариации насыщенности определяют на начальном этапе.
11. Способ по любому из пп. 1-10, содержащий этап, на котором выполняют одно или более моделирований для получения упомянутых смоделированных данных.
12. Способ по любому из пп. 1-11, в котором упомянутый способ выполняют как часть итеративной методики адаптации по истории, содержащей выполнение инверсного процесса для минимизации рассчитанных несовпадений.
13. Способ по п. 12, в котором этап вычисления несовпадения содержит этап, на котором измеряют разности между соответствующими наблюдаемыми и смоделированными параметрами расположения фронта, и измерение разности между соответствующими наблюдаемыми и смоделированными параметрами интенсивности вытеснения.
14. Способ по п. 13, в котором измерение разности между соответствующими наблюдаемыми и смоделированными параметрами интенсивности вытеснения выполняют на по-ячейковой основе.
15. Способ по п. 13 или 14, в котором изменение разности между соответствующими наблюдаемыми и смоделированными параметрами расположения фронта содержит этап, на котором измеряют расстояние каждой ячейки до фронта.
16. Способ по п. 13, 14 или 15, в котором измеренную разность между соответствующими наблюдаемыми и смоделированными параметрами расположения фронта и измеренную разность между соответствующими наблюдаемыми и смоделированными параметрами интенсивности вытеснения комбинируют с относительным весом перед упомянутым вычислением несовпадения.
17. Способ по любому из пп. 12-16, в котором упомянутая методика адаптации по истории использует методику ансамблевого фильтра Кальмана.
18. Способ по любому из пп. 1-17, в котором преобразование выполняют для каждого рассматриваемого флюида.
19. Способ по п. 18, в котором преобразование выполняют для каждого флюида отдельно.
20. Способ по п. 18, в котором преобразование выполняют одновременно для каждого флюида, или любого подмножества упомянутых флюидов.
21. Способ по любому пп. 1-20, дополнительно содержащий этап, на котором используют результаты упомянутого способа для облегчения извлечения углеводородов из пласта-коллектора.
22. Компьютерная программа, содержащая компьютерно-читаемые инструкции, которые, при выполнении на подходящем компьютерном устройстве, заставляют компьютерное устройство выполнять способ по любому из пп. 1-21.
23. Компьютерный программный носитель, содержащий компьютерную программу по п. 22.
24. Устройство, выполненное с возможностью выполнения этапов способа по любому из пп. 1-21.