蘇明金許光泉李佩全劉澤功嚴(yán)家平黃向菁宗云峰

(1.安徽理工大學(xué)地球與環(huán)境學(xué)院,安徽省淮南市,232001; 2.淮南礦業(yè)集團(tuán),安徽省淮南市,232001)

★地質(zhì)與勘探 ★

克立格插值法在礦井生產(chǎn)中的應(yīng)用＊

蘇明金1許光泉1李佩全2劉澤功1嚴(yán)家平1黃向菁2宗云峰2

(1.安徽理工大學(xué)地球與環(huán)境學(xué)院,安徽省淮南市,232001; 2.淮南礦業(yè)集團(tuán),安徽省淮南市,232001)

依據(jù)克立格插值法在張集煤礦某工作面基巖面標(biāo)高實(shí)際控制應(yīng)用情況,分析了影響控制精度的各因素,為勘探各階段合理布置鉆孔提供理論指導(dǎo)。同時采用交叉驗(yàn)證與采樣驗(yàn)證相結(jié)合的方法,對不同控制方案的控制精度進(jìn)行評價,確定合理的評價方案。

礦井測量 克立格插值法 控制誤差因素 工作面基巖面 標(biāo)高控制

AbstractAccording to the application of Kriging interpolation method on the bedrock surface elevation of actual control of the situation in a working face of Zhangji Coal Mine,factors that control the final precision are analyzed,thus providing theoretical guidance for the various stages of the exploration and rational arrangement of drilling.At the same time,a combined utilization of cross-validation and sampling verification methods is carried out in evaluating the accuracy of different control programs and a reasonable evaluation program is determined.

Key wordsmine survey,Kriging interpolation method,factor of error control,working face bedrock surface,elevation control

隨著計算機(jī)在礦井地質(zhì)生產(chǎn)中的廣泛應(yīng)用,很多地質(zhì)計算方法、計算原理被集成在各個軟件包中,地質(zhì)工作者所要做的就是簡單的數(shù)據(jù)輸入與計算結(jié)果的輸出,而忽略了各個軟件的應(yīng)用前提條件、產(chǎn)生誤差因素以及計算結(jié)果精度。本文通過克立格插值法在張集煤礦生產(chǎn)中的應(yīng)用實(shí)例分析了控制誤差因素,并對計算精度在不同計算方案下進(jìn)行交叉驗(yàn)證,說明通過設(shè)計合理的施工方案對生產(chǎn)成本與生產(chǎn)精度控制起著重要作用。

1 克立格插值法

1.1 對采樣點(diǎn)屬性值賦予權(quán)重求估計值

n個采樣點(diǎn)得到n個屬性值(i=1,2,3,…n),通過屬性值來估計某一點(diǎn)上的未知值z＊,可以通過對各個屬性值賦予不同的權(quán)重wi后累加得到:

式中:ui為第i采樣點(diǎn)屬性值。

wi、wj表示點(diǎn)i、j的權(quán)重,如果該點(diǎn)上沒有采樣,則沒法知道真正的誤差大小,所以 (2)式不可解。

1.2 通過距離相關(guān)量求權(quán)重

空間統(tǒng)計方法認(rèn)為估計值與采樣點(diǎn)屬性值的相似程度是通過點(diǎn)對的距離相關(guān)量來度量的。點(diǎn)對距離相關(guān)量只與采樣點(diǎn)間的相互距離有關(guān),而與它們的絕對位置無關(guān)。地質(zhì)工作中,采樣點(diǎn)以及估計點(diǎn)的位置坐標(biāo)是已知的,因而可得到采樣點(diǎn)之間、采樣點(diǎn)與估計點(diǎn)之間的距離相關(guān)量。距離相關(guān)量c(ui,uj)表示采樣點(diǎn) i和 j的距離相關(guān)關(guān)系, c(ui,z)表示采樣點(diǎn) i和估計點(diǎn)z的距離相關(guān)關(guān)系?？肆⒏袼惴ㄖ惺褂镁嚯x相關(guān)量 c(ui,uj)和c(ui,z)來代替 (2)式中的未知部分,得到:

從公式 (4)可以看出,通過采樣點(diǎn)之間、采樣點(diǎn)與估計點(diǎn)之間的距離相關(guān)量可以得到各個采樣點(diǎn)對估計點(diǎn)貢獻(xiàn)的權(quán)重wi。這樣可以通過 (1)式求得采樣區(qū)內(nèi)任一估計點(diǎn)的估計值。

1.3 通過變差函數(shù)求距離相關(guān)量

