江西贛州梓山地區(qū)富硒土壤重金屬元素安全性評價

郄海滿,文幫勇,王繼強,周強強,劉冰權(quán)

(江西省地質(zhì)調(diào)查研究院,南昌 330030)

江西贛州梓山地區(qū)富硒土壤重金屬元素安全性評價

郄海滿,文幫勇,王繼強,周強強,劉冰權(quán)

(江西省地質(zhì)調(diào)查研究院,南昌 330030)

在江西梓山地區(qū)采集土壤樣712件、水稻樣46件,分析該區(qū)土壤、水稻中As、Cd等重金屬元素含量及其在土壤—水稻系統(tǒng)中的安全性。研究區(qū)土壤環(huán)境質(zhì)量總體優(yōu)良,清潔和較清潔土壤圖斑面積為11 883.06 hm2,占全區(qū)面積的96.82%;輕微污染的土壤圖斑面積為368.66 hm2,占全區(qū)面積的3.00%;重金屬污染的土壤圖斑主要分布于富硒土壤區(qū)。水稻中As等重金屬元素含量較低,均符合國家食品安全標(biāo)準(zhǔn)。水稻對As等元素的吸收富集能力受其土壤總量及其有效性制約,而重金屬元素的生物有效性主要受重金屬總量與TFe2O3的控制。研究區(qū)土壤脫硅富鐵鋁的特征表明,該區(qū)土壤As等重金屬元素的生物有效性較低,可開發(fā)富硒水稻。

重金屬元素;富硒土壤區(qū);安全性評價;遷移轉(zhuǎn)化

江西省地質(zhì)調(diào)查研究院經(jīng)過10余年的多目標(biāo)區(qū)域地球化學(xué)調(diào)查及專題研究,在江西省內(nèi)發(fā)現(xiàn)大面積的富硒土壤資源,并將之成果應(yīng)用化。2014年,位于豐城市面積約50 km2的“中國生態(tài)硒谷”產(chǎn)值達25億元,取得了顯著的經(jīng)濟效益。但富硒土壤往往也富集重金屬,江西省已發(fā)現(xiàn)的富硒土壤中Cd元素超標(biāo)率達25.59%。重金屬污染問題一直是研究熱點,目前已開展了富硒土壤區(qū)重金屬元素含量、形態(tài)特征及其生態(tài)效應(yīng),Se與重金屬元素的交互關(guān)系等研究[1-10],但是關(guān)于重金屬超標(biāo)富硒土壤資源研究及評價鮮有報道。本文對江西省贛州梓山地區(qū)富硒土壤中的As、Cd等重金屬元素進行安全性評價,并研究As、Cd等重金屬元素在土壤—水稻系統(tǒng)中的遷移轉(zhuǎn)化規(guī)律,為江西省煤系地層區(qū)富硒土壤資源開發(fā)提供參考。

1 研究區(qū)概況

研究區(qū)位于江西省贛州市于都縣東北部約16 km處,東西寬約16.5 km,南北長約19 km,面積約122 km2。地理坐標(biāo)為: 115°28′16″～115°38′09″ E,25°54′42″～26°04′00″ N。研究區(qū)交通便利,國道G323、贛龍鐵路與省道S218橫穿該區(qū)。研究區(qū)屬亞熱帶季風(fēng)氣候,雨量充沛,四季分明,氣候溫和,年平均氣溫為19.7 ℃,年平均降雨量達1 507 mm。

江西省地質(zhì)調(diào)查研究院于2011～2015年開展“江西興國—尋烏地區(qū)1∶250 000多目標(biāo)區(qū)域地球化學(xué)調(diào)查”時首次發(fā)現(xiàn)該區(qū)表層土壤中存在較大面積的富硒土壤。2013年綜合考查贛州市富硒土壤空間分布、土地利用情況及交通樞紐位置,在于都縣東部梓山地區(qū)進行1∶5 0000土壤環(huán)境地質(zhì)調(diào)查評價。該區(qū)以低丘和崗地為主,廣布紅壤及水稻土,耕地主要種植水稻等糧食作物,園地主要種植臍橙、油茶等經(jīng)濟作物。該區(qū)出露震旦系、泥盆系、石炭系、二疊系、侏羅系、白堊系和第四系地層。該區(qū)中部和北東部石炭系、二疊系地層較發(fā)育,并有煤層出露。富硒土壤主要分布于研究區(qū)中部和北東部地區(qū)。

