⚡ electrical engineering

Seismic Stability of Karst Sinkholes under high peak ground acceleration: a three dimensional numerical study in Barangay Dontogan, Baguio City, Philippines

이 3차원 수치 연구는 높은 최대 지반 가속도 조건 하에서 바기오 시 돈토간 마을에 위치한 44개의 검증된 싱크홀에 대한 지진 안정성을 평가하며, 얕은 공동들은 0.40 g에서 안정적으로 유지되는 반면, 높이 대 직경 비율이 1.00 이하인 공동들의 경우 0.80 g의 설계 요구가 붕괴 및 임계 상태로의 분기를 유발한다는 점을 밝혀내어 기초 안전성 평가를 위한 4단계 안정성 임계치 매트릭스를 제안한다.

원저자: Deryl Jhon Lucas

게시일 2026-09-15

📖 3 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Deryl Jhon Lucas

원본 논문은 CC BY 4.0 (https://creativecommons.org/licenses/by/4.0/) 라이선스로 제공됩니다. ✨ 이것은 아래 논문에 대한 AI 생성 설명입니다. 저자가 작성하거나 승인한 것이 아닙니다. 기술적 정확성을 위해서는 원본 논문을 참조하세요. 전체 면책 조항 읽기

필리핀의 험준한 고지대에서 우리가 발을 딛고 있는 지표면은 항상 견고한 것만은 아닙니다. 수천 년에 걸쳐 고대의 빗물이 석회암을 서서히 녹여낸 곳에서는 지하에 숨겨진 동굴과 빈 공간의 네트워크가 형성됩니다. 이러한 빈 공간의 천장이 상부의 무게를 지탱하기에 너무 얇아지면, 지표면이 갑자기 무너지며 싱크홀이 발생합니다. 이것들은 단순히 깊은 구덩이가 아니라, 예고 없이 집과 도로를 삼켜버릴 수 있는 갑작스럽고 파괴적인 붕괴입니다. 위험은 지진이 발생할 때 증폭됩니다. 지진은 단순히 건물을 흔드는 것이 아니라 지표면 자체를 요동치게 하며, 이는 간신히 안정적인 상태를 유지하던 암석 천장을 붕괴의 끝으로 밀어 넣을 수 있는 강력하고 리드미컬한 힘을 가합니다. 엔지니어와 도시 계획가들에게 있어, 특히 지반이 이미 불안정한 것으로 알려진 지역에서 그 천장이 언제 안전하고 언제 붕괴 직전인지를 정확히 아는 것은 오랫동안 난제였습니다.

서늘한 기후와 가파른 지형으로 유명한 산악 휴양지인 바기오 시(Baguio)에서 이 문제는 매우 심각합니다. 이 도시는 균열이 있는 석회암 지층 위에 위치해 있으며, 특히 돈토간(Dontoggan) 마을은 이 지역에서 기록된 싱크홀이 가장 많이 집중되어 있습니다. 지역 건축 법규는 지진을 견딜 수 있도록 구조물을 설계하는 엄격한 규칙을 제공하지만, 숨겨진 지하 공동 바로 위에 위치한 건물들에 대한 구체적인 지침은 제공하지 않습니다. 이 공백을 메우기 위해, 유니버시티 오브 더 코디예라(University of the Cordilleras)의 한 연구자는 싱크홀 공동이 대규모 지진의 강한 흔들림을 받을 때 어떤 일이 발생하는지 시뮬레이션하기 위해 첨단 컴퓨터 모델링을 활용했습니다. 그들의 목표는 단순히 지표면이 갈라지는 것을 보는 것이 아니라, 안정적인 천장이 붕괴하는 천장으로 변하는 정확한 임계점을 결정하고, 그 발견을 공동체를 보호하기 위한 실질적인 도구로 변환하는 것이었습니다.

연구자는 공동의 기하학적 구조, 구체적으로는 공동을 덮고 있는 토양과 암석의 두께와 공동 자체의 너비 사이의 관계에 집중했습니다. 그들은 돈토간의 실제 조건을 모방하여, 흙 층과 절리된 석회암 층 아래에 묻힌 구형 구멍의 3차원 디지털 모델을 제작했습니다. 그런 다음 이 모델에 시뮬레이션된 지진력을 가했으며, 흔들림의 강도는 이 지역의 높은 지진 요구 수준인 0.80g까지 높였습니다. 이는 중력 가속도와 거의 맞먹는 힘을 나타내는 지면 가속도의 척도입니다. 연구자는 다양한 크기의 공동과 덮개 두께에 걸쳐 시뮬레이션을 실행함으로써, 흔들림의 응력에 따라 암석과 토양이 어떻게 반응하는지 관찰할 수 있었습니다.

결과는 명확하면서도 다소 경각심을 불러일으키는 패턴을 보여주었습니다. 흔들림이 중간 정도일 때는 구멍의 크기에 상관없이 지표면이 안정적으로 유지되었습니다. 그러나 시뮬레이션된 지진의 강도가 높아짐에 따라 천장의 안정성은 급격히 떨어졌습니다. 가장 결정적인 발견은 가장 높은 설계 수준의 흔들림이 발생했을 때 나타났는데, 연구자는 특정 기하학적 임계점을 식서했습니다. 만약 공동 위의 암석과 토양 덮개가 공동의 너비보다 얇다면, 시스템은 치명적으로 불안정해졌습니다. 이러한 시나리오에서 컴퓨터 모델은 암석 천장이 하중을 견디지 못하고 붕괴로 이어진다는 것을 보여주었습니다. 시뮬레이션 결과, 돈토간의 특정 조건에서 공동보다 얇은 천장은 주요 지진 발생 시 붕괴할 가능성이 높은 반면, 공동보다 현저히 두꺼운 천장은 안전한 것으로 나타났습니다.

이 복잡한 연구 결과를 일상적인 의사결정에 유용하게 사용하기 위해, 연구자는 계산된 안전 마진에 따라 위험을 네 가지 범주로 나누는 간단한 분류 체계를 개발했습니다. '안정(stable)' 등급은 표준 콘크리트 슬래브나 단일 말뚝과 같은 일반적인 기초로 충분함을 의미합니다. '경계(marginal)' 등급은 공동 아래의 단단한 암반 깊숙이 건물을 고정하는 그라우팅 마이크로 파일과 같은 더 강한 지지력이 필요함을 시사합니다. '위험(critical)' 등급은 공동 전체를 가로지르는 래프트 형태의 연결된 말뚝과 같은 특수 솔루션이 요구될 만큼 위험이 높음을 나타냅니다. 마지막으로 '붕괴(collapse)' 등급은 해당 지점이 건축하기에 너무 위험하며, 유일한 안전한 선택지는 부지를 피하거나 공동 훨씬 아래의 단단한 암반을 통해 건물의 무게를 전달하는 것임을 의미합니다.

이 새로운 프레임워크는 현재 돈토간에 매핑된 44개의 검증된 싱크홀에 적용되었습니다. 분석 결과, 건물의 안전은 주변 지역의 안정성에 근거하여 가정될 수 없으며, 모든 개별 부지는 반드시 개별적으로 평가되어야 함이 드러났습니다. 오늘날 안전해 보이는 부지라도 대규모 지진이 발생한다면, 특히 지하 공동의 기하학적 구조가 위험 또는 붕괴 구역에 해당한다면 붕괴 직전의 상태일 수 있습니다. 이 연구는 이 지역에서 정적인 안정성에 의존하는 전통적인 방식이 불충분하다는 점을 확인시켜 주며, 지진의 동적인 힘이 결정적인 요소임을 입증합니다. 이 위험을 평가하는 명확하고 정량적인 방법을 제공함으로써, 이 연구는 시 관계자와 엔지니어들에게 미래의 재난을 예방하고, 새로운 구조물이 진정으로 지표면이 버틸 수 있는 곳에만 세워질 수 있도록 하는 중요한 도구를 제공합니다.

기술 요약: 높은 최대 지반 가속도 하에서의 카르스트 싱크홀의 지진 안정성

문제 정의
필리핀 바기오 시의 Barangay Dontogan은 카르스트 침하의 높은 위험에 직면해 있으며, 129개의 지역 인벤토리 중 가장 많은 44개의 검증된 싱크홀이 집중되어 있다. 이 싱크홀들은 개구부(aperture) 크기가 1.0에서 50.0미터에 이르며, 절리가 발달한 Kennon 석회암층 위에 위치한다. 국가 구조 설계 기준인 필리핀 구조 설계 기준(NSCP)은 Zone 4(근원지 설계 요구량 약 0.80 g ~ 0.83 g)에 대한 지진 설계 규정을 제공하지만, 지진 하중 하에서의 지하 공동(cavity) 안정성을 평가하기 위한 구체적인 기준은 결여되어 있다. 기존의 암반 분류 체계와 결정론적 해석 방법은 공동을 포함하는 절리 암반의 3차원 비선형 거동을 반영하지 못하며, 이는 이 고위험 지역 내 주거 개발을 위한 인허가 단계의 의사결정 과정에서 규제 공백을 야기한다.

연구 방법론
본 연구는 PLAXIS 3D를 이용한 3차원 유한요소법(FEM)을 사용하여 구형 카르스트 공동의 지진 응답을 모델링한다.

모델 기하 구조: 2층 구조의 층서 프로파일을 구축하였다: Mohr–Coulomb 모래/자갈 성분의 실트-클레이 상부층이 일반화된 Hoek–Brown 절리 Kennon 석회암층 위에 놓인 형태이다. 구형 공동은 석회암-상부층 경계면이 공동의 중간 높이를 통과하도록 배치되었다.
파라미터 스윕(Parametric Sweep): 5가지 높이-직경(H/D) 비율(0.50, 0.75, 1.00, 1.50, 2.25)과 3가지 최대 지반 가속도(PGA) 단계(0.40 g, 0.60 g, 0.80 g)를 교차하여 15개의 시나리오 매트릭스를 생성하였다. 0.80 g 단계는 해당 지역의 NSCP 2015 설계 요구량을 나타낸다.
지진 입력: 분석에는 1990년 Mw 7.8 Luzon 지진의 시계열 가속도 기록을 사용하였으며, 이를 목표 PGA로 스케일링하였다. 본 연구는 ASCE 7-22의 최소 기준인 11개의 기록 대신 단일 기록을 사용한 한계를 인정하며, 결과가 특정 기록에 종속적임을 명시한다.
분석 절차: 시뮬레이션은 4단계 파이프라인을 따른다: (1) 단계별 시공 및 정적 평형; (2) 동적 시계열 하중; (3) 진동 후 잔류 재평형; (4) 강도 감소법(SRM)을 통한 잔류 극한 평형 상태에서의 동적 안전율( $\Sigma M_{sf}$ ) 추출.
검증: 모델은 격자 독립성(mesh independence)에 대해 검증되었으며, 기존의 3D 극한 분석 연구들과 벤치마킹되었다. 안정성-붕괴 전이 구간을 확정하기 위해 중간 PGA(0.55 g, 0.65 g, 0.75 g)에서의 이산 이벤트 검증 실행이 수행되었다.

주요 결과
본 연구는 다양한 지진 요구 조건에 따른 싱크홀 구성의 안정성을 정량화하였다:

0.40 g 단계: 5가지 H/D 구성 모두 안정(Stable) 상태를 유지하였다 ( $\Sigma M_{sf}$ 범위: 1.642–1.683). 안정성 여유는 기하학적 구조에 약하게 반응하였다.
0.60 g 단계: 안정성이 크게 저하되었다. $\Sigma M_{sf}$ 대역이 1.007–1.210으로 수렴하였으며, 이는 0.40 g 단계와 비교했을 때 H/D에 따른 편차가 4.95배 증가했음을 의미한다. 대부분의 구성이 한계(Marginal) (1.10–1.30) 또는 위험(Critical) (1.00–1.10) 구역에 진입하였다.
0.80 g 단계 (설계 기준): 스윕 결과 두 개의 뚜렷한 분기로 갈라졌다:
- 붕괴 분기(Collapse Branch): 두 가지 구성이 실패하였다. H/D = 0.50은 동적 붕괴를 일으켰다 ( $\Sigma M_{sf}$ = 0.9313). H/D = 1.00은 정적 1단계(사전 붕괴 발생)에서 실패하였으며 ( $\Sigma M_{sf}$ = 0.9092), 이는 격자 독립성이 확인되었다.
- 위험 분기(Critical Branch): 세 가지 구성(H/D = 0.75, 1.50, 2.25)은 불안정하지만 안정적인 상태를 유지하였으며, $\Sigma M_{sf}$ 값은 1.001에서 1.006 사이였다.
식별된 임계치: 본 연구는 $\Sigma M_{sf}$ = 1.00의 기계적 임계치와 0.80 g 설계 요구량에서의 임계 기하학적 범위 H/D $\le$ 1.00을 식별하였다. 안정성에서 붕로로의 전이는 0.75 g와 0.80 g 사이에서 발생한다.

기여 및 산출물
본 연구의 주요 기여는 Barangay Dontogan에 특화된 **안정성 임계치 매트릭스(Stability Threshold Matrix)**를 개발한 것이다. 이 도구는 추상적인 임계치( $\Sigma M_{sf}$ = 1.00 및 H/D = 1.00)를 네 가지 실행 가능한 판정과 기초 권고 사항으로 변환한다:

안정 (Stable, $\Sigma M_{sf} \ge$ 1.30): 일반적인 얕은 기초 또는 단일 말뚝 기초.
한계 (Marginal, 1.10 $\le$ $\Sigma M_{sf}$ < 1.30): 그라우팅된 마이크로 파일 기초.
위험 (Critical, 1.00 $\le$ $\Sigma M_{sf}$ < 1.10): 마이크로 파일 래프트 기초 또는 건설 회피.
붕괴 (Collapse, $\Sigma M_{sf}$ < 1.00): 건설 회피 또는 공동을 통한 공학적 하중 전달.

이 매트릭스는 44개의 검증된 싱크홀 인벤토리에 적용되었으며, 이를 통해 정성적 구역 설정(zoning)을 넘어 부지별 footprint 단위의 스크리닝 도구를 제공하였다.

의의 및 주장
본 논문은 필리핀의 카르스트 싱크홀에 대해 최초의 정량적, 지진 특화적 의사결정 지원 도구를 제공함으로써 중요한 규제 공백을 메운다고 주장한다. 본 연구는 정적 상부 하중보다 지진 요구량이 해당 지역의 얕은 카르스트 지붕 안정성을 결정짓는 핵심 요소임을 입증하였다. 즉, 정적 조건에서는 안정적으로 보이는 구성이라도 높은 PGA 하에서는 실패할 수 있다. 연구는 0.80 g 설계 요구량에서의 안정성 여유가 매우 얇으므로, 주변 지반의 역사적 안정성에 의존하기보다는 모든 제안된 건축물 부지에 대해 개별적인 수치적 스크리닝이 필요함을 강조한다. 저자는 안정성 임계치 매트릭스를 향후 NSCP 개정판의 모듈로 제안하며, MGB-CAR의 정성적 취약성 지도와 구조 엔지니어 및 건축 행정관의 정량적 요구 사항 사이의 가교 역할을 한다고 본다. 저자는 결과의 일반화 가능성에 대해서는 신중한 태도를 유지하며, 본 임계치들이 특정 기록에 종속적이고 Barangay Dontogan의 특정 지질에 맞춰 보정된 것이므로 다른 카르스트 환경에서는 재보정이 필요함을 명시한다.

유사한 논문