回采工作面地質構造預測研究

郝治朝  

（同煤集團地質勘測處，山西，大同，037003）

摘要：本文引入空間插值方法將回采區(qū)域格網化，并以此格網為基礎，結合生產實際揭露的煤層底板標高，借助數學分析方法，預測回采工作面前方的地質構造。

關鍵字： 空間插值  格網  構造  預測

1 概述

由于礦井小斷層構造發(fā)育規(guī)律認識不清，探查手段落后，構造預測準確性較差，給煤礦生產和建設帶來了很多困難。同時，因對小斷層構造認識不清而導致底板突水，甚至淹沒礦井的現象也時有發(fā)生。因此，必須對礦井小斷層構造井進行探查和預測，并盡量提高預測的準確性。利用數學地質中統計分析的方法，通過對已被揭露的小斷層構造的分析研究，經數學計算，建立一種數學模型，用它來擬合礦井小斷層構造的出現和展布規(guī)律，可以達到對礦山小斷層構造定量的認識。

隨著GIS技術的不斷發(fā)展以及研究中對空間高質量數據的要求，空間數據插值應用越來越廣。通常大家通過曲面擬合來描述煤層底板，這在實際應用中受到各種地質條件的限制，本文采用了更為靈活的格網插值方式，對于斷層影響下的不連續(xù)曲面達到了更符合客觀的擬合。以此為基礎，可以借助數據分析方法更為科學地預測回采工作面前方的地質構造。

2 空間格網

2.1三角格網

所謂三角形網格化也就是對平面點集的計算其三角剖分，具體實現過程是將平面點集各點用互不相交的直線段連接起來，構成一個三角形網格。在這種情況下，三角網中各三角形頂點仍為原始數據點。在追蹤等值線時，以其中的每個三角形為單元，進行等值點的計算及等值線的追蹤。

一般而言，三角形網格構成方式有多種，不同的三角化方式有不同的標準。目前被普遍接受和廣泛采用的三角剖分有最小權三角剖分、貪婪三角剖分和Delaunay三角剖分。最小權三角剖分是指該離散點集可以夠成的所有三角剖分中的邊總長最小的三角剖分；而貪婪三角剖分是指一次加入三角剖分的一條邊，直至產生要求的邊數構成的三角剖分；Delaunay三角剖分需要滿足這樣特性即在三角剖分時獲得的任意三角形的外接圓內不包含原始離散點集中的其它點。

繪制等值線時，為了獲得較好的等值點內插精度，總希望最臨近的數據點構成三角形，即三角網中三角形的邊長之和最小。最小權三角剖分滿足這一準則，但其求解算法問題至今尚未解決。而貪婪三角剖分是常以局部情況為基礎作最優(yōu)選擇而不追求最優(yōu)解，它只希望得到較為滿意解的方法，而不考慮各種可能的整體情況。因此從理論上來講Delaunay三角剖分被認為是可實現的最優(yōu)剖分算法。

2.2矩形格網

所謂的規(guī)則格網即依據一定的插值原理將原始離散數據點內插為按一定的行距和列距分布在相互垂直的兩組平行線的交點上的點數據集。其插入原理是：假設有n個測量數據點(Z₁，Z₂…Zn)，其中Z_i=(x_i，y_i，z_i)，把這些數據點投影到XOY平面上，得到一個相應的點集，記為{Z_i(x_i，y_i) | i=1，2…n}，進而以此點集為依據在XOY平面上按照選定的插值方法建立一定數量行列的規(guī)則矩形網格，每一個網格節(jié)點即為一個插入點。

對于矩形格網建立采用的插值方法已比較成熟，包括反距離加權、克里金、趨勢面法、改進謝德別法、移動平均法等。雖然矩形網格法繪制等值線圖具有算法簡單的優(yōu)點，但它要求原始數據點分布較為均勻或均勻。當計算格網中左上部分節(jié)點特征值時，由于參考數據較為鄰近結果也較為精確，而計算格網右下方向部分節(jié)點特征值時，由于原始數據點稀疏，參考點與當前計算點相關度較低導致結果可能與實際相差較大。

2.3兩種格網比較

