王曉亮,趙元藝,柳建平,路璐,楊永強(qiáng),初娜

1) 中國(guó)地質(zhì)大學(xué)地球科學(xué)與資源學(xué)院，北京，100083； 2) 中國(guó)地質(zhì)科學(xué)院礦產(chǎn)資源研究所,北京,100037；3) 北京中色資源環(huán)境工程有限公司,北京,100070

內(nèi)容提要：采用連續(xù)提取法對(duì)德興銅礦大塢河流域的土壤樣品中Cd元素的形態(tài)進(jìn)行分析，并基于緩變型地球化學(xué)災(zāi)害模型，對(duì)大塢河流域土壤中Cd的環(huán)境地球化學(xué)特征進(jìn)行研究。結(jié)果顯示，大塢河上游的祝家村附近和中游的德銅醫(yī)院附近的土壤中Cd元素活動(dòng)性污染總濃度(total concentration of active specie of Cd,以下記為T(mén)CASCd)隨著污染物的可釋放總量(total releseable content of the pollutant of Cd,以下記為T(mén)RCPCd)增加而增加,二者之間呈線性關(guān)系。雖然大塢河上—中游的樣品中TRCPCd范圍均在國(guó)家三級(jí)土壤標(biāo)準(zhǔn)之下，但隨著該地區(qū)土壤中Cd元素的不斷輸入，這種情況下所造成的污染相對(duì)于緩變型地球化學(xué)災(zāi)害造成的污染，應(yīng)當(dāng)更加引起重視。在大塢河下游的下沽口村附近土壤中,Cd的各形態(tài)數(shù)據(jù)特征符合緩變型地球化學(xué)災(zāi)害特征，且有TRCPCd向TCASCd釋放轉(zhuǎn)化的趨勢(shì)。Cd元素的緩變型地球化學(xué)災(zāi)害模型為y=1.12167E-04x3-6.25584E-02x2+1.15298E+01x-6.07452E+02。爆發(fā)臨界點(diǎn)為T(mén)RCPCd =185.90ng/g，有約36%的樣品超過(guò)爆發(fā)臨界點(diǎn)，具有爆發(fā)緩變型地球化學(xué)災(zāi)害的可能性。這種認(rèn)識(shí)對(duì)大塢河流域土壤中Cd元素污染的預(yù)警和對(duì)生態(tài)環(huán)境的改善具有重要意義。

1 德興銅礦環(huán)境地質(zhì)特征

德興銅礦位于江西省上饒地區(qū)德興縣泗洲鎮(zhèn),距德興縣城22 km,是世界上銅資源量在800萬(wàn)噸以上的8個(gè)斑巖型礦床之一,為目前亞洲生產(chǎn)規(guī)模最大的銅礦。德興銅礦包括富家塢、銅廠和朱砂紅3個(gè)礦床,3處共探明可開(kāi)采的礦石總儲(chǔ)量為16億噸,含銅900萬(wàn)噸,含鉬29萬(wàn)噸,還有十余種價(jià)值較高的伴生元素。礦區(qū)是一個(gè)以黃銅礦為主的多元素共生硫化物礦床,共生元素有Cu、S、Fe、Pb、Zn、Mo、Ag、Au、Co、Bi、As、Mn、Si、Al、Ti、Ca、Mg等,礦石儲(chǔ)量大,Cu品位雖低但較均勻,礦石類(lèi)型簡(jiǎn)單,礦體埋藏較淺,在溫?zé)釟夂驐l件下,由于強(qiáng)烈的風(fēng)化和淋溶作用,土體下面形成厚達(dá)數(shù)米至數(shù)十米的紅色風(fēng)化殼,Ca、Mg、K、Na等金屬元素有明顯遷移,而Al、Fe明顯富集,構(gòu)成鐵鋁型風(fēng)化殼。大塢河河流全長(zhǎng)14 km,匯水面積34 km2。流經(jīng)礦區(qū)腹地,在礦區(qū)北部匯入樂(lè)安河,樂(lè)安河最終匯入鄱陽(yáng)湖。大塢河是樂(lè)安河的支流之一,發(fā)源于德興銅礦區(qū)東端官帽山,流貫礦區(qū),流域匯水面積34 km2以上,有明顯的山區(qū)河流特征,河流全長(zhǎng)為14 km,經(jīng)張家坂、沽口村流入樂(lè)安河。大塢河匯集了礦床、巖隙、廢石堆滲出的酸性水,流域內(nèi)環(huán)境污染問(wèn)題已引起研究者的不斷關(guān)注(趙元藝等,2003;初娜等,2006;趙元藝等,2006，2009)。

