硅质页岩中SiO2的成因主要包括陆源、生物、热液三种[29]。根据前人研究,陆源输入通常与Al元素的含量相关,而热液的参与则影响Fe、Mn元素的富集[30⁃31],因此常用Al-Fe-Mn三角图和Al/（Al+Fe+Mn）比值来判断硅质页岩的成因。Al/（Al+Fe+Mn）比值通常反映热液对沉积物的贡献量,热液物质的贡献量越多,Fe、Mn元素越富集,比值就越小。通常纯热液沉积物的Al/（Al+Fe+Mn）比值近似0.01,受热液影响的沉积物Al/（Al+Fe+Mn）比值通常小于0.35,当比值大于0.4时,表明主要受陆源输入的影响,而纯生物作用形成的硅质岩该比值接近0.6[32⁃33]。例如东太平洋洋中脊纯热液沉积物的Al/（Al+Fe+Mn）比值低至0.01,半远洋环境下纯生物成因的Kamiaso硅质岩的比值为0.60[32]。研究区五峰组—龙马溪组Al/（Al+Fe+Mn）比值介于0.42~0.78,表明该套页岩硅质成分受热液影响较小,生物成因和陆源输入可能是该套页岩SiO2形成的主要因素。据样品测试点在Al-Fe-Mn三角图上的投点显示,研究区样品主要落在生物成因区内,仅有一个样品点落在生物成因区外,靠近热液成因区,进一步表明研究区五峰组—龙马溪组页岩硅质成因以生物和陆源输入为主（图6）,但也不能排除热液的影响；另外研究区五峰组—龙马溪组页岩镜下可见海绵骨针、腕足类化石以及放射虫等硅质生物化石（图7）。

图  6  昭通太阳区块五峰组—龙马溪组页岩Al⁃Fe⁃Mn三角图解（底图据文献[32]）

Figure 6.  Al⁃Fe⁃Mn ternary diagram of the Wufeng Formation⁃Longmaxi Formation shale, Taiyang block of Zhaotong (base map is from reference [32])

图  7  昭通太阳区块五峰组—龙马溪组页岩典型硅质生物照片

Figure 7.  Typical siliceous biological photographs of the Wufeng Formation⁃Longmaxi Formation shale, Taiyang block of Zhaotong

Si/Al比值也可用来判断硅质成因。研究区五峰组—龙马溪组页岩Si/Al比值较高,Holdaway et al.[34]提出过量硅是指高于正常碎屑沉积环境下的SiO2含量,通常反映生物成因部分的硅质,其计算公式为：Si过量=Si样品-[（Si/Al）背景×Al样品],（Si/Al）背景采用平均页岩比值3.11。研究区五峰组—龙马溪组页岩Si过量介于0.29%~23.73%,平均值为8.52%。Si-Al相关图显示,五峰组—龙马溪组页岩与Barnett页岩类似,在伊利石Si/Al线上为过量硅部分,反映硅质的生物成因部分[29,35]；而部分样品点落入伊利石Si/Al线下区域（图8）,反映五峰组—龙马溪组页岩硅质的陆源碎屑来源。通过硅质总量减去过量硅得到陆源硅的含量,五峰组—龙马溪组页岩陆源硅含量介于6.94%~43.23%,平均值为20.11%,自下而上陆源硅含量逐渐增多,生物硅含量逐渐减少（图3）,表明在五峰组—龙马溪组沉积时期陆源输入不断增多。焦石坝地区五峰组—龙马溪组下段页岩样品过量硅含量介于19.79%~55.31%,平均值为33.4%[30⁃31],川南长宁地区五峰组—龙马溪组页岩的生物硅含量多介于0.24%~61.23%,五峰组—龙马溪组底部平均值可达47.52%[36⁃37]。而研究区陆源硅总体偏低,推测原因为昭通示范区靠近黔中古陆、康滇古陆,导致陆源硅输入较多,生物硅相对较少。

图  8  昭通太阳区块五峰组—龙马溪组页岩及Barnett页岩Si-Al相关图（底图据文献[35]）

