滇东弥勒盆地褐煤的矿物质特征：内陆盆地泥炭堆积过程中海水输入的新证据

2022-01-19刘晶晶代世峰SPIROBaruch宋晓林魏建朋

煤炭学报 2021年12期

刘晶晶，代世峰,， SPIRO Baruch F，宋晓林，魏建朋

(1. 中国矿业大学国际煤地质学研究中心，江苏徐州 221116；2. 中国矿业大学(北京) 煤炭资源与安全开采国家重点实验室，北京 100083；3. 伦敦自然历史博物馆地球科学部，英国伦敦 SW7 5BD；4. 云南云投镇雄矿业能源开发有限公司，云南昭通 657000)

在沉积过程中受到海水影响的煤通常含有高含量的硫。继WHITE和THIESSEN[1]在研究美国Illinois煤盆地时，首次发现海水与沉积环境之间的关系后，已有很多学者对海相成因的煤进行了研究[1-6]。WILLIAMS和KEITH[7]在研究Lower Kittanning煤时提出，上覆海相地层的煤通常比正常的煤含有较高含量的硫，并且煤中硫的含量在很大程度上取决于受海水影响的程度[2,8-9]。例如，内蒙古乌达煤田的太原组的9号和10号煤，在其泥炭堆积时受到了海水的影响，全硫质量分数分别为3.46%和3.42%，受海水影响中等的12号和13号煤层，全硫质量分数分别为2.29%和0.88%[5]。煤中高质量分数的硫会导致一系列的环境问题。在煤燃烧过程中，煤中的硫有可能随着烟气的冷却形成SO3，然后转变为H2SO4，造成对燃烧炉的腐蚀[10]。同时，燃烧炉中释放的硫氧化物也是酸雨的主要来源[2]。煤中硫的主要载体矿物(黄铁矿)通常也富含其他有害元素砷(As)和汞(Hg)等[11]，在煤的利用过程中也会对环境和人体健康造成危害。另一方面，海相环境中形成的高硫煤也可能富集战略性关键元素，如铀(U)、钒(V)和硒(Se)，特别是稀土元素(Rare Earth Elements和Yittrium，REY)[12-13]，这些高度富集关键元素的煤可作为这些元素提取的潜在来源。基于上述2个方面，海相成因的高硫煤已在世界范围内引起学者的关注[2-6]。

中国西南地区，如广西合山、贵州贵定和紫云、云南砚山，晚二叠世海相碳酸盐岩层中沉积形成的煤通常称为超高有机硫煤(有机硫质量分数4%～11%)，煤中的有机硫主要来源于海底喷流[12,14]，海水的入侵不是主导因素[15-16]。受热液影响而导致煤中硫的富集，在其他地区也有报道。例如，TURNER 和RICHARDSON[17]在研究英格兰Northumberland煤田Westphalian A 和B煤时，发现煤中部分次生硫主要来源于富硫的热液，而该热液的供给源远离海水影响的区域。

成煤过程中受海水影响的煤层，通常具有特定的地质和地球化学特征。CHOU[9]提出，一般情况下，海相泥炭中黄铁矿硫平均占全硫的14%，比淡水中形成的煤层中黄铁矿质量分数高出一个数量级。莓球状黄铁矿通常形成于海相泥炭的早期成岩阶段[18-19]。我国西南地区晚二叠世的高硫煤或超高有机硫煤一般有海相碳酸盐顶板或沉积于局限碳酸盐台地[12,14-16,20-21]。煤中的有孔虫及腕足动物碎片、腹足类、介形虫和微生物化石的富集也指示了海相沉积环境[22]。GOODARZI和SWAINE[23]在研究澳大利亚和加拿大煤时，提出用B元素的质量分数评估受海水的影响程度：海相沉积的煤中B的质量分数>110×10-6，尽管其质量分数有时也受其他因素(如热液流体)的影响。另外，一些对沉积环境敏感的微量元素的质量分数以及它们的比值也可用来判别海水或淡水沉积环境。例如，受海水影响的煤中黄铁矿经常富集As，Se，Hg，Tl，Pb，Mo，Cd，Ni，Sb和Cu等元素[24]。SPIRO等[25]和DAI等[26]认为可以通过煤中同生的并且是海水来源的Sr的同位素，判断泥炭的堆积环境。

