姜?jiǎng)佥x，朱龍海*，胡日軍，張 偉，劉愛江

（1.中國(guó)海洋大學(xué) 海底科學(xué)與探測(cè)技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，山東 青島 266100；2.中國(guó)海洋大學(xué)海洋地球科學(xué)學(xué)院，山東 青島 266100；3.青島海洋地質(zhì)工程勘察院，山東 青島 266071）

萊州灣為渤海三大海灣之一，為一半封閉性海灣，由于其位置的獨(dú)特性，生物資源豐富，經(jīng)濟(jì)發(fā)達(dá)，生態(tài)環(huán)境相對(duì)較脆弱。近年來，海岸資源開發(fā)利用的強(qiáng)度、密度大幅增加，產(chǎn)生了自然岸線資源縮減、海灣屬性弱化、海灣污染加劇等一系列問題［1］。

隨著萊州灣周邊的開發(fā)建設(shè)，科研工作者在這一區(qū)域進(jìn)行了大量的研究。利用ECOM模型研究了黃河口改道對(duì)萊州灣潮流場(chǎng)造成的影響［2］；利用衛(wèi)星遙感圖像和觀測(cè)資料，得出底質(zhì)中細(xì)顆粒物質(zhì)含量較高的海域容易在潮流作用下再懸浮，有助于了解懸浮物在渤海中的運(yùn)輸狀態(tài)［3］；最近50年來，萊州灣西、南海岸在地貌演變上存在顯著的差異性，這種差異主要與自然環(huán)境條件和海洋沉積動(dòng)力過程有關(guān)［4］。前人對(duì)萊州灣的水動(dòng)力條件、泥沙輸運(yùn)有了大量的研究，但萊州灣水動(dòng)力變化的研究相對(duì)匱乏。本文在前人水動(dòng)力研究的基礎(chǔ)上通過對(duì)比萊州灣2003和2013年兩個(gè)時(shí)間段的潮流場(chǎng)、納潮量和水交換變化，分析圍填海工程建設(shè)對(duì)水動(dòng)力的影響。

1 研究區(qū)概況

萊州灣西起現(xiàn)代黃河口，東迄屺姆島高角，本文研究區(qū)域?yàn)辄S河口至刁龍嘴所圍海域（見圖1），為萊州灣的內(nèi)灣。萊州灣內(nèi)灣水深較淺，最大水深15m，大部分區(qū)域水深在10m以內(nèi)；萊州灣的潮汐主要受黃河口外半日無潮點(diǎn)的影響，潮汐類型屬不規(guī)則的混合半日潮海區(qū)；受海灣形狀影響，萊州灣的強(qiáng)浪向和常浪向以北向浪（NE-NW）為主［5］。2003年內(nèi)灣的海域面積約4 376.24km2，近年來由于大規(guī)模填海造地、圍海養(yǎng)殖、鹽田圍海、港口開發(fā)等開發(fā)活動(dòng)，海灣面積有所減小。

圖1 地理位置圖Fig.1 Geographical position

2 資料與方法

本文運(yùn)用MIKE21模擬軟件，采用數(shù)值模擬的方法，分別模擬了2003和2013年萊州灣的潮流場(chǎng)、納潮量、水交換情況，通過對(duì)比分析，研究近十年來圍填海工程建設(shè)對(duì)萊州灣水動(dòng)力影響。

2.1 水深地形和岸線

水深數(shù)據(jù)主要采用中國(guó)航海圖書出版社出版發(fā)行的萊州灣1∶15萬(wàn)海圖，同時(shí)結(jié)合近岸港口建設(shè)實(shí)測(cè)的大比例尺水深地形測(cè)量數(shù)據(jù)。

研究區(qū)海岸線采用平均大潮高潮線，本文收集了萊州灣2003年5月Landsat-7ETM＋影像數(shù)據(jù)和2013年5月landsat8影像數(shù)據(jù)，其空間分辨率均可達(dá)到15m。其中，自然岸線區(qū)從海圖中提取岸線，人工岸線區(qū)通過兩期衛(wèi)星圖像配準(zhǔn)后解譯。由于人工海岸是由水泥和石塊構(gòu)筑，在圖像上具有較高的光譜反射率，與光譜反射率很低的海域區(qū)分明顯，能夠準(zhǔn)確的確定其海岸線［6－7］。

2.2 潮位和潮流

潮位資料收集了濰河口和彌河口2個(gè)潮汐站歷史觀測(cè)資料經(jīng)調(diào)和分析后，選用 M2、S2、K1、O14個(gè)分潮的調(diào)和常數(shù)預(yù)報(bào)出大潮期的潮位與計(jì)算結(jié)果進(jìn)行驗(yàn)證；潮流資料利用長(zhǎng)江委水文局長(zhǎng)江下游水文水資源勘測(cè)局2007年3月6—8日（大潮期）進(jìn)行的2個(gè)站位（C3、C4）和交通運(yùn)輸部天津水運(yùn)工程科學(xué)研究所2011年3月12—13日進(jìn)行的2個(gè)站位（C1、C2）海流觀測(cè)資料進(jìn)行潮流驗(yàn)證，驗(yàn)證點(diǎn)位置如圖2所示。

圖2 網(wǎng)格分布及驗(yàn)證點(diǎn)位置圖Fig.2 Grid distribution and verification point

2.3 數(shù)值模型

數(shù)值模擬采用MIKE21模型，該模型適用于模擬河口、海灣、海岸及海洋的水流、波浪、泥沙及環(huán)境，國(guó)內(nèi)外眾多學(xué)者采用該模型進(jìn)行海灣研究。模型采用非結(jié)構(gòu)網(wǎng)格剖分計(jì)算域，標(biāo)準(zhǔn)Galerkin有限元法進(jìn)行水平空間離散，在時(shí)間上采用顯式迎風(fēng)差分格式離散動(dòng)量方程和輸運(yùn)方程［8］。論文采用模型中水動(dòng)力模塊和對(duì)流擴(kuò)散模塊，進(jìn)行潮流場(chǎng)、水交換和納潮量模擬。模型采用潮汐驅(qū)動(dòng)的方式，開邊界潮位由北黃海及渤海大模型提取，采用三角網(wǎng)格，計(jì)算域?yàn)辄S河口和萊州港連線及岸線所圍海域，最小空間步長(zhǎng)約10m，由于萊州灣水深較淺，垂向上不分層。模擬了K1、O1、M2、S2 4個(gè)主要分潮驅(qū)動(dòng)下潮流場(chǎng)狀況，經(jīng)潮位、潮流驗(yàn)證，模擬結(jié)果與實(shí)測(cè)結(jié)果吻合度較好。計(jì)算域及網(wǎng)格劃分見圖2。

2.4 計(jì)算方法

2.4.1 納潮量計(jì)算方法 納潮量是一個(gè)水域可以接納的潮水的體積，是表征半封閉海灣生存能力的重要指標(biāo)。納潮量的計(jì)算采用傳統(tǒng)的算法，假定潮灘為坡度均勻的斜面，進(jìn)行大潮期納潮量的推算［9］。

其中：H2－H1為潮差；S1為高潮時(shí)海灣的面積；S2為低潮時(shí)海灣的面積。

2.4.2 水交換計(jì)算方法 假定在口門附近海域有同樣的溶解態(tài)保守污染物，使水體中污染物質(zhì)的濃度達(dá)到某一個(gè)特定的量值Wc，且初始時(shí)刻均勻分布在附近水域內(nèi)，假設(shè)水邊界入流時(shí)給定這種物質(zhì)在開邊界的濃度為0，計(jì)算出不同時(shí)刻此種物質(zhì)的濃度值W，進(jìn)而通過以下公式計(jì)算不同時(shí)刻被區(qū)外海水置換的比率（也即水體交換率），來比較規(guī)劃方案實(shí)施后附近海域的水體凈化能力［10－11］。

式中：n為水體交換率，WC為原有的污染物濃度值，W為每一時(shí)刻水體中的污染物濃度值。

3 萊州灣變化

通過對(duì)2003和2013年萊州灣遙感影像對(duì)比發(fā)現(xiàn)，2003—2013年間萊州灣西側(cè)東營(yíng)市海域圍填海工程較少，中部濰坊市海域開展了較大規(guī)模的圍填海，特別是濰坊港垂直岸線垂直岸邊防波堤向?yàn)硟?nèi)延伸長(zhǎng)度達(dá)17.3km，大幅改變了萊州灣的水動(dòng)力環(huán)境；東側(cè)萊州市海域圍填海工程主要位于岸線附近，對(duì)水動(dòng)力環(huán)境影響較小。至2013年萊州灣內(nèi)灣的面積減小至4 053.13km2，面積減小了323.11km2，萊州灣減小的面積占整個(gè)海灣面積的7.38%。

4 結(jié)果

4.1 模型驗(yàn)證