2 樣品與實(shí)驗(yàn)

2.1 樣品采集

取樣時(shí)間是2004年4月,共采集樣品22件。對(duì)大塢河上、中、下游分別取樣,其中上游6件,中游2件,下游14件。在大塢河上游西岸分別采集了兩個(gè)土壤剖面,下游地區(qū)同樣采集土壤剖面(圖1)。

采樣時(shí)盡量避開(kāi)受居民生活環(huán)境污染的土壤,選擇具有代表性的采樣點(diǎn)。應(yīng)用GPS定位儀確定采樣點(diǎn)的經(jīng)緯度,并記錄所取樣品周?chē)沫h(huán)境狀況和土壤狀況。采集完的樣品，現(xiàn)場(chǎng)需要手工挑出其中的草根、礫石等雜物。然后盛裝在透氣性良好的樣品袋中。

2.2 樣品分析

所采樣品送至中國(guó)地質(zhì)科學(xué)院地球物理地球化學(xué)勘查研究所進(jìn)行分析測(cè)試，用508型原子吸收光度計(jì)對(duì)樣品中Cd元素含量進(jìn)行分析。采用連續(xù)提取法對(duì)樣品中Cd的水溶態(tài)、吸附態(tài)、碳酸鹽態(tài)、有機(jī)結(jié)合態(tài)、硫化物態(tài)、硅酸鹽態(tài)等形態(tài)進(jìn)行分析。

從表1中我們可以得到，Cd元素吸附態(tài)的相對(duì)偏差RE為30.9%，估計(jì)也是由于含量少的原因，但是Cd元素的主要形態(tài)有機(jī)態(tài)和硅酸鹽態(tài)的相對(duì)偏差在5%左右，因此不會(huì)造成大的影響。

從表2中可以看到，Cd大部分樣品的σ(大塢河流域土壤Cd元素總量與形態(tài)含量之和的相對(duì)誤差)在10%～20%范圍內(nèi)，有個(gè)別的樣品σ值偏大，可能由于樣品中Cd的含量很低造成的誤差。

表 1 Cd 土壤元素形態(tài)精密度比較Table 1 The precision of Cd element speciations in soil samples

表 2 大塢河流域土壤Cd元素總量與形態(tài)含量之和的相對(duì)誤差表Table 2 The compariton between Cd element contents and sum of the speciations

注：Σ是Cd元素各形態(tài)含量之和， σ=(Cd元素總量-Σ)/總量·100%。

圖 1 德興銅礦大塢河土壤樣品分布示意圖Fig.1 Sketch map showing the distribution of samples of soil along the Dawu River， Dexing orefield

3 大塢河流域土壤中Cd的元素含量特征

文中Cd元素的水溶態(tài)、吸附態(tài)、碳酸鹽態(tài)、硫化物態(tài)、有機(jī)結(jié)合態(tài)、硅酸鹽態(tài)分別以Cd1、Cd2、Cd3、Cd4、Cd5、Cd6表示。研究區(qū)Cd元素污染物的可釋放總量(total releseable content of the pollutant of Cd,以下記為T(mén)RCPCd)(即Cd元素各形態(tài)含量之和)TRCPCd=Cd1+Cd2+Cd3+Cd4+Cd5+Cd6。由于Cd的硅酸鹽態(tài)性質(zhì)穩(wěn)定，在自然界正常條件下不易釋放，能長(zhǎng)期穩(wěn)定的存在于沉積物中，不易為植物吸收，故在整個(gè)土壤生態(tài)系統(tǒng)中對(duì)食物鏈影響較小(李宇慶等，2004)。所以研究區(qū)Cd元素活動(dòng)性污染總濃度(total concentration of active specie of Cd,以下記為T(mén)CASCd)TCASCd=TRCPCd-Cd6。這六種結(jié)合態(tài)按照其在水中的溶解度的從大到小排列為：Cd1>Cd2>Cd3>Cd4>Cd5>Cd6，其中排在后面的元素形態(tài)常?？梢赞D(zhuǎn)化成靠前的形態(tài)(陳明等，2005)，例如Cd5可以轉(zhuǎn)化成Cd4、Cd3、Cd2、Cd1。

