保姆级教程：用QSWAT+3.10.6从DEM到出流量曲线，水文模拟避坑指南

零基础实战：QSWAT+3.10.6水文模拟全流程避坑指南

水文模拟是环境工程和水利领域的重要工具，而QSWAT作为QGIS的插件，为初学者提供了相对友好的操作界面。本文将带你从零开始，一步步完成从DEM数据处理到最终流量曲线生成的完整流程，特别针对Windows平台上的常见错误进行预警，确保你的首次水文模拟之旅顺畅无阻。

1. 环境准备与数据获取

1.1 软件安装注意事项

在开始之前，我们需要确保所有必要的软件正确安装。以下是关键组件及其安装要点：

• QGIS 3.10.6：建议从官网下载QGIS-OSGeo4W-3.10.6-1-Setup-x86_64.exe安装包
• SWATPlus Tools：需要单独安装swatplustools-installer-1.2.3.exe
• QSWAT插件：将在QGIS中通过插件管理器安装

注意：所有安装路径必须使用纯英文，避免任何中文字符。这是后续操作能否成功的关键前提。

1.2 基础数据准备

水文模拟需要三类核心数据：

建议提前创建专用项目文件夹，结构如下：

水文项目/ ├── raw_data/ # 存放原始下载数据 ├── processed/ # 存放处理后的中间文件 └── output/ # 存放最终结果

2. DEM数据处理实战

2.1 数据导入与初步检查

打开QGIS后，按以下步骤导入DEM：

1. 点击菜单栏的Layer→Add Layer→Add Raster Layer
2. 浏览到DEM文件位置（通常为.tif格式）
3. 点击Add导入，然后Close对话框

导入后，右键点击图层选择Properties，检查以下信息：

• 像素大小（Pixel Size）是否符合预期
• 坐标系（CRS）是否明确
• 数据范围（Extent）是否覆盖研究区域

2.2 精确裁剪与投影转换

流域模拟通常不需要处理整个DEM文件，裁剪可以显著提升效率：

# 伪代码展示裁剪逻辑 def clip_dem(original_dem, study_area): # study_area应为稍大于目标流域的矩形范围 clipped_dem = original_dem.clip_by_extent(study_area) return clipped_dem

实际操作步骤：

1. 点击Raster→Extraction→Clip raster by extent
2. 在Clipping extent栏手动输入经纬度范围（建议比目标流域大10-20%）
3. 指定输出文件位置（建议保存在processed文件夹）
4. 点击Run执行裁剪

投影转换是关键步骤，常见问题包括：

• 错误1：选择了不合适的投影导致面积计算失真
• 错误2：忘记取消勾选"No projection"选项
• 错误3：输出文件名包含空格或特殊字符

推荐使用UTM投影，具体zone号可通过在线工具查询。转换完成后，务必删除原始图层以保持工程整洁。

3. 土地利用与土壤数据处理

3.1 数据获取与标准化

从SWAT官网获取数据时需注意：

• 土地利用数据（landuse map）和土壤数据（soil map）应覆盖研究区域
• 下载后检查文件完整性，特别是.tif和配套的.dbf文件
• 数据时间范围应与研究时段匹配

3.2 投影统一化处理

所有输入数据必须使用相同的投影系统：

1. 对土地利用数据重复DEM的投影转换步骤
2. 对土壤数据同样处理
3. 使用Raster→Miscellaneous→Build Overviews优化显示性能

提示：可以使用Raster→Alignment→Align Rasters工具确保所有栅格严格对齐，避免后续HRU划分出现问题。

4. QSWAT插件配置与流域划分

4.1 插件安装与项目创建

在QGIS中安装QSWAT插件：

1. 点击Plugins→Manage and Install Plugins
2. 搜索"QSWAT"，勾选QSWATPlus3_64
3. 安装完成后，工具栏会出现"S+"图标

创建新项目时的关键设置：

• 项目路径必须全英文，建议使用简短路径如C:/swat_projects/test1
• 初始设置保持默认，后续可调整
• 项目名称应具有描述性，如yangtze_river_2023

4.2 流域自动划分技巧

使用处理好的DEM进行流域划分：

1. 点击Delineate Watershed按钮
2. 选择处理后的DEM文件
3. 调整Channel/Stream threshold参数（初次尝试可保持默认）
4. 点击Create streams生成河网

常见问题解决方案：

• 河网过于密集：增大threshold值
• 无河网生成：检查DEM是否有洼地，需先填洼
• 出水口定位不准：放大到足够比例尺再标记

# 伪代码展示河网生成算法逻辑 def create_streams(dem, threshold): flow_direction = calculate_flow_direction(dem) flow_accumulation = calculate_flow_accumulation(flow_direction) streams = flow_accumulation > threshold return streams

5. HRU划分与模型运行

5.1 土地利用与土壤数据整合

创建HRUs是模拟精度的关键：

1. 在Create HRUs界面加载处理好的土地利用和土壤数据
2. 选择正确的数据表：
   • Landuse table:global_landuses
   • Soil table:global_soils
   • Usersoil table:global_usersoil
3. 设置坡度带（Slope bands），通常10个足够
4. 勾选Generate FullHRUs shapefile以便后续分析

5.2 气象数据准备

对于没有实测气象数据的情况：

1. 从SWAT+官网下载swatplus_wgn.sqlite全球气象数据库
2. 在Weather Generator中选择wgn_cfsr_world表
3. 根据研究区域经纬度自动匹配最近的气象站点

注意：如果使用本地气象数据，需确保格式符合SWAT+要求，包括日降水量、最高最低温度等必要字段。

6. 模型运行与结果可视化

6.1 参数设置与模型执行

在运行前检查关键设置：

• 模拟时间段是否正确
• 输出频率（日/月/年）是否符合需求
• 必要变量是否已选择输出

运行过程可能耗时较长，取决于：

• 流域面积大小
• HRU数量多少
• 计算机性能

6.2 流量曲线生成与分析

成功运行后，可视化出流量曲线：

1. 在Visualise界面选择channel_sd_day表
2. 筛选目标时间段
3. 变量选择flo_out（出流量）
4. 点击Add plot生成曲线

典型问题排查：

• 无结果显示：检查是否选择了正确的输出单元（subbasin编号）
• 曲线异常平坦：可能气象数据输入有问题
• 负值出现：检查模型参数是否合理

# 示例：使用Matplotlib绘制类似结果 import matplotlib.pyplot as plt import pandas as pd data = pd.read_csv('output/channel_sd_day.csv') plt.plot(data['date'], data['flo_out']) plt.xlabel('Date') plt.ylabel('Flow (m3/s)') plt.title('Daily Streamflow Output') plt.show()

7. 常见错误与解决方案

在实际教学中，我发现以下几个高频错误点值得特别注意：

中文路径问题：尽管反复强调，仍有约40%的初学者会犯此错误。症状包括模型无法运行、结果文件缺失等。解决方法只有一种 - 彻底检查所有文件路径，确保从项目文件夹到每个输入文件都不含中文。

投影不一致：表现为HRU划分失败或流域边界异常。必须确保DEM、土地利用和土壤数据使用完全相同的投影系统，包括相同的中央经线等参数。

气象数据缺失：特别是在中国西部地区，全球气象数据库可能没有足够站点。这时可以考虑使用邻近站点数据，或手动输入本地气象观测。

计算机内存不足：处理大流域时���建议先在小区域测试整套流程，确认无误后再扩展。也可以调整HRU划分参数，减少最终HRU数量。