CN116244958A - Method for researching mining and supplementing balance of seawater invasion underground water in river basin - Google Patents
Method for researching mining and supplementing balance of seawater invasion underground water in river basin Download PDFInfo
- Publication number
- CN116244958A CN116244958A CN202310257229.2A CN202310257229A CN116244958A CN 116244958 A CN116244958 A CN 116244958A CN 202310257229 A CN202310257229 A CN 202310257229A CN 116244958 A CN116244958 A CN 116244958A
- Authority
- CN
- China
- Prior art keywords
- water
- river basin
- river
- flow
- groundwater
- 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
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 105
- 239000013535 sea water Substances 0.000 title claims abstract description 42
- 230000009545 invasion Effects 0.000 title claims abstract description 21
- 238000000034 method Methods 0.000 title claims abstract description 19
- 238000005065 mining Methods 0.000 title claims description 13
- 230000001502 supplementing effect Effects 0.000 title claims description 10
- 239000003673 groundwater Substances 0.000 claims abstract description 49
- 238000001556 precipitation Methods 0.000 claims abstract description 21
- 230000008595 infiltration Effects 0.000 claims description 24
- 238000001764 infiltration Methods 0.000 claims description 24
- 238000001704 evaporation Methods 0.000 claims description 17
- 230000008020 evaporation Effects 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 15
- 230000035699 permeability Effects 0.000 claims description 15
- 230000009189 diving Effects 0.000 claims description 9
- 238000003973 irrigation Methods 0.000 claims description 9
- 230000002262 irrigation Effects 0.000 claims description 9
- 239000005442 atmospheric precipitation Substances 0.000 claims description 6
- 238000009933 burial Methods 0.000 claims description 6
- 239000002609 medium Substances 0.000 claims description 6
- 239000012736 aqueous medium Substances 0.000 claims description 3
- 238000007599 discharging Methods 0.000 claims description 3
- 238000009826 distribution Methods 0.000 claims description 3
- 230000005484 gravity Effects 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 claims description 3
- 238000005325 percolation Methods 0.000 claims description 3
- 238000003860 storage Methods 0.000 claims description 3
- 230000000007 visual effect Effects 0.000 claims description 3
- 230000029142 excretion Effects 0.000 claims description 2
- 239000010865 sewage Substances 0.000 claims 1
- 230000002265 prevention Effects 0.000 abstract description 6
- 239000013505 freshwater Substances 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000006378 damage Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000002939 deleterious effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- 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
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
-
- 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—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Entrepreneurship & Innovation (AREA)
- Marketing (AREA)
- General Business, Economics & Management (AREA)
- Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Evolutionary Computation (AREA)
- Computer Hardware Design (AREA)
- Primary Health Care (AREA)
- Geometry (AREA)
- Educational Administration (AREA)
- General Health & Medical Sciences (AREA)
- Water Supply & Treatment (AREA)
- Public Health (AREA)
- General Engineering & Computer Science (AREA)
- Sewage (AREA)
Abstract
The invention discloses a method for researching the balance of seawater invasion groundwater in a river basin, which comprises the following steps: s1: calculating the sink-source item coefficient in the stream domain; s2: defining a research scope and boundaries in a flow domain, and establishing a hydrogeologic conceptual model in the flow domain: s3: and inputting the sink source item coefficients in the river basin into the hydrogeologic conceptual model to obtain underground water balance data in the river basin, and evaluating whether the river basin is invaded by seawater according to the underground water balance data. According to the method, the relation between the submerged depth and the precipitation change is researched by analyzing the dynamic change trend of the groundwater level and the precipitation of a research area in a flow area, and the influence of the precipitation on the submerged depth is analyzed; by establishing a hydrogeologic conceptual model, an effective method is provided for researching and analyzing seawater invasion prevention and control measures at home and abroad, and effective data is provided for providing seawater invasion prevention and control measures suitable for a research area.
Description
Technical Field
The invention relates to the technical field of seawater invasion research, in particular to a method for researching balance of seawater invasion groundwater mining and compensation in a river basin.
Background
Seawater invasion refers to a deleterious hydrogeologic effect of seawater moving toward the continental aquifer due to the destruction of the natural equilibrium conditions of seawater and groundwater by excess groundwater extraction in the seashore area. From the above definition, the nature of seawater intrusion is: the coastal area is characterized in that the groundwater is greatly reduced due to the fact that groundwater is artificially and excessively extracted, the original dynamic balance between seawater and fresh water is destroyed, and therefore a salty-fresh water interface is pushed to the land. In recent years, the development and utilization of groundwater resources on a large scale are started, the exploitation amount is greatly increased, and the development and utilization of groundwater play a great supporting and guaranteeing role for the rapid development of economy and society. Meanwhile, the problem of super mining caused by unreasonable development and utilization of underground water is also highlighted, and sustainable utilization of water resources and sustainable development of economy and society are seriously threatened. Therefore, the groundwater super-mining area is controlled and treated in time, the water resource management work of the groundwater super-mining area is enhanced, the seawater invasion condition of the groundwater in the coastal area needs to be researched, and a basis is provided for the groundwater fresh water treatment.
Disclosure of Invention
Aiming at the defects in the prior art, the invention provides a research method for the balance of the groundwater production and the compensation of the seawater invasion in the river basin, which provides basis for the seawater invasion in the river basin.
In order to achieve the aim of the invention, the invention adopts the following technical scheme:
the method for researching the balance of seawater invasion groundwater mining and supplementing in the river basin comprises the following steps:
s1: calculating a sink item coefficient in a flow field, wherein the sink item coefficient comprises atmospheric precipitation infiltration supply quantity, river leakage supply quantity, agricultural irrigation infiltration quantity, reservoir infiltration quantity, diving evaporation quantity, underground water exploitation quantity and lateral outflow quantity;
s2: defining a research scope and boundaries in a flow domain, and establishing a hydrogeologic conceptual model in the flow domain:
wherein omega is a seepage area, epsilon is a source and sink coefficient of an aquifer, h is a groundwater level, h 0 K is the initial water level distribution of the underground water system x 、K y Horizontal permeability coefficients, K, in the x-direction and y-direction, respectively, of the aqueous medium z Is the vertical z permeability coefficient of the water-containing medium, S is the water storage rate of the water-containing layer below the free surface, mu is the gravity water supply rate of the submerged water-containing layer, Γ 0 Is the upper boundary of the seepage area, namely the free surface of the groundwater, p is the evaporation and precipitation seepage intensity of the diving surface, h 1 For a known boundary water level value Γ 1 Is the known water level boundary Γ 2 As a flow boundary of the percolation region,
is the normal direction of the boundary surface, +.>
The permeability coefficient in the normal direction of the boundary surface is q is the flow rate in the unit area of the flow rate boundary, the inflow is positive, the outflow is negative, the water is isolated, and t is the study time;
s3: and inputting the sink source item coefficients in the river basin into the hydrogeologic conceptual model to obtain underground water balance data in the river basin, and evaluating whether the river basin is invaded by seawater according to the underground water balance data.
Further, the method for calculating the atmospheric precipitation infiltration replenishment quantity in the step S1 comprises the following steps:
wherein Q is i k Supplementing the amount of precipitation infiltration of the kth month of the position i in the river basin, alpha i For precipitation infiltration coefficient corresponding to position i, P i k Precipitation for the kth month of position i, F i The calculated area for position i;
river leakage supply quantity Q R The calculation method of (1) is as follows:
wherein h is r Is the water level of the river, h is the groundwater level, W is the river width, M is the thickness of the bed bottom layer, K s Is the permeability coefficient of the bed bottom layer, L is the river length, C R Hydraulic conductivity coefficient for river seepage;
the infiltration amount of agricultural irrigation is 50% of the average agricultural irrigation amount in the river basin;
reservoir infiltration quantity Q Reservoir The calculation method of (1) is as follows: q (Q) Reservoir =v·f·Δt, where V is the seepage velocity of the reservoir in the flow, F is the area of the blowby region of the reservoir in the flow, Δt is the seepage time of the reservoir in the flow;
the calculation method of the diving evaporation capacity comprises the following steps:
wherein E is g For the evaporation intensity of groundwater E 0 The evaporation intensity of the water surface, d 0 The evaporation limit burial depth of the underground water is d, the water level burial depth of the underground water is n, the empirical coefficient is n=1;
the underground water exploitation quantity comprises agricultural exploitation quantity, life exploitation quantity and town self-provided well exploitation quantity, and the agricultural exploitation quantity, the life exploitation quantity and the town self-provided well exploitation quantity are obtained by counting the agricultural exploitation quantity, the life exploitation quantity and the town self-provided well exploitation quantity in a flow area;
lateral outflow Q b The calculation method of (1) is as follows:
wherein C is b Is the hydraulic conductivity coefficient of an external water source and a river basin, h b H' is the water level of the external water source and is the water level in the flow field; l x W is the cross-sectional area of the outside water source facing the basin, K is the permeability coefficient of the medium between the outside water source and the basin, and D is the distance between the outside water source and the basin.
Further, step S3 includes:
s31: using visual software GMS as a platform for simulating underground water flow in a river basin, and simulating the underground water flow in the river basin according to a hydrogeological conceptual model;
s32: dividing the river basin into a plurality of rows and columns of regular rectangular grids on a plane according to the aquifer structure, boundary conditions and underground water flow field characteristics in the river basin;
s33: substituting the sink item coefficients into the hydrogeologic conceptual model for calculation to obtain a groundwater balance table in the river basin, and calculating the balance value by using the total supply quantity and the total excretion quantity of the river basin: balance = total drainage-total replenishment; comparing the equalization amount with an equalization amount threshold:
if the equilibrium amount is larger than the equilibrium amount threshold value, determining excessive consumption of the groundwater in the flow area, wherein a pressure difference is formed between the groundwater and the seawater, and the groundwater is invaded by the seawater;
if the equilibrium amount is less than or equal to the equilibrium amount threshold value, determining that the underground water in the river basin is sufficient, forming pressure difference between the underground water and the seawater, and discharging the excessive underground water into the sea, wherein the river basin is not invaded by the seawater.
The beneficial effects of the invention are as follows: according to the method, the relation between the submerged depth and the precipitation change is researched by analyzing the dynamic change trend of the groundwater level and the precipitation of a research area in a flow area, and the influence of the precipitation on the submerged depth is analyzed; by establishing a hydrogeological conceptual model, researching and analyzing the groundwater mining and supplementing balance of a research area by utilizing a three-dimensional groundwater flow model; based on the underground water flow numerical model, the development trend of seawater invasion under different mining conditions can be effectively predicted, an effective method is provided for research and analysis of seawater invasion prevention and control measures at home and abroad, and effective data is provided for providing proper seawater invasion prevention and control measures in a research area.
Drawings
Fig. 1 is a schematic diagram of a horizontal mesh subdivision of a sow river basin.
Fig. 2 is a schematic view of a vertical mesh subdivision of a sow river basin.
Detailed Description
The following description of the embodiments of the present invention is provided to facilitate understanding of the present invention by those skilled in the art, but it should be understood that the present invention is not limited to the scope of the embodiments, and all the inventions which make use of the inventive concept are protected by the spirit and scope of the present invention as defined and defined in the appended claims to those skilled in the art.
The method for researching the balance of seawater invasion groundwater mining and supplementing in the river basin comprises the following steps:
s1: and calculating a sink item coefficient in the flow field, wherein the sink item coefficient comprises atmospheric precipitation infiltration supply quantity, river leakage supply quantity, agricultural irrigation infiltration quantity, reservoir infiltration quantity, diving evaporation quantity, underground water exploitation quantity and lateral outflow quantity.
The method for calculating the atmospheric precipitation infiltration replenishment quantity comprises the following steps:
wherein Q is i k Supplementing the amount of precipitation infiltration of the kth month of the position i in the river basin, alpha i For precipitation infiltration coefficient corresponding to position i, P i k Precipitation for the kth month of position i, F i The calculated area for location i.
River leakage supply quantity Q R The calculation method of (1) is as follows:
wherein h is r Is the water level of the river, h is the groundwater level, W is the river width, M is the thickness of the bed bottom layer, K s Is the permeability coefficient of the bed bottom layer, L is the river length, C R Hydraulic conductivity coefficient for river seepage;
the infiltration amount of agricultural irrigation is 50% of the average agricultural irrigation amount in the river basin;
reservoir infiltration quantity Q Reservoir The calculation method of (1) is as follows: q (Q) Reservoir =v·f·Δt, where V is the seepage velocity of the reservoir in the flow, F is the area of the blowby region of the reservoir in the flow, Δt is the seepage time of the reservoir in the flow;
the calculation method of the diving evaporation capacity comprises the following steps:
wherein E is g For the evaporation intensity of groundwater E 0 The evaporation intensity of the water surface, d 0 The evaporation limit burial depth of the underground water is d, the water level burial depth of the underground water is n, the empirical coefficient is n=1;
the underground water exploitation quantity comprises agricultural exploitation quantity, life exploitation quantity and town self-provided well exploitation quantity, and the agricultural exploitation quantity, the life exploitation quantity and the town self-provided well exploitation quantity are obtained by counting the agricultural exploitation quantity, the life exploitation quantity and the town self-provided well exploitation quantity in a flow area;
lateral outflow Q b The calculation method of (1) is as follows:
wherein C is b Is the hydraulic conductivity coefficient of an external water source and a river basin, h b H' is the water level of the external water source and is the water level in the flow field; l is L x W is the horizontal direction facing the external water source in the river basinThe cross-sectional area, K, is the permeability coefficient of the medium between the external water source and the basin, and D is the distance between the external water source and the basin.
S2: defining a research scope and boundaries in a flow domain, and establishing a hydrogeologic conceptual model in the flow domain:
wherein omega is a seepage area, epsilon is a source and sink coefficient of an aquifer, h is a groundwater level, h 0 K is the initial water level distribution of the underground water system x 、K y Horizontal permeability coefficients, K, in the x-direction and y-direction, respectively, of the aqueous medium z Is the vertical z permeability coefficient of the water-containing medium, S is the water storage rate of the water-containing layer below the free surface, mu is the gravity water supply rate of the submerged water-containing layer, Γ 0 Is the upper boundary of the seepage area, namely the free surface of the groundwater, p is the evaporation and precipitation seepage intensity of the diving surface, h 1 For a known boundary water level value Γ 1 Is the known water level boundary Γ 2 As a flow boundary of the percolation region,
is the normal direction of the boundary surface, +.>
The permeability coefficient in the normal direction of the boundary surface is q is the flow rate in the unit area of the flow rate boundary, the inflow is positive, the outflow is negative, the water is isolated, and t is the study time; />
S3: and inputting the sink source item coefficients in the river basin into the hydrogeologic conceptual model to obtain underground water balance data in the river basin, and evaluating whether the river basin is invaded by seawater according to the underground water balance data.
The step S3 comprises the following steps:
s31: using visual software GMS as a platform for simulating underground water flow in a river basin, and simulating the underground water flow in the river basin according to a hydrogeological conceptual model;
in this embodiment, taking a sow river basin in a peninsula in shandong as an example, as shown in fig. 1, a horizontal mesh subdivision schematic diagram of the sow river basin is simulated.
S32: dividing the river basin into a plurality of rows and columns of regular rectangular grids on a plane according to the aquifer structure, boundary conditions and underground water flow field characteristics in the river basin; as shown in fig. 2, a schematic view of the vertical mesh subdivision of a sow river basin is shown.
S33: substituting the sink source term coefficients into the hydrogeologic conceptual model for calculation to obtain a groundwater balance table in the flow domain, as shown in the following table 1,
table 1 sow river basin 2014 groundwater year balance table
The equilibrium amount is calculated by using the total supply amount and the total discharge amount of the basin: balance = total drainage-total replenishment; comparing the equalization amount with an equalization amount threshold:
if the equilibrium amount is larger than the equilibrium amount threshold value, determining excessive consumption of the groundwater in the flow area, wherein a pressure difference is formed between the groundwater and the seawater, and the groundwater is invaded by the seawater;
if the equilibrium amount is less than or equal to the equilibrium amount threshold value, determining that the underground water in the river basin is sufficient, forming pressure difference between the underground water and the seawater, and discharging the excessive underground water into the sea, wherein the river basin is not invaded by the seawater.
According to the method, the relation between the submerged depth and the precipitation change is researched by analyzing the dynamic change trend of the groundwater level and the precipitation of a research area in a flow area, and the influence of the precipitation on the submerged depth is analyzed; by establishing a hydrogeological conceptual model, researching and analyzing the groundwater mining and supplementing balance of a research area by utilizing a three-dimensional groundwater flow model; based on the underground water flow numerical model, the development trend of seawater invasion under different mining conditions can be effectively predicted, an effective method is provided for research and analysis of seawater invasion prevention and control measures at home and abroad, and effective data is provided for providing proper seawater invasion prevention and control measures in a research area.
Claims (3)
1. A method for researching the balance of seawater invasion groundwater mining and supplementing in a river basin is characterized by comprising the following steps:
s1: calculating a sink item coefficient in a flow field, wherein the sink item coefficient comprises atmospheric precipitation infiltration supply quantity, river leakage supply quantity, agricultural irrigation infiltration quantity, reservoir infiltration quantity, diving evaporation quantity, underground water exploitation quantity and lateral outflow quantity;
s2: defining a research scope and boundaries in a flow domain, and establishing a hydrogeologic conceptual model in the flow domain:
wherein omega is a seepage area, epsilon is a source and sink coefficient of an aquifer, h is a groundwater level, h 0 K is the initial water level distribution of the underground water system x 、K y Horizontal permeability coefficients, K, in the x-direction and y-direction, respectively, of the aqueous medium z Is the vertical z permeability coefficient of the water-containing medium, S is the water storage rate of the water-containing layer below the free surface, mu is the gravity water supply rate of the submerged water-containing layer, Γ 0 Is the upper boundary of the seepage area, namely the free surface of the groundwater, p is the evaporation and precipitation seepage intensity of the diving surface, h 1 For a known boundary water level value Γ 1 Is the known water level boundary Γ 2 As a flow boundary of the percolation region,
is the normal direction of the boundary surface, +.>
The permeability coefficient in the normal direction of the boundary surface is q is the flow rate in the unit area of the flow rate boundary, the inflow is positive, the outflow is negative, the water is isolated, and t is the study time;
s3: and inputting the sink source item coefficients in the river basin into the hydrogeologic conceptual model to obtain underground water balance data in the river basin, and evaluating whether the river basin is invaded by seawater according to the underground water balance data.
2. The method for balancing seawater intrusion into groundwater production and compensation according to claim 1, wherein the method for calculating the infiltration replenishment amount of atmospheric precipitation in step S1 is as follows:
wherein Q is i k Supplementing the amount of precipitation infiltration of the kth month of the position i in the river basin, alpha i For precipitation infiltration coefficient corresponding to position i, P i k Precipitation for the kth month of position i, F i The calculated area for position i;
the river seepage supply quantity Q R The calculation method of (1) is as follows:
wherein h is r Is the water level of the river, h is the groundwater level, W is the river width, M is the thickness of the bed bottom layer, K s Is the permeability coefficient of the bed bottom layer, L is the river length, C R Hydraulic conductivity coefficient for river seepage;
the infiltration amount of the agricultural irrigation is 50% of the average agricultural irrigation amount in the river basin;
the reservoir infiltration quantity Q Reservoir The calculation method of (1) is as follows: q (Q) Reservoir =v·f·Δt, where V is the seepage velocity of the reservoir in the flow, F is the area of the blowby region of the reservoir in the flow, Δt is the seepage time of the reservoir in the flow;
the calculation method of the diving evaporation capacity comprises the following steps:
wherein E is g Is the groundEvaporation intensity of sewage, E 0 The evaporation intensity of the water surface, d 0 The evaporation limit burial depth of the underground water is d, the water level burial depth of the underground water is n, the empirical coefficient is n=1;
the underground water exploitation amount comprises agricultural exploitation amount, life exploitation amount and town self-provided well exploitation amount, and the agricultural exploitation amount, the life exploitation amount and the town self-provided well exploitation amount are obtained by counting the agricultural exploitation amount, the life exploitation amount and the town self-provided well exploitation amount in a flow area;
the lateral outflow Q b The calculation method of (1) is as follows:
wherein C is b Is the hydraulic conductivity coefficient of an external water source and a river basin, h b H' is the water level of the external water source and is the water level in the flow field; l x W is the cross-sectional area of the outside water source facing the basin, K is the permeability coefficient of the medium between the outside water source and the basin, and D is the distance between the outside water source and the basin.
3. The method for studying the balance of seawater invasion groundwater production and compensation in a river basin according to claim 1, wherein the step S3 comprises:
s31: using visual software GMS as a platform for simulating underground water flow in a river basin, and simulating the underground water flow in the river basin according to a hydrogeological conceptual model;
s32: dividing the river basin into a plurality of rows and columns of regular rectangular grids on a plane according to the aquifer structure, boundary conditions and underground water flow field characteristics in the river basin;
s33: substituting the sink item coefficients into the hydrogeologic conceptual model for calculation to obtain a groundwater balance table in the river basin, and calculating the balance value by using the total supply quantity and the total excretion quantity of the river basin: balance = total drainage-total replenishment; comparing the equalization amount with an equalization amount threshold:
if the equilibrium amount is larger than the equilibrium amount threshold value, determining excessive consumption of the groundwater in the flow area, wherein a pressure difference is formed between the groundwater and the seawater, and the groundwater is invaded by the seawater;
if the equilibrium amount is less than or equal to the equilibrium amount threshold value, determining that the underground water in the river basin is sufficient, forming pressure difference between the underground water and the seawater, and discharging the excessive underground water into the sea, wherein the river basin is not invaded by the seawater.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310257229.2A CN116244958A (en) | 2023-03-13 | 2023-03-13 | Method for researching mining and supplementing balance of seawater invasion underground water in river basin |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310257229.2A CN116244958A (en) | 2023-03-13 | 2023-03-13 | Method for researching mining and supplementing balance of seawater invasion underground water in river basin |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116244958A true CN116244958A (en) | 2023-06-09 |
Family
ID=86626033
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310257229.2A Pending CN116244958A (en) | 2023-03-13 | 2023-03-13 | Method for researching mining and supplementing balance of seawater invasion underground water in river basin |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116244958A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116993030A (en) * | 2023-09-27 | 2023-11-03 | 长江水利委员会水文局 | Reservoir pressure salty taste adjustment method and system under variable conditions |
| CN118229040A (en) * | 2024-05-24 | 2024-06-21 | 水利部水利水电规划设计总院 | River channel infiltration compensation evaluation management method based on water balance |
-
2023
- 2023-03-13 CN CN202310257229.2A patent/CN116244958A/en active Pending
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116993030A (en) * | 2023-09-27 | 2023-11-03 | 长江水利委员会水文局 | Reservoir pressure salty taste adjustment method and system under variable conditions |
| CN116993030B (en) * | 2023-09-27 | 2023-12-08 | 长江水利委员会水文局 | Reservoir pressure salty taste adjustment method and system under variable conditions |
| CN118229040A (en) * | 2024-05-24 | 2024-06-21 | 水利部水利水电规划设计总院 | River channel infiltration compensation evaluation management method based on water balance |
| CN118229040B (en) * | 2024-05-24 | 2024-08-06 | 水利部水利水电规划设计总院 | River channel infiltration compensation evaluation management method based on water balance |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116244958A (en) | Method for researching mining and supplementing balance of seawater invasion underground water in river basin | |
| Bobba | Numerical modelling of salt-water intrusion due to human activities and sea-level change in the Godavari Delta, India | |
| Oude Essink | Salt water intrusion in a three-dimensional groundwater system in the Netherlands: a numerical study | |
| Vandenbruwaene et al. | Sedimentation and response to sea-level rise of a restored marsh with reduced tidal exchange: Comparison with a natural tidal marsh | |
| Da Lio et al. | Saltwater contamination in the managed low-lying farmland of the Venice coast, Italy: An assessment of vulnerability | |
| CN108446489A (en) | Measurement method and its processing unit of the Coastline Changes to wetland influence on groundwater | |
| Charlesworth et al. | The Burdekin Delta—Australia’s oldest artificial recharge scheme | |
| Zhang et al. | An overview of ecohydrology of the Yellow River delta wetland | |
| Luternauer et al. | Salt marshes | |
| Nofal et al. | Groundwater: Extraction versus recharge; vulnerability assessment | |
| Geriesh et al. | Geoenvironmental impact assessment of El-Salam Canal on the surrounding soil and groundwater flow regime, Northwestern Sinai, Egypt | |
| CN116226994A (en) | Underground dam height optimization method for preventing seawater intrusion and improving self-cleaning of aquifer | |
| Sana et al. | Feasibility study of using treated wastewater to mitigate seawater intrusion along northern coast of Oman | |
| Zhang et al. | Spatial distribution characteristics of soil water-salt gradients in the ecological buffer zone of arid zone lakes and their influencing factors | |
| Cross et al. | The East Kent Chalk Aquifer during the 1988‐92 Drought | |
| Don et al. | Groundwater resources and management for paddy field irrigation and associated environmental problems in an alluvial coastal lowland plain | |
| Hamdan et al. | Estimation of seepage loss from Nasser Lake to the adjacent Nubian sandstone aquifer, southern Egypt | |
| Subetto et al. | Palaeorunoff from Lake Ladoga to the Baltic Sea during the Holocene | |
| Saliev et al. | Analysis of changes in moisture transport parameters in soils under waterlogged conditions | |
| CN114215000B (en) | Wetland ecological restoration method based on water system reconstruction | |
| Khayyun et al. | Modelling of Groundwater Quality of Tigris River Reach-in Baghdad-Iraq Using Groundwater Modeling System Software | |
| Tsapleva et al. | Construction of a mathematical model for the extraction of mineral raw materials | |
| Lecca et al. | Modeling seawater intrusion in the Korba aquifer (Tunisia) | |
| Kelleners | Effluent salinity of pipe drains and tube-wells: a case study from the Indus plain | |
| Masterson et al. | Hydrogeology and Ground-water Resources of the Coastal Aquifers of Southeastern Massachusetts |
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 |