陳亞萍

（蘭州資源環(huán)境職業(yè)技術(shù)大學(xué)，甘肅蘭州 730022）

巖溶地貌又稱喀斯特地貌，在我國(guó)分布比較廣泛，是以水對(duì)碳酸鹽巖、石灰?guī)r、石膏等可溶性巖石的化學(xué)性溶蝕為主，以流水的沖刷、潛蝕為輔而形成的地質(zhì)容貌。由于煤田地下巖溶發(fā)育地區(qū)巖土地質(zhì)結(jié)構(gòu)比較脆弱，再加上人類地質(zhì)活動(dòng)向地下巖溶發(fā)育地區(qū)巖土地質(zhì)結(jié)構(gòu)施加壓力，導(dǎo)致煤田巖溶塌陷發(fā)生。通常情況下煤田巖溶塌陷具有立即性、隨機(jī)性、可預(yù)測(cè)性等特點(diǎn)，巖溶塌陷的發(fā)生會(huì)影響到煤田開采工程進(jìn)度和成本，同時(shí)也會(huì)影響到施工人員人身安全，因此對(duì)煤田巖溶塌陷預(yù)測(cè)是非常有必要的。現(xiàn)行的預(yù)測(cè)方法有基于層次分析法的煤田巖溶塌陷預(yù)測(cè)方法，該方法是以層次分析法為理論依據(jù)，根據(jù)煤田巖溶監(jiān)測(cè)數(shù)據(jù)對(duì)巖溶塌陷發(fā)生概率進(jìn)行分析。由于傳統(tǒng)方法選取的預(yù)測(cè)因子不夠全面，需要對(duì)預(yù)測(cè)因子進(jìn)行逐層分析才能得出最終的預(yù)測(cè)結(jié)果，預(yù)測(cè)過程比較復(fù)雜，導(dǎo)致傳統(tǒng)方法準(zhǔn)確率較低，已經(jīng)無法滿足煤田巖溶塌陷預(yù)測(cè)在精度方面的需求，為此提出基于多點(diǎn)統(tǒng)計(jì)的煤田巖溶塌陷預(yù)測(cè)方法研究。

1 煤田巖溶塌陷形成機(jī)理

由于煤田巖溶發(fā)育地區(qū)地質(zhì)環(huán)境不同，誘發(fā)巖溶塌陷的因素不同，會(huì)使煤田巖溶發(fā)育地區(qū)產(chǎn)生不同的受力情況，從而導(dǎo)致不同程度的巖溶塌陷發(fā)生。但是通常情況下，煤田巖溶塌陷的發(fā)生與生成主要是受到一種主導(dǎo)機(jī)制影響，多種力學(xué)機(jī)制相互耦合的結(jié)果。圖1為煤田巖溶塌陷形成機(jī)理示意圖。

圖1 煤田巖溶塌陷形成機(jī)理示意圖

如上圖所示，煤田巖土覆蓋層下覆含有溶洞或者裂縫的基巖，并且煤田地下水含水層位于該覆蓋層處，煤田地下水地質(zhì)結(jié)構(gòu)與巖土覆蓋層巖溶發(fā)育地質(zhì)結(jié)構(gòu)相互重合，由此構(gòu)成了煤田巖溶塌陷的基本形成條件。在煤田巖土覆蓋層處的溶洞或者裂縫，受到周圍地下水的潛蝕作用影響，溶洞和裂縫面積逐漸增大[1]。

2 基于多點(diǎn)統(tǒng)計(jì)的煤田巖溶塌陷預(yù)測(cè)方法設(shè)計(jì)

2.1 煤田巖溶塌陷預(yù)測(cè)因子選取與分析

對(duì)煤田巖溶塌陷形成機(jī)理分析可以得出，煤田巖溶塌陷影響因素主要包括地下水潛蝕、崩解、真空吸蝕[2-3]，因此選取該三個(gè)因素作為煤田巖溶塌陷預(yù)測(cè)因子。根據(jù)太沙基平衡原理，導(dǎo)致煤田巖溶塌陷橫向潛蝕作用形成的水力梯度計(jì)算如下：

