麗水36-1終端至交付點輸氣管線工程海床沖刷研究

孫國民,韓志遠,李文丹,郎一鳴,劉志剛,龐啟秀

（1.海洋石油工程股份有限公司，天津529228；2.交通運輸部天津水運工程科學研究所工程泥沙交通行業重點實驗室,天津300456）

麗水36-1終端至交付點輸氣管線工程海床沖刷研究

孫國民1,韓志遠2,李文丹2,郎一鳴1,劉志剛1,龐啟秀2

（1.海洋石油工程股份有限公司，天津529228；2.交通運輸部天津水運工程科學研究所工程泥沙交通行業重點實驗室,天津300456）

利用已有大量實測資料對工程海域的水動力泥沙環境及工程海域特別是海底管線路由位置的海床穩定性進行了研究。建立了工程海域考慮波浪作用的二維潮流泥沙數學模型，對考慮周圍近期及遠期工程建設影響下的灘面沖淤進行了計算，并給出了最終的沖淤分布。

輸氣管線；海床沖刷；數值模擬；潮流；泥沙；麗水

浙江溫州地區是我國經濟發展較快的地區，對天然氣的需求增長快速，為滿足其對天然氣日益增長的需求，中海油集團擬在洞頭縣霓嶼島上建設麗水36-1氣田陸上終端，將霓嶼島終端分離出的天然氣通過海底管道輸送至溫州龍灣分輸站。輸氣管道設計管徑12英寸，路由全長約28.8 km，在霓嶼島南山的終端站入海，向西經甌江南口水道海域至龍灣海濱城附近登陸（圖1）。海底管道建設和運行過程中，其安全性直接影響到海洋油氣資源開發和海洋環境保護。大量工程實踐中，海床不穩定特別是局部沖刷將會導致海底管道懸空，是影響海底管道安全的重要因素［1-4］。影響海床穩定性的主要因素包括海域水動力特征、海域泥沙環境、海床沉積特性、人類活動引起的海床邊界條件等［5-9］。

擬建海底管道位于甌江南口海域，該海域為強潮海區，潮差大、潮流急、島嶼林立、淺灘密布、灘槽交錯、地形復雜，且受甌江徑流及波浪潮流共同作用，水動力泥沙環境較為復雜。本區域由于經濟快速發展，人類活動影響頻繁（圖1）：已建工程有靈霓大堤、溫州淺灘一期工程、藍田作業區進港航道；在建工程有溫州淺灘二期促淤堤；擬建工程有甌飛邊灘圍墾工程、甬臺溫高速復線大橋、甬臺溫管線等。這些工程建成后將導致該區域邊界條件發生大幅變化，因而該區域的水動力泥沙環境及海床沖淤趨勢也將隨之發生變化，這也必將影響該區域海床的穩定性，對海底管道的安全產生影響。因此在海底管道建設之前需充分掌握工程海域的海床沖淤特性及演變趨勢。

甌江口海域已有研究主要著眼于港口開發、航道整治、淺灘圍墾促淤等，且集中在甌江北口水道、溫州淺灘等水域［10-14］，目前對于南口水道海床穩定性及其對海底管道的影響鮮有研究。本文在充分研究甌江南口水道海域水動力泥沙環境的基礎上，通過建立考慮波浪作用的二維潮流泥沙數學模型，研究相關工程實施后對海底管道沿程海床沖淤變化情況，為麗水36-1海底輸氣管道的設計提供基礎依據。

圖1 工程示意圖Fig.1 Sketch of project position

1 自然條件分析

1.1 河口地貌

甌江是浙江省第二大河，發源于浙南慶元縣仙霞嶺，在溫州市龍灣區和樂清市交界處入海。甌江為源短流小的山溪性河流，多年平均入海徑流總量為148億m3，多年平均入海懸沙總量為205.1萬t。每年4～9月為洪季，10月～翌年3月份為枯季。甌江入海口被靈昆島分為南口、北口，河口東側為洞頭列島，排列多個島嶼。甌江河口為強潮河口，在徑流和潮流共同作用下，甌江口灘槽沖淤變化較為劇烈，發育有多個江心淺灘和島嶼如七都島、靈昆島等。甌江南口水域北側為靈昆島及溫州淺灘，東側為洞頭列島，西側為甌飛邊灘，東南側接霓嶼深槽，南側為溫州灣開闊海域。

1.2 風和波浪

據洞頭氣象站1971～2001年風況資料，本地區全年常風向為N—NE向，頻率為54.7%，全年平均風速為3.8 m/s；強風向為SSW向，最大風速為32.0 m/s。全年中夏季多為SW向大風，春秋季節多為偏S向或偏N向大風，以偏N向大風為主，冬季盛行N—NE向大風。年平均≥6級風日數為37 d，年平均≥7級風日數為8.5 d。該地區每年5～11月是臺風影響期，其中7～8月份為臺風盛行期，有嚴重影響或在本地區登陸的大約每兩年一次。甌江南口水域被溫州陸地、靈昆島—靈霓大堤、洞頭列島圍成一個僅向南開敞的海灣，因此NNE、NE向入射波浪幾乎全部被阻擋，E向波浪的影響也很小；只有ESE向、SE向、SSE向及S向波浪對南口水道的影響最直接。據南口水道附近2001～2004年實測資料，小于0.5 m波高頻率為91.9%，因此工程海域波浪作用較小。

1.3 潮汐和潮流

1.4 底質和含沙量

甌江南口附近表層沉積物為砂質沉積物，中值粒徑大于0.15mm；靈昆島至霓嶼附近底質以粉砂質粘土為主，局部有粘土質粉砂分布，中值粒徑基本小于0.004mm，僅霓嶼西北側及甌飛淺灘南側在0.004～0.008mm范圍內。工程海域海床性質屬于淤泥質范疇。甌江南口附近含沙量較高，大潮漲、落潮實測平均含沙量分別為2.51 kg/m3和2.40 kg/m3，南口水道南部含沙量較低，大潮漲、落潮實測平均含沙量介于0.23～0.40 kg/m3。可見，甌江南口水域含沙量較高，且由河口向外海含沙量逐漸降低。工程水域的主要泥沙來源為甌江河流來沙、外海來沙及淺灘局部再搬運的泥沙。

1.5 海床沖淤演變趨勢

根據前人研究成果，1999年之前，甌江南口水域海床略呈淤積趨勢［15］。據1999年之后海底管道沿程地形沖淤分析可知（圖2）：①1999～2005年間，南口水道深槽內沖淤相間，整體呈沖刷趨勢，最大沖刷深度超過0.6m；小霓嶼西側呈淤積趨勢，淤積幅度基本在0～0.4m，局部超過0.6m；小霓嶼南側呈沖刷趨勢，沖刷幅度基本在0～0.2m；霓嶼南側呈淤積趨勢，淤積幅度基本在0.2～0.6 m。②2005～2011年間，南口水道北段呈淤積趨勢，淤積幅度超過0.4m；淺灘一期西側及南側為藍田航道開挖段，水深普遍增加超過4m；小霓嶼西側及南側淺灘沖刷趨勢，沖刷幅度基本在0～0.4m；霓嶼南側沖淤相間，沖淤變化幅度基本在-0.2～0.4m。

整體而言，2002年以后，隨著靈霓大堤（2003～2006年修建）、淺灘一期南圍堤和東圍堤的建成（2006～2009年修建）、藍田作業區航道建設（2010年1月竣工）及淺灘二期圍涂促淤堤開工（2010年6月開工，工期2 a），該水域水動力泥沙環境發生了變化。溫州淺灘水域由于水動力環境減弱，呈持續淤積狀態；而南口水道北段沖淤基本平衡，南段局部沖刷幅度較大，整體呈基本穩定略有沖刷的狀態；霓嶼深槽水域淺灘地形呈沖刷趨勢，但近年來沖刷趨勢有所減緩。

圖2 輸氣管線路由沖淤分布Fig.2 Distribution of siltation and erosion in position of gas transmission pipeline

2 二維潮流泥沙數值模擬

2.1 模型建立

考慮波浪作用的二維潮流泥沙數學模型基本方程及數值方法參見文獻［15］。

根據工程的特點，本次計算所建立的數學模型東邊界至121°35′經線，北邊界至樂清灣頂，南邊界至南麂島南側。計算域的東西距離為100 km，南北距離120 km。模型開邊界由Chinatide提供［16］，甌江河口徑流采用流量控制。

