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

基于离散单元法碎石垫层隔震性能的宏细观分析

张昕 戴国亮 栾阳 张倩 贾其军

徐伟进,李雪婧,谢卓娟,吕悦军,高战武,2021. 中国海域及邻区地震时间分布特征研究. 震灾防御技术,16(1):39−50. doi:10.11899/zzfy20210105. doi: 10.11899/zzfy20210105
引用本文: 张昕,戴国亮,栾阳,张倩,贾其军,2023. 基于离散单元法碎石垫层隔震性能的宏细观分析. 震灾防御技术,18(3):595−603. doi:10.11899/zzfy20230316. doi: 10.11899/zzfy20230316
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: Zhang Xin, Dai Guoliang, Luan Yang, Zhang Qian, Jia Qijun. Macro and Micro Analysis of Seismic Isolation Performance of Gravel Cushion Based on Discrete Element Method[J]. Technology for Earthquake Disaster Prevention, 2023, 18(3): 595-603. doi: 10.11899/zzfy20230316

基于离散单元法碎石垫层隔震性能的宏细观分析

doi: 10.11899/zzfy20230316
基金项目: 中交集团重点科研课题(2016-zjjt-24)
详细信息
    作者简介:

    张昕,男,生于1975年。高级工程师。主要从事工程项目管理及岩土工程抗震方面研究工作。E-mail:zhangx@crbc.com

    通讯作者:

    贾其军,男,生于1976年。博士,高级工程师。主要从事岩土工程研究工作。E-mail:18513671198@163.com

