Quaternary Activity Characteristic of the Kaifeng Fault
-
摘要: 开封断裂是郑州-开封断裂带的重要组成部分，也是黄河中下游平原的一条重要控制性断裂。利用浅层地震勘探发现，开封断裂为总体走向EW，倾向N的正断层；断裂分为东、西2支，断裂活动性具有明显的分段特征，东支活动性较西支强。钻孔联合剖面揭示开封断裂东支断裂上断点埋深27～35 m，根据地层年代学结果，其最新活动时间为晚更新世。研究成果为开封市城乡规划、重大工程选址和地震区划工作提供了重要的基础数据。Abstract: Kaifeng fault is an important segment of Zhengzhou-Kaifeng fault zone,and it is also an important controlling fault in the middle and lower reaches of the Yellow River Plain. It is found that Kaifeng fault is a normal fault with overall strike EW and dip N by using shallow seismic exploration. The fault can be divided into East and West branches, and has obvious differentiated segmented characteristics. The activity of the East branch is stronger than that of the West Branch. The borehole joint geological profile reveals that the buried depth of the breakpoint on the fault is between 27 m and 35 m. Combining with the chronological results, the latest age of activity is the late Pleistocene. The results can provide important basic data for Kaifeng urban and rural planning, major project site selection and seismic zoning.1) 2河南省平原第四纪地质图及地质图说明书（1/50万）. 郑州：河南省地质局. 开封市活动断层探测与地震危险性评价. 郑州：河南省地震局.
-
引言
地震时间分布模型是进行地震预测和地震危险性分析的重要理论基础，分析海域地震时间分布特征是进行海域地震危险性分析的重要研究内容。地震的发生在时间上符合泊松模型是当前概率地震危险性分析和地震预测中广泛接受的基本假设。泊松模型表示地震发生是相互独立的，即本次地震发生在时间上不受上次地震的影响，也不会对下次地震造成影响。地震时间独立分布模型包含泊松模型、双泊松模型、分段泊松模型、更新模型以及特征地震模型等（Cornell，1968；Gardner等，1974；Schwartz等，1984；胡聿贤，1990；Console等，2003），其中泊松模型在地震学和工程地震学中应用广泛。地震学家们针对地震时间分布是否符合泊松模型做了大量统计检验工作（Gardner等，1974；Bufe等，2005；Michael，2011；Shearer等，2012；Parsons等，2012），目前中国、美国等国家在地震区划和地震危险性分析中均采用泊松模型（高孟潭，1996；潘华等，2013；Petersen等，2014）。
相对于地震的时间独立模型，还有用于描述地震时间分布的时间相依模型，时间相依模型依据Reid（1910）提出的弹性回跳理论。弹性回跳理论认为地震发生是断层应变能渐进且连续的积累过程，当该断层出现突然的应变释放时就会发生构造地震，地震将大大减轻地壳的应变，在稳定的构造力作用下，压力慢慢地重新累积变大，最终发生下次地震。按照弹性回跳理论，断层上本次地震的发生受上次地震的影响，同样对下次地震产生影响，地震是时间相依的。近几十年来，科学家们使用不同地区的地震目录对地震时间相依性进行了实证研究，发现在许多情况下地震时间分布并不符合泊松模型，而时间相依模型能够更好地描述地震的时间分布特征（Utsu，1984；Nishenko等，1987；Ogata，1991；Tripathi，2006；Sharma等，2010；Ellsworth等，2015）。科学家们提出了采用伽马模型（Gamma）、对数正态模型（Lognormal）、威布尔模型（Weibul）以及布朗过程时间模型（Brownian Passage Time，BPT）描述地震时间分布特征（Utsu，1984；Matthews等，2002；Tripathi，2006；Pasari等，2015，2018；Bajaj等，2019），这些模型为地震预测和地震危险性分析提供了重要的理论支撑。
本研究中，以中国海域统一地震目录为基础资料，以泊松模型（指数分布）、伽马模型、对数正态模型、威布尔模型以及布朗过程时间模型为目标模型，回归了各模型参数，并通过AIC、BIC判定准则和K-S检验选择能够描述海域地震活动时间分布特征的最优模型。采用扩散熵分析（Diffusion Entropy Analysis，DEA）和标准差分析（Standard Deviation Analysis，SDA）方法对中国海域地震时间丛集特征和时间相关性进行研究。
1. 数据、模型与方法
1.1 数据
海域地震事件记录是进行海域地震活动特征分析、海域地震危险性分析以及海域地震灾害预测等的重要基础数据，近几十年来，随着地震观测技术的提高和更多地震台网的布设，我国及周边国家、地区均记录了丰富的海域地震事件。地震学家收集了我国海域及周边国家、地区的地震台网记录，经过分析和整理，编制了我国海域统一地震目录，共计61285条地震事件（Xie等，2021），该地震目录为进行我国海域地震危险性分析提供了重要的基础资料。本研究使用的海域地震目录是Xie等（2021）分析整理的中国海域统一地震目录（图1（a））。Xie等（2021）还对海域地震目录做了震级转换和完整性分析等工作，本研究中直接采用其海域地震时间完整性分析结果。使用Gardner等（1974）得出的时空窗法删除了地震目录中的余震事件，该方法被广泛应用于地震区划图编制和地震活动性分析中的余震删除（Shearer等，2012；Daub等，2012；Petersen等，2014）。删除余震后，还剩29450条主震事件，使用删除余震后的主震目录研究地震时间分布特征。
图 1 中国海域及邻区地震分布Figure 1. Earthquake distribution in China sea areas and adjacent areas注：（a）图中红色圆圈为主震，蓝色圆圈为余震；（b）图为近海大陆架各地震带范围，其中，①为华北平原地震带；②为郯庐地震带；③为长江下游-南黄海地震带；④为朝鲜地震带；⑤为东海地震统计区；⑥为华南沿海地震带；⑦为台湾西部地震带；⑧为台湾-马尼拉地震带；⑨为南海地震统计区；⑩为琉球海沟地震带根据新划分的中国海域及领区地震带（图1（b）），以地震带为单元研究了海域地震时间分布特征。海域及领区地震带是进行海域地震活动性参数计算的基本单元，研究海域地震带地震时间分布特征对海域地震活动性参数的确定具有重要科学意义。
1.2 模型
本研究中选择泊松模型（指数分布）、伽马模型、对数正态模型、威布尔模型以及布朗过程时间模型作为参考分布模型，通过统计检验选择能够描述海域地震事件分布特征的最优模型。
（1）指数模型（指数分布）
指数分布是无记忆分布，可由指数分布在极端情况下推导得出，地震学家们通过检验地震时间间隔是否符合指数分布判断地震是否符合指数分布模型。指数分布概率密度函数为：
$$f\left(x \right) = \frac{1}{\mu }{{\rm{e}}^{\frac{{ - x}}{\mu }}}$$ (1) 式中，
$ \mu $ 为均值。（2）伽马（Gamma）模型
对于连续随机变量
$x$ ，具有参数$\alpha $ 和$\lambda $ 的伽马分布概率密度函数为：$$f\left(x \right) = \frac{{{x^{\left({\lambda - 1} \right)}}{{{\rm{e}} }^{\left({{{ - x} / \lambda }} \right)}}}}{{{\lambda ^\alpha }\Gamma \alpha }}$$ (3) 式中，
$\Gamma \alpha = \displaystyle\int_0^\infty {{t^{\alpha - 1}}{{{\rm{e}}}^{ - t}}{\rm{d}}t}$ ，为伽马函数，$\alpha $ 和$\lambda $ 分别为尺度参数和形状参数。（3）对数正态（Log-Normal）模型
在地震学中，对数正态分布通常用于描述大地震时间相依特征。对数正态分布和正态分布关系密切，如果
$x$ 的自然对数服从正态分布，则大小为$n$ 的随机变量$x$ 对数正态分布概率密度函数为：$$f\left(x \right) = \frac{1}{{\sigma x\sqrt {2{\text{π}} } }}{{\rm{e}} ^{\frac{{ - {{\left({\ln \left(x \right) - \mu } \right)}^2}}}{{2{\sigma ^2}}}}}$$ (4) 式中，
$\mu $ 和$\sigma $ 分别为$\ln x$ 均值和标准差。（4）威布尔（Weibull）模型
由Weibull等（1951）提出，Kagan（1997）认为该模型可以用于评估地震的复发概率。Weibull概率密度函数
$f\left(x \right)$ 表达式如下：$$f\left(x \right) = \frac{\beta }{\alpha }{\left({{x / \alpha }} \right)^{\beta - 1}}{{\rm{e}}^{ - {{\left({{x / \alpha }} \right)}^\beta }}}$$ (2) 式中，
$\alpha $ 为尺度参数，$\beta $ 为形状参数，当$\beta $ =1时，威布尔分布退化为指数分布。（5）布朗过程时间（BPT）模型
Matthews等（2002）基于Reid（1910）的弹性回跳理论，结合多个物理参数提出了地震复发的布朗过程时间模型。该模型假设发震断层以恒定的加载速率进行加载，采用标准布朗运动来描述累积加载过程。在实际使用中，模型均值和方差可直接采用特征地震序列计算得到。BPT模型在理论上反映了地震孕育和发生的内在物理机制，是近年来地震学家们广泛使用的模型，概率密度函数为：
$${f_{{\rm{BPT}}}}\left(t \right) = \sqrt {\frac{\mu }{{2{\text{π}} {\alpha ^2}{t^3}}}} {\rm{e}} \left({ - \frac{{{{\left({t - \mu } \right)}^2}}}{{2{\alpha ^2}\mu t}}} \right)$$ (5) 式中，μ为复发间隔均值，α=σ/μ，σ为复发间隔标准差。
布朗过程时间模型累积概率分布函数为：
$$F\left(t \right) = P\left({T \leqslant t} \right) = \int_0^t {{f_{{\rm{BPT}}}}\left(\tau \right)} {\rm{d}}\tau = \varPhi \left[ {{u_1}\left(t \right)} \right] + {{\rm{e}}^{\frac{2}{{{\alpha ^2}}}}}\varPhi \left[ { - {u_2}\left(t \right)} \right]$$ (6) 式中，
${u_1}\left(t \right) = {\alpha ^{ - 1}}\left({{t^{{1 / 2}}}{\mu ^{{{ - 1} / 2}}} - {t^{{{ - 1} / 2}}}{\mu ^{{1 / 2}}}} \right)$ ，${u_2}\left(t \right) = {\alpha ^{ - 1}}\left({{t^{{1 / 2}}}{\mu ^{{{ - 1} / 2}}} + {t^{{{ - 1} / 2}}}{\mu ^{{1 / 2}}}} \right)$ ，$\varPhi \left(t \right) = \dfrac{1}{{\sqrt {2{\text{π}} } }}\displaystyle\int_{ - \infty }^t {{{\rm{e}}^{ - {{{x^2}} / 2}}}{\rm{d}}x}$ ，$\varPhi$ 为标准正态分布函数累积概率函数。1.3 方法
采用扩散熵分析和标准差分析方法研究中国海域地震的时间丛集特征和时间相关性。DEA法是Scafetta等（2004）在计算青少年生育现象时首先采用的，该方法在是否为高斯过程的情况下还是在变异数为无穷大的情况下，都能得到正确的标度参数。此后，该方法被广泛应用于分析不同地区地震的时间分布特征（Jiménez等，2006；Zhou等，2016）。Mega等（2003）将DEA应用于解释1976—2002年加州地区不同地震丛集在时间上的分布关系，得出加州地区大尺度时间上的地震丛集和下个地震丛集具有相关性。Tsai等（2008）用该方法分析了台湾地区地震与地震之间的相关性，发现台湾地区地震时间概率密度函数的标度值为0.83。
DEA方法简述如下：设定单位时间
$\Delta t$ ，将地震序列分成m个单位时间序列ξ i，设定门槛值k，当第i个单位时间内地震个数大于k时，ξ i=1，否则ξ i=0。以t为窗口宽度，移动窗口，每次移动1个单位时间，由n=0到n=m−t开始移动，共生成m−t+1道轨迹，在第n道轨迹，窗口内所有ξ i总和为：$${x_n}\left(t \right) = \int_{n\Delta t}^{n\Delta t + t} {\xi \left({{t^{'}}} \right)} d{t^{'}}$$ (7) 通过上述分析可将由ξ i组成的单一轨迹时间序列变换成不同轨迹在另一变量x的辅助空间和扩散时间t的关系，以
$\varepsilon $ 为间隔（本研究中$\varepsilon $ =1），将x分成多段，计算出现在第j个段落的轨迹数${N_j}\left(t \right)$ ，可得在时间t时落在第j段内轨迹的概率为：$${p_j}\left(t \right) = \frac{{{N_j}\left(t \right)}}{{m - t + 1}}$$ (8) 计算宣农熵（Shannon，1948）：
$$S\left(t \right) = - \int_{ - \infty }^{ + \infty } {{p_j}\left({x,t} \right)\ln \left[ {{p_j}\left({x,t} \right)} \right]} {\rm{d}}x$$ (9) 标度不变性是许多复杂系统的重要性质，并且在许多时间序列的分析中被广泛应用（Mandelbrot，1982）。在对x的概率密度函数进行渐进时间变换后，可得到标度关系（Zinn-Justin，2002）：
$$p\left({x,t} \right) = \frac{1}{{{t^\delta }}}F\left({\frac{x}{{{t^\delta }}}} \right)$$ (10) 式中，
$\delta $ 为概率密度函数标度参数，$F\left({\dfrac{x}{{{t^\delta }}}} \right)$ 为统计分布函数。将式（10）带入式（9）可得：
$$S\left(t \right) = A + \delta \ln \left(t \right)$$ (11) 式中，
$A = - \displaystyle\int_{ - \infty }^{ + \infty } {{\rm{d}}yF\left(y \right)} \ln \left[ {F\left(y \right)} \right]$ ，$\delta $ 在0.5到1之间变化，当$\delta $ 为0.5时，说明时间序列具有正常扩散，事件之间在时间上不具相关性。2. 时间分布模型结果分析
以中国海域地震目录为基本输入，采用极大似然估计方法，回归了上文介绍的各模型参数。为选出最优模型，采用赤池信息准则（Akaike Information Criterion，AIC）、贝叶斯信息准则（Bayesian Information Criterion，BIC）以及K-S检验判断模型拟合优度。AIC和BIC的定义分别为：
$${\rm{AIC}} = 2k - 2\ln \left(L \right)$$ (12) $${\rm{BIC}} = k\ln \left(n \right) - 2\ln \left(L \right)$$ (13) 式中，k为模型参数个数；L为模型似然函数；n为样本量。训练模型时，增加参数数量相当于增加模型复杂度，会增大似然函数，也会导致过拟合现象，针对该问题，AIC和BIC均引入与模型参数个数相关的惩罚项，BIC的惩罚项较AIC的惩罚项大，考虑了样本数量，样本数量过多时，可有效防止模型精度过高造成的模型复杂度过高。从一组可供选择的模型中选择最佳模型时，通常选择AIC或BIC最小的模型。
表1为根据整个中国海域及邻区地震目录回归得到的各模型参数值和模型拟合优度判定准则参数值。对于5.0级以上地震，从表1中可以看出，指数分布、威布尔分布以及伽马分布的AIC、BIC以及K-S检验值相近，其中威布尔分布具有最小的AIC和BIC，指数分布具有最小的K-S检验值（K-S检验值表示地震经验分布与理论模型的差异，该值越小表示越接近），这意味M≥5级地震不拒绝时间独立性的假设，但也具有一定的时间相关性。对于6.0和7.0级以上的地震，根据BIC和K-S检验值，可判定指数分布能够最优地表达海域地震的时间分布特征。由图2的累积分布函数曲线可以看出，指数分布、威布尔分布以及伽马分布与实际数据吻合较好，对数正态分布和布朗过程时间分布与实际数据吻合较差。由上述分析可知，对于整个中国海域及邻区，指数分布能够最佳描述地震的时间分布特征，威布尔分布和伽马分布与实际观测数据吻合良好，但具有2个参数，模型相对复杂。
表 1 中国海域及邻区模型参数及AIC、BIC和K-S检验值Table 1. Model parameters and AIC, BIC and K-S test results for earthquakes in China Sea and adjacent areas震级 模型 模型参数 −lnL AIC BIC K-S检值验 指数分布 μ=23.2760 ― 2770.4772 5542.9545 5547.4587 0.0366 威布尔分布 α=22.1757 β=0.9038 2764.3739 5532.7479 5541.7565 0.0384 1976年以来M≥5 对数正态分布 μ=2.4560 σ=1.5526 2882.3087 5768.6174 5777.6260 0.1067 伽马分布 α=0.8514 λ=27.3383 2764.4759 5532.9518 5541.9604 0.0391 布朗过程时间分布 μ=23.2760 α=0.0329 4548.7329 9101.4659 9110.4745 0.1667 指数分布 μ=188.8710 ― 1422.9627 2847.9253 2851.3547 0.0471 威布尔分布 α=180.4296 β=0.9047 1421.0141 2846.0283 2852.8870 0.0476 1900年以来M≥6 对数正态分布 μ=4.5348 σ=1.7015 1478.6316 2961.2631 2968.1218 0.1144 伽马分布 α=0.8355 λ=226.0487 1420.3845 2844.7689 2851.6276 0.0541 布朗过程时间分布 μ=188.8710 α=0.0874 2152.2080 4308.4161 4315.2748 0.4251 指数分布 μ=920.4309 ― 352.1179 706.2358 708.0424 0.0857 1900年以来M≥7 威布尔分布 α=869.2773 β=0.8697 351.4152 706.8304 710.4438 0.1214 对数正态分布 μ=6.0154 σ=12.3052 372.1317 748.2635 751.8768 0.2007 伽马分布 α=0.7408 λ=1242.5236 350.6158 705.2316 708.8449 0.1290 布朗过程时间分布 μ=920.4309 α=0.0344 545.7260 1095.4521 1099.0654 0.5404 
图 2 中国海域及邻区不同震级地震的时间间隔累积经验分布函数及其对应的5个模型累积分布函数Figure 2. The cumulative empirical distribution function of the time interval of M≥5, 6, 7 earthquakes in the Sea area of China and adjacent areas and the cumulative distribution function of the corresponding five models表2为根据华南沿海地震带地震目录回归得到的各模型参数值和模型拟合优度判定准则参数值。从AIC、BIC和K-S检验值上看，对于M≥4级的地震，指数分布、威布尔分布和伽马分布均表现优良，模型间差异较小，K-S检验接受地震时间间隔分布符合这3种分布的假设；对于M≥5级的地震，指数分布被判定为最优模型，威布尔分布和伽马分布也有良好表现；对于M≥6的地震，K-S检验显示所有模型均不拒绝时间独立性的假设，所有模型的AIC和BIC判定准则值均非常接近，特别是对于M≥4和M≥5级地震表现较差的对数正态分布和布朗过程时间分布，对于M≥6级的地震表现优良，拟合优度好于威布尔分布和伽马分布。这意味着在发震构造特征相似的地震带内，大地震时间分布模型可用对数正态分布和布朗过程时间分布描述，这与国际上的研究结果相似。
表 2 华南沿海地震带模型参数及AIC、BIC和K-S检验值Table 2. Model parameters and AIC, BIC and K-S test results for earthquakes in South China Coastal seismic zone震级 分布模型 模型参数值 −lnL AIC BIC K-S检验值 1970年以来M≥4 指数分布 μ=124.1940 ― 820.8802 1643.7603 1646.7091 0.0524 威布尔分布 α=119.3586 β=0.9161 819.9501 1643.9002 1649.7977 0.0428 对数正态分布 μ=4.1325 σ=11.6783 855.7540 1715.5080 1721.4055 0.1243 伽马分布 α=0.8537 λ=145.4822 819.6564 1643.3129 1649.2104 0.0481 布朗过程时间分布 μ=124.1940 α=0.1551 1205.4927 2414.9854 2420.8829 0.8448 1900年以来M≥5 指数分布 μ=1095.4519 ― 295.9601 593.9202 595.5312 0.1106 威布尔分布 α=1070.9469 β=0.9473 295.8694 595.7388 598.9607 0.1293 对数正态分布 μ=6.3317 σ=11.5744 303.5748 611.1497 614.3715 0.2180 伽马分布 α=0.8788 λ=1246.5401 295.7488 595.4976 598.7195 0.1382 布朗过程时间分布 μ=1095.4519 α=0.20633 327.2432 658.4864 661.7083 0.5460 1900年以来M≥6 指数分布 μ=3647.1044 ― 82.8152 167.6304 167.8276 0.1553 威布尔分布 α=3876.6133 β=1.1812 82.6192 169.2384 169.6329 0.1896 对数正态分布 μ=7.8006 σ=10.9423 82.4707 168.9414 169.3359 0.1443 伽马分布 α=1.3897 λ=2624.4039 82.5429 169.0858 169.4803 0.1870 布朗过程时间分布 μ=3647.1044 α=0.5123.8113 82.5444 169.0887 169.4832 0.1601 表3为根据台湾西部地震带的地震目录回归得到的各模型参数值和模型拟合优度判定准则结果。从AIC、BIC和K-S检验值上看，对于M≥4和M≥5级的地震，指数分布为最优分布，威布尔分布和伽马分布的拟合优度较好；对于M≥6.0级的地震，正态对数分布和布朗过程时间分布拟合优度最好。这再次说明对于大地震，采用对数正态分布、布朗过程时间分布等模型来描述其时间分布特征更合适。
表 3 台湾西部地震带模型参数及AIC、BIC和K-S检验值Table 3. Model parameters and AIC, BIC and K-S test results for earthquakes in the western Taiwan earthquake zone震级 分布模型 模型参数值 −lnL AIC BIC K-S检验值 1970年以来M≥4 指数分布 μ=63.9192 ― 1413.1877 2828.3755 2831.9886 0.0321 威布尔分布 α=62.5341 β=0.9513 1412.6234 2829.2469 2836.4731 0.0413 对数正态分布 μ=3.5179 σ=1.4155 1447.9090 2899.8181 2907.0443 0.0801 伽马分布 α=0.9123 λ=70.0631 1412.4117 2828.8235 2836.0498 0.0396 布朗过程时间分布 μ=63.9192 α=2.0285 1737.7414 3479.4829 3486.7091 0.5304 1900年以来M≥5 指数分布 μ=360.3387 ― 812.6713 1627.3425 1630.1132 0.0691 韦伯分布 α=349.6849 β=0.9394 812.2484 1628.4968 1634.0381 0.0627 对数正态分布 μ=5.2615 σ=1.4189 829.5743 1663.1486 1668.6900 0.1093 伽马分布 α=0.9307 λ=387.1736 812.4684 1628.9368 1634.4782 0.0650 布朗过程时间分布 μ=360.3387 α=2.3725 1047.7444 2099.4887 2105.0301 0.7232 1900年以来M≥6 指数分布 μ=1540.5848 ― 225.1778 452.3555 453.6514 0.1196 韦伯分布 α=1528.7487 β=0.9856 225.1718 454.3435 456.9352 0.1177 对数正态分布 μ=6.8444 σ=0.9483 221.6868 447.3737 449.9653 0.0680 伽马分布 α=1.1471 λ=1343.0297 225.0235 454.0470 456.6387 0.1230 布朗过程时间分布 μ=1540.5848 α=0.4539 221.4897 446.9793 449.5710 0.0680 表4为根据台湾-马尼拉海沟地震带地震目录回归得到的各模型参数值和模型拟合优度判定准则参数值。从AIC、BIC和K-S检验值上看，对于M≥5级的地震，威布尔分布拟合优度最好，伽马分布与威布尔分布较接近；对于M≥6级的地震，指数分布为最优分布，威布尔分布和伽马分布表现较好；对于M≥7级的地震，综合判定指数分布为最优模型，威布尔分布和伽马分布具有较好的拟合优度。从地震事件分布模型上可以看出，在该地震带内5级以上地震具有较强的时间相依性，6级以上地震可接受指数分布假设，但也具有一定的时间相关性。这是由于台湾-马尼拉海沟地震带地震活动非常强烈，经常发生8级左右的地震，大震的发生对稍小震级地震具有一定影响，且震级越小所受影响越大。
表 4 台湾-马尼拉海沟地震带模型参数及AIC、BIC和K-S检验值Table 4. Model parameters and AIC, BIC and K-S test results for earthquakes in Taiwan - Manila trench seismic zone震级 分布模型 模型参数值 −lnL AIC BIC K-S检验值 1900年以来M≥5 指数分布 μ=198.7702 ― 1365.3964 2732.7928 2736.1727 0.1098 威布尔分布 α=172.2580 β=0.7822 1352.9227 2709.8454 2716.6052 0.0247 对数正态分布 μ=4.4145 σ=1.6354 1372.5919 2749.1837 2755.9435 0.0912 伽马分布 α=0.6900 λ=288.0880 1353.9309 2711.8617 2718.6215 0.0348 布朗过程时间分布 μ=198.7702 α=5.1720 1566.5258 3137.0516 3143.8114 0.4935 1900年以来M≥6 指数分布 μ=530.8835 ― 589.2379 1180.4759 1182.8703 0.0686 威布尔分布 α=525.5586 β=0.9762 589.2005 1182.4010 1187.1899 0.0720 对数正态分布 μ=5.6565 σ=1.3571 597.8469 1199.6937 1204.4826 0.1442 伽马分布 α=0.9408 λ=564.3035 589.1379 1182.2757 1187.0646 0.0745 布朗过程时间分布 μ=530.8835 α=1.9082 623.7724 1251.5448 1256.3337 0.3173 1900年以来M≥7 指数分布 μ=2534.7483 ― 132.5677 267.1355 267.8435 0.1465 威布尔分布 α=2598.8333 β=1.0804 132.5075 269.0151 270.4312 0.1351 对数正态分布 μ=7.2616 σ=1.4120 135.4003 274.8006 276.2167 0.1871 伽马分布 α=1.0014 λ=2531.1015 132.5677 269.1355 270.5516 0.1464 布朗过程时间分布 μ=2534.7483 α=0.8039 139.9211 283.8422 285.2583 0.3793 本文计算了琉球海沟地震带各模型参数值和模型拟合优度判定准则参数值，得到了与上述地震带相似的结论。特别地，在琉球海沟地震带，对于M≥7级地震，正态对数分布和布朗过程时间分布拟合优度最好。这再次说明了对数正态分布和布朗过程时间分布在描述大地震时间分布特征的方面具有优越性，同时意味着琉球海沟地震带内7级以上地震具有一定的时间相关性。
根据上述分析可知，对于震级相对较小的地震（M<6），指数分布、伽马分布以及威布尔分布均能较好地描述地震时间分布特征。在大的区域范围内（如整个海域），震级相对较大的地震（M>6）可采用指数分布描述其时间分布特征。在较小的区域范围内（如地震带），大地震时间间隔可能更符合对数正态分布和布朗过程时间分布。
3. 地震时间丛集性和相关性分析
地震在时间上呈现丛集现象是相对于地震在时间上是完全随机分布（指数分布）而言的，是指地震在一段时间内集中成组发生，而后在一段时间内表现出相对频度降低的特征，然后又集中成组发生的现象。地震时间丛集一直是地震学家关心的重要问题，其对地震预测和地震危险性分析具有重要影响。围绕地震在时间上是否丛集，地震学家们做了许多探索和研究，有的研究表明全球地震在时间上遵从指数分布（Michael，2011；Shearer等，2012；Daub等，2012；Parsons等，2012），有的研究认为地震呈现丛集特征，但每个丛集之间相互独立（Mega等，2003），还有的研究表明地震之间存在相互影响且呈现丛集特征（Scafetta等，2004；Bufe，2005；Kulkarni等，2013；Salditch，2020）。
本研究从地震时间间隔的变异系数(Coefficient of Variation)入手，采用扩散熵分析和标准差分析方法对中国海域地震时间丛集特征进行研究。
一套地震时间序列中，若存在丛集性，则在统计量上表现为地震之间的时间间隔离散性变大（相对于泊松模型），变异系数是比较2组数据（本研究为实际观测数据和泊松模型数据）离散程度大小的理想参数，其定义为标准差与均值的比值。对于实际地震记录，连续地震之间的时间间隔是确定的，因此可以计算出1个变异系数。对于指数分布，其变异系数的理论值为1，但由于受样本量的影响，符合泊松模型的随机变量的变异系数不可能完全为1，而是呈现以1为均值的分布。我们采用蒙特卡洛方法模拟1000套和实际观测数据样本相同的符合指数分布的随变量，计算每1套模拟数据的变异系数，将其分布与实际观测数据的变异系数进行比较。图3为整个海域M≥5、M≥6、M≥7级地震变异系数与泊松模型变异系数的对比，可以看出在1%的置信度下，M≥5地震被拒绝接受为指数分布的假设，其变异系数即使在1%置信度的极端情况下也大于指数分布的变异系数，这意味着M≥5地震在时间上具有丛集性。M≥6、M≥7级地震的变异系数在5%置信度下接受为指数分布的假设，这也意味着海域地震中主要是M5～M6级地震体现出时间丛集性。
图 3 海域地震时间间隔变异系数与泊松模型变异系数对比Figure 3. Comparison between the variation coefficient of sea area earthquake time interval and the variation coefficient of Poisson model注：红线为实际地震记录计算的变异系数，直方图为蒙特卡洛模拟1000次泊松过程的变异系数分布采用扩散熵分析和标准差分析方法分析地震的时间相关性。扩散熵分析法的标度参数为δ，标准差分析法的标度参数为H（Hurst exponent）。当随机变量的δ=H=0.5时，表示随机变量符合指数分布，变量之间无相关性，实际情况下若计算的标度值接近0.5±0.05，则认为变量无相关性；当H>δ>0.5时，则认为变量之间具有相关性且具有长期记忆性，即在某一丛集期内，一次地震的发生促使下一次地震的发生；当H<0.5时，则认为地震之间具有负相互作用，即一次地震的发生使下次地震延迟了，在实际情况下，由于地震样本量较少，即使符合指数分布，计算的标度值偏小，这已被蒙特卡洛模拟结果证明，因此若计算的标度值小于0.5，则认为变量符合指数分布。
图4为整个海域地震的扩散熵分析法和标准差分析法分析结果，可以看出对于M≥5级以上地震，标度参数均大于0.5，这意味着海域5级以上地震之间具有相关性；对于M≥6和M≥7地震，标度参数则非常接近0.5，表现出完全随机性，说明地震之间无相关性。这意味海域地区M5～M6地震受更大地震（M≥6、M≥7）活动的影响，大地震发生后促使了中小地震的发生。
图 4 整个海域地震的扩散熵和标准差随时间变化及其标度值Figure 4. Diffusion entropy and standard deviation of earthquakes as the function of time in the whole sea area and the scale values针对不同地震带计算不同震级地震目录的扩散熵分析法和标准差分析法标度参数（表5），对于华南沿海地震带，M≥4级以上地震表现出时间相关性，M≥5、M≥6级以上地震可认为符合指数分布，这说明该地震带M4～M5级地震受M≥5、M≥6级以上地震的影响。在长江下游-南黄海地震带，5级以上地震表现为时间丛集，且受更大地震（M≥6级以上）的控制。对于台湾西部地震带，M≥5、M≥6级地震均具有时间相关性，这说明该地震带内M=5、M=6级地震相互促进或受更大震级地震的控制。在台湾-马尼拉海沟地震带，M≥5地震呈现出时间丛集性，M≥6、M≥7级地震则符合指数分布，这说明了M5～M6级地震受更大地震的影响。在琉球海沟地震带内，M≥7级地震表现出相互促进作用，在地震预测和地震危险性分析中需要重点考虑。
表 5 各地震带不同起始震级地震目录的扩散熵分析法和标准差分析法标度值Table 5. DEA and SDA scale values of earthquake catalogs with different initial magnitudes in each seismic zone地震带 M≥4 M≥5 M≥6 M≥7 δ H δ H δ H δ H 华南沿海地震带 0.56 0.58 0.37 0.46 0.38 0.46 ― ― 长江下游-南黄海地震带 0.44 0.44 0.54 0.58 0.42 0.51 ― ― 台湾西部地震带 0.47 0.47 0.65 0.74 0.56 0.64 ― ― 马尼拉海沟地震带 ― ― 0.59 0.77 0.50 0.57 0.45 0.46 琉球海沟地震带 ― ― ― ― 0.50 0.51 0.56 0.63 以往的研究中，科学家们关注大地震之间的时间相关性（Bufe，2005；Michael，2011；Shearer等，2012；Daub等，2012；Kulkarni等，2013；Salditch，2020），本研究表明，即使在删除余震和前震后，地震目录中的大地震与中小地震之间也具有时间相关性，这意味着大地震的发生促进了中小地震的发生，使地震活动在统计特征上表现出长期记忆性，中小地震在时间轴上呈现丛集特征。
4. 结论与讨论
由上述分析可知，在大的区域范围内（如整个海域区域），起始震级稍小的地震（M≤5、M≤6）时间分布多与指数分布吻合更好，威布尔分布和伽马分布的拟合优度也较好；起始震级大的地震（M≥6、M≥7）则完全符合指数分布。在一个较小区域范围内（如地震带），指数分布、威布尔分布和伽马分布均能较好描述震级较小地震的时间分布特征，对于起始震级大（M≥6、M≥7）的地震，对数正态分布和布朗过程时间分布有时具有较好的拟合优度。这与其他学者的研究结果相似（Matthews等，2002；Bajaj等，2019）。由于威布尔分布和伽马分布在特殊情况下可退化为指数分布，因此，在较大区域内，在遵从科学性的前提下，本着便于使用的原则，泊松模型可选择作为进行地震危险性分析的基本模型。
由于大地震发生率对于地震危险性具有重要的意义，在较小区域内（如某潜在震源区或某条断层），若历史地震目录较为丰富，可尝试采用对数正态模型或布朗过程时间模型计算高震级地震的发生概率和发生频度，这为地震危险性分析中地震发生率的计算提供了新途径，从而使地震危险性计算结果更加科学合理。目前，也有一些学者和机构采用对数正态分布模型或布朗过程时间模型来计算地震危险性（Hebden等，2009；Working Group on California Earthquake Probabilities，2013；Petersen等，2014），并取得了新的认识，为地震灾害的预测预防提供了重要支撑。
对海域地震时间丛集性和相关性进行分析，结果表明，在某区域内稍小震级的地震易受更大地震的影响，在时间上呈现丛集性，而震级大的地震则更多地表现出完全随机性（符合指数分布）。一方面这是由于目前没有有效的方法区分主震和余震，即使删除余震，地震目录中仍然包含数量相当的余震（中小震）。另一方面，在大地震发生后，应力扩散和调整需要较长时间，表现为中小震级地震的发生（Shearer等，2012），在统计特征上表现为丛集性和长期记忆性（即H>0.5），因此，在对某地区进行地震危险性分析时应充分考虑该地区最近一次发生的强震对未来中小地震的影响。本研究还发现在地震活动非常强烈的琉球海沟地震带，M≥7级地震表现出正向的相关性，即某次7级以上地震的发生会促使下次7级以上地震的发生，需要引起特别注意，这一发现对于理解地震孕育发生机理具有一定科学意义。
-
图 3 ZK1和ZK4地震反射时间剖面及其解译成果
Figure 3. Seismic reflection time profile and geological interpretation for survey line ZK1 and ZK4
图 4 ZK2测线地震反射时间剖面及其解译成果
Figure 4. Seismic reflection time profile and geologic interpretation for survey line ZK2
图 5 ZK3测线地震反射时间剖面及其解译成果
Figure 5. Seismic reflection time profile and geologic interpretation for survey line ZK3
表 1 观测系统参数表
Table 1. Parameter table of the observation system
线号 道距/m 炮距/m 偏移距/m 接收道数 覆盖次数 采样间隔/ms 前放增益/db 记录时间/s 震源吨位/t ZK1 5 20 0 240 30 0.5 12 3 26 ZK2 5 20 0 260 32 0.5 12 3 26 ZK3 4 16 0 320 40 0.5 12 2 9 ZK4 5 20 0 260 32 0.5 12 3 26 
下载: 导出CSV
表 2 勘察断点参数表
Table 2. Parameter table of the upper offset point of fault
资料来源 探测断裂 上断点埋深 断距特征 断裂活动性判断 收集石油测线 开封断裂 浅层数据被切除，
可识别最大埋深300 m新近系底部断距150 m PreQ 古近系底界 2 000 m ZK1 开封断裂西支 95 m 3～5 m $ {\mathrm{Q}}_{\mathrm{P}}^{2} $ ZK4 开封断裂西支 139 m 5～8 m $ {\mathrm{Q}}_{\mathrm{P}}^{2} $ ZK2 开封断裂东支 50 m 2～3 m $ {\mathrm{Q}}_{\mathrm{P}}^{3} $ ZK3 开封断裂东支 45～50 m 1～3 m $ {\mathrm{Q}}_{\mathrm{P}}^{3} $
下载: 导出CSV
表 3 钻孔联合剖面主要地层层断距
Table 3. Respective displacement of the mark layer in composite geological profile
主要标志层 地层岩性 下盘埋深/m 上盘埋深/m 断距/m 地层①泥炭层 灰黑色泥碳层 Z2:15.29～15.95
Z4:15.45～15.74
Z8:15.81～16.21
Z9:15.67～16.42Z1:15.37～16.55
Z7:15.49～15.81
Z5:15.97～16.04
Z10:15.47～15.81— 地层②底界面 上部为棕灰、棕黄色细砂，
下部棕灰、棕黄色黏土-粉土，
界线明显Z1:35.29
Z6:35.53
Z8:35.84
Z9:35.74
Z11:35.23
Z12:35.84Z1:36.36
Z3:36.78
Z5:36.55
Z7:37.42
Z10:36.990.75 地层③泥炭夹层 黑色泥炭层 Z2:40.84～41.47
Z4:40.20～40.77
Z6:39.92～40.37
Z8:40.93～42.35
Z9:40.14～42.51
Z11:40.19～40.34Z5:42.23～42.38
Z7:42.03～42.24
Z10:42.96～43.261 地层③底界面 上部为黄棕色、灰棕色粉砂，下部为深棕色、黄棕色粉质黏土，界线清晰 Z2:50.20
Z3:50.81
Z4:50.26
Z11:50.74Z7:52.24
Z5:51.80
Z10:51.971.7 地层④底界面 上部深棕色粉砂，下部为深棕色黏土分界面，界线清晰 Z2:60.37
Z4:60.51
Z9:62.48Z5:65.58
Z7:65.61
Z1:66.443.1 地层⑤底界面 上部为棕黄色、浅棕色粉砂
下部深棕色黏土层，界线清晰Z2:77.48
Z3:79.20
Z4:78.27Z1:88.64
Z5:86.987.1
下载: 导出CSV
-
[1] 柴炽章, 孟广魁, 杜鹏等, 2006. 隐伏活动断层的多层次综合探测——以银川隐伏活动断层为例. 地震地质, 28（4）: 536—546 doi: 10.3969/j.issn.0253-4967.2006.04.002Chai C. Z. , Meng G. K. , Du P. , et al. , 2006. Comprehensive multi-level exploration of buried active fault: an example of Yinchuan buried active fault. Seismology and Geology, 28(4): 536—546. (in Chinese) doi: 10.3969/j.issn.0253-4967.2006.04.002 [2] 邓起东, 徐锡伟, 张先康等, 2003. 城市活动断裂探测的方法和技术. 地学前缘, 10（1）: 93—104 doi: 10.3321/j.issn:1005-2321.2003.01.012Deng Q. D. , Xu X. W. , Zhang X. K. , et al. , 2003. Methods and techniques for surveying and prospecting active faults in urban areas. Earth Science Frontiers, 10(1): 93—104. (in Chinese) doi: 10.3321/j.issn:1005-2321.2003.01.012 [3] 侯江飞, 邢磊, 张扬等, 2021. 新乡—商丘断裂延津段浅部地层结构特征研究. 工程地球物理学报, 18（4）: 486—494 doi: 10.3969/j.issn.1672-7940.2021.04.0011Hou J. F. , Xing L. , Zhang Y. , et al. , 2021. The shallow structural characteristics of the Yanjin section of Xinxiang-Shangqiu fault. Chinese Journal of Engineering Geophysics, 18(4): 486—494. (in Chinese) doi: 10.3969/j.issn.1672-7940.2021.04.0011 [4] 雷启云, 柴炽章, 郑文俊等, 2014. 钻探揭示的黄河断裂北段活动性和滑动速率. 地震地质, 36（2）: 464—477 doi: 10.3969/j.issn.0253-4967.2014.02.015Lei Q. Y. , Chai C. Z. , Zheng W. J. , et al. , 2014. Activity and slip rate of the northern section of Yellow River fault revealed by drilling. Seismology and Geology, 36(2): 464—477. (in Chinese) doi: 10.3969/j.issn.0253-4967.2014.02.015 [5] 李峰, 环文林, 2021. 粤东滨海断裂带第四纪活动特征研究. 震灾防御技术, 16（1）: 19—28 doi: 10.11899/zzfy20210103Li F. , Huan W. L. , 2021. Quaternary activity characteristics of the coastal fault zone in eastern Guangdong. Technology for Earthquake Disaster Prevention, 16(1): 19—28. (in Chinese) doi: 10.11899/zzfy20210103 [6] 刘保金, 柴炽章, 酆少英等, 2008. 第四纪沉积区断层及其上断点探测的地震方法技术——以银川隐伏活动断层为例. 地球物理学报, 51（5）: 1475—1483 doi: 10.3321/j.issn:0001-5733.2008.05.021Liu B. J. , Chai C. Z. , Feng S. Y. , et al. , 2008. Seismic exploration method for buried fault and its up-breakpoint in Quaternary sediment area-An example of Yinchuan buried active fault. Chinese Journal of Geophysics, 51(5): 1475—1483. (in Chinese) doi: 10.3321/j.issn:0001-5733.2008.05.021 [7] 宋佳佳, 田晓峰, 王帅军等, 2021. 南华北块体及邻区地质结构构造研究进展. 地球物理学进展, 36（3）: 940—952 doi: 10.6038/pg2021EE0286Song J. J. , Tian X. F. , Wang S. J. , et al. , 2021. Research progress of geological structure in the South North China and adjacent areas. Progress in Geophysics, 36(3): 940—952. (in Chinese) doi: 10.6038/pg2021EE0286 [8] 王定一, 汤锡元, 陈乃明, 1994. 开封坳陷构造特征、形成演化与油气远景. 石油学报, 15（2）: 39—47 doi: 10.3321/j.issn:0253-2697.1994.02.005Wang D. Y. , Tang X. Y. , Chen N. M. , 1994. Structural features, evolution and oil and gas prospects of Kaifeng depression. Acta Petrolei Sinica, 15(2): 39—47. (in Chinese) doi: 10.3321/j.issn:0253-2697.1994.02.005 [9] 王晶, 李德文, 马保起等, 2021. 黄河郑州—济南段河型变化特征及其与隐伏断层活动的关系. 地质论评, 67（6）: 1605—1618Wang J. , Li D. W. , Ma B. Q. , et al. , 2021. River pattern anomalies of Zhengzhou—Jinan segment of Yellow River, North of China, and its relationship with activity of subtle faults. Geological Review, 67(6): 1605—1618. (in Chinese) [10] 王学潮, 向宏发, 2001. 聊城—兰考断裂综合研究及黄河下游河道稳定性分析. 郑州: 黄河水利出版社, 1—163. [11] 王志铄, 马兴全, 2018. 郑州-开封断裂新生代活动特征. 地震地质, 40（3）: 511—522 doi: 10.3969/j.issn.0253-4967.2018.03.001Wang Z. S. , Ma X. Q. , 2018. The activity characteristics of Zhengzhou-Kaifeng fault during kainozoic. Seismology and Geology, 40(3): 511—522. (in Chinese) doi: 10.3969/j.issn.0253-4967.2018.03.001 [12] 闫震鹏, 2009. 地质构造对黄河下游悬河稳定性的控制作用. 人民黄河, 31（7）: 16—17, 19 doi: 10.3969/j.issn.1000-1379.2009.07.007Yan Z. P. , 2009. Control function of Geologic structure to the stability of elevated section of the Lower Yellow River. Yellow River, 31(7): 16—17, 19. (in Chinese) doi: 10.3969/j.issn.1000-1379.2009.07.007 [13] Liang K. , Sun C. B. , Ma B. Q. , et al. , 2018. Investigation of the Yellow River buried fault in the Wuhai basin, northwestern ordos block, China, using deep/shallow seismic reflection and drilling techniques. Journal of Asian Earth Sciences, 163: 54—69. doi: 10.1016/j.jseaes.2018.05.025 [14] Yu Z. Y. , Pan H. , Xi H. , et al. , 2019. Late Quaternary paleoseismicity of the Xiadian fault in the North China Plain with implications for earthquake potential. Journal of Asian Earth Sciences, 184: 103997. doi: 10.1016/j.jseaes.2019.103997 期刊类型引用(1)
1. 陆东华，陈鲲，魏佳彤. 台湾地区浅层地震活动性模型初步研究. 震灾防御技术. 2024(03): 494-503 . 
本站查看其他类型引用(2)
-
下载:
点击查看大图
计量
- 文章访问数: 261
- HTML全文浏览量: 66
- PDF下载量: 63
- 被引次数: 3
