本發(fā)明涉及水土保持學和恢復生態(tài)學領域，特別涉及一種礦山排土場邊坡土壤侵蝕速率估算方法。
背景技術：
開發(fā)建設項目導致的水土流失是人類生產(chǎn)建設活動過程中擾動地表和地下巖層、堆置廢棄物、構筑人工邊坡以及排放各種有毒有害物質而造成的水土資源和土地生產(chǎn)力的破壞和損失，是一種典型的人為加速侵蝕。礦山排土場由采礦區(qū)巖土剝離、運輸、堆墊而形成，原來的土體和礦層的上覆巖層經(jīng)過劇烈的擾動混合后，以松散堆積體狀態(tài)堆置在內外排土場，形成人工巨大松散堆積地貌景觀，平臺-陡坡是其基本地貌單元，其水土流失分布是點、線、片、帶、塊、面中的一種或多種形式的組合。排土場土壤侵蝕過程復雜，主要表現(xiàn)在土壤侵蝕部位集中、類型多樣、強度劇烈。排土場堆置厚度各部位不等，顆粒組成差異大，自然固結率不同，在自重營力及入滲水流的作用下，造成不同部位壓縮沉降速率不同，即所謂的“非均勻沉降”。
目前，主要通過定性和定量的方法測算礦山排土場工程堆積體土壤侵蝕速率。定量研究方法主要有徑流小區(qū)法、侵蝕測針法、元素示蹤法、模型模擬法。近年來，建設工程項目區(qū)土壤侵蝕的研究方法主要采用人工降雨、放水沖刷或二者組合以及天然降水條件下小尺度上的短期定位觀測或模擬實驗為主。水土流失量預測的方法有數(shù)學模型法(如USLE模型)、類比分析法(如毗鄰項目類比)以及典型調查推算法(如類似對象調查推算法)。
目前，針對無實測數(shù)據(jù)地區(qū)，如何快速精確計算礦山排土場邊坡土壤侵蝕模數(shù)尚無可靠的方法。
技術實現(xiàn)要素：
本發(fā)明提供了一種礦山排土場邊坡土壤侵蝕速率估算方法，以解決現(xiàn)有技術中無法快速精確計算礦山排土場邊坡土壤侵蝕模數(shù)的問題。
本發(fā)明中的礦山排土場邊坡土壤侵蝕速率估算方法，包括：步驟1，利用無人機獲取礦山排土場影像數(shù)據(jù)；步驟2，利用RTK-GPS獲取排土場的地面控制點準確的空間位置信息；步驟3，優(yōu)選地可利用Agisoft Photoscan professional 1.1.2軟件，根據(jù)所述GPS位置信息及所述影像數(shù)據(jù)，生成排土場點云數(shù)據(jù)和正射影像數(shù)據(jù)；步驟4，對所述正射影像數(shù)據(jù)進行目視解釋以生成侵蝕溝邊界，根據(jù)所述侵蝕溝邊界篩選出位于侵蝕溝內部的高程點數(shù)據(jù)并根據(jù)高程點數(shù)據(jù)生成侵蝕溝DEM；步驟5，根據(jù)所述侵蝕溝邊界得到緩沖區(qū)邊界，篩選位于所述緩沖區(qū)邊界內的高程點數(shù)據(jù)并根據(jù)高程點數(shù)據(jù)生成緩沖區(qū)DEM；步驟6，利用柵格計算器工具計算所述緩沖區(qū)DEM與所述侵蝕溝DEM的差從而得到侵蝕溝平均深度DEM；步驟7，根據(jù)所述侵蝕溝平均深度DEM及侵蝕溝面積計算得到侵蝕溝的侵蝕土壤體積；步驟8，根據(jù)下式計算得到土壤侵蝕模數(shù)：
其中，Ms為土壤侵蝕模數(shù)，單位t/(km2·a)；Ws為年侵蝕總量，單位t，且Ws為所述侵蝕土壤體積與土壤容重的乘積；F為侵蝕面積，單位為km2；T為侵蝕時限，單位為a。
優(yōu)選地，所述步驟1還包括：規(guī)劃無人機航跡以使影像間的重疊度約為40％；利用無人機獲取排土場清晰無云的影像，且該影像包含GPS信息。
優(yōu)選地，所述GPS位置信息采用RTK-GPS獲得。
優(yōu)選地，所述地面控制點選自地面標志點。
優(yōu)選地，所述地面標志點為道路交叉點、和/或水渠交叉點、和/或房屋轉角。
優(yōu)選地，所述侵蝕溝DEM、緩沖區(qū)DEM及侵蝕溝平均深度DEM采用ArcGIS軟件得到。
本發(fā)明適合于無長期監(jiān)測數(shù)據(jù)地區(qū)，可通過無人機影像和高精度DEM數(shù)據(jù)之間的結合，快速推算出排土場邊坡土壤侵蝕速率，并產(chǎn)生一定的經(jīng)濟、社會和生態(tài)效益。
具體實施方式
本發(fā)明涉及一種快速估算礦山排土場邊坡土壤侵蝕速率方法，具體涉及利用無人機影像生成高密度點云和正射影像，監(jiān)測礦山邊坡微地形變化的方法，用于估算礦山排土場邊坡土壤侵蝕速率，其通過使用無人機遙感影像，獲取研究區(qū)高程密集點云，生成同一時期較為精確的兩組DEM數(shù)據(jù)，通過目視解譯，提取侵蝕溝邊界，以此來計算土壤侵蝕模數(shù)的方法。
在一個實施例中，本發(fā)明中的礦山排土場邊坡土壤侵蝕速率估算方法，包括以下步驟：
步驟1，利用無人機獲取礦山排土場影像數(shù)據(jù)。針對礦山排土場坡陡、堆積物松散，人員難以全面實地勘察的特點，利用無人機合理規(guī)劃航跡，獲取排土場清晰無云的影像，影像需要包含GPS信息，影像間重疊度應在40％左右，便于后期影像拼接及數(shù)據(jù)生成。
步驟2，采用RTK-GPS獲取排土場的地面控制點的GPS位置信息。例如，使用RTK-GPS在排土場測得地面實際控制點坐標，控制點一般選擇在道路交叉點、水渠交叉點、房屋轉角點等在影像上易分辨且較精細的特征點。
步驟3，根據(jù)所述GPS位置信息及所述影像數(shù)據(jù)，生成點云數(shù)據(jù)和正射影像數(shù)據(jù)。例如，可使用Agisoft Photoscan軟件加載無人機拍攝的照片，導入RTK-GPS測得的控制點，軟件自動將照片對齊，生成高密度點云數(shù)據(jù)和正射影像數(shù)據(jù)。
步驟4，對所述正射影像數(shù)據(jù)進行目視解釋以生成侵蝕溝邊界，根據(jù)所述侵蝕溝邊界篩選出位于侵蝕溝內部的高程點數(shù)據(jù)并根據(jù)高程點數(shù)據(jù)通過空間插值生成侵蝕溝DEM；
步驟5，根據(jù)對所述侵蝕溝邊界進行緩沖區(qū)分析得到緩沖區(qū)邊界，分別篩選位于所述緩沖區(qū)邊界內的高程點數(shù)據(jù)并根據(jù)高程點數(shù)據(jù)通過空間插值生成緩沖區(qū)DEM；
步驟6，利用柵格計算器工具計算所述緩沖區(qū)DEM與所述侵蝕溝DEM的差從而得到侵蝕溝平均深度DEM；
步驟7，根據(jù)所述侵蝕溝平均深度DEM及侵蝕溝面積計算得到侵蝕溝的侵蝕土壤體積V(m3)，其中：
Vi＝Si*H_meani
其中，Vi為各侵蝕溝侵蝕土壤體積(m3)，Si為各侵蝕溝面積(m2)，H_meani為各侵蝕溝的侵蝕溝平均深度DEM(m)；
整個研究區(qū)的土壤侵蝕總體積為：
V總＝V1+V2+V3+…+Vn+…+Vi。
步驟8，根據(jù)下式計算得到土壤侵蝕模數(shù)(單位面積上每年侵蝕土壤的平均重量)：
其中，Ms為土壤侵蝕模數(shù)，單位t/(km2·a)；Ws為年侵蝕總量，單位t，且Ws為所述侵蝕土壤體積與土壤容重的乘積，即Ws＝τ*V總，其中，τ為土壤
容重1.5t/m3；F為侵蝕面積，單位為km2；T為侵蝕時限，單位為a(年)。
由于采用了上述技術方案，本發(fā)明適合于無長期監(jiān)測數(shù)據(jù)地區(qū)，可通過無人機影像和高精度DEM數(shù)據(jù)之間的結合，快速推算出排土場邊坡土壤侵蝕速率，并產(chǎn)生一定的經(jīng)濟、社會和生態(tài)效益。
優(yōu)選地，所述侵蝕溝DEM、緩沖區(qū)DEM及侵蝕溝平均深度DEM采用ArcGIS軟件得到。
本發(fā)明的有益效果是：
1、科學性
以水土保持學和土壤學原理為依據(jù)，針對礦山排土場邊坡無長時間監(jiān)測數(shù)據(jù)的情況，提出了應用無人機獲取影像，獲取高密度點云數(shù)據(jù)，通過人工目視解譯，篩選地面高程數(shù)據(jù)點，對比了應用不同插值方法生成DEM數(shù)據(jù)的精度，研究結果表明該方法科學有效。
2、可操作性強
本發(fā)明只需要使用無人機獲取研究區(qū)影像及地面控制點坐標，通過軟件即可完成影像自動配準，生成高密度點云數(shù)據(jù)和正射影像。通過人工目視解譯和數(shù)據(jù)空間分析，即可計算無長時間監(jiān)測數(shù)據(jù)地區(qū)的土壤侵蝕速率，操作簡單，具有較強的可操作性。
現(xiàn)有技術的土壤侵蝕監(jiān)測方法費時、費力，測量誤差具有不確定性，而本發(fā)明可利用無人機快速獲取高精度影像，并采用數(shù)字高程模型(DEM)快速獲取地形數(shù)據(jù)，這樣，可通過無人機遙感監(jiān)測技術快速獲取高精度的地表影像，提取高密度點云，生成正射影像和高精度DEM數(shù)據(jù)，從而能夠有效監(jiān)測地形變化。