콘크리트 내부로의 염화물 침투는 어떤 확산 거동을 보이며 어떤 법칙을 따르나요?

콘크리트 내부로의 염화물 침투는 시간에 따라 농도 구배가 변하는 비정상 상태(Non-steady state) 확산 거동을 보이며, 1차원 반무한체 조건에서 Fick의 제2법칙을 따릅니다.

염화물 확산 계수 역해석과 Fick의 제2법칙을 활용한 콘크리트 잔존 수명 예측 엔지니어링

2026년 9월 10일 읽는 시간 약 11분

해안가 교량, 항만 케이슨, 해양 플랜트 구조물의 수명을 갉아먹는 가장 지배적인 열화 요인은 비래염분과 해수에 포함된 염화물 이온(Cl⁻)입니다. 콘크리트 타설 초기에는 피복두께가 철근을 보호하지만, 세월이 흐르며 침투한 염화물이 철근 위치에서 임계 부식 한계 농도(Critical Chloride Threshold)에 도달하면 부동태 피막이 파괴되고 급격한 팽창 균열이 발생합니다.

설계 단계에서는 실내 촉진 시험(NT BUILD 492 등)으로 얻은 이론적 확산 계수를 적용하지만, 수십 년간 공용 중인 구조물의 실제 잔존 수명을 정밀 예측하려면 실측 염화물 프로파일 데이터를 바탕으로 확산 계수를 ‘역해석(Back-analysis)’해야 합니다.

본 글에서는 Fick의 제2법칙 지배 방정식, 오차함수를 활용한 비선형 커브 피팅 역해석 절차, 그리고 시간 의존성(Aging Factor)을 반영한 잔존 수명 평가식을 정리합니다.

1. Fick의 제2법칙과 오차함수(erf) 기반 1차원 확산 모델

콘크리트 내부로의 염화물 침투는 시간에 따라 농도 구배가 변하는 비정상 상태(Non-steady state) 확산 거동을 보입니다. 1차원 반무한체(Semi-infinite medium) 조건에서 지배 방정식은 Fick의 제2법칙을 따릅니다.

[Fick의 제2법칙 편미분 방정식]

∂C / ∂t = D · (∂²C / ∂x²)

C : 깊이 x와 시간 t에서의 염화물 농도 (kg/m³ 또는 시멘트 질량 대비 %)

D : 겉보기 염화물 확산 계수 (Apparent Chloride Diffusion Coefficient, m²/s 또는 cm²/year)

x : 콘크리트 표면으로부터의 침투 깊이 (m 또는 mm)

t : 염화물 노출 경과 시간 (초 또는 년)

다음과 같은 현장 경계 조건을 적용합니다.

초기 조건 (t = 0): C(x, 0) = C₀ (초기 콘크리트에 포함되어 있던 잔류 염화물 농도)
경계 조건 (x = 0): C(0, t) = Cs (표면 염화물 농도, 대기 비래염분이나 해수 접촉에 의해 일정하게 유지된다고 가정)

위 조건을 만족하는 해석적 해(Analytical Solution)는 가우스 오차함수(Error Function, erf)를 통해 다음과 같이 유도됩니다.

[Fick의 제2법칙 닫힌 형태의 해]

C(x, t) = C₀ + (Cs – C₀) · [ 1 – erf( x / (2 · √(D · t)) ) ]

erf(z) : 가우스 오차함수 (0에서 z까지 정규분포 적분값)

[ 1 – erf(z) ] = erfc(z) : 여오차함수 (Complementary Error Function)

따라서 식을 간결하게 쓰면 다음과 같습니다.

C(x, t) = C₀ + (Cs – C₀) · erfc( x / (2 · √(D · t)) )

2. 현장 코어 시편을 활용한 확산 계수(D) 역해석(Back-analysis) 절차

공용 중인 구조물에서 실제 확산 계수를 구하기 위해 실측 프로파일 데이터를 비선형 회귀분석으로 역추적합니다.

[염화물 확산 계수 역해석 4단계 흐름]

1. 현장 코어 채취 (직경 100mm 코어 드릴링)
                ↓
