胥　巖，潘　毅，陳永平

（河海大學港口海岸與近海工程學院，江蘇　南京　210098）

綜合

老虎石公園海灘養(yǎng)護岸線演變模擬及預測

胥巖，潘毅*，陳永平

（河海大學港口海岸與近海工程學院，江蘇南京210098）

在實測岸線資料的基礎上對老虎石公園海灘養(yǎng)護工程的養(yǎng)灘效果進行了分析，然后基于GENESIS數(shù)值模型構建了岸線演變模型，采用經(jīng)驗證的模型對工程后3~5 a內(nèi)的岸線演變趨勢進行了預測和分析。結果表明，工程結束后，東段海灘岸線后退，5 a后逐漸達到平衡；中段海灘逐漸淤長，5 a后海灘寬度增長40 m左右；西段海灘持續(xù)遭受侵蝕。

GENESIS模型；人工養(yǎng)灘；岸線演變

1　概述

20世紀中葉以前，人們通常通過硬性海岸防護建筑（潛堤、防波堤和海墻等）來抵御海岸的風浪侵蝕。多年的工程應用和研究結果表明，采用硬性海岸防護建筑會暫時抵御風浪的侵蝕，但會改變當?shù)厮畡恿l件，引起泥沙橫向運動失衡，離岸輸運增強。長遠看來，反而加重了海灘的侵蝕程度。為了克服海岸硬防護的不足，人們開始傾向于海岸軟防護。一般說來，若當?shù)赜谐浞值难a沙來源，人工養(yǎng)灘（Beach nourishment）是最佳且應優(yōu)先考慮的海岸軟防護方式[1]。

1922年美國紐約的柯尼島公共岸灘工程是現(xiàn)存記錄中最早的人工養(yǎng)灘工程，截至1996年，美國的養(yǎng)灘地點已超390處，養(yǎng)灘工程已超過1 300 個[2]。美國大規(guī)模的人工養(yǎng)灘工程在保護海岸免受侵蝕、降低風暴潮災害損失的同時，還改善了海岸環(huán)境，促進了旅游業(yè)的發(fā)展[2]。歐洲人工養(yǎng)灘工程的規(guī)模與美國接近[3]，工程大都采用了補沙與丁壩和離岸堤相結合的措施，并定期地對養(yǎng)灘工程進行后續(xù)補沙[4]。我國大陸的人工養(yǎng)灘工作最早開展于20世紀70年代，如青島、茂名等地，但規(guī)模較小，缺乏保護措施，每年的補沙都侵蝕的很快[5]。1990年香港淺水灣工程是我國最早的較大規(guī)模的人工養(yǎng)灘工程[2]。

老虎石公園位于中國著名的避暑和療養(yǎng)勝地——北戴河風景區(qū)的中心。由于諸多原因，老虎石海灘的海洋資源和生態(tài)環(huán)境均受到了一定程度的損害。為了遏制海灘坡度變陡、砂質(zhì)粗化等問題，有效改善、恢復海岸帶旅游功能，于2013 年3—4月實施了老虎石海灘養(yǎng)護工程，將海岸線向海推進了約30 m[6]。本文采用GENESIS（Generalized Model for Simulating Shoreline change）數(shù)值模擬的方法，在驗證模型結果與實測資料基本吻合的前提下，對老虎石海灘養(yǎng)護工程的養(yǎng)護效果進行了分析、預測和評價，為以后進一步的海灘養(yǎng)護工程提供參考和依據(jù)。

2　工程簡介

老虎石海灘長約820 m，西起原旅游碼頭，東至浪濤嶼浴場，為基巖岬灣海岸和砂質(zhì)海岸的混合海岸。海灘養(yǎng)護工程方案見圖1，海灘分3段進行灘肩補沙，補沙粒徑在0.42～0.62 mm之間，老虎石岬頭前方東西兩側各建1個潛壩。西側沙壩長250 m，壩頂高程為-0.5 m，距離岸線約250 m；東側沙壩長440 m，壩頂高程為-0.8 m，距岸約195 m；壩底寬50 m，壩頂寬30 m[7]。

