CN113407897B - 基于多源互济供水模式的分布式水循环模型的设计方法 - Google Patents
基于多源互济供水模式的分布式水循环模型的设计方法 Download PDFInfo
- Publication number
- CN113407897B CN113407897B CN202110946366.8A CN202110946366A CN113407897B CN 113407897 B CN113407897 B CN 113407897B CN 202110946366 A CN202110946366 A CN 202110946366A CN 113407897 B CN113407897 B CN 113407897B
- Authority
- CN
- China
- Prior art keywords
- water
- source
- module
- water supply
- hru
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 title claims abstract description 1116
- 238000000034 method Methods 0.000 title claims abstract description 76
- 238000013461 design Methods 0.000 title claims abstract description 17
- 238000004088 simulation Methods 0.000 claims abstract description 57
- 230000008569 process Effects 0.000 claims abstract description 53
- 230000006870 function Effects 0.000 claims abstract description 24
- 239000003621 irrigation water Substances 0.000 claims abstract description 18
- 238000003973 irrigation Methods 0.000 claims description 71
- 230000002262 irrigation Effects 0.000 claims description 65
- 238000004364 calculation method Methods 0.000 claims description 41
- 239000002689 soil Substances 0.000 claims description 38
- 238000012546 transfer Methods 0.000 claims description 30
- 239000003673 groundwater Substances 0.000 claims description 29
- 239000010865 sewage Substances 0.000 claims description 23
- 150000001875 compounds Chemical class 0.000 claims description 18
- 238000001704 evaporation Methods 0.000 claims description 16
- 239000002352 surface water Substances 0.000 claims description 15
- 230000008020 evaporation Effects 0.000 claims description 14
- 238000001764 infiltration Methods 0.000 claims description 11
- 230000008595 infiltration Effects 0.000 claims description 11
- 238000004519 manufacturing process Methods 0.000 claims description 11
- 238000012384 transportation and delivery Methods 0.000 claims description 11
- 238000009826 distribution Methods 0.000 claims description 10
- VHKQAUGSGQUCDU-UHFFFAOYSA-N (6'-methoxy-3-oxospiro[2-benzofuran-1,9'-xanthene]-3'-yl) 2-(pyridin-2-yldisulfanyl)benzoate Chemical group C=1C(OC)=CC=C(C2(C3=CC=CC=C3C(=O)O2)C2=CC=3)C=1OC2=CC=3OC(=O)C1=CC=CC=C1SSC1=CC=CC=N1 VHKQAUGSGQUCDU-UHFFFAOYSA-N 0.000 claims description 9
- 238000010276 construction Methods 0.000 claims description 9
- 239000008400 supply water Substances 0.000 claims description 9
- 230000009189 diving Effects 0.000 claims description 7
- 239000008235 industrial water Substances 0.000 claims description 7
- 239000013589 supplement Substances 0.000 claims description 6
- 238000011438 discrete method Methods 0.000 claims description 4
- 230000003993 interaction Effects 0.000 claims description 4
- 238000005065 mining Methods 0.000 claims description 4
- 101100492584 Caenorhabditis elegans ast-1 gene Proteins 0.000 claims description 3
- 230000000295 complement effect Effects 0.000 claims description 3
- 238000005520 cutting process Methods 0.000 claims description 2
- 230000000694 effects Effects 0.000 abstract description 14
- 238000011161 development Methods 0.000 abstract description 10
- 239000008239 natural water Substances 0.000 abstract description 6
- 230000007246 mechanism Effects 0.000 abstract description 3
- 230000010354 integration Effects 0.000 abstract 1
- 238000007493 shaping process Methods 0.000 description 128
- 108091006146 Channels Proteins 0.000 description 52
- 230000008878 coupling Effects 0.000 description 9
- 238000010168 coupling process Methods 0.000 description 9
- 238000005859 coupling reaction Methods 0.000 description 9
- 230000018109 developmental process Effects 0.000 description 9
- 238000007726 management method Methods 0.000 description 9
- 230000004048 modification Effects 0.000 description 8
- 238000012986 modification Methods 0.000 description 8
- 238000003860 storage Methods 0.000 description 6
- 230000008859 change Effects 0.000 description 5
- 241000209140 Triticum Species 0.000 description 4
- 235000021307 Triticum Nutrition 0.000 description 4
- 238000011160 research Methods 0.000 description 4
- 230000006872 improvement Effects 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 241000272525 Anas platyrhynchos Species 0.000 description 2
- 101000997295 Drosophila melanogaster Potassium voltage-gated channel protein Shal Proteins 0.000 description 2
- 101001135472 Rattus norvegicus Potassium voltage-gated channel subfamily D member 2 Proteins 0.000 description 2
- 230000032683 aging Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 2
- 230000008901 benefit Effects 0.000 description 2
- 230000002457 bidirectional effect Effects 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000003344 environmental pollutant Substances 0.000 description 2
- 230000004720 fertilization Effects 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 238000005192 partition Methods 0.000 description 2
- 231100000719 pollutant Toxicity 0.000 description 2
- 238000001556 precipitation Methods 0.000 description 2
- 238000004804 winding Methods 0.000 description 2
- 102000010637 Aquaporins Human genes 0.000 description 1
- 108010063290 Aquaporins Proteins 0.000 description 1
- 235000017060 Arachis glabrata Nutrition 0.000 description 1
- 244000105624 Arachis hypogaea Species 0.000 description 1
- 235000010777 Arachis hypogaea Nutrition 0.000 description 1
- 235000018262 Arachis monticola Nutrition 0.000 description 1
- 229920000742 Cotton Polymers 0.000 description 1
- 244000241257 Cucumis melo Species 0.000 description 1
- 235000015510 Cucumis melo subsp melo Nutrition 0.000 description 1
- 241000196324 Embryophyta Species 0.000 description 1
- 241000208152 Geranium Species 0.000 description 1
- 240000007594 Oryza sativa Species 0.000 description 1
- 235000007164 Oryza sativa Nutrition 0.000 description 1
- 244000000231 Sesamum indicum Species 0.000 description 1
- 235000003434 Sesamum indicum Nutrition 0.000 description 1
- 229910000831 Steel Inorganic materials 0.000 description 1
- IXKSXJFAGXLQOQ-XISFHERQSA-N WHWLQLKPGQPMY Chemical compound C([C@@H](C(=O)N[C@@H](CC=1C2=CC=CC=C2NC=1)C(=O)N[C@@H](CC(C)C)C(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H](CC(C)C)C(=O)N1CCC[C@H]1C(=O)NCC(=O)N[C@@H](CCC(N)=O)C(=O)N[C@@H](CC(O)=O)C(=O)N1CCC[C@H]1C(=O)N[C@@H](CCSC)C(=O)N[C@@H](CC=1C=CC(O)=CC=1)C(O)=O)NC(=O)[C@@H](N)CC=1C2=CC=CC=C2NC=1)C1=CNC=N1 IXKSXJFAGXLQOQ-XISFHERQSA-N 0.000 description 1
- 240000008042 Zea mays Species 0.000 description 1
- 235000005824 Zea mays ssp. parviglumis Nutrition 0.000 description 1
- 235000002017 Zea mays subsp mays Nutrition 0.000 description 1
- FJJCIZWZNKZHII-UHFFFAOYSA-N [4,6-bis(cyanoamino)-1,3,5-triazin-2-yl]cyanamide Chemical compound N#CNC1=NC(NC#N)=NC(NC#N)=N1 FJJCIZWZNKZHII-UHFFFAOYSA-N 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 239000002131 composite material Substances 0.000 description 1
- 235000005822 corn Nutrition 0.000 description 1
- 238000003967 crop rotation Methods 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000011049 filling Methods 0.000 description 1
- ZZUFCTLCJUWOSV-UHFFFAOYSA-N furosemide Chemical compound C1=C(Cl)C(S(=O)(=O)N)=CC(C(O)=O)=C1NCC1=CC=CO1 ZZUFCTLCJUWOSV-UHFFFAOYSA-N 0.000 description 1
- 238000003306 harvesting Methods 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000002844 melting Methods 0.000 description 1
- 230000008018 melting Effects 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 235000015097 nutrients Nutrition 0.000 description 1
- 235000020232 peanut Nutrition 0.000 description 1
- 230000035515 penetration Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 230000008929 regeneration Effects 0.000 description 1
- 238000011069 regeneration method Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000013468 resource allocation Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 235000009566 rice Nutrition 0.000 description 1
- 229920006395 saturated elastomer Polymers 0.000 description 1
- 238000000638 solvent extraction Methods 0.000 description 1
- 238000009331 sowing Methods 0.000 description 1
- 239000010959 steel Substances 0.000 description 1
- 238000000547 structure data Methods 0.000 description 1
- 238000000859 sublimation Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 235000013311 vegetables Nutrition 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Images
Classifications
-
- G01V20/00—
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/04—Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0631—Resource planning, allocation, distributing or scheduling for enterprises or organisations
- G06Q10/06315—Needs-based resource requirements planning or analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0637—Strategic management or analysis, e.g. setting a goal or target of an organisation; Planning actions based on goals; Analysis or evaluation of effectiveness of goals
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
- G06Q50/06—Electricity, gas or water supply
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/60—Analysis
- G01V2210/66—Subsurface modeling
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F8/00—Arrangements for software engineering
- G06F8/20—Software design
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Economics (AREA)
- Physics & Mathematics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- General Business, Economics & Management (AREA)
- Tourism & Hospitality (AREA)
- Marketing (AREA)
- Operations Research (AREA)
- Game Theory and Decision Science (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Mathematical Physics (AREA)
- Public Health (AREA)
- Pure & Applied Mathematics (AREA)
- Water Supply & Treatment (AREA)
- Computational Mathematics (AREA)
- Primary Health Care (AREA)
- General Health & Medical Sciences (AREA)
- Mathematical Analysis (AREA)
- Data Mining & Analysis (AREA)
- Mathematical Optimization (AREA)
- Software Systems (AREA)
- Algebra (AREA)
- Databases & Information Systems (AREA)
- General Engineering & Computer Science (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geophysics (AREA)
Abstract
本发明提供一种基于多源互济供水模式的分布式水循环模型的设计方法，包括：步骤S1，嵌套式HRU划分；步骤S2，HRU属性设计；步骤S3，多源互济供水模块设计；步骤S4，SWAT模型改进。本发明以SWAT模型为基础，开发一种分布式的自然‑人工水循环动态互馈模拟模型，使其具有了自然水循环与人工侧枝水循环动态互馈、水资源开发利用与调控一体化模拟功能，可完成基于城市多水源供水和多水源灌溉供水模式的流域自然‑人工水循环过程模拟，不仅反映水循环过程对人工取用水的影响，同时还反映水资源开发利用与调控对水循环过程的实时干预作用，从而为深入了解强人类活动影响下的流域水循环机制，合理开发和利用水资源提供科学的参考依据。
Description
技术领域
本发明涉及水文模型的设计技术领域，具体是一种基于多源互济供水模式的分布式水循环模型的设计方法。
背景技术
随着人类生产力的提高和技术的进步，人类活动对水循环的影响明显加大，流域水循环演变呈现出复杂的自然-人工复合特性。传统的水文模型多是以产汇流机制为核心开发的，一般以地表径流模拟为主，对资源开发利用的处理较为粗略，且不能对人工取用水活动及调控决策进行考虑，难以适应这种基于强人类扰动过程的精细化模拟需求。众多来自水文及水文相关学科的专家学者开始尝试在传统的水文模型框架上寻求新的突破，研究人类扰动下的水循环演变规律，分布式自然-人工水循环模型应运而生。
将分布式水文模型与集总式水资源配置模型耦合是构建自然-人工分布式水循环模型的最常用方式之一。分布式水文模型以精确模拟水循环过程为目标，能够补充水资源配置模型无法给出的水循环过程，提供其所需要的实时水资源边界情况；水资源配置模型能够处理人类调控下的水资源开发利用过程，给出取、用、耗、排水数据，有效提高水文模型的模拟精度。二者的耦合可以取长补短，充分发挥各自优势，实现自然-人工水循环精确模拟。如文献一（赵勇.广义水资源合理配置研究[D]. 中国水利水电科学研究院, 2006.）中，将平原区分布式水循环模型(PDWCM)与广义水资源合理配置模型(GWRAM)耦合，采用分解和聚合的方式进行信息交互，实现了自然-人工水循环全过程模拟；在文献二（张洪刚,熊莹,邴建平,李明新.NAM模型与水资源配置模型耦合研究[J].人民长江,2008(17):15-17.）中，利用水文模型NAM与水资源配置模型的耦合，研究强人类活动影响下的汉江流域水资源响应；在文献三（张守平. 基于流域水循环及其伴生过程的水量水质联合配置研究[D].中国水利水电科学研究院,2015.）中，采用集总式水资源配置模型模拟人工“取水-用水-耗水-排水”过程，将以三级区套地市为单位的水资源配置分区数据传输到以子流域高程带划分的计算单元上，实现了与分布式水循环模型WEP-L的耦合，进行了渭河流域自然-人工水循环及伴生过程模拟。
但是，这些模型的构建多以分布式水文模型与集总式水资源配置模型松散耦合的方式进行，也即是采用公用参数或者输出文件单向传递的方式实现模型的耦合，忽略自然水循环与社会水循环过程之间的动态互馈作用，对人类活动用水过程与天然水文过程之间的相互影响反映不足，不利于水资源的精细化模拟与调控。
发明内容
针对当前分布式自然-人工水循环模型在自然-人工水循环耦合模拟及水资源开发利用与调控方面的不足，本发明基于自然-人工水循环理论，以SWAT模型为工具，开发一种分布式自然-人工水循环动态互馈模拟模型，通过添加多源互济供水模块，使SWAT模型具有了自然-人工水循环之间的动态互馈模拟、水资源开发利用与调控等功能，在模型运行中，时刻保持自然水循环与人工侧枝水循环之间的动态互馈关系，不仅反映水循环过程对人工取用水的影响，同时还反映经济社会活动对水循环过程的实时干预作用。
本发明采用如下技术方案实现：
一种基于多源互济供水模式的分布式水循环模型的设计方法，包括如下步骤：
步骤S1：嵌套式HRU划分：采用基于“流域-水资源区-行政区-灌区-土地利用-土壤-坡度”多元属性的嵌套式坡面离散方法进行HRU划分，划分后的HRU具有相应的属性；
步骤S2：HRU属性设计：构建HRU属性识别模块，所述HRU属性识别模块用于识别HRU的属性；
步骤S3：设计多源互济供水模块，所述多源互济供水模块用于调用HRU属性识别模块，以识别HRU的属性，根据识别的属性判定土地利用类型、对应的供水源及供水源的供水优先次序，根据供水源的供水优先次序调用相应的水源模块进行取水；
步骤S4：SWAT模型改进：将多源互济供水模块与SWAT模型中修改的模块进行衔接以实现数据实时交互，所述HRU根据输入的需水信息、供水水源类型信息、供用水优先序规则信息、水利工程信息和由SWAT模型提供的天然水文条件信息进行逐日水资源分配与调控，输出逐日“供水-用水-耗水-排水”的人工侧枝水循环信息，将人工侧枝水循环信息传递给SWAT模型，SWAT模型根据人工侧枝水循环信息模拟刻画天然水循环过程，实时模拟水文变化对水资源开发利用的影响，以及模拟人工用水对下一阶段的水资源和供水的影响，为人工侧枝水循环提供实时天然水文边界条件信息，实现“自然-人工”水循环动态互馈式模拟。
进一步的，所述步骤S1包括：
（1）采用ArcGIS对DEM提取流域河网，划分天然子流域；
（2）在天然子流域上叠加土地利用信息、土壤类型信息和坡度信息，划分出天然HRU；
（3）将天然HRU套上行政区、水资源区边界，将天然HRU进一步切割；
（4）根据灌区的分布，将灌区与天然HRU进行叠加，最终完成HRU划分。这样，每一个HRU都具有子流域属性、水资源区属性、行政区属性、灌区属性、土地利用类型属性和土壤类型属性。
进一步的，所述步骤S2包括：
（1）HRU属性识别模块构建
构建HRU属性识别模块，用于读取指定的HRU的属性，所述指定的HRU的属性包括子流域属性、水资源区属性、行政区属性、灌区属性；
（2）HRU属性识别模块的调用
将构建的HRU属性识别模块放在SWAT模型中的主模块中，便于对HRU属性识别模块的调用。
进一步的，所述步骤S3包括：
（1）水源代码信息文件设计：
所述水源代码信息文件用于读取指定的供水源信息，设置的供水源类型包括外调水、水库水、河道水、浅层地下水、深层地下水、坑塘水6种水源类型，水源代码信息文件通过程序指令读入；
（2）供水优先次序信息文件设计：
所述供水优先次序信息文件用于读取水源供水优先次序信息，指定供水源供水优先次序，供水优先次序信息文件通过程序指令读入；
（3）取水控制量信息文件设计：
所述取水控制量信息文件用于读取地表水供水控制总量和地下水开采控制总量信息，识别HRU所属行政区的年度地表水供水控制总量和地下水开采控制总量，用于后续的水源的取水量计算。
（4）多源互济供水计算设计，计算流程如下：
1）HRU类型识别
首先识别HRU的土地利用类型，如果是居工地，则程序进入城市供水计算流程；如果是灌溉用地，则程序进入灌溉用水计算流程；如果是其他土地利用类型，程序结束；
2）供水源及其优先序识别
通过识别HRU的取水源识别代码，调用相应的水源模块，通过水源代码信息文件和供水优先次序文件的读取，识别每一个HRU的供水源个数、类型、取水先后顺序，根据取水源识别代码依次调用各水源模块；
3）多源互济供水计算
根据HRU的供水次序寻求各水源进行取水，直至满足该HRU的当日的生活、工业用水需求、农业灌溉用水需求、或者到最后一个供水源供水完毕为止。
进一步的，所述多源互济供水计算的步骤包括：
①明确目标HRU设定的当日需水量WD；
②明确目标HRU的供水水源个数k、水源代码及其供水优先次序，其中，k≤30；
③按次序调用水源模块，计算水源的取水量，水源模块包括河道取水模块（rchuse模块）、水库取水模块（res模块）、watuse模块（包含地下水、坑塘水和外调水的取水计算），水源取水量取决于当日HRU的需水量和水源的可供水量，可供水量则取决于水源的可利用水量、取水工程（如引水渠道、供水管道、电井等）供水能力和取水控制量，计算公式如下：
式中，i表示HRU序号；j表示水源的供水优先次序号；WSP表示水源的当日实际取水量，m3；WD表示HRU的当日需水量，m3；Wsc表示水源的当日可供水量，m3；WF表示取水工程的供水能力，m3；WA表示水源的当日可利用水量，m3；WMX表示年度取水控制量，m3，WUM表示年用水控制总量，m3；WSM表示年度地表取水控制总量，m3；WGM表示年度地下水开采控制总量，m3；
④对供水优先次序为1的水源，则优先从该水源取水，如果水源的可供水量Wsc1＞WD，那么该水源供水量为WSP1=WD，供水程序结束，HRU的水源总供水量为WSP= WSP1；否则，WSP1=Wsc1，并且HRU的需水量变为Wf =WD-Wsc1，程序将继续寻找下一级供水源；
⑤对供水优先次序为j（j=2,…,k-1；k≤30）的水源，则优先从该水源取水。如果该水源的当日可供水量Wscj＞Wf，那么该水源供水量为WSPj=Wf，程序结束，HRU的水源总供水量为WSP=WSP+WSPj；否则，WSPj=Wscj，并且HRU的需水量变成Wf=Wf-Wscj，程序将继续寻找下一级供水源；
⑥对供水优先次序为k（k≤30）的水源，如果该水源的当日可供水量Wsck＞Wf，那么该水源供水量为WSPk=Wf，程序结束，HRU的水源总供水量为WSP=∑WSP i ；否则，WSPk=Wsck，并且HRU的需水量变成Wf=Wf-Wsck，程序结束。
进一步的，所述步骤S4中修改SWAT模型具体包括
（1）模块屏蔽
屏蔽各水源模块，即rchuse 模块、res模块、watuse模块、irr_rch 模块、irr_res模块和irrsub模块，不再采取单一水源取水模式；将以上模块放在多源互济供水模块中进行调用；
（2）修改rchuse 模块、res模块
分别添加代码，用参数WSPi（i=1,2）分别替换掉rchuse 模块、res模块中的参数waterrch、wuresn，实现多源互济供水模块Multi_sc与rchuse 模块、res模块的衔接和调用；
（3）修改watuse模块
①由于该模块中没有外调水供水功能，通过修改程序增加外调水取水以及外调水量限制约束，供水量不能超过调水限制总量，计算公式表达如下：
其中，waterout (i,j)表示第i天、第j个HRU的调水利用量，m3；WX5表示调水限制总量，m3；
②在watuse模块中添加代码，用参数WSPi（i=3,4,5,6）分别替换掉watuse模块中的参数watershal、waterdeep、waterout和waterpnd，实现多源互济供水模块Multi_sc与watuse模块的衔接与调用；
（4）添加管网输水模拟功能
城市供水管网在向用水户输送水量的过程中，会因管网材质老化或破损因素产生“跑、冒、滴、漏”输水损失，因此有必要考虑城市供水损失，在rchuse 模块、res模块、watuse模块中添加如下公式的计算程序：
式中，pip表示供水管网漏损率；
（5）修改irrsub模块
由于模块中没有坑塘灌溉，添加坑塘灌溉功能；完善外调水灌溉功能，通过增加供水量限制约束，灌溉取水量不能超过调水限制总量：
其中，wirrout (i,j)表示第i天、第j个HRU的调水灌溉利用量，m3；WX5表示调水限制总量，m3，wirrpnt (i,j)表示第i天、第j个HRU的坑塘灌溉利用量，m3；WX6表示坑塘可供水量，m3；
（6）添加灌溉渠道模拟功能
①SWAT模型并不考虑灌溉渠道的输水损失对水循环的影响，灌溉渗漏视为系统损失，通过修改irr_rch模块、irr_res模块、irrsub模块的源代码，增加渠系输水过程模拟，包括渠道水量损失和渠道退水，其中，渠道水量损失包括渠道内的水分蒸发损失和渠道渗漏损失两部分，主要计算公式如下：
式中，ET can 表示渠系蒸发损失量，mm；IRR can 表示进入渠道的灌溉水量，mm；Ls can 表示渠系渗漏损失量，mm；Surp can 表示渠系退水量，mm；
表示渠系水有效利用系数；
表示渠系蒸发系数；
表示渠系渗漏系数；
②渗漏损失水量进入上层土壤补给土壤水，参与土壤水分循环，通过修改percmain模块相关代码，添加渗漏损失的计算程序，计算公式如下：
式中，
表示第t天第一层土壤的土壤含水量，mm；
表示第t-1天第一层土壤的土壤含水量，mm；inf pcp 表示降雨入渗量，mm；inf irr 表示灌溉入渗量，mm；inf wet 表示湖库湿地渗漏量，mm；
（7）修改gwmod模块代码
供水管网漏失水量进入地下含水层，补给潜水，通过修改gwmod模块中的地下水补给代码，实现对管网漏损水量的模拟，计算公式如下：
式中，
表示第t天的潜水补给量，mm；
表示第t-1天的潜水补给量，mm；prc表示补给地下水的土壤渗漏水量，mm；GW_DELAY表示地下水补给延迟系数，mm；Area表示HRU的面积，m2；
（8）subbasin模块修改
在subbasin模块中添加多源互济供水模块，便于SWAT模型运行时的调用，通过读取指定的水源类型、水源个数、水源代码、取水量、取水时间，进行年内的水源动态互补供水运算，实现SWAT模型运行时的多水源联合供水模拟；
（9）surface模块修改
渠系退水是指渠道排泄超引灌溉水，渠道退水直接进入河道，参与河道汇流计算，通过修改surface模块相关代码，将叠加地表产流量，计算公式如下：
式中，
表示渠道退水后的径流量，mm；
表示渠道退水前的径流量，mm；
（10）点源模块修改
城市生活和工业污水产生后，通过排水管网系统直接排入河流，或者输送至污水处理厂进行处理，处理达标后的污水一部分直接排放进入河道，一部分通过深度处理达到再生水利用标准后回用于绿化、生活杂用以及生产；SWAT模型采用其点源模块进行模拟，所述点源模块包括recday模块和recmon模块，在recday模块和recmon模块中修改相关代码，将污染排放参数WDR分别替换参数floday、flomon，相关计算公式如下：
其中，WDR为城市污水产生量，m3；WP为污水排放量，m3；r表示用水消耗率；v表示污水处理厂的污水处理率；re表示再生水利用率；
（11）main主模块修改
将构建的readattr模块放在SWAT模型中的main主模块中，便于对readattr模块的调用。
与现有技术相比，本发明具有如下有益效果：
（1）模型采用嵌套式坡面离散方法，很好地处理了天然HRU边界、行政区边界、灌区边界互不重合的问题，既能体现传统分布式水循环模型单元划分特点，也满足了水资源区管理和行政区管理相结合的需求；
（2）嵌入的多源互济供水模块实现社会水循环的多水源、多行业联合配置，具有刻画多水源（包括河道、水库、地下水、外调水、坑塘水等）、多工程（包括蓄水工程、引水工程、提水工程、调水工程等）的功能，将水资源系统中的各种水源与用水户之间的拓扑关系、水分在降雨-产流过程以及各社会生产部门中的传递关系进行客观、清晰的描述，真实体现人类活动对水循环过程的影响；
（3）通过添加多源互济供水模块和修改SWAT模型的相关模块，构建了分布式自然-人工水循环动态互馈模拟模型，弥补了以往的水循环模型不能很好地反映人类活动用水过程与天然水文过程之间的相互影响的不足，并且模型在进行长序列模拟时可以根据生活生产取用水的需要以及水源的水量情况，灵活取用不同的水源，更加适用于强人类活动地区。该模型具有了自然-人工水循环之间的动态互馈模拟、水资源开发利用与调控等功能，实现了自然水循环与人工侧枝水循环之间的双向耦合模拟，充分反映复杂水循环系统的动态性、互馈性特征，可以作为区域自然-人工水循环互馈模拟以及水资源系统精细化管理的有力支撑工具。
附图说明
图1为本发明实施例基于多源互济供水模式的分布式水循环模型的设计方法的流程图；
图2为本发明实施例HRU划分流程图；
图3为本发明实施例多源互济供水模块的运行流程图；
图4为本发明实施例水源供水优先计算流程图；
图5为本发明实施例新店铺站的实测与模拟月径流过程对比；
图6为2016水平年白河流域自然-人工水循环转化关系图；
图7为2016水平年白河流域各水源年内供水过程图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明的一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本发明保护的范围。
本发明实施例提供一种基于多源互济供水模式的分布式水循环模型的设计方法，主要步骤包括嵌套式HRU划分、HRU属性设计、多源互济供水模块设计和SWAT模型改进四部分内容。多源互济供水模块设计就是设计并添加相关模块，增加SWAT模型的城乡多水源供水和多水源灌溉模拟功能，完善人工侧枝水循环，进而构建自然-人工分布式水循环模型。SWAT模型改进，即是修改SWAT模型相关模块，便于水源类型、水源个数、水源代码、取水优先序信息的读取以及多源互济模块的调用。所涉及的模块主要包括rchuse模块、res模块、watuse模块、irr_rch模块、irr_res模块、irrsub模块、subbasin模块、percmain模块、gwmod模块、surface模块、recday模块、recmon模块、main模块等。
本发明所述方法包括以下步骤：
步骤S1：嵌套式HRU划分；
步骤S2：HRU属性设计；
步骤S3：多源互济供水模块设计；
步骤S4：SWAT模型改进。
在所述的步骤S1中，采取基于“流域-水资源区-行政区-灌区-土地利用-土壤-坡度”多元属性的嵌套式坡面离散方法进行HRU划分。其主要步骤（图2）主要包括：（1）天然子流域划分；（2）天然HRU划分；（3）天然HRU叠加管理单元；（4）天然HRU叠加灌区。
（1）采用ArcGIS中的 Arc Hydro Tools 工具进行天然子流域划分。首先对 DEM底图进行填洼等修正操作，依据设定的生成子流域临界集水面积CSA（Critical SourceArea）阈值要求，确定格网水流流向，识别流域分水线，并依据流域坡度、坡向、坡长等河网特征参数划分天然子流域；
（2）按照子流域内的土地利用类型、土壤类型和坡度类型进行叠加，划分出若干个天然HRU。 HRU是指下垫面特征相对单一和均匀的区域，在这个区域中的下垫面具有相同的水文特性，每个HRU内只有一种植被、一种土壤和一种坡度类型的组合。
（3）采用选取HRU叠加地市/县行政区和水资源区的嵌套方式。采用GIS的“EditorToolbar”功能，将管理单元边界GIS图直接叠加到上一步所划分的天然HRU图上，在天然HRU与管理区划交界地区，依照边界线将天然HRU一分为二，并将划分的HRU附上管理单元属性；
（4）采用GIS的“Editor Toolbar”功能，在天然HRU与灌区交界地区，依照边界线将天然HRU一分为二，并将划分的HRU附上灌区属性，最终完成HRU的划分。最后，每一个HRU都具有子流域属性、水资源区属性、行政区属性、灌区属性、土地利用类型属性和土壤类型属性。
在所述步骤S2中，HRU属性的设计包括如下内容：
（1）HRU属性识别模块readattr构建：
构建HRU属性识别模块readattr，用于读取指定的HRU属性。设置的HRU属性分别是子流域属性、水资源区属性、行政区属性、灌区属性。输入文件格式见表1。
（2）HRU属性识别模块的调用
将构建的readattr模块（HRU属性识别模块）放在SWAT模型中的main主模块中，便于对readattr模块的调用。
在所述步骤S3中，多源互济供水模块Multi_sc的设计包括如下内容：（1）水源代码信息文件设计；（2）供水优先次序信息文件设计；（3）多源互济供水计算设计。
（1）水源代码信息文件设计：
用于读取指定的供水源信息。设置的供水源类型包括外调水、水库水、河道水、浅层地下水、深层地下水、坑塘水6种水源类型，并且设定单个HRU的同类型水源最多为5个，即每个HRU最多可以设置30个供水源。水源代码以输入文件读入的方式进行指定，并针对不同的HRU类型，居工地HRU以行政区为单位进行指定，同一行政区内的HRU供水源一致；灌溉用地HRU以灌区为单位进行指定，同一灌区内的HRU供水源一致。输入文件格式分别见表2和表3。
（2）供水优先次序信息文件设计：
用于读取指定的水源供水优先次序信息。每一个水源可以为多个HRU利用，每一个HRU也可以从多个供水源取水。根据每一个HRU的供水源类型和个数，设定供水次序。指定供水源的供水优先次序，供水优先次序以输入文件读入的方式进行指定，并针对不同的HRU类型，居工地HRU以行政区为单位进行指定，同一行政区内的HRU水源供水顺序一致；灌溉用地HRU以灌区为单位进行指定，同一灌区内的HRU水源供水顺序一致。优先次序以1到k（k为HRU的供水水源个数，k≤30）进行排序，1为最高优先序，对未指定的水源，序号设置为0。输入文件格式见表4和表5。
（3）取水控制量信息文件设计：
用于读取各行政区各年度的供水控制总量信息。每个HRU根据自身的行政区属性，识别所属行政区的年度地表水供水控制总量和地下水开采控制总量，用于后续的水源的取水量计算。输入文件格式见表6。
（4）多源互济供水计算设计：
多源互济供水模块Multi_sc的运行流程（图3）主要包括：1）HRU类型识别；2）供水源及其优先序识别；3）多源互济供水计算。
1）HRU类型识别
首先识别HRU的土地利用类型。如果是建设用地（城镇用地、农村用地），则程序进入城乡供水计算流程；如果是灌溉用地，则程序进入灌溉用水计算流程；如果是其他土地利用类型，则该HRU无供水要求，程序退出。
2）供水源及其优先序识别
多源互济供水模块通过识别HRU的取水源识别代码，调用相应的水源模块，河道、水库、浅层地下水、深层地下水、外调水和坑塘水的取水源识别代码分别是1、2、3、4、5、6。通过水源代码信息文件和供水优先次序文件的读取，识别每一个HRU的供水源个数、类型、取水先后顺序，根据取水源识别代码依次调用各水源模块。
3）多源互济供水计算
多源互济供水模块根据HRU的供水次序寻求各水源进行取水，直至满足该HRU的当日的生活、工业用水需求（建设用地HRU）或者农业灌溉用水需求（农业用地HRU）或者到最后一个供水源供水完毕为止。供水优先计算的步骤（图4）如下：
①明确目标HRU设定的当日需水量WD；
②明确目标HRU的供水水源个数k（k≤30）、水源代码及其供水优先次序；
③按次序调用水源模块，计算水源的取水量。水源模块包括rchuse模块（河道取水模块）、res模块（水库取水模块）、watuse模块（地下水、坑塘水的取水计算）。水源取水量取决于当日HRU的需水量和水源的可供水量。可供水量则取决于水源的可利用水量、取水工程（如引水渠道、供水管道、电井等）供水能力和取水控制量。计算公式如下：
式中，i表示HRU序号；j表示水源的供水优先次序号；WSP表示水源的当日实际取水量，m3；WD表示HRU的当日需水量，m3；Wsc表示水源的当日可供水量，m3；WF表示取水工程的供水能力，m3；WA表示水源的当日可利用水量，m3；WMX表示年度取水控制量，m3。WUM表示年用水控制总量，m3；WSM表示年度地表取水控制总量，m3；WGM表示年度地下水开采控制总量，m3。
④对供水优先次序为1的水源，则优先从该水源取水。如果水源的可供水量Wsc1＞WD，那么该水源供水量为WSP1=WD，供水程序结束，HRU的水源总供水量为WSP= WSP1；否则，WSP1=Wsc1，并且HRU的需水量变为Wf =WD-Wsc1，程序将继续寻找下一级供水源；
⑤对供水优先次序为j（j=2,…,k-1；k≤30）的水源，如该水源的当日可供水量Wscj＞Wf，那么该水源供水量为WSPj=Wf，程序结束，HRU的水源总供水量为WSP=WSP+WSPj；否则，WSPj=Wscj，并且HRU的需水量变成Wf=Wf-Wscj，程序将继续寻找下一级供水源；
⑥对供水优先次序为k（k≤30）的水源，如果该水源的当日可供水量Wsck＞Wf，那么该水源供水量为WSPk=Wf，程序结束，HRU的水源总供水量为WSP=∑WSP i ；否则，WSPk=Wsck，并且HRU的需水量变成Wf=Wf-Wsck，程序结束。
在所述步骤S4中，SWAT模型改进主要是将多源互济供水模块与SWAT模型中修改的模块进行衔接以实现数据实时交互，所述HRU根据输入的需水信息、供水水源类型信息、供用水优先序规则信息、水利工程信息和由SWAT模型提供的天然水文条件信息进行逐日水资源分配与调控，输出逐日“供水-用水-耗水-排水”的人工侧枝水循环信息，将人工侧枝水循环信息传递给SWAT模型。其中修改SWAT模型中的相关模块，便于水源类型、水源个数、水源代码、取水优先序信息的读取以及多源互济模块的调用。修改的模块主要有rchuse模块、res模块、watuse模块、irr_rch模块、irr_res模块、irrsub模块、subbasin模块、surface模块、gwmod模块、recday模块、recmon模块和main模块等。各模块主要功能如下：
rchuse 模块：河道供水模块，从指定的河道中取水，用于生活工业用水。
res模块：水库供水模块，从指定的水库中取水，用于生活工业用水。
watuse模块：地下水、坑塘供水模块，从指定的子流域内的浅层地下含水层、深层地下含水层和坑塘中取水，用于生活工业用水。
irr_rch 模块：河道灌溉模块，从指定的河道中取水，用于农田灌溉用水。
irr_res模块：水库灌溉模块，从指定的水库中取水，用于农田灌溉用水。
irrsub模块：地下水、外调水灌溉模块，从指定的子流域内的浅层地下含水层、深层地下含水层和外源水源中取水，用于农田灌溉。
percmain模块：渗漏损失模块，根据降雨入渗、灌溉入渗、坑塘入渗计算结果，模拟土壤水变化。
gwmod模块：地下水模块，根据地表入渗结果，模拟浅层、深层地下水变化以及地下径流。
subbasin模块：子流域模块，主要包括蒸发模拟、产流模拟、地下水模拟、植被生长模拟、农田管理模拟、营养物质迁移转化模拟等内容，通过对水文、水质过程相关模块的调用实现。
surface模块：地表径流模拟模块，模拟冠层截流、积雪融雪、地表产流等主要地表水文过程。
recday模块：日点源模块，通过输入逐日污染物排放量，模拟逐日点源变化过程。
recmon模块：月点源模拟，通过输入逐月污染物排放量，模拟逐月点源变化过程。
main模块：SWAT模型主模块，主要包括读入输入文件，分配数组大小、初始化参数、水文过程模拟等内容。
修改SWAT模型中的相关模块，具体包括：
（1）模块屏蔽
屏蔽水源模块，即rchuse 模块、res模块、watuse模块、irr_rch 模块、irr_res模块和irrsub模块，不再采取单一水源取水模式。将以上模块放在多源互济供水模块中进行调用。
（2）修改rchuse 模块、res模块
添加代码，将参数WSP i （i=1,2）分别替换掉rchuse 模块、res模块中的参数waterrch、wuresn，实现多源互济供水模块Multi_sc与rchuse 模块、res模块的衔接和调用；
（3）修改watuse模块
①由于该模块中没有外调水供水功能，通过修改程序增加外调水取水以及外调水量限制约束，供水量不能超过调水限制总量，计算公式表达如下：
其中，
waterout (i,j)表示第i天、第j个HRU的调水利用量，m3；WX5表示调水限制总量，m3。
②在watuse模块中添加代码，用参数WSP i （i=3,4,5,6）分别替换掉watuse模块中的参数watershal、waterdeep、waterout和waterpnd，实现多源互济供水模块Multi_sc与watuse模块的衔接与调用。
（4）添加管网输水模拟功能
城市供水管网在向用水户输送水量的过程中，会因管网材质老化或破损等因素产生“跑、冒、滴、漏”等输水损失。因此有必要考虑城市供水损失。在rchuse 模块、res模块、watuse模块中添加如下公式的计算程序代码：
式中，pip表示供水管网漏损率。
（5）修改irrsub模块
由于模块中没有坑塘灌溉，添加坑塘灌溉功能；完善外调水灌溉功能，通过增加供水量限制约束，灌溉取水量不能超过调水限制总量。
其中，wirrout (i,j)表示第i天、第j个HRU的调水灌溉利用量，m3；WX5表示调水限制总量，m3。wirrpnt (i,j)表示第i天、第j个HRU的坑塘灌溉利用量，m3；WX6表示坑塘可供水量，m3。
（6）添加灌溉渠道模拟功能
①SWAT模型并不考虑灌溉渠道的输水损失对水循环的影响，灌溉渗漏视为系统损失。通过修改irr_rch模块、irr_res模块、irrsub模块的源代码，增加渠系输水过程模拟，包括渠道水量损失和渠道退水。其中，渠道水量损失包括渠道内的水分蒸发损失和渠道渗漏损失两部分。主要计算公式如下：
式中，ET can 表示渠系蒸发损失量，mm；IRR can 表示进入渠道的灌溉水量，mm；Ls can 表示渠系渗漏损失量，mm；Surp can 表示渠系退水量，mm；
表示渠系水有效利用系数；
表示渠系蒸发系数；
表示渠系渗漏系数。
②渗漏损失水量进入上层土壤补给土壤水，参与土壤水分循环。通过修改percmain模块相关代码，添加渗漏损失计算程序。计算公式如下：
式中，
表示第t天第一层土壤的土壤含水量，mm；
表示第t-1天第一层土壤的土壤含水量，mm；inf pcp 表示降雨入渗量，mm；inf irr 表示灌溉入渗量，mm；inf wet 表示湖库湿地渗漏量，mm；
（7）修改gwmod模块代码
供水管网漏失水量进入地下含水层，补给潜水。通过修改gwmod模块中的地下水补给代码，实现对管网漏损水量的模拟。计算公式如下：
式中，
表示第t天的潜水补给量，mm；
表示第t-1天的潜水补给量，mm；prc表示补给地下水的土壤渗漏水量，mm；GW_DELAY表示地下水补给延迟系数，mm；Area表示HRU的面积，m2；
（8）subbasin模块修改
在subbasin模块中添加多源互济供水模块，便于SWAT模型运行时的调用。通过读取指定的水源类型、水源个数、水源代码、取水量、取水时间，进行年内的水源动态互补供水运算，实现SWAT模型运行时的多水源联合供水模拟。
（9）surface模块修改
渠系退水是指渠道排泄超引灌溉水，渠道退水直接进入河道，参与河道汇流计算。通过修改surface模块相关代码，将叠加地表产流量。计算公式如下：
式中，
表示渠道退水后的径流量，mm；
表示渠道退水前的径流量，mm；
（10）点源模块修改
城市生活和工业污水产生后，通过排水管网系统直接排入河流，或者输送至污水处理厂进行处理。处理达标后的污水一部分直接排放进入河道，一部分通过深度处理达到再生水利用标准后回用于绿化、生活杂用以及生产等用途。SWAT模型采用其点源模块（recday模块和recmon模块）进行模拟。在recday模块和recmon模块中修改相关代码，将污染排放参数WDR分别替换参数floday、flomon。相关计算公式如下：
其中，WDR为城市污水产生量，m3；WP为污水排放量，m3；r表示用水消耗率；v表示污水处理厂的污水处理率；re表示再生水利用率。
（11）main主模块修改
将构建的readattr模块放在SWAT模型中的main主模块中，便于对readattr模块的调用。
用Fortran语言在Windows平台应用程序开发环境Visual Studio 2012中根据以上结构对SWAT源代码进行修改，实现了对多源互济供水模块的调用，具有了城市多水源供水和多水源农业灌溉模拟功能，从而构建了以SWAT模型为基础的、具有多源互济功能的分布式自然-人工水循环双向互馈模拟模型。
本发明选取白河流域作为实施案例，白河流域自然-人工水循环模拟过程介绍如下：
1、研究区概况
白河流域地处汉江流域中部，北起伏牛山地区，西与老鹳河为邻，东靠唐河流域，南抵汉江，大部分地区位于南阳盆地内，流域总面积为1.23万km2。地形呈北高南低态势，水系主要由白河及其支流湍河、刁河、严陵河等，多年平均水资源量20.8亿m3。流域内有河南省南阳市卧龙区、宛城区、新野县、湖北省襄阳市襄州区等3个地级市、共12个县级行政单元。白河流域境内有鸭河口灌区、引丹灌区等大型灌区4个，赵湾水库灌区、高丘水库灌区、陡坡水库灌区等中型灌区18个。为防洪抗旱、保障农业灌溉用水，流域内建有大型水库（鸭河口水库）1座，赵湾水库、虎山水库、陡坡水库、官沟水库等中型水库十几座。
2、基础数据收集
模型构建所需要的数据有空间数据，包括DEM（90m×90m），土地利用图（1:10万）、土壤分布图（1:100万）、行政区分布图、灌区分布图、以及水系图等；白河流域内及邻近的4个气象站的1990~2016年气象数据，包括日降水量、日最高及最低气温、相对湿度、日照时数风速等气象要素；水利工程数据，主要包括水库位置、死库容、调节库容、总库容等数据，灌溉渠道过水能力、电机井日提水能力等数据；新店铺水文站1991~2016年的逐月径流数据，用于模型的率定和验证；农业灌溉用水数据，通过查询2006~2016年南阳市水资源公报、襄阳市水资源公报等资料，获得历年农业灌溉用水数据；灌溉面积与种植结构数据，包括小麦、水稻、花生、油菜、芝麻、棉花、蔬菜、瓜果等各类作物种植面积。此外，还有播种、灌溉、施肥、收割等作物管理措施有关的时间、次数、单次（灌溉、施肥）数量等基本信息等。
3、建模过程
针对白河流域水资源分布、水利工程特点及各行业取用水情况，模型构建步骤如下：
（1）HRU划分：基于DEM数据、土地利用类型图、土壤类型图划分天然子流域和HRU，共获得34个天然子流域和507个天然HRU；在此基础上，根据行政区划及灌区分布情况，将天然HRU细分成1027个HRU，每个HRU均具有子流域属性、水资源区属性、行政区属性、灌区属性、土地利用类型属性和土壤类型属性。
（2）农业种植与管理信息输入：包括各行政区县的种植作物类型、作物种植面积与灌溉面积、轮作制度与灌溉制度等。
（3）供水规则信息输入：主要包括水源供水对象、分区供水原则、水源供水优先级、用水户分水原则、供水水源类型与个数、行业用水优先级等信息；此外，需要输入的信息还有水库、渠道、电机井等水利工程信息，渠系水有效利用系数、田间水有效利用系数等、耗水率等信息。
（4）气象数据信息输入：选择建立好的降水、气温、风速、辐射和相对湿度数据输入，再把所有属性数据及水库数据输入，模型构建完成后即可开始运行。
4、参数率定与模型校验
通过上述过程和分析，确定模型的主要参数值，参数调整的最终值见表7所示。通过对模型参数的调整，模型的径流模拟与实测径流过程对比结果见表8和图5所示。在率定期（1995~2005年），新店铺站月径流模拟值与实测值的相关系数为0.792，纳什效率系数为0.756；在验证期（2006~2016年），新店铺站月径流模拟值与实测值的相关系数为0.643，纳什效率系数为0.635。可以看出，月径流模拟值与实测值流量过程线拟合程度较好，模型的模拟精度达到了要求值。
对全部位于流域内的行政区（卧龙区、宛城区、镇平县）的模拟结果的偏差百分比见表9所示。可以看出，2010~2016年三个县（区）的供用水量偏差很小，其生活用水量、工业用水量、农业用水量、地表水供水量以及地下水供水量与实际水量的误差均在10%以内。可见，模型的模拟结果很好地反映了白河流域实际用水和供水情况。
5、模拟结果分析
通过对白河流域2016年水循环模拟汇总，获得流域水资源循环转化关系见图6所示。2016年全流域总降水量95.47亿m3，地表产流量22.73亿m3，年末土壤水蓄变量8.62亿m3、地表水蓄变量1.25亿m3、地下水蓄变量为0.65亿m3。全流域耗水总量为70.61亿m3，其中土壤蒸发29.90亿m3，植被散发33.09亿m3，截留蒸发6.25亿m3，积雪升华0.003亿m3，水面蒸发1.17亿m3，生活工业消耗0.45亿m3。考虑外调水利用，全流域的经济社会总用水量为9.53亿m3，其中，地表水利用量6.21亿m3（含外调水4.47亿m3），地下水利用量3.32亿m3；人工总耗水量为5.34亿m3，人工排水量为1.71亿m3，总出境水量20.30亿m3。
汇总2016年白河流域各水源供水情况如图7所示。2016年，白河流域水资源利用以地下水和外调水（丹江口水库供水）为主，全年地下水开发利用量占流域用水总量的34.83%（3.32亿m3），外调水利用量占流域用水总量的46.92%（4.47亿m3），河道水利用量最小，仅有0.18亿m3，占流域用水总量的1.95%。全年的用水峰值集中在3月和8月，3月份是冬小麦用水的关键的时期，外调水利用量为0.86亿m3，其次为蓄水供水（0.42亿m3），由于地表水资源量相对充足，地下水仅用于补充供水，补充供水量为0.26亿m3。8月份为玉米生长用水的关键的时期，外调水利用量为0.77亿m3，由于蓄水资源量不足，供水量仅为0.2亿m3，地下水起到了重要的补充供水作用，补充供水量为0.75亿m3，在数量上几乎与外调水供水量相等。此外，在10-11月、1月出现用水小高峰，主要是由于种植冬小麦后的灌溉用水。由于这段时间处于枯水季节，地表水资源量有限，供水主要由丹江口水库供水（外调水）和地下水补充灌溉。
表1
| HRU代码 | 子流域属性 | 水资源区属性 | 行政区属性 | 灌区属性 |
| 字段格式 | 整形 | 整形 | 整形 | 整形 |
| 说明 | 子流域代码 | 水资源区代码 | 行政区代码 | 灌区代码 |
| 内容 | HRU所属子流域代码 | HRU所属水资源区代码 | HRU所属的行政区代码 | HRU所属的灌区代码 |
表2
| 参数名称 | 数据类型 | 说明 | 内容 |
| cntyID | 整形 | 行政区代码 | 行政区序号 |
| Riv1 | 整形 | 1号河道水源 | 河道所在的子流域代码 |
| Riv2 | 整形 | 2号河道水源 | 河道所在的子流域代码 |
| Riv3 | 整形 | 3号河道水源 | 河道所在的子流域代码 |
| Riv4 | 整形 | 4号河道水源 | 河道所在的子流域代码 |
| Riv5 | 整形 | 5号河道水源 | 河道所在的子流域代码 |
| Res1 | 整形 | 1号水库水源 | 水库代码 |
| Res2 | 整形 | 2号水库水源 | 水库代码 |
| Res3 | 整形 | 3号水库水源 | 水库代码 |
| Res4 | 整形 | 4号水库水源 | 水库代码 |
| Res5 | 整形 | 5号水库水源 | 水库代码 |
| Shal1 | 整形 | 1号浅水层水源 | 浅水层所在的子流域代码 |
| Shal2 | 整形 | 2号浅水层水源 | 浅水层所在的子流域代码 |
| Shal3 | 整形 | 3号浅层水水源 | 浅水层所在的子流域代码 |
| Shal4 | 整形 | 4号浅水层水源 | 浅水层所在的子流域代码 |
| Shal5 | 整形 | 5号浅水层水源 | 浅水层所在的子流域代码 |
| Deep1 | 整形 | 1号深水层水源 | 深水层所在的子流域代码 |
| Deep2 | 整形 | 2号深水层水源 | 深水层所在的子流域代码 |
| Deep3 | 整形 | 3号深层水水源 | 深水层所在的子流域代码 |
| Deep4 | 整形 | 4号深水层水源 | 深水层所在的子流域代码 |
| Deep5 | 整形 | 5号深水层水源 | 深水层所在的子流域代码 |
| Out1 | 整形 | 1号坑塘水源 | 坑塘所在的子流域代码 |
| Out2 | 整形 | 2号坑塘水源 | 坑塘所在的子流域代码 |
| Out3 | 整形 | 3号坑塘水源 | 坑塘所在的子流域代码 |
| Out4 | 整形 | 4号坑塘水源 | 坑塘所在的子流域代码 |
| Out5 | 整形 | 5号坑塘水源 | 坑塘所在的子流域代码 |
| Pnd1 | 整形 | 1号外流域水源 | 外源水代码 |
| Pnd2 | 整形 | 2号外流域水源 | 外源水代码 |
| Pnd3 | 整形 | 3号外流域水源 | 外源水代码 |
| Pnd4 | 整形 | 4号外流域水源 | 外源水代码 |
| Pnd5 | 整形 | 5号外流域水源 | 外源水代码 |
表2中若无供水，代码以0表示。
表3
| 参数名称 | 数据类型 | 说明 | 内容 |
| irrID | 整形 | 灌区代码 | 灌区序号 |
| Riv1 | 整形 | 1号河道水源 | 河道所在的子流域代码 |
| Riv2 | 整形 | 2号河道水源 | 河道所在的子流域代码 |
| Riv3 | 整形 | 3号河道水源 | 河道所在的子流域代码 |
| Riv4 | 整形 | 4号河道水源 | 河道所在的子流域代码 |
| Riv5 | 整形 | 5号河道水源 | 河道所在的子流域代码 |
| Res1 | 整形 | 1号水库水源 | 水库代码 |
| Res2 | 整形 | 2号水库水源 | 水库代码 |
| Res3 | 整形 | 3号水库水源 | 水库代码 |
| Res4 | 整形 | 4号水库水源 | 水库代码 |
| Res5 | 整形 | 5号水库水源 | 水库代码 |
| Shal1 | 整形 | 1号浅水层水源 | 浅水层所在的子流域代码 |
| Shal2 | 整形 | 2号浅水层水源 | 浅水层所在的子流域代码 |
| Shal3 | 整形 | 3号浅层水水源 | 浅水层所在的子流域代码 |
| Shal4 | 整形 | 4号浅水层水源 | 浅水层所在的子流域代码 |
| Shal5 | 整形 | 5号浅水层水源 | 浅水层所在的子流域代码 |
| Deep1 | 整形 | 1号深水层水源 | 深水层所在的子流域代码 |
| Deep2 | 整形 | 2号深水层水源 | 深水层所在的子流域代码 |
| Deep3 | 整形 | 3号深层水水源 | 深水层所在的子流域代码 |
| Deep4 | 整形 | 4号深水层水源 | 深水层所在的子流域代码 |
| Deep5 | 整形 | 5号深水层水源 | 深水层所在的子流域代码 |
| Out1 | 整形 | 1号坑塘水源 | 坑塘所在的子流域代码 |
| Out2 | 整形 | 2号坑塘水源 | 坑塘所在的子流域代码 |
| Out3 | 整形 | 3号坑塘水源 | 坑塘所在的子流域代码 |
| Out4 | 整形 | 4号坑塘水源 | 坑塘所在的子流域代码 |
| Out5 | 整形 | 5号坑塘水源 | 坑塘所在的子流域代码 |
| Pnd1 | 整形 | 1号外流域水源 | 外源水代码 |
| Pnd2 | 整形 | 2号外流域水源 | 外源水代码 |
| Pnd3 | 整形 | 3号外流域水源 | 外源水代码 |
| Pnd4 | 整形 | 4号外流域水源 | 外源水代码 |
| Pnd5 | 整形 | 5号外流域水源 | 外源水代码 |
表3中若无供水，代码以0表示。
表4
| 参数名称 | 数据类型 | 说明 | 内容 |
| cntyID | 整形 | 行政区代码 | 行政区序号 |
| sup_Riv1 | 整形 | 1号河道水源的取水序号 | 供水序号代码 |
| sup_Riv2 | 整形 | 2号河道水源的取水序号 | 供水序号代码 |
| sup_Riv3 | 整形 | 3号河道水源的取水序号 | 供水序号代码 |
| sup_Riv4 | 整形 | 4号河道水源的取水序号 | 供水序号代码 |
| sup_Riv5 | 整形 | 5号河道水源的取水序号 | 供水序号代码 |
| sup_Res1 | 整形 | 1号水库水源的取水序号 | 供水序号代码 |
| sup_Res2 | 整形 | 2号水库水源的取水序号 | 供水序号代码 |
| sup_Res3 | 整形 | 3号水库水源的取水序号 | 供水序号代码 |
| sup_Res4 | 整形 | 4号水库水源的取水序号 | 供水序号代码 |
| sup_Res5 | 整形 | 5号水库水源的取水序号 | 供水序号代码 |
| sup_Shal1 | 整形 | 1号浅水层水源的取水序号 | 供水序号代码 |
| sup_Shal2 | 整形 | 2号浅水层水源的取水序号 | 供水序号代码 |
| sup_Shal3 | 整形 | 3号浅层水水源的取水序号 | 供水序号代码 |
| sup_Shal4 | 整形 | 4号浅水层水源的取水序号 | 供水序号代码 |
| sup_Shal5 | 整形 | 5号浅水层水源的取水序号 | 供水序号代码 |
| sup_Deep1 | 整形 | 1号深水层水源的取水序号 | 供水序号代码 |
| sup_Deep2 | 整形 | 2号深水层水源的取水序号 | 供水序号代码 |
| sup_Deep3 | 整形 | 3号深层水水源的取水序号 | 供水序号代码 |
| sup_Deep4 | 整形 | 4号深水层水源的取水序号 | 供水序号代码 |
| sup_Deep5 | 整形 | 5号深水层水源的取水序号 | 供水序号代码 |
| sup_Out1 | 整形 | 1号坑塘水源的取水序号 | 供水序号代码 |
| sup_Out2 | 整形 | 2号坑塘水源的取水序号 | 供水序号代码 |
| sup_Out3 | 整形 | 3号坑塘水源的取水序号 | 供水序号代码 |
| sup_Out4 | 整形 | 4号坑塘水源的取水序号 | 供水序号代码 |
| sup_Out5 | 整形 | 5号坑塘水源的取水序号 | 供水序号代码 |
| sup_Pnd1 | 整形 | 1号外流域水源的取水序号 | 供水序号代码 |
| sup_Pnd2 | 整形 | 2号外流域水源的取水序号 | 供水序号代码 |
| sup_Pnd3 | 整形 | 3号外流域水源的取水序号 | 供水序号代码 |
| sup_Pnd4 | 整形 | 4号外流域水源的取水序号 | 供水序号代码 |
| sup_Pnd5 | 整形 | 5号外流域水源的取水序号 | 供水序号代码 |
表4中若无供水，代码以0表示。
表5
| 参数名称 | 数据类型 | 说明 | 内容 |
| irrID | 整形 | 灌区代码 | 灌区序号 |
| irr_Riv1 | 整形 | 1号河道水源的取水序号 | 供水序号代码 |
| irr_Riv2 | 整形 | 2号河道水源的取水序号 | 供水序号代码 |
| irr_Riv3 | 整形 | 3号河道水源的取水序号 | 供水序号代码 |
| irr_Riv4 | 整形 | 4号河道水源的取水序号 | 供水序号代码 |
| irr_Riv5 | 整形 | 5号河道水源的取水序号 | 供水序号代码 |
| irr_Res1 | 整形 | 1号水库水源的取水序号 | 供水序号代码 |
| irr_Res2 | 整形 | 2号水库水源的取水序号 | 供水序号代码 |
| irr_Res3 | 整形 | 3号水库水源的取水序号 | 供水序号代码 |
| irr_Res4 | 整形 | 4号水库水源的取水序号 | 供水序号代码 |
| irr_Res5 | 整形 | 5号水库水源的取水序号 | 供水序号代码 |
| irr_Shal1 | 整形 | 1号浅水层水源的取水序号 | 供水序号代码 |
| irr_Shal2 | 整形 | 2号浅水层水源的取水序号 | 供水序号代码 |
| irr_Shal3 | 整形 | 3号浅层水水源的取水序号 | 供水序号代码 |
| irr_Shal4 | 整形 | 4号浅水层水源的取水序号 | 供水序号代码 |
| irr_Shal5 | 整形 | 5号浅水层水源的取水序号 | 供水序号代码 |
| irr_Deep1 | 整形 | 1号深水层水源的取水序号 | 供水序号代码 |
| irr_Deep2 | 整形 | 2号深水层水源的取水序号 | 供水序号代码 |
| irr_Deep3 | 整形 | 3号深层水水源的取水序号 | 供水序号代码 |
| irr_Deep4 | 整形 | 4号深水层水源的取水序号 | 供水序号代码 |
| irr_Deep5 | 整形 | 5号深水层水源的取水序号 | 供水序号代码 |
| irr_Out1 | 整形 | 1号坑塘水源的取水序号 | 供水序号代码 |
| irr_Out2 | 整形 | 2号坑塘水源的取水序号 | 供水序号代码 |
| irr_Out3 | 整形 | 3号坑塘水源的取水序号 | 供水序号代码 |
| irr_Out4 | 整形 | 4号坑塘水源的取水序号 | 供水序号代码 |
| irr_Out5 | 整形 | 5号坑塘水源的取水序号 | 供水序号代码 |
| irr_Pnt1 | 整形 | 1号外流域水源的取水序号 | 供水序号代码 |
| irr_ Pnt 2 | 整形 | 2号外流域水源的取水序号 | 供水序号代码 |
| irr_ Pnt 3 | 整形 | 3号外流域水源的取水序号 | 供水序号代码 |
| irr_ Pnt 4 | 整形 | 4号外流域水源的取水序号 | 供水序号代码 |
| irr_ Pnt 5 | 整形 | 5号外流域水源的取水序号 | 供水序号代码 |
表5中若无供水，代码以0表示。
表6
| 参数名称 | 数据类型 | 说明 | 内容 |
| Cnty_surf(1,1) | 实型 | 1号行政区第1年的地表水取水量 | 取水量 |
| Cnty_surf(1,2) | 实型 | 1号行政区第2年的地表水取水量 | 取水量 |
| : | : | : | : |
| Cnty_surf(M,N) | 实型 | M号行政区第N年的地表水取水量 | 取水量 |
| Cnty_gw(1,1) | 实型 | 1号行政区第1年的地下水取水量 | 取水量 |
| Cnty_gw(1,2) | 实型 | 1号行政区第2年的地下水取水量 | 取水量 |
| : | : | : | : |
| Cnty_gw(M,N) | 实型 | M号行政区第N年的地下水取水量 | 取水量 |
表7
| 参数名称 | 调参方式 | 物理意义 | 调参值 |
| CN2 | r | 湿润条件II下的初始SCS径流曲线数 | 1.37 |
| GWQMN | v | 浅层含水层产生“基流”的阈值深度 | 1098 |
| GW_DELAY | v | 地下水补给延迟系数 | 10.4 |
| ALPHA_BF | v | 基流alpha因子 | 0.38 |
| ESCO | v | 土壤蒸发补偿系数 | 0.49 |
| EPCO | v | 植物吸收补偿因子 | 0.29 |
| GW_REVAP | v | 浅层地下水再蒸发系数 | 0.10 |
| REVAPMN | v | 浅层地下水再蒸发阈值 | 522 |
| <i>α</i> | v | 渠系输水损失蒸发比 | 0.09 |
| SOL_K | r | 土壤饱和水力传导度 | 0.52 |
| RCHRG_DP | v | 深含水层渗透比 | 0.15 |
| <i>β</i> | v | 渠系输水损失入渗比 | 0.50 |
| SOL_AWC | r | 土壤有效含水量 | 1.11 |
| <i>pip</i> | v | 管网漏损率 | 0.10 |
| <i>φ</i> | v | 渠系水有效利用系数 | 0.56~0.69 |
| <i>ω</i> | v | 田间水有效利用系数 | 0.95 |
| <i>r</i> | v | 居工地耗水率 | 0.2~0.45 |
| <i>v</i> | v | 污水处理率 | 1.0 |
| <i>re</i> | v | 再生水利用率 | 0 |
表7中，v表示调参值替代原参数值；r表示原参数值乘以调参值。
表8
表9
本发明以SWAT模型为基础，开发一种分布式自然-人工水循环模型，使其具有了自然-人工水循环动态互馈及水资源开发利用与调控一体化模拟功能，可完成基于城市多水源供水和多水源灌溉供水模式的流域自然-人工水循环过程模拟，在模型运行中，能够时刻保持自然水循环与人工侧枝水循环之间的动态互馈关系，不仅反映水循环过程对人工取用水的影响，同时还反映水资源开发利用与调控对水循环过程的实时干预作用，从而为深入了解强人类活动影响下的流域水循环机制，合理开发和利用水资源提供科学的参考依据。
以上所述，仅为本发明的具体实施方式，但本发明的保护范围并不局限于此，任何属于本技术领域的技术人员在本发明揭露的技术范围内，可轻易想到的变化或替换，都应涵盖在本发明的保护范围之内。因此，本发明的保护范围应该以权利要求的保护范围为准。
Claims (4)
1.一种基于多源互济供水模式的分布式水循环模型的设计方法，其特征在于，包括如下步骤：
步骤S1：采用基于“流域-水资源区-行政区-灌区-土地利用-土壤-坡度”多元属性的嵌套式坡面离散方法进行HRU划分，划分后的HRU具有相应的属性；
步骤S2：构建HRU属性识别模块，所述HRU属性识别模块用于识别HRU的属性；
步骤S3：设计多源互济供水模块，所述多源互济供水模块用于调用HRU属性识别模块，以识别HRU的属性，根据识别的属性判定土地利用类型、对应的供水源及供水源的供水优先次序，根据供水源的供水优先次序调用相应的水源模块进行取水；
步骤S4：将多源互济供水模块与SWAT模型中修改的模块进行衔接以实现数据实时交互，所述HRU根据输入的需水信息、供水水源类型信息、供用水优先序规则信息、水利工程信息和由SWAT模型提供的天然水文条件信息进行逐日水资源分配与调控，输出逐日“供水-用水-耗水-排水”的人工侧枝水循环信息，将人工侧枝水循环信息传递给SWAT模型；
所述步骤S3包括：
设计水源代码信息文件，所述水源代码信息文件用于读取指定的供水源信息，设置的供水源类型包括外调水、水库水、河道水、浅层地下水、深层地下水、坑塘水6种水源类型，水源代码信息文件通过程序指令读入；
设计供水优先次序信息文件，所述供水优先次序信息文件用于读取水源供水优先次序信息，指定供水源供水优先次序，供水优先次序信息文件通过程序指令读入；
设计取水控制量信息文件，所述取水控制量信息文件用于读取供水控制总量信息，识别HRU所属行政区的年度地表水供水控制总量和地下水开采控制总量，用于后续的水源的取水量计算；
设计多源互济供水计算流程，具体计算流程如下：
首先识别HRU的土地利用类型，如果是建设用地，则程序进入城乡供水计算流程；如果是农业用地，则程序进入灌溉用水计算流程；如果是其他土地利用类型，程序结束；
通过识别HRU的取水源识别代码，调用相应的水源模块，通过水源代码信息文件和供水优先次序文件的读取，识别每一个HRU的供水源个数、类型、取水先后顺序，根据取水源识别代码依次调用各水源模块；
根据HRU的供水次序寻求各水源进行取水，直至满足该HRU的当日的生活、工业用水需求、农业灌溉用水需求、或者到最后一个供水源供水完毕为止；
所述多源互济供水计算的步骤包括：
明确目标HRU设定的当日需水量WD；
明确目标HRU的供水水源个数k、水源代码及其供水优先次序，其中，k≤30；
按次序调用各水源模块，计算水源的取水量，水源模块包括rchuse模块、res模块、watuse模块、irr_rch 模块、irr_res模块和irrsub模块，水源取水量取决于当日HRU的需水量和水源的可供水量，可供水量则取决于水源的可利用水量、取水工程供水能力和取水控制量，计算公式如下：
式中，i表示HRU序号；j表示水源的供水优先次序号；WSP表示水源的当日实际取水量，单位是m3；WD表示HRU的当日需水量，单位是m3；Wsc表示水源的当日可供水量，单位是m3；WF表示取水工程的供水能力，单位是m3；WA表示水源的当日可利用水量，单位是m3；WMX表示年度取水控制量，单位是m3，WUM表示年用水控制总量，单位是m3；WSM表示年度地表取水控制总量，单位是m3；WGM表示年度地下水开采控制总量，单位是m3；
对供水优先次序为1的水源，则优先从该水源取水，如果水源的可供水量Wsc1＞WD，那么该水源供水量为WSP1=WD，供水程序结束，HRU的水源总供水量为WSP= WSP1；否则，WSP1=Wsc1，并且HRU的需水量变为Wf =WD-Wsc1，程序将继续寻找下一级供水源；
对供水优先次序为j的水源，j=2,…,k-1；k≤30，如果该水源的当日可供水量Wscj＞Wf，那么该水源供水量为WSPj=Wf，程序结束，HRU的水源总供水量为WSP=WSP+WSPj；否则，WSPj=Wscj，并且HRU的需水量变成Wf=Wf-Wscj，程序将继续寻找下一级供水源；
对供水优先次序为k的水源，k≤30，如果该水源的当日可供水量Wsck＞Wf，那么该水源供水量为WSPk=Wf，程序结束，HRU的水源总供水量为WSP=∑WSP i ；否则，WSPk=Wsck，并且HRU的需水量变成Wf=Wf-Wsck，程序结束。
2.根据权利要求1所述一种基于多源互济供水模式的分布式水循环模型的设计方法，其特征在于，所述步骤S1包括：
采用ArcGIS对DEM提取流域河网，划分天然子流域；
在天然子流域上叠加土地利用信息、土壤类型信息和坡度信息，划分出天然HRU；
将天然HRU套上行政区、水资源区边界，将天然HRU进一步切割；
根据灌区的分布，将灌区与天然HRU进行叠加，最终完成HRU划分，这样，每一个HRU都具有子流域属性、水资源区属性、行政区属性、灌区属性、土地利用类型属性和土壤类型属性。
3.根据权利要求2所述一种基于多源互济供水模式的分布式水循环模型的设计方法，其特征在于，所述步骤S2包括：
构建HRU属性识别模块，用于读取指定的HRU的属性，所述指定的HRU的属性包括子流域属性、水资源区属性、行政区属性、灌区属性；
将构建的HRU属性识别模块放在SWAT模型中的主模块中，便于对HRU属性识别模块的调用。
4.根据权利要求1所述基于多源互济供水模式的分布式水循环模型的设计方法，其特征在于，所述步骤S4中修改SWAT模型中的相关模块具体包括：
屏蔽rchuse 模块、res模块、watuse模块、irr_rch 模块、irr_res模块和irrsub模块，不再采取单一水源取水模式；将以上模块放进多源互济供水模块中进行调用；
添加相关代码，将参数WSP i 分别替换掉rchuse 模块、res模块中的参数waterrch、wuresn，实现多源互济供水模块Multi_sc与rchuse 模块、res模块的衔接和调用，其中i=1,2；
修改相关程序，增加外调水取水和外调水量约束功能，使供水量不能超过调水限制总量，计算公式表达如下：
其中，waterout (i,j)表示第i天、第j个HRU的调水利用量，单位是m3；WX5表示调水限制总量，单位是m3；
在watuse模块中添加代码，用参数WSP i 分别替换掉watuse模块中的参数watershal、waterdeep、waterout和waterpnd，实现多源互济供水模块Multi_sc与watuse模块的衔接与调用，其中i=3,4,5,6；
在rchuse 模块、res模块、watuse模块中添加如下公式的计算程序：
式中，pip表示供水管网漏损率；
添加坑塘灌溉功能，并完善外调水灌溉功能，增加供水量限制约束，使灌溉取水量不能超过调水限制总量：
其中，wirrout (i,j)表示第i天、第j个HRU的调水灌溉利用量，单位是m3；MX5表示调水限制总量，单位是m3，wirrpnt (i,j)表示第i天、第j个HRU的坑塘灌溉利用量，单位是m3；MX6表示坑塘可供水量，单位是m3；
修改irr_rch模块、irr_res模块、irrsub模块的源代码，增加渠系输水过程模拟，包括渠道水量损失和渠道退水，其中，渠道水量损失包括渠道内的水分蒸发损失和渠道渗漏损失两部分，主要计算公式如下：
式中，ET can 表示渠系蒸发损失量，单位是mm；IRR can 表示进入渠道的灌溉水量，单位是mm；Ls can 表示渠系渗漏损失量，单位是mm；Surp can 表示渠系退水量，单位是mm；
表示渠系水有效利用系数；
表示渠系蒸发系数；
表示渠系渗漏系数；
修改percmain模块相关代码，添加渗漏损失的计算程序，计算公式如下：
式中，
表示第t天第一层土壤的土壤含水量，单位是mm；
表示第t-1天第一层土壤的土壤含水量，单位是mm；inf pcp 表示降雨入渗量，单位是mm；inf irr 表示灌溉入渗量，单位是mm；inf wet 表示湖库湿地渗漏量，单位是mm；
修改gwmod模块中的地下水补给代码，实现对管网漏损水量的模拟，计算公式如下：
式中，
表示第t天的潜水补给量，单位是mm；
表示第t-1天的潜水补给量，单位是mm；prc表示补给地下水的土壤渗漏水量，单位是mm；GW_DELAY表示地下水补给延迟系数，单位是mm；Area表示HRU的面积，单位是m2；
在subbasin模块中添加多源互济供水模块，通过读取指定的水源类型、水源个数、水源代码、取水量、取水时间，进行年内的水源动态互补供水运算，实现SWAT模型运行时的多水源联合供水模拟；
通过修改surface模块相关代码，将渠系退水叠加地表产流量，参与河道汇流计算，相关计算公式如下：
式中，
表示渠道退水后的径流量，单位是mm；
表示渠道退水前的径流量，单位是mm；
点源模块包括recday模块和recmon模块，在recday模块和recmon模块中修改相关代码，将污染排放参数WDR分别替换参数floday、flomon，相关计算公式如下：
其中，WDR为城市污水产生量，单位是m3；WP为污水排放量，单位是m3；r表示用水消耗率；v表示污水处理厂的污水处理率；re表示再生水利用率；
将构建的HRU属性识别模块放在SWAT模型中的main主模块中，便于对HRU属性识别模块的调用。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110946366.8A CN113407897B (zh) | 2021-08-18 | 2021-08-18 | 基于多源互济供水模式的分布式水循环模型的设计方法 |
| US17/867,221 US20230099257A1 (en) | 2021-08-18 | 2022-07-18 | Design method for distributed hydrological cycle model based on multi-source complementary water supply mode |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110946366.8A CN113407897B (zh) | 2021-08-18 | 2021-08-18 | 基于多源互济供水模式的分布式水循环模型的设计方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113407897A CN113407897A (zh) | 2021-09-17 |
| CN113407897B true CN113407897B (zh) | 2021-10-29 |
Family
ID=77688848
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110946366.8A Active CN113407897B (zh) | 2021-08-18 | 2021-08-18 | 基于多源互济供水模式的分布式水循环模型的设计方法 |
Country Status (2)
| Country | Link |
|---|---|
| US (1) | US20230099257A1 (zh) |
| CN (1) | CN113407897B (zh) |
Families Citing this family (11)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN113837666B (zh) * | 2021-11-23 | 2022-02-11 | 长江水利委员会长江科学院 | 一种基于改进作物水分指数的农业干旱评估方法 |
| CN115062447B (zh) * | 2022-04-25 | 2023-03-10 | 中国科学院地理科学与资源研究所 | 一种基于swat模型模拟水文环境的方法及系统 |
| CN114673214B (zh) * | 2022-04-28 | 2023-08-18 | 科源智水(北京)科技有限公司 | 多水厂智能联动供水控制系统及其控制方法 |
| CN116202576B (zh) * | 2023-05-06 | 2023-09-26 | 广东广宇科技发展有限公司 | 一种城市水利的河湖动态管理的监测方法及系统 |
| CN116663838A (zh) * | 2023-06-08 | 2023-08-29 | 长江水利委员会长江科学院 | 基于城市建成区的分布式水文模型、构建及供水模拟方法 |
| CN116402410B (zh) * | 2023-06-08 | 2023-08-25 | 长江水利委员会长江科学院 | 一种分布式水量水质配置方法及系统 |
| CN116822877B (zh) * | 2023-06-29 | 2024-03-22 | 浪潮智慧科技有限公司 | 一种基于知识图谱的水资源分配方法、设备及介质 |
| CN116562475B (zh) * | 2023-07-11 | 2023-09-26 | 北京师范大学 | 一种获取城市水资源消耗和水资源韧性优化参数的方法和系统 |
| CN117114347B (zh) * | 2023-09-12 | 2024-04-26 | 深圳市环境科学研究院 | 一种潮汐河流生态补水优化配置方法 |
| CN117237677B (zh) * | 2023-11-15 | 2024-03-26 | 南京信息工程大学 | 一种基于深度学习的强降水空间整体相似度的降水预报订正方法 |
| CN117893179A (zh) * | 2024-03-18 | 2024-04-16 | 四川嘉源生态发展有限责任公司 | 基于数字孪生的水利电力工程建造管控方法及系统 |
Family Cites Families (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20200026735A1 (en) * | 2018-07-19 | 2020-01-23 | Bank Of Montreal | Systems and methods for digitization services |
| CN111899126B (zh) * | 2020-08-10 | 2023-12-15 | 中国水利水电科学研究院 | 基于水循环模拟的三条红线控制指标划分方法 |
| CN112199456B (zh) * | 2020-09-17 | 2022-03-25 | 西南科技大学 | 一种基于供给型水文生态系统服务的流域水资源管理方法 |
| CN112765800B (zh) * | 2021-01-12 | 2022-10-21 | 长江水利委员会长江科学院 | 一种分布式水资源配置模型的设计方法 |
-
2021
- 2021-08-18 CN CN202110946366.8A patent/CN113407897B/zh active Active
-
2022
- 2022-07-18 US US17/867,221 patent/US20230099257A1/en active Pending
Also Published As
| Publication number | Publication date |
|---|---|
| CN113407897A (zh) | 2021-09-17 |
| US20230099257A1 (en) | 2023-03-30 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113407897B (zh) | 基于多源互济供水模式的分布式水循环模型的设计方法 | |
| CN112765800B (zh) | 一种分布式水资源配置模型的设计方法 | |
| CN103544533B (zh) | 一种调水工程水资源优化配置动态模型构建系统 | |
| CN105528734A (zh) | 一种多层水资源管理的水量分配方法 | |
| Multsch et al. | Improving irrigation efficiency will be insufficient to meet future water demand in the Nile Basin | |
| CN108108556B (zh) | 一种基于耗散-汇合结构的灌区水循环模型构建方法 | |
| CN108764573A (zh) | 一种面向内陆干旱区的水资源多维均衡配置系统 | |
| Omar | Evaluation of actions for better water supply and demand management in Fayoum, Egypt using RIBASIM | |
| CN111062537A (zh) | 临河灌区节水量计算方法与系统 | |
| Sivakumar | Management policy of water table in dry zone of Sri Lanka to subsidise the pain of non rice crop cultivators for the food productivity improvement. | |
| Ojeda-Bustamante et al. | Using spatial information systems to improve water management in Mexico | |
| Sun et al. | Response of groundwater to different water resource allocation patterns in the Sanjiang Plain, Northeast China | |
| CN116402410B (zh) | 一种分布式水量水质配置方法及系统 | |
| CN116956407A (zh) | 一种污染物的拦截与消纳方法 | |
| Srivastava et al. | GIS based integrated modelling framework for agricultural canal system simulation and management in Indo-Gangetic plains of India | |
| Wang et al. | An integrated model for simulating regional water resources based on total evapotranspiration control approach | |
| Sidhu et al. | A study of availability and utilization of water resources in Punjab | |
| CN111738554A (zh) | 基于耗水量的山水林田湖草各系统均衡的水资源配置方法 | |
| CN115496622B (zh) | Mikebasin水资源配置模型文件快速生成方法 | |
| Nguyen | Harmonizing multi-sectorial water management with minimum flow requirements in an anthropogenically impacted river basin. The case of Vu Gia–Thu Bon, Central Viet Nam | |
| Yaozhou et al. | Pricing of irrigation water in China | |
| Fipps et al. | Implementation of a district management system in the Lower Rio Grande Valley of Texas | |
| Boels et al. | Analysis of water management in Sirsa District in Haryana: model testing and application | |
| Esnaashari et al. | Modeling the effects of changes made in Sistan's water resources system using WEAP | |
| Zehui | ECOLOGICAL RESTORATION BASED ON A REGIONAL POLYCENTRIC GOVERNANCE MODEL--REGENERATION OF THE KAREZ WATER SYSTEM IN TURPAN REGION. |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |