CN116402410A

CN116402410A - Distributed water quantity and water quality configuration method and system

Info

Publication number: CN116402410A
Application number: CN202310671504.5A
Authority: CN
Inventors: 杨明智; 许继军; 桑连海; 王冬; 景唤; 定培中; 宋权
Original assignee: Changjiang River Scientific Research Institute Changjiang Water Resources Commission
Current assignee: Changjiang River Scientific Research Institute Changjiang Water Resources Commission
Priority date: 2023-06-08
Filing date: 2023-06-08
Publication date: 2023-07-07
Anticipated expiration: 2043-06-08
Also published as: CN116402410B

Abstract

The invention discloses a distributed water quantity and water quality configuration method and system, and relates to the technical field of water quality and water quantity distribution, wherein the method comprises the following steps: dividing the configuration area into a sub-drainage basin, a sub-region and a computing unit by adopting a superposition type unit dividing method; determining daily water quantity water quality configuration data of the calculation unit, so as to determine the daily water quantity water quality configuration data, the monthly water quantity water quality configuration data and the annual water quantity water quality configuration data of the calculation unit and the subareas and configuration areas; the determination process of the daily water quantity and water quality configuration data of any computing unit comprises the following steps: determining water consumption parameters and land utilization types of water departments and water supply parameters of water sources; and obtaining daily water quantity and water quality configuration data of each calculation unit according to the water consumption parameters, the land utilization types and the water supply parameters. The invention realizes the fine management of regional water quantity and water quality.

Description

Distributed water quantity and water quality configuration method and system

Technical Field

The invention relates to the technical field of water quality and water quantity distribution, in particular to a distributed water quantity and water quality configuration method and system.

Background

Most of the traditional water resource allocation belongs to a centralized allocation mode, and administrative partitions, water resource partitions or nested partitions of the administrative partitions and the water resource partitions are generally adopted as water resource allocation partitions. The spatial differences of natural geographic conditions (such as precipitation, terrain, land utilization and the like) and economic and social parameters (such as population density, industrial structure and total quantity, crop types and planting areas thereof, irrigation areas, point source positions and emission intensity thereof and the like) in the subareas are ignored in space, and the configuration result based on the spatial differences is only the average value of the subareas, so that the problem of water resources at the corner positions is easily ignored. The lumped configuration mode generally adopts time steps of month, ten days and the like in time, and ignores the feedback effect of dynamic changes of water quantity and water quality on water resource configuration. These are detrimental to the fine management of water resources.

Along with the development of the distributed hydrologic model, a dynamic water resource allocation method is provided, namely, a lumped water resource allocation model is coupled with a traditional distributed hydrologic model, water resource boundary conditions are provided by the distributed hydrologic model, the lumped allocation model provides calculation results of water taking, water supply, water consumption, water drainage and pollution discharge, and the model coupling is realized by means of output file transmission or parameter unidirectional transmission, so that the dynamic water resource allocation is completed. For example, zhao Yong and the like perform information interaction in a decomposition and aggregation mode, and perform unidirectional coupling on a plain area distributed water circulation model (PDMCM) and a generalized water resource reasonable allocation model (GWRAM), so that drainage basin generalized water resource allocation is realized. Gu Angwen and the like realize loose coupling of a distributed hydrological model and a lumped water resource allocation model through information transfer between a WEP model and an RPWAS model, and develop a comprehensive water quality simulation platform with macroscopic economy, water resources, ecology and environment functions for the water quantity and the water quality of a river basin. Zhang Shouping with the water resource partition divided by the three-level water resource region nested city range as a basic configuration unit, transmitting data to a sub-basin elevation zone of a hydrologic model, realizing the coupling of the water resource configuration model and the WEP-L model, adding a QUAL2E module, and constructing a water quantity and water quality joint configuration model based on water circulation and associated processes. What has been proposed is a distributed water resource allocation platform DTVGM-WEAR based on a binary water circulation process based on a distributed hydrological model DTVGM. However, the method belongs to a loose coupling mode, ignores the dynamic interaction between natural water circulation and social water circulation, cannot well reflect the interaction between the water using process and the natural hydrologic process of human activities, easily causes deviation of water resource allocation results, and is not beneficial to fine management of water resources.

In recent years, some experts propose a tight coupling method of a water resource allocation model and a distributed hydrological model, a water resource allocation module is embedded in the structural design of the hydrological model or a hydrological module is embedded in the design of the water resource allocation model, two-way coupling of the two is completed, and accurate simulation and dynamic allocation of a water resource complex system are finally realized. For example, sang Xuefeng and the like construct a conceptual distributed water resource comprehensive simulation and allocation model WAS by expanding the hydrologic cycle simulation function of the water resource allocation model, and realize the real-time linkage mutual feedback simulation of the natural-social water cycle process; yang Mingzhi and the like expand the water resource allocation simulation function by improving the SWAT model, develop a distributed water resource allocation model SWAT-WARM based on water circulation, and realize the integrated simulation of water circulation and water resource allocation in a river basin. However, the WAS has simplified consideration of hydrologic cycle and accompanying pollution process, and SWAT-WARM does not consider the influence of point source and surface source pollution on water resource allocation, so that the requirement for fine management of regional water quantity and water quality is still insufficient.

Disclosure of Invention

The invention aims to provide a distributed water quantity and water quality configuration method and system, which realize the fine management of regional water quantity and water quality.

In order to achieve the above object, the present invention provides the following solutions:

a distributed water volume water quality configuration method, the method comprising:

dividing the whole configuration area into a plurality of sub-domains;

dividing the sub-drainage basin into a plurality of hydrologic simulation units according to the underlying surface information of the sub-drainage basin; the underlying information includes land use type and soil type; the land use types include construction land and farmland land;

dividing the hydrologic simulation unit into a plurality of calculation units according to the range of the subareas in the configuration area and the range of the hydrologic simulation unit;

determining daily water quantity and water quality configuration data of each calculation unit, and determining the daily water quantity and water quality configuration data of the subareas, the daily water quantity and water quality configuration data of the sub-watershed and the daily water quantity and water quality configuration data of the configuration area based on the daily water quantity and water quality configuration data of all calculation units by adopting a superposition method; the daily water quantity and water quality configuration data comprises: daily water intake of the water source, daily water demand, daily water consumption, daily water discharge and daily pollution discharge of each water department in the calculation unit;

determining the lunar water quality configuration data of the computing units based on the solar water quality configuration data of the computing units by adopting a superposition method, and determining the lunar water quality configuration data of the subareas, the lunar water quality configuration data of the sub-basins and the lunar water quality configuration data of the configuration areas based on the lunar water quality configuration data of all the computing units by adopting the superposition method; the lunar water quantity and water quality configuration data comprises: the month water intake of the water source, the month water demand, month water consumption, month water discharge and month pollution discharge of each water department in the calculation unit;

Determining annual water quantity water quality configuration data of the computing units based on the monthly water quantity water quality configuration data of the computing units by adopting a superposition method, and determining annual water quantity water quality configuration data of the subareas, annual water quantity water quality configuration data of the sub-basins and annual water quantity water quality configuration data of the configuration areas based on the annual water quantity water quality configuration data of all the computing units by adopting the superposition method; the annual water quantity and water quality configuration data comprises: annual water intake of the water source, annual water demand, annual water consumption, annual water discharge and annual pollution discharge of each water department in the calculation unit;

wherein, the determination process of the daily water quantity and water quality configuration data of any computing unit comprises the following steps:

determining water consumption parameters and land utilization types of all water consumption departments of the computing unit; the water consumption parameters include: the water demand, the water sequence and the upper limit of the total water amount of the allocation area, wherein the land utilization type comprises construction land and farmland land;

determining water supply parameters of each water source; the water supply parameters include: water supply sequence, water quality, water supply availability, water supply availability, surface water resource utilization control amount, and groundwater resource utilization control amount;

