柳克松

（中煤科工集團(tuán)沈陽(yáng)設(shè)計(jì)研究院有限公司，遼寧 沈陽(yáng) 110000）

近些年礦區(qū)周邊土地利用狀況與礦區(qū)生態(tài)環(huán)境持續(xù)惡化，礦區(qū)生態(tài)保護(hù)迫在眉睫，礦區(qū)生態(tài)恢復(fù)與治理已成為當(dāng)前研究的熱點(diǎn)[1]。礦區(qū)的水土保持度是衡量礦區(qū)生態(tài)恢復(fù)情況的主要指標(biāo)之一[2]，土方量的精確計(jì)算是水土保持監(jiān)測(cè)的重要手段。在傳統(tǒng)的土方量計(jì)算方法中，散點(diǎn)法適用于地形起伏小且高程點(diǎn)分布較為均勻地區(qū)[3-4]；方格網(wǎng)法對(duì)高程的測(cè)量有著更高要求。隨著傾斜攝影測(cè)量的飛速發(fā)展，數(shù)字高程模型應(yīng)用在土方量計(jì)算中優(yōu)勢(shì)明顯[5-6]，但是其對(duì)于礦區(qū)低山丘陵區(qū)域的復(fù)雜地形下土方量計(jì)算未能提出完整的處理方法。針對(duì)具有低山丘陵區(qū)域特征的礦區(qū)水土保持監(jiān)測(cè)困難的問(wèn)題，基于無(wú)人機(jī)傾斜攝影測(cè)量—數(shù)字高程模型技術(shù)對(duì)礦區(qū)進(jìn)行土方量多次計(jì)算，得到礦區(qū)水土流失比例，引入土壤流失方程，準(zhǔn)確評(píng)估水土保持率，結(jié)合礦區(qū)地質(zhì)土壤分布特性，為礦區(qū)生態(tài)恢復(fù)評(píng)估提供技術(shù)支撐。

1 原理與方法

1.1 水土保持監(jiān)測(cè)基本原理

本文水土保持監(jiān)測(cè)的原理以土方量多次計(jì)算為基礎(chǔ)，結(jié)合水土流失模型完成水土保持研究。土方計(jì)算的基礎(chǔ)數(shù)據(jù)是傾斜攝影測(cè)量獲得的DEM數(shù)據(jù)，它以離散的數(shù)字進(jìn)行表達(dá)，關(guān)系式表示如下：

Z=f（x，y）i=1,2,3...,n（1）式中：x、y分別為離散變量，表示地表網(wǎng)格某點(diǎn)的坐標(biāo)；Z表示該點(diǎn)對(duì)應(yīng)的高程。土方量計(jì)算即為一個(gè)二重積分的應(yīng)用，基于柵格相減的方法，并將所有格網(wǎng)的單元計(jì)算的量進(jìn)行統(tǒng)計(jì)，實(shí)現(xiàn)低山丘陵區(qū)域的土方量計(jì)算。

為研究水土保持率，引入通用土壤流失方程A=R×K×LS。此方程未將植被覆蓋度與特定保護(hù)措施作為影響因子。A是模型預(yù)測(cè)的年土壤侵蝕量；R是降雨侵蝕力因子，即為降水產(chǎn)生的徑流對(duì)土壤造成侵蝕的動(dòng)力指標(biāo)。以月平均降雨量和年平均降雨量為基數(shù)，R可表示為如公式（2）所示[7]，其中Pi為月平均降雨量，P為年平均降雨量。

式中：LS為地形因子，L即坡長(zhǎng)因子，是指標(biāo)準(zhǔn)化到22.13 m 坡長(zhǎng)上的土壤侵蝕量；S即坡度因子，是指標(biāo)準(zhǔn)化到5.14°坡度下的土壤侵蝕量。在實(shí)際的研究中小尺度研究區(qū)直接利用實(shí)測(cè)地形數(shù)據(jù)，而大尺度研究區(qū)會(huì)利用DEM 數(shù)據(jù)提取信息來(lái)計(jì)算LS因子。本文研究區(qū)為大尺度區(qū)域，利用無(wú)人機(jī)傾斜攝影測(cè)量技術(shù)生成的DEM 作為計(jì)算地形因子的基礎(chǔ)數(shù)據(jù)。K是土壤可蝕性因子，其反映了土壤對(duì)侵蝕應(yīng)力分離和搬運(yùn)作用的敏感性。

1.2 技術(shù)路線

基于無(wú)人機(jī)航空傾斜攝影測(cè)量技術(shù)獲取基礎(chǔ)影像數(shù)據(jù)，采用GPS RTK 技術(shù)進(jìn)行野外像控點(diǎn)測(cè)量，運(yùn)用ContextCapture 軟件進(jìn)行空三加密等處理生成DEM 和DOM。進(jìn)而將獲取的DEM 導(dǎo)入ArcGIS 中，進(jìn)行原地貌高程統(tǒng)計(jì)，推演規(guī)劃設(shè)計(jì)地貌，計(jì)算其與原地貌高程差值，獲得礦區(qū)的土方量。收集降雨數(shù)據(jù)、計(jì)算LS因子與土壤可蝕性因子，進(jìn)而計(jì)算土壤侵蝕量，與DEM 土方量計(jì)算相互印證。技術(shù)路線圖如圖1。

2 實(shí)驗(yàn)與分析

2.1 實(shí)驗(yàn)區(qū)概況

礦區(qū)位于遼寧省本溪市東南郊16 km，南芬區(qū)北9 km 處思山嶺村北側(cè)一帶，礦區(qū)面積3.74 km2。地形以低山丘陵區(qū)域?yàn)橹?，礦區(qū)廢石場(chǎng)所在區(qū)域土層較薄，施工前地表土未剝離，未有其他相關(guān)水土保持措施布設(shè)。

