基于WetSpa分布式水文模型的洪水多维度研究：以南方典型流域为例

基于WetSpa分布式水文模型的洪水多维度研究：以南方典型流域为例一、引言1.1研究背景与意义洪水灾害作为一种极具破坏力的自然灾害，始终对人类的生存与发展构成严重威胁。在全球范围内，洪水的发生频率与影响范围呈上升趋势，给社会经济、生态环境以及人类生活带来了巨大冲击。20世纪90年代，我国洪灾造成的直接经济损失约12000亿元人民币，仅1998年就高达2600亿元人民币，水灾损失占国民生产总值（GNP）的比例在1%-4%之间，远高于美国、日本等发达国家。洪水不仅会导致人员伤亡，还会对水利、交通、电力、通信等基础设施造成严重破坏，使农田、城镇、村庄大面积受淹，农作物减产甚至绝收，耕地失去耕作条件，企业停工停产，正常的生活生产秩序被打乱。同时，洪水还可能引发河流改道、生态破坏、水环境污染等问题，进一步危及人民群众的健康和生存环境。面对洪水灾害带来的严峻挑战，深入开展洪水研究显得尤为重要。准确模拟洪水过程，深入分析土地利用变化对洪水的影响，以及科学推求设计洪水，对于制定有效的防洪减灾措施、保障人民生命财产安全、促进社会经济可持续发展具有不可替代的作用。随着计算机技术和“3S”（GIS、RS、GPS）技术的飞速发展，分布式水文模型应运而生，为洪水研究提供了更为先进和有效的手段。WetSpa分布式水文模型作为一种基于物理过程的模型，能够充分考虑流域下垫面条件的空间变异性，真实地反映降雨径流的形成机理和产汇流过程。它通过对流域内土壤、植被、大气间的水汽传输及能量交换进行细致模拟，能够更准确地预测洪水的发生和发展。将WetSpa模型应用于洪水研究，可以为洪水管理、土地合理规划以及水利工程科学管理提供坚实的理论依据和技术支持，从而有效提高防洪减灾能力，减少洪水灾害带来的损失。1.2国内外研究现状在国外，WetSpa模型的应用和研究起步较早，已在多个流域开展了相关研究，并取得了一定成果。在欧洲的一些流域，研究者利用WetSpa模型对不同土地利用类型下的水文响应进行了模拟分析，结果表明该模型能够较好地捕捉土地利用变化对径流的影响，为流域土地利用规划和水资源管理提供了科学依据。在北美洲的部分地区，研究人员运用WetSpa模型进行洪水模拟和预测，通过与实测数据的对比验证，发现该模型在洪水过程模拟方面具有较高的精度，能够为防洪决策提供可靠的支持。在国内，WetSpa模型的应用研究也逐渐受到关注。舒晓娟等人引进WetSpa模型，以广东省流溪河水库流域为研究对象，利用流域的DEM、土地利用、土壤等数字化资料，建立了流溪河水库流域的WetSpa模型框架。通过对1980-2004年流域出口断面的流量资料进行分析，选取20场有代表性的大中小洪水，其中15场用于率定参数，5场用于验证模拟结果，研究结果表明WetSpa模型适用于流溪河水库入库洪水的模拟。此外，还有学者利用WetSpa模型研究了广东省流溪河水库流域8种土地利用情景下的2次典型洪水过程，定量分析了该流域植树造林的洪水响应，发现将灌木改造成林地和退耕还林都能起到减少洪量和削减洪峰的作用，且土地利用/植被覆盖变化对小洪水的影响比大洪水更大。尽管国内外在利用WetSpa模型研究洪水方面取得了一定进展，但仍存在一些空白与不足。一方面，不同地区的下垫面条件和气候特征差异较大，WetSpa模型在不同区域的适用性和模拟精度还需要进一步验证和优化。另一方面，目前对于土地利用变化与洪水响应之间的复杂关系研究还不够深入，缺乏全面系统的认识。此外，在设计洪水推求方面，虽然基于WetSpa模型的方法已经得到应用，但与传统方法的对比分析还不够充分，如何结合多种方法提高设计洪水推求的准确性和可靠性，仍有待进一步研究。1.3研究内容与方法本研究基于WetSpa分布式水文模型，主要开展以下几个方面的研究：洪水模拟研究：收集目标流域的地形、土地利用、土壤、气象等数据，利用WetSpa模型构建流域水文模型框架。通过对历史洪水事件的模拟，率定和验证模型参数，评估模型在研究流域的适用性和模拟精度，分析洪水的产汇流过程和时空分布特征。土地利用变化对洪水的响应研究：根据研究流域的土地利用现状和未来发展趋势，设定不同的土地利用情景。运用WetSpa模型模拟不同情景下的洪水过程，定量分析土地利用变化对洪峰流量、洪量、峰现时间等洪水特征参数的影响，揭示土地利用变化与洪水响应之间的内在关系。设计洪水推求研究：分别采用基于流量资料的方法、传统的推理公式法以及基于WetSpa模型的方法推求设计洪水洪峰流量。对三种方法计算的结果进行对比分析，探讨不同方法的优缺点和适用条件，为设计洪水的合理推求提供参考依据。在研究方法上，本研究综合运用多种技术手段和方法。利用地理信息系统（GIS）技术对流域数据进行处理和分析，为WetSpa模型的构建提供基础数据支持；采用E-FAST全局参数敏感性分析法对模型参数进行敏感性分析，确定对模型输出影响较大的关键参数；运用PEST独立参数自动率定程序与人工试错法相结合的方式对模型参数进行优化率定，提高模型的模拟精度；通过模型选择准则对模拟结果进行评价和分析，确保研究结果的可靠性和科学性。二、WetSpa分布式水文模型概述2.1WetSpa模型基本原理2.1.1模型结构WetSpa模型是一种基于单元格的分布式物理水文模型，其结构设计基于对流域复杂水文过程的深入理解。该模型将流域离散化为多个大小相同的单元格，每个单元格都被视为一个独立的水文响应单元，通过对这些单元格的模拟来反映整个流域的水文特征。在垂直方向上，流域被概化为大气、植物冠层、植物根系区、水汽传输带（非饱和层）及饱和带5层。这种分层结构能够细致地描述水分在不同层面的运动和转化过程。降雨首先与植物冠层相互作用，部分降雨被冠层拦截，满足拦截容量后降落到地面的雨水被称为穿透雨。穿透雨在地面上会经历不同的过程，一部分在降雨初期填满地面上的洼地后产生地面径流，另一部分则下渗到土壤中。进入土壤的水又依据不同的土壤含水量，或以土壤水的形式存储在根系区，或沿水平方向流动形成壤中流。土壤水向下运动被简化为一维垂直流，包括地表下渗、深层下渗、非饱和区毛管水及地下水补给。其中，地表下渗是指水分从地表进入土壤表层的过程，其速率受到土壤质地、初始含水量、降雨强度等因素的影响；深层下渗则是土壤水在重力作用下向更深层土壤的运动；非饱和区毛管水的运动则与土壤孔隙结构和水分含量的梯度有关，它使得水分在非饱和土壤中能够进行再分配；当土壤水到达饱和带时，就会形成地下水补给，为地下径流提供水源。在每个单元格中，WetSpa模型通过建立水量平衡方程和能量平衡方程来模拟水文过程。水量平衡方程考虑了降雨、截留、蒸发、下渗、地表径流、壤中流和地下径流等各个分量之间的关系，确保在任何时刻单元格内的水量保持平衡。能量平衡方程则考虑了太阳辐射、感热通量、潜热通量和土壤热通量等能量项，用于计算蒸发和蒸腾等与能量相关的水文过程。通过对这些方程的求解，可以得到每个单元格在不同时刻的水文状态变量，如土壤含水量、地下水位、径流流量等。这种基于物理过程的模拟方式，使得WetSpa模型能够真实地反映流域内的水文循环过程，为洪水研究提供了坚实的理论基础。2.1.2产汇流机制降雨产流：降雨产流是洪水形成的首要环节，WetSpa模型采用超渗产流和蓄满产流相结合的机制来模拟这一过程。当降雨发生时，首先满足植物冠层的截留，截留量与植被类型、叶面积指数等因素密切相关。截留后的降雨到达地面，一部分会填补地面的洼坑，形成填洼量。填洼量的大小取决于地面的粗糙度和地形条件。当降雨强度超过土壤的下渗能力时，就会产生超渗地面径流；而当土壤含水量达到田间持水量后，后续降雨则会形成蓄满地面径流。土壤下渗能力是一个关键参数，它受到土壤质地、结构、初始含水量以及降雨特性等多种因素的影响。在WetSpa模型中，通常采用霍顿下渗公式或格林-安普特下渗公式来计算下渗率。霍顿下渗公式假设下渗率随时间呈指数衰减，即f=f_c+(f_0-f_c)e^{-kt}，其中f为t时刻的下渗率，f_0为初始下渗率，f_c为稳定下渗率，k为衰减系数。格林-安普特下渗公式则基于土壤水动力学原理，考虑了土壤吸力和饱和导水率等因素，能更准确地描述下渗过程，但计算相对复杂。在实际应用中，根据流域的具体情况选择合适的下渗公式。对于土壤质地较为均匀、下渗特性变化较小的流域，霍顿下渗公式可能更为适用；而对于土壤特性空间变异性较大的流域，格林-安普特下渗公式能更好地反映下渗过程的复杂性。通过准确模拟降雨产流过程，WetSpa模型可以为后续的汇流计算提供可靠的输入。坡面流与河道汇流演算：坡面流是指降雨产生的地表径流在坡面上的流动过程。在WetSpa模型中，坡面流采用运动波方程进行演算。运动波方程基于水流的连续性方程和动量方程推导而来，它假设水流的惯性力和重力相比可以忽略不计，主要考虑了水流的沿程阻力。通过求解运动波方程，可以得到坡面流的流速和流量随时间和空间的变化。坡面流的流速受到坡面坡度、糙率以及水深等因素的影响。一般来说，坡面坡度越大，流速越快；糙率越大，流速越小。在实际应用中，糙率的取值需要根据坡面的植被覆盖情况、土壤类型等因素进行确定。河道汇流演算是将坡面流汇集到河道后，对河道内水流的运动进行模拟。WetSpa模型采用马斯京根法或扩散波方程进行河道汇流演算。马斯京根法是一种基于线性水库概念的河道汇流方法，它通过建立河段的水量平衡方程和蓄泄关系方程，将上游断面的流量过程转化为下游断面的流量过程。扩散波方程则考虑了水流的惯性力和扩散作用，能更准确地描述河道内水流的非恒定流特性。在选择河道汇流演算方法时，需要综合考虑河道的形态、水力特性以及数据的可获取性等因素。对于河道形态较为规则、水流变化相对平缓的情况，马斯京根法通常能够满足计算精度要求；而对于河道情况复杂、水流变化剧烈的流域，扩散波方程可能更为合适。通过合理的坡面流与河道汇流演算，WetSpa模型能够准确地模拟洪水在流域内的传播过程，为洪水预报和防洪决策提供重要的依据。2.2WetSpa模型特点与优势WetSpa模型具有诸多显著特点与优势，使其在洪水研究领域具有重要的应用价值。首先，该模型能够充分反映流域地理特征的空间变异性。通过将流域划分为多个单元格，每个单元格可以赋予不同的地形、土壤、土地利用等属性，从而精确地考虑流域下垫面条件的差异对水文过程的影响。在地形复杂的山区，不同坡度和坡向的区域会导致降雨截留、下渗和坡面流等过程的不同，WetSpa模型能够通过对单元格属性的设置，准确地模拟这些空间变化，而传统的集总式水文模型则难以做到这一点。其次，WetSpa模型在模拟洪水要素的空间分布方面表现出色。它可以输出每个单元格的径流、土壤含水量等信息，从而直观地展示洪水在流域内的发生和发展过程。研究人员可以通过这些详细的空间分布数据，分析洪水的传播路径、淹没范围以及不同区域的洪水风险程度，为防洪减灾措施的制定提供科学依据。例如，在城市规划中，可以根据WetSpa模型模拟的洪水风险图，合理布局建筑物和基础设施，避免在高风险区域进行建设，从而降低洪水灾害带来的损失。再者，WetSpa模型基于物理过程的模拟方式，使其具有坚实的理论基础。它考虑了土壤、植被、大气间的水汽传输及能量交换等复杂过程，能够真实地反映自然界降雨径流的形成机理。与一些基于经验或统计关系的模型相比，WetSpa模型在不同的气候和下垫面条件下具有更好的通用性和适应性。在不同的气候区域，无论是湿润地区还是干旱地区，WetSpa模型都可以通过合理设置参数，准确地模拟当地的水文过程，为水资源管理和洪水研究提供可靠的工具。此外，WetSpa模型与地理信息系统（GIS）技术的紧密结合，进一步增强了其数据处理和分析能力。通过GIS技术，可以方便地获取和处理流域的地形、土地利用、土壤等数字化资料，并将这些数据导入到WetSpa模型中，实现模型的快速构建和参数化。同时，利用GIS的空间分析功能，可以对模型的模拟结果进行可视化展示和分析，使研究人员能够更直观地理解水文过程的空间分布特征和变化规律。通过制作洪水过程的动态地图，能够清晰地展示洪水在不同时刻的演进情况，为防洪指挥决策提供有力的支持。2.3WetSpa模型应用范围与局限性WetSpa模型在水文研究领域具有广泛的应用范围，但也存在一定的局限性。从应用范围来看，WetSpa模型适用于多种类型的流域，尤其是地形、土壤和土地利用等下垫面条件空间变异性较大的流域。在山区流域，由于地形起伏较大，不同区域的水文过程差异明显，WetSpa模型能够通过对单元格的细致划分和属性设置，准确地模拟降雨径流过程，为山区的水资源管理和洪水防治提供重要依据。在城市化程度较高的流域，土地利用类型复杂多样，人类活动对水文循环的影响显著，WetSpa模型可以考虑这些因素，分析城市化对洪水的影响，为城市防洪规划提供科学指导。此外，WetSpa模型还适用于对流域水资源进行长期模拟和预测，以及研究气候变化对水文过程的影响等方面。通过设置不同的气候变化情景，模拟未来气候条件下的水文过程变化，为应对气候变化提供决策支持。然而，WetSpa模型在应用过程中也面临一些局限性。首先，该模型对数据的要求较高，需要大量的地形、土壤、土地利用、气象等数据来准确描述流域的下垫面条件和气象特征。在实际应用中，获取这些数据可能存在困难，尤其是在数据稀缺的地区，数据的质量和完整性也可能影响模型的模拟精度。在一些偏远山区或经济欠发达地区，可能缺乏长期的气象观测数据和详细的土壤调查资料，这会给WetSpa模型的应用带来一定的阻碍。其次，WetSpa模型的参数众多，参数率定是一个复杂且耗时的过程。不同的参数对模型输出的影响程度不同，需要通过敏感性分析来确定关键参数，并采用合适的方法进行率定。参数率定过程需要大量的实测数据进行验证和校准，而实测数据的不足或误差也会影响参数率定的准确性，进而影响模型的模拟效果。在一些流域，由于缺乏足够的实测流量数据，难以对模型参数进行有效的率定，导致模型的模拟精度难以保证。此外，WetSpa模型在模拟一些特殊的水文过程时可能存在一定的局限性。在模拟岩溶地区的水文过程时，由于岩溶地貌的特殊性，存在大量的地下溶洞和管道，水流运动复杂，传统的WetSpa模型难以准确描述这些特殊的水文现象，需要对模型进行改进或结合其他专门的岩溶水文模型进行模拟。三、基于WetSpa模型的洪水模拟案例分析3.1研究区域选择与数据收集3.1.1研究区域概况本研究选取广东省流溪河水库流域作为主要研究区域。流溪河水库位于广东省中部从化市东北，因建于流溪河上游而得名，建成于1958年。该流域控制面积达539平方千米，库容为3.76亿立方米，大坝高78米，长255.5米，水电站装机容量4.2万千瓦，与神岗大凹坝配合，可灌溉农田45万亩，在防洪、灌溉、发电、水产养殖等方面发挥着重要作用。流溪河水库流域属于亚热带季风气候，气候温暖湿润，雨量充沛。多年平均降雨量约为2000毫米，降雨主要集中在4-9月，约占全年降雨量的80%以上，这期间极易发生洪水灾害。流域内地形以山地和丘陵为主，地势北高南低，海拔高度在50-1200米之间，地形起伏较大，坡度变化明显。这种地形条件导致流域内的水文过程复杂，不同区域的降雨截留、下渗和坡面流等过程存在显著差异。流域内的植被类型丰富，主要包括常绿阔叶林、针叶林、针阔混交林以及灌木林等。植被覆盖率较高，达到70%以上，植被的存在对流域的水文循环有着重要影响，它不仅可以截留降雨，减少地表径流的产生，还能增加土壤的入渗能力，调节土壤水分的动态变化。此外，流域内的土壤类型多样，主要有红壤、黄壤、赤红壤等，土壤质地和结构的差异也会影响水分的下渗和储存。3.1.2数据收集与预处理DEM数据：从相关测绘部门获取了研究流域1:5万比例尺的数字高程模型（DEM）数据。该数据精度较高，能够准确反映流域的地形地貌特征。在获取数据后，首先利用ArcGIS软件对DEM数据进行预处理，包括去除噪声、填补空洞、平滑处理等操作，以提高数据的质量。通过这些处理，确保DEM数据能够准确地用于后续的流域水系提取、地形分析以及WetSpa模型的构建。利用DEM数据提取流域的水系网络，确定河流的流向和汇流关系，为洪水模拟中的汇流演算提供基础数据。通过对DEM数据的分析，获取流域的坡度、坡向等地形参数，这些参数对于理解流域的水文过程，如坡面流的运动方向和速度等具有重要意义。土地利用数据：土地利用数据来源于最新的遥感影像解译成果。通过对高分辨率遥感影像的解译和分类，得到了研究流域的土地利用类型图，包括耕地、林地、草地、水域、建设用地等。为了使土地利用数据能够与WetSpa模型的要求相匹配，对解译得到的土地利用类型进行了重新编码和分类。将一些相似的土地利用类型进行合并，确保模型能够准确识别和处理不同的土地利用条件。同时，对土地利用数据进行了空间校正和投影转换，使其与DEM数据和其他数据具有相同的坐标系统，以便于数据的集成和分析。土壤数据：土壤数据主要包括土壤质地、土壤孔隙度、土壤饱和导水率等信息。这些数据通过收集相关的土壤调查资料和野外采样分析获得。在收集到土壤数据后，根据土壤类型和空间分布特征，将其转换为WetSpa模型所需的参数格式。利用地理信息系统（GIS）的空间分析功能，将土壤数据与DEM数据和土地利用数据进行叠加分析，确定不同区域的土壤参数，为模型模拟土壤水分运动和下渗过程提供准确的数据支持。降雨径流数据：降雨数据来自流域内及周边的多个雨量站的实测记录，时间跨度为1980-2004年。在收集降雨数据时，对数据的完整性和准确性进行了严格的检查和验证，剔除了明显错误和缺失的数据。对于部分缺失的数据，采用空间插值方法进行补充，确保降雨数据能够完整地反映流域内的降雨时空分布特征。径流数据则来源于流域出口断面的实测流量资料，同样对其进行了质量控制和预处理。为了满足WetSpa模型的模拟需求，将降雨和径流数据按照一定的时间步长进行整理和统计，生成模型所需的输入文件。在数据预处理过程中，还对降雨和径流数据进行了一致性检验，确保数据在时间和空间上的连续性和合理性。3.2WetSpa模型构建与参数率定3.2.1模型构建流程本研究以ArcView3.2为操作平台，基于数字化资料构建WetSpa模型框架，具体步骤如下：流域离散化：利用DEM数据，通过ArcView的水文分析模块，将流溪河水库流域划分为多个大小相同的单元格。根据流域的地形复杂程度和研究精度要求，确定单元格的大小为500米×500米。这种尺度的划分既能较好地反映流域下垫面条件的空间变异性，又能在保证模拟精度的前提下，控制计算量在可接受范围内。在划分单元格的过程中，充分考虑了地形的起伏变化，确保每个单元格都能准确地代表其所在区域的地形特征。对于地形变化剧烈的区域，适当减小单元格的尺寸，以提高模拟的精度；而在地形相对平坦的区域，则适当增大单元格的尺寸，以减少数据量和计算时间。参数赋值：根据收集到的土地利用、土壤等数据，为每个单元格赋予相应的参数值。对于土地利用数据，根据不同的土地利用类型，确定其对应的植被截留参数、糙率参数等。林地的植被截留能力较强，糙率较大，而耕地的植被截留能力相对较弱，糙率较小。对于土壤数据，根据土壤质地、孔隙度等信息，确定土壤的下渗参数、饱和导水率等。砂土的下渗能力较强，饱和导水率较大，而黏土的下渗能力较弱，饱和导水率较小。通过准确地为每个单元格赋值，使得模型能够真实地反映流域内不同区域的水文特性。边界条件设定：明确流域的边界条件，包括流域的进口和出口。在进口处，根据实测的降雨数据和气象条件，设定降雨输入和蒸发边界条件。考虑到流域内可能存在的水汽输送和降水的时空变化，采用动态的降雨输入边界条件，以更准确地模拟实际的降雨过程。在出口处，根据实测的径流数据，设定流量输出边界条件。确保模型能够准确地模拟流域内的水量平衡和径流过程。同时，还考虑了流域与周边地区的水分交换和能量交换，设定了相应的边界条件，以提高模型的模拟精度。模型集成与调试：将上述步骤得到的单元格参数、边界条件等信息集成到WetSpa模型中，构建完整的模型框架。在模型集成过程中，仔细检查各个参数和边界条件的设置是否正确，确保模型的逻辑一致性。完成模型构建后，进行初步的调试和运行，检查模型是否能够正常运行，输出的结果是否合理。通过对模型的调试，及时发现并解决可能存在的问题，如参数设置不合理、数据格式不匹配等，为后续的参数率定和模拟分析奠定基础。3.2.2参数敏感性分析采用E-FAST（ExtendedFourierAmplitudeSensitivityTest）全局参数敏感性分析法对WetSpa模型的参数进行敏感性分析。E-FAST方法是一种基于傅里叶变换的全局敏感性分析方法，它能够同时考虑多个参数之间的相互作用，全面地评估参数对模型输出的影响。在进行敏感性分析时，首先确定需要分析的参数范围和变化步长。根据WetSpa模型的特点和研究区域的实际情况，选择了包括土壤下渗参数、植被截留参数、坡面糙率参数、河道糙率参数等在内的15个关键参数。对每个参数在其合理的取值范围内进行均匀采样，生成一系列的参数组合。将这些参数组合输入到WetSpa模型中进行模拟，得到相应的模型输出结果，如径流流量、土壤含水量等。通过对模型输出结果的分析，计算每个参数的敏感性指数，敏感性指数越大，说明该参数对模型输出的影响越大。根据敏感性分析的结果，确定了对模型输出影响较大的关键参数，包括土壤饱和导水率、植被最大截留量、坡面糙率等。这些关键参数的准确确定对于提高模型的模拟精度至关重要，在后续的参数率定过程中，将重点对这些参数进行优化调整。3.2.3参数率定方法与结果采用人工试错法与PEST（ParameterESTimation）程序联合对WetSpa模型进行参数率定。人工试错法是一种基于经验和对模型原理理解的传统参数率定方法，它通过不断地调整参数值，观察模型模拟结果与实测数据的差异，逐步逼近最优参数值。虽然人工试错法主观性较强，且对于复杂模型和大量参数时效率较低，但它能够充分利用率定者的经验和知识，对模型参数进行初步的筛选和调整。PEST是一种独立于模型的非线性参数估计和不确定性分析的综合软件，它通过优化算法，能够在尽量少的模型运行次数下估计参数。在参数率定过程中，首先利用人工试错法对模型参数进行初步调整，使模型模拟结果与实测数据在趋势上基本一致。然后，将初步调整后的参数作为PEST程序的初始值，利用PEST程序进行自动优化率定。PEST程序通过不断地迭代计算，调整参数值，使得模型模拟结果与实测数据的误差最小化。在对白盆珠水库流域1992-2001年日径流进行连续模拟研究时，经过多次人工试错和PEST程序的优化，最终得到了一组较为满意的参数值。将率定后的参数代入WetSpa模型进行模拟，结果表明，模型模拟的日径流过程与实测数据具有较好的一致性，相关系数达到0.85以上，确定性系数也在0.75以上，说明通过人工试错法与PEST程序联合率定得到的参数能够较好地反映流域的水文特性，提高了模型的模拟精度。3.3洪水模拟结果与验证3.3.1模拟结果展示次洪模拟结果：利用率定后的WetSpa模型对研究流域的不同场次次洪进行模拟，得到了各场次洪水的流量过程线。以1998年的一场典型洪水为例，模拟结果显示，模型能够较好地捕捉到洪水的涨水和落水过程。在涨水阶段，模型模拟的流量迅速上升，与实测数据的变化趋势基本一致，准确地反映了降雨产流后地表径流快速汇集的过程。在落水阶段，模拟流量逐渐下降，也与实测数据相符，体现了洪水消退的过程。通过对比模拟流量过程线和实测流量过程线，可以直观地看到模型在模拟次洪过程中的准确性和可靠性。连续径流模拟结果：对流域的连续径流进行模拟，得到了长期的径流变化过程。以1992-2001年的连续日径流模拟结果为例，模型模拟的径流过程能够较好地反映出流域内的径流年内变化和年际变化特征。在年内变化方面，模拟径流在雨季（4-9月）明显增加，与实际降雨情况相吻合，体现了降雨对径流的直接影响。在年际变化方面，模拟径流的总体趋势与实测数据一致，能够反映出不同年份间气候条件和下垫面条件对径流的综合影响。通过连续径流模拟结果，可以更全面地了解流域的水文循环特征和水资源变化规律。3.3.2模拟结果验证与评价验证指标选择：为了准确验证模拟结果的准确性，选用了多种评价指标，包括相关系数（R）、确定性系数（DC）、洪峰相对误差（REP）、洪量相对误差（REE）等。相关系数用于衡量模拟值与实测值之间的线性相关程度，其值越接近1，说明两者的相关性越好；确定性系数反映了模型模拟结果对实测数据的解释能力，取值范围在0-1之间，值越大表示模型的模拟效果越好；洪峰相对误差和洪量相对误差分别用于评估模拟洪峰流量和洪量与实测值的相对偏差，误差越小说明模拟结果越准确。不同流域模拟结果评价：在流溪河水库流域，基于E-FAST参数敏感性分析的模拟，除了峰现时间不够理想外，其余检验指标表明洪水模拟的总体效果较好。相关系数达到0.8以上，确定性系数在0.7左右，洪峰相对误差控制在15%以内，洪量相对误差在10%以内，说明模型能够较好地模拟该流域的洪水过程。在祯江流域，WetSpa仅通过一场洪水进行参数率定，另外3场洪水进行检验，除了低水部分拟合不够好外，洪峰、洪量、洪水过程的高水部分拟合较好。相关系数在0.7-0.8之间，确定性系数在0.6-0.7之间，洪峰和洪量相对误差在20%以内，表明模型在该流域也具有一定的适用性。在黄龙带水库流域，虽然只收集到一场数据可靠的降雨洪水资料，但利用该场洪水进行参数率定后，模拟效果总体较好。相关系数和确定性系数均在0.7以上，洪峰和洪量相对误差在15%以内，说明模型能够在数据有限的情况下，对该流域的洪水进行有效的模拟。在白盆珠水库流域，对1992-2001年日径流的连续模拟效果较好，相关系数和确定性系数都在0.8以上，洪峰和洪量相对误差较小，表明模型能够准确地模拟该流域的长期径流变化过程。综合不同流域的模拟结果验证与评价，WetSpa分布式水文模型在我国南方温润地区具有良好的水文拟合性能，能够为洪水研究和防洪决策提供可靠的支持。四、土地利用变化对洪水的影响——基于WetSpa模型的分析4.1土地利用变化情景设定4.1.1研究区域土地利用现状流溪河水库流域当前的土地利用类型丰富多样，主要包括林地、耕地、草地、水域、建设用地等。其中，林地占据了流域面积的较大比例，约为60%，这主要得益于流域内的山地和丘陵地形，适宜树木生长，形成了以常绿阔叶林、针叶林和针阔混交林为主的植被覆盖。林地不仅在保持水土、涵养水源方面发挥着关键作用，还对调节气候、维护生物多样性具有重要意义。耕地面积占流域面积的20%左右，主要分布在地势相对平坦的河谷和山间盆地。这些区域土壤肥沃，灌溉条件便利，有利于农作物的种植。耕地主要种植水稻、蔬菜等农作物，为当地居民提供了重要的粮食和农产品来源。草地面积占比约为10%，多分布在山坡和一些土壤条件相对较差的区域。草地的存在对于防止土壤侵蚀、保持土壤肥力具有一定的作用，同时也为畜牧业的发展提供了一定的饲料资源。水域面积约占流域面积的5%，主要包括河流、湖泊和水库等。流溪河作为流域内的主要河流，为周边地区提供了重要的水资源，支持着农业灌溉、居民生活用水和工业用水等。流溪河水库则在防洪、灌溉、发电等方面发挥着重要的综合效益。建设用地面积占比相对较小，约为5%，主要集中在城镇和乡村地区。随着经济的发展和城市化进程的推进，建设用地面积呈现出逐渐增加的趋势，这对流域的土地利用格局和生态环境产生了一定的影响。通过对土地利用现状的分析，可以看出流域内的土地利用类型分布与地形、土壤、水源等自然条件密切相关。不同土地利用类型的空间分布特征，直接影响着流域的水文循环过程和生态系统功能。在山地和丘陵地区，林地的大面积覆盖有效地减少了地表径流的产生，增加了土壤的入渗能力，从而起到了涵养水源、调节径流的作用。而在河谷和平原地区，耕地和建设用地的分布则对地表径流和地下水的补给产生了不同程度的影响。4.1.2未来土地利用变化情景设定根据流溪河水库流域的发展趋势和限制条件，设定了以下三种土地利用变化情景：自然增长情景：假设在未来一段时间内，流域的土地利用变化主要遵循自然发展的趋势，不考虑大规模的人为干预。在这种情景下，随着城市化进程的持续推进，建设用地将逐渐增加，主要通过占用部分耕地和少量林地来实现扩张。预计到2030年，建设用地面积将增加5%左右，而耕地面积将相应减少3%，林地面积减少1%。同时，由于生态保护意识的逐渐提高，部分退化的草地可能会得到一定程度的恢复和改良，草地面积将略有增加，约增长1%。水域面积则基本保持稳定，仅因自然因素的影响可能会有微小的波动。生态保护情景：该情景以加强生态保护为主要目标，加大对林地和草地的保护力度，严格限制建设用地的扩张，并积极开展退耕还林、还草等生态工程。预计到2030年，林地面积将增加10%，主要通过将部分坡度较大的耕地和退化的草地进行还林还草来实现。耕地面积将减少8%，其中大部分转化为林地，少量转化为草地。建设用地面积将得到严格控制，仅允许在现有城镇和乡村的基础上进行适度的优化和改造，面积基本保持不变。草地面积将增加5%，通过对退化草地的修复和生态治理，提高草地的质量和生态功能。水域面积同样保持相对稳定。经济开发情景：此情景侧重于经济的快速发展，大力推进工业化和城市化进程，加大对土地资源的开发利用力度。预计到2030年，建设用地面积将大幅增加20%，主要通过大规模占用耕地和林地来满足经济发展对土地的需求。耕地面积将减少15%，其中大部分被建设用地占用，少量因生态环境恶化而撂荒。林地面积将减少12%，除了被建设用地占用外，还有部分因过度砍伐和森林火灾等原因而减少。草地面积将减少5%，主要是由于建设用地的扩张和不合理的放牧等活动导致草地退化。水域面积可能会因填湖造地等开发活动而略有减少，约减少1%。4.2不同情景下洪水响应模拟4.2.1模拟方法与过程运用已构建并率定好的WetSpa模型，对上述三种土地利用变化情景下的典型洪水过程进行模拟。首先，根据不同情景设定的土地利用变化情况，对WetSpa模型中的土地利用参数进行相应的调整。在自然增长情景下，按照建设用地增加、耕地和林地减少的比例，重新分配各单元格的土地利用类型，并更新与之相关的植被截留、糙率等参数。在生态保护情景中，依据林地和草地增加、耕地减少的设定，调整模型参数，以反映生态保护措施对土地利用和水文过程的影响。对于经济开发情景，根据建设用地大幅增加、耕地和林地大量减少的情况，对模型参数进行全面更新。在模拟过程中，保持其他输入条件不变，如气象数据（降雨、蒸发等）和地形数据等，以确保模拟结果能够准确反映土地利用变化对洪水响应的影响。将调整后的模型参数输入到WetSpa模型中，以小时为时间步长，对典型洪水过程进行连续模拟。在模拟过程中，模型根据设定的参数和输入数据，计算每个单元格在不同时刻的降雨截留、下渗、地表径流、壤中流和地下径流等水文过程，并通过坡面流和河道汇流演算，得到流域出口断面的洪水流量过程线。在模拟完成后，对模拟结果进行整理和分析，提取洪峰流量、洪量、峰现时间等关键洪水要素，以便进行后续的结果分析。4.2.2模拟结果分析洪峰流量变化：在自然增长情景下，模拟结果显示洪峰流量略有增加，增幅约为5%。这是由于建设用地的增加导致不透水面积扩大，地表径流的产生速度加快，使得洪峰流量有所上升。在生态保护情景中，洪峰流量明显减小，降幅达到15%左右。林地和草地面积的增加，增强了植被对降雨的截留和土壤的入渗能力，有效地削减了地表径流，从而降低了洪峰流量。而在经济开发情景下，洪峰流量大幅增加，增幅高达30%。大量的耕地和林地被建设用地占用，使得流域的下垫面条件发生了显著变化，不透水面积急剧扩大，地表径流迅速汇集，导致洪峰流量急剧上升。洪量变化：自然增长情景下，洪量增加约8%，主要是因为建设用地的扩张减少了土壤的蓄水能力，使得更多的降雨转化为地表径流，从而增加了洪量。生态保护情景中，洪量减少约20%，林地和草地的增加不仅增加了植被截留和土壤入渗，还延长了洪水的汇流时间，使得更多的雨水能够被储存和缓慢释放，从而减少了洪量。在经济开发情景下，洪量大幅增加约40%，建设用地的大量增加破坏了原有的水文循环平衡，导致地表径流大量增加，洪量显著增大。峰现时间变化：自然增长情景下，峰现时间略有提前，提前约1-2小时，这是由于地表径流的增加使得洪水的传播速度加快。生态保护情景中，峰现时间有所推迟，推迟约3-4小时，植被和土壤对洪水的调节作用使得洪水的汇流过程减缓。经济开发情景下，峰现时间明显提前，提前约5-6小时，大量不透水面积的存在使得洪水能够迅速汇集和传播，导致峰现时间大幅提前。综合不同情景下的模拟结果可以看出，土地利用变化对洪水的影响十分显著。生态保护情景下的土地利用方式调整能够有效地削减洪峰、减少洪量并推迟峰现时间，对防洪减灾具有积极作用；而经济开发情景下的大规模土地开发则会加剧洪水灾害的风险，增加洪峰流量、洪量并提前峰现时间。自然增长情景下的土地利用变化对洪水的影响介于两者之间，但也呈现出洪峰和洪量增加、峰现时间提前的趋势，需要引起关注。4.3土地利用变化与洪水关系探讨通过对不同土地利用变化情景下洪水响应的模拟分析，可以总结出土地利用变化对洪水影响的规律。建设用地的增加会导致不透水面积扩大，减少土壤的入渗和蓄水能力，使得地表径流迅速增加，从而增大洪峰流量和洪量，提前峰现时间，加剧洪水灾害的风险。而林地和草地面积的增加，则能够增强植被对降雨的截留和土壤的入渗能力，延长洪水的汇流时间，有效地削减洪峰、减少洪量并推迟峰现时间，对洪水起到调节和缓冲作用，降低洪水灾害的影响。基于上述规律，为了实现合理的土地利用规划，降低洪水灾害风险，提出以下建议：加强生态保护与修复：加大对林地和草地的保护力度，严格限制对森林和草原的破坏。积极开展退耕还林、还草工程，增加植被覆盖面积，提高流域的生态系统服务功能。在山区和丘陵地区，应优先保护和恢复森林植被，充分发挥森林在涵养水源、保持水土、调节径流方面的重要作用。通过植树造林、封山育林等措施，提高森林覆盖率，改善生态环境，从而有效减轻洪水灾害的威胁。优化建设用地布局：在城市化和工业化进程中，应科学合理地规划建设用地，避免盲目扩张。尽量减少对耕地和林地的占用，充分利用闲置土地和低效建设用地进行开发。在城市建设中，推广海绵城市理念，增加城市绿地、湿地和透水地面的面积，提高城市的雨水吸纳和蓄渗能力，减少地表径流的产生，降低城市内涝的风险。合理布局城市排水系统，确保雨水能够及时排出，避免积水造成的灾害。推进农业可持续发展：在耕地利用方面，应推广可持续的农业生产方式，减少对土壤的破坏和污染。采用节水灌溉技术，合理施肥和使用农药，保护土壤的结构和肥力，提高土壤的入渗能力。通过建设梯田、等高种植等措施，减少坡面径流和土壤侵蚀，保护农田生态环境。同时，适度调整农业产业结构，减少高耗水作物的种植面积，发展节水型农业，降低农业用水对水资源的压力，从而减少因农业活动对洪水产生的不利影响。五、基于WetSpa模型的设计洪水推求5.1设计洪水推求方法概述设计洪水推求是水利工程规划设计中的关键环节，其结果直接关系到工程的安全与效益。目前，常见的设计洪水推求方法主要有流量资料法、推理公式法和WetSpa模型法，它们各自基于不同的原理和数据基础。流量资料法是一种较为传统且常用的方法，其原理基于概率论和数理统计。该方法假设洪水系列是独立同分布的随机变量，通过对实测洪水流量资料进行分析，运用频率计算的方法来推求设计洪水。在实际操作中，首先需要收集和整理研究流域的实测洪水流量数据，对这些数据进行“三性”审查，即可靠性、一致性和代表性审查。可靠性审查主要是检查数据的测量和记录是否准确可靠，避免因观测误差或数据处理不当导致的错误；一致性审查则关注洪水资料产生的流域产流和汇流条件在观测期内是否基本相同，若发生较大变化，需将资料还原到天然状态；代表性审查要求实测洪水年数一般不少于20-30年，且有特大洪水加入，以确保样本能够充分反映总体的特征。经过审查后，采用年最大值法从资料中逐年选取一个最大流量和固定时段的最大洪水总量，组成洪峰流量和洪量系列。然后，根据概率论中的频率计算方法，如矩法、适线法等，确定洪水系列的统计参数，如均值、变差系数和偏态系数等，进而推求出不同频率的设计洪峰流量和设计洪量。这种方法适用于有较长系列实测洪水流量资料的流域，其优点是直接利用实测数据，计算结果具有一定的可靠性和直观性。然而，该方法也存在局限性，当实测资料系列较短或代表性不足时，计算结果的不确定性会增加，难以准确反映稀遇洪水的特征。推理公式法是一种基于流域水文参数的经验性洪水推算方法，主要适用于小流域或缺乏实测洪水资料的地区。该方法的原理是利用流域的地形地貌、土壤性质、植被覆盖等特征，结合降雨资料，通过建立经验公式来推求设计洪水。其核心在于将流域的产流和汇流过程进行简化和概化，以流域面积、降雨强度、汇流时间等参数为基础，建立起设计洪峰流量与这些参数之间的关系。推理公式法通常采用的计算公式为Q=C\cdotI\cdotA，其中Q为设计洪峰流量，A为流域面积，C为径流系数，I为设计暴雨强度。在实际应用中，首先需要确定流域的面积、主河道长度、河道平均比降等流域特征参数，以及损失参数、汇流参数等产汇流参数，还有暴雨参数，如t=1小时的平均雨强、暴雨指数等。这些参数的确定通常依赖于流域的地形测量数据、土壤调查资料、降雨观测数据以及相关的水文手册和图集。通过将这些参数代入推理公式，可以计算出设计洪峰流量。推理公式法的优点是计算相对简单，所需数据相对较少，能够在资料有限的情况下对设计洪水进行初步估算。但由于其基于经验公式，对流域的概化程度较高，计算精度相对较低，且对于流域特征和降雨条件的变化较为敏感，在应用时需要结合实际情况进行合理的调整和验证。WetSpa模型法是基于分布式水文模型的设计洪水推求方法，它充分考虑了流域下垫面条件的空间变异性和水文过程的物理机制。WetSpa模型将流域离散化为多个单元格，每个单元格都被视为一个独立的水文响应单元，通过对这些单元格的模拟来反映整个流域的水文特征。在垂直方向上，流域被概化为大气、植物冠层、植物根系区、水汽传输带（非饱和层）及饱和带5层，通过建立水量平衡方程和能量平衡方程来模拟水文过程。在推求设计洪水时，首先利用研究流域的DEM、土地利用、土壤等数字化资料，结合气象数据，构建WetSpa模型框架。然后，通过对模型参数的率定和验证，使模型能够准确地模拟流域的水文过程。在模型运行过程中，根据设定的设计暴雨情景，模型可以模拟出相应的洪水过程，包括地表径流、壤中流和地下径流的产生和汇流过程，从而得到设计洪水的洪峰流量、洪量和洪水过程线等信息。这种方法的优势在于能够更真实地反映流域的实际情况，考虑了地形、土地利用、土壤等因素对洪水的影响，模拟结果具有较高的精度和可靠性。但该方法对数据的要求较高，需要大量的地形、土壤、土地利用、气象等数据来准确描述流域的下垫面条件和气象特征，且模型参数众多，参数率定过程复杂，计算量较大。5.2基于WetSpa模型的设计洪水计算5.2.1模型输入与设置数据输入：地形数据：采用高分辨率的数字高程模型（DEM）数据，其精度直接影响到对流域地形特征的刻画。本研究使用的是1:5万比例尺的DEM数据，通过ArcGIS软件进行预处理，包括去除噪声、填补空洞、平滑处理等操作，以确保数据的准确性。利用DEM数据提取流域的水系网络，确定河流的流向和汇流关系，同时获取流域的坡度、坡向等地形参数。这些地形信息对于理解流域的水文过程至关重要，坡度影响坡面流的速度和方向，坡向则与太阳辐射、降水分布等因素相关，进而影响流域的产汇流过程。土地利用数据：来源于最新的遥感影像解译成果，经过重新编码和分类，将土地利用类型划分为耕地、林地、草地、水域、建设用地等。为使数据与WetSpa模型要求匹配，对其进行空间校正和投影转换，使其与DEM数据和其他数据具有相同的坐标系统。土地利用类型直接影响植被截留、糙率等水文参数，林地的植被截留能力较强，糙率较大，而建设用地的不透水性质会导致地表径流迅速增加。土壤数据：通过收集相关的土壤调查资料和野外采样分析获得，包括土壤质地、土壤孔隙度、土壤饱和导水率等信息。利用GIS的空间分析功能，将土壤数据与DEM数据和土地利用数据进行叠加分析，确定不同区域的土壤参数，为模型模拟土壤水分运动和下渗过程提供准确的数据支持。不同土壤质地的下渗能力和蓄水能力差异显著，砂土的下渗能力较强，而黏土的蓄水能力相对较高。气象数据：包括降雨、蒸发、气温等数据。降雨数据来自流域内及周边多个雨量站的实测记录，时间跨度为1980-2004年，对数据进行完整性和准确性检查，剔除错误和缺失数据，对于缺失数据采用空间插值方法补充。蒸发数据可通过气象站的观测资料或基于能量平衡原理的估算方法获得，气温数据同样来自气象站记录。这些气象数据是驱动WetSpa模型的重要输入，直接影响到降雨产流和蒸发蒸腾等水文过程。参数设置：模型结构参数：根据研究流域的实际情况，确定模型的单元格大小。本研究中单元格大小设定为500米×500米，这一尺度既能较好地反映流域下垫面条件的空间变异性，又能在保证模拟精度的前提下控制计算量。同时，明确流域的边界条件，包括流域的进口和出口，在进口处设定降雨输入和蒸发边界条件，在出口处设定流量输出边界条件。水文过程参数：对于降雨产流过程，设置植物冠层截留参数、地面填洼参数、土壤下渗参数等。植物冠层截留参数根据植被类型和叶面积指数确定，不同植被类型的截留能力不同；地面填洼参数与地面粗糙度和地形条件相关；土壤下渗参数采用霍顿下渗公式或格林-安普特下渗公式进行计算，根据土壤质地、初始含水量等因素确定公式中的参数值。在坡面流和河道汇流演算过程中，设置坡面糙率参数和河道糙率参数，坡面糙率根据坡面的植被覆盖情况和土壤类型确定，河道糙率则与河道的形态、河床材料等因素有关。5.2.2计算过程与结果计算过程：模型初始化：将输入的数据和设置的参数加载到WetSpa模型中，对模型进行初始化，确定模拟的起始时间和结束时间，设置时间步长，本研究中时间步长设定为1小时，以满足对洪水过程动态模拟的需求。降雨产流计算：模型根据输入的降雨数据，首先计算植物冠层截留量，当降雨超过截留容量后，剩余降雨到达地面。地面降雨一部分用于填补洼坑，形成填洼量，当降雨强度超过土壤下渗能力时，产生超渗地面径流；当土壤含水量达到田间持水量后，产生蓄满地面径流。模型根据设定的下渗公式计算下渗率，确定下渗到土壤中的水量，同时考虑土壤水分在根系区的存储和沿水平方向的壤中流运动。坡面流与河道汇流演算：产生的地表径流和壤中流在坡面和河道中进行汇流演算。坡面流采用运动波方程进行演算，根据坡面坡度、糙率和水深等因素计算坡面流的流速和流量。河道汇流演算采用马斯京根法，通过建立河段的水量平衡方程和蓄泄关系方程，将上游断面的流量过程转化为下游断面的流量过程。在汇流演算过程中，模型考虑了河道的糙率、河宽、水深等因素对水流运动的影响。洪水过程模拟：通过不断迭代计算，模型模拟出整个洪水过程中流域内各单元格的水文状态变量，如土壤含水量、地下水位、径流流量等，并最终得到流域出口断面的洪水流量过程线。在模拟过程中，模型实时更新各单元格的水文状态，考虑不同单元格之间的水分交换和能量交换，以准确反映洪水在流域内的传播和演变。计算结果：洪峰流量：通过模型模拟得到不同频率设计洪水的洪峰流量。以某一特定频率（如50年一遇）为例，模拟得到的洪峰流量为[X]立方米每秒。与历史实测洪水资料对比，该洪峰流量在历史洪水序列中处于较高水平，反映了该频率洪水的峰值特征。洪量：计算得到不同时段（如1日、3日、7日等）的设计洪量。以3日设计洪量为例，模拟结果为[Y]立方米，与历史洪水过程中相应时段的洪量对比分析，能够了解不同频率洪水在不同时段内的水量积累情况，为防洪工程的规划和调度提供重要参考。洪水过程线：模型输出的洪水过程线直观地展示了洪水流量随时间的变化情况。通过对洪水过程线的分析，可以清晰地看到洪水的涨水和落水过程，以及峰现时间等关键信息。对于某一设计洪水，模型模拟的洪水过程线显示，洪水在[具体时间]达到洪峰，随后逐渐消退，整个洪水过程持续时间为[Z]小时。这些结果为评估洪水对流域的影响范围和程度提供了重要依据，有助于制定合理的防洪措施。5.3不同方法计算结果比较与分析流量资料法：利用历史实测洪水流量资料，经过“三性”审查和频率计算，得到设计洪峰流量为[流量资料法洪峰流量值]。该方法基于实际观测数据，具有一定的可靠性，但由于实测资料系列的局限性，对于稀遇洪水的推求存在一定的不确定性。当实测资料系列较短时，样本的代表性不足，可能导致频率计算结果偏离真实值，使得设计洪峰流量的估计存在偏差。推理公式法：根据流域的地形地貌、土壤性质、植被覆盖等特征，结合降雨资料，通过推理公式计算得到设计洪峰流量为[推理公式法洪峰流量值]。该方法计算相对简单，适用于小流域或资料缺乏的地区，但由于其基于经验公式，对流域的概化程度较高，计算结果的精度相对较低。在实际应用中，推理公式法对流域特征参数和降雨参数的变化较为敏感，参数的微小误差可能导致计算结果的较大偏差。WetSpa模型法：基于WetSpa模型模拟得到的设计洪峰流量为[WetSpa模型法洪峰流量值]。该方法充分考虑了流域下垫面条件的空间变异性和水文过程的物理机制，能够更真实地反映流域的实际情况，模拟结果具有较高的精度和可靠性。然而，该方法对数据的要求较高，需要大量的地形、土壤、土地利用、气象等数据，且模型参数众多，参数率定过程复杂。在数据获取困难或数据质量不高的情况下，模型的模拟精度可能受到影响。结果比较：将三种方法计算得到的设计洪峰流量进行对比，发现WetSpa模型法计算结果与流量资料法较为接近，这表明WetSpa模型能够较好地模拟流域的洪水过程，其结果具有一定的可靠性。而推理公式法计算结果与其他两种方法存在一定差异，主要原因在于推理公式法对流域的简化和概化程度较高，未能充分考虑流域下垫面条件的空间变异性和水文过程的复杂性。WetSpa模型法优势与不足：优势：考虑空间变异性：能够充分考虑流域地形、土地利用、土壤等下垫面条件的空间变异性，通过对每个单元格的模拟，更准确地反映流域内不同区域的水文过程差异，从而提高设计洪水计算的精度。在地形复杂的山区，不同坡度和坡向的区域会导致降雨截留、下渗和坡面流等过程的不同，WetSpa模型能够准确捕捉这些差异，为洪水模拟提供更详细的信息。物理机制明确：基于物理过程的模拟方式，考虑了土壤、植被、大气间的水汽传输及能量交换等复杂过程，具有坚实的理论基础，能够更真实地反映自然界降雨径流的形成机理，使模拟结果更具可信度。与基于经验或统计关系的模型相比，WetSpa模型在不同的气候和下垫面条件下具有更好的通用性和适应性。洪水要素空间分布展示：可以输出每个单元格的径流、土壤含水量等信息，直观展示洪水在流域内的发生和发展过程，为分析洪水的传播路径、淹没范围以及不同区域的洪水风险程度提供了详细的数据支