最新湖水污染分析模型

1、摘要在兩種情況下分析湖水中的污染物,分別建立模型即理論模型和實際模 型.理論模型是根據(jù)伊利湖和安大概湖各自的污染物流入流出的關系建立污染 物量關于時間的差分方程：伊利湖的污染物總量an+1 0.62an安大概湖的污染 物總量 bn 6129.0323 0.62n 7020.3360 0.87n 192.3077bn 在 n 時趨于一個定值192.3077,這個定值就是安大概湖系統(tǒng)的平衡值；當n 35時 bn 245.95安大概湖的污染程度減少到目前水平的10%;當3 n 1是系統(tǒng)的污染物的量是一直增加的,當20 n 3系統(tǒng)的污染物量急劇減少,大約從n 40 開始系統(tǒng)的污染物量幾乎保持不變.實際

2、模型中首先根據(jù)湖水的實際更新情況 重新確定湖水流入和流出占湖水總量的百分數(shù),又由于湖水中污染物的濃度時 刻變化,所以用時間微元的方法對實際污染物流出的比例進行修正.分析鋁廠排放的污染物時,鋁廠排放的污染物是赤泥,根據(jù)赤泥的物化性質(zhì)利用重力沉 降原理求得赤泥顆粒從湖面沉降到湖底的時間t,把一年分成多份t同時將鋁廠每年向湖水中排放的污染物量 25單位按t分成多份,每一個單位時間鋁廠排 放到湖里的污染物量是q 0.3單位,那么安大概湖的湖水中將始終保持有0.3單位的赤泥,其余的赤泥都將在湖底沉積.綜合安大概湖中赤泥和伊利湖流入的 污染物的情況預測了未來十年內(nèi)的情況.模型中重力沉降原理指出顆粒的直徑

3、影響沉降速度間接影響赤泥的排出量直徑越小排出量越大,同時直徑是最可能實現(xiàn)改良的因素.在直徑小于20um時赤泥的排出量急劇增加.為減少安大概 湖的污染盡量把顆粒直徑做小.、問題分析專業(yè)技術資料伊利湖的湖水每年有38%的更新,湖水的更新引起湖內(nèi)污染物量的變化假設流入伊利湖的湖水是不含有污染物的,而流出伊利湖的湖水又將攜帶污染 物,那么伊利湖是一個沒有污染物注入只有污染物排除的系統(tǒng) ,污染物的量逐 漸減少,根據(jù)污染物排除的情況獲得伊利湖污染物量隨時間變化的關系.安大 略湖湖水每年更新13%,在認為流入湖的水僅來自伊利湖的情況下,安大概湖 的污染物源只有伊利湖和鋁廠,根據(jù)每一年的湖水更新量和鋁廠的排放

4、量求一 二階差分方程,用來描述安大概湖內(nèi)污染物伊利湖流入的和鋁廠排放的和 時間變化的情況.但是這些都只是很簡單的處理問題,考慮實際情況時.要有 很多因素影響著伊利湖和安大概湖的污染物總量,湖水的更新不僅僅包括湖水 的流入和流出還有將于和蒸發(fā)等；鋁廠排放的污染物在湖水中的性質(zhì)與伊利湖 流入的污染物不同；湖水內(nèi)污染物的量時刻變化等實際情況建立更加合理的模 型.三、理論模型建立1.模型假設1 .污染物在湖水中均勻分布,每一處污染物濃度相同2 .對理論模型假設每年安大概湖的湖水補充僅來自伊利湖的流入.3湖水的更新量中湖水的增加僅來自湖水的流入量,減少量僅是湖水的流出即忽略湖水的蒸發(fā)和降雨等.4.安大概

