Open-pit mine slope optimization with the morphology and spatial support effect
-
摘要:
为探究露天煤矿到界边坡的坡面形态优化设计方法，以国家能源集团神华新疆黑山露天矿南帮边坡为工程实例，应用简化Bishop法，基于边坡断面形态效应理论，设计了边坡在二维视角下不同凹凸程度的坡面形态并进行了稳定性评价与对比，得到较优设计；基于强度折减法，应用数值模拟方法仿真了较优坡面形态在边坡三维支挡作用下的形变过程，通过分析形变特性，进一步提高边坡坡面的凹凸程度得到了最优坡面形态。结果表明：二维设计条件下，受局部边坡稳定性的限制，在黑山露天矿南帮边坡+2 385 m水平以下自底周界向上依次布置两组组合台阶，可得凹凸程度较优的坡面形态，坡面的整体稳定性系数为1.344，局部稳定性系数为1.204。三维设计条件下，较优形态在内排土场与横采工作帮的三维支挡作用影响下，边坡局部与整体的稳定性均得到改善，存在进一步形态优化的空间。可在较优形态基础上，通过增设一组组合台阶的方式局部加陡+2 430 m水平以下边坡区段以取得最优坡面形态。最优形态在数值模拟边坡形变过程中，进入临界失稳前变速阶段后呈现的形变速度相对缓慢，模拟结果稳定，最终强度折减系数为1.20。最优形态在地表界与底周界不变的前提下，保障黑山露天矿边坡安全的同时可节省首采区岩石剥离量约384万m3。
Abstract:In order to explore the optimal design method of slope morphology of open-pit coal mine, the south slope of Heishan Open-pit Mine was taken as an engineering example. Based on the simplified Bishop method and the theory of slope morphology effect, the slope with different concave and convex degrees under the two-dimensional perspective was designed, and the stability evaluations are carried out to obtain the better slope morphology. Based on the strength reduction method, the numerical simulation method is used to simulate the deformation process of the better slope morphology under the action of three-dimensional retaining. By analyzing the deformation characteristics, the optimal slope morphology is obtained by further improving the concave and convex degree of the slope surface. The results show that the slope morphology with better concave and convex degree under the two-dimensional design conditions is obtained. Due to the limitation of local slope stability, only two groups of composite steps can be arranged from the bottom to the top below the level of + 2385 m in the southern slope. The overall stability coefficient of the slope is 1.344, and the local stability coefficient is 1.204. Under the condition of three-dimensional design, the local and overall stability of the slope is improved under the influence of the support from the internal dump and the horizontal mining working side, and there is further morphological optimization space. On the basis of the better shape, the optimal slope shape can be obtained by adding a group of combined steps to local steepen the slope section below + 2430 m level. In the process of numerical simulation , the deformation speed presented is relatively slow after entering the variable speed stage before critical instabilityand the final strength reduction coefficient is 1.20. Under the premise that the surface boundary and bottom boundary are unchanged, the optimal slope morphology can ensure the slope safety and save about 3.84 million cubic meters of rock stripping in the first mining area.
-
0. 引 言
煤炭截割产尘、煤与瓦斯突出等矿山灾害问题时刻困扰着煤矿的安全高效生产[1-2]。煤矿开采会产生大量粉尘[3]，由粉尘引起的尘肺病是中国最严重的职业病，约占职业病病例总数的80%[4-5]。此外，煤与瓦斯突出时，井下工作面和巷道内的气流会受到严重破坏，并可能引发瓦斯爆炸等次生灾害[6]。煤层注水作为一种抑制煤与瓦斯突出[7]、降低煤尘生产量的有效方法，它可以湿润煤体内的原生煤尘，降低煤尘的产生量，抑制煤层瓦斯的解析。但是，注水润湿、瓦斯渗流与煤岩内部孔裂隙的发育情况有直接关联[8]，煤是一种多孔介质，孔径跨度范围大且孔隙形态多样、结构复杂[9]，这些孔隙裂隙所构成的空间是煤层中唯一存储和输运流体的空间[10]，因此，煤的孔隙裂隙结构一直是煤矿安全领域研究的重点。煤是一种应力敏感的有机岩石，在地应力作用下，煤内部裂隙、孔隙等细观结构将发生变化，而变化的孔隙通道对煤体注水、润湿有着至关重要的作用[11]，煤矿中煤柱处于单轴压缩状态，当煤柱内裂隙发育后漏风，漏粉导致采空区自燃发火，因此，探明在单轴压缩作用下煤体裂隙孔隙细观结构演化过程，掌握煤体细观结构改变的空间展布形态，对认清煤体破坏过程中的裂纹扩展及其分布，具有理论研究价值和实际工程意义。
为了更好理解煤岩内部细观结构，以往的研究人员不仅分析了煤岩变形破坏过程中微观结构对破坏形式产生的影响以及渗透率的变化[12-13]，而且在单轴压缩条件下研究了裂纹的微损伤和细观损伤的演化特征以及微裂隙初始密度与煤岩单轴抗压强度之间的联系[14-15]。人们对于煤体内部微结构的研究也越来越深入[16-19]。随着CT、核磁共振、扫描电镜等细观成像技术的发展，通过细观尺度的研究来揭示宏观结构特征演化的内在原因，尤其CT技术应用较为广泛。CT扫描技术利用X射线对岩石进行全方位快速且无损扫描并重构三维结构特征，能呈现样品中孔裂隙空间分布与连通性，精细定量表征煤岩孔裂隙形态、分布以及矿物分布，因此该方法被广泛应用于表征煤岩微米-毫米级连通微裂隙特征[20-21]。目前，诸多学者应用CT扫描技术对煤岩的精细结构进行系统评价。研究主要包括：评价特定煤岩的三维孔裂隙空间分布特征、动态监测煤岩的孔裂隙系统随外界条件的变化趋势、及探究不同图像处理技术在煤岩评价上的优缺点等[22-24]。现代运用CT对煤体结构的分析主要是通过CT扫描观测吸附瓦斯煤体细观结构的受载破坏过程，分析孔隙结构对煤体受载变形的影响，结合分型理论对煤岩细观结构进行分析[25-31]。研究人员运用CT可以定量表征煤岩孔隙度、孔喉特征、配位数、连通性、通道长度、喉道尺寸和孔隙体积的测量分布之间的相关性等参数[32-34]。以往的CT扫描都是状态的扫描，不能实时监测样品的变化情况，原位加载技术使样品可以在保持负载的情况下进行扫描，提高扫描样品时孔裂隙的精确度[35]。原位CT可以在加载样本时逐层捕获详细的损坏状态并扫描测量裂纹[36]，其对于研究煤岩细观结构有着重要意义，通过它能够获得样品内部细观结构随着应力加载的动态演化规律，进而定量、定性揭示其内部裂隙孔隙时空展布特征及裂隙的萌生、扩展、贯通与承载能力的关联[37]。
尽管CT扫描技术对前人探究煤岩微观结构奠定了基础，但仅运用CT并不能全面动态地研究煤岩微观结构随着轴向压力变化的规律，而原位加载系统则为实时监测煤样随着轴向压力变化的规律特点提供了支持，但利用CT扫描技术结合原位加载系统对煤岩多尺度孔裂隙评价的研究相对较少。本文基于原位加载系统并结合NanoVoxel-3502EX射线三维显微镜对煤样进行单轴压缩，通过AVIZO进行煤样三维重构，探究煤岩细观结构的时空演化特征。
1. 试验及煤样
1.1 试验装置
采用NanoVoxel-3502E X射线三维显微镜和原位加载系统进行试验，X射线三维扫描成像原理如图1所示。原位加载系统载荷力范围是0～5 kN。NanoVoxel-3502E X射线三维显微镜是具有超高分辨率的无损伤三维显微成像设备。X射线三维显微镜，采用独特的X光光学显微成像技术，利用不同角度的X射线透视图像，结合三维数字重构技术，提供样品内部复杂结构的高分辨率三维数字图像，对样品内部的微观结构进行微米尺度上的数字化三维表征，以及对构成样品的物质属性进行分析。
![]()
图 1 X射线三维扫描成像原理及二维切片展示Figure 1. X-ray three-dimensional scanning imaging principle and two-dimensional section presentation1.2 煤 样
本试验煤样取自中国山东省济宁市兴隆庄煤矿，属于烟煤，经工业分析表明，煤样的灰分9.4%，水分2.8%，挥发分28.4%，固定碳59.4%(图2)。为便于扫描将煤样试件磨制成长方体，试件尺寸为：6 mm×6 mm×12 mm，用岩石断面打磨机打磨煤样试件端面，使端面不平整度误差不大于0.02 mm。扫描前用热缩管包裹煤样试件，防止煤样试件破坏后崩落，影响后期试样的CT扫描。
1.3 试验过程
将煤样试件固定于样品台，通过NanoVoxel-3502E X射线三维显微镜和原位加载系统进行煤样试件的扫描。为使煤样与加载台接触，初期原位加载台以1.0 mm/min的速度进行移动，并设定预加轴向力为20 N。单轴压缩实验过程中采用轴向位移控制方式进行加载，位移加载速度变为0.1 mm/min，轴向位移每加载0.1 mm进行一次CT扫描，共计进行了5次扫描。试验过程参数见表1。
表 1 位移-轴向力关系与三维扫描参数Table 1. Displacement-axial force telationship and 3D scanning parameters加载速率/
(mm·min−1)轴向位移/
mm轴向力/
N扫描电流/
μA扫描电压/
kV1.0 0.146 20 100 130 0.1 0.246 431 0.346 732 0.446 1100 0.546 492 试验路径如图3所示。获得数据后结合AVIZO将煤样试件三维重构，其中图像尺寸为1200 pix×1300 pix×1300 pix、像素分辨率为8 μm。将煤样外层包裹的热缩管裁切后得到816个切片。
2. CT扫描数据处理情况
将五个状态下的孔裂隙进行骨架化，从图4中可以清楚地看出孔裂隙（包括贯通裂隙）随时间和应力增长而逐渐发育成熟的过程：20 N应力状态下微小孔裂隙出现；431 N下孤立孔裂隙逐渐密集，发育出贯通裂隙，并且垂直X-Y平面；732 N下的孔裂隙相互勾通，局部孔裂隙被压实，贯通煤样试件的裂隙逐渐发育，裂隙网络开始显现；峰值1100 N轴向应力的作用下使贯通煤样试件的裂隙发育地更加成熟；最终在峰后492 N下可以看出残余荷载使微小孔裂隙部分消失。基于孔裂隙骨架化可以获得孔裂隙配位数的情况。
3. 煤岩细观情况分析
煤样试件置于NanoVoxel-3502E X射线三维显微镜设备中在各状态下进行扫描，将所得数据运用Avizo处理后得到各状态下的球形度、配位数、曲率、开度一系列参数。
3.1 煤岩球形度情况分析
将处理所得球形度数据进行统计分析得到图5。球形度反映了煤岩内部孔裂隙的形状，形貌上越接近球的颗粒，其球形度越接近于1，球形度越偏离1，则其孔裂隙形状越不规则，越狭长。
通过图5可知在时间和应力的作用下，孔隙由20～431 N先逐渐增多至18万个，431～732 N逐渐减少至7万个，732～492 N（峰后）再次逐渐增多至11.5万个。煤样孔隙在单轴压缩过程初期先发育增多，而孔裂隙在732 N时开始被压实压密，导致孔裂隙减少9万个，但随着单轴压缩过程的推进，孔隙逐渐增长，部分孔裂隙虽然在732 N被压实压密，但由于轴向力的增大以及残余荷载的作用，孔裂隙的增多是必然趋势。20 N孔隙球形度在0.5～2.5、2.9范围内，431N孔隙球形度在0.5～2.5、2.9范围内，集中分布在1.0～2.5和2.9，球形度分布范围没有明显变化，732～492 N（峰后）孔隙球形度集中在1.5～2.5、2.9范围内，球形度集中分布范围此时开始变窄，部分孔隙被融合并入裂隙网络中，使球形度开始明显发生变化。煤样在单轴压缩过程初期开始中出现微小孔裂隙，从431 N开始孔隙球形度集中分布范围变窄，部分孔裂隙逐渐融合进入裂隙网络，部分孔裂隙被压实压密，峰值和峰后仍然存在微小孔裂隙，但球形度更多的是集中在1.5～2.0。
3.2 煤岩配位数情况分析
对煤岩内部孔裂隙进行骨架化后得到配位数情况，配位数是反映煤岩内部孔裂隙连通程度的重要参数，配位数越大，储层性质越好。
随着轴向力增加，节点与配位数之间的变化情况如图6所示。在20 N轴向力下，节点有1933个，配位数最大值是20，配位数集中分布范围是0～5；施加至431 N的轴向力作用后，节点增加至4.5万个，配位数最高是98，配位数集中分布范围是0～10；当轴向力由431 N增至732 N，节点个数减少了2.6万个，配位数最高是204，配位数集中分布范围是0～10；当轴向力由732 N增至1100 N（峰值），节点增加了4.9万个，配位数最高是405，配位数集中分布范围是0～15；当轴向力由1100 N（峰值）降至492 N（峰后），此时节点数为8.3万个，配位数最高值为319，配位数集中分布范围为0～15。
配位数分布范围随着轴向力的增大而逐渐增大，在20～431 N时存在一部分节点配位数高于集中分布范围，而自732 N后配位数只有1或2个高于集中分布范围，即732 N后配位数集中率更高，同时也说明煤样内部连通性更好。
3.3 煤岩裂隙曲率和开度情况分析
为进一步分析煤样内部裂隙发生演化规律，运用AVIZO对产生贯通裂隙的四个状态（431 N、732 N、1100 N、492 N）进行网格划分，如图7所示，展示了431 N的贯通裂隙网格划分。网格划分完成后得到4个状态的贯通裂隙的曲率和开度的数据，如图8、图9所示。
曲率在一定程度上控制了裂缝发育的密度、方向、宽度和深度。曲率事故几何体不平坦程度的一种衡量，地层曲率用来表征地层受构造应力挤压后的弯曲程度，曲率越大说明岩石所受的应力越大，地层的弯曲程度越大，裂缝裂隙也就越发育。图8中看出4个状态下贯通裂隙的曲率在裂隙网络的边缘普遍较高，即裂隙网络边缘为裂隙发育带，其中732 N的曲率红色范围最多，1100 N（峰值）和492 N（峰后）的曲率绿色范围居多，说明732 N是裂隙网络发育的巅峰，1100 N（峰值）和492 N（峰后）裂隙网络的发育已接近饱和。四个状态下的开度展示如图9所示。在整个发育过程中，431 N和492 N（峰后）的裂隙开度较小，只有小区域的裂隙开度较大，732 N的裂隙逐渐发育，当1100 N（峰值）时，裂隙开度开始出现较大的区域。
4. 结 论
1）煤样在单轴压缩过程初期开始中出现微小孔裂隙，孔隙个数增多，从431 N开始孔隙球形度集中分布范围变窄，部分孔裂隙逐渐融合进入裂隙网络，部分孔裂隙被压实压密，732 N时部分孔隙被压实压密，此时孔隙个数减少9万个，但随着单轴压缩过程的推进，孔隙个数仍然逐渐增长，部分孔裂隙虽然在732 N被压实压密，但由于轴向力的增大以及残余荷载的作用，孔裂隙的增多是必然趋势，峰值和峰后仍然存在微小孔裂隙，但球形度更多的是集中在1.5～2.0。
2）配位数分布范围随着轴向力的增大而逐渐增大，在20～431 N时存在一部分节点配位数高于集中分布范围，而自732 N后配位数只有1或2个高于集中分布范围，即732 N后配位数集中率更高，同时也说明煤样内部连通性更好。
3）4个状态下贯通裂隙的曲率在裂隙网络的边缘普遍较高，即裂隙网络边缘为裂隙发育带，其中732 N的曲率红色范围最多，裂隙开度逐渐变大，732 N是裂隙网络发育的巅峰。1100 N（峰值）和492 N（峰后）的曲率绿色范围居多，1100 N（峰值）时，裂隙开度开始出现较大的区域，492 N（峰后）的裂隙较窄。
-
表 1 岩土体物理力学指标
Table 1 Geotechnical physical mechanics index
岩层 黏聚力/kPa 摩擦角/(°) 容重/(kN•m−3) 弹性模量/GPa 泊松比 中砂岩 88 38.9 26.1 3.38 0.29 煤 45 29 14.1 0.8 0.33 粉砂岩 82 40 26.6 5.48 0.28 
下载: 导出CSV
-
[1] 方玉树. 边坡稳定性分析条分法最小解研究[J]. 岩土工程学报,2008,30(3):331−335. FANG Yushu. The lowest solution of slice method for slope stability analysis[J]. Chinese Journal of Geotechnical Engineering,2008,30(3):331−335.
[2] 张福海,梁 喆,陈 雷. 土质边坡稳定性分析向心条分法[J]. 防灾减灾工程学报,2016,36(4):544−550. ZHANG Fuhua,LIANG Zhe,CHEN Lei. Arc-to-center slice method for soil slope stability analysis[J]. Journal of Disaster Prevention and Mitigation,2016,36(4):544−550.
[3] 王 振,叶晓明,刘永新. 考虑滑坡渐进破坏的改进简布条分法[J]. 岩土力学,2018,39(2):675−682. WANG Zhen,YE Xiaoming,LIU Yongxin. Improved Janbu slices method considering progressive destruction in landslide[J]. Rock and Soil Mechanics,2018,39(2):675−682.
[4] 张晓平,吴顺川,张志增,等. 含软弱夹层土样变形破坏过程细观数值模拟及分析[J]. 岩土力学,2008,29(5):1200−1205. doi: 10.3969/j.issn.1000-7598.2008.05.010 ZHANG Xiaoping,WU Shunchuan,ZHANG Zhizeng,et al. Numerical simulation and analysis of failure process of soil with weak intercalated layer[J]. Rock and Soil Mechanics,2008,29(5):1200−1205. doi: 10.3969/j.issn.1000-7598.2008.05.010
[5] XU Tao,XU Qiang,MAO Lindeng,et al. A numerical analysis of rock creep-induced slide: A case study from Jiweishan Mountain, China[J]. Environmental Earth Sciences,2014,72(6):2111−2128. doi: 10.1007/s12665-014-3119-7
[6] HOEK E, BRAY J W, 卢世宗. 岩石边坡工程[M]. 北京: 冶金出版社, 1983: 223−232. HOEK E, BRAY J W, LU Shizong. Rock slope engineering [M]. Beijing: Metallurgical Press, 1983: 223−232.
[7] 卢坤林,朱大勇. 坡面形态对边坡稳定性影响的理论与试验研究[J]. 岩石力学与工程学报,2014,33(1):35−42. LU Kunlin,ZHU Dayong. Theoretical and experimental study of effects of slope topography on its stability[J]. Chinese Journal of Rock Mechanics and Engineering,2014,33(1):35−42.
[8] 马 力,李克民,韩 流,等. 平直与水平凹形边坡力学结构及稳定性对比研究[J]. 采矿与安全工程学报,2016,33(4):728−733. MA Li,LI Kemin,HAN Liu,et al. Comparative study on stability and mechanical structure of straight slope and horizontal concave slope[J]. Journal of Mining & Safety Engineering,2016,33(4):728−733.
[9] 王 东,张晓敏,曹兰柱,等. 凹边坡稳定性三维极限平衡法研究与应用[J]. 中国安全科学学报,2016,26(7):91−95. WANG Dong,ZHANG Xiaomin,CAO Lanzhu,et al. 3D limit-equilibrium method for calculating concave slope stability and its application[J]. China Safety Science Journal,2016,26(7):91−95.
[10] 田 宇,张 波,刘博文,等. 露天矿多软弱夹层边坡变形机理与形态优化研究[J]. 煤炭科学技术,2021,49(6):189−194. TIAN Yu,ZHANG Bo,LIU Bowen,et al. Study on deformation mechanism and shape optimization of multi-weak interlayer slope in open-pit mine[J]. Coal Science and Technology,2021,49(6):189−194.
[11] 佘长超,陈 涛,韩 流,等. 冻融循环下排土场软岩边坡抗剪特性及时效性研究[J]. 煤炭科学技术,2021,49(10):50−57. SHE Changchao,CHEN Tao,HAN Liu,et al. Research on shear characteristics and timeliness of soft rock slope of dump under freezing-thawing cycle[J]. Coal Science and Technology,2021,49(10):50−57.
[12] 孙 江,田 光,韩 流,等. 内排阻水对非工作帮边坡稳定性影响规律研究[J]. 煤炭科学技术,2021,49(12):186−191. doi: 10.3969/j.issn.0253-2336.2021.12.mtkxjs202112023 SUN Jiang,TIAN Guang,HAN Liu,et al. Study on influence of water blocking of internal dumping site on stability of non-working slope[J]. Coal Science and Technology,2021,49(12):186−191. doi: 10.3969/j.issn.0253-2336.2021.12.mtkxjs202112023
[13] 王 珍,曹兰柱. 斜交断层影响下露天煤矿顺倾层状边坡稳定性研究[J]. 煤炭科学技术,2020,48(5):69−74. WANG Zhen,CAO Lanzhu. Study on dip bedded slope stability of open - pit coal mine under influence of oblique fault[J]. Coal Science and Technology,2020,48(5):69−74.
[14] 王 东,张 禹,曹兰柱,等. 露天矿边坡稳定性的断面形态效应研究[J]. 中国安全科学学报,2018,28(1):105−111. WANG Dong,ZHANG Yu,CAO Lanzhu. et al. Study of section shape effect on stability of open-pit slope[J]. China Safety Science Journal,2018,28(1):105−111.
[15] Alan W BISHOP. The use of the slip circle in the stability analysis of earth slopes[J]. Geotechnique,1955,5(1):7−17. doi: 10.1680/geot.1955.5.1.7
[16] 李广贺,王 东,曹兰柱. 横采内排追踪压帮高度对端帮稳定性影响研究[J]. 安全与环境学报,2019,19(1):191−195. LI Guanghe WANG Dong,CAO Lanzhu. Influence of the covering height on the stability of the end slope in the condition of the criss-cross-way exploitation with inner dumping[J]. Journal of Safety and Environment,2019,19(1):191−195.
[17] 张 禹,吕文伟. 含多层级软弱互层的露天矿边坡稳定性优化[J]. 露天采矿技术,2020,35(6):80−83. ZHANG Yu,LYU Wenwei. Optimization of slope stability of open-pit mine with multi-level weak interlayer[J]. Opencast Mining Technology,2020,35(6):80−83.
[18] 赵红泽,邵志奔,杜海瑞,等. 露天矿背斜区内排土场边坡基底处理方式优选研究[J]. 煤炭科学技术,2021,49(2):77−84. ZHAO Hongze,SHAO Zhiben,DU Hairui,et al. Study on optimization of treatment methods for inclined basement of inner dump slope in anticline area of open-pit mine[J]. Coal Science and Technology,2021,49(2):77−84.
[19] 张 萌,刘宏磊,乔文光,等. 排土场稳定性影响因素显著性研究[J]. 中国安全科学学报,2021,31(2):69−75. doi: 10.16265/j.cnki.issn1003-3033.2021.02.010 ZHANG Meng,LIU Honglei,QIAO Wenguang,et al. Significant research on influencing factors of dump stability[J]. China Safety Science Journal,2021,31(2):69−75. doi: 10.16265/j.cnki.issn1003-3033.2021.02.010
[20] 曹兰柱,黄浩轩,王 东,等. 胜利东二露天矿南帮滑坡区压煤回采边坡稳定性研究[J]. 煤炭科学技术,2019,47(10):81−87. doi: 10.13199/j.cnki.cst.2019.10.009 CAO Lanzhu,HUANG Haoxuan,WANG Dong,et al. Research on slope stability of coal recovery on landslide area underlying southern slope in Shengli No. 2 East Open-pit Mine[J]. Coal Science and Technology,2019,47(10):81−87. doi: 10.13199/j.cnki.cst.2019.10.009
-
期刊类型引用(5)
1. 詹召伟，葛俊岭，徐振，刘涛. 基于高分子增透剂的煤层注水增渗作用分析. 能源与环保. 2024(05): 84-90 . 
百度学术
2. 赵琛，赵炯博. 高压环境下油页岩物理性质及演化规律研究. 能源与环保. 2024(11): 269-274 . 百度学术
3. 张村，贾胜，王永乐，赵毅鑫，陈彦宏，王方田. 煤样CT扫描重构研究进展：原理、方法及应用. 煤炭学报. 2024(S2): 800-820 . 百度学术
4. 李刚，张亚兵，杨庆贺，邹军鹏，才天，刘航，赵艺鸣. 基于Real-ESRGAN的岩石CT图像超分辨率重建. 工矿自动化. 2023(11): 84-91 . 百度学术
5. 李小超，李东涛，尚学锋，刘涛. 基于高分子改性的煤层注水增透剂研发及性能研究. 中国矿业. 2023(12): 177-183 . 百度学术
其他类型引用(3)
计量
- 文章访问数: 126
- HTML全文浏览量: 13
- PDF下载量: 44
- 被引次数: 8
