反距離加權(quán)平均插值法

1、反距離加權(quán)插值法在污染場(chǎng)地評(píng)價(jià)中的應(yīng)用3陽(yáng)文銳王如松33黃錦樓陳展李鋒中國(guó)科學(xué)院生態(tài)環(huán)境研究 中心城市與區(qū)域生態(tài)國(guó)家重點(diǎn)實(shí)驗(yàn)室北京100085摘 要確定污染場(chǎng)地的土壤修復(fù)范圍在實(shí)際操作中存在很大的難 度.本文以北方某廢棄農(nóng)藥廠為例結(jié)合專家判定采樣法與網(wǎng) 格采樣法采用地統(tǒng)計(jì)學(xué)方法中的反距離加權(quán) 空間插值法分兩種情形將大于土壤環(huán)境基準(zhǔn)值的區(qū)域劃定為治理范圍.結(jié)果表明由于國(guó)家土壤環(huán)境二級(jí)標(biāo)準(zhǔn)主要適用于農(nóng)業(yè)生產(chǎn)的 土地依據(jù)反距離加權(quán) 空間插值后以此標(biāo)準(zhǔn)確定的治理范圍 較大而基于健康風(fēng)險(xiǎn)水平閾值所確定的修復(fù)范圍和修復(fù)成 本較小較合理且經(jīng)濟(jì)可行采用風(fēng)險(xiǎn)評(píng)價(jià)的 反距離加權(quán)插值 法確定污染場(chǎng)地的污染范圍為今后

2、的場(chǎng)地評(píng)價(jià)和土壤修復(fù) 提供了思路.關(guān)鍵詞場(chǎng)地評(píng)價(jià)尼梅羅污染指數(shù)風(fēng)險(xiǎn)評(píng)價(jià)反距離加權(quán)插值土壤修復(fù)文章編號(hào)1001 -9332200709-2013-06 中圖分類號(hào) X171.5X826 文 獻(xiàn)標(biāo)識(shí)碼AApplicationofinversedistanceweightedinterpolationmethodinco ntaminatedsiteassess2ment.YANGWen2ruiWANGRu2songHU ANGJin2louCHENZhanLIFengStateKeyLa2boratoryofUrbanan dRegionalEcologyResearchCenterforEco2

3、EnvironmentalScienc esChineseAcademyofSciencesBeijing100085China.2Chin .J .App I.Ecol.2007189:2013-2018.Abstract:Therearemanydifficultiesin ascertainingtheareathatneedstoberemediedincontami2natedsite.T hispaperintegratedtheexperts ' judgmentswithgridsasthesamplings trategyinana2bandonedpesticide

4、plantinnorthernChinaandapplied ageostatisticalmethodinversedistanceweightedinterpolationIDWt osimulatethepollutants2andriskspatialdistribution.Basedonsoilen vironmentalbenchmarktwoscenariosweredesignedforascertainin gthepollutedarea.Theresultsshowedthattheareaneededtoberemedi edanddeterminedbyIDWwas

5、somewhatlargeraccordingtothenatio nalenvironmental2ndstandardforagriculturalsoilsbuttheareaandc ostneededforremedyweresmallerandmoreeconomicandreasonabl ebasedonthehealthriskthresh2oldlevelandIDWmethod.Acombina tionofIDWandhealthriskassessmentinascertainingpollu2tedarea wouldbeanewwayforsiteassessme

6、ntandsoilremediationinthefutur e.Keywords:siteassessmentNemerowpollutionindexNPIriskasses smentinversedistanceweightedinterpolationsoilremediation.3 國(guó) 家自然科學(xué)基金重點(diǎn)資助項(xiàng)目 70433001.33 通訊作者 .E2mail:2006211230 收稿 2007206216 接 受 .1 引 言當(dāng)前中國(guó)面臨高速的城市化發(fā)展城市有限的 土地空間成為城市發(fā)展的主要限制因子之一.為創(chuàng)建和諧與健康的城市環(huán)境安全合理的土地利用方式成為保障城市居 民健康

7、的關(guān)鍵性要素 .過(guò)去的 50 年由于歷史條件的限制眾多 污染型化工企業(yè)在城市內(nèi)部發(fā)展占用了大量的土地對(duì)當(dāng)?shù)氐淖匀画h(huán)境造成了嚴(yán)重的污染其中一些重金屬、有毒有害化 學(xué)品和持久性有機(jī)污染物POPs隨著土地的重新利用進(jìn)入水體和暴露在大氣中嚴(yán)重威脅人群健康2城市工業(yè)區(qū)土壤環(huán)境污染的問(wèn)題已經(jīng)引起學(xué)者的廣泛關(guān)注11-13.由化工企業(yè)造成的環(huán)境問(wèn)題特別是土地污染問(wèn)題必須得到有效的治理以保 障城市居民的健康北方某大城市為解決這些老化工企業(yè)遺 留下來(lái)的場(chǎng)地污染問(wèn)題已經(jīng)責(zé)令市內(nèi)十幾家大型化工生產(chǎn) 企業(yè)搬出城市外部這些污染型化工企業(yè)遺留下來(lái)的場(chǎng)地污 染問(wèn)題已經(jīng)引起當(dāng)?shù)丨h(huán)保部門(mén)的高度重視.由于污染物在土壤中的遷移不同于

8、空氣和水體介質(zhì)污染物質(zhì)遷移緩慢對(duì)于 大面積的化工企業(yè)土壤污染而言搬遷過(guò)程中對(duì)于場(chǎng)地土壤 的干擾使得土壤中的這些污染物質(zhì)運(yùn)移具有很大的不確定 性給環(huán)境評(píng)價(jià)帶來(lái)了巨大的難度如何在經(jīng)濟(jì)花費(fèi)最小的條 件下使污染的治理達(dá)到土地利用方式的安應(yīng)用生態(tài)學(xué)報(bào) 2007年9月第18卷第9期ChineseJournalofAppliedEcologySep.2007189:2013-2018 全水 平確保污染場(chǎng)地經(jīng)過(guò)評(píng)價(jià)修復(fù)后土地再利用的潛在健康風(fēng) 險(xiǎn)達(dá)到安全水平是當(dāng)前場(chǎng)地評(píng)價(jià)中的一大難題.目前有關(guān)場(chǎng)地評(píng)價(jià)的研究工作大部分停留在以環(huán)境介質(zhì)的環(huán)境質(zhì)量標(biāo) 準(zhǔn)為參考采用逐步采樣的方法評(píng)價(jià)場(chǎng)地是否存在污染其成 本過(guò)高.本文

9、嘗試應(yīng)用 反距離加權(quán) 空間插值法進(jìn)行場(chǎng)地評(píng)價(jià) 提出了基于風(fēng)險(xiǎn)評(píng)價(jià)的土地治理措施為污染場(chǎng)地修復(fù)提供了思路 .2 研究地區(qū)與研究方法 211 自然概況研究場(chǎng)地為 北方某農(nóng)藥廠廠區(qū)面積 4 萬(wàn) m2 經(jīng)過(guò)幾十年的發(fā)展后遺留下 來(lái)的污染物引起周邊居民的不安 20 世紀(jì) 80 年代初該農(nóng)藥廠 廠區(qū)用地轉(zhuǎn)化為某油漆廠生產(chǎn)用地 2005 年全廠拆遷搬出該 地作為未來(lái)房地產(chǎn)開(kāi)發(fā)用途 .場(chǎng)地土壤 pH715915 為粘質(zhì)粉土 有刺激性氣味 .212 樣品采集根據(jù)該廠的歷史生產(chǎn)沿革經(jīng)過(guò) 與原廠區(qū)工作人員以及專家討論從廠區(qū)的平面圖上判別出 可能產(chǎn)生污染物的車(chē)間設(shè)置重點(diǎn)采樣范圍.以整個(gè)廠區(qū)的廠界為分析范圍劃分12mK

10、 12m單元網(wǎng)格以網(wǎng)格交叉點(diǎn)為采樣 點(diǎn)對(duì)重點(diǎn)采樣區(qū)域根據(jù)實(shí)地情況進(jìn)行疏密采樣.采樣工具為工程勘察院的 SH30 型鉆機(jī)每取 1 個(gè)樣本鉆頭均作清理以免 土壤樣本交叉污染 .采樣深度分別為 20、80、180 和 380cm 為客觀真實(shí)地反映場(chǎng)地污染情況在污染嚴(yán)重的點(diǎn)位適當(dāng)進(jìn) 行深度采樣分別增加 500 和 600cm 兩個(gè)土樣 .樣品放置于內(nèi) 襯聚乙烯塑料袋的棉布專用樣品袋 .213 測(cè)定方法由于研究 地塊為農(nóng)藥廠在場(chǎng)地評(píng)價(jià)時(shí)針對(duì)該廠的主要農(nóng)藥產(chǎn)品進(jìn)行 分析本研究主要以六六六 HCH 和滴滴涕 DDT 兩類物質(zhì)作為 討論對(duì)象 .樣品置于陰涼通風(fēng)處自然風(fēng)干剔除樣品中植物根 系、有機(jī)殘?jiān)约翱梢?jiàn)侵

