地震活动具有复杂性，用传统的线性动力学理论来刻画复杂系统已显得越来越局限，而用非线性科学的理论去描述和研究复杂系统已成为近年来的热点。R/S时间序列分析方法作为非线性时间序列分析方法之一自1965年由Hurst提出以来，已被广泛应用于研究分子布朗运动(FBM)和自然现象的分形特征，R/S时间序列分析方法在地震学及地震预报中也有了一定的应用^[1-4]。这些研究都各有特点，各有侧重，而关于地震震级变化趋势的分析大多是基于统计学，用平均值、均方差等统计特征量描述某地区过去不同震级地震发生的总体情况和涨落等特征^[5-6]。大地震对人类的威胁极大，研究地震震级的变化趋势对地震研究和安全研究工作都有参考作用。笔者应用赫斯特指数(Hurst exponent)方法分析全球7.0级以上地震震级的变化趋势，旨在探讨地震震级变化趋势异常表现形态的客观性。

1 数据来源及震级总体趋势

从文献[7]取得全球7.0级以上地震震级数据，数据时期为1898年4月15日7点12分的7.3级地震到2011年1月3日4点20分18秒的智利中部7.1级地震，共计2 194个震级数据。数据总体趋势见图 1，从左至右是1898年往2011年方向，X为地震发生顺序值，Y为相应震级，趋势线为Y=-0.000 07X+7.399，斜率为-0.000 07，即全球7.0级以上地震震级趋势有10万分之7的下降趋势，与“地球自转长期减慢”^[8]相符合(地球自转影响地质构造应力及于影响地震震级)。

2 Hurst指数及相关概念

H.E.Hurst(1900—1978)在研究尼罗河水库水流量和贮存能力的关系时，发现用有偏的随机游走(分形布朗运动)能更好地描述水库的长期存贮能力，并用重标极差(R/S)分析方法建立Hurst指数(H)，作为判断时间序列数据遵从随机游走还是有偏的随机游走过程的指标。

2.1 重标极差法

重标极差法(rescaled range analysis，R/S)由Hurst提出，后经Mandelbrot补充和完善^[9-13]，已经广泛地应用于判断时序的分形特征及长程相关性，区分随机与非随机系统，并确定趋势性的持续和强度。基本分析原理如下。

设震级时间序列为X=X₁, X₂，…，X_n，则X序列的均值和标准差为：

$ \overline X = \frac{{\left( {{X_1} + \cdots + {X_n}} \right)}}{n}, $

(1)

$ S=\sqrt {\frac{1}{{n - 1}}\sum\limits_{i = 1}^n {{{\left( {{X_i} - \overline X} \right)}^2}}, } i = 1, 2, \cdots, n。$

(2)

对数据进行标准化处理得序列Z，Z=0：

$ {Z_i} = \frac{{{X_i} - \overline X}}{S}, i = 1, 2, \cdots, n, $

(3)

由Y₁=Z₁-Z_i, i=2, …, n；Y₂=Z₂-Z_i, i=3, …, n; …Y_n-1=Z_n-1-Z_n得(n-1)个时间数据序列Y_i，i=1, 2, …, n-1，计算极差序列R₁, …, R_n-1：

$ {R_i} = \max \left( {{Y_i}} \right) - \min \left( {{Y_i}} \right), i = 1, 2, \cdots, n - 1。$

(4)

计算Y_i的标准差S_i，样本数为n_i，并计算重标极差R/S，得${\left( {\frac{R}{S}} \right)_i}$, i=1, 2, …, n-1。

不同时间尺度对应不同的重标度极差(R/S)_n。Mandelbrot证明log(R/S)_n与log(n)之间存在线性关系：

$ \log {\left( {\frac{R}{S}} \right)_n} = H \cdot \log \left( n \right) + \log \left( b \right)。$

(5)

以log(n)为自变量、log(R/S)_n为因变量作双对数坐标散点图，采用最小二乘估计拟合直线，直线斜率H值即为Hurst指数，b为常数，表征系统固有特性。

2.2 V统计量

V统计量是由Hurst提出，后经Peters改进，用于测量序列“长期记忆”的长度，其定义式为^[11-14]：

$ {V_n} = \frac{{\left( {\frac{{{R_n}}}{{{S_n}}}} \right)}}{{\sqrt n }}。$

(6)

其特征为：1) 独立随机过程的时间序列，V_n关于log(n)的曲线是一条水平直线。2) 如果序列具有状态持续性，即当H＞0.5时，V_n关于log(n)的曲线向上倾斜。3) 如果序列具有逆状态持续性，即当H＜0.5时，V_n关于log(n)的曲线向下倾斜。4) 当V_n关于log(n)的曲线形状改变时，就产生了突变，长期记忆消失。用V统计量曲线可以直观地看出一个时间序列某一时刻的值对以后值的影响时间的界限。

2.3 相关系数C和分形维D

为了度量现在对未来的影响，Mandelbrot于1968年引入一个相关性度量指标^[11-14]：

$ C = {2^{\left( {2H - 1} \right)}} - 1。$

(7)

时间序列分形维D与Hurst指数有以下关系^[15-16]：

$ D = 2 - H。$

(8)

时间序列的分形维在1和2之间。随机游动的累积变化的分形维数是1.50，一条线的分形维数是1.0，一个几何平面的分形维数是2.0。

1) 当C=0，H=0.5，D=1.5，是布朗运动的Hurst指数和维数，变量之间是相互独立的，相应的相关系数是零，现在不会影响未来，因此时间序列是随机的。

2) 当-0.5 ≤C＜0，0 ≤H＜0.5，D＞1.5，变量之间是负相关的，称之为反持续性的、遍历性的或均值回复的时间序列。如果某一时刻，序列向上(下)，那么下一刻，它可能反转向下(上)。由于频繁出现逆转，因此反持续性的时间序列具有比随机噪声更剧烈的波动性。

