基于GA优化BP神经网络的城市暴雨灾害风险评估方法及系统

技术领域

本发明涉及城市暴雨灾害风险评估的技术领域，尤其是涉及一种基于GA优化BP神经网络的城市暴雨灾害风险评估方法及系统。

背景技术

当前，我国每年由于自然灾害受到的人员伤亡和经济损失数目巨大，其中暴雨引发的灾害及衍生灾害在城市中十分普遍，因此迫切需要构建一个针对城市区域的暴雨灾害风险评估模型，对潜在可能发生的风险进行评估，既可以为城市土地资源利用提供参考，也对城市防灾减灾发挥着重要作用。

暴雨洪涝灾害风险评估是对其风险发生的强度和形式进行定量评定和估计。要进行风险评估，首先必须存在风险源，即存在自然灾变；第二，必须有风险承载体(承灾体)，即人类社会，自然灾害是自然力作用于承灾体的结果。因此，暴雨洪涝风险评估实际上是评估暴雨洪涝对承灾体的负面影响。

关于暴雨洪涝灾害风险评估，国内外学者做了大量研究，认为灾害的形成是承载体脆弱性、致灾因子和暴露度等方面综合作用的结果。暴雨洪涝灾害风险评估应综合考虑致灾因子、承灾体和防灾能力等因素，构建评估模型开展风险评估，其主要包括暴雨灾害危险性、承灾体暴露性、承灾体脆弱性以及综合风险分析，风险等级划分及其风险应对措施等内容。

暴雨洪涝灾害所造成的损失是由众多灾害影响因素相互作用的结果，而这些影响因素无法用精确的数学模型来描述。目前所采用的评估方法本身存在诸多不足，如评价机理不够直观、部分评价方法由于其建立的数学基础本身的原因，影响评估结果的可信度等。同时，防范暴雨洪涝灾害迫切需要实时动态的风险评估，而当前灾害风险实时评估的工作相对匮乏。

发明内容

本发明的目的是提供一种基于GA优化BP神经网络的城市暴雨灾害风险评估方法及系统，其具有提高评估全面性和准确性的效果。

本发明的上述发明目的一是通过以下技术方案得以实现的：

一种基于GA优化BP神经网络的城市暴雨灾害风险评估方法，包括以下步骤：

建立包括致灾因子危险性、孕灾环境敏感性、承灾体的脆弱性、防灾抗灾能力的暴雨灾害风险评估体系；

基于k-means聚类历史灾损数据生成风险等级标签，所述风险等级标签包括有极高风险、高风险、中风险和低风险四个标签；

根据所述暴雨灾害风险评估体系和风险等级标签，构建GA优化神经网络暴雨灾害风险评估模型；

将实时降雨量输入所述暴雨灾害风险评估模型中，获得特定时段内的风险等级标签。

本发明进一步设置为：所述致灾因子危险性由等效日降雨量衡量，所述等效日降雨量在某一时段的转换关系为：

D＝εPT

其中，D为等效日降雨量，P为T时长的降雨量，ε和β为转化系数；

所述致灾因子的危险性的量化计算公式为：

其中，W为暴雨致灾因子危险性，p′为累计降雨量阈值，当等效日降雨量D小于p′时，致灾因子危险性为0，当等效日降雨量超过p′时，致灾因子危险性呈指数型增长。

本发明进一步设置为：所述孕灾环境敏感性包括地形、地势、水系的评价指标，所述地形从高程数据获取，所述地势通过某一格点的高程与周围八个格点的高程标准差来计算，所述水系通过城市中河流面积和长度表示；

所述承灾体的脆弱性包括人口数量、人均GDP、路网密度三个评价指标；

所述防灾抗灾能力包括防汛路段长度和防汛点数量的两个评价指标；

将所述地形、地势、水系、人口数量、人均GDP、路网密度、防汛路段长度和防汛点数量八个评价指标无量纲化处理。

本发明进一步设置为：所述基于k-means聚类历史灾损数据生成风险等级标签，包括以下步骤：

获取历史受灾数据，统计受灾人数、直接经济损失和农作物受灾面积；

将所述受灾人数、直接经济损失和农作物受灾面积采用min-max标准化方法进行无量纲处理；

