确定污染源的位置

《确定污染源的位置》由会员上传分享，免费在线阅读，更多相关内容在应用文档-天天文库。

1、5污染源的判定在问题（2）中，通过考虑不同的条件，建立了三个模型，首先单独作分析，然后综合三个模型得出结论，具体情况分析如下：1、差分方程反演模型的建立与求解1.1模型I的建立问题II研究和分析的是长江干流近一年多来受高锰酸钾盐指数及氨氮的污染源情况，建立该模型主要考虑以下几个因素：一、浓度是随时间、距离变化的量。二、长江中的水流量及水中污染物在消耗减少，同时也在补充，支流以及沿干流的排污都是其“外援”。三、长江自身具有自净能力，由于降解系数介于0.1~0.5之间，我们考虑取其值0.2。四、考虑河流在受污染时的扩散能力，其扩散系数为。根据以上考虑的因素及流体力学的基本理论，我们建立一个反应扩散

2、方程来反映长江流域污染物浓度的变化情况：5-1-1其中，为污染物质在长江内浓度；为该段河流水流速度；为外界输入的污染物浓度，根据现有数据资料，我们不妨先假设，在任意一个分段内为常数，即整个是一个分段的常数函数。为扩散系数，为降解系数。1.2差分方程反演模型对污染源的判定根据每一段已经给出的数据，我们确定每一小段的水流速度。首先考虑整条长江流域的平均流速　　5-1-2其中12是第个干流站点的水流速度。假设考虑权重，可以用于表示任一小段的平均水流速度，并且和整条长江的平均水流速度联系到起来，表示如下：　　5-1-3其中，。由于沿长江往下，水流速度在逐渐减少，在上游速度比较大时可以用后一个点的流速决

3、定该流水段的水流速度，到下游时可以将权重逐渐转移到了上，由此我们取，其中3215是干流上站点之间的总间隔。表示从上游往下权重的变化情况。依据上述方程5-1-3可以依次解出。用下表列出结果（表5-1-1）：表5-1-1各段的平均流速：第一段第二段第三段第四段第五段第六段第1月2.56751.53780.90.919321.01431.1第2月2.4261.38470.83470.938631.11.1第3月2.6261.67841.23471.33861.51.5第4月2.8261.87841.43471.51.51431.6第5月2.59681.77211.41.41931.52851.7第6

4、月4.88773.34771.76941.93862.28553.4第7月2.53381.92521.53471.61.61431.7第8月2.13381.33780.73470.819320.90.9第9月1.96751.12520.73470.80.80.8第10月1.67531.03150.53470.619320.70.7第11月1.29220.718920.469390.619320.70.7第12月1.19220.665760.469390.638630.80.8第13月1.4630.825220.43470.538630.714270.8从资料[1]中我们查得，扩散系数取。对微分

5、方程两边积分，积分区间为（），其中。12其中。对上式进行处理，含积分项利用梯形公式，含偏导数项利用一阶向前差商，得到的表达式：5-1-4其中，。当，时，时间偏导的差分用向后差商，得到的结果为：5-1-5当，时，扩散项对距离偏导的差分用向后差商处理，得到结果如下：5-1-61.3模型的结果与分析根据表达式5-1-4，5-1-5，5-1-6，利用Matlab编程代入数据分别求解出污染物高锰酸钾盐在各段流域的的值（程序见附录2程序1），由于表示的是长江水流中污染物的浓度，所以总是不小于零的数，故将其中负数归零处理（其它关于的数值同样处理），列入下表(表5-1-2)：12表5-1-2高锰酸钾盐的值钾盐

6、第一段第二段第三段第四段第五段第六段第1月0.871920.339570.7850.564790.526740.17498第2月0.169950.669680.78130.578860.653240.18947第3月0.672740.824140.650680.4135900.47018第4月0.836330.605761.0490.2570100.38149第5月00.771460.989130.20470.0860610.30163第6月0.159310.366510.92202000.83506第7月0.464330.955710.607950.0875170.576890.33663第

7、8月0.293070.371780.585910.440010.43340.32784第9月0.229630.49910.785280.442480.163450.43812第10月0.268560.398611.04470.440880.58720.49527第11月0.380780.399250.659980.612690.41070.4394第12月0.388420.421370.81066