港區建設對入海河流泥沙淤積的影響分析

金子嵩

（遼寧省水利水電勘測設計研究院有限責任公司,遼寧 沈陽 110006）

1 緒論

河流入海口附近河段受潮汐作用明顯,其泥沙淤積特性受所在海域的潮流、波浪、泥沙等特性的綜合作用,因此不能僅僅局限于河道內及入海口,而需要采用適當的潮流泥沙數學模型對受影響的整片海域進行計算。

本文以河北省某港區工程項目為例,建立二維潮流泥沙數學模型,分析港區附近兩條入海河流——陡河和沙河河口區域工程前后的泥沙淤積特性,以及工程建設對河道泥沙淤積的影響。

2 工程概況

某港區工程利用沙河口至陡河口之間岸線形成環抱式港池,港池內東側為突堤結合挖入式港池的形式,西側為順岸布置形式。港內主要布置干散貨泊位區、通用泊位區及臨港工業區、倉儲物流區。沙河河口修建了擋沙堤,與港口圍墾線平行。港區平面布置見圖1。

圖1 港區平面布置圖

3 泥沙來源分析

3.1 河流來沙

豐南港區附近有陡河與沙河兩條小型河流入海,其來水來沙量均很小（表1）,其中陡河全長120 km,在陡河村附近入海；沙河全長108 km,在黑沿子村與黑沿子排干交匯后入海,沙河口位于陡河口以東約10 km處。建閘前兩河口多年平均入海年徑流量和泥沙量分別為1.05 億m3和15 萬t,入海泥沙較少。目前兩河口上游都已經建閘,來水來沙量更是微乎其微,對本工程影響不大。

表1 沙河與陡河多年水沙特征值統計表

3.2 淺灘泥沙

本海區海底平緩,一般為1/1000 ~1/2000 ,灘面主要由細顆粒泥沙組成,在風浪和潮流共同作用下灘地泥沙極易被掀起,并隨漲、落潮流挾帶輸移。大風浪是造成近岸泥沙懸移運動的主要因素。在小風或無風天水體含沙量很低,小于0.1 kg/m3；而隨著風級的增大,其含沙量逐漸增高,近岸最高可測到2 kg/m3~4 kg/m3。

由此看來,淺灘泥沙是影響陡河、沙河河口泥沙的主要來源,要分析河口泥沙淤積特性,需要建立適當的潮流泥沙數學模型。

4 數學模型簡介

本項目的潮流計算采用由丹麥DHI水工所開發的Mike21 軟件的三角形網格水動力模塊（HD模塊）[1],該模塊在國內外工程項目中應用廣泛,模塊計算成果得到大量驗證,受到業內人士普遍認同。MIKE21 軟件的水動力學模塊（HD模塊）基于水流運動控制方程——淺水方程進行計算[2]。

波浪引起的輻射應力作為驅動力影響著潮流流場,同時波浪也是影響泥沙輸運和淤積的一個重要因素。本項目的波浪采用SWAN模型來模擬。SWAN模型基于波作用譜平衡方程進行計算[3]。

泥沙運動采用竇國仁等（1995）基于波流共同作用下挾沙力概念的平面二維泥沙數學模型[4],其表達式如下：

式中：h為水深；t為時間坐標；x和y為水平坐標；S為沿深度平均的含沙量；S*為波流共同作用下的挾沙能力；u和v分別為沿x方向和y方向的流速；為沉降幾率或恢復飽和系數；w為泥沙沉速。根據波流挾沙的原理,S*可近似為：

其中,S*c和S*W分別為潮流和波浪作用下的挾沙能力。潮流作用下的挾沙能力可表示為

式中：βc為根據實驗或者現場資料確定的系數；ρS和ρ分別為泥沙與水的密度；c為謝才系數；V為垂向平均流速。

對于波浪作用下的挾沙能力,根據實際波能演化原理,修正為如下形式：

其中,fw為床面摩阻系數,Hrms為均方根波高,T為波浪周期,k為波數,g為重力加速度,DB2為由于波浪破碎引起的波能耗散, β1與β2為系數。

5 模型布置

對港區范圍內海域建立潮流泥沙模型進行計算,本文采用大、中、小三重模型以嵌套方式進行計算,分別由大模型為中模型、中模型為小模型提供邊界條件。圖2 顯示了大、中、小模型計算范圍。大模型包含整個渤海海域。本工程緊鄰周邊若干其他港口,為充分考慮周邊港口建設對本工程的影響,中模型范圍包含了整個渤海灣,東西向最長約129 km,南北向最長約147 km。小模型為工程局部范圍。

圖2 模型計算范圍示意圖

