IGUG 그래프 이론 기반 3차원 중력 데이터 역전 패키지

초록

**
본 논문은 중력 데이터의 3차원 역전을 위해 개발된 오픈소스 MATLAB 패키지 IGUG를 소개한다. 점질량들의 좌표와 총 질량을 모델 파라미터로 삼고, 이들을 완전 가중 그래프의 정점에 대응시켜 최소 신장 트리를 구한다. 균등거리 함수로 점들의 공간 분포를 제어하고, 정규화 파라미터 λ을 통해 데이터 적합도와 안정화 항의 균형을 맞춘다. 유전 알고리즘을 이용한 전역 최적화와 새로운 λ 선택 진단 방법을 제시하며, 합성 및 실제 Mobrun 광체 데이터에 대한 실험 결과를 통해 알고리즘의 효율성과 신뢰성을 검증한다.

**

상세 분석

**
IGUG는 Bijani et al.(2015)의 그래프 이론 기반 중력 역전 아이디어를 MATLAB 환경에 구현한 확장형 도구이다. 핵심 아이디어는 지하의 균일한 물체를 다수의 단순 점질량으로 근사하고, 각 점질량의 3차원 좌표와 전체 질량을 추정 파라미터로 설정하는 것이다. 점질량들 사이의 유클리드 거리를 가중치로 하는 완전 그래프를 구성한 뒤, Kruskal 알고리즘을 적용해 최소 신장 트리(MST)를 구한다. MST의 모든 간선 길이가 가능한 한 동일하도록 하는 균등거리 함수 Θ(p)를 정규화 항으로 도입함으로써, 점들이 과도하게 군집하거나 비정상적으로 퍼지는 것을 방지하고, 물체의 전반적인 형태를 균일하게 재구성한다.

목적함수 Γ(q)=Φ(q)+λΘ(p)는 데이터 적합도 Φ(q)와 안정화 항 Θ(p)를 λ라는 정규화 파라미터로 가중합한다. Φ(q)는 관측 중력값과 모델이 예측한 중력값 사이의 가중 제곱 오차이며, 가중 행렬 W_d는 각 관측점의 잡음 분산(σ_i^2)을 역으로 이용해 정의한다. 이를 통해 잡음이 큰 데이터에 낮은 가중치를 부여하고, 잡음이 작은 데이터에 높은 신뢰도를 부여한다.

전역 최적화는 파생 정보를 필요로 하지 않는 유전 알고리즘(GA)으로 수행한다. 초기 개체군은 파라미터 범위(좌표와 질량)에 대한 물리적 제약을 만족하도록 무작위 생성되며, 각 개체는 목적함수 값을 피트니스로 갖는다. 선택, 교차, 돌연변이 연산을 반복하면서 최적 해 q_opt를 탐색한다. GA는 지역 최소에 빠지는 위험을 감소시키면서도 복잡한 비선형 탐색 공간을 효율적으로 탐색한다.

λ 선택을 위한 기존 L-curve 방법 대신, 저자들은 여러 λ 값에 대해 제한된 횟수의 GA 실행을 수행하고, Θ(p)와 Φ(q) 수렴 곡선의 선형 회귀 분석을 통해 최적 λ를 추정한다. 회귀의 결정계수(R²)와 잔차 분석을 이용해 λ가 과소·과대 정규화되는 경우를 자동으로 감지한다. 이는 계산 비용을 크게 절감하면서도 사용자에게 직관적인 λ 선택 기준을 제공한다.

실험에서는 두 개의 합성 모델(단일 구형 체와 복합 구조)과 캐나다 퀘벡 주의 Mobrun 광체 데이터를 사용하였다. 합성 실험에서 제시된 물리적 경계조건(좌표와 질량 범위)은 실제 지질 정보를 반영했으며, 결과는 원래 모델과 높은 일치도를 보였다. 실제 데이터에서는 IGUG가 기존 선형 방법보다 더 간결하고 해석 가능한 점질량 골격을 도출했으며, 추정된 총 질량과 위치는 기존 지질 조사와 일관되었다.

전반적으로 IGUG는 그래프 이론과 진화 연산을 결합한 새로운 역전 프레임워크를 제공한다. 균등거리 정규화는 물리적 균일성을 강제함으로써 비선형 역전 문제의 불안정성을 완화하고, GA 기반 전역 탐색은 초기값 의존성을 최소화한다. 또한, λ 진단 절차는 사용자가 정규화 파라미터를 선택할 때 경험적 판단에 의존하지 않도록 돕는다. 이러한 장점은 지구물리학, 광물 탐사, 환경 지질학 등 다양한 분야에서 3차원 중력 역전 도구로 활용될 가능성을 높인다.

**