文章提出了一种耦合水动力与形态动力学的泥沙输移模型（HSWEMED），解决了传统浅水模型在非平衡条件下预测精度不足的问题。该模型通过Legendre多项式展开捕捉垂向流速结构，直接计算近床流速和剪切应力，并耦合床载与悬移质输移方程。创新性地引入双曲正则化技术保证模型稳定性，推导了包含侵蚀-沉积过程的Exner方程闭合系统。数值实验表明，模型在溃坝问题中较经典模型更精确地复现了自由水面演变、河床变形和悬沙浓度分布，尤其在PVC颗粒等轻质泥沙案例中显著提升了预测准确性。
传统浅水泥沙输移模型依赖深度平均流速和经验对数剖面估算近床剪切应力，难以准确描述非稳态流中的垂向结构。文章针对此局限，提出融合矩方法的多过程耦合框架，核心创新包括：
基于非均质Navier-Stokes方程，通过垂向标度变换和Galerkin投影导出矩方程。质量守恒方程表示为：
\[ \partial_t h + \partial_x(h u_m) = F_b \]
其中床面交换率 $ F_b $ 关联侵蚀-沉积过程。动量方程新增三项关键耦合项：
悬移质浓度方程与Exner方程构成闭合系统：
\[ \partial_t(hc) + \partial_x(h c u_m) = (1-\psi) F_b \]
\[ \partial_t h_b + \frac{1}{1-\psi} \partial_x Q_b = -F_b \]
侵蚀-沉积模型采用Garcia-Parker公式，沉降速度由颗粒雷诺数动态计算：
\[ \omega_o = \sqrt{ \left( \frac{13.95\nu_w}{d_s} \right)^2 + 1.09 \rho_w \left( \frac{\rho_s}{\rho_w}-1 \right) g d_s } - \frac{13.95\nu_w}{d_s} \]
通过正则化技术解决高维矩模型特征值虚根问题。定义双曲矩阵 $ \tilde{A}_{ED} $ 满足：
\[ \partial_t W + \tilde{A}_{ED} \partial_x W = S(W) \]
其特征多项式分解为：
\[ \chi_{\tilde{A}_{ED}}(\lambda) = (u_m - \lambda) \left[ -\lambda ((\lambda - u_m)^2 - gh - \alpha_1^2) + gh(\delta_h + \lambda \delta_q + c\delta_c + 2\alpha_1 \delta_q) \right] \cdot \chi_{A_2}(\lambda - u_m) \]
该形式确保所有特征值为实数，证明模型具备物理相容的有限传播速度。
对比HSWEMED、经典SWEED和HSWEM模型：
针对三种初始配置（平坦/非连续河床、干湿界面）和两种床质（PVC/砂）：
HSWEMED模型通过矩方法有效统一了床载与悬移质输移机理，其优势体现在：
未来研究将扩展非恒定浓度剖面建模，并开发适应干湿前沿的高阶数值格式。