在本次研究中，Cd元素各形態(tài)含量如表3所示，在大塢河上游地區(qū)，Cd元素TRCPCd的總量范圍在122.4～277.8ng/g之間，其平均值為198.9ng/g。Cd元素的有效態(tài)含量TCASCd的范圍在44.8～193.8ng/g之間。其平均值為116.2 ng/g。在大塢河中游地區(qū)，Cd元素TRCPCd的總量范圍在187.4～243.4ng/g之間，其平均值為215.3 ng/g。 Cd元素的有效態(tài)含量TCASCd的范圍在127.4～161.2ng/g 之間，其平均值為144.3 ng/g。在大塢河下游地區(qū)，Cd元素TRCPCd的總量范圍在109.4～293.3ng/g之間，其平均值為180.3 ng/g。Cd元素的有效態(tài)含量TCASCd的范圍在61～222.3ng/g之間，其平均值為98.8 ng/g。

表 3 大塢河流域Cd元素各形態(tài)含量(ng/g)Table 3 All kinds of the speciations of Cd element in Dawu River

注：Cd1為Cd的水溶態(tài)含量; Cd2為Cd的吸附態(tài)含量; Cd3為Cd的碳酸鹽態(tài)含量; Cd4為Cd的硫化物態(tài)含量; Cd5為Cd的有機(jī)物態(tài)含量; Cd6為Cd的硅酸鹽態(tài)含量?！茷镃d元素各形態(tài)含量之和。

4 討論

4.1 大塢河上—中游土壤中Cd元素污染特征

通過(guò)對(duì)大塢河上—中游地區(qū)土壤中Cd元素各形態(tài)含量以及對(duì)大塢河上—中游地區(qū)環(huán)境系統(tǒng)中的污染物可釋放總量(TRCPCd)與環(huán)境系統(tǒng)中的活動(dòng)性污染物總濃度(TCASCd)的關(guān)系圖(圖2)分析可知TRCPCd與TCASCd之間呈線性關(guān)系，而緩變型地球化學(xué)災(zāi)害模型中TRCP存在一個(gè)“緩沖區(qū)”(羅杰等，2012)即TCAS在TRCP“緩沖區(qū)”含量范圍內(nèi)將相對(duì)平穩(wěn)，沒(méi)有明顯的變化。然而，在大塢河上—中游地區(qū)土壤中Cd元素的污染情況卻與之相反，該區(qū)域的土壤中的Cd的活動(dòng)性污染物總濃度(TCASCd)隨著Cd的污染物可釋放總量(TRCPCd)的增加而增加，雖然大塢河上—中游的樣品中TRCPCd范圍在122.4～277.8ng/g之間，均在國(guó)家三級(jí)土壤標(biāo)準(zhǔn)(1000ng/g)之下，但隨著該地區(qū)土壤中Cd元素的不斷輸入，這種情況下所造成的污染相對(duì)于緩變型地球化學(xué)災(zāi)害造成的污染情況，應(yīng)當(dāng)更加引起重視。

4.2 大塢河下游Cd的緩變型地球化學(xué)災(zāi)害的識(shí)別

緩變型地球化學(xué)災(zāi)害是對(duì)“化學(xué)定時(shí)炸彈”(Stigliani,1988; Stigliani et al.,1991; Xie Xuejin,1993)的擴(kuò)展和科學(xué)化。其內(nèi)涵和外延比“化學(xué)定時(shí)炸彈”要豐富得多。其過(guò)程可以用一個(gè)具有普遍意義的非線性方程來(lái)表示。污染物的形態(tài)轉(zhuǎn)化是緩變型地球化學(xué)災(zāi)害發(fā)生的重要機(jī)理，“污染物的可釋放總量”和“活動(dòng)性污染物的總濃度”是個(gè)重要的新概念與“環(huán)境容量”一起構(gòu)成研究緩變型地球化學(xué)災(zāi)害的基本工具(陳明等,2005a;陳明等,2005b)。

