CN108957537B - 倾斜地震波作用下场地卓越周期的计算方法 - Google Patents

倾斜地震波作用下场地卓越周期的计算方法 Download PDF

Info

Publication number
CN108957537B
CN108957537B CN201810638788.7A CN201810638788A CN108957537B CN 108957537 B CN108957537 B CN 108957537B CN 201810638788 A CN201810638788 A CN 201810638788A CN 108957537 B CN108957537 B CN 108957537B
Authority
CN
China
Prior art keywords
wave
waves
incident
seismic
layer
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201810638788.7A
Other languages
English (en)
Other versions
CN108957537A (zh
Inventor
范留明
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Xian University of Technology
Original Assignee
Xian University of Technology
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Xian University of Technology filed Critical Xian University of Technology
Priority to CN201810638788.7A priority Critical patent/CN108957537B/zh
Publication of CN108957537A publication Critical patent/CN108957537A/zh
Application granted granted Critical
Publication of CN108957537B publication Critical patent/CN108957537B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/30—Analysis
    • G01V1/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开的倾斜地震波作用下场地卓越周期的计算方法，具体操作步骤如下：首先确定场地和地震波的各参数以及地震脉冲响应持时和时间采样间隔；计算入射至地面处产生的反射系数；然后，计算纵波和横波斜入射至各地层界面i处的反射系数和透射系数，再计算纵波和横波在各个土层界面的脉冲响应得到入射至地表面处的脉冲响应；最后，计算地面水平方向的脉冲响应时程并将其离散傅里叶变换获得脉冲地震动的傅里叶谱，从而确定场地卓越周期。本发明公开的倾斜地震波作用下场地卓越周期的计算方法，解决了现有计算方法通产将地震波看成垂直向上的平面剪切波，导致计算的场地卓越周期不准确的问题，该方法不受土层层数的限制，计算精度高且速度快。

Description

倾斜地震波作用下场地卓越周期的计算方法
技术领域
本发明属于岩土工程勘察技术领域，涉及一种倾斜地震波作用下场地卓 越周期的计算方法。
背景技术
当建筑结构的自振周期与场地自振周期接近或一致时，两者构成的系统 在地震过程中产生共振作用，由于共振作用而被放大的场地自振周期称为卓 越周期。为了避免地面建筑物与工程场地在地震时产生共振，首先需要确定 场地的卓越周期。在应用波速法计算场地的卓越周期时，假定覆盖土层及其 下卧基岩为力学性质沿竖向成层变化、沿横向均匀无限延伸的水平成层场 地，并将入射波看作是垂直向上的平面剪切波。事实上，地震波通常以某个 角度传播到地面，将地震波视为垂向传播不仅与实际情况明显不符，而且斜入射地震波对基岩土体及其地面建(构)筑物的作用与垂直入射地震波的情 形也明显不同，因此考虑地震波入射方式对场地卓越周期的影响是非常必要 的。
发明内容
本发明的目的是提供一种倾斜地震波作用下场地卓越周期的计算方法， 解决了现有计算方法通常将地震波看成垂直向上的平面剪切波，导致计算的 场地卓越周期不准确的问题，本发明的方法能够考虑地震波入射角度对卓越 周期的影响，为确定实际工程设计提供更全面的场地自振特性参数。
本发明所采用的技术方案是，倾斜地震波作用下场地卓越周期的计算方 法，具体操作步骤如下：
步骤1.根据岩土工程的勘查结果，确定场地土层数n、地层的密度ρi、 纵波速度vpi、横波速度vsi以及地层厚度hi，其中i＝1，2，......，n；
步骤2.根据需要计算的精度要求，确定地震脉冲响应持时tl和时间采 样间隔Δt，并计算采样点数nt；
步骤3.根据实际工程勘察资料选择地震横波入射角θ，计算地震波的射 线参数pw；
步骤4.分别计算地震纵波和横波在场地各土层中的传播时间以及其入 射至地面处产生的反射系数R1 pp+、R1 ps+、R1 sp+和R1 ss+；
步骤5.分别计算地震纵波和横波斜入射至各地层界面i处的反射系数
Figure BDA0001702104110000021
和透射系数
Figure BDA0001702104110000022
步骤6.根据步骤5得到的地震纵波和横波斜入射至各地层界面i处的 反射系数
Figure BDA0001702104110000023
和透射系数
Figure BDA0001702104110000024
分别计算纵波和横波在各个土 层界面的脉冲响应
Figure BDA0001702104110000025
然后得到纵波和横波入射至地表面处的脉冲 响应
Figure BDA0001702104110000026
结合步骤4得到的结果合成地面水平方向的脉冲响应时程
Figure BDA0001702104110000027
步骤7.将步骤6得到脉冲响应时程
Figure BDA0001702104110000028
通过进行离散Fourier变换(DFT) 获得脉冲地震动的Fourier谱F(u)；
步骤8.通过分析步骤7得到的脉冲地震动的Fourier谱F(u)，确定场地 卓越周期，即在Fourier谱中出现的若干个不连续尖峰状谱中，选取尖峰极 大值点所对应的周期作为场地卓越周期，其中，最大卓越周期为基本卓越周 期。
本发明的其他特点还在于，
步骤2中地震脉冲响应持时tl和时间采样间隔Δt的确定依据为：Δt< 0.1tsmin，tl>10 tsmax，其中，tsmin＝min{nts1，nts2，......，ntsi，......}，tsmax＝nts1+ nts2+......+ntsi+......，ntsi＝hi/vsi/Δt；
采样点数nt的计算方法为：nt＝tl/Δt。
步骤3中地震波射线参数pw计算方法为：pw＝sinθ/vsn，且pwvpmax<1， 其中，vpmax＝max{vp1，vp2，......，vpi，......}。
步骤4中地震纵波和横波在工程场地各土层中的传播时间通过公式1计 算得到：
Figure BDA0001702104110000031
其中，hi是地层i的厚度，vpi、vsi分别是地震波在地层i中的纵波速度、 横波速度；
步骤4中地震纵波和横波斜入射至地面处产生的反射系数R1 pp+、R1 ps+、 R1 sp+和R1 ss+分别通过公式2和3计算得到：
Figure BDA0001702104110000032
Figure BDA0001702104110000033
其中，R1 pp+、R1 ps+分别是纵波入射至地面产生的同类波，记为P波、转 换波，记为S波的反射系数；R1 sp+、R1 ss+分别是横波入射至自由面处产生的P波、S波的反射系数，上标“pp”表示入射P波产生的同类波，“ps”表示 入射P波产生的转换波，“sp”表示入射S波产生的转换波，“ss”表示入射S 波产生的同类波；“+”表示地震波从下向上入射，“-”表示地震波从上向下 入射；vp1、vs1分别是地层1的纵波速度和横波速度；α1、β1分别是纵波、 横波传播至地面处的入射角和反射角。
步骤5中地震纵波从第i层入射至第i-1层时在界面i处的反射系数
Figure BDA0001702104110000044
Figure BDA0001702104110000045
和透射系数
Figure BDA0001702104110000046
通过公式4计算得到：
Figure BDA0001702104110000041
步骤5中地震横波从第i层入射至第i-1层时在界面i处的反射系数
Figure BDA0001702104110000047
Figure BDA0001702104110000048
和透射系数
Figure BDA0001702104110000049
通过公式5计算得到：
Figure BDA0001702104110000042
地震纵波从第i-1层入射至第i层时在界面i处的反射系数
Figure BDA00017021041100000410
和 透射系数
Figure BDA00017021041100000411
通过公式6计算得到：
Figure BDA0001702104110000043
地震横波从第i-1层入射至第i层时在界面i处的反射系数
Figure BDA00017021041100000412
和透射系数
Figure BDA00017021041100000511
通过公式7计算得到：
Figure BDA0001702104110000051
公式4-7中，i是地层i-1与地层i之间的地层界面；vpi、vsi分别是地层 i的纵波速度、横波速度，vpi-1、vsi-1分别是地层i-1的纵波速度、横波速度； αi-1、αi分别是纵波在界面i处的入射角或反射角、透射角；βi-1、βi是横波在 界面i处的入射角或反射角、透射角；i＝2，3，......，n。
步骤6中纵波和横波在各个土层界面的脉冲响应
Figure BDA0001702104110000052
和
Figure BDA0001702104110000053
通过公式8 计算得到：
Figure BDA0001702104110000054
其中，
Figure BDA0001702104110000055
和
Figure BDA0001702104110000056
分别表示从第i层入射至第i-1层时界面i处的纵波和横 波的响应；
Figure BDA0001702104110000057
和
Figure BDA0001702104110000058
分别是从第i-1层入射至第i层时界面i处纵波和横 波的响应；i＝2，......，n-1；j＝1，2，......，nt；
当i＝1时，纵波和横波入射至地表面处的脉冲响应
Figure BDA0001702104110000059
和
Figure BDA00017021041100000510
如公式9所 示：
Figure BDA0001702104110000061
其中，
Figure BDA0001702104110000062
和
Figure BDA0001702104110000063
分别是从第1层入射至地表面处的纵波和横波；
Figure BDA0001702104110000064
和
Figure BDA0001702104110000065
分别是从空气介质入射至地表面处的纵波和横波，j＝1，2，......，nt；
当i＝n时，公式8表示为公式10：
Figure BDA0001702104110000066
其中，
Figure BDA0001702104110000067
和
Figure BDA0001702104110000068
分别是从第n层入射至第n-1层时界面n处的纵波和 横波响应；
Figure BDA0001702104110000069
和
Figure BDA00017021041100000610
分别是从第n-1层入射至第n层时界面n处的纵波和 横波响应；j＝1，2，......，nt。
步骤6中合成地面水平方向的脉冲响应时程
Figure BDA00017021041100000612
通过公式11得到：
Figure BDA00017021041100000611
其中，j＝1，2，...，nt。
本发明的有益效果是，首次采用定量计算方法揭示了场地卓越周期与入 射方式的关系，由该方法得到的卓越周期不仅包含了周期信息，而且还包含 了入射方式和入射波型方面的信息，可为实际工程的抗震设计提供更全面的 场地自振特性参数，而由传统方法测定的卓越周期只是SV波垂直入射时的 一个计算特例。
附图说明
图1是本发明的一种倾斜地震波作用下场地卓越周期的计算方法的计算 流程图；
图2是本发明的一种倾斜地震波作用下场地卓越周期的计算方法中场地 反应模型；
图3是本发明的一种倾斜地震波作用下场地卓越周期的计算方法中倾斜 入射波作用下的地面反射，(a)P波入射，(b)SV波入射；
图4是本发明的一种倾斜地震波作用下场地卓越周期的计算方法中地震 波从第i层入射至第i-1层时在界面i处的反射波与透射波的示意图，(a)P 波入射，(b)SV波入射；
图5是本发明的一种倾斜地震波作用下场地卓越周期的计算方法中地震 波从第i-1层入射至第i层时在界面i处的反射波与透射波的示意图，(a)P 波入射，(b)SV波入射；
图6是本发明的一种倾斜地震波作用下场地卓越周期的计算方法中地震 波入射至地面的P波时程
Figure BDA0001702104110000071
图7是本发明的一种倾斜地震波作用下场地卓越周期的计算方法中地震 波入射至地面的S波时程
Figure BDA0001702104110000072
图8是本发明的一种倾斜地震波作用下场地卓越周期的计算方法中地震 波脉冲响应水平分量时程
Figure BDA0001702104110000073
图9是本发明的一种倾斜地震波作用下场地卓越周期的计算方法中地震 脉冲响应水平分量时程
Figure BDA0001702104110000074
的Fourier谱。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明的倾斜地震波作用下场地卓越周期的计算方法，具体操作步骤如 下：
步骤1.根据岩土工程的勘查结果，确定场地土层数n、地层的密度ρi、 纵波速度vpi、横波速度vsi以及地层厚度hi，其中，i＝1，2，......，n；
步骤2.根据需要计算的精度要求，确定地震脉冲响应持时tl和时间采 样间隔Δt，并计算采样点数nt；
步骤3.依据实际工程需要选择地震波的横波入射角θ，计算地震波的射 线参数pw；
步骤4.分别计算地震纵波和横波在工程场地各土层中的传播时间以及 地震纵波和横波斜入射至地面处产生的反射系数R1 pp+、R1 ps+、R1 sp+和R1 ss+；
步骤5.分别计算地震纵波和横波斜入射至各地层界面i处的反射系数
Figure BDA0001702104110000087
和透射系数
Figure BDA0001702104110000088
步骤6.根据步骤5得到的地震纵波和横波斜入射至各地层界面i处的 反射系数
Figure BDA0001702104110000085
和透射系数
Figure BDA0001702104110000086
分别计算纵波和横波在各个土 层界面的脉冲响应
Figure BDA0001702104110000081
然后得到纵波和横波入射至地表面处的脉冲 响应
Figure BDA0001702104110000082
结合步骤4得到的合成地面水平方向的脉冲响应时程
Figure BDA0001702104110000083
步骤7.将步骤6得到脉冲响应时程
Figure BDA0001702104110000084
通过进行离散Fourier变换(DFT) 获得脉冲地震动的Fourier谱F(u)；
步骤8.通过分析步骤7得到的脉冲地震动的Fourier谱F(u)，确定场地 卓越周期，即在Fourier谱中出现的若干个不连续尖峰状谱中，选取尖峰极 大值点所对应的周期作为场地卓越周期，其中最大卓越周期为基本卓越周 期。
步骤2中地震脉冲响应持时tl和时间采样间隔Δt的确定依据为：Δt< 0.1tsmin，tl>10 tsmax，采样点数nt的计算方法为：nt＝tl/Δt；
其中，tsmin＝min{nts1，nts2，......，ntsi，......}，tsmax＝nts1+nts2+......+ ntsi+......，ntsi＝hi/vsi/Δt。
步骤3中地震波的射线参数pw的计算方法为：pw＝sinθ/vsn，且pwvpmax <1，其中，vpmax＝max{vp1，vp2，......，vpi，......}。
步骤4中地震纵波和横波在工程场地各土层中的传播时间样点数通过公 式1计算得到：
Figure BDA0001702104110000091
其中，hi、vpi、vsi分别是地层i的厚度、纵波速度、横波速度；pw是射 线参数，由步骤3中计算方法确定。
步骤4中地震纵波和横波斜入射至地面处产生的反射系数
Figure BDA0001702104110000098
Figure BDA0001702104110000096
和
Figure BDA0001702104110000097
分别通过公式2和3计算得到：
Figure BDA0001702104110000092
Figure BDA0001702104110000093
其中，
Figure BDA0001702104110000095
分别是纵波入射至地面产生的同类波，记为P波、转 换波，记为S波的反射系数；
Figure BDA0001702104110000094
分别是横波入射至自由面处产生的 P波、S波的反射系数，上标“pp”表示入射P波产生的同类波，“ps”表示 入射P波产生的转换波，“sp”表示入射S波产生的转换波，“ss”表示入射S 波产生的同类波；“+”表示地震波从下向上入射，“-”表示地震波从上向下 入射；vp1、vs1分别是地层1的纵波速度和横波速度；α1、β1分别是纵波、 横波传播至地面处的入射角和反射角。
步骤5中地震纵波从第i层入射至第i-1层时在界面i处的反射系数
Figure BDA0001702104110000106
Figure BDA0001702104110000104
和透射系数
Figure BDA0001702104110000105
通过公式4计算得到：
Figure BDA0001702104110000101
步骤5中地震横波从第i层入射至第i-1层时在界面i处的反射系数
Figure BDA0001702104110000107
Figure BDA0001702104110000108
和透射系数
Figure BDA0001702104110000109
通过公式5计算得到：
Figure BDA0001702104110000102
地震纵波从第i-1层入射至第i层时在界面i处的反射系数
Figure BDA00017021041100001010
和 透射系数
Figure BDA00017021041100001011
通过公式6计算得到：
Figure BDA0001702104110000103
地震横波从第i-1层入射至第i层时在界面i处的反射系数
Figure BDA00017021041100001012
和透射系数
Figure BDA00017021041100001013
通过公式7计算得到：
Figure BDA0001702104110000111
公式4-7中，i是地层i-1与地层i之间的地层界面；vpi、vsi分别是地层 i的纵波速度、横波速度，vpi-1、vsi-1分别是地层i-1的纵波速度、横波速度； αi-1、αi分别是纵波在界面i处的入射角或反射角、透射角；βi-1、βi是横波在 界面i处的入射角或反射角)、透射角；i＝2，3，......，n。
步骤6中纵波和横波在各个土层界面的脉冲响应
Figure BDA0001702104110000112
和
Figure BDA0001702104110000113
通过公式8 计算得到：
Figure BDA0001702104110000114
其中，
Figure BDA0001702104110000115
和
Figure BDA0001702104110000116
分别表示从第i层入射至第i-1层时界面i处的纵波和横 波的响应；
Figure BDA0001702104110000117
和
Figure BDA0001702104110000118
别是从第i-1层入射至第i层时界面i处纵波和横波 的响应；i＝2，......，n-1；j＝1，2，......，nt；
当i＝1时，纵波和横波入射至地表面处的脉冲响应
Figure BDA0001702104110000119
和
Figure BDA00017021041100001110
如公式9所 示：
Figure BDA00017021041100001111
其中，
Figure BDA0001702104110000121
和
Figure BDA0001702104110000122
分别是从第1层入射至地表面处的纵波和横波；
Figure BDA0001702104110000123
和
Figure BDA0001702104110000124
分别是从空气介质入射至地表面处的纵波和横波，j＝1，2，......，nt；
当i＝n时，公式8表示为公式10：
Figure BDA0001702104110000125
其中，
Figure BDA0001702104110000126
和
Figure BDA0001702104110000127
分别是从第n层入射至第n-1层时界面n处的纵波和 横波响应；
Figure BDA0001702104110000128
和
Figure BDA0001702104110000129
分别是从第n-1层入射至第n层时界面n处的纵波和 横波响应；j＝1，2，......，nt。
步骤6中合成地面水平方向的脉冲响应时程
Figure BDA00017021041100001210
通过公式11得到：
Figure BDA00017021041100001211
其中，j＝1，2，...，nt。
本发明的倾斜地震波作用下场地卓越周期的计算方法主要是将基岩入 射至场地的地震波作为输入信号，地面震动作为输出信号，把地震场地看作 地震信号滤波器，当输入信号为脉冲波时的地面震动称之为脉冲响应，其傅 里叶频谱函数就是滤波器的频率特性函数，据此得到的函数周期就是场地的 卓越周期。由该方法得到的卓越周期不仅包含了周期信息，而且还包含了入 射方式和入射波型方面的信息，可为实际工程的抗震设计提供更全面的场地 自振特性参数，而由传统方法测定的卓越周期只是SV波垂直入射时的一个计算特例。
具体实施例如下：
实施例1
某二层场地的地层参数如表1所示，下面以SV波20°入射角为例，详 细介绍在倾斜地震波作用下场地卓越周期的确定过程，具体步骤如下：
(1)确定地层参数，如表1所示；
(2)确定脉冲响应持时采样点数，本例中，取脉冲响应持时tl＝5s，时 间采样间隔△t＝0.001s，采样点数nt＝5000；
(3)计算射线参数pw，选择SV入射角θ＝20°，地层层数n＝20，vs2＝ 500 m/s，则pw＝sin20/500＝0.00068404 s/m；
(4)计算地震波在第1层单程传播的时间：
Figure BDA0001702104110000131
Figure BDA0001702104110000132
(5)计算地震波入射至地面的反射系数；根据斯奈尔定理，计算入射 角α1和β1，α1＝arc sin(vp1sin20°/vs2)＝arc sin(346sin20°/500)＝13.6905°；β1＝ arc sin(vs1sin20°/vs2)＝arc sin(200sin20°/500)＝7.8632°，将入射角α1＝ 13.6905°，β1＝7.8632°，vp1＝346 m/s，vs1＝200 m/s，代入上述公式2和3， 计算可得：R1 pp+＝-0.9140、R1 ps +＝0.5286、R1 sp+＝-0.3115、R1 ss+＝-0.9140； 根据上述公式4～7，计算可得：R2 pp+＝-0.2932、R2 ps+＝0.5855、T2 pp+＝1.3701、 T2 ps+＝-0.6507、R2 sp+＝-0.3943、R2 ss+＝-0.2460、T2 sp +＝0.3800、T2 ss+＝1.5086、 R2 pp-＝0.5326、R2 ps-＝0.2749、T2 pp-＝0.4401、T2 ps-＝-0.1812、R2 sp-＝-0.1620、 R2 ss-＝0.4853、T2 sp-＝0.1232、T2 ss-＝0.4241；
(6)根据上述公式8～10计算得到
Figure BDA0001702104110000133
j＝1，2，......，5000，如 图6所示；根据
Figure BDA0001702104110000141
合成地面水平方向的脉冲响应时程
Figure BDA0001702104110000142
j＝1，2，......， 5000，如图7所示；
(7)通过离散傅里叶变换(DFT)获得地震脉冲响应的Fourier谱，如 图8所示；
(8)根据地震脉冲响应的Fourier谱的峰值极大值所对应的时间，确定 场地卓越周期，前三阶卓越周期分别是0.968 s、0.321 s和0.206 s，基本卓 越周期为0.968 s；选取0°、5°、10°、15°、25°、30°、33°入射角，分别重复 上述步骤1～10，可以得到不同入射角条件下的场地卓越周期，计算结果如 表2所示。
表1场地1地层参数
Figure BDA0001702104110000143
表2场地1卓越周期
Figure BDA0001702104110000144
实施例2
某7层场地的地层参数如表3所示，选取0°、5°、10°、15°、25°、30°、 32°入射角，分别按照步骤1～8计算，可以得到不同入射角条件下的场地卓 越周期，计算结果如表4所示。
表3场地2地层参数
Figure BDA0001702104110000151
表4场地2卓越周期
Figure BDA0001702104110000152
实施例3
某11层场地的地层参数如表5所示，选取0°、5°、10°、15°、25°、30°、 35°入射角，分别按照步骤1～8计算，可以得到不同入射角条件下的场地卓 越周期，计算结果如表6所示。
实施例4
某9层场地的地层参数如表7所示，选取0°、5°、10°、15°、24°入射角， 分别按照步骤1～8计算，可以得到不同入射角条件下的场地卓越周期，计 算结果如表8所示。
表5场地3地层参数
Figure BDA0001702104110000161
表6场地3卓越周期
Figure BDA0001702104110000162
实施例5
某10层场地的地层参数如表9所示，选取0°、5°、9°入射角，分别按照 步骤1～10计算，可以得到不同入射角条件下的场地卓越周期，计算结果如 表10所示。
表7场地4地层参数
Figure BDA0001702104110000171
表8场地4卓越周期
Figure BDA0001702104110000172
表9场地5地层参数
Figure BDA0001702104110000181
表10场地5卓越周期
Figure BDA0001702104110000182

