借助Surfer軟件進行水下地形演變監(jiān)測分析

余學淵

摘 要：長江干線航道水下地形復雜多變，沖刷和回淤幅度此起彼伏，本文簡要敘述利用水深測量數(shù)據(jù)對長江南京以下12.5米深水航道水下地形演變監(jiān)測分析的方法，及制定航道水深監(jiān)測周期的意義，著重介紹Surfer曲面測繪軟件的等高線圖示法、三維圖示法和斷面圖示法及其應用。

關鍵詞：河床演變；XYZ坐標；數(shù)據(jù)白化；礙航淤積速度

中圖分類號：U697 文獻標識碼：A 文章編號：1006—7973（2017）05-0029-02

1 航道演變監(jiān)測分析的意義

經(jīng)過多年自然演變和人工治理，長江南京以下12.5米深水航道上下游便捷條件逐漸趨于穩(wěn)定，多汊并存的河勢格局基本形成。但上游來水、來沙條件的不穩(wěn)定性和因河勢不同造成的沿程阻力多變決定沿程及底沙輸移的不均勻性，航道水下地形演變監(jiān)測、回淤規(guī)律研究顯得尤為重要。動態(tài)跟蹤掌握航道演變趨勢與規(guī)律，可以為航道維護決策提供科學依據(jù)；積累觀測資料和分析成果，能夠為航道建設提供技術支持；深入研究掌握航道演變規(guī)律，創(chuàng)新航道養(yǎng)護方式，能夠增強航道維護管理的科學性和系統(tǒng)性，同時不斷提高航道養(yǎng)護管理質(zhì)量和水平。

2 監(jiān)測分析的方法和應用

由于長江干線航道的來水來沙條件瞬息多變，河床周界條件因地而異，河床演變的形式及過程極其復雜，現(xiàn)階段要進行精確的定量計算，尚有不少困難，但可借助某些手段對河床演變進行定性分析或定量估算?，F(xiàn)階段常用的幾種分析途經(jīng)如下：①利用條件相似河段的資料進行類比分析；②運用模型試驗的基本理論，通過河工模型試驗，對河床演變進行預測；③運用泥沙運動基本規(guī)律及河床演變基本原理，對河床變形進行理論計算；④天然河道實測資料分析。上面四種方法中，天然河道實測資料分析方法，是最基本、最常用的方法。

根據(jù)長年實測的水深測量數(shù)據(jù)，借助制圖軟件和地理信息系統(tǒng)軟件的數(shù)據(jù)分析功能進行航道演變規(guī)律分析應用，而最直觀的分析方法就是采用圖示法，其中最常用的為等高線圖示法、三維圖示法和斷面圖示法。下面以Surfer 12為例，簡要分析如何應用水深點測量數(shù)據(jù)進行航道水下地形演變的分析。

2.1 采集原始測量數(shù)據(jù)

收集不同時期獲得的以XYZ坐標離散點形式的水深點數(shù)據(jù)，所謂的不同時期，可以是枯水季、洪水季、同年或者不同年等，根據(jù)分析目的不同自行決定。通常不采用格柵或者內(nèi)插的水深點，目的是保證當前位置水深值的準確性，也是航道水下地形演變分析可靠性的關鍵。如果采用的軟件不具備投影改正和基準轉換功能，則需要收集經(jīng)過投影改正和基準轉換處理后的水深測量數(shù)據(jù)，并且不同時期的數(shù)據(jù)宜為同一基準面和坐標系。

2.2 格網(wǎng)化處理并繪制成圖

在數(shù)據(jù)格網(wǎng)化的時候，格網(wǎng)化算法有12種之多。為了如實反應地貌特征，盡可能選擇克里金法（Kriging）和反距離加權插值法（Inverse Distance to a Power）。因為這兩種方法在插值點與取樣點重合時，插值點的值就是樣本點的值，而其它方法不能體現(xiàn)樣本點的真實值。當所有其他參數(shù)相等的條件下，采用反距離加權插值法，越靠近結點（計算出的規(guī)則點）的原始數(shù)據(jù)點，對計算該結點的Z值貢獻越大。把不同時期采集水深點測量數(shù)據(jù)，通過格網(wǎng)化處理產(chǎn)生不同時期的規(guī)則矩形網(wǎng)格的三維數(shù)組X，Y，Z文件數(shù)據(jù)（*.grd），并繪制等高線圖（Contour Map）。直觀的等高線圖需要進行詳細設置，比如主要等高線、輔助等高線，注記、色塊填充及坐標系等的設置。為提高工作效率，減少重復工作量，在利用不同期測量數(shù)據(jù)繪制圖形時，我們可以采用同一等深線配置文件（*.LVL文件）。在高線圖和線框圖中可以直觀分辨出異?；蛘咄蛔兊臄?shù)據(jù)，通過檢查原始測量數(shù)據(jù)，配合粘貼圖，使用網(wǎng)格結點編輯器變更網(wǎng)格文件（*.grd）內(nèi)網(wǎng)格結點的Z值，并重新生成網(wǎng)格文件和各類圖元。

2.3 邊界文件及數(shù)據(jù)白化處理

