Aims & Scope

주입 지역에서 Paaratte 층은 Skull Creek Mudstone 층을 덮고 있는 삼각주 환경의 퇴적체이다(그림 2, 4). 감마선 검층과 시추코어 주상도에서 삼각주 퇴적체는 최하부에서 상부로 가면서 전반적으로 상향조립화하는 것이 관찰된다. 이는 광역적으로 분포하는 이암층 위에 해안선이 바다쪽으로 이동하며 전진하는 삼각주(prograding delta) 환경의 발달을 지시한다.

CRC-2 공의 주입 대상 지층에서 Paaratte 층은 크게 5 개의 퇴적상으로 구성된다(그림 4). Paaratte 층의 하부는 삼각주전면(delta front), 말단사주(distal mouth bar) 환경에서 퇴적된 이질암이 우세하고, 상부로 가면서 근위사주(proximal mouth bar), 분류하도(distributary channel) 환경에서 퇴적된 사질암이 우세하게 발달한다(그림 4). 주입 대상 지층은 상부에 말단사주 환경에서 퇴적된 이암층에 의해 피복된다(그림 4). 이는 퇴적 중심이 다른 곳으로 이동하였거나, 갑작스런 구조적 침강에 따라 침식 기준면이 깊어지면서 일어난 해침에 의한 것으로 해석된다. 해침이 일어나면서 퇴적된 이암층은 광역적인 덮개층을 형성하였다.

Fig. 4. 
Well log correlation of reservoir and seal intervals in the Paaratte Formation. For well location, see Figure 1.

대상 지층의 저장능력과 주입능력을 평가하기 위해 CRC-1 코어시료에서 측정한 공극률과 투과도, 물리검층으로부터 계산한 공극률과 투과도를 검토하였다(그림 5, 6; 표 1). 저류층의 공극률은 평균 25∼30%로 매우 양호하며, 투과도도 수 다르시(Darcy)에 이를 정도로 양호하다. 공극률과 투과도는 퇴적상별로 큰 차이를 보인다(그림 5; 표 1). 저류층을 구성하는 주요 퇴적상인 근위사주와 분류하도의 시료는 공극률과 투과도가 뚜렷하게 높다. 반면에 덮개층의 주요 퇴적상인 말단사주, 삼각주전면 시료는 매우 낮게 측정되었다. 속성작용에 의해 침전된 교질물(cement section) 퇴적상은 공극률과 투과도가 현저하게 낮다. CRC-1 코어시료 분석에 의하면 대개 공극률이 20% 이하이고, 투과도가 4 mD 이하인 경우 비저류암으로 구분된다(표 2).

Fig. 5. 
Cross plot of porosity versus permeability for different types of facies, (a) proximal mouth bar, (b) distributary channel, (c) distal mouth bar, and (d) cement section in the Paaratte Formation.

CRC-2 공의 주입 구간을 포함하는 저류층의 물리검층 자료로부터 계산한 공극률과 투과도 값은 코어에서 측정한 값과 대개 유사하며, 퇴적상에 따른 분포 양상도 매우 유사하다(그림 6). 투과도 값은 공극률과 상관성이 매우 높으며, 퇴적상에 따라 분포 범위가 한정된다.

Fig. 6. 
Cross plot of porosity versus permeability derived from well log.

주입공을 포함하는 모델 영역의 총 면적은 325,000 m² 이다. 단위격자(cell)의 크기는 x, y 방향으로 15 m × 15 m 이고, 남쪽 경계단층의 방향(북서-남동)을 고려하여 약 70°로 기울어지게 배열하였다. 단위격자의 방향을 주요 지질구조인 단층 방향과 일치하게 배열하여 지질구조에 의한 이방성을 반영하고, 모델 경계의 바깥에 놓이는 노드(node)의 개수를 최소화할 수 있다(Anderson and Woessner, 1992).

