Effect of 3-D step topography on ground motion
-
摘要: 运用谱元法模拟了不同类型三维台阶地形对点源地震动的影响, 发现在中小地震中, 大型山系边缘的三维台阶地形对地表地震波具有一定的影响. 台阶斜坡与上水平地层交线和下水平地层交线的影响作用各不相同. 通过地表接收到的地震动描述了各种不同地形的波形快照特点, 并且通过绘制的台阶斜坡与上水平地层交线及下水平地层交线的测线图, 发现波形相对于水平地形有明显的后续波. 同时发现, 台阶地形上下测线相对于水平地形具有放大作用.下测线波形放大倍数基本在1.1—1.5, 上测线放大倍数在0.9—1之间. 通过描述山体边缘台阶结构在点源地震中产生的影响, 定量计算了三维台阶不同角度各个位置相对于水平地形的放大倍数, 为以后工程研究提供一个依据.Abstract: This paper simulates the different types of three-dimensional step topography effects on ground motion due to a point source by using spectral element method. Level topography and several different 3-D step topographies are simulated, which shows that the 3-D step topographies have an important impact on seismic waves. The effect from the intersecting lines of step slopes with upper level is different from that with the lower level. Moreover, the characteristics of various topographies are described using snapshots of the ground motions; the measuring line diagrams of the intersecting lines of step slope with upper level and those with the lower level demonstrate that the waveform diagrams of 3-D step topography have obvious scattering waves relative to level topography. Also it shows that the values on the upper and lower measuring lines on 3-D step topography are amplified relative to level topography. On the waveforms of lower measuring lines the amplification ranges from 1.1 to 1.5, and that is from 0.9 to 1 on the waveforms of upper measuring lines. These studies describe the impact of step topography at the edge of the mountain due to a point-source earthquake and give the respective positions’magnification of the 3-D step topography relative to level topography, providing a basis for forward engineering studies.
-
引言
正演模拟是研究地震波在地球介质中传播规律的重要工具．目前波动方程类的数值模拟方法主要有有限差分法、虚谱法和有限元法等(董渊等，2003)．有限差分法(裴正林，2004)的主要优点是计算速度快，占用内存小；缺点是会产生数值频散，而且传统的规则网格有限差分法精度较低．交错网格有限差分法(Virieux，1986；董良国等，2000)则可以大大减少数值频散、提高计算精度，是实际应用中常用的一种数值模拟方法．
地球介质中广泛存在着各向异性，横向各向同性(transversely isotropic media，简写为TI)介质是最普遍的一种，包括VTI(vertical transversely isotropic media)，HTI(horizon taltransversely isotropic media)和TTI(tilted transversely isotropic media)介质．基于各向同性假设的正演模拟方法不足以准确地描述实际地球介质中地震波的传播．在以往的研究中，考虑各向异性的弹性波动方程正演模拟取得了一定的进展(裴正林，2004)．该方法是用三分量矢量场来描述弹性波波场，其正演结果能够获得比标量波场更丰富的波场信息．但其存在一定的缺陷(何兵寿，张会星，2006)：①各向异性参数复杂，物理意义不明确，在实际生产中难以获得；②多分量数据占用资源大，计算效率低；③由于多分量地震勘探成本高，目前野外采集仍然以纵波地震资料为主，多分量弹性波理论的应用受到资料不足的限制．
为克服各向异性弹性波动方程在实际应用中的局限，Alkhalifah(2000)从VTI介质的精确qP-qSV波频散关系出发，假设垂向剪切波速度vS0=0，推导了VTI介质四阶拟声波方程．为降低四阶偏微分方程计算的复杂性，Zhou等(2006)和Du等(2008)从声学近似的频散关系出发，通过引入不同的辅助波场函数，分别推导了两种VTI介质二阶耦合拟声波方程．为结合交错网格技术实现精确的有限差分数值模拟，Hestholm(2009)推导了一阶偏微分形式的VTI介质拟声波方程．Duveneck等(2008)从声学近似的胡克定律出发，也推导了二阶形式和一阶应力-速度形式的VTI介质拟声波方程．与Alkhalifah(2000)方程一样，几种拟声波方程均存在低速度、低振幅的qSV人为干扰波问题(Grechka et al，2004)，而这种干扰波的存在会影响正演模拟和偏移成像的最终结果．
本文以VTI介质为例，首先从声学近似的胡克定律和频散关系两个不同角度出发，研究了3种不同形式的VTI介质一阶拟声波方程，并引入旋转坐标系，将其推广到TTI介质中；然后采用交错网格有限差分方法求解其关于时间和空间的一阶偏微分算子，实现了各向异性拟声波数值模拟，并通过数值结果的对比分析验证了3种一阶拟声波方程在运动学和动力学上的等价性；最后应用VTI介质拟声波逆时偏移算法实现了各向异性HESS模型的准确成像．
1. 各向异性介质一阶拟声波方程推导
1.1 从声学近似的胡克定律角度
首先对VTI介质胡克定律进行经典的声假设近似(Alkhalifah，2000).令垂向剪切波速度vS0=0，得到用Thomsen参数表示的VTI介质声假设近似胡克定律：
式中：σij和εkl分别为应力张量和应变张量(下标1，2，3分别对应x，y，z方向)；ρ为密度；vP0为垂向纵波速度；ε和δ为Thomsen各向异性参数(Thomsen，1986)．
波动方程1.式(1)两端对时间t求导，并结合运动方程ρ∂vi/∂t=∂σij/∂xj(下标ij=1，2，3分别对应x，y，z方向)，可得到一阶应力-速度形式的VTI介质拟声波方程(Duvenecket al，2008)：
式中：u，v，w分别为沿x，y，z方向的速度分量；p=σ11=σ22和q=σ33分别为水平应力分量和垂直应力分量；vPx为对称平面内的qP波速度，v2Px=v2P0(1+2ε)；vPn为qP波的正常时差速度，v2Pn=v2P0(1+2δ)．
1.2 从声学近似的qP-qSV波频散关系角度
从声假设近似的三维VTI介质频散关系出发，
式中：w为角频率；kx，ky，kz为空间波数；vP0，vPx，vPn的物理意义同式(2)．由于式(3)中相应的时空域四阶偏微分方程中存在四阶空间和时空混合偏导，不易直接求解，Alkhalifah(2000)引入辅助中间变量p(x，y，z，t)，令p=∂F2/∂t2，将其转换为二阶偏微分形式的波动方程:
基于同样的声近似频散关系，Du等(2008)引入新的辅助波场，得到二阶耦合形式的VTI介质拟声波方程:
式中，p=p(x，y，z，t)为伪应力波场，q=q(x，y，z，t)是为简化计算而引入的辅助波场．
波动方程2.Hestholm(2009)在式(4)的基础上，通过引入3个中间变量κ，ψ，ζ，得到一阶偏微分形式的VTI介质拟声波方程：
式中：u，v，w分别为x，y，z方向上的质点振动速度；η=(ε-δ)/(1+2δ)为非椭圆参数．
波动方程3.本文在式(5)的基础上，引入波场p和q的伪速度分量up，vp，wp和uq，vq，wq，并分别沿x，y，z3个方向进行分解，得到新的VTI介质一阶应力-速度拟声波方程：
式中，up，vp，wp和uq，vq，wq分别表示波场p和q在x，y，z方向的伪速度分量．
1.3 从VTI介质推广至TTI介质
由VTI介质推广至TTI介质只是增加了代数复杂性，两者在物理意义上并没有本质的差异．因此，可以通过引入对称轴倾角θ和方位角φ，将垂直坐标系X(x，y，z)旋转至倾斜坐标系X′(x′，y′，z′)中．引入三维坐标系的旋转关系(Cheng，Kang，2014)：
式中，旋转矩阵A依赖于倾角θ和方位角φ，即
由此可得旋转坐标系与垂直坐标系下的一阶偏导算子关系为
式中AT为旋转矩阵A的转置．将VTI介质一阶拟声波方程中相应的偏导算子通过坐标旋转后即可得到TTI介质形式．
2. VTI介质一阶拟声波方程的数值解法
2.1 交错网格高阶有限差分格式的构造
交错网格技术的基本思想是计算两个相邻网格点之间的偏导数(董良国等，2000)．本文以式(7)的二维形式为例，利用交错网格对其应力和速度分量进行差分离散(图 1)，得到时间二阶和空间高阶的交错网格有限差分格式(式(2)和式(6)的差分格式同理可得)：
式中：Δt为时间步长；Δx和Δz分别为x和z方向上的空间步长；i，j为空间离散点号；k为时间离散点号；pk+1/2i，j表示(k+1/2)时刻的pi，j值；2N为差分方程的精度；C(N)n为交错网格一阶差分系数．
![]()
图 1 交错网格及波场分量位置示意图Figure 1. Schematic diagram of staggered-grids and the location of wavefield2.2 完全匹配层边界条件
完全匹配层(perfectlymatchedlayer，简写为PML)边界条件是Berenger(1994)提出的一种高效吸收边界条件，被认为是目前吸收效果最好的边界条件而得到广泛应用．本文以式(7)的二维形式为例，依据PML方程分裂思路，推导相应的时间域PML控制方程如下：
式中：d(x)，d(z)为x和z方向的阻尼因子；px，pz和qx，qz分别为波场p，q在x和z方向上的分裂算子．
3. qSV人为干扰波噪音及压制方法
通过数值求解各向异性介质拟声波方程，可以比较准确地描述各向异性介质中qP波的运动学特征．然而在qP波传播的同时，其内部还存在一种菱形状的qSV人为干扰波(Grechkaet al，2004)．它的存在不仅会产生噪声问题，而且由于其速度值相对较低，往往还会引起较为严重的数值频散，进而影响正演模拟结果和降低偏移成像精度，因此必须对其进行压制或消除．
3.1 人为干扰波的产生机制
为求解各向异性介质拟声波方程，利用平面波解F(x，z，t)=A(t)exp[i(kxx+kzz)]作为试验解，将其代入式(4)的二维形式可得到两组复数解：
其中
分析可知，该两组复数解按指数增加或衰减的关键取决于指数项的符号．Alkhalifah(2000)研究发现，在实际介质中a1恒为负数，此时式(13)中F1的解对应qP波，其描述的是qP波的传播．当a2为正值，即η=(ε-δ)/(1+2δ)＜0时，式(13)中F2的解呈指数增长，这会导致数值结果不稳定．而当a2为负值，即η＞ 0时，F2对应的解会产生一组“人造的”地震波，其波速小于qP波，实际介质中这种人为干扰波并不存在．其产生的根源是基于声近似的拟声波方程并不是纯qP波方程，虽然将沿对称轴方向的剪切波相速度设为零，但在其它传播方向上的值并不为零，更为重要的是qSV波的群速度值也是非零的．
3.2 人为干扰波的压制方法
由人为干扰波的产生机制可知，当η=0(即ε=δ)时，式(4)退化为各向同性介质纵波方程，此时式(13)中F2的解变为一个与时间无关的函数，即其对应的qSV波不会在各向同性或椭圆各向异性(ε=δ)介质中传播．因此，在激发震源的周围设置一个环状的各向同性或椭圆各向异性盒子，即震源环(Duvenecket al，2008)，则可以有效地压制人为干扰波的传播．此外，也可以采用投影滤波方法来进一步压制人为干扰波(Zhang，Zhang，2009)．Xu和Zhou(2014)通过引入一个辅助标量算子S，提出一种纯的拟声波方程，以更好地解决复杂各向异性介质中产生的qSV人为干扰波问题．
4. 模型试算
4.1 各向异性均匀介质模型
采用空间十阶和时间二阶的交错网格有限差分方法对3种一阶拟声波方程式(2)、(6)、(7)进行数值模拟．均匀VTI介质模型参数为vP0=3000m/s，ε=0.24，δ=0.10，ρ=1g/cm3．网格大小为301×301，网格尺寸为Δx=Δz=10m，时间采样间隔Δt=1ms，雷克子波震源主频为30Hz，位于模型中心处．
图 2为VTI介质弹性波方程与3种一阶拟声波方程正演模拟得到的波场快照结果．由图 2b，c，d可知，外层传播较快的波是本文主要研究的qP波，且3种拟声波方程正演模拟均存在qSV人为干扰波(内层传播较慢的菱形状波形)．与图 2a对比可知，3种一阶拟声波方程均能有效地描述VTI介质中qP波的运动学特征，因此三者在运动学上是等价的．需注意，本文后面给出的各种波动方程的波场快照和地震记录均表示与图 2相同的分量．
![]()
图 2 波场快照对比(t=400 ms)(a)VTI介质弹性波方程正演得到的垂直应力分量波场；(b)用式(2)正演得到的波场q；(c)用式(6)正演得到的波场p；(d)用式(7)正演得到的波场pFigure 2. Comparison of wavefield snapshots (t=400 ms)(a)Wavefield of vertical stress component modeled by elastic wave equation in VTI media;(b)Wavefield q modeled by equation(2);(c)Wavefield p modeled by equation(6);(d) Wavefield p modeled by equation(7)图 3为应用公式(12)中的PML边界条件对qP波(图 3a)和qSV波(图 3b)的吸收效果．可以看出，PML边界条件不仅可以有效地衰减qP波，对qSV波也有较好的吸收效果，并未因为低速度的人为干扰波而在模型边界处引起较大的数值频散．
![]()
图 3 应用PML边界条件后的波场快照(a)qP波，t=600 ms；(b)qSV波，t=2 000 msFigure 3. Wavefield snapshots with PML boundary condition(a)qP wave，t=600 ms；(b)qSV wave, t=2 000 ms图 4a，b给出了设置震源环前后式(7)所对应的正演模拟结果．在震源附近设置半径为30m的小区域震源环，可以有效地压制拟声波正演模拟中产生的qSV人为干扰波．同时，利用旋转坐标系，可以得到TTI拟声波正演模拟对应的波场快照(图 4c，d)．
![]()
图 4 设置震源环前后各向异性介质拟声波正演的波场快照对比(t=400 ms)(a)VTI介质(设置震源环前)；(b)VTI介质(设置震源环后)；(c)TTI介质(设置震源环前，θ=45°)；(d)TTI介质(设置震源环后，θ=45°)Figure 4. Wavefield snapshots comparison of anisotropic pseudo-acoustic wave modeling without and with source box (t=400 ms)(a)VTI media (without source box)；(b)VTI media (with source box)；(c)TTI media (without source box, θ=45°)；(d)TTI media (with source box, θ=45°)4.2 TTI介质层状模型
4.2.1 波场特征分析
图 5给出了TTI介质层状模型及其各向异性参数．采用3种一阶拟声波方程进行正演模拟，其观测系统如下：模型网格大小为301×151，网格间距为10m×10m；选用震源主频为30Hz的雷克子波，并将炮点置于模型中(1500m，10m)位置处；采用中间放炮、两边接收的方式，共301道，道间距为10m；时间采样间隔Δt=1ms，记录时间为1.5s．
图 6分别显示了3种VTI一阶拟声波方法(图 6a-f)、TTI一阶拟声波方法(图 6g)、VTI弹性波方法(图 6h)和各向同性声波方法(图 6i)正演模拟得到的合成地震记录．从图 6a，b，c中地面位置为1200m处抽取一个近偏移距地震道记录，得到如图 7a所示的单道地震信号对比图．同样，从图 6d，e，f中相应位置处抽取一个地震道，得到的结果如图 7b所示．需注意，图 7中所有的单道地震信号对比图均是切除直达波后的结果，其振幅值表示相对大小，并且只取400—1200ms的旅行时范围，主要研究两个反射层位的有效信号．
![]()
图 6 单炮地震记录对比(a,d)用式(2)正演得到的炮记录；(b,e)用式(6)正演得到的炮记录；(c,f)用式(7)正演得到的炮记录；(g)TTI介质拟声波方程正演得到的炮记录；(h)VTI介质弹性波方程正演得到的炮记录；(i)各向同性声波方程正演得到的炮记录.其中，(a)-(c)为设置震源环前的结果；(d)-(f)为设置震源环后的结果Figure 6. Comparison of single shot seismic records(a,d) Shot records modeled by equation (2)；(b,e) Shot records modeled by equation (6)；(c,f) Shot records modeled by equation (7);(g) Shot records modeled by TTI pseudo-acoustic wave equation；(h)Shot records modeled by VTI elastic wave equation；(i) Shot records modeled by isotropic acoustic wave equation. Figs.(a)-(c) are results without source box, and Figs.(d)-(f) are results with source box![]()
图 7 单道地震信号对比(a)从图6a，b，c中近偏移距(地面位置为1 200 m)处抽取的单道信号对比；(b,c)从图6d e，f中近偏移距(地面位置为1 200 m)和远偏移距(地面位置为200 m)处抽取的单道信号对比；(d,e)从图6f，h中近偏移距(地面位置为1 200 m)和远偏移距(地面位置为200 m)处抽取的单道信号对比；(f,g) 从图6f，g，i中近偏移距(地面位置为1 200 m)和远偏移距(地面位置为200 m)处抽取的单道信号对比Figure 7. Comparison of single trace seismic signalsFig.(a) is single trace signals comparison of near offset (distance=1 200 m) extracted from Figs.6a,b,c; Fig.(b) and (c) are single trace signals comparison of near offset (distance=1 200 m) and far offset (distance=200 m) extracted from Figs.6d,e,f; Figs.(d) and (e) are single trace signals comparison between near offset (distance=1 200 m) and far offset (distance=200 m) extracted from Figs.6f, h; Figs.(f) and (g) are single trace signals comparison between near offset(distance=1 200 m) and far offset (distance=200 m) extracted from Figs.6f, g, i图 7a中由于低速度的qSV人为干扰波造成了数值频散，设置震源环后，人为干扰波得到了很好地压制，且几乎不会对qP波产生影响，如图 7b所示．由图 7b还可以看出，3种一阶拟声波正演模拟方法具有很好的一致性，各层同相轴、旅行时和振幅均吻合较好，三者在运动学和动力学上是等价的．类似地，图 7c是从图 6d,e,f中远偏移距(地面位置为200m)处抽取的单道信号对比结果，也能得到相同的结论．
图 7d是从图 6f,h中近偏移距(地面位置为1200m)处抽取的单道地震记录．对比图 7d中蓝线与绿线可知，在两个反射层上，VTI介质一阶拟声波方程与弹性波方程的qP波反射旅行时吻合较好，但振幅值有所差异．因此，本文推导的一阶拟声波方程能够比较准确地模拟VTI介质中qP波的运动学(旅行时和速度等)特征，而在动力学(振幅等)特征上也是VTI弹性波方程的有效近似．类似地，图 7e中远偏移距(地面位置为200m)处抽取的地震道信号对比结果也能得到相同的结论．
图 7f,g分别是从图 6f,g,i中近偏移距(地面位置为1200m)和远偏移距(地面位置为200m)处抽取的单道地震记录．通过对比可知，由于考虑各向异性介质(VTI和TTI)的影响，各层的反射旅行时和振幅均会发生变化，且往往造成旅行时和振幅减小，特别是远偏移距(地面位置为200m)处，这种影响更明显．同时可以看出，模型第二层中受到TTI介质的倾角影响后，振幅和旅行时会进一步发生改变．因此，在实际应用中，对大倾角构造进行处理时必须要考虑到TI介质的倾角因素．
4.2.2 计算效率分析
在波动方程类的正演模拟和逆时偏移中，数值计算效率也是值得关注的．本文将3种VTI一阶拟声波方程、各向同性声波方程及VTI弹性波方程的计算效率进行了对比，其中偏微分算子的有限差分实现占据了波动方程类数值模拟的主要计算量．表 1展示了二维情形下几种波动方程各自的计算量，以及对于同一层状模型(图 5)正演模拟所需的平均计算时间(取10次结果的平均)．虽然计算机硬件不同可能会导致计算时间的差异，但是在同等条件下得到的计算时间还是具有一定的参考价值．可以看出：VTI一阶拟声波方程(2)的计算时间最少，与各向同性声波方程相当；而另外两种VTI一阶拟声波方程(6)和(7)的计算时间相对多一些，但仍然少于VTI弹性波方程的计算时间．因此，利用VTI介质拟声波方程不仅能减少计算所需的内存，还能有效地提高计算效率．
表 1 几种波动方程的计算效率对比Table 1. Computational efficiency comparison of the several wave equations方程类型 计算量 平均计算时间/s VTI一阶拟声波方程(2) ∂p/∂x,∂q/∂z,∂u/∂x,∂w/∂z 72 VTI一阶拟声波方程(6) ∂p/∂x,∂p/∂z,∂u/∂x,∂k/∂z,∂w/∂z,∂ζ/∂z 105 VTI一阶拟声波方程(7) ∂p/∂x,∂p/∂z,∂q/∂x,∂q/∂z,∂up/∂x,∂wq/∂z 103 各向同性声波方程 ∂p/∂x,∂p/∂z,∂u/∂x,∂w/∂z 70 VTI弹性波方程 ∂τxx/∂x,∂τxz/∂z,∂τxz/∂x,∂τzz/∂z,∂vx/∂x,∂vx/∂z,
∂vz/∂x,∂vz/∂z130 4.3 VTI介质HESS模型
为进一步验证本文方法的准确性，首先用本文的VTI一阶拟声波方程正演模拟方法得到HESS模型的合成地震记录，然后对其进行逆时偏移成像．HESS模型的网格大小为601×501，网格间距为10 m×10 m，其各向异性参数如图 8所示．正演模拟的观测系统为：选用震源主频为30 Hz的雷克子波，采用中间放炮、两边接收的方式；第一个炮点位于模型的(1 500 m，10 m)位置处，总共31炮，每炮301道接收；炮间距为100 m，道间距为10 m；时间采样间隔Δt=1 ms，总采样时间为4 s．分析对比两种偏移结果可知：由于各向异性因素的影响，常规的各向同性声波逆时偏移方法成像效果较差，特别是陡倾角断层及各向异性较强的区域(图 9a)；而利用VTI介质拟声波逆时偏移方法则能够准确地刻画地下各向异性构造和陡倾角断层，如图 9中红色方框所示区域．
![]()
图 9 各向同性(a)和各向异性(b)介质逆时偏移剖面Figure 9. Fig.9 Reverse-time migration sections for isotropic (a) and anisotropic (b) media5. 讨论与结论
从声学近似的VTI介质胡克定律和频散关系两种思路出发，对比分析了3种不同形式的VTI介质一阶拟声波方程，并利用交错网格高阶有限差分法实现了VTI介质一阶拟声波数值模拟，然后通过旋转坐标系将其推广到TTI介质中．相比于四阶和二阶形式，一阶拟声波方程形式简单，易于结合交错网格技术得到精确的正演模拟结果.波场快照验证了本文推导的PML边界条件能够达到理想的边界反射吸收效果．最后，从VTI介质拟声波方程的解析解出发，简单讨论了其固有的qSV人为干扰波的产生机制，并通过设置震源环的方法实现了对干扰波噪音的有效压制．
数值结果表明，3种一阶拟声波方程在运动学和动力学上是等价的，它们均能比较准确地模拟各向异性介质中qP波的运动学(旅行时和速度等)特征，这便于研究qP波的传播规律．结合计算效率分析可知，相对于各向异性弹性波正演模拟，拟声波正演模拟不仅节省了计算内存、效率更高，同时克服了前者的弹性参数意义不明确、波场分离困难等缺陷．设置小区域的震源环能够有效地压制人为干扰波，且不会对有效qP波造成较大影响．由于各向异性因素的影响，地震记录的波场特征(旅行时和振幅等)会发生一定变化，特别是远偏移距处这种影响更为明显，进而会影响偏移和反演的精度．因此，在实际大偏移距、宽方位地震资料的应用中，各向异性因素不能忽略．
本文采用震源环方法虽能从震源处实现对人为干扰波的有效压制，但其难以处理波传播过程中散射源新产生的干扰噪音．因此，人为干扰波的压制方法还有待完善，各向异性介质中纯qP波模式的波动方程及其正演模拟方法也将是今后的研究重点．
衷心感谢审稿专家和编辑部对本文提出的宝贵意见和建议；感谢中国石油大学(华东)Leon课题组的帮助．
-
![]()
图 1 三类地形模型 (a) 水平地形； (b) 简单的二维阶梯地形； (c) 有尖凸的阶梯地形； (d) 有尖凹的阶梯地形； (e) 有圆凸的阶梯地形； (f) 中间有圆凹的阶梯地形
Figure 1. The three types of terrain models (a) The level terrain； (b) The terrain with a simple step； (c) The step terrain with a sharp convex； (d) The step terrain with a sharp concave； (e) The step terrain with a rounded convex; (f) The step terrain with a rounded concave
![]()
图 2 对应图1(a)—(f) 6种地形的地表x分量速度7—9 s快照
Figure 2. Snapshots of x-component velocity waveforms on the surface position corresponding to the six types of terrains in Figs.1(a)-(f) ranging from 7 s to 9 s
![]()
图 3 台阶下测线 蓝色、 绿色、 黄色、 红色、 黑色分别 对应图1中(b)—(f)的地形
Figure 3. The lower lines of steps Blue, green, yellow, red, black lines correspond to the terrains in Figs.1(b)-(f), respectively
![]()
图 4 不同角度台阶下测线 绿色、 红色分别对应图1中(c)、 (d)不同角度的地形
Figure 4. The lower lines of steps with different angles Green and red lines correspond to the terrains with different angles in Figs.1(c) and 1(d), respectively
![]()
图 5 各个模型地表接收到的狓方向速度峰值图 (a)水平地形;(b)夹角为60°的尖凸地形;(c)夹角为90°的尖凸地形;(d)夹角为120°的尖凸地形; (e)简单台阶地形;(f)夹角为120°的尖凹地形;(g)夹角为90°的尖凹地形;(h)夹角为60°的尖凹地形
Figure 5. The x-componet peak velocities observed on the surface of each model (a)Thelevel terrain;(b)The terrain with 60° sharp convex;(c) Theterrain with 90° sharp convex;(d) The terrain with 120° sharpconvex;(e) Simple two-dimensional step terrain;(f) The terrain with 120°sharp concave;(g) The terrain with 90° sharp concave;(h) The terrain with 60° sharp concave
![]()
图 6 台阶端线x分量波形图(左侧图为上测线波形， 右侧图为下测线波形) (a) 简单二维台阶地形波形； (b) 尖凸台阶地形波形； (c) 尖凹台阶地形波形
Figure 6. The x-component velocity waveforms (the left panels are waveforms measured on upper line， and the right ones are measured on lower line) (a) The diagram of the simple two-dimensional step terrain； (b) The diagram of the step terrain with sharp convex； (c) The diagram of the step terrain with the sharp concave
![]()
图 7 震源在上测线(左侧图)和下测线(右侧图)投影点的x分量速度波形图 图(a)， (b)， (c)分别为图6(a), (b), (c)中震源在上下测线投影点的速度波形
Figure 7. The x-component velocity waveformｓ of the projected point of the source on the upper (left panels) and lower (right panels) lines Figs.(a), (b), (c) separately correspond to the velocity waveforms projected on the point from the source to the survey line in Figs.6(a), (b), (c)
![]()
图 8 震源在上测线(左侧图)和下测线(右侧图)投影点的x分量加速度波形图 图(a)， (b)， (c)分别为图6(a)， (b)， (c)中震源在上下测线投影点的加速度波形
Figure 8. The x-component acceleration waveforms of the projected point of the source on the upper (left panels) and lower (right panels) lines Figs.(a), (b), (c) separately correspond to the velocity waveforms projected on the point from the source to the survey line in Figs.6(a), (b), (c)
![]()
图 9 表1和表2中测点位置
Figure 9. Distribution of the measuring points listed in Tables 1 and 2
表 1 三维台阶地形与水平地形速度峰值比较(放大倍数)
Table 1 Comparison of the peak velocities of three-dimensional step terrain with those of the horizontal terrain (magnification)
下载: 导出CSV
表 2 三维台阶地形与水平地形加速度峰值比较(放大倍数)
Table 2 Comparison of the peak accelerations of three-dimensional step terrain with those of the horizontal terrain (magnification)
下载: 导出CSV
-
薄景山, 李秀领, 李山有. 2003. 场地条件对地震动影响研究的若干进展[J]. 世界地震工程, 19(2): 11-15. Bo J S, Li X L, Li S Y. 2003. Some progress of study on the effect of site conditions on ground motion[J]. World Earthquake Engineering, 19(2): 11-15 (in Chinese).
陈丙午. 1982. 不规则地形对地震动及震害的影响[J]. 地震工程与工程振动, 2(1): 12-19. Chen B W. 1982. Effect of irregular topography on earthquake motion and earthquake damage[J]. Earthquake Engineering and Engineering Vibration, 2(1): 12-19 (in Chinese).
董良国. 2005. 复杂地表条件下地震波传播数值模拟[J]. 勘探地球物理进展, 28(3): 187-194. Dong G L. 2005. Numerlcal simulation of seismic wave propagation under complex near surface conditions[J]. Progress in Exploration Geophysics, 28(3): 187-194 (in Chinese).
胡元鑫, 刘新荣, 罗建华, 张梁, 葛华. 2011. 汶川震区地震动三维地形效应的谱元法模拟[J]. 兰州大学学报: 自然科学版, 47(4): 393-402. Hu Y X, Liu X R, Luo J H, Zhang L, Ge H. 2011. Simulation of three-dimensional topographic effects on seismic ground motion in Wenchuan earthquake region based upon the spectral-element method[J]. Journal of Lanzhou University: Natural Sciences, 47(4): 393-402 (in Chinese).
景立平, 卓旭炀, 王祥建. 2005. 复杂场地对地震波传播的影响[J]. 地震工程与工程振动, 25(6): 16-23. Jing L P, Zhuo X Y, Wang X J. 2005. Effect of complex site on seismic wave propagation[J]. Earthquake Engineering and Engineering Vibration, 25(6): 16-23 (in Chinese).
刘晶波. 1996. 局部不规则地形对地震地面运动的影响[J]. 地震学报, 18(2): 239-245. Liu J B.1996. The effect of local irregular topography on seismic ground motion[J]. Acta Seismologica Sinica, 18(2): 239-245.
袁晓铭, 廖振鹏. 1996. 任意圆弧形凸起地形对平面SH波的散射[J]. 地震工程与工程振动, 16(2): 1-12. Yuan X M, Liao Z P. 1996. Scattering of plane SH waves by a cylingrical hill of circular-arc cross-section[J]. Earthquake Engineering and Engineering Vibration, 16(2): 1-12 (in Chinese).
荣棉水, 李小军. 2007. 局部地形对出平面运动谱特性的影响分析[J]. 中国地震, 23(2): 147-156. Rong M S, Li X J. 2007. The effect analysis of topography on the spectrum property of anti-plane movement[J]. Earthquake Research in China, 23(2): 147-156 (in Chinese).
荣棉水, 李小军, 吕悦军, 尤红兵. 2009. 平台地形对地震地面运动特征周期值的影响[J]. 中国地震, 25(2): 178-185. Rong M S, Li X J, Lu Y J, You H B. 2009. Effect of protruding topography on characteristic period of ground motion[J]. Earthquake Research in China, 25(2): 178-185 (in Chinese).
王海云, 谢礼立. 2010. 自贡市西山公园地形对地震动的影响[J]. 地球物理学报, 53(7): 1631-1638. Wang H Y, Xie L L. 2010. Effects of topography on ground motion in the Xishan park, Zigong city[J]. Chinese Journal of Geophysics, 53(7): 1631-1638 (in Chinese).
王伟. 2011. 地震动的山体效应[D]. 哈尔滨: 中国地震局工程力学研究所: 1-173. Wang W. 2011. Effect of Hill Topography on Ground Motion[D]. Harbin: Institute of Engineering Mechanics, China Earthquake Administration: 1-173 (in Chinese).
周红. 2007. 两种地震波传播方法的研究[D]. 北京: 北京大学: 1-47. Zhou H. 2007. The Study of Two Kinds of Seismic Wave Propagation Methods[D]. Beijing: Peking University: 1-47 (in Chinese).
周红, 陈晓非. 2007. 凹陷地形对Rayleigh面波传播影响的研究[J]. 地球物理学报, 50(4): 1182-1189. Zhou H, Chen X F. 2007. A study on the effect of depressed topography on Rayleigh surface wave[J]. Chinese Journal of Geophysics, 50(4): 1182-1189 (in Chinese).
周红, 高孟潭, 愈言祥. 2010. SH波地形效应特征的研究[J]. 地球物理学进展, 25(3): 775-782. Zhou H, Gao M T, Yu Y X. 2010. A study of topographical effect on SH waves[J]. Progress in Geophysics, 25(3): 775-782 (in Chinese).
周正华, 张艳梅, 孙平善, 杨柏坡. 2003. 断层对震害影响的研究[J]. 自然灾害学报, 12(4): 20-24. Zhou Z H, Zhang Y M, Sun P S, Yang B P. 2003. Study on effect of fault on seismic damage[J]. Journal of Natural Disasters, 12(4): 20-24 (in Chinese).
Boore D M. 1972. A note on the effect of simple topography on seismic SH waves[J]. Bull Seismol Soc Am, 62(1): 275-284.
Jih R S, Mclaughlin K L, Der Z A. 1988. Free-boundary conditions of arbitrary topography in a two-dimensional explicit elastic finite-difference scheme[J]. Geohphysics, 53(8): 1045-1055.
Komatitsch D, Vilotte J P. 1998. The spectral element method: An efficient tool to simulation the seismic response of 2D and 3D geological structure[J]. Bull Seismol Soc Am, 88(2): 368-392.
Komititsch D, Tromp J. 1999. Introduction to the spectral element method for three-dimensional seismic wave propagation[J]. Geophys J Int, 139(3): 806-822.
Seriani G. 1998. 3-D large-scale wave propagation modeling by spectral element method on Gray T3E multiprocessor[J]. Comp Meth Appl Mech Eng, 164(1): 235-247.
Tessmer E, Kosloff D, Behle A. 1992. Elastic wave-propagation simulation in the presence of surface-topograghy[J]. Geophys J Int, 108(2): 621-632.
Tessmer E, Kosloff D. 1994. 3-D elastic modeling with surface topography by a Chebychev spectral method[J]. Geophysics, 59(3): 464-473.
Wang Y B, Takenaka H, Furumura T. 2001. Modeling seismic wave propagation in a two-dimensional cylindrical whole-earth model using pseudo-spectral method[J]. Geophys J Int, 145(3): 689-708.
-
期刊类型引用(4)
1. 陈利新，王胜雷，万效国，苏洲，马兵山. 哈拉哈塘地区共轭走滑断裂差异特征及演化. 西南石油大学学报(自然科学版). 2024(04): 19-37 . 
百度学术
2. 曾庆才，王清华，曾同生，陈胜，张凯. 基于测井修正的TTI伪弹性波逆时偏移优化方法. 西南石油大学学报(自然科学版). 2024(04): 38-50 . 百度学术
3. 秦宁. VTI介质角度叠加逆时偏移. 石油地球物理勘探. 2019(02): 341-347+237-238 . 百度学术
4. 李振春，杨富森，王小丹. 基于LS-RSGFD方法优化的横向各向同性(TI)介质一阶qP波高精度数值模拟. 地球物理学报. 2016(04): 1477-1490 . 百度学术
其他类型引用(6)
计量
- 文章访问数: 556
- HTML全文浏览量: 252
- PDF下载量: 22
- 被引次数: 10