2.2 網格劃分

為擬合復雜岸線和橋墩等細致建筑物邊界，數學模型中采用無結構三角形網格對計算域進行剖分。圖3給出了計算所采用的網格示意圖，管線所在路由位置進行了局部加密。用于最終計算的模型約39 813個網格節點，其中最小網格尺度近1～2 m。大橋橋墩群為透水建筑，一方面由于樁基阻力的影響，流速將減小，另一方面，又因樁柱體存在而使過水斷面縮小，流速將增加，此外，由于樁的存在，水流中還會形成旋渦，樁群內的水流非常復雜，因此，在有樁基的潮流數值模擬中，應該考慮橋墩的影響，本文通過局部網格加密將橋墩處理為陸地。圖4為橋墩局部網格示意圖。

圖3 計算網格Fig.3 Computational mesh

圖4 橋墩網格示意圖Fig.4 Sketch of computational mesh of pier

2.3 模型驗證

為了驗證模型的合理性，采用2006年10月、2009年2月和2010年10月3次甌江口水文全潮測驗（每次測驗均為大、中、小3個潮型）共9個潮次的資料進行了潮位、流速流向和含沙量的驗證。驗證結果表明，計算值與實測值吻合良好。

海底地形沖淤變化的驗證是本次計算的基礎，考慮了大、中、小潮、各向向岸波分級波浪（根據計算域內南麂島國家海洋局長期測波站1994年全年波浪資料結合本工程的特點，選擇ESE向和由SE向、SSE向、S向加權所得的代表波浪，并對其波高、波向和波周期進行了驗證）及徑流的組合。本文采用溫州淺灘附近水域最新地形資料對模型進行灘面沖淤驗證；采用藍田航道2010年1月～2011年7月兩次實測水深測量資料對航道內淤積進行了驗證；另外，根據淺灘一期西圍堤西側及南側沖刷趨勢，本模型進一步對其沖刷范圍及深度進行了驗證，驗證結果同該區域2005～2010年間的地形沖淤趨勢的分析結果基本一致。

3 海床沖淤預測

3.1 計算方案

由前言可知，麗水36-1輸氣管線所在的溫州南口水域工程頻繁，在建工程有溫州淺灘二期圍涂促淤堤；擬建工程有甌飛邊灘圍墾工程、甬臺溫高速復線大橋、甬臺溫管線。大規模的工程建設將會對管線路由位置海底地形產生一定的沖淤影響，因此，有必要對周邊工程建設后工程海域地形進行模擬計算。周邊工程分近期和遠期兩種工況考慮。近期：甌飛邊灘不實施+淺灘二期考慮實堤形式+甬臺溫高速復線大橋實施；遠期：甌飛邊灘實施+淺灘二期考慮實堤形式+甬臺溫高速復線大橋實施。

3.2 工程后潮流場

采用經驗證的潮流數學模型對近期和遠期兩種工況下的潮流場進行了計算，有以下結論：

（1）工程所在海域漲、落潮流受島嶼及陸地岸線限制，基本呈往復流運動趨勢。漲潮時，外海潮流向岸灘及甌江河口方向運動；落潮時，甌江的徑流及納潮水體下泄，大部分水流經大門島南側北口水道，少部分經南口水道與大門島北側水道入海。

（2）近期和遠期工況實施后，并沒有改變大范圍海域潮流總體運動規律，僅工程局部海域潮流有所變化。近期工程實施后，淺灘二期西側局部水流流速增加，大小霓嶼島南側水域流速減小。遠期工程實施后，工程區大部分水域流速呈減小趨勢，甌飛邊灘與淺灘二期工程之間水域束窄，流速呈增加趨勢。

圖5 漲落急流場圖(原型)Fig.5 Flow field at flood and ebb strength of tide(prototype)

圖6 漲落急流場圖（方案）Fig.6 Flow field at flood and ebb strength of tide(scheme)

3.3 輸氣管線路由位置沖刷研究

