技术领域
[0001] 本发明属于地震技术领域，尤其涉及一种宽频带有限震源模型构建方法。
相关背景技术
[0002] 合理的地震动输入是结构抗震计算可靠性的必要前提。对于大坝、核电站等重要的结构，需要采用地震动时程作为输入。鉴于真实的地震动记录数量较少，尚不能满足工程需求，因此需要采用人工模拟地震动作为补充。
[0003] 目前，大坝场址地震动的合成方法多采用基于随机方法的随机有限震源法。该优点是能够考虑地震动的统计特性，对于模拟高频地震动效果较好。然而，采用随机方法无法考虑地震动三分量的差异，且难以计入局部地形起伏对地震动模拟的影响，并且对传播路径及场地的描述过于简化，无法反映其真实的特性。若采用从震源出发的直接数值模拟方法，如有限元法，谱元法等，则能较好的解决这些问题。但由于反演方法和数据频带的限制，目前仅能获得历史地震较为粗糙的震源特性。若想获得工程场址可用的宽频带地震动，必须构造能够反映地震破裂细节的震源模型，以生成更加合理的地震动高频成分。
具体实施方式
[0070] 下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
[0071] 实施例
[0072] 本发明提出了一种能够表征真实地震高频成分的宽频带有限震源模型构建方法，该方法通过多段不同的高斯子函数，将地震发生时的滑动速率模型分解为能够辐射高频能量的初始破裂过程和尾部持续破裂过程，近似模拟真实地震发生时的震源破裂过程，采用该方法，不仅能够约束震源错动的总持时，又能分别反映滑动过程中高频和低频能量的辐射水平，因此能够生成满足工程需求的宽频带地震动。如图1所示，一种宽频带有限震源模型构建方法，其实现方法如下：
[0073] S1、确定震源参数；
[0074] 本实施例中，根据场址近场区地震构造评价及断裂活动性特征确定发震断层基本信息。主要包括可能发震的矩震级Mw，断层的埋深htop、走向φ、滑动角λ、倾角δ，以及断层的长L和宽W，并将断层破裂面划分为500m*500m的子源。以及获取断层附近地壳结构，包括密度，p波及s波波速。
[0075] S2、根据震源参数，通过计算宽频带滑动量模型波数域低频分布和高频成分，构建宽频带滑动量模型，其实现方法如下：
[0076] S201、基于历史地震的反演位错模型或凹凸体模型，构建震源滑动量基础模型；
[0077] S202、根据震源滑动量基础模型以及震源参数，计算得到宽频带滑动量模型波数域低频分布；
[0078] S203、按照k2衰减规律计算宽频带滑动量模型波数域高频成分；
[0079] S204、将步骤S202与步骤S203的结果进行叠加，得到宽频带滑动量模型。
[0080] 本实施例中，确定宽频带滑动量模型包括确定震源滑动量基础模型以及补充滑动量模型高频细节。
[0081] 本实施例中，确定宽频带滑动量模型包括确定震源滑动量基础模型，如下：
[0082] 对于曾经发生过大震的断层，可采用历史地震的反演位错模型在子源上插值获得震源滑动量基础模型。除此之外，基础模型还可采用凹凸体模型。
[0083] 采用凹凸体模型时，凹凸体的总面积Sa可按下式计算：
[0084]
[0085] 其中，β表示剪切波速，vR表示断层破裂传播速度，S＝WL表示破裂的总面积，M0和分别表示标量地震矩和加速度震源谱幅值，可由如下经验公式计算：
[0086] M0＝101.5×Mw+16.05
[0087]
[0088] 断层总的平均位错可由标量地震矩M0估计：
[0089]
[0090] 其中，μ表示剪切模量。
[0091] 凹凸体的平均位错 分布与凹凸体数量Na有关：
[0092]
[0093] 凹凸体的数量可由活动断层分区确定，一般为1～3个，平均为2.6个。
[0094] 本实施例中，补充滑动量模型高频细节，如下：
[0095] (1)由基础模型计算得到宽频带滑动量模型波数域低频分布。
[0096] 若滑动量基础模型在断层平面的空间分布为uB(x,y)，则可采用二维傅立叶变换计算uB(x,y)在波数域的分布UB(ks,kd)。其中，x,y分别为沿走向和倾向的坐标，ks,kd分别为沿走向和倾向的波数。宽频带滑动量模型的低频成分可由滑动量基础模型经低通滤波得到：
[0097] UL(ks,kd)＝UB(ks,kd)F
[0098]
[0099] 其中，cs,cd分别表示沿走向和倾向的拐角波数，若基本模型采用凹凸体模型，cs,cd可分别取凹凸体沿走向或倾向最小尺寸的一半，此时滤波器指数N为1；若采用反演模型，cs,cd可分别取反演模型沿走向和倾向的最小子源分辨率，此时滤波器指数N为4。
[0100] (2)按照k2衰减规律计算宽频带滑动量模型的高频成分。
[0101] 假定滑动量模型在波数域服从k2衰减规律，则宽频带滑动量模型的高频成分可按下式确定：
[0102]
[0103] 其中，θ表示随机相位，在[0,2π]间均匀随机取值，as,ad分别表示沿走向和倾向的相关长度：
[0104] log as＝0.5Mw‑2.5
[0105]
[0106] A2(ks,kd)表示Von Karmen相关函数的波数功率谱密度，可由下式计算：
[0107]
[0108] 其中， 表示归一化的波数，H表示Hurst指数因子，H＝0.75。
[0109] 最终宽频带滑动量模型波数谱由高频部分和低频部分叠加得到：
[0110] U(ks,kd)＝UL(ks,kd)+UH(ks,kd)
[0111] 其中，U(ks,kd)表示宽频带滑动量模型。
[0112] 将U(ks,kd)进行二维逆傅里叶变换，可得最终的宽频带滑动量模型空间分布u(x,y)。如图2所示，其给出了某个情景地震采用反演模型并补充震源细节后得到的宽频带滑动量模型，及破裂传播前锋到达时间(等值线，间隔为2s)。
[0113] S3、根据宽频带滑动量模型，计算得到子源总上升时间模型；
[0114] 本实施例中，得到源宽频带滑动量模型后，可据此计算子源总上升时间模型。整个断层面的子源平均上升时间τAve可根据标量地震矩M0估计：
[0115]
[0116] 其中，aτ表示与δ相关的系数，δ表示倾角：
[0117]
[0118] 当δ在45°～60°之间时，aτ由线性插值确定。
[0119] 断层面上第i个子源的总上升时间τi与子源滑动量ui以及子源所在深度z有关：
[0120]
[0121] 当子源所在深度z在5～8km之间时，τi由线性插值确定。上式中，kτ表示系数，使得所有子源的总上升时间平均值与τAve相同。如图3所示，其给出了子源总上升时间的空间分布。
[0122] S4、根据宽频带滑动量模型，计算子源破裂前锋分布模型；
[0123] 本实施例中，断层破裂的传播速度Vr由断层的剪切波速和深度确定：
[0124]
[0125] 当子源所在深度z在5～8km之间时，Vr由线性插值确定。第i个子源破裂前锋的平均到达时间Ti0可通过解断层平面上的二维波动方程获得，其最终到达时间TiF与子源滑动量有关：
[0126]
[0127] 其中，uAve表示所有子源滑动量的均值，ui表示子源滑动量，但应大于0.05uAve，umax表示子源滑动量最大值，Δt表示与标量地震矩M0相关的调整系数：
[0128]
[0129] 如图2所示，其中的等值线给出了断层破裂前锋到达时间的分布。
[0130] S5、根据子源总上升时间模型、子源破裂前锋分布模型以及宽频带滑动量模型对子源破裂过程进行拆分，其实现方法如下：
[0131] S501、通过拟合历史地震记录，确定震源破裂过程中脉冲时间比例以及子过程脉冲能量占比；
[0132] S502、将子源总上升时间乘以每个子过程的脉冲时间比例，计算子过程的持续时间，得到子过程上升时间分布；
[0133] S503、结合子过程的持续时间和子源破裂前锋分布模型，计算每个子过程的开始时间，得到每个子过程的破裂前锋分布；
[0134] S504、将子源滑动量乘以子过程脉冲能量占比，得到子过程的滑动量分布，其中，若拆分前的震源共有n个子源，每个子源拆分为m个子过程，则拆分后有n*m个子过程，i＝1,2,...,n，i表示第i个子源。
[0135] 本实施例中，子源的破裂过程可拆分为破裂前锋初至时的快速破裂及后续的缓慢破裂过程。则第i个子源的滑动速率sr,i(t)可由m'个高斯模型叠加得到：
[0136]
[0137] 其中，Ci,j表示每个子函数的能量占比， a表示震源衰减率，Ai,j表示时间比例因子，τi,j和t0i,j分别表示第i个子源、第j个子函数的持续时间和开始时间；可用子源总上升时间τi表示：
[0138]
[0139]
[0140] 上式中，Ci,j表示每个子函数的能量占比， 及 均为系数，可通过拟合历史地震记录或拟合动力学滑动速率模型获得。最终，子源每个子函数的破裂前锋到达时间TiF,j和滑动量ui,j可由下式计算：
[0141]
[0142] ui,j＝Ci,jui
[0143] 其中，ui表示子源滑动量。
[0144] 本实施例中，如图4所示，其给出了采用两段高斯子函数拟合Liu et al.(2006)滑动速率模型的示意图，图中G1和G2为归一化的高斯模型。
[0145] S6、根据拆分结果，得到宽频带有限震源模型，其具体为：
[0146] 根据子过程上升时间分布、子过程的破裂前锋分布以及子过程的滑动量分布，得到宽频带有限震源模型。
[0147] 本实施例中，通过步骤S1到步骤S6获得宽频带有限震源模型后，可采用谱元法模拟震源破裂时的波动传播过程。
[0148] 下面对本发明作进一步的说明。
[0149] A、某次地震矩震级为Mw4.5，震源深度8km，震源总上升时间2.4s。考虑其震级较小，采用点源模型。模拟时，分别采用本方法提出的多段高斯子函数滑动速率模型和Liu et al.2006提出的滑动速率模型分别模拟其震源破裂过程，获得距离震源水平距离6.5km的某基岩台站的模拟结果。模拟和实际记录的水平向加速度反应谱RotD50值对比如图5所示。RotD50为两个水平向记录反应谱在水平向所有方位角分量的50％分位数。从图中可以看出，采用本方法提出的多段高斯子函数滑动速率模型获得的模拟地震动有更加丰富的高频成分，且与实际记录更加接近。
[0150] B、某次地震距震级为Mw 6.8，由于震级较大，采用上述第一步～第六步的有限震源面源模型模拟该地震破裂过程，获得距离震源水平距离～7km的河谷基岩台站模拟结果。模拟和实际记录的水平向地震动加速度反应谱RotD50值对比如图6所示，该台站三向地震动时程模拟和实际记录对比如图7所示。
[0151] 从图6中可以看出，采用本方法提出的基于多段高斯子函数滑动速率模型的宽频带有限震源模型能够模拟丰富的高频成分，且其加速度反应谱与实际记录较为接近。同时，从图7中可以看出，该方法能够模拟河谷中三方向地震动的差异，且三向地震动无论加速度、速度、位移均与实际地震记录形状相近。