And determining the daily water quantity of the corresponding water department according to the water consumption parameter, the land utilization type and the water supply parameter of each water supply source, thereby obtaining daily water quantity and water quality configuration data of each calculation unit.

According to the specific embodiment provided by the invention, the invention discloses the following technical effects:

the invention discloses a distributed water quantity and water quality configuration method and system, wherein the method comprises the following steps: dividing the configuration area into a sub-drainage basin, a sub-region, a hydrological simulation unit and a calculation unit according to the hierarchy; determining daily water quantity and water quality configuration data of each computing unit, thereby determining daily water quantity and water quality configuration data of the subareas, the sub-waterbasins and the configuration areas; determining the monthly water quantity water quality configuration data of the calculation unit based on the daily water quantity water quality configuration data of the calculation unit by adopting a superposition method, thereby determining the monthly water quantity water quality configuration data of the subareas, the sub-drainage basins and the configuration areas; determining annual water quantity water quality configuration data of the computing unit based on the monthly water quantity water quality configuration data of the computing unit by adopting a superposition method, thereby determining annual water quantity water quality configuration data of the subareas, sub-waterbasins and configuration areas; wherein, the determination process of the daily water quantity and water quality configuration data of any computing unit comprises the following steps: determining water consumption parameters and land utilization types of all water consumption departments of the computing unit; determining water supply parameters of each water source; and determining the daily water quantity of the corresponding water department according to the water consumption parameter, the land utilization type and the water supply parameter of each water supply source, thereby obtaining daily water quantity and water quality configuration data of each calculation unit. The water quantity and water quality configuration method considers the configuration of different time-space scales such as year-month-day and river basin-administrative area-irrigation area-water function area-calculation unit and the like, and realizes the fine management of regional water quantity and water quality.

Drawings

In order to more clearly illustrate the embodiments of the present invention or the technical solutions of the prior art, the drawings that are needed in the embodiments will be briefly described below, it being obvious that the drawings in the following description are only some embodiments of the present invention, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art.

FIG. 1 is a schematic flow chart of a distributed water volume and water quality configuration method provided by an embodiment of the invention;

fig. 2 is a schematic structural diagram of a distributed water volume and water quality configuration system according to an embodiment of the present invention.

Detailed Description

The following description of the embodiments of the present invention will be made clearly and completely with reference to the accompanying drawings, in which it is apparent that the embodiments described are only some embodiments of the present invention, but not all embodiments. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.

The invention aims to provide a distributed water quantity and water quality configuration method and system, and aims to realize the fine management of regional water quantity and water quality.

In order that the above-recited objects, features and advantages of the present invention will become more readily apparent, a more particular description of the invention will be rendered by reference to the appended drawings and appended detailed description.

Fig. 1 is a schematic flow chart of a distributed water quantity and water quality configuration method according to an embodiment of the invention. As shown in fig. 1, the distributed water volume and water quality configuration method in this embodiment includes:

step 101: the entire configuration area is divided into a plurality of sub-domains.

Step 102: dividing the sub-watershed into a plurality of hydrologic simulation units according to the underlying surface information of the sub-watershed; the underlying information includes land use type and soil type; land use types include construction land and farmland land.

Specifically, land use types also include other land types in addition to construction land and farmland land.

Step 103: dividing the hydrologic simulation unit into a plurality of calculation units according to the range of the subareas in the configuration area and the range of the hydrologic simulation unit; the subareas comprise county-level administrative areas, water function areas and irrigation areas; the calculation units are divided into town calculation units, rural calculation units, farmland calculation units, woodland calculation units, grassland calculation units and water area calculation units.

Step 104: the daily water quantity and water quality configuration data of each calculation unit are determined, and a superposition method is adopted to determine the daily water quantity and water quality configuration data of the subareas, the daily water quantity and water quality configuration data of the sub-basins and the daily water quantity and water quality configuration data of the configuration areas based on the daily water quantity and water quality configuration data of all calculation units.

Wherein, daily water yield quality of water configuration data includes: daily water intake of the water source, daily water demand, daily water consumption, daily water discharge and daily pollution discharge of each water department in the calculation unit; the water department includes: urban resident water sector, rural resident water sector, industrial water sector, building industry water sector, service industry water sector, livestock and poultry breeding water sector, urban ecological water sector, rural ecological water sector and irrigation water sector.

Step 105: the method comprises the steps of determining the lunar water quantity water quality configuration data of a calculation unit based on the solar water quantity water quality configuration data of the calculation unit by adopting a superposition method, and determining the lunar water quantity water quality configuration data of a subarea, the lunar water quantity water quality configuration data of a sub-basin and the lunar water quantity water quality configuration data of a configuration area based on the lunar water quantity water quality configuration data of all the calculation units by adopting the superposition method.

Wherein, the lunar water quantity and water quality configuration data comprises: the month water intake of the water source and the month water demand, month water consumption, month water discharge and month pollution discharge of each water using department in the calculating unit.

Step 106: the method comprises the steps of determining annual water quantity water quality configuration data of a calculation unit based on monthly water quantity water quality configuration data of the calculation unit by adopting a superposition method, and determining annual water quantity water quality configuration data of a subarea, annual water quantity water quality configuration data of a sub-basin and annual water quantity water quality configuration data of a configuration area based on annual water quantity water quality configuration data of all the calculation units by adopting the superposition method.

Wherein, annual water volume water quality configuration data includes: annual water intake of the water source, annual water demand, annual water consumption, annual water discharge and annual pollution discharge of each water department in the calculation unit.

determining water consumption parameters and land utilization types of all water consumption departments of the computing unit; the water consumption parameters include: the water demand, the water use sequence and the upper limit of the total water use amount of the allocation area, and land use types include construction land and farmland land.

Determining water supply parameters of each water source; the water source includes: river course, reservoir, shallow aquifer, deep aquifer, pit pool, external water source, regeneration water factory, rain collecting pool and desalination water factory, water supply parameters include: water supply sequence, water quality, water supply availability, water supply availability, surface water resource utilization control and groundwater resource utilization control.

As an alternative embodiment, determining the daily water amount of the corresponding water department according to the water consumption parameter, the land utilization type and the water supply parameter of each water supply source, thereby obtaining daily water amount and water quality configuration data of each calculation unit, specifically comprising:

and determining a configuration mode according to the land utilization type of the computing unit.

When the land use type is a construction land, the calculation units comprise urban resident water departments, rural resident water departments, industrial water departments, building industry water departments, service industry water departments, livestock and poultry raising water departments, urban ecological water departments and rural ecological water departments, and the water quantity and the water quality of each calculation unit are configured according to water consumption parameters and water supply parameters of each water supply source by adopting a first configuration mode.

When the land use type is farmland land, the calculation units comprise irrigation water departments, and the water quantity and quality configuration is carried out on each calculation unit according to the water consumption parameters and the water supply parameters of each water supply source by adopting a second configuration mode.

As an alternative embodiment, a first configuration mode is adopted, and the water quantity and quality configuration is carried out on each computing unit according to the water consumption parameter and the water supply parameter of each water supply source, which specifically comprises the following steps:

and carrying out water quantity and quality configuration on each calculation unit in turn according to the priority of each calculation unit.

When the current calculation unit is configured with water quantity and water quality, each water department is configured with water quantity and water quality in sequence according to the water use sequence.

When the current water utilization department carries out water quantity and quality configuration, each water source is sequentially called according to the water supply sequence.

When the current water source is called, judging whether the water quality of the current water source meets the preset water quality condition or not, and obtaining a first judging result.

If the first judgment result is negative, replacing the current water source with the next water source, and returning to judge whether the water quality of the current water source meets the preset water quality condition to obtain the first judgment result.

