파동 산란을 이용한 해저 및 거칠기 동시 추정

초록

본 논문은 무한 차원 베이지안 프레임워크를 도입해 장파(종파) 산란 데이터를 활용, 해저 표면의 형태와 그 거칠기를 동시에 추정한다. 해저를 통계적 등방성(Stationary, Isotropic) 랜덤 함수로 가정하고, 분수 미분가능성(Fractional Differentiability) 지표를 거칠기 파라미터로 정의한다. 제안된 알고리즘은 사후분포를 효율적으로 탐색해 해저 프로파일과 거칠기 추정값을 제공하며, 불확실성 정량화도 수행한다. 수치 실험을 통해 대규모 해저 탐사에 적용 가능한 정확도와 안정성을 입증하였다.

상세 분석

이 연구는 해저 탐사를 위한 역산 문제를 무한 차원 베이지안 관점에서 재구성함으로써 기존의 유한 차원 접근법이 갖는 차원 저주와 수치 불안정성을 근본적으로 회피한다. 핵심 가정은 해저 표면이 통계적 등방성을 띤 랜덤 함수라는 점이다. 등방성 가정은 지역적인 파동 측정만으로도 전역적인 스펙트럼 특성을 추정할 수 있게 하며, 이는 기존의 주파수별 파워 스펙트럼 추정 방식보다 더 일반적인 프레임워크를 제공한다.

거칠기 파라미터는 전통적인 스펙트럼 지수나 자기상관 길이 대신, 함수의 분수 차수(α)로 정의된다. 분수 차수는 파워 스펙트럼의 고주파 감쇠율을 직접적으로 나타내며, 노이즈와 불완전한 측정에 대해 강인한 특성을 가진다. 논문은 이 분수 차수를 베이지안 사전분포에 포함시켜, 사후분포가 측정 데이터에 대해 Lipschitz 연속성을 유지함을 수학적으로 증명한다.

전방 모델은 2차원 영역에서의 스칼라 파동 방정식을 사용해 장파(p‑wave)의 전파를 기술한다. 물과 암석의 물성은 각각 깊이에 따라 선형적으로 변하는 밀도와 라메 파라미터로 모델링하고, 해저 인터페이스 h(x)를 통해 물성과 파동 전파 속도가 구분된다. 경계조건은 해수면에서의 Neumann 조건과 개방 해양을 모사하는 비반사 경계조건을 적용해 실제 해양 환경을 근사한다.

수치 해석은 FEM–FD 결합 방식으로 구현된다. 공간을 유한요소로 이산화하고 시간은 중앙 차분을 사용해 안정적인 전진 시뮬레이션을 수행한다. 이때 해저 인터페이스는 레벨셋 형태로 파라미터화되어, 인터페이스 자체가 무한 차원 베이지안 변수로 취급된다.

베이지안 사후분포 탐색은 사전분포(등방성 가우시안 랜덤 필드)와 관측 모델(선형/비선형 연산자 G) 그리고 가우시안 노이즈를 결합한 형태이며, 마르코프 체인 몬테 카를로(MCMC) 혹은 변분 베이지안 방법을 변형한 사전‑후처리 알고리즘을 사용한다. 논문은 특히 사후분포의 지역적 Lipschitz 연속성을 이용해 효율적인 사전‑후처리 샘플링 전략을 설계하고, 이를 통해 추정된 해저 프로파일과 거칠기 파라미터의 불확실성 영역을 시각화한다.

실험에서는 합성 데이터와 실제 해저 시뮬레이션 데이터를 이용해 두 가지 시나리오(사전 분포 내 샘플, 사전 분포 외 샘플)를 검증한다. 결과는 전통적인 최적화 기반 역산 방법에 비해 평균 제곱 오차가 현저히 낮으며, 특히 거칠기 추정에서 높은 정확도와 낮은 분산을 보인다. 또한 불확실성 분석을 통해 데이터가 부족하거나 노이즈가 큰 영역에서 추정값이 크게 편향될 수 있음을 확인하고, 이를 기반으로 측정 설계(센서 배치, 주파수 선택) 최적화에 활용할 수 있음을 제시한다.

전반적으로 이 논문은 무한 차원 베이지안 프레임워크와 분수 미분가능성 개념을 결합해 해저 탐사의 핵심 난제인 “다중 해저 구성 → 유사 측정” 문제를 효과적으로 해결한다는 점에서 학술적·실용적 의의가 크다.