克立格算法中使用變差函數(shù)求得距離相關(guān)量。變差函數(shù)定義為:假設(shè)空間點(diǎn)只在一維 x軸上變化,就把區(qū)域化變化量 z(x)在 x,x+h兩點(diǎn)處的值之差的方差之半定義為z(x)在 x方向上的變差函數(shù),記為 r(x,h)。

式中 E為數(shù)學(xué)期望,可以看出,變差函數(shù)通過求某一方向上距離不同的點(diǎn)對所帶信息間差平方的平均,實(shí)際上是要表達(dá)在一個方向上距某點(diǎn)一定距離位置上的信息與這個點(diǎn)上信息之間的關(guān)系,其自變量是一個矢量。

2 應(yīng)用實(shí)例

2.1 項(xiàng)目概況

張集煤礦位于淮南礦區(qū)的西部,為新生界松散層所覆蓋的較厚隱伏式煤田,上覆新生界地層厚度為70～550 m,從上至下大致可分為3個含水層(組)、相應(yīng)隔水層 (組)以及底部礫石層,其中底部礫石層直接覆蓋在煤系地層上。煤系地層中含煤多層,其傾角為10～20°,煤層穩(wěn)定,目前主采煤層厚為4 m左右,且構(gòu)造相對簡單。原設(shè)計單位在礦井設(shè)計時,留設(shè)了80 m垂高的防水煤巖柱,積壓了淺部煤炭資源。通過對該礦區(qū)地質(zhì)條件進(jìn)一步的認(rèn)識,開展了“張集煤礦13-1#、11-1#煤層合理煤巖柱留設(shè)可行性研究”工作,研究影響上提工作面的新生界下部松散層巖性結(jié)構(gòu)、厚度、富水性、基巖風(fēng)化帶特征、覆巖結(jié)構(gòu)性質(zhì)與類型、回采時“二帶”高度的破壞規(guī)律,從而為張集煤礦1100、1300工作面提高上限的安全開采提供基礎(chǔ)性論證。防水煤巖柱留設(shè)高度是指煤層頂板至基巖風(fēng)化帶頂面 (即松散層沉積物沿沉積前的古地形露頭面)之間的垂直距離,因此對各巖層頂?shù)装鍢?biāo)高精確控制至關(guān)重要。巖層面標(biāo)高控制是指利用現(xiàn)有的巖層面控制點(diǎn)標(biāo)高來描述巖層展布狀況。

圖1 鉆孔分布、基巖面標(biāo)高等值線—矢量綜合圖/m

2.2 克立格插值法應(yīng)用

將克立格插值法在張集煤礦 13-1#(11-1#)煤層的1300(1100)工作面基巖面標(biāo)高控制中進(jìn)行了應(yīng)用。在詳勘設(shè)計施工中,上提工作面及其附近11個鉆孔位置、及由這11個鉆孔數(shù)據(jù)經(jīng)克立格插值計算得到基巖面標(biāo)高等值線—矢量圖如圖1所示。由圖1可以看出,基巖面標(biāo)高各向異性顯著,南北為主梯度方向 (傾向)。擬工作面附近等值線彎曲變形與整體趨勢不協(xié)調(diào),說明該局部范圍基巖面標(biāo)高差異較大,同時擬工作面內(nèi)控制鉆孔太少,對標(biāo)高的估計可能存在較大誤差,因此后期又補(bǔ)勘4個鉆孔,即抽水1、抽水2、覆巖1、覆巖2。所加4個孔的位置及加密鉆孔后等值線圖如圖2所示,覆巖孔2布置的位置是為了加強(qiáng)主梯度方向的控制。

鉆孔加密后基巖面標(biāo)高等值線較之加密前擬工作面附近彎曲更為明顯,說明鉆孔加密前對該局部范圍基巖面標(biāo)高估計存在較大誤差,驗(yàn)證了鉆孔加密控制的必要性。

圖2 鉆孔分布、加密控制后基巖面標(biāo)高等值線圖/m

在鉆孔未加密的條件下由克立格插值法計算出的4個補(bǔ)勘鉆孔相應(yīng)位置的基巖面標(biāo)高與實(shí)際測量取得的標(biāo)高之間的誤差如表1所示。

表1 補(bǔ)勘鉆孔位置實(shí)測基巖面標(biāo)高與插值求得基巖面標(biāo)高比較

