劉鵬霞， 劉子洲， 宋 新， 顧艷鎮(zhèn)， 翟方國， 李培良

(1.中國海洋大學海洋與大氣學院，山東 青島 266100； 2.山東省濰坊市海洋環(huán)境監(jiān)測中心站，山東 濰坊 261100； 3.浙江大學海洋學院，浙江 舟山 316021)

海洋牧場指在一定的海域內(nèi)，采用規(guī)模化漁業(yè)設施和系統(tǒng)化管理，對魚、蝦、藻、貝等海洋生物資源進行有目的、有計劃的海上放養(yǎng)。海洋牧場對于推進“藍色糧倉”建設有重要意義，是提升經(jīng)濟社會效益的重要途徑[1]。目前及以后一段時期是我國海洋漁業(yè)乃至整個海洋產(chǎn)業(yè)轉型升級的關鍵階段。魯海豐海洋牧場位于青島市黃島區(qū)，屬南黃海西北海域。探究魯海豐海洋牧場的流場分布及變化特征，實現(xiàn)對海洋牧場的可測、可控化管理，對提高我國海洋牧場的漁情預測和預報減災能力具有指導意義。

魯海豐海洋牧場區(qū)域的流場受到鄰近海域甚至整個南黃海的動力環(huán)境影響，針對該海域流系特征，前人已經(jīng)做過不少研究：趙保仁[2]根據(jù)青島近海某浮標站的海流觀測數(shù)據(jù)分析了局地風對當?shù)睾Ａ鞯挠绊?，指出余流隨風變化的反應敏感。湯毓祥等[3]采用一種二維非線性潮汐模式對南黃海的潮流和潮余流進行模擬，指出該海域M2潮流的運動規(guī)律受到旋轉潮波系統(tǒng)的控制，在山東半島南岸近海和膠州灣一帶，潮流為左旋。湯毓祥等[4]根據(jù)中韓黃海水循環(huán)動力學合作調(diào)查所獲資料指出，初春時節(jié)來自山東沿岸的低鹽水，自成山頭近海區(qū)域大致沿50 m等深線向南部延伸，至青島外海在36°5′N, 122°E形成低溫中心，即青島冷水團。Tang等[5]根據(jù)2003年青島外海的實測海流資料分析指出，在該海域潮流作用顯著，其中半日潮流占主，全日潮流次之，潮流垂向均勻且多是逆時針旋轉；余流有3～5 d的波動周期，主要源于對經(jīng)向風的響應。Pang等[6]利用POM模式得出在青島以東的山東南部近海存在反氣旋式海水運動的結論。張志欣[7]根據(jù)2006—2007年青島近岸水域的海床基測流站數(shù)據(jù)，指出青島外海水域余流是風和山東半島地形共同作用的結果，夏季該水域表層余流呈現(xiàn)西北向，并與南側相鄰海域的余流構成氣旋式環(huán)流，與夏季青島外海冷水團環(huán)流一致。

前人對于黃海海域的水動力環(huán)境做了大量的研究工作并取得了許多重要成果。但是，由于密集的垂釣和拖網(wǎng)捕魚活動使得南黃海長期觀測的海流資料缺乏[8]；其次，以往的調(diào)查多集中在冬、夏兩季，導致研究資料具有較大的時、空局限性。本文旨在總結前人研究成果的基礎上，利用魯海豐海洋牧場站點的實測海流資料，對該海域的流場結構進行深入的分析，進一步掌握該海洋牧場的水動力環(huán)境特征，為其建設提供指導。

1 數(shù)據(jù)和方法

1.1 數(shù)據(jù)

本文采用坐底式觀測平臺在魯海豐海洋牧場站點連續(xù)觀測的流矢量資料，站點位置為120°12′51.69″E, 35°55′41.52″N，站點的平均水深為12.4 m，地形特征如圖1所示。流矢量資料由AWAC浪龍采集所得，采樣時間范圍為2016年3月15日—6月30日，原始采樣頻率為10和20 min。設置垂直測量單元層數(shù)為20層，每層間隔2 m，盲區(qū)為0.51 m，傳感器具有3個波束，自海底向上發(fā)射，根據(jù)層數(shù)返回相應的流速數(shù)據(jù)。考慮到觀測點的水深及儀器在近表層的觀測誤差，本文選擇以海底為基準向上5層的流矢量數(shù)據(jù)進行分析。在進行數(shù)據(jù)分析前，對原始觀測數(shù)據(jù)進行質量控制，剔除奇異值，通過線性插值補齊缺省值并對每層數(shù)據(jù)進行小時平均，得到觀測期間每層2 592個數(shù)據(jù)。

(紅色五角星為魯海豐海洋牧場位置。Red star is the location of LHF Marine Ranch.)

