地震损失评估建模实战从哥伦比亚7.4级强震、188人遇难看烈度衰减与损失预测当地时间8月10日上午7点34分哥伦比亚西部发生7.4级强震震中位于首都波哥大以西约240公里是哥伦比亚过去十年来录得的最强地震。据中新社8月12日报道哥伦比亚全国首府协会8月11日发布的最新统计显示此次强震已造成188人死亡、1677人受伤考卡山谷省首府卡利市等多地建筑受损巴拿马、厄瓜多尔、委内瑞拉等邻国均有震感联合国秘书长古特雷斯也通过发言人发声表示哀悼。每一次强震之后都会有一批人试图回答同一个问题地震造成的损失到底有多大要预测损失先要预测破坏要预测破坏先要理解震动本身是如何从震中向外衰减的。这一整条链就是地震损失评估建模的主线。一、两个容易混的概念震级与烈度建模之前先把两个词分清。震级衡量的是地震释放能量的大小是一个物理量全球只有一个数值烈度衡量的是地震在某地造成破坏的强弱和距离、场地条件有关同一个地震在不同地方烈度完全不同。概念衡量什么一个地震有几个值影响因素震级震源释放的能量1个震源本身烈度地表的破坏程度多个每个地点一个距离、震源深度、场地土质、建筑质量这次哥伦比亚地震中国地震台网测定的震级为7.5级、震源深度80公里哥伦比亚地质局测定为7.4级——不同机构测出微小差异是正常的因为台网分布和算法不同。而烈度则是卡利市区和波哥大市区各有一个数值前者离震中更近烈度更高。建模的第一个任务就是把震中附近震级高这种定性描述变成离震中X公里处烈度大约是多少的定量关系。二、烈度衰减模型地震的能量账本烈度衰减的核心规律可以用一句话概括离震中越远震动越弱但衰减不是线性的。经典的做法是用经验公式拟合最常见的形式是震级项减去距离项以国内常用的烈度衰减模型为例某地烈度大致等于一个与震级有关的常数项加上震级与对数距离的线性组合。用中文说就是烈度 基础值随震级增大而增大 距离修正离震中越远衰减越多且用对数刻度描述衰减速度再考虑震源深度和场地条件的修正项。为什么要用对数距离而不是直接用公里数因为地震动的衰减是先快后慢震中附近几百米烈度掉得飞快到了上百公里外每多100公里衰减得就不那么明显了。对数函数恰好能刻画这种早期快速衰减的特征这也是地震学教材里最常见的模型形式。建模时通常分三步。第一步收集历史地震的烈度观测点数据震级、距离、观测烈度第二步用最小二乘回归估计模型参数第三步代入目标地震的震级和位置输出各区域的预测烈度。拟合优度、残差分布是必须检查的环节——如果残差随距离呈现系统性偏差说明模型形式选错了。三、从烈度到损失一条完整的评估链有了区域烈度下一步是把烈度翻译成损失。损失评估的标准链条是烈度→破坏等级→经济损失每一环都是一个子模型。第一环烈度到破坏等级。按我国现行规范建筑破坏分为基本完好、轻微破坏、中等破坏、严重破坏、倒塌五档。对同一类建筑不同烈度下处于各破坏等级的比例不一样——例如某烈度下砖混结构中等破坏及以上的占比可能达三成而钢结构可能只有一成。这组比例通常来自震害调查统计是评估的核心参数。第二环破坏等级到经济损失。每栋建筑的修复或重建成本按破坏等级折算轻微破坏取造价的很小比例倒塌取全额。把区域内每栋建筑的造价×折损比例加起来就是直接经济损失。第三环人员伤亡估算。常用方法是建筑倒塌率→被困人数→伤亡率的链式估算结合人口密度与发生时间白天还是凌晨人员是否在室内输出伤亡的期望区间。评估环节输入输出常用方法烈度预测震级、位置、深度区域烈度分布烈度衰减经验公式回归破坏评估烈度、建筑类型各破坏等级比例易损性矩阵查表损失计算破坏比例、造价经济损失加权求和伤亡估算倒塌率、人口密度伤亡区间链式估算情景分析四、这次哥伦比亚地震模型能告诉我们什么把哥伦比亚这次地震代入评估框架可以检验每个环节的价值与局限。烈度环节7.4级、震源深度80公里、震中距卡利市较近——按衰减模型推算震中附近烈度应达到较高的破坏烈度等级而距离较远的波哥大震动明显减弱这也与邻国有震感、波哥大未现大规模伤亡报道的公开信息一致。损失环节伤亡高度集中在震中附近的考卡山谷省说明距离衰减在现实中真实上演——这就是衰减模型最直接的应用在地震发生后数小时内不需要逐村实地调查先用模型圈出高烈度区指导救援力量优先投放。188人遇难、1677人受伤的官方统计则是检验模型预测的标答。但模型也有明显的局限。第一伤亡数字随时间更新最初统计111人随后升至188人外媒后续报道已超200人模型给的是期望值现实充满离散的偶然性第二建筑质量、土质条件这些场地效应难以精确参数化同一烈度下不同城市的损失可以差出好几倍第三数据匮乏时经验公式的外推风险很大不能用别国的模型参数直接套用。五、国赛怎么用这类题建模启示与延伸地震损失评估是国赛非常喜欢出题的灾害类背景它的建模启示可以归纳为三点。第一链条式建模。一个复杂问题拆成烈度→破坏→损失→伤亡四个环节每个环节单独建、再串联比一个大而全的模型清晰得多也更容易拿分。第二参数获取的合理性。灾害建模中参数从哪来是评委必问的问题合理做法是优先引用公开的震害统计、规范表格实在没有就明确说明假设并做敏感性分析展示参数变化对结论的影响。第三模型的落地表达。最终成果要能指导行动输出一张高烈度区分布图、一张建议救援优先级表比一个抽象的预测数值有用得多。评分标准里模型的应用价值往往比模型的理论深度更值钱。灾害面前模型不能阻止地震但可以让救援更快、损失更小一点——这就是这类建模题最朴素的意义。六、给参赛者的三个实战提醒把这套框架落到赛场上还有三个实战提醒都是容易丢分的地方。提醒一烈度表要引用权威依据。我国的烈度评估通常参照国家标准里对烈度与破坏现象的描述表如房屋破坏程度、人的感觉、器物反应等建模时务必注明烈度档位的划分依据来自何处哪怕引用的是教材或规范条文——评委看到依据明确四个字比看到任何花哨算法都放心。提醒二数据来源分三层交代。灾害类赛题的数据往往有三种来源公开的历史震害统计、题目附件、自己合理假设的参数。写论文时把这三层分开标注参数假设部分主动做敏感性分析——哪个参数变化对结论影响最大就重点讨论哪个这比列一堆公式更能体现建模能力。提醒三把结论翻译成行动语言。同样是预测烈度为八度好的表述是建议将高烈度区的救援优先级设为第一档并在24小时内完成区域人员排查。评委评分表里成果的应用价值一栏吃的就是这种把数字翻译成决策的功夫。建得出模型是及格用得上模型才是优秀。记住这三条地震类赛题就不只是会做而是做得漂亮。