Scientific journal
Fundamental research
ISSN 1812-7339
"Перечень" ВАК
ИФ РИНЦ = 1,674

APPLICATION OF THE NONLINEAR LAW OF FILTRATION IN NATURAL POROUS RESERVOIRS IN HYDRODYNAMIC MODELS

Strekalov A.V. 1 Sarancha A.V. 2
1 Federal state budget higher educational institution «Tyumen State Oil and Gas University»
2 Federal state budget higher educational institution «Tyumen State Oil and Gas University»
This article considers the problem of mathematical physics arising in the theoretical description of the process the energy of a system of pore channels on the flow filtering them through the media. It is known that in conditions of the most primitive accounting of viscous forces within the pore friction, such processes are expressed in the law of filtration Darcy. However, its applicability is very limited range of speeds and pressure gradients, both from the critically low and high values. Thus, the applicability of Darcy’s law arise non-linear hydraulic effects associated, with capillary forces and rheological properties of environments, and with the output parameters of the flow outside the laminar flow conditions. The article demonstrates several techniques of mathematical description of the process inside porous viscous friction for use in the preparation of the system of equations of hydrodynamic models of oil and gas productive formations.
hydrodynamic model
nonlinear filtration law
1. Grachev S.I., Strekalov A.V. Opyt v reshenii zadach modelirovanija i optimizacii razrabotki mestorozhdenij nefti i gaza // Vestnik CKR Rosnedra. 2012. no. 2.
2. Svidetelstvo no. 2002611864 o registracii programmy dlja JeVM. Kompleks universalnogo modelirovanija tehnicheskih gidravlicheskih sistem podderzhanija plastovogo davlenija (HydraSym). M., 2002.
3. Semenov A.A., Dmitriev N.M, Kadet V.V., Mihajlov N.N. Jeffekt asimmetrii pri filtracii v anizotropnyh poristyh sredah // Tehnologii nefti i gaza. 2007. no. 1(48). рр. 52–55.
4. Semenov A.A. O variantah nelinejnyh zakonov filtracii dlja poristyh sred, projavljajushhih v zakone Darsi ortotropnye filtracionnye svojstva // Tezisy dokladov 6-j vserossijskoj konferencii molodyh uchenyh, specialistov i studentov po problemam gazovoj promyshlennosti Rossii. M.: Neft i gaz, 2005.
5. Strekalov A.V. Matematicheskie modeli gidravlicheskih sistem dlja upravlenija sistemami podderzhanija plastovogo davlenija. Tjumen: OAO Tjumenskij dom pechati, 2007. 664 р.

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

Для удобства дальнейшего описания и вследствие явного подобия закона Дарси закону Ома в случае соблюдения линейного закона фильтрации будем понимать под гидравлическим сопротивлением ячейки i пласта между определенной парой граней отношение

strekalov01.wmf [Па?с/м3], (1)

где S – ось, вдоль которой рассматривается течение в элементе; ?ps = pc – pgran – перепад давления между центром ячейки и одной из граней; qg – объемный расход флюида через грань g.

Таким образом, замыкающее отношение, соответствующее линейному закону течения Дарси вдоль оси S в ячейке, будет

strekalov02.wmf (2)

где S – ось – направление от центра ячейки к ее грани gran, 0 = X, 1 = Y, 2 = Z; Ls – длина ячейки вдоль оси S, м; Fgran – площадь грани (средняя площадь фильтрации), м2;

Сопротивление согласно (1) и (2):

strekalov03.wmf (3)

где k(s) – абсолютная проницаемость породы ячейки вдоль оси S, м2.

Направления фильтрации ТС будем условно полагать только вдоль осей X, Y и Z.

Наличие множества компонентов текучих сред при фильтрации учитывается согласно относительным проницаемостям. В этом случае сопротивление фильтрации компонента Ф вдоль оси S ячейки будет (рис. 1):

strekalov04.wmf (4)