2 樣品采集及分析測試

2.1 樣品采集

共采集46件水稻樣品,712件土壤樣品(包括水稻根系土樣品46件)(圖1),其中水田、旱地、蔬菜地樣點垂直向下均勻采集0～20 cm土壤,種植果樹、油茶等的園地及適宜發(fā)展為園地的林地樣點垂直向下均勻采集0～60 cm土壤;水田樣品采集時間為作物收割后的農(nóng)閑期,采集的各樣點土壤經(jīng)手工掰碎,挑出根系、秸稈、石塊以及蟲體等雜物,樣品采樣時避開溝渠、林帶、路邊、舊房基、糞堆及微地形高低不平、無代表性的地段。土壤、植物樣品采集、切割、洗滌和保存等環(huán)節(jié)嚴(yán)格按照《區(qū)域生態(tài)地球化學(xué)評價技術(shù)要求》執(zhí)行。

圖1 研究區(qū)樣點位置圖Fig. 1 Map show sampling locations in the study area

2.2 樣品采集

采集的土壤元素全量分析在國土資源部武漢礦產(chǎn)資源監(jiān)督檢測中心完成,水稻樣品元素分析及土壤樣品形態(tài)分析在國土資源部合肥礦產(chǎn)資源監(jiān)督檢測中心完成。土壤樣品主要分析Se、As、Cd、Cr、Hg、Pb、Cu、Zn、Ni等元素,植物樣品主要分析Se、As、Cd、Cr、Hg、Pb等元素，不同樣品的分析方法及檢出限見表1。所有測試結(jié)果均滿足DD2005-3《生態(tài)地球化學(xué)評價樣品分析技術(shù)要求》規(guī)定的精密度和準(zhǔn)確度要求,對重復(fù)樣相對偏差(RD=|A1-A2|/(A1+A2)×100)的統(tǒng)計分析表明,樣品的合格率均>85%,檢測結(jié)果可信。

表1 樣品分析方法及檢出限/(mg·kg-1)

3 土壤中As等重金屬元素含量特征及環(huán)境質(zhì)量評價

3.1 土壤中As等重金屬元素含量特征

對富硒土壤資源進行開發(fā),除土壤富硒外,生產(chǎn)基地的土壤環(huán)境質(zhì)量還應(yīng)達到國家土壤環(huán)境質(zhì)量GB15618-1995二級及以上標(biāo)準(zhǔn)。研究區(qū)土壤中As等重金屬含量元素測試結(jié)果如表2所示,其平均含量由高到低依次為Zn>Cr>Pb>Cu>Ni>As>Cd>Hg。

3.2 環(huán)境質(zhì)量評價

參照國家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)GB15618-1995二級標(biāo)準(zhǔn)限量值,采用單項污染指數(shù)法對研究區(qū)土壤進行單指標(biāo)土壤環(huán)境地球化學(xué)等級劃分評價,結(jié)果見表2。研究區(qū)土壤中As、Cd等8種重金屬元素的單項環(huán)境質(zhì)量等級多為一級,而受污染的三級、四級及五級分布圖斑面積較少,如土壤As一級(清潔區(qū))土壤圖斑面積為12 004.35 hm2,占全區(qū)面積的97.8%;二級(較清潔區(qū))土壤圖斑面積為22.17 hm2,占全區(qū)面積的0.18%;三級(輕度污染區(qū))土壤圖斑面積為244.52 hm2,占全區(qū)面積的1.99%;四級(中度污染區(qū))土壤圖斑面積為1.93 hm2,占全區(qū)面積的0.02%;五級(重度污染區(qū))土壤在該區(qū)幾乎沒有分布。

在單指標(biāo)土壤環(huán)境地球化學(xué)等級劃分的基礎(chǔ)上,采用一票否決法進行土壤環(huán)境地球化學(xué)綜合等

表2 土壤中As等元重金屬素含量及環(huán)境質(zhì)量等級統(tǒng)計

注：含量單位為mg/kg,面積單位為hm2,比例單位為%，n=712。

