📊 statistics

A Log-Gaussian Cox process framework for zone-free probabilistic seismic hazard assessment in Sumatra

이 논문은 임의적인 지진원 구역 경계를 제거함으로써 현재의 인도네시아 설계 기준보다 유의미하게 높은 설계 가속도를 산출하는 것을 입증하며, 텍토닉 공변량과 ETAS 여진 모델을 결합한 로그-가우스 콕스 과정을 통해 수마트라를 위한 구역 미지정 확률론적 지진 위험도 평가 프레임워크를 제안한다.

원저자: Irwan Endrayanto Aluicius, Nanang Susyanto, Wiwit Suryanto, Dwi Ertiningsih, Fajar Adi-Kusumo

게시일 2026-09-03

📖 4 분 읽기☕ 가벼운 읽기

CC BY 4.0

원저자: Irwan Endrayanto Aluicius, Nanang Susyanto, Wiwit Suryanto, Dwi Ertiningsih, Fajar Adi-Kusumo

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

지진은 지구상에서 가장 파괴적인 힘 중 하나이며, 수십 년 동안 엔지니어들은 건물이 지진을 견딜 수 있도록 설계하기 위해 특정 방법에 의존해 왔습니다. 확률론적 지진 위험 분석(probabilistic seismic hazard analysis)이라고 알려진 이 방법은 특정 장소에서 장기간에 걸쳐 지면 흔들림이 얼마나 강할지를 예측하려고 시근합니다. 전통적으로 과학자들은 지역을 지도 위에 상자를 그리는 것처럼 크고 고정된 구역으로 나누는 방식으로 접근했습니다. 각 상자 내부에서는 지진 위험이 어디서나 동일하며, 지진이 서로 독립적이고 무작위로 발생한다고 가정합니다. 그러나 거대한 지각판이 충돌하고 맞물려 돌아가는 곳에서는 이러한 가정이 종종 실패합니다. 지하에서 쌓이는 응력은 그려진 상자의 가장자리에서 갑자기 멈추지 않고, 지형을 따라 매끄럽게 흐릅니다. 더욱이, 거대한 지진이 발생하면 그것은 단독으로 일어나지 않습니다. 그것은 단순한 무작위성과는 다른 패턴을 따르는 길고 밀집된 여진의 연쇄 반응을 유발합니다.

인도네시아 수마트라섬에 초점을 맞춘 새로운 연구에서, 연구진은 이러한 경직된 상자 형태를 넘어서는 방법을 제안했습니다. 지진 위험을 미리 정의된 구역 안에 강제로 끼워 넣는 대신, 그들은 지면의 물리적 실체로부터 위험 지도가 스스로 나타나게 하는 유연한 수학적 프레임워크를 사용했습니다. 수십 년간의 지진 기록을 분석하고 이를 지각판이 얼마나 단단하게 맞물려 있는지에 대한 측정값과 결합함으로써, 연구팀은 연속적이고 흐르는 듯한 지진 위험 지도를 만들어냈습니다. 이 접근 방식은 위험이 이전 생각보다 특정 구역에 훨씬 더 집중되어 있음을 드러냈는데, 특히 오랫동안 조용했던 지역이지만 실제로는 엄청난 에너지를 축적하고 있는 지역에서 그러했습니다. 이 결과는 현재 해당 지역의 건축 법규가, 특히 지면이 가장 단단하게 맞물려 있고 미끄러질 준비를 하고 있는 지역에서 발생할 수 있는 흔들림을 상당히 과소평가하고 있을 수 있음을 시사합니다.

연구진은 지구에서 가장 활동적이고 위험한 판 경계 중 하나에 위치한 길고 좁은 섬인 수마트라에 이 새로운 방법을 적용했습니다. 이곳에서는 거대한 인도-호주 판이 순다 판 아래로 파고들며 강렬한 압력을 만들어냅니다. 연구팀은 1998년부터 2024년 사이에 기록된 5,000개 이상의 지진 목록으로 시작했습니다. 진정한 배경 위험을 이해하기 위해, 그들은 먼저 본진과 그 뒤를 따르는 여진을 구분해야 했는데, 이는 마치 본 사건과 그 메아리를 구별하는 것과 같습니다. 일단 본진을 분리한 후, 그들은 물리적 요인들이 어디에서 이러한 지진을 일으키는지 조사했습니다. 그들은 두 가지 핵심 동인을 발견했는데, 하나는 섬을 가로지르는 주요 균열인 그레이트 수마트라 단층(Great Sumatran Fault)으로부터의 거리였고, 다른 하나는 두 지각판이 경계에서 얼마나 단단하게 붙어 있는지를 나타내는 척도인 '지간 결합(inter-seismic coupling)'이었습니다.

