一种山洪灾害人口损失评估模型的构建方法及其应用

摘要

本发明提供一种山洪灾害人口损失评估模型的构建方法，该方法包括获取待评估区域地形参数、土壤参数、土地资源参数、区域水文特征参数和人口分布数据；获取该待评估区域单次历史山洪灾害的气象参数和死亡人口分布数据；基于以上参数确定该次历史山洪灾害的最大淹没深度分布、最大流速分布和最大冲击力分布，计算得到最大分布乘积数据；根据最大分布乘积数据的数值大小将待评估区域划分为多个山洪影响区；确定各个山洪影响区的人口死亡率；根据各个山洪影响区的最大分布乘积数据和人口死亡率，通过拟合构建该待评估区域山洪灾害人口损失评估模型。该模型兼顾考虑了水深、流速和冲击力三种因素共同对山洪灾害导致的人口死亡的影响，精度较高。

说明书

技术领域

本发明涉及山洪灾害的风险评估技术领域。更具体地，涉及一种山洪灾害人口损失评估模型的构建方法。

背景技术

山洪灾害(Flash Floods)指山区小流域由降雨引起的突发的、暴涨暴落的洪水及由山洪诱发的泥石流、滑坡等对国民经济和人民生命财产造成损失的灾害(张平仓等，2009)。人类“依水而居”的特性使其成为山洪灾害最敏感、最脆弱的承灾体，山洪灾害影响人口占所有灾害影响人口的52％(Jonkman，2005)。随着人口和财富的高度集中，使得人口受山洪灾害影响与日俱增、风险加剧(Hanson et al，2011；Peduzzi et al，2012)。过去40年世界人口增长了87％，而山洪易发区的人口增长了114％(Velasquez et al，2012)。IPCC第五次评估报告指出，今后强降水等极端气候事件的出现将更加频繁，未来山洪灾害人口风险将进一步增大(IPCC，2013)。因此，如何从科学意义上认识山洪灾害对人口的影响，定量评估山洪灾害人口损失，成为山洪灾害风险评估研究亟待解决的重要议题。