2. 건식 다이아몬드 밀링 (표면부터 2~5mm 두께 단위로 분말 샘플링)
                ↓
3. 전위차 적정법 (각 깊이별 전염화물/수용성 염화물 농도 Ci 정량 분석)
                ↓
4. 비선형 최소자승법(NLS) 피팅 (Cs 및 유효 확산 계수 Dapp 동시 역산)

비선형 회귀분석(Non-linear Least Squares) 원리

실측 깊이(x_i)에서 시험으로 얻은 염화물 농도(C_meas(x_i))와 모델 공식에 의해 계산된 농도(C_calc(x_i)) 사이의 잔차 제곱합(RSS, Residual Sum of Squares)을 목적함수로 둡니다.

[목적함수 최소화식]

RSS = Σ [ C_meas(x_i) – C_calc(x_i, t_exposure, Cs, D) ]² → 최소화(Minimize)

엑셀의 ‘해답 찾기(Solver)’ 또는 파이썬(SciPy curve_fit), MATLAB을 이용해 RSS를 최소화하는 최적의 표면 염화물 농도(Cs)와 겉보기 확산 계수(D) 두 변수를 동시에 역산해 냅니다.

[염화물 농도 프로파일 피팅 곡선 개념도]

  농도 C(x)
    ▲
 Cs ● (표면 농도)
    │ ╲
    │   ╲  ● 실측 데이터 포인트 (실제 코어 분석값)
    │     ╲
    │─────── 모델 피팅 곡선: C(x) = Cs · erfc( x / 2√(Dt) )
    │         ╲
    │           ╲
 C₀ ┼─────────────●───────●───────●─── (초기 농도 C₀)
    └─────────────────────────────────▶ 깊이 x (mm)
                       ▲
                    피복두께 위치

3. 확산 계수의 시간 의존성: 에이징 팩터(Aging Factor, m)의 적용

Fick의 단순 제2법칙은 확산 계수(D)가 시간에 따라 변하지 않는 상수라고 가정하지만, 실제 현장 콘크리트는 시간이 지날수록 시멘트와 혼화재(플라이애시, 고로슬래그)의 지속적인 수화 반응으로 공극 구조가 치밀해집니다.

따라서 확산 계수(D)는 시간이 흐를수록 급격히 감소하며, 이를 반영하지 않으면 잔존 수명이 극단적으로 과소평가(조기 파괴로 잘못 예측)됩니다. Life-365 모델 및 국내 콘크리트학회 규준에서는 다음과 같은 시간 의존성 수식을 채택합니다.

[시간 의존적 확산 계수 모델식]

D(t) = D_ref · (t_ref / t)^m

D(t) : 재령 t 시점에서의 염화물 확산 계수

D_ref : 기준 재령(통상 t_ref = 28일)에서 측정한 기준 확산 계수

t : 실제 공용 경과 시간

m : 연령 계수 (Aging Factor, 바인더 종류에 따른 감쇄 지수)

바인더 종류별 연령 계수(m) 적용 기준치

보통 포틀랜드 시멘트 (OPC 단독 배합): m = 0.20 ~ 0.25 (공극 치밀화 효과가 미미)
플라이애시 (FA 20~30% 치환 배합): m = 0.50 ~ 0.60 (장기 포졸란 반응으로 확산 계수가 지속 감소)
고로슬래그 미분말 (GGBS 40~60% 치환 배합): m = 0.60 ~ 0.70 (치밀한 C-S-H 겔 형성 및 염화물 고정화 능력 우수)

[실무 계산 주의: 감쇄 한계 시간(t_real)]

수화 반응에 의한 공극 치밀화는 무한정 지속되지 않습니다. 보통 공용 25~30년 시점에 도달하면 조직 발달이 수렴하므로, t > 30년 이후에는 D(t)를 30년 시점의 고정값으로 유지하는 것이 실무 표준입니다.

4. 철근 부식 개시 시점 및 잔존 수명(Residual Service Life) 정량 산정

