蘇明金許光泉李佩全劉澤功嚴(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)。同時(shí)采用交叉驗(yàn)證與采樣驗(yàn)證相結(jié)合的方法,對(duì)不同控制方案的控制精度進(jìn)行評(píng)價(jià),確定合理的評(píng)價(jià)方案。

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

1 克立格插值法

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

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

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

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

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

空間統(tǒng)計(jì)方法認(rèn)為估計(jì)值與采樣點(diǎn)屬性值的相似程度是通過(guò)點(diǎn)對(duì)的距離相關(guān)量來(lái)度量的。點(diǎn)對(duì)距離相關(guān)量只與采樣點(diǎn)間的相互距離有關(guān),而與它們的絕對(duì)位置無(wú)關(guān)。地質(zhì)工作中,采樣點(diǎn)以及估計(jì)點(diǎn)的位置坐標(biāo)是已知的,因而可得到采樣點(diǎn)之間、采樣點(diǎn)與估計(jì)點(diǎn)之間的距離相關(guān)量。距離相關(guān)量c(ui,uj)表示采樣點(diǎn) i和 j的距離相關(guān)關(guān)系, c(ui,z)表示采樣點(diǎn) i和估計(jì)點(diǎn)z的距離相關(guān)關(guān)系。克立格算法中使用距離相關(guān)量 c(ui,uj)和c(ui,z)來(lái)代替 (2)式中的未知部分,得到:

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

1.3 通過(guò)變差函數(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ù)通過(guò)求某一方向上距離不同的點(diǎn)對(duì)所帶信息間差平方的平均,實(shí)際上是要表達(dá)在一個(gè)方向上距某點(diǎn)一定距離位置上的信息與這個(gè)點(diǎn)上信息之間的關(guān)系,其自變量是一個(gè)矢量。

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

2.1 項(xiàng)目概況

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

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

3 結(jié)論

通過(guò)對(duì)克立格插值法原理以及對(duì)作為此基礎(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ù)目的控制鉆孔,在主變異方向相對(duì)多布置點(diǎn),能夠提高控制精度。

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

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

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

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

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

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

[6]萬(wàn)麗.基于變差函數(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)目 (編號(hào):2007BAK28B06)