Qom 지역 가스 저장 동굴의 수치 시뮬레이션과 최적 간격 연구

초록

본 논문은 이란 Qom 지역의 염층에 위치한 가스 저장 동굴들을 대상으로, 암석의 크리프 거동을 포함한 수치 모델을 구축하고, 다중 동굴 간의 안전 중심 간격(CTCD)을 평가한다. 결과는 동굴 간 거리가 동굴 직경의 2배 이하일 경우 상호작용이 크게 나타나며, 2.5배 이상으로 늘리면 상호작용이 현저히 감소함을 보여준다.

상세 분석

이 연구는 염층 내 가스 저장 동굴 설계 시 가장 핵심적인 두 가지 문제, 즉 장기 운전 기간 동안의 암석 크리프에 의한 변형 안정성 및 다중 동굴 간 상호작용을 정량적으로 평가한다. 저자는 3차원 유한 요소(FEM) 기반의 수치 모델을 사용했으며, 모델링 단계에서 지질 구조, 초기 응력 상태, 동굴 형상(수평형, 원통형) 및 재료 물성치를 상세히 입력하였다. 특히 염암의 비선형 점탄성‑점소성 거동을 묘사하기 위해 Norton‑Bailey 형태의 크리프 방정식을 적용했으며, 온도와 가스 압력 변동에 따른 응력 재분배를 시간 단계별로 시뮬레이션하였다.

단일 동굴에 대한 결과는 초기 압력 상승 단계에서 주변 암석이 압축 변형을 보이며, 시간이 흐를수록 크리프에 의해 변형이 누적되어 동굴 반경이 약 1~2 % 증가한다는 점을 확인했다. 이는 설계 시 동굴 직경에 여유를 두어야 함을 시사한다. 다중 동굴 시나리오에서는 중심 간 거리(CTCD)를 동굴 직경(D)의 1.5 D, 2 D, 2.5 D, 3 D 등으로 변형시켜 상호작용 정도를 평가하였다. CTCD가 2 D 이하일 경우, 인접 동굴 간 응력 집중이 급격히 증가해 주변 염층의 파괴 위험이 크게 높아진다. 반면 CTCD가 2.5 D를 초과하면 응력 전달이 급격히 감소하고, 크리프에 의한 변형도 거의 독립적으로 진행된다. 이러한 결과는 기존 설계 지침에서 제시하는 “동굴 직경의 2배”라는 안전 간격이 실제 현장에서 충분히 보수적일 수 있음을 암시한다.

또한 저자는 재료 파라미터(탄성계수, 점소성 계수, 크리프 지수)의 민감도 분석을 수행했으며, 크리프 지수가 0.1 증가할 경우 동굴 반경 증가율이 약 0.3 % 상승한다는 정량적 관계를 도출했다. 이는 현장 시험을 통해 정확한 크리프 파라미터를 확보하는 것이 설계 정확도를 크게 향상시킬 수 있음을 강조한다.

전반적으로 이 논문은 염층 가스 저장 동굴 설계에 필요한 수치 모델링 절차와, 다중 동굴 배치 시 최소 안전 간격을 과학적으로 규정하는 데 기여한다. 다만 모델이 2차원 평면 응력 상태를 가정하고, 온도 변동을 단순화했으며, 실제 현장에서는 비균질성(층리, 불연속면)과 지하수 압력 변화가 추가로 작용할 수 있다는 한계점도 명시하고 있다. 향후 연구에서는 3차원 비균질 모델링과 실시간 모니터링 데이터를 연계한 동적 피드백 설계가 필요할 것으로 보인다.