在實際航道工程測量工作中，我們采集的數(shù)據(jù)量通常是比較大，采集區(qū)域和我們分析所需的區(qū)域不一，這時就要使用軟件中的邊界和白化處理功能。Surfer 對由離散的XYZ 數(shù)據(jù)文件生成網(wǎng)格文件時，是根據(jù)原始X、Y的取值范圍和所選用的數(shù)學模型，自動生成一個矩形網(wǎng)格。根據(jù)航道形狀把航道范圍線外（線內(nèi)）的數(shù)據(jù)剔除，生成一個由航道邊線封閉多邊形定義的區(qū)域，按照surfer的.BLN文件格式存儲在一個或多個以.BLN 為擴展名的邊界文件中，通過邊界文件及白化處理使航道邊線封閉多邊形定義的區(qū)域外（區(qū)域內(nèi)）的區(qū)域被剔除形成空白。不同時期的測量數(shù)據(jù)處理都采用同一邊界文件可以大大減少重復工作量，提高工作效率，方便數(shù)據(jù)比較分析。以包含邊界文件的新網(wǎng)格文件（*.grd）重新生成各類圖，可以看到，等高線圖上航道邊線范圍外的等高線被消除掉；在線框圖上航道邊線范圍外成為平行于基點平面的區(qū)域。下圖為長江南京以下12.5米深水航道某水道2016 年8月份兩次測圖獲得的離散高程數(shù)據(jù)，插值生成規(guī)則的格網(wǎng)數(shù)據(jù)并經(jīng)邊界文件空白處理后繪制的等高線圖。圖1 為其中兩次測量的等高線圖。圖2 為水深數(shù)據(jù)生成的曲面圖和線框圖。從等高線圖上等高線的疏密和色塊的變化，以及三維圖的空間特征可以直觀地看出航道水下地形演變趨勢。

2.4 回淤沖刷量及礙航淤積速度的比較和分析

不管是回淤沖刷量比較還是礙航淤積速度的分析，通常是使用多組數(shù)據(jù)進行比較分析。不同測量數(shù)據(jù)文件生成具有相同X、Y坐標的多個文件，并按比較分析的順序要求，利用Surfer軟件“數(shù)學”功能，相互之間以C=B-A生成一個新的輸出網(wǎng)格文件，其中A和B是輸入的網(wǎng)格文件，C是輸出網(wǎng)格文件。由于具有相同X、Y坐標的網(wǎng)格結點，這樣就從輸入網(wǎng)格文件A和B的相應結點上計算出Z的差值（Z`），計算的結果（X，Y，Z`）將保存到具有相同坐標的網(wǎng)格文件C內(nèi)。利用B-A函數(shù)得出的C1、C2、C3……，計算其體積和表面積，可以得到回淤沖刷量，結合兩次數(shù)據(jù)采集周期可以計算出礙航淤積速度。Surfer 軟件還可通過“殘差”功能計算出任一坐標位置的高程與設計航道底標高之間的差值。通過剖面生成功能在等值圖上按一定方向和距離切畫出航道剖面線，剖面生成的航道剖面數(shù)據(jù)文件能被多種二維繪圖軟件用來生成航道剖面圖。圖3為長江南京以下12.5米深水航道洪水季某一水道某一斷面相鄰四個月測量結果疊加的比較。若采用固定斷面測量的數(shù)據(jù)進行比較分析可得出更精確的結果，圖4為長江南京以下12.5米深水航道某一水道2016年8月到12月的測量數(shù)據(jù)分析計算的結果，清晰呈現(xiàn)出每個時期礙航淤積速度的變化。

2.5 航道水下地形變化的演示

為了更加直觀地反映航道水下地形變化，用Surfer 建立各個時期航道水下地形圖件，特別是三維立體顯示與渲染圖，并通過對圖形的迭合和圖形的堆壘，以及對繪制好的圖進行圖形整飾（例如修改、標注、渲染等），作為幻燈片制成演示文稿，可以方便、清晰、直觀地進行各個時期航道水下地形變化的演示（如圖5），結合比色刻度尺，圖上整個區(qū)域水深變化一目了然。

3 結束語

長江河道呈現(xiàn)多種形態(tài)分布，河床演變錯綜復雜，不同河道灘槽格局穩(wěn)定性差異較大。長江河床演變除了受徑流和潮流作用以外，還有其它一些因素的影響，如橋渡、港口碼頭、取水工程、護岸工程等等。諸多因素的分析后，再由此及彼、由表及里采用多種途徑方法結合起來進行綜合分析，探明河床演變的基本規(guī)律及主要影響因素，預估河床演變的發(fā)展趨勢，為制定合理的長江航道維護和整治工程方案提供科學依據(jù)。

參考文獻：

[1]GB 12327-1998. 海道測量規(guī)范[S]. 北京：中國標準出版社，1999.

[2]胡友元，黃杏元. 計算機地圖制圖[M]. 北京：測繪出版社，1991.

[3]鄔倫. 地理信息系統(tǒng)—原理、方法和應用[M]. 北京：科學出版社，2001.

[4]郭文偉，張鐵軍，袁世中等. 水運工程測量手冊[M].北京：人民交通出版社，2001.