CN108615090B - 基于蒙特卡洛仿真的河湖水质预测及富营养化风险评估方法 - Google Patents
基于蒙特卡洛仿真的河湖水质预测及富营养化风险评估方法 Download PDFInfo
- Publication number
- CN108615090B CN108615090B CN201810269206.2A CN201810269206A CN108615090B CN 108615090 B CN108615090 B CN 108615090B CN 201810269206 A CN201810269206 A CN 201810269206A CN 108615090 B CN108615090 B CN 108615090B
- Authority
- CN
- China
- Prior art keywords
- water quality
- prediction
- quality index
- index
- chl
- 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.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 155
- 238000000034 method Methods 0.000 title claims abstract description 65
- 238000012851 eutrophication Methods 0.000 title claims abstract description 27
- 238000000342 Monte Carlo simulation Methods 0.000 title claims abstract description 25
- 238000012502 risk assessment Methods 0.000 title claims abstract description 16
- 235000016709 nutrition Nutrition 0.000 claims abstract description 31
- 230000008569 process Effects 0.000 claims abstract description 26
- 230000035764 nutrition Effects 0.000 claims abstract description 24
- 238000004088 simulation Methods 0.000 claims description 19
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 claims description 13
- 125000001477 organic nitrogen group Chemical group 0.000 claims description 12
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims description 11
- 229910052760 oxygen Inorganic materials 0.000 claims description 11
- 239000001301 oxygen Substances 0.000 claims description 11
- OAICVXFJPJFONN-UHFFFAOYSA-N Phosphorus Chemical compound [P] OAICVXFJPJFONN-UHFFFAOYSA-N 0.000 claims description 10
- 229910052698 phosphorus Inorganic materials 0.000 claims description 10
- 239000011574 phosphorus Substances 0.000 claims description 10
- 239000002245 particle Substances 0.000 claims description 9
- MMDJDBSEMBIJBB-UHFFFAOYSA-N [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] Chemical compound [O-][N+]([O-])=O.[O-][N+]([O-])=O.[O-][N+]([O-])=O.[NH6+3] MMDJDBSEMBIJBB-UHFFFAOYSA-N 0.000 claims description 8
- XKMRRTOUMJRJIA-UHFFFAOYSA-N ammonia nh3 Chemical compound N.N XKMRRTOUMJRJIA-UHFFFAOYSA-N 0.000 claims description 8
- NBIIXXVUZAFLBC-UHFFFAOYSA-K phosphate Chemical compound [O-]P([O-])([O-])=O NBIIXXVUZAFLBC-UHFFFAOYSA-K 0.000 claims description 8
- 229910052757 nitrogen Inorganic materials 0.000 claims description 6
- 238000010606 normalization Methods 0.000 claims description 5
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 claims description 4
- 229910052799 carbon Inorganic materials 0.000 claims description 4
- ATNHDLDRLWWWCB-AENOIHSZSA-M chlorophyll a Chemical compound C1([C@@H](C(=O)OC)C(=O)C2=C3C)=C2N2C3=CC(C(CC)=C3C)=[N+]4C3=CC3=C(C=C)C(C)=C5N3[Mg-2]42[N+]2=C1[C@@H](CCC(=O)OC\C=C(/C)CCC[C@H](C)CCC[C@H](C)CCCC(C)C)[C@H](C)C2=C5 ATNHDLDRLWWWCB-AENOIHSZSA-M 0.000 claims description 4
- 239000011159 matrix material Substances 0.000 claims description 4
- 229930002875 chlorophyll Natural products 0.000 claims description 3
- 235000019804 chlorophyll Nutrition 0.000 claims description 3
- 125000004432 carbon atom Chemical group C* 0.000 claims description 2
- -1 phytoplankton Chemical compound 0.000 claims description 2
- 235000003784 poor nutrition Nutrition 0.000 claims description 2
- 230000007246 mechanism Effects 0.000 abstract description 7
- 230000007613 environmental effect Effects 0.000 abstract description 2
- 238000010276 construction Methods 0.000 description 3
- 238000010586 diagram Methods 0.000 description 3
- 239000013505 freshwater Substances 0.000 description 3
- 238000009792 diffusion process Methods 0.000 description 2
- 235000015097 nutrients Nutrition 0.000 description 2
- 150000003839 salts Chemical class 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000001228 trophic effect Effects 0.000 description 2
- 238000003911 water pollution Methods 0.000 description 2
- 238000012935 Averaging Methods 0.000 description 1
- 229910002651 NO3 Inorganic materials 0.000 description 1
- 229930002868 chlorophyll a Natural products 0.000 description 1
- 238000005094 computer simulation Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000003912 environmental pollution Methods 0.000 description 1
- 238000011156 evaluation Methods 0.000 description 1
- 238000013178 mathematical model Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000004172 nitrogen cycle Methods 0.000 description 1
- 235000003715 nutritional status Nutrition 0.000 description 1
- 235000008935 nutritious Nutrition 0.000 description 1
- 238000004175 phosphorus cycle Methods 0.000 description 1
- 230000002265 prevention Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 239000000779 smoke Substances 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0635—Risk analysis of enterprise or organisation activities
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- 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
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Operations Research (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Educational Administration (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种基于蒙特卡洛仿真的河湖水质预测及富营养化风险评估方法,属于环境工程技术领域。首先基于水质指标演化机理,建立8个水质指标的水质动力学模型。然后,利用蒙特卡罗仿真进行河湖水质预测,在已知水质动力学模型水质指标和模型参数的先验分布基础上,利用蒙特卡罗仿真预测水质指标的演化过程,获得未来时刻水质指标取值的概率分布,实现水质预测。最后,构造综合营养状态指数,结合水质指标预测结果,计算综合营养状态指数的概率分布和处于不同营养程度的概率,实现富营养化风险评估。本发明能够有效实现水质预测和富营养化分析,且考虑更加全面、准确,克服了单值预测结果带来的偶然性。
Description
技术领域
本发明设计一种基于蒙特卡洛仿真的河湖水质预测方法及富营养化风险评估方法,属于环境工程技术领域。
背景技术
随着现代社会的飞速发展,水体富营养化已成为一个全球性的重大水环境问题。目前,在世界范围内,受富营养化影响的沿海地区有500多个。此外,许多世界淡水湖泊、溪流和水库遭受富营养化。相关研究表明,我国当前已有80%的湖泊出现富营养化现象。在我国除一些人烟稀少的淡水湖区域以及还保持着原始状态的湖泊之外,其他湖泊营养盐水平的浓度已经超过了湖泊富营养化的标准浓度。尤其是我国的几大淡水湖的营养盐浓度都已经超过了磷、氮富营养化的发生浓度,总氮浓度甚至超过10倍之上。当下,我国的太湖以及巢湖都已经完全变为富营养化状态。因此,水环境问题已迫不容缓。水质预测能有效的掌握水质情况,预防大规模水污染事件的发生,对水污染防治具有现实意义。所以实现合理有效的水质预测,并在此基础上评估富营养化风险是解决水环境问题的重要思路。
已有的水质预测方法主要分为两类:基于机理的水质预测方法和数据驱动的水质预测方法。(1)基于机理的水质预测方法。这类方法主要是通过分析水体内部的机理演化过程来建立水质模型,再基于该水质模型对水质指标进行预测。常用的机理模型包括:QUAL2E、WASP和EFDC等。QUAL2E是一个一维水质模型,适用于模拟完全混合的枝状河流水质;它假定存在主流输送机理,即假定平流与扩散混合都是沿着河流的主流向,而在河流的横向与垂向上水质组分是完全均匀混合的,允许河流沿程有多个污染源、取水口以及支流汇入;其基本方程是一维平流—扩散质量迁移方程。WASP(The water quality analysissimulation program)是美国环保局(EPA)开发并推荐使用的用于各类水体水质模拟与分析的软件;它可用于模拟水动力学、河流一维不稳定流、湖泊和河口三维不稳定流、常规污染物和有毒污染物在水中的迁移和转化规律;其基本原理是平移—扩散方程。EFDC(TheEnvironmental Fluid Dynamics Code)模型是涵盖一维到三维地表水水质的数学模型;它可实现河流、湖泊、水库、湿地系统、河口和海洋等水体的水动力学、水质模拟、污染物迁移和泥沙输移等,是一个多参数有限差分模型。(2)数据驱动的水质预测方法。这类方法利用数据分析模型和工具对水质指标的实测数据进行建模,再通过此模型对水质指标进行预测。常用的数据驱动模型包括:人工神经网络、时间序列模型和支持向量机等,该类模型均是利用不同时刻的水质指标观测数据建立表征水质演化过程的数据驱动模型,并在此基础上实现水质过程预测。两类方法各具有一定优势,然而已有方法大都得到的是水质指标的取值变化预测(单值预测),其结果通常具有一定的偶然性,精确度和可信度难以令人信服,在此基础上的富营养化评价也会存在偶然性和不确定性。
发明内容
本发明为了弥补单值预测方法带来的不足,结合水质动力学模型,提出了一种基于蒙特卡罗仿真的河湖水质预测及富营养化风险评估方法,为水质预测和富营养化评估的不精确性问题提供一种新思路和可行解决方案。
本发明基于蒙特卡洛仿真的河湖水质预测方法及富营养化风险评估方法,通过以下步骤实现:
步骤一:基于水质指标,DO、Phyt、BOD、OP、PO4、ON、NH3-N和NO3-N,建立水质动力学模型;DO、Phyt、BOD、OP、PO4、ON、NH3-N和NO3-N分别为溶解氧、浮游植物、碳生化需氧量、有机磷、正磷酸盐、有机氮、氨氮和硝酸盐氮。
步骤三:选取Chl_a、总磷TP和总氮TN三个水质指标,计算三个水质指标的营养状态指数,通过权值归一化得出综合营养状态指数TLI;对TLI做分级处理,并计算得出TLI的概率分布和处于不同营养程度的概率;上述Chl_a为Phyt中叶绿素的浓度;TP为OP与PO4;TN为ON、NH3-N与NH3-N。
步骤二:利用蒙特卡罗仿真预测水质指标的演化过程,获得各水质指标在不同时间点上取值的概率分布。蒙特卡罗仿真的基本原理是当问题或对象本身具有概率特征时,可以用计算机模拟的方法产生抽样结果,根据抽样计算统计量或者参数的值;随着模拟次数的增多,可以通过对各次统计量或参数的估计值求平均的方法得到稳定结论。在基于蒙特卡罗仿真的预测过程中,水质指标演化过程会被多次模拟。对每一次模拟,每一个水质指标在每一个预测时间点的数值均不相同,这些数值同时具有规律性和随机性。规律性体现在预测值是基于确定的水质动力学模型、水质指标和模型参数初值而产生的;随机性体现在每一个水质指标在每一个预测时间点受噪声影响而得到的预测值不尽相同。当模拟次数充分大时,结合概率统计规律,可以获得各水质指标在不同时间点上取值的概率分布。
步骤三:富营养化风险评估方法的构建
选取Chl_a、总磷TP和总氮TN三个水质指标,计算三个水质指标的营养状态指数,通过权值归一化得出综合营养状态指数TLI;对TLI做分级处理,并计算得出TLI的概率分布和处于不同营养程度的概率。上述Chl_a为Phyt中叶绿素的浓度;TP为OP与PO4;TN为ON、NH3-N与NH3-N。
本发明的优点在于:
1)本发明基于蒙特卡洛仿真的河湖水质预测及富营养化风险评估方法,利用蒙特卡罗仿真预测8个水质指标的演化过程,获得8个水质指标取值的概率分布,实现水质预测。弥补了单一定值的水质预测的不精确性,降低了预测的偶然性。
2)本发明基于蒙特卡洛仿真的河湖水质预测及富营养化风险评估方法,通过构造综合营养状态指数,结合水质指标预测结果,计算综合营养状态指数的概率分布和处于不同营养程度的概率,实现富营养化风险评估。提高了水体的富营养化风险评估的准确度。
附图说明
图1是本发明整体流程图;
图2是水体生态系统演化过程图;
图3a是本发明所涉及的水质指标中溶解氧DO演化曲线图;
图3b是本发明所涉及的水质指标中浮游植物Phyt演化曲线图;
图3c是本发明所涉及的水质指标中碳生化需氧量BOD演化曲线图;
图3d是本发明所涉及的水质指标中有机磷OP演化曲线图;
图3e是本发明所涉及的水质指标中正磷酸盐PO4演化曲线图;
图3f是本发明所涉及的水质指标中有机氮ON演化曲线图;
图3g是本发明所涉及的水质指标中氨氮NH3-N演化曲线图;
图3h是本发明所涉及的水质指标中硝酸盐氮NO3-N演化曲线图;
图4a是本发明所涉及的水质指标中不同时刻下溶解氧DO的概率密度函数图;
图4b是本发明所涉及的水质指标中不同时刻下浮游植物Phyt的概率密度函数图;
图4c是本发明所涉及的水质指标中不同时刻下碳生化需氧量BOD的概率密度函数图;
图4d是本发明所涉及的水质指标中不同时刻下有机磷OP的概率密度函数图;
图4e是本发明所涉及的水质指标中不同时刻下正磷酸盐PO4的概率密度函数图;
图4f是本发明所涉及的水质指标中不同时刻下有机氮ON的概率密度函数图;
图4g是本发明所涉及的水质指标中不同时刻下氨氮NH3-N的概率密度函数图;
图4h是本发明所涉及的水质指标中不同时刻下硝酸盐氮NO3-N的概率密度函数图;
图5是应用本发明方法获得的不同时刻下综合营养状态指数TLI概率密度函数图;
图6a是应用本发明方法,预测天数为100天的贫营养概率曲线图;
图6b是应用本发明方法,预测天数为100天的中营养概率曲线图。
具体实施方式
下面将结合附图和实施例对本发明作进一步的详细说明。
本发明是基于蒙特卡洛仿真的河湖水质预测方法及富营养化风险评估方法,如图1所示,具体通过如下步骤实现:
步骤一:建立水质动力学模型。
水体生态系统演化过程,如图2所示,水体生态系统主要分为以下四个过程:溶解氧平衡过程、浮游植物动力学过程、磷循环和氮循环。水体生态系统主要由八个水质指标描述其演化过程:溶解氧(DO)、浮游植物(Phyt)、碳生化需氧量(BOD)、有机磷(OP)、正磷酸盐(PO4)、有机氮(ON)、氨氮(NH3-N)和硝酸盐氮(NO3-N)。基于上述水质指标演化机理,针对水体生态系统的富营养化过程建立水质动力学模型,其数学模型为上述8个水质指标方程组成的一组非线性连续微分方程组。
其中,分别表示水质指标DO、Phyt、BOD、OP、PO4、ON、NH3-N和NO3-N的浓度。Phyt的浓度以叶绿素a(Chl_a)的浓度为准,参数{k1,k2,…,k13}为模型参数,其含义和单位如表1所示。
表1水质动力学模型参数
上式中,f(x,θ)为水质动力学模型函数。
步骤二:基于蒙特卡洛仿真的河湖水质预测方法的构建。
将式(2)写成离散化形式如下:
上式中,x∈R8×1,θ∈R13×1,f:R8×1×R13×1→R8×1;水质指标的过程噪声模型参数过程噪声其中,Σx、Σθ为相应的噪声协方差阵;h是预测步长,h≤1day,t为时间点,t=0,h,2h,…,T,T为预测天数。
在基于蒙特卡罗仿真的预测过程中,水质指标演化过程会被多次模拟。对每一次模拟,每一个水质指标在每一个预测时间点的数值均不相同,这些数值同时具有规律性和随机性。规律性体现在预测值是基于确定的水质动力学模型、水质指标和模型参数初值而产生的;随机性体现在每一个水质指标在每一个预测时间点受噪声影响而得到的预测值不尽相同。当模拟次数充分大时,结合概率统计规律,可以获得各水质指标在不同时间点上取值的概率分布。基于蒙特卡罗仿真的水质预测的具体流程步骤为:
A、令水质指标的初值为x0,模型参数的初值为θ0,蒙特卡罗仿真粒子数为N,预测步长为h,初始时刻t=0。
x0 (i)=x0,θ0 (i)=θ0 (4)
C、对于时刻t=h,2h,…,T,基于t-1时刻的水质指标和模型参数值{xt-1 (i),θt-1 (i)},利用式(3)进行单步预测,得到t时刻的水质指标和模型参数预测值{xt (i),θt (i)}。
D:若t<T,t←t+h,返回进行步骤C,否则进行步骤E。
上式中,P(xt|(x0,θ0))为基于(x0,θ0)的水质指标概率密度函数,P(θt|(x0,θ0))为基于(x0,θ0)的模型参数概率密度函数,δ为狄拉克函数。
步骤三:富营养化风险评估方法的构建。
由于水体生态系统主要由八个水质指标描述其演化过程:溶解氧(DO)、浮游植物(Phyt)、碳生化需氧量(BOD)、有机磷(OP)、正磷酸盐(PO4)、有机氮(ON)、氨氮(NH3-N)和硝酸盐氮(NO3-N)。选取Chl_a、总磷(TP,TP=OP+PO4)和总氮(TN,TN=ON+NH3-N+NH3-N)三个水质指标,计算三个水质指标的营养状态指数,通过权值归一化得出综合营养状态指数TLI;对TLI做分级处理,并计算得出TLI的概率分布和处于不同营养程度的概率。具体流程步骤为:
a、对于时刻t=h,2h,…,T,粒子x(i),i∈[1,N]:
(1)对于Chl_a,CPhyt,t∈xt (i),有
TLIt(Chl_a)=10(2.500+1.086ln(CPhyt,t)) (6)
对于TP,COP,t、CPO4,t∈xt (i),有
TLIt(TP)=10(9.436+1.624ln(COP,t+CPO4,t)) (7)
(2)以Chl_a为基准参数,则Chl_a、TP、TN的权重分别为:
上式中,则w1、r11分别是Chl_a的权重和与自身的相关系数,w2、r12分别是TP的权重和与基准参数Chl_a的相关系数,w3、r13分别是TN的权重与基准参数Chl_a的相关系数。
(3)结合各水质指标的营养状态指数和权重,计算TLI:
TLIt=w1·TLIt(Chl_a)+w2·TLIt(TP)+w3·TLIt(TN) (10)
b、:对TLI做分级处理。TLI的值越大营养程度越重,以0—100分级,TLI<30为贫营养(Oligotrophic),30<TLI<50为中营养(Mesotrophic),TLI>50为富营养(Eutropher)。
c、计算TLI的概率密度函数:
上式中,为t时刻预测结果的TLIt属于区间[Em,Em+1]的概率,B[Em,Em+1]为TLIt属于区间[Em,Em+1]的个数,区间[Em,Em+1]为TLIt的分级区间,Em与Em+1分别为指TLIt的第m个分级区间的上限和下限,即当属于第m个分级区间时,TLIt Em<TLIt<Em+1。
实施例1:
步骤一:建立水质动力学模型
基于水质指标演化机理,建立水质动力学模型,模型方程如式(1)所示。
步骤二:基于蒙特卡洛仿真的河湖水质预测方法的构建
水质指标初值:x0=[6.5,0.05019,3,0.07,0.06,0.35,0.36,0.1]T,
水质指标过程噪声协方差阵:Σx=diag{0.1,0.001,0.01,0.001,0.001,0.01,0.01,0.01},
模型参数初值:θ0=[0.1,1.208,0.01,0.01,0.15,0.13,0.0009,0.001,0.025,0.075,0.9,0.35,1.02,1.3]T,
模型参数过程噪声协方差阵:
Σθ=diag{0.1,0.01,0.001,0.1,0.001,0.001,0.0001,0.0001,0.0001,0.01,0.001,0.01,0.1}。
预测步长h=1/24(单位:day)、粒子数N=1000和仿真时长T=100(单位:day)。
按照式(1)~(4)和上述参数数值计算8个水质指标演化结果,结果如图3a~3h所示。按照式(5)计算不同时刻下的8个水质指标的概率密度函数,结果如图4a~4h所示。
步骤三:富营养化风险评估方法的构建
基于上述水质指标演化结果,构造综合营养状态指数,计算综合营养状态指数的概率分布和处于不同营养程度的概率,实现富营养化风险评估。选取Chl_a(即Phyt)、总磷(TP)和总氮(TN)三个水质指标,计算三个水质指标的营养状态指数,通过权值归一化得出TLI;对TLI做分级处理,并计算得出TLI的概率分布和处于不同营养程度的概率。
中国湖泊(水库)的基准参数Chl_a与其他参数之间的相关系数r1j如下表2所示。
表2中国湖泊(水库)的基准参数Chl_a与其他参数之间的相关关系rij
按照式(6)~(11)和表2计算综合营养状态指数的概率密度函数,结果如图5所示,按照式(12)计算水体处于不同营养状态的概率,结果如图6a、6b所示。
Claims (1)
1.基于蒙特卡洛仿真的河湖水质预测方法及富营养化风险评估方法,其特征在于:通过以下步骤实现:
步骤一:基于水质指标,DO、Phyt、BOD、OP、PO4、ON、NH3-N和NO3-N,建立水质动力学模型;DO、Phyt、BOD、OP、PO4、ON、NH3-N和NO3-N分别为溶解氧、浮游植物、碳生化需氧量、有机磷、正磷酸盐、有机氮、氨氮和硝酸盐氮;
步骤二:利用蒙特卡罗仿真预测水质指标的演化过程,获得各水质指标在不同时间点上取值的概率分布;具体为:
将式(2)写成离散化形式如下:
上式中,x∈R8×1,θ∈R13×1,f:R8×1×R13×1→R8×1;水质指标的过程噪声模型参数过程噪声其中,Σx、Σθ为相应的噪声协方差阵;h是预测步长,h≤1day,t为时间点,t=0,h,2h,…,T,T为预测天数;
在基于蒙特卡罗仿真的预测过程中,水质指标演化过程会被多次模拟;对每一次模拟,每一个水质指标在每一个预测时间点的数值均不相同,这些数值同时具有规律性和随机性;规律性体现在预测值是基于确定的水质动力学模型、水质指标和模型参数初值而产生的;随机性体现在每一个水质指标在每一个预测时间点受噪声影响而得到的预测值不尽相同;当模拟次数充分大时,结合概率统计规律,可以获得各水质指标在不同时间点上取值的概率分布;基于蒙特卡罗仿真的水质预测的具体流程步骤为:
A、令水质指标的初值为x0,模型参数的初值为θ0,蒙特卡罗仿真粒子数为N,预测步长为h,初始时刻t=0;
x0 (i)=x0,θ0 (i)=θ0 (2)
C、对于时刻t=h,2h,…,T,基于t-1时刻的水质指标和模型参数值{xt-1 (i),θt-1 (i)},利用式(1)进行单步预测,得到t时刻的水质指标和模型参数预测值{xt (i),θt (i)};
D:若t<T,t←t+h,返回进行步骤C,否则进行步骤E;
上式中,P(xt|(x0,θ0))为基于(x0,θ0)的水质指标概率密度函数,P(θt|(x0,θ0))为基于(x0,θ0)的模型参数概率密度函数,δ为狄拉克函数;
步骤三:选取Chl_a、总磷TP和总氮TN三个水质指标,计算三个水质指标的营养状态指数,通过权值归一化得出综合营养状态指数TLI;对TLI做分级处理,并计算得出TLI的概率分布和处于不同营养程度的概率;上述Chl_a为Phyt中叶绿素的浓度;TP为OP与PO4;TN为ON、NH3-N与NH3-N;具体为:
a、对于时刻t=h,2h,…,T,粒子x(i),i∈[1,N]:
(1)对于Chl_a,CPhyt,t∈xt (i),有
TLIt(Chl_a)=10(2.500+1.086ln(CPhyt,t)) (3)
(2)以Chl_a为基准参数,则Chl_a、TP、TN的权重分别为:
上式中,则w1、r11分别是Chl_a的权重和与自身的相关系数,w2、r12分别是TP的权重和与基准参数Chl_a的相关系数,w3、r13分别是TN的权重与基准参数Chl_a的相关系数;
(3)结合各水质指标的营养状态指数和权重,计算TLI:
TLIt=w1·TLIt(Chl_a)+w2·TLIt(TP)+w3·TLIt(TN) (7)
b、:对TLI做分级处理;TLI的值越大营养程度越重,以0—100分级,TLI<30为贫营养,30<TLI<50为中营养,TLI>50为富营养;
c、计算TLI的概率密度函数:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810269206.2A CN108615090B (zh) | 2018-03-29 | 2018-03-29 | 基于蒙特卡洛仿真的河湖水质预测及富营养化风险评估方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810269206.2A CN108615090B (zh) | 2018-03-29 | 2018-03-29 | 基于蒙特卡洛仿真的河湖水质预测及富营养化风险评估方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108615090A CN108615090A (zh) | 2018-10-02 |
CN108615090B true CN108615090B (zh) | 2020-07-17 |
Family
ID=63659205
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810269206.2A Active CN108615090B (zh) | 2018-03-29 | 2018-03-29 | 基于蒙特卡洛仿真的河湖水质预测及富营养化风险评估方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108615090B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109272171A (zh) * | 2018-10-16 | 2019-01-25 | 天津大学 | 一种预测暴雨期二级河道排水对干流水质动态影响的方法 |
CN109614570B (zh) * | 2018-11-15 | 2023-04-18 | 北京英视睿达科技股份有限公司 | 预测断面水质参数数据的方法及装置 |
CN110889201B (zh) * | 2019-11-05 | 2020-07-03 | 北京工商大学 | 基于混杂动态系统的水质机理建模及健康风险评估方法 |
CN110956365B (zh) * | 2019-11-11 | 2020-07-17 | 北京工商大学 | 基于混杂贝叶斯网络的小麦粉供应链的菌落总数动态风险评估方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103425853A (zh) * | 2012-05-16 | 2013-12-04 | 上海市环境科学研究院 | 一种湖泊富营养化控制标准定值及评估方法 |
CN103810537A (zh) * | 2014-02-12 | 2014-05-21 | 南京大学 | 一种基于水质模型的区域环境风险评估方法 |
KR20140140361A (ko) * | 2013-05-29 | 2014-12-09 | 창원대학교 산학협력단 | 앙상블 유출 예측기법을 적용한 하천 수질예측 시스템 |
CN105095997A (zh) * | 2015-07-30 | 2015-11-25 | 浙江大学 | 基于蒙特卡洛和层次分析法的突发性水污染事故预警方法 |
CN107688701A (zh) * | 2017-08-23 | 2018-02-13 | 北京工商大学 | 基于wasp模型的水质软测量及水体富营养化评价方法 |
-
2018
- 2018-03-29 CN CN201810269206.2A patent/CN108615090B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103425853A (zh) * | 2012-05-16 | 2013-12-04 | 上海市环境科学研究院 | 一种湖泊富营养化控制标准定值及评估方法 |
KR20140140361A (ko) * | 2013-05-29 | 2014-12-09 | 창원대학교 산학협력단 | 앙상블 유출 예측기법을 적용한 하천 수질예측 시스템 |
CN103810537A (zh) * | 2014-02-12 | 2014-05-21 | 南京大学 | 一种基于水质模型的区域环境风险评估方法 |
CN105095997A (zh) * | 2015-07-30 | 2015-11-25 | 浙江大学 | 基于蒙特卡洛和层次分析法的突发性水污染事故预警方法 |
CN107688701A (zh) * | 2017-08-23 | 2018-02-13 | 北京工商大学 | 基于wasp模型的水质软测量及水体富营养化评价方法 |
Non-Patent Citations (1)
Title |
---|
蒙特卡洛算法在水动力水质模型(DYRESM-CAEDYM)参数优选中的应用;翟海涛 等;《生态科学》;20140131;第33卷(第1期);第38-45页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108615090A (zh) | 2018-10-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108615090B (zh) | 基于蒙特卡洛仿真的河湖水质预测及富营养化风险评估方法 | |
Cox | A review of currently available in-stream water-quality models and their applicability for simulating dissolved oxygen in lowland rivers | |
CN105913411B (zh) | 一种基于因子定权模型的湖泊水质评价预测系统及方法 | |
Wakatsuki et al. | Rates of weathering and soil formation | |
Sandin et al. | Local, landscape and regional factors structuring benthic macroinvertebrate assemblages in Swedish streams | |
Laine | Adaptive MCMC methods with applications in environmental and geophysical models | |
CN107688701B (zh) | 基于wasp模型的水质软测量及水体富营养化评价方法 | |
CN105631537A (zh) | 基于气象服务平台的空气质量预报业务系统 | |
Zhang et al. | Integration of a benthic sediment diagenesis module into the two dimensional hydrodynamic and water quality model–CE-QUAL-W2 | |
CN114240196A (zh) | 一种基于水动力-水质-生态模型的湖泊多水源调控方法 | |
Xu et al. | Parameter uncertainty and sensitivity analysis of the three Gorges Reservoir and Xiangxi River EFDC model | |
CN115062874B (zh) | 一种水体污染物监测预测分析方法及系统 | |
Memiş et al. | An application of fuzzy parameterized fuzzy soft matrices in data classification | |
CN109033648B (zh) | 基于果蝇优化算法的水质机理建模及水质预测方法 | |
Osidele et al. | Identification of model structure for aquatic ecosystems using regionalized sensitivity analysis | |
Tsiaras et al. | Model simulations of a mesocosm experiment investigating the response of a low nutrient low chlorophyll (LNLC) marine ecosystem to atmospheric deposition events | |
Wang et al. | A two-dimensional numerical model for eutrophication in Baiyangdian Lake | |
Hinrichsen et al. | Biogeochemical and physical controls of nitrogen fluxes in a highly dynamic marine ecosystem—model and network flow analysis of the Baltic Sea | |
Pauer et al. | The impact of alternative nutrient kinetics and computational grid size on model predicted primary production and hypoxic area in the northern Gulf of Mexico | |
Makarova et al. | Estimating mercury footprint in the regions of the Russian Federation | |
Paredes et al. | Two mathematical approaches to study the phosphorus eutrophication of a wetland in Puerto Rico | |
AL KINDI et al. | APPLICATION OF A ONE-DIMENSIONAL STEADY STATE MODEL FOR SIMULATION THE WATER QUALITY IN TIGIRS RIVER | |
Mouaouia et al. | Mathematical modelling of nitrogen removal in horizontal subsurface flow constructed wetland | |
Lowrance | Water Quality: Modeling | |
CN116562176B (zh) | 一种基于耦合神经网络与水文物理模型的径流模拟方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |