论文部分内容阅读
海底地下水排泄和海水入侵是两个相互作用、相互影响的处于动态平衡的过程,为了进一步了解稳定状态下均质各向异性潜水含水层中海水与地下水之间的相互作用,本文利用数值模拟法重点研究了渗透系数、各向异性比、纵横向弥散度之比和海底坡度对含水层中盐度分布和海底地下水排泄量的影响。其中,渗透系数的变化范围为10-9m/s~1m/s(弱透层到含水层、粘土层到砾石层),各向异性比的取值范围为1~100,纵横向弥散度之比的取值范围为1~500,海底坡度的变化范围为0.01~?2(弧度)。模型区域内陆边界为给定水头边界、盐度为零,区域底部为隔水边界,上边界(地表)为不考虑降雨和入渗的隔水边界,海向边界的压力水头为作用于其上的海水水柱的压强,当海水流入含水层时为给定浓度35g/L,当水流流出含水层时为零弥散边界条件。采用二维数值代码MARUN来求数值解,空间离散采用三角形有限单元法,依据海底坡度的不同,所考察区域的海向延伸长度不等,结点数在16605到41123间变化,三角形单元数在32320到80160之间变化,网格的Peclet数均小于2,Courant数小于0.95,计算过程采用Picard迭代法,压力水头迭代收敛标准为10-5m。研究结果表明:当zx L?KK?10-6m2/s时(其中zK是垂直渗透系数、L?是纵向弥散度),盐度的空间分布几乎独立于水平渗透系数xK,且此时地下淡水排泄量fQ和海水再循环量sQ均与水平渗透系数xK呈线性关系,同时用一个半解析方法验证了这一系统性结论;各向异性比增大时,咸淡水界面向海向移动、并使海水再循环量sQ减小,消除了相关文献对此的争议;给定纵向弥散度后、增大纵横向弥散度之比会使混合带变倾斜并降低海水再循环量sQ;当海底坡度?从0.01增大到?2时,咸淡水混合带会向内陆方向移动,当?<?4时,海底坡度?的增大会使海水再循环量sQ有轻微的增加,当?>?4时,海底坡度?的增大会降低海水再循环量sQ。这些研究结果填补了前人在研究稳定态海水-地下水相互作用时时尚存的空白,进一步对影响SGD的因素进行了挖掘与分析,具有重要的意义。