為了解觀測點及鄰近海域的地形狀況，本文采用美國國家地球物理數(shù)據(jù)中心(The National Geophysical Data Center)發(fā)布的全球海底地形數(shù)據(jù)etopo1和etopo5(https://www.ngdc.noaa.gov/mgg/global/global.html)，etopo1數(shù)據(jù)網(wǎng)格精度為1′×1′, etopo5數(shù)據(jù)網(wǎng)格精度為1/12°×1/12°。

為探究觀測點的流矢量與海面風的相關關系，本文采用歐洲中尺度天氣預報中心(ERA-Interim)發(fā)布的10 m處風矢量數(shù)據(jù)(http://apps.ecmwf.int/datasets/data/interim-full-daily/levtype=sfc/), 選取時間范圍為2016年3月15日—6月30日，水平分辨率為0.125°×0.125°，采樣間隔為6 h。將ERA-Interim再分析的海面10 m處風場資料進行雙線性插值獲得觀測點處的風矢量數(shù)據(jù)，并進行日平均處理。

為探究鄰近海域流場對觀測點流場的作用，本文采用HYCOM(Hybrid Coordinate Ocean Model)再分析的日平均表層流矢量資料(https://hycom.org/dataserver/gofs-3pt0/analysis), 選取范圍為20°N～42°N，115°E～135°E，選取的時間范圍為2016年3月15日—6月30日，數(shù)據(jù)水平分辨率為1/12°×1/12°。

1.2 方法

本文通過功率譜分析、調(diào)和分析、旋轉譜分析、相關分析等方法對觀測點的平均流、潮流和余流特征進行分析和研究。

利用T_TIDE MATLAB程序包對實測海流資料進行潮流調(diào)和分析[9]，得到各分潮的頻率、長軸、短軸、傾角、初相位等要素信息。根據(jù)中國《港口工程技術規(guī)范》(中華人民共和國交通部，1988)的規(guī)定，潮流類型特征值可以按照下式計算得到：

(1)

其中UK1、UO1、UM2分別為K1,O1,M2的最大分潮流，即潮流橢圓的長半軸。

2 結果

2.1 平均流

鑒于當?shù)仄骄顬?2.4 m, AWAC自海底向上獲得流矢量資料，故分別取第1、3、5層作為觀測點的底、中、表層。U、V別表示緯向流和經(jīng)向流，在U、V方向上分別計算各層流的相關系數(shù)(見表1)，結果均大于0.8，說明上下層水體的運動具有良好的一致性。

表1 各層流速相關系數(shù)

2.1.1 垂向平均流 對1～5層的實測流資料進行垂向平均，然后分別對U分量和V分量做功率譜分析，得到垂向平均流的功率譜(見圖2)。觀測點的經(jīng)向流和緯向流均表現(xiàn)出12.35 h的顯著周期，并通過了紅噪聲置信度檢驗，說明觀測點的流具有顯著的半日周期波動特征。另外，U、V分量流速還呈現(xiàn)出幾個更高頻的周期，但是這種波動周期是否顯著，還需進一步驗證。

圖2 垂向平均流功率譜分析

2.1.2 觀測期間平均流 對實測海流資料進行48 h低通濾波去除潮流的影響，再對得到的余流場進行時間平均并計算標準差，得到整個觀測期間的平均流和標準差(見圖3)。結果顯示，觀測期間U、V方向的平均流速均在10 cm/s以內(nèi)，U分量表現(xiàn)為西向流，V分量表現(xiàn)為北向流。由表及底，西向流先增大后減小，至底層達到最小；北向流隨深度緩慢增大，至底層趨于穩(wěn)定且達到最大。U分量的標準差隨深度減小，減小的速度先快后慢，說明西向流在表層變化較大，至中層和底層趨于穩(wěn)定；V分量的標準差在底層以上與U分量的標準差變化一致，但在底層突然增大?？傮w上，由表及底U分量減弱，V分量增強，且流速逐漸趨于平穩(wěn)。

2.1.3 月平均流 對實測海流資料進行48 h低通濾波去除潮流的影響，再對得到的余流場進行月平均處理并分別計算標準差，得到每月的平均流和標準差的垂向分布圖(見圖4)。整體看來，月平均流在垂向上的分布情況與總的平均流相差不大。其中4月的U分量最強，V分量最弱，且V分量在表層出現(xiàn)流速轉向；6月的U分量最弱，V分量最強，U分量隨深度逐漸減小，V分量隨深度先減小后增大而后減小，變化較大；3和5月呈現(xiàn)中等強度的U向流和V向流。每月的標準差與上述總的標準差變化基本一致，U分量的標準差在3月最小，在5月最大；V分量的標準差在6月最小，在其它月份量值相差無幾，另外，在5月V分量的標準差至底層突然增大，表現(xiàn)出流動的不穩(wěn)定性。

2.2 調(diào)和分析和旋轉譜分析

2.2.1 潮流調(diào)和分析 對實測海流資料進行潮流調(diào)和分析，得到O1、K1、M2、S2四個最顯著分潮的潮流橢圓參數(shù)(見表2)?？傮w看來，觀測點處M2分潮具有最大的分潮流速，其次為S2分潮，O1和K1分潮相對較小，表現(xiàn)出半日潮族和全日潮族的特征。潮流類型判別標準表明，A≤0.5為正規(guī)半日潮流，計算得表、中、底層的A分別為0.27, 0.21, 0.21，均小于0.5，故觀測點的潮流類型為正規(guī)半日潮。

圖3 平均流速和標準差垂向分布圖

(實線表示U分量；虛線表示V分量。The solid line represents the U component and the dotted line represents the V component.)

表2 潮流橢圓參數(shù)

短軸的正負表示潮流橢圓的旋轉方向，正值表示潮流橢圓左向旋轉(逆時針)，負值表示潮流橢圓右向旋轉(順時針)[10]。畫出各層的潮流橢圓(見圖5)，可以看出，當?shù)爻绷鞅憩F(xiàn)為以東北向和西南向流為主導的往復流，與當?shù)貣|北-西南向的岸線一致。其中，M2(分潮)為第一大主要分潮，在各層占有絕對優(yōu)勢，由表到底幾乎不發(fā)生偏轉，只是強度略有變化；S2分潮為第二大分潮，由表及底傾角變化不大，強度相對來說有明顯減??；O1和K1分潮在各層都很小，由表及底強度明顯變小且旋轉明顯，在表層時均為右旋的橢圓，而到中、底層變?yōu)樽笮?，這與Tang等[5]的結論基本一致。

