• ISSN 1673-5722
  • CN 11-5429/P

开封断裂第四纪活动特征

孙杰 张予川 徐婉君 李涛 张扬

徐伟进,李雪婧,谢卓娟,吕悦军,高战武,2021. 中国海域及邻区地震时间分布特征研究. 震灾防御技术,16(1):39−50. doi:10.11899/zzfy20210105. doi: 10.11899/zzfy20210105
引用本文: 孙杰,张予川,徐婉君,李涛,张扬,2022. 开封断裂第四纪活动特征. 震灾防御技术,17(2):234−241. doi:10.11899/zzfy20220204. doi: 10.11899/zzfy20220204
Xu Weijin, Li Xuejing, Xie Zhuojuan, Lv Yuejun, Gao Zhanwu. Temporal Distribution Characteristics of Earthquakes in the China Sea and Adjacent Areas[J]. Technology for Earthquake Disaster Prevention, 2021, 16(1): 39-50. doi: 10.11899/zzfy20210105
Citation: Sun Jie, Zhang Yuchuan, Xu Wanjun, Li Tao, Zhang Yang. Quaternary Activity Characteristic of the Kaifeng Fault[J]. Technology for Earthquake Disaster Prevention, 2022, 17(2): 234-241. doi: 10.11899/zzfy20220204

开封断裂第四纪活动特征

doi: 10.11899/zzfy20220204
基金项目: 中央公益性科研院所基本科研业务专项(IGCEA1902);中国地震局地震科技星火计划项目(XH20037);开封市活动断层探测与地震危险性评价服务项目;河南省地震构造探查工程项目
详细信息
    作者简介:

    孙杰,男,生于1984年。硕士,高级工程师。主要从事地震地质相关研究。 E-mail: 33837077@qq.com

  • 2河南省平原第四纪地质图及地质图说明书(1/50万). 郑州:河南省地质局. 开封市活动断层探测与地震危险性评价. 郑州:河南省地震局.

