采用水文分析法對林區道路涵洞選址的研究

朱磊，陳偉，鄭云峰，金久宏

(1.國家林業局華東林業調查規劃設計院，杭州310019；2.杭州市蕭山區林業局，杭州 310019)

采用水文分析法對林區道路涵洞選址的研究

朱磊1，陳偉1，鄭云峰1，金久宏2

(1.國家林業局華東林業調查規劃設計院，杭州310019；2.杭州市蕭山區林業局，杭州 310019)

為了對林區道路的涵洞進行科學選址，本文以杭州市蕭山區為研究區域，以DEM數據為研究基礎，采用GIS水文分析法對DEM進行填充洼地、計算流向、計算匯流累積量、閾值篩選、疊加求交等操作，得到涵洞選址大致位置，從而進行小范圍詳細勘察，得到合理的涵洞選址 。與傳統作業模式相比，該方法能充分依據GIS理論和計算機技術科學合理地選點，從而大大減少了外業勘查設計人員的工作量，不失為一種新的探索模式。

水文分析；林區道路；涵洞；DEM

0 引言

建設林區道路是促進森林旅游及山區經濟發展的有效途徑[1]。林區道路建設對當地生態環境破壞較大，建設不當會引起山洪、泥石流、山體滑坡及水土流失等地質災害，因此防水排水工程至關重要，如何為林區道路的涵洞科學選址成為研究熱點。

水文分析法基于數字高程模型(DEM)建立水系模型，對流域水文特征和模擬地表水文過程進行研究，并對未來的地表水文情況作出估計[2-3]。該方法廣泛應用于區域規劃、農林、災害預測、道路設計等行業和領域[4-7]。近年來，國內外學者對水文分析法及其應用做了大量的研究。Nishio M等[8]對DEM進行水文分析，模擬洪水水文過程，制定水文過程圖，從而精確劃分洪水區；Jaseim E S等[9]對DEM進行水文分析，提取伊拉克基爾庫克省內的水路信息，精度較高；馬勤學等[10]基于水文分析法提取水系網絡，并根據水系分級，合理設置水土保持防護措施；張國珍、安祺等[11-13]對10m分辨率的DEM進行水系提取，提取的結果與地形圖數據做對比，發現精度很高。Walker Jeffrey P、Kenward T[14-16]等研究表明，水文信息提取精度與DEM分辨率有很大關系，分辨率過低不能真實反映地形地貌特征，分辨率過高，易產生零碎斑塊，影響精度。林區道路涵洞選址的精度取決于水文信息的提取精度，而水文信息的精度與DEM的分辨率和處理方法有關。本文將水文分析法方法應用到林區道路涵洞選址中，通過林區道路與水文、水系交叉求交集得到涵洞大致位置，然后人工對重點區域進行詳細勘察設計，此方法可大大減少外業工作量，提高工作效率，不失為一種新的探索模式。

本文以5m分辨率DEM為基礎，首先，對DEM進行填洼操作、流向計算、流量計算，得到累積流量柵格圖；然后，通過設置閾值，進行柵格計算，得到符合精度要求的柵格水系圖；最后，把柵格數據轉換為矢量數據，與林區道路求交集，得到涵洞布置點。

1 數據及處理

1.1 研究區概況

浙江省杭州市蕭山區地處錢塘江南岸，為杭州市南大門，介于東經120°04′22″～120°43′46″，北緯29°50′54″～ 30°23′47″，土地總面積1 414.17 km2，全境東西寬約57.2 km，南北長約59.4 km(如圖1所示)。屬于浙東低山丘陵區北部至浙北平原南部的過渡地區，地勢南高北低，自西南向東北傾斜，中部略呈低洼。全境可分為4個不同特征的地貌區：低山區、丘陵區、山間河谷平原區、錢塘江河口平原區。

圖1 研究區域Fig.1 Study Area

蕭山區江河眾多，水網密布，水系發達，水域面積約241.4 km2，占全區土地總面積的17.7%。本區主要江、河、湖泊，按地理位置和河流特征，可分為北部、中部、南部3個自然區域，相應形成3個既自成體系又相互關聯的水系，統屬錢塘江水系。境內森林物種較多，因人類活動頻繁，原生植被基本已消失，局部保留天然次生林。主要森林植被類型有暖性針葉林、常綠闊葉林、落葉闊葉林、常綠落葉闊葉混交林、針闊混交林、竹林、經濟林和灌木林等。

1.2 數據來源與處理

研究以1∶2 000基礎地形圖為數據源，坐標系為國家2 000大地坐標，等高距為2 mm。該數據精度較高，能精細表達地形地貌。采用Arcgis對地形圖處理流程如下：①提取等高線，對等高線進行范圍裁剪；②對等高線進行空間插值操作，生成5 m分辨率DEM；③對DEM進行空間投影，投影方式為高斯克呂格，3°帶，40帶號，中央經線為120°。由此得到的5 m分辨率DEM。