圖1　工程概況圖（2005年航拍圖）Fig.1　Overview of beach nourishment (Aerial Photo taken in 2005)

2.1地形資料

整個北戴河海域內(nèi)，金山嘴附近水深較深、地形較陡，其它海域地形平緩，呈現(xiàn)自西北向東南緩傾。海底地形主要特征可以概括為：除近岸處坡度較陡外，整體上坡度較緩，0～2 m等深線海域坡度較大，介于1.18%～1.72%，2～5 m等深線坡度為0.32%～0.34%，5～10 m等深線之間坡度最小，不足0.1%，介于0.094%～0.098%之間[8]。

2.2波浪資料

波浪測站位于距老虎石海灘約750 m處，波浪資料為2011年6月—2014年6月測得的波浪數(shù)據(jù)，波浪玫瑰圖如圖2所示。該海域的波浪總體較小，測量期間絕大多數(shù)波浪的波高在0.9 m以下，其中，0.6 m以下的波浪占全部波浪的2/3以上。

圖2　波浪頻率玫瑰圖（%）Fig.2　Rose diagram of wave height(%)

3　實測岸線資料分析

為了掌握工程后的岸灘演變趨勢，為后續(xù)的海灘維護提供參考資料，在工程前后，進行了海灘剖面監(jiān)控。根據(jù)各時期剖面數(shù)據(jù)繪制的岸線位置變遷圖見圖3。

圖3　實測岸線位置變遷圖Fig.3　Change of observed shorelines

從圖3可以看出，工程前岸線與底圖所示岸線非常接近，故可以認為在工程前，岸線保持一種相對穩(wěn)定的狀態(tài)，自2005年—2012年變動不大。工程后18個月內(nèi)，在非常明顯的自東向西的沿岸輸沙的作用下，工程填沙發(fā)生了重分布，其中老虎石東側尤為明顯。根據(jù)海灘侵淤特征，工程區(qū)海灘可分為3段來討論：

1）東段海灘的填沙向海突出，在自東向西的沿岸輸沙的作用下，向中段海灘輸移，繼而岸線后退。

2）中段海灘內(nèi)凹于平直的岸線，并有老虎石陸連島阻隔了該段海灘向西的泥沙輸運，故在自東向西的沿岸輸沙作用下，該段海灘攔截了絕大多數(shù)來自于東段海灘的泥沙，海灘面積迅速增長。

3）西段海灘受到老虎石的掩蔽作用但并未形成較大的掩蔽區(qū)，再加上下游泥沙缺乏有效的阻擋，向北戴河西海灘流失，故西段海灘在工程后岸線位置不斷后退。

4　數(shù)值模擬

本文采用GENESIS數(shù)值模擬方法對老虎石海灘養(yǎng)護工程的效果進行預測和評價。GENESIS[9]是基于一線模型所開發(fā)的模擬砂質(zhì)海岸岸線長期演變的系統(tǒng)，在大量的工程應用中逐步完善并成熟起來，現(xiàn)在國際上已廣泛的應用于預測岸線的長期演變及岸線對海岸建筑物和人工養(yǎng)灘等的響應。GENESIS模型的控制方程為：

式中：x為沿岸線方向；y為垂直岸線方向；DB為海水所能到達的最大高程（灘肩高度）；DC為存在沿岸輸沙的最大深度（封閉水深）；q=qs+q0為橫向輸沙率（離岸輸沙+向岸輸沙）；Q為沿岸輸沙率，通過下式計算：式中：H為有效波高；Cg為波群速度；下標b表示波浪破碎時的參數(shù)；θbs為破波角；a1和a2為無量綱參數(shù)，定義如下：

