CN110244355A - 一种基于震源断层模型的脉冲地震动模拟方法 - Google Patents
一种基于震源断层模型的脉冲地震动模拟方法 Download PDFInfo
- Publication number
- CN110244355A CN110244355A CN201910677592.3A CN201910677592A CN110244355A CN 110244355 A CN110244355 A CN 110244355A CN 201910677592 A CN201910677592 A CN 201910677592A CN 110244355 A CN110244355 A CN 110244355A
- Authority
- CN
- China
- Prior art keywords
- model
- pulse
- frequency domain
- response
- earthquake
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 34
- 230000000694 effects Effects 0.000 claims abstract description 30
- 238000001228 spectrum Methods 0.000 claims description 58
- 230000001133 acceleration Effects 0.000 claims description 24
- 238000004364 calculation method Methods 0.000 claims description 21
- 230000005540 biological transmission Effects 0.000 claims description 15
- 238000003325 tomography Methods 0.000 claims description 14
- 238000006073 displacement reaction Methods 0.000 claims description 9
- 238000012544 monitoring process Methods 0.000 claims description 9
- 238000000611 regression analysis Methods 0.000 claims description 6
- 230000003321 amplification Effects 0.000 claims description 3
- 238000012512 characterization method Methods 0.000 claims description 3
- 238000002591 computed tomography Methods 0.000 claims description 3
- 238000012937 correction Methods 0.000 claims description 3
- 239000000463 material Substances 0.000 claims description 3
- 238000012821 model calculation Methods 0.000 claims description 3
- 238000003199 nucleic acid amplification method Methods 0.000 claims description 3
- 238000007619 statistical method Methods 0.000 claims description 3
- 230000009466 transformation Effects 0.000 claims description 3
- 210000001367 artery Anatomy 0.000 claims description 2
- 210000003462 vein Anatomy 0.000 claims description 2
- 230000010354 integration Effects 0.000 claims 1
- 238000011160 research Methods 0.000 abstract description 8
- 238000013461 design Methods 0.000 abstract description 6
- 230000015572 biosynthetic process Effects 0.000 abstract description 5
- 238000003786 synthesis reaction Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 3
- 238000004445 quantitative analysis Methods 0.000 abstract description 2
- 238000009825 accumulation Methods 0.000 description 2
- 230000006378 damage Effects 0.000 description 2
- 230000007812 deficiency Effects 0.000 description 2
- 208000004860 Blast Crisis Diseases 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000004615 ingredient Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000001308 synthesis method Methods 0.000 description 1
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/24—Recording seismic data
-
- 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/282—Application of seismic models, synthetic seismograms
-
- 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
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
本发明公开了一种基于震源断层模型的脉冲地震动模拟方法。本发明基于脉冲地震动数据库，结合脉冲地震动衰减规律和建立的山体地形效应模型，借助地震动模拟技术可在研究区域内按照栅格生产脉冲地震动时程，使用脉冲定量化方法对脉冲地震动进行定量表达，提出脉冲地震动的人工重现方法，为相关研究领域提供丰富合理的人工脉冲地震动，奠定研究基础，为抗震设计规范的完善和震后灾害的评估等提供数据支持。
Description
技术领域
本发明涉及基于震源断层模型的脉冲地震动模拟技术领域，具体涉及一种基于震源断层模型的脉冲地震动模拟方法。
背景技术
相关资料表明，除部分低概率、摧毁性的飓风、洪灾外，地震是破坏力最强的自然灾害，以2016年为例，仅为数不多的自然巨灾就造成了1750亿美元的损失，其中飓风与地震是最主要的损失来源。我国是地震频发的国家，全境分布有约20条地震带，强震分布广，地震破坏严重。
地震诱发滑坡的近断层分布规律及与断层走向的特定关系说明近断层地震动的某些特性对滑坡的启动和分布发育起着非常重要的作用。因此，近断层地震动成为近二十多年来地震学和地震工程学中的研究热点。近断层地震动对边坡、结构的破坏作用主要在于其速度时程上具有明显脉冲，能够在短时间内释放大量能量，进而使得边坡、结构产生强烈振动，最终发生失稳破坏。国内外很多学者对近断层脉冲地震动相关科学问题进行了研究，如速度脉冲对边坡隧道、桥梁、塔架等的破坏作用，为抗震设计规范的完善和震后灾害的评估等提供了一定参考。
目前相关研究领域最为棘手的问题之一就是实测近断层地震动记录匮乏。通常采用合成地震动来弥补这一不足，显然合成地震动的合理性是结构抗震设计须考虑的首要问题。目前，主要有4种方法以等效替代近断层地震动：①等效脉冲模型，该方法未考虑地震动高频成分的影响，并不能真实反映近断层地震动的特征；②远场地震动叠加简单脉冲模型，因近断层地震动高频成分与远场地震动存在显著差异，以远场地震动代替近场地震动高频成分并不合理；③地震动记录分解法，尽管其可准确反映地震动不同成分对结构响应的影响，但在抗震设计中难以利用该方法进行地震动合成；④无脉冲近场地震动叠加简单脉冲模型，该方法以无脉冲近场地震动代替基于震源断层模型的脉冲地震动高频成分，但未考虑两者间差异。
发明内容
针对现有技术中的上述不足，本发明提供的一种基于震源断层模型的脉冲地震动模拟方法解决了基于震源断层模型的脉冲地震动模拟结果不准确的问题。
为了达到上述发明目的，本发明采用的技术方案为：一种基于震源断层模型的脉冲地震动模拟方法，包括以下步骤：
S1、创建目标断层，按照需求将目标断层在纵向上将目标断层划分为nl个纵向单元，在横向上将目标断层划分为nw个横向单元，令i和j的初始值均为1，i为纵向单元的序数，j为横向单元的序数；
S2、提取目标断层所有单元的地质信息；
所述地质信息包括单元滑动位移；
S3、对脉冲地震动数据库中的所有记录进行统计和回归分析，得到脉冲衰减路径效应模型；
S4、基于多尺度潜在地震危险性区划成果或设定具体发震断层和震级的特定地震，按照区域内历史地震事件进行统计分析，得到脉冲震源效应模型；
S5、通过单元滑动位移计算脉冲震源效应模型的频域值，并通过脉冲衰减路径效应模型的频域值、脉冲震源效应模型的频域值和基于震源断层模型的脉冲地震动响应模型的频域值计算断层单元的频域值；
S6、利用傅里叶逆变换方法将断层单元的频域值转化为时域值，并将时域值进行累加计算加速度值；
S7、判断j是否小于nw，若是则进入步骤S8，否则令j加1，并返回步骤S1；
S8、判断i是否小于nl，若是则输出每个断层单元的加速度值作为加速度时程，通过加速度时程模拟脉冲地震动，否则令i加1，并返回步骤S1。
进一步地：所述步骤S5中基于震源断层模型的脉冲地震动响应模型的建立方法为：
S51、初始化脉冲地震记录数据序号s＝1，模型序号n＝1；
S52、利用FLAC3D软件建立m个模型，设置模型参数，并在模型上设置监测点；
所述模型包括山体模型和自由场模型；
S53、在模型n中输入脉冲地震记录数据f(t)，分别计算各监测点及对应自由场的加速度时程，n＝1,2,3…m；
S54、将加速度时程经过傅里叶变换，得到傅里叶谱；
S55、通过傅里叶谱计算各监测点的传递函数；
S56、计算脉冲地震记录数据fs(t)的傅里叶谱F(ws)；
S57、利用传递函数和脉冲地震记录数据fs(t)的傅里叶谱F(ws)计算频域响应；
S58、将频域响应经过傅里叶变化得到时域响应，并将时域响应转化为反应谱；
所述反应谱包括山体反应谱和对应自由场反应谱；
S59、通过山体反应谱和对应自由场反应谱的比值计算山体反应谱比值；
S510、判断s是否等于脉冲地震记录数据的总数，若是则进入步骤S511，否则，令s加1，并返回步骤S56；
S511、判断n是否等于模型总数m，若是则进入步骤S512，否则，令n加1，并返回步骤S53；
S512、通过频谱关系对山体反应谱比值进行回归分析，建立基于震源断层模型的脉冲地震动响应模型。
进一步地：所述步骤S52中的模型参数包括材料粘聚力和内摩擦角。
进一步地：所述步骤S55中传递函数的计算公式为：
H(w,zr)＝G(w,zr)/F(w)
上式中，H(w,zr)为传递函数，G(w,zr)和F(w)均为傅里叶谱。
进一步地：所述步骤S57中频域响应的计算公式为：
G(ws,zr)＝H(w,zr)*F(ws)
上式中，G(ws,zr)为频域响应，H(w,zr)为传递函数，F(ws)为脉冲地震记录数据fs(t)对应的傅里叶谱。
进一步地：所述步骤S512中的频谱关系为：
y＝c0+c1x1+c2x2+c3x3+…+ckxk+ε
上式中，y为山体反应谱比值，x1、x2、x3、…、xk均为变量，c0、c1、c2、c3、…、ck均为变量系数，ε为常数。
进一步地：所述步骤S512中的基于震源断层模型的脉冲地震动响应模型具体为：
ln(Am1)＝c0+c1ln(T)+c2[ln(T)]2+c3θ+c4A0
ln(ym,s)＝c1Mwm,s+c2ln(rm,s)+c3xm,s+c4ln(Tm+csc)+c5[ln(Tm+csc)]2
+c6ψ+c7ψ2+c8A0m+c9θm+c10ln(Tm)+ck+ζ
上式中，ln(ym,s)为模型m在脉冲地震记录s下的反应谱放大比，Am1为谱比平均值，T为自振周期，θ为山体倾角，A0为山顶宽度，rm,s为震级与震源距的函数，ln(rm,s)＝ln(xm,s+exp(c1+c2Mwm,s))，xm,s为模型m在脉冲地震记录s下的震源距，Mwm,s为模型m在脉冲地震记录s下的地震震级，Tm为模型m的自振周期，A0m为模型m的山顶宽度，θm为模型m的山体倾角，csc为控制近场反应谱的常数，ψ为测点与山体高度比，ζ为模型残差，包括由不同脉冲地震记录引起的残差及模型参数引起的残差。
进一步地：所述自振周期T的计算公式为：
上式中，H为山体高度，fα为修正系数，fα＝1.5309-1/(0.69α+1.8836)，α为表征山体倾斜程度的参数，α＝2Hcotθ/A0，Cs为剪切波速。
进一步地：所述步骤S5中断层单元的频域值的计算公式为：
上式中，Aij(f)为断层单元的频域值，AL(f)为山体地形效应模型的频域值，为脉冲衰减路径效应模型的频域值，为脉冲震源效应模型的频域值；
脉冲震源效应模型的频域值的计算公式为：
上式中，k为系数，Sij为单元滑动位移。
进一步地：所述步骤S6中时域值的计算公式为：
上式中，aij(t)为时域值，Aij(f)为断层单元的频域值，f为输入脉冲地震记录数据，t为时间；
所述加速度值的计算公式为：
上式中，a(t)为加速度值，Δtij为各相邻单元间的时间间隔。。
本发明的有益效果为：本发明基于脉冲地震动数据库，结合脉冲地震动衰减规律和建立的山体地形效应模型，借助地震动模拟技术可在研究区域内按照栅格生产脉冲地震动时程，使用脉冲定量化方法对脉冲地震动进行定量表达，提出脉冲地震动的人工重现方法，为相关研究领域提供丰富合理的人工脉冲地震动，奠定研究基础，为抗震设计规范的完善和震后灾害的评估等提供数据支持。
附图说明
图1为本发明流程图；
图2为建立基于震源断层模型的脉冲地震动响应模型的方法流程图。
具体实施方式
下面对本发明的具体实施方式进行描述，以便于本技术领域的技术人员理解本发明，但应该清楚，本发明不限于具体实施方式的范围，对本技术领域的普通技术人员来讲，只要各种变化在所附的权利要求限定和确定的本发明的精神和范围内，这些变化是显而易见的，一切利用本发明构思的发明创造均在保护之列。
如图1所示，一种基于震源断层模型的脉冲地震动模拟方法，其特征在于，包括以下步骤：
S1、创建目标断层，按照需求将目标断层在纵向上将目标断层划分为nl个纵向单元，在横向上将目标断层划分为nw个横向单元，令i和j的初始值均为1，i为纵向单元的序数，j为横向单元的序数；
S2、提取目标断层所有单元的地质信息；
所述地质信息包括单元滑动位移；
S3、对脉冲地震动数据库中的所有记录进行统计和回归分析，得到脉冲衰减路径效应模型；
S4、基于多尺度潜在地震危险性区划成果或设定具体发震断层和震级的特定地震，按照区域内历史地震事件进行统计分析，得到脉冲震源效应模型；
S5、通过单元滑动位移计算脉冲震源效应模型的频域值，并通过脉冲衰减路径效应模型的频域值、脉冲震源效应模型的频域值和基于震源断层模型的脉冲地震动响应模型的频域值计算断层单元的频域值；
如图2所示，基于震源断层模型的脉冲地震动响应模型的建立方法为：
S51、初始化脉冲地震记录数据序号s＝1，模型序号n＝1；
S52、利用FLAC3D软件建立m个模型，设置模型参数，并在模型上设置监测点；模型参数包括材料粘聚力和内摩擦角。
所述模型包括山体模型和自由场模型；
S53、在模型n中输入脉冲地震记录数据f(t)，分别计算各监测点及对应自由场的加速度时程，n＝1,2,3…m；
S54、将加速度时程经过傅里叶变换，得到傅里叶谱；
S55、通过傅里叶谱计算各监测点的传递函数；
传递函数的计算公式为：
H(w,zr)＝G(w,zr)/F(w)
上式中，H(w,zr)为传递函数，G(w,zr)和F(w)均为傅里叶谱。
S56、计算脉冲地震记录数据fs(t)的傅里叶谱F(ws)；
S57、利用传递函数和脉冲地震记录数据fs(t)的傅里叶谱F(ws)计算频域响应；
频域响应的计算公式为：
G(ws,zr)＝H(w,zr)*F(ws)
上式中，G(ws,zr)为频域响应，H(w,zr)为传递函数，F(ws)为脉冲地震记录数据fs(t)对应的傅里叶谱。
S58、将频域响应经过傅里叶变化得到时域响应，并将时域响应转化为反应谱；
所述反应谱包括山体反应谱和对应自由场反应谱；
S59、通过山体反应谱和对应自由场反应谱的比值计算山体反应谱比值；
S510、判断s是否等于脉冲地震记录数据的总数，若是则进入步骤S511，否则，令s加1，并返回步骤S56；
S511、判断n是否等于模型总数m，若是则进入步骤S512，否则，令n加1，并返回步骤S53；
S512、通过频谱关系对山体反应谱比值进行回归分析，建立基于震源断层模型的脉冲地震动响应模型。
频谱关系为：
y＝c0+c1x1+c2x2+c3x3+…+ckxk+ε
上式中，y为山体反应谱比值，x1、x2、x3、…、xk均为变量，c0、c1、c2、c3、…、ck均为变量系数，ε为常数。
基于震源断层模型的脉冲地震动响应模型具体为：
ln(Am1)＝c0+c1ln(T)+c2[ln(T)]2+c3θ+c4A0
ln(ym,s)＝c1Mwm,s+c2ln(rm,s)+c3xm,s+c4ln(Tm+csc)+c5[ln(Tm+csc)]2
+c6ψ+c7ψ2+c8A0m+c9θm+c10ln(Tm)+ck+ζ
上式中，ln(ym,s)为模型m在脉冲地震记录s下的反应谱放大比，Am1为谱比平均值，T为自振周期，θ为山体倾角，A0为山顶宽度，rm,s为震级与震源距的函数，ln(rm,s)＝ln(xm,s+exp(c1+c2Mwm,s))，xm,s为模型m在脉冲地震记录s下的震源距，Mwm,s为模型m在脉冲地震记录s下的地震震级，Tm为模型m的自振周期，A0m为模型m的山顶宽度，θm为模型m的山体倾角，csc为控制近场反应谱的常数，ψ为测点与山体高度比，ζ为模型残差，包括由不同脉冲地震记录引起的残差及模型参数引起的残差。
自振周期T的计算公式为：
上式中，H为山体高度，fα为修正系数，fα＝1.5309-1/(0.69α+1.8836)，α为表征山体倾斜程度的参数，α＝2Hcotθ/A0，Cs为剪切波速。
断层单元的频域值的计算公式为：
上式中，Aij(f)为断层单元的频域值，AL(f)为山体地形效应模型的频域值，为脉冲衰减路径效应模型的频域值，为脉冲震源效应模型的频域值；
脉冲震源效应模型的频域值的计算公式为：
上式中，k为系数，Sij为单元滑动位移。
S6、利用傅里叶逆变换方法将断层单元的频域值转化为时域值，并将时域值进行累加计算加速度值；
时域值的计算公式为：
上式中，aij(t)为时域值，Aij(f)为断层单元的频域值，f为输入脉冲地震记录数据，t为时间；
所述加速度值的计算公式为：
上式中，a(t)为加速度值，Δtij为各相邻单元间的时间间隔。
S7、判断j是否小于nw，若是则进入步骤S8，否则令j加1，并返回步骤S1；
S8、判断i是否小于nl，若是则输出每个断层单元的加速度值作为加速度时程，通过加速度时程模拟脉冲地震动，否则令i加1，并返回步骤S1。
本发明基于脉冲地震动数据库，结合脉冲地震动衰减规律和建立的山体地形效应模型，借助地震动模拟技术可在研究区域内按照栅格生产脉冲地震动时程，使用脉冲定量化方法对脉冲地震动进行定量表达，提出脉冲地震动的人工重现方法，为相关研究领域提供丰富合理的人工脉冲地震动，奠定研究基础，为抗震设计规范的完善和震后灾害的评估等提供数据支持。
Claims (10)
1.一种基于震源断层模型的脉冲地震动模拟方法，其特征在于，包括以下步骤：
S1、创建目标断层，按照需求将目标断层在纵向上将目标断层划分为nl个纵向单元，在横向上将目标断层划分为nw个横向单元，令i和j的初始值均为1，i为纵向单元的序数，j为横向单元的序数；
S2、提取目标断层所有单元的地质信息；
所述地质信息包括单元滑动位移；
S3、对脉冲地震动数据库中的所有记录进行统计和回归分析，得到脉冲衰减路径效应模型；
S4、基于多尺度潜在地震危险性区划成果或设定具体发震断层和震级的特定地震，按照区域内历史地震事件进行统计分析，得到脉冲震源效应模型；
S5、通过单元滑动位移计算脉冲震源效应模型的频域值，并通过脉冲衰减路径效应模型的频域值、脉冲震源效应模型的频域值和脉冲地震动响应模型的频域值计算断层单元的频域值；
S6、利用傅里叶逆变换方法将断层单元的频域值转化为时域值，并将时域值进行累加计算加速度值；
S7、判断j是否小于nw，若是则进入步骤S8，否则令j加1，并返回步骤S1；
S8、判断i是否小于nl，若是则输出每个断层单元的加速度值作为加速度时程，通过加速度时程模拟脉冲地震动，否则令i加1，并返回步骤S1。
2.根据权利要求1所述的基于震源断层模型的脉冲地震动模拟方法，其特征在于，所述步骤S5中基于震源断层模型的脉冲地震动响应模型的建立方法为：
S51、初始化脉冲地震记录数据序号s＝1，模型序号n＝1；
S52、利用FLAC3D软件建立m个模型，设置模型参数，并在模型上设置监测点；
所述模型包括山体模型和自由场模型；
S53、在模型n中输入脉冲地震记录数据f(t)，分别计算各监测点及对应自由场的加速度时程，n＝1,2,3…m；
S54、将加速度时程经过傅里叶变换，得到傅里叶谱；
S55、通过傅里叶谱计算各监测点的传递函数；
S56、计算脉冲地震记录数据fs(t)的傅里叶谱F(ws)；
S57、利用传递函数和脉冲地震记录数据fs(t)的傅里叶谱F(ws)计算频域响应；
S58、将频域响应经过傅里叶变化得到时域响应，并将时域响应转化为反应谱；
所述反应谱包括山体反应谱和对应自由场反应谱；
S59、通过山体反应谱和对应自由场反应谱的比值计算山体反应谱比值；
S510、判断s是否等于脉冲地震记录数据的总数，若是则进入步骤S511，否则，令s加1，并返回步骤S56；
S511、判断n是否等于模型总数m，若是则进入步骤S512，否则，令n加1，并返回步骤S53；
S512、通过频谱关系对山体反应谱比值进行回归分析，建立基于震源断层模型的脉冲地震动响应模型。
3.根据权利要求2所述的基于震源断层模型的脉冲地震动模拟方法，其特征在于，所述步骤S52中的模型参数包括材料粘聚力和内摩擦角。
4.根据权利要求2所述的基于震源断层模型的脉冲地震动模拟方法，其特征在于，所述步骤S55中传递函数的计算公式为：
H(w,zr)＝G(w,zr)/F(w)
上式中，H(w,zr)为传递函数，G(w,zr)和F(w)均为傅里叶谱。
5.根据权利要求2所述的基于震源断层模型的脉冲地震动模拟方法，其特征在于，所述步骤S57中频域响应的计算公式为：
G(ws,zr)＝H(w,zr)*F(ws)
上式中，G(ws,zr)为频域响应，H(w,zr)为传递函数，F(ws)为脉冲地震记录数据fs(t)对应的傅里叶谱。
6.根据权利要求2所述的基于震源断层模型的脉冲地震动模拟方法，其特征在于，所述步骤S512中的频谱关系为：
y＝c0+c1x1+c2x2+c3x3+…+ckxk+ε
上式中，y为山体反应谱比值，x1、x2、x3、…、xk均为变量，c0、c1、c2、c3、…、ck均为变量系数，ε为常数。
7.根据权利要求6所述的基于震源断层模型的脉冲地震动模拟方法，其特征在于，所述步骤S512中的基于震源断层模型的脉冲地震动响应模型具体为：
ln(Am1)＝c0+c1ln(T)+c2[ln(T)]2+c3θ+c4A0
ln(ym,s)＝c1Mwm,s+c2ln(rm,s)+c3xm,s+c4ln(Tm+csc)+c5[ln(Tm+csc)]2
+c6ψ+c7ψ2+c8A0m+c9θm+c10ln(Tm)+ck+ζ
上式中，ln(ym,s)为模型m在脉冲地震记录s下的反应谱放大比，Am1为谱比平均值，T为自振周期，θ为山体倾角，A0为山顶宽度，rm,s为震级与震源距的函数，ln(rm,s)＝ln(xm,s+exp(c1+c2Mwm,s))，xm,s为模型m在脉冲地震记录s下的震源距，Mwm,s为模型m在脉冲地震记录s下的地震震级，Tm为模型m的自振周期，A0m为模型m的山顶宽度，θm为模型m的山体倾角，csc为控制近场反应谱的常数，ψ为测点与山体高度比，ζ为模型残差，包括由不同脉冲地震记录引起的残差及模型参数引起的残差。
8.根据权利要求7所述的基于震源断层模型的脉冲地震动模拟方法，其特征在于，所述自振周期T的计算公式为：
上式中，H为山体高度，fα为修正系数，fα＝1.5309-1/(0.69α+1.8836)，α为表征山体倾斜程度的参数，α＝2Hcotθ/A0，Cs为剪切波速。
9.根据权利要求1所述的基于震源断层模型的脉冲地震动模拟方法，其特征在于，所述步骤S5中断层单元的频域值的计算公式为：
上式中，Aij(f)为断层单元的频域值，AL(f)为山体地形效应模型的频域值，为脉冲衰减路径效应模型的频域值，为脉冲震源效应模型的频域值；
脉冲震源效应模型的频域值的计算公式为：
上式中，k为系数，Sij为单元滑动位移。
10.根据权利要求1所述的基于震源断层模型的脉冲地震动模拟方法，其特征在于，所述步骤S6中时域值的计算公式为：
上式中，aij(t)为时域值，Aij(f)为断层单元的频域值，f为输入脉冲地震记录数据，t为时间；
所述加速度值的计算公式为：
上式中，a(t)为加速度值，Δtij为各相邻单元间的时间间隔。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910677592.3A CN110244355B (zh) | 2019-07-25 | 2019-07-25 | 一种基于震源断层模型的脉冲地震动模拟方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201910677592.3A CN110244355B (zh) | 2019-07-25 | 2019-07-25 | 一种基于震源断层模型的脉冲地震动模拟方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN110244355A true CN110244355A (zh) | 2019-09-17 |
| CN110244355B CN110244355B (zh) | 2021-06-08 |
Family
ID=67893478
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201910677592.3A Active CN110244355B (zh) | 2019-07-25 | 2019-07-25 | 一种基于震源断层模型的脉冲地震动模拟方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN110244355B (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115236732A (zh) * | 2022-07-19 | 2022-10-25 | 福州大学 | 一种基于相位谱重构的近断层地震动拟合方法 |
| CN115826053A (zh) * | 2022-11-30 | 2023-03-21 | 西南交通大学 | 近走滑活动断层岩土工程灾变评价地震动确定方法 |
| CN115903017A (zh) * | 2022-10-17 | 2023-04-04 | 西南交通大学 | 一种近断层脉冲型地震动参数化方法 |
| CN116609823A (zh) * | 2023-04-13 | 2023-08-18 | 西南交通大学 | 考虑断层和地震动特性的概率性地震滑坡危险性评价方法 |
Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1144771A (ja) * | 1997-07-28 | 1999-02-16 | Tokyu Constr Co Ltd | 地震動予測方法及び地震動予測のプログラムを記録した記録媒体 |
| CN1493005A (zh) * | 2001-01-25 | 2004-04-28 | 维斯特恩格科地震控股有限公司 | 海洋地震数据处理方法和地震测量方法 |
| WO2007037721A1 (fr) * | 2005-09-28 | 2007-04-05 | Schlumberger Canada Limited | Procede de monitorage sismique d'une faille hydraulique d'une couche productrice (et variantes) |
| CN103440369A (zh) * | 2013-08-14 | 2013-12-11 | 中国石油集团东方地球物理勘探有限责任公司 | 脉冲震源的优化方法及装置 |
| CN103926621A (zh) * | 2014-05-06 | 2014-07-16 | 中国地震局工程力学研究所 | 两阶段残差分析建立地震动衰减关系的方法 |
| WO2014205162A1 (en) * | 2013-06-21 | 2014-12-24 | Schlumberger Canada Limited | Determining change in permeability caused by a hydraulic fracture in reservoirs |
| CN104462814A (zh) * | 2014-12-03 | 2015-03-25 | 沈阳工业大学 | 一种近地表地震动模拟的网格分级方法 |
| CN106461803A (zh) * | 2014-03-20 | 2017-02-22 | 施蓝姆伯格技术公司 | 从时间分布式引爆气枪阵列数据重建脉冲震源地震数据 |
| CN106646602A (zh) * | 2016-12-30 | 2017-05-10 | 中科宇图科技股份有限公司 | 一种基于多种震源模型的震后快速生成震动图的方法 |
| CN107133414A (zh) * | 2017-05-18 | 2017-09-05 | 中国科学院武汉岩土力学研究所 | 一种块体地震动力响应的分析方法 |
| CN107480391A (zh) * | 2017-08-24 | 2017-12-15 | 中铁二院工程集团有限责任公司 | 基于数据驱动的近断层非平稳地震动模拟方法 |
| CN107807388A (zh) * | 2017-11-02 | 2018-03-16 | 辽宁工程技术大学 | 一种基于多普勒效应的地震断层滑动速度计算方法 |
| CN108169794A (zh) * | 2017-12-14 | 2018-06-15 | 中铁二院工程集团有限责任公司 | 近断层人工地震波与规范反应谱匹配的调整方法 |
| CN108828661A (zh) * | 2018-03-06 | 2018-11-16 | 西安理工大学 | 基于地震脉冲响应谱测定场地卓越周期的方法 |
| CN109188524A (zh) * | 2018-10-26 | 2019-01-11 | 辽宁工程技术大学 | 考虑介质品质因子的地震中场地反应计算方法 |
| CN109375253A (zh) * | 2018-12-13 | 2019-02-22 | 中国地震局地球物理研究所 | 基于全部发震构造最大可信地震的地震动参数评价方法 |
| CN109375252A (zh) * | 2018-12-13 | 2019-02-22 | 中国地震局地球物理研究所 | 考虑不同发震构造最大可信地震的地震动参数评价方法 |
| CN109444879A (zh) * | 2018-10-19 | 2019-03-08 | 西南交通大学 | 一种DInSAR近断层同震形变场提取方法 |
-
2019
- 2019-07-25 CN CN201910677592.3A patent/CN110244355B/zh active Active
Patent Citations (18)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JPH1144771A (ja) * | 1997-07-28 | 1999-02-16 | Tokyu Constr Co Ltd | 地震動予測方法及び地震動予測のプログラムを記録した記録媒体 |
| CN1493005A (zh) * | 2001-01-25 | 2004-04-28 | 维斯特恩格科地震控股有限公司 | 海洋地震数据处理方法和地震测量方法 |
| WO2007037721A1 (fr) * | 2005-09-28 | 2007-04-05 | Schlumberger Canada Limited | Procede de monitorage sismique d'une faille hydraulique d'une couche productrice (et variantes) |
| WO2014205162A1 (en) * | 2013-06-21 | 2014-12-24 | Schlumberger Canada Limited | Determining change in permeability caused by a hydraulic fracture in reservoirs |
| CN103440369A (zh) * | 2013-08-14 | 2013-12-11 | 中国石油集团东方地球物理勘探有限责任公司 | 脉冲震源的优化方法及装置 |
| CN106461803A (zh) * | 2014-03-20 | 2017-02-22 | 施蓝姆伯格技术公司 | 从时间分布式引爆气枪阵列数据重建脉冲震源地震数据 |
| CN103926621A (zh) * | 2014-05-06 | 2014-07-16 | 中国地震局工程力学研究所 | 两阶段残差分析建立地震动衰减关系的方法 |
| CN104462814A (zh) * | 2014-12-03 | 2015-03-25 | 沈阳工业大学 | 一种近地表地震动模拟的网格分级方法 |
| CN106646602A (zh) * | 2016-12-30 | 2017-05-10 | 中科宇图科技股份有限公司 | 一种基于多种震源模型的震后快速生成震动图的方法 |
| CN107133414A (zh) * | 2017-05-18 | 2017-09-05 | 中国科学院武汉岩土力学研究所 | 一种块体地震动力响应的分析方法 |
| CN107480391A (zh) * | 2017-08-24 | 2017-12-15 | 中铁二院工程集团有限责任公司 | 基于数据驱动的近断层非平稳地震动模拟方法 |
| CN107807388A (zh) * | 2017-11-02 | 2018-03-16 | 辽宁工程技术大学 | 一种基于多普勒效应的地震断层滑动速度计算方法 |
| CN108169794A (zh) * | 2017-12-14 | 2018-06-15 | 中铁二院工程集团有限责任公司 | 近断层人工地震波与规范反应谱匹配的调整方法 |
| CN108828661A (zh) * | 2018-03-06 | 2018-11-16 | 西安理工大学 | 基于地震脉冲响应谱测定场地卓越周期的方法 |
| CN109444879A (zh) * | 2018-10-19 | 2019-03-08 | 西南交通大学 | 一种DInSAR近断层同震形变场提取方法 |
| CN109188524A (zh) * | 2018-10-26 | 2019-01-11 | 辽宁工程技术大学 | 考虑介质品质因子的地震中场地反应计算方法 |
| CN109375253A (zh) * | 2018-12-13 | 2019-02-22 | 中国地震局地球物理研究所 | 基于全部发震构造最大可信地震的地震动参数评价方法 |
| CN109375252A (zh) * | 2018-12-13 | 2019-02-22 | 中国地震局地球物理研究所 | 考虑不同发震构造最大可信地震的地震动参数评价方法 |
Non-Patent Citations (7)
| Title |
|---|
| ZHIWANG CHANG 等: ""An improved energy‐based approach for selecting pulse‐like ground motions"", 《EARTHQUAKE ENGINEERING & STRUCTURAL DYNAMICS》 * |
| 任叶飞 等: ""运用广义反演法研究汶川地震动效应"", 《土木工程学报》 * |
| 宋健 等: ""近断层速度脉冲地震动对边坡滑移的影响分析"", 《岩石力学与工程学报》 * |
| 王宏伟 等: ""地震动震源、路径及场地效应分析的广义反演方法研究进展"", 《地球物理学进展》 * |
| 王海云 等: ""近断层强地震动预测中的有限断层震源模型"", 《地球科学——中国地质大学学报》 * |
| 罗全波 等: ""近断层速度脉冲与震源机制的关系浅析"", 《震灾防御技术》 * |
| 高平军: ""梯形山体的地形效应模型研究"", 《中国硕士学位论文全文数据库 工程科技II辑》 * |
Cited By (8)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN115236732A (zh) * | 2022-07-19 | 2022-10-25 | 福州大学 | 一种基于相位谱重构的近断层地震动拟合方法 |
| CN115236732B (zh) * | 2022-07-19 | 2024-06-04 | 福州大学 | 一种基于相位谱重构的近断层地震动拟合方法 |
| CN115903017A (zh) * | 2022-10-17 | 2023-04-04 | 西南交通大学 | 一种近断层脉冲型地震动参数化方法 |
| CN115903017B (zh) * | 2022-10-17 | 2024-03-15 | 西南交通大学 | 一种近断层脉冲型地震动参数化方法 |
| CN115826053A (zh) * | 2022-11-30 | 2023-03-21 | 西南交通大学 | 近走滑活动断层岩土工程灾变评价地震动确定方法 |
| CN115826053B (zh) * | 2022-11-30 | 2023-08-18 | 西南交通大学 | 近走滑活动断层岩土工程灾变评价地震动确定方法 |
| CN116609823A (zh) * | 2023-04-13 | 2023-08-18 | 西南交通大学 | 考虑断层和地震动特性的概率性地震滑坡危险性评价方法 |
| CN116609823B (zh) * | 2023-04-13 | 2024-01-30 | 西南交通大学 | 考虑断层和地震动特性的概率性地震滑坡危险性评价方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN110244355B (zh) | 2021-06-08 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110244355A (zh) | 一种基于震源断层模型的脉冲地震动模拟方法 | |
| Régnier et al. | International benchmark on numerical simulations for 1D, nonlinear site response (PRENOLIN): Verification phase based on canonical cases | |
| Koketsu et al. | Finite-element simulation of seismic ground motion with a voxel mesh | |
| CN103221842A (zh) | 用于利用相位外推的数据反演的系统和方法 | |
| Bolt | Seismic input motions for nonlinear structural analysis | |
| CN108983285A (zh) | 一种基于矩张量的多种震源波场模拟方法及装置 | |
| Ferraro et al. | Site effects evaluation in Catania (Italy) by means 0f 1-D numerical analysis | |
| Sgambi et al. | Monte Carlo simulation for seismic analysis of a long span suspension bridge | |
| Castro et al. | Stochastic simulation of strong-motion records from the 26 September 1997 (M w 6), Umbria-Marche (central Italy) earthquake | |
| Zollo et al. | Modelling of ground acceleration in the near source range: the case of 1976, Friuli earthquake (M= 6.5), northern Italy | |
| Wu et al. | Seismic monitoring of super high‐rise building using ambient noise with dense seismic array | |
| Tang et al. | The Hualien large-scale seismic test for soil-structure interaction research | |
| Jiang et al. | Seismology-based hybrid ground motion prediction models of PGA in Sichuan, China | |
| Albrecht et al. | Feedback mechanisms controlling Antarctic glacial-cycle dynamics simulated with a coupled ice sheet–solid Earth model | |
| CN113281488B (zh) | 一种反倾岩质边坡模型试验的模拟地应力试验装置及方法 | |
| Bandini et al. | The shaking table with laminar box of the University of Messina | |
| Yan | Seismic deformation analysis of earth dams: a simplified method | |
| Atkinson et al. | Reassessing the New Madrid seismic zone | |
| Lemos et al. | Numerical and analytical modelling of the seismic behaviour of shallow tunnels | |
| Dashdondog et al. | Determining 1D velocity model from local earthquake data in the South Hangay region, central Mongolia | |
| Carnevale et al. | Generation of non-synchronous accelerograms for evaluate the seismic bridge response, including local site amplification. | |
| THESSALONIKI | Copyright© 2018 A. Kiratzi, Z. Roumelioti, Ch. Benetatos, N. Theodulidis, A. Savvaidis, A. Panou, IN Tziavos, P. Savvaidis, Th. Hatzigogos, S. | |
| Datta | Effect of non-linear soil structure interaction on the seismic response of critical infrastructure facilities | |
| Dong et al. | Soil Dynamic Constitutive Considering Post-Liquefaction Deformation and Reversible Pore-Water Pressure. Sustainability 2022, 14, 16512 | |
| Pradono et al. | Strong Ground Motion by September 30, 2009 Pariaman Earthquake and Damage To Large Scale Buildings |
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 |