第３６卷 第１１期海 洋 学 报
Ｖｏｌ．３６，Ｎｏ．１１２０１４年１１月
ＡＣＴＡＯＣＥＡＮＯＬＯＧＩＣＡＳＩＮＩＣＡ
Ｎｏｖｅｍｂｅｒ２０１４
陈泾，朱建荣．长江河口青草沙水库盐水入侵来源［Ｊ］．海洋学报，２０１４，３６（１１）：１３１－１４１，ｄｏｉ：１０．３９６９／ｊ．ｉｓｓｎ．０２５３-４１９３．２０１４．１１．０１５
ＣｈｅｎＪｉｎｇ，ＺｈｕＪｉａｎｒｏｎｇ．ＳｏｕｒｃｅｓｆｏｒｓａｌｔｗａｔｅｒｉｎｔｒｕｓｉｏｎａｔｔｈｅｗａｔｅｒｉｎｔａｋｅｏｆＱｉｎｇｃａｏｓｈａＲｅｓｅｒｖｏｉｒｉｎｔｈｅＣｈａｎｇｊｉａｎｇＥｓｔｕａｒｙ［Ｊ］．Ａｃ-ｔａＯｃｅａｎｏｌｏｇｉｃａＳｉｎｉｃａ（ｉｎＣｈｉｎｅｓｅ），２０１４，３６（１１）：１３１－１４１，ｄｏｉ：１０．３９６９／ｊ．
ｉｓｓｎ．０２５３-４１９３．２０１４．１１．０１５长江河口青草沙水库盐水入侵来源
陈泾１，朱建荣１*
（１．
华东师范大学河口海岸学国家重点实验室，上海２０００６２）收稿日期：２０１３-０８-３０；修订日期：２０１４-０６-０５。
基金项目：国家自然科学基金项目（４１１７６０７１）；上海市科学技术委员会重大项目（１２２３１２０３１０１）；水利部公益性项目（２０１２０１０６８）。
作者简介：陈泾（１９８８－），男，福建省龙海市人，从事河口海岸动力学研究。
Ｅ-ｍａｉｌ：ｃｈｅｎｊｉｎｇｆｊｘｍ＠ｓｉｎａ．ｃｏｍ*通信作者：朱建荣，男，教授，从事河口海洋动力学研究。
Ｅ-ｍａｉｌ：ｊｒｚｈｕ＠ｓｋｌｅｃ．
ｅｃｎｕ．ｅｄｕ．ｃｎ摘要：应用改进的三维数值模式ＥＣＯM -ｓｉ，从模式计算的盐度和流向的变化过程、涨憩和落憩时刻盐度等值线和淡水区域的变化，分析在一般动力条件下青草沙水库取水口盐水入侵来源。
计算结果表明，
小潮后中潮、大潮、大潮后中潮和小潮期间北支倒灌占青草沙水库取水口表层盐水入侵比例分别为６９．５％、８９．３％、９８．５％和９９．５％，占底层盐水入侵比例分别为３４．９％、８８．９％、９８．５％和９９．５％。
除了小潮后中潮期间底层盐水入侵来源主要来自下游外海（占６５．１％），青草沙水库取水口表层和底层盐水入侵来源主要来自北支盐水倒灌，尤其是大潮后中潮和小潮期间几乎全部来自北支盐水倒灌。
关键词：青草沙水库；盐水入侵；来源；数值计算
中图分类号：Ｐ７３１．１２
文献标志码：Ａ
文章编号：０２５３-４１９３（２０１４）１１-０１３１-１１
1 引言
长江河口地区，经济发达，人口密集，需要大量优质原水。
２０１０年前上海的用水主要取自黄浦江，水量不足，水质较差，为典型的水质性缺水城市。
长江河口水量充沛，水质优良，要从根本上解决上海的用水难题，必须从长江河口取水。
为了解决用水难题，上海市在长江口建成了大型河口江心水库———青草沙水库。
青草沙水库位于南北港分汊口附近、长兴岛西北侧（见图１），水域面积达到６６．１５ｋｍ２，于２０１０年建成并开始向上海供水，承担了上海市约５０％的原水供应，规划供水规模为７．１９×１０６ｍ３／ｄ，
受益人口超过１０００万人［１］。
青草沙水库的兴建，大大缓解了上海供水紧张的局面，为上海城市发展和工农业生
产提供有利的保障。
但从长江河口取水面临的主要问题是枯季盐水入侵，如果取水口盐度大于饮用水标
准０．４５（实用盐度单位，按国际惯例单位一般不标注）就不能取水。
长江河口枯季经常面临盐水入侵，最大的特点是南支除受外海盐水入侵外，还受上游北支盐水倒灌的
入侵［２］。
已有的大量观测和研究表明，潮汐和径流量是影响盐水入侵的主要原因，另外还受风应力［３—４］、口外陆架环流［５］和河势变化［２］等影响。
以往对长江河口盐水入侵来源的研究，沈焕庭等［２］指出潮汐河口
外海盐水入侵导致的盐度周日变化一般规律是盐度最高、最低值分别出现在涨憩、落憩附近，且周日变幅大。
长江口南支、南北港受北支盐水倒灌后，改变了盐度周日变化规律，
具体反映为盐度周日变幅小，周日的盐度峰值和谷值出现在落憩和涨憩附近。
茅志
昌等［６］通过分析现场观测资料，指出青草沙水源地的
盐水源自外海盐水入侵和北支盐水倒灌，
其中以受北支倒灌盐水团过境作用为主。
顾玉亮等［７］根据观测
资料分析，指出影响南支水源地的盐水来源有两个，
即北支盐水倒灌和南北港外海盐水入侵。
乐勤等［
１］根据多年监测数据指出，影响青草沙水库水域氯化物变化的咸潮入侵源有３个，即北支咸潮倒灌、南港咸
潮入侵和北港咸潮入侵，
其中最主要的入侵因素是北支咸潮倒灌。
以往对长江河口盐水入侵的研究，主要针对陈行水库盐水入侵来源。
陈行水库于１９９２年建
成启用，运营时间较早，对其研究较多。
青草沙水库
于２０１０年建成启用，运营时间较晚，对盐水入侵来源研究相对较少。
开展青草沙水库盐水入侵来源的研究，定量给出不同潮型下北港外海盐水入侵和上游北支盐水倒灌在水库取水口盐水来源的百分比，可为水库避咸蓄淡提供科技依据。
本文应用长江河口盐水入侵三维数值模式，在模式验证的基础上，开展青草沙水库盐水入侵来源的研究。
图１ 长江河口形势图和测站位置
Ｆｉｇ．１ MａｐｏｆｔｈｅＣｈａｎｇｊｉａｎｇＥｓｔｕａｒｙａｎｄｔｈｅｄｉｓｔｒｉｂｕｔｉｏｎｏｆｏｂｓｅｒｖａｔｉｏｎｓｔａｔｉｏｎｓ
三角形为堡镇水文站、圆点为船测站Ｂ和Ｃ，正方形为浮筒测站Ａ、Ｄ和水库取水口
ＴｈｅｔｒｉａｎｇｌｅｒｅｐｒｅｓｅｎｔｓＢａｏｚｈｅｎｈｙｄｒｏｍｅｔｒｉｃｓｔａｔｉｏｎ，ｂｌａｃｋｄｏｔｓｒｅｐｒｅｓｅｎｔｂｏａｔｏｂｓｅｒｖａｔｉｏｎｓｔａｔｉｏｎｓＢａｎｄＣ，ａｎｄｔｈｅｓｑｕａｒｅｓ
ｒｅｐｒｅｓｅｎｔｂｕｏｙｓｔａｔｉｏｎｓａｎｄｔｈｅｗａｔｅｒｉｎｔａｋｅｓｉｔｅｏｆｔｈｅｒｅｓｅｒｖｏｉｒ
２ 三维数值模式的设置和验证
２·1 模式设置
本文应用改进的三维数值模式ＥＣＯM -ｓｉ［８—１０］
，该模式长期应用于长江河口地区水动力过程和盐水入侵等方面的研究，并取得诸多成果［３，１１—１４］。
模式采用水平曲线非正交网格，范围包括整个长江河口、杭州湾和邻近海区，上游边界设在长江枯季潮区界大通，外海开边界东边到１２４．５°Ｅ附近，
北边到３３°Ｎ附近，南边到２８°Ｎ附近（见图２）。
对长江河口区域，包括南北支分汊口和深水航道工程区域的网格进行局部加密，并且较好地拟合了岸线和导堤。
口内网格分辨率为１００～５００ｍ不等，口外网格较疏，分辨率最大为１０ｋｍ左右。
垂向采用σ
坐标均匀分为１０层，时间步长取６０ｓ。
长江河口区
域浅滩较多，
模式运用干湿判别法实现潮滩移动边界的模拟，临界水深取０．２ｍ。
模式地形采用２０１０
年岸线及水深资料。
外海开边界由潮位驱动，考虑１６个分潮（M ２，Ｓ２，Ｎ２，Ｋ２，Ｋ１，Ｏ１，Ｐ１，Ｑ１，MＵ２，
ＮＵ２，Ｔ２，Ｌ２，２Ｎ２，Ｊ１，M １，ＯＯ１），由各分潮调和常数合成得到。
分潮调和常数由全球潮汐数值模式ＮＡＯＴＩＤＥ计算的结果得到（ｈｔｔｐ：／／ｗｗｗ．ｍｉｚ．
ｎａｏ．ａｃ．ｊｐ／）。
初始水位和流速取零。
关于初始盐度场，从大通到徐六泾６２０ｋｍ的河道，均为淡水，
盐度为０，从徐六泾到口门由枯季各月多次实测资料插值得到，在长江口外由《渤海黄海东海海洋图集
（水文）》各月分布图数字化得到［１５］，考虑海表面风
应力的作用。
２
３１海洋学报 ３６卷
图２ 模式计算区域和网格（ａ），南北支分汊口及其附近水域网格（ｂ），口门深水航道及其邻近水域网格（ｃ）Ｆｉｇ．２ Ｄｏｍａｉｎａｎｄｍｅｓｈｆｏｒｎｕｍｅｒｉｃａｌｍｏｄｅｌ（ａ），ｅｎｌａｒｇｅｄｖｉｅｗｓｏｆｔｈｅｍｏｄｅｌｍｅｓｈａｒｏｕｎｄｔｈｅｂｉｆｕｒｃａｔｉｏｎｏｆｔｈｅＮｏｒｔｈＢｒａｎｃｈａｎｄｔｈｅＳｏｕｔｈＢｒａｎｃｈ（ｂ），ａｎｄｔｈｅＤｅｅｐWａｔｅｒｗａｙＰｒｏｊｅｃｔ（ｃ）
２·２ 模式验证
改进的ＥＣＯM-ｓｉ模式已在长江河口进行了大量的验证，模式计算的水位、流速、流向和盐度与实测资料吻合良好［３，１０，１４］。
本文采用２０１１年１２月２４日至２０１２年１月１３日长江口观测资料对模式作进一步验证，测站的分布图见图１。
模式从２０１１年１１月１日起算，径流量采用大通水文站每日的实测值，风场采用ＱＳＣＡＴ／ＮＣＥＰ提供的时间分辨率为６ｈ、空间分辨率为０．５°×０．５°的数值产品。
图３和图４分别为船测站点Ｂ和Ｃ表底层流速、流向和盐度观测值和模拟计算值随时间变化，测点Ｂ位于南北支分汊口附近，测点Ｃ位于北港北汊口。
测点Ｂ实测最大流速表层约为１．４ｍ／ｓ，底层约为１ｍ／ｓ，因底摩擦作用表层流速大于底层流速，模式计算流速大于实测值。
落潮流历时大于涨潮流历时，流向同实测流向一致。
观测表层盐度小于底层盐度，计算盐度与观测盐度较为接近。
测点Ｃ模式计算流速、流向和盐度与实测值吻合良好，涨憩时刻盐度达到最大，落憩时刻盐度达到最小，表明盐水入侵源自下游外海。
同样，落潮流历时大于涨潮流历时，底层盐度大于表层盐度。
图５为浮筒和水文站观测和模式计算表层盐度随时间变化。
浮筒测点Ａ位于北支口，盐度在２５～３０之间波动，模式计算值略高于实测值。
浮筒测点Ｄ位于北港拦门沙区域，该处存在盐度锋面，盐度随时间变化大，模式计算值略大于实测值。
浮筒测点青草沙水库取水口在大部分时间盐度低于饮用水标准０．４５，在２０１２年１月６－１０日出现大的波动和峰值，模式很好地模拟出了这个过程。
在堡镇水文站，在２０１１年１２月２４至２０１２年１月４日期间，盐度在０．４上下变化，波动很小。
在１月５－１３日之间盐度出现大的波动和峰值，实测盐度最大值达到９．０，同样模式很好地再现了这个过程。
从上面的验证结果可以看出模式计算的流速、流向和盐度同实测资料吻合良好，可以较好地模拟长江河口水动力和盐水入侵过程。
２·3 数值试验
本文研究青草沙水库盐水入侵来源，影响水库取水口盐水入侵的因素有径流量、潮汐和风应力等。
先设计一个径流量和风况为冬季１－２月一般状况的控制数值试验，模拟长江河口和水库取水口盐水入侵状况。
数值模式计算时段为盐水入侵严重的１－２月，模式从１月１日开始运行，至２月底结束，输出２月的计算结果作分析和比较。
径流量取大通１９５０年以来１月和２月平均值，分别为１．１１×１０４ｍ３／ｓ和１．２×１０４ｍ３／ｓ。
风场取ＮＣＥＰ多年半月平均风场，外海开边界考虑１６个主要分潮。
为了解在不同动力条件下青草沙水库盐水入侵状况和来源，设计在北支上段封堵的数值试验，将没有北支盐水倒灌情况下计算结果与控制试验结果比
３３１
１１期 陈泾等：长江河口青草沙水库盐水入侵来源
图３ 船只测点Ｂ处流速（ａ、ｂ）、流向（ｃ、ｄ）和盐度（ｅ、ｆ）
随时间变化Ｆｉｇ．３ Ｖａｒｉａｔｉｏｎｓｏｆｃｕｒｒｅｎｔｓｐｅｅｄ（ａ、ｂ），ｄｉｒｅｃｔｉｏｎ（ｃ、ｄ）ａｎｄｓａｌｉｎｉｔｙ（ｅ、ｆ）ｆｒｏｍ１７：３０ｏｎＪａｎｕａｒｙ１１，２０１２ｔｏ
７：００ｏｎＪａｎｕａｒｙ１３，２０１２ａｔｔｈｅｓｔａｔｉｏｎＢ
ａ、ｃ、ｅ为表层，ｂ、ｄ、ｆ为底层；测定时间：２０１２年１月１１日１７：３０至１月１３日０７：００；
点为实测值，实线为模拟值ａ、ｃ、ｅｉｓｓｕｒｆａｃｅ，ｂ、ｄ、ｆｉｓｂｏｔｔｏｍ；ｄｏｔｓｒｅｐｒｅｓｅｎｔｏｂｓｅｒｖｅｄｄａｔａ，ａｎｄｂｌａｃｋｌｉｎｅｓｒｅｐｒｅｓｅｎｔｍｏｄｅｌ
ｒｅｓｕｌｔｓ
图４ 船只测点Ｃ处流速（ａ、ｂ）、流向（ｃ、ｄ）和盐度（ｅ、ｆ）
随时间变化Ｆｉｇ．４ Ｖａｒｉａｔｉｏｎｓｏｆｃｕｒｒｅｎｔｓｐｅｅｄ（ａ、ｂ），ｄｉｒｅｃｔｉｏｎ（ｃ、ｄ）ａｎｄｓａｌｉｎｉｔｙ（ｅ、ｆ）ｆｒｏｍ１２：３０ｏｎＪａｎｕａｒｙ１０，２０１２ｔｏ
７：
３０ｏｎＪａｎｕａｒｙ１３，２０１２ａｔｔｈｅｓｔａｔｉｏｎＣａ、ｃ、ｅ为表层，ｂ、ｄ、ｆ为底层；测定时间：２０１２年１月１０日１２：００至１月１３日０７：３０；
点为实测值，实线为模拟值ａ、ｃ、ｅｉｓｓｕｒｆａｃｅ，ｂ、ｄ、ｆｉｓｂｏｔｔｏｍ；ｄｏｔｓｒｅｐｒｅｓｅｎｔｏｂｓｅｒｖｅｄｄａｔａ，ａｎｄｂｌａｃｋｌｉｎｅｓｒｅｐｒｅｓｅｎｔｍｏｄｅｌｒｅｓｕｌｔｓ
４
３１海洋学报 ３６卷
图５ 浮筒和水文站测点２０１１年１２月２４日０∶００至２０１２年１月１３日２４∶００表层盐度随时间变化
Ｆｉｇ．５ Ｖａｒｉａｔｉｏｎｓｏｆｓｕｒｆａｃｅｓａｌｉｎｉｔｙｏｆｔｈｅｏｂｓｅｒｖｅｄｄａｔａａｎｄｍｏｄｅｌｒｅｓｕｌｔｓａｔｔｈｅｂｕｏｙｓｔａｔｉｏｎｓａｎｄｈｙｄｒｏｍｅｔｒｉｃｓｔａｔｉｏｎｓｆｒｏｍ００∶００ｏｎＤｅｃｅｍｂｅｒ２４，２０１１ｔｏ２４∶００ｏｎＪａｎｕａｒｙ１３，２０１２．
ａ．浮筒Ａ，ｂ．浮筒Ｂ，ｃ．浮筒青草沙水库取水口，ｄ．堡镇水文站
ａ．ｂｕｏｙｓｔａｔｉｏｎＡ，ｂ．ｂｕｏｙｓｔａｔｉｏｎＢ，ｃ．ｂｕｏｙｓｔａｔｉｏｎａｔｔｈｅｗａｔｅｒｉｎｔａｋｅｓｏｆＱｉｎｇｃａｏｓｈａｒｅｓｅｒｖｏｉｒ，ｄ．Ｂａｏｚｈｅｎｈｙｄｒｏｍｅｔｒｉｃｓｔａｔｉｏｎ
较，定量分析青草沙水库盐水入侵来源。
3 结果与分析
3．1 控制试验
图６为青草沙水库取水口水位、流速、流向和表层、底层盐度随时间变化，潮汐为半日潮，日不等现象明显，半月大小潮变化显著。
大潮最高潮位约２．１ｍ，最低潮位约－１．１ｍ（黄海８５基面）。
水流为往复流、表层最大流速约１３０ｃｍ／ｓ。
取水口大部分时间盐度低于０．４５，２月２４日至３月２日盐度大于０．４５，水库不宜取水。
在２月４－８日小潮后中潮期间底层盐度出现大的峰值，最大值达到４．８，量值远大于表层盐度。
峰值出现在涨憩时刻，表明取水口底层盐水入侵来自北港外海。
除了小潮后中潮期间表底层盐度差异显著，其他时段差异微小，表明垂向混合均匀。
为清楚地表示不同潮型期间取水口盐水来源，将２月１１－２６日表层流向和盐度放大（见图７）。
在大潮后中潮期间１１－１４日（结合图６中水位过程），前２天落憩时刻盐度达到峰值，涨憩时刻达到谷值，表明盐水入侵来自上游，即北支盐水倒灌。
后２天盐度大于０．４５，但波动很小，表明倒灌盐水团正处于取水口水域，水平混合均匀。
在小潮期间１５－１８日，绝大部
５３１
１１期 陈泾等：长江河口青草沙水库盐水入侵来源
分盐度小于０．４５，水库能取淡水，盐度波动小，但峰值出现在涨憩时刻，表明盐水入侵来自下游外海。
但该低盐水源自北支倒灌、小潮期间移至取水口下游，故本质上还是源自上游北支倒灌。
在小潮后中潮期间１９－２２日，盐度小于０．４５，盐度峰值出现在涨憩，盐水入侵来自下游外海。
在大潮期间２３－２６日，盐度随时间逐渐升高，盐度峰值出现在落憩时刻，盐度入侵来自上游北支倒灌。
因此，大潮和大潮后中潮期，水库取水口盐水入侵来自上游北支盐水倒灌，小潮和小潮后中潮期来自下游外海盐水入侵。
图６ 控制数值试验计算的青草沙水库取水口２月１日至３月２日水位、表层流速和流向、盐度随时间变化Ｆｉｇ．６ Ｔｅｍｐｏｒａｌｖａｒｉａｔｉｏｎｏｆｗａｔｅｒｌｅｖｅｌ，ｓｕｒｆａｃｅｃｕｒｒｅｎｔｓｐｅｅｄａｎｄｄｉｒｅｃｔｉｏｎ，ａｎｄｓａｌｉｎｉｔｙｆｒｏｍＦｅｂｒｕａｒｙ１ｔｏMａｒｃｈ２ａｔ
ｔｈｅｗａｔｅｒｉｎｔａｋｅｏｆＱｉｎｇｃａｏｓｈａＲｅｓｅｒｖｏｉｒ，ｃａｌｃｕｌａｔｅｄｂｙｔｈｅｃｏｎｔｒｏｌｅｘｐｅｒｉｍｅｎｔ
表层黑色实线，底层红色虚线，黑色虚线为盐度０．４５；ｔ１和ｔ２为小潮后中潮期间涨憩和落憩时刻，ｔ３和ｔ４为大潮期间涨憩和落憩
时刻，对应下面盐度平面分布输出时刻
Ｂｌａｃｋｌｉｎｅａｎｄｒｅｄｄａｓｈｅｄｌｉｎｅｒｅｐｒｅｓｅｎｔｓｓｕｒｆａｃｅｄａｔａａｎｄｂｏｔｔｏｍｄａｔａ，ｒｅｓｐｅｃｔｉｖｅｌｙ，ａｎｄｔｈｅｓａｌｉｎｉｔｙｒｅｐｒｅｓｅｎｔｅｄｂｙｔｈｅｂｌａｃｋｄａｓｈｅｄｌｉｎｅｉｓ
０．４５．ｔ１ａｎｄｔ２ｒｅｐｒｅｓｅｎｔｆｌｏｏｄｓｌａｃｋａｎｄｅｂｂｓｌａｃｋｉｎｍｏｄｅｒａｔｅｔｉｄｅｆｏｌｌｏｗｉｎｇｎｅａｐｔｉｄｅ，ｒｅｓｐｅｃｔｉｖｅｌｙ，ａｎｄｔ３ａｎｄｔ４ｒｅｐｒｅｓｅｎｔｆｌｏｏｄｓｌａｃｋ
ａｎｄｅｂｂｓｌａｃｋｉｎｓｐｒｉｎｇｔｉｄｅ，ｒｅｓｐｅｃｔｉｖｅｌｙ，ａｎｄｔｈｅｓｅｆｏｕｒｔｉｍｅｐｏｉｎｔｓｃｏｒｒｅｓｐｏｎｄｔｏｔｈｅｏｕｔｐｕｔｔｉｍｅｏｆｔｈｅｆｏｌｌｏｗｉｎｇｈｏｒｉｚｏｎｔａｌｓａｌｉｎｉｔｙｄｉｓｔｒｉ-
ｂｕｔｉｏｎｉｎｔｈｅｆｉｇｕｒｅｓ
２月４－７日小潮后中潮期间涨憩时刻（ｔ１）、落憩
时刻（ｔ２）表层和底层盐度平面图可见图８，在青草沙
水库附近北港水域盐度从口外向口内递减，北港盐度
锋面明显。
在涨憩时刻，表层等盐度线０．４５正在经
过取水口，上游为盐度低于０．４５的淡水，下游为盐
水，底层北港北汊高盐水入侵北港十分明显，盐度显
著升高，取水口下游附近出现等盐度线３，表明时刻
底层盐水入侵强烈，这与图５给出的小潮后中期期间
底层高盐度峰值是一致的。
在落憩时刻，等盐度线
０．４５向下游移动，取水口附近淡水面积增大，上游倒
灌的盐水距取水口尚远。
故从取水口涨憩时刻盐度
增大、淡水面积减少，落憩时刻盐度减小、淡水面积增
加的盐度平面分布，可得出小潮后中潮期间取水口盐
水入侵来自下游外海。
６３１海洋学报 ３６卷
图７ 局部时段放大的青草沙取水口表层流向和盐度随时间变化（２月１１－２６日）Ｆｉｇ．７ Ａｚｏｏｍｅｄｖｉｅｗｏｆｔｅｍｐｏｒａｌｖａｒｉａｔｉｏｎｏｆｓｕｒｆａｃｅｓａｌｉｎｉｔｙａｎｄｃｕｒｒｅｎｔｄｉｒｅｃｔｉｏｎａｔｔｈｅ
ｗａｔｅｒｉｎｔａｋｅｏｆＱｉｎｇｃａｏｓｈａＲｅｓｅｒｖｏｉｒｆｒｏｍＦｅｂｒｕａｒｙ１１ｔｏ
２６
图８ 小潮后中潮期间涨憩时刻（ｔ１）表层（ａ）、底层（ｂ），落憩时刻（ｔ２）表层（ｃ）、底层（ｄ）
盐度平面分布Ｆｉｇ．８ Ｓａｌｉｎｉｔｙｄｉｓｔｒｉｂｕｔｉｏｎｓａｔｓｕｒｆａｃｅ（ａ）ａｎｄｂｏｔｔｏｍ（ｂ）ｌａｙｅｒａｔｆｌｏｏｄｓｌａｃｋ（ｔ１），ａｎｄａｔｓｕｒｆａｃｅ（ｃ）ａｎｄｂｏｔｔｏｍ（ｄ）
ａｔｅｂｂｓｌａｃｋ（ｔ２）ｉｎｍｏｄｅｒａｔｅｔｉｄｅｆｏｌｌｏｗｉｎｇｎｅａｐｔｉｄｅ
红色虚线代表０．４５等盐度线
Ｔｈｅｉｓｏｈａｌｉｎｅｂｙｔｈｅｒｅｄｄａｓｈｅｄｌｉｎｅｓｉｓ０．４５
７
３１１１期 陈泾等：
长江河口青草沙水库盐水入侵来源
在大潮期间２月２５日涨憩和落憩时刻（
图９），北支倒灌进入南支的盐度高于１．０盐水位于南支上段北侧，盐度低于０．４５的淡水位于横沙小港北侧附近水域，在涨憩和落憩时刻取水口盐度均大于０．４５，
不宜取水。
从取水口盐度变化看，涨憩时刻下游淡水区域趋近，落憩时刻下游淡水区域远离，而上游高盐水趋近，所以大潮期间取水口盐水入侵源自上游北支倒灌。
这与上面从盐度变化过线分析结果一致。
图９ 大潮期间涨憩（ｔ３）表层（ａ）、底层（ｂ），落憩（ｔ４）表层（ｃ）、底层（ｄ）
盐度平面分布Ｆｉｇ．９ Ｓａｌｉｎｉｔｙｄｉｓｔｒｉｂｕｔｉｏｎｓａｔｓｕｒｆａｃｅ（ａ）ａｎｄｂｏｔｔｏｍ（ｂ）ｌａｙｅｒａｔｆｌｏｏｄｓｌａｃｋ（ｔ３），ａｎｄａｔｓｕｒｆａｃｅ（ｃ）ａｎｄｂｏｔｔｏｍ（ｄ）
ａｔｅｂｂｓｌａｃｋ（ｔ４）ｄｕｒｉｎｇｓｐｒｉｎｇｔｉｄｅ
3·２ 北支盐水倒灌影响
在控制实验的基础上，封堵北支上段，完全消除北支倒灌的影响，模拟和对比分析青草沙水库盐水入侵。
图１０为封堵北支上段后青草沙水库取水口表层和底层盐度随时间变化过程，对比未封堵时情况（见图６），除了２月４－７日底层盐度仍出现峰值、量值基本一致外，其他时段盐度大幅下降。
这表明取水口表层盐水入侵来自北支盐水倒灌，底层小潮后中潮来自下游外海，其他时段主要来自上游北支盐水倒灌。
为定量分析北支倒灌在青草沙水库取水口附近
水域盐水入侵中的比例，
我们计算了沿北港横断面不同潮型期间涨潮和落潮总的平均盐通量（为绝对值之和平均，见表１），在控制数值试验中小潮后中潮、大潮、大潮后中潮和小潮期间平均盐通量分别为２０．６６、２０．７６、３４．１１和１５．２８ｔ／ｓ，
封堵北支上段后分别变为１２．０８、２．８０、０．１１和０．００ｔ／ｓ。
通过北港横断面的盐通量由来自上游北支盐水倒灌和下游外海盐水入侵共同引起，那么封堵北支上段后盐通量的减少部分可认为是由北支倒灌产生，从整个断面角度看北支倒灌占盐水入侵的百分比在小潮后中潮、大潮、大潮后中
８
３１海洋学报 ３６卷
潮和小潮期间分别为４１．５％、８６．５％、９９．７％和
１００％，
也即除了在小潮后中潮期间青草沙水库取水口盐水入侵５８．５％来自下游外海，
其他潮型主要来自北支盐水倒灌，尤其是大潮后中潮和小潮几乎全部来自北支盐水倒灌。
图１０ 封堵北支上段情况下青草沙水库取水口２月１日至３月２日表层盐度（黑线）、底层盐度（红线）随时间变化Ｆｉｇ．１０ Ｖａｒｉａｔｉｏｎｏｆｓｕｒｆａｃｅ（ｂｌａｃｋｌｉｎｅ）ａｎｄｂｏｔｔｏｍ（ｒｅｄｄａｓｈｅｄｌｉｎｅ）ｓａｌｉｎｉｔｙａｔｔｈｅｗａｔｅｒｉｎｔａｋｅｏｆＱｉｎｇｃａｏｓｈａＲｅｓｅｒｖｏｉｒｆｒｏｍ
Ｆｅｂｒｕａｒｙ１ｔｏMａｒｃｈ２，
ａｆｔｅｒｔｈｅＮｏｒｔｈＢｒａｎｃｈｗａｓｂｌｏｃｋｅｄ 从青草沙水库取水口表层和底层盐度变化看，
小潮后中潮、大潮、大潮后中潮和小潮期间北支倒灌占表层盐水入侵比例分别为６９．５％、８９．３％、９８．５％和
９９．５％，占底层盐水入侵比例分别为３４．９％、８８．９％、９８．５％和９９．５％（
表２），也就是说青草沙水库取水口表层盐水入侵来源主要来自北支盐水倒灌，尤其是大潮后中潮和小潮期间几乎全部来自北支盐水倒灌。
底层盐水入侵来源小潮后中潮期间主要来自下游外海，占６５．１％，其他潮型主要来自北支盐水倒灌，尤其是大潮后中潮和小潮期间几乎全部来自北支盐水倒灌。
长江河口盐水入侵最大特点是北支盐水倒灌，严重影响南支盐水入侵和水源地。
因为青草沙水库取水口距离北港口门较远，北支盐水倒灌成为了除小潮后中潮期间外盐水入侵外的主要来源。
表1 北港断面不同潮型期间平均盐通量（t ／s
）T a b ．1 A v e r a g e d s a l i n i t y
f l u x e s i nd i f f e r e n t t i d e p a t t e r n s a t t h e c r o s s s e c t i o n o f t h eN o r t hC h a n n e l 小潮后中潮
大潮大潮后中潮小潮控制数值试验２０．６６２０．７６３４．１１１５．２８封堵北支试验１２．０８２．８００．１１０．００北支倒灌百分比
４１．５％
８６．５％
９９．７％
１００％
表２ 青草沙水库取水口不同潮型期间平均盐度
T a b ．２ A v e r a g e d s a l i n i t y i nd i f f e r e n t t i d e p a t t e r n s a t t h ew a t e r i n t a k e o f t h e Q i n g
c a o s h aR e s e r v o i r 小潮后中潮大潮
大潮后中潮
小潮
表层
底层表层底层表层底层表层底层控制数值试验０．２０００．３７２０．３４６０．３５００．６０９０．６１４０．３８７０．３８９封堵北支试验０．０６１０．２４２０．０３７０．０３９０．００９０．００９０．００２０．００２北支倒灌百分比
６９．５％
３４．９％
８９．３％
８８．９％
９８．５％
９８．５％
９９．５％
９９．５％
９
３１１１期 陈泾等：
长江河口青草沙水库盐水入侵来源
４ 结论
应用改进的三维数值模式ＥＣＯM-ｓｉ，考虑径流量、潮汐和风应力等的作用，研究青草沙水库取水口盐水入侵的来源。
采用２０１１年１２月至２０１２年１月长江口观测资料对模式作流速、流向和盐度的验证，模式计算结果与实测值吻合良好。
设计控制数值试验，根据盐度和流向的变化过程、涨憩和落憩时刻盐度等值线和淡水区域的变化，分析在一般动力条件下青草沙水库取水口盐水入侵来源。
设计封堵北支上段数值试验，计算不同潮型下通过北港横断面的涨潮和落潮平均盐通量和取水口表层和底层平均盐度，并与控制数值试验结果比较，定量得出盐水入侵来源的百分比。
从通过北港断面盐通量计算结果看，青草沙水库取水口附近水域盐水入侵中北支倒灌所占百分比在小潮后中潮、大潮、大潮后中潮和小潮期间分别为４１．５％、８６．５％、９９．７％和１００％。
从青草沙水库取水口表层和底层盐度变化看，小潮后中潮、大潮、大潮后中潮和小潮期间北支倒灌占表层盐水入侵比例分别为６９．５％、８９．３％、９８．５％和９９．５％，占底层盐水入侵比例分别为３４．９％、８８．９％、９８．５％和９９．５％。
表层盐水入侵来源主要来自北支盐水倒灌，底层盐水入侵来源小潮后中潮期间主要来自下游外海，占６５．１％。
除了小潮后中潮，其他潮型表层和底层盐水入侵来源主要来自北支盐水倒灌，尤其是大潮后中潮和小潮期间几乎全部来自北支盐水倒灌。
需要指出的是本文结果是在一般动力因子情况下模式结算的结果，若冬季出现偏北大风，北港的盐水入侵会大幅加强，期间盐水入侵来源可能主要来自下游外海。
参考文献：
［１］ 乐勤，关许为，刘小梅，等．青草沙水库取水口选址与取水方式研究［Ｊ］．给水排水，２００９（２）：４６－５１．
［２］ 沈焕庭，茅志昌，朱建荣．长江河口盐水入侵［M］．北京：海洋出版社，２００３：１６．
［３］ ＬｉＬ，ＺｈｕＪＲ，WｕＨ．ＩｍｐａｃｔｓｏｆｗｉｎｄｓｔｒｅｓｓｏｎｓａｌｔｗａｔｅｒｉｎｔｒｕｓｉｏｎｉｎｔｈｅＹａｎｇｔｚｅＥｓｔｕａｒｙ［Ｊ］．ＳｃｉｅｎｃｅＣｈｉｎａＥａｒｔｈＳｃｉｅｎｃｅｓ，２０１２，５５（７）：１１７８－１１９２．
［４］ 朱建荣，傅利辉，吴辉．风应力和科氏力对长江河口没冒沙淡水带的影响［Ｊ］．华东师范大学学报：自然科学版，２００８（６）：１－８．
［５］ 项印玉，朱建荣，吴辉．冬季陆架环流对长江河口盐水入侵的影响［Ｊ］．自然科学进展，２００９，１９（２）：１９２－２０２．
［６］ 茅志昌，沈焕庭，肖成献．长江口北支盐水倒灌南支对青草沙水源地的影响［Ｊ］．海洋与湖沼，２００１，３２（１）：５８－６６．
［７］ 顾玉亮，吴守培，乐勤．北支盐水入侵对长江口水源地影响研究［Ｊ］．人民长江，２００３，３４（４）：１－４．
［８］ 朱建荣．海洋数值计算方法和数值模式［M］．北京：海洋出版社，２００３：７２－１２６．
［９］ 朱建荣，朱首贤．ＥＣＯM模式的改进及在长江河口、杭州湾及邻近海区的应用［Ｊ］．海洋与湖沼，２００３，３４（４）：３６４－３７４．
［１０］ WｕＨ，ＺｈｕＪ．Ａｄｖｅｃｔｉｏｎｓｃｈｅｍｅｗｉｔｈ３ｒｄｈｉｇｈ-ｏｒｄｅｒｓｐａｔｉａｌｉｎｔｅｒｐｏｌａｔｉｏｎａｔｔｈｅｍｉｄｄｌｅｔｅｍｐｏｒａｌｌｅｖｅｌａｎｄｉｔｓａｐｐｌｉｃａｔｉｏｎｔｏｓａｌｔｗａｔｅｒｉｎｔｒｕｓｉｏｎｉｎｔｈｅＣｈａｎｇｊｉａｎｇＥｓｔｕａｒｙ［Ｊ］．ＯｃｅａｎMｏｄｅｌｌｉｎｇ，２０１０，３３（１）：３３－５１．
［１１］ WｕＨ，ＺｈｕＪ，ＣｈｅｎＢ，ｅｔａｌ．ＱｕａｎｔｉｔａｔｉｖｅｒｅｌａｔｉｏｎｓｈｉｐｏｆｒｕｎｏｆｆａｎｄｔｉｄｅｔｏｓａｌｔｗａｔｅｒｓｐｉｌｌｉｎｇｏｖｅｒｆｒｏｍｔｈｅＮｏｒｔｈＢｒａｎｃｈｉｎｔｈｅＣｈａｎｇｊｉａｎｇＥｓｔｕａｒ-ｙ：Ａｎｕｍｅｒｉｃａｌｓｔｕｄｙ［Ｊ］．Ｅｓｔｕａｒｉｎｅ，ＣｏａｓｔａｌａｎｄＳｈｅｌｆＳｃｉｅｎｃｅ，２００６，６９（１）：１２５－１３２．
［１２］ ＣｈｅｎＢ，ＺｈｕＪ，ＦｕＬ．ＦｏｒｍａｔｉｏｎｍｅｃｈａｎｉｓｍｏｆｆｒｅｓｈｗａｔｅｒｚｏｎｅａｒｏｕｎｄｔｈｅMｅｉｍａｏＳａｎｄｂａｎｋｉｎｔｈｅＣｈａｎｇｊｉａｎｇｅｓｔｕａｒｙ［Ｊ］．ＣｈｉｎｅｓｅＪｏｕｒｎａｌｏｆＯ-ｃｅａｎｏｌｏｇｙａｎｄＬｉｍｎｏｌｏｇｙ，２０１０，２８：１３２９－１３３９．
［１３］ ＬｉＬ，ＺｈｕＪ，WｕＨ，ｅｔａｌ．ＡｎｕｍｅｒｉｃａｌｓｔｕｄｙｏｎｗａｔｅｒｄｉｖｅｒｓｉｏｎｒａｔｉｏｏｆｔｈｅＣｈａｎｇｊｉａｎｇ（Ｙａｎｇｔｚｅ）ｅｓｔｕａｒｙｉｎｄｒｙｓｅａｓｏｎ［Ｊ］．ＣｈｉｎｅｓｅＪｏｕｒｎａｌｏｆＯ-ｃｅａｎｏｌｏｇｙａｎｄＬｉｍｎｏｌｏｇｙ，２０１０，２８：７００－７１２．
［１４］ WｕＨ，ＺｈｕＪ，ＨｏＣｈｏｉＢ．ＬｉｎｋｓｂｅｔｗｅｅｎｓａｌｔｗａｔｅｒｉｎｔｒｕｓｉｏｎａｎｄｓｕｂｔｉｄａｌｃｉｒｃｕｌａｔｉｏｎｉｎｔｈｅＣｈａｎｇｊｉａｎｇＥｓｔｕａｒｙ：Ａｍｏｄｅｌ-ｇｕｉｄｅｄｓｔｕｄｙ［Ｊ］．Ｃｏｎｔｉ-ｎｅｎｔａｌＳｈｅｌｆＲｅｓｅａｒｃｈ，２０１０，３０（１７）：１８９１－１９０５．
［１５］ 海洋图集编委会．渤海黄海东海海洋图集（水文）［M］．北京：海洋出版社，１９９２：１３－１６８．
０４１海洋学报 ３６卷
S o u r c e s f o r s a l t w a t e r i n t r u s i o n a t t h ew a t e r i n t a k e o f Q i n g
c a o s h a R e s e r v o i r i n t h eC h a n g j i a n g E s t u a r y
ＣｈｅｎＪｉｎｇ１，ＺｈｕＪｉａｎｒｏｎｇ１
（１．S t a t eK e y L a b o r a t o r y o f E s t u a r i n e a n dC o a s t a lR e s e a r c h ，E a s t C h i n aN o r m a lU n i v e r s i t y ，S h a n g
h a i ２０００６２，C h i n a ）A b s t r a c t ：Ｂａｓｅｄｏｎｔｅｍｐｏｒａｌｖａｒｉａｔｉｏｎｏｆｔｈｅｓｉｍｕｌａｔｅｄｓａｌｉｎｉｔｙａｎｄｃｕｒｒｅｎｔｄｉｒｅｃｔｉｏｎ，ｃｈａｎｇｅｏｆｔｈｅｉｓｏｈａｌｉｎｅａｎｄｆｒｅｓｈｗａｔｅｒｚｏｎｅａｔｆｌｏｏｄｓｌａｃｋａｎｄｅｂｂｓｌａｃｋ，ｔｈｉｓｓｔｕｄｙｕｓｅｄｔｈｅｍｏｄｉｆｉｅｄｎｕｍｅｒｉｃａｌ３-ＤｍｏｄｅｌＥＣＯM -ｓｉｔｏａｎａｌｙｚｅｓｏｕｒｃｅｓｏｆｓａｌｔｗａｔｅｒｉｎｔｒｕｓｉｏｎａｔｔｈｅｗａｔｅｒｉｎｔａｋｅｏｆＱｉｎｇｃａｏｓｈａＲｅｓｅｒｖｏｉｒｕｎｄｅｒｎｏｒｍａｌｄｙｎａｍｉｃｃｏｎｄｉｔｉｏｎ．Ｔｈｅｍｏｄｅｌｅｄｒｅｓｕｌｔｓｓｈｏｗｅｄｔｈａｔｓａｌｔｗａｔｅｒ-ｓｐｉｌｌ-ｏｖｅｒｆｒｏｍｔｈｅＮｏｒｔｈＢｒａｎｃｈｉｎｔｏｔｈｅＳｏｕｔｈＢｒａｎｃｈ（ＳＳＯ）ａｃｃｏｕｎｔｓｆｏｒ６９．５％，８９．３％，９８．５％ａｎｄ９９．５％ｏｆｔｈｅｓｕｒｆａｃｅｓａｌｔｗａｔｅｒｉｎｔｒｕｓｉｏｎａｔｔｈｅｗａｔｅｒｉｎｔａｋｅｏｆｔｈｅＱｉｎｇｃａｏｓｈａＲｅｓｅｒ-ｖｏｉｒｄｕｒｉｎｇｍｏｄｅｒａｔｅｔｉｄｅｆｏｌｌｏｗｉｎｇｎｅａｐｔｉｄｅ，ｓｐｒｉｎｇｔｉｄｅ，ｍｏｄｅｒａｔｅｔｉｄｅｆｏｌｌｏｗｉｎｇｓｐｒｉｎｇｔｉｄｅａｎｄｎｅａｐｔｉｄｅ，ｒｅｓｐｅｃ-
ｔｉｖｅｌｙ，
ａｎｄ３４．９％，８８．９％，９８．５％，９９．５％ｏｆｔｈｅｂｏｔｔｏｍｓａｌｔｗａｔｅｒｉｎｔｒｕｓｉｏｎａｔｔｈｅｗａｔｅｒｉｎｔａｋｅｏｆｔｈｅｒｅｓｅｒｖｏｉｒｄｕｒｉｎｇｔｈｅｓａｍｅｆｏｕｒ-ｔｙｐｅｔｉｄｅｓｍｅｎｔｉｏｎｅｄａｂｏｖｅ．Ｅｘｃｅｐｔｆｏｒｔｈｅｂｏｔｔｏｍｓａｌｔｗａｔｅｒｉｎｔｒｕｓｉｏｎｍａｉｎｌｙｆｒｏｍｄｏｗｎ-ｓｔｒｅａｍｏｐｅｎｓｅａ（ａｃｃｏｕｎｔｓｆｏｒ６５．１％）ｄｕｒｉｎｇｍｏｄｅｒａｔｅｔｉｄｅｆｏｌｌｏｗｉｎｇｎｅａｐｔｉｄｅ，ｔｈｅｓｏｕｒｃｅｏｆｔｈｅｓａｌｔｗａｔｅｒｉｎｔｒｕ-ｓｉｏｎａｔｓｕｒｆａｃｅａｎｄｂｏｔｔｏｍｌａｙｅｒａｔｔｈｅｗａｔｅｒｉｎｔａｋｅｏｆｔｈｅｒｅｓｅｒｖｏｉｒｉｓｍａｉｎｌｙｆｒｏｍｔｈｅＳＳＯ，ｅｓｐｅｃｉａｌｌｙｄｕｒｉｎｇｍｏｄ-ｅｒａｔｅｔｉｄｅｆｏｌｌｏｗｉｎｇｓｐｒｉｎｇａｎｄｎｅａｐｔｉｄｅ．
K e y w
o r d s ：ＱｉｎｇｃａｏｓｈａＲｅｓｅｒｖｏｉｒ；ｓａｌｔｗａｔｅｒｉｎｔｒｕｓｉｏｎ；ｓｏｕｒｃｅ；ｎｕｍｅｒｉｃａｌｓｉｍｕｌａｔｉｏｎ１４１１１期 陈泾等：长江河口青草沙水库盐水入侵来源
长江河口青草沙水库盐水入侵来源
作者：陈泾， 朱建荣， Chen Jing， Zhu Jianrong
作者单位：华东师范大学 河口海岸学国家重点实验室,上海,200062刊名：
海洋学报（中文版）
英文刊名：Acta Oceanologica Sinica
年，卷(期)：2014(11)
本文链接：/Periodical_hyxb201411015.aspx。
