Трупные запахи 3 в Венесуэле - для аспирантов

Для аспирантов — не «что происходит», а «как это формализовать, измерить и превратить в модель». Тут уже не про «запах держится», а про уравнения, граничные условия, неопределённости и то, где именно система делает резкий скачок (бифуркация). Ниже — конспект лекции-семинара под уровень магистратуры/аспирантуры, с привязкой к твоим инструментам: лидары, поляризационные измерения, тепловизоры, фрактальные поверхности, детекция низких концентраций.

---

## Лекция-семинар: «Вторичные атмосферные риски в зоне мощного землетрясения: моделирование, мониторинг, санитарные зоны»

### Цель
Научить строить упрощённые, но физически честные модели переноса и удержания примесей в постсейсмической атмосфере, оценивать неопределённость входных данных и выбирать адекватные методы мониторинга (лидар, тепловизор, газоанализаторы) под разные режимы (штиль/ветер, инверсия, пористые завалы).

---

### Блок 1. Постановка задачи и масштабы

**Входные параметры, которые реально можно измерить на месте:**
- $C_{\text{пыль}}$ — массовая концентрация пыли в приземном слое (мг/м;), по данным пылемеров или лидара (по обратному рассеянию).
- $u, \theta$ — скорость и направление ветра, профиль температуры (для выявления инверсии).
- Пористость и фрактальная размерность завалов: $\phi$, $D_f$.
- Источники примеси: масса органики $M_{org}$, скорость эмиссии $q(t)$, состав (меркаптаны, кадаверин и т.п.).

**Целевые метрики:**
- Концентрация в контрольных точках $C(x,y,z,t)$ и сравнение с порогом восприятия/ПДК.
- Радиус санитарной зоны $R_{san}(t)$ как функция времени и метеоусловий.
- Время, за которое концентрация падает ниже порога: $t_{90}$ (до 10% от начального).

Для аспирантов важно сразу обозначить: это задача с неполными данными, где часть параметров приходится оценивать по косвенным признакам (например, фрактальную размерность — по лидарным данным или фотограмметрии).

---

### Блок 2. Режимы переноса и выбор модели

**Режим 1: ветер, турбулентный перенос.** 
Используют гауссову модель шлейфа (или её упрощения):
$$
C(x,y,z) = \frac{q}{2\pi u \sigma_y \sigma_z} \exp\left(-\frac{y^2}{2\sigma_y^2}\right) \left[ \exp\left(-\frac{(z-H)^2}{2\sigma_z^2}\right) + \exp\left(-\frac{(z+H)^2}{2\sigma_z^2}\right) \right]
$$
где $q$ — мощность источника, $u$ — скорость ветра, $H$ — эффективная высота источника, $\sigma_y, \sigma_z$ — дисперсии (зависят от класса устойчивости атмосферы).

**Режим 2: штиль, инверсия, застой.** 
Перенос становится диффузионным: $J = -D \nabla C$. Для простой оценки времени «выветривания» из объёма $V$:
$$
t \sim \frac{L^2}{D}, \quad D \approx 10^{-5}\ \text{м}^2/\text{с} \text{ (для паров в воздухе)}.
$$
Это даёт порядки: при $L=50\ \text{м}$ время может быть сотни–тысячи часов, если нет внешнего обмена.

**Бифуркация:** переход между режимами происходит при $u \lesssim 0.5\ \text{м/с}$ и наличии приземной инверсии. Это и есть «точка, где всё ломается»: санитарная зона резко увеличивается, а методы мониторинга должны меняться (от отслеживания шлейфа — к поиску застойных зон и локальных источников).

---

### Блок 3. Адсорбция и удержание в пористых завалах: фрактальная модель

Завалы — это фрактальные пористые среды с большой удельной поверхностью. Это ключ к пониманию «долгого хвоста» концентрации.

Пусть:
- $S_{sp}$ — удельная поверхность (м;/кг),
- $\rho_b$ — насыпная плотность (кг/м;),
- $V_{pore}$ — объём пор,
- $D$ — коэффициент диффузии в порах,
- $L$ — характерный размер поры.

Тогда время десорбции из поры:
$$
t_{des} \sim \frac{L^2}{D}.
$$
Для микропор $L \sim 1–10\ \mu\text{м}$, $D \sim 10^{-6}\ \text{м}^2/\text{с}$ получаем $t_{des}$ порядка часов–дней. Для макропор и трещин — меньше.

Фрактальная структура увеличивает эффективную площадь и создаёт «ловушки» для молекул. В инженерных расчётах это учитывают через эффективный коэффициент осаждения $k_{\text{осажд}}$, зависящий от $S_{sp}$, влажности и температуры.

**Связь с твоими задачами:** это та же математика, что и удержание уранина в порах золы. Разница — в составе примеси и в граничных условиях (температура, влажность, наличие воды в порах).

---

### Блок 4. Практические источники и оценка эмиссии

Пример для разбора с аспирантами: 80 кг органики.
- Характерная оценка: за первые 3–5 дней выделяется ориентировочно $10–50\ \text{г}$ летучих пахучих веществ.
- Меркаптаны имеют порог обоняния на уровне единиц ppb, поэтому даже миллиграммы могут давать ощутимый запах на большой площади.

