DESIGN METHOD FOR DISTRIBUTED HYDROLOGICAL CYCLE MODEL BASED ON MULTI-SOURCE COMPLEMENTARY WATER SUPPLY MODE

Abstract

The present disclosure provides a design method for a distributed hydrological cycle model based on a multi-source complementary water supply mode, the method including the following steps: S1, nested hydrological response unit (HRU) division; S2, HRU attribute design; S3, design of a multi-source complementary water supply module; and S4, improvement on a SWAT model. Based on the Soil and Water Assessment Tool (SWAT) model, the present disclosure develops a distributed natural-artificial hydrological dynamic reciprocal simulation model. The model is endowed with the functions of simulating dynamic reciprocation of natural water cycle and artificial water cycle, and integration of development, utilization and regulation of water resources, thereby simulating a natural-artificial hydrological cycle based on modes of urban multi-source water supply and multi-source irrigation water supply.

Claims

1. A design method for a distributed hydrological cycle model based on a multi-source complementary water supply mode, comprising the following steps: Step S1: conducting hydrological response unit (HRU) division by adopting a nested slope discretization method based on attributes of “basin, water resources region, administrative region, irrigation area, land use, soil, slope”, wherein HRUs obtained after division each have corresponding attributes; Step S2: constructing an HRU attribute recognition module, wherein the HRU attribute recognition module is configured to recognize attributes of an HRU; Step S3: designing a multi-source complementary water supply module, wherein the multi-source complementary water supply module is configured to invoke the HRU attribute recognition module to recognize the attributes of each HRU, determining a land use type, a corresponding water source and a water supply priority of the water source according to the recognized attributes, and invoking, by the water supply priority of the water source, a corresponding water source module to conduct water withdrawal; and Step S4: connecting the multi-source complementary water supply module with modules modified in a Soil and Water Assessment Tool (SWAT) model to realize real-time data exchange, wherein the HRU allocates and regulates water resources according to input information about water demand, types of water sources, rules of water supply priority, and water conservancy projects, and information about natural hydrological conditions which is provided by the SWAT model, and outputs and transfers information about an artificial hydrological cycle regarding daily “water supply, water use, water consumption, water drainage” to the SWAT model, wherein step S3 specifically comprises: designing a water source code information file, wherein the water source code information file is configured to read designated water source information, six types of water sources are set and comprise transferred water, reservoir water, urban river water, shallow groundwater, deep groundwater and pit-pond water, and the water source code information file is read by program instructions; designing a water supply priority information file, wherein the water supply priority information file is configured to read information about water supply priority, and specify a water supply priority of a water source, and is read by program instructions; designing a water withdrawal control information file, wherein the water withdrawal control information file is configured to read information about water supply control volume, and recognize an annual surface water supply control volume and an annual groundwater exploitation control volume of an administrative region to which the HRU belongs for the subsequent calculation of water withdrawal volume of water sources; and designing a calculation process for multi-source complementary water supply, wherein the specific calculation process is as follows: first, recognizing the land use type of the HRU, wherein if it is construction land, a program enters a calculation process for urban and rural water supply; if it is agricultural land, the program enters a calculation process for irrigation water; and if it is other land use type, the program ends; invoking a corresponding water source module by recognizing a water withdrawal source identification code of the HRU, recognizing the number, type and water withdrawal sequence of water sources of each HRU by reading the water source code information file and water supply priority information file, and invoking each water source module in turn according to the water withdrawal source identification code; and seeking water sources and conducting water withdrawal from each water source according to a water supply sequence of the HRU until the HRU’s daily demand for domestic water, industrial water, and agricultural irrigation water is satisfied, or until the last water source finishes water supply; and the calculation for multi-source complementary water supply comprises the following steps: specifying a daily water demand WD set by a target HRU; specifying the number k, water source codes and water supply priority of water sources of the target HRU, wherein k≤30; invoking the water source modules in sequence to calculate a water withdrawal volume of a water source, wherein the water source modules comprise a rchuse module, a res module, watuse module, an irr rch module, an irr res module and an irrsub module, the water withdrawal volume of the water source depends on a daily water demand of the HRU and an available water supply of the water source, while the available water supply depends on an accessible water volume of the water source, the water supply capacity of a water withdrawal project and the water withdrawal control volume, wherein the calculation formulas are as follows: $WS{P}_{ij}=min\left(W{D}_i-{\displaystyle \underset{k=1}{\overset{j-1}{.Math.}}WS{P}_k,Ws{c}_{ij}}\right)$ $Ws{c}_{ij}=min\left(W{A}_{ij},W{F}_i,WM{X}_{ij}-{\displaystyle \underset{k=1}{\overset{j-1}{.Math.}}WS{P}_k}\right)$ $WM{X}_{ij}=\left\{\begin{array}{l}min\left(WUM-{\displaystyle \underset{m=1}{\overset{i-1}{.Math.}}{\displaystyle \underset{k=1}{\overset{j}{.Math.}}WS{P}_k,WSM-{\displaystyle \underset{m=1}{\overset{i-1}{.Math.}}{\displaystyle \underset{k=1}{\overset{j}{.Math.}}WS{P}_k}}}}\right)\\ min\left(WUM-{\displaystyle \underset{m=1}{\overset{i-1}{.Math.}}{\displaystyle \underset{k=1}{\overset{j}{.Math.}}WS{P}_k,WGM-{\displaystyle \underset{m=1}{\overset{i-1}{.Math.}}{\displaystyle \underset{k=1}{\overset{j}{.Math.}}WS{P}_k}}}}\right)\end{array}\right)$ wherein, i denotes a sequence number of an HRU; j denotes a water supply priority number of a water source; WSP denotes an actual daily water withdrawal (m.sup.3) of a water source; WD indicates a daily water demand (m.sup.3) of an HRU; Wsc indicates a daily available water supply (m.sup.3) of a water source; WF indicates the water supply capacity (m.sup.3) of a water withdrawal project; WA indicates a daily accessible water volume of a water source (m.sup.3) ; WMX denotes an annual water withdrawal control volume (m.sup.3) , and WUM denotes an annual water consumption control volume (m.sup.3) ; WSM denotes an annual surface water withdrawal control volume (m.sup.3) ; and WGM denotes an annual groundwater exploitation control volume (m.sup.3) ; wherein for a water source with a water supply priority of 1, priority is given to water withdrawal from the water source; if the available water supply of the water source is Wsc.sub.1>WD, then the water supply of the water source is WSP.sub.1=WD, a water supply program ends, and the total water supply of the water source of the HRU is WSP= WSP.sub.1; otherwise, WSP.sub.1=Wsc.sub.1, and the water demand of the HRU changes to Wf =WD-Wsc.sub.1, and the program will continue to seek the next grade of water source; for a water source with a water supply priority of j, j=2, ..., k-1; k≤30, if the daily available water supply of the water source is Wsc.sub.j>Wf, then the water supply of the water source is WSP.sub.j=Wf, the programends, and the total water supply of the water source the HRU is SP=WSP+WSP.sub.j; otherwise, WSP.sub.j=Wsc.sub.j, the water demand of the HRU changes to Wf=Wf-Wsc.sub.j, and the program will continue to seek the next grade of water source; and for a water source with a water supply priority of k, k≤30, if the daily available water supply of the water source is Wsc.sub.k>Wf, then the water supply of the water source is WSP.sub.k=Wf, the program ends, and the total water supply of the HRU is WSP=ΣWSP.sub.i; otherwise, WSP.sub.k=Wsc.sub.k, the water demand of the HRU changes to Wf=Wf-Wsc.sub.k, and the program ends.

