由于河口盐水入侵灾害的日益突出,该问题的研究受到了人们的关注。通过建立水流、盐度数学模型,进行河口盐度输移的数值模拟,对河口盐度的分布特征以及其受河口各种水动力因素的影响进行探讨分析,将有利于河口盐水入侵的防治。本文跟踪了国内外河口盐水入侵问题研究的最新动向,采用数值模拟方法对河口盐水入侵展开了系列研究,内容涉及一维、一、二维联解及三维的水流、盐度数值模拟、模型建立和运行中一些关键问题的处理、径流量、潮差、波浪和风对河口盐水入侵的影响分析、以及珠江河口伶仃洋的盐水入侵模拟。论文采用多套实测资料对建立的模型进行了验证,合理地再现与分析了河口区的盐水入侵现象。论文的主要内容可以归纳为以下几个方面:1.系统归纳与评述了河口盐水入侵的危害、复杂性及各种影响因素,总结了盐水入侵的主要研究方法和国内外盐水入侵问题的研究现状。2.针对网河区水系复杂交错、河口区较宽广的不同特点以及实际工程研究的不同需要,分别建立了一维和一、二维联解水流、盐度模型。结合具体的工程实例,进行了模型的验证及工程方案对盐水入侵影响的评价。此外,还对盐水入侵强度与径流量、潮差的相关关系进行了探讨。3.推导了水平正交曲线变换和垂向σ变换后的三维非线性斜压流场、盐度场数学模式,给出了离散求解过程,对一些模型参数和边界条件作了处理,并首次将波浪和风因素纳入河口盐度数值模拟系统。并且,对常用的干湿判别动边界处理法进行了改进,通过计算水域边界邻近网格边长的系数修正,更准确地定义了计算区域,反映出了水域边界移动的连续性。利用弯曲分汊河段—洲头嘴的实测资料,进行了三维水流计算和动边界处理的验证。4.建立了理想河口径流、潮流、温度、波浪和风综合影响下的三维水流、盐度数学模型,通过对相同边界、初始条件下,有浪和无浪、有风和无风情况下的河口盐水入侵进行模拟、对比分析,探讨了波浪和风对理想河口盐水入侵的影响。5.建立伶仃洋多因素综合作用下三维水流、盐度数学模型,利用大量实测数据对模型进行了较全面的验证,并对伶仃洋的三维盐度分布特征以及波浪和风对其盐度输移的影响进行了研究。此外,探讨了盐度初始场与循环计算的问题。论文的创新性体现在以上2~5所述的内容当中,通过本文的研究,可以得到以下有价值的结论:1.从本文中各个实例的计算结果可以看出,文中所建立的模型能正确反映出工程方案对盐水入侵的影响及河口盐度三维空间分布特征,采用一维、一、二维联解及三维水流、盐度数值模拟方法研究河口上游三角洲及河口区的盐水入侵问题是确实可行的。此外,动边界处理方法的改进、以及盐度初始场与循环计算问题的讨论,对于提高数值模拟的准确度和效率也是非常有利的。2.径流和潮流是盐水入侵的两个主要影响因素。盐水入侵距离和各个横断面含盐度均随着上游径流量的增加而减小,近似于线性变化关系。盐水入侵距离和各横断面含盐度均随着下游潮差的增加(涨潮量增大)而增加,其中盐水入侵距离与潮差呈线性变化关系,而横断面含盐度与潮差之间的变化近似于单调递增的曲线关系。3.理想河口和珠江河口伶仃洋的计算充分表明了波浪对河口盐度输移的影响。在本文所加波浪场作用下,河口盐水入侵前沿锋向河口上移,盐度等值线增加,总体而言,波浪越大时其对盐度输移的影响也越大,并且考虑波浪后表层盐度的变化大于底层盐度的变化。此外,也说明了波浪对盐度输移有两方面的影响:一是改变了水流流动,间接影响盐水入侵;二是直接地影响盐度的紊动扩散过程。4.风对河口区盐水入侵也存在一定的影响。在与涨潮流同向的风场作用下,受表面应力作用,河口区整体(表层至底层)盐度均有所增加,但由于实际河口情况复杂、风的存在同时会影响到深槽与浅滩的分流,因而,局部可能出现不一致的情况。较宽阔的外海区,其表层盐度增加,底层由于垂向风生环流的作用,也会出现盐度减小的情况。总之,从网河区一维、网河区和河口区一、二维联解到河口三维水流、盐度数学模型,从建模前期资料的整理到模型的求解、参数的确定及具体实例的计算、模型的验证,从径流、潮差到波浪和风各因素对河口盐水入侵影响的分析,本文的工作涉及到河口盐水入侵数值模拟研究的诸多方面。本文对盐度模拟研究中的一些重要问题进行了分析。本文给出的模型与经验为河口盐水入侵的进一步研究提供了有益的参考。