[发明专利]一种基于广义熵理论的水环境风险管理动态智能综合分析方法在审

专利信息
申请号: 201510400221.2 申请日: 2015-07-09
公开(公告)号: CN104992245A 公开(公告)日: 2015-10-21
发明(设计)人: 赵君;徐进超;朱节中;杭庆丰;曹莉莉;陈雅倩;战建;陈秀英 申请(专利权)人: 南京信息工程大学
主分类号: G06Q10/04 分类号: G06Q10/04;G06Q50/06
代理公司: 南京瑞弘专利商标事务所(普通合伙) 32249 代理人: 杨晓玲
地址: 210044 *** 国省代码: 江苏;32
权利要求书: 查看更多 说明书: 查看更多
摘要: 发明公开了一种基于广义熵理论的水环境风险管理动态智能综合分析方法。本方法利用水利科学、智能科学、信息科学与风险理论等集成的途径,从整体结构和组成部分出发,对水环境风险管理中水量、水质、水生态等风险因素进行探讨,提出利用广义分布函数及其广义熵理论统一描述、物理解析系统中各种不确定性信息,构建基于遗传算法、模糊层次分析、投影寻踪的水环境风险管理广义熵动态智能综合分析方法。采用本方法对水环境风险进行分析得到的模拟与实测预警结果相符,准确率较高,这说明本方法能够在一定程度上合理地利用指标对水环境风险进行管理,可以应用于水环境评估方面的研究,对于水环境的科学预警具有重要意义。
搜索关键词: 一种 基于 广义 理论 水环境 风险 管理 动态 智能 综合分析 方法
【主权项】:
一种基于广义熵理论的水环境风险管理动态智能综合分析方法,其特征在于,包括如下步骤:步骤A,建立筛选子模型:首先,应用广义熵理论建立水环境风险评价指标体系,再从数据库直接获取水文专家对所述水环境风险评价指标体系中两两评价指标间重要性的判断,构建模糊互补判断矩阵,计算所述互补判断矩阵的评价指标权重;然后选取评价指标权重中均值最大的Mj个评价指标,再从所述Mj个评价指标中选取标准差最小的Nj个评价指标,组成最终的水环境风险评价指标体系;步骤B,建立预测子模型:首先根据研究区域来选择预测子模型,然后针对研究区域,应用预测子模型模拟计算水环境风险评价指标值;步骤C,建立评价子模型:根据步骤A建立的水环境风险评价指标体系,将步骤B得出的各评价指标的时空立体数据进行无量纲化处理,无量纲化公式为:<mrow><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>=</mo><mrow><mo>(</mo><msubsup><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow><mn>0</mn></msubsup><mo>-</mo><msubsup><mi>x</mi><mrow><mi>j</mi><mi>min</mi></mrow><mn>0</mn></msubsup><mo>)</mo></mrow><mo>/</mo><mrow><mo>(</mo><msubsup><mi>x</mi><mrow><mi>j</mi><mi>m</mi><mi>a</mi><mi>x</mi></mrow><mn>0</mn></msubsup><mo>-</mo><msubsup><mi>x</mi><mrow><mi>j</mi><mi>min</mi></mrow><mn>0</mn></msubsup><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>1</mn><mo>)</mo></mrow></mrow>式中:为第i个样本的第j个评价指标值,i=1,2,...,n,j=a,2,...,m;n为样本个数,m为评价指标个数;分别为第j个评价指标的样本最大值和最小值;设a为m维单位投影方向向量,其分量为a1,a2,...,am,则xij的一维投影特征值zi描述为:<mrow><msub><mi>z</mi><mi>i</mi></msub><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>j</mi><mo>=</mo><mn>1</mn></mrow><mi>m</mi></munderover><msub><mi>a</mi><mi>j</mi></msub><msub><mi>x</mi><mrow><mi>i</mi><mi>j</mi></mrow></msub><mo>,</mo><mrow><mo>(</mo><mi>i</mi><mo>=</mo><mn>1</mn><mo>,</mo><mn>2</mn><mo>,</mo><mn>...</mn><mo>,</mo><mi>n</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>2</mn><mo>)</mo></mrow></mrow>z=(z1,z2,...,zi,...,zn)为投影特征值集合;构造投影指标:设s(zi,zk)=|zi‑zk|,(k=1,2,...,n)为任意两投影特征值间的绝对值距离;将待聚类样本zi分为N(2≤N<n)类,用Θh(h=1,2,...,N)表示第h类样本投影特征值集合,即:Θh={zi|d(Ah‑zi)≤d(At‑zi)},(t=1,2,...,N,t≠h)   (3)式中:d(Ah‑zi)=|zi‑Ah|,d(At‑zi)=|zi‑At|,Ah和At分别为第h类和第t类的初始聚核;类同样本的邻近程度用类同聚集度dd(a)表示为:<mrow><msub><mi>d</mi><mi>d</mi></msub><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><munderover><mo>&Sigma;</mo><mrow><mi>h</mi><mo>=</mo><mn>1</mn></mrow><mi>N</mi></munderover><msub><mi>d</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>4</mn><mo>)</mo></mrow></mrow>式中:<mrow><msub><mi>d</mi><mi>h</mi></msub><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>z</mi><mi>i</mi></msub><mo>,</mo><msub><mi>z</mi><mi>k</mi></msub><mo>&Element;</mo><msub><mi>&Theta;</mi><mi>h</mi></msub></mrow></munder><mi>s</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>i</mi></msub><mo>,</mo><msub><mi>z</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>,</mo></mrow>样本间的离散程度用类异分散度表示为:<mrow><msub><mi>s</mi><mi>s</mi></msub><mrow><mo>(</mo><mi>a</mi><mo>)</mo></mrow><mo>=</mo><munder><mo>&Sigma;</mo><mrow><msub><mi>z</mi><mi>i</mi></msub><mo>,</mo><msub><mi>z</mi><mi>k</mi></msub><mo>&Element;</mo><mi>Z</mi></mrow></munder><mi>s</mi><mrow><mo>(</mo><msub><mi>z</mi><mi>i</mi></msub><mo>,</mo><msub><mi>z</mi><mi>k</mi></msub><mo>)</mo></mrow><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>5</mn><mo>)</mo></mrow></mrow>最后,根据动态聚类构建的投影指标表示为:QQ(a)=ss(a)‑dd(a)   (6)采用加速遗传算法对目标函数公式(6)优化求解,得到最优投影方向向量;步骤D,将步骤C得到的最优投影方向向量代入各年评价指标数据中得到水环境风险评价的投影特征值;步骤E,建立预警子模型:当步骤D中计算出的水环境风险评价的投影特征值高于2时,水环境系统处于预警状态,系统开始报警;运用对比判断法划分警度,所述运用对比判断法划分警度的具体方法为:假设E最优值是Mb,最劣值是为Ma,则隶属度公式是:<mrow><msub><mi>R</mi><mi>x</mi></msub><mo>=</mo><mfenced open = '{' close = ''><mtable><mtr><mtd><mn>1</mn></mtd><mtd><mrow><mi>E</mi><mo>&le;</mo><msub><mi>M</mi><mi>b</mi></msub></mrow></mtd></mtr><mtr><mtd><mrow><mo>(</mo><mi>E</mi><mo>-</mo><msub><mi>M</mi><mi>a</mi></msub><mo>)</mo><mo>/</mo><mo>(</mo><msub><mi>M</mi><mi>b</mi></msub><mo>-</mo><msub><mi>M</mi><mi>a</mi></msub><mo>)</mo></mrow></mtd><mtd><mrow><msub><mi>M</mi><mi>b</mi></msub><mo>&lt;</mo><mi>E</mi><mo>&lt;</mo><msub><mi>M</mi><mi>a</mi></msub></mrow></mtd></mtr><mtr><mtd><mn>0</mn></mtd><mtd><mrow><mi>E</mi><mo>&GreaterEqual;</mo><msub><mi>M</mi><mi>a</mi></msub></mrow></mtd></mtr></mtable></mfenced><mo>-</mo><mo>-</mo><mo>-</mo><mrow><mo>(</mo><mn>7</mn><mo>)</mo></mrow></mrow>式中:Rx是预警对象的隶属度,预警对象的警限对应的隶属度设为0.80,根据隶属度将警度划分为:轻警区间是Rx∈[0.60,0.80);中警区间是Rx∈[0.40,0.60);重警区间是Rx∈[0.20,0.40);巨警区间是Rx∈[0,0.20)。
下载完整专利技术内容需要扣除积分,VIP会员可以免费下载。

该专利技术资料仅供研究查看技术是否侵权等信息,商用须获得专利权人授权。该专利全部权利属于南京信息工程大学,未经南京信息工程大学许可,擅自商用是侵权行为。如果您想购买此专利、获得商业授权和技术合作,请联系【客服

本文链接:http://www.vipzhuanli.com/patent/201510400221.2/,转载请声明来源钻瓜专利网。

×

专利文献下载

说明:

1、专利原文基于中国国家知识产权局专利说明书;

2、支持发明专利 、实用新型专利、外观设计专利(升级中);

3、专利数据每周两次同步更新,支持Adobe PDF格式;

4、内容包括专利技术的结构示意图流程工艺图技术构造图

5、已全新升级为极速版,下载速度显著提升!欢迎使用!

请您登陆后,进行下载,点击【登陆】 【注册】

关于我们 寻求报道 投稿须知 广告合作 版权声明 网站地图 友情链接 企业标识 联系我们

钻瓜专利网在线咨询

周一至周五 9:00-18:00

咨询在线客服咨询在线客服
tel code