土壤重金屬復合污染整體模型研究

李瑩等

摘 要 在重金屬遷移建模研究成果的基礎上，考慮土壤對重金屬的吸附特性，推導出土壤重金屬復合污染物遷移轉化整體模型.為求解該模型，分別引進Langmuir等溫競爭吸附模型與土壤水分遷移模型.最后，以重金屬Zn和Cd為例，數值模擬其在土壤中的垂向遷移過程，并利用MATLAB畫出相關圖形.

關鍵詞 土壤；重金屬；復合污染

中圖分類號 O213.9 文獻標識碼 A

Study on the Heavy Metal

Combined Pollution Model in Soil

LI Ying1，WANG Zhizhong1，XIAO Guoguang2

（1.College of Mathematics and Statistics， Central South University， Changsha，Hunan 410083，China；

2. Changsha Research Institute of Mining and Metallurgy， Changsha，Hunan 410012， China）

Abstract Based on the previous achievements and considering the soils adsorption characteristics of heavy metal， this paper proposed an overall migration model of heavy metal compound contaminants in soil. In order to solve the problem， the Langmuir isotherm competitive adsorption model and water migration model in soil were applied. Finally， an numerical simulation example of Zn and Cd was shown as heavy metals vertical migration in soil and displayed by MATLAB.

Key words soil； heavy metal； combined pollution

1 引 言

土壤是人類賴以生存的主要自然資源之一，也是生態系統的重要組成部分[1].近年來，隨著工業的發展，特別是采礦與冶煉的迅速發展，使得重金屬對土壤的污染越來越嚴重[2-4].土壤重金屬污染是指人為因素有意或無意的將對人類本身和其他生命有害的物質施加到土壤中，使其某種成分的含量明顯高于原有含量，并引起土壤環境質量惡化的現象[5].因其能沿食物鏈富集、放大，較低濃度的重金屬便可能損害自然生態系統和人體健康，所以如何有效地減少和修復重金屬污染，緩解其對人體、生物和生態環境的危害，是當今環境領域所面臨的難題之一[6，7].現實環境中的重金屬污染多為多元素共存所造成的復合污染， 而且各種污染元素之間還存在著某種形式的聯合作用.因此， 研究如何較好的表征重金屬的復合污染的綜合效應， 實現對重金屬復合污染的預測、評價和防治，具有重大的現實意義[8].由于土壤及地下水重金屬污染具有隱蔽性以及地下重金屬運移具有很多的不確定性，并且人們對重金屬在地下水中的彌散機理和各種化學反應過程都了解不夠，目前對重金屬在土壤中遷移轉化的整體模擬不多，大多數是對這一系統的分散的零碎的研究，并沒有將這些過程聯系起來綜合研究，而且其多考慮單一重金屬污染，對重金屬復合污染數學模型的研究很少.

本文綜合前人所建立的模型， 將重金屬在土壤中的遷移轉化作為一個整體進行研究，并加入復合污染因素，推導出土壤重金屬復合污染整體模型.對模型進行數值模擬時，為了模型準確性 ，對Langmuir模型和土壤水分遷移模型進行改進，用以說明土壤重金屬復合污染整體模型中參數之間的關系，而不是簡單的以常數代替參數.最后以重金屬Zn和Cd為例，顯示模型的應用步驟.

2 重金屬復合污染遷移轉化模型

土壤由泥沙顆粒、孔隙水和孔隙氣組成，即重金屬在土壤中由固相、液相和氣相組成.地表表層土體是非飽和土，地表深層土體是含地下水的土，為飽和土.從地表向下，土壤一般分為兩個帶，上面是非飽和帶，下面是飽和帶.

與重金屬在水體中的遷移轉化模型類似[9，10]，目前更適合建立重金屬復合污染遷移轉化的整體模型，而不是分相模型（分別就固相，液相，氣相顆粒建立模型）.基于這一結論和土壤中遷移建模的成果，簡化的結果和概念在本文中將推廣和延伸，應用于推導重金屬在土壤中遷移轉化三維方程中[11].

在土壤中取一個控制體dx，dy，dz，假定它由泥沙顆粒空間、孔隙水空間和空隙氣體空間三部分組成.土壤中的重力水通過孔隙向下流失，假定其空隙大小沒有改變，即孔隙度水空間是常數.忽略生物化學作用傳質，不考慮生物體對重金屬的富集，忽略孔隙中的泥沙運動.對于泥沙顆粒濃度的變化，吸附解吸是生化反應的主要機理.

為了推導方便，首先定義下列參數：

6 結 論

