【正文】 654917xy??? ??兩個(gè)點(diǎn)，對(duì)比前面的模型，誤差 為 ‰ Zn分布圖： 16 可得 138479589xy??? ??一個(gè)點(diǎn)，對(duì)比前面的模型，誤差 為 ‰ 由此可看出，模型的準(zhǔn)確性 稍差 。 （刮風(fēng)，下雨，沙塵 暴等）也會(huì)影響重金屬的擴(kuò)散。 c如果假設(shè)重金屬污染物在土壤深度 L之內(nèi)是均勻變化的，且在下部邊界沒 有污染物的流失，而土壤中污染物的揮發(fā)速率受土壤上部空氣去處揮發(fā)出得污染物速率的控制，即相對(duì)于流失而言，擴(kuò)散的影響可以忽略；如果空氣停止流動(dòng)，則污染物的揮發(fā)量為零。 [3] 朱旭，李煥琴，籍萬新， MATLAB 軟件與基礎(chǔ)數(shù)學(xué)實(shí)驗(yàn)，西安：西安交通大學(xué)出版 社， 2022。 八 、附錄 附錄一 ① 城區(qū)地形示意圖如下： 附錄二 ① 一區(qū)各元素的分布圖： 21 111 112 ② 二區(qū)各元素分布圖： 211 22 212 ③ 三區(qū)各元素分布圖： 311 312 ④ 四區(qū)各元素分布圖： 23 411 412 ⑤ 五區(qū)各元素分布圖： 511 24 512 附錄三 ① As在城區(qū)情況分布圖： 121 ② Cd在城區(qū)情況分布圖： 25 221 ③ Cr在城區(qū)情況分布圖： 321 ④ Cu在城區(qū)情況分布圖： 26 421 ⑤ Hg在城區(qū)情況分布圖： 521 ⑥ Ni在城區(qū)情況分布圖： 27 621 ⑦ Pb在城區(qū)情況分布圖： 721 ⑧ Zn在城區(qū)情況分布圖： 28 821 附錄四 %1區(qū) As分布圖 x=data(1:44,2)。)。 [xi,yi]=meshgrid(ti,tu)。 z=data(1:44,8)。o39。cubic39。 tu=0:10:18449。 y=data(1:44,3)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 zi=griddata(x,y,z,xi,yi,39。 ti=0:10:26453。) %2區(qū) As分布圖 x=data(45:80,2)。)。 [xi,yi]=meshgrid(ti,tu)。 z=data(45:80,8)。o39。cubic39。 tu=0:10:17980。 y=data(45:80,3)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 zi=griddata(x,y,z,xi,yi,39。 ti=0:10:22674。) %3區(qū) As分布圖 x=data(81:146,2)。)。 [xi,yi]=meshgrid(ti,tu)。 z=data(81:146,8)。o39。cubic39。 tu=0:10:18397。 y=data(81:146,3)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 zi=griddata(x,y,z,xi,yi,39。 ti=0:10:28654。) %4區(qū) As分布圖 x=data(147:284,2)。)。 [xi,yi]=meshgrid(ti,tu)。 z=data(147:284,8)。o39。cubic39。 tu=0:10:17949。 y=data(147:284,3)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 zi=griddata(x,y,z,xi,yi,39。 ti=0:10:26424。) %5區(qū) As分布圖 x=data(285:319,2)。)。 [xi,yi]=meshgrid(ti,tu)。 z=data(285:319,8)。o39。cubic39。 tu=0:10:17133。 y=data(285:319,3)。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 zi=griddata(x,y,z,xi,yi,39。 ti=0:10:20582。) %城區(qū)地形圖 x=data(1:319,2)。)。 zi=griddata(x,y,z,xi,yi,39。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 y=data(285:319,3)。 tu=0:10:17133。cubic39。o39。 z=data(285:319,10)。 [xi,yi]=meshgrid(ti,tu)。)。) %5區(qū) Cr分布圖 x=data(285:319,2)。 ti=0:10:20582。 zi=griddata(x,y,z,xi,yi,39。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 y=data(147:284,3)。 tu=0:10:17949。cubic39。o39。 z=data(147:284,10)。 [xi,yi]=meshgrid(ti,tu)。)。) 36 %4區(qū) Cr分布圖 x=data(147:284,2)。 ti=0:10:26424。 zi=griddata(x,y,z,xi,yi,39。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 y=data(81:146,3)。 35 tu=0:10:18397。cubic39。o39。 z=data(81:146,10)。 [xi,yi]=meshgrid(ti,tu)。)。) %3區(qū) Cr分布圖 x=data(81:146,2)。 ti=0:10:28654。 zi=griddata(x,y,z,xi,yi,39。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 y=data(45:80,3)。 tu=0:10:17980。cubic39。o39。 z=data(45:80,10)。 [xi,yi]=meshgrid(ti,tu)。)。) %2區(qū) Cr分布圖 x=data(45:80,2)。 ti=0:10:22674。 zi=griddata(x,y,z,xi,yi,39。 mesh(xi,yi,zi) hold on plot3(x,y,z,39。 y=data(1:44,3)。 tu=0:10:18449。cubic39。o39。 z=data(1:44,10)。 [xi,yi]=meshgrid(ti,tu)。)。) %1區(qū) Cr分布圖 x=data(1:44,2)。 ti=0:10:26453。 zi=griddata(x,y,z,xi,yi,39。 [9]國家環(huán)境保護(hù)局，土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)，北京：中國環(huán)境出版社， 1995。 七 、參考文獻(xiàn) [1] 鄭彤 ,陳春云，環(huán)境系統(tǒng)數(shù)學(xué)模型，北京：化學(xué)工業(yè)出版社， 2022。 a邊界條件 當(dāng) t=0且 0 ZL??時(shí)， C= 0C ； 當(dāng) Z=0且 t0時(shí)， C=0； 當(dāng) Z=L時(shí) ， =0CZ?? 此時(shí)，擴(kuò)散方程的解為： 17 22( 2 1 )014 ( 1 ) ( 2 1 ) ( )( , ) c o s2 1 2n DnnC n L ZC Z t enL? ?? ? ???? ? ?? ?? 揮發(fā)量為： 2201[1 2 ( 1 ) e x p ( ) ]nnDCf n L D tDt? ??? ? ? ?? 式中， L為土壤層的厚度；其余符號(hào)同前。 ⑵ 加入的信息 鑒于以上缺點(diǎn)，我們認(rèn)為，還應(yīng)收集以下信息來完善模型的準(zhǔn)確性， 。 ： 我們使用 matlab軟件的 griddata命令下的 39。22 : , , , 。 ( , ) ( 0 , 0)C x y x y?? 邊界條件為： 0( 0 , 0 , ) 。在四類區(qū)中， Cd、 Pb、 Cu、 Hg、 Zn 的污染貢獻(xiàn)普遍很大，由于四類區(qū)為主干道路區(qū)，其污染原因最為可能是汽車尾氣的排放，以及汽車輪胎的磨損。 等級(jí) P 污染程度 污染水平 1 ≤ 安全 清潔 2 ～ 警戒線 尚清潔 3 ～ 輕污染 超過背景值，視輕污染 4 ～ 中污染 土壤開始受中度污染 5 ≥ 重污染 污染已相當(dāng)嚴(yán)重 表 113 從 上 圖 1 及表 113，可以得出三類區(qū) 的內(nèi)梅 羅污染指數(shù) 剛過警戒線，尚清潔；一類區(qū)和五類區(qū)的內(nèi)梅羅污染指數(shù)介于 ～ ，屬于輕污染；四類區(qū)的內(nèi)梅羅污染指數(shù)介于 ～ ，屬中污染 ； 二 類區(qū)的內(nèi)梅羅污染指數(shù)高于 ，屬于重污染。 ， 利用單因子污染指數(shù)法和內(nèi)梅羅綜合指數(shù)法來評(píng)價(jià)城區(qū)內(nèi)不同區(qū)域金屬的污染程度。 a 各元素在 各區(qū)的情況分布圖，附圖如下： （見附錄二） b 各元素在城區(qū)情況分布圖： （見附錄三） 第二部分： 問題一的解答 4 （一） 問題 1 模型的建立和求解 8 種重金屬元素在該城區(qū)的空間分布，我們基于對(duì)附件一及附件二中海量數(shù)據(jù)的分析，以及查證到的信息，利用 matlab 軟件，繪制出了 8種重金屬元素在該城區(qū)的不同區(qū)域的分布情況圖（ 見附錄一 ），以及 某元素在全城區(qū)的分布情況圖（ 見附錄二 ）。 (2) 通過數(shù)據(jù)分析，說明重金屬污染的主要原因。對(duì)城市土壤地質(zhì)環(huán)境異常的查證，以及 應(yīng)用查證獲得的海量數(shù)據(jù)資料開展城市環(huán)境質(zhì)量評(píng)價(jià)，研究人類活動(dòng)影響下城市地質(zhì)環(huán)境的演變模式，日益成為人們關(guān)注的焦點(diǎn)。 本文 使用 excel 軟件 對(duì)數(shù)據(jù)進(jìn)行了處理 與分析 ， 得 出城區(qū)污染的主要 污染物是 Hg 和 Pb。 我們知道，抄襲別人的成果是違反競賽規(guī)則的 , 如果引用別人的成果或其他公開的資料（包括網(wǎng)上查到的資料），必須按照規(guī)定的參考文獻(xiàn)的表述方式在正文引用 處 和參考文獻(xiàn)中明確列出