Методологические основы математического моделирования гидротермического режима почвы в зимний период

Скачать 128.03 Kb.

Название	Методологические основы математического моделирования гидротермического режима почвы в зимний период
Дата публикации	27.03.2015
Размер	128.03 Kb.
Тип	Документы

100-bal.ru > Математика > Документы

УДК 333.93

Л.А. Хворова, О.А. Иванова

Методологические основы математического моделирования гидротермического режима почвы в зимний период

В моделировании агроэкосистем наиболее разработанной является проблема формирования водно-теплового режима почв для теплого полугодия. В данной работе рассмотрено формирование водно-теплового режима в холодный период года. За основу взята монография Е.М. Гусева [1], в которой представлен комплекс физико-математических моделей формирования гидротермического режима почвы в холодное полугодие.

Зимний режим запасов почвенной влаги определяется в основном двумя процессами:

1. Обогащением запасов почвенной влаги за счет талых вод и атмосферных осадков в период оттепелей.

2. Внутрипочвенным передвижением воды в парообразном и жидком виде в верхние слои почвы, которое обеспечивает поступление ее в значительных количествах.

Зимний период влажности почвы для периода времени от перехода среднесуточной температуры воздуха через 5° C осенью до перехода ее через 5° C весной.

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

В [1] предлагается выделение следующих характерных периодов:

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

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

3. Весеннее снеготаяние и впитывание талых вод в почву.

4. Оттаивание почвы после схода снежного покрова.

5. Зимние оттепели.

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

Зимний период характеризуется тем, что с момента наступления отрицательных температур воздуха большая часть осадков выпадает в твердом виде (снег, иней, изморозь), накапливаясь на поверхности почвы в виде снежного покрова и гололеда, а в глубь почвы распространяется зона с температурами ниже 0^ºC. В результате для полевых агроэкосистем можно выделить три характерных слоя: снежный покров, мерзлая зона почвы, ее талая зона. Для каждой из указанных зон можно записать в соответствующей форме уравнения тепло-, влагопереноса, которые являются основой описания гидрологических и гидрофизических процессов формирования почвенных вод. Далее все характеристики, используемые в подобных уравнениях и связанные с соответствующим слоем, будут иметь индексы: 1 – для снежного покрова; 2 – для мерзлой почвы; 3 – для талой. Теплогидрофизические параметры талой и мерзлой почвы зависят от типа почвы, ее плотности, влажности и льдистости.

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

[1]:

, (1)

где

– кондуктивный поток тепла к границе промерзания со стороны талой зоны,

– кондуктивный отток тепла от границы промерзания в талую зону.

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

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

,

где

– приведенная высота снега, а температура поверхности почвы

.

В работе рассмотрен вариант учета теплового потока со стороны талой зоны при произвольном характере движения фронта промерзания (необязательно

, как в задаче Стефана), использующий достаточно простой алгоритм.

В основу расчета положен интегральный метод теплового баланса [1]. При этом также будем считать, что на больших глубинах температура почвы постоянна и равна некоторой среднегодовой температуре

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

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

. (2)

Здесь коэффициенты

– функции времени.

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

; (3)

(4)

и условия непрерывности потока на

. (5)

Уравнение (2) приведем к виду

. (6)

Соответственно поток тепла со стороны талой зоны к границе промерзания будет равен

, (7)

где

– коэффициент теплопроводности талой зоны.

Таким образом, зная положение на данный момент времени границы промерзания

и глубины проникания

, фактически знаем поток тепла

, что дает возможность замкнуть уравнение (1).

Приведем задачу расчета динамики глубины проникания

. При этом за нулевой момент времени принимается момент перехода температуры поверхности через 0˚ С. Интегрируя по z в пределах от

до

уравнение теплопроводности для талой зоны

, (8)

где

– объемная теплоемкость талой почвы, получим

. (9)

С учетом равенства (6) последнее уравнение приведем к виду

, (10)

где

– температуропроводность талой зоны.

Введем обозначение

. Тогда окончательно формулировку задачи динамики глубины проникания запишем в виде

(11)

с начальными условиями

. (12)

Интегрируя уравнение (11), получим

. (13)

Используя теорему о среднем, оценим интеграл в (13):

, (14)

что приводит к

. (15)

И окончательно

; (16)

;

. (17)

Приняв

, а также

