Моделирование миграции подземных вод

РАСЧЕТЫ ИЗОБРАЖЕНИИ ВЫХОДНЫХ КОНЦЕНТРАЦИЙ ПО ЛАПЛАСУ—КАРСОНУ

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

Vx (v2—v,) e ptl'\...... . . + (Vn

(П.1>

Где n — общее число отрезков ло­маной линии, причем обычно можно - разбивать такие кривые на четы­ре— шесть отрезков (п = 4—6).

Можно также использовать об­щие формулы механических квад­ратур [13] со сдвигом начала рас­чета на время to. Тогда при расче­тах по пятиточечной формуле полу­чим

(П. 2)

Где tk = to-\-JktP, причем значения tk и Ak принимаются по данным табл. 1. При таких расчетах рекомендуется задавать значения tp в ин­тервалах от 0,5^о ДО 5, соответствующего с—0,5; при tp>i0,5 фор­мулой (П.2) пользоваться не следует.

При опыте с трассерным «пакетом», для которого типичная вы­ходная кривая представлена на рис. 36, расчеты изображений ре­комендуется проводить по формуле (П.2) при максимальном зна­чении tp=0,4tmax (где t„,ах — максимальное время достоверной фиксации выходной концентрации).

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

С = см ехр [— а, (t — tM) —Ozit — tuf ],

P

-pK

V*e

РАСЧЕТЫ ИЗОБРАЖЕНИИ ВЫХОДНЫХ КОНЦЕНТРАЦИЙ ПО ЛАПЛАСУ—КАРСОНУ

Рис. 61. График изменения кон­центрации во времени. Скорости изменения концентраций: Pl-C|/(*|—<!>; 02= (C2—Ct)/(І2—1\)

Параметры которой аг и а2 определяются по графику, построенно­му в координатах 0= (/—/м)_1^см/с и (t—tM), где опытные точки должны лечь на прямую линию, дающую значения aj=0(O) и

Таблица 1

Значения Ан и h в формуле (П.2)

K

1

2

3

4

5

H

0,335

1,13

2,4

4,17

6,49

Ah

0,403

0,332

0,138

0,0316

0,004

Таблица 2

Значения коэффициентов Ai(xu)

Величины при различных значениях

І

0,2

0,3

0,4

0,5

0,6

0,8

1,0

1

—2,173

—2,118

-.1,663

—1,074

—0,457

0,695

1,631

2

5,247

4,185

2,554

0,828

—0,820

—3,639

—5,710

3

—2,465

—1,178

0,267

1,624

2,817

4,636

5,712

4

0,392

0,108

—0,164

—0,385

-0,547

-0,695

—0,632

А2=2,ЗА9/А (t—t,.,), причем см и /„ - максимальная концентрация и время ее прохождения. Расчет изображения при этом проводится по формуле

Ср =-1- ехр (- £ + а,4 - а242) Erfc,

Где <7=а,-2ааМ-р.

Обратный переход от изображения по Лапласу—Карсону F(p) к его оригиналу f(t) для монотонных выпуклых кривых JI. К. Гох - берг рекомендует [13] проводить по формуле

ЛО 4-С4)F(Pl), (П. З)

Где Ai(xk)—коэффициенты, которые находят из табл. 2.

Для расчета по формуле (П. З) прежде всего выбирают макси­мальное расчетное время процесса ^пах, по которому затем рассчи­тывают четыре параметра преобразования pi, причем p\ — 8/tmax, p2==4/^max, Рз—2/fmax, p4=l/*max. ПрИ ЭТИХ ЗНаЧЄНИЯХ pi ВЫЧИСЛЯ­ЮТ изображение F(pi). Далее для каждого хк из табл. 2 берутся коэффициенты Аі(хк) и по уравнению (П. З) находятся значения f(t) При t—Xktmax. ТаКИМ ОбрЭЗОМ ПОЛуЧЭЮТ ЗНЭЧеНИЯ f (t) В ИНТер - ВаЛе ВреМеНИ 0,2^гпах t ^гпах - Если необходимо получить значе­ния f(t) в другом интервале времени, то соответственно следует изменить задаваемое значение /max.

[1] Маркс К., Энгельс Ф. Поли. собр. соч., т. 20, с. Б87.

[2] Раздел составлен по материалам работ [5, 9, 16].

[3] Здесь и далее буквой (г) обозначено газообразное состояние вещества, (к) — твердое и (в) — жидкое.

[4] Сорбционная емкость N связывается с емкостью поглощения ЕП соотно­шением М=уа ЕП, где Yn — объемная масса породы.

[5] Нередко используется несколько иная запись закона Фика для пористой среды', при которой площадь поперечного сечения задается только в пределах сечения порового пространства: vd~—Dit*n0grad с(поЬи*—Ви)< Однако, сохра­няя используемые в теории фильтрации положения механики сплошной среды, будем в дальнейшем считать предпочтительной форму записи (4.3).

[6] Энгельс Ф. Диалектика природы. М., Политиздат, 1982, с. 193.

[7] Выражения (4.7)—(4.10) для Dt приведены здесь в несколько трансфор­мированном виде, исходящем из представления закона Фика в форме (4.3) и - принимая в качестве кинематической характеристики скорость фильтрации. В за­рубежной литературе обычно используется коэффициент гидродисперсии Кіж - Dijn.

[8] + (то*/то) mj Сопоставляя выражения (4,18) и (4.15), находим соотношение mQ* 1

