Что такое перколяционный способ

Перколяция

В физике и химии явлением перколяции (от лат. percōlāre , просачиваться, протекать) называется явление протекания или не протекания жидкостей через пористые материалы, электричества через смесь проводящих и непроводящих частиц и другие подобные процессы. Теория перколяции находит применение в описании разнообразных систем и явлений, в том числе таких, как распространение эпидемий и надежность компьютерных сетей.

Некоторые примеры задач, которые решаются через теорию перколяции:

• Сколько надо добавить медных опилок в ящик с песком, чтобы смесь начала проводить ток?
• Какой процент людей должен быть восприимчив к болезни, чтобы стала возможна эпидемия?

Содержание

Описание

Явление перколяции (или протекания среды) определяется:

1. Средой, в которой наблюдается это явление;
2. Внешним источником, который обеспечивает протекание в этой среде;
3. Способом протекания среды, который зависит от внешнего источника.

Пример

В качестве простейшего примера можно рассмотреть модель протекания (например электрического пробоя) в двумерной квадратной решетке, состоящей из узлов, которые могут быть проводящими или непроводящими. В начальный момент времени все узлы сетки являются непроводящими. Со временем источник заменяет непроводящие узлы на проводящие, и число проводящих узлов постепенно растет. При этом узлы замещаются случайным образом, то есть выбор любого из узлов для замещения является равновероятным для всей поверхности решетки.

Перколяцией называют момент появления такого состояния решетки, при котором существует хотя бы один непрерывный путь через соседние проводящие узлы от одного до противоположного края. Очевидно, что с ростом числа проводящих узлов, этот момент наступит раньше, чем вся поверхность решетки будет состоять исключительно из проводящих узлов.

Обозначим непроводящее и проводящее состояние узлов нулями и единицами соответственно. В двумерном случае среде будет соответствовать бинарная матрица. Последовательность замены нулей матрицы на единицы будет соответствовать источнику протекания.

В начальный момент времени матрица состоит полностью из непроводящих элементов:

0	0	0	0
0	0	0	0
0	0	0	0
0	0	0	0

При воздействии внешнего источника в матрице начинают добавляться проводящие элементы, однако поначалу их недостаточно для перколяции:

0	0	0	1
1	0	0	0
0	0	1	0
0	0	1	0

По мере увеличения числа проводящих узлов наступает такой критический момент, когда происходит перколяция, как показано ниже:

0	0	0	1
1	1	0	0
0	1	1	0
0	0	1	1

Видно, что от левой к правой границе последней матрицы имеется цепочка элементов, которая обеспечивает протекание тока по проводящим узлам (единицам), непрерывно следующим друг за другом.

Перколяция может наблюдаться как в решетках, так и других геометрических конструкциях, в том числе непрерывных, состоящих из большого числа подобных элементов или непрерывных областей соответственно, которые могут находиться в одном из двух состояний. Соответствующие математические модели называются решеточными или континуальными.

В качестве примера перколяции в непрерывной среде может выступать прохождение жидкости через объемный пористый образец (например, воды через губку из пеноообразующего материала), в котором происходит постепенное надувание пузырьков до тех пор, пока их размеров не станет достаточно для просачивания жидкости от одного края образца до другого.

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

Порог протекания

Совокупность элементов, по которым происходит протекание, называется перколяционным кластером. Будучи по своей природе связным случайным графом, в зависимости от конкретной реализации он может иметь различную форму. Поэтому принято характеризовать его общий размер. Порогом протекания называется минимальная концентрация, при которой возникает протекание.

Ввиду случайного характера переключений состояний элементов среды, в конечной системе чётко определенного порога (размера критического кластера) не существует, а имеется так называемая критическая область значений , в которую попадают значения порога перколяции, полученные в результате различных случайных реализаций. С увеличением размеров системы область сужается в точку.

Источник

Перколяция — Percolation

В физике , химии и науки о материалах , перколяции (от латинского percolare , «к фильтру» или «струйкой через») относится к движению и фильтрации жидкостей через пористые материалы. Это описывается законом Дарси . С тех пор были разработаны более широкие приложения, которые охватывают связность многих систем, смоделированных в виде решеток или графов, аналогично связности компонентов решетки в задаче фильтрации, которая модулирует способность к просачиванию.

СОДЕРЖАНИЕ

