王美华

(中化地质矿山总局 浙江地质勘查院，浙江 杭州 310002)

0 引 言

硒(Se)是人体必不可少的微量元素之一[1]，直接参与机体代谢，具有抗癌、抗辐射、增强人体免疫力以及抗氧化的作用，同时还对镉、汞、砷等有害重金属产生很强拮抗作用，有效阻挡重金属对人体健康造成的潜在危害[2-3]，被称为“生命之火”。大量研究表明，人体缺硒会引起一系列的疾病，常导致身体免疫力下降，引发关节炎、克山病等疾病，严重的会引起该区癌症病例的集中出现，如四川盐亭、江苏扬中等地为癌症高发病区[4-7]。中国土壤中的硒含量总体并不高，有一条大致从东北向西南直至青藏高原方向典型的缺硒带，其中30%为严重缺硒地区[8]。人体摄入食物中硒的主要途径来自农产品，而农产品植物中的硒主要来自土壤中的可溶性硒酸盐、亚硒酸盐等，通过植物根部吸收，进入植物。因此，富硒土壤是天然富硒农产品生产的物质基础，利用和开发天然的富硒土壤已越来越受到各地政府的重视。前人研究发现控制和影响土壤硒含量的因素很多，其中成土母质与富硒土壤关系密切[9-10]，寒武系荷塘组黑色岩系是浙西富硒土壤的主要来源[11]，高硒地质背景对土壤硒有明显的控制作用，其他土壤类型、土地利用方式等亦较大地影响土壤硒迁移，但每个地区影响程度又存在差异[12]。

通过常山县土地质量地球化学调查，在白石镇、何家乡、紫港街道中部、招贤镇东南部圈定了富硒土壤面积10 685.34 hm2[13]。据刘道荣等[12]相关研究，该县的土壤硒富集主要来源于自然背景，寒武系荷塘组黑色岩系是富硒土壤的主要来源。浙西石煤资源丰富，随着国家对自然生态环境的日趋重视，大量石煤矿山相继停采，并历史遗留大量废弃石煤矿山。前人对富硒土壤主要以区域面上调查研究为主[12-17]，针对典型石煤矿山周边耕地的相关研究较少。为系统研究石煤矿山赋矿岩层寒武系荷塘组黑色岩系与富硒土壤的相互关系，对浙江省常山县辉埠地区典型矿山周边耕地进行表层土壤、农产品等样品采集和相关野外调查，研究了石煤矿山周边耕地的土壤硒含量与分布特征，并对主要影响因素进行了讨论，以期为富硒土地资源开发利用、种植结构调整与农田土壤污染修复等提供技术支持。

1 材料和方法

1.1 研究区概况

图1 研究区地质简图Fig.1 Geological sketch map of the study area

研究区位于浙西山地丘陵区，面积9.42 km2。地势总体北西为山地，中间为河谷平原，一条北西向常流大溪自北西往南东流出，为区内主要灌溉水系。出露地层由老至新为南华系、寒武系、奥陶系、志留系、石炭系、白垩系及第四系，以第四系地层为主，岩性主要为冲洪积物；次为寒武系地层，岩性主要为碳酸盐岩、碳质硅质岩和石煤层；奥陶系主要以泥岩和砂岩为主。北西向断裂为该区主要断裂构造，无侵入岩出露(图1)。土壤类型主要以红壤、水稻土为主，次为石灰岩土和粗骨土，其中潴育型水稻土分布于研究区中部，为水田主要土壤类型。土地利用类型主要有林地(有林地、灌木林地和其他林地)、水田、园地(果园、茶园、其他园地)、城镇及工矿用地、旱地等。耕地(水田和旱地)集中于中部，为本次重点研究对象，面积为3.25 km2。主要农作物为水稻、甘蔗、茭白等。研究区成土母质类型主要有洪冲积物、碳硅质岩类风化物、灰岩类风化物及砂砾岩类风化物。

石煤矿山位于北西向断裂北东角山坡上，其赋矿层位属下寒武统荷塘组黑色含矿岩系，20世纪50年代开采，2010年左右停采。目前已形成2个大废弃采矿坑，采场宕面裸露面积较大，长约200 m，宽约100 m，高10～60 m。宕底堆积大量石煤、含碳硅质岩废渣。矿山停采后未进行规范的生态恢复治理，自然复绿效果一般，堆渣裸露面积较大。

1.2 样品采集与分析测试

