地震安全性评价课程第10章

SSEISMICSAFETYCHAPTER10地震安全性评价课程地震危险性的概率分析概率地震危险性分析（PSHA）的基本原理、参数模型与工程应用CHAPTER10适用专业土木工程·地质工程CONTENTS本章内容概览01引言与基本概念第3–7页：从确定性方法的局限出发，阐明概率分析的核心思想04地震动预测方程第19–22页：GMPE的表达形式、选择论证02地震区带与潜在震源区第8–12页：讲解三级层级源区模型的划分原则05概率计算与不确定性第23–27页：全概率积分公式、逻辑树方法03地震活动性参数确定第13–18页：涵盖G-R关系拟合、完备性校正06结果表述与工程应用第28–31页：危险性曲线解读、统一反应谱生成地震安全性评价·概率地震危险性分析CHAPTER01引言与基本概念从确定性到概率分析的范式转变PROBABILISTICANALYSISSEISMICSAFETYEVALUATION为什么需要概率分析确定性方法仅考虑"最大潜在地震"，无法回答"该地震在工程寿命期内发生的可能性有多大"；概率分析将地震发生的时间、空间、强度均视为随机变量，给出不同超越概率下的地震动参数，为基于风险的抗震设计提供定量依据。确定性方法的局限确定性方法仅考虑最大潜在地震置于最近位置，给出的地震动值是"可能的最大值"而非"可能发生的值"，无法回答工程寿命期内遭遇该水平的概率有多大。规范设防要求现代抗震设计规范（如GB50011）已全面采用基于概率的设防目标，要求提供50年超越概率63%、10%、2%三水准地震动参数，只有PSHA能系统给出这些值。全概率整合概率分析通过全概率公式整合区域内所有震源、所有震级和距离的贡献，避免了确定性方法对单一控制震源的过度依赖，结果更具统计稳健性。风险决策支持PSHA输出的危险性曲线可直接用于风险决策：业主可根据可接受的风险水平选择对应的设防水准，实现安全与经济的最优平衡。THEORY概率地震危险性分析的核心思想PSHA将地震危险性表达为"场地处地震动参数超过某一给定值的年超越概率"，通过全概率公式对区域内所有潜在震源、所有可能震级和距离进行积分叠加，最终得到一条完整的危险性曲线。概率定义的完整性：PSHA将地震危险性定义为"场地处地震动参数A超过给定值a的年超越概率P(A≥a)"，这是一个完整的概率分布而非单一数值全概率公式的应用：对所有潜在震源i、所有可能震级M和距离R，将ν(M,R)·P(A≥a|M,R)积分叠加得到总超越率λ(a)关键参数的导出：ν(M,R)是震级-距离联合发生率，由G-R关系和空间分布函数导出；P(A≥a|M,R)由地震动预测方程的对数正态分布给出危险性曲线的意义：最终得到的λ(a)-a曲线即危险性曲线，横轴为地震动值，纵轴为年超越率或其倒数重现期，是工程设计取值的直接依据MethodologyPSHA的基本流程与技术框架完整PSHA包含四个核心环节：潜在震源区划分与建模、地震活动性参数确定、地震动预测方程选取、全概率积分计算与不确定性校正；每个环节的输入质量直接决定最终结果的可靠性。01潜在震源区划分与几何建模划分地震区、地震带和潜在震源区，建立几何模型；这一步决定了后续参数统计的空间单元和积分范围02地震活动性参数确定确定每个统计区的G-R参数(a,b)、空间分布函数f(x,y)和各潜在震源区的震级上限Mu，构成地震活动性模型03地震动预测方程选取选择适用于场址区域的地震动预测方程(GMPE)，并论证其在震级、距离和构造环境方面的适用性04全概率积分计算与不确定性校正执行全概率积分计算，得到均值危险性曲线；通过逻辑树处理认知不确定性，输出多分位结果供工程选用S地震安全性评价TERMINOLOGY关键术语辨析超越概率指在指定时段内地震动参数≥某阈值的概率；重现期是该超越概率的倒数；年超越率是单位时间内事件发生的平均次数；三者互为换算关系，是理解PSHA输出结果的基础概念。01超越概率P(T,t)在t年内地震动参数≥某阈值的概率，泊松假设下P=1-exp(-λt)，其中λ为年超越率。02重现期Tr超越概率的倒数概念，Tr=1/λ表示该水平地震动平均多少年出现一次；注意重现期不等于"每隔Tr年必然发生一次"。03年超越率λ单位时间内地震动参数超过给定值的期望次数，是PSHA的直接输出量；多条危险性曲线即λ随阈值变化的函数。04三者换算示例50年超越概率10%对应λ=-ln(0.9)/50≈0.0021/年，重现期Tr≈475年；这是我国抗震规范"中震"设防基准的概率含义。CHAPTER02地震区、带与潜在震源区划分概率分析的几何建模基础SEISMICZONATION地震区与地震带的划分原则地震区依据地震活动空间分布的分区性和构造背景的相似性划定；地震带在地震区内按成带性和活动构造一致性进一步细分；两者共同构成地震活动性参数统计的基本单元。地震区——最高级分区单元：依据地震活动空间分布的分区性和深部构造背景的相似性划定；同一地震区内地震成因机制基本一致。地震带——区内细分单元：在地震区内按地震活动的成带性和活动构造带的一致性进一步细分；同一条带上的地震往往受同一组活动断裂控制。多源资料综合判定：划分时需综合考虑历史地震分布、活动构造图、地球物理场特征和地壳结构资料；边界通常沿构造单元界线或地震活动梯度带划定。参数统计的基本单元：地震区和地震带不仅是描述性分区，更是G-R参数统计和空间分布函数构建的基本单元，划分不合理会导致参数估计偏差。SEISMICSOURCEMODELPROBABILISTICANALYSIS潜在震源区的三级层级结构GB17741-2025规定概率分析采用三层级源区模型：地震统计区（提供G-R关系与空间分布函数）、背景潜在震源区（反映弥散地震活动）、构造潜在震源区（对应具体发震构造），三者嵌套互补。01地震统计区

