구면 웨이브렛 기반 전지구 단층 모델 해석과 희소성

초록

본 논문은 구면 웨이브렛을 이용해 전지구 지진단층 모델을 다중해상도로 표현하고, 파형 데이터와 결합한 역문제에서 ℓ₁ 정규화를 통해 희소한 해를 얻는 방법을 제시한다. 큐브드 스피어를 이용한 6개의 청크와 경계 처리를 포함한 두 종류의 이산 웨이브렛 변환을 설계하고, 실제 모델에 적용해 계수의 스케일별 통계와 희소성을 분석한다. 합성 실험을 통해 제안 방법이 기존 방식보다 높은 해상도와 계산 효율성을 제공함을 입증한다.

상세 분석

이 연구는 전지구 지진단층 역문제에 적용 가능한 구면 웨이브렛 기반 프레임워크를 체계적으로 구축한다. 먼저, Ronchi 등(1996)의 큐브드 스피어를 6개의 청크로 분할하고, 각 청크 내부를 직교 좌표(ξ, η)로 매핑한다. 이 좌표계는 구면 전체에 거의 균일한 공간 해상도를 제공하지만 청크 경계에서 비연속성이 발생한다는 점이 핵심 난점이다. 이를 해결하기 위해 저자들은 “interval wavelet” 기법을 도입, 경계 근처에서는 특수 필터와 사전조건(preconditioner)을 적용해 다중해상도 기반의 다중스케일 분석이 가능한 이산 웨이브렛 변환을 구현한다. 두 가지 변환은 (1) Daubechies의 4‑tap 정규 직교 웨이브렛(D4)과 (2) Cohen‑Daubechies‑Feauveau(CDF) 4‑2 바이오쏘리얼 웨이브렛이다. D4는 대칭성이 없지만 완전 직교성을 제공하고, CDF는 분석·합성 쌍을 통해 대칭성을 유지하면서도 높은 재구성 정확도를 보인다.

논문은 두 개의 공개된 전지구 속도 모델(예: S20RTS, GyPSum)에서 각 청크와 깊이 층별로 웨이브렛 계수를 계산하고, 계수 절댓값을 임계값으로 차단(threshold)하여 재구성 오류와 압축 비율을 평가한다. 결과는 고스케일(큰 파장) 계수가 전체 에너지의 대부분을 차지하고, 저스케일(작은 파장) 계수는 희소하게 분포함을 보여준다. 특히, CDF 기반 변환은 경계 인접 영역에서 발생하는 인공적인 진동을 최소화하며, D4 변환은 경계 근처에서 약간의 아티팩트를 보인다.

역문제 해결 단계에서는 데이터 적합도 ℓ₂ 노름과 웨이브렛 계수 ℓ₁ 노름의 가중합을 최소화하는 최적화 문제를 설정한다. 이때, 선형 연산자는 전파 민감도 커널(유한 주파수 커널)과 웨이브렛 변환 행렬의 곱으로 구성되며, 빠른 행렬-벡터 곱을 위해 병렬 FFT와 멀티코어 GPU 구현이 가능하도록 설계되었다. 최적화 알고리즘은 ISTA/FISTA와 같은 가속된 경사 하강법을 사용해 수백 번의 반복만에 수렴한다. 합성 실험에서는 잡음이 포함된 가상의 데이터에 대해 전통적인 구면 조화 기반 역문제와 비교했을 때, 제안된 웨이브렛 기반 방법이 동일한 데이터 적합도에서 약 30 % 이상의 파라미터 수 감소와 더 선명한 구조 복원을 달성했다.

이러한 결과는 전지구 지진단층 모델링에서 공간적 해상도를 위치에 따라 가변적으로 조정하고, 물리적으로 의미 있는 구조를 최소한의 자유도(희소성)로 표현할 수 있음을 시사한다. 또한, 웨이브렛 기반 파라미터화는 정규화 선택에 따라 모델의 스무딩 정도를 직관적으로 제어할 수 있어, 지구 물리학적 해석과 데이터 기반 모델 검증 사이의 간극을 메우는 유용한 도구가 된다.