利用極化率參數(shù)監(jiān)測煤礦滯后突水的可行性

劉志新， 薛國強， 張小楷,4

1 中國礦業(yè)大學(xué)資源與地球科學(xué)學(xué)院， 徐州 221116 2 中國科學(xué)院地質(zhì)與地球物理研究所， 中國科學(xué)院礦產(chǎn)資源研究重點實驗室， 北京 100029 3 中國科學(xué)院大學(xué)地球與行星科學(xué)學(xué)院， 北京 100049 4 貴州省地礦局第二工程勘察院， 遵義 563000

0 引言

在煤炭開采過程中，煤層底板的完整性會遭到一定程度的破壞，在裂隙發(fā)育區(qū)域，往往導(dǎo)致工作面底板突水事故的發(fā)生.特別是采空區(qū)，由于煤炭開采活動便底板產(chǎn)生一定深度的破壞區(qū)域，同時也會使原有隱伏構(gòu)造的裂隙逐步擴大，如果隱伏構(gòu)造與底板下承壓含水層連通，承壓水會沿著隱伏構(gòu)造裂隙涌入工作面采空區(qū)，形成采空區(qū)底板滯后突水災(zāi)害(劉樹才，2008；胡彥博，2020)，如圖1所示.如果工作面底板存在隱伏斷層、陷落柱等構(gòu)造時，煤層開采活動可能誘發(fā)的滯后突水將更加嚴(yán)重.

圖1 采空區(qū)滯后突水示意圖Fig.1 Schematic diagram of lagging water inrush

滯后突水是礦井生產(chǎn)中經(jīng)常遇到的特殊水害類型，因具有一定的隱蔽性、突發(fā)性，對礦井生產(chǎn)的危害極大.近幾年來，煤礦工作面滯后出水發(fā)生較為頻繁，華北石炭-二疊系煤層隨著煤炭資源的高效產(chǎn)出，礦井逐漸向深部延伸，地層中廣泛發(fā)育的巖溶及高壓奧灰?guī)r溶水，均嚴(yán)重威脅安全開采，其中以滯后性突水事故造成的人員傷亡和經(jīng)濟損失最為巨大(虎維岳，2010；虎維岳和尹尚先，2010).工作面滯后突水嚴(yán)重威脅著礦井安全生產(chǎn)，輕者使突水工作面關(guān)閉停產(chǎn)，嚴(yán)重影響礦井正常生產(chǎn)；重者會造成淹井事故，危及財產(chǎn)和工人生命安全.

地球物理方法在與環(huán)境領(lǐng)域勘查中發(fā)揮著重要的作用(Xue et al., 2007；底青云等，2020).通常采用地成電磁方法對煤礦含水結(jié)構(gòu)體進行的探測或監(jiān)測(Xue et al., 2013; Chen et al., 2015, 2019a,b).目前，在采煤工作面回采之前需要進行礦井地球物理勘探，采明工作面煤層內(nèi)部、兩側(cè)及頂、底板采動影范圍內(nèi)是否存在構(gòu)造異常或?qū)Шǖ?，常用的礦井電磁法技術(shù)主要有：礦井瞬變電磁，無線電波透視，礦井直流電法，礦井音頻電透視技術(shù)等(薛國強等,2007,2008,2021; Xue et al., 2019).但是，對于規(guī)模尺度較小的隱伏構(gòu)造，由于探測技術(shù)與資料解釋精度限制，傳統(tǒng)的礦井地球物理方法不能對其進行有效識別，可能導(dǎo)致后期采空區(qū)滯后突水事故的發(fā)生.