Macro and Micro Analysis of Seismic Isolation Performance of Gravel Cushion Based on Discrete Element Method

  • 摘要: 为研究非接触桩箱复合基础中碎石垫层隔震性能,采用颗粒流软件PFC3D对碎石垫层和沉箱进行模拟分析,从宏细观多角度对碎石垫层隔震性能进行分析。研究结果表明,竖向压力、垫层厚度对垫层隔震效果的影响较大;地震动加载过程中,垫层底部颗粒水平相对位移较大,垫层中上部颗粒水平相对位移较小,颗粒配位数及垫层孔隙率与碎石垫层隔震效果关系密切。
  • 地震时间分布模型是进行地震预测和地震危险性分析的重要理论基础,分析海域地震时间分布特征是进行海域地震危险性分析的重要研究内容。地震的发生在时间上符合泊松模型是当前概率地震危险性分析和地震预测中广泛接受的基本假设。泊松模型表示地震发生是相互独立的,即本次地震发生在时间上不受上次地震的影响,也不会对下次地震造成影响。地震时间独立分布模型包含泊松模型、双泊松模型、分段泊松模型、更新模型以及特征地震模型等(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.  Non-contact pile box composite foundation

    图  2  垫层数值模型

    Figure  2.  Numerical simulation model of bedding layer

    图  3  输入地震动加速度时程曲线

    Figure  3.  Input seismic recording acceleration timehistory curves

    图  4  非接触桩箱复合基础水平承载力变化曲线

    Figure  4.  Horizontal bearing capacity change curves of non-contact pile box composite foundation

    图  5  非接触桩箱复合基础水平荷载-位移关系曲线

    Figure  5.  Horizontal load-displacement relationship curves of non-contact pile box composite foundation

    图  6  不同时刻垫层内颗粒平均位移

    Figure  6.  Average displacement of particles of bedding layer under different moments

    图  7  不同时刻垫层内颗粒平均加速度

    Figure  7.  Average acceleration of particles of bedding layer under different moments

    图  8  不同时刻垫层内颗粒平均位移

    Figure  8.  Average particle displacement of bedding layer under different moments

    图  9  不同时刻垫层内颗粒平均加速度

    Figure  9.  Average acceleration of particles of bedding layer under different moments

    图  10  不同时刻垫层内颗粒平均位移

    Figure  10.  Average particle displacement of particles of bedding layer under different moments

    图  11  不同时刻垫层内颗粒平均加速度

    Figure  11.  Average acceleration of particles of bedding layer under different moments

    图  12  垫层平均孔隙率随时间变化曲线

    Figure  12.  Average porosity curves of bedding layer

    图  13  垫层平均配位数随时间变化曲线

    Figure  13.  Average matting coordination number of bedding layer

    图  14  输入地震动加速度时程曲线

    Figure  14.  Input acceleration timehistory curve

    图  15  垫层不同高度处平均孔隙率随时间变化曲线

    Figure  15.  Mean porosity curves at different heights of bedding layer

    图  16  垫层不同高度处颗粒平均配位数随时间变化曲线

    Figure  16.  Curve of mean particle coordination number at different heights of bedding layer

    表  1  垫层材料参数

    Table  1.   Bedding material parameters

    参数数值
    粒径/mm5.2,3.1~7.0,0.1~9.9
    厚度/mm20.0,30.0,40.0,60.0
    阻尼比0.055
    黏聚力/kPa0
    重度/(kN·m−316.51
    颗粒间摩擦系数0.57
    侧限压缩模量/MPa60.1
    下载: 导出CSV

    表  2  隔震率对比

    Table  2.   Comparison of isolation ratio

    方案编号垫层厚度/mm压力/kPa输入加速度/(m·s−2输出加速度/(m·s−2数值模拟减震率/%文献减震率/%误差/%
    1-12003.53.441.7050.657.111.4
    1-22008.33.402.3530.932.13.8
    1-33003.53.411.5255.436.8433.5
    1-43008.33.411.7847.844.117.7
    下载: 导出CSV
  • 董学武, 周世忠, 2004. 希腊里翁-安蒂里翁大桥的设计与施工. 世界桥梁, 32(4): 1—4

    Dong X. W. , Zhou S. Z. , 2004. Design and construction of Rion-Antirion bridge in Greece. World Bridges, 32(4): 1—4. (in Chinese)
    范立础, 王君杰, 2001. 桥梁抗震设计规范的现状与发展趋势. 地震工程与工程振动, 21(2): 70—77 doi: 10.3969/j.issn.1000-1301.2001.02.013

    Fan L. C. , Wang J. J. , 2001. Design code for earthquake-resistance of bridges: current situation and trend. Earthquake Engineering and Engineering Vibration, 21(2): 70—77. (in Chinese) doi: 10.3969/j.issn.1000-1301.2001.02.013
    李志强, 栾阳, 张倩等, 2018. 碎石垫层隔震性能离散元数值分析. 公路, 63(1): 158—164.
    刘宗华, 2012. 土耳其Izmit海湾大桥主塔基础设计与施工构思. 施工技术, 41(17): 11—13, 20

    Liu Z. H. , 2012. The conception of main tower foundation design and construction for Izmit bay bridge in turkey. Construction Technology, 41(17): 11—13, 20. (in Chinese)
    吕伟华, 缪林昌, 2013. 刚性桩复合地基桩土应力比计算方法. 东南大学学报(自然科学版), 43(3): 624—628

    Lv W. H. , Miao L. C. , 2013. Calculation method of pile-soil stress ratio of rigid pile composite foundation. Journal of Southeast University (Natural Science Edition), 43(3): 624—628. (in Chinese)
    南文文, 2015. 沉箱-垫层-桩复合深水基础水平承载性能研究. 南京: 东南大学, 43—51

    Nan W. W., 2015. Experimental study on bearing performance of caisson-cushion-pile deepwater composite foundation under lateral loads. Nanjing: Southeast University, 43—51. (in Chinese)
    陶景晖, 梁书亭, 龚维明等, 2009. 高层建筑刚性桩复合地基承载受力性状研究. 东南大学学报(自然科学版), 39(S2): 238—245

    Tao J. H. , Liang S. T. , Gong W. M. , et al. , 2009. Study on bearing behavior of composite foundation with rigid pile for high-rise buildings. Journal of Southeast University (Natural Science Edition), 39(S2): 238—245. (in Chinese)
    魏磊, 2013. 村镇建筑基础下碎石垫层隔震性能试验研究. 西安: 西安建筑科技大学, 21—35

    Wei L., 2013. An experimental research on seismic reduced function gravel layer under rural buildings foundation. Xi'an: Xi'an University of Architecture and Technology, 21—35. (in Chinese)
    姚国伟, 2011. 复合地基褥垫层减震研究. 广州: 广州大学, 38—41

    Yao G. W., 2011. Study on seismic reduced function of composite foundation cushion. Guangzhou: Guangzhou University, 38—41. (in Chinese)
    赵少伟, 窦远明, 郭蓉等, 2005. 基础下砂垫层隔震性能振动台试验研究. 河北工业大学学报, 34(3): 92—97 doi: 10.14081/j.cnki.hgdxb.2005.03.019

    Zhao S. W. , Dou Y. M. , Guo R. , et al. , 2005. An experimental study of isolating properties of sand cushion under the foundation by shaking table. Journal of Hebei University of Technology, 34(3): 92—97. (in Chinese) doi: 10.14081/j.cnki.hgdxb.2005.03.019
    Anastasopoulos I. , Gazetas G. , Loli M. , et al. , 2010. Soil failure can be used for seismic protection of structures. Bulletin of Earthquake Engineering, 8(2): 309—326. doi: 10.1007/s10518-009-9145-2
    Biesiadecki G. L., Dobry R., Leventis G. E., et al., 2004. Rion- Antirion Bridge foundations: A blend of design and construction innovation. In: Fifth International Conference on Case Histories in Geotechnical Engineering. New York: University of Missouri, 16.
    Combault J. , Pecker A. , Teyssandier J. P. , et al. , 2005. Rion-antirion bridge, Greece - concept, design, and construction. Structural Engineering International, 15(1): 22. doi: 10.2749/101686605777963387
    Dobry R. , Pecker A. , Mavroeidis G. P. , et al. , 2003. Damping/global energy balance in FE model of bridge foundation lateral response. Soil Dynamics and Earthquake Engineering, 23(6): 483—495. doi: 10.1016/S0267-7261(03)00050-2
    Lyngs J. H., Kasper T., Bertelsen K. S., 2013. Modelling of soil-structure interaction for seismic analyses of the Izmit Bay Bridge. In: Proceedings of the 18 th International Conference on Soil Mechanics and Geotechnical Engineering. Paris: ICSMGE, 763—768.
    Pfauth T. , Poling J. , 2017. Strain measurement in a specimen subjected to out-of-plane movement: Using an open-source digital image correlation-based tool. The Journal of Purdue Undergraduate Research, 7(1): 7.
    Yang D., Dobry R., Peck R. B., 2001. Foundation-soil-inclusion interaction modelling for rion-antirion bridge seismic analysis. In: Fourth International Conference on Recent Advances in Geotechnical Earthquake Engineering and Soil Dynamics. New York: University of Missouri, 620.
  • 期刊类型引用(1)

    1. 陆东华,陈鲲,魏佳彤. 台湾地区浅层地震活动性模型初步研究. 震灾防御技术. 2024(03): 494-503 . 本站查看

    其他类型引用(2)

  • 加载中
图(16) / 表(2)
计量
  • 文章访问数:  89
  • HTML全文浏览量:  77
  • PDF下载量:  8
  • 被引次数: 3
出版历程
  • 收稿日期:  2022-03-24
  • 刊出日期:  2023-08-31

目录

/

返回文章
返回