Тема 2. дифференциальные уравнения

Геофильтрации и геомиграции

Теоретическая часть

Емкостные свойства горных пород

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

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

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

Тема 2. дифференциальные уравнения - student2.ru = Wo- Wн; (2.1)

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

Тема 2. дифференциальные уравнения - student2.ru ' = Wo- Wе. (2.2)

Напомним, что под полной влагоемкостью (Wo)понимают максимально возможное количество воды (все виды воды) при полном заполнении свободного пространства в минеральном скелете горной породы, а под наименьшей влагоемкостью (водоудерживающая способность горной породы) (Wн) понимают влажность, определяемая максимальным количеством прочносвязанной, осмотической, капиллярно-разобщенной и капиллярно-подвешенной воды.

Таким образом, величина водоотдачи зависит от свойств самой горной породы (общая пустотность, размеры пустот, структура "свободного" пространства, минеральный состав и др.), определяющих суммарное количество различных видов связанной воды (Wн) и объем воды при полном насыщении (Wo), и в этом смысле может рассматриваться как объективное свойство (параметр) данной породы. В определенной мере водоотдача зависит также от скорости осушения горной породы (скорости снижения уровня полного насыщения), поскольку эта скорость может определять (при прочих равных условиях) формирование различных объемов капиллярной воды (капиллярно-подвешенная вода).

Величина недостатка насыщения ( Тема 2. дифференциальные уравнения - student2.ru ') в реальных условиях зависит от залегания и изменения естественной влажности (We), поэтому этот показатель не может рассматриваться как гидрогеологический параметр горной породы, а является комплексной характеристикой ее свойств в определенных условиях залегания .

В заключении рассмотрения фильтрационных и емкостных параметров горных пород приведем классификацию горных пород по проницаемости (табл.2.1).

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

Запишем теперь выражение для относительного изменения объема воды V0 в рассматриваемом объеме породы Vп, имея в виду, что V0 = n Vп:

Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru d H, (2.3)

где Тема 2. дифференциальные уравнения - student2.ru - коэффициент упругоемкости пластахарактеризует объем жидкости, который может быть получен с единицы объема горной породы за счет упругих свойств как горной породы, так и воды при снижении напоров на 1 м.

Таблица 2.1

Классификация горных пород по проницаемости

Группа пород Проницаемость (Кф, м/сут) Гравитационная емкость (водоотдача), m Типы горных пород
Очень высокопроницаемые более 100 0.2 - 0.3 0.1 - 0.15 Гравийно-галечниковые образования, интенсивно закарстованные породы (известняки, гипсы)
Хорошо (высоко) проницаемые 10-100 0.2 - 0.25 0.1 - 0.15 Крупнозернистые и грубозернистые пески, интенсивно трещиноватые скальные породы (песчаники, граниты, базальты )
Проницаемые 0.1 - 10 0.15 - 0.2 0.05 - 0.1 Разнозернистые, тонко- и мелкозернистые глинистые пески, слаботрещиноватые и закарстованные горные породы
Слабопроницаемые 10-6 - 10-2 0.01 - 0.03 Суглинки, глины песчанистые или микротрещиноватые, очень слаботрещиноватые скальные породы, незакарсто- ванные известняки
Практически непроницаемые меньше 10-6 менее 0.01 Плотные нетрещиноватые глины, практически нетрещиноватые скальные породы

Численно он равен

Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru , (2.4)

где g0 - объемная масса жидкости; е - коэффициент пористости, е = n/(1-n); Eв= модуль Юнга для воды (для чистой воды Ев= 2 Тема 2. дифференциальные уравнения - student2.ru 103МПа); ас - коэффициент сжимаемости горных пород.

Коэффициент упругоемкости Тема 2. дифференциальные уравнения - student2.ru имеет размерность метр в минус первой степени. В целом значения коэффициентов упругоемкости относительно невелики (табл.2.2).

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

Таблица 2.2

Значения коэффициент упругоемкости для различных типов пород

Порода Коэффициент упругоемкости Тема 2. дифференциальные уравнения - student2.ru-1)
Пески (0.5 Тема 2. дифференциальные уравнения - student2.ru 5) Тема 2. дифференциальные уравнения - student2.ru 10-4
Супеси и суглинки 10-4-10-3
Трещиноватые породы 10-5-10-6