11、入體用木質(zhì)工具碾碎并用瑪瑙研缽 研磨分別過(guò) 10、60 和 100 目尼龍篩 8.依據(jù) GB/T14550 93 采用氣相色譜法測(cè)定 HCH 和 DDT 含量 .數(shù)據(jù)采用 Excel 處理 后采用 Arc2GIS910 軟件進(jìn)行插值分析 .214 評(píng)價(jià)方法 21411 污染評(píng)價(jià) 以國(guó)家土壤環(huán)境質(zhì)量二級(jí)標(biāo)準(zhǔn) GB15618 1995 為標(biāo)準(zhǔn)采用國(guó)內(nèi)外普遍采用的內(nèi)梅羅綜合指數(shù)法 9.PCiSi2 平 均CiSi2最大式中:P為土壤中污染物i的環(huán)境質(zhì)量指數(shù) PW1 為未受污染P1為已污染P值越大污染越嚴(yán)重 Ci為土壤中污 染物 i 的實(shí)測(cè)濃度的統(tǒng)計(jì)平均值 mg? 9? 9kg-1Si 為污染物 i

12、的評(píng)價(jià)標(biāo)準(zhǔn) mg? 9?9kg-1HCH 為 0 1 05mg? 9? 9kg-1 DDT 為 0 1 05mg? 9? 9kg-1.21412 風(fēng)險(xiǎn)評(píng)價(jià) 以 EPA 健康風(fēng)險(xiǎn)評(píng)價(jià)導(dǎo) 則為依據(jù)在對(duì)人體健康進(jìn)行風(fēng)險(xiǎn)評(píng)價(jià)時(shí)主要考慮的攝入途 徑為經(jīng)口攝入 .污染物日攝入量計(jì)算公式4:ADDCX IR >EDX 10-6/BW MAT式中：ADD為污染物日均攝入 量 mg? 9? 9kg-1 ? 9? 9d-1C 為土體中污染物的濃度 mg? 9?9kg-1IR 為日土壤攝入量 mg? 9? 9d-1ED 為人體暴露 的時(shí)間 dBW 為人群平均體質(zhì)量 kg10-6 為轉(zhuǎn)換系數(shù) AT 為人 均壽命

13、d.致癌性污染物按不同年齡段進(jìn)行計(jì)算.致癌性風(fēng)險(xiǎn)度 3:RADD> Q 式中 :R 為致癌性污染物所引起的健康風(fēng)險(xiǎn) a-1ADD 為致癌性污染物的日平均攝入劑量mg? 9?9kg-1 ? 9? 9d-1 按照壽命周期為 70 歲 Q 為污染物的 致癌斜率因子 kg? 9? 9d? 9? 9mg-1 參考 EPA 標(biāo)準(zhǔn)經(jīng)口 HCH 為 118kg?9? 9d-1? 9?9mg-1DDT 為0143kg?9? 9d-1? 9?9mg-1 我國(guó)國(guó)家工業(yè)企業(yè)土壤環(huán)境質(zhì) 量風(fēng)險(xiǎn)評(píng)價(jià)基準(zhǔn) HJ/T25 1 999 規(guī)定的標(biāo)準(zhǔn) :HCH 為 113kg? 9? 9d?9? 9mg-1DDT 為 014

14、3kg? 9? 9d? 9? 9mg-1.風(fēng)險(xiǎn)計(jì)算中按照人體在不同年齡段攝入的累積量計(jì)算表1.表1各人群的土壤日攝入量參數(shù)Tab.1Dailyintakeparameterofsoilfordifferentpeoplegroups 年齡段 Age日攝入量Dailyintakemg ? 9? 9d-1年暴露時(shí)間Yearlyexposuretimed ? 9? 9a-1 平均體質(zhì)量Averagebodymasskg 一般天氣 Ordinaryclimate 兒童Children1630021020 少年 Teenagers71820012045 職業(yè)成年人 Adults196010012060