Quaternary Activity Characteristic of the Kaifeng Fault

  • 摘要: 开封断裂是郑州-开封断裂带的重要组成部分,也是黄河中下游平原的一条重要控制性断裂。利用浅层地震勘探发现,开封断裂为总体走向EW,倾向N的正断层;断裂分为东、西2支,断裂活动性具有明显的分段特征,东支活动性较西支强。钻孔联合剖面揭示开封断裂东支断裂上断点埋深27~35 m,根据地层年代学结果,其最新活动时间为晚更新世。研究成果为开封市城乡规划、重大工程选址和地震区划工作提供了重要的基础数据。
    1)  2河南省平原第四纪地质图及地质图说明书(1/50万). 郑州:河南省地质局. 开封市活动断层探测与地震危险性评价. 郑州:河南省地震局.
  • 地震时间分布模型是进行地震预测和地震危险性分析的重要理论基础,分析海域地震时间分布特征是进行海域地震危险性分析的重要研究内容。地震的发生在时间上符合泊松模型是当前概率地震危险性分析和地震预测中广泛接受的基本假设。泊松模型表示地震发生是相互独立的,即本次地震发生在时间上不受上次地震的影响,也不会对下次地震造成影响。地震时间独立分布模型包含泊松模型、双泊松模型、分段泊松模型、更新模型以及特征地震模型等(Cornell,1968Gardner等,1974Schwartz等,1984胡聿贤,1990Console等,2003),其中泊松模型在地震学和工程地震学中应用广泛。地震学家们针对地震时间分布是否符合泊松模型做了大量统计检验工作(Gardner等,1974Bufe等,2005Michael,2011Shearer等,2012Parsons等,2012),目前中国、美国等国家在地震区划和地震危险性分析中均采用泊松模型(高孟潭,1996潘华等,2013Petersen等,2014)。

    相对于地震的时间独立模型,还有用于描述地震时间分布的时间相依模型,时间相依模型依据Reid(1910)提出的弹性回跳理论。弹性回跳理论认为地震发生是断层应变能渐进且连续的积累过程,当该断层出现突然的应变释放时就会发生构造地震,地震将大大减轻地壳的应变,在稳定的构造力作用下,压力慢慢地重新累积变大,最终发生下次地震。按照弹性回跳理论,断层上本次地震的发生受上次地震的影响,同样对下次地震产生影响,地震是时间相依的。近几十年来,科学家们使用不同地区的地震目录对地震时间相依性进行了实证研究,发现在许多情况下地震时间分布并不符合泊松模型,而时间相依模型能够更好地描述地震的时间分布特征(Utsu,1984Nishenko等,1987Ogata,1991Tripathi,2006Sharma等,2010Ellsworth等,2015)。科学家们提出了采用伽马模型(Gamma)、对数正态模型(Lognormal)、威布尔模型(Weibul)以及布朗过程时间模型(Brownian Passage Time,BPT)描述地震时间分布特征(Utsu,1984Matthews等,2002Tripathi,2006Pasari等,20152018Bajaj等,2019),这些模型为地震预测和地震危险性分析提供了重要的理论支撑。

    本研究中,以中国海域统一地震目录为基础资料,以泊松模型(指数分布)、伽马模型、对数正态模型、威布尔模型以及布朗过程时间模型为目标模型,回归了各模型参数,并通过AIC、BIC判定准则和K-S检验选择能够描述海域地震活动时间分布特征的最优模型。采用扩散熵分析(Diffusion Entropy Analysis,DEA)和标准差分析(Standard Deviation Analysis,SDA)方法对中国海域地震时间丛集特征和时间相关性进行研究。

    海域地震事件记录是进行海域地震活动特征分析、海域地震危险性分析以及海域地震灾害预测等的重要基础数据,近几十年来,随着地震观测技术的提高和更多地震台网的布设,我国及周边国家、地区均记录了丰富的海域地震事件。地震学家收集了我国海域及周边国家、地区的地震台网记录,经过分析和整理,编制了我国海域统一地震目录,共计61285条地震事件(Xie等,2021),该地震目录为进行我国海域地震危险性分析提供了重要的基础资料。本研究使用的海域地震目录是Xie等(2021)分析整理的中国海域统一地震目录(图1(a))。Xie等(2021)还对海域地震目录做了震级转换和完整性分析等工作,本研究中直接采用其海域地震时间完整性分析结果。使用Gardner等(1974)得出的时空窗法删除了地震目录中的余震事件,该方法被广泛应用于地震区划图编制和地震活动性分析中的余震删除(Shearer等,2012Daub等,2012Petersen等,2014)。删除余震后,还剩29450条主震事件,使用删除余震后的主震目录研究地震时间分布特征。

    图 1  中国海域及邻区地震分布
    Figure 1.  Earthquake distribution in China sea areas and adjacent areas
    :(a)图中红色圆圈为主震,蓝色圆圈为余震;(b)图为近海大陆架各地震带范围,其中,①为华北平原地震带;②为郯庐地震带;③为长江下游-南黄海地震带;④为朝鲜地震带;⑤为东海地震统计区;⑥为华南沿海地震带;⑦为台湾西部地震带;⑧为台湾-马尼拉地震带;⑨为南海地震统计区;⑩为琉球海沟地震带

    根据新划分的中国海域及领区地震带(图1(b)),以地震带为单元研究了海域地震时间分布特征。海域及领区地震带是进行海域地震活动性参数计算的基本单元,研究海域地震带地震时间分布特征对海域地震活动性参数的确定具有重要科学意义。

    本研究中选择泊松模型(指数分布)、伽马模型、对数正态模型、威布尔模型以及布朗过程时间模型作为参考分布模型,通过统计检验选择能够描述海域地震事件分布特征的最优模型。

    (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$为标准正态分布函数累积概率函数。

    采用扩散熵分析和标准差分析方法研究中国海域地震的时间丛集特征和时间相关性。DEA法是Scafetta等(2004)在计算青少年生育现象时首先采用的,该方法在是否为高斯过程的情况下还是在变异数为无穷大的情况下,都能得到正确的标度参数。此后,该方法被广泛应用于分析不同地区地震的时间分布特征(Jiménez等,2006Zhou等,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=mt开始移动,共生成mt+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时,说明时间序列具有正常扩散,事件之间在时间上不具相关性。

    以中国海域地震目录为基本输入,采用极大似然估计方法,回归了上文介绍的各模型参数。为选出最优模型,采用赤池信息准则(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
    震级模型模型参数−lnLAICBICK-S检值验
    指数分布μ=23.27602770.47725542.95455547.45870.0366
    威布尔分布α=22.1757β=0.90382764.37395532.74795541.75650.0384
    1976年以来M≥5对数正态分布μ=2.4560σ=1.55262882.30875768.61745777.62600.1067
    伽马分布α=0.8514λ=27.33832764.47595532.95185541.96040.0391
    布朗过程时间分布μ=23.2760α=0.03294548.73299101.46599110.47450.1667
    指数分布μ=188.87101422.96272847.92532851.35470.0471
    威布尔分布α=180.4296β=0.90471421.01412846.02832852.88700.0476
    1900年以来M≥6对数正态分布μ=4.5348σ=1.70151478.63162961.26312968.12180.1144
    伽马分布α=0.8355λ=226.04871420.38452844.76892851.62760.0541
    布朗过程时间分布μ=188.8710α=0.08742152.20804308.41614315.27480.4251
    指数分布μ=920.4309352.1179706.2358708.04240.0857
    1900年以来M≥7威布尔分布α=869.2773β=0.8697351.4152706.8304710.44380.1214
    对数正态分布μ=6.0154σ=12.3052372.1317748.2635751.87680.2007
    伽马分布α=0.7408λ=1242.5236350.6158705.2316708.84490.1290
    布朗过程时间分布μ=920.4309α=0.0344545.72601095.45211099.06540.5404
    下载: 导出CSV 
    | 显示表格
    图 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
    震级分布模型模型参数值−lnLAICBICK-S检验值
    1970年以来M≥4指数分布μ=124.1940820.88021643.76031646.70910.0524
    威布尔分布α=119.3586β=0.9161819.95011643.90021649.79770.0428
    对数正态分布μ=4.1325σ=11.6783855.75401715.50801721.40550.1243
    伽马分布α=0.8537λ=145.4822819.65641643.31291649.21040.0481
    布朗过程时间分布μ=124.1940α=0.15511205.49272414.98542420.88290.8448
    1900年以来M≥5指数分布μ=1095.4519295.9601593.9202595.53120.1106
    威布尔分布α=1070.9469β=0.9473295.8694595.7388598.96070.1293
    对数正态分布μ=6.3317σ=11.5744303.5748611.1497614.37150.2180
    伽马分布α=0.8788λ=1246.5401295.7488595.4976598.71950.1382
    布朗过程时间分布μ=1095.4519α=0.20633327.2432658.4864661.70830.5460
    1900年以来M≥6指数分布μ=3647.104482.8152167.6304167.82760.1553
    威布尔分布α=3876.6133β=1.181282.6192169.2384169.63290.1896
    对数正态分布μ=7.8006σ=10.942382.4707168.9414169.33590.1443
    伽马分布α=1.3897λ=2624.403982.5429169.0858169.48030.1870
    布朗过程时间分布μ=3647.1044α=0.5123.811382.5444169.0887169.48320.1601
    下载: 导出CSV 
    | 显示表格

    表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
    震级分布模型模型参数值−lnLAICBICK-S检验值
    1970年以来M≥4指数分布μ=63.91921413.18772828.37552831.98860.0321
    威布尔分布α=62.5341β=0.95131412.62342829.24692836.47310.0413
    对数正态分布μ=3.5179σ=1.41551447.90902899.81812907.04430.0801
    伽马分布α=0.9123λ=70.06311412.41172828.82352836.04980.0396
    布朗过程时间分布μ=63.9192α=2.02851737.74143479.48293486.70910.5304
    1900年以来M≥5指数分布μ=360.3387812.67131627.34251630.11320.0691
    韦伯分布α=349.6849β=0.9394812.24841628.49681634.03810.0627
    对数正态分布μ=5.2615σ=1.4189829.57431663.14861668.69000.1093
    伽马分布α=0.9307λ=387.1736812.46841628.93681634.47820.0650
    布朗过程时间分布μ=360.3387α=2.37251047.74442099.48872105.03010.7232
    1900年以来M≥6指数分布μ=1540.5848225.1778452.3555453.65140.1196
    韦伯分布α=1528.7487β=0.9856225.1718454.3435456.93520.1177
    对数正态分布μ=6.8444σ=0.9483221.6868447.3737449.96530.0680
    伽马分布α=1.1471λ=1343.0297225.0235454.0470456.63870.1230
    布朗过程时间分布μ=1540.5848α=0.4539221.4897446.9793449.57100.0680
    下载: 导出CSV 
    | 显示表格

    表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
    震级分布模型模型参数值−lnLAICBICK-S检验值
    1900年以来M≥5指数分布μ=198.77021365.39642732.79282736.17270.1098
    威布尔分布α=172.2580β=0.78221352.92272709.84542716.60520.0247
    对数正态分布μ=4.4145σ=1.63541372.59192749.18372755.94350.0912
    伽马分布α=0.6900λ=288.08801353.93092711.86172718.62150.0348
    布朗过程时间分布μ=198.7702α=5.17201566.52583137.05163143.81140.4935
    1900年以来M≥6指数分布μ=530.8835589.23791180.47591182.87030.0686
    威布尔分布α=525.5586β=0.9762589.20051182.40101187.18990.0720
    对数正态分布μ=5.6565σ=1.3571597.84691199.69371204.48260.1442
    伽马分布α=0.9408λ=564.3035589.13791182.27571187.06460.0745
    布朗过程时间分布μ=530.8835α=1.9082623.77241251.54481256.33370.3173
    1900年以来M≥7指数分布μ=2534.7483132.5677267.1355267.84350.1465
    威布尔分布α=2598.8333β=1.0804132.5075269.0151270.43120.1351
    对数正态分布μ=7.2616σ=1.4120135.4003274.8006276.21670.1871
    伽马分布α=1.0014λ=2531.1015132.5677269.1355270.55160.1464
    布朗过程时间分布μ=2534.7483α=0.8039139.9211283.8422285.25830.3793
    下载: 导出CSV 
    | 显示表格

    本文计算了琉球海沟地震带各模型参数值和模型拟合优度判定准则参数值,得到了与上述地震带相似的结论。特别地,在琉球海沟地震带,对于M≥7级地震,正态对数分布和布朗过程时间分布拟合优度最好。这再次说明了对数正态分布和布朗过程时间分布在描述大地震时间分布特征的方面具有优越性,同时意味着琉球海沟地震带内7级以上地震具有一定的时间相关性。

    根据上述分析可知,对于震级相对较小的地震(M<6),指数分布、伽马分布以及威布尔分布均能较好地描述地震时间分布特征。在大的区域范围内(如整个海域),震级相对较大的地震(M>6)可采用指数分布描述其时间分布特征。在较小的区域范围内(如地震带),大地震时间间隔可能更符合对数正态分布和布朗过程时间分布。

    地震在时间上呈现丛集现象是相对于地震在时间上是完全随机分布(指数分布)而言的,是指地震在一段时间内集中成组发生,而后在一段时间内表现出相对频度降低的特征,然后又集中成组发生的现象。地震时间丛集一直是地震学家关心的重要问题,其对地震预测和地震危险性分析具有重要影响。围绕地震在时间上是否丛集,地震学家们做了许多探索和研究,有的研究表明全球地震在时间上遵从指数分布(Michael,2011Shearer等,2012Daub等,2012Parsons等,2012),有的研究认为地震呈现丛集特征,但每个丛集之间相互独立(Mega等,2003),还有的研究表明地震之间存在相互影响且呈现丛集特征(Scafetta等,2004Bufe,2005Kulkarni等,2013Salditch,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≥4M≥5M≥6M≥7
    δHδHδHδH
    华南沿海地震带0.560.580.370.460.380.46
    长江下游-南黄海地震带0.440.440.540.580.420.51
    台湾西部地震带0.470.470.650.740.560.64
    马尼拉海沟地震带0.590.770.500.570.450.46
    琉球海沟地震带0.500.510.560.63
    下载: 导出CSV 
    | 显示表格

    以往的研究中,科学家们关注大地震之间的时间相关性(Bufe,2005Michael,2011Shearer等,2012Daub等,2012Kulkarni等,2013Salditch,2020),本研究表明,即使在删除余震和前震后,地震目录中的大地震与中小地震之间也具有时间相关性,这意味着大地震的发生促进了中小地震的发生,使地震活动在统计特征上表现出长期记忆性,中小地震在时间轴上呈现丛集特征。

    由上述分析可知,在大的区域范围内(如整个海域区域),起始震级稍小的地震(M≤5、M≤6)时间分布多与指数分布吻合更好,威布尔分布和伽马分布的拟合优度也较好;起始震级大的地震(M≥6、M≥7)则完全符合指数分布。在一个较小区域范围内(如地震带),指数分布、威布尔分布和伽马分布均能较好描述震级较小地震的时间分布特征,对于起始震级大(M≥6、M≥7)的地震,对数正态分布和布朗过程时间分布有时具有较好的拟合优度。这与其他学者的研究结果相似(Matthews等,2002Bajaj等,2019)。由于威布尔分布和伽马分布在特殊情况下可退化为指数分布,因此,在较大区域内,在遵从科学性的前提下,本着便于使用的原则,泊松模型可选择作为进行地震危险性分析的基本模型。

    由于大地震发生率对于地震危险性具有重要的意义,在较小区域内(如某潜在震源区或某条断层),若历史地震目录较为丰富,可尝试采用对数正态模型或布朗过程时间模型计算高震级地震的发生概率和发生频度,这为地震危险性分析中地震发生率的计算提供了新途径,从而使地震危险性计算结果更加科学合理。目前,也有一些学者和机构采用对数正态分布模型或布朗过程时间模型来计算地震危险性(Hebden等,2009Working Group on California Earthquake Probabilities,2013Petersen等,2014),并取得了新的认识,为地震灾害的预测预防提供了重要支撑。

    对海域地震时间丛集性和相关性进行分析,结果表明,在某区域内稍小震级的地震易受更大地震的影响,在时间上呈现丛集性,而震级大的地震则更多地表现出完全随机性(符合指数分布)。一方面这是由于目前没有有效的方法区分主震和余震,即使删除余震,地震目录中仍然包含数量相当的余震(中小震)。另一方面,在大地震发生后,应力扩散和调整需要较长时间,表现为中小震级地震的发生(Shearer等,2012),在统计特征上表现为丛集性和长期记忆性(即H>0.5),因此,在对某地区进行地震危险性分析时应充分考虑该地区最近一次发生的强震对未来中小地震的影响。本研究还发现在地震活动非常强烈的琉球海沟地震带,M≥7级地震表现出正向的相关性,即某次7级以上地震的发生会促使下次7级以上地震的发生,需要引起特别注意,这一发现对于理解地震孕育发生机理具有一定科学意义。

  • 图  1  区域地质构造图

    Figure  1.  Regional tectonic map

    图  2  地震勘探测线位置图

    Figure  2.  Location map of seismic exploration line

    图  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

    图  6  LH2钻孔联合探测地质剖面

    Figure  6.  Composite geological profile of LH2

    表  1  观测系统参数表

    Table  1.   Parameter table of the observation system

    线号道距/m炮距/m偏移距/m接收道数覆盖次数采样间隔/ms前放增益/db记录时间/s震源吨位/t
    ZK15200240300.512326
    ZK25200260320.512326
    ZK34160320400.51229
    ZK45200260320.512326
    下载: 导出CSV

    表  2  勘察断点参数表

    Table  2.   Parameter table of the upper offset point of fault

    资料来源探测断裂上断点埋深断距特征断裂活动性判断
    收集石油测线开封断裂浅层数据被切除,
    可识别最大埋深300 m
    新近系底部断距150 mPreQ
    古近系底界 2 000 m
    ZK1开封断裂西支95 m3~5 m$ {\mathrm{Q}}_{\mathrm{P}}^{2} $
    ZK4开封断裂西支139 m5~8 m$ {\mathrm{Q}}_{\mathrm{P}}^{2} $
    ZK2开封断裂东支50 m2~3 m$ {\mathrm{Q}}_{\mathrm{P}}^{3} $
    ZK3开封断裂东支45~50 m1~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.42
    Z1: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.84
    Z1:36.36
    Z3:36.78
    Z5:36.55
    Z7:37.42
    Z10:36.99

    0.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.34
    Z5:42.23~42.38
    Z7:42.03~42.24
    Z10:42.96~43.26
    1
    地层③底界面 上部为黄棕色、灰棕色粉砂,下部为深棕色、黄棕色粉质黏土,界线清晰 Z2:50.20
    Z3:50.81
    Z4:50.26
    Z11:50.74
    Z7:52.24
    Z5:51.80
    Z10:51.97
    1.7
    地层④底界面 上部深棕色粉砂,下部为深棕色黏土分界面,界线清晰 Z2:60.37
    Z4:60.51
    Z9:62.48
    Z5:65.58
    Z7:65.61
    Z1:66.44
    3.1
    地层⑤底界面 上部为棕黄色、浅棕色粉砂
    下部深棕色黏土层,界线清晰
    Z2:77.48
    Z3:79.20
    Z4:78.27
    Z1:88.64
    Z5:86.98
    7.1
    下载: 导出CSV
  • [1] 柴炽章, 孟广魁, 杜鹏等, 2006. 隐伏活动断层的多层次综合探测——以银川隐伏活动断层为例. 地震地质, 28(4): 536—546 doi: 10.3969/j.issn.0253-4967.2006.04.002

    Chai 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.012

    Deng 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.0011

    Hou 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.015

    Lei 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/zzfy20210103

    Li 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.021

    Liu 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/pg2021EE0286

    Song 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.005

    Wang 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—1618

    Wang 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.001

    Wang 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.007

    Yan 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)

  • 加载中
图(6) / 表(3)
计量
  • 文章访问数:  261
  • HTML全文浏览量:  66
  • PDF下载量:  63
  • 被引次数: 3
出版历程
  • 收稿日期:  2022-02-22
  • 刊出日期:  2022-06-30

目录

/

返回文章
返回