式中：ρ表示海水密度，取1.03×103kg/cm3；ρs為沙的密度，取2.65×103kg/cm3；p為海灘的孔隙率，砂質(zhì)海岸一般取0.4[9]；tan β為海灘平均坡度；1.416為有效波高向均方根波高的轉化系數(shù)；K1、K2為控制沿岸輸沙率大小的調(diào)整參數(shù)，可通過調(diào)整K1、K2值對模型進行調(diào)試和標定。

波浪由外海傳到近岸區(qū)的波浪參數(shù)（如：波高、周期、波向、譜的改變等）采用STWAVE模塊（STeady-state spectral WAVE）進行計算。STWAVE[9]是一個基于波作用平衡方程的有限差分模型，控制方程為：

式中：C為波速；Cg為波群速度；α為波浪正交方向（垂直于波峰線的方向）；μ為波向線方向（能量傳播的方向）；ωr為相對角頻率；ωa為絕對角頻率；E為波能密度；S為能量的源、匯項。

4.1模型設置

模擬區(qū)域計算網(wǎng)格的x軸大致垂直于岸線方向（x軸與正北方向夾角0°），長度為2 000 m（向海側的水深取至5.35 m）；y軸與x軸垂直，長度為1 000 m（包括了整個中海灘在內(nèi)）；網(wǎng)格間距設定為10 m。

模型中涉及的老虎石海灘地形條件較為復雜，為了通過數(shù)值模擬最大化的還原實際條件，模型對老虎石海灘的海岸情況進行了詳細的處理：海灘上的礁石概化為海堤（阻止海灘進一步后退），近岸的礁石概化為離岸堤或離岸潛堤（削減波浪），伸出海岸的礁石概化為丁壩（引起波浪繞射、阻止沿岸輸沙）。根據(jù)實測資料發(fā)現(xiàn)，工程后18個月內(nèi)，人工水下沙壩剖面是穩(wěn)定的，并顯示發(fā)生明顯變形，因此，在模型中概化為離岸潛堤。

4.2模型率定和驗證

為了更準確地預測海灘岸線演變的趨勢，本次建模先對養(yǎng)灘前的岸線演變結果進行模型率定，再對養(yǎng)灘后的岸線演變結果進行模型驗證。

4.2.1模型率定

在岸線長期演變的結果下，認為工程前的岸線形狀處在一個相對比較穩(wěn)定的狀態(tài)。模擬工程前岸線10 a內(nèi)的演變，繪制于圖4?？梢钥吹皆跓o工程狀態(tài)下，10 a間岸線演變模擬結果與初始岸線非常接近，可以認為該模型能夠反映老虎石浴場及周圍海灘的水動力和地貌動力演變規(guī)律。確定下來的GENESIS模型參數(shù)設置如下：閉合水深DC=7 m，灘肩高度DB=2 m，中值粒徑D50= 0.5 mm（采用了客沙粒徑），經(jīng)驗參數(shù)K1=0.8，經(jīng)驗參數(shù)K2=0.96。驅(qū)動模型的波浪資料為工程區(qū)域外海的波浪測站2011年6月—2014年6月的測量結果。

圖4　工程前岸線變遷模擬Fig.4　Simulation of shoreline changes before beach nourishment

4.2.2模型驗證

通過對工程后岸線演變的模擬，對建立的模型進行驗證。選定2013年9月30日和2014年12月23日兩個驗證時段。

首先將實測的工程岸線演變繪制于圖5（a），將模擬的工程后的岸線演變繪制于圖5（b）。如圖5（a）所示，工程后東段海灘在自東向西的沿岸輸沙作用下持續(xù)侵蝕、中段海灘受老虎石陸連島阻擋作用影響持續(xù)淤積、西段海灘因缺乏有效的阻擋持續(xù)侵蝕。這一趨勢在圖5（b）所示的模擬結果中同樣體現(xiàn)出來。

圖5　模擬與實測岸線演變過程比較Fig.5　Comparison of shoreline location between numerical predictions and observations

