不同加气和粒径条件下浑水滴灌滴头堵塞特性研究

　　摘要：为探明浑水滴灌过程中，水中加气对滴头堵塞的影响，以齿形迷宫流道滴头为研究对象，采用周期性间歇灌水测试方法，用最大粒径小于 0.1mm 的泥沙配置了 5 种不同的浑水，运用激光粒度分析仪和场发射扫描电镜等方法研究了滴头堵塞状况。结果表明：加气和泥沙颗粒级配对滴头堵塞具有极显著影响(P<0.01);加气提升了毛管内泥沙输移能力、促进大泥沙颗粒排出，减小淤积泥沙中值粒径，且泥沙最大粒径越小，加气对泥沙运移的影响越小，毛管淤积物质量、淤积泥沙中值粒径较未加气处理减少 8.75%～31.92%、8.59%～35.64%;泥沙粒径为 0.075～0.1mm 时，滴头流量下降最快，下降了 70.3%，泥沙为粒径 0～0.1mm 时，滴头流量下降最慢，下降了 40.0%;加气增大了水流紊动程度，促进浑水中大颗粒泥沙在流道内的输移;加气加剧了浑水中小颗粒泥沙在流道入口处黏附，加速了流道入口堵塞，滴头 Dra 和 Cu 比未加气处理低 9%～18.7%和 16.2%～36.4%，这是造成加气加速滴头堵塞的主要原因。建议进行毛管冲洗，降低流道入口堵塞风险，以提高加气滴灌滴头的抗堵塞性。

　　牛文全; 张二信; 吕畅; 孙军; 董爱红; 邬梦龙， 农业机械学报 发表时间：2021-10-15

　　关键词：迷宫流道;加气;泥沙粒径;堵塞

　　0 引 言

　　滴灌过程中，滴头附近作物根区土壤处于低氧状态，易造成植物根系缺氧，严重时会影响作物根系的正常生理活动。加气滴灌具有协同调节作物根区水、肥、气状况，改善根系生长环境，缓解作物根部缺氧问题，提高作物水肥利用效率、果实产量和品质等优点，被广泛应用于蔬菜、大田作物和果树种植上[1-3]。地下滴灌是加气灌溉的最佳方法，但滴头狭窄的流道易被水中的颗粒、化学沉淀、溶解盐、微生物和其他杂质堵塞 [4]，众多学者试图通过过滤设备配置、滴头流道优化及毛管冲洗等方法解决堵塞问题[5-7]，但该难题至今仍没有得到很好的解决。

　　灌溉水源的复杂性决定了滴头堵塞类型的多样性，一般有物理、化学和生物堵塞 3 种类型。灌溉水中细小泥沙颗粒(粒径小于 0.1mm)易在流道内发生聚积、沉降、黏附等行为造成滴头发生物理堵塞 [8-9]，文献[10]发现粒径范围为 0.031～0.038mm 的泥沙在流道内易发生团聚行为形成团聚体，黏附于流道壁面不易随水流输出，文献[11]研究表明，毛管中粒径范围为 0.10～0.15mm 的泥沙主要以推移质的形式进入灌水器，且易引发滴头发生突然完全堵塞。此外，水力特性的改变会影响悬浮颗粒运移过程中所受的黏滞阻力、摩擦力等，从而对滴头抗堵塞性产生影响。文献[12]研究发现，加气能够改变水流水力特性，增大水流流速，文献[13]发现，微气泡破碎时产生巨大的能量能够加速周围颗粒的运移，文献[14]研究表明，加气能够减小壁面的摩擦阻力，降低水流能量损失。目前，对于细小颗粒造成的滴头堵塞问题，多采用多级过滤器组合模式来减缓滴头堵塞，但这些措施极大增加了滴灌成本和能源损耗[15]。倘若明确细颗粒泥沙在加气条件下的输移、淤积规律，进而采取合理的管控措施来缓解滴头堵塞，减轻过滤系统负荷，对提高滴灌系统运行效率具有重要意义。

　　目前，针对加气条件下细小泥沙对滴头堵塞的影响研究较少，且加气对滴头堵塞机理的影响尚不明确。为此，本文通过研究加气对不同粒径浑水滴灌滴头堵塞的影响，对比分析加气前后滴头堵塞物结构及毛管淤积泥沙粒径的变化，以探明加气滴灌堵塞机理，为加气滴灌系统提高滴头抗堵塞能力提供依据。