где ?j(?j) – функция относительной фазовой проницаемости (ОФП) компонента j от его объемной насыщенности – ?j, д.е.; ?j – динамическая вязкость компонента j, Па?с; ?ji(?i) – функция относительного изменения ОФП компонента j от насыщенности компонента i, д.е.; N – количество рассматриваемых компонентов, шт.

pic_29.tif

Рис. 1. Схема элемента гидродинамической модели (ГДМ)

Функция fg(qg) отражает закон фильтрации, иначе говоря, закон вязкостного внутрипорового трения, который также может быть получен через функцию градиента давления от скорости фильтрации:

strekalov05.wmf (5)

strekalov06.wmf (6)

где ? = qg/Fg – скорость фильтрации через грань gran, м/с; qg – объемный расход среды, м3/с; ? – кинематическая вязкость, м2/с; kф – фазовая проницаемость, м2; ? – плотность флюида, кг/м3; ?0(?) – функция линейного закона внутрипорового трения; ? – объемная насыщенность, д.е.; m – текущая (так как она изменчива при изменении внутрипорового давления) пористость, д.е.

Описанные отношения (1)–(6) соответствуют линейному закону фильтрации.

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

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

1. Общий случай – нелинейные произвольные законы фильтрации. Нелинейный закон фильтрации (I приоритет) по функции корректировки градиента давления от числа Рейнольдса – ?I(?):

strekalov07.wmf (7)

где Rd(Re) – функция, зависящая от числа Рейнольдса для пористых сред и показывающая увеличение или уменьшение сопротивления при изменении Re. Функция Rd(Re) равна 1 д.е., если линейный закон фильтрации соблюдается.

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

Здесь strekalov08.wmf – число Рейнольдса, д.е., а strekalov09.wmf – истинная скорость фильтрации, м/с.

2. Общий случай – нелинейные произвольные законы фильтрации. Нелинейный закон фильтрации (II приоритет) учитывается в функции корректировки градиента давления в зависимости от истинной скорости ?II(?). Данная функция задается в исходных данных формулой или табулирована с последующей интерполяцией:

strekalov10.wmf (8)

где ?(?ист) – зависящая от истинной скорости фильтрации и показывает увеличение или уменьшение сопротивления при изменении ?ист. Функция равна единице, если закон соответствует закону Дарси.

3. Частный случай – нелинейный составной закон фильтрации (III приоритет) – ?III(?).

Для расчета зависимости градиента давления от скорости фильтрации или наоборот требуется разделить закон фильтрации ?III(?) на две области определения:

– первая соответствует линейному закону фильтрации

strekalov11.wmf

т.е. strekalov12.wmf (9)

– вторая соответствует нелинейному закону фильтрации

strekalov13.wmf

при strekalov14.wmf (10)

strekalov15.wmf

при strekalov16.wmf (11)

где ? – коэффициент, равный 1 Па/м; A и B – коэффициенты, отыскиваемые в процессе настройки ?III(?); ? – показатель нелинейности, д.е., который задается в исходных данных.

Переход от линейного к нелинейному закону фильтрации определяется числом Рейнольдса по Щелкачев:

strekalov17.wmf (12)

Критическая скорость и соответствующий ей критический градиент давления

strekalov18.wmf

strekalov19.wmf (13)

где Reкр – критическое число Рейнольдса (задается в исходных данных), д.е.

Коэффициенты A и B находятся решением системы уравнений:

strekalov20.wmf (14)

strekalov21.wmf (15)

На рис. 2 показаны примеры составного закона фильтрации.

pic_30.tif

Рис. 2. Составной закон фильтрации для различных коэффициентов j при a = 1

Градиент давления при известной скорости фильтрации для данного (составного) закона фильтрации можно определить решением уравнения

strekalov22.wmf (16)

strekalov23.wmf (17)

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

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

Неизвестный расход флюида – qg через грань – g можно рассчитать, решив уравнение вида

strekalov24.wmf (18)

при известном перепаде давления ?pg.

pic_31.tif

Рис. 3. Схема расчета перетока между ячейками

