Processing math: 100%
  • 全国中文核心期刊
  • 中国科技核心期刊
  • 美国工程索引(EI)收录期刊
  • Scopus数据库收录期刊

基于二阶功准则的NPR锚索支护效果研究

胡杰, 何满潮, 李兆华, 冯吉利, 肖勇杰

胡杰, 何满潮, 李兆华, 冯吉利, 肖勇杰. 基于二阶功准则的NPR锚索支护效果研究[J]. 岩土工程学报, 2021, 43(10): 1870-1877. DOI: 10.11779/CJGE202110013
引用本文: 胡杰, 何满潮, 李兆华, 冯吉利, 肖勇杰. 基于二阶功准则的NPR锚索支护效果研究[J]. 岩土工程学报, 2021, 43(10): 1870-1877. DOI: 10.11779/CJGE202110013
HU Jie, HE Man-chao, LI Zhao-hua, FENG Ji-li, XIAO Yong-jie. Anchorage effect of NPR cable based on second-order work criterion[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(10): 1870-1877. DOI: 10.11779/CJGE202110013
Citation: HU Jie, HE Man-chao, LI Zhao-hua, FENG Ji-li, XIAO Yong-jie. Anchorage effect of NPR cable based on second-order work criterion[J]. Chinese Journal of Geotechnical Engineering, 2021, 43(10): 1870-1877. DOI: 10.11779/CJGE202110013

基于二阶功准则的NPR锚索支护效果研究  English Version

基金项目: 

国家重大研发计划重点专项项目 2016YFC0600901

国家自然科学基金项目 41172116

国家自然科学基金项目 U1261212

中国矿业大学(北京)深部岩土力学与地下工程国家重点实验室开放基金项目 SKLGD UEK1718

详细信息
    作者简介:

    胡杰(1992— ),男,博士研究生,主要从事岩土工程数值模拟方面研究工作。E-mail:hi_hujie@163.com

    通讯作者:

    李兆华, E-mail:472697969@qq.com

  • 中图分类号: TU457

Anchorage effect of NPR cable based on second-order work criterion

  • 摘要: NPR锚索因其具有恒定工作阻力和较大变形量的特性,已广泛应用于巷道的加固和监测预警中。为了研究NPR锚索在露天矿山边坡中的支护效果,首先基于有限差分软件FLAC3D中的锚索单元开发出NPR锚索单元,并通过数值拉伸试验验证了NPR锚索比普通锚索具有更优异的吸能特性;之后引入二阶功作为边坡稳定性评价准则,并将其植入显式动力学算法中,详细讨论了其合理性和可靠性;最后,以浩尧尔忽洞金矿西采场北帮开挖诱发的滑坡为研究对象,对现有普通矿用锚索和NPR锚索支护体系下的边坡稳定性进行对比分析。根据数值模拟结果可知,NPR锚索支护体系下的边坡稳定性明显优于普通矿用锚索。
    Abstract: Due to the characteristics of constant working resistance and steady large deformation, the NPR cable has been widely applied in monitoring and reinforcement of a large number of roadways. In order to study the anchorage effect of NPR cables in open-pit slope, firstly, the NPR cable element is developed based on the cable element in the finite difference software FLAC3D, and it is proved by the tensile tests that the NPR cable has better energy absorbing capacity than the traditional one. Secondly, the second-order work criterion is then implemented, and the rationality and advantages of the second order work criterion as a safety factor in explicit numerical algorithm are discussed. Finally, the rock slope induced by the excavation at the north side of the west stope in Haoyaoerhudong Opencast Gold Ore is taken as the research object, and the stability of rock slope under the existing traditional mine cable support system and NPR cable support system is compared and analyzed. According to the numerical results, the global second-order work and kinetic energy of each loading step have the opposite variation tendency, and the slope stability of NPR cable support system is obviously better than that of the traditional mine cable.
  • 高压实膨润土作为高放废物深地质处置库首选缓冲/回填材料,具有维持结构稳定、延缓地下水入渗、阻滞核素迁移和传导核衰变热等重要缓冲功能(图 1)。在处置库长达上万年的服役期间,高压实膨润土将长期受到围岩地下水及其化学成分的侵入作用,导致矿物相态和微观结构不断演化,缓冲性能逐渐衰减,最终威胁处置库的长期稳定与运营安全。因此,深入剖析膨润土与孔隙溶液的相互作用机制是研究处置库复杂近场环境下膨润土缓冲性能长期演化的关键。

    图  1  高放废物深地质处置库及多重屏障缓冲示意图
    Figure  1.  Schematic diagram of deep geological repository and multiple barrier buffer for high-level waste

    为此,众多学者针对高压实膨润土在不同孔隙溶液条件下缓冲性能的演化规律开展了大量试验研究,并建立孔隙溶液对膨润土膨胀力、渗透系数等水力-力学特性和吸附量、扩散系数等化学阻滞特性影响的理论模型[1-4]。化-水-力膨胀模型方面,Guimarães等[5]以阳离子交换量作为孔隙溶液的化学行为描述,对巴塞罗那膨胀模型(BExM)[6]进行了耦合改进;Lu等[2]、Dominijanni等[7]分别以渗透吸力和迁移特性参数变化作为孔隙溶液的化学行为描述,改进了有效应力的表达。Bennethum等[8]引入化学势,利用自由能和热力学第二定律建立了膨胀多孔介质的本构框架;徐永福[9]基于膨润土吸水膨胀特性,用分形表面吸附理论建立了膨润土膨胀变形的统一理论。渗透模型方面,Ruan等[10]对Kozeny-Carmen模型进行双孔修正,提高了在压实膨润土中的适用性;Komine[11]基于阳离子交换量提出了多组分阳离子在膨润土中等效渗透系数的计算模型。吸附扩散模型方面,Torstenfelt等[12]将菲克第二定律对多孔介质的扩散模型进行了吸附耦合修正;Muurinen等[13]排除层间孔隙影响,以有效孔隙代替总孔隙进行修正,将孔隙溶液对吸附扩散的化学行为简化为对有效孔隙和线性吸附分配系数的影响。

    然而,现有研究表明,膨润土晶层胀缩和扩散双电层厚度受孔隙溶液中阳离子交换的影响[14],且交换难易程度与其干密度和水化约束条件有关[15-16],并表现出对不同价态阳离子的选择性[17-18]。由于膨润土对阳离子的响应明显大于阴离子,因此,以吸力或阳离子交换作为膨润土在溶液中有效应力和水化响应的修正依据,忽略了缓冲性能随溶液扩散的演化过程、各参数互相影响和溶液类型的差异性,也缺乏描述溶液作用下孔隙结构变化的过渡模型和化学吸附模型,尤其是微观吸附模型、细观孔隙迁移模型和宏观力学体变模型在尺度上存在极大差异。至今尚未建立起跨尺度和溶液作用描述的缓冲性能模型。

    事实上,高压实膨润土在溶液中的响应十分复杂,既有晶层胀缩、集合体裂解、扩散双电层水化、阳离子交换等不改变矿物晶胞结构的物理作用[19-20];也有产生矿物晶胞相变、化学胶结的化学作用[20-21]。本文在阐述溶液对高压实膨润土缓冲性能影响的基础上,总结分析了与孔隙溶液的物理作用机制和仍待研究的重点方向,以期为高放废物深地质处置库工程屏障在近场复杂环境下长期缓冲性能演化评估、设计建造和运营管理提供依据。

    膨润土是以层状蒙脱石为主要矿物的高膨胀性黏土。在蒙脱石晶体结构中,硅氧四面体中的四价Si可被三价Al同晶置换[22],铝氧八面体中的三价Al又可被二价Mg,Ca,Fe等同晶置换,使得每个单位晶胞有0.2~0.6的永久负电荷[23]。为了平衡晶层负电荷,晶层间吸附有大量可交换金属阳离子和极性水分子。蒙脱石特殊结构为吸水膨胀和离子吸附提供了强劲动力;其膨胀挤密和吸附黏滞效应也使膨润土具有高愈合性、低渗透性、强阻滞性等重要特性。

    高压实膨润土可分成单元晶层,层叠体(准晶体)和集合(聚)体3个结构层次,如图 2所示。其中,层叠体是由数个至数百个单元晶层层状堆积形成,其数量与蒙脱石层间的阳离子类型、饱和度和孔隙溶液有关;集合体是由若干个层叠体聚集组成,尺度一般在微米级别。同时,这3种基本单元也构成了高压实膨润土的3种主要结构层次的孔隙,即晶层间孔隙、层叠体间孔隙和集合体间孔隙[24]

    图  2  高压实膨润土微观结构组成[24]
    Figure  2.  Microstructural composition of highly compacted bentonite[24]

    通常,层叠体结构内晶层间可吸附1~4层极性水分子。这些水分子一部分被强吸附在晶层表面形成氢原子向内的定向排列,另一部分则以水合阳离子形成层内溶液[25]。这些定向排列的强吸附水通常密度大于1.0 g/cm3,黏滞性比自由水大100倍左右,活动度、扩散性和氢键作用强度远比自由水弱,一定程度上表现出固体性质[26],使得水分子难以穿越晶层。因此,层叠体作为基本微观结构单元,也是与孔隙溶液直接作用的基本受力单元。层叠体内晶层膨胀和层叠体外扩散双电层膨胀共同构成了膨润土的总体膨胀。

    膨润土工程屏障中的孔隙溶液主要由近场围岩地下水侵入形成,其化学成分十分复杂,与处置库围岩和场地环境密切相关。比利时、瑞士和法国为代表的黏土岩型处置库主要考虑黏土岩孔隙水,如法国处置库围岩地下水TDS约为5700 mg/L,主要阳离子中Na+含量高达1035 mg/L,Ca2+含量为280 mg/L[27];日本沉积岩型处置库主要考虑海水的侵入,其地下水TDS高达42400 mg/L,主要阳离子中Na+含量高达11040 mg/L,Ca2+含量为440 mg/L[28];西班牙、韩国和中国的花岗岩型处置库主要为花岗岩裂隙水,以中国北山地下水为例,其TDS值介于1793.3~3594.5 mg/L,pH介于7.5~8.5,化学相为Cl·SO4-Na型,其中Na+含量高达494~1036 mg/L,Ca2+含量为90.2~183.0 mg/L [29]

    同时,膨润土中的孔隙水长期遭受衰变热、气体运移、围岩应力、膨润土吸水膨胀等复杂多相多场耦合作用后,会因干湿循环作用而产生盐分浓缩积聚,形成浓度更高、化学成分更复杂的孔隙溶液[20]。此外,混凝土衰解和膨润土溶蚀转变等化学作用也会对孔隙溶液的组成产生重要影响[30]

    处置库中拼砌的高压实膨润土屏障服役后,将长期处于非饱和状态,并在吸力作用下从围岩中吸收水分,发生膨胀愈合形成机械屏障和化学屏障阻滞核素迁移。然而,缓冲屏障建设时不可避免地在块体间、块体与围岩间留下大量施工接缝,成为地下水及化学成分入渗的优势通道,使部分膨润土很快达到饱和状态。

    处置库服役初期,高压实膨润土主要通过水汽吸湿产生水化膨胀,典型的膨胀力时程曲线见图 3 [31]。当吸力大于149.6 MPa时,膨润土处于初始水化阶段,主要发生蒙脱石晶层膨胀,膨胀力大小取决于晶层间阳离子类型,不受干密度(初始孔隙率)影响。随着水化的持续进行,吸力逐渐降低,膨胀力大小与干密度呈现正相关,增长趋势因干密度不同而表现出显著差异。当吸力低于38 MPa时,膨胀力因层叠体大规模裂解而发生累积“楔”力释放和孔隙坍塌而表现出减小,并因层叠体表面积增多使得扩散双电层总膨胀增大而再次增长,呈现出“先增加—后减小—再增加”的双峰演化特征[32]

    图  3  不同吸力下压实膨润土膨胀力时程曲线[31]
    Figure  3.  Time-swelling pressure curves with increase of suction[31]

    此外,膨润土吸湿过程中形成的晶层水化厚度、层叠体裂解后的晶层数与水化吸力或含水率对应。Saiyouri等[25-26]测试了不同吸力下MX80-Na基膨润土的晶层衍射数据,并基于Pons等[33]统计得到层间吸附水分子层数和层叠体内晶片层数随吸力变化规律(图 4)。结果表明,吸力在7~50 MPa时,层间含有2层水分子;吸力小于0.06 MPa时,层间最多含有4层水分子。相应地,当吸力从50 MPa降至7 MPa时,层叠体内晶层从350层裂解至150层;吸力降至3 MPa时,晶层数量骤减至3层。可见,蒙脱石层叠体在吸力低于50 MPa时开始大规模裂解,在吸力低于7 MPa的高饱和阶段,层叠体已基本裂解分散。

    图  4  吸力对蒙脱石水化层间距和层叠体厚度的影响[25-26]
    Figure  4.  Effects of suction on interlayer distance and laminate thickness of montmorillonite[25-26]

    根据水化膨胀主导机制差异,通常将膨润土由非饱和向完全饱和的水化膨胀过程分为3个阶段(图 5)。①第1阶段,蒙脱石层间吸水膨胀;②第2阶段,层叠体膨胀裂解,并填充层叠体内孔隙;③第3阶段,集合体逐步崩解、大孔隙坍塌并被裂解的薄层叠体重新构建,此时,薄层叠体迅速增加,并在薄层外表面快速形成扩散双电层继续水化膨胀[34]

    图  5  膨润土水化膨胀机制
    Figure  5.  Hydration swelling mechanism of bentonite

    溶液中高压实膨润土水化膨胀同样经历膨胀力快速发展的晶层膨胀、膨胀力不断调整的层叠体裂解和膨胀力再增加的扩散双电层膨胀3个阶段。溶液入渗初期,由于水分运移速率大于离子,当离子进入膨润土时,膨润土已水化膨胀。因此,当离子依次通过大小孔隙时,势必打破不同尺度孔隙溶液的吸力平衡,引起土颗粒表面阳离子竞争吸附和孔隙结构不断调整,尤其以溶液浓度对膨胀力演化影响最为显著。如图 6所示[30],NaCl浓度越高,晶层膨胀力越小,孔隙塌陷时膨胀力下降越大,扩散双电层膨胀越小。更多研究表明[35-39],溶液浓度增加对最终膨胀力既有促进作用,又有抑制作用,并与离子类型有关。此外,孔隙溶液引起的膨润土化学体变也因阳离子类型不同而呈现差异性(图 7)。Na基膨润土在NaCl溶液的盐—淡循环中,化学体变具有可逆性;而在KCl或CaCl2溶液的盐—淡循环中却呈现出不可逆收缩变形[40-41]。可见,溶液离子对高压实膨润土膨胀的影响机制是多因素共同作用的结果,物理作用的可逆性也会因干密度和水化约束条件不同而不同。

    图  6  溶液作用下膨胀力三阶段演化特征[30]
    Figure  6.  Three stages of swelling pressure hydrated with salt solutions[30]
    图  7  溶液对膨润土最终膨胀力的影响[35-39]
    Figure  7.  Effects of ion species and concentration on final swelling pressure of compacted bentonite[35-39]

    孔隙溶液对晶层膨胀的作用与渗透吸力有关。蒙脱石晶层表面定向排列的强吸附水分子一定程度上表现出固体性质[26],使得水分子难以通过晶层。当层叠体外溶液的溶质吸力小于层内溶质吸力时,层间吸附阳离子和水分子不能自由扩散至层外,仅能通过层叠体外水分子向内迁移达到吸力平衡。当晶层受到约束,限制水分子向内迁移时,层叠体将产生向内的渗透力Oπ,大小与吸力差成正比。当晶层在约束力Pc作用下达到平衡时,膨胀力Ps等于约束力和渗透力之和,即Ps = Pc + Oπ,如图 8所示[42]

    图  8  层叠体颗粒受力平衡示意图
    Figure  8.  Schematic diagram of stress balance of layered particles

    图 8可知,膨润土最大膨胀力和膨胀变形并不是发生在浓度为零的去离子水环境,而是蒙脱石晶层内外溶质吸力差降为零的临界吸力溶液中[36, 38]。当孔隙溶质吸力低于临界吸力时,孔隙溶液浓度的增加会促进晶层膨胀。进一步分析图 79可知[15, 43],临界吸力的大小与约束力(如恒体积、恒定荷载或自由膨胀等水化边界条件)有关,约束力越大,层间膨胀量越小,此时层间溶质吸力越大,对应的临界吸力也越大。恒体积条件下临界吸力对应的浓度小于0.05 mol/L,而自由膨胀条件下普遍小于0.02 mol/L。

    图  9  离子类型和浓度对膨润土膨胀指数的影响[15, 43]
    Figure  9.  Effects of ion species and concentration on free swell index of bentonite[15, 44]

    当层叠体外溶液的溶质吸力大于临界吸力时,层叠体内水分子可向层外运移、层叠体外离子也向层内扩散以达到吸力平衡。因此,孔隙溶液对晶层膨胀的抑制包括阳离子交换和层间排水,分别对应渗透固结和渗透诱导固结两种体变机制。

    阳离子内迁交换对晶层的影响与层间溶质吸力密切相关。发生阳离子交换时,晶层更倾向选择高价态或强吸附性的低水合吸力离子,以降低整体水化势能。如图 10(a)[45],层间相同正价浓度中,阳离子水化吸力大小为Na+ > K+ > Ca2+,相应地,膨润土膨胀性能相对大小为Na基膨润土>K基膨润土>Ca基膨润土。此外,层外阳离子内迁还与层叠体约束力密切相关。约束力增加在提高临界浓度的同时,也抑制了晶层膨胀间距和扩散孔径,增加阳离子交换难度,从而使层叠体内外吸力差更倾向于通过层间排水达到平衡。在恒体积膨胀时,相同浓度条件下,阳离子溶质吸力越大的溶液对膨胀力的抑制作用越强,高价态阳离子的抑制作用程度低于低价态阳离子,且水合离子半径越小,影响越大(图 11[16],各阳离子影响程度依次为:Na+ > K+ > Ca2+ > Mg2+ > Al3+ > Fe3+。恒定荷载条件下(图 12),土体膨胀使孔隙增大,当K+进入层叠体,因其离子半径与晶层中硅氧四面体六联环空间相近,形成K-O配位多面体而具有更强的交换优势,表现出比Na+更强的抑制作用[15]。自由膨胀条件下(图 9),层叠体不受约束,晶层完全膨胀,孔隙进一步增大,高价态Ca2+也能进入层叠体内,从而表现出比Na+更强的抑制作用[15, 43]

    图  10  溶液浓度与渗透吸力的关系[45]
    Figure  10.  Relationship between concentration of different solutions and osmotic suction[45]
    图  11  离子类型对高压实膨润土膨胀力的影响[16]
    Figure  11.  Effects of ion species on final swelling pressure of compacted bentonite[16]
    图  12  离子类型对高压实膨润土膨胀变形的影响[15]
    Figure  12.  Effects of ion species on final swelling deformation of compacted bentonite[15]

    膨润土在孔隙溶液作用下,层叠体表面将产生扩散双电层膨胀。扩散双电层厚度可为[14]

    1K = ε0εRkBT2n0e2v2 (1)

    式中:1/K为双电层厚度;ε0为真空介电常数,取8.854×10−12 C2·J−1·m−1εR为孔隙溶液相对介电常数;kB为玻尔兹曼常数(1.38×10−23 J·K−1);e为单位电荷带电量(1.602×10−19 C);n0为电解质浓度;v为离子价态。

    由此可知,扩散双电层厚度主要与孔隙溶液浓度和离子价有关,浓度越大,离子价越高,扩散双电层厚度越小。如图 12所示,相同浓度Na+溶液中,SO42-比Cl-对膨胀抑制作用更强[15]。而且低价阳离子的抑制作用强于高价阳离子[35, 37],但这一现象与粉末状膨润土的自由膨胀刚好相反[15, 43]。可见,高压实膨润土中晶层膨胀量大于扩散双电层膨胀量,且高价阳离子很难进入晶层微孔,多处于层叠体大孔中。因此,孔隙溶液对扩散双电层膨胀影响较晶层膨胀更为显著。进一步研究表明,膨润土完全水化后,集合体裂解产生的层叠体厚度也不均一。土体内存在层间微孔(约2 nm)、薄层叠体间小孔(10~150 nm,峰值孔径20 nm)和厚层叠体间大孔(大于150 nm)[24]。孔隙溶液浓度相同时,离子类型对膨润土膨胀力和饱和渗透系数大小的影响规律依次为Na+ > K+ > Ca2+ > Mg2+ > Al3+ > Fe3+[16]。这一现象表明,离子很难在微孔中运移,主要通道为层叠体间孔隙,且受水化双电层阻滞效应的影响,实际有效孔隙更小(图 13)。实际上,扩散双电层中的水分子和离子越靠近层叠体表面,吸附固定作用越强[45],提供给溶液运移的有效孔隙直径更小。由此可见,孔隙溶液对晶层和扩散双电层膨胀的抑制作用,显著增加了有效孔隙通道并减小弯曲度,使膨润土渗透性能(图 14)和核素离子扩散性能(图 15)明显提高[16, 46]。高压实膨润土力学缓冲性能和化学缓冲性能也因其它离子的加速侵入而衰减。

    图  13  Stern理论扩散双电层构造[45]
    Figure  13.  Diffusion double-layer structure of Stern theory[45]
    图  14  离子类型及浓度对膨润土饱和渗透系数的影响[46]
    Figure  14.  Effects of ion species and concentrations on saturation permeability coefficient of compacted bentonite[46]
    图  15  离子强度对膨润土中核素离子扩散的影响[47]
    Figure  15.  Effects of ionic strength on diffusion of nuclides[47]

    膨润土对溶液离子吸附作用包括表面络合(化学)吸附和静电(物理)吸附,分别对应两类表面活性位点。如图 16所示,表面络合吸附主要是离子与蒙脱石矿物表面硅/铝原子形成羟基络合位点(Si/Al―OH)形成共价键结合的内层配合物,化学反应式为

    nSi/AlOH+Mn+(Si/AlO)nM+nH+ (2)
    图  16  蒙脱石表面吸附作用示意图[49]
    Figure  16.  Schematic diagram of surface adsorption of montmorillonite[49]

    式中:≡Si/AlOH 为羟基络合位点,Mn+为吸附阳离子。

    静电吸附主要是水合离子与矿物表面在范德华力作用下形成外层配合物,也叫离子交换吸附。发生在层间的静电吸附(离子交换反应)可表示为

    nXNa+Mn+(X)nM+nNa+ (3)

    式中:X为极性负电荷位点。

    表面络合吸附与吸附离子形成的共价键比静电吸附形成的范德华键的亲和力高出一个或几个数量级[48-49]。因此,静电吸附的离子在环境条件改变时容易发生解吸附和置换,而络合吸附的离子不易解吸附。

    溶液中H+或OH-的浓度不仅影响膨润土表面吸附位点的活性,也会改变离子水解结合形式,其在表面吸附位点的反应机制见表 1[50]。通过平衡常数可知,酸性环境中,H+率先与蒙脱石矿物边缘的铝羟基(Al―OH)发生质子化反应而带正电;酸性更强时,H+又可作为竞争阳离子吸附占据阳离子交换位点,影响膨润土对目标离子的吸附活性。膨润土易与异性电荷离子产生静电吸附,而与同性电荷离子发生络合吸附。当膨润土表面由负电荷变为正电荷时,静电吸附发生吸附反转,表现出对阴离子的吸附活性,并抑制对阳离子的吸附作用。碱性环境中,OH-对蒙脱石表面羟基络合位点(Si/Al―OH)的去质子化反应(表 1),会使原本羟基络合位点带负电,进入准络合状态,增强络合活性。如图 17所示[51-55],核素离子吸附率在强酸环境时,均因强竞争吸附和质子化而随pH减小而迅速降低。

    表  1  Na基膨润土表面活性位点的反应机制[51]
    Table  1.  Reactions of active sites on Na-bentonite[51]
    表面活性位点反应 双层模型参数(lgK) 非静电模型参数(lgK)
    XNaX+Na+ -1.58 -1.58
    XNa+H+fXH+Na+ 1.50 2.10
    AlOH+H+fAlOH + 2 6.15 5.83
    AlOH+OHfAlO + H2O -9.27 -7.02
    SiOH+OHfSiO + H2O -9.06 -8.75
    下载: 导出CSV 
    | 显示表格
    图  17  pH对GMZ-Na基膨润土吸附核素的影响[52-56]
    Figure  17.  Effects of pH on adsorption of nuclides by GMZ-Na bentonite[52-56]

    除对膨润土吸附位点的活性影响外,pH对核素水合离子形态的影响也极为复杂。核素离子由于价态高,通常存在多级水解反应,离子形态随pH升高会结合多个OH-而减小正电荷,从而提高被吸附量(图 17),也可能因带负电而减小吸附量,更有如U(Ⅵ)存在多种同价异构离子而表现出特殊吸附特性(图 18[53]。此外,碱性环境下核素离子的沉淀作用也会导致核素离子去除率迅速升高的重要原因之一。膨润土对核素离子静电吸附和络合吸附的主次关系与孔隙溶液pH密切相关。以阳离子形态存在的核素离子在弱酸溶液中主要以静电吸附为主,在强酸和碱性溶液中主要以络合吸附为主[50, 56]

    图  18  不同温度中U(Ⅵ)离子形态分布随pH的变化[49]
    Figure  18.  Species distribution of U(Ⅵ) ions with pH change under different temperatures[49]

    GMZ-Na基膨润土表面活性位点的分布密度见表 2,其中蒙脱石矿物表面阳离子交换位点的分布密度约为羟基络合位点的5倍(络合位点仅占17.27%)[50]。两类吸附位点对离子吸附响应有着显著的活性差别,Na+,Ca2+,Mg2+,K+等金属离子水解能力较强,与羟基的络合能力弱,吸附主要发生在阳离子交换位点,对羟基络合位点影响较小。存在多种离子时,有限吸附位点对目标离子的吸附量会因竞争吸附而降低[17-18, 57],并表现出选择吸附性。

    表  2  GMZ-Na基膨润土表面活性位点的分布密度[50]
    Table  2.  Distribution density of active sites on surface of GMZ-Na bentonite [50]
    活性位点类型 密度
    蒙脱石层间阳离子交换位点,(X 1.35×10-5 mol/m2, (66.5±7 cmol/kg)
    铝氧八面体边缘羟基位点,(AlOH 9.39×10-7 mol/m2, (4.61 cmol/kg)
    硅氧四面体边缘羟基位点,(SiOH 1.88×10-6 mol/m2, (9.23 cmol/kg)
    下载: 导出CSV 
    | 显示表格

    研究表明,离子在黏土矿物表面的竞争吸附优势主要与其带电性、一级水解常数、离子半径、离子强度和体系pH等因素有关。对相同正价电荷的金属离子,一级水解常数越小,水合离子半径越小,与膨润土竞争吸附亲和力越强[58]。而核素离子与背景离子的竞争吸附机制会因该核素离子水解价态的变化而出现不同的促进或抑制结果。如图 19所示,酸性环境下,U(Ⅵ)以阳离子形式存在,离子强度增加会导致竞争优势减弱而吸附减少[59];碱性环境下,U(Ⅵ)离子正价减小反而促进吸附,当水解为阴离子时,则抑制吸附[54]。因U(Ⅵ)水解特殊性,阴离子和阳离子都会对吸附产生较大影响(图 1820),其中K+影响大于Na+和Li+,HCO3影响则因pH变化而变化。离子类型方面,Ca2+会显著降低膨润土对U(Ⅵ)吸附容量,而Na+,K+,Mg2+对吸附影响较小;阴离子中,CO32−对吸附影响最大,HCO3次之,而SO42−,NO3对吸附影响较小(图 21[60]

    图  19  离子强度对GMZ-Na基膨润土吸附U(Ⅵ)的影响[54]
    Figure  19.  Effects of ionic strength on adsorption of U(Ⅵ) ions by GMZ-Na bentonite[54]
    图  20  阴阳离子类型对GMZ-Na基膨润土吸附U(Ⅵ)的影响[54]
    Figure  20.  Effects of cationic and anion species on adsorption of U(Ⅵ) by GMZ-Na bentonite[54]
    图  21  背景离子类型对U(Ⅵ)吸附分配系数的影响[60]
    Figure  21.  Effects of ion species on equilibrium constants of U(Ⅵ) [60]

    孔隙溶液对高压实膨润土的物理作用机制包括晶层膨胀、扩散双电层膨胀和吸附作用,相应的内在影响机理如下所示:

    (1)层叠体是膨润土与孔隙溶液直接发生吸力平衡和吸附平衡的基本单元,且与水化约束力有关。孔隙溶液在低浓度时能够促进晶层膨胀,在高浓度时抑制晶层膨胀,取决于与临界吸力的差;孔隙溶液对层叠体外扩散双电层膨胀具有抑制作用。

    (2)层叠体间孔隙是离子运移的主要通道,孔隙溶液对层叠体膨胀的力-化作用伴随着层叠体裂解和孔隙坍塌调整。孔隙溶液对扩散双电层的抑制作用是导致有效孔隙通道扩大、渗透性和扩散性增大的主要原因。

    (3)膨润土对溶液离子吸附作用包括表面络合作用的化学吸附和静电作用的物理吸附。孔隙溶液对吸附作用的影响与其pH、离子成分及浓度有关。pH对膨润土表面活性位点和核素水合离子带电性质的影响是引起吸附特性变化的主要原因;背景离子的竞争吸附作用致使膨润土对核素离子的吸附量显著减少。

    由于高压实膨润土缓冲性能受蒙脱石矿物层间阳离子、孔隙溶液成分、干密度和孔隙结构的影响十分复杂。现有溶液对膨润土缓冲性能的量化表征研究仍十分有限,国内外学者仅通过考虑化学作用效应对现有模型进行修正,无法从膨润土多尺度结构与孔隙溶液相互作用机理出发构建微宏观跨尺度的缓冲性能预测模型。因此,笔者认为以下工作仍需进一步研究:

    (1)开展蒙脱石层叠体与孔隙溶液相互作用的系统性试验,研究不同孔隙溶液、干密度和水化约束条件下,蒙脱石层叠体膨胀裂解、吸力平衡和吸附平衡的作用特性,量化孔隙溶液对膨润土水化膨胀的影响程度,构建与宏观本构模型的参数关联。

    (2)开展压实膨润土中不同尺度孔隙的等效量化研究,探究不同尺度孔隙对压实膨润土膨胀、渗透、扩散等模型的溶液响应特性,确定孔隙的参数选用。

    (3)开展约束条件下膨润土在动态、多组分溶质中的竞争吸附研究,构建约束条件下高压实膨润土的竞争吸附模型。

    (4)在此基础上,提出考虑化学效应的高压实膨润土缓冲性能演化长期评价方法,为中国深地质处置缓冲/回填材料的选取及设计提供理论和技术支撑。

  • 图  1   NPR锚索结构示意图

    Figure  1.   Components of NPR cables

    图  2   NPR锚索单元实现流程

    Figure  2.   Implementation of NPR Cable element

    图  3   锚索拉伸试验模型

    Figure  3.   Tensile test model of cables

    图  4   锚索拉伸荷载-变形量曲线

    Figure  4.   Tensile force-axial deformation curves of traditional and NPR cables

    图  5   边坡研究区域

    Figure  5.   Research area of rock slope

    图  6   锚杆(索)布置图

    Figure  6.   Distribution of bolts and cables

    图  7   岩质边坡三维地质模型及锚杆(索)布置图

    Figure  7.   Geometric model for rock slope in 3D and distribution of bolts and cables

    图  8   最大剪应变云图

    Figure  8.   Contours of maximum shear strain

    图  9   4种工况下整体二阶功和动能对比图

    Figure  9.   Comparison between global second-order work and kinetic energy of rock slope under 4 working conditions

    图  10   传统矿用锚索及NPR锚索加固下边坡的整体二阶功对比图

    Figure  10.   Comparison of global second-order work of rock slope reinforced by traditional and NPR cables

    表  1   锚杆(索)物理力学参数

    Table  1   Physical and mechanical parameters for bolts and cables

    锚固类别强度/kN杨氏模量/GPa长度/m预应力/kN横截面/m2恒阻力/kN变形量/mm
    NPR锚索1500205277500.22388501000
    矿用锚索1500205277500.2238
    矿用锚杆800205800.373×10-3
    下载: 导出CSV

    表  2   边坡岩体物理力学参数

    Table  2   Physical and mechanical parameters of involved rockmass

    类型ρ /(kg·cm-3)E/GPaν c/MPaφ /(°)TS/MPa
    Ⅲ级片岩2.74.540.250.40350.12
    Ⅳ级片岩2.54.110.300.15280.10
    Ⅳa级片岩2.54.110.300.15280.10
    片岩破碎带2.42.500.350.02180.10
    开挖区2.54.110.300.15280.10
    下载: 导出CSV
  • [1] 何满潮. 露天矿高边坡工程[M]. 北京: 煤炭工业出版社, 1991.

    HE Man-chao. High Slope Engineering in Open-Pit Mine[M]. Beijing: China Coal Industry Publishing House, 1991. (in Chinese)

    [2]

    HE M C, GONG W L, WANG J, et al. Development of a novel energy-absorbing bolt with extraordinarily large elongation and constant resistance[J]. International Journal of Rock Mechanics and Mining Sciences, 2014, 67: 29-42. doi: 10.1016/j.ijrmms.2014.01.007

    [3] 吕谦, 陶志刚, 李兆华, 等. 恒阻大变形锚索弹塑性力学分析[J]. 岩石力学与工程学报, 2018, 37(4): 792-800. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201804002.htm

    LÜ Qian, TAO Zhi-gang, LI Zhao-hua, et al. Elasto-plastic analysis of large deformation cables[J]. Chinese Journal of Rock Mechanics and Engineering,2018, 37(4): 792-800. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201804002.htm

    [4]

    CAI Y, ESAKI T, JIANG Y. A rock bolt and rock mass interaction model[J]. International Journal of Rock Mechanics & Mining Sciences, 2004, 41(7): 1055-1067.

    [5]

    LI C C. A new energy-absorbing bolt for rock support in high stress rock masses[J]. International Journal of Rock Mechanics and Mining Sciences, 2010, 47(3): 396-404. doi: 10.1016/j.ijrmms.2010.01.005

    [6]

    ORTLEPP W D, READ J J. Yieldable rock bolts for shock loading and grouted bolts for faster rock stabilization[J]. Mining Engineering, 1970, 60(3): 12-17.

    [7]

    ST-PIERRE L, HASSANI F P, RADZISZEWSKI P H, et al. Development of a dynamic model for a cone bolt[J]. International Journal of Rock Mechanics & Mining Sciences, 2009, 46(1): 107-114.

    [8]

    CHARETTE F, PLOUFFE M. A new rock bolt concept for under ground excavations under high stress conditions[C]//Proceedings of the 6thInternational Symposium on Ground Support in Mining and Civil Engineering Construction, 2008, Johannesburg.

    [9]

    LI C C. Field observations of rock bolts in high stress rock masses[J]. Rock Mechanics and Rock Engineering, 2010, 43(4): 491-496. doi: 10.1007/s00603-009-0067-8

    [10]

    HE M, SOUSA L R E. Experiments on rock burst and its control[C]//Australasian Ground Control in Mining Conference, 2014, Sidney.

    [11] 何满潮. 滑坡地质灾害远程监测预报系统及其工程应用[J]. 岩石力学与工程学报, 2009, 28(6): 1081-1090. doi: 10.3321/j.issn:1000-6915.2009.06.001

    HE Man-chao. Real-time remote monitoring and forecasting system for geological disasters of landslides and its engineering application[J]. Chinese Journal of Rock Mechanics and Engineering, 2009, 28(6): 1081-1090. (in Chinese) doi: 10.3321/j.issn:1000-6915.2009.06.001

    [12]

    HE M, TAO Z, ZHANG B. Application of remote monitoring technology in landslides in the Luoshan mining area[J]. Mining Science and Technology, 2009, 19(5): 609-614. doi: 10.3969/j.issn.2095-2686.2009.05.011

    [13] 何满潮, 李晨, 宫伟力, 等. NPR锚杆/索支护原理及大变形控制技术[J]. 岩石力学与工程学报, 2016, 35(8): 1513-1529. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201608001.htm

    HE Man-chao, LI Chen, GONG Wei-li, et al. Support principles of NPR bolts/cables and control techniques of large deformation[J]. Chinese Journal of Rock Mechanics and Engineering, 2016, 35(8): 1513-1529. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201608001.htm

    [14] 何满潮, 郭志飚. 恒阻大变形锚杆力学特性及其工程应用[J]. 岩石力学与工程学报, 2014, 33(7): 1297-1308. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201407001.htm

    HE Man-chao, GUO Zhi-biao. Mechanical property and engineering application of anchor bolt with constant resistance and large deformation[J]. Chinese Journal of Rock Mechanics and Engineering, 2014, 33(7): 1297-1308. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201407001.htm

    [15]

    SUN X M, ZHANG Y, WANG D, et al. Mechanical properties and supporting effect of CRLD bolts under static pull test conditions[J]. International Journal of Minerals, Metallurgy, and Materials, 2017, 24(1): 1-9. doi: 10.1007/s12613-017-1372-y

    [16]

    HE M, LI C, GONG W, et al. Dynamic tests for a Constant-Resistance-Large-Deformation bolt using a modified SHTB system[J]. Tunnelling & Underground Space Technology, 2017, 64(5): 103-116.

    [17]

    Itasca . FLAC3D User Manual, Version 6.0[J]. Itasca Consulting Group Inc., USA2017.

    [18]

    LI Z H, JIANG Y J, TAO Z G, et al. Monitoring prediction of a rockslide in an open-pit mine and numerical analysis using a material instability criterion[J]. Bulletin of Engineering Geology and the Environment, 2019, 78(3): 2041-2053. doi: 10.1007/s10064-017-1224-z

    [19] 陶志刚, 李海鹏, 孙光林, 等. 基于恒阻大变形锚索的滑坡监测预警系统研发及应用[J]. 岩土力学, 2015, 36(10): 3032-3040. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201510046.htm

    TAO Zhi-gang, LI Hai-peng, SUN Guang-lin, et al. Development of monitoring and early warning system for landslides based on constant resistance and large deformation anchor cable and its application[J]. Rock & Soil Mechanics, 2015, 36(10): 3032-3040. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201510046.htm

    [20]

    ADHIKARY D P, DYSKIN A V. Modelling of progressive and instantaneous failures of foliated rock slopes[J]. Rock Mechanics and Rock Engineering, 2007, 40(4): 349-362. doi: 10.1007/s00603-006-0085-8

    [21] 杨有成, 李群, 陈新泽, 等. 对强度折减法若干问题的讨论[J]. 岩土力学, 2008, 29(4): 1103-1106. doi: 10.3969/j.issn.1000-7598.2008.04.047

    YANG You-cheng, LI Qun, CHEN Xin-ze, et al. Discussion on strength reduction using FLAC[J]. Rock & Soil Mechanics, 2008, 29(4): 1103-1106. (in Chinese) doi: 10.3969/j.issn.1000-7598.2008.04.047

    [22]

    PRUNIER F, CHOMETTE B, BRUN M, et al. Designing geotechnical structures with a proper stability criterion as a safety factor[J]. Computers and Geotechnics, 2016, 71(5): 98-114.

    [23]

    HILL R. A general theory of uniqueness and stability in elastic-plastic solids[J]. Journal of the Mechanics and Physics of Solids, 1958, 6(3): 236-249. doi: 10.1016/0022-5096(58)90029-2

    [24]

    MERRIEN-SOUKATCHOFF V, DURIEZ J, GASC-BARBIER M, et al. Mechanical Stability Analyses of Fractured Rock Slopes[M]//Rockfall Engineering, New York: John Wiley & Sons, Inc, 2013: 67-112.

    [25]

    NICOT F, LERBET J, DARVE F. Second-order work criterion: from material point to boundary value problems[J]. Acta Mechanica, 2017, 228(7): 2483-2498. doi: 10.1007/s00707-017-1844-1

    [26]

    NICOT F, SIBILLE L, DONZE F, et al. From microscopic to macroscopic second-order work in granular assemblies[J]. Mechanics of Materials, 2007, 39(7): 664-684. doi: 10.1016/j.mechmat.2006.10.003

    [27]

    NICOT F, HADDA N, BOURRIER F, et al. Inertia effects as a possible missing link between micro and macro second-order work in granular media[J]. International Journal of Solids and Structures, 2012, 49(10): 1252-1258. doi: 10.1016/j.ijsolstr.2012.02.005

    [28]

    NICOT F, DARVE F, KHOA H D V. Bifurcation and second-order work in geomaterials[J]. International Journal for Numerical and Analytical Methods in Geomechanics, 2007, 31(8): 1007-1032.

    [29] 黄茂松, 曲勰, 吕玺琳. 基于状态相关本构模型的松砂静态液化失稳数值分析[J]. 岩石力学与工程学报, 2014, 33(7): 1479-1487. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201407021.htm

    HUANG Mao-song, QU Xie, LÜ Xi-lin. Instability and static liquefaction analysis of loose sands with a state-dependent constitutive model[J]. Chinese Journal of Rock Mechanics & Engineering, 2014, 33(7): 1479-1487. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201407021.htm

    [30] 吕玺琳, 赖海波, 黄茂松. 饱和土体静态液化失稳理论预测[J]. 岩土力学, 2014(5): 1329-1333, 1339. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201405018.htm

    LÜ Xi-lin, LAI Hai-bo, HUANG Mao-song. Theoretically predicting instability of static liquefaction of saturated soils[J]. Rock & Soil Mechanics, 2014(5): 1329-1333, 1339. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201405018.htm

    [31] 吕玺琳, 钱建固, 黄茂松. 不排水加载条件下K0固结饱和砂土失稳预测[J]. 岩土工程学报, 2015, 37(6): 1010-1015. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201506007.htm

    LÜ Xi-lin, QIAN Jian-gu, HUANG Mao-song. Prediction of instability of K0-consolidated saturated sands under undrained loading conditions[J]. Chinese Journal of Geotechnical Engineering, 2015, 37(6): 1010-1015. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201506007.htm

    [32] 刘洋, 樊猛, 晏洲毅. 常偏应力剪切条件下砂土失稳模式的离散元模拟[J]. 岩土工程学报, 2020, 42(3): 467-475. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202003011.htm

    LIU Yang, FAN Meng, YAN Zhou-yi. DEM simulation of instability mode in sand under constant shear drained conditions[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(3): 467-475. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC202003011.htm

    [33]

    LI Z, DUFOUR F, DARVE F. Hydro-elasto-plastic modelling with a solid/fluid transition[J]. Computers & Geotechnics, 2016, 75(2): 69-79.

    [34]

    LI Z H, JIANG Y J, LV Q, et al. Consistent modeling of a catastrophic flowslide at the Shenzhen landfill using a hydro-elasto-plastic model with solid-fluid transition[J]. Acta Geotechnica, 2018, 13(6): 1451-1466.

    [35] 王立忠, 舒恒. Hill稳定条件在有限元法计算地基承载力中的应用[J]. 岩石力学与工程学报, 2010, 29(增刊1): 3122-3131. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2010S1079.htm

    WANG Li-zhong, SHU Heng. Application of hill's stability condition to bearing capacity computation of foundation with finite element method[J]. Chinese Journal of Rock Mechanics & Engineering, 2010, 29(7): 3122-3131. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX2010S1079.htm

    [36]

    NICOT F, DAOUADJI A, LAOUAFA F, et al. Second-order work, kinetic energy and diffuse failure in granular materials[J]. Granular Matter, 2011, 13(1): 19-28.

    [37]

    NICOT F, DARVE F. Failure in rate-independent granular materials as a bifurcation toward a dynamic regime[J]. International Journal of Plasticity, 2012, 29(1): 136-154.

  • 期刊类型引用(14)

    1. 姜淑印,李向阳,杨超,尹磊建,王佳奇,朱利勇. 考虑析水效应的PPGF浆液扩散规律与抗分散特征. 金属矿山. 2025(04): 43-53 . 百度学术
    2. 蔡跃辉. 动水注浆堵漏技术研究现状与发展情况. 科技创新与应用. 2025(13): 177-180 . 百度学术
    3. 李海燕,夏茂哲,张锟,张波,孙怀凤,赵国东,韩俊飞,刘功杰,贺恩磊. 岩溶凹陷式露天矿山大流量涌水治理技术. 煤炭科学技术. 2024(01): 267-279 . 百度学术
    4. 林久卿,牛昊,刘致延,李晓亮,王彦哲,李召峰,陈经棚. 水泥基矽土注浆材料抗海水侵蚀性能研究. 防灾减灾工程学报. 2024(03): 551-559 . 百度学术
    5. 付贵永,肖杨,史金权,周航,刘汉龙. 干湿循环下EICP联合黄原胶加固钙质粉土劣化特性试验研究. 岩土工程学报. 2024(11): 2341-2351 . 本站查看
    6. 陈亮,孙晨,王雷雨,邵晓妹,胡靖宇. 引水隧洞超前预处理灌浆材料研究与应用进展. 南水北调与水利科技(中英文). 2024(06): 1181-1188 . 百度学术
    7. 雷华阳,施福硕,刘旭,崔溦. 砂性地层中植物胶改性泥浆性质及渗透成膜试验研究. 岩土工程学报. 2023(02): 394-401 . 本站查看
    8. 张胜杰,王鸥,王天亮,王林,刘松松. 黄原胶及瓜尔胶改良尾矿砂强度特性及微观机制. 工程地质学报. 2023(02): 441-448 . 百度学术
    9. 周中,邓卓湘,鄢海涛,张俊杰. 岩溶区隧道新型绿色注浆材料试验研究. 铁道工程学报. 2023(07): 63-68 . 百度学术
    10. 吴龙骥,吴志军,翁磊. 聚丙烯酸酯改性水泥注浆材料力学性能与微观结构研究. 力学与实践. 2023(05): 999-1009 . 百度学术
    11. 夏冲,李传贵,冯啸,赵宏魁,张思峰,武剑峰. 水泥粉煤灰-改性水玻璃注浆材料试验研究与应用. 山东大学学报(工学版). 2022(01): 66-73+84 . 百度学术
    12. 付宏渊,查焕奕,潘浩强,曾铃,刘杰. 生物聚合物改良预崩解炭质泥岩水稳性及冲刷试验研究. 中南大学学报(自然科学版). 2022(07): 2633-2644 . 百度学术
    13. 张昊,胡相明,王伟,梁运涛,王兆喜,刘金举,白光星,赵艳云,吴明跃. 黄原胶和氧化镁改性黏土-水泥基新型喷涂堵漏风材料的制备及特征. 煤炭学报. 2021(06): 1768-1780 . 百度学术
    14. 康正斌,李小强,巩越. 强渗透涌水地层注浆新材料的配制与工程特性研究. 新型建筑材料. 2021(12): 19-23 . 百度学术

    其他类型引用(16)

图(10)  /  表(2)
计量
  • 文章访问数:  344
  • HTML全文浏览量:  49
  • PDF下载量:  121
  • 被引次数: 30
出版历程
  • 收稿日期:  2020-08-10
  • 网络出版日期:  2022-12-02
  • 刊出日期:  2021-09-30

目录

/

返回文章
返回