5、湖的湖水流出速度恒定.2.符號說明專業(yè)技術資料n :代表年份取值從1開始,初始的那年為第1年.an：第n年時的伊利湖的污染物總量bn：第n年時的安大概湖的污染物總量3,差分方程的建立和求解伊利湖每年的湖水更新量是38%,流入的的水流中沒有污染物,流出的湖水攜帶污染物.安大概湖的湖水補充僅來自伊利湖,每年伊利湖湖水總量的38%經(jīng)尼亞加拉河流入安大概湖,安大概湖的污染物消失都來自湖水的流出貢獻.安大概湖的污染物總量用函數(shù)b (n)表示,伊利湖的污染物總量用函數(shù)a(n)表示.根據(jù)伊利湖水的污染物總量的變化可得如下方程：an+20.62烝+1an+10.62ann 1,2,3.根據(jù)安大概湖的污染物情況

6、列出可以表示污染物變化關系的函數(shù)式25bn 2 0.38an 1 bn 1 0.13bn+1bn 1 0.38an bn 0.13bn 25n 1,2,3.把an和an 1消去得到的二階差分方程如下bn 2-1.49bn 1 0.62 0.87bn25 0.38可得上式的齊次方程：2 1.490.620.87 010.62 ,20.87.專業(yè)技術資料齊次方程的通解為：bn Ci 0.62n C2 0.87n 設 1 式 的 通解為： bn Ci 0.62n C2 0.87n b式代入1式是滿足：bn+2 1.49bn1 0.62 0.87bn 0式代入1式滿足bn+2 1.49bn i 0.6

7、2 0.87bn 25 0.38.根據(jù)這樣二個規(guī)那么 可求出 b 的方程：b 1.49b 0.62 0.87b 25 0.38,解得 b=192.3077,貝U 1 式的通解3變?yōu)閎nC1 0.62n C2 0.87n 192.3077.第一年安大概湖的污染物總量是2500單位即b1 2500 ,第二年污染物總量增加至3150單位 b2 3150 ;通過 MATLAB 編程可求得 C16129.0323, C2 7020.3360 ,所以二階差分方程的通解可以代表安大概湖的污染物總量的方程為：bn 6129.0323 0.62n 7020.3360 0.87n 192.307744 .系統(tǒng)的平

8、衡值安大概湖污染物總量 bn 6129.0323 0.62n 7020.3360 0.87n 192.3077 當 n從1開始增長下去,6129.0323 0.62 n 7.2.336.0.87n逐漸趨近于0所以bn 隨著n的增長逐漸趨近于192.3077.通過方程可以證實系統(tǒng)存在平衡值是192.3077,隨著時間的推移安大概湖內(nèi)的污染物的流入和排除水平相同,湖內(nèi)的污染物總量保持在192.3077單位.其實,通過實際情況的分析亦可得到相同的 結論,安大概湖的污染物僅來自鋁廠和伊利湖,鋁廠每年的排入量恒定25 單位；伊利湖的污染物每年流出總量中的一局部,而且沒有外界的補充,這樣時間的增長使伊利湖

9、內(nèi)的污染物總量越來越少趨近于零,伊利湖流入安大概湖專業(yè)技術資料的污染物自然減少至零從某一年開始安大概湖的污染物僅來自于鋁廠,并且流入和流出的水平相同,所以說每年鋁廠流入安大概湖的25單位的污染物的同時,安大概湖也要通過每年13%湖水的更新中把25單位的污染物排走,污染物再湖水中平均分布,鋁廠流入的25單位的污染物相當于安大概湖污染物平衡時的污染物總量的13%.根據(jù)這個關系我們可以分析出安大概湖的污染物平衡值是192.3077.通過實際分析出來的平衡值和方程顯示的平衡值相統(tǒng)一.5 .污染物減少到目前水平10%的時間和系統(tǒng)的長期表現(xiàn)目前污染物總量為2500單位,減少到10%即250單位左右.對污染

