一种基于水动力数值模拟的感潮河网潮位预报方法技术

本发明专利技术提供了一种感潮河网内潮位预报的方法，综合运用了计算流体力学(CFD)方法和水位预报方法，解决了感潮河网内水位受多汊道分流及外海潮汐影响而难以准确预报的问题，提高感潮河网潮位预报的适用范围和精度，属于水动力模拟和水文预报技术领域。主要包括以下步骤，1)建立感潮河网区及上游河道一维水动力数学模型；2)利用调和分析预测河口控制站逐时潮位过程并利用预报前一日实测数据进行预报值校正；3)然后采用加权移动平均方法推求上游控制站点流量过程；4)利用河网区及上游干流区一维水动力数学模型求解感潮河网区各断面潮位过程，实现河网区任意断面潮位预报。

【技术实现步骤摘要】

本专利技术涉及水动力模拟和水文预报
，特别涉及一种基于水动力数值模拟的感潮河网潮位预报方法。
技术介绍
随着河口及三角洲经济建设的不断发展，对水文预报工作提出了更高要求，如何快速准确地进行潮位预报，成为当前急需解决的问题。防汛决策、水利枢纽设计参考、闸泵群启闭方案优化等对感潮网河区的潮位预报需求日益增加。感潮河网的水位变化过程受上游径流、下游潮汐、河汊分流、区间入流及降水等多种因素共同影响，其水文特性既有内陆河流的特性，表现出洪、中、枯水期的季节变化，又有海洋的某些特征，表现出大、中、小潮的变化和潮汐引起水位涨落与流速流向的变化，此外，还受陆地气候的影响，上游径流不仅年内有明显的变化，而且在年际间变化也很不相同；受河道地形的影响也较大，情况十分复杂，感潮河网内潮位预报是长期以来河口及三角洲治理开发中的一个难点问题。传统的潮位预报多为高水位预报，一般采用相应水位法，即由某站水位、传输距离推求另一站水位，其预报有效期受传输时间左右，无法满足长期、中期预报需求；另有分离过程法将预报站水位分离成上、下游水位相关的两部分分别进行预报，其不足之处与相关法相同，主要是预报时间较为受限；神经网络法预报水位则主要存在可能由于局部极值难收敛、振幅不大时较精准、高水位预报误差大的问题。为满足三角洲对潮位预报需求，发展一种基于数值模拟的能实现感潮河网内任意水域潮位快速预报的方法是十分必要的。
技术实现思路
本专利技术的目的在于克服河网地区现有潮位预报方法误差较大、适用范围小的不足，提供一种基于水动力数学模型的感潮河网任意水域潮位快速预报及实时校正方法，提高感潮河网潮位预报的适用范围和精度。本专利技术的目的通过下述技术方案实现：一种基于水动力数值模拟的感潮河网潮位预报方法，包括以下步骤：步骤S1、根据预报位置选择合适的研究范围，确定河口下边界及上游控制站点与边界；步骤S2、对研究范围内河道进行断面剖分和地形插值；步骤S3、建立上游控制站至下游河口边界之间的网河区的一维非恒定流水动力数学模型；步骤S4、采用实测水文数据对一维非恒定流水动力数学模型进行率定和验证；步骤S5、根据河口控制站点调和参数进行潮位过程预报计算：步骤S6、根据实测数据对河口潮位边界进行校正；步骤S7、采用加权移动平均方法推求上游控制站点流量过程；步骤S8、根据实时实测数据对上游流量边界进行校正；步骤S9、输入预报上游边界流量和下游河口边界潮位过程，进行网河区的一维非恒定流水动力数学模型计算，求得网河区内各断面预潮位过程。进一步地，所述步骤S3、建立上游控制站至下游河口边界之间的网河区的一维非恒定流水动力数学模型具体过程如下：步骤S31、采用圣维南方程组作为一维潮流数学模型，方程如下：连续方程：动量方程：将相关支流汇入或流出点作为网河区汊点，汊口点水流按以下连接条件控制：流量连接条件：水位连接条件：Zi,j＝Zm,n＝……＝Zl,k；式中，Z为断面平均水位；Q，A，R分别为断面流量、过水面积、水力 半径；B为水面宽度；X,T为距离和时间；G为重力加速度；谢才系数C＝R1/6/n，n为糙率；步骤S32、采用四点加权Preissmann固定网格隐式差分格式将方程离散，采用一维河网三级联解算法求解方程，具体求解方法为：进行网格划分，以S代表流量Q和水位Z，则S在Δx河段、Δt时段内的加权平均量及相应偏导数可分别表示为：式中：θ为加权系数，θ一般取0.5～1.0。进一步地，所述步骤S5、根据河口控制站点调和参数进行潮位过程预报计算具体过程如下：步骤S51、有分潮参数的直接根据预报时间计算潮位过程；步骤S52、无分潮参数，但有长于一年的实测历史数据的，调和分析各分潮参数，然后重复步骤S51；步骤S53、分潮参数和历史数据均无则参考临近站点简化处理。进一步地，所述步骤S51、有分潮参数的直接根据预报时间计算潮位过程的具体计算公式如下：式中：为t时刻的潮位预报值，m；A0为从某基准面算起的平均海平面高度，m；fi为第i个分潮的交点因子；Hi为第i个分潮的振幅，m；σi为分潮角速度；(v0+u)i为分潮的天文初相角；gi为分潮迟角；i为分潮编号；N为分潮总数；Hi、gi亦称为调和常数。进一步地，所述步骤S6、根据实测数据对河口潮位边界进行修正具体过程如下：步骤S61、在[t0-26,t0]范围内搜索判别出高高潮位ZHPRE、低低潮位ZLPRE，计算潮幅ZAMPRE，ZAMPRE＝ZHPRE+ZLPRE；步骤S62、同样在[t0-26,t0]范围内搜索出实测潮位过程中的高高潮位ZHMESU、低低潮位ZLMESU，计算潮幅ZAMMESU，ZAMMESU＝ZHMESU+ZLMESU；步骤S63、根据步骤S61及步骤S62结果对预报潮位中轴线进行调整，使得ZAMPRE＝ZAMMESU；步骤S64、计算中轴线偏移量DIF_ZZERO，DIF_ZZERO＝(ZAMPRE-ZAMMESU)/2.0；步骤S65、对于预报潮位过程值ZPRE逐一进行中轴线偏移校正，ZPRE＝ZPRE-DIF_ZZERO；步骤S66、计算潮幅偏差DIF_ZAM，DIF_ZAM＝ZAMPRE/ZAMMESU；进一步地，所述步骤S7中加权移动平均方法的公式如下： Q p r e e t = Q A t * 2 + Q p r e e t - 1 ]]>式中：为t时刻预报流量，m3/s；为第t时刻的流量趋势率的修正量，m3/s；为第t-1时刻的预报流量，m3/s。进一步地，所述流量趋势率的修正量的求解步骤如下：步骤S71、对于每一段预报期(t0～tn)内的首、末趋势率：式中：分别为预报期内首时刻t0、末时刻tn对应流量趋势率；QT0-1和QT0-2分别为预报前1、前2小时的实测流量；步骤S72、按线性插值法计算每一段预报期内首末时刻以外其它各时刻趋势率，式中：为预报期内首末时刻以外其它各时刻趋势率；T为预报期总时长。本专利技术相对于现有技术具有如下的优点及效果：本专利技术公开了一种基于水动力数学模型的感潮河网任意水域潮位快速预报及实时校正方法，综合运用了计算流体力学(CFD)方法和水位预报方法，解决了感潮河网内水位受多汊道分流及外海潮汐影响而难以准确预报的问题，提高感潮河网潮位预报的适用范围和精度。附图说明图1是本专利技术公开的的一种基于水动力数值模拟的感潮河网潮位预报方法 的流程步骤图；图2是一维河网水动力模型断面剖分示意图；图3(a)是本专利技术在中顺大围调度系统中一维河网水动力模型研究范围图；图3(b)是本专利技术在中顺大围调度系统中中顺大围区域的研究范围图；图4(a)是98.6洪水期三沙口站水位过程验证成果图；图4(b)是98.6洪水期勒竹站水位过程验证成果图；图4(c)是05.1枯水期小榄站水位过程验证成果图；图4(d)是05.1枯水期三沙口站水位过程验证本文档来自技高网...

【技术保护点】
一种基于水动力数值模拟的感潮河网潮位预报方法，其特征在于，包括以下步骤：步骤S1、根据预报位置选择合适的研究范围，确定河口下边界及上游控制站点与边界；步骤S2、对研究范围内河道进行断面剖分和地形插值；步骤S3、建立上游控制站至下游河口边界之间的网河区的一维非恒定流水动力数学模型；步骤S4、采用实测水文数据对一维非恒定流水动力数学模型进行率定和验证；步骤S5、根据河口控制站点调和参数进行潮位过程预报计算：步骤S6、根据实测数据对河口潮位边界进行校正；步骤S7、采用加权移动平均方法推求上游控制站点流量过程；步骤S8、根据实时实测数据对上游流量边界进行校正；步骤S9、输入预报上游边界流量和下游河口边界潮位过程，进行网河区的一维非恒定流水动力数学模型计算，求得网河区内各断面预潮位过程。

【技术特征摘要】
1.一种基于水动力数值模拟的感潮河网潮位预报方法，其特征在于，包括以下步骤：步骤S1、根据预报位置选择合适的研究范围，确定河口下边界及上游控制站点与边界；步骤S2、对研究范围内河道进行断面剖分和地形插值；步骤S3、建立上游控制站至下游河口边界之间的网河区的一维非恒定流水动力数学模型；步骤S4、采用实测水文数据对一维非恒定流水动力数学模型进行率定和验证；步骤S5、根据河口控制站点调和参数进行潮位过程预报计算：步骤S6、根据实测数据对河口潮位边界进行校正；步骤S7、采用加权移动平均方法推求上游控制站点流量过程；步骤S8、根据实时实测数据对上游流量边界进行校正；步骤S9、输入预报上游边界流量和下游河口边界潮位过程，进行网河区的一维非恒定流水动力数学模型计算，求得网河区内各断面预潮位过程。2.根据权利要求1所述的一种基于水动力数值模拟的感潮河网潮位预报方法，其特征在于，所述步骤S3、建立上游控制站至下游河口边界之间的网河区的一维非恒定流水动力数学模型具体过程如下：步骤S31、采用圣维南方程组作为一维潮流数学模型，方程如下：连续方程：动量方程：将相关支流汇入或流出点作为网河区汊点，汊口点水流按以下连接条件控制：流量连接条件：水位连接条件：Zi,j＝Zm,n＝……＝Zl,k；式中，Z为断面平均水位；Q，A，R分别为断面流量、过水面积、水力半径；B为水面宽度；X,T为距离和时间；G为重力加速度；谢才系数C＝R1/6/n，n为糙率；步骤S32、采用四点加权Preissmann固定网格隐式差分格式将方程离散，采用一维河网三级联解算法求解方程，具体求解方法为：进行网格划分，以S代表流量Q和水位Z，则S在Δx河段、Δt时段内的加权平均量及相应偏导数可分别表示为： ∂ S ∂ t = S i + 1 n + 1 + S i n + 1 - S i + 1 n - S i n 2 Δ t ∂ S ∂ x = θ S i + 1 n + 1 + S i n + 1 Δx i + ( 1 - θ ) S i + 1 n ...

【专利技术属性】
技术研发人员：杨莉玲，胡晓张，杨芳，杨志伟，余顺超，刘壮添，
申请(专利权)人：珠江水利委员会珠江水利科学研究院，
类型：发明
国别省市：广东;44