2. The design method for a distributed hydrological cycle model based on a multi-source complementary water supply mode according to claim 1, wherein step S1 comprises: extracting a river network of a basin from a DEM using ArcGIS to conduct division to obtain natural sub-basins; superimposing land use information, soil type information and slope information on the natural sub-basins to conduct division to obtain natural HRUs; setting boundaries of an administrative region and a water resources region for the natural HRU to further divide the natural HRUs; and superimposing irrigation areas with the natural HRUs according to the distribution of the irrigation areas to finally complete HRU division, wherein each HRU has a sub-basin attribute, a water resources region attribute, an administrative region attribute, an irrigation area attribute, a land use type attribute and a soil type attribute.

3. The design method for a distributed hydrological cycle model based on a multi-source complementary water supply mode according to claim 2, wherein step S2 comprises: constructing the HRU attribute recognition module which is configured to read specified HRU attributes, wherein the specified HRU attributes comprise a sub-basin attribute, a water resources region attribute, an administrative region attribute, and an irrigation area attribute; and putting the constructed HRU attribute recognition module in a main module in the SWAT model to facilitate invocation of the HRU attribute recognition module.

4-5. (canceled)

6. The design method for a distributed hydrological cycle model based on a multi-source complementary water supply mode according to claim 1, wherein modification for the relevant modules in the SWAT model in step S4 specifically comprises: shielding the rchuse module, the res module, the watuse module, the irr_rch module, the irr_res module and the irrsub module, and forgoing adopting a single water source withdrawal mode; and putting the foregoing modules into the multi-source complementary water supply module for invoking; adding relevant codes, and replacing parameters waterrch and wuresn in the rchuse module and the res module with parameter WSP.sub.i, respectively to achieve connection of the multi-source complementary water supply module Multi_sc with the rchuse module and the res module as well as invoking, wherein i =1,2; modifying relevant programs to add functions of transferred water withdrawal and transferred water volume restriction so as to control water supply within a total transferred water limit, wherein a calculation formula is expressed as follows: ${\displaystyle \underset{i}{.Math.}{\displaystyle \underset{j}{.Math.}\text{waterout}\left(i,j\right)\le MX5}}$ wherein waterout (i,j) denotes transferred water consumption (m.sup.3) of the jth HRU on the ith day; and WX5 denotes total transferred water limit (m.sup.3) ; adding codes in the watuse module, and replacing parameters watershal, waterdeep, waterout, and waterpnd in the watuse module with parameter WSP.sub.i, respectively to achieve connection of the multi-source complementary water supply module Multi_sc with the watuse module as well as invoking, wherein i=3,4,5,6; adding a calculation program of the following formula in the rchuse module, the res module, and the watuse module: $\text{WSP}=\text{WSP}\cdot \left(1-pip\right)$ wherein, pip denotes a leakage rate of water supply pipe network; adding a pit-pond irrigation function, completing a transferred water irrigation function, and imposing water supply restriction to control an irrigation water withdrawal within the total transferred water limit: ${\displaystyle \underset{i}{.Math.}{\displaystyle \underset{j}{.Math.}\text{wirrout}\left(i,j\right)\le MX5}}$ ${\displaystyle \underset{i}{.Math.}{\displaystyle \underset{j}{.Math.}wirrpnt\left(i,j\right)\le MX6}}$ wherein, wirrout (i,j) denotes transferred water irrigation consumption (m.sup.3) of the jth HRU on the ith day; WX5 denotes total transferred water limit (m.sup.3) , and wirrpnt (i,j) denotes pit-pond irrigation consumption (m.sup.3) of the jth HRU on the ith day; and WX6 denotes pit-pond available water supply (m.sup.3) ; modifying source codes of the irr_rch module, the irr_res module, and the irrsub module to add simulation on a channel system delivery process comprising channel water loss and channel recession, wherein the channel water loss comprises two parts of channel water evaporation loss and channel leakage loss, and the main calculation formulas are as follows: $E{T}_{can}=IR{R}_{can}\cdot \left(1-\phi \right)\cdot \alpha$ $L{s}_{can}=IR{R}_{can}\cdot \left(1-\phi \right)\cdot \beta$ $Sur{p}_{can}=IR{R}_{can}\cdot \left(1-\phi \right)\cdot \left(1-\alpha -\beta \right)$ wherein, ET.sub.can denotes a channel system evaporation loss (mm); IRR.sub.can denotes an irrigation water volume (mm) entering a channel; Ls.sub.can denotes a channel system leakage loss (mm) ; Surp.sub.can denotes a channel system recession volume (mm); φ denotes an effective utilization coefficient of channel system water; α denotes a channel system evaporation coefficient; and β denotes a channel system leakage coefficient; adding a calculation program for leakage loss by modifying relevant codes of a percmain module, wherein a calculation formula is as follows: $Ws{l}_{yr1,t}=Ws{l}_{yr1,t-1}+in{f}_{pcp}+in{f}_{irr}+in{f}_{wet}+L{s}_{can}$ wherein, .sup.Wslyr1,t+1 denotes soil water content (mm) of a first layer of soil on the t-th day; .sup.Wslyr1,t denotes soil water content (mm) of a first layer of soil on the (t-1) th day; inf.sub.pcp denotes precipitation infiltration capacity (mm); inf.sub.irr denotes irrigation infiltration capacity (mm); and inf.sub.wet denotes lake and reservoir wetland infiltration capacity (mm); modifying groundwater recharge codes in a gwmod module to achieve simulation on water leakage of a pipe network, wherein a calculation formula is as follows: $\begin{array}{l}r{h}_{\epsilon }=\left(1-\exp \left(-1/GW\_DELAY\right)\right)\cdot \left(prc+WSP\cdot pip/Area\right)\\ +\exp \left(-1/GW\_DELAY\right)\cdot r{h}_{L-1}\end{array}$ wherein, .sup.rht denotes groundwater recharge capacity (mm) on the t-th day; .sup.rht-1 denotes groundwater recharge capacity (mm) on the (t-1) th day; prc denotes soil water leakage (mm) of recharged groundwater; GW_DELAY denotes groundwater recharge delay coefficient (mm) ; and Area denotes the area (m.sup.2) of an HRU; adding the multi-source complementary water supply module in a subbasin module, and conducting in-year dynamic complementary water supply operation on water sources by reading specified type, number, water source codes, water withdrawal volume, and water withdrawal time of water sources to achieve multi-source combined water supply simulation during the running of the SWAT model; superimposing channel system recession with earth surface runoff by modifying relevant codes of a surface module to participate in calculation of flow concentration of river channels, wherein a calculation formula is as follows: $sur{f}_t=sur{f}_0+Sur{p}_{can}$ wherein, .sup.surft denotes runoff (mm) after channel recession; and .sup.surf0 denotes runoff (mm) before channel recession; wherein a point source module comprises a recday module and a recmon module, wherein relevant codes are modified in the recday module and the recmon module, a pollution discharge parameter WDR is used to replace parameters floday and flomon, respectively, and the calculation formulas are as follows: $\text{WP}\text{=}\text{WSP}\cdot \left(1-\text{r}\right)$ $\text{WDR}\text{=}\text{WP}\cdot \left(1-v\right)+\text{WP}\cdot \text{v}\cdot \left(1-re\right)$ wherein WDR denotes urban sewage output (m.sup.3) ; WP denotes sewage discharge (m.sup.3) ; r denotes a water consumption rate; v denotes a sewage disposal rate of a sewage disposal plant; and re denotes a reclaimed water utilization rate; and putting the constructed HRU attribute recognition module in the main module in the SWAT model to facilitate invocation of the HRU attribute recognition module.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0071] FIG. 1 is a flowchart illustrating a design method for a distributed hydrological cycle model based on a multi-source complementary water supply mode according to an embodiment of the present disclosure;

[0072] FIG. 2 is a flowchart illustrating HRU division according to an embodiment of the present disclosure;

[0073] FIG. 3 is a running flowchart of a multi-source complementary water supply module according to an embodiment of the present disclosure;

[0074] FIG. 4 is a flowchart illustrating calculation of water supply priority according to an embodiment of the present disclosure;