　　1材料与方法

　　1.1 试验材料

　　据调查[16]，宁夏段黄河引水渠含沙量为 0.61～1.1g/L，故试验配置浑水含沙量为 1g/L，试验泥沙采自渭河漫滩，由于目前滴灌设备中多采用 150 目的过滤器，即通过过滤器的最大颗粒粒径为 0.1mm[17]，所以本研究选取泥沙最大粒径为 0.1mm，同时根据泥沙颗粒在水中的跟随性[17]、泥沙颗粒对堵塞的危险性[18]、以及已有研究确定的滴头堵塞敏感粒径范围[10,19]，将采集的泥沙经自然风干后过筛，分为 5 种粒径范围，配置 5 种不同粒径级配的浑水，泥沙颗粒机械组成如表 1 所示。试验用水为杨凌自来水，pH 值为 7.83～8.03，悬浮物质量浓度为 106～152mg/L，电导率为 217.3～372.6µS/cm，总硬度为 71～81mg/L，细菌菌体浓度小于 1CFU/mL，符合灌水标准。

　　依据文献[20]的研究，内镶贴片式滴头更适用于引黄灌区灌溉，故本试验选用内镶贴片式滴头(图 1)进行浑水测试，滴灌带参数为：额定流量 2.5L/h，毛管外径 16mm，壁厚 0.2mm，流道深(Z)0.63mm，流道宽(W)0.60mm，流道长 18.71mm，齿间距(L) 1.40mm，齿间角 37°。

　　1.2 试验装置

　　试验装置主要由抗堵塞测试平台和加气装置组成，如图 2 所示。抗堵塞测试平台由蓄水箱(长度、宽度、高度均为 60cm)、水泵(功率 0.75kW，扬程 60m)、过滤器(120 目网式过滤器)、搅拌机(功率 0.75kW)、压力表(量程 0.25MPa，精度 0.001MPa)、控制阀门、支管以及待测试毛管等组成。毛管布置于测试架上(长度 6m、宽度 0.5m、高度 1m)，共 3 个测试架，每个架子上铺设 3 根 6m 长滴灌带，每根滴灌带有 22 个滴头，滴头间距为 25cm。加气装置为 FRGW-10 型水肥气耦合机，(河南丰润环保科技有限公司生产)，经前期测定，最优加气时长为 5min，加气压力 0.45～0.46MPa，加气后形成乳白色水气混合液(水中微气泡粒径为 15µm，溶氧量为 8.36～8.45mg/L)。

　　1.3 试验方案和方法

　　试验分为清水试验和浑水试验两部分。

　　1.3.1 清水试验

　　清水试验包括加气和未加气清水测试，测试压力为 0.02～0.12MPa，压力间隔为 0.02MPa。测试前对测试平台进行冲洗，以保证试验数据的可靠性。测试方法：打开测试平台阀门，将压力表的示数调节至测定压力值，待示数稳定后，将烧杯(1000mL)置于毛管下方测定每个滴头的流量，测定时长为 10min，每组重复 3 次，取均值，清水流量计算公式为：

　　1.3.2 浑水试验

　　选用了 1 种齿形迷宫流道滴头的滴灌带，分别在加气和不加气条件下进行短周期灌水试验，测试不同粒径浑水滴灌滴头堵塞状况，浑水泥沙粒径级配如表 1 所示，浑水试验共计 10 组。

　　1.3.3 试验过程

　　试验在西北农林科技大学水力学大厅内进行，试验开展时间为 2021 年 4—6 月，试验测试时，先按照 1.3.1 节进行清水试验，后进行短周期间歇性浑水试验。灌水压力控制在0.1MPa，单次灌水持续时间为 1.5h，相邻灌水间歇时间为 0.5h。每次灌水结束前 10min 压力表的示数稳定在 0.1MPa 后测定毛管各滴头流量，重复 3 次取平均值，当累计灌水 44 次后结束测试。取下毛管置于通风处晾干，剖开毛管收集毛管淤积物和滴头堵塞物，并用 S4800 型场发射扫描电镜和 APA2000 型激光粒度仪对滴头堵塞物结构、毛管淤积物粒径分别进行观测。为保证浑水的均匀性，灌水期间搅拌机持续工作。每个处理灌水结束后，更换新的毛管并对测试平台进行冲洗。

