Experimental research on visible seepage of sandstone fissure using digital image-based method
-
摘要: 自主研发的可视化裂隙渗流试验装置,应用颜色示踪的图像数字化处理技术,在不同渗透压和法向应力分级加卸载耦合作用下,进行空间垂直角度的砂岩裂隙可视化渗流试验。通过毫秒级分帧技术捕捉裂隙渗流状态并对渗流面积色域分区,基于数字化图像技术进行二值化处理,通过图层二次叠加渲染,自识别渗流面积参数,建立渗透压与法向应力作用下的渗流面积扩散规律。研究裂隙渗流优势路径各断面的过流宽度变化特征,提出可视化渗流速度求取公式,指出渗流速度沿流径变化为骤增、骤降、均匀损耗3个阶段,建立最小断面宽度对应的峰值流速与渗透压和法向应力耦合作用下的函数模型。基于实测数据建立主渗流路径中法向应力、渗透压与雷诺数的幂函数关系,标定渗流状态的惯性作用区与黏性作用区的转变临界点,构建可视化裂隙渗流雷诺数预测模型,对裂隙渗流中真实渗流路径识别、流速矢量实时变化、流体状态判定等关键科学问题提出了新的研究理论及方法。Abstract: The visual fracture seepage test devices independently developed apply the color tracer image digital processing technology to carry out the visual seepage tests on the sandstone fractures with vertical spatial angle under the coupling effects of loading and unloading with different osmotic pressure and normal stress gradations. By means of the millisecond frame division technology, the seepage state of the crack is captured and the color domain of the seepage area is partitioned. Based on the digital image technology, the binary processing is carried out. The parameters of the seepage area are identified by double overlay of layers, and the diffusion law of the seepage area under the action of osmotic pressure and normal stress is established. The characteristics of the change of the flow width of each section along the dominant seepage path of the fracture are studied, and the formula for calculating the visual seepage velocity is put forward. It is pointed out that the seepage velocity changes along the flow diameter into three stages: sudden increase, sudden decrease and uniform loss. A functional model for the peak velocity corresponding to the minimum section width is established under the coupling action of osmotic pressure and normal stress. Based on the measured data of the main seepage path, the power function for the stress normal, osmotic pressure and Reynolds number is established, the change point of seepage state in the inertia function area and viscous effect area is calibrated, and the prediction model for Reynolds number of visualized fissure flows is formulated. The new theories and research method for the key scientific issues such as identification of real seepage in fissure seepage path, real-time change of velocity vector and determination of fluid state are put forward.
-
0. 引言
泥水平衡盾构广泛应用于水下隧道建设中，其采用拌制的泥水悬浮液（泥浆）作为支护材料，通过泥浆渗透形成不透水泥膜，将泥水压力转化为有效支护力从而平衡开挖面水土压力[1]。盾构机在进行带压开舱作业时，泥膜还需要有较好的闭气保压能力。而在高渗透地层，由于地层孔隙大，普通的膨润土泥浆大量流失，泥膜质量差，往往无法满足闭气要求甚至无法成膜。
目前普遍采用向泥浆中添加粗颗粒的方法提高成膜成功率和泥膜质量，常用的粗颗粒材料有天然河砂、锯末和蛭石等，已在Zimmerberg Base Tunnel [2]、南京和燕路隧道[3]等工程中成功应用。核桃壳是一种生物质粗颗粒材料，其主要成分为木质素、纤维素和半纤维素[4]，常应用于污水处理、砂浆改性等领域，在长期浸泡和酸碱盐浸渍条件下均有较好的耐久性[5]。由于核桃壳颗粒具有一定强度且吸水膨胀后可产生柔性变形，因此可同时承担骨架和填充变形作用[6]，因而也是钻井液最常用的堵漏材料之一。近年来，由于核桃壳在满足工程经济性的基础上，拥有比砂子更好的地层适应性和泥浆匹配性，含核桃壳泥浆在高渗透地层的盾构施工中也开始应用。
明确泥浆流变性、渗透成膜特性是施工中相关参数设定的基础。在泥浆流变性方面，盾构泥浆被认为是一种塑性流体或屈服假塑性流体[7-8]，泥浆中固相含量的增加对浆液的流变行为有显著影响，主要表现为黏度增加和流变曲线形状改变[9-10]。在泥浆渗透性方面，国内外众多学者通过室内泥浆渗透试验对泥浆渗透成膜进行全过程分析，把泥浆在地层中的渗透分为液相滤失和固相滞留两部分，提出了渗滤模型[11]，并对颗粒的堵漏机制[12]、堆积形态[13]、颗粒与地层孔隙的匹配关系[14]开展了进一步研究。在带压开舱时，粗颗粒对泥膜闭气性也有明显影响，粗颗粒含量的适当增加可改善地层内部网架结构的稳定性，从而提高泥膜的闭气性能[15]。但以上研究主要针对含砂泥浆，对含核桃壳泥浆各项性质认识不足，导致施工中泵送压力、泥浆配比、闭气压力等参数的设定多依靠经验。
本文拟采用核桃壳作为盾构泥浆粗颗粒材料，研究添加核桃壳后泥浆流变性变化，同时进行高渗透地层下的泥浆渗透和泥膜闭气试验，对比含核桃壳泥浆和含砂泥浆的渗透成膜特性和泥膜闭气性，以期为高渗透地层的带压开舱泥浆配制提供施工指导和建议。
1. 试验材料与方法
1.1 试验材料
（1）地层材料
为了模拟高渗透地层的孔隙特征，泥浆渗透及泥膜闭气试验所用地层采用粒径范围为0.5~2 mm的连续级配砂土，地层平均孔径为0.205 mm。控制每次试验的地层干密度为1.441 g/cm3，地层孔隙率为0.461，采用常水头渗透试验测得地层渗透系数为9.04×10-3 m/s。
（2）泥浆材料
试验中采用的泥浆材料主要包括纯净水、钠基膨润土、羧甲基纤维素钠（CMC）、核桃壳和天然河砂。试验用基浆按1︰8的膨水比配置，并添加质量分数为1‰的CMC作为增黏剂。其中，钠基膨润土采用0.075 mm标准筛做筛分处理，以消除杂质的影响，粗粒材料粒径的选取依据Min等[16]提出地层匹配性准则。试验采用的核桃壳由废弃物中回收而来，通过清洗、破碎、分拣、筛分等工序，制备成不同大小的颗粒。泥浆配制完成后，充分膨化24 h后对其进行基本性质测定。试验采用相同配比基浆（LB0/ST0），添加不同种类、含量的粗颗粒进行泥浆配制，研究泥浆的流变性、渗透性和泥膜闭气性，试验组别安排与泥浆基本性质参数如表 1所示。其中，编号ST0-ST4的含核桃壳泥浆同时进行流变性测试和泥浆渗透及泥膜闭气试验，编号ST5-ST8的含砂泥浆仅进行泥浆渗透及泥膜闭气试验。
表 1 泥浆基本性质参数Table 1. Basic parameters and properties of slurry泥浆编号 粗颗粒材料 添加量/(g·L-1) 密度/(g·cm-3) 漏斗黏度/s 胶体率/% ST0 — 0 1.075 35 83.5 ST1 0.15~0.2 mm核桃壳 10 1.082 55 87.0 ST2 20 1.086 61 90.0 ST3 30 1.090 65 92.0 ST4 40 1.092 68 88.0 ST5 0.15~0.2 mm砂 10 1.082 49 81.0 ST6 20 1.093 52 87.0 ST7 30 1.100 57 90.0 ST8 40 1.106 51 88.4 1.2 试验装置
（1）NXS-11型旋转黏度计
泥浆流变性测试试验采用NXS-11型旋转黏度计，采用控制剪切速率的方式进行流变测试。黏度计共5组测量系统，可测量剪切速率范围为15.5~996.1 s-1，共15级转速可调。
（2）泥浆渗透及泥膜闭气试验系统
采用自制的泥浆渗透及泥膜闭气试验系统，由加压系统、渗透柱和渗透流量采集系统3部分组成，加压系统包括空压机、调压装置和压力表；渗透柱为内径10 cm高70 cm的有机玻璃柱；渗透流量采集系统包括高精度天平、滤液收集器及配套的数据采集软件，如图 1所示。
![]()
图 1 泥浆渗透及泥膜闭气试验系统Figure 1. Apparatus for measuring volume of slurry permeation and airtightness of filter cakes1.3 试验步骤
（1）泥浆流变性测试试验
试验前，根据试验泥浆的预估黏度选择合适的测量系统，本试验选用测量系统A（外筒内径4 cm，内筒外径3.846 cm）。安装仪器，外筒中倒入规定量待测浆液并按照顺序安装内、外筒。待仪器各部分安装完毕后将转子置于工作状态，转动旋钮调整转速，待指针摆动稳定后（±0.5格）读取对应的刻度盘读数，按下式计算得到某一剪切速率所对应的泥浆剪切应力及表观黏度：
τ=N⋅z，η=γτ。} (1) 式中：τ为剪切应力（Pa）；N为刻度盘读数；z为转子常数，取0.2767 Pa/格；η为表观黏度（Pa·s）；γ为剪切速（s-1）。
试验中采用两种不同的剪切速率加载方式。其中，加载曲线1用于测试泥浆在剪切速率分级增加和降低条件下的流动性变化，分别可得到流变上行曲线（即流变曲线）和下行曲线；加载曲线2用于测试泥浆在恒定剪切速率下长时间剪切的流动性变化，剪切速率加载曲线如图 2所示。
（2）泥浆渗透及泥膜闭气试验
首先进行地层填装，底部滤层采用6~8 mm白色砾石，高度为5 cm，试验地层采用粒径为1~2 mm的天然河砂，高度26 cm，自下向上进行地层饱和，完成后缓慢注入制备好的泥浆1000 ml，静止5 min后开始泥浆渗透试验。
试验采用分级加压，共50，100，150，200 kPa共4级压力，每级压力下保持3 min。待每级压力稳定后打开渗透柱底部的出水阀门，数据采集系统自动采集流量，采集间隔为1 s。试验结束后缓慢对渗透柱进行泄压，打开渗透柱侧壁的阀门将剩余泥浆排出，观察泥膜形态并测量泥膜厚度。
对泥浆渗透试验中形成稳定泥膜的泥浆组进行泥膜闭气试验，采用分级加压的方式测定泥膜的耐高气压能力，以30 kPa为增量逐级增压至泥膜破坏，每级压力下保持1 min。试验结束后取出泥膜，观察泥膜的破坏形态并进行电镜扫描。
2. 试验结果分析
2.1 核桃壳颗粒对盾构泥浆流变性影响分析
（1）泥浆流变曲线分析
对ST0~ST4共5组含核桃壳颗粒泥浆进行流变性测试，得到宽剪切范围下泥浆的剪切应力-剪切速率（τ-γ）曲线（上行曲线）和表观黏度-剪切速率（η-γ）曲线，其中表观黏度为剪切应力与其相对应的剪切速率之比。泥浆的τ-γ曲线如图 3所示，当核桃壳含量为0 g/L时，泥浆表现为典型的屈服假塑性流体。加入核桃壳后τ-γ曲线形态产生改变，低剪切速率时剪切应力增长较快，中剪切速率处出现剪切应力平台。
泥浆的η-γ曲线如图 4所示，在中低剪切速率段表观黏度迅速降低，并随剪切速率的增大逐渐趋于一定值，即表现为明显的剪切稀释性。加入核桃壳前后曲线形态并无明显变化，但核桃壳的添加使泥浆在低剪切速率下的初始黏度有明显提升，黏度降幅增大，剪切稀释特性更为显著。
由于τ-γ曲线中段剪切应力平台的出现，现有的流变模型均不能完整反映添加核桃壳后泥浆的流变行为。为准确表达含核桃壳泥浆的流变特点，基于泥浆在不同剪切速率下表现出的不同流动特性，将其分为中低剪切速率下的屈服假塑性区和高剪切速率下的塑性区，分别采用不同的流变模型进行拟合，以保持较高的拟合相关程度。
由于泥浆中胶体颗粒间的相互作用，泥浆内部形成较为稳定的空间网架结构。在低剪切速率区，泥浆的空间网架结构逐渐破坏，泥浆表观黏度迅速降低。达到某一剪切速率后，泥浆的网架结构基本被破坏，此时泥浆黏度的主要由流体内部摩擦决定，这种摩擦贡献逐渐趋于定值，泥浆达到稳定流动状态，表观黏度不再变化。基浆中以胶体颗粒为主，固体粗颗粒较少，泥浆可以在较低的剪切速率下达到稳定流动。加入核桃壳后，粗颗粒间的碰撞和无序运动延长了泥浆的不稳定流动段。随着剪切速率的继续增大，在高剪切速率区，粗颗粒运动产生的摩擦贡献逐渐趋于稳定，泥浆也达到稳定流动状态，黏度不再随剪切速率变化。
在中低剪切速率区，泥浆的流动符合屈服假塑性流体的特征，采用幂律模型对该段进行拟合，拟合相关度达到了0.83以上。其流变控制方程为τ=Kγn，其中K为稠度系数，反映泥浆的流动阻力，若K值过大，则泥浆泵送困难，过小则泥浆携渣能力不足；n为流动指数，反映流体的流动行为偏离牛顿行为的程度（当n＜1时，流体表现为剪切稀释性），在施工过程中，一般希望泥浆有较低的n值，以便于泥浆流动与土渣携带，拟合结果参数见表 2。核桃壳颗粒的加入在降低泥浆流动指数的同时提升了泥浆的稠度，对工程中泥浆的泵送性和携渣能力均有正向贡献。
表 2 中低剪切速率区幂律模型拟合结果Table 2. Fitting results by power law model at low and medium shear rates核桃壳含量/(g·L-1) 幂律模型τ=Kγn n K/(Pa·s) R2 0 0.270 0.995 0.993 10 0.166 2.454 0.841 20 0.142 2.736 0.843 30 0.168 2.915 0.832 40 0.157 2.163 0.883 采用Bingham模型对泥浆稳定流动后高剪切速率下的塑性区进行拟合，拟合相关度达到了0.98以上。其流变控制方程为τ=τ0+μγ，其中τ0为泥浆的屈服应力，反映了泥浆内部各组分之间相互作用力的大小，亦即形成三维空间网架结构的能力；μ为泥浆的塑性黏度，该黏度不随剪切速率变化，与幂律模型中的稠度系数K具有相似的物理意义，反映泥浆稳定流动时的流动阻力，拟合参数结果如表 3所示。
表 3 高剪切速率区Bingham模型拟合结果Table 3. Fitting results by Bingham at high shear rates核桃壳含量/(g·L-1) Bingham模型τ=τ0+μγ τ0/Pa μ/(Pa·s) R2 0 2.024 0.0038 0.995 10 3.102 0.0059 0.998 20 3.444 0.0064 0.989 30 3.601 0.0067 0.982 40 3.360 0.0072 0.999 泥浆屈服应力τ0和塑性黏度μ随核桃壳含量的变化曲线如图 5所示，核桃壳颗粒的添加显著提高了泥浆的屈服应力，这表示核桃壳颗粒较好地嵌入了泥浆的三维网架结构，并提高了泥浆三维空间网架结构的稳定性。另外，随着核桃壳添加量的增加，总表面积增大，泥浆各组分间内摩擦力随之增加，宏观表现为泥浆塑性黏度的增大。
（2）泥浆的触变性分析
泥浆触变性是泥浆经剪切后变稀而静置后变稠的特性。膨润土泥浆作为典型的触变性非牛顿流体，主要特征为在剪切速率恒定的条件下，泥浆的剪切应力逐渐降低并趋于稳定，该值的大小和稳定时间与泥浆自身结构和剪切速率有关。该性能对于描述长时间输送后泥浆流变特性的变化具有重要意义。
经过计算，以盾构施工常用的泥浆输送系统为例，管径0.35 m、流量1000 m3/h的泥浆泵送设备的有效剪切速率约为124.75 s-1，因此采用旋转黏度计的第7级转速（105.1 s-1），在该转速下测试泥浆长时间剪切的流变性衰减，如图 6所示。虽然各组泥浆的初始黏度不同，但经过恒定剪切速率的长时间作用后，剪切应力趋于一个较为统一的稳定值，达到这一稳定值的时间也无较大区别。因此，核桃壳颗粒的加入对盾构泥浆固定速率下长时间剪切后的流变性并无明显改变。
![]()
图 6 固定剪切速率下泥浆剪切应力衰变曲线Figure 6. Decay curves of shear stress of slurry at constant shear rate触变性流体在流变曲线上的表现为剪切速率上行和下行曲线不重合，形成一个滞后圈，可采用滞后圈面积表示流体触变性[17]。触变滞后圈包含的面积被定义为使材料三维网格或凝胶结构被破坏所需要的能量，其量纲为能量/体积。
核桃壳添加量为0，30 g/L的泥浆滞后圈如图 7所示，添加核桃壳前后泥浆的滞后圈形态有较大差异，主要体现在上行曲线形态的改变和下行曲线斜率的不同。
通过积分计算各组泥浆滞后圈面积如图 8所示，加入核桃壳后泥浆滞后圈面积显著增加，随着核桃壳添加量的增加，泥浆滞后圈面积先增加后减小，在添加量为30 g/L时达到峰值。这表示核桃壳的加入提升了泥浆内部网架结构和凝胶结构强度，破坏该结构需要更多的能量，核桃壳添加量为30 g/L的泥浆具有的凝胶强度最高，触变性最好。
2.2 含核桃壳泥浆渗透特性分析
对添加核桃壳颗粒（ST1-4）和天然河砂（ST5-8）的九组泥浆开展高渗透砂层下的泥浆渗透试验，试验中的总渗透流量变化曲线分别如图 9，10所示。九组泥浆中，除了未添加粗颗粒的基浆外，均形成了稳定的泥膜+渗透带型低渗透区。由于高渗透砂层的地层孔径过大，基浆的液相和固相颗粒均穿过地层孔隙造成滤失，因此无法形成稳定的低渗透区；加入合适粒径粗颗粒后，一部分粗颗粒在地层内部滞留堵塞降低了地层渗透系数形成渗透带，另一部分粗颗粒在地层表面不断堆积形成泥膜，该过程中泥浆渗透流量逐渐降低并趋于稳定，这标志着不透水泥膜的形成。
![]()
图 9 含核桃壳泥浆渗透流量曲线Figure 9. Curves of discharged water of slurry containing walnut shell图 11为泥浆渗透成膜的总渗透流量变化图，随着粗颗粒含量的增加，泥浆的成膜滤失量呈现出先减小后增加的趋势，且含核桃壳泥浆的成膜渗透流量大于传统含砂泥浆。另外，试验中发现，含核桃壳泥浆生成的泥膜厚度（约为8~10 mm）大于含砂泥浆（约为5~6 mm）。由于试验中的渗透方向为竖向渗透，砂颗粒比重较大，可以较快地沉积于地层表面，而核桃壳颗粒与泥浆比重相近，均匀分散悬浮于泥浆中，随泥浆缓慢渗透进入地层堵塞地层孔隙，因此在竖向渗透试验中含核桃壳泥浆的成膜时间较长，渗透流量较大。但在实际盾构施工中，泥浆为横向渗透，将不会出现此种渗透流量的差别。反而由于砂颗粒密度大，多聚集于开挖舱底部，上部由于缺乏堵漏材料而产生大量漏失，常出现泥膜质量不均匀甚至无法成膜的问题。
![]()
图 11 泥浆总渗透流量变化曲线Figure 11. Relationship between water discharge and coarse content of slurry2.3 含核桃壳泥浆与含砂泥浆的泥膜耐高气压能力对比
在带压开舱时，泥水舱的空气压力经泥膜转化为有效应力用以平衡地层的水土压力，从而维持开挖面稳定。泥膜在一定气压下不被气体击穿并保持这种压力转化功能的能力即为泥膜的闭气性能。
对ST1-4、ST6-8共7组含粗颗粒泥浆渗透形成的泥膜进行耐高气压能力测试试验，试验过程中的渗透流量变化曲线分别如图 12，13所示。
![]()
图 12 含核桃壳泥膜闭气渗透流量曲线Figure 12. Curves of discharged water of slurry containing walnut shell渗透流量曲线在不同位置出现突变点，即泥膜的进气点。由于泥膜内部实际由许多孔隙连通形成毛管，研究中通常将此现象简化为多孔介质毛管模型进行解释：进气点之前，外气压与泥膜孔隙毛管中液体的表面张力形成平衡，此时渗透流量几乎为零，气舱压力无明显变化；随着气压的增大，毛管液体弯液面不断“下凹”，以获得更大的垂向张力来抵抗气压。当弯液面变为半球形时达到极限状态（即进气点），外气压瞬间完成对泥膜孔隙液体的驱替，渗透流量瞬间增大，泥膜表现为圆孔状击穿破坏，如图 14所示。
采用闭气值作为泥膜耐高气压能力的评价指标，闭气值是泥膜在压气条件下能维持闭气性能的最大压力值，闭气值越大，泥膜质量越好，在试验中取为泥膜进气点的前一级气压。泥膜闭气值随粗颗粒种类和含量的变化如表 4所示。对于同种粗颗粒材料，随着添加量的增加，泥膜闭气值逐渐增大，在30 g/L分别达到最大值120，90 kPa。随着含量的继续增加，过多的粗颗粒表面吸附了泥浆中大量的水，抑制膨润土颗粒水化的同时颗粒平均水化膜厚度变薄，泥膜各组分间的连接紧密性降低，宏观表现为泥膜耐高气压能力的降低。
表 4 泥膜闭气值汇总表Table 4. Threshold pressures of airtightness of filter cakes添加量/(g·L-1) 泥膜闭气值/kPa 核桃壳 砂 10 30 — 20 90 30 30 150 90 40 60 30 对于不同的粗颗粒材料，在相同质量含量条件下，含核桃壳泥浆的泥膜闭气值均高于含砂泥浆。两种泥浆均在粗颗粒添加量为30 g/L达到泥膜闭气的最优值，且含核桃壳泥浆的泥膜最优闭气值较含砂泥浆高30 kPa。因此，含核桃壳泥浆形成的泥膜具有更好的耐高气压能力。
3. 讨论
3.1 核桃壳对泥浆性质影响机理分析
盾构泥浆是一种多相分散体系，主要由膨润土或黏土等胶粒和分散介质（水）组成，胶体颗粒与分散介质的界面性质对泥浆的性质起关键作用[18]。表面负电性的膨润土颗粒和极性水分子通过物理吸附作用建立联系，形成以膨润土胶体颗粒为“节点”，以水化阳离子为“桥”的三维网架结构，在静止状态下泥浆保持均匀分散和稳定性。核桃壳的主要成分为纤维素、半纤维素、木质素等，且基本单元中含有糖类物质，富含极性羟基，因而具有较强的亲水性[18-19]，可以通过化学连接和物理吸附作用吸引水化阳离子形成水化膜。因此，大量的膨润土胶体颗粒可以于核桃壳表面建立有效连接并将其包裹，使其较好地嵌入到泥浆的三维网架结构中，提高了结构的稳定性，宏观表现为泥浆胶体率和屈服应力的提高。通过电镜扫描试验对含核桃壳泥浆形成的泥膜进行细观分析，如图 15所示。颗粒表面与连接处均被叠状连接的膨润土胶粒包裹，形成致密的结构，颗粒界线轮廓已不明显。
![]()
图 15 含30 g/L核桃壳泥膜细观分析图Figure 15. Microscopic analysis of slurry containing walnut shell of 30 g/L3.2 核桃壳改善泥膜特性作用机理分析
泥浆在向地层中渗透的过程中，大部分液相滤失、固相则堆积成为泥膜。泥膜的基本单元体为骨架状结构，骨架间的孔隙是泥膜进气的“突破点”。含粗颗粒泥浆形成的泥膜主要以粗颗粒作为骨架结构，骨架孔隙中由膨润土颗粒和CMC大分子进行填充。由于粗颗粒刚度较大且稳定性强，泥膜高气压破坏的击穿圆孔多位于粗颗粒之间的孔隙处，因此，填充泥浆的凝胶强度和粗颗粒与填充物间的界面结合强度是泥膜耐高气压能力的关键影响因素。
传统泥浆常用的堵漏材料，如天然河砂，为刚性材料，颗粒之间为硬接触，如图 16（a）所示。而核桃壳颗粒吸水后表面会发生软化，可产生一定柔性变形和挤压，颗粒间为软接触，如图 16（b），因此颗粒间的孔隙明显小于刚性颗粒，孔隙中可平衡外气压的表面张力越大，泥膜越不容易被气体击穿，闭气值越大。
![]()
图 16 颗粒接触方式及孔隙大小对比Figure 16. Contact modes of particles and correlation of pore sizes除此之外，结合2.1节中对含核桃壳泥浆流变性的测试结果，核桃壳粗颗粒的加入可提高泥浆黏度和在高剪切速率下的凝胶强度，触变性也有明显改善，这也使泥膜中颗粒间的黏结性更强、泥膜的三维空间网架结构更为稳定，不易被高气压冲破。
4. 结论
针对传统含砂泥浆在高渗透砂层带压开舱时遇到成膜困难及闭气性不足的问题，本文选用轻质有机颗粒核桃壳作为泥浆堵漏材料，进行流变性测试试验、泥浆渗透成膜及泥膜闭气试验，对含核桃壳泥浆的流变性和渗透成膜特性进行研究并分析其作用机理，主要得到4点结论。
（1）核桃壳可以增强泥浆的分散稳定性。核桃壳颗粒由于表面水化和吸附作用可以较好地嵌入泥浆的三维网架结构，提高胶体结构稳定性，泥浆黏度和胶体率均有所提升。
（2）核桃壳颗粒的添加会影响泥浆的流变性，在提升泥浆屈服应力和黏度的同时保持了较好的流动性。含核桃壳泥浆的流变曲线形态发生改变，可以分为中低剪切速率下的屈服假塑性区和高剪切速率下的塑性区，分别采用幂律模型和Bingham模型拟合的相关性较高。
（3）相同添加量下，含核桃壳泥浆较传统含砂泥浆的泥膜耐高气压能力更强。与砂子间的硬接触相比，核桃壳颗粒的软接触方式可以减少颗粒间孔隙尺寸、且泥浆流变性和触变性均有一定程度改善，使泥膜更加致密，耐高气压能力更强。
（4）盾构施工中，在泥浆性质方面，应尽量选用稳定性高、泵送性强、凝胶强度高的泥浆，在泥膜质量方面，应尽量选择闭气值高的泥膜。综合考虑各试验结果，采用核桃壳替代传统砂子可改善泥浆性质和泥膜质量，在渗透系数为9.04×10-3 m/s的高渗透性砂层中，推荐的添加量为30 g/L。
-
![]()
图 7 分帧、分区、特征点提取图像示例
Figure 7. Example of image extraction by frame, partition and feature points
![]()
图 8 渗流面积临界Ⅰ值求取示意图
Figure 8. Schematic diagram of critical Ⅰ value of contact area of fracture surface
![]()
图 10 识别盲区与误区
Figure 10. Identification blind and misunderstanding area of seepage area
![]()
图 15 渗流面积与法向应力、渗透压的关系
Figure 15. Relationship among seepage area and normal stress and osmotic pressure
![]()
图 19 法向应力0.1 MPa时各断面速度关系
Figure 19. Relationship between speed section at normal stress of 0.1 MPa
![]()
图 20 法向应力0.6 MPa时各断面速度关系图
Figure 20. Relationship between speed sections at normal stress of 0.6 MPa
![]()
图 22 渗透压与峰值流速vmax的影响关系
Figure 22. Relationship between osmotic pressure and maximum velocity
![]()
图 24 0.1 MPa不同渗透压下渗流路径长度与雷诺数的关系
Figure 24. Reynolds number of flow path under different osmotic pressures at normal stress of 0.1 MPa
表 1 渗流面积分类与RGB值对应表
Table 1 Classification of seepage area and values of RGB
区域命名 渗流面积分类 颜色表观 RGB数值大小反映颜色的偏向 临界 归一化 R G B Ⅰ值 Ⅰ值 1区域 有效渗流面积 液体扩散区 70~105 40~85 35~85 0.32 0.11~0.32 2区域 临界渗流面积 液体微过渡区 101~130 51~90 51~90 — 0.32~0.342 3区域 非渗流面积 液体未扩散区 131~160 90~120 90~120 0.342 0.342~0.918 
下载: 导出CSV
表 2 黑白高亮二值化图
Table 2 Black and white highlighted binarization diagram
法向应力/MPa 渗透压/MPa 0.1 0.2 0.3 0.4 0.5 0.6 0.1 
0.6 
下载: 导出CSV
表 3 裂隙固有渗流参数δ值
Table 3 Values of intrinsic seepage parameter of crack
法向应力/MPa 0.1 0.2 0.3 0.4 0.5 0.6 δ值 49.73 47.17 45.93 41.90 39.17 23.43
下载: 导出CSV
表 4 主惯性与主黏性临界点
Table 4 Critical points of main inertia and main viscosity
临界点坐标 渗透压/MPa 0.1 0.2 0.3 0.4 0.5 0.6 x（X） 8.33 7.52 7.65 7.93 8.08 8.52 y（Re） 0.71 1.15 1.18 1.74 1.86 2.27
下载: 导出CSV
表 5 雷诺数预测值的平均相对误差
Table 5 Average relative errors of predicted values of Reynolds number
法向应力/MPa 0.1 0.2 0.3 0.4 0.5 0.6 平均相对误差ε 0.2659 0.1442 0.1811 0.1872 0.2059 0.1667
下载: 导出CSV
-
[1] LOUIS C. Rock Hydraulics in Rock Mechanics[M]. New York: Springer-Verlag, 1974.
[2] 速宝玉, 詹美礼, 赵坚. 仿天然岩体裂隙渗流的试验研究[J]. 岩土工程学报, 1995, 17(5): 19-24. doi: 10.3321/j.issn:1000-4548.1995.05.004 SU Bao-yu, ZHAN Mei-li, ZHAO Jian. Study on fracture seepage in the imitative nature rock[J]. Chinese Journal of Geotechnical Engineering, 1995, 17(5): 19-24. (in Chinese) doi: 10.3321/j.issn:1000-4548.1995.05.004
[3] 许光祥, 张永兴, 哈秋舲. 粗糙裂隙渗流的超立方和次立方定律及其试验研究[J]. 水利学报, 2003, 34(3): 74-79. doi: 10.3321/j.issn:0559-9350.2003.03.014 XU Guang-xiang, ZHANG Yong-xing, HA Qiu-ling. Super- cubic and sub-cubic law of rough fracture seepage and its experimental study[J]. Journal of Hydraulic Engineering, 2003, 34(3): 74-79. (in Chinese) doi: 10.3321/j.issn:0559-9350.2003.03.014
[4] SINGH K K, SINGH D N, RANJITH P G. Laboratory simulation of flow through single fractured granite[J]. Rock Mechanics and Rock Engineering, 2015, 48(3): 987-1000. doi: 10.1007/s00603-014-0630-9
[5] BABADAGLI T, REN X, DEVELI K. Effects of fractal surface roughness and lithology on single and multiphase flow in a single fracture:an experimental investigation[J]. International Journal of Multiphase Flow, 2015, 68: 40-58. doi: 10.1016/j.ijmultiphaseflow.2014.10.004
[6] BRUSH D J, THOMSON N R. Fluid flow in synthetic roughwalled fractures:Navier-Stokes, Stokes, and local cubic law simulations[J]. Water Resources Research, 2003, 39(4): 1085-1100.
[7] OLSSON R, BARTON N. An improved model for hydro mechanical coupling during shearing of rock joints[J]. International Journal of Rock Mechanics and Mining Sciences, 2001, 38(3): 317-329. doi: 10.1016/S1365-1609(00)00079-4
[8] 胡昱, 源新, 刘光廷, 等. 多轴应力作用下砂砾岩单裂隙渗流规律试验研究[J]. 地下空间与工程学报, 2007, 3(6): 1009-1013. doi: 10.3969/j.issn.1673-0836.2007.06.006 HU Yu, YUAN Xin, LIU Guang-ting, et al. Experiment research on the laws of seepage in calcirudite rock within single fracture under multiaxial stresses[J]. Chinese Journal of Underground Space and Engineering, 2007, 3(6): 1009-1013. (in Chinese) doi: 10.3969/j.issn.1673-0836.2007.06.006
[9] 刘才华, 陈从新. 三轴应力作用下岩石单裂隙的渗流特性[J]. 自然科学进展, 2007, 17(7): 989-994. doi: 10.3321/j.issn:1002-008X.2007.07.022 LIU Cai-hua, CHEN Cong-xin. Seepage characteristics of single fracture under triaxial stress[J]. Advances in Natural Science, 2007, 17(7): 989-994. (in Chinese) doi: 10.3321/j.issn:1002-008X.2007.07.022
[10] 王志良, 申林方, 徐则民, 等. 岩体裂隙面粗糙度对其渗流特性的影响研究[J]. 岩土工程学报, 2016, 38(7): 1262-1268. https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201607013.htm WANG Zhi-liang, SHEN Lin-fang, XU Ze-min, et al. Influence of roughness of rock fracture on seepage characteristics[J]. Chinese Journal of Geotechnical Engineering, 2016, 38(7): 1262-1268. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YTGC201607013.htm
[11] QIAN J, ZHAN H, LUO S, et al. Experimental evidence of scale-dependent hydraulic conductivity for fully developed turbulent flow in a single fracture[J]. Journal of Hydrology, 2007, 339(3/4): 206-215.
[12] RANJITH P G, DARLINGTON W. Nonlinear single-phase flow in real rock joints[J]. Water Resources Research, 2007, 43(9): 146-156.
[13] ZHANG Z, NEMCIK J. Fluid flow regimes and nonlinear flow characteristics in deformable rock fractures[J]. Journal of Hydrology, 2013, 477(16): 139-151.
[14] 胡少华, 周佳庆, 陈益峰, 等. 岩石粗糙裂隙非线性渗流特性试验研究[J]. 地下空间与工程学报, 2017, 13(1): 48-56. https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201701008.htm HU Shao-hua, ZHOU Jia-qing, CHEN Yi-feng, et al. Laboratory research on nonlinear flow behavior of rough fractures[J]. Chinese Journal of Underground Space and Engineering, 2017, 13(1): 48-56. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-BASE201701008.htm
[15] 鞠杨, 谢和平, 郑泽民, 等. 基于3D打印技术的岩体复杂结构与应力场的可视化方法[J]. 科学通报, 2014, 59(32): 3109-3119. https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201432002.htm JU Yang, XIE He-ping, ZHENG Ze-ming, et al. Visualization of the complex structure and stress field inside rock by means of 3D printing technology[J]. Chinese Science Bull, 2014, 59(32): 3109-3119. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-KXTB201432002.htm
[16] 刘建军, 汪尧, 宋睿, 等. 基于透明岩土材料的可视化渗流实验及其应用前景[J]. 地球科学, 2017, 42(8): 1287-1295. https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX201708004.htm LIU Jian-jun, WANG Yao, SONG Rui, et al. Visual seepage experiment based on transparent rock-soil material and its application prospect[J]. Earth Science, 2017, 42(8): 1287-1295. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-DQKX201708004.htm
[17] 盛金昌, 刘继山, 赵坚. 基于图像数字化技术的裂隙岩体非稳态渗流分析[J]. 岩石力学与工程学报, 2006, 25(7): 1402-1407. https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200607015.htm SHENG Jin-chang, LIU Ji-shan, ZHAO Jian. Analysis of transient fluid flow in fractured rock masses with digital image-based method[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(7): 1402-1407. (in Chinese) https://www.cnki.com.cn/Article/CJFDTOTAL-YSLX200607015.htm
[18] FOXRW , MCDONALD A T, PRITCHARD P J. Introduction to Fluid Mechanics[M]. 6th ed. Hoboken: John Wiley and Sons, 2004: 348.
计量
- 文章访问数: 316
- HTML全文浏览量: 42
- PDF下载量: 264