过去的研究通常认为，沉积于海相环境的煤具有海相碳酸盐顶板[14-16, 20-21]。但是，对于陆相环境中形成的煤或无海相碳酸盐顶板的煤，很难考虑到其在泥炭堆积期受到海水输入的影响。因此，通过分析煤的地球化学和矿物学特征，可以为此类型煤的形成提供是否受到海水影响的证据，并为煤盆地的区域构造演化历史提供煤地质学依据。本文的研究区为滇东新近纪弥勒盆地，长期被认为沉积于陆相沉积环境[27]。在对该盆地煤中锶同位素的研究基础上[25]，笔者通过分析弥勒盆地中新世煤的元素地球化学和矿物学特征，发现了海水入侵的新证据，以期进一步揭示该盆地泥炭堆积时所经历的地质过程。

1 研究区地质背景

云南位于印度板块之东，扬子地块之西南，同时又处于特提斯构造域与环太平洋构造域的交接复合部位。因此，云南区域构造域与印度板块、扬子地块和太平洋板块的活动状态以及3者之间的相互作用密切相关[27]。本次研究区为云南省弥勒县(图1，修改自SCHOENBOHM等[28])境内的山心村煤矿，北距省会昆明143 km，南距红河州府蒙自148 km，隶属于弥勒盆地南端的跨竹矿区。山心村煤矿东西平均宽约1 km，南北长约0.78 km，面积约0.78 km2。矿区含煤地层为中新世小龙潭组，厚度为40～150 m，平均厚度为100 m，主要由泥岩、砂质泥岩，以及少量的粉砂岩、炭质泥岩、煤层和非海相的碳酸盐组成。小龙潭组下伏地层为始新统—渐新统木花果组，其不整合覆于中三叠世个旧组之上，厚度为150～530 m，上部和底部为棕红色砂岩，中部为深灰色微现暗红色厚层状砾岩，钙质胶结，与下伏地层呈不整合接触；个旧组地层厚度为300～2 500 m，为浅海碳酸盐岩相沉积，主要由灰岩及少量白云质灰岩和白云岩组成，岩溶发育，为盆地矿坑充水的主要含水层。小龙潭组上覆地层为上新世师宗渡组，平均厚度为156 m，主要由砂砾岩、粉砂岩和泥岩组成。第四系沉积厚度为0～15 m，主要由黄褐色黏土、细砂及砾石组成。

图1 哀牢山—红河断裂带及研究区地理位置(修改自SCHOENBOHM等[28])

本次采样区为弥勒盆地跨竹矿区，采集了含煤地层底部的M1煤层，共采集样品61件，包括顶板、底板和2个夹矸样品，以及57个煤样品(表1)。从煤层顶部至底部，依次命名为ML-1～ML-61，顶、底板和夹矸分别附后缀-R，-F，-P加以表示(表1)。

表1 滇东弥勒煤中的矿物(石膏和黄铁矿，灰基，%)、主要元素的氧化物(CaO，MgO，P2O5，煤基，%)和稀土元素的质量分数(μg/g)，沉积环境敏感元素质量分数比(Ca/Mg，Sr/Ba，Th/U和V/Ni)以及锶同位素比

续表

2 分析方法

按照ASTM Standards D3174-11[31]测定煤中的灰分，按照ASTM Standards D3177-02[32]测定煤中全硫的质量分数。煤中元素C，H，N的质量分数由元素分析仪(Elementar，Vario MACRO)测得。

用X射线荧光光谱仪(ARL ADVANT′XP+)测定煤的氧化物及岩石样品中的常量元素的氧化物(CaO， MgO和P2O5)的质量分数。用电感耦合等离子质谱仪(Thermo Fisher，X series II ICP-MS)测定样品中的微量元素(Sr，Ba，Th，U，V，Ni和稀土元素)的质量分数。待测样品的预处理过程包括：称样、加酸、微波、消解、赶酸、提取和定容。对于煤和岩石(顶、底板和夹矸)样品，每50 mg样品，分别添加5 mL 65% HNO3，2 mL 40% HF和2 mL 65% HNO3，5 mL 40% HF 作为消解试剂，ICP-MS测试方法的检出限不高于0.02 ng/mL。