典型緩變型地球化學(xué)災(zāi)害的演化過(guò)程是具有多重套合結(jié)構(gòu)特性的非線性過(guò)程,可以劃分為3個(gè)演化階段,每個(gè)階段之間各內(nèi)蘊(yùn)一個(gè)具有特定數(shù)學(xué)特征的臨界點(diǎn)。圖3中,X坐標(biāo)表示環(huán)境系統(tǒng)的污染物可釋放總量(Total releasable content of the pollutant)TRCP(C),Y坐標(biāo)表示環(huán)境系統(tǒng)中的活動(dòng)性污染物總濃度(Total concentration of active specie)TCAS(Q)(陳明等, 2005),兩條虛線分別表示一階和二階導(dǎo)數(shù)的圖形。隨著TRCP的增長(zhǎng),TCAS的增長(zhǎng)趨勢(shì)發(fā)生變化,當(dāng)TRCP的增量為ΔC時(shí),TCAS增長(zhǎng)了ΔQ1,隨著污染物濃度的累積,同樣的ΔC的增長(zhǎng),TCAS增長(zhǎng)了ΔQ2, ΔQ2>>ΔQ1,即TCAS與TRCP的關(guān)系是非線性的,可以用多項(xiàng)式表示如下:

Q=a0+a1C1+a2C2+a3C3+……

圖 2 德興銅礦大塢河上—中游土壤中TRCPCd與TCASCd含量關(guān)系圖Fig.2 The relation between TRCPCd and TCASCd in the soil of upper—middle stream of the Dawu River， Dexing Copper orefield

在一個(gè)演化周期內(nèi),該多項(xiàng)式的最高次數(shù)一般為3式中一階、二階導(dǎo)數(shù)為零處分別代表緩變型地球化學(xué)災(zāi)害爆發(fā)的臨界點(diǎn)、爆發(fā)點(diǎn)。具有特定數(shù)學(xué)特征的臨界點(diǎn)包括:爆發(fā)臨界點(diǎn),Q′=Q″=0時(shí),曲線左側(cè)向下凹右側(cè)向上凹;爆發(fā)點(diǎn),Q′取極大值、Q″=0時(shí),曲線向上凹;積累臨界點(diǎn),Q″取極小值時(shí),曲線左側(cè)向上凹右側(cè)向下凹(陳明等，2005；袁峰等，2008；陶春軍等，2009)。

圖 3 污染物緩變型地球化學(xué)災(zāi)害數(shù)學(xué)模型(陳明等，2005)Fig.3 The mathematical model of delayed geochemical hazard of pollutant(Chen Ming et al., 2005)

通過(guò)對(duì)大塢河下游地區(qū)土壤中Cd元素各形態(tài)含量的分析，以及對(duì)下游地區(qū)環(huán)境系統(tǒng)中的污染物可釋放總量(TRCP)與環(huán)境系統(tǒng)中的活動(dòng)性污染物總濃度(TCAS)的關(guān)系圖的對(duì)比研究發(fā)現(xiàn)：大塢河下游地區(qū)下沽口村附近的土壤中Cd 元素符合緩變型地球化學(xué)災(zāi)害特征，且有TCRPCd向TCASCd釋放轉(zhuǎn)化的趨勢(shì)。Cd元素的緩變型地球化學(xué)災(zāi)害數(shù)學(xué)模型為:

y=1.12167E-04x3-6.25584E-02x2+

1.15298E+01x-6.07452E+02

求其二階導(dǎo)數(shù)y〞并使y〞=0可得TRCPCd=185.90ng/g和TCASCd=185.90ng/g。當(dāng)TRCPCd=185.9ng/g時(shí)曲線左側(cè)向下凹,右側(cè)向上凹,說(shuō)明TRCPCd釋放向TCASCd轉(zhuǎn)化的速度加速,該點(diǎn)為緩變型地球化學(xué)災(zāi)害的爆發(fā)臨界點(diǎn)。當(dāng)TRCPCd的濃度在161.2 ～185.9ng/g之間時(shí),方程的斜率幾乎為零,也就是說(shuō),TRCPCd濃度的增加不會(huì)導(dǎo)致“有效態(tài)的”Cd的增加。當(dāng)TRCPCd的濃度在187.3 ～293.3ng/g之間時(shí),曲線向上凹，在TRCPCd=293.3ng/g時(shí)方程的y′=極大,也就是達(dá)到了緩變型地球化學(xué)災(zāi)害的爆發(fā)點(diǎn),此時(shí)緩變型地球化學(xué)災(zāi)害爆發(fā)達(dá)到最劇烈的階段。統(tǒng)計(jì)顯示，在大塢河下游的下沽口村附近的土壤采樣剖面中約有36%的樣品超過(guò)爆發(fā)臨界點(diǎn)，具有爆發(fā)緩變型地球化學(xué)災(zāi)害的可能性。