式中，j表示導(dǎo)致煤田巖溶塌陷橫向潛蝕作用形成的水力梯度；G表示煤田巖土覆蓋層土顆粒比重；n表示煤田巖土覆蓋層土體的孔隙度[4]。

由達(dá)西定律可知，地下水縱向潛蝕與地下水流速和水力坡度相關(guān)，其計(jì)算公式如下：

式中，jz表示煤田巖溶縱向潛蝕力；k表示煤田巖土覆蓋層滲透系數(shù)；s表示煤田巖土覆蓋層地下水流速；d表示煤田巖土覆蓋層縱向水力坡度[5]。根據(jù)煤田巖溶覆蓋層地下水原始流速，計(jì)算其煤田巖溶縱向潛蝕動(dòng)能變化量：

式中，ΔE表示煤田巖溶縱向潛蝕動(dòng)能變化量；m表示煤田巖土覆蓋層地下水質(zhì)量；v1表示煤田巖土覆蓋層地下水最大流速；v2表示煤田巖土覆蓋層地下水原始流速[6]。

2.2 煤田巖溶塌陷預(yù)測(cè)因子數(shù)據(jù)采集

對(duì)于煤田巖溶塌陷的預(yù)測(cè)是以上文選取的預(yù)測(cè)因子的真實(shí)數(shù)據(jù)為依據(jù)，此次根據(jù)實(shí)際需求，利用可控源音頻大地電磁法對(duì)預(yù)測(cè)因子數(shù)據(jù)進(jìn)行采集[9]。利用地下煤田巖溶區(qū)域內(nèi)地下水、巖土等不同介質(zhì)的導(dǎo)電率，測(cè)量一次場(chǎng)電位與磁場(chǎng)強(qiáng)度變化，根據(jù)卡格尼阿理論，計(jì)算出煤田巖溶區(qū)域內(nèi)不同介質(zhì)的電阻率，其計(jì)算公式如下：

式中，emin表示電磁信號(hào)最小收發(fā)距離；emax表示磁信號(hào)最大收發(fā)距離；I表示V8多功能電磁儀供電電流；dl表示V8多功能電磁儀供電偶極長(zhǎng)度；β表示V8多功能電磁儀測(cè)線角度；Rmin表示在給定噪聲下V8多功能電磁儀磁探頭能觀測(cè)到的最小電場(chǎng)。根據(jù)上述公式計(jì)算結(jié)果設(shè)定V8多功能電磁儀電磁信號(hào)收發(fā)距離參數(shù)，為后續(xù)基于多點(diǎn)統(tǒng)計(jì)分析煤田巖溶塌陷概率提供數(shù)據(jù)依據(jù)。

2.3 基于多點(diǎn)統(tǒng)計(jì)估算煤田巖溶塌陷概率

以采集到的預(yù)測(cè)因子數(shù)據(jù)作為數(shù)據(jù)依據(jù)，利用多點(diǎn)統(tǒng)計(jì)方法對(duì)預(yù)測(cè)因子進(jìn)行分析，將煤田巖溶塌陷作為多點(diǎn)統(tǒng)計(jì)的“數(shù)據(jù)事件”，公式表示如下：

式中，d(u)表示多點(diǎn)統(tǒng)計(jì)數(shù)據(jù)事件，即煤田巖溶塌陷；i(ua)表示數(shù)據(jù)模板上向量所在位置的狀態(tài)值；SKa表示數(shù)據(jù)模板上的向量屬性；a表示可取的狀態(tài)值數(shù)量。

2.4 預(yù)估煤田巖溶塌陷程度及危險(xiǎn)性

此次結(jié)合實(shí)際需求，對(duì)于描述煤田巖溶塌陷程度設(shè)定了一級(jí)、二級(jí)、三級(jí)、四級(jí)五級(jí)五個(gè)等級(jí)，多點(diǎn)統(tǒng)計(jì)分析結(jié)果在0%～20%，發(fā)生煤田巖溶塌陷可能性非常小，塌陷程度為一級(jí)，所帶來的危險(xiǎn)性也非常低；多點(diǎn)統(tǒng)計(jì)分析結(jié)果在20%～40%，發(fā)生煤田巖溶塌陷可能性比較小，塌陷程度為二級(jí)，所帶來的危險(xiǎn)性也比較低；多點(diǎn)統(tǒng)計(jì)分析結(jié)果在40%～60%。