이러한 물리적 측정값을 사용하여, 연구팀은 지진 위험을 별개의 구역 모음이 아닌 매끄럽고 보이지 않는 표면으로 취급하는 모델을 구축했습니다. 이 모델에서 위험은 판이 가장 단단하게 맞물려 있고 단층이 가장 가까운 곳에서 가장 높습니다. 이러한 접근 방식은 전통적인 방법이 놓쳤던 패턴을 명확히 보여주었습니다. 예를 들어, 모델은 멘타와이 제도(Mentawai Islands) 해안의 '지진 공백(seismic gap)'을 명확히 식와냈습니다. 이 지역은 2세기 넘게 거대한 파열을 경험하지 못했지만, 모델은 이곳의 판이 매우 단단하게 맞물려 있어 향후 거대 지진이 발생할 위험이 매우 높다고 예측합니다. 이 모델은 과거의 지진 위치에 의존하여 구역을 정의하지 않기 때문에, 물리적 조건(예: 판의 단단한 맞물림)만 적절하다면 오랫동안 조용했던 지역에서도 높은 위험을 예측할 수 있습니다.

연구진이 이 새로운 연속적 위험 지도를 사용하여 건축 설계의 표준 벤치마크인 475년 재현 주기에 대한 예상 지면 가속도를 계산했을 때, 현재의 국가 건축 법규와 매우 다르다는 놀라운 결과를 얻었습니다. 멘타와이 공백 바로 위에 위치한 파당(Padang) 시의 경우, 새 모델은 현재 법규가 요구하는 것보다 거의 3배 더 높은 지면 가속도를 예측했습니다. 반다 아체(Banda Aceh)와 같은 다른 도시들의 경우 차이는 더 작았지만 여전히 유의미했습니다. 연구진은 이러한 차이가 두 가지 구체적인 모델링 선택에서 비롯되었음을 밝혀냈습니다. 즉, 가장 큰 규모의 지진을 제외하지 않고 가능한 모든 규모의 범위를 고려했으며, 지면 흔들림이 어떻게 행동하는지에 대한 전체 통계적 변동성을 포함했는데, 이는 평균 추정치보다 더 높은 정점을 만들어내는 경향이 있습니다.

또한 연구진은 이 새로운 지도가 미래의 지진을 예측할 수 있는지 테스트했습니다. 그들은 1998년부터 2015년까지의 데이터로 모델을 학습시킨 후, 2016년부터 2024년까지의 지진 패턴을 예측하도록 했습니다. 모델은 지진이 발생할 일반적인 위치를 성공적으로 예측했으며, 이는 그들이 사용한 물리적 요인들인 판의 맞물림과 단층 거리가 지면이 흔들릴 가능성이 높은 곳을 가이드하는 신뢰할 만한 지표임을 확인시켜 주었습니다. 그러나 모델은 실제 발생한 것보다 더 많은 총 지진 횟수를 예측했습니다. 연구진은 학습 기간에 거대한 2004년 지진 이후의 이례적으로 높은 활동기가 포함되어 배경 발생률을 일시적으로 부풀렸기 때문이라고 설명했습니다. 지역이 정상적인 속도로 돌아감에 따라 예측되는 발생률은 감소할 것으로 예상되지만, 위험이 가장 높은 공간적 패턴은 안정적으로 유지됩니다.

이러한 발견의 함의는 지진 취약 지역에 대한 우리의 안전 관념에 심오한 영향을 미칩니다. 지도 위에 상자를 그리는 전통적인 방식은 고위험 구역 바로 바깥쪽 지역에 가짜 안도감을 주거나, 물리적으로 파열 준비가 된 조용한 지역에서의 위험을 놓칠 수 있습니다. 지구의 물리 법칙과 함께 흐르는 프레임워크를 사용함으로써, 연구진은 수마트라의 위험이 균일하지도 정적이지도 않다는 것을 보여주었습니다. 가장 큰 위험은 멘타와이 공백에 있으며, 이곳은 수 세기 동안 조용했지만 현재 판 경계에서 가장 단단하게 맞물려 있는 곳입니다. 이 연구는 건축 법규가 이러한 연속적이고 물리 기반의 현실을 반영하도록 재평가되어야 하며, 이를 통해 구조물들이 그 아래 지면의 실제 잠재력을 견딜 수 있도록 설계되어야 함을 시사합니다. 이 작업은 단순히 새로운 지도를 제공하는 것이 아니라, 지진의 위험이 고정된 경계가 아닌 매끄럽고 변화하는 풍경임을 인식하는, 지구를 바라보는 새로운 방식을 제시합니다.