[0075] FIG. 5 shows comparison between actual measurement of a Xindianpu station and a simulated monthly runoff according to an embodiment of the present disclosure;

[0076] FIG. 6 is a relation diagram illustrating transformation of natural-artificial hydrological cycle of Baihe Basin in Level Year 2016; and

[0077] FIG. 7 is a diagram illustrating in-year water supply process of various water resources of Baihe Basin in Level Year 2016.

DETAILED DESCRIPTION OF THE EMBODIMENTS

[0078] In order to make the objectives, technical solutions and advantages of the examples of the present disclosure clearer, the following clearly and completely describes the technical solutions in the examples of the present disclosure with reference to accompanying drawings in the examples of the present disclosure. Apparently, the described examples are some rather than all of the examples of the present disclosure. All other embodiments obtained by those of ordinary skill in the art based on the embodiments of the present disclosure without making inventive efforts shall fall within the scope of protection of the present disclosure.

[0079] Embodiments of the present disclosure provide a design method for a distributed hydrological cycle model based on a multi-source complementary water supply mode, the method mainly including the following four steps: nested HRU division; HRU attribute design, design of a multi-source complementary water supply module; and improvement on a SWAT model. The design of a multi-source complementary water supply module refers to the process of designing and adding relevant modules to add the functions of urban and rural multi-source water supply and multi-source irrigation simulation to the SWAT model, and improve artificial hydrological cycle, thus building a natural-artificial distributed hydrological cycle model. Improvement on the SWAT model refers to the process of modifying the relevant modules of the SWAT model to facilitate the reading of information about types, number, water source codes, and water withdrawal priority of water sources and the invocation of multi-source complementary module. Modules involved mainly include rchuse module, res module, watuse module, irr_rch module, irr_res module, irrsub module, subbasin module, percmain module, gwmod module, surface module, recday module, recmon module, main module, etc.

[0080] The method of the present includes the following steps: [0081] Step S1: nested HRU division; [0082] Step S2: HRU attribute design; [0083] Step S3: design of a multi-source complementary water supply module; and [0084] Step S4: improvement on a SWAT model.

[0085] In step S1: HRU division is conducted by adopting a nested slope discretization method based on attributes of “basin, water resources region, administrative region, irrigation area, land use, soil and slope”. As shown in FIG. 2, step S1 mainly includes: (1) division of natural sub-basins; (2) division of natural HRUs; (3) superimposition of natural HRUs with an administrative unit; and (4) superimposition of natural HRUs with irrigation areas. [0086] (1) Division of natural sub-basins is conducted by adopting Arc Hydro Tools in ArcGIS. Firstly, the correction operations such as depression filling are carried out on the DEM base map, and according to the set threshold requirements for generating Critical SourceArea (CSA) of sub-basins, the flow direction of a grid is determined, the waterline of the basin is identified, and division of the natural sub-basins is achieved based on the characteristic parameters of a river network such as slope, aspect and slope length of the basin; [0087] (2) Superimposition is conducted according to the land use type, soil type and slope type in sub-basins to conduct division to obtain a plurality of natural HRUs. HRU refers to a region where underlying surfaces have relatively monotonous and uniform features, that is, underlying surfaces in this region have identical hydrological characteristics, and there is a merely a combination of one kind of plantation cover, one kind of soil and one kind of slope type within each HRU. [0088] (3) A nest mode of superimposing HRUs with a prefecture-level / county-level administrative region and water resources region is adopted. By adopting the “EditorToolbar” function of GIS, a boundary GIS map of an administrative unit is directly superimposed onto a natural HRU map as obtained through division in the previous step. In a boundary region between the natural HRU and the management division, the natural HRU is divided into two parts according to the boundary line, and the HRUs obtained after division are endowed with attributes of the administrative unit; [0089] (4) By adopting the “EditorToolbar” function of GIS, in a boundary region between the natural HRU and an irrigation area, the natural HRU is divided into two parts according to the boundary line, then the HRUs obtained after division are endowed with attributes of the irrigation area, thus completing division the HRU. Finally, each HRU has a sub-basin attribute, a water resources region attribute, an administrative region attribute, an irrigation area attribute, a land use type attribute and a soil type attribute.

[0090] Step S2 of HRU attribute design includes: [0091] (1) construction of an HRU attribute recognition module readattr, which specifically includes: [0092] construct the HRU attribute recognition module readattr which is configured to read specified HRU attributes. Set HRU attributes are a sub-basin attribute, a water resources region attribute, an administrative region attribute, and an irrigation area attribute, respectively. The format of an input document can be found in Table 1. [0093] (2) invoking the HRU attribute recognition module, which specifically includes [0094] putting the constructed readattr module (HRU attribute recognition module) in the main module in the SWAT model to facilitate invocation of the readattr module.

