Method and device for determining water infiltration depth of double-layer soil and storage medium

Abstract

The invention discloses a method and device for determining the water infiltration depth of double-layer soil and a storage medium. The method comprises the steps that the water infiltration duration is collected; in response to the situation that the collected water infiltration duration is smaller than or equal to the predetermined duration when water infiltrates into the interface of the upper layer soil and the lower layer soil, it is indicated that water does not infiltrate into the lower layer soil, and the water infiltration depth of the double-layer soil is determined according to the relational expression of the water infiltration duration in the upper layer soil and the water infiltration depth; in response to the situation that the collected water infiltration duration is greater than the predetermined duration of water infiltration into the interface of the upper and lower layers of soil, indicating that the water has infiltrated into the lower layer of soil, determining a relational expression between the water infiltration duration and the water infiltration depth in the lower layer of soil according to the infiltration pressure at the interface of the upper and lower layers of soil, and determining the water infiltration depth of the double-layer soil in combination with a relational expression between the water infiltration duration and the water infiltration depth in the upper-layer soil. The infiltration depth determined by the method is more accurate and better fits engineering practice.

Description

Technical Field

[0001] The invention relates to the technical field of hydrology and geotechnical engineering, and in particular to a method, a device and a storage medium for determining the water infiltration depth of a double-layer soil. Background Art

[0002] Determining the depth of water infiltration in the soil is of great engineering significance. Above the current water infiltration depth, the soil is saturated; below the current water infiltration depth, the soil is unsaturated; there are large differences between saturated soil and unsaturated soil in terms of deadweight, water content, shear strength, etc. In irrigation projects, the depth of water infiltration in the soil is closely related to plant growth. In geotechnical engineering, when a retaining wall or foundation pit encounters rainfall, the depth of rainwater infiltration in the soil is related to the stability of the retaining wall or foundation pit. Similar problems also occur in projects such as embankments, cuttings, and slopes.

[0003] Currently known methods for calculating soil moisture infiltration depth all assume that the soil is a homogeneous single-layer soil. Although this assumption simplifies the problem, in actual engineering, most soils are layered. When water infiltrates into different soils, its travel speed will change. Therefore, the existing method of calculating soil moisture infiltration depth assuming that the soil is a homogeneous single-layer soil is bound to be inconsistent with actual engineering practices, resulting in technical problems such as large differences between the calculated results and the actual situation. Summary of the invention

[0004] The purpose of the present invention is to overcome the deficiencies in the prior art and to provide a method, device and storage medium for determining the water infiltration depth of a double-layer soil, so that the determined infiltration depth is more accurate and more in line with engineering practice.

[0005] To achieve the above object, the present invention is implemented by adopting the following technical solutions:

[0006] In a first aspect, the present invention provides a method for determining the water infiltration depth of a double-layer soil, the method comprising the following steps:

[0007] Collect water infiltration time;

[0008] In response to the collected water infiltration time being less than or equal to the predetermined time for water to infiltrate to the interface between the upper and lower soil layers, indicating that water has not infiltrated into the lower soil layer, the water infiltration depth of the double-layer soil is determined according to the relationship between the water infiltration time in the upper soil layer and the water infiltration depth;

[0009] In response to the collected water infiltration time being greater than the predetermined time for water to infiltrate to the interface between the upper and lower soil layers, indicating that water has infiltrated into the lower soil layer, the relationship between the water infiltration time and the water infiltration depth in the lower soil layer is determined according to the infiltration pressure at the interface between the upper and lower soil layers, and the water infiltration depth of the double-layer soil is determined in combination with the relationship between the water infiltration time and the water infiltration depth in the upper soil layer.

[0010] In combination with the first aspect, further, a method for determining the relationship between the water infiltration time and the water infiltration depth in the upper soil layer includes:

[0011] The relationship between water infiltration and water infiltration depth is determined according to the saturated water content and initial water content of the upper soil:

[0012] (a)

[0013] Where I is the infiltration of water into the upper soil layer, z is the infiltration depth of water into the upper soil layer, θ s1 ,θ i1 are the saturated water content and initial water content of the upper soil layer, respectively;

[0014] Based on formula (a), the upper soil water infiltration amount I is derived from the water infiltration time t to determine the upper soil surface water infiltration rate i:

[0015] (b)

[0016] According to the surface water thickness H of the upper soil layer and the saturated permeability coefficient K of the upper soil layer, s1 and initial water content θ i1 The corresponding matrix suction s f1 Reconstruct the expression of the water infiltration rate i of the upper soil surface:

[0017] (c)

[0018] Since the water infiltration rate of the same layer of homogeneous soil is equal, combining equation (b) and equation (c), we get:

[0019] (d)

[0020] By integrating formula (d) over time, we can obtain the relationship between the water infiltration time t and the water infiltration depth z of the upper soil layer:

[0021] (e).

[0022] In combination with the first aspect, further, the method for determining the time taken for the water to infiltrate the interface between the upper and lower soil layers includes:

[0023] Substituting the upper soil thickness z1 as the upper soil water infiltration depth z into formula (e), the time t1 for water to infiltrate to the interface between the upper and lower soil layers is calculated as:

[0024] (j).

[0025] In combination with the first aspect, further, the relationship between the water infiltration time and the water infiltration depth in the lower soil layer determined according to the infiltration pressure at the interface between the upper and lower soil layers includes:

[0026] If the infiltration water volume at the interface between the upper and lower soil layers is less than or equal to the infiltrative water volume of the lower soil layer, the infiltration pressure at the interface between the upper and lower soil layers is determined to be less than or equal to zero. The infiltration rate i2 of the lower soil surface is calculated based on the relevant parameters of the upper soil layer:

[0027] (k)

[0028] Where H is the thickness of water on the upper soil surface; K s1 、s f1 are the saturated permeability coefficient of the upper soil and the matrix suction corresponding to the initial water content; z1 is the thickness of the upper soil;

[0029] According to the saturated water content of the lower soil θ s2 , initial water contentθ i2 , and the water infiltration time t2 in the lower soil, the expression of the infiltration rate i2 on the lower soil surface is reconstructed:

[0030] (l)

[0031] Where z2 is the water infiltration depth in the lower soil layer;

[0032] Since the water infiltration rate of the same layer of homogeneous soil is equal, the combination of equation (k) and equation (l) yields:

[0033] (m)

[0034] After integrating the formula (m) over time, the relationship between the water infiltration time t2 and the water infiltration depth z2 in the lower soil layer is obtained as follows:

[0035] (n).

[0036] In combination with the first aspect, further, the relationship between the water infiltration time and the water infiltration depth in the lower soil layer determined according to the infiltration pressure at the interface between the upper and lower soil layers includes:

[0037] If the infiltration water volume at the interface between the upper and lower soil layers is greater than the infiltrative water volume of the lower soil layer, it is determined that the infiltration pressure at the interface between the upper and lower soil layers is greater than zero. The surface infiltration rate i2 of the lower soil layer is calculated based on the predetermined infiltration pressure at the interface between the upper and lower soil layers and the relevant parameters of the lower soil layer. The expression is:

[0038] (p)

[0039] Where H' is the infiltration pressure at the interface between the upper and lower soil layers; z2 is the water infiltration depth in the lower soil layer; K s2 、s f2 are the saturated permeability coefficient of the underlying soil and the matrix suction corresponding to the initial water content;

[0040] According to the saturated water content of the lower soil θ s2 , initial water contentθ i2 , and the water infiltration time t2 in the lower soil, the expression of the infiltration rate i2 on the lower soil surface is reconstructed as:

[0041] (q)

[0042] Since the water infiltration rate of the same layer of homogeneous soil is equal, the combination of equation (p) and equation (q) yields:

[0043] (r)

[0044] By integrating both sides of formula (r) over time, we can obtain the relationship between the water infiltration time t2 and the water infiltration depth z2 in the lower soil layer:

[0045] (s).

[0046] In combination with the first aspect, further, the method for determining the infiltration pressure at the interface between the upper and lower soil layers includes:

[0047] Assuming that the water infiltration depth in the lower soil layer is equal to the thickness z1 of the upper soil layer, the positions of the upper and lower soil layers are interchanged, and all other parameters remain unchanged. Based on the equal energy lost in the infiltration of the upper and lower soil layers, the time required to reach the infiltration depth is also equal. Therefore, the equivalent equation is established as follows:

[0048] (w)

[0049] Where H´ is the infiltration pressure at the interface between the upper and lower soil layers; θ s1 ,θ i1 are the saturated water content and initial water content of the upper soil layer, respectively; K s1 、s f1are the saturated permeability coefficient of the upper soil and the matrix suction corresponding to the initial water content; θ s2 ,θ i2 are the saturated water content and initial water content of the lower soil layer, respectively; K s2 、s f2 are the saturated permeability coefficient of the lower soil and the matrix suction corresponding to the initial water content; H is the thickness of the surface water in the upper soil layer;

[0050] The infiltration pressure H´ at the interface between the upper and lower soil layers is calculated using formula (w).

[0051] In a second aspect, the present invention provides a device for determining the water infiltration depth of a double-layer soil, the device comprising:

[0052] Collection module: used to collect water infiltration duration;

[0053] The first determination module: in response to the collected water infiltration time being less than or equal to the predetermined time for water to infiltrate to the interface between the upper and lower soil layers, indicating that water has not infiltrated into the lower soil layer, the water infiltration depth of the double-layer soil is determined according to the relationship between the water infiltration time in the upper soil layer and the water infiltration depth;

[0054] The second determination module: in response to the collected water infiltration time being greater than the predetermined time for water to infiltrate to the interface between the upper and lower soil layers, indicating that the water has infiltrated into the lower soil layer, the relationship between the water infiltration time and the water infiltration depth in the lower soil layer is determined according to the infiltration pressure at the interface between the upper and lower soil layers, and the water infiltration depth of the double-layer soil is determined in combination with the relationship between the water infiltration time and the water infiltration depth in the upper soil layer.

[0055] In a third aspect, the present invention provides a computer device, comprising a processor and a memory connected to the processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the steps of the method described in any one of the first aspects are performed.

[0056] In a fourth aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, characterized in that when the program is executed by a processor, the steps of any one of the methods described in the first aspect are implemented.

[0057] In a fifth aspect, the present invention provides a computer program product, comprising a computer program / instruction, which, when executed by a processor, implements the steps of any one of the methods described in the first aspect.

[0058] Compared with the prior art, the present invention has the following beneficial effects:

[0059] The technical solution provided in the present application determines whether water has infiltrated into the lower soil layer according to the water infiltration time. If the water has not infiltrated into the lower soil layer, the water infiltration depth is determined according to the water infiltration situation of a single homogeneous soil. If the water has infiltrated into the lower soil layer, the water infiltration depth is determined according to the water infiltration situation of a double-layer soil. For the water infiltration situation of the double-layer soil, the infiltration pressure at the interface between the upper and lower soil layers is further considered, and the influencing factor of the different water infiltration speeds in different soils is incorporated into the water infiltration depth calculation, so that the finally determined water infiltration depth is more accurate and more in line with engineering practice. BRIEF DESCRIPTION OF THE DRAWINGS

[0060] In order to more clearly illustrate the technical solutions in the embodiments of the present invention, the drawings required for use in the description of the embodiments will be briefly introduced below. Obviously, the drawings in the following description are only embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on the provided drawings without paying creative work.

[0061] Figure 1 It is a flow chart of a method for determining water infiltration depth of double-layer soil provided in the first embodiment of the present invention;

[0062] Figure 2 is a schematic diagram of double-layer soil water infiltration provided by an embodiment of the present invention;

[0063] Figure 3 is a schematic diagram of calculating infiltration pressure provided by an embodiment of the present invention;

[0064] Figure 4 It is a block diagram of the internal structure of a computer device provided in Embodiment 3 of the present invention.

[0065] Description of reference numerals:

[0066] 1. Upper soil layer; 2. Lower soil layer; 3. Surface water. DETAILED DESCRIPTION

[0067] The technical solution of the present invention is described in detail below through the accompanying drawings and specific embodiments. It should be understood that the embodiments of the present application and the specific features in the embodiments are detailed descriptions of the technical solution of the present application, rather than limitations on the technical solution of the present application. In the absence of conflict, the embodiments of the present application and the technical features in the embodiments can be combined with each other.

[0068] Embodiment 1:

[0069] This embodiment provides a method for determining the water infiltration depth of a double-layer soil. Figure 1 The method of this implementation specifically includes the following steps:

[0070] Step S01: Collecting water infiltration time;

[0071] like Figure 2 As shown, it is a schematic diagram of water infiltration in a double-layer soil. Considering that the double-layer soil is not a single homogeneous soil, the water infiltration speed is different in different soils. Therefore, for the double-layer soil, it should first be determined whether the water has infiltrated from the upper soil 1 to the lower soil 2, and then the water infiltration depth at the current moment is determined for the two situations of water not infiltrating into the lower soil 2 and water infiltrating into the lower soil 2. Compared with the prior art that assumes that the soil is a homogeneous single-layer soil, this embodiment takes into account the characteristics of soil stratification in engineering practice, so that the final calculation result is more in line with engineering practice.

[0072] In an embodiment of the present invention, the time it takes for water to infiltrate to the interface between the upper and lower soil layers is determined in advance based on the thickness of the upper soil layer and related parameters, and then the collected water infiltration time is compared with the determined time it takes for water to infiltrate to the interface between the upper and lower soil layers, and then it is determined whether the water has infiltrated into the lower soil layer 2 based on the comparison result.

[0073] Step S02: Determine whether the water has infiltrated into the lower soil 2 according to the water infiltration time, and determine the water infiltration depth for the two situations of water not infiltrating into the lower soil 2 and water infiltrating into the lower soil 2, as follows:

[0074] Case 1: In response to the collected water infiltration time being less than or equal to the predetermined time for water to infiltrate to the interface between the upper and lower soil layers, indicating that water has not infiltrated into the lower soil layer, the water infiltration depth of the double-layer soil is determined according to the relationship between the water infiltration time in the upper soil layer and the water infiltration depth;

[0075] Scenario 2: In response to the collected water infiltration time being greater than the predetermined time for water to infiltrate to the interface between the upper and lower soil layers, indicating that water has infiltrated into the lower soil 2, the relationship between the water infiltration time and the water infiltration depth in the lower soil is determined according to the infiltration pressure at the interface between the upper and lower soil layers, and the relationship between the water infiltration time and the water infiltration depth in the upper soil 1 is combined to determine the water infiltration depth of the double-layer soil.

[0076] In different soils, the infiltration speed of water is different. In this embodiment, the water infiltration depth is determined according to the double-layer soil water infiltration situation when the water has infiltrated into the lower soil 2. For the double-layer soil water infiltration situation, the infiltration pressure at the interface between the upper and lower soil layers is further considered, and the influencing factor of the different infiltration speeds of water in different soils is incorporated into the water infiltration depth calculation, so that the finally determined water infiltration depth is more accurate and more in line with the actual project.

[0077] For the aforementioned situation 1, the method for determining the relationship between the water infiltration time and the water infiltration depth in the upper soil 1 specifically includes the following steps:

[0078] Step S02-01: Determine the relationship between water infiltration and water infiltration depth according to the saturated water content and initial water content of the upper soil 1:

[0079] (a)

[0080] Where I is the water infiltration of the upper soil layer, z is the water infiltration depth of the upper soil layer, θ s1 ,θ i1 are the saturated water content and initial water content of the upper soil 1, respectively;

[0081] Step S02-02: Based on formula (a), the water infiltration amount is differentiated with respect to the water infiltration time t to determine the water infiltration rate i on the upper soil surface:

[0082] (b)

[0083] Step S02-03: Based on the thickness H of the surface water 3 of the upper soil layer and the saturated permeability coefficient K of the upper soil layer 1 s1 and initial water content θ i1 The corresponding matrix suction s f1 Reconstruct the expression of water infiltration rate i on the surface of upper soil 1:

[0084] (c)

[0085] Step S02-04: The water infiltration rate of the same layer of homogeneous soil is equal. Therefore, by combining equations (b) and (c), we get:

[0086] (d)

[0087] Step S02-05: Integrate the time by equation (d) to obtain the relationship between the water infiltration time t and the water infiltration depth z of the upper soil layer 1:

[0088] (e).

[0089] Formula (e) is the relationship between the water infiltration depth z and the water infiltration time t when the water has not infiltrated into the lower soil 2. In the actual calculation process, the collected water infiltration time can be substituted into formula (e) to obtain the water infiltration depth in the upper soil 1, which is the final double-layer soil water infiltration depth. As an equivalent replacement, the function graph of the water infiltration depth z and the water infiltration time t can also be plotted according to formula (e), and the water infiltration depth at any time can be determined according to the function graph.

[0090] As an embodiment of the present invention, the time length for the water to infiltrate to the interface between the upper and lower soil layers described in step S01 can also be determined by formula (e). Specifically, the thickness z1 of the upper soil layer 1 can be substituted into formula (e) as the water infiltration depth z of the upper soil layer 1 to calculate the time length t1 for the water to infiltrate to the interface between the upper and lower soil layers:

[0091] (j)

[0092] Combining formula (j) we can determine:

[0093] When the infiltration time t≤t1, the water does not enter the lower soil 2, which belongs to the water infiltration problem of single-layer homogeneous soil;

[0094] When the infiltration time t>t1, water enters the lower soil layer 2, which belongs to the double-layer soil water infiltration problem.

[0095] For the aforementioned situation 2, when water infiltrates into the lower soil 2, its infiltration time in the upper soil 1 can be determined by formula (j). However, for the infiltration depth of water in the lower soil 2, considering the different characteristics of the upper and lower soil layers, the water infiltration depth may be different for the same water infiltration time. For this reason, according to the actual engineering situation, this embodiment divides the double-layer soil into two situations: the upper layer has strong water permeability and the lower layer has weak water permeability (i.e., "strong upper and weak lower type"), and the upper layer has weak water permeability and the lower layer has strong water permeability (i.e., "weak upper and strong lower type"), and determines the water infiltration depth respectively.

[0096] Continue reading Figure 2 , when the amount of infiltrated water at the interface between the upper and lower soil layers is less than or equal to the infiltrative water of the lower soil layer, the double-layer soil is of the "upper weak and lower strong type". When water infiltrates into the lower soil layer, its infiltration time in the upper soil layer can be determined by formula (j). For the "upper weak and lower strong" double-layer soil, the water infiltration of the lower soil layer is pressureless infiltration. At this time, the infiltration pressure at the interface between the upper and lower soil layers is less than or equal to zero, that is, the water infiltrating at the interface between the upper and lower soil layers depends on the infiltration rate of the upper soil surface. Therefore, the infiltration rate of the lower soil layer at the interface is equal to the infiltration rate of the soil surface. Therefore, let z=z1 in formula (c) to calculate the infiltration rate i2 of the lower soil surface, and its expression is:

[0097] (k)

[0098] From formula (b), we can know that the saturated water content of the lower soil θ s2 , initial water contentθ i2 , and the water infiltration time t2 in the lower soil, the expression of the infiltration rate i2 on the lower soil surface is reconstructed:

[0099] (l)

[0100] Where, z2 is the water infiltration depth in the lower soil 2; t2 and z2 are variables and are unknown parameters;

[0101] The water infiltration rate of the same layer of homogeneous soil is equal, so combining equations (k) and (l), we get

[0102] (m)

[0103] The right side of the equal sign of formula (m) is a constant. After integrating formula (m) over time, the relationship between the water infiltration time t2 and the water infiltration depth z2 in the lower soil layer is obtained:

[0104] (n)

[0105] There exists t=t1+t2, so adding equation (j) and equation (n) yields:

[0106] (o)

[0107] In formula (o), except for t and z2, which are to be determined, the other parameters are known. The water infiltration duration can be obtained through collection and determined by user input. Substituting it into formula (o) can determine the lower soil infiltration depth z2 at the corresponding moment. It is also possible to draw a t~z2 graph based on formula (o), and the corresponding lower soil infiltration depth z2 at any time can be obtained. Add z2 to the upper soil thickness z1 to obtain the total water infiltration depth z, which is the double-layer soil water infiltration depth that needs to be determined in the end.

[0108] On the contrary, Figure 2 As shown in the figure, when the amount of infiltrated water at the interface between the upper and lower soil layers is greater than the amount of infiltrated water in the lower soil layer 2, it is considered that the permeability of the upper soil layer 1 is greater than that of the lower soil layer 2, and the double-layer soil is of the "strong upper and weak lower type". Since the amount of infiltrated water is greater than the amount of infiltrated water, the water infiltration of the lower soil layer is pressurized infiltration. At this time, the infiltration pressure at the interface between the upper and lower soil layers is greater than zero, that is, the water infiltrated at the interface between the upper and lower layers depends on the infiltration pressure of the contact surface between the upper and lower soil layers. Assuming that the infiltration pressure at the interface between the upper and lower soil layers is H', the expression of the infiltration rate i2 of the lower soil surface can be obtained:

[0109] (p)

[0110] In the formula, K s2 、s f2 are the saturated permeability coefficient of the underlying soil and the matrix suction corresponding to the initial water content, respectively.

[0111] It should be noted that the infiltration pressure H' at the interface between the upper and lower soil layers is different from the surface water thickness H. H represents the pressure on the soil surface, and H' is the infiltration pressure at the contact surface between the upper and lower soil layers. The two are related but not the same. H' is unknown and is an undetermined parameter.

[0112] According to formula (b), the upper soil parameters are replaced by the lower soil parameters, and the saturated water content θ of the lower soil 2 is s2 , initial water contentθ i2 , and the water infiltration time t2 in the lower soil, the expression of the infiltration rate i2 on the lower soil surface is reconstructed as:

[0113] (q)

[0114] The water infiltration rate of the same layer of homogeneous soil is equal, so combining equations (p) and (q), we get

[0115] (r)

[0116] By integrating both sides of formula (r) over time, we can obtain the relationship between the water infiltration time t2 and the water infiltration depth z2 in the lower soil:

[0117] (s)

[0118] Adding equation (j) and equation (s), we get:

[0119] (t)

[0120] In formula (t), except for t, H' and z2, the other parameters are known. If we want to obtain the t~z2 image, we also need to solve H'. Figure 3 , and explain the process of solving H´ in detail.

[0121] like Figure 3 As shown in the figure, the principle diagram for solving the infiltration pressure H' is given. Figure 3 The left side (a) is a schematic diagram before the upper soil layer and the lower soil layer are exchanged, and the right side (b) is a schematic diagram after the upper soil layer and the lower soil layer are exchanged. First, assume that the infiltration depth of water in the upper soil layer 1 and the lower soil layer 2 is z1. By comparing Figure 3 From (a) and (b), we can see that except for the position swap of the upper soil 1 and the lower soil 2, all other parameters remain unchanged. Therefore, the energy lost by water during infiltration is equal, that is, the time required to reach the infiltration depth is equal. Based on this, an equivalent equation is established to determine H´.

[0122] Figure 3 The double-layer soil shown in (a) is of the "strong upper and weak lower" type. According to formula (t), and taking z2=z1 in formula (t), we can get

[0123] (u)

[0124] Figure 3 In (b), the upper soil 1 and the lower soil 2 are interchanged, and the soil is converted to the "upper weak and lower strong type". According to formula (o), z2=z1 is taken in formula (o). Because the upper and lower soil layers are interchanged, the formula (o) should be θ s1 and θ s2 Interchange, θ i1 and θ i2 Swap, K s1 and K s2 Swap, s f1 Replace with s f2 , so we can get

[0125] (v)

[0126] The time t in equations (u) and (v) is equal, that is,

[0127] (w)

[0128] H' can be calculated using formula (w). Substituting the calculated H' into formula (t) can obtain the "strong on top and weak on bottom" double-layer soil t~z2 image, that is, the corresponding lower soil infiltration depth z2 at any time can be obtained. Adding z2 to the upper soil thickness z1 is the total infiltration depth z.

[0129] In summary, the method provided in the embodiment of the present invention can determine t1 according to the calculation parameters by formula (j):

[0130] When the infiltration time t≤t1, the water does not infiltrate into the lower soil layer, and the infiltration depth z is determined by formula (e).

[0131] When the infiltration time t>t1, water enters the lower soil layer, and the infiltration pressure H´ at the contact surface of the upper and lower soil layers is determined by formula (w);

[0132] When H´>0, it is a "strong on top and weak on bottom" type. The formula (t) is used to determine the water infiltration depth z2 of the lower soil layer. The total infiltration depth z is obtained by adding z2 to the thickness of the upper soil layer z1.

[0133] When H´≤0, it is a "weak on top and strong on the bottom" type. Use formula (o) to determine the water infiltration depth z2 of the lower soil layer. Add z2 to the thickness of the upper soil layer z1 to get the total infiltration depth z.

[0134] In combination with specific experimental data, the water infiltration depth is determined using the method provided by the embodiment of the present invention:

[0135] A double-layer soil layer, the upper layer is silt, thickness z1=2m, parameter θs1 =0.40,θ i1 =0.20, K s1 =0.04m / h,s f1 =1.4m; the lower layer is clay soil, the parameter is θ s2 =0.45,θ i2 =0.18, K s2 =0.018m / h,s f2 =2.2m; depth of water accumulation on the soil surface H=0.05m.

[0136] According to the above parameters, it can be known from formula (j) that t1=3.72h.

[0137] When the infiltration time t≤3.72h, the water does not enter the lower soil 2, and the relationship between the infiltration depth z and the time t can be calculated by formula (e):

[0138] (x)

[0139] When the infiltration time t>3.72h, the water has entered the lower soil layer. The infiltration pressure H´ at the contact surface of the upper and lower soil layers is determined by formula (w) to be -1.43m, which is less than 0. Therefore, it is pressureless infiltration, which is a "weak upper and strong lower type".

[0140] The relationship between the depth of water infiltration z2 in the lower soil and time t is determined by formula (o):

[0141] (y)

[0142] In formula (y), , so formula (y) can also be written as

[0143] (z)

[0144] According to formula (x) and (z), the water infiltration depth at any time can be calculated. For example, when t=1h, it is less than 3.72h, and the infiltration depth z is 0.91m when t=5h, it is greater than 3.72h, and the infiltration depth z is 2.32m when t=5h.

[0145] It should be noted that Figure 1 Only the logical sequence of the method described in this embodiment is shown. In other possible embodiments of the present invention, different logical sequences may be used without conflict. Figure 1 The steps shown or described are completed in the order shown. The method provided in this embodiment can be applied to a terminal and can be performed by a double-layer soil water infiltration depth determination device, which can be implemented by software and / or hardware, and the device can be integrated in a terminal, for example: any smart phone, tablet computer or computer device with communication function.

[0146] Embodiment 2:

[0147] An embodiment of the present invention provides a device for determining water infiltration depth of a double-layer soil, the device comprising:

[0148] Collection module: used to collect water infiltration duration;

[0149] The first determination module: in response to the collected water infiltration time being less than or equal to the predetermined time for water to infiltrate to the interface between the upper and lower soil layers, indicating that water has not infiltrated into the lower soil layer, the water infiltration depth of the double-layer soil is determined according to the relationship between the water infiltration time in the upper soil layer and the water infiltration depth;

[0150] The second determination module: in response to the collected water infiltration time being greater than the predetermined time for water to infiltrate to the interface between the upper and lower soil layers, indicating that the water has infiltrated into the lower soil layer, the relationship between the water infiltration time and the water infiltration depth in the lower soil layer is determined according to the infiltration pressure at the interface between the upper and lower soil layers, and the water infiltration depth of the double-layer soil is determined in combination with the relationship between the water infiltration time and the water infiltration depth in the upper soil layer.

[0151] The device provided in the embodiment of the present invention can execute the method provided in the first embodiment, and has the corresponding functional modules and beneficial effects of the execution method.

[0152] Embodiment three:

[0153] The embodiment of the present invention further provides a computer device, which may be a server, and its internal structure diagram may be as shown in FIG. Figure 4 As shown. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. Among them, the processor, the memory and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store data acquired and generated in the method of the robot autonomously entering the packaging container. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, the method described in the aforementioned embodiment 1 is implemented.

[0154] Those skilled in the art will understand that Figure 4The structure shown in the figure is only a block diagram of a part of the structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine certain components, or have a different arrangement of components.

[0155] The computer device provided in the embodiment of the present invention can execute the method provided in any embodiment of the present invention, and has the corresponding functional modules and beneficial effects of the execution method.

[0156] Embodiment 4:

[0157] An embodiment of the present invention provides a computer-readable storage medium on which a computer program is stored. When the program is executed by a processor, the steps of the method described in the first embodiment are implemented.

[0158] Embodiment five:

[0159] An embodiment of the present invention provides a computer program product, including a computer program / instruction, which implements the steps of the method described in the first embodiment when executed by a processor.

[0160] Those skilled in the art will appreciate that the embodiments of the present application may be provided as methods, systems, or computer program products. Therefore, the present application may adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware. Moreover, the present application may adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program codes.

[0161] The present application is described with reference to the flowcharts and / or block diagrams of the methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or box in the flowchart and / or block diagram, as well as the combination of the processes and / or boxes in the flowchart and / or block diagram, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing device to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing device generate instructions for implementing the processes in the flowchart and / or block diagram. Figure 1 A process or multiple processes and / or boxes Figure 1 A device that provides the functions specified in a block or multiple blocks.

[0162] These computer program instructions may also be stored in a computer-readable memory capable of directing a computer or other programmable data processing device to operate in a specific manner, so that the instructions stored in the computer-readable memory produce an article of manufacture comprising an instruction device, which implements the process Figure 1 A process or multiple processes and / or boxes Figure 1 A function specified in one or more boxes.

[0163] These computer program instructions can also be loaded onto a computer or other programmable data processing device so that a series of operating steps are executed on the computer or other programmable device to produce a computer-implemented process, thereby providing instructions for implementing the process. Figure 1 A process or multiple processes and / or boxes Figure 1 The steps for the functions specified in one or more boxes.

[0164] The above is only a preferred embodiment of the present invention. It should be pointed out that for ordinary technicians in this technical field, several improvements and modifications can be made without departing from the technical principles of the present invention. These improvements and modifications should also be regarded as the scope of protection of the present invention.

Claims

1. A method for determining the water infiltration depth of a double-layer soil, characterized in that: The method comprises the following steps: Collect water infiltration time; In response to the collected water infiltration time being less than or equal to the predetermined time for water to infiltrate to the interface between the upper and lower soil layers, indicating that water has not infiltrated into the lower soil layer, the water infiltration depth of the double-layer soil is determined according to the relationship between the water infiltration time in the upper soil layer and the water infiltration depth; In response to the collected water infiltration time being greater than the predetermined time for water to infiltrate to the interface between the upper and lower soil layers, indicating that water has infiltrated into the lower soil layer, the relationship between the water infiltration time and the water infiltration depth in the lower soil layer is determined according to the infiltration pressure at the interface between the upper and lower soil layers, and the water infiltration depth of the double-layer soil is determined in combination with the relationship between the water infiltration time and the water infiltration depth in the upper soil layer.

2. The method for determining the double-layer soil water infiltration depth according to claim 1, characterized in that: The method for determining the relationship between the water infiltration time and the water infiltration depth in the upper soil layer includes: The relationship between water infiltration and water infiltration depth is determined according to the saturated water content and initial water content of the upper soil: （a） Where I is the infiltration of water into the upper soil layer, z is the infiltration depth of water into the upper soil layer, θ s1 ,θ i1 are the saturated water content and initial water content of the upper soil layer, respectively; Based on formula (a), the upper soil water infiltration amount I is derived from the water infiltration time t to determine the upper soil surface water infiltration rate i: （b） According to the surface water thickness H of the upper soil layer and the saturated permeability coefficient K of the upper soil layer, s1 and initial water content θ i1 The corresponding matrix suction s f1 Reconstruct the expression of the water infiltration rate i of the upper soil surface: （c） Since the water infiltration rate of the same layer of homogeneous soil is equal, combining equation (b) and equation (c), we get: （d） By integrating formula (d) over time, we can obtain the relationship between the water infiltration time t and the water infiltration depth z of the upper soil layer: （e）。 3. The method for determining the double-layer soil water infiltration depth according to claim 2, characterized in that: The method for determining the time length of water infiltration to the interface between the upper and lower soil layers comprises: Substituting the upper soil thickness z1 as the upper soil water infiltration depth z into formula (e), the time t1 for water to infiltrate to the interface between the upper and lower soil layers is calculated as: （j）。 4. The method for determining the double-layer soil water infiltration depth according to claim 1, characterized in that: The relationship between the water infiltration time and the water infiltration depth in the lower soil layer determined according to the infiltration pressure at the interface between the upper and lower soil layers includes: If the infiltration water volume at the interface between the upper and lower soil layers is less than or equal to the infiltrative water volume of the lower soil layer, the infiltration pressure at the interface between the upper and lower soil layers is determined to be less than or equal to zero. Based on the relevant parameters of the upper soil layer, the infiltration rate i2 of the lower soil surface is calculated as: （k） Where H is the thickness of water on the upper soil surface; K s1 、s f1 are the saturated permeability coefficient of the upper soil and the matrix suction corresponding to the initial water content; z1 is the thickness of the upper soil; According to the saturated water content of the lower soil θ s2 , initial water contentθ i2 , and the water infiltration time t2 in the lower soil, the expression of the infiltration rate i2 on the lower soil surface is reconstructed: （l） Where z2 is the water infiltration depth in the lower soil layer; Since the water infiltration rate of the same layer of homogeneous soil is equal, the combination of equation (k) and equation (l) yields: （m） After integrating the formula (m) over time, the relationship between the water infiltration time t2 and the water infiltration depth z2 in the lower soil layer is obtained as follows: （n）。 5. The method for determining the double-layer soil water infiltration depth according to claim 1 or 4, characterized in that: The relationship between the water infiltration time and the water infiltration depth in the lower soil layer determined according to the infiltration pressure at the interface between the upper and lower soil layers includes: If the infiltration water volume at the interface between the upper and lower soil layers is greater than the infiltrative water volume of the lower soil layer, it is determined that the infiltration pressure at the interface between the upper and lower soil layers is greater than zero. The surface infiltration rate i2 of the lower soil layer is calculated based on the predetermined infiltration pressure at the interface between the upper and lower soil layers and the relevant parameters of the lower soil layer. The expression is: （p） Where H' is the infiltration pressure at the interface between the upper and lower soil layers; z2 is the water infiltration depth in the lower soil layer; K s2 、s f2 are the saturated permeability coefficient of the underlying soil and the matrix suction corresponding to the initial water content; According to the saturated water content of the lower soil θ s2 , initial water contentθ i2 , and the water infiltration time t2 in the lower soil, the expression of the infiltration rate i2 on the lower soil surface is reconstructed as: （q） Since the water infiltration rate of the same layer of homogeneous soil is equal, the combination of equation (p) and equation (q) yields: （r） By integrating both sides of formula (r) over time, we can obtain the relationship between the water infiltration time t2 and the water infiltration depth z2 in the lower soil layer: （s）。 6. The method for determining the double-layer soil water infiltration depth according to claim 1, characterized in that: The method for determining the infiltration pressure at the interface between the upper and lower soil layers comprises: Assuming that the water infiltration depth in the lower soil layer is equal to the thickness z1 of the upper soil layer, the positions of the upper and lower soil layers are interchanged, and all other parameters remain unchanged. Based on the equal energy lost in the infiltration of the upper and lower soil layers, the time required to reach the infiltration depth is also equal. Therefore, the equivalent equation is established as follows: （w） Where H´ is the infiltration pressure at the interface between the upper and lower soil layers; θ s1 ,θ i1 are the saturated water content and initial water content of the upper soil layer, respectively; K s1 、s f1 are the saturated permeability coefficient of the upper soil and the matrix suction corresponding to the initial water content; θ s2 ,θ i2 are the saturated water content and initial water content of the lower soil layer, respectively; K s2 、s f2 are the saturated permeability coefficient of the lower soil and the matrix suction corresponding to the initial water content; H is the thickness of the surface water in the upper soil layer; The infiltration pressure H´ at the interface between the upper and lower soil layers is calculated using formula (w).

7. A device for determining the water infiltration depth of a double-layer soil, characterized in that: The device comprises: Collection module: used to collect water infiltration duration; The first determination module: in response to the collected water infiltration time being less than or equal to the predetermined time for water to infiltrate to the interface between the upper and lower soil layers, indicating that water has not infiltrated into the lower soil layer, the water infiltration depth of the double-layer soil is determined according to the relationship between the water infiltration time in the upper soil layer and the water infiltration depth; The second determination module: in response to the collected water infiltration time being greater than the predetermined time for water to infiltrate to the interface between the upper and lower soil layers, indicating that the water has infiltrated into the lower soil layer, the relationship between the water infiltration time and the water infiltration depth in the lower soil layer is determined according to the infiltration pressure at the interface between the upper and lower soil layers, and the water infiltration depth of the double-layer soil is determined in combination with the relationship between the water infiltration time and the water infiltration depth in the upper soil layer.

8. A computer device, characterized in that: The method comprises a processor and a memory connected to the processor, wherein a computer program is stored in the memory, and when the computer program is executed by the processor, the steps of the method according to any one of claims 1 to 6 are executed.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the steps of the method described in any one of claims 1 to 6 are implemented.

10. A computer program product comprising a computer program / instructions, characterized in that When the computer program / instructions are executed by a processor, the steps of the method described in any one of claims 1 to 6 are implemented.

Images