최근 정부의 수소경제 활성화 로드맵 발표(’19.1)와 함께 적극적인 정책 추진의 결과로 수소산업의 발전이 급속하게 이루어지고 있다. 수소에너지는 차세대 청정연료로서 수송·발전·건물·산업 등의 원료로 활용되고 있으며, 관련 인프라 확충과 기술개발로 인해 급격한 성장세를 보이고 있다. 특히, 1998년 국가 G7 사업 및 차세대 자동차 개발 사업을 기점으로 수소차 개발을 시작하여 교통 인프라 보급이 확대되고 있으며 2040년까지 수소차 290만 대, 수소충전소 1,200개소를 목표로 하고 있다. 이러한 수소차 보급 및 시장 활성화를 위해 수소충전소의 확대가 필요하나 이전까지 국내에서는 충전 인프라에 대한 신뢰성 있는 위험성평가 프로그램이 부재하였기 때문에 안전성 확인에 어려움이 있었다.

본 연구에서는 이를 위해 수소충전소 전용 위험성평가 프로그램인 Hy-KoRAM을 활용하여 제조식 수소충전소에 대한 정량적 위험성평가를 수행하였다.

수소충전소는 수소의 공급 방식에 따라 저장식과 제조식으로 나누어진다. 제조식 수소충전소는 수소의 원료가 되는 가스를 개질기와 같은 추출장치로 수소를 직접 생산하여 수소차에 공급하는 방식으로서 충전소 건설비용이 높아 국내에서는 주로 저장식 수소충전소 방식을 이용하고 있다. 하지만 수소충전차량의 수요가 점점 증가함에 따라 수소 공급 요구량이 증대되면서 여러 지역에서 제조식 수소충전소 설치사업을 확대 실시하고 있다. 따라서 본 연구를 통해 제조식 수소충전소에 대한 피해 범위 및 영향을 분석하여 안전성을 제고하고자 한다.

본 연구는 제조식 수소충전시스템의 피해범위와 영향을 분석하기 위하여 정량적 위험성평가 프로그램인 Hy-KoRAM을 활용하였다. Hy-KoRAM은 가스 누출에 따른 피해 범위 예측 및 사고 빈도를 고려하는 프로그램으로서 department of energy (DOE)에서 개발한 HyRAM^1,2)에 정성적·정량적 모듈을 보완하여 만든 위험성평가 프로그램이다. 현재, 국내에서는 Phast/Safeti, EFFECTS/RISKCURVES 등 주로 석유화학플랜트의 위험성을 평가하는 해외 범용 프로그램을 사용하고 있다. 하지만 해외 범용 프로그램은 국내환경 부합에 제한적이며 사용되는 고장률 데이터베이스에 따라 평가 결과의 편차가 크다.

따라서 본 연구에서는 최근 국내에서 국제공동연구를 통해 개발한 수소환경에 부합하는 위험성평가 프로그램인 Hy-KoRAM을 활용하여 제조식 수소충전소의 피해범위와 영향을 분석하고자 한다.

Hy-KoRAM은 수소충전소 사고 발생 시 피해반경을 기반으로 개인적 위험도 산출 및 주변 주민과 환경을 고려한 사회적 위험도(frequency-number curve) 산출이 가능하다. 또한, 수소 누출 및 제트화염 실증 시험을 기반으로 한 피해 영향 모델이 프로그램에 적용되었으며, 고장률데이터의 통계처리 방법을 이용해 수소설비에 적합하도록 failure frequency가 보정된다. Hy-KoRAM의 누출 제트 및 제트화염 모델 식은 HyRAM과의 산출방식과 동일하고 개인적위험도와 사회적위험도 산출 기능이 추가 되었으며 각각의 모델식은 아래와 같다^1,3).

Jet release model equations

υ : velocity, ρ : density, cl : centerline, amb : ambient
B : characteristic half-width
λ : ratio of density spreading relative to velocity
r : perpendicular to the stream-wise direction
Y : mass fraction of hydrogen

Jet flame model equation

Srad : total emitted radiative power
C* : normalized radiative heat flux
X_rad : radiant fraction
m˙fuel: mass flow rate of fuel
∆H_c : heat of combustion (118.83 MJ/kg)
τ_f : flame residence time
f_s : mass fraction of hydrogen in a stoichiometric mixture
ρ_f , W_f and L_vis : flame density, width and visible length
ρ_j , d_j and u_j : density, diameter and velocity of jet at exit

Individual risk

IR_x,y : Total individual risk from death at geographic location “x,y”
IR_x,y,i : Individual risk from death at geographic location “x,y” of accident case “i”
n : Total number of accident result cases considered in the analysis
f_i : Frequency of accident case “i” in frequency analysis
P_f,i : Probability of accident case “i” causing death at location “x,y”

사회적 위험도를 산출하는 F-N 곡선을 계산하기 위한 과정은 개인적 위험도 산출 방식과 동일하며, 따라서 각 사고 결과 사례의 영향을 받는 사람의 수를 예측하기 위해서는 해당 정보를 인구 자료와 결합하여야 한다.

Societal risk

N_i : Number of deaths caused in accident case “i”
P_x,y : Number of people at location “x,y”

본 연구는 제조식 수소충전소에 대해 국내에서 개발한 위험성 평가프로그램인 Hy-KoRAM을 이용하여 정량적 위험성평가를 수행하였다. 각 시설의 위험요인을 분석하기 위하여 국제표준(IEC 60079-10-1)에서 제시하고 있는 2차 누출 조건을 적용하였으며, 수소충전소 사고에 따른 피해 범위 및 영향을 확인하였다. 연구결과를 요약하면 다음과 같다.

1) 피해영향분석에서 2차 압축기의 피해영향이 가장 큰 것으로 나타났으며 압축기와 압축가스 설비가 방호벽 내의 공동 공간에 배치되어 있는 상태이기 때문에 한 설비에서 사고 발생 시 다른 설비로 사고가 전이 될 수 있어 격벽 등의 방호조치로 안전성을 강화하는 것이 권장된다.

2) 본 연구의 대상인 제조식 수소충전소의 개인적 위험도 분석 결과 작업자에 대한 허용 불가 영역(>1×10^-3) 및 작업장 밖에서의 일반 대중에 대한 허용 불가 영역(>1×10^-4)은 나타나지 않았으며, 충전소 주변의 인구현황 분포를 적용한 사회적 위험도 평가 결과 ALARP 영역 및 완전 허용 가능 구역 내에 분포하는 것으로 확인되었다.