[0095] The design of a multi-source complementary water supply module in step S3 includes: (1) designing a water source code information file; (2) designing a water supply priority information file; and (3) designing calculation for multi-source complementary water supply. [0096] (1) designing a water source code information file: [0097] the water source code information file is configured to read information about a specified water source. Six types of water sources are set and include transferred water, reservoir water, urban river water, shallow groundwater, deep groundwater and pit-pond water, and it is specified that every single HRU has at most 5 water sources of the same kind, that is, at most 30 water sources can be disposed for each HRU. Water source codes are specified in a way of reading an input document based on different types of HRUs; for residential and industrial land HRUs, specification is conducted with the administrative region as a unit, and HRUs within the same administrative region have consistent water sources; and for irrigation area HRUs, specification is conducted with the administrative region as a unit, and HRUs within the same irrigation area have consistent water sources. The format of the input document can be found in Table 2 and Table 3, respectively. [0098] (2) designing a water supply priority information file: [0099] the water supply priority information file is configured to read information about a specified water supply priority of a water source. Each water source is accessible to a plurality of HRUs, and meanwhile, each HRU can withdraw water from a plurality of water sources. A water supply sequence is set according to the type and number of water sources for each HRU. The water supply priority of water sources is specified in a way of reading an input document based on different types of HRUs; for residential and industrial land HRUs, specification is conducted with the administrative region as a unit, and HRUs within the same administrative region have consistent water supply priority; and for irrigation area HRUs, specification is conducted with the administrative region as a unit, and HRUs within the same irrigation area have consistent water supply priority. The priority is set from 1 to k (k denotes the number of water sources for HRU, k≤30), 1 denotes the highest priority, and the sequence number is set to 0 for unspecified water sources. The format of the input document can be found in Table 4 and Table 5. [0100] (3) designing a water withdrawal control information file: [0101] which is configured to read information about annual water supply control volume of each administrative region. Each HRU recognizes an annual surface water supply control volume and an annual groundwater exploitation control volume of an administrative region to which the HRU belongs based on attributes of the administrative region for the subsequent calculation of water withdrawal volume of water sources. The format of an input document can be found in Table 6. [0102] (4) designing calculation for multi-source complementary water supply: [0103] as shown in FIG. 3, the running process of the multi-source complementary water supply module Multi _sc mainly includes: 1) recognizing types of HRUs; 2) recognizing the water source and priority thereof; and 3) calculating multi-source complementary water supply. [0104] 1) recognizing types of HRUs, which specifically includes [0105] first, recognizing the land use type of an HRU, where if it is construction land (urban land and rural land), a program enters a calculation process for urban and rural water supply; if it is irrigation land, the program enters a calculation process for irrigation water; and if it is other land use type, the HRU has no requirement for water supply, and the program ends. [0106] 2) recognizing the water source and priority thereof, which specifically includes [0107] using the multi-source complementary water supply mode to revoke corresponding water source modules by recognizing water withdrawal source identification codes of an HRU, where the water withdrawal source identification codes for urban river water, reservoir water, shallow groundwater, deep groundwater, transferred water and pit-pond water are 1, 2, 3, 4, 5, and 6, respectively. The number, type and water withdrawal sequence of water sources of each HRU are recognized by reading the water source code information file and water supply priority information file, and each water source module is revoked in turn according to the water withdrawal source identification code. [0108] 3) calculating multi-source complementary water supply, which specifically includes [0109] using the multi-source complementary water supply module to seek water sources and conducting water withdrawal from each water source according to a water supply sequence of the HRU until the HRU’s daily demand for domestic water (construction land HRU), industrial water, and agricultural irrigation water (agricultural land HRU) is satisfies or until the last water source finishes water supply. As shown in FIG. 4, calculation of water supply priority includes: [0110] 1) specifying a daily water demand WD set by a target HRU; [0111] 2) recognizing the number k (k≤30), water source codes and water supply priority of water sources of the target HRU; and [0112] 3) invoking water source modules in sequence to calculate the water withdrawal of the water sources. The water source modules include a rchuse module, a res module and a watuse module. The water withdrawal volume of the water source depends on a daily water demand of the HRU and an available water supply of the water source. The available water supply depends on an accessible water volume of the water source, the water supply capacity of a water withdrawal project (such as water diversion channels, water supply pipelines, and power-driven wells) and the water withdrawal control volume. The calculation formulas are as follows: where, i denotes a sequence number of an HRU; j denotes a water supply priority number of a water source; WSP denotes an actual daily water withdrawal volume (m.sup.3) of a water source; WD indicates a daily water demand (m.sup.3) of an HRU; Wsc indicates a daily available water supply (m.sup.3) of a water source; WF indicates the water supply capacity (m.sup.3) of a water withdrawal project; WA indicates a daily accessible water volume of a water source (m.sup.3); WMX denotes an annual water withdrawal control volume (m.sup.3), and WUM denotes an annual water consumption control volume (m.sup.3); WSM denotes an annual surface water withdrawal control volume (m.sup.3); and WGM denotes an annual groundwater exploitation control volume (m.sup.3). [0113] 4) For a water source with a water supply priority of 1, priority is given to water withdrawal from the water source. If the available water supply of the water source is Wsc.sub.1>WD, then the water supply of the water source is WSP.sub.1=WD, a water supply program ends, and the water supply of the water source of the HRU is WSP= WSP.sub.1; otherwise, WSP.sub.1=Wsc.sub.1, and the water demand of the HRU changes to Wf =WD-Wsci, and the program will continue to seek the next grade of water source; [0114] 5) for a water source with a water supply priority of j (j=2, ..., k-1; k≤30), if the daily available water supply of the water source is Wsc.sub.j>Wf, then the water supply of the water source is WSP.sub.j=Wf, the program ends, and the total water supply of the water source the HRU is SP=WSP+WSP.sub.j; otherwise, WSP.sub.j=Wsc.sub.j, the water demand of the HRU changes to Wf=Wf-Wsc.sub.j, and the program will continue to seek the next grade of water source; and [0115] 6) for a water source with a water supply priority of k (k≤30), if the daily available water supply of the water source is Wsc.sub.k>Wf, then the water supply of the water source is WSP.sub.k=Wf, the program ends, and the total water supply of the HRU is WSP=.Math.WSP.sub.i; otherwise, WSP.sub.k=Wsc.sub.k, the water demand of the HRU changes to Wf=Wf-Wsc.sub.k, and the program ends.

[0116] In step S4, improvement on the SWAT model mainly involves connecting the multi-source complementary water supply module with modified modules in a SWAT model to realize real-time data exchange, where the HRU allocates and regulates water resources according to input information about water demand, types of water sources, rules of water supply priority, and water conservancy projects, and information about natural hydrological conditions which is provided by the SWAT model, and outputs and transfers information about an artificial hydrological cycle regarding daily “water supply, water use, water consumption, water drainage” to the SWAT model. The relevant modules of the SWAT model are modified to facilitate the reading of information about types, number, water source codes, and water withdrawal priority of water sources and the invocation of multi-source complementary module. The modified modules mainly include rchuse module, res module, watuse module, irr_rch module, irr_res module, irrsub module, subbasin module, surface module, gwmod module, recday module, recmon module, main module, etc. The main functions of each module are as follows:

[0117] rchuse module: the river channel water supply module, which is configured to conduct water withdrawal from a designated river channel for domestic and industrial water.

[0118] res module: a reservoir water supply module, which is configured to conduct water withdrawal from a designated reservoir for domestic and industrial water.

[0119] watuse module: a groundwater and pit-pond water supply module, which is configured to conduct water withdrawal from a shallow underground aquifer, a deep underground aquifer and pit-ponds in a designated sub-basin for domestic and industrial water.

[0120] irr_rch module: a river channel irrigation module, which is configured to conduct water withdrawal from a designated river channel for farmland irrigation.

[0121] irr_res: a reservoir irrigation module, which is configured to conduct water withdrawal from a designated reservoir for farmland irrigation.

[0122] irrsub module: a groundwater and transferred water irrigation module, which is configured to conduct water withdrawal from a shallow underground aquifer, a deep underground aquifer and an external water source in a designated sub-basin.

[0123] percmain module: a leakage loss module, which is configured to simulate moisture change in soil according to the calculation results of precipitation infiltration, irrigation infiltration and pit-pond infiltration.

[0124] gwmod module: a groundwater module, which is configured to simulate the variation in shallow groundwater and deep groundwater according to the results of surface water infiltration.

[0125] subbasin module: a sub-basin module, which mainly conducts evaporation simulation, runoff producing simulation, groundwater simulation, vegetation growth simulation, farmland management simulation, nutrient migration and transformation simulation, etc., which is realized by invoking modules related to the hydrological and water-quality process.

[0126] surface module: an earth surface runoff simulation module, which is configured to simulate hydrological processes of primary surface water such as canopy interception, accumulated snow and snowmelt, and earth surface runoff.

[0127] recday module: a daily point source module, which is configured to simulate the changing process of daily point sources by inputting daily pollutant emissions.

[0128] recmon module: a monthly point-source module, which is configured to simulate the changing process of monthly point sources by inputting monthly pollutant emissions.

[0129] main module; a main module of the SWAT model, which is mainly configured to read input documents, allocating array size, initializing parameters, simulating the hydrological process and so on.

[0130] Modification of relevant modules in the SWAT model specifically includes:

Shielding Modules

[0131] Shield the water source modules, namely the rchuse module, the res module, the watuse module, the irr_rch module, the irr_res module and the irrsub module, and forgo adopting a single water source withdrawal mode; and put the foregoing modules into the multi-source complementary water supply module for invoking.

Modifying the Rchuse Module and the Res Module

[0132] Add codes separately, and replace parameters waterrch and wuresn in the rchuse module and the res module with parameter WSP.sub.i (i=1,2) respectively to achieve connection of the multi-source complementary water supply module Multi_sc with the rchuse module and the res module well as invoking;

Modifying the Watuse Module

[0133] 1) Modify programs to add functions of transferred water withdrawal and transferred water volume restriction so as to control water supply within a total transferred water limit given that the module does not have the function of transferred water supply, where a calculation formula is expressed as follows:

[00023]${\displaystyle \underset{i}{.Math.}{\displaystyle \underset{j}{.Math.}\text{waterout}\left(i,j\right)\le MX5}}$

where waterout (i,j) denotes transferred water consumption (m.sup.3) of the jth HRU on the ith day; and WX5 denotes total transferred water limit (m.sup.3).

[0134] 2) Add codes in the watuse module, and replacing parameters watershal, waterdeep, waterout, and waterpnd in the watuse module with parameter WSP.sub.i (i=3,4,5,6), respectively to achieve connection of the multi-source complementary water supply module Multi_sc with the watuse module as well as invoking.

Adding a Function of Simulating Water Delivery of a Pipe Network

