Enrichment characteristics of strategic metal elements in Late Paleozoic coal-bearing sequences in Renjiazhuang Mine Field
-
摘要:
煤系是战略性金属矿产增储的重要来源。宁东煤田煤炭资源丰富，部分地区煤中富集战略性金属。以宁东煤田任家庄井田煤系为研究对象，系统采集主煤层剖面煤样、夹矸及顶(底)板样，采用电感耦合等离子体质谱法(ICP-MS)、X-射线荧光光谱(XRF)和X射线衍射(XRD)进行了系统的矿物学和地球化学测试分析，查明常量元素和微量元素分布及富集情况，讨论战略性金属元素富集成因并进行物源分析。结果表明，任家庄井田石炭－二叠纪煤系富集战略性金属矿产，太原组和山西组煤中战略性金属元素整体含量水平较高，以Li-Ga-Zr(Hf)-Nb(Ta)-Th(U)-Pb及REY共富集为特征，且顶底板及夹矸更为富集。其中：5号煤层Li、Ba高度富集，Rb、Zr富集；9号煤层Li高度富集，Zr、Hf、Th富集；9上煤层微量元素含量整体相对较低，Cr富集。5号煤层和9号煤层战略性金属元素富集程度较9上煤层更好。5号煤层和9号煤层灰基Li、Ga均达到了开发利用品位，5号煤层灰基REY也达到了开发利用品位，具有良好开发利用前景。基于元素赋存特征、沉积环境、构造演化的综合分析，认为任家庄井田战略性金属元素富集主要受陆源碎屑供给影响，由鄂尔多斯盆地西缘北部阴山古陆与西北部阿拉善地块蚀源区的长英质−中性岩碎屑物质经古水流搬运进入泥炭沼泽，在特定沉积环境及地球化学环境下通过有机质与无机质相互作用而富集，并经多期次构造演化控制形成现今赋存分布状态及富集特征。
Abstract:Coal-bearing sequences are important source for increasing strategic metal mineral reserves. There are abundant coal resources in Ningdong Coalfield,and strategic metals enrichment have been found in several coal-bearing sequences. The study selected coal-bearing sequences of Renjiazhuang Mine Field in Ningdong Coal Field as the research object, the coal samples, gangue, roof and floor plate of the main coal seam section were collected systematically. The mineralogical and geochemical tests were carried out by inductively coupled plasma mass spectrometry (ICP-MS), X-ray fluorescence spectroscopy and X-ray diffraction, and the distribution characteristics of major and trace elements were ascertained. The causes of strategic metal element enrichment and their provenance were analyzed.The results showed that the Permo-Carboniferous coal-bearing of late Paleozoic are the enrichment beds of strategic metal minerals in RenJiazhuang Mine Field. The enrichment of strategic metal elements in coal is common in Taiyuan and Shanxi formations, and the overall content of trace elements is high, which is characterized by Li-Ga-Zr(Hf)-Nb(Ta)-Th(U)-Pb,rare earth element and Yttrium(REY), co-enrichment. Moreover, strategic metal elements are more enriched in coal-bearing near the roof,floor and the gangue.Among them, Li and Ba are highly enriched, Rb and Zr are enriched in No.5 coal seam. Li is highly enriched in No.9 coal seam, while Zr, Hf and Th are enriched. The content of strategic metal elements in upper No.9 coal seam is relatively low and Cr is enriched. The concentration of trace elements in No.5 and No.9 coal is better than that in upper No.9 coal seam. Li and Ga of the ash foundation of No. 5 and No. 9 coal seam have reached the industrial grade, and REY of the ash foundation of No. 5 coal seam has also reached the industrial grade, which presenting a very promising potential for further exploration.Based on the comprehensive analysis of element occurrence characteristics, sedimentary environment and tectonic evolution, it is concluded that the strategic metal element enrichment in Renjiazhuang Mine Field is mainly influenced by the clastic supplying of the erosion source area, and the felsitic-neutral rock detrites from the northern Yinshan ancient landthe and the northwestern Alxa block in the western margin of Ordos Basin were transported into the peat swamp by palaeocurrent. It was enriched through interaction between organic matter and inorganic matter in specific sedimentary environment and geochemical environment, and formed the present distribution state and enrichment characteristics under the control of multi-stage tectonic evolution.
-
Keywords:
- Ordos Basin /
- Ningdong Coalfield /
- strategic metal element /
- permo-carboniferous /
- coal-bearing strata /
- provenance
-
0. 引 言
国家经济的飞速发展与能源的大量消耗密不可分。我国“相对富煤、少油、缺气”的资源特征，决定了今后相当长一段时间煤炭资源仍将成为我国重要的能源支持[1]。煤岩是一种复杂的多孔介质，主要是由基质和裂隙组成[2]。煤岩的孔隙结构决定煤层气吸附、解吸规律，而裂隙是煤层气的有效运移通道，决定着煤层气的运移和产出[3]。然而，煤储层裂隙微观结构复杂，往往表现出非均质性，严重影响煤储层裂隙的渗流特性。与此同时，煤基质将产生膨胀、收缩变形，导致了裂隙变形，而煤岩渗透特性受煤岩裂隙表面粗糙度、应力压缩及滑脱效应等影响，会有不同程度的变化[4]。因此，研究在多种因素综合作用下煤层气在储层的运移规律，对提高煤层气采收率有重要的指导意义与参考价值。
将物性试验与分形特征结合是目前对煤岩孔裂隙结构进行定量分析的一种重要方法[5]。XIE等[6]在对岩石力学问题的研究中将分形维数与岩石力学相关理论相结合，提出分形－岩石力学理论，完善了之前采用常规方法进行研究的缺点；与此同时，王世芳等[7]基于分形理论和达西定律，提出了牛顿流体在粗糙裂缝型多孔介质中渗透率的分形模型，揭示了影响粗糙裂缝型多孔介质渗透率的物理机理；基于此，赵明凯等[8]通过分形理论及流体力学基本原理定量表征煤岩裂隙表面粗糙微元及裂隙开度的分布规律，提出不可忽视裂隙粗糙度对渗流的影响；在粗糙微通道的孔隙形态和内壁粗糙微元分形研究基础上，王登科等[9]推导出粗糙微纳米孔隙中的瓦斯气体传输方程。
理论上，气测和液测煤岩绝对渗透率一致，然而实际分析测试过程中普遍存在着气测煤岩和液测煤岩绝对渗透率不同的问题，其主要原因是气体存在滑脱效应，且随着煤层绝对渗透率降低，这种气体滑脱效应变得更加显著[10]。ZHU等[11]采用含瓦斯煤热‒流‒固耦合三轴伺服渗流装置对煤样应变及气相渗透率进行了观测，建议在对渗透率与有效应力特性进行研究时充分考虑气体滑脱效应；ZHOU等[12]分析了克林肯伯格系数与裂隙宽度之间的关系，并建立了具有有效应力和气体吸附效应的裂隙宽度理论模型；魏建平等[13]通过实验分析与理论推导，建立了考虑气体动力黏度，压缩因子与滑脱效应的煤层瓦斯渗透率方法。因此，在煤层气的抽采过程中，煤岩裂隙表面粗糙度和滑脱效应都会对煤岩的渗透率产生一定的影响。
综上所述，相关学者分别基于气体压力及其引起的吸附变形对煤岩渗透特性演化机制进行了诸多研究，但考虑煤岩裂隙表面相对粗糙度与滑脱效应等多因素共同影响的煤岩表观渗透率规律研究鲜有报道。鉴于此，笔者选取贵州六盘水矿区煤样，进行扫描电镜和氮气吸附试验，分析了几个矿区的煤岩孔裂隙特征及渗流规律。并基于分形理论定量表征煤岩裂隙表面相对粗糙度，建立了考虑裂隙表面相对粗糙度、吸附变形及滑脱效应的表观渗透率模型。基于渗透率测试结果，通过匹配试验结果与理论结果，进一步探讨并量化分析煤岩裂隙表面粗糙程度及滑脱系数与渗透率的响应机理。研究实际储层中的煤岩渗流规律，旨在为煤层气抽采和煤层增透过程中提供一定的理论支撑。
1. 试验样品采集与制备
本试验煤样分别取自贵州六盘水的金佳煤矿、马依煤矿、东李煤矿17号、湘桥煤矿、打牛厂煤矿、月亮田煤矿。对现场采集的新鲜煤样进行封存，并运送至试验室进行进一步研究，取用0.18～0.25 mm（60～80目）的干燥煤样依据标准GB/T212—2008《煤的工业分析方法》对煤样的水分（Mad）、灰分（Ad）、挥发分（Vdaf）的成分进行测定，工业分析测试结果见表1。
表 1 煤样工业分析测试结果Table 1. Industrial analysis of coal samples煤样编号 采样位置 Mad/% Ad/% Vdaf/% 煤阶 JJ 金佳煤矿 2.26 27.50 11.87 低挥发分烟煤 MY 马依煤矿 1.35 17.44 12.29 低挥发分烟煤 DL 东李煤矿17号 1.27 3.79 20.64 中挥发分烟煤 XQ 湘桥煤矿 1.05 6.97 29.44 中高挥发分烟煤 DNC 打牛厂煤矿 2.44 9.06 30.88 中高挥发分烟煤 YLT 月亮田煤矿 1.32 15.97 32.98 中高挥发分烟煤 由于煤的挥发分反映了煤的变质情况，煤样按照由高至低的变质程度进行排序。由表1可知，JJ煤样变质程度最高，MY次之，YLT煤矿变质程度最低，均为变质程度较高的烟煤。
2. 煤岩孔裂隙结构特征
煤岩是由基质和裂隙组成的双孔隙、非均质性较强的多孔介质。为探究煤岩非均质性对煤岩渗流的影响，本节先通过扫描电镜试验探究裂隙表面形 态特征，再利用低温氮气吸附试验结果结合分形模 型和Hurst理论，观测表征煤岩的裂隙粗糙程度和孔 裂隙分布特征。
2.1 煤岩孔裂隙结构特征
2.1.1 扫描电镜试验
将煤样制成长×宽×高为1 cm×1 cm×1 cm的煤块，为保证效果，将煤样进行喷金处理。采用德国ZEISS Gemini SEM 300型自带能谱扫描电镜，保持煤样表面水平，进行扫描电镜（SEM）试验，得到煤岩表面平整度和裂隙结构在不同放大倍数时的成像结果，并在此基础上对煤岩微观裂隙结构的分布规律进行了分析。为探究煤岩裂隙表面特征，笔者选取煤样放大500倍和2 000倍的SEM图像进行分析（如图1所示）。由图1可知，SEM图像直观展示出不同地区煤岩裂隙形态特征，JJ、MY、DL、XQ、DNC、YLT煤样裂隙均较为发育，裂缝宽度在1～10 μm，裂隙表面非均质性较强。观察不同地区煤岩裂隙表面形态特征可知，JJ煤样（如图1a）表面凹凸不平，极大地增加了表面积，裂隙比较发育。MY煤样（如图1b）表面有一定量的煤屑存在，矿物颗粒之间充填比较完全，胶结程度较好，裂隙形状较发育。DL煤样（如图1c）表面有少量碎屑分布，在不同方向上裂隙发育良好。XQ煤样（如图1d）表面分散有少量碎屑，裂隙较宽，明显看到其裂隙表面不光滑。DNC煤样（如图1e）表面不平坦，裂隙发育较宽。YLT煤样（如图1f）表面孔裂隙发育，裂隙宽度是最宽的，明显看出其裂隙表面粗糙。通过对SEM的煤岩裂隙形态特征的观察，裂隙表面粗糙度也许会影响煤层气流动速率，进而导致煤岩渗透率降低。
![]()
图 1 基于SEM的煤岩裂隙形态特征Figure 1. Morphological characteristics of coal fractures based on SEM2.1.2 低温氮气吸附试验
前一小节通过SEM试验发现煤岩表面具有明显的非均质性，于是该小节利用低温氮气吸附试验分析煤岩孔隙微观结构特征。低温氮气吸附试验采用美国Micromeritics公司ASAP2020型比表面微孔分析仪进行分析。利用低温氮气吸附试验分析煤岩孔隙微观结构特征。煤岩低温氮气吸附试验结果如图2和表2所示，当压力较大时，低温氮气试验中会发生毛细凝聚，造成凝聚与蒸发过程中的相对压力不同而形成滞后环。滞后环形状不一，反映了孔的类型及孔径分布。根据最新的IUPAC分类方案[14]，将滞后环分为H1、H2、H3、H4和H5这5种类型。如图2所示，试验表明6种类型的孔均产生脱附回线。从滞后环的形状来看，煤样的吸附‒脱附等温线均为H2(b)&H3复合型，MY、DL、DNC煤样的滞后环较小，表明其孔隙率较低，存在楔形半封闭孔隙。
![]()
图 2 煤样低温氮气吸附试验等温线Figure 2. Low temperature nitrogen adsorption isotherm of coal samples表 2 低温氮气吸附法测定孔隙结构参数Table 2. Analysis results of low- temperature nitrogen adsorption isotherms煤样 BET比表面积/
(m2·g‒1)总孔容/
(cm3·g‒1)平均孔径/
nmJJ 1.374 0.003 001 8.737 MY 1.205 0.004 770 15.830 DL 0.786 0.003 082 15.682 XQ 1.357 0.005 199 15.328 DNC 2.266 0.008 404 14.618 YLT 0.970 0.004 445 18.322 2.2 煤岩分形特征
分形维数能较好地反映煤岩孔裂隙的发育程度，为了探讨分形维数与煤岩孔裂隙特征及渗流特性之间的联系，根据低温液氮试验资料，采用FHH模型对贵州六盘水矿区煤岩孔裂隙分形维数进行计算。煤岩孔隙表面分形维数介于2～3，愈靠近2说明煤岩孔隙表面愈平滑，靠近3说明煤岩孔隙表面粗糙度愈大。FHH[15]模型计算表达式为
$$ \ln \left( {\frac{V}{{{V_{\text{m}}}}}} \right) = A\ln \left[ {\ln \left( {\frac{{{p_{\text{0}}}}}{p}} \right)} \right] + C $$ (1) 式中：V为平衡压力为p时吸附的气体分子体积，cm3/g；Vm为单分子层饱和吸附量，cm3/g；A为双对数曲线的斜率，与煤的孔隙分形维数D及吸附过程相关；p0为液氮温度时氮的饱和蒸气压，MPa；p为氮气分压，MPa；C为与气体吸附和凝结有关的常数；基于分形FHH模型，前人认为lnV和ln[ln(p0/p)]关系式的斜率A与D存在明确关系，并根据D=A+3计算得到煤岩分形维数，其计算结果如图3和表3所示。
表 3 基于FHH模型的孔隙分形维数Table 3. Pore fractal dimension based on FHH model煤样编号 A R2 D JJ ‒0.378 0.997 2.622 ‒0.645 0.993 2.355 MY ‒0.205 0.964 2.795 ‒0.858 0.991 2.142 DL ‒0.338 0.993 2.664 ‒0.747 0.985 2.253 XQ ‒0.384 0.998 2.616 ‒0.673 0.993 2.327 DNC ‒0.294 0.967 2.706 ‒0.779 0.992 2.221 YLT ‒0.300 0.998 2.700 ‒0.940 0.986 2.060 研究表明，D表征的为孔隙表面的非均质性特征。如图3所示，低压段分形维数越大，则孔隙内部表面的形状越不规则，光滑度越差；而高压段分形维数主要反映孔隙结构的特征，包括孔径分布与大小，分形维数越大，则孔径分布越分散，孔径越小[16]。6个煤样均为高压区分形维数（p>0.5）<低压区分形维数（p<0.5）。这表明，煤样的中大孔隙（d>2 nm）与微孔孔隙（d<2 nm）相比，孔隙结构更复杂。
2.3 煤岩粗糙度的定量描述与表征
近年来，多重分形理论在各领域都得到了国内外学者的广泛运用。利用低温氮气吸附试验对不同变质程度煤样进行孔隙结构的测定，而多重分形奇异谱α～f(α)和广义维数谱q～Dq这2种算法等价于描述多重分形特征[17]，可用于对煤样的氮气吸附分布数据进行分析。图4所示是煤样广义分形维数谱q～Dq，其中Dq为与q相关的广义分形维数，q表示矩阵的阶数（‒∞<q<+∞）。在本研究中，q的范围为−10～10，对应特征参数见表4。D1、D2和D3分别表示容量维度、信息维度和相关维度。研究发现，当q=2时，表示如下[18]
![]()
图 4 煤样孔径分布广义分形维数谱曲线Figure 4. Generalized fractal dimension spectrum plot of coal sample pore size distribution表 4 煤样孔径分布广义分形维数谱参数Table 4. Parameters of generalized fractal dimension spectrum of pore size distribution of coal samples样品编号 D‒10 D0 D1 D2 H D10 D‒10‒D10 D‒10‒D0 D0‒D10 JJ 1.482 1 0.740 0.536 0.768 0.329 1.474 0.689 0.787 MY 1.687 1 0.587 0.455 0.681 0.213 1.152 0.482 0.671 DL 1.774 1 0.558 0.333 0.667 0.195 1.724 0.919 0.805 XQ 1.600 1 0.552 0.324 0.662 0.189 1.517 0.706 0.811 DNC 1.919 1 0.669 0.452 0.723 0.272 1.328 0.600 0.728 YLT 1.706 1 0.646 0.414 0.707 0.245 1.290 0.535 0.755 $$ H = {{\left( {{D_2} + 1} \right)} \mathord{\left/ {\vphantom {{\left( {{D_2} + 1} \right)} 2}} \right. } 2} $$ (2) 式中：H为Hurst指数（0.5<H<1），H可用于表征粗糙度，指定空间孔隙的连通性。
全部样品的D0值为1，说明每个盒子中都存在孔隙，这与一维分布的欧式维数相对应。q～Dq谱能够有效地表征各阶段孔径分布特征的非均质性与复杂程度，从而揭示局部差异特征。如图4所示，全部煤样谱线都呈现反“S”型单调递减的趋势，说明孔径分布呈现多重分形的特点，这与μ(q,ε)函数的分析结果一致。就煤样谱宽变化而言，其与煤样变质程度并无明显联系。Hurst指数可用于反映孔隙连通性情况，由表4可见，JJ煤样Hurst指数最大，孔隙间联通性良好。
3. 裂隙粗糙度对煤岩渗透率的影响
3.1 建立考虑煤岩裂隙表面粗糙度和气体滑脱效应的渗透率模型
在第2章节中，通过扫描电镜试验观测到孔裂隙结构表面存在粗糙微元，同时将低温氮气吸附试验与FHH模型结合，定量表征煤岩非均质性。研究表明煤岩孔裂隙结构发育，存在分形特征且非均质性较强，并利用Hurst指数定量表征煤岩粗糙度表明其对煤岩渗透率的重要性。然而不少研究学者发现煤岩渗流过程中会发生滑脱效应，且其影响不可忽略，因此，笔者利用煤岩裂隙表面粗糙度修正滑脱系数，定量表征滑脱系数对煤岩渗透率的影响，并建立煤岩表观渗透率模型。
3.1.1 考虑煤岩裂隙表面粗糙度、应力及气体吸附的煤岩裂隙开度模型
基于扫描电镜试验分析结果，假设裂隙面两侧的粗糙程度相同，为简化计算，将裂隙通道简化为平板，粗糙微元假设为半球体结构（如图5所示）。并且根据分形标度定律，裂隙与粗糙微元由结构参数分形表征。
基于假设条件、粗糙裂隙平行板简化模型及王登科等[9]对粗糙微通道的分形表征研究基础，建立考虑煤岩裂隙表面粗糙度的相对粗糙度模型。裂隙表面半球体粗糙微元的有效平均高度可表示为
$$ \overline{h_{\mathrm{e}}}=\frac{V}{S\mathrm{_t}}=\frac{\varphi a_{\max}}{3}\frac{2-D_{\text{a}}}{3-D_{\text{a}}}\frac{1-\varepsilon^{3-D_{\text{a}}}}{1-\varphi} $$ (3) 式中：$ \overline{h_{\mathrm{e}}} $为裂隙半球体粗糙微元的等效平均高度；V为裂隙的单侧表面粗糙微元的体积；St为粗糙裂隙单侧表面的总面积；φ为在1个分形集合中，全部半球体底面总面积与分形集合单元总面积之比；amax为粗糙微元底边的最大直径；Da为粗糙微元底边直径分布分形维数，二维情况下，0≤Da<2；ε为将裂隙面粗糙微元的底边最小直径amin和最大直径amax之比，ε=amin/amax。
煤岩裂隙表面相对粗糙度ξ为粗糙微元的等效高度与裂隙开度b的比值，即
$$ \xi = \frac{{2\overline {{h_{\text{e}}}} }}{b} $$ (4) 赵明凯等[8]研究发现研究域内裂隙面的粗糙度几乎相同。由于裂隙半球体粗糙微元的高度h为粗糙微元底边直径a的一半，即h=a/2，结合式（3）和式（4），可以得到裂隙开度与粗糙微元底面直径的关系
$$ \frac{{{b_{\min }}}}{b} = \frac{{{{\overline {{h_e}} }_{{\text{min}}}}}}{{\overline {{h_e}} }} = \frac{{{{\left( {{a_{\max }}} \right)}_{\min }}}}{{{a_{\max }}}} $$ (5) 式中：bmin为裂隙开度的最小值；$ \overline h_{{\mathrm{e}}\text{min}} $为裂隙半球体粗糙微元等效平均高度的最小值；(amax)min为开度最小的裂隙的粗糙元底边最大长度。
将式（3）和式（5）代入式（4），由此可得煤岩裂隙表面相对粗糙度的表达式
$$ \xi = \frac{{{\text{4}}\varphi {{\left( {{h_{{\text{max}}}}} \right)}_{\min }}}}{{3{b_{\min }}}}\frac{{2 - {D_{\text{a}}}}}{{3 - {D_{\text{a}}}}}\frac{{1 - {\varepsilon ^{3 - {D_{_{\text{a}}}}}}}}{{1 - \varphi }} $$ (6) 式中：(hmax)min为开度最小的裂隙的粗糙微元最大高度。
笔者将裂隙开度b与粗糙微元等效高度$ \overline {{h_{\text{e}}}} $的差值作为流体在裂隙通道内流动的初始裂隙开度b0。仅考虑煤岩裂隙表面相对粗糙度的初始裂隙开度可表式为
$$ {b_{\text{0}}} = b - 2\overline {{h_e}} = b\left( {1 - \xi } \right) $$ (7) 结合式（7）和ZHOU等[12]对渗透率模型研究，有效应力、裂隙表面相对粗糙度及气体吸附效应综合作用下的有效裂隙开度为
$$ {b_{{\text{eff}}}} = \left. {\left\{ \begin{gathered} {b_{\text{0}}} - \frac{{{b_{\text{0}}}\alpha }}{{{\phi _{\text{0}}}K}}\left[ {\left( {\sigma - {\sigma _{\text{0}}}} \right) - \left( {P - {P_{\text{0}}}} \right)} \right] \\ - {b_{\text{0}}}f\left( {\frac{{{S_{{\text{max}}}}P}}{{{P_{\text{L}}} + P}} - \frac{{{S_{{\text{max}}}}{P_{\text{0}}}}}{{{P_{\text{L}}} + {P_{\text{0}}}}}} \right) \\ \end{gathered} \right.} \right\} $$ (8) 式中：beff为考虑应力和裂隙表面相对粗糙度及气体吸附效应综合作用下的有效裂隙开度；α为Biot的系数；σ和P分别为应力和孔隙压力；σ0和P0分别为初始平均应力和初始孔隙压力；ϕ0为煤样初始孔隙度；K为煤岩体积模量；f为内膨胀系数，表示微裂隙因基质吸附膨胀而产生的变形与基质膨胀变形的比值；Smax为吸附引起的最大应变；PL为测量应变为0.5基质块的吸附诱导应变时的孔隙压力。由于Kf=ϕ×K/α，Kf为裂隙的体积模量，ϕ为煤样孔隙度；Cf=1/Kf，Cf为裂隙压缩系数，即Cf=α/ϕ×K。将其代入式（8），可得
$$ {b_{{\text{eff}}}} = \left\{ \begin{gathered} {b_{\text{0}}} - {b_{\text{0}}}{C_{\text{f}}}\left[ {\left( {\sigma - {\sigma _{\text{0}}}} \right) - \left( {P - {P_{\text{0}}}} \right)} \right] \\ - {b_{\text{0}}}f\left( {\frac{{{S_{{\text{max}}}}P}}{{{P_{\text{L}}} + P}} - \frac{{{S_{{\text{max}}}}{P_{\text{0}}}}}{{{P_{\text{L}}} + {P_{\text{0}}}}}} \right) \\ \end{gathered} \right\} $$ (9) 3.1.2 考虑滑脱效应的煤岩表观渗透率模型
假定裂隙内流体流动形态是层流，根据立方定律，可以计算单条微裂缝中气体初始流量q为
$$ q = \frac{{lb_{^0}^3}}{{12\mu }}\frac{{\Delta p}}{{{L_{\text{t}}}}} $$ (10) 式中：l为裂隙长度；Δp为气体压力差；μ为流体黏度；煤岩中裂隙发育，具有分形特征，当考虑裂隙弯曲度时，煤岩裂隙实际渗流通道长度Lt可表示为$ {L_{\text{t}}}\left( b \right) = {b^{1 - {D_{\text{T}}}}}L_0^{{D_{\text{T}}}} $，L0为气体沿流动方向直线长度，其中Lt(b)≥L0；DT为裂隙开度分布的迂曲度分形维数，二维情况下，1<DT<2，三维为1<DT<3。通过煤岩横截面积A0的总流量为
$$ Q = - \int_{{b_{{\text{min}}}}}^{{b_{{\text{max}}}}} {q(b){\mathrm{d}}N(b)} $$ (11) 式中：bmin/bmax<1。由于粗糙微元底边长度在无限小范围内的元素总数为：$ - {\text{d}}N = {D_{\text{f}}}b_{{\text{max}}}^{{D_{\text{f}}}}{b^{ - \left( {{D_{\text{f}}} + 1} \right)}}{\text{d}}b $，将其和ROBERTSON等[19]推导出的煤的孔隙度模型ϕ=3b/l代入式（11），可得
$$ Q = \frac{{\Delta p{D_{\text{f}}}}}{{4\mu \phi L_{\text{0}}^{{D_{\text{T}}}}}}\frac{{{{\left( {1 - \xi } \right)}^3}b_{\max }^{3 + {D_{\text{T}}}}}}{{\left( {3 + {D_T} - {D_{\text{f}}}} \right)}} $$ (12) 根据分形裂隙模型，横截面上的流体通道面积之和为
$$ A = \frac{{3{D_{\text{f}}}}}{{2 - {D_{\text{f}}}}}\frac{{1 - \phi }}{{{\phi ^2}}}\left( {1 - \xi } \right)b_{\max }^2 $$ (13) 根据达西定律，并结合式（12）和式（13），可得粗糙裂隙的初始渗透率表达式为
$$ {k_0} = \frac{{Q\mu {L_0}}}{{A\Delta P}} = \frac{{\phi L_0^{1 - {D_{\text{T}}}}\left( {2 - {D_{\text{f}}}} \right){{\left( {1 - \xi } \right)}^2}b_{\max }^{1 + {D_{\text{T}}}}}}{{12\left( {3 + {D_{\text{T}}} - {D_{\text{f}}}} \right)\left( {1 - \phi } \right)}} $$ (14) 裂隙宽度随煤岩的膨胀或收缩而改变，裂隙宽度对煤岩渗透率有重要影响。因此，绝对渗透率可推导为[12]
$$ {k_\infty } = {k_0}{\left( {{{{b_{{\text{eff}}}}} \mathord{\left/ {\vphantom {{{b_{{\text{eff}}}}} {{b_0}}}} \right. } {{b_0}}}} \right)^3} $$ (15) 将式（9）和式（14）代入式（15）中，得到绝对渗透率为
$$ \begin{gathered} {k_\infty } = \frac{{L_0^{1 - {D_{\text{T}}}}{{\left( {1 - \xi } \right)}^2}b_{\max }^{1 + {D_{\text{T}}}}\left( {2 - {D_{\text{f}}}} \right)}}{{12\left( {3 + {D_{\text{T}}} - {D_{\text{f}}}} \right)}}\frac{\phi }{{1 - \phi }} \times \\ {\left\{ \begin{gathered} 1 - {C_{\text{f}}}\left[ {\left( {\sigma - {\sigma _{\text{0}}}} \right) - \left( {P - {P_0}} \right)} \right] \\ - \frac{{3f}}{{{\phi _0}}}\left( {\frac{{{S_{\max }}P}}{{{P_L} + P}} - \frac{{{S_{\max }}{P_0}}}{{{P_L} + {P_0}}}} \right) \\ \end{gathered} \right\}^3} \end{gathered} $$ (16) 克林伯格方程考虑到气体滑移效应能进一步有效预测致密多孔介质中气体的表观渗透率[20]，并将渗透率与压力的联系起来，克林伯格方程表示为
$$ {k_{{\text{app}}}} = {k_\infty }\left( {1 + {{{B_{{\text{eff}}}}} \mathord{\left/ {\vphantom {{{B_{{\text{eff}}}}} P}} \right. } P}} \right) $$ (17) 式中：Kapp为表观渗透率；Beff为考虑有效应力、裂隙表面相对粗糙度和气体吸附的综合作用下的滑脱系数，可用$ {B_{{\text{eff}}}} = {{\left( {16 {\mathrm{c}}\mu \sqrt {{{2RT} \mathord{\left/ {\vphantom {{2RT} {\left( {\pi M} \right)}}} \right. } {\left( {\pi M} \right)}}} } \right)} \mathord{\left/ {\vphantom {{\left( {16c\mu \sqrt {{{2RT} \mathord{\left/ {\vphantom {{2RT} {\left( {\pi M} \right)}}} \right. } {\left( {\pi M} \right)}}} } \right)} {{b_{{\text{eff}}}}}}} \right. } {{b_{{\text{eff}}}}}} $来确定[21]；c为常数（通常取0.9)，M为气体分子量；R和T分别为通用气体常数和绝对温度；Beff受气体本身性质的影响，与裂隙开度密切相关。于是，考虑应力、裂隙表面相对粗糙度和气体吸附的综合作用下的滑脱系数与裂隙开度的关系式为
$$ {{{B_{{\text{eff}}}}} \mathord{\left/ {\vphantom {{{B_{{\text{eff}}}}} {{B_{\text{0}}}}}} \right. } {{B_{\text{0}}}}} = {{{b_{\text{0}}}} \mathord{\left/ {\vphantom {{{b_{\text{0}}}} {{b_{{\text{eff}}}}}}} \right. } {{b_{{\text{eff}}}}}} $$ (18) 式中：B0为初始状态（不受应力和气体吸附效应综合作用）的滑脱系数。
将式（16）、式（18）代入式（17）可得，考虑有效应力、裂隙表面相对粗糙度、气体吸附及滑脱效应综合作用下的煤岩表观渗透率模型为
$$ \begin{gathered} {k_{{\text{app}}}} = \frac{{L_0^{1 - {D_{\text{T}}}}\left( {2 - {D_{\text{f}}}} \right){{\left( {1 - \xi } \right)}^2}b_{\max }^{1 + {D_{\text{T}}}}}}{{12\left( {3 + {D_{\text{T}}} - {D_{\text{f}}}} \right)}}\frac{\phi }{{1 - \phi }} \\ {\left\{ \begin{gathered} 1 - {C_{\text{f}}}\left[ {\left( {\sigma - {\sigma _{\text{0}}}} \right) - \left( {P - {P_0}} \right)} \right] \\ - \frac{{3f}}{{{\phi _0}}}\left( {\frac{{{S_{\max }}P}}{{{P_L} + P}} - \frac{{{S_{\max }}{P_0}}}{{{P_L} + {P_0}}}} \right) \\ \end{gathered} \right\}^3} \times \\ \left\{ {1 + {{{B_0}} \mathord{\left/ {\vphantom {{{B_0}} {\left\{ {P\left\{ \begin{gathered} 1 - {C_{\text{f}}}\left[ {\left( {\sigma - {\sigma _{\text{0}}}} \right) - \left( {P - {P_{\text{0}}}} \right)} \right] \\ - \frac{{3f}}{{{\phi _0}}}\left( {\frac{{{S_{{\text{max}}}}P}}{{{P_{\text{L}}} + P}} - \frac{{{S_{{\text{max}}}}{P_{\text{0}}}}}{{{P_{\text{L}}} + {P_{\text{0}}}}}} \right) \\ \end{gathered} \right\}} \right\}}}} \right. } {\left\{ {P\left\{ \begin{gathered} 1 - {C_{\text{f}}}\left[ {\left( {\sigma - {\sigma _{\text{0}}}} \right) - \left( {P - {P_{\text{0}}}} \right)} \right] \\ - \frac{{3f}}{{{\phi _0}}}\left( {\frac{{{S_{{\text{max}}}}P}}{{{P_{\text{L}}} + P}} - \frac{{{S_{{\text{max}}}}{P_{\text{0}}}}}{{{P_{\text{L}}} + {P_{\text{0}}}}}} \right) \\ \end{gathered} \right\}} \right\}}}} \right\} \\ \end{gathered} $$ (19) 3.2 模型验证与分析
此前，作者团队针对贵州六盘水地区的煤样开展了恒定外应力条件下的煤岩渗透率测试，并建立了考虑煤岩裂隙分形特征的渗透率模型，分析了煤岩孔裂隙下分形维数和渗透率的耦合机理[22]。然而，该研究尚未考虑煤岩裂隙表面粗糙度的影响。在本文中，基于分形理论定量表征煤岩裂隙表面相对粗糙度，建立了考虑裂隙表面相对粗糙度、吸附应力及滑脱效应的表观渗透率模型。基于渗透率测试结果，通过匹配试验结果与理论结果，进一步探讨并量化分析煤岩裂隙表面相对粗糙度及滑脱系数与渗透率的响应机理。试验条件为恒定有效应力，存在(σ‒σ0)‒(P‒P0)=0，煤岩内气体运移通道大致为一条直线，DT=1。为了验证考虑裂隙表面相对粗糙度的煤岩表观渗透率模型的合理性，将粗糙微元最小与最大底边长度的比值ε、粗糙微元的最大高度与裂隙开度的比值hmax/b、φ代入式（6）中，粗糙微元底边长度分布的分形维数Da的值可通过式Da=dE−(lnφ/lnε)获得，可得煤岩裂隙表面相对粗糙度ξ。进一步将Smax、PL、f、Df、bmax（见表5）代入式（19）得到模型计算结果与试验数据对比如图6所示，球形实点为试验数据、曲线为拟合数据。
表 5 模型验证参数Table 5. Model validation parameters由图6可知，所建表观渗透率模型与试验数据具有相同的变化趋势，表明所建立的表观渗透率模型可较好地表征裂隙的表观渗透率随孔隙压力的变化关系。同时可知煤岩表观渗透率随孔隙压力的增大呈先减小并逐渐趋于平缓的变化趋势。究其原因：随孔隙压力增大前期，滑脱效应占主导，由于滑脱效应导致裂隙表面存在一定附加流量而导致煤岩表观渗透率在孔隙压力增大前期较大；随孔隙压力增大后期，气体吸附效应逐渐占主导，而由于煤层气逐渐吸附于裂隙壁面占据裂隙通道，裂隙通道减小而增大渗流阻力，煤岩表观渗透率逐渐减小。然而煤岩表观渗透率逐渐趋于平缓的变化趋势，究其原因：气体吸附效应随孔隙压力的增大，气体吸附量逐渐达到饱和，基质膨胀变形量逐渐减小而对裂隙影响减小并逐渐趋于可忽略状态。此外，通过与WANG[22]等计算值的比较，发现所提模型计算值与其的变化趋势相似以及计算值与其在同一数量级，侧面证明笔者所建模型的可靠性与合理性。
![]()
图 6 煤岩表观渗透率与孔隙压力的演化关系Figure 6. Evolution of apparent permeability and pore pressure in coal3.3 煤岩裂隙表面相对粗糙度与滑脱系数的动态演化
由式（18）可知，滑脱系数主要与煤岩裂隙表面相对粗糙度、初始滑脱系数和有效裂隙开度beff及裂隙开度相关。为探究煤岩裂隙表面相对粗糙度对滑脱系数的影响，将裂隙表面相对粗糙度ξ分别设为0.004，0.006，0.008，0.01，得到不同孔隙压力下滑脱系数Beff与ξ的演化规律，如图7所示。
![]()
图 7 不同孔隙压力下相对粗糙度与滑脱系数的变化关系Figure 7. Relationship between relative roughness and slippage coefficient under different pore pressures由图7可知，滑脱系数随煤岩裂隙表面相对粗糙度的增大而增大，究其原因：一方面，滑脱效应是指由于煤层气运移过程中与裂隙壁面碰撞而促使壁面产生的附加流速，相对粗糙度越大，气体分子与壁 面碰撞越多，滑脱效应增大。即对气体分子产生相应的动力，促使煤层气分子运移速度增大，滑脱效应增大。另一方面，裂隙表面相对粗糙度增大会导致煤岩有效裂隙开度减小，进而导致煤层气体分子与壁面的碰撞频率增大从而导致滑脱效应增大。
3.4 滑脱系数与煤岩表观渗透率的响应机制
滑脱效应是指气体运移过程中与裂隙壁面频繁碰撞而在裂隙壁面产生的附加流速，导致煤层气运移速度增大而促使煤岩表观渗透率增大。为探究煤岩裂隙表面相对粗糙度对滑脱系数的影响，煤岩内气体运移通道可以近似于直线，即DT=1，结合式（14）、式（15）、式（17）和式（18）可得煤岩粗糙度与滑脱系数的关系式
$$ \begin{gathered} {k_{{\text{app}}}} = \frac{{\left( {2 - {D_{\text{f}}}} \right){{\left( {1 - \xi } \right)}^2}b_{\max }^2}}{{12\left( {4 - {D_{\text{f}}}} \right)}}\frac{\phi }{{1 - \phi }} \times \\ {\left( {\frac{{{b_{{\text{eff}}}}}}{{{b_{\text{0}}}}}} \right)^3} \times \left( {1 + \frac{{{B_{\text{0}}}}}{{P\left( {{{{b_{{\text{eff}}}}} \mathord{\left/ {\vphantom {{{b_{{\text{eff}}}}} {{b_{\text{0}}}}}} \right. } {{b_{\text{0}}}}}} \right)}}} \right) \\ \end{gathered} $$ (20) 图8所示为不同初始状态下滑脱系数B0下，粗糙裂隙的表观渗透率随孔隙压力的变化关系图。由图8可知，相同初始滑脱系数条件下，煤岩表观渗透率均随孔隙压力增大而减小并趋于平缓。孔隙压力增大阻碍气体通过煤岩渗流通道，且孔隙压力越小时，阻碍作用越明显。而相同孔隙压力条件下，随初始滑脱系数B0的增大，煤岩表观渗透率随之增大。此时，滑脱系数增大，气体分子碰撞煤壁的概率增大，即滑脱效应增强，滑脱效应对煤岩渗透起促进作用，所以，煤岩表观渗透率也随之增大。滑脱系数越大，其作用效应越强。
![]()
图 8 不同B0条件下煤岩表观渗透率与孔隙压力的关系Figure 8. Relationship between apparent permeability and pore pressure of coal with different B0 values4. 结 论
1）将扫描电镜与低温氮气吸附试验相结合，对煤岩孔裂隙进行联合表征，通过扫描电镜试验观测到孔裂隙结构表面存在粗糙微元，同时低温氮气吸附试验与FHH模型结合结果表明，煤岩孔裂隙结构发育，存在分形特征且非均质性较强。
2）粗糙微元分形维数的含义表示煤岩裂隙表面相对粗糙度的粗糙程度，相对粗糙度会随粗糙微元分形维数的增大而增大。当粗糙微元分形维数恒定时，相对粗糙度随粗糙微元的尺度比例ε的增大而减小，煤岩裂隙开度b越小，裂隙表面相对粗糙度对裂隙开度越敏感。
3）随着煤岩裂隙表面相对粗糙度增加，滑脱系数也随之增加。一方面，相对粗糙度越大，气体分子 与壁面碰撞越多，导致滑脱效应增大。另一方面相对粗糙度增大会导致煤岩有效裂隙开度减小，进而导致煤层气体分子与壁面的碰撞频率增大从而导致滑脱效应增大。
4）在同一初始滑脱系数下，煤岩表观渗透率随孔隙压力升高而减小且趋于平缓。孔隙压力增大阻碍气体通过煤岩渗流通道，且孔隙压力越小时，阻碍作用越明显。相同孔隙压力条件下，随初始滑脱系数B0的增大，煤岩表观渗透率随之增大。滑脱效应对煤岩渗透起促进作用，滑脱系数越大，其作用效应越强。
-
![]()
![]()
图 2 工业组分及全硫含量垂向分布
Figure 2. Vertical distribution of proximate components and total sulfur content
表 1 煤样工业分析和全硫分析
Table 1 Proximate analysis and total sulfur content of coal samples
煤层号 样品号 工业分析/% 元素分析St,d/% Mad Ad Vdaf 5煤 RJZ5-1 0.6 36 40 1.87 RJZ5-4 0.3 15 37 0.96 RJZ5-6 0.6 32 35 0.53 RJZ5-7 0.6 24 35 3.63 RJZ5-8 1.3 34 37 3.67 RJZ5-9 0.9 17 33 1.01 RJZ5-10 0.6 33 36 0.59 RJZ5-11 0.5 23 35 1.49 RJZ5-12 1.6 20 38 0.91 RJZ5-13 0.3 12 33 0.77 RJZ5-14 0.4 41 41 0.71 9上煤 RJZ9上-1 0.1 5 37 2.14 RJZ9上-2 2.3 49 47 25.16 RJZ9上-3 0.9 3 38 2.83 RJZ9上-4 1.3 4 37 2.78 RJZ9上-5 0.2 3 42 3.15 RJZ9上-6 1.0 4 44 3.82 9煤 RJZ9-1 0.6 26 39 5.25 RJZ9-2 0.8 14 41 3.34 RJZ9-3 0.6 49 44 1.82 RJZ9-6 0.7 38 46 1.93 RJZ9-7 1.8 5 38 3.05 RJZ9-8 0.3 5 40 3.09 RJZ9-10 0.8 46 33 1.08 RJZ9-11 0.8 30 31 1.25 注：下标ad为空气干燥基；d为干燥基；daf为干燥无灰基；St,d为全硫。 
下载: 导出CSV
表 2 样品中常量元素氧化物质量分数（全煤基）和烧失量
Table 2 Mass fraction and Loss on ignition of major element oxides in samples (on whole coal basis)
% 样品 质量分数 烧失量 SiO2 TiO2 Al2O3 Fe2O3 MnO MgO CaO Na2O K2O P2O5 RJZ5-1 24.69 0.815 7.48 2.01 0.004 0.148 0.126 0.050 0.716 0.015 0.24 RJZ5-2 61.79 0.762 16.92 4.95 0.009 0.449 0.062 0.123 2.341 0.030 0.20 RJZ5-3 43.24 0.800 20.53 8.05 0.005 0.393 0.154 0.162 1.763 0.038 0.28 RJZ5-4 7.67 0.198 6.47 0.28 0.001 0.030 0.164 0.017 0.021 0.008 0.32 RJZ5-5 42.20 1.294 36.43 0.25 0.001 0.081 0.065 0.162 0.122 0.025 0.10 RJZ5-6 16.81 0.328 14.16 0.17 0.001 0.064 0.074 0.038 0.093 0.024 0.08 RJZ5-7 10.32 0.259 8.78 4.15 0.003 0.089 0.182 0.026 0.022 0.021 0.20 RJZ5-8 15.22 0.591 13.06 4.18 0.003 0.150 0.282 0.037 0.082 0.026 0.28 RJZ5-9 7.19 0.376 7.02 0.50 0.003 0.148 0.323 0.037 0.014 0.410 0.68 RJZ5-10 16.68 1.027 14.56 0.32 0.001 0.063 0.073 0.036 0.092 0.034 0.14 RJZ5-11 10.74 0.303 9.71 1.39 0.002 0.067 0.315 0.032 0.021 0.222 0.36 RJZ5-12 8.11 0.278 7.15 0.22 0.009 0.970 1.754 0.070 0.012 0.109 0.97 RJZ5-13 6.05 0.091 5.37 0.09 0.001 0.053 0.110 0.023 0.008 0.020 0.50 RJZ5-14 20.60 0.943 18.00 0.72 0.002 0.098 0.111 0.045 0.107 0.030 0.22 RJZ9上-顶 22.69 0.369 11.18 24.00 0.122 0.339 2.925 0.096 1.056 0.073 0.46 RJZ9上-1 2.93 0.211 1.40 0.08 0.001 0.033 0.178 0.033 0.012 0.002 0.52 RJZ9上-2 2.01 0.095 1.44 45.13 0.009 0.064 0.201 0.157 0.015 0.005 0.04 RJZ9上-3 1.11 0.031 0.97 0.07 0.002 0.168 0.305 0.127 0.009 0.002 1.34 RJZ9上-4 1.60 0.101 1.32 0.13 0.003 0.154 0.310 0.073 0.006 0.002 1.38 RJZ9上-5 0.58 0.027 0.55 0.11 0.002 0.378 0.577 0.049 0.004 0.001 1.38 RJZ9上-6 0.33 0.012 0.35 0.51 0.003 0.560 0.848 0.058 0.007 0.002 0.70 RJZ9上-底 94.29 0.119 2.71 0.70 0.006 0.099 0.089 0.099 0.564 0.016 0.11 RJZ9-顶 66.28 1.243 18.76 1.37 0.005 0.209 0.073 0.100 2.111 0.056 0.13 RJZ9-1 10.46 0.324 7.00 4.14 0.010 0.801 1.196 0.055 0.385 0.014 0.88 RJZ9-2 7.91 0.131 3.21 1.04 0.003 0.302 0.435 0.080 0.151 0.006 1.44 RJZ9-3 24.74 1.018 20.98 0.47 0.004 0.294 0.397 0.103 0.098 0.039 0.55 RJZ9-4 28.32 0.950 24.52 0.18 0.004 0.398 0.594 0.101 0.067 0.030 0.57 RJZ9-5 43.96 0.710 37.39 0.07 0.002 0.083 0.075 0.125 0.116 0.017 0.20 RJZ9-6 19.47 0.292 16.74 0.13 0.003 0.312 0.418 0.118 0.038 0.012 0.76 RJZ9-7 1.74 0.032 1.52 0.14 0.003 0.335 0.504 0.158 0.009 0.023 1.68 RJZ9-8 1.77 0.034 1.57 0.34 0.004 0.329 0.524 0.094 0.006 0.029 0.51 RJZ9-9 26.80 0.657 23.24 4.59 0.004 0.134 0.146 0.090 0.078 0.020 0.11 RJZ9-10 22.73 0.651 20.57 0.44 0.001 0.074 0.161 0.110 0.032 0.536 0.23 RJZ9-11 15.51 0.342 13.50 0.20 0.001 0.045 0.087 0.045 0.024 0.040 0.26
下载: 导出CSV
表 3 样品中微量元素质量分数（全煤基）
Table 3 Concentration of trace elements in samples (on whole coal basis)
μg/g 样品 质量分数 Li Be Sc V Cr Co Ni Cu Zn Ga Rb Sr Y Zr Nb Sn Cs Ba La RJZ5-1 48.6 9.5 20.2 57.2 59.8 17.4 29.3 18.7 14.2 37.8 28.4 261.7 31.0 218.5 13.2 4.2 2.4 135.7 23.8 RJZ5-2 53.0 2.6 11.6 97.7 202.4 29.9 51.2 13.7 59.0 23.9 118.8 68.2 27.8 216.5 16.3 4.0 9.4 336.2 35.5 RJZ5-3 91.6 3.8 13.3 126.3 162.5 21.7 73.5 21.3 108.6 26.6 92.4 112.4 31.5 246.4 19.6 5.1 8.8 7 148.7 44.0 RJZ5-4 42.9 7.5 6.6 20 8.0 0.9 3.1 11.5 6.1 14.6 1.0 23 22.8 132.9 6.5 1.0 0.1 50.3 9.2 RJZ5-5 469.0 3.2 5.4 16.9 6.3 0.4 2.9 5.8 8.7 39.8 6.2 53.5 11.7 154.7 47.4 2.5 0.3 29.2 17.2 RJZ5-6 123.5 3.5 7.9 26.4 8.3 0.5 3.6 10.2 4.4 21.1 5.6 31 21.9 277.1 15.6 2.8 0.6 27.5 32.0 RJZ5-7 75.8 3.3 8.3 21.5 66.5 5.2 11.9 18.3 7.2 13.1 1.1 73.2 19.5 145.7 9 1.8 0.1 34.8 12.3 RJZ5-8 123.4 3.1 9.4 22.4 13.9 9.4 12.4 11.2 7.5 18.8 4.7 45.2 21.7 171.7 20.1 2.4 0.5 21.9 25.2 RJZ5-9 76.5 2.7 5.3 23.0 6.4 1.3 2.4 8.9 7.2 15.7 0.6 3213 19.7 88.9 7.1 1.1 0 107.4 228.5 RJZ5-10 157.1 2.3 10 31.3 12.0 0.5 2.8 13.7 6.8 17.8 4.9 60.7 22.4 181.8 24.7 3.4 0.4 32.4 15.7 RJZ5-11 116.4 1.8 6.5 16.1 30.8 1.4 2.5 11.6 5.4 18.6 1.2 544.9 17.7 134.1 8.6 1.1 0.1 35.0 32.7 RJZ5-12 84.4 1.4 6.9 26.0 8.4 1.6 2.1 9.0 5.4 16.0 0.5 273.0 23.6 278.6 10.2 1.3 0 50.6 12.1 RJZ5-13 52.3 1.5 4.9 12.1 4.4 1.1 1.4 8.9 3.7 10.5 0.4 120.7 13.2 71.2 3.6 0.5 0 27.5 13.7 RJZ5-14 175.5 2.1 6.8 20 13.4 0.9 4.7 10.6 6.9 14.5 5.9 102.1 17.0 208.3 23.9 2.9 0.7 34.0 13.9 RJZ9上-顶 41.4 1.3 4.5 234.9 131.8 13.5 152.3 19.5 154.9 11.7 45.6 319.4 9.9 92.8 7.4 2.2 3.3 135.0 50.5 RJZ9上-1 4.3 4.2 4.5 34.8 16.6 0.6 6.7 15 4.5 8.5 0.5 38.6 8.2 149.9 5.5 0.6 0 22.5 2.9 RJZ9上-2 5.9 1.1 0.9 6.8 481.2 1.3 7.6 7.2 8.2 1.6 1.7 36.2 3.6 20 2.4 0.4 0.1 30.2 6.0 RJZ9上-3 4.4 1.6 1.0 3.2 9.0 0.8 1.1 4.1 2.8 4.4 0.3 51.1 3.5 21.3 1.0 0.4 0 23.3 3.8 RJZ9上-4 6.4 0.7 1.4 8.2 5.6 0.7 1.4 10 2.8 2.6 0.2 59.8 3.6 36.3 1.9 0.3 0 21.2 4.1 RJZ9上-5 9.8 0.4 1.0 2.6 8.4 0.7 1.5 2.2 0.8 1.6 0.1 6.8 2.4 18.5 1.1 0.2 0 3.6 1.6 RJZ9上-6 2.1 0.5 2.5 2.2 5.0 1.3 2.0 4.8 4.1 20.6 0.2 116.9 4.0 6.0 0.4 3.5 0 21.8 0.9 RJZ9上-底 3.5 0.1 0.8 9.4 594 1.7 7.8 3.1 7.8 2.9 13.8 41.9 4.0 102.0 2.4 0.3 0.2 22.8 7.4 RJZ9-顶 111.9 1.2 9.6 60.5 415.9 15.6 36.6 14.9 85.0 19.5 62.6 154.7 31.7 689.8 26.7 1.9 2.6 429.5 53.3 RJZ9-1 53.3 1.3 10.8 72.8 53.6 9.7 24.2 12.6 39.5 18.8 11.8 212.2 12.4 220.7 8.6 1.2 0.5 111.3 12.9 RJZ9-2 18.8 1.2 3.8 10.2 15.4 1.8 4.5 6.1 8.8 13.6 4.3 116.9 10.5 56.4 3.5 1.2 0.2 52.8 8.0 RJZ9-3 287.6 1.3 13.0 52.9 27.8 1.1 5.8 52.4 7.8 26.7 4.3 131.8 24.2 345.5 31.5 5.1 0.2 73.0 45.8 RJZ9-4 221.8 1.6 18.7 40.3 19.5 0.6 3.4 33.0 8.7 20.5 3.0 144.5 29.6 261 27.2 5.4 0.2 77.8 35.0 RJZ9-5 218.3 1.2 6.8 10.5 6.5 0.5 1.6 7.6 7.1 35.4 6.1 40.9 8.0 175.1 34.5 2.7 0.6 52.8 8.0 RJZ9-6 120.1 1.3 8.2 20.4 6.8 1.9 2.6 13.8 3.8 27.4 1.9 117.4 16.5 370.9 13.7 2.1 0.2 30 8.2 RJZ9-7 19.1 0.8 1.8 2.9 2.6 0.6 1.1 7.9 2.8 4.0 0.4 158.5 7.0 18.0 0.7 1.0 0 26.2 5.9 RJZ9-8 4.6 1.5 3.1 9.5 9.5 2.9 5.0 2.1 1.3 4.3 5.1 15.6 4.2 34.3 2.1 0.5 0.4 17.3 2.1 RJZ9-9 289 1.6 7.4 17.9 81.8 1.3 5.3 10.2 8.7 25.3 3.8 79.5 16.7 177.5 18.5 4.2 0.3 56.6 8.3 RJZ9-10 347.8 5.4 15.7 21.8 14.3 0.2 1.5 19.5 6.7 16.9 1.6 2 395.2 54.7 240.6 17.8 4.4 0.1 234.1 235.5 RJZ9-11 227.7 3.1 9.6 12.9 10 0.3 1.7 28.6 6.3 9.5 1.1 164.4 29.9 128.1 9.6 2.3 0.1 60 44.7
下载: 导出CSV
续表3 样品 质量分数 Ce Pr Nd Sm Eu Gd Tb Dy Ho Er Tm Yb Lu Hf Ta Tl Pb Th U RJZ5-1 40.3 4.5 15.9 3.5 0.6 3.5 0.7 4.6 1.0 3.3 0.5 3.6 0.5 4.8 0.9 1.3 24.3 12.0 3.7 RJZ5-2 68.3 7.6 25.9 4.8 0.8 3.9 0.7 4.7 0.9 2.8 0.4 2.7 0.4 5.9 1.2 6.4 27.2 11.9 4.0 RJZ5-3 84.7 9.3 31.7 5.9 1.8 4.9 0.8 5.3 1.1 3.2 0.5 3.2 0.5 6.5 1.3 3.7 68.6 17.2 4.8 RJZ5-4 22.4 2.8 10.7 2.6 0.4 2.5 0.5 3.2 0.7 2.4 0.4 2.5 0.4 3.5 0.4 0.2 8.7 7.6 1.8 RJZ5-5 27.6 2.6 7.6 1.3 0.2 1.2 0.3 1.9 0.4 1.1 0.2 1.1 0.2 5.6 1.7 0.1 24.6 12.7 4.9 RJZ5-6 68.2 7.6 26.5 5.2 0.7 3.8 0.6 3.9 0.7 2.1 0.3 2.0 0.3 8.2 1.6 0.1 33.3 26.9 3.3 RJZ5-7 26.2 3.1 11.7 2.8 0.5 2.6 0.5 3.4 0.7 2.1 0.3 2.0 0.3 4.2 0.7 2.8 47.3 12.1 3.3 RJZ5-8 53.0 6.2 21.9 4.4 0.7 3.7 0.6 4.1 0.8 2.4 0.3 2.3 0.3 5.0 1.3 1.1 60.2 18.5 5.5 RJZ5-9 330.8 31.1 99.6 16.7 2.6 11.3 1.2 4.8 0.7 1.9 0.2 1.6 0.2 2.2 0.4 0.1 10.5 6.9 2.8 RJZ5-10 38.6 4.8 17.5 3.9 0.7 3.2 0.6 3.9 0.8 2.4 0.4 2.3 0.3 5.0 1.6 0.1 26.3 13.5 6.3 RJZ5-11 64.6 7.2 25.5 4.6 0.9 3.9 0.6 3.6 0.6 1.7 0.2 1.6 0.2 3.5 0.5 0.4 14.7 8.3 2.9 RJZ5-12 23.2 2.7 10.2 2.6 0.6 3.3 0.6 4.0 0.8 2.2 0.3 1.8 0.3 7.2 0.5 0.1 15.2 12.0 2.9 RJZ5-13 23.8 2.6 9.7 2.1 0.4 1.9 0.3 2.2 0.4 1.3 0.2 1.2 0.2 1.8 0.1 0.1 8.4 4.2 1.7 RJZ5-14 27.7 3.5 13.1 2.6 0.5 2.4 0.5 2.8 0.5 1.7 0.3 1.7 0.2 6.2 1.8 0.3 26.5 23.8 5.0 RJZ9上-顶 105 11.0 35.1 4.8 0.8 2.6 0.4 1.9 0.3 1.1 0.2 1.1 0.2 2.7 0.5 12.0 21.4 9.7 25.2 RJZ9上-1 6.7 0.8 2.9 0.7 0.2 0.8 0.2 1.2 0.3 0.9 0.1 1 0.2 2.3 0.3 0 4.4 3.9 22.0 RJZ9上-2 13.9 1.6 5.7 1.1 0.2 0.9 0.1 0.6 0.1 0.3 0 0.3 0 0.5 0.3 0.5 2.7 5.2 2.8 RJZ9上-3 8.0 0.9 3.2 0.6 0.1 0.5 0.1 0.5 0.1 0.3 0 0.3 0 0.6 0 0 1.6 0.6 0.5 RJZ9上-4 8.4 1.0 3.5 0.7 0.1 0.6 0.1 0.6 0.1 0.4 0.1 0.3 0 0.9 0.1 0 3.1 2.2 0.8 RJZ9上-5 3.3 0.4 1.5 0.4 0.1 0.3 0.1 0.4 0.1 0.3 0 0.3 0 0.5 0.1 0.4 6.0 1.6 0.4 RJZ9上-6 2.0 0.3 1.1 0.3 0.1 0.4 0.1 0.7 0.1 0.4 0.1 0.3 0 0.1 0 0 0.2 0.3 0.3 RJZ9上-底 14.1 1.6 5.8 1.1 0.2 0.9 0.1 0.7 0.2 0.4 0.1 0.5 0.1 2.7 0.2 0.1 4.7 2.1 0.7 RJZ9-顶 101.9 12.4 44.0 7.6 1.3 5.8 0.9 5.4 1.0 3.3 0.5 3.3 0.5 18.1 1.7 0.5 21.1 16.9 7.5 RJZ9-1 25.9 3.1 11.7 2.4 0.5 2.2 0.4 2.3 0.5 1.4 0.2 1.2 0.2 5.2 0.4 0.5 8.3 6.3 18.2 RJZ9-2 16.4 2.0 7.3 1.5 0.3 1.5 0.3 1.6 0.3 1.0 0.1 0.9 0.1 1.6 0.3 0.1 5.1 3.8 1.6 RJZ9-3 93.1 10.4 35.7 6.7 1.2 5.6 1.0 5.6 1.0 2.7 0.4 2.5 0.4 9.6 2.4 0.1 50 38.7 12.8 RJZ9-4 81.8 9.6 33.9 7.2 1.3 6.1 1.1 6.8 1.2 3.4 0.5 3.1 0.4 9.2 2.2 0.1 76.7 53.8 11.3 RJZ9-5 17.5 1.8 6.3 1.3 0.2 1.2 0.2 1.5 0.3 0.8 0.1 0.8 0.1 6.8 2.3 0.1 22.8 17.2 7.3 RJZ9-6 19.4 2.4 9.3 2.3 0.4 2.3 0.5 3.2 0.6 1.8 0.3 1.7 0.3 9.4 1.1 0.1 35.5 19.8 10.6 RJZ9-7 13.2 1.6 6.3 1.3 0.2 1.2 0.2 1.2 0.2 0.7 0.1 0.7 0.1 0.5 0.1 0 2.1 1.5 0.5 RJZ9-8 3.9 0.4 1.6 0.3 0.1 0.4 0.1 0.6 0.1 0.5 0.1 0.5 0.1 0.7 0.1 0.4 3.3 1.8 0.5 RJZ9-9 18.8 2.2 8.5 2.1 0.4 2.2 0.5 2.9 0.6 1.6 0.3 1.5 0.2 6.0 1.5 0.6 42.2 19.5 5.5 RJZ9-10 366.2 35.7 115 21.0 3.4 16.1 2.3 12.0 2.0 5.4 0.7 4.1 0.6 7.2 1.3 0.1 39.7 27.2 5.3 RJZ9-11 79.2 8.2 27.3 5.1 0.9 4.5 0.8 5.1 1.0 2.9 0.4 2.3 0.3 3.8 0.9 0.1 34.8 17.1 3.6
下载: 导出CSV
表 4 煤中微量元素富集程度统计
Table 4 Enrichment degree of trace elements in coal
煤层 高度富集 富集 轻度富集 正常范围 亏损 5 Li、Ba Rb、Zr Be、Sc、V、Cr、Ni、Zn、Ga、Sr、
Y、Nb、Sn、Cs、La、Ce、Pr、Nd、
Sm、Er、Yb、Hf、Ta、Pb、ThCo、Cu、Eu、Gd、Tb、Dy、Ho、Tm、Lu、Tl、U — 9上 — Cr Tl、U Li、Zr、Hf、Th、Sc、Ga、Sr、Y、Nb、
Sn、La、Ce、Pr、Nd、Sm、Er、Ta、
Pb、Be、V、Co、Ni、Cu、Zn、Rb、YbCs、Ba、Eu、Gd、Tb、
Dy、Ho、Tm、Lu9 Li Zr、Hf、
ThSc、Cr、Ga、Sr、Y、Nb、Sn、La、
Ce、 Pr、Nd、Sm、Tb、Er、Ta、Pb、UBe、V、Co、Ni、Cu、Zn、Rb、Ba、Eu、
Gd、Dy、Ho、Tm、Yb、LuCs、Tl
下载: 导出CSV
-
[1] 宁树正, 邓小利, 李聪聪, 等. 中国煤中金属元素矿产资源[M]. 北京: 科学出版社, 2019. [2] 宁树正,邓小利,李聪聪,等. 中国煤中金属元素矿产资源研究现状与展望[J]. 煤炭学报,2017,42(9):2214−2225. doi: 10.13225/j.cnki.jccs.2017.0683 NING Shuzheng,DENG Xiaoli,LI Congcong,et al. Research status and prospect of metal element mineral resources in China[J]. Journal of China Coal Society,2017,42(9):2214−2225. doi: 10.13225/j.cnki.jccs.2017.0683
[3] 代世峰,赵 蕾,魏 强,等. 中国煤系中关键金属资源: 富集类型与分布[J]. 科学通报,2020,65:3715−3729. DAI Shifeng,ZHAO Lei,WEI Qiang,et al. Resources of critical metals in coal-bearing sequences in China: Enrichment types and distribution (in Chinese)[J]. Chin Sci Bull,2020,65:3715−3729.
[4] 魏 强,代世峰. 中国煤型锗矿床中的关键金属和有害元素: 赋存特征与富集成因[J]. 煤炭学报,2020,45(1):296−303. WEI Qiang,DAI Shifeng. Critical metals and hazardous elements in the coal-hosted germanium ore deposits of China: Occurrence characteristics and enrichment causes[J]. Journal of China Coal Society,2020,45(1):296−303.
[5] 代世峰,刘池洋,赵 蕾,等. 煤系中战略性金属矿产资源: 意义和挑战[J]. 煤炭学报,2022,47(5):1743−1749. DAI Shifeng,LIU Chiyang,ZHAO Lei,et al. Strategic metal resources in coal-bearing strata: Significance and challenges[J]. Journal of China Coal Society,2022,47(5):1743−1749.
[6] 李 晶,王 园,袁 伟,等. 陕西典型矿区煤系关键金属元素富集特征及其成因机制[J]. 煤炭学报,2022,47(5):1808−1821. doi: 10.13225/j.cnki.jccs.mj22.0225 LI Jing,WANG Yuan,YUAN Wei,et al. Enrichment characteristics and mechanism of critical metal elements enriched in coals from typical coalfields in Shaanxi Province[J]. Journal of China Coal Society,2022,47(5):1808−1821. doi: 10.13225/j.cnki.jccs.mj22.0225
[7] 赵存良. 鄂尔多斯盆地与煤伴生多金属元素的分布规律和富集机理[D]. 北京: 中国矿业大学(北京), 2015. ZHAO Cunliang. Distribution and enrichment mechanism of multi-metallic elements associated with coal in Ordos Basin [D]. Beijing: China University of Mining & Technology-Beijing, 2015.
[8] 秦国红,邓丽君,刘 亢,等. 鄂尔多斯盆地西缘煤中稀土元素特征[J]. 煤田地质与勘探,2016,44(6):8−14. doi: 10.3969/j.issn.1001-1986.2016.06.002 QIN Guohong,DENG Lijun,LIU Kang,et al. Characteristic of rare earth elements in coal in western margin of Ordos basin[J]. Coal Geology & Exploration,2016,44(6):8−14. doi: 10.3969/j.issn.1001-1986.2016.06.002
[9] 秦国红. 鄂尔多斯盆地晚古生代煤中微量元素富集特征与成因类型[D]. 北京: 中国矿业大学(北京), 2020. QIN Guohong. Enrichment characteristics and genetic types of trace elements in the Late Paleozoic Coal from Ordos Basin[D]. Beijing: China University of Mining & Technology-Beijing, 2020.
[10] 石志祥,马家亮,高 政,等. 宁东煤田中部侏罗纪煤的矿物学特征[J]. 煤炭学报,2017,42(6):1535−1543. SHI Zhixiang,MA Jialiang,GAO Zheng,et al. Mineralogical of the Jurassic coals in the middle of Ningdong Coalfield[J]. Journal of China Coal Society,2017,42(6):1535−1543.
[11] 曹代勇, 魏迎春, 等. 鄂尔多斯盆地煤系矿产赋存规律与资源评价[M]. 北京: 科学出版社, 2019. [12] 吴 蒙,秦云虎,王晓青,等. 宁东地区晚古生代煤中硫的地化特征及其对有害元素富集的影响[J]. 煤炭学报,2020,45(S2):932−942. WU Meng,QIN Yunhu,WANG Xiaoqing,et al. Geochemical characteristics of sulfur and its impact on accumulation of hazardous trace elements in late Paleozoic coal from Ningdong area[J]. Journal of China Coal Society,2020,45(S2):932−942.
[13] 吴 亮,何 伟. 宁夏宁东煤田灵武矿区煤中砷含量及赋存形态探讨[J]. 煤炭与化工,2021,44(4):57−60,117. WU Liang,HE Wei. Discussion on arsenic content and occurrence form in Ningdong coalfield in Ningxia[J]. Coal and Chemical Industry,2021,44(4):57−60,117.
[14] 何 伟,吴 亮,魏向成,等. 宁东煤田中侏罗统延安组稀有稀散稀土元素地球化学特征及其对沉积环境的指示意义[J]. 岩矿测试,2022,41(6):962−977. HE Wei,WU Liang,WEI Xiangcheng,et al. Geochemical characteristics of three kinds of rare elements in middle jurassic yan’an formation of ningdong coalfield and their indicative significance to sedimentary environment[J]. Rockand Mineral Analysis,2022,41(6):962−977.
[15] 刘 亢. 鄂尔多斯盆地西缘煤系矿产资源共生组合特征研究[D]. 北京: 中国矿业大学(北京), 2016. LIU Kang. Combination characters of coal series mineral resources in the west margin of ordos basin[D]. Beijing: China University of Mining & Technology-Beijing, 2016.
[16] GB/T 482—2008, 煤层煤样采取方法[S]. [17] 王双明. 鄂尔多斯盆地叠合演化及构造对成煤作用的控制[J]. 地学前缘, 2017, 24(2): 54-63. WANG Shuangming. Ordos Basin superposed evolution and structural controls of coal forming activities. Earth Science Frontiers, 2017, 24(2): 54-63.
[18] GB/T212—2008, 煤的工业分析方法[S]. [19] GB/T214—2007, 煤中全硫的测定方法[S]. [20] QUEROL X,WHATELEY M K G,FERNANDEZ-TURIEL L,et al. Geological controls on the mineralogy and geochemistry of the Bey pazari lignite, central Anatolia, Turkey[J]. International Journal of Coal Geology,1997,33(3):255−271. doi: 10.1016/S0166-5162(96)00044-4
[21] DAI S, SEREDIN V V, WARD C R, et al. Enrichment of U-Se-Mo - Re - V in coals preserved within marine carbonate successions: geochemical and mineralogical data from the Late Permian Guiding coalfield, Guizhou, China[J]. Mineralium Deposita, 2015, 50(2): 159-186.
[22] KETRIS M P,YUDOVICH Y E. Estimations of Clarkes for Carbonaceous biolithes: World averages for trace element contents in black shales and coals[J]. International Journal of Coal Geology,2009,78(2):135−148. doi: 10.1016/j.coal.2009.01.002
[23] 宁树正, 朱士飞, 曹代勇, 等. 2017. 含煤岩系矿产综合勘查评价研究(2015-2016年度)[R]. 北京: 中国煤炭地质总局. [24] 代世峰,任徳贻,周义平等. 煤型稀有金属矿床: 成因类型、赋存状态和利用评价[J]. 煤炭学报,2014,39(8):1707−1715. DAI Shifeng,REN Deyi,ZHOU Yiping,et al. Coal-hosted rare metal deposits: Genetic types, modes of occurrence, and utilization evaluation[J]. Journal of China Coal Society,2014,39(8):1707−1715.
[25] 熊 科, 何田田. 2019. 宁夏宁鲁煤电有限责任公司任家庄煤矿九号煤煤质补充勘探报告[R]. 保定: 中国煤炭地质总局一七三勘探队. [26] 田和明,代世峰,李大华,等. 重庆南川晚二叠世凝灰岩的元素地球化学特征[J]. 地质论评,2014,60(1):169−177. doi: 10.16509/j.georeview.2014.01.019 TIAN Heming,DAI Shifeng,LI Dahua,et al. Geochemical Features of the Late Permian Tuff in Nanchuan District, Chongqing, Southwestern China[J]. Geological Review,2014,60(1):169−177. doi: 10.16509/j.georeview.2014.01.019
[27] 陈全红,李文厚,胡孝林,等. 鄂尔多斯盆地晚古生代沉积岩源区构造背景及物源分析[J]. 地质学报,2012,86(7):1150−1162. doi: 10.3969/j.issn.0001-5717.2012.07.011 CHEN Quanhong,LI Wenhou,HU Xiaolin,et al. Tectonic setting and provenance analysis of late paleozoic sedimentary rocks in the Ordos Basin[J]. Acta Geologica Sinica,2012,86(7):1150−1162. doi: 10.3969/j.issn.0001-5717.2012.07.011
[28] WINCHESTER J A,FLOYD P A. Geochemical discrimination of different magma series and their differentiation products using immobile elements[J]. Chemical Geology,1977,20:325−343. doi: 10.1016/0009-2541(77)90057-2
[29] 曹代勇,徐 浩,刘 亢,等. 鄂尔多斯盆地西缘煤田构造演化及其控制因素[J]. 地质科学,2015,50(2):410−427. doi: 10.3969/j.issn.0563-5020.2015.02.005 CAO Daiyong,XU Hao,LIU Kang,et al. Coalfield tectonic evolution and its controlling factors at the western margin of Ordos Basin[J]. Chinese Journal of Geology,2015,50(2):410−427. doi: 10.3969/j.issn.0563-5020.2015.02.005
[30] 李 俊, 赵红格, 汪 建, 等. 鄂尔多斯盆地西缘中部三叠系延长组碎屑岩沉积环境及物源示踪[J/OL]. 沉积学报, 1−23[2023-02-14]. https://doi.org/10.14027/j.issn.1000-0550.2022.120. LI Jun, ZHAO Hongge, WANG Jian, et al. Sedimentary Environment and Provenance Tracing of Clastic Rocks from the Triassic Yanchang Formation in the Western Margin of the Ordos Basin[J/OL]. Acta Sedimentologica Sinica, 1−23[2023-02-14]
计量
- 文章访问数: 117
- HTML全文浏览量: 12
- PDF下载量: 43
