Stress-dependent properties of physicochemical interaction of unsaturated saline clay
-
摘要: 为了揭示物理化学作用对含盐溶液非饱和黏土压缩行为的影响及其与应力水平的依赖特性，对孔隙含蒸馏水、氯化钠和硫酸钠溶液及控制基质吸力条件下的试样进行了一维压缩试验。然后，根据试验结果测定了不同条件下黏土的压缩指数、次固结系数和屈服应力，并标定了它们随基质吸力和渗透吸力的变化规律。此外，对不同基质吸力和渗透吸力下非饱和盐渍黏土主、次固结行为展开了深入分析，明确了非饱和盐渍黏土物理化学作用的应力依赖特征，并对非饱和盐渍黏土的LC屈服行为进行了探索。结果表明：不同物理化学力下非饱和盐渍黏土的次固结系数与压缩指数的比值以及屈服应力可以使用基质吸力和渗透吸力来统一描述；由塑性加载区压缩曲线斜率以及特征参量、与竖向应力的相关性指出非饱和盐渍黏土物理化学作用与应力水平息息相关；此外，化学-水力-力学耦合作用下非饱和盐渍黏土的LC屈服曲线是由MLC屈服曲线和OLC屈服曲线组成的一条光滑曲线。Abstract: To reveal the impact of physicochemical effects on the compressive behaviors of unsaturated clay containing salt solution and its dependence properties on stress level, one-dimensional compression tests are performed on the specimens with pores containing distilled water, sodium chloride solution, sodium sulfate solution and controlled matric suction conditions. Then, the compression index, secondary compression coefficient and yield stress of clay under different conditions are measured according to the test results, and their variation laws with the matric suction and osmotic suction are calibrated. Furthermore, the stress-dependent characteristics of physicochemical action are clarified through an in-depth analysis of the primary and secondary consolidation behaviors of unsaturated saline clay under different matric suctions and osmotic suctions, and the LC yielding behaviors of unsaturated saline clay are explored. The results show that the ratio of the secondary compression coefficient to the compression index and the yield stress of unsaturated saline clay at different physicochemical forces can be described uniformly using the osmotic suction and matric suction. From the slope of compression curve in the plastic loading zone and the correlation between characteristic parameter, and the vertical stress, it is noted that the physicochemical action of unsaturated saline clay is closely related to the stress level. Moreover, the LC yield curve of unsaturated saline clay is a smooth curve composed of MLC yield curve and OLC yield curve under chemo- hydro-mechanical coupling.
-
0. 引言
随着中国煤炭资源开采大范围向深部转移，冲击矿压（也称“冲击地压”）灾害日益凸显，严重影响煤矿安全生产[1]。冲击矿压主要是由采矿活动引起的煤岩结构局部应力集中引起，大量研究表明，在应力异常集中区（如断层、褶曲、煤柱区等），更容易发生冲击矿压[2]。
此外，相变构造区似乎也会增加煤层的异常应力集中[3]。国内外学者针对相变构造区诱发冲击矿压的研究，大多集中在煤层厚度变化诱发冲击矿压方面，针对坚硬覆岩厚度变化导致应力异常及诱冲效应的研究较少。孙振武[4]根据地应力场测量和矿压观测结果，发现在局部煤层厚度变化的区域往往会发生构造应力异常现象；Álvarez-Fernández等[5]运用数值模拟分析了分段崩落法开采煤层时，由于煤层厚度变薄导致在某超前巷道应力显著增大、变形较大的具体情况；Zhu等[6]通过现场观测和数值分析，研究了煤厚变化区初始应力场及采动应力场的分布规律；南存全等[7]研究发现，当工作面接近煤厚变异区时，煤厚变异区对工作面超前支承压力分布特征影响显著；王勇等[8]分析了不同煤层厚度变化条件下应力及能量分布规律；赵同彬等[9]研究了煤层厚度变化对超前支承压力及能量演化规律的影响，揭示了煤厚变异区冲击矿压发生的力学机制。
关于坚硬覆岩厚度变化与煤层厚度变化诱发冲击矿压机理的问题，其实质都是在厚度变化区产生局部应力集中，但是覆岩厚度变化不像煤厚变化能随着工作面的推进直观表现出来，只能从部分钻孔了解覆岩的分布情况。因此，坚硬覆岩厚度变化诱发冲击矿压的问题很难引起研究人员的注意。
针对现有研究成果的局限性，本文针对中国内蒙深部矿区坚硬覆岩厚度变化区煤层开采诱发冲击矿压机理问题，采用弹性力学理论对坚硬覆岩厚度变化区的构造应力特征进行力学分析，试图揭示覆岩厚度变化区应力异常分布的成因机制；并利用FLAC3D数值模拟研究覆岩厚度变化区工作面超前支承压力分布特征以及煤壁前方的能量演化规律，进一步探讨坚硬覆岩厚度变化区煤层开采诱发冲击矿压的机理。
1. 覆岩厚度变化区煤层应力异常形成机制
当煤岩体的应力不超过弹性范围时，最适宜用弹性力学的方法解析煤岩地下工程的力学行为[10]。原岩应力条件下，煤层及顶底板未受采掘扰动影响，假设煤岩体均未超过其弹性范围，可将覆岩厚度变化区简化为覆岩1、覆岩2两个组合的弹性元件[9]。其整体力学属性由两个弹性元件串联、并联组合表示，模型如图 1所示。
由图 1可知，在覆岩厚度变化区，区域一、区域二、区域三的弹性模量均为覆岩1、覆岩2的弹性模量串联组成，其等效弹性模量E为
1E=H1H⋅ER1+H2H⋅ER2， (1) 式中：H1，H2分别为区域一覆岩1、覆岩2的厚度；ER1，ER2分别为覆岩1、覆岩2的弹性模量。
根据煤岩体应力应变关系可知：
σ=Eε。 (2) 本模型内区域一、区域二、区域三相互并联，故区域一、区域二、区域三的应变εa，εb，εc相等；同时串联部分应力相等，联立式（1），式（2）可得覆岩变化区下方煤层应力关系为
σaσc=E1E3=H′2⋅ER1+H′1⋅ER2H2⋅ER1+H1⋅ER2， (3) 式中：σa，σc为区域一、区域三覆岩厚度变化区下方的应力，H′1，H′2为区域三覆岩1、覆岩2的厚度。
为了定性分析覆岩厚度变化区下方的应力分布情况，假设区域一覆岩1和覆岩2的厚度H1和H2保持不变，将式（3）优化为ER2/ER1和H1'/H1的函数：
σaσc=[HH1 + H′1H1⋅(ER2ER1−1)]/[HH1 + (ER2ER1−1)]。 (4) 对于坚硬顶板赋存条件ER1＞ER2，区域一坚硬覆岩厚度H1大于区域三坚硬覆岩厚度H′1。由式（4）可知，区域一所受的应力σa大于区域三所受的应力σc，即覆岩变厚区的构造应力比覆岩变薄区大。
为研究ER2/ER1与H′1/H1变化对σa/σc的影响，控制式（4）中的自变量，将H/H1设为定值；同时ER2/ER1与H′1/H1的数值仅代表覆岩性质及厚度变化的一个趋势，不做实际情况下精确数值的讨论。图 2（a）为覆岩厚度比值一定的情况下，构造应力比值与覆岩弹性模量比值关系图，H1'/H1一定，σa/σc随ER2/ER1增加呈反比例关系，即两覆岩性质相差越大，构造应力变化越大。图 2（b）为覆岩性质一定的情况下，构造应力比值与覆岩厚度变化关系图，当ER2/ER1一定，σa/σc随H1'/H1增加呈线性递减关系，说明坚硬覆岩厚度变化越大，构造应力变化越大。
![]()
图 2 构造应力比值与覆岩性质/厚度变化关系图Figure 2. Relationship between tectonic stress ratio and roof properties or variation of roof thickness2. 覆岩厚度变化区煤层开采应力与能量分布特征
2.1 覆岩厚度变化区煤层开采数值模型建立
为研究坚硬覆岩厚度及性质变化对煤层应力分布特征的影响，本文以某矿煤层基本顶的厚度及性质变化为例，利用FLAC3D建立覆岩厚度变化区煤层开采数值计算模型。模型尺寸x×y×z为300 m×400 m ×150 m，覆岩变化区为y轴中部100 m范围，其它模型参数如图 3所示。初始应力状态下，假设煤岩均处于线弹性状态，模型本构设置为Elastic弹性模型[11]。
许家林等[12]基于岩层控制关键层理论指出，只有关键层上部的载荷可以简化为均布载荷。因此，本文在模型顶部设置一层20 m厚的关键层；模拟煤层埋深约700 m，在模型顶部施加15 MPa的均布载荷，重力加速度为10 m/s2；底部固定x，y，z方向的位移；x，y方向为固定边界，并在x，y方向施加水平渐变应力，x，y方向应力系数分别设为1.0和1.5。
模型材料参数以某矿实际的煤岩力学参数为基准，对材料参数进行简化处理[13]，如表 1所示。模拟方案主要为：方案一：模拟H1=60 m，H′1=30 m，ER2=10 GPa，ER1为15，20，25，30，35 GPa 5种情况下，覆岩性质变化对煤层初始应力分布特征的影响；方案二：模拟ER1=30 GPa，ER2=10 GPa，H1=60 m，H1'为10，20，30，40，50 m 5种情况下，覆岩厚度变化对煤层初始应力分布特征的影响。
表 1 模型选用参数Table 1. Model parameters序号 岩性 厚度/m 密度/(kg·m-3) 弹性模量/GPa 泊松比 1 细粒砂岩 20 2400 20 0.20 2 中粒砂岩 30* 2500 30* 0.15 砂质泥岩 30* 2200 10 0.25 3 砂质泥岩 10 2200 10 0.25 4 细粒砂岩 20 2400 20 0.20 5 煤层 6 1400 5 0.30 6 砂质泥岩 14 2200 10 0.25 7 细粒砂岩 20 2400 20 0.20 注：表中带“*”参数为方案一、方案二中的可变量。 2.2 覆岩厚度变化区煤层初始应力分布特征
图 4所示为不同方案下煤层内初始应力分布规律，由图 4（a）可知，由覆岩较薄区至覆岩较厚区，煤层内初始应力由均布状态逐渐减小，然后随着覆岩厚度的增加而增加，当处于覆岩较厚区时，初始应力又逐渐减小，最终恢复至均布状态；并且最大（最小）应力值随覆岩弹性模量比的增加而增加（减小）。由图 4（b）可知，初始应力分布规律与方案一相似，最大（最小）应力值随覆岩厚度变化率的增加而增加（减小）。
由于在模型顶部施加了15 MPa的均布载荷，因此，坚硬覆岩性质及其厚度变化对煤层初始应力的影响程度为煤层内的应力值减模型顶部均布载荷之后的比值，如图 5所示。方案一中H′1/H1不变，覆岩弹性模量ER2/ER1由10/35增加至10/15时，煤层内应力变化比值呈反比例函数由1.70降到1.21；方案二中ER2/ER1不变，覆岩厚度情况H′1/H1由1/6增加至5/6，煤层内应力变化比值呈线性关系由2.14递减至1.22。这与理论分析式（4）所得规律一致。
![]()
图 5 覆岩性质及厚度变化对应力变化率的影响Figure 5. Influences of properties and thickness of overburden on change of stress rate2.3 覆岩厚度变化区煤层采动应力分布特征
煤岩体所受应力在采动应力影响下会超过其峰值应力，煤岩体进入全应力应变曲线的峰后阶段，而应变软化模型能很好的反映煤岩体的峰后破坏阶段[14]。因此，将工作面回采过程中的煤岩体本构模型设置为应变软化模型，煤岩体应变软化参数的设置参考文献[14]，如表 2所示，其余模型参数见表 1。
表 2 煤岩体应变软化参数Table 2. Strain-softening parameters of coal and rock塑性应变 中粒砂岩 砂质砂岩 黏聚力/MPa 摩擦角/(°) 剪胀角/(°) 黏聚力/MPa 摩擦角/(°) 剪胀角/(°) 0 10.00 36 18 5.00 28 12 1×10-4 6.00 32 10 3.00 24 6 2×10-4 4.00 28 5 2.00 22 2 5×10-4 2.00 24 0 1.00 20 0 1 2.00 24 0 1.00 20 0 塑性应变 细粒砂岩 煤层 黏聚力/MPa 摩擦角/(°) 剪胀角/(°) 黏聚力/MPa 摩擦角/(°) 剪胀角/(°) 0 8.00 32 16 4.00 25 10 1×10-4 5.00 28 8 2.00 22 5 2×10-4 3.00 24 4 1.00 20 2 5×10-4 1.50 22 0 0.50 18 0 1 1.50 22 0 0.50 18 0 为研究覆岩厚度变化对工作面回采的影响，采用的开挖模型为：H1=H1'=60 m，覆岩厚度不变；H1= 2H′1=60 m，覆岩厚度发生改变。煤层回采方式包括沿覆岩由厚向薄推进和沿覆岩由薄向厚推进，如图 6所示。为消除边界效应对工作面回采的影响，在工作面两侧留50 m的边界煤柱。此外，工作面每开挖20 m，计算一次平衡。
工作面回采过程中，煤岩体内应力重新分布并在工作面超前区域形成应力集中区。为研究覆岩厚度变化对工作面回采超前支承压力的影响，在煤层中部，x=150 m，z=37 m处布置一条测线，将工作面沿覆岩由薄向厚回采、由厚向薄回采与沿覆岩厚度不变回采的超前支承压力进行对比分析。定义应力突变系数为工作面沿覆岩厚度变化推进的超前支承压力与沿覆岩厚度不变推进的超前支承压力之比。限于篇幅，本文选用超前工作面10，20，30，40，50 m处的支承压力进行对比分析，如图 7所示。
![]()
图 7 不同推进方式工作面超前支承压力突变规律Figure 7. Abrupt change of abutment pressure in advance of working faces with different mining modes由图 7（a）可以看出，当工作面沿覆岩由薄向厚推进时，应力突变系数呈先减小后升高再减小的趋势。工作面在覆岩较厚区，应力突变系数最大，平均为1.042；在覆岩较薄区，应力突变系数最小为0.947，同时超前距离越近，应力突变系数越接近1，表明距离回采位置越近，受超前支承压力影响越大，受覆岩厚度变化影响较小。
由图 7（b）可得，当工作面沿覆岩由厚向薄推进时，超前工作面的应力突变系数有先增加后降低再升高的趋势。工作面在覆岩较厚区，应力突变系数最大约为1.039；在覆岩较薄区，应力突变系数最小为0.943，同时超前距离越近，受超前支承压力影响越大，受覆岩厚度变化影响较小。
通过分析可知，应力突变系数最大的区域，均在覆岩较厚区至覆岩变化区范围，当工作面推进至此区域时，超前支承压力与突变的构造应力进行叠加，增加了冲击矿压发生的可能性。因此，工作面在坚硬覆岩厚度变化区推进时，在覆岩较厚区至覆岩变化区应力集中系数最大，当工作面推进至此区域时，因重点关注此区域的应力异常情况，必要时进行卸压处理。
2.4 覆岩厚度变化区煤层开采能量分布特征
谢和平等[15]指出，处于三向应力状态下的工程岩体常常面临卸载破坏的危险，并且煤岩体单元在卸载过程中释放的弹性能密度为
WE=σ21+σ22+σ23−2ν⋅σ1⋅σ2−2ν⋅σ2⋅σ3−2ν⋅σ1⋅σ32E， (5) 式中：σ1，σ2，σ3为最大、中间、最小主应力；E为弹性模量；ν为泊松比。
为研究覆岩厚度变化对工作面回采煤壁前方积聚弹性应变能演化规律的影响，由式（5）得出了工作面回采过程中模型内储存的弹性能密度。为了更直观地反映工作面回采过程中煤层内的能量聚积情况，将工作面回采过程中每一步的弹性能密度云图进行切片（切片取z=37 m）、组合处理，如图 8所示。
![]()
图 8 不同推进方式下煤层弹性能密度切片云图Figure 8. Section cloud diagram of elastic energy density in different mining modes工作面沿覆岩由薄向厚推进时，除了在工作面开挖附近出现了能量聚积区，还在覆岩较厚区出现第二个能量聚积区；随工作面继续向覆岩厚度变化区推进，能量聚集区范围增加，冲击危险性增加。工作面沿覆岩由厚向薄推进时，仅在工作面前方出现一个能量聚集区，该能量聚集区在覆岩较厚区至覆岩变化区范围较大，冲击危险性更高；随工作面继续推进，能量聚集区范围逐渐减小。
3. 覆岩厚度变化区煤层开采诱发冲击矿压机理
冲击矿压发生的“动静载叠加诱冲原理”表示，采掘空间周围煤岩体中的静载荷与矿震形成的动载荷叠加超过煤岩体发生冲击的最小载荷时，就容易诱发冲击矿压[16-17]，即
σs+σd⩾σbmin， (6) 式中：σs为煤岩体静载应力，σd为矿震动载，σbmin为冲击矿压的临界应力。
根据覆岩厚度变化区煤层开采过程中超前支承压力分布特征及动静载叠加诱冲原理，得出了不同回采方式下覆岩厚度变化区煤层开采诱发冲击矿压的机理，如图 9所示。
![]()
图 9 冲击矿压动静载叠加诱发机理示意图Figure 9. Schematic diagram of induced mechanism of dynamic and static load superposition under rock burst pressure当工作面沿覆岩由薄向厚推进时，超前支承压力由变薄区向变厚区传递，煤层中的静载应力相比覆岩厚度不变时有先减小后升高的趋势，在覆岩较厚区会形成应力的第二峰值区。当工作面前方覆岩破断诱发动载时，动载荷与静载荷相互叠加，可能会出现两个冲击矿压潜在发生区，即：工作面前方超前支承压力峰值区和应力第二峰值区与动载荷叠加诱冲区。其中第二峰值区与动载荷叠加诱冲区位于工作面前方，产生的冲击能量将主要向工作面两巷或者工作面临空面释放，冲击矿压影响范围增加。
当工作面沿覆岩由厚向薄推进时，超前支承压力由变厚区向变薄区传递，煤层中的静载应力相比覆岩厚度不变时先升高后减小，仅在超前支承压力影响范围与覆岩较厚区形成一个峰值应力；与动载荷叠加时，形成一个冲击矿压潜在发生区，但是该区域相对覆岩厚度不变时影响范围更大、冲击能量更高。
因此，工作面在覆岩较厚区至覆岩厚度变化区回采时，煤壁前方超前支承压力与突变的构造应力进行叠加，积聚的弹性能增加，更容易诱发冲击矿压。冲击矿压防治的主要工作是降低覆岩变化区第二峰值应力，将应力集中区向煤岩体深部转移。
4. 覆岩厚度变化区工程案例分析
4.1 工作面沿覆岩由薄向厚推进案例分析
（1）工程概况
伊泰某矿103工作面为矿井一采区的第2个回采工作面，东南侧为实体煤；西北侧为101采空区，两者留有两个30 m宽的双煤柱；东北邻南翼辅运大巷；西南方向开切眼为实体煤柱。工作面设计推进长度2480 m，宽度210 m，主采3-1煤层，煤层埋深平均706.3 m，煤层厚度约6.7 m，倾角1°~3°，平均倾角2°，为近水平煤层。
统计103工作面附近的钻孔，并绘制出覆岩厚度变化等值线图，发现在3-1煤层上方近200 m范围主要存在一层较厚的坚硬顶板，该覆岩位于煤层上方20 m左右，层厚在60~110 m；经计算可知，该岩层为工作面覆岩结构的关键层。随着工作面的推进，该覆岩由薄变厚，如图 10所示。
![]()
图 10 伊泰某矿103工作面覆岩厚度等值线及微震平面定位图Figure 10. Contours of roof thickness and distribution map of microseismic plane location of working face 103 in Yitai Mine（2）案例分析
103工作面自2018年3月开始回采以来，大能量事件频发，103辅运顺槽出现了不同程度的动力显现，造成巷道多次破坏和设备损坏，破坏以巷道底板的瞬间开裂、底鼓为主。通过分析微震监测数据，可判断煤岩体的应力分布状态及覆岩运移规律[18]。统计分析103工作面回采期间的微震数据与坚硬覆岩厚度变化的关系，以103临空顺槽做剖面，如图 11所示。
![]()
图 11 沿覆岩由薄向厚推进覆岩厚度变化与微震关系图Figure 11. Relationship between variation of roof thickness and microseismicity when working face advancing from thin to thick along roof103工作面回采期间，位于工作面前方的微震事件频次占比87.7%，微震能量占比91.0%，表明103工作面回采期间微震事件与工作面前方煤岩体应力分布状态关系较大。同时，微震频次、能量的分布情况与覆岩厚度变化存在一定的对应关系，工作面沿覆岩由薄向厚推进时，工作面推进位置前方微震能量、频次均处于较大的水平，工作面冲击危险性较高，此时对应于冲击矿压的潜在发生区Ⅰ，主要由工作面采动引起的超前支承压力造成；在坚硬覆岩较厚区域，微震频次及＜103 J的微震总能量存在极大值，此时对应于冲击矿压的潜在发生区Ⅱ，主要受超前支承压力与突变的构造应力叠加影响。
4.2 工作面沿覆岩由厚向薄推进案例分析
（1）工程概况
中煤某矿3102工作面为矿井第2个回采工作面，主采3-1煤层，工作面走向长5540 m，倾向宽300 m，东侧为实体煤，西侧与3101采空区相邻，南侧为井田边界，北侧为3-1煤辅助回风巷，区段煤柱宽度35 m。该工作面煤层埋深690~700 m，煤层厚度4.35~5.47 m，倾角1°~4°，平均倾角2°，为近水平煤层。
统计3102工作面附近的钻孔，并绘制出覆岩厚度变化等值线图，发现在3-1煤层上方近200 m范围主要存在一层厚度变化较大的坚硬顶板，该覆岩位于煤层上方40 m左右，层厚在40~60 m。随着工作面的推进，该覆岩由厚变薄，如图 12所示。
![]()
图 12 中煤某矿3102工作面覆岩厚度等值线及微震平面定位图Figure 12. Contours of roof thickness and distribution map of microseismic plane location of working face3102 in China Coal Mine（2）案例分析
3102工作面自2017年8月开始回采，大能量事件频发，3102回风巷出现了不同程度的帮鼓、顶板下沉、底鼓、单体损坏等冲击显现情形。统计分析3102工作面回采期间所有的微震数据与坚硬厚度变化的关系，以3102临空巷做剖面，见图 13。
![]()
图 13 沿覆岩由厚向薄推进覆岩厚度变化与微震关系图Figure 13. Relationship between variation of roof thickness and microseismicity when working face advancing from thick to thin along roof3102工作面回采期间，微震事件大多位于工作面前方，其中微震事件频次占比62.3%，微震能量占比87.1%。工作面沿覆岩由厚向薄推进时，在坚硬覆岩较厚区至覆岩厚度变化区，微震总频次分布较多，极大值分布在坚硬覆岩厚度较大的区域；并且小于103 J的微震总能量与微震频次具有相似的分布规律。表明在覆岩较厚区至覆岩厚度变化区的冲击危险性较高，此时对应于工作面沿覆岩由厚向薄推进时的冲击矿压潜在发生区，主要受超前支承压力与突变的构造应力叠加影响。
以上分析可知，工作面沿覆岩厚度变化区回采期间，微震事件大多分布在工作面前方，与煤岩体应力分布状态关系较大；同时，微震频次及能量小于103 J的微震总能量与坚硬覆岩厚度变化相关性较好。工作面受超前支承压力与突变的构造应力叠加影响，静载应力较大，煤岩体产生较多的微裂隙，导致坚硬覆岩厚度变化区至覆岩较厚区微震事件频发。这与前文研究成果相符。
5. 结论
（1）基于弹性力学理论分析了坚硬覆岩厚度变化区的力学机制，结果表明，覆岩变厚区的构造应力比覆岩变薄区大，覆岩厚度变化越大或覆岩之间弹性模量相差越大，覆岩厚度变化区构造应力变化越大。
（2）采用数值模拟方法研究了坚硬覆岩厚度及物理力学性质变化对煤层构造应力分布特征的影响，结果表明，自覆岩较薄区到覆岩较厚区，构造应力均有先减小后增大再减小的趋势，覆岩厚度变化越大或覆岩物理力学性质差异越大，应力变化比值越大。
（3）工作面在覆岩厚度变化区开采时，受超前支承压力与突变的构造应力叠加影响，覆岩厚度变化区至覆岩较厚区应力集中程度较大，该区域积聚的弹性能主要向工作面前方巷道释放，冲击危险性更大；冲击矿压防治思路是降低覆岩变化形成的第二应力峰值区。
（4）现场案例表明，工作面沿覆岩由薄向厚推进或沿覆岩由厚向薄推进时，在坚硬覆岩厚度变化区及变化区向较厚区过渡时微震能量、频次较高，冲击矿压危险上升，与理论与模拟分析较为吻合。
-
![]()
图 1 非饱和盐渍黏土的一维压缩试验结果
Figure 1. One-dimensional compression test results of unsaturated saline clay
![]()
图 2 非饱和盐渍黏土在压缩过程中体积应变的变化情况
Figure 2. Change in volume strain of unsaturated saline clay during compression process
![]()
图 3 不同含水和含盐条件下试样的吸力
Figure 3. Suctions of samples under different water and salt contents
![]()
图 4 非饱和盐渍黏土物理化学作用应力依赖性的示意图
Figure 4. Schematic of stress dependence of physicochemical interaction of unsaturated saline clay
![]()
图 5 压缩指数随基质吸力和渗透吸力的变化规律
Figure 5. Variation of compression index with matric suction and osmotic suction
![]()
图 6 非饱和盐渍黏土物理化学作用的应力依赖性分析
Figure 6. Analysis of stress dependence of physicochemical interaction of unsaturated saline clay
![]()
图 7 不同应力水平下与渗透吸力和基质吸力的关系
Figure 7. Relationship among , osmotic suction and matric suction at different stress levels
![]()
图 9 屈服应力随基质吸力和渗透吸力的变化规律
Figure 9. Variation of yield stress with matric suction and osmotic suction
表 1 试验用土中易溶盐的离子种类和含量
Table 1 Ion species and contents of soluble salts in experimental soil
阳离子/% 阴离子/% 总含盐量/% pH值 0.0113 0.4470 0.0481 0.0546 — 0.0235 0.3222 0.0912 1.0 7.47 
下载: 导出CSV
表 2 不同基质吸力和渗透吸力下的次固结系数
Table 2 Secondary consolidation coefficients under different conditions of matric suction and osmotic suction
s = 0 kPa s = 50 kPa s = 100 kPa s = 200 kPa π = 457 kPa π = 457 kPa π = 5630 kPa π = 5902 kPa π = 457 kPa π = 457 kPa π = 5630 kPa π = 5902 kPa σ/kPa Ca×102 σ/kPa Ca×102 σ/kPa Ca×102 σ/kPa Ca×102 σ/kPa Ca×102 σ/kPa Ca×102 σ/kPa Ca×102 σ/kPa Ca×102 100 0.1165 100 0.1165 100 0.0764 97 0.0872 101 0.1156 97 0.1161 99 0.0792 97 0.1126 200 0.1685 200 0.1222 200 0.1650 201 0.1676 194 0.1160 194 0.1141 200 0.1451 201 0.1540 300 0.2186 300 0.1318 300 0.1555 291 0.1657 291 0.1220 291 0.1140 300 0.1669 291 0.1673 400 0.2184 401 0.1755 400 0.2141 388 0.2339 387 0.1625 388 0.1159 400 0.1256 387 0.1675 600 0.2672 601 0.1959 600 0.2644 581 0.2626 581 0.1739 581 0.1645 599 0.2098 581 0.2201 800 0.2911 802 0.2645 800 0.3076 775 0.3529 775 0.2040 775 0.1744 800 0.2511 774 0.2722 1000 0.3623 1003 0.3075 1001 0.3952 969 0.4083 969 0.2915 969 0.2510 1000 0.3028 968 0.3217
下载: 导出CSV
-
[1] 张莎莎, 谢山杰, 杨晓华, 等. 火山灰改良粗粒硫酸盐渍土路基填料及其作用机理研究[J]. 岩土工程学报, 2019, 41(3): 588-594. doi: 10.11779/CJGE201903023 ZHANG Shasha, XIE Shanjie, YANG Xiaohua, et al. Action mechanism of coarse particle sulfate soil subgrade modified by volcanic ash[J]. Chinese Journal of Geotechnical Engineering, 2019, 41(3): 588-594. (in Chinese) doi: 10.11779/CJGE201903023
[2] ZHANG Y, YE W M, CHEN B, et al. Desiccation of NaCl-contaminated soil of earthen heritages in the Site of Yar City, northwest China[J]. Applied Clay Science, 2016, 124-125: 1-10. doi: 10.1016/j.clay.2016.01.047
[3] LU N, GODT J W, WU D T. A closed‐form equation for effective stress in unsaturated soil[J]. Water Resources Research, 2010, 46(5): W05515.
[4] YAO CQ, CHEN P, MA TT, et al. Physicochemical effect on shear strength characteristics of clayey soils based on ring-shear experiment[J]. Canadian Geotechnical Journal, 2020, 57(12): 1820-1831. doi: 10.1139/cgj-2019-0513
[5] 陈云敏. 环境土工基本理论及工程应用[J]. 岩土工程学报, 2014, 36(1): 1-46. doi: 10.11779/CJGE201401001 CHEN Yunmin. A fundamental theory of environmental geotechnics and its application[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(1): 1-46. (in Chinese) doi: 10.11779/CJGE201401001
[6] THYAGARAJ T, DAS A P. Physico-chemical effects on collapse behaviour of compacted red soil[J]. Géotechnique, 2017, 67(7): 559-571. doi: 10.1680/jgeot.15.P.240
[7] 孙重初. 酸液对红粘土物理力学性质的影响[J]. 岩土工程学报, 1989, 11(4): 89-93. http://cge.nhri.cn/cn/article/id/9296 SUN Zhongchu. The effect of acidizing fluid on the physico mechanical properties of red clay[J]. Chinese Journal of Geotechnical Engineering, 1989, 11(4): 89-93. (in Chinese) http://cge.nhri.cn/cn/article/id/9296
[8] 陈声凯, 凌建明, 罗志刚. 路基土回弹模量应力依赖性分析及预估模型[J]. 土木工程学报, 2007, 40(6): 95-99, 104. https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC200706018.htm CHEN Shengkai, LING Jianming, LUO Zhigang. Stress-dependent characteristics and prediction model of the resilient modulus of subgrade soils[J]. China Civil Engineering Journal, 2007, 40(6): 95-99, 104. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-TMGC200706018.htm
[9] 弋晓明, 李术才, 王松根, 等. 非饱和粉土回弹模量的应力依赖性与水敏感性耦合分析[J]. 山东大学学报(工学版), 2013, 43(2): 84-88, 110. https://www.cnki.com.cn/Article/CJFDTOTAL-SDGY201302015.htm YI Xiaoming, LI Shucai, WANG Songgen, et al. Coupling analysis of stress dependence and water sensitivity for the resilient modulus of unsaturated silt soil[J]. Journal of Shandong University (Engineering Science), 2013, 43(2): 84-88, 110. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-SDGY201302015.htm
[10] SUN D A, SHENG D C, CUI H B, et al. A density-dependent elastoplastic hydro-mechanical model for unsaturated compacted soils[J]. International Journal for Numerical & Analytical Methods in Geomechanics, 2010, 31(11): 1257-1279.
[11] XU X T, JIAN W B, WU N S, et al. Void ratio-dependent water retention model for a deformable residual clay[J]. International Journal of Geomechanics, 2020, 20(8): 04020131-1. doi: 10.1061/(ASCE)GM.1943-5622.0001773
[12] CUI Y J, DELAGE P. Yielding and plastic behaviour of an unsaturated compacted silt[J]. Géotechnique, 1996, 46(2): 291-311. doi: 10.1680/geot.1996.46.2.291
[13] RAMPINO C, MANCUSO C, VINALE F. Experimental behaviour and modelling of an unsaturated compacted soil[J]. Canadian Geotechnical Journal, 2000, 37(4): 748-763. doi: 10.1139/t00-004
[14] 高登辉, 陈正汉, 郭楠, 等. 干密度和基质吸力对重塑非饱和黄土变形与强度特性的影响[J]. 岩石力学与工程学报, 2017, 36(3): 736-744. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201703023.htm GUO Denghui, CHEN Zhnghan, HUO Nan, et al. The influence of dry density and matric suction on the deformation and the strength characteristics of the remolded unsaturated loess soils[J]. Chinese Journal of Rock Mechanics and Engineering, 2017, 36(3): 736-744. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX201703023.htm
[15] DI MAIO C, SANTOLI L, SCHIAVONE P. Volume change behaviour of clays: the influence of mineral composition, pore fluid composition and stress state[J]. Mechanics of Materials, 2004, 36(5/6): 435-451.
[16] 张彤炜, 邓永锋, 刘松玉, 等. 渗透吸力对重塑黏土的压缩和渗透特性影响的试验研究[J]. 岩土工程学报, 2014, 36(12): 2260-2266. doi: 10.11779/CJGE201412014 ZHANG Tongwei, DENG Yongfeng, LIU Songyu, et al. Experimental investigation of osmotic suction effect on hydro-mechanical behaviour of remolded clay[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(12): 2260-2266. (in Chinese) doi: 10.11779/CJGE201412014
[17] WITTEVEEN P, FERRARI A, LALOUI L. An experimental and constitutive investigation on the chemo-mechanical behaviour of a clay[J]. Géotechnique, 2013, 63(3): 244-255. doi: 10.1680/geot.SIP13.P.027
[18] 颜荣涛, 赵续月, 于明波, 等. 盐溶液饱和黏土的等向压缩特性[J]. 岩土力学, 2018, 39(1): 129-138. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201801017.htm YAN Rongtao, ZHAO Xuyue, YU Mingbo, et al. Isotropic compression characteristics of clayey soil saturated by salty solution[J]. Rock and Soil Mechanics, 2018, 39(1): 129-138. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX201801017.htm
[19] 周凤玺, 王立业, 赖远明. 饱和盐渍土的一维蠕变试验与模型研究[J]. 岩土工程学报, 2020, 42(1): 142-149. doi: 10.11779/CJGE202001016 ZHOU Fengxi, WANG Liye, LAI Yuanming. One-dimensional creep tests and model studies on saturated saline soil[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(1): 142-149. (in Chinese) doi: 10.11779/CJGE202001016
[20] 王立业, 周凤玺, 秦虎. 饱和盐渍土分数阶蠕变模型及试验研究[J]. 岩土力学, 2020, 41(2): 543-551. https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202002026.htm WANG Liye, ZHOU Fengxi, QIN Hu. Fractional creep model and experimental study of saturated saline soil[J]. Rock and Soil Mechanics, 2020, 41(2): 543-551. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTLX202002026.htm
[21] ALONSO E E, GENS A, JOSA A. A constitutive model for partially saturated soils[J]. Géotechnique, 1990, 40(3): 405-430. doi: 10.1680/geot.1990.40.3.405
[22] 陈正汉, 扈胜霞, 孙树国, 等. 非饱和土固结仪和直剪仪的研制及应用[J]. 岩土工程学报, 2004, 26(2): 161-166. http://cge.nhri.cn/cn/article/id/11362 CHEN Zhenghan, HU Shengxia, SUN Shuguo, et al. Development & application of consolidation apparatus and direct shear apparatus for unsaturated soils[J]. Chinese Journal of Geotechnical Engineering, 2004, 26(2): 161-166. (in Chinese) http://cge.nhri.cn/cn/article/id/11362
[23] 陈正汉. 重塑非饱和黄土的变形、强度、屈服和水量变化特性[J]. 岩土工程学报, 1999, 21(1): 82-90. http://cge.nhri.cn/cn/article/id/10260 CHEN Zhenghan. Deformation, strength, yield and moisture change of a remolded unsaturated loess[J]. Chinese Journal of Geotechnical Engineering, 1999, 21(1): 82-90. (in Chinese) http://cge.nhri.cn/cn/article/id/10260
[24] THYAGARAJ T, THOMAS S R, DAS A P. Physico -chemical effects on shrinkage behavior of compacted expansive clay[J]. International Journal of Geomechanics, 2017, 17(2): 06016013. doi: 10.1061/(ASCE)GM.1943-5622.0000698
[25] LEONG E, RAHARDJO H, HE L. Factors affecting the filter paper method for total and matric suction measurements[J]. Geotechnical Testing Journal, 2002, 25(3): 8198.
[26] SHENG D C. Review of fundamental principles in modelling unsaturated soil behaviour[J]. Computers and Geotechnics, 2011, 38(6): 757-776. doi: 10.1016/j.compgeo.2011.05.002
[27] RAMPINO C, MANCUSO C, VINALE F. Laboratory testing on an unsaturated soil: equipment, procedures, and first experimental results[J]. Canadian Geotechnical Journal, 1999, 36(1): 1-12. doi: 10.1139/t98-093
[28] LU N, LIKOS W J. Suction stress characteristic curve for unsaturated soil[J]. Journal of Geotechnical & Geoenvironmental Engineering, 2006, 132(2): 131-142.
[29] 周凤玺, 王立业, 赖远明. 饱和盐渍土渗透吸力的回顾及研究[J]. 岩土工程学报, 2020, 42(7): 1199-1210. doi: 10.11779/CJGE202007003 ZHOU Fengxi, WANG Liye, LAI Yuanming. Review and research on osmotic suction of saturated saline soils[J]. Chinese Journal of Geotechnical Engineering, 2020, 42(7): 1199-1210. (in Chinese) doi: 10.11779/CJGE202007003
[30] PHILIP J R. Unitary approach to capillary condensation and adsorption[J]. Journal of Chemical Physics, 1977, 66(11): 5069-5075. doi: 10.1063/1.433814
[31] LU N, ZHANG C. Soil Sorptive Potential: Concept, Theory, and Verification[J]. Journal of Geotechnical and Geoenvironmental Engineering, 2018, 145(4): 04019006.
[32] 赵建斌, 申俊敏, 张艳聪. 考虑土体应力依赖性的水泥混凝土路面基层力学响应研究[J]. 公路, 2014, 59(12): 54-58. https://www.cnki.com.cn/Article/CJFDTOTAL-GLGL201412012.htm ZHAO Jianbin, SHEN Junmin, ZHANG Yancong. Research on mechanical responses of the base for cement concrete pavement considering stress dependence of subgrade[J]. Highway, 2014, 59(12): 54-58. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-GLGL201412012.htm
[33] KONG L M, YAO Y P, QI J L. Modeling the combined effect of time and temperature on normally consolidated and overconsolidated clays[J]. Acta Geotechnica, 2020, 15(3): 2451-2471.
[34] MESRI G, GODLEWSKI P M. Time and stress-compressibility interrelationship[J]. Journal of the Geotechnical Engineering Division, 1977, 103(5): 417-430. doi: 10.1061/AJGEB6.0000421
-
期刊类型引用(2)
1. 周斌，张弘，林春秀，付冬平. 降雨条件下基坑坡度和平台宽度对稳定性的影响. 广东土木与建筑. 2025(02): 28-31 . 
百度学术
2. 樊鸿志，王光进，蓝蓉，刘明生. 露天矿山岩质边坡软弱夹层赋存状态影响边坡稳定性规律研究. 重庆大学学报. 2024(04): 22-33 . 百度学术
其他类型引用(2)
-
其他相关附件
计量
- 文章访问数: 331
- HTML全文浏览量: 36
- PDF下载量: 82
- 被引次数: 4