탄성파 자료를 이용하여 저류층의 하부 경계면에 해당되는 반사면(Top Skull Creek Mudstone)을 대비하고 분포를 맵핑하여 지질격자 모델의 수평 경계면으로 삽입하였다(그림 7). 지질격자 모델에서 저류층과 덮개층은 CRC-1, CRC-2, Naylor-1 공과 인근의 시추공으로부터 계산한 두께를 이용하여 각각 삽입하였다(그림 7). 주입구간을 포함하는 저류층의 두께는 약 27 m, 상부 덮개층의 두께는 약 3 m 이다.

Fig. 7. 
Simplified structural grid containing multi-layer reservoir and single-layer seal. Note that each layer follows the lower boundary of reservoir.

저류층 내에 발달하는 층형(stratal pattern)을 묘사하기 위해 레이어를 생성하였다. 레이어는 0.5 m 두께로 일정하고, 하부 경계면과 평행하게 설정하였다(그림 7). 해침 퇴적계의 경우 내부에 발달하는 지층들이 하부 경계면을 따라 윗걸침(onlap)하며 채우는 층서적 관계가 있다. 지질격자 모델의 레이어를 상부 경계면과 평행하게 차례대로 배열하면 해침 퇴적계의 위걸침하는 층서적 관계를 모델에 반영할 수 있다. 반면에 고해수면 퇴적계는 상부 경계면을 따라 지층이 상향종식(toplap or truncation)하기 때문에 하부 경계면과 평행하게 레이어를 구성하는 것이 적절하다. 주입구간의 저류층은 삼각주 환경에서 전진구축하며 퇴적된 사암층이 우세하기 때문에 하부 경계면과 평행하게 레이어를 생성하여 저류층내 층서적 관계를 묘사하였다(그림 7).

수직 레이어의 두께는 저류층내 분포하는 퇴적상 두께의 비균질성을 고려하여 결정해야 한다. 사암 퇴적상(proximal mouth bar, distributary channel)은 수직 방향으로 누적되어 나타나기 때문에 수 m 두께로 비교적 두껍고, 이암 퇴적상(distal mouth bar, delta front, cement section)은 수 cm 혹은 수 m 두께로 두께가 다양하다. 이암 퇴적상이 저류층 내에 협재할 경우 주입 이산화탄소의 흐름을 방해하기 때문에 이들의 수직적 분포를 묘사할 수 있도록 수직 레이어의 두께를 산정하였다. 시추공에서 구분한 퇴적상에 대해 수직 베리오그램 분석을 수행하여 각 퇴적상의 수직 상관거리를 유추하였고, 코어에서 관찰한 각 퇴적상의 수직적 두께 변화와 주기성 등을 고려하여 0.5 m 두께의 레이어를 생성하였다. 덮개층은 말단사주 환경에서 퇴적된 이암으로 연속성이 좋고, 균질한 지층이기 때문에 수직 레이어로 세분화하지 않았다. 지질격자 모델의 전체 단위격자의 개수는 118,440 이다.

세 개 시추공 위치에서 지질격자 모델에 퇴적상, 공극률, 투과도 값을 할당하기 위해 업스케일을 수행하였다. 코어에서 구분한 교질물 퇴적상은 수직 레이어 두께인 0.5 m 보다 얇은 층으로 존재할 경우 업스케일한 셀에 표현되지 못한다. 이를 보완하기 위해 교질물 퇴적상은 업스케일 과정에서 다른 퇴적상들보다 높은 가중치를 주었다. 삼각주전면 퇴적상은 저류층내에 아주 제한적으로 나타나기 때문에 업스케일 과정에서 제외하고 모델에 반영하지 않았다. 업스케일한 퇴적상은 코어에서 구분한 원래 퇴적상과 유사한 비율로 나타난다(그림 8).

Fig. 8. 
Histogram showing proportion of upscaled and original facies. The upscaled facies preserve the proportions of the original faceis. For facies codes, see Table 1.

