绘制环境中化学元素的区域分布地图对于确定和理解大规模地球化学过程以及记录空间和时间地球化学变化至关重要。同时，地球化学制图在土地利用规划和政策制定方面有许多实际应用，且在了解化学元素在地方病发生和分布中具有重要作用^[1-4]。1984年，国际原子能机构(IAEA)首次提出了“全球一张地球化学图”的概念；1988—1997 年，国际地球化学填图计划和全球地球化学基准计划相继开展^[5]。此外，各种具有特定目标的地球化学调查活动被广泛开展，如EuroGeoSurveys农业和牧场土壤地球化学制图项目(GEMAS)，该项目旨在获取欧洲农牧区土壤的金属含量数据^[6-7]；在中国东部第四纪平原开展了为环境调节和农业实践提供数据的多用途生态地球化学制图项目^[8]；巴伦支海生态地球化学项目^[9]；北美土壤景观项目等^[10]。这些调查通过地球化学制图以获取环境中有害微量元素的含量、分布和影响因素等信息。

煤炭是保证全球经济增长的主要能源之一，未来很长一段时间将继续保持其在全球能源生产中的主导地位。中国是全球最大的煤炭生产国，2020 年原煤产量为39亿t，约占全球煤炭产量的51%。由于煤炭工业规模庞大，其开发利用过程造成的环境污染问题备受关注^[11-12]，重点问题之一是煤中有害微量元素对环境构成的潜在风险^[13-15]。煤中有害微量元素向地下水、空气和土壤的迁移已被广泛报道^[16-18]。中国已经对有害微量元素开展了大量研究，一般集中在小规模矿区或局部地区。在全国范围内，较为完整的中国煤中有害微量元素含量分布图少有报道。笔者利用ArcGIS技术绘制中国含煤区中Be、Co、Cu、Mo、Th、Zn的含量分布图，以期为煤中微量元素和地球化学过程研究提供科学信息，并为煤炭环境管理提供直观有效的参考。

1   研究方法

1.1   数据来源

