2. 西南石油大学 土木与建筑学院, 四川 成都 610500;
3. 中国科学院 岩土力学与工程国家重点实验室, 湖北 武汉 430071;
4. 东华理工大学 地球科学学院, 江西 南昌 330013
2. School of Civil Engineering and Architecture, Southwest Petroleum University, Chengdu 610500, China;
3. State Key Laboratory of Geomechanics and Geotechnical Engineering, Chinese Academy of Sciences, Wuhan 430071, China;
4. School of Earth Sciences, East China University of Technology, Nanchang 330013, China
岩体质量评价是随着各种工程的开展而发展起来的, 在工程建设的各个阶段中具有重要的作用.在高放废物地质处置工程中, 处置库围岩作为隔绝放射性核素和生物圈的最后一道屏障, 是影响地下处置库长期稳定性和安全性的关键因素[1].因此围岩岩体质量研究显得尤为重要.
近几十年, 随着岩土工程的发展, 国内外很多专家对不同工程领域的岩体质量评价进行了深入研究, 并提出相应的评价方法.例如传统的RQD法、RMR法、Q系统分类法、国标《工程岩体分级标准》水利水电工程围岩工程地质分类HC法等围岩分级和岩体质量评价方法.随着科学技术的不断进步, 一些系统理论被引入到围岩分级和岩体质量评价中, 如可拓学理论、分形理论、未确知测度理论、粗糙集和人工神经网络.
针对高放废物地质处置工程, 文献[2]根据高放废物地质处置工程不同阶段的特点建立与之相适应的岩体工程质量评价体系; McEwen等[3]提出了RSC指标进行岩体质量评价; 宗自华等[1]开展了高放废物地质处置岩体评价技术及其应用研究, 初步获得了甘肃北山预选区岩体结构面分维特征参数, 并采用Q法、RMR等方法对预选区岩体质量进行了分析; 陈亮等[4]以Q系统为基础提出了中国高放废物地质处置工程岩体质量评价方法.总体而言, 与高放废物地质处置工程不同阶段的特点相对应的岩体质量评价技术体系尚未形成.
目前我国高放废物地质处置工程还处在场址比选阶段, 此阶段主要进行预选区地表岩体节理调查工作.基于地表露头节理调查统计结果, 本文提出了利用JSR指标进行岩体质量评价, 并应用到芨芨槽、算井子和新场地段岩体质量分级中, 为高放废物地质处置场址比选提供依据, 同时讨论了JSR指标与传统Q指标的关系,并修正了JSR指标的分级范围.
1 节理调查方法和内容在野外节理调查中, 首先通过实际踏勘, 在研究范围内选择发育良好的露头, 用GPS确定露头的位置, 然后采用测线法进行节理调查.图 1给出了节理与测线的关系, 将节理分为Ⅰ(相交型)、Ⅱ(延长线相交型)和Ⅲ(不相交型)三种类型, 并测量露头内每一条节理几何参数.为了准确快捷测量, 布置测线时, 应尽量使测线与所有节理迹线相交.主要调查节理产状、间距(密度)和迹长等几何参数[5].
岩体质量评价是指利用岩体质量指标评定岩体的工程地质性质.由于各类工程考虑的侧重点不同, 故岩体质量评价指标和参评因素也不相同.目前国内外各类岩体质量评价分级所考虑的主要因素大致可概括为岩体完整性、结构面条件和环境因素[6].
随着高放废物处置工程的不断发展, 其相应的工程岩体质量评价方法和围岩分级研究也越来越受重视.目前我国高放废物处置研究正处在场址筛选阶段, 因此有必要建立适合于我国的处置库场址评价方法.
围岩作为处置库的天然屏障, 是高放废物与人类生存环境隔离的关键.花岗岩作为高放废物处置的首选围岩, 通过实验得出:花岗岩具有孔隙度小、含水率极低和吸水性差且单轴抗压强度Rc多在170 MPa左右, 大于60 MPa属于坚硬岩石[7-9].但花岗岩中含有一些节理裂隙, 这些裂隙是核素伴随地下水运移的主要通道, 这样就容易对生物圈造成破坏.而高放废物处置场址比选工程主要考虑两方面的问题:核素迁移和围岩稳定.由于花岗岩属于硬岩, 故围岩稳定性基本可以保证; 而花岗岩岩体中节理分布特征成为主要研究方向.因此岩体完整性成为高放废物处置地表岩体质量评价中的主要影响因素.
鉴于岩体完整性对高放废物处置地表岩体质量起决定性作用, 根据国际岩石力学学会的建议[5], 岩体完整性可以用岩石块体大小来定量描述.目前岩石块体大小用节理间距或Jv(节理体密度)来计算.
但是在实际工程中, 尤其是当节理组数较多时, 不可能看到一个露头全部的岩体块体构成情况.而地表岩体的块体大小主要由节理间距、组数和延伸情况等决定.图 2给出了节理间距、组数和延伸情况对岩石块体大小的影响.鉴于上述原因, 针对高放废物处置场址比选阶段, 利用JSR指标进行节理发育程度和地表岩体质量评价:
(1) |
式中:Da为全部露头的平均迹线中点面密度; L为全部露头的平均迹长; Wn为节理优势组数量.由式(1) 可知, JSR的3个参数可以通过地表节理调查统计获得, 根据它们各自的评分值就可以算出JSR值的大小, 从而进行岩体质量分级[10].表 1给出了JSR值对应的节理发育程度和地表岩体质量等级.
高放废物处置甘肃北山预选区新场、芨芨槽和算井子岩体位于甘肃省西北部, 属于塔里木板块北缘北山地区南带.研究区内岩体主要为花岗岩类, 岩性主要是花岗闪长岩和二长花岗岩且呈大面积岩基广泛发育.花岗岩的岩体露头出露新鲜, 几乎没有蚀变, 地貌特征均为强烈夷平和平缓起伏的石漠或低丘[11].
图 3给出了地表岩体节理发育情况.可以看出:在甘肃北山预选区地表岩体发育大量的Ⅳ级结构面, 即节理.地表岩体大量发育的构造节理形成雨水和地下水连通的通道, 对高放废物处置工程长期稳定性具有一定的影响.因此在地表调查构造结构面的基础上, 如何对预选区岩体进行比选是工程的关键问题.本文将通过对新场、芨芨槽和算井子岩体质量评价来确定处置库建设场址.
通过数学与概率统计法对新场、芨芨槽和算井子地段地表野外地质调查获得的节理几何特征进行统计处理, 对3个参数进行评分, 最后计算每个露头点的JSR值并绘制等值线图.以算井子地段为例, 露头点N2节理组数为2组, 相应的评分值为3;节理中点面密度统计值0.02 m-2, 相应的评分值为4;节理迹长为23.09 m, 相应的评分值为15.故露头点N2的JSR值为180.同时得出表 2算井子地段地表露头点对应的JSR值.根据表 2的结果绘制JSR指标等值线图, 如图 4所示.
根据地表岩体JSR指标等值线图统计各级JSR值所占面积分数并计算整个地段的JSR指标的期望值,由表 3可知:芨芨槽地段地表岩体JSR期望值小, 说明芨芨槽地段地表岩体节理发育程度轻微、岩体坚硬、完整性较好且基本稳定, 岩体质量等级为Ⅱ级; 而算井子地段地表岩体节理发育程度中等, 岩体坚硬但比较破碎,稳定性较差, 岩体质量等级为Ⅲ级.因此可以认为芨芨槽地段地表岩体质量较好, 更适合作为高放废物处置库建造预选场址.
在高放废物地质处置场址比选阶段, 针对同一场址, 利用不同的方法评价岩体质量的结果可能存在差异, 因此以算井子地表岩体为研究对象, 建立JSR指标与传统Q指标之间的定量对应关系并修正其分级范围.图 5给出了JSR与Q的定量关系, 在相关系数为0. 65的情况下, 关系式为
(2) |
利用式(2) 和Q指标分级结果对JSR指标进行修正, 具体结果为:JSR值区间[0~25], [25~141], [141~373], [373~489]和[489~720]相对应的岩体质量等级分别为Ⅰ, Ⅱ, Ⅲ, Ⅳ和Ⅴ.
5 结论1) 采用测线法调查地表露头节理几何特征, 利用概率论和数理统计方法进行节理几何特征统计, 然后提出了基于JSR指标进行地表岩体质量评价方法.
2) 利用JSR指标对甘肃北山地表岩体进行节理发育程度和地表岩体质量评价, 结果表明, 芨芨槽地段地表岩体质量较好, 更适合作为高放废物处置库建造预选场址.
3) 研究了JSR指标与传统Q指标的关系并修正了JSR指标的分级范围, 即修正后的分级区间为[0~25], [25~141], [141~373], [373~489]和[489~720].
4) 在我国高放废物地质处置北山预选区地表岩体质量评价中, 因JSR指标既考虑了影响工程的主要因素又容易获取所需参数而被使用.但是JSR指标在其他工程如隧道开挖或地下洞室中并不适用.因此有必要对JSR指标进行进一步的研究以应用到其他工程岩体质量评价中.
[1] |
宗自华, 陈亮, 陈强, 等.
高放废物处置北山预选区地表岩体结构与岩体质量相关性及其应用研究[J]. 岩土力学, 2013(sup1): 279–284, 292.
( Zong Zi-hua, Chen Liang, Chen Qiang, et al. Study of relationship between rock structure and rock quality and its application to site evaluation of Beishan HLW repository[J]. Rock and Soil Mechanics, 2013(sup1): 279–284, 292. ) |
[2] | Hagros A.Host rock classification (HRC) system for nuclear waste disposal in crystalline bedrock[D].Helsinki:University of Helsinki, 2006. |
[3] | McEwen T, Kapyaho A, Hella P, et al.Rock suitability classification RSC 2012[R].Helsinki:Posiva Oy, 2012. |
[4] | Chen L, Wang J, Zong Z H, et al. A new rock mass classification system QHLW for high-level radioactive waste disposal[J]. Engineering Geology, 2015, 190: 33–51. DOI:10.1016/j.enggeo.2015.02.006 |
[5] | Barton N. Suggested methods for the quantitative description of discontinuities in rock masses:ISRM[J]. International Journal of Rock Mechanics and Mining Sciences and Geomechanics Abstracts, 1978, 15(6): 319–368. DOI:10.1016/0148-9062(78)91472-9 |
[6] | Tzamos S, Sofianos A I. A correlation of four rock mass classification systems through their fabric indices[J]. International Journal of Rock Mechanics and Mining Sciences, 2007, 44(4): 477–495. DOI:10.1016/j.ijrmms.2006.08.003 |
[7] |
魏翔, 杨春和.
含陡倾角软弱结构面的花岗岩力学特性及声发射特征[J]. 工程科学学报, 2016, 38(12): 1667–1673.
( Wei Xiang, Yang Chun-he. Experimental study of mechanics features and acoustic emission characteristics for granite of weak structure plane with steep dip angle[J]. Chinese Journal of Engineering, 2016, 38(12): 1667–1673. ) |
[8] | Chen L, Wang C P, Liu J F, et al. Effects of temperature and stress on the time-dependent behavior of Beishan granite[J]. International Journal of Rock Mechanics and Mining Sciences, 2017, 93(3): 316–323. |
[9] | Zhao X G, Cai M, Wang J, et al. Strength comparison between cylindrical and prism specimens of Beishan granite under uniaxial compression[J]. International Journal of Rock Mechanics and Mining Sciences, 2015, 76(1): 10–17. |
[10] |
雷光伟. 岩体节理几何特征研究及三维网络模拟[D]. 重庆: 重庆大学, 2015.
( Lei Guang-wei.Study on rock mass joints geometric feature and three-dimension network simulation[D].Chongqing:Chongqing University, 2015. http://cdmd.cnki.com.cn/Article/CDMD-10611-1015971004.htm ) |
[11] | Zhao X G, Wang J, Qin X H, et al. In-situ stress measurements and regional stress field assessment in the Xinjiang candidate area for China's HLW disposal[J]. Engineering Geology, 2015, 197: 42–56. DOI:10.1016/j.enggeo.2015.08.015 |