級評價,即每個評價單元的土壤環(huán)境地球化學(xué)綜合等級等同于單指標(biāo)劃分出環(huán)境等級最差的等級。結(jié)果表明,研究區(qū)土壤環(huán)境質(zhì)量總體優(yōu)良,一級清潔區(qū)(未受任何重金屬污染)的土壤圖斑面積為11 318.55 hm2,占研究區(qū)總面積的92.22%;二級較清潔區(qū)的土壤圖斑面積為564.51 hm2,占研究區(qū)總面積的4.60%;三級輕度污染區(qū)的土壤圖斑面積為368.66 hm2,占研究區(qū)總面積的3.00%;中度污染區(qū)和重度污染區(qū)的土壤圖斑面積分別占研究區(qū)總面積的0.11%和0.06%。

研究區(qū)土壤Se含量為0.08～4.18 mg/kg,平均含量為0.29 mg/kg,富硒樣品共154件,富硒率為21.63%,富硒土壤區(qū)主要分布于梅江與貢水之間的低丘及崗地,而河流兩岸的農(nóng)田富硒土壤分布面積較少。對研究區(qū)富硒土壤及非富硒土壤樣品中As等重金屬元素含量進行統(tǒng)計(表3),可知,富硒土壤區(qū)As等8種重金屬元素含量均高于非富硒土壤區(qū)的重金屬元素含量,尤其富硒土壤區(qū)的As、Cr、Ni等元素含量分別是非富硒土壤區(qū)As、Cr、Ni等元素含量的2.48倍、2.22倍和2.25倍。參照國家土壤環(huán)境質(zhì)量標(biāo)準(zhǔn)GB15618-1995二級標(biāo)準(zhǔn)限量值,與非富硒土壤相比,富硒土壤的As、Cd、Cr、Cu、Hg、Ni等6種重金屬元素超標(biāo)的樣品數(shù)和超標(biāo)率較高,如富硒土壤樣品中Cd元素超標(biāo)20件,超標(biāo)率為12.99%;而非富硒土壤樣品中Cd元素超標(biāo)12件,超標(biāo)率為2.15%。這可能由于煤系地層中除富含Se元素外,也富含As、Cd等重金屬元素。巖石自然風(fēng)化及人類礦產(chǎn)開采活動為梓山地區(qū)土壤提供豐富的Se外,也帶來了Cd、Ni、Cr等重金屬元素。

表3 研究區(qū)土壤富硒區(qū)及非富硒區(qū)重金屬元素含量統(tǒng)計/(mg·kg-1)

4 水稻中As等重金屬元素安全性評價及富集特征

水稻中As等重金屬元素含量見表4。由表4可知,稻米中Cd、Hg、Pb、Cr等重金屬元素含量達標(biāo),其中1件稻米樣品As元素含量為0.568 mg/kg,略微超標(biāo)(0.5 mg/kg)。參照國家食品衛(wèi)生限量值標(biāo)準(zhǔn)(表4),稻米中As等重金屬元素含量均達標(biāo),食品合格率為100%。

表4 梓山地區(qū)稻米As等元素含量/(mg·kg-1)

注：*數(shù)據(jù)來源NY5115-2008;**數(shù)據(jù)來源NY861-2004；n=46。

稻米中Se元素平均含量為0.087 mg/kg,含量范圍為0.039～0.509 mg/kg。參照國家富硒稻谷標(biāo)準(zhǔn)硒含量區(qū)間0.04～0.3 mg/kg[11],該地區(qū)有45件稻谷樣品達到富硒標(biāo)準(zhǔn),富硒率為97.83%。

假定水稻中重金屬元素均來自土壤,統(tǒng)計梓山地區(qū)稻米的生物富集系數(shù)(表4),結(jié)果表明稻米對不同重金屬元素的富集能力差異較大,其中稻米對Cd和Se的富集能力較強,富集系數(shù)分別為42.65%和30.68%,稻米對As和Hg的富集系數(shù)分別為5.16%和3.56%,而稻米對Pb、Cr的富集能力較弱,均<1%。

研究區(qū)46件水稻根系土中,富硒樣品為5件,足硒樣品為36件，乏硒樣品為5件,說明在不同Se含量(富硒、足硒、乏硒)的農(nóng)田中均可產(chǎn)富硒水稻。水稻根系土中有2件樣品Cd元素含量超標(biāo),1件Hg元素含量超標(biāo),即在不同重金屬元素含量(清潔、超標(biāo))的農(nóng)田中均可產(chǎn)符合國家食品衛(wèi)生標(biāo)準(zhǔn)的安全水稻,說明水稻對Se和As等重金屬元素的富集能力不僅與土壤中重金屬元素總量有關(guān),還與其有效含量或生物有效性相關(guān)[12-15]。

