基于果蝇优化算法的水质机理建模及水质预测方法技术

本发明专利技术公开了一种基于果蝇优化算法的水质机理建模及水质预测方法，属于环境工程技术领域。本发明专利技术首先结合水质演化机理，提出了一种基于果蝇优化算法的湖库水质机理建模及水质预测方法。在已知水质测量数据的基础上，利用果蝇优化算法估计水质机理模型的未知参数。在此基础上，利用蒙特卡洛仿真预测水质演化过程，获得水质指标未来时刻取值的概率分布，实现水质预测。本发明专利技术能建立准确的水质机理模型，基于果蝇优化算法的水质机理建模方法能够准确估计水质机理模型的未知参数，与已有方法相比提高了优化精度和速度；水质预测方法能有效实现水质预测，且考虑更加全面、准确，克服了单值预测结果带来的偶然性。

【技术实现步骤摘要】
基于果蝇优化算法的水质机理建模及水质预测方法
本专利技术设计一种基于果蝇优化算法的水质机理建模及水质预测方法，属于环境工程

技术介绍
水是我们星球上生命的基础，但这种宝贵的资源越来越受到威胁。国际湖泊环境委员会与联合国环境规划署合作开展了题为“世界湖泊状况调查”的项目。该项目收集了全球217个湖泊详细数据。通过这个项目，可以确定六个主要的环境问题都对水质有重大影响，富营养化问题就是其中之一。目前，全球不同地区的湖泊富营养化百分比分别为：亚洲54％，欧洲53％，非洲28％；北美洲48％；南美洲41％，南非60％以上。在中国，据《2016年中国环境公报》，108个监测营养状态的湖泊(水库)中，贫营养的湖泊(水库)有10个，中营养的湖泊(水库)有73个，轻度富营养的湖泊(水库)为20个，中度富营养的湖泊(水库)为5个。因此，水环境问题已刻不容缓。水质机理建模能表现水体内部环境变化，水质预测能有效的掌握水质情况，二者结合能预防大规模水污染事件的发生，对水污染防治具有现实意义，这是解决水环境问题的重要思路。已有的水质建模方法主要分为两类：基于机理的水质建模方法和数据驱动的水质建模方法。(1)基于机理的水质建模方法。这类方法主要是通过分析水体内部的机理演化过程来建立水质模型。已有的经典水质机理模型包括：MIKE、WASP和EFDC等。MIKE是丹麦水动力研究所(DHI)开发，用于模拟河网、河口、滩涂等地区的水质情况，同时也能模拟水体环境中垂向变化常被忽略的湖泊、河口、海岸地区等水质情况，它拥有一维动态模型和二维动态模型，其中二维动态模型能处理三维空间。WASP(Thewaterqualityanalysissimulationprogram)是美国环保局(EPA)开发并推荐使用的用于各类水体水质模拟与分析的软件；它可用于模拟水动力学、河流一维不稳定流、湖泊和河口三维不稳定流、常规污染物和有毒污染物在水中的迁移和转化规律；其基本原理是平移—扩散方程。EFDC(TheEnvironmentalFluidDynamicsCode)模型是涵盖一维到三维地表水水质的数学模型；它可实现河流、湖泊、水库、湿地系统、河口和海洋等水体的水动力学、水质模拟、污染物迁移和泥沙输移等，是一个多参数有限差分模型。(2)数据驱动的水质建模方法。这类方法利用数据分析模型和工具对水质指标的实测数据进行建模。常用的数据驱动模型包括：人工神经网络、时间序列模型和支持向量机等，该类模型均是利用不同时刻的水质指标观测数据建立表征水质演化过程的数据驱动模型。两类方法各具有一定优势，然而基于数据驱动的水质建模效果取决于数据量的数量和质量，且无法合理解释水质指标演化过程。目前，已有的经典水质机理模型通常包含大量未知参数，且仅给定经验取值范围，在此基础上预测水质变化难以满足精度要求。
技术实现思路
本专利技术为了准确计算水质机理模型中未知参数，提出基于果蝇优化算法(Fruitflyoptimizationalgorithm，FOA)的水质机理建模方法，同时，本专利技术为了弥补单值预测方法带来的不足，结合蒙特卡洛仿真提出一种水质预测方法，为水质预测的不精确性问题提供一种新思路和可行解决方案。本专利技术提供一种基于果蝇优化算法的水质机理建模及水质预测方法，该方法具体步骤如下：步骤一：建立水质机理模型。基于水质指标，DO、Phyt、BOD、OP、PO4、ON、NH3-N和NO3-N，建立水质机理模型；DO、Phyt、BOD、OP、PO4、ON、NH3-N和NO3-N分别为溶解氧、浮游植物、碳生化需氧量、有机磷、正磷酸盐、有机氮、氨氮和硝酸盐氮；步骤二：基于果蝇优化算法的水质机理建模方法。果蝇优化算法是一种基于果蝇觅食行为来寻求全局优化的新方法，属于仿生智能优化算法。果蝇是一种感官知觉优于其他物种的昆虫，尤其是嗅觉和视觉。果蝇的嗅觉器官能很好的搜集漂浮在空气中的各种气味，甚至能够嗅到40公里以外的食物源。然后，飞到食物位置附近后亦可使用敏锐的视觉发现食物和同伴聚集的位置，并且向该方向飞去。果蝇优化算法在寻优过程中具有一定的随机性，为了确保果蝇群体向着正确的方向飞行，该算法引入了气味浓度判定值和气味浓度判定函数。果蝇优化算法具有计算过程简单、全局收敛性强、执行的时间短等优点，可应用于求解最优解，在求解数学函数极值、数据挖掘、广义回归神经网络参数优化等领域已有初步成果。步骤三：基于蒙特卡洛仿真的河湖水质预测。蒙特卡洛仿真的基本原理是当问题或对象本身具有概率特征时，可以用计算机模拟的方法产生抽样结果，根据抽样结果计算统计量或者参数的值；随着模拟次数的增多，可以通过对各次统计量或参数的估计值求平均的方法得到稳定结论。在基于蒙特卡洛仿真的预测过程中，水质指标演化过程会被多次模拟。对每一次模拟，每一个水质指标在每一个预测时间点的数值均不相同，这些数值同时具有规律性和随机性。规律性体现在预测值是基于确定的水质动力学模型、水质指标和模型参数初值而产生的；随机性体现在每一个水质指标在每一个预测时间点受噪声影响而得到的预测值不尽相同。当模拟次数充分大时，结合概率统计规律，可以获得各水质指标在不同时间点上取值的概率分布，实现水质预测。本专利技术的优点在于：(1)本专利技术利用果蝇优化算法估计水质机理模型的未知参数。克服了已有的经典水质机理模型通常包含大量未知参数，且仅给定经验取值范围的问题，同时与粒子群优化算法和遗传算法对比，果蝇优化算法的寻优速度和精度均优于粒子群优化算法和遗传算法。通过该果蝇优化算法方法获得水质机理建模准确度更高。(2)本专利技术利用蒙特卡洛仿真预测8个水质指标的演化过程，获得8个水质指标取值的概率分布，实现水质预测。弥补了单一定值的水质预测的不精确性，降低了预测的偶然性。附图说明图1是本专利技术基于果蝇优化算法(FOA)的水质机理建模的8个水质指标演化曲线、基于粒子群优化算法(PSO)的水质机理建模的8个水质指标演化曲线、基于遗传算法(GA)的8个水质指标机理建模的水质演化曲线和8个水质指标真实值曲线图。图2是本专利技术本专利技术基于果蝇优化算法(FOA)的水质机理建模的水质演化误差、基于粒子群优化算法(PSO)的水质机理建模的水质演化误差和基于遗传算法(GA)的8个水质指标机理建模的水质演化误差减小曲线图。图3是本专利技术基于蒙特卡洛仿真的水质预测方法的8个水质指标演化曲线图。图4是本专利技术不同时刻下的8个水质指标的概率密度函数曲线图。具体实施方式下面将结合附图和实施例对本专利技术作进一步的详细说明。本专利技术是一种基于果蝇优化算法的水质机理建模及水质预测方法，首先结合水质演化机理，提出了一种湖库水质机理建模及水质预测方法；在已知水质测量数据的基础上，利用果蝇优化算法估计水质机理模型的未知参数；在此基础上，利用蒙特卡洛仿真预测水质演化过程，获得水质指标未来时刻取值的概率分布，实现水质预测。本专利技术提供的所述基于果蝇优化算法的水质机理建模及水质预测方法，通过如下步骤实现：步骤一：建立水质机理模型。水体生态系统主要分为以下四个过程：溶解氧平衡过程、浮游植物动力学过程、磷循环和氮循环。水体生态系统主要由八个水质指标描述其演化过程：溶解氧(DO)、浮游植物(Phyt)、碳生化需氧量(BOD)、有机磷(OP)、本文档来自技高网...

