中国传媒大学2010 学年第一学期数学建模与数学实验课程数学建模与数学实验题目Pristine湖污染问题的建模与求解学生姓名学号班级学生所属学院任课教师教师所属学院成绩Pristine湖污染问题的建模与求解摘要本文讨论了湖水污染浓度变化趋势的预测问题。
通过分析水流输入输出湖泊的过程，建立了湖水污染浓度随时间变化的含参变量的微分方程模型，在河水污染浓度恒定和自然净化速率呈线性关系的情况下，求得其精确解，带入具体数据得到结论：在PCA声称的河水污染浓度下，湖的环境不会恶化；在工作人员实地测得的河水浓度下，湖的环境将会恶化。
同时建立了计算机模拟模型，带入具体数值，运用时间步长法来仿真模拟了在湖水污染浓度稳定以前湖水每天的变化情况，输出自PCA建厂以来每年的湖水污染浓度，得到与微分方程模型相同的结论。
在全停产和半停产时，通过前面的两个模型可以计算湖水污染浓度在自然净化影响下的恢复到净化指标所需的年限。
并可得到结论：在半停产状态下，在选定的自然净化速率常数的约束下，只有当河水污染浓度降至原来的3.15%（自然净化速率呈线性关系），4.7%（自然净化速率呈指数关系），才有可能使河水在100年内恢复至0.001mol/l，然后给出整改建议。
一、问题重述Pure河是流入Pristine湖的唯一河流。
50年前PCA公司在此河旁建起一个生产设施并投入运行。
PCA将为处理的湖水排入河中，导致Pristine湖被污染。
PCA公司声称：已排放的废水的标准多年从未改变切不会对湖的环境有影响。
现已知：Pristine湖的湖容量为L，流入（流出）的水流速度为L/年。
PCA公司声称河水污染浓度仅为0.001mol/L，自工厂以来没有改变过。
讨论下列问题：（1）建立数学模型用PCA提供的公开数据判断湖的环境是否会恶化；（2）以目前湖水污染浓度0.03mol/L，和河水污染浓度0.05mol/L为新数据判断湖的环境是否会恶化；二、模型的合理假设和符号系统2.1 模型的合理假设（1）降水量和增发量相等；（2）湖中流入量和流出量相等且一直未变；（3）污水量远小于河水注入量，且污水与河水混合均匀；（4）湖水混合均匀，且流入污水的扩散速度无限大；（5）湖内除Pure河外，无其他污染源；2.2 符号系统：河水污染浓度mol/L；：湖水污染物浓度mol/L；：湖泊容量L；：自然净化速率mol/（L。
年）：流入（流出）的水流速度L/年；：从PCA建厂至考察时刻的时间段。
三、问题的分析3.1 问题分析：对于问题中几个词语的理解：1.是否会恶化——湖的环境恶化即湖水污染浓度大于0.001mol/L，要判断其是否会恶化，则需计算在某一污染物积累速度（分析影响此速度的因素）下，湖水能达到的最大污染浓度和其变化趋势，以及湖水经几年超过0.001mol/L，经过几年达到最大污染浓度。
2.自然净化——自然净化是独立的生态系统进行自我调节的方式之一，是在空气，阳光，水和细菌的参与下，进行包括物理沉降，化学反应和生物转化三大方面的活动，其最终作用是将污染物转化为无害物质，从而净化生态系统。
当河水输入湖泊并均匀混合之后，影响湖水污染浓度的唯一因素便是自然净化速度。
湖水污染问题水流的动态流程图：此问题中，我们考察对象是湖水污染浓度的变化趋势：1.在整改之前，其增加的趋势，超过净化指标0.001ml/L（即湖水恶化）的可能性和时限；2.在全停产（无污染物输入）和半停产的情况下，其降低的趋势，达到净化指标的时限。
3.在整改之后，其增加的趋势，未定与净化指标之下某一水平的时限、在前假设条件的基础之上，湖水容量不变，出河水外无其他的污染源，故我们可将湖泊作为一个封闭的生态系统，其简化的湖水被污染的动态过程为：受污河输入湖泊，河水与湖水均匀混合，受污河水进行自我的净化，湖水数出湖泊。
湖中污染物的量直接决定了湖水污染浓度，而污染物的量受到以下两方面因素的影响：1.河水的污染浓度及其流入速度（根据已知此速度不变），2.湖水的自然净化速度，前者使其增加，后者使其减少（负增加）。
问题一、二、三的实质都是要分析污染浓度的变化趋势，其去表便在于前一因素的不同。
问题一中，河水污染浓度不变，恒为0.001mol/L；问题二中，河水污染浓度可能会变化，受PCA效益的影响而按一定规律波动；问题三中，在全停产或半停产的情况下，和硕污染浓度为0或减为问题二中的一部分。
后一因素（自然净化速度）在三个问题中的作用都是相同的。
根据微积分的知识可知，在适当短的时间段之内，通过建立微分方程，可以将连续的过程离散化，从而可以得到湖水污染浓度与时间之间的关系式。
利用时间步长发，缩小步长值（从年到月到天），并与微分方程所得的精确解做出比较。
四、模型建立与求解问题一：根据PCA的公开申明和所提供数据，可认为：河水污染浓度恒为0.001mol/L。
从存在自然净化和不存在自然净化两个方面考虑：（1）．在不考虑自然净化的情况下：由于假设湖中流入量和流出量相等，而在经过与湖水均匀混合后，流出湖水污染浓度明显减小，故流出污染物的量小于流出污染物的量，污染物将在湖中沉积，从而使湖水污染浓度增加，当其增加至于输入的河水污染浓度相等时，河水污染浓度达到最大，并稳定在这一数值，不在增加。
建立湖水污染浓度随时间变化的微分方程模型：设在极短时间，湖水污染浓度增加，在将湖水被污染这一连续动态过程简化为离散的瞬间静止状态（如问题分析中所述）之后，根据湖中剩余量=输入量—输出量，我们可以列微分方程如下：化简可得： (1)带入L，L/年和mol/L的数据，我们可得到和的关系式如下：通过此关系式我们可知，当时，湖水污染浓度将趋近与0.001mol/L，即湖的环境不会恶化。
建立计算机模拟模型：湖水污染浓度的变化时有湖中污染物随时间的积累而引起的，这个逐步积累的过程我们可以用计算机进行仿真模拟，其实质为完成一个循环累加的过程，并可改变时间步长，如一年一年的累积，一月一月的累积，一天一天的累积，从而使我们的模拟值逐步精确，可与微分方程求得精确解比较，分析误差。
为提高模拟结果的精确性和运算的效率，我们采用了逐天累加，数出年污染浓度的方式。
模拟程序见附件一，从后面的分析中我们可知：在此河水污染浓度恒定，无自然净化的最简单的情况下建立的模型是以后问题的基础，后面的问题只是改变条件或数据，其实质是不变的，故我们在此程序中加入了多个选择语句，在不同德条件或数据下执行不同的命令，从而用一个程序解决全部的模拟问题。
模拟所得的数据如下表1：年份 1 3 5 8 101.73 4.35 6.14 7.82 8.51年份15 20 30 40 509.43 9.78 9.97 9.995 9.999根据模拟数据所作的湖水污染浓度变化趋势图如下图1：模拟所得的数据显示：湖水污染浓度将稳定于0.001mol/L，这与分析微分方程所得的结果是相符的。
（2）在考虑自然净化的情况下：A.首先，我们应该了解什么事自然净化。
根据资料显示：湖水中的污染物可分为有机污染物和无机污染物两大类，在多种环境因素（阳光，空气，水，水中生物，水中化学物质，重力等）的作用下，通过物理沉降，化学反应和生物转化一系列复杂的活动，它们的量会发生变化。
有机物在水环境中的迁移转化过程如图2：可以看出整个过程是相当复杂的，不仅过程多，而且在相同的过程中，不同的物质有着不同的结果。
此问题中没有明确给出输入湖中的污染物种类，也没有对湖泊环境做任何描述，无疑给问题的解决增加了极大的难度。
为是问题简化，我们从一般的情况出发，假设：湖中所进行的反应均为一级反应，有机物的存在不会对环境参数造成改变，环境固定（自然净化速率恒定），考虑湖泊生态系统中起主要作用的几种过程作简要分析：1.物理沉降。
不同物质有着不同的沉降速度常数，沉降速率为；2.挥发。
不同物质有着不同的挥发速率常数，在有机物在水体上的大气中的分压为0的条件下，挥发速率为(为水体深度)3.水解反应。
不同物质有着不同的水解速率常数，在一级反应的条件下，水解速率为；4.生物降解反应。
不同物质有着不同的降解速率常数，在一级反应的条件下，生物降解速率；B.在最简模拟程序的基础上，从每天的积累量中减去每天的自然净化量时间步长，重复循环可得到逐年湖水污染浓度的值，设。
越大，湖水污染浓度将越快稳定于一个更小的浓度值。
模拟所得的数据如下表2：年份 1 3 5 8 101.587 3.437 4.315 4.878 4.984年份15 20 30 40 505.087 5.103 5.106 5.106 5.106根据模拟数据所作的湖水污染浓度变化趋势图如下图3：C.进一步考虑自然净化速率的影响。
更贴近于实际的情况是，自然净化速率与湖水污染浓度成指数关系，随着的增加而增加，但增加的速率会逐渐减小，用关系式表达即可为：，A,B为与环境和污染物种类有关的常数。
将此关系式带入微分方程（2），得到一个新的微分方程，此方程无解。
但是我们可以通过计算机模拟求得数值解。
A的大小影响着最终稳定浓度，B的大小影响着达到最终稳定浓度的快慢。
在前面达到线性关系的基础上，在=0/。
001mol/L处的值大小确定A至数量级，并稳定年限定在10—20年间，从而得到一组估计值：模拟所得数据如下表3：年份 1 3 5 8 101.4932.9713.539 3.816 3.871年份15 20 30 40 503.903 3.906 3.906 3.906 3.906根据模拟所得数据所作的湖水污染浓度变化趋势图如下图4：观察数据可知：无论自然净化速率与污染浓度成线性关系还是指数关系，当河水污染浓度=0.001mol/L时，都将稳定于一个小于0.001mol/L的值，也即：湖的环境不会恶化。
问题二：根据实际情况，在此我们只考虑存在自然净化的情况。
（1）.河水污染浓度恒定为0.05mol/L:与问题一（2）相同，只是=0.05mol/L，沿用问题一（2）的微分方程模型和差分模拟模型，我们可以得到以下结果：A.微分方程模型（自然净化速率与湖水污染浓度成线性关系）：当时间时，，即湖水污染浓度稳定于一个与K有关的值。
当=0.05mol/L时，如要稳定于0.001mol/L则K值为。
根据我们设定的可计算，湖水最终污染浓度=0.0317mol/L，超过净化指标0.001mol/L，故在此条件下，湖的环境将会恶化。
B.差分模拟模型：在与成线性关系时，得到模拟数据如下表4：年份 1 3 5 8 100.82 1.88 2.46 2.88 3.01年份15 20 30 40 503.13 3.16 3.17 3.17 3.17在与成指数关系时，得到模拟数据如下表5：年份 1 3 5 8 100.81 2.01 2.81 3.59 4.02年份15 20 30 40 504.32 4.48 4.56 4.58 4.58观察以上数据可知：由于河水污染没得过大，河水污染没得从第一年起就超过了净化指标，并逐年增加，是湖的环境恶化。