Я защищался по акустотермографии, методу, который позволяет услышать (!) температуру. Расскажу, что это. И спрошу помощи у понимающих в радиоэлектронике.
![Рис. 1 Поработали на тренажере – бицепс нагрелся и зашумел Рис. 1 Поработали на тренажере – бицепс нагрелся и зашумел](https://habrastorage.org/getpro/habr/upload_files/ea6/242/f44/ea6242f44c7b4c735e55e28037e73349.jpg)
Температура – это мера средней кинетический энергии движения атомов и молекул. Электронщики знают о Джонсоновском тепловом шуме резистора, так вот в воде распространяется этой же природы механический тепловой шум в ультразвуковом диапазоне. Энергии молекул имеют некоторое статистическое распределение. И случаются флуктуации, при которых локальное повышение плотности, т.е. возмущение, приводит к возбуждению механической волны. Эта механическая волна – ультразвук. Таким образом вода и мягкие ткани человека наполнены идущим во всех направлениях ультразвуком, несущим информацию о температурном распределении.
Спектр теплового излучения соответствует формуле Рэлея–Джинса – спектральная плотность излучения с ростом частоты растет, т.к. в любой области пространства с ростом частоты растет количество мод стоячих колебаний (степеней свободы), а в соответствии с законом равнораспределения энергии на каждую такую моду должна приходиться равная энергия.
![Рис. 2 В оптике и радиоэлектронике зависимость спектральной плотности от частоты растет, а затем падает.. В акустике только растет. Рис. 2 В оптике и радиоэлектронике зависимость спектральной плотности от частоты растет, а затем падает.. В акустике только растет.](https://habrastorage.org/getpro/habr/upload_files/bed/184/cc7/bed184cc7e6932589ab2651d45268978.png)
В электромагнитной теории этот рост приводит к парадоксу «ультрафиолетовой катастрофы», который разрешается квантованием энергии волны и переходом зависимости Релея-Джинса в зависимость Планка. В акустике такого перехода не происходит, т.к. на достаточно больших частотах длина волны оказывается меньше атомарных расстояний и достаточно высокочастотных волн просто нет.
Рост интенсивности теплового шума с частотой мешает акустикам, слушающим что-либо еще кроме теплового шума, т.к. тепловой шум заглушает остальные сигналы. Нам же, термоакустикам, помогает. Однако и нам не интересно слушать слишком высокие частоты, т.к. с ростом частоты растет и коэффициент поглощения. Наиболее высокочастотные тепловые шумы – это и есть тепловое движение атомов и молекул, которое несет информацию только о состоянии ближнего слоя. Низкие же частоты проходят с незначительным затуханием большие расстояния, перенося информацию о тепловом распределении в среде на большой глубине. Поэтому мы обычно слушаем на частоте около 2 МГц.
Как акустический шум слушать? – пьезопластинкой с напылёнными электродами. Обычная пьезопластинка, которая при подаче переменного напряжения излучает ультразвук, будет микрофоном, если напряжение с нее снимать.
![Рис. 3 Обычный пьезоизлучатель, он же приемник УЗ-волны, на частоте равной частоте механического резонанса керамической пластинки. Рис. 3 Обычный пьезоизлучатель, он же приемник УЗ-волны, на частоте равной частоте механического резонанса керамической пластинки.](https://habrastorage.org/getpro/habr/upload_files/bc1/b38/d9d/bc1b38d9de21fe3495c4feb5a65c9af9.jpg)
Микрофон получается нелинейным, снимает только частоты, с которыми пластинка входит в резонанс. Спектр теплового шума колоколообразный вокруг резонансной частоты пластинки. Это не значит, что спектр теплового шума на самом деле колоколообразный, это значит, что мы умеем только так слушать.
Спектр теплового излучения сплошной или шумоподобный. И спектр непериодического сигнала тоже шумоподобный. Это значит, что воздействие шумоподобной волны аналогично непериодическому сигналу, как если бы время от времени по пластинке ударяли. И под такими ударами пластинка колеблется с частотой собственных свободных колебаний.
С одной стороны пластинки должен быть воздух, с другой вода, иначе энергия волны как придет на пластинку, так и дальше в воду уйдет. А нужно, чтобы из спектра приходящей волны вырезались частоты близкие к резонансной, энергия этих волн пластинку раскачивала, сигнал был по возможности сильным и можно было измерять его среднеквадратичную амплитуду, пропорциональную температуре излучающей среды.
Но измерять этот синусоидальный сигнал с нановольтовой амплитудой непросто. Наши коллеги из ИПФ РАН подобные датчики делают. Я пытался датчик повторить под свои нужды, но не получается. Как измерять амплитуду синусоидального сигнала размером в нановольты? Может быть АМ-радиомикросхемой на частоте 2МГц? Или усилителем подобным вот такому? Есть опыт и умение делать УЗ-датчики, программировать микроконтроллеры, разводить платы, использовать АЦП, обрабатывать данные. Но нет понимания, как и какие использовать чувствительные компактные усилители. Если кто понимает, как это делается, пожалуйста, помогите разъяснением или сотрудничеством (e-mail: iva2000@gmail.com).
Если волна пришла на датчик под прямым углом – она воздействует на всю пластину одновременно, раскачивая ее. Волна, падающая на пластину под косым углом, в некоторых областях пластинки создает положительное давление, в других в тот же момент отрицательное, и суммарное воздействие оказывает слабее. Из-за этого эффекта пластинка слушает только волны, приходящие прямо на нее и почти не слышит волны, приходящие под углом. И с ростом частоты диаграмма чувствительности сжимается. То есть пластинка слушает в узком направлении, имеет узкую диаграмму направленности, работает как направленный микрофон. И это можно использовать!
![Рис. 4 Распределение чувствительности УЗ-приемника от направления – теория и эксперимент. Рис. 4 Распределение чувствительности УЗ-приемника от направления – теория и эксперимент.](https://habrastorage.org/getpro/habr/upload_files/f29/31d/30f/f2931d30f7d36a28a401caf0427620d9.png)
Я смоделировал самодельную нейросетку, угадывающую положение источника волны. Комплект из десяти датчиков дает низкое разрешение по глубине, но на порядок лучшее по линии, вдоль которой стоят датчики.
![Рис. 4 Восстановление положения и размеров источника сигнала по данным десяти датчиков (компьютерный эксперимент) Рис. 4 Восстановление положения и размеров источника сигнала по данным десяти датчиков (компьютерный эксперимент)](https://habrastorage.org/getpro/habr/upload_files/022/884/e65/022884e6535438d76949a8f179fcbc7f.png)
Затем поставили в аквариум пачку датчиков, двигали перед ними вытянутым горячим объектом и «слушали» как он двигается. Результаты натурного эксперимент впечатляют меньше компьютерного моделирования – нужно использовать побольше датчиков поменьше размером, а больше нет.
![Рис. 7 Экспериментальное восстановление образа, перемещающегося вытянутого нагретого тела (натурный эксперимент) Рис. 7 Экспериментальное восстановление образа, перемещающегося вытянутого нагретого тела (натурный эксперимент)](https://habrastorage.org/getpro/habr/upload_files/9d3/5b6/010/9d35b6010033b6dd9ca289195a6869a9.png)
А вот как звучат тающий лед, и пузырьки в газировке. Также увлекательно как слушать пение китов в океане, только в микромасштабе.
![Рис. 5 Нетепловые акустические сигналы в воде
Рис. 5 Нетепловые акустические сигналы в воде](https://habrastorage.org/getpro/habr/upload_files/33f/214/164/33f2141644310b0c8e1b845fb33a12fb.png)
Моделируя операции с термообляцией использовали говяжью печень с погруженным в нее жалом паяльника. Измеряли распределение температур термодатчиками, и пытались его же восстановить акустотермометром. Включали и выключали паяльник и смотрели как сигналы растут и снижаются, и видели, что чем дальше «луч» от паяльника, тем меньше амплитуда изменений. Слушали печень и с одного направления, и с разных.
![Рис. 6 Паяльник в говяжьей печени как модельный источник тепла Рис. 6 Паяльник в говяжьей печени как модельный источник тепла](https://habrastorage.org/getpro/habr/upload_files/6ba/829/124/6ba82912493fe4c63c0d379679f7e739.png)
А еще можно использовать корреляционный прием! Строить корреляционную функцию сигналов с двух пересекающихся датчиков, сдвигая сигналы по времени. И получать сигнал отличный от нуля только из-за тех «событий», волна от которых приходит на оба датчика с соответствующей задержкой. Остальные события корреляционная функция обнулит. И получится что-то вроде фазированной антенной решетки, которая умеет менять направления приема, не поворачиваясь физически.
Расчеты и эксперименты показывают, что за 10…50 с измеряется температура в объеме 1 см3 с точностью 0,5…1 K на глубине до 5…10 см. Акустической термографией можно контролировать глубинную температуру тканей во время медицинских процедур, и диагностировать заболевания, сопровождающиеся глубинным нагревом.
Недавно исследовали грудную клетку 59-летнего человека, перенесшего COVID-19 с подтвержденным на томографе двусторонним поражением легких. Во время болезни работать не было возможности, обследовали сразу после окончания карантина. Интегральная температура тканей за 10 дней выросла с 32,2 ± 0,07 до 33,0 ± 0,03°С. Выход температур на плато может свидетельствовать о завершении реабилитации.
![Рис. 7 Изменения акустояркостной температуры торса (1), со стороны груди (2) и со спины (3) во время реабилитации после COVID-19. Рис. 7 Изменения акустояркостной температуры торса (1), со стороны груди (2) и со спины (3) во время реабилитации после COVID-19.](https://habrastorage.org/getpro/habr/upload_files/677/9f6/0d9/6779f60d99c256652f97d5dee3852637.png)
Коэффициент корреляции между средней акустояркостной температурой и температурой комнаты равен –0,14, т. е. корреляция отсутствует. Контрольная инфракрасная термометрия изменений температуры кожных покровов не показала. Идеальное исследование должно включать МРТ как контрольный метод, но МРТ – дорого. А от предложения для контроля вводить в пациента толстыми полыми иглами термопары участники решительно отказались.
И тем не менее, качественный результат есть. Легочные патологи ухудшают кровоснабжение легких и бронхов, что может снизить температуру тканей. И термометрируя легкие после выписки можно наблюдать как восстанавливается их кровенаполнение.
Дальнейшая работа, очевидно, связана с увеличением объема выборки.
Работа поддержана РФФИ: грант 20–02–00759 «Восстановление с помощью многоканального акустотермометра пространственного распределения меняющейся во времени глубинной температуры тела человека при локальной гипертермии».
Комментарии (5)
abbyyit
14.01.2022 12:18+1измерять амплитуду синусоидального сигнала размером в нановольты? Может быть АМ-радиомикросхемой на частоте 2МГц?
цифровые и аналоговые синхронные детекторы отлично выполняют измерение нановольтовых и фемтоамперных сигналов
синхронное детектирование АМ сигналов привлекает сейчас внимание как специалистов, так и радиолюбителей(
iva2000 Автор
14.01.2022 13:57Очень интересно, внимательно читаю ваши ссылки!
abbyyit
14.01.2022 14:16не факт, что получится аж так оптимизировать)) понадобится приличная фантазия
даже профильщики к акустотермографу приделывали нейросети
Динамическое восстановление глубинной температуры методом акустотермографии с помощью нейронных сетей (с. 768-772 (2013)
может проще будет запилить ацп+датчик температуры
или юзать модуль температуры запьястья
ru-bezh.ru/blog/36861-tochnost-distanczionnogo-opredeleniya-temperaturyi-tela-chelovek
Radisto
Что-то вроде акустико-эмиссионного контроля, но для человека?
iva2000 Автор
Этимологически это действительно акустико-эмиссионный контроль. Прямо слово в слово.
Но акустико-эмиссионный контроль проводят под калиброванной нагрузкой, слушая как система из-за нагрузки трещит. А при акустотермии вообще никак не воздействуют на систему, только слушают как она шумит сама.