在 20 世紀(jì) 90 年代中期就研發(fā)了煤層底板突水征兆監(jiān)測設(shè)備，實現(xiàn)底板位移、鉆孔水壓和彈性波等的監(jiān)測，開展了水害監(jiān)測預(yù)警相關(guān)理論、硬件和軟件等方面的研究(胡耀青等，2000；姜福興等,2008),先后研發(fā)了松散層突水監(jiān)測預(yù)警系統(tǒng)、底板突水災(zāi)害預(yù)警系統(tǒng)等(喬偉等,2020).針對華北型煤田突水類型及特點，提出了監(jiān)測預(yù)警指標(biāo)體系和礦井突水監(jiān)測預(yù)警模式(劉德民,2015；劉德民等，2019)，即根據(jù)監(jiān)測指標(biāo)將水害預(yù)警劃分為物理預(yù)警模式和化學(xué)預(yù)警模式兩種主要類型，并制定了隔水層電阻率參數(shù)動態(tài)全面監(jiān)測和關(guān)鍵部位單點多參數(shù)監(jiān)測的現(xiàn)場實時監(jiān)測方式.另外，以底板“下三帶”理論為基礎(chǔ)，提出了集多頻連續(xù)電法充水水源監(jiān)測、“井-地-孔”聯(lián)合微震采動底板破壞帶監(jiān)測以及監(jiān)測大數(shù)據(jù)智能預(yù)警為一體的煤層底板突水三維監(jiān)測與智能預(yù)警技術(shù)思路(靳德武等，2020).目前，基于電磁法的礦井突水預(yù)測(楊天鴻等,2007)，主要是在已有巷道空間內(nèi)進行電位或電阻率參數(shù)連續(xù)觀測，實現(xiàn)礦井突水過程的實時監(jiān)測，同時可實現(xiàn)對水質(zhì)進行快速檢測識別，判斷水源類型.本論文提出的基于極化率參數(shù)的突水過程監(jiān)測系統(tǒng)，將接收電極或電磁傳感器通過鉆孔布置在采空區(qū)范圍內(nèi)，也可利用下抽巷空間向兩側(cè)鉆孔，采用行列式方式提前布置電極或電磁傳感器，實現(xiàn)對采空區(qū)滯后突水過程的時實動態(tài)監(jiān)測.

論文主要采用三維激發(fā)極化法裝置形式，利用理論分析、數(shù)值模擬和物理模型實驗的方法，研究了采空區(qū)底板下陷落柱構(gòu)造在不同突水階段極化率參數(shù)的變化規(guī)律，論證了利用極化率參數(shù)監(jiān)測煤礦滯后突水的可行性.

1 極化率參數(shù)實時監(jiān)測技術(shù)

1.1 激發(fā)極化效應(yīng)

激發(fā)極化法是利用地殼中巖、礦石不同激電效應(yīng)差異為基礎(chǔ)，外部人為供直流電(時間域)或交流電(頻率域)，通過研究和分析其激電場分布規(guī)律進行找礦和解決地質(zhì)問題的一個電法勘探方法，本文正是基于時間域激電法，在礦井地質(zhì)條件下進行相關(guān)研究分析.

激發(fā)極化效應(yīng)(激電效應(yīng))是在外電流場作用下，巖、礦石與周圍的水溶液產(chǎn)生物理-電化學(xué)反應(yīng)，引起地質(zhì)介質(zhì)中出現(xiàn)電荷分離，產(chǎn)生附加“過電位”的一種物理化學(xué)現(xiàn)象.見圖2，A和B為直流供電點，可通過MN測量總電位差ΔU(t)(一次電位ΔU1與二次電位差ΔU2(t)之和)和二次電位差ΔU2.極化體被極化是一個先快后慢的過程，經(jīng)過Ts秒后趨于飽和狀態(tài).斷電后，ΔU1已不存在，但極化引起的電位差ΔU2(t)先快后慢的逐漸減小，直至消失.

圖2 激發(fā)極化現(xiàn)象示意圖Fig.2 Schematic diagram of induced polarization phenomenon

一般為了得到激發(fā)極化的全過程，需要充電的時間T足夠長(一般幾分鐘)，然后斷電，應(yīng)測量供電時和斷電后電位差變化的全過程.但在實際中，這樣費時費力，成本不劃算，所以一般只取若干個時間點測量的電位差間接求得某些參數(shù)來描述激發(fā)極化特性.常用的參數(shù)是視極化率ηs：

(1)

1.2 極化率參數(shù)正演模擬算法