15、老人 Elders607020015050 揚(yáng)塵天氣 Dustyclimate四類人群Alltypes-100020-參數(shù)參考場(chǎng)地環(huán)境 評(píng)價(jià)指南考慮到北方揚(yáng)塵天氣作了適當(dāng)?shù)男薷腡heco2efficientswererefer redtoGuidestoField ' sEnvironmentalA ssessmentandmodi2fiedforthedustweatherinNorthChina.41O2 應(yīng) 用生態(tài)學(xué)報(bào)18卷圖1不同深度土壤尼梅羅污染指數(shù)空間分布Fig.1SpatialdistributionofNPIindifferentsoildepths.A:20cmB:8

16、0cmC:180cmD:380cmE:500cm.下同 Thesamebelow.關(guān)于風(fēng)險(xiǎn)評(píng)價(jià)的標(biāo)準(zhǔn)美國(guó)環(huán)保局根據(jù)動(dòng)物的毒理學(xué)試驗(yàn)建議致癌化學(xué)物的健康風(fēng)險(xiǎn)可接受水平值為10-6? 9? 9a-110.美國(guó)和歐洲許多國(guó)家對(duì)污染場(chǎng)地進(jìn)行風(fēng)險(xiǎn)評(píng)價(jià)時(shí)亦采用該值對(duì)人體健康進(jìn)行風(fēng)險(xiǎn)評(píng)價(jià).我國(guó)目前在風(fēng)險(xiǎn)評(píng)價(jià)基準(zhǔn)值方面的研究尚顯不足沒(méi)有相關(guān)的標(biāo)準(zhǔn)與參考規(guī)范 本研究以10-6? 9? 9a-1作為 風(fēng)險(xiǎn)評(píng)價(jià)標(biāo)準(zhǔn).21413空間分析 空間插值法被廣泛用于資源管理、災(zāi)害管理和生態(tài)環(huán)境治理中57應(yīng)用較多的有臨近法、反距離權(quán)重法、克里格法、樣條函數(shù)法、趨勢(shì)面法和多元回 歸等一系列模型方法空間異質(zhì)性是空間插值研究的隱含前

17、 提即要素的非均勻空間分布才需要空間插值空間相關(guān)性則 是空間插值研究的基礎(chǔ)16.反距離權(quán)重空間 插值法是一種加權(quán)平均51029期陽(yáng)文銳等：反距離加權(quán)插值法在污染場(chǎng)地評(píng)價(jià)中的應(yīng)用 圖2 不同土壤層六六六和DDT對(duì)兒童健康風(fēng)險(xiǎn)空間分布圖Fig.2 Spatialdistributionofthechildren' shealthriskbyHCHandDDTindifferentsoildepths 10-66102 應(yīng)用生態(tài)學(xué)報(bào)18卷內(nèi)插法該方法認(rèn)為任何一個(gè)觀測(cè)值都對(duì)鄰近的區(qū)域有 影響且影響的大小隨距離的增大而減小10.采用ArcGIS910軟件中的空間分析模塊IDW空間插值工具對(duì)污染評(píng)

18、價(jià)和風(fēng) 險(xiǎn)評(píng)價(jià)的結(jié)果進(jìn)行空間插值.3結(jié)果與分析311污染綜合 指數(shù)空間插值根據(jù) 反距離加權(quán) 法空間插值的結(jié)果對(duì)不同土 壤層綜合污染指數(shù)的 反距離加權(quán) 法空間插值圖1.從圖1可以 看出該廠區(qū)兩種主要污染物在土壤中的分布范圍隨著土壤 深度的增加而呈現(xiàn)逐漸縮小的趨勢(shì).污染主要集中在表層0180cm在380cm 土層污染物濃度急劇下降表明此層以下土 壤污染逐步減輕.根據(jù)該污染指數(shù)空間插值圖及國(guó)家土壤環(huán) 境質(zhì)量二級(jí)標(biāo)準(zhǔn) GB15618 1995對(duì)需要治理的區(qū)域進(jìn)行判斷.結(jié)果表明廠區(qū) 0180cm 的土壤全部受到了不同程度的污染 180500cm 部分土壤需要進(jìn)行治理治理土方量約為 167200m3.31