本文所使用的分析与制图数据来源于中国煤中微量元素数据库(Trace Elements in Coal of China Database，TECC)，包括1 167个Be、1 315个Co、1 406个Cu、1 191个Mo、1 247个Th和1 390个Zn样品数据(图1)。TECC是第一个旨在管理中国煤中微量元素数据的数据库系统，收录了中国主要煤矿的地理坐标、微量元素含量、工业分析等数据单元^[19-20]。这些数据主要收集自Web of Science、中国知网 (CNKI)文献数据库和相关专著。中国含煤区域的矢量边界绘制依据为开放源地理空间基金会(OSGeo)−中国中心(http://www.osgeo.cn)网站的“中国煤炭资源分布在线地图(1∶3200万)”^[21]。

图  1  中国煤炭样品中 Be、Co、Cu、Mo、Th 和 Zn 的含量频率分布

Figure  1.  Concentration frequency distribution of Be, Co, Cu, Mo, Th and Zn in Chinese coals

下载: 全尺寸图片 幻灯片

1.2   数据处理与成图

采用箱线图分析获得原始含量数据集的四分位数和异常值。四分位数可反映数据的离散程度，计算过程如下：首先对给定的数据进行排序，排序后的数据范围为X_i, i=1,2,···, n。随后，第1个四分位数(Q1)、第2个四分位数(Q2)和第3个四分位数 (Q3) 都计算为$ ({1 - {g}}){X_{i}} + g{X_{{i} + {\text{1}}}} $，其中i是等式的整数部分，g是等式的小数部分；(n+1)/4，(n+1)/2和3(n+1)/4分别是用于计算 Q1、Q2和Q3的公式。四分位距(IQR) 计算式为IQR=Q3−Q1。相应地，上限和下限分别使用Q3+1.5IQR和Q1−1.5IQR计算。结合实际情况，将下限指定为0，即超出[0,Q3+1.5(Q3−Q1)]范围的数据被定义为异常值。

各元素的含量数据经过对数(lg)转换后，呈现正态分布或近似正态分布(图2)。利用ArcGIS 9.0软件将lg转换后的含量数据和地理坐标整合到地理信息系统中，并运用反距离权重方法绘制中国煤中微量元素含量的空间分布图。地图中的颜色等级根据样本数据的占比进行分类，即：0~10%、10%~25%、25%~50%、50%~75%、75%~90%和大于90%。异常点不参与制图，但其空间位置在地图上标出。

下载: 全尺寸图片 幻灯片

  2  煤中 Be、Co、Cu、Mo、Th和Zn含量数据转换箱线图

  2.  Box plot representing the lg-transformed concentration distribution of Be, Co, Cu, Mo, Th, and Zn in Chinese coals

下载: 全尺寸图片 幻灯片

2   结果与讨论

2.1   煤中有害微量元素的含量特征

中国煤中Be、Co、Cu、Mo、Th、Zn含量均表现出较大的跨度，且数据分布呈正偏性，不符合正态分布(图1)。大部分数据点聚集在相对较低的值上，标准偏差(σ)和变异系数(cv)较高(表1)。

表  1  样品数据统计指标

Table  1.  Statistical indicators of sample data

指标		Be	Co	Cu	Mo	Th	Zn	备注
样品数量		1167	1315	1406	1191	1247	1390	原始数据
标准差σ		31.17	17.10	37.44	47.67	18.79	58.69
变异系数cv		4.10	1.88	1.31	3.94	2.07	1.34
含量范围和算术平均值/(mg·kg−1)		0~542/7.61	0~250/9.11	0~559/28.66	0~665/12.10	0~581/9.08	0~865/43.70
Q3 (25%)/(mg·kg−1)		1	1.90	8.10	0.90	2.58	12.85
Q2 (50%)/(mg·kg−1)		2	4.30	16.71	1.85	5.69	24.76
Q1 (75%)/(mg·kg−1)		3.63	10.25	36.77	4.91	12.13	50.47
样品数量		1002	1190	1306	984	1162	1255	剔除异 常值后 的数据
标准差σ		1.65	5.05	18.05	2.12	6.27	24.49
变异系数cv		0.78	0.91	0.85	0.97	0.85	0.82
含量范围和算术平均值/(mg·kg−1)		0.05~7.51/2.10	0.04~22.73/5.53	0.06~78.5/21.36	0.01~10.6/2.19	0.04~26.1/7.35	0.38~106/30.02
元素平均含量 /(mg·kg−1)	Dai Shifeng等[22]	2.11	7.08	17.5	3.08	5.84	41.4
	任德贻等[23]	2.13	7.05	18.35	3.11	5.81	42.18
	唐修义等[24]	1.90	7.00	13.00	4.00	6.00	38.00
	白向飞等[25]	1.75	10.62	17.87	2.70	5.88	—

下载: 导出CSV 

| 显示表格

基于箱线图分析排除异常含量值后，识别出中国煤中Be、Co、Cu、Mo、Th、Zn的一般含量范围分别为0.05~7.51、0.04~22.73、0.06~78.50、0.01~10.60、0.04~26.10、0.38~106.00 mg/kg，分别占样本总量的86%、90%、93%、83%、93%和 90%。排除异常值后，σ、cv、算术均值显著降低。这意味着即使分析中包含相对较少的异常值，统计结果也会受到显著影响(表1)。在不考虑异常值的情况下，中国煤中 Be、Co、Cu、Mo、Th和Zn含量的算术平均值与Dai Shifeng^[22]、任德贻^[23]、唐修义^[24]、白向飞^[25]等的计算结果接近。因此，在区域和全国范围内评价煤中微量元素的平均含量水平时，有必要筛选出异常值，以避免造成较大的分析误差。

根据本文计算结果，我国煤中Mo的平均含量为2.19 mg/kg，与全世界煤中Mo的平均水平(2.2 mg/kg)^[26]相当；此外，相比于世界煤炭中Be(1.6 mg/kg)、Co(5.1 mg/kg)、Cu(16 mg/kg)、Th(3.3 mg/kg)和Zn (23 mg/kg)的平均含量^[26]，我国煤中Be、Co、Cu、Th、Zn的平均含量更为富集，平均含量分别为2.10、5.53、21.36、7.35和30.02 mg/kg。

2.2   煤中有害微量元素含量的空间分布

根据聚煤作用特征及聚煤盆地的演化形成，我国煤炭资源分布被划分为5个赋煤区：东北赋煤区、华北赋煤区、西北赋煤区、华南赋煤区以及滇藏赋煤区(图3)。就成煤时代而言，晚石炭世−早二叠世(C₂−P₁)、晚二叠世(P₃)、晚三叠世(T₃)、早−中侏罗世(J_1-2)、晚侏罗世和早白垩纪 (J₃−K₁)，以及古近纪和新近纪 (E−N)是我国6个主要煤炭成煤期。其中，可采煤炭储量主要为J_1-2、C₂−P₁和J₃−K₁煤，分别占煤炭总储量的39.6%、38.1%和12.1%；P₃、T₃和E−N煤的占比相对较小，分别为7.5%、0.4%和2.3%^[22]。就地域而言，C₂−P₁煤主要分布在华北地区，J_1-2煤主要分布在中国西北地区，东北地区以J₃−K₁、E−N煤为主，南部地区以P₃、T₃煤为主，少量的E−N和T₃煤分布在西藏−云南西部地区。

图  3  中国含煤区和成煤时代分布

Figure  3.  Distribution of coal-bearing areas and coal-forming periods in China

下载: 全尺寸图片 幻灯片

本文绘制了中国煤中Be、Co、Cu、Mo、Th和Zn含量的空间分布(图4—图9)。以地域进行统计而不区分成煤时代的情况下，根据本文计算得到的我国煤中Be、Co、Cu、Mo、Th和Zn的平均含量作为衡量基准(表1)，各元素含量的空间分布情况如下：① 对于Be元素， 安徽、四川、黑龙江、重庆、云南、湖北等地煤中的平均含量较高，高于全国平均水平；山西、新疆、湖南、宁夏等地平均含量在全国平均水平附近；其余省份/自治区的含量低于全国平均水平(图10)。② 对于Co元素，高于全国平均含量的煤主要分布在四川、云南、吉林、贵州、河南、宁夏、广西、重庆等地；安徽、黑龙江、山东、湖南、新疆等地的煤接近全国平均水平；其余地区的含量相对较低。③ 对于Cu元素，云南、贵州、四川、吉林、河南、湖北等地煤中的含量较高；湖南、重庆、陕西、安徽、山西、河北等地处于中等水平；新疆和辽宁地区煤中Co元素较低。④ 对于Mo元素，安徽、辽宁、云南、宁夏、广西等地煤中含量较高；陕西、河北、河南等地的含量低于全国平均含量；其余地区处于中间水平。⑤ 对于Th元素，广西、云南、内蒙古、四川、江西、吉林、重庆等地的含量较高，且平均含量基本相当；河北、山西、陕西、河南、贵州、安徽、山东、湖南、黑龙江等地的含量依次降低；宁夏、辽宁、湖北和新疆等地煤中Th含量明显低于全国平均含量。⑥ 对于Zn元素，高于全国平均含量的煤主要分布在江西、云南、四川、重庆、广西等地；河北、新疆、山东、黑龙江等地煤中Zn含量远低于全国平均含量；其余地区处于中间水平且含量较为接近。

图  4  中国煤中Be的空间分布

注：图例中色柱两侧数值含义：≤7.51 mg/kg的样品占比88%；≤3.60 mg/kg的样品占比75%；≤2.00 mg/kg的样品占比50%；≤1.00 mg/kg的样品占比25%；≤0.51 mg/kg的样品占比10%；图5—图9含义相同

Figure  4.  Spatial distribution map of Be in Chinese coals

下载: 全尺寸图片 幻灯片

图  5  中国煤中Co的空间分布

Figure  5.  Spatial distribution map of Co in Chinese coals

下载: 全尺寸图片 幻灯片

图  6  中国煤中Cu的空间分布

Figure  6.  Spatial distribution map of Cu in Chinese coals

下载: 全尺寸图片 幻灯片

图  7  中国煤中Mo的空间分布

Figure  7.  Spatial distribution map of Mo in Chinese coals

下载: 全尺寸图片 幻灯片

图  8  中国煤中Th的空间分布

Figure  8.  Spatial distribution map of Th in Chinese coals

下载: 全尺寸图片 幻灯片

图  9  中国煤中Zn的空间分布

Figure  9.  Spatial distribution map of Zn in Chinese coals

下载: 全尺寸图片 幻灯片

图  10  不同省(自治区、直辖市)煤中Be、Co、Cu、Mo、Th和Zn的平均含量

Figure  10.  Average concentration of Be, Co, Cu, Mo, Th and Zn in coals in different regions

下载: 全尺寸图片 幻灯片

2.3   煤中有害微量元素含量空间分布格局的成因

任德贻等^[27-28]将煤中有害微量元素的富集成因类型初步划分为5类，包括陆源富集型、沉积−生物作用富集型、岩浆−热液作用富集型、深大断裂−热液作用富集型和地下水作用富集型。Dai Shifeng等^[22]通过对中国煤中微量元素深入分析，将其富集成因类型划分为物源区控制型、海洋环境控制型、热液控制型(包括岩浆控制型、低温热液控制型和海底喷气控制型)、地下水控制型和火山灰控制型。刘桂建等^[29]结合前人研究结论，将煤中微量元素富集成因按照时间进程分为原生因素(物源区母岩、成煤植物、沼泽水介质、气候、古土壤)、次生因素(岩浆热液、煤化程度)和后生因素(构造运动、风化作用、地下水作用、围岩物质交换)3种类型。由此可知，中国煤中有害微量元素含量空间分布格局的形成是长时期、多阶段、多因素综合作用的结果。Cao Qingyi等^[30]通过对比中国煤、上地壳和深层土壤(>1 m)中有害微量元素含量特征，指出有害微量元素平均含量在三者中的丰度具有相同的变化规律，并且针对同一元素，其在3种介质平面空间上的含量分布特征也是类似的，表明煤、岩石和土壤中微量元素的丰度分布同步受到岩石圈内元素地球化学循环的影响。此外，煤中有害微量元素的异常富集点与断层带分布具有空间关联性，并且异常富集点位置往往伴随着热液作用的痕迹^[30]，这指示断裂构造和热液作用可能作为煤中有害微量元素的运移通道和物质来源，热液作用可作为煤中有害微量元素异常富集的典型特征。

3   结 论

a. 基于TECC数据库，利用ArcGIS技术绘制了中国煤中Be、Co、Cu、Mo、Th、Zn含量的空间分布地图。这些地图有助于快速了解元素的空间分布信息，可以为煤中有害微量元素的潜在环境风险研究和煤炭环境管理提供参考。中国煤中有害微量元素含量的空间分布极不均匀。煤中Be、Co、Cu、Mo、Th、Zn的含量主要分布在0.05~7.51、0.04~22.73、0.06~78.50、0.01~10.60、0.04~26.10和0.38~106.00 mg/kg，平均含量水平分别为2.10、5.53、21.36、2.19、7.35和30.02 mg/kg。

b. 煤中有害微量元素含量的空间分布格局的形成是长时期、多阶段、多因素综合作用的结果。典型影响因素包括物源区母岩、热液作用、水运移作用等。煤中有害微量元素含量分布受到岩石圈元素地球化学循环的影响。热液作用是煤中有害微量元素异常富集的典型特征。

c. 尽管TECC数据库中已包含了丰富的中国煤中有害微量元素的数据，但尚不能完全覆盖我国所有的含煤地区。来自于内蒙古、山西、陕西、安徽、贵州、广西等地区的样本数据较多，而来自于西北地区和滇藏地区的数据相对较少。今后仍需要进行持续补充分析测试工作。此外，本文中目标元素含量的空间成图范围采用的是含煤区范围，而非煤田范围。这是由于煤田范围较为分散，不利于空间数据展示。因此，今后需要对空间数据模型做进一步优化，提高空间数据的展示精度，使元素含量的空间分布更加符合实际情况。