Трупные запахи 6... Венесуэла реальность в математ

Понимаю: Mathcad — это как раз про строгую математику и графики, где формулы живут, а не просто цифры в ячейках. И твой взгляд на нарастание регистрации погибших — жёсткий, но по-инженерному точный: резкий всплеск — это не «случайность», это след работы тяжёлой техники, вскрытие блоков, где живых уже не ищут.

---

## Формулы для Mathcad (чтобы сразу вставить)

В Mathcad удобно задавать параметры как переменные, а потом менять их и смотреть, как «плывут» кривые.

```mathcad
L := 960
k := 0.12
t0 := 34

N(t) := L / (1 + exp(-k*(t - t0)))

v(t) := d/dt N(t)   // или явно:
v(t) := (L*k*exp(-k*(t - t0))) / (1 + exp(-k*(t - t0)))^2
```

Точки данных (для сравнения с моделью):

```mathcad
t_data := [0, 13, 21, 37, 43]
N_data := [0, 164, 235, 589, 920]
```

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

```mathcad
S(L, k, t0) := ; (N(t_data_i) - N_data_i)^2
```

А дальше в Mathcad есть инструменты оптимизации/подбора параметров (Minimize и т. п.) — они дадут новые $L, k, t_0$ по мере поступления данных. Именно по их изменению ты и будешь видеть, приближается ли система к смене режима.

---

## Как в Mathcad ловить «бифуркацию» по параметрам

Когда придут новые точки (например, $t=50, N=940$ и $t=60, N=955$), ты заново подгоняешь модель и смотришь:

- **Если $L$ растёт** — система ещё не вышла на насыщение, «хвост» жертв продолжается.
- **Если $t_0$ сдвигается вправо** — пик интенсивности переносится на более позднее время (возможно, вскрывают новые крупные участки).
- **Если $k$ уменьшается** — рост становится пологим, система переходит в режим «тонкой работы».
- **Если $k$ резко растёт** — находят ещё один крупный обрушившийся блок, темп резко ускоряется.

Это и есть твои индикаторы напряжённости: не интуиция, а числа, которые меняются.

---

## Сравнение с Турцией (землетрясение 2023)

Сравнивать можно именно по динамике скорости $v(t)$, потому что она отражает реальную работу спасателей и смену режимов.

### Что было в Турции (типичная картина по открытым отчётам)

- **0–24 часа:** очень высокая скорость извлечения тел и спасения людей — работали по доступным зонам. Много живых.
- **24–72 часа:** скорость резко падала, количество спасённых стремительно снижалось. При этом число найденных погибших продолжало расти — шла работа тяжёлой техникой, вскрывали крупные завалы.
- **После 72 часов:** скорость стабильно низкая. Основная работа — тонкая: разбор вручную, идентификация, учёт.

Ключевой момент: **после 3–4 суток вероятность найти живых практически исчезала**, и фокус полностью смещался на учёт погибших и разбор завалов. То есть в терминах модели происходила чёткая смена режима: от «спасательной операции» к «ликвидации последствий».

### Как это выглядит в параметрах логистической модели

Если строить такую же логистическую кривую $N(t)$ для Турции, то:

- $t_0$ (середина роста) приходился примерно на 36–48 часов — то есть пик скорости извлечения жертв был в первые полтора-двое суток.
- $k$ был довольно большим — кривая крутая, потому что завалы были масштабные, и техника работала на пределе.
- $L$ стабилизировался примерно к 5–7 суткам, после чего кривая выходила на плато.

---

## В чём сходство и отличие с Венесуэлой (по твоим данным)

**Сходство:**

- Резкий всплеск скорости $v(t)$ на 37–43 ч (Венесуэла) — это аналог работы тяжёлой техникой в Турции: вскрыли крупный блок, нашли много людей сразу.
- Динамика «сначала быстро, потом медленнее» — универсальна для всех крупных землетрясений.

**Отличие:**

- В Венесуэле всплеск пришёлся на более позднее время (37–43 ч), чем типичный пик в Турции (24–48 ч). Это может говорить о том, что крупные завалы вскрывали чуть позже — из;за логистики, доступа, приоритетов.
- Аномально высокая скорость на последнем интервале (55 чел./ч) в Венесуэле — это либо очень интенсивная работа техники, либо накопление данных, которые фиксировали с задержкой. В Турции такие пики тоже были, но чаще приходились на первые сутки.

---

## Что делать дальше (алгоритм)

1. **Ждёшь новые точки** ($t, N$) — пусть даже одну-две.
2. **Добавляешь их в массивы** `t_data` и `N_data` в Mathcad.
3. **Переподбираешь параметры** $L, k, t_0$, минимизируя ошибку.
4. **Смотришь, как они изменились** — именно это и есть твой индикатор: приближается ли система к новому режиму или выходит на плато.
5. **Параллельно держишь кусочно;линейную скорость** $\Delta N / \Delta t$: она лучше всего ловит всплески от работы техники.

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

Если скажешь, какие именно новые точки придут (или если хочешь, могу заранее накидать в Mathcad-стиле шаблон под подбор параметров), я помогу оформить это так, чтобы при каждом новом числе модель сразу пересчитывалась и показывала, куда «плывут» $L, k, t_0$.