通過以上對三角網和規(guī)則格網的介紹可知，就對原始數據點的要求而言，三角網格法不僅適合于規(guī)則分布的原始數據點，也適合于不規(guī)則、甚至畸形分布的原始數據點。而矩形網格法要求采集的原始離散數據點分布較為均勻，被研究的空間變量連續(xù)分布，不如三角網格法靈活。就格網化效率而言，在實際網格化過程中，由于一般采集的原始離散數據點分布不均勻，因此使用三角網格法可以省去將原始離散數據點變換成矩形網格節(jié)點數據的處理過程，格網化效率較高。就插值準確度而言，由于三角網格的網格節(jié)點即是原始數據點，在格網點的插值過程中也是直接依據原始數據計算，因此在原始數據點分布比較均勻的情況下可以生成格網點分布規(guī)律且數值精度較高的格網。更重要的是相對于規(guī)則矩形格網，三角網構建較為靈活，可以方便地將斷層線等線性約束條件嵌入到三角網中，對繪制支持斷層等地質構造有重要意義。

3 格網化中斷層處理

3.1 生成delaunay三角網

根據逐點插入法及優(yōu)化算法的思想，生成Delaunay三角網的基本步驟如下：

（1）檢索數據集查找到數據點的最小與最大橫縱坐標分別記錄為xmin、xmax、ymin、ymax，從而確定包納所有的邊界點和散點的超級三角形頂點坐標，并將該三角形作為Delaunay三角網的第一條記錄。

超級三角形頂點的記錄順序是順時針的，而且采用逐點插入法生成的Delaunay三角網中所有的三角形將均以順時針方向記錄頂點。這在斷層識別中有重要意義。

（2）從離散數據中取一點，遍歷已有三角網，判斷該離散點是否位于當前三角形外接圓內。如果是，將當前三角形各邊放入優(yōu)化隊列，并從三角網內刪除該三角形；當遍歷完畢，如果散點不在任何一個三角形中，則表明點落在了計算區(qū)域之外，應刪除此點而不添加到三角網中。反之，按如下方法進行優(yōu)化：

（a）從優(yōu)化隊列中取出一條邊，開始優(yōu)化此邊；

（b）如果此邊屬于邊界邊，則此邊不用優(yōu)化；

（c）如果此邊所在的正在被優(yōu)化的三角形的外接形的除此邊外的另一個頂點。在此情況下應將當前的優(yōu)化邊從優(yōu)化隊列中刪除，同時將該三角形的另外的兩邊加入到優(yōu)化隊列中刪除后送入優(yōu)化后的隊列，便于建立新的三角形；

（d）重復（a）、（b）、（c）步，直到優(yōu)化隊列為空。

（3）重復執(zhí)行2，直到所有的離散點處理完畢，得到要繪制的Delaunay三角網。

在步驟2執(zhí)行過程中，最重要的問題是判斷點與三角形外接圓的關系。解決方法如圖3.1所示，當判斷點Q與三角形ABC外接圓關系時，首先計算三角形的外接圓圓心O點坐標，進而比較OQ距離與半徑R的大小。

圖3.1點與三邊形外接圓關系示意圖

3.2嵌入斷層

嵌入斷層即將斷層線作為約束條件加入已生成的三角網，實現所有斷層線線段均作為三角網中某三角形一邊，以便于計算等值點和追蹤等值線。計算方法如下：

（1）遍歷Delaunay三角網，檢索并記錄作為三角形邊的所有斷層線線段，對比之下得到斷層線段中非三角網某邊的斷層線線段并為記錄；

（2）計算線段中點，并將各中點逐次插入原三角網中，并將該中點按其位置記入斷層線；

（3）執(zhí)行1、2直到斷層線中各線段均為Delaunay三角網的一邊；

3.3斷層線節(jié)點賦值

本文斷層線節(jié)點的賦值思想是將斷層線節(jié)點特征值與其所處的三角形直接關聯，即在包含斷層節(jié)點的三角形中直接記錄斷層線節(jié)點值。此時，以三角網為基礎追蹤等值線時可以采用通用方法完成包含斷層線的等值線追蹤。如在圖3.2中紅色的線表示斷層線，紅色的點代表斷層線節(jié)點，而黑色的點表示其它數據點。當為斷層線節(jié)點B點賦值時，需在△ABE、△ABJ、△BCE、△BCK、△BKJ中分別記錄B點特征值。