随机选择四个聚类的质心，对每个数据点求其到质心的距离，将每个数据点的类划分为离该数据点最近的质心，形成四个簇，根据分类好的簇，在每个簇内重新计算质心，重复迭代更新质心，直到达到稳定状态；

将分类后的数据，根据受灾情况分别对应到特大灾、大灾、中灾和小灾四个受灾等级，将受灾等级映射为对应的风险等级标签，特大灾对应为极高风险，大灾对应为高风险，中灾对应中风险，小灾对应低风险。

本发明进一步设置为：所述构建GA优化神经网络暴雨灾害风险评估模型，包括以下步骤：

构建包括输入层、隐含层和输出层的BP神经网络；

获取样本数据，将样本数据分为训练数据和测试数据；

使用遗传算法优化BP神经网络的权重和阈值。

本发明进一步设置为：所述输入层包括有致灾因子危险性、地形、地势、水系、人口数量、人均GDP、路网密度、防汛路段长度和防汛点数量9个神经元节点；

所述隐含层为1层，所述隐含层包括有19个神经元节点；

所述输出层包括有特大灾、大灾、中灾和小灾4个神经元节点。

本发明进一步设置为：所述使用遗传算法优化BP神经网络的权重和阈值，包括以下步骤：

步骤一：种群初始化，采用二进制编码方式对个体进行编码，连接所有所述权值和阈值的编码即为一个个体编码，当所述权值和阈值均为n位二进制编码时，所述个体的二进制编码包括有171*n位输入层与隐含层连接权值编码、19*n位隐含层阈值编码、76*n位隐含层与输出层连接权值编码、4*n位输出层编码，所述个体的二进制编码供270*n位；

步骤二：计算个体适应度，包括以下公式：

其中，minE(n)为目标函数，Fit(n)为适应度函数，N为训练集中样本个数，y

步骤三：选择运算；

步骤四：交叉运算；

步骤五：变异运算；

步骤六：当均方差E小于预设值ε或达到最大迭代次数时，则输出最佳个体，解码获得最佳神经网络权值和阈值，否则，解码得到权值和阈值，将权值和阈值赋值给新建的BP神经网络，使用训练样本训练网络，使用测试样本测试网络，再计算测试误差，返回步骤二继续执行。

本发明的上述发明目的二是通过以下技术方案得以实现的：

一种基于GA优化BP神经网络的城市暴雨灾害风险评估系统，包括：

暴雨灾害风险评估体系建立模块，用于建立包括致灾因子危险性、孕灾环境敏感性、承灾体的脆弱性、防灾抗灾能力的暴雨灾害风险评估体系；

风险等级标签生成模块，用于基于k-means聚类历史灾损数据生成包括有极高风险、高风险、中风险和低风险四个标签的风险等级标签；

暴雨灾害风险评估模型建立模块，用于根据所述暴雨灾害风险评估体系和风险等级标签，构建GA优化神经网络暴雨灾害风险评估模型；

结果输出模块，用于将实时降雨量输入所述暴雨灾害风险评估模型中，获得特定时段内的风险等级标签。

综上所述，本发明的有益技术效果为：

基于GA优化BP神经网络对暴雨灾害风险进行评估，相比于依赖层次分析法、专家打分法确定各暴雨风险评估前的权重，使用神经网络模型训练评估指标和风险等级的关系更具有科学性。该方法不依赖于主观赋值，有着很强大的非线性映射能力、自学习能力和自适应性，可以有效的从暴雨灾害历史灾损数据训练得到更加高效精准的评估指标和风险等级之间的关系模型，进而得到实时暴雨灾害风险评估产品。此外，GA优化的神经网络模型克服了BP神经网络局部极小化的不足，对BP神经网络的权值和阈值进行优化。

附图说明

图1是本发明实施例一的整体流程图；

图2是本发明实施例一中BP神经网络的拓扑结构图。

具体实施方式

实施例一

本发明公开了一种基于GA(Genetic Algorithm，遗传算法)优化BP(backpropagation，反向传播)神经网络的城市暴雨灾害风险评估方法，参照图1，包括以下步骤：

S1、建立包括致灾因子危险性、孕灾环境敏感性、承灾体的脆弱性、防灾抗灾能力的暴雨灾害风险评估体系。

