CN108280553B - 基于gis-神经网络集成的山洪灾害风险区划及预测方法 - Google Patents
基于gis-神经网络集成的山洪灾害风险区划及预测方法 Download PDFInfo
- Publication number
- CN108280553B CN108280553B CN201810157583.7A CN201810157583A CN108280553B CN 108280553 B CN108280553 B CN 108280553B CN 201810157583 A CN201810157583 A CN 201810157583A CN 108280553 B CN108280553 B CN 108280553B
- Authority
- CN
- China
- Prior art keywords
- risk
- disaster
- layer
- mountain torrent
- index
- 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
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/02—Neural networks
- G06N3/08—Learning methods
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Strategic Management (AREA)
- Economics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Game Theory and Decision Science (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Tourism & Hospitality (AREA)
- General Business, Economics & Management (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Computation (AREA)
- Artificial Intelligence (AREA)
- Biomedical Technology (AREA)
- Biophysics (AREA)
- Computational Linguistics (AREA)
- Data Mining & Analysis (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明涉及一种基于GIS‑神经网络集成的山洪灾害风险区划及预测方法,包括:S1.利用关联规则挖掘山洪灾害中风险因子与风险等级之间的关联关系,辨识风险因子,构建定量化的山洪灾害风险评价指标体系;S2.采用层次分析法确定危险性和易损性指标体系及其权重,生成各要素图层;S3.利用ArcGIS将山洪灾害危险性和易损性分布图层叠加得到山洪灾害风险分布图;S4.采用ISO最大似然法聚类及自下而上区域合并与自上而下定性分析结合的方法,形成山洪灾害风险区划;S5.利用Elman神经网络分析评价指标与风险等级、灾情数据之间的非线性关系,构建山洪灾害风险评价及损失预估模型。本发明解决了变化环境下山洪灾害评估中的空间尺度不确定性问题,可广泛用于山洪灾害风险评估。
Description
技术领域
本发明涉及山洪灾害防治领域,更具体地,涉及一种基于GIS-神经网络集成的山洪灾害风险区划及预测方法。
背景技术
我国是一个多山的国家,山丘区面积约占全国陆地面积的2/3。复杂的地形地质条件、暴雨多发的气候特征、密集的人口分布和人类活动的影响,导致山洪灾害发生频繁。我国山洪灾害呈现出影响范围广、发生频率高、灾害历时短、危害破坏程度大、地域差异显著的整体特征。据《全国山洪灾害防治规划报告》数据统计,我国山丘区流域面积在100km2以上的山溪河流约5万条,其中70%因受降雨、地形及人类活动影响会发生山洪灾害。据官方统计数据表明,我国每年因山洪灾害直接或间接导致人民群众生命财产的损失约占自然灾害造成的总损失的40%左右。目前,我国山洪灾害的研究已经涉及成因、空间分布特征、灾害损失评估、风险评价与制图等各方面。风险评估与管理逐渐也成为国际上倡导和推广的减灾防灾有效途径之一。目前,山洪灾情评估工作得到了来自地学工作者、工程专家和各级政府部门的高度重视,并逐渐成为国际性的研究项目。特别是在山洪风险评估方面的表现尤为突出。
变化环境下山洪灾害的形成条件及灾情分布具有空间尺度上的非线性特点,无论是在风险评价指标体系,还是在风险评估数学模型中,都存在不同空间尺度范围的不确定性,这对准确评估山洪灾害风险带来了很大的不确定性。应用时间序列分析、多元回归分析、延神经元网络、人工神经网络、模糊数学方法、混沌优化算法、灰色聚类分析等数理统计方法、地理信息系统的空间分析方法等方法构建山洪灾害评估体系,进行山洪灾害定量评估具有很大前景性。
发明内容
本发明的目的在于克服现有技术的不足,提供一种基于GIS-神经网络集成的山洪灾害风险区划及预测方法,构建了山洪灾害风险评价及损失预估模型,解决了变化环境下山洪灾害评估中的空间尺度不确定性问题。
为解决上述问题,本发明提供的技术方案为:一种基于GIS-神经网络集成的山洪灾害风险区划及预测方法,其中,包括以下步骤:
S1.利用关联规则挖掘山洪灾害中风险因子与风险等级之间的关联关系,辨识风险因子,构建定量化的山洪灾害风险评价指标体系;
本发明利用关联规则挖掘山洪灾害中风险因子与风险等级之间的关联关系,构建了定量化的山洪灾害风险评价指标体系。在初步构建的风险因子数据库,用项集A来表示,如式(1)所示,基于历史灾情分析的损失风险等级用项集B来表示,如式(2)所示,由此构建用于山洪灾害的“风险因子-风险等级”不确定性分析的事物集。
A={A1,A2,A3,...,Am} (1)
B={B1,B2,B3,...,Bn} (2)
其中,A为山洪灾害的环境条件因素项集,A1,A2,A3,…,Am分别对应于某个山洪灾害案例的环境条件因素状态的组合;B为山洪灾害的风险等级项集,B1,B2,B3,…,Bn等表征了山洪灾害的风险特征。
关联规则可以用X=>Y来表示,其中X、Y是项集,I是所有项集的总和,X∈I,Y∈I,X∩Y=Φ,X=>Y表示在数据库的事件中,包含项集X的事件同时也包含项集Y的可能性。然后,在遍历所有项集的基础上,利用式(3)和式(4)进行关联规则分析,挖掘不同灾害环境组合条件项与不同灾害等级项之间的强关联规则,如图1所示。进一步地,探究其间的物理成因,与关联规则的分析结果进行验证。进而对各风险因子在山洪灾害风险评价中的作用大小及其不确定性进行定量分析,构建相应的风险评价指标体系,为研究更加符合山洪灾害成因机理和客观规律的风险评价模型提供科学合理的物理基础。
其中,式(3)为支持度的计算式,表示在数据库中项集X和项集Y同时出现的概率,DX∪Y为包含X且包含Y的元组数,D表示总元组数;式(4)为置信度的计算式,表示在出现项集X的事件中,项集Y也出现的概率,DX∪Y为包含X且包含Y的元组数,DX表示包含X的元组数。
S2.采用层次分析法确定危险性和易损性指标体系及其权重,生成各要素图层;
本发明利用层次分析法对山洪灾害风险评价指标体系进行定权,其主要原理是将与决策相关的元素分解成目标层(A层),准则层(C层),指标层或决策层(C层)三个层面;将各层次的指标两两进行比较,对其予以优劣评判;并利用评判结果建立矩阵模型来计算各指标的权重系数,即得到在某一准则下各个指标相对重要性的度量,具体计算步骤为:
S21.建立层次递阶结构
层次分析的结构模型大致分为三层,从最高层到最底层依次是目标层,准则层和变量层。目标层一般情况下是唯一的,研究所要达到的目标;准则层是为实现总目标所采取的措施;变量层是影响措施的因素
S22.构建判断矩阵
判断矩阵是针对上一层次中某元素而言,评定该层次中各元素间相对重要程度的判断。两两指标之间的相对重要程度的度量尺度一般采用9分法。划分标准如表1所示:
表1 AHP法程度赋值及意义
在山洪灾害发生的危险性因子中,降雨是引发山洪灾害的主导因素,其影响力最大,地形地貌、河网水系等下垫面条件又为山洪灾害的发生提供了孕灾环境,但影响力较降雨小;而在山洪灾害易损性分析的各指标中,因为人口是山洪灾害中最重要的保护对象,城镇化率是防洪能力的体现,所以人口和城镇化率指标是最重要的指标,地均GDP、一产比重反映了资产暴露状况和下垫面土壤持水能力影响防洪能力,对山洪灾害易损性评价有重大影响,由此建立各层次判断矩阵见表2。
表2各层次判断矩阵
S23.计算判断矩阵的最大特征向量与特征根;
最后,计算判断矩阵的最大特征值λmax:
S24.计算判断矩阵一致性指标,并检验其一致性;引入一致性比率CR:
式中,n为判断矩阵阶数,CI为一致性指标,CR为随机一致性比率,RI为随机一致性指标;
由于不同阶数的矩阵两两判断的比例测度达到一致性的难度不一样,用RI来对不同阶数矩阵一致性指标CI进行修正。随机一致性指标的统计平均值如表3所示。
表3平均随机一致性指标
若CR<0.1,判断矩阵具有很好一致性,判断合理;
若CR=0.1,判断矩阵具有较好一致性,判断较合理;
若CR>0.1,判断矩阵不符合一致性原则,需重新调整;
根据以上步骤,计算得到判断矩阵的最大特征向量、特征根、一致性检验值如表4所示:
表4最大特征向量、特征根、一致性检验值
S25.通过加权综合计算得到准则层C层指标的相对权重;通过步骤S21至S24的计算,得到一层指标相对于上一层指标的权重,再计算指标层各指标相对于目标层的总权重,其计算公式为:
则,指标层对目标层的总权重为:
根据以上方法,构造出通过一致性检验的判断矩阵,并确定准则层和指标层各指标的相对权重。
S3.利用ArcGIS将山洪灾害危险性和易损性分布图层叠加得到山洪灾害风险分布图;
S31.依照层次分析法建立的指标体系和权重对各指标层的指标进行汇总,得出各栅格的山洪风险值;风险值越高,意味着面临的山洪风险越大。
S32.利用ArcGIS栅格计算器通过以下公式将不同指标图层叠加起来,得到危险性、易损性和风险性栅格分布图,其中,
风险性计算公式为:D=f(H,V)=ωHH+ωVV
式中,hi、vi为各个栅格单元经过标准化处理后危险性指标和易损性指标的取值;ω为危险性和易损性各指标权重
S4.采用ISO最大似然法聚类以及自下而上区域合并与自上而下定性分析结合的方法,形成山洪灾害风险区划;
ISO最大似然法聚类:ISO聚类是以空间相似性为原则的非监督分类方法,主要运算思想是:输入N个样本,预选Nc个初始聚类中心和各参数指标:预期的聚类中心数目K、每一聚类域中最少的样本数θN、聚类域中样本距离分布的标准差θS、两聚类中心间的最小距离θc、一次迭代运算中可以合并的聚类中心的最多对数L、迭代运算的次数I;当某两类聚类中心距离θc小于阀值则归为一类,否则为不同类;若某类样本数θN小于阀值,则将其取消。ISO聚类不需要提供监督样本,自动进行合并和分裂,进而得到较为合理的聚类效果。将其与最大似然法结合,能得到与实际条件相吻合的区划分类结果。在本次区划评价中主要通过ArcGIS中ISO聚类工具和最大似然法分类工具实现。
“自上而下”分区与“自下而上”合并相结合:参考《全国山洪灾害风险区划》,以“自上而下”分区与“自下而上”合并相结合的方法为主要思路进行区域划分。区划原则应与区划方法相对应,本次评价所遵循的原则包括:主导因素与综合因素相结合原则、区域单元内部相对一致性原则、以人为本的经济社会分析原则、区域共轭性原则、定性分析与定量研究相结合原则,而“自上而下”区划方法是为相对一致性原则而设计的,“自下而上”区划方法是为区域共轭性原则而设计的。
以上两种方法都是自然灾害区划乃至自然区划中最通用的方法。“自上而下”方法进行区划时,要掌握宏观格局,根据某些区划指标,首先进行最高级别单位的划分然后依次将已划分出的高级单位再划分成低四级的单位,一直划分到最低级区划单位为止;“自下而上”方法则恰恰相反,它通过对最小图斑的指标的分析,首先合并出最低级的区划单位,然后再在低级区划单位的基础上,逐步合并出较高级别的单位,直到得出最高级别的区划单位为止。
S5.利用Elman神经网络分析评价指标与风险等级、灾情数据之间的非线性关系,构建山洪灾害风险评价及损失预估模型。
Elman神经网络是Elman于1990年提出的一种典型的局部回归网络。除了拥有一般的神经网络所具有的输入层、隐含层、输出层之外,还有一个特殊的承接层,用于构成局部反馈,可以改善BP神经网络训练时间较长,和易于陷入局部最优的状况。承接层选择线性函数为传递函数,由于多出一个延迟单元,承接层能够对过去状态进行记录,这种记录在下一次也会作为输入单元,这种动态记忆功能,让Elman神经网络在时间序列预测上有较好的表现。
在本发明中,若Elman神经网络的输入u(k-1)∈Rr,输出y(k)∈Rm,隐含层输出x(k)∈Rn,承接层输出xc(k)∈Rn,Elman网络的数学模型为:
x(k)=f(W1xc(k)+W2u(k-1))
xc(k)=αxc(k-1)+x(k-1)
y(k)=g(W3x(k))
其中,W1∈Rn×nW2∈Rn×r、W3∈Rm×n分别为承接层到隐含层、输入层到隐含层、隐含层到输出层的连接权矩阵;f(·)和g(·)分别为隐含层神经元和输出层神经元的激发函数组成的非线性向量函数;
进一步地,所述的S5步骤中利用Elman神经网络,构建山洪灾害风险评价及损失预测模型包括以下步骤:
S51.数据准备,将各评价指标原始数据作为输入数据,利用ArcGIS自然间断法将各指标数据分为5个灾害风险等级并作为风险评价Ⅰ,将基于历史灾情评价的损失风险评估结果作为损失评价集Ⅱ,将风险评价集Ⅰ和损失评价集Ⅱ作为输出数据,随机选取训练样本和验证样本,训练样本和验证样本的比例为3:1。
S52.数据预处理,采用premnmx函数对输入和输出数据进行归一化处理至[-1,1]区间内;神经网络的输入数据大多存在着多维、多纲、多源的特点,若将数据直接输入网络中,可能需要花费较长的训练时间,且当数值差别很大时,大数值往往会更多的影响网络的学习,大大降低网络学习速度更可能导致最终无法收敛。因此,有必要对输入数据进行合理的预处理。
S53.Elman神经网络构建,网络结构关系决定着神经网络解决问题的能力和效率,但在隐含层层数和节点数的选择上,当前还未存在确切的理论和方法。一般情况下,单隐含层的神经网络能够逼近任意的非线性函数,因此本发明选择单隐含层的神经网络即三层Elman神经网络来进行模拟灾害风险评价和灾害损失预估。在隐含层节点的选择上,通过训练数据不断调试来确定。隐含层传输函数选择非线性的双曲正切函数(tanh或tansig),其公式为
输出层使用线性传输函数purelin,训练函数选择带有动量项的自适应学习算法traingdx;训练函数的选择上,主要考虑以下三种方式:带有动量项的梯度下降法(traingdm)、带有动量项的自适应学习算法(traingdx)和Levenberg-Marquardt算法(trainlm)。其中,traingdm和traingdx是基于基本梯度下降法(traing)发展而来的方法,在收敛速度上,traing<traingdm<traingdx;对于中等规模的网络,trainlm是达到收敛需要训练次数最少的方法,但其在运算过程中内存占用较大。
S54.Elman神经网络训练,根据S3步骤构建的神经网络,分别对灾害风险评价和灾害损失预估利用训练输入数据进行调试,traingdx训练方法设定训练步数为10000步,期望训练精度根据不同的数据输入,选择0.008-0.05不等,通过比较模拟值和实际值的拟合程度,确定模型隐含层节点;
S55.Elman神经网络预测,将验证输入数据分别输入到训练好的灾害风险评价和灾害损失预估模型中进行预测,分别得到山洪灾害风险评价结果和山洪灾害损失预估结果,应用上述训练好的网络,将各指标图层栅格值作为输入数据,模拟灾害风险及损失,利用ArcGIS输出山洪灾害风险评价等级图和灾害损失预测分布图;
S56.误差分析与泛化验证,通过误差函数评价模拟效果,并更换训练数据和验证数据,对模型的泛化能力进行验证。
与现有技术相比,有益效果是:本发明提供的一种基于GIS-神经网络集成的山洪灾害风险区划及预测方法,利用关联规则构建了定量化的山洪灾害风险评价指标体系,减少了构建指标体系过程中主观随意性的影响;集成AHP-GIS分析方法,分析山洪灾害危险性和易损性水平,进行风险评价,得到山洪灾害风险等级分布图;采用ISO最大似然法聚类及自上而下定性分析方法,形成山洪灾害风险区划;最后利用Elman神经网络分析评价指标与风险等级、灾情数据之间的非线性关系,构建山洪灾害风险评价及损失预估模型,解决了变化环境下山洪灾害评估中的空间尺度不确定性问题,本发明可用于山洪灾害风险研究分析,为山洪灾害的空间预测和防治提供科学依据。
附图说明
图1为本发明方法整体流程图。
图2为本发明关联规则辨识山洪灾害驱动因子概化图。
图3为本发明基于神经网络的山洪灾害风险评估模型。
图4为本发明实施例中广东省山洪灾害风险等级分布图。
图5为本发明实施例中广东省山洪灾害风险分析聚类结果。
图6为本发明实施例中广东省山洪灾害风险区划。
图7为本发明实施例中选取的历史山洪灾害情况验证结果。
具体实施方式
如图1至图3所示,一种基于GIS-神经网络集成的山洪灾害风险区划及预测方法,其中,包括以下步骤:
步骤1:本发明选择广东省作为实验区域。广东省位于我国大陆南部,地形地貌类型复杂多样,且以山地丘陵为主,海拔500m以上的山地占31.7%,丘陵占28.5%。加之地处亚热带季风区,暴雨频发,受特殊自然地理环境、日趋极端的灾害性天气以及山丘区人类经济社会活动等多种自然和人为因素的影响,导致广东省山洪灾害发生频繁。因此,为进一步做好山洪灾害的防御工作,本发明提出一种山洪灾害风险区划方法,并将其用于广东省的山洪灾害风险区划,为广东省山洪灾害的空间预测和防治提供科学依据。
步骤2:首先利用层次分析法对山洪灾害风险评价指标体系进行定权,如表5、表6所示,将各层次的指标两两进行比较,对其予以优劣评判,并利用评判结果建立矩阵模型来计算各指标的权重系数,即得到在某一准则下各个指标相对重要性的度量。本发明认为广东省山洪灾害危险性受降雨因子、地形因子(地貌、坡度)、河网密度影响,其物理意义在于:区域暴雨、地形、水文水系条件是影响山洪过程分布格局最为重要的3项因子,且只有3项因子均达到较高水平,才有可能发育规模-频率水平较高的山洪过程。本发明认为反映人口集中程度的人口密度和城镇化率指标应赋予高权重,地均GDP和第一产业产值比重反映资产的集中程度,权重设置可以相对低些。
表5广东省山洪灾害危险性评价指标体系及权重
步骤3:山洪灾害风险是山洪灾害危险性、易损性的结合体,是在山洪灾害危险性、易损性等要素评估结果的基础上进行评价的,用山洪灾害综合风险指数来表征其等级。本发明利用层次分析法确定风险性评价指标权重后,将危险性、易损性的栅格图叠加后,得到广东省山洪灾害风险等级分布图,如图4所示。
步骤4:利用ArcGIS软件对广东省山洪灾害风险评价结果进行ISO最大似然法聚类分析,得到广东省山洪灾害风险分析5类聚类结果,如图5所示。
步骤5:为了保证区划单元的准确性和科学性,在单纯地采用聚类分析后,需将自下而上区域合并工作与自上而下定性分析紧密结合,从而得到最终的风险区划单元。根据全国山洪灾害风险三级区划,广东省划分为滇东南粤桂山地中高风险区和粤桂低山平原中风险区2个三级区划单元。在三级区划框架内,按照区划的集中连片原则以及区域共轭性原则,同时参阅《广东省地貌区划》、《广东省气候区划》等区划成果,进一步划分28个山洪灾害风险等级相对一致的四级区划单元,形成广东省山洪灾害风险四级区划,如图6所示。
步骤6:将本发明评价结果与历史山洪灾害分布情况进行对比,结果显示风险等级较高的区域均是历史山洪灾害发生次数与损失较大的区域,如图7所示。如清远市的清新、英德、佛冈、阳山、连山、连南;韶关市的乐昌;河源市的连平、东源、紫金;梅州市的五华、兴宁;惠州市的龙门、博罗;阳江市的阳春;茂名市的信宜、高州、电白、化州等地。根据初步对比验证结果,认为本次广东省山洪灾害风险区划结果较为符合实际情况,可为山洪灾害的空间预测和防治提供科学依据。
山洪灾害风险区划作为一项防洪减灾非工程措施,是洪水灾害风险分析系统中必不可少的组成部分。应用关联规则、人工神经网络以及地理信息系统等方法构建山洪灾害评估体系,进行山洪灾害定量评估,实施山洪灾害风险区划,是山洪灾害风险研究的主要发展方向。本发明利用关联规则构建了定量化的山洪灾害风险评价指标体系,减少了构建指标体系过程中主观随意性的影响;集成AHP-GIS分析方法,分析山洪灾害危险性和易损性水平,进行风险评价,得到山洪灾害风险等级分布图;采用ISO最大似然法聚类及自上而下定性分析方法,形成山洪灾害风险区划;最后利用Elman神经网络分析评价指标与风险等级、灾情数据之间的非线性关系,构建山洪灾害风险评价及损失预估模型,解决了变化环境下山洪灾害评估中的空间尺度不确定性问题。因此本发明可用于山洪灾害风险研究分析,并可为山洪灾害的空间预测和防治提供科学依据。
显然,本发明的上述实施例仅仅是为清楚地说明本发明所作的举例,而并非是对本发明的实施方式的限定。对于所属领域的普通技术人员来说,在上述说明的基础上还可以做出其它不同形式的变化或变动。这里无需也无法对所有的实施方式予以穷举。凡在本发明的精神和原则之内所作的任何修改、等同替换和改进等,均应包含在本发明权利要求的保护范围之内。
Claims (2)
1.一种基于GIS-神经网络集成的山洪灾害风险区划及预测方法,其特征在于,包括以下步骤:
S1.利用关联规则挖掘山洪灾害中风险因子与风险等级之间的关联关系,辨识风险因子,构建定量化的山洪灾害风险评价指标体系;
首先构建山洪灾害风险指标因子数据库,用项集A来表示,如式①所示,将山洪灾害的风险等级用项集B来表示,如式②所示,同时,考虑径流量在暴雨山洪灾害发生过程中所起表征作用,以网格化后流域的典型断面为研究对象,将河段径流量变化幅度划分阈值,用项集C来表示,如式③所示,由此构建用于山洪灾害“风险因子-径流量变化-灾害等级”不确定性分析的事物集;
A={A1,A2,A3,...,Am} ①
B={B1,B2,B3,...,Bn} ②
C={C1,C2,C3,...,Cn} ③
其中,A为山洪灾害的环境条件因素项集,A1,A2,A3,…,Am分别对应于某个山洪灾害案例的环境条件因素状态的组合;B为山洪灾害的风险等级项集,B1,B2,B3,…,Bn表征了山洪灾害的风险特征;C为典型断面径流量的变化幅度项集,C1,C2,C3,…,Cn表征了山洪灾害案例中径流量的变化特征;
关联规则用X=>Y来表示,其中X、Y是项集,I是所有项集的总和,X∈I,Y∈I,X∩Y=Φ,X=>Y表示在数据库的事件中,包含项集X的事件同时也包含项集Y的可能性;然后,在遍历所有项集的基础上,利用式④和式⑤进行关联规则分析,挖掘不同灾害环境组合条件项与不同灾害等级项之间的强关联规则;进一步地,探究其间的物理成因,与关联规则的分析结果进行验证;进而对各暴雨山洪灾害风险因子在径流量变化、暴雨山洪灾害风险评价中的作用大小及其不确定性进行定量分析,构建相应的风险评价指标体系,为研究更加符合暴雨山洪灾害成因机理和客观规律的风险评价模型提供科学合理的物理基础;
其中,式④为支持度的计算式,表示在数据库中项集X和项集Y同时出现的概率,DX∪Y为包含X且包含Y的元组数,D表示总元组数;式⑤为置信度的计算式,表示在出现项集X的事件中,项集Y也出现的概率,DX∪Y为包含X且包含Y的元组数,DX表示包含X的元组数;
S2.采用层次分析法确定危险性和易损性指标体系及其权重,生成各要素图层;所述的S2步骤具体包括:
S21.建立层次递阶结构,层次分析的结构模型分为三层,从最高层到最底层依次是目标层,准则层和变量层;
S22.构建判断矩阵;判断矩阵是针对上一层次中某元素而言,评定该层次中各元素间相对重要程度的判断;
S23.计算判断矩阵的最大特征向量与特征根;
最后,计算判断矩阵的最大特征值λmax:
S24.计算判断矩阵一致性指标,并检验其一致性;引入一致性比率CR:
式中,n为判断矩阵阶数,CI为一致性指标,CR为随机一致性比率,RI为随机一致性指标;
若CR<0.1,判断矩阵具有很好一致性,判断合理;
若CR=0.1,判断矩阵具有较好一致性,判断较合理;
若CR>0.1,判断矩阵不符合一致性原则,需重新调整
S25.通过加权综合计算得到准则层C层指标的相对权重;通过步骤S21至S24的计算,得到一层指标相对于上一层指标的权重,再计算指标层各指标相对于目标层的总权重,其计算公式为:
则,指标层对目标层的总权重为:
S3.利用ArcGIS将山洪灾害危险性和易损性分布图层叠加得到山洪灾害风险分布图;
S4.采用ISO最大似然法聚类以及自下而上区域合并与自上而下定性分析结合的方法,形成山洪灾害风险区划;
S5.利用Elman神经网络分析评价指标与风险等级、灾情数据之间的非线性关系,构建山洪灾害风险评价及损失预估模型;所述的S5步骤中Elman神经网络的输入u(k-1)∈Rr,输出y(k)∈Rm,隐含层输出x(k)∈Rn,承接层输出xc(k)∈Rn,Elman网络的数学模型为:
x(k)=f(W1xc(k)+W2u(k-1))
xc(k)=αxc(k-1)+x(k-1)
y(k)=g(W3x(k))
其中,W1∈Rn×n、W2∈Rn×r、W3∈Rm×n分别为承接层到隐含层、输入层到隐含层、隐含层到输出层的连接权矩阵;f(·)和g(·)分别为隐含层神经元和输出层神经元的激发函数组成的非线性向量函数;
所述的S5步骤中利用Elman神经网络,构建山洪灾害风险评价及损失预测模型包括以下步骤:
S51.数据准备,将各评价指标原始数据作为输入数据,利用ArcGIS自然间断法将各指标数据分为5个灾害风险等级并作为风险评价Ⅰ,将基于历史灾情评价的损失风险评估结果作为损失评价集Ⅱ,将风险评价集Ⅰ和损失评价集Ⅱ作为输出数据,随机选取训练样本和验证样本,训练样本和验证样本的比例为3:1;
S52.数据预处理,采用premnmx函数对输入和输出数据进行归一化处理至[-1,1]区间内;
S53.Elman神经网络构建,选择单隐含层的神经网络即三层Elman神经网络来进行模拟灾害风险评价和灾害损失预估,其中,
隐含层传输函数选择非线性的双曲正切函数tanh或tansig,其公式为
输出层使用线性传输函数purelin,训练函数选择带有动量项的自适应学习算法traingdx;
S54.Elman神经网络训练,根据S3步骤构建的神经网络,分别对灾害风险评价和灾害损失预估利用训练输入数据进行调试,traingdx训练方法设定训练步数为10000步,期望训练精度根据不同的数据输入,选择0.008-0.05不等,通过比较模拟值和实际值的拟合程度,确定模型隐含层节点;
S55.Elman神经网络预测,将验证输入数据分别输入到训练好的灾害风险评价和灾害损失预估模型中进行预测,分别得到山洪灾害风险评价结果和山洪灾害损失预估结果,应用上述训练好的网络,将各指标图层栅格值作为输入数据,模拟灾害风险及损失,利用ArcGIS输出山洪灾害风险评价等级图和灾害损失预测分布图;
S56.误差分析与泛化验证,通过误差函数评价模拟效果,并更换训练数据和验证数据,对模型的泛化能力进行验证。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810157583.7A CN108280553B (zh) | 2018-02-24 | 2018-02-24 | 基于gis-神经网络集成的山洪灾害风险区划及预测方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810157583.7A CN108280553B (zh) | 2018-02-24 | 2018-02-24 | 基于gis-神经网络集成的山洪灾害风险区划及预测方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108280553A CN108280553A (zh) | 2018-07-13 |
CN108280553B true CN108280553B (zh) | 2020-10-02 |
Family
ID=62808699
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810157583.7A Active CN108280553B (zh) | 2018-02-24 | 2018-02-24 | 基于gis-神经网络集成的山洪灾害风险区划及预测方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108280553B (zh) |
Families Citing this family (31)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109086887A (zh) * | 2018-08-06 | 2018-12-25 | 贵州大学 | 深度rbf神经网络与基于熵权的ahp结合的预警方法 |
CN109359738A (zh) * | 2018-10-19 | 2019-02-19 | 西南交通大学 | 一种基于qpso-bp神经网络的滑坡危险性评估方法 |
CN109740924B (zh) * | 2018-12-29 | 2022-12-06 | 西安电子科技大学 | 融合属性信息网络和矩阵分解的物品评分预测方法 |
CN109872062A (zh) * | 2019-02-02 | 2019-06-11 | 湖南大学 | 基于三维指标框架的县域尺度水土保持防控区划方法及系统 |
CN109784779B (zh) * | 2019-03-04 | 2021-06-18 | 广州杰赛科技股份有限公司 | 财务风险预测方法、装置及存储介质 |
CN109933921B (zh) * | 2019-03-18 | 2020-06-09 | 中国科学院地理科学与资源研究所 | 一种滚石灾害风险评估方法、装置、系统及存储介质 |
CN110377680B (zh) * | 2019-07-11 | 2020-04-28 | 中国水利水电科学研究院 | 基于网络爬虫和语义识别的山洪灾害数据库构建及更新的方法 |
CN110633858A (zh) * | 2019-09-18 | 2019-12-31 | 国网湖南省电力有限公司 | 一种用于输电线路的暴雨地质灾害的聚类预警方法及系统 |
CN110738355B (zh) * | 2019-09-19 | 2023-07-04 | 河源职业技术学院 | 一种基于神经网络的城市内涝预测方法 |
CN111008734A (zh) * | 2019-11-26 | 2020-04-14 | 华南理工大学 | 洪水危险性的区划和时空演变规律的研究方法 |
CN111079999A (zh) * | 2019-12-04 | 2020-04-28 | 中国地质大学(武汉) | 一种基于cnn和svm的洪水灾害易发性预测方法 |
CN111047099B (zh) * | 2019-12-16 | 2020-08-21 | 杭州鲁尔物联科技有限公司 | 一种区域性山洪风险预测方法及系统 |
CN111798644A (zh) * | 2020-07-03 | 2020-10-20 | 成都理工大学 | 一种便携式山洪灾害预警方法与系统 |
CN112036680A (zh) * | 2020-07-08 | 2020-12-04 | 贵州电网有限责任公司 | 一种基于深度神经网络的电网灾情应急演练管理系统 |
CN112001565A (zh) * | 2020-09-08 | 2020-11-27 | 清华大学合肥公共安全研究院 | 基于Softmax回归模型的地震灾害损失预测与评估的方法及系统 |
CN112308292B (zh) * | 2020-10-10 | 2023-01-20 | 广东工业大学 | 一种火灾风险等级分布图的绘制方法 |
CN112966856A (zh) * | 2021-02-10 | 2021-06-15 | 四川水利职业技术学院 | 一种山洪风险预测方法及预测系统 |
CN112966926B (zh) * | 2021-03-02 | 2022-04-22 | 河海大学 | 一种基于集成学习的洪水敏感性风险评估方法 |
CN113191582B (zh) * | 2021-03-15 | 2022-09-06 | 西南石油大学 | 一种基于gis与机器学习的道路山洪易损性评价方法 |
CN113487123B (zh) * | 2021-05-11 | 2022-05-27 | 长江水利委员会长江科学院 | 高光谱监测与gis耦合山洪灾害动态风险评估方法 |
CN113177737A (zh) * | 2021-05-26 | 2021-07-27 | 南京恩瑞特实业有限公司 | 基于ga优化bp神经网络的城市暴雨灾害风险评估方法及系统 |
CN113313384B (zh) * | 2021-05-28 | 2022-10-25 | 华南理工大学 | 一种融合弹性的城市洪涝灾害风险评估方法 |
CN113283817B (zh) * | 2021-07-20 | 2021-09-21 | 光谷技术有限公司 | 一种灾害评估方法及系统 |
CN113554345B (zh) * | 2021-08-10 | 2023-04-18 | 山西省地震局 | 地震滑坡灾害区域及灾害链风险评估方法 |
CN113869438B (zh) * | 2021-09-30 | 2023-04-07 | 自然资源部第三海洋研究所 | 一种基于k-means的风暴潮灾害重点防御区选划方法 |
CN114186780B (zh) * | 2021-11-04 | 2022-07-22 | 河海大学 | 一种基于机器学习的山洪灾害区划方法 |
CN114723218B (zh) * | 2022-02-21 | 2024-04-02 | 浙江海洋大学 | 一种基于信息量-神经网络的油气管道地质灾害评价方法 |
CN114595581A (zh) * | 2022-03-11 | 2022-06-07 | 重庆地质矿产研究院 | 基于影响因素动态权重的区域地质灾害危险性概率模型 |
CN116090839B (zh) * | 2023-04-07 | 2023-06-16 | 水利部交通运输部国家能源局南京水利科学研究院 | 水资源耦合系统多重风险分析与评估方法及系统 |
CN116433032B (zh) * | 2023-04-26 | 2024-04-09 | 中国农业科学院农业环境与可持续发展研究所 | 基于网络爬虫方式的智能评估方法 |
CN117610945B (zh) * | 2024-01-24 | 2024-04-05 | 四川省华地建设工程有限责任公司 | 基于人工智能的复杂山区地质灾害风险评价方法及系统 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101726750A (zh) * | 2009-11-30 | 2010-06-09 | 中国矿业大学(北京) | 地裂缝评价的“三图法” |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CA2854132A1 (en) * | 2011-10-31 | 2013-05-10 | Insurance Bureau Of Canada | System and method for predicting and preventing flooding |
-
2018
- 2018-02-24 CN CN201810157583.7A patent/CN108280553B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101726750A (zh) * | 2009-11-30 | 2010-06-09 | 中国矿业大学(北京) | 地裂缝评价的“三图法” |
Non-Patent Citations (4)
Title |
---|
Elman递归神经网络在车牌字符识别中的应用;杨晓燕 等;《中北大学学报(自然科学版)》;20121015;第33卷(第5期);第3-4节 * |
基于GIS的喜马拉雅中部地区地质灾害危险性评价研究;李晓乐;《中国优秀硕士学位论文全文数据库 基础科学辑》;20130115;第四,五章 * |
基于GIS的闽江上游滑坡敏感性评价研究;陈齐;《中国优秀硕士学位论文全文数据库 基础科学辑》;20170515;第5.3节 * |
模糊综合评价和改进神经网络相结合的电力建设项目合同风险评价;李泽红等;《科技和产业》;20090325;第9卷(第03期);第1-3节 * |
Also Published As
Publication number | Publication date |
---|---|
CN108280553A (zh) | 2018-07-13 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108280553B (zh) | 基于gis-神经网络集成的山洪灾害风险区划及预测方法 | |
CN109858647B (zh) | 一种耦合gis和gbdt算法的区域洪涝灾害风险评价与预估方法 | |
CN106650767B (zh) | 基于聚类分析和实时校正的洪水预报方法 | |
Zou et al. | Comprehensive flood risk assessment based on set pair analysis-variable fuzzy sets model and fuzzy AHP | |
CN112819207B (zh) | 基于相似性度量的地质灾害空间预测方法、系统及存储介质 | |
CN111582386A (zh) | 一种基于随机森林的地质灾害多灾种综合风险评价方法 | |
Minglei et al. | Classified real-time flood forecasting by coupling fuzzy clustering and neural network | |
CN111665575B (zh) | 一种基于统计动力的中长期降雨分级耦合预报方法及系统 | |
CN113408917B (zh) | 自然灾害风险评估方法、装置、计算机设备及存储介质 | |
CN106934237A (zh) | 雷达抗干扰效能评估可信性度量实现方法 | |
CN105843829A (zh) | 一种基于分层模型的大数据可信性度量方法 | |
CN115239127A (zh) | 生态脆弱性评价方法、计算机设备、存储介质和验证方法 | |
CN116596303A (zh) | 一种干旱风险评估与区划方法、系统、介质、设备及终端 | |
CN115829326A (zh) | 基于优化组合赋权模型的山区公路建设风险评价方法 | |
CN113191582B (zh) | 一种基于gis与机器学习的道路山洪易损性评价方法 | |
Li et al. | Drought assessment based on data fusion and deep learning | |
CN117540303A (zh) | 基于交叉半监督机器学习算法的滑坡易发性评估方法及系统 | |
CN117350146A (zh) | 一种基于ga-bp神经网络的排水管网健康性评价方法 | |
CN112434887A (zh) | 一种结合网络核密度估计和svm的供水管网风险预测方法 | |
CN115829209A (zh) | 基于碳路径的绿色智能仓库环保质量分析方法及装置 | |
Jin et al. | Integration of an improved transformer with physical models for the spatiotemporal simulation of urban flooding depths | |
CN115906669A (zh) | 一种顾及负样本选取策略的密集残差网络滑坡易发性评价方法 | |
CN115829150A (zh) | 一种积水预测系统 | |
CN115358587A (zh) | 区域多部门协同的基础设施规划方法及系统 | |
Owolabi et al. | Development of priority index assessment model for road pavements in Nigeria |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |