一种湖泊富营养化状态发展趋势的预测方法技术

本发明专利技术公开了一种湖泊富营养化状态发展趋势的预测方法，属于环境保护和水资源可持续利用技术领域,其技术方案要点是首先选择与湖泊富营养化关系密切的工业废水排放量、生活污水排放量、农业面源污染负荷比和工业废水达标排放率、污水处理率等作为富营养化影响因素；然后建立富营养化状态预测S型曲线指数公式；进而结合富营养化指标监测数据和影响因素数据，应用粒子群算法对公式中的参数进行优化，得出适用于富营养化发展趋势的预测指数公式。达到了能适用于不同地区的湖泊富营养化状态发展趋势预测的技术效果,有效的保护湖泊生态系统。

【技术实现步骤摘要】
一种湖泊富营养化状态发展趋势的预测方法
本专利技术属于环境保护和水资源可持续利用
，更具体地说，它涉及一种湖泊富营养化状态发展趋势的预测方法。
技术介绍
湖泊富营养化是当今全球范围内最严重的水质问题之一。国务院在京发布《2011年中国环境状况公报》，公报中指出我国湖泊(水库)富营养化问题突出，监测的26个湖泊(水库)中，处于富营养化状态的湖泊(水库)占53.8％。环保部周生贤部长曾提出，国家在深入推进重点流域水污染防治的同时，已着手优先保护水质良好和生态脆弱的江河湖泊。近20年来，随着经济、社会的发展和人口增长，我国长江中、下游湖泊的富营养化进程已呈现逐渐加快的态势。不少湖泊水质恶化，可利用功能丧失，严重的制约了该地区人类社会和经济的发展。为了解湖泊富营养化状况和未来发展趋势，需要在现有的富营养化指标监测数据、污染物排放量及影响湖泊富营养化的社会、经济等因素之间的相关分析基础上，建立描述湖泊富营养化状态与影响因素之间数量关系的预测模型。如公开号为CN106295121A的中国专利公开了一种景观湖泊贝叶斯水质风险预测方法,通过联合使用动态模型方程和贝叶斯统计方法建立了景观湖泊富营养化风险概率预测模型。湖泊富营养化预测模型通常分为动态模型和稳态模型。动态模型要求具有较多的调查和观测资料，计算复杂，且模型适应度差；稳态模型需要的资料较少，模型结构简单，计算结果是以年为时间尺度的平均状态，预测的结果与实际情况吻合较好，应用较广。其中具有代表性的稳态模型有Vollenweider模型、Dillon模型、Canfield模型及合田键模型。但是上述模型实际上只是单一的负荷模型，只能对磷的浓度进行预测，不能对水体中包含有Chla、TN和CODMn等多指标的湖泊富营养化状态进行综合预测。
技术实现思路
针对现有技术存在的不足，本专利技术的目的在于提供湖泊富营养化状态发展趋势的预测方法，其具有适用于不同地区的湖泊富营养化状态发展趋势预测的优势。为实现上述目的，本专利技术所提供的技术方案具体步骤如下：1.构建预测模型的评价因子采用五个参数作为预测湖泊富营养化状态发展趋势的影响因子：湖区工业废水排放量Q1，生活污水排放量Q2，农业面源有机污染负荷比η3，工业废水达标排放率η1和污水处理率η2等影响指标。2.湖泊富营养化的S型曲线预测指数公式建立如式(1)表示的湖泊富营养化的S型曲线预测指数公式：式中，式中，Qt1和ηt1分别为第t年工业废水排放量和工业废水达标排放率；Qt2和ηt2分别为第t年生活污水排放量和污水处理率；Q01和η01分别为所选择的“基准年”的工业废水排放量和工业废水达标排放率；Q02和η02分别为所选择的“基准年”的生活污水排放量和污水处理率；ηt3和η03分别为第t年和“基准年”农业面源污染负荷量所占比例。3.湖泊富营养化评价指标指数值计算(Ⅰ).单指标富营养化状态分指数值计算富营养化预测评价的指标j的营养状态分指数计算公式如式(3)所示。式中，cj为指标j的测定值；cjmax和cjmin分别为取指标j的最高营养级(极富)和最低营养级(极贫)标准限值。与富营养化关系密切的10项富营养化指标(以下简称指标)的cjmin和cjmax如表1所示。表110项湖泊富营养化评价指标的cjmin和cjmax及指数分级标准(Ⅱ).综合富营养化状态指数值计算富营养化预测评价的指标总体的营养状态指数计算公式如式(4)所示。式中，Wj为指标j的归一化权值，多数情况下，各指标可视作等权的，故Wj＝1/n；n为参与评价的评价指数数目。Ij为由式(3)计算得到的指标j的营养状态分指数。10项指标6级标准浓度限值(cjk)及由式(3)计算得到的单指标j的营养状态分指数Ij如表1所示。将表1中的10项指标的各级标准分指数值代入式(4)，可得出富营养状态的6级分级标准与综合指数值EI的对应关系。4.应用粒子群算法优化参数应用粒子群算法优化确定式(1)中的参数α和βi(i＝1,2,3)，设计优化目标函数：式中，EIt和EIt0分别为由式(1)计算得到的湖泊富营养化预测指数值和由幂函数加和型富营养化评价综合指数式(4)计算得到的湖泊富营养化评价综合指数值；N为建模所用年份总数。粒子群算法的迭代公式如式(6)、(7)所示。vid(t+1)＝w·vid(t)+c1·r1(pid-xid(t))+c2·r2(pgd-xgd(t))(6)xid(t+1)＝xid(t)+vid(t+1)(7)式中，w为惯性权值；c1和c2为加速系数；r1和r2为两个在[0，1]内变化的随机数。5.趋势预测应用优化后的富营养状态S型预测指数公式进行富营养化趋势预测。综上所述，本专利技术具有以下有益效果：1、本方法适应于多指标的富营养化预测评价，同时能提供尽可能优化的计算公式，有效的保护湖泊生态系统；2、本方法选择对富营养化影响密切的湖区工业废水排放量、生活污水排放量、农业有机面源污染负荷比及工业废水达标排放率和污水处理率作为富营养化预测的影响因子，建立湖泊富营养化预测指数公式，并应用粒子群算法优化参数，得到富营养化预测指数公式，并设定湖区未来的工业废水排放量、生活污水排放量、农业面源污染负荷比和工业废水达标排放率及污水处理率等影响因素的变化，通过预测指数公式,能够合理预测未来湖泊富营养状态的变化趋势。附图说明图1是1987～2005年湖区富营养化综合评价指数值的曲线图；图2是湖区富营养化状态变化趋势(2005～2035年)曲线图。具体实施方式下面以长江中下游地区某一湖泊为实例进一步说明本专利技术，该湖泊富营养化状态发展趋势的预测方法的具体步骤如下。1.构建预测模型的评价因子对长江中、下游湖泊，由于湖区大、工业生产发达、农业蓄水量达，人口密度大，采用对湖泊富营养化进程起重要影响的有湖区工业废水排放量Q1、生活污水排放量Q2、农业面源有机污染负荷比η3及工业废水达标排放率η1和污水处理率η2等影响指标。1999～2004年某湖区工业废水排放量Qt1、生活污水排放量Qt2、农业面源污染负荷比ηt3及工业废水达标排放率ηt1和污水处理率ηt2如表2所示。依据公式，计算出各项指标预测排放量。表2湖区工业废水、生活污水、面源污染负荷比、工业废水达标排放率、污水处理率及Rti(i＝1,2,3)2.湖泊富营养化的S型曲线预测指数公式以选取的指标作为湖泊富营养化预测模型的因子，建立如式(1)表示的湖泊富营养化的S型曲线预测指数公式：式中，式中，Qt1和ηt1分别为第t年工业废水排放量和工业废水达标排放率；Qt2和ηt2分别为第t年生活污水排放量和污水处理率；Q01和η01分别为所选择的“基准年”的工业废水排放量和工业废水达标排放率；Q02和η02分别为所选择的“基准年”的生活污水排放量和污水处理率；ηt3和η03分别为第t年和“基准年”农业面源污染负荷量所占比例。3.湖泊富营养化评价指标指数值计算依据式(3)计算得到的湖区3项富营养状态分指数值Ij(j＝1,2,3)及由式(4)计算得到的富营养状态综合指数值EIt如表3所示。表3湖区富营养状态分指数值Ij及综合指数值EIt(1999～2005)由式(4)计算得到的1987～2005年某湖区富营养化状态综合评价指数值如图1所示。4.应用粒子群算本文档来自技高网...

【技术保护点】
1.一种湖泊富营养化状态发展趋势的预测方法，其特征在于，包括如下步骤：(1)构建预测模型的评价因子；(2)建立湖泊富营养化的S型曲线预测指数公式；(3)湖泊富营养化评价指标指数值计算；(4)应用粒子群算法优化参数；(5)预测湖泊富营养化状态发展趋势。

【技术特征摘要】
1.一种湖泊富营养化状态发展趋势的预测方法，其特征在于，包括如下步骤：(1)构建预测模型的评价因子；(2)建立湖泊富营养化的S型曲线预测指数公式；(3)湖泊富营养化评价指标指数值计算；(4)应用粒子群算法优化参数；(5)预测湖泊富营养化状态发展趋势。2.根据权利要求1所述的湖泊富营养化状态发展趋势的预测方法，其特征在于，步骤(1)中采用五个参数作为预测湖泊富营养化状态发展趋势的影响因子：湖区工业废水排放量Q1，生活污水排放量Q2，农业面源有机污染负荷比η3，工业废水达标排放率η1和污水处理率η2。3.根据权利要求2所述的湖泊富营养化状态发展趋势的预测方法，其特征在于，步骤(2)中的湖泊富营养化的S型曲线预测指数公式如下：式中，式中，Qt1和ηt1分别为第t年工业废水排放量和工业废水达标排放率；Qt2和ηt2分别为第t年生活污水排放量和污水处理率；Q01和η01分别为所选择的“基准年”的工业废水排放量和工业废水达标排放率；Q02和η02分别为所选择的“基准年”的生活污水排放量和污水处理率；ηt3和η03分别为第t年和“基准年”农业面源污染负荷量所占比例。4.根据权利要求3所述的湖泊富营养化状态发展趋势的预测方法，其特征在于，步骤(3)包括(Ⅰ)单指标富营养化状态分指数值计算和(Ⅱ)综合富营养化状态指数值计算。5.根据权利要求4所述的湖泊富营养化状态发展趋势的预测方法，其特征在于，指标j的营养状态分指数计算公式如式(3)所示：式中，...

【专利技术属性】
技术研发人员：余春雪，杨志峰，李祚泳，尹心安，
申请(专利权)人：东莞理工学院，
类型：发明
国别省市：广东,44