2차원 장거리 상호작용 버지지 노프 모델을 이용한 지진 시뮬레이션 연구

초록

본 논문은 2차원 버지지-노프(BK) 스프링‑블록 모델에 탄성체 인접 효과를 반영한 1/r³ 형태의 장거리 상호작용을 도입하고, 이를 수치 시뮬레이션으로 조사한다. 장거리 모델은 단거리(최근접 이웃) 모델에 비해 시공간 상관이 약해지며, 규모‑분포는 B≈0.55의 고정된 지그프리히‑리히터(B‑값)를 보이는 ‘준임계’ 거동을 넓은 파라미터 구간에서 유지한다. 또한 평균 응력 강하가 규모와 무관하게 일정함을 확인해 실제 관측과 일치시킨다. 큰 사건은 전조 현상(전진진동)과 도넛형 소진 구역을 동반하지만, 사후 진동은 거의 나타나지 않는다.

상세 분석

본 연구는 기존의 2차원 버지지‑노프(BK) 모델이 가정하던 최근접 이웃 상호작용을 탈피하여, 실제 지진 단층이 주변 3차원 탄성체와 연계된 물리적 상황을 반영한다. 저자들은 탄성 이론을 정적 근사(파동 전파 속도가 파열 전파 속도보다 빠른 경우) 하에 풀어, 블록 간 유효 스프링 상수 K(r)∝1/r³ 형태의 장거리 상호작용을 도출하였다. 이는 거리 r이 멀어질수록 힘이 급격히 감소하지만, 최근접 모델에 비해 더 넓은 범위의 블록이 서로 영향을 주고받게 함으로써 시스템의 차원성과 상호작용 범위가 통계적 임계 현상에 미치는 영향을 탐구할 수 있게 한다.

수치 구현에서는 차원 없는 형태로 정규화된 방정식(식 2, 4)을 사용하고, 속도 약화 마찰 법칙 φ(·)을 카를슨‑랑거식(σ, α 파라미터)으로 채택하였다. σ는 정지 마찰에서 동적 마찰으로의 즉각적인 강하량을, α는 속도 증가에 따른 마찰 감소율을 제어한다. 로딩 속도 ν는 실질적으로 0에 가깝게 설정해 지진 발생 중 외부 구동이 없는 조건을 모사하였다. 시뮬레이션은 160×80(또는 60×60) 격자에 대해 10⁵~10⁷개의 이벤트를 생성했으며, 초기 10⁵ 이벤트는 이른 수렴을 위해 폐기하였다.

주요 결과는 규모‑분포 R(µ)에서 나타난다. α가 01 구간에서는 작은 규모(µ≈2 이하)만 발생하지만, α가 1을 초과하면 ‘준임계’ 구간이 넓어져 로그-선형 구간이 B≈0.55(즉, GR 법칙의 b값≈0.85)로 유지된다. 이는 단거리 모델에서 B값이 파라미터에 민감하게 변하고, 특정 α와 l 값에서만 임계 거동을 보이던 것과 대조적이다. 따라서 장거리 상호작용은 규모 분포의 보편성을 강화하고, 실제 지진에서 관측되는 B≈0.50.6과 일치한다는 점에서 중요한 의미를 가진다.

또한 평균 응력 강하 Δσ̄가 규모와 무관하게 일정함을 확인하였다. 이는 관측된 ‘규모 독립적 응력 강하’ 현상을 재현한 것으로, 단거리 모델에서는 큰 사건일수록 Δσ̄가 증가하는 경향을 보였다. 이는 장거리 상호작용이 전체 시스템에 걸친 응력 재분배를 효율적으로 수행함을 시사한다.

시공간 상관 분석에서는 큰 사건 전후에 전조 진동(foreshocks)이 집중되고, 사건 중심부에 도넛형 소진(quiescence) 구역이 형성되는 패턴이 관찰되었다. 반면 사후 진동(aftershocks)은 거의 없으며, 상관 함수 C(t)도 급격히 0에 수렴한다. 이는 장거리 모델이 ‘주요 사건 후 응력 재분배가 빠르게 평형에 도달’한다는 물리적 해석을 가능하게 하며, 실제 지진대에서 관측되는 ‘주요 사건 후 잔여 응력 감소’와 일맥상통한다. 반면 단거리 모델은 강한 잔여 응력과 풍부한 aftershock 시퀀스를 보여, 두 모델 간 차이가 뚜렷하다.

마지막으로 1차원 장거리 모델(1/r² 상호작용)과의 비교를 부록에 제시함으로써 차원 효과를 강조하였다. 1D에서는 여전히 파라미터 튜닝이 필요하지만, 2D 장거리 모델은 차원과 상호작용 범위가 결합되어 보다 넓은 파라미터 공간에서 임계 거동을 보인다. 이는 통계물리학에서 차원과 상호작용 범위가 임계 지수에 미치는 영향과 일치한다.

전반적으로, 장거리 1/r³ 상호작용을 도입한 2D BK 모델은 실제 지진 통계(GR 법칙, 규모 독립적 응력 강하, 전조‑후조 패턴)를 보다 자연스럽게 재현하며, 단거리 모델이 갖는 인공적인 파라미터 민감성을 극복한다는 점에서 지진 물리학 모델링에 중요한 진전을 제공한다.