兩個模擬時段的模擬岸線與實測岸線位置比較繪制于圖6（a）和圖6（b）?？梢钥吹?，在兩個選定的驗證時段，模擬結果與實測結果均符合較好。綜合養(yǎng)灘前和養(yǎng)灘后模型的驗證結果，本次建立的岸線演變數(shù)學模型能夠正確地反應老虎石及鄰近岬灣海灘的演變特征，可以用于工程后岸線演變趨勢的預測。

圖6　模擬與實測岸線不同時期位置比較Fig.6　Comparison between numerically predicted shoreline changes and observed shoreline changes in different periods

4.3模型的預測

利用建立的數(shù)值模型，對老虎石及鄰近岬灣海岸的岸線演變進行預測，模擬時長分別為3 a 和5 a。模擬結果繪制于圖7。

圖7　岸線演變趨勢預測Fig.7　Prediction of shoreline evolution trend

由圖7，工程后3~5 a的時間內(nèi)，老虎石東側海灘岸線形狀變化緩慢，已基本穩(wěn)定。東段海灘在自東向西的沿岸輸沙的作用下，岸線后退，但后退速度變緩；中段海灘堆積了絕大部分來自于東段海灘的泥沙，海灘面積增長，工程5 a后岸線已經(jīng)較完工時向海推進了40 m左右，受淤長后岸線走向的影響，淤長速度會變緩；西段海灘的填沙持續(xù)流失，部分海灘已經(jīng)后退到養(yǎng)灘前位置。

5　結語

本文在實測岸線資料的基礎上對老虎石公園海灘養(yǎng)護工程的養(yǎng)灘效果進行了分析，采用GENESIS數(shù)值模擬的方法構建了岸線演變模型，根據(jù)實測資料和模擬結果的比較對模型進行了率定和驗證，利用該模型預測和分析了工程后3~5 a內(nèi)的岸線演變趨勢，該結果為海灘后續(xù)養(yǎng)護提供了參考依據(jù)。

海灘剖面監(jiān)控結果表明，工程后18個月內(nèi)，在非常明顯的自東向西的沿岸輸砂的作用下，工程填砂發(fā)生了重分布。其中東段海灘岸線后退，中段海灘面積迅速增長，西段海灘泥沙岸線持續(xù)后退。

采用GENESIS數(shù)值模擬的方法構建的岸線演變模型，無工程狀態(tài)下10 a內(nèi)的岸線演變模擬結果與現(xiàn)有岸線非常接近，符合岸線演變的長期趨勢。在養(yǎng)灘后選定的兩個驗證時段，模擬結果與實測結果符合均較好。證明了該模型能夠反映老虎石及周圍海灘的水動力和地貌動力演變規(guī)律，可以用于工程后岸線演變趨勢的預測。

利用構建的岸線演變模型，預測了未來3~5 a內(nèi)，老虎石海灘的岸線演變。根據(jù)模型結果，東段海灘的岸線在沿岸輸砂作用下后退，隨著岸線后退，該段海灘的形狀漸趨平直，后退速度也逐漸減慢，至工程5 a后逐漸達到平衡，海灘的寬度仍明顯寬于養(yǎng)灘前。中段海灘上游接收來自東段海灘的來沙，下游受到老虎石陸連島的阻擋，因此工程后海灘逐漸淤長，工程5 a后海灘寬度相比工程完工時寬40 m左右。西段海灘上游來沙受到老虎石的阻隔，而下游有部分泥沙越過礁石向北戴河西海灘流失，故該段海灘同樣遭受侵蝕，預計直到老虎石的掩蔽作用和礁石的擋砂作用能夠與泥沙的流失相平衡時才能逐漸穩(wěn)定。

由于海灘不同位置的侵淤，海灘整體寬度不均衡，易造成侵蝕部分岸線后退而淤積部分坡度過緩，因此建議對侵蝕部分岸段采取保沙措施或后續(xù)補沙。

[1]鄒志利.海岸動力學[M].北京：人民交通出版社，2009：231-237.