建立一个完备的暴雨灾害风险评估指标体系对预测和评估灾害风险的发生至关重要，如果指标构建不全面，会导致使评估结果与实际产生较大偏差。如果将自然灾害风险从系统角度去理解，它的构成首先应该包括风险源。风险源不但从根本上决定某种自然灾害是否存在，而且还决定了该种风险的大小。我们通常采用致灾因子危险性来表述。在暴雨风险灾害系统中，将降雨强度作为致灾因子，采用模型转化到同一维度来量化致灾因子的危险性。但是风险源存在不一定意味着风险存在，风险是相对于承受灾害的行为主体的，在城市区域中，通常将人口和地均DGP作为承灾主体。此外由于交通对一个城市的人口流动和经济发展也起着至关重要的作用，故把路网密度也作为一项重要的承灾体。此外风险源的变异强度与孕灾环境有关，例如，同样的降水强度，地势低洼的地方更容易出现洪涝灾害。拟选择地势、地形和水系指标来分析孕灾环境。最后，在人类能力范围的防汛措施可大大降低暴雨灾害风险，在风险评估时亦需要将当前防灾减灾措施作为指标考虑进来。

S2、基于k-means聚类历史灾损数据生成风险等级标签，风险等级标签包括有极高风险、高风险、中风险和低风险四个标签。

S3、根据暴雨灾害风险评估体系和风险等级标签，构建GA优化神经网络暴雨灾害风险评估模型。城市暴雨灾害风险评估指标与灾害风险等级之间具有高度不确定的非线性关系，采用主观赋权法，缺乏一定的科学性，受制于专家的主观因素，故采用BP神经网络模型，将输入的风险评估指标与输出的风险级别进行非线性映射。同时由于BP神经网络存在初始值敏感，易陷于局部最优等缺点，故采用GA算法对权值和阈值进行优化，构建GA优化的BP神经网络模型。

S4、将实时降雨量输入暴雨灾害风险评估模型中，获得特定时段内的风险等级标签。依据BP神经网络模型，即可通过输入的灾害风险评估系统的各个指标，得到相应的可能发生的灾害的风险等级。在孕灾环境、承灾体及防灾抗灾能力不变的情况下，实时输入降雨量，转化为等效日降雨量并量化为致灾因子危险度后，即可对应输出该时段该地区的暴雨灾害风险等级。

在步骤S1中，致灾因子危险性由等效日降雨量衡量，等效日降雨量在某一时段的转换关系为：

D＝εPT

其中，D为等效日降雨量，P为T时长的降雨量，ε和β为转化系数；

致灾因子的危险性的量化计算公式为：

其中，W为暴雨致灾因子危险性，p′为累计降雨量阈值，当等效日降雨量D小于p′时，致灾因子危险性为0，当等效日降雨量超过p′时，致灾因子危险性呈指数型增长。

孕灾环境敏感性在本实施例中主要考虑地形、地势、水系三个方面，地形从高程数据获取，地势通过某一格点的高程与周围八个格点的高程标准差来计算，水系通过城市中河流面积和长度来量化。

承灾体的脆弱性是指当一旦发生灾变，对灾害承受主体可能遭到的损失的度量，在城市中的暴雨灾害承载体主要考虑人口数量、人均GDP、路网密度三个评价指标。

防灾抗灾能力包括防汛路段长度和防汛点数量的两个评价指标。

将地形、地势、水系、人口数量、人均GDP、路网密度、防汛路段长度和防汛点数量八个评价指标无量纲化处理。

在步骤S2中，基于k-means聚类历史灾损数据生成风险等级标签，包括以下步骤：

S21、获取历史受灾数据，统计受灾人数、直接经济损失和农作物受灾面积。

S22、将受灾人数、直接经济损失和农作物受灾面积采用min-max标准化方法进行无量纲处理。

S23、随机选择四个聚类的质心，对每个数据点求其到质心的距离，将每个数据点的类划分为离该数据点最近的质心，形成四个簇，根据分类好的簇，在每个簇内重新计算质心，重复迭代更新质心，直到达到稳定状态。

S24、将分类后的数据，根据受灾情况分别对应到特大灾、大灾、中灾和小灾四个受灾等级，将受灾等级映射为对应的风险等级标签，特大灾对应为极高风险，大灾对应为高风险，中灾对应中风险，小灾对应低风险。

