CN105740607B - 一种溃坝洪水致生命损失的计算方法 - Google Patents
一种溃坝洪水致生命损失的计算方法 Download PDFInfo
- Publication number
- CN105740607B CN105740607B CN201610047890.0A CN201610047890A CN105740607B CN 105740607 B CN105740607 B CN 105740607B CN 201610047890 A CN201610047890 A CN 201610047890A CN 105740607 B CN105740607 B CN 105740607B
- Authority
- CN
- China
- Prior art keywords
- dam
- human loss
- submodule
- flood
- human
- 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.)
- Active
Links
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16Z—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS, NOT OTHERWISE PROVIDED FOR
- G16Z99/00—Subject matter not provided for in other main groups of this subclass
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种溃坝洪水致生命损失的计算方法,其特征在于,包括:生命损失数据选择模块、生命损失影响因素模块、生命损失子模块以及耦合模型这4个递进计算模块。本发明提供的一种溃坝洪水致生命损失的计算方法,筛选了对溃坝洪水致生命损失的重要影响因素,建立了致灾因子、孕灾环境、承灾体和救灾能力4个生命损失子模块,并建立了生命损失耦合模型;本发明可应用于中国各类溃坝条件下的生命损失和生命损失率计算;可对已溃坝水库进行生命损失灾后评价;可对未溃坝水库进行生命损失风险预测;为灾后重建以及灾前防范预测提供了科学的计算方法。
Description
技术领域
本发明涉及一种溃坝洪水致生命损失的计算方法,特别是溃坝引起的洪水与受灾群体生命损失数量之间的关系以及计算方法,属于水灾害技术领域。
背景技术
中国有2500多年的筑坝史,是人类筑坝历史最悠久的国家之一,也是当今世界拥有水库数量最多的国家。中国目前已建成水库97721座,总库容达8298亿m3,其中大型水库687座,中型水库3774座,小型水库93260座。水库在蓄水、发电、防洪、航运、养殖等方面发挥着巨大的效用,但是水库一旦溃坝,就会对人类生命及自然环境造成灾难性的破坏。1954年至2013年共60年间,中国发生溃坝事故3 544起,1954~1960年,年均溃坝50座,1971~1980年共溃坝2031座,年均溃坝达到峰值203座,70年代以后,年均溃坝数迅速下降,90年代年均溃坝达23座,2011~2013年,年均溃坝降至4座。溃坝事故数量虽大幅下降,却持续剥夺人类的生命,1993年青海沟后水库溃坝,导致320人丧生;2010年吉林大河水库溃坝,导致31人遇难。
中国水库溃坝生命损失研究刚起步发展,大多在国外方法上进行系数修正,但国外方法中使用的溃坝事故几乎不含中国案例,适用性较差,计算精度低。本发明立足中国国情,加深对生命损失的研究,通过对历史溃坝事故的分析,优选生命损失影响因素、建立生命损失子模块以及耦合模型,得到一套完整的溃坝洪水致生命损失的计算方法。
发明内容
本发明所要解决的技术问题是,提供一种溃坝洪水致生命损失的计算方法,该方法可应用于中国各类溃坝条件下的生命损失和生命损失率计算,并能对溃坝水库施予生命损失灾后评价,以及对未溃坝水库(现有水库或规划水库)进行生命损失风险预测,为我国溃坝洪水导致的生命损失原因分析以及生命损失风险预测提供了可靠的数据和方法支撑,为灾后重建以及灾前防范预测提供了科学的计算方法。
为解决上述技术问题,本发明采用的技术方案为:
一种溃坝洪水致生命损失的计算方法,其特征在于,包括:生命损失数据选择模块、生命损失影响因素模块、生命损失子模块以及耦合模型这4个递进计算模块;
所述生命损失数据选择模块以中国溃坝事故为数据材料,采用分层抽样法,均衡考虑水平区域层、垂直海拔层和时间年代层,进行三维时空抽样,并综合考虑各类溃坝条件,优选溃坝案例;
所述生命损失影响因素模块选取的生命损失影响因素包含溃坝洪水严重程度SF、溃坝发生时间TB、建筑物损毁ND、受灾区与大坝平均距离DD、风险人口PR、对溃坝理解程度UB、预警时间TW、撤离条件EC、溃坝方式MB、库容SR和溃坝时天气WE;
根据灾害系统三要素——致灾因子DCF、孕灾环境DPE及承灾体DAB,同时考虑救灾能力DRA对溃坝生命损失有着重要影响,将生命损失影响因素分为4类;
第一类:致灾因子DCF为诱发灾害的因素,包括SF、MB、和SR;
SF的计算方法为用计算断面的溃坝峰值流量Qtop与溃坝洪水引起的最大水面宽度Wmax之商作为SF的估计值:SF≈Qtop/Wmax;
从单因素角度考虑,一般SF越大,则生命损失L越大;SF<0.5m2/s,对生命损失影响可以忽略不计;0.5m2/s≤SF<4.6m2/s,SF为低严重性;4.6m2/s≤SF<12.0m2/s,SF为中严重性;SF≥12.0m2/s,SF为高严重性;
MB表示溃坝方式,从单因素角度考虑,一般MB的频率越高,溃坝事故可能性越大,L也越大;
SR表示库容,从单因素角度考虑,一般SR越大,溃坝越严重,L也越大;
第二类:孕灾环境DPE为形成灾害的外部环境,包括TB、WB、ND和DD;
TB表示溃坝发生时间,从单因素角度考虑,从白天工作期逐步进入凌晨睡眠期,人类活动逐步减少,警觉性下降,L的危险性也逐渐增加;
WE表示溃坝时天气,从单因素角度考虑,WE越恶劣,L取值越大;
ND表示建筑物损毁,从单因素角度考虑,ND越多,对人类威胁性越大,L越多;
DD表示受灾区与大坝平均距离,从单因素角度考虑,DD越近,L越多;
第三类:承灾体DAB为受到损害的对象,包括PR和UB;
PR表示风险人口,为下泄洪水或溃坝洪水淹没范围内的人数,PR作为L的基数,对L影响巨大,从单因素角度考虑,PR越多,L相应也越多;
UB表示对溃坝理解程度,UB受当地文化教育程度以及当地溃坝洪水灾害宣传情况影响,UB分为清楚和模糊两类,从单因素角度考虑,UB越模糊,L越多;
第四类:救灾能力DRA为受灾对象的自救或被救的条件,包括TW和EC;
TW表示预警时间,从单因素角度考虑,TW越短,L相应也越多;
EC表示撤离条件,EC受到地形环境的影响,对受灾人员救援疏散起着重要作用,将EC分为好、中、差三类,从单因素角度考虑,EC越差,L越多;
根据第一类~第四类建立致灾因子子模块M1,孕灾环境子模块M2,承灾体子模块M3,救灾能力子模块M4;结合熵权法建立所述生命损失子模块,同时分析生命损失影响因素IFL对所述生命损失子模块的贡献率,从而筛选重要的IFL:
溃坝洪水致生命损失的计算方法的步骤如下:
S01,建立判断矩阵:
式(1)中,xij为IFL初始值,有m个溃坝事件,n个IFL,n在M1~M4中分别取3、4、2、2;
S02,M1~M4子模块标准化处理:
rij=xij/max{xij} (2)
式(2)中,rij代表xij的标准化值;
S03,根据熵的定义,计算M1~M4子模块的信息熵ej:
式(3)中,k=1/lnm,熵值反映了信息的无序化程度,可以用来度量信息量的大小,式(3)中IFL比重为pij:
S04,IFL权重wj计算和筛选:
IFL的效用价值越高,则对M1~M4子模块的重要性就越大,该IFL的权重也就越大;当wj<5%,表明IFL的效用价值太小,删除此IFL,只保留wj≥5%的IFL;
筛选留下SF,TB,ND,DD,PR,TW,EC,SR和WE共8个重要IFL,得到W1=(0.203,0.766),W2=(0.073,0.636,0.254),W3=(0.966),W4=(0.701,0.299);
S05,子模块Mi:
式(6)中Mi代表M1~M4子模块中的任一个;
S06,对fL与子模块Mi进行多元非线性回归分析,建立L耦合模型:
式(7)中,fL代表生命损失率;M1代表致灾因子子模块;M2代表孕灾环境子模块;M3代表承灾体子模块;M4代表救灾能力子模块;
fL=L/PR (8)
式(8)中,L为生命损失,指受到溃坝洪水淹没而遇害的死亡人数;
由式(7)、式(8)可推导出L为:
所述水平区域层包括华东地区、华南地区、华中地区、华北地区、西北地区、西南地区和东北地区;所述垂直海拔层包括库面海拔最低至30m,最高达3100m;所述时间年代层包括从20世纪50年代至21世纪10年代,每个年代选取2座溃坝水库。
所述MB包括4种,按照溃坝频率从高到低分别为漫坝、坝体质量差、管理不当和其他类;所述坝体质量差包括坝体、输水洞和溢洪道由于质量问题而发生渗漏、堵塞;所述管理不当包括超蓄、维护运用不良和无人管理;所述其他类包括白蚁建巢、临危扒口和原因不详。
所述SR包括大一型水库、大二型水库、中型水库、小一型水库和小二型水库;所述大一型水库总库容大于10亿立方米;所述大二型水库总库容在1亿立方米以上而小于10亿立方米;所述中型水库总库容大于或等于1000万立方米以上而小于1亿立方米;所述小一型水库总库容大于或等于100万立方米而小于1000万立方米;所述小二型水库总库容在10万立方米以上而小于100万立方米。
所述TB分为3个时段,白天工作期6:00-18:00、夜晚休息期18:00-22:00;和凌晨睡眠期22:00-次日6:00。
所述WB分为5个等级:一级:暴雨、暴雪、台风、雾霾;二级:大雨、大雪、大风;三级:中雨、中雪;四级:小雨、阵雨、小雪;五级:晴天、多云及阴天。
所述TW分为5段;0-15min、15-30min、30-45min、45-60min和大于60min。
留一交叉分析验证法,对fL与子模块Mi结果进行率定验证得到式(7)。本发明提供的一种溃坝洪水致生命损失的计算方法,筛选了对溃坝洪水致生命损失的重要影响因素—溃坝洪水严重程度、溃坝发生时间、建筑物损毁、受灾区与大坝平均距离、风险人口、预警时间、撤离条件和库容,建立了致灾因子、孕灾环境、承灾体和救灾能力4个生命损失子模块,并建立了生命损失耦合模型;本发明可应用于中国各类溃坝条件下的生命损失和生命损失率计算;可对已溃坝水库进行生命损失灾后评价;可对未溃坝水库(现有水库或规划水库)进行生命损失风险预测;为灾后重建以及灾前防范预测提供了科学的计算方法。
说明书附图
图1为本发明的结构示意图。
具体实施方式
下面结合附图对本发明作更进一步的说明。
如图1所示,一种溃坝洪水致生命损失的计算方法,其特征在于,包括溃坝事故分层抽样(即生命损失数据选择模块),生命损失影响因素(即生命损失影响因素模块)、生命损失子模块以及耦合模型4个递进计算软件。
一种溃坝洪水致生命损失的计算方法,其特征在于,所述方法的建立依次含有以下步骤:
1.以中国溃坝事故为数据材料,采用分层抽样法,均衡考虑水平区域、垂直海拔和时间年代3维时空抽样,并综合考虑各类溃坝条件,优选溃坝案例。
水平区域层:分为全国七大区域:华东地区、华南地区、华中地区、华北地区、西北地区、西南地区和东北地区;
垂直海拔层分布:库面海拔最低至30m,最高达3100m;
时间年代层:从20世纪50年代至21世纪10年代,每个年代选取2座溃坝水库。
2.考虑因素的功能性、可获取性、可比性、完整性、非重叠性以及定量和定性相结合的6大原则,选取的生命损失影响因素包含溃坝洪水严重程度SF、溃坝发生时间TB、建筑物损毁ND、受灾区与大坝平均距离DD、风险人口PR、对溃坝理解程度UB、预警时间TW、撤离条件EC、溃坝方式MB、库容SR和溃坝时天气WE,总共11项生命损失影响因素IFL。
3.根据灾害系统三要素——致灾因子、孕灾环境及承灾体,同时考虑救灾能力对溃坝生命损失有着重要影响,将生命损失影响因素分为4类,并依据生命损失影响因素的不同性质进行划分。
致灾因子DCF为诱发灾害的因素,包含了SF、MB、和SR共3项。
SF表示溃坝洪水严重程度,对L有着重要直接作用。SF一般不能通过调查直接获取其数值,可用计算断面的溃坝峰值流量Qtop与溃坝洪水引起的最大水面宽度Wmax之商作为SF的估计值:
SF≈Qtop/Wmax
从单因素角度考虑,一般SF越大,则L越大。SF<0.5m2/s,对生命损失影响可以忽略不计;0.5m2/s≤SF<4.6m2/s,SF为低严重性;4.6m2/s≤SF<12.0m2/s,SF为中严重性;SF≥12.0m2/s,SF为高严重性。
MB表示溃坝方式,溃坝的主要方式为4种,按照溃坝频率从高到低分别为漫坝、坝体质量差(坝体、输水洞、溢洪道等质量问题发生渗漏、堵塞等)、管理不当(超蓄、维护运用不良、无人管理等)和其他类(白蚁建巢、临危扒口、原因不详等)。从单因素角度考虑,一般MB的频率越高,溃坝事故可能性越大,L也越大。
SR表示库容。根据国家《防洪标准》(GB50201-94),水库工程分为五个等别:大(一)型水库:总库容大于10亿立方米;大(二)型水库:总库容在1亿立方米以上而小于10亿立方米;中型水库:总库容大于或等于1000万立方米以上而小于1亿立方米;小(一)型水库:总库容大于或等于100万立方米而小于1000万立方米;小(二)型水库:总库容在10万立方米以上而小于100万立方米。从单因素角度考虑,一般SR越大,溃坝越严重,L也越大。
孕灾环境DPE为形成灾害的外部环境,包含了TB、WB、ND和DD共4项。
TB表示溃坝发生时间。按照常人作息规律将DPE分为3个时段,白天工作期(6:00-18:00);夜晚休息期(18:00-22:00);凌晨睡眠期(22:00-次日6:00)。从单因素角度考虑,从白天工作期逐步进入凌晨睡眠期,人类活动逐步减少,警觉性下降,L的危险性也逐渐增加。
WE表示溃坝时天气。溃坝发生时,若遇到极端恶劣天气,将不利于逃离或救援,相比较若无风雨将利于逃离和救援。本文按照天气恶劣程度将WB分为5个等级:一级:暴雨、暴雪、台风、雾霾等极端天气;二级:大雨、大雪、大风等;三级:中雨、中雪等;四级:小雨、阵雨、小雪等;五级:晴天、多云及阴天。从单因素角度考虑,WB越恶劣,L取值越大。
ND表示建筑物损毁。洪水冲毁的建筑物越多,相应避难处越少,且水流中携带建筑物的残骸硬物越多,居民被建筑物碰伤砸死的可能性增加。从单因素角度考虑,ND越多,对人类威胁性越大,L越多。
DD表示受灾区与大坝平均距离。受灾处与大坝溃口的距离影响着洪水严重程度,也个关系着居民反应逃离时间,间接作用着L大小。从单因素角度考虑,DD越近,L越多。
承灾体DAB为受到损害的对象,包含了PR和UB共2项。
其中,PR表示风险人口,为下泄洪水或溃坝洪水淹没范围内的人数。PR作为L的基数,对L影响巨大。从单因素角度考虑,PR越多,L相应也越多。
UB表示对溃坝理解程度。UB受当地文化教育程度以及当地溃坝洪水灾害宣传情况影响。本发明将UB分为清楚和模糊两类,从单因素角度考虑,UB越模糊,L越多。
救灾能力DRA为受灾对象的自救或被救的条件,包含了TW和EC共2项。
TW表示预警时间。溃坝洪水多为突发性事件,争取更多TW会大大降低L,因而TW就是挽救生命至关重要的黄金时间。本发明按照TW长短分为5段,0-15min,15-30min,30-45min,45-60min,>60min。从单因素角度考虑,TW越短,L相应也越多。
EC表示撤离条件。EC受到地形环境的影响,对受灾人员救援疏散起着重要作用。将EC分为好、中、差三类,从单因素角度考虑,EC越差,L越多。
4.建立4个生命损失子模块:致灾因子子模块M1,孕灾环境子模块M2,承灾体子模块M3,救灾能力子模块M4。结合熵权法建立L子模块,同时分析IFL对子模块的贡献率,从而筛选重要的IFL:
建立判断矩阵:
其中,xij为IFL初始值,有m个溃坝事件,n个IFL,n在4次独立的子模块计算下分别取3、4、2、2。
子模块标准化处理:
rij=xij/max{xij}
根据熵的定义,计算子模块的信息熵ej:
其中,k=1/lnm。熵值反映了信息的无序化程度,可以用来度量信息量的大小。式中IFL比重pij:
IFL权重wj计算和筛选:
IFL的效用价值越高,则对子模块的重要性就越大,该IFL的权重也就越大。当wj<5%,表明IFL的效用价值太小,删除此IFL,只保留wj≥5%的IFL。筛选留下SF,TB,ND,DD,PR,TW,EC,SR和WE共8个重要IFL,得到W1=(0.203,0.766),W2=(0.073,0.636,0.254),W3=(0.966),W4=(0.701,0.299)。
子模块Mi:
5.对fL与子模块进行多元非线性回归分析,建立L耦合模型:
其中fL生命损失率,公式中子模块的函数形式F(M)需要通过子模块计算以及多元非线性回归分析进一步确定。本发明采用留一交叉分析验证法,对fL与子模块Mi结果进行率定验证。
因为fL定义为:
fL=L/PR
其中L为生命损失,指受到溃坝洪水淹没而遇害的死亡人数;PR为风险人口,指下泄洪水或溃坝洪水淹没范围内的人数,则可计算L:
本发明的计算方法依次含有以下步骤:
输入待计算的溃坝水库生命损失影响因素信息:SF,TB,ND,DD,PR,TW,EC,SR和WE;
代入生命损失子模块得到子模块计算结果。
代入耦合模型计算,得出相应的生命损失率fL和生命损失L。
本方法建立步骤如下:1、分层抽样优选水库溃坝事件;2、生命损失影响因素IFL选取;3、生命损失影响因素分类;4、建立致灾因子、孕灾环境、承灾体和救灾能力4个生命损失子模块,并筛选出重要的生命损失影响因素;5、多元非线性回归耦合子模块,形成溃坝洪水致生命损失的计算方法。方法计算步骤如下:1、输入待计算的溃坝水库生命损失影响因素信息;2、生命损失子模块计算;3、耦合模型计算,得出相应的生命损失率fL和生命损失L。本发明为我国溃坝洪水导致的生命损失结果计算,生命损失原因分析以及生命损失风险预测提供了可靠的数据和方法支撑,为灾后重建以及灾前防范预测提供了科学的计算方法。
本实施例计算过程如下:
1.以已溃决横江水库为例进行溃坝生命损失计算和原因分析,水库于1970年9月15日溃坝,输入待计算的溃坝水库生命损失影响因素信息:SF,TB,ND,DD,PR,TW,EC,SR和WE,通过子模块计算得到4个子模块结果分别为0.784,0.673,0.756,0.431,表明横江水库的救灾能力较好,但致灾因子、孕灾环境以及承灾体条件相对较差。其中致灾因子库容较大,溃坝瞬时泄水量多,建筑物质量差是加重生命损失的主要原因。通过耦合模型计算得到生命损失率fL和生命损失L分别为0.023和1155人,实际的生命损失率fL0和生命损失L0为0.019和941人,生命损失率fL和生命损失L的相对误差分别为21.1%和22.7%,利用本方法计算的结果与真实值接近。生命损失计算结果分析表明,应重点降低致灾因子风险性,并增强建筑物抗洪能力,完善灾后建设。
2.以规划建造的云山水库为例进行溃坝生命损失风险预测,通过子模块计算得到4个子模块结果分别为0.274,0.337,0.672,0.900,子模块计算结果显示云山水库溃坝致灾因子、孕灾环境相对较好,而承灾体情况不佳,救灾能力最差。其中预警时间不充分,撤离条件差将会加重生命损失风险性。通过耦合模型计算预测生命损失率fL和生命损失L分别为0.0243和27人,生命损失结果预测表明,应重点提高水库的救灾能力,延长预警时间,提前规划完备的应急预案实施方案,加强灾害防范和宣传力度,完善撤离条件,以降低溃坝生命损失风险。
以上所述仅是本发明的优选实施方式,应当指出:对于本技术领域的普通技术人员来说,在不脱离本发明原理的前提下,还可以做出若干改进和润饰,这些改进和润饰也应视为本发明的保护范围。
Claims (8)
1.一种溃坝洪水致生命损失的计算方法,其特征在于,包括:生命损失数据选择模块、生命损失影响因素模块、生命损失子模块以及耦合模型这4个递进计算模块;
所述生命损失数据选择模块以中国溃坝事故为数据材料,采用分层抽样法,均衡考虑水平区域层、垂直海拔层和时间年代层,进行三维时空抽样,并综合考虑各类溃坝条件,优选溃坝案例;
所述生命损失影响因素模块选取的生命损失影响因素包含溃坝洪水严重程度SF、溃坝发生时间TB、建筑物损毁ND、受灾区与大坝平均距离DD、风险人口PR、对溃坝理解程度UB、预警时间TW、撤离条件EC、溃坝方式MB、库容SR和溃坝时天气WE;
根据灾害系统三要素——致灾因子DCF、孕灾环境DPE及承灾体DAB,同时考虑救灾能力DRA对溃坝生命损失有着重要影响,将生命损失影响因素分为4类;
第一类:致灾因子DCF为诱发灾害的因素,包括SF、MB、和SR;
SF的计算方法为用计算断面的溃坝峰值流量Qtop与溃坝洪水引起的最大水面宽度Wmax之商作为SF的估计值:SF≈Qtop/Wmax;
从单因素角度考虑,SF越大,则生命损失L越大;SF<0.5m2/s,对生命损失影响可以忽略不计;0.5m2/s≤SF<4.6m2/s,SF为低严重性;4.6m2/s≤SF<12.0m2/s,SF为中严重性;SF≥12.0m2/s,SF为高严重性;
MB表示溃坝方式,从单因素角度考虑,MB的频率越高,溃坝事故可能性越大,L也越大;
SR表示库容,从单因素角度考虑,SR越大,溃坝越严重,L也越大;
第二类:孕灾环境DPE为形成灾害的外部环境,包括TB、WB、ND和DD;
TB表示溃坝发生时间,从单因素角度考虑,从白天工作期逐步进入凌晨睡眠期,人类活动逐步减少,警觉性下降,L的危险性也逐渐增加;
WE表示溃坝时天气,从单因素角度考虑,WE越恶劣,L取值越大;
ND表示建筑物损毁,从单因素角度考虑,ND越多,对人类威胁性越大,L越多;
DD表示受灾区与大坝平均距离,从单因素角度考虑,DD越近,L越多;
第三类:承灾体DAB为受到损害的对象,包括PR和UB;
PR表示风险人口,为下泄洪水或溃坝洪水淹没范围内的人数,PR作为L的基数,对L影响巨大,从单因素角度考虑,PR越多,L相应也越多;
UB表示对溃坝理解程度,UB受当地文化教育程度以及当地溃坝洪水灾害宣传情况影响,UB分为清楚和模糊两类,从单因素角度考虑,UB越模糊,L越多;
第四类:救灾能力DRA为受灾对象的自救或被救的条件,包括TW和EC;
TW表示预警时间,从单因素角度考虑,TW越短,L相应也越多;
EC表示撤离条件,EC受到地形环境的影响,对受灾人员救援疏散起着重要作用,将EC分为好、中、差三类,从单因素角度考虑,EC越差,L越多;
根据第一类~第四类建立致灾因子子模块M1,孕灾环境子模块M2,承灾体子模块M3,救灾能力子模块M4;结合熵权法建立所述生命损失子模块,同时分析生命损失影响因素IFL对所述生命损失子模块的贡献率,从而筛选重要的IFL:
溃坝洪水致生命损失的计算方法的步骤如下:
S01,建立判断矩阵:
式(1)中,xij为IFL初始值,有m个溃坝事件,n个IFL,n在M1~M4中分别取3、4、2、2;
S02,M1~M4子模块标准化处理:
rij=xij/max{xij} (2)
式(2)中,rij代表xij的标准化值;
S03,根据熵的定义,计算M1~M4子模块的信息熵ej:
式(3)中,k=1/lnm,熵值反映了信息的无序化程度,可以用来度量信息量的大小,式(3)中IFL比重为pij:
S04,IFL权重wj计算和筛选:
IFL的效用价值越高,则对M1~M4子模块的重要性就越大,该IFL的权重也就越大;当wj<5%,表明IFL的效用价值太小,删除此IFL,只保留wj≥5%的IFL;
筛选留下SF,TB,ND,DD,PR,TW,EC,SR和WE共8个重要IFL,得到W1=(0.203,0.766),W2=(0.073,0.636,0.254),W3=(0.966),W4=(0.701,0.299);
S05,子模块Mi:
式(6)中Mi代表M1~M4子模块中的任一个;
S06,对fL与子模块Mi进行多元非线性回归分析,建立L耦合模型:
式(7)中,fL代表生命损失率;M1代表致灾因子子模块;M2代表孕灾环境子模块;M3代表承灾体子模块;M4代表救灾能力子模块;
fL=L/PR (8)
式(8)中,L为生命损失,指受到溃坝洪水淹没而遇害的死亡人数;
由式(7)、式(8)可推导出L为:
2.根据权利要求1所述的一种溃坝洪水致生命损失的计算方法,其特征在于:所述水平区域层包括华东地区、华南地区、华中地区、华北地区、西北地区、西南地区和东北地区;所述垂直海拔层包括库面海拔最低至30m,最高达3100m;所述时间年代层包括从20世纪50年代至21世纪10年代,每个年代选取2座溃坝水库。
3.根据权利要求1所述的一种溃坝洪水致生命损失的计算方法,其特征在于:所述MB包括4种,按照溃坝频率从高到低分别为漫坝、坝体质量差、管理不当和其他类;所述坝体质量差包括坝体、输水洞和溢洪道由于质量问题而发生渗漏、堵塞;所述管理不当包括超蓄、维护运用不良和无人管理;所述其他类包括白蚁建巢、临危扒口和原因不详。
4.根据权利要求1所述的一种溃坝洪水致生命损失的计算方法,其特征在于:所述SR包括大一型水库、大二型水库、中型水库、小一型水库和小二型水库;所述大一型水库总库容大于10亿立方米;所述大二型水库总库容在1亿立方米以上而小于10亿立方米;所述中型水库总库容大于或等于1000万立方米以上而小于1亿立方米;所述小一型水库总库容大于或等于100万立方米而小于1000万立方米;所述小二型水库总库容在10万立方米以上而小于100万立方米。
5.根据权利要求1所述的一种溃坝洪水致生命损失的计算方法,其特征在于:所述TB分为3个时段,白天工作期6:00-18:00、夜晚休息期18:00-22:00;和凌晨睡眠期22:00-次日6:00。
6.根据权利要求1所述的一种溃坝洪水致生命损失的计算方法,其特征在于:所述WB分为5个等级:一级:暴雨、暴雪、台风、雾霾;二级:大雨、大雪、大风;三级:中雨、中雪;四级:小雨、阵雨、小雪;五级:晴天、多云及阴天。
7.根据权利要求1所述的一种溃坝洪水致生命损失的计算方法,其特征在于:所述TW分为5段;0-15min、15-30min、30-45min、45-60min和大于60min。
8.根据权利要求1所述的一种溃坝洪水致生命损失的计算方法,其特征在于:采用留一交叉分析验证法,对fL与子模块Mi结果进行率定验证,得到式(7)。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610047890.0A CN105740607B (zh) | 2016-01-25 | 2016-01-25 | 一种溃坝洪水致生命损失的计算方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610047890.0A CN105740607B (zh) | 2016-01-25 | 2016-01-25 | 一种溃坝洪水致生命损失的计算方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105740607A CN105740607A (zh) | 2016-07-06 |
CN105740607B true CN105740607B (zh) | 2018-07-27 |
Family
ID=56247578
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610047890.0A Active CN105740607B (zh) | 2016-01-25 | 2016-01-25 | 一种溃坝洪水致生命损失的计算方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105740607B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108038273B (zh) * | 2017-11-24 | 2021-12-21 | 西安理工大学 | 一种瞬间溃坝坝址流量过程线的快速预测方法 |
CN110909413A (zh) * | 2019-12-09 | 2020-03-24 | 山东大学 | 一种混凝土面板堆石坝水平位移预测精确方法及系统 |
CN115034543B (zh) * | 2022-04-14 | 2024-01-05 | 水利部交通运输部国家能源局南京水利科学研究院 | 一种溃坝生命损失动态估算方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3886751A (en) * | 1973-11-12 | 1975-06-03 | Jimenez Labora Mauricio Porraz | Aquatic construction module and method of forming thereof |
US5059065A (en) * | 1991-01-25 | 1991-10-22 | David Doolaege | Apparatus and a method for joining water structure sections or the like |
CN101520324A (zh) * | 2009-03-24 | 2009-09-02 | 中国安全生产科学研究院 | 一种尾矿库溃坝灾害位移定量预警方法 |
CN104991999A (zh) * | 2015-06-17 | 2015-10-21 | 大连理工大学 | 一种基于二维sph的溃坝洪水演进模拟方法 |
-
2016
- 2016-01-25 CN CN201610047890.0A patent/CN105740607B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US3886751A (en) * | 1973-11-12 | 1975-06-03 | Jimenez Labora Mauricio Porraz | Aquatic construction module and method of forming thereof |
US5059065A (en) * | 1991-01-25 | 1991-10-22 | David Doolaege | Apparatus and a method for joining water structure sections or the like |
CN101520324A (zh) * | 2009-03-24 | 2009-09-02 | 中国安全生产科学研究院 | 一种尾矿库溃坝灾害位移定量预警方法 |
CN104991999A (zh) * | 2015-06-17 | 2015-10-21 | 大连理工大学 | 一种基于二维sph的溃坝洪水演进模拟方法 |
Non-Patent Citations (4)
Title |
---|
《A procedure for estimating loss of life caused by dam failure》;GRAHAM W J.;《Sedimentation & River Hydraulics》;19991231;第6卷(第5期);第1-43页; * |
《基于支持向量机的溃坝生命损失》;王志军等;;《水力发电》;20050131;第34卷(第1期);第67-70页; * |
《基于灰色关联分析法的溃坝》;侯保灯;《水力发电》;20121031;第38卷(第10期);第76-80页; * |
《溃坝生命损失分析方法研究》;周克发;《中国优秀博硕士学位论文全文数据库 (硕士)工程科技Ⅱ辑》;20070115(第1期);第C037-11页; * |
Also Published As
Publication number | Publication date |
---|---|
CN105740607A (zh) | 2016-07-06 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104898183B (zh) | 暴雨城市内涝建模评估方法 | |
CN104408900B (zh) | 基于动态优化的神经网络洪水预警装置及其方法 | |
Odufuwa Bashir et al. | Floods of fury in Nigerian cities | |
KR101315121B1 (ko) | 지아이에스 기반의 홍수방어 시스템 및 그 방법 | |
CN109272189A (zh) | 一种基于链式结构的城市尺度洪涝灾害风险评估方法 | |
CN105740607B (zh) | 一种溃坝洪水致生命损失的计算方法 | |
Bashir et al. | Floods of fury in Nigerian cities | |
CN111178712B (zh) | 一种用于引调水工程事故后援救预警方法及系统 | |
CN113935609B (zh) | 一种长河段多溃口凌汛洪水淹没联合风险评估方法 | |
Arnalds et al. | Avalanche hazard zoning in Iceland based on individual risk | |
Karamouz et al. | Uncertainty-based flood resiliency evaluation of wastewater treatment plants | |
CN114492233B (zh) | 一种基于webGIS平台及考虑综合利用要求的流域水模拟方法 | |
CN115809800A (zh) | 洪涝灾害风险评估方法 | |
Adhikari et al. | Achieving residential coastal communities resilient to tropical cyclones and climate change | |
CN117540931B (zh) | 一种自然灾害动态综合风险研判及预警方法 | |
CN105590034A (zh) | 防汛减灾能力评价模型 | |
CN109784720A (zh) | 一种台风灾害下基于时空网格关联的配电网风险评估方法 | |
Fang et al. | AN INTEGRATED ONE‐DIMENSIONAL AND TWO‐DIMENSIONAL URBAN STORMWATER FLOOD SIMULATION MODEL 1 | |
CN115293241A (zh) | 基于多源数据融合的河道崩岸预警方法及装置 | |
Yaoxian | Chinese experience with post-natural-disaster reconstruction | |
Cai et al. | Management of flood disasters in the Jianghan Plain, China | |
Balsells et al. | Flood resilience assessment of New Orleans neighborhood over time | |
Zhang et al. | Mortality during Hurricane Sandy: the effects of waterfront flood protection on Staten Island, New York | |
Sterlacchini | Vulnerability Assessment: concepts, definitions and methods | |
Stoyanova et al. | Flood hazard mapping using two-dimensional hydraulic modeling results |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |