强震巨灾后死亡人数早期预测方法.docVIP

强震巨灾后死亡人数早期预测方法【摘要】目的 利用强震巨灾后人员死亡时增量 K 值随时间 t 变化的函数关系，研究早期评估死亡总人数的方法。方法 选取中国国际救援队提供的持续性信息保障工作以来积累的强震巨灾的10次典型震例，利用统计学软件SPSS 17.0对数据进行计算，并建立数学模型。结果 人员死亡时增量 K 值最大值（ Kmax）所对应的时间点为 T，为8～12 h，2 T 时间内的曲线呈近似正态分布，对正态分布曲线进行积分求和估计线下面积即2T时间内的上报死亡人数W，即对W=∫02Tf（t）dt函数计算。W与死亡总数 M 正相关（ P 0.05），时增量 K 值在0≤ t ≤2 T 时间内变化的函数关系为f（t）=Kmax e — π K max2（t—T）2，具体结果见表1。 2.2 10个震例在据震后0～2 T 时间内估计的死亡人数（ W ） 2 T 时间内 K 值曲线下的面积为死亡统计上报的人数（ W ），计算2 T 时间内的死亡上报人数，即为计算函数f（x）的线下面积，正态分布的线下面积利用微积分公式可得：W=∫02Tf（t）dt。依据此函数，本研究计算出10个强震巨灾震例2 T 时间内的死亡上报总人数（ W ），结果见表1。 从表1可知，10个震例中， T 时即据发震时间最短4.8 h，最长30.8 h，（12.94±8.18） h。另外，死亡总人数越多，预测的2 T 时间内死亡人数也越多（震例2、3、4、7）。但预测的2 T 时间内死亡人数比最后统计的实际死亡人数相差较大（震例7相差4倍多）。 2.3 预测值评价 死亡总数（ M ）是应急救援工作结束后权威部门的统计数，与2 T 时间内的死亡上报数值呈一定的相关关系。 T 时上报死亡人数是实际统计的数值，2 T 时间内预测的死亡人数（估计的死亡人数 W ）， W 是根据 T 时和 Kmax值对死亡人数的函数式推导的数值，见表2。 从表2可知，预测的2 T 时间内死亡人数和 T 时实际上报死亡人数与最后权威部门统计的死亡总数呈高度正相关，并且均有统计学意义。 2.4 模型模拟 由 T 时上报的死亡人数预测2 T 时的死亡上报人数，由于2 T 上报人数与死亡总人数存在相关关系，本研究对10个强震巨灾震例，在震后2 T 时间内预测的死亡人数与最后权威部门统计的死亡总人数做模拟回归方程，共建立9个模拟模型，模型模拟结果和参数估计结果见表3。 表3的9个模拟模型中，模型模拟的 R 2值拟合效果优势：幂函数指数函数或复合函数三次方函数二次方函数或线性函数S曲线函数对数函数逆模型函数。因此，选择幂函数为模型模拟函数， R 2值为0.88，由该函数的参数可得死亡总人数 M 的模拟函数式， M = W 1.23×0.194， W 为预测的2 T 时间内死亡人数。利用该函数式可在 T 时间点估算某震例的死亡总人数，从而达到早期判断死亡总人数的目的。 3 讨论 地震巨灾发生后，人员死亡主要由建筑物倒塌破坏，使人员被压、被困造成。地震次生灾害包括火灾、水灾和危险物品等因素对人员死亡也有重要影响[6—8]。这些因素众多，在地震灾难中表现得错综复杂，具有高度不确定性，由此导致由发生原因正向估算震后死亡人数异常困难。我国从“九五”开始，开展了城市震害预测和防御对策研究项目，结合地理信息系统（GIS）等工具进一步估计出人员死亡在预测单元内的大概分布，为灾难救援工作提供了科学依据。这种人员死亡评估对把握地震总体损失有重要作用，但许多导致人员死亡的不确定因素如地震发生的强度、时间、地点、建筑物抗震强度、各种次生灾害等因素，使得各个震例都有自己的特性，用公式或模型理论计算的结果很难符合每一个震例的真实情况。 　　1990年尹之潜等[9]基于建筑物倒塌导致的死伤人数估计模型，精度不够；而2009年吴新燕等[10]发表的地震报道死亡人数随时间变化的修正指数模型，时效性不强，起不到早期预测的目的。2011年李晓杰[11]利用修正或确定地震影响场的各种方法建立人员损失的回归模型，李帆等[3]和宁宝坤等[4]根据震后1到2 d得到的死亡人数报道和震级，利用力学关系式估算地震死亡人数的分布，虽然具有动态性和区域性，但可操作性不强。本文从10个典型的强震巨灾的上报数据可知， T 时的平均时间是震后10 h左右，故本方法在震后8～12 h即可推测死伤总人数。由于受伤总数与死亡总数呈一定的函数关系式，这将在下一步进行研究。并且震后2 T 时间即16～24 h可知死亡人员的80%，所以比较准确，加上后面的幂函数推断，精度增加。 该方法尚有不足，人员死亡时增量 K 值反映实时信息统计的快慢、准确与否。实时信息统计又受很多因素影响，如统计和救援力量薄弱，上报数据可能不准确，