近十年来，众多水文模型(Pilgrim et al.,1988；Vélezand Francés,2005；Braud et al.,2010)和水动力模型(Bates et al.,2006；Schumann et al.,2009；Murillo and García-Navarro,2010；Lacasta et al.,2014)被应用于山洪灾害的危险性评估中，但针对特定时空尺度下的山洪灾害事件，由于缺少有效的降雨与流量等气象水文数据，加之其迟滞时间短，使得已有水文、水动力模型对山洪进行模拟时表现出高度不确定性(Segura-Beltran et al.，2016)。因此，需结合详实的灾后实地调研来进行山洪建模，以便更好地理解山洪灾害致灾过程。如在野外实地调查的基础上，研究者利用LIQUID水文模型平台对2002 年发生在法国南部Gard地区的山洪灾害进行了降水-径流过程模拟，并着重分析了降雨特征及土壤湿度变化对山洪洪峰流量的影响(Manus et al.，2008； Sandrineet al.，2010；Isabelle Braud，2010)。Roux(2011)在充分考虑坡度因素的影响下，构建了MARINE分布式水文模型用于山洪灾害洪峰流量的预测。Miao(2016)利用GBHM(Geomorphology-based Hydrological Model)模型(GBHM)完成了在中国湿润、半湿润、半干旱山区，无资料支持小流域 (493-1601km

脆弱性曲线是最直观且实用的人口损失评估模型。其构建方法有两种，一是通过多个历史灾情数据；二是通过单次历史事件数据。第一种方法，如 Brown and Graham(1988)基于1950年以来发生在美国和欧洲的24例大型溃坝洪水和洪水灾情数据，构建了“可出逃时间-人口死亡率”脆弱性曲线。针对欧洲洪水，研究者们认为人口损失是洪水特征(水深、流速等)、位置特征(人在建筑物内外、房屋特征等)、人口特征(年龄、健康状况等)综合作用的结果，由此提出了人口损失评估的概念模型，并结合专家经验，提出了人口受伤数量和死亡数量的计算公式(Reiter，2001；Ramsbottom et al.，2003； Penning-Roswell etal.，2005)。Boyd et al.(2005)基于发生在美国的7次历史飓风山洪的灾情数据，构建了“水深-人口死亡率”脆弱性曲线。美国FEMA (Federal Emergency Management Agency)构建的HAZUS系统，提供了基于历史案例构建的“水深-人口死亡率”脆弱性曲线(HAZUS,2006)。Jonkman(2007) 基于1953年荷兰洪水事件、1953年英国山洪事件及1959年日本山洪事件3 次事件的灾情数据，构建了山洪流速较大区域的“水深-人口死亡率”曲线。Brazdovaet al.(2014)针对中欧地区的洪水灾害案例，构建了人口死亡数量与物质损失、洪水管理因素(包括对洪水危险性的认知、洪水经验、洪水管理相关文件、随水文预报的反应、对洪水预警的反应、撤离/应急行动等8个指标)的多元回归模型。而基于单次历史事件构建脆弱性曲线法，则多用于历史灾情数据相对匮乏的研究中。如基于1953年发生在荷兰的洪水事件，Duiser (1989)通过官方新闻报道获取了此次事件的人口死亡及水动力特征数据，构建了“水深-人口死亡率”脆弱性曲线；Waarts(1992)通过进一步完善Duiser(1989)的数据库，在充分考虑预警、出逃、流速、建筑物崩塌因素的基础上，构建了“水深-人口死亡率”脆弱性曲线；Vrouwenvelder and Steenhuis(1997) 在Waarts(1992)构建曲线的基础上，在自变量中加入水位上升速率变量，构建了“水深-水位上升速率-人口死亡率”的三维脆弱性曲线模型；Jonkman (2001)则在前人研究的基础上，经充分分析人员死亡原因后，分别构建了“水深-人口死亡率”曲线和“流速-人口死亡率”曲线。Boyd(2005)分析了1965年9月发生在美国新奥尔良地区飓风洪水事件人员死亡情况，并根据有限的灾情数据，构建了线性“水深-人口死亡率”曲线。Jonkman(2007) 利用SOBEK-1D2D水动力模型，模拟了美国新奥尔良地区2005年卡特里娜飓风影响的洪水灾害事件的水深、流速、水深*流速、水位上升速率、洪水到达时间等洪水特征指标的空间分布情况，并结合人口死亡数据，构建了“水深-人口死亡率”、“水位上升速率-人口死亡率”曲线。

总体来看，前述利用单次洪水历史事件构建的“水深-人口死亡率”曲线、“流速-人口死亡率”曲线以及“水位上升速率-人口死亡率”曲线，都只能反映单一因素如水深、流速或水位上升速率对人口死亡率的影响。但是山洪灾害过程中，水深、流速、冲击力在造成人口死亡时，是共同作用的，因此，上述方法得到的模型用于评估山洪灾害人口损失时，存在精度不高的问题。

因此，需要提供一种相对准确地山洪灾害人口损失评估的方法。

发明内容

本发明的一个目的在于提供一种山洪灾害人口损失评估模型的构建方法，该方法包括：

获取待评估区域地形参数、土壤参数、土地资源参数、区域水文特征参数和人口分布数据；

获取该待评估区域单次历史山洪灾害数据，包括气象参数和山洪死亡人口分布数据；

基于所述地形参数、土壤参数、土地资源参数、区域水文特征参数和气象参数，确定该次历史山洪灾害的最大淹没深度分布、最大流速分布和最大冲击力分布；

计算该次历史山洪灾害的最大淹没深度分布、最大流速分布和最大冲击力分布的乘积，得到最大分布乘积数据；

根据最大分布乘积数据的数值大小将待评估区域划分为多个山洪影响区；

分别确定各个山洪影响区的人口总数、死亡人口数和人口死亡率；

根据各个山洪影响区的最大分布乘积数据和人口死亡率，通过拟合构建该待评估区域山洪灾害人口损失评估模型。

优选地，该方法进一步包括根据所述地形参数、土地资源参数、土壤参数、区域水文特征参数、气象参数，基于HEC-HMS半分布式水文模型，确定逐时径流数据；

根据确定的逐时径流数据，基于FLO-2D水动力模型，确定该次历史山洪灾害的最大淹没深度分布、最大流速分布和最大冲击力分布。

优选地，该方法进一步包括根据所述地形参数、土地资源参数、土壤参数、区域水文特征参数、气象参数确定产流数据、汇流数据和河道演进数据，由得到的产流数据、汇流数据和河道演进数据，确定逐时径流数据。

优选地，所述地形参数包括有坡度；所述土地资源参数包括有土地利用类型；所述土壤参数包括土壤类型、土壤质地、土壤前期湿润状况中的一项或多项；所述区域水文特征参数包括流域单元面积、河道断面宽度、河道断面水深、河底比降、河段长度、初始基流量中的一项或多项；所述气象参数包括总降雨量、降雨时长中的一项或两项。

优选地，该方法进一步包括根据坡度、土地利用类型、土壤类型、土壤前期湿润状况、土壤质地、总降雨量，采用SCS-CN法，确定产流数据。

优选地，该方法进一步包括根据坡度、土地利用类型、土壤类型、土壤前期湿润状况、土壤质地、流域单元坡度、流域单元面积和降雨时长，采用 SCS单位线法，确定汇流数据。

优选地，该方法进一步包括根据河道断面宽度、河道断面水深、河底比降、河段长度、初始基流量，采用马斯京根法，确定河道演进数据。

优选地，该方法进一步包括根据所述待评估区域的建筑物分布及其建筑面积数据和人均住房面积，确定人口分布数据。

本发明第二个目的在于提供一种特定气象条件下山洪灾害人口损失评估方法，该方法包括：

利用如权利要求1-8任一项定义的方法，构建待评估区域的山洪灾害人口损失评估模型；

获取该待评估区域的预报或实时气象条件的气象参数；

基于该待评估区域的地形参数、土地资源参数、土壤参数、区域水文特征参数和所述气象参数，确定本次气象条件的最大淹没深度分布、最大流速分布和最大冲击力分布；

计算本次气象条件的最大淹没深度分布、最大流速分布和最大冲击力分布的乘积，得到本次最大分布乘积数据；

根据计算得到的本次最大分布乘积数据、所述待评估区的人口分布数据和所述待评估区域的山洪灾害人口损失评估模型，评估本次气象条件的人口死亡数。

本发明的有益效果如下：

本发明提供的山洪灾害人口损失评估模型的构建方法是以单次山洪历史事件进行的，可用于历史灾情数据相对匮乏的地区，得到的人口脆弱性曲线是最直观且实用的人口损失评估模型。根据本发明提供的山洪灾害人口损失评估模型的构建方法得到的人口损失评估模型，兼顾考虑了水深、流速和冲击力三种因素共同对山洪灾害导致的人口死亡的影响，精度较高。利用该模型，可以根据待评估区域地形参数、土地资源参数、土壤参数、区域水文特征参数、人口分布数据以及预报或实时气象条件的气象参数，得到本次气象条件的最大淹没深度分布、最大流速分布和最大冲击力分布，进而评估本次气象条件的人口死亡数，为开展相应的预防、救援措施做准备。

附图说明

下面结合附图对本发明的具体实施方式作进一步详细的说明。

图1示出本发明的山洪灾害人口损失评估模型的构建流程图。

图2示出实施例中台头沟流域位置及其地形特征。

图3示出实施例中由逐时降雨数据产生的逐时径流数据，即流量过程线。

图4示出实施例中台头沟流域的最大分布乘积数据。

图5示出实施例中台头沟流域中的死亡失踪人口分布点位图。

图6示出实施例中台头沟流域居民住宅分布及人口分布(左侧为居民住宅分布；右侧为人口分布)。

图7示出实施例中台头沟流域的山洪灾害的人口脆弱性曲线。

图8示出实施例中台头沟流域的DEM数据。

图9示出实施例中台头沟流域的坡度。

图10示出实施例中台头沟流域的土地利用类型。

图11示出实施例中台头沟流域的土壤类型。

图12示出实施例中台头沟流域的土壤质地(粉砂土)的分布。

图13示出实施例中台头沟流域的土壤质地(砂土)的分布。

图14示出实施例中台头沟流域的土壤质地(粘土)的分布。

具体实施方式

为了更清楚地说明本发明，下面结合优选实施例和附图对本发明做进一步的说明。附图中相似的部件以相同的附图标记进行表示。本领域技术人员应当理解，下面所具体描述的内容是说明性的而非限制性的，不应以此限制本发明的保护范围。

下面将结合图1具体说明本发明的山洪灾害人口损失评估模型的构建方法的流程图。

需要说明的是，本发明中的死亡人口指的是可实际确认山洪灾害中死亡的人口以及山洪灾害中失踪的人口。

图1的S110步骤中，获取待评估区域地形参数、土地资源参数、土壤参数、区域水文特征参数和人口分布数据；获取该待评估区域单次历史山洪灾害数据，包括气象参数和山洪死亡人口分布数据；

本领域技术人员公知，影响山洪灾害的发生的因素有地形参数，如坡度；土地资源参数，如土地利用类型；土壤参数，如土壤类型、土壤质地、土壤前期湿润状况等；区域水文特征参数，如流域单元面积、河道断面宽度、河道断面水深、河底比降、河段长度、初始基流量；以及气象参数，如总降雨量和降雨时长。

其中，土地利用类型、土壤类型、土壤质地、初始基流量可通过查询相关资料得到；河道断面水深则可以通过实地调研得到；河道断面宽度则可以结合遥感图像，经过简单计算得到；坡度、流域单元面积、河段长度、河底比降则可以经过通过数字高程数据(DEM数据)经过计算得到；总降雨量、降雨时长、土地前期湿润状况则可以根据气象站点的逐时降雨数据，经过处理得到。

人口分布数据则可通过遥感影像统计建筑分布及其建筑面积，结合人均建筑面积，得到人口分布数据。

历史山洪灾害中的山洪死亡人口分布数据则可通过实际走访获得待评估区域内死亡人口的具体位置，获得其分布情况。

S120步骤中，根据地形参数、土地资源参数、土壤参数、区域水文特征参数和气象参数，确定该次历史山洪灾害的最大淹没深度分布、最大流速分布和最大冲击力分布，计算该次历史山洪灾害的最大淹没深度分布、最大流速分布和最大冲击力分布的乘积，得到最大分布乘积数据。

在具体的实施方式中，根据地形参数、土地资源参数、土壤参数、区域水文特征参数、气象参数，基于HEC-HMS半分布式水文模型，可以确定逐时径流数据；由逐时径流数据、曼宁粗糙系数和FLO-2D水动力模型确定该次历史山洪灾害的最大淹没深度分布数据、最大流速分布数据和最大冲击力分布数据。

其中，HEC-HMS半分布式水文模型是对山洪灾害事件的降雨-径流过程进行重建，而FLO-2D水动力模型则是对山洪演进过程进行重建。

利用HEC-HMS半分布式水文模型进行降雨-径流过程重建时，主要进行了产流计算、汇流计算及河道演进模拟，最终得到了逐时径流数据，即流量过程线。

在具体的实施过程中，在获得坡度、土地利用类型、土壤类型、土壤前期湿润状况、土壤质地、总降雨量的基础上，产流计算过程可采用SCS-CN 法。该方法对输入数据的要求较低、模拟精度高且能全面反映流域下垫面自然条件对可收集雨水资源量的影响。具体的计算过程如下：

I

式(1)-(4)中，P

在具体的实施过程中，在获得坡度、土地利用类型、土壤类型、土壤前期湿润状况、土壤质地、流域单元坡度、流域单元面积和净降雨历时的情况下，采用SCS单位线法，确定汇流数据。SCS单位线是一个三角形单位线，由降雨开始到洪峰出现的时间T

T

式(5)-(9)中，Q

在具体的实施过程中，考虑到更好的与FLO-2D模型进行耦合，选择马斯京根法确定河道演进数据，需要提前测量得到河道断面宽度、河道断面水深、河底比降、河段长度、初始基流量和曼宁粗糙系数。具体计算公式如下：

S(t)＝KQ′(t)＝K(xI(t)+(1-x)Q(t)) (11)

式(10)-(11)中，K为槽蓄系数，即稳定流情况下河段传播时间，具有时间因次，单位为h；x为无因次的流量比重因子，无量纲；K、x统称为演进系数； S(t)为河段t时刻的槽蓄水量，单位为m

Q

式(12)-(18)中，I

根据上述计算过程得到产流、汇流以及河道演进过程，确定由逐时降雨数据产生的逐时径流数据，即流量过程线。将流量过程线引入FLO-2D模型，将 HEC-HMS半分布式水文模型和FLO-2D模型两个模型进行耦合，在相同的连续方程控制下，FLO-2D模型则侧重流体的运动，故而用流体运动时的动量守恒方程(如式(20)-(21)所示)替换HEC-HMS模型中侧重计算径流槽蓄控制方程，并结合流域内的曼宁粗糙系数n，完成整个山洪事件的实时演进过程模拟，所得结果包括该次历史山洪灾害中的最大淹没深度分布数据、最大流速分布数据、最大冲击力分布数据。具体计算过程如下：

式(19)-(21)中，h是水深，单位为m；Sy是摩擦坡降，无量纲；i为河底比降，无量纲；g是重力加速度，单位为m/s

步骤S130中，接着根据该次历史山洪灾害的最大分布乘积数据的数值大小，将待评估区域划分为多个山洪影响区。

步骤S140中，根据各个山洪影响区的人口分布数据，确定各个山洪影响区的总人口；结合各个山洪影响区的死亡人口，得到各个山洪影响区的人口死亡率。

在具体的实施方式中，可通过待评估区域内的遥感影像数据，得到待评估区域内的建筑物分布及建筑面积信息；然后结合当地的人均住房面积数据，计算每栋建筑物内的居住人口，进而绘制待评估区域内的人口分布图；接着根据待评估区域内最大分布乘积数据，将待评估区域划分为多个山洪影响区；由各个山洪影响区的人口分布数据，确定各个山洪影响区的总人口；再根据各个山洪影响区的总人口和死亡人口，得到该次山洪灾害中各个山洪影响区的人口死亡率。

在优选的实施方式中，确定建筑物分布及建筑面积信息时，可通过实地调研对获取的结果进行校正，校正内容包括(1)建筑物分布校正，去除无人居住建筑物，(2)建筑面积校正，部分富裕家庭的住宅为二层楼房，其建筑面积需在人工目视解译获取结果的基础上翻倍，(3)建筑物类型校正，只保留居民住宅类型。

在优选的实施方式中，结合当地的人均住房面积数据，得到待评估区域内的人口分布图后，可通过比对统计得到的灾害影响区内的总人口数和实际人口数，验证所绘制的人口分布图的准确性，并进行人口分布的抽样调查，完成对所绘制的人口分布图的校正工作并获得最终的较精准的人口分布图。

图1步骤S150中，根据各个山洪影响区的最大分布乘积数据和山洪影响区的人口死亡率建立山洪灾害人口损失评估模型。

本发明第二个目的在于提供一种特定气象条件下山洪灾害人口损失评估方法，包括以下步骤：

利用本发明第一个目的提供的山洪灾害人口脆弱性模型构建方法，构建待评估区域的山洪灾害人口损失评估模型；

获取该待评估区域的预报或实时气象条件的气象参数；

基于该待评估区域的地形参数、土壤资源参数、土壤参数、区域水文特征参数和所述气象参数，确定本次气象条件的最大淹没深度分布、最大流速分布和最大冲击力分布；

计算本次气象条件的最大淹没深度分布、最大流速分布和最大冲击力分布的乘积，得到本次最大分布乘积数据；

根据计算得到的本次最大分布乘积数据、待评估区的人口分布数据和所述待评估区域的山洪灾害人口损失评估模型，评估本次气象条件的人口死亡数量。

在具体的实施过程中，由于各个待评估区域的地形参数、土地资源参数、土壤参数、区域水文特征参数等不同，且人口分布在密度、居住习惯等方面亦不同，因此，针对不同的待评估区域，都需要先根据单次历史山洪灾害事件建立起山洪灾害人口损失评估模型。

此外，针对某一待评估区域，其地形参数、土地资源参数、土壤参数、区域水文特征参数以及人口分布在一定时间内基本是不改变的，可通过数据查询或计算得到。但是每次山洪灾害中的气象参数之间的差异非常大，且只能够在山洪灾害发生前的一段有限时间内获得，因此在获得预报或实时气象条件的气象参数时，便可得到该待评估区域在本次气象条件下的最大淹没深度分布、最大流速分布和最大冲击力分布；进而得到本次气象条件的最大分布乘积数据；

根据计算得到的本次气象条件最大分布乘积数据、待评估区的人口分布数据和已得到的所述待评估区域的山洪灾害人口损失评估模型，便可以评估本次气象条件的人口死亡数，为开展相应的预防、救援措施提供参考，以有效保护国家和个人生命财产安全。

下面将结合本发明优选实施例具体说明根本发明的山洪灾害人口损失评估模型的构建方法。

实施例

下面将以我国河北省井陉县台头沟流域洪灾为例，具体说明本发明的山洪灾害人口损失评估模型的构建方法。在该实施例中，台头沟流域即为待评估区域。

台头沟流域地处太行山东麓，河北省西陲井陉县，与山西省接壤(如图2 所示)。其流域面积为31km

流域位于井陉县，据统计，流域从1989年至今，先后发生过9次山洪灾害事件，均造成严重损失，以2016年7月19日发生的特大暴雨山洪灾害的损失最为惨重，即“7·19”特大暴雨山洪事件。

该次山洪灾害共造成台头沟流域26人死亡失踪，1000余间房屋不同程度毁坏，1500余亩农田受灾，200余辆各色车辆被冲走或淹坏，50余家商店被洗劫一空，物质席卷无数，企业全部停产，水、电、路、通信设施全部瘫痪。

降雨-径流过程重建

关于该次山洪灾害事件的人口损失评估模型的构建过程中，首先获取该次山洪灾害中的地形参数、土地资源参数、土壤参数、区域水文特征参数、气象参数、人口分布数据和山洪死亡人口分布数据。

通过河北地理信息局、地理国情监测云平台等，搜集或经进一步计算，获取台头沟的地形参数、土地资源参数、土壤参数以及区域水文特征参数，具体如下：DEM数据如图8所示，经过GeoHEC-HMS软件处理，可得坡度数据如图9所示，以及流域单元面积、河底比降、河段长度；土地利用类型如图10所示，土壤类型如图11所示，土壤质地如图12-14所示。

而气象参数则可以通过流域周边气象监测站点获得，此次“7·19”特大暴雨山洪事件中的降水过程三个阶段：

第一阶段18日23时到19日16时，共历时17小时，是降水的开始阶段。主要天气影响系统是高空槽和地面倒槽。降雨前期，亚洲中高纬为两脊一槽型，低槽位于河套西部，东部的大陆高压和副热带高压叠加，并不断加强，成阻挡形势。从低槽动态看，低槽东移缓慢，且不断加深。槽前强劲的西南气流向河北省源源不断地输送水汽。此时低层850hPa在河北省西部太行山附近有一南北向的切变线，且有东南气流箱河北省输送水汽；而高层200hPa河北省上空为明显的疏散气流，地面上河北省南部受低压倒槽影响，倒槽的辐合线正位于河北省西南部太行山附近。高空辐散、低层辐合有利于垂直运动的加强。分析邢台探空亦可见，低层湿层厚，有利降水；有下暖湿、上干冷的不稳定层结条件，有一定的不稳定能量，且抬升凝结高度很低，极易触发不稳定能量释放；低层偏东风持续并不断加强，地形抬升作用有利于太行山前降水加强；以上形势为河北省西部太行山暴雨的有利形势。

第二阶段从19日17时到21日0时，历时32小时，是“7.19”过程的主要降雨阶段，阶段中降雨的强度、强降雨的范围以及降雨的累计雨量均是本次过程中最强、最大的阶段。该阶段主要影响系统是高空低涡和地面气旋。影响第一阶段降水的高空低槽缓慢东移并不断加深，19日20时在山西南部切出低涡，低涡不断发展，并向东北方向移动。低涡生成初期，虽然中心未在河北省，但其北部的偏东气流不断加强，强劲的偏东气流在太行山迎风坡强迫抬升，使河北省西部太行山山前的强降雨持续并增强。随着低涡向东北向移动，与其对应的地面气旋20日02时进入河北省东南部并缓慢向北移动，气旋本身的动力条件以及气旋前部强劲的东南气流的水汽输送，使河北省东部平原出现另一强降雨中心，该降雨中心位于气旋中心附近暖锋北侧，并随着气旋的移动向北移动。

第三阶段从21日1时到14时，历时14小时，是“7.19”过程的尾声阶段，该阶段高空低涡减弱向东北方向移动，地面气旋减弱消失，河北省东北地区有降水，但降雨区域、降水强度都明显低于第二阶段，降水趋于结束。

根据上述降水过程及流域附近气象监测站点数据，可以得到“7.19”特大暴雨山洪事件中的总降雨量为611mm，降雨时长为49h。

结合上述收集得到的参数，可知，台头沟流域的其他相关参数，如反映降雨前流域特征的一个综合参数CN

将上述参数输入式(5)-(18)可得台头沟暴雨山洪事件期间由逐时降雨数据产生的逐时径流数据，即流量过程线，如图3所示。据图3可知此次降雨过程在19日晚10时左右降雨量最大，对应的径流数据也在几乎同一时间达到流量的最大值，即洪峰流量值。据野外实地调查，采访亲历者得知，此次山洪事件在19日晚10时左右降雨最为猛烈，产生的洪流最为强劲，水位最高处可达8米，间接验证模型模拟效果的精确性。

山洪演进过程重建

由于HEC-HMS模型的输出结果，仅为流量过程线，无法提供表征此次事件致灾强度指标的实时空间分布情况，为此，将经HEC-HMS半分布式水文模型模拟估算的流量过程线引入到FLO-2D模型，将两个模型进行耦合，在相同的连续方程控制(如式(19)所示)下，FLO-2D模型则侧重流体的运动，故而用流体运动时的动量守恒方程(如式(20)-(21)所示)，替换HEC-HMS 模型中侧重计算径流槽蓄控制方程，并结合流域内的曼宁粗糙系数，完成整个山洪事件的实时演进过程模拟，得到最大淹没深度×最大流速×最大冲击力的分布数据，如图4所示。

死亡失踪人口调查利用前期从中国民政府减灾中心获得的初步死亡失踪人口数据，研究人员亲赴研究区对死亡失踪人口展开详细调查，调查结果如表1和图5所示。

表1死亡失踪人员概况

建筑物提取及人口分布制图

研究区内建筑物的提取主要通过以下步骤获得：首先，获取台头沟特大暴雨山洪灾害事件发生前台头沟流域遥感影像数据；接着，通过人工目视解译的方式获取流域内建筑物分布及建筑面积信息；然后，通过实地调研对获取结果进行校正，校正内容包括(1)建筑物分布校正，去除无人居住建筑物，(2) 建筑面积校正，部分富裕家庭的住宅为二层楼房，其建筑面积需在人工目视解译获取结果的基础上翻倍，(3)建筑物类型校正，由于灾害发生时间为夜间，居民全部在自己住宅内，故本研究仅关心居民住宅，并将其提取出来，最终获取研究区内居民住宅的分布及建筑面积数据，此处仅以研究区内的台头村为例进行展示(如图6左侧所示)。

在获取研究区精准的建筑物分布数据的基础上，绘制研究区人口分布图，具体步骤如下：首先，通过《河北省经济统计年鉴》(2017)获取2016年河北省农村人均住房面积数据；接着，结合已获得的房屋分布及其建筑面积数据，计算每栋建筑物内的居住人口数，并统计研究区内的总人口数；然后，将统计获得的研究区内总人数与实地调查获得的总人数据进行比对，验证所绘制的人口分布图的准确性，并进行人口分布的抽样调查，完成对所绘制的人口分布图的校正工作并获得最终的较精准的人口分布图，此处仅以研究区内的台头村为例进行展示(如图6右侧所示)。

人口损失评估模型构建

本发明中的人口损失评估模型是指由该次历史山洪灾害的最大淹没深度分布、最大流速分布和最大冲击力分布的乘积，即最大分布乘积数据和人口死亡率建立起来的，其中，人口死亡率是指不同山洪影响区内的死亡和失踪人口总数与该区域内总人口数的比值。

人口损失评估模型的构建过程，具体为：(1)根据最大分布乘积数据的数值大小将台头沟流域等分成31个山洪影响区，即0-30万Nm/s、30～60万Nm/s、60～90万Nm/s、90～120万Nm/s、…、870～900万Nm/s、>900万Nm/s；(2) 根据各个山洪影响区的人口分布数据，获得各个山洪影响区各自的总人口数和各个山洪影响区各自的死亡人口数；进而获得各个山洪影响区的人口死亡率；(3)将各个山洪影响区的最大分布乘积数据与人口死亡率进行匹配，最终拟合构建得到山洪灾害人口损失评估模型。

为使最大分布乘积数据与人口死亡率匹配过程中更加精确，本发明还采用了结合死亡失踪人口的位置信息，提取各个位置的最大分布乘积数据，并计算出各个山洪影响区的死亡人口位置的最大分布乘积数据均值。然后再将各个山洪影响区计算得到的均值与人口死亡率进行匹配，拟合得到山洪灾害人口损失评估模型。

经过上述构建过程， 本发明得到的“7·19”特大暴雨山洪事件的人口损失评估模型如表2所示，其中M

表2山洪灾害人口损失评估模型

显然，本发明的上述实施例仅仅是为清楚地说明本发明所作的举例，而并非是对本发明的实施方式的限定，对于所属领域的普通技术人员来说，在上述说明的基础上还可以做出其它不同形式的变化或变动，这里无法对所有的实施方式予以穷举，凡是属于本发明的技术方案所引伸出的显而易见的变化或变动仍处于本发明的保护范围之列。