CN117037932A - Method for calculating farmland supporting ditch water withdrawal amount and identifying water age distribution based on isotopes - Google Patents
Method for calculating farmland supporting ditch water withdrawal amount and identifying water age distribution based on isotopes Download PDFInfo
- Publication number
- CN117037932A CN117037932A CN202310919136.1A CN202310919136A CN117037932A CN 117037932 A CN117037932 A CN 117037932A CN 202310919136 A CN202310919136 A CN 202310919136A CN 117037932 A CN117037932 A CN 117037932A
- Authority
- CN
- China
- Prior art keywords
- water
- farmland
- branch
- ditch
- branch ditch
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 144
- 238000000034 method Methods 0.000 title claims abstract description 33
- 239000003621 irrigation water Substances 0.000 claims abstract description 53
- QVGXLLKOCUKJST-UHFFFAOYSA-N atomic oxygen Chemical compound [O] QVGXLLKOCUKJST-UHFFFAOYSA-N 0.000 claims abstract description 20
- 229910052760 oxygen Inorganic materials 0.000 claims abstract description 20
- 239000001301 oxygen Substances 0.000 claims abstract description 20
- 239000003673 groundwater Substances 0.000 claims abstract description 18
- 230000000694 effects Effects 0.000 claims abstract description 13
- 229910052704 radon Inorganic materials 0.000 claims abstract description 10
- SYUHGPGVQRZVTB-UHFFFAOYSA-N radon atom Chemical compound [Rn] SYUHGPGVQRZVTB-UHFFFAOYSA-N 0.000 claims abstract description 10
- 238000012544 monitoring process Methods 0.000 claims description 9
- 238000004364 calculation method Methods 0.000 claims description 7
- 238000005194 fractionation Methods 0.000 claims description 6
- 238000002156 mixing Methods 0.000 claims description 5
- 238000012360 testing method Methods 0.000 claims description 5
- XOFYZVNMUHMLCC-ZPOLXVRWSA-N prednisone Chemical compound O=C1C=C[C@]2(C)[C@H]3C(=O)C[C@](C)([C@@](CC4)(O)C(=O)CO)[C@@H]4[C@@H]3CCC2=C1 XOFYZVNMUHMLCC-ZPOLXVRWSA-N 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 239000000203 mixture Substances 0.000 claims description 2
- 230000007547 defect Effects 0.000 abstract description 2
- 230000001502 supplementing effect Effects 0.000 abstract description 2
- 230000008595 infiltration Effects 0.000 abstract 7
- 238000001764 infiltration Methods 0.000 abstract 7
- 230000002262 irrigation Effects 0.000 description 6
- 238000003973 irrigation Methods 0.000 description 6
- 230000014759 maintenance of location Effects 0.000 description 4
- 238000002791 soaking Methods 0.000 description 3
- 235000015097 nutrients Nutrition 0.000 description 2
- 239000000575 pesticide Substances 0.000 description 2
- 230000001373 regressive effect Effects 0.000 description 2
- 239000002689 soil Substances 0.000 description 2
- 239000000126 substance Substances 0.000 description 2
- FAPWRFPIFSIZLT-UHFFFAOYSA-M Sodium chloride Chemical compound [Na+].[Cl-] FAPWRFPIFSIZLT-UHFFFAOYSA-M 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000003657 drainage water Substances 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000001788 irregular Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012856 packing Methods 0.000 description 1
- 239000012466 permeate Substances 0.000 description 1
- 238000007789 sealing Methods 0.000 description 1
- 239000011780 sodium chloride Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Mining
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
Abstract
The invention discloses a method for calculating the water withdrawal amount of a branch ditch of a farmland and identifying the water age distribution based on isotopes, which is used for respectively calculating the underwater infiltration rate of irrigation water into the underground of each farmland unit and the lateral infiltration rate of irrigation water into the branch ditch according to the oxygen isotope content and radon isotope activity, and calculating the water quantity of irrigation water lateral infiltration into the branch ditch of each farmland unit, the total groundwater quantity for supplementing the water withdrawn from the branch ditch, the proportion of the water withdrawn from the branch ditch of the lower infiltration into the branch ditch, the proportion of the water withdrawn from the branch ditch of the lateral infiltration into the branch ditch and the proportion of the groundwater entering the branch ditch into the branch ditch of the branch ditch, so that the defects that the prior art cannot accurately track the underwater infiltration rate of irrigation water of the farmland, the proportion of the lateral infiltration discharge into the branch ditch and the groundwater supply effect in the branch ditch withdrawal process are overcome. According to the invention, the age distribution condition of the farmland-branch ditch regression water body is obtained by drawing a gamma curve.
Description
Technical Field
The invention belongs to the field of agricultural hydraulic engineering, relates to a monitoring method for researching a drainage process of a farmland irrigation ditch, and particularly relates to a technology for analyzing the amount of returned water, water source components and water ages of farmland regression water based on multiple isotopes.
Background
Most areas still adopt the irrigation mode of large water flood irrigation, but the drainage system of farmland ditch is not matched and is not sound, and the irrigation is often caused to have or not to have drainage, hold and not leak or drain, irrigate seepage and get into groundwater, destroyed the balance of the regional groundwater volume of irrigating, the groundwater level rises sharply, leads to the saline and alkaline buried in the air-packing area to accumulate to the surface soil rapidly, causes the regional soil secondary salinization of filling. On the other hand, the farmland is short in water-withdrawal duration, the ditch flood peak is formed rapidly, the water quantity is huge, the water-withdrawal contains a large amount of surplus nutrient substances and residual pesticides, and the excessive nutrient substances and the residual pesticides are directly discharged into the ditch step by step to enter the river, so that the water body in the whole area is easily polluted by the non-point source.
The agricultural activities in the areas are not uniform in development time, irregular in agricultural space distribution and large in difference of agricultural water irrigation measures, so that the farmland water-draining process is complex. The ditches in the space area are crisscrossed, the geographical boundary and the occurrence position of the concrete field-water withdrawal are difficult to identify and determine, and the ditch has larger space heterogeneity. The flow production and converging processes are overlapped in time, and meanwhile, the randomness and uncertainty of the water withdrawal time sequence process are increased due to the influence of the adjustment of the rainfall and irrigation and drainage system. The existing technical method is difficult to effectively identify the farmland water-withdrawal process, and the farmland water-withdrawal quantity and quality monitoring has the problems of high monitoring cost, high monitoring difficulty, low monitoring precision and the like.
Disclosure of Invention
The invention aims at solving the problems that the prior art cannot accurately track the proportion of the farmland irrigation water seepage and side seepage discharged into a ditch, the supply effect of underground water in the ditch water withdrawal process and the retention time (water age) distribution of the return water in the ditch, and provides a method for calculating the farmland ditch water withdrawal amount and identifying the water age distribution based on isotopes.
The above object of the present invention is achieved by the following technical solutions:
the method for calculating the farmland supporting ditch water withdrawal amount based on isotopes comprises the following steps:
step 1, analyzing position arrangement and quantity information of farmland and ditch units in a region to be researched, and determining effective flow areas of each farmland unit to each branch ditch;
step 2, collecting water samples of irrigation water, farmland water and branch ditch units, testing oxygen isotope content of the water samples, and calculating the ratio of the irrigation water of each farmland unit to the underground;
step 3, monitoring the radon isotope activities of irrigation water, farmland steeping water and branch ditch units, and calculating the side seepage rate of farmland irrigation water in each farmland unit into the branch ditch by combining the elevation difference between the steeping water surface of each farmland unit and the adjacent water surface of each branch ditch;
and 4, calculating the total amount of underground water in the drainage amount of each branch ditch and the drainage amount of the side seepage of the irrigation water of the farmland unit into the branch ditches by using an end member mixing model, and calculating the influence degree of the underwater seepage of the irrigation water of the farmland unit into the branch ditches, the side seepage into the branch ditches and the regional underground water on the drainage amount of the branch ditches.
Further, the effective flow area from each farmland unit m to each branch ditch n in step 1Wherein M is the number m=1 of the farmland unit, …, M, N is the number of the branch ditch, n=1, …, N, a m Indicating the area of each farmland unit->Represents the adjacent length of each branch n to each farmland unit m>
Further, in the step 2, the clean sample bottle is used for collecting irrigation water in the area, the water for soaking the farmland of each farmland unit m and the drainage samples in each branch ditch n, and the samples are filtered and then stored at a low temperature (not higher than 4 ℃) in a sealing way. The content of oxygen isotopes in the detected sample is respectively recorded as delta I 、δ m 、δ n The ratio f1 of the irrigation water of each farmland unit m to the underground m The calculation formula is as follows:
wherein delta E The average water vapor oxygen isotope content above the regional farmland steeping water is represented by solving for:
ε + ＝α + -1
α + ( 18 O)＝exp[-7.685/10 -3 +6.7123/(273.15+T)-1666.4/(273.15+T) 2 +350410/(273.15+T) 3 ])
ε K ( 18 O)＝0.0277(1-h)
wherein,represents the average oxygen isotope content epsilon of the steeping water and the branch furrows of each farmland unit + Is the equilibrium separation coefficient of isotope, alpha + Is a function of the regional air temperature T, epsilon K Is the dynamic fractionation coefficient of the isotope, h is the relative humidity of the atmosphere at the 2m position of the region, delta A Represents the average oxygen isotope content epsilon in the atmosphere of the region K ( 18 O) represents the kinetic fractionation coefficient of the oxygen isotope.
Further, in the step 3, the radon isotope activity in the water sample is collected by directly testing irrigation water, field water soaking water of each farmland unit m and drainage water in each branch ditch n on site by using a radon isotope detector 222 Rn) are respectively denoted as R I 、R m 、R n The irrigation water side of each farmland unit m seeps into the branch ditch nThe calculation is as follows:
wherein,for the height difference between the water surface of the farmland unit m and the water surface of the branch ditch n which is directly adjacent to the farmland unit m,/->Is the elevation difference between the ground of the farmland unit m and the bottom surface of the directly adjacent branch ditch n.
Further, in step 4, farm Tian Shanyuan m is irrigatedThe water quantity of the water side seeping into the branch groove nThe method comprises the following steps:
wherein IR m Irrigation water quantity for the farmland unit m;
for branch channel n, the water amount of the branch channel is collected into the upper bucket channelThe method comprises the following steps:
wherein,for the side seepage of farmland irrigation water of a farmland unit m into the water quantity of a branch ditch n, GW n The total amount of groundwater to be supplied to the branch ditch n;
based on the isotope mass balance end member mixture model, GW can be obtained n 、
For the branch ditch n, each farmland unit m irrigates and seeps into the branch ditch under water to occupy the water yield of the branch ditch nProportion of->Calculated by the following formula:
irrigation water side seepage entering branch ditch of farmland unit m occupies the amount of water returned by branch ditch nProportion of->Calculated by the following formula:
the groundwater in the region enters the branch ditch to occupy the amount of the water returned by the branch ditch nProportion of->Calculated by the following formula:
the residence time distribution of farmland-branch ditch regression water return is analyzed based on the method, namely the water age distribution of farmland-branch ditch regression water return is analyzed, and the specific operation is as follows:
constructing a distribution gamma curve h of water detention time, namely water age tau, in the branch ditch n n (τ)：
Wherein,the distance from the farmland unit m to the branch ditch n to the upper bucket ditch is the distance from the farmland unit m to the branch ditch n; />The average value of the water retention time in the branch channel n is shown.
The invention has the beneficial effects that:
according to the method, the ratio of the irrigation water side seepage of each farmland unit into the underground is calculated according to the oxygen isotope content, the ratio of the irrigation water side seepage of each farmland unit into the branch ditch is calculated according to the radon isotope activity, the water quantity of the irrigation water side seepage of each farmland unit into the branch ditch is calculated by combining the irrigation water quantity of the farmland unit, the total groundwater quantity of the water to be supplied to the branch ditch for water withdrawal, the ratio of the irrigation water side seepage of each farmland unit into the branch ditch for the water withdrawal of the branch ditch and the ratio of the groundwater entering the area into the branch ditch for the water withdrawal of the branch ditch are calculated by an end member mixing model with the isotope mass balance, so that the influence degree of the irrigation water side seepage of the farmland unit into the branch ditch, the side seepage of the branch ditch and the area groundwater on the water withdrawal of the branch ditch can not be accurately tracked, and the defects of the prior art in the groundwater supply effect in the branch ditch water withdrawal process can be solved.
According to the invention, the residence time distribution of farmland-branch ditch regression water is obtained by drawing a gamma curve, namely the water age distribution condition of farmland-branch ditch regression water.
Drawings
FIG. 1 is a technical flow chart of the method of the present invention.
FIG. 2 is a schematic diagram showing the arrangement of the farmland and trench units in the embodiment.
FIG. 3 shows the delta of irrigation water, farmland water and a gutter unit according to the embodiment 18 Box plot of O.
FIG. 4 shows irrigation water, farmland water and a gutter unit according to an embodiment 222 Box plot of Rn activity.
FIG. 5 is a graph showing the average influence ratio of the groundwater in the underground side-seepage branch ditches and the area to the water withdrawal amount of each branch ditches in the embodiment.
FIG. 6 is a graph showing a gamma curve of the retention time (age τ) of water in the rural ditches in the example.
Detailed Description
The following description of the embodiments of the present invention will be made more apparent and fully hereinafter with reference to the accompanying drawings, in which some, but not all embodiments of the invention are shown. All other embodiments, which can be made by one of ordinary skill in the art without undue burden on the person of ordinary skill in the art based on embodiments of the present invention, are within the scope of the present invention.
Example 1
Taking Ningxia Pingroo advancing farm as an example, corresponding investigation and monitoring are carried out on farmland supporting ditches water return condition in 2023 and 5 months.
As shown in fig. 1, the identification method of the farmland branch ditch water withdrawal amount, the water source component and the water body age based on isotopes comprises the following steps:
step 1, analyzing the position arrangement and the quantity information of farmland and ditch units in the area, and determining the effective flow area of each farmland unit to each branch ditch, wherein the method comprises the following specific operations:
1) Determining the number of regional branches n=14 and the length L of each branch n N represents the number of the branch groove (n=1, …, 14).
2) Dividing farmland units into regions, determining the number as M=13, and determining the area A of each farmland unit m M is denoted as the number of farmland units (m=1, …, 13).
3) Determining adjacent lengths of each branch ditch n and each farmland unit m according to the positions and arrangement information of the farmland and each ditch unitThereby calculating the effective flow area of each farmland unit m to each branch ditch n>As shown in FIG. 2, the green and yellow areas represent farmland, and the blue areas represent ditches.
And 2, collecting water samples of irrigation water, farmland water and branch ditch units, testing the oxygen isotope content of the water samples, and calculating the ratio of the irrigation water of each farmland unit to the underground.
The specific operation is as follows:
the irrigation water in the collecting area of the clean sample bottle, the soaking water in each farmland unit m and the drainage sample in each branch ditch n are filtered and stored at a low temperature (not higher than 4 ℃). Determination of the oxygen isotope content (delta) in sample Water Using a liquid Water isotope Analyzer (PICARROL 2130-i, picaro, USA) 18 O). Respectively marked as delta I 、δ m 、δ n As shown in fig. 3.
Proportion f1 of the irrigation water of each farmland unit m to the underground m The calculation formula is as follows:
wherein delta E The average water vapor oxygen isotope content above the regional farmland steeping water is represented by solving for:
ε + ＝α + -1
α + ( 18 O)＝exp[-7.685/10 -3 +6.7123/(273.15+T)-1666.4/(273.15+T) 2 +350410/(273.15+T) 3 ])
ε K ( 18 O)＝0.0277(1-h)
wherein,represents the average oxygen isotope content epsilon of the steeping water and the branch furrows of each farmland unit + Is the equilibrium separation coefficient of isotope, alpha + Is a function of the regional air temperature T (DEG C), delta A Represents the average oxygen isotope content epsilon in the atmosphere of the region K ( 18 O) represents the kinetic fractionation coefficient of the oxygen isotope, ε K Is the kinetic fractionation coefficient of the isotope, h is the relative humidity of the atmosphere at region 2m, and is measured to be t=22.5 ℃, h=32.2%.
F1 in the present embodiment m The average value is 45.22%, the minimum value is 22.68%, and the maximum value is 65.72%.
Step 3, monitoring radon isotope activities of irrigation water, farmland steeping water and branch ditch units, and calculating the side seepage rate of farmland irrigation water in each farmland unit into the branch ditch by combining the elevation difference of the steeping water surface of each farmland unit and the directly adjacent water surface of each branch ditch, wherein the specific operation is as follows:
directly testing irrigation water, field water of each farmland unit m and radon isotope activity in drainage collection water sample in each branch ditch n on site by using radon isotope detectors (RAD 7 and RAD AQUA water fittings) 222 Rn) are respectively denoted as R I 、R m 、R n As shown in fig. 4.
Proportion of each farmland unit m to side-seep irrigation water into the branch ditch nThe calculation is as follows:
wherein,for the height difference between the water surface of the farmland unit m and the water surface of the branch ditch n which is directly adjacent to the farmland unit m,/->Is the elevation difference between the ground of the farmland unit m and the bottom surface of the directly adjacent branch ditch n. Calculated (i.e. the case of->The average value is 27.73%, the minimum value is 11.28%, and the maximum value is 51.16%.
And 4, analyzing different water source components in the drainage process of each branch ditch by using an end member mixing model, and respectively determining the influence degree of the farmland unit irrigation water seepage into the underground, the side seepage into the branch ditches and the regional groundwater on the drainage quantity of the branch ditches. As shown in fig. 5, the contribution of 13 farmland unit irrigation water to the water withdrawal average water quantity of 14 branch ditches is shown as follows:
for branch channel n, the water amount of the branch channel is collected into the upper bucket channelThe method comprises the following steps:
wherein,for the side seepage of farmland irrigation water of a farmland unit m into the water quantity of a branch ditch n, GW n The total amount of groundwater to be supplied to the branch trench n.
For the followingThen there is
Wherein IR m The irrigation water quantity of the farmland unit m.
The GW can be obtained by considering the end member mixing model of isotope mass balance n 、
Thus, for the branch trench n, each farmland unit n irrigates water to permeate into the branch trench n for water withdrawalProportion of->It is possible to calculate:
irrigation water side seepage entering branch ditch of farmland unit m occupies the amount of water returned by branch ditch nProportion of->It is possible to calculate:
groundwater pair supplementing branch ditch n water-withdrawal quantityProportion of->It is possible to calculate:
and 5, analyzing residence time (water age) distribution of farmland-branch ditch regression water withdrawal, wherein the method comprises the following specific operations:
constructing a gamma curve h of water retention time (water age tau) distribution in the branch ditch n n (τ)：
Wherein,the distance from the farmland unit m to the branch ditch n to the upper bucket ditch is the distance from the farmland unit m to the branch ditch n; />The average value of the water retention time in the branch channel n is shown.
h n (τ) can reflect the distribution of the retention time of farmland-by-gutter regressive water (water age τ), as shown in FIG. 6, showing the cumulative integral of the retention time of farmland-by-gutter regressive water (water age τ) of the by-gutters 1-7Cloth, branch grooves 1-7And (3) days.
Claims (6)
1. The method for calculating the farmland supporting ditch water withdrawal amount based on the isotopes is characterized by comprising the following steps of:
step 1, analyzing position arrangement and quantity information of farmland and ditch units in a region to be researched, and determining effective flow areas of each farmland unit to each branch ditch;
step 2, collecting water samples of irrigation water, farmland water and branch ditch units, testing oxygen isotope content of the water samples, and calculating the ratio of the irrigation water of each farmland unit to the underground;
step 3, monitoring the radon isotope activities of irrigation water, farmland steeping water and branch ditch units, and calculating the side seepage rate of farmland irrigation water in each farmland unit into the branch ditch by combining the elevation difference between the steeping water surface of each farmland unit and the adjacent water surface of each branch ditch;
and 4, calculating the total amount of underground water in the drainage amount of each branch ditch and the drainage amount of the side seepage of the irrigation water of the farmland unit into the branch ditches by using an end member mixing model, and calculating the influence degree of the underwater seepage of the irrigation water of the farmland unit into the branch ditches, the side seepage into the branch ditches and the regional underground water on the drainage amount of the branch ditches.
2. The method for calculating the water withdrawal amount of a branch ditch of a farmland based on isotopes according to claim 1, characterized in that the effective flow area from each farmland unit m to each branch ditch n in step 1Wherein M is the number m=1 of the farmland unit, …, M, N is the number of the branch ditch, n=1, …, N, a m Indicating the area of each farmland unit->Represents the adjacent length of each branch n to each farmland unit m>
3. The method for calculating the water withdrawal amount of a branch ditch of a farmland based on isotopes according to claim 2, wherein the oxygen isotope content in the water samples of the irrigation water, the farmland steeping water and the branch ditch unit in the step 2 is respectively recorded as delta I 、δ m 、δ n The ratio f1 of the irrigation water of each farmland unit m to the underground m The calculation formula is as follows:
wherein delta E The average water vapor oxygen isotope content above the regional farmland steeping water is represented by solving for:
ε + ＝α + -1
α + ( 18 O)＝exp[-7.685/10 -3 +6.7123/(273.15+T)-1666.4/(273.15+T) 2 +350410/(273.15+T) 3 ])
ε K ( 18 O)＝0.0277(1-h)
wherein,represents the average oxygen isotope content epsilon of the steeping water and the branch furrows of each farmland unit + Is the equilibrium separation coefficient of isotope, alpha + Is a function of the regional air temperature T, epsilon K Is the dynamic fractionation of isotopesCoefficient h is the relative humidity of the atmosphere at region 2m, delta A Represents the average oxygen isotope content epsilon in the atmosphere of the region K ( 18 O) represents the kinetic fractionation coefficient of the oxygen isotope.
4. The method for calculating the water withdrawal amount of a branch ditch of a farmland based on isotopes according to claim 3, wherein the activities of radon isotopes in the water samples of the irrigation water, the farmland steeping water and the branch ditch unit in the step 3 are respectively marked as R I 、R m 、R n The irrigation water side of each farmland unit m seeps into the branch ditch nThe calculation is as follows:
wherein,for the height difference between the water surface of the farmland unit m and the water surface of the branch ditch n which is directly adjacent to the farmland unit m,/->Is the elevation difference between the ground of the farmland unit m and the bottom surface of the directly adjacent branch ditch n.
5. The method for isotope-based calculation of the water withdrawal amount from a branch trench in a farmland according to claim 4, wherein in step 4, the irrigation water from the farmland unit m side infiltrates into the branch trench n in an amount of waterThe method comprises the following steps:
wherein IR m Irrigation water quantity for the farmland unit m;
for branch channel n, the water amount of the branch channel is collected into the upper bucket channelThe method comprises the following steps:
wherein,for the side seepage of farmland irrigation water of a farmland unit m into the water quantity of a branch ditch n, GW n The total amount of groundwater to be supplied to the branch ditch n;
based on the isotope mass balance end member mixture model, GW can be obtained n 、
For the branch ditch n, each farmland unit m irrigates and seeps into the branch ditch under water to occupy the water yield of the branch ditch nProportion of->Calculated by the following formula:
irrigation water side seepage entering branch ditch of farmland unit m occupies the amount of water returned by branch ditch nProportion of->Calculated by the following formula:
the groundwater in the region enters the branch ditch to occupy the amount of the water returned by the branch ditch nProportion of->Calculated by the following formula:
6. the method for calculating the water withdrawal amount of the farmland branch ditch based on isotopes according to claim 5, further comprising the step 5 of analyzing the residence time distribution of the farmland-branch ditch regression water, namely the water age distribution of the farmland-branch ditch regression water, specifically comprising the following steps:
constructing a distribution gamma curve h of water detention time, namely water age tau, in the branch ditch n n (τ)：
Wherein,the distance from the farmland unit m to the branch ditch n to the upper bucket ditch is the distance from the farmland unit m to the branch ditch n; />The average value of the water retention time in the branch channel n is shown; alpha is a distance parameter and beta is a water quantity variation parameter.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310919136.1A CN117037932B (en) | 2023-07-25 | 2023-07-25 | Method for calculating farmland supporting ditch water withdrawal amount and identifying water age distribution based on isotopes |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310919136.1A CN117037932B (en) | 2023-07-25 | 2023-07-25 | Method for calculating farmland supporting ditch water withdrawal amount and identifying water age distribution based on isotopes |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN117037932A true CN117037932A (en) | 2023-11-10 |
| CN117037932B CN117037932B (en) | 2024-03-19 |
Family
ID=88621787
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310919136.1A Active CN117037932B (en) | 2023-07-25 | 2023-07-25 | Method for calculating farmland supporting ditch water withdrawal amount and identifying water age distribution based on isotopes |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN117037932B (en) |
Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080052011A1 (en) * | 2006-07-29 | 2008-02-28 | Yongdong Wang | Determination of chemical composition and isotope distribution with mass spectrometry |
| CN102160495A (en) * | 2010-12-23 | 2011-08-24 | 西北农林科技大学 | Method for eliminating and degrading non-point source pollutants in field in situ |
| US20130118994A1 (en) * | 2011-03-21 | 2013-05-16 | R360 Environmental Solutions | System and Methods for Wastewater and Produced Water Cleaning and Reclamation |
| CN107764960A (en) * | 2017-09-15 | 2018-03-06 | 国家电网公司 | Pumped Storage Plant Underground Powerhouse infiltration source method is surveyed based on isotope tracer technique |
| CN109470668A (en) * | 2018-11-14 | 2019-03-15 | 华东理工大学 | A kind of combination nitrogen oxygen isotope and three-dimensional fluorescence carry out the method that pollutant is traced to the source |
| CN110674985A (en) * | 2019-09-20 | 2020-01-10 | 北京建筑大学 | Urban resident domestic water consumption prediction method and application thereof |
| WO2020155510A1 (en) * | 2019-01-28 | 2020-08-06 | 湖南恒凯环保科技投资有限公司 | Soil conditioner used for restoring mercury-contaminated agricultural land, and preparation method and application thereof |
| CN114496108A (en) * | 2022-01-21 | 2022-05-13 | 中国地质调查局武汉地质调查中心(中南地质科技创新中心) | Water phosphate tracing method based on phosphate oxygen isotope and water chemistry characteristics |
| CN115754248A (en) * | 2022-12-09 | 2023-03-07 | 中国科学院南京土壤研究所 | Diving-soil moisture movement tracing experimental device and method under freezing condition based on radon isotope |
| CN116254321A (en) * | 2022-12-16 | 2023-06-13 | 中国农业科学院农业资源与农业区划研究所 | Accurate quantification method for individual nitrogen utilization of soil microorganisms |
| CN116298184A (en) * | 2022-12-07 | 2023-06-23 | 中国科学院南京土壤研究所 | Tracing method for farmland soil evapotranspiration water vapor migration process by utilizing oxyhydrogen isotope |
-
2023
- 2023-07-25 CN CN202310919136.1A patent/CN117037932B/en active Active
Patent Citations (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20080052011A1 (en) * | 2006-07-29 | 2008-02-28 | Yongdong Wang | Determination of chemical composition and isotope distribution with mass spectrometry |
| CN102160495A (en) * | 2010-12-23 | 2011-08-24 | 西北农林科技大学 | Method for eliminating and degrading non-point source pollutants in field in situ |
| US20130118994A1 (en) * | 2011-03-21 | 2013-05-16 | R360 Environmental Solutions | System and Methods for Wastewater and Produced Water Cleaning and Reclamation |
| CN107764960A (en) * | 2017-09-15 | 2018-03-06 | 国家电网公司 | Pumped Storage Plant Underground Powerhouse infiltration source method is surveyed based on isotope tracer technique |
| CN109470668A (en) * | 2018-11-14 | 2019-03-15 | 华东理工大学 | A kind of combination nitrogen oxygen isotope and three-dimensional fluorescence carry out the method that pollutant is traced to the source |
| WO2020155510A1 (en) * | 2019-01-28 | 2020-08-06 | 湖南恒凯环保科技投资有限公司 | Soil conditioner used for restoring mercury-contaminated agricultural land, and preparation method and application thereof |
| CN110674985A (en) * | 2019-09-20 | 2020-01-10 | 北京建筑大学 | Urban resident domestic water consumption prediction method and application thereof |
| CN114496108A (en) * | 2022-01-21 | 2022-05-13 | 中国地质调查局武汉地质调查中心(中南地质科技创新中心) | Water phosphate tracing method based on phosphate oxygen isotope and water chemistry characteristics |
| CN116298184A (en) * | 2022-12-07 | 2023-06-23 | 中国科学院南京土壤研究所 | Tracing method for farmland soil evapotranspiration water vapor migration process by utilizing oxyhydrogen isotope |
| CN115754248A (en) * | 2022-12-09 | 2023-03-07 | 中国科学院南京土壤研究所 | Diving-soil moisture movement tracing experimental device and method under freezing condition based on radon isotope |
| CN116254321A (en) * | 2022-12-16 | 2023-06-13 | 中国农业科学院农业资源与农业区划研究所 | Accurate quantification method for individual nitrogen utilization of soil microorganisms |
Non-Patent Citations (2)
| Title |
|---|
| 杨丽慧等: "基于排水沟的宁夏银北灌区农田退水污染演变研究", 《中国生态农业学报》, vol. 23, no. 12, 31 December 2015 (2015-12-31), pages 1580 - 1587 * |
| 章斌等: "用氯和氢氧同位素揭示洋戴河平原地下水的形成演化规律", 环境科学学报, no. 11, 6 November 2013 (2013-11-06), pages 50 - 57 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN117037932B (en) | 2024-03-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Mahannopkul et al. | Influences of root concentration and suction on Chrysopogon zizanioides reinforcement of soil | |
| CN102169617B (en) | Early warning method of mud-rock flow used in regions short of rainfall data | |
| Fu et al. | Importance of subsurface fluxes of water, nitrogen and phosphorus from rice paddy fields relative to surface runoff | |
| Klammler et al. | Determining water and nitrogen balances for beneficial management practices using lysimeters at Wagna test site (Austria) | |
| CN105334311A (en) | Multi-layer soil water fertilizer migration and transformation parameters testing device and method | |
| CN104239729B (en) | A kind of agricultural non -point pollution measuring method and device based on water circulation | |
| CN106770928A (en) | The assay method of the lower WATER LEVEL CHANGES threshold of ecology fragility mining area vegetation constraint | |
| Rowe et al. | Disposition of rainfall in two mountain areas of California | |
| Maestre-Valero et al. | Effects of drip irrigation systems on the recovery of dissolved oxygen from hypoxic water | |
| Khanzada et al. | Groundwater uptake and sustainability of Acacia and Prosopis plantations in Southern Pakistan | |
| CN111473819B (en) | Surface runoff monitoring and analyzing method | |
| Klocke et al. | Design, installation, and performance of percolation lysimeters for water quality sampling | |
| CN117037932B (en) | Method for calculating farmland supporting ditch water withdrawal amount and identifying water age distribution based on isotopes | |
| Woodbury et al. | Effectiveness of a passive feedlot runoff control system using a vegetative treatment area for nitrogen control | |
| Broughton | The performance of subsurface drainage systems on two Saint Lawrence lowland soils. | |
| Faulkner et al. | Tracer movement through paired vegetative treatment areas receiving silage bunker runoff | |
| CN111983190A (en) | Soil erosion monitoring system for depression and use method thereof | |
| Wu et al. | Reviews on agricultural nonpoint source pollution monitoring techniques and methods | |
| Feyen et al. | Runoff processes in catchments with a small scale topography | |
| Wierenga | Influence of trickle and surface irrigation on return flow quality | |
| CN214471140U (en) | Device for measuring underground water level of pit and automatically following underground water level change of field | |
| Tong et al. | Long-term remote monitoring the effects of rainfall and agricultural non-point sources pollution on the surface waters quality in lake Taihu | |
| Trudgill et al. | Application of simple soil-water flow models to the transfer of surface-applied solutes to streams | |
| CN115754240A (en) | Method for monitoring non-point source pollution of farmland in hilly and mountainous regions | |
| Shahrajabian et al. | A Review of Lysimeter Studies and Experiments by Considering Agricultural Production |
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 |