圖 4 德興銅礦大塢河下游土壤中TRCPCd與TCASCd關(guān)系圖Fig.4 The relation between TRCPCd and TCASCd in the soil of the lower stream of the Dawu River,Dexing Copper orefield

4.3 大塢河下游Cd的緩變型地球化學(xué)災(zāi)害成因分析

土壤的性質(zhì)易受到母質(zhì)、氣候、地貌和生物等因素影響，其環(huán)境容量的影響因素包括有機(jī)質(zhì)、土壤質(zhì)地、土壤陽(yáng)離子交換量、pH、Eh 等?。其中在 pH 值等于或小于 4.2 的地區(qū)，其森林土壤的酸化以及 Al 的濾除，引起了這些地區(qū)森林大面積的消亡(Hauhs et al.，1986; Schulze et al.，1989)。另外，微生物的富集作用以及重金屬高度污染地區(qū)食肉動(dòng)物的生物活化作用等也會(huì)對(duì)動(dòng)物群落產(chǎn)生巨大的影響(Coughrey et al.，1977；Bengtsson et al.，1984)。

對(duì)于大塢河下游，本文主要考慮的是研究區(qū)酸性降水、礦山酸性排水等對(duì)研究區(qū)土壤的影響。根據(jù)前人研究(吳丹等，2006)，江西省為酸雨污染較嚴(yán)重地區(qū)，又由于低品位礦石經(jīng)人工酸性淋濾萃取銅等多種金屬產(chǎn)生的酸性廢水和重金屬，經(jīng)滲漏和排放的方式排入大塢河中(初娜等，2007，2008)，這就導(dǎo)致了研究區(qū)多酸性土壤。在下游地區(qū)，由于在中游地區(qū)從大山選廠和泗洲選廠排進(jìn)的堿性廢水中含有大量的選礦廠選礦時(shí)用的石灰，故大塢河中碳酸鹽較高，導(dǎo)致下游土壤中碳酸鹽態(tài)含量很高。在下游地區(qū)碳酸鹽態(tài)含量雖然明顯低于硅酸鹽態(tài)，但仍占10%～30%。又根據(jù)碳酸鹽態(tài)的重金屬一般是由于沉淀或共沉淀作用與CaCO3反應(yīng)而沉淀下來(lái)，用弱酸就可將此態(tài)的重金屬溶解出來(lái)(黨志等，2001)。實(shí)驗(yàn)證明，碳酸鹽態(tài)重金屬對(duì)pH變化反應(yīng)強(qiáng)烈，pH的大小決定其形態(tài)的遷移轉(zhuǎn)化(莫爭(zhēng)等，2002；Urasa et al.,1996)。那么隨著研究區(qū)土壤的pH值降低，Cd的碳酸鹽態(tài)逐漸被活化。土壤受到長(zhǎng)期持續(xù)Cd元素輸入積累，最終超過(guò)環(huán)境容量，這將是研究區(qū)爆發(fā)Cd元素的緩變型地球化學(xué)災(zāi)害的主要原因。

4.4 研究區(qū)土壤中Cd元素污染預(yù)警

大塢河小白菜和水稻中Cd元素含量高于背景值,其中Cd元素含量處于國(guó)家蔬菜標(biāo)準(zhǔn)的臨界值(初娜等，2008),因此Cd碳酸鹽態(tài)含量高的地區(qū)，不宜種植農(nóng)作物，因?yàn)檗r(nóng)作物在酸性土壤可能較容易吸收碳酸鹽態(tài)的Cd而造成污染。同時(shí)，累積于土壤表層的Cd由于降水作用，可溶態(tài)部分隨水流動(dòng)很可能發(fā)生遷移，進(jìn)入界面土壤和附近的河流或湖泊，造成二次污染(王云等，1995)。所以對(duì)大塢河流域的土壤中Cd元素的污染進(jìn)行預(yù)警將具有重要意義。