用X射线粉末衍射仪，结合商用软件SiroquantTM，对煤中的矿物(黄铁矿)进行定性和定量。用X射线荧光光谱仪(ARL ADVANT′XP+)测定煤的氧化物及岩石样品中的常量元素的氧化物(CaO，MgO和P2O5)。用电感耦合等离子质谱仪(Thermo Fisher，X series II ICP-MS)测定样品中的微量元素(Sr，Ba，Th，U，V，Ni和稀土元素)质量分数。用带能谱的场发射扫描电镜(FEI QuantaTM650 FEG)观察和分析煤中的矿物形态和显微结构，并分析样品中部分元素的分布特征。

在同位素比值87Sr/86Sr测试前，将已破碎至200目(0.075 mm)的煤粉末样置于超纯水中，充分搅拌以溶解煤中的石膏。过滤以分离固体物和石膏溶解液。采用Triton Plus质谱分析仪测试锶的同位素，锶的同位素质量分馏用87Sr/86Sr=8.375 209校正[33]。国际标准样品NBS987测试的平均值87Sr/86Sr=0.710 250±0.000 010。

3 结果

3.1 原煤中的矿物组成

由于有机质的影响，所研究的57个煤样中仅有11个煤样品的XRD分析结果理想(表1)。弥勒原煤中的矿物主要为高岭石、伊利石、石膏、黄铁矿和锐钛矿。在灰分高于40%的样品中有金红石和绿泥石。部分样品中含有少量纤磷钙铝石、绿泥石和方解石。石膏的质量分数变化范围较大(为1.80%～70.10%，灰基)，均值约为30%，个别煤分层样品(ML-49)中质量分数高达70%(表1)，远高于其他地区煤中石膏的质量分数[34-35]。下部煤分层中石膏质量分数比上部煤分层中的高。

图2 弥勒煤中石膏和黄铁矿的扫描电子显微镜背散射图像

3.2 元素地球化学特征

3.2.1常量和微量元素

研究区煤中硫的质量分数变化为0.19%～3.78%，均值为0.78%(干燥基)(表1)。高硫含量(>1%)的分层位于剖面的底部和上部。2个连续的煤分层ML-48和ML-49中全硫的质量分数分别为3.78%和3.60%。在这些样品中，有机硫和硫酸盐硫是硫的主要赋存形式，分别占全硫质量分数的40%～60%和25%～42%。

煤中部分常量元素氧化物、微量元素的含量以及锶同位素比值87Sr/86Sr见表1。CaO的质量分数为1.01%～4.76%，均值为2.97%(煤基)，高于中国煤中CaO质量分数的均值(1.23%[38])。MgO的质量分数为0.28%～1.14%，均值为0.69%，高于中国煤中MgO质量分数的均值(0.22%[38])。

微量元素Sr的质量分数为61.73～1 011.28 μg/g，均值为204.4 μg/g，高于世界低阶煤中Sr的平均质量分数(120 μg/g[39])。Ba的质量分数均值为137.01 μg/g，与世界低阶煤中的均值(150 μg/g[39])相当。微量元素Th和U的质量分数均值分别为2.30，1.78 μg/g，均低于世界低阶煤的均值(分别为3.30,2.90 μg/g[39])。

3.2.2 稀土元素

稀土元素包括镧系元素(La，Ce，Pr，Nd，Sm，Eu，Gd，Tb，Dy，Ho，Er，Tm，Yb和Lu)和Y，共15个元素。按照SEREDIN和DAI[40]提出的稀土元素分类方法，将稀土元素分为轻稀土(LREY:La，Ce，Pr，Nd，Sm)、中稀土(MREY:Eu，Gd，Tb，Dy和Y)和重稀土(HREY:Ho，Er，Tm，Yb和Lu)。利用上地壳中稀土元素的质量分数对研究区样品中稀土元素的质量分数进行标准化后，将样品中稀土元素的富集类型划分为3种：轻稀土富集型(LaN/LuN>1)、中稀土富集型(LaN/SmN<1，GdN/LuN>1)和重稀土富集型(LaN/LuN<1)。据BAU等[41]和DAI等[42]提出的计算方法，对Ce，Eu，La，Gd和Y的异常按照式(1)～(4)计算。

