Method for Predicting Soil Hydraulic Properties Based on Upward Infiltration Experiment of Cutting Ring Soil Sample

Abstract

The present disclosure relates to a method for predicting soil hydraulic properties based on an upward infiltration experiment of a cutting ring soil sample. According to the new design approach, based on an inverse process of one-dimensional water upward infiltration analytical solution of homogeneous dry soils under a boundary condition of a constant pressure, by recording the data regarding cumulative infiltration amounts and infiltration times during the upward infiltration process of the cutting ring soil sample and a time for the wetting front to reach the upper surface of the cutting ring soil sample, a saturated hydraulic conductivity of the cutting ring soil sample is measured, and the soil hydraulic properties are quickly obtained. The present disclosure overcomes the defects that the methods in the prior art for measuring soil hydraulic properties are time consuming, low in efficiency, and not suitable for rapid acquisition of large area in the field, avoids the problems of multiple solutions and non-convergence of parameters, and improves the efficiency of acquiring soil hydraulic properties, thereby providing an effective means for the investigation of soil hydraulic properties in a large area in the field.

Claims

1. A method for predicting soil hydraulic properties based on an upward infiltration experiment of a cutting ring soil sample, comprising the following steps: step A: during a one-dimensional water upward infiltration experiment of homogeneous dry soil for the cutting ring soil sample under a boundary condition of a constant pressure h.sub.p, recording a change of a cumulative infiltration amount I of the cutting ring soil sample with infiltration time t, and an infiltration time t.sub.e for a wetting front in soil to reach an upper surface of the cutting ring soil sample, measuring a saturated hydraulic conductivity K.sub.s of the cutting ring soil sample, and proceeding to step B; step B: according to the change of the cumulative infiltration amount I of the cutting ring soil sample over the infiltration time t, obtaining a cumulative infiltration amount I.sub.te, and an advancing distance z.sub.fte of the wetting front in soil at a time t.sub.e, as well as a ratio of z.sub.fte to I.sub.te, namely, an average soil moisture content increase U in a wetting region of the cutting ring soil sample, and proceeding to step C; step C: obtaining a shape coefficient a of a soil moisture profile of the cutting ring soil sample based on a saturated soil moisture content θ.sub.s, an initial soil moisture content θ.sub.i and a residual soil moisture content θ.sub.r of the cutting ring soil sample, as well as the average soil moisture content increase U in the wetting region of the cutting ring soil sample, and proceeding to step D; and step D: according to the shape coefficient a of the soil moisture profile of the cutting ring soil sample, obtaining a soil pore size distribution index n of the cutting ring soil sample.

2. The method for predicting soil hydraulic properties based on an upward infiltration experiment of a cutting ring soil sample according to claim 1, wherein after step D is executed, the method further comprises the following steps: step E: obtaining an effective saturated hydraulic conductivity K.sub.e of the cutting ring soil sample according to the soil pore size distribution index n of the cutting ring soil sample, the saturated hydraulic conductivity K.sub.s of the cutting ring soil sample, the saturated soil moisture content θ.sub.s, the initial soil moisture content θ.sub.i, and the residual soil moisture content θ.sub.r of the cutting ring soil sample, and the shape coefficient a of the soil moisture profile of the cutting ring soil sample, and proceeding to step F; and step F: obtaining a sorptivity s of the cutting ring soil sample according to the effective saturated hydraulic conductivity K.sub.e of the cutting ring soil sample, and the change of the cumulative infiltration amount I of the cutting ring soil sample over the infiltration time t.

3. The method for predicting soil hydraulic properties based on an upward infiltration experiment of a cutting ring soil sample according to claim 2, wherein in step E, by substituting the soil pore size distribution index n of the cutting ring soil sample, the saturated hydraulic conductivity K.sub.s of the cutting ring soil sample, the saturated soil moisture content θ.sub.s, the initial soil moisture content θ.sub.i, and the residual soil moisture content θ.sub.r of the cutting ring soil sample and the shape coefficient a of the soil moisture profile of the cutting ring soil sample into the following formula: $K_e=\frac{\left(a+1\right)\left({\theta }_s-{\theta }_i\right)}{\left(\mathrm{am}/n+1\right)\left[{\theta }_s-{\theta }_i-a\left({\theta }_i-{\theta }_r\right)\right]}{K}_s$ the effective saturated hydraulic conductivity K.sub.e of the cutting ring soil sample is obtained, wherein m=3n+2.

4. The method for predicting soil hydraulic properties based on an upward infiltration experiment of a cutting ring soil sample according to claim 2, wherein in step F, by substituting the effective saturated hydraulic conductivity K.sub.e of the cutting ring soil sample, and the change of the cumulative infiltration amount I of the cutting ring soil sample over the infiltration time t into the following formula: $t=-\frac{I}{K_e}-\frac{s^2}{2K_e^2}\ln \left(1-\frac{2K_e}{s^2}I\right)$ the sorptivity s of the cutting ring soil sample is obtained through a fitting operation.

5. The method for predicting soil hydraulic properties based on an upward infiltration experiment of a cutting ring soil sample according to claim 2, wherein after step F is executed, the method further comprises the following step: step G: by substituting the sorptivity s of the cutting ring soil sample, the soil pore size distribution index n of the cutting ring soil sample, the saturated soil moisture content θ.sub.s of the cutting ring soil sample, the initial soil moisture content θ.sub.i of the cutting ring soil sample, the saturated hydraulic conductivity K.sub.s of the cutting ring soil sample and the constant lower boundary pressure h.sub.p into the following formula: $h_d=\{\begin{array}{cc}\frac{\left(3n+1\right)s^2}{2K_s\left({\theta }_s-{\theta }_i\right)}& h_p\le -h_d\\ \frac{m-1}{m}\left(\frac{s^2}{2K_s\left({\theta }_s-{\theta }_i\right)}-h_p\right)& h_p>-h_d\end{array}$ an air-entry suction value h.sub.d of the cutting ring soil sample is obtained, wherein m=3n+2.

6. The method for predicting soil hydraulic properties based on an upward infiltration experiment of a cutting ring soil sample according to claim 1, wherein in step C, by substituting the saturated soil moisture content θ.sub.s, the initial soil moisture content θ.sub.i and the residual soil moisture content θ.sub.r of the cutting ring soil sample, as well as the average soil moisture content increase U in the wetting region of the cutting ring soil sample into the following formula: $a=\frac{{\theta }_s-{\theta }_i-U}{{\theta }_i+U-{\theta }_r}$ the shape coefficient a of the soil moisture profile of the cutting ring soil sample is obtained.

7. The method for predicting soil hydraulic properties based on an upward infiltration experiment of a cutting ring soil sample according to claim 1, wherein in step D, according to the shape coefficient a of the soil moisture profile of the cutting ring soil sample, and in combination with the following formula: $a=\frac{n\left(2n+3\right)}{\left(2n+2\right)\left(2n+1\right)}$ the soil pore size distribution index n of the cutting ring soil sample is obtained using a fixed-point iteration method.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0027] FIG. 1 is a flowchart of a method for predicting soil hydraulic properties based on an upward infiltration experiment of a cutting ring soil sample according to the present disclosure; and

[0028] FIGS. 2A, 2B, and 2C show comparison between a cumulative infiltration amount curve simulated by an upward infiltration analytical model and a numerical simulation curve for three kinds of soil, where FIG. 2A shows a comparison diagram for the cumulative infiltration of a clay soil, FIG. 2B shows a comparison diagram for the cumulative infiltration of a silty soil, and FIG. 2C shows a comparison diagram for the cumulative infiltration of a sandy soil.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0029] The specific implementation of the present disclosure will be further described in detail below with reference to the drawings.

[0030] The present disclosure provides a method for predicting soil hydraulic properties based on an upward infiltration experiment of a cutting ring soil sample. Firstly, Richards Equation (1931) and boundary conditions thereof are analyzed as follows:

[00006]$\begin{array}{cc}\frac{\partial \theta }{\partial t}=\frac{\partial }{\partial z}\left[K\left(h\right)\frac{\partial h}{\partial z}+K\left(h\right)\right]& \left(1\right)\end{array}$$\theta \left(z,0\right)={\theta }_i$$\theta \left(0,t\right)={\theta }_s$$\theta \left(\infty ,t\right)={\theta }_i$

[0031] The Brooks-Corey model is analyzed to describe the relationships among unsaturated hydraulic conductivity, soil matrix potential and soil moisture content as follows:

[00007]$\begin{array}{cc}S=\{\begin{array}{cc}\frac{\theta -{\theta }_r}{{\theta }_s-{\theta }_r}={.Math.\frac{h_d}{h}.Math.}^n,& h<-h_d\\ 1,& h\ge -h_d\end{array}& \left(2\right)\end{array}$$\begin{array}{cc}K\left(h\right)=K_s{.Math.\frac{h_d}{h}.Math.}^m={K_s\left(\frac{\theta -{\theta }_r}{{\theta }_s-{\theta }_r}\right)}^{\frac{m}{n}}& \left(3\right)\end{array}$

[0032] Where h denotes a soil matrix potential (cm).

[0033] Through the principle of least action, flux-concentration assumption and mean value theorem of integrals, a profile equation describing soil water upward infiltration is obtained, where the flux-concentration assumption and mean value theorem of integrals are used to simplify the functional form of the profile equation, and the flux-concentration assumption is expressed as follows (Philip, 1973):

[00008]$\begin{array}{cc}F\left(\theta \right)=\frac{J_w+K_i}{J_{w0}+K_i}=\frac{S-S_i}{1-S_i}& \left(4\right)\end{array}$

[0034] Where J.sub.wo denotes a water flux density (cm min.sup.−1) at a water inlet, J.sub.w denotes a soil water flux density (cm min.sup.−1), and K.sub.i denotes an unsaturated hydraulic conductivity (cm min.sup.−1) corresponding to an initial soil moisture content; the meanings of other parameters are denoted above.

[0035] Based on the principle of mass conservation and profile equation, complete analytical models (6)-(8) describing one-dimensional water upward infiltration of homogeneous dry soil under a boundary condition of a constant pressure are derived.

[0036] The principle of mass conservation is used to obtain the functional relationship between the cumulative infiltration amount I and the infiltration time t, and specifically, the principle of mass conservation is physically expressed as follows:

[00009]$\begin{array}{cc}{J}_{w0}=\frac{\mathrm{dI}}{\mathrm{dt}}& \left(5\right)\end{array}$$\begin{array}{cc}S=\frac{\theta -{\theta }_r}{{\theta }_s-{\theta }_r}={\left(1-\left(1-S_i^{\frac{1}{a}}\right)\frac{z}{z_f}\right)}^a& \left(6\right)\end{array}$$\begin{array}{cc}t=-\frac{I}{K_e}-\frac{s^2}{2K_e^2}\ln \left(1-\frac{2K_e}{s^2}I\right)& \left(7\right)\end{array}$$\begin{array}{cc}I={\mathrm{Uz}}_f& \left(8\right)\end{array}$$\mathrm{Where},$$\begin{array}{cc}U=\frac{{\theta }_s-{\theta }_i-a\left({\theta }_i-{\theta }_r\right)}{a+1}& \left(9\right)\end{array}$$\begin{array}{cc}{s}^2=2K_e\left({\theta }_s-{\theta }_i\right)\left(h_p-H_f\right)& \left(10\right)\end{array}$$\begin{array}{cc}a=\frac{n\left(2n+3\right)}{\left(2n+2\right)\left(2n+1\right)}& \left(11\right)\end{array}$$\begin{array}{cc}{K}_e=\frac{\left(a+1\right)\left({\theta }_s-{\theta }_i\right)}{\left(\mathrm{am}/n+1\right)\left[{\theta }_s-{\theta }_{𝔦}-a\left({\theta }_i-{\theta }_r\right)\right]}{K}_s& \left(12\right)\end{array}$$\begin{array}{cc}{H}_f=\{\begin{array}{cc}-\frac{\left(\mathrm{am}/n+1\right)\left[{\theta }_s-{\theta }_i-a\left({\theta }_i-{\theta }_r\right)\right]}{\left(a+1\right)\left({\theta }_s-{\theta }_i\right)\left(m-1\right)}{h}_d& h_p\le -h_d\\ h_p-\frac{\left(\mathrm{am}/n+1\right)\left[{\theta }_s-{\theta }_i-a\left({\theta }_i-{\theta }_r\right)\right]}{\left(a+1\right)\left({\theta }_s-{\theta }_i\right)}\left(\frac{m}{m-1}{h}_d+h_p\right)& h_p>-h_d\end{array}& \left(13\right)\end{array}$

[0037] Based on the foregoing formulas (6)-(13), the present disclosure provides a method for predicting soil hydraulic properties based on an upward infiltration experiment of a cutting ring soil sample. In practical application, as shown in. FIG. 1, step A to step G are executed as follows to predict soil hydraulic properties of the cutting ring soil sample.

[0038] step A: During a one-dimensional water upward infiltration experiment of a homogeneous dry soil for the cutting ring soil sample under a boundary condition of a constant pressure h.sub.p, record the change of the cumulative infiltration amount I of the cutting ring soil sample over infiltration time t, and a infiltration time t.sub.e for the wetting front in soil to reach the upper surface of the cutting ring soil sample, measure a saturated hydraulic conductivity K.sub.s of the cutting ring soil sample, and proceed to step B.

[0039] step B: According to the change of the cumulative infiltration amount I of the cutting ring soil sample over the infiltration time t, obtain a cumulative infiltration amount I.sub.te, and an advancing distance z.sub.fte of the wetting front in soil at the time t.sub.e, as well as a ratio of z.sub.fte to I.sub.te, ie, an average soil moisture content increase U in the wetting region of the cutting ring soil sample, and proceed to step C.

[0040] step C: By substituting the saturated soil moisture content θ.sub.s, the initial soil moisture content θ.sub.i and the residual soil moisture content θ.sub.r of the cutting ring soil sample, as well as the average soil moisture content increase U in the wetting region of the cutting ring soil sample into the following formula:

[00010]$\begin{array}{cc}a=\frac{{\theta }_s-{\theta }_i-U}{{\theta }_i+U-{\theta }_r}& \left(14\right)\end{array}$

[0041] obtain the shape coefficient a of the soil moisture profile of the cutting ring soil sample, and proceed to step D.

[0042] step D: According to the shape coefficient a of the soil moisture profile of the cutting ring soil sample, and in combination with the following formula:

[00011]$\begin{array}{cc}a=\frac{n\left(2n+3\right)}{\left(2n+2\right)\left(2n+1\right)}& \left(15\right)\end{array}$

[0043] obtain the soil pore size distribution index n of the cutting ring soil sample using a fixed-point iteration method, and proceed to step E.

[0044] step E: By substituting the soil pore size distribution index n of the cutting ring soil sample, the saturated hydraulic conductivity K.sub.s of the cutting ring soil sample, the saturated soil moisture content 0,, the initial soil moisture content θ.sub.i, and the residual soil moisture content θ.sub.r of the cutting ring soil sample, and the shape coefficient a of the soil moisture profile of the cutting ring soil sample into the following formula:

[00012]$\begin{array}{cc}{K}_e=\frac{\left(a+1\right)\left({\theta }_s-{\theta }_i\right)}{\left(\mathrm{am}/n+1\right)\left[{\theta }_s-{\theta }_i-a\left({\theta }_i-{\theta }_r\right)\right]}{K}_s& \left(16\right)\end{array}$

[0045] the effective saturated hydraulic conductivity K.sub.e of the cutting ring soil sample is obtained (m=3n+2), and proceed to step F.

[0046] step F: By substituting the effective saturated hydraulic conductivity K.sub.e of the cutting ring soil sample, and the change of the cumulative infiltration amount I of the cutting ring soil sample over the infiltration time t into the following formula:

[00013]$\begin{array}{cc}t=-\frac{I}{K_e}-\frac{s^2}{2K_e^2}\ln \left(1-\frac{2K_e}{s^2}I\right)& \left(17\right)\end{array}$

[0047] the sorptivity s of the cutting ring soil sample is obtained through a fitting operation. Then proceed to step G.

[0048] step G: By substituting the sorptivity s of the cutting ring soil sample, the soil pore size distribution index n of the cutting ring soil sample, the saturated soil moisture content θ.sub.s of the cutting ring soil sample, the initial soil moisture content θ.sub.i of the cutting ring soil sample, the saturated hydraulic conductivity K.sub.s of the cutting ring soil sample and the constant lower boundary pressure h.sub.p into the following formula:

[00014]$\begin{array}{cc}{h}_d=\{\begin{array}{cc}\frac{\left(3n+1\right)s^2}{2K_s\left({\theta }_s-{\theta }_i\right)}& h_p\le -h_d\\ \frac{m-1}{m}\left(\frac{s^2}{2K_s\left({\theta }_s-{\theta }_i\right)}-h_p\right)& h_p>-h_d\end{array}& \left(18\right)\end{array}$

[0049] an air-entry suction value h.sub.d of the cutting ring soil sample is obtained, where m=3n+2.

[0050] In order to verify the accuracy of the upward infiltration analytical model, FIGS. 2A, 2B and 2C respectively show the difference between the cumulative infiltration amount curve simulated by analytical solution and the numerical solution of formula (1) for three typical soil textures (see Table 1 for soil hydraulic properties). It can be seen from the curve diagrams that the analytical solution results are almost consistent with the numerical simulation results, the relative errors between them are 2.17%, 1.12% and 1.56% respectively, and determination coefficients are 0.996, 0.999 and 0.998 respectively. Therefore, the upward infiltration analytical model can accurately simulate the upward infiltration for different soil textures, which provides a theoretical basis for obtaining corresponding hydraulic parameters.

[0051] When the present disclosure is applied to practice, the following three typical soil textures (see Table 1 for soil hydraulic properties) and soil sample test cases with different ring cutter heights are used to further illustrate the method for quickly estimating soil hydraulic properties based on the water upward infiltration experiment.

TABLE-US-00001 TABLE 1 θ.sub.r θ.sub.s θ.sub.i h.sub.d s Soil texture cm.sup.3 cm.sup.−3 cm n m cm min.sup.−0.5 Sandy soil 0.041 0.453 0.065 14.66 0.322 2.966 0.4759 Silty soil 0.015 0.501 0.02  20.75 0.211 2.633 0.3742 clay soil 0.09  0.475 0.121 37.31 0.131 2.393 0.1413 Note: the meanings of the parameters in this Table are denoted above.

[0052] (1) In case where a ring cutter has a height of 5 cm, and a soil sample is in a texture of clay soil, set a lower boundary pressure to −37.31 cm in a one-dimensional water upward infiltration experiment.

[0053] Based on the above calculation parameters, the calculation process of soil hydraulics properties is as follows.

[0054] Step A: During a one-dimensional water upward infiltration experiment for a clay soil sample under a boundary condition of a constant pressure, record data about the cumulative infiltration amount and infiltration time, a time t.sub.e=7800 s for the wetting front to reach the upper surface of the soil sample, and measure a saturated hydraulic conductivity K.sub.s=0.001 cm min.sup.−1 of the soil sample.

[0055] Steps B-D: According to the time for the wetting front to reach the upper surface of the soil sample, query the data about the cumulative infiltration amount and infiltration time to obtain a cumulative infiltration amount I.sub.te=1.536 cm and an advancing distance z.sub.fte=5 cm of the wetting front, obtain U=0.307 according to formula (8), obtain α=0.139 according to formula (14), and obtain the soil pore size distribution index (0.120) of the soil sample according to formula (15).

[0056] Steps E-G: Obtain K.sub.e=0.00031 cm min.sup.−1 according to formula (16), for a known value of K.sub.e, fit time series data of the cumulative infiltration amount according to formula (17) to obtain a sorptivity of 0.1371 cm min.sup.−0.5, and obtain an air-entry suction value of 36.98 cm according to formula (18).

[0057] (2) In case where a ring cutter has a height of 5 cm, and a soil sample is in a texture of silty soil, set a lower boundary pressure to −20.75 cm in a one-dimensional water upward infiltration experiment.

[0058] Based on the above calculation parameters, the calculation process of soil hydraulics properties is as follows.

[0059] Step A: During a one-dimensional water upward infiltration experiment for a silty soil sample under a boundary condition of a constant pressure, record data about the cumulative infiltration amount and infiltration time, a time t.sub.e=1913 s for the wetting front to reach the upper surface of the soil sample, and measure a saturated hydraulic conductivity K.sub.s=0.011 cm min.sup.−1 of the soil sample.

[0060] Steps B-D: According to the time for the wetting front to reach the upper surface of the soil sample, query the data about the cumulative infiltration amount and infiltration time to obtain a cumulative infiltration amount I.sub.te=1.981 cm and an advancing distance z.sub.fte=5 cm of the wetting front, obtain U=0.396 according to formula (8), obtain α=0.193 according to formula (14), and obtain the soil pore size distribution index (0.186) of the soil sample according to formula (15).

[0061] Steps E-G: Obtain K.sub.e=0.0037 cm min.sup.−1 according to formula (16), for a known value of K.sub.e, fit time series data of the cumulative infiltration amount according to formula (17) to obtain a sorptivity of 0.3593 cm min.sup.−0.5, and obtain an air-entry suction value of 19.58 cm according to formula (18).

[0062] (3) In case where a ring cutter has a height of 5 cm, and a soil sample is in a texture of sandy soil, set a lower boundary pressure to −14.66 cm in a one-dimensional water upward infiltration experiment.

[0063] Based on the above calculation parameters, the calculation process of soil hydraulics properties is as follows.

[0064] Step A: During a one-dimensional water upward infiltration experiment for a silty soil sample under a boundary condition of a constant pressure, record data about the cumulative infiltration amount and infiltration time, a time t.sub.e=694 s for the wetting front to reach the upper surface of the soil sample, and measure a saturated hydraulic conductivity Ks=0.043 cm min.sup.−1 of the soil sample.

[0065] Steps B-D: According to the time for the wetting front to reach the upper surface of the soil sample, query the data about the cumulative infiltration amount and infiltration time to obtain a cumulative infiltration amount I.sub.te=1.522 cm and an advancing distance z.sub.fte=5 cm of the wetting front, obtain U=0.304 according to formula (8), obtain α=0.255 according to formula (14), and obtain the soil pore size distribution index (0.290) of the soil sample according to formula (15).

[0066] Steps E-G: Obtain K.sub.e=0.0156 cm min.sup.−1 according to formula (16), for a known value of K.sub.e, fit time series data of the cumulative infiltration amount according to formula (17) to obtain a sorptivity of 0.4714 cm min.sup.−0.5, and obtain an air-entry suction value of 12.94 cm according to formula (18).

[0067] (4) In case where a ring cutter has a height of 10 cm, and a soil sample is in a texture of clay soil, set a lower boundary pressure to −37.31 cm in a one-dimensional water upward infiltration experiment.

[0068] Based on the above calculation parameters, the calculation process of soil hydraulics properties is as follows.

[0069] Step A: During a one-dimensional water infiltration experiment for a clay soil sample under a boundary condition of a constant pressure, record data about the cumulative infiltration amount and infiltration time, a time t.sub.e=539 min for the wetting front to reach the upper surface of the soil sample, and measure a saturated hydraulic conductivity K.sub.s=0.001 cm min.sup.−1 of the soil sample.

[0070] Steps B-D: According to the time for the wetting front to reach the upper surface of the soil sample, query the data about the cumulative infiltration amount and infiltration time to obtain a cumulative infiltration amount I.sub.e=3.059 cm and an advancing distance z.sub.fte=10 cm of the wetting front, obtain. U=0.306 according to formula (8), obtain α=0.143 according to formula (14), and obtain the soil pore size distribution index (0.123) of the soil sample according to formula (15).

[0071] Steps E-G: Obtain K.sub.e=0.00031 cm min.sup.−1 according to formula (16), for a known value of K.sub.e, fit time series data of the cumulative infiltration amount according to formula (17) to obtain a sorptivity of 0.1371 cm min.sup.−0.5, and obtain an air-entry suction value of 36.36 cm according to formula (18).

[0072] (5) In case where a ring cutter has a height of 10 cm, and a soil sample is in a texture of silty soil, set a lower boundary pressure to −20.75 cm in a one-dimensional water upward infiltration experiment.

[0073] Based on the above calculation parameters, the calculation process of soil hydraulics properties is as follows.

[0074] Step A: During a one-dimensional water upward infiltration experiment for a silty soil sample under a boundary condition of a constant pressure, record data about the cumulative infiltration amount and infiltration time, a time t.sub.e=138 min for the wetting front to reach the upper surface of the soil sample, and measure a saturated hydraulic conductivity Ks=0.011 cm min.sup.−1 of the soil sample.

[0075] Steps B-D: According to the time for the wetting front to reach the upper surface of the soil sample, query the data about the cumulative infiltration amount and infiltration time to obtain a cumulative infiltration amount I.sub.te=3.926 cm and an advancing distance z.sub.fte=10 cm of the wetting front, obtain U=0.393 according to formula (8), obtain α=0.204 according to formula (14), and obtain the soil pore size distribution index (0.203) of the soil sample according to formula (15).

[0076] Steps E-G: Obtain K.sub.e=0.0038 cm min.sup.−1 according to formula (16), for a known value of K.sub.e, fit time series data of the cumulative infiltration amount according to formula (17) to obtain a sorptivity of 0.3638 cm min.sup.−0.5, and obtain an air-entry suction value of 19.77 cm according to formula (18).

[0077] (6) in case where a ring cutter has a height of 10 cm, and a soil sample is in a texture of sandy soil set, a lower boundary pressure to −14.66 cm in a one-dimensional water upward infiltration experiment.

[0078] Based on the above calculation parameters, the calculation process of soil hydraulics properties is as follows.

[0079] Step A: During a one-dimensional water upward infiltration experiment for a silty soil sample under a boundary condition of a constant pressure, record data about the cumulative infiltration amount and infiltration time, a time t.sub.e=54 min for the wetting front to reach the upper surface of the soil sample, and measure a saturated hydraulic conductivity Ks=0.043 cm min.sup.−1 of the soil sample.

[0080] Steps B-D: According to the time for the wetting front to reach the upper surface of the soil sample, query the data about the cumulative infiltration amount and infiltration time to obtain a cumulative infiltration amount I.sub.te=2.992 cm and an advancing distance z.sub.fte=10 cm of the wetting front, obtain U=0.299 according to formula (8), obtain α=0.275 according to formula (14), and obtain the soil pore size distribution index (0.333) of the soil sample according to formula (15).

[0081] Steps E-G: Obtain K.sub.e=0.0161 cm min.sup.−1 according to formula (16), for a known value of K.sub.e, fit time series data of the cumulative infiltration amount according to formula (17) to obtain a sorptivity of 0.4758 cm min.sup.−0.5, and obtain an air-entry suction value of 13.51 cm according to formula (18).

[0082] From the above six specific calculation results (Table 2 shows the comparison between the soil hydraulic properties calculated based on the present disclosure and tested soil parameters), it can be seen that there is little difference between the soil hydraulic properties calculated by the present disclosure and the tested soil parameters. For the soil pore size distribution indexes obtained from the two methods, a maximum error between them is only 0.032, which occurs in the sandy soil with the ring cutter having a height of 5 cm; for the sorptivity obtained from the two methods, they are almost the same; for the air-entry suction values, the maximum error between them is only 1.72 cm, which occurs in the sandy soil with the ring cutter having a height of 5 cm. Therefore, based on the one-dimensional water upward infiltration experiment of homogeneous dry soil under a boundary condition of a constant pressure, the present disclosure can accurately predict the hydraulic properties of cutting ring soil sample, especially with fine soil texture.

TABLE-US-00002 TABLE 2 Soil n″ s s′ s″ h.sub.d h.sub.d′ h.sub.d″ texture n n′ % cm min.sup.−0.5 % cm % Sandy soil 0.322 0.290 0.333 0.4759 0.4714 0.4758 14.66 12.94 13.51 Silty soil 0.211 0.186 0.203 0.3742 0.3593 0.3628 20.75 19.58 19.77 Clay soil 0.131 0.120 0.123 0.1413 0.1371 0.1371 37.31 36.98 36.36 Note: the meanings of the parameters in this Table are denoted above, where n′, s′ and h.sub.d′ denote estimation parameters for a soil sample with a ring cutter height of 5 cm, and n″, s″ and h.sub.d″ denote estimation parameters for a soil sample with a ring cutter height of 10 cm.

[0083] According to the above technical solution, a method for predicting soil hydraulic properties based on an upward infiltration experiment of cutting ring soil sample is designed. According to the new design approach, based on an inverse process of one-dimensional water upward infiltration of homogeneous dry soil under a boundary condition of a constant pressure, by recording the data regarding a cumulative infiltration amount and a infiltration time during the upward infiltration process of the cutting ring soil sample and a time for the wetting front to reach the upper surface of the cutting ring soil sample, a saturated hydraulic conductivity of the cutting ring soil sample is measured, and the soil hydraulic properties are quickly obtained. The present disclosure overcomes the defects that the methods in the prior art for measuring soil hydraulic properties are time consuming, low in efficiency, and not suitable for rapid acquisition of large area in the field, avoids the problems of multiple solutions and non-convergence of parameters, and improves the efficiency of acquiring soil hydraulic properties, thereby providing an effective means for the investigation of soil hydraulic properties in a large area in the field.

[0084] Although the embodiments of the present disclosure are described in detail above in conjunction with the drawings, the present disclosure is not limited to the above-described embodiments, and various changes may be made without departing from the spirit of the present disclosure within the knowledge of those skilled in the art.