Python-инструмент для оценки нейтронной дозы вне поля при протонной терапии

Протонная терапия считается одним из самых точных методов лечения рака, однако у нее есть скрытый недостаток: при облучении опухоли протоны взаимодействуют с веществом, порождая вторичные нейтроны. Эти нейтроны рассеиваются по процедурной комнате и создают дополнительную дозу облучения для здоровых тканей, что может увеличивать риск развития вторичных злокачественных новообразований. До недавнего времени точный расчет такой «внепольной» нейтронной дозы был сложной задачей, требующей длительного компьютерного моделирования. Теперь группа исследователей из Медицинского центра Университета Радбауд в Нидерландах создала простой в использовании Python-калькулятор, который решает эту проблему, позволяя быстро оценивать нейтронную дозу в любой точке процедурной комнаты на основе данных плана лечения.

Новый инструмент для оценки нейтронной дозы при протонной терапии

Разработанный инструмент, описанный в статье, опубликованной в журнале Physics in Medicine & Biology, представляет собой скрипт на Python, который на основе данных из плана лечения в формате DICOM RT рассчитывает эквивалентную дозу нейтронов. Калькулятор использует аналитические модели, полученные из результатов моделирования методом Монте-Карло, и требует всего несколько секунд для выполнения расчета. Это делает его пригодным для клинического использования — например, для быстрой оценки риска для конкретного пациента или для оптимизации расположения оборудования в комнате. Авторы протестировали инструмент на реальных клинических случаях и показали, что его точность сопоставима с результатами полномасштабного моделирования Монте-Карло: погрешность составляет менее 20% для большинства точек, что приемлемо для дозиметрических оценок. Инструмент доступен для скачивания на GitHub под лицензией с открытым исходным кодом, что позволяет любому медицинскому физику адаптировать его под свои нужды.

Какую проблему решает Python-калькулятор нейтронной дозы?

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

Как работает Python-калькулятор?

Инструмент принимает на вход файл плана лечения в формате DICOM RT, который содержит информацию о геометрии пучка, его энергии, положении пациента и параметрах облучения. На основе этих данных калькулятор вычисляет распределение нейтронного поля в комнате, используя предварительно рассчитанные коэффициенты, полученные из моделирования Монте-Карло для типичных конфигураций. Пользователь может задать координаты интересующей точки, например место расположения гонад пациента или хрусталика глаза, и получить оценку эквивалентной дозы нейтронов. Результат выводится в виде числа и графического отображения на схеме комнаты. Для проверки точности авторы сравнили результаты калькулятора с независимыми измерениями на реальном аппарате — расхождение не превысило 15% для большинства точек.

Каковы технические детали создания инструмента?

Исследователи провели обширное моделирование Монте-Карло с помощью кода Geant4 для типичного аппарата протонной терапии с пассивным рассеянием. Они варьировали энергию протонов от 70 до 230 МэВ, размер поля, глубину мишени и другие параметры. На основе полученных данных были построены аналитические функции, аппроксимирующие зависимость нейтронной дозы от расстояния до изоцентра, угла и энергии пучка. Эти функции и легли в основу Python-калькулятора. Такой подход позволил сократить время расчета с часов до секунд, сохранив при этом приемлемую точность.

Кому пригодится новый инструмент и как его использовать?

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

Перспективы развития и внедрения

Авторы планируют расширить возможности калькулятора для других типов ускорителей, в частности с активным сканированием пучка, и включить расчет нейтронной дозы для персонала. Также ведется работа по интеграции инструмента в коммерческие системы планирования лечения. Ожидается, что в ближайшие год-два подобные калькуляторы станут стандартной частью рабочего процесса протонной терапии. Открытый исходный код способствует широкому внедрению и дальнейшему совершенствованию метода.

Итог

Разработанный Python-инструмент представляет собой важный шаг к повышению безопасности протонной терапии. Он позволяет быстро и точно оценивать внепольную нейтронную дозу, что критически важно для минимизации риска вторичных раков. Следить за этой темой стоит всем, кто интересуется современными технологиями лучевой терапии и их влиянием на здоровье пациентов.