鄭淑倩,楊青萍,劉犀力,駱旭佳
(1.浙江華東測繪地理信息有限公司,浙江 杭州 310014)
多波束測深系統(tǒng)是由單頻測深系統(tǒng)發(fā)展而來的,一種具有高效率、高精度和高分辨率的水下地形測量新技術(shù)[1]。單波束測深過程采用單點連續(xù)測量方法,測線上的測深數(shù)據(jù)十分密集,而在兩測線之間沒有數(shù)據(jù),是一種非覆蓋型測量[2],從而導(dǎo)致最終采集的數(shù)據(jù)不具代表性以及漏采集部分特征點。多波束則將測深技術(shù)從原先的點、線狀推展到面狀,測量范圍大、速度快且精度高。多波束測深一次的數(shù)據(jù)采集量較大,通常可達(dá)百萬級,無法直接用于構(gòu)建地形,因此如何做好多波束點云數(shù)據(jù)預(yù)處理工作是關(guān)鍵。目前,國內(nèi)關(guān)于該項技術(shù)的研究主要集中在杭州海洋二所、大連艦艇學(xué)院以及武漢大學(xué)測繪學(xué)院[3];但就多波束測深數(shù)據(jù)的加工、集中管理方面,還沒有一套現(xiàn)成的具有商品化水平的多波束數(shù)據(jù)后處理成圖軟件和數(shù)據(jù)管理軟件能為國內(nèi)用戶提供滿意的解決方案[4]。
多波束數(shù)據(jù)處理技術(shù)的滯后直接阻礙了其在水下工程應(yīng)用中的推廣。目前,沖刷淤積監(jiān)測評估大多仍采用斷面計算法和等高線容積法,前者以每隔一定距離的斷面特征曲線為基礎(chǔ)數(shù)據(jù),通過比較不同時期的曲線變化確定各斷面位置沖刷淤積面積,將相鄰斷面構(gòu)成的區(qū)域視為梯形體,再分別對各區(qū)域沖刷體積、淤積體積求和作為整個研究范圍的沖刷淤積量;后者則按不同高程面將水體細(xì)分成n層梯形體,計算得到每層體積,再將特征水位對應(yīng)高程以下的每層體積求和作為各特征水位對應(yīng)的庫容,然后通過前后兩次庫容求差得到對應(yīng)的沖刷淤積量,若為正值則認(rèn)為泥沙淤積,反之則泥沙被沖刷。楊立晉[5]等運用等高線容積法計算了黃羊泉水庫2007~2010年的泥沙淤積量,結(jié)果準(zhǔn)確度有所提高,具有較強(qiáng)的參考意義。但上述研究方法中均存在一定弊端,斷面計算法較適用于拐彎較少、叉河較少的線形區(qū)域,結(jié)果精度直接取決于斷面布設(shè)的密度;等高線容積法只能獲取整個庫區(qū)的沖刷量或淤積量,無法分別估算,且當(dāng)測量的等高線不閉合時,計算工作量巨大,難以在短時間內(nèi)完成。事實上,若能采用有效的手段對多波束采集的數(shù)據(jù)進(jìn)行預(yù)處理,保證特征點密度合理、位置分布均勻,以此構(gòu)建的水下地形既直觀又具有較高的精度。同時,利用GIS強(qiáng)大的空間分析和空間數(shù)據(jù)庫功能,提供了其他研究手段難以解決的水下地形沖淤演變、空間分布和計算結(jié)果可視化等工具[6],為研究水下地形的變化規(guī)律以及做好水利設(shè)施的安全維護(hù)提供了重要技術(shù)保證。
觀音巖水電站位于云南省麗江市華坪縣與四川省攀枝花市交界的金沙江中游河段,為金沙江中游河段規(guī)劃的8個梯級電站的最末一個梯級,上游與魯?shù)乩娬鞠嚆暯印k娬緣沃肺挥谔翂魏雍涌谏嫌卧颇?、四川兩省交界的河段上,河段長約3 km。壩址內(nèi),河道呈弓背型,向西凸出,河流由北西向流入后轉(zhuǎn)為北東向流出(圖1)。
本文以2013年、2014年壩址區(qū)多波束水下地形掃測結(jié)果為實驗數(shù)據(jù),采用丹麥Reason公司生產(chǎn)的SeaBat 7125 SV2 全功能高效版多波束測深儀,其主要性能參數(shù)見表1[7]。
圖1 研究區(qū)概況
表1 SeaBat性能參數(shù)說明表
預(yù)處理是對所有多波束測量數(shù)據(jù)進(jìn)行的初步整理,主要包括剔除異常點、稀釋冗余點和平滑等步驟,其目的是在保證數(shù)據(jù)精度的前提下,最大程度地稀釋點云密度,使最終構(gòu)建的DEM能真實反映水下地形空間變化特征(圖2)。
1)格網(wǎng)劃分。將實驗區(qū)按照指定的格網(wǎng)密度進(jìn)行劃分,判斷各點所在格網(wǎng)的索引號。在此基礎(chǔ)上,以格網(wǎng)內(nèi)所有點為操作對象,采用多線程并行技術(shù)遍歷所有格網(wǎng)執(zhí)行后續(xù)操作。一方面,各格網(wǎng)內(nèi)的點數(shù)據(jù)量是有限的,一般計算機(jī)硬件環(huán)境均可支持程序的正常運行;另一方面,多線程并行技術(shù)可使多個格網(wǎng)同時執(zhí)行相應(yīng)的運算流程,互不干擾,大大提高了數(shù)據(jù)處理效率。
2)異常點剔除。大數(shù)據(jù)集通常滿足正態(tài)分布規(guī)則,離期望值μ越近,異??赡苄栽叫?;反之,則異??赡苄栽酱蟆R虼?,程序?qū)⒙湓谥眯艆^(qū)間以外的點視為異常點并予以剔除。
3)點云密度稀釋。點密度過大將直接導(dǎo)致DEM無法生成或程序運行時間過長,生成的結(jié)果不夠連續(xù),視覺效果差。如何在滿足精度要求的前提下,最大限度地稀釋點云,本文采用(X,Y,Z)三維坐標(biāo)來標(biāo)識各點,并依據(jù)集聚統(tǒng)計分析方法建立決策對所有點進(jìn)行分類,按其相似性將所有點的空間分布分割或合并成一群。對于落在同一群的點可按一定比例進(jìn)行隨機(jī)抽樣,最終保留的點空間分布均勻且代表性強(qiáng),可作為構(gòu)建水下DEM的基礎(chǔ)數(shù)據(jù)。
經(jīng)測試,該算法可在32 s內(nèi)完成400萬點云數(shù)據(jù)的預(yù)處理工作。
圖2 多波束數(shù)據(jù)預(yù)處理流程圖
壩前淤積與壩后沖刷監(jiān)測管理信息系統(tǒng)是一款集成了多波束數(shù)據(jù)成圖,水下三維場景構(gòu)建與展示,庫容、沖刷淤積等重要信息自動分析等多種功能的綜合管理軟件。系統(tǒng)分為4個功能區(qū)(圖3):
1)數(shù)據(jù)管理模塊。該模塊主要包括多波束數(shù)據(jù)入庫和水下地形成圖。經(jīng)過預(yù)處理的多波束數(shù)據(jù)可由系統(tǒng)自動導(dǎo)入SQL Server數(shù)據(jù)庫,并按照“一期一表”的原則進(jìn)行管理。水下地形成圖包括調(diào)用數(shù)據(jù)表、生成點shp文件、構(gòu)建TIN以及輸出DEM等步驟。系統(tǒng)繼承了“一鍵式”的概念,將上述步驟整合成一個功能,大大簡化了操作流程。
2)基本操作模塊。該模塊主要負(fù)責(zé)圖形數(shù)據(jù)的二三維瀏覽和屬性查詢等。
3)二維分析模塊。該模塊提供對各類數(shù)據(jù)的空間分析操作,以獲取可靠有用的專題信息,包括緩沖區(qū)分析、疊加分析等。其中,剖面分析是最具特色的,程序支持任何斷面位置高程、坡度等信息的快速獲取,并可對同一位置不同時刻的高程(坡度)變化曲線進(jìn)行比對,以便掌握水下地形的變化過程。
4)三維展示與分析模塊。DEM數(shù)據(jù)可依據(jù)高程因子,配以適當(dāng)?shù)目鋸埾禂?shù),達(dá)到最佳的三維顯示效果。系統(tǒng)支持水下三維地形沿指定線路漫游,用戶可全方位多角度地瀏覽其地理變化特征。此外,該模塊提供了強(qiáng)大的沖淤分析功能,不僅可快速準(zhǔn)確地獲取沖刷淤積量,更能直觀地瀏覽沖刷區(qū)、淤積區(qū)的空間分布情況。
圖3 系統(tǒng)功能模塊
系統(tǒng)已投入觀音巖水電站試用,本文以水電站樞紐區(qū)2013年、2014年水下多波束點云為數(shù)據(jù)源,完成了水下地形建模(圖4)。經(jīng)驗證,模型精度達(dá)到三維仿真要求,并可用于沖刷淤積、庫容等重要信息的自動提取。
圖4 水下地形與大壩疊加的三維顯示效果
本文采用地形計算法分析了水電站樞紐區(qū)2013~2014年水下沖刷淤積變化情況(沖刷淤積量=后一次生成的DEM-前一次生成的DEM),若單元格為正值,則表明泥沙淤積,反之為泥沙被沖刷。從統(tǒng)計結(jié)果來看,研究區(qū)以泥沙淤積為主,淤積量達(dá)96.25 萬 m3,占總面積的12.49%,且集中分布在地勢較低的庫區(qū)底部(圖5);沖刷量則較少,僅有5.06 萬 m3,占總面積的2.88%。
圖5 沖刷淤積區(qū)空間分布圖
庫容曲線是表示水位與其相應(yīng)庫容關(guān)系的曲線,以水位為縱坐標(biāo),庫容為橫坐標(biāo)繪制而成,是水庫規(guī)劃設(shè)計和管理調(diào)度的重要依據(jù)。目前,在調(diào)洪演算時發(fā)現(xiàn)我國絕大部分水庫存在水量不平衡問題,其主要原因為庫容曲線精度不高,而水下DEM的精度直接影響了庫容統(tǒng)計結(jié)果的可靠性。本文以預(yù)處理后多波束點云構(gòu)建的DEM為基礎(chǔ)數(shù)據(jù),通過調(diào)用ArcEngine提供的庫容計算接口,分別統(tǒng)計了觀音巖水電站樞紐區(qū)2013年、2014年投影面積以及庫容隨水位的變化情況(圖6、表2)。
圖6 水位—面積、水位—庫容關(guān)系曲線圖
從統(tǒng)計結(jié)果來看,低水位區(qū)域庫容隨水位的變化較緩慢,而高水位區(qū)域庫容隨水位的增長呈直線上升趨勢。從時間尺度來看,2014年的庫容與投影面積比2013年均有所減少,證明了研究區(qū)在2013~2014年以泥沙淤積為主要現(xiàn)象這一論點。
表2 投影面積、庫容隨水位變化表(部分結(jié)果)
本文提出了一套切實可行的多波束數(shù)據(jù)預(yù)處理與應(yīng)用方法,在稀釋點密度時綜合考慮了平面位置和高程信息,使得最終保留的點空間分布均勻,以保證構(gòu)建的DEM能最真實地還原水下地形特征。該方法在觀音巖水電站樞紐區(qū)進(jìn)行了試用,并分析了2013~2014 年樞紐區(qū)水下地形變化過程。從運行結(jié)果來看,自動提取的沖刷淤積、庫容等重要信息是可靠的,且該系統(tǒng)具有工作效率高、工作強(qiáng)度小等優(yōu)點。
[1] 陳堅雄.多波束水下地形測量系統(tǒng)及其在水利建設(shè)方面的應(yīng)用[J].人民珠江,2005(4):21-23
[2] 張偉.多波束測深系統(tǒng)在水下地形測量中的應(yīng)用研究[D].北京:中國地質(zhì)大學(xué),2009:4
[3] 劉雁春.海洋測深空間結(jié)構(gòu)及其數(shù)據(jù)處理[J].測繪學(xué)報,2001,30(2):186
[4] 宋玲玲.多波束測量數(shù)據(jù)預(yù)處理研究[D].南京:南京航空航天大學(xué),2007:4
[5] 楊立晉,李曉飛.水庫泥沙淤積測量方法及淤積庫容計算的分析[J].甘肅水利水電技術(shù),2011,47(6):22-24
[6] 白世彪,王建,閭國年,等.水下地形沖淤變化可視化計算方法[J].工程勘察,2006(11):54-56
[7] 湯小云,吳金才.RESON SeaBat 7125多波束測深系統(tǒng)在水下測量中的應(yīng)用[J].科技視界,2014(34):133-135