[0135] In order to resolve loss in urban water supply caused when water “escapes, overflows, drips or leaks” due to aging or breakage of a pipe network, add a calculation program code of the following formula in the rchuse module, the res module, and the watuse module:

[00024]$WSP=WSP\cdot \left(1-pip\right)$

where, pip denotes a leakage rate of water supply pipe network.

Modifying the Irrsub Module

[0136] In terms of the absence of a pit-pond irrigation function in the module, add the pit-pond irrigation function; completing a transferred water irrigation function; and impose water supply restriction to control an irrigation water withdrawal within the total transferred water limit.

[00025]${\displaystyle \underset{i}{.Math.}{\displaystyle \underset{j}{.Math.}\text{wirrout}\left(i,j\right)\le MX5}}$

[00026]${\displaystyle \underset{i}{.Math.}{\displaystyle \underset{j}{.Math.}\text{wirrpnt}\left(i,j\right)\le MX6}}$

where wirrout (i,j) denotes transferred water irrigation consumption (m.sup.3) of the jth HRU on the ith day; and WX5 denotes total transferred water limit (m.sup.3). wirrpnt (i,j) denotes pit-pond irrigation consumption (m.sup.3) of the jth HRU on the ith day; and WX6 denotes pit-pond available water supply (m.sup.3).

Adding a Function of Simulating an Irrigation Channel

[0137] 1) Given that the SWAT model does not consider the influence of water delivery loss of an irrigation channel on a hydrological cycle, irrigation-related leakage is deemed as system loss, add modify the irr_rch module, the irr_res module, and the irrsub module to add simulation on a channel system delivery process comprising channel water loss and channel recession, where the channel water loss includes two parts of channel water evaporation loss and channel leakage loss, and the main calculation formulas are as follows:

[00027]$E{T}_{can}=IR{R}_{can}\cdot \left(1-\phi \right)\cdot \alpha$

[00028]$L{s}_{can}=IR{R}_{can}\cdot \left(1-\phi \right)\cdot \beta$

[00029]$Surp{}_{can}=IR{R}_{can}\cdot \left(1-\phi \right)\cdot \left(1-\alpha -\beta \right)$

where, ET.sub.can denotes a channel system evaporation loss (mm); IRR.sub.can denotes an irrigation water volume (mm) entering a channel; Ls.sub.can denotes a channel system leakage loss (mm); Surp.sub.can denotes a channel system recession volume (mm); φ denotes an effective utilization coefficient of channel system water; α denotes a channel system evaporation coefficient; and β denotes a channel system leakage coefficient.

[0138] 2) Drive water from leakage loss to enter upper soil as replenished soil water and participate in the soil moisture cycle, and add a calculation program for leakage loss by modifying relevant codes of a percmain module, and the calculation formula is as follows:

[00030]$Ws{l}_{yr1,t}=Ws{l}_{yr1,t-1}+{\inf }_{pcp}+{\inf }_{irr}+{\inf }_{wet}+L{s}_{can}$

where,

[00031]$Ws{l}_{yr1,t+1}$

denotes soil water content (mm) of a first layer of soil on the t-th day;

[00032]$Ws{l}_{yr1,t}$

denotes soil water content (mm) of a first layer of soil on the (t-1)th day; inf.sub.pcp denotes precipitation infiltration capacity (mm); inf.sub.irr denotes irrigation infiltration capacity (mm); and inf.sub.wet denotes lake and reservoir wetland infiltration capacity (mm);

Modifying Codes of the Gwmod Module

[0139] Drive water loss from leakage of a water supply pipe network to enter an underground aquifer as replenished groundwater, and modify groundwater recharge codes in a gwmod module to achieve simulation on water leakage of a pipe network, and the calculation formula is as follows:

[00033]$\begin{array}{c}rh{}_{\epsilon }=\left(1-\exp \left(-1/GW\_DELAY\right)\right)\cdot \left(prc+WSP\cdot pip/Area\right)\\ +\exp \left(-1/GW\_DELAY\right)\cdot r{h}_{L-1}\end{array}$

where,

[00034]$r{h}_t$

denotes groundwater recharge capacity (mm) on the t-th day;

[00035]$r{h}_{t-1}$

denotes groundwater recharge capacity (mm) on the (t-1)th day; prc denotes soil water leakage (mm) of recharged groundwater; GW_DELAY denotes groundwater recharge delay coefficient (mm); and Area denotes the area (m.sup.2) of an HRU;

Modifying the Subbasin Module

[0140] Add the multi-source complementary water supply module in a subbasin module to facilitate revoking during the running of the SWAT model, and conduct in-year dynamic complementary water supply operation on water sources by reading specified type, number, water source codes, water withdrawal volume, and water withdrawal time of water sources to achieve multi-source combined water supply simulation during the running of the SWAT model.

Modifying the Surface Module

[0141] Channel system recession refers to irrigation water discharged and overdrawn from the channel, which directly enters the river channel and participates in the calculation of flow concentration of river channel, and superimpose the earth surface runoff by modifying the relevant codes of the surface module, and the calculation formula is as follows:

[00036]$sur{f}_t=sur{f}_0+Sur{p}_{can}$

where,

[00037]$sur{f}_t$

denotes runoff (mm) after channel recession; and

[00038]$sur{f}_0$

denotes runoff (mm) before channel recession.

Modifying the Point Source Module

[0142] After urban domestic and industrial sewage is produced, it is directly discharged into the river through the drainage pipe network system, or transported to the sewage disposal plant for disposal. Some of the up-to-standard sewage obtained after disposal is directly discharged into the river, and some of the up-to-standard reclaimed water obtained after advanced disposal is reused for greening, domestic miscellaneous use, production and other purposes. The SWAT model conducts simulation using a point source module (recday module or recmon module). Modify relevant codes in the recday module and recmon module, and use a pollution discharge parameter WDR to replace parameters floday and flomon, respectively, and the calculation formulas are as follows:

[00039]$WP=WSP\cdot \left(1-r\right)$

[00040]$WDR=WP\cdot \left(1-v\right)+WP\cdot v\cdot \left(1-re\right)$

where WDR denotes urban sewage output (m.sup.3); WP denotes sewage discharge (m.sup.3); r denotes a water consumption rate; v denotes a sewage disposal rate of a sewage disposal plant; and re denotes a reclaimed water utilization rate.

Modifying the Main Module

[0143] Put the constructed readattr module in the main module in the SWAT model to facilitate invocation of the readattr module.

[0144] Modify SWAT source codes in the Windows platform application development environment Visual Studio 2012 using Fortran language based on the foregoing structures, which achieves revocation of a multi-source complementary water supply module with the functions of urban multi-source supply and multi-source agricultural irrigation. Therefore, a SWAT model based distributed model for simulating bidirectional reciprocation of natural-artificial hydrological cycles with the function of multi-source complementation is constructed.

[0145] The present disclosure selects Baihe Basin as an embodiment, and the simulation process of the natural-artificial hydrological cycle in Baihe Basin is described as follows:

1. Overview of Study Area

[0146] Suited in the middle part of the Hanjiang River basin, Baihe Basin runs from Funiu Mountain in the north to the Hanjiang River in the south, bordering Laoguan River in the west and Tang River Basin in the east. With most of regions located in the Nanyang Basin, Baihe Basin covers a total area of 12,300 km.sup.2. Baihe Basin exhibits a tendency of elevating from south to north. Its river systems mainly cover Baihe and its tributaries including Tuan River, the turbulent River, Diao River and Yanling River. The mean annual gross amount of water resources reaches 2.08 billion m.sup.3. Within the basin, there are 3 prefecture-level cities including Nanyang City, Henan Province, as well as Xiangyang City, Hubei Province, and a total of 12 county-level administrative units including Wolong District and Wancheng District, Xinye County and Xiangzhou District. Baihe Basin covers 4 large-sized irrigation areas such as Yahekou irrigation area and Yindan irrigation area, and 18 medium-sized irrigation areas such as Zhaowan reservoir irrigation area, Gaoqiu reservoir irrigation area and Doupo reservoir irrigation area. In order to achieve flood control and drought relief and ensure accessibility of agricultural irrigation water, there are one large-sized reservoir (Yahekou Reservoir) and more than a dozen medium-sized reservoirs such as Zhaowan Reservoir, Hushan Reservoir, Doupo Reservoir and Guangou Reservoir.