Figure 8.  Correlation diagram of the Wufeng Formation⁃Longmaxi Formation shale and Barnett shale Si⁃Al, Taiyang block of Zhaotong (base map is from reference [35])

陆源碎屑物质的输入对有机质的富集起到关键性作用[38]。通常陆源碎屑的输入会带来丰富的营养物质,使得水体中的藻类大量繁殖,由此促进有机质的富集；当陆源碎屑输入过多时,会对有机质产生稀释作用,从而抑制有机质的富集[38⁃39]。Al、Ti等元素在沉积过程中不易受风化作用和成岩作用的影响,常被用来指示陆源物质输入[39]。Murry et al.[40]提出利用公式：陆源物质质量分数（%）＝ω（Ti样品）/ω（TiPAAS）×100来计算陆源输入量。经计算,研究区五峰组—龙马溪组页岩陆源输入量介于10.47%~53.49%,平均值为35.0%,自下而上表现为先减小再增加的特征（图9）。陆源输入的变化与沉积时的古环境密切相关[1⁃2]。在晚奥陶世末期存在短暂的冰期,导致全球海平面下降,陆源输入增多；早志留世早期冰川发生大规模融化,海平面快速上升,发生海侵,陆源输入迅速减少,随着早志留世中后期海退的发生,水体逐渐变浅,陆源碎屑输入逐渐增多[37⁃39],这与U/Th、V/Cr所指示的沉积环境变化（由缺氧环境转变为富氧环境）相符（图3）,随着陆源碎屑物质不断沉积,海水中的沉积物趋于稳定[41]。

图  9  昭通太阳区块五峰组—龙马溪组页岩陆源输入箱型图

Figure 9.  Histogram of terrestrial source input of the Wufeng Formation⁃Longmaxi Formation shale, Taiyang block of Zhaotong

研究区五峰组—龙马溪组页岩TOC含量在龙一1亚段底部1、2小层达最大,向上逐渐减小。据陆源输入与TOC交会图显示,龙一1亚段底部1、2小层的TOC与陆源输入呈较好的正相关关系,R2分别为0.40、0.57（图10）,表明在1、2小层沉积时期陆源碎屑的输入带来了丰富的营养物质,有利于生物繁盛,进而促进有机质的富集,在五峰组、龙一1亚段3、4小层及龙一2亚段陆源碎屑的输入较多,对有机质产生一定的稀释作用,因此在一定程度上抑制了有机质的富集。

图  10  昭通太阳区块五峰组—龙马溪组页岩TOC、陆源输入相关关系图

Figure 10.  Correlation diagram of total organic carbon (TOC) and terrestrial source input of the Wufeng Formation⁃Longmaxi Formation shale, Taiyang block of Zhaotong

稀土元素在水体中主要赋存在细粒碎屑物质或悬浮物中,这些载体在水体中停留时间的长短决定了REE的分异程度,因此REE分异程度是研究沉积速率的重要手段之一。在沉积过程中轻稀土元素通常会被吸附在黏土矿物或者有机质中优先沉积,而重稀土元素常常在水体中形成稳定络合物,进而导致轻重稀土元素强烈分异[42]。La和Yb分别指示轻稀土元素和重稀土元素。因此,可以利用LaN/YbN值和稀土元素标准化配分模式曲线的斜率来体现沉积速率相对大小[42⁃43]。当沉积物在水体中的停留时间越长,沉积速率越慢时,稀土元素的分解程度越彻底,有机质等物质发生反应的程度越高,REE分异程度越强,稀土元素配分模式曲线斜率越大,LaN/YbN值明显大于1或小于1,即偏离1的程度越大；反之,沉积物与水体中的物质交换少,REE分异程度弱,曲线越平缓,分异程度越小,LaN/YbN值越接近1 [43]。研究区五峰组—龙马溪组LaN/YbN值介于3.72~14.42,平均为9.19,自下而上LaN/YbN值呈由低—高—低—高的变化趋势,反映五峰组—龙马溪组页岩的沉积速率降低—增加—降低的变化过程,与陆源输入呈大致相同的变化趋势（图3）。沉积速率的变化与前文所述的沉积环境变化相符合。五峰组沉积时期陆源输入较多,沉积速率较快；龙马溪组沉积早期陆源输入减少,沉积物在水体中长时间滞留,沉积速率减慢,早志留世中后期陆源碎屑输入不断增多,沉积速率增大,随着沉积物稳定沉降[38⁃39],陆源输入减小,沉积速率也再次减小。研究区五峰组—龙马溪组稀土元素标准化配分模式曲线较为平缓（图5）,整体沉积速率较快,表明其沉积时期距离物源区较近。