If the first judgment result is that the current water source belongs to surface water, judging whether the current surface water resource utilization control amount is larger than the available water amount of the current water source or if the first judgment result is that the current water source belongs to underground water, judging whether the current underground water resource utilization control amount is larger than the available water amount of the current water source, and obtaining a second judgment result.

If the second judgment result is that the current water source belongs to surface water, judging whether the underground water resource utilization control amount is larger than the available water amount of the current water source or if the second judgment result is that the current water source belongs to underground water, judging whether the surface water resource utilization control amount is larger than the available water amount of the current water source, and obtaining a third judgment result.

If the third judgment result is that the current water source belongs to surface water, replacing the current water source with the next water source belonging to the ground water type or if the third judgment result is that the current water source belongs to ground water, replacing the current water source with the next water source belonging to the surface water type.

If the third judgment result is negative, the current water use department is replaced by the next water use department, and the water source is returned to the current water use department, and when the water quantity and the water quality of the current water use department are configured, all water sources are sequentially called according to the water supply sequence.

If the second judgment result is yes, updating the current water resource utilization control amount of the current water resource according to the current water resource utilization control amount and the water supply amount of the current water source, judging whether the upper limit of the current water total amount of the configuration area is larger than the water supply amount of the current water source, and obtaining a fourth judgment result.

If the fourth judgment result is negative, the current water use department is replaced by the next water use department, and the water source is returned to the current water use department, and when the water quantity and the water quality of the current water use department are configured, all water sources are sequentially called according to the water supply sequence.

If the fourth judgment result is yes, updating the current water consumption total amount upper limit of the configuration area according to the current water consumption total amount upper limit of the configuration area and the current water source water supply amount, judging whether the available water amount of the current water source is larger than the water demand of the current water use department, and obtaining a fifth judgment result.

If the fifth judgment result is yes, replacing the water supply amount of the current water source with the water demand of the current water department, replacing the water demand of the current water department with zero, replacing the current water department with the next water department, and returning to 'when the current water department carries out water quantity and water quality configuration, sequentially calling the water sources according to the water supply sequence'.

If the fifth judgment result is negative, replacing the water supply amount of the current water source with the available water amount of the current water source, updating the water demand of the current water department, judging whether the next water source exists, and obtaining a sixth judgment result.

If the sixth judgment result is yes, replacing the current water source with the current water source and replacing the current water source with the next water source, and returning to judge whether the water quality of the current water source meets the preset water quality condition or not to obtain a first judgment result.

If the sixth judgment result is negative, judging whether the next water use department exists, and obtaining a seventh judgment result.

If the seventh judgment result is yes, the current water department is replaced by the next water department, and the water source is returned to the step of calling each water source in sequence according to the water supply sequence when the current water department carries out water quantity and water quality configuration.

If the seventh judgment result is no, the daily water quantity and water quality configuration data of the calculation unit are obtained when the land utilization type is the construction land, the current calculation unit is replaced by the next calculation unit, and the water quantity and water quality configuration is returned to each water department in turn according to the water use sequence when the current calculation unit is configured, until the daily water quantity and water quality configuration data of all calculation units are obtained.

As an optional implementation manner, a second configuration manner is adopted, and the water quantity and quality configuration is carried out on each calculation unit according to the water consumption parameter and the water supply parameter of each water supply source, which specifically comprises the following steps:

When the current computing unit is configured with water quantity and water quality, each water source is sequentially called according to the water supply sequence to configure the water quantity and water quality of the irrigation water department.

When the current water source is called, judging whether the water quality of the current water source meets the preset water quality condition or not, and obtaining an eighth judging result.

If the eighth judgment result is no, replacing the current water source with the next water source, and returning to judge whether the water quality of the current water source meets the preset water quality condition to obtain the eighth judgment result.

If the eighth judgment result is that the current water source belongs to surface water, judging whether the current surface water resource utilization control amount is larger than the available water amount of the current water source or if the eighth judgment result is that the current water source belongs to underground water, judging whether the current underground water resource utilization control amount is larger than the available water amount of the current water source, and obtaining a ninth judgment result.

If the ninth judgment result is that the current water source belongs to surface water, judging whether the underground water resource utilization control amount is larger than the available water amount of the current water source or if the ninth judgment result is that the current water source belongs to underground water, judging whether the surface water resource utilization control amount is larger than the available water amount of the current water source, and obtaining a tenth judgment result.

If the tenth judgment result is that the current water source belongs to surface water, replacing the current water source with the next water source belonging to the ground water type or if the tenth judgment result is that the current water source belongs to ground water, replacing the current water source with the next water source belonging to the surface water type.

If the tenth judgment result is negative, the current calculation unit is replaced by the next calculation unit, and the method returns to the step of sequentially calling each water source to carry out water quantity and water quality configuration on the irrigation water department according to the water supply sequence when the current calculation unit carries out water quantity and water quality configuration.

If the ninth judgment result is yes, updating the current water resource utilization control amount of the current water resource according to the current water resource utilization control amount and the water supply amount of the current water source, and judging whether the upper limit of the current water total amount of the configuration area is larger than the water supply amount of the current water source or not, so as to obtain an eleventh judgment result.

If the eleventh judgment result is no, replacing the current calculation unit with the next calculation unit, and returning to the step of sequentially calling each water source to perform water quantity and water quality configuration on the irrigation water department according to the water supply sequence when the current calculation unit performs water quantity and water quality configuration.

If the eleventh judgment result is yes, updating the current water consumption total amount upper limit of the configuration area according to the current water consumption total amount upper limit of the configuration area and the current water source water supply amount, judging whether the available water amount of the current water source is larger than the water demand of the current water use department, and obtaining a twelfth judgment result.

If the twelfth judgment result is yes, replacing the water supply amount of the current water source with the water demand of the irrigation water department, replacing the water demand of the irrigation water department with zero, replacing the current calculation unit with the next calculation unit, and returning to the step of sequentially calling each water source to carry out water quantity and water quality configuration on the irrigation water department according to the water supply sequence when the current calculation unit carries out water quantity and water quality configuration.

If the twelfth judgment result is negative, replacing the water supply amount of the current water source with the available water amount of the current water source, updating the water demand of the irrigation water department, judging whether the next water source exists, and obtaining the thirteenth judgment result.

If the thirteenth judgment result is yes, replacing the current water source with the current water source and replacing the current water source with the next water source, and returning to judge whether the water quality of the current water source meets the preset water quality condition or not to obtain an eighth judgment result.

If the thirteenth judgment result is negative, the daily water quantity and water quality configuration data of the current calculation unit are obtained, the current calculation unit is replaced by the next calculation unit, and the water quantity and water quality configuration is returned to each water department in turn according to the water sequence when the current calculation unit is configured, until the daily water quantity and water quality configuration data of all calculation units are obtained.

Specifically, in order to implement the above method, the process of performing model construction is as follows:

step1: the computing unit divides the design.

The calculation unit division design refers to further dividing a hydrologic simulation unit of the hydrologic model on the basis of slope surface dispersion of the hydrologic model by comprehensively considering distribution conditions of irrigation areas, administrative areas and water function areas to form the calculation unit.

(1) Hydrologic simulation unit division: the distributed hydrologic model divides a research area (namely a configuration area) into a plurality of sub-watersheds according to the topographic information, and then divides the sub-watersheds into a plurality of hydrologic simulation units according to the underlying surface information in each sub-watershed.

(2) Dividing a computing unit: according to the distribution of the county administrative district range and the irrigation district range in the research district, the hydrologic simulation unit is overlapped with the water function district, the county administrative district and the irrigation district, and the hydrologic simulation unit is divided into a plurality of calculation units. It is possible to specify which basin, which administrative district, which water function district, which irrigation district, which land use type and soil type a certain computing unit belongs to. The calculation unit can be divided into various types such as town calculation unit, rural calculation unit, farmland calculation unit, forest calculation unit, grassland calculation unit, water area calculation unit, etc. according to land utilization attribute.