4.1.1 潮位驗(yàn)證 利用濰河口、彌河口的潮位觀測(cè)資料與計(jì)算結(jié)果進(jìn)行驗(yàn)證。潮位驗(yàn)證曲線如圖3a、b所示。潮位驗(yàn)證結(jié)果表明，對(duì)應(yīng)觀測(cè)點(diǎn)上模擬得到的潮位值與實(shí)測(cè)潮位基本吻合。

4.1.2 潮流驗(yàn)證 利用潮流模型模擬萊州灣2007年3月6—8日的潮流場(chǎng)，提取C3、C4站位的流速流向與實(shí)測(cè)潮流進(jìn)行對(duì)比，潮流驗(yàn)證曲線如圖4c～d所示。然后，在不改變曼寧系數(shù)等參數(shù)的情況下，利用潮流模型分別模擬萊州灣2011年3月12—13日和2014年5月28—29日的潮流場(chǎng)，提取C1、C2站位的流速流向與2011年實(shí)測(cè)潮流進(jìn)行對(duì)比，潮流驗(yàn)證曲線如圖4a～b所；提取C5、C6站位的流速流向與2014年實(shí)測(cè)潮流進(jìn)行對(duì)比，潮流驗(yàn)證曲線如圖4e～f所示。

潮流驗(yàn)證結(jié)果表明，對(duì)應(yīng)觀測(cè)點(diǎn)上模擬得到的潮流流速流向與實(shí)測(cè)潮流基本吻合，能夠較好地反映規(guī)劃周邊海域潮流狀況。

4.2 數(shù)值模擬結(jié)果

本文分別對(duì)萊州灣2003和2013年潮流場(chǎng)、納潮量和水交換情況進(jìn)行數(shù)值模擬。

4.2.1 潮流場(chǎng)模擬結(jié)果 本文中給出研究區(qū)潮流場(chǎng)潮時(shí)變化是以濰坊港中港區(qū)口門處的潮汐變化為參照時(shí)間。

2003 年潮流場(chǎng)的模擬結(jié)果：圖5a是漲急時(shí)潮流場(chǎng)，萊州灣內(nèi)的潮流由灣外向?yàn)硟?nèi)流動(dòng)，潮流流速一般在30～50cm/s之間。北部海域潮流流向?yàn)槟衔飨?，萊州淺灘處流速達(dá)到最大值，為106cm/s。繞過萊州淺灘，潮流轉(zhuǎn)向南流動(dòng)，進(jìn)入太平灣后，向偏東方向流動(dòng)。萊州淺灘以西至黃河口之間海域，潮流由東北往西南流，至灣頂附近逐漸轉(zhuǎn)為近西向流，在黃河口南側(cè)近岸海域潮流由東南往西北流。圖5b為落急時(shí)潮流場(chǎng)，潮流分布情況與漲急時(shí)潮流場(chǎng)相似，只是流向與漲急時(shí)潮流場(chǎng)大致相反。潮流整體由萊州灣內(nèi)向?yàn)惩饬鲃?dòng)，潮流流速一般在30～50cm/s之間，在萊州淺灘處流速最大值可達(dá)119cm/s。

圖3 潮位驗(yàn)證曲線Fig.3 Tidal level validation curve

圖4 潮流驗(yàn)證曲線Fig.4 Tidal current validation curve

圖5 a 漲急時(shí)潮流場(chǎng)Fig.5a Tidal current field during flood tide

圖5 b 落急時(shí)潮流場(chǎng)Fig.5b Tidal current field during ebb tide

2013 年潮流場(chǎng)的模擬結(jié)果：2013年潮流場(chǎng)與2003年潮流趨勢(shì)基本一致，漲急時(shí)萊州灣內(nèi)的潮流由灣外向?yàn)硟?nèi)流動(dòng)，落急時(shí)潮流整體由萊州灣內(nèi)向?yàn)惩饬鲃?dòng)，海灣外部流速潮流場(chǎng)變化不大，靠近岸邊區(qū)域整體呈不同程度的減小。

4.2.2 納潮量計(jì)算結(jié)果 本文分別模擬了2003和2013年萊州灣大潮期和小潮期的納潮量，納潮量模擬結(jié)果見表1，2003年萊州灣大潮期納潮量為63.94×108m3，小潮期納潮量為45.44×108m3，納潮量平均為54.69×108m3；2013年萊州灣大潮期納潮量為60.08×108m3，小潮期納潮量為43.22×108m3，納潮量平均為51.65×108m3；大潮期納潮量減小了約6.05%，小潮期減小了4.90%，平均減小了5.57%。由對(duì)比可知，十年來萊州灣的納潮量發(fā)生變化，呈現(xiàn)減小的趨勢(shì)。

表1 2003—2013年萊州灣納潮量Table 1 Storage capacity of tidal water in Laizhou Bay from 2003to 2013

4.2.3 萊州灣水體交換預(yù)測(cè)結(jié)果 本文分別模擬了2003和2013年萊州灣15和30d的水交換情況，由模擬結(jié)果可知：水交換15d，水交換率大于10%的面積均大于63%，水交換率大于30%的面積均大于40%，水交換大于50%的面積占整個(gè)海灣面積25%左右；水交換30天，萊州灣內(nèi)所有海域水交換率基本上均大于10%，水交換率大于30%的面積均大于58%，水交換大于50%的面積占整個(gè)海灣面積30%左右。2003和2013年萊州灣水交換率情況見表2。

5 討論

5.1 圍填海工程對(duì)潮流場(chǎng)的影響

為了研究2003—2013年間圍填海工程建設(shè)對(duì)萊州灣潮流影響，選取大潮時(shí)漲急時(shí)和落急時(shí)潮流場(chǎng)進(jìn)行流速對(duì)比，漲急時(shí)流速變化見圖6a、落急時(shí)流速變化見圖6b。

表2 2003年和2013年萊州灣水交換率情況Table 2 The rate of water exchange in Laizhou Bay form 2003to 2013

由對(duì)比結(jié)果可知，近十年圍填海工程的建設(shè)對(duì)萊州灣潮流場(chǎng)產(chǎn)生了一定的影響，其中對(duì)近岸區(qū)域和萊州灣中部影響較大，主要是由于近岸工程的建設(shè)改變了近岸區(qū)的水動(dòng)力環(huán)境，濰坊港防波堤向萊州灣內(nèi)延伸較長(zhǎng)，對(duì)萊州灣水動(dòng)力影響較大。圍填海工程使萊州灣流速整體減小，減小較大的區(qū)域集中于濰坊港東北側(cè)，流速減小超過4cm/s（變化率約10%）的區(qū)域漲急和落急時(shí)平均面積達(dá)760km2，約占整個(gè)海灣面積的24.3%；流速增大區(qū)的面積和增加量均較小，流速增大大于4cm/s區(qū)域僅占海灣面積的4.0%。根據(jù)流向?qū)Ρ冉Y(jié)果，流向變化主要位于靠近岸邊區(qū)域，灣內(nèi)流向變化一般小于5°，流向變化整體較小。

5.2 圍填海工程對(duì)納潮量的影響

對(duì)比2003年和2013年納潮量結(jié)果可知，2003—2013年萊州灣大量圍填海工程的建設(shè)導(dǎo)致了萊州灣納潮量平均減小了5.57%。分析萊州灣的納潮量減小的原因主要有以下2個(gè)方面：圍填海工程的建設(shè)直接占用了海域面積，影響海灣納潮的空間；圍填海工程建設(shè)使萊州灣內(nèi)的水動(dòng)力減弱，潮差有所減小，根據(jù)數(shù)值模擬結(jié)果，大潮期萊州灣潮差減小了約15.3mm，小潮期潮差減小了5.4mm，潮差的減小導(dǎo)致了納潮量的減小。

5.3 圍填海工程對(duì)海灣水交換率的影響

統(tǒng)計(jì)結(jié)果表明，2013年較2003年各交換率各時(shí)間段的水交換均有不同程度的減小；交換15d相同交換率面積減小量在0.87%～1.30%之間，30d的相同交換率面積減小量在0.55%～7.02%之間，2013年水交換面積大于10%的減小量為7.02%，減小較大是由于圍填海區(qū)直接占用了2003年水交換率大于10%的區(qū)域。2003—2013年圍填海工程一定程度上減小了萊州灣的水交換率，萊州灣的自凈能力有所減弱。由于圍填海工程主要沿岸邊建設(shè)，圍填海工程大部分位于灣頂、沒有改變口門的寬度，圍填海工程建設(shè)對(duì)萊州灣水交換影響相對(duì)較小。

5.4 對(duì)比分析水動(dòng)力變化

2000—2010 年圍填海工程導(dǎo)致錦州灣海域面積減小了31.0%，特征點(diǎn)流速變化量在7.99%～15.13%之間，納潮量減少了11.92%，水交換率減少了17.51%［12］。與圍填海工程建設(shè)較多的錦州灣相比，萊州灣的水動(dòng)力變化不明顯；但十年來海灣面積減小7.38%，萊州灣水動(dòng)力的變化遠(yuǎn)超自然演變速率，圍填海工程對(duì)水動(dòng)力的影響不容忽視。近年來，廈門灣等開展了綜合整治工程，水動(dòng)力條件有了較大改善［13］。相關(guān)部門應(yīng)科學(xué)編制萊州灣圍填海規(guī)劃，實(shí)現(xiàn)萊州灣海域的可持續(xù)發(fā)展［14］。