研究區(qū)水稻Se和As、Cd等重金屬元素的富集系數(shù)與其根系土相應(yīng)元素含量呈負(fù)相關(guān)(圖2),隨著水稻根系土As、Cd等重金屬元素含量上升,水稻對As、Cd等重金屬元素的生物富集能力下降,如根系土Cd元素含量由0.20 mg/kg上升至0.30 mg/kg時,水稻對Cd元素的富集系數(shù)以乘冪指數(shù)方程形式由約50.00%快速下降至10.00%左右。

圖2 研究區(qū)水稻As等元素富集系數(shù)與其根系土元素含量關(guān)系散點圖Fig. 2 Scatter plots showing enrichment coefficients of elements such as As in rice and contents of elements in root soils in the study area

5 土壤As等重金屬元素生物有效性的影響因素

土壤重金屬元素的有效量反映了植物可吸收利用的含量水平,其有效性(元素有效量與其總量比值)是影響植物對重金屬元素富集能力的重要因素。重金屬元素的離子交換態(tài)、水溶態(tài)、碳酸鹽態(tài)可直接被植物吸收,對植物危害大,而重金屬其他形態(tài)很難直接被植物吸收,對植物的危害較小。本文將植物可直接吸收利用的水溶態(tài)、離子交換態(tài)和碳酸鹽態(tài)三者之和稱為活動態(tài)(可利用態(tài)),將腐植酸態(tài)、鐵錳氧化態(tài)、強有機結(jié)合態(tài)及殘渣態(tài)四者之和稱為穩(wěn)定態(tài)[9-10]。

研究區(qū)表層土壤As等重金屬元素形態(tài)百分含量統(tǒng)計見表5。由表5可知,土壤中As等重金屬元素各形態(tài)含量變化較大,但多數(shù)元素(除Cd)穩(wěn)定態(tài)百分含量占優(yōu)勢,范圍為82.39%～97.42%。

表層土壤可被植物直接吸收的活動態(tài)Cd百分含量最高,為43.29%,其次為活動態(tài)Pb,百分含量為17.61%?；顒討B(tài)Ni、活動態(tài)Cu、活動態(tài)Zn、活動態(tài)As、活動態(tài)Se和活動態(tài)Cr的百分含量分別為10.37%、9.47%、5.04%、4.88%、2.74%和2.58%。

表5 表層土壤As等重金屬元素形態(tài)百分含量統(tǒng)計/%(n=13)

將As等重金屬元素活動態(tài)百分含量與總量、有機質(zhì)、pH、TFe2O3作相關(guān)性矩陣分析(表6)。可知,土壤As、Cr、Cu、Zn、Ni和Se活動態(tài)百分含量與其總量呈較好的負(fù)相關(guān),其相關(guān)系數(shù)分別為-0.77、-0.71、-0.72、-0.77、-0.56和-0.67,隨土壤元素總量增加,其可被植物吸收利用的活動態(tài)百分含量降低(圖3)。對比圖2和圖3可知,水稻As等重金屬元素的富集系數(shù)、土壤As等重金屬元素的生物有效性均與土壤中相應(yīng)元素具有類似的變化趨勢,說明土壤中As等重金屬元素的生物有效性影響水稻對As等重金屬元素的吸收。

表6 土壤中As等重金屬元素活動態(tài)百分含量與其總量、有機質(zhì)、pH及TFe2O3等的相關(guān)矩陣(n=13)

活動態(tài)Cd、Zn與土壤有機質(zhì)呈較好的正相關(guān),相關(guān)系數(shù)分別為0.66和0.61,其他元素與土壤有機質(zhì)沒有顯著的相關(guān)性。

活動態(tài)Cu、Pb與土壤pH呈較好的負(fù)相關(guān),相關(guān)系數(shù)分別為-0.51和-0.56,而其他元素與土壤pH沒有顯著的相關(guān)性。研究區(qū)土壤呈酸性,如形態(tài)樣品土壤pH為4.62～5.76,變化幅度較小且酸化較嚴(yán)重,可能導(dǎo)致多數(shù)土壤中的重金屬元素形態(tài)含量無顯著變化。