В течение последних десятилетий теория перколяции , математическое исследование перколяции, принесла новое понимание и новые методы в широкий круг вопросов в физике, материаловедении, сложных сетях , эпидемиологии и других областях. Например, в геологии под перколяцией понимается фильтрация воды через почву и проницаемые породы. Вода поступает перезарядки в грунтовые воды в грунтовых вод и водоносных слоев . В местах, где инфильтрационные бассейны или септические дренажные поля планируются для удаления значительного количества воды, необходимо заранее провести испытание на просачивание, чтобы определить, будет ли предполагаемая конструкция успешной или нет. В двумерной квадратной решетке перколяция определяется следующим образом. Сайт «занят» с вероятностью p или «пуст» (в этом случае его края удаляются) с вероятностью 1 — p; соответствующая проблема называется перколяцией сайтов, см. рис. 2.

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

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

Источник

Перколяция

Из Википедии — свободной энциклопедии

В физике и химии явлением перколяции (от лат. percōlāre — просачиваться, протекать) называется явление протекания или непротекания жидкостей через пористые материалы. Теория перколяции находит применение в описании разнообразных систем и явлений, в том числе таких, как распространение эпидемий, надежность компьютерных сетей и протекание электричества через смесь проводящих и непроводящих частиц.

Некоторые примеры задач, которые решаются через теорию перколяции:

• Сколько надо добавить медных опилок в ящик с песком, чтобы смесь начала проводить ток?
• Какая доля людей должна быть восприимчива к болезни, чтобы стала возможна эпидемия?

Изначально термин перколяция возник в гидрологии и обозначал воду, просачивающуюся через почву. Просачивание происходит только при превышении критической емкости, то есть максимальной влагоемкости почвы. Если это происходит по всей толщине ненасыщенной зоны, образуются новые грунтовые воды.

В фармацевтической технологии перколяция используется для получения активных ингредиентов из растений. Здесь обычно теплый растворитель, такой как вода или спирт, пропускается через растения или части растений; в лабораторной практике в основном с помощью насадки Сокслета. Хорошо известным примером перколяции является также приготовление фильтрованного кофе. Если несколько перколяторов подключаются последовательно (на практике до пяти), это также называется повторным соединением.

При производстве спиртных напитков перколяция означает извлечение эфирных веществ путем «вытеснения». Для этого используются цилиндрические сосуды из меди или керамики, сужающиеся к носику. Измельченные вещества помещают в перколятор между двумя ситами, а сверху наливают спирт с крепостью 30-60 %. Топливо диффундирует и накапливается в течение нескольких дней вместе с растворимыми в спирте веществами, содержащимися в растениях. Через три-шесть дней предварительный перколят медленно сливают. В то же время щелок продолжает поступать через приток до тех пор, пока не будут извлечены компоненты частей растения. Затем следует постперколяция, при которой делается попытка выжать оставшееся топливо путем добавления воды.

Источник

Перколяционный метод

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

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

Общие принципы современного подхода в подземной гидромеханике

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

Транспортные или кондуктивные свойства пористых сред определяются, прежде всего, наличием проводящих поровых каналов. Наиболее простая и удобная модель, позволяющая описывать взаимодействие этих каналов и, как результат, формирование проницаемости пористой среды как макрообъекта – пространственная решетка проводящих капилляров. При этом естественно считать, что радиусы капилляров в такой решетке распределены в соответствии с реальной порометрической кривой – функцией плотности распределения поровых каналов по радиусам f(r).

Перколяционная модель, позволяющая описать проводимость неоднородной среды, если известна функция распределения элементов структуры по величине их собственной проводимости, развивалась в работах автора, например [1, 2]. Она базируется на обобщении модели бесконечного кластера (БК) Шкловского – де-Жена (рис. 1) на случай, когда решетка содержит проводящие элементы, имеющие распределение по величине их собственной проводимости. Исходя из представлений о структуре БК, приходим к задаче определения проводимости его скелета, ответственного за транспортные свойства БК и, соответственно, пористой среды в целом.

Рис. 1. Схематичное представление структуры БК, положенное в основу модели Шкловского — де-Жена без учета извилистости (фрактальности) составляющих его проводящих путей.

Проводящие каналы, по которым осуществляется течение флюида, представляют собой цепочки гидравлически связанных между собой поровых каналов (капилляров) (рис. 2) различного радиуса r (проводимость σ

r λ , где λ определяется характером процесса переноса). Проводимость цепочки будет определяться самым тонким капилляром, поэтому его радиус естественно считать основной характеристикой такой цепочки. Используя его, введем понятие «r-цепочки» – будем так именовать цепочки капилляров, в которых минимальный радиус cоставляющих ее капилляров лежит в диапазоне .