2.2 基礎(chǔ)數(shù)據(jù)處理與土方計(jì)算

2.2.1 場(chǎng)地影像采集

本文通過(guò)對(duì)實(shí)測(cè)區(qū)域的實(shí)地勘察，結(jié)合《低空數(shù)字?jǐn)z影航空攝影測(cè)量外業(yè)規(guī)范》要求，對(duì)飛行區(qū)域與航線規(guī)劃如下：采用飛馬D2000 搭配DOP-3000 鏡頭對(duì)礦區(qū)進(jìn)行傾斜攝影測(cè)量，在飛行區(qū)域內(nèi)共規(guī)劃航線1 條，設(shè)置仿地飛行高度為120 m，地面采樣距離為3 cm，飛行路線的旁向重疊為70%，航向重疊為80%，本航線共獲取的相片數(shù)為27 000余張。本實(shí)驗(yàn)在實(shí)驗(yàn)區(qū)域內(nèi)均勻設(shè)置地面像控點(diǎn)121 個(gè)，每個(gè)像控點(diǎn)間距為10 m，在地形復(fù)雜區(qū)域處適當(dāng)加密像控點(diǎn)。研究區(qū)域如圖2。

圖 1 技術(shù)路線圖

圖2 研究區(qū)域

2.2.2 三維模型處理

（1）在Contextcapture 中導(dǎo)入無(wú)人機(jī)航測(cè)相片并進(jìn)行空中三角測(cè)量后，輸入地面相控坐標(biāo)數(shù)據(jù)結(jié)合傾斜相片進(jìn)行刺點(diǎn)工作，完成礦區(qū)絕對(duì)三維模型的生成；（2）數(shù)字高程模型生成與研究區(qū)域土方計(jì)算。在礦區(qū)三維模型中選取地面高程點(diǎn)，起算面高程300 m，填方3 208 332.678 2 m3，挖方168 039.007 3 m3，并由DEM 計(jì)算分析得到坡度分析如圖3。

圖3 坡度計(jì)算分析圖

針對(duì)潛在水土流失區(qū)域大多位于坡度較陡區(qū)域，因此結(jié)合坡度分析圖中坡度范圍在12°～25°、25°～35°、35°～90°內(nèi)水土保持度較差，故在三維模型中提取坡長(zhǎng)參數(shù)作為水土保持度分析因素之一。提取坡長(zhǎng)區(qū)域如圖4，坡長(zhǎng)見(jiàn)表1。

表1 坡長(zhǎng)統(tǒng)計(jì)表

圖 4 坡長(zhǎng)圖

2.3 水土保持度分析

項(xiàng)目區(qū)主要土壤類(lèi)型為棕壤土、草甸土、水稻土三種土類(lèi)，工程占地范圍內(nèi)土壤類(lèi)型以棕壤土為主，土層較薄，厚度在10～30 cm 不等。項(xiàng)目區(qū)處于長(zhǎng)白植物區(qū)系與華北植物區(qū)系的交匯區(qū)，地帶性植被為溫帶針闊葉混交林，林草覆蓋率為72.60%。土壤侵蝕類(lèi)型區(qū)為東北黑土區(qū)-長(zhǎng)白山山地丘陵區(qū)，土壤侵蝕類(lèi)型主要以水力侵蝕為主。結(jié)合上述坡長(zhǎng)、坡度數(shù)據(jù)及將年平均降雨量等因子換算到同一投影坐標(biāo)系下，像元大小為30 m×30 m，基于ArcGIS10.2 柵格運(yùn)算，根據(jù)水利部頒布的《土壤侵蝕分類(lèi)分級(jí)標(biāo)準(zhǔn)》可將水土流失分為6 個(gè)等級(jí)，即微度侵蝕為小于500 t/km2，輕度侵蝕為500～2500 t/km2，中 度 侵 蝕 為2500～5000 t/km2，強(qiáng) 度 侵 蝕5000～8000 t/km2， 極 強(qiáng) 度 侵 蝕 為8000～15 000 t/km2，劇烈侵蝕為大于等于15 000 t/km2。對(duì)其進(jìn)行柵格分類(lèi)，得到各個(gè)侵蝕等級(jí)像元數(shù)見(jiàn)表2。

表2 土壤侵蝕強(qiáng)度分級(jí)柵格統(tǒng)計(jì)

項(xiàng)目區(qū)域平均土壤侵蝕模數(shù)為：占比最大的等級(jí)為輕度侵蝕占比50.40%，其次是微度侵蝕占比21.94%，中度侵蝕以上占比達(dá)到27.66%，整個(gè)研究區(qū)侵蝕強(qiáng)度較大。

3 結(jié)論

本文提出的方法針對(duì)礦區(qū)水土保持研究特點(diǎn)為：結(jié)合坡長(zhǎng)、坡度數(shù)據(jù)及將年平均降雨量等因子換算到同投影坐標(biāo)系下，計(jì)算各等級(jí)土壤侵蝕等級(jí)面積占比，結(jié)果依次為微度侵蝕占比21.94%、輕度侵蝕占比50.40%、中度侵蝕占比12.90%、強(qiáng)度侵蝕占比9.04%、極強(qiáng)度侵蝕占比5.31%、劇烈侵蝕占比0.40%。從各等級(jí)占比中可得礦區(qū)水土保持度較差。此外無(wú)人機(jī)影像不僅詳細(xì)展現(xiàn)礦區(qū)內(nèi)的地表地物，還可通過(guò)DEM 計(jì)算礦區(qū)土方量，在多視角下為礦山生態(tài)監(jiān)測(cè)提供直接的信息支撐。