Claims (4)

1.倾斜地震波作用下场地卓越周期的计算方法，其特征在于，具体操作步骤如下：
步骤1.根据岩土工程的勘查结果，确定场地地层数n、地层的密度ρi、纵波速度vpi、横波速度vsi以及地层厚度hi，其中i＝1，2，......，n；
步骤2.根据需要计算的精度要求，确定地震脉冲响应持时tl和时间采样间隔Δt，并计算采样点数nt；
步骤3.根据实际工程勘察资料选择地震横波入射角θ，计算地震波的射线参数pw；
步骤4.分别计算地震纵波和横波在场地各地层中的传播时间以及其入射至地面处产生的反射系数R1 pp+、R1 ps+、R1 sp+和R1 ss+；
步骤5.分别计算地震纵波和横波斜入射至各地层界面i处的反射系数
Figure FDA0002275349920000011
和透射系数
Figure FDA0002275349920000012
Figure FDA0002275349920000013
当i＝2，......，n时，具体如下：
地震纵波从第i层入射至第i-1层时在界面i处的反射系数
Figure FDA0002275349920000014
和透射系数
Figure FDA0002275349920000015
通过公式(4)计算得到：
Figure FDA0002275349920000016
地震横波从第i层入射至第i-1层时在界面i处的反射系数
Figure FDA0002275349920000017
和透射系数
Figure FDA0002275349920000018
通过公式(5)计算得到：
Figure FDA0002275349920000021
地震纵波从第i-1层入射至第i层时在界面i处的反射系数
Figure FDA0002275349920000022
和透射系数
Figure FDA0002275349920000023
通过公式(6)计算得到：
Figure FDA0002275349920000024
地震横波从第i-1层入射至第i层时在界面i处的反射系数
Figure FDA0002275349920000025
和透射系数
Figure FDA0002275349920000026
通过公式(7)计算得到：
Figure FDA0002275349920000027
公式(4)～公式(7)中，i是地层与地层界面的序号，i取i＝2，......，n；vpi、vsi分别是地层i的纵波速度、横波速度，vpi-1、vsi-1分别是地层i-1的纵波速度、横波速度；
公式(4)中，αi、βi、αi-1、βi-1分别是当纵波从第i层入射至第i-1层时，在界面i处的纵波反射角、横波反射角、纵波透射角、横波透射角；
公式(5)中，αi、βi、αi-1、βi-1分别是当横波从第i层入射至第i-1层时，在界面i处的纵波反射角、横波反射角、纵波透射角、横波透射角；
公式(6)中，αi、βi、αi-1、βi-1分别是当纵波从第i-1层入射至第i层时，在界面i处的纵波透射角、横波透射角、纵波反射角、横波反射角；
公式(7)中，αi、βi、αi-1、βi-1分别是当横波从第i-1层入射至第i层时，在界面i处的纵波透射角、横波透射角、纵波反射角、横波反射角；
步骤6.根据步骤4和步骤5得到的反射系数和透射系数，计算得到各个地层界面的脉冲响应
Figure FDA0002275349920000031
以及地面脉冲响应
Figure FDA0002275349920000032
和
Figure FDA0002275349920000033
其中，i取i＝2，......，n；时间样点序号j＝1，2，......，nt；根据R1 pp+、R1 ps+、R1 sp+、R1 ss+、
Figure FDA0002275349920000034
计算得到地面水平方向的脉冲响应时程
Figure FDA0002275349920000035
具体过程如下：
纵波和横波在各个地层界面的脉冲响应
Figure FDA0002275349920000036
和
Figure FDA0002275349920000037
通过公式(8)计算得到：
Figure FDA0002275349920000038
其中，
Figure FDA0002275349920000039
和
Figure FDA00022753499200000310
分别表示从第i层入射至第i-1层时界面i处的纵波的脉冲响应和横波的脉冲响应；
Figure FDA00022753499200000311
和
Figure FDA00022753499200000312
分别是从第i-1层入射至第i层时界面i处纵波的脉冲响应和横波的脉冲响应；i取i＝2，......，n-1；j＝1，2，......，nt；
当i＝1时，纵波和横波入射至地面处的脉冲响应
Figure FDA00022753499200000313
和
Figure FDA00022753499200000314
如公式(9)所示：
Figure FDA0002275349920000041
其中，
Figure FDA0002275349920000042
和
Figure FDA0002275349920000043
分别是从第1层入射至地面处的纵波的脉冲响应和横波的脉冲响应；
Figure FDA0002275349920000044
和
Figure FDA0002275349920000045
分别是从空气介质入射至地面处的纵波的脉冲响应和横波的脉冲响应，j＝1，2，......，nt；
当i＝n时，公式(8)表示为公式(10)：
Figure FDA0002275349920000046
其中，
Figure FDA0002275349920000047
和
Figure FDA0002275349920000048
分别是从第n层入射至第n-1层时界面n处的纵波的脉冲响应和横波的脉冲响应；
Figure FDA0002275349920000049
和
Figure FDA00022753499200000410
分别是从第n-1层入射至第n层时界面n处的纵波的脉冲响应和横波的脉冲响应；j＝1，2，......，nt；
步骤6中合成地面水平方向的脉冲响应时程
Figure FDA00022753499200000411
通过公式(11)得到：
Figure FDA00022753499200000412
其中，j＝1，2，...，nt；
步骤7.将步骤6得到脉冲响应时程
Figure FDA00022753499200000413
通过进行离散Fourier变换(DFT)获得脉冲地震动的Fourier谱F(u)；
步骤8.通过分析步骤7得到的脉冲地震动的Fourier谱F(u)，确定场地卓越周期，即在Fourier谱中出现的若干个不连续尖峰状谱中，选取尖峰极大值点所对应的周期作为场地卓越周期，其中，最大卓越周期为基本卓越周期。
2.如权利要求1所述的倾斜地震波作用下场地卓越周期的计算方法，其特征在于，所述步骤2中地震脉冲响应持时tl和时间采样间隔Δt的确定依据为：Δt<0.1tsmin，tl>10tsmax，其中，tsmin＝min{nts1，nts2，......，ntsi，......}，tsmax＝nts1+nts2+......+ntsi+......，ntsi＝hi/vsi/Δt；
采样点数nt的计算方法为：nt＝tl/Δt。
3.如权利要求1所述的倾斜地震波作用下场地卓越周期的计算方法，其特征在于，所述步骤3中地震波的射线参数pw的计算方法为：pw＝sinθ/vsn，且pwvpmax<1，其中，vpmax＝max{vp1，vp2，......，vpi，......}。
4.如权利要求1所述的倾斜地震波作用下场地卓越周期的计算方法，其特征在于，所述步骤4中地震纵波和横波在工程场地各地层中的传播时间通过公式1计算得到：
Figure FDA0002275349920000051
其中，hi是地层i的厚度，vpi、vsi分别是地震波在地层i中的纵波速度、横波速度；
所述步骤4中地震纵波和横波斜入射至地面处产生的反射系数R1 pp+、R1 ps+、R1 sp+和R1 ss+分别通过公式2和3计算得到：
Figure FDA0002275349920000052
Figure FDA0002275349920000061
其中，R1 pp+、R1 ps+分别是纵波入射至地面引起的反射P波的反射系数、反射S波的反射系数；R1 sp+、R1 ss+分别是横波入射至地面处产生的反射P波的反射系数、反射S波的反射系数；上标“pp”表示入射P波产生的同类波，“ps”表示入射P波产生的转换波，“sp”表示入射S波产生的转换波，“ss”表示入射S波产生的同类波；“+”表示地震波从下向上入射，“-”表示地震波从上向下入射；vp1、vs1分别是第一层中的纵波速度和横波速度；α1、β1分别是地面处纵波反射P波的反射角、横波反射S波的反射角。
CN201810638788.7A 2018-06-20 2018-06-20 倾斜地震波作用下场地卓越周期的计算方法 Active CN108957537B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810638788.7A CN108957537B (zh) 2018-06-20 2018-06-20 倾斜地震波作用下场地卓越周期的计算方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810638788.7A CN108957537B (zh) 2018-06-20 2018-06-20 倾斜地震波作用下场地卓越周期的计算方法

