秦美茵，孫秀峰，張飛龍,鐘小陽

(1.廣東省水利水電科學(xué)研究院，廣東 廣州 510635;(2.中山大學(xué) 土木工程學(xué)院，廣東 珠海 519000)

1 概述

廣東省河流眾多，以珠江流域及獨流入海的韓江流域和粵東沿海、粵西沿海諸河為主，集水面積占全省面積的99.8%，其余屬于長江流域的鄱陽湖和洞庭湖水系。全省流域面積在100 km2以上的各級干支流共614條。由于歷史原因，許多地區(qū)河湖管理范圍邊界不清，侵占河湖、破壞河湖問題時有發(fā)生，嚴(yán)重影響河湖生態(tài)空間管控。依法劃定河湖管理范圍，明確管理邊界線，是加強河湖管理的基礎(chǔ)性工作，也是《水法》《防洪法》《河道管理條例》等法律法規(guī)作出的規(guī)定，更是中央全面推行河長制湖長制明確的任務(wù)要求。據(jù)測算，肇慶市德慶縣50 km2以上河道總長為359 km，其中無堤河段總長為238 km，占比66%，在河道管理范圍劃定工作中，無堤河段易受地形、岸坡植被、建筑物等影響，是劃界工作的重難點。

針對河道管理范圍劃定任務(wù)，國內(nèi)學(xué)者開展了大量的研究工作。在技術(shù)應(yīng)用方面，歐陽顯良等[1]探索無人機航測技術(shù)在中小河流劃界中的應(yīng)用，認(rèn)為該技術(shù)可獲得測區(qū)立體像對、高分辨率三維點云數(shù)據(jù)和正射影像圖，可較清晰地識別堤腳線等河道管理范圍的起算線；夏祖?zhèn)サ萚2]結(jié)合正射影像光譜信息豐富、地面分辨率高的特點，提出一種融合激光雷達點云和正射影像的管理與保護范圍線提取方法；殷丹[3]提出將ArcGIS應(yīng)用于有堤河段和有規(guī)劃設(shè)計資料無堤河段管理范圍劃界工作中；李青林等[4]提出應(yīng)用奧維互動地圖提高河湖和水利工程管理范圍劃界工作效率。但以上技術(shù)應(yīng)用研究多集中于有堤河段管理范圍劃定，而山區(qū)中小河流大部分為無堤河段，針對無規(guī)劃設(shè)計資料的無堤河段，曹均昌[5]以紅柳河和新橋水庫為例就劃界過程中洪水分析計算進行探討；趙巨偉[6]針對無堤河段劃界，提出基于洪水分析計算的技術(shù)流程，并結(jié)合遼寧省丹東市草河劃界案例，剖析了無堤段劃界的實施過程和技術(shù)要點；羅姍姍等[7]以貴州省某河道劃界案例，在開展測量工作的基礎(chǔ)上分析了無堤河段劃界技術(shù)流程；以上無堤防河段管理范圍劃定研究多著重于洪水分析部分，鮮有將上述技術(shù)結(jié)合應(yīng)用，劃界成果精度及生產(chǎn)效率有待提高，本文以德慶縣悅城河為例，結(jié)合無人機航測、無人船斷面測量技術(shù)，剖析無堤防河段管理范圍劃定的實施流程和技術(shù)方法，提出無堤河段設(shè)計水位與岸坡交線自動化快速生成技術(shù)方案，以供技術(shù)人員交流和參考。

2 流域概況

悅誠河屬西江一級支流，發(fā)源于德慶禾叉頂，流經(jīng)古有、莫村、永豐、播植、悅城等鄉(xiāng)鎮(zhèn)，于悅城鎮(zhèn)注入西江。悅城河全長為82 km，集雨面積為877 km2，其中有43.4 km2由廣寧縣流入，9.97 km由高要市流入，總落差為798 m，平均坡降為1.62‰，多年平均流量為1.1 m3/s；流經(jīng)德慶縣河長為66.61 km，河口位于德慶縣悅城鎮(zhèn)。流域水系見圖1所示。

圖1 悅城河流域水系示意

3 劃界技術(shù)流程

3.1 資料收集與分析

對劃界河道開展調(diào)研工作，收集河道相關(guān)規(guī)劃設(shè)計資料、水利工程設(shè)計報告、河道所在流域水文站系列資料等，并對所采用的資料進行合理性分析。

3.2 無人機航測

