[发明专利]一种考虑汊点面积的通用河网水流数值模拟方法

摘要

本发明公开了一种考虑汊点面积的通用河网水流数值模拟方法。获取河网平面几何形状数据和各河段横断面数据，建立各河段和汊点的拓扑关系以及各汊点处的水位面积关系。河网采用有限体积单元离散，水力要素值存储在单元中心。通过建立河网汊点处水量平衡方程，求解汊点处的水位值，该水位值将为与该汊点相连接的河段提供边界条件，汊点提供的水位边界条件与各河段的外边界条件配合后，每条河段的首末端均具备了边界定解条件，再通过Godunov显格式有限体积法来分别模拟各河段水流运动过程。该方法既能处理平原区的缓流流态，也能很好处理山丘区河网的复杂流态和大梯度水面的模拟问题，具有很好的通用性，能够有效弥补现有平原河网水流模拟方法的不足。

主权项

1.一种考虑汊点面积的通用河网水流数值模拟方法，其特征在于：通过建立河网汊点处水量平衡方程，求解汊点处的水位值，为与该汊点相连接的各条河段提供边界条件，将该边界条件与各河段外边界条件配合，通过Godunov显格式有限体积法来模拟各河段水流运动过程，包括以下具体步骤：1)获取河网平面几何形状和各河段横断面数据，各横断面间的间距不大于1km；2)对各河段和汊点进行编号，建立汊点和各河段的对应关系；采用有限体积单元离散各河段，水力变量存储在各单元中心点；根据河网平面几何形状和首末端单元的断面形状建立各汊点处的水位面积关系Z‑A；3)模型初始化：给各河段一维河道单元设置糙率和水力条件初始值；根据各河段首末单元的初始水位获取汊点处初始水位；4)确定计算时间步长dt：一维河道模型采用显格式的有限体积法求解，该时间步长获取受到CFL条件限制；获取当前时刻各河段外边界条件，当前时刻以t时刻表示；5)采用完整的一维Saint‑Venant方程组描述河道水流运动，采用HLL显格式计算通过单元界面处的数值通量，通过t时刻各单元的水力初始值、外边界条件以及汊点初始水位计算获得各河段t+dt时刻各单元的水力要素值；6)建立汊点处水量平衡方程，并根据各汊点水量平衡方程获得t+dt时刻各汊点的水位值；采用如公式(3)的各汊点水量平衡方程计算t+dt时刻各汊点的水位值：式中：为第k个汊点在t时刻的水面面积，为第k个汊点在t时刻的水位，为需要求解的第k个汊点在t+dt时刻的水位，n为与第k个汊点相连接的河段个数，为与汊点k相连接的第i个河段在t时刻流入/流出汊点的流量，流入为正，流出为负；7)令t＝t+dt，重复步骤4)～6)，直到计算结束。