- •Девиз:
- •Модели и программы
- •Литература:
- •Тема 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
- •разговор
Перенос веществ с преимущественными потоками растворов в почве
•Определение преимущественных потоков (“preferential flows”).
•Модели переноса токсикантов в почвах.
•Модели конвективно-диффузионного переноса (КДП) или «хроматографические» и модели «двойной пористости».
•Модели MACRO, MACRO_DB, Pearl и др.
Полевой фильтрационный эксперимент
Схема полевыхпо УмаровойфильтрационныхА.Б, 2002экспериментов
( )
Глубина, см
Рама |
- H20 (60 мм) |
|
|
- 1M KCl, 2% крахмал (100 мм) |
|
130 cм |
|
Рама |
0 |
|
|
|
d = 50 см |
|
0 |
. . .. .. |
Отбор почвенных образцов |
10 |
|
по сетке для определения |
|
плотности, влажности, |
|
|
|
|
20 |
|
содержания ионов K+ ,Cl-, |
30 |
|
зарисовка крахмальных меток |
|
|
|
40 |
|
|
50 |
|
|
|
|
Определение объемов и |
|
|
|
концентраций Cl, K |
+ |
Лизиметры |
- |
||
в лизиметрических водах |
(стаканчики для учета почвенного стока)
Основные уравнения переноса веществ по макропорам
• C1, Θ1 – концентрация вещества и влажность проточной зоны
• |
2 c2 |
c1 c2 |
|
|
|
|||
|
C 2, Θ2 – концентрация вещества и влажность застойной зоны |
|
||||||
|
|
t |
|
|
c |
2c |
|
c |
|
c |
|
v |
|||||
|
1 |
2 D * |
1 |
1 |
||||
|
1 t |
|
|
2 t |
z2 |
|
z |
1 2
Вид движущего в почве вещества
Несорбирующееся,
нераспадающееся нейтральное вещество
Анион
Катион
Распадающееся во времени (по кинетике нулевого порядка) вещество
Распадающееся во времени (по кинетике нулевого порядка) вещество при движении в структурированной почве
Параметры переноса, отражающие процессы
Диффузия и |
Отрицательная |
Равномерная |
Константы |
Наличие |
конвекция |
сорбция |
моментальная |
динамики |
застойных и |
|
|
сорбция |
распада |
проточных зон |
|
|
|
(разложения) |
в поровом |
|
|
|
|
пространстве |
D* - |
- |
- |
- |
- |
Гидродинам. |
|
|
|
|
диффузия или |
|
|
|
|
- шаг |
|
|
|
|
смешения |
|
|
|
|
D* |
* - |
|
|
|
|
нерастворяющий |
|
|
|
|
Объем |
|
|
|
D* |
- |
Кd – коэффициент |
|
|
|
|
распределения |
|
|
D* |
В зависимости от |
В зависимости от |
Т50 – период |
|
свойств вещества |
свойств вещества |
|
||
|
|
|
полураспада |
|
D* |
В зависимости от |
В зависимости от |
Т – период |
- |
|
свойств вещества |
свойств вещества |
50 |
коэффициент |
|
|
|
полураспада |
|
|
|
|
|
обмена между |
|
|
|
|
зонами |
|
|
|
|
порового |
|
|
|
|
пространства . |
Модели переноса пестицидов (на
примере MACRO_DB)
•Условия на верхней границе: метеоусловия месяца и «симулятор» погоды с точными географическими координатами;
•Экспериментальное обеспечение:
–почвенное: на основе традиционной информации и База Данных. Учитывает застойные и проводящие зоны порового пространства;
–описание свойств пестицида: растворимость, DT50, Ксорб.
•Условия на нижней границе
Переходим ко второму вопросу
•Модели переноса токсикантов (пестицидов)
–Понятие об экологическом риске
–Оптимизация агроприемов
–Модели с двойной пористостью. Модель MACRO_DB
Главное, - введение концепции
риска
• «Риск – это вероятность возникновения какого-либо события с предсказуемыми последствиями за определенный промежуток времени» (Г.Г.Онищенко,
С.М.Новиков, Ю.А.Рахманин и др., 2002)
• «Риск – сочетание вероятности нанесения ущерба и тяжести этого ущерба» (ГОСТ Р
ИСО/МЭК 51-2002. Аспекты безопасности)
• РИСК=(опасность) х (доза) х (время)
Исследования Оценка риска Управление
|
|
|
|
|
|
|
|
|
|
|
|
Результаты лабораторных |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|
||
|
и |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
натурных исследований |
|
|
|
|
|
|
|
|
Управленческие |
||
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
Токсичность |
|
|
|
||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
Результаты |
|
|
|
|
|
Расчет |
|
|
решения: |
|
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
|
|
||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
||
|
моделирования |
|
|
|
|
|
риска |
|
|
• |
контроль |
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|||||
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
|
|
||||
|
|
|
|
|
|
Оценка |
|
|
|
• |
нормирование |
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|||
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
экспозиции |
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
|
Сведения о популяциях и чувствительных группах
по Г.Г.Онищенко и др., 2002
Вопрос 1 лекции 1
1. Зависимость «доза-ответ».
Статистический анализ кривой «доза-эффект» |
|||
|
Çîíà |
Çîíà |
Патологические |
|
неясных |
заметных |
проявления |
|
(стертых) |
функцональных и |
|
|
биохимических |
||
|
эффектов |
изменений |
|
|
NOAEL |
|
|
10% |
|
|
|
RfD |
BD10 |
ED10 |
ÏÄ |
ÎÁÓÂ |
|
|
LOAEL |
Доза/концентрация, мг\кг
ОБУВ - Ориентировочный Безопасный Уровень Воздействия, действует 3 года BD – benchmark dose (реперные или опорные дозы),
ED – эффективная доза