중력 단층 촬영을 위한 타원체 모델 활용

초록

본 논문은 매장체를 수직 바 형태로 분할해 직접 중력 문제를 해결하고, 역문제에서는 균일한 회전 타원체(구면체)로 근사한다. 변환된 구면체 중력식과 제약 조건을 적용한 티크노프 정규화를 통해 해의 유일성과 안정성을 확보한다. 초기 깊이·질량 추정에 사용되는 Bulakh 알고리즘을 개선하고, 두 개·다섯 개 체의 수치 사례를 통해 방법의 유효성을 검증한다.

상세 분석

이 연구는 중력 탐사에서 흔히 마주치는 비선형·불안정한 역문제를 해결하기 위해 두 단계 접근법을 제시한다. 첫 번째 단계는 직접 문제로, 매장된 물체를 동일한 단면적을 가진 수직 바(바렐)들로 세분화한다. 각 바에 대한 중력 기여를 선형적으로 합산함으로써 복잡한 지오메트리를 간단히 표현하고, 계산 효율성을 크게 높인다. 두 번째 단계는 역문제로, 전체 매장을 하나의 균일한 회전 타원체(구면체)로 모델링한다. 기존에 알려진 구면체의 z‑성분 중력식은 적분 형태가 복잡해 실용적 적용에 제약이 있었으나, 저자들은 이를 폐쇄형 식으로 변형하고, 매개변수(반경, 장축, 밀도·질량)를 명시적으로 나타냈다.

역문제 해석에서는 티크노프 정규화(정규화된 최소제곱) 방식을 채택하고, 물리적·지질학적 제약(예: 깊이와 장축 비율의 범위)을 추가한다. 이러한 제약은 해 공간을 제한해 ill‑posed 문제를 유일하고 안정적인 해로 전환한다. 특히, 정규화 파라미터 선택을 자동화하기 위해 L‑곡선 분석과 교차 검증을 병행한다.

초기 추정값을 제공하는 Bulakh 알고리즘은 기존에 깊이와 질량을 단순히 비례 관계로 추정했으나, 여기서는 매장체의 비대칭성을 반영하도록 수정하였다. 구면체의 장축·단축 비율을 사전 추정하고, 이를 기반으로 깊이·질량 초기값을 도출한다. 이렇게 얻은 초기값은 비선형 최적화 과정에서 수렴 속도를 크게 향상시킨다.

수치 실험에서는 두 개와 다섯 개의 독립적인 매장체를 가정하고, 각각에 대해 직접·역문제 시뮬레이션을 수행했다. 결과는 모델 파라미터가 실제값과 5 % 이내의 오차로 복원되었으며, 잡음이 섞인 데이터에서도 안정적인 재구성이 가능함을 보여준다. 특히, 다중 매장체 상황에서 서로 간의 중력 신호 간섭을 효과적으로 분리해낼 수 있었다.

이 논문의 핵심 기여는 (1) 직접 문제를 바 기반으로 단순화한 계산 효율성, (2) 구면체 모델을 이용한 역문제의 해석적 식 도출, (3) 제약 기반 티크노프 정규화로 해의 유일성·안정성 확보, (4) Bulakh 알고리즘의 개선을 통한 초기값 최적화이며, 이는 중력 단층 촬영(gravimetric tomography) 분야에 새로운 해석 프레임워크를 제공한다.