膠州灣大橋建設前后灣內泥沙沖淤數值模擬

張永強，張菀君，遲萬清，邊淑華，曹成林，胡澤建，劉建強

(自然資源部第一海洋研究所，山東 青島 266061)

當前，我國主要經濟帶主要分布于東部沿海地區，跨海大橋的建設為我國經濟建設做出了重大的貢獻。為了加強東部沿海城市與島嶼之間的聯系，近年來我國跨海大橋的建設取得了空前的成就，但存在的問題是，跨海大橋的建設對研究海域的海流流態、泥沙沖淤及生物的生態環境等造成了一定的影響，因此對所屬海域大橋建成前后的環境評估至關重要。

王晨陽等(2010)基于無結構網格的港珠澳大橋所在海區平面二維潮流數學模型，分析了港珠澳大橋工程周圍海域的潮流動力影響[1]。何杰等(2012)采用平面二維潮流數學模型模擬了港珠澳大橋對珠江口水域水動力的影響發現大橋工程對珠江口水域的水動力分布格局基本不產生影響[2]。盛天航(2015)模擬了順河跨海橋工程對河道行洪的影響，指出跨海大橋的建立是導致河道大渦旋的直接成因，河道道行洪所所造成的河床的沖刷，經過計算后得出其最大沖刷的深度滿滿足對特殊大橋的基底埋深的安全值[3]。劉波等(2016)基于GOCI數據對杭州灣跨海大橋兩側水域懸浮泥沙濃度空間分異規律進行了研究，指出杭州灣大橋兩側懸浮泥沙濃度呈現一定的梯度特征[4]。石曉雨等(2018)基于區域海洋水動力模型ECOM的模擬結果分析了近年來膠州灣典型水文氣象要素的變化特征及建橋前后水動力環境等的變化對冬季冰情的影響，指出跨海大橋建設從多方面影響膠州灣北部海冰的生消[5]。

1 研究區域與方法

1.1 研究區概況

膠州灣位于中國黃海中部(36.06°～36.25°N，120.10°～120.37°E)，為一淺水海灣，總體上呈簸箕形，直傾斜在灣口區又轉而向東傾斜，灣內平均水深約7 m，最大水深64 m[6]。

1.2 海流觀測

2013年11月19日至2013年12月12日期間，在研究海域分別進行了大潮期間海流多船定點同步連續觀測，站位分布見圖1。

圖1 膠州灣海流觀測站位示意圖Fig.1 Position of ocean current observation stations in Jiaozhou Bay

1.3 研究區泥沙沖淤數值模擬

本研究運用MIKE21軟件從而達到有效的設計條件和參數[7-8]，并將沉積物分析、洋流觀測等方法與水深地形和工程地質數據相結合來模擬潮流作用下研究區周圍海底地形的演變。以2013年11月19日水文觀測期間的潮型的開邊界為計算條件，計算時采用海流觀測期間的岸線，并對當時的水文觀測資料進行了驗證。為了減少其他膠州灣沿岸工程對計算結果的影響，在考慮膠州灣大橋建設前計算情景時以2013年11月19日水文觀測期間的岸線及水深作為計算岸線和水深，與建設預測情景相比增加了橋墩，并且根據實際尺寸將橋墩設計為模型中的防滲結構。

1.3.1 泥沙運動控制方程 泥沙控制方程為：

(1)

1.3.2 沉積物沉積和侵蝕的計算公式 ①沉積速率。根據以前的研究結果提出的方法來計算沉積速率(SD)，公式如下：

SD=ws·cb·pd

(2)

式(2)中：ws為沉降速度；cb為底層懸浮泥沙濃度；pd為淤積速率的斜坡函數。沉降速度計算公式如下：

(3)

式(3)中：k,γ為系數；c為體積濃度；ws,r為沉降速度系數；cgel為泥沙絮凝點；ws,n為組分能量常數。

Pd的計算公式如下：

(4)

