1. 引言
随着社会经济的快速发展，人类活动对水文水资源系统的干预越来越强烈，人类活动与自然水文水资源系统之间的相互关系日趋紧密，动摇了传统水文水资源学科理论和方法的科学基础，基于历史长系列数据的水文模拟分析和水资源管理面临严峻的挑战。国际地圈–生物圈研究计划(IGBP)和全球水系统计划(GWSP)已将气候变化和人类活动对水循环产生的影响作为重点研究问题之一；联合国教科文组织的国际水文计划第七阶段(IHP-VII)任务中，把气候变化对水文循环和水资源的影响以及提高变化环境下水资源管理能力作为研究的重点领域[1] 。国际水文科学协会(IAHS)2013年正式启动Panta Rhei (2013~2022)科学计划[2] ，主题是变化中的水文与社会(Change in Hydrology and Society)，该计划提出了6个科学问题，即：① 人们认知水文变化的主要差距在哪？② 水文系统变化如何通过水文过程驱动与自然、社会系统相互作用和反馈？③ 什么是水文系统变化的外部驱动力和内部系统特征？④ 怎样运用水文和社会耦合系统新知识来提升预测能力，包括预测结果的不确定性计算、可预测性评价？⑤ 如何增强水文监测和数据分析能力、以便预测和管理水文过程变化？⑥ 怎样考虑水文自然变化与人类活动影响的相互反馈和不确定性，支撑经济社会发展适应这些变化条件？通过对上述6个科学问题进行研究探讨，以实现科学认知、预测评估、实际应用三大目标，构成水文科学研究的整体平台，应对处于正在变化中的环境系统和经济社会发展的需求。
我国人多水少、水资源时空分布严重不均，洪涝灾害、水资源短缺、水污染严重、水生态恶化等问题十分突出，已成为制约经济社会可持续发展的主要瓶颈。国务院2012年发布《关于实行最严格水资源管理制度的意见》[3] ，确立水资源开发利用控制、用水效率控制和水功能区限制纳污“三条红线”，实施包括用水总量控制、用水效率控制制度、水功能区限制纳污制度和水资源管理责任和考核制度在内的“四项制度”，从制度上推动经济社会发展与水资源、水环境承载能力相适应[4] 。
可以看出，国际国内社会近年来对如何应对未来可能发生的环境变化，改变现有水文模拟预测理论方法和资源管理方式，以及提出相关适应性对策等方面的研究给予了很大的关注。因此，针对水文循环和水资源系统动态演化，开展以水循环模拟为基础的区域动态水资源安全评价、水资源可持续利用和适应性管理等重大需求的科学研究，也必然成为当前水文水资源科学研究的前沿与热点问题。
2. 国内外研究现状与进展
变化环境下水文水资源与社会耦合系统的研究是一个新的多学科交叉的研究领域。它包括水文系统变化规律研究、水文过程与自然、社会系统相互作用和反馈，水文和社会耦合系统的建立，水资源优化配置与管理等等。与其相关的研究包括变化环境下水文水资源分析、计算与模拟的理论和方法，水文模型与气候模型的耦合可能性与挑战，水文水资源和社会耦合系统的概念与模型，等。下面首先综述国内外研究现状与进展。
2.1. 变化环境下水文分析与模拟
除地球系统自然条件的变化外，流域或区域内经济快速发展、人类活动影响加剧、水利工程的兴建、地下水的开发、林业植被、矿山开采、坡地治理、河道采砂等人类活动导致下垫面条件不断改变，使流域降雨径流等关系发生了巨大变化。张建云等[5] 研究发现：近50年我国六大江河(长江、黄河、珠江、松花江、海河、淮河)的实测径流量都呈下降趋势，水文特征值发生了显著改变，径流下降幅度最大的是海河流域，径流组成及年内分配均发生了一定程度的变化。流域水文预测预报通过建立降雨径流的数学关系定量模拟径流过程，研究和探讨降雨径流关系是否稳定是水文预测预报非常关键的一步，目前主要的分析手段有统计学方法和建立流域水文模型两种。统计学方法所包含的分析方法较广泛，其原理通常是假设研究变量为随机变量，通过计算样本的均值、标准差、方差、径流系数、变异系数和显著性检验以及年际变化绝对比率、不均匀系数、年径流和年降水距平百分率的变化等指标来描述降雨径流关系[6] [7] 。任立良等[8] 对中国北方的黄河、海河、辽河、松花江流域的长时间系列观测资料进行分析，认为这些流域的径流都存在不同程度减少，其中上游水库拦蓄和人类取用水增加是造成径流减少的主要原因，计算发现1980s~1990s相同量级的降水量产生的径流量比1950s~1960s减少了20%~50%。冯平等[9] 对潘家口水库1956~2003年的入库径流变化趋势及影响因素进行了分析，通过数据分析，认为下垫面变化是导致径流减少的主要原因。李慧赟等[10] 采用模糊聚类的方法对1942~2005年年降雨序列聚类分析，研究不同时间分布类型的降雨下人类活动对流域水文效应的影响。结果表明，人类活动的影响是造成流域径流减少的主要原因，并且降雨量越小，年内分布越分散，流域径流受人类活动影响越剧烈。李致家等[11] 以湿润、半湿润、半干旱以及干旱地区的15个典型流域为研究对象，分别建立其P + Pa~Rs相关关系，然后对其中3个嵌套流域的P + Pa~Rs相关关系进行分析与验证，比较各流域嵌套前后的P + Pa~Rs相关关系变化，对不同自然地理条件下降雨径流相关关系的区域性规律进行探讨，李二辉等[12] 应用Mann-Kendall秩次相关检验、流量历时曲线法、双累积曲线法等方法对黄河干流陕县站和河口镇站1919~2010年径流量演变过程进行了分析。结果表明：区域面平均降水量趋势性变化不显著，而上游及中游年径流量自1985年以来呈显著减少趋势，中游径流量的降幅高于上游，并定量分析了降水和人类活动对径流量的影响，上游和中游人类活动对径流量减少的影响程度分别占88.1%和84.9%，水利水保工程、生产生活用水等人类活动是引起黄河径流量减少的主要因素。
流域水文模型在水文预报和解决生产实际问题中起着重要的作用，如防洪减灾、水资源可持续利用、水环境和生态系统保护等。流域水文模型还是分析研究气候变化和人类活动对洪水，水资源和水环境影响的有效工具，该方法用人类活动影响之前的资料率定流域水文模型参数，用率定后的模型来推求人类活动影响之后的流域水文过程，并与实测资料进行对比。陈利群和刘昌明[13] 利用2个分布式水文模型(SWAT和VIC)模拟了黄河源区土地覆被变化与气候波动对径流的影响，并区分两者的影响贡献率。王忠静等[14] 分析了1957~2000年期间永定河官厅水库以上流域的还原水量与入库水量，并采用GBHM模型模拟了1990和2000年两期土地利用下的流域径流，分析了下垫面变化对径流的影响。权全等[15] 运用DRAINMOD水文模型模拟了1988~2007年西安市城市土地利用变化对降雨径流产生的影响，探讨了降低城市地表径流、增加雨水入渗的调控措施，结果表明土地利用类型、土壤入渗能力以及降雨强度是影响城区降雨–径流关系的重要因素。史晓亮等[16] 以滦河流域为研究区，利用1985和2000年土地利用数据，结合SWAT分布式水文模型定量评价了流域土地利用/覆被变化的水文效应，并分析了流域地表径流变化与主要景观类型的响应关系。贺瑞敏等[17] 利用VIC模型，在海河流域选取了6个典型流域资料来率定VIC模型的参数，通过模型参数移植技术，建立了全流域的径流模拟平台，根据假定的气候变化情景，分析了海河流域河川径流对气候变化的响应机理。许崇育等[18] 全面回顾总结了水文模型发展的历史，取得的成就和存在的问题，重点讨论了环境变化下水文模拟的几个主要问题：①水文过程和降雨–径流关系的非稳定性问题；②水文模型在无资料流(区)域的应用问题；③水文模型在变尺度(时间和空间)情况下存在的问题；④水文模型与气候模型耦合的可能性问题。在讨论问题的同时提出了作者独到的观点和想法。
许多水文学家对如何甄别水文序列的变化和探讨其变化原因做了大量的研究工作，也得出了不同的结论[19] -[21] 。Zhang等[22] 初步研究和比较了中国九大流域的降雨和流量近50年(1956~2005)的变化趋势以及降雨径流关系的变化规律。结果表明，在中国北方干旱和半干旱地区(黄河，海河，塔里木河，松辽河)，降雨和径流明显呈现出不一致的变化趋势，导致了降雨径流关系的非平稳性。而在中国南方湿润地区(长江，珠江，岷江，及西南诸河)，虽然降雨和径流都有不同程度的变化，但降雨和径流的变化趋势基本一致，说明南方诸河的降雨径流关系趋于稳定。
2.2. 变化环境下水文模型结构和参数随时间和空间转移的可能性
在水文序列和降雨径流关系稳定的前提下建立的水文模型，被广泛用来模拟气候变化和土地利用变化情景下的水文响应。一般的计算方法和步骤是：①通过历史资料对模型参数进行率定；②利用不同的方法和模型建立未来气候和土地利用变化情景；③将第二步产生的情景，作为第一步率定好模型的输入信息和资料，模拟生成变化环境下的径流和其它水量平衡项。现有研究中存在的核心问题是水文模型的结构和参数是否具备在变化环境下的转移应用能力，而大量的研究表明：大多数的水文模型并不具备这个功能[23] -[26] 。主要表现在即使是非常相似的模型，而且其在率定期模拟结果相差无几，一旦把变化环境下非稳定性水文序列作为模型的输入，其结果相差巨大[27] 。
如何才能得到变化环境情景下较为客观的水文响应结果？最简单的也是最常用的方法就是采用多个水文模型的结果加以平均，也就是简单的集合预测方法。该方法虽然能在一定程度上降低单一模型的不确定性，但由于它不能确定不同模型的优劣，导致了它们的平均值仍可能有很大的不确定性，而且这个不确定性的程度是未知的。Young [28] 发现模型参数随率定时段的不同而发生显著变化。也就是说，如果用不同时期资料率定的模型参数，模拟环境变化情景下的水文响应将得出不同的结果。Merz等[29] 在奥地利选择273个流域，用41年的资料对HBV模型每5年率定一次，结果发现不管是模型的效率系数还是模型参数值都随时间发生明显变化。由于未来的情景是不可测的，环境变化下水文预测预报的结果也无法检验。国际水文科学协会(IAHS)前主席Klemes [30] 提出了利用历史资料检验水文模型模拟变化环境下水文响应的方法，即differential-split sample检验，不能通过这个检验的模型不适合用来做这方面的应用。Xu [31] [32] 采用WASMOD模型开展了类似检验的研究，但由于流域个数和序列长度有限，还有待对其它模型在不同区域的应用效果作进一步的验证。Li等[33] 利用differential-split sample检验法研究了DWBM和SIMHD模型参数在非稳定气候条件下移植可能性，结论是不同气候条件下率定的参数在交叉应用时将产生显著误差。Refsgaard等[34] 比较了几种在气候变化情景下水文模型参数转移的方法，得出通过differential-split sample检验的模型要优于其它模型的结论。
无(缺)资料流(区)域模型参数是需要通过移用有资料流域的参数来实现的。因为每个流域的气候和下垫面条件不完全相同，所以无(缺)资料流(区)域水文模拟也属于变化环境下的水文模拟问题。预测同一个流域土地利用变化下的水文响应可以认为是其中一种特殊的情景。无(缺)资料流(区)域的水文模拟经历自有水文模型以来的多年探索，特别是国际水文科学协会(IAHS)2002年启动的无(缺)资料流域的水文预测(Prediction in Ungauged Basin，简称PUB) (2003~2012)研究计划。该计划实施10余年来，取得了一些重要的理论和应用成果[35] [36] 。无(缺)资料流(区)域的模型参数估算目前常见的方法是多元回归法[37] -[39] 、流域空间邻近法和流域物理特性相似法[40] -[42] 、全局平均法[43] [44] 和多目标区域率定法[45] [46] 。
2.3. 水文模型和气候模型耦合的问题
为延长水文预报的预见期，20世纪90年代以来，越来越多的水文气象学者将数值天气预报模式与水文模型进行耦合研究，根据水文模型的输出结果是否对数值天气预报模式形成反馈，可将耦合方法分成两类：单向耦合和双向耦合。Lin等[47] 采用双向耦合(MC2/CLASS*/GUH)系统成功地重现了一场暴雨洪水过程线，显示了大气水文耦合模式预报暴雨洪水的能力，并能提供较长的预见期；Seuffert等[48] 将中尺度数值天气预报系统LM和陆面水文模型TOPLATS双向耦合，用TOPLATS取代了LM中土壤模块，克服了未考虑土壤侧向水流及汇流过程的不足，在同一个交互界面上实现了两个模式的在线耦合，提高了大气模式对陆气之间的能量通量和降水的预报能力。Anderson等[49] 采用单向的方法耦合了中尺度预报模式MM5与HEC-HMS水文模型，并对Calaveras流域3日的1小时连续流量过程进行了预报，同时用图示法比较了MM5降水预报与实测值的差别，指出其降雨预报精度还有待提高；Karsten等[50] 将5个数值天气预报模式和一个陆面分布式水文模型WaSiM-ETH分别作了单向耦合，研究耦合模式在洪水预报中的应用，他们指出耦合预报的结果因数值预报模式和挑选场次洪水不同而不同，预报降水的位置和雨量即使只有相对较小的误差，也将导致流量预报的巨大误差，耦合模式的改进将主要依靠大气模式的改进；Walter Collischonn等[51] 用中尺度数值预报模式ARPS模拟预报了连续167天的3种不同空间分辨率下的降水资料，以此驱动分布式水文模型模拟预报径流过程，并与零降雨输入、完美预报输入(即实测降雨)、实测径流过程三种情况进行了比较，但未对ARPS的降水预报进行评定；陆桂华等[52] 利用加拿大中尺度数值模式MC2模拟的降水驱动半分布式新安江模型，建立了降雨预报评价标准，进行了产汇流计算，结果显示耦合模型能成功的预报研究场次的暴雨和洪水过程，且能有效的延长洪水预见期。郭生练等[53] 建立了中尺度数值预报模式MM5和水文模型VIC的单向耦合预报系统，并采用汉江流域汛期水文资料对模型预报效果进行了验证。高冰等[54] 采用新一代中尺度数值天气模式WRF和分布式水文模型GBHM，开展了基于气象水文耦合的三峡入库洪水预报研究，结果表明，基于数值天气预报的三峡入库洪水预报具有一定精度，并且能够显著延长洪水预报的预见期。
2.4. 水资源优化配置模型与水文模型的耦合
20世纪90年代以来，为推动水资源系统规划的发展与应用，国外开发了一系列水资源系统模拟软件，如：丹麦水利与环境研究所(DHI)的Mikebasin、美国Brigham Young大学与陆军工程兵团共同开发的WMS、奥地利环境软件与服务公司开发的流域综合管理软件Waterware、以美国农业部(USDA)为主开发的流域水资源模拟模型Aquarius、澳大利亚研制的ICMS (Interactive Component Modeling System)水资源系统管理模型等，水资源配置对象从单纯的水量配置逐步发展为考虑水质因素的水资源配置。与此同时，国内水资源配置研究始于20世纪60年代以水库优化调度为手段的水资源分配研究，经过一系列的国家科技攻关、自然科学基金等项目的研究以及生产实践的推动，水资源配置理论体系逐步丰富与完善。黄强等[55] 以西安市水资源系统为对象，建立多水源、多目标水资源配置模型，提出了相应的求解方法，获得了系统有限水量在不同用水部门的优化分配方案。陈晓宏等[56] 以大系统分解协调理论作为技术支持，运用逐步宽容约束法及递阶分析法，建立东江流域水资源优化调配的实用模型，并对该流域特枯年的水资源量进行了优化配置和供需平衡分析。王浩等[57] 在“黄淮海水资源合理配置研究”中，首次提出水资源“三次平衡”的配置思想，系统地阐述了基于流域水资源可持续利用的系统配置方法，其核心内容是在国民经济用水过程和流域水循环转化过程两个层面上分析水量亏缺态势，并在统一的用水竞争模式下研究流域之间的水资源配置问题。刘建林等[58] 以系统分析理论为基础，对跨流域多水源多目标多工程调水所涉及的水资源问题进行研究，建立了南水北调东线工程联合调水仿真模型，为实际应用提供了决策平台。刘丙军等[59] 针对南方地区水资源利用特点，以东江流域为例，建立了以水资源合理利用为目标，综合考虑防洪、供水、河道生态等约束的水资源合理配置多目标分析模型，并运用大系统“分解协调”原理提出了相应算法。王德智等[60] 针对动态规划在供水库群优化调度中存在“维数灾”，且难以获得真正最优解的问题，将改进的加速遗传算法应用到水库群的水资源优化配置之中。金菊良等[61] 为确定区域水资源合理配置各评价指标的权重、有效处理区域水资源分配过程中的模糊性和随机性，提出了基于信息熵原理的模糊模式识别模型与基于加速遗传算法的模糊层次分析法相耦合，进行各子区域水资源量合理配置的新模型(EFPR-FAHP)。辛芳芳和梁川[62] 在多目标线性规划求解中引入模糊数学中的隶属度概念，建立以目标隶属度最大为目标的数学模型，应用于都江堰灌区水资源合理配置中。陈卫宾等[63] 提出基于记忆梯度混合遗传算法用于灌区水资源优化配置。耿福明等[64] 根据水资源价值、边际机会成本、最优化数学等理论，基于水资源净效益最大化的水资源优化配置思想，以水资源承载力、用水部门生存条件、用水公平性及可持续发展约束等水资源开发利用约束为准则，建立了基于水资源净效益最大化的水资源优化配置模型，并将该模型应用于南水北调受益区。王福林和吴丹[65] 基于区域产业结构发展趋势及其演变规律，以保障区域生活用水需求与生态环境用水需求为前提，在区域水资源可供给总量与产业发展用水需求量的约束条件下，建立基于水资源优化配置的区域产业结构动态演化模型。陈春林、郑垂勇等(2010)针对传统权重确定方法的不足，提出了一种新的确定权重的方法(ANP方法)，将该方法应用于南水北调工程受水区水资源优化配置。陈强等[67] 利用改进的SWAT模型多水源灌溉模块，将水资源配置模型的农业灌溉用水分布到改进后的SWAT模型中，实现二模型的松散式耦合，这种分布式水文模型与配置模型的耦合为流域水资源管理提供了一种新的途径。张礼兵等[68] 采用对数平均迪氏指数法(LMDI)建立了工业用水恒等式，将工业用水影响效应划分为工业节水效应、经济规模效应、产业结构效应和用水强度效应，对安徽省2003~2011年总体及逐年的工业行业用水变化分别进行因素分解计算和分析。左其亭等[69] 从和谐论的观点出发，把和谐平衡定义为“利益相关者考虑各自利益和总体和谐目标而呈现的一种相对静止且相关者各方暂时都能接受的平衡状态”，提出了基于和谐论的和谐平衡计算方法，确定了实现河南省水资源与经济社会和谐平衡的条件。张守平等[70] [71] 构建了供需平衡、耗水平衡和基于水资源优化配置的水质模拟系统，提出基于水功能区纳污能力的污染物总量分配优化模型，分析分质供水的实现方式以及不同模块之间的耦合关系。李昱等[72] 根据水库群的拓扑结构和共同供水目标的分布情况，分层构建了两个虚拟聚合水库及其联合调度图；并根据水库自身属性以及供水任务，联合使用分水比例法和补偿调节法两种方法分配共同供水任务；同时设置补给限制线以避免过度补偿，以此构建水库群联合调度模型，并采用改进遗传算法对模型进行优化求解，得到相应联合调度规则。
综上所述，水资源配置研究主要表现在[73] -[75] ：① 配置研究的水源范围上，从最初的地表水量分配为主，发展到地表水、地下水联合调度配置，再到对常规水源和非常规水源的统一调配，从一次性水资源到再生性水资源的配置；②配置目标上，从单一的供水效益最大化发展到对流域水资源管理目标的多属性识别；③ 配置对象上，从供水量的配置发展到对取水量与耗水量的统一分配，从水量分配到考虑水量水质的联合分配；④ 配置指导思想上，从“以需定供”和“以供定需”到基于宏观经济的水资源配置，再到面向可持续发展配置；⑤ 配置方式上，从单一的工程措施扩展到采用经济、协商等多种手段并用的配置方式；⑥ 配置模型构建上，从“还原论”为基础的组合模型方法到整”体模型方法，系统动力学、大系统分解协调、多目标规划与决策技术等都引入到模型的构建和分析中；⑦ 配置模型求解技术上，根据所构建模型特点，已从线性规划、动态规划发展到多目标求解技术以及各种启发式智能算法，如模拟退火算法(SA)、蚁群算法(ACO)、遗传算法(GA)、粒子群算法(PSO)、人工免疫算法、人工鱼群算法、和声算法、多目标广义演化免疫算法(MEMOIA)、第二代非支配排序遗传算法(NSGA-II)等。
2.5. 流域水循环“自然–社会”二元变化的驱动机理与模拟
随着社会经济的快速发展，人类活动对流域水文循环过程和水资源演变规律的影响逐渐加剧，引发了一系列严重的水问题和水危机，其主要原因是流域水循环的驱动力、循环结构、循环参数在人类活动的影响下发生了二元演变效应，呈现出明显的“自然–社会”二元特性，并导致了径流性水资源衰减、水体污染和环境污染、天然生态退化等各种不良后果[76] [77] 。因此，流域水文循环过程的模拟不仅仅要研究变化环境下的流域降雨径流关系，更需要研究人类活动影响下的流域水文效应和人类活动对水循环要素变异的驱动机理。王浩等[78] 归纳了流域水循环“自然–社会”二元特性主要表现在：一是循环驱动力的二元化，即流域水循环的内在动力已由过去一元自然驱动演变为现在的“自然–社会”二元驱动；二是循环结构的二元化，即人类聚集区的水循环过程往往由自然循环和人工侧支循环耦合而成，两大循环之间保持动力关系，通量之间此消彼涨；三是水资源服务功能的二元化，即水分在其循环转化过程中，同时支撑了同等重要的经济社会系统和生态环境系统。贾仰文等[79] [80] 针对高强度人类活动作用下海河流域水循环的“自然–社会”二元特性，开发了流域二元水循环模型，该模型由分布式流域水循环模型(WEP)、水资源合理配置模型(ROWAS)和多目标决策分析模型(DAMOS)3个模型耦合而成，针对各模型的优势与不足，提出了两层耦合的技术路线，以实现统筹考虑水资源、宏观经济与生态环境的流域水资源综合管理分析的功能。刘家宏等[81] 以海河流域为原型对象，分析了二元水循环的关键要素、平衡方程以及演变历程，应用建立的二元水循环模式理论，研究了海河流域二元水循环十项关键要素的演化规律。王浩等[82] 在全面剖析流域水循环演变“自然–社会”二元特性的基础上，构建海河流域二元水循环及其伴生过程综合模拟平台，并以此为工具通过基于气候模式预测结果与水量控制条件的情景设置，对水资源、水生态和水环境的未来演变情势进行了预估。刘佳嘉等[83] 提出一种基于分布式水文模型的分解气候变化和人类活动对流域水循环过程影响的方法，用于对多个因素进行影响量分解。该方法能综合考虑不同因素之间的相互作用，使得单因素影响量之和恰好等于所有因素综合影响量，并从理论和应用两个方面对此进行论证。总体来看，通过一个时期的探索，“自然–社会”二元水循环的耦合研究已形成了诸多的知识点[84] 。
3. 主要存在的科学技术问题
3.1. 降雨径流的非稳定性关系
目前几乎所有模型的应用都是建立在稳定降雨径流关系假设的基础上，针对非稳定的降雨径流关系，建立和应用水文模型来模拟降雨径流关系的专门研究还不够深入。已有的研究在处理非稳定方面存在很大的缺陷，大多数都是针对其中的一个方面，缺乏系统性的理论和技术方法研究。因此，变化环境下的降雨径流非稳定性关系模拟所面临问题和挑战很多，包括模型参数和结构随时间(变化的气候)转移的问题；模型参数和结构随空间(不同的流域)转移的问题；上述两种情况的组合，即变化环境下无资料流域水文响应的模拟；水文模型和气候模型耦合的可能性等。环境变化下带来的这些问题，导致了现有水文模型存在着重大的理论缺陷，如果不及时解决和突破，水文模型应用结果的可靠性和精确性就值得商榷和怀疑。
国内外许多水文学家对如何甄别水文序列的变化和探讨其变化原因做了大量的研究工作，但对于如何甄别各水文序列相互关系的稳定性，分析各水文过程之间关系的稳定性以及非稳定原因和模拟等方面，还缺少成熟的理论和方法[85] 。
3.2. 水文模型参数的时空尺度转化
变化环境下水文模型参数转移问题到现在还没有很好的解决方法，主要表现为不同的研究不能得出统一的结论，有时甚至是相反的结论。原因是因为所用的模型及流域的类型和数量各不相同，大部分研究都表示在该领域还需要做大量深入的研究。特别是模型参数和结构同时随时间(不同的时期和时间尺度)和空间(不同的区域和空间尺度)转移的问题。该问题已成为当前和今后水文模型研究的关键和前沿问题之一，迫切需要在充分论证和比较已有方法的基础上，创建一些新的理论和方法，该领域的研究应该引起全世界水文学家的高度重视[86] [87] 。
3.3. 水文模型与气候模型的耦合
国内外的水文气象学家为水文气候模型的耦合做了大量的研究工作并取得了一些进展，但真正成功耦合(双向耦合并完全按照气候模型的时间和空间尺度运行)仍局限于气候模型与陆面过程模型的耦合。而陆面过程模型并不是传统意义上的水文模型，它本来就是气候模型的一部分或者专门为气候模型开发的以能量平衡为核心的垂直过程模式。传统的水文模型以水量平衡为基础、以土壤含水量的时空变化为核心，具有较高的模拟精度。但传统水文模型的基本结构、时空尺度和无能量平衡模块等特点决定了其在近期内还不能与气候模型双向耦合。该领域的研究重点和难点在于除了为水文模型增加能量平衡模块外，更重要的是研究传统水文模型的结构和参数在空间尺度(由流域到区域甚至全球)和时间尺度(由日到小时再到15分钟(气候模型时间步长))的转换可能性和转换规律。在变化环境下的水文水资源影响研究中，水文模型参数随时间步长的变化以及如何变化的问题，一直没有得到重视、研究和解决。针对不同时间步长，水文模型结构和参数如何演绎需要进一步深入研究。Littlewood等[88] 针对水文模型的参数随时间步长变化的问题，做了些探讨和研究，结论是水文模型参数的率定值随模拟时间步长变化显著。然而目前大部分的水文模型都没有进行水文模型结构和参数随时间步长变化的比较研究。
3.4. 水文模型与水资源优化配置模型的耦合
国内外在水文模拟预测和水资源配置模型进行了较为深入的研究，但由于变化环境下水文机理的复杂性、流域的多样性、资料的限制和模型的局限性等种种原因，许多科学和技术难题并未解决，尚未将水资源系统动态演化和系统内各单元之间相互联系之间进行耦合，也未考虑动态演化条件下优化配置模型求解算法。因此，十分有必要开展降雨径流非稳定性关系的模拟理论与方法研究，建立一套客观公正的方法体系来应对变化环境下的水文模拟预报、水资源供需动态模拟和水资源优化配置研究，这不仅会丰富发展水文水资源学科理论，而且也是研究变化环境下水资源系统适应性及其对策的基础[89] [90] 。
3.5. 水文水资源与社会耦合系统
目前的研究对于“自然–社会”耦合系统的内在联系、驱动机制、演进规律和反馈机制的系统认知还没有形成，亟待加大研究和创新力度，建立和完善“自然–社会”相互作用的二元水循环科学体系，通过研究水资源随人类社会经济系统发展而演化的科学机理，实现对社会水循环系统有效科学的调控和水资源可持续利用，保障经济社会的可持续发展。
4. 研究内容以及关键科学问题
4.1. 变化环境下非稳定性降雨径流关系及识别方法
不同历史时期的人类活动和气候变化是导致降雨径流关系非稳定性的主要原因。而水文模型通过确定参数来模拟降雨径流关系，通过分析降雨径流序列的变化同模型参数的变化是否存在相关，通过模型参数的变化来反映降雨径流关系非稳定性的变化规律，同时通过模型参数和变化环境的关系来诊断降雨径流关系非稳定性的成因，建立由水文模型参数的稳定性来识别降雨径流关系稳定性的方法。但是，由于水文机理的复杂性、流域的多样性、资料的限制和模型的局限性等原因，模型参数的时间和空间移用尚有很多未解决的难题。现有模型的率定，流量资料是唯一的标准，但是在模型中大部分参数同流量资料并无直接关系，用流量资料率定模型无法确保所有的参数的合理性。通过分析模型“异参同效”的原因、模型参数间的相互关系，以及模型参数和下垫面的关系，分析得到不同模型的参数在时间和空间的演绎规律及成因。
4.2. 考虑时间变化的降雨径流产汇流机理
流域的产流，主要决定于非饱和带地下水运动的机理、特性与运动规律。现有水文模型的基本理论依据是将非饱和带地下水运动基本微分方程组中的连续方程简化为水量平衡方程，仅仅考虑了质量守恒而未考虑能量守恒，从数学上讲，是不能求解的，但是该法将时间视为常数，从静态去研究流域产流，即不考虑时间变量，因而可以求解。考虑时间变化的降雨径流产汇流计算由非饱和带地下水运动基本微分方程推导流域产流方程，导出地表、地下径流产流时刻的计算公式，并以此建立一个将下渗能力曲线、不稳定雨强(随时间而变的)和前期土壤含水量有机地联系起来的流域产流计算方法，从而解决利用下渗能力曲线计算流域产流量的途径。
4.3. 基于时间变化产汇流机理的分布式水文模型
变化环境影响研究中，不同的模型存在时空尺度差异，水文模型时间步长较大，气候模型时间步长小。现有的水文模型是静态产流过程，无法满足小步长的计算要求。因此，探讨水文模型产汇流机制的改进，研究水文模型结构和参数随时间步长的变化规律，以满足所建立的水文模型能适应不同时间步长的计算要求。并同区域气候模型进行耦合，分析变化环境下的水文过程响应。
4.4. 基于水联网技术的水资源优化配置模型
水联网技术包括物理水网(现实的河湖连通及供用水通道系统)、虚拟水网(物理水循环通路及其边界的信息化表达)和市场水网(水资源供需的市场信息、优化调配机制及交互反馈)。物理水网以水资源分区为基础，根据蓄水、引水、提水和调水工程等之间的水力联系，虚拟水网是构建面向对象的水循环通路及其边界的信息模型和信息服务，市场水网是在水资源需求量预测模型的基础上，建立不同部门、不同区域、不同时段用水主体之间相互联系的定量表示，并与虚拟水网、物理水网相耦合，在不同层次、不同区域上建立水资源的供用水关系模型，从而构建出既考虑当代的发展，又要照顾到后代可持续发展的需要；既要照顾到发达地区的发展现实，又要求发达地区的今后发展不应以损害欠发达地区的可持续发展能力为代价；既要追求以提高水资源总体配置效率为中心的优化配置模式，又要注意效益在全体社会成员之间的公平分配的水资源配置模型。
4.5. 流域水文水资源与社会系统耦合概念模型及应用
流域水文水资源与社会系统耦合概念模型能对水文气象信息、土地利用变化及各种水资源信息等在流域进行空间分析管理，模拟变化环境下水文响应过程，从流域/区域角度研究水文气象变化、土地利用格局变化、水资源开发利用等过程对流域水文过程的综合影响。流域水文水资源与社会系统可实现统筹考虑水文水资源、宏观经济与生态环境的流域水资源综合管理分析功能的目标。不同的模型存在时空尺度差异，如水文模型结构和参数的时间特性，能适应气候模型或者气象模型时间步长小的特点，同样的问题也存在于多目标决策模型与水资源优化配置模型之间、分布式水文模型与水资源优化配置模型之间。因此流域水文水资源与社会系统耦合概念模型模拟的关键问题为以基于可变下渗能力曲线的分布式水文模型为核心嵌套模型之间的尺度转换问题。
5. 结论
以国际水文科学协会Panta Rhei (2013~2022)科学计划中提出的科学问题为研究方向，结合我国实行最严格水资源管理制度的试点流域为研究对象，开展水文水资源与社会耦合系统研究，分析探讨流域水文模拟预测和水资源优化配置与管理中存在的科学技术问题，研究提出水文水资源与社会耦合系统的新理论和新方法，为流域防洪减灾及经济社会可持续发展提供技术支撑，为国际水文科学计划做出我国应有的贡献，对提高我国水文水资源学科的水平和国际影响力，具有重大的理论价值和现实意义。
探索并揭示流域水文水资源与社会耦合系统的外部驱动因素和内部变化规律以及相互作用和反馈机制。具体的科学目标包括：①识别降雨径流的非稳定性关系，分析现有水文模型随时间和空间变化的演绎规律及内在成因；②通过对非饱和带地下水运动的机理与特性分析，建立能反映时间特性的降雨径流产汇流计算方程，实现变化环境下的水文过程响应；③研究和实现考虑了大气–植被–土壤之间物理交换过程的分布式水文模型，反映土壤、植被、大气中水热状态变化和水热传输；④考虑社会经济发展变化，建立能适应社会经济动态演化的水资源优化配置模型；⑤建立流域水文水资源与社会系统耦合概念模型，以流域为目标区域，实现水资源合理配置和科学管理的目标。
基金项目
国家自然科学基金重点项目(51539009)资助。