덮개층에 해당되는 구간은 말단사주 퇴적상 하나로 구성된다. 주입 저류층의 퇴적상 분포 모델링은 세 개 시추공에서 업스케일한 퇴적상을 입력자료로 활용하였고, 각 퇴적상의 수직 베리오그램 분석 결과를 기반으로 순차지표시뮬레이션을 이용하여 수행되었다. 즉, 시추공 위치에서 업스케일된 퇴적상은 원래 퇴적상과 비율이 거의 동일하고, 시추공 사이에서의 퇴적상 분포는 순차지표시뮬레이션 기법에 의해 예측되었다. 퇴적상 분포의 이방성(major and minor direction)과 연속되는 상관거리(range)는 코어, 물리검층 자료로부터 수평 베리오그램 분석하여 시뮬레이션에 반영할 수 있다. 연구지역은 시추공 수가 부족하여 각 퇴적상의 수평 베리오그램 분석이 어렵기 때문에 노두에 드러난 퇴적층에서 분석한 결과를 간접적으로 활용하였다(표 3). Paaratte 층과 퇴적환경이 유사한 고기의 지층과 현생의 퇴적환경에 대한 연구결과를 검토하여 각 퇴적상 분포의 수평 상관거리 및 이방성을 다음과 같이 가정하였다.

삼각주 환경의 분류하도 사암은 두께와 폭의 비율이 1:50 정도로 알려져 있다(Tye, 2004). 코어에서 하도 사암의 최대 두께가 대개 4 m 임을 고려할 때 하도의 폭이 200 m 이하일 것으로 유추할 수 있다. 따라서, 분류하도의 길이는 최대 800 m 정도 연장성을 가지며 분포할 것으로 여겨진다(Gani and Bhattacharya, 2007).

말단사주와 삼각주전면 퇴적상은 마하캄 삼각주와 같은 현생 환경에서 수 십 km 규모로 분포하기 때문에 모델링 영역내에서 연속적으로 분포할 가능성이 높다. 실제로 덮개층을 구성하는 말단사주 퇴적상은 모든 시추공에서 연속적으로 나타나기 때문에 이방성을 고려하지 않았다.

저류층 내에 분포하는 교질물 퇴적상은 속성작용에 의해 2차적으로 형성된 산물이기 때문에 연속성과 규모를 예측하기 어렵다. Paaratte 층과 유사한 퇴적환경에서 퇴적된 암상의 노두 관찰에 의하면 교질물 퇴적상은 수평 방향으로 약 30 m, 수직 방향으로 2.5 m의 상관거리를 가지며(Dutton et al., 2002), 지하투과레이다(Ground Penetrating Radar) 자료를 이용한 3차원 분포 예측에 의하면, 최대 두께가 1.0∼1.4 m 이다(Lee et al., 2007). 시추코어 관찰에서 수평 상관거리에 대한 분석이 불가능하기 때문에 노두 분석 결과를 바탕으로 수평 상관거리를 30 m로 가정하였다. 시추코어에서 교질물 퇴적상의 두께와 출현빈도가 매우 다양하여 수직 상관거리는 1 m로 작게 가정하였으며, 이방성은 고려하지 않았다.

Paaratte 삼각주 퇴적시스템의 발달은 퇴적 동시성 성장 단층의 발달과 관련이 있다. 북서-남동 방향의 성장 단층은 삼각주 퇴적시스템이 단층 방향과 수직인 북동-남서 방향으로 발달하였음을 지시한다. Paaratte 층에 누적되어 있는 삼각주 퇴적체는 주로 북동-남서 방향으로 전진과 후퇴를 반복하면서 형성되었을 것이다. 따라서 분류하도, 근위사주 퇴적상의 분포는 주로 북동-남서 방향으로 우세하게 발달하였을 것으로 추측된다.