Коэффициент упругой водоотдачи Тема 2. дифференциальные уравнения - student2.ru * связан с коэффициентом упругоемкости следующим образом:

Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru Тема 2. дифференциальные уравнения - student2.ru m, (2.5)

где m - мощность водоносного горизонта, т.к. в соответствии с таблицей значения Тема 2. дифференциальные уравнения - student2.ru незначительны, то для безнапорных горизонтов, в которых при снижении уровней проявляются и гравитационные и упругие запасы, последними обычно пренебрегают, т.е. считают

Тема 2. дифференциальные уравнения - student2.ru << Тема 2. дифференциальные уравнения - student2.ru . (2.6)

Дифференциальное уравнение упругого режима фильтрации

Упругий режим геофильтрации учитывает сжимаемость пласта и воды, т.е. уравнения состояния имеют вид

- для воды (34)

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

С учетом вышесказанного дифференциальное уравнение упругого режима фильтрации можно представить согласно 2.7

Тема 2. дифференциальные уравнения - student2.ru = m Тема 2. дифференциальные уравнения - student2.ru Тема 2. дифференциальные уравнения - student2.ru (2.7)

Для однородного изотропного пласта

Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru Тема 2. дифференциальные уравнения - student2.ru , (2.8)

где

a* = Тема 2. дифференциальные уравнения - student2.ru . (2.9)

Уравнение вида (2.8), называется уравнением Фурье и широко исследовано в теории теплопроводности.

Таким образом, в отличии от жесткого режима движения в напорном пласте, уравнения упругого режима напорной фильтрации содержат производную по времени, т.е. в этом случае Н=Н (x, y, z, t), и движение является нестационарным. Физически это означает, что по мере уменьшения напоров во времени в водоносном пласте постепенно срабатываются его упругие запасы. Высвобождающиеся при этом объемы воды "вносятся" в общий баланс фильтрационного потока в водоносном пласте. Реакция от возмущения напоров на границе или в какой-либо области рассматриваемого пласта распространяется от границы (области) по пласту постепенно, причем скорость распространения тем больше, чем выше проницаемость и чем меньше упругоемкость горной породы. Следовательно, отношение a* = Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru является показателем скорости изменения напора (гидростатического давления) в пласте. Соответственно величина а* получила название коэффициента пьезопроводности ; ее размерность м2/сут.

Плановая фильтрация в изолированном напорном пласте

Уравнение напорной фильтрации для жесткого режима:

Тема 2. дифференциальные уравнения - student2.ru = 0, (2.10)

где Tx = kxm и Ty = kym - коэффициенты водопроводимости пласта в направлениях осей x и y (совпадающих в общем случае с направлением главных осей анизотропии).

При Tx = Ty = T = const получаем:

Тема 2. дифференциальные уравнения - student2.ru = 0. (2.11)

В упругом режиме имеем

Тема 2. дифференциальные уравнения - student2.ru = m* Тема 2. дифференциальные уравнения - student2.ru . (2.12)

где m* = m Тема 2. дифференциальные уравнения - student2.ru m - коэффициент упругой водоотдачи пласта

Для однородного изотропного пласта постоянной мощности:

Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru Тема 2. дифференциальные уравнения - student2.ru , (2.13)

где а* - коэффициент пьезопроводности, который в данном случае может быть выражен в виде:

а* = T / m*. (2.14)

Плановая напорная фильтрация при наличии перетекания

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

Тема 2. дифференциальные уравнения - student2.ru

Рис.2.1 Схема к выводу уравнения неразрывности в пласте с перетеканием

При расчетах подобных водоносных систем принимаются следующие предположения, известные как предпосылки перетекания (предпосылки Мятиева-Гиринского):

1) движение в водоносных пластах является плановым (линии тока параллельны напластованию);