由表1可知:覆巖1、覆巖2所在位置的估計標(biāo)高與實(shí)測標(biāo)高存在較大誤差,其中覆巖2所在的位置基巖面標(biāo)高估計誤差為19.30 m,已經(jīng)超過《建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開采規(guī)程》定義的保護(hù)層厚度12.0 m(H保=3× 4.0 m=12.0 m,其中3為系數(shù),4 m為采高),這說明在鉆孔未加密控制情況下的任何防水煤巖柱留設(shè)設(shè)計都沒有實(shí)際意義。在基巖面露頭附近因構(gòu)造、風(fēng)化作用強(qiáng)烈,地形起伏變化較大,在鉆孔控制不夠或布置不合理的情況下存在較大誤差是經(jīng)常遇見的。由克立格插值法原理及變差函數(shù)性質(zhì)結(jié)合其在標(biāo)高控制中應(yīng)用分析可知,控制誤差因素有以下兩個方面。

(1)區(qū)域構(gòu)造情況,古地貌風(fēng)化剝蝕程度。這是因?yàn)閿?shù)據(jù)的空間變異總是可以區(qū)分為受大尺度區(qū)域因素控制的規(guī)律變化或方向性變化和受小尺度局部因素制約的隨機(jī)變化。

(2)待估值點(diǎn)周邊控制點(diǎn)的疏密程度及布點(diǎn)方式。克立格插值法是通過統(tǒng)計采樣點(diǎn)之間距離、自動設(shè)置基本步長 (滯后距),進(jìn)行變差函數(shù)模型擬合的。同時克立格插值法是加權(quán)平均,權(quán)值大小由點(diǎn)對之間距離決定的。

分析控制誤差因素為施工中進(jìn)行合理布孔提供理論依據(jù)。在設(shè)計布孔方案時,認(rèn)真分析鉆孔所處區(qū)域構(gòu)造、地層產(chǎn)狀,盡可能使鉆孔布置區(qū)地質(zhì)環(huán)境一致,以減小誤差;在地層走向與傾向上合理控制鉆孔密度,在消減誤差因素影響的同時達(dá)到經(jīng)濟(jì)原則。

2.3 通過交叉驗(yàn)證,確定勘探精度

通過以上分析,在理解影響克立格插值法精度的因素基礎(chǔ)上,運(yùn)用交叉驗(yàn)證評價克立格插值法的精度。所謂交叉驗(yàn)證就是假設(shè)研究變量為 z(x),依次刪去采樣點(diǎn) xi(i=1,2,3,…n),其它點(diǎn)的屬性值保持不變,利用剩下的n-1個點(diǎn)的屬性值,插值計算被刪除的采樣點(diǎn)上的屬性值 z＊(x),并對n個插值計算結(jié)果與實(shí)際的結(jié)果作比較,進(jìn)行誤差的統(tǒng)計學(xué)分析。其理論基礎(chǔ)是:在研究區(qū)域滿足二階平穩(wěn)假設(shè)前提下,若控制點(diǎn)能夠從整體控制采區(qū)屬性,剔除個別屬性值,通過克立格插值法計算得出的相應(yīng)估計值應(yīng)接近真實(shí)值?？紤]到工作面是在原有的設(shè)計基礎(chǔ)上進(jìn)行上提的,所以只需選取上提工作面風(fēng)巷附近抽水1、覆巖1、十一1、抽水2及二十65等鉆孔進(jìn)行交叉驗(yàn)證。現(xiàn)采用兩種方案對基巖面標(biāo)高進(jìn)行控制并通過交叉驗(yàn)證分析控制精度,通過比較選擇合理控制方案。

(1)方案1。不考慮誤差產(chǎn)生因素,利用現(xiàn)有的15個鉆孔資料,對基巖面標(biāo)高進(jìn)行控制,并對風(fēng)巷附近各鉆孔標(biāo)高進(jìn)行交叉驗(yàn)證以確定控制精度,見表2。

提高開采上限只需對煤層露頭附近基巖標(biāo)高進(jìn)行局部控制,主要研究由隨機(jī)變化引起的誤差。由距離相關(guān)性可知,相距較遠(yuǎn)的兩點(diǎn)之間幾乎沒有聯(lián)系,卻增加模型擬合的難度,同時弱化了鄰近鉆孔的作用 (尺度效應(yīng))?；谶@種考慮,設(shè)計了標(biāo)高控制方案2,并對風(fēng)巷附近鉆孔進(jìn)行交叉驗(yàn)證。

表2 方案1交叉驗(yàn)證產(chǎn)生的誤差統(tǒng)計

(2)方案2?？紤]誤差產(chǎn)生因素,在現(xiàn)有的15個鉆孔中剔除162、117、二十32、165、112等較遠(yuǎn)鉆孔,利用余下的10個鉆孔資料對基巖面標(biāo)高進(jìn)行控制并對風(fēng)巷附近各鉆孔進(jìn)行交叉驗(yàn)證以確定控制精度,見表3。

