二维潮流泥沙数值模拟在福姜沙河段深水航道治理工程中的应用.pdfVIP

第十三届中国海洋(岸)工程学术讨论会论文集 483 二维潮流泥沙数值模拟在福姜沙河段 深水航道治理工程中的应用 闻云呈，夏云峰 (南京水利科学研究院，江苏南京210029) 摘要：通过对长江下游河口段福姜沙分汊河段水沙特性分析，在把握本河段河床演变规律的基础上，建立二维非恒定潮流 非平衡非均匀输沙数学模型，经水动力和河床冲淤验证，模拟了天然汉道的潮流传播和潮流运动以及泥沙输移引起的河床变 化。利用本模型研究了潮汐分汉河道航道整治工程对各汉道水动力条件的影响及河床冲淤变化，分析工程所处的潮流界变动 区，不同径流潮汐共同作用下各汉道的输水输沙特性，以及各汊道航道的变化规律。 关键词：长江下游航道整治工程；冲刷；淤积；潮流泥沙数学模型 1概述 随着长江三角洲经济迅猛发展以及长江口深水航道整治工程的实施，5万吨级海轮将由长江口直达南 京。根据交通部交规划发文“关于长江口深水航道向上延伸有关干线航道发展规划的批复”文件 精神，长江口深水航道要上延至南京，到2020年全河段实现12．5m的深水航道维护尺度。长江南京～浏 河口航道目前有多个碍航浅滩，航道尺度标准偏低，其中福姜沙已经成为长江下游重点碍航航道，严重制 约了长江下游黄金水道所发挥的效益。因此对福姜沙水道海轮深水航道碍航浅滩段的整治迫在眉睫。由于 本河段受上游径流及下游潮汐影响，水流条件复杂，滩槽冲淤变化、汊道兴衰将影响航道的维护和发展。 拟通过二维潮流泥沙数学模型计算，分析各种计算水文泥沙条件下比选整治方案实施后航槽的冲淤变化。 2二维潮流数学模型建立 对于水平尺度远大于垂直尺度的河流，湖泊和海洋等水域，其物理量沿水深方向的变化相对于水平方 向变化小的多时，可将各物理量沿水深积分，得水深积分二维水流泥沙运动方程。且假定：①水流为均质 不可压流体；②水流为常密度并服从静水压力分布：P=P。+店(f—Z)；③刚盖假定，忽略垂向加速度： ④床面切应力采用二次形式，紊动切应力采用Boussinesq假定；⑤由流速沿垂线分布不均匀在积分时产生 的修正系数为1．0。 2．1笛卡儿坐标系二维水深积分水流运动基本方程 2．1．1水流运动基本方程 连续方程： (1) af ———兰—一J_ !也±兰幽+!险±￡幽：o 动量方程： at dx 0Y (2) 瓦栅面w 扣妻+Vauay=未(屹妻)+孙势g豢+刍一刍+声2瓦卜忑／J+面l屹面j—g素+茄一声+∥ (3) 一枷一¨一2万lVe万J+瓦I叱瓦j_g贵+。声一茄√“ 式中：M，v为水深平均流速在x,y方向分量，H为水深，日=忽+6 转角速度， 缈为纬度； y。为有效黏性系数： y。=l，，+V，y，为紊动黏性系数； ％=一“√：甭，％=pC,fV√丽分别为底部切应力在x,y方向分量；c，为底部摩擦系数： n为河底糙率系数。z。，彳掣分别为表面风应力在x,y方向分量：彳。=,ok。W，WI， c，=n29／H％ 丁掣=肚。w，lWI，W=、厂孺其中五=。为系数，本文计算中暂不考虑风应力的影响，令彳。，z掣为零。 第十三届中国海洋(岸)工程学术讨论会论文集 2．2笛卡儿坐标系二维泥沙输运方程 2．2．1悬沙不平衡输运方程 匀在积分时产生的修正系数： z 1 一 二IMlsdz