중력·자기 데이터의 Lp‑노름 공동역전 및 교차기울기 제약 통합 방법

초록

**
본 논문은 중력과 자기 데이터의 공동역전을 위해 Lp‑노름(0 ≤ p ≤ 2) 정규화와 교차기울기(cross‑gradient) 제약을 결합한 통합 프레임워크를 제시한다. 모델 파라미터와 그 기울기에 대한 L0, L1, L2 정규화를 선택적으로 적용할 수 있으며, 깊이 가중치와 하드 제약을 포함한 가중치 행렬을 손쉽게 삽입한다. 비선형 목적함수를 선형화하고 IRLS‑CG 알고리즘으로 반복 해결함으로써, 부드러운, 희소한, 혹은 블록형(단단한) 구조를 물리적으로 의미 있는 형태로 복원한다. 3‑D 합성 모델 실험을 통해 독립 역전 대비 정확도와 구조 일관성이 크게 향상됨을 보였다.

**

상세 분석

**
이 연구는 지구물리학에서 흔히 사용되는 중력·자기 데이터의 공동역전 문제를 일반화된 Lp‑노름 정규화와 교차기울기 제약이라는 두 축을 통해 통합적으로 해결한다는 점에서 의미가 크다. 먼저, 저자는 기존 단일 데이터 역전에서 사용되는 L0(희소성), L1(총변동), L2(스무딩) 정규화를 하나의 가중치 행렬 W 에 포함시켜, 파라미터 자체 혹은 그 기울기에 적용할 수 있는 유연한 프레임워크를 구축한다. 특히, (4)식에서 제시된 ε‑파라미터를 통해 정규화 강도를 조절함으로써, ε가 작을 때는 희소하고 블록형 구조를, ε가 클 때는 부드러운 연속 구조를 얻을 수 있다.

교차기울기 제약은 두 모델 m₁(밀도)와 m₂(자화율) 사이의 구조적 일관성을 강제한다. t(m)=∇m₁×∇m₂=0 조건은 두 기울기가 동일하거나 반대 방향이 되도록 하여, 경계면이 동일한 위치에 존재한다는 지질학적 가정을 수학적으로 구현한다. 이는 물리적 파라미터 간의 경험적 관계를 필요로 하지 않으며, 구조적 정보를 서로 교환하도록 만든다.

목표함수 P(α,λ) 은 데이터 적합도, 정규화 항, 교차기울기 항을 각각 가중치 α, λ 로 조절한다. 비선형성을 완화하기 위해 IRLS(Iteratively Re‑Weighted Least Squares)와 CG(Conjugate Gradient) 결합 방식을 채택했으며, 매 반복마다 선형화된 시스템을 효율적으로 풀어 계산 비용을 크게 낮춘다. 또한, 깊이 가중치 W_depth 와 하드 제약 W_h 를 포함한 블록 대각 행렬 W 을 도입해, 깊이 의존성 감소와 사전 지질 정보 반영을 동시에 구현한다.

실험에서는 두 개의 3‑D 합성 모델(하나는 부드러운 구조, 다른 하나는 블록형 구조)을 사용해 독립 역전과 비교하였다. 결과는 공동역전이 데이터 적합도(χ²)와 구조 재현 정확도(RMSE) 모두에서 우수했으며, 특히 경계면 위치와 형태가 정확히 일치하는 것을 확인했다. 이는 교차기울기 제약이 두 데이터 세트 간의 구조적 정보를 효과적으로 공유함을 증명한다.

이 프레임워크는 정규화 형태와 가중치 매개변수를 자유롭게 조합할 수 있어, 탐사 목적에 따라 ‘희소·블록형’, ‘부드러운’, 혹은 ‘혼합형’ 모델을 손쉽게 생성할 수 있다. 또한, 기존 단일 데이터 역전 코드에 최소한의 수정만으로 적용 가능하므로 실무 적용성도 높다. 향후 연구에서는 실제 현장 데이터와 비선형 전자기 효과(예: 잔류자화) 등을 포함한 확장과, 자동 파라미터 선택(α, λ, ε) 알고리즘 개발이 기대된다.

**