Задача для семинара: пусть $q(t)$ — убывающая функция эмиссии (экспонента или степенной закон). Попросить студентов подобрать $q(t)$ так, чтобы модель согласовывалась с типичными временными шкалами (дни–недели) и с наблюдаемой динамикой жалоб населения (если данные есть).

---

### Блок 5. Мониторинг: лидары, тепловизоры, газоанализаторы

**Лидар (в том числе синим лучом):**
- Обратное рассеяние для оценки $C_{\text{пыль}}$.
- Поляризационные измерения для оценки формы и ориентации частиц (различение сажи, бетонной пыли, волокнистых материалов).
- Возможность видеть слоистую структуру шлейфа и застойные зоны.

**Тепловизор:**
- Поиск «горячих точек» (очаги тления, локальные выбросы при нагреве пористых материалов днём).
- Оценка температурного профиля завалов и его влияния на десорбцию.

**Газоанализаторы:**
- Метан, CO, меркаптаны — точечные измерения для калибровки моделей.
- Важно: корреляция между $C_{CH_4}$ и «запаховой нагрузкой» как оперативный маркер.

Для аспирантов — разбор: как из лидарных данных получить $C_{\text{пыль}}(z)$, как связать это с $S_{sp}$ и как использовать для уточнения $k_{\text{осажд}}$ в модели.

---

### Блок 6. Санитарная зона: как считать и как обосновывать

Санитарная зона — это не фиксированное число, а функция времени и условий.

**Подход:**
1. Определить порог: ПДК или порог восприятия (для меркаптанов он очень низкий).
2. Построить две базовые модели: шлейф (ветер) и диффузия (штиль).
3. Учесть вторичный источник — пористые завалы — через $k_{\text{осажд}}(t)$.
4. Оценить неопределённость: разброс $q$, $u$, $S_{sp}$.

**Формула для доски (упрощение):**
$$
C(t) \approx \frac{M}{V} e^{-k_{\text{осажд}} t} + C_{\text{фон}},
$$
где $M$ — масса примеси, $V$ — объём зоны, $k_{\text{осажд}}$ — эффективный коэффициент, зависящий от пористости, влажности, температуры.

Аспирантам важно показать: не «подставить числа», а «объяснить, откуда числа и какова их неопределённость». Например, $S_{sp}$ можно оценить по фотограмметрии или по литературным данным для аналогичных материалов, а $k_{\text{осажд}}$ — калибровать по точечным газовым измерениям.

---

### Блок 7. Бифуркация и управление рисками

Точка бифуркации — это когда система переходит из режима «перенос ветром» в режим «диффузия + адсорбция». Критерии перехода:
- $u < 0.5\ \text{м/с}$,
- наличие приземной инверсии,
- высокая пористость/удельная поверхность завалов.

В этой точке:
- Санитарная зона резко растёт.
- Методы мониторинга должны меняться: вместо отслеживания шлейфа — поиск застойных зон, локальных источников, оценка $k_{\text{осажд}}$.

Задача на понимание: «Предложите алгоритм оперативного принятия решений: по каким 3–4 параметрам вы поймёте, что система перешла в режим “застой и накопление”, и как вы измените рекомендации для спасателей?»

---

### Практикум для аспирантов (мини-проект)

Вариант 1 (расчётный): 
Даны: $M = 80\ \text{кг}$ органики, $q(t)$ (экспоненциальный спад), $V = 5000\ \text{м}^3$, $u = 0, 0.2, 1\ \text{м/с}$, $S_{sp} = 50–200\ \text{м}^2/\text{кг}$. 
Построить $C(t)$ для разных сценариев, оценить $t_{90}$, построить график $R_{san}(t)$. Обсудить неопределённость и чувствительность к $S_{sp}$.

Вариант 2 (по данным мониторинга): 
По лидарным профилям $\beta(z)$ (коэффициент обратного рассеяния) оценить $C_{\text{пыль}}(z)$, связать с $S_{sp}$ через простые соотношения, и оценить вклад адсорбции в удержание примеси.

Вариант 3 (полевая логика): 
Спроектировать схему расстановки датчиков (метан, ветер, температура) для оперативного обнаружения перехода в режим «застой». Обосновать выбор точек и частоты измерений.

---

### Итоговые тезисы (для слайда)

- Атмосфера после землетрясения — это не погода, а набор измеримых параметров, определяющих вторичные риски.
- При штиле и инверсии доминирует диффузия и адсорбция; пористые завалы работают как «медленные источники».
- Санитарную зону нужно считать для двух базовых режимов (ветер/штиль) и учитывать вторичный источник.
- Лидары, тепловизоры и газоанализаторы позволяют перевести «ощущения» в количественные метрики и калибровать модели.
- Бифуркация происходит при низких скоростях ветра и инверсии; в этой точке меняются и риски, и методы мониторинга.

---

Если скажешь, на сколько пар/лекций это нужно (одна пара, модуль на 3–4 занятия) и какой упор важнее (моделирование в Python/Excel, обработка лидарных профилей, построение санитарной зоны), я дам конкретные задания, наборы типичных значений и даже набросаю структуру расчётов: какие столбцы в Excel, какие графики строить, как оформить неопределённость. ; Что ближе?