河西走廊中段绿洲是我国重要的农业生产区，也是国内单季粮食产量较高的区域之一。近35年以来，人为开垦使绿洲面积增加了1.5倍^[1-2]。由于不同的开垦年限以及人为耕作措施的影响，形成了不同类型的绿洲农田土壤，其特殊的土壤母质，气候条件和管理方式对土壤孔隙结构和水分运动过程产生着深刻的影响^[3-5]。荒漠绿洲农田土壤由于受到干湿交替、冻融作用和土壤特性等非生物因素和动物洞穴、植物根系等生物因素的影响，大孔隙广泛存在^[3]。荒漠绿洲区土壤水分运动的研究以往多关注均匀流^[4]，也有相关学者通过染色示踪技术表明了荒漠绿洲地区存在优先流现象，发现大孔隙流现象是荒漠绿洲区土壤水分运动的重要形式之一^[5-6]。土壤大孔隙虽然只占土壤体积的0.1%~5%，但却可传导90%的水流通量，对土壤水分运动起着关键作用^[7]。绿洲农田灌溉水分利用效率提高、养分流失和地下水污染与土壤大孔隙引起的优先流现象有着密切关系^[8]。大孔隙作为土壤优先流通道，将水分传至土壤深层，影响着荒漠绿洲区的生态水文过程^[9-11]。

目前对大孔隙的研究方法主要有渗透仪法、染色示踪法、CT扫描技术以及溶质穿透曲线法等^[12-14]。染色示踪法观测简单直观，适合在农田进行土壤大孔隙流现象研究。Zhang等^[15]将染色技术应用到水稻田土壤裂隙的研究，李胜龙等^[16]利用染色技术研究了稻田－田埂过渡区的优先流特征。通过染色可以观察到土壤大孔隙流形态特征，但难以获得土壤三维孔隙结构特征^[17]，缺乏对土壤大孔隙的定量化研究。近年来CT扫描技术被国内外学者广泛应用在土壤大孔隙结构研究中^[18-21]，CT扫描是一种无损的3D成像技术，能够在不破坏土壤结构的情况下，获取土壤大孔隙结构^[22]，进一步分析得到土壤大孔隙结构的特征参数。Hu等^[18]利用CT扫描技术研究了青海湖流域不同土地利用类型的大孔隙结构，结果表明土壤类型是影响土壤大孔隙特征的重要因素之一。

目前关于土壤大孔隙特征的研究多集中在水分条件较好的草地、林地、水稻田^{[9, 11, 14]}等，而对干旱区绿洲农田的大孔隙特征的研究相对较少，特别是绿洲农田的土壤大孔隙结构特征如何？大孔隙的产生对绿洲农田水分入渗的影响如何？基于上述认识，本文将CT扫描和原位染色示踪技术相结合，对干旱区典型绿洲农田的大孔隙特征进行了较为全面深入的探究，旨在揭示大孔隙对荒漠绿洲农田水分入渗过程的影响，提高荒漠绿洲区地表水与地下水转换规律的认识。

1 材料与方法 1.1 研究区概况

研究区地处甘肃省河西走廊中段(99°51′~100°30′E，38°57′~39°42′N，海拔1 374 m)。该区是典型的干旱荒漠气候，年平均气温7.6℃，年平均日照时数3 051h，年平均风速3.2 m·s^-1，多年平均降水量为117 mm，蒸发量2 388 mm，7-9月降水量约占全年降水量的60%。研究区主要景观类型有外围荒漠、荒漠绿洲过渡带和中心绿洲区。黑河中游是西北干旱区典型的绿洲灌溉农业，农业种植主要依赖于黑河地表水和地下水灌溉，以种植玉米为主。近几十年来，荒漠边缘被逐步开垦，形成了不同开垦年限的绿洲农田^[23]，土壤类型为灌漠土。其中新开垦的边缘绿洲农田土壤有机质低，保水性能差，土壤为砂土；经过百年尺度耕作的老绿洲农田，由于人为施肥灌溉，其土壤肥力较高，土壤质量好；分布在北山山前的绿洲，受地层条件和山洪洪水所携带泥沙淤积，开垦形成的山前绿洲农田，土壤耕作层主要是黏土层，土壤为黏土(图 1)。

1.2 试验设计

2019年7-10月在研究区选取了三种典型农田：老绿洲农田、边缘绿洲农田、山前绿洲农田开展试验。利用高30 cm，内径19 cm的PVC管在三个样地进行原状土柱的采集，每个样地取3个重复，共计9个土柱。在原状土柱转运过程中注意对样品的保护，确保其土壤结构的完整性，为避免土壤水分的蒸发，采样结束后利用保鲜膜将土柱的两端密封。采集的原状土柱在临泽县人民医院利用CT机(Optima CT520，GE，USA)进行扫描，激发电压为100kvp，电流为100mA，扫描图像的水平分辨率为0.469 mm，垂直分辨率为0.625 mm，由于受到CT扫描分辨率的限制，本研究中绿洲农田土壤大孔隙孔径大于0.47 mm。为避免作物根系对土壤结构的影响，原状土柱的采集和入渗试验均在玉米农田的行间空地上进行。

在采集原状土柱的样点附近，开展单环染色入渗试验，单环直径20 cm，设计灌水量为100 mm，与当地实际灌水定额保持一致。入渗试验中，采用亮蓝溶液作为染色剂，浓度为4 g·L^-1，利用马氏瓶供水，保持定水头5 cm。入渗过程中记录马氏瓶水位刻度随时间变化(实验刚开始时，每3 min记录一次水面位置，30 min以后每10 min记录一次)，在每个样地进行3个重复。入渗完成24 h后，入渗湿润锋基本稳定，垂直开挖土壤剖面，剖面深度为40 cm。将剖面修理平整后，用相机对染色剖面进行拍照，为避免阳光直射，采取遮光拍照获取染色剖面图像。并在0~30 cm土层每隔10cm采集土壤样品，测定土壤容重、饱和含水量、总孔隙度和机械组成等基本性质。采用环刀法测定土壤容重、土壤总孔隙度和饱和含水量等，土壤机械组成采用激光粒度仪(Mastersizer 2000)测定。

1.3 数据处理

土壤染色图像经过Photoshop CS6进行裁剪处理、几何校正，处理好的图像在Image-J软件中进行图像二值化处理，然后利用ArcMap10计算土壤剖面染色特征参数(最大入渗深度、优先流长度指数、土壤染色面积比)。染色面积比是指土壤染色区域面积与所选取的整个剖面面积的比值。基质流入渗深度为染色面积比 > 80%以上的土层深度^[23]。优先流长度指数是整个土壤染色剖面上下层染色面积比的绝对差值，计算公式为$ \mathrm{L}=\sum_{i=1}^{n}|\operatorname{Dc}(i+1)-\operatorname{Dc}(i)|$，将土壤染色剖面划分成n层，Dc(i)表示第i层的染色面积比^[24]。

为消除边缘效应的影响，CT扫描后的原始图像去除上下各20张图像，并选择中间12 cm×12 cm的感兴趣区域进行分析。采用目视解译分析方法，利用Image-J对图像进行阈值分割，对照每个土柱CT扫描的原始图像，进行对比选取合理的阈值，进行二值化阈值分割处理。利用Image-J中的3D viewer模块对图像进行三维结构重建，构建土壤大孔隙三维结构图，利用Image-J中的BoneJ插件构建三维骨架结构，其中孔隙数、孔隙度、大孔隙面积密度、平均孔隙体积、分形维数、分支密度、连接点密度、平均孔隙直径、连通性密度均可通过Image-J直接计算得出^[25-26]。

利用Excel对数据进行计算处理，利用SPSS 20.0对数据进行单因素方差分析、显著性检验和相关性分析，利用Origin 2018进行制图。

2 结果 2.1 绿洲农田土壤基本物理性质

表 1显示了绿洲农田土壤的基本物理性质，0~30 cm土层的三种绿洲农田砂粒含量存在显著性差异(P < 0.05)，但不同土层间的差异并不显著，老绿洲、边缘绿洲农田的砂粒含量分别为山前绿洲的2.58倍、3.43倍。根据美国农业部土壤质地分类标准，山前绿洲农田质地为粉黏壤土，老绿洲农田质地为壤质砂土，边缘绿洲农田质地为砂土。边缘绿洲和山前绿洲农田的土壤容重显著大于老绿洲农田(P < 0.05)，但老绿洲农田的总孔隙度最大，为45%。

2.2 绿洲农田土壤大孔隙特征

由图 2可知，三种绿洲农田的土壤大孔隙大多集中于农田的耕作层(0~20 cm)。老绿洲农田的孔隙发育好，土壤大孔隙结构复杂，土壤孔隙的连续性较好；老绿洲农田土壤的大孔隙度明显高于边缘绿洲农田和山前绿洲农田，且在20 cm土层以下，仍有一定数量的大孔隙分布。边缘绿洲农田大孔隙相对较少，且其连通性性低，在10~20 cm土层中，有少量的连续性孔隙分布，可能是由腐烂的作物根系形成的。不同类型的绿洲农田土壤大孔隙结构呈现明显的异质性特征。

经过Image-J软件处理计算得到不同绿洲农田土壤大孔隙特征参数(表 2)，土壤大孔隙度表现为：老绿洲农田(0.54%) > 山前绿洲农田(0.22%) > 边缘绿洲农田(0.14%)。老绿洲农田的大孔隙度、大孔隙数、分形维数、分支密度、连接点密度、连通性密度和山前绿洲农田、边缘绿洲农田之间存在显著性差异(P < 0.05)。土壤大孔隙分形维数与土壤结构密切相关，老绿洲农田的分形维数特征值显著高于其他两种农田(P < 0.05)。山前绿洲、老绿洲、边缘绿洲的土壤大孔隙平均直径分别为3.27 mm、3.14 mm、3.58 mm，不同绿洲农田的平均土壤大孔隙体积、平均孔隙直径差异不显著(P > 0.05)。较低的分支密度、连接点密度和连通性密度意味着土壤大孔隙连通性较差，而老绿洲农田的分支密度、连接点密度以及连通性指数最大，因此其大孔隙连通性最好。

如图 3所示，三种绿洲农田的大孔隙面积密度均表现出随土壤深度增加而波动性减小，并在不同深度处存在峰值现象。边缘绿洲的大孔隙面积密度在10~15 cm土层深度处呈现波动性峰值，可能是由于该层次土壤中存在较大裂缝引起的。山前绿洲农田和边缘绿洲农田20 cm土层以下的大孔隙面积密度较小，而老绿洲农田的大孔隙面积密度与上层差异较小，表明在20 cm土层下依旧存在着一定数量的大孔隙。老绿洲农田与边缘绿洲农田、山前绿洲农田的平均大孔隙面积密度存在显著性差异(P < 0.05)。

2.3 绿洲农田土壤水分入渗特征

在土壤水分入渗过程中，通常用初始入渗速率、稳定入渗速率、累计入渗量等指标来衡量土壤的入渗性能。由图 4可知，在土壤水分入渗过程中，边缘绿洲的累计入渗量大于老绿洲和山前绿洲，山前绿洲、老绿洲、边缘绿洲90 min累计入渗量分别为26.2 mm、33.9 mm、58.3 mm，表明山前绿洲农田的入渗速率最慢，入渗性能低于其他两种农田。

三种绿洲农田的入渗指标存在差异(图 5)，稳定入渗速率表现为边缘绿洲农田(0.48 mm·min^-1) > 老绿洲农田(0.28 mm·min^-1) > 山前绿洲农田(0.16 mm·min^-1)，有显著性差异(P < 0.05)。初始入渗速率表现为边缘绿洲农田(3.02 mm·min^-1) > 山前绿洲农田(2.47 mm·min^-1) > 老绿洲农田(1.99 mm·min^-1)，边缘绿洲和山前绿洲的初始入渗速率显著大于老绿洲农田(P < 0.05)。边缘绿洲农田的入渗能力较高，其稳定入渗速率为老绿洲的1.71倍、山前绿洲的3倍，这主要是由于土壤机械组成决定的(表 1)。

2.4 绿洲农田土壤水分运动的染色特征

如图 6所示，三种绿洲农田的土壤染色特征存在差异，边缘绿洲农田的入渗以基质流为主，基质流入渗深度占其总入渗深度比例高，达到79%，入渗深度浅；而老绿洲最大入渗深度大，土壤水分分布的空间异质性强，优先流现象明显。由图 7可知，三种农田的染色面积比随土层深度的增加而波动下降，总染色面积比为：老绿洲农田(56%)显著大于边缘绿洲和山前绿洲(P < 0.05)。最大入渗深度为：老绿洲(33 cm) > 山前绿洲(21.8 cm) > 边缘绿洲(18.4 cm)。老绿洲、山前绿洲、边缘绿洲的优先流长度指数分别为310、200、174，优先流长度指数越大表示其土壤水分运动的优先流现象明显。老绿洲农田的优先流长度指数和最大入渗深度大于山前绿洲和边缘绿洲，土壤水分运动以大孔隙流为主。

2.5 入渗特征的影响因素

由表 3可知，干旱区绿洲农田稳定入渗速率与砂粒含量呈显著正相关(P < 0.01)，砂粒含量越高，其稳定入渗速率越大，与黏粒、粉粒含量呈显著负相关(P < 0.01)。最大入渗深度、优先流长度指数与土壤大孔隙度、分形维数、连接点密度、连通性密度呈显著正相关(P < 0.01)，土壤大孔隙度越高，且连通性越好，其土壤水分运动的优先流现象越明显。

3 讨论

干旱区绿洲农田受到灌溉的干湿交替、作物根系和动物活动等因素的影响，绿洲农田的土壤大孔隙特征存在着异质性^[23]。绿洲农田的土壤大孔隙大多集中于0~20 cm耕作层，这与Zhang等^[25]对红壤区水稻田研究结果一致，在20 cm以下为犁底层，受耕作影响小，土壤容重较大，土壤的大孔隙度低。老绿洲的土壤大孔隙结构特征较边缘绿洲、山前绿洲复杂，且土壤大孔隙度、分形维数、连接点密度显著高于边缘绿洲和山前绿洲。农事活动耕作对土壤大孔隙结构发育有一定的影响，蔡太义等^[27]利用CT技术对煤矿区不同复垦年限的土壤孔隙结构的研究表明，随着复垦年限的延长，土壤大孔隙度、孔隙数量总体呈现逐渐增加的趋势。老绿洲农田由于开垦时间长，在灌溉、人为耕作和施肥等因素的影响下^[24]，其黏粉粒含量要高于边缘绿洲农田，土壤质量较高，土壤孔隙发育良好，其孔隙结构也更加复杂。边缘绿洲农田开垦时间较短，砂粒含量大于80%，容重较大，土壤有机质含量低，发育程度较差，土壤大孔隙结构相对简单。山前绿洲农田由于其耕作层土壤主要是洪水冲积物，黏粉粒含量较高，达到75%，容重较大，但开垦时间短，约10年，其孔隙发育较差，大孔隙度与边缘绿洲农田差异不显著。绿洲农田土壤大孔隙结构发育与土壤质地、农事耕作活动密切相关，绿洲化过程改变了土壤的理化性质和水力传导度，土壤大孔隙度增加、孔隙结构逐步形成^{[3, 11]}。

土壤质地是影响土壤水分入渗过程的重要因素^[28]，解文艳和樊贵盛^[29]对潇河灌区农田入渗能力的研究表明，土壤质地对土壤入渗能力的影响显著。边缘绿洲农田土壤质地为砂土，老绿洲农田为壤砂土，山前绿洲农田为粉黏壤土(表 1)，在整个入渗过程中，边缘绿洲农田的入渗能力最强。由稳定入渗速率和累计入渗量结果表明，边缘绿洲农田显著大于老绿洲、山前绿洲，主要是因为其砂粒含量较高。本研究与谭丽丽^[30]的研究结果一致，其通过对滨州农田入渗的研究结果表明，土壤的稳渗速率、累积入渗量与黏粒含量呈显著负相关；吕刚等^[31]在辽西北风沙地的研究也表明稳定入渗率与土壤中粉粒含量呈负相关。同时，土壤水分入渗的性能与孔隙度、容重表现出很强的相关性^[4]，魏恒等^[32]对塔里木河上游绿洲不同植被类型土壤水分入渗的研究表明，土地利用类型对入渗的影响较大。在干旱区绿洲农田，土壤稳定入渗速率主要受土壤质地的影响。

土壤大孔隙是形成土壤优先流的主要原因之一^[24]，王发等^[33]研究表明耕作方式改变了土壤结构，形成了不同特征的土壤大孔隙，进而影响着土壤优先流运动特征。边缘绿洲和山前绿洲的土壤大孔隙度低，大孔隙连通性差，入渗以基质流为主，平均最大入渗深度在20 cm左右，优先流长度指数较低，大孔隙流现象不明显。而老绿洲农田的大孔隙度、分形维数、分支密度、连接点密度等大孔隙特征参数优于边缘绿洲和山前绿洲，大孔隙连通性好，土壤大孔隙结构复杂；最大入渗深度和优先流长度指数显著大于山前绿洲和边缘绿洲，大孔隙特征参数与优先流程度表现出很强的相关性，更易产生大孔隙流现象。郭会荣等^[34]通过溶质穿透实验表明，大孔隙度越高，优先流对溶质穿透曲线的贡献越大。土壤大孔隙为土壤水分的运动提供了优先流通道，使其能够绕过部分土壤基质向土壤深处运动，不同的土壤大孔隙结构特征影响着绿洲农田的水分入渗过程。干旱区绿洲农田受到农事活动措施(开垦年限、耕作和灌溉等)及自然因素(土壤母质、干湿冻融交替等)的综合作用^[35]，改变了土壤理化性质、水力传导度及其土壤结构的演变过程，形成了不同土壤类型的绿洲农田，土壤大孔隙度增加、孔隙结构形成，影响着绿洲农田土壤水分入渗的途径和过程。

4 结论

干旱区绿洲农田的土壤大孔隙特征及土壤水分入渗特征存在差异。老绿洲农田大孔隙度、大孔隙数、分形维数、连接点密度、连通性密度显著大于山前绿洲、边缘绿洲(P < 0.05)。三种绿洲农田的稳定入渗速率存在着显著性差异(P < 0.05)，稳定入渗速率和累计入渗量表现为：边缘绿洲农田 > 老绿洲农田 > 山前绿洲农田。老绿洲农田的优先流长度指数和最大入渗深度大于山前绿洲和边缘绿洲，其土壤水分运动优先流现象明显。绿洲农田的土壤水分入渗特征与土壤质地、土壤大孔隙特征有较强的相关性，土壤质地对稳定入渗率影响较大，而大孔隙在土壤水分运移过程发挥着重要作用。绿洲化过程改变了土壤的理化性质和入渗特性，老绿洲土壤大孔隙度增加、孔隙结构形成，孔隙连通性好，更易产生优先流，土壤水分运动的空间异质性程度高。