由于各元素稳定性的差异,部分稳定性较强微量及稀土元素能够较好地保留源区的地球化学信息,因此常被用来表征物源类型[1,44]。通常物质来源稳定的沉积物稀土元素配分模式较为相似[45⁃46]。例如何佳伟等[2]研究发现盐津牛寨剖面龙马溪组各页岩样品之间的配分模式差异不大,配分模式曲线间近于平行,源岩主要为沉积岩及花岗岩,来源基本一致。对研究区五峰组—龙马溪组页岩样品进行球粒陨石标准化[47],据稀土元素配分模式图显示,五峰组和龙马溪组的稀土元素配分模式基本一致,呈较为平缓的右倾特征,五峰组的配分模式曲线整体波动幅度较小,龙马溪组大多数样品配分模式曲线基本平行,个别样品出现异常,稀土元素含量整体偏高（图5）,说明可能存在混合物源；此外,配分模式图中还显示研究区五峰组—龙马溪组页岩存在Eu负异常和Ce弱负异常,与上地壳的稀土元素配分模式一致,反映研究区五峰组—龙马溪组的源岩来自上地壳[1,44,46]。

通常情况下花岗岩的δEu值多小于0.90,玄武岩、中性斜长石等δEu值大于0.90[48⁃49]。研究区五峰组—龙马溪组δEu介于0.42~0.85,均小于0.9,指示其源岩主要为花岗岩。Roser et al.[50]最早提出了依据主量元素含量计算的砂、泥岩物源区判别函数（表4）,判别函数图（F1-F2）显示,研究区五峰组—龙马溪组页岩样品主要落在石英质沉积岩源区、酸性火山岩源区,少部分样品落在中性火山岩源区（图11a）。微量元素Hf-La/Th图解中[51],样品点主要落在长英质物源以及与基性岩混合物源区域,部分样品落在被动大陆边缘附近（图11b）。随着岩浆不断演化,La、Th等元素不断富集,Sc、Co等元素逐渐亏损,因此酸性火成岩通常显示有更高的La/Sc值和更低的Co/Th值[49,52⁃53]。研究区样品La/Sc比值介于2.57~9.74,Co/Th值介于0.10~1.93,在La/Sc-Co/Th图解中主要分布在长英质火成岩附近[53],少数样品点散落在长英质火成岩和花岗岩中间区域（图11c）,判断其原始物质可能来自长英质火成岩和花岗岩。ΣREE-La/Yb图解常能够有效判别源岩属性[46,54],该图解显示研究区五峰组—龙马溪组样品主要落在沉积岩和花岗岩交汇区域,少数落在碱性玄武岩、花岗岩交汇区（图11d）,反映五峰组—龙马溪组页岩的源岩主要为沉积岩类和花岗岩类。综合上述判断研究区五峰组—龙马溪组可能存在混合物源,其源岩主要来自花岗岩和沉积岩。

表 4  砂、泥岩物源区主量元素判别函数的变量及其系数[50]

Table 4.  Variables and their coefficients of discriminant function for provenance of sandstone and mudstone[50]

图  11  昭通太阳区块五峰组—龙马溪组页岩的物源区岩石属性判别交会图

Figure 11.  Discrimination diagram of the Wufeng Formation⁃Longmaxi Formation shale, Taiyang block of Zhaotong