對劃界河段進行無人機航測，航測作業(yè)開展前，首先通過現(xiàn)場調(diào)研及資料收集情況確定航測范圍；然后利用TY-3垂起型航測無人機搭載LIDAR航拍系統(tǒng)，根據(jù)預(yù)先確定的航攝范圍，執(zhí)行航測外業(yè)作業(yè)和像控布設(shè)；通過內(nèi)業(yè)數(shù)據(jù)處理計算，獲得劃界河段三維點云數(shù)據(jù)、1：2 000數(shù)字高程模型、1：2 000正射影像數(shù)據(jù)等工作底圖。

3.3 無人船河道斷面測量

采用基于網(wǎng)絡(luò)RTK技術(shù)的無人船測量系統(tǒng)進行水下地形測量，由千尋CORS、船載GNSS和便攜式電腦組成導(dǎo)航定位系統(tǒng)進行定位導(dǎo)航，采用船載測深儀自動采集水深，后借助專業(yè)數(shù)據(jù)處理軟件和繪圖工具進行數(shù)據(jù)內(nèi)業(yè)處理，獲取1：500河道斷面數(shù)據(jù)。

3.4 設(shè)計洪水分析計算

基于河道斷面測量成果，采用水文比擬法、綜合單位線法和新推理公式法3種方法計算設(shè)計洪水?，F(xiàn)有設(shè)計資料可靠的河段，直接采用現(xiàn)有洪峰流量成果，或采用水文比擬法計算設(shè)計洪水；無可靠設(shè)計資料的河段采用綜合單位線法和新推理公式法分別進行計算，并協(xié)調(diào)兩種方法算得的洪峰流量值相差不超過20%，然后選用綜合單位線法計算成果[8-9]。

3.5 確定管理區(qū)邊界

利用已有設(shè)計水位成果或建立河道一維數(shù)學(xué)模型推算劃界河段設(shè)計水位；利用Arcgis生成設(shè)計水位與岸坡交線并結(jié)合人工目視解譯獲得無堤河段管理范圍基準(zhǔn)線；然后根據(jù)劃界標(biāo)準(zhǔn)，分河段緩沖分析得到管理區(qū)邊界。

4 技術(shù)和方法

4.1 無人機航測

1)外業(yè)航拍

外業(yè)航拍需對攝區(qū)的地點及設(shè)施進行現(xiàn)場勘察和調(diào)研，選擇攝區(qū)最有利的氣象條件，并要盡可能的避免或減少地表植被和其他覆蓋物對攝影和測圖的不良影響，確保航攝像片能夠真實地顯現(xiàn)地面細(xì)部。在避開航攝范圍內(nèi)高壓電力線和軍民航空器的前提下，保證航攝飛行路線的直線性，并把項目區(qū)分成若干測段，每一測段再分為若干航帶。在航攝飛行中盡量保持飛機姿態(tài)的平穩(wěn)，轉(zhuǎn)彎半徑要大，飛機傾斜角不得大于15°，以防止GPS信號失鎖。

2)像控測量

在完成一定區(qū)域范圍的航攝，解算出影像數(shù)據(jù)后，即開始選取控制點、刺點的工作。采用千尋網(wǎng)連續(xù)運行參考站系統(tǒng)(CORS)作為測區(qū)控制系統(tǒng)參考，利用1個或若干個固定的、連續(xù)運行的GPS參考站，實時自動地提供經(jīng)過檢驗的GPS觀測值、各種改正數(shù)、狀態(tài)信息。通過RTK獲得測區(qū)的平面坐標(biāo)信息，即國家2 000坐標(biāo)系統(tǒng)。采用RTK獲得大地高程數(shù)據(jù)，通過獲得的RTK數(shù)據(jù)經(jīng)過廣東省地理信息中心的后期結(jié)算，獲得測區(qū)準(zhǔn)確的國家85高程系統(tǒng)下的像控點數(shù)據(jù)。

3)內(nèi)業(yè)處理

對外業(yè)航拍及像控測量數(shù)據(jù)進行內(nèi)業(yè)計算處理，獲得劃界河段三維點云數(shù)據(jù)、1：2 000數(shù)字高程模型(DEM)、1：2 000正射影像(DOM)數(shù)據(jù)等工作底圖，悅城河某河段DEM及DOM成果示意分別見圖2～3。

圖2 悅城河某河段DEM成果示意

圖3 悅城河某河段DOM成果示意

4.2 無人船河道斷面測量

1)導(dǎo)航定位