圖3 給出了中模型以及小模型的計算網格。本模型由三角形非結構化網格構成,其中港區內航道、建筑物及陡河、沙河河道內網格進行加密。

圖3 中、小模型網格示意圖

6 模型驗證

需對工程海域含沙量場進行模擬,率定泥沙模型參數。含沙量的驗證采用2010 年6 月最新實測含沙量資料。圖4 為大潮各測點實測含沙量與計算值的比較情況。

圖4 2010年6月大潮含沙量驗證（0時刻為13日15時）

通過比較分析,各測站的含沙量在連續變化過程中均與實測值接近,其中含沙量平均誤差不超過30%,測點的驗證結果符合要求。因此,本次研究中所建立的泥沙運動模型能夠合理地反映工程區附近海域的含沙量分布,可進一步用于地形沖淤計算。

7 計算成果及分析

在對工程海域波浪場模擬的基礎上,將計算出的波浪動力提供給潮流模型和泥沙模型,對工程海域泥沙運動和地形沖淤進行模擬[5]。潮流計算采用2010 年6 月10 日12 時~2010 年6 月25 日12 時計算潮過程,這其中包含了大、中和小潮連續作用半個月。代表波浪與潮流組合占全年計算時段的23%,純潮流作用時間占77%。兩種動力條件組合得到的回淤計算結果為年回淤計算結果。圖5 和圖6 分別反映了工程前后河口附近的波浪場,圖7 ~圖12 給出了工程前后河口的含沙量場分布圖。

圖5 工程前河口附近波浪場圖

圖6 工程后河口附近波浪場圖

圖7 工程前大潮潮流與代表波浪共同作用下河口附近含沙量場圖

圖8 工程前大潮潮流與代表波浪共同作用下陡河河口含沙量場局部放大圖

圖9 工程前大潮潮流與代表波浪共同作用下沙河河口含沙量場局部放大圖

圖10 工程后大潮潮流與代表波浪共同作用下河口附近含沙量場圖

圖11 工程后大潮潮流與代表波浪共同作用下陡河河口含沙量場局部放大圖

圖12 工程后大潮潮流與代表波浪共同作用下沙河河口含沙量場局部放大圖

從圖中可以看出,工程后陡河與沙河河口及附近海域含沙量明顯較工程前減小,港區以外海域則變化不明顯,這主要是由于港區工程及東側擋沙堤阻擋泥沙輸運所致。

根據本研究建立的針對河口地區的泥沙數學模型知,在只考慮外海來沙的條件下,陡河入海河口段工程前年淤積量在92 萬 m3,工程后年淤積量在57 萬m3；沙河入海河口段工程前年淤積量在179 萬m3,工程后年淤積量在97 萬m3。

綜合考慮泥沙數學模型研究成果,可得到以下結論：（1）在只考慮外海來沙的條件下,由于工程區域的掩護作用,工程建設后引起西側陡河口泥沙淤積減少；（2）本工程為了減小對兩側河口的影響,在沙河河口東側建設擋沙堤,距圍墾堤東線1 km,與圍墾堤線平行。擋沙堤走向與漲落潮流流向一致,對周邊流場的影響較小,比較含沙量分布圖可知,擋沙堤可有效攔截東側淺灘泥沙,圍墾線與擋沙堤之間的含沙量大幅減小,引起沙河河口段的泥沙淤積減小。

兩側河口涉及到當地防洪、漁業等民生問題,需引起足夠重視。雖然計算結果表明,港區建設方案對兩側河口泥沙淤積的直接影響不大,但周邊也在進行大規模圍填,可能造成河口區動力條件進一步減弱,對河口健康的維持不利,需加強監測分析。港區圍填施工過程中,可能會造成含沙量的增加,致使河口淤塞,需及時分析采取相應對策。建議港區開發建設單位對兩側河口進行疏浚,一方面可提高河口防洪航運效益,抵消開發方案的不利影響,另一方面疏浚土可作為港區回填料。

8 結論

本研究采用河北省某港區項目整體潮流泥沙數學模型,對港區附近海域及陡河、沙河河口工程前后的泥沙淤積進行計算。水力學計算模型的采用及概化處理等對類似入海河流泥沙分析具有較為普遍的參考價值。得出結論如下：

（1）陡河、沙河河口上游都已建閘,河流來沙很小,河口來沙以外海來沙為主。

（2）從泥沙淤積計算成果看,由于工程區域的掩護作用,工程建設后引起西側陡河口泥沙淤積減少。沙河河口擋沙堤可有效攔截東側淺灘泥沙,引起沙河河口段的泥沙淤積減小。