表3 方案2交叉驗(yàn)證產(chǎn)生的誤差統(tǒng)計

由表2與表3可知:覆巖1孔的誤差均較大。抽水1、覆巖1、十一1鉆孔在同一走向方向,其中覆巖1孔所處地勢低洼 (自動擬合變差函數(shù)模型,默認(rèn)該3點(diǎn)應(yīng)在同一水平),同時主梯度方向上沒有相對較近的控制鉆孔,這必將導(dǎo)致覆巖1孔估計誤差較大。為評價上述兩種方案哪種對標(biāo)高估計更為合理,在剔除變異較大的覆巖1孔 (不具有統(tǒng)計意義),對兩種方案交叉驗(yàn)證產(chǎn)生的誤差進(jìn)行了統(tǒng)計分析,如表4所示。剔除覆巖1孔并不影響統(tǒng)計結(jié)果的可信度,在方案1中 (或方案2),假設(shè)控制鉆孔中沒有抽水1孔,現(xiàn)在抽水1孔所在位置進(jìn)行取樣驗(yàn)證,其實(shí)質(zhì)就是方案1(或方案2)中對抽水1孔的交叉驗(yàn)證。可依次對十一1、抽水2、二十65孔作同樣處理,這種處理具有交叉驗(yàn)證法的經(jīng)濟(jì)優(yōu)點(diǎn)。同時避免了因個別主控制鉆孔在交叉驗(yàn)證中產(chǎn)生較大的誤差使統(tǒng)計分析失去意義。

由表4看出,方案2各統(tǒng)計量均優(yōu)于方案1,這是由克立格插值法自身結(jié)構(gòu)及鉆孔所處地質(zhì)環(huán)境決定的。本次評價中,采用方案2進(jìn)行基巖面標(biāo)高控制。由表3和表4可知,誤差已經(jīng)能夠滿足生產(chǎn)要求,且此誤差系剔除采樣點(diǎn)后插值計算取得的,大于實(shí)際控制誤差。

表4 方案1與方案2誤差比較 m

3 結(jié)論

通過對克立格插值法原理以及對作為此基礎(chǔ)進(jìn)行估值的變差函數(shù)的介紹,結(jié)合在張集煤礦1300 (1100)工作面基巖面標(biāo)高實(shí)際應(yīng)用情況,分析了影響控制精度因素有3點(diǎn):

(1)區(qū)域構(gòu)造情況,古地貌風(fēng)化剝蝕程度。區(qū)域構(gòu)造越復(fù)雜,古地貌風(fēng)化剝蝕作用越強(qiáng)烈,一定的控制精度需要的控制鉆孔越多。

(2)控制鉆孔的布置方式及疏密程度。一定數(shù)目的控制鉆孔,在主變異方向相對多布置點(diǎn),能夠提高控制精度。

(3)局部控制時,剔除較遠(yuǎn)的相關(guān)性差的鉆孔,能夠提高控制精度 (尺度效應(yīng))。

[1]侯景儒,尹鎮(zhèn)南,李維明等.實(shí)用地質(zhì)統(tǒng)計學(xué)[M].北京:地質(zhì)出版社,1998

[2]張仁鐸.空間變異理論及應(yīng)用 [M].北京:科學(xué)出版社,2005

[3]許光泉,胡友彪,涂敏等.松散含水體下合理安全煤巖柱高度留設(shè)回顧與探討 [J].煤炭科學(xué)技術(shù),2003 (10)

[4]國家煤炭工業(yè)局.建筑物、水體、鐵路及主要井巷煤柱留設(shè)與壓煤開采規(guī)程 [M].北京:煤炭工業(yè)出版社,2000

[5]趙鵬大.定量地學(xué)方法及應(yīng)用 [M].北京:高等教育出版社,2004

[6]萬麗.基于變差函數(shù)的空間異質(zhì)性定量分析 [J].理論新探,2006

(責(zé)任編輯 張毅玲)

Application of Kriging interpolation in production in a coal mine

Su Mingjin1,Xu Guangquan1,Li Peiquan2,Liu Zegong1, Yan Jiaping1,Huang Xiangjing2,Zong Yunfeng2
(1.School of Earth and Environment Engineering,Anhui University of Science and Technology,Huainan,Anhui province 232001,China; 2.Huainan Mining Group,Huainan,Anhui province 232001,China)

TD175.6

B

蘇明金 (1982-),男,安徽理工大學(xué)在讀碩士研究生。

＊科技部支撐項(xiàng)目 (編號:2007BAK28B06)