[논문 리뷰] Scaling in a cellular automaton model of earthquake faults

Marian Anghel, W. Klein외 2명|arXiv (Cornell University)|2000. 02. 29.

earthquake and tectonic studies참고 문헌 3인용 수 5

한 줄 요약

이 논문은 지진 단층 모델이 서로 다른 물리적 메커니즘에 의해 이끌리는 세 가지 구분되는 스케일링 영역을 보임을 제안한다: 소규모 사건은 스핀odal 점 근처에서의 핵형성에 기인하고, 중간 규모 사건은 정지된 핵형성에 기인하며, 대규모 '브레이크아웃' 사건은 정지된 클러스터의 융합에 기인한다. 장거리 응력 전달을 포함하는 세포자동기 모델을 사용하여, 저자들은 스핀odal 유사 임계 변동에 의해 스케일링 행동이 발생함을 보여주며, 이는 1.5와 2.0의 비율 법칙을 따르지만, 가장 큰 사건은 평균장 가정을 초월한 집합적 비평형 역학으로 인해 이격됨을 밝혀냈다.

ABSTRACT

We present theoretical arguments and simulation data indicating that the scaling of earthquake events in models of faults with long-range stress transfer is composed of at least three distinct regions. These regions correspond to three classes of earthquakes with different underlying physical mechanisms. In addition to the events that exhibit scaling, there are larger ``breakout'' events that are not on the scaling plot. We discuss the interpretation of these events as fluctuations in the vicinity of a spinodal critical point.

연구 동기 및 목표

• 장거리 응력 전달을 포함하는 단층 모델에서 지진 스케일링의 물리적 기원을 이해하기 위해.
• 모든 스케일링 사건이 동일한 기초 메커니즘을 공유하는지 여부를 규명하기 위해.
• 작은 사건에서 관찰된 비율 법칙에서 벗어나는 대규모 '브레이크아웃' 사건의 성격을 조사하기 위해.
• 임계점 근처에서 지진 역학이 스핀odal 유사 특이점에 의해 지배되며, 평형 통계역학과 유사하다는 가설을 검증하기 위해.
• 핵형성, 정지 성장, 융합 유도 폭주 파손에 의해 유도되는 사건을 구분하기 위해.

제안 방법

• 2차원 세포자동기 모델은 스프링으로 연결된 블록을 포함하는 단층 시스템을 시뮬레이션하며, 공간적으로 감쇠하는 그린 함수를 통한 장거리 응력 전달을 포함한다.
• 각 블록의 응력은 인접 블록과 로더 플레이트의 기여를 합산하여 계산되며, 응력이 임계값을 초과하면 파손이 발생한다.
• 파손 후 응력은 재분배되며, 비율 K_C/K는 이웃 블록 간에 공유되고, K_L/K는 소산되며, 에너지 방출과 단층의 회복을 모델링한다.
• 지진 모멘트는 사건 동안 파손된 모든 블록의 총 스립으로 정의되며, 사건은 파손 클러스터의 최소 응력에 따라 분류된다.
• 이론적 분석은 랑주아ン 방정식과 안장점 근사법을 사용하여 핵형성과 정지 성장을 모델링하며, 이를 임계 현상과 침투 이론과 연결한다.
• 수치 시뮬레이션은 응력 임계값과 클러스터 크기를 기준으로 사건을 세 가지 유형으로 분류한다: 파손 핵형성 (σ_min ∈ [0.90, 0.99]), 정지 핵형성 (σ_min < 0.90), 브레이크아웃 사건 (크고 비스케일링 성격의 사건).

실험 결과

연구 질문

• RQ1 지진 단층 모델에서 관찰된 세 가지 명백한 스케일링 영역을 뒷받침하는 물리적 메커니즘은 무엇인가?
• RQ2 왜 가장 큰 지진 사건은 작은 사건에서 관찰된 비율 법칙에서 벗어나는가?
• RQ3 핵형성과 정지 성장 과정은 모델에서 관찰된 스케일링 행동에 어떻게 기여하는가?
• RQ4 스케일링 행동이 스핀odal 임계성과 평균장 통계역학으로 얼마나 잘 설명될 수 있는가?
• RQ5 스케일링 도표의 대규모 사건 쪽에서 관찰되는 데이터 산란의 근본 원인은 무엇인가?

주요 결과

• 스케일링 영역은 두 가지 구분된 비율 법칙 영역으로 구성되어 있다: 파손 핵형성 사건에서 기인하는 기울기 ≈1.5의 영역과 정지 핵형성 사건에서 기인하는 기울기 ≈2.0의 영역.
• 파손 핵형성 사건의 평균 크기는 약 1000이며, 이는 예측된 스케일링 1/s^{3/2}와 일치한다.
• 정지 핵형성 사건의 평균 크기는 약 10^4이며, 그 수밀도는 1/s^2로 스케일링되어 이론적 예측과 일치한다.
• 스케일링 도표의 대규모 사건 쪽에서 관찰되는 데이터 산란은 세 번째 유형의 사건—'브레이크아웃' 사건—에 기인하며, 이는 정지된 핵형성 클러스터의 융합으로 발생한다.
• 브레이크아웃 사건은 평균장 스핀odal 이론으로는 설명되지 않으며, 감쇠하는 고응력 클러스터에서의 응력 방출에 의해 유도된 비평형 폭주 파손으로 인해 발생한다.
• 모델은 스핀odal 점 근처에서 시스템의 행동이 평형 통계역학으로 잘 기술되며, 이소프린드 모델과 유체에서 관찰되는 것과 유사한 임계 변동에 의해 스케일링이 지배됨을 확인한다.