土壤水分特征曲线的压力膜仪检测实验报告

土壤水分特征曲线的压力膜仪检测实验报告一、实验材料与方法1.1供试土壤样品采集与预处理本次实验供试土壤采自华北平原典型农田耕作层0~20cm深度区域，土壤类型为潮土。采样过程采用五点混合法，共采集原状土柱5个，同时采集扰动土样20kg，装入密封袋后带回实验室处理。原状土柱采用PVC环刀（容积100cm³，高度5cm）采集，采样过程保持土壤结构完整，避免挤压扰动；扰动土样剔除根系、砾石等杂质后，自然风干、研磨后过2mm筛，密封保存备用。实验前对土壤基本理化性质进行测定：采用比重瓶法测定土壤颗粒密度为2.68g·cm⁻³，环刀法测定田间原状土容重为1.35g·cm⁻³，质地分类为粉质壤土，其中砂粒（2~0.02mm）含量28.7%、粉粒（0.02~0.002mm）含量56.2%、黏粒（<0.002mm）含量15.1%，有机质含量12.3g·kg⁻¹，pH值为8.1。1.2实验仪器与工作原理本次实验采用美国SoilmoistureEquipmentCorp生产的1500F1型压力膜仪，仪器主要由压力腔室、调压系统、陶土板、排水收集装置四部分组成。核心工作原理为：将饱和土壤样品放置在具有特定进气值的陶土板上，密封压力腔后通入压缩氮气，通过调压系统精确控制腔室内压力值，土壤中超过基质吸力对应的水分会透过陶土板排出，直至土壤水势与外部压力达到平衡，此时测定土壤含水量，即可得到该吸力下的土壤水分特征点。实验配套陶土板进气值分别为1bar（适用于0~1bar吸力范围）和15bar（适用于1~15bar吸力范围），压力控制精度为±0.01bar，含水量测定采用精度0.001g的电子分析天平，烘干温度控制采用105℃恒温鼓风干燥箱。1.3实验设计与操作流程实验设置原状土与扰动土两组处理，每组处理设置3个重复。具体操作流程如下：（1）样品饱和：将环刀盛装的土壤样品（原状土直接使用，扰动土按容重1.35g·cm⁻³装填）放置于盛有去离子水的托盘内，水面高度低于环刀上沿1cm，让土壤通过毛细作用自下而上饱和，饱和时间为48h，饱和后取出擦干环刀表面水分，称重记录饱和重量。（2）低压段（0~1bar）测定：将饱和样品放置于1bar陶土板上，陶土板预先浸泡24h排除内部气泡。密封压力腔室后，依次设置压力梯度为0.1、0.2、0.3、0.5、0.8、1.0bar，每个压力梯度下平衡时间为24~48h，直至排水口无连续水滴流出视为达到平衡。平衡后快速取出样品，擦干表面水分后称重，记录重量后放回压力腔继续下一个梯度测定。（3）高压段（1~15bar）测定：低压段测定完成后，更换15bar陶土板，依次设置压力梯度为2、3、5、7、10、15bar，每个压力梯度平衡时间为48~72h，达到平衡后同样称重记录重量。（4）烘干称重：全部压力梯度测定完成后，将土壤样品放入105℃干燥箱内烘干至恒重，称量干土重量，计算每个压力梯度下的土壤质量含水量与容积含水量。1.4数据处理与模型拟合土壤容积含水量计算公式为：θ_v=(W_w-W_d)/(V×ρ_w)×100%其中θ_v为容积含水量（cm³·cm⁻³），W_w为湿土重量（g），W_d为干土重量（g），V为环刀容积（cm³），ρ_w为水的密度（1g·cm⁻³）。采用vanGenuchten模型对实验测定的水分特征曲线进行拟合，模型表达式为：θ(h)=θ_r+(θ_s-θ_r)/[1+(α|h|)^n]^(1-1/n)其中θ(h)为基质吸力h对应的容积含水量，θ_s为饱和含水量，θ_r为残余含水量，α为与进气值相关的拟合参数，n为与孔隙分布均匀性相关的形状参数。拟合过程采用OriginPro2023软件的非线性拟合模块完成，拟合精度采用决定系数R²和均方根误差RMSE评价。二、实验结果与分析2.1土壤水分特征曲线基本变化规律实验测定的原状土与扰动土水分特征曲线变化趋势一致，均表现为随着基质吸力升高，土壤含水量逐渐降低，且降低速率呈现“快-慢-稳”的特征。在低吸力段（<0.3bar），土壤含水量随吸力升高快速下降，原状土容积含水量从饱和时的0.462cm³·cm⁻³降至0.318cm³·cm⁻³，降幅达31.2%；扰动土饱和含水量为0.458cm³·cm⁻³，降至0.297cm³·cm⁻³，降幅为35.2%。该阶段排出的水分主要是土壤中的大孔隙重力水，说明低吸力段水分变化主要受土壤大孔隙结构控制。中吸力段（0.3~5bar），土壤含水量下降速率明显放缓，原状土含水量从0.318cm³·cm⁻³降至0.187cm³·cm⁻³，降幅为41.2%；扰动土含水量从0.297cm³·cm⁻³降至0.156cm³·cm⁻³，降幅为47.5%。该阶段土壤中毛管孔隙中的水分逐步排出，是土壤有效水的主要分布区间，不同处理的含水量差异逐渐显现，原状土含水量始终高于扰动土，说明原状土的毛管孔隙连通性更好，持水能力更强。高吸力段（5~15bar），土壤含水量下降趋于平缓，基质吸力从5bar升高到15bar时，原状土含水量仅从0.187cm³·cm⁻³降至0.122cm³·cm⁻³，降幅为34.8%；扰动土从0.156cm³·cm⁻³降至0.098cm³·cm⁻³，降幅为37.2%。该阶段排出的水分主要是受土壤颗粒吸附的束缚水，水分移动难度大，15bar对应的含水量通常被认为是土壤凋萎系数，本次实验测定原状土凋萎系数为0.122cm³·cm⁻³，扰动土为0.098cm³·cm⁻³，表明原状土能够为作物提供更多的有效水分。2.2原状土与扰动土水分特征曲线差异分析对比两组处理的测定结果，发现原状土与扰动土的水分特征曲线存在显著差异，主要体现在三个方面：一是饱和含水量差异，原状土饱和含水量为0.462cm³·cm⁻³，略高于扰动土的0.458cm³·cm⁻³，这是由于原状土中存在部分较大的生物孔隙和团聚体间孔隙，而扰动土经过研磨装填后，大孔隙被破坏，小孔隙比例增加，整体孔隙度略有下降。二是进气值差异，原状土的进气值约为0.08bar，而扰动土的进气值约为0.12bar，说明扰动土的平均孔隙直径更小，需要更高的压力才能将孔隙中的空气压入。进气值的差异直接导致低吸力段水分释放速率不同，在0.1bar吸力下，原状土含水量降低了0.078cm³·cm⁻³，而扰动土仅降低了0.059cm³·cm⁻³，进一步验证了原状土大孔隙比例更高的特征。三是有效水含量差异，将田间持水量（0.3bar对应的含水量）与凋萎系数（15bar对应的含水量）的差值作为有效水含量，计算得到原状土有效水含量为0.196cm³·cm⁻³，扰动土为0.199cm³·cm⁻³，二者差异不大，但分布区间不同：原状土在低吸力段（<1bar）的有效水占比为62.3%，而扰动土仅为54.7%，说明原状土中的有效水更易被作物吸收利用，而扰动土中更多有效水存在于较小的孔隙中，作物吸收难度更大。2.3vanGenuchten模型拟合结果采用vanGenuchten模型对两组处理的实验数据进行拟合，拟合参数结果见表1。|处理类型|θ_s(cm³·cm⁻³)|θ_r(cm³·cm⁻³)|α(bar⁻¹)|n|R²|RMSE||----------|----------------|----------------|-----------|---|----|------||原状土|0.465|0.087|0.042|1.46|0.992|0.0087||扰动土|0.460|0.065|0.035|1.38|0.989|0.0094|拟合结果显示，两组处理的决定系数R²均大于0.98，均方根误差RMSE均小于0.01，说明vanGenuchten模型能够很好地描述本次实验测定的土壤水分特征曲线。从拟合参数来看，原状土的残余含水量θ_r高于扰动土，这是因为原状土中团聚体结构更稳定，包裹的束缚水含量更高；α值原状土大于扰动土，与实验测定的进气值结果一致，反映了原状土孔隙更大的特征；n值原状土更高，说明原状土的孔隙分布更不均匀，大孔隙和小孔隙的差异更显著。2.4实验误差来源分析本次实验的平行样相对误差均控制在5%以内，符合土壤物理实验精度要求，但仍存在部分系统误差，主要来源包括：一是平衡时间误差，不同质地土壤在不同压力下的平衡时间存在差异，本次实验采用固定平衡时间，部分黏粒含量较高的样品可能未达到完全平衡，导致测定的含水量略高于实际值，尤其在高吸力段误差更为明显。二是样品蒸发损失，在取出样品称重过程中，土壤表面水分会快速蒸发，尤其在高吸力段样品含水量较低时，蒸发损失带来的相对误差更大，实验过程中采用快速称重、每次称重时间控制在30s以内的方式降低该误差，估算蒸发损失带来的误差约为1.2%。三是陶土板滞后效应，陶土板本身的水分特征曲线存在滞后现象，即湿润和干燥过程的水分传导率不同，本次实验采用的是干燥过程测定，若进行湿润过程测定，结果会存在一定差异。四是容重变化影响，扰动土在多次加压过程中可能发生轻微压实，导致容重升高，孔隙结构改变，而原状土结构更稳定，容重变化较小，这也是两组处理差异的原因之一。三、实验结论与应用参考3.1实验可靠性评价本次实验采用压力膜仪测定的土壤水分特征曲线重复性好，模型拟合精度高，能够准确反映供试潮土的持水特性。测定得到的原状土田间持水量为0.318cm³·cm⁻³，凋萎系数为0.122cm³·cm⁻³，有效水含量为0.196cm³·cm⁻³，与华北平原潮土的已有研究结果基本一致，说明实验操作规范，结果可靠。对比原状土与扰动土的结果可以发现，土壤结构扰动会显著改变水分特征曲线形态，尤其在低吸力段差异更为明显，因此在进行农田水文模拟、灌溉制度制定等应用时，应优先采用原状土测定结果，避免使用扰动土带来的误差。3.2农业生产应用参考本次实验结果可为华北平原潮土区的农田水分管理提供数据支撑：一是灌溉阈值确定，当土壤含水量降低至0.18cm³·cm⁻³（对应基质吸力约5bar）时，作物吸水难度显著增加，应及时进行灌溉，避免发生旱害；二是灌溉定额计算，根据有效水含量计算，0~20cm耕作层土壤最大可储存有效水量约为39.2mm，因此单次灌溉定额控制在30~40mm较为合理，既能够满足作物需求，又可以避免深层渗漏损失；三是土壤改良效果评价，通过对比改良前后的水分特征曲线变化，可以定量评价秸秆还田、有机肥施用等措施对土壤持水能力的提升效果。3.3实验方法优化建议为进一步提高压力膜仪测定精度，建议在后续实验中进行以下优化：一是