底板突水激發(fā)極化法的三維有限元正演模擬是利用“等效電阻率法”進行替換求解模擬.第一步，利用軟件數(shù)值模擬得到一次場ΔU1(不含激電效應(yīng))；第二步，根據(jù)“等效電阻率法”替換相應(yīng)地質(zhì)體的電阻率，用軟件再次數(shù)值模擬計算得到總場ΔU(含有激電效應(yīng))；第三步，根據(jù)軟件兩次模擬計算的得到的結(jié)果，用總場減去一次場得到相應(yīng)的二次場ΔU2，然后由視極化率公式得到模擬視極化率.

激發(fā)極化法的視極化率參數(shù)正演計算是利用電阻率來求解的等效電阻率法.由于極化率η是二次場電位差ΔU2與總場電位差ΔU的比值：

(2)

因此，設(shè)大地電阻率為ρ，極化率為η，未產(chǎn)生極化效應(yīng)的一次場電位差為：

(3)

設(shè)K為裝置系數(shù)：

(4)

電阻率為：

(5)

當(dāng)供電時間使巖礦石達(dá)到飽和時，測量電極MN之間的總場電位為：

(6)

其中ρ*為等效電阻率，由式(6)可得：

(7)

根據(jù)極化率η的計算式，可得η、ρ和ρ*三者間的關(guān)系：

(8)

于是有：

(9)

利用ρ*代替ρ來求含激電效應(yīng)總場極化總場電位差為：

(10)

將ρ*代替ρ來求含激電效應(yīng)總場的方法叫“等效電阻率法”，由式(3)、(10)得二次場電位差為：

(11)

兩式相除，得極化率：

(12)

上面通過等效電阻率法計算得到的極化率η即為激發(fā)極化法三維正演模擬的視極化率.

1.3 監(jiān)測系統(tǒng)設(shè)計

采空區(qū)滯后突水觀測系統(tǒng)是采用三維激發(fā)極化法工作方式，觀測系統(tǒng)采用單極-單極工作方式，即供電點A放置于采煤切眼中間位置，接收電極(電磁傳感器)可通過有線或無線方式傳輸采集的電位數(shù)據(jù)，記錄點位置為接收電極正下方.如果主采煤層下方存在瓦斯抽放巷道，可在下方巷道內(nèi)通過鉆孔向兩側(cè)按照行列式方式布置電磁傳感器，具體位置如圖中紅點所示.如果是單一主采煤層，可利用采煤切眼兩端巷道空間向采空區(qū)布置2～4個水平鉆孔，將電磁傳感器通過鉆孔放置在指定位置.本文中要分析兩個水平鉆孔Z1、Z2內(nèi)的電磁傳感器對應(yīng)極化率電性參數(shù)變化規(guī)律，如圖3所示.通過兩個鉆孔中電磁傳感器處極化率參數(shù)的變化規(guī)律及對應(yīng)水平位置，實現(xiàn)對底板下隱伏突水構(gòu)造進行實時動態(tài)監(jiān)測的目的.

圖3 底板監(jiān)測系統(tǒng)示意圖Fig.3 Schematic diagram of monitoring system for bottom plate

2 極化率參數(shù)實時監(jiān)測數(shù)值模擬分析

2.1 模型建立

(1)整個地電模型為長方體，長寬高分別為400 m、300 m、300 m；根據(jù)煤礦采空區(qū)頂?shù)装鍘r層分布，分別設(shè)置三層，底部為底板，厚度190 m，中部為煤層和采空區(qū)，厚度為10 m；上部為煤層頂板，厚度為100 m(圖4)，各參數(shù)設(shè)置見表1.

表1 正演模型參數(shù)表Table 1 Forward model parameters table

