技术领域
[0001] 本
发明涉及近断层地震动评估领域，特别是一种海底近断层地震动的拟合方法。
背景技术
[0002] 近断层地震动特性的研究表明，近断层地震动具有区别于远场地震动的两个显著特征：一个是速度时程包含长周期、大峰值的脉冲；另一个是位移时程中可能包含永久地面位移引起的阶跃型脉冲。因此，位于近断层附近的结构会产生更为严重的地震破坏。
[0003] 由于目前已有的近断层地震动记录十分有限，为解决近断层地震动记录十分有限的问题，有关学者提出近断层地震动的拟合方法，基本思路为：高频成分的
加速度时程与低频成分的脉冲成分进行
叠加。以上人工近断层地震动拟合方法的出现，使得工程师能更准确地评估位于近断层附近的建筑结构的抗震性能，避免了以往仅以远场地震动来评估抗震性能而带来的危险低估。
[0004] 然而，已有的近断层地震动拟合方法是基于陆地场地近断层地震动记录而提出的拟合方法。由于
海水层的存在，处于海域环境的地震传播规律显著不同于陆地场地：(1)
海水层不仅会直接影响
地震波在近海场地中的传播；(2)海水层还会增大近海场地土层的
饱和度和孔隙水压
力，进而影响地震动的场地放大效应。因此，目前针对陆地场地的近断层地震动拟合方法不能适用于拟合海底近断层地震动。
[0005] 随着我国沿海经济的飞速发展，在海域、滨海地区的地震
断裂带附近或穿越断裂带修建大跨交通工程以及大型海洋工程等生命线工程的需求越来越大。例如刚刚建好的海文大桥即受活动断层影响十分严重。另一方面，目前已有的海底近断层地震动记录还较少，无法满目前海洋工程抗震分析的需求。
[0006] 目前尚无适用于海底近断层地震动的拟合方法，因此，亟需提出适用于海底近断层地震动的拟合方法，以便为跨海结构的抗震分析提供合适的地震动输入，以提高位于近断层附件的跨海结构抗震安全性。
发明内容
[0007] 有鉴于此，本发明的目的是提供一种海底近断层地震动的拟合方法，为对海底近断层结构的抗震性能分析提供更为准确合理的地震动输入。
[0008] 本发明采用以下方案实现：一种海底近断层地震动的拟合方法，包括以下步骤：
[0009] 步骤S1：计算近海场地的地震动传递函数H(ω)；
[0010] 步骤S2：拟合海底地震动的高频成分，得到近海域近断层地震动的高频加速度时程分量a3(t)；
[0011] 步骤S3：拟合海底地震动的低频脉冲成分，得到地震动的低频脉冲型加速度时程a5(t)；
[0012] 步骤S4：将高频加速度时程分量a3(t)与低频脉冲型加速度时程a5(t)叠加，得到海底近断层地震动。
[0013] 进一步地，所述步骤S1具体包括以下步骤：
[0014] 步骤S11：求基岩的动力
刚度矩阵[SR]；
[0015] 步骤S12：求海
底土层的动力刚度矩阵[SL]；
[0016] 步骤S13：求海水层的动力刚度矩阵[SW]；
[0017] 步骤S14：通过海水层动力刚度矩阵[SW]、海底土层动力刚度矩阵[SL]和基岩动力刚度矩阵[SR]组合得到海底场地的动力刚度矩阵；通过场地的动力刚度矩阵、场地位移向量和场地荷载向量确定场地动力平衡方程，求解方程得到海底层底部位移和基岩自由表面处位移之比即近海场地传递函数H(ω)。
[0018] 进一步地，所述步骤S11具体包括以下内容：
[0019] 根据场地中基岩的参数取值，计算基岩中地震波的传播速度；
[0020] 基岩中S波的传播速度表示为：
[0021]
[0022] 式中，G为基岩的
剪切模量，ρ为基岩的
密度；
[0023] 基岩中P波的传播速度表示为：
[0024]
[0025] 式中，G为基岩的剪切模量，ρ为基岩的密度，υ为基岩的泊松比；
[0026] 考虑了
迟滞阻尼的S波
波速计算公式如下：
[0027]
[0028] 式中，Vp为S波波速，ξ为阻尼比；
[0029] 考虑了迟滞阻尼的P波波速计算公式如下：
[0030]
[0031] 式中，Vp为P波波速，ξ为阻尼比；
[0032] 提供基岩中P波的入射
角度，则通过
相位速度的计算公式计算出相位速度；确定了相位速度，通过相位速度的计算公式陆续推导出其他介质中不同地震波的入射角度；
[0033] 相位速度的计算公式如下：
[0034]
[0035] 式中，c为相位速度； 和 分别为考虑了迟滞阻尼的P波和S波的波速；lx和mx分别为P波和S波入射角度的余弦值；上标W、L和R分别表示海水层、海底土层和基岩；
[0036] 利用相位速度、基岩中地震波传播速度和基岩的参数取值，得到基岩的动力刚度矩阵；
[0037] 平面外SH波作用下基岩的动力刚度矩阵如下：
[0038]
[0039] 式中,GR为基岩剪切模量，ξR为基岩阻尼比，tR为平面外SH波在基岩处入射角度的正切值，ω为角
频率，c为相位速度；
[0040] 平面内P波和SV波作用下基岩的动力刚度矩阵如下：
[0041]
[0042] 其中，
[0043]
[0044]
[0045]
[0046] 式中， 为考虑了迟滞阻尼的基岩剪切模量，t'R为平面内SV波在基岩处入射角度的正切值，s'R为平面内P波在基岩处入射角度的正切值，ω为角频率，c为相位速度。
[0047] 进一步地，所述步骤S12具体包括以下内容：
[0048] 通过近海场地的海底土层饱和度计算含水土体的泊松比，含水土体的泊松比表达式如下：
[0049]
[0050] 式中，υ'、G分别为土骨架的泊松比和剪切模量，α、M分别为与土颗粒和孔隙
流体压缩性相关的系数，通过土体孔隙率n和土体饱和度Sr计算；
[0051] 用含水土体泊松比计算海底土层P波传播速度，海底土层P波传播速度表达式如下：
[0052]
[0053] 式中，G为土骨架的剪切模量，α、M分别为与土颗粒和孔隙流体压缩性相关的系数，λ为土骨架的
拉梅常数，ρ＝(1-n)ρs+nρf为含水土体的密度，ρs和ρf分别为土颗粒和孔隙流体的密度；
[0054] 根据场地中海底土层的参数取值，计算海底土层中S波的传播速度，S波波速计算公式表示为：
[0055]
[0056] 式中，G为土骨架的剪切模量，ρ为含水土体的密度；
[0057] 考虑了迟滞阻尼的S波波速计算公式如下：
[0058]
[0059] 式中，Vp为S波波速，ξ为阻尼比；
[0060] 考虑了迟滞阻尼的P波波速计算公式如下：
[0061]
[0062] 式中，Vp为P波波速，ξ为阻尼比；
[0063] 通过相位速度计算公式，求出各自的入射角；相位速度的计算公式如下：
[0064]
[0065] 式中，c为相位速度； 和 分别为考虑了迟滞阻尼的P波和S波的波速；lx和mx分别为P波和S波入射角度的余弦值；上标W、L和R分别表示海水层、海底土层和基岩；
[0066] 利用海底土层中地震波的入射角和海底土层的参数取值，得出海底土层的动力刚度矩阵；某一土层的动力刚度矩阵如下：
[0067]
[0068] 其中，每一个kL与角频率、相位速度、地震波入射角度及场地参数相关；
[0069] 将每一个土层的动力刚度矩阵组合得到海底整个土层的动力刚度矩阵。
[0070] 进一步地，所述步骤S13具体包括以下内容：根据海水层的参数取值计算海水层中P波传播速度，海水层中P波传播速度表达式为：
[0071]
[0072] 式中，K为海
水体积模量，ρ为海水密度；
[0073] 通过公式 求得的相位速度求出海水层P波入射角，式中，c为相位速度； 和 分别为考虑了迟滞阻尼的P波和S波的波速；lx和mx分别为P波和S波入射角度的余弦值；上标W、L和R分别表示海水层、海底土层和基岩；
[0074] 再结合海水层深度和海水层阻尼比，得出海水层的动力刚度矩阵，见公式(20)。
[0075]
[0076] 其中，每一个kW与角频率、相位速度、地震波入射角度及地震波传播速度、海水层剪切模量、海水层厚度、海水层阻尼比等参数相关。
[0077] 进一步地，所述步骤S14具体包括以下内容：
[0078] 在平面内P波和SV波共同作用下的场地动力平衡方程表示为：
[0079] [SP-SV]{uP-SV}＝{PP-SV} (21.1)
[0080] 式中，[SP-SV]为海底场地平面内方向的动力刚度矩阵，{uP-SV}为P波和SV波共同作用下的场地位移，{PP-SV}为P波和SV波共同作用下的荷载向量；
[0081] 在平面外SH波作用下的场地动力平衡方程表示为：
[0082] [SSH]{uSH}＝{PSH} (21.2)
[0083] 式中，[SSH]为海底场地平面外方向的动力刚度矩阵，{uSH}为SH波作用下的场地位移，{PSH}为SH波作用下的荷载向量；
[0084] 通过公式(21.1)和公式(21.2)的两个场地动力平衡方程，求得近海场地地震动的三个方向的传递函数；近海场地地震动传递函数一般表达式如下：
[0085]
[0086] 式中，ub为海水层底部位移幅值，uo为基岩自由表面处位移幅值。
[0087] 进一步地，所述步骤S2的具体包括以下步骤：
[0088] 步骤S21：计算海底土层表面的幅值谱；
[0089] 采用修正Tajimi-Kanai
功率谱密度模型对基岩自由表面处三个方向的地震动功率谱密度函数进行模拟，表达式如下：
[0090]
[0091] 式中，Sbr(ω)为基岩自由表面地震动功率谱密度函数；|Hp(ω)|为高通滤波函数；Sg(ω)为Tajimi-Kanai功率谱密度函数；ωf和ξf分别为高通滤波函数的中心频率和阻尼比；ωg和ξg分别Tajimi-Kanai功率谱密度函数的中心频率和阻尼比；
[0092] 由基岩自由表面功率谱Sbr(ω)计算基岩自由表面傅里叶幅值谱Ar(ω)，表达式如下：
[0093] Abr(ω)＝[4Sbr(ω)Δω]1/2 (24)
[0094] 式中，Δω为频率间隔，Δω＝2π×
采样频率/FFT长度；
[0095] 确定基岩自由表面处幅值谱后，通过近海场地传递函数H(ω)，求得近海场地海底土层表面的幅值谱As1(ω)表达式如下：
[0096] As1(ω)＝|H(ω)|Abr(ω) (25)
[0097] 式中，Abr(ω)、As1(ω)分别表示基岩自由表面、近海场地海底土层处地震动的傅里叶幅值谱，H(ω)表示计算得出的近海场地地震动传递函数；
[0098] 步骤S22：通过matlab的rand函数生成[0,2π]间均匀分布的随机相位谱φ1(ω)；
[0099] 步骤S23：基于式(26)进行傅里叶逆变换，得到高频加速度时程a1(t)。
[0100]
[0101] 步骤S24：对求得的高频加速度时程a1(t)乘以强度包络线g(t)，得到非平稳加速度时程a2(t)，表达式如下：
[0102] a2(t)＝a1(t)g(t) (27.1)
[0103] 强度包络线由上升时间段、延续时间段和下降时间段的包络线组成。延续时间段包络线的线形是幅值为1的水平线，上升时间段和下降时间段包络线的线形可以是直线、抛物线或指数曲线。
[0104] 上升时间段直线、抛物线或指数曲线的包络线函数的表达式分别为：
[0105] 直线：
[0106] 抛物线：
[0107] 指数曲线：
[0108] 式中，t为时间变量；T1为上升时间长度；c1为上升时间段指数曲线包络线的线形参数。
[0109] 下降时间段直线、抛物线或指数曲线的包络线函数的表达式分别为：
[0110] 直线：
[0111] 抛物线：
[0112] 指数曲线：
[0113] 式中，t为时间变量；T2为0时刻至延续时间段结束时刻的时间长度；c2为下降时间段指数曲线包络线的线形参数；T为地震波总的时间长度。
[0114] 步骤S25：对非平稳加速度时程a2(t)进行傅里叶变换，得到幅值谱As2(ω)和相位谱φ2(ω)。将低于临界频率的幅值设为0，得到调整后的幅值谱As3(ω)。基于式(27.8)进行傅里叶逆变换，得到海域近断层地震动的高频加速度时程分量a3(t)。
[0115]
[0116] 进一步地，所述步骤S3的具体包括以下步骤：
[0117] 步骤S31：通过等效速度脉冲模型生成等效速度脉冲时程，对等效速度脉冲时程求导数得到等效加速度时程：
[0118] 等效速度脉冲模型如下：
[0119] v(t)＝vp·ω(t)·cos[2πfp(t-t1)] 0≤t≤T (28)
[0120] 式中，vp为速度脉冲的峰值，fp是速度脉冲的频率，速度脉冲的周期是Tp＝1/fp，-2πfPt1是速度脉冲的相位角，T是速度时程的持续时间，速度时程的包络函数ω(t)按下式取值：
[0121]
[0122] 式中，t0是包络函数的峰值发生时刻，γ表示峰值衰减速率；
[0123] 根据速度脉冲周期的统计公式，得到速度脉冲的周期TP和频率fP,统计公式如下：
[0124] logTp＝-2.9+0.5Mw (30)
[0126] fP＝1/TP (31)
[0127] 根据速度脉冲峰值的统计公式，得到速度脉冲的峰值vP，统计公式如下：
[0128] log10vp＝-2.22+0.69Mw-0.58log10R (32)
[0129] 式中，Mw为矩震级,R为场地到断层的最短距离；
[0130] 根据具体的脉冲型记录确定形状参数γ、t0和t1的取值；将参数TP、fP、vP、γ、t0和t1代入等效速度脉冲模型，由等效速度脉冲模型生成等效速度脉冲时程；对等效速度时程求导得到等效加速度时程；
[0131] 步骤S32：初步得到低频脉冲型加速度时程：
[0132] 将等效加速度时程进行傅里叶变换，得到低频加速度时程的幅值谱As4(ω)和相位谱φ4(ω)；通过近海场地传递函数H(ω)，求得海底土层表面的幅值谱如下：
[0133] As5(ω)＝|H(ω)|As4(ω) (33)
[0134] 式中，As4(ω)、As5(ω)分别表示基岩自由表面、近海场地海底土层处地震动的傅里叶幅值谱，H(ω)表示计算得出的近海场地地震动传递函数；
[0135] 经过式(34)的傅里叶逆变换，最终得到地震动的低频脉冲型加速度时程a4(t)；
[0136]
[0137] 进一步地，所述步骤S4的具体包括以下步骤：
[0138] 步骤S41：重新生成低频脉冲型加速度时程；
[0139] 定义低频脉冲加速度峰值的到达时刻与高频加速度峰值的到达时刻相同；根据这一定义，在时间轴上平行移动低频加速度时程，使得低频加速度时程的峰值时刻与高频加速度时程的峰值时刻相重合；
[0140] 在时间轴上平移加速度时程a4(t)，得到低频脉冲成分的加速度时程a5(t)；
[0141] 步骤S42：将平移后的低频脉冲型加速度时程分量a5(t)与高频加速度时程分量a3(t)进行叠加，得到近海域近断层地震动加速度时程a(t)，公式如下:
[0142] a(t)＝a5(t)+a3(t) (35)
[0143] 与
现有技术相比，本发明具有以下有益效果：
[0144] 本发明能够切实地解决目前跨海结构抗震分析缺乏的近断层地震动记录的
瓶颈，为位于海底断层附近的海洋结构抗震性能分析提供更为准确合理的地震动输入，有助于推动和提高海洋结构的抗震性能评估。
附图说明
具体实施方式
[0146] 下面结合附图及实施例对本发明做进一步说明。
[0147] 应该指出，以下详细说明都是例示性的，旨在对本
申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属技术领域的普通技术人员通常理解的相同含义。
[0148] 需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本
说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。
[0149] 如图1所示，本实施例提供一种海底近断层地震动的拟合方法，包括以下步骤：
[0150] 步骤S1：计算近海场地的地震动传递函数H(ω)；
[0151] (1)计算海水层的动力刚度矩阵
[0152] 考虑海水层的存在会影响地震动的场地放大效应，根据海水层的参数取值计算海水层中P波传播速度。
[0153] 通过已求得的相位速度求出海水层P波入射角，再结合海水层深度和海水层阻尼比，可以得出海水层的动力刚度矩阵。
[0154] (2)计算海底土层的动力刚度矩阵
[0155] 考虑到海底多孔土层的存在会影响地震动的场地放大效应，通过近海场地的海底土层饱和度计算含水土体的泊松比，进一步用含水土体泊松比计算海底土层P波传播速度；
[0156] 根据近海场地海底土层参数，计算海底土层中SH波和SV波的传播速度，通过对应的相位速度求出各自的入射角；
[0157] 利用海底土层中地震波的入射角和海底土层的参数取值，可以得出海底土层的动力刚度矩阵。
[0158] (3)计算基岩的动力刚度矩阵
[0159] 根据场地中基岩的参数取值，计算基岩中地震波的传播速度；
[0160] 假设基岩中地震波的入射角度，通过地震波的传播速度计算出相位速度；
[0161] 利用相位速度、基岩中地震波传播速度和基岩的参数取值，可以得出基岩的动力刚度矩阵。
[0162] (4)计算近海场地传递函数
[0163] 海底场地的动力刚度矩阵可通过以上求得的基岩动力刚度矩阵、海底土层动力刚度矩阵和海水层动力刚度矩阵组合得出，然后采用场地动力平衡方程，可求得近海场地传递函数。
[0164] 步骤S2：拟合海底地震动的高频成分，得到近海域近断层地震动的高频加速度时程分量a3(t)；
[0165] (1)计算海底土层表面的幅值谱
[0166] 采用功率谱密度函数模型对基岩自由表面处地震动功率谱进行模拟，通过功率谱与傅里叶幅值谱的关系式求出基岩自由表面处傅里叶幅值谱。
[0167] 确定基岩自由表面处幅值谱后，通过近海场地传递函数，求得海底土层表面的幅值谱。
[0168] (2)由随机
相位差谱模型计算随机相位谱。
[0169] (3)基于求得的幅值谱和相位谱，采用傅里叶逆变换，得到海底近断层地震动的高频加速度时程分量，并乘上包络函数，得到非平稳的高频加速度时程。
[0170] 步骤S3：拟合海底地震动的低频脉冲成分，得到地震动的低频脉冲型加速度时程a5(t)；
[0171] (1)通过等效速度脉冲模型生成等效速度脉冲时程，对等效速度脉冲时程求导数得到等效加速度时程。
[0172] (2)考虑海水层对低频成分的影响
[0173] 对等效加速度时程进行傅里叶变换，得到低频加速度时程的相位谱和基岩处幅值谱。
[0174] 基于基岩处幅值谱与近海场地的传递函数，求得海底土层表面的幅值谱。
[0175] 基于傅里叶变换得到的相位谱与海底土层表面的幅值谱，经过傅里叶逆变换，得到低频脉冲型加速度时程。
[0176] 步骤S4：将高频加速度时程分量a3(t)与低频脉冲型加速度时程a5(t)叠加，得到海底近断层地震动。
[0177] (1)重新生成低频脉冲型加速度时程
[0178] 在时间轴上平移低频加速度时程，使得低频加速度时程的峰值时刻与高频加速度时程的峰值时刻相重合。重新生成速度脉冲加速度时程，即得到海底近断层地震动的低频脉冲型加速度时程分量。
[0179] (2)将低频脉冲型加速度时程分量与高频加速度时程分量进行叠加，得到海底近断层地震动加速度模拟时程。
[0180] 在本实施例中，所述步骤S1具体包括以下步骤：
[0181] 步骤S11：求基岩的动力刚度矩阵[SR]；
[0182] 步骤S12：求海底土层的动力刚度矩阵[SL]；
[0183] 步骤S13：求海水层的动力刚度矩阵[SW]；
[0184] 步骤S14：通过海水层动力刚度矩阵[SW]、海底土层动力刚度矩阵[SL]和基岩动力刚度矩阵[SR]组合得到海底场地的动力刚度矩阵；通过场地的动力刚度矩阵、场地位移向量和场地荷载向量确定场地动力平衡方程，求解方程得海底层底部位移和基岩自由表面处位移之比即近海场地传递函数H(ω)。
[0185] 在本实施例中，所述步骤S11具体包括以下内容：
[0186] 根据场地中基岩的参数取值，计算基岩中地震波的传播速度；
[0187] 基岩中S波的传播速度表示为：
[0188]
[0189] 式中，G为基岩的剪切模量，ρ为基岩的密度；
[0190] 基岩中P波的传播速度表示为：
[0191]
[0192] 式中，G为基岩的剪切模量，ρ为基岩的密度，υ为基岩的泊松比；
[0193] 考虑了迟滞阻尼的S波波速计算公式如下：
[0194]
[0195] 式中，Vp为S波波速，ξ为阻尼比；
[0196] 考虑了迟滞阻尼的P波波速计算公式如下：
[0197]
[0198] 式中，Vp为P波波速，ξ为阻尼比；
[0199] 提供基岩中P波的入射角度，则通过相位速度的计算公式计算出相位速度；确定了相位速度，通过相位速度的计算公式陆续推导出其他介质中不同地震波的入射角度；
[0200] 相位速度的计算公式如下：
[0201]
[0202] 式中，c为相位速度； 和 分别为考虑了迟滞阻尼的P波和S波的波速；lx和mx分别为P波和S波入射角度的余弦值；上标W、L和R分别表示海水层、海底土层和基岩；
[0203] 利用相位速度、基岩中地震波传播速度和基岩的参数取值，得到基岩的动力刚度矩阵；
[0204] 平面外SH波作用下基岩的动力刚度矩阵如下：
[0205]
[0206] 式中,GR为基岩剪切模量，ξR为基岩阻尼比，tR为平面外SH波在基岩处入射角度的正切值，ω为角频率，c为相位速度；
[0207] 平面内P波和SV波作用下基岩的动力刚度矩阵如下：
[0208]
[0209] 其中，
[0210]
[0211]
[0212]
[0213] 式中， 为考虑了迟滞阻尼的基岩剪切模量，t'R为平面内SV波在基岩处入射角度的正切值，s'R为平面内P波在基岩处入射角度的正切值，ω为角频率，c为相位速度。
[0214] 在本实施例中，所述步骤S12具体包括以下内容：通过近海场地的海底土层饱和度计算含水土体的泊松比，含水土体的泊松比表达式如下：
[0215]
[0216] 式中，υ'、G分别为土骨架的泊松比和剪切模量，α、M分别为与土颗粒和孔隙流体压缩性相关的系数，通过土体孔隙率n和土体饱和度Sr计算；
[0217] 用含水土体泊松比计算海底土层P波传播速度，海底土层P波传播速度表达式如下：
[0218]
[0219] 式中，G为土骨架的剪切模量，α、M分别为与土颗粒和孔隙流体压缩性相关的系数，λ为土骨架的拉梅常数，ρ＝(1-n)ρs+nρf为含水土体的密度，ρs和ρf分别为土颗粒和孔隙流体的密度；
[0220] 根据场地中海底土层的参数取值，计算海底土层中S波的传播速度，S波波速计算公式表示为：
[0221]
[0222] 式中，G为土骨架的剪切模量，ρ为含水土体的密度；
[0223] 考虑了迟滞阻尼的S波波速计算公式如下：
[0224]
[0225] 式中，Vp为S波波速，ξ为阻尼比；
[0226] 考虑了迟滞阻尼的P波波速计算公式如下：
[0227]
[0228] 式中，Vp为P波波速，ξ为阻尼比；
[0229] 通过相位速度计算公式，求出各自的入射角；相位速度的计算公式如下：
[0230]
[0231] 式中，c为相位速度； 和 分别为考虑了迟滞阻尼的P波和S波的波速；lx和mx分别为P波和S波入射角度的余弦值；上标W、L和R分别表示海水层、海底土层和基岩；
[0232] 利用海底土层中地震波的入射角和海底土层的参数取值，得出海底土层的动力刚度矩阵。某一土层的动力刚度矩阵如下：
[0233]
[0234] 其中，每一个kL与角频率、相位速度、地震波入射角度及场地参数相关；
[0235] 将每一个土层的动力刚度矩阵组合得到海底整个土层的动力刚度矩阵。
[0236] 在本实施例中，所述步骤S13具体包括以下内容：根据海水层的参数取值计算海水层中P波传播速度，海水层中P波传播速度表达式为：
[0237]
[0238] 式中，K为海水体积模量，ρ为海水密度；
[0239] 通过公式(5) 求得的相位速度求出海水层P波入射角，式中，c为相位速度； 和 分别为考虑了迟滞阻尼的P波和S波的波速；lx和mx分别为P波和S波入射角度的余弦值；上标W、L和R分别表示海水层、海底土层和基岩；
[0240] 再结合海水层深度和海水层阻尼比，得出海水层的动力刚度矩阵，见公式(20)。
[0241]
[0242] 其中，每一个kW与角频率、相位速度、地震波入射角度及地震波传播速度、海水层剪切模量、海水层厚度、海水层阻尼比等参数相关。
[0243] 在本实施例中，所述步骤S14具体包括以下内容：
[0244] 在平面内P波和SV波共同作用下的场地动力平衡方程表示为：
[0245] [SP-SV]{uP-SV}＝{PP-SV} (21.1)
[0246] 式中，[SP-SV]为海底场地平面内方向的动力刚度矩阵，{uP-SV}为P波和SV波共同作用下的场地位移，{PP-SV}为P波和SV波共同作用下的荷载向量；
[0247] 在平面外SH波作用下的场地动力平衡方程表示为：
[0248] [SSH]{uSH}＝{PSH} (21.2)
[0249] 式中，[SSH]为海底场地平面外方向的动力刚度矩阵，{uSH}为SH波作用下的场地位移，{PSH}为SH波作用下的荷载向量；
[0250] 通过公式(21.1)和公式(21.2)的两个场地动力平衡方程，求得近海场地地震动的三个方向的传递函数；近海场地地震动传递函数一般表达式如下：
[0251]
[0252] 式中，ub为海水层底部位移幅值，uo为基岩自由表面处位移幅值。
[0253] 在本实施例中，所述步骤S2的具体包括以下步骤：
[0254] 步骤S21：计算海底土层表面的幅值谱；
[0255] 采用修正Tajimi-Kanai功率谱密度模型对基岩自由表面处三个方向的地震动功率谱密度函数进行模拟，表达式如下：
[0256]
[0257] 式中，Sbr(ω)为基岩自由表面地震动功率谱密度函数；|Hp(ω)|为高通滤波函数；Sg(ω)为Tajimi-Kanai功率谱密度函数；ωf和ξf分别为高通滤波函数的中心频率和阻尼比；ωg和ξg分别Tajimi-Kanai功率谱密度函数的中心频率和阻尼比；
[0258] 由基岩自由表面功率谱Sbr(ω)计算基岩自由表面傅里叶幅值谱Ar(ω)，表达式如下：
[0259] Abr(ω)＝[4Sbr(ω)Δω]1/2 (24)
[0260] 式中，Δω为频率间隔，Δω＝2π×
采样频率/FFT长度；
[0261] 确定基岩自由表面处幅值谱后，通过近海场地传递函数H(ω)，求得近海场地海底土层表面的幅值谱As1(ω)表达式如下：
[0262] As1(ω)＝|H(ω)|Abr(ω) (25)
[0263] 式中，Abr(ω)、As1(ω)分别表示基岩自由表面、近海场地海底土层处地震动的傅里叶幅值谱，H(ω)表示计算得出的近海场地地震动传递函数；
[0264] 步骤S22：通过matlab的rand函数生成[0,2π]间均匀分布的随机相位谱φ1(ω)；
[0265] 步骤S23：基于式(26)进行傅里叶逆变换，得到高频加速度时程a1(t)。
[0266]
[0267] 步骤S24：对求得的高频加速度时程a1(t)乘以强度包络线g(t)，得到非平稳加速度时程a2(t)，表达式如下：
[0268] a2(t)＝a1(t)g(t) (27.1)
[0269] 强度包络线由上升时间段、延续时间段和下降时间段的包络线组成。延续时间段包络线的线形是幅值为1的水平线，上升时间段和下降时间段包络线的线形可以是直线、抛物线或指数曲线。
[0270] 上升时间段直线、抛物线或指数曲线的包络线函数的表达式分别为：
[0271] 直线：
[0272] 抛物线：
[0273] 指数曲线：
[0274] 式中，t为时间变量；T1为上升时间长度；c1为上升时间段指数曲线包络线的线形参数。
[0275] 下降时间段直线、抛物线或指数曲线的包络线函数的表达式分别为：
[0276] 直线：
[0277] 抛物线：
[0278] 指数曲线： 式中，t为时间变量；T2为0时刻至延续时间段结束时刻的时间长度；c2为下降时间段指数曲线包络线的线形参数；
T为地震波总的时间长度。
[0279] 步骤S25：对非平稳加速度时程a2(t)进行傅里叶变换，得到幅值谱As2(ω)和相位谱φ2(ω)。将低于临界频率的幅值设为0，得到调整后的幅值谱As3(ω)。基于式(27.8)进行傅里叶逆变换，得到海域近断层地震动的高频加速度时程分量a3(t)。
[0280]
[0281] 在本实施例中，所述步骤S3的具体包括以下步骤：
[0282] 步骤S31：通过等效速度脉冲模型生成等效速度脉冲时程，对等效速度脉冲时程求导数得到等效加速度时程：
[0283] 模拟低频脉冲型加速度时程分量时选取等效速度脉冲模型，因为该模型采用连续的单一函数形式表达速度时程，极大地方便了速度脉冲的曲线拟合以及其参数的确定，等效速度脉冲模型如下：
[0284] v(t)＝vp·ω(t)·cos[2πfp(t-t1)] 0≤t≤T (28)
[0285] 式中，vp为速度脉冲的峰值，fp是速度脉冲的频率，速度脉冲的周期是Tp＝1/fp，-2πfPt1是速度脉冲的相位角，T是速度时程的持续时间，速度时程的包络函数ω(t)按下式取值：
[0286]
[0287] 式中，t0是包络函数的峰值发生时刻，γ表示峰值衰减速率；
[0288] 根据速度脉冲周期的统计公式，得到速度脉冲的周期TP和频率fP,统计公式如下：
[0289] logTp＝-2.9+0.5Mw (30)
[0290] 式中，Mw为矩震级；
[0291] fP＝1/TP (31)
[0292] 根据速度脉冲峰值的统计公式，得到速度脉冲的峰值vP，统计公式如下：
[0293] log10vp＝-2.22+0.69Mw-0.58log10R (32)
[0294] 式中，Mw为矩震级,R为场地到断层的最短距离；
[0295] 根据具体的脉冲型记录确定形状参数γ、t0和t1的取值；将参数Tp、fp、vp、γ、t0和t1代入等效速度脉冲模型，由等效速度脉冲模型生成等效速度脉冲时程；对等效速度时程求导得到等效加速度时程；
[0296] 步骤S32：初步得到低频脉冲型加速度时程：
[0297] 将等效加速度时程进行傅里叶变换，得到低频加速度时程的幅值谱As4(ω)和相位谱φ4(ω)；通过近海场地传递函数H(ω)，求得海底土层表面的幅值谱如下：
[0298] As5(ω)＝|H(ω)|As4(ω) (33)
[0299] 式中，As4(ω)、As5(ω)分别表示基岩自由表面、近海场地海底土层处地震动的傅里叶幅值谱，H(ω)表示计算得出的近海场地地震动传递函数；
[0300] 经过式(34)的傅里叶逆变换，最终得到地震动的低频脉冲型加速度时程a4(t)；
[0301]
[0302] 在本实施例中，所述步骤S4的具体包括以下步骤：
[0303] 步骤S41：重新生成低频脉冲型加速度时程；
[0304] 定义低频脉冲加速度峰值的到达时刻与高频加速度峰值的到达时刻相同；根据这一定义，在时间轴上平行移动低频加速度时程，使得低频加速度时程的峰值时刻与高频加速度时程的峰值时刻相重合；
[0305] 在时间轴上平移加速度时程a4(t)，得到低频脉冲成分的加速度时程a5(t)；
[0306] 步骤S42：将平移后的低频脉冲型加速度时程分量a5(t)与高频加速度时程分量a3(t)进行叠加，得到近海域近断层地震动加速度时程a(t)，公式如下:
[0307] a(t)＝a5(t)+a3(t) (35)
[0308] 以上所述仅为本发明的较佳实施例，凡依本发明申请
专利范围所做的均等变化与修饰，皆应属本发明的涵盖范围。