解決斷層線節(jié)點賦值的關鍵問題在于搜索合理的參考點以計算特征值，由于斷層兩側數據的無關性，當為某個斷層線節(jié)點賦值時需要分別以該斷層線兩側數據點集作為參考計算特征值，并將所得的特征值分別賦予其所在三角形。求解斷層線與等值線交點關鍵是為了顯示等值線分布趨勢而精度要求并非十分精確，因此本文中斷層線節(jié)點特征值計算時選擇參考點為與其直接相關的離散數據點或矩形格網節(jié)點，即在三角網中與其共存于同一三角形中的其它數據點。如圖3.2中，當計算斷層線節(jié)點A的特征值時， D、E、G、I、J五點為參考點。

圖3.2 含斷層的約束三角網示意圖

獲得參考數據點集后，需要按其與斷層位置關系進行分類。為了便于為三角形中斷層線節(jié)點賦值，本文采用的方法是對當前斷層線節(jié)點為其頂點的三角形進行分類，進而求解每組三角形中斷層線節(jié)點的特征值。如圖3.2 中，當計算B點特征值時，首先檢索出與包含B點的三角形分別為△ABE、△ABJ、△BCE、△BCK、△BKJ。 進而按其與斷層線的位置關系將△ABE和△BCE分為一組，△ABJ、△BCK、△BKJ劃為另一組。在第一組中特征值計算參考的數據點為J、K點， 而第二組中參考的數據點為E點。獲取參考數據點后，即可采用一定插值方法（如克里金、反距離加權）計算特征值并賦予對應的三角形頂點。

本文分類算法的關鍵思想是當兩三角形的公用頂點為非斷層線節(jié)點時，此三角形即劃為一類。該方法的正確性用反證法即可證明，若兩三角形的公用頂點為非斷層線節(jié)點時，且這兩個位于斷層線兩側，則其中一個三角形的該非斷層線節(jié)點與其第三個頂點的連線必與斷層線相交，這與三角網的基本特性矛盾。應用此方法需要注意一種特殊情況如圖3.3所示，同樣紅色的線表示斷層線，紅色的點代表斷層線節(jié)點，而黑色的點表示其它數據點。當計算斷層線節(jié)點B點時，需△ABE、△ABJ、△BCE、△BCK、△BKJ分組，其中△ABE為獨立為一組且該三角形頂點均為斷層線節(jié)點，此時參考的數據點為與該三角形公用當前斷層線節(jié)點對邊的三角形中的非斷層線節(jié)點，在圖3.3中所示為D點。

圖3.3 斷層線節(jié)點特征值計算特殊情況圖示

總結斷層線節(jié)點計算步驟如下：

（1）檢索Delaunay三角網，將所有包含斷層線節(jié)點的三角形記錄成集；

（2）在斷層線節(jié)點中取出一點；

（3）檢索1中的三角形集，記錄并存儲與當前斷層線節(jié)點所處的三角形；

（4）將與相前斷層線節(jié)點相關的三角形按位于斷層線方位進行分類并檢索獲取每一類的參考點集；

（5）分別計算兩側的斷層節(jié)點特征值，并賦予相應的三角形；

（6）重復執(zhí)行2-5，直到所有的斷層線節(jié)點處理完畢。

4 地質構造預測方法

4.1小型褶曲的預測

首先通過整個或局部格網確立一個平面，進而分析實際揭露各點煤層底板標高到該平面的距離在某一方向（縱向或橫向）上的數據分布的相關關系，如果呈現出增加到峰值然后降低或減少到谷底然后增加，且近似對稱分布，則是褶曲出現的標志。

平面法線式參照方程：

Xcosα+ycosβ+zocsγ-p0=0

平面的確定要求與煤層產狀面擬合較好，否則預測準確率較低。因此可以選用局部格網作為擬合平面確立依據，以增加擬合精度。與此同時，需確定一個判斷指標以評價其擬合準確度。

