플루리 가우시안 기반 실시간 지질 구역 업데이트 기법

초록

본 논문은 플루리‑가우시안 시뮬레이션과 Gibbs 샘플러, 다중 데이터 동화가 가능한 Ensemble Kalman Filter‑MDA, 그리고 회전 기반 반복 가우시안화(RBIG)를 결합해 지질 구역(범주형 변수)과 연계된 등급 변수들을 near‑real‑time으로 업데이트하는 방법을 제안한다. 합성 데이터와 남호주 IOCG 매장량의 실제 드릴링 데이터를 이용한 두 사례 연구에서, 제안 알고리즘이 기존 방법에 비해 구역 경계 재현성, 등급 예측 정확도, 그리고 지질학적 현실성을 크게 향상시킴을 확인하였다.

상세 분석

이 연구는 지질 구역 모델링에서 가장 큰 난제 중 하나인 범주형 변수의 실시간 동화를 가능하게 하는 새로운 프레임워크를 제시한다. 핵심 아이디어는 기존 플루리‑가우시안 시뮬레이션(PGS)에서 사용되는 두 개 이상의 가우시안 랜덤 필드(GRF)를, Gibbs 샘플러를 통해 관측된 구역 데이터에 조건화하고, 이를 EnKF‑MDA(Ensemble Kalman Filter with Multiple Data Assimilations)와 RBIG(Rotation‑based Iterative Gaussianisation)로 연속형 변수와 결합하여 동시에 업데이트한다. 구체적으로, 1) 사전 PGS 실현을 생성하고, 2) 관측된 구역 정보를 이용해 Gibbs 샘플링으로 GRF 값을 수정한 뒤, 3) EnKF‑MDA를 적용해 다중 동화 단계에서 측정 오차를 인플레이션하여 모델‑관측 차이를 최소화한다. 여기서 RBIG는 다변량 비정규성을 가우시안 형태로 변환해 EnKF의 가우시안 가정에 부합하도록 하며, 변환·역변환 과정을 반복함으로써 복잡한 상관 구조를 보존한다. 또한, 최적의 절단 임계값을 Optuna 기반 하이퍼파라미터 튜닝으로 자동 도출하여 구역 경계의 정확성을 높인다. 합성 사례에서는 사전 실현과 실제 관측을 비교했을 때, 구역 재현률이 85 % 이상으로 향상되었으며, 등급 RMSE가 평균 30 % 감소하였다. 실제 IOCG 매장량 사례에서는 40 000개의 드릴 샘플을 활용해 5개의 구역과 3개의 연관 등급(Au, Cu, U)을 동시에 업데이트했으며, 업데이트 후 도메인 비율과 등급 통계가 관측값과 거의 일치했다. 특히, 하드 경계가 존재하는 HEM‑VOLC 전이 구역에서 기존 SIS 기반 방법이 과도하게 경계를 흐리게 만든 반면, 제안 방법은 명확한 경계와 현실적인 지질 구조를 유지했다. 계산 효율성 측면에서도, Gibbs 샘플링 1 000회와 EnKF‑MDA 4~5회 동화를 수행하는 전체 파이프라인이 30분 이내에 완료돼 현장 의사결정에 실시간 적용이 가능함을 보여준다. 그러나 두 개의 GRF에 제한된 모델링, 고차원 다변량 변환 시 RBIG 반복 횟수 증가에 따른 연산 부하, 그리고 관측 데이터의 공간적 편중이 결과에 미치는 영향 등 몇 가지 한계점도 논의된다. 전반적으로 이 연구는 범주형 지질 정보와 연속형 등급 데이터를 통합적으로 다루는 최초의 실시간 동화 프레임워크를 제공하며, 향후 복합 매장량 모델링 및 자동화된 자원 관리 시스템에 적용될 잠재력이 크다.