기술 요약: 수마트라의 구역 자유(Zone-Free) 확률론적 지진 재해 평가를 위한 로그-가우시안 콕스 과정 프레임워크

문제 제lement
Cornell (1968)의 프레임워크를 따르는 전통적인 확률론적 지진 재해 분석(PSHA)은 수마트라와 같은 섭입대에서 정당화하기 어려운 두 가지 주요 가정에 의존한다: (1) 지진 발생 영역이 공간적으로 균일한 배경 발생률을 가진 면적 소스 구역으로 분할된다는 점과, (2) 지진 이벤트가 시간적으로 독립적으로 발생한다는 점이다. 실제로 지각 하중은 불연속적인 단계가 아니라 공간적으로 매끄럽게 변화하며, 대규모 본진은 시공간적으로 강력하게 군집된 여진 시퀀스를 생성한다. 저자들은 이러한 가정들이 지진 활동의 물리적 동인을 가릴 수 있다고 주장하며, 공간적 소스 항을 물리적 공변량으로부터 유도된 적합된 확률적 강도로 사용하는 "구역 자유(zone-free)" 대안의 필요성을 제기한다.

방법론
본 연구는 배경 지진 활동을 위한 **로그-가우시안 콕스 과정(Log-Gaussian Cox Process, LGCP)**과 트리거된 군집 현상을 위한 자기 흥분형 에피데믹 유형 여진 시퀀스(Epidemic-Type Aftershock Sequence, ETAS) 모델을 결합한 프레임워크를 제안한다.

데이터 및 전처리:
- 카탈로그: 1998–2024년 사이의 수마트라( $95^\circ\text{E}–108^\circ\text{E}$ , $6^\circ\text{S}–6^\circ\text{N}$ ) 내 $M \ge 4.5$ 이벤트 5,356개를 사용한 작업 데이터셋이 사용되었다.
- 디클러스터링(Declustering): 배경 지진 활동을 격리하기 위해 Gardner–Knopoff (1974) 시공간 윈도우 방법이 적용되었으며, 이를 통해 카탈로그는 693개의 본진으로 축소되었다. 이 단계는 2004–2005년 메가스러스트 여진 시퀀스가 배경 발생률을 10배 가까이 부풀리는 것을 방지하기 위해 필수적이었다.
배경 강도 모델링 (LGCP):
- 배경 강도 $\lambda(s)$ 는 $\log \lambda(s) = \mu(s) + Z(s)$ 인 이중 확률 과정(doubly stochastic process)으로 모델링된다.
- 공변량(Covariates): 평균 추세 $\mu(s)$ $μ (s)$ 는 두 가지 물리적으로 해석 가능한 지각 공변량의 로그-선형 함수이다:
  - 지각 간 결합(Inter-seismic coupling, $\phi$ ): 지질 구조 역산으로부터 도출된 순다 메가스러스트 인터페이스의 결합 비율로, 여전히 잠겨 있는 멘타와이(Mentawai) 지진 갭(seismic gap)에서 정점을 찍는다.
  - 단층 거리( $D_{fault}$ ): 그레이트 수마트라 단층(GSF)으로부터의 수직 거리.
- 잠재장(Latent Field): $Z(s)$ 는 공변량으로 설명되지 않는 잔여 미세 규모 변동성을 포착하기 위해 지수 공분산 구조를 가진 가우시안 랜덤 필드(GRF)이다.
- 추정: 파라미터는 $60 \times 60$ 그리드 상에서 포아송 로그-선형 최대 우도 추정법(MLE)을 통해 추정되었다.
트리거 성분:
- 트리거된 군집 성분은 별도로 기술되는 자기 흥분형 ETAS 모델에 의해 설명되며, 오모리-우츠(Omori–Utsu) 시간 커널을 사용한다. 이 모델은 공변량 기반 배경 모델이 설명하지 못하는 이벤트 간 트리거링을 고려한다.
재해 통합:
- 적합된 LGCP 강도 표면은 코넬 재해 적분식에서 직접 공간적 소스 항으로 사용되어 구역 구획의 필요성을 제거한다.
- 연간 초과 발생률인 최대 지반 가속도(PGA)는 구텐베르크-리히터(Gutenberg–Richter) 규모 분포와 Youngs 등(1997)의 섭입-인터페이스 지반 운동 예측 방정식(GMPE)을 사용하여 계산된다.
검증 및 비교:
- 모델 선택: LGCP를 물리적 공변량이 없는 통계적 베이스라인인 네이만-스콧 토마스(Neyman–Scott Thomas) 클러스터 과정과 비교하였으며, AIC(Akaike Information Criterion)와 비균질 L-함수(inhomogeneous L-function)를 사용하였다.
- 예측 검증: 공간적 예측 성능을 평가하기 위해 시간적 홀드-아웃 테스트(1998–2015년 학습, 2016–2024년 테스트)와 지역 지진 발생 가능성 모델(RELM) 테스트 스위트를 사용하였다.

