
- •Девиз:
- •Модели и программы
- •Литература:
- •Тема 1: «Экспериментальное обеспечение моделей»
- •Мы изучаем модели:
- •Физическое обоснование. Уравнение неразрывности
- •Вывод уравнения неразрывности
- •Расчетная схема
- •Основные уравнения
- •Основное уравнение переноса в конечно-разностной форме
- •Сеточная схема расчета (используется при численном моделировании методом конечных элементов, -метод Галеркина)
- •Таким образом, экспериментальным обеспечением моделей влагопереноса
- •Начальные условия
- •Модель движения воды в почве
- •Граничные условия
- •Продолжение темы 1
- •Нижнее граничное условие
- •Это было о потоках влаги, о динамике влажности почвы
- •Пункт 2 темы: Перенос тепла в почве qТ dTdz
- •Уравнение теплопереноса
- •Основное уравнение теплопереноса
- •Расчет переноса тепла внутри почвы
- •Уравнение теплопереноса в
- •Условия для применения расчетной схемы
- •Расчет переноса тепла внутри почвы
- •Уравнение теплопереноса в
- •Условия для применения расчетной схемы
- •Условия для применения расчетной схемы
- •стоп
- •Педотрансферные функции
- •«Педотрансферные функции превращают данные, которые у нас есть, в данные, которые нам нужны!»
- •Определяют педотрансферные функции
- •1-й способ (точечный)
- •2-й способ (параметрический)
- •Педотрасферные функции
- •Нейросетевые алгоритмы
- •Математическая модель нейрона
- •«Педотрансферные функции превращают данные, которые у нас есть, в данные, которые нам нужны!»
- •Тема 2. Аппроксимация эмпирических зависимостей
- •Виды функций
- •Убывающие функции
- •С одним максимумом
- •V.С несколькими экстремумами
- •Определение
- •Следующий вопрос темы 2
- •Вхимии почв: уравнения сорбции
- •Физический смысл параметров
- •3.Определение параметров аппроксимации
- •Задача – найти минимум S, подобрав соответствующие значения b1 и b2
- •Определение
- •Расчет среднеквадратической
- •Метод сканирования
- •СТАТИСТИКИ ДЛЯ АНАЛИЗА ПАРАМЕТРОВ
- •Анализ ошибок (погрешностей)
- •Практическое задание 1.
- •Зависимость Y от X
- •«Выдача» после аппроксимации
- •Гистограмма погрешностей
- •График реальных значений от расчетных
- •График зависимости погрешностей от расчетной величины
- •Анализ на систематические погрешности
- •Сравнение параметров
- •Статистическое доказательство значимости различий параметров
- •Физический смысл сравнения параметров (пример с термофильными микроорганизмами)
- •Правила подбора функций
- •Практическое занятие №1
- •Модели и программы
- •Тема 2. Аппроксимация эмпирических зависимостей
- •Виды функций
- •Убывающие функции
- •С одним максимумом
- •V.С несколькими экстремумами
- •Определение
- •Следующий вопрос темы 2
- •Вхимии почв: уравнения сорбции
- •Физический смысл параметров
- •3.Определение параметров аппроксимации
- •Задача – найти минимум S, подобрав соответствующие значения b1 и b2
- •Определение
- •Расчет среднеквадратической
- •Метод сканирования
- •СТАТИСТИКИ ДЛЯ АНАЛИЗА ПАРАМЕТРОВ
- •Анализ ошибок (погрешностей)
- •Практическое задание 1.
- •Зависимость Y от X
- •«Выдача» после аппроксимации
- •Гистограмма погрешностей
- •График реальных значений от расчетных
- •График зависимости погрешностей от расчетной величины
- •Анализ на систематические погрешности
- •Сравнение параметров
- •Статистическое доказательство значимости различий параметров
- •Физический смысл сравнения параметров (пример с термофильными микроорганизмами)
- •Правила подбора функций
- •Практическое занятие №1
- •Задание после аппроксимации
- •Заключительные вопросы (!)
- ••ЗАДАНИЕ ПО ПАРАМЕТРИЗАЦИИ ЭКСПЕРИМЕНТАЛЬНЫХ ЗАВИСИМОСТЕЙ (АППРОКСИМАЦИЯ ЭКСПЕРИМЕНТАЛЬНЫХ ДАННЫХ)
- •Получить параметры аппроксимации для 2-й экспериментальной зависимости. Сравнить достоверность отличий полученных параметров по
- •Пример работы с программой STATISTICA
- •Аппроксимация в программе STATISTICA
- •Сравнение функций для описания экспериментальных данных и выбор лучшей
- •Сравнение моделей. Критерий Вильямса-Клюта
- •Из раздела: Аппроксимация
- •Пример расчета критерия Вильямса-Клюта
- •4. Функции, наиболее употребительные в почвоведении
- •Вхимии почв: уравнения сорбции
- •Физический смысл параметров
- •А если мы не знаем вида? Или функция очень сложная? С несколькими экстремумами?
- •Пример Харитоновой
- •4. Функции, наиболее употребительные в почвоведении
- •Вхимии почв: уравнения сорбции
- •Физический смысл параметров
- •В прошлый раз… Заканчиваем тему «Аппроксимация экспериментальных данных»
- •Впрограмме STATISTICA
- •Сравнение функций для описания экспериментальных данных и выбор лучшей
- •Сравнение моделей. Критерий Вильямса-Клюта
- •Для одних и тех же данных можно
- •• Критери й Вильямса-Клута (по Рыбалко
- •Пример расчета критерия Вильямса-Клюта
- •4. Функции, наиболее употребительные в почвоведении
- •Вхимии почв: уравнения сорбции
- •В биологии почв, экологии: физический смысл параметров гауссиады
- •Раздел почвоведения, использвание
- •А если мы не знаем вида? Или функция очень сложная? С несколькими экстремумами?
- •Общие правила подбора функций для экспериментальных данных
- •Новая тема
- •Работа с физически обоснованными моделями Модель водо- и солеобмена почв HYDRUS
- •Структура моделей
- •Типы наиболее употребительных моделей
- •2 типа моделей:
- •Составляющие водного баланса почв
- •Водный баланс почвы (одномерный)
- •Балансовые модели основаны на использовании гидрологических
- •Устройство балансовых моделей
- •Имитационные модели ……
- •HYDRUS
- •Что будем считать?
- •«Первая задача с прогнозной моделью HYDRUS»
- •Основные расчетные процессы модели HYDRUS
- •Geometry information
- •Заказ
- •Time information
- •Формирование препроцессора
- •Условия сеточной схемы
- •Параметры ван Генухтена и педотрансферные функции
- •NB! Запомнить, записать.
- •Граничные условия
- •Метеоусловия (условия на
- •Графический редактор Слои (материал 1 и 2)
- •Графический редактор (начальные условия: давление влаги по профилю)
- •Постпроцессор
- •Об умении читать выдачи
- •Пример работы с программой STATISTICA
- •Тема 3 «МОДЕЛИ РАВНОВЕСИЯ И ПЕРЕНОСА ВЕЩЕСТВ В ПОЧВАХ»
- •Ионные равновесия в растворах (гомогенный обмен)
- •Основные уравнения
- •Ионные равновесия при наличии твердой фазы (гетерогенный обмен)
- •3. Модели переноса веществ
- •Конвективный переносC веществ
- •Движение в почве с учетом
- •Явление гидродинамической дисперсии (по Holzbecher,1998)
- •Гидродинамическая дисперсия и
- •Конвективно-диффузионное уравнение
- •Определения
- •Конвективно-диффузионное уравнение
- •Это конвективно-дифф.уравнение, но без учета процессов, которые могут приводить к появлению дополнительных количеств
- •Конвективно-диффузионное уравнение
- •1.Ионный обмен
- •Итак, о слагаемом в конвективно-
- •Нерастворяющий объем (объем выталкивания, exclusive water)
- •1.Ионный обмен
- •Сорбция/десорбция кинетическая
- •Кинетика нулевого порядка
- •Кинетическая сорбция (разложение) веществ по кинетике 1-го порядка
- •Кинетики десорбции/разложения
- •Аналогично сорбции/десорбции описываются и процессы разложения, распада, роста (микроорганизмов) и пр.
- •Мы разобрались с описаниями всех процессов,
- •Основные параметры переноса солей
- •ВЫВОД:
- •Перенос веществ с преимущественными потоками растворов в почве
- •Полевой фильтрационный эксперимент
- •Основные уравнения переноса веществ по макропорам
- •Вид движущего в почве вещества
- •Модели переноса пестицидов (на
- •Переходим ко второму вопросу
- •Главное, - введение концепции
- •Исследования Оценка риска Управление
- •Вопрос 1 лекции 1
- •Статистический анализ кривой «доза-эффект»
- •Установление реперной концентрации (ВС5) для формальдегида
- •ПДК, ОБУВ и риски. Определения
- •Пример с агротехнологиями
- •ВЫВОД:
- •Полевой фильтрационный эксперимент
- •Основные уравнения переноса веществ по макропорам
- •Вид движущего в почве вещества
- •Модели переноса пестицидов (на
- •Конвективно-диффузионное уравнение
- •1.Ионный обмен
- •1.Ионный обмен
- •Сорбция/десорбция кинетическая
- •Кинетика нулевого порядка
- •Кинетическая сорбция (разложение) веществ по кинетике 1-го порядка
- •Кинетики десорбции/разложения
- •Аналогично сорбции/десорбции описываются и процессы разложения, распада, роста (микроорганизмов) и пр.
- •Основные параметры переноса солей
- •ВЫВОД:
- •Перенос веществ с преимущественными потоками растворов в почве
- •Полевой фильтрационный эксперимент
- •Основные уравнения переноса веществ по макропорам
- •Вид движущего в почве вещества
- •Модели переноса пестицидов (на
- •Переходим ко второму вопросу
- •Главное, - введение концепции
- •Исследования Оценка риска Управление
- •Вопрос 1 лекции 1
- •Статистический анализ кривой «доза-эффект»
- •Установление реперной концентрации (ВС5) для формальдегида
- •Вид движущего в почве вещества
- •Модели переноса пестицидов (на
- •Блок метеоусловий
- •Симулятор погоды
- •Сравнение моделей. Критерий Сайерта
- •Основные параметры модели
- •Схема передвижения веществ к корню растения
- •Основные уравнения модели
- •Определение чувствительности модели
- •Анализ чувствительности модели
- •Анализ чувствительности модели к параметрам ван Генухтена
- •Оценка
- •Тема IV
- •разговор

Вхимии почв: уравнения сорбции
•Степенная функция
Q mCn |
Уравнение Фрейндлиха |
• Логистическая
A A |
|
KС |
|
|
л |
Уравнение Ленгмюра |
|
|
KС |
||
1 |
Аналогичные уравнения Моно, |
||
|
|
л |
Михаэлиса-Ментен и др. |
|
|
|
•Гауссиада (агрохимия, биология почв, статистика и др.)

В биологии почв, экологии: физический смысл параметров гауссиады
• В уравнении гауссиады
|
|
|
|
x b2 2 |
|
|
|
b1 |
exp |
|
|
|
|
|
|
b3 |
|||||||
|
|
|
|
|
|
||
|
|
|
|
|
b1 – обилие вида
b2 – биологический оптимум
b3 – толерантность (мера экологической амплитуды)

Раздел почвоведения, использвание
Накопление почвенного органического вещества.
Почвоведение Основная гидрофизическая
характеристика.
Почвенная гидрология Основная гидрофизическая
характеристика.
Почвенная гидрология Моментальная адсорбция
Химия почв Моментальная адсорбция
Химия почв Поглощение веществ растением.
Агрохимия Кинетика 1-го порядка.
Химия, физика почв Кинетика 2-го порядка.
Кинетика нулевого порядка
Логарифмическая зависимость урожая от химических веществ (удобрений)
Функции в почвоведении |
|||||||||||||||||||||
Функция |
Вид |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Автор |
|
Название |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Уравнение Костычева-Иенни по |
||
Показательная |
|
|
y t |
y0 |
1 |
|
e |
kt |
|||||||||||||
|
|
|
|
|
накоплению органического |
||||||||||||||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
1 |
вещества почв |
|
|
|
|
|
|
|
|
|
|
|
n |
|
|
|
|
|
|
||||
Степенная |
|
i r |
|
|
PдляР |
|
|
|
|
|
|
Brooks and Corey, 1964 |
|||||||||
|
|
|
|
к c |
|
|
|
|
|
|
|
|
к с |
|
|||||||
|
Se |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
s r |
|
|
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|||||
|
|
|
1дляР |
|
|
|
|
|
|
|
|||||||||||
|
|
|
|
|
|
|
|
|
|
||||||||||||
|
|
|
|
|
|
|
|
|
|
к с |
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
Логистическая |
|
|
|
|
|
|
|
|
|
1 |
|
|
|
|
|
|
|
m |
|
|
ван Генухтеном (van Genuchten, |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
дляPк c 1980 |
|||||
|
|
i |
r |
|
|
|
|
|
|
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
n |
||||||||||
|
Se |
|
|
|
|
1 |
( P |
|
|
) |
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
|
|
|
|||||||||||
|
|
s |
r |
|
|
|
|
|
к с |
|
|
|
|
|
|
||||||
|
|
|
|
|
|
|
|
|
|
|
0 |
|
|
|
|||||||
|
|
|
|
|
|
дляР |
|
|
|
|
|
||||||||||
Степенная |
|
|
|
|
|
sк с |
|
|
|
|
|
|
|
|
|
Уравнение Фрейндлиха |
Логит функция
Логит функция
Экспоненциальн
ая Логит функция
Линейная
Q mCn
A A KС
л
1 KС
л
Jr Jmax Kmc cmin
1 Kmc
Ci C0 exp( k1t)
Ci |
|
|
C0 |
|
|
|
1 |
k2 C0 |
t |
||||
|
Уравнение Ленгмюра.
Уравнение Михаэлиса-Ментен
Логарифмичес
кая
Ci C0 k0t
гдеlgА –Aмаксимальныйy lg A b(урожай,x c) у
– урожай при внесении того или иного химического вещества в почву в количестве х, и b – фактор действия этого вещества,
Уравнение Макса Айдхарда Альфреда Мичерлиха

А если мы не знаем вида? Или функция очень сложная? С несколькими экстремумами?
•Например, кривая сорбции паров воды.
•Это S-образная функция. У нас не было функций такого вида. Что делать в этом случае?
•Надо использовать принцип пошаговой регрессии, начиная с полинома.

Общие правила подбора функций для экспериментальных данных
• Определите вид функции (выбор из VШ типов:
возрастающие, убывающие, с одним максимумом, минимумом….).
•Задайте функцию с наименьшим числом параметров по правилу: n=b+3, где n -число параметров, b – число пар «функция- аргумент».
•Испытайте другой тип функции и по критерию Вильямса-Клюта докажите различия и выберите лучшую.
•Попытайтесь выяснить и оценить физическое значение параметров.

Новая тема
«Использование физически обоснованных моделей
•Балансовые и имитационные (с сеточной схемой расчета) модели;
•Препроцессор и постпроцессор;
•Формулировка задачи; пространственные и временные границы моделируемого процесса;
•Задание почвенных свойств и характеристик. Графический редактор;
•Постпроцессор. Выдачи модели.

Работа с физически обоснованными моделями Модель водо- и солеобмена почв HYDRUS

Структура моделей
•Препроцесс
ор
•Пространственные
данные
•Временные
•Задание выдач
•Почвенные данные
•ОГХ
•Кф - фильтрация
•Метео-
•Пост
проце
ссор
•Выдачи

Типы наиболее употребительных моделей
Физически обоснованные модели
|
|
Имитационные |
|
балансовые |
|||
|
(с сеточной схемой расчета) |
||
|
|

2 типа моделей:
балансовые и имитационные (механистические, «сеточные» модели, физически обоснованные).
1.Устройство моделей. Их характерные признаки.
2.Препроцессор моделей
3.Недостатки каждого из типов.