Ricker子波在非一致性地震合成记录中的应用以及在软土层中的地震响应.pdfVIP

2利用Ricker子波合成地震记录 利用Ricker子波，分别合成在弹性半空间内的岩层和软土层条件下的非一致性地面震 动记录[4|。岩层的横波速度为600m／s，泊松比为0．33，密度为2200kg／m3。软土层的横波 速度为200m／s，泊松比为0．33，密度为1800 kg／m3。为了减少其他外界条件的影响，假定 岩层和软土层的阻尼忽略不计。 2．1入射波动方程 拉梅一纳维微分方程描述了线性弹性体的位移场。在弹性体均质性假设初始条件以及无 体力作用下，此方程可以转化为边界积分方程(1)： f “ Z = “ Z≠ f ／L z r dfdr一 ／L z￡ t工、 f、， Zr drdr ／L ．|、 、， ／L ．J‘ 、， X *池 、， X k k *浓 巧 哟 ，¨川v ，¨川掣 式(1)中，￡盘和Ui。k表示t时刻在点X处荷载作用在边界点车处，得到的应力和位移 积分项，由附近边界点车的几何形状所决定。 对于入射波的散射波动方程(1)，不能直接应用到全部波场域，因为它不能满足发散条 件。能够应用在入射波和散射波场的叠加中，公式(1)可以用以下形式表示： 钒“引)=j’J．“孟(z㈠e，r)勺(z，-r)drdrx ‘0 (2) 一m孟(z^e，r)哟(z，r)drdlr。+娩(z，￡) 在公式(2)中，蹦ii为边界入射位移向量。一般来讲，在边界r上的各点位移和应力 是已知的，公式(2)可用于确定未知的边界值。对于数值计算，边界积分方程(2)在时 间、空间上离散后求解。 在边界条件下对公式(2)求积分只能得到函数在空间上的变化。对于特殊几何形状的 边界条件，这些函数不能被有效地积分。因此，边界11分解为常量、线性或二次等参数边 界单元。公式(2)可表示为以下形式： 5”) 钒“i(拿，tN)=∑∑lT{f—m圳(z，拿)dF／k(z)“5”)』)“i 121 m21} L N， u5f—m川(z，e)dnk(‘)￡5m’ (3) =∑∑I 扛lm-1置 公式(3)中tk和Uk分别为应力和位移核函数，可以求出基本解的时间积分。公式 以用矩阵符号表示如下： ∥￡Ⅳ=一“N+EN—uiN： 372 EN=芝：P∥一+1一泸∥一+1 ，(4) ：鬲 J。 通过与时间域的三维边界元素法计算结果的对比得到验证【5 2．2非一致性地面震动记录 地面位移以及对应的傅立叶频谱，Ricker子波的入射方向为90。，波前垂直于波的传播方向。 在本研究中考虑了入射方向分别为30。、60。和90。这3种情况的地震记录的合成。假定在一 有限元法和边界元法，计算观测点处地表的位移和加速度响应。 0．08 0．1 O．04 0 2