采用經過驗證的二維潮流泥沙數學模型對近期和遠期兩種工況實施后灘面沖淤達到平衡情況進行了計算（圖7和圖8）。由圖可見：（1）由擬建工程位置可知，輸氣管線與高速復線大橋相交叉，為避免橋墩局部沖刷，管線路由于大橋橋墩中間穿過，管線兩側橋墩間距為50m。工程實施后，橋孔位置流速增加，流速的增加，勢必引起管線路由位置沖刷，經計算，其沖刷深度約為1.1m。（2）近期工程實施后：小霓嶼—大霓嶼附近水域為淤積狀態；小霓嶼西側為沖刷狀態；電纜與輸氣管線交匯位置處為2.0m以下沖刷狀態。（3）遠期工程實施后：小霓嶼—大霓嶼附近水域為淤積狀態；小霓嶼西側水域為沖刷狀態（輸氣管線位置沖刷小于0.5 m）；電纜與輸氣管線交匯位置處水域為0.5～2.0m沖刷狀態；甌飛邊灘與二期工程之間水域為2.0～3.0m沖刷狀態。

圖7 近期工程實施后灘面沖淤Fig.7 Distribution of siltation and erosion after recent project

圖8 遠期工程實施后灘面沖淤Fig.8 Distribution of siltation and erosion after future project

4 結論

本文在分析工程海域的水動力泥沙環境及海床穩定性的基礎上，通過建立二維潮流泥沙數學模型，計算了近期及遠期工程建設影響下的海床灘面沖淤變化，給出了灘面最終平衡狀態下的沖刷深度，得出主要結論如下：（1）甌江南口海域為強潮海區，流速大、含沙量高。海床底質以粘土質粉砂為主，屬淤泥質海岸性質。（2）近期溫州淺灘呈淤積狀態，南口水道呈基本穩定并略有沖刷的狀態，霓嶼深槽水域呈緩慢沖刷趨勢。（3）甬臺溫高速復線大橋實施后，輸氣管線路由所在橋孔位置沖刷深度約為1.1m。（4）近期工程實施后，小霓嶼西側沖刷深度介于2.0～3.0m，電纜與輸氣管線交匯位置沖刷深度小于2.0m。（5）遠期工程實施后，小霓嶼西側水域沖刷深度小于0.5 m，電纜與輸氣管線交匯位置處水域沖刷深度介于0.5～2.0 m，甌飛邊灘與淺灘二期工程之間水域沖刷深度介于2.0～3.0m。

［1］喻國良,陳琴琴,李艷紅.海底管道防沖刷保護技術的發展現狀與趨勢［J］.水利水電技術,2007,38（11）:30-33. YU G L,CHEN Q Q,LIY H.Status and tendency of development of scour-prevention technique for submarine pipeline［J］. Water Resources and Hydropower Engineering,2007,38（11）:30-33.

［2］許文兵.平湖油氣田海底輸油管線臨時修復方法［J］.中國海洋平臺,2003（2）:39-42. XUW B.The Temporary Repairing Method for Sub-sea Oil Pipeline of Ping-hu Oil and Gas Field［J］.China Offshore Platform, 2003（2）:39-42.

［3］苗文成.埕島海域海底管道隱患分析及治理［J］.石油工程建設,2004（3）:48-50. MIAOW C.Analysis and Settlement of Hidden Trouble of Pipeline on Chengdao Seabed［J］.Petroleum Engineering Construction, 2004（3）:48-50.

［4］王利金,劉錦昆.埕島油田海底管道沖刷懸空機理及對策［J］.油氣儲運,2004（1）:44-48. WANG L J，LIU JK.The Scouring Span and its Control Methods for Subsea Pipeline in Chengdao Oilfield［J］.Oil&Gas Storage and Transportion,2004（1）:44-48.

［5］曹祖德,肖輝.潮流作用下淤泥質海床沖淤演變預測及應用［J］.水道港口,2009,30（1）:1-8. CAO Z D,XIAO H.Seabed evolution prediction and application on muddy coast under action of tidal current［J］.Journal of Waterway and Harbor,2009,30（1）:1-8.

［6］黃世昌,張舒羽,余炯.杭州灣灰鱉洋海域海床演變趨勢研究［J］.泥沙研究,2005（1）:46-52. HUANG S C,ZHANG S Y,YU J.Investigation on seabed evolution of Huibieyang zone in Hangzhou Bay［J］.Journal ofSediment Research,2005（1）:46-52.