样品采集于2020年11月，采集的样品主要有表层土壤样、土壤垂向剖面样、地表水样、岩石样及农产品甘蔗样(图2)。表层土壤采样点全部布设于耕地，共采集144件土壤样品(不含8件重复样)，采样密度约为44件/km2，其中水田点位数占比83%，旱地点位数占比17%。表层土壤样由1个主坑和4个子样坑0～20 cm表层土壤组合而成，采样时避开田埂、明显人为污染等区域；土壤垂向剖面布设于矿山附近水田，共施工2条，剖面深度约1.2 m，采用浅钻(洛阳铲)开挖方式，按耕作层、潴育层、母质层等土壤发生层采集样品，每个剖面采集2～4个样品。每个样品质量不少于1 000 g，经自然风干后过20目筛并送样检测；岩石样采自石煤矿山荷塘组地层新鲜基岩，由多点采集新鲜基岩组合而成，共采集5件；农产品甘蔗以0.1～0.2 hm2为采样单元，按对角线法多点取样，等量混匀组成1个混合样品，样品总质量为1 000～2 000 g(鲜质量)，共采集甘蔗样20件；另外在矿山矿坑、小溪和灌溉水渠采集地表水样3件。

图2 研究区土地利用现状及采样点分布Fig.2 Status of land use and distribution of sampling points in the study area

样品测试均由中化地质华东分析测试研究中心完成，土壤样品分析指标有pH值、有机质、重金属(As、Cd、Cr、Cu、Pb、Zn、Ni、Hg)、Se等元素；岩样分析指标有重金属、Se等元素。pH测试采用玻璃电极法(Glass electrode method)，检出限为0.1；有机质测试采用容量法，检出限为0.1%；Cr、Cu、Pb、Zn测试采用X荧光光谱法(X-ray fluorescence spectrometry)，检出限分别为5 mg/kg、1 mg/kg、2 mg/kg、3 mg/kg；Cd和Ni测试采用电感耦合等离子体质谱法(Inductively coupled plasma mass spectrometry)，检出限分别为0.02 mg/kg、1 mg/kg；Hg测试采用冷蒸气-原子荧光光谱法(Cold vapor atomic fluorescence spectrometry)，检出限为0.5 μg/kg；Se、As测试采用氢化物发生-原子荧光光谱法(Hydride generation-Atomic fluorescence spectrometry),检出限分别为0.01 mg/kg、1 mg/kg；水样按《地下水质量标准》(GB/T 14848—2017)[18]等相关规范要求进行地表水重金属含量测定。测试时采用标准物质监控分析质量，经检查，所有元素相对偏差均小于20%，样品中各检测指标的准确度、精密度及报出率均符合《土地质量地球化学评价规范》(DZ/T 0295—2016)[19]的要求，数据可靠。

农产品按《生态地球化学评价样品分析技术要求》(试行)(DD2005—03)[20]生物样品分析方法要求，测试Se、As、Hg、Cr、Cu、Zn、Ni、Cd、Pb等指标。其中Cr、Cu、Zn、Ni、Cd、Pb测试采用电感耦合等离子体质谱法，检出限分别为0.2 mg/kg、1 mg/kg、1 mg/kg、0.1 mg/kg、0.03 mg/kg、0.1 mg/kg；Se、As、Hg测试采用原子荧光光谱法，检出限分别为0.005 mg/kg、0.3 mg/kg、0.005 mg/kg。另随机抽取5%样品送外检分析，原样与送检样相对偏差合格率达100%。

2 富硒土壤地球化学特征

2.1 表层土壤硒含量特征

对研究区144件表层土壤样品硒含量进行统计，由于样品数据不服从正态分布，故取循环剔除3倍离差后的算术平均值作为背景值，其结果见表1。经统计样品Se含量范围为0.18～3.71 mg/kg，平均0.87 mg/kg，背景值为0.72 mg/kg，均高于浙江省土壤平均值[21]和世界表层土壤平均值[22]。变异系数为0.66，属中等变异[23]，说明硒元素分布空间变异较大。由全硒含量的频率分布直方图(图3)可见，数据呈右偏型正态分布，样品硒含量集中分布于0.47～1.34 mg/kg之间，按0.40～3.00 mg/kg范围作为富硒土壤标准进行判定[24]，占样品总数的89.58%(129件样品)达到富硒土壤标准，说明该区富硒土壤资源丰富，分布面广。