2) в разделяющем слое линии тока перпендикулярны к напластованию; физически эта предпосылка вполне объяснима: вода стремится пройти участок с большим сопротивлением (водоупорный слой) по кратчайшему пути.

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

Ip = Тема 2. дифференциальные уравнения - student2.ru .

Замечание. Одновременно мы тем самым предположили, что режим фильтрации в разделяющем слое является жестким, т.е. в этом слое мгновенно устанавливается распределение напоров в соответствии с напорами на его кровли Н' и подошве Н. Иначе говоря, мы пренебрегли упругими запасами воды в разделяющем слое.

Следовательно, через верхнюю грань столбика поступает дополнительное питание, равное

eп = kp Тема 2. дифференциальные уравнения - student2.ru

на единицу площади пласта в единицу времени; величина eп входит в качестве дополнительного члена в уравнение неразрывности, так что вместо (41) получаем:

Тема 2. дифференциальные уравнения - student2.ru + Тема 2. дифференциальные уравнения - student2.ru - r eп = 0 .

Соответственно преобразуются и дифференциальные уравнения фильтрации. Так, вместо результирующего уравнения (2.13) получаем

Тема 2. дифференциальные уравнения - student2.ru + Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru Тема 2. дифференциальные уравнения - student2.ru , (2.15)

где B = Тема 2. дифференциальные уравнения - student2.ru - так называемый параметр (фактор) перетекания, имеющий размерность длины. Чем меньше величина В, тем интенсивнее, при прочих равных условиях, идет перетекание.

Конвективный перенос в подземных водах

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

Vд = Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru (2.16)

где na – так называемая активная пористость (трещиноватость).

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

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

Тема 2. дифференциальные уравнения - student2.ru

Конвективный перенос сорбируемого вещества

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

При достаточно длительном взаимодействии раствора и породы зависимость между концентрацией компонента в растворе c и на твердой фазе N оказывается для данной температуры однозначной

N0 = f( c), (2.17)

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

Для мало концентрированных растворов достаточно характерна линейная изотерма Генри:

N0( с ) = kг с, (2.18)

где безразмерная константа kг определяющая отношение концентраций No и c, называется коэффициентом распределения для данной пары “раствор – порода”. Для несорбируемых компонентов kг = 0, а для хорошо сорбируемых – величина kгизмеряется сотнями и даже тысячами относительных единиц.

С другой стороны, ясно, что сорбция протекает во времени. Полная сорбционная емкость породы реализуется не сразу. Количество сорбируемого компонента, усвоенного единицей объема горной породы N(c, t), достигает своего предельного значения N0 (с) только после достаточно длительного соприкосновения минеральных зерен с внутрипоровым раствором, имеющим концентрацию с. Это обстоятельство приближенно отражается уравнением кинетики сорбции, связывающим линейно интенсивность сорбции с отклонением от равновесной концентрации:

Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru [ c - cc ( N ) ], (2.19)

