Mathematical modeling of the corpse cooling under conditions of varying ambient temperature

Full Text

ОБОСНОВАНИЕ

Основным условием корректности определения дав­ности наступления смерти (ДНС) с помощью существующих математических моделей охлаждения трупа является постоянство температуры внешней среды [1]. Однако на практике данное требование выполняется лишь в случаях наступления смерти в закрытых помещениях. Вне помещений температура воздуха всегда подвергается суточным колебаниям, отражающимся на процессе охлаждения мёртвого тела. Использование известных математических моделей в подобных условиях будет сопровождаться дополнительными неконтролируемыми и весьма существенными ошибками определения ДНС [1, 2].

Фактически единственным методом математического моделирования охлаждения трупа в условиях изменяющейся температуры окружающей среды является базирующийся на двойной экспоненциальной модели Marshall–Hoare весьма трудоёмкий метод пошагового моделирования, в соответствии с которым посмертный период разбивается на серию коротких отрезков, в каждом из которых температура окружающей среды условно считается постоянной [3]. Решений, обеспечивающих простое непосредственное аналитическое термометрическое определение ДНС в условиях изменяющейся температуры окружающей среды, пока не предложено. Между тем математическое моделирование охлаждения трупа в данных условиях вполне возможно на основе закона охлаждения Ньютона–Рихмана, хорошо подходящего для описания остывания поверхностных тканей трупа [4]. В условиях постоянной температуры окружающей среды применительно к охлаждению трупа данный физический закон имеет вид дифференциального уравнения

$\frac{dT}{dt}=-k(T-T_a)$,                        (1)

где Т — температура диагностической зоны трупа, °С; t — ДНС, ч; k — коэффициент пропорциональности (постоянная охлаждения); Т_а — температура окружающей среды, °С.

Решением (1) является лежащее в основе известной уже более 150 лет простой экспоненциальной модели определения ДНС выражение [5, 6]:

 $T=T_a+(T_0-T_a) e^{-kt}$.

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

Цель исследования — разработка математической модели определения ДНС на основе закона охлаждения Ньютона–Рихмана в условиях изменяющейся температуры окружающей среды.

МАТЕРИАЛ И МЕТОДЫ ИССЛЕДОВАНИЯ

Методологический дизайн исследования представляет собой математическое моделирование охлаждения трупа в условиях изменяющейся температуры окружающей среды, осуществлённое на базе закона Ньютона–Рихмана и ориентированное на решение задачи термометрического определения ДНС. В качестве диаг­ностической зоны трупа использована подмышечная впадина, температура которой в момент смерти человека принята за 36,6°С. Вычислительные процедуры производили с использованием приложений Microsoft Excel пакета Office 2016 и Statistica (StatSoft) версии 7.0. Математические выражения выводили вручную, а также с использованием приложения Wolfram|Alpha. Код программы для ЭВМ составляли на языке программирования C# с использованием приложения Microsoft Visual Studio 2019.

РЕЗУЛЬТАТЫ

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

Тем не менее поставленная задача может быть решена в рамках обратного метода математического моделирования охлаждения трупа, в соответствии с которым за начальную температуру трупа принимается температура в диагностической зоне на момент повторной термометрии, а весь процесс охлаждения воспроизводится в обратном направлении (от температуры трупа в ходе его осмотра к температуре человека в момент его смерти). Для постоянной температуры внешней среды обратный метод моделирования на основе закона Ньютона–Рихмана (1) следует записать с противоположным знаком:

$\frac{dT}{dt}=k(T-T_a)$.                             (2)

В условиях изменяющейся температуры внешней среды показатель Т_а в выражении (2) нужно заменить соответствующей аппроксимирующей её функцией.

Обычно суточные колебания температуры воздуха в период до двух суток до осмотра трупа на месте его обнаружения удаётся аппроксимировать каким-либо представителем семейства полиномов 1–4‑й степеней:

$T_a=T_{a0}+{\beta }_{1}t+...+{\beta }_n t^n$,

где Т_а0 — температура окружающей среды на момент повторной термометрии трупа, °С; β₁, … β_n — коэффициенты полиномиальной аппроксимации; n — степень полинома; t — период времени до повторной термометрии трупа, ч.

Тогда модель обратного воспроизведения охлаждения трупа в соответствии с законом Ньютона–Рихмана описывается уравнением

$\frac{dT}{dt}=k(T-T_{a0}-{\beta }_{1}t-...-{\beta }_n t^n)$.       (3)

Общим решением (3) для полиномиальной аппроксимации 4‑й степени является функция

 $T=T_{a0}+{\beta }_1(t+\frac{1}{k})+{\beta }_2(t^2+\frac{2t}{k}+\frac{2}{k^2})+{\beta }_3(t^3+\frac{3t^2}{k}+\frac{6t}{k^2}+\frac{6}{k^3})+{\beta }_4(t^4+\frac{4t^3}{k}+\frac{12t^2}{k^2}+\frac{24t}{k^3}+\frac{24}{k^4})+(T_0+T_{a0}-\frac{{\beta }_1}{k}-\frac{2{\beta }_2}{k^2}-\frac{6{\beta }_3}{k^3}-\frac{24{\beta }_4}{k^4}) e^{kt}$ (4)

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

$T_1=T_{a2}+{\beta }_1(∆t+\frac{1}{k})+{\beta }_2(t^2+\frac{2∆t}{k}+\frac{2}{k^2})+{\beta }_3(∆t^3+\frac{3∆t^2}{k}+\frac{6∆t}{k^2}+\frac{6}{k^3})+{\beta }_4(∆t^4+\frac{4∆t^3}{k}+\frac{12∆t^2}{k^2}+\frac{24∆t}{k^3}+\frac{24}{k^4})+(T_2-T_{a2}-\frac{{\beta }_1}{k}-\frac{2{\beta }_2}{k^2}-\frac{6{\beta }_3}{k^3}-\frac{24{\beta }_4}{k^4}) e^{k∆t}$ (5)

где Т₁ и Т₂ — соответственно результаты первой и повторной термометрии трупа, °С; Та₂ — температура окружающей среды на момент повторной термометрии трупа, °С; Δt — временной интервал между первой и повторной термометриями трупа, ч.

Отсюда ДНС относительно момента первой термометрии трупа для указанной степени полиномиальной аппроксимации является решением неявной функции $36.6=T_{a1}+{\beta }_1(t+\frac{1}{k})+{\beta }_2(t^2+\frac{2t}{k}+\frac{2}{k^2})+{\beta }_3(t^3+\frac{3t^2}{k}+\frac{6t}{k^2}+\frac{6}{k^3})+{\beta }_4(t^4+\frac{4t^3}{k}+\frac{12t^2}{k^2}+\frac{24t}{k^3}+\frac{24}{k^4})+(T_1-T_{a1}-\frac{{\beta }_1}{k}-\frac{2{\beta }_2}{k^2}-\frac{6{\beta }_3}{k^3}-\frac{24{\beta }_4}{k^4}) e^{kt}$ (6)

Аналогичные решения для полиномиальных аппроксимаций меньших степеней получаются путём приравнивания отсутствующих регрессионных коэффициентов к нулю с последующим упрощением функций (4–6).

Пример 1

Осмотр трупа, обнаруженного на улице, начат в 02.00. В это время температура воздуха на уровне трупа составила 19,5°С. Температура в подмышечной впадине трупа на начало осмотра равнялась 32,7°С, а при повторном измерении через 1 ч после первого — 31,7°С. Температура воздуха на уровне трупа на момент его повторной термометрии составила 18,4°С. В предшествующие сутки в месте нахождения трупа имела место немонотонная динамика температуры воздуха, представленная в таблице.