表1 研究区土壤硒含量特征(n=144)

图3 表层土壤Se含量频率分布直方图Fig.3 Frequency distribution histogram of Se content in surface soil

根据谭见安对我国生态景观划分界限值[24]和《土地质量地球化学评价规范》(DZ/T 0295—2016)[19]对表层土壤硒含量和硒效应进行分级(表2)，结果显示：该区无缺硒(缺乏)和潜在缺硒(边缘)样品，足硒样品(适量)有13件，占比9.03%；富硒样品(高)129件，占比89.58%；硒中毒样品(过剩)2件，占比1.39%，表明该区土壤硒效应分级以富硒为主，硒含量分级主要为高。

表2 研究区土壤全硒含量及硒效应分级结果(n=144)

2.2 土壤垂向剖面硒含量特征

为了解土壤硒含量的深部变化情况，分析了2条位于矿山南部附近水田的富硒土壤垂向剖面，变化特征见图4。由图4可知，土壤垂向剖面各采样点硒含量总体较高，都达到富硒土壤标准。其中CSP16土壤垂向剖面的硒含量变化范围为0.38～1.06 mg/kg，平均值为0.77 mg/kg，在0～40 cm深度范围内硒含量相对富集，其平均值为1.06 mg/kg；CSP17土壤垂向剖面硒含量变化范围为0.54～2.24 mg/kg，平均值为1.39 mg/kg，在0～18 cm深度范围内硒相对富集，硒含量为2.24 mg/kg。2条土壤垂向剖面都呈现硒含量随深部增加而降低的规律性变化。表层土壤的硒含量均明显高于母质层，表明表层土壤硒除继承于母质层之外，还可能受到次生富集作用的影响。

图4 土壤剖面硒含量垂向变化Fig.4 Vertical change of selenium content in the soil profile

2.3 富硒土壤空间分布

从研究区土壤Se地球化学图(图5)可知，表层土壤硒高值区位于北东部大梗村、山背村一带，呈北东走向，总体与北东部石煤矿山走向一致，硒含量均在0.85 mg/kg以上，以矿山附近200 m范围内最高，含量最高达3.71 mg/kg；离矿山2.3 km南西部西坞村稍低，说明该区耕地内土壤硒含量较高与石煤矿山有直接关系，推测其形成与成土母质碳硅质岩类风化物有密切关系。采用土壤硒评价等级标准[19]，对相应采样地块图斑进行Se含量赋值，无样点分布耕地图斑按临近点位Se含量数据赋值，绘制研究区耕地表层土壤Se丰缺评价图(图6)。发现耕地区内(总面积3.25 km2)富硒土壤位于研究区中部的大梗村一带，面积2.84 km2，占耕地总面积的87.38%；硒含量适量土壤位于西坞村一带，面积0.36 km2，占耕地总面积的11.08%；少量硒含量过剩土壤位于石煤矿山附近，面积0.047 km2，占耕地总面积的1.45%。

2.4 土壤重金属含量特征

为查明富硒土壤重金属污染状况，统计了研究区As、Cd、Cu、Pb、Zn、Ni、Hg、Cr 8种土壤重金属含量特征，其结果见表3。由表3可知，Cd、Cu、As、Zn、Ni 5种元素中位数含量都小于均值，且相差较大，说明此5种元素数据中偏大数较多，中位数更能代表此5种元素的一般含量，但8种重金属元素平均含量和中位数都超过中国土壤背景值[25]。Cd、Cu、As、Zn、Ni 5种重金属元素含量采用中位数，Pb、Hg、Cr 3种重金属元素含量采用平均值，同中国土壤背景值比较，重金属元素超标倍数由大至小依次为Cd、Hg、Cu、As、Zn、Ni、Pb、Cr，说明该区土壤中8种重金属元素均有较高富集，其中Cd富集程度最大。

由变异系数统计可知，Cd、As、Ni、Cu、Zn变异系数都在1以上，属强变异；Hg、Pb、Cr元素变异系数在0.14～0.79之间，为中等变异。Cd变异系数最大，表明Cd元素分布的空间变异最大。

按《GB15618—2018土壤环境质量农用地土壤污染风险管控标准(试行)》[26]，计算重金属点位超标率(超过风险筛选值)，结果表明该区土壤重金属超标率由大到小依次为Cd>As>Cu>Ni>Zn>Hg，Pb、Cr不超标，其中Cd点位超标率达66.67%。说明该区土壤重金属污染总体比较严重。

