20xx年數(shù)學(xué)建模a題優(yōu)秀論文-城市表層土壤重金屬污染分析-預(yù)覽頁
【正文】 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）： 重金屬原始含量數(shù)據(jù)和海拔的相關(guān)系數(shù)矩陣 注： 利用 SPSS軟件 做 出了八種重金屬在不同海拔內(nèi)的分布圖（見附錄 ）。 （二）根據(jù)空間分布圖、區(qū)域散點(diǎn)圖和主要重金屬含量土壤單項(xiàng)污染的指數(shù)進(jìn)行分析： 對于 Cr 和 Ni ，在來源上關(guān)聯(lián)較密切，該市表層土壤 Cr 和 Ni 基本未污染，只有個(gè)別點(diǎn)富集程度較高，污染達(dá)到中度污染，該富集中心的位置主要分布在生活區(qū)周邊和 交通 區(qū)周邊，這可能是由于生活廢水的排放和交通源汽車尾氣的排放等原因造成的。 對于 Hg , 其高含量點(diǎn)主要分布在交通繁忙的 交通 區(qū)周邊和工業(yè)區(qū)周邊 , Hg 污染的一個(gè)主要原因是由于燃煤造成的，無論是工業(yè)用煤還是居民用煤，而且燃燒方式落后。 綜上所述，可以認(rèn)為工業(yè)區(qū)、 交通 區(qū)和生活區(qū)的活動是造成該城區(qū)表層土壤重金屬污染的主要原因。 模型的優(yōu)缺點(diǎn) 優(yōu)點(diǎn)： （ 1） 解決問題一的第一小問時(shí)，我們用 MATLAB 對原始數(shù)據(jù)進(jìn)行差值擬合。 （ 2） 解決問題二時(shí)，我們用 SPSS對各種重金屬元素濃度和海拔做因子分析，得出各種元素濃度和海拔相關(guān)性，各元素濃度和海拔呈現(xiàn)負(fù)相關(guān)，正好驗(yàn)證了第一問中求得的山區(qū)各重金屬濃度最低，污染程度最輕這一結(jié)果。 模型建立 為了更好地研究城市地質(zhì)環(huán)境的演變模式，需要做以下幾個(gè)工作。 （ 1）高斯?jié)舛饶Ｐ徒? 本題假設(shè)污染物濃度符合正態(tài)分布，建立高斯?jié)舛确植寄Ｐ停? ? ? ? ? ???????? ????? 22222223212,zyxzyutxzyxeQzyxC ??????? （ 11） 其中 ? ?zyxC , 是空間任意一點(diǎn) ? ?zyx , 的濃度； 0Q 是污染源重金屬產(chǎn)生率；amp。 （ 2）高斯模型的修正 因?yàn)楦沙练e和腐蝕是影響重金屬元素在土壤中的衰減的重要因素，故引入其濃度隨時(shí)間的變化關(guān)系： uteCC ??? 0 （ 12） 其中 v 是衰減系數(shù) 則某時(shí)刻下空間點(diǎn) ? ?zyx , 處的重金屬濃度為： ? ? ? ? utzyutxzyxeeDQzyxC zyx ????????? ?????? 22222223212,???????, （ 13） 六．模型的評價(jià)與推廣 對于問題一所建模型，我們在求各區(qū)污染程度的時(shí)候，僅僅拿各區(qū)重金屬的平均濃度進(jìn)行分析，得出的結(jié)果只能反映各區(qū)的平均污染等級，不能反映各區(qū)在不同位置的污染等級。 對于問題三所建模型，我們只考慮到海拔與空間距離，因此擬合的效果只能差強(qiáng)人意。 [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 分布圖
公司管理相關(guān)推薦
鄂ICP備17016276號-1