퇴적상 모델링 결과, 분류하도와 근위사주 퇴적상과 같은 저류암은 북동-남서 방향으로 발달하는 이방성과 연속성이 적절하게 묘사되었고, 교질물 퇴적상은 불연속으로 분포한다(그림 9a). 말단사주 퇴적상은 저류층 하부 구간에 우세하게 분포하고 상부로 가면서 분류하도와 근위사주 퇴적상이 우세한 경향이 나타난다(그림 10a). 이러한 경향은 물리검층과 코어자료 분석에서 인지된 상향 조립화 경향이 모델에 반영되어 있음을 보여준다. 퇴적상 모델에서의 가장 큰 불확실성은 교질물 퇴적상의 분포 양상이다. 모델에서 입력한 수평 상관거리가 작고, 이방성을 고려하지 않았기 때문에 교질물 퇴적상의 분포가 단속적이다 (그림 9, 10). 주입 지층의 주요 저류암인 분류하도와 근위사주는 투과도가 매우 양호하기 때문에, 주입 이산화탄소의 흐름은 교질물 퇴적상의 횡적 연속성과 방향성에 의해 주로 영향을 받을 것으로 여겨진다. 속성작용과 탄산염 교질물 침전과 용해 기작을 이해하고 분포 특성을 파악하기 위해서는 고퇴적환경의 수직적, 수평적 변화에 대한 해석이 필요하다.

Fig. 9. 
(a) Facies model built by variogram-based, Sequential Indicator Simulation method, (b) Porosity model guided by the facies distribution (used as secondary data), (c) Permeability model co-krigged with the porosity. Sequential Gaussian Simulation method was used to model porosity and permeability.

Fig. 10. 
Cross-sectional views of (a) facies, (b) porosity, and (c) permeability models.

시추공 사이의 공간에서 공극률 분포를 예측할 만한 직접적인 자료가 불충분하고, 공극률의 공간적 분포를 간접적으로 지시할 수 있는 3차원 탄성파 속성(attribute) 분석 결과가 불량하기 때문에 세 개 시추공에서 업스케일한 공극률 자료를 이용하여 공간적 분포를 예측하였다. 물리검층 자료에서 계산한 공극률은 각 퇴적상에 따라 분포가 상이하다(그림 6). 즉, 공극률의 공간적 분포는 저류층내 퇴적상 분포 경향을 반영한다. 따라서, 공극률 분포를 예측하기 위해 퇴적상 분포 모델링 결과를 반영하는 조건부 시뮬레이션을 수행하였고, 순차가우시안 시뮬레이션 기법을 이용하였다.

순차가우시안 시뮬레이션 기법은 정규분포를 가지는 자료에 대해 지구통계적 분포 예측이 유용하다. 시추공 위치에서 공극률 값들이 정규분포를 가지도록 변환(normal score transformation)하고 순차가우시안 시뮬레이션을 이용하여 모델링을 수행한 후 다시 자료를 역변환(back transformation)하여 공극률 모델 결과를 얻었다. 순차가우시안 시뮬레이션은 시추공 사이의 공간에서 공극률을 무작위로 계산하지만 시추공 위치에서 공극률 값이 가지는 통계학적 특성은 최종 모델의 전체 공극률의 통계학적 특성과 일치한다. 시추공 개수가 충분하지 않기 때문에 수평 베리오그램 분석을 통한 공극률 분포의 이방성과 상관거리는 반영하지 않았다. 그럼에도 불구하고, 각 퇴적상별로 공극률의 분포를 예측하였기 때문에 최종 공극률 모델은 퇴적상의 공간적 비균질성을 잘 반영하고 있다(그림 9b, 10b). 즉, 고공극률은 근위사주, 분류하도 퇴적상과 같이 사암이 우세한 퇴적상에서 뚜렷하게 나타난다.

투과도 모델링은 순차가우시안 시뮬레이션 기법을 사용하였다. 투과도 변화와 공극률 변화의 상관성이 매우 높기 때문에 선행된 공극률 모델링 결과를 공동크리깅(co-kriging)으로 반영하였다. 또한, 공극률 모델링과 마찬가지로 퇴적상 분포 모델을 조건부 자료로 활용하였다. 따라서, 투과도의 공간적 분포는 불균질한 퇴적상의 분포 경향을 반영하고, 예측된 공극률의 분포 경향과 유사하다(그림 9c, 10c).

