[发明专利]一种陡坎河床洪水冲刷演变模拟方法及系统

说明书

技术领域

本发明涉及模拟仿真领域，尤其涉及一种陡坎河床洪水冲刷演变模拟方法及系统。

背景技术

强烈的地壳运动或河道采砂等人为因素可导致河床局部塌陷或抬升，形成水流落差较大的陡坎。陡坎导致水流流态的强烈变化，反过来，水流强烈搬运床面泥沙导致陡坎河段河床的快速演变。洪水期，这样的演变更加剧烈，直接或间接导致桥梁等跨河工程的毁坏，给生产和生活带来安全隐患。2008年我国汶川地震形成的堰塞坝、引起2002年某铁路大桥坍塌的河床大规模采砂坑都是具有灾害隐患的陡坎河床。洪水作用下，陡坎河床泥沙冲刷演变，致使大桥基础变浅，承载力不足而坍塌。研究陡坎河段河床演变的数值模拟方法，可用于河床演变实时预报，对涉河工程防灾预报非常有意义。本发明特别针对因大规模采沙、地震诱发的堰塞坝形成的陡坎河床在洪水过程的冲刷与淤积问题，采用一维非恒定水动力方程、推移质泥沙输运方程及河床变形方程，考虑了间断点流体数值计算格式，并将有限体积法与有限差分法混合使用，提高了计算速度和计算精度，最终还通过计算机实现了数值模拟。本发明可用于特殊地形的河床演变预测和跨河交通设施防洪的安全管理中，实现随洪水过程实时掌握河床冲淤演变状况，避免因洪水冲刷而导致的工程设施毁坏或灾害。

现有技术中陡坎河床演变的研究现状主要有以下问题：1)在目前的跨河工程的安全条例中，仅根据工程规模的级别来规定工程与河床陡坎的安全距离，未全面考虑河道中的水流和地形条件引起的河床演变所带来的隐患；2)陡坎河床的冲淤演变定性分析较多，而少见用水、沙运动及河床变形理论进行定量和高精度数值模拟、预测的方法；3)陡坎附近的水流是大梯度跨临界流动，存在间断，传统的水流计算数值方法——差分法、有限体积法在计算陡坎间断点处容易出现数值振荡，这是由于水流的控制方程要求计算变量在空间点处处连续，而陡坎附近流态不满足该要求；4)任何一种数值格式难以同时满足数值计算的如下要求：处处连续、时空高阶精度和计算的高效率。

发明内容

本发明的实施例提供了一种陡坎河床洪水冲刷演变模拟方法，针对上述问题1)和2)，提供一种陡坎河床洪水冲刷演变的数值模拟方法，能快速、实时地定量预测洪水冲刷，给工程安全防护提供重要信息，避免灾害的发生或减轻灾害损失；针对上述问题3)，提供处理陡坎河床间断流的水流方程混合数值离散方法，提供处理间断面流动的守恒格式，和流体物质识别方法；针对上述问题4)，提供了混合离散格式数值解的时、空数值重构法，使方程组的解同时满足时间和空间二阶精度，且较仅用有限体积法计算效率高；为了本发明的实施，编制了C++语言计算程序。

本发明提供了如下方案，包括：

采用空间离散分项混合法对水流向量方程：进行处理；

其中，变量通量源项其中，为水面梯度，y为水位，y＝h+z_b；h为水深；q＝uh，即单宽流量；u为水深平均流速；S₀和S_f分别为底部坡度和摩阻坡度；g为重力加速度；

采用Godunov-Riemann格式对进行空间离散分项混合法处理的水流向量方程求解，获取关于变量U的解；

采用差分法离散求解源项S中的水面梯度并针对水面梯度获取离散后的水面梯度项，对离散后的水面梯度项进行间断点的空间数值重构，进一步求解右界面水位临界值和左界面水位临界值

根据推移质泥沙输运方程求解单宽输沙率q_s，

对河床变形方程：采用后差分格式的有限差分法离散；得到方程：再根据所述单宽输沙率q_s，和修正后的变量U,求解河床面高程z_b；其中，Δt表示时间步长，Δx表示空间步长，λ为河床孔隙率。

根据本发明的上述方法，所述采用空间离散分项混合法对水流向量方程：进行处理，包括：对惯性与对流项采用有限体积法离散，源项S采用中心差分格式的有限差分法离散。

根据本发明的上述方法，包括：基于所述对惯性与对流项采用有限体积法离散，源项S采用中心差分格式的有限差分法离散获得的方程，在间断点采用Godunov同时间层分段常数函数法，在间断界面的通量F用Riemann算子HLL求解。

根据本发明的上述方法，所述获取关于变量U的解，包括：对采用Godunov-Riemann格式所求的关于变量U的解在间断点上采用时空数值重构，修正所述关于变量U的解的精度。

根据本发明的上述方法，所述修正所述关于变量U的解的精度，包括：