采用透水网格法的研究伶仃洋鸡抱沙抛石促淤作用.pdfVIP

采用透水网格法研究 伶仃洋鸡抱沙抛石促淤作用木 包芸姚冬静 (中山大学近岸海洋中心，广州510275) 摘要：本文采用自主开发的二维河口海洋水沙特性数值计算模式ccost-2d，建立了河口抛石线模 型。通过透水网格法模拟抛石促淤和围垦，以鸡抱沙为例对伶仃洋围垦前后水动力变化和沉积速 率进行分析，并对自然演变和人类活动影响强度进行量化。结果表明抛石后洪季的淤积厚度是自 然演变的1～2倍。透水网格法可以有效的模拟人工抛石促淤演变过程。 关键词：透水网格法；抛石促淤；伶仃洋；C．c．OSt模型 伶仃洋河口湾是珠江最大的河口湾，是珠江主要的人海口，汇集了珠江八大口门中的东 四El门的来水来沙，是陆路、水路的结点。特别是伶仃洋地处广州、香港之间，是广州和黄 埔港及内河航道等的主要出海水道。其滩槽演变趋势与珠江三角洲防洪排涝，黄浦港和伶仃 洋出海航道的维护、整治以及海涂资源的开发利用等关系较为密切，直接影响这些地区未来 经济发展。近几十年来伶仃洋周围地区经济快速发展，人类活动频繁，各种近海工程对伶仃 洋的环境演变产生了巨大的影响。以往对近海工程多采用物模对比等方法进行研究。本文采 用自主开发的ccost-2d模型，用透水网格法模拟伶仃洋近几十年地形演变过程中最重要的人 类活动——抛石围垦，并以鸡抱沙为例分析抛石前后伶仃洋水沙特性变化。 1 ccost．2d模式及抛石线模型 1．1水动力模式 及荷兰的DELFT模式和丹麦模式等。这些模式在珠江河口的水动力计算中各有优缺点。ADI 方法和汉堡模式是在矩型网格上进行，不利于珠江河口复杂地形的模拟。POM模式为显式 格式，不方便后期的长周期模型的开发。ECOM模式是在POM模式的基础上发展起来的， 水位计算改用半隐式计算，使计算效率有所提高。但ECOM模式计算中科氏力项存在弱不 稳定性，在计算的编程上存在着所有网格点都进行计算求解，而后再把水点的物理量乘l， 陆点的乘。的缺陷，计算中有大量的浪费。还有许多环境演变的模式，程序规模都非常大， 只适合应用，不便进行新的研究的再开发和改造。 ·973课题资助项目(2002CB412407)，国家基金资助项目 109 在以上总结现有的多个海洋模式的基础上，我们开发了自己的二维水动力模式CCOSt．2d。 该模式的水动力模块有以下特点： (1)采用平面正交曲线网格， (2)采用平面有限体积法求解水动力方程组， (3)通过旋转矩阵克服科氏力产生的计算弱不稳定问题， (4)水位计算类似于ECOM采用隐式技术， (5)采用冻结网格法模拟滩地出没水面的动边界技术， (6)对水点记数，只计算求解水点的物理量。 Ccost-2d模型中也包含了泥沙输运和河床变形的模拟。 1．2透水网格法模拟抛石促淤 抛石的目的是消弱水势、束夹水流以稳定河槽，拦沙淤涨，加速滩地的淤积。抛石堤的 存在，并非完全隔断水沙的运动，为此我们借鉴空气动力学中飞机透气壁的模拟思路，采用 透水网格法来模拟现实中抛石堤的作用。 抛石堤的宽度与其长度相比可忽略不计，模拟过程中，将原来仅起流场划分作用的计算 网格线，在抛石线处加上透水率的限制，使流过该网格边界的水、沙被阻而又不完全被阻断， 减少通过流量，使该透水网格线起到模拟抛石的作用，如图1。 透水网格线 计算网格线 图1透水网格线与计算网格线示意图 在透水网格线上设定透水率盯(0盯I)，透水率是水位高低的函数。当盯=0时， 该网格线不透水；当or=l，该网格线成为原来的计算网格线，不阻碍水沙的流动；当0仃1 时，流过该网格线的水、沙被阻，通过流量减少，从而起到模拟抛石的作用。透水网格法关 键问题是透水率盯的确定，透水率仃的大小由多次数值模拟对比确定。 2模型计算结果分析 模型计算了上世纪80年代伶仃洋鸡抱沙围垦过程中抛石堤在促淤工程中的作用。模型 计算区域如图2所示。为了能更好的对比人类抛石下的冲淤与自然冲淤的区别，本文在同样