花崗巖礦儲量變化監(jiān)測研究

盧 偉

（廣東省地質局第六地質大隊，廣東江門529000）

1 概述

隨著我國城鎮(zhèn)建設的飛速發(fā)展，工程建造的原料需求量也不斷增大。其中的砂石料，作為各地方高速路、高鐵、建構筑物等基礎建設中不可或缺的物料，其需求量也在急劇增大。在如此大的市場需求環(huán)境下，石料的出廠價相比往年出現(xiàn)了大幅上漲，其利潤空間之大，導致各地開設了許多采石場，由此進入采石賺錢的時代。當然，為了貫徹資源開采的合理化與可持續(xù)化發(fā)展，政府定期要對采石場進行資源儲量核實工作，同時遞交儲量核實報告，其目的是調查礦山開采現(xiàn)狀，估算礦山保有資源量，為下步開發(fā)提供地質依據(jù)。在這樣的發(fā)展背景下，如何運用新的測量技術來快速有效地對儲量變化總量進行計算，成為目前一個研究的熱點與難點。為此，本文依托于廣東省成熟的CORS系統(tǒng)、無人機成像系統(tǒng)、免棱鏡全站儀系統(tǒng)相結合，實現(xiàn)了快速準確的計算土石方量，為政府決策與部署提供了強有力的數(shù)據(jù)支持。

2 數(shù)據(jù)獲取與處理

2.1 圖根控制測量

廣東省連續(xù)運行衛(wèi)星定位服務系統(tǒng)（簡稱GD?CORS）于2006年11月開始試運行，至今已有十三年多，實時坐標定位服務高效穩(wěn)定，能夠提供平面毫米級定位精度，高程厘米級精度。廣東省似大地水準面精化模型的外符合精度為±0.048m，完全符合1∶2000地形圖測量對高程精度的要求。基于成熟的GDCORS數(shù)據(jù)以及傳統(tǒng)的RTK技術需要收集礦區(qū)附近的控制點數(shù)據(jù)，架設基準站并進行參數(shù)計算，攜帶設備多，花費很多時間與精力的缺點，甚至有些礦區(qū)及周邊無法收集到控制點數(shù)據(jù)，這就給實地采集測量數(shù)據(jù)帶來很多麻煩。因此完全可以基于GDCORS的GNSS-RTK（型號為中海達TS7）直接采集測站點和后視點坐標及高程。為檢驗GDCORS網(wǎng)絡RTK定位精度，選取江門市某區(qū)8個D級GPS控制點及四等水準點進行精度統(tǒng)計（每個點均利用腳架對中采用平滑觀測），GPS網(wǎng)圖如圖1所示，統(tǒng)計精度如表1、表2所示。

圖1 GPS網(wǎng)圖

由表1數(shù)據(jù)可知，點位精度較差均小于3cm，完全滿足圖根控制點測圖的要求。較差均值1.67cm說明該GDCORS系統(tǒng)提供數(shù)據(jù)穩(wěn)定可靠。由表格2數(shù)據(jù)可知，8個點位較差均小于5cm說明GDCORS提供高程數(shù)據(jù)滿足大比例尺測圖要求。綜上證明了基于GD?CORS的網(wǎng)絡RTK測量可以達到厘米級精度，符合礦山動態(tài)監(jiān)測精度要求。為了使礦山石方量計算工作能更加高效、經(jīng)濟，充分利用GDCORS是相當有必要的。

表1 D級GPS控制點與網(wǎng)絡RTK平面成果比對表

2.2 確定地形修測范圍

表2 D等水準點與網(wǎng)絡RTK高程比對

由于大部分采石場開采量大，一般分臺階開采，很難確定監(jiān)測周期間地形變化范圍。為此，本項目采用大疆精靈4無人機進行拍照，每次在固定地點起飛至固定高度，同角度拍攝照片與上一次監(jiān)測照片進行對比可以很輕松地確定需要修測的范圍。表格3為該無人機主要技術參數(shù)。

2.3 全站儀采集地形數(shù)據(jù)

表3 大疆精靈4無人機主要技術參數(shù)

本項目采用中緯ZOOM 35 PRO全站儀進行野外數(shù)據(jù)采集，測角精度為2″，免棱鏡測量模式標稱為1000m(在礦山觀測背景下親測可達到500～700m)，測距精度為3mm+2ppm。由于礦山按照開發(fā)利用方案分臺階開采，人工立標尺不僅效率低下而且安全隱患特別大。因此找好有利測站點，利用全站儀免棱鏡測量功能可以做到一個測站搞定全部修測。如表4所示為全站儀主要技術參數(shù)。

表4 全站儀主要技術指標

3 工程實例

3.1 項目概況