19、2 風(fēng)險(xiǎn)評(píng)價(jià)空間插值針對(duì)不同土層土壤對(duì)不 同人群的健康風(fēng)險(xiǎn)進(jìn)行空間插值以兒童為例圖 2.從污染物對(duì) 不同人體健康風(fēng)險(xiǎn)概率進(jìn)行空間插值的結(jié)果來(lái)看在所列舉 的 4 類人群中兩種污染物對(duì)兒童的健康風(fēng)險(xiǎn)較大因此選用兒 童風(fēng)險(xiǎn)概率為基準(zhǔn)所制定的土壤修復(fù)標(biāo)準(zhǔn)是比較嚴(yán)格的.按照該兒童風(fēng)險(xiǎn)概率計(jì)算的結(jié)果推導(dǎo)出該廠區(qū)兩種污染物需 要治理的濃度閾值即風(fēng)險(xiǎn)評(píng)價(jià)閾值以該值為基準(zhǔn)確定土壤 治理范圍 .通過(guò)軟件 ArcGIS 對(duì)各圖層進(jìn)行疊加確定該廠區(qū)的 修復(fù)范圍以及各層需要治理的深度圖 3 其需要治理的土方量 約為 130400m3.313 兩種治理方案的比較基于綜合污染指 數(shù)的土壤修復(fù)范圍比基于風(fēng)險(xiǎn)評(píng)價(jià)的土壤修復(fù)范圍

20、要大這 是因?yàn)榍罢叩脑u(píng)價(jià)標(biāo)準(zhǔn)為國(guó)家土壤環(huán)境質(zhì)量二級(jí)標(biāo)準(zhǔn) GB156181995 由于該場(chǎng)地將作為城市居民區(qū)的建設(shè)用地 該標(biāo)準(zhǔn)對(duì)于居民區(qū)而言污染物在土壤中的濃度限制較為嚴(yán) 格 015mg? 9? 9kg-1 使修復(fù)成本較高 .后者通過(guò)對(duì)特定人群進(jìn) 行風(fēng)險(xiǎn)假設(shè)和風(fēng)險(xiǎn)評(píng)價(jià)確定修復(fù)范圍其參數(shù)在設(shè)定時(shí)考慮 到了不同人群因?yàn)椴煌顒?dòng)方式的污染土壤攝入量基于受 害人群的風(fēng)險(xiǎn)閾值推斷土壤污染物治理濃度閾值 HCH 為 2119mg? 9? 9kg-1DDT 為 814mg? 9? 9kg-1. 而國(guó)家工業(yè)企業(yè) 土壤環(huán)境質(zhì)量風(fēng)險(xiǎn)評(píng)價(jià)基準(zhǔn) HJ/T25 1999 的土壤基準(zhǔn)中HCH 為 51mg? 9? 9kg-

21、1DDT 為 195mg? 9? 9kg-1 故其治理范 圍比前者要小也更加經(jīng)濟(jì)合理.而對(duì)污染場(chǎng)地進(jìn)行風(fēng)險(xiǎn)評(píng)價(jià) 就是為了制定污染物的治理濃度閾值6并選擇基于風(fēng)險(xiǎn)概率計(jì)算的土壤治理范圍圖3基于兒童健康風(fēng)險(xiǎn)的土壤治理范 圍 Fig.3Areaneededforremedybasedonchildrenhealthriskassessment.17: 采樣點(diǎn)Samplingsites.10500cm20180cm380500cm30380cm40180cm3 80500cm50380cm60180cm70180cm.4 討 論目前雖然對(duì) 具體污染物的分析監(jiān)測(cè)技術(shù)較成熟對(duì)污染物在土壤中的遷 移和轉(zhuǎn)化機(jī)

22、理的研究較多但對(duì)于城市中受到人為干擾嚴(yán)重 的土壤特別是化工企業(yè)在搬遷過(guò)程中遺留下來(lái)的污染問(wèn)題 在確定治理范圍和治理量上存在很大的不確定性給環(huán)境治 理帶來(lái)很大難度在區(qū)域環(huán)境治理中反距離權(quán)重空間插值法已被應(yīng)用于對(duì)污染物分布范圍的模擬14-15但在局部土壤治理中的應(yīng)用還未見(jiàn)報(bào)道.本研究對(duì) 反距離加權(quán) 空間插值法在 場(chǎng)地土壤污染治理中的應(yīng)用進(jìn)行了嘗試但尚存在以下問(wèn)題需要進(jìn)一步改進(jìn)和完善:1反距離加權(quán)空間插法一種局部插值 法其假設(shè)前提是未知值的點(diǎn)受較近控制點(diǎn)的影響比較遠(yuǎn)控 制點(diǎn)的影響更大1該方法要求對(duì)受影響的局部控制點(diǎn)有清楚 的認(rèn)識(shí)且其結(jié)論直接受到采樣點(diǎn)數(shù)值的影響采樣點(diǎn)越多對(duì) 局部的真實(shí)反映越強(qiáng)利用 反