주요 결과

모델 성능: 공변량 기반 LGCP는 우도비 개선도 $2\Delta\ell = 45.7$ ( $\chi^2_2, p < 10^{-4}$ ) 및 $\Delta$ AIC 643.7을 기록하며 토마스 클러스터 베이스라인을 결정적으로 능가하였다. LGCP는 사전 정의된 경계 없이도 관측된 재해 구배(hazard gradient), 즉 멘타와이 갭을 포함한 고위험 지역을 성공적으로 재현하였다.
공변량 유의성: 두 공변량 모두 매우 유의미하였다 (Cox–Snell pseudo- $R^2 \approx 0.61$ ). 결합 계수는 양수였고(높은 결합 $\rightarrow$ 높은 발생률), 단층 거리 계수는 음수였다(GSF에 가까울수록 $\rightarrow$ 높은 발생률).
재해 추정치: 본 프레임워크는 인도네시아 국가 설계 기준(SNI 1726:2019)보다 상당히 높은 475년 재현 주기 설계 가속도를 산출하였다. LGCP 유도 PGA와 코드 PGA의 비율은 다음과 같다:
- 파당(Padang): 2.7배 (1.79 g)
- 페칸바루(Pekanbaru): 2.5배
- 메단(Medan): 2.1배
- 반다 아체(Banda Aceh): 1.5배
공간적 패턴: 가장 높은 재해는 북서쪽과 파당 회랑을 따라 집중되었으며, 이는 높은 결합력을 가진 멘타와이 갭과 직접적으로 일치한다. 반면, 반다 아체는 2004년 파열 이후 낮은 잔여 결합력으로 인해 낮은 재해를 보였다.
검증: 시간적 홀드-아웃 테스트는 LGCP의 공간적 패턴이 미학습 데이터에서도 보존됨을 확인하였다 (L-test 퀀타일 $\gamma = 0.21$ ). 그러나 절대적인 발생률은 테스트 기간 동안 3.3배 과대 예측되었는데, 이는 학습 창에 2004년 이후 상승한 여진 활동이 포함되었기 때문으로 분석된다.

의의 및 주장
본 논문은 적합된 공변량 기반 콕스 강도와 코넬 재해 적분 사이의 직접적인 수학적 결합을 입증했다고 주장하며, 이는 섭입대 환경에서 이전에 입증된 바 없는 조합이다. 주요 의의는 사전 정의된 구역이 아닌 연속적인 확률적 소스 모델로부터 설계 지반 운동을 도출할 수 있다는 점에 있다.

저자들은 본 프레임워크가 다음을 수행함을 강조한다:

물리적 구배 재현: 입력된 파열 이력이나 구역 경계 없이도, 오직 물리적 공변량(결합 및 단층 거리)만을 사용하여 멘타와이 지진 갭과 수마트라의 재해 구배를 포착한다.
보수적 가이드 제공: 국가 코드보다 높은 재해 값을 산출한 것은 두 가지 구체적이고 방어 가능한 모델링 선택에 기인한다: (a) 카탈로그와 일치하는 규모 완결성 하한( $M_c=4.5$ )을 사용하여 큰 이벤트에 더 높은 가중치를 부여한 점, (b) GMPE에 절단되지 않은 무작위 변동성(aleatory variability) 항을 적용한 점.
시간 척도 구분: 저자들은 LGCP가 지각 하중이 발생률을 결정하는 장기적 공학 설계(수십 년/수백 년)에 적합한 반면, 자기 흥행형 ETAS 모델은 트리거된 군집 현상이 지배적인 단기 운영 예측(수일/수개월)에 더 적합하다고 명시한다.

결론적으로, 본 논문은 구역 자유 LGCP 프레임워크가 특히 지각 하중이 매끄럽게 변화하거나 역사적 파열 데이터가 부족한 지역(예: 잠겨 있는 멘타와이 갭)에서 전통적인 면적 소스 PSHA에 대한 지질학적으로 타당하고 재현 가능한 대안을 제공한다고 주장한다.

기술 요약: 수마트라의 구역 자유(Zone-Free) 확률론적 지진 재해 평가를 위한 로그-가우시안 콕스 과정 프레임워크

유사한 논문