張寶通

（中鐵工程設計咨詢集團有限公司濟南設計院，山東 濟南 250000）

礦山開采對地質環(huán)境與生態(tài)環(huán)境的影響較大，為了合理地開采礦山，應對礦山開采區(qū)周圍的地質進行動態(tài)監(jiān)測。通過了解礦山地質周圍的環(huán)境問題，監(jiān)測地質環(huán)境的變化，對礦山開采區(qū)進行合理地評估與預測，最終達到合理化開采，維護礦山周圍的生態(tài)環(huán)境穩(wěn)定發(fā)展[1]。

礦山開采區(qū)的活動受多種因素的影響，例如巖石、土壤、地質作用等，與生物圈之間密切聯(lián)系。開采礦山時，周圍巖體受到外力作用導致平衡狀態(tài)出現(xiàn)失衡現(xiàn)象。在巖體調整平衡過程中，巖土層在一定時間內會出現(xiàn)位移、變形的情況，從而形成不均勻沉降、塌陷坑等斷裂，長時間地斷裂，在開采區(qū)上部會形成較多裂縫，容易導致山體崩塌，形成地質災害[2]。礦山開采對水源也存在一定的影響，開采區(qū)上部斷裂破壞了水資源的平衡，使地下水滲透，水位下降，嚴重時甚至會引發(fā)水資源干枯[3]。礦山開采區(qū)在開采中形成的礦坑水、工業(yè)廢水等自然排放到河流中，致使開采區(qū)附近的水體與土壤污染嚴重。礦山開采產生的廢水中化學物質較多，流入河流土壤中，造成該區(qū)域的地方病源，存在嚴重的安全隱患?；诖耍岢隽说V山開采區(qū)周圍不均勻沉降監(jiān)測分析方法，對礦山開采區(qū)的地質變化進行監(jiān)測，減少因開采礦區(qū)導致的一系列問題。

1 礦山開采區(qū)周圍不均勻沉降監(jiān)測分析方法設計

針對礦山開采區(qū)地表沉降速度大和非線性變化的特征，設計了礦山開采區(qū)周圍不均勻沉降監(jiān)測分析方法。通過構建短基線礦山地表沉降監(jiān)測模型、對不均勻沉降監(jiān)測數(shù)據(jù)的預處理、分析監(jiān)測地表移動變形規(guī)律，最終實現(xiàn)對礦區(qū)地質變化的監(jiān)測，具體步驟如下。

1.1 構建短基線礦山地表沉降監(jiān)測模型

短基線礦山地表沉降監(jiān)測模型的構建依據(jù)是在短時間基線與短空間基線原則上，選擇沉降干涉數(shù)據(jù)，根據(jù)空間失相原理所構建而成。由于大氣波動變化較大，礦山開采區(qū)周圍的地形相位與形變相位之間的相位精度變化較明顯，應當對兩個相位進行干涉圖處理。使用D-InSAR干涉處理獲取短基線內地表沉降的時空變化，在礦山開采區(qū)選定好相應的位置后，獲取相差分干涉數(shù)據(jù)。由于SAR礦區(qū)影響受到失相干的影響，在選擇區(qū)域時，兼顧相干性與相位梯度，建立地表相差干涉像對組合結構圖，如圖1所示。

圖1 短基線礦山地表沉降監(jiān)測結構

根據(jù)圖1可知，確定礦山開采實驗區(qū)，當?shù)乇沓两禎M足短基線原則時，對數(shù)據(jù)對進行配準操作，有DEM數(shù)據(jù)生成干涉圖，為差分失相處理提供基礎。

1.2 不均勻沉降監(jiān)測數(shù)據(jù)預處理

基于上述構建的短基線礦山地表沉降監(jiān)測模型，對礦山開采區(qū)的不均勻沉降監(jiān)測數(shù)據(jù)進行預處理。根據(jù)礦區(qū)的開采方式不同，礦區(qū)地表的坑口沉降也不同。

各級斷裂受地質構造分割影響，開采過程中的連續(xù)性較差。受長期的地表擾動狀態(tài)影響，地表沉降與變形的持續(xù)時間較長，導致某些斷層出現(xiàn)活化現(xiàn)象，加劇了礦區(qū)地表沉陷災害。根據(jù)礦區(qū)地表的面積，對獲取的波段變化數(shù)據(jù)進行計算，得到4＊4弧秒的DEM數(shù)據(jù)，在空間分辨率為110m時，每隔4.5弧秒采集一次礦區(qū)地表變化數(shù)據(jù)。采集的數(shù)據(jù)水平參考系使用WAG65，垂直參考系選擇BMN74水平面。在礦區(qū)地質環(huán)境相對平坦時，DEM數(shù)據(jù)對于監(jiān)測礦區(qū)地表的不均勻沉降結果更加精準。