23、距離加權(quán)插值法 所確定的污染范圍就越準(zhǔn)確.在71029期陽(yáng)文銳等：反距離加權(quán)插值法 在污染場(chǎng)地評(píng)價(jià)中的應(yīng)用 場(chǎng)地評(píng)價(jià)中由于受到采樣和分析成本的限制采樣點(diǎn)數(shù)目有 限其結(jié)論可能不能完全真實(shí)地反映場(chǎng)地的實(shí)際狀況.2采樣點(diǎn)的位置對(duì)空間插值有著重要影響17.采樣點(diǎn)均勻分布能更好地反映研究要素在空間上的分布特征.但在實(shí)際場(chǎng)地評(píng)價(jià)中特別是在拆遷后的場(chǎng)地中可能遇到人為地質(zhì)條件如地下掩 埋的建筑物基底的阻礙導(dǎo)致采樣點(diǎn)位置的挪移從而造成采 樣點(diǎn)空間分布不均勻?qū)Y(jié)果造成一定影響.3采樣網(wǎng)格的大小影響著空間插值的結(jié)果.由于廠區(qū)搬遷或者其他人為活動(dòng)對(duì) 于土壤的干擾污染物質(zhì)分布的不確定性大在局部的場(chǎng)地評(píng) 價(jià)中采取多大的網(wǎng)

24、格進(jìn)行插值分析可以盡可能的減少這些 不確定性需要在今后的研究中進(jìn)一步探索參考文獻(xiàn)1ChangKT.2002.Trans.ChenJ2F 陳健飛.2003.1 ntroductiontoGeographiclnformationSystems.Bei2jing:SciencePress:245-254inChinese2 ChenH2H 陳鴻漢 ChenH2W諶宏偉HeJ2T何江濤 etal.2006.Health2basedriskassessmentofcontaminatedsites:Princ iplesandmethods.EarthScienceFrontiers 地學(xué)前緣 131

25、:216-223inChinese3ChenYCMaHW.2006.Modelcomparisonforriskassessment:Acasestudyofcontaminatedgroundwater.Chemosphere63:751-7614 EPA.1992.GuidelinesforExposureAssessment.Wash2ingtonDC. FederalRegister57104:22888-229385 GuoX2D 郭旭東 FUB2J 傅伯杰 MaK2M 馬克明 etal.2000.Spatialvariabilityofsoilnutrientsbasedongeo

26、statisticsco mbinedwithGIS AcasestudyinZunhuaCityofHebeiProvince.Ch ineseJournalofAppliedEcology 應(yīng)用生態(tài)學(xué)報(bào) 114:557-563inChinese6 JiangL 姜 林 WangY 王 巖 .2004.GuidestoField 'sEnvironmentalAessment.Beijing:Chi naEnvironmentalSciencePress:42-67inChinese7 JiangY 姜 勇 ZhangY2G 張玉革 LiangW2J 梁文舉 etal.2003.Sp

27、atialvariabilityofsoilnutri2entsincultivatedsurfacesoi lofSujiatunDistrictShen2yangCity.ChineseJournalofAppliedEcol ogy 應(yīng)用生態(tài)學(xué)報(bào) 1410:1673-1676inChinese8 LuR2K 魯如 坤.2000.AnalyticalMethodofSoilAg2ro2chemistry.Beijing:Chin eseAgriculturalScienceandTechnologyPress.inChinese9 LuY2S 陸雍森.1999.EnvironmentalA

28、ssessment.Shanghai:TongjiUniversityP ress:531-557inChi2nese10 QuC2X 曲宸緒 JiangY 姜 勇 WuY2P 武燕萍etal.2006.Studyonthedevelopmentofachoro2plethatlasoncancerm ortalityusingtheinversedistanceweightinterpolationinthe1990 hinaJournalofEp2idemiology 中華流行病學(xué)雜志273:230-233inChinese11 WcisoEIovenDKucharskiRetal.200

29、2.Humanhealthriskassessmen tcasestudy:AnabandonedmentalsmeltersiteinPoland.Chemospher e47:507-51512 WuS2T 吳雙桃 WuX2F 吳曉芙 HuY2L 胡曰 利 etal.2004.StudiesonsoilpollutionaroundPb2Znsmeltingfactoryan dheavymetalshyperaccumula2tors.EcologyandEnvironment 生 態(tài)環(huán)境 132:156-157inChinese13 ZhaoQ2N 趙沁娜 XuQ2X 徐 啟新 YangK