2차원 버리게 노프 모델에서 관측된 Weibull‑log Weibull 전이와 지진 발생 간격 통계

초록

2차원 Burridge‑Knopoff 모델로 생성한 인공 지진 카탈로그를 분석한 결과, 연속 지진 사이의 발생 간격(Inter‑occurrence time)은 Weibull 분포와 log‑Weibull 분포의 선형 결합으로 잘 설명된다. 진폭 임계값을 높이면 log‑Weibull 영역에서 Weibull 영역으로 전이하는 ‘Weibull‑log Weibull 전이’ 현상이 나타나며, 이 전이는 마찰계수와 판 경도와 같은 모델 파라미터에 민감하게 반응한다. 연구는 이 모델이 실제 지진 데이터의 복합 통계적 특성을 재현할 수 있는 물리적 프레임워크임을 제시한다.

상세 분석

본 논문은 2차원 Burridge‑Knopoff (BK) 모델을 이용해 대규모 합성 지진 카탈로그를 생성하고, 사건 간 발생 간격인 inter‑occurrence time (IOT)의 확률분포를 정량적으로 분석하였다. 기존 연구에서는 IOT가 단일 지수형 혹은 파워‑law 형태로 기술된 경우가 많았지만, 저자들은 Weibull 분포와 log‑Weibull 분포의 가중합 형태가 데이터에 가장 적합함을 확인하였다. Weibull 분포는 꼬리가 짧고 급격히 감소하는 형태로, 큰 진폭을 가진 이벤트 사이의 간격을 잘 포착한다. 반면 log‑Weibull 분포는 로그 변환 후 Weibull 형태를 띠어, 작은 진폭 이벤트가 다수 존재하는 영역에서 긴 꼬리를 제공한다.

모델 파라미터인 마찰계수 α와 판의 탄성계수 kx, ky를 변화시켜 다양한 물리적 조건을 모사했으며, 그 결과 IOT 분포는 α와 k값에 따라 가중치 비율이 크게 달라졌다. 특히, 진폭 임계값 Mth를 높이면 log‑Weibull 성분의 가중치가 급격히 감소하고 Weibull 성분이 우세해지는 ‘Weibull‑log Weibull 전이’가 관찰되었다. 이는 큰 지진이 발생할수록 사건 간 간격이 보다 규칙적이고, 통계적으로 Weibull 형태에 가까워진다는 물리적 해석을 가능하게 한다.

또한, 전이 현상이 단순히 통계적 피팅에 의한 것이 아니라, 마찰 비선형성 및 판 경도에 의해 발생하는 응력 누적과 해제 메커니즘과 직접 연결된다는 점을 강조한다. 마찰계수가 클수록 응력 해소가 급격히 일어나 Weibull 지수 β가 증가하고, 판 경도가 높을수록 시스템이 보다 강직하게 동작해 로그 변환 효과가 감소한다. 이러한 결과는 실제 지진 관측 데이터에서 보고된 Weibull‑like 및 log‑Weibull‑like 특성을 모델이 물리적으로 재현할 수 있음을 시사한다.

마지막으로, 저자들은 기존의 단일 분포 모델이 갖는 한계점을 지적하고, 복합 분포 접근법이 지진 발생 간격의 다중 스케일 특성을 포착하는 데 유리함을 입증하였다. 향후 연구에서는 3차원 모델 확장, 비균질 마찰 특성 도입, 그리고 실제 관측 카탈로그와의 정량적 비교를 통해 전이 메커니즘을 더욱 정교화할 필요가 있다.