지역 규모 열전도도 분석을 통한 지열 설계·유역 시뮬레이션 정확도 향상

초록

본 연구는 독일 몰라세 분지의 주요 암석군 400여 개 시료에 대해 열전도도, 밀도, 공극률, 음속을 측정하고, 단순 통계 → 교차플롯 → 전산역학적 역전산의 3단계 분석법을 제시한다. 단계별로 불확실성을 0.8 W/(m·K)에서 0.07 W/(m·K)까지 감소시켰으며, 혼합법칙과 광물 비율을 0.3 W/(m·K) 수준으로 정밀 추정했다. 결과는 지열 시추 위험 감소와 유역 규모 열·유체 흐름 모델링에 바로 적용 가능하다.

상세 분석

이 논문은 지열 개발과 유역 규모 열·유체 흐름 시뮬레이션에 필수적인 암석 물성값의 불확실성을 체계적으로 감소시키는 방법론을 제시한다. 첫 번째 단계는 전형적인 단변량 기술통계로, 평균값과 사분위수를 산출해 초기 불확실성 범위를 제시한다. 여기서는 암석군을 세부 지층으로 나누지 않을 경우 열전도도 불확실성이 0.8 W/(m·K)까지 확대될 수 있음을 보여준다. 반면, 동일 지층 내에서도 석회암·점판암·사암 등 물성적으로 유사한 하위 집단으로 구분하면 불확실성이 0.1 W/(m·K) 수준으로 크게 감소한다. 두 번째 단계는 교차플롯과 ‘quick‑look’ 모델을 활용한다. Wyllie 1956식(음속), 밀도 혼합법칙, Woodside‑Messmer 1961식(열전도도) 등을 기반으로 3‑성분(석회석, 셰일, 유체) 시스템을 가정하고, 실험 데이터가 이론 삼각형 안에 들어오는지를 검증한다. 이 과정에서 셰일의 음속·열전도도와 같은 엔드멤버 값을 실측 데이터에 맞게 보정해야 함을 강조한다. 마지막 단계는 선택된 시료군에 대해 전산역학적 역전산을 수행해 혼합법칙의 파라미터와 광물 비율을 동시에 추정한다. 역전산 결과, 석회석 매트릭스의 실제 열전도도가 문헌값(3.6 W/(m·K))보다 낮아 3.1 W/(m·K) 정도이며, 이를 반영하면 전체 암석의 열전도도 예측 오차가 0.07 W/(m·K)까지 감소한다. 또한, 음속과 밀도에 대한 예측 오차는 각각 70 m/s, 8 kg/m³ 수준으로, 현장 로그 데이터와 직접 연계하기에 충분히 정확하다. 이러한 3단계 접근법은 데이터 양과 품질에 따라 유연하게 적용 가능하며, 지역 규모 물성 데이터베이스 구축 시 효율적인 불확실성 관리 도구가 된다.