(2)陷落柱形態(tài)設(shè)置為圓柱體，半徑為30 m，極化率為8%，其頂部距離采空區(qū)底板高度h(圖4).在模擬計算過程中h值分別為100 m、80 m、60 m、40 m、20 m 5個高度變化模擬水位上升.陷落柱中心點到煤層與采空區(qū)邊界(采煤工作面)橫向距離S，計算是分別取值130 m、100 m、70 m、40 m，模擬含稅水陷落柱靠近采煤工作面；陷落柱中心偏離采空區(qū)中心線的縱向距離為L，取值0 m、40 m、80 m、110 m，模擬縱向含水陷落柱不同位置；以S和L的變化模擬陷落柱的不同位置對電位和極化率影響及其規(guī)律研究.

圖4 模型剖面示意圖Fig.4 Model profile diagram

(3)如圖3所示，供電采用單點電流源A，電流大小10 A.理論上可分析采空區(qū)范圍內(nèi)任一測點的電性參數(shù)變化規(guī)律.考慮工作的可行性，主要研究Z1和Z2兩個鉆孔內(nèi)各測點的參數(shù)變化情況.Z1、Z2是在工作面兩側(cè)巷道端向采空區(qū)內(nèi)布置兩個水平交叉鉆孔.

2.2 底板突水過程監(jiān)測模擬

2.2.1S=100 m，L=0 m，不同高度h引起的視極化率分布

圖5中，A為供電電流源，假設(shè)水位逐漸上升，水位到煤層底板的距離h值分別設(shè)定為100 m、80 m、60 m、40 m、20 m，不同水位情況下的極化率動態(tài)參數(shù)變化分布如圖5a—e所示.從圖5中可以看出，陷落柱異常體投影在縱向和橫向上位于采空區(qū)中心部位，陷落柱異常隨著h降低(水位上升)逐漸出現(xiàn)并越加明顯，極化率異常與陷落柱投影位置靠近.當(dāng)h=100 m時極化率異?；静伙@示，當(dāng)h小于80 m時極化率隨著水位上升，異常響應(yīng)越來越明顯.

圖5 A3供電h變化底板極化率監(jiān)測等值線圖Fig.5 Monitoring isoline map of goaf floor polarizability of A3 power supply with h changed

分別提取不同水位情況下(水位到煤層底板的距離h)的異常極化率極值，獲得極化率極大值隨水位h變化曲線圖(圖6)，h值從100 m減小到80 m時(水位升高20 m)，視極化率參數(shù)值增加了6.0%，曲線斜率為0.268；h值從80 m減小到60 m時，視極化率參數(shù)值增加了11.04%，曲線斜率為0.522；h值從60 m減小到40 m時，視極化率參數(shù)值增加了19.5%，曲線斜率為1.024；h值從40 m減小到20 m時，視極化率參數(shù)值增加了34.1%，曲線斜率為2.142.極大值隨水位高度增加而變化，水位上升，極大值增長率和曲線斜率越來越大，底板極化率變化愈加明顯.

圖6 極化率極大值隨水位h變化曲線圖Fig.6 Curve of maximum polarizability variation with water level h

各模型Z1測線和Z2測線上異常極化率變化圖7所示，曲線總體呈現(xiàn)中間高兩側(cè)低的開口狀，極化率值隨著h減小而升高，在兩條鉆孔方向異常極值位置與陷落柱異常體位置相同，可通過兩條鉆孔測線不同方向的異常極值位置圈定含水陷落柱異常體的位置.

圖7 模型Z1測線和Z2側(cè)線上異常極化率變化Fig.7 Abnormal polarizability changes on Z1 and Z2 lateral lines of the model

2.2.2h=40 m，L=0 m，不同橫向距離S引起的視極化率分布

圖8為A3為供電電流源，異常陷落柱逐漸靠近采煤工作面，橫向距離S分別取130 m、100 m、70 m、40 m變化的動態(tài)極化率參數(shù)變化分布圖.隨著陷落柱逐漸靠近采煤工作面，極化率異常響應(yīng)越加明顯，極化率異常中心逐步向采煤工作面移動，位置偏向異常體在底板投影位置附近偏離供電電極一側(cè)，且隨著橫向距離S變小，極化率異常越來越明顯.

圖8 A供電橫向距離S變化采空區(qū)底板極化率監(jiān)測等值線圖Fig.8 Monitoring isoline map of goaf floor polarizability of A power supply with S changed

