【正文】 模型進行評價。評價每個因素對地質(zhì)環(huán)境質(zhì)量的影響，以原始資料為前提，結(jié)合前人研究資料 并參照實際情況予以綜合判定。地質(zhì)災害易發(fā)程度反映滑坡、崩塌、泥石流、地面塌陷的易發(fā)程度 。層次分析法 [4]是一種定性和定量相結(jié)合的評價方法。具體計算步驟如下 : 1)建立層次結(jié)構(gòu)模型 。若滿足一致性要求，即可對底層可操作要素歸一化得到要素的組合權(quán)重，否則重新組建兩兩比較判斷矩陣。其他數(shù)值類似。通過計算礙出的各指標權(quán)重如表 3。 ix ： 評價因子的性狀數(shù)據(jù)取值 。 建議 與展望 污染物來源的鑒別，在微觀上對人們認識污染元素在土壤中的化學行為及其與作物和環(huán)境的關系等方面的機理提供重要的證據(jù)，宏觀上是對環(huán)境質(zhì)量現(xiàn)狀、污染程度進行正確評價和對污染源進行準確、有效治理的前提。又如在進行剖面分析時，采集深層樣品將會比較耗費與耗時等等。此外，與人類活 21 動有關的微量元素可以通過很多途徑進入土壤，它們在土壤中的行為和趨向依據(jù)它們的來源和形態(tài)有所不同。 缺點： 1）沒有考慮城市大氣干濕沉降和城市土壤重金屬長期累積的監(jiān)測，以及隨季節(jié)變化和年際變化的研究分析； 2）由于忽略了所給數(shù)據(jù)中涉及到重金屬污染物的濃度隨海拔高度而產(chǎn)生的變化，又由于計算機模擬帶有一定的隨機性，以致得到模型的三個指標不是很讓人滿意。y=A(:,2)。v439。 % request output from contour clabel(cs) % add labels identifying heightsxlabel(39。) figure,meshc(X,Y,Z)%三維曲面下面帶有等值線 figure,mesh(X,Y,Z)%三維曲面 30 附錄 2： 城區(qū)表層土壤 8 種重金屬元素含量空間分布圖： As Cd 31 Cr Cu 32 Hg Ni 33 Pb Zn 。),ylabel(39。%插值 pcolor(X,Y,Z)。 scatter(x,y,5,z)%散點圖 figure [X,Y,Z]=griddata(x,y,z,linspace(min(x),max(x),200)39。 22 參考文獻 [1] 劉紹貴 .張?zhí)伊?.王興祥 .潘劍軍 , 南昌市城郊表層土壤重金屬污染特征研究[J].土壤通報， 2020,4（ 2） [2] 敬久旺 .趙玉紅 .張涪平 .李成芳 .曹湊貴 .劉宗磊， 藏中礦區(qū)表層土壤重金屬污染評價 [J].貴州農(nóng)業(yè)科學， 2020,39（ 7）： 126~128 [3] 王政權(quán)，地統(tǒng)計學在生態(tài)學中的應用 [M],北京：科學出版社， 1999 [4] 田福林， 受體模型應用于典型 環(huán)境介質(zhì)中多環(huán)芳烴、二嚷英和多氯聯(lián)苯的來源解析研究 [J],大連理工大學博士論文， 2020 [5] 王昌全，成都平原城市化土壤重（類）金屬演變及環(huán)境效應研究 [D]，重慶：西南農(nóng)業(yè)大學， 2020 [6] 李紹生， 永城礦區(qū)土壤重金屬污染評價研究 [J]，河南科學， 2020,29（ 4） [7] 劉衍君 .湯慶新 .白振華等， 基于地質(zhì)累積與內(nèi)梅羅指數(shù)的耕地重金屬污染研究 [J].中國農(nóng)學通報， 2020,25（ 20）： 174178. [8] 楊忠平，長春市城市重金屬污染的生態(tài)地球化學特征及其來源解析 [J]，吉林大學博士學位 論文， 2020 附件 附錄 1： 表層土壤中 Cr 元素含量的空間分布的 Matlab 程序： A=[ 74 781 1373 731 1321 1791 0 1787 1049 2127 1647 2728 2883 3617 2383 3692 2708 2295 2933 1767 4233 895 4043 1895 2427 3971 3526 4357 5062 4339 4777 4897 5868 4904 6534 5641 5481 6004 4592 4603 2486 5999 3299 6018 3573 6213 23 4741 6434 5375 8643 5635 7965 5394 8631 5291 7349 4742 7293 4948 7293 5567 6782 7004 6226 7304 5230 7048 4600 8180 4496 9328 4311 9090 5365 8049 5439 8077 6401 8017 7210 6869 7286 7056 8348 7747 8260 8457 8991 9460 8311 9062 7639 9319 6799 10631 6472 10685 5528 10643 4472 11702 4480 11730 5532 11482 6354 10700 8184 10630 8774 11678 8618 11902 7709 13244 7056 12746 8450 12855 8945 13797 9621 14325 8666 15467 8658 12442 4329 13093 4339 13920 5354 14844 5519 24 16569 6055 16387 6609 16061 7352 15658 7594 14298 7418 14177 6684 15092 6936 12778 5799 17044 10691 17087 11933 17075 12924 17962 12823 18413 11721 19007 11488 18738 10921 17814 10707 18134 10046 17198 9810 17144 9081 18393 9183 19767 8810 21006 8819 21091 9482 22846 9149 23664 9790 22304 10527 21418 10721 21439 11383 20554 11228 20201 10774 21072 10404 20215 9951 18993 12371 19968 12961 21766 12348 22674 12173 22535 11293 25221 5795 26453 5577 26416 6508 27816 5581 25361 6423 24065 7353 25998 7032 25 27177 7771 26424 8639 26073 8807 24631 9422 24702 9522 25461 9834 24813 10799 26086 11094 26015 12078 27700 11609 27696 11621 27346 13331 26591 13715 27823 14737 27232 14482 24580 13319 24153 12450 22965 13535 23198 13523 24685 14278 28654 8755 24003 15286 21684 13101 22193 12185 17079 5894 15255 5110 15007 5535 3518 2571 3469 2308 3762 2170 3927 2110 4153 2299 3267 793 4684 1364 5495 1205 5664 1653 5541 2093 5451 2757 4020 2990 4026 3913 5101 4080 5438 3994 5382 3012 5314 2060 26 5503 1127 5636 133 6605 374 7093 1381 7100 2449 6837 3490 7906 3978 8045 3052 8394 2035 8403 1075 8079 0 9663 1288 9469 2286 9178 3299 9095 3975 10225 3821 10210 2789 10340 1764 11557 1581 11415 2585 11649 3515 12734 4015 12696 3024 12400 2060 12591 1063 13765 1353 13694 2357 13855 3345 14862 2524 14896 1603 15387 729 15810 2307 16032 3061 15801 3966 15087 3512 16872 2798 17734 3629 16823 4207 17008 4775 17203 6218 17005 7212 16947 7487 16301 8299 17904 8287 27 18303 7385 18438 6539 18556 5588 18954 4874 18012 4414 19072 8519 20282 8590 21475 8540 21450 7555 20261 7586 19569 7348 19411 6934 19501 6091 20582 6548 19909 5300 21018 5764 22176 5492 23359 5325 23238 6502 22624 4818 21703 6591 5006 8846 5734 9659 6395 10443 7405 10981 8446 11200 7612 11938 7912 12840 8866 13143 9296 13102 9475 12020 9212 11305 8629 12086 7776 10613 8622 10638 9237 9872 8307 9726 7106 9467 6423 8831 7458 8920 8904 8868 10547 9591 10398 10360 10395 11203 28 11529 11243 11563 10298 11646 9381 12641 9560 14000 8970 14207 9980 14065 10987 12734 10344 12727 7691 14173 11941 15467 12080 15140 11101 15198 10100 15248