ZOU Zhi-li.Coastal hydrodynamics[M].Beijing:China Communications Press,2009:231-237.

[2]張振克.美國東海岸海灘養(yǎng)護工程對中國砂質(zhì)海灘旅游資源開發(fā)與保護的啟示[J].海洋地質(zhì)動態(tài)，2002，18（3）：23-27.

ZHANG Zhen-ke.Discussions on development and protection of touring resource of Chinese sandy beach from beach nourishment engineering done for American east coast[J].Marine Geology Letters,2002,18(3):23-27.

[3]季小梅，張永戰(zhàn)，朱大奎.人工海灘研究進展[J].海洋地質(zhì)動態(tài)，2006，22（7）：21-25.

JI Xiao-mei,ZHANG Yong-zhan,ZHU Da-kui.Research development of artificial beach[J].Marine Geology Letters,2006,22(7): 21-25.

[4]HANSON H,BRAMPTON A,CAPOBIANCOET M.Beach nourishment projects,practices,and objectives:a European overview [J].Coastal Engineering,2002,47:81-111.

[5]莊振業(yè)，王永紅，包敏，等.海灘養(yǎng)護過程和工程技術[J].中國海洋大學學報：自然科學版，2009，39（5）：1 019-1 024.

ZHUANG Zhen-ye,WANG Yong-hong,BAO Min,et al.Beach nourishment process and engineering technology[J].Periodical of Ocean University of China:Natural Science,2009,39(5):1 019-1 024.

[6]楊燕雄，張甲波，劉松濤.秦皇島海灘養(yǎng)護工程的實踐與方法[J].海洋地質(zhì)前沿，2014，30（3）：1-15.

YANG Yan-xiong,ZHANG Jia-bo,LIU Song-tao.What we have learnt from the beach nourishment project in Qinhuangdao[J]. Marine Geology Frontiers,2014,30(3):1-15.

[7]冒小丹，匡翠萍，潘毅.潮流對老虎石海灘養(yǎng)護工程響應的數(shù)值預測[J].海洋地質(zhì)前沿，2013，29（2）：71-78.

MAO Xiao-dan,KUANG Cui-ping,PAN Yi.Numerical study on hydrodynamics responses to nourishment project of Laohushi beach [J].Marine Geology Frontiers,2013,29(2):71-78.

[8]潘毅.數(shù)值模型研究秦皇島海灘養(yǎng)護工程[M].上海：同濟大學土木工程學院，2009.

PAN Yi.Study on Qinhuangdao beach nourishment project with numerical model[M].Shanghai:College of Civil Engineering of Tongji University,2009.

[9]HANSON H,KRAUS N C.GENESIS:Generalized model for simulating shoreline change[R].Washington DC:US Army Corps of Engineers,1989.

Simulating and forecasting of shoreline evolution of Laohushi Park

XU Yan,PAN Yi*,CHEN Yong-ping
（College of Harbor,Coastal and Offshore Engineering,Hohai University,Nanjing,Jiangsu 210098,China）

This paper analyzes the effect of the beach nourishment project in Laohushi Park based on the field observation of the shoreline.The shoreline evolution model was built using the GENESIS(Generalized Model for Simulating Shoreline change)method,and then the model was calibrated and verified according to the comparison between simulation results and field data.Finally,this model was applied to predict and analyze the tendency of shoreline evolution in 3 to 5 years'time. Results show that the eastern beach shoreline will retreat until gradually achieving balance after 5 years.Middle beachwillsilt gradually,increasing its width by 40 meters.The western beach remains constantly eroded.

GENESIS;beach nourishment;shoreline change

U652.4；P753

A

2095-7874（2016）05-0001-05

10.7640/zggwjs201605001

2015-11-30

2016-01-02

國家自然科學基金（51309092）；江蘇省自然科學基金（BK20130833）

胥巖（1988— ），女，黑龍江綏化人，碩士研究生，主要從事海岸防災研究。

潘毅，E-mail：panyi21hhu@gmail.com