海底热液对黑色页岩沉积具有重要影响[48,55]。前人常用Ba/Sr、Co/Zn值来识别海底热液活动。通常正常海相沉积物的Ba/Sr值小于1,受海底热液影响的沉积物Ba/Sr值大于1[56]。研究区五峰组—龙马溪组页岩的Ba/Sr值介于1.16~17.50,普遍大于1,表明其沉积时受到海底热液影响。热液成因的沉积物中Co/Zn值约为0.15,铁锰结核等正常沉积形成的岩石中Co/Zn值一般在2.5左右[27]。研究区五峰组—龙马溪组页岩的Co/Zn值介于0.01~0.40,平均为0.12,接近0.15,同样反映存在热液活动。Zn⁃Ni⁃Co三角图也是判断热液影响的方法之一[57⁃58],该图解显示样品主要落在受热液影响的沉积区及其附近,进一步反映热液活动的存在（图12a）。通常可以利用ΣREE/Fe值来反映受热液影响的沉积物距离热液活动中心的远近,该比值越小,距离热液活动中心就越近,反之则越远[59⁃60]。如东太平洋洋隆附近沉积物的ΣREE/Fe值为6.2×10-4,而距离洋隆802 km之外的沉积区该比值约为28.8×10-4[59]。研究区五峰组页岩的ΣREE/Fe值介于57.32×10-4~141.91×10-4,平均为96.77×10-4,龙马溪组页岩的ΣREE/Fe值介于14.90×10-4~563.30×10-4,平均为73.85×10-4,说明研究区五峰组—龙马溪组页岩沉积时离热液活动中心较远。Alexander et al.[61]提出Eu/Sm-Sm/Yb二元图解用于判断海水与海底热液的贡献量。研究区五峰组—龙马溪组页岩接近于水成铁锰质地壳与海水（图12b）,表明在原始溶液中海底热液的占比尚不足0.1%。此外海底热液常表现出强烈的Eu正异常,学者们利用Eu正异常程度判别高温热液的贡献量[62],当样品的δEu值大于1,并且数值越大,说明样品正异常程度越高,高温热液的贡献量也就越大。研究区样品均表现为Eu负异常,进一步说明原始溶液中海底热液组分较少,对页岩沉积影响极小,推测其为远离热液活动中心的远端沉积[27]。前文的Al-Fe-Mn三角图解（图6）主要是针对硅质成因进行判别,该图解没有样品点落入热水成因区中,并不表示无热液影响。通过以上研究证实研究区五峰组—龙马溪组页岩的硅质来源受到了热液的影响。但由于热液组分含量极低,从而对于Fe、Mn元素的响应不大,导致无样品点落入热水成因区,热液对于硅质成分的贡献极小,因此这与上文黑色页岩受热液影响的结论并不相悖[29]。

图  12  昭通太阳区块五峰组—龙马溪组页岩物源判别图解

Figure 12.  Provenance discrimination diagram of the Wufeng Formation⁃Longmaxi Formation shale, Taiyang block of Zhaotong

伴随着沉积作用及构造活动的发生,不同地质背景下的陆源碎屑岩中地球化学特征会存在差异,部分微量元素和稀土元素在沉积过程中比较稳定,其元素特征和比值范围能够较好地保留物源区的相关信息,常被用来研究物源区母岩属性及构造背景[46,48⁃49,63]。Bhatia et al.[63]研究发现不同构造背景下的杂砂岩的REE分布存在一定差异,并总结出规律作为不同构造背景的识别依据（表5）,通过对研究区La、Ce、∑REE及La/Yb平均值等数值进行计算,再与前人研究总结出的各背景下的范围值进行比较,依据对比结果综合判断出最为接近的构造背景。此外利用主量元素判别不同构造背景已得到广泛应用,常用的主量元素指标有SiO2、K2O、Na2O和Al2O3等[24,63]。

表 5  昭通太阳区块五峰组—龙马溪组页岩样品与不同构造背景沉积盆地杂砂岩的REE特征对比[63]

Table 5.  Comparison of rare earth element characteristics of mixed sandstones from the Wufeng Formation⁃Longmaxi Formation shale samples from Taiyang block of Zhaotong with heterogeneous sandstones in sedimentary basins with different tectonic backgrounds[63]

利用Sc/Cr-La/Y图解对研究区五峰组—龙马溪组样品数据进行投点,结果显示研究区样品投点较为集中,大部分样品落在被动大陆边缘内,少数样品落在活动大陆边缘和大陆岛弧内及其附近（图13a）。利用前人建立的SiO2/Al2O3-K2O/Na2O双变量判别图版进行投点[63],,所有样品集中落在被动大陆边缘范围内（图13b）,说明研究区构造背景为被动大陆边缘。通常处于被动大陆边缘的沉积物稀土元素具有轻稀土富集,Eu负异常的特征,而活动大陆边缘的沉积物多以重稀土元素富集、无Eu亏损为特征,且多为分异程度较低的火成岩[43,48,64]。研究区样品的LREE/HREE平均值为9.11,δEu平均值为0.69,指示构造背景为被动大陆边缘。将研究区样品稀土元素含量与各背景的范围值对比[63],结果表明研究区La、Ce、∑REE及La/Yb平均值分别为44.79 µg/g、90.68 µg/g、225.98 µg/g、14.73（表5）,均与被动大陆边缘接近。综合以上研究结果,判断研究区五峰组—龙马溪组页岩构造背景主要为被动大陆边缘。

图  13  昭通太阳区块五峰组—龙马溪组构造背景判别图解

Figure 13.  Tectonic setting discrimination diagrams of the Wufeng Formation⁃Longmaxi Formation, Taiyang block of Zhaotong

研究区主要位于上扬子地块东南部,而扬子地块是前寒武纪华南板块的重要组成部分,在早古生代与华北板块、塔里木板块发生分离[22,29]。寒武纪末期,华夏板块与扬子地块受广西造山运动的影响开始发生汇聚,扬子板块东南部不断抬升,扬子地区进入被动大陆边缘演化阶段。扬子地区奥陶世先以碳酸盐台地沉积为主,自晚奥陶世起,以碳酸盐岩—碎屑岩混合沉积为主。在奥陶纪—志留纪之交发生都匀运动,华南板块向北持续俯冲,与华北板块、滇缅板块等发生碰撞,使得扬子板块内部发生挤压变形,川中隆起、黔中—滇东隆起、雪峰隆起等开始形成,由此控制了晚古生代的构造演化,扬子西缘整体仍在持续抬升,因此仍处于被动大陆边缘环境,与上述研究结论一致[22,64]。奥陶纪末期冈瓦纳大陆南极冰盖开始扩张,扬子地台发生海退,早志留世早期,随着华夏板块的不断扩张以及华南板块大部分地区的持续隆升,康滇古陆、黔中隆起以及雪峰隆起等古陆及水下高地开始形成,并且该时期伴随着冰川的大规模融化,由于各个古隆起的阻隔,上扬子地区处于半封闭的沉积环境,广泛沉积一套黑色页岩,即龙马溪组页岩[3,8,22]。

研究区西侧的康滇古陆位于扬子地台西南缘。康滇古陆自显生宙以来长期暴露地表[65],前寒武纪以前,康滇古陆逐渐开始形成岛弧褶皱带,晋宁运动导致岛弧褶皱发生转变,转变成早震旦世的安第斯型山弧,并与扬子大陆统一[66⁃67]。晚震旦世康滇构造带基本处于稳定状态,中奥陶世—二叠纪,地幔物质上涌,上地壳发生上隆,康滇古陆区域也大面积隆升。据前人研究发现,康滇古陆部分区域（如泸定—石棉一线）主要出露新元古代地台盖层,岩性主要为一套中—酸性火成岩—火山碎屑岩,在其内部发现有大量同源花岗岩侵入体[2,29,67],与前文研究得出的花岗岩源岩属性结论一致。结合五峰组—龙马溪组页岩沉积时构造背景,在晚奥陶世—早志留世黔中隆起紧邻研究区开始发育,与龙马溪组开始沉积时间基本同步[2,68],据此综合推断研究区五峰组—龙马溪组页岩的源岩可能来自康滇古陆新元古地台盖层和黔中隆起。