10、物總量和時間繪二維圖形如圖1所示.通過MATLAB計算由表1可知道n 34時污染物總量b34 253.95, n 35時b35 245.95.由這些數(shù)據(jù)可以知道現(xiàn)在作為第一年,要等到第35年污染物總量才能低于現(xiàn)在水平的 10%.表1.年份-污染物總量單位年份污染物量污染物減少量年份污染物總量污染物減少量年份污染物總量污染物減少量1250035245.958.0269192.7780.0723150 .-65036238.9786.9770192.7170.0633354.50-204.537232.916.0771192.6640.0543308.545.9138227.6335.287219

11、2.6170.05專業(yè)技術資料953129.88178.7139223.0414.5973192.5770.0462888.37241.5140219.0453.9974192.5420.0372624.92263.4541215.5693.4875192.5120.0382362.64262.2842212.5453.0176192.4850.0392113.95248.6943209.9142.6577192.4620.02101884.88229.0744207.6252.2878192.4420.02111677.70207.1845205.6342.0079192.4240.0181

12、21492.57185.1346203.9011.7380192.4090.015131328.48164.0947202.3941.5181192.3960.013專業(yè)技術資料141183.84144.6448201.0831.3182192.3840.012151056.85126.9949199.9421.1483192.3740.0116945.636111.2150198.9500.9984192.3660.00817848.43497.2051198.0860.8685192.3580.00818763.59084.8452197.3350.7386192.3510.0071968

13、9.60473.9953196.6810.6887192.3460.00520625.13063.9754196.1130.5788192.3410.00521568.97156.1655195.6180.4989192.3360.00522520.07248.956195.1870.4390192.3320.00423477.5042.5757194.8130.3791192.3290.003專業(yè)技術資料424440.45437.0558194.4870.3392192.3260.00325408.21132.2459194.2040.2893192.3240.00226380.15328.

14、0660193.9570.2594192.3220.00227355.73924.4161193.7430.2195192.3200.00228334.49721.2462193.5560.1996192.3180.00229316.01518.4863193.3940.1697192.3170.00130299.93416.0964193.2530.1498192.3160.00131285.94413.9965193.1300.1299192.3150.00132273.7712.1766193.0230.11100192.3140.00133263.1810.5967192.9300.0

15、9專業(yè)技術資料34253.9689.2168192.8490.07把上表的數(shù)據(jù)繪制成二維圖形可以直觀的看出污染物總量和污染物減少量隨時間變化的趨勢,從圖1中可以明顯的看出系統(tǒng)的長期情況安大概湖的 污染物總量僅在前三年的范圍內(nèi)值增長的,從第四年開始污染物總量的值開始下降,從第四年開始的前20年下降速度非常的快,圖2也說明了安大概湖系統(tǒng)的污染物量逐漸趨于穩(wěn)定,從第四十年開始污染物變化量幾乎為 00圖1.污染物總量與時間的關系圖2污染物變化量與時間關系實際模型建立符號說明：A實際情況下第n年時的伊利湖的污染物總量專業(yè)技術資料B實際情況下第n年時的安大概湖的污染物總量.:修正后的伊利湖污染物每年排除的

16、百分數(shù),可直接表示污染物變化情況.:修正后的伊利湖每年排到安大概湖污染物的百分數(shù),可直接表示伊利湖排到安大概湖的污染物情況.:修正后的安大概湖污染物每年排放的百分數(shù),可直接表示污染物變化情況.模型假設1 .鋁廠向安大概湖排放的污染物中只有赤泥,而且赤泥顆粒在湖水中的運動和 分布相互獨立.2 .對于安大概湖除伊利湖的流入含有污染物外,降雨和其他支流都不含有污染 物.3 .由于安大概湖的平均水深較深,在污染物赤泥沉降到湖底以后,只能在 湖底層附近移動,不能浮起.4 .設污染物赤泥在湖中沉積過程中以勻速下沉,其速度符合重力沉降原 理.5 .湖水的密度為正常水的密度即為1 103kg/m3.6 .湖底