廣東省江門市露天花崗巖礦山數(shù)量多，本單位承接該市兩區(qū)共計18個露天礦山的動態(tài)監(jiān)測工作。這就務必有一套高效率、低成本且安全行之有效的解決方案。本次項目選擇了監(jiān)測項目里礦區(qū)面積最大、年開采量最多并且是該省知名花崗巖礦山的泰盛石場有限公司建筑用花崗巖礦區(qū)（以下簡稱泰盛石場）作為研究對象，討論礦山動態(tài)監(jiān)測的工作方法。研究區(qū)概況敘述如下。泰盛石場行政隸屬該市該區(qū)某鎮(zhèn)管轄，礦區(qū)面積0.4519km2，開采標高：+360～+35m，開采礦種：建筑用花崗巖，開采方式：露天開采，生產(chǎn)規(guī)模：200×104m3/a。

3.2 室內數(shù)據(jù)整理成圖及當期開采儲量計算

3.2.1 確地監(jiān)測范圍

該花崗巖礦山作為露天礦每月采掘量很大，初步估計有幾十萬立方米。本文選擇一年的時間差為基準，計算該礦區(qū)石方量的變化。確定計算范圍根據(jù)大疆無人機航飛確定，如圖2所示范圍確定。

3.2.2 計算石方變化量

本次的處理過程是以該花崗巖礦山的地形地質圖為例，該地形地質圖的比例尺為1∶2000，數(shù)據(jù)的坐標系是1980年西安坐標系，基準面采用1985國家高程基準。用南方CASS 10.1打開2019年1月3日和2020年1月5日該礦的地形地質圖進行比較。對一年內有所變化的區(qū)域進行初步的判斷，并用多段線進行框選，通過CASS工具欄中的工具進行數(shù)據(jù)批量刪減，刪減掉指定區(qū)域外多余的數(shù)據(jù)，實現(xiàn)開采范圍的初步選取。通過CASS工具欄里的“建立DTM”功能，進行DTM的創(chuàng)建，在生成DTM時可由圖面高程點生成，并且考慮在構建三角網(wǎng)過程中陡坎與地性線的影響，之后對得到了三角網(wǎng)進行細部修改，將不合理三角網(wǎng)點刪除，進行局部構建成網(wǎng)。三角網(wǎng)如圖3所示。最后采用南方CASS 10.1自動生成剖面圖，利用平行斷面法計算。計算挖方量如圖4所示。

圖2 確定監(jiān)測范圍

圖3 點數(shù)據(jù)導入CASS中顯示

3.2.3 數(shù)據(jù)比較與分析

數(shù)據(jù)a：通過采石場提供真實出貨量，該數(shù)值作為真值。

數(shù)據(jù)b：它的兩期數(shù)據(jù)均以地形地質圖為基礎，是通過本文提出的快速監(jiān)測方法計算得到的儲量變化總量。

數(shù)據(jù)c：它是以傳統(tǒng)RTK外業(yè)采集模式，利用CASS計算石方量變化值。

由圖5所示數(shù)據(jù)分析可知，數(shù)據(jù)b是采用本文提出的方法快速計算的石方量變化值與真值之間差距微小。數(shù)據(jù)c的獲取是按照傳統(tǒng)外業(yè)RTK模式，外業(yè)采集高程通過CASS軟件計算而來。實際作業(yè)中有一定的危險性。數(shù)據(jù)b數(shù)值為610.79×104m3；數(shù)據(jù)c數(shù)值為661.91×104m3，由(c-b)/b=8.34%，可以看出數(shù)據(jù)b與數(shù)據(jù)c的值之間相差較大，基于本文提出的方法計算的數(shù)據(jù)比傳統(tǒng)監(jiān)測的數(shù)據(jù)要少了8%，并且數(shù)據(jù)b數(shù)值更接近真值。該組數(shù)據(jù)證明基于本文提出的方法計算的石方儲量變化結果要比傳統(tǒng)監(jiān)測結果更為全面，更能夠反映出采石場的真實開采量。

圖4 石方變化量

圖5 數(shù)據(jù)比較圖

4 結語

本文通過使用廣東省CORS系統(tǒng)實時提供平面坐標及高程，且在該地區(qū)信號穩(wěn)定，固定速度快，精度可行，能省時省成本省人力。利用拍照型無人機快速確定修測范圍，可避免重復測繪，節(jié)省時間，并獲得直觀的礦山俯視圖給礦山主管部門參考。長距離免棱鏡全站儀的使用，在外業(yè)觀測中節(jié)省人力物力，提高工作效率。綜合以上三種設備融合一套計算土方量變化的系統(tǒng)。該系統(tǒng)能夠準確及時計算出石方儲量變化數(shù)值。通過某礦的實測數(shù)據(jù)計算結果與采石場提供真實數(shù)據(jù)進行比較，證明了該系統(tǒng)的可靠性與可行性。本文研究是基于露天采礦領域的嘗試，希望這種監(jiān)測方式可以逐漸走向成熟，并在不久的將來被生產(chǎn)監(jiān)測單位所認可。