뇌종양 광음향 영상 시뮬레이션 모델

초록

본 논문은 MRI 기반 디지털 뇌 팬텀을 이용해 광음향 영상 재구성을 전산적으로 평가하는 방법을 제시한다. 유한요소법으로 광전달을 모델링하고, 펄스 광음향 원리를 통해 초기 압력장을 계산한다. 선형화된 1차 연속 방정식과 k‑space 알고리즘을 사용해 탄성 두개골을 포함한 전파를 시뮬레이션하며, 무작위 매질 접근법으로 산란, 파워법칙으로 흡수를 구현한다. 최종적으로 k‑space 시간역전법으로 압력 신호를 재구성해 3차원 뇌종양 광음향 단층 영상을 얻는다. 모델의 수렴성과 실용성을 검증하였다.

상세 분석

이 연구는 광음향 영상(PAI)의 뇌 적용 가능성을 검증하기 위해 전 단계 물리 현상을 통합적으로 모델링한 점이 가장 큰 강점이다. 첫 번째로, MR 영상을 기반으로 한 이미지 세분화 과정을 통해 실제 해부학적 구조를 반영한 디지털 팬텀을 구축하였다. 이는 조직별 광학적 흡수계수와 산란계수를 정확히 할당할 수 있게 하여, 광전달 시뮬레이션의 현실성을 크게 높인다. 두 번째로, 유한요소법(FEM)을 이용한 광전달 모델은 복잡한 뇌 조직의 이방성 및 비균질성을 고려할 수 있다. 특히, 혈액 내 헤모글로빈 흡수 스펙트럼을 반영함으로써 종양 혈관 신생과 같은 병리학적 특징을 대비할 수 있다.

압력 생성 단계에서는 펄스 레이저에 의해 순간적인 열 팽창이 일어나 초기 압력 파동을 만든다는 전통적인 광음향 원리를 적용했으며, 이때 열-기계 변환 효율을 조직별 비열용량과 열전도도와 연계해 정량화하였다. 이어서, 선형화된 연속 방정식과 운동 방정식을 1차 형태로 결합한 연속 파동 방정식을 사용해 압력 전파를 기술하였다. 여기서 핵심은 k‑space 방법을 적용한 점이다. k‑space 알고리즘은 고주파 성분의 정확한 전파와 수치적 안정성을 동시에 제공하며, 특히 탄성 매질(두개골)에서의 전파를 다루기 위해 기존의 액체 전용 k‑space를 확장한 것이 눈에 띈다.

두개골은 고강성 탄성체이며, 초음파가 통과할 때 산란과 흡수가 크게 발생한다. 이를 모델링하기 위해 무작위 매질(random medium) 접근법을 도입해 공간적 변동성을 가진 탄성 계수를 부여하고, 파워법칙(power law) 형태의 흡수 모델을 적용하였다. 파워법칙은 조직별 주파수 의존성 흡수를 현실적으로 재현하며, 특히 고주파 손실이 심한 두개골을 정확히 묘사한다.

마지막으로, 전방 시뮬레이션으로 얻어진 압력 데이터를 k‑space 시간역전(time reversal) 기법에 투입해 재구성하였다. 시간역전은 경계 조건과 매질 이질성을 자연스럽게 반영하면서 역방향 파동을 전파시켜 원래의 초기 압력 분포를 복원한다. 3차원 전역 재구성을 수행함으로써 종양의 위치, 형태, 그리고 혈관 분포까지 시각화할 수 있었다.

전체 모델은 각 단계별 수치적 수렴성을 검증했으며, 파라미터 스윕을 통해 격자 해상도, 시간 스텝, 그리고 무작위 매질의 통계적 특성이 최종 영상 품질에 미치는 영향을 정량화하였다. 결과적으로, 제안된 통합 시뮬레이션 프레임워크는 실제 실험적 PAI 시스템 설계와 최적화, 그리고 임상 적용 전 가상 검증 단계에 매우 유용한 도구가 될 것으로 기대된다.