Step2: and (5) national economy module design.

The national economy module design refers to the design of an exponential national economy development prediction program, reads economic and social development data and change rate data thereof, predicts population and national economy total amount in a designated year, and spreads the population and national economy total amount to a calculation unit of the fine configuration module.

And the national economy module obtains national economy data such as population number, livestock and poultry breeding number, GDP total amount and the like which take county administrative regions as units, and population growth rate, urbanization rate, livestock and poultry breeding number growth rate, GDP growth rate and three-time industrial structure data by reading the input social economy data file. Future population and economic and social data are predicted based on the specified predicted year.

The urban and rural population quantity prediction calculation formula is as follows:

P _t =P ₀ ·(1+g _p ) ^yr 。

P _urban =P _t ·r _urban 。

P _rural =P _t -P _urban 。

wherein P is _t To predict the population total for year t, the unit is ten thousand; p (P) ₀ The population total amount of the reference annual county level administrative district is ten thousand; g _p Annual population growth rate for county-level administrative areas; yr is the predicted years; p (P) _urban 、P _rural The urban population number and the rural population number of the county-level administrative areas are respectively ten thousands of people; r is (r) _urban Is the town ratio of the county-level administrative district.

The prediction formula of the livestock and poultry breeding quantity is as follows:

AN _t =AN ₀ ·(1+g _an ) ^yr 。

Wherein AN _t The total amount of livestock and poultry in the predicted year t of the county administrative district is ten thousand units; AN (Access network) ₀ The total amount of livestock and poultry in standard year is ten thousand; g _an The annual average growth rate of livestock and poultry in county administrative areas is increased.

The GDP total and three industrial increment value prediction formulas are as follows:

GDP _t =GDP ₀ ·(1+g _g ) ^yr 。

GDP _a =GDP _t ·e _a (a=1,2,3)。

。

GDP _ind =GDP ₂ ·e _ind 。

GDP _con =GDP ₂ -GDP _ind 。

wherein, GDP _t Predicting the total amount of GDP (in terms of GDP invariance) of year t for a county-level administrative area in hundred million yuan; GDP (GDP) ₀ The total GDP amount is the unit of hundred million yuan for the reference annual county level administrative district; g _g The GDP of the county administrative district is increased in annual average speed; GDP (GDP) _a (a=1, 2, 3) is the increment value of the a industry in the forecast year t in hundred million yuan in the county administrative district; e, e _a (a=1, 2, 3) is the percentage of the incremental value of industry a in the predicted year t to the total GDP in the county administrative district; GDP (GDP) _ind 、GDP _con Each of county-level administrative areasIndustrial increment value and building industry increment value in hundred million yuan in the predicted year t; e, e _ind The industrial increment value of the predicted year t for the county administrative area accounts for the percentage of the second industrial increment value; GDP (GDP) ₂ The increment value of the year t is predicted for the second industry of the county-level administrative district, and the unit is hundred million yuan; a=2 represents a first industry, a=2 represents a second industry, and a=3 represents a third industry.

And the national economy module performs space distribution on the predicted economic and social data of the county administrative district according to the urban and rural land utilization type distribution and the area size in the administrative district. The town population, the industry added value, the building industry added value and the service industry added value are distributed on the town calculating unit, and the rural population and the livestock and poultry breeding quantity are distributed on the rural calculating unit.

For town calculating units, the spread formula is:

。

wherein, data _urban Urban national economy total data (including urban population total, industrial added value total, building added value total and service added value total) of county-level administrative areas; data _urban,p National economy data of a p-th town calculating unit in a county-level administrative area; area (Area) _urban,p Calculating the area of a unit for the p-th town in the county-level administrative area; t is the total number of town calculation units in the county level administrative district.

For rural computing units, the spread formula is:

。

wherein, data _rural Rural national economy total data (including rural population total and livestock and poultry raising total) of county-level administrative areas; data _rural,q National economy data of a q-th rural calculation unit in a county-level administrative area; area (Area) _rural,q Calculating the area of a unit for the q-th rural area in the county-level administrative area; f is a rural computing unit in a county-level administrative districtTotal number.

Step3: and (5) refining the design of the configuration module.

The fine configuration module design refers to the fine configuration program design of water quantity and water quality, and comprises a water demand module, a water quantity and water quality configuration module, a water consumption and drainage module and a pollution discharge module.

(1) Water requiring module

The water demand module is connected with the national economy module by receiving various national economy data on the computing unit, obtains the water consumption intensity and water consumption efficiency data of urban and rural resident domestic water quota, the water consumption in ten thousand yuan industry increment value, the water delivery pipe network leakage loss rate, the irrigation water effective utilization coefficient and the like by reading the input water consumption intensity data file, and simulates and calculates the daily water demand of each water department. The calculation formula is as follows:

WD _urban =P _urban ·d _urban /1000/(1-a _pip )。

WD _rural =P _rural ·d _rural /1000。

Wherein WD _urban The daily water demand of urban residents in the urban calculation unit is expressed as m ³ /d；d _urban The unit is L/(d. Person) for the urban resident domestic water intensity in the urban calculation unit; a, a _pip The leakage rate of the water delivery pipe network; WD _rural The daily water demand of rural residents as rural calculation units is expressed as m ³ /d；d _rural The unit is L/(d. Person) for the domestic water intensity of village residents in the rural calculation unit.

The daily water demand for livestock and poultry cultivation is calculated by the following formula: WD _an =P _an ·d _an /1000。

Wherein WD _an The daily water demand of livestock and poultry in rural calculation units is m ³ /d；P _an The total amount of the livestock and the poultry is ten thousand; d, d _an The water intensity of livestock and poultry in rural calculation units is L/(d.head).

The water demand of industry, construction industry and service industry calculates daily water demand according to the industry increment value and the water quota for ten thousand yuan industry increment value,the calculation formulas of water demand in three industries of industry, construction industry and service industry can be as follows: WD _e =d _e ·GDP _e /(1-a _pip )/num。

Wherein WD _e Daily water demand in m for the industrial or construction or service industry of town calculating units ³ /d；d _e The unit of the industry or the construction industry or the service industry is ten thousand yuan industry to increase the water quota for the value added, m ³ The unit is ten thousand yuan; GDP (GDP) _e The value of daily industry is increased for the industry or the building industry or the service industry, and the unit is ten thousand yuan/d; num is the total number of days of the year.

Agricultural irrigation water is calculated according to the area of an irrigation calculation unit, irrigation time and irrigation quantity input into the distributed hydrologic model by adopting the following formula: WD _irr =Area _irr ·I _irr ·10/b _irr 。

Wherein WD _irr Agricultural irrigation daily water demand as farmland calculating unit with unit of m ³ /d；I _irr The unit is mm for the single irrigation water quantity; area (Area) _irr The unit is ha which is the area of the irrigation calculation unit; b _irr The coefficient is effectively utilized for irrigation water.

(2) Water quantity and water quality configuration module

The water quantity and water quality configuration module is connected with the water demand module by receiving the water demand data of each department on the calculation unit, outputs the water consumption data of the calculation unit and transmits the water consumption data to the water consumption and drainage module. The water quantity and water quality configuration module comprises a construction land configuration sub-module and an irrigation land configuration sub-module. The program calls a construction land configuration sub-module by identifying the land utilization type of the computing unit and if the land is a construction land; if the irrigation land is the farmland land, the program calls an irrigation land configuration program; if the land use is of other land use types, the program ends and the next calculation unit cycle is entered.

1) Construction land allocation submodule

The program calculation flow in the construction land configuration submodule comprises the steps of department code and water use sequence information identification, water source code and water supply sequence information identification, water source module calling, water source quality condition identification, water resource total amount utilization limit identification, water use total amount limit identification, water source water taking calculation and the like. The module outputs daily water consumption data of all water departments and transmits the daily water consumption data to the water consumption and drainage module. The operation flow of the construction land configuration submodule is as follows:

(1) Department code and water sequence information identification

First, the set department codes and water use priorities of each construction land are identified, and the water source water supply calculation cycles are sequentially carried out. For town computing units, the water sector includes town resident life, town ecology, construction, industry, and service. For rural computing units, the water use departments comprise rural resident life, rural ecology and livestock and poultry cultivation. The department water sequence may be set according to the actual demand of each computing unit.

(2) Water source code and water supply sequence information identification

The program identifies the number of water sources, the water source codes and the water supply priority order set by each water department, and calls each water source module in sequence. The water supply source type comprises 9 water sources such as river channels, reservoirs, shallow aquifers, deep aquifers, pits, external water sources, regenerated water plants, rain collecting ponds, desalinated water plants and the like. The number of each type of water source is not limited, and is denoted by I, J, K, L, M, N, X, Y, Z.

(3) Water source module call

The program reads the water source identification codes of the corresponding water sources sequentially according to the water supply priority sequence from 1 to o (o=i+j+k+l+m+n+x+y+z) (the water source water supply priority sequence is set according to the need and is an input item of a model (see (1)). For example, the water supply source of the computing unit is river channel 1, river channel 5, reservoir 2 and reservoir 6. The water supply priority sequence of the computing unit can be set such that reservoir 2> river channel 1> reservoir 6> river channel 5). And when the water source identification code is 1, the river channel module is called, when the water source identification code is 2, the reservoir module is called, when the water source identification code is 3 and 4, the underground water module is called, and when the water source identification code is 5, 6, 7, 8 and 9, the external water regulation module, the pit pool module, the reclaimed water module, the rainwater collecting module and the desalted water module are respectively called.

(4) Water source quality condition identification

After the water source module is called, the program calculates the water quality concentration of the water source according to the pollution load quantity and the water resource quantity result simulated by the distributed hydrologic model, and identifies whether the water quality of the water supply source meets the water taking requirement according to the input water quality requirement file of the water supply source. The water supply source water quality requirement file contains a water quality concentration limit value for each water source, and if the water quality concentration of a certain water source exceeds the water quality concentration limit value on a certain day, water cannot be taken from the water source. If the water resource utilization total amount limit identification process meets the requirement, the program enters the water resource utilization total amount limit identification process; otherwise, the program searches for a water source with the secondary water supply priority; if the water supply source is not available, the program judges whether the computing unit has the next water use department, and if so, the program enters the next water use department for circulation; if not, the program enters the next calculation unit loop.

(5) Water resource utilization total amount limiting identification

If the water quality condition of the water supply source is proper, the program enters the water resource utilization total amount limit identification. Program sets the control quantity Lim of surface water resource utilization _s And groundwater resource utilization control amount Lim _g Two constraints (for a particular program, lim _s And Lim _g Is a fixed value). For a water source with water supply priority J (j=1, 2, …, o; o=i+j+k+l+m+n+x+y+z), when the water intake source belongs to a surface water source, if Lim is present _s ≥W _j （W _j =min(WD,Wsc _j ) Wherein W is _j Water supply capacity for water source with priority j; wsc _j The total water quantity of the water source with the priority order of j; WD is the water demand of the sector. ) The control amount of the surface water resource utilization is reduced to Lim _s -W _j And the program enters a water total amount limiting and identifying link; if there is Lim _s <W _j And Lim _g ≥W _j The program searches for a secondary water supply priority water source until a water source belonging to the groundwater type is found, and then enters a water total amount limiting and identifying link; if there is Lim _s <W _j And Lim _g ＜W _j The program goes to the next water department cycle. When the water intake source belongs to the underground water source, if Lim exists _g ≥W _j The groundwater resource utilization control amount is reduced to Lim _g -W _j And the program enters a water total amount limiting and identifying link; if there is Lim _g <W _j And Lim _s ≥W _j The program searches for a secondary water supply priority water source until a water source belonging to the surface water type is found, and then enters a water total amount limiting and identifying link; if there is Lim _g <W _j And Lim _s ＜W _j The program goes to the next water department cycle.

(6) Water total limit identification

The program sets a limit for total water consumption in the area ₃ (for a particular program, lim) ₃ Fixed value). For water sources with water supply priority J (j=1, 2, …, o; o=i+j+k+l+m+n+x+y+z), if Lim ₃ ≥W _j The area water control amount is reduced to Lim ₃ -W _j And the program enters a water source water taking calculation link; if Lim ₃ <W _j The program goes to the next water department cycle.

(7) Water source water intake calculation

For a water source with water supply priority J (j=1, 2, …, o; o=i+j+k+l+m+n+x+y+z), if the available water amount Wsc of the water source is available _j > WD, then the water supply to the water source is WSP _j WD, and the water demand of this water sector is reduced to 0, the program goes to the next water sector cycle; otherwise, the water source supplies water to WSP _j =Wsc _j And the water demand of the water department is reduced to WD-Wsc _j The program will find the secondary water supply priority source and continue the water supply calculation until the water demand of the water department is met or until the last water supply source is finished. Calculating the water consumption of the water department

. The program goes to the next water department cycle.

2) Irrigation land configuration submodule

Because the irrigation land only has one water using department of the irrigation land, the water using departments and the water using sequence thereof do not need to be identified, and therefore, the program design is only less than the construction land configuration submodule by one department cycle. The program calculation flow in the irrigation land configuration submodule comprises the steps of water source code and water supply sequence information identification, water source module calling, water source quality condition identification, water resource total utilization amount limit identification, water use amount limit identification, water source water taking calculation and the like. The module outputs irrigation water consumption data and transmits the irrigation water consumption data to the water consumption and drainage module. The operation flow of the irrigation land configuration submodule is as follows:

(1) Water source code and water supply sequence information identification

The program identifies the number of water sources, the water source codes and the water supply priority order set by each irrigation calculation unit, and calls each water source module in sequence. The water supply source type comprises 9 water sources such as river channels, reservoirs, shallow aquifers, deep aquifers, pits, external water sources, regenerated water plants, rain collecting ponds, desalinated water plants and the like. The number of each type of water source is not limited, and is denoted by I, J, K, L, M, N, X, Y, Z.

(2) Water source module call

The program reads the water source identification codes of the corresponding water sources in order of water supply priority from 1 to o (o=i+j+k+l+m+n+x+y+z). And when the water source identification code is 1, the river channel module is called, when the water source identification code is 2, the reservoir module is called, when the water source identification code is 3 and 4, the underground water module is called, and when the water source identification code is 5, 6, 7, 8 and 9, the external water regulation module, the pit pool module, the reclaimed water module, the rainwater collecting module and the desalted water module are respectively called.

(3) Water source quality condition identification

After the water source module is called, the program calculates the water quality concentration of the water source according to the simulated pollution load and the water resource quantity result, and identifies whether the water quality of the water supply source meets the irrigation water taking requirement according to the input water quality requirement file. If the water resource utilization total amount limit identification process meets the requirement, the program enters the water resource utilization total amount limit identification process; otherwise, the program searches for a water source with the secondary water supply priority; if there is no water supply, the program proceeds to the next calculation unit cycle.

(4) Water resource utilization total amount limiting identification

