1. 引言
葫芦岛市位于辽东湾西南部沿海地区,北依松岭山脉,东邻锦州,西接山海关,南临渤海辽东湾,近海水域和海岸带空间广阔,海洋资源丰富,海洋经济一直是全市国民经济的重要组成部分。海岸线长度为303.49 km,包括自然岸线104.65 km、人工岸线198.84 km。其中,天角山海岸带是葫芦岛港以东,其东部和西部分别分布有区域稀缺的自然砂质海滩,中部的基岩岸线和砾石滩是区域重要的自然岸线资源;此外,也是区域十分重要的山地–基岩岸线–海岸湿地–滩涂生态系统 [1] 。该海域为一个半封闭的浅水海湾,湾口朝向东,湾口水深达10 m以上,平均水深3.5 m [2] 。
近年来,辽宁海岸带围填海工程项目不断开发,受周边开发、人为破坏、海浪侵蚀等因素影响严重,人工痕迹愈加明显,面临外来的威胁压力逐渐加大,沙滩、滨海湿地、基岩岸线、滩涂等生态功能退化明显,并有进一步恶化趋势,亟需保护和修复;在带来经济效益的同时,海湾面积大量减少,1991~2015年各类围填海面积共增加约50.67 km2,长期以来且工程区岸线侵蚀情况严重,道路、防浪墙、自然山体等破损明显,防灾能力亟需提高。人工围填海改变了海湾的潮流系统,海底淤积冲刷改变剧烈,生态功能退化,环境恶化、资源衰退,造成海湾资源严重破坏 [3] 。近年来我国一些地区开展了以恢复滨海湿地为核心目标的海岸带生态修复工程 [4] 。通过因地制宜开展海岸带保护修复工程,在遵循自然规律和生态系统特征的情况下,宜滩则滩、宜荒则荒,促进海岸带区域生态、减灾协同增效。
因此开展天角山海域的潮汐潮流泥沙特征研究对协调资源的开发利用和环境保护尤为重要。该海区最近的实测资料分析尚为2010年前的实测资料 [5] [6] [7] ,人们缺乏对该海域潮汐潮流动力状况的最新认识。因此,本文拟利用2020年冬季的最新观测资料对葫芦岛天角山近海海域的潮汐、潮流及泥沙状况进行数据统计分析,更新人们对葫芦岛天角山近海海域的水文环境状况的认识。这对在该海域开展近海工程、航道利用、海域使用、水产养殖和环境保护等具有重要的科学价值和现实意义。
2. 资料来源
数据主要来源于天津市地质研究和海洋地质中心于2020年冬季在葫芦岛天角山近海海域的水文调查项目。项目共布设了潮位观测站1个和海流观测站9个,观测时间及采用的仪器见表1,观测位置见图1。
![](Images/Table_Tmp.jpg)
Table 1. Information of hydrological observation in Tianjiao Mountain sea area
表1. 天角山海域水文观测信息表
本次海洋水文调查观测内容为潮汐、海流及悬浮泥沙特征。观测指标为水位、26小时连续同步流速、流向要素,并在整点时刻进行悬沙水样采集。
![](//html.hanspub.org/file/6-2830236x8_hanspub.png?20230907095239553)
Figure 1. The distribution of hydrological observation station in Tianjiao Mountain sea area
图1. 天角山海域水文观测站位图
3. 数据处理计算方法
参考文献 [8] 方法﹐分别对潮位、各站各层海流资料进行准调和分析,应用差比关系法计算太阴主要半日分潮M2、太阳主要半日分潮S2、太阴赤纬半日分潮O1、太阳赤纬半日分潮K1、太阴浅水1/4分潮M4、太阴–太阳1/4分潮MS4。6个分潮的潮位、潮流调和常数、椭圆要素和余流。
式子中式中的
、
、
、
分别为主太阴日分潮、太阴太阳赤纬日分潮、主太阴半日分潮和主太阳半日分潮。
根据《海港水文规范》(JTs145-2-2013) [9] 潮汐性质可按下式计算标准判别:
(1)
式子中式中的
、
、
、
分别为主太阴日分潮、太阴太阳赤纬日分潮、主太阴半日分潮和主太阳半日分潮的平均振幅。
当F ≤ 0.5时为正规半日潮;
当0.5 < F ≤ 2.0时为不正规半日潮;
当2.0 < F ≤ 4.0时为不正规全日潮;
当4.0 < F时为正规全日潮;
潮流按其性质可分为正规的、非正规的半日潮流或全日潮流。
判别标准为全日潮流振幅之和(
)与主太阴半日分潮流振幅(
)之比K值:
K ≤ 0.5为正规半日潮流;
0.5 < K ≤ 2.0为非正规半日潮流;
2.0 < K ≤ 4.0为非正规全日潮流;
4.0 < K正规全日潮流;
悬浮泥沙采用重量法进行测试分析,计算公式如下所示:
(2)
式中:SPM为悬浮体浓度;
为带样品滤膜质量;
为水样滤膜质量;
为空白校正滤膜校正值;V为过滤样品水的体积。
4. 结果分析
4.1. 潮汐特征
根据实测潮位数据进行统计分析,获得潮位特征值(表2),最高潮位为5.21 m,最低潮位为0.67 m,平均高潮位为4.12 m,平均低潮位为1.92 m,平均海平面为3.05 m,最大潮差4.21 m,最小潮差为1.22 m,平均潮差2.01 m,平均涨潮历时为6.04 h,平均落潮历时为6.18 h。根据依据葫芦岛验潮站1969~1988年实测潮位资料,统计出本站潮汐特征值与本次调查结果比较而言,平均潮差有所减小但相差不大,最大潮差相对增大,最小潮差相对减小。涨落潮历时相差不大。随着地形和岸线的变化,天角山近海海域的潮汐特征有所变化。
![](Images/Table_Tmp.jpg)
Table 2. Characteristics value of short-term tide station
表2. 短期验潮站潮位特征表
另外,采用最小二乘法对调查的潮位数据进行调和分析 [10] [11] ,计算获得12个天文潮分潮的潮汐调和常数振幅和迟角,结果见表3。
![](Images/Table_Tmp.jpg)
Table 3. Tide harmonic calculation of short-term tide station
表3. 短期验潮站潮汐调和计算表
采用公式(1)计算F值,潮汐示型系数为0.66,潮型系数大于0.5且小于2.0,因此天角山海域的潮汐类型属于不正规半日混合潮。
4.2. 潮流特征
4.2.1. 潮流特征分析
根据本次调查资料的统计结果,为直观分析海流特征,由原始数据提取得到了各站各层(底、中、表)的海流矢量图(图2)。
通过该海区实测海流矢量图,可以得到如下特征。由图2不难发现该海域外海的3个海洋水文站(S3、S6和S9)流速明显大于湾内6个海洋水文站位,由海到陆流速有逐渐降低的趋势。另外该海区海流具有明显的往复流特征,流向主要受到岸线、港口和地形等影响。在外海开阔海域,涨潮流向基本一致,流向主要呈N-E走势;受岸线及堤坝的阻挡和导流的作用的影响,落潮时S3站位部分落潮流流入湾内,S9主要落潮流主要来源于湾内。在湾内,涨潮流整体呈现为N-W走势,落潮流主要呈现为S-E走势。受到防波堤挑流影响,S2涨潮流以N向为主,落潮流以S向为主。可见落潮过程中外海潮流部分潮水由湾内北侧海域进入锦州湾口,锦州湾内水流从南侧湾口流出天角山湾内。
根据天角山验潮站在水文观测期间的水位变化,为进一步说明本区潮流的时空变化,绘成潮位–潮流关系图(图3),从图中可以看出涨潮时刻潮流线均大于水位线,落潮时潮流线均小于水位线。在平潮前后时刻流速接近于零,随着涨落潮开始(水位线切线为零)流速逐渐增加,至涨落潮最大(水位线切线最大)潮流也达到最大值,然后减小,到潮位线达到最大(小),涨落潮停止,潮流也接近于零。
![](//html.hanspub.org/file/6-2830236x24_hanspub.png?20230907095239553)
Figure 2. Distribution of average current vertically at observed stations during springtides
图2. 大潮期观测站垂线平均海流矢量图
![](//html.hanspub.org/file/6-2830236x25_hanspub.png?20230907095239553)
Figure 3. Relationship of average current vertically and tide at observed stations during spring tides
图3. 大潮期观测站垂线平均潮位–潮流关系图
4.2.2. 潮流调和分析
潮流调和分析的目的在于推算分潮的调和常数及分潮的椭圆要素,用以分析潮流性质。分析结果表明,主太阴半日分潮流M2是本海区的优势分潮流。因此,各测站M2分潮流的椭圆长轴走向决定了本海区潮流的主流向。根据公式计算得到潮流示型系数,结果见表4,各站层潮流类型判别系数在均小于0.5。此该海域总体为正规半日潮流区。
![](Images/Table_Tmp.jpg)
Table 4. The type parameter of tidal current at observed stations
表4. 观测站潮流类型系数
潮流的运动形式分旋转流和往复流,通常以椭圆率K的绝对值大小来判断,当
时,潮流椭圆成圆形,各方向流速相等,为纯旋转流;当
时,潮流椭圆为一直线,海水在直线上往返流动,为典型往复流。K值通常在0~1之间,
值越大,旋转流的形式越显著,
值越小,往复流的形式越显著 [12] 。分析结果表明,主太阴半日分潮流M2是本海区的优势分潮流,各测站M2分潮流的椭圆长轴走向决定了本海区潮流的主流向。以M2分潮流的椭圆率来对潮流运动形式作近似分析,各分潮流椭圆率计算结果见表5。
![](Images/Table_Tmp.jpg)
Table 5. Statistical table of sub-tide M2 ellipticity at observed stations
表5. 观测站M2分潮流椭圆率统计表
由表5可知,M2分潮流的
值在0.01~0.64之间,跨度较大,其中S2、S5站位M2均在0.5左右,表现为一定的旋转流,主要原因为受到湾内地形地貌的影响,其余水文站分潮K值较小(0.15左右),主要表现往复流。故因此总体来看该海域为明显的往复流特征。
4.2.3. 余流特征
该海区余流包括冲淡水流及风海流,也包括潮汐引起的长周期或定常的流动。由表6可见,大潮期最大余流流速为12.1 cm/s,最小余流流速为1.4 cm/s。余流较大站位是S3、S6和S9站位,余流最小的站位为S7,可见外海的余流流速明显高于湾内其他站位。余流流向方向不难发现除S7站位外,其他站位流向为偏北方向,主要原因天角山海域冬季盛行北风有关。
![](Images/Table_Tmp.jpg)
Table 6. The residual current at observed stations during the spring tide
表6. 大潮期各站余流及流向统计表
4.3. 悬沙特征
大潮期剖面平均悬沙含量为22.9 mg/dm3 (14.1~31.2 mg/dm3),其中S2站和S3站位悬沙浓度较高,S1站和S7站浓度较低,悬沙浓度见图4。
![](//html.hanspub.org/file/6-2830236x31_hanspub.png?20230907095239553)
Figure 4. Distribution of the average sediment content vertically at observed stations during the spring tide
图4. 各站位大潮沙含量垂向平均分布图
空间分布上各水层含沙浓度的分布具有由近岸到离岸逐渐升高的趋势,近岸站位S4是潮湖通道,流速较大,携沙量较强;S1、S5、S7和S8站位属于港湾半封闭海区内,水动力较弱,含沙量较低;S3、S6和S9站位属于远海站位,水动力较强,含沙量较大;S2属于航道区,船只活动活跃区,含沙量高基本属于人类活动影响作用。
5. 结论
通过2020年冬季天角山近海海域1个站位的潮位资料和9个站位的海流泥沙资料进行了统计分析,研究发现因人类活动海洋水文特征发生了一定的改变,具体表现为:
1) 调查海域潮汐性质属于不正规半潮(0.25 < F = 0.48 < 1.5),近年来围填海工程,由于地形和岸线的变化造成天角山近海海域的潮汐特征有所变化,平均潮差相对减小,最大潮差相对增大,最小潮差相对减小。
2) 大潮期实测最大流速分别为69 cm/s。垂向上流速从表层到底层流速逐渐减小,空间分布上外海流速明显大于湾内。湾内张潮流为西北向,落潮流为西南流,外海涨潮流为东北向,落潮流为西南向,流速有所增大,主要原因为湾内本身岸线及堤坝的阻挡和导流的作用的影响。
3) 调查海域上为正规半日潮流,潮流运动形式总体表现为往复流;各站、层余流流速皆较小,绝大部分余流均小于10 cm/s,空间分布上外海余流大于湾内,主要原因为冬季北向季风影响。
4) 大潮期剖面平均含沙量为22.9 mg/dm3,悬沙浓度空间分布上存在由近岸到离岸逐渐升高的趋势。悬沙浓度和人类活动及潮流大小存在一定相关性。
NOTES
*第一作者。