啟動測量軟件，按軟件提示設(shè)定記錄間隔、船速檔位、通訊ip、船只吃水等，由導(dǎo)航人員引導(dǎo)測量船至測區(qū)后開始作業(yè)，計算機實時采集定位、水深等數(shù)據(jù)，顯示到圖形界面。測量采用手動自動相結(jié)合的測量模式，測量數(shù)據(jù)的記錄間隔按1 s/點，實現(xiàn)導(dǎo)航、定位、數(shù)據(jù)采集自動化。作業(yè)過程中測量人員實時觀察RTK流動站的衛(wèi)星信號鎖定、固定解情況。

2)外業(yè)測量

水深外業(yè)測量采用船載測深儀自動采集，測量精度0.1%h。無人船的測量數(shù)據(jù)由無人船傳輸天線通過網(wǎng)橋傳輸?shù)降孛婀ぷ髡局校瑹o人船在工作中可以通過地面工作站和遙控器進行控制。測深過程中，測量人員通過電臺實時查看測深數(shù)據(jù)，及時調(diào)整航速，使測深數(shù)據(jù)保持穩(wěn)定、可靠[10]。

3)內(nèi)業(yè)處理

數(shù)據(jù)內(nèi)業(yè)處理分為數(shù)據(jù)預(yù)處理、數(shù)據(jù)改正和數(shù)據(jù)檢查，主要是借助專業(yè)數(shù)據(jù)處理軟件和繪圖工具進行數(shù)據(jù)處理。其中數(shù)據(jù)預(yù)處理主要包括定位中心與測深中心的偏心改正、水深點采集時間延時檢測、對定位數(shù)據(jù)進行數(shù)據(jù)異常值檢測及修正以及水深數(shù)據(jù)異常值檢測及修正；數(shù)據(jù)改正主要包括水深數(shù)據(jù)的吃水改正、聲速改正和水位改正；數(shù)據(jù)檢查主要包括制作圖件，并采用等深線圖示法及三維顯示法對圖件中存在的殘留異常值進行修正，最后輸出成果、制作成果圖件，斷面測量成果見圖4。

圖4 無人船河道斷面測量成果示意

4.3 設(shè)計洪水分析計算

現(xiàn)有設(shè)計資料可靠的河段，直接采用現(xiàn)有洪峰流量成果，或采用水文比擬法計算設(shè)計洪水；無設(shè)計資料的河段，采用綜合單位線法和推理公式法分別計算劃界河段各頻率設(shè)計洪水成果，并協(xié)調(diào)兩種方法算得的洪峰流量值相差不超過20%，然后選用綜合單位線法計算成果，悅城河設(shè)計洪水計算結(jié)果見表1。

表1 悅城河設(shè)計洪水計算結(jié)果

4.4 設(shè)計水位計算

獲得設(shè)計洪水計算結(jié)果后，建立河道一維數(shù)學(xué)模型推算悅城河設(shè)計水位[11-12]，計算結(jié)果見表2。

表2 悅城河設(shè)計水位計算結(jié)果

4.5 設(shè)計水位與岸坡交線自動化快速生成技術(shù)

利用Arcgis將設(shè)計水位計算結(jié)果與無人機航測獲取的高分辨率DEM數(shù)據(jù)進行疊加分析，獲取劃界河段設(shè)計水位與岸坡交線。具體過程分析如下：

1)獲取設(shè)計水位數(shù)據(jù)后，利用Arcgis對水位數(shù)據(jù)沿河道中心線進行二維插值，獲得二維水位數(shù)據(jù)。

2)將二維插值后的水位數(shù)據(jù)生成不規(guī)則三角網(wǎng)(TIN)數(shù)據(jù)集文件，具體操作步驟為：打開ArcToolbox工具→3DAnalyst工具→數(shù)據(jù)管理→TIN→創(chuàng)建TIN。

3)將不規(guī)則三角網(wǎng)(TIN)數(shù)據(jù)集文件轉(zhuǎn)換成柵格數(shù)據(jù)文件，具體操作步驟為：打開ArcToolbox工具→3DAnalyst工具→轉(zhuǎn)換→由TIN轉(zhuǎn)出→TIN轉(zhuǎn)柵格。完成柵格數(shù)據(jù)轉(zhuǎn)換，得到設(shè)計水位柵格圖層，設(shè)計水位柵格數(shù)據(jù)示意見圖5。

圖5 設(shè)計水位柵格數(shù)據(jù)示意