圖5 表、中、底層潮流橢圓

2.2.2 流矢量的旋轉譜分析 流速數(shù)據(jù)是在地球慣性坐標系中測量的，流速矢量具有旋轉性，隨時間做順時針或逆時針方向旋轉，單純地把流速矢量分解成東、北兩個正交分量并分別進行功率譜分析和調(diào)和分析的做法是存在缺陷的[11]。為了更好地反映驅動力的作用和地球自轉對海水的影響，本文對流速矢量進行旋轉譜分析，結果如圖6所示。

(實線代表逆時針旋轉譜，虛線表示順時針旋轉譜。The solid line represents the counterclockwise rotary spectra and the dotted line represents the clockwise rotary spectra.)

通過圖6可以看出觀測點處半日周期潮流占主，全日潮流次之。半日周期(12.4 h)潮流的逆時針部分在整個水體厚度上均略大于順時針的部分，說明當?shù)氐陌肴粘绷鞅憩F(xiàn)為一種逆時針旋轉的近乎直線的往復運動。全日周期潮流的順時針部分在表層略大于逆時針部分，而在中層和底層相反，說明表層全日潮流(24 h)為順時針旋轉，隨著深度增加改變旋轉方向，且全日潮流也表現(xiàn)出往復流的特點，這與調(diào)和分析的結果完全一致。

另外，觀測點處的慣性流為順時針方向，由慣性運動的周期公式T=2π/f，其中f為觀測點所在緯度的科氏參數(shù)，計算知當?shù)貞T性流周期約為20.5 h，與全日潮流周期接近。分段選取流矢量資料做旋轉譜分析，得到各層順時針旋轉的全日潮流的周期均為24 h左右，即未表現(xiàn)出顯著的慣性周期特征，說明總體上觀測點表現(xiàn)出的日周期波動只由全日潮流引起，不受慣性流影響。

2.3 余流分析

2.3.1 觀測點余流 余流是指從實測海流中剔除周期性流動之后剩余的水體流動[12]。為了探究當?shù)赜嗔鞯慕Y構特征及其對海面風的響應，從原始的觀測海流中減去經(jīng)調(diào)和分析得到的周期性潮流成分，然后再進行24 h平均即得到觀測點處的日平均的余流時間序列，海面10 m處風速和觀測點處各層余流的時間序列如下圖7所示。

通過余流的時間序列可以看出觀測點的余流主要呈現(xiàn)西北向，且由表及底余流減小。其中，表層余流最強，方向變化劇烈，在觀測期間表現(xiàn)出順時針旋轉的特征，中層和底層余流較弱且流向平穩(wěn)，以西北向和北向流為主。結合風速的時間序列可以看出，在觀測期間當?shù)匾阅巷L為主，表層余流對風的響應敏感，余流隨時間的變化與風基本一致，風較強時余流加強，風較弱時余流較弱，風轉向時余流也隨之轉向，計算風場與表層余流的相關系數(shù)，得到緯向風與緯向流的相關系數(shù)為0.68，經(jīng)向風與經(jīng)向流的相關系數(shù)為0.64(置信水平超過95%)，而對比中層和底層余流與海面風的時間序列發(fā)現(xiàn)中層和底層余流對海面風的響應差。因此魯海豐海洋牧場的表層余流主要受風場驅動，而中層和表層余流受風影響較小。

