단일면 접근으로 구현하는 가상 음향 및 복합 초점 함수

초록

본 논문은 전면 전체가 아닌 한쪽 면만 접근 가능한 복합 매질에서, 시간역전 대신 복합 초점 함수를 이용해 내부에 가상 음원·수신기를 생성하는 방법을 제시한다. Marchenko 방법을 기반으로 한 단일면 초점 설계와 대칭화 과정을 통해 실제 물리적 소스·레시버 없이도 매질 내부의 복잡한 전파·다중산란 경로를 재현한다. 초음파 물리 모델과 실제 지진 반사 데이터를 이용한 실험 결과가 이를 입증한다.

상세 분석

이 연구는 전통적인 시간역전(Time‑Reversal, TR) 기법이 전체 경계면에서 측정·재발사된 파동을 필요로 하는 한계를 극복하고자 한다. TR은 손실이 없는 매질에서 파동 방정식이 시간에 대해 대칭이라는 사실을 이용해, 내부에 존재하는 점음원의 응답을 경계에서 기록한 뒤 이를 시간역전시켜 다시 방출함으로써 가상 음원을 재현한다. 그러나 매질이 비균질이고 임피던스 대비가 크게 변하는 경우, 특히 접근 가능한 경계가 한쪽 면에 국한될 때는 TR이 정확히 수렴하지 않는다. 이때 발생하는 ‘ghost focus’와 비등방성 방사 패턴은 실제 응용에서 큰 장애가 된다.

저자들은 이러한 문제를 해결하기 위해 ‘단일면 초점 함수(single‑sided focusing function)’ F(x,s,t)를 도입한다. F는 경계면 S₀(접근 가능한 면)에서만 정의되며, Marchenko 방정식의 다차원 확장을 이용해 반사 데이터 G(x₀,x,t)로부터 직접 추정한다. 핵심 아이디어는 (1) 직접 도착 파동을 배경 속도 모델을 이용해 초기 초점 함수의 근사값으로 사용하고, (2) 식 (2)와 (3)에서 제시된 대칭화 연산을 통해 시간에 대해 반대칭인 인공 아티팩트를 소거한다는 것이다. 대칭화된 형태는 G(r,s,t)+G(r,s,−t)=Symmetrize∫_{S₀}G(r,x,t)*F(x,s,t)dx와 같이, 전체 경계가 아니라 접근 가능한 면만을 적분함에도 불구하고 전체 TR과 동일한 좌변을 얻는다.

수학적으로는 Huygens 원리를 재해석하여, 경계면의 ‘2차 소스’가 실제 파동이 아니라 기록된 반사 응답과 초점 함수의 컨볼루션으로 구성된다는 점을 강조한다. 또한, 반사 데이터가 스칼라 파동(음향)뿐 아니라 손실이 없는 전자기·탄성 파동에도 적용 가능함을 이론적으로 제시한다. 손실이 존재하는 경우에는 사전 손실 보정(pre‑compensation) 절차를 통해 Marchenko 알고리즘을 그대로 사용할 수 있다.

실험 부분에서는 두 가지 사례를 제시한다. 첫 번째는 실험실에서 제작한 3차원 실리콘 겔·밀랍 복합 매질이다. 물탱크 위에 초음파 트랜스듀서를 배치하고, 대각선 경로를 따라 301 × 10⁶개의 반사 트레이스를 수집하였다. 기록된 G(x₀,x,t)를 2 D Marchenko 반복 과정에 투입해 초점 함수 F를 계산하고, 식 (4)를 통해 내부 가상 수신기 G(r,x,t)를 복원한 뒤, 다시 식 (3)에 대입해 가상 음원‑수신기 응답 G(r,s,t)+G(r,s,−t)를 재구성하였다. 결과는 레이어 경계에서의 다중산란, 전파 속도 차이에 따른 굴절 및 반사 현상을 정확히 재현했으며, 유한 apertures와 2 D 모델링 한계로 인한 소규모 왜곡만을 보였다.

두 번째 사례는 1994년 노르웨이 Vøring Basin에서 수집된 고전적인 지진 반사 데이터이다. 여기서는 2 D 배경 속도 모델만을 사용해 초기 초점 함수를 추정하고, 동일한 Marchenko‑대칭화 절차를 적용하였다. 결과적으로 지진 파동의 1차 및 다중 산란 이벤트가 가상 내부 응답에 그대로 나타났으며, 기존 지진 이미지와 일관된 구조를 보여준다.

본 논문의 핵심 기여는 (1) 단일면 접근만으로도 전면 TR과 동등한 가상 음향 응답을 얻을 수 있는 수학적 프레임워크를 제공한 점, (2) Marchenko 방법을 실험적·현장 데이터에 적용해 실제 매질의 복잡한 구조 정보를 사전 모델 없이도 복원한 점, (3) 의료 초음파, 포토어쿠스틱, 지진 모니터링 등 다양한 분야에 바로 활용 가능한 범용성을 제시한 점이다. 한계점으로는 (i) 고주파 대역에서의 evanescent 파동 무시, (ii) 2 D 구현이 3 D 매질에 적용될 때 발생하는 근사오차, (iii) 초점 함수 추정에 필요한 충분한 반사 데이터 확보가 필요하다는 점을 들 수 있다. 향후 연구에서는 3 D Marchenko 알고리즘 고도화, 손실 보정 자동화, 그리고 벡터 파동(전기·탄성)으로의 확장을 목표로 할 수 있다.