【技术保护点】
1.基于果蝇优化算法的水质机理建模及水质预测方法，其特征在于：具体步骤如下，步骤一：建立水质机理模型；基于水质指标，DO、Phyt、BOD、OP、PO4、ON、NH3‑N和NO3‑N，建立水质机理模型；DO、Phyt、BOD、OP、PO4、ON、NH3‑N和NO3‑N分别为溶解氧、浮游植物、碳生化需氧量、有机磷、正磷酸盐、有机氮、氨氮和硝酸盐氮；水质机理模型的数学模型为上述8个水质指标方程组成的一组非线性连续微分方程组：

【技术特征摘要】
1.基于果蝇优化算法的水质机理建模及水质预测方法，其特征在于：具体步骤如下，步骤一：建立水质机理模型；基于水质指标，DO、Phyt、BOD、OP、PO4、ON、NH3-N和NO3-N，建立水质机理模型；DO、Phyt、BOD、OP、PO4、ON、NH3-N和NO3-N分别为溶解氧、浮游植物、碳生化需氧量、有机磷、正磷酸盐、有机氮、氨氮和硝酸盐氮；水质机理模型的数学模型为上述8个水质指标方程组成的一组非线性连续微分方程组：其中，分别表示水质指标DO、Phyt、BOD、OP、PO4、ON、NH3-N和NO3-N的浓度，分别表示所述水质指标浓度的微分；Phyt的浓度由叶绿素a的浓度表征；参数{k1,k2,…,k13}为水质机理模型的未知的模型参数；令θ＝(k1,k2,k3,k4,k5,k6,k7,k8,k9,k10,k11,k12,k13)Τ且假设水质机理模型的模型参数是慢变的，则式(1)表示为：上式中，f(x,θ)为式(1)给出的水质机理模型；步骤二：基于果蝇优化算法的水质机理建模方法；首先，将式(2)写成离散化形式如下：上式中，xt为t时刻的水质指标，xt-1为t-1时刻的水质指标；θt为t时刻的模型参数，θt-1为t-1时刻的模型参数；水质指标的过程噪声ωx,t-1为t-1时刻的水质指标的过程噪声，模型参数的过程噪声ωθ,t-1为t-1时刻的模型参数的过程噪声，其中Σx、Σθ为相应的过程噪声协方差阵；h是预测步长，h≤1day，t为离散时间点；将果蝇优化算法寻优过程中的模型参数代入味道浓度函数中，获得该组模型参数的味道浓度，在果蝇种群中，对比每个果蝇的味道浓度，并找出最优的味道浓度和该果蝇位置，其余果蝇利用视觉飞往所述的最佳位置，最后在不断迭代过程中获得最优参数解；基于果蝇优化算法的水质机理建模方法包括初始化和迭代过程两个部分，具体步骤如下：(1)初始化：设定模型参数初值设定数据长度为ND,t1＝1,2,…,ND；计算目标函数初值,保存味道浓度Smell(0)为最佳味道浓度BestSmell；设定果蝇种群个体数为MaxGroup；设定果蝇种群个体初始位置均为设定最大迭代次数MaxGeneration和味道浓度阈值ObjSmell；(2)迭代过程：forgen＝1,2,…,MaxGeneration(2.1)种群寻优过程：forgro＝1,2,…,MaxGroup种群中每个果蝇利用嗅觉搜寻食物的随机距离与方向，每一个模型参...

【专利技术属性】
技术研发人员：赵峙尧，王小艺，周宇琴，许继平，王立，于家斌，
申请(专利权)人：北京工商大学，
类型：发明
国别省市：北京,11