(1)

(2)

(3)

(4)

其中，下标N为用上地壳的稀土元素含量进行标准化后的值；*为标准化后稀土元素含量的期望值，异常值大于1和小于1分别表示稀土元素的正异常和负异常，等于1表示无异常。

图3 弥勒煤及顶底板、夹矸中稀土元素的分布特征

3.2.3 石膏中锶同位素87Sr/86Sr特征

由于Sr与Ca具有相似的离子半径，石膏的矿物晶格中可能有部分Ca被Sr替换。因此，同沉积的自生石膏中Sr的同位素值87Sr/86Sr能够反映泥炭堆积时期溶液的性质。弥勒盆地13个煤样中石膏的87Sr/86Sr比值变化为0.708 350～0.708 591(表1)，它在煤层剖面上的变化如图4(a) 所示。从煤层顶部至底部，可将该数据点依次分为A，B和C三组，并且整体上呈递减趋势(图4(a))。

4 讨论：泥炭堆积受到海水影响的证据

4.1 矿物学证据

高含量的同生黄铁矿通常是受海水影响的煤的显著特征之一[1,3,7,45-50]，虽然在煤层形成的不同阶段热液成因的黄铁矿也经常被发现。在不同海相沉积物中，莓球状黄铁矿的形成机制相似，均是由四方硫铁矿(FeS0.9)经硫复铁矿(Fe3S4)转变形成。而黄铁矿单个晶体颗粒则一般由FeS与聚硫化物作用形成[51-53]。在上述黄铁矿的2种形成机制中，海水中的硫酸盐是各作用过程中所需的H2S的主要来源。因此，弥勒煤中丰富的同生黄铁矿，表明在泥炭堆积过程中受到了海水的影响。

弥勒煤中夹矸和底板中石膏的质量分数均小于3%，远低于煤层中石膏的质量分数。沉积岩中的石膏为蒸发作用的产物，湖盆和海盆中的卤水在干旱的气候条件下，经蒸发和浓缩作用后结晶而成。煤中的石膏通常被认为是次生矿物[36-37]。SPIRO等[25]通过对弥勒煤中石膏赋存状态的研究，认为其为同生沉积，来源可能为海水，经泥炭失水、压实后结晶而形成。在底部的煤层中(ML-48，ML-49，ML-50)，石膏的质量分数远高于顶部和中部煤层(表1)，表明底部煤层受海水影响程度较大，表明随着泥炭的堆积，海水的输入逐渐减少。

4.2 元素地球化学证据

除了上述黄铁矿外，煤中有机硫的性质能够为其沉积环境提供一定的信息[54]。成煤过程中海水的输入为泥炭提供了充足的硫酸盐作为硫源，因此，煤中丰富的黄铁矿和有机硫指示了成煤的海相沉积环境[3,7-9,55]。CHOU[55]在研究伊利诺伊盆地赫林煤时，发现海水向盆地的输入导致形成煤中高含量的硫及微量元素B，Mo和U，并终止了泥炭的堆积。在弥勒煤中，底部的煤分层中全硫的质量分数高于顶部和中部的煤分层，从底部至顶部硫的质量分数逐渐降低(图5)，推测是海水逐渐退出煤盆地所致。在硫质量分数高于1%的煤样中，有机硫和硫酸盐硫分别占全硫的50%和40%，也表明了其具有海相成因的特点[9]。

图5 弥勒煤的灰分Aad(%)，全硫St(%)、P2O5和CaO的质量分数(%)，以及Ca/Mg，Sr/Ba，Th/U，V/Ni和Y/Ho质量分数比在剖面上的变化

由于海水中的Mg2+质量分数远远高于Ca2+[56]，Ca/Mg质量分数比可用作沉积环境的指示剂，反映沉积物水体的盐度。水体的盐度越高，Ca/Mg质量分数比越低[57-58]。本次研究的弥勒煤中Ca/Mg质量分数比为5.26(表1)，与欧洲莱茵河附近受海水影响的褐煤中Ca/Mg质量分数比6[59]相当，也表明弥勒盆地的煤层在形成过程中受到了海水的影响。