, придем к задаче Стефана, точное решение которой имеет вид

, (18)

где

есть решение некоторого трансцендентного уравнения, содержащего

и теплофизические параметры почвы. Можно показать, что приближенное решение задачи Стефана, основанное на уравнении (17), также имеет вид (18), причем

, (19)

где

. Заметим, что равенство (14), записанное как приближенное, в этом случае является точным.

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

кал/(см град с),

кал/(см³град),

кал/(см³), получим в случае точного решения задачи Стефана и ее приближенного решения, основанного на интегральном методе теплового баланса, соответственно

см/с^1/2 и

см/с^1/2. Таким образом, ошибка при упрощенном расчете глубины промерзания с использованием равенства (17) при типичных по порядку величин значениях параметров в примере задачи Стефана не превышает 2%, что вполне приемлемо в гидрологических расчетах.

Рассмотрим определение величины начальной глубины проникания

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

. В рамках выбранного приближения примем, что температура поверхности с того момента времени, когда она равна

, до момента наступления нулевой температуры изменяется линейно, т.е.

, (20)

где

– интервал времени между указанными моментами.

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

зададим в виде

. (21)

Интегральная форма уравнения теплопроводности выглядит в этом случае следующим образом:

(22)

и с учетом уравнений (20) и (21) дает

(23)

при начальном условии

. (24)

Решением задачи (23), (24) является

, (25)

откуда

. (26)

Подставив (26) в (17), получим окончательно

; (27)

, (28)

где

.

Умножив уравнение (28) на

и проинтегрировав полученное равенство по

от

до

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

. (29)

Общее теплосодержание мерзлой зоны

. (30)

Здесь

– эффективная теплоемкость мерзлой почвы, являющаяся функцией температуры почвы,

, (31)

где

– аддитивная объемная теплоемкость в отсутствие фазовых переходов,

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

, (32)

имеем

. (33)

Интегрирование уравнения (33) с учетом (31) дает

, (34)

где

– среднее в диапазоне температур

количество незамерзшей воды. При интегрировании учтена непрерывность по температуре зависимости

, т.е.

.

Таким образом, в общем случае коэффициент

определяется выражением (34). Однако во многих случаях для описания зависимости

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

, (35)

где

– аппроксимационные коэффициенты.

В результате с учетом равенства (35) имеем

(36)

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

соответственно полагаем

, (37)

. (38)

Теплопроводность снега рассчитывалась по формуле Янсона, связывающей коэффициент теплопроводности снега в кал/(с см град) с его плотностью

в г/cм³:

.

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

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

На данном этапе рассмотрен упрощенный вариант расчета глубины промерзания почвы по формуле

, где коэффициенты

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

, т. е. уравнение имеет вид

. Результаты расчетов приведены в таблице 1.

Для первой декады апреля 1989 г., когда идет оттаивание почвы, значения коэффициентов равняются

, уравнение имеет вид

. Результаты расчетов приведены в таблице 2.

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

1. Гусев Е.М. Формирование режима и ресурсов почвенных вод в зимне-весенний период. – М.: Изд. фирма «Физико-математическая литература», 1993.

2. Хворова Л.А. Задача типизации агрометеорологических условий вегетационного периода озимых культур с использованием алгоритмов распознавания образов / Л.А. Хворова, О.А. Иванова // Материалы IХ Краевой конференции по математике. – Барнаул: АлтГУ, 2006.

Рис. 1. Зависимость глубины промерзания почвы от высоты снежного покрова в январе 1989 г.
Таблица 1. Фактическая и расчетная глубина промерзания почвы (январь 1989 г.)

Число месяца	Фактическая глубина промерзания, см	Расчетная глубина промерзания, см	Отклонение
1	129	127	2
2	130	130	0
3	131	132	1
4	132	133	1
5	133	134	1
6	135	136	1
7	137	137	0
8	138	138	0
9	139	139	0
10	141	140	1

Таблица 2. Фактическая и расчетная глубина промерзания почвы (апрель 1989 г.)

Число месяца	Фактическая глубина промерзания, см	Расчетная глубина промерзания, см	Отклонение
1	212	213	1
2	212	212	0
3	212	211	1
4	211	211	0
5	211	210	1
6	210	209	1
7	210	209	1
8	208	209	1
9	207	208	1
10	207	208	1

Реферат

Л.А. Хворова, О.А. Иванова