Publications (2)

Publication Number Publication Date
CN108957537A CN108957537A (zh) 2018-12-07
CN108957537B true CN108957537B (zh) 2020-03-31

Family

ID=64490663

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810638788.7A Active CN108957537B (zh) 2018-06-20 2018-06-20 倾斜地震波作用下场地卓越周期的计算方法

Country Status (1)

Country Link
CN (1) CN108957537B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108828661B (zh) * 2018-03-06 2019-10-25 西安理工大学 基于地震脉冲响应谱测定场地卓越周期的方法
CN110609323B (zh) * 2019-09-24 2022-02-22 西安理工大学 一种基于地层信息的场地卓越周期的精简计算方法
CN112444875B (zh) * 2020-10-28 2023-10-03 西安理工大学 一种获取场地卓越周期的精确解的方法
CN113820748B (zh) * 2021-08-23 2024-03-26 西安理工大学 一种地震sv波发生全反射时的地面滑行波时程获取方法

Family Cites Families (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN2079760U (zh) * 1990-08-14 1991-06-26 福建省建筑设计院 一种地脉动波频谱测量仪
JP3859587B2 (ja) * 2001-12-12 2006-12-20 東京瓦斯株式会社 配管耐震性評価方法および配管耐震性評価装置
US7131238B2 (en) * 2003-07-21 2006-11-07 Fm Global Technologies, Llc Method of testing seismic braces
CN101354444B (zh) * 2007-07-25 2011-02-09 中国石油天然气集团公司 一种确定地层岩性和孔隙流体的方法
CN103792059A (zh) * 2012-10-29 2014-05-14 同济大学 利用多点振动台模拟地下结构非一直激励的分段式模型箱
US10677944B2 (en) * 2014-08-23 2020-06-09 Trimble Inc. Earthquake and displacement early warning system
CN104807607B (zh) * 2015-04-23 2017-10-24 成都畅达通检测技术股份有限公司 基于激励状态时响应谱特性的构筑物病害检测系统及其检测方法
CN105421609B (zh) * 2016-01-14 2017-07-28 同济大学 半主动变刚度变阻尼混合阻尼器
CN107119959A (zh) * 2017-06-30 2017-09-01 安徽理工大学 组合式单摆式调谐质量阻尼器结构系统及施工方法
CN108089225B (zh) * 2017-12-02 2019-08-13 桂林理工大学 一种基于单台站初至p波的震级快速估算方法

Also Published As

Publication number Publication date
CN108957537A (zh) 2018-12-07

Similar Documents

Publication Publication Date Title
CN108957537B (zh) 倾斜地震波作用下场地卓越周期的计算方法
CN105308479B (zh) 通过与偏移距相关的弹性fwi的多参数反演
Sabra et al. Emergence rate of the time-domain Green’s function from the ambient noise cross-correlation function
WO2017024702A1 (zh) 一种射线弹性参数的反演系统
US9075163B2 (en) Interferometric seismic data processing
US11650342B2 (en) Method for exploring passive source seismic frequency resonance
US20050254343A1 (en) Methods for processing dispersive acoustic waveforms
US20050152220A1 (en) Seismic data processing
WO2010118624A1 (zh) 地震波速度场构建中井约束横向可变h-v曲线构建方法
CN104820218B (zh) 一种基于频域自相关的浅海海底单参数反演方法
CN111025387B (zh) 一种页岩储层的叠前地震多参数反演方法
Duncan et al. Characteristics of sound propagation in shallow water over an elastic seabed with a thin cap-rock layer
EP2073041A1 (en) Method to estimate a seismic ray parameter for a seismogram
CN112327358B (zh) 一种粘滞性介质中声波地震数据正演模拟方法
CN109946742B (zh) 一种TTI介质中纯qP波地震数据模拟方法
Bodet et al. Elasticity profile of an unconsolidated granular medium inferred from guided waves: Toward acoustic monitoring of analogue models
CN106896409A (zh) 一种基于波动方程边值反演的变深度缆鬼波压制方法
Duroux et al. Extracting guided waves from cross-correlations of elastic diffuse fields: Applications to remote structural health monitoring
Marsset et al. Deep-towed high resolution seismic imaging II: Determination of P-wave velocity distribution
Du et al. Observation and inversion of very-low-frequency seismo-acoustic fields in the South China Sea
CN112540403A (zh) 基于模拟地震软弱层的应力计算及动力响应阶段划分方法
CN110376646B (zh) 一种基于曲坐标系纵横波解构方程的弹性棱柱波逆时偏移成像方法
CN112558159A (zh) 一种声波测井初至拾取方法
JP2020177011A5 (ja) 構造解析装置
Li et al. Virtual head waves in ocean ambient noise: Theory and modeling

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant