Simulation method, device and equipment for coupling migration of aeration zone, medium and product

Abstract

The invention provides an aeration zone coupling migration simulation method, device and equipment, a medium and a product. The method comprises the following steps: acquiring atmosphere information of an arid and semi-arid region and basic physical information of a soil body within a preset time; determining a moisture motion equation, a heat motion equation and a solute transport equation according to the basic physical information, and constructing a coupling transport model; calculating and inverting hydraulic parameters of the soil body according to the basic physical information; determining an initial condition and a boundary condition of the coupling migration model according to the atmosphere information and the basic physical information; and solving the coupling migration model by utilizing hydraulic parameters, initial conditions and boundary conditions of the soil body, and simulating a coupling migration process of water, heat and solute in the aeration zone within a preset time. The device is used for simulating the aeration zone water, heat and solute coupling migration process under the evaporation condition of the arid and semi-arid region within the preset time period, and quantitatively analyzing the influence mode of the solute on aeration zone soil moisture movement.

Description

Technical Field

[0001] The present invention relates to the field of environmental science and technology, and in particular to a simulation method, device, equipment, medium and product for coupled migration of aeration zones. Background Art

[0002] With the expansion of human living area and the progress of production activities, soil salinization is becoming more and more serious. The salt entering the soil not only changes the soil structure, but also affects the movement and distribution of water in the soil. It not only affects the movement of soil liquid water through the osmotic potential generated in the soil, but also affects the movement of soil water vapor by affecting the relative humidity of water vapor in the soil.

[0003] At present, the study of groundwater seepage in the vadose zone often ignores the influence of solutes on the movement of liquid water and water vapor. At the same time, due to the complexity of the model equations, the influence of temperature and osmotic potential on the relative humidity of soil water vapor, the influence of solute concentration on film water thickness and surface tension, and the influence of temperature and matrix potential on osmotic potential are usually ignored in existing studies. All of these make the prediction of groundwater movement in arid and semi-arid salinized areas inaccurate.

[0004] Therefore, there is currently a lack of a simulation method for the coupled migration of water, heat, and solutes in the vadose zone with strong system correlation and clear logical structure, which can be used to simulate and analyze the coupled migration process of water, heat, and solutes in the vadose zone. Summary of the invention

[0005] In view of the above problems, the present invention is proposed to provide a method, device, equipment, medium and product for simulating the coupled migration of water, heat and solute in the vadose zone that overcomes the above problems or at least partially solves the above problems.

[0006] To achieve the above and other related purposes, the present invention provides a method for simulating coupled migration of aeration zones, the method comprising:

[0007] Obtaining atmospheric information and basic physical information of soil in arid and semi-arid areas within a preset time; wherein the atmospheric information includes air temperature, air humidity, solar radiation and water vapor saturation; the basic physical information includes soil moisture content, soil temperature, soil matrix potential and soil conductivity;

[0008] Determine the water motion equation, the heat motion equation and the solute transport equation according to the basic physical information, and construct a coupled transport model using the water motion equation, the heat motion equation and the solute transport equation;

[0009] Calculating and inverting soil hydraulic parameters according to the basic physical information; and determining initial conditions and boundary conditions of the coupled transport model according to the atmospheric information and the basic physical information;

[0010] The coupled transport model is solved by using the soil hydraulic parameters, the initial conditions and the boundary conditions, and the coupled transport process of water, heat and solute in the vadose zone within a preset time is simulated.

[0011] Optionally, the calculating and inverting soil hydraulic parameters according to the basic physical information includes:

[0012] Draw a soil moisture characteristic curve according to the soil moisture content and the soil matrix potential in the basic physical information, and determine a relationship equation between the soil moisture content and the soil matrix potential according to the soil moisture characteristic curve;

[0013] The hydraulic conductivity of liquid water in the vadose zone is inverted and calculated using the relational equation.

[0014] Optionally, determining the initial conditions of the coupled transport model according to the atmospheric information and basic physical information includes:

[0015] The soil moisture content and the soil temperature in the basic physical information are extracted, and the soil moisture content and the soil temperature are set as initial conditions of the coupled transport model.

[0016] Optionally, determining the boundary conditions of the coupled transport model according to the atmospheric information and basic physical information includes:

[0017] The soil temperature in the basic physical information is extracted, and the soil temperature is set as a boundary condition of the temperature field in the coupled transport model.

[0018] Optionally, determining the boundary conditions of the coupled transport model according to the atmospheric information and the basic physical information includes:

[0019] The evaporation efficiency of the soil is calculated according to the atmospheric information and the basic physical information, and the evaporation rate is set as the atmospheric boundary condition of the coupled transport model.

[0020] Optionally, solving the coupled transport model by using the soil hydraulic parameters, the initial conditions and the boundary conditions includes:

[0021] The coupled transport model is numerically simulated by using the soil hydraulic parameters, the initial conditions and the boundary conditions to obtain a temperature simulation value and a water content simulation value.

[0022] In a second aspect, the present invention further provides a simulation device for coupled migration of aeration zones, the device comprising:

[0023] An acquisition module is used to acquire atmospheric information and basic physical information of soil in arid and semi-arid areas within a preset time; wherein the atmospheric information includes air temperature, air humidity, solar radiation and water vapor saturation; and the basic physical information includes soil moisture content, soil temperature, soil matrix potential and soil conductivity;

[0024] A construction module is used to determine the water motion equation, the heat motion equation and the solute transport equation by using the basic physical information, and to construct a coupled transport model by using the water motion equation, the heat motion equation and the solute transport equation;

[0025] A determination module, used to calculate and invert soil hydraulic parameters according to the basic physical information; and to determine the initial conditions and boundary conditions of the coupled transport model according to the atmospheric information and the basic physical information;

[0026] The calculation module is used to solve the coupled transport model through the soil hydraulic parameters, the initial conditions and the boundary conditions, and simulate the coupled transport process of water, heat and solute in the vadose zone within a preset time.

[0027] In a third aspect, the present invention provides an electronic device, comprising: a memory and a processor; the memory is used to store a computer program; the processor is used to execute the computer program stored in the memory, so that the electronic device performs the steps of the simulation method of coupled migration of the vadose zone as described above.

[0028] In a fourth aspect, the present invention provides a computer-readable storage medium having a computer program stored thereon, which, when executed by an electronic device, implements the steps of the method for simulating coupled migration of aeration zones as described above.

[0029] In a fifth aspect, the present invention provides a computer program product, which includes computer program code. When the computer program code runs on a computer, the computer implements the steps of the simulation method of coupled migration of vadose zones as described above.

[0030] The above one or more technical solutions provided by the present invention may have the following advantages or at least achieve the following technical effects:

[0031] The simulation method, device and equipment, medium and product of coupled migration of vadose zone of the present invention, the method comprises: obtaining atmospheric information and basic physical information of soil in arid and semi-arid areas within a preset time; wherein the atmospheric information comprises air temperature, air humidity, solar radiation and water vapor saturation; the basic physical information comprises soil moisture content, soil temperature, soil matrix potential and soil conductivity; according to the basic physical information, determining the water motion equation, heat motion equation and solute migration equation, and constructing a coupled migration model by using the water motion equation, heat motion equation and solute migration equation; according to the basic physical information, calculating and inverting soil hydraulic parameters; and according to the atmospheric information and the basic physical information, determining the initial conditions and boundary conditions of the coupled migration model; using the soil hydraulic parameters, the initial conditions and the boundary conditions, solving the coupled migration model, and simulating the coupled migration process of water, heat and solute in the vadose zone within the preset time.

[0032] Therefore, the present invention is used to simulate the coupled migration process of water, heat and solute in the vadose zone under evaporation conditions in arid and semi-arid areas within a preset time period, and quantitatively analyze the influence pattern of solute on soil moisture movement in the vadose zone, which can improve the accuracy of groundwater movement prediction in arid and semi-arid areas. BRIEF DESCRIPTION OF THE DRAWINGS

[0033] Figure 1 It is a schematic diagram showing the principle of a method for simulating coupled migration of water, heat and solute in a vadose zone according to one embodiment of the present invention;

[0034] Figure 2 Shown is a schematic flow chart of a method for simulating coupled migration of water, heat and solute in a vadose zone according to an embodiment of the present invention;

[0035] Figure 3 A schematic diagram showing a water content simulation result in one embodiment of the present invention;

[0036] Figure 4 Shown is a schematic diagram of a temperature simulation result in one embodiment of the present invention;

[0037] Figure 5 It is a schematic diagram of the functional modules of a simulation device for coupled transport of water, heat and solute in the vadose zone according to one embodiment of the present invention;

[0038] Figure 6 Shown is a schematic diagram of an electronic device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0039] The following describes the embodiments of the present invention by specific examples, and those skilled in the art can easily understand other advantages and effects of the present invention from the contents disclosed in this specification. The present invention can also be implemented or applied through other different specific embodiments, and the details in this specification can also be modified or changed in various ways based on different viewpoints and applications without departing from the spirit of the present invention. It should be noted that the following embodiments and features in the embodiments can be combined with each other without conflict.

[0040] It should be noted that the illustrations provided in the following embodiments are only schematic illustrations of the basic concept of the present invention, and thus the drawings only show components related to the present invention rather than being drawn according to the number, shape and size of components in actual implementation. In actual implementation, the type, quantity and proportion of each component may be changed arbitrarily, and the component layout may also be more complicated.

[0041] In the following description, numerous details are discussed to provide a more thorough explanation of the embodiments of the present invention. However, it is obvious to those skilled in the art that the embodiments of the present invention can be implemented without these specific details. In other embodiments, well-known structures and devices are shown in the form of block diagrams rather than in detail to avoid making the embodiments of the present invention difficult to understand.

[0042] The terms "first", "second", etc. in the specification and claims of the disclosed embodiments and the above-mentioned drawings are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the terms used in this way can be interchanged where appropriate, so that the embodiments of the disclosed embodiments described herein. In addition, the terms "including" and "having" and any variations thereof are intended to cover non-exclusive inclusions.

[0043] Unless otherwise stated, the term "plurality" means two or more.

[0044] In the embodiment of the present disclosure, the character " / " indicates that the preceding and following objects are in an "or" relationship. For example, A / B indicates: A or B.

[0045] The term "and / or" is a description of the association relationship between objects, indicating that three relationships can exist. For example, A and / or B means: A or B, or, A and B.

[0046] For ease of understanding, the following first explains the relevant nouns or terms used in the embodiments of the present invention:

[0047] <1> , the vadose zone, also known as the unsaturated zone or unsaturated soil layer, refers to the area below the surface and above the groundwater level.

[0048] The migration of soil moisture, heat and solute in the vadose zone is a mutually influencing and coupled process. Based on the interaction relationship and influence mode between the three physical fields, the three-field coupling control equations (the water movement equation, the heat movement equation and the solute migration equation) were established.

[0049] <2> , coupled transport model, which describes the interaction and migration of water, heat and solute in the soil; it includes the water motion equation, heat motion equation and solute transport equation.

[0050] <3> , inversion method, which is used to determine unknown input parameters from known results (such as observational data).

[0051] <4> , Soil matrix potential, also known as soil matrix suction or soil water potential, is a physical quantity that describes the energy state of water in the soil. It is an important driving force for soil water movement and reflects the ability of soil particles to adsorb water. Soil matrix potential is a negative value, indicating that the potential energy of water in the soil is lower than the potential energy of water in the atmosphere.

[0052] <5> , Soil Water Characteristic Curve (SWCC), also known as water release curve or water retention curve, describes the relationship between soil moisture content and soil matrix potential.

[0053] In a saturated state, the soil matrix potential is close to zero and the soil moisture content is at its highest. As the soil dries, the soil matrix lowers (becomes more negative) and the soil moisture content gradually decreases.

[0054] <6> , hydraulic conductivity, also known as permeability, is a physical quantity that describes the water flow capacity in soil, rock or other porous media; it is one of the key parameters controlling liquid water flux. By measuring this parameter, the movement of soil moisture can be better predicted and controlled.

[0055] The technical solutions in the embodiments of the present invention will be described in detail below in conjunction with the accompanying drawings in the embodiments of the present invention.

[0056] See also Figure 1 , Figure 1 Schematic diagram showing the principle of the simulation method for the coupled transport of water, heat and solute in the vadose zone; see Figure 1As shown in the figure, in the coupled transport model (involving the water movement equation, heat movement equation and solute transport equation), the water movement equation (including liquid water movement and water vapor movement), heat movement equation and solute transport equation and numerical method can be used to describe and simulate to accurately predict the migration and distribution of water, heat and solute in the soil. For example, when water flows in the soil, the heat of the solute in the water also moves with it. This phenomenon is called convective heat transfer. The faster the water flows, the faster the heat transfer.

[0057] See also Figure 2 An embodiment of the present invention provides a method for simulating coupled migration of water, heat, and solute in a vadose zone. The method may include the following steps S10 to S40:

[0058] Step S10, obtaining atmospheric information and basic physical information of soil in arid and semi-arid areas within a preset time; wherein the atmospheric information includes air temperature, air humidity, solar radiation and water vapor saturation; and the basic physical information includes soil moisture content, soil temperature, soil matrix potential and soil electrical conductivity.

[0059] The preset time can be set by researchers according to actual research needs. In this embodiment, the preset time is preferably set to 150 hours.

[0060] Atmospheric information includes but is not limited to: air temperature, air humidity, solar radiation, and water vapor saturation.

[0061] Basic physical information includes but is not limited to: soil moisture content, soil temperature, soil matrix potential, and soil electrical conductivity.

[0062] In a specific implementation, atmospheric information of arid and semi-arid areas within a preset time can be obtained through atmospheric sensors, and basic information of soil in arid and semi-arid areas within a preset time can be collected through soil sensors; wherein the atmospheric information may include air temperature, air humidity, solar radiation and water vapor saturation; basic physical information may include soil moisture content, soil temperature, soil matrix potential and soil conductivity.