.cth то

Т0 т0

Из которого численным путем устанавливается, что при m0< 1 с точностью до 10 % получается постоянное значение т0*/т0—0,33. Таким образом, для блоков

[9] Под плотностью переноса понимается количество мигранта, переносимого в единицу времени через единичную площадь общего поперечного сечения по­тока (рассматриваемого с позиций квазиоднородной сплошной среды).

[10] Описание такой модели дано в основном по материалам исследований, проведенных в ГДР. Обзор исследований в этом направлении в СССР и других странах приведен в работе [1].

[11] дР дх2

£ = ^ . (9.5)

Соответственно величину п можно определить по времени ім прохождения максимума концентрации трассера, преобразуя фор­мулу (8.12) применительно к условиям радиального потока:

168

Моделирование миграции подземных вод

Составление теоретических гидрогеомиграционных моделей

Представим в свободном виде общие теоретические модели мигра­ционных процессов, реализующие рассмотренные выше концеп­ции о различных гидрогеомиграционных (геомиграционных) про­цессах на системных (кибернетических) позициях. Базовым для построения гидрогеомиграционных моделей явля­ется балансовое уравнение …

Опыты в фильтрующих образцах

Фильтрующие образцы (нарушенного и ненарушенного строения) имеют сравнительно небольшие размеры (обычно 10—20 см). На таких образцах целесообразно определять параметры обменных лроцессов между твердой и жидкой фазами, с учетом физико-хи­мических превращений …

Биогенные преобразования химического состава

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

Как с нами связаться:

Украина:
г.Александрия
тел. +38 05235 7 41 13 Завод
тел./факс +38 05235  77193 Бухгалтерия
+38 067 561 22 71 — гл. менеджер (продажи всего оборудования)
+38 067 2650755 - продажа всего оборудования
+38 050 457 13 30 — Рашид - продажи всего оборудования
e-mail: msd@inbox.ru
msd@msd.com.ua
Скайп: msd-alexandriya

Схема проезда к производственному офису:
Схема проезда к МСД

Представительство МСД в Киеве: 044 228 67 86
Дистрибьютор в Турции
и странам Закавказья
линий по производству ПСВ,
термоблоков и легких бетонов
ооо "Компания Интер Кор" Тбилиси
+995 32 230 87 83
Теймураз Микадзе
+90 536 322 1424 Турция
info@intercor.co
+995(570) 10 87 83

Оперативная связь

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