4.2隱伏斷層預測

煤強度（P）相對而言降低和節(jié)理發(fā)育強度（J）相對增大，是采掘前方出現斷層的兩個重要征兆。

據此預測步驟如下：

1）對已知斷層的觀測數據建立經驗公式

Sq=K1H+b1

Sj=K2H+b2

2）連續(xù)監(jiān)測

當P≥A或J≤B時，視為正常征兆；

當P＜A或J＞B時，則為異常征兆，此時記錄Pm1和Jn1點在井下的位置或三維坐標，特別要注意Pm2、Pm3和Jn1、Jn2的觀測。

3）動態(tài)預測

f₀=(P≥A)∩(J≤B)

f₁= P_m2∩J_n1∪P_m1∩J_n2

f₂= P_m1∩J_n3∪P_m2∩J_n2

f₃= P_m2∩J_n3∪P_m3∩J_n2

f₄= P_m3∩J_n3

f真值表

4.3半隱伏及完全隱伏斷層發(fā)育規(guī)模預測

半隱伏及完全隱伏斷層發(fā)育規(guī)模的預測依據如下公式：

α——巖層真傾角；

β——斷層真傾角；

θ——斷層傾向與剖面方位夾角；

ω——巖層走向與斷層走向的夾角；

φ——巖層傾向與剖面方位夾角；

ρ——斷層走向線與斷煤交線在斷層面內的夾角；

γ——斷層面綜合擦痕伏角。   

5 地質構造預測實例

同煤集團同發(fā)東周窯礦地質構造極其復雜，在首采面順槽掘進的過程中多次遇到斷層、陷落柱等地質構造，對正常的生產建設造成很大影響，甚至一度因掘全巖而建設進度緩慢。本文研究以此礦資料為基礎。

格網化依據的數據來源于兩方面，一是勘探過程中的鉆孔資料，二是來源于巷道掘進過程中實測的煤層底板標高。首先以此基礎數據生成經格網化后生成格網，進而在此基礎上，借助上述數據分析方法預測掘進、回采前方的地質構造，收到了較好的效果。

表5.1 工作面兩巷揭露標高

 

如上表所示為工作面兩巷揭露的部分煤層底板標高，以此數據為基礎生成格網，進而繪制煤層底板等高線如圖5.1所示。在巷道掘進過程中揭露的煤層底板標高與理論值差異如圖5.2所示，進而可以預測該區(qū)域斷層的存在。

 

圖5.1 以格網為基礎繪制的底板等高線圖

圖5.2 掘進揭露與理論擬合的差異

6 結論

本文首先提出了采用格網化的方式取代曲面擬合來描述煤層底板面，進而研究了斷層等地質構造存在條件下的格網化賦值方式。并以此格網為基礎，借助數據分析方法，預測回采工作面前方褶曲、斷層等構造的存在及其規(guī)模，為地質構造定量化描述和預測提供了依據。為完成上述成果，本文在一些方面做了探索。如以斷層線結點賦值的方式，考慮了格網化過程中已知斷層對數值分布的影響，為進一步分析提供了更為可靠的基礎格網。本文在褶曲的預測中，基于局部格網選擇擬合平面，避免了以整個區(qū)域擬合平面造成的誤差較大的問題。

參考文獻

1.崔洪慶．煤礦生產地質保障技術[M]．北京：煤炭工業(yè)出版社，2005

2.喬美萍，周榮福等．基于GIS的空間插值法在礦井地質構造預測中的應用研究[J]．礦業(yè)工程，2010(2)

3.李志華，竇林名，張衛(wèi)東．礦井地質構造預測與應用[J].煤炭工程．2005(2)

4王國華．煤礦中小型地質構造預測方法[J]．水力采煤與管道運輸．2006(3)

 

Research on geological structures in working face of coal mine

Abstract:This paper introduced spatial interpolation method into regional grid transformation, and then based on fruiting grid,combined with the floor level of production practice and using mathematical analysis method, to predict the geological structure of front mining face.

 

Key words:Spatial interpolation;Grid;geological structure; prediction

附件下載