式(4)中：τb為底床剪切力；τcd為臨界淤積剪切應力。

②泥沙濃度的分布。泥沙濃度分布有2種計算方法，其一是Teeter公式，如下：

(5)

其二如下：

(6)

式(6)中：C為懸浮泥沙濃度；R為Rouse參數，計算公式如下：

(7)

式(7)中：ε為擴散系數；z為垂向笛卡爾坐標；Ca為深度基準面處的懸浮泥沙濃度；a為深度基準面。

③底床侵蝕。底床侵蝕根據底床的密度，有2種計算侵蝕的方法，一是壓實固結底床侵蝕的計算公式，如下：

(8)

式(8)中：SE為侵蝕率；E為侵蝕系數；τce為臨界沖刷剪切應力；n為侵蝕能力。

二是軟、部分固結底床侵蝕計算公式，如下：

SE=Eexp[α(τb-τce)1/2],τb>τce

(9)

式(9)中：α為參考系數。

1.3.3 輸入參數確定 水深地形：以研究海域最新實測地形并結合航保部公布的海圖數字化的數據提供模型計算水深；計算岸線：與水動數模型計算岸線給定方式一致；沉積物的類型以及粒度參數特征：我們依據研究區最近和歷史沉積物的表層調查資料，2013年8月研究區沉積物調查結果表明，在膠州灣灣頂發育了大片潮灘，研究區的潮灘以及海底的沉積物的類別主要有粉砂T、砂質粉砂ST及粉砂質砂TS，中值粒徑在1.8～7.8 φ之間(圖2)；在大橋中部、2 m水深以深的海域，沉積物顆粒較細，主要成分為粉砂，中值粒徑在5.9～6.3 φ之間，沉積物分布較均勻；在大橋的東西兩側接陸端、1 m水深以深的海域沉積物分布相對偏粗，主要是砂質粉砂，中值粒徑在3.7～4.9 φ之間；其他區域包括滄口水道海域，沉積物分布比較均勻，主要成分是粉砂，顆粒較細，中值粒徑在5.2～7.8 φ之間。

圖2 中值粒徑φ值分布Fig.2 Distribution map of median particle size (φ)灰色線為等深線，藍色線為中值粒徑等值線，紅色線為膠州灣大橋。

2 結果與討論

2.1 研究區潮流場數值模擬以及對模型的驗證

本研究中數值模型采用的是二維平面潮流數值模型—Mike21-FM模型[9-11]，將模擬結果與研究區各站位的實測海流資料進行比較(圖3、4)，得出二者的趨勢基本一致，準確地說明了大橋建成前后的潮流場分布情況。

圖3 L04站潮位驗證Fig.3 Verification chart of tide level at L04 station

圖4 各站位潮流驗證Fig.4 Verification chart of tides at each station

2.2 膠州灣大橋建設前后潮流計算結果比較分析

圖5 落急時刻膠州灣大橋施工前后流場比較圖Fig.5 Current fields during the ebb before and after the bridge construction工程前流矢量圖為紅色，工程后流矢量圖為黑色。

由膠州灣大橋建設前后落急時流場比較圖(圖5～9)可以看出，落急時流向改變較大的海域主要是在膠州灣大橋橋墩附近。總體特征是：相對于施工前，膠州灣跨海大橋建成后的流向在東部海區域沿逆時針方向偏轉，而在西部海區域則沿順時針方向偏轉。流速改變表現為在膠州灣大橋連接陸側海域中，在大橋的北側，大橋建設后流速較大橋建設前變小；在大橋南側，大橋建設后流速較建設前增大;在大橋東側通航孔處，兩個最大跨度橋墩之間，工程后流速較工程前增大幅度較大，可達0.06 m/s以上；東西兩側的海域中，工程后的流速相對于工程前變小，下降幅度約為0.02 m/s。從各個橋墩附近的海域中可以看出，在橋墩南北兩側，工程后流速相對于工程前變小。在大橋東西兩側的海域中，施工后流速相比工程前變大，但變化幅度不大。落急時在膠州灣大橋與紅島連接線的東西兩側海域，工程后的流速相對于工程前變大，變大幅度在0.02 m/s左右。在大橋中部通航孔附近，跨度最大的橋墩之間的海域，工程后的流速相對于工程前變大，而在東西兩側海域工程后流速相對于工程前變小，變小幅度在0.02 m/s左右。在西側通航孔附近，跨度最大的橋墩之間海域，工程后流速相對于工程前變大，而在東西兩側海域工程后流速相對于工程前變小，變小幅度在0.04 m/s左右。