沉积物中Sr和Ba的质量分数与水体的盐度有关[60-63]。全球现代海洋中Sr和Ba的平均质量浓度分别为8.0，0.013 mg/L，河流中平均质量浓度分别为0.07，0.02 mg/L[64]。因此Sr/Ba质量分数比可用来判别沉积环境的介质溶液盐度[65-66]。对于煤的沉积环境，该值可用以识别其受到淡水、海水或者两者的共同影响[25-26]。

弥勒样品(煤、夹矸、顶底板)中Sr/Ba质量分数比与Ad的相关性(图6)表明，低灰分(<40%，多数低于22%)样品中的Sr/Ba质量分数比>1，而高灰分样品(44%～80%)中Sr/Ba质量分数比<1，并且Sr/Ba质量分数比与灰分呈一定的负相关性。Sr/Ba质量分数比=1为煤层剖面上指示海水/淡水影响的临界值(图6)。因此，Sr/Ba质量分数比>1的含石膏的煤层(Ad<40%)受到了海水的影响，而Sr/Ba质量分数比<1的底板、顶板和夹矸(Ad>40%)形成主要受淡水的影响。由于海水中Sr的质量浓度比Ba的质量浓度高出100倍[61]，因此，受淡水影响的沉积环境中，只要有少量海水的输入，就会导致该沉积环境介质溶液具有海水87Sr/86Sr特征[25-26]。此外，在煤层剖面上，Sr/Ba质量分数比和CaO的质量分数在下部煤层中较高，随后至煤层上部显著降低(图5)，反映了海水输入影响减弱，而淡水输入影响增加的过程。

图6 弥勒煤中Sr/Ba质量分数比与煤的灰分(%)的关系(修改自SPIRO等[25])

煤中具有高质量分数P2O5的层位(位于1 394～1 546 cm和385～513 cm处)也具有高的CaO，Ca/Mg和Sr/Ba质量分数比(图5)，因此高质量分数的P2O5与海水的强烈影响有关[67]。

Th/U比值可用作沉积环境的指示剂[68]。通常情况下，陆相风化产物中Th/U质量分数比较高(>7)，海相黑色页岩和灰岩中该元素质量分数较低(<2)[69]。BOUKA[70]指出，陆相沉积物中Th/U质量分数比>7，海相沉积物中Th/U质量分数比<7。因此，受海水输入影响较大煤中的Th/U质量分数比应低于未受海水影响或影响较小的煤层。在弥勒煤中，Th/U质量分数比为0.05～4.34，均值为1.43(表1)。从剖面底部至顶部，随着海水作用的逐渐减弱，Th/U质量分数比也呈现递增的趋势。在煤层深度1 300～1 400 cm处(图5)，Th/U质量分数比骤然降低，与P2O5的质量分数以及Ca/Mg，Sr/Ba的质量分数比在该处出现的拐点一致，表明在该成煤阶段，海水对煤层的影响减弱的程度最为显著。同时，由于海水中相对富集可溶性氧化铀的化合物[68]，因此在个别煤分层(ML-6，ML -11，ML -59)中U相对富集(最高为14 μg/g)，也是泥炭层中海水输入的结果。

与Ni相比，V更倾向富集于贫氧的海水中[71-72]，因此，海水中的V/Ni质量分数比高于淡水[73]，并且，沉积岩中的V/Ni质量分数比可用于氧化-还原沉积环境的判别[26,71,74-80]。在弥勒煤层剖面上，V/Ni质量分数比从下部至上部逐渐降低，也指示了沉积环境由还原性海相环境向氧化性的淡水环境的转变。