　　1.4 评价指标和方法

　　本文采用平均相对流量 Dra 和克里斯森均匀度系数 Cu 来评价滴头整体的堵塞程度[21]，一般以实测滴头浑水流量占清水流量的百分比小于 75%，滴灌系统灌水均匀系数小于 80% 时判定滴头发生堵塞[22-23]。堵塞率计算公式为： = 100% c n n   (2)式中 n-滴头总数，个 nc-发生堵塞的滴头个数，个

　　1.5 数据处理

　　采用 SPSS25.0 软件进行统计分析，采用主效应方差分析泥沙粒径和加气对滴头堵塞的影响，用 LSD(最小显著差异)法检测各处理组间相对流量及堵塞物的差异性，利用 Origin2021 和 CAD2014 软件进行绘图。

　　2 结果与分析

　　2.1 加气对滴头水力性能的影响

　　将不同压力下测定的滴头流量代入式(1)，通过曲线拟合得到水力特性曲线如图 3 所示，可看出加气增大了滴头的流量系数，降低了流态指数，加气处理的流量系数 kd、流态指数 x 为 0.28 和 0.47，未加气的为 0.24 和 0.50。流量系数 kd越小说明流量波动越小，而流态指数 x 越小说明滴头流量对进口压力越不敏感 [24]。因此，本试验中加气增大了水流波动，降低了流量与压力之间的敏感度。对比同一压力下不同处理的滴头流量，可发现加气增大了滴头流量。

　　2.2 加气对滴头流量及均匀度的影响

　　滴头 Dra 和 Cu 随灌水次数的变化过程如图 4 所示。不同粒径浑水加气处理的滴头 Dra 和 Cu 下降速率大于未加气处理，其中，D2 浑水 Dra 下降速率最快，易造成滴头堵塞，D1 浑水 Dra 下降速率最慢，不易造成滴头堵塞。试验结束时，加气处理的滴头 Dra 和 Cu 比未加气处理低 9%～18.7%和 16.2%～36.4%。

　　灌水初期，加气处理的滴头 Dra 和 Cu 随时间增长呈缓慢下降趋势，而未加气处理的滴头 Dra和 Cu在 95%附近波动。当灌水 15 次时，加气处理的 D3、D4 浑水最先引发滴头发生堵塞，流量比未加气处理分别低 17.4%和 22.0%，加气处理的滴头 Dra 和 Cu 下降速度大于未加气处理。灌水 30 次后，5 种浑水加气处理的滴头 Dra 均小于 75%，滴头发生严重堵塞，而未加气处理中只有 D2 和 D4 浑水的滴头发生了堵塞。5 种浑水加气处理的滴头 Dra 分别比未加气处理低 11.3%、27.0%、11.3%、11.9%和 31.4%，D5 浑水加气与未加气处理的差异性最大。灌水结束时，D1、D2、D3、D4 和 D5 浑水加气处理的滴头 Dra分别比未加气处理低 18.6%、18.4、9%、8.9%和 18.7%，Cu 分别低 36.3%、26.8%、16.2%、28.5%和 36.4%。

　　灌水结束时，滴头流量方差和显著性分析结果表明(表 2、3)，泥沙粒径和加气对滴头堵塞具有极显著性影响(P<0.01)，且交互作用也极显著，表明加气加剧了滴头堵塞，且与泥沙级配密切相关(表 2)。D1 浑水加气处理的滴头 Dra、Cu与 D3、D4、D5 浑水处理差异性显著(表 3)，整个灌水过程中该浑水处理的滴头 Dra、Cu 均保持在最高的水平，表明 D1 浑水相比其他粒径浑水造成滴头堵塞的程度较轻，且各滴头堵塞程度差异性较小。

　　2.3 滴头相对流量和灌水均匀度的动态变化