最外层单元提供区域G-R关系和空间分布函数；其范围通常覆盖一个完整的地震带或构造区，面积数千至数万平方公里。02背景潜在震源区

弥散地震活动位于统计区内、构造源之外，用于容纳未被明确发震构造解释的弥散地震活动；其震级上限一般低于构造源0.5-1.0级。03构造潜在震源区

发震构造对应具体的发震构造（断层或褶皱），边界依据断层几何、历史地震分布和古地震资料划定；每个构造源有独立的震级上限。三层嵌套互补

建模校验三层嵌套确保模型既反映已知构造的控制作用，又不遗漏背景地震活动；实际建模时需检查所有观测地震是否都被某一源区覆盖。METHODOLOGY构造潜在震源区的边界确定方法构造潜在震源区以发震构造的几何产状、历史地震空间分布和古地震资料为依据划定边界；需考虑断层倾角、破裂宽度及定位不确定性，通常沿断层走向外扩一定缓冲带作为概率扩散范围。01走向方向沿断层走向方向，边界通常延伸至历史地震序列端点或断层几何不连续处（弯折、阶跃）；若无明确端点，可按断层全长外推10%-20%02垂直走向方向需考虑断层倾角和破裂宽度；高倾角断层缓冲带较窄（3-5km），低倾角逆冲断层可达10-20km以涵盖上盘效应03深度方向破裂下限通常取脆性-韧性转换带深度（大陆区约15-20km）或历史地震最大震源深度；浅源与深源地震可能需分层建模04不确定性处理边界划定中的不确定性应通过逻辑树处理：设置"窄边界""基准边界""宽边界"三个分支，分别赋予权重，体现几何建模的认知不确定性SEISMICSOURCEMODELBACKGROUNDSOURCE背景潜在震源区的作用与划定背景源用于容纳未被明确发震构造解释的中小地震活动，其震级上限一般低于构造源；划定范围覆盖地震统计区内构造源之外的区域，确保所有观测到的地震活动在模型中都有归属。01核心作用——容纳弥散地震背景源的核心作用是容纳"弥散地震"——那些无法归因于已知发震构造、但在统计上显著存在的中小地震活动。02划定范围划定范围为地震统计区扣除所有构造源后的剩余区域；若构造源密集重叠，背景源可能被分割为多个不连通子区，各自独立统计。03震级上限Mu,bg通常取该区域内历史最大弥散地震震级加0.5级作为裕度；若无历史记录，可参照邻近类似构造区的经验值。04忽略背景源的风险忽略背景源会导致低估中小地震的危险性贡献，尤其在构造源稀疏但背景活动频繁的地区（如华北平原内部），背景源贡献可达总量的30%以上。CHAPTER03地震活动性参数的确定G-R关系、空间分布函数与时间模型SEISMICITYMODEL震级-频度关系（G-R关系）lgN=a-bM描述震级M以上地震的年累积次数N；b值反映大小地震比例，通常在0.6-1.2之间；a值反映整体活动水平；两者由历史地震目录经完备性校正后回归得到，是PSHA最核心的输入参数。01公式释义—lgN=a-bM是地震学最基本的经验关系之一：N为震级≥M的地震年累积次数，a表征整体活动水平，b表征大小地震比例02b值的物理意义—b值通常在0.6-1.2之间；b偏低意味着大地震占比偏高（如俯冲带megathrust），b偏高则小震占主导（如火山或地热区）03参数回归方法—a值与b值须由完备地震目录经最小二乘或最大似然法回归得到；使用不完备数据会系统性低估b值、高估大地震发生率u04截断震级Mu—G-R关系仅在Mmin≤M≤Mu范围内有效；超出Mu后发生率截断为零，这一截断对长周期危险性影响极大，需审慎确定MuSEISMICITYCATALOGCOMPLETENESS地震目录完备性校正历史地震记录随时间推移逐渐完整，早期小震遗漏严重；需分时段确定最小完备震级Mc，仅使用Mc以上的数据进行G-R拟合；常用方法包括MAXC法、MBASS法和EMR法，选择不当会显著偏估b值。问题根源历史地震目录在早期仅记录破坏性大震，小震严重遗漏；直接使用全时段数据拟合G-R会低估b值、高估大震风险核心方法完备性校正的核心是分时段确定最小完备震级Mc(t)：例如1900年前Mc=6.0，1900–1960年Mc=5.0，1960年后Mc=3.0常用估计方法MAXC（最大曲率法）、MBASS（分段拟合转折点）和EMR（整个震级范围拟合法）；建议多种方法交叉验证拟合与评估校正后的G-R拟合仅使用各时段Mc以上的数据；拟合优度可用χ²检验或AIC准则评估，残差过大时应重新审视分区合理性SSeismicSafetyMETHODOLOGYSpatialDistribution空间分布函数的构建空间分布函数f(x,y)描述地震统计区内不同位置发生地震的相对概率密度；可基于历史地震核密度估计、活动构造权重或平滑模型构建；归一化后在整个统计区上积分为1，是将区域发生率分配到各潜在震源区的关键桥梁。核心定义：空间分布函数f(x,y)将地震统计区的总发生率按空间位置分配：∫∫f(x,y)dxdy=1，f(x,y)≥0；它是连接区域G-R关系与各潜在震源区的关键桥梁。核密度估计法（KDE）：最常用方法，以每个历史地震为中心放置二维高斯核，带宽h控制平滑程度。h过小则噪声大，过大则抹平构造信号，需依据区域地震目录密度合理选取。活动构造权重：在已知断层附近人为提高f值，远离断层处降低；权重系数需通过敏感性分析校准，避免主观过强导致空间分布失真。稀少地区处理：对于地震稀少的地区，纯数据驱动的KDE不稳定，可采用平滑模型或借用邻近相似构造区的空间分布模式作为先验，以增强估计稳健性。SEISMICITYMODELPSHA·TEMPORAL地震发生的时间模型标准PSHA假设地震发生服从泊松过程，即事件独立、平稳、无记忆；该假设使年超越概率P=1-exp(-λt)成立；对于具有准周期特征的大震，可采用时间相依模型修正，但需充分的古地震证据支持。01泊松过程假设标准PSHA假设地震服从平稳泊松过程：事件独立、发生率恒定、无记忆性；这使得年超越概率P=1−exp(−λt)严格成立。02假设适用范围泊松假设对中小地震通常合理，但对具有准周期复发特征的大震（如板块边界特征地震）可能低估短期风险或高估长期风险。03时间相依模型时间相依模型（如BPT、对数正态复发模型）引入距上次大震的时间间隔作为条件变量，可更准确描述周期性复发行为。04数据与保守性采用时间相依模型需充分的古地震或历史重复地震证据支持；若数据不足强行使用，反而引入更大不确定性，此时泊松模型仍是保守稳妥的选择。地震安全性评价17SEISMICPARAMETERS震级上限Mu的确定方法震级上限是潜在震源区可能发生的地震最大震级；构造源依据断层长度-震级经验公式、古地震和历史最大震综合判定；背景源取统计区内弥散地震的最大观测震级加适当裕度；Mu对长周期危险性影响显著。01构造源Mu估算首选断层长度-震级经验公式（如Wells&Coppersmith1994）估算，再与历史最大震和古地震特征震级对比取高值。02断层分段与联合破裂若断层分段明确，需判断Mu对应单段破裂还是多段联合破裂；联合破裂情景的Mu更高但频率更低，应作为逻辑树分支处理。03背景源Mu确定取统计区内弥散地震的最大观测震级加0.5级裕度；若无观测记录，参照全球类似构造背景的经验上限（通常5.5–6.0级）。04长周期敏感性与审慎原则Mu对长周期（T≥2s）危险性影响远大于短周期；敏感性分析显示Mu变化0.5级可使Sa(3s)变化30%–50%，因此需特别审慎。CHAPTER04地震动预测方程连接震源与场地的桥梁SeismicGroundMotionGMPE地震动预测方程的表达形式GMPE以数学函数或表格形式给出地震动参数随震级、距离变化的期望值与标准差；需反映高频饱和特性，必要时包含断层类型项；反应谱周期点在对数坐标下均匀分布且不少于20个，最长周期不小于6秒。标准形式—GMPE的标准形式为lgY=c1+c2·M+c3·lgR+c4·R+εσ，其中Y为地震动参数，M为震级，R为距离，ε为标准正态随机变量扩展参数项—现代GMPE还包含断层类型项（走滑/逆冲/正断）、场地条件项（Vs30分类）和盆地深度项等，以更精细地捕捉物理效应反应谱周期—反应谱GMPE需在0.01–6s（或更长）范围内给出20个以上周期点的系数；周期点在对数坐标下近似均匀分布，确保插值精度高频饱和—高频段（T<0.1s）应体现震级饱和和距离饱和特性，避免大震近场预测值无限增长；这通常通过非线性M-R耦合项实现METHODOLOGYGMPE的选择与适用性论证强震动数据充足的地区应采用本地统计建立的GMPE；数据不足时可用类比法引入邻近地区方程并做校正；选用时需验证其对场址所在构造环境、震级范围和距离范围的适用性，不得外推使用。Priority优先选用本地回归方程优先选用基于本地强震动数据回归的GMPE；中国西部地区可选用俞言祥等(2021)或陈鲲等(2022)的最新方程Fallback全球数据库类比与校正本地数据不足时，可从NGA-West2、ESM等全球数据库中选择构造环境相似的方程，并用少量本地记录做残差校正（host-to-target调整）Validation适用性三维论证适用性论证需覆盖三方面：震级范围（方程适用范围是否包含Mu）、距离范围（是否涵盖近场至远场）、构造类型（走滑/逆冲是否与场址一致）Constraint禁止外推与不确定性说明不得将方程外推到标定范围之外；若场址条件超出所有可用方程的适用范围，应在报告中明确说明并给出不确定性放大因子DISTANCEMETRICS距离度量与近场处理GMPE中的距离可选震源距、断层距或投影距，必须与方程定义一致；近场（≤10km）强震动数据稀缺，需基于有限观测单独评定或采用近场饱和模型；距离计算的准确性直接影响短周期段的危险性估计。01距离定义的严格匹配GMPE中的距离定义必须严格匹配：Repi（震中距）、Rhyp（震源距）、Rrup（断层距）、Rjb（Joyner-Boore投影距）不可混用。02断层距Rrup的优势与局限断层距在现代GMPE中最常用，因其最能反映有限断层辐射场的几何衰减；但计算需已知断层面几何参数，对小震或弥散地震常退化为点源距离。03近场数据稀缺性近场（R≤10km）强震动记录稀缺且变异性大；GMPE在此区间的外推可靠性低，需基于有限的近场观测单独评定或采用近场饱和模型约束。04GB17741-2025近场评定要求当最近距离≤10km时，GB17741-2025要求基于近场强震动数据统计评定地震动参数，而非单纯依赖GMPE外推，以避免低估或高估近场危险性。CHAPTER05概率危险性计算与不确定性校正全概率积分与认知不确定性的处理05PROBABILISTICSEISMICHAZARD全概率积分的计算公式年超越率λ(a)=Σ∫∫ν(M,R)·P(A≥a|M,R)·f(M)f(R)dMdR，对所有震源、震级和距离积分；实际计算离散化为多重求和，需保证震级步长≤0.1、距离步长足够密以收敛。λ(a)=Σi∫∫νi(M,R)·P(A≥a|M,R)·f(M)·f(R)dMdR总超越率公式：λ(a)=Σi∫∫νi(M,R)·P(A≥a|M,R)dRdM，对所有n个震源求和联合发生率：νi(M,R)=Ni(M)·fi(R|M)，其中N(M)由G-R关系给出，f(R|M)由空间分布函数和源区几何导出超越概率：P(A≥a|M,R)=1−Φ[(lga−lgYmean(M,R))/σ]，Φ为标准正态CDF，lgYmean由GMPE给出，σ为对数标准差物理含义：将所有可能的(M,R)组合下"发生该级别地震的概率"乘以"该地震产生超限地面运动的概率"，再累加得到总超越率S地震安全性评价NUMERICALINTEGRATION数值积分的实现要点震级离散化从Mmin到Mu，步长0.05-0.1；距离按对数或线性分段，近场加密；每个(M,R)组合调用GMPE获取均值和对数标准差σ，再查正态分布表得超越概率；总计算量与震源数×网格数成正比。01震级离散化步长—ΔM取0.05–0.1；步长过大会遗漏G-R曲线的非线性段，尤其在高震级端截断附近；建议做步长收敛测试以确认结果稳定02距离离散化策略—近场加密（ΔR=1–2km）、远场放宽（ΔR=10–20km）或对数等间距；总距离点数通常50–200个以保证积分精度03GMPE调用与概率计算—每个(M,R)网格点调用GMPE计算均值和对数标准差，再查正态分布表或用erfc函数计算超越概率；向量化运算可大幅提速04结果验证—计算完成后需验证：①λ(a)单调递减；②a→0时λ(a)趋近总发生率；③a极大时λ(a)→0；若不满足则检查积分边界或GMPE调用错误S地震安全性评价METHODOLOGYEpistemicUncertainty认知不确定性与逻辑树方法GMPE选择、b值、Mu等均存在认知不确定性；逻辑树将每个不确定因素设为分支节点，赋予权重，遍历所有路径计算加权平均危险性曲线；分支过多时需采样近似，避免组合爆炸。认知不确定性的本质源于知识不完整（如GMPE选哪个、b值取多少、Mu定多大），不同于自然变异性的偶然不确定性，原则上可通过更多研究减小。逻辑树的基本原理将每个不确定因素设为一个节点，每个节点分出若干互斥分支并赋予权重（总和为1）；遍历所有路径得到加权平均危险性曲线。典型逻辑树结构包含3–5个节点：GMPE选择（3–4条）、b值（±0.05~0.1）、Mu（±0.3~0.5级）、空间分布模型（KDEvs构造加权）等。分支爆炸与采样近似分支总数超过数百条时全枚举计算量过大；可采用蒙特卡洛采样或拉丁超立方抽样近似，采样数通常1000–5000次即可收敛到均值±5%。UNCERTAINTYPSHA·LECTURE27偶然不确定性与分位数表达GMPE的对数标准差σ代表同一(M,R)下地震动的自然变异性，属偶然不确定性；PSHA结果通常给出均值危险性曲线，也可给出16%、84%分位曲线以表征置信区间；工程应用常取均值或保守分位。对数标准差σ量化偶然不确定性偶然不确定性由GMPE的对数标准差σ量化，代表相同(M,R)下不同地震事件的天然变异性；σ通常为0.5–0.8（以10为底对数单位）。均值曲线与置信范围的缺口PSHA均值曲线已内含σ的影响（通过对数正态积分），但无法反映结果本身的置信范围；为此需输出分位数曲线。16%与84%分位曲线的工程含义16%分位曲线表示"有84%把握实际危险性高于此值"，84%分位反之；工程保守设计常取84%分位或均值+1σ作为设防基准。分位数的两种计算方法分位数可通过逻辑树各路径结果的排序统计得到，也可在单次计算中对σ做偏移模拟；两种方法等价但后者计算效率更高。CHAPTER06结果表述与工程应用从危险性曲线到设计地震动参数SEISMICHAZARDHAZARDCURVE危险性曲线的解读与应用横轴为地震动参数（PGA或Sa(T)），纵轴为年超越率或其倒数重现期；工程设计通常取50年超越概率10%（重现期475年）和2%（2475年）对应的地震动值作为设防基准；曲线形态反映场地受远近震控制的程度。01坐标系与曲线走向横轴为PGA或Sa(T)，纵轴为年超越率λ（对数坐标）或重现期Tr=1/λ；曲线从左上方（高频小震）向右下方（低频大震）单调递减。小震高频→大震低频02设防基准取值50年超越概率10%对应λ=0.0021/年（Tr=475年），为我国"中震"设防基准超越概率2%对应λ=0.0004/年（Tr=2475年），为"大震"基准475年中震重现期2475年大震重现期03曲线斜率与构造源陡降：说明小震主导、大震罕见，危险性随震级快速衰减平缓：大震贡献显著，通常意味着近场存在高Mu构造源陡降型平缓型04工程取值注意事项同一场地不同周期的危险性曲线形态可能差异很大——短周期受近场小震控制，长周期受远场大震控制，因此不可对不同周期的谱值做简单缩放，须逐周期独立取值。禁止跨周期简单缩放设防谱值SEISMICANALYSISPSHA·UHRS统一危险性反应谱的生成对多个周期分别做PSHA得到各自的重现期-谱值关系，在同一超越概率下截取各周期谱值，连成场地相关反应谱；该谱反映了场址特有的地震环境，不同于规范标准化谱形，可直接用于结构抗震验算。01逐周期执行PSHA对每个感兴趣的结构周期Ti分别执行PSHA，得到λ–Sa(Ti)曲线；在同一目标超