Таблица. Реальная динамика температуры воздуха и её линейная и квадратичная полиномиальные аппроксимации 

Table. Real dynamics of air temperature and its linear and quadratic polynomial approximations

Несмотря на нелинейный характер температурного тренда в целом, динамика температуры воздуха на отрезке от 18.00 предыдущих суток до осмотра трупа могла быть хорошо аппроксимирована линейной функцией

$T_a=18.7615+0.8107t$,

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

$T=T_{a0}+{\beta }_1(t+\frac{1}{k})+(T_0-T_{a0}-\frac{\beta }{k}) e^{kt}$.

Коэффициент охлаждения k находим, решив уравнение

$T_1=T_{a2}+{\beta }_1(∆t+\frac{1}{k})+(T_2-T_{a2}-\frac{\beta }{k}) e^{k∆t}$,

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

$32,7=18,4+0,8107+\frac{0,8107}{k}+(31.7-18.4-\frac{0.8107}{k}) e^k$.

Отсюда k = 0,074661. ДНС определяем по формуле

$36.6=T_{a1}+{\beta }_1(t+\frac{1}{k})+(T_1-T_{a1}-\frac{{\beta }_1}{k}) e^{kt}$,

в соответствии с которой

$36.6=19.5+0.8107+\frac{0.8107}{0.074661}+(32.7-19.5-\frac{0.8107}{0.074661}) e^{0.074661t}$,

а ДНС равно 3,85 ч.

Динамику температуры в данном примере можно также с хорошей степенью приближения аппроксимировать квадратичным полиномом

$T_a=18.2583+1.3247t-0.0508t^2$,

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

В этом случае коэффициент охлаждения k найдём, решив неявную функцию

$T_1=T_{a2}+{\beta }_1(t+\frac{1}{k})+{\beta }_2(∆t^2+\frac{2∆t}{k}+\frac{2}{k^2})+(T_2-T_{a2}-\frac{{\beta }_1}{k}-\frac{2{\beta }_2}{k^2})e^{k∆t}$,

а ДНС — решив неявную функцию

$36.6=T_{a1}+{\beta }_1(t+\frac{1}{k})+{\beta }_2(t^2+\frac{2t}{k}+\frac{2}{k^2})+(T_1-T_{a1}-\frac{{\beta }_1}{k}-\frac{2{\beta }_2}{k^2})e^{kt}$.

После подстановки значений итерационным методом находим, что k = 0,07601, t = 4,02 ч.

Полученный результат близок таковому в модели с линейным температурным трендом окружающей среды. Подобное соответствие закономерно на приближенно линейном отрезке времени от осмотра трупа до 18.00 предыдущих суток.

Степень полинома, который следует использовать для аппроксимации динамики температуры окружающей среды перед обнаружением трупа, зависит от предполагаемой ДНС: чем больше ДНС, тем больше должна быть степень полиномиальной аппроксимации. В частности, полином 1‑й степени можно применять при ДНС не более 6–9 ч, когда температурный тренд окружающей среды имеет линейный характер. При увеличении ДНС до 20 ч температурный тренд окружающей среды, как правило, имеет форму унимодальной кривой, которую целесообразно аппроксимировать кубическим полиномом. При ДНС от 20 до 30 ч динамика температуры окружающей среды принимает S‑образную форму с наличием одного локального максимума и одного локального минимума, аппроксимация которой требует уже использования кубического полинома. Наконец, при ДНС более 30 ч температурный тренд окружающей среды приобретает вид двугорбой кривой, адекватная аппроксимация которой становится возможной лишь при использовании полинома 4‑й степени. При этом хорошие аппроксимации можно получить путём оценивания анализируемых температурных кривых методом наименьших квадратов. Принципы обратного математического моделирования охлаждения трупа для полиномиальных аппроксимаций 3‑й и 4‑й степени аналогичны таковым для линейной и квадратной полиномиальных аппроксимаций и отличаются от последних лишь повышением сложности выполняемых математических операций.

Пример 2

На начало осмотра трупа, обнаруженного на улице, температура воздуха на уровне трупа составила 16,4°С. Температура в подмышечной впадине трупа на начало осмотра равнялась 20,83°С, а при повторном измерении через 1 ч после первого — 20,57°С. Температура воздуха на уровне трупа на момент его повторной термометрии составила 16,0°С. В предшествующие двое суток в месте нахождения трупа имела место сложная динамика температуры воздуха с наличием двух локальных максимумов (рис. 1).

Рис. 1. Модельные кривые охлаждения трупа и температуры окружающей среды, рассчитанные на основе полиномиальной аппроксимации 4-й степени.

Динамика температуры воздуха в предшествующие осмотру трупа двое суток могла быть адекватно аппроксимирована полиномом 4‑й степени:

$T_a=14.5732+2.6364t-0.31266t^2+0.01179t^3-0.00014t^4$.

Коэффициент охлаждения k находим, решив неявную функцию

$16+\frac{2.6364}{k}-\frac{2*0.31266}{k^2}+\frac{6*0.01179}{k^3}-\frac{24*0.00014}{k^4}+e^k(20.57-16-\frac{2.6364}{k}+\frac{2*0.31266}{k^2}-\frac{6*0.01179}{k^3}+\frac{24*0.00014}{k^4})=20.83$

Итерационным методом находим, что k = 0,074371. Для определения ДНС необходимо найти переменную t в уравнении

$36.6=16.4+2.6364(t+\frac{1}{0.074371})-0.31266(t^2+\frac{2t}{0.074371}+\frac{2}{0.{074371}^2})+0.01179(t^3+\frac{3t^2}{0.074371}+\frac{6t}{0.{074371}^2}+\frac{6}{0.{074371}^2})-$

$-0.00014(t^4+\frac{4t^3}{0.074371}+\frac{12t^2}{0.{074371}^2}+\frac{24t}{0.{074371}^3}+\frac{24}{0.{074371}^4})+e^{0.074371t}(20.83-16.4-\frac{2.6364}{0.74371}+\frac{2*0.31266}{0.{074371}^2}-\frac{6*0.01179}{0.{074371}^3}+\frac{24*0.00014}{0.{074371}^4})$

решением которого является значение ДНС, равное 42,1 ч.

При расчётах на основе зарегистрированного в момент первой термометрии трупа значения температуры воздуха без учёта её предшествующих изменений (т. е. принимая температуру среды при первой термометрии трупа за константу) получаем ДНС равную 25,1 ч. Столь значительная ошибка объясняется тем, что за 10 ч до осмотра трупа температура окружающей среды сначала сравнялась, а затем превысила температуру в диагностической зоне трупа (см. рис. 1). Из-за этого температура трупа стала расти, а затем сразу после нового падения температуры окружающей среды вновь начала уменьшаться, после чего подверглась динамической термометрии в ходе осмотра трупа.

Приведённый пример показывает, что метод обратного математического моделирования охлаждения трупа учитывает практически любые изменения температуры окружающей среды, в том числе и такие, когда температура воздуха превышает температуру трупа. Однако корректное применение закона Ньютона–Рихмана при определении ДНС требует правильного выбора диагностической зоны в трупе, т. к. данный закон не предназначен для описания процесса охлаждения глубоких тканей трупа. Следовательно, изложенные методы математического моделирования не подходят к таким видам посмертной термометрии, как ректальная, печёночная и краниоэнцефальная. Зато перспективными в этом отношении являются тимпаническая и внутриглазная разновидность термометрии [7, 8]. Указанные диагностические зоны по своим теплообменным свойствам относятся к поверхностным тканям трупа, но в отличие от кожи в силу своего анатомического расположения в гораздо меньшей степени подвержены влиянию различных случайных факторов, включая инсоляцию и иные погодные условия, особенности ложа трупа, наличие и характер одежды. Для начала практического использования указанных диагностических зон в рамках предложенной математической модели достаточно стандартизации измерительных процедур и определения средних температур указанных зон на момент смерти человека.

Детерминистский характер простой экспоненциальной модели охлаждения трупа позволяет рассматривать определение ДНС как косвенное измерение, результат которого связан с результатами прямых измерений температуры трупа и окружающей среды и интервала времени между первой и повторной термометрией дифференцируемой нелинейной зависимостью. Указанное свойство позволяет вычислять предельные абсолютные погрешности определения ДНС непосредственно, без неизбежной для остальных моделей проверки на трупном материале. Математическая модель оценки погрешностей определения ДНС на основе закона охлаждения Ньютона–Рихмана и примеры её реализации будут охарактеризованы нами в отдельном сообщении.

Для избавления пользователя от выполнения трудоёмких математических операций изложенный алгоритм определения ДНС и её погрешностей был реализован в формате программы для ЭВМ Warm Bodies NRN (свидетельство о государственной регистрации № 2021611972). Данная программа предназначена для определения ДНС на основе закона охлаждения Ньютона–Рихмана как в условиях постоянной, так и изменяющейся температуры внешней среды. Диагностической точкой в рамках данного физического закона могут служить только поверхностные ткани трупа. Программа поддерживает любые полиномиальные аппроксимации температурного тренда внешней среды 1–4‑й степени. При постоянной температуре окружающей среды приложение находит ДНС аналитически по формуле H. Rainy [6], являющейся следствием уравнения (1). В условиях изменяющейся температуры внешней среды программа определяет ДНС, решая неявные функции (5) и (6), предварительно сгруппировав все слагаемые и множители в правых частях уравнений и приравняв их к нулю. Поиск корней неявных функций приложение реализует методом касательных Ньютона, который имеет высокую скорость сходимости, и с его помощью приближенное значение корня с точностью до 6 знаков находится всего за 3–15 итераций. Затем программа рассчитывает предельные абсолютные погрешности определения ДНС с учётом инструментальных погрешностей средств измерения температуры и времени на основе математической модели косвенного измерения. Стартовое окно программы приведено на рис. 2.

Рис. 2. Стартовое окно приложения Warm Bodies NRN с нередактируемым выпадающим списком вариантов наличия или отсутствия изменений температуры внешней среды за предшествующие двое суток перед обнаружением трупа и текстовыми полями для ввода данных динамической термометрии трупа и внешней среды, а также инструментальных погрешностей измерительных средств.

ЗАКЛЮЧЕНИЕ

Разработан метод математического моделирования охлаждения трупа, позволяющий объективно определять ДНС в условиях изменяющейся температуры внешней среды. Метод основан на законе охлаждения Ньютона–Рихмана и реализуется путём двукратной термометрии поверхностных тканей трупа и внешней среды. Разработанная математическая модель определения ДНС и её предельных абсолютных погрешностей реализована в формате прикладной программы Warm Bodies NRN.

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

About the authors

German V. Nedugov

Author for correspondence.
Email: nedugovh@mail.ru
ORCID iD: 0000-0002-7380-3766
SPIN-code: 3828-8091
Scopus Author ID: 25947646500

Cand. Sci. (Med.), Assistant Professor

References

Supplementary files

Supplementary Files

Action

1. JATS XML

Download

2. Fig. 1. Model curves of cadaver cooling and ambient temperature calculated on the basis of a polynomial approximation of the 4th degree.

Download (756KB)

Indexing metadata

3. Fig. 2. Start window of the Warm Bodies NRN application with a noneditable drop-down list of options for the presence or absence of changes in the ambient temperature for the previous two days before the discovery of the corpse and text fields for entering data on the dynamic thermometry of the corpse and external environment and errors of measuring instruments.

Download (651KB)

Indexing metadata