步骤S3中，构建GA优化神经网络暴雨灾害风险评估模型，包括以下步骤：

S31、构建包括输入层、隐含层和输出层的BP神经网络。其中，根据暴雨灾害风险评估体系可知，从致灾因子、孕灾环境、承灾体和防灾抗灾能力四个方面分别选取了1、3、3、2个指标。故本文将降雨强度、地形、地势、水系、人口、地均GDP、路网密度、防汛点数量、防汛路段长度作为BP网络的输入层，因此输入层神经元节点数目为9个。

由于过多隐含层层数会导致训练时间增长，计算量巨大，同时也会导致在训练过程中更容易陷入局部最优，使得局部最小误差增加，基于上述考虑，本模型隐含层层数选择1层。隐含层神经元数目依据经验公式来确定，隐含层神经元个数n

输出层的节点选择由评估结果来确定。在模型的训练阶段，期望输出值是从历史灾损数据中获得，并根据上文提到的分类方法将历史灾情数据定级为特大灾、大灾、中灾和小灾四个等级，对应的输出应分别是(1,0,0，0)、(0,1,0,0)、(0，0,1,0)、(0,0,0,1),显然输出层的神经元节点数为4。

S32、获取样本数据，将样本数据分为训练数据和测试数据。训练数据用于网络训练，通过不断修正权值和阈值，来降低测试样本的测试误差。

S33、使用遗传算法优化BP神经网络的权重和阈值。优化步骤如下：

S331、种群初始化，参照图2，采用二进制编码方式对个体进行编码，连接所有权值和阈值的编码即为一个个体编码，由于该模型中输入层神经元个数为9个，隐含层神经元为19个，输出层神经元为4个，所以共有权值9*19+19*4＝247个，阈值19+4＝23个。当权值和阈值均为10位二进制编码时，个体的二进制编码包括有9*19*n位输入层与隐含层连接权值编码、19*n位隐含层阈值编码、19*4*n位隐含层与输出层连接权值编码、4*n位输出层编码，个体的二进制编码供270*n位。若每个权值和阈值使用10位的二进制编码，则每个个体的二进制编码长度为2700，其中前1710位是输入层与隐含层连接权值编码，1711-1900位是隐含层阈值编码，1901-2660位是隐含层与输出层连接权值编码，2661-2700位则为输出层编码。

S332、为了使BP网络预测时，预测值与期望值的均值误差尽可能小，选取网络在训练样本中产生的均值误差作为学习的目标函数，并利用目标函数的倒数作为适应度函数。计算个体适应度，包括以下公式：

其中，minE(n)为目标函数，Fit(n)为适应度函数，N为训练集中样本个数，y

S333、选择运算，采用随机遍历抽样，为现有技术。

S334、交叉运算，采用最简单的单点交叉算子，为现有技术。

S335、变异运算，变异以一定概率产生变异基因数，用随机方法选出发生变异的基因。如果所选的基因的编码为1，则将其变为0；反之则变为1。

S336、当均方差E小于预设值ε或达到最大迭代次数时，则输出最佳个体，解码获得最佳神经网络权值和阈值，否则，解码得到权值和阈值，将权值和阈值赋值给新建的BP神经网络，使用训练样本训练网络，使用测试样本测试网络，再计算测试误差，返回步骤二继续执行。

实施例二

本发明公开了一种基于GA优化BP神经网络的城市暴雨灾害风险评估系统，包括：

暴雨灾害风险评估体系建立模块，用于建立包括致灾因子危险性、孕灾环境敏感性、承灾体的脆弱性、防灾抗灾能力的暴雨灾害风险评估体系；

风险等级标签生成模块，用于基于k-means聚类历史灾损数据生成包括有极高风险、高风险、中风险和低风险四个标签的风险等级标签；

暴雨灾害风险评估模型建立模块，用于根据暴雨灾害风险评估体系和风险等级标签，构建GA优化神经网络暴雨灾害风险评估模型；

结果输出模块，用于将实时降雨量输入暴雨灾害风险评估模型中，获得特定时段内的风险等级标签。

本具体实施方式的实施例均为本发明的较佳实施例，并非依此限制本发明的保护范围，故：凡依本发明的结构、形状、原理所做的等效变化，均应涵盖于本发明的保护范围之内。