중소 규모 중력 측정망을 이용한 에오트보스 행렬 요소 추정 방법

초록

본 논문은 전통적으로 토션 밸런스로 측정하던 에오트보스 행렬(Eötvös matrix) 요소를, 물리적 지표면상의 인접한 세 지점에서의 중력값만으로 추정할 수 있는 새로운 중력계 기반 방법을 제시한다. 이론적 유도와 수치 예시를 통해 제안된 방법이 충분히 정확함을 검증한다.

상세 분석

에오트보스 행렬은 중력 그라디언트(2차 미분)와 직접 연관된 3×3 대칭 행렬로, 수직 편향(deflection of the vertical), 중력 이상(gravity anomaly), 그리고 지오이드 고도 계산 등에 핵심적인 역할을 한다. 전통적으로는 토션 밸런스(torsion balance)라는 고감도 장치를 이용해 현장별로 직접 측정했으며, 이는 장비 비용·운용 난이도·측정 시간 측면에서 제한적이었다. 저자들은 이러한 제약을 극복하고자, “극소 네트워크”라 불리는 세 점의 중력값만으로 행렬 요소를 역산하는 방식을 고안했다.

핵심 아이디어는 중력장 g(x,y,z)의 2차 테일러 전개를 이용해, 기준점 P₀에서의 중력값 g₀와 인접점 P₁, P₂의 차이를 행렬 요소와 거리 벡터의 내적 형태로 표현하는 것이다. 구체적으로, Δg_i ≈ Γ·r_i (i=1,2) 로 나타내며, 여기서 Γ는 에오트보스 행렬, r_i는 P₀에서 Pi까지의 위치 벡터이다. 세 점의 Δg와 r을 알면, 6개의 독립적인 행렬 성분(대칭성으로 인해 6개) 중 5개는 선형 연립방정식으로 풀 수 있다. 남은 하나는 중력의 라플라시안(∇²g=0, 자유 공간 가정) 혹은 추가 측정점(예: 수직 방향으로의 작은 이동)으로 보강한다.

수식 전개 과정에서 저자들은 다음과 같은 가정을 명시한다. 첫째, 측정점 간 거리가 수십 미터 이하로, 지구 곡률과 대규모 지질 구조의 변화를 무시할 수 있다. 둘째, 중력장에 대한 3차 이상 항은 무시하고, 2차 항만을 고려한다. 셋째, 중력 측정 오차는 통계적으로 독립이며, 표준편차는 10⁻⁶ m/s² 수준이라고 가정한다. 이러한 전제 하에, 행렬 요소를 구하는 선형 시스템은 조건수가 낮아 수치적으로 안정적이며, 최소자승법을 적용해 오차를 최소화한다.

논문은 또한 오차 전파 분석을 수행한다. 위치 오차(Δr)와 중력 오차(Δg)가 행렬 요소에 미치는 영향을 편미분을 통해 도출하고, 실험적 시뮬레이션으로 0.5 % 이하의 상대 오차를 얻었다. 특히, E₁₁, E₂₂, E₁₂ 같은 수평 성분은 중력값 차이의 직접적인 비례관계로 인해 높은 정확도를 보였으며, 수직 성분(E₁₃, E₂₃, E₃₃)은 추가적인 수직 이동 측정이 없을 경우 약간의 편향이 발생한다는 점을 지적한다.

수치 예시에서는 가상의 지형 모델을 설정하고, 기준점에서 20 m·30 m·10 m 거리의 세 점에 대해 중력값을 가상 측정한다. 이 데이터를 이용해 행렬을 역산한 결과, 실제 행렬과 평균 상대 오차 0.38 %를 기록했다. 이는 토션 밸런스 기반 측정과 비교했을 때, 비용·시간 측면에서 큰 이점을 제공하면서도 실용적인 정확도를 유지함을 의미한다.

결론적으로, 저자들은 “극소 중력 네트워크” 접근법이 기존 장비 의존성을 크게 낮추고, 현장 조사나 위성 기반 중력 모델의 지역 보정에 유용하게 활용될 수 있음을 제시한다. 향후 연구에서는 더 큰 네트워크, 비대칭 배치, 그리고 위성 중력 데이터와의 결합을 통해 정확도와 적용 범위를 확대할 여지가 있다.