4)應(yīng)用設(shè)計水位柵格數(shù)據(jù)與相應(yīng)河道高分辨DEM數(shù)據(jù)相減，獲得水深數(shù)據(jù)，具體操作步驟為：打開ArcToolbox工具→3DAnalyst工具→柵格計算→減。完成設(shè)計水位柵格圖層與相應(yīng)河道高分辨DEM數(shù)據(jù)處理，得到水深數(shù)據(jù)文件，河道高分辨DEM數(shù)據(jù)示意見圖6，水深數(shù)據(jù)示意見圖7。

圖6 河道高分辨DEM數(shù)據(jù)示意

圖7 水深數(shù)據(jù)示意

5)對水深數(shù)據(jù)進行重分類，便可獲得設(shè)計水位與岸坡交線，具體操作步驟為：將水深圖層載入ArcMap→右鍵圖層屬性→符號系統(tǒng)→已分類→分類。對水深數(shù)據(jù)進行重分類，數(shù)值小于零表示地面高程高于設(shè)計水位的區(qū)域，數(shù)值大于零表示設(shè)計洪水淹沒范圍，設(shè)置區(qū)間-0.5～0作為緩沖區(qū)，此數(shù)值可根據(jù)數(shù)據(jù)分辨率設(shè)定。將無物理意義區(qū)間-247.38～-0.5設(shè)置為無顏色區(qū)域，緩沖區(qū)-0.5～0與設(shè)計洪水淹沒范圍區(qū)0～141.14分別設(shè)置不同顏色，便可獲得設(shè)計水位與岸坡清晰交線，同時緩沖區(qū)的設(shè)置亦可直觀表現(xiàn)地勢的陡緩。水位岸坡交線成果圖見圖8。

圖8 水深數(shù)據(jù)重分類后獲取水位岸坡交線成果示意

獲得設(shè)計水位與岸坡交線圖層后，基于高分辨率無人機航攝正射影像和現(xiàn)場調(diào)研情況結(jié)合人工目視解譯繪制設(shè)計水位與岸坡交線獲得無堤河段管理范圍基準(zhǔn)線；然后根據(jù)劃界標(biāo)準(zhǔn)，分河段緩沖區(qū)分析得到管理區(qū)邊界。根據(jù)《廣東省河湖管理范圍劃定技術(shù)指引(試行)》，無堤河段劃界分為無堤山區(qū)河段和無堤平原河段，無堤山區(qū)河段以設(shè)計水位與岸坡交線作為管理區(qū)邊界，無堤平原河段以設(shè)計水位與岸坡交線為管理范圍基準(zhǔn)線，根據(jù)劃界標(biāo)準(zhǔn)，結(jié)合護堤地寬度要求外延一定距離作為管理區(qū)邊界。

5 結(jié)語

本文通過無人機航測獲取劃界河段三維點云數(shù)據(jù)、高分辨率數(shù)字高程模型、正射影像數(shù)據(jù)等工作底圖，通過無人船河道斷面測量獲取1：500河道斷面數(shù)據(jù)，采用水文比擬法、綜合單位線法和新推理公式法3種方法分析計算設(shè)計洪水，然后建立河道一維數(shù)學(xué)模型推算劃界河段設(shè)計水位；最后利用Arcgis將設(shè)計水位計算結(jié)果與無人機航測獲取的高分辨率DEM數(shù)據(jù)進行疊加分析，實現(xiàn)設(shè)計水位與岸坡交線的自動化快速生成。

在運用無人機航測獲取的地形數(shù)據(jù)確定設(shè)計洪水位的過程中，受岸坡植被、建筑物等影響，存在生成的設(shè)計水位與岸坡交線與實際地形誤差，對于河岸建筑物較多的敏感河段應(yīng)結(jié)合現(xiàn)場調(diào)研，RTK定點測量等方式進行修正，提高劃界成果的精度。另外，應(yīng)用設(shè)計水位確定河道管理范圍較適用于山區(qū)河道或者地勢較高的平原河道，對于地勢較低的平原河道，計算的淹沒范圍普遍偏大，而中小河流一般河寬較窄，應(yīng)用設(shè)計水位與岸坡交線作為河道管理區(qū)邊界易導(dǎo)致管理范圍偏大。

文章提出無堤河段設(shè)計水位與岸坡交線自動化生成的技術(shù)方案，讀者后續(xù)通過Python語言或Arcgis自定義工具箱可簡化操作步驟，進一步提高生產(chǎn)效率。