广州小虎二桥行洪影响数学模型计算.pdfVIP

广州小虎二桥行洪影响数学模型计算 郑国栋 黄东 张政 黄本胜 （广东省水利水电科学研究院，广州，510630） 摘 要：建立较大范围的一维河网潮流数学模型嵌套局部平面二维数学模型，运用有限差分、有限元法、Pressimann 隐格式、三级联解、Galerkin 法及Delaunay 三角元剖分等方法和技术，对不概化的桥墩直接模拟，并探讨性给 出局部损失系数与糙率的经验转换关系式，模型亦对此式进行了验证。 关键词：河网 三级联解 Galerkin 法 Delaunay 剖分 局部损失系数 糙率 1 引言 根据广州市城市总体规划的广州黄沙快速干线工程将通过小虎二桥横跨小虎沥水道中部，小 虎沥水道为广州片网河的重要行洪、纳潮通道之一，小虎二桥有13 排桥墩布置在小虎沥河道内， 桥墩及其基础将占用所在河段的行洪、纳潮过流面积，会对小虎沥水道的行洪、纳潮态势造成一 定程度的影响。珠江三角洲河网的特点是水道纵横交错、相互贯通，在上游径流及下游南海潮汐 动力的共同作用下，水流形态不仅十分复杂，而且互相关联。在局部河道上兴建工程不仅会改变 该河道的水流形态，还可能对附近其他河道造成相应影响。河网数学模型对于模拟河网区工程对 行洪纳潮的影响具有很强的优势，但这仅限于宏观一维计算情况，到二维河网计算则受计算机容 量和计算速度的限制，难于大范围应用，三维计算则更加困难了。所以目前较通用的做法是宏观 上采用一维河网数学模型计算，结合局部工程段应用二、三维数学模型（或物理模型试验）对水 流流态进行较精细的模拟。本文采用一维河网数学模型嵌套局部平面二维有限元数学模型模拟广 州小虎二桥的行洪影响。求解过程中用到了河网的三级联解、二维有限元数学模型的 Galerkin 法及 Delaunay 三角元剖分等技术。本文还对局部工程引起的局部水头损失与沿程损失的转换关 系，给出了局部损失系数与糙率的经验转换关系式，模型亦对此式进行了验证。 2 数学模型 2.1 一维河网数学模型 2.1.1 基本方程 数学模型的基本方程采用圣维南方程组： ∂Q ∂z ∂x + B T ∂t q l (2.1.1)  2 ∂Q ∂Q 2 ∂z 2 ∂A n u Q ∂t + 2 u ∂x + ( gA − Bu ) ∂x − u ∂x + g R 3 4 0 (2.1.2) 式中：Q 为流量；Z 为水位；R 为水力半径；u 为流速；ql 为旁侧入流；n 为糙率系数，可用谢才 公式计算；BT 为包括主河道泄流宽度和仅起调蓄作用的附加宽度；B 为过流河宽；A 为过水面积； g 为重力加速度；x、t 为空间和时间坐标。 补充河网节点关系式（2.2）（或（2.3））、（2.4）即构成一维网河数值计算的数学模型， 通过运用三级联解法离散求解（2.1）～（2.4），即可实现河网模拟。对于节点应满足下列流量 连接条件和动力连接条件： n dw ∑Qi (2.2) i 1 dt 192 式中，Q 为节点过流量；i 表示汇集于同一节点的各河道断面的编号；w 为节点蓄量。 若节点为无蓄量的几何点，则w=0。因此： n ∑Qi 0 (2.3) i 1 动力连接条件取为： Z =Z +1 k=1，2，3ΛK-1 (2.4)