對(duì)于大塢河上—中游地區(qū)，即祝家村和德銅醫(yī)院和附近土壤中Cd元素在上游地區(qū)隨深度變化其碳酸鹽態(tài)含量也呈現(xiàn)一定的規(guī)律性。在樣品4-18T04-1到4-18T04-4中可發(fā)現(xiàn)隨深度的增加，碳酸鹽態(tài)含量所占百分比由37% 到18.4%(表3)。又因?yàn)楸韺油寥乐蠧d元素的遷移轉(zhuǎn)化能力比深層土壤中Cd強(qiáng)，土壤酸性條件增強(qiáng)時(shí)，表層土壤中的Cd將造成更大的危害，以往的研究成果中也曾有同樣的結(jié)論(周啟星，1992)。由于在大塢河上—中游地區(qū)的地層中碳酸鹽化十分廣泛，多在礦床的淺部發(fā)育。在銅廠南山采場(chǎng)和銅廠北山采坑處礦物中CaO和CO2的量占到主導(dǎo)。其中銅廠和朱砂紅兩礦床以白云石化和含鐵白云石化為主(朱訓(xùn)等，1983)。故在上—中游的重金屬碳酸鹽態(tài)可能是重金屬受本身礦床礦物影響所致。根據(jù)表3可知：Cd元素在上游和中游地區(qū)土壤中碳酸鹽態(tài)占主導(dǎo)，在樣品中含量最高達(dá)到103ng/g,所占的百分比在30%左右。一旦該區(qū)域土壤理化性質(zhì)改變，大塢河的上游與中游地區(qū)土壤中的Cd的碳酸鹽態(tài)將會(huì)在弱酸的環(huán)境下轉(zhuǎn)化成Cd的水溶態(tài)和吸附態(tài)，這將導(dǎo)致這一區(qū)域中的Cd將被種植物直接吸收，造成更加嚴(yán)重的污染。

根據(jù)研究區(qū)Cd元素各形態(tài)的含量特征以及大塢河下游土壤中Cd元素污染物可釋放總量與Cd元素活動(dòng)性污染總濃度之間的關(guān)系(圖4)可知：在大塢河下游地區(qū)，即在下沽口一帶，Cd元素具有爆發(fā)緩變型地球化學(xué)災(zāi)害的可能性。也就是說(shuō)大塢河下游地區(qū)的下沽口村一帶，可以判定為緩變型地球化學(xué)災(zāi)害的敏感區(qū)域，應(yīng)加以監(jiān)測(cè)，盡快確定緩變型地球化學(xué)災(zāi)害爆發(fā)時(shí)間。

5 結(jié)論

(1)大塢河上游的祝家村附近和中游的德銅醫(yī)院附近的土壤中Cd元素活動(dòng)性污染總濃度隨著污染物的可釋放總量的增加而增加，二者之間呈線性關(guān)系。雖然大塢河上—中游TRCPCd范圍均在國(guó)家三級(jí)土壤標(biāo)準(zhǔn)之下，但隨著該地區(qū)土壤中Cd元素的不斷輸入，這種情況下所造成的污染相對(duì)于緩變型地球化學(xué)災(zāi)害造成的污染，應(yīng)當(dāng)更加引起重視。

(2)大塢河下游地區(qū)即下沽口村一帶的土壤中，Cd 元素符合緩變型地球化學(xué)災(zāi)害特征。Cd元素的緩變型地球化學(xué)災(zāi)害數(shù)學(xué)模型為y=1.12167E-04x3-6.25584E-02x2+1.15298E+01x-6.07452E+02。當(dāng)TRCPCd=185.9ng/g時(shí)，該點(diǎn)為緩變型地球化學(xué)災(zāi)害的爆發(fā)臨界點(diǎn)。在TRCPCd=293.3ng/g時(shí)，達(dá)到了緩變型地球化學(xué)災(zāi)害的爆發(fā)點(diǎn)。在大塢河下游約有36%的樣品超過(guò)爆發(fā)臨界點(diǎn)，具有爆發(fā)緩變型地球化學(xué)災(zāi)害的可能性。

致謝: 本項(xiàng)研究工作得到陳毓川院士、李家熙研究員、張彥英前院長(zhǎng)的指導(dǎo)。 野外工作期間得到德興銅礦地測(cè)部、環(huán)保部等部門(mén)的支持和孫信芽高級(jí)工程師的指導(dǎo)，相關(guān)樣品處理得到了中國(guó)地質(zhì)科學(xué)院地球物理地球化學(xué)勘查研究所中心實(shí)驗(yàn)室張勤主任的指導(dǎo)，在此一并表示感謝。

注釋/Note

? 楊善謀.2010.銅陵金屬礦集區(qū)土壤中Cu ,Cd元素污染評(píng)價(jià)及其換變型地球化學(xué)災(zāi)害研究. 導(dǎo)師:袁峰. 合肥： 合肥工業(yè)大學(xué)碩士學(xué)位論文,58～59.