為了提高基線數(shù)據(jù)的精度，精密的軌道數(shù)據(jù)處理是關鍵環(huán)節(jié)。為了避免出現(xiàn)數(shù)據(jù)缺失的情況，將DEM模擬干涉相位進行轉化，轉化為接近SAR圖像數(shù)據(jù)的分辨率，得到UTM坐標下的SAR圖像數(shù)據(jù)。在DEM數(shù)據(jù)干涉相位中，會產生因DEM數(shù)據(jù)精度的地形變化誤差，通過計算可知，地形誤差對干涉相位的影響為：

其中，B表示地表沉降垂直基線；Δh表示沉降高程誤差；R表示地面點到地表沉降的距離；λ表示信號變動波長；θ表示地表沉降入射角。其他因素對相位影響忽略不計時，沉降高程誤差引起的地表形變誤差較大，應當采取有效的方法減小地表沉降的高程誤差。

1.3 分析監(jiān)測地表移動變形規(guī)律

根據(jù)監(jiān)測的礦山開采區(qū)的數(shù)據(jù)結果，沿礦山開采區(qū)的任意方向做剖面動態(tài)下沉分析。利用地表變化奇異值分解能夠得到礦區(qū)的平均沉降速率，對各時間段的速度進行記錄，獲取速度變化規(guī)律，根據(jù)時間的變化求得基于時間序列的沉降結果。距離礦山開采區(qū)域較近的區(qū)域，其不均勻沉降的速率較快；遠離礦山開采區(qū)的位置，不均勻沉降的速率較慢。并且對于監(jiān)測到的礦山地表的變形規(guī)律中，不僅存在不均勻沉降的形式，地表還存在抬升形變的可能，出現(xiàn)抬升現(xiàn)象的區(qū)域大多集中分布在30～40mm/a范圍內，稱為礦區(qū)中的漏斗位置。

礦山開采區(qū)周圍的不均勻沉降還受到降雨量的影響，一般在夏季，降水量較大時，地表會出現(xiàn)大面積的沉降形變，據(jù)分析，最大的沉降值可達345mm。通過監(jiān)測礦山開采區(qū)地表變化的數(shù)據(jù)，最終得到地表移動變形規(guī)律。

2 對比實驗

選擇唐山市某礦山開采區(qū)作為研究對象，采用DORMT軌道測量數(shù)據(jù)，以單視影像數(shù)據(jù)為源數(shù)據(jù)。根據(jù)短基線時序SAR方法，處理覆蓋唐山市礦區(qū)地面沉降的SAR數(shù)據(jù)。依據(jù)實驗區(qū)失相相位的散射特點，為了將監(jiān)測過程操作簡便，采用閾值法確定礦山開采區(qū)位置。地表沉降離差指數(shù)為相干目標，對干涉圖的失相不均勻沉降基于高程不同進行分析。在分解干涉對組合時，采用短基線方式，設置空間基線為350m，時間基線閾值取值范圍介于200～250d，根據(jù)地表沉降的入射角與垂直基線，計算出相應的不均勻沉降高程差。將本文提出的礦山開采區(qū)周圍不均勻沉降監(jiān)測分析方法與傳統(tǒng)的雷達測量地面點分析方法的沉降高程誤差結果進行對比，如表1所示。

表1 兩種分析方法的地表沉降高程誤差對比

根據(jù)表1可知，本文提出的礦山開采區(qū)周圍不均勻沉降監(jiān)測分析方法監(jiān)測的沉降高程差較傳統(tǒng)方法相比誤差較小，更接近實際結果。

因此，本文提出的分析方法更適用于礦山開采區(qū)周圍不均勻沉降監(jiān)測，較精確的地表沉降高程誤差有利于礦山的合理開采，保證礦山周圍地質環(huán)境的平衡。

3 結束語

本文研究了礦山開采區(qū)周圍不均勻沉降監(jiān)測分析方法，通過構建短基線礦山地表沉降監(jiān)測模型、不均勻沉降監(jiān)測數(shù)據(jù)的預處理、分析監(jiān)測地表移動變形規(guī)律，最終完成監(jiān)測。此分析方法能夠降低模擬監(jiān)測結果的誤差，使監(jiān)測值更加接近于實際值。

但由于SAR數(shù)據(jù)獲取有限，故監(jiān)測結果的準確率還有待提高。在未來的研究中，可以根據(jù)礦山開采區(qū)的點位時序數(shù)據(jù)分析，對這一問題進行改善。