土壤As、Cd、Cr、Zn和Se活動態(tài)百分含量與土壤TFe2O3呈較好的負(fù)相關(guān),其相關(guān)系數(shù)分別為-0.75、-0.50、-0.56、-0.59和-0.51,可能由于梓山地區(qū)自然形成的土壤具有脫硅富鐵鋁的特征,而鐵鋁氧化物形成的Fe(OH)3膠體或粘土礦物等將吸附或固化土壤中的As等元素,降低了這些元素的遷移活動能力和生物可利用性。

圖3 研究區(qū)土壤中As等重金屬元素活動態(tài)百分含量與其總量散點圖Fig. 3 Scatter plots showing the relationship between percentage compositions of heavy elements and the total amount in soil in the study area

總體而言,土壤中As等重金屬元素的生物有效性主要受其總量及TFe2O3含量的制約,隨總量、TFe2O3含量的增高,重金屬元素更傾向于以穩(wěn)定態(tài)形式賦存,活動態(tài)百分含量降低,重金屬的生物有效性也隨之降低,水稻對重金屬元素的富集能力隨之下降。土壤中雖存在重金屬超標(biāo)現(xiàn)象,但其生物有效性較低,不利于重金屬元素向水稻中聚集,水稻中的重金屬元素含量較低。水稻高硒低重金屬元素的特征說明,研究區(qū)存在重金屬輕微超標(biāo)的農(nóng)田可以開發(fā)富硒大米。

6 結(jié)論及建議

(1)研究區(qū)土壤環(huán)境質(zhì)量總體優(yōu)良,清潔和較清潔的土壤圖斑面積為11 883.06 hm2,占全區(qū)面積的96.82%。研究區(qū)土壤局部存在As、Cd、Cr、Cu、Hg、Ni等重金屬元素輕微超標(biāo),且主要分布于富硒土壤區(qū)。

(2)46件水稻樣品中As、Cd、Cr、Pb、Hg等元素含量較低,符合國家食品安全水平。水稻中Se含量達0.087 mg/kg,富硒率高達97.83%。水稻富硒貧As等重金屬元素是開發(fā)水稻富硒農(nóng)產(chǎn)品的前提。

(3)水稻對As等重金屬元素的吸收富集能力受土壤中相應(yīng)的As等重金屬元素總量及其有效性制約,重金屬元素的生物有效性主要受其總量與TFe2O3的控制。梓山地區(qū)土壤As等重金屬元素的生物有效性較低,不利于水稻對As等重金屬元素的吸收。在重金屬元素含量輕微超標(biāo)的高背景區(qū),可生產(chǎn)出符合國家食品安全的稻谷,也為煤系地層區(qū)的富硒土壤資源開發(fā)提供參考。

(4)因不同農(nóng)作物對重金屬元素的敏感性不同,本次研究僅采集水稻1種介質(zhì),下一步需采集其他農(nóng)作物樣品進一步探討富硒農(nóng)產(chǎn)品的安全性，同時需增加水稻根、莖、葉中重金屬元素含量測試,為土壤—重金屬元素的遷移提供生物學(xué)依據(jù)。

[1] 戎秋濤,翁煥新. 環(huán)境地球化學(xué)[M]. 北京:地質(zhì)出版社,1990:231-233.

[2] 文幫勇,張濤亮,李西周,等.江西龍南地區(qū)富硒土壤資源開發(fā)可行性研究[J].中國地質(zhì),2014,41(4):256-263.

[3] 宋明義,劉建新,黃春雷.浙北富硒土壤地球化學(xué)特征與生物學(xué)效應(yīng)[J].廣東微量元素科學(xué),2012, 19(3):61-68.

[4] 宋明義,馮雪外,周濤發(fā),等.浙江典型富硒區(qū)硒與重金屬的形態(tài)分析[J].現(xiàn)代地質(zhì),2008,22(6):960-965.

[5] 牟明輝,石楊程,張曉晴,等.恩施富硒茶園土壤重金屬和氟含量及風(fēng)險評價[J].河南農(nóng)業(yè)科學(xué),2016,45(5):571-579.

