CN111861133A - 一种山洪灾害防治能力评价方法 - Google Patents
一种山洪灾害防治能力评价方法 Download PDFInfo
- Publication number
- CN111861133A CN111861133A CN202010585403.2A CN202010585403A CN111861133A CN 111861133 A CN111861133 A CN 111861133A CN 202010585403 A CN202010585403 A CN 202010585403A CN 111861133 A CN111861133 A CN 111861133A
- Authority
- CN
- China
- Prior art keywords
- weight
- index
- evaluation
- consistency
- indexes
- 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.)
- Pending
Links
- 238000011156 evaluation Methods 0.000 title claims abstract description 56
- 230000002265 prevention Effects 0.000 title claims abstract description 33
- 238000000034 method Methods 0.000 claims abstract description 58
- 238000013278 delphi method Methods 0.000 claims abstract description 6
- 238000004458 analytical method Methods 0.000 claims abstract description 5
- 238000013210 evaluation model Methods 0.000 claims abstract description 5
- 238000012821 model calculation Methods 0.000 claims abstract description 4
- 238000011002 quantification Methods 0.000 claims abstract description 4
- 239000011159 matrix material Substances 0.000 claims description 54
- 238000004364 calculation method Methods 0.000 claims description 47
- 239000013598 vector Substances 0.000 claims description 20
- 230000008569 process Effects 0.000 claims description 14
- 239000002131 composite material Substances 0.000 claims description 13
- 238000012360 testing method Methods 0.000 claims description 9
- 238000010606 normalization Methods 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 5
- 238000006243 chemical reaction Methods 0.000 claims description 3
- 230000000694 effects Effects 0.000 abstract description 6
- 238000007726 management method Methods 0.000 abstract description 4
- 238000013139 quantization Methods 0.000 description 3
- 239000000470 constituent Substances 0.000 description 2
- 230000007123 defense Effects 0.000 description 2
- 238000012544 monitoring process Methods 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000012854 evaluation process Methods 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 230000008520 organization Effects 0.000 description 1
- 238000004451 qualitative analysis Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000012163 sequencing technique Methods 0.000 description 1
- 239000000243 solution Substances 0.000 description 1
- 239000012086 standard solution Substances 0.000 description 1
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/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
-
- 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
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- 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)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Entrepreneurship & Innovation (AREA)
- General Physics & Mathematics (AREA)
- Marketing (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Theoretical Computer Science (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种山洪灾害防治能力评价方法,步骤1,采集待分析区域的数据作为初始数据,确定评价指标;步骤2,对于定量指标直接进行赋值;对于定性指标利用德尔菲法进行指标量化后进行赋值;步骤3,利用主观赋权采用层次分析法、客观赋权采用熵权法来确定各评价指标的权重,利用主观赋权法和客观赋权法这两者结合来确定综合权重;步骤4,建立山洪灾害防治能力评估模型,依据模型计算结果进行分析评价。本发明通过建立评价体系,综合考虑山洪灾害的风险性和山洪灾害抵御能力,对山洪灾害的防治能力进行评价,有利于提升政府的灾害管理水平,提高山丘区洪灾的预防与治理能力,也可以更直观地了解山洪灾害防治措施发挥的作用。
Description
技术领域
本发明涉及应急防灾领域,具体一种山洪灾害防治能力评价指标体系与山洪灾害防治能力评价方法。
背景技术
山洪灾害是指由于降雨在山丘区引发的洪水及由山洪诱发的泥石流、滑坡等对国民经济和人民生命财产造成损失的灾害。复杂的地形地质条件、暴雨多发的气候特征、密集的人口分布和人类活动的影响,导致我国山洪灾害发生频繁。山区的人类活动集中,城镇和重大工程沿河谷展布,人口、经济密集区与山洪灾害危险区高度重叠,山洪灾害不仅使经济受到严重损失而且还造成了大量的人员伤亡。现阶段我国已初步建立了山洪灾害防治体系,并取得了较好的成效,但随着山丘区社会经济的快速发展与群众对山洪预警公共服务的迫切需求,防治工作又被提出了新的更高要求,山洪灾害一旦发生将会使山丘区受到更加严重的迫害,所以山洪灾害防治措施需要不断补充建设。
现在很多学者利用GIS和RS等技术对研究区的危险性和易损性进行定性、定量的分析,对山洪灾害风险度的评价、山洪灾害风险区划及相应的风险图的编制已经比较成熟,但在山洪灾害防治能力的评价过程中存在模糊性与不相容问题,采用模糊物元方法建立模型可以很好解决这个问题,因此建立该评价指标体系具有现实意义。
发明内容
本发明的目的在于提供一种山洪灾害防治能力评价方法,旨在建立山洪灾害防治能力评价模型,为水行政主管部门提供定量依据,提升政府的灾害管理水平,提高山丘区洪灾的预防与治理能力。
本发明的一种山洪灾害防治能力评价方法,包括以下步骤:
步骤1,采集待分析区域的数据作为初始数据,确定评价指标,包括28个定量指标和7个定性指标,并对指标进行分类:第一类是通过查询项目相关资料或已公布的统计数据获取的统计类情景评价指标;第二类是采用水文水利计算、损失评估方法获取的计算类评价指标;第三类是需由决策者做出主观判断的定性指标;
步骤2,对于定量指标直接进行赋值;对于定性指标利用德尔菲法进行指标量化后进行赋值;
步骤3,利用主观赋权采用层次分析法、客观赋权采用熵权法来确定各评价指标的权重,利用主观赋权法和客观赋权法这两者结合来确定综合权重;
步骤4,建立山洪灾害防治能力评估模型,具体包括以下过程:
步骤4-1,构建n个事物的m维复合模糊物元Rmn,表达式为:
式中,Rmn为n个事物的m个模糊特征的复合物元,Mn为第n个事物,Cj为第j个特征向量;xmn为第n个事物的第m个模糊特征对应的模糊量值;
步骤4-3,构建差平方复合模糊物元RΔ,表达式如下:
步骤4-4,依据欧式贴近度公式计算模糊算子并构建欧式贴近度复合模糊物元RρH,表达式如下:
依据模型计算结果进行分析评价。
所述步骤2中利用德尔菲法进行指标量化的步骤具体包括以下处理:
对定性指标量化;
N位专家对量化值打分;
如果打分结果是较好,则根据打分结果修改正量化值;否则,重新打分。
所述步骤3中确定综合权重的步骤具体包括以下处理:
3.1、层次分析法主观赋权,其具体的操作步骤如下:
(311)建立层次结构:将问题所含的因素按照相互关系进行分层,各因素分别归入不同的层次结构,以框架结构说明各层次之间的从属关系;
(312)构造判断矩阵:通过专家打分,把某一层中的要素与高一层次因素之间的相对程度用矩阵表示出来;
(313)计算层次单排序权重:层次单排序的任务归结为计算判断矩阵的特征根和特征向量
其满足下式:
Aw=λmaxW
其中,λmax为判断矩阵A的最大特征根,W为对应于λmax的正规化特征向量,W的分量Wi为对应元素单排序的权重值;
(314)计算层次单排序一致性指标CI层次单排序、层次单排序随机一致性比例CR层次单排序,公式如下:
其中,λmax为判断矩阵A的最大特征根、N是矩阵维数;
依据层次单排序一致性指标CI层次单排序、层次单排序随机一致性比例C层次单排序调整判断矩阵,直至判断矩阵是一致性矩阵;
计算各个指标在整个指标体系中的总排序,计算公式如下:
Wxi=Wx×Wx,xi
其中,Wx为一级指标x在整个指标体系中的权重,Wx,xi为二级指标xi在x层次上的权重;
为进行一致性检验,计算层次总排序一致性指标CI、层次总排序平均一致性指标RI、层次总排序随机一致性比例CR,公式如下:
CR=CI/RI
式中:Waj表示A层次中的元素Aj的层次总排序权重;CI为层次总排序一致性指标;CIj为与Aj对应下一层次中判断矩阵的随机一致性指标;CR为层次总排序随机一致性比例;同样,当CR<0.1时,则认为层次总排序的计算结果具有令人满意的一致性;否则就需对本层次的各判断矩阵进行调整,从而使层次总排序具有令人满意的一致性;
3.2利用熵权法进行客观赋权,其计算步骤如下:
(321)构建n个区域m个评估指标的判断矩阵:
Ri×j=(xij)(i=1,2,…,m;j=1,2,…,n);
(322)判断矩阵归一化处理,得到归一化矩阵B:
其中,maxxi、minxj分别为同一指标下不同对象中最满意者或最不满意者;
(323)根据熵的定义,利用n个评估区域m个评估指标构造的判断矩阵可以确实指标的熵Hj为:
(324)指标差异性系Gi
Gi=1-Hi;
(325)熵权重Wi确定:
3.3、采用Kendall协和系数来检验主观权重与客观权重的一致性:首先把确定的权重向量转换成按评价值的高低顺序向量,其计算步骤如下:
(331)假设H0:n种权重计算方法计算结果不具有一致性;假设H1:n种权重计算方法计算结果具有一致性,给定显著水平α;
(332)将i组权重记作Wi=(W1 i,W2 i,…,Wm i),转换成高低排序值:
Pij(1≤pij≤m)(i=1,2,…,m)
其中,m为指标数;
(333)计算Kendall协和系数检验统计量,公式如下:
其中,n为权重组数;
(334)依据显著水平α,查阅Kendall协和系数检验的临界表,确定临界值;当K≤Kα时,接收H0,拒绝H1,即n中权重计算方法所赋权值无一致性;K≧Kα,拒绝H0,接收H1,即n中权重计算方法所赋权值具有一致性,是合理的;
3.4、基于离差平方和的综合权重计算方法进行综和权重计算,计算公式如下:
W=ρ1W1+ρ2W2+…+ρlWl
其中,W为综合权重向量,ρ1,ρ2,…,ρl为不同方法计算所得权重的组合系数;且权重系数满足:
与现有技术相比,本发明的有益效果如下:
通过建立评价体系,综合考虑山洪灾害的风险性和山洪灾害抵御能力,对山洪灾害的防治能力进行评价,有利于提升政府的灾害管理水平,提高山丘区洪灾的预防与治理能力,也可以更直观地了解山洪灾害防治措施发挥的作用,发现项目中存在的不足与问题,为下阶段地山洪防治建设提供借鉴与依据,促进山洪灾害防治长效机制的建立,丰富山洪灾害防治研究体系。
附图说明
图1是本发明的一种山洪灾害防治能力评价方法整体流程图;
图2是利用德尔菲法对定性指标值进行指标量化流程图;
图3是确定综合权重流程图。
具体实施方式
以下结合具体实施例以及附图,对本发明的技术方案进行说明。
如图1所示,是本发明的一种山洪灾害防治能力评价方法流程框图,如图1所示:
步骤1,根据待评估单元,选取合理的指标,构建评价体系。采集待分析区域的相关数据作为初始数据,包括森林覆盖率、年平均暴雨量、24小时最大降雨量、危险区暴露人口占比、危险区房屋价值、山洪灾害管理机构密度等35个指标,其中包括28个定量指标和7个定性指标。指标可分为三类,第一类是可以通过查询项目相关资料或已公布的统计数据获取的统计类情景评价指标;第二类是采用水文水利计算、损失评估等方法获取的计算类评价指标;第三类是无法通过统计或计算得到,需由决策者做出主观判断的定性指标。35个评估指标共分为7类,作为模型的7个准则层,分别是山洪灾害危险性、山洪灾害易损性、灾害预防能力、监测预警能力、群测群防能力、工程防御能力和应急救援能力本实施例选取云南省7个灾害频发且有已建山洪灾害防治项目的县(市)作为评估单元,根据当地情况,最终筛选出包括森林覆盖率、年平均暴雨量、危险区暴露人口占比、水雨情监测站点密度等21项指标(见表1)来构建山洪灾害防治能力指标体系,同时收集相应指标在山洪灾害风险评估中的标准值;
表1
步骤2,对于定量指标直接进行赋值;对于定性指标利用德尔菲法进行指标量化后进行赋值,本例选取7个样本县2014年数据作为样本数据使用;
步骤3,确定各评价指标的权重,即使用层次分析法进行主观赋权、使用熵权法客观赋权这两者结合,确定综合权重;其中主观赋权法采用层次分析法主观赋权,客观赋权法采用熵权法客观赋权;
3.1层次分析法主观赋权
层次分析法采用一种多目标、多标准的两两比较的方法,列出评估指标集中各个评估指标的先后顺序与权重系数。
具体的操作步骤如下:
(1)建立层次结构。将问题所含的因素,按照相互关系进行分层,各因素可分别归入不同的层次结构,以框架结构说明各层次之间的从属关系。
(2)构造判断矩阵。通过专家打分,把某一层中的要素,与高一层次因素之间的相对程度用矩阵表示出来,以便比较本层次的各因素与某一因素之间的相对程度。假设B层次中的组成元素B1,B2,B3…,Bn和上一层次A中的组成元素Ak有关系,则可以运用判断矩阵将其关系表示出来,如表2所示,为判断矩阵的形式。
表2
不同因子间的相互比较采用1-9标度方法进行打分,不同重要程度分别赋予不同分值,如表3所示,为标度值含义表。
表3
(3)计算层次单排序权重。层次单排序是对于上层次中的某元素而言,确定本层次与之有联系的元素重要性次序的权重值,是本层次所有元素对上一层次而言的重要性排序。层次单排序的任务就可以归结为计算判断矩阵的特征跟和特征向量的问题,及对于判断矩阵A,计算满足Aw=λmaxWλmax为A的最大特征根,W为对应于λmax的正规化特征向量,W的分量Wi就是对应元素单排序的权重值。
当一致性指标CI=0时,说明判断矩阵是完全一致性矩阵,一致性指标CI数值越小,则判断矩阵一致性就越好。一般而言1阶或2阶的判断矩阵就是完全一致性矩阵。若判断矩阵为2阶以上,则矩阵的一致性指标CI与同阶的平均一致性指标RI比值,称为判断矩阵的随机一致性比例CR。一般只有当CR<0.1时认为判断矩阵是令人满意的一致性矩阵;否则就需要继续调整判断矩阵,直到达到令人满意为止。
RI与N的随机一致性查算表结果如表4所示:
表4
(3)计算层次总排序权重。利用统一层次中所有层次单排序的结果,就可以计算针对上一层次而言的本层次所有元素的重要性权重值,这就是层次总排序。层次总排序需要从上到下逐层进行,各个指标在整个指标体系中的总排序计算方法如下:
Wxi=Wx×Wx,xi
式中:Wx是一级指标x在整个指标体系中的权重,Wx,xi是二级指标xi在x层次上的权重,x=A,B,C,D,E,…(分别代表一级指标)。为进行一致性检验,需计算下列指标:
CR=CI/RI
式中:Waj判断矩阵A中的元素aj的层次总排序权重;CI为层次总排序一致性指标;CIj为与aj对应下一层次中判断矩阵的随机一致性指标;CR为层次总排序随机一致性比例。同样,当CR<0.1时,则认为层次总排序的计算结果具有令人满意的一致性;柔则就需对本层次的各判断矩阵进行调整,从而使层次总排序具有令人满意的一致性。
3.2熵权法客观赋权
熵权法属于客观赋权法,通过已知评估指标数值构造判断矩阵,根据各项指标值的变异程度,来确定权重,从而得到较为客观的评估结果。该方法有效避免权重确定过程中主观意念的影响,使评估结果更加符合客观实际。其计算步骤如下:
(1)构建n个区域m个评估指标的判断矩阵:
Ri×j=(xij)(i=1,2,…,m;j=1,2,…,n)。
(2)判断矩阵归一化处理,得到归一化矩阵B:
maxxi、minxj分别为同一指标下不同对象中最满意者或最不满意者。
(3)根据熵的定义,利用n个评估区域m个评估指标构造的判断矩阵可以确实指标的熵Hj为:
(4)指标差异性系Gi
Gi=1-Hi
(5)熵权重Wi确定:
由以上公式可知,熵权重就是根据判断矩阵中指标量值的差异性来确定的。单项指标数值差异性越大,指标权重越大,评估指标对评估对象的作用就越大。
3.3、主客观权重一致性检验
本发明采用Kendall协和系数来检验主观权重与客观权重的一致性。
Kendall协和系数计算法属于一种秩相关的分析方法。在检验之前,首先要把两种方法确定的权重向量转换成按评价值的高低顺序向量,其计算步骤如下:
(1)设H0表示权重计算方法计算结果不具有一致性;设H1表示权重计算方法计算结果具有一致性。给定显著水平α。
(2)将i组权重记作Wi=(W1 i,W2 i,…,Wm i),转换成高低排序值:
Pij(1≤pij≤m)(i=1,2,…,m)
式中:n为权重组数目;m为指标数目。
(3)计算Kendall协和系数检验统计量:
(4)依据显著水平α,查阅Kendall协和系数检验的临界表,确定临界值。当K≤Kα时,接收H0,拒绝H1,即权重计算方法所赋权值无一致性;K≧Kα,拒绝H0,接收H1,即权重计算方法所赋权值具有一致性,是合理的。
3.4、基于离差平方和的综和权重计算
本文选取基于离差平方和的综合权重计算方法,该方法能融合主客观方法的优点,原理简单易理解,计算方便快捷。其计算公式如下:
W=ρ1W1+ρ2W2+…+ρlWl
其中,W为综合权重向量,ρ1,ρ2,…,ρl为不同方法计算所得权重的组合系数。且权重系数满足:
步骤4,建立山洪灾害防治能力评估模型,具体包括以下过程:
步骤4-1,理论中所描述的事物M及其特征向量C和量值x组成物元R=(M,C,x)或R=(M,C,C(M)),n个事物的m维模糊物元组合到一起,构建n个事物的m维复合模糊物元Rmn,表达式如下:
式中,Rmn为n个事物的m个模糊特征的复合物元,Mn为第n个事物,Cj为第j个特征向量;xmn为第n个事物的第m个模糊特征对应的模糊量值;
步骤4-3,构建差平方复合模糊物元RΔ,表达式如下:
步骤4-4,依据欧式贴近度公式计算模糊算子并构建欧式贴近度复合模糊物元RρH,表达式如下:
依据模型计算结果进行分析评价。
表5
以上为对本申请所提供的一种山洪灾害防治能力评价指标体系与评价方法所做的介绍。以上实施步骤仅用于帮助理解本发明的具体方法与核心思想,并不用以限制本发明。对于本领域的一些研究人员,依据本申请的思想,在具体实施方式及应用范围上均会有改变之处。
如图2所示,是利用德尔菲法进行指标量化流程图。该流程中:
对定性指标量化;
N位专家对量化值打分(N≥6);
如果打分结果是较好,则根据打分结果修改正量化值;否则,重新打分。
如图3所示,是确定综合权重流程图。该流程中:
对层次分析法主观赋权、熵权法客观赋权分别进行主客观权重一致性检验;基于离差平方和的综和权重计算;确定综合权重。
Claims (3)
1.一种山洪灾害防治能力评价方法,其特征在于,该方法包括以下步骤:
步骤1,采集待分析区域的数据作为初始数据,确定评价指标,包括28个定量指标和7个定性指标,并对指标进行分类:第一类是通过查询项目相关资料或已公布的统计数据获取的统计类情景评价指标;第二类是采用水文水利计算、损失评估方法获取的计算类评价指标;第三类是需由决策者做出主观判断的定性指标;
步骤2,对于定量指标直接进行赋值;对于定性指标利用德尔菲法进行指标量化后进行赋值;
步骤3,利用主观赋权采用层次分析法、客观赋权采用熵权法来确定各评价指标的权重,利用主观赋权法和客观赋权法这两者结合来确定综合权重;
步骤4,建立山洪灾害防治能力评估模型,具体包括以下过程:
步骤4-1,构建n个事物的m维复合模糊物元Rmn,表达式为:
式中,Rmn为n个事物的m个模糊特征的复合物元,Mn为第n个事物,Cj为第j个特征向量;xmn为第n个事物的第m个模糊特征对应的模糊量值;
步骤4-3,构建差平方复合模糊物元RΔ,表达式如下:
步骤4-4,依据欧式贴近度公式计算模糊算子并构建欧式贴近度复合模糊物元RρH,表达式如下:
依据模型计算结果进行分析评价。
2.如权利要求1所述的一种山洪灾害防治能力评价方法,其特征在于,所述步骤2中利用德尔菲法进行指标量化的步骤具体包括以下处理:
对定性指标量化;
N位专家对量化值打分;
如果打分结果是较好,则根据打分结果修改正量化值;否则,重新打分。
3.如权利要求1所述的一种山洪灾害防治能力评价方法,其特征在于,所述步骤3中确定综合权重的步骤具体包括以下处理:
3.1、层次分析法主观赋权,其具体的操作步骤如下:
(311)建立层次结构:将问题所含的因素按照相互关系进行分层,各因素分别归入不同的层次结构,以框架结构说明各层次之间的从属关系;
(312)构造判断矩阵:通过专家打分,把某一层中的要素与高一层次因素之间的相对程度用矩阵表示出来;
(313)计算层次单排序权重:层次单排序的任务归结为计算判断矩阵的特征根和特征向量
其满足下式:
Aw=λmaxW
其中,λmax为判断矩阵A的最大特征根,W为对应于λmax的正规化特征向量,W的分量Wi为对应元素单排序的权重值;
(314)计算层次单排序一致性指标CI层次单排序、层次单排序随机一致性比例CR层次单排序,公式如下:
其中,λmax为判断矩阵A的最大特征根、N是矩阵维数;
依据层次单排序一致性指标CI层次单排序、层次单排序随机一致性比例C层次单排序调整判断矩阵,直至判断矩阵是一致性矩阵;
计算各个指标在整个指标体系中的总排序,计算公式如下:
Wxi=Wx×Wx,xi
其中,Wx为一级指标x在整个指标体系中的权重,Wx,xi为二级指标xi在x层次上的权重;
为进行一致性检验,计算层次总排序一致性指标CI、层次总排序平均一致性指标RI、层次总排序随机一致性比例CR,公式如下:
CR=CI/RI
其中,Waj为判断矩阵A中的元素aj的层次总排序权重;
依据次总排序一致性指标CI、层次总排序平均一致性指标RI、层次总排序随机一致性比例CR对本层次的各判断矩阵进行调整,从而使层次总排序具有一致性;
3.2、利用熵权法进行客观赋权,其计算步骤如下:
(321)构建n个区域m个评估指标的判断矩阵:
Ri×j=(xij)(i=1,2,…,m;j=1,2,…,n);
(322)判断矩阵归一化处理,得到归一化矩阵B:
其中,maxxi、minxj分别为同一指标下不同对象中最满意者或最不满意者;
(323)根据熵的定义,利用n个评估区域m个评估指标构造的判断矩阵可以确实指标的熵Hj为:
(324)指标差异性系Gi
Gi=1-Hi;
(325)熵权重Wi确定:
3.3、采用Kendall协和系数来检验主观权重与客观权重的一致性:首先把确定的权重向量转换成按评价值的高低顺序向量,计算步骤如下:
(331)设H0表示权重计算方法计算结果不具有一致性;设H1表示权重计算方法计算结果具有一致性,给定显著水平α;
(332)将i组权重记作Wi=(W1 i,W2 i,…,Wm i),转换成高低排序值:
Pij(1≤pij≤m)(i=1,2,…,m)
其中,m为指标数;
(333)计算Kendall协和系数检验统计量,公式如下:
其中,n为权重组数;
(334)依据显著水平α,查阅Kendall协和系数检验的临界表,确定临界值;当K≤Kα时,接收H0,拒绝H1,即n中权重计算方法所赋权值无一致性;K≧Kα,拒绝H0,接收H1,即n中权重计算方法所赋权值具有一致性,是合理的;
3.4、基于离差平方和的综合权重计算方法进行综和权重计算,计算公式如下:
W=ρ1W1+ρ2W2+…+ρlWl
其中,W为综合权重向量,ρ1,ρ2,…,ρl为不同方法计算所得权重的组合系数;且权重系数满足:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010585403.2A CN111861133A (zh) | 2020-06-23 | 2020-06-23 | 一种山洪灾害防治能力评价方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202010585403.2A CN111861133A (zh) | 2020-06-23 | 2020-06-23 | 一种山洪灾害防治能力评价方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN111861133A true CN111861133A (zh) | 2020-10-30 |
Family
ID=72989353
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202010585403.2A Pending CN111861133A (zh) | 2020-06-23 | 2020-06-23 | 一种山洪灾害防治能力评价方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN111861133A (zh) |
Cited By (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112580078A (zh) * | 2020-12-22 | 2021-03-30 | 江苏方天电力技术有限公司 | 一种针对电力数据指纹的评估方法 |
CN112633679A (zh) * | 2020-12-21 | 2021-04-09 | 贵州电网有限责任公司电力科学研究院 | 信息质量量化方法、装置、计算机设备和存储介质 |
CN112862266A (zh) * | 2021-01-19 | 2021-05-28 | 昆明理工大学 | 一种基于模糊决策的煤矿事故应急预案评估模型构建方法 |
CN113112124A (zh) * | 2021-03-22 | 2021-07-13 | 西安理工大学 | 一种淤地坝系风险评价方法 |
CN113222376A (zh) * | 2021-04-29 | 2021-08-06 | 北京建筑大学 | 文物自然灾害风险监测技术适宜性评价方法和装置 |
CN113408928A (zh) * | 2021-06-30 | 2021-09-17 | 广东电网有限责任公司 | 一种配电自动化终端评价方法及系统 |
CN113496359A (zh) * | 2021-07-13 | 2021-10-12 | 天津大学 | 一种多组合主客观均匀优化赋权风险评价方法 |
CN113837621A (zh) * | 2021-09-27 | 2021-12-24 | 中国电子科技集团公司第二十八研究所 | 一种机场飞行区安全风险综合评价方法 |
CN113869354A (zh) * | 2021-08-16 | 2021-12-31 | 天津大学 | 一种凌汛期堤防危险性综合评价方法 |
CN113934904A (zh) * | 2021-09-24 | 2022-01-14 | 广东烟草惠州市有限责任公司 | 一种基于rfm模型的卷烟零售客户价值评价方法 |
CN114202257A (zh) * | 2022-02-18 | 2022-03-18 | 中国科学院地理科学与资源研究所 | 自然灾害综合减灾能力评估方法及装置 |
CN116798633A (zh) * | 2023-08-22 | 2023-09-22 | 北京大学人民医院 | 创伤数据安全风险评估系统构建方法和电子设备 |
-
2020
- 2020-06-23 CN CN202010585403.2A patent/CN111861133A/zh active Pending
Non-Patent Citations (3)
Title |
---|
杨通通: "山洪灾害防治能力与效果综合评估模型及应用研究", 《中国优秀硕士学位论文全文数据库》 * |
潘汀超 等: "组合赋权_模糊聚类算法的改进及其在洪灾风险评价的应用", 《南水北调与水利科技(中英文)》 * |
苑希民 等: "基于熵权法的模糊物元模型在山洪防御能力评估中的应用", 《水利水电技术》 * |
Cited By (17)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112633679A (zh) * | 2020-12-21 | 2021-04-09 | 贵州电网有限责任公司电力科学研究院 | 信息质量量化方法、装置、计算机设备和存储介质 |
CN112580078A (zh) * | 2020-12-22 | 2021-03-30 | 江苏方天电力技术有限公司 | 一种针对电力数据指纹的评估方法 |
CN112580078B (zh) * | 2020-12-22 | 2022-11-08 | 江苏方天电力技术有限公司 | 一种针对电力数据指纹的评估方法 |
CN112862266B (zh) * | 2021-01-19 | 2022-08-09 | 昆明理工大学 | 一种基于模糊决策的煤矿事故应急预案评估模型构建方法 |
CN112862266A (zh) * | 2021-01-19 | 2021-05-28 | 昆明理工大学 | 一种基于模糊决策的煤矿事故应急预案评估模型构建方法 |
CN113112124A (zh) * | 2021-03-22 | 2021-07-13 | 西安理工大学 | 一种淤地坝系风险评价方法 |
CN113222376A (zh) * | 2021-04-29 | 2021-08-06 | 北京建筑大学 | 文物自然灾害风险监测技术适宜性评价方法和装置 |
CN113408928A (zh) * | 2021-06-30 | 2021-09-17 | 广东电网有限责任公司 | 一种配电自动化终端评价方法及系统 |
CN113496359B (zh) * | 2021-07-13 | 2022-03-22 | 天津大学 | 一种多组合主客观均匀优化赋权风险评价方法 |
CN113496359A (zh) * | 2021-07-13 | 2021-10-12 | 天津大学 | 一种多组合主客观均匀优化赋权风险评价方法 |
CN113869354A (zh) * | 2021-08-16 | 2021-12-31 | 天津大学 | 一种凌汛期堤防危险性综合评价方法 |
CN113869354B (zh) * | 2021-08-16 | 2024-05-07 | 天津大学 | 一种凌汛期堤防危险性综合评价方法 |
CN113934904A (zh) * | 2021-09-24 | 2022-01-14 | 广东烟草惠州市有限责任公司 | 一种基于rfm模型的卷烟零售客户价值评价方法 |
CN113837621A (zh) * | 2021-09-27 | 2021-12-24 | 中国电子科技集团公司第二十八研究所 | 一种机场飞行区安全风险综合评价方法 |
CN114202257A (zh) * | 2022-02-18 | 2022-03-18 | 中国科学院地理科学与资源研究所 | 自然灾害综合减灾能力评估方法及装置 |
CN116798633A (zh) * | 2023-08-22 | 2023-09-22 | 北京大学人民医院 | 创伤数据安全风险评估系统构建方法和电子设备 |
CN116798633B (zh) * | 2023-08-22 | 2023-11-21 | 北京大学人民医院 | 创伤数据安全风险评估系统构建方法和电子设备 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN111861133A (zh) | 一种山洪灾害防治能力评价方法 | |
CN106651211B (zh) | 一种不同尺度区域洪水灾害风险评估的方法 | |
CN108710984B (zh) | 一种矿山地质环境综合评价方法及系统 | |
Yang et al. | Multiple flood vulnerability assessment approach based on fuzzy comprehensive evaluation method and coordinated development degree model | |
CN108280553B (zh) | 基于gis-神经网络集成的山洪灾害风险区划及预测方法 | |
CN107544253B (zh) | 基于改进模糊熵权法的大型导弹装备退役安全控制方法 | |
CN112287018A (zh) | 一种台风灾害下10kV杆塔受损风险评估方法及系统 | |
CN113191644A (zh) | 一种基于自发-引发风险评价模型城市内涝风险评价方法 | |
CN113033997A (zh) | 基于改进集对分析的城市水质等级确定方法、装置及介质 | |
CN111047043A (zh) | 一种基于云模型对堰塞坝风险等级模糊综合评价分析方法 | |
CN114881396A (zh) | 一种基于ahp和topsis的隧道坍塌风险评估方法 | |
CN113128893A (zh) | 区域防旱减灾的评估方法和装置 | |
CN108269016A (zh) | 一种基于信息扩散的小流域山洪灾害风险分析方法 | |
CN113191582B (zh) | 一种基于gis与机器学习的道路山洪易损性评价方法 | |
CN113505978A (zh) | 一种不同形式城市社区的防灾功能评估方法和装置 | |
CN107169289A (zh) | 一种基于可拓最优组合赋权法的滑坡危险性评价方法 | |
CN111401677A (zh) | 生态水利工程的综合评价方法 | |
CN113516396A (zh) | 基于模糊层次分析的城市动态风险指数评估方法 | |
CN113657732A (zh) | 一种基于层次分析法的灾度指数计算方法和设备 | |
CN115829326A (zh) | 基于优化组合赋权模型的山区公路建设风险评价方法 | |
CN111953543A (zh) | 一种基于pca-ahp的量子通信网络可靠性状况的评估方法 | |
CN111062588A (zh) | 基于层次分析的区域电网设备强降雨灾害影响评估方法 | |
CN116757368A (zh) | 一种基于gis的应急避难场所评估及选址方法和系统 | |
CN114444910A (zh) | 一种面向电力物联网的边缘网络系统健康度评估方法 | |
CN108805471A (zh) | 基于复合系统作用关系分析的水资源承载能力评价方法 |
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 | ||
WD01 | Invention patent application deemed withdrawn after publication |
Application publication date: 20201030 |
|
WD01 | Invention patent application deemed withdrawn after publication |