CN113656756A - Method for calculating critical burial depth of boundary groundwater between oasis and transition zone of inland river arid region - Google Patents
Method for calculating critical burial depth of boundary groundwater between oasis and transition zone of inland river arid region Download PDFInfo
- Publication number
- CN113656756A CN113656756A CN202110985114.6A CN202110985114A CN113656756A CN 113656756 A CN113656756 A CN 113656756A CN 202110985114 A CN202110985114 A CN 202110985114A CN 113656756 A CN113656756 A CN 113656756A
- Authority
- CN
- China
- Prior art keywords
- diving
- oasis
- boundary
- water
- transition zone
- 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
Links
- 230000007704 transition Effects 0.000 title claims abstract description 89
- 238000009933 burial Methods 0.000 title claims abstract description 37
- 239000003673 groundwater Substances 0.000 title claims abstract description 36
- 238000000034 method Methods 0.000 title claims abstract description 36
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 151
- 230000009189 diving Effects 0.000 claims abstract description 118
- 239000002689 soil Substances 0.000 claims abstract description 69
- 230000009471 action Effects 0.000 claims abstract description 27
- 238000011160 research Methods 0.000 claims abstract description 22
- 230000008020 evaporation Effects 0.000 claims abstract description 14
- 238000001704 evaporation Methods 0.000 claims abstract description 14
- 238000001556 precipitation Methods 0.000 claims abstract description 12
- 238000004364 calculation method Methods 0.000 claims description 24
- 230000005484 gravity Effects 0.000 claims description 9
- 229920006395 saturated elastomer Polymers 0.000 claims description 9
- 238000005070 sampling Methods 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 5
- 238000013508 migration Methods 0.000 claims description 4
- 230000005012 migration Effects 0.000 claims description 4
- 239000011148 porous material Substances 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 3
- 239000002245 particle Substances 0.000 claims description 3
- 238000011161 development Methods 0.000 abstract description 6
- 230000007246 mechanism Effects 0.000 abstract description 5
- 230000008569 process Effects 0.000 abstract description 3
- 230000002262 irrigation Effects 0.000 description 10
- 238000003973 irrigation Methods 0.000 description 10
- 238000010586 diagram Methods 0.000 description 8
- 238000005273 aeration Methods 0.000 description 4
- 238000011835 investigation Methods 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000031068 symbiosis, encompassing mutualism through parasitism Effects 0.000 description 2
- 206010003694 Atrophy Diseases 0.000 description 1
- 102100021022 Gastrin Human genes 0.000 description 1
- 108010052343 Gastrins Proteins 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000037444 atrophy Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- AOXOCDRNSPFDPE-UKEONUMOSA-N chembl413654 Chemical compound C([C@H](C(=O)NCC(=O)N[C@H](CC=1C2=CC=CC=C2NC=1)C(=O)N[C@H](CCSC)C(=O)N[C@H](CC(O)=O)C(=O)N[C@H](CC=1C=CC=CC=1)C(N)=O)NC(=O)[C@@H](C)NC(=O)[C@@H](CCC(O)=O)NC(=O)[C@@H](CCC(O)=O)NC(=O)[C@@H](CCC(O)=O)NC(=O)[C@H](CCC(O)=O)NC(=O)[C@H](CCC(O)=O)NC(=O)[C@H](CC(C)C)NC(=O)[C@H](CC=1C2=CC=CC=C2NC=1)NC(=O)[C@H]1N(CCC1)C(=O)CNC(=O)[C@@H](N)CCC(O)=O)C1=CC=C(O)C=C1 AOXOCDRNSPFDPE-UKEONUMOSA-N 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 208000035475 disorder Diseases 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000013589 supplement Substances 0.000 description 1
- 238000012795 verification Methods 0.000 description 1
- 239000011800 void material Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
- G06F17/13—Differential equations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/15—Correlation function computation including computation of convolution operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Mathematical Optimization (AREA)
- Pure & Applied Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Theoretical Computer Science (AREA)
- General Engineering & Computer Science (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Computing Systems (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Evolutionary Biology (AREA)
- Probability & Statistics with Applications (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to a method for calculating critical buried depth of underground water at boundary of oasis and transition zones of inland river arid regions, which comprises the following steps: determining soil type and precipitation related parameters and the like; determining the boundary moisture driving conditions of oasis and transition zones in the arid regions of inland rivers; determining a theoretical formula of an internal water distribution curve of the diving affected layer; calculating the thickness of the diving influence layer; determining a water distribution curve of a submerged evaporation affected zone in a research area; determining the thickness of the usable layer of the boundary vegetation of the oasis and the transition zone; and determining the critical buried depth of the boundary groundwater of the oasis and the transition zone. The method for calculating the critical burial depth of the oasis and transition zone boundary underground water in the inland river arid region has clear physical process and action mechanism, has universality for determining a suitable development region and scale control of the oasis, and provides reliable theoretical and technical support for guiding reasonable development and utilization of underground water resources in the inland river arid region and reasonable control of the oasis.
Description
Technical Field
The invention relates to an ecological and environment-friendly hydrological and water conservancy calculation method, in particular to a calculation method for critical buried depth of underground water at the boundary of oasis and transition zones of inland rivers in arid regions.
Background
The method has the advantages that the reasonable critical burial depth of the groundwater at the boundary of the oasis and the transition zone of the inland river is determined, the significance for determining the boundary line of the oasis and the transition zone is great, and the technical difficulty for determining the appropriate development scale of the oasis in the arid zone from the perspective of natural hydrological conditions needs to be broken through urgently.
The mountainous area rainfall runoff production flow in the arid area of the inland river in the northwest area, the runoff consumption in the plain area, the underground water submerged flow field formed by the runoff at the mountain outlet supports the plain area to form an ecological pattern of oasis-transition zone-desert symbiosis, and the transition zone is a natural ecological barrier of the oasis. On one hand, the area of the transition zone is directly occupied by a large amount of irrigation regions in the oasis in disorder expansion, on the other hand, irrigation causes water resources to be concentrated towards the interior of the oasis, the water level of underground water is reduced, vegetation of the transition zone is degraded, and the desert approaches the oasis. Research shows that the total artificial oasis of the northwest region from 70 to 90 years is increased by 1.53 km26800km of natural oasis atrophy2The area of the staggered transition zone of the oasis desert is reduced by 4.4 km2The desertification area is enlarged by 3.6 km2。
At present, aiming at the problems that the oasis and the transition zone are lack of clear definition and are difficult to distinguish in practice, although expert scholars divide the oasis and the transition zone, the critical burial depth at the boundary of the oasis and the transition zone is not quantitatively calculated. The disordered expansion of the oasis can be effectively controlled by researching critical hydrological conditions of the boundaries of the oasis and the transition zone, and a reasonable area is kept, so that the ecological safety of the transition zone of the natural barrier of the oasis is ensured.
Disclosure of Invention
The invention aims to solve the technical problem that the oasis and transition zones are lack of clear definition at present and are difficult to distinguish in practice, and provides a method for calculating the critical burial depth of the groundwater at the boundary of the oasis and the transition zones in the inland river arid regions, which comprises the following steps:
the method comprises the following steps: sampling soil water in an experimental area at multiple points, and collecting soil and meteorological related parameters;
step two: determining the boundary moisture driving conditions of oasis and transition zones in the arid regions of inland rivers;
step three: determining an internal water distribution curve of a diving influence layer formed by diving evaporation in an arid region of an inland river;
step four: calculating the thickness h of the diving influence layer0;
Step five: determining parameters of a water distribution curve of a submerged evaporation affected zone in an experimental area;
step six: determining the thickness m of the overlapping part of the oasis and the transition zone boundary diving influence layer and the root action layer;
step seven: and determining the critical buried depth of the boundary groundwater of the oasis and the transition zone.
Wherein, soil and meteorological relevant parameter in step one includes: and (4) carrying out multipoint sampling and averaging to obtain the soil moisture content, the soil particle size, the soil porosity and the precipitation parameters above the submerging surface.
And the determination method of the boundary water driving condition of the oasis and the transition zone in the inland river arid region in the second step comprises the following steps: according to the hydrologic cycle characteristics of the arid region of the inland river, analyzing the ecological circle structure under the driving of both rainfall and diving in the plain region, wherein the relation between the rainfall supply P and the diving supply Q in the ecological water demand of the oasis vegetation is that P/Q is less than or equal to 1; the relationship between the rainfall supply P and the diving supply Q in the ecological water demand of the vegetation in the transition zone is that P/Q is more than or equal to 1; therefore, the critical hydrological condition of the ecological water demand of the boundary vegetation of the oasis and the transition zone is determined to be that the rainfall replenishment quantity P is equal to the diving replenishment quantity Q, and the calculation formula is as follows:
P=Q。
the method for determining the distribution curve in the third step comprises the following steps: the driving force of water migration in the diving influence layer is capillary force, the impedance is gravity, and a characteristic function omega (h) of soil capillary water movement conforms to an inverted S-shaped curve and is expressed by the following formula:
wherein h isD=h0/2;
In the formula: omega is the water content of the soil inside the diving affected layer; omega0The initial water content of the soil; omegasThe saturated water content of the soil is obtained; h is the height from any place of the diving influence layer to the diving surface; alpha is the capillary force intensity coefficient and is related to the soil grain diameter and the pore characteristics; h is0The thickness of the diving influence layer; h isDThe height from the turning point of the water distribution curve of the diving affected layer to the diving surface.
Wherein, the thickness h of the diving affected layer in the fourth step0The calculation method comprises the following steps: the diving affected zone is formed by diving under the action of capillary force, so that the theoretical formula of the maximum capillary water rise height can be used for calculation:
in the formula: t is the temperature of soil water; σ is the surface tension of the soil water, ρ is the density of the soil water, g is the acceleration of gravity.
And in the fifth step, the parameter determination method for the water distribution curve of the submerged evaporation affected zone in the experimental area comprises the following steps: the measured soil water parameters are utilized to solve the parameters of the water distribution curve of the diving affected layer in the research area, and the initial water content of the research area is 5 percent, the saturated water content is 38 percent, the diving affected layer is 1.29m, and the capillary force intensity coefficient is 3.168, so that the water distribution curve of the diving affected layer in the research area is determined as follows:
and in the sixth step, the method for determining the thickness m of the overlapping part of the oasis and transition zone boundary diving influence layer and the root action layer comprises the following steps: and integrating the water distribution curve of the diving influence layer to obtain the water quantity of the diving supply vegetation, and determining the thickness m of the overlapping part of the oasis and the transition zone boundary diving influence layer and the root action layer by using the known precipitation quantity P in a reverse calculation way according to the following calculation formula by utilizing the boundary hydrological condition:
and in the seventh step, the method for determining the critical burial depth of the boundary groundwater of the oasis and the transition zone comprises the following steps: the thickness of the overlapping part of the thickness of the diving influence layer and the thickness of the heel system action layer subtracted by the sum of the thickness of the diving influence layer and the thickness of the heel system action layer is the critical burial depth h of the oasis and the boundary underground water of the transition zone1:
h1=h0+D-m;
In the formula: h is1The critical buried depth of the boundary groundwater of the oasis and the transition zone; h is0The thickness of the diving influence layer; d is the thickness of the root action layer.
The implementation of the invention has the following beneficial effects: according to the invention, by analyzing a transition zone moisture driving mechanism and a submerged evaporation surface vegetation replenishment mechanism, the ecological water demand relationship of vegetation at the boundary of the oasis and the transition zone is determined, and the critical burial depth of underground water at the boundary of the oasis and the transition zone in the inland river arid region is defined; deducing a soil water distribution curve formed in the aeration zone by the submerged evaporation, and determining the water quantity of the groundwater absorbed and utilized by the vegetation at the boundary of the oasis and the transition zone; determining the critical burial depth of the groundwater in the oasis and the transition zone by utilizing the ecological water requirement source relation in the transition zone water driving mechanism; the method has clear physical process and action mechanism, has universality for determining a suitable development area and scale control of the oasis, and provides reliable theoretical and technical support for guiding reasonable development and utilization of underground water resources in arid areas of inland rivers and reasonable control of the oasis.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, it is obvious that the drawings in the following description are only some embodiments of the present invention, and for those skilled in the art, other drawings can be obtained according to the drawings without creative efforts.
FIG. 1 is a schematic diagram of the boundary water driving conditions of oasis and transition zones in the inland river arid regions;
FIG. 2 is a water distribution diagram of a diving affected zone;
FIG. 3 is a schematic diagram illustrating calculation of critical burial depth of groundwater at the boundary between oasis and transition zones;
FIG. 4 is a schematic flow chart of the method for calculating critical burial depth of groundwater at the boundary of oasis and transition zones of inland river arid regions;
FIG. 5 is a schematic diagram of a groundwater burial depth survey point in a Roche irrigation area.
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 present application, 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.
The invention provides a method for calculating the critical burial depth of underground water at the boundary of oasis and transition zones of inland river arid regions, which has the following principle:
and (3) the mountainous area rainfall runoff production flow in the arid area of the inland river in the northwest region, the runoff consumption in the plain area, and the groundwater submerged field formed by the runoff at the mountain outlet supports the plain area to form an ecological pattern of oasis-transition zone-desert symbiosis. The ecological water requirement of oasis vegetation in the plain area mainly depends on the diving supply amount, and the rainfall supply amount is auxiliary; the ecological water demand of the transition zone vegetation mainly depends on the rainfall supply amount, and the diving supply amount is indispensable. Since precipitation can be directly obtained by observation, it is mainly the diving supplement that is the difficulty of calculation.
The diving in the arid region does not directly act on the ground surface vegetation, but forms soil water in the aeration zone through diving evaporation, and the soil water is absorbed and utilized by the ground surface vegetation, so the distribution rule of the soil water formed by the diving evaporation in the aeration zone and the characteristic that the root system of the vegetation absorbs water are researched, the water quantity of the diving replenishing ground surface vegetation can be quantitatively determined, and the critical burial depth of underground water at the boundary of the oasis and the transition zone is determined, which is the key for determining the essential difference between the oasis and the transition zone and is also the key for guiding the reasonable development scale of the oasis.
The method for calculating the critical burial depth of the boundary groundwater between the oasis and the transition zone of the inland river arid region comprises the following steps:
please refer to fig. 1, fig. 1 is a schematic diagram of the boundary moisture driving conditions of oasis and transition zones in the inland river arid regions. According to the ecological circle structure theory of the arid region of the inland river, the relationship between the rainfall supply P and the diving supply Q in the ecological water demand of the oasis vegetation is that P/Q is less than or equal to 1; the relationship between the rainfall supply P and the diving supply Q in the ecological water demand of the vegetation in the transition zone is that P/Q is more than or equal to 1. Therefore, the critical hydrological condition of the ecological water demand of the boundary vegetation of the oasis and the transition zone is determined to be that the rainfall replenishment quantity P is equal to the diving replenishment quantity Q, and the calculation formula is as follows:
P=Q;
referring to fig. 2, fig. 2 is a water distribution diagram of the diving affected layer. The mechanical characteristics of distribution of soil water in an aeration zone formed by diving evaporation are that the driving force of water migration is capillary force, the impedance is gravity, a characteristic function omega (h) of a water distribution curve in a diving influence layer accords with the expression of an available S curve, and the formula expression is as follows:
wherein h isD=h0/2;
In the formula: omega is the water content of the soil inside the diving affected layer; omega0The initial water content of the soil; omegasThe saturated water content of the soil is obtained; h is the height from any place of the diving influence layer to the diving surface, and the maximum range height h is equal to the maximum rising height of the capillary water; alpha is the capillary force intensity coefficient and is related to the soil grain diameter and the pore characteristics; h is0The thickness of the diving influence layer; h isDThe height from the turning point of the water distribution curve of the diving affected layer to the diving surface.
Calculating the thickness h of the diving influence layer0The steps of (1): the formation of the diving influence layer is mainly under the action of capillary force when diving, so the theoretical formula of the maximum capillary water rise height can be used for calculating:
in the formula: t is the temperature of soil water; σ is the surface tension of the soil water, ρ is the density of the soil water, g is the acceleration of gravity.
The measured soil water parameters are utilized to solve the parameters of the water distribution curve of the diving affected layer in the research area, and the initial water content of the research area is 5 percent, the saturated water content is 38 percent, the diving affected layer is 1.29m, and the capillary force intensity coefficient is 3.168, so that the water distribution curve of the diving affected layer in the research area is determined as follows:
referring to fig. 3, fig. 3 is a schematic diagram illustrating calculation of critical burial depth of groundwater at the boundary between oasis and transition zones. And integrating the water distribution curve of the diving influence layer to obtain the water quantity of the diving supply vegetation, and determining the thickness m of the overlapping part of the oasis and the transition zone boundary diving influence layer and the root action layer by using the known precipitation quantity P in a reverse calculation way according to the following calculation formula by utilizing the boundary hydrological condition:
subtracting the thickness of the overlapped part of the diving influence layer thickness and the heel system action layer thickness from the sum of the diving influence layer thickness and the heel system action layer thickness, namely the critical burial depth of the oasis and the transition zone boundary underground water:
h1=h0+D-m;
in the formula h1The critical buried depth of the boundary groundwater of the oasis and the transition zone; h is0The thickness of the diving influence layer; d is the thickness of the root action layer.
Referring to fig. 4, fig. 4 is a schematic flow chart of the method for calculating the critical burial depth of groundwater at the boundary between oasis and transition zones of inland river arid regions. The method for calculating the critical burial depth of the boundary groundwater of the oasis and the transition zone of the inland river arid region comprises the following steps of:
the method comprises the following steps: carrying out multi-point sampling on soil water in an experimental area, and collecting soil and meteorological related parameters: and (4) carrying out multipoint sampling and averaging to obtain the soil moisture content, the soil particle size, the soil porosity and the precipitation parameters above the submerging surface.
Step two: determining the boundary moisture driving conditions of oasis and transition zones of inland river arid regions: according to the hydrologic cycle characteristics of the arid region of the inland river, analyzing the ecological circle structure under the driving of both rainfall and diving in the plain region, wherein the relation between the rainfall supply P and the diving supply Q in the ecological water demand of the oasis vegetation is that P/Q is less than or equal to 1; the relationship between the rainfall supply P and the diving supply Q in the ecological water demand of the vegetation in the transition zone is that P/Q is more than or equal to 1; therefore, the critical hydrological condition of the ecological water demand of the boundary vegetation of the oasis and the transition zone is determined to be that the rainfall replenishment quantity P is equal to the diving replenishment quantity Q, and the calculation formula is as follows:
P=Q。
step three: determining the inner water distribution curve of a diving influence layer formed by diving evaporation in an arid region of an inland river: the driving force of water migration in the diving influence layer is capillary force, the impedance is gravity, and a characteristic function omega (h) of soil capillary water movement conforms to an inverted S-shaped curve and is expressed by the following formula:
wherein h isD=h0/2;
In the formula: omega is the water content of the soil inside the diving affected layer; omega0The initial water content of the soil; omegasThe saturated water content of the soil is obtained; h is the height from any place of the diving influence layer to the diving surface; alpha is the capillary force intensity coefficient and is related to the soil grain diameter and the pore characteristics; h is0The thickness of the diving influence layer; h isDThe height from the turning point of the water distribution curve of the diving affected layer to the diving surface.
Step four: calculating the thickness h of the diving influence layer0: the diving affected zone is formed by diving under the action of capillary force, so that the theoretical formula of the maximum capillary water rise height can be used for calculation:
in the formula: t is the temperature of soil water; σ is the surface tension of the soil water, ρ is the density of the soil water, g is the acceleration of gravity.
Step five: determining parameters for obtaining a water distribution curve of a submerged evaporation affected zone in an experimental area: the measured soil water parameters are utilized to solve the parameters of the water distribution curve of the diving affected layer in the research area, and the initial water content of the research area is 5 percent, the saturated water content is 38 percent, the diving affected layer is 1.29m, and the capillary force intensity coefficient is 3.168, so that the water distribution curve of the diving affected layer in the research area is determined as follows:
step six: determining the thickness m of the overlapping part of the oasis and transition zone boundary diving influence layer and the root action layer: and integrating the water distribution curve of the diving influence layer to obtain the water quantity of the diving supply vegetation, and determining the thickness m of the overlapping part of the oasis and the transition zone boundary diving influence layer and the root action layer by using the known precipitation quantity P in a reverse calculation way according to the following calculation formula by utilizing the boundary hydrological condition:
step seven: determining the critical buried depth of the boundary groundwater of the oasis and the transition zone: the thickness of the overlapping part of the thickness of the diving influence layer and the thickness of the heel system action layer subtracted by the sum of the thickness of the diving influence layer and the thickness of the heel system action layer is the critical burial depth h of the oasis and the boundary underground water of the transition zone1:
h1=h0+D-m;
In the formula: h is1The critical buried depth of the boundary groundwater of the oasis and the transition zone; h is0The thickness of the diving influence layer; d is the thickness of the root action layer.
The specific application embodiment is as follows:
and selecting a Luo city irrigation area in Gastrin county of Zhangye city with serious salinization at the lower reaches in the black river basin as a case area for research. And calculating the critical burial depth of the groundwater at the boundary of the oasis and the transition zone of the Lucheng irrigated area by adopting a constructed theoretical formula, and verifying and analyzing the calculation result by combining with the actual observation data in the field of the groundwater.
According to a classification system adopted by national soil species, statistical analysis is carried out on the soil in different regions of a research area to obtain that the main soil type of the Roche irrigation area is aeolian sandy soil, the aeolian sandy soil type is mainly loamy sandy soil, the effective grain diameter d is 0.1mm, and the void degree n is 42%. And combining with precipitation information statistics of national weather stations, obtaining the average precipitation of 78mm in 1956-2018 years in the research area.
According to the sampling condition of 33 points of field sampling, the temperature range of the soil before the irrigation period in the research area is 10-30 ℃, and the minimum temperature of the soil is 10 ℃. Studies have shown that the surface tension of water is primarily dependent on temperatureThe influence is in inverse proportion in a certain range. Taking the surface tension of 74 multiplied by 10 corresponding to the soil temperature of 10 DEG C-3N/m, the maximum value of the surface tension of water at 10 ℃. The density of water is 1 x 103kg/m3The change in density of water is small. Therefore, neglecting the density change of water with temperature, the gravity acceleration g takes 9.8N/kg.
substituting the related parameters into a formula, and calculating to obtain the salinization critical underground water level of 1.29 m.
The measured soil water and the like are utilized to solve the parameters of the water distribution curve of the diving affected layer in the research area, and the initial water content of the research area is 5 percent, the saturated water content is 38 percent, the diving affected layer is 1.29m, and the capillary force intensity coefficient is 3.168, so that the water distribution curve of the diving affected layer in the research area is determined as follows:
when the precipitation is 78mm, the corresponding m is 0.81, h0The root system D of the vegetation passing through the boundary of the wild oasis and the transition zone is 4-5 m, the average value is 4.5m, and the critical buried depth of underground water at the boundary of the oasis and the transition zone is 4.98 m.
Field verification based on actual measurement underground water burial depth data:
the area of field investigation is about 600km with Luocheng irrigation area as main body2The investigation points are distributed on both sides of the main flow of the black river.
By investigating the groundwater burial depths of different regions of an oasis and a transition zone of a research region, 68 point data are obtained in total, please refer to fig. 5, and fig. 5 is a schematic diagram of groundwater burial depth investigation points of a levy irrigation region.
Through analysis, the buried depths of underground water in regions with serious salinization in oasis of a Luo city irrigation area are all within 1m, and the buried depths of underground water in regions with weak salinization are 1-2 m and are mainly distributed near river main flows and in cultivated lands; the buried water level of underground water in oasis of a Luo city irrigation area is about 2-6 m; the underground water burial depth of the transition zone area at the periphery of the oasis is about 6-13 m, and the underground water burial depths corresponding to different vegetation communities are different.
The buried depth of the oasis and the underground water in the area is analyzed to be between 2 and 6m, the buried depth of the underground water in the transition zone area is 6 to 13m, the buried depth of the underground water at the boundary of the oasis and the transition zone is about 5 to 6m, and the buried depth basically accords with the result of 4.98m calculated by a deduced theoretical formula, so that the calculation formula is proved to have certain reliability.
Finally, it should be noted that the above is only for illustrating the technical solution of the present invention and not for limiting, although the present invention has been described in detail with reference to the preferred arrangement, it should be understood by those skilled in the art that the technical solution of the present invention (such as various terrains and soil types, derivation process of application of formula, sequence of steps, etc.) may be modified or equivalently replaced without departing from the spirit and scope of the technical solution of the present invention.
Claims (8)
1. A method for calculating the critical burial depth of groundwater at the boundary of an oasis and a transition zone of an inland river arid region is characterized by comprising the following steps of:
the method comprises the following steps: sampling soil water in an experimental area at multiple points, and collecting soil and meteorological related parameters;
step two: determining the boundary moisture driving conditions of oasis and transition zones in the arid regions of inland rivers;
step three: determining an internal water distribution curve of a diving influence layer formed by diving evaporation in an arid region of an inland river;
step four: calculating the thickness h of the diving influence layer0;
Step five: determining parameters of a water distribution curve of a submerged evaporation affected zone in an experimental area;
step six: determining the thickness m of the overlapping part of the oasis and the transition zone boundary diving influence layer and the root action layer;
step seven: and determining the critical buried depth of the boundary groundwater of the oasis and the transition zone.
2. The method for calculating the critical burial depth of the inland river arid area oasis and transition zone boundary groundwater according to claim 1, wherein the soil and weather-related parameters in the first step comprise: and (4) carrying out multipoint sampling and averaging to obtain the soil moisture content, the soil particle size, the soil porosity and the precipitation parameters above the submerging surface.
3. The method for calculating the critical burial depth of the groundwater at the boundary of the oasis and the transition zone of the inland river arid region according to claim 2, wherein the method for determining the driving condition of the water at the boundary of the oasis and the transition zone of the inland river arid region in the second step comprises the following steps: according to the hydrologic cycle characteristics of the arid region of the inland river, analyzing the ecological circle structure under the driving of both rainfall and diving in the plain region, wherein the relation between the rainfall supply P and the diving supply Q in the ecological water demand of the oasis vegetation is that P/Q is less than or equal to 1; the relationship between the rainfall supply P and the diving supply Q in the ecological water demand of the vegetation in the transition zone is that P/Q is more than or equal to 1; therefore, the critical hydrological condition of the ecological water demand of the boundary vegetation of the oasis and the transition zone is determined to be that the rainfall replenishment quantity P is equal to the diving replenishment quantity Q, and the calculation formula is as follows:
P=Q。
4. the method for calculating the critical burial depth of the inland river arid area oasis and transition zone boundary groundwater according to claim 3, wherein the distribution curve determination method in the third step is as follows: the driving force of water migration in the diving influence layer is capillary force, the impedance is gravity, and a characteristic function omega (h) of soil capillary water movement conforms to an inverted S-shaped curve and is expressed by the following formula:
hD=h0/2;
in the formula: omega is the water content of the soil inside the diving affected layer; omega0The initial water content of the soil; omegasThe saturated water content of the soil is obtained; h is influence of divingThe height of any place of the layer from the diving surface; alpha is the capillary force intensity coefficient and is related to the soil grain diameter and the pore characteristics; h is0The thickness of the diving influence layer; h isDThe height from the turning point of the water distribution curve of the diving affected layer to the diving surface.
5. The method for calculating the critical burial depth of the groundwater at the boundary of the oasis and the transition zone of the inland river arid region according to claim 4, wherein the thickness h of the diving influence layer in the fourth step is0The calculation method comprises the following steps: the diving affected zone is formed by diving under the action of capillary force, so that the theoretical formula of the maximum capillary water rise height can be used for calculation:
in the formula: t is the temperature of soil water; σ is the surface tension of the soil water, ρ is the density of the soil water, g is the acceleration of gravity.
6. The method for calculating the critical burial depth of the boundary groundwater between the oasis and the transition zone in the inland river arid region according to claim 5, wherein in the fifth step, the parameter determination method for the water distribution curve of the submerged evaporation affected zone in the experimental region comprises the following steps: the measured soil water parameters are utilized to solve the parameters of the water distribution curve of the diving affected layer in the research area, and the initial water content of the research area is 5 percent, the saturated water content is 38 percent, the diving affected layer is 1.29m, and the capillary force intensity coefficient is 3.168, so that the water distribution curve of the diving affected layer in the research area is determined as follows:
7. the method for calculating the critical burial depth of the inland river arid area oasis and transition zone boundary groundwater according to claim 6, wherein the method for determining the thickness m of the overlapping part of the oasis and transition zone boundary diving influence layer and the root action layer in the sixth step is as follows: and integrating the water distribution curve of the diving influence layer to obtain the water quantity of the diving supply vegetation, and determining the thickness m of the overlapping part of the oasis and the transition zone boundary diving influence layer and the root action layer by using the known precipitation quantity P in a reverse calculation way according to the following calculation formula by utilizing the boundary hydrological condition:
8. the method for calculating the critical burial depth of the oasis and the transition zone boundary groundwater in the inland river arid region according to claim 7, wherein the method for determining the critical burial depth of the oasis and the transition zone boundary groundwater in the step seven comprises the following steps: the thickness of the overlapping part of the thickness of the diving influence layer and the thickness of the heel system action layer subtracted by the sum of the thickness of the diving influence layer and the thickness of the heel system action layer is the critical burial depth h of the oasis and the boundary underground water of the transition zone1:
h1=h0+D-m;
In the formula: h is1The critical buried depth of the boundary groundwater of the oasis and the transition zone; h is0The thickness of the diving influence layer; d is the thickness of the root action layer.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110985114.6A CN113656756B (en) | 2021-08-26 | Method for calculating critical buried depth of groundwater at boundary of oasis and transition zone in arid region of inland river |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110985114.6A CN113656756B (en) | 2021-08-26 | Method for calculating critical buried depth of groundwater at boundary of oasis and transition zone in arid region of inland river |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113656756A true CN113656756A (en) | 2021-11-16 |
CN113656756B CN113656756B (en) | 2024-06-04 |
Family
ID=
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115658822A (en) * | 2022-09-23 | 2023-01-31 | 中国水利水电科学研究院 | Computing method for determining reasonable ecological pattern of farming-pasturing staggered zone in semiarid region |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020144953A1 (en) * | 1995-05-05 | 2002-10-10 | Kerfoot William B. | Gas-gas-water treatment system for groundwater and soil remediation |
CN101385433A (en) * | 2008-09-19 | 2009-03-18 | 中国科学院新疆生态与地理研究所 | Groundwater buried depth abductive technique of desert plant seedling root system growth |
JP2012044913A (en) * | 2010-08-26 | 2012-03-08 | Nishizawa:Kk | Leaf vegetable harvester |
CN103093114A (en) * | 2013-02-05 | 2013-05-08 | 河海大学 | Distributed-type river basin water deficit calculating method based on terrain and soil characteristics |
CN104094693A (en) * | 2014-06-24 | 2014-10-15 | 天津泰达绿化集团有限公司 | Greening method for controlling and discharging salt stepwise in coastal saline area |
CN104751011A (en) * | 2015-04-23 | 2015-07-01 | 中国水利水电科学研究院 | Critical burial depth computing method for groundwater recharging vegetation |
CN106770928A (en) * | 2016-11-11 | 2017-05-31 | 陕西省地质环境监测总站 | The assay method of the lower WATER LEVEL CHANGES threshold of ecology fragility mining area vegetation constraint |
CN109061105A (en) * | 2018-08-02 | 2018-12-21 | 中国水利水电科学研究院 | A kind of calculation method of the critical underground water buried depth of the soil salinization |
CN110059980A (en) * | 2019-04-28 | 2019-07-26 | 安徽农业大学 | A kind of controllable groundwater level depth crop water Sensitivity Index calculation method |
CN110503330A (en) * | 2019-08-20 | 2019-11-26 | 中国科学院新疆生态与地理研究所 | A kind of Arid Region oasis water resource monitoring method and system |
CN110851969A (en) * | 2019-11-01 | 2020-02-28 | 中国人民解放军63653部队 | Hydrogeological feature and regional groundwater circulation simulation method |
CN112085409A (en) * | 2020-09-18 | 2020-12-15 | 中国水利水电科学研究院 | Underground water level regulation and control method for synergistically relieving desertification and secondary salinization of arid region |
CN112115410A (en) * | 2020-09-18 | 2020-12-22 | 中国水利水电科学研究院 | Method for estimating water resource quantity complementary potential of desertification and salinization in arid region |
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20020144953A1 (en) * | 1995-05-05 | 2002-10-10 | Kerfoot William B. | Gas-gas-water treatment system for groundwater and soil remediation |
CN101385433A (en) * | 2008-09-19 | 2009-03-18 | 中国科学院新疆生态与地理研究所 | Groundwater buried depth abductive technique of desert plant seedling root system growth |
JP2012044913A (en) * | 2010-08-26 | 2012-03-08 | Nishizawa:Kk | Leaf vegetable harvester |
CN103093114A (en) * | 2013-02-05 | 2013-05-08 | 河海大学 | Distributed-type river basin water deficit calculating method based on terrain and soil characteristics |
CN104094693A (en) * | 2014-06-24 | 2014-10-15 | 天津泰达绿化集团有限公司 | Greening method for controlling and discharging salt stepwise in coastal saline area |
CN104751011A (en) * | 2015-04-23 | 2015-07-01 | 中国水利水电科学研究院 | Critical burial depth computing method for groundwater recharging vegetation |
CN106770928A (en) * | 2016-11-11 | 2017-05-31 | 陕西省地质环境监测总站 | The assay method of the lower WATER LEVEL CHANGES threshold of ecology fragility mining area vegetation constraint |
CN109061105A (en) * | 2018-08-02 | 2018-12-21 | 中国水利水电科学研究院 | A kind of calculation method of the critical underground water buried depth of the soil salinization |
CN110059980A (en) * | 2019-04-28 | 2019-07-26 | 安徽农业大学 | A kind of controllable groundwater level depth crop water Sensitivity Index calculation method |
CN110503330A (en) * | 2019-08-20 | 2019-11-26 | 中国科学院新疆生态与地理研究所 | A kind of Arid Region oasis water resource monitoring method and system |
CN110851969A (en) * | 2019-11-01 | 2020-02-28 | 中国人民解放军63653部队 | Hydrogeological feature and regional groundwater circulation simulation method |
CN112085409A (en) * | 2020-09-18 | 2020-12-15 | 中国水利水电科学研究院 | Underground water level regulation and control method for synergistically relieving desertification and secondary salinization of arid region |
CN112115410A (en) * | 2020-09-18 | 2020-12-22 | 中国水利水电科学研究院 | Method for estimating water resource quantity complementary potential of desertification and salinization in arid region |
Non-Patent Citations (6)
Title |
---|
张越等: "河套灌区地下水时空变化趋势成因分析", 人民黄河, vol. 38, no. 11, 18 November 2016 (2016-11-18), pages 44 * |
汪勇等: "基于生态圈层结构稳定的地下水位计算与调控", 水科学进展, vol. 32, no. 4, 15 March 2021 (2021-03-15), pages 597 * |
王水献;吴彬;杨鹏年;李霞;董新光;: "焉耆盆地绿洲灌区生态安全下的地下水埋深合理界定", 资源科学, vol. 33, no. 03, 15 March 2011 (2011-03-15), pages 422 * |
玛依拉·麦麦提艾力: "吐鲁番地区土地利用/覆盖变化与生态经济耦合关系研究", 中国优秀硕士学位论文全文数据库工程科技Ⅰ辑, no. 10, 15 October 2013 (2013-10-15), pages 027 - 71 * |
陈敏建等: "西辽河平原地下水补给植被的临界埋深", 水科学进展, vol. 30, no. 1, 26 December 2018 (2018-12-26), pages 24 * |
鲁建荣等: "鄱阳湖典型洲滩湿地植物根系对水分垂向通量的影响", 中国环境科学, vol. 40, no. 5, 20 May 2020 (2020-05-20), pages 2180 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115658822A (en) * | 2022-09-23 | 2023-01-31 | 中国水利水电科学研究院 | Computing method for determining reasonable ecological pattern of farming-pasturing staggered zone in semiarid region |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Ma et al. | Hydrological connectivity from glaciers to rivers in the Qinghai–Tibet Plateau: roles of suprapermafrost and subpermafrost groundwater | |
Raidla et al. | A δ18O isoscape for the shallow groundwater in the Baltic Artesian Basin | |
Kumar | Assessment of groundwater potential | |
Li et al. | Variation of groundwater hydrochemical characteristics in the plain area of the Tarim Basin, Xinjiang Region, China | |
Petrič et al. | Effects of the vadose zone on groundwater flow and solute transport characteristics in mountainous karst aquifers–the case of the Javorniki–Snežnik massif (SW Slovenia) | |
Delsman | Saline groundwater-Surface water interaction in coastal lowlands | |
Liu et al. | Sediment yield modeling using SWAT model: case of Changjiang River Basin | |
Odhiambo et al. | Spatial and isotopic analysis of watershed soil loss and reservoir sediment accumulation rates in Lake Anna, Virginia, USA | |
Xi et al. | Groundwater storage changes and estimation of stream lateral seepage to groundwater in desert riparian forest region | |
CN113656756A (en) | Method for calculating critical burial depth of boundary groundwater between oasis and transition zone of inland river arid region | |
Gao et al. | Vulnerability of eco-hydrological environment in the Yellow River Delta wetland | |
CN113656756B (en) | Method for calculating critical buried depth of groundwater at boundary of oasis and transition zone in arid region of inland river | |
Adar | Quantification of aquifer recharge distribution using environmental isotopes and regional hydrochemistry | |
Islaih et al. | Impact of climate change on flash floods using hydrological modelling and GIS: Case study Zarqa Ma’in area | |
TOKA et al. | Effects of the Vadose Zone on Groundwater Flow and Solute Transport Characteristics in Mountainous Karst Aquifers–the Case of the Javorniki–Snežnik Massif (SW Slovenia) | |
Kim et al. | Suitability analysis of saline water intrusion monitoring wells near a waterway based on a numerical model and electric conductivity time series | |
Duranel et al. | Modelling the hydrological interactions between a fissured granite aquifer and a valley mire in the Massif Central, France | |
Várallyay et al. | NIolel for the IEstimation of Water (and Solute)" I ransport from the Groundwater to Overlying Soil Horizons | |
Kawade et al. | Carbon Dioxide Emission from a Reservoir in India | |
Reghu et al. | Demarcation of groundwater prospective zones in humid tropical river basin: a geospatial approach | |
Wang et al. | Scale effects of eroded sediment transport in Wujiang River Basin, Guizhou Province, China | |
Corniello et al. | Hydrogeological map of the Monti Picentini Regional Park (southern Italy) at 1: 50,000 scale. | |
Mjemah et al. | Hydrogeological mapping and estimation of potential evapotranspiration and recharge rate of Quaternary sand aquifers in Dar-es-Salaam, Tanzania | |
Xiaoyong et al. | The influence of change of land-use type on hydrological cycle | |
Ni et al. | Based on GMS management of shallow groundwater resource in Ningjin, China |
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 |