圖6 a 2003—2013年漲急時(shí)潮流流速變化Fig，6a Tidal current velocity variation during flood tide from 2003to 2013

圖6 b 2003—2013年落急時(shí)潮流流速變化圖Fig.6b Tidal current velocity variation during ebb tide from 2003to 2013

6 結(jié)論

（1）2003—2013年萊州灣大面積的圍填海建設(shè)導(dǎo)致萊州灣海域面積減小了7.38%。

（2）圍填海工程使萊州灣流速整體減小，流速減小超過4cm/s（變化率約10%）約占整個(gè)海灣面積的24.3%，集中在濰坊港周邊；流速增加大于4cm/s區(qū)域僅占海灣面積的4.0%，灣內(nèi)流向變化一般小于5°，流向變化整體較小。

（3）圍填海工程的建設(shè)使萊州灣的納潮量減小，大小潮期平均減小了5.57%。

（4）由于圍填海工程主要沿岸邊建設(shè)且沒有改變口門的寬度，圍填海工程建設(shè)對(duì)萊州灣水交換影響相對(duì)較小。與2003年相比，2013年15天水體交換率減小0.87%～1.30%，30天水體交換率減小0.55%～7.02%。

（5）2003—2013年圍填海工程的建設(shè)使萊州灣潮流場(chǎng)、納潮量、水交換不同程度的減小，海灣屬性弱化，將會(huì)對(duì)萊州灣的資源和環(huán)境帶來一系列不良影響。

［1］吳云凱.萊州灣海洋環(huán)境變化趨勢(shì)及管理措施研究［J］.海洋開發(fā)與管理，2011，9：90-92.

［2］張永強(qiáng)，遲萬(wàn)清，胡澤建，等.黃河清水溝流路大嘴的形成對(duì)萊州灣潮流場(chǎng)影響的數(shù)值研究［J］.海洋科學(xué)進(jìn)展，2010，28（2）：149-156.

［3］江文勝，王厚杰.萊州灣懸浮泥沙分布形態(tài)及其與底質(zhì)分布的關(guān)系［J］.海洋與湖沼，2005，36（2）：97-103.

［4］李蒙蒙，王慶，張安定，等.最近50年來萊州灣西—南部淤泥質(zhì)海岸地貌演變研究［J］.海洋通報(bào)，2013，32（2）：141-151.

［5］中國(guó)海灣志編篆委員會(huì).中國(guó)海灣志（第三分冊(cè)）［M］.北京：海洋出版社，1991.

［6］馬小峰，趙冬至，張豐收，等.海岸線衛(wèi)星遙感提取方法研究進(jìn)展［J］.遙感技術(shù)與應(yīng)用，2007，24（4）：575-579.

［7］Palmer P J，Tilton J C.Refining image segmentation by integration of edge and region data［J］.IEEE Trans Geo sci Remote Sens，1995，33：605-615.

［8］DHI Water ＆Environment.MIKE21User Guide［M］.Denmark：DHI Water ＆Environment，2002.

［9］伍成成，鄭西來，林國(guó)慶.雙臺(tái)子河口趕潮段納潮量研究［J］.水資源與水工程學(xué)報(bào)，2010，21（4）：105-110.

［10］羅鋒，廖光洪，楊成浩，等.樂清灣水交換特征研究［J］.海洋學(xué)研究，2011，29（2）：79-88.

［11］Wang Cui，Zhang Xue-qing，Sun Ying-lan.Numerical Simulation of Water Exchange Characteristics of the Jiaozhou Bay Based on A Three-Dimensional Lagrandian Model［J］.China Ocean Engineering，2009，23（2）：277-290.

［12］劉明，席小慧，等.錦州灣圍填海工程對(duì)海灣水交換能力的影響［J］.大連海洋大學(xué)學(xué)報(bào)，2013，2（28）：110-114.

［13］陸榮華.圍填海工程對(duì)廈門灣水動(dòng)力環(huán)境的累積影響研究［D］.廈門：國(guó)家海洋局第三海洋研究所，2010.

［14］狄乾斌，韓增林.大連市圍填?；顒?dòng)的影響及對(duì)策研究［J］.海洋開發(fā)與管理，2008，10：122-126.