RU2670174C1 - Способ многомасштабного моделирования нелинейных процессов подземной гидродинамики - Google Patents
Способ многомасштабного моделирования нелинейных процессов подземной гидродинамики Download PDFInfo
- Publication number
- RU2670174C1 RU2670174C1 RU2017144198A RU2017144198A RU2670174C1 RU 2670174 C1 RU2670174 C1 RU 2670174C1 RU 2017144198 A RU2017144198 A RU 2017144198A RU 2017144198 A RU2017144198 A RU 2017144198A RU 2670174 C1 RU2670174 C1 RU 2670174C1
- Authority
- RU
- Russia
- Prior art keywords
- phase
- physical properties
- phases
- imb
- equations
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 40
- 230000008569 process Effects 0.000 title claims abstract description 16
- 238000006073 displacement reaction Methods 0.000 claims abstract description 37
- 230000000704 physical effect Effects 0.000 claims abstract description 24
- 238000004364 calculation method Methods 0.000 claims abstract description 23
- 229930195733 hydrocarbon Natural products 0.000 claims abstract description 10
- 150000002430 hydrocarbons Chemical class 0.000 claims abstract description 10
- 230000000694 effects Effects 0.000 claims abstract description 9
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 7
- 239000004215 Carbon black (E152) Substances 0.000 claims abstract description 5
- 238000001914 filtration Methods 0.000 claims description 15
- 230000035699 permeability Effects 0.000 claims description 10
- 230000037230 mobility Effects 0.000 claims description 7
- 239000013598 vector Substances 0.000 claims description 5
- 239000006185 dispersion Substances 0.000 claims description 4
- 230000010349 pulsation Effects 0.000 claims description 3
- 239000000126 substance Substances 0.000 abstract 1
- 239000012071 phase Substances 0.000 description 39
- 239000007788 liquid Substances 0.000 description 14
- 238000012935 Averaging Methods 0.000 description 8
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 7
- 238000011161 development Methods 0.000 description 5
- 238000005755 formation reaction Methods 0.000 description 5
- 238000010276 construction Methods 0.000 description 4
- 239000012530 fluid Substances 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000004422 calculation algorithm Methods 0.000 description 3
- 230000007423 decrease Effects 0.000 description 3
- 238000009792 diffusion process Methods 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 230000014509 gene expression Effects 0.000 description 3
- 239000000203 mixture Substances 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 230000002028 premature Effects 0.000 description 2
- 238000011084 recovery Methods 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000000739 chaotic effect Effects 0.000 description 1
- 239000003795 chemical substances by application Substances 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000002708 enhancing effect Effects 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- HJUFTIJOISQSKQ-UHFFFAOYSA-N fenoxycarb Chemical compound C1=CC(OCCNC(=O)OCC)=CC=C1OC1=CC=CC=C1 HJUFTIJOISQSKQ-UHFFFAOYSA-N 0.000 description 1
- 230000004907 flux Effects 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 239000007791 liquid phase Substances 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 230000000149 penetrating effect Effects 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06G—ANALOGUE COMPUTERS
- G06G7/00—Devices in which the computing operation is performed by varying electric or magnetic quantities
- G06G7/48—Analogue computers for specific processes, systems or devices, e.g. simulators
Landscapes
- Engineering & Computer Science (AREA)
- Computer Hardware Design (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- General Physics & Mathematics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Изобретение относится к способам многомасштабного моделирования нелинейных процессов подземной гидродинамики. Сущность: разбивают исследуемую геологическую структуру на крупномасштабные блоки (КМБ). Разбивают КМБ на более мелкомасштабные блоки (ММБ). Вычисляют значения физических свойств для ММБ. Вычисляют значения физических свойств для КМБ с учетом дополнительных корректирующих членов уравнений, коэффициенты для которых получают из вычислений для ММБ. При этом в качестве физического свойства, определяемого на ММБ, определяют неустойчивость на фронте вытеснения углеводорода из пласта. Технический результат: повышение точности расчетов и сокращение времени их проведения. 4 з.п. ф-лы, 6 ил.
Description
Изобретение относится к области информационных технологий, предназначенных для специализированной обработки данных, в частности к способу вычислительного моделирования процессов в подземной гидродинамике и к разработке аппаратно-программных средств предсказательного моделирования в подземной гидродинамике.
Существуют коммерческие гидродинамические симуляторы, не включающие в себя возможность учета мелкомасштабной (подсеточной) неустойчивости:
1) ECLIPSE, https://www.software.slb.com/products/eclipse.
2) STARS CMG, http://www.cmgl.ca/stars.
3) Техсхема, http://tecscheme.org/.
4) Т-навигатор, http://rfdyn.ru/technology/.
Программного пакета, позволяющего проводить моделирование нелинейных процессов вытеснения углеводородов с учетом возникающей неустойчивости фронта на мелком масштабе, на данный момент нет.
Методы математического моделирования процессов фильтрации на основе параметров, полученных из геологических исследований в скважинах и лабораторных исследований свойств жидкостей и породы, описываются во многих патентах (RU 2148169 С1, 17.09.1998; RU 2166619 С1, 03.09.1999; RU 2166630 С1, 03.09.1999; RU 2183268 С2, 14.07.2000; RU 2522494 С1, 10.05.2013).
Одним из недостатков известных методов является то, что моделирование происходит на макроуровне без учета возможных процессов на микроуровне, а именно, неустойчивости, возникающей на фронте вытеснения углеводорода из пласта, что оказывает существенное влияние на результат моделирования. Впоследствии это может привести к выбору неправильной стратегии разработки месторождения и к существенным финансовым потерям организации, использующей представленные методы моделирования.
Учет физических свойств на микроуровне производится в методе, описанном в патенте RU 2428739 С2, 16.05.2007. Суть метода заключается в разбиении исследуемой геологической структуры на крупномасштабные блоки (КМБ), а затем разбиение КМБ на более мелкомасштабные блоки (ММБ) и вычисление значения физических свойств для ММБ после чего путем ремасштабирования осуществляют вычисление значения физических свойств для КМБ.
Недостатком известного метода является то, что неустойчивость, возникающая на фронте вытеснения, в данной модели не учитывается как на макроуровне, так и на микроуровне, из-за чего снижается точность расчетов и увеличивается время их проведения при построении реальной модели прогнозирования протекания процессов в нефтесодержащих пластах.
Задачей является повышение точности расчетов и сокращение времени их проведения при построении реальной модели прогнозирования протекания процессов в нефтесодержащих пластах.
Технический результат заключается в изменении модели на макроуровне за счет дополнительных потоков в уравнениях в зависимости от физических свойств, полученных при исследовании процесса на микроуровне.
Указанные задача и технический результат достигаются использованием способа многомасштабного моделирования нелинейных процессов подземной гидродинамики, включающим разбиение исследуемой геологической структуры на крупномасштабные блоки (КМБ), разбиение КМБ на более мелкомасштабные блоки (ММБ), вычисление значения физических свойств для ММБ после чего вычисление значения физических свойств для КМБ, при этом, вычисления значений физических свойств для КМБ осуществляют с учетом дополнительных корректирующих членов уравнений, коэффициенты для которых получены из вычислений для ММБ, а в качестве физического свойства, определяемого на ММБ, определяют неустойчивость на фронте вытеснения углеводорода из пласта.
А также тем, что для вычисления физических свойств для ММБ формируют диапазон возможных безразмерных параметров и контрольных точек в фазовом пространстве этих параметров внутри этого диапазона, проводят различные варианты расчетов на основе указанных параметров и формируют интерполяционные зависимости результатов от введенных параметров, с помощью которых определяют коэффициенты дополнительных корректирующих членов уравнений для определения значений физических свойств для КМБ.
А также тем, что вычисление физических свойств для ММБ осуществляют из соотношений:
где ρ - плотность, φ - пористость, s1, s2 - насыщенности фаз, - относительные проницаемости фаз, μ1, μ2 - динамические вязкости фаз, K - абсолютная проницаемость, v1, v2 - вектора скорости фильтрации фаз, р1, р2 - давление в фазах, рC - капиллярное давление.
А также тем, что вычисление физических свойств для КМБ осуществляют из соотношений:
j=1, 2, где m1, m2 - относительные мобильности фаз, v - суммарная объемная скорость фильтрации, ρ1, ρ2 - плотности фаз, φ - пористость, s1, s2 - насыщенности фаз, , - относительные проницаемости фаз, μ1, μ2 - динамические вязкости фаз, v1, v2 - вектора скорости фильтрации фаз; штрихом обозначены пульсационные составляющие переменных, чертой сверху обозначены средние величины.
А также тем, что дополнительные корректирующие члены уравнений определяют из соотношений:
где М - отношение вязкостей, Ре - число Пекле, Δ - характерный размер расчетной ячейки глобальной задачи, d - характерная дисперсия фаз, возникающая за счет капиллярных эффектов, m1, m2 - относительные мобильности фаз, K - абсолютная проницаемость, v - суммарная объемная скорость фильтрации, рC - капиллярное давление, μ1, μ2 - динамические вязкости фаз, s1, s2 - насыщенности фаз.
Все существующие методы объединяет то, что для моделирования на макро- и микро- уровнях используются одинаковые системы уравнений. Значительным отличием изобретения от известных существующих методов является то, что в патентуемом способе для проведения апскейлинга в дифференциальные уравнения на макроуровне добавляются дополнительные корректирующие члены, коэффициенты для которых получены из исследований на микроуровне. В зависимости от физических свойств, полученных при исследовании процесса на микроуровне, меняется модель на макроуровне за счет учета дополнительных потоков в уравнениях. Такой подход дает возможность учесть неустойчивость, не прибегая к мелкому разбиению, что также сопровождается экономией расчетного времени.
Коэффициент извлечения нефти по большинству отечественных месторождений не превосходит 30%, а в большом числе случаев опускается ниже 20%, при этом, трудно извлекаемые запасы нефти в России составляют 60-70%. Один из методов нефтедобычи предполагает вытеснение углеводорода, при этом, на фронте вытеснения развивается неустойчивость, существенно сказывающаяся на качестве нефтедобычи.
Изобретение позволяет создать специализированный программный комплекс для предсказательного моделирования неустойчивого вытеснения нефти из природных коллекторов на малом масштабе, работающего как самостоятельно, так и допускающего интегрирование в любой коммерческий симулятор подземной гидродинамики.
Особенностью специализированного программного комплекса является то, что он позволит осуществлять многомасштабное моделирование нелинейных процессов вытеснения углеводородов с учетом возникающей неустойчивости фронта и даст возможность оценки качества вытеснения в условиях неустойчивых режимов с разрешением на уровне 1 см, и в то же время без уменьшения общей области моделирования, а также без увеличения времени вычислений. Полученные результаты могут привести к увеличению нефтеотдачи на 3-10%.
Патентуемый способ позволит учитывать эффекты с разрешением на несколько порядков выше, но без увеличения времени расчета.
Фиг. 1 иллюстрирует двумерную картину развития неустойчивости фронта вытеснения жидкости при перепаде вязкостей 1:100 (вытесняющей к вытесняемой жидкости). Изображена насыщенность вытесняющей фазы (минимальная - синий, максимальная - пурпурный). Под каждой картинкой подписано безразмерное время; единице в этой шкале соответствует время, за которое бы при устойчивом фронте вытеснения вся жидкость в области была бы замещена. Видно, что при неустойчивом фронте прорыв вытесняющей жидкости осуществляется в 5 раз раньше этого момента, но объем вытеснения значительно уменьшен.
Фиг. 2 иллюстрирует вид поверхности, разделяющей вытесняемую и вытесняющую жидкости, с двух точек зрения в последовательные моменты времени при трехмерном моделировании. Соотношение вязкостей вытесняющей и вытесняемой жидкости 1:100. Под каждой картинкой подписано безразмерное время; единице в этой шкале соответствует время, за которое бы при устойчивом фронте вытеснения вся жидкость в области была бы замещена. Прорыв происходит раньше, чем при двумерном моделировании, результаты которого показаны на предыдущей фигуре.
Фиг. 3 иллюстрирует общий вид уравнений глобальной модели вытеснения. Обозначено: ρ - плотность, φ - пористость, s1, s2 - насыщености фаз, - относительные проницаемости фаз, μ1, μ2 - динамические вязкости фаз, K - абсолютная проницаемость, мобильности фаз, v1, v2 - вектора скорости фильтрации фаз, р1, р2 - давление в фазах, рС - капиллярное давление.
Фиг. 4 иллюстрирует дополнительный член, возникающий при осреднении модели вытеснения в уравнении динамики массы фаз, и его структура. Обозначено: mj - относительная мобильность фазы, v - суммарная объемная скорость фильтрации, штрихом обозначены пульсационные составляющие переменных, чертой сверху обозначены средние величины.
Фиг. 5 иллюстрирует моделирование дополнительного потока фазы, возникающего за счет неустойчивости вытеснения на микроуровне. Обозначено: М - отношение вязкостей, Ре - число Пекле, Δ - характерный размер расчетной ячейки глобальной задачи, d - характерная дисперсия фаз, возникающая за счет капиллярных эффектов.
Фиг. 6 иллюстрирует блок-схему работы алгоритма, с включением расчета на двух уровнях, глобальном и локальном.
Один из основных методов нефтедобычи предполагает вытеснение углеводорода менее вязкой жидкостью (водой), при этом, на фронте вытеснения развивается неустойчивость, существенно сказывающаяся на качестве нефтедобычи.
Задача вытеснения углеводородов осложняется тем, что в практических случаях ее приходится решать на значительном пространственном масштабе. В то же время из-за использования менее вязкой воды для вытеснения более вязкой смеси углеводородов (пластовой нефти) этот процесс является неустойчивым: фронт вытеснения, даже если он был вначале прямым, искажается за счет пальцеобразных выступов, по которым менее вязкая вода прорывается к добывающей скважине, оставляя позади себя невытесненную нефть. После достижение одним из таких выступов добывающей скважины картина вытеснения резко меняется - большая часть подаваемой воды будет направляться в канал, созданный лидирующим пальцем. В скважину начинает закачиваться уже не нефть, а ее смесь с вытесняющей водой, причем процент воды будет быстро увеличиваться. Оставшаяся же в пласте нефть после такого прорыва будет практически неподвижна. На фиг. 1 приведена картина развития неустойчивости фронта вытеснения жидкости при перепаде вязкостей 1:100; маловязкая жидкость вытесняет более вязкую, видно образование каналов, по которым в основном и течет вытесняющая жидкость, оставляя позади вытесняемую.
Такая неустойчивость многомасштабна, и дело ее моделирования осложняется тем, что наиболее быстро растут мелкомасштабные «пальцы». Основная проблема прямого моделирования неустойчивости в том, что разрешение мелкомасштабных потоков на практической сетке нереальна даже с использованием современной вычислительной техники. Например, при обработке малого участка месторождения размером 10×10×1 км с ячейками размером 1×1×1 м потребуется 100 млрд. ячеек. Такая крупная задача находится на пределе возможности современных вычислительных комплексов. Но разрешить неустойчивости размером по ширине менее 3 ячеек численная схема не способна. Тем самым пальцеобразные выступы вытесняющего агента шириной менее 3 м, проникающие в вытесняемую нефть, на такой сетке не только не будут обнаружены, но решение вообще не будет учитывать их влияние, хотя бы и скрытое, на процесс вытеснения. Уменьшение же линейного размера ячейки в 10 раз увеличивает число ячеек двумерной сетки в 100, а трехмерной в 1000 раз. В то же время опыт показывает, что значение имеет пальцеобразная неустойчивость фронта вытеснения для размеров порядка 1 см. Весьма важно и то, что реальный фронт вытеснения имеет существенно трехмерную структуру, что видно из развития трехмерной картины вытеснения, изображенной на фиг. 2, где в двух проекциях изображен вид поверхности, разделяющей вытесняемую и вытесняющую жидкости. Картина эта показывает, что прорыв маловязкой жидкости в случае неустойчивости фронта вытеснения проходит по каналам, имеющим форму, близкую к цилиндрической, и тем самым объем невытесненной жидкости при трехмерном моделировании выше, чем при двумерном. Это говорит о том, что для решения задачи, позволяющей учитывать мелкую структуру вытеснения, требуется именно трехмерное моделирование с мелким разрешением. Измельчение исходной сетки до такого масштаба делает решение задачи вытеснения если и не невозможным, то не применимым на практике.
В то же время учесть влияние мелкомасштабной (подсеточной) неустойчивости при решении глобальной задачи на грубой сетке возможно. Для этого с уравнениями глобальной задачи поступают так же, как поступают с уравнениями Навье-Стокса при моделировании эффектов подсеточной турбулентности: вводится временной масштаб осреднения параметров (в данном случае насыщенности фаз, их скорости и давления). Затем решение представляется как сумма среднего на этом временном масштабе и осциллирующего, после чего уравнения осредняются.
В случае уравнений Навье-Стокса такое осреднение называется осреднением по Рейнольдсу, оно приводит к тому, что в осредненных уравнениях появляются новые члены, прежде всего в конвективных составляющих уравнений импульса. Эти новые члены, дополнительные потоки импульса, требуют дополнительного эмпирического моделирования, которое в различных моделях турбулентности проводится различным способом, но самым простым и, тем не менее, результативным во многих случаях является способ алгебраический.
Уравнения динамики вытеснения приведены на фиг. 3. Система содержит уравнение динамики массы каждой фазы (3.1), закон Дарси для каждой фазы (3.2), связь между давлениями в фазах через капиллярное давление (3.3) и дополнительные алгебраические соотношения (3.4). Уравнения глобальной системы приведены в размерном виде.
В случае подобной обработки уравнений динамики насыщенности новые члены также возникнут, и они будут соответствовать дополнительным потокам массы фаз. Эти потоки учитывают влияние мелкомасштабной неустойчивости, которая в данной модели скрыта на фоне среднего решения, и, тем не менее, выражается в резком усилении фактических фильтрационных потоков фаз. Собственно, именно эти потоки и ответственны за возможный преждевременный прорыв воды к добывающей скважине; преждевременным он будет относительно времени прорыва, полученного на грубой (глобальной) сетке без учета дополнительных членов. Возникающий в результате осреднения член содержится в уравнениях динамики массы потоков, и его структура приведена на фиг. 4. Уравнение (4.1) - это уравнение динамики массы фаз после осреднения, содержащее дополнительный член, выражение (4.2) определяет структуру зависимости потока от насыщенности, которая подвергается дополнительному осреднению, выражения (4.3) и (4.4) определяют относительную мобильность фаз и суммарную скорость фильтрации, соответственно.
Данные потоки моделируются алгебраически, будучи частично разделены на конвективную и диффузионную составляющие. Дополнительные потоковые члены должны обнуляться в том случае, когда вытеснение устойчиво на микроуровне, и тогда, когда насыщенность какой-либо фазы уменьшается до остаточной. Они зависят алгебраически от небольшого набора безразмерных параметров, регулирующих неустойчивость вытеснения. В частности, члены, усиливающие неустойчивость, должны в общем случае увеличиваться с ростом отношения вязкости вытесняемой жидкости к вязкости вытесняющей, и уменьшаться с ростом диффузии, возникающей вследствие капиллярных эффектов, и формирующей вместе с характерной скоростью фильтрации жидких фаз безразмерное число Пекле. На фиг. 5 приведена форма моделирования дополнительных потоков (5.1) и перечислены основные безразмерные параметры (5.2), от которых в большей или меньшей степени зависят дополнительные фильтрационные потоки фаз. Выражение (5.3) показывает, каким образом оценивается характерная диффузия, возникающая за счет капиллярных эффектов; черта сверху означает осреднение этих данных по глобальной ячейке. Форма самой зависимости дополнительного скаляра F в конвективном члене и тензора D в дисперсионном члене, моделирующих дополнительны поток, от указанных параметров, определяется характером неустойчивости и подлежит нахождению путем решения локальной задачи.
Патентуемый способ предполагает создание программного комплекса, предназначенного прежде всего для определения параметров алгебраической модели эффектов мелкомасштабной неустойчивости. Это достигается путем решения двух задач - локальной (мелкомасштабной) и глобальной (задачи с крупными расчетными ячейками). Алгоритм состоит из нескольких блоков: локальное построение сетки, локальное построение начального состояния, локальный расчет (последовательность шагов по времени), постпроцессор локальных расчетов (собственно формирование потоков, моделирующих неустойчивость на микроуровне); глобального построения сетки, глобального построения начального состояния, глобального расчета (шагов по времени). Дополнительный блок - постпроцессор, обрабатывающий данные глобального моделирования. Все расчетные блоки можно подразделить на две крупных группы: локальную (расчет на микроуровне и получение дополнительных потоков) и глобальную (расчет на макроуровне с применением членов, моделирующих особенности течения на микроуровне).
Последовательность работы алгоритма, в том случае, когда решение локальной задачи интегрировано в глобальную, изображена в виде блок-схемы на фиг. 6.
Подробности работы таковы.
Вначале строится глобальная сетка под цели расчета крупномасштабной задачи вытеснения.
Затем для работы локального блока строится сетка, моделирующая типичную ячейку глобальной сетки как рабочую область.
Формируется диапазон возможных безразмерных параметров локальной задачи и контрольные точки в фазовом пространстве этих параметров внутри этого диапазона. Параллельно на локальной сетке проводятся расчеты для различных вариантов (строятся начальные условия, проводятся шаги по времени, формируется результат в виде профиля насыщенности и потоков компонент на входе и выходе). Для этого используются модуль построения начального состояния локальной задачи и модуль ее решения.
Затем модуль постпроцессора локальной задачи обрабатывает полученные результаты, формируя интерполяционную зависимость результатов локальной задачи от введенных параметров. В зависимости от полученных данных такая зависимость представляется либо простым алгебраическим соотношением, в котором следует лишь определить коэффициенты, либо (в общем случае) табличными данными, которые могут быть интерполированы локально различными методами: сплайн с привязкой к точным значениям (для регулярного расположения полученных данных), сплайн без таковой привязки с минимальным средне-квадратичным отклонением (для менее регулярного их расположения), использование функций радиального базиса (для случая хаотично или неравномерно расположенных данных в фазовом пространстве большой размерности).
Затем полученные данные используются при решении глобальной задачи, в модулях которой - построения начального состояния и решения - вводятся дополнительные потоки, сформированные моделированием неустойчивости на микроуровне.
В качестве решателя глобальной задачи может быть использован как собственный решатель, так и любой существующий решатель для задач подземной гидродинамики.
Применение изобретения позволит:
1. Создать программное обеспечение для предсказательного моделирования неустойчивого вытеснения нефти из природных коллекторов;
2. Провести компьютерный анализ и автоматизированную обработку результатов численных экспериментов;
3. Разработать функциональные формы и таблицы коэффициентов для описания динамики эволюции зоны взаимного проникновения жидкостей и ее структуры;
4. Создать приложение, совместимое с коммерческими гидродинамическими симуляторами и улучшающего качество моделирования вытеснения;
Claims (11)
1. Способ многомасштабного моделирования нелинейных процессов подземной гидродинамики, включающий разбиение исследуемой геологической структуры на крупномасштабные блоки (КМБ), разбиение КМБ на более мелкомасштабные блоки (ММБ), вычисление значения физических свойств для ММБ, последующее вычисление значения физических свойств для КМБ, отличающийся тем, что вычисления значений физических свойств для КМБ осуществляют с учетом дополнительных корректирующих членов уравнений, коэффициенты для которых получены из вычислений для ММБ, а в качестве физического свойства, определяемого на ММБ, определяют неустойчивость на фронте вытеснения углеводорода из пласта.
2. Способ многомасштабного моделирования по п.1, отличающийся тем, что для вычисления физических свойств для ММБ формируют диапазон возможных безразмерных параметров и контрольных точек в фазовом пространстве этих параметров внутри этого диапазона, проводят различные варианты расчетов на основе указанных параметров и формируют интерполяционные зависимости результатов от введенных параметров, с помощью которых определяют коэффициенты дополнительных корректирующих членов уравнений для определения значений физических свойств для КМБ.
3. Способ многомасштабного моделирования по п.1, отличающийся тем, что вычисление физических свойств для ММБ осуществляется из соотношений
4. Способ многомасштабного моделирования по п.1, отличающийся тем, что вычисление физических свойств для КМБ осуществляется из соотношений
, , , j=1, 2, где m1, m2 - относительные мобильности фаз, v - суммарная объемная скорость фильтрации, ρ1, ρ2 - плотности фаз, φ - пористость, s1, s2 - насыщенности фаз, , - относительные проницаемости фаз, μ1, μ2 - динамические вязкости фаз, v1, v2 - векторы скорости фильтрации фаз; штрихом обозначены пульсационные составляющие переменных, чертой сверху обозначены средние величины.
5. Способ многомасштабного моделирования по п.1, отличающийся тем, что дополнительные корректирующие члены уравнений определяются из соотношений
где М - отношение вязкостей, Ре - число Пекле, Δ - характерный размер расчетной ячейки глобальной задачи, d - характерная дисперсия фаз, возникающая за счет капиллярных эффектов, m1, m2 - относительные мобильности фаз, K - абсолютная проницаемость, v - суммарная объемная скорость фильтрации, рC - капиллярное давление, μ1, μ2 - динамические вязкости фаз, s1, s2 - насыщенности фаз.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2017144198A RU2670174C1 (ru) | 2017-12-18 | 2017-12-18 | Способ многомасштабного моделирования нелинейных процессов подземной гидродинамики |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
RU2017144198A RU2670174C1 (ru) | 2017-12-18 | 2017-12-18 | Способ многомасштабного моделирования нелинейных процессов подземной гидродинамики |
Publications (1)
Publication Number | Publication Date |
---|---|
RU2670174C1 true RU2670174C1 (ru) | 2018-10-18 |
Family
ID=63862460
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
RU2017144198A RU2670174C1 (ru) | 2017-12-18 | 2017-12-18 | Способ многомасштабного моделирования нелинейных процессов подземной гидродинамики |
Country Status (1)
Country | Link |
---|---|
RU (1) | RU2670174C1 (ru) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2703359C1 (ru) * | 2018-12-13 | 2019-10-16 | Общество с ограниченной ответственностью (ООО) "ЛУКОЙЛ-ПЕРМЬ" | Инженерный симулятор процесса добычи и транспортировки продукции скважин |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2428739C2 (ru) * | 2006-07-07 | 2011-09-10 | Эксонмобил Апстрим Рисерч Компани | Укрупнение сетки для моделей коллекторов путем повторного использования расчетов потока, полученных на основе геологических моделей |
WO2015153506A1 (en) * | 2014-03-31 | 2015-10-08 | Ingrain, Inc. | Digital rock physics-based trend determination and usage for upscaling |
RU2573739C2 (ru) * | 2011-02-28 | 2016-01-27 | Шлюмбергер Текнолоджи Б.В. | Многомасштабное цифровое моделирование породы для моделирования пласта |
-
2017
- 2017-12-18 RU RU2017144198A patent/RU2670174C1/ru active
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2428739C2 (ru) * | 2006-07-07 | 2011-09-10 | Эксонмобил Апстрим Рисерч Компани | Укрупнение сетки для моделей коллекторов путем повторного использования расчетов потока, полученных на основе геологических моделей |
RU2573739C2 (ru) * | 2011-02-28 | 2016-01-27 | Шлюмбергер Текнолоджи Б.В. | Многомасштабное цифровое моделирование породы для моделирования пласта |
WO2015153506A1 (en) * | 2014-03-31 | 2015-10-08 | Ingrain, Inc. | Digital rock physics-based trend determination and usage for upscaling |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
RU2703359C1 (ru) * | 2018-12-13 | 2019-10-16 | Общество с ограниченной ответственностью (ООО) "ЛУКОЙЛ-ПЕРМЬ" | Инженерный симулятор процесса добычи и транспортировки продукции скважин |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Geiger et al. | A novel multirate dual-porosity model for improved simulation of fractured and multiporosity reservoirs | |
De Rosis et al. | A comparison between the interpolated bounce-back scheme and the immersed boundary method to treat solid boundary conditions for laminar flows in the lattice Boltzmann framework | |
Langlo et al. | Macrodispersion for two-phase, immiscible flow in porous media | |
US20160341850A1 (en) | Geomechanical and geophysical computational model for oil and gas stimulation and production | |
Wang et al. | Discrete fracture modeling using Centroidal Voronoi grid for simulation of shale gas plays with coupled nonlinear physics | |
AU2015208052A1 (en) | Method for transient quasi three-dimensional simulation of multiphase fluid flow in pipelines | |
WO2012071090A1 (en) | Variable discretization method for flow simulation on complex geological models | |
NO334387B1 (no) | Multiskala, endelig volum-metode for bruk ved strømningssimulering under jordoverflaten | |
CN102224502A (zh) | 用于多孔介质中的流动的迭代多尺度方法 | |
US9418184B2 (en) | Determining flow through a fracture junction in a complex fracture network | |
EP3096252A2 (en) | Adaptive multiscale multi-fidelity reservoir simulation | |
WO2015048760A1 (en) | Interface point method modeling of the steam-assisted gravity drainage production of oil | |
Ravshanov et al. | An approximate analytical solution of the problem of fluid filtration in the multilayer porous medium | |
Gimenez et al. | An extended validation of the last generation of particle finite element method for free surface flows | |
Weijermars et al. | High-resolution streamline models of flow in fractured porous media using discrete fractures: Implications for upscaling of permeability anisotropy | |
Hamon et al. | Analysis of hybrid upwinding for fully-implicit simulation of three-phase flow with gravity | |
Chen et al. | Lattice Boltzmann simulation of particle motion in binary immiscible fluids | |
Douillet-Grellier et al. | Simulations of intermittent two-phase flows in pipes using smoothed particle hydrodynamics | |
RU2670174C1 (ru) | Способ многомасштабного моделирования нелинейных процессов подземной гидродинамики | |
US10699034B2 (en) | Flow transition technology | |
Mazo et al. | Petroleum reservoir simulation using super element method | |
Brito-Parada et al. | A finite element formulation to model the flow of flotation foams | |
Van Doren et al. | Controllability and observability in two-phase porous media flow | |
Matthäi et al. | A high‐order TVD transport method for hybrid meshes on complex geological geometry | |
Imankulov et al. | Design of HPC system for analysis the gel-polymer flooding of oil fields |