If the water quality condition of the water supply source is proper, the program enters the water resource utilization total amount limit identification. Program sets the control quantity Lim of surface water resource utilization _s And groundwater resource utilization control quantity Limg. For a water source with water supply priority J (j=1, 2, …, o; o=i+j+k+l+m+n+x+y+z), when the water intake source belongs to a surface water source, if Lim is present _s ≥W _j （W _j =min (WD, wscj)), the surface water resource utilization control amount is reduced to Lim _s Wj and the program goes to the water total limit identification link; if there is Lim _s <W _j And Lim _g ≥W _j The program searches for a secondary water supply priority water source until a water source belonging to the groundwater type is found, and then enters a water total amount limiting and identifying link; if there is Lim _s <W _j And Lim _g ＜W _j The program enters the next calculation unit loop. When the water intake source belongs to the underground water source, if Lim exists _g ≥W _j The groundwater resource utilization control amount is reduced to Lim _g -W _j And the program enters a water total amount limiting and identifying link; if there is Lim _g <W _j And Lim _s ≥W _j The program searches for a secondary water supply priority water source until a water source belonging to the surface water type is found, and then enters a water total amount limiting and identifying link; if there is Lim _g <W _j And Lim _s ＜W _j The program enters the next calculation unit loop.

(5) Water total limit identification

The program sets a limit for total water consumption in the area ₃ . For water sources with water supply priority J (j=1, 2, …, o; o=i+j+k+l+m+n+x+y+z), if Lim ₃ ≥W _j The area water control amount is reduced to Lim ₃ -Wj, and the program enters a water source water intake calculation link; if Lim ₃ <W _j The program enters the next calculation unit loop.

(6) Water source water intake calculation

For water sources with water supply priority J (j=1, 2, …, o; o=i+j+k+l+m+n+x+y+z), if availableWater amount Wsc _j > WD, then the water supply to the water source is WSP _j =wd, and the irrigation water demand is reduced to 0, the program goes to the next water department cycle; otherwise, the water source supplies water to WSP _j =Wsc _j And the irrigation water demand is reduced to WD-Wsc _j The program will find the secondary water supply priority source and continue the water supply calculation until the irrigation water demand of the calculation unit is met or until the last water supply source is supplied. Calculating the irrigation water consumption

. The program goes to the next calculation unit loop.

(3) Water consumption and drainage module

The water consumption and drainage module receives daily water consumption data of water departments of the calculation unit output by the water consumption and quality configuration module, and obtains data such as urban and rural resident life water consumption rate, ecological environment water consumption rate, industrial water consumption rate, building water consumption rate, service industry water consumption rate, livestock water consumption rate and the like through reading input water consumption rate data files, calculates water consumption and drainage respectively, and transmits the water consumption and drainage to the pollution discharge module. The calculation formula is as follows:

WC _i =WU _i ·f _i 。

WR _i =WU _i -WC _i 。

Wherein WC _i For calculating the daily water consumption of the ith water department on the unit, the unit is m ³ /d；WU _i For calculating daily water consumption of ith water department (including urban and rural resident life, urban and rural ecology, industry, building industry, service industry, livestock and poultry raising and the like) on a unit, the unit is m ³ /d；WR _i For calculating the daily displacement of the ith water department on the unit, the unit is m ³ /d；f _i The water consumption rate (including urban and rural resident life water consumption rate, ecological environment water consumption rate, industrial water consumption rate, building water consumption rate, service water consumption rate, livestock water consumption rate and the like) of the ith water consumption department on the unit is calculated.

The water consumption and the water discharge of the farm irrigation are obtained by simulation of a soil moisture process through a distributed hydrological model, and the calculation formula is as follows:

W _aird =WU _irr ·(1-rto)。

W _qird =WU _irr ·rto。

wherein W is _aird For calculating the daily irrigation water consumption on the unit, the unit is m ³ /d；WU _irr For calculating the daily irrigation water consumption of the unit, the unit is m ³ /d；W _qird For calculating the irrigation displacement on the unit, the unit is m ³ /d; rto is the irrigation production rate on the calculation unit.

(4) Pollution discharge module

The drainage module receives drainage data output by the drainage consumption module, calculates point source pollution (urban living pollution and industrial pollution) and rural non-point source pollution (rural living pollution, livestock and poultry cultivation pollution and farmland pollution), and transmits the point source pollution and the rural non-point source pollution to the river channel module and the sub-river basin module of the distributed hydrological model respectively. The pollutants mainly comprise ammonia nitrogen, organic nitrogen, nitrate nitrogen, total phosphorus, biochemical oxygen demand and the like.

The calculation formula of the urban point source emission is as follows: NP (NP) _urban =(WR _urban +WR _in )·[(1-r)·c ₀ +r·c _d ]/1000。

Wherein NP _urban The unit of the pollution discharge amount is kg/d for the point source of the town; WR (WR) _urban And WR _in The unit of the discharge amount of the urban domestic sewage and the discharge amount of the industrial sewage is m ³ /d; r is the sewage treatment rate of a sewage treatment plant; c ₀ And c _d The concentration of the sewage pollutants is respectively the untreated sewage and the treated sewage, and the unit is mg/L.

The rural non-point source pollution amount expression is as follows: NP (NP) _rural =(P _rural ·m _rural +P _an ·m _an )·10。

Wherein NP _rural The unit of the non-point source pollution emission is kg/d in rural areas; m is m _rural And m _an The emission intensity of the living pollution of rural residents and the emission intensity of the livestock and poultry breeding pollutants are respectively in g/(d. People) or g/(d. People).

Step4: and modifying the distributed hydrologic model.

Distributed hydrologic model modification refers to code modification of related modules in the distributed hydrologic model. The modified content includes:

(1) Water source water taking program design

In order to realize the connection of the water source module of the distributed hydrologic model and the water quantity and quality configuration module, a calculation program with the following formula is added in the water source modules such as a river channel module, a reservoir module, a groundwater module, a pit module, an external water regulating module and the like, so as to finish the daily dynamic change simulation of the water intake of each water source in the configuration process.

WD1=WD0–min(WSP,Wsc)。

Wherein WD1 is the modified water demand of the department, and the unit is m ³ The method comprises the steps of carrying out a first treatment on the surface of the WD0 is the water demand of the department before modification, and the unit is m ³ The method comprises the steps of carrying out a first treatment on the surface of the WSP is water source and water supply quantity, and the unit is m ³ The method comprises the steps of carrying out a first treatment on the surface of the Wsc is the available water quantity of the water source, and the unit is m ³ 。

(2) Water intake restriction program design

(1) And adding the surface water supply limiting and restricting related codes into the river channel module, the reservoir module, the pit module, the reclaimed water module, the rain collecting module, the desalted water module and other modules, so that the surface water intake cannot exceed the total surface water utilization limiting amount. The related expression is:

。

wherein i1 is a day number; t1 is the number of simulated days; j1 is the calculation unit sequence number; m is the total number of calculation units; WSP (Wireless sensor protocol) _Si1j1 The unit is m, which is the actual water intake of the surface water resource on the i1 th day of the j1 st calculation unit ³ ；Lim _S The unit of the control quantity is m for the utilization of surface water resources ³ 。

(2) And adding an underground water supply quantity limiting constraint related code into the underground water module, so that the underground water exploitation quantity cannot exceed the total quantity of the underground water utilization limit. The related expression is:

。

wherein WSP _Gi1j1 For the j1 st meterCalculating the actual water intake of the groundwater resource on the i1 th day of the unit, wherein the unit is m ³ ；Lim _G The unit of the control quantity is m for the utilization of groundwater resources ³ 。

(3) Design of external water regulating program

And adding related codes into the external water regulating module to perfect the external water regulating and supplying function, so that the external water regulating consumption cannot exceed the total water regulating limit. The related expression is:

。

Wherein wout is _i1j1 The water supply amount of the external regulating water source of the jth 1 computing unit on the ith 1 day is expressed as m ³ ；out _MX To regulate the total amount of water, the unit is m ³ 。

(4) Unconventional water source programming

And designing an unconventional water source module such as a rain collecting module, a desalted water module, a regenerated water module and the like, adding a water supply calculation related code, and realizing daily dynamic change simulation of the water intake of unconventional water. The correlation formula is:

。

wherein wucw is _i1j1 The water supply amount of the j1 st calculation unit on the i1 st day is the unconventional water (rainwater, desalted water and regenerated water) with the unit of m ³ ；ucw _MX The water quantity is available for unconventional water, and the unit is m ³ 。