由于水体中稀土元素Y和Ho的清除速率不同，特别是在海水中，Ho会被海洋颗粒物优先清除，导致海洋沉积物中Y/Ho质量分数比显著高于硅质碎屑沉积物[81-83]。因此，Y/Ho质量分数比可用于判断硅质碎屑沉积物(Y/Ho质量分数比25～30)和海相沉积物(Y/Ho质量分数比60～70)[26,81,83]。弥勒煤中Y/Ho质量分数比为18.67～64.00(表1)，表明泥炭堆积过程中既有硅质碎屑的供给，又有海水的影响。Y/Ho质量分数比最大值位于煤层深度1 300～1 400 cm处(ML-34煤分层)，随后比值显著降低(图5)，指示陆源碎屑供给加强、海水作用的减弱，与本文中所讨论的其他参数异常特征一致。

在稀土元素方面，弥勒煤显示出La和Gd的正异常特征。La正异常是现代海水典型的稀土元素特征之一[84]，弥勒煤中La的正异常可能指示了海水的影响。中国大部分煤通常具有较弱的Gd负异常特征[42]。煤中的Gd正异常主要受海水、热液和其他水体的影响[42]。本次研究的煤中尚未发现热液流体的存在以及受热液流体形成的矿物，因此该煤中Gd正异常可能是由于海水的输入所致。

4.3 锶同位素证据

HODELL等[85]测定了深海钻探项目(Deep Sea Drilling Project，DSDP)3个钻孔(519，588和607钻孔)261件保存完好的浮游有孔虫化石样品的87Sr/86Sr比值，并结合磁地层学、氧同位素地层学、生物地层学和稳定同位素地层学等手段精确测定的地层年代，绘制了24 Ma至今海水中87Sr/86Sr比值的变化曲线(图4(b))。弥勒煤中石膏的87Sr/86Sr比值与中新世早期海水中87Sr/86Sr变化一致，且呈递增的变化趋势，表明形成石膏的母质溶液来源于海水。煤中A，B和C组87Sr/86Sr比值的均值在HODELL等[85]绘制的曲线上对应的地质年代分别为18.6，20.5，21.7 Ma(图4(b))，表明研究区泥炭的堆积始于早中新世21.7 Ma，终止于18.6 Ma，整个沉积过程持续了约3 Ma。

弥勒盆地煤中的87Sr/86Sr比值对应的地质年龄(21.5 ～18.5 Ma)恰好处于特提斯海域的关闭和哀牢山—红河断裂带碰撞变形区形成的时间21 Ma[86]或22～17 Ma[87]时间段内。同时，同一地质历史时期内，研究区煤中的87Sr/86Sr比值与海水中87Sr/86Sr比值一致，表明本次研究的石膏中的锶同位素值并未受到其他因素的影响(如与盆地基底三叠系个旧组灰岩之间发生水-岩反应)。因此，该锶同位素比值能够反映出研究区气候和沉积环境的变化。

4.4 研究区海水的输入机制

新生代以来，印度板块持续向北推移，欧亚板块则被动向南推挤。由于云南位于印度板块东侧，因此其仅显示向南的滑移推挤。印度板块北北东向的推挤，一方面使青藏地块东部向东和南东向扩展移动；另一方面，由于其西北端相对东北端向北运动缓慢，使印度板块本身发生逆时针方向旋转，产生向东的侧向挤压，形成哀牢山—红河左行走滑断裂带[88]。另外，印度板块东侧地块向南滑移的过程中，滇中地区表现出比两侧地块更大的滑移速度，导致沿东侧鲜水河、小江断裂带上发生明显的左行位移，形成当前的鲜水河—小江断裂系，并与哀牢山—红河断裂带交汇。本文研究区即位于该断裂系内[28]。与此同时，太平洋板块的俯冲方向由北北西转为北西西，促使原北东、北北东方向的断裂(如富源—弥勒、宣威—寻甸断裂带、弥勒—师宗断裂带)，转为右旋剪切运动，并伴随着哀牢山—红河断裂带北部区域(包括本研究区弥勒盆地)的抬升。

弥勒盆地为山间内陆盆地，其与哀牢山—红河断裂带中的特提斯海域之间的水系连接可能为变形断裂带。泥炭堆积期，长期处于多阶段活动期的断裂带变形区，为海水输入至泥炭沼泽提供了通道。哀牢山—红河断裂带北部区域的抬升，使弥勒盆地煤层形成时受到海水的影响逐渐减弱。