3) 当0＜C≤1，0.5＜H≤1，D＜1.5，是持续性或分形时间序列，波动比较平缓，其变量之间不是相互独立的，而是正相关的。时间序列具有持久性或趋势增强性的特性，如果某一时刻，序列向上(向下)，那么，下一时刻，它可能继续向上(下)。H越接近1，这种持续性越强。这种正相关性，即长期记忆结构，使得时间序列在时间现实方面显示自相似性，即在不同时间尺度上有类似的统计特性。

2.4 Hurst指数判断依据

Hurst指数的取值范围为0＜H≤1，根据H的不同取值，即可区分系统属性，判断时间序列的相关性质^[9-13]。

1) 当0≤H＜0.5，系统是反持久的或遍历性的，经常被称为均值回复的。一个系统如果在前一个期间是向上(下)走的，则在下一个期间多半向下(上)走, 这种反持久性行为的强度依赖于H距离零有多近，距离零越近，C=2^(2H-1)-1中的C就越接近于-0.5，负相关，这种时间序列具有比随机序列更强的突变性或易变性，它是由频繁出现的逆转构成的。

2) 当H=0.5，标志着序列是随机的，事件是随机的和不相关的。

3) 当0.5＜H≤1，是一个持久的或趋势性增强的序列。如果序列在前一个期间是向上(下)走的，则在下一个期间将继续向上(下)走，趋势明显。趋势增强行为的强度或持久性，随H接近于1或C=2^(2H-1)-1的相关性而增加。H越接近0.5，其噪声就越大，趋势也就越不确定。持久性序列是分形布朗运动或有偏随机游动，偏倚的强度依赖于H比0.5大多少。

随机的打乱序列次序计算Hurst指数。

1) 如果序列是一个真正独立的序列，没有长期记忆效应或相关性，那么，Hurst指数性质保持不变，打乱序列次序对于序列的定性方面没有影响。

2) 如果有长期记忆效应，那么，打乱序列次序将破坏系统的结构，计算的Hurst指数H值有明显差别或接近0.5。

3 计算与分析 3.1 计算变量数据

根据式(1)~(6)对2011年1月3日到1898年4月15日全球7.0级以上地震震级数据进行R/S分析，重标极差R/S、样本数n、V等值如表 1。

3.2 计算Hurst指数及分形维D和相关系数C 3.2.1 计算Hurst指数

根据表 1数据绘制log(R/S)与log(n)图，见图 2，得趋势线

$ \log {\left( {\frac{R}{S}} \right)_n} = 0.118 \cdot \log \left( n \right) + 0.328。$

(9)

全球7.0级以上地震震级序列具有以下特征：1)H=0.118＜0.5，是反持久的或遍历性的均值回复过程。2)D=1.882，是分形维，是一个有更多的逆转的系统。3) 相关系数C=2^(2H-1)-1=-0.411，为负相关，且较接近-0.5，是反持续性的时间序列，具有比随机噪声更剧烈的波动性。

3.2.2 打乱序列次序后计算Hurst指数

打乱序列次序的log(R/S)与log(n)图，见图 3，得趋势线

$ \log {\left( {\frac{R}{S}} \right)_n} = 0.050 \cdot \log \left( n \right) + 0.556。$

(10)

由图 3可知：1)Hurst指数H=0.050＜0.5，是反持久的或遍历性的均值回复过程。2)D=1.950，是分形维，是一个有更多的逆转的系统。3) 相关系数C=2^(2H-1)-1=-0.464，为负相关，且较接近-0.5，是反持续性的时间序列，具有比随机噪声更剧烈的波动性。

根据Hurst指数判据，打乱序列次序后的Hurst指数性质保持不变，说明全球7.0级以上地震震级序列是一个真正独立的序列，没有长期记忆效应或相关性，打乱序列次序对于序列的定性方面没有影响。

3.3 V统计量曲线图

根据表 1中V及log(n)数据，做V_n关于log(n)的曲线图如图 4，得趋势线

$ {V_n} = - 0.239\log \left( n \right) + 0.877。$

(11)

曲线是向下倾斜的，全球7.0级以上地震震级序列具有逆状态持续性。

4 结语

应用基于R/S的Hurst指数分析法，探讨了全球7.0级以上地震震级序列的变化趋势，得到了一些有意义的结论，可供地震研究和安全研究工作参考。

1) 地球自转减慢导致的构造应力减弱是地震震级下降的重要原因，地震震级长期情况下有下降趋势，且是长期的极缓慢的过程。

2) 依据R/S分析法计算的Hurst指数可知，全球7.0级以上地震震级时间序列的hurst指数均在(0，0.5)范围之内，表明该序列过程是分形布朗运动，粉红噪声(反持续性)即均值回复过程。H=0.118较接近于0，因此该时间序列具有较强的反持久性。相关系数C=-0.411，为负相关，且较接近于-0.5，该时间序列具有比随机序列更强的突变性或易变性，且由频繁出现的逆转构成。

3) 打乱序列秩序后H=0.050，更接近于0，仍是分形布朗运动状况，粉红噪声(反持续性)即均值回复过程属性不变。V_n关于log(n)的曲线向下倾斜，表明全球7.0级以上地震震级序列具有逆状态持续性，Hurst指数分析得出的分形布朗运动、粉红噪声(反持续性)即均值回复过程是有效的。

通过对全球7.0级以上地震震级序列的长期相关性与持久性的分析表明，发生7.0级以上地震长期情况下有极缓慢下降趋势，而这一趋势又具有极强的反持久性，因此短期内表现出频繁逆转。