API 5L X65 배관의 신뢰도 평가를 위한 파열압력 분포 추정에 관한 연구
A Study on the Estimating Burst Pressure Distributions for Reliability Assessment of API 5L X65 Pipes
Article information
Trans Abstract
Purpose
The purpose of this paper is to present a probability distribution of the burst pressure of API 5L X65 pipes for the reliability assessment of corroded gas pipelines.
Methods
Corrosion is a major cause of weakening the residual strength of the pipe. The mean residual strength on the corrosion defect can be obtained using the burst pressure code. However, in order to obtain the pipe reliability, a probability distribution of the burst pressure should be provided. This study is concerned with estimating the burst pressure distribution using Monte Carlo simulation. A response surface method is employed to represent the distribution parameter as a model of the corrosion defect size.
Results
The experimental results suggest that the normal or Weibull distribution should be suitable as the probability distribution of the burst pressure. In particular, it was shown that the probability distribution parameters can be well predicted by using the depth and length of the corrosion defect.
Conclusion
Given a corrosion defect on the pipe, its corresponding burst pressure distribution can be provided at instant. Subsequently, a reliability assessment of the pipe is conducted as well.
1. 서 론
배관 표면에 발생하는 부식은 잔류강도를 약화시켜 배관을 고장에 이르게 하는 주요 원인 중 하나이다. 배관의 고장은 잔류강도 즉 배관의 파열압력이 운영압력보다 작을 때 발생한다고 간주한다. 따라서 배관의 신뢰도는 파열압력이 운영압력보다 클 확률로 정의된다. 신뢰도를 구하기 위해서는 파열압력과 운영압력에 대한 확률분포가 필요하다. 보통 운영압력의 확률분포는 공급관리소에서 측정된 데이터를 이용하여 쉽게 얻을 수 있지만, 반면에 파열압력에 대한 확률분포를 구하는 것은 매우 복잡하다. 본 논문은 파열압력의 확률분포를 추정하는 방안을 다루고자 한다.
배관 매설기간이 증가하면서 부식은 배관 고장에 중대한 원인으로 부각됨에 따라 배관의 신뢰도 및 안전성 확보에 대한 관심이 점차 커지고 있다. 이미 오래전부터 매설배관을 운영해 온 서유럽에서는, 배관 고장의 약 27%가 부식에 기인한 것으로 요약되고 있다(Dundulis et al. 2016). 부식으로 인한 고장은 특히 2009년 이후 크게 증가하고 있는데 이는 배관 사용연수가 진행됨에 따라 부식의 영향이 점차 커지고 있음을 의미한다. 또한 세계적으로도 배관부식에 따른 직접비용은 연간 약 1.4조 달러 수준으로서 이는 전세계 GNP의 약 2%에 해당된다(Ossai 2015). 부식결함이 발생한 지점에서의 배관 신뢰도 평가기술은 장기 사용 배관의 건전성 관리 및 안전성 확보에 토대가 된다.
부식손상 배관의 파열압력을 계산하는 방법은 많은 연구에서 찾아볼 수 있다. 파열압력 계산을 위한 대표적인 코드로는 PCORRC, ASME Modified B31G, DNV99, Shell92 등을 들 수 있다(Caleyo et al. 2002). 하지만 이들 계산 코드는 파열압력의 평균을 얻기 위한 것이고 파열압력의 산포에 대한 연구는 상대적으로 미흡한 실정이다. 배관의 신뢰도 평가와 수명예측에는 파열압력의 산포에 대한 정보가 필요하고 특히 확률을 계산하기 위해서는 파열압력의 확률분포를 구해야 한다. 본 논문은 몬테카를로 시뮬레이션을 이용하여 파열압력의 확률분포를 추정하는 절차를 제안하고자 한다. 또한 본 연구는 부식결함의 크기가 주어질 때 확률분포의 파라미터를 추정할 수 있도록 반응표면 실험계획을 이용한 모델링 방안을 제시한다. 국내 가스배관으로 널리 사용 중인 API 5L X65 등급의 배관을 대상으로 제안된 방법의 적용과정을 예시한다. 본 연구의 결과는 부식이 검출된 배관의 신뢰도 평가 및 유지보수 의사결정에 횔용될 수 있을 것으로 기대된다. 본 논문의 구성은 다음과 같다. 2장은 PCORRC 코드를 이용한 파열압력 계산과 부하-강도 모델을 이용한 신뢰도 평가절차에 대해 기술한다. 3장에서는 파열압력 분포를 추정하기 위해 몬테카를로 시뮬레이션을 이용하는 방안을 설명하고 반응표면분석법을 통해서 파라미터 추정 모델을 개발하는 과정을 제시한다. 마지막으로 4장은 연구결과의 요약 및 결론을 수록하고 향후 연구방향에 대해 전망한다.
2. 배관 파열압력 및 신뢰도 평가
2.1 부식결함과 배관의 파열압력
배관의 부식은 장기사용 배관의 위험요인이다. 배관 표면에 발생한 부식결함을 식별하고 시간에 따른 부식성장을 모니터링하는 것은 장기사용 배관의 건전성 관리에 중요한 역할을 차지한다(Kim et al. 2017). 사용 중인 배관의 상태를 점검하기 위해 라인내 검사(In-line Inspection)를 정기적으로 실시한다. 이 때 배관 내부로 투입되는 측정 장비를 피그(Pig)라고 부른다. 라인내 검사는 크게 Magnetic Flux Leakage(MFL) 방식과 Ultrasonic 방식으로 구분되는데 다음 그림은 MFL 검사로 얻은 3차원 결함정보를 배관두께에 대한 2차원 변수로 변환하는 과정을 개략적으로 나타낸 것이다. 이렇게 얻어진 결과를 토대로 부식결함의 발생 위치, 형상, 크기 등을 기록한 파일을 피깅(Pigging) 데이터셋이라고 부른다. 일반적으로 MFL 방식에 따른 피깅 데이터는 10% 내외의 측정 및 계산 오차를 포함하는 것으로 알려져 있다(Cosham et al. 2007).
그림 1에 나타낸 바와 같이, 부식결함의 크기는 깊이 d와 길이 L로 기록된다. w는 배관의 두께를 의미한다. 캐나다 표준협회가 제정한 석유 및 가스 배관에 관한 표준규격 CSA Z662-2007에 따르면, 부식과 관련된 배관 고장은 Small Leak, Large Leak, Rupture 등 3가지로 구분된다. 이는 아래 그림 2에서 보는 바와 같이 부식결함의 깊이와 길이를 기준으로 한 것이다.
Small Leak는 길이는 매우 짧지만 깊이가 큰 핀홀 형태의 부식결함으로 인해 가스가 누설되는 것을 말한다. Large Leak는 국부적인 부식결함으로 인해 배관의 잔류강도가 약화되면서 가스가 누설되는 것이고, Rupture는 배관 전역에서 부식이 광범위하게 진행되어 심각한 두께 및 강도 손실과 함께 파단이 일어나는 것을 뜻한다. 그림 안에 g1, g2, g3라 표기한 것은 이들 3가지 고장에 대응하는 한계상태함수이다. 그림을 보면, g2가 대부분의 영역을 포괄하고 있음을 알 수 있는 데 바로 배관의 파열(Burst)이 여기에 해당된다.
파열압력이란 파열이 발생하는 압력을 말하는 데 배관 내부의 압력이 이를 초과할 경우 배관의 파열이 발생한다는 의미이다. 파열압력은 배관의 강도, 직경, 두께 등의 함수이고 특히 부식 깊이와 길이에 의해 큰 영향을 받는다. 이에 따른 파열압력 계산식으로 널리 이용되는 코드는 PCORRC, ASME Modified B31G, DNV-99 등이다. 각각의 특징이나 차이점에 대해서는 Caleyo et al.(2002)의 연구를 참고하기 바란다. 본 논문은 선행연구인 Kim et al. (2019)에서와 같이 국내에서 널리 채택되고 있는 PCORRC를 사용한다. 이때 파열압력 Pb의 계산식은 다음과 같다.
단 U, D, w는 배관의 인장강도, 직경, 두께를 각각 나타낸다. 그리고 d와 L은 부식결함의 깊이와 길이를 의미한다. 부식 깊이와 길이는 시간에 따라 커지게 되므로 파열압력은 시간에 따라 감소하는 형태가 된다. 특히 d가 w로 접근하면 파열압력은 0으로 수렴하게 됨을 알 수 있다. 반면 d가 충분히 작을 때에는 L은 파열압력에 그리 큰 영향을 미치지 않는다.
2.2 파열압력을 이용한 배관 신뢰도 평가
배관의 신뢰도 평가를 위해서는 우선 배관의 고장을 정의하는 한계상태함수를 결정해야 한다. 본 논문에서는 CSA Z662-2007에서 제시한 g2를 채택하고자 한다. 이는 파열압력 Pb과 운영압력 Pop의 차이로서 g2 = Pb - Pop와 같이 쓸 수 있다. 전술한 바와 같이, 운영압력이 파열압력을 초과하면 배관의 고장이 발생한다는 의미이다. 파열압력과 운영압력을 확률변수라고 하면 배관의 신뢰도 R은 파열압력이 운영압력보다 클 확률 또는 한계상태함수가 0보다 클 확률로 정의되어지며 아래 식으로 나타낼 수 있다.
따라서 신뢰도를 구하기 위해서는 Pb와 Pop의 결합 확률분포가 필요하다. 보통은 Pb와 Pop는 서로 독립으로 간주될 수 있으므로 개별 확률분포가 주어지면 강도-부하 모형을 이용하여 어렵지 않게 신뢰도를 계산할 수 있다. 예를 들어 확률변수 Pb와 Pop가 서로 독립이며 정규분포를 따른다고 가정하자. 그리고 Pb의 평균과 표준편차를 각각 μb와 σb로 놓고 Pop의 평균과 표준편차를 각각 σop와 σop라 놓으면 식 (2)의 신뢰도는 다음과 같이 나타낼 수 있다.
단 Φ(·)는 표준정규분포의 분포함수를 나타내고 특히
물론 이러한 신뢰도 평가는 파열압력과 운영압력이 서로 독립인 정규분포를 따를 때 가능하다. 하지만 정규분포를 따른다는 보장이 없고 심지어 파열압력에 대한 확률분포는 알 수 없는 경우가 많다. 사실, 운영압력에 대한 확률분포는 배관운영회사에서 측정을 통해 용이하게 확보할 수 있지만 현재 사용하고 있는 배관의 파열압력은 측정 자체가 불가능하거나 식 (1)로도 볼 수 있듯이 많은 변수의 영향과 불확실성 하에 있기 때문에 확률분포를 수리적으로 유도하는 것은 매우 복잡하다. 배관 신뢰도를 구하고자 할 때 등장하는 이러한 문제점을 해결하기 위해 근사해법 또는 Monte Carlo Simulation(MCS)이 적용되어 왔다.
근사해법으로는 First Order Second Moment(FOSM), First Order Reliability Method(FORM), Second Order Reliability Method(SORM) 등을 들 수 있다. 이들은 적용하기 간편하지만 정규분포를 가정한다는 점과 식 (1)의 파라미터 설정에 민감하다는 점이 한계이다. 특히 이들은 분산을 추정할 뿐이지 분포 자체를 추정하는 방법은 아니라는 점도 상기해야 한다. MCS를 이용하면 식 (1)을 이용하여 많은 표본을 발생시켜서 경험적 분포를 구하거나 기존 확률분포에 적합시키는 것이 가능하다. 예를 들어, Caleyo et al. (2009)는 MCS를 이용하여 부식 깊이와 성장속도에 대한 확률분포를 추정하였다. Hasan et al. (2012)는 MCS를 이용하여 파열압력의 분산을 추정하고 신뢰도 지수를 계산하였다. Qian et al. (2013)은 한계상태함수를 평가하고 고장확률을 계산하는 데 MCS를 이용하였다. 이처럼 MCS는 배관 신뢰도 문제를 다룰 때 다각적으로 적용되고 있지만 파열압력의 확률분포를 직접적으로 다룬 연구는 상대적으로 미흡한 실정이다.
본 논문은 배관의 신뢰도 평가에 요구되는 파열압력의 분포추정을 위해 MCS를 이용하는 방안을 제안한다. 특히 부식결함의 깊이와 길이 수준에 따라 다양한 확률분포가 나타날 수 있기 때문에 반응표면 실험계획법을 이용하여 분포 파라미터 모델을 구축하고 활용하는 과정을 API 5L X65 배관을 대상으로 예시한다. 아래 그림 3은 본 연구의 내용을 토대로 API 5L X65 배관의 신뢰도를 평가하는 절차를 보여주고 있다.
3. 파열압력의 분포 추정과 파라미터 모델링
3.1 파열압력 계산을 위한 설정과 Monte Carlo simulation
API 5L X65 등급 배관은 국내 천연가스 공급을 담당하고 있는 주요 배관 중 하나이며 직경은 508mm 또는 762mm, 두께는 10.3mm에서 17.5mm까지 다양한 제품군으로 이루어져 있다. 본 연구는 17.5mm 배관을 대상으로 파열압력 분포의 추정 문제를 다룬다. 그림 3에서 나타낸 바와 같이, 파열압력을 계산하기 위해서는 인장강도, 직경, 두께, 부식 깊이, 부식 길이 등을 입력해야 한다. X65 배관의 경우 각 변수들의 공칭 값은 표 1과 같이 주어진다.
각 변수의 확률분포는 CSA Z662-2007이 권고한 내용을 따르기로 한다. 표 1에서 볼 수 있듯이, 기초 변수의 확률분포는 서로 독립인 정규분포를 가정한다. 다만 부식결함의 크기에 대해서는 정규분포 외에도 와이블분포나 대수정규분포의 적용이 가능하다. MCS 수행을 위해서는 확률분포의 평균과 표준편차가 필요하다. 평균은 표 1의 공칭 값을 그대로 사용하지만 표준편차는 관련 규격이나 사내 운영데이터를 이용해서 결정해야 한다. 배관공학 분야에서는 표준편차 대신 변동계수(COV)를 많이 이용하는 데, CSA Z662-2007이 제시한 변동계수 설정기준에 따르면 인장강도는 4%, 배관직경은 0.06%, 배관두께는 1% 정도가 적절할 것으로 판단된다. 배관 시스템의 신뢰도 기법에 관한 국제규격인 ISO 16708-2006(E)에서도 인장강도 1.5-6.0%, 배관직경 0.1% 미만, 배관두께 0.5%-2.0%로 변동계수를 다룰 것을 권고하고 있으므로 대체로 비슷한 수준이다. 다만 부식결함의 경우는 매설 지역과 기간에 따라 변동이 크기 때문에 현장에서 수집된 데이터에 근거하여 변동계수를 설정하는 것이 타당할 것이다. 예를 들어, Qian et al.(2013)은 부식결함의 깊이와 길이에 대해 변동계수를 각각 10%로 채택하였다. 이보다는 큰 값으로 Zhang et al.(2019)의 연구에서는 깊이를 약 40%, 길이를 약 20%로 변동계수를 적용하였다. 본 연구에서는 10km 이상 구간에서 수집된 6가지 피깅 데이터셋을 검토한 결과, 깊이와 깊이의 COV가 각각 15-48%, 14-51% 수준으로 얻어졌다. 이를 근거로 할 때 부식결함의 COV는 대략 10-50% 범위에서 선택하되 해당 구간의 매설기간이나 환경적 특성을 고려하는 것이 바람직할 것으로 사료된다. 참고로, 그림 4는 어떤 3개의 구간에서 검출된 부식결함의 깊이를 상자 그림으로 간단하게 나타낸 것이다. 대상 구간에 따라서 변동에 차이가 나타날 수 있음을 예시하고 있다.
그림 5는 어떤 부식결함에 대해서 깊이와 길이의 COV를 각각 50%로 하였을 때 MCS를 통해 얻은 파열압력의 분포를 보여준다. 런 수가 늘어남에 따라 분포의 모습이 안정화되고 있다. 대체로 정규분포에 가까운 모습이고 런 수가 200 정도이면 큰 변화를 나타내지 않으며 분포를 추정하는 데 충분할 것으로 판단된다. 이러한 관측을 근거로 본 연구에서는 MCS 런 수로 200회로 설정한다.
그림 5의 부식결함은 Depth to Thickness Ratio와 Length to Diameter Ratio를 각각 30%, 15%로 적용한 것이다. 이후로는 이를 (d/t)=0.3, (L/D)=0.15로 표기하기로 한다. 여기서 이 조건을 선택한 이유는 그림 2에 나타낸 고장 발생 주요 영역의 중앙 지점에 해당되기 때문이다. 전술한 바와 같이 (d/t)와 (L/D)는 매설 지역과 기간에 따라 다양하게 변동하기 때문에 여러 조합에서 파열압력의 변화를 관찰할 필요가 있다. 본 연구는 반응표면 실험계획법을 이용하여 (d/t)와 (L/D)의 대표 조합을 구하고 파열압력 분포의 파라미터를 모델링하고자 한다.
3.2 파라미터 모델 구축을 위한 반응표면 실험계획법
반응표면 실험계획법은 공정이나 제품의 최적 설계조합을 찾고자 품질공학 분야에서 널리 이용되는 기법으로, Kim et al. (2016)의 리뷰논문을 통해 품질 분야의 활용 현황을 찾아볼 수 있다. 반응표면 실험계획법은 Central Composite Design과 Box-Behnken Design으로 나뉘는 데 본 실험에서는 설계인자가 X1 = (d/t)와 X2 = (L/D), 2개이므로 전자를 채택한다. 기존의 피깅 데이터셋을 검토했을 때 X1과 X2의 실험영역은 각각 0.1에서 0.5, 0.05에서 0.25로 정하면 적절할 것으로 판단된다. 이 경우, 중심점은 (0.30, 0.15)가 되고 축점은 (0.5828, 0), (0.0172, 0), (0, 0.2914), (0, 0.0086)이 된다. 중심점에서의 반복을 5회로 할 때 전체 실험점은 아래의 표 2와 같다.
표에는 X1과 X2의 실험조건과 함께 이들을 표준화한 U1과 U2의 실험조건도 포함되어 있다. 즉 U1과 U2는 U1 = (X1-0.30)/0.20과 U2 = (X2-0.15)/0.10와 같이 구한 것이다. 또한 표 2에는 각 실험조건에서 MCS를 이용하여 생성된 200개의 파열압력 표본을 이용해서 얻어진 표본평균 Y1과 표본표준편차 Y2도 수록하고 있다. 앞서 언급한 바와 같이 파열압력은 X1과 X2의 COV에 크게 의존하기 때문에 이를 10%, 30%, 50%로 변화시켜가면서 표본을 얻도록 하였다. 예를 들어, 부식 깊이와 길이의 COV를 각각 30%로 설정한 후 가장 최악으로 예상되는 실험조건 6번과 중심점에 대응되는 실험조건 9번에 대한 실험결과를 그래프로 나타내면 아래 그림 6과 같다. 그림에는 정규분포, 대수정규분포, 감마분포, 와이블분포를 이용한 적합 결과를 동시에 보여주고 있다. 두 실험결과 모두 정규분포로 잘 묘사되는 것으로 판단된다. 한편 와이블분포 적합은 분포의 치우침을 보다 적극적으로 수용하는 모습을 보인다.
분석 결과, Y1과 Y2에 대한 반응표면 회귀식은 다음과 같이 구할 수 있다.
어떤 부식결함에 대해서 그 깊이와 크기 즉 (X1, X2)이 주어지면 상기의 회귀식을 이용하여 해당 파열압력의 분포를 추정할 수 있게 된다. 여기에 운영압력의 분포가 정해지면 배관의 신뢰도를 계산할 수 있다. 그리고 부식성장 예측치가 주어질 때에는 식 (4)와 (5)를 이용하여 시간에 따른 파열압력의 분포를 구할 수 있으므로 배관의 잔여수명 예측이 가능해진다.
3.3 적용절차 예시
제안된 방법의 적용절차를 예시하기 위해 동일한 테스트 구간에서 하나의 부식결함이 검출되었다고 가정하자. 이 결함의 깊이와 길이가 (X1, X2) = (0.25, 0.30)라고 할 때 먼저 평균과 표준편차를 추정한다. 이를 위해 식 (4)와 (5)를 이용하면 (
이 결과를 고장확률로 나타내면
단, Pop와 D는 내부 운영압력(bar)과 배관 직경(m)을 각각 의미한다. 식 (6)에 Pop = 6.86Mpa = 68.6bar와 D = 0.762m를 대입하면 안전등급 4일 때의 허용한계는 plim = 1.647×10-7과 같이 얻어진다. 따라서 상기의 결함은 매우 안전한 상태에 있다는 결론이 가능하다.
신뢰도 계산에 이어 잔여수명 평가를 수행하려면 우선 부식성장 예측치가 필요하다. 부식성장을 다루기 위해 선형 모델이나 멱함수 모델이 널리 이용된다. 선형 모델이 부식크기를 과대추정하는 경향이 있기 때문에 멱함수 모델이 현실적으로 타당하지만 적절한 파라미터를 추정하는 것이 어렵다는 단점이다(Caleyo et al. 2009, Zhang et al. 2019). 여기서는 예시를 위해 부식이 시간에 선형적으로 매년 0.5%씩 증가하는 것으로 가정하고 잔여수명을 평가하기로 한다.
부식성장 예측치가 해마다 주어지면 식 (4)와 (5)를 이용하여 (Y1, Y2)을 추정, 신뢰도를 구할 수 있다. 신뢰도 계산값이 1에 가깝기 때문에 고장확률에 로그를 취해서 나타내는 것이 편리하다. 아래의 그림 7은 시간 경과에 따른 고장확률이 점차 증가하고 있음을 보여주고 있다. 고장확률은 약 19년 후에 ISO 허용한계를 초과할 것으로 예상되므로 배관의 잔여수명은 19년이라고 예측할 수 있다. 이 시점에서 예상되는 부식크기는 (X1, X2)= (0.345, 0.395)이고 파열압력의 평균과 표준편차는 (
4. 결론
본 논문은 부식결함이 존재할 때 배관의 신뢰도 평가를 위한 파열압력 분포의 추정방법을 Monte Carlo 시뮬레이션과 반응표면 실험계획법을 통해 제시하였다. 배관의 매설 시기와 구간에 따라 다양한 부식결함이 나타나게 되고 피깅 데이터의 규모도 점차 증가하고 있기 때문에 효율적인 신뢰도 평가절차가 절실하게 요청되며 이는 배관의 유지보수 의사결정에도 중요한 역할을 차지한다. 또한 신뢰도 평가를 위해 파열압력의 분포를 직접 다룬 연구가 상대적으로 미흡하고 근사해법은 분포가 아니라 분산을 구해서 신뢰도를 평가하는 방식이라는 점에서 본 연구의 의의를 찾을 수 있다.
배관의 파열압력은 인장강도, 직경, 두께, 부식결함 등 여러 변수의 영향을 받기 때문에 확률분포를 이론적으로 구하는 것은 거의 불가능하다. 이에 본 논문은 Monte Carlo 시뮬레이션을 이용하였으며 피깅 데이터로부터 부식 깊이와 길이가 주어지면 그에 상응하는 파열압력의 분포 파라미터를 추정할 수 있도록 반응표면 회귀식을 개발하였다. 본 연구의 파열압력 분포추정 모델은 부식성장 예측치가 주어졌을 때 잔여수명 평가를 효율적으로 수행할 수 있으므로 배관 유지보수 실무에서도 유용할 것으로 기대된다.
본 연구는 API 5L X65 17.5mm 배관을 대상으로 하였으며 PCORRC 파열압력 계산식을 이용하여 이루어졌다. 향후 X42 등 현장에서 많이 쓰이는 타 배관에 대해서도 적용이 확대될 수 있도록 연구내용이 보완되어야 하겠고, 잔여수명 예측에 필요한 고장확률 허용한계를 거리당 연간 평균값으로 환산하는 기준도 마련되어야 한다. 끝으로, 제안된 방법은 부식결함이 하나로 주어졌을 때 적용가능하다는 점에서 향후에는 다수의 부식결함을 다룰 수 있는 신뢰도 및 잔여수명 평가방법 연구가 필요하다.