邢利英　張國珍　孫三祥　饒小波

摘要：采用Landweber迭代建立了一維河流水體污染物非恒定輸運逆問題模型，通過某時刻實測值來反演污染物初始濃度分布，應(yīng)用Landweber迭代算法反演了一維純擴散和對流一擴散過程兩個算例。結(jié)果表明：Landweber迭代能準(zhǔn)確重構(gòu)污染物的初始分布，并且該方法具有良好的抗干擾性；當(dāng)測量誤差較大時，該算法仍能較好地反演污染物初始濃度，即Landweber達代可以較精確地反演一些非連續(xù)的污染物初始濃度，可以較好地解決一些連續(xù)污染物初值和污染源項的識別問題。

關(guān)鍵詞：純擴散方程；對流-擴散方程；Landweber迭代；污染物；初值

水污染模型是環(huán)境動力學(xué)中重要的研究課題。一些有機物質(zhì)、有毒有害物質(zhì)、重金屬的輸運過程符合對流一擴散方程，關(guān)于對流一擴散方程的反問題，有許多不同的算法，比如基于馬爾科夫鏈的貝葉斯法（B-MCMC）、Tikhonov正則法（TRM）、遺傳算法（GAM）和時空全域配點法（MQCM）等。Hazart A，等應(yīng)用B-MCMC重構(gòu)了點污染源的位置、釋放時間和釋放量：李功勝等應(yīng)用一種最佳迭代正則化算法反演污染源項，數(shù)值計算結(jié)果表明，在測量數(shù)據(jù)存在干擾時，這種算法仍可以精確、穩(wěn)定地反演污染源項；徐波等提出應(yīng)用GAM來確定源項，將源項識別問題轉(zhuǎn)變?yōu)閮?yōu)化問題，從任意的初始猜測值開始，通過交叉變異算子較為精確地尋找全局最優(yōu)解：關(guān)于初始條件的反演問題，李子等采用時空全域配點法求解一維污染物對流一擴散過程的源項識別問題。

上述數(shù)值方法反演速度比較快，能夠滿足工程的要求，但不可避免地存在一些問題：B-MCMC需要估計一些參數(shù)的先驗概率分布，若這些估計不準(zhǔn)確，則可能導(dǎo)致較大的誤差：TRM需要一個較好的近似解，不斷驗算正則化參數(shù)，以便得到一個較為合適的重構(gòu)值；針對GAM反演源項，需要編制復(fù)雜的程序，包括選擇算子、交叉算子、變異算子，這3個算子又包含交叉概率和變異概率等參數(shù)，這些參數(shù)對于模擬結(jié)果有較大影響。目前來說，這些參數(shù)的選擇主要依靠經(jīng)驗，若選擇不合適，則可能需要更多的訓(xùn)練時間搜尋最佳結(jié)果。Landweber迭代不存在以上問題，汪繼文等應(yīng)用Landweber迭代重構(gòu)了基于時間和空問變量的熱傳導(dǎo)方程的源項。基于上述分析，本文采用Landweber迭代重構(gòu)非恒定河流污染物輸運的初值。