где с – текущая средняя концентрация порового раствора; cc ( N ) – концентрация раствора непосредственно на контакте с поверхностью минерального скелета (связанная с текущей величиной N(t) (2.17); Тема 2. дифференциальные уравнения - student2.ru – константа скорости массообмена, имеющая обычно порядок десятков – сотен суток в минус первой степени. Уточним, что речь здесь идет, конечно, о внутрипоровой кинетике сорбции, т.е. о скорости обмена веществом между минеральными зернами и смежными с ними порами.

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

Рассмотрим в таком варианте процесс поршневого вытеснения впределах фиксированной трубки тока (фильтрационный поток – стационарный с расходом Q) . Взаимодействующие жидкости считаем однокомпонентными, так что количество солей однозначно характеризуется концентрацией: для вытесняемой (пресной) воды она равна с0, для вытесняющей (соленой) – с0. На фронте вытеснения имеет место сорбция, протекающая мгновенно. Составим баланс соли для элемента трубки (рис.) длиной dl, отвечающей перемещению фронта вытеснения за время dt.

Тема 2. дифференциальные уравнения - student2.ru

Количество соли, поступившей слева, равно с0 Q dt, а поток соли вправо равен с0 Q dt. Накопление соли в поровом пространстве элемента dlравно (с0 – с0 ) n w dl(w – усредненное поперечное сечение трубки), а накопление соли на минеральном скелете за счет сорбции, согласно зависимости (2.17), равно (kгc0– kгc0) w dl. Уравнение баланса имеет вид,

( c0 - c0) Q dt = (с0 – с0 ) n w dl + kг (c0–c0) w dl. (2.20)

Отсюда

Тема 2. дифференциальные уравнения - student2.ru , (2.21)

где vп – скорость перемещения фронта вытеснения. При отсутствии сорбции vп = vд, где средняя действительная скорость движения жидкости vдопределяется формулой (2.16).

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

nэ = n + kг (2.22)

которая формально выражает общую емкость водоносных пород по отношению к рассматриваемому сорбируемому компоненту. Подчеркнем, что при больших kг величина n э>>1.

Диффузионно - дисперсионные процессы.

Молекулярная диффузия и гидродисперсия

При наличии градиента концентрации возникает самопроизвольный перенос вещества под действием молекулярных сил. Общий поток вещества Qc в направлении Lчерез расчетное сечение wl при этом ориентирован в сторону уменьшения концентрации с и определяется законом Фика:

Тема 2. дифференциальные уравнения - student2.ru = - Dм wl Тема 2. дифференциальные уравнения - student2.ru , (2.23)

где Dм– коэффициент молекулярной диффузии в пористой среде, зависящий от ее структурных особенностей и физических свойств поровых растворов (его характерные значения для песчано-глинистых пород имеют порядок 10-5 м2/сут).

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

В конечном счете внешним проявлением диффузии во всех случаях оказывается "размыв" границы между вытесняемой и вытесняющей жидкостями, образование вблизи фронта вытеснения переходной зоны Lп пределах которой концентрация вещества изменяется постепенно от с0до с0 (рис.), в отличие от резкой границы, характерной для чисто конвективного переноса (т.е. для схемы поршневого вытеснения). По такому внешнему проявлению с молекулярной диффузией ассоциируется еще один механизм переноса вещества – механическая дисперсия.

Тема 2. дифференциальные уравнения - student2.ru

Однако природа этого явления существенно иная: механическая дисперсия обусловлена неоднородностью поля действительных скоростей. Если, например, все окрашенные частицы воды на (рис) начинают одновременно двигаться в фильтрационном потоке, имеющем среднюю действительную скорость vд, от поперечного сечения АА в направлении ВВ, то за время t=L/vд одни из этих частиц пройдут путь существенно больший, чем L, а другие – меньший. В результате вместо резкого контакта окрашенной и неокрашенной жидкости, характерного для начального момента (сечение АА), вблизи сечения ВВэтот контакт окажется размытым – вследствие рассеяния окрашенных частиц фильтрационным потоком. Согласно экспериментам, подобное рассеяние, называемое фильтрационной диффузией, или механической дисперсией, также подчиняется закону Фика (2.23), однако при замене коэффициента молекулярной диффузии Dм на другой параметр – коэффициент механической дисперсии Dд.

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

Dд = Тема 2. дифференциальные уравнения - student2.ru v, (2.24)

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

Для однородных песков параметр Тема 2. дифференциальные уравнения - student2.ru имеет порядок миллиметров, а для трещиноватых пород – порядок расстояний между основными водопроводящими трещинами (поэтому нередки случаи, когда Тема 2. дифференциальные уравнения - student2.ru измеряется метрами, а то и десятками метров.

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

Аналогия в характере проявления молекулярной диффузии и механической дисперсии способствовала их объединению через суммарный коэффициент:

D = Dм + Dд = Dм + Тема 2. дифференциальные уравнения - student2.ru v, (2.25)

который далее для краткости будет именоваться коэффициентом дисперсии.

Соответственно, суммарное рассеяние вещества, вызываемое гидродисперсией, или просто дисперсией, определяется законом Фика в виде:

Тема 2. дифференциальные уравнения - student2.ru = - Dм wl Тема 2. дифференциальные уравнения - student2.ru , (2.26)

Феноменологическая закономерность (2.26) является, таким образом, некоторым усредненным детерминированным отражением процесса, имеющего на самом деле ярко выраженный вероятностный характер. Благодаря подобным усреднениям, мы отказываемся от стохастического анализа процессов массопереноса на микроуровне, связанного с необходимостью использования статистической теории фильтрации. Важно отметить, что дисперсия может идти не только вдоль направления фильтрационного потока (как мы до сих пор предполагали), но и в перпендикулярном к нему направлении. Такая – поперечная – дисперсия также подчиняется закону вида (2.26), но с иным коэффициентом D', имеющим, подобно (2.25), вид:

D' = Dм +v, (2.27)

геометрический параметр Тема 2. дифференциальные уравнения - student2.ru при поперечной дисперсии примерно на порядок меньше, чем при продольной, хотя опытного материала по этому поводу очень мало.

Значение поперечной дисперсии проще всего проиллюстрировать, запуская меченые частицы жидкости не по всей ширине, а лишь на ограниченном участке фильтрационного потока. Поперечная дисперсия приводит к распространению меченых частиц за пределами зоны, ограниченной крайними траекториями фильтрации меченой жидкости, к появлению здесь переходных зон и к уменьшению скоростей продвижения фронта вытеснения в продольном направлении. Все это хорошо иллюстрируется (рис.), полученным по конкретному эксперименту. Слагаемые в формулах (2.25, 2.27) отражают интенсивность молекулярно-диффузионного и механического дисперсионного массопереноса: доля каждого из них в общем балансе массового потока определяется скоростью фильтрации v и структурой горных пород ( для однородных пород - характерным размером минеральных зерен dз).Так, из результатов лабораторных экспериментов следует, что при vdз, порядка 0,0001 м2/сут и менее доминирующей является молекулярная диффузия, а при параметре vdз измеряемом сотыми – десятыми долями квадратных метров в сутки, резко преобладает уже механическая. Поэтому, с определенной долью погрешности, можно говорить, что в породах глинистого состава обычно можно учитывать только молекулярную диффузию, в трещиноватых породах - только механическую дисперсию, но в остальных случаях это требует дополнительного обоснования.

Конвективно-дисперсионный перенос в однородных водоносных пластах

При выводе фундаментального решения задачи конвективно-дисперсионного переноса в однородных водоносных пластах мы ограничимся рассмотрением жидкостей с постоянной плотностью (и вязкостью). Это позволяет значительно упростить анализ миграционных процессов: «фильтрационная часть» задачи может решаться заранее, отдельно,– независимо от задачи переноса; поэтому фильтрационное поле считается нами заданным.

Фундаментальное решение

Тема 2. дифференциальные уравнения - student2.ru

Рассмотрим совместное проявление конвекции и дисперсии в однородном пласте, сложенном гомогенными водоносными породами – пористыми или «чисто» трещиноватыми (рис.а). Фильтрационный поток считаем одномерным (плоскопараллельным) и стационарным. Исходная концентрация вещества повсеместно равна c0. В момент t = 0концентрация на левой границе принимает постоянное значение c0 и соленые воды начинают перемещаться по пласту в направлении оси X. Поток соли qcчерез произвольное сечение пласта обусловлен конвекцией ( Тема 2. дифференциальные уравнения - student2.ru = cq, где c(x,t) – текущее значение концентрации; q – удельный расход фильтрационного потока) и дисперсией ( Тема 2. дифференциальные уравнения - student2.ru = – Dm Тема 2. дифференциальные уравнения - student2.ru – согласно закону Фика):

qс = Тема 2. дифференциальные уравнения - student2.ru + Тема 2. дифференциальные уравнения - student2.ru = cq – Dm Тема 2. дифференциальные уравнения - student2.ru (2.28)

Составим уравнение неразрывности для элемента dx:

qc dt - [ qc - Тема 2. дифференциальные уравнения - student2.ru ] dt = Тема 2. дифференциальные уравнения - student2.ru ( n m dx c) dt, (2.29)

где справа записано приращение количества соли в элементе dx(с объемом порового пространства, равным nmdx, причем при наличии сорбции величина nзаменяется на nэсогласно (2.22) за время dt, а слева – разность между количеством соли, поступившей в этот элемент и вытекшей из него за то же время dt. Подставляя сюда выражение (2.28) для qc приходим к дифференциальному уравнению конвективно-дисперснонного переноса относительно неизвестной концентрации c(х, t)

n Тема 2. дифференциальные уравнения - student2.ru + v Тема 2. дифференциальные уравнения - student2.ru = D Тема 2. дифференциальные уравнения - student2.ru . (2.30)

Начальные и граничные условия имеют вид:

Тема 2. дифференциальные уравнения - student2.ru (2.31)

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

Тема 2. дифференциальные уравнения - student2.ru = ( c - c0 ) / ( c0 + c0)

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

n Тема 2. дифференциальные уравнения - student2.ru + v Тема 2. дифференциальные уравнения - student2.ru = D Тема 2. дифференциальные уравнения - student2.ru ; (2.32)

Тема 2. дифференциальные уравнения - student2.ru (2.33)

Решение этой задачи, формально напоминающей фундаментальную задачу фильтрации , можно представить в виде

Тема 2. дифференциальные уравнения - student2.ru = 0,5 ( erfc Тема 2. дифференциальные уравнения - student2.ru + e Тема 2. дифференциальные уравнения - student2.ru erfc Тема 2. дифференциальные уравнения - student2.ru ), (2.34)

где

Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru ; Тема 2. дифференциальные уравнения - student2.ru = Тема 2. дифференциальные уравнения - student2.ru ; Тема 2. дифференциальные уравнения - student2.ru ,

а erfc(z) – табличная функция.

При достаточно больших q вторым членом в правой части (2.34) можно пренебречь, внося при этом погрешность е порядка: e=0,3/ Тема 2. дифференциальные уравнения - student2.ru . Тогда решение (2.34) принимает вид:

Тема 2. дифференциальные уравнения - student2.ru = 0,5 erfc Тема 2. дифференциальные уравнения - student2.ru . (2.35)

Графическое представление решения (2.35) дается на (рис б). На графике выделяются три зоны: I – зона вытесняющего раствора, III – зона вытесняемого раствора и II – переходная зона ', в пределах которой относительная концентрация меняется от значений, близких к 1, до значений, близких к 0.

На рисунке пунктиром показан также график изменения концентрации при поршневом вытеснении, когда дисперсия отсутствует (D = 0). Понятно, что точка xп, согласно (2.16), отвечает условию:

xп = Тема 2. дифференциальные уравнения - student2.ru

и, следовательно, в случае, описываемом решением (2.35),

Тема 2. дифференциальные уравнения - student2.ru ( xп ) = 0,5 erfc Тема 2. дифференциальные уравнения - student2.ru 0 = 0,5, (2.36)

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

Ориентировочно можно полагать, что дисперсией допустимо пренебречь при

Тема 2. дифференциальные уравнения - student2.ru < e, (2.37)

где е – некоторое малое число. Например, полагая в=1 %, приходим к условию:

Тема 2. дифференциальные уравнения - student2.ru < 0,0001 xп. (2.38)

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

Литература

1. Мироненко, В.А. Динамика подземных вод. [Текст] / В.А. Мироненко - М.:Изд-во МГУ, 1996.-519 с.

4. Миpоненко, В.А. Теоpия и методы интеpпpетации опытно-фильтpационных pабот. [Текст] / В.А.Миpоненко, В.М.Шестаков - М.: Hедpа, 1978. - 325 с.

5. Шестаков, В.М. Динамика подземных вод. [Текст] / В.М. Шестаков -М.:Изд-во МГУ,1979. – 418 с.

Практическая часть

Наши рекомендации