　　从滴头 Cu和 Dra拟合结果(图 5)可以看出，滴头 Dra和 Cu具有明显的协同性，两者均随着灌水次数的增加而同步减小，D1～D4 浑水加气处理的 Cu和 Dra拟合斜率为 0.312～ 0.752，明显小于未加气处理的 0.334～0.915，而 D5 浑水则相反。斜率越小表明该滴头越易发生突然堵塞，因此，加气加剧了 D1～D4 浑水滴灌滴头突然堵塞的发生，减缓了 D5 浑水滴灌滴头突然堵塞的发生。加气与未加气处理斜率的差异性随泥沙最大粒径的减小而增大(D2～D5 浑水加气处理的拟合斜率分别是未加气处理的 0.822、0.84、0.87、1.46 倍)。这是由于滴头流道尺寸较小，无论加气与否，易发生沉降的大颗粒泥沙均易引发滴头发生突然堵塞，而对于小颗粒而言，未加气时流道内沉积的泥沙易随水流冲出，滴头 Dra 和 Cu 下降速率缓慢，而加气时泥沙黏附加剧造成滴头 Dra和 Cu下降较快(图 4)，因此，泥沙最大粒径越小，加气处理的滴头 Dra和 Cu协同性较未加气处理差异性更大。

　　2.4 泥沙粒径对滴头堵塞率的影响

　　不同处理的滴头堵塞率统计结果如图 6 (图中不同大写字母表示不同粒径之间差异显著，不同小写字母表示加气与未加气处理之间差异显著(P<0.05)，下同)所示。加气处理的滴头堵塞率显著高于未加气处理(P<0.05)，D1～D5 浑水未加气处理的堵塞率为 23.9%～ 46.3%，加气处理的堵塞率为 37.3%～61.2%，比未加气处理增大了 13.4～14.9 个百分点。无论是否加气，5 种浑水灌溉滴头堵塞率从大到小均依次为 D2、D4、D5、D3、D1，表明滴头堵塞的程度主要由泥沙机械组成所决定。加气处理下 D1 浑水的堵塞率与 D2～D4 浑水的堵塞率差异性显著(P<0.05)，且其堵塞率最低，表明 D1 浑水加气灌溉条件下，滴头经过长时间的灌溉后仍具有较强的出流能力，这与图 4 中 D1 浑水加气处理的滴头 Dra 和 Cu 随时间变化表现出的规律一致。

　　2.5 加气对毛管泥沙淤积量的影响

　　不同处理的毛管淤积泥沙参数如表 4、5 所示。加气显著降低了毛管淤积物质量、减小了淤积泥沙中值粒径、增大了淤积物比表面积，且各处理间差异性显著(P<0.05)。加气处理的毛管淤积物质量随着泥沙最大粒径的减小而增大，未加气处理无明显规律，加气处理的毛管淤积物较未加气降低 8.75%～31.92%。泥沙中值粒径随着泥沙最大粒径的减小而减小，比表面积则相反，D1～D5 浑水加气处理的中值粒径为未加气处理的 0.64、0.88、 0.89、0.91、0.93 倍，较未加气减小 8.59%～35.64%，比表面积为未加气处理的 1.59、 2.33、1.23、1.19、1.01 倍，较未加气增大 6.01%～132.81%。泥沙中值粒径越大，表明沉积的大颗粒泥沙越多，比表面积越大表明沉积的细颗粒泥沙越多。因此，加气促进了毛管内泥沙运移和大颗粒的排出，且其对泥沙运移的影响作用随泥沙最大粒径的减小而减弱。

　　2.6 加气对滴头堵塞物和流道堵塞位置的影响

　　加气显著影响滴头堵塞物质量(P<0.05)，减少了 D1、D2 和 D3 浑水处理的滴头堵塞物质量，增大了 D4 和 D5 浑水处理的滴头堵塞物质量，D1、D2 和 D3 浑水加气处理的滴头堵塞物较未加气处理低 8.92%～51.00%，D4 和 D5 浑水则高 5%～29.8%(图 7)。加气增大了流道入口堵塞风险，降低了流道完全堵塞的风险(表 6)，5 种浑水加气处理的流道完全堵塞分别比未加气处理低 0.6%、7.6%、6.7%、5.0%和 4.4%，流道入口堵塞分别高 3.3%、 10.3%、10.0%、6%和 2.9%。加气增大了水流紊动性，促进了流道内大颗粒泥沙输移，小颗粒泥沙在气泡桥力作用下易发生黏结不易被水流冲出，因此加气降低了 D1～D3 浑水处理的堵塞物，增大了 D4～D5 处理的堵塞物。

　　为确定滴头堵塞参数之间的关系，计算了相对流量、灌水均匀度、堵塞率、毛管淤积物质量(M)、滴头堵塞物质量(m)、比表面积(SSA)、泥沙中值粒径(SMS)、滴头完全堵塞占比(AB)、滴头入口堵塞占比(EB)、滴头中部堵塞占比(CB)和滴头出口堵塞占比(BO)各参数之间的 Spearman 相关系数(图 8，图中**表示 P<0.01，*表示 P<0.05)。Dra 和 Cu 与毛管淤积物质量、滴头堵塞率相关系数最大，说明毛管淤积物质量和滴头堵塞率也是影响滴头抗堵塞性能的重要因素。Dra和 Cu与毛管淤积物质量呈显著正相关关系(r=0.8～ 1.0, P<0.05)，与滴头堵塞率呈显著负相关关系(r=-0.8～-1.0，P<0.05)。毛管淤积泥沙对滴头堵塞位置有极显著的影响(P<0.01)，毛管淤积泥沙中值粒径越小、泥沙比表面积越大，越易引发滴头流道发生堵塞。

　　3 讨 论

　　泥沙颗粒相互碰撞后发生凝聚、沉降行为堵塞流道被认为是造成滴头堵塞的根本原因 [8]。文献[25-26]发现水中加气能够改变毛管内水动力特性，降低壁面对泥沙颗粒的阻力，从而有效提升输沙效率，减缓泥沙沉降。文献[27]研究表明，微纳米气泡桥接过程会产生很强的长程疏水作用力，可提高粗颗粒与气泡的粘附能力，减小脱落概率，促进水流对泥沙颗粒的运移，这与本研究结果一致。此外，本研究还发现加气对泥沙输移的促进作用随着泥沙最大粒径的减小而逐渐减弱(表 4)。由于沉积的大颗粒在毛管运移时主要受内壁摩擦阻力的影响[28]，当泥沙颗粒沉积时会与毛管内壁发生碰撞导致沉积平均速度反向，而大颗粒的反向速度大于小颗粒[29]，从而降低了大颗粒的沉积、减小了毛管内壁摩擦阻力。此外，加气增大了水流波动，提高了水流流速(图 3)，进一步促进了大颗粒的运移，而细小颗粒(粒径小于 0.031mm)在其负电荷作用下易在毛管内壁发生絮凝堆积，堆积物黏附力较强不易随水流冲动，增大了内壁摩擦阻力[30-31]，减缓了泥沙在毛管内的运移。同时，细颗粒泥沙在微气泡的桥接作用下易形成絮团，增加了泥沙与气泡的接触面积，从而提升了泥沙颗粒碰撞概率、增大了颗粒与气泡的黏附概率，加速了细颗粒泥沙在毛管内壁的黏附，有效提高了毛管絮团泥沙密度[32-33]，因此，加气减缓了大颗粒泥沙沉积，而增大小颗粒沉积(表 5)。文献[34]研究表明，微气泡能够吸附水中悬浮颗粒，增大大颗粒的沉降，减缓小颗粒的沉降。这与本研究结果存在差异，文献[34]试验研究泥沙为微细粉尘颗粒，其黏性颗粒占比较高，大颗粒发生絮凝沉降后形成密实度较高、体积较大的沉积物，沉积物粘附性较强且形状不易发生破碎。而本研究中使用的河沙孔隙率较大，黏粒占比较低，发生沉降的泥沙颗粒之间黏结力较弱，较易随水流以悬移质形式向前运移[35]。因此，加气对泥沙输移的影响与泥沙性质、粒径等密切相关。

　　此外，文献[36]研究结果表明，泥沙颗粒与气泡之间存在的长程疏水力在颗粒间的团聚、气泡与颗粒的碰撞吸附过程中发挥了重要的作用，而泥沙颗粒大小是影响长程疏水力的 最 主 要 因 素 。 本 研 究 发 现 不 同 粒 径 浑 水 灌 溉 时 ， 加 气 减 少 了 D1(0 ～ 0.1mm) 、 D2(0.075~0.1mm)和D3(0.05~0.075mm)浑水处理的滴头堵塞物，增大了D4(0.031～0.05mm) 和D5(0～0.031mm)浑水处理的滴头堵塞物(图7)。这是因为加气增大浑水中大颗粒泥沙的运移的同时，导致小颗粒泥沙被滞留于毛管和流道内(表5)，随着灌水时间的增加，沉积的细小颗粒逐渐增多，而细小颗粒在微气泡的絮凝促进作用下逐渐发生聚团[37]，当随水流进入流道时，流道入口处流速骤减，使得一部分泥沙在此沉积，随着泥沙粒径减小，黏性泥沙占比逐渐增多，流道入口处沉积黏附泥沙增多，形成粗糙的壁面粘附层造成滴头堵塞 [8]。而D1～D3浑水机械组成主要以沙粒为主，水流对该粒径段(0.05～0.1mm)泥沙大颗粒输移大于小颗粒的沉降，从而减少了流道淤积物质量，而D4～D5浑水(0～0.05mm)则相反。因此，加气加剧了小颗粒在流道入口黏附，加速了流道入口堵塞(表6)是造成滴头堵塞的主要原因(图4)。文献[28]研究表明，通过冲洗可以减少细颗粒泥沙沉积，降低絮凝率，抑制堵塞物的形成从而降低流道入口堵塞的风险。建议进行毛管冲洗以减缓流道入口堵塞风险。

　　从滴头堵塞物结构观测图可以看出(图9)，加气条件下的堵塞物结构较为紧密，而未加气较为松散，这是由于加气增大了细小颗粒沉积造成的。5种浑水中，D1、D2、D3浑水大颗粒泥沙占比较高(表1)，易在重力作用下沉降形成堆积体，引起滴头堵塞。同时由于微气泡吸附泥沙颗粒的负电荷，增大了微气泡的界面电位和颗粒间的非键能，促进了细小颗粒在大颗粒表面黏附[38-39]，减小了颗粒间孔隙(图9)，增大了堆积体密实度，滴头堵塞后堆积体不易在水流剪切作用下发生破坏，因此，加气处理的滴头堵塞更为严重(图4)。而 D4、D5浑水小颗粒泥沙占比较高，细小颗粒在微气泡的桥力作用、颗粒间的吸附水膜促进作用下形成团聚体的机率相比其他粒径泥沙较大[38,40]，同时加气水流紊动性较大提升了颗粒间碰撞频率，增大了团聚体的密实度，从而加剧了滴头堵塞。D5浑水加气处理的堵塞物团聚体体积虽然最大，但形成的团聚体缺少大颗粒骨架，沉积泥沙颗粒易被水流夹带流出流道[41]，这也是该粒径浑水流量下降较为缓慢的主要原因。D2(0.075～0.1mm)浑水加气处理的堵塞物中基础骨架颗粒较多，且颗粒间孔隙率较小密实度高，不易随水流带出流道，从而表现出D2浑水加气处理的流量下降速率最快(图4)。D1浑水沙粒占比小于D2浑水，黏粒占比小于D3和D4浑水，因此，其沉积泥沙堆积和絮凝作用较为缓和，从而表现出D1 浑水流量下降最慢(图4)。

　　4 结 论

　　(1)加气和泥沙粒对滴头堵塞具有极显著影响(P<0.01)，加气加剧了滴头堵塞。泥沙粒径为 0.075～0.1mm 的浑水流量下降速率最快，粒径为 0～0.1mm 的浑水流量下降速率最慢，与未加气处理相比，加气处理滴头流量和均匀度低 9%～18.7%和 16.2%～36.4%。

　　(2)加气增大了水流波动，促进了毛管内泥沙运移和大颗粒的排出，减小淤积泥沙中值粒径，且加气对泥沙运移的影响作用随泥沙最大粒径的减小而减弱。与未加气处理相比，加气处理毛管淤积物质量降低 8.75%～31.92%，淤积泥沙中值粒径减小 8.59%～ 35.64%。

　　(3)加气对流道淤积物影响显著(P<0.05)，加气对 0.05~0.1mm 粒径浑水中大颗粒的运移作用大于小颗粒的沉降作用，减少了流道淤积物，而对 0～0.05mm 粒径浑水作用则相反：加气加剧浑水中细小颗粒在流道入口黏附，是造成滴头堵塞的主要原因。