[0063] Step S20, determining the water motion equation, the heat motion equation and the solute transport equation according to the basic physical information, and constructing a coupled transport model using the water motion equation, the heat motion equation and the solute transport equation.

[0064] Among them, the water movement equation can be used to describe the water movement in unsaturated and saturated soils; it can be expressed as:

[0065]

[0066] Among them, θ land θ v are the volume content of liquid water and the volume content of water vapor (m 3 m -3 );q L is the liquid water flux (ms -1 );q v is the water vapor flux (ms -1 ), t is time (s), and z is the position coordinate (m).

[0067] Among them, the liquid water flux (q L ), the volume of liquid water passing through a unit area per unit time; it can be expressed as:

[0068]

[0069] Where K is the hydraulic conductivity of liquid water in the vadose zone (ms -1 ); γ is the surface tension of water (gs -2 ); γ0 is the surface tension of water at 25°C (gs -2 );h l is the matrix potential (m); σ is the reduction factor (dimensionless); v is the number of ions per molecule of the ionized solute (dimensionless); g is the gravitational acceleration (=9.81ms -2 ); R is the ideal gas constant (=8.314 Jmol -1 K -1 ), φ is the osmotic potential coefficient (dimensionless).

[0070] Among them, the water vapor flux (q v ), refers to the volume of water vapor passing through a unit area per unit time; it can be expressed as:

[0071]

[0072] Where D is the water vapor diffusion coefficient in soil; H rl is the relative humidity affected by matrix potential (dimensionless); H ro is the relative humidity affected by osmotic potential (dimensionless).

[0073] Among them, the water vapor volume content (θ v ) can be expressed as:

[0074]

[0075] Among them, H r is the relative humidity of water vapor in the soil (dimensionless); ρ sv is the saturated water vapor density (kg m -3 );θ a is the volume content of gas in soil pores (m3 m -3 );ρ l is the density of liquid water (kg m -3 ).

[0076] The heat equation of motion can describe the process of heat transfer in the soil, including heat conduction and convection; it can be expressed as:

[0077]

[0078] Where, T is temperature (K); C s , C v , C w are the volumetric heat capacities of soil particles, liquid water, and water vapor (J m -3 K -1 ), T0 is an arbitrary reference temperature, θ n is the volume fraction of the solid phase in the soil (m 3 m -3 ); L0 is the volumetric latent heat of evaporation of liquid water in the soil at time T0 (J m -3 ); λ is the apparent thermal conductivity of soil (W m -1 K -1 ).

[0079] The solute transport equation can describe the transport of solutes in soil moisture; it can be expressed as:

[0080]

[0081] Where, C is the solute concentration (mol / kg); q s is the salt flux (kg m -2 s -1 ).

[0082] Among them, the salt flux (q s ), refers to the total amount of salt passing through a unit area per unit time; it can be expressed as:

[0083]

[0084] Among them, D i is the molecular diffusion coefficient (m 2 s -1 );D m is the mechanical diffusion coefficient (m 2 s -1 ).

[0085] The coupled transport model can integrate the water movement equation, heat movement equation and solute transport equation into a unified framework to simulate the interaction and migration of water, heat and solute in soil.

[0086] In the specific implementation, the water movement equation, heat movement equation and solute transport equation can be determined according to the basic physical information of the soil; then the water movement equation, heat movement equation and solute transport equation can be integrated into a unified framework to construct the corresponding coupled transport model, so that the coupled transport model can be used to simulate the interaction and migration of water, heat and solute in the soil.

[0087] Step S30, calculating and inverting soil hydraulic parameters according to the basic physical information; and determining initial conditions and boundary conditions of the coupled transport model according to the atmospheric information and the basic physical information.

[0088] The soil hydraulic parameters may be a series of key coefficients describing the water movement characteristics of the soil; they may include: water characteristic curve, hydraulic conductivity, saturated hydraulic conductivity, etc.

[0089] Initial conditions refer to the state of each point in the coupled transport model at the beginning of the simulation.

[0090] As an example, the initial water content and soil temperature etc. can be set for the coupled transport model.

[0091] Boundary conditions define the state on the boundary of the coupled transport model, which can be set by researchers according to actual research needs. In this embodiment, the boundary conditions may include: boundary conditions of the temperature field and (upper and lower) boundary conditions of the seepage field.

[0092] The temperature field boundary condition is a boundary condition that specifies the temperature value on the boundary; it can be used to simulate the heat exchange between the soil surface and the atmosphere.

[0093] As an example, the first type of boundary condition of the temperature field can be set according to the measured temperature value of the soil surface.