Учитывая возможность выбора четырех вариантов закона фильтрации из ?0(?), ?I(?), ?II(?) и ?III(?) в дальнейшем будем полагать закон вязкостного внутрипорового трения заданным для каждой ячейки в виде функций, связывающих градиент давления и скорости фильтрации ?i(?), где i – глобальный или локальный индекс ячейки. Если далее индекс не указан, то значит, имеется в виду общий случай.

pic_32.tif

Рис. 4. Общий вид зависимости закона фильтрации

Как видно из графика (рис. 4), strekalov25.wmf – критический градиент начала фильтрации, который складывается из градиента начала течения для флюида (задается в «PVT» свойствах флюида [2]) и из градиента начала фильтрации порового пространства пласта, который задается в свойствах каждого пласта.

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

strekalov26.wmf (19)

где ?(?) – закон фильтрации, который может быть одним из четырех приоритетов: линейный и нелинейные – I, II, и III.

Известно, что большинство фильтрационно-емкостных свойств (ФЕС) динамичны и изменяются в зависимости от внутрипорового давления.

Так, например, при изменении пористости, вызванном изменением внутрипорового давления, изменяется и абсолютная проницаемость, а следовательно, и все фазовые проницаемости (рис. 5). В общем виде ?m(?m) – это функция относительного приращения проницаемости от относительного приращения пористости.

pic_33.tif

Рис. 5. Зависимость относительного изменения (множителя) проницаемости от относительного приращения пористости

Здесь величина strekalov27.wmf – относительное приращение пористости, д.е.; ?k – относительное изменение проницаемости при относительном приращении пористости, д.е.

То есть итоговой абсолютной проницаемостью в любом из законов фильтрации будет

strekalov28.wmf (20)

Также возможно предусмотреть корректировку абсолютной проницаемости от внутрипорового давления напрямую через функцию ?p(?p) (рис. 6).

pic_34.tif

Рис. 6. Зависимость относительного изменения (множителя) проницаемости от относительного приращения давления

Здесь подобно предыдущему примеру strekalov29.wmf – относительное приращение давления, д.е.; p0 – начальное пластовое давление, которому соответствует начальная абсолютная проницаемость – k0; ?k – относительное изменение проницаемости при относительном приращении давления, д.е.

Данная функция необходима для учета изменения проницаемости от силы воздействия на пристеночную область – там, где ТС контактирует со стенками поровых каналов.

В этом случае итоговой абсолютной проницаемостью будет

strekalov30.wmf (21)

Полученная проницаемость в зависимости от выбранного закона фильтрации подставляется в соответствующее математическое описание (7)–(9).

Рассмотренное здесь математическое описание законов фильтрации интегрировано в ГДМ месторождений нефти и газа и используется в расчетном комплексе Hydra’Sym [2].

Результаты ряда лабораторных исследований авторов данной статьи и других ученых [3, 4] подтверждают необходимость учета в комплексной ГДМ месторождения нелинейных законов фильтрации, учитывающих уникальные особенности пористой среды в сочетании с физическими свойствами ТС.

В большей части случаев рассмотренного в статье математического описания достаточно для адаптации ГДМ посредством прямого переноса закона фильтрации из лабораторных исследований керна в модель через (7) и (8). Такие закономерности в целом сходны с видом функции на рис. 4 и отражают критический градиент давления начала течения, повышенные сопротивления на низких и высоких скоростях течения, линейный закон фильтрации на некоторых средних скоростях течения.

Рецензенты:

Грачев С.И., д.т.н., профессор, заведующий кафедрой «Разработка и эксплуатация нефтяных и газовых месторождений», Институт геологии и нефтегазодобычи, ФГБОУ «Тюменский государственный нефтегазовый университет», г. Тюмень;

Сохошко С.К., д.т.н., заведующий кафедрой «Моделирование и управление процессами нефтегазодобычи», Институт геологии и нефтегазодобычи, ФГБОУ «Тюменский государственный нефтегазовый университет», г. Тюмень.