由膠州灣大橋建設前后漲急時流場比較圖(圖10～14)可以看出，漲急時流向改變較大的海域主要是在膠州灣大橋橋墩附近，總體特征是：施工后，與施工前相比流向發生了逆時針偏轉。流速改變表現為：大橋東側通航孔附近，兩個最大跨度墩之間的海域，工程后流速相對于工程前變大，變大幅度在0.06 m/s左右；在各個橋墩的南北兩側海域，工程后流速均相對于工程前變小，而東西兩側工程后流速相對于工程前變大，但存在流速改變的海域范圍不大。膠州灣大橋與紅島連接線附近海域，工程后流速相對于工程前變小，變小幅度在0.02～0.14 m/s之間。在大橋中部通航孔附近，跨度最大的橋墩之間海域，工程后流速相對于工程前變大，而東西兩側海域工程后流速相對于工程前變小，變小幅度在0.02 m/s左右。在西側通航孔附近，跨度最大的橋墩之間海域，工程后流速相對于工程前變大，而東西兩側海域工程后流速相對于工程前變小，變小幅度在0.04 m/s左右。

圖6 落急時刻膠州灣大橋施工前后流速比較等值線Fig.6 Current velocity contours during the ebb before and after the bridge construction紅色區域表示流速增加，藍色區域表示流速減少；框1、2、3分別為圖7、8、9展示區域。

圖8 落急時膠州灣大橋中部通航口附近流速比較等值線Fig.8 Current velocity contours during the ebb near the central navigation port紅色區域表示流速增加，藍色區域表示流速減少。

圖9 落急時膠州灣大橋西側通航口附近流速比較等值線Fig.9 Current velocity contours during the ebb near the west navigation port紅色區域表示流速增加，藍色區域表示流速減少。

圖10 漲急時刻膠州灣大橋工程前后流場比較Fig.10 Current fields during the flood before and after the bridge construction工程前流矢量圖為紅色，工程后流矢量圖為黑色。

圖11 漲急時刻膠州灣大橋工程前后流速比較等值線Fig.11 Current velocity contours during the flood before and after the bridge construction紅色區域表示流速增加，藍色區域表示流速減少；框1、2、3分別為圖12、13、14展示區域。

圖12 漲急時膠州灣大橋東側通航口附近流速比較等值線Fig.12 Current velocity contours during the flood near the east navigation port紅色區域表示流速增加，藍色區域表示流速減少。

圖13 漲急時膠州灣大橋中部通航口附近流速比較等值線Fig.13 Current velocity contours during the flood near the central navigation port紅色區域表示流速增加，藍色區域表示流速減少。

圖14 漲急時膠州灣大橋西側通航口附近流速比較等值線Fig.14 Current velocity contours during the flood near the west navigation port紅色區域表示流速增加，藍色區域表示流速減少。

2.3 膠州灣大橋建設后沖淤計算結果

由膠州灣大橋建設后年沖淤變化圖(圖15)可知，在滄口水道與紅島連接線之間的南側海域，工程后處于沖刷區，沖刷強度最大可達0.015 m/a；在大橋連接青島側的岬灣處表現為淤積環境，淤積強度為0.050 m/a左右；在中部通航孔與紅島連接線之間海域表現為沖淤平衡區，淤積強度在-0.004～0.004 m/a之間，在西側通航孔附近表現為沖刷區，沖刷強度可達0.013 m/a左右。

圖15 膠州灣大橋建設后年沖淤變化Fig.15 Annual scouring and silting changes after the bridge construction

在膠州灣大橋建設后東側通航孔附近，由于大橋建設后的阻流效應，在跨度最大的兩個橋墩之間表現為明顯的沖刷區，沖刷強度最大可達0.030 m/a；其他小型橋墩表現為：在南北向的兩個小橋墩之間為淤積區，淤積強度為0.050 m/a以上，在東西向的平行橋墩之間則為沖刷區，沖刷強度在0.015 m/a左右。

膠洲灣大橋與紅島連接線附近向陸側海域明顯表現為淤積區，在橋墩附近淤積強度為0.080 m/a左右，在橋墩之間淤積強度為0.004 m/a左右；在南側海域主要表現為沖刷區，沖刷強度為0.006 m/a。

在中部通航孔附近，除跨度最大的兩個橋墩之間表現為沖刷區(年沖刷強度在0.017 m/a左右)外，其他大部分海域表現為淤積區，且橋墩附近海域淤積強度可達0.050 m/a以上。

在西側通航孔附近，跨度最大的兩個橋墩之間和通航孔北側海域表現為沖刷區，沖刷強度分別在0.020 m/a左右和0.009 m/a左右，通航孔南側海域則表現為沖淤平衡區。

2.4 膠州灣大橋建設前沖淤計算結果

由膠州灣大橋建設前年沖淤變化圖(圖16)可知，大橋建設前后研究海域在大范圍上的沖淤趨勢上基本一致。在東側通航孔附近，部分海域表現為沖刷區，沖刷強度在0.007 m/a左右。

在大橋與紅島連接線之間的南側海域表現為沖刷區，沖刷強度在0.007 m/a左右，北側海域表現為弱沖刷區，沖刷強度在0.001 m/a左右。

在中部通航孔附近，大部分海域表現為淤積區，其中跨度最大的兩個橋墩之間淤積強度在0.030 m/a左右，其他區域淤積強度可達0.010 m/a左右；在中部通航孔與紅島連接線之間海域表現為沖淤平衡區。

在西側通航孔附近海域表現為沖刷區，沖刷強度在0.007 m/a左右，大橋北側的沖刷趨勢較大橋南側稍強；在大橋連接青島側的岬灣處表現為淤積區。

3 結論

由膠州灣大橋建設后各典型海域如東側通航孔、紅島連接線附近、中部通航孔、西側通航孔等的潮流場數值模擬結果可知，在通航孔處跨度最大的橋墩之間，流速略高于周圍的海域。與周圍海域相比，小橋墩南側和北側的流速相對較小，東西兩側海域的速度大于周圍海域的速度。在連接線附近海域表現為：落急時橋墩南側海域流速較周邊海域小，漲急時橋墩北側海域流速較周邊海域小。

通過對膠州灣大橋建設前的流場進行數值模擬可知，膠州灣大橋的建設對膠州灣整體的流態影響較小，漲落急等特征時刻流場分布在膠州灣大橋建設前后基本一致。

膠州灣大橋的建設對膠州灣灣內水動力的改變主要在橋墩建設附近海域，主要表現為在漲落急特征時刻，橋墩南北兩側的流速比工程前要小，橋墩東西兩側的流速比工程前要大。特別在跨度較大的通航孔附近海域，橋墩東西兩側流速與工程前相比增大幅度較大。

圖16 膠州灣大橋建設前年沖淤變化Fig.16 Annual scouring and silting changes before the bridge construction

通過對大橋建設前后沖淤數值進行模擬，結果可知，工程前后工程附近沖淤趨勢基本一致，只是在大橋橋墩附近出現的沖淤變化較大的情況，表現為在橋墩的南北兩側出現較大橋建設前淤積加強區，在橋墩的東西兩側出現沖刷加強區，在通航孔橋墩跨度較大的海域，沖刷強度增大比較明顯。