저류층 하부에 넓게 분포하는 말단사주 퇴적상의 분포로 인해 저류층 하부는 투과도가 낮고, 상부로 가면서 투과도가 높아진다(그림 10c). 저류층 상부 구간에서 이산화탄소를 주입할 경우 분류하도 퇴적상의 투과도가 높기 때문에 이들의 수직 및 수평적 연속성에 따라 이산화탄소의 이동 양상이 달라질 것으로 여겨진다. 분류하도 퇴적상은 대개 근위사주 퇴적상과 교호하며, 두 퇴적상의 투과도는 모두 양호하기 때문에 분류하도 퇴적상의 연속성과 이방성이 불확실하더라도 이산화탄소의 이동에 큰 영향은 없을 것으로 여겨진다. 반면에 낮은 투과도를 가지는 교질물 퇴적상의 분포가 이산화탄소 흐름을 지연시키거나 방해하는 주요한 요인이 될 것이다.

CRC-2 공 주입 지층을 구성하는 퇴적상을 대상으로 수은주입모세관압(mercury injection capillary pressure, MICP) 시험 방법으로 침투압(entry pressure)을 측정하였다(Daniel et al., 2011). 이 때, 이산화탄소가 주입되는 구간의 압력과 온도, 그리고 이미 공극을 점유하고 있는 물의 염도를 고려하였고, 덮개암이라는 가정 하에 이산화탄소가 투과하기 위해서 요구되는 압력을 나타내는 이산화탄소 칼럼(column)의 높이를 측정하였다.

이산화탄소 칼럼의 높이는 교질물 퇴적상에서 1.2∼197 m의 넓은 범위를 구성하는 반면, 저류암의 경우 1.2∼10.3 m로 상대적으로 좁은 범위를 구성한다. 이러한 두 개의 다른 범위를 바탕으로 앞의 정적 모델링에서 선정한 다섯 개 퇴적상에 대한 침투압을 유추하였다.

상대투과도에 대한 모델은 석유분야에서 흔히 쓰는 Brooks-Corey의 모델(Brooks and Corey, 1964)을 따랐다. 어느 한 유체의 상대투과도는 그 유체의 실질적인 포화도에 의한 함수로 습윤상(wetting phase) 유체의 유효포화도(Swe)는 다음 식과 같이 정의된다.

수식 1에서 보는 봐와 같이 침투압과 Brooks-Corey 모델 상수 λ, 그리고 퇴적상에 따른 물의 잔류포화도(residual saturation) 등이 앞서 소개한 Daniel 등이 측정한 이산화탄소의 컬럼 높이에서 유추하여 표 4와 5에 정리하였다. 다섯 개 퇴적상들을 대상으로 덮개암의 효율과 상대투과도를 Brooks-Corey 모델로 표현함으로써 다상유체의 특성이 규명되었다.

이산화탄소 주입과 거동 시뮬레이션을 위해 CRC-2 공의 주입구간에 150일 동안 159 톤/일의 이산화탄소를 주입하는 것으로 가정하였다. 이산화탄소 주입 전의 초기조건은 저류층 내에 염수가 공극을 점유하고 있고, 주어진 심도에서 정수압(hydrostatic) 조건을 만족하는 것으로 가정하였다. 이는 대수층에 존재하는 물이 외부의 영향이 없을 경우 움직이지 않고 정적으로 존재하는 상태를 가정한다. 그리고 경계조건의 경우에는 가장 적은 수의 노드만을 정해서 초기조건과 마찬가지로 정수압 조건을 물의 압력과 상응하는 모세관압을 계산해서 설정하였다.

앞서 가정한 동적 모델링의 초기조건과 경계조건, 약 24,000 톤 정도의 이산화탄소를 주입하기 위한 주입조건과 주입구간, 그리고 정적 모델링에서 구현한 다섯 개의 퇴적상에 대한 Brooks-Corey 모델의 재해석을 완성하면, 예비적인 이산화탄소 주입에 대한 최소한의 모델링 조건이 구성된다. 이에 기초해서 총 92,205개의 노드와 총 86,039개의 요소로 이루어진 동적모델링 도메인을 구성했다. 0.5일 시간간격으로 150일 동안 이산화탄소를 주입하는 시뮬레이션을 수행하였다.

