151324 (594679), страница 4
Текст из файла (страница 4)
| (1.4.31) |
| (1.4.32) |
| (1.4.33) |
| (1.4.34) |
При этом плотность загрязнителя, входящая в (1.4.27) – (1.4.29), также будет разлагаться по параметру асимптотического разложения , причём это разложение производится независимо от разложения (1.4.26), хотя и по тому же принципу.
-
Математическая постановка задачи теплопереноса в нулевом приближении
Из (1.4.29) для коэффициентов при (нулевое приближение) получим
, тогда
. Таким образом, в нулевом приближении температура загрязнителя является функцией только от r и t. Из условий сопряжения (1.4.30)
. Следовательно, температура загрязнителя в каждом вертикальном сечении одинакова по всей высоте несущего пласта
. Приравнивая коэффициенты при
к нулю в уравнении (1.4.29), получим
| (1.4.35) |
Сумму первых двух слагаемых в правой части этого уравнения, не зависящую от z, обозначим через
| (1.4.36) |
Тогда
| (1.4.37) |
следовательно,
| (1.4.38) |
При z = 1, воспользовавшись (1.4.30)
| (1.4.39) |
при z = – 1
| (1.4.40) |
Вычитая и складывая два последних уравнения, получим для функций и
следующие выражения:
| (1.4.41) |
| (1.4.42) |
Проинтегрировав (1.4.38), получим
| (1.4.43) |
здесь функция, не зависящая от z, значение которой предстоит найти.
Подставив выражение из (1.4.41) в (1.4.36), получим для нулевого приближения уравнение гиперболического типа со следами производных из внешних областей
| (1.4.44) |
Окончательная постановка задачи в нулевом приближении наряду с (1.4.44) включает также уравнения для окружающих сред, начальные, граничные условия и условия сопряжения
| (1.4.45) |
| (1.4.46) |
| (1.4.47) |
| (1.4.48) |
| (1.4.49) |
| (1.4.50) |
Последнее слагаемое в правой части уравнения (1.4.44) устанавливает изменение температуры за счёт энергии, выделяющейся при радиоактивном распаде. Отметим, что температурное поле в нулевом приближении определяется не значениями плотностей радиоактивного загрязнителя в точках, а усреднёнными значениями по вертикальной координате в интервале пласта. Как будет показано ниже, усреднённое таким образом значение плотности совпадает с нулевым приближением соответствующей задачи массопереноса (см. пункт 1.5.3).
Для определения в нулевом приближении поля температур в среде, как следует из (1.4.44) – (1.4.50), необходимо задание функции плотности радиоактивного загрязнителя. Постановка этой задачи осуществлена в пункте 1.5, а её решению посвящена глава II.
-
Постановка задачи теплопереноса в первом приближении
Уравнения (1.4.27), (1.4.28) для коэффициентов при (первое приближение) принимают вид
| (1.4.51) |
| (1.4.52) |
Для коэффициентов при в (1.4.29)
| (1.4.53) |
Условия сопряжения, начальные и граничные условия
| (1.4.54) |
| (1.4.55) |
| (1.4.56) |
| (1.4.57) |
| (1.4.58) |
Решение отыскивается в виде квадратного многочлена относительно z (1.4.43), где
и
определяются как (1.4.41), (1.4.42), а значение
предстоит найти.
Уравнения (1.4.51) – (1.4.58) определяют постановку задачи теплообмена в первом приближении. Здесь также зависит от плотности загрязнителя, что обусловливается выражениями для
,
.
1.5. Задача массопереноса
1.5.1. Математическая постановка задачи массопереноса и её обезразмеривание
Геометрия задачи массопереноса практически ничем не отличается от температурной задачи и представлена на рис. 1.2.
Рис. 1.2. Геометрия задачи массопереноса |
Математическая постановка задачи массопереноса для всех областей включает уравнение диффузии с учётом радиоактивного распада в покрывающем
| (1.5.1) |
и подстилающем
| (1.5.2) |
пластах, а также уравнение конвективной диффузии с учётом радиоактивного распада в пористом пласте
| (1.5.3) |
При этом граничные условия включают в себя равенства плотностей и потоков растворённого вещества на границах раздела пластов
| (1.5.4) |
| (1.5.5) |
Плотность загрязнителя в скважине, радиус которой мы считаем малым по сравнению с расстояниями до точки наблюдения, равна , т.е.
| (1.5.6) |
В начальный момент времени полагаем плотность загрязнителя равной нулю
| (1.5.7) |
Кроме того, на бесконечности выполняются условия регулярности
| (1.5.8) |
Перейдём к безразмерным координатам (1.4.8). При этом получим следующую постановку задачи: для покрывающего пласта
| (1.5.9) |
для пористого пласта
| (1.5.10) |
для подстилающего пласта
| (1.5.11) |
При этом во втором слагаемом в левой части уравнения (1.5.9) появляется отношение коэффициента диффузии к коэффициенту температуропроводности
| (1.5.12) |
величина которого оказывается порядка ÷
.
Вновь, как и в задаче теплопереноса, последнее слагаемое в левой части уравнения (1.5.10) содержит сомножитель Рd который при существующих объёмах закачки имеет порядок 102, так что конвективная составляющая (вдоль координаты r) для поля концентраций оказывается много значимей, чем диффузионная составляющая. Поэтому в уравнениях (1.5.9) – (1.5.11) пренебрежём молекулярной диффузией вдоль оси r.
Вводя обозначения
| (1.5.13) |
выпишем окончательно интересующие нас уравнения:
| (1.5.14) |
| (1.5.15) |
| (1.5.16) |
Условия сопряжения, граничные и начальные условия при этом принимают вид
| (1.5.17) |
| (1.5.18) |
| (1.5.19) |
| (1.5.20) |
| (1.5.21) |
Уравнения (1.5.14) – (1.5.21) определяет математическую постановку задачи массопереноса.
1.5.2. Разложение задачи массопереноса по асимптотическому параметру
Рассмотрим более общую задачу, получающуюся введением произвольного асимптотического параметра путём формальной замены коэффициента диффузии
на частное
. В соответствии с принятыми обозначениями это отвечает следующим заменам:
,
. Задача (1.5.14) – (1.5.16) становится, таким образом, частным случаем (при
) более общей задачи, содержащей параметр асимптотического разложения как в уравнении для пласта, так и в условиях сопряжения:
| (1.5.22) |
| (1.5.23) |