콘크리트 구조물의 내구적 한계 상태는 일반적으로 철근 위치(깊이 = 피복두께 cover)에서 염화물 농도가 부식 발생 임계 농도(C_th)에 도달하는 순간으로 정의합니다.

부식 개시 수명(T_initiation) 유도 공식

단순 확산 모델을 기준으로 잔존 부식 개시 수명(t_corr)을 역산하면 다음과 같습니다.

[부식 개시 수명 산정식]

T_init = (cover)² / [ 4 · D_app · ( erfc⁻¹( (C_th – C₀) / (Cs – C₀) ) )² ]

T_init : 타설 후 부식 개시까지 걸리는 총 소요 연수 (년)

cover : 설계 철근 피복두께 (mm 또는 m)

C_th : 철근 부식 임계 염화물 농도 (국내 기준: 보통 시멘트 중량 대비 0.40 wt%, 또는 콘크리트 단위 용적당 1.2 kg/m³)

erfc⁻¹ : 역여오차함수 (Inverse Complementary Error Function)

최종적인 현시점 기준 잔존 수명(Remaining Life)은 다음과 같이 결정됩니다.

잔존 수명(Remaining Life) = T_init – 공용 경과 년수(t_current)

만약 잔존 수명이 음수(-)로 나온다면, 해당 구조물은 이미 철근 위치에서 임계 부식 농도를 초과하여 활성 녹 발생 및 단면 손실 단계에 진입했음을 의미합니다.

5. 실무 잔존 수명 평가 계산 워크스루 예시

구조물 현황: 공용 10년 경과(t = 10년)된 서해안 교각, 실측 평균 피복두께 = 50 mm
코어 역해석 결과: 표면 농도(Cs) = 3.5 kg/m³, 역해석 평균 확산 계수(D_app) = 1.2 × 10⁻¹² m²/s (약 37.8 mm²/year)
내구성 기준치: 초기 염화물(C₀) = 0.0 kg/m³, 부식 임계 농도(C_th) = 1.2 kg/m³

Step 1. 여오차함수 비 계산

C_th / Cs = 1.2 / 3.5 ≈ 0.343

erfc(z) = 0.343을 만족하는 인자 z를 표준 오차함수 수치표(또는 엑셀 ERFPRECISE)에서 검색하면:

z ≈ 0.67

Step 2. 인자 z 공식 대입

z = cover / ( 2 · √(D_app · T_init) ) = 0.67

0.67 = 50 / ( 2 · √(37.8 · T_init) )

√(37.8 · T_init) = 50 / (2 × 0.67) = 50 / 1.34 ≈ 37.31

37.8 · T_init = (37.31)² ≈ 1,392

T_init = 1,392 / 37.8 ≈ 36.8년

Step 3. 잔존 수명 도출

부식 개시 시점(약 37년)에서 현재 공용 연수(10년)를 차감합니다.

잔존 수명 = 36.8 – 10 = 약 26.8년

즉, 별도의 단면 보수나 표면 발수/도장 코팅을 하지 않아도 향후 약 27년간은 철근 부식 균열 없이 건전성을 유지할 것으로 판정할 수 있습니다.

현장 품질 엔지니어링 체크 요약

최외단 2mm 분말 데이터 배제:
표면 최외각 0~2mm 영역은 대기 강우에 의한 염화물 세척(Washing), 탄산화에 의한 고정 염화물 유리화, 표면 건습 반복으로 인해 순수 확산 거동에서 벗어납니다. 커브 피팅 시 표면 첫 번째 샘플(0~2mm)은 이상치로 제외하고 역해석을 돌려야 정확한 확산 계수를 얻을 수 있습니다.
단순 전염화물 vs 수용성 유리 염화물 구분:
C-S-H 겔이나 프리델염에 고정된 염화물은 부식을 일으키지 않습니다. 엄밀한 잔존 수명 평가를 위해서는 산가용성(전염화물) 프로파일로 확산 계수를 역해석하되, 임계 한계치 도달 여부는 물가용성(유리 염화이온) 시험 결과와 반드시 교차 검증해야 불필요한 과대 보수를 방지할 수 있습니다.