Introduction

석유는 시추를 통해 채굴하는 전통 자원과 전통 자원과 다른 방식으로 채굴하는 비전통 자원으로 나뉜다. 비전통 자원의 대표적인 예로는 오일샌드(oil sand), 셰일오일(shale oil), 셰일가스(shale gas), 석탄층메탄(coalbed methane) 등이 있다. 비전통 자원은 석유, 가스 생산에서 오일은 9%, 가스는 27% 이상을 차지하고 있으며 비중은 커지고 있다⁽¹⁾.

오일샌드는 비전통 자원의 하나로 비투멘(bitumen)과 물, 점토나 모래 등이 섞인 혼합물 형태로서 지하에 존재한다. 지하에 매장되어 있는 오일샌드의 경우 SAGD(steam assisted gravity drainage) 공정을 통해 오일로 생산된다. 비투멘은 지하에 10000 cP(10 kg/mㆍs) 이상의 고점성 상태로 존재하여 유동성이 낮아 추출이 매우 어렵다. 따라서 고온의 증기(steam)을 주입하여 비투멘을 추출한다⁽²⁾.

고온의 증기를 주입하면 비투멘의 온도가 올라가며 점도는 낮아진다. 점도가 낮아진 비투멘은 광구의 하단으로 모이며 이러한 점도가 낮아진 비투멘을 유정(oil well)의 밑부분에 설치된 파이프를 통해 추출된다.

추출된 비투멘은 SAGD 초기공정에서 3상 유수분리기(FWKO, free water knock out)를 통해 물과 오일로 분리된다. 물과 오일간의 밀도 차이로 분리되며, 수평형(horizontal type)과 수직형(vertical type)이 있다. 수직형의 경우 모래와 같은 고체 입자의 분리에 이점이 있으며, 수평형의 경우 오일 액적(droplet)이 유수분리기 내부 유동과 반대되지 않아 분리가 쉽다는 장점이 있다.

유수분리기에서 물과 오일을 분리하는 이유로는 파이프라인에서의 부하 방지와 유수분리기 이후 공정에서의 설비 소형화를 통한 비용 절감, 가열 비용 절감 등이 있다.

유수분리기의 성능 향상을 위한 많은 연구들이 진행되어왔다. 그러나 캐나다의 경우 겨울에는 기온이 매우 낮은 극한지라는 점에서 극한지 환경 조건을 고려하였을 때 실험적 연구보다는 주로 수치해석적인 연구가 수행되어왔다. Lu 등⁽³⁾은 유수분리기의 내부 배플(baffle)이 분리효율에 미치는 영향에 대하여 연구하였다. Kang 등⁽⁴⁾은 배플의 배치 간격에 따른 유수분리기의 성능에 대하여 연구하였으며, Kwon 등⁽⁵⁾은 유수분리기의 방해판 개수에 따른 성능과 가동을 위한 최적의 watercut에 대하여 실험적으로 연구하였다. Kang 등⁽⁶⁾은 유수분리기의 입구에 베인 형상을 부여하여 분리효율에 대하여 연구하였으며, Jung 등⁽⁷⁾은 유수분리기의 위어 각도, 체류시간, 물과 오일의 비(water-to-oil ratio)를 변수로 설정하여 다양한 조건에서 유수분리기의 성능에 대하여 연구하였다.

본 연구에서는 위어 상단부 형상의 곡면을 조절하여 유수분리기의 성능에 대하여 수치해석적으로 연구하였으며, 또한 오일샌드를 채굴하는 광구마다 점도 차이가 있다는 점을 고려하여 점도에 따른 유수분리기의 성능에 대하여 수치해석적으로 연구를 진행하였다.

본 연구에서는 ANSYS사의 상용 코드 Fluent Ver 21.1을 사용하였다. 수치해석에서 다상 유동 모델로 VOF(volume of fluid)모델을 사용하였다. VOF 모델은 표면 추적 기법이며, 하나의 운동량 방정식을 각 상에 적용하며, 격자에서는 각 상의 부피만을 고려하여 계산을 한다. VOF모델의 지배 방정식은 다음과 같다.

운동량 방정식:

체적분율 보존 방정식:

해석에서 표면장력을 고려하는 경우 표면장력의 영향은 운동량 방정식에 추가적인 항을 고려하여 아래의 식으로 표현한다.

난류 모델로 realizable k-ε모델을 사용하였다. realizable k-ε모델은 소산(dissipation)에 대한 식과 eddy viscosity에 대한 새로운 식을 적용하여 기존의 k-ε모델보다 높은 정확도를 가진다는 장점이 있으며 다음 식으로 표현된다⁽¹⁰⁾.

해석은 과도상태(transient)로 진행하였으며, 계산시간단계(time step)는 0.001, 반복횟수(iteration)는 82400, 압력은 10 bar, 수렴조건은 10^-3로 설정하였다. 경계조건은 Table 3에 나타냈다.

유수분리기의 성능을 평가하는 인자로는 분리효율(η)로 설정하였는데, 오일 출구에서의 오일 유량(Q˙o)과 혼합물 유량(Q˙m=Q˙w+Q˙o)의 비로 정의하였다.

수치해석 연구의 정확성을 위해 격자의존시험(grid dependency test)을 수행하였고, 약 500000개로 설정하여 수치해석 실험을 진행하였다(Fig. 3 참조).

Fig. 3 
Grid dependency test

본 연구에서는 유수분리기의 위어 상단부 형상과 오일의 점도 변화에 따른 분리효율에 대해 수치해석적으로 분석하였다. 위어 상단에 곡면 형상을 부여하여, 오일 출구에서 오일과 혼합물의 비, 즉 분리효율을 통해 성능을 평가하였다. 오일의 점도 변화에 따른 유수분리기의 경우 분리효율과 오일 출구에서의 오일 유량을 통해 성능을 평가하였으며 결과를 요약하면 다음과 같다.

• 1) 유동 방향으로 곡면 형상을 준 경우 분리효율은 오일 출구에서 오일 유량과 혼합물 유량의 비로 계산하였으며, 곡률반지름 30 mm, 50 mm에서 각각 10.9% 증가한 91.89%, 11.59% 증가한 92.58%로 유동 방향으로 곡률 반지름이 증가할수록 분리효율이 증가하는 것을 확인하였다. 유동 반대 방향의 경우 반지름 30 mm에서 분리효율이 1.9% 감소한 79.09%를 기록하여 분리효율이 유동 반대 방향에서는 감소하는 것을 확인하였다. 따라서 유수분리기 위어 상단부에 곡률 반지름을 유동 방향으로 주어 유수분리기의 성능을 증가시킬 수 있을 것으로 판단된다.
• 2) 오일의 점도를 100 cP, 1000 cP, 10000 cP로 총 3가지 케이스로 나누어 수치해석연구를 진행하였다. 분리효율은 유출 유량과 유입 유량의 비로 계산하였으며, 100 cP는 54.37%, 1000 cP는 43.85%, 10000 cP는 12.86%로 측정되었다. 오일 유량의 경우 100 cP는 0.0093 kg/s, 1000 cP는 0.0075 kg/s, 10000 cP는 0.0022 kg/s로 측정되었다. 점도가 높을수록 유동성이 떨어져 위어 부근에서 오일이 흐르지 않으며 정체되어 생산되는 오일 유량이 적어지는 것을 확인하였다. 따라서 유수분리기의 생산 유량 측면에서 성능 향상을 위해 오일의 유동성 확보가 중요한 것으로 판단된다.