CN108828661A - 基于地震脉冲响应谱测定场地卓越周期的方法 - Google Patents

基于地震脉冲响应谱测定场地卓越周期的方法 Download PDF

Info

Publication number
CN108828661A
CN108828661A CN201810184529.1A CN201810184529A CN108828661A CN 108828661 A CN108828661 A CN 108828661A CN 201810184529 A CN201810184529 A CN 201810184529A CN 108828661 A CN108828661 A CN 108828661A
Authority
CN
China
Prior art keywords
seismic
time
pulse
predominant period
pulse response
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.)
Granted
Application number
CN201810184529.1A
Other languages
English (en)
Other versions
CN108828661B (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 CN201810184529.1A priority Critical patent/CN108828661B/zh
Publication of CN108828661A publication Critical patent/CN108828661A/zh
Application granted granted Critical
Publication of CN108828661B publication Critical patent/CN108828661B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/307—Analysis for determining seismic attributes, e.g. amplitude, instantaneous phase or frequency, reflection strength or polarity
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/63—Seismic attributes, e.g. amplitude, polarity, instant phase

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

本发明公开的基于地震脉冲响应谱测定场地卓越周期的方法，具体操作步骤如下：首先确定场地土层数、地层密度、剪切速度以及地层厚度、地震脉冲响应持时和时间采样间隔；然后计算地震脉冲响应持时采样点数、地震波在各土层中传播时间采样点数、地震波在地层界面处的反射系数并确定脉冲地震动响应时程；最后将脉冲地震动响应时程通过离散傅里叶变换获得脉冲地震动的傅里叶谱，通过脉冲地震动的Fourier谱确定场地卓越周期。本发明基于地震脉冲响应谱测定场地卓越周期的方法，解决了现有技术需要使用实际地震记录的缺点，以及处理复杂工程场地时计算结果不准确、求解过程困难等问题，本发明的方法不受土层层数的限制，计算精度高且速度快。

Description

基于地震脉冲响应谱测定场地卓越周期的方法
技术领域
本发明属于岩土工程勘察技术领域，涉及一种基于地震脉冲响应谱测定场地卓越周期的方法。
背景技术
卓越周期是指与地基土层发生共振作用的地震波的一个谐波分量。目前，确定场地卓越周期的方法大致分为直接测定法和波速法两大类。直接测定法是通过对强震记录或地脉动记录进行Fourier分析确定卓越周期，但这类方法需要有良好的强震或地脉动记录，且容易受场地局部条件的限制或人为因素的干扰；波速法则是应用理论方法、简化方法或数值方法计算得到卓越周期，应用这类方法处理复杂工程场地时，可能面临计算结果不准确或求解过程困难等问题。
发明内容
本发明的目的是提供一种基于地震脉冲响应谱测定场地卓越周期的方法，解决了现有技术需要使用实际地震记录的缺点，以及处理复杂工程场地时计算结果不准确、求解过程困难的问题。
本发明所采用的技术方案是，基于地震脉冲响应谱测定场地卓越周期的方法，具体操作步骤如下：
步骤1.根据岩土工程的勘查结果，确定场地土层数n、地层的密度ρi、剪切速度vi以及地层厚度hi，其中，i＝1，2，......，n；
步骤2.根据需要计算的精度要求，确定地震脉冲响应持时tl和时间采样间隔Δt；
步骤3.根据步骤2得到的地震脉冲响应持时tl和时间采样间隔Δt，计算地震脉冲响应持时采样点数nt；
步骤4.结合步骤1和步骤2确定的参数，计算地震波在各土层中传播时间采样点数nti和地震波在地层界面处的反射系数Ri，其中，i＝1，2，......，n；
步骤5.根据步骤4得到的各土层中传播时间采样点数nti和地震波在地层界面处的反射系数Ri，计算各个土层界面的次级地震波时程ui d和
步骤6.通过步骤5得到的界面的次级地震波时程u0 d，通过函数关系转换得到脉冲地震动响应时程u0；
步骤7.将步骤6得到的脉冲地震动响应时程u0通过进行离散Fourier变换(DFT)获得脉冲地震动的Fourier谱F(u)；
步骤8.通过分析步骤7得到的脉冲地震动的Fourier谱F(u)，确定场地卓越周期，即在Fourier谱中出现的若干个不连续尖峰状谱中，选取尖峰极大值点所对应的周期作为场地卓越周期，其中，最大卓越周期为基本卓越周期。
本发明的其他特点还在于，
步骤2中地震脉冲响应持时tl和时间采样间隔Δt的确定依据为：
Δt<0.1tmin，tl>10tmax.，其中，tmin＝min{nt1，nt2，......，nti，......}，tmax＝nt1+nt2+......+nti+......，其中i＝1，2，......，n。
步骤3中地震脉冲响应持时采样点数nt的计算方法如公式1：
nt＝tl/Δt (1)。
步骤4中地震波在各土层中传播时间采样点数nti通过公式2计算得到：
nti＝hi/vi/Δt (2)
所述地震波在地层界面处的反射系数Ri通过公式3计算得到：
Ri＝(ρi+1vi+1-ρivi)/(ρi+1vi+1+ρivi) (3)
其中，ρi为第i地层的密度，vi为第i地层的剪切速度，ρi+1为第i+1地层的密度，vi+1为第i+1地层的剪切速度，hi为第i地层的地层厚度，Δt为时间采样间隔。
步骤5中各个土层界面的次级地震波时程ui d和通过递推公式4-9计算得到：
u0 d(j)＝-u1 p(j-nt1) (4)
u1 d(j)＝(1+R1)u0 d(j-nt1)-R1u2 p(j-nt2) (5)
u1 p(j)＝R1u0 d(j-nt1)+(1-R1)u2 p(j-nt2) (6)
………………………………
ui d(j)＝(1+Ri)ui-1 d(j-nti)-Riui+1 p(j-nti+1) (7)
ui p(j)＝Riui-1 d(j-nti)+(1-Ri)ui+1 p(j-nti+1) (8)
………………………………
un p(j)＝Rnun d(j-ntn) (9)
其中，j＝1，2，......，nt，un p(1)＝1，u0 d，u1 d，u1 p，……，ui d，……，un p是关于j的函数，分别表示地震作用下各土层之间的界面0，界面1，……，界面i，……，界面n产生的次级地震波时程；上标d，p分别表示向下，向上传播的次级波；i表示界面编号，j表示时间样点序号。
步骤6中脉冲地震动时程u0通过公式10的函数关系得到：
u0＝2u0 d(j) (10)
其中，j＝1，2，...，nt。
本发明的有益效果是，基于地震脉冲响应谱测定场地卓越周期的方法，解决了现有技术需要使用实际地震记录的缺点，以及处理复杂工程场地时计算结果不准确、求解过程困难的问题。根据大地滤波作用，将场地土视为地震波信号滤波器，并以从基岩入射至场地的地震波为输入，入射地震波引起的地面震动为输出。当输入脉冲波时，此时输出的地震动信号称为地震脉冲响应，其傅立叶(Fourier)谱即是场地滤波器的频率特性函数，通过对地震脉冲响应谱分析可确定场地的卓越周期。
附图说明
图1是本发明的基于地震脉冲响应谱测定场地卓越周期的方法中场地卓越周期的计算流程图；
图2是本发明的基于地震脉冲响应谱测定场地卓越周期的方法中场地反应模型；
图3是本发明的基于地震脉冲响应谱测定场地卓越周期的方法的实施例1的场地脉冲响应时程；
图4是本发明的基于地震脉冲响应谱测定场地卓越周期的方法中实施例1的场地脉冲响应谱；
图5是本发明的基于地震脉冲响应谱测定场地卓越周期的方法中实施例2的场地脉冲响应时程；
图6是本发明的基于地震脉冲响应谱测定场地卓越周期的方法中实施例2的场地脉冲响应谱；
图7是本发明的基于地震脉冲响应谱测定场地卓越周期的方法中实施例3的场地脉冲响应时程；
图8是本发明的基于地震脉冲响应谱测定场地卓越周期的方法中实施例3的场地脉冲响应谱；
图9是本发明的基于地震脉冲响应谱测定场地卓越周期的方法中实施例4的场地脉冲响应时程；
图10是本发明的基于地震脉冲响应谱测定场地卓越周期的方法中实施例4的场地脉冲响应谱；
图11是本发明的基于地震脉冲响应谱测定场地卓越周期的方法中实施例5的场地脉冲响应时程；
图12是本发明的基于地震脉冲响应谱测定场地卓越周期的方法中实施例5的场地脉冲响应谱。
具体实施方式
下面结合附图和具体实施方式对本发明进行详细说明。
本发明的基于地震脉冲响应谱测定场地卓越周期的方法，具体操作步骤如下：
步骤1.根据岩土工程的勘查结果，确定场地土层数n、地层的密度ρi、剪切速度vi以及地层厚度hi，其中，i＝1，2，......，n；
步骤2.根据需要计算的精度要求，确定地震脉冲响应持时tl和时间采样间隔Δt。理论上，Δt取值越小，tl取值越大，计算精度越高；
步骤3.根据步骤2得到的地震脉冲响应持时tl和时间采样间隔Δt，计算地震脉冲响应持时采样点数nt；
步骤4.结合步骤1和步骤2确定的参数，计算地震波在各土层中传播时间采样点数nti和地震波在地层界面处的反射系数Ri，其中，i＝1，2，......，n；
步骤5.根据步骤4得到的各土层中传播时间采样点数nti和地震波在地层界面处的反射系数Ri，计算各个土层的界面次级地震波时程ui d和其中，ui d，分别是地震作用下由土层界面i产生的向下、向上传播的次级地震波。
步骤6.通过步骤5得到的界面次级地震波时程u0 d，通过函数关系转换得到脉冲地震动响应时程u0；
步骤7.将步骤6得到的脉冲地震动响应时程u0通过进行离散Fourier变换(DFT)获得脉冲地震动的Fourier谱F(u)；
步骤8.通过分析步骤7得到的脉冲地震动的Fourier谱F(u)，确定场地卓越周期，即在Fourier谱中出现的若干个不连续尖峰状谱中，选取尖峰极大值点所对应的周期作为场地卓越周期，其中，最大卓越周期为基本卓越周期。
步骤2中时间采样间隔Δt的确定依据为：
在保证计算精度的前提下，Δt<0.1tmin，tl>10tmax，其中，tmin＝min{nt1，nt2，......，nti，......}，tmax＝nt1+nt2+......+nti+......，其中i＝1，2，......，n。
步骤3中地震脉冲响应持时采样点数nt的计算方法如公式1：
nt＝tl/Δt (1)。
步骤4中地震波在各土层中传播时间的采样点数nti通过公式2计算得到：
nti＝hi/vi/Δt (2)
所述地震波在地层界面处的反射系数Ri通过公式3计算得到：
Ri＝(ρi+1vi+1-ρivi)/(ρi+1vi+1+ρivi) (3)
其中，ρi为第i地层的密度，vi为第i地层的剪切速度，ρi+1为第i+1地层的密度，vi+1为第i+1地层的剪切速度，hi为第i地层的地层厚度，Δt为时间采样间隔。
步骤5中各个土层界面次级地震波时程ui d和通过递推公式4-9计算得到：
u0 d(j)＝-u1 p(j-nt1) (4)
u1 d(j)＝(1+R1)u0 d(j-nt1)-R1u2 p(j-nt2) (5)
u1 p(j)＝R1u0 d(j-nt1)+(1-R1)u2 p(j-nt2) (6)
………………………………
ui d(j)＝(1+Ri)ui-1 d(j-nti)-Riui+1 p(j-nti+1) (7)
ui p(j)＝Riui-1 d(j-nti)+(1-Ri)ui+1 p(j-nti+1) (8)
………………………………
un p(j)＝Rnun d(j-ntn) (9)
其中，j＝1，2，......，nt，un p(1)＝1，u0 d，u1 d，u1 p，……，ui d，……，un p是关于j的函数，分别表示地震作用下各土层之间的界面0，界面1，……，界面i，……，界面n产生的次级地震波时程；上标d，p分别表示向下，向上传播的次级波；i表示界面编号，j表示时间样点序号。
步骤6中脉冲地震动时程u0通过公式10的函数关系得到：
u0＝2u0 d(j) (10)
其中，j＝1，2，...，nt。
本发明的基于地震脉冲响应谱测定场地卓越周期的方法是综合利用直接测定法和波速法两者的优点，分两步完成对场地卓越周期的测定，第一步是利用场地土的厚度、密度、波速参数，计算得到脉冲震源作用下的人造地震动时程，克服了缺失场地实际地震记录的不足；第二步是对人造地震动时程进行Fourier分析得到卓越周期，在保证计算精度的前提下避免了复杂的计算过程。此方法不受土层层数的限制，具有计算精度高、计算速度快、方便直观判断等特点。
实施例1
某三层土层场地，土厚20m，有两类土，各个土层参数见表1所示，
表1场地地层
场地卓越周期计算如下：
(1)确定地层参数；
n＝3；
ρ1＝2041、v1＝200、h1＝4；
ρ2＝2143、v2＝300、h2＝4；
ρ3＝2041、v3＝200、h3＝12；
(2)确定时间参数；
tl＝5、Δt＝0.0001
(3)计算地震脉冲响应持时采样点数；
nt＝50000
(4)计算地震波在各层中传播时间采样点数和反射系数
nt1＝200、nt2＝133、nt3＝600；R1＝0.2233、R2＝-0.2233、R3＝0.85
(5)根据递推公式4-9，计算得到各个土层界面的所激发的次级反射波时程u0 d(j)，u1 d(j)，u1 p(j)，u2 d(j)，u2 p(j)，u3 p(j)，其中j＝1，2，......，50000；
(6)通过步骤5得到的u0 d(j)，根据公式u0＝2u0 d(j)计算得到地震动时程u0，结果如图3所示；
(7)将得到的脉冲地震动记录数据通过DFT获得地震脉冲响应的Fourier谱，计算结果如图4所示；
(8)根据地震脉冲响应的Fourier谱的峰值极大值所对应的时间，确定场地卓越周期，由图4可得，场地前4阶的卓越周期分别为0.393s，0.119s，0.074s和0.056s。
为了验证此项发明的效果，与其他方法得到的结果进行了对比，详见表2所示，由表2可见，各种方法计算所得到的误差均不超过1％。
表2场地卓越周期Ti(秒)
表3场地地层
实施例2
某2层土场地，土厚20m，各个土层参数见表3所示。按照上述实施例1步骤计算所得的地脉冲响应时程及其响应谱分别如图5与图6所示，场地前4阶的卓越周期分别为0.402s，0.131s，0.076s和0.053s，与实际结果一致。
实施例3
一个5层土场地，各个土层参数见表4所示，按照上述实施例1步骤计算所得的地脉冲响应时程及其响应谱分别如图7与图8所示，基本卓越周期为0.21s，与地球脉动法所得场地卓越周期0.25s接近。
实施例4
一个地层较多的工程场地，共12层，土厚126m，各个土层参数见表5所示，基本卓越周期为1.316s。按照上述实施例1步骤计算所得的地脉冲响应时程及其响应谱分别如图9与图10所示。由图10得场地前4阶的卓越周期分别为1.319s，0.590s，0.361s和0.259s，基本卓越周期(1.319s)与实际场地结果(1.316s)吻合良好。
实施例5
一个非常复杂的多层场地，共20层，土厚250m，各个土层参数见表6所示，按照上述实施例1步骤计算所得的地脉冲响应时程及其响应谱分别如图11与图12所示。计算得到场地前4阶的卓越周期分别为2.032s，0.845s，0.536s和0.392s。
表4场地地层
表5场地地层
表6场地地层

Claims (6)

1.基于地震脉冲响应谱测定场地卓越周期的方法，其特征在于，具体操作步骤如下：
步骤1.根据岩土工程的勘查结果，确定场地土层数n、地层的密度ρi、剪切速度vi以及地层厚度hi，其中，i＝1，2，......，n；
步骤2.根据需要计算的精度要求，确定地震脉冲响应持时tl和时间采样间隔Δt；
步骤3.根据步骤2得到的地震脉冲响应持时tl和时间采样间隔Δt，计算地震脉冲响应持时采样点数nt；
步骤4.结合步骤1和步骤2确定的参数，计算地震波在各土层中传播时间采样点数nti和地震波在地层界面处的反射系数Ri，其中，i＝1，2，......，n；
步骤5.根据步骤4得到的各土层中传播时间采样点数nti和地震波在地层界面处的反射系数Ri，计算各个土层界面的次级地震波时程ui d和
步骤6.通过步骤5得到的界面的次级地震波时程u0 d，通过函数关系转换得到脉冲地震动响应时程u0；
步骤7.将步骤6得到的脉冲地震动响应时程u0通过进行离散Fourier变换(DFT)获得脉冲地震动的Fourier谱F(u)；
步骤8.通过分析步骤7得到的脉冲地震动的Fourier谱F(u)，确定场地卓越周期，即在Fourier谱中出现的若干个不连续尖峰状谱中，选取尖峰极大值点所对应的周期作为场地卓越周期，其中，最大卓越周期为基本卓越周期。
2.如权利要求1所述的基于地震脉冲响应谱测定场地卓越周期的方法，其特征在于，所述步骤2中地震脉冲响应持时tl和时间采样间隔Δt的确定依据为：
Δt<0.1tmin，tl>10tmax.，其中，tmin＝min{nt1，nt2，......，nti，......}，tmax＝nt1+nt2+......+nti+......，其中i＝1，2，......，n。
3.如权利要求1所述的基于地震脉冲响应谱测定场地卓越周期的方法，其特征在于，所述步骤3中地震脉冲响应持时采样点数nt的计算方法如公式1：
nt＝tl/Δt (1)。
4.如权利要求1所述的基于地震脉冲响应谱测定场地卓越周期的方法，其特征在于，所述步骤4中地震波在各土层中传播时间采样点数nti通过公式2计算得到：
nti＝hi/vi/Δt (2)
所述地震波在地层界面处的反射系数Ri通过公式3计算得到：
Ri＝(ρi+1vi+1-ρivi)/(ρi+1vi+1+ρivi) (3)
其中，ρi为第i地层的密度，vi为第i地层的剪切速度，ρi+1为第i+1地层的密度，vi+1为第i+1地层的剪切速度，hi为第i地层的地层厚度，Δt为时间采样间隔。
5.如权利要求1所述的基于地震脉冲响应谱测定场地卓越周期的方法，其特征在于，所述步骤5中各个土层界面的次级地震波时程ui d和通过递推公式4-9计算得到：
u0 d(j)＝-u1 p(j-nt1) (4)
u1 d(j)＝(1+R1)u0 d(j-nt1)-R1u2 p(j-nt2) (5)
u1 p(j)＝R1u0 d(j-nt1)+(1-R1)u2 p(j-nt2) (6)
………………………………
ui d(j)＝(1+Ri)ui-1 d(j-nti)-Riui+1 p(j-nti+1) (7)
ui p(j)＝Riui-1 d(j-nti)+(1-Ri)ui+1 p(j-nti+1) (8)
………………………………
un p(j)＝Rnun d(j-ntn) (9)
其中，j＝1，2，......，nt，un p(1)＝1，u0 d，u1 d，u1 p，……，ui d，是关于j的函数，分别表示地震作用下各土层之间的界面0，界面1，……，界面i，……，界面n产生的次级地震波时程；上标d，p分别表示向下，向上传播的次级波；i表示界面编号，j表示时间样点序号。
6.如权利要求5所述的基于地震脉冲响应谱测定场地卓越周期的方法，其特征在于，所述步骤6中脉冲地震动时程u0通过公式10的函数关系得到：
u0＝2u0 d(j) (10)
其中，j＝1，2，...，nt。
CN201810184529.1A 2018-03-06 2018-03-06 基于地震脉冲响应谱测定场地卓越周期的方法 Active CN108828661B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201810184529.1A CN108828661B (zh) 2018-03-06 2018-03-06 基于地震脉冲响应谱测定场地卓越周期的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201810184529.1A CN108828661B (zh) 2018-03-06 2018-03-06 基于地震脉冲响应谱测定场地卓越周期的方法

Publications (2)

Publication Number Publication Date
CN108828661A true CN108828661A (zh) 2018-11-16
CN108828661B CN108828661B (zh) 2019-10-25

Family

ID=64154031

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201810184529.1A Active CN108828661B (zh) 2018-03-06 2018-03-06 基于地震脉冲响应谱测定场地卓越周期的方法

Country Status (1)

Country Link
CN (1) CN108828661B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110244355A (zh) * 2019-07-25 2019-09-17 西南交通大学 一种基于震源断层模型的脉冲地震动模拟方法
CN110609323A (zh) * 2019-09-24 2019-12-24 西安理工大学 一种基于地层信息的场地卓越周期的精简计算方法
CN112444875A (zh) * 2020-10-28 2021-03-05 西安理工大学 一种获取场地卓越周期的精确解的方法
CN115903017A (zh) * 2022-10-17 2023-04-04 西南交通大学 一种近断层脉冲型地震动参数化方法

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050120894A1 (en) * 2003-11-14 2005-06-09 Fischer & Krecke Gmbh & Co. Method and apparatus for shifting a doctor blade
US20130286779A1 (en) * 2012-04-30 2013-10-31 Conocophillips Company Quasi-impulsive displacement source
CN105421609A (zh) * 2016-01-14 2016-03-23 同济大学 半主动变刚度变阻尼混合阻尼器
CN107119959A (zh) * 2017-06-30 2017-09-01 安徽理工大学 组合式单摆式调谐质量阻尼器结构系统及施工方法
CN108957537A (zh) * 2018-06-20 2018-12-07 西安理工大学 倾斜地震波作用下场地卓越周期的计算方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20050120894A1 (en) * 2003-11-14 2005-06-09 Fischer & Krecke Gmbh & Co. Method and apparatus for shifting a doctor blade
US20130286779A1 (en) * 2012-04-30 2013-10-31 Conocophillips Company Quasi-impulsive displacement source
CN105421609A (zh) * 2016-01-14 2016-03-23 同济大学 半主动变刚度变阻尼混合阻尼器
CN107119959A (zh) * 2017-06-30 2017-09-01 安徽理工大学 组合式单摆式调谐质量阻尼器结构系统及施工方法
CN108957537A (zh) * 2018-06-20 2018-12-07 西安理工大学 倾斜地震波作用下场地卓越周期的计算方法

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
王家全，等: "场地土卓越周期变化的试验分析", 《广西大学学报（自然科学版）》 *
范留明，等: "考虑基岩影响的场地自振特性", 《西安理工大学学报》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110244355A (zh) * 2019-07-25 2019-09-17 西南交通大学 一种基于震源断层模型的脉冲地震动模拟方法
CN110244355B (zh) * 2019-07-25 2021-06-08 西南交通大学 一种基于震源断层模型的脉冲地震动模拟方法
CN110609323A (zh) * 2019-09-24 2019-12-24 西安理工大学 一种基于地层信息的场地卓越周期的精简计算方法
CN110609323B (zh) * 2019-09-24 2022-02-22 西安理工大学 一种基于地层信息的场地卓越周期的精简计算方法
CN112444875A (zh) * 2020-10-28 2021-03-05 西安理工大学 一种获取场地卓越周期的精确解的方法
CN112444875B (zh) * 2020-10-28 2023-10-03 西安理工大学 一种获取场地卓越周期的精确解的方法
CN115903017A (zh) * 2022-10-17 2023-04-04 西南交通大学 一种近断层脉冲型地震动参数化方法
CN115903017B (zh) * 2022-10-17 2024-03-15 西南交通大学 一种近断层脉冲型地震动参数化方法

Also Published As

Publication number Publication date
CN108828661B (zh) 2019-10-25

Similar Documents

Publication Publication Date Title
CN108828661B (zh) 基于地震脉冲响应谱测定场地卓越周期的方法
CN102373923B (zh) 一种储层识别方法
WO2017024523A1 (zh) 一种射线弹性参数的反演方法
US11650342B2 (en) Method for exploring passive source seismic frequency resonance
CN104730579B (zh) 一种基于表层横波速度反演的纵横波联合静校正方法
CN104931363B (zh) 节理岩体变形模量测试方法
CN103454686A (zh) 基于地层切片的小尺度沉积相进行储层预测的方法及系统
CN104570110B (zh) 一种基于纵横波匹配的多分量资料联合速度分析方法
CN110133715A (zh) 一种基于初至时差和波形叠加的微地震震源定位方法
CN105719433B (zh) 一种基于孔内地震波的超前预报方法
CN109085663A (zh) 一种致密砂岩储层层理缝识别方法
US20230324577A1 (en) Real-Time Microseismic Magnitude Calculation Method and Device Based on Deep Learning
CN105259581B (zh) 一种地震资料时深转换方法
CN106324682A (zh) 一种应用于冻土层地区的表层结构调查方法
CN103399345B (zh) 一种潜山裂缝分布的勘测方法与装置
CN107506556A (zh) 一种确定新鲜完整岩块声波纵波速度值的简便方法
CN113009576B (zh) 一种基于本征模态滤波的储层固有频率原位非接触式检测方法
CN108919340A (zh) 一种单孔法波速测试方法、装置及终端设备
CN103628866A (zh) 获取地层的地层压力系数的方法
CN102914797A (zh) 一种获得地层各向异性系数的方法及装置
CN106324661A (zh) 一种利用天然源面波进行岩体风化分带的方法
CN207232405U (zh) 一种高速公路隧道中夹岩隔墙的微震信息监测系统
CN113031070B (zh) 一种深度域合成地震记录的制作方法
CN102455440B (zh) Vti介质中地震波的走时计算方法
CN106094051B (zh) 一种花岗岩地层划分对比方法

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