提取圖中極化率極大值，以不同橫向距離S為橫坐標(biāo)，對應(yīng)的異常極化率極大值以S，視極化率極大值為縱坐標(biāo)得到圖9.當(dāng)S從130 m到100 m時，極化率極大值增長為6.2%，兩點曲線斜率為0.771；當(dāng)S從100 m到70 m時，極化率極大值增長5.6%，曲線斜率為0.743；當(dāng)S從70 m到40 m時，極化率極大值增長1.3%，曲線斜率為0.186.隨著異常體靠近煤層，當(dāng)S從130 m逐漸靠近到S=40 m時，以30 m為步長，視極化率極大值增長率和曲線斜率由大變小.

圖9 極化率極大值隨橫向距離S變化曲線圖Fig.9 Curve of maximum polarizability variation with transverse distances S

隨著橫向距離增大(從40～130 m)，由于異常體極化效應(yīng)較強，異常體對采空區(qū)底板極化率影響越來越小，導(dǎo)致底板總體極化率減小.當(dāng)陷落柱靠近采空區(qū)與煤層邊界時，靠近供電電源，極化率值相應(yīng)增大.各模型Z1側(cè)線和Z2側(cè)線上異常極化率變化如圖10所示，兩條鉆孔側(cè)線方向曲線規(guī)律相同，極化率曲線總體呈現(xiàn)中間高兩側(cè)低的開口狀，極化率的值隨著S增大而逐漸減小，極大值對應(yīng)鉆孔坐標(biāo)位置遠(yuǎn)離巷道方向移動，位置接近異常體中心.

2.2.3h=40 m，S=100 m，不同縱向距離L引起的視極化率分布

圖11為異常體投影在橫向上處于采空區(qū)底板中部，縱向上距離L取0 m、30 m、60 m、90 m變化的動態(tài)極化率參數(shù)變化分布圖.圖中顯示極化率異常中心區(qū)域位于異常體投影遠(yuǎn)離供電電源一側(cè)，隨著陷落柱縱向距離L增大，異常位置向采空區(qū)上邊界移動，異常體位置和極化率異常位置同步移動，隨著L的增大，異常極值小于6%的幅度在減小.

圖11 A3供電縱向距離L變化采空區(qū)底板極化率監(jiān)測等值線圖Fig.11 Monitoring isoline map of goaf floor polarizability of A3 power supply with L changed

從等值線圖中看不出明顯參數(shù)值變化，提取圖中極化率極大值，以縱向距離L為橫坐標(biāo)，異常極化率極大值為縱坐標(biāo)，得到圖12.當(dāng)L從0 m到30 m時，異常極化率極大值增長-0.12%，曲線斜率為-0.033；當(dāng)L從30 m到60 m時，異常極化率極大值增長-1.0%，曲線斜率為-0.26；當(dāng)L從60 m到90 m時，異常極化率極大值增長-5.1%，曲線斜率為-1.327.極化率異常極大值保持減小，增長率越來越小，總體值得變化幅度很小.

圖12 極化率極大值隨縱向距離L變化曲線圖Fig.12 Curve of maximum polarizability variation with longitunal distances L

A供電隨著縱向距離L增大(從0～90 m)，由于異常體極化率高，極化效應(yīng)影響采空區(qū)底板極化率分布，水位沒有變化的情況下，縱向距離L變化對底板極化率影響較小.所以當(dāng)陷落柱靠近采空區(qū)邊界時，極化率值變化不大.各模型Z1側(cè)線和Z2側(cè)線上異常極化率變化如圖13所示，Z1方向上異常極值位置位于異常體在采空區(qū)中的投影位置遠(yuǎn)離巷道一側(cè)，極大值隨著L增大呈現(xiàn)先增大后減小的趨勢，位置逐漸向右側(cè)移動.Z2方向上極化率曲線隨著L增大而減小，極值位置向著巷道移動.根據(jù)兩條鉆孔方向上極值位置做平行于采空區(qū)上邊界的直線能大致圈定異常體的位置.

