Omori's Law Applied to Mining-Induced Seismicity and Re-entry Protocol Development

被引:21
作者
Vallejos, J. A. [1 ]
McKinnon, S. D. [1 ]
机构
[1] Queens Univ, Dept Min Engn, Kingston, ON K7L 3N6, Canada
关键词
Omori's law; Mining-induced seismicity; Aftershock sequences; Mine safety; Ground control; Re-entry protocol; AFTERSHOCK SEQUENCES; EARTHQUAKE CATALOGS; COMPLETENESS; MAGNITUDE; FORMULA; MODELS;
D O I
10.1007/s00024-009-0010-7
中图分类号
P3 [地球物理学]; P59 [地球化学];
学科分类号
0708 ; 070902 ;
摘要
This paper describes a detailed study of the Modified Omori's law n(t) = K/(c + t) (p) applied to 163 mining-induced aftershock sequences from four different mine environments in Ontario, Canada. We demonstrate, using a rigorous statistical analysis, that this equation can be adequately used to describe the decay rate of mining-induced aftershock sequences. The parameters K, p and c are estimated using a uniform method that employs the maximum likelihood procedure and the Anderson-Darling statistic. To estimate consistent decay parameters, the method considers only the time interval that satisfies power-law behavior. The p value differs from sequence to sequence, with most (98%) ranging from 0.4 to 1.6. The parameter K can be satisfactorily expressed by: K = kappa N (1), where kappa is an activity ratio and N (1) is the measured number of events occurring during the first hour after the principal event. The average kappa values are in a well-defined range. Theoretically kappa a parts per thousand currency sign 0.8, and empirically kappa a [0.3-0.5]. These two findings enable us to develop a real-time event rate re-entry protocol 1 h after the principal event. Despite the fact that the Omori formula is temporally self-similar, we found a characteristic time T (MC) at the maximum curvature point, which is a function of Omori's law parameters. For a time sequence obeying an Omori process, T (MC) marks the transition from highest to lowest event rate change. Using solely the aftershock decay rate, therefore, we recommend T (MC) as a preliminary estimate of the time at which it may be considered appropriate to re-enter an area affected by a blast or large event. We found that T (MC) can be estimated without specifying a p value by the expression: T (MC) = a N (1) (b) , where a and b are two parameters dependent on local conditions. Both parameters presented well-constrained empirical ranges for the sites analyzed: a a [0.3-0.5] and b a [0.5-0.7]. These findings provide concise and well-justified guidelines for event rate re-entry protocol development.
引用
收藏
页码:91 / 106
页数:16
相关论文
共 37 条
[1]   NEW LOOK AT STATISTICAL-MODEL IDENTIFICATION [J].
AKAIKE, H .
IEEE TRANSACTIONS ON AUTOMATIC CONTROL, 1974, AC19 (06) :716-723
[2]   A TEST OF GOODNESS OF FIT [J].
ANDERSON, TW ;
DARLING, DA .
JOURNAL OF THE AMERICAN STATISTICAL ASSOCIATION, 1954, 49 (268) :765-769
[3]  
Board M, 2001, UNDERGROUND MINING METHODS, P483
[5]  
DIETERICH JH, 1986, EARTHQUAKE SOURCE ME, V37, P37
[6]   IDENTIFICATION OF AFTERSHOCKS OF DEEP EARTHQUAKES BY A NEW RATIOS METHOD [J].
FROHLICH, C ;
DAVIS, S .
GEOPHYSICAL RESEARCH LETTERS, 1985, 12 (10) :713-716
[7]  
GROSS SJ, 1994, B SEISMOL SOC AM, V84, P1571
[8]  
HABERMANN RE, 1984, B SEISMOL SOC AM, V74, P1743
[9]  
HEAL D, 2005, SEISMIC HAZARD AGNIC
[10]   Subcritical and supercritical regimes in epidemic models of earthquake aftershocks [J].
Helmstetter, A ;
Sornette, D .
JOURNAL OF GEOPHYSICAL RESEARCH-SOLID EARTH, 2002, 107 (B10)