图6 研究区耕地表层土壤Se丰缺评价图Fig.6 Evaluation map of Se abundance and deficiency in the surface soil of cultivated land in the study area

3 农作物中硒和重金属含量特征

石煤矿山附近耕地以大面积种植甘蔗为主，本次采集20件甘蔗样品，Se含量分析结果见表4。结果表明，甘蔗新鲜样品的Se含量范围为0.003～0.037 mg/kg，均值为0.017 mg/kg；Cd含量范围为0.02～0.51 mg/kg，均值为0.15 mg/kg。按照《富硒食品硒含量分类标准》(DB 36/T566—2017)[27]中水果食品富硒标准0.01～0.05 mg/kg对农产品进行评价，研究区甘蔗富硒率为65%，另外按《食品安全国家标准 食品中污染物限量》(GB 2762—2017)[28]评价标准，对水果类有标准要求的Cd、Pb进行重金属评价，结果表明：该区农产品甘蔗重金属超标严重，20件甘蔗样品有8件Cd超标，超标率为40%，另有1件样品Pb含量超标，农产品质量安全存在较大风险。

表3 石煤矿山周边表层土壤重金属含量特征(n=144)

表4 研究区甘蔗硒含量和重金属含量特征(n=20)

4 土壤硒含量的影响因素

4.1 地质背景

李家熙等[29]研究认为高硒岩石和煤层是高硒土壤中硒的主要来源，石煤层和含碳页岩中的硒含量相对较高，是形成富硒土壤的主要因素。刘道荣等[12]在常山地区调查发现：常山县西阳山组岩石硒含量均值最高，达0.90 mg/kg；其次为荷塘组硅质岩、碳质页岩夹石煤层，岩石中硒含量均值达0.59 mg/kg；休宁组砂岩岩石硒含量为0.11 mg/kg，均高于中国东部上地壳值(0.05 mg/kg)[30]。据本次矿山岩性采样，寒武系下统荷塘组地层中的碳质页岩、硅质岩岩石硒含量均值达4.39 mg/kg，富集系数达87.8。该区土壤富硒高值区与荷塘组地层分布在空间上基本一致，说明该地层是富硒土壤硒含量的主要来源，其次为碳酸盐岩地层。

统计不同地层对应的土壤硒含量(表5)进行研究分析可知，各地层土壤硒含量差异明显。本次采样点涉及南华系下统休宁组、寒武系下统荷塘组、寒武系中—上统、奥陶系中—下统、奥陶系上统三衢山组及第四系等地层，不同地层单元对应土壤硒含量由高至低依次为1h>2-3>O3s>Nh1x>Q>O1-2，土壤硒含量平均值最高为寒武系下统荷塘组，达1.31 mg/kg；其次为寒武系中—下统和奥陶系上统三衢山组地层，由于该区成土母质灰岩类风化物硒含量较高，成土作用后形成的土壤硒含量均值亦相对较高，分别为1.24 mg/kg、1.19 mg/kg；而后为南华系下统休宁组砂岩地层分布区，土壤硒含量均值为0.94 mg/kg；含量较少的为第四系和奥陶系中—下统地层，土壤硒含量分别为0.78 mg/kg和0.64 mg/kg。

采用本次采集的第四系土壤硒含量均值0.78 mg/kg、荷塘组地层硒含量均值4.39 mg/kg，和刘道荣等[12]在常山地区调查研究中得出的中—上寒武统地层硒含量均值0.90 mg/kg、下奥陶统印渚埠组地层硒含量均值0.06 mg/kg分别作为第四系、下寒武统荷塘组、中—上寒武统和下奥陶统印渚埠组地层硒平均含量界线，编制A-A’剖面土壤Se含量和地层Se含量对应分布图(图7)。从图7可知，荷塘组地层硒平均含量高于对应土壤硒含量，中—上寒武统地层和印渚埠组地层硒平均含量总体低于对应土壤硒含量，主要地层硒平均含量与其所对应的土壤硒含量总体上有较好的对应关系，荷塘组黑色岩系是影响该区富硒土壤的主要因素。

表5 研究区不同地层分布区土壤硒含量

图7 A-A’剖面土壤Se含量和地层Se含量对应分布图Fig.7 Corresponding distribution map of soil Se content and stratum Se content in section A-A’

4.2 土地利用类型