Скелет БК будут составлять r-цепочки всего диапазона области определения функции , причем их количество для каждого будет свое и, вообще говоря, заранее неизвестно. Поэтому необходимо построить некоторый алгоритм определения количества образовавшихся r-цепочек, их проводимостей и последующего суммирования с целью нахождения суммарной проводимости БК. Именно данный алгоритм и является главным элементом перколяционной модели.

Представляемый здесь подход основан на введении определенной систематизации или иерархии r-цепочек, которая позволяет реализовать указанную выше схему суммирования.

Рис. 2. Схема формирования проводящих r-цепочек.

Учет данного обстоятельства позволяет получить окончательное выражение для средней проводимости единицы длины элемента порового 3D пространства (моделирующей его пространственной решетки капилляров)

Численный коэффициент g (порядка единицы), зависящий от типа решетки, призван скорректировать тот факт, что перетоки между проводящими параллельными цепочками не были учтены. При получении выражения (6) учтено численное значение показателя извилистости, связанного с фрактальным поведением элементов склета БК

Проиллюстрируем возможности развитого математического аппарата на примере исследования ряда явлений, возникающих при течении пластовых флюидов к коллекторах.

Эффективная вязкость минерализованной воды при течении в пористой среде.

Чем больше отношение площади поверхности соприкосновения флюида с твердым телом к объему движущегося флюида, тем сильнее течение отличается от объемного или так называемого «течения капельной жидкости».

В ряде экспериментов по исследованию течения флюидов в пористых средах [3, 4, 5] отмечалось превышение эффективной, то есть наблюдаемой в эксперименте, вязкости по сравнению с ее классическим значением для капельной жидкости. Причем особенности течения в микроканалах зависят от природы жидкости и геометрических параметров канала.

В работе [6] построен модель, позволяющая рассчитывать эффективную вязкость при фильтрационном течении электролитов и проведена её верификация путём сопоставления с экспериментальными данными.

Пусть жидкость, текущая по решетке под влиянием приложенного градиента давления, представляет собой раствор симметричного электролита, обладающий диэлектрической проницаемостью и вязкостью . При этом на стенках капилляров существует дзета-потенциал .

Для вывода зависимоcти макропараметров течения от микрохарактеристик среды получим вначале соотношение между внешним перепадом давлений в среде и расходом жидкости по каждому проводящему пути в такой системе.

В качестве первого шага рассмотрим течение по отдельному каналу радиуса а (рис.3).

Рис. 3. Схема электрокинетического течения в капилляре.

Массовая сила, действующая в такой системе, будет определяться взаимодействием потенциала протекания с зарядом ионного облака объемной плотностью . Таким образом, установившееся течение в капилляре будет описываться уравнением Навье-Стокса следующего вида

Здесь – модифицированная функция Бесселя первого рода нулевого порядка, обычно называют безразмерным электрокинетическим радиусом (отношение радиуса канала к толщине ДЭС).

При установившемся электрокинетическом течении в канале ток проводимости равен току протекания, следовательно, общий ток в (10) равен нулю. Учитывая также, что в соответствии с принципом

определяющее эффективную вязкость электролита в процессе его течения в пористой среде.

На рис. 4 отражены результаты расчёта величины в случае, когда максимальный радиус порометрической кривой был постоянным, а минимальный радиус варьировался. Видно, что влияние ДЭС на течение жидкости возрастает при увеличении доли капилляров с меньшим радиусом.

Установлено, что при наличии достаточного числа каналов с радиусом, сравнимым с толщиной ДЭС, скорость фильтрации электролита значительно (в разы) меньше предсказываемой теорией из расчётов, не учитывающих влияние ДЭС, что интерпретируется как возрастание вязкости (эффект электровязкости). Результаты расчетов согласуются с данными проведенного авторами эксперимента.

Соответственно, при проектировании разработки необходимо учитывать, что величина вязкости капельной жидкости может существенно отличаться от её эффективной вязкости при течении в пористой среде.

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

В процессе разработки месторождения такими методами, как циклическое заводнение и смена направления фильтрационных потоков, вытеснение нефти водой сменяется вытеснением воды нефтью и наоборот. Такое изменение в характере движения влияет на зависимости относительных фазовых проницаемостей (ОФП) от водонасыщенности. Это явление называется гистерезисом ОФП при дренаже и пропитке. Пренебрежение этим явлением при гидродинамическом моделировании нефтяных месторождений приводит к существенным ошибкам.

Наличие гистерезиса ОФП отмечалось в ряде экспериментальных работ 7.

