Моделирование процессов тепло- и массопереноса при закачке радиоактивных растворов в глубокозалегающие пласты (стр. 9 из 26)

1.5.4. Математическая постановка задачи массообмена в первом приближении

Уравнения (1.5.31), (1.5.32) для коэффициентов первого приближения принимают вид

(1.5.58)

(1.5.59)

Коэффициенты при

в уравнении (1.5.33) дают

(1.5.60)

Начальные, граничные условия и условия сопряжения

(1.5.61)

(1.5.62)

(1.5.63)

(1.5.64)

Причем, решение

отыскивается в форме квадратного многочлена относительно z(1.5.45), где

задаются выражениями (1.5.48) и (1.5.49), а

неизвестно. Для его определения перепишем (1.5.60) в виде

(1.5.65)

где оператор

(1.5.66)

введён для более компактной записи получающихся соотношений и удобства преобразований. Отметим, что из (1.5.42) следует

(1.5.67)

Учитывая (1.5.45), (1.5.65), а также линейность оператора

, получим

(1.5.68)

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

(1.5.69)

используя которое определим выражения для следов производных на границах сопряжения (1.5.63) через вспомогательные функции, не зависящие от вертикальной координаты z

(1.5.70)

(1.5.71)

Умножая левую и правую части (1.5.71) на

и вычитая полученное из (1.5.70), приходим к уравнению для определения функции

, входящей в квадратичное представление первого коэффициента разложения

(1.5.72)

Уравнение для определения первого коэффициента разложения получается путем подстановки (1.5.68), (1.5.72), (1.5.48), (1.5.49) в (1.5.60)

(1.5.73)

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

(1.5.74)

(1.5.75)

Начальные условия, условия сопряжения и граничные условия

(1.5.76)

(1.5.77)

(1.5.78)

(1.5.79)

Уравнения (1.5.73) – (1.5.79) представляют собой математическую постановку задачи массопереноса для коэффициентов первого приближения.

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

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

1.5.5. Дополнительное интегральное условие для первого приближения

Усредним равенство (1.5.15) по z в пределах несущего пласта согласно

(1.5.80)

Последовательно для каждого слагаемого

(1.5.81)

(1.5.82)

(1.5.83)

Окончательно, после усреднения, получим следующую постановку задачи осреднённого по несущему пласту поля плотностей загрязнителя

(1.5.84)

(1.5.85)

(1.5.86)

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

(1.5.87)

(1.5.88)

(1.5.89)

(1.5.90)

Полученная задача совпадает с задачей (1.5.51) – (1.5.57) для нулевого приближения плотности загрязнителя. В силу единственности решения следует, что