本文首先推導了土壤中重金屬復合污染整體模型，為了對該模型進行數值模擬，分別引進Langmuir等溫競爭吸附模型與土壤水分遷移模型，得出固相濃度Si與液相濃度Ci的關系以及液相滲流速度q與含水率θ的數值表達式，并利用Matlab畫出含水率隨時間和垂向距離變化圖.最后運用所得模型對重金屬Zn，Cd進行數值模擬，利用Matlab分別畫出Zn和Cd的液相濃度隨垂向距離和時間t的變化圖.從圖看出：深度增加，含水率減少，土壤中Zn，Cd濃度減少；Zn，Cd共存時吸附量明顯小于單一離子存在時的吸附量；土壤的吸附使重金屬的遷移明顯滯后；土壤中Zn的競爭吸附能力大于Cd的.endprint

參考文獻

[1] 林青.土壤中重金屬Cu、Zn、Pb吸附及遷移的實驗研究[D].青島：青島大學環境科學系，2008.

[2] 陳麗蓉，陳小羅，羅志輝，易曉明，江國防.土壤中重金屬污染物鉻的遷移轉化及治理[J].科技創新導報，2011，3（33）：122-123.

[3] S CHOTPANTARAT， S K ONG， C SUTTHIRAT，et al. Competitive modeling of sorption and transport of Pb2+，Ni2+，Mn2+ and Zn2+ under binary and multimetal system in lateritic soil columns[J]. Geoderma， 2012，190（1-3）：278-287.

[4] Uta KRAUS， Jens WIEGAND. Longterm effects of the Aznalc llar mine spill—heavy metal content and mobility in soils and sediments of the Guadiamar river valley（SW Spain） [J]. Science of the Total Environment， 2006，367（2-3）：855-871.

[5] 周東美，王慎強，陳懷滿.土壤中有機污染物的交互作用[J].土壤與環境，2000，9（2）：143-145.

[6] 韋朝陽，陳同斌.重金屬超富集植物及植物修復技術研究進展[J].生態學報，2001，21（7）：196-1203.

[7] 李戰，李坤.重金屬污染的危害與修復[J].現代農業科技，2010，7（16）：268-270.

[8] 孟昭福，薛澄澤，張增強，等.土壤中重金屬復合污染的表征[J].AgroEnvironmental Protection，1999，17（2）：87-91.

[9] 黃歲樑.沖積河流重金屬污染物遷移轉化數值模擬控制方程及其物理意義[J] 中國科學：技術科學， 2010，40（5） ： 515-524.

[10]黃歲樑.泥沙運動在重金屬污染物遷移轉化中的作用[D].北京：中國水利水電科學研究院，1993.

[11]紀書華.多孔介質中重金屬反應性運移的數值模擬研究[D].青島：青島大學化學與環境工程學院，2009.

[12]張宏，張敬華.生物吸附的熱力學平衡模型和動力學模型綜述[J].天中學刊，2009，24（5）：19-22.

[13]楊亞提，張一平.銅在土壤—溶液界面吸附—解吸特性的研究[J].西北農業學報，1998，7（4）：82-85.

[14]楊德軍.基于SPAC系統的土壤水動力學模型研究[D].杭州：浙江大學土木系，2009.

[15]張圣微，雷玉平，鄭力，等.基于GIS的土壤水分運動模型[J].土壤通報，37（6）：1006-1070.

[16]趙彥旭，張虎元，呂擎峰，等.壓實黃土非飽和滲透系數試驗研究[J].巖土力學，2010，31（6）：1809-1812.

[17]王偉，傅臣家，李光永，等.超水質對土壤水分擴散率的影響[J].山東農業大學學報：自然科學版，2011，42（4）：627-631.

[18]張耀峰.特征有限元及其在土壤水分運動數值模擬中的應用[D]. 西安：西安理工大學理學院，2009.endprint

參考文獻

[1] 林青.土壤中重金屬Cu、Zn、Pb吸附及遷移的實驗研究[D].青島：青島大學環境科學系，2008.

[2] 陳麗蓉，陳小羅，羅志輝，易曉明，江國防.土壤中重金屬污染物鉻的遷移轉化及治理[J].科技創新導報，2011，3（33）：122-123.

[3] S CHOTPANTARAT， S K ONG， C SUTTHIRAT，et al. Competitive modeling of sorption and transport of Pb2+，Ni2+，Mn2+ and Zn2+ under binary and multimetal system in lateritic soil columns[J]. Geoderma， 2012，190（1-3）：278-287.

[4] Uta KRAUS， Jens WIEGAND. Longterm effects of the Aznalc llar mine spill—heavy metal content and mobility in soils and sediments of the Guadiamar river valley（SW Spain） [J]. Science of the Total Environment， 2006，367（2-3）：855-871.

[5] 周東美，王慎強，陳懷滿.土壤中有機污染物的交互作用[J].土壤與環境，2000，9（2）：143-145.

[6] 韋朝陽，陳同斌.重金屬超富集植物及植物修復技術研究進展[J].生態學報，2001，21（7）：196-1203.

[7] 李戰，李坤.重金屬污染的危害與修復[J].現代農業科技，2010，7（16）：268-270.

[8] 孟昭福，薛澄澤，張增強，等.土壤中重金屬復合污染的表征[J].AgroEnvironmental Protection，1999，17（2）：87-91.

[9] 黃歲樑.沖積河流重金屬污染物遷移轉化數值模擬控制方程及其物理意義[J] 中國科學：技術科學， 2010，40（5） ： 515-524.

[10]黃歲樑.泥沙運動在重金屬污染物遷移轉化中的作用[D].北京：中國水利水電科學研究院，1993.

[11]紀書華.多孔介質中重金屬反應性運移的數值模擬研究[D].青島：青島大學化學與環境工程學院，2009.

[12]張宏，張敬華.生物吸附的熱力學平衡模型和動力學模型綜述[J].天中學刊，2009，24（5）：19-22.

[13]楊亞提，張一平.銅在土壤—溶液界面吸附—解吸特性的研究[J].西北農業學報，1998，7（4）：82-85.

[14]楊德軍.基于SPAC系統的土壤水動力學模型研究[D].杭州：浙江大學土木系，2009.

[15]張圣微，雷玉平，鄭力，等.基于GIS的土壤水分運動模型[J].土壤通報，37（6）：1006-1070.

[16]趙彥旭，張虎元，呂擎峰，等.壓實黃土非飽和滲透系數試驗研究[J].巖土力學，2010，31（6）：1809-1812.

[17]王偉，傅臣家，李光永，等.超水質對土壤水分擴散率的影響[J].山東農業大學學報：自然科學版，2011，42（4）：627-631.

[18]張耀峰.特征有限元及其在土壤水分運動數值模擬中的應用[D]. 西安：西安理工大學理學院，2009.endprint

參考文獻

[1] 林青.土壤中重金屬Cu、Zn、Pb吸附及遷移的實驗研究[D].青島：青島大學環境科學系，2008.

[2] 陳麗蓉，陳小羅，羅志輝，易曉明，江國防.土壤中重金屬污染物鉻的遷移轉化及治理[J].科技創新導報，2011，3（33）：122-123.

[3] S CHOTPANTARAT， S K ONG， C SUTTHIRAT，et al. Competitive modeling of sorption and transport of Pb2+，Ni2+，Mn2+ and Zn2+ under binary and multimetal system in lateritic soil columns[J]. Geoderma， 2012，190（1-3）：278-287.

[4] Uta KRAUS， Jens WIEGAND. Longterm effects of the Aznalc llar mine spill—heavy metal content and mobility in soils and sediments of the Guadiamar river valley（SW Spain） [J]. Science of the Total Environment， 2006，367（2-3）：855-871.

[5] 周東美，王慎強，陳懷滿.土壤中有機污染物的交互作用[J].土壤與環境，2000，9（2）：143-145.

[6] 韋朝陽，陳同斌.重金屬超富集植物及植物修復技術研究進展[J].生態學報，2001，21（7）：196-1203.

[7] 李戰，李坤.重金屬污染的危害與修復[J].現代農業科技，2010，7（16）：268-270.

[8] 孟昭福，薛澄澤，張增強，等.土壤中重金屬復合污染的表征[J].AgroEnvironmental Protection，1999，17（2）：87-91.

[9] 黃歲樑.沖積河流重金屬污染物遷移轉化數值模擬控制方程及其物理意義[J] 中國科學：技術科學， 2010，40（5） ： 515-524.

[10]黃歲樑.泥沙運動在重金屬污染物遷移轉化中的作用[D].北京：中國水利水電科學研究院，1993.

[11]紀書華.多孔介質中重金屬反應性運移的數值模擬研究[D].青島：青島大學化學與環境工程學院，2009.

[12]張宏，張敬華.生物吸附的熱力學平衡模型和動力學模型綜述[J].天中學刊，2009，24（5）：19-22.

[13]楊亞提，張一平.銅在土壤—溶液界面吸附—解吸特性的研究[J].西北農業學報，1998，7（4）：82-85.

[14]楊德軍.基于SPAC系統的土壤水動力學模型研究[D].杭州：浙江大學土木系，2009.

[15]張圣微，雷玉平，鄭力，等.基于GIS的土壤水分運動模型[J].土壤通報，37（6）：1006-1070.

[16]趙彥旭，張虎元，呂擎峰，等.壓實黃土非飽和滲透系數試驗研究[J].巖土力學，2010，31（6）：1809-1812.

[17]王偉，傅臣家，李光永，等.超水質對土壤水分擴散率的影響[J].山東農業大學學報：自然科學版，2011，42（4）：627-631.

[18]張耀峰.特征有限元及其在土壤水分運動數值模擬中的應用[D]. 西安：西安理工大學理學院，2009.endprint