圖7 海面風矢量與余流矢量圖

2.3.2 觀測點鄰近海域余流 為了進一步探究觀測點的余流結構，結合觀測期間平均的HYCOM表層流場資料(HYCOM再分析的流資料已剔除潮流)，得到觀測點與鄰近海域的平均余流場對比圖(見圖8)。基于HYCOM資料得到的余流場，發(fā)現(xiàn)在青島近岸海域呈現(xiàn)東北向流，與青島近岸西南-東北向的地形基本一致，與Tang等[5]對站點C的觀測結果一致。同時，在觀測點的東南部存在一個反氣旋式環(huán)流，該環(huán)流的中心大約位于121°30′E, 35°30′N，與青島冷水團位置相近，該結果與Pang等[6]用POM模式得出的在青島以東的山東南部近海存在反氣旋式的海水運動的結論相符，證明了結合HYCOM再分析資料來探究觀測點的流場特征是可行的。

觀測期間平均海面10 m處風場(見圖8)表明，在整個研究區(qū)域盛行東南風和南風，其中觀測點附近為較強的東南風。觀測點離岸較近，其平均流在各層均表現(xiàn)為西北向，呈現(xiàn)出向岸流的特點，且由表及底，西向流減弱，北向流增強，方向右旋。為探究實測流表現(xiàn)出的右旋特征，利用有限深海Ekman理論的解析解進行驗證計算，發(fā)現(xiàn)計算得到的Ekman流方向為東北向，由表及底右旋，這與實測流方向并不一致，且量值遠小于實測流的量值。因此觀測點處的平均流表現(xiàn)出的右旋特征受到多種因素的影響，不能簡單描述該流動為有限深海Ekman流。

(黑色箭頭表示HYCOM表層流場;綠色箭頭表示鄰近海域10 m處平均風場;紅色五角星表示魯海豐海洋牧場的位置;藍色箭頭表示表層平均流;紅色表示中層流;黃色表示底層流。The black arrows represent the surface current field of HYCOM data；The green arrows represent the mean wind field at the adjacent sea area；The red star represents the location of LHFMarine Ranch；The blue arrow represents the surface mean flow；The red arrow represents the middle mean flow；The yellow arrow represents the bottom mean flow.)

將觀測的表層海流做天平均，然后與鄰近海域相同時間長度的HYCOM表層流場資料做相關性分析，得到通過95%置信水平的觀測點的日平均流與鄰近海域日平均流的相關系數(shù)等值線圖(見圖9)。在觀測期間內(nèi)，觀測點表層平均流的U分量與渤海萊州灣，黃海西部近岸海域和黃海中部海域相關性較高，表層平均流的V分量與渤海的渤海灣、萊州灣以及山東半島南部近岸海域相關性較好，表明觀測點海流具有的低頻時間變化特征與上述海域較一致,因此觀測點的余流也受到山東半島沿岸流的影響。

(a.U分量; b.V分量。a. U component; b.V component.)

3 結論

通過對魯海豐海洋牧場觀測的3個半月的流矢量資料進行分析并結合鄰近海域的海流資料進行研究，得出如下結論：

(1) 觀測點各層的水體運動具有良好的一致性，時間平均的垂向流顯示，在整個水體厚度上，西北向流占優(yōu)，呈現(xiàn)出向岸流的特點，且西向流隨深度減小，北向流隨深度增大。

(2) 月平均流在垂向上的分布情況與總的平均流相差不大，其中4月份的U分量最強，V分量最弱，且V分量在表層出現(xiàn)流速轉向；6月份的U分量最弱，V分量最強，且U分量隨深度逐漸減小，V分量隨深度先減小后增大而后減小，變化較大；3和5月份呈現(xiàn)中等強度的U向流和V向流。

(3) 功率譜分析表明，觀測點的流具有顯著的半日周期特征；潮流調(diào)和分析表明，觀測點的潮流類型為正規(guī)半日潮，M2為主要分潮流，在各層占有絕對優(yōu)勢，各層潮流均顯示出旋轉流性質，旋轉方向以逆時針為主；旋轉譜的分析結果與調(diào)和分析結果一致，總體上觀測點的日周期波動只由全日潮流引起，不受當?shù)貞T性流影響。

(4) 觀測點表層余流對風的響應敏感，緯向風與表層緯向余流的相關系數(shù)為0.68，經(jīng)向風與表層經(jīng)向余流的相關系數(shù)為0.64；結合HYCOM再分析的表層流場資料發(fā)現(xiàn)觀測點的日平均流與鄰近的渤黃海近岸海域變化較一致，說明觀測點的平均余流受到海面風的驅動和山東半島沿岸流的影響。