CN112268595A - Method for monitoring regional average groundwater level - Google Patents
Method for monitoring regional average groundwater level Download PDFInfo
- Publication number
- CN112268595A CN112268595A CN202011186789.6A CN202011186789A CN112268595A CN 112268595 A CN112268595 A CN 112268595A CN 202011186789 A CN202011186789 A CN 202011186789A CN 112268595 A CN112268595 A CN 112268595A
- Authority
- CN
- China
- Prior art keywords
- level
- groundwater
- area
- groundwater level
- calculating
- 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.)
- Pending
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01F—MEASURING VOLUME, VOLUME FLOW, MASS FLOW OR LIQUID LEVEL; METERING BY VOLUME
- G01F23/00—Indicating or measuring liquid level or level of fluent solid material, e.g. indicating in terms of volume or indicating by means of an alarm
-
- 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
-
- 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
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- Pure & Applied Mathematics (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Computational Mathematics (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Operations Research (AREA)
- Fluid Mechanics (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
The invention relates to the field of groundwater resource monitoring, and discloses a method for monitoring regional average groundwater level, which comprises dividing a research region of water level to be monitored into regional calculation units of different levels; collecting information and basic information of a monitoring well in a water level research area to be monitored; calculating the groundwater level and groundwater level variation of the final stage area calculating unit to obtain the groundwater level and groundwater level variation of all the final stage area calculating units; sequentially and iteratively calculating the groundwater level and the groundwater level variation of each stage of regional calculation unit according to the groundwater level and the groundwater level variation of all final stage regional calculation units; and obtaining the groundwater level and the groundwater level variation of the water level research area to be monitored. According to the invention, the regional average groundwater level and the groundwater level change value can be rapidly and accurately calculated according to the monitoring data of the regional groundwater level monitoring well, and a basis and a support are provided for management work such as groundwater treatment effect evaluation and assessment.
Description
Technical Field
The invention relates to the field of underground water resource monitoring, in particular to a method for monitoring the average underground water level in a region.
Background
Groundwater has important resource supply and ecological maintenance functions, and has irreplaceable effects. The development and utilization scale of underground water is continuously enlarged, and a large amount of underground water is excessively extracted in some regions, so that the aquifer is nearly drained, and a series of serious environmental and geological problems such as ground settlement, ground collapse and ground cracks, seawater (salt water) invasion, land desertification and the like are caused.
The underground water development and utilization mainly adopts a shaft sinking and pumping mode, and takes the occurrence characteristics of underground water into consideration, so that the motor-pumped wells have the characteristics of wide distribution, scattered positions, numerous quantities and the like, and the difficulty in monitoring, metering, monitoring and managing the underground water exploitation amount is higher. In actual work, only a few underground water mining methods are monitored by instruments, but most of the underground water, especially agricultural water, is almost not provided with metering equipment, the underground water mining amount is mainly based on a statistical estimation method, and the precision is difficult to meet the requirements of supervision and examination. When the groundwater level supervision and examination is implemented, the regional average water level and the regional average water level change need to be calculated according to the groundwater level of the monitoring well. The current main calculation methods comprise a kriging interpolation method, an isoline weighting method and the like. The contour weighting method is to draw the planar distribution condition of the regional groundwater by a group of groundwater level contours so as to further obtain the regional average groundwater level. The method can visually reflect the distribution condition of the groundwater level in the space, and is particularly suitable for areas with obvious terrain changes. The method has the disadvantages that more water level data of the monitoring well are needed for drawing the equal water level line, a time-interval contour map is not needed to be drawn again, and the workload is large. The kriging interpolation method is to obtain water level values at undetermined points of other unmonitored water levels by monitoring well water level values through the kriging interpolation method. The method has complex calculation steps and large interpolation calculation amount, is generally used in space calculation software such as GIS and the like, has high requirement on the technical level of managers, and is difficult to popularize in practical application.
Therefore, in view of the urgent need of groundwater level monitoring, a simple and convenient method for calculating the average groundwater level in the region is urgently needed for effectively, accurately and quickly calculating the average groundwater level and the water level amplitude.
Disclosure of Invention
The present invention provides a method for monitoring the average groundwater level in a region, thereby solving the above problems of the prior art.
A method for monitoring the average groundwater level in an area comprises the following steps:
s1) dividing the water level research area to be monitored into N-level area calculation units with different levels, wherein the Nth-level area calculation unit is a final-level area calculation unit;
s2) collecting information of a monitoring well in the research area of the water level to be monitored and basic information of the research area of the water level to be monitored;
s3) calculating the groundwater level and groundwater level variation of the final stage area calculating units according to the information of the monitoring wells in the water level research area to be monitored and the basic information of the water level research area to be monitored, and obtaining the groundwater level and groundwater level variation of all the final stage area calculating units;
s4) sequentially calculating the groundwater level and the groundwater level variation of the N-1 level region calculating unit, the groundwater level and the groundwater level variation of the N-2 level region calculating unit, …, the groundwater level and the groundwater level variation of the first level region calculating unit in an iterative manner according to the groundwater level and the groundwater level variation of all the last level region calculating units;
s5) obtaining the groundwater level and the groundwater level variation of the research area of the water level to be monitored.
Further, in step S1), the N-level different-level region calculating units include a primary region calculating unit, the primary region calculating unit includes a plurality of secondary region calculating units, …, an N-2 level region calculating unit includes a plurality of N-1 level region calculating units, and an N-1 level region calculating unit includes a plurality of N-level region calculating units.
The invention divides the research area of the water level to be monitored into calculation units of different levels, such as a primary unit, a secondary unit, a tertiary unit and the like according to a bottom-up method. Taking the most common administrative districts as an example for further explanation, the groundwater level change value of the provincial administrative district is calculated, the provincial administrative district can be divided into a provincial district calculation unit (first level), a city district calculation unit (second level) and a county district calculation unit (third level) step by step according to the administrative district level, and if necessary, the provincial district can be further divided into smaller calculation units such as towns, villages and villages.
Further, in step S2), the information of the monitoring wells in the research area of the water level to be monitored includes the number of the monitoring wells in the research area, the positions of the monitoring wells, the levels of the monitoring wells and the groundwater levels measured by the monitoring wells at different time nodes; the basic information of the water level research area to be monitored comprises the area of the research area, the area of each level of area calculation unit, the number of monitoring wells in each level of area calculation unit and the distribution positions of the monitoring wells.
The data collected according to the purposes of the invention mainly comprise two types: firstly, monitoring well information of a water level research area to be monitored is collected, wherein the monitoring well information comprises the number, the position and the position of the monitoring well in a research area, the underground water level (buried depth) of a corresponding time node and the like. If the average water level (burial depth) change of the water level research area to be monitored is calculated, for the same monitoring well in the research area, the groundwater monitoring water level (burial depth) data at the previous moment is subtracted from the groundwater monitoring water level (burial depth) data at the later moment, and the groundwater level (burial depth) change value of each monitoring well in the research area is obtained. And secondly, collecting basic information of the water level research area to be monitored, wherein the basic information comprises the area of the water level research area to be monitored, the area of each level of area calculation unit, and the quantity and distribution position conditions of the monitoring wells in each level of area calculation unit.
Further, in step S3), the method for calculating the groundwater level and the groundwater level variation of each final stage area calculating unit according to the information of the monitoring well in the research area of the water level to be monitored and the basic information of the research area of the water level to be monitored to obtain the groundwater level and the groundwater level variation of all the final stage area calculating units includes the following steps:
s31), the total number of the final-stage area calculation units is v, and the measured underground water level and the underground water level variation of all monitoring wells in each final-stage area calculation unit are obtained;
s32) calculating the groundwater level of the jth final stage area calculating unit according to the groundwater levels measured by all monitoring wells in the jth final stage area calculating unit
And groundwater level change amount of jth final stage area calculation unit
Groundwater level of the jth last stage region calculation unit
A groundwater level variation amount of the jth last stage region calculation unit
Underground water level of the tth monitoring well in the jth N-level area calculation unit,
calculating the groundwater level variation of the tth monitoring well in the jth N-stage area calculating unit, wherein t is more than or equal to 1 and less than or equal to N, N is the total number of the monitoring wells in the jth last-stage area calculating unit, and j is more than or equal to 1 and less than or equal to v;
s33) obtaining the groundwater level and the groundwater level variation amount of all the final stage area calculation units.
Further, in step S4), sequentially and iteratively calculating the groundwater level and groundwater level variation of the N-1 stage region calculation unit, the groundwater level and groundwater level variation of the N-2 stage region calculation unit, …, and the groundwater level and groundwater level variation of the one stage region calculation unit based on the groundwater level and groundwater level variation of all the last stage region calculation units, the method includes the following steps:
s41) recording the total number of h-1 level region calculating units as w, recording the total number of h-level region calculating units in the e-1 level region calculating unit as m, wherein the k-th h-level region calculating unit in the m h-level region calculating units comprises r h +1 level region calculating units; e is more than or equal to 1 and less than or equal to w; k is more than or equal to 1 and less than or equal to m; h is more than or equal to 1 and less than or equal to N-1;
s42) acquiring the area of each h + 1-level area calculation unit according to the information of the monitoring well in the water level research area to be monitored and the basic information of the water level research area to be monitored;
s43) calculating the groundwater level of the kth level area calculating unit
Calculating the groundwater level variation of the kth h level area calculation unit
Calculating the area of a unit for the ith h + 1-level region;
s44) sequentially calculating the groundwater level and the groundwater level variation of all the h-level region calculating units to obtain a groundwater level set { D ] of the h-level region calculating unitsh,1、Dh,2、…、Dh,k、…、Dh,mSet of groundwater level variation { Δ D } and level-h regional computing unitsh,1、ΔDh,2、…、ΔDh,k、…、ΔDh,m}；
S45) calculating groundwater level of the e-th h-1 level area calculating unit according to groundwater levels and groundwater level variation of all h level area calculating units
Groundwater level variation of the e-th h-1 level area calculation unit
The area of the cell is calculated for the kth h level region.
The method comprises the steps of firstly, starting from a divided final-stage area computing unit, and obtaining the groundwater level and groundwater level variation of all final-stage area computing units by an arithmetic mean method according to monitoring data of monitoring wells in all final-stage computing units. Then, according to the groundwater level (buried depth) of each final stage region calculation unit, the variation value of the groundwater level (buried depth) and the final stage region calculation unit area which are obtained by calculation, and the subordination relation with the previous-stage area calculation unit is obtained by calculating the underground water level (buried depth) and the underground water level (buried depth) change value of each previous-stage calculation unit by adopting an area weighted average method, continuously iterating, calculating the underground water level (buried depth) and the underground water level (buried depth) change value of each previous-stage unit by adopting the area weighted average method and obtaining the underground water level (buried depth) and the underground water level (buried depth) change value of each previous-stage unit by using the underground water level (buried depth) and the underground water level (buried depth) change value of the current-stage calculation unit and the subordination relation with the previous-stage unit, and finally calculating the underground water level (buried depth) and the underground water level (buried depth) change value of the water level research area to be monitored.
The invention has the beneficial effects that: the method can simply and conveniently calculate the groundwater level of the region, is simple and easy to operate, requires few calculation tools and data, and greatly reduces the calculation cost compared with the methods of a kriging interpolation method, an isoline weighting method and the like in the prior art. The method has the advantages of high calculation speed, high calculation precision and more comprehensive and scientific calculation result, not only considers the current situation of underground water monitoring, but also considers regional area distribution and weight, has better applicability, and is suitable for administrative divisions such as provinces, cities, counties, villages and the like, and is also suitable for water resource drainage basin divisions such as water resource first-level districts, second-level districts, third-level districts and the like.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings used in the embodiments are briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on the drawings without creative efforts.
Fig. 1 is a schematic structural diagram of a method for monitoring the average groundwater level in a region according to a first embodiment of the present invention.
Fig. 2 is a distribution diagram of a water level research area to be monitored and an underground water monitoring well provided in the first embodiment.
Fig. 3 is a partition diagram of a secondary computing unit in the research area of the water level to be monitored according to the first embodiment.
Fig. 4 is a division diagram of a three-level region calculation unit in the second-level unit d according to the first embodiment.
Fig. 5 is a schematic diagram of the distribution of groundwater monitoring wells of the county level unit r in the first embodiment.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to the accompanying drawings. It should be understood that the detailed description and specific examples, while indicating the invention, are intended for purposes of illustration only and are not intended to limit the scope of the invention. It is noted that the terms "comprises" and "comprising," and any variations thereof, in the description and claims of the present invention and the above-described drawings are intended to cover a non-exclusive inclusion, such that a process, method, system, article, or apparatus that comprises a list of elements is not necessarily limited to those elements, but may include other elements not expressly listed or inherent to such process, method, article, or apparatus.
Embodiment one, a method for monitoring regional average groundwater level, as shown in fig. 1, includes the following steps:
s1) dividing the water level research area to be monitored into N-level area calculation units with different levels, wherein the Nth-level area calculation unit is a final-level area calculation unit;
s2) collecting information of a monitoring well in the research area of the water level to be monitored and basic information of the research area of the water level to be monitored;
s3) calculating the groundwater level and groundwater level variation of the final stage area calculating units according to the information of the monitoring wells in the water level research area to be monitored and the basic information of the water level research area to be monitored, and obtaining the groundwater level and groundwater level variation of all the final stage area calculating units;
s4) sequentially calculating the groundwater level and the groundwater level variation of the N-1 level region calculating unit, the groundwater level and the groundwater level variation of the N-2 level region calculating unit, …, the groundwater level and the groundwater level variation of the first level region calculating unit in an iterative manner according to the groundwater level and the groundwater level variation of all the last level region calculating units;
s5) obtaining the groundwater level and the groundwater level variation of the research area of the water level to be monitored.
In step S1), the N-level different-level region calculating units include a primary region calculating unit, the primary region calculating unit includes a plurality of secondary region calculating units, …, an N-2 level region calculating unit includes a plurality of N-1 level region calculating units, and an N-1 level region calculating unit includes a plurality of N-level region calculating units.
The invention divides the research area of the water level to be monitored into calculation units of different levels, such as a primary unit, a secondary unit, a tertiary unit and the like according to a bottom-up method. Taking the most common administrative districts as an example for further explanation, the groundwater level change value of the provincial administrative district is calculated, the provincial administrative district can be divided into a provincial district calculation unit (first level), a city district calculation unit (second level) and a county district calculation unit (third level) step by step according to the administrative district level, and if necessary, the provincial district can be further divided into smaller calculation units such as towns, villages and villages.
In the step S2), the information of the monitoring wells in the water level research area to be monitored comprises the number of the monitoring wells in the research area, the positions of the monitoring wells and the groundwater levels measured by the monitoring wells at different time nodes; the basic information of the water level research area to be monitored comprises the area of the research area, the area of each level of area calculation units, the number of monitoring wells in each level of area calculation units and the distribution positions of the monitoring wells.
The data collected according to the purposes of the invention mainly comprise two types: firstly, monitoring well information of a water level research area to be monitored is collected, wherein the monitoring well information comprises the number, the position and the position of the monitoring well in a research area, the underground water level (buried depth) of a corresponding time node and the like. If the average water level (burial depth) change of the water level research area to be monitored is calculated, for the same monitoring well in the research area, the groundwater monitoring water level (burial depth) data at the previous moment is subtracted from the groundwater monitoring water level (burial depth) data at the later moment, and the groundwater level (burial depth) change value of each monitoring well in the research area is obtained. And secondly, collecting basic information of the water level research area to be monitored, wherein the basic information comprises the area of the water level research area to be monitored, the area of each level of area calculation unit, and the quantity and distribution position conditions of the monitoring wells in each level of area calculation unit.
Step S3), the groundwater level and groundwater level variation of each final stage area computing unit are computed according to the information of the monitoring well in the research area of the water level to be monitored and the basic information of the research area of the water level to be monitored, and the groundwater level and groundwater level variation of all the final stage area computing units are obtained, and the method comprises the following steps:
s31), the total number of the final-stage area calculation units is v, and the measured underground water level and the underground water level variation of all monitoring wells in each final-stage area calculation unit are obtained;
s32) calculating the groundwater level of the jth final stage area calculating unit according to the groundwater levels measured by all monitoring wells in the jth final stage area calculating unit
And groundwater level change amount of jth final stage area calculation unit
Groundwater level of the jth last stage region calculation unit
A groundwater level variation amount of the jth last stage region calculation unit
Underground water level of the tth monitoring well in the jth N-level area calculation unit,
calculating the groundwater level variation of the tth monitoring well in the jth N-stage area calculating unit, wherein t is more than or equal to 1 and less than or equal to N, N is the total number of the monitoring wells in the jth last-stage area calculating unit, and j is more than or equal to 1 and less than or equal to v;
s33) obtaining the groundwater level and the groundwater level variation amount of all the final stage area calculation units.
Step S4), sequentially and iteratively calculating the groundwater level and the groundwater level variation of the N-1-level region calculating unit, the groundwater level and the groundwater level variation of the N-2-level region calculating unit, …, the groundwater level and the groundwater level variation of the first-level region calculating unit according to the groundwater level and the groundwater level variation of all the last-level region calculating units, including the following steps:
s41) recording the total number of h-1 level region calculating units as w, recording the total number of h-level region calculating units in the e-1 level region calculating unit as m, wherein the k-th h-level region calculating unit in the m h-level region calculating units comprises r h +1 level region calculating units; e is more than or equal to 1 and less than or equal to w; k is more than or equal to 1 and less than or equal to m; h is more than or equal to 1 and less than or equal to N-1;
s42) acquiring the area of each h + 1-level area calculation unit according to the information of the monitoring well in the water level research area to be monitored and the basic information of the water level research area to be monitored;
s43) calculating the groundwater level of the kth level area calculating unit
Calculating the groundwater level variation of the kth h level area calculation unit
Calculating the area of a unit for the ith h + 1-level region;
s44) sequentially calculating the groundwater level and the groundwater level variation of all the h-level region calculating units to obtain a groundwater level set { D ] of the h-level region calculating unitsh,1、Dh,2、…、Dh,k、…、Dh,mSet of groundwater level variation { Δ D } and level-h regional computing unitsh,1、ΔDh,2、…、ΔDh,k、…、ΔDh,m}；
S45) calculating groundwater level of the e-th h-1 level area calculating unit according to groundwater levels and groundwater level variation of all h level area calculating units
Groundwater level variation of the e-th h-1 level area calculation unit
The area of the cell is calculated for the kth h level region.
The method comprises the steps of firstly, starting from a divided final-stage area computing unit, and obtaining the groundwater level and groundwater level variation of all final-stage area computing units by an arithmetic mean method according to monitoring data of monitoring wells in all final-stage computing units. Then, according to the groundwater level (buried depth) of each final stage region calculation unit, the variation value of the groundwater level (buried depth) and the final stage region calculation unit area which are obtained by calculation, and the subordination relation with the previous-stage area calculation unit is obtained by calculating the underground water level (buried depth) and the underground water level (buried depth) change value of each previous-stage calculation unit by adopting an area weighted average method, continuously iterating, calculating the underground water level (buried depth) and the underground water level (buried depth) change value of each previous-stage unit by adopting the area weighted average method and obtaining the underground water level (buried depth) and the underground water level (buried depth) change value of each previous-stage unit by using the underground water level (buried depth) and the underground water level (buried depth) change value of the current-stage calculation unit and the subordination relation with the previous-stage unit, and finally calculating the underground water level (buried depth) and the underground water level (buried depth) change value of the water level research area to be monitored.
In the embodiment, aiming at the urgent need of groundwater level monitoring, in order to effectively, accurately and quickly calculate the average groundwater level or groundwater level change in a region, the change of groundwater level burial depth in 2019 of a plain region of a part of northern regions is selected as a representative example to be analyzed, the area data of the administrative regions are from relevant statistics yearbooks, and the groundwater level burial depth data are from actually measured data of a hydrological department. According to three-level administrative divisions of provincial, local and county, as shown in fig. 2 and 3, a research area 1 of water level to be monitored is provided with a plurality of monitoring wells 2, and the research area 1 of water level to be monitored is divided into four secondary units (regional units) and 61 tertiary units (county units). The four secondary units (ground level units) are respectively a ground level unit a, a ground level unit b, a ground level unit c and a ground level unit d, the ground level unit d is divided into 16 county level units (see fig. 4), the ground level unit c is divided into 11 county level units, the ground level unit b is divided into 18 county level units, and the ground level unit a is divided into 16 county level units.
The present embodiment collects basic data required for consolidation calculations, including study area and monitoring well information. The area of the research area can be obtained according to the area of the administrative area issued by the relevant department; and for the data of the underground water monitoring wells, the number, the positions and the monitoring horizons of the monitoring wells in the research area and each computing unit are arranged, and the underground water level burial depth at the beginning and the end of the year and the like are arranged. When the underground water level buried depth change of each monitoring well is calculated, the underground water level buried depth change value delta dw of each monitoring well is obtained by subtracting the underground water level buried depth change at the beginning of the year from the underground water level buried depth change at the end of the yeart2-dt1Delta dw represents the underground water level buried depth change value at the end of 2019 and at the beginning of the year of a certain monitoring well, dt2Indicates the buried depth of underground water level at the end of 2019 years of the monitoring well, dt1And the underground water level burial depth of the monitoring well 2019 in the early years is shown.
The embodiment starts from county level unit (final level area calculation unit) and buries according to the underground water level of each monitoring well and each monitoring well in each county level unitAnd (4) calculating the average value of the buried depth change calculation numbers of the groundwater level of the monitoring well as the buried depth change of the county level unit groundwater. Taking the first county level unit in the district level unit d as an example (see fig. 5), the county level administrative district has 25 underground water monitoring wells in total, and the underground water level buried depth change value of the county level unit is
ΔdwiThe underground water level buried depth change value of a certain underground water monitoring well in a county level unit.
And then, calculating the groundwater level burial depth change of each level unit by adopting an area weighted average method according to the calculated groundwater level burial depth change and unit area of each county level unit and the subordination relation with the upper level unit. Taking a grade unit d as an example, 16 county-level units are totally arranged in the grade unit, and the underground water level buried depth change value of the grade unit d is calculated
ΔdS＝(ΔdX1×SX1+ΔdX2×SX2+…+ΔdX15×SX15+…+ΔdX16×SX16)/(SX1+SX2+…SX15+
SX16)，ΔdXiThe underground water level buried depth change value of the ith county level unit in the ground level unit d; sXiIs the area of the ith county level cell in the ground level cell d.
According to the underground water buried depth change and the unit area of each ground level unit obtained by calculation, calculating the underground water buried depth change value of the research area (namely, provincial level unit) by adopting an area weighted average method
Δd＝(ΔdS1×Ss1+ΔdS2×SS2+ΔdS3×Ss3+ΔdS4×SS4)/(SS1+SS2+SS3+SS4)，ΔdSiIs the groundwater level burial depth variation value of the ith ground level unit, SSiIs the area of the ith level cell.
In the embodiment, a county administrative district is used as a final-stage regional calculation unit, and in actual application, units with smaller towns, villages and villages can be selected as the final-stage calculation unit according to the data base and the working precision requirement.
The simple and convenient regional average groundwater level calculation method provided by the invention can be used for quickly and accurately calculating the regional average groundwater level and the water level change value of the regional average groundwater level relative to a certain moment according to the monitoring data of the regional groundwater level monitoring well, so that a basis and a support are provided for management work such as groundwater treatment effect evaluation and assessment.
By adopting the technical scheme disclosed by the invention, the following beneficial effects are obtained:
the method can simply and conveniently calculate the groundwater level of the region, is simple and easy to operate, needs few calculation tools and data, and greatly reduces the calculation cost compared with the methods of a kriging interpolation method, an isoline weighting method and the like in the prior art. The method has the advantages of high calculation speed, high calculation precision and more comprehensive and scientific calculation result, not only considers the current situation of underground water monitoring, but also considers the regional area distribution and the weight, has better applicability, and is not only suitable for administrative divisions such as provinces, cities, counties, villages and the like, but also suitable for water resource drainage basin divisions such as water resource first-level districts, second-level districts, third-level districts and the like.
The foregoing is only a preferred embodiment of the present invention, and it should be noted that it will be apparent to those skilled in the art that various modifications and improvements can be made without departing from the principle of the present invention, and such modifications and improvements should also be considered within the scope of the present invention.
Claims (5)
1. A method for monitoring the average groundwater level in an area is characterized by comprising the following steps:
s1) dividing the water level research area to be monitored into N-level area calculation units with different levels, wherein the Nth-level area calculation unit is a final-level area calculation unit;
s2) collecting information of a monitoring well in the research area of the water level to be monitored and basic information of the research area of the water level to be monitored;
s3) calculating the groundwater level and groundwater level variation of the final stage area calculating units according to the information of the monitoring wells in the water level research area to be monitored and the basic information of the water level research area to be monitored, and obtaining the groundwater level and groundwater level variation of all the final stage area calculating units;
s4) sequentially and iteratively calculating the groundwater level and the groundwater level variation of the N-1 level region calculating unit, the groundwater level and the groundwater level variation of the N-2 level region calculating unit, …, the groundwater level and the groundwater level variation of the first level region calculating unit according to the groundwater level and the groundwater level variation of all the last level region calculating units;
s5) obtaining the groundwater level and the groundwater level variation of the research area of the water level to be monitored.
2. The method of claim 1, wherein in the step S1), the N-level different-level region calculating units include a primary region calculating unit, the primary region calculating unit includes a plurality of secondary region calculating units, …, an N-2 level region calculating unit includes a plurality of N-1 level region calculating units, and an N-1 level region calculating unit includes a plurality of N-level region calculating units.
3. The method for monitoring the average groundwater level according to claim 1 or 2, wherein in step S2), the information of the monitoring wells in the research area of the water level to be monitored comprises the number of the monitoring wells in the research area, the positions of the monitoring wells, the levels of the monitoring wells and the groundwater levels measured by the monitoring wells at different time nodes; the basic information of the water level research area to be monitored comprises the area of the research area, the area of each level of area calculation unit, the number of monitoring wells in each level of area calculation unit and the distribution positions of the monitoring wells.
4. The method for monitoring the average groundwater level according to claim 1, wherein in step S3), the groundwater level and the groundwater level variation of each final stage area computing unit are computed according to the information of the monitoring well in the research area of the water level to be monitored and the basic information of the research area of the water level to be monitored, and the groundwater level variation of all the final stage area computing units are obtained, comprising the steps of:
s31), obtaining the groundwater level and groundwater level variation quantity measured by all monitoring wells in each final stage area calculation unit, wherein the total number of the final stage area calculation units is v;
s32) calculating the groundwater level of the jth final stage area calculating unit according to the groundwater levels measured by all monitoring wells in the jth final stage area calculating unit
And groundwater level variation amount of jth last-stage area calculation unit
Groundwater level of the jth last stage region calculation unit
A groundwater level variation amount of the jth last stage region calculation unit
Underground water level of the tth monitoring well in the jth N-level area calculation unit,
calculating the groundwater level variation of the tth monitoring well in the jth N-stage area calculating unit, wherein t is more than or equal to 1 and less than or equal to N, N is the total number of the monitoring wells in the jth last-stage area calculating unit, and j is more than or equal to 1 and less than or equal to v;
s33) obtaining the groundwater level and the groundwater level variation amount of all the final stage area calculation units.
5. The method of claim 4, wherein the step S4) of sequentially and iteratively calculating the groundwater level and the groundwater level variation of the N-1 stage region calculation unit, the groundwater level and the groundwater level variation of the N-2 stage region calculation unit, …, the groundwater level and the groundwater level variation of the one stage region calculation unit according to the groundwater level and the groundwater level variation of all the last stage region calculation units comprises the steps of:
s41) recording the total number of h-1 level region calculating units as w, recording the total number of h-level region calculating units in the e-1 level region calculating unit as m, wherein the k-th h-level region calculating unit in the m h-level region calculating units comprises r h +1 level region calculating units; e is more than or equal to 1 and less than or equal to w; k is more than or equal to 1 and less than or equal to m; h is more than or equal to 1 and less than or equal to N-1;
s42) acquiring the area of each h + 1-level area calculation unit according to the information of the monitoring well in the water level research area to be monitored and the basic information of the water level research area to be monitored;
s43) calculating the groundwater level of the kth level area calculating unit
Calculating the groundwater level variation of the kth h level area calculation unit
Calculating the area of a unit for the ith h + 1-level region;
s44) sequentially calculating the groundwater level and the groundwater level variation of all the h-level region calculating units to obtain a groundwater level set { D ] of the h-level region calculating unitsh,1、Dh,2、…、Dh,k、…、Dh,mSet of groundwater level variation { Δ D } and level-h regional computing unitsh,1、ΔDh,2、…、ΔDh,k、…、ΔDh,m}；
S45) calculating the e-th h-1-level area according to the groundwater levels and the groundwater level variation of all the h-level area calculating unitsGroundwater level of a computing unit
Groundwater level variation of e-th h-1 level area calculation unit
The area of the cell is calculated for the kth h level region.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011186789.6A CN112268595A (en) | 2020-10-29 | 2020-10-29 | Method for monitoring regional average groundwater level |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202011186789.6A CN112268595A (en) | 2020-10-29 | 2020-10-29 | Method for monitoring regional average groundwater level |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN112268595A true CN112268595A (en) | 2021-01-26 |
Family
ID=74345237
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202011186789.6A Pending CN112268595A (en) | 2020-10-29 | 2020-10-29 | Method for monitoring regional average groundwater level |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN112268595A (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113780882A (en) * | 2021-09-26 | 2021-12-10 | 中国水利水电科学研究院 | Groundwater level amplitude evaluation system after groundwater super-mining treatment |
| CN114593792A (en) * | 2022-03-29 | 2022-06-07 | 中国水利水电科学研究院 | Underground water level monitoring method and device and storage medium |
| CN115687883A (en) * | 2022-11-18 | 2023-02-03 | 水利部信息中心 | Shallow groundwater storage variable calculation method based on kilometer grids |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2008118249A (en) * | 2008-05-07 | 2009-11-20 | Государственное образовательное учреждение высшего профессионального образования "Самарский государственный университет путей сооб | METHOD FOR DETERMINING MEASUREMENT OF GROUND WATER LEVEL |
| CN107038258A (en) * | 2017-05-18 | 2017-08-11 | 中国地质环境监测院 | Groundwater monitoring data acquisition release management system |
| CN110019604A (en) * | 2017-10-31 | 2019-07-16 | 于盛 | A kind of groundwater monitoring early warning system |
| CN110031596A (en) * | 2019-03-14 | 2019-07-19 | 中国科学院新疆生态与地理研究所 | A kind of real-time dynamic monitoring method for inland river basin in arid zone underground water |
-
2020
- 2020-10-29 CN CN202011186789.6A patent/CN112268595A/en active Pending
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| RU2008118249A (en) * | 2008-05-07 | 2009-11-20 | Государственное образовательное учреждение высшего профессионального образования "Самарский государственный университет путей сооб | METHOD FOR DETERMINING MEASUREMENT OF GROUND WATER LEVEL |
| CN107038258A (en) * | 2017-05-18 | 2017-08-11 | 中国地质环境监测院 | Groundwater monitoring data acquisition release management system |
| CN110019604A (en) * | 2017-10-31 | 2019-07-16 | 于盛 | A kind of groundwater monitoring early warning system |
| CN110031596A (en) * | 2019-03-14 | 2019-07-19 | 中国科学院新疆生态与地理研究所 | A kind of real-time dynamic monitoring method for inland river basin in arid zone underground water |
Non-Patent Citations (1)
| Title |
|---|
| 廖梓龙等: "地下水控制性临界水位及确定方法", 《水利水电技术》 * |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113780882A (en) * | 2021-09-26 | 2021-12-10 | 中国水利水电科学研究院 | Groundwater level amplitude evaluation system after groundwater super-mining treatment |
| CN114593792A (en) * | 2022-03-29 | 2022-06-07 | 中国水利水电科学研究院 | Underground water level monitoring method and device and storage medium |
| CN115687883A (en) * | 2022-11-18 | 2023-02-03 | 水利部信息中心 | Shallow groundwater storage variable calculation method based on kilometer grids |
| CN115687883B (en) * | 2022-11-18 | 2024-04-30 | 水利部信息中心 | Shallow groundwater storage variable calculation method based on kilometer grid |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN112268595A (en) | Method for monitoring regional average groundwater level | |
| Singh et al. | Analysis of drivers of trends in groundwater levels under rice–wheat ecosystem in Haryana, India | |
| CN110232737B (en) | Urban catchment area division method | |
| Sayl et al. | Optimization of area–volume–elevation curve using GIS–SRTM method for rainwater harvesting in arid areas | |
| Kumar et al. | Sensitivity analysis of the GIUH based Clark model for a catchment | |
| Nagarajan et al. | Remote sensing-and GIS-based runoff modeling with the effect of land-use changes (a case study of Cochin corporation) | |
| Dadhich et al. | Assessment and prediction of groundwater using geospatial and ANN modeling | |
| CN113240257A (en) | Territorial space partitioning method and device based on minimum cumulative resistance model | |
| Sahu et al. | Density-based spatial clustering of application with noise approach for regionalisation and its effect on hierarchical clustering | |
| Peterson et al. | Joint estimation of gross recharge, groundwater usage, and hydraulic properties within HydroSight | |
| Hwang et al. | Quantifying dam‐induced fluctuations in streamflow frequencies across the Colorado River Basin | |
| Barthel et al. | Integrated assessment of groundwater resources in the Oueme basin, Benin, West Africa | |
| Suryanarayana et al. | Groundwater-level assessment and prediction using realistic pumping and recharge rates for semi-arid coastal regions: a case study of Visakhapatnam city, India. | |
| Hagras | Runoff modeling using SCS-CN and GIS approach in the Tayiba Valley Basin, Abu Zenima area, South-west Sinai, Egypt | |
| Xue et al. | Determination of hydraulic conductivity and its spatial variability in the Jianghan Plain using a multi-format, multi-method approach | |
| Michalek et al. | Evaluation of CMIP6 HighResMIP for hydrologic modeling of annual maximum discharge in Iowa | |
| Zhang et al. | Coastline changes in mainland China from 2000 to 2015 | |
| Fang et al. | Non-identical models for seasonal flood frequency analysis | |
| Wright et al. | Identification of hydrologic models, optimized parameters, and rainfall inputs consistent with in situ streamflow and rainfall and remotely sensed soil moisture | |
| CN103150616A (en) | Digging method for distribution characteristics of prehistoric settlement site | |
| Wang et al. | Probabilistic modeling of the rainfall severity and height for locating the surface artificial recharge structure | |
| CN111369102B (en) | Method and device for extracting waterlogging risk points | |
| Gayen | Field-based spatio-temporal monitoring of hydrograph network stations to predict the long-term behavioral pattern of groundwater regime and its implications in India: A review | |
| Singh et al. | Hydrologic modeling with GIS | |
| Cloete et al. | Revision of regional maximum flood (RMF) estimation in Namibia |
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 | ||
| RJ01 | Rejection of invention patent application after publication | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20210126 |