［7］樂培九，朱玉德，程小兵，等.清水沖刷河床調整過程試驗研究［J］.水道港口,2007,28（1）:23-29. YUE P J,ZHU Y D,CHENG X B,et al.Study on the process of river-bed adjustment caused by clear water scouring［J］.Journal of Waterway and Harbor,2007,28（1）:23-29.

［8］楊華,趙沖久,侯志強,等.粉沙質海岸波浪和潮流作用下泥沙垂線分布研究［J］.水道港口,2006,27（3）:142-146. YANG H,ZHAO C J,HOU ZQ,etal.Study on sediment vertical distribution under the wave and current of silt coast［J］.Journal of Waterway and Harbor,2006,27（3）:142-146.

［9］吳以喜,麥苗,趙洪波.沖淤演變分析研究及其在海岸工程中應用［C］//中國海洋工程學會.第十四屆中國海洋（岸）工程學術討論會論文集：上冊.北京：海洋出版社，2009：50-54.

［10］李孟國,蔣厚武,吳以喜.溫州淺灘圍涂工程對港口航道的影響研究［J］.港工技術,2006（1）:1-4. LIM G,JIANG H W,WU Y X.Study on the influence of the reclamation of the Wenzhou Shoal on Wenzhou Harbor and its navigational channel［J］.Port Engineering Technology，2006（1）:1-4.

［11］李孟國,趙洪波,鄭敬云.甌江口航道治理工程水沙環境問題研究［J］.水利發電學報,2012,31（3）:107-113,118. LIM G,ZHAO H B,ZHENG JY.Study on tidal current and sediment environment problems in entrance channel regulation project in the Oujiang Estuary［J］.Journal of Hydroelectric Engineering,2012,31（3）:107-113,118.

［12］葉建國,莊曉榮.樂清灣港池開挖后驟淤可能性分析［J］.海洋工程,2007,25（1）:114-119. YE JG,ZHUANG X R.Analysis on possibility of silting up after harbor basin excavation in Yue-qing Bay［J］.The Ocean Engineering,2007,25（1）:114-119.

［13］曾長新,葉建國.樂清灣港區淺水深用可行性分析［J］.中國水運,2007（2）:44-45. ZENG C X,YE JG.Feasibility analysis of utilization of shallow water for deepwater port purposes in Yueqing Bay［J］.China Water Transport,2007（2）:44-45.

［14］李孟國,李文丹.甌江南口封堵方案研究［J］.中國港灣建設,2009（3）:19-22. LIM G,LIW D.Study on schemes of stopping flow in South Branch of Oujiang Estuary［J］.China Harbour Engineering，2009（3）:19-22.

［15］韓志遠,李文丹,龐啟秀.浙江麗水36-1終端至交付點輸氣管線工程海底管道沖刷研究［R］.天津：交通部天津水運工程科學研究所,2012.

［16］李孟國,鄭敬云.中國海域潮汐預報軟件Chinatide的應用［J］.水道港口,2007,28（1）:65-68. LIM G,ZHENG JY.Introduction to Chinatide software for tide prediction in China seas［J］.Journal of Waterway and Harbor, 2007,28（1）:65-68.

Study on sea bed erosion of gas transmission line from terminal to delivery point in Lishui 36-1

SUN Guo-m in1,HAN Zhi-yuan2,LIW en-dan2,LANG Yi-m ing1,LIU Zhi-gang1,PANG Qi-xiu2
(1.Offshore Oil Engineering Co.,Ltd.,Tianjin 529228,China; 2.Tianjin Research Institute for Water Transport Engineering,Key Laboratory of Engineering Sediment,Ministry of Transport,Tianjin 300456,China)

Based on a large number of measured data,seabed stability of project area especially in position of gas transmission line and hydrodynamic and sediment environment were studied in this paper.Then,a mathematical model of tidal flow and sediment was set up to calculate the siltation and erosion of seabed,and the distribution of siltation and erosion was given out.

gas transmission line;seabed erosion;numerical simulation;tidal flow;sediment;Lishui

TV 142；O 242.1

A

1005-8443（2013）04-0317-06

2013-01-07；

2013-03-28

孫國民（1966-），男，天津市人，高級工程師，主要從事海底管道設計研究工作。

Biography：SUN Guo-min(1966-),male,senior engineer.