土地利用类型是土壤硒含量的影响因素之一。本次采集的样品都位于水田和旱地，统计这两种不同类型的土壤硒含量可知(表6)，水田和旱地土壤Se含量平均值都在0.9 mg/kg左右。对水田和旱地两种土地利用类型硒含量进行差异性T检验，其两者的差异概率P值为0.57(>0.05)，说明两者差异性不显著。

据刘道荣等[31]对浙西球川富硒区耕地土壤硒含量研究发现，主要由寒武系下统荷塘组黑色岩系形成富硒土壤，其中水田潴育型水稻土因透水、透气，受长期耕作影响，易使土壤活性态硒流失，造成土壤硒含量减少，而旱地土壤中硒更多继承成土母质，致使旱地平均硒含量显著高于水田。黄钊等[17]对云南省新平县哀牢山地区富硒土壤研究中发现该区旱地平均硒含量高于水田。本研究中未发现旱地与水田土壤硒平均含量存在显著差异，推测其原因是本次旱地采样点分布于水田四周，点位各地层分布不均，大部分旱地采样点位于地质背景硒含量较低的奥陶系中—下统、第四系等地层，致使旱地硒含量不高。

表6 不同土地利用类型土壤Se含量变化

4.3 土壤类型

不同成土母质、气候、地形高度、地理位置等因素产生不同土壤类型，因此土壤类型亦是影响土壤硒含量产生差异的重要因素。根据第二次土壤普查资料，该区土壤类型主要分为水稻土、红壤、石灰岩土、粗骨土4类，水稻土主要分布于中部，分布面积最大，红壤分布于外围，石灰岩土主要分布于区内东部山背村一带，粗骨土分布较少，主要集中于北部石煤矿山附近。Se含量按不同土壤类型统计结果见表7。区内硒含量平均值由大到小依次为粗骨土、石灰岩土、水稻土、红壤。粗骨土硒含量均值最高，达1.15 mg/kg，其主要原因是粗骨土成土母质为荷塘组的碳质页岩风化物，硒含量较高的成土母质是导致粗骨土土壤硒含量最高的直接原因。石灰岩土成土母质主要为灰岩类风化物，灰岩类岩性硒含量较高亦导致石灰岩土硒含量较高。刘道荣等[12]对浙西常山地区研究发现常山县红壤硒含量大于水稻土硒含量。该区水稻土硒含量高于红壤推测是由于北部硒含量较高的碳质硅质岩等黑色岩系、碳酸盐岩等原岩，经风化、淋滤作用，其土壤硒不断补充南部水田所致，而本次红壤采样样点较少，且大部分位于硒含量背景值较低的泥岩和砂岩分布区，亦是造成红壤硒含量较低的主要原因。

图8 土壤有机质含量与硒含量相互关系Fig.8 The relationship between soil organic matter content and selenium content

表7 不同土壤类型Se含量变化

4.4 pH、有机质对土壤硒含量影响及土壤硒含量与其他理化指标的相关性

表层土壤中pH和有机质等理化指标对土壤硒含量有着重要的影响，亦是控制有效硒的主要因素。土壤有机质含量越高，对硒吸附能力越强。研究区土壤有机质含量为0.64%～7.28%，平均2.51%，不同土壤类型有机质含量排列为石灰岩土>水稻土>红壤>粗骨土。从表8可知，有机质与硒含量相关系数达0.65，相关性较好。根据图8，无论水田和旱地都表现出土壤有机质与硒的极显著相关性(P<0.01)。这是因为富含有机质的土壤能增多土壤团粒结构，增大土壤比表面积，从而导致土壤吸附硒能力加强，有机质含量越高土壤吸附硒能力越强。统计表明pH值与土壤硒含量不相关；pH值按pH<5.5、5.5≤pH<6.5、6.5≤pH<7、pH≥7与对应土壤硒含量分别进行统计亦不相关。前人对云南哀牢山地区研究发现，富硒土壤区土壤硒含量与pH呈显著负相关[17]，临安东部土壤硒与pH呈正相关[32]，四川省沐川县北部土壤硒与pH不相关[33]，这可能是地区不同导致的差异所致。

此外，从表8可知，土壤Se含量与Cd、As、Hg、Cu、Pb、Zn、Cr、Ni亲硫重金属元素含量呈极显著的正相关关系(P<0.01)，而位于北部的石煤矿山赋矿层位为寒武系荷塘组，以富集V、Ni、Cr等元素而区别于其他地层(表9)，进一步说明土壤Se含量与荷塘组地层关系密切，对该地层形成成土母质具有一定继承性，这与浙西富硒土壤硒来源的研究相似[11]。