다섯 개의 퇴적상으로 구성된 비균질한 지질 모델을 기반으로 이산화탄소를 주입한 경우, 주입구간을 중심으로 분포한 퇴적상들의 비균질한 물성을 반영하여 이산화탄소 플룸이 확장되는 점을 관찰할 수 있다. 이를 3차원으로 도시하기 위해서 이산화탄소 포화도 분포를 그림 11에 도시했다. 이 경우에는 물성에 따라 다상유체유동에서 보여주는 단속적인(discrete) 구분이 일어남을 잘 관찰할 수 있다. 즉, 비균질한 지질 모델에서 이산화탄소 플룸은 주입 구간을 따라 규모와 형태가 다양하고, 횡적 이동의 정도가 상이하며, 포화도의 차이가 관찰된다. 이산화탄소 포화도는 주입 구간의 하부에서 상대적으로 높고, 상부로 가면서 작아진다. 이산화탄소 플룸은 전반적으로 하부에서 넓게 퍼지고 상부에서 분포 범위가 좁은 피라미드형으로 발달한다(그림 11). 레이어 규모에서 일어나는 현상을 살펴보면, 투과도가 좋은 레이어에 주입된 이산화탄소는 횡적 이동이 일어나기 시작한다. 상부에 불투과성 레이어가 존재할 경우, 불투과성 레이어의 모세관압이 주입 이산화탄소의 상태압(phase pressure)보다 크기 때문에 이산화탄소의 수직 이동이 제한되며, 하부 레이어에서 이산화탄소는 계속 누적되고 수평방향으로 넓게 퍼지게 된다. 주입이 계속되는 동안 이산화탄소의 상태압이 침투압을 극복하게 되면 상부의 불투과성 레이어로 침투가 일어나고, 다시 수평 방향으로의 횡적 이동이 일어나는 현상이 반복적으로 나타난다. 이러한 현상은 저류층 내에 분포하는 불투과성 퇴적상의 연속성과 분포에 영향을 받는 것으로 해석된다. 주입 조건에서 이산화탄소 플룸은 덮개암에 의해 분리되어 저장되는 모습을 관찰할 수 있다.

Fig. 11. 
Simulated CO₂ saturation and migration for heterogeneous facies model. Simulation time is 50 days. (a) CO₂ plume with saturation larger than 0.01. (b) Cross-sectional view of CO₂ saturation. (c) Three dimensional view of CO₂ saturation. (d) CO₂ plume sealed by seal formation.

저류암과 비저류암으로 구분한 단순한 지질 모델에 이산화탄소를 주입한 경우, 주입구간을 중심으로 분포한 저류암 혹은 비저류암의 물성을 반영하여 이산화탄소의 플룸이 확장되는 점을 관찰할 수 있다. 이를 3차원으로 관찰하기 위해서 이산화탄소 포화도 분포도를 그림 12에 도시했다. 단순한 지질 모델에서 이산화탄소 플룸은 규모, 형태, 이동 정도, 포화도가 비교적 단조롭고 균질한 것이 특징이다. 전반적으로 이산화탄소 플룸은 수평 이동 보다 밀도차에 의한 수직 이동이 일어나 실린더형으로 발달한다(그림 12). 비균질한 지질 모델에 대한 이산화탄소 거동 시뮬레이션과 다른 차이를 보이는 것은 염수층내에 분포하는 불투과성 퇴적상의 분포가 단조롭고 흐름을 방해하는 효과가 낮기 때문이다. 현재 시뮬레이션 조건하에서 단순한 지질 모델의 덮개층은 주입 이산화탄소를 충분히 막지 못하고 누출이 일어난다. 이는 실제 누출 가능성을 지시하기 보다 균질한 물성의 지질 모델에서 이산화탄소 주입에 따른 동적 물성의 변화가 단순화되어 나타난 결과로 해석된다.

Fig. 12. 
Simulated CO₂ saturation and migration for homogeneous facies model. Simulation time is 50 days. (a) CO₂ plume with saturation larger than 0.01. (b) Cross-sectional view of CO₂ saturation. (c) Three dimensional view of CO₂ saturation. (d) Breakthrough of CO₂ plume into seal formation.