1. 引言
曹妃甸地区的地势低洼,地面标高较低,这使得该地区容易受到海洋的影响。在晚更新世晚期,滦河的冲洪积扇逐渐形成了一个三角洲,河流和海洋的共同作用导致了大量的沉积物在这一区域堆积。这个地区是由河流和海洋共同作用形成的冲海积低平原区。曹妃甸填海造陆区位于渤海湾西北部,是我国规模最大的填海造陆工程 [1] [2] [3] 。曹妃甸岛所处的区域属于暖温带半湿润滨海大陆性季风气候区,四季分明。该地区的河流主要有纳潮河、电力郎道河、电厂排水渠、石化河和通岛河等。曹妃甸海城潮汐属于不规则的半日潮,潮流速度较小,属于弱流区 [4] [5] 。褚宏宪等通过海洋地球物理方法发现曹妃甸沙岛建成后海底地形地貌发生了较大改变,进而影响该海域海洋水动力变化 [6] 。总体上该海区海洋水文实测资料较少,人们缺乏对该海域潮汐潮流动力状况的最新认识,因此开展该海域的现状水动力调查,对潮汐、潮流进行统计分析具有重要意义,为曹妃甸港的规划、运营维护和未来发展提供科学依据。
2. 资料来源
数据主要来源于天津市地质研究和海洋地质中心于2023年春季在唐山市曹妃甸近海海域的水文调查项目。项目共布设了潮位观测站1个和海流观测站3个,观测时间及采用的仪器见表1,观测位置见图1。
本次海洋水文调查观测内容为潮汐和海流特征。观测指标为水位、大小潮26小时连续同步流速、流向要素。
3. 数据处理计算方法
参考文献 [7] [8] [9] [10] 方法﹐分别对潮位、各站各层海流资料进行准调和分析,应用差比关系法计算太阴主要半日分潮M2、太阳主要半日分潮S2、太阴赤纬半日分潮O1、太阳赤纬半日分潮Kl、太阴浅水1/4分潮M4、太阴-太阳1/4分潮MS4。根据《海港水文规范》(JTS145-2-2013) [11] 潮汐性质可按下式计算标准判别:
Table 1. Information table of hydrological observation in Caofeidian sea area
表1. 曹妃甸海域水文观测信息表
Figure 1. The distribution of hydrological observation station in Caofeidian sea area
图1. 曹妃甸海域水文观测站位图
式子中式中的HO1、HK1、HM2、HS2分别为主太阴日分潮、太阴太阳赤纬日分潮、主太阴半日分潮和主太阳半日分潮的平均振幅。
当F ≤ 0.5时为正规半日潮;
当0.5 < F ≤ 2.0时为不正规半日潮;
当2.0 < F ≤ 4.0时为不正规全日潮;
当4.0 < F时为正规全日潮;
潮流按其性质可分为正规的、非正规的半日潮流或全日潮流。
判别标准为全日潮流振幅之和(WO1 + Wk1)与主太阴半日分潮流振幅(WM2)之比K值:
K ≤ 0.5为正规半日潮流;
0.5 < K ≤ 2.0为非正规半日潮流;
2.0 < K ≤ 4.0为非正规全日潮流;
4.0 < K正规全日潮流。
4. 结果分析
4.1. 潮汐特征
根据实测潮位数据(图2)进行统计分析,获得潮位特征值(表2),实测期间潮汐示性系数0.62,属不正规半日潮,最高潮位为0.98 m,最低潮位为−1.52 m,平均高潮位及平均低潮位分别为0.78 m和−0.68 m,一周平均潮位为0.07 m,实测期间最大潮差达到2.30 m,最小潮差为0.53 m,平均潮差为1.80 m,涨落潮历时验潮站表现为落潮历时略大于涨潮历时的特征,平均落潮历时长于平均涨潮历时,差值为1小时6分钟。
Figure 2. Tidal course diagram of tidal station in half a month
图2. 半个月内潮位站潮位过程线图
Table 2. Table of tidal level characteristics of short-term tidal stations
表2. 短期验潮站潮位特征表
4.2. 潮流特征
4.2.1. 实测潮流分析
为分析测区实测流况的基本特征,根据测验成果中的《潮流观测记录报表》,进行了各测站分层最大涨、落流速(流向)、垂线平均最大流速(流向)以及垂线平均流速出现频率和涨、落潮流历时等各项统计。大小潮汛期各测站垂线平均潮流见图3和图4,潮流与潮位关系见图5。
Figure 3. Distribution of average current vertically at observed stations during springtides
图3. 大潮垂线平均流矢图
Figure 4. Distribution of average current vertically at observed stations during neap tide
图4. 小潮垂线平均流矢图
流速上分析:大潮汛,测区海域各测站垂线平均最大落潮流均大于涨潮流,落、涨潮差值在0.03~0.19 m/s之间。小潮汛与大潮汛有所不同,S1、S3测站表现为最大涨潮流大于落潮流,涨、落潮差值在0.05~0.14 m/s之间;S2测站则表现为最大落潮流大于涨潮流,落、涨潮差值为0.04 m/s。测区流速平面分布有以下规律,涨、落潮最大流速均出现在S3测站,其次为S2测站,相对较小的流速则出现在S1测站。各测站垂线平均流速仅S3测站出现1.0节流速以上的测次,频率达44.23%,各测站垂线平均流速出现0.5节流速以上的测次,出现频率最大的也为S3测站,达76.92%,其次为S2测站为38.46%,出现频率最小为是S1测站,仅为11.54%。
流向上分析:S1测站,流向频次出现在16方位上的频率较为分散,在各方位上的频率为1.92%~17.31%。S2测站,出现频率较高的为NNE向,频率为23.08%,其余方位上频率都较小,说明S2测站流向也较为分散。相较S1和S2测站,S3测站流向较为集中,涨潮方向大多集中在WNW~NW向,频率为37.80%;落潮方向大多集中在ESE方向,频率为36.54%。垂线平均流速对应的流向分布上,总体上看,测区各测站潮流运动方式基本都呈往复流的特征。
Figure 5. Relationship of average current vertically and tide at observed stations during spring tides
图5. 测站大潮垂线平均潮流与潮位关系图
潮位与潮流之间的运动关系,是海洋中潮汐升降与潮流涨落之间的相位关系。这一关系,通常用来判别潮波的属性,亦可用来掌握潮流的特征时刻(涨急、涨憩、落急、落憩)与潮位的特征时刻(低潮、高潮、半潮面(低潮至高潮、或高潮至低潮))之间的相位匹配情况。若转流时间(垂线平均流速接近或为零)发生在高潮与低潮的中间时刻,则潮流运动具有前进波性质,若转流时间发生在高潮和低潮时,则潮流运动具有驻波性质。以垂线S3站大潮为例,通过潮位、垂线平均流速。由此可以看出:潮流的“转流”时刻与潮位过程的高、低潮时刻相差1~2小时;潮流的涨急、落急时刻也与潮位过程的半潮面时刻有时间差;这两项基本特征说明测区海域的潮位、潮流运动关系呈现一定的前进波特征。
4.2.2. 潮流调和分析
1) 潮流类型
潮流类型由主太阴日分潮流(O1)与太阴太阳赤纬日分潮流(K1)的椭圆长半轴之和与主太阴半日分潮流(M2)的椭圆长半轴之比值,即(WK1 + WO1)/WM2来确定。同时,为了考察测区浅海分流的大小与作用,往往又将四分之一日主要浅海分潮流M4与主要半日分潮流M2的椭圆长半轴之比作为判据,进行分析。各测站垂线平均海流类型判据(WK1 + WO1)/WM2和WM4/WM2的计算结果如表3中所列。
Table 3. Statistical table of tidal current type criterion for the average vertical line of each station
表3. 各测站垂线平均的潮流类型判据统计表
由表可知:各测站垂线平均潮流类型判据(WK1 + WO1)/WM2之比,在0.15~0.38之间,小于0.50,故半日潮流的特征明显;其中S1、S2测站WM4/WM2之比大于0.04,为0.08,因此,严格地说S1和S2测站潮流性质为非正规半日潮海流的类型,浅水效应较为显著;S3测站WM4/WM2之比为0.02,小于0.04,因此S3测站属正规半日潮海流类型。
2) 潮流类型
潮流运动形式一般可分为旋转流和往复流两种。在半日潮流占主导地位的测区,潮流运动可用M2分潮流的椭圆率K值来表述,一般认为当K值大于0.25时,潮流表现出较强的旋转性,潮流运动可定义为旋转流,当K值小于0.25时,潮流方向主要集中在涨、落潮流的两个方向上,潮流运动可定义为往复流。其中,潮流的旋转方向则以K值的正负来表征,正值为逆时针的左旋,负值为顺时针的右旋。
Table 4. Statistical table of the average M2 subcurrent K value of the vertical lines of each measuring station
表4. 各测站垂线平均的M2分潮流K值的统计表
由上表4可知:S1和S2测站垂线平均的M2分潮流椭圆率(K值)均大于0.25,表明测区海域潮流运动形式具有显著的旋转流特征;S3测站垂线平均的M2分潮流椭圆率(K值)为0.00,S3测站潮流运动形式具有显著的往复流特征。
3) 天文潮可能最大流速
在潮流椭圆要素计算的基础上按《海港水文规范》 [11] 中的有关公式,对水文测验资料各测站的表层、0.6H、底层及垂线平均可能最大流速进行计算;现将天文潮可能最大流速计算结果,统一列于表5中以供比较。
Table 5. Statistical table of the calculation results of the possible maximum velocity of astronomical tide at each station
表5. 各测站天文潮可能最大流速计算结果的统计表
由表可知:各测站主要层及垂线平均的可能最大流速均表现为涨潮流的方向,从而印证了整个测区涨潮流速大于落潮流速的特征;其次,在计算的天文潮可能最大流速均大于实测最大流速,在流向上与实测最大流速的方向较为一致,在分布上亦与实测流况相似,从而表明了计算结果的合理性。
4) 余流特征
余流乃指消除周期性潮流后的一种相对稳定的流动。然而由于受分析方法和计算资料序列的限制,表6列出的余流值可能包含了部分未被分离的潮流成份,但仍可由此获取某些统计性的认识,余流的大小及方向对物质输送的研究具有重要意义。大、小潮余流矢量见图6。
Table 6. The residual current at observed stations
表6. 各站余流及流向统计表
水文测验期间,各测站主要层次最大余流为0.13 m/s,发生在小潮汛的D2测垂线平均最大余流为0.10 m/s,发生在大潮汛D3测站。各测站余流的量值在垂向分布上基本与实测流一致,表现为表层较大,0.6H层次之,底层较小的特征。
Figure 6. Residual vector diagram of large and small tides
图6. 大、小潮余流矢量图
5. 结论
通过对2023年春季曹妃甸近海海域1个站位的潮位资料和3个站位的海流资料进行了统计分析,研究发现该海域海洋水文的一些特征,具体表现为:
1) 调查海域潮位实测期间潮汐示性系数0.62,属不正规半日潮,最高潮位为0.98 m,最低潮位为−1.52 m,平均高潮位及平均低潮位分别为0.78 m和−0.68 m,一周平均潮位为0.07 m,实测期间最大潮差达到2.30 m,最小潮差为0.53 m,平均潮差为1.80 m,涨落潮历时表现为落潮历时略大于涨潮历时的特征,平均落潮历时长于平均涨潮历时,差值为1小时6分钟。
2) 最大涨潮流速极值为0.89 m/s,出现在S3测站大潮汛;最大落潮流速极值为1.12 m/s,亦出现在S3测站大潮汛。各测站的最大涨落潮流速,随着潮汛的演变,总体上具有良好的变化规律,即大潮汛较大,小潮汛较小。
3) 各测站垂线平均潮流类型在0.15~0.38之间,半日潮流的特征明显。其中S1、S2测站潮流性质为非正规半日潮海流的类型,浅水效应较为显著,测区海域潮流运动形式具有显著的旋转流特征;S3测站属正规半日潮海流类型,具有显著的往复流特征。
水文测验期间,天文潮可能最大流速均大于实测最大流速,在流向上与实测最大流速的方向较为一致,在分布上亦与实测流况相似。各测站余流的量值在垂向分布上基本与实测流一致,表现为表层较大,中层次之,底层较小的特征。