圖13 模型Z1測線和Z2側(cè)線上異常電位變化Fig.13 Abnormal potential changes on Z1 and Z2 lateral lines of the models

3 物理模型實驗

3.1 實驗監(jiān)測方案設(shè)計

實驗?zāi)Ｐ烷L寬高大小分別為380 cm×300 cm×180 cm，分為三部分：上部是覆蓋層(黏土)30 cm；中部是弱含水層(細(xì)砂)130 cm，在其內(nèi)部存在導(dǎo)水陷落柱(充填礫石)大小為100 cm×100 cm×130 cm，導(dǎo)水陷落柱位于模型中部；下部為含水層(粗砂)10 cm和底部隔水層(黏土)10 cm，圖14是物理模型的設(shè)計圖，通過從下部含水層和供水管道陷落柱異常體中持續(xù)供水，右側(cè)安裝水位監(jiān)測管監(jiān)測水位變化.

圖14 物理模型實驗平臺設(shè)計圖Fig.14 Physical model experiment design

實驗儀器選用法國IRIS公司生產(chǎn)的SYSCAL電法儀，其中SYSCAL Pro單元是專門針對高效電法勘探所設(shè)計研發(fā)生產(chǎn)的電阻率&IP測量計.SYSCAL Pro單元能夠測量一次場電壓和衰減電壓曲線值，從而可以計算出電阻率和極化率.

電極布設(shè)：根據(jù)觀測系統(tǒng)的設(shè)計方案，本次實驗電極布置方式采用三維激發(fā)極化法觀測布極方式，觀測系統(tǒng)采用單-單極方式，供電電極A位于36號點，另一供電電極B裝置于室外20 m以外.供電電流為1 A.電極在模型上表面以間隔40 cm距離均勻布置，根據(jù)實驗?zāi)Ｐ痛笮≡O(shè)置為7排9列，供63個電極.其中邊界電極與模型邊界距離都為30 cm.在模型的對角位置上布置兩條斜向測線Z1和Z2，主要分析兩條斜向測線Z1和Z2位置上的測點極化率曲線變化規(guī)律，如圖15所示.

實驗步驟如下：

(1)根據(jù)電極布設(shè)方案在模型頂面布設(shè)電極，如圖15所示.

圖15 實驗電極布設(shè)圖Fig.15 Experimental electrode layout

(2)向模型底部含水層充水，從水位監(jiān)測管觀察水位的高度.

(3)通過觀察水位監(jiān)測管，水位從80 cm開始第一次監(jiān)測數(shù)據(jù)采集，上升過程中由于含水層水的滲流，降低水的注速度，保持水位不變直到完成數(shù)據(jù)采集.以水位上升20 cm為一個周期，重復(fù)上述采集過程完成80 cm、100 cm、120 cm、140 cm水位數(shù)據(jù)采集工作.

(4)完成數(shù)據(jù)采集工作之后，結(jié)束數(shù)據(jù)采集工作，整理儀器設(shè)備，進行數(shù)據(jù)反演處理與解釋.

3.2 實驗數(shù)據(jù)分析

首先對實驗數(shù)據(jù)采進行三維成像，可以得出隨著陷落柱構(gòu)造內(nèi)容水位上升，極化率逐漸增大.由于模型頂面由干燥的黏土組成，當(dāng)水位較低的時候，水位對頂面的影響較小，隨著水位上升，模型含水層水面上升，模型整體極化率逐漸升高.為了更清晰反應(yīng)極化率參數(shù)隨水位變化規(guī)律，提取了三維立體數(shù)據(jù)中深度為0.7 m處的水平切片數(shù)據(jù)，得到-0.7 m平面極化率等值線分布情況，如圖16所示.在圖16a是水位80 cm時-0.7 m深度截面極化率分布等值線圖，從圖中可以看出中部陷落柱位置出現(xiàn)一個低極化率異常區(qū)域，左側(cè)出現(xiàn)一個高極化率異常區(qū)域和右側(cè)靠近邊緣出現(xiàn)一個低極化率異常區(qū)域.由于前期實驗?zāi)Ｐ瓦M行了第一次物理模型實驗，排水時間不夠?qū)е履Ｐ腿鹾畬涌紫端扛?，反之陷落柱?nèi)部孔隙直徑較大，殘留的水含量相對含水層低.加之此時模型底板供水水位為80 cm，-0.7 m界面上水位未達(dá)到，因此陷落柱異常體極化率異常顯示為低極化率異常.左側(cè)高極化率異常是由于模型建好之后先進行了污染物模擬實驗所殘留的污染物吸附于細(xì)砂顆粒中引起.右側(cè)設(shè)置了排水管，相對左側(cè)排水效果更好，孔隙水含量低，導(dǎo)致了右側(cè)低極化率異常.

圖16 物理模型-0.7 m平面極化率分布等值線圖Fig.16 Physical model -0.7 m plane isoline map of polarizability distribution

圖16b、c分別為水位100 cm和120 cm時-0.7 m深度截面極化率分布等值線圖，相對于圖16a，圖16b、c陷落柱位置及右側(cè)低極化率異常區(qū)域減小，左側(cè)高極化率區(qū)域增大，等值線密度增大.此時水位到達(dá)截面-0.6 m深度，超過了-0.7 m深度截面，由于陷落柱水位上升，其與四周含水層保持接觸，水由于水壓較高向四周擴散，導(dǎo)致深度為-0.7 m的截面極化率截面顯示為高低極化率異常均存在的情況.圖16d為水位為140 cm時-0.7 m深度截面極化率分布等值線圖，與圖16a、b、c相比，圖16d低極化率異常區(qū)域幾乎不存在，只有右側(cè)存在很小區(qū)域低極化率異常，極化率等值線圖中心陷落柱異常體位置極化率由于水位上升而極化率升高.物理模型實驗表明極化率參數(shù)在監(jiān)測底板陷落柱突水過程的變化規(guī)律與數(shù)值模擬結(jié)果相似，能較好的反應(yīng)底板隱伏構(gòu)造內(nèi)水位上升情況，證明該方法對監(jiān)測采空區(qū)底板下隱伏構(gòu)造滯后突水的有效性.

圖17為深度-0.7 m時，與數(shù)值模擬相同的兩條鉆孔測線Z1和Z2在不同供水水位情況下極化率曲線.從圖中可以看出，在Z1和Z2中心坐標(biāo)處，兩條測線都顯示中部低極化率異常，Z1左側(cè)高極化率異常，Z2右側(cè)的高極化率異常.隨著水位上升，各水位對應(yīng)的極化率曲線極大值和極小值差值逐漸減小.這是由于含水層本省含水，當(dāng)陷落柱從不含水到含水時，極化率升高，同時水的極化率值不會無限升高，所以高極化率異常和低極化率異常區(qū)別越來越不明顯.與數(shù)值模擬鉆孔測線所顯示規(guī)律一致，隨著水位升高，陷落柱異常體極化率隨之升高.

圖17 斜向鉆孔測線Z1和Z2極化率曲線Fig.17 Polarizability curves of oblique borehole survey lines Z1 and Z2

4 結(jié)論

(1) 當(dāng)水位上升，異常體響應(yīng)越來越明顯，極化率參數(shù)的極大值增大，通過極大值點的位置能確定陷落柱異常體在采空區(qū)得平面位置.兩條鉆孔內(nèi)極化率曲線異常極大值位置能圈定異常體位置，同時能反應(yīng)水位的變化情況.當(dāng)利用極化率參數(shù)對底板突水情況進行監(jiān)測的時候，根據(jù)不同的測線所觀測異常疊加的方法來圈定異常體的位置，并依據(jù)異常的值的大小變化監(jiān)測其水位變化.

(2)在物理模型實驗平臺采用三維直流電法觀測方式，以三維極化率立體圖作為對比.物理模型監(jiān)測結(jié)果與數(shù)值模擬結(jié)果具有很高的一致性，結(jié)果表明極化率參數(shù)對采空區(qū)底板突水監(jiān)測技術(shù)的有效行和可行性.