2 選址方法

2.1 選址原則

(1)減少對生態環境的破壞和降低地質災害發生的風險。依據地形地貌的水文特征合理布設涵洞，做到充分利用自然地形。

(2)合理設置匯流累積量閾值，均勻布設涵洞。匯流累積量小的地方通過設置排水溝，引導匯集附近涵洞排水。

(3)林區道路建設中應處理好技術指標選用與生態保護的關系。道路設計時要充分考慮道路建設對森林、景觀、野生動植物、生物多樣性和其他自然資源保護的影響。

2.2 選址技術路線

涵洞選址首先是對DEM進行預處理，得到高精度基礎數據；其次合理設置閾值，并結合地形圖的水系數據進行對比修正；然后由水文分析法得到的水系圖與林區道路進行疊加求交集，得到涵洞初步選址點；最后再與林相圖疊加，綜合考慮，篩選出最優結果。如圖2所示。

圖2 技術流程圖Fig.2 Flow Chart of the Methods

(1)填充洼地：在【ArcToolbox】中，執行命令[SpatialAnalyst工具]→[水文分析]→[填洼]。DEM四周較高區域圍繞的洼地，影響水的流向，無法進行有效的水文分析。洼地產生主要有兩方面原因：一是在對矢量等高線進行空間插值的時候，產生局部誤差，這種誤差不可避免；二是山地林區存在山塘水池、采石場、泥水沖刷溝壑等真實地形地貌。如圖3所示，填洼前，有部分區域顯示 “椒鹽現象”，實為洼地，填洼后，該區域洼地消失。對DEM進行填洼后，得到質量較高的DEM。

圖3 填充洼地效果Fig.3 The effect of filling depression

(2)計算流向：在【ArcToolbox】中，執行命令[SpatialAnalyst工具]→[水文分析]→[流向]。對DEM進行填洼后，DEM表面比較平滑，滿足水流動力條件。計算水流方向是創建從每個像元到其最陡下坡相鄰點的流向的柵格單元，根據像素點與四周等距像素點的高程差來決定水流方向。水流流向的精度跟填洼質量密切相關。本文流向分為1(東)、2(東南)、4(南)、8(西南)、16(西)、32(西北)、64(北)、128(東北)八個方向，如圖4所示。

圖4 水流流向圖Fig.4 Flow direction

(3)計算匯流累積量：在【ArcToolbox】中，執行命令[SpatialAnalyst工具]→[水文分析]→[流量]。對流向柵格的每個像元累計流量進行統計，流量累積將基于流入輸出柵格中每個像元的像元數。高流量的輸出像元是集中流動區域，可用于標識河道。流量為零的輸出像元是局部地形高點，可用于識別山脊，如圖5所示，白色線條為匯流累積量較大區域。

圖5 流量分布Fig.5 The distribution map of flow accumulation

(4)提取水文特征：在【ArcToolbox】中，執行命令[Spatial Analyst 工具]→[地圖代數]→[柵格計算器]，在[地圖代數表達式]中輸入公式：Con(Flow Accumulation1>600，1)。上步得到的匯流流量柵格數據是像元值連續的柵格數據，需要設置合理閾值，提取滿足精度要求的水文特征線。根據經驗閾值設置為600，與實際地形比較吻合。如圖6所示，細黑色為基礎地形圖水系，粗黑色部分為提取水文特征線，疊加分析，吻合較好。

(5)柵格水系矢量化：對上步得到的水文特征線進行矢量化后得到矢量水文特征線，由于柵格數據轉為矢量易產生鋸齒狀，為此，需要進行平滑處理。

(6)矢量水文特征線與林區道路求交集即可得到涵洞坐標點。

圖6 水文特征線與水系對比Fig.6 Comparison between hydrological characteristic line and river system

3 結果分析與討論

3.1 涵洞初始布設結果

由圖7可以看到涵洞布設在水文特征線與林區道路交叉點處，在此處設置涵洞可以解決以下幾個問題：①充分利用地形地貌特征，減少土方施工；②積水順流而下，自然通過涵洞，減少泥水沖刷，有效保護了生態環境；③借助于GIS空間分析方法，可大大減少施工勘查的工作量，有效輔助涵洞的科學布設。

圖7 涵洞選址結果Fig.7 Culver location results

3.2 需要注意的問題

(1)林區道路應根據沿線地面水和地下水等情況，設置必要的邊溝、截水溝、排水溝和盲溝等排水設施，與涵洞融會貫通。

(2)涵洞選址正確的關鍵是匯流累積量閾值的篩選，篩選不當會生成一些偽水系，所以建議選址后，進行抽樣驗證。

(3)本文主要是研究涵洞選址，但需要根據現地情況進行合理修正。比如有的地方水流量較大，溝壑跨度較寬，這時可能需要修建橋梁。

4 結束語

