土壤下渗湿润峰吸力和降雨初损量计算方法
【技术领域】
[0001] 本公开涉及水文学降雨产流模型的土壤参数确定，特别涉及一种基于双环入渗的 土壤下渗湿润峰吸力和降雨初损量计算方法。
【背景技术】
[0002] 我国是暴雨洪水最频发的国家之一，对洪水灾害进行及时准确的预报预警对于降 低灾害损失十分重要。洪水灾害预报预警的关键问题之一是准确模拟场次降雨产流过程。 Green-Ampt模型是物理机制清晰的产流模型，然而其湿润峰固定吸力参数难以确定，严重 限制了该模型的应用。传统的双环入渗实验仅仅能够测试土壤饱和导水率，也不能确定 Green-Ampt模型中的湿润峰固定吸力参数，难以满足Green-Ampt模型应用于流域水文学对 土壤参数的需求。因此，亟需一种可靠实用的方法，能够快速计算Green-Ampt模型中的湿润 峰固定吸力参数和降雨初损量。
【发明内容】
[0003] 为了弥补现有技术的缺失，本公开提供一种土壤下渗湿润峰吸力和降雨初损量计 算方法。
[0004] 本公开采用以下解决方案：
[0005] -种土壤下渗湿润峰吸力和降雨初损量计算方法，包括以下步骤：
[0006] 步骤1:通过双环入渗实验，测量并记录供水水量Q随入渗时间t的变化，根据公式 (1)计算实时累积入渗量h: / Q
[0007] η =- S (:1)
[0008] 其中，S表示采用的双环入渗仪的内环面积；
[0009] 步骤2:基于所述入渗时间t和实时累积入渗量h，根据公式(2)拟合参数A和参数B:
[0010] ^ - h - B ^ LN (1 -i-) A B (:2)
[0011] 步骤3:根据公式⑶计算湿润峰吸力hf:
[0012] hf = B-h〇 (3)
[0013] 其中，ho表示固定入渗水头高度；
[0014] 步骤4:根据公式(4)计算土壤饱和导水率K:
[0015] 塞二一^-- A (4)
[0016] 其中，θ〇表示土壤初期含水量，0S表示土壤孔隙度；
[0017] 步骤5:根据公式(5)计算降雨初损量Ia:
[0018] 4 = τ)Κν Ρ-Κ (3、
[0019] 其中，Ρ表示降雨强度。
[0020] 优选地，通过环刀取样测量所述土壤孔隙度0S。
[0021] 优选地，通过土壤含水量测试测定所述土壤初期含水量θ〇。
[0022] 优选地，通过最小二乘法拟合所述参数Α和参数Β。
[0023] 优选地，所述供水水量Q是供水马氏瓶的水量。
[0024] 本公开的有益效果是基于双环入渗计算土壤下渗湿润峰吸力和降雨初损量，能够 快速准确地确定Green-Ampt模型所涉及的土壤参数，进而进行洪水灾害预报预警。
【附图说明】
[0025] 通过结合附图对本公开示例性实施例进行更详细的描述，本公开的上述以及其它 目的、特征和优势将变得更加明显，其中，在本公开示例性实施例方式中，相同的参考标号 通常代表相同部件。
[0026] 图1示出了根据示例性实施例的土壤下渗湿润峰吸力和降雨初损量计算方法的流 程图；
[0027]图2示出了根据示例性实施例的入渗量与入渗时间的曲线图。
【具体实施方式】
[0028] 下面将参照附图更详细地描述本公开的优选实施例。虽然附图中显示了本公开的 优选实施例，然而应该理解，可以以各种形式实现本公开而不应被这里阐述的实施例所限 制。相反，提供这些实施例是为了使本公开更加透彻和完整，并且能够将本公开的范围完整 地传达给本领域的技术人员。
[0029] 图1示出了根据示例性实施例的土壤下渗湿润峰吸力和降雨初损量计算方法的流 程图，该方法包括以下步骤：
[0030] 步骤1，通过双环入渗实验，测量并记录供水水量Q随入渗时间t的变化，根据以下 公式(1)计算实时累积入渗量h:
[0031] h = - S (1)
[0032] 其中，S表示采用的双环入渗仪的内环面积。
[0033] 在实际操作中，优选采用双环入渗仪进行双环入渗实验。双环入渗仪是本领域中 常用的测量水渗入土壤的渗透速度的设备。标准的双环入渗仪包括直径不同的外环和内 环。垂直渗透水流向边缘时，外环可以起到隔离的作用。由于内环中的水是垂直流动的，因 而测量仅限于在内环中进行。
[0034] 作为一种优选方案，采用马氏瓶作为供水设备，因此供水水量Q是指供水马氏瓶的 水量，当采用其他供水设备时，供水水量Q是指所采用的供水设备的水量。
[0035] 步骤2,基于步骤1中的入渗时间t和实时累积入渗量h，根据以下公式(2)，拟合参 数A和参数B:
[0036] - = k - B * LN (} H-) ^ B (2)
[0037] 在实际实施过程中，可以通过最小二乘法编制计算机程序或者利用通用数据处理 软件(例如TableCurve 2D)，拟合参数A和参数B。
[0038] 步骤3,根据以下公式(3)计算湿润峰吸力hf:
[0039] hf = B-ho (3)
[0040] 其中，ho表示固定入渗水头高度。
[0041] 步骤4,根据以下公式(4)计算土壤饱和导水率K:
[0042] K = -~- Λ. (4)
[0043]其中，θ〇表示土壤初期含水量，可通过土壤含水量测试来测定，0S表示土壤孔隙度， 可通过环刀取样来测定。
[0044]步骤5,根据以下公式(5)计算降雨初损量Ia:
[0045] 4 = Μ6>v _ 怂)-^―： p- κ (5)
[0046] 其中，P表示降雨强度。根据实际应用，P可以通过雨量计测量，也可以是设计的降 雨强度。
[0047] 应用示例
[0048] 本示例提供一种基于野外双环入渗实验的土壤下渗湿润峰吸力和降雨初损量计 算方法，包括以下步骤：
[0049] 步骤1:进行双环入渗实验，以1秒为采样周期，利用挂式电子秤连接电脑串口测量 并记录供水马氏瓶的水量随入渗时间t的变化，根据公式（1)计算实时累积入渗量h，得到的 实时累积入渗量h和入渗时间t的曲线图如图2所示。
[0050] 步骤2:基于步骤1中的实时累积入渗量h和入渗时间t的数据，根据公式(2)，通过 最小二乘法拟合公式(2)中的参数A和参数B，得到A = 32.45s/cm，B = 25.91cm，R2 = 0.9997， 其中R2表示拟合确定系数。
[0051 ] 步骤3 :测量得到固定入渗水头高度ho = 10cm，根据公式(3)计算湿润峰吸力hf = 15.91cm〇
[0052] 步骤4:通过土壤含水量测试测量得到土壤初期含水量θ〇 = 〇,通过环刀取样测得 的土壤孔隙度0S = 36%，根据公式(4)计算土壤饱和导水率Κ = 399.22mm/hr。
[0053]步骤5:根据公式(5)计算降雨初损量： 399
[0054] / =15.91 *(0.36-0)--二二一 a V-399.22。
[0055] 上述技术方案只是本发明的一种实施例，对于本领域内的技术人员而言，在本发 明公开了应用方法和原理的基础上，很容易做出各种类型的改进或变形，而不仅限于本发 明上述具体实施例所描述的方法，因此前面描述的方式只是优选的，而并不具有限制性的 意义。
【主权项】
1. 一种±壤下渗湿润峰吸力和降雨初损量计算方法，包括W下步骤： 步骤1:通过双环入渗实验，测量并记录供水水量Q随入渗时间t的变化，根据公式（1)计 算实时累积入渗量h:(1) 其中，S表示采用的双环入渗仪的内环面积； 步骤2:基于所述入渗时间t和实时累积入渗量h，根据公式(2)拟合参数A和参数B:(2) 步骤3:根据公式(3)计算湿润峰吸力hf: hf = B-ho (3) 其中，ho表示固定入渗水头高度； 步骤4:根据公式(4)计算±壤饱和导水率K:(4) 其中，θ〇表示±壤初期含水量，0S表示±壤孔隙度； 步骤5:根据公式巧)计算降雨初损量la:(5) 其中，P表示降雨强度。2. 根据权利要求1所述的±壤下渗湿润峰吸力和降雨初损量计算方法，其中通过环刀 取样测量所述±壤孔隙度0S。3. 根据权利要求1所述的±壤下渗湿润峰吸力和降雨初损量计算方法，其中通过±壤 含水量测试测定所述±壤初期含水量9〇。4. 根据权利要求1所述的±壤下渗湿润峰吸力和降雨初损量计算方法，其中通过最小 二乘法拟合所述参数A和参数B。5. 根据权利要求1所述的±壤下渗湿润峰吸力和降雨初损量计算方法，其中所述供水 水量Q是供水马氏瓶的水量。
【专利摘要】一种土壤下渗湿润峰吸力和降雨初损量计算方法，包括以下步骤：通过双环入渗实验，测量并记录供水水量随入渗时间的变化，计算实时累积入渗量；基于入渗时间和实时累积入渗量，拟合必要参数；基于拟合参数和固定入渗水头高度，计算湿润峰吸力；基于拟合参数、土壤初期含水量和土壤孔隙度，计算土壤饱和导水率；基于湿润峰吸力、土壤饱和导水率、降雨强度、土壤初期含水量和土壤孔隙度，计算降雨初损量。
【IPC分类】G01N15/08
【公开号】CN105547957
【申请号】CN201510917405
【发明人】李军, 杨默远, 王中根, 刘昌明
【申请人】中国科学院地理科学与资源研究所
【公开日】2016年5月4日
【申请日】2015年12月10日