2. Basic Data Collection

[0147] Data required for the construction of the model include spatial data, covering DEM (90 m×90 m), land use map (1:100000), soil distribution map (1:1000000), administrative region distribution map, irrigation area distribution map, and water system map, etc; meteorological data from 1990 through 2016 of four meteorological stations in and adjoining the Baihe Basin, including daily precipitation, daily maximum and minimum temperature, relative humidity, sunshine duration, wind speed, and other meteorological elements; water conservancy project data, mainly including reservoir location, inactive reservoir capacity, regulated reservoir capacity, total reservoir capacity, discharge capacity of irrigation channels, and daily water pumping capacity of motor-driven wells; monthly runoff data of Xindianpu Hydrological Station from 1991 through 2016 used for model calibration and validation; data about agricultural irrigation water over the years, which is obtainable by querying Nanyang City Water Resources Bulletin and Xiangyang City Water Resources Bulletin during 2006 - 2016; and data about irrigation areas and plantation structures, including plantation areas of wheat, rice, peanuts, rape, sesame, cotton, vegetables, melons and fruits and other crops. In addition, it is necessary to collect basic information about the time, frequency and quantity of single (irrigation, fertilization) related to crop management measures such as sowing, irrigation, fertilization and harvesting.

3. Modeling Process

[0148] In combination with the water resources distribution, characteristics of water conservancy project, and water withdrawal and utilization of various industries in the Baihe Basin, model building includes the following steps:

[0149] (1) HRU division: obtain, by division, a total of 34 natural sub-basins and 507 natural HRUs according to DEM data, land-use type map and soil type map; and on this basis, subdivide natural HRUs into 1,027 HRUs according to administrative division and distribution of irrigation areas, the HRUs each having a sub-basin attribute, a water resources region attribute, an administrative region attribute, an irrigation area attribute, a land use type attribute and a soil type attribute.

[0150] (2) Input of information about agricultural plantation and management, including crop types, crop planting area and irrigation area, crop rotation system and irrigation system in administrative regions and counties.

[0151] (3) Input of information about water supply rules, mainly including information such as water supply objects, regional water supply principle, water supply priority, water distribution principle of water users, type and number of water sources, and industry water use priority; in addition, the information that needs to be input also includes information about water conservancy projects such as reservoirs, channels, motor-driven wells, effective utilization coefficient of channel system water, effective utilization coefficient of field water, water consumption rate and so on.

[0152] (4) Input of meteorological data information: input data about established precipitation, air temperature, wind speed, radiation and relative humidity, afterwards, input all attribute data and reservoir data, and then start running after the construction of the model is completed.

4. Parameter Calibration and Model Validation

[0153] Through the above process and analysis, the main parameter values of the model are determined, and the final values after parameter adjustment are shown in Table 7. Through adjustment in the parameters of the model, the results of comparison between the simulated and measured runoff of the model are shown in Table 8 and FIG. 5. In the calibration period (1995-2005), the correlation coefficient between monthly runoff simulation value and measured value in Xindianpu station is 0.792, and the Nash-Sutcliffe efficiency coefficient is 0.756; in the validation period (2006-2016), the correlation coefficient between monthly runoff simulation value and measured value in Xindianpu station is 0.643, and the Nash-Sutcliffe efficiency coefficient is 0.635. The results show that the fitting degree of the monthly runoff simulated value and measured value is high, and the simulation accuracy of the model reaches the required value.

[0154] The deviation percentages of the simulation results of all administrative regions (Wolong District, Wancheng District, and Zhenping County) located in the basin are shown in Table 9. As can be seen from this table, the water supply volumes of the three counties (districts) from 2010 through 2016 have seen a slight deviation, in which deviations of domestic water consumption, industrial water consumption, agricultural water consumption, surface water supply and groundwater supply from actual water volumes are all within 10%. As can be seen, simulation results of the model well reflect actual water consumption and supply in the Baihe Basin.

5. Analysis of Simulation Results

[0155] Through the simulation and summary of hydrological cycle in Baihe Basin in 2016, the relationship of cycle transformation in water resources in the Basin is shown in FIG. 6. In 2016, the whole basin saw 9,547 million m.sup.3 of total precipitation, 2,273 million m.sup.3 of earth surface runoff, 862 million m.sup.3 of year-end soil water storage variable, 125 million m.sup.3 of surface water storage variable, and 65 million m.sup.3 of groundwater storage variable. The whole basin saw 7,061 million m.sup.3 of total water consumption, including 2,990 million m.sup.3 of soil evaporation, 3,309 million m.sup.3 of vegetation transpiration, 625 million m.sup.3 of interception evaporation, 0.3 million m.sup.3 of snow sublimation, 117 million m.sup.3 of water surface evaporation, and 45 million m.sup.3 of domestic and industrial consumption. When the consumption of transferred water is considered, the whole basin saw 953 million m.sup.3 of total economic and social water consumption, including 621 million m.sup.3 of surface water consumption (including 447 million m.sup.3 of transferred water) and 332 million m.sup.3 of groundwater consumption; and moreover, the whole basin saw 534 million m3 of total artificial water consumption, 171 million m.sup.3 of artificial displacement, and 2,030 million m.sup.3 of total 9.

[0156] The water supply of water sources in the Baihe Basin in 2016 is summarized in FIG. 7. In 2016, water sources in Baihe Basin mainly centered around groundwater and transferred water (supplied by Danjiangkou Reservoir). Throughout the year, the development and consumption of groundwater accounted for 34.83% (332 million m.sup.3) of the total water consumption in the Basin, the consumption of transferred water accounted for 46.92% (447 million m.sup.3) of the total water consumption in the Basin, and the consumption of urban river water came last, which only accounted for 1.95% (18 million m.sup.3) of the total water consumption in the Basin. Throughout the year, peak of water consumption appeared in March and August. March is the key period of water consumption for winter wheat, and during this period, the consumption of transferred water is 86 million m.sup.3, followed by water storage and water supply (42 million m.sup.3); given that the surface water resources are relatively abundant, groundwater is only used for supplementary water supply, and the supplementary water supply is 26 million m.sup.3. August is the key period of water consumption for maize growth, and during this period, the consumption of transferred water is 77 million m.sup.3, water storage resources are inadequate and thus can only supply 20 million m.sup.3 of water, while groundwater plays an important role in supplementary water supply, and can supply 75 million m.sup.3 of water, which is almost equal to the supply of transferred water. In addition, the small peak of water consumption appeared from October to November and in January, mainly due to the irrigation water after winter wheat was planted. During this period, there are limited surface water resources, water is mainly supplied from Danjiangkou Reservoir (transferred water) and groundwater for supplementary irrigation.

TABLE-US-00001 HRU code Sub-basin attribute Water resources region attribute Administrative region attribute Irrigation area attribute Field format Shaping Shaping Shaping Shaping Note Code of sub-basin Code of water resources region Code of administrative region Code of irrigation area Content Code of sub-basin to which an HRU belongs Code of water resources region to which an HRU belongs Code of administrative region to which an HRU belongs Code of irrigation area to which an HRU belongs