本文以杭州市蕭山區為研究區域，采用GIS水文分析法來對蕭山林區的道路涵洞進行選址。采用地形圖精度較高，生成DEM也比較精細。并對DEM進行填洼、計算流向、計算匯流累積量等步驟得到水文特征線，最終水文特征線與林區道路進行疊加求交，得到林區道路涵洞的坐標點。

(1)該方法能充分依據GIS理論和計算機技術科學合理地選點，對于提高涵洞勘查選址精確度有一定的輔助作用。

(2)通過科學計算得到較為精確的重點區域，幫助外業人員重點作業，大大減少野外工作強度和周期。

[1]白秀萍，陳紹志，何友均，等.國外林區道路發展現狀及啟示[J].世界林業研究，2015，28(1)：85-91.

[2]胡婕，李崇巍，崔鐵軍.基于DEM地表水文特征和匯水子流域的提取研究[J].環境科學與管理，2011，36(9)：131-135.

[3]王云，梁明，汪桂生.基于ArcGIS的流域水文特征分析[J].西安科技大學學報，2012，32(5)：581-585.

[4]李俊超，馬倩，陶鈞.基于ArcGIS的水文流域分析及應用[J].地理空間信息，2012，10(6)：121-123.

[5]Torterotot J B，Perrier C，Bergeron N E，et al.Influence of Forest Road Culverts and Waterfalls on the Fine-Scale Distribution of Brook Trout Genetic Diversity in a Boreal Watershed[J].Transactions of the American Fisheries Society，2014，1436(6)：1577-1591.

[7]Sadeghi S H.Analysis of Temporal Variation of Hydrological Processes in Experimental Small Plots in a Forest Watershed[J].2015.

[8]Nishio M，Mori M.Hydrologic Analysis of a Flood Based on a New Digital Elevation Model[J].ISPRS-International Archives of the Photogrammetry.Remote Sensing and Spatial Information Sciences，2015，XL-7/W4：127-134.

[9]Jaseim E S，Mohammed F G.Waterways Finding in the Province of Kirkuk-Iraq Based on Hydrological Analysis of Digital Elevation Model[J].International Journal of Scientific & Engineering Research，2015，6(1).

[10]馬勤學，李紅旭，于洋，等.ArcGIS水文分析模塊在干旱區煤田滅火區生態修復中的應用研究[J].新疆環境保護，2013，35(3)：33-37.

[11]張國珍，嚴恩萍，洪奕豐，等.基于DEM的東江湖風景區水文分析研究[J].中國農學通報，2013，29(2)：172-177.

[12]安祺，楊霏，陳丹，等.基于Hydrology的山區1：1萬DEM水系提取研究[J].重慶工商大學學報(自然科學版)，2012，29(11)：87-92.

[13]關穎慧，鄭粉莉，王彬，等.基于DEM的黑龍江賓州河流域水系提取試驗研究[J].水土保持通報，2012，32(1)：127-131.

[14]Saud M M A.Hydrological Analysis of DEMs[M].Flood Control Management for the City and Surroundings of Jeddah，Saudi Arabia.2015.

[15]Gregoretti C，Degetto M，Bernard M.Comparative Analysis of the Differences Between Using LiDAR and Contour-based DEMs for Hydrological Modeling of Runoff Generating Debris Flows in the Dolomites[J].Frontiers in Earth Science，2015，3：21.

[16]顧用紅，舒振文，張火青.DEM在流域水文特征分析中的應用[J].人民珠江，2001，22(4)：5-6.

Research on Hydrological Analysis for Locating Forest Road Culvert

Zhu Lei1，Chen Wei1，Zheng Yunfeng1，Jin Jiuhong2

(1.East China Forest Inventory and Planning Institute，State Forestry Administration，Hangzhou 310019；2.Xiaoshan District Forestry Bureau of Hangzhou Municipality，Hangzhou 311203)

In order to scientifically locating forest road culvert，this paper takes Xiaoshan District of Hangzhou as the study area，using GIS hydrological analysis for filling DEM depressions，calculating flow direction and accumulation，selecting threshold，overlapping and intersecting the maps，etc.to get the location of culverts.After a small range of detailed investigation，a reasonable location of culverts can be obtained.Compared with the traditional work，this method can select the locations scientifically and reasonably based on GIS and computer technology.It can greatly reduce the workload of field survey，which can be a new exploration method.

Hydrological analysis；forest road；culvert；DEM

2016-10-28

國家自然科學基金資助課題(61201177)；國家林業公益性行業科研專項(201504323)

朱磊，碩士，高級工程師。研究方向：智慧林業應用。

*通信作者：陳偉，博士研究生。研究方向：智慧林業與資源環境遙感。E-mail：lajiao.1225@163.com

朱磊，陳偉，鄭云峰,等. 水文分析法對林區道路涵洞選址研究[J].森林工程，2017，33(3)：85-88.

U 445.7

A

1001-005X(2017)03-0085-04