Disclosure of Invention
Aiming at the serious problem of underground water super-mining in North China, the nation provides comprehensive treatment measures of 'one reduction and one increase', and the significance of saving water (namely 'one reduction') and increasing multi-channel water source supply (namely 'one increase') as important ways for realizing underground water pressure mining is great. Considering the spatial imbalance of the regional underground water pressure mining target and the accessibility of measures such as water conservation and water increase, the invention provides a system for analyzing the contribution rate balance of the alternative water source based on the underground water pressure mining target.
In the existing balance analysis of the contribution rate of the alternative water source, the defects of contribution rate of water saving and water increasing to underground water pressure mining cannot be described at the same time, and the invention takes water shortage as a virtual water source for balancing the gap between the water saving, water increasing and underground water supply variation; by improving an optimized alternative water source contribution rate balance analysis system, the alternative water source replacement relation (the replacement relation is described by the contribution rate) between the water saving amount, the water increasing amount, the water shortage amount and the groundwater supply amount (the groundwater supply amount is the difference between the actual mining amount and the groundwater pressure mining target) under different scenes in the region is quantitatively analyzed, and reference is provided for water resource management work such as water distribution.
The invention provides a system for analyzing the contribution rate balance of a substitute water source based on an underground water pressure mining target, which comprises:
the first module is used for predicting the water demand corresponding to each combination scheme and determining an underground water pressure mining target;
the combination scheme comprises a scheme 1, a scheme 2 and a scheme 3, and is defined as follows:
scheme 1 refers to: the current annual underground water exploitation amount is taken as water supply constraint, and the general water saving is taken as water demand input;
scheme 2 refers to: the current annual underground water exploitation amount is used as water supply restriction, and the reinforced water saving is used as water demand input;
scheme 3 refers to: the situation scheme that the underground water mining amount based on the underground water pressure mining target is taken as water supply restriction, and the strengthened water saving is taken as water demand input;
the second module is used for receiving the water demand and the underground water pressure collection target output by the first module and obtaining the water supply amount and the water shortage amount corresponding to each optimized combination scheme by using the GAMS optimization module;
the GAMS optimization module calculates the water shortage by adopting an objective function shown in formula (1), constraint conditions shown in formulas (2) to (3) and formula (4):
Gk,i≤LIMITk,i (2)
wherein t represents a time period number, and tl is the total time period number; m represents the number of the alternative water source types, and ml represents the total number of the alternative water source types; i represents different area units, and il is the total number of the area units; k denotes a combination scheme number;
δi,m,trepresenting the marginal benefit of the water supply of the ith unit mth type water source in the tth time period; sk,i,mIndicates the water supply amount of the m-th type water source of the i-th zone unit in the k-th scheme, Gk,iIndicating groundwater feed of the i-th zone unit in the k-th scheme; LIMITk,iRepresenting the groundwater mining target of the ith zone unit in the kth scheme; qk,iRepresenting the water demand of the unit of the ith area in the kth scheme; n is a radical ofk,iIndicating the water shortage of the unit of the ith area in the kth scheme; sk,i,m、Gk,i、Nk,iOptimizing the output of the module for GAMS;
the third module is used for receiving the water supply amount and the water shortage amount output by the second module, judging whether the regional unit lacks water under each combination scheme, and calculating the contribution rate of the regional unit according to the water shortage situation;
when the regional unit does not have a water shortage scene, respectively calculating the water saving contribution rate and the multi-source water increasing contribution rate by adopting an equation (5) and an equation (6); when the regional unit has a water shortage scene, calculating a water saving contribution rate, a multi-source water increasing contribution rate and a water shortage contribution rate by adopting an equation (5), an equation (6) and an equation (7) respectively;
the water saving contribution rate
Is calculated as follows:
the alternative water source contribution rate
Is calculated as follows:
the water shortage contribution rate
Is calculated as follows:
wherein Q is1,i、Q2,iWater demand, G, of the i-th zone unit in the schemes 1 and 2, respectively1,i、G2,i、G3,iGroundwater feed, S, for the i-th zone unit in scheme 1, scheme 2, and scheme 3, respectively1,i,m、S2,i,mThe water supply amount of the mth type water source of the ith area unit in the scheme 1 and the scheme 2 is respectively; r1,i,m、R2,i,mThe water supply amount of the mth type water source of the ith area unit in the scheme 1 and the scheme 2 is respectively;
Q1,i、Q2,iwater demand predicted for the first module, G1,i、G2,i、G3,i、S1,i,m、S2,i,m、R1,i,m、R2,i,mBoth are GAMS optimized output.
Furthermore, in the first module, the water demand corresponding to the combination scheme is predicted by using a water demand quota method.
Furthermore, in the first module, under the situation scheme that general water saving is taken as water demand input, the water demand adopts the product of the current actual water usage quota and the planned socioeconomic index; under the scenario scheme of taking the intensified water conservation as the water demand input, the water demand adopts the product of the planned annual water demand quota and the planned socioeconomic index.
Furthermore, in the first module, the underground water pressure mining target is obtained from data compiled by a water administration department.
Compared with the prior art, the invention has the following characteristics:
(1) at present, no alternative water source contribution rate balance analysis system takes water shortage as a virtual water source for balancing the gap between water conservation, water increase and groundwater supply variation. The invention introduces a virtual water source, and can show the influence degree of multiple ways (water saving, multi-source water increasing and water shortage) on the groundwater supply variation.
(2) The method is simple and easy to implement, the complex multi-path alternative water source displacement relation is condensed into a paradigm, a decision maker can conveniently identify the action rule of each path on the underground water pressure mining, and decision support is provided for practice management work such as water distribution and the like.
Detailed Description
The following further describes the technical solutions of the present invention with reference to specific embodiments, so that those skilled in the art can better understand the present invention and can implement the present invention.
The system of the present invention will be described in detail below.
Scheme identification.
Based on the assumption that the effect generated when two major elements of the demand side and the supply side change simultaneously is equal to the effect accumulation generated when a single element changes, when a demand scheme (general water conservation W1 and intensified water conservation W2) and a groundwater supply constraint scheme (current mining L1 and pressing mining target L2), a pairwise combination scheme is formed: the relationships between the variation of groundwater supply and the water saving amount, the multi-source water supply amount and the water shortage amount under the two schemes are analyzed respectively when W1& L1, W2& L1, W2& L1 and W2& L2 are adopted. It is worth noting here that the W1& L2 solution does not actually exist, so the combination of W1& L2 and W1& L1 is not considered.
The scheme W1& L1 is a scenario scheme that the current annual underground water exploitation amount is used as water supply restriction, and general water saving is used as water demand input, and is called a general water saving & current excess mining scheme for short. The W2& L1 protocol refers to: the current situation scheme that the underground water exploitation amount is taken as water supply constraint and the intensified water saving is taken as water demand input is called 'intensified water saving & current situation excess mining scheme' for short. The W2& L2 protocol refers to: the scenario scheme that the underground water exploitation amount based on the underground water pressure exploitation target is taken as water supply constraint and the intensified water conservation is taken as water demand input is called as an intensified water conservation and pressure exploitation scheme for short.
For convenience of description, the W1& L1 scheme, the W2& L1 scheme, and the W2& L2 scheme will be hereinafter abbreviated as K1, K2, and K3 schemes, respectively.
And (II) predicting water requirement and determining an underground water pressure mining target.
The water demand forecasting method is a conventional method, and can generally adopt a water demand quota method for forecasting, the method is the product of the water demand quota and social economic indexes, the social economic indexes can adopt data in data such as ' statistical yearbook ', urban general planning ' and the like, and the water demand quota is determined after comprehensive analysis according to the local water demand quota standard, the current annual actual water demand quota and the planned annual water demand quota. For the invention, the water demand of the 'general water-saving scheme' is determined by the product of the current actual water demand quota and the planned socioeconomic index, and the water demand of the 'reinforced water-saving scheme' is determined by the product of the planned water demand quota and the planned socioeconomic index. The present invention uses the following table 1 to predict water demand for K1, K2, and K3 protocols, respectively. And for each scheme, respectively predicting the agricultural water demand, the industrial water demand, the domestic water demand and the ecological environment water supplement of each regional unit, and summing to obtain the water demand of the regional unit.
TABLE 1 Water demand rationing method forecast scheme water demand
The underground water pressure mining target can be obtained through data compiled by a water administrative department, and is generally obtained through files compiled by a local water administrative department, such as an underground water super mining comprehensive treatment action scheme, and if the files clearly suggest the underground water pressure mining target in the area, the underground water pressure mining target is directly adopted. And if the underground water exploitability is provided in the file, determining an underground water pressure exploitation target according to the difference value between the actual underground water exploitation amount and the underground water exploitability, wherein the actual underground water exploitation amount is from statistical data of water resource bulletin.
And (III) designing a balance analysis system for alternative water source contribution rate.
The alternative water source contribution rate balance analysis system is developed based on GAMS general modeling software, and realizes a balance analysis function by utilizing an optimization target and a linear solving algorithm. Various optimization algorithms are built in the GAMS, the method mainly adopts an LP algorithm (linear optimization algorithm), and the GAMS automatically optimizes output variables after input variables and constraint conditions are introduced into a GAMS interface. The method takes water demand Q and an underground water pressure mining target LIMIT as input variables, takes underground water supply not more than the underground water pressure mining target as a constraint condition (see formula 1), takes a formula (2) as an objective function, takes water supply S and water shortage N as output variables, and substitutes the output variables into a contribution rate calculation module in an overlapping mode.
The meaning of the objective function adopted by the optimization algorithm of the invention is as follows: the marginal benefits of different water sources supplying different zone units vary, for example: the benefits of external water supply for agriculture of a unit of a certain area and industrial use of a unit of another area are different, and besides considering the cost problem and the balanced water distribution problem, the optimization target of the objective function is to hope to obtain the maximum water supply benefit.
Gk,i≤LIMITk,i (1)
Formulae (1) to (5):
k denotes the combination scheme number, defining: k is 1 to represent general water saving and current super mining scheme, k is 2 to represent enhanced water saving and current super mining scheme, and k is 3 to represent enhanced water saving and pressure mining scheme;
i represents the unit numbers of different areas, and administrative units such as counties, districts and the like are generally adopted; il is the total number of area units;
m represents the number of the alternative water source types, and ml represents the total number of the alternative water source types; in the invention, the alternative water source mainly comprises 3 types of local surface water, external water and unconventional water, namely ml is 3;
t represents a time period number, and tl is the total time period number;
Gk,ithe groundwater supply quantity of the ith area unit in the kth scheme is represented and is an output variable after GAMS optimization;
LIMITk,irepresenting the underground water pressure mining target of the ith area unit in the kth scheme, wherein the underground water pressure mining target is the input quantity of an optimization algorithm;
δi,m,tindicating the marginal benefit of the water supply of the mth type water source of the ith area unit in the t time period; marginal benefit is calculated by adopting a marginal benefit function, specifically, a Cobb-Douglas production function can be adopted for solving, or the economic hypothesis which is the biggest in total profit is deduced based on the marginal cost equal to the marginal benefit;
Sk,i,mthe water supply amount of the mth type water source of the ith area unit in the kth scheme is represented and is the output amount after GAMS optimization;
Qk,irepresenting the water demand of the ith area unit in the kth scheme, which is the input quantity of the optimization algorithm;
Nk,ithe water shortage of the i-th area unit in the k scheme is represented and is the output quantity after GAMS optimization;
representing the alternative water source contribution, x corresponds to the contribution type,
respectively representing a water saving contribution rate calculation function, a multi-source water increasing contribution rate calculation function and a water shortage contribution rate calculation function.
The alternative water source contribution rate calculation module provided by the invention is nested in GAMS modeling software, takes the water supply amount S as an input variable and takes the contribution rate
Is a module that outputs variables.
The water source replacement and water saving are two important ways of underground hydraulic mining, and under the condition that no water shortage exists in the regional unit, the sum of the water saving contribution rate and the multi-source water increasing contribution rate is 1; under the condition that the regional unit has water shortage, considering the water shortage as a virtual water source, the water saving contribution rate
Multi-source water-increasing contribution rate
Water deficit contribution rate
The sum is 1. In the optimization algorithm, the water shortage N of the area unit under each combination scheme can be calculated by using a formula (4)
k,iAnd judging whether the area unit has water shortage under the combination scheme according to the calculated water shortage. When water shortage N
k,i0, water deficit; when water shortage N
k,iGreater than 0, indicating no water shortage; n is a radical of
k,iIf the GAMS is not more than 0, the implicit condition that the water supply does not exceed the water demand is considered during GAMS optimization. Contribution rate water-saving contribution rate and multi-source water increasing under water shortage situationThe contribution rate and the water shortage contribution rate have a relationship shown in formula (6):
in formula (6):
representing the water saving contribution rate of the ith zone unit;
(ii) an alternative water source contribution rate for the mth type of water source for the ith zone unit;
is the water shortage contribution rate of the unit of the ith area.
Water shortage or non-water shortage situation, water saving contribution rate
All calculated by equation (7):
in formula (7):
Q1,i、Q2,irespectively "general water conservation&Current super-mining scheme 'and' intensified water-saving&The water demand of the unit in the ith area in the current super-sampling scheme;
G1,i、G2,irespectively "general water conservation&Current super mining scheme 'and' mining intensified water saving&The groundwater supply amount of the i-th area unit in the current super scheme is the output amount after GAMS optimization;
S1,i,m、S2,i,mrespectively "general water conservation&Current super-mining scheme 'and' intensified water-saving&In the current super-sampling scheme, the water supply amount of the mth type water source of the ith area unit is the output amount after GAMS optimization.
Lack of water or notWater shortage situation, multi-source water increase contribution rate
All calculated by equation (8):
in formula (8):
G3,ito' intensify water conservation&The groundwater supply amount of the i-th area unit in the pressure mining scheme is the output amount after GAMS optimization;
R1,i,m、R2,i,mrespectively "intensified water saving&Current super-mining scheme 'and' intensified water-saving&And the water supply amount of the mth type water source of the ith area unit in the pressure collection scheme is the output amount after GAMS optimization.
Water deficit contribution rate
Is calculated as follows:
the fourth step: and (4) carrying out balance analysis on the contribution rate of the alternative water source.
The calculated contribution rate value can represent the influence of water saving and water increasing on the underground water mining of the regional unit. When in use
The larger the water is, the more obvious the effect of water saving on underground hydraulic mining of the regional unit is;
the larger the water source is, the more obvious the effect of increasing m types of water sources on underground hydraulic mining of the regional unit is;
the larger the area, the more water-saving the area unit,The water increasing way has no obvious effect on underground hydraulic pressure mining and has large difficulty in pressure mining. Therefore, the outputted contribution rate can be used as a reference for water resource management decision, and particularly, the decision reference can be provided as an alternative water source replacement path possibly involved in underground water pressure mining.
The technical solution provided by the present invention is not limited by the above embodiments, and all technical solutions formed by utilizing the structure and the mode of the present invention through conversion and substitution are within the protection scope of the present invention.