구면 파동과 스레피안 트리 변환을 이용한 지구물리 데이터 역산

초록

본 논문은 구면 위의 파동함수와 스레피안 기반 다중해상도 사전을 제시한다. 저자들은 큐브드‑스피어(cubed‑sphere) 위에 정의된 Haar·D4·D6 파동 변환을 구현하고, 지형·지진파속도 모델에 적용해 높은 압축률과 낮은 재구성 오차를 보였다. 또한, 영역‑집중성을 갖는 스레피안 함수를 이진 트리 구조로 계층화한 “스레피안 트리”를 제안하여, 대역제한·공간집중·다중해상도 특성을 동시에 만족하는 사전을 구축한다. 실험 결과, 두 방법 모두 전역 지구물리 역문제에서 희소성 기반 정규화와 결합했을 때 해의 정확도와 계산 효율을 크게 향상시킨다.

상세 분석

이 논문은 구면 데이터의 효율적 표현을 위해 두 가지 혁신적인 수학적 도구를 제시한다. 첫 번째는 큐브드‑스피어 좌표계를 이용한 구면 파동 변환이다. 기존의 구면 파동함수는 복잡한 구면 조화와 달리 지역적 지원과 다중해상도 특성을 제공하지만, 구현이 까다롭다. 저자들은 Ronchi 등(1996)이 제안한 6개의 겹치지 않거나 겹치는 “청크”(chunk)로 구면을 평면에 매핑한 뒤, 기존의 1‑D 및 2‑D 직교·비직교 파동 변환(예: D2 Haar, D4, D6)을 그대로 적용한다. 변환 전후에 경계 필터와 사전조건(preconditioning)을 도입해 경계 효과와 다항식 소거 문제를 해결한다. 실험에서는 북미 청크에 대한 지형 데이터와 670 km 깊이의 전성·전단 파동속도 모델을 변환하고, 85 %·95 % 하드 임계값을 적용해 계수의 85 % 이상을 제거하면서 재구성 오차를 5 % 이하로 유지함을 보였다. 특히 D4 파동이 블록 현상 없이 부드러운 재구성을 제공한다는 점이 강조된다. 두 번째 제안은 “스레피안 트리” 사전이다. 전통적인 스레피안 함수는 대역제한(L)과 공간집중(Region R)이라는 두 제약을 동시에 만족하지만, 전체 대역을 차지하는 전역 함수 집합으로 구성돼 다중해상도 특성이 부족하다. 저자들은 영역 R을 이진 트리 방식으로 반복적으로 절반씩 분할하고, 각 하위 영역에 대해 대역제한 L을 유지하면서 상위 n개의 스레피안 함수를 선택한다. 노드 용량 n과 트리 높이 H는 샤논 수(N)와 목표 집중도에 따라 수식적으로 결정된다(예: H≈log₂(|R|·(L+1)²/4π·n)). 이렇게 구성된 사전 D_R,L,n은 총 n·(2^{H+1}−1)개의 함수로, 서로 다른 스케일과 위치에서 거의 직교(낮은 코히어런스)를 보이며, 희소 표현과 ℓ₁ 정규화 기반 역문제에 적합하다. 실험에서는 아프리카 대륙을 대상으로 트리를 구축하고, 각 노드의 첫 번째 스레피안 함수를 시각화해 지역별 집중성을 확인한다. 또한, 트리 기반 사전이 기존 순수 스레피안 기반 역산보다 잔차 오차를 현저히 감소시키는 것을 보고한다. 전체적으로, 파동 변환은 데이터 압축과 잡음 억제에, 스레피안 트리는 지역적 해상도 조절과 희소성 촉진에 강점을 보이며, 두 기법을 결합하면 전역·지역·다중스케일을 모두 포괄하는 강력한 파라미터화 체계를 제공한다.