TABLE-US-00002 Parameter name Data type Note Content cntyID Shaping Code of administrative region Sequence number of administrative region Riv1 Shaping No. 1 urban river water source Code of sub-basin where a river channel lies Riv2 Shaping No. 2 urban river water source Code of sub-basin where a river channel lies Riv3 Shaping No. 3 urban river water source Code of sub-basin where a river channel lies Riv4 Shaping No. 4 urban river water source Code of sub-basin where a river channel lies Riv5 Shaping No. 5 urban river water source Code of sub-basin where a river channel lies Res1 Shaping No. 1 reservoir water source Code of reservoir Res2 Shaping No. 2 reservoir water source Code of reservoir Res3 Shaping No. 3 reservoir water source Code of reservoir Res4 Shaping No. 4 reservoir water source Code of reservoir Res5 Shaping No. 5 reservoir water source Code of reservoir Shall Shaping No. 1 shallow water source Code of sub-basin where a shallow water layer lies Shal2 Shaping No. 2 shallow water source Code of sub-basin where a shallow water layer lies Shal3 Shaping No. 3 shallow water source Code of sub-basin where a shallow water layer lies Shal4 Shaping No. 4 shallow water source Code of sub-basin where a shallow water layer lies Shal5 Shaping No. 5 shallow water source Code of sub-basin where a shallow water layer lies Deep1 Shaping No. 1 deep water source Code of sub-basin where a deep water layer lies Deep2 Shaping No. 2 deep water source Code of sub-basin where a deep water layer lies Deep3 Shaping No. 3 deep water source Code of sub-basin where a deep water layer lies Deep4 Shaping No. 4 deep water source Code of sub-basin where a deep water layer lies Deep5 Shaping No. 5 deep water source Code of sub-basin where a deep water layer lies Out1 Shaping No. 1 pit-pond water source Code of sub-basin where a pit-pond lies Out2 Shaping No. 2 pit-pond water source Code of sub-basin where a pit-pond lies Out3 Shaping No. 3 pit-pond water source Code of sub-basin where a pit-pond lies Out4 Shaping No. 4 pit-pond water source Code of sub-basin where a pit-pond lies Out5 Shaping No. 5 pit-pond water source Code of sub-basin where a pit-pond lies Pnd1 Shaping No. 1 external-basin water source Code of external water Pnd2 Shaping No. 2 external-basin water source Code of external water Pnd3 Shaping No. 3 external-basin water source Code of external water Pnd4 Shaping No. 4 external-basin water source Code of external water Pnd5 Shaping No. 5 external-basin water source Code of external water

[0157] If there is no water supply in Table 2, the code is represented by 0.

TABLE-US-00003 Parameter name Data type Note Content irrID Shaping Code of irrigation area Sequence number of irrigation area Riv1 Shaping No. 1 urban river water source Code of sub-basin where a river channel lies Riv2 Shaping No. 2 urban river water source Code of sub-basin where a river channel lies Riv3 Shaping No. 3 urban river water source Code of sub-basin where a river channel lies Riv4 Shaping No. 4 urban river water source Code of sub-basin where a river channel lies Riv5 Shaping No. 5 urban river water source Code of sub-basin where a river channel lies Res1 Shaping No. 1 reservoir water source Code of reservoir Res2 Shaping No. 2 reservoir water source Code of reservoir Res3 Shaping No. 3 reservoir water source Code of reservoir Res4 Shaping No. 4 reservoir water source Code of reservoir Res5 Shaping No. 5 reservoir water source Code of reservoir Shal1 Shaping No. 1 shallow water source Code of sub-basin where a shallow water layer lies Shal2 Shaping No. 2 shallow water source Code of sub-basin where a shallow water layer lies Shal3 Shaping No. 3 shallow water source Code of sub-basin where a shallow water layer lies Shal4 Shaping No. 4 shallow water source Code of sub-basin where a shallow water layer lies Shal5 Shaping No. 5 shallow water source Code of sub-basin where a shallow water layer lies Deep1 Shaping No. 1 deep water source Code of sub-basin where a deep water layer lies Deep2 Shaping No. 2 deep water source Code of sub-basin where a deep water layer lies Deep3 Shaping No. 3 deep water source Code of sub-basin where a deep water layer lies Deep4 Shaping No. 4 deep water source Code of sub-basin where a deep water layer lies Deep5 Shaping No. 5 deep water source Code of sub-basin where a deep water layer lies Out1 Shaping No. 1 pit-pond water source Code of sub-basin where a pit-pond lies Out2 Shaping No. 2 pit-pond water source Code of sub-basin where a pit-pond lies Out3 Shaping No. 3 pit-pond water source Code of sub-basin where a pit-pond lies Out4 Shaping No. 4 pit-pond water source Code of sub-basin where a pit-pond lies Out5 Shaping No. 5 pit-pond water source Code of sub-basin where a pit-pond lies Pnd1 Shaping No. 1 external-basin water source Code of external water Pnd2 Shaping No. 2 external-basin water source Code of external water Pnd3 Shaping No. 3 external-basin water source Code of external water Pnd4 Shaping No. 4 external-basin water source Code of external water Pnd5 Shaping No. 5 external-basin water source Code of external water

[0158] If there is no water supply in Table 3, the code is represented by 0.

TABLE-US-00004 Parameter name Data type Note Content cntyID Shaping Code of administrative region Sequence number of administrative region sup_Riv1 Shaping Water withdrawal sequence number of No. 1 urban river water source Code of water supply sequence number sup_Riv2 Shaping Water withdrawal sequence number of No. 2 urban river water source Code of water supply sequence number sup_Riv3 Shaping Water withdrawal sequence number of No. 3 urban river water source Code of water supply sequence number sup_Riv4 Shaping Water withdrawal sequence number of No. 4 urban river water source Code of water supply sequence number sup_Riv5 Shaping Water withdrawal sequence number of No. 5 urban river water source Code of water supply sequence number sup_Res1 Shaping Water withdrawal sequence number of No. 1 reservoir water source Code of water supply sequence number sup_Res2 Shaping Water withdrawal sequence number of No. 2 reservoir water source Code of water supply sequence number sup_Res3 Shaping Water withdrawal sequence number of No. 3 reservoir water source Code of water supply sequence number sup_Res4 Shaping Water withdrawal sequence number of No. 4 reservoir water source Code of water supply sequence number sup_Res5 Shaping Water withdrawal sequence number of No. 5 reservoir water source Code of water supply sequence number sup_Shal1 Shaping Water withdrawal sequence number of No. 1 shallow water source Code of water supply sequence number sup_Shal2 Shaping Water withdrawal sequence number of No. 2 shallow water source Code of water supply sequence number sup_Shal3 Shaping Water withdrawal sequence number of No. 3 shallow water source Code of water supply sequence number sup_Shal4 Shaping Water withdrawal sequence number of No. 4 shallow water source Code of water supply sequence number sup_Shal5 Shaping Water withdrawal sequence number of No. 5 shallow water source Code of water supply sequence number sup_Deep1 Shaping Water withdrawal sequence number of No. 1 deep water source Code of water supply sequence number sup_Deep2 Shaping Water withdrawal sequence number of No. 2 deep water source Code of water supply sequence number sup_Deep3 Shaping Water withdrawal sequence number of No. 3 deep water source Code of water supply sequence number sup_Deep4 Shaping Water withdrawal sequence number of No. 4 deep water source Code of water supply sequence number sup_Deep5 Shaping Water withdrawal sequence number of No. 5 deep water source Code of water supply sequence number sup_Out1 Shaping Water withdrawal sequence number of No. 1 pit-pond water source Code of water supply sequence number sup_Out2 Shaping Water withdrawal sequence number of No. 2 pit-pond water source Code of water supply sequence number sup_Out3 Shaping Water withdrawal sequence number of No. 3 pit-pond water source Code of water supply sequence number sup_Out4 Shaping Water withdrawal sequence number of No. 4 pit-pond water source Code of water supply sequence number sup_Out5 Shaping Water withdrawal sequence number of No. 5 pit-pond water source Code of water supply sequence number sup_Pnd1 Shaping Water withdrawal sequence number of No. 1 external-basin water source Code of water supply sequence number sup_Pnd2 Shaping Water withdrawal sequence number of No. 2 external-basin water source Code of water supply sequence number sup_Pnd3 Shaping Water withdrawal sequence number of No. 3 external-basin water source Code of water supply sequence number sup_Pnd4 Shaping Water withdrawal sequence number of No. 4 external-basin water source Code of water supply sequence number sup_Pnd5 Shaping Water withdrawal sequence number of No. 5 external-basin water source Code of water supply sequence number

