Sensitivity Analysis of Shield Tunnel-Induced Ground Settlement to Stratigraphic Spatial Variability: From Random Field Simulation to Causal Decoupling

기술 요약: 지층의 공간적 변동성에 따른 쉴드 터널 유발 지표 침하의 민감도 분석

문제 정의
복잡한 지질 환경, 특히 불규칙한 지층 경계와 뚜렷한 강성 대비를 특징으로 하는 풍화암반 형성층에서의 쉴드 터널 굴착은 지반 침하 예측에 상당한 불확실성을 초래한다. 기존의 결정론적 해석은 대개 균질한 토양 특성을 가정하거나 평균값에 의존하며, 이는 지층 순서의 이산적 변동성과 지층 내 지반 매개변수의 연속적 공간 변동성이라는 지질학적 불확실성의 이중적 성격을 포착하지 못한다. 또한, 기존의 민감도 분석은 전체 토사 프로파일을 하나의 통합된 시스템으로 취급하여, 굴착 지층(변형의 원인)과 상부 덮개층(변형 전달의 경로) 사이의 근본적인 역학적 구분을 간과하는 경우가 많다. 본 연구는 특히 심천 간선 철도(Shenzhen Intercity Railway)의 맥락 내에서, 이러한 공간적 변동성으로부터 발생하는 침하 불확실성을 정량화하고, 설명하며, 예측할 필요성을 다룬다.

연구 방법론
저자들은 랜덤 필드 시뮬레이션, 제어된 수치 실험, 그리고 기호 회귀(symbolic regression)를 결합한 통합 프레임워크를 개발하였다:

1. 수치 모델링 및 검증: 직경 8.0m, 평균 매설 깊이 20m의 EPB(토압식) 쉴드 터널을 모사하는 PLAXIS 3D 기반의 3차원 유한 요소 모델을 구축하였다. 모델에는 토양에 대한 Mohr–Coullet 구성 모델과 구조 부재에 대한 선형 탄성 모델을 사용하였다. 이 모델은 현장 계측 데이터 및 Peck 경험식과 비교하여 검증되었으며, 풍화암반에서의 지표 침하 트러프를 높은 정확도로 재현함을 입증하였다.
2. 이중 불확실성 결합 프레임워크: 지질학적 공간 변동성을 시뮬레이션하기 위해, 연구는 두 가지 모듈 접근 방식을 채택하였다:
   • 지층 시뮬레이션: 결합 마르코프 체인(Coupled Markov Chains, CMC)을 사용하여 11개의 독립적인 지질 구현체를 생성함으로써, 지층 순서(예: 변동하는 풍화 경계면)의 이산적 불확실성을 포착하였다.
   • 매개변수 할당: 각 구현체 내에서 지수 자기상관 함수를 갖는 조건부 랜덤 필드(Conditional Random Fields, CRF)를 적용하여 변형 계수($E$), 점착력($c$), 내부 마찰각($\phi$)의 공간적 가변 값을 할당함으로써 연속적인 매개변수 변동성을 포착하였다.
3. 제어된 매개변수 실험: 지질 매개변수의 효과를 분리하기 위해 37회의 일변수 변화법(one-variable-at-a-time, OVAT) 수치 실험을 수행하였다. 연구는 굴착 지층과 상부 덮개층 각각에 대해 $E$, $c$, $\phi$를 독립적으로 변화시켰으며, 층의 개수를 조절하여 최대 지표 침하($S_{max}$)에 대한 인과적 기여도를 격리하였다.
4. 기호 회귀 및 해석: 매개변수 공간을 보다 균일하게 커버하기 위해 라틴 하이퍼큐브 샘플링(Latin Hypercube Sampling, LHS)을 사용하여 데이터셋을 확장하였다. 이후 기호 회귀를 사용하여 명시적이고 해석 가능한 예측 방정식을 도출하였다. 마지막으로, 결과 모델에 SHAP(SHapley Additive exPlanations) 분석을 적용하여 민감도 순위를 검증하고 물리적 해석 가능성을 확보하였다.