And adding the rain collecting module, the desalted water module and the regenerated water module into the refined configuration module for calling.

(5) Water quality calculation program design

And adding a calculation program with the following formula into the river channel module, the reservoir module, the pit module and the underground water module to calculate the water quality concentration of the water body in real time. The calculation formula is as follows: c=q×1000/W.

Wherein c is the water quality concentration of the water source, and the unit is mg/L; q is pollution load of a water source, and the unit is kg; w is the water resource amount of the water source, and the unit is m ³ 。

Step5: and (5) a module embedding design.

The module embedding design means to realize the embedding of national economy module, fine configuration module and distributed hydrologic model. The main content comprises:

(1) Parameters such as population number, three times of industry increment value, livestock number and the like in the national economy module are added into a water demand calculation program in the fine configuration module, so that the connection of the national economy module and the fine configuration module is realized.

(2) And shielding source codes of water source modules such as a river channel module, a reservoir module, a groundwater module, a pit module and the like in the distributed hydrological model, moving the modules into a fine configuration module, and calling each water source module by the fine configuration module.

(3) Relevant codes for calling the fine configuration module are added in the sub-basin module, so that the embedding of the distributed hydrologic model to the fine configuration module is realized, and the fine configuration simulation of water quantity and water quality is carried out when the hydrologic model runs.

Step6: and (5) superposing and calculating the design.

The superposition calculation design refers to space superposition calculation and daily-month-year-by-year time superposition calculation of a calculation unit, a water function area, an irrigation area, a administrative area, a sub-drainage basin and a full drainage basin, so that a fine configuration function is embodied. The main content comprises:

(1) Adding a daily scale superposition module in a daily circulation program of the distributed hydrologic model, and accumulating daily configuration results of a computing unit with the same spatial attribute to obtain daily scale water quantity and quality configuration results of a water functional area, a irrigated area, a administrative area, a sub-drainage area and a full-drainage area, wherein the daily scale water quantity and quality configuration results comprise daily water intake of each water source, daily water demand of each department, daily water consumption, daily water drainage and daily pollution discharge. The superposition formula is: w1 (X, i 1) =w0 (X, i 1) +m (j 1, i 1).

Wherein X is a water function area or a irrigated area or a administrative area or a sub-drainage basin or a drainage basin code; j1 is the calculation unit sequence number; w1 (X, i 1) is the day configuration result of the superimposed X on the i1 th day; w0 (X, i 1) is the day configuration result of X before superposition on the i1 th day; m (j 1, i 1) is the daily configuration result of the j1 st calculation unit on the i1 st day, and the daily configuration result comprises the daily water intake of a water source, the daily water demand of each department, the daily water intake, the daily water consumption, the daily water drainage and the daily pollution discharge.

(2) And adding a month scale superposition code into the day scale superposition module, and accumulating day-by-day configuration results of the computing units with the same spatial attribute to obtain month scale water quantity and water quality configuration results of the water functional area, the irrigation area, the administrative area, the sub-drainage area and the full drainage area. If the calculation unit j1 has an X attribute, the superposition formula is: m1 (X, mon) =m0 (X, mon) +w1 (X, i 1).

Wherein mon is a month number; j1 is the calculation unit sequence number; m1 (X, mon) is the month configuration result of the superposed X in the mon month; m0 (X, mon) is the month configuration result of X before superposition at mon month.

(3) And adding an annual scale superposition module in an annual circulation program of the distributed hydrologic model, and accumulating month-to-month configuration results of the computing units with the same spatial attribute to obtain annual scale water quantity and water quality configuration results of the water functional area, the irrigation area, the administrative area, the sub-watershed and the full watershed. The superposition formula is: n1 (X, yr 1) =n0 (X, yr 1) +m1 (X, mon).

Wherein yr1 is a year number; n1 (X, yr) is the year configuration result of superimposed X in yr1 year; n0 (X, yr 1) is the annual configuration result of X in yr1 before superposition, the annual configuration result including annual water intake of each water source, annual water demand of each department, annual water consumption, annual water discharge and annual pollution discharge.

Step7: and (5) storing and regressing the design.

The storage and regression design refers to the storage of simulation results of parameters which need to participate in the next cycle, and preparation is made for the simulation of the next cycle; and carrying out regression setting on the simulation results of the parameters which do not participate in the next cycle, and preventing cycle accumulation. The main content comprises:

(1) The storage design is mainly aimed at day scale parameters participating in the simulation of the next hydrologic cycle and associated process, including point source day drainage and pollution discharge. In order to connect the fine configuration module with the point source module, the design program adopts the wrong day transmission parameters to store data, and arranges the point source day drainage and the pollution discharge of the current day at the lower partIs discharged every day. The transfer formula is as follows:

。

wherein Point is the daily or pollution discharge of the Point source; point1 is the wrong day delivery parameter; j2 is a point source code; lst is the last day of the year.

(2) The regression design is mainly to carry out zero setting on parameters which do not participate in the next cycle.

(1) The daily scale configuration parameters which do not participate in the next daily cycle comprise daily water intake of a daily water source, daily water demand of departments, water consumption and the like. And uniformly zeroing the parameters after participating in month superposition calculation, and releasing the storage space.

(2) The configuration parameters of the month scale which do not participate in the next month circulation comprise the water intake of a month water source, the water consumption of a department month, the water consumption, the pollution discharge and the like. The parameters are uniformly zeroed after participating in the annual scale superposition calculation, and the accumulation of the month parameter values is prevented.

Thus, the construction of the model is completed.

More specifically, the above model will be described by taking a specific region as an example:

the invention has the beneficial effects that:

1. space refinement: based on the traditional water resource partition and administrative partition, the natural underlying surface condition in the partition and the space heterogeneity characteristic of the water used by the economy and society are fully considered, the calculation unit is further refined and configured, and the space downscaling method is adopted to simulate the space distribution conditions of regional water demand, water consumption, drainage, water shortage and pollution loads under the water resource configuration scheme.

2. Time refinement: the daily scale calculation step length is adopted to dynamically describe complex physicochemical processes such as regional water resource change, point source and non-point source pollutant migration and conversion, and the influence of water resource allocation on the water resource quantity, quality and spatial distribution of the water resource in the next period is reflected in real time.

3. Refining elements:

(1) The national economy development and the economic society water simulation are refined: in the water resource allocation process, the regional population distribution, the industrial layout and the spatial differences of water demand, water consumption and pollution emission are fully considered, and the regional economic and social development and the water resource utilization condition are truly reflected.

(2) Water source water supply simulation refinement: the multi-water source supply, multi-process regulation and control and multi-user dissipation functions of the water resource system are described, and the topological relation among water taking sources, water supply projects and water departments and the transfer relation of water quantity and water quality in the water resource system can be clearly and finely described.

(3) Water quantity and water quality simulation refinement: the influence of space-time distribution of water resources on water resource allocation is considered, and the influence of water supply according to different qualities is also considered, so that the requirements of different water departments on water quality of water sources are met, and the water resources with different water quality types are reasonably utilized.

Fig. 2 is a schematic structural diagram of a water quality joint configuration system based on distributed water quantity provided by the embodiment of the invention. As shown in fig. 2, the water quality joint configuration system based on distributed water quantity in this embodiment includes:

the sub-basin dividing module 201 is configured to divide the entire configuration area into a plurality of sub-basins.

The hydrological simulation unit dividing module 202 is configured to divide the sub-basin into a plurality of hydrological simulation units according to the underlying surface information of the sub-basin; the underlying information includes land use type and soil type; land use types include construction land and farmland land.

A calculation unit dividing module 203, configured to divide the hydrologic simulation unit into a plurality of calculation units according to the range of the subregion within the configuration region and the range of the hydrologic simulation unit; the subareas comprise county-level administrative areas, water function areas and irrigation areas; the calculation units are divided into town calculation units, rural calculation units, farmland calculation units, woodland calculation units, grassland calculation units and water area calculation units.