[6] 李慧,魏昌華,鮑征宇,等.恩施富硒茶葉中Se含量與對應(yīng)土壤中Se及重金屬元素As、Cd、U的關(guān)系[J].地質(zhì)科技情報,2011,30(3):103-107.

[7] 耿建梅,王文斌,溫翠萍,等.海南稻田土壤硒與重金屬的含量、分布及其安全性[J].生態(tài)學(xué)報,2012,32(11):3477-3486.

[8] 肖鈺杰.富硒土壤重金屬在農(nóng)作物中的遷移轉(zhuǎn)化特征及硒的影響效應(yīng)[D].海口:海南大學(xué),2014:43-54.

[9] 楊明桂,王光輝,徐梅桂,等.江西省及鄰區(qū)濱太平洋構(gòu)造活動的基本特征[J].華東地質(zhì),2016,37(1):10-18.

[10] 賈十軍.安徽省富硒土壤評價標(biāo)準(zhǔn)及富硒土壤成因淺析[J].資源調(diào)查與環(huán)境,2013,34(2):133-138.

[11] 中華人民共和國國家質(zhì)量監(jiān)督檢驗檢疫總局,中國國家標(biāo)準(zhǔn)化管理委員會. GB/T富硒稻谷[M].北京：中國標(biāo)準(zhǔn)出版社,2009.

[12] 雷鳴,廖柏寒,秦普豐,等. 土壤重金屬化學(xué)形態(tài)的生物可利用性評價[J]. 生態(tài)環(huán)境學(xué)報,2007,16(5):1151-1556.

[13] 雷鳴,廖柏寒,秦普豐,等. 礦區(qū)污染土壤Pb、Cd、Cu和Zn的形態(tài)分布及其生物活性的研究[J]. 生態(tài)環(huán)境學(xué)報,2007,16(3):807-811.

[14] 文幫勇,黃錦,張濤亮,等. 贛南大余地區(qū)農(nóng)田土壤As、Cd元素化學(xué)形態(tài)含量分析[J].華東地質(zhì),2015,36(4):298-306.

[15] 陳學(xué)誠,李生志,董文庚,等.江西省大余污灌田中鎘的形態(tài)分布及治理研究[J]. 河北輕化工學(xué)院學(xué)報,1989(Z):67-72.

Safetyevaluationofheavymetalcontentsinselenium-richsoilintheZishanarea,Ganzhou,JiangxiProvince

QIE Hai-man, WEN Bang-yong,WANG Ji-qiang,ZHOU Qiang-qiang,LIU Bing-quan

(GeologicalSurveyofJiangxiProvince,Nanchang330030,China)

712 soil samples and 46 rice samples were collected in the Zishan area of Jiangxi Province to test the contents of heavy metals (such as As, Cd, etc) in the soil and rice and finally evaluate safety of land environment in the Zishan area. The test results show that soil environment quality of the study area is generally good. The clean soil area covers an area of 11 883.06 hm2,which takes up 96.82% of the total area. The slightly polluted soil area has an area of 368.66 hm2,which takes up 3.00% of the total area. The soil polluted by heavy metals is mainly distributed in Se-rich soil areas. The contents of heavy metals such as arsenic in rice samples are low and meet the national food safety standards. Absorption and enrichment ability of rice to heavy metal elements such as As is controlled by soil amount and validity. While bio-availability of heavy metals is mainly controlled by the total heave metal amount and TFe2O3. Soil characteristic of desiliconization and ferrallitic enrichment indicates that the biological effectiveness of heavy metals such as As is low in the area, which can be used to grow Se-rich rice.

heavy metals; Se-rich soil area; safety evaluation; migration and transformation

X131.1

：A

：2096-1871(2017)03-234-07

10.16788/j.hddz.32-1865/P.2017.03.010

2016-07-18

：2016-11-23責(zé)任編輯：譚桂麗

中國地質(zhì)調(diào)查局“江西興國-尋烏地區(qū)1∶250 000多目標(biāo)區(qū)域地球化學(xué)調(diào)查(編號：J【2011】01-03-07)”資助。

郄海滿，1984年生，女，工程師，主要從事環(huán)境地球化學(xué)及農(nóng)業(yè)地質(zhì)研究。

文幫勇，1984年生，男，工程師，主要從事環(huán)境地球化學(xué)及農(nóng)業(yè)地質(zhì)研究。