Математическая модель данного явления построена и исследована в работах 11.

В качестве основных механизмов возникновения гистерезиса были учтены следующие явления. В процессе дренажа активные компоненты нефти адсорбируются на поверхности породообразующих минералов [7], что приводит к формированию нефтяной нанопленки, не меняющей вид f(r), т.к. толщина пленки на несколько порядков меньше радиуса самых мелких капилляров. Это приводит к гидрофобизации поверхности капилляров, что рассматривается в качестве первого механизма возникновения гистерезиса. Перемешивание флюидов в процессе дренажа может приводить к изменению реологических свойств обеих жидкостей. Таким образом, второй механизм заключается в изменении реологических свойств нефти и воды.

В результате прохождения нефти через изначально гидрофильную пористую среду поверхностные свойства некоторой части капилляров меняются. Для описания движения флюидов в такой среде используется перколяционная модель среды с микрогетерогенной смачиваемостью [1], в которой вводятся параметры κ и α: κ – доля капилляров с неизмененными свойствами поверхности (тип 1), (1-κ) – доля капилляров с измененными свойствами поверхности (тип 2), α=cosθ₂/cosθ₁, где θ₁ и θ₂ – углы смачивания в капиллярах 1-го и 2-го типа соответственно. При этом поверхностные свойства могут поменяться двумя принципиально разными способами (рис. 2): либо уменьшится «степень гидрофильности» поверхности капилляров (0 2+ или Mg 2+ ;

· Закачиваемая вода имеет пониженную концентрацию солей по сравнению с пластовой водой.

Расчет влияния закачки в коллектор воды с ионным составом и минерализацией, отличной от состава и минерализации пластовой воды, на эффективность заводнения представлен в работе [15].

При этом был проведен анализ изменения фильтрационно-ёмкостных свойств (ФЕС) глиносодержащего коллектора в процессе смены гидрохимического режима разработки пласта.

Оказалось, что пористость может изменяться в пределах 15-20% в зависимости от концентрации солей и ионного состава раствора и обменного комплекса глин.

Расчеты показывают, что вследствие усадки глин, при взаимодействии Na-бентонита с раствором хлорида калия происходит рост фазовой проницаемости по воде и снижение по нефти. Вследствие набухания, при взаимодействии K-бентонита с раствором хлорида натрия, происходит снижение фазовой проницаемости по воде и рост по нефти. Результаты теоретических расчетов подтверждаются экспериментальными исследованиями поведения ОФП в глиносодержащих пластах при фильтрации различных растворов электролитов.

В случае разбухания глин коллектор переходит в разряд тонкопоровых, поэтому необходим учет капиллярных сил. В связи с этим фильтрация флюидов рассматривалась в рамках модели Раппопорта-Лиса, а полная система дифференциальных уравнений в частных производных для определения концентрации солей в электролите C и водонасыщенности S имеет вид

где S_r – доля адсорбированной воды в породе (знак «+» перед S_r при набухании глин, «–» в случае усадки), F(S,C) – функция Бакли-Леверетта, q – удельный суммарный расход фаз, P_k(S,C) – функция капиллярного давления, J(S) – функция Леверетта, χ – поверхностное натяжение фаз, M_i(S,C) – подвижность фазы (i = 1 соответствует водной фазе), w_i – скорость фильтрации фазы, ΔP_i – разность между давлением на контуре и на стенке скважины в i-й фазе, R_k – радиус контура пласта, r_w – радиус скважины, K(C) – абсолютная проницаемость, f(ρ(C)) – функция плотности распределения капилляров по радиусам ρ, γ´ – коэффициент, отражающий наличие перетоков между проводящими цепочками (порядка единицы), ν´ – индекс радиуса корреляции, ρ_c – критический радиус капилляров, l – средняя длина капилляров, D * – коэффициент, включающий в себя диффузию и гидродинамическую дисперсию, R´(C) – скорость гидратации. Действием гравитационных сил пренебрегаем.

Видно, что в сравнении со стандартной закачкой пластовой воды, закачка хлорнатриевого раствора дает прирост нефтеотдачи на 6%, а закачка слабоминерализованной воды дает прирост на 15%. Кроме того, смена гидрохимического режима значительно интенсифицирует разработку пласта. Так, нефтеизвлечение 60%, достигнутое закачкой пластовой воды за 77 лет, при закачке хлорнатриевого раствора достигается на 44 года раньше, а при закачке слабоминерализованной воды раньше на 52 года.

Циклическая закачка раствора электролита представляется более экономически выгодной по сравнению с постоянной.

Видно, что при частоте 0,01 год -1 концентрация электролита в пласте незначительна, и, следовательно, фазовая проницаемость по нефти низкая, по воде – высокая. Кроме того, цикл закачки пластовой воды пришелся на период после прорыва воды к скважине, что усилило темп обводнения. Увеличение частоты приводит к тому, что после прорыва воды и на конец разработки пришелся цикл закачки электролита, что способствовало довытеснению нефти. При высоких частотах закачки (ω ≥ 1 год -1 ) амплитуда колебания концентрации электролита в пласте падает, диапазон изменения ФЕС снижается, и, как следствие, снижается эффективность нефтевытеснения.

Проведенный анализ ионообменных процессов при течении минерализованных жидкостей в глиносодержащей пористой среде показал, что изменение объема глинистой компоненты приводит к изменению пористости на 15-20%. В результате объем пор, содержащих воду, уменьшается, и последняя «выдавливается» в нефтесодержащую часть порового пространства. Это приводит к дополнительному вытеснению нефти по сравнению со стандартным процессом заводнения пластовой водой.

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

Закачка в глиносодержащий пласт растворов солей высокоактивных металлов лития и натрия увеличивает нефтеотдачу на 10 и 8% соответственно (рис. 11). Применение циклической закачки позволяет вдвое снизить объем реагентов – солей лития и натрия – при потерях в нефтеотдаче по сравнению с постоянной закачкой на 2 и 1% соответственно.

1. V.I.Selyakov, V.V.Kadet. Percolation Models for Transport in Porous Media. With Applications to Reservoir Engineering. Kluwer Academic Publishers. Dordrecht/ Boston/ London, 1996.

2. V.V.Kadet, A.A. Maximenko. Determination of relative permeabilities using the network models of porous media. Journal of Petroleum Science and Engineering, Vol. 28 Issued 3, pp. 145-152, November 2000.

3. Tuckermann, D. B., and Pease, R. F. W. High-performance Heat Sinking for VLSI // IEEE Electron. Device Lett. , 1981. № 2(5). P. 126.

4. Pfahler, J. N., Liquid Transport in Micron and Submicron Size Channels, Ph.D. thesis, Department of Mechanical Engeneering and Applied Mechanics, Univ. of Pennsylvania, 1992.

5. Peng, X. F., Peterson, G. P., and Wang, B. X. Frictional flow characteristics of water flowing through rectangular microchannels // Exp. Heat Transfer 7, 1994. P. 249.

6. V.V.Kadet, A.S.Koryuzlov. Effective Viscosity of Mineralized Water Flowing in a Porous Medium: Theory and Experiment. Theoretical Foundations of Chemical Engineerin, 2008, Vol. 42, No. 6, pp. 899-904.

7. Variations in bounding and scanning relative permeability curves with different carbonate rock types / Dernaika M.R., Basioni M.A., Dawoud A., Kalam M.Z., Skjæveland S.M. // Paper SPE 162265 presented at the 2012 Abu Dhabi International Petroleum Conference and Exhibition, Abu Dhabi, UAE, Nov. 11-14.

8. Wei J. Z., Lile O. B. Influence of wettability and saturation sequence on relative permeability hysteresis in unconsolidated porous media // [Electron. resource]. http://www.onepetro.org/mslib/app/Preview. do?paperNumber=00 025282&societyCode=SPE.

9. Braun E. M., Holland R. F. Relative permeability hysteresis: laboratory measurements and a conceptual model // SPE Reservoir Engineering., 1995, V. 10, № 3, рр.222-228.

10. V.V Kadet, A.M. Galechyan. Percolation Modeling of Relative Permeability Hysteresis. Journal of Petroleum Science and Engineering, vol. 119, pp. 139-148, July 2014.

11. Кадет В.В., Галечан А.М. Учет реологии флюидов в гидрофобизационной модели гистерезиса относительных фазовых проницаемостей. ПМТФ, 2017, Т. 58, N.6. С.58 — 68.

12. Tang G. Q., Morrow N. R. Influence of brine composition and fines migration on crude oil/brine/rock interactions and oil recovery // Journal of Petroleum Science and Engineering. 1999. V. 24.P.99-111.

13. Lager A., Webb K.J., Black C.J.J., Singleton M., Sorbie K.S. Low salinity recovery – an experimental investigation // Petrophysics. 2008. V. 49. P. 28-35.

14. Aboulghasem Kazemi Nia Korrani, Gary R. Coupled Geochemical-Based Modeling of LSW // SPE 169115-MS.

Источник