表8 土壤Se含量同其他理化指标相关性

表9 研究区荷塘组地层部分重金属元素含量(wB/(mg/kg))

4.5 外源输入对土壤硒含量影响

硒的不同来源组分可利用概率累积曲线的拐点进行初步判断[35]。从图9可知，该区表层土壤硒含量概率累积曲线根据最小二乘法线性拟合，分离出A、B两条拟合直线，且相关系数接近于1，表明Se元素含有A、B两种组分的物质来源，A组分为自然背景输入，拐点处累积贡献百分比达80%，推测受碳质硅质岩、碳酸盐岩等形成成土母质的影响，B组分直线斜率较低，样点较少，累积贡献百分比达20%，推测是受外源输入组分的影响。据刘道荣等[12]对整个常山县土壤硒的累积曲线分析，自然背景影响接近100%，外源输入组分比例总体偏少，而该研究区外源输入组分比例总体偏高。通过调查，该区附近除北部有过大面积石煤矿山开采外，未发现有大面积人工施用硒肥。目前矿山虽已停采，但复绿效果一般，宕面裸露面积仍旧较大。据野外石煤矿山矿坑及周边3处地表水采样分析(表10)，结果表明矿坑内地表水pH呈强酸性，Cd、Zn、Cu、Se含量明显高于周边地表水。对照《GB5084—2005农田灌溉水质标准》[36]，矿坑水Cd含量超标29.7倍、Zn超标3.1倍、Cu超标1.2倍、Se超标8.6倍，指示随雨水淋滤和地表水径流，矿山地表水中该4种元素由北西往南东随水流作用对周边灌溉水的较强参与作用。说明土壤Se含量除主要来源于自然背景外，还与大气沉降、矿山富Se废水汇入灌溉水，引入水田有直接关联。

图9 土壤硒含量概率累积分布Fig.9 Probabilitic cumulative distribution of soil seleniumcontent

表10 研究区地表水样品分析结果

5 结 论

(1)石煤矿山周边表层土壤Se含量的整体水平较高，Se含量高值区总体与北东向石煤矿山赋存位置一致，样品Se含量集中于0.47～1.34 mg/kg之间，均值为0.87 mg/kg，高于浙江省平均值，89.58%采样点位达到富硒土壤标准要求。在垂向上，土壤Se表现为表面富集，即从地表至深部土壤硒含量逐步降低。土壤硒效应总体以富Se为主，但存在土壤重金属超标风险。

(2)研究区甘蔗Se含量范围为0.003～0.037 mg/kg，均值为0.017 mg/kg，65%甘蔗样的Se含量达到富硒标准，但存在Cd等重金属元素超标，Cd超标率达40%，研究区农产品存在较大的质量安全风险。

(3)寒武系下统荷塘组碳质硅质岩及寒武系碳酸盐岩岩石的硒含量较高，其对应土壤Se含量由高至低依次为1h>2-3> O3s> Nh1x> Q>O1-2。土壤硒与地质背景关系密切，成土母质是土壤Se的主要来源；此外土壤硒还受土壤类型等因素的制约，呈现出明显不同的含量特征。不同土壤类型中，粗骨土硒含量均值最高，平均值达1.15 mg/kg。

(4)土壤中有机质高的土壤更易富硒，水田和旱地土壤硒含量与有机质均呈极显著正相关，表层土壤Se与pH无相关性，与Cd、As、Hg、Cu、Pb、Zn、Cr、Ni等重金属元素呈极显著的正相关关系，表明Se主要继承于成土母质。

(5)研究区土壤硒含量主要受控于自然背景，寒武系下统荷塘组碳质硅质岩和碳酸盐岩风化物是形成富硒土壤的主要因素，小部分与以往矿山开采引起的富硒废水输入等人类工程活动外源组分影响有关。

致谢：中化地质矿山总局浙江地质勘查院宋元青、袁野等同志参与了项目的野外工作，论文编写过程中得到浙江省地质调查院宋明义教授级高级工程师、魏迎春老师和浙江省林业科学研究院李贺鹏博士等的指导，审稿专家在该文修改过程中提出宝贵意见，在此一并表示衷心感谢。