[0094] The upper boundary condition of the seepage field can be set as the atmospheric boundary condition, which is the second type of boundary condition expressed as the evaporation rate of the soil; it can be used to simulate the water exchange between the soil surface and the atmosphere, such as evaporation.

[0095] The evaporation rate (E) of soil can be expressed as:

[0096]

[0097] Among them, ρ sa is the water vapor density in saturated air (kg m -3 );r a is the air resistance coefficient (sm -1 );r sis the soil resistance coefficient (sm -1 );H ra is the relative humidity of water vapor in the air (m 3 m -3 ).

[0098] The lower boundary condition of the seepage field can be set as a constant water level boundary condition.

[0099] It should be noted that the boundary conditions in the coupled transport model are usually combined with the boundary conditions of soil moisture and solute transport to achieve a comprehensive simulation of the interaction between moisture, heat and solutes in the soil system.

[0100] In a specific implementation, the hydraulic parameters of the soil can be calculated using the inversion method based on the basic physical information; and the initial conditions and boundary conditions of the coupled transport model can be determined based on the atmospheric information and the basic physical information.

[0101] Step S40, using the soil hydraulic parameters, the initial conditions and the boundary conditions, the coupled transport model is solved, and the coupled transport process of water, heat and solute in the vadose zone within a preset time is simulated.

[0102] In a specific implementation, the coupled transport model can be solved by a fully coupled numerical method (e.g., finite element method) according to the hydraulic parameters, initial conditions, and boundary conditions of the soil, and the simulated values ​​of the moisture content and temperature can be obtained (see formula (2), formula (3), and formula (5) in step S20 above); and the interaction and migration of water, heat, and solutes in the soil can be simulated through the changes in moisture content and temperature.

[0103] Further, in one embodiment, step S30 may include the following sub-step A10:

[0104] Sub-step A10: extracting soil moisture content and soil temperature from the basic physical information, and setting the soil moisture content and soil temperature as initial conditions of the coupled transport model.

[0105] The initial conditions can be determined by:

[0106] ① Experimental data: Initial conditions (e.g., soil temperature and soil moisture content) can be obtained through field or laboratory tests.

[0107] ②Historical data: It can be the use of previously simulated or observed data as initial conditions (such as soil temperature and soil moisture content).

[0108] In a specific implementation, the soil moisture content and soil temperature in the basic physical information can be extracted, and the soil temperature and soil moisture content can be set as the initial conditions of the coupled transport model.

[0109] Furthermore, in one embodiment, step S30 may further include the following sub-step A20:

[0110] Sub-step A20: extracting soil temperature from the basic physical information and setting the soil temperature as a boundary condition of the temperature field in the coupled transport model.

[0111] In a specific implementation, the soil temperature (ie, the measured temperature value) can be extracted from the basic physical information, and then the measured temperature value can be set as the boundary condition of the temperature field in the coupled transport model to simulate the heat exchange between the soil surface and the atmosphere.

[0112] Furthermore, in one embodiment, step S30 may also include the following sub-step A30:

[0113] Sub-step A30, calculating the evaporation efficiency of the soil according to the atmospheric information and the basic physical information, and setting the evaporation rate as the atmospheric boundary condition of the coupled transport model.

[0114] In a specific implementation, the evaporation rate of the soil can be calculated using atmospheric information and basic physical information (see formula (5) in step S30 above), and then the evaporation rate of the soil is set as the atmospheric boundary condition of the coupled transport model to simulate the water exchange between the soil surface and the atmosphere.

[0115] In this embodiment, atmospheric information and basic physical information of soil in arid and semi-arid areas within a preset time are obtained; based on the basic physical information, the water movement equation, heat movement equation and solute transport equation are determined, and a coupled transport model is constructed using the water movement equation, heat movement equation and solute transport equation; based on the basic physical information, soil hydraulic parameters are calculated and inverted; and based on the atmospheric information and the basic physical information, the initial conditions and boundary conditions of the coupled transport model are determined; the coupled transport model is solved using the soil hydraulic parameters, the initial conditions and the boundary conditions, and the coupled transport process of water, heat and solute in the vadose zone within a preset time is simulated. The present invention is used to simulate the coupled transport process of water, heat and solute in the vadose zone under evaporation conditions in arid and semi-arid areas within a preset time period, quantitatively analyze the influence mode of solute on the movement of soil water in the vadose zone, and improve the accuracy of groundwater movement prediction in arid and semi-arid areas.

[0116] Based on the above embodiments, a second embodiment of the solar irradiance prediction method of the present invention is proposed. In this embodiment, step S30 may further include the following sub-steps S301 to S302:

[0117] Sub-step S301, drawing a soil moisture characteristic curve according to the soil moisture content and the soil matrix potential in the basic physical information, and determining a relationship equation between the soil moisture content and the soil matrix potential according to the soil moisture characteristic curve.

[0118] Among them, the relationship equation between soil matrix potential and soil moisture content can be expressed as:

[0119]

[0120] Among them, θ r is the residual soil moisture content (cm 3 cm -3 );θ s is the soil saturated moisture content (cm 3 cm -3 ); α is a parameter related to the inverse of the air intake value of the soil moisture characteristic curve (1 / cm); n and m are empirical parameters (dimensionless) reflecting the soil pore distribution, m = 1-1 / n. In this embodiment, θ r =0.027,θ s =0.304, α = 0.02368 cm -1 , n=5.42.

[0121] In a specific implementation, a soil moisture characteristic curve can be drawn according to the soil moisture content and the soil matrix potential; and then according to the soil moisture characteristic curve, the relationship between the soil moisture content and the soil matrix potential can be fitted into a relationship equation.

[0122] Sub-step S302, using the relationship equation, invert and calculate the hydraulic conductivity of liquid water in the aeration zone.

[0123] Among them, the hydraulic conductivity coefficient (K) is a physical quantity that describes the water flow capacity in the soil and is one of the key parameters controlling the liquid water flux; it can be expressed as:

[0124]

[0125] Among them, K s is the saturated hydraulic conductivity (ms -1 );S e is the effective saturation (dimensionless); l is an empirical parameter (dimensionless). In this embodiment, K s =1.9×10 -5 , l=0.5.

[0126] In a specific implementation, the relationship equation between soil moisture content and soil matrix potential can be used to determine the parameters α and n, and calculate the hydraulic conductivity coefficient of liquid water in the vadose zone.

[0127] In this embodiment, a soil moisture characteristic curve is drawn according to the soil moisture content and soil matrix potential in the basic physical information, and the relationship equation between the soil moisture content and the soil matrix potential is determined according to the soil moisture characteristic curve; the hydraulic conductivity coefficient of the liquid water in the vadose zone is calculated and inverted using the relationship equation. Thus, by calculating the hydraulic conductivity coefficient of the liquid water in the vadose zone, the liquid water flux in the soil moisture movement can be determined, and the movement of soil moisture can be better predicted.

[0128] Based on the above embodiments, a third embodiment of the solar irradiance prediction method of the present invention is proposed. In this embodiment, step S40 may include the following sub-step S401:

[0129] Sub-step S401, using the soil hydraulic parameters, the initial conditions and the boundary conditions, numerical simulation calculations are performed on the coupled transport model to obtain temperature simulation values ​​and water content simulation values.

[0130] The temperature simulation value may be the soil temperature obtained by fully coupling and numerically solving the coupled transport model (including the water transport equation, the heat transport equation, and the solute transport equation).

[0131] The simulated value of moisture content may be the soil moisture content obtained by performing a fully coupled numerical solution to a coupled transport model (including: a water movement equation, a heat movement equation, and a solute transport equation).

[0132] In a specific implementation, the soil hydraulic parameters, initial conditions and boundary conditions can be used to solve the water motion equation, heat motion equation and solute transport equation by a fully coupled numerical method, and the temperature simulation value (see formula (5) in step S20 above) and the water content simulation value (see formula (2) and formula (3) in step S20 above) can be calculated. Then, the temperature simulation value and the water content simulation value are analyzed to evaluate the distribution and migration process of water, capacity and solute in the soil.

[0133] As an example, in order to verify the simulation effect of the coupled transport model, the simulated values ​​(temperature simulation values ​​and moisture content simulation values) and the measured values ​​(soil temperature and moisture content) are compared and fitted. The comparison results are shown in Figure 3 and Figure 4 , Figure 3 Shown is a schematic diagram of the water content simulation results. Figure 4 It shows a schematic diagram of the temperature simulation results; it can be seen that the simulated values ​​of the coupled transport model (soil temperature and moisture content) and the measured values ​​have a good degree of fit, which better reflects the change pattern of soil moisture and salinity within the preset time. The simulation results are reliable and can be used for simulation and analysis of the coupled transport process of water, heat and solute in the vadose zone.

[0134] Furthermore, the liquid water flux, water vapor flux and salt flux can be calculated based on the temperature simulation value, water content simulation value and hydraulic conductivity coefficient obtained. Thus, the migration process of liquid water, water vapor and salt in the soil can be simulated through the coupled transport model.

[0135] In this embodiment, the coupled transport model is numerically simulated using the soil hydraulic parameters, the initial conditions, and the boundary conditions to obtain the temperature simulation value and the water content simulation value. Thus, the dynamic behavior of water, heat, and solute migration in the soil is reflected through the changes in the temperature simulation value and the water content simulation value, thereby realizing the simulation and analysis of the coupled transport process of water, heat, and solute in the vadose zone using the coupled transport model.

