吉倫特河口灣
簡介
吉倫特河口灣 Gironde 法國西南部亞奎丹大區吉倫特省比斯開灣一河口灣。由加倫河(Garonne River)以及多爾多涅河(Dordogne River)匯成,長約72公里(45哩)。儘管有沙洲和大潮,仍通行遠洋輪。沿灣波亞克(Pauillac)的貝克達貝斯邁(Bec d'Ambesmay)有煉油廠,濱海勒韋爾東(Le Verdun-sur-Mer)為客運港。
針對法國吉倫特河口灣的大尺度地貌動力學數值模擬研究
前言
一直以來,由於技術條件的限制,在複雜河口灣條件下對泥沙遷移現象進行數值模擬研究經常採取簡化條件的方法以降低對計算資源的需求。而近幾年來,隨着數值模擬方法的不斷優化以及大規模並行計算的使用,現階段的技術能力已經能夠在有限的計算資源條件下儘可能的擴大模擬範圍以涵蓋各種可能的初始和邊界條件。而在本次我們將介紹一個利用法國電力自由水力學軟件TELEMAC-MASCARET對法國吉倫特河口灣進行大尺度地貌動力學和泥沙遷移的數值模擬案例。
研究背景
本次研究的目的為初步建立一個基於真實條件的大尺度地貌動力學模型以對大型河口灣地區的泥沙遷移以及中期(5年)河床演變進行相對準確的預測。該模型最終可被用於測試各類解決方案以防止河口灣地區的一些關鍵設施遭到泥沙遷移的侵蝕和堵塞等影響。
此次研究的區域是位於法國西部的吉倫特(Gironde)河口灣。作為法國乃至歐洲西部沿海最大的河口灣,吉倫特是多爾多涅河(Dordogne)和加龍河(Garonne)的匯合點以及入海口,其寬度從3到11公里不等,而縱深有將近80公里,總面積達到635平方公里。受到兩條主幹河流的影響,其累計流量一般在50到2000 立方米/秒不等,在洪澇季節甚至可超過5000立方米/秒。其複雜的河口條件以及巨大的面積使得該地區在數值模擬方面需要大量的計算資源,因此對該地區泥沙遷移以及地貌動力學方面的研究直到最近都一直難以推進。
本次的研究將涵蓋整個吉倫特河口灣,在其中將選取多個測量點以在模型輸出結果和觀測數據之間做出一個初步的比對。研究工具是法國電力自主開發的自由水力學軟件TELEMAC-MASCARET。其中為平衡模型精度和計算資源要求,主要利用2D模型來進行研究。
模型建立
水力學模型
水力學部分將由TELEMAC-MASCARET系統中的二維水力學模塊TELEMAC-2D完成建模。模型的計算範圍將覆蓋整個吉倫特河口,包括上遊河段、中央匯流區和入海口。TELEMAC-2D中所使用的非結構有限元網格可以更加適應大尺度建模的需要,同時在重點區域(例如在本案例中的中央匯流區)可以進行較為靈活的網格加密以反映該區複雜的物理現象。另外深度平均淺水控制方程的使用也可以節省更多的計算資源。
泥沙遷移模型
泥沙遷移部分的建模由同樣包含於TELEMAC-MASCARET系統的Sisyphe模塊完成。該模塊的核心控制方程為埃克納方程(Exner equation),其廣泛用於反應泥沙等沉積物的質量守恆。在泥沙遷移的計算過程中Sisyphe會將整個泥沙部分分成河床底積以及懸移質,前者可通過半經驗公式Meyer-Peter Muller 來估計,而後者將通過解算另一個遷移方程以獲得深度平均泥沙含量。
模型耦合
鑑於泥沙遷移現象主要發生在河流近床部分,模型耦合過程中河流底部的速度將被作為主要參數之一參與泥沙遷移的計算。而考慮到基於深度平均淺水方程的TELEMAC-2D的輸出為深度平均流速,因此需要特殊的校正處理。主要方式是通過在解算過程中加入一個基於空間變化的位於[0,1]區間之內的校正係數來人為降低河流的遷移速率,當然這種處理會使模型在數值解算過程中產生不穩定性,因此在水平移流解算過程中加入一個隱式的沉降速率項以避免計算崩潰。事實上由於這種方法已被作為一種常用工具集成在TELEMAC-MASCARET系統當中,因此實際操作起來非常方便。
模擬結果
大尺度水力學模型的驗證
首先要驗證大尺度下水力學模型的有效性。驗證地點選取入海口附近的Verdon(點1)站點以及中央地區的Pauillac(點4)站點。
大尺度泥沙遷移模型的驗證
1. 粒度分布驗證
由於泥沙在遷移過程中涉及多個粒度的混合,我們將首先驗證粒度分布的計算和觀測數據,在網格密集的中部地區粒度分布的計算結果與觀測數據非常契合。
2. 中期河床演變
在中期河床演變部分,我們將展示吉倫特河口中央匯流區域的河床演變情況。時間方面選取1995-2000年5年的變化情況,計算結果與觀測數據的比對如下圖所示:
左圖為1995年-2000年5年內河床演變的觀測數據,而右圖為針對該時間段的計算結果。由此可見計算結果在定性與定量方面都比較契合觀測數據。
研究結論
本次研究利用TELEMAC-2D與Sisyphe的耦合完成了對法國吉倫特河口灣的初步模擬。現階段模型對泥沙遷移和中期(5年)河床演變的模擬在定性層面上能夠與觀測數據達到很好的一致性。而定量方面,中央匯流區的模擬結果相對較好,而入海口和上游部分由於計算資源限制模擬效果相對有限。鑑於吉倫特河口灣巨大的面積,複雜的地理,水文形勢以及泥沙成分,達到完全契合的定量分析相對較難。因此對於本次研究來說,對大尺度河口灣的泥沙遷移和中期河床演變模擬已初步取得一定的效果。為獲得更加精細的模擬結果,在將來的工作中需要進一步完善各類參數的設置,以最終達成對大尺度河口灣更加契合實際情況的模擬效果。