Study on deformation characteristics of room-and-pillar mined-out area based on physical simulation test
-
摘要:
房柱式采空区具有采出率较小、留设煤柱较大的特点，其覆岩移动规律及地表变形规律与长壁式采空区差别较大，在煤柱失效情况下，地表会产生较大变形，威胁地面建构筑物的安全。为研究浅埋缓倾房柱式采空区覆岩和松散层的变形发育机制，基于PIV图像处理技术和模块化组装思路，在室内建立物理模型开展物理模拟试验，反演了房柱式采空区覆岩和松散层的变形发育过程，分析了房柱式开采和煤柱失稳条件下覆岩和松散层的变形机制。研究结果表明：① 房柱式开采条件下，浅埋缓倾采空区覆岩变形可划分为3个区域：“条带状”变形区域、“梯形状”变形区域和“倒漏斗型”变形区域。松散层中的变形呈“梯形状”发育。② 煤柱失稳条件下，覆岩变形规律保持不变，松散层变形特征由“梯形状”发育为“漏斗型”，且随着失稳煤柱的增加，松散层变形增大，“漏斗型”变形区域由中心向两端逐渐扩展。③ 房柱式开采结束后，采空区顶板受力可简化为承受均布荷载的简支梁，在煤柱持续失稳的情况下，覆岩和松散层的静载将导致采空区顶板达到极限平衡状态，顶板破坏模式为拉破断，在拉应力作用下产生裂缝，但顶板并未完全断裂，采空区顶板形成一种带裂缝的悬臂梁结构。浅埋缓倾房柱式采空区煤柱失稳对地表变形影响剧烈，松散层变形与覆岩变形基本处于同一量级。
Abstract:The room-and-pillar mined-out area has the characteristics of small recovery rate and large coal pillar, and its overburden movement law and surface deformation law are quite different from those of long-wall mined-out area. In the case of coal pillar failure, the surface will produce large deformation, threatening the safety of surface structures. In order to study the deformation and development mechanism of overlying rock and unconsolidated formation in shallow and gently tilted room-and-pillar mined-out area, based on PIV image processing technology and modular assembly idea, a physical model was established in the laboratory to carry out physical simulation experiments, inversion of the deformation and development process of overlying rock and unconsolidated formation in room-and-pillar mined-out area, and analysis of the deformation mechanism of overlying rock and unconsolidated formation under the condition of room-and-pillar mining and coal pillar instability. The research results show that: ① the overburden deformation of shallow buried and gently mined goaf can be divided into three areas: “banded” deformation area, “trapezoidal” deformation area and “inverted funnel” deformation area. The deformation in the unconsolidated formation is “trapezoidal”. ② Under the instability condition of coal pillar, the overburden deformation law remains unchanged, and the deformation characteristics of unconsolidated formation develop from “trapezoidal” to “funnel”. With the increase of unstable coal pillar, the deformation of unconsolidated formation increases, and the “funnel” deformation area gradually expands from the center to both ends. ③ After the end of room-and-pillar mining, the stress of the goaf roof can be simplified to simply supported beams bearing uniform load. In the case of continuous instability of coal pillar, the static load of overlying rock and unconsolidated formation will lead to the goaf roof reaching the ultimate equilibrium state. The failure mode of the roof is tensile breaking, and cracks are generated under the action of tensile stress, but the roof is not completely broken. The goaf roof forms a cantilever beam structure with cracks. The instability of coal pillar in shallow and gently tilting room-and-pillar mined-out area has a strong influence on the surface deformation, and the deformation of unconsolidated formation and overlying rock are in the same order of magnitude.
-
0. 引 言
在20世纪90年代和21世纪初，中国煤炭主产区的村办集体煤矿和私人煤矿多以房柱式采煤法开采煤炭资源，导致中国老矿区遗留了“星罗密布”的房柱式小窑煤矿采空区。房柱式采空区具有采出率较小、留设煤柱较大的特点，其覆岩移动规律及地表变形规律与长壁式采空区差别较大[1-2]。在地下水、煤柱风化等作用下，煤柱失稳引起房柱式采空区顶板失稳塌落，进而引起地表变形对上方建构筑物的稳定性造成影响，甚至引起矿震导致地面的大面积塌陷[3-7]。煤柱失稳曾引发我国房柱式采空区发生多次大面积塌陷，引起的矿震对上方建构筑物和人员生命安全造成了极大的威胁。2005年11月6日，中国发生了53 000 m2的灾难性地面沉降，造成38人受伤，37人死亡。自2009年9月至2015年底，中国某矿区共监测到矿震13次，震级2.1～3.3级[8]。
采空区覆岩的变形机制[9-11]直接决定了地表沉陷的发育机制和发育程度。覆岩厚度的不同导致其破坏机制也存在较大差异，厚岩层的主要破坏机制为压剪破坏，剪切裂纹占主导地位[12-15]，而薄岩层的破坏为基岩的弯曲为小挠度问题，根据煤矿实测数据结合现场技术人员经验，基岩层位移达到0.1 m，就十分危险，而基岩层的平均厚度小于20 m，挠度与厚度之比远小于1/5。根据断裂力学对裂纹扩张的的描述[16-18]，可根据裂纹形成的力学机制将裂纹划分为张开型裂纹、滑移型裂纹和复合型裂纹。鉴于此，分析采空区覆岩运动规律时，确定覆岩的变形机制具有十分重要的现实意义。基岩的沉降变形为上部松散层的运动创造了空间，松散层的抗拉、抗剪及抗弯能力都很低，基岩塑性变形会引起松散层的滑移破坏、拉裂破坏和剪切破坏，在充分采动的作用下，覆盖层最终会形成倒梯形破坏形态[19-20]。
综上所述，目前对房柱式采空区变形失稳的研究主要集中在煤柱失稳方面，而对采空区覆岩和松散层变形机制的研究主要是针对长壁式采空区。房柱式采空区上覆地层的失稳变形机制还未引起足够的重视。鉴于此，笔者以山西晋城某房柱式煤矿采空区为研究对象，探讨房柱式采空区上覆地层的失稳变形机制，以期为类似工程防灾减灾和环境保护提供一定的借鉴。
1. 研究概况
近年来，随着基础设施的持续建设，大量的公路、铁路、建筑物不可避免地位于房柱式采空区之上。某高速公路在建设过程中，路线走廊带内存在大量的房柱式小窑采空区。由于房柱式采空区开采煤层较浅、回采率较小、留设煤柱较大等特点，煤柱的失效会引起小窑采空区的失稳变形，对上方公路的建设和安全运营造成较大的威胁[21]。
研究区位于山西省晋城市的沁水煤田，煤系地层主要二叠系下统山西组（P1s），煤层产状近水平向，在20世纪90年代经历了大面积的房柱式开采，采深约50 m，煤柱留设宽度约7～10 m，回采率在30%～50%，平均采厚5 m。该处地层典型断面从上至下依次为：
1）第四系中更新统离石组（Q2l）黏土，厚度36 m。
2）二叠系上统石千峰组（P2sh）紫红色泥岩，厚度13.5 m。
3）二叠系下统下石盒子组（P1x）粉砂岩，厚度3 m。
4）二叠系下统山西组（P1s）3号煤层，厚度5 m。
5）采空区底板为二叠系下统山西组（P1s）泥岩。
本文以该典型地层断面为基础，开展室内物理模拟试验，探讨房柱式采空区上覆地层的失稳变形机制。
2. 物理模拟试验设计
本次物理模拟试验装置由定制模型箱、模型材料、监测设备组成。
2.1 定制模型箱
选择几何相似比为1∶100的平面应力物理模型。试验所用模型箱内部尺寸为长150 cm、宽10 cm、高100 cm。模型箱采用三角钢焊制边框，侧面采用定制钢化玻璃做表面，除背面外，其余三侧固定。背面下部20 cm高度的钢化玻璃固定，在25 cm处角铁上打孔，做限位螺栓，用于支撑上部活动玻璃。下部20 cm填实作为采空区底板，20～25 cm部分为煤层，如图1所示。
2.2 模型材料
本次模型试验的材料主要分3部分：
1）煤层及采空区：煤层部位采用43块宽3.5 cm、高5 cm的木块模拟煤层，木块长度为20 cm，便于抽取以模拟煤层开采。
2）上覆岩层：选取河沙和重晶石粉为骨料，石膏作为胶结剂，运用相似原理，通过室内试验确定配比。
泥岩：水∶重晶石粉∶河砂∶石膏质量比为 1∶1∶1∶1，制作长×宽×高=20 cm×10 cm×1 cm的块体，对块体进行分层交错码砌，模拟采空区上覆泥岩。
粉砂岩：水∶重晶石粉∶河砂∶石膏质量比为 0.8∶1∶1∶1.2，制作长×宽×高=50 cm×10 cm×1 cm的块体，对块体进行分层对称码砌，模拟采空区顶板粉砂岩。
3）松散层：采用5～10 mm的砂砾模拟松散层。
2.3 监测系统及内容
在模型试验开始前，对模拟岩层的石膏板进行处理，采用3D打印机制作小型滚轮器具，在滚轮上裹上凸颗粒橡胶膜，蘸印油后在板的侧面均匀涂黑点，以便于粒子图像识别，对岩层的变形进行监测。
模型试验过程中，利用专业相机对物理模型的正面进行实时拍照，如图2所示。试验结束后利用PIV配套软件对照片进行后处理，计算同一粒子不同阶段的位移变化，就可以得到松散层和覆岩的位移场变化信息。本次监测测量的空间分辨率为0.02 mm。
2.4 试验方案
根据该处采空区特征，并结合房柱式开采煤柱留设宽度经验，本次物理模拟试验先完成开挖步序1—5，模拟不同回采率采空区；然后依次完成开挖步序6—9，模拟煤柱失稳，具体开挖步序和位置如图3和表1所示，每个开挖步序按木块编号从大至小依次开挖。
表 1 开挖步序与模拟工况对照表Table 1. Excavation sequence and simulation condition comparison table开挖步序 模拟工况 1 采出率30% 2 采出率35% 3 采出率40% 4 采出率45% 5 采出率50% 6 1号煤柱失稳 7 2号煤柱失稳 8 3号煤柱失稳 9 4号煤柱失稳 为保证上覆地层的充分变形，每一个开挖步序结束后停放1 h，保证位移充分，然后拍照记录。若遇上覆地层明显移动，则也停下放置1 h，再次拍照记录。
3. 试验结果分析
3.1 不同回采率房柱式采空区上覆地层变形特征
根据试验结果，得到不同回采率情况下采空区覆岩和松散层的变形特征，如图4所示。
![]()
图 4 不同回采率采空区覆岩和松散层变形特征注：图中不同颜色区域为2条竖向变形等值线数值的中间区域，等值线数值标度单位为mm。Figure 4. Deformation characteristics of overlying rock and unconsolidated formation in mined-out area with different mining rates房柱式开采结束后，由于开采的卸载作用引起覆岩和上覆松散层的变形，但覆岩整体处于稳定状态，未发生垮落，能对上覆松散层起到较好的支撑作用。
覆岩整体发生变形，其主要变形部位集中在采空区中部顶板处，松散层在覆岩移动影响下也发生变形，变形区域发育至地表。当煤层回采率为30%时，基本顶处最大变形为0.55 mm。随着回采率的增大，覆岩整体变形增大，覆岩变形由中部基本顶逐渐向两侧扩展，变形区域逐渐增大，且变形量级也不断增大，松散层变形区域也逐渐增大。
由图4可以看出，采空区中部基本顶覆岩变形呈“条带状”发育，远离变形中心部位的覆岩变形呈“梯形状”发育。由于覆岩较薄，靠近变形中心部位的覆岩变形呈“梯形状”发育，远离变形中心部位的覆岩变形呈“倒漏斗型”发育。但随着回采率的增大，覆岩变形区域由中心向两端扩展，“梯形状”变形区域逐渐转变为“倒漏斗型”。
松散层属于松散介质，其抗拉、抗剪及抗弯能力都很低，覆岩变形引起的松散层变形呈“梯形状”，随着回采率的增加，梯形面积不断增大。由图4可以看出，松散层中的变形呈“梯形状”发育，随着回采率的增大，梯形面积逐渐增大，变形区域不断由中心向两端扩展，当回采率达到50%时，地面沉降变形超过0.2 mm，梯形顶边消失。
采空区覆岩变形可划分为3个区域：“条带状”变形区域、“梯形状”变形区域和“倒漏斗型”变形区域。其中，“条带状”变形区域在整个采空区中部的基本顶发育，变形最为剧烈，对覆岩的“梯形状”变形区域的大小和变形量级起控制作用；“倒漏斗型”变形区域范围最大，且变形量随着距离“条带状”变形区域越远，衰减幅度越大。
3.2 煤柱失稳条件下上覆地层变形特征
由图5煤柱失稳条件下采空区覆岩及松散层变形特征可知，煤柱失稳条件下，松散层变形特征由“梯形状”发育为“漏斗型”，随着失稳煤柱的增加，松散层变形增大，“漏斗型”变形区域由中心向两端逐渐扩展。覆岩变形规律与房柱式开采煤柱未失稳条件下相似，但“条带状”变形区域、“梯形状”变形区域和“倒漏斗型”变形区域变形量显著增大。
![]()
图 5 煤柱失稳条件下采空区覆岩及松散层变形特征Figure 5. Deformation characteristics of overlying rock and unconsolidated formation in mined-out area under instability of coal pillar在房柱式采煤结束后，顶板已经发生弯曲变形，最大变形位置为模型中部煤层基本顶处，基本顶最大下沉量为0.65 mm。在1号煤柱失稳的情况下，基本顶变形继续增大，最大下沉量达到1.5 mm，覆岩开始产生离层裂缝。随着煤柱的持续失稳，基本顶变形继续增大，离层裂缝持续发育，基本顶中部出现裂缝。
在2号煤柱失稳后，基本顶形成“悬臂梁”支撑结构，稳定的支撑住了上覆岩体荷载，而覆岩整体形成拱形稳定结构；随着3号煤柱和4号煤柱的失稳，顶板和覆岩虽然下沉量持续增大，但“悬臂梁”结构保持稳定，拱形结构由下至上逐渐扩展，拱脚位置由中间逐渐向两边扩展，“悬臂梁”结构支撑住了上部岩体的重量，未发生垮塌。
4. 分析与讨论
4.1 松散层覆岩联动破坏过程分析
根据本次物理模拟试验的结果，将浅埋缓倾房柱式采空区变形失稳条件下的松散层覆岩破坏联动演化过程描述为：① 煤层开采前，松散层自重产生的静载全部作用在基岩上，处于静力平衡状态；② 房柱式开采结束后，原先的静力平衡被打破，在采空区中部的基本顶处发生应力集中，采空区覆岩变形由应力集中区域向外可划分为“条带状”、“梯形状”和“倒漏斗型”3个区域，并引起上覆松散层的同步协调沉降，松散层中的变形呈“梯形状”发育。③ 在煤柱失稳的条件下，松散层变形由“梯形状”发育为“漏斗型”，由于松散层自身性质无法形成稳定承载结构，其运移区域内的大范围介质仍将以载荷形式作用在基岩上，从而进一步加剧了采空区覆岩变形的发展；④ 覆岩的加剧变形，反过来又致使松散层进一步沉降运移，“漏斗型”破坏区域进一步扩大，达到新的平衡状态；⑤ 煤柱持续失稳，直至顶板逐渐变成一种带裂缝的悬臂梁结构，达到极限平衡状态。
综上可知，松散层引起的静载是覆岩变形的根本动力，而松散层与覆岩力学特性及两者的联合演化是形成其特殊变形机制的根本原因。
4.2 房柱式采空区顶板岩层破断机制分析
浅埋缓倾煤层房柱式开采结束后，顶板受力可简化为承受均布荷载q的简支梁，如图6所示，设该段采空区顶板的长度为l，则均布荷载的合力为$ F_{\mathrm{R}}=q l $，均布荷载q作用下简支梁的挠曲轴方程为：
$$ w=\frac{q x}{24 E I}\left(2 l x^2-x^3-l^3\right) $$ (1) 式中：EI为梁的刚度。
由式（1）可知，简支梁跨度中点截面C的挠度最大，其值为：
$$ \delta_{\max }=\frac{5 q l^4}{384 E I} $$ (2) 在1号煤柱失稳情况下，离层裂隙开始发育，顶板跨度翻倍，$ \delta_{\max } $与$ l^{4} $呈正比，$ \delta_{\max } $增大16倍，顶板挠度最大值出现在煤柱1上方。
随着2号煤柱失稳，顶板断裂，可认为上覆地层荷载由两段悬臂梁支撑。由于离层裂隙增大，悬臂梁部分与上覆岩层脱离，其受力可以简化为在距离支点a处存在一个集中荷载，如图7所示，该种情况下，悬臂梁的挠曲轴方程为：
$$ w=\frac{F x^2}{6 E I}(x-3 a),(0 \leqslant x \leqslant a) $$ (3) $$ w=\frac{F a^2}{6 E I}(a-3 x),(a \leqslant x \leqslant l) $$ (4) 在上覆地层荷载作用下，挠度最大处为B点（图7），即变形主要发生在断裂处：
$$ w_B=-\frac{F a^2}{6 E I}(3 l-a) $$ (5) 此时，顶板达到极限平衡状态。
5. 结 论
1）浅埋缓倾煤层在房柱式开采条件下，采空区覆岩变形可划分为3个区域：“条带状”变形区域、“梯形状”变形区域和“倒漏斗型”变形区域。其中，“条带状”变形区域在整个采空区中部顶板处发育，变形最为剧烈，对覆岩的“梯形状”变形区域的大小和变形量级起控制作用；“倒漏斗型”变形区域范围最大，且变形量随着距离“条带状”变形区域越远，衰减幅度越大。随着回采率的增大，覆岩变形区域由中心向两端扩展，“梯形状”变形区域逐渐转变为“倒漏斗型”。松散层中的变形呈“梯形状”发育，随着回采率的增大，梯形面积逐渐增大，变形区域不断由中心向两端扩展。
2）煤柱失稳条件下，覆岩变形规律保持不变，松散层变形特征由“梯形状”发育为“漏斗型”，且随着失稳煤柱的增加，松散层变形增大，“漏斗型”变形区域由中心向两端逐渐扩展。
3）房柱式开采结束后，采空区顶板受力可简化为承受均布荷载的简支梁，在煤柱持续失稳的情况下，覆岩和松散层的静载将导致顶板达到极限平衡状态，顶板破坏模式为拉破断，在拉应力作用下产生裂缝，但顶板并未完全断裂，采空区顶板形成一种带裂缝的悬臂梁结构。
-
![]()
图 4 不同回采率采空区覆岩和松散层变形特征
注：图中不同颜色区域为2条竖向变形等值线数值的中间区域，等值线数值标度单位为mm。
Figure 4. Deformation characteristics of overlying rock and unconsolidated formation in mined-out area with different mining rates
![]()
图 5 煤柱失稳条件下采空区覆岩及松散层变形特征
Figure 5. Deformation characteristics of overlying rock and unconsolidated formation in mined-out area under instability of coal pillar
表 1 开挖步序与模拟工况对照表
Table 1 Excavation sequence and simulation condition comparison table
开挖步序 模拟工况 1 采出率30% 2 采出率35% 3 采出率40% 4 采出率45% 5 采出率50% 6 1号煤柱失稳 7 2号煤柱失稳 8 3号煤柱失稳 9 4号煤柱失稳 
下载: 导出CSV
-
[1] SMITH G J,ROSENBAUM M S. Recent underground investigations of abandoned chalk mine workings beneath Norwich City,Norfolk[J]. Engineering Geology,1993,36(1-2):67−78. doi: 10.1016/0013-7952(93)90019-9
[2] MILLER R D,STEEPLES D W,SCHULTE L. Shallow seismic reflection study of salt dissolution well field near Hutchinson,KS[J]. International Journal of Rock Mechanics and Mining Sciences & Geomechanics Abstracts,1994,31(2):1−21.
[3] WANG X R,YANG T H,GUAN K,et al. Stability Evaluation of multi-pillar and roof system based on instability theory[J]. Rock Mechanics and Rock Engineering,2022,55:1461−1480. doi: 10.1007/s00603-021-02712-8
[4] LUO R,LI G Y,CHEN L,et al. Ground subsidence induced by pillar deterioration in abandoned mine districts[J]. Journal Of Central South University,2020,27:2160−2172. doi: 10.1007/s11771-020-4438-3
[5] HENRY S. HASEGAWA,ROBERT J. WETMILLER and DON. J. GENDZWILL. Induced Seismicity in Mines in Canada--An Overview[J]. PAGEOPH,1989,129(3/4):423-453.
[6] BELL F G,BRUYN I A D. Subsidence problems due to abandoned pillar workings in coal seams[J]. Bulletin of Engineering Geology and the Environment,1999,57:225−237 doi: 10.1007/s100640050040
[7] ROGER J M,YOGINDER P C,THOMAS R. Subsidence prediction in shallow room and pillar mines[J]. International Journal of Mining and Geological Engineering,1986,4:39−46 doi: 10.1007/BF01553755
[8] 朱卫兵,许家林,陈璐,等. 浅埋近距离煤层开采房式煤柱群动态失稳致灾机制[J]. 煤炭学报,2019,44(2):358−366. ZHU Weibing,XU Jialin,CHEN Lu,et al. Mechanism of disaster induced by dynamic instability of coal pillar group in room-and-pillar mining of shallow and close coal seams[J]. Journal of China Coal Society,2019,44(2):358−366.
[9] DOU L M,HE X Q,HE H,et al. Spatial structure evolution of overlying strata and inducing mechanism of rockburst in coal mine[J]. Transactions of Nonferrous Metals Society of China,2014,24(4):1255−1261. doi: 10.1016/S1003-6326(14)63187-3
[10] CHLEBOWSKI D,ZORYCHTA A,BURTAN Z,et al. Location of opening-out cross-cuts in relation to rockburst hazard conditions along the face in the light of model testing[J]. Archives of Mining Sciences,2017,62(1):121−130. doi: 10.1515/amsc-2017-0009
[11] JU Y,WANG Y,SU C,et al. Numerical analysis of the dynamic evolution of mining-induced stresses and fractures in multilayered rock strata using continuum-based discrete element methods[J]. International Journal of Rock Mechanics and Mining Sciences,2019,113:191−210. doi: 10.1016/j.ijrmms.2018.11.014
[12] CHEN J,ZHOU L,XIA B,et al. Numerical Investigation of 3D distribution of mining-induced fractures in response to longwall mining[J]. Natural Resources Research,2021,30(1):889−916. doi: 10.1007/s11053-020-09759-4
[13] GAO F,STEAD D,COGGAN J. Evaluation of coal longwall caving characteristics using an innovative UDEC Trigon approach[J]. Computers and Geotechnics,2014,55:448−460. doi: 10.1016/j.compgeo.2013.09.020
[14] MAJDI A,HASSANI F P,NASIRI M Y. Prediction of the height of destressed zone above the mined panel roof in longwall coal mining[J]. International Journal of Coal Geology,2012,98:62−72. doi: 10.1016/j.coal.2012.04.005
[15] KELLY M,LUO X,CRAIG S. Integrating tools for longwall geomechanics assessment[J]. International Journal of Rock Mechanics and Mining Sciences,2002,39(5):661−676. doi: 10.1016/S1365-1609(02)00063-1
[16] LEI B,ZUO J,LIU H,et al. Experimental and numerical investigation on shale fracture behavior with different bedding properties[J]. Engineering Fracture Mechanics,2021,247:107639. doi: 10.1016/j.engfracmech.2021.107639
[17] ZUO J,SHEN J. The Hoek-Brown Failure criterion-From theory to application[M]. Singapore:Springer,2020.
[18] JI S,WANG Z,KARLOVšEK J. Analytical study of subcritical crack growth under mode I loading to estimate the roof durability in underground excavation[J]. International Journal of Mining Science and Technology,2022,32(2):1−15.
[19] 左建平,孙运江,钱鸣高. 厚松散层覆岩移动机理及“类双曲线”模型[J]. 煤炭学报,2017,42(6):1372−1379. ZUO Jianping,SUN Yunjiang,QIAN Minggao. Movement mechanism and analogous hyperbola model of overlying strata with thick alluvium[J]. Journal of China Coal Society,2017,42(6):1372−1379.
[20] 张广超,陶广哲,孟祥军,等. 巨厚松散层下软弱覆岩破坏规律[J]. 煤炭学报,2022,47(11):3998−4010. ZHANG Guangchao,TAO Guangzhe,MENG Xiangjun,et al. Failure law of weak overburden stratum underlying extra-thick alluvium[J]. Journal of China Coal Society,2022,47(11):3998−4010.
[21] 曹树刚,勾攀峰,樊克恭. 采煤学[M]. 北京:煤炭工业出版社,2017. CAO Shugang,GOU Panfeng,FAN Kegong. Coal mining[M]. Beijing:China Coal Industry Press,2017.
计量
- 文章访问数: 26
- HTML全文浏览量: 2
- PDF下载量: 26