[0136] Based on the same inventive concept, the fourth embodiment of the present invention also provides a simulation device for coupled migration of vadose zone corresponding to the simulation method for coupled migration of vadose zone in the previous embodiment. Since the principle of solving the problem by the device in the fourth embodiment of the present invention is similar to the simulation method for coupled migration of vadose zone in the previous embodiment of the present invention, the implementation of the device can refer to the implementation of the method, and the repeated parts will not be repeated. Figure 5 The simulation device for coupled migration of the aeration zone of the present invention may include:

[0137] The acquisition module 10 is used to acquire the atmospheric information and basic physical information of the soil in the arid and semi-arid areas within a preset time; wherein the atmospheric information includes air temperature, air humidity, solar radiation and water vapor saturation; the basic physical information includes soil moisture content, soil temperature, soil matrix potential and soil conductivity;

[0138] A construction module 20 is used to determine the water motion equation, the heat motion equation and the solute transport equation by using the basic physical information, and to construct a coupled transport model by using the water motion equation, the heat motion equation and the solute transport equation;

[0139] A determination module 30 is used to calculate and invert soil hydraulic parameters according to the basic physical information; and to determine initial conditions and boundary conditions of the coupled transport model according to the atmospheric information and the basic physical information;

[0140] The calculation module 40 is used to solve the coupled transport model through the soil hydraulic parameters, the initial conditions and the boundary conditions, and simulate the coupled transport process of water, heat and solute in the vadose zone within a preset time.

[0141] In addition, the present invention also provides a computer-readable storage medium, on which a computer program is stored, and when the program is executed by a processor, the above-mentioned simulation method of coupled migration of the vadose zone is implemented.

[0142] Figure 6 is a schematic block diagram of an electronic device provided in an embodiment of the present application. Figure 6 As shown, the electronic device includes: at least one processor 401, a memory 402, at least one network interface 403 and a user interface 405. The various components in the electronic device are coupled together through a bus system 404. It can be understood that the bus system 404 is used to realize the connection and communication between these components. In addition to the data bus, the bus system 404 also includes a power bus, a control bus and a status signal bus. However, for the sake of clarity, Figure 6 In the specification, various buses are labeled as bus systems.

[0143] The user interface 405 may include a display, a keyboard, a mouse, a trackball, a click gun, keys, buttons, a touch pad or a touch screen.

[0144] It is understood that the memory 402 can be a volatile memory or a non-volatile memory, and can also include both volatile and non-volatile memories. Among them, the non-volatile memory can be a read-only memory (ROM), a programmable read-only memory (PROM), which is used as an external cache. By way of example but not limitation, many forms of RAM are available, such as static random access memory (SRAM), synchronous static random access memory (SSRAM). The memory described in the embodiments of the present invention is intended to include but is not limited to these and any other suitable categories of memory.

[0145] The memory 402 in the embodiment of the present invention is used to store various categories of data to support the operation of the electronic device 400. Examples of these data include: any executable program for operating on the electronic device 400, such as an operating system 4021 and an application 4022; the operating system 4021 includes various system programs, such as a framework layer, a core library layer, a driver layer, etc., for implementing various basic services and processing hardware-based tasks. The application 4022 may include various applications, such as a media player (Media Player), a browser (Browser), etc., for implementing various application services. The method for determining the true height of a web page provided in the embodiment of the present invention may be included in the application 4022.

[0146] The method disclosed in the above embodiment of the present invention can be applied to the processor 401, or implemented by the processor 401. The processor 401 may be an integrated circuit chip with signal processing capabilities. In the implementation process, each step of the above method can be completed by the hardware integrated logic circuit in the processor 401 or the instruction in the form of software. The above processor 401 may be a general processor, a digital signal processor (DSP, Digital Signal Processor), or other programmable logic devices, discrete gates or transistor logic devices, discrete hardware components, etc. The processor 401 can implement or execute the methods, steps and logic block diagrams disclosed in the embodiment of the present invention. The general processor 401 can be a microprocessor or any conventional processor, etc. In combination with the steps of the method for determining the real height of a web page provided in the embodiment of the present invention, it can be directly embodied as a hardware decoding processor to perform, or it can be performed by a combination of hardware and software modules in the decoding processor. The software module can be located in a storage medium, which is located in a memory, and the processor reads the information in the memory and completes the steps of the above method in combination with its hardware.

[0147] In an exemplary embodiment, the electronic device 400 may be implemented by one or more application specific integrated circuits (ASIC), DSP, programmable logic device (PLD), complex programmable logic device (CPLD) to execute the aforementioned method.

