本发明涉及水域纳污技术领域，具体涉及一种基于最小二乘法的水库纳污能力预警方法。
背景技术：
水域纳污能力分析是水污染总量控制的基础，其定量核算对保护水环境具有重要的意义。现行水库的纳污能力即水域污染纳污能力污染负荷计算方法采用《水域纳污能力计算规程》gb/t25173—2001，存在以下问题：在水库纳污能力计算之前，需要进行水文资料、水质资料、排污口资料、水下地形资料等内容的调查，耗时长，且数据容易缺失和存在误差；模型计算法中的均匀混合模型、非均匀混合模型和分成模型需要大量数据的检验才能确定参数，存在一定的偶然；模型的求解及计算采用条件简化和假定的方式，存在主观判断，最终的结果无法客观的体现出水库的纳污能力。在不能准确计算水库的实时纳污能力的条件下，需要对水库水质进行实时监测并预测，以供管理者辅助决策，及时采取措施，以防水库水质恶劣变化。
技术实现要素：
为解决上述现有技术存在的问题，本发明提出一种基于最小二乘法的水库纳污能力预警方法，通过最小二乘法对未来12小时的水库水质参数进行滚动预测，及时报警提示。
为达到上述目的，本发明的技术方案是：一种基于最小二乘法的水库纳污能力预警方法，包括以下步骤：
步骤s1：在水库选取若干个水质监测点；
步骤s2：在水质监测点对高、中、低三个水层的水质进行监测，利用水质传感器在线监测水质参数，所述水质参数包括：温度、ph值、浊度、含氧量、电导率、高锰酸盐、氨氮、总磷总氮、叶绿素、蓝绿藻；
步骤s3：计算水质监测点所在的区域水质参数和整个水库的水质参数：
smn＝(dmn+emn+fmn)/3
其中，m＝1,2,...,m,m为水质监测点数量，n＝1,2,...n,n为监测的水质参数的种类数，smn为第m个水质监测点所在区域的第n种水质参数数据，dmn、emn、fmn分别为第m个水质监测点高、中、低层的第n种水质参数数据，sn为水库的第n种水质参数数据；
步骤s4：根据水库的历史水质参数数据，利用最小二乘法曲线拟合计算未来12小时的水质参数；
步骤s5：根据水库水质等级要求设定每种水质参数预警值，当预测时间段内出现预警值则报警提示。
进一步地，所述步骤s1按如下方法选取水质监测点：将水库的水面拟化成若干个面积相等的多边形，利用多边形最大内切圆算法得到圆心，所述圆心所在地即选取为水库的水质监测点。
进一步地，所述步骤s4具体包括：
步骤s41：将监测到的预测时刻前240小时的水库水质参数按时间和水质参数数值在直角坐标系中形成若干点pi(xi,yi)，x轴为时间，y轴为水质参数数值，得到y＝f(x)的曲线,i＝1,2,...i,i为历史数据个数；
步骤s42：计算y＝f(x)的拟合曲线
设拟合曲线的多项式为：
y＝a0+a1x+…+akxk
令点pi(xi,yi)到拟合曲线的偏差的平方和最小：
通过多项式的分解计算得到a0,a1,...,ak，从而得到拟合曲线
与现有技术相比，本发明具有以下有益效果：
(1)数据来源于水质实时在线监测，无需完整的水文资料或水下地形资料等；
(2)对监测数据进行数据优化后进行函数拟合，数据客观性和实时性较高。
附图说明
图1是本发明一种基于最小二乘法的水库纳污能力预警方法的流程示意图。
具体实施方式
下面结合附图及实施例对本发明做进一步说明。
如图1所示，一种基于最小二乘法的水库纳污能力预警方法，包括以下步骤：
步骤s1：将水库的水面拟化成10个面积相等的多边形，利用多边形最大内切圆算法得到圆心，所述圆心所在地即选取为水库的水质监测点。
步骤s2：在水质监测点对高、中、低三个水层的水质进行监测，利用水质传感器在线监测水质参数，所述水质参数包括：温度、ph值、浊度、含氧量、电导率、高锰酸盐、氨氮、总磷总氮、叶绿素、蓝绿藻；
步骤s3：计算水质监测点所在的区域水质参数和整个水库的水质参数：
smn＝(dmn+emn+fmn)/3
其中，m＝1,2,...,m,m,为水质监测点数量，m＝10，n＝1,2,...n,n为监测的水质参数的种类数，smn为第m个水质监测点所在区域的第n种水质参数数据，dmn、emn、fmn分别为第m个水质监测点高、中、低层的第n种水质参数数据，sn为水库的第n种水质参数数据；
步骤s4：根据水库的历史水质参数数据，利用最小二乘法曲线拟合计算未来12小时的水质参数；
步骤s4具体包括：
步骤s41：将监测到的预测时刻前240小时的水库水质参数按时间和水质参数数值在直角坐标系中形成若干点pi(xi,yi)，x轴为时间，y轴为水质参数数值，得到y＝f(x)的曲线,i＝1,2,...i,i为历史数据个数；在本实施例中，以1小时为时间间隔，每种水质参数有240个数据；还可以缩短时间间隔，选取更多的历史数据；
步骤s42：计算y＝f(x)的拟合曲线
设拟合曲线的多项式为：
y＝a0+a1x+…+akxk
令点pi(xi,yi)到拟合曲线的偏差的平方和最小：
通过多项式的分解计算得到a0,a1,...,ak，从而得到拟合曲线
步骤s5：按照水库水质等级要求，根据《中华人民共和国地表水环境质量标准》，对照《地表水环境质量标准基本项目标准限值》设定每种水质参数预警值，当预测时间段内出现预警值则报警提示。
以上所述仅为本发明的较佳实施例，凡依本发明申请专利范围所做的均等变化与修饰，皆应属本发明的涵盖范围。