3 實(shí)驗(yàn)論證分析

以某煤田巖溶為實(shí)驗(yàn)對(duì)象，該煤田巖溶區(qū)域面積為1 643.16m2，溶洞數(shù)量為25個(gè)，裂縫數(shù)量為43條，巖溶巖土結(jié)構(gòu)為蜂窩狀、支架狀，巖土成分為蒙脫石、石英石，熔巖發(fā)育區(qū)域經(jīng)常發(fā)生塌陷，利用設(shè)計(jì)方法與傳統(tǒng)方法對(duì)該煤田巖溶塌陷進(jìn)行預(yù)測(cè)。實(shí)驗(yàn)準(zhǔn)備了五臺(tái)V8多功能電磁儀，布置了10條測(cè)線，100個(gè)測(cè)點(diǎn)，測(cè)線方位為110°，測(cè)線長(zhǎng)度為1 250m，測(cè)點(diǎn)間距為45m，V8多功能電磁儀供電偶極位于測(cè)線以南大約5.5km 的地方，供電偶極極長(zhǎng)為985m，根據(jù)式（8）計(jì)算出電磁信號(hào)最小收發(fā)距離為154m，電磁信號(hào)最大收發(fā)距離975m，結(jié)合實(shí)際情況將電磁信號(hào)收發(fā)間距參數(shù)設(shè)定為462m。將該煤田巖溶發(fā)育區(qū)域分成8個(gè)測(cè)區(qū)，利用V8多功能電磁儀測(cè)量的相關(guān)數(shù)據(jù)總量為364.52GB，具體如表2所示。

根據(jù)表1對(duì)預(yù)測(cè)因子與子因子進(jìn)行計(jì)算和分析，并建立數(shù)據(jù)模板與二維訓(xùn)練圖像，對(duì)每個(gè)測(cè)區(qū)煤田巖溶塌陷概率進(jìn)行計(jì)算，最終確定塌陷程度和危險(xiǎn)性。實(shí)驗(yàn)對(duì)每個(gè)測(cè)區(qū)進(jìn)行巖溶塌陷預(yù)測(cè)，使用OHSF 軟件對(duì)兩種方法預(yù)測(cè)結(jié)果與實(shí)際情況進(jìn)行對(duì)比分析，計(jì)算出兩種方法準(zhǔn)確率，將其作為檢驗(yàn)兩種方法可行性的指標(biāo)，具體如表2所示。

表1 煤田巖溶觀測(cè)數(shù)據(jù)表

表2 兩種方法準(zhǔn)確率對(duì)比（%）

從表2可以看出，設(shè)計(jì)方法預(yù)測(cè)準(zhǔn)確率均在99%以上，預(yù)測(cè)結(jié)果基本與實(shí)際情況相符，遠(yuǎn)高于傳統(tǒng)方法，因此實(shí)驗(yàn)結(jié)果證明了設(shè)計(jì)方法具有較高的預(yù)測(cè)精度，在該方面優(yōu)于傳統(tǒng)方法，更適用于煤田巖溶塌陷預(yù)測(cè)。

4 結(jié)束語

煤田巖溶塌陷預(yù)測(cè)是確保煤田開采工作安全有序進(jìn)行的重要保障，此次針對(duì)傳統(tǒng)方法在實(shí)踐中存在的問題，利用多點(diǎn)統(tǒng)計(jì)方法提出了一個(gè)新的媒體巖溶塌陷預(yù)測(cè)思路，并利用實(shí)驗(yàn)論證了該思路的可行性與可靠性，此次研究對(duì)提高煤田巖溶塌陷預(yù)測(cè)精度具有重要的現(xiàn)實(shí)意義。