《ISPRS International Journal of Geo-Information》：A GIS-Based Spatiotemporal Digital Twin-Oriented Framework for a Dammed River Shoreline: Methods, Validation, and Multi-Epoch Analysis
编辑推荐：
数字孪生（Digital Twins）在地理研究中日益被采用为环境系统的动态表示；然而，它们在受调控河流岸线中的应用仍然有限，特别是在必须将测深变化（bathymetric change）、水文变率（hydrological variability）和岸线状态
数字孪生（Digital Twins）在地理研究中日益被采用为环境系统的动态表示；然而，它们在受调控河流岸线中的应用仍然有限，特别是在必须将测深变化（bathymetric change）、水文变率（hydrological variability）和岸线状态动力学（shoreline-state dynamics）整合到一个基于GIS（地理信息系统）的框架中的情况下。研究人员开发并验证了一个基于GIS的时空数字孪生导向框架（GIS-based spatiotemporal digital twin-oriented framework），用于波兰W?oc?awek大坝下游受大坝影响的岸线。该框架将2008–2011年期间采集的四次测深调查（bathymetric surveys）、水位记录、机载激光扫描（airborne laser scanning）数据和三维水工基础设施整合到一个统一的地理数据库（geodatabase）中，旨在实现动态岸线状态重建、多时期分析和环境监测。一个关键的方法要素是将水位作为动态参考面（dynamic reference surface），从而能够为观测和情景水文条件下的淹没区和暴露区进行自动划定。重建的测深表面被组织为一个多维栅格数据集（multidimensional raster dataset），以时间作为显式分析维度，支持可重复的变化检测、横断面解释和趋势分析。为了将框架扩展到纯粹的回顾性分析之外，通过将IMGW（波兰水文与气象研究所）API（应用程序接口）的操作水位观测数据摄入地理数据库，实现了一个近实时水文更新组件（near-real-time hydrological updating component）。基于趋势的2011年预测验证结果：R2 = 0.967，RMSE（均方根误差）= 0.44 m，MAE（平均绝对误差）= 0.28 m，偏差（bias）= ?0.06 m。所提出的框架为在水文条件变化情况下对受调控河流岸线动力学进行时空建模和监测提供了一个可转移的地理空间基础。
## 论文解读：基于GIS的时空数字孪生框架用于筑坝河流岸线动态分析
### 研究背景与问题
数字孪生（Digital Twin, DT）在地理研究中正从静态三维复制转向动态、可更新的分析环境，尤其在环境系统模拟中日益受关注。然而，在受调控河流（regulated river）岸线领域，DT的应用仍十分有限，尤其是在需要将测深变化（bathymetric change）、水文变率（hydrological variability）和岸线状态动力学（shoreline-state dynamics）整合至单一GIS（地理信息系统）框架时。波兰维斯瓦河下游的W?oc?awek大坝（建于1970年）下游河段，受水力调控、沉积物不连续及水工基础设施共同影响，呈现高度非均匀的侵蚀-堆积格局和河床不稳定性。此前研究主要以离散战役式评估（campaign-based assessment）报道，缺乏将测深、水文、地形及基础设施几何形态整合于统一时空分析模型中的系统框架。因此，本研究旨在开发并验证一个基于GIS的时空DT导向框架，实现动态岸线状态重建、多时期分析及验证，并发表在《ISPRS International Journal of Geo-Information》。
### 主要关键技术方法
研究人员采用以下关键方法：
1. **数据整合与统一地理数据库**：集成2008–2011年的四次单波束回声测深调查（bathymetric surveys）、W?oc?awek水电站水位记录、国家机载激光扫描（Airborne Laser Scanning, ALS）数字地形/表面模型（Digital Terrain Model/Surface Model, DTM/DSM）以及三维水工基础设施模型，所有数据统一至PL-2000坐标系（EPSG:2177）和PL-EVRF2007-NH高程基准。
2. **动态水位参考面（Dynamic Water-Level Referencing）**：将水位作为动态参考面，用于自动划分每个测深时期的水面以上（暴露区）和水面以下（淹没区）边界，而非固定高程阈值。
3. **多维栅格结构（Multidimensional Raster Structure）**：重建的测深表面组织为以时间为显式维度的多维栅格数据集（2 m分辨率），支持逐像元比较、趋势分析和短时预测。
4. **近实时水文更新组件（Near-Real-Time Hydrological Updating）**：在ArcGIS Pro 3.7中通过Python 3.13.13和ArcPy 3.7实现，从IMGW（波兰水文与气象研究所）公开API获取W?oc?awek水位站的操作观测数据，经预处理后存储至地理数据库的水文参考状态表（如hydro_current_state），驱动岸线状态重新计算。
5. **测深表面重建与插值选择**：对每个测深时期原始点数据执行Hampel滤波后，采用薄板样条（Thin Plate Spline, TPS）插值（经4000点基准测试，TPS的均方根误差RMSE=0.529 m优于反距离加权IDW，与普通克里金Ordinary Kriging相近）。
6. **趋势预测扩展**：基于2008–2010年三个有序栅格进行逐像元时间趋势外推，生成2011年预测表面，并与实际2011年表面验证。
### 研究结果
**4.1 数字孪生表示与多时期状态组织**：框架将四个测深时期转换为有序空间状态序列，每个时期关联其水位条件，共同存储于统一三维场景。多维栅格结构使各时期测深表面作为独立时间切片，支持直接比较和趋势分析，实现岸线状态基于测深与水位关系而非单纯岸线线提取的动态重建。
**4.2 水动力形态与运行分析支持**
- **4.2.1 侵蚀与沉积空间分布**：通过逐时期差异图（2009–2008、2010–2009、2011–2010、2011–2008）揭示研究区内侵蚀/沉积呈交替格局，最强变化集中在中央活跃带和水工设施附近。体积统计显示2010–2009时期侵蚀体积最大（约70,000 m3），2009–2008时期沉积最大（约50,000 m3）；受影响面积方面，2010–2009侵蚀区约220,000 m2，2011–2010沉积区约210,000 m2。
- **4.2.2 水位参考岸线状态**：以2011年两个记录水位（43.5 m和45.8 m，PL-EVRF2007-NH）为例，2.3 m水位差导致暴露和淹没范围显著变化，表明低坡度近岸区域岸线状态对水位敏感。
- **4.2.3 形态变化横断面表达**：沿两条代表性断面（近坝段Profile 1、中央段Profile 2）提取多时期床面高程曲线和趋势预测2011剖面。Profile 1显示上游段垂直变化达2–3 m，下游段稳定；Profile 2呈现更复杂起伏。趋势预测剖面整体拟合良好，但在床面复杂区域存在局部偏差。
**4.3 自动化、可重复性与预测扩展**：实现了可重复的GIS工作流，将数据获取、预处理、插值、水位赋值、验证和输出生成整合，包括近实时水文更新模块（从IMGW API摄入当前水位）。基于2008–2010序列的逐像元趋势预测2011年床面，再现了主要空间模式，但中央活跃带和设施附近偏差较大。
**4.4 预测验证与残差分析**：残差图显示大部分区域差异小，最大残差集中在中央活跃带和水工设施附近。定量验证基于81,477个有效栅格单元：R2 = 0.967，RMSE = 0.44 m，MAE = 0.28 m，偏差 = ?0.06 m，表明预测表面与实测高度一致，但高精度部分反映床面空间连续性而非完全预测能力。
### 讨论与结论
讨论指出，床面调整空间持续但时间不均，与W?oc?awek河段先前研究一致。水位作为主要分析变量而非附属元数据是核心发现，使DT对水文条件变化响应。多维栅格结构体现分析核心，但预测应视作短时探索性延伸而非确定性模型。不确定性需显式处理（如通过暴露/淹没/过渡区分类）。框架在描述、预测和决策支持功能上具有条件可转移性（需多次测深、统一参考、水位记录等）。局限性包括测深时期有限（仅四期）、未纳入分布式水力学变量，未来应扩展自动化水文摄入、延长时序并集成水力-沉积过程。
**研究结论翻译**：本研究开发并验证了一个基于GIS的时空框架，用于W?oc?awek大坝下游受调控河流环境的动态岸线状态重建。该框架作为数字孪生导向的地球空间系统实现，将多时期测深调查、水位记录、ALS派生地形背景和三维水工基础设施整合于统一地理空间环境。所提出的方法超越了传统战役式评估，转向在受调控水文条件下河床变化和岸线状态变化的显式时间模型。结果表明，水位应作为主要分析变量而非附属描述信息。将测深表面与同期水位关联改善了岸线解释，特别是在低坡度近岸区域，而多维时序栅格结构支持可重复变化检测、横断面分析和趋势评估。为将框架延伸至纯粹回顾性分析之外，通过将IMGW API的操作观测数据摄入地理数据库实现了近实时水文更新组件，为岸线状态重建中所用水文参考状态提供了实用的物理-虚拟同步路径，但测深组件仍为基于战役式。所采用工作流适用于环境尺度形态动力学分析。预测扩展再现了2011年实测床面的主要空间格局，但较大残差仍集中在复杂形态区域和水工设施附近重复调整区域。因此，该预测步骤应视为短时探索性概念验证而非完整形态动力学预报模型。总体而言，所提出框架不仅作为可视化环境，而且作为结构化的、部分可更新的分析模型，将形态学、水文学、地形和基础设施在显式时间GIS框架中连接。其可转移性依赖于协调的测深、水文和地形数据的可用性。在这一更广泛视角下，该框架将描述、预测和决策支持能力整合于单一时空环境，用于受调控河流岸线监测。随着持续发展，特别是在自动化水文更新、更长时间序列和显式定量不确定性建模方面，可能演变为更接近运行的数字孪生导向系统，用于监测受调控河流岸线。