[0148] In summary, the present invention is used to simulate the coupled migration process of water, heat and solute in the vadose zone under evaporation conditions in arid and semi-arid areas within a preset time period, and quantitatively analyze the influence pattern of solute on the movement of soil moisture in the vadose zone, which can improve the accuracy of groundwater movement prediction in arid and semi-arid areas.

[0149] The above embodiments are merely illustrative of the principles and effects of the present invention, and are not intended to limit the present invention. Anyone familiar with the art may modify or alter the above embodiments without departing from the spirit and scope of the present invention. Therefore, all equivalent modifications or alterations made by a person of ordinary skill in the art without departing from the spirit and technical concept disclosed by the present invention shall still be covered by the claims of the present invention.

Claims

1. A method for simulating coupled migration of aeration zones, characterized in that: The method comprises: Obtaining atmospheric information and basic physical information of soil in arid and semi-arid areas within a preset time; wherein the atmospheric information includes air temperature, air humidity, solar radiation and water vapor saturation; the basic physical information includes soil moisture content, soil temperature, soil matrix potential and soil conductivity; Determine the water motion equation, the heat motion equation and the solute transport equation according to the basic physical information, and construct a coupled transport model using the water motion equation, the heat motion equation and the solute transport equation; Calculating and inverting soil hydraulic parameters according to the basic physical information; and determining initial conditions and boundary conditions of the coupled transport model according to the atmospheric information and the basic physical information; The coupled transport model is solved by using the soil hydraulic parameters, the initial conditions and the boundary conditions, and the coupled transport process of water, heat and solute in the vadose zone within a preset time is simulated.

2. The method according to claim 1, characterized in that The calculation and inversion of soil hydraulic parameters according to the basic physical information includes: Draw a soil moisture characteristic curve according to the soil moisture content and the soil matrix potential in the basic physical information, and determine a relationship equation between the soil moisture content and the soil matrix potential according to the soil moisture characteristic curve; The hydraulic conductivity of liquid water in the vadose zone is inverted and calculated using the relational equation.

3. The method according to claim 1, characterized in that Determining the initial conditions of the coupled transport model according to the atmospheric information and basic physical information includes: The soil moisture content and the soil temperature in the basic physical information are extracted, and the soil moisture content and the soil temperature are set as initial conditions of the coupled transport model.

4. The method according to claim 1, characterized in that: Determining the boundary conditions of the coupled transport model according to the atmospheric information and basic physical information includes: The soil temperature in the basic physical information is extracted, and the soil temperature is set as a boundary condition of the temperature field in the coupled transport model.

5. The method according to claim 1, characterized in that Determining the boundary conditions of the coupled transport model according to the atmospheric information and the basic physical information includes: The evaporation efficiency of the soil is calculated according to the atmospheric information and the basic physical information, and the evaporation rate is set as the atmospheric boundary condition of the coupled transport model.

6. The method according to claim 1, characterized in that The step of solving the coupled transport model by using the soil hydraulic parameters, the initial conditions and the boundary conditions includes: The coupled transport model is numerically simulated by using the soil hydraulic parameters, the initial conditions and the boundary conditions to obtain a temperature simulation value and a water content simulation value.

7. A simulation device for coupled migration of vadose zones, characterized in that: The device comprises: An acquisition module is used to acquire atmospheric information and basic physical information of soil in arid and semi-arid areas within a preset time; wherein the atmospheric information includes air temperature, air humidity, solar radiation and water vapor saturation; and the basic physical information includes soil moisture content, soil temperature, soil matrix potential and soil conductivity; A construction module is used to determine the water motion equation, the heat motion equation and the solute transport equation by using the basic physical information, and to construct a coupled transport model by using the water motion equation, the heat motion equation and the solute transport equation; A determination module, used to calculate and invert soil hydraulic parameters according to the basic physical information; and to determine the initial conditions and boundary conditions of the coupled transport model according to the atmospheric information and the basic physical information; The calculation module is used to solve the coupled transport model through the soil hydraulic parameters, the initial conditions and the boundary conditions, and simulate the coupled transport process of water, heat and solute in the vadose zone within a preset time.

8. An electronic device, characterized in that: The electronic device comprises: a memory and a processor, wherein the memory is used to store a computer program, and the processor is used to execute the computer program stored in the memory, so that the electronic device performs the steps of the method as claimed in any one of claims 1 to 6.

9. A computer-readable storage medium, characterized in that: The computer-readable storage medium stores a program, and when the program is run, it is used to implement the steps of the method according to any one of claims 1 to 6.

10. A computer program product, characterized in that The computer program product includes computer program codes, and when the computer program codes are executed on a computer, the computer is enabled to implement the steps of the method according to any one of claims 1 to 6.

Images