【正文】 途徑（只分析傳播特征，未具體考慮其傳播途徑）； （ 5）考慮到海拔高度與濃度的關(guān)系比較復(fù)雜，根據(jù)數(shù)據(jù)分析，相同水平距離的兩點(diǎn)之間的海拔差距與濃度的關(guān)系具有不確定性，所以 簡(jiǎn)化了模型； （ 6）由于地域差異等各種因素的影響，沒有統(tǒng)一的評(píng)價(jià)標(biāo)準(zhǔn)，沒有對(duì)未知點(diǎn)進(jìn)行估計(jì)模擬后再計(jì)算。 土壤重金屬含量屬于空間連續(xù)的變量，具有區(qū)域化變量特征。半方差函數(shù)是地統(tǒng)計(jì)學(xué)中的主要工具，一方面是利用半方差函數(shù)對(duì)參數(shù)的空間分布進(jìn)行結(jié)構(gòu)分析和變異性分析，另一方面是應(yīng)用結(jié)構(gòu)分析的結(jié)果和克里格法進(jìn)行估值?？梢跃C合考慮各種水環(huán)境因素對(duì)其遷移過程的影響，水質(zhì)影響最直接的溶解態(tài)重金屬為研究對(duì)象，不考慮重金屬的懸浮態(tài)和底泥態(tài)，借助 BODDO 水質(zhì)模型中有機(jī)降解的概念，以綜合沉降系數(shù) k 來分析懸浮態(tài)和底泥態(tài)重金屬的沉降和再懸浮關(guān)系（沉降為正，懸浮為負(fù)）。s 為污染源。土壤環(huán)境質(zhì)量綜合評(píng)價(jià)方法探討 [J].山東科技大學(xué)學(xué)報(bào)（自然科學(xué)版）， 2020,24（ 4）： 113118. [8] 張朝生，章中，何建邦 .長(zhǎng)江水系沉積物重金屬含量空間分布特征研究 ——地統(tǒng)計(jì)學(xué)方法 [J].地理學(xué)報(bào)， 1997,52(2):184192. 。 考慮到不同地區(qū)的地質(zhì)因素不同，因此我們?cè)谖廴疚餄舛忍荻葦?shù)學(xué)模型的基礎(chǔ)上進(jìn)行改進(jìn)，設(shè)置一個(gè)地質(zhì)因子 γ，不同地質(zhì)重金屬擴(kuò)散的速度不同，所以 γ的取值也不同。2xc c cE kc st x x?? ? ?? ? ? ?? ? ? 19 其中 :c 為溶解態(tài)重金屬濃度， t 為時(shí)間， u 為縱向水流流速， Ex 為縱向彌散系數(shù)， h。 為更好地研究城市地質(zhì)環(huán)境演變模式，還需要該城區(qū)的水文分布情況的信息，與其相關(guān)的不同時(shí)間段的采樣數(shù)據(jù)，該城區(qū)地質(zhì)情況及深層土壤采樣數(shù)據(jù)。并運(yùn)用 GIS 技術(shù)對(duì)重金屬進(jìn)行空間分布的分析。這種空間異質(zhì)性包括土壤的水分特征、物理性質(zhì)、化學(xué)性質(zhì)、土壤重金屬及其他元素的性質(zhì)隨著空間位置不同而發(fā)生的變異，正是因?yàn)檫@種高度復(fù)雜的空間異質(zhì)性，使得對(duì)土壤變化的空間動(dòng)態(tài)的研究變得十分困難。 檢驗(yàn)結(jié)果：求得 Pb 、 Cu 和 Cd 元素污染源的坐標(biāo)與圖對(duì)比的結(jié)果相符。從圖上可以看出這 些 k 值較大的點(diǎn)對(duì)應(yīng)的采樣點(diǎn)序號(hào)集中在 五 個(gè)區(qū)域附近。 檢驗(yàn)結(jié)果：求得 Cr 和 Ni元素污染源的坐標(biāo)與圖對(duì)比的結(jié)果相符。 將這 50個(gè)匯總結(jié)果數(shù)據(jù) 用 excel 直觀表示（圖 5）如下： 圖 5 重金屬元素 Cr 和 Ni的采樣點(diǎn)篩選結(jié)果 圖中橫坐標(biāo)表示采樣點(diǎn)序號(hào)；縱坐標(biāo)表示采樣點(diǎn)數(shù)目。 結(jié)果檢驗(yàn)： 首先用附件給采樣點(diǎn)的坐標(biāo)參數(shù)及對(duì)應(yīng)的 Zn 濃度值，以濃度代替高度，用matlab 繪制出位置與濃度的類等高線關(guān)系圖（圖 6）： 0 50 100 150 200161116212631364146采樣點(diǎn)序號(hào) Zn 15 圖中橫坐標(biāo)、縱坐標(biāo)分別代表采樣點(diǎn)的 x 與 y 值，不同顏色區(qū)域代表不同濃度（其中紅色代表 Zn 元素的濃度最高）。 將這 50 個(gè)匯總結(jié)果數(shù)據(jù) 用 excel 直觀表示（圖 9）如下： 14 圖 9 重金屬元素 Zn 的采樣點(diǎn)篩選結(jié)果 圖中橫坐標(biāo)表示采樣點(diǎn)序號(hào)；縱坐標(biāo)表示采樣點(diǎn)數(shù)目。 結(jié)果檢驗(yàn)： 首先用附件給采樣點(diǎn)的坐標(biāo)參數(shù)及對(duì)應(yīng)的 Hg 濃度值，以濃度代替高度，用matlab 繪制出位置與濃度的類等高線關(guān)系圖（圖 6）： 圖中橫坐標(biāo)、縱坐標(biāo)分別代表采樣點(diǎn)的 x 與 y 值，不同顏色區(qū)域代表不同濃度（其中紅色代表 Hg 元素的濃度最高）。 將這 50 個(gè)匯總結(jié)果數(shù)據(jù) 用 excel直觀表示（圖 7）如下： 圖 7 重金屬元素 Hg 的采樣點(diǎn)篩選結(jié)果 圖中橫坐標(biāo)表示采樣點(diǎn)序號(hào)；縱坐標(biāo)表示采樣點(diǎn)數(shù)目。 結(jié)果檢驗(yàn)： 首先用附件給采樣點(diǎn)的坐標(biāo)參數(shù)及對(duì)應(yīng)的 As 濃度值，以濃度代替高度，用matlab 繪制出位置與濃度的類等高線關(guān)系圖（圖 6）： 0 50 100 150 200 2501713192531374349采樣點(diǎn)序號(hào) As 12 圖中橫坐標(biāo)、縱坐標(biāo)分別代表采樣點(diǎn)的 x 與 y 值，不同顏色區(qū)域代表不同濃度（其中紅色代表 As 元素的濃度最高）。 得到集中區(qū)域在序號(hào) 27， 84， 178 這三個(gè)采樣點(diǎn)附近。因此，在污染源附近的采樣點(diǎn)之間的污染物濃度梯度值 k 較大（即單位距離上的污染物濃度變化值較大） 。 對(duì)與 Pb 、 Cu 、 Cd 三 種元素，由于 Pb 與其它三種元素的相關(guān)性都較好，這表明 Pb 的濃度變化與其它三種元素變化都比較接近，所以它的含量變化可以倚重其它三種元素。所以對(duì)該城區(qū)的 8 種重金屬元素含量之間進(jìn)行相關(guān)性分析可以知道其來源是否相同。 農(nóng)業(yè)上化肥、農(nóng)藥農(nóng)用化學(xué)品的施用是土壤重金屬污染的主要途徑之一 ， 重金屬污染來自化肥、農(nóng)藥農(nóng)用化學(xué)品，如磷肥、各類復(fù)合肥及不少殺蟲劑、除草劑中，都不同程度含有重金屬（ Cd、 Hg、 Pb、 As 等） 。 汽車輪胎磨損以及煤 燃燒產(chǎn)生的粉塵、煙塵中均含有鋅及化合物 與 農(nóng)業(yè)上用灌溉 用水 含鋅污染物 ， 這也導(dǎo)致了 主干道區(qū)和生活區(qū)鋅的含量 較高。 該城區(qū)內(nèi)工業(yè)區(qū)、主干道區(qū)表層土壤中 Cu 的含量相對(duì)較高 .由于工農(nóng)業(yè)生產(chǎn)的快速發(fā)展， Cu 的用量不斷增加，含 Cu 的污染物排放越來越多 [3].主干道區(qū)機(jī)動(dòng)車輛含 Cu 零件的磨損也是該表層土壤中 Cu 含量增多的主要因素之一 . 工業(yè)區(qū)和主干道區(qū) Hg 的含量明顯高于其它區(qū)域。 Pb 在土壤中有積累效應(yīng)，使得工業(yè)區(qū) Pb的含量高于其他區(qū)域 。 鉻化合物 是 電鍍、染料、制藥、皮革、顏料等制造企業(yè) 的主要原料，這些企業(yè) 排放 的 廢棄物及生活垃圾 中的這些產(chǎn)品都會(huì)導(dǎo)致 Cr 污染 。 ? :公園綠地區(qū)； ? ：主干道路區(qū)； ：山區(qū)； ：生活區(qū)； ? ：工業(yè)區(qū) 圖 2 功能區(qū)等高線分布圖 圖 2 是圖 1 在平面上投影得到的等高線圖 。對(duì)于問題三，這 是本文章的難點(diǎn)，也是亮點(diǎn)。 (4) 分析你所建立模型的優(yōu)缺點(diǎn)，為更好地 研究城市地質(zhì)環(huán)境的演變模式 ，還應(yīng)收集什么信息？有了這些信息，如何建立模型解決問題？ 的 假設(shè) 與符號(hào)說明 模型的假設(shè) ，且無重大地殼變動(dòng)致使土壤地質(zhì)環(huán)境改變； ； （ 0~10 厘米深度）進(jìn)行取樣的深度 相同 ； ； ，未受人類活動(dòng)的任何影響； 符合說明 P 表示單 因子評(píng)價(jià)指數(shù) iP 表示 土壤中污染物 i 的指數(shù) ip 表示 土壤中 污染物 i 的實(shí)測(cè)含量 p 表示 綜合評(píng)價(jià)指數(shù) 3 maxip 表示 最嚴(yán)重污染物的單因子指數(shù) r 表示相關(guān)性檢驗(yàn)系數(shù) 1C 表示 Cr 和 Ni 的綜合指標(biāo) 2C 表示 Cu 、 Pb 和 Cd 的綜合指標(biāo) h 表示某個(gè)采樣點(diǎn)的海拔高度 k 表示任意兩個(gè)采樣點(diǎn)的重金屬濃度梯度 ,ijcc 分別表示采樣點(diǎn) i 和采樣點(diǎn) j 的重金屬濃度 ,xy 分別表示某個(gè)采樣點(diǎn)在某一坐標(biāo)系下的橫坐標(biāo)與縱坐標(biāo) 問題一要求繪制出 8 種主要重金屬元素在該城區(qū)的空間分布圖，考慮到里面含有的因素（坐標(biāo)，海拔 h ， 功能區(qū) ， 金屬濃度）過多，繪制圖形比較困難所以可以將這工作分為幾部分：第一部分，繪制地形圖；第二部分，繪制 x， y， h 與功能區(qū)的圖；第三部分，繪制各功能區(qū)不同金屬的濃度圖。 附件 1 列出了采 樣點(diǎn)的 位置、海拔高度及其所屬功能區(qū)等信息，附件 2 列出了 8 種主要重金屬 元素 在采樣點(diǎn)處 的濃度，附件 3 列 出了 8 種主要重金屬元素的背景值。 按照功能劃分，城區(qū)一般可分為 生活區(qū)、 工業(yè)區(qū)、 山區(qū)、 主干道路 區(qū)及公園綠地區(qū)等， 分別記 為 1 類區(qū)、 2 類區(qū)、??、 5 類 區(qū)， 不同的區(qū)