CN109375253A - 基于全部发震构造最大可信地震的地震动参数评价方法 - Google Patents
基于全部发震构造最大可信地震的地震动参数评价方法 Download PDFInfo
- Publication number
- CN109375253A CN109375253A CN201811528301.6A CN201811528301A CN109375253A CN 109375253 A CN109375253 A CN 109375253A CN 201811528301 A CN201811528301 A CN 201811528301A CN 109375253 A CN109375253 A CN 109375253A
- Authority
- CN
- China
- Prior art keywords
- earthquake
- site
- parameter
- whole
- scheme
- 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
- 238000011156 evaluation Methods 0.000 title claims abstract description 19
- 230000001133 acceleration Effects 0.000 claims abstract description 51
- 238000000034 method Methods 0.000 claims abstract description 35
- 238000001228 spectrum Methods 0.000 claims abstract description 35
- 230000015572 biosynthetic process Effects 0.000 claims abstract description 33
- 238000003786 synthesis reaction Methods 0.000 claims abstract description 33
- 230000004044 response Effects 0.000 claims abstract description 25
- 238000004364 calculation method Methods 0.000 claims abstract description 14
- 238000004458 analytical method Methods 0.000 claims abstract description 13
- 238000012916 structural analysis Methods 0.000 claims abstract description 6
- 238000013461 design Methods 0.000 claims description 16
- 238000009826 distribution Methods 0.000 claims description 11
- 238000004088 simulation Methods 0.000 claims description 8
- 230000000694 effects Effects 0.000 claims description 7
- 238000007619 statistical method Methods 0.000 claims description 7
- 238000010835 comparative analysis Methods 0.000 claims description 3
- 238000011835 investigation Methods 0.000 claims description 3
- 230000007246 mechanism Effects 0.000 claims description 3
- 230000008569 process Effects 0.000 description 5
- 238000003325 tomography Methods 0.000 description 5
- 230000019771 cognition Effects 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 230000008901 benefit Effects 0.000 description 2
- 230000035939 shock Effects 0.000 description 2
- 230000003321 amplification Effects 0.000 description 1
- 230000002146 bilateral effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 230000006854 communication Effects 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000009510 drug design Methods 0.000 description 1
- BTCSSZJGUNDROE-UHFFFAOYSA-N gamma-aminobutyric acid Chemical compound NCCCC(O)=O BTCSSZJGUNDROE-UHFFFAOYSA-N 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003199 nucleic acid amplification method 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/01—Measuring or predicting earthquakes
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01M—TESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
- G01M7/00—Vibration-testing of structures; Shock-testing of structures
-
- 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
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/10—Aspects of acoustic signal generation or detection
- G01V2210/12—Signal generation
- G01V2210/121—Active source
-
- 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/66—Subsurface modeling
- G01V2210/667—Determining confidence or uncertainty in parameters
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请涉及一种基于全部发震构造最大可信地震的地震动参数评价方法，包括：基于区域构造和区域地震地质资料，确定地震动模拟震源设置方案，并对其进行不确定性分析；确定地震动模拟方法，并对模拟结果影响较大的参数取值进行分析；考虑全部震源影响，对地震动模拟方案进行设计；采用随机有限断层法进行场址地震动模拟，获得全部方案的加速度时程、加速度峰值和不同周期的加速度反应谱值；对考虑全部震源参与设计的每个方案的地震动模拟结果进行分析，并对场址地震动参数进行综合评价；对全部模拟方案的凹凸体模型和随机模型计算得到的参数统计值进行分析，取凹凸体模型和随机模型计算结果的较大值作为场址的参数；根据场址地震动参数，并考虑结构特性，对用于结构动力分析的加速度时程进行模拟。
Description
技术领域
本申请涉及一种基于全部发震构造最大可信地震的地震动参数评价方法，适用于地震工程的技术领域。
背景技术
最大可信地震是指目前已知或假定的地壳构造框架下，有理由相信的已确定的断层或在地理上界定的地震构造区内发生的最大地震，即代表一个合理设计的结构所能承受地震动的最高水平，是抗御超越设计地震动的最大能力。对于重大工程，不仅要求能够满足目前抗震设计中最大设计地震的单一水准要求，还需要满足最大可信地震的要求。
对于重大工程，对其能产生重要影响的一般为离场址较近的大断裂，在确定重要结构的地震输入时，还需要解决以下几个问题：1.将近断裂大震震源作为地表的一个点源看待，明显是不合理的；2.对于地震地质环境较为复杂的场址，对场址有影响的近场大断裂可能不止一个；3.地震从震源-介质-地表传播过程是一个非常复杂的过程，对这个过程的模拟还存在着一些随机不确定性，以及在发震过程中对某些参数的认知不确定等问题。因此，要获得发震构造最大可信地震的场址地震动参数，首先需要合理解决这些问题。
发明内容
本方法将采用随机有限断层法计算场址地震动参数，来解决地震动模拟中复杂震源参数设置与计算效率的问题；采用多方案和多重权重系数的设计来解决随机不确定性和认知不确定性等问题，从而建立了一种适合于工程应用的基于发震构造最大可信地震的、考虑对场址有影响的全部震源的地震动参数评价方法。
本申请涉及一种基于全部发震构造最大可信地震的地震动参数评价方法，包括以下步骤：
(1)基于区域构造和区域地震地质资料，确定地震动模拟震源设置方案，并对其进行不确定性分析；
(2)确定地震动模拟方法，并对模拟结果影响较大的参数取值进行分析；
(3)考虑全部震源影响，对地震动模拟方案进行设计；
(4)采用随机有限断层法进行场址的地震动模拟，获得全部方案的加速度时程、加速度峰值和不同周期的加速度反应谱值；
(5)对考虑全部震源参与设计的每个方案的地震动模拟结果进行分析，并对场址地震动参数进行综合评价；
(6)对全部模拟方案的凹凸体模型和随机模型计算得到的参数统计值进行分析，取凹凸体模型和随机模型计算结果的较大值作为场址的参数；
(7)根据场址地震动参数，并考虑结构特性，对用于结构动力分析的加速度时程进行模拟。
优选地，步骤(3)包括：
考虑全部震源对场址的影响，采用凹凸模型设置地震动模拟方案，并计算每种方案总体的权重系数；
采用随机模型设置用于对比分析的基于全部震源的地震动模拟方案，并计算每种方案的权重系数；
优选地，步骤(5)包括：
分析全部震源参与设计的每个方案计算得到结果，将每个方案的凹凸体模型和随机模型计算得到的地震动加速度峰值和不同周期的加速度谱值，分别按照从小到大进行排列，并根据其对应的加权系数，形成场址地震动参数排列表；
基于全部震源设计方案计算得到的场址地震动参数排列表，对全部方案凹凸体模型计算得到的地震动加速度峰值和不同周期点加速度反应谱值进行统计分析，给出最小值、50％分位数值、均值、84％分位数值、95％分位数值及最大值；
基于全部震源设计方案计算得到的场址地震动参数排列表，对全部方案随机模型计算得到的地震动加速度峰值和不同周期点加速度反应谱值进行统计分析，给出最小值、50％分位数值、均值、84％分位数值、95％分位数值及最大值。
优选地，步骤(1)中还包括：
(1)基于研究区域构造背景，确定对场址有重要影响的发震构造，并完成以下分析：
A.依据发震构造最新活动时代及地震震中分布特点，确定发震构造的震级上限；
B.依据地震、地质资料、野外调查、并结合区域深部构造背景，确定发震构造的位置，对不能确定明确位置的断裂，采用多个位置方案，并按确定位置的资料的可信度及对场址的影响水平分析，设置加权系数；
C.依据活动断裂的活动性(活动时代、活动速率)及大震复发周期等，对影响场址的不同发震断裂设置加权系数；
(2)依据研究区域地震、地质、小震精定位、震源机制解等资料，对发震断层的参数进行分析，确定地震动模拟震源方案，对于具有不确定性的参数，设置多个模拟方案，并设置各个方案的加权系数。
优选地，步骤(7)中的模拟方法为：
A.以场址地震活动环境为基础，挑选与场址加速度反应谱符合较好的实际地震加速度时程为种子时程，然后基于目标反应谱进行调整；
B.以场址地震活动环境为基础，建立地震动频率非平稳特性参数，人工拟合得到符合于目标反应谱的地震加速度时程；或者
C.从模拟得到的加速度时程中挑选出与目标反应谱符合较好的时程。
根据本申请的一种基于全部发震构造最大可信地震的地震动参数评价方法，其具有以下技术优势：
1.本方法通过对每个震源设置加权系数，来考虑未来一定时间内不同震源的发震概率及其对评价结果的影响。
2.本方法所采用的地震动模拟方法，不仅可以描述更加符合真实情况的断层破裂源、考虑地震波路径衰减、场地衰减等因素，而且计算效率高，便于工程应用。
3.本方法采用多方案设计，能够对未知地震发生过程中的随机不确定性，以及在发震过程中对某些参数的认知不确定进行充分考虑。
4.本方法通过在每个分析节点上设置不同的权重系数来合理描述参数不确定的影响。
5.本方法能够给出多风险水平的地震动参数评价结果，便于实际工程按照结构的重要性进行决策。
附图说明
图1显示了根据本申请的考虑全部震源的地震动参数评价方法的示意图。
图2显示了对场址有重要影响的断裂位置的示例。
图3显示了随机有限断层法原理的示意图。
图4显示了应力降分布频度图。
图5显示了κ值随震中距离的分布图。
图6显示了一个震源地震动模拟凹凸体模型多方案设计的示意图。
图7显示了一个震源地震动模拟随机模型多方案设计的示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明白，下文中将结合附图对本发明的实施例进行详细说明。需要说明的是，在不冲突的情况下，本申请中的实施例及实施例中的特征可以相互任意组合。
根据本申请的基于全部发震构造最大可信地震的地震动参数评价方法，包括以下步骤：
(1)基于区域构造和区域地震地质资料，确定地震动模拟震源设置方案，并对其进行不确定性分析；
(2)确定地震动模拟方法，并对模拟结果影响较大的参数取值进行分析；
(3)考虑全部震源影响，对地震动模拟方案进行设计；
(4)采用随机有限断层法进行场址的地震动模拟，获得全部方案的加速度时程、加速度峰值和不同周期的加速度反应谱值；
(5)对考虑全部震源参与设计的每个方案的地震动模拟结果进行分析，并对场址地震动参数进行综合评价；
(6)对全部模拟方案的凹凸体模型和随机模型计算得到的参数统计值进行分析，取凹凸体模型和随机模型计算结果的较大值作为场址的参数；
(7)根据场址地震动参数，并考虑结构特性，对用于结构动力分析的加速度时程进行模拟。
下面依次对上述各步骤进行具体的说明。
步骤(1)：基于区域构造和区域地震地质资料，确定地震动模拟震源设置方案，并对其进行不确定性分析
地震动是由震源位置发生的破裂产生的地震波在地表引起的地表运动。震源设置影响到地震波传播过程中的能量衰减、断层破裂面的展布等因素。由于不同地区的地震地质条件、发震构造背景存在显著的差异，强震的发震位置具有较强的随机性，在对设定地震进行地震动模拟时，首先需要确定可能的发震断层，以及位于发震断层上的震源位置，进一步利用地震地质资料构建震源的有限断层模型来进行地震动模拟。因此，对设定地震的震源设置无疑是地震动模拟中一个重要的考虑因素。具体地，步骤(1)中还可以包括以下内容：
1.基于研究区域构造背景，确定对场址有重要影响的发震构造，并完成以下分析：
A.依据发震构造最新活动时代及地震震中分布特点，确定发震构造的震级上限；
B.依据地震、地质资料、野外调查、并结合区域深部构造背景，确定发震构造的位置，对不能确定明确位置的断裂，采用多个位置方案，并按确定位置的资料的可信度及对场址的影响水平分析，设置加权系数。例如，如图2所示，对场址有重要影响的断裂有3个，其中对断裂1设置了两个不同的位置；
C.依据活动断裂的活动性(例如，活动时代、活动速率等)及大震复发周期等，对影响场址的不同发震断裂设置加权系数。
2.依据研究区域地震、地质、小震精定位、震源机制解等资料，对发震断层的参数进行分析，确定地震动模拟震源方案，对于具有不确定性的参数，设置多个模拟方案，并设置各个方案的加权系数，如表1所示。不确定性的标准主要是基于对未知地震发生过程中的随机不确定性，以及在发震过程中对某些参数的认知不确定来决定。
下面以表1中的“破裂位置”参数为例，来说明如何确定多方案还是单方案。破裂位置的多方案是发震过程可能存在的多种随机可能性决定，而发震位置的确定基本对应于破裂方向的确定。比如如果破裂位置在断层中部，其破裂方向为双侧破裂；如果破裂位置在断层的端点，则破裂方向为单侧破裂。所以发震位置取多方案，而破裂方向取单方案。
表1震源参数设置方案
| 序号 | 参数名称 | 参数设置 | 说明 | 权重 |
| 1 | 断层尺度 | 可按断层性质及震级上限估计 | 多方案 | ■ |
| 2 | 走向 | 按地质资料及地震反演结果取值 | 单方案 | |
| 3 | 倾向 | 按地质资料及地震反演结果取值 | 单方案 | |
| 4 | 倾角 | 按地质资料及地震反演结果取值 | 多方案 | ■ |
| 5 | 运动性质 | 正断层、逆断层、走滑断层 | 单方案 | |
| 6 | 震源深度 | 按地震资料及其反演结果取值 | 单方案 | |
| 7 | 破裂位置 | 按照断层两端和中间位置设置 | 多方案 | ■ |
| 8 | 破裂方向 | 依据破裂位置设置单向或双向 | 单方案 | |
| 9 | 凹凸体模型 | 依据地震规模大小设置 | 多方案 | ■ |
步骤(2)：确定地震动模拟方法，并对模拟结果影响较大的参数取值进行分析
场址地震动模拟采用随机有限断层法进行，该方法属于一种半经验半理论的地震动模拟方法，不仅可以设置比较符合实际情况的震源参数，考虑距离衰减和场地效应等影响，而且计算工作量小，因此比较适合于工程应用。
如图3所示，该方法的主要思想是将发震断层划分成一系列子断层，然后将每个子断层视为点源，应用随机点源的方法，得到每个子源对于场址的影响。计算每个子源在场址产生的地震动，首先是在频域上综合考虑震源、传播路径以及场地效应的影响，拟合地震动傅立叶谱，然后转换到时域上，生成地震动时程。如果设一个地震矩为M0的点源产生的傅立叶谱为FA(M0,f,R)，它可以表示为震源谱S(M0,f)、距离衰减项P(R,f)、场地效应影响项G(f)以及地震动类型因子I(f)之间的乘积，即
FA(M0,f,R)＝S(M0,f)·P(R,f)·G(f)·I(f) (1)
式中：f是频率，R是场地距震中的距离，M0为地震矩。
对模拟结果影响较大的参数取值进行分析，主要包括对应力降和高频衰减因子kappa(κ)值进行分析。
应力降：该参数是影响地震动模拟结果的最重要的参数，因此应选用多种方法确定其取值及分布。可选择研究区内记录到近震震级ML≥4.0以上地震进行该区域应力降的反演，结合该区域其它方法得到的应力降的取值，获得应力降的取值范围及分布情况，并根据其分布情况，确定取值水平及各取值水平的加权系数。
如图2所示的研究区域，其应力降的主要取值范围集中在20bar-60bar之间，其频度的分布图如图4所示。根据Brune的震源谱模型，利用应力降分布计算出平均震源谱，并求出对应的平均应力降。在实际模拟中，为了充分考虑到可能发生的最危险情况，在应力降取值时对研究区域内的应力降取值范围按照平均应力降乘以12.5％的放大倍数，并按照10bar的间隔设置5个水平，即(25,35,45,55,65)bar进行分析，根据每个应力降取值与平均应力降的偏离程度，以公式(2)计算出其对应的权重系数。
权重系数＝1/exp(|取值-均值|/均值) (2)
高频衰减因子kappa(κ)值：利用研究区域内的测震记录，经过零线校正和滤波后，求出其加速度谱，在半对数坐标中计算S波的傅立叶谱，选择其高频衰减的起始频率和终止频率进行计算，将各个台站的计算结果通过一阶近似得到0km处的κ值，即为目标场地与传播路径无关的κ值。
图5显示了κ值随震中距离的分布图。如图5所示的研究区内，高频衰减系数κ可设为三个取值水平0.0246,0.0277和0.0308，即高频衰减参数在计算中通过拟合结果可以求出均值和误差，这里以均值(mean)和加减一倍方差(delta)来进行取值，即κ＝均值±方差，相应的加权系数分别取为(方差/均值，方差/均值，1-2*方差/均值)。
由于采用随机有限断层法进行地震动模拟时，采用的相位是随机的，所以对每个方案需要产生的地震加速度时程的样本数量需要进行分析。按照样本均值稳定，方差较小的原则，经分析最后选取的样本数量为30条，其相应的权重系数为1/30。
步骤(3)：考虑全部震源影响，对地震动模拟方案进行设计；
根据步骤(1)和步骤(2)对每个参数设置的方案和相应的加权系数，设计地震动模拟方案。本申请按照考虑全部震源的方式进行考虑。
A.考虑全部震源对场址的影响，采用凹凸体模型设置地震动模拟方案，如图6所示，并根据公式(5)，计算每种方案总体的权重系数，如表2所示。
表2全部震源地震动模拟凹凸体模型方案权重系数
B.考虑全部震源对场址的影响，采用随机模型设置用于对比分析地震动模拟方案，如图7所示，并根据公式(6)计算每种方案的权重系数，如表3所示。
表3全部震源地震动模拟凹凸体模型方案权重系数
步骤(4)：采用随机有限断层法进行场址的地震动模拟，获得全部方案的加速度时程、加速度峰值和不同周期的加速度反应谱值
按照步骤(3)所设置的地震动模拟方案，采用随机有限断层法进行场址地震动模拟，获得全部方案的加速度时程、加速度峰值(PGA)以及不同周期的加速度反应谱值，周期点个数按实际需要采样。
步骤(5)：对考虑全部震源参与设计的每个方案的地震动模拟结果进行分析，并对场址地震动参数进行综合评价
分析全部震源参与设计的每个方案计算得到结果，将每个方案的凹凸体模型和随机模型计算得到的地震动加速度峰值(PGA)和不同周期(例如，0.03s-10s)的加速度谱值(Sa)，分别按照从小到大进行排列，并根据表2和表3列出其对应的加权系数，形成一系列场址地震动参数排列表，如表4和表5分别给出了M个方案凹凸体模型计算得到的PGA和相应于周期0.03s的谱值Sa03的排列表，表6和表7分别给出了M个方案随机模型计算得到的PGA和相应于周期0.03s谱值Sa03的排列表，其中pi＜pj＜L＜pm，Sa03_i＜Sa03_j＜L＜Sa03_m。
表4全部震源计算得到的场址地震动PGA排列表(凹凸体模型)
表5全部震源计算得到的场址地震动Sa03排列表(凹凸体模型)
表6全部震源计算得到的场址地震动PGA排列表(随机模型)
表7全部震源计算得到的场址地震动Sa03排列表(随机模型)
然后，根据上述场址地震动参数排列表，对考虑全部震源设计方案的计算结果进行统计分析：
A.基于全部震源设计方案计算得到的场址地震动参数排列表(如表4和表5)，对全部方案凹凸体模型计算得到的PGA值和不同周期点加速度反应谱值进行统计分析，给出最小值(Min)、50％分位数值(Q50)、均值(Mean)、84％分位数值(Q84)、95％分位数值(Q95)及最大值(Max)，如表8所示。
表8全部震源凹凸体模型的地震动参数统计值
B.基于全部震源设计方案计算得到的场址地震动参数排列表(如表6和表7)，对全部方案随机模型计算得到的PGA值和不同周期点加速度反应谱的谱值进行统计分析，给出最小值(Min)、50％分位数值(Q50)、均值(Mean)、84％分位数值(Q84)、95％分位数值(Q95)及最大值(Max)，如表9所示。
表9全部震源随机模型的地震动参数统计值
步骤(6)：对全部模拟方案的凹凸体模型和随机模型计算得到的参数统计值进行分析，考虑凹凸体模型和随机模型计算结果的较大值作为场址地震动参数；
步骤(7)：根据场址地震动参数，并考虑结构特性，对用于结构动力分析的加速度时程进行模拟。
基于步骤6得到的场址地震动参数，确定考虑全部发震构造最大可信地震影响的场址地震动参数(PGA，加速度反应谱等)。考虑全部发震构造最大可信地震影响的场址地震动参数可采用84％分位数的统计值。
根据结构特性，对用于结构动力分析的加速度时程按照3种方案进行模拟：
A：以场址地震活动环境为基础，挑选与场址加速度反应谱符合较好的实际地震加速度时程为种子时程，然后基于目标反应谱进行调整；
B：以场址地震活动环境为基础，建立地震动频率非平稳特性参数，人工拟合得到符合于目标反应谱的地震加速度时程；
C：从模拟得到的加速度时程中挑选出与目标反应谱符合较好的时程。
根据本申请的基于全部发震构造最大可信地震的地震动参数评价方法，其与考虑不同发震构造最大可信地震的地震动参数评价方法，在分析思路上都是采用多个方案来进行。区别点在于，不同发震构造是对每个震源单独进行方案设计，然后按照工程需要取其中一个震源或全部震源的结果的外包线作为评价结果；全部发震构造是对场址有影响的震源一起进行分析，可通过考虑未来一定时间内不同震源的发震概率来设置每个震源的加权系数，对全部方案的结果进行评价。两者在评价方法上存在显著的差异。
虽然本发明所揭露的实施方式如上，但所述的内容只是为了便于理解本发明而采用的实施方式，并非用以限定本发明。任何本发明所属技术领域内的技术人员，在不脱离本发明所揭露的精神和范围的前提下，可以在实施的形式上及细节上作任何的修改与变化，但本发明的专利保护范围，仍须以所附的权利要求书所界定的范围为准。
Claims (4)
1.一种基于全部发震构造最大可信地震的地震动参数评价方法，其特征在于，包括以下步骤：
(1)基于区域构造和区域地震地质资料，确定地震动模拟震源设置方案，并对其进行不确定性分析；
(2)确定地震动模拟方法，并对模拟结果影响较大的参数取值进行分析；
(3)考虑全部震源影响，对地震动模拟方案进行设计；
(4)采用随机有限断层法进行场址的地震动模拟，获得全部方案的加速度时程、加速度峰值和不同周期的加速度反应谱值；
(5)对考虑全部震源参与设计的每个方案的地震动模拟结果进行分析，并对场址地震动参数进行综合评价；
(6)对全部模拟方案的凹凸体模型和随机模型计算得到的参数统计值进行分析，取凹凸体模型和随机模型计算结果的较大值作为场址的参数；
(7)根据场址地震动参数，并考虑结构特性，对用于结构动力分析的加速度时程进行模拟。
2.根据权利要求1所述的地震动参数评价方法，其特征在于，
步骤(3)包括：
考虑全部震源对场址的影响，采用凹凸模型设置地震动模拟方案，并计算每种方案总体的权重系数；
采用随机模型设置用于对比分析的基于全部震源的地震动模拟方案，并计算每种方案的权重系数；
步骤(5)包括：
分析全部震源参与设计的每个方案计算得到结果，将每个方案的凹凸体模型和随机模型计算得到的地震动加速度峰值和不同周期的加速度谱值，分别按照从小到大进行排列，并根据其对应的加权系数，形成场址地震动参数排列表；
基于全部震源设计方案计算得到的场址地震动参数排列表，对全部方案凹凸体模型计算得到的地震动加速度峰值和不同周期点加速度反应谱值进行统计分析，给出最小值、50％分位数值、均值、84％分位数值、95％分位数值及最大值；
基于全部震源设计方案计算得到的场址地震动参数排列表，对全部方案随机模型计算得到的地震动加速度峰值和不同周期点加速度反应谱值进行统计分析，给出最小值、50％分位数值、均值、84％分位数值、95％分位数值及最大值。
3.根据权利要求1或2所述的地震动参数评价方法，其特征在于，步骤(1)中还包括：
(1)基于研究区域构造背景，确定对场址有重要影响的发震构造，并完成以下分析：
A.依据发震构造最新活动时代及地震震中分布特点，确定发震构造的震级上限；
B.依据地震、地质资料、野外调查、并结合区域深部构造背景，确定发震构造的位置，对不能确定明确位置的断裂，采用多个位置方案，并按确定位置的资料的可信度及对场址的影响水平分析，设置加权系数；
C.依据活动断裂的活动性及大震复发周期，对影响场址的不同发震断裂设置加权系数；
(2)依据研究区域地震、地质、小震精定位、震源机制解等资料，对发震断层的参数进行分析，确定地震动模拟震源方案，对于具有不确定性的参数，设置多个模拟方案，并设置各个方案的加权系数。
4.根据权利要求1或2或3所述的地震动参数评价方法，其特征在于，步骤(7)中的模拟方法为：
A.以场址地震活动环境为基础，挑选与场址加速度反应谱符合较好的实际地震加速度时程为种子时程，然后基于目标反应谱进行调整；
B.以场址地震活动环境为基础，建立地震动频率非平稳特性参数，人工拟合得到符合于目标反应谱的地震加速度时程；或者
C.从模拟得到的加速度时程中挑选出与目标反应谱符合较好的时程。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811528301.6A CN109375253B (zh) | 2018-12-13 | 2018-12-13 | 基于全部发震构造最大可信地震的地震动参数评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201811528301.6A CN109375253B (zh) | 2018-12-13 | 2018-12-13 | 基于全部发震构造最大可信地震的地震动参数评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN109375253A true CN109375253A (zh) | 2019-02-22 |
| CN109375253B CN109375253B (zh) | 2020-01-21 |
Family
ID=65374636
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201811528301.6A Active CN109375253B (zh) | 2018-12-13 | 2018-12-13 | 基于全部发震构造最大可信地震的地震动参数评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN109375253B (zh) |
Cited By (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110244355A (zh) * | 2019-07-25 | 2019-09-17 | 西南交通大学 | 一种基于震源断层模型的脉冲地震动模拟方法 |
| CN111159947A (zh) * | 2019-12-30 | 2020-05-15 | 中国水利水电科学研究院 | 高拱坝坝肩加固方法及装置 |
| CN111551994A (zh) * | 2020-05-11 | 2020-08-18 | 中国地震局地球物理研究所 | 一种匹配多目标的频率非平稳地震动拟合方法 |
| CN112882093A (zh) * | 2021-01-18 | 2021-06-01 | 中国测绘科学研究院 | 一种计算弹性地球内部同震变形的方法和系统 |
| CN113486507A (zh) * | 2021-06-28 | 2021-10-08 | 中国地震局工程力学研究所 | 一种地震动时程图确定方法、装置、电子设备及存储介质 |
| CN116127247A (zh) * | 2023-02-14 | 2023-05-16 | 中国地震局地球物理研究所 | 一种耦合多种震源模型的概率危险性解析计算方法 |
| CN116609832A (zh) * | 2023-07-18 | 2023-08-18 | 北京建筑大学 | 一种震源断层位错分布的确定方法和确定装置 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140022075A1 (en) * | 2011-02-26 | 2014-01-23 | Michael John Price | Customizable policy engine |
| CN107422376A (zh) * | 2017-07-03 | 2017-12-01 | 河海大学 | 一种非平稳非高斯地震动时程模拟方法 |
-
2018
- 2018-12-13 CN CN201811528301.6A patent/CN109375253B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20140022075A1 (en) * | 2011-02-26 | 2014-01-23 | Michael John Price | Customizable policy engine |
| CN107422376A (zh) * | 2017-07-03 | 2017-12-01 | 河海大学 | 一种非平稳非高斯地震动时程模拟方法 |
Non-Patent Citations (3)
| Title |
|---|
| 张翠然等: "基于汶川和芦山地震的沙牌坝址地震动输入研究", 《水力发电学报》 * |
| 张翠然等: "沙牌坝址的最大可信地震研究", 《水利学报》 * |
| 陈厚群等: "2008年汶川地震中沙牌拱坝的震情检验和分析", 《水利水电技术》 * |
Cited By (13)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN110244355B (zh) * | 2019-07-25 | 2021-06-08 | 西南交通大学 | 一种基于震源断层模型的脉冲地震动模拟方法 |
| CN110244355A (zh) * | 2019-07-25 | 2019-09-17 | 西南交通大学 | 一种基于震源断层模型的脉冲地震动模拟方法 |
| CN111159947B (zh) * | 2019-12-30 | 2023-04-07 | 中国水利水电科学研究院 | 高拱坝坝肩加固方法及装置 |
| CN111159947A (zh) * | 2019-12-30 | 2020-05-15 | 中国水利水电科学研究院 | 高拱坝坝肩加固方法及装置 |
| CN111551994A (zh) * | 2020-05-11 | 2020-08-18 | 中国地震局地球物理研究所 | 一种匹配多目标的频率非平稳地震动拟合方法 |
| CN112882093A (zh) * | 2021-01-18 | 2021-06-01 | 中国测绘科学研究院 | 一种计算弹性地球内部同震变形的方法和系统 |
| CN112882093B (zh) * | 2021-01-18 | 2024-03-05 | 中国测绘科学研究院 | 一种计算弹性地球内部同震变形的方法和系统 |
| CN113486507B (zh) * | 2021-06-28 | 2022-09-13 | 中国地震局工程力学研究所 | 一种地震动时程图确定方法、装置、电子设备及存储介质 |
| CN113486507A (zh) * | 2021-06-28 | 2021-10-08 | 中国地震局工程力学研究所 | 一种地震动时程图确定方法、装置、电子设备及存储介质 |
| CN116127247A (zh) * | 2023-02-14 | 2023-05-16 | 中国地震局地球物理研究所 | 一种耦合多种震源模型的概率危险性解析计算方法 |
| CN116127247B (zh) * | 2023-02-14 | 2023-08-18 | 中国地震局地球物理研究所 | 一种耦合多种震源模型的概率危险性解析计算方法 |
| CN116609832A (zh) * | 2023-07-18 | 2023-08-18 | 北京建筑大学 | 一种震源断层位错分布的确定方法和确定装置 |
| CN116609832B (zh) * | 2023-07-18 | 2023-10-27 | 北京建筑大学 | 一种震源断层位错分布的确定方法和确定装置 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN109375253B (zh) | 2020-01-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109375253A (zh) | 基于全部发震构造最大可信地震的地震动参数评价方法 | |
| CN109375252A (zh) | 考虑不同发震构造最大可信地震的地震动参数评价方法 | |
| Chaljub et al. | Quantitative comparison of four numerical predictions of 3D ground motion in the Grenoble Valley, France | |
| CN106814391A (zh) | 基于菲涅尔体层析反演的地面微地震事件定位方法 | |
| CN105093319B (zh) | 基于三维地震数据的地面微地震静校正方法 | |
| CN108254780A (zh) | 一种微地震定位及各向异性速度结构层析成像方法 | |
| CN107894618B (zh) | 一种基于模型平滑算法的全波形反演梯度预处理方法 | |
| CN105116444B (zh) | 一种地面微地震监测各向异性速度模型 | |
| CN105386756A (zh) | 一种应用应变量计算脆性地层孔隙度的方法 | |
| CN103149586A (zh) | 一种倾斜层状粘弹性介质中波场正演模拟方法 | |
| CN105301639A (zh) | 基于vsp旅行时双加权层析反演速度场的方法及其装置 | |
| CN106646645A (zh) | 一种新的重力正演加速方法 | |
| CN104360396B (zh) | 一种海上井间tti介质三种初至波走时层析成像方法 | |
| CN103913768A (zh) | 基于地震波资料对地表中浅层进行建模的方法及装置 | |
| CN104166159A (zh) | 四维微地震监测的裂缝形态处理方法和系统 | |
| CN105717538B (zh) | 起伏地表地震数据偏移基准面转换方法及装置 | |
| Bydlon et al. | Combining dynamic rupture simulations with ground‐motion data to characterize seismic hazard from Mw 3 to 5.8 earthquakes in Oklahoma and Kansas | |
| US20240134080A1 (en) | Method and System for Real-Time Calculating a Microseismic Focal Mechanism Based on Deep Learning | |
| Tao et al. | Rupture directivity and hanging wall effect in near field strong ground motion simulation | |
| Magrin | Multi-scale seismic hazard scenarios | |
| Cao et al. | Slip‐correlated corner frequency for stochastic finite‐fault modeling of ground motion | |
| CN114089419B (zh) | 优化的变网格地震正演模拟方法 | |
| Edwards et al. | Ground motion prediction equations | |
| Ghodrati Amiri et al. | Ground motion prediction equations (GMPEs) for elastic response spectra in the Iranian plateau using gene expression programming (GEP) | |
| CN103837894A (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 |