基于差分方程的湖水污染问题
- 格式：docx
- 大小：1.03 MB
- 文档页数：15
在对于整个大环境的探究过程中，我们暂定认为，该水域由一条流进湖的pure河中间的Pristine湖和一条流出湖的组成。对于该水域的动态变化过程可以认为，该水域的污染来源是PCA公司的污染投放，该水域对于污染的净化能力分为其湖水的自我净化能力和流出该湖的清洁能力。
根据微积分的知识可知，在适当短的时间段之内，通过建立微分方程，可以将连续的过程离散化，从而可以得到湖水污染浓度与时间之间的关系式。
5.3对于净化标准划定的分析
根据题意已知，当该水域的污染程度达到0.001mol/l，才达到需要治理水质的要求，这是一个基本红线。
5.4对于污染源划定的基本分析
污染源的划定，我们假设其作为含有某系列物质的溶液，其主要分为含有不同污染能力的固态物质若等，对于题三的整改部分，我们可以对公司提供一定的建议。鉴于公司不能直接切实有效对污染源的量做到减少，只有在排放的部分中做一定改善。对于不同污染能力的污染源物质，可以采取固态物质粉碎颗粒化，对于溶液进一步稀释。
问题：
1.在花费时间和经费去测试之前，建立数学模型用PCA提供的公开数据判断湖的环境是否会恶化。
2.派出野外工作人员测得目前湖水污染浓度为0.03mol/l，再测得河水污染浓度为0.05mol/l。以新数据为依据考虑湖水污染问题的数学模型。
3.现在假设你是环保局的所聘请的高级顾问，请向你的雇主提交一份报告.内容包括:
观察以上各式可发现：在一定条件下，各个速率与湖中污染程度呈线性关系，又因为某种化学物的总消失速率等于各过程消失速率之和，故我们可以认为：自然净化速率 （ 为正常数，与环境及污染种类有关）
如上所述，我们可得上述表达式的改进。
（3）
把该式求解得
（4）
当时间t= 时， ，即湖水污染浓度稳定于一个与K有关的值。
的大小影响着 达到最终稳定浓度， 的大小影响着 达到最终稳定浓度的快慢。在前面达到线性关系的基础上，在 =0.001 的c值大小确定A至数量级，并稳定年限在10-20年间，从而得到一股估计值。
= ， 1
3
5
8
10
1.493
2.971
3.539
3.816
3.871
15
关键字：时间步长法微分方程模型logistic模型
二问题重述
Pure河是流入Pristine湖的唯一河流。50年前PCA公司在此河旁建起一个生产设施并投入运行。PCA将为处理的废水排入河中,导致了Pristine湖被污染。PCA公司声称:已排放的废水的标准多年从未改变且不会对湖的环境有影响。
Pristine湖的湖容量为 ，流入(流出)的水流速度为 年。PCA声称河水污染浓度仅为0.001mol/l，自工厂开工以来没有改变过。
同时建立了计算机模拟模型，带入具体数值，运用时间步长法来仿真模拟了在湖水污染浓度稳定以前湖水每天的变化情况，输出自PCA建厂以来每年的湖水污染浓度变化，得到与微分方程模型相同的结论。
在全停产和半停产时，通过前面的两个模型可以计算湖水污染浓度在自然进化影响下的恢复到自然净华指标所需年限。并可得到结论：在半停产条件下，在选定的自然净化速率常数的约束下，只有当河水污染浓度降至原来的3.15%（自然净化速率成线性关系），4.7%（自然净化速率呈指数关系），才有可能使河水在100年内恢复至0.001 ，然后给出整改建议。
20
30
40
50
3.903
3.906
3.906
3.906
3.906
表3
根据模拟所得数据所做的湖水污染浓度变化趋势如下图6所示：
图6
观察数据可知：无论自然净化速率c与湖水污染浓度 呈线性关系还是指数关系，当湖水污染浓度 =0.001 时， 都将稳定于一个小于0.001 的值，也即：湖的环境不会恶化。
问题二：
模拟程序见附件一，从后面的分析中我们可知：在此河水污染浓度恒定，无自然净化的最简单的情况下建立的模型是以后问题的基础，后面的问题只是改变条件或数据，其实质是不变的，故我们在此程序中加入了多个选择语句，在不同的条件或数据下执行不同的命令，从而用一个程序解决了全部的模拟问题。
模拟所得数据如下表1所示：
年份
根据实际情况，我们在此只考虑存在自然净化的情况。
5.在整个大环境下，河水和湖水在污染源的浓度不一致。
6.在时间允许的范围内，假设整个湖水容量不发生变化。
7.污染源属于什么态的物质
四符号说明和名词解释
河水污染浓度
湖水污染物浓度
V
湖泊容量
L
c
自然净化速率
流入（流出）的水流速度
t
从PCA建厂至考察时刻的时间段
天
五问题分析
5.1对于水域基本情况分析
由于题目假设可得，我们对于整个水质的大环境假设为一个容积不变的大容器。基于在该流域里，降水量和蒸发量保持相等，湖水流入量和流出量相等且一直不变，所以我们假设该水域环境下，其容积不发生变化。
综上所述，对于该题，我们得到如下基本分析。在前假设条件的基础之上，湖水容量不变，出河水外无其他的污染源，故我们可将湖泊作为一个封闭的生态系统，其简化的湖水被污染的动态过程为：受污河水输入湖水中，河水与湖水均匀混合，受污水质在湖水区域内进行自我净化，湖水从输出端流出湖泊。湖中污染物的量直接决定了湖水污染浓度，而污染物的量受到以下两方面因素的影响：1.河水的污染浓度及其流入速度（根据已知此速度不变），
2.湖水的自我净化速度，前者使其增加，后者使其减少。
问题一、二、三的实质都是要分析污染浓度的变化趋势，而关键就在于前一因素的不同。问题一中，河水污染浓度不变，恒为0.001mol/L；问题二中，河水污染浓度可能会变化，受PCA效益的影响而按一定规律波动；问题三中，在全停产或半停产的情况下，和硕污染浓度为0或减为问题二中的一部分。后一因素（自然净化速度）在三个问题中的作用都是相同的。
根据所得数据如下表2所示
年份
1
3
5
8
10
1.587
3.437
4.315
4.878
4.984
年份
15
20
30
40
50
5.087
5.103
5.106
5.106
5.106
表2
根据模拟数据所作的湖水污染浓度变化趋势如图5所示，
图5
C进一步考虑自然净化速率的影响。更贴切于实际的情况是，自然净化速率c与湖水污染程度 成指数关系，c随的增加而增加，但增加的速率会逐步减小，用关系式表达即可为： ， 为环境和污染程度种类有关的常数。将此关系式代入微分方程（3），得到一个新的微分方程，此方程无解。但是我们可以通过此关系式，建立计算机模拟求得数值解。
建立湖水污染浓度随时间变化的微分方程模型：设在极短时间dt，湖水污染浓度增加d ，在将湖水被污染这一连续动态过程简化为离散的瞬间静止状态（如问题分析中所述）之后，根湖中剩余量=输入量—输出量，我们可以列微分方程如下
我们做如下规定：
（1）
把该式子求解得， （2）
将相关数据带入化简可得：可得 和t的关系式如下：
图3
图4
可以看出整个过程是相当的复杂，不仅过程多，而且在相同的过程中，不同的物质有着不同的结果。此问题中没有明确给出输出湖中的污染物的种类，也没有对湖泊环境做任何描述，无疑给问题的解决增加了极大的难度。为使问题简化，我们从一般情况出发，假设：湖中所进行的反应均为一级反应，有机物的存在不会对环境参数造成改变，环境固定（自然净化速率恒定），考虑湖泊生态系统中起主要作用的几种过程作简要分析。
针对污染源的不同污染和净化的情况，我们把该水域分为两个部分，湖水和河水。这两部分水的成分只在湖水的连接区域保持联系。
鉴于河水的流动性，水质轮新的速度较快，所以我们可以做这样的定义：其中湖水作为净化的场所，而对于河水而言，其功能上就是作为PC A公司投放污染源的场所。所以我们考虑对于污染源进入该流域内的速度，根据题意直接所得的相关表示形式。
1.物理沉降。不同物质有着不同的沉降速度常数，沉降速率为
2.挥发。不同物质有着不同的挥发速率常数 ，在有机物在水体上的大气中的
分压为0的条件下，挥发速率为 (Z为水体深度)
3.水解反应。不同物质有着不同的水解速率常数 ，在一级反应的条件下，水解速率为 ；
4.生物降解反应。不同物质有着不同的降解速率常数 ，在一级反应的条件下，生物降解速率 ；
5.2对于题设条件下对于水质净化作用的分析
而对于该水域内自带的水质净化作用，我们初步定义，对于水质的净化作用需要考虑该水域的面积和深度（即容积的大小）等因素。对于该题而言，该流域的净化作用能力主要体现在流出过程中的排污能力和多种形式的自我净化作用能力。
针对水质的自我净化作用和污染作用这两个过程，我们考虑是一个相互抗衡的阶段，根据题设已知污染源进入水质的速度恒定，我们考虑得，水质的净化作用速度收到该水域环境下的污染浓度限制，这样的过程是一个负反馈的过程。由基本假设得到，只有当达到污染源的浓度超过某个阶段，也就意味着水质的净化作用相应有一定幅度地减少。而针对在无污染的情况下，水质的净化作用随着水质的逐渐改善不断提高。
1
3
5
8
10
1.73
4.35
6.14
7.82
8.51
年份
15
20
30
40
50
9.43
9.78
9.97
9.995
9.999
表1
根据模拟数据所作的湖水污染浓度变化趋势如下图2所示
图2
（1）考虑湖水的自我净化作用
自我净化作用的原因：
根据资料显示：湖水中的污染物可分为有机污染物和无机污染物两大类，在多种环境因素（阳光，空气，水，水中生物，水中化学物质，重力等）的作用下，通过物理沉降，化学反应和生物转化一系列复杂的活动，它们的量会发生变化。有机污染物在水环境中的迁移过程过程如下图3，4所示。
根据微积分的知识可知，在适当短的时间段之内，通过建立微分方程，可以将连续的过程离散化，从而可以得到湖水污染浓度与时间之间的关系式。
5.3对于净化标准划定的分析
根据题意已知，当该水域的污染程度达到0.001mol/l，才达到需要治理水质的要求，这是一个基本红线。
5.4对于污染源划定的基本分析
污染源的划定，我们假设其作为含有某系列物质的溶液，其主要分为含有不同污染能力的固态物质若等，对于题三的整改部分，我们可以对公司提供一定的建议。鉴于公司不能直接切实有效对污染源的量做到减少，只有在排放的部分中做一定改善。对于不同污染能力的污染源物质，可以采取固态物质粉碎颗粒化，对于溶液进一步稀释。
问题：
1.在花费时间和经费去测试之前，建立数学模型用PCA提供的公开数据判断湖的环境是否会恶化。
2.派出野外工作人员测得目前湖水污染浓度为0.03mol/l，再测得河水污染浓度为0.05mol/l。以新数据为依据考虑湖水污染问题的数学模型。
3.现在假设你是环保局的所聘请的高级顾问，请向你的雇主提交一份报告.内容包括:
观察以上各式可发现：在一定条件下，各个速率与湖中污染程度呈线性关系，又因为某种化学物的总消失速率等于各过程消失速率之和，故我们可以认为：自然净化速率 （ 为正常数，与环境及污染种类有关）
如上所述，我们可得上述表达式的改进。
（3）
把该式求解得
（4）
当时间t= 时， ，即湖水污染浓度稳定于一个与K有关的值。
的大小影响着 达到最终稳定浓度， 的大小影响着 达到最终稳定浓度的快慢。在前面达到线性关系的基础上，在 =0.001 的c值大小确定A至数量级，并稳定年限在10-20年间，从而得到一股估计值。
= ， 1
3
5
8
10
1.493
2.971
3.539
3.816
3.871
15
关键字：时间步长法微分方程模型logistic模型
二问题重述
Pure河是流入Pristine湖的唯一河流。50年前PCA公司在此河旁建起一个生产设施并投入运行。PCA将为处理的废水排入河中,导致了Pristine湖被污染。PCA公司声称:已排放的废水的标准多年从未改变且不会对湖的环境有影响。
Pristine湖的湖容量为 ，流入(流出)的水流速度为 年。PCA声称河水污染浓度仅为0.001mol/l，自工厂开工以来没有改变过。
同时建立了计算机模拟模型，带入具体数值，运用时间步长法来仿真模拟了在湖水污染浓度稳定以前湖水每天的变化情况，输出自PCA建厂以来每年的湖水污染浓度变化，得到与微分方程模型相同的结论。
在全停产和半停产时，通过前面的两个模型可以计算湖水污染浓度在自然进化影响下的恢复到自然净华指标所需年限。并可得到结论：在半停产条件下，在选定的自然净化速率常数的约束下，只有当河水污染浓度降至原来的3.15%（自然净化速率成线性关系），4.7%（自然净化速率呈指数关系），才有可能使河水在100年内恢复至0.001 ，然后给出整改建议。
20
30
40
50
3.903
3.906
3.906
3.906
3.906
表3
根据模拟所得数据所做的湖水污染浓度变化趋势如下图6所示：
图6
观察数据可知：无论自然净化速率c与湖水污染浓度 呈线性关系还是指数关系，当湖水污染浓度 =0.001 时， 都将稳定于一个小于0.001 的值，也即：湖的环境不会恶化。
问题二：
模拟程序见附件一，从后面的分析中我们可知：在此河水污染浓度恒定，无自然净化的最简单的情况下建立的模型是以后问题的基础，后面的问题只是改变条件或数据，其实质是不变的，故我们在此程序中加入了多个选择语句，在不同的条件或数据下执行不同的命令，从而用一个程序解决了全部的模拟问题。
模拟所得数据如下表1所示：
年份
根据实际情况，我们在此只考虑存在自然净化的情况。
5.在整个大环境下，河水和湖水在污染源的浓度不一致。
6.在时间允许的范围内，假设整个湖水容量不发生变化。
7.污染源属于什么态的物质
四符号说明和名词解释
河水污染浓度
湖水污染物浓度
V
湖泊容量
L
c
自然净化速率
流入（流出）的水流速度
t
从PCA建厂至考察时刻的时间段
天
五问题分析
5.1对于水域基本情况分析
由于题目假设可得，我们对于整个水质的大环境假设为一个容积不变的大容器。基于在该流域里，降水量和蒸发量保持相等，湖水流入量和流出量相等且一直不变，所以我们假设该水域环境下，其容积不发生变化。
综上所述，对于该题，我们得到如下基本分析。在前假设条件的基础之上，湖水容量不变，出河水外无其他的污染源，故我们可将湖泊作为一个封闭的生态系统，其简化的湖水被污染的动态过程为：受污河水输入湖水中，河水与湖水均匀混合，受污水质在湖水区域内进行自我净化，湖水从输出端流出湖泊。湖中污染物的量直接决定了湖水污染浓度，而污染物的量受到以下两方面因素的影响：1.河水的污染浓度及其流入速度（根据已知此速度不变），
2.湖水的自我净化速度，前者使其增加，后者使其减少。
问题一、二、三的实质都是要分析污染浓度的变化趋势，而关键就在于前一因素的不同。问题一中，河水污染浓度不变，恒为0.001mol/L；问题二中，河水污染浓度可能会变化，受PCA效益的影响而按一定规律波动；问题三中，在全停产或半停产的情况下，和硕污染浓度为0或减为问题二中的一部分。后一因素（自然净化速度）在三个问题中的作用都是相同的。
根据所得数据如下表2所示
年份
1
3
5
8
10
1.587
3.437
4.315
4.878
4.984
年份
15
20
30
40
50
5.087
5.103
5.106
5.106
5.106
表2
根据模拟数据所作的湖水污染浓度变化趋势如图5所示，
图5
C进一步考虑自然净化速率的影响。更贴切于实际的情况是，自然净化速率c与湖水污染程度 成指数关系，c随的增加而增加，但增加的速率会逐步减小，用关系式表达即可为： ， 为环境和污染程度种类有关的常数。将此关系式代入微分方程（3），得到一个新的微分方程，此方程无解。但是我们可以通过此关系式，建立计算机模拟求得数值解。
建立湖水污染浓度随时间变化的微分方程模型：设在极短时间dt，湖水污染浓度增加d ，在将湖水被污染这一连续动态过程简化为离散的瞬间静止状态（如问题分析中所述）之后，根湖中剩余量=输入量—输出量，我们可以列微分方程如下
我们做如下规定：
（1）
把该式子求解得， （2）
将相关数据带入化简可得：可得 和t的关系式如下：
图3
图4
可以看出整个过程是相当的复杂，不仅过程多，而且在相同的过程中，不同的物质有着不同的结果。此问题中没有明确给出输出湖中的污染物的种类，也没有对湖泊环境做任何描述，无疑给问题的解决增加了极大的难度。为使问题简化，我们从一般情况出发，假设：湖中所进行的反应均为一级反应，有机物的存在不会对环境参数造成改变，环境固定（自然净化速率恒定），考虑湖泊生态系统中起主要作用的几种过程作简要分析。
针对污染源的不同污染和净化的情况，我们把该水域分为两个部分，湖水和河水。这两部分水的成分只在湖水的连接区域保持联系。
鉴于河水的流动性，水质轮新的速度较快，所以我们可以做这样的定义：其中湖水作为净化的场所，而对于河水而言，其功能上就是作为PC A公司投放污染源的场所。所以我们考虑对于污染源进入该流域内的速度，根据题意直接所得的相关表示形式。
1.物理沉降。不同物质有着不同的沉降速度常数，沉降速率为
2.挥发。不同物质有着不同的挥发速率常数 ，在有机物在水体上的大气中的
分压为0的条件下，挥发速率为 (Z为水体深度)
3.水解反应。不同物质有着不同的水解速率常数 ，在一级反应的条件下，水解速率为 ；
4.生物降解反应。不同物质有着不同的降解速率常数 ，在一级反应的条件下，生物降解速率 ；
5.2对于题设条件下对于水质净化作用的分析
而对于该水域内自带的水质净化作用，我们初步定义，对于水质的净化作用需要考虑该水域的面积和深度（即容积的大小）等因素。对于该题而言，该流域的净化作用能力主要体现在流出过程中的排污能力和多种形式的自我净化作用能力。
针对水质的自我净化作用和污染作用这两个过程，我们考虑是一个相互抗衡的阶段，根据题设已知污染源进入水质的速度恒定，我们考虑得，水质的净化作用速度收到该水域环境下的污染浓度限制，这样的过程是一个负反馈的过程。由基本假设得到，只有当达到污染源的浓度超过某个阶段，也就意味着水质的净化作用相应有一定幅度地减少。而针对在无污染的情况下，水质的净化作用随着水质的逐渐改善不断提高。
1
3
5
8
10
1.73
4.35
6.14
7.82
8.51
年份
15
20
30
40
50
9.43
9.78
9.97
9.995
9.999
表1
根据模拟数据所作的湖水污染浓度变化趋势如下图2所示
图2
（1）考虑湖水的自我净化作用
自我净化作用的原因：
根据资料显示：湖水中的污染物可分为有机污染物和无机污染物两大类，在多种环境因素（阳光，空气，水，水中生物，水中化学物质，重力等）的作用下，通过物理沉降，化学反应和生物转化一系列复杂的活动，它们的量会发生变化。有机污染物在水环境中的迁移过程过程如下图3，4所示。
