Введение 2
Постановка задачи 4
1 Система 5
1.1 Существующие решения 5
1.2 Требования к системе 6
1.2.1 Требования к характеристикам компьютера, обрабатывающего акустические данные 6
1.2.2 Требования к входным данным 6
1.2.3 Требования к процессу обработки 8
2 Прототип системы для оценивания характеристик 10
2.1 Компонент «Моделирование волны» 11
2.1.1 Обзор: способы моделирования 11
2.1.2 Детали реализации 13
2.1.3 Векторизация 14
2.2 Компонент «Определение характеристик» 15
2.2.1 Оптимизационная задача 15
2.2.2 Детали реализации 16
2.3 Пользовательский интерфейс 19
3 Экспериментальная апробация 21
3.1 Эксперимент с данными с отсутствием
объектов в водной среде 22
3.2 Эксперимент с данными по объекту в водной среде . 22
Заключение 25
Список литературы 26
Решение обратных задач по акустическим данным — один из известных методов медицинского обследования, сейсморазведки, дефектоскопии.
Ультразвуковое исследование (УЗИ) внутренних органов является самым известным применением такого рода задач. Если рассмотреть самый обыкновенный УЗИ-аппарат, то его датчик, излучая ультразвук, принимает отраженный ультразвук, обрабатывает его, зная как и в каких средах он распространяется, и выдает изображение. Этот метод намного быстрее других, на приеме у врача результаты УЗИ можно получить в тот же день.
Для сейсморазведки этот метод является важным, так как позволяет существенно сократить затраты, например, для бурения скважин, где на счету стоит очень много: разведка местности, оценка месторождений, и непосредственно бурение. Если хотя бы на одном этапе уже на месте обнаружится ошибка, будет потеряно очень много денег, не говоря уже об угрозах для жизни, которые могут возникнуть во время бурения. В процессе сейсморазведки применяются разные способы: с применением взрывчатых и невзрывчатых веществ. В независимости от того, что именно используется, применятся это для того, чтобы в толще породы образовались (возбудились) упругие волны. Эти самые волны фиксируются раставлен- ными в нужных местах датчиками, таким образом получаются сейсмограммы. В дальнейшем остается только обработать эти сейсмограммы, чтобы найти полезные ископаемые, например, газ, нефть, металлы, минералы и т.д.
Также этот метод очень важен в дефектоскопии. Возможность узнать заранее о том, что те или иные деталь, резервуар, балон и т.д. неисправны или содержат внутри дефекты и микротрещены, может существенно помочь в своевременной замене этих изделий, а также обеспечит безопасность для производства и персонала
При обработке акустических данных требуется найти звуковые характеристики объектов в исследуемой области, например внутри детали для дефектоскопии или органа для УЗИ. Эти характеристики дожны бать такими, чтобы расчётные данные как можно лучше приближали экспериментальные. Значения таких характеристик представитляется в виде распределения скорости звука в среде.
Поиск распределения включает в себя такие этапы: с текущим полем скорости звука моделируется распространение звуковых волн на протяжении определенного количества времени; во время моделирования волн с записывающих датчиков на протяжении этого же времени снимаются диаграммы, представляющие собой значения амплитуд звукового сигнала; после моделирования и снятия диаграмм путем варьирования характеристик, ставится задача уменьшить разницу между извлеченными диаграммами и экспериментальными, представленная минимизацией функционала ошибки между всеми амплитудами записанных диаграмм. Описанный функционал оптимизируется с помощью метода наименьших квадратов, часто применяемый в задачах статистической регрессии и методах оптимизации. В качестве основного метода решения задачи о наименьших квадратов выбран метод Левенберга - Марквардта...
Получены следующие результаты:
• сформулированы требования к системе;
• разработана архитектура системы;
• разработан прототип системы с применением
- векторизации через библиотеку PyTorch,
- CUDA-ядер Nvidia,
- SPSA;
• прототип применен к акустическим данным с отсутствием объектов в водной среде и одному объекту в водной среде.
Система применена в международном проекте университета с Ху-ачжунским университетом науки и технологии города Ухань.
[1] Галямина В. С. Вычисление времени прихода («Time Of Flight») акустических сигналов. - Конференция СПИСОК2019, 2019.
[2] Гонта К. А. Технология восстановления особых областей на основе данных акустической томографии. Санкт-Петербург, 2019.
[3] Леонова А. В. Система для расчета скоростей звука в особых областях по данным УЗИ - томографии. Санкт-Петербург, 2020.
[4] Сенин И.И. Рандомизированный алгоритм при обработке данных ультразвуковых исследований. - Стохастическая оптимизация в информатике, 2016.
[5] Турсунова М.Б. Рандомизированный алгоритм при обработке данных ультразвуковых исследований. - Конференция СПИСОК2019, 2019.
[6] Bohme C., Holmberg A., Nilsson Lind M. Numerical Analysis of the Two Dimensional Wave Equation: Using Weighted Finite Differences for Homogeneous and Hetrogeneous Media. - 2020.
[7] Duric N. et al. Detection of breast cancer with ultrasound tomography: First results with the Computed Ultrasound Risk Evaluation (CURE) prototype //Medical physics. - 2007. - Т. 34. - №. 2. - С. 773-785.
[8] Erofeeva V. et al. Detection of specific areas and densities for ultrasound tomography //Cybernetics and Physics. - 2019. - Т. 8. - №. 3. - С. 121-127.
[9] Granichin O., Amelina N. Simultaneous perturbation stochastic approximation for tracking under unknown but bounded disturbances //IEEE Transactions on Automatic Control. - 2014. - Т. 60. - №. 6. - С. 1653-1658.
[10] Hutchison D. W. On an efficient distribution of perturbations for simulation optimization using simultaneous perturbation stochastic approximation //Proceedings of IASTED International Conference. - 2002. - С. 4-6.
[11] Jovanovic I. Inverse problems in acoustic tomography: theory and applications : дис. - Verlag nicht ermittelbar, 2008.
[12] Li C. et al. In vivo breast sound-speed imaging with ultrasound tomography //Ultrasound in medicine biology. - 2009. - Т. 35. - №. 10. - С. 1615-1628.
[13] Matthews T. P. et al. Regularized dual averaging image reconstruction for full-wave ultrasound computed tomography //IEEE transactions on ultrasonics, ferroelectrics, and frequency control. - 2017. - Т. 64. - №.
5. - С. 811-825.
[14] Peng J., Liu J., Wei H. A compressed matrix sequence method for solving normal equations of bundle adjustment //Machine Vision and Applications. - 2021. - Т. 32. - №. 4. - С. 1-10.
[15] QT Ultrasound. - https://www.qtimaging.com/publications/ (дата обращения: 1 мая 2023)...29