CN106295121B - 景观湖泊贝叶斯水质风险预测方法 - Google Patents

景观湖泊贝叶斯水质风险预测方法 Download PDF

Info

Publication number
CN106295121B
CN106295121B CN201610594828.3A CN201610594828A CN106295121B CN 106295121 B CN106295121 B CN 106295121B CN 201610594828 A CN201610594828 A CN 201610594828A CN 106295121 B CN106295121 B CN 106295121B
Authority
CN
China
Prior art keywords
landscape
water
water quality
risk
impoundments
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.)
Expired - Fee Related
Application number
CN201610594828.3A
Other languages
English (en)
Other versions
CN106295121A (zh
Inventor
成龙
彭森
杨立焜
吴卿
赵新华
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Tianjin University
Original Assignee
Tianjin University
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Tianjin University filed Critical Tianjin University
Priority to CN201610594828.3A priority Critical patent/CN106295121B/zh
Publication of CN106295121A publication Critical patent/CN106295121A/zh
Application granted granted Critical
Publication of CN106295121B publication Critical patent/CN106295121B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
    • G06N7/00—Computing arrangements based on specific mathematical models
    • G06N7/01—Probabilistic graphical models, e.g. probabilistic networks
    • G—PHYSICS
    • G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
    • G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
    • G16C20/70—Machine learning, data mining or chemometrics
    • 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
    • Y02A20/00—Water conservation; Efficient water supply; Efficient water use
    • Y02A20/152—Water filtration

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Software Systems (AREA)
  • Artificial Intelligence (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Evolutionary Computation (AREA)
  • Computing Systems (AREA)
  • General Health & Medical Sciences (AREA)
  • Probability & Statistics with Applications (AREA)
  • Crystallography & Structural Chemistry (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Medical Informatics (AREA)
  • Health & Medical Sciences (AREA)
  • Databases & Information Systems (AREA)
  • Chemical & Material Sciences (AREA)
  • Algebra (AREA)
  • Computational Mathematics (AREA)
  • Computer Vision & Pattern Recognition (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Mathematical Physics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

本发明的景观湖泊富营养化风险预测方法，通过联合使用动态模型方程和贝叶斯统计方法建立了景观湖泊富营养化风险概率预测模型。经验证，该模型预测精度良好(精确度高于80％)，从而使研究人员在面临景观湖泊入水水质变化时，可以及时、准确地预测景观湖水富营养化风险。通过比较不同入水下景观湖水质预测分布情况，指出氮和磷营养盐浓度是影响景观湖泊富营养化风险的主要因素，结合风险计算，初步拟定了不同入水景观湖富营养风险等级的应对措施，使景观湖泊在面临富营养风险时，可以采取合理措施保证景观湖泊水质，降低富营养化引发的藻类爆发。本发明对景观湖泊水质管理也可以提供一定参考。

Description

景观湖泊贝叶斯水质风险预测方法
技术领域
本发明涉及一种景观湖泊富营养化风险预测方法建立，特别是涉及一种可以为景观湖泊入水水质提供依据，使用贝叶斯方法对景观湖泊富营养化风险进行预测的方法。
背景技术
富营养化通常作为大量污染物排入水体导致有机物质加速产生的过程(尤其藻类)富营养化在全世界江河湖海中越来越严重并导致水环境恶化，尤其景观湖水体尤为严重。为了防止富营养化导致的浮游植物大爆发，对景观湖水体的评估和管理非常重要。在景观湖水体评估管理中机理模型用于解释生态过程，预测景观湖水环境因外来营养物质引起的改变，最后为政策制定做参考。然而传统模型运算忽略了模型参数的不确定性，将模型参数作为固定值进行校验，传统模型进行预测容易导致精度不高的现象发生。所以为了模型预测运算实现更加好的精度，必须考虑模型结构和参数中的不确定性因素影响。这样的模型才能容易的被决策制定者和政策规划人员利用。
发明内容
为了解决现有技术中存在的问题，本发明提供一种景观湖泊贝叶斯水质风险预测方法，克服现有技术中景观湖水体评估管理中机理模型预测不准确的问题。
本发明的技术方案是：一种景观湖泊富营养化风险贝叶斯预测模型方法，包括如下步骤：
(1)进行数据监测和采集；
(2)通过通讯系统将数据传输到服务器上；
(3)通过使用景观湖泊富营养化动态模型和贝叶斯统计方法建立景观湖泊富营养化风险不确定性预测模型；
(4)读取服务器数据库中存储的相应监测数据，并将监测数据输入水质模型，对景观湖水体富营养风险进行预测、分析和决策支持。
对于进水水质有周期性波动的景观湖，在建立模型时，选取的样本监测数据要具备周期性变化的特征，即选取一个完整水质变化周期内的样本监测数据进行建模，以保证模型模拟结果的稳定可靠。
所述步骤(1)须监测和采集的指标包括：氨氮、硝态氮、磷酸盐、叶绿素和溶解氧浓度。
所述的步骤(4)对景观湖水体富营养风险进行预测、分析和决策支持包括以下几步骤：
(1)根据不同入水浓度分布下监测到的进水氨氮、硝态氮、磷酸盐、叶绿素和溶解氧浓度来判断对景观湖泊富营养风险概率；将入水浓度情况，带入景观湖泊富营养化风险不确定性预测模型，预测未来一段时间内水质变化情况，获取模拟均值和上下95％置信值；
(2)以《地表水环境质量标准》水质指标限值为依据，将估计的整个预测值在标准内各类水质分布的概率作为预警强度划分的依据，通过使用指标限值的量△作为预警强度划分的依据；
(3)将95％上下置信值和均值分别作预警强度计算，获得一个预警强度范围。
所述步骤(2)指标限值的量△是指：
若区间估计的上限值大于相应的指标限值，而区间估计的下限值小于指标限值，则：
其中，l1和l2分别为标准内各类水质上下限值，l为模型预测出来相应水质指标浓度；α为参数：当预测值为I类水体数值为0，II类水体数值为20，III类水体数值为40，IV类水体数值为60，V类水体数值为80。
若区间估计的下限值大于V类水质指标限值，则：△＝100％。
本发明的有益效果为：本发明的景观湖泊富营养化风险预测方法和系统，通过联合使用动态模型方程和贝叶斯统计方法建立了景观湖泊富营养化风险概率预测模型。经验证，该模型预测精度良好(精确度高于80％)，从而使研究人员在面临景观湖泊入水水质变化时，可以及时、准确地预测景观湖水富营养化风险。通过比较不同入水下景观湖水质预测分布情况，指出氮和磷营养盐浓度是影响景观湖泊富营养化风险的主要因素，结合风险计算，拟定了不同入水景观湖富营养风险等级的应对措施，使景观湖泊在面临富营养风险时，可以采取合理措施保证景观湖泊水质，降低富营养化引发的藻类爆发。本发明可以为景观湖泊水质管理提供参考。
附图说明
图1是本发明的景观水水质在线预警方法的原理结构示意图；
图2是本发明的景观水水质预测预警部分的运行流程图。
具体实施方式
下面结合实施例和附图对本发明的景观湖富营养化风险预测方法做出详细说明。
本发明的景观湖富营养化风险预测方法，包括如下步骤：
1)进行数据监测和采集，是由若干个遥测远传终端组成，用于采集景观湖进、出水监测点的相关理化指标数据。须监测和采集的指标包括：氨氮、硝态氮、磷酸盐、叶绿素和溶解氧浓度；
2)通过通讯系统将数据传输到服务器上；该数据库中按时序存储了在线监测到的理化指标数据以及工艺基本情况。建立服务器用于控制数据监测和采集系统对景观湖进、出水监测点理化指标进行循环数据采集，将采集数据存储到数据库中并进行操作与分析；服务器可随时操作数据库中的数据，如进行存储/读取，以图表方式显示数据等。
2.景观湖泊贝叶斯水质风险预测方法，包括：
(1)建立景观湖泊水质模型主要水质模型，主要模型方程见下：
完整模型公式见说明书附录1。
(2)所述的对景观湖水质进行贝叶斯预测需要通过使用贝叶斯统计方法使用马尔科夫链-蒙托卡洛斯(MCMC)算法实现。马尔科夫链-蒙托卡洛斯(MCMC)公式为：
式中m为变量数，n是实测值数量。∑δj为模型误差，yj和fj(θ,x,y0)分别为模型计算值和模型模拟值；p(θ)是前验参数密度分布，p(σ2)基于基于前验密度的反伽马分布。将(6)和(7)联立，获得方程(8)。
(3)使用水质模型方程和贝叶斯统计方法对模型进行拟合，结合方程(8)，将入水氨氮、硝态氮、磷酸盐、溶解氧浓度和景观湖水氨氮、硝态氮、磷酸盐、叶绿素、溶解氧浓度输入到模型。读取服务器数据库中存储的相应监测数据，并将监测数据输入水质模型，对景观湖水体富营养风险进行预测；对于进水水质有周期性波动的景观湖，在建模时，选取的样本监测数据要具备周期性变化的特征，即选取一个完整水质变化周期内的样本监测数据进行建模，以保证模型模拟结果的稳定可靠。在实际使用中，本发明中所构建的富营养化风险预测模型可以每年更新一次建模数据(即每年重新建模一次)，而如果景观湖泊的年际水质相对稳定，模型更新周期可以设定得更长。最后通过模型计算，得到后验参数分布，获取模型求解结果(平均值和上下95％置信线)，以及检查模型模拟准确度。
3.根据不同入水情况下浓度的分布情况，来判断景观湖水质风险预警；
(1)所述的根据不同入水浓度分布下监测到的进水总氮、总磷、氨氮、硝态氮、磷酸盐、溶解氧浓度来判断对景观湖泊富营养风险概率。将入水浓度情况，带入1中所求解的模型，预测未来一段时间内水质变化情况，获取模拟均值和上下95％置信值。
(2)以获得模拟数据为基础，由《地表水环境质量标准》(GB3838-2002)可知，为了满足景观水体的观赏功能和生态服务功能，要求景观水水质至少能够满足V类水质的标准。所以，以GB3838-2002水质指标限值为依据，将估计的整个预测值在标准内各类水质分布的概率作为预警强度划分的依据，通过使用指标限值的量△作为预警强度划分的依据。若区间估计的上限值大于相应的指标限值，而区间估计的下限值小于指标限值，则：
其中，l1和l2分别为标准内各类水质上下限值，l为模型预测出来相应水质指标浓度；α为参数：当预测值为I类水体数值为0，II类水体数值为20，III类水体数值为40，IV类水体数值为60，V类水体数值为80。
若区间估计的下限值大于V类水质指标限值，则：Δ＝100％；
所以，预警依据和预警强度的对照表如下表1：
表1景观水体水质预警依据与预警强度对照表
(3)将95％上下置信值和均值分别作预警强度计算，可以获得一个预警强度范围。针对不同的预警强度范围，水质管理决策者就应该采取相应的处理措施。例如，对于预警强度范围为无警，这说明景观水的水质能够满足I类水体的标准；对于强度范围为轻警，可以对景观水进行日常的监测，但需要引起注意；对于强度范围为中警，就需要对景观水体加强监测，并引起较大的注意；对于范围内有重警，有很大可能景观水体已经无法满足水质标准，需要进行更为严格和迅速的水质治理。
尽管上面结合附图对本发明进行了描述，但是本发明并不局限于上述的具体实施方式，上述的具体实施方式仅仅是示意性的，并不是限制性的，本领域的普通技术人员在本发明的启示下，在不脱离本发明宗旨和权利要求所保护的范围情况下，还可以做出很多形式，这些均属于本发明的保护范围之内。
附录1
表2模型数学公式表达

Claims (3)

1.一种景观湖泊富营养化风险贝叶斯模型预测方法，其特征在于，包括如下步骤：
(1)进行数据监测和采集；
(2)通过通讯系统将数据传输到服务器上；
(3)通过联合景观湖泊富营养化动态模型方程和贝叶斯方法建立景观湖泊富营养化风险不确定性预测模型；
(4)读取服务器数据库中存储的相应监测数据，并将监测数据输入水质模型，对景观湖水体富营养风险进行预测、分析和决策支持，包括以下几步骤：
(a)根据不同入水浓度分布下监测到的进水氨氮、硝态氮、磷酸盐、叶绿素和溶解氧浓度来判断对景观湖泊富营养风险概率；将入水浓度情况，带入景观湖泊富营养化风险不确定性预测模型，预测未来一段时间内水质变化情况，获取模拟均值和上下95％置信值；
(b)以《地表水环境质量标准》水质指标限值为依据，将估计的整个预测值在标准内各类水质分布的概率作为预警强度划分的依据，通过使用指标限值的量△作为预警强度划分的依据；其中指标限值的量△是指：
若区间估计的上限值大于相应的指标限值，而区间估计的下限值小于指标限值，则：
其中，l1和l2分别为标准内各类水质上下限值，l为模型预测出来相应水质指标浓度；α为参数：当预测值为I类水体数值为0，II类水体数值为20，III类水体数值为40，IV类水体数值为60，V类水体数值为80；
若区间估计的下限值大于V类水质指标限值，则：△＝100％；
(c)将95％上下置信值和均值分别作预警强度计算，获得一个预警强度范围。
2.根据权利要求1所述景观湖泊富营养化风险贝叶斯模型预测方法，其特征在于，对于进水水质有周期性波动的景观湖，在建立模型时，选取的样本监测数据要具备周期性变化的特征，即选取一个完整水质变化周期内的样本监测数据进行建模，以保证模型模拟结果的稳定可靠。
3.根据权利要求1所述的景观湖泊富营养化风险贝叶斯模型预测方法，其特征在于，所述步骤(1)须监测和采集的指标包括：氨氮、硝态氮、磷酸盐、叶绿素和溶解氧浓度。
CN201610594828.3A 2016-07-21 2016-07-21 景观湖泊贝叶斯水质风险预测方法 Expired - Fee Related CN106295121B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610594828.3A CN106295121B (zh) 2016-07-21 2016-07-21 景观湖泊贝叶斯水质风险预测方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610594828.3A CN106295121B (zh) 2016-07-21 2016-07-21 景观湖泊贝叶斯水质风险预测方法

Publications (2)

Publication Number Publication Date
CN106295121A CN106295121A (zh) 2017-01-04
CN106295121B true CN106295121B (zh) 2019-01-18

Family

ID=57652861

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610594828.3A Expired - Fee Related CN106295121B (zh) 2016-07-21 2016-07-21 景观湖泊贝叶斯水质风险预测方法

Country Status (1)

Country Link
CN (1) CN106295121B (zh)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108152468A (zh) * 2018-01-26 2018-06-12 中冶华天工程技术有限公司 一种基于大数据的水质健康评估方法和装置
CN108537439A (zh) * 2018-04-09 2018-09-14 中国科学院遥感与数字地球研究所 一种流域多尺度景观格局与水质指标关系研究方法
CN108830425A (zh) * 2018-06-26 2018-11-16 北京高能时代环境技术股份有限公司 水库水质预测系统及方法
CN109460869B (zh) * 2018-11-06 2021-03-12 北京英视睿达科技有限公司 水环境预警方法和装置
CN110070144A (zh) * 2019-04-30 2019-07-30 云南师范大学 一种湖泊水质预测方法及系统
CN110379463B (zh) * 2019-06-05 2021-10-26 山东大学 基于机器学习的海洋藻类成因分析及浓度预测方法和系统
CN116400035B (zh) * 2023-04-11 2023-12-15 江苏省秦淮河水利工程管理处 一种湖泊生态监测管理方法与系统

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102807301A (zh) * 2012-07-31 2012-12-05 天津大学 再生水厂出水水质在线监测与实时预测系统及其控制方法
CN104155423A (zh) * 2014-08-08 2014-11-19 天津大学 一种变权组合的景观水水质预警方法

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102807301A (zh) * 2012-07-31 2012-12-05 天津大学 再生水厂出水水质在线监测与实时预测系统及其控制方法
CN104155423A (zh) * 2014-08-08 2014-11-19 天津大学 一种变权组合的景观水水质预警方法

Non-Patent Citations (7)

* Cited by examiner, † Cited by third party
Title
A case study of an enhanced eutrophication model with stoichiometric zooplankton growth sub-model calibrated by Bayesian method;Likun Yang et al;《Environmental Science & Pollution Research》;20160118;第23卷(第9期);8398-8409
Eutrophication risk assessment using Bayesian calibration of process-based models Application to a mesotrophic lake;George B.Arhonditsis et al;《Ecological Modelling》;20070712;215-229
Integration of Bayesian analysis for eutrophication prediction and assessment in a landscape lake;Likun Yang et al;《《Environmental Monitoring & Assessment》;20141203;第187卷(第1期);1-19
Nutrient dynamics and the eutrophication of shallow lakes Kasumigaura (Japan), Donghu (PR China), and Okeechobee (USA);K.EHavens et al;《Environment Pollution》;20110228;第111卷(第2期);263-272
Predicting the Frequency of Water Quality Standard Violations Using Bayesian Calibration of Eutrophication Models;Weitao Zhang et al;《Journal of Great Lakes Research》;20081231;第34卷(第4期);698-720
三峡库区环境水质监测管理信息系统的设计和研究;陈益;《中国优秀硕士学位论文全文数据库（电子期刊）信息科技辑》;20110415;第2011年卷(第4期);第4-5章
基于区间型贝叶斯的湖泊水质评价模型;祝慧娜;《监测与评价》;20160115;第33卷(第2期);130-134

Also Published As

Publication number Publication date
CN106295121A (zh) 2017-01-04

Similar Documents

Publication Publication Date Title
CN106295121B (zh) 景观湖泊贝叶斯水质风险预测方法
CN108171259A (zh) 一种基于NB-IoT的智慧末梢水质监测系统及方法
CN107607897B (zh) 一种电压监测仪在线校验及预测装置和方法
CN115774953B (zh) 一种基于数据处理的污染时空风险监管评估系统及方法
CN113610381B (zh) 一种基于5g网络的水质远程实时监测系统
CN106354940B (zh) 基于wasp水质模型不确定性输入的景观水水质模拟与预警方法
CN111523727B (zh) 基于不确定过程的考虑恢复效应的电池剩余寿命预测方法
JP6038011B2 (ja) ダム下流河川流量予測システム
CN103226737B (zh) 基于卡尔曼滤波与灰色预测的化工异常工况趋势预测方法
CN102509020B (zh) 基于传感器网络的复杂环境下多目标信息融合方法
CN110244387A (zh) 一种基于大气可降水量预测降雨天气的方法、装置、设备以及存储介质
US11480956B2 (en) Computing an explainable event horizon estimate
Badyalina et al. Streamflow estimation at ungauged basin using modified group method of data handling
RU2480825C2 (ru) Способ прогнозирования времени наступления и уровня паводков
CN116187204B (zh) 一种水域遥感反演水环境污染的风险决策方法
CN113962253A (zh) 基于深度小波极限学习机的轴承剩余寿命预测方法及系统
CN109215821A (zh) 核电站蒸发器主给水流量冗余测量通道比较方法及系统
Chenyu et al. Coupling of adjoint-based Markov/CCMT predictive analytics with data assimilation for real-time risk scenario forecasting of industrial digital process control systems
CN102122004A (zh) 铁路防灾风预测方法
CN110517731A (zh) 基因检测质量监控数据处理方法和系统
CN113887058B (zh) 一种考虑距海岸线距离和风速影响的氯离子沉积率预测方法
CN107679020A (zh) 一种直接预测海冰初、终冰日的方法
CN114624411A (zh) 一种水源地取水口附近水域的水质预警方法
Kabi et al. Leveraging IoT and Machine Learning for Improved Monitoring of Water Resources-A Case Study of the Upper Ewaso Nyiro River
Hamidi et al. Prediction of MEMS-based INS error using interval type-2 fuzzy logic system in INS/GPS integration

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20190118

Termination date: 20210721

CF01 Termination of patent right due to non-payment of annual fee