[0159] If there is no water supply in Table 4, the code is represented by 0.

TABLE-US-00005 Parameter name Data type Note Content irrID Shaping Code of irrigation area Sequence number of irrigation area irr_Riv1 Shaping Water withdrawal sequence number of No. 1 urban river water source Code of water supply sequence number irr_Riv2 Shaping Water withdrawal sequence number of No. 2 urban river water source Code of water supply sequence number irr_Riv3 Shaping Water withdrawal sequence number of No. 3 urban river water source Code of water supply sequence number irr_Riv4 Shaping Water withdrawal sequence number of No. 4 urban river water source Code of water supply sequence number irr_Riv5 Shaping Water withdrawal sequence number of No. 5 urban river water source Code of water supply sequence number irr_Res1 Shaping water withdrawal sequence number of No. 1 reservoir water source Code of water supply sequence number irr_Res2 Shaping water withdrawal sequence number of No. 2 reservoir water source Code of water supply sequence number irr_Res3 Shaping water withdrawal sequence number of No. 3 reservoir water source Code of water supply sequence number irr_Res4 Shaping water withdrawal sequence number of No. 4 reservoir water source Code of water supply sequence number irr_Res5 Shaping water withdrawal sequence number of No. 5 reservoir water source Code of water supply sequence number irr_Shal1 Shaping water withdrawal sequence number of No. 1 shallow water source Code of water supply sequence number irr_Shal2 Shaping water withdrawal sequence number of No. 2 shallow water source Code of water supply sequence number irr_Shal3 Shaping water withdrawal sequence number of No. 3 shallow water source Code of water supply sequence number irr_Shal4 Shaping water withdrawal sequence number of No. 4 shallow water source Code of water supply sequence number irr_Shal5 Shaping water withdrawal sequence number of No. 5 shallow water source Code of water supply sequence number irr_Deep1 Shaping water withdrawal sequence number of No. 1 deep water source Code of water supply sequence number irr_Deep2 Shaping water withdrawal sequence number of No. 2 deep water source Code of water supply sequence number irr_Deep3 Shaping water withdrawal sequence number of No. 3 deep water source Code of water supply sequence number irr_Deep4 Shaping water withdrawal sequence number of No. 4 deep water source Code of water supply sequence number irr_Deep5 Shaping water withdrawal sequence number of No. 5 deep water source Code of water supply sequence number irr_Out1 Shaping water withdrawal sequence number of No. 1 pit-pond water source Code of water supply sequence number irr_Out2 Shaping water withdrawal sequence number of No. 2 pit-pond water source Code of water supply sequence number irr_Out3 Shaping water withdrawal sequence number of No. 3 pit-pond water source Code of water supply sequence number irr_Out4 Shaping water withdrawal sequence number of No. 4 pit-pond water source Code of water supply sequence number irr_Out5 Shaping water withdrawal sequence number of No. 5 pit-pond water source Code of water supply sequence number irr_Pnt1 Shaping water withdrawal sequence number of No. 1 external-basin water source Code of water supply sequence number irr_Pnt2 Shaping water withdrawal sequence number of No. 2 external-basin water source Code of water supply sequence number irr_Pnt3 Shaping water withdrawal sequence number of No. 3 external-basin water source Code of water supply sequence number irr_Pnt4 Shaping water withdrawal sequence number of No. 4 external-basin water source Code of water supply sequence number irr_Pnt5 Shaping water withdrawal sequence number of No. 5 external-basin water source Code of water supply sequence number

[0160] If there is no water supply in Table 5, the code is represented by 0.

TABLE-US-00006 Parameter name Data type Note Content Cnty_surf(1,1) Example Surface water withdrawal of No. 1 administrative region in the 1st year Water withdrawal Cnty_surf(1,2) Example Surface water withdrawal of No. 1 administrative region in the 2nd year Water withdrawal : : : : Cnty _surf(M,N) Example Surface water withdrawal of No. M administrative region in Nth year Water withdrawal Cnty_gw(1,1) Example Groundwater withdrawal of No. 1 administrative region in 1st year Water withdrawal Cnty_gw(1,2) Example Groundwater withdrawal of No. 1 administrative region in the 2nd year Water withdrawal Cnty_gw(M,N) Example Groundwater withdrawal of No. M administrative region in Nth year Water withdrawal

TABLE-US-00007 Parameter name Way of parameter adjustment Physical significance Parameter adjustment value CN2 r Initial SCS Runoff Curve Number under humid condition II 1.37 GWQMN v Threshold depth of “base flow” produced by shallow aquifer 1098 GW DELAY v Groundwater recharge delay coefficient 10.4 ALPHA BF v Baseflow alpha factor 0.38 ESCO v Soil evaporation compensation coefficient 0.49 EPCO v Plant absorption compensation factor 0.29 GW_REVAP v Shallow groundwater reevaporation coefficient 0.10 REVAPMN v Shallow groundwater reevaporation threshold 522 α v Evaporation ratio of channel system water loss 0.09 SOL K r Saturated hydraulic conductivity in soil 0.52 RCHRG DP v Permeability ratio of deep aquifer 0.15 β v Infiltration ratio of channel system water loss 0.50 SOL AWC r Available soil moisture 1.11 pip v Leakage rate of pipe network 0.10 φ v Effective utilization coefficient of channel system water 0.56~0.69 ω v Effective utilization coefficient of field water 0.95 r v Water consumption rate of residential and industrial land 0.2~0.45 v v Sewage disposal rate 1.0 re v Reclaimed water utilization rate 0

[0161] In Table 7, v indicates that the parameter adjustment value replaces the original parameter value; and r indicates that the original parameter value is multiplied by the parameter adjustment value.

TABLE-US-00008 Hydrological station R.sup.2 E.sub.ns Calibration period Validation period Calibration period Validation period Xindianpu 0.792 0.643 0.756 0.635

TABLE-US-00009 Region Water consumption Water supply Domestic water Industrial water Agricultural water Total water consumption Surface water supply Groundwater supply Wolong District 0.54% 0.72% 3.41% 1.62% 0.01% 1.57% Wancheng District -0.04% -0.02% 1.15% 0.59% 0 1.01% Zhenping County 0.01% 0 0 0 0.01% 0

[0162] Based on the SWAT model, the present disclosure develops a distributed natural-artificial hydrological model cycle. The distributed natural-artificial hydrological cycle model is endowed with the functions of simulating dynamic reciprocation between natural-artificial hydrological cycles, and integrating development, utilization and regulation of water resources, thereby simulating the process of natural-artificial hydrological cycles of a basin based on modes of urban multi-source water supply and multi-source irrigation water supply. During the running of the model, the dynamic reciprocation between natural hydrological cycle and artificial hydrological cycle can be maintained all the time. In this way, not only the influence of the hydrological cycle on artificial water withdrawal and use is reflected, but also the real-time intervention effect on development, utilization and regulation of water resources on the hydrological cycle is reflected. Thus, it provides a scientific reference basis for in-depth understanding of the hydrological cycle mechanism of a basin under the influence of high-intensity human activities, as well as rational development and utilization of water resources.

[0163] The above described are merely specific implementations of the present disclosure, and the protection scope of the present disclosure is not limited thereto. Any modification or replacement easily conceived by those skilled in the art within the technical scope of the present disclosure should fall within the protection scope of the present disclosure. Therefore, the protection scope of the present disclosure should be subject to the protection scope of the claims.