20xx年數學建模a題優(yōu)秀論文-城市表層土壤重金屬污染分析-預覽頁
【正文】 1 .307** .436** Sig. (2tailed) .000 .000 .000 .000 .066 .000 .000 N 319 319 319 319 319 319 319 319 Pb Pearson Correlation .290** .660** .383** .520** .298** .307** 1 .494** Sig. (2tailed) .000 .000 .000 .000 .000 .000 .000 N 319 319 319 319 319 319 319 319 Zn Pearson Correlation .247** .431** .424** .387** .196** .436** .494** 1 Sig. (2tailed) .000 .000 .000 .000 .000 .000 .000 N 319 319 319 319 319 319 319 319 **. Correlation is significant at the level (2tailed). 圖（ 11）： 重金屬原始含量數據和海拔的相關系數矩陣 注： 利用 SPSS軟件 做 出了八種重金屬在不同海拔內的分布圖（見附錄 ）。 （二）根據空間分布圖、區(qū)域散點圖和主要重金屬含量土壤單項污染的指數進行分析： 對于 Cr 和 Ni ，在來源上關聯較密切，該市表層土壤 Cr 和 Ni 基本未污染，只有個別點富集程度較高，污染達到中度污染，該富集中心的位置主要分布在生活區(qū)周邊和 交通 區(qū)周邊，這可能是由于生活廢水的排放和交通源汽車尾氣的排放等原因造成的。 對于 Hg , 其高含量點主要分布在交通繁忙的 交通 區(qū)周邊和工業(yè)區(qū)周邊 , Hg 污染的一個主要原因是由于燃煤造成的，無論是工業(yè)用煤還是居民用煤，而且燃燒方式落后。 綜上所述，可以認為工業(yè)區(qū)、 交通 區(qū)和生活區(qū)的活動是造成該城區(qū)表層土壤重金屬污染的主要原因。 模型的優(yōu)缺點 優(yōu)點： （ 1） 解決問題一的第一小問時，我們用 MATLAB 對原始數據進行差值擬合。 （ 2） 解決問題二時，我們用 SPSS對各種重金屬元素濃度和海拔做因子分析，得出各種元素濃度和海拔相關性，各元素濃度和海拔呈現負相關，正好驗證了第一問中求得的山區(qū)各重金屬濃度最低，污染程度最輕這一結果。 模型建立 為了更好地研究城市地質環(huán)境的演變模式，需要做以下幾個工作。 （ 1）高斯?jié)舛饶Ｐ徒? 本題假設污染物濃度符合正態(tài)分布，建立高斯?jié)舛确植寄Ｐ停? ? ? ? ? ???????? ????? 22222223212,zyxzyutxzyxeQzyxC ??????? （ 11） 其中 ? ?zyxC , 是空間任意一點 ? ?zyx , 的濃度； 0Q 是污染源重金屬產生率；amp。 （ 2）高斯模型的修正 因為干沉積和腐蝕是影響重金屬元素在土壤中的衰減的重要因素，故引入其濃度隨時間的變化關系： uteCC ??? 0 （ 12） 其中 v 是衰減系數 則某時刻下空間點 ? ?zyx , 處的重金屬濃度為： ? ? ? ? utzyutxzyxeeDQzyxC zyx ????????? ?????? 22222223212,???????, （ 13） 六．模型的評價與推廣 對于問題一所建模型，我們在求各區(qū)污染程度的時候，僅僅拿各區(qū)重金屬的平均濃度進行分析，得出的結果只能反映各區(qū)的平均污染等級，不能反映各區(qū)在不同位置的污染等級。 對于問題三所建模型，我們只考慮到海拔與空間距離，因此擬合的效果只能差強人意。 [3] 景麗潔，王敏 ，《 不同類型土壤對重金屬的吸附特性 》 , 《生態(tài)環(huán)境》 ,2020，第 17 卷第 1 期 ： 245248。A4:E32239。z=A(:,4)。)。) figure,contourf(X,Y,Z) figure,contour(X,Y,Z) title(39。,1,39。y=A(:,3)。 y=y(x1(:))。 x=A(:,2)。 y=y(x2(:))。 x=A(:,2)。 y=y(x3(:))。 x=A(:,2)。 y=y(x4(:))。 x=A(:,2)。 y=y(x5(:))。功能區(qū)分布 39。工業(yè)區(qū) 39。交通區(qū) 39。F:\A\cumcm2020A184。)。 nx=linspace(min(x),max(x),100)。v439。 hold on d=data x=d(:,1)。 ny=linspace(min(y),max(y),100)。)。string39。userdata39。 p4=[ ]。 p2=[ ]。 b=[ 130/1000 31 35 31 69 ]。 w4=p4./b。 p1z=sqrt((p1p^2+p1m^2)/2) p2p=sum(w2)/8。 p3z=sqrt((p3p^2+p3m^2)/2) p4p=sum(w4)/8。 p5z=sqrt((p5p^2+p5m^2)/2) As分布圖 Cd分布圖 0100200300400500600700As (μg/g) Cd(ng/g)Cr (μg/g) Cu (μg/g) Hg(ng/g)Ni (μg/g) Pb (μg/g) Zn (μg/g) 生活區(qū) 工業(yè)區(qū) 山區(qū) 交通區(qū) 公園綠地區(qū) Cr分布圖 Cu分布圖 Hg 分布圖 Ni 分布圖 Pb 分布圖 Zn 分布圖
公司管理相關推薦
鄂ICP備17016276號-1