湖水污染问题的数学建模与求解

中国传媒大学2010 学年第一学期数学建模与数学实验课程

数学建模与数学实验

题目Pristine湖污染问题的建模与求解

学生姓名

学号

班级

学生所属学院

任课教师

教师所属学院

成绩

Pristine湖污染问题的建模与求解

摘要

本文讨论了湖水污染浓度变化趋势的预测问题。

通过分析水流输入输出湖泊的过程，建立了湖水污染浓度随时间变化的含参变量的微分方程模型，在河水污染浓度恒定和自然净化速率呈线性关系的情况下，求得其精确解，带入具体数据得到结论：在PCA声称的河水污染浓度下，湖的环境不会恶化；在工作人员实地测得的河水浓度下，湖的环境将会恶化。

同时建立了计算机模拟模型，带入具体数值，运用时间步长法来仿真模拟了在湖水污染浓度稳定以前湖水每天的变化情况，输出自PCA建厂以来每年的湖水污染浓度，得到与微分方程模型相同的结论。

在全停产和半停产时，通过前面的两个模型可以计算湖水污染浓度在自然净化影响下的恢复到净化指标所需的年限。并可得到结论：在半停产状态下，在选定的自然净化速率常数的约束下，只有当河水污染浓度降至原来的%（自然净化速率呈线性关系），%（自然净化速率呈指数关系），才有可能使河水在100年内恢复至l，然后给出整改建议。

一、问题重述

Pure河是流入Pristine湖的唯一河流。50年前PCA公司在此河旁建起一个生产设施并投入运行。PCA将为处理的湖水排入河中，导致Pristine湖被污染。PCA公司声称：已排放的废水的标准多年从未改变切不会对湖的环境有影响。

10L，流入（流出）的水流速度为149.1L/现已知：Pristine湖的湖容量为15

年。PCA公司声称河水污染浓度仅为L，自工厂以来没有改变过。

讨论下列问题：

（1）建立数学模型用PCA提供的公开数据判断湖的环境是否会恶化；

（2）以目前湖水污染浓度L，和河水污染浓度L为新数据判断湖的环境是否会恶化；

二、模型的合理假设和符号系统

模型的合理假设

（1）降水量和增发量相等；

（2）湖中流入量和流出量相等且一直未变；

（3）污水量远小于河水注入量，且污水与河水混合均匀；

（4）湖水混合均匀，且流入污水的扩散速度无限大；

（5）湖内除Pure河外，无其他污染源；

符号系统

0ρ：河水污染浓度mol/L ；

ρ：湖水污染物浓度mol/L ； V ：湖泊容量15

10L ；

c ：自然净化速率mol/（L 。年）

μ：流入（流出）的水流速度14

9.1L/年；

t ：从PCA 建厂至考察时刻的时间段。

三、问题的分析

问题分析：

对于问题中几个词语的理解：

1.是否会恶化——湖的环境恶化即湖水污染浓度大于L ，要判断其是否会恶化，则需计算在某一污染物积累速度（分析影响此速度的因素）下，湖水能达到的最大污染浓度和其变化趋势，以及湖水经几年超过L ，经过几年达到最大污染浓度。

2.自然净化——自然净化是独立的生态系统进行自我调节的方式之一，是在空气，阳光，水和细菌的参与下，进行包括物理沉降，化学反应和生物转化三大方面的活动，其最终作用是将污染物转化为无害物质，从而净化生态系统。当河水输入湖泊并均匀混合之后，影响湖水污染浓度的唯一因素便是自然净化速度。 湖水污染问题水流的动态流程图：

此问题中，我们考察对象是湖水污染浓度的变化趋势：

1.在整改之前，其增加的趋势，超过净化指标L （即湖水恶化）的可能性和时限；

2.在全停产（无污染物输入）和半停产的情况下，其降低的趋势，达到净化指标的时限。

3.在整改之后，其增加的趋势，未定与净化指标之下某一水平的时限、

在前假设条件的基础之上，湖水容量不变，出河水外无其他的污染源，故我们可将湖泊作为一个封闭的生态系统，其简化的湖水被污染的动态过程为：受污

河输入湖泊，河水与湖水均匀混合，受污河水进行自我的净化，湖水数出湖泊。湖中污染物的量直接决定了湖水污染浓度，而污染物的量受到以下两方面因素的影响：1.河水的污染浓度及其流入速度（根据已知此速度不变），2.湖水的自然净化速度，前者使其增加，后者使其减少（负增加）。问题一、二、三的实质都是要分析污染浓度的变化趋势，其去表便在于前一因素的不同。问题一中，河水污染浓度不变，恒为L ；问题二中，河水污染浓度可能会变化，受PCA 效益的影响而按一定规律波动；问题三中，在全停产或半停产的情况下，和硕污染浓度为0或减为问题二中的一部分。后一因素（自然净化速度）在三个问题中的作用都是相同的。

根据微积分的知识可知，在适当短的时间段之内，通过建立微分方程，可以将连续的过程离散化，从而可以得到湖水污染浓度与时间之间的关系式。

利用时间步长发，缩小步长值（从年到月到天），并与微分方程所得的精确解做出比较。

四、模型建立与求解

问题一：

根据PCA 的公开申明和所提供数据，可认为：河水污染浓度恒为L 。从存在自然净化和不存在自然净化两个方面考虑： （1）．在不考虑自然净化的情况下：

由于假设湖中流入量和流出量相等，而在经过与湖水均匀混合后，流出湖水污染浓度明显减小，故流出污染物的量小于流出污染物的量，污染物将在湖中沉积，从而使湖水污染浓度增加，当其增加至于输入的河水污染浓度相等时，河水污染浓度达到最大，并稳定在这一数值，不在增加。

建立湖水污染浓度随时间变化的微分方程模型：设在极短时间dt ，湖水污染浓度

增加d ρ，在将湖水被污染这一连续动态过程简化为离散的瞬间静止状态（如问题分析中所述）之后，根据湖中剩余量=输入量—输出量，我们可以列微分方程如下：

0dt V

d V dt ρμρρρ

μ⨯⨯+⨯=++⨯

化简可得：0()d dt V ρμ

ρρ=⨯- (1)

带入15

10V =L ，141.910μ=⨯L/年和00.001ρ=mol/L 的数据，我们可得到ρ和t 的关系式如下：

0.190.190.001(1)t t

e e ρ-⨯⨯=⨯⨯- 通过此关系式我们可知，当t =∞时，湖水污染浓度将趋近与L ，即湖的环境不会