码头工程建设对长江中游典型河段水流影响的数值模拟研究

　　摘要：为了研究码头工程建设对荆江公安河段河道水流的影响，通过建立平面二维水流数学模型，计算分析了工程建设前后附近河段的水位、流速变化情况。结果表明：工程建设对其所在河段的水位和流速影响较小，不会对工程河段河势、行洪带来较大影响。建立的二维水流数学模型可准确模拟工程河段水流运动特性，可为研究工程前后河道水流、水位变化提供较为可靠的依据。

　　关键词：水位;流速;数学模型;码头工程;公安河段;长江中游

　　1概况

　　1.1河段概况

　　长江是我国第一大河流，其干支流通航里程已超过6.5万km，占全国内河的52%，水运量占80%，是沟通我国东、中、西部地区的水路运输大动脉，在流域经济社会发展中具有极其重要的地位，素有“黄金水道”之称。目前，长江干线已成为世界上运输最繁忙、运量最大的通航河流。与其他交通运输方式相比，长江水运具有占地少、成本低、能耗小、污染轻、运能大、效益高等优势。随着水运经济的发展，长江沿岸码头工程的数量越来越多，码头对附近水域产生的阻水作用很有可能对该河段的河势及行洪产生影响。因此，研究码头建设对其所在河段水流的影响极为必要。本文以荆江公安河段左岸某码头工程为例，采用平面二维水流数学模型，对工程建设前后的河道水流流速、水位的变化情况进行了模拟计算，为该河段工程建设提供技术支撑。

　　长江公安河段上起观音寺，下止新厂，全长约54km，由陡湖堤河弯、郝穴河弯及过渡段组成，河型为微弯分汊。据历史资料记载，白元朝至1756年，陡湖堤至郝穴河段为单一的微弯河段，从观音寺到马家寨，河道紧靠现在的荆江大堤。1830年后，随着陡湖堤河弯南岸不断崩塌，河面不断展宽，逐渐形成突起洲。郝穴河弯中的采石、白沙、新淤、新泥和白脚等5个江心洲将河道分为左、右两汊，左汊为主河槽，右汊即黄水套。1852年起，郝穴一带开始實施护岸，由于历年泥沙的不断淤长，采石洲等5个江心洲逐渐联成一个整体，成为现在的南五洲。近30a来，公安河段河势演变基本平稳，岸线、深泓及滩槽位置相对稳定。随着三峡水库蓄水运用，大量泥沙被拦截在库内，下泄泥沙量大幅度减少，公安河段冲刷强度普遍增[1-2]，公安河段河势见图1。

　　1.2工程概况

　　码头工程位于荆州市江陵县郝穴镇城区上游，长江左岸的郝穴作业区，对应荆江大堤桩号为鄂江左712+492-712+960。上距在建的荆州长江公铁大桥约4.5km、距沙市水文站约45km，下距江陵县城约3km。

　　码头工程水工建筑物主要为4个3000吨级货运泊位，包括高桩直立式码头1座，浮码头2座(见图2)。为了减少码头对荆江大堤的不利影响，码头采用顺岸布置，码头前沿线与上游中航粮油综合码头前沿线呈折线，夹角为174°，布置丁21-22m等高线附近，与水流方向基本一致。泊位总长468m。

　　高桩码头平台长220m，宽25m，顶标高40.5m。码头平台通过2座引桥和1条皮带机廊道分别与大堤和后方陆域廊道相接，上游1号固定引桥长58.4m，宽12m，下游2号固定引桥长75m，宽12m，2座引桥与大堤平交后，通过下堤道路与规划道路连接。皮带机廊道布置在1号固定引桥的上游侧，采用多跨钢引桥结构跨堤，并与后方陆域廊道对接。

　　浮码头每个泊位从江侧往大堤方向均设置1艘75.0mx16.0m(长×宽)钢质趸船、2座钢撑杆系统、1座60mx4.5m(长×宽)活动钢引桥和1座转运楼墩台，墩台顶高40.77m。两个转运楼墩台之间各通过1座钢引桥连接到2号转运楼墩台，墩台顶高程47.70m，2号转运楼墩台通过固定钢引桥与后方陆域连接，引桥全长360.74m，宽8m。

　　2模型原理

　　本次数学模型基于正交曲线坐标系下沿水深平均的平面二维水流运动方程组，微分方程的数值离散采用有限体积法[3-8]。定解条件包括初始条件与边界条件，边界条件为上游给定垂线平均流速沿河宽的分布，下游出水位沿河宽的分布。对于岸边界，则采用水流无滑移条件，即取岸边水流流速为0。在计算时，由计算开始时刻上、下边界的水位确定初始条件，初始流速取0，随着计算的进行，初始条件的偏差将逐渐得到修正，对最终计算结果的精度不会产生影响[9]。

　　3模型建立与验证

　　3.1模型建立

　　模型选取观音寺附近至茅林口全长约49.15km的河段作为二维数学模型的计算区域(见图1)。计算河段的地形资料采用2016年10月实测地形资料。

　　计算网格采用正交曲线网格，在工程河段布置600x120个网格，其中ξ方向网格数600个，平均断面间距8lm;η方向网格数120个，平均网格宽度约33m。同时，为使网格大小与工程实际尺寸一致，对工程附近的局部河段进行了网格加密处理。

　　3.2模型验证

　　采用2014年12月该河段实测的水面线与垂线平均流速资料进行水位与流速验证。水位验证成果见表1，流速验证成果见图3。水位、流速计算值与实测值基本一致，水位误差大部分在5cm以内，流速误差在0.1m/s之内，符合JST/T231-4-2018《内河航道与港口水流泥沙模拟技术规程》要求。因此，模型中相关参数的取值是合理的，可以模拟码头修建对河道水位与垂线平均流速的影响。

　　4工程影响计算分析

　　4.1工程概化

　　码头工程主要阻水建筑物为桩柱、墩台，在防洪设计洪水时还有现浇承台及部分引桥桥面。为在模型中反映码头工程对河道水流的影响，在网格剖分时，尽可能在码头附近对网格进行局部加密，同时采用局部地形修正与局部糙率调整来进行概化处理。

　　推荐阅读：《水资源保护》杂志是河海大学和环境水利研究会合办的科学技术期刊。该刊以宣传贯彻国家水资源保护方针政策，反映国内外水资源保护先进科学技术，促进我国水资源保护和管理以及水环境保护工作为宗旨。