The daily water quality configuration module 204 is configured to determine daily water quality configuration data of each calculation unit, and determine daily water quality configuration data of the sub-area, daily water quality configuration data of the sub-basin and daily water quality configuration data of the configuration area based on the daily water quality configuration data of all calculation units by adopting a superposition method; the daily water quantity and water quality configuration data comprise: daily water intake of the water source, daily water demand, daily water consumption, daily water discharge and daily pollution discharge of each water department in the calculation unit; the water department includes: urban resident water sector, rural resident water sector, industrial water sector, building industry water sector, service industry water sector, livestock and poultry breeding water sector, urban ecological water sector, rural ecological water sector and irrigation water sector.

A lunar water quality configuration module 205, configured to determine lunar water quality configuration data of the computing unit based on solar water quality configuration data of the computing unit by adopting a superposition method, and determine lunar water quality configuration data of the sub-area, lunar water quality configuration data of the sub-basin and lunar water quality configuration data of the configuration area based on lunar water quality configuration data of all the computing units by adopting a superposition method; the lunar water quantity and water quality configuration data comprise: the month water intake of the water source and the month water demand, month water consumption, month water discharge and month pollution discharge of each water using department in the calculating unit.

The annual water quantity water quality configuration module 206 is configured to determine annual water quantity water quality configuration data of the calculation unit based on monthly water quantity water quality configuration data of the calculation unit by adopting a superposition method, and determine annual water quantity water quality configuration data of the sub-zone, annual water quantity water quality configuration data of the sub-basin and annual water quantity water quality configuration data of the configuration region based on annual water quantity water quality configuration data of all the calculation units by adopting a superposition method; the annual water quantity and water quality configuration data comprise: annual water intake of the water source, annual water demand, annual water consumption, annual water discharge and annual pollution discharge of each water department in the calculation unit.

In terms of the determination process of the daily water volume and water quality configuration data for any computing unit, the daily water volume and water quality configuration module 204 specifically includes:

the water consumption parameter and land utilization type determining submodule is used for determining water consumption parameters and land utilization types of all water consumption departments of the computing unit; the water consumption parameters include: the water demand, the water use sequence and the upper limit of the total water use amount of the allocation area, and land use types include construction land and farmland land.

The water supply parameter determination submodule is used for determining water supply parameters of all water sources; the water source includes: river course, reservoir, shallow aquifer, deep aquifer, pit pool, external water source, regeneration water factory, rain collecting pool and desalination water factory, water supply parameters include: water supply sequence, water quality, water supply availability, water supply availability, surface water resource utilization control and groundwater resource utilization control.

The daily water quantity and water quality configuration submodule is used for determining the daily water quantity of the corresponding water department according to the water consumption parameter, the land utilization type and the water supply parameter of each water supply source so as to obtain daily water quantity and water quality configuration data of each calculation unit.

As an alternative embodiment, the solar water volume and water quality configuration submodule specifically includes:

The configuration mode determining unit is used for determining the configuration mode according to the land utilization type of the calculating unit.

The first configuration unit is used for carrying out water quantity and water quality configuration on each calculation unit according to water consumption parameters and water supply parameters of each water supply source by adopting a first configuration mode when the land utilization type is a construction land, wherein the calculation unit comprises a town resident water department, a rural resident water department, an industrial water department, a building industry water department, a service industry water department, a livestock and poultry cultivation water department, a town ecological water department and a rural ecological water department.

And the second configuration unit is used for carrying out water quantity and water quality configuration on each calculation unit according to the water consumption parameter and the water supply parameter of each water supply source by adopting a second configuration mode when the land utilization type is farmland land.

As an optional implementation manner, the first configuration unit specifically includes:

the first configuration subunit is used for sequentially configuring the water quantity and the water quality of each calculation unit according to the priority of each calculation unit.

As an optional implementation manner, the second configuration unit specifically includes:

and the second configuration subunit is used for sequentially configuring the water quantity and the water quality of each calculation unit according to the priority of each calculation unit.

In the present specification, each embodiment is described in a progressive manner, and each embodiment is mainly described in a different point from other embodiments, and identical and similar parts between the embodiments are all enough to refer to each other. For the system disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant points refer to the description of the method section.

The principles and embodiments of the present invention have been described herein with reference to specific examples, the description of which is intended only to assist in understanding the methods of the present invention and the core ideas thereof; also, it is within the scope of the present invention to be modified by those of ordinary skill in the art in light of the present teachings. In view of the foregoing, this description should not be construed as limiting the invention.

Claims

1. A distributed water volume and quality configuration method, the method comprising:

dividing the whole configuration area into a plurality of sub-domains;

2. A distributed water volume water quality configuration system, the system comprising:

the sub-drainage basin dividing module is used for dividing the whole configuration area into a plurality of sub-drainage basins;

the hydrologic simulation unit dividing module is used for dividing the sub-watershed into a plurality of hydrologic simulation units according to the underlying surface information of the sub-watershed; the underlying information includes land use type and soil type; the land use types include construction land and farmland land;

the calculation unit dividing module is used for dividing the hydrologic simulation unit into a plurality of calculation units according to the range of the subareas in the configuration area and the range of the hydrologic simulation unit;

the daily water quantity and water quality configuration module is used for determining daily water quantity and water quality configuration data of each calculation unit, and determining daily water quantity and water quality configuration data of the subareas, daily water quantity and water quality configuration data of the sub-watershed and daily water quantity and water quality configuration data of the configuration area by adopting a superposition method based on the daily water quantity and water quality configuration data of all calculation units; the daily water quantity and water quality configuration data comprises: daily water intake of the water source, daily water demand, daily water consumption, daily water discharge and daily pollution discharge of each water department in the calculation unit;

The lunar water quality configuration module is used for determining lunar water quality configuration data of the computing units based on solar water quality configuration data of the computing units by adopting a superposition method, and determining lunar water quality configuration data of the subareas, lunar water quality configuration data of the sub-basins and lunar water quality configuration data of the configuration areas based on lunar water quality configuration data of all the computing units by adopting a superposition method; the lunar water quantity and water quality configuration data comprises: the month water intake of the water source, the month water demand, month water consumption, month water discharge and month pollution discharge of each water department in the calculation unit;

the annual water quantity water quality configuration module is used for determining annual water quantity water quality configuration data of the calculation units based on monthly water quantity water quality configuration data of the calculation units by adopting a superposition method, and determining annual water quantity water quality configuration data of the subareas, annual water quantity water quality configuration data of the sub-watershed and annual water quantity water quality configuration data of the configuration area based on annual water quantity water quality configuration data of all calculation units by adopting a superposition method; the annual water quantity and water quality configuration data comprises: annual water intake of the water source, annual water demand, annual water consumption, annual water discharge and annual pollution discharge of each water department in the calculation unit;

In the aspect of determining daily water quantity and water quality configuration data of any computing unit, the daily water quantity and water quality configuration module specifically comprises:

the water consumption parameter and land utilization type determining submodule is used for determining water consumption parameters and land utilization types of all water consumption departments of the computing unit; the water consumption parameters include: the water demand, the water sequence and the upper limit of the total water amount of the allocation area, wherein the land utilization type comprises construction land and farmland land;

the water supply parameter determination submodule is used for determining water supply parameters of all water sources; the water supply parameters include: water supply sequence, water quality, water supply availability, water supply availability, surface water resource utilization control amount, and groundwater resource utilization control amount;

and the daily water quantity and water quality configuration submodule is used for determining the daily water quantity of the corresponding water department according to the water consumption parameter, the land utilization type and the water supply parameter of each water supply source so as to obtain daily water quantity and water quality configuration data of each calculation unit.