주요 결과

• 공간적 변동성의 영향: 11개의 랜덤 필드 구현체는 $-11.25$ mm에서 $-14.63$ mm 사이의 최대 지표 침하를 생성하였으며, 그 변동 폭은 $3.38$ mm였다. 이는 지질학적 공간 변동성 자체가 실질적인 침하 불확실성을 유발할 수 있음을 확인시켜 주며, 평균값에 기반한 결정론적 모델은 이를 과소평가할 수 있음을 보여준다.
• 민감도 순위: 제어된 실험을 통해 명확한 민감도 계층 구조가 확립되었다: $E_{exc}$(굴착 지층 변형 계수) $> E_{upper}$(상부 덮개층 변형 계수) $> c_{exc}$(굴착 지층 점착력).
  • $E_{exc}$ (35–75 MPa)의 변화는 $3.72$ mm의 침하 변화를 일으켰다.
  • $E_{upper}$ (5–30 MPa)의 변화는 $2.57$ mm의 변화를 일으켰다.
  • $c_{exc}$ (25–48 kPa)의 변화는 $1.09$ mm의 변화를 일으켰다.
  • 내부 마찰각($\phi$)과 지층 복잡성(층 구성)은 침하 변화가 $0.18$ mm 미만으로 나타나 무시할 만한 영향을 미쳤다.
• 역학적 역할: 본 연구는 굴착 지층이 변형의 원인(source) 역할을 하며, 이곳에서는 강성과 더불어 강도가 변형 생성을 지배한다는 것을 확인하였다. 반면, 상부 덮개층은 주로 변형 전달 경로(pathway) 역할을 하며, 여기서는 강성이 유일한 지배 요인이다. 상부 덮개층의 강도 매개변수($c$ 및 $\phi$)는 지표 침하에 측정 가능한 영향을 미치지 않았다.
• 예측 모델: $E_{exc}$, $E_{upper}$, $c_{exc}$를 포함하는 도출된 기호 회귀 방정식은 $R^2 = 0.806$, MAE $= 0.265$ mm, RMSE $= 0.454$ mm를 달성하였다. SHAP 분석을 통한 모델 검증 결과, 민감도 순위($E_{exc} \gg E_{upper} > c_{exc}$)가 재현되어 모델의 물리적 일관성을 확인하였다.

의의 및 주장
본 논문은 쉴드 터널 침하 분석을 단순한 불확실성 특성화를 넘어 메커니즘 기반의 인과적 분리 접근법으로 진전시킨다고 주장한다. 주요 기여는 다음과 같다:

1. 방법론적 혁신: 지층 순서의 불확실성(CMC를 통해)과 지층 내 매개변수 변동성(CRF를 통해)을 동시에 모델링하는 이중 불확실성 프레임워크를 개발함으로써, 결정론적 또는 단일 소스 확률 모델보다 풍화암반 형성을 더욱 현실적으로 표현하였다.
2. 역학적 통찰: 굴착 지층과 상부 덮개층을 명시적으로 분리함으로써, 두 지층이 근본적으로 다른 역학적 역할을 수행함을 밝혀냈다. 이는 토사 프로파일을 단일 단위로 취급하는 기존의 전역 민감도 분석에 도전하며, 매개변수 평가에 있어 더욱 물리적으로 의미 있는 근거를 제공한다.
3. 실무적 유용성: 변형 계수가 점착력이나 내부 마찰각보다 침하를 제어하는 지배적인 매개변수임을 입증하였다. 또한, "블랙박스" 형태의 머신러닝 모델에 대한 투명한 대안으로서 해석 가능한 기호 회귀 방정식을 제안하였다. 이 도구는 물리적 해석 가능성을 희생하지 않으면서도 신속한 예비 침하 예측, 지질 매개변수 스크리닝 및 리스크 평가를 용이하게 한다.

저자들은 본 연구를 복잡한 확률론적 시뮬레이션과 실무적인 공학적 의사결정 사이의 간극을 메우는, 물리 정보 기반의 예측 모델링을 향한 단계로 규정하고 있다.