CN112161897A - Negative pressure boundary control method for large-scale soil water movement experiment system - Google Patents

Negative pressure boundary control method for large-scale soil water movement experiment system Download PDF

Info

Publication number
CN112161897A
CN112161897A CN202011031105.5A CN202011031105A CN112161897A CN 112161897 A CN112161897 A CN 112161897A CN 202011031105 A CN202011031105 A CN 202011031105A CN 112161897 A CN112161897 A CN 112161897A
Authority
CN
China
Prior art keywords
negative pressure
soil
boundary
vacuum
air chamber
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011031105.5A
Other languages
Chinese (zh)
Other versions
CN112161897B (en
Inventor
王康
李立
刘涛
李余隆
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Wuhan University WHU
Original Assignee
Wuhan University WHU
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Wuhan University WHU filed Critical Wuhan University WHU
Priority to CN202011031105.5A priority Critical patent/CN112161897B/en
Publication of CN112161897A publication Critical patent/CN112161897A/en
Application granted granted Critical
Publication of CN112161897B publication Critical patent/CN112161897B/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01NINVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
    • G01N13/00Investigating surface or boundary effects, e.g. wetting power; Investigating diffusion effects; Analysing materials by determining surface, boundary, or diffusion effects
    • G01N13/04Investigating osmotic effects
    • FMECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
    • F04POSITIVE - DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS FOR LIQUIDS OR ELASTIC FLUIDS
    • F04BPOSITIVE-DISPLACEMENT MACHINES FOR LIQUIDS; PUMPS
    • F04B37/00Pumps having pertinent characteristics not provided for in, or of interest apart from, groups F04B25/00 - F04B35/00
    • F04B37/10Pumps having pertinent characteristics not provided for in, or of interest apart from, groups F04B25/00 - F04B35/00 for special use
    • F04B37/14Pumps having pertinent characteristics not provided for in, or of interest apart from, groups F04B25/00 - F04B35/00 for special use to obtain high vacuum

Landscapes

  • Engineering & Computer Science (AREA)
  • Mechanical Engineering (AREA)
  • General Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Chemical & Material Sciences (AREA)
  • Analytical Chemistry (AREA)
  • Biochemistry (AREA)
  • General Health & Medical Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Immunology (AREA)
  • Pathology (AREA)
  • Investigation Of Foundation Soil And Reinforcement Of Foundation Soil By Compacting Or Drainage (AREA)

Abstract

The invention discloses a negative pressure boundary control method for a large-scale soil water movement experiment system, which comprises experiment soil and an experiment soil negative pressure measuring sensor; negative pressure boundary control device is in including setting up vacuum chamber, the setting of experiment soil below are in the pottery of vacuum chamber top contacts the membrane, negative pressure chamber vacuum can be adjusted, pottery contacts membrane connection experiment soil and vacuum chamber and is used for carrying out the negative pressure transmission by vacuum chamber to soil. Under the experimental condition of simulating the soil water movement by the control boundary, for the first class boundary condition (Dirichlet boundary) or the second class boundary condition (Newman boundary), flux is set based on the soil negative pressure measured by the soil negative pressure sensor buried in the boundary layer and the boundary layer, and the control of the first class or the second class boundary condition is realized by controlling and adjusting the vacuum air chamber negative pressure in the whole experimental period based on the functional relation of the vacuum air chamber negative pressure-boundary flux-soil negative pressure determined in the step one.

Description

Negative pressure boundary control method for large-scale soil water movement experiment system
Technical Field
The invention belongs to the technical field of soil-water boundary condition control, and particularly relates to a negative pressure boundary control method for a large-scale soil-water movement experiment system.
Background
Large soil water movement experiment systems such as an lysimeter, a slope flow experiment system including a soil seepage experiment and the like inevitably need to cut off soil by setting boundaries, boundary condition changes have obvious influence on soil water flow, and flux and processes measured by different boundaries (such as a seepage surface boundary, a gravity drainage boundary and a low-permeability water blocking boundary) show obvious differences. The boundary condition control is one of the unsolved core problems of the large-scale soil water movement experiment system for carrying out in-situ simulation experiments and control experiments.
The current common method for solving the boundary control problem is to enlarge the scale of the soil water movement experiment system, for example, the experiment system with large burial depth (more than 4m) is adopted, although the method reduces the influence of the boundary condition on the measurement result to a certain extent, the size of the corresponding experiment system is very large, and very high requirements are provided for monitoring the soil water movement and the accompanying and accompanying processes thereof. More importantly, the problem is not solved fundamentally in mechanism.
In addition, large soil water movement experimental systems typically have boundary conditions set to either watertight boundaries or positive pressure boundaries. The unsaturated soil water movement is divided into bounded and unbounded conditions, and the bounded unsaturated soil water movement is determined when the boundary is a saturated soil layer (underground water) or an impermeable interlayer; when the unsaturated zone is infinitely deep (or the groundwater burial depth is quite large) and the supply of soil water in the unsaturated zone to groundwater (or the supply of groundwater to the aeration zone) is not researched, the research is often focused on the soil water movement and the accompanying and accompanying processes within a certain depth range, for the unbounded soil water movement problem, an experimental area needs to be selected within a soil layer range with a certain depth, and the boundary is a negative pressure boundary in the case. However, all the existing large-scale soil water experimental devices cannot meet the requirements for developing the experiment for controlling the negative pressure boundary.
Disclosure of Invention
Aiming at the problems in the background art, the invention provides a negative pressure boundary control method for a large-scale soil water movement experiment system, establishes a method for realizing the negative pressure boundary control of the large-scale soil water movement experiment system through air chamber vacuum negative pressure control, and solves the key theory and technical problems of the reliability improvement of the large-scale soil water movement experiment.
In order to solve the technical problems, the invention adopts the following technical scheme:
a negative pressure boundary control method for a large-scale soil water movement experiment system comprises a control device and a control method for carrying out an in-situ simulation experiment by using the control device and controlling a negative pressure boundary under the condition of a boundary simulation soil water movement experiment, wherein the control device comprises:
the experimental soil is used for developing an in-situ simulation soil water flow experiment of the large-scale soil water movement system;
the experimental soil negative pressure measuring sensor is buried at the boundary position of the experimental soil and is used for measuring the negative pressure at the boundary position of the experimental soil;
negative pressure boundary control device, it is used for realizing normal position soil negative pressure boundary reduction, including setting up vacuum chamber, the setting of experiment soil below are in the pottery of vacuum chamber top contacts the membrane, negative pressure chamber vacuum can be adjusted, pottery contacts membrane connection experiment soil and vacuum chamber and is used for carrying out the negative pressure transmission to soil by vacuum chamber.
The invention also provides a negative pressure boundary control method of the large-scale soil water movement experiment system, which utilizes the control device to carry out the in-situ simulation experiment and control the negative pressure boundary under the condition of the boundary simulation soil water movement experiment and comprises the following steps:
s1, establishing a relation of vacuum air chamber negative pressure-soil water flux through parameter calibration of a ceramic contact membrane;
s2, measuring the negative pressure at the in-situ boundary position through an in-situ soil negative pressure measuring sensor embedded at the in-situ soil boundary position;
s3, according to the negative pressure measurement result at the original position of the boundary of the soil in the step S2, adopting the negative pressure-boundary flux-negative pressure function relation of the vacuum air chamber determined in the step S1, and adjusting the negative pressure of the vacuum air chamber to enable the experimental soil negative pressure at the boundary position of the experimental system to be consistent with the in-situ lower boundary soil negative pressure;
s4, in the experiment process, continuously measuring the negative pressure change of the in-situ soil boundary through the in-situ soil negative pressure measuring sensor embedded in the in-situ soil boundary position, and dynamically regulating and controlling the simulation experiment boundary through the step S3 to meet the condition that the simulation experiment boundary condition is consistent with the in-situ condition;
s5, collecting and discharging a water body entering a vacuum air chamber from a soil boundary;
s6, under the experimental condition of simulating the soil water movement by the control boundary, for the first kind of boundary condition, namely Dirichlet boundary or the second kind of boundary condition, namely Newman boundary, on the basis of the soil negative pressure measured by the experimental soil negative pressure measuring sensor for burying the experimental soil boundary layer and the set flux of the boundary layer, the control of the first kind or second kind of boundary condition is realized by controlling and adjusting the negative pressure of the vacuum air chamber in the whole experimental period on the basis of the functional relation of the vacuum air chamber negative pressure-boundary flux-soil negative pressure determined in the step S1.
Further, under the conditions of the in-situ simulation experiment and the control boundary experiment, the relationship of vacuum chamber negative pressure, soil negative pressure and boundary water flux under the condition that the ceramic contact membrane established by negative pressure boundary regulation and control is subjected to negative pressure transmission is as follows:
Figure BDA0002703712920000031
wherein t is a time dimension, Q is a boundary water flux continuously flowing in the experimental soil-ceramic contact membrane-vacuum air chamber, and delta h is the negative pressure difference between the experimental soil and the vacuum air chamber;
Figure BDA0002703712920000032
wherein R reflects the relative permeability of the ceramic contact membrane, K (h) is the soil hydraulic conductivity, and K is a function of the water potential h under unsaturated conditionsp,lpRespectively the permeability coefficient of the ceramic contact membrane, namely the water flux and the thickness under unit negative pressure gradient; b is the area of the negative pressure regulating boundary surface, a1Are first order coefficients. Further, in step S1, the ceramic membrane layer has a permeability coefficient KpAnd a first order coefficient a1To a desired rateThe parameters are determined by controlling the vacuum air chamber to form different negative pressures and measuring the corresponding experimental soil boundary negative pressure and flux, namely the second type of boundary condition controls the flux Q set according to the boundary and the experimental soil negative pressure measured by the experimental soil negative pressure measuring sensor buried in the experimental soil, determines the vacuum air chamber negative pressure and regulates and controls, realizes the boundary constant flux control in the whole experimental period, and fits the ceramic contact membrane permeability coefficient K according to the test resultpAnd a first order coefficient a1
In step S1, in the continuum of soil-ceramic contact membrane-vacuum chamber, the negative pressure generated by the vacuum chamber is transmitted to the soil through the ceramic contact membrane, forming a continuous negative pressure transmission process, in which the soil water motion equation is:
Figure BDA0002703712920000033
at the boundary position, unsaturated flux formed in the soil penetrates through the ceramic contact membrane to form continuous flow, and under balanced conditions, the flow equation in the ceramic contact membrane can be expressed as follows:
Figure BDA0002703712920000034
or
Figure BDA0002703712920000035
In the formula: z and t are the position and time dimensions, respectively, Kp,lpThe permeability coefficient and the thickness of the ceramic contact membrane are respectively; h, negative pressure; c and K are the soil water content and the hydraulic conductivity respectively and are a function of h, namely c (h) and K (h) are functions of water potential h under unsaturated conditions, and c (h) can be used as a constant (average value in a variation range of h) under the continuous flowing condition of the soil-ceramic contact membrane.
Under the condition of continuous flow, the Lappas transformation is carried out on the formula (1),
Figure BDA0002703712920000041
where p is the Laplace change constant, the general solution is:
Figure BDA0002703712920000042
wherein c is1And c2For the integration constant, Δ h is the negative pressure (water potential) difference between the soil boundary and the ceramic contact membrane, and the derivative of equation (3) on the position dimension z is obtained:
Figure BDA0002703712920000043
and (3) substituting the Laplace transform of the equation (2) into the equation (4) to obtain:
Figure BDA0002703712920000044
due to the interface position of the ceramic contact membrane and the air chamber
Figure BDA0002703712920000045
Therefore, the following equations (5) and (6) show that:
Figure BDA0002703712920000046
thus determining the integration constant c1=c2=c
Substituting the formula (7) into the formula (5) to obtain:
Figure BDA0002703712920000047
performing inverse Laplace transform on the formula (8), and obtaining the following result according to a residue theorem:
Figure BDA0002703712920000048
wherein, according to the relationship between the hyperbolic function and the complex variable trigonometric function:
Figure BDA0002703712920000049
chx-cosix, where x is a variable, i is an imaginary number,
Figure BDA00027037129200000410
order:
Figure BDA00027037129200000411
then (9) is:
Figure BDA0002703712920000051
then further according to the decomposition theorem, get:
Figure BDA0002703712920000052
wherein,
Figure BDA0002703712920000053
d is the diffusivity of soil water, and is substituted into the formula (11) to obtain:
Figure BDA0002703712920000054
wherein, negative pressure regulation and control boundary surface area is B, then soil boundary and ceramic contact membrane continuous water flux Q is:
Figure BDA0002703712920000055
that is, when the difference between the negative pressure at the soil boundary and the negative pressure at the air chamber is Δ h, the boundary flux Q at the time t is
Figure BDA0002703712920000056
When in use
Figure BDA0002703712920000057
The first term of the series can be taken (with an error of no more than 5%).
Formula (14) can be represented as
Figure BDA0002703712920000058
Thus, the boundary layer water flux Q generated under the condition that the soil boundary and the ceramic contact membrane continuously flow and the negative pressure difference between the soil and the air chamber is delta h is determined. (15) In the formula, because
Figure BDA0002703712920000059
The parameters comprise soil hydrodynamic parameters, including hydraulic conductivity K (h) and diffusivity D; ceramic contact membrane parameters including the permeability coefficient (water flux per negative pressure gradient) K of the ceramic contact membranepAnd a thickness lpAnd a first order coefficient a1. Wherein the soil hydrodynamic parameters can be directly measured by a conventional method, and the thickness l of the ceramic contact filmpDirectly measured by geometric method. Permeability coefficient K of ceramic membrane layerpAnd a first order coefficient a1Are parameters that need to be calibrated. As shown in fig. 3, the parameter calibration device performs parameter calibration by controlling the air chamber to form different negative pressures and measuring the corresponding soil boundary negative pressure and flux.
According to the technical scheme, the second type (Newman boundary) boundary control can determine the negative pressure of the vacuum air chamber according to the flux Q set by the set boundary and the soil negative pressure measured by the negative pressure sensor buried in the soil according to the (16) and regulate and control, and the constant flux control of the boundary is realized during the whole experiment. Wherein formula (16) is directly obtainable from formula (15):
Figure BDA0002703712920000061
under the experimental condition of simulating the soil water movement by the control boundary, for the first kind of boundary condition, namely Dirichlet boundary or the second kind of boundary condition, namely Newman boundary, based on the soil negative pressure measured by the experimental soil negative pressure measuring sensor for burying the experimental soil boundary layer and the set flux of the boundary layer, based on the functional relation of the vacuum air chamber negative pressure-boundary flux-soil negative pressure determined in the step S1, the control of the first kind or second kind of boundary condition is realized by controlling and adjusting the vacuum air chamber negative pressure in the whole experimental period.
Furthermore, under the condition that the experimental soil body is the viscous soil, an experimental soil negative pressure measuring sensor based on the soil vacuum degree galvanic couple induction as the measuring principle is adopted.
Further, ceramic contact membrane below is equipped with tipping bucket formula flow measuring device and soil water receiving flask in proper order among the vacuum air chamber, and it is used for collecting soil water among the experiment soil, be equipped with vacuum air chamber negative pressure measuring transducer in the vacuum air chamber and be used for measuring negative pressure in the vacuum air chamber, the vacuum air chamber lateral wall is equipped with the connection the connecting device of vacuum air chamber and vacuum pump regulates and control vacuum chamber internal vacuum through the vacuum pump and realizes the regulation and control of air chamber negative pressure, the vacuum pump is connected with negative pressure control regulator.
Compared with the prior art, the invention has the beneficial effects that:
the invention provides a method for realizing negative pressure boundary control of a large-scale soil water movement experiment system by air chamber vacuum negative pressure control, which is completely based on a soil negative pressure-flux-air chamber negative pressure function relation in a soil-ceramic contact membrane-vacuum air chamber continuous flow process, solves the key theory and technical problem of reliability improvement of the large-scale soil water movement experiment, and avoids the problem that high-requirement monitoring on soil water movement and accompanying processes thereof in the boundary control problem in the prior art by enlarging the scale of the soil water movement experiment system.
Drawings
FIG. 1 is a schematic view of a negative pressure boundary regulating device of a large soil water movement experiment system according to the present invention;
FIG. 2 is a negative pressure boundary regulating device of the large soil water movement experiment system of the invention;
FIG. 3 is a graphical representation of the parameter rate of the boundary negative pressure regulating device of an example of the present invention;
the system comprises 1-experimental soil, 2-negative pressure boundary control device, 3-experimental soil negative pressure measuring sensor, 4-negative pressure control regulator, 5-vacuum pump, 6-ceramic contact membrane, 7-tipping bucket type flow measuring device, 8-soil water collecting bottle, 9-vacuum air chamber, 10-vacuum air chamber negative pressure measuring sensor and 11-connecting device.
Detailed Description
The technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It is to be noted that the experimental methods described in the following embodiments are all conventional methods unless otherwise specified, and the reagents and materials, if not otherwise specified, are commercially available; in the description of the present invention, the terms "lateral", "longitudinal", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer", etc., indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, are only for convenience in describing the present invention and simplifying the description, and do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus, should not be construed as limiting the present invention.
Furthermore, the terms "horizontal", "vertical", "overhang" and the like do not imply that the components are required to be absolutely horizontal or overhang, but may be slightly inclined. For example, "horizontal" merely means that the direction is more horizontal than "vertical" and does not mean that the structure must be perfectly horizontal, but may be slightly inclined.
In the description of the present application, it is further noted that, unless expressly stated or limited otherwise, the terms "disposed," "mounted," "connected," and "connected" are to be construed broadly, e.g., as meaning either a fixed connection, a removable connection, or an integral connection; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meaning of the above terms in the present application can be understood in a specific case by those of ordinary skill in the art.
The present invention will be further explained with reference to the accompanying drawings and embodiments, where the embodiment provides a negative pressure boundary control method for a large soil water movement experiment system, including: the system comprises experimental soil, an experimental soil negative pressure measuring sensor and a negative pressure boundary control device, wherein the experimental soil is used for carrying out an in-situ simulation soil water flow experiment of a large-scale soil water movement system; the experimental soil negative pressure measuring sensor is buried at the experimental soil boundary position and used for measuring the negative pressure at the experimental soil boundary position; the edge pressing boundary control device comprises a vacuum air chamber arranged below the experimental soil and a ceramic contact membrane arranged at the top end of the vacuum air chamber, the vacuum degree of the negative pressure air chamber can be adjusted, and the ceramic contact membrane is connected with the experimental soil and the vacuum air chamber and is used for transmitting negative pressure to the soil from the vacuum air chamber.
Further preferred embodiment in the ceramic contact membrane below be equipped with tipping bucket formula flow measuring device and soil water receiving flask in proper order in the vacuum air chamber, it is used for collecting soil water in the experiment soil, be equipped with vacuum air chamber negative pressure measuring transducer in the vacuum air chamber and be used for measuring negative pressure in the vacuum air chamber, the vacuum air chamber lateral wall is equipped with the connection vacuum air chamber and vacuum pump's connecting device, adjusts and controls the regulation and control that vacuum chamber internal vacuum degree realized the air chamber negative pressure through the vacuum pump, the vacuum pump is connected with negative pressure control regulator.
When a large soil water movement system is adopted to carry out an in-situ simulation soil water flow experiment, the reduction of an in-situ boundary is realized through the following steps:
s1, establishing a relation of vacuum air chamber negative pressure-soil water flux through parameter calibration of a ceramic contact membrane;
s2, measuring the negative pressure at the in-situ boundary position through an in-situ soil negative pressure measuring sensor embedded at the in-situ soil boundary position;
s3, according to the negative pressure measurement result at the original position of the boundary of the soil in the step S2, adopting the negative pressure-boundary flux-negative pressure function relation of the vacuum air chamber determined in the step S1, and adjusting the negative pressure of the vacuum air chamber to enable the experimental soil negative pressure at the boundary position of the experimental system to be consistent with the in-situ lower boundary soil negative pressure;
s4, in the experiment process, continuously measuring the negative pressure change of the in-situ soil boundary through the in-situ soil negative pressure measuring sensor embedded in the in-situ soil boundary position, and dynamically regulating and controlling the simulation experiment boundary through the step S3 to meet the condition that the simulation experiment boundary condition is consistent with the in-situ condition;
s5, collecting and discharging a water body entering a vacuum air chamber from a soil boundary;
s6, under the experimental condition of simulating the soil water movement by the control boundary, for the first kind of boundary condition, namely Dirichlet boundary or the second kind of boundary condition, namely Newman boundary, on the basis of the soil negative pressure measured by the experimental soil negative pressure measuring sensor for burying the experimental soil boundary layer and the set flux of the boundary layer, the control of the first kind or second kind of boundary condition is realized by controlling and adjusting the negative pressure of the vacuum air chamber in the whole experimental period on the basis of the functional relation of the vacuum air chamber negative pressure-boundary flux-soil negative pressure determined in the step S1.
According to the technical scheme, the functional relation of the negative pressure of the vacuum air chamber, the negative pressure of the soil and the flow flux of the boundary water under the ceramic contact membrane negative pressure transfer condition is established in the step S1 as follows:
in the continuum formed by the soil-ceramic contact membrane-vacuum air chamber, the negative pressure formed by the vacuum air chamber is transmitted to the soil through the ceramic contact membrane to form a continuous negative pressure transmission process, and in the continuous process, the soil water motion equation is as follows:
Figure BDA0002703712920000091
at the boundary position, unsaturated flux formed in the soil penetrates through the ceramic contact membrane to form continuous flow, and under balanced conditions, the flow equation in the ceramic contact membrane can be expressed as follows:
Figure BDA0002703712920000092
or
Figure BDA0002703712920000093
In the formula: z and t are the position and time dimensions, respectively, Kp,lpThe permeability coefficient and the thickness of the ceramic contact membrane are respectively; h, negative pressure; c and K are the soil water content and the hydraulic conductivity respectively and are a function of h, namely c (h) and K (h) are functions of water potential h under unsaturated conditions, and c (h) can be used as a constant (average value in a variation range of h) under the continuous flowing condition of the soil-ceramic contact membrane.
Under the condition of continuous flow, the Lappas transformation is carried out on the formula (1),
Figure BDA0002703712920000094
where p is the Laplace change constant, the general solution is:
Figure BDA0002703712920000095
wherein c is1And c2For the integration constant, Δ h is the negative pressure (water potential) difference between the soil boundary and the ceramic contact membrane, and the derivative of equation (3) on the position dimension z is obtained:
Figure BDA0002703712920000096
and (3) substituting the Laplace transform of the equation (2) into the equation (4) to obtain:
Figure BDA0002703712920000097
due to the interface position of the ceramic contact membrane and the air chamber
Figure BDA0002703712920000101
Therefore, the following equations (5) and (6) show that:
Figure BDA0002703712920000102
thus determining the integration constant c1=c2=c
Substituting the formula (7) into the formula (5) to obtain:
Figure BDA0002703712920000103
performing inverse Laplace transform on the formula (8), and obtaining the following result according to a residue theorem:
Figure BDA0002703712920000104
wherein, according to the relationship between the hyperbolic function and the complex variable trigonometric function:
Figure BDA0002703712920000105
chx-cosix, where x is a variable, i is an imaginary number,
Figure BDA0002703712920000106
order:
Figure BDA0002703712920000107
then (9) is:
Figure BDA0002703712920000108
then further according to the decomposition theorem, get:
Figure BDA0002703712920000109
wherein,
Figure BDA00027037129200001010
d is the diffusivity of soil water, and is substituted into the formula (11) to obtain:
Figure BDA00027037129200001011
wherein, negative pressure regulation and control boundary surface area is B, then soil boundary and ceramic contact membrane continuous water flux Q is:
Figure BDA0002703712920000111
that is, when the difference between the negative pressure at the soil boundary and the negative pressure at the air chamber is Δ h, the boundary flux Q at the time t is
Figure BDA0002703712920000112
When in use
Figure BDA0002703712920000113
The first term of the series can be taken (with an error of no more than 5%).
Formula (14) can be represented as
Figure BDA0002703712920000114
In this way, the soil boundary and the ceramic are determinedAnd under the condition of continuous flow of the contact membrane, the water flux Q of the boundary layer occurs under the condition that the negative pressure difference between the soil and the air chamber is delta h. (15) In the formula, because
Figure BDA0002703712920000115
The parameters comprise soil hydrodynamic parameters, including hydraulic conductivity K (h) and diffusivity D; ceramic contact membrane parameters including the permeability coefficient (water flux per negative pressure gradient) K of the ceramic contact membranepAnd a thickness lpAnd a first order coefficient a1. Wherein the soil hydrodynamic parameters can be directly measured by a conventional method, and the thickness l of the ceramic contact filmpDirectly measured by geometric method. In this example, k (h) is 4.6 × 10 according to the experimental soil properties-6exp(0.0478h)m/s,D=2.18×10-2m2S; thickness l of ceramic contact filmp0.24mm, and the negative pressure regulating boundary surface area B of 78.5cm2. Permeability coefficient K of ceramic membrane layerpAnd a first order coefficient a1The parameters needing calibration are determined by setting calibration experiments under the condition of water potential difference deltah of different fluxes Q. As shown in fig. 3, the parameter calibration device performs parameter calibration by controlling the air chamber to form different negative pressures and measuring the corresponding soil boundary negative pressure and flux.
According to the technical scheme, the second type (Newman boundary) boundary control can determine the negative pressure of the vacuum air chamber according to the flux Q set by the set boundary and the soil negative pressure measured by the negative pressure sensor buried in the soil according to the (16) and regulate and control, and the constant flux control of the boundary is realized during the whole experiment. Wherein formula (16) is directly obtainable from formula (15):
Figure BDA0002703712920000116
under the experimental condition of simulating the soil water movement by controlling the boundary, for the first kind of boundary condition, namely Dirichlet boundary or the second kind of boundary condition, namely Newman boundary, based on the soil negative pressure measured by the experimental soil negative pressure measuring sensor for burying the experimental soil boundary layer and the boundary layer set flux, based on the functional relation of the vacuum air chamber negative pressure-boundary flux-soil negative pressure determined in the step S1, the control of the first kind or the second kind of boundary condition is realized by controlling and adjusting the vacuum air chamber negative pressure in the whole experimental period, as shown in FIG. 3, the relation of the vacuum air chamber negative pressure and the soil boundary negative pressure under the same condition of two groups of control Q fluxes is adopted, and the measurement results of all delta h and Q satisfy the control error of the formula (16)
Figure BDA0002703712920000121
Minimum, determine the permeability coefficient K of the ceramic contact membranepAnd first order coefficient a1Are respectively 7.87 multiplied by 10-8m/s and 0.0554.
According to the above technical scheme, the soil negative pressure sensor for measuring the soil negative pressure needs to have sufficient sensitivity, the soil negative pressure sensor using the potential energy balance principle as the measurement principle can generate larger control error due to the balance time process under the condition of cohesive soil quality, and the soil negative pressure sensor using the soil vacuum degree galvanic couple induction as the measurement principle, such as the soil negative pressure sensor
Figure BDA0002703712920000122
A sensor.
The foregoing examples are provided for illustration and description of the invention only and are not intended to limit the invention to the scope of the described examples. Furthermore, it will be understood by those skilled in the art that the present invention is not limited to the embodiments described above, and that many variations and modifications may be made in accordance with the teachings of the present invention, all of which fall within the scope of the invention as claimed.

Claims (5)

1. A negative pressure boundary control method for a large-scale soil water movement experiment system is characterized by comprising a control device and a control method for carrying out an in-situ simulation experiment by using the control device and controlling a negative pressure boundary under the condition of a boundary simulation soil water movement experiment, wherein the control device comprises:
the experimental soil is used for developing an in-situ simulation soil water flow experiment of the large-scale soil water movement system;
the experimental soil negative pressure measuring sensor is buried at the boundary position of the experimental soil and is used for measuring the negative pressure at the boundary position of the experimental soil;
the negative pressure boundary control device is used for realizing in-situ soil negative pressure boundary reduction and comprises a vacuum air chamber arranged below the experimental soil and a ceramic contact membrane arranged at the top end of the vacuum air chamber, wherein the vacuum degree of the negative pressure air chamber can be adjusted, and the ceramic contact membrane is connected with the experimental soil and the vacuum air chamber and is used for transmitting negative pressure from the vacuum air chamber to the soil;
the control method specifically comprises the following steps:
s1, establishing a relation of vacuum air chamber negative pressure-soil water flux through parameter calibration of a ceramic contact membrane;
s2, measuring the negative pressure at the in-situ boundary position through an in-situ soil negative pressure measuring sensor embedded at the in-situ soil boundary position;
s3, according to the negative pressure measurement result at the original position of the boundary of the soil in the step S2, adopting the negative pressure-boundary flux-negative pressure function relation of the vacuum air chamber determined in the step S1, and adjusting the negative pressure of the vacuum air chamber to enable the experimental soil negative pressure at the boundary position of the experimental system to be consistent with the in-situ lower boundary soil negative pressure;
s4, in the experiment process, continuously measuring the negative pressure change of the in-situ soil boundary through the in-situ soil negative pressure measuring sensor embedded in the in-situ soil boundary position, and dynamically regulating and controlling the simulation experiment boundary through the step S3 to meet the condition that the simulation experiment boundary condition is consistent with the in-situ condition;
s5, collecting and discharging a water body entering a vacuum air chamber from a soil boundary;
s6, under the experimental condition of simulating the soil water movement by the control boundary, for the first kind of boundary condition, namely Dirichlet boundary or the second kind of boundary condition, namely Newman boundary, on the basis of the soil negative pressure measured by the experimental soil negative pressure measuring sensor for burying the experimental soil boundary layer and the set flux of the boundary layer, the control of the first kind or second kind of boundary condition is realized by controlling and adjusting the negative pressure of the vacuum air chamber in the whole experimental period on the basis of the functional relation of the vacuum air chamber negative pressure-boundary flux-soil negative pressure determined in the step S1.
2. The negative pressure boundary control method for the large-scale soil water movement experiment system according to claim 1, wherein under the conditions of the in-situ simulation experiment and the control boundary experiment, the relationship of vacuum air chamber negative pressure-soil negative pressure-boundary water flux under the condition that the ceramic contact membrane established by negative pressure boundary regulation and control is subjected to negative pressure transmission is as follows:
Figure FDA0002703712910000021
wherein t is a time dimension, Q is a boundary water flux continuously flowing in the experimental soil-ceramic contact membrane-vacuum air chamber, and delta h is the negative pressure difference between the experimental soil and the vacuum air chamber;
Figure FDA0002703712910000022
wherein R reflects the relative permeability of the ceramic contact membrane, K (h) is the soil hydraulic conductivity, and K is a function of the water potential h under unsaturated conditionsp,lpRespectively the permeability coefficient of the ceramic contact membrane, namely the water flux and the thickness under unit negative pressure gradient; b is the area of the negative pressure regulating boundary surface, a1Are first order coefficients.
3. The method as claimed in claim 2, wherein in step S1, the ceramic membrane layer has a permeability coefficient KpAnd a first order coefficient a1The parameters to be calibrated are parameter calibration by a method of controlling the vacuum air chamber to form different negative pressures and measuring the corresponding experimental soil boundary negative pressure and flux, namely, the second type of boundary condition controls the flux Q set according to the boundary and the experimental soil negative pressure measured by the experimental soil negative pressure measuring sensor buried in the experimental soil, determines the vacuum air chamber negative pressure and regulates and controls, thereby realizing the constant flux control of the boundary during the whole experiment,fitting the permeability coefficient K of the ceramic contact membrane according to the test resultpAnd a first order coefficient a1
4. The method as claimed in claim 1, wherein a negative pressure measuring sensor for soil vacuum degree galvanic couple induction is used for measuring the negative pressure of the experimental soil under the condition that the experimental soil is viscous soil.
5. The negative pressure boundary control method of the large-scale soil water movement experiment system according to claim 1, characterized in that: ceramic contact membrane below is equipped with tipping bucket formula flow measuring device and soil water receiving flask in proper order among the vacuum chamber, and it is used for collecting soil water among the experiment soil, be equipped with vacuum chamber negative pressure measuring transducer in the vacuum chamber and be used for measuring negative pressure in the vacuum chamber, the vacuum chamber lateral wall is equipped with the connection the connecting device of vacuum chamber and vacuum pump, the regulation and control of vacuum chamber negative pressure is realized to vacuum chamber interior vacuum through vacuum pump regulation and control, the vacuum pump is connected with negative pressure control regulator.
CN202011031105.5A 2020-09-27 2020-09-27 Negative pressure boundary control method for large-scale soil water movement experiment system Expired - Fee Related CN112161897B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011031105.5A CN112161897B (en) 2020-09-27 2020-09-27 Negative pressure boundary control method for large-scale soil water movement experiment system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011031105.5A CN112161897B (en) 2020-09-27 2020-09-27 Negative pressure boundary control method for large-scale soil water movement experiment system

Publications (2)

Publication Number Publication Date
CN112161897A true CN112161897A (en) 2021-01-01
CN112161897B CN112161897B (en) 2021-10-22

Family

ID=73864280

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011031105.5A Expired - Fee Related CN112161897B (en) 2020-09-27 2020-09-27 Negative pressure boundary control method for large-scale soil water movement experiment system

Country Status (1)

Country Link
CN (1) CN112161897B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113945316A (en) * 2021-08-31 2022-01-18 河海大学 Negative pressure type under-film blind ditch and soil body air pressure distribution field test method

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102062773A (en) * 2010-11-05 2011-05-18 中国科学院新疆生态与地理研究所 Method for simulating soil solute transport in laboratory
US20120103077A1 (en) * 2010-10-29 2012-05-03 Solum, Inc. Microsampling Nutrient Measurement
CN107167576A (en) * 2017-05-09 2017-09-15 武汉大学 Unsaturated soil disorderly scattered kinetic coefficient assay method and device
CN107421868A (en) * 2017-07-19 2017-12-01 武汉大学 The assay method of unsaturated soil moisture diffusance based on the response of Complete heart block flow process
CN108801839A (en) * 2018-06-05 2018-11-13 武汉大学 Large-scale lysimeter intelligent monitor system
CN109270244A (en) * 2018-08-24 2019-01-25 河海大学 Solute transfer analogue observation device in soil under a kind of rainfall

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120103077A1 (en) * 2010-10-29 2012-05-03 Solum, Inc. Microsampling Nutrient Measurement
CN102062773A (en) * 2010-11-05 2011-05-18 中国科学院新疆生态与地理研究所 Method for simulating soil solute transport in laboratory
CN107167576A (en) * 2017-05-09 2017-09-15 武汉大学 Unsaturated soil disorderly scattered kinetic coefficient assay method and device
CN107421868A (en) * 2017-07-19 2017-12-01 武汉大学 The assay method of unsaturated soil moisture diffusance based on the response of Complete heart block flow process
CN108801839A (en) * 2018-06-05 2018-11-13 武汉大学 Large-scale lysimeter intelligent monitor system
CN109270244A (en) * 2018-08-24 2019-01-25 河海大学 Solute transfer analogue observation device in soil under a kind of rainfall

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
KANG WANG 等: "Characterizing heterogeneous soil water flow and solute transport using information measures", 《JOURNAL OF HYDROLOGY》 *
周云成 等: "变水头边界条件下土壤水分运动数值模拟", 《沈阳农业大学学报》 *
盛丰 等: "土壤非均匀水流运动与溶质运移的两区-两阶段模型", 《水利学报》 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113945316A (en) * 2021-08-31 2022-01-18 河海大学 Negative pressure type under-film blind ditch and soil body air pressure distribution field test method
CN113945316B (en) * 2021-08-31 2022-08-02 河海大学 Negative pressure type under-film blind ditch and soil body air pressure distribution field test method

Also Published As

Publication number Publication date
CN112161897B (en) 2021-10-22

Similar Documents

Publication Publication Date Title
Selker et al. Fingered flow in two dimensions: 1. Measurement of matric potential
US10712253B2 (en) Simulation device for interaction between deep reservoir rock and fluid in basin and method for using same
CN112161897B (en) Negative pressure boundary control method for large-scale soil water movement experiment system
JP5038516B2 (en) Groundwater surface runoff velocity flow rate measuring device and its monitoring device
EP3332227B1 (en) Mems flow sensor
CN105758495B (en) A kind of air pressure sensing liquidometer and its measurement method
JP2015517673A (en) Pressure transmitter with hydrogen getter
TW593727B (en) Combination differential and absolute pressure transducer for load lock control
CN110542537B (en) Experimental device for simulating influence of tides on underground water level and application method thereof
CN104458528A (en) Method and device for measuring field slope earth surface infiltration performance and permeability coefficient
Van Grinsven et al. Automated in situ measurement of unsaturated soil water flux
CN104535257B (en) A kind of silicon pressure drag temperature-compensating appraisal procedure
Perera et al. Surface-groundwater flow numerical model for barrier beach with exfiltration incorporated bottom boundary layer model
CN111855970A (en) Aeration zone moisture infiltration detection simulation device
Zhou et al. A grid-free, nonlinear shallow-water model with moving boundary
CN211270012U (en) Pressure regulating device and cooking utensil
CN109975350B (en) Method for determining heat conductivity coefficient of soil mineral
Lentz et al. An automated vacuum extraction control system for soil water percolation samplers
CN208223565U (en) A kind of big flow open channel measuring automatically flow system
Charuchittipan et al. Turbulent kinetic energy dissipation in the surface layer
CN105424543B (en) A kind of soil gas conductivity stable state measuring method
CN209148026U (en) A kind of pore characteristics experimental provision
CN206832309U (en) A kind of bailing test weir case of ultrasonic measuring automatically flow amount
Li et al. Calibration method for Manning's roughness coefficient for a river flume model
Chaudhry et al. A study of some unsteady free surface problems in porous media

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
CF01 Termination of patent right due to non-payment of annual fee
CF01 Termination of patent right due to non-payment of annual fee

Granted publication date: 20211022