17、初始時已有的赤泥的量為00安大概湖是安大概省西端中最小的五大湖,平均深度86米283英尺,流入專業(yè)技術資料安大概湖的水中79%是伊利湖通過安大概省西端尼亞加拉河流入的,剩下的流動支流來自安大概省西端盆地(14%)和降水(7%)0 93%的湖水從安大概省西端 的圣勞倫斯河流出,其余的7%是蒸發(fā)失去.安大概湖的總蓄水量是1688立方 千米,伊利湖的總蓄水量是1638立方千米.安大概湖的污染物來源主要是伊利 湖湖水的流入和鋁廠廢物的排放.鋁廠排放的廢物主要是赤泥,赤泥是從鋁土礦 中提煉氧化鋁后排出的工業(yè)固體廢物.赤泥的化學成分為：FezO3、Si.?、CaO AI2O3、MgO、Na2O >

18、K2O > TiO2;赤泥的物 理性質(zhì)是：顆粒直徑 0. 080. 25mm .1 .伊利湖的污染物模型對于伊利湖,伊利湖的污染物總量沒有外界的補充,污染物隨著伊利湖湖水的流出逐漸減少.每年的湖水更新量是38%,其中補充的水量包括降雨(7%)和河流流入,他們都不包含污染物,排除的水量包括河流的流出和蒸發(fā)其占所有的更新量的(7%);所以每年湖水通過 尼亞加拉排除的更新量是用 38% (1 7%) 35%.由于污染物的量隨時間連續(xù)變化,所以不能用35%水流量的更新直接表示污染物的排除,對35%作如下修正.將一年的時間分為m份,每一份時間內(nèi)流出的水流量占總水流量的0.31/m,那么可以知道每一

19、份時間過后湖中剩余的污染物量dm：d1 d0 0.35d0/m = (1 0.35/m) d0d2 d1 0.35dl. m = (1 0.35. m)2dodm dm 1 0.35dmi.m = (1 0.35. m)md0當m趨向于無窮時dm e0.35d0,即dm 0.705C0所以35%修正后的結果為專業(yè)技術資料1 0.70529.5% o每年伊利湖中污染物的量減少 29.5%.所以伊利湖中污染物的量可以表示為：An nA0n=1,2,3 2 .安大概湖污染物量對于安大概湖的情況污染物的來源只有伊利湖的流入和鋁廠的排放,研究安大概湖的污染物情況時把伊利湖流入的污染物和鋁廠排放的污染物分

20、開考慮 分別建立模型.2.1源自伊利湖的污染物通過網(wǎng)上資料總結得到以下湖泊河流以及降水蒸發(fā)之間的關系見圖3用出一專業(yè)技術資料圖3.伊利湖與安大列湖的關系每年伊利湖流入安大概湖的水流量占安大概湖湖水更新局部(13%)的79%,這是由于安大概湖的79%的水是來自伊利湖的.根據(jù)這個關系可以知道 每年從伊利湖流入安大概湖的水流量 Q=13%*79%*1688=175.552 立方千米, 這局部占伊利湖湖水總量的10.7%也就是說每年伊利湖的污染物通過這 10.7% 湖水更新流入安大概湖.由于伊利湖中湖水的流出是連續(xù)的,湖水中的污染物 的量也是連續(xù)變化的,因此在求解污染物每年流出的量的百分數(shù)時不能直接使

21、 用10.7%,對10.7%修正,修正方法同伊利湖模型中的相同.將一年的時間分 為m份,每一份時間內(nèi)流出的水流量占總水流量的0.107/m,那么可以知道每一份時間過后湖中剩余的污染物量Pm：P1 P0 0.107p0/m = (1 0.107/m) p02p2 p1 0.107Pl/m = (1 0.107/m) P0Pm Pm 1 0.107Pm1；m = (1 0.1071)%當m趨向于無窮時Pm e 0.107p0 ,即pm 0.8985p01 0.89850.1015所以10.7%經(jīng)修正后的結果是10.15% ,可以直接描述每年伊利湖中排放到安大概湖的污染物量.安大概湖的湖水每年更新1

22、3%,而更新的湖水93%來自湖水的流出,7%是 湖水的蒸發(fā),湖水的蒸發(fā)對污染物的排除是沒有奉獻的.所以每年安大概湖的 湖水通過圣勞倫斯河流出更新是 13%*93%=12.09% .安大概湖污染物的量也是專業(yè)技術資料隨著時間連續(xù)變化的,通過上面提到的百分數(shù)的修正q0 0.1209q0/m = (1 0.1209/m) q()2q2q1 0.1209q1; m = (1 0.1209； m) q0qm 1 0.1209qm 1/m = (1 0.1209 m)mq0當m趨向于無窮時qm e020%.,即qm 0.8861q那么每年安大概湖排除的污染物的百分比修正后的結果是.1 0.886111.3

23、9%根據(jù)安大概湖, , 0污染物流入和排除關系可得到差分方程：Bn 2 0.1015An1 Bn 1 0.1139Bn+1Bn 1 0.1015An Bn 0.1139BnA n+1(1+Bn 1(1) Bn+1(1) 0可得上式的齊次方程：方程的通解為：Bn& (1)C2 n安大概湖第一年的污染物總量是2500單位,第二年3150單位.從鋁廠污染物的計算中可以知道,由于鋁廠排放的赤泥在湖水中溶解的十分的少絕大局部沉到湖底,我們假設安大概湖水體中的污染物都是伊利湖排放來的.通過B二2500和B2 3150可以求解出C1 =4605.99 , C2=-5360.59安大概湖中來自伊利湖的

24、污染物隨時間變化的方程：Bn 4605.99 (1)n-5360.59 n近十年的情況如圖2所示.通過圖4可以觀察到,那些源于伊利湖的污染物在安大略湖中從第三年開始呈 線性下降趨勢.由于伊利湖內(nèi)的污染物逐漸減少,排入專業(yè)技術資料到安大概湖的污染物自然也減少,照這個趨勢一致走下去,伊利湖和安大概湖內(nèi)來自伊利湖的污染物都將減少至 00圖4.安大概湖中來自伊利湖的污染物變化情況對應上圖的的未來三年的數(shù)據(jù)如表2所示.在第五年污染物的總量又回到了初始的水平.表2.未來三年由伊利湖流入安大概湖的污染物的預測年份第一年第二年第三年第四年第五年污染物量250035103067279925042.2鋁廠排放的赤

25、泥考慮安大概湖的污染物實際情況,實際情況中污染物除了伊利湖的湖水流入攜帶污染物外,僅有鋁廠向安大概湖排放鋁廠的廢物.通過查詢資料可知鋁赤泥的化學物理組成決定了廠像安大概湖排放的廢物是赤泥,因此只需要知道赤泥在湖水中的存在情況就 可以描述鋁廠的排放對安大概湖污染物量的影響專業(yè)技術資料赤泥不能像一般污染物一樣溶解在湖水中,并在湖水中擴散,隨湖水的流出而排除湖水系統(tǒng).實際上赤泥以顆粒的形式在湖水中向下沉降,絕大多數(shù)的赤泥會沉降到湖底,不能夠隨著水流排除到圣勞倫斯河中,為驗證這一情況,下面 做出一些討論.符號說明t:赤泥顆粒從湖頂沉積到湖底的時間.:湖水的密度水的密度.Ut:污染物顆粒在水中的平均速度

26、.d ：赤泥顆粒的直徑0.16mm.S ：赤泥的密度.:水的粘滯系數(shù).g :重力加速度.h:安大概湖的平均水深.q :赤泥顆粒從湖面沉到湖底的時間內(nèi)鋁廠向湖內(nèi)排放的赤泥的量Q :赤泥每年隨湖水排除的污染物的量.Qn：第n年時湖底積累的赤泥的量.Wn：第n年安大概湖總的污染物量.Re-沉降雷諾數(shù).專業(yè)技術資料赤泥顆粒在水中運動時的速度非常的小,所以忽略到加速段的影響,認為赤泥顆粒在水中沉降的速度是勻速的重力沉降原理模型：當K>69.1時為湍流區(qū)Newton牛頓區(qū)下限段；當2.62<K<69.1時為過渡區(qū)Allen阿倫區(qū)；當K<2.62時為滯留區(qū)Stokes斯托克斯區(qū)算得K

27、1.6 10數(shù)群K值判斷法：K4 1 103(2.8 103 1 103) 9.8:(0.839 103)24.71根據(jù)K =4.71 ,確定赤泥顆粒的沉降在重力沉降原理模型中的過渡區(qū)那么有RetduUt 0.27d s -Ret0.6-4 、求得Ret 4.27ut 2.24 10 m/s安大概湖的平均水深為86米,那么赤泥從湖面沉降到湖底的平均時間,h 86t 4ut2.24 10 43.84 105(s)將t的單位換算成年：t一一 _ 53.84 10365 24 60 600.012( year)把t作為一個單位時間,那么一年一共由83.3個單位時間組成；鋁廠一年向 安大概湖龐芳的污染

28、物的總量為25單位,那么每一個時間t內(nèi)排放到安大概湖25的赤泥的量為q q,q 0.3單位.由于每一個單位時間是顆粒從湖面沉 83.3q 0.3 單降到湖底的時間,每一個單位時間鋁廠排放到湖里的污染物量是專業(yè)技術資料位,那么湖水里將始終保持有0.3單位的赤泥,隨著時間推移,排入到安大概湖的赤泥越來越多,湖水中始終保持0.3單位的量,更多的赤泥就在湖底沉積.沉 在湖底的赤泥在沒有外界清楚地情況下降一直在湖底積累,湖水中含有的0.3單 位的赤泥可以隨著湖水的流出而局部排除.安大概湖的湖水每年更新13%,其中93%是安大概湖湖水的流出,其余是湖水的蒸發(fā),湖水的蒸發(fā)對污染物的排除沒有作用,所以在湖水中

29、污染物始終保持0.3單位的情況下每年有Q q 13% 93%, Q 0.036單位,每年鋁廠排放的25單位污染物中只有0.036被排除,剩下的24.964單位都將沉積到湖底.到第n年時湖底積累的赤泥Qn(n 1)(25 Q) n =1,2,3 安大概湖的污染物量由伊利湖的流入和鋁廠的排放組成,安大概湖污染物總量與時間的關系函數(shù)式：Wn Qn Bn.未來十年的開展情況如圖5,對未來的十 年預測如表3中所示.We lutal 3irciTj!仃 pQlLrlant hiJ dTiturt n'附RWE fem entFigure 5.實際安大概湖污染物情況專業(yè)技術資料Time1234567

30、8910Thetotal2500.3175.3116.2873.2604.2350.2124.1925.1750.1599.amountof0099084082pollutantThe total024.9649.9274.8999.85124.8149.7174.7199.7224.6amount of red4262417mudThe total2500315030672779.2504.2226.1974.1750.1551.1374.amount of9206315pollutant fromErie表3.實際情況下安大概湖的污染物情況預覽未來十年赤泥在湖底逐年累積,流出局部非常的少而由伊利湖流入的污染物逐年減少,總污染物量的減少是伊利湖污染物減少的作用.模型結果的分析不同的數(shù)群判斷法值K對應不同的區(qū)域,K d s 2的“,其中當K>69.1時為湍流區(qū)Newton牛頓區(qū)下限段；當2.62<K<69.1時為過渡區(qū)Allen阿倫區(qū)；當K<2.62時為滯留區(qū)Stokes斯托克斯區(qū).專業(yè)技術資料換算成顆粒直徑為：當時污染物沉降處于滯留區(qū),當d 2.36時