杜賢軍 牛浩 李明陽(yáng) 韓冬

摘要：為了探明多煤層采空區(qū)分布，利用地震波頻率諧振勘探法及可控源音頻大地電磁測(cè)深法進(jìn)行綜合勘查。本文闡述了地震波頻率諧振勘探法及可控源音頻大地電磁法的方法原理以及采空區(qū)的地球物理特征，結(jié)合實(shí)際案例，從波阻抗斷面圖及視電阻率斷面圖，分析了采空區(qū)的異常特征并圈定采空區(qū)范圍，后經(jīng)鉆探驗(yàn)證，推斷結(jié)果與實(shí)際情況基本吻合，并對(duì)兩種勘探方法的優(yōu)缺點(diǎn)進(jìn)行了分析。本次研究將地震波頻率諧振勘探技術(shù)應(yīng)用于多煤層采空區(qū)的探測(cè)取得了良好的效果，對(duì)多煤層采空區(qū)的探測(cè)有很重要的借鑒和指導(dǎo)意義。

關(guān)鍵詞：采空區(qū)；波阻抗；視電阻率；地震波頻率諧振勘探;可控源音頻大地電磁法

中圖分類(lèi)號(hào)：P631????文獻(xiàn)標(biāo)識(shí)碼：A????doi：10.12128/j.issn.1672-6979.2024.02.005

引文格式：杜賢軍，牛浩，李明陽(yáng)，等.綜合物探對(duì)多煤層采空區(qū)綜合勘查對(duì)比研究［J］.山東國(guó)土資源，2024，40（2）：2732. DU Xianjun， NIU Hao， LI Mingyang， et al. Comprehensive Geophysical Exploration and Comparative Study of Multi-coal Strata Goaf Area［J］.Shandong Land and Resources，2024，40（2）：2732.

0?引言

隨著我國(guó)工業(yè)進(jìn)程的不斷推進(jìn)，對(duì)于能源的需求量不斷加大。煤炭是我國(guó)今后一段時(shí)間的主要能源，但我國(guó)煤炭出產(chǎn)地煤層大多以多層結(jié)構(gòu)賦存。我國(guó)很多煤礦由于歷史原因建礦初期采煤方法不正規(guī)，為了出煤快，早見(jiàn)效益，首先開(kāi)采的是主要大巷附近的煤層，而且沒(méi)有遵循正規(guī)采煤順序，形成多層采空區(qū)，為煤礦后期生產(chǎn)帶來(lái)不少安全隱患。隨著時(shí)代的發(fā)展，煤礦對(duì)安全生產(chǎn)意識(shí)的重視，查明歷史采空區(qū)的范圍及富水性尤為重要。對(duì)于多層采空區(qū)勘查是當(dāng)今物探勘查研究的熱點(diǎn)，常規(guī)電法及電磁法橫、縱分辨低且易受地表干擾，不易對(duì)多層采空區(qū)進(jìn)行解譯。而地震波頻率諧振能夠同時(shí)對(duì)多層采空區(qū)具有良好的分辨能力。本次采用地震波頻率諧振勘探技術(shù)（SFRT技術(shù)）和可控源音頻大地電磁法（CSAMT）進(jìn)行勘探，分析研究勘查區(qū)內(nèi)采空區(qū)的地球物理特征，進(jìn)而查明其空間分布情況，分析兩種方法在多層采空區(qū)勘查的優(yōu)缺點(diǎn)，為多層采空區(qū)勘查提供一種方案。

1?方法原理

1.1?SFRT勘探原理

地震波頻率諧振勘探技術(shù)（SFRT技術(shù)）［112］是近幾年來(lái)在傳統(tǒng)地震勘探基礎(chǔ)上發(fā)展起來(lái)的一種新的地球物理勘探方法，可以進(jìn)行無(wú)源方式施工，具有勘探深度大，分辨率高、施工相對(duì)方便的特點(diǎn)。其基本原理是根據(jù)地下不同地層及地質(zhì)體具有各自不同的頻率諧振的特性進(jìn)行勘探，通過(guò)提取地層的波阻抗參數(shù)（波速和密度），來(lái)分辨目標(biāo)地質(zhì)體的勘探方法。地下地層及地質(zhì)體均具有其固有的頻率，其頻率特性由其形狀、組分以及屬性（波速和密度）不同而不同，自然界賦予了不同物體固有頻率。物體會(huì)對(duì)不同的振動(dòng)，做出不同的對(duì)應(yīng)的反應(yīng)，當(dāng)物體的自身固有頻率與其受到的振動(dòng)頻率相一致時(shí)，物體將發(fā)生類(lèi)似于共振現(xiàn)象，放大振動(dòng)的幅度，此時(shí)，可以說(shuō)物體產(chǎn)生了共振（諧振），這些響應(yīng)中包含了地下介質(zhì)的信息。通過(guò)分析這些響應(yīng)中的特征頻率、幅度和相位等參數(shù)，可以建立地下模型，進(jìn)而推測(cè)其物性參數(shù)、巖層分布、儲(chǔ)層性質(zhì)等地質(zhì)信息。地震波頻率諧振勘探技術(shù)主要用于礦產(chǎn)勘查、地質(zhì)災(zāi)害勘查、油氣資源勘查、井壓裂勘查以及采空區(qū)勘查中，成為地質(zhì)勘查中一種主要方法。

1.2?CSAMT法原理

CSAMT法是在AMT法及MT法基礎(chǔ)上改進(jìn)的一種人工控制場(chǎng)源的電磁測(cè)深法。該方法通過(guò)改變發(fā)射端的頻率來(lái)達(dá)到探測(cè)不同深度的目的。該方法具有生產(chǎn)效率高，抗干擾能力強(qiáng)，勘探深度可達(dá)到2km，水平方向分辨率高且不受高阻層屏蔽等優(yōu)點(diǎn)。近年來(lái)，廣泛用于地?zé)峥辈榧安煽諈^(qū)勘查中［1320］。

2?勘查區(qū)地質(zhì)及地球物理特征

2.1?地質(zhì)概況

勘查區(qū)位于構(gòu)造侵蝕堆積地形的Ⅱ級(jí)超河漫灘階地之上，地勢(shì)較為平坦，地面標(biāo)高485～491m。區(qū)內(nèi)以往地質(zhì)勘查成果鉆孔揭露，區(qū)內(nèi)地層由老至新主要有中生界下白堊統(tǒng)九佛堂組、阜新組及新生界第四系。九佛堂組主要巖性：底部為一層厚層安山質(zhì)礫巖夾砂巖、頁(yè)巖，下部為薄層凝灰質(zhì)粉砂巖、紙片狀頁(yè)巖、凝灰質(zhì)頁(yè)巖，中、上部為砂巖、粉砂巖與頁(yè)巖、粉砂質(zhì)頁(yè)巖、泥巖不等厚互層；阜新組（為區(qū)內(nèi)主要含煤地層）主要巖性為煤層及細(xì)砂巖；第四系主要巖性為礫石層及粘土，厚0～80 m（圖1）。

勘查區(qū)含有6個(gè)煤組11個(gè)可采煤層，即42、4、52、5、62、61A、61、6、60、71及7煤層本次主要針對(duì)4、52、5、62煤層進(jìn)行勘查。

勘查區(qū)煤礦初期由于采煤方法不正規(guī)，采用的是房柱式采煤方式，開(kāi)采的是主要大巷附近的煤層，而且沒(méi)有遵循由上而下的正規(guī)采煤順序，主要開(kāi)采4煤、52煤和5煤煤炭資源。

勘查區(qū)地層呈單斜構(gòu)造，地層走向SN、傾向E、傾角8°～18°，一般為12°，地層平緩?？辈閰^(qū)外部有3條斷層，勘查區(qū)內(nèi)斷裂及褶曲構(gòu)造不發(fā)育。

2.2?地球物理特征

2.2.1?波阻抗特征

黏土的波阻抗一般小于1 000 kg/（m2·s）；玄武巖、礫石層的波阻抗一般在1 800～4 000 kg/（m2·s）；砂巖、泥巖、粉砂巖的波阻抗一般在1800～4 000 kg/（m2·s）之間；煤層、煤層采空區(qū)的波阻抗范圍為500～2000kg/（m2·s），與圍巖存在明顯的波阻抗差異。

煤層采空區(qū)在地震上的地球物理特征主要取決于采空區(qū)與圍巖之間存在波阻抗差異。當(dāng)采空區(qū)形成后，頂板的巖層出現(xiàn)跨落、破碎并伴隨裂隙產(chǎn)生，且采空冒落區(qū)域通常充填碎石、空氣和水，造成與周?chē)膸r層密度和速度發(fā)生變化，與原有巖層的密度和速度存在明顯的差異。由于斷裂帶和彎曲帶巖石上的應(yīng)力的變化，這兩個(gè)區(qū)域的巖石大多是松散的，其波阻抗在一定程度上降低。煤系地層與軟弱巖石交錯(cuò)，與正常完整巖石相比，呈現(xiàn)低密度及低速度異常，其波阻抗異常相對(duì)明顯［16］。

2.2.2?電性特征

第四系：巖性為礫石層及黏土，電阻率一般在50～500Ω·m;阜新組煤層及細(xì)砂巖，電阻率一般在100～1 000Ω·m;九佛堂組泥巖，電阻率一般在500～2000Ω·m。

當(dāng)煤層采動(dòng)形成采空區(qū)后，改變了圍巖原有的電性特征，電阻率值與非采空區(qū)有明顯的差別，當(dāng)煤層采空區(qū)充水時(shí)形成低阻反映，從而為電法勘探尋找采空區(qū)提供了地球物理依據(jù)。

3?數(shù)據(jù)采集

本次物探采用CSAMT法和SFRT技術(shù)兩種方法。地震波頻率諧振勘探進(jìn)行面積性勘探，網(wǎng)度為40 m×20 m，布置測(cè)線16條；CSAMT法主要用于采空區(qū)富水性進(jìn)行評(píng)價(jià)，網(wǎng)度為80 m×20 m，布置測(cè)線8條。其中L1、L3、L5、L7、L9、L11、L13、L15線位置同時(shí)施工了地震波頻率諧振勘探和CSAMT法施工。

SFRT技術(shù)采用北京派特森FRTIMGⅠ型采集站，本采集站內(nèi)置頻帶寬度為0.05～150 Hz的檢波器；施工前對(duì)道間幅度進(jìn)行一致性試驗(yàn)，保證各道間幅度一致性及道間相位符合規(guī)范要求［16］，點(diǎn)距20 m，采集時(shí)長(zhǎng)為18 min。

CSAMT法勘探采用儀器V8多功能電法工作站，野外采用赤道偶極裝置。工作頻率1～9 600 Hz，共選擇了44個(gè)避開(kāi)50Hz工頻及其倍頻的頻點(diǎn)，并且頻率點(diǎn)分布均勻，每個(gè)頻點(diǎn)采集時(shí)間為60 s，收發(fā)距r為8～9 km，發(fā)射偶極距AB為1 600 m；高頻最小發(fā)射電流為9 A，低頻最大發(fā)射電流為18 A；接收端偶極距為20 m，點(diǎn)距20 m。

4?資料處理

4.1?SFRT數(shù)據(jù)處理

本次地震波頻率諧振主要經(jīng)過(guò)以下處理步驟：

（1）波場(chǎng)整形及波場(chǎng)一致性校正：對(duì)采集的原始數(shù)據(jù)進(jìn)行濾波處理。

（2）波場(chǎng)疊加及波場(chǎng)平衡：主要通過(guò)傅里葉變換將時(shí)間序列數(shù)據(jù)變換為頻率域數(shù)據(jù)并進(jìn)行多次疊加。

（3）波場(chǎng)頻率分析及波場(chǎng)分離：對(duì)多次疊加的各分量頻率數(shù)據(jù)進(jìn)行噪音處理及磨光處理。

（4）深度與波阻抗率反演：進(jìn)行諧振遴選和深度分析，最后獲得測(cè)線的波阻抗斷面圖。

4.2?CSAMT法數(shù)據(jù)處理

本次CSAMT法數(shù)據(jù)處理包括以下幾個(gè)步驟：

（1）對(duì)觀測(cè)的數(shù)據(jù)進(jìn)行預(yù)處理（包括數(shù)據(jù)解編、信號(hào)的回放檢查等）。

（2）將各場(chǎng)量的時(shí)間域信號(hào)轉(zhuǎn)換為頻率域測(cè)深曲線。

（3）張量阻抗的性質(zhì)及計(jì)算。

（4）進(jìn)行場(chǎng)源校正以及靜態(tài)位移校正，對(duì)得到的數(shù)據(jù)進(jìn)行二維反演生成視電阻率斷面圖。

5?資料解釋

5.1?SFRT資料解釋

圖2b為地震波頻率諧振視波阻抗成果圖。在淺部，波阻抗小于1000kg/（m2·s），為第四系黏土的反映；在標(biāo)高450m以淺，存在一層波阻抗高值條帶，推斷為礫石層波阻抗的反映；在標(biāo)高430～440 m，波阻抗突變帶，且波阻抗為相對(duì)低值，結(jié)合已知資料為第四系底界面的反映。在標(biāo)高430 m以深，波阻抗等值線呈層狀分布，與實(shí)際地層反映一致，呈單斜特征，整體地層傾角約為13°。圖中存在數(shù)條高低相間的低波阻抗特征，為煤層低速、低密度的反映，結(jié)合已知資料，推斷分別為4煤層、52煤層、5煤層、62煤層的反映，煤層之間的相對(duì)高波阻抗區(qū)域主要為砂巖、泥巖的反映。其中在標(biāo)高220～280 m，樁號(hào)71 800～71 970 m位置，5煤層賦存位置，該異常區(qū)波阻抗值在1 200 kg/（m2·s）左右，呈現(xiàn)相對(duì)低值，推斷該處煤層采空且塌陷，且該塌陷發(fā)育的裂隙影響至到52煤層；樁號(hào)71 860～72 140 m，52煤層賦存位置，波阻抗呈現(xiàn)相對(duì)低值，推斷為52煤層采空區(qū)的影響；另在樁號(hào)71 850～72 140 m，4煤層賦存位置，存在相對(duì)低波阻抗區(qū)域，推斷為4煤層采空區(qū)的影響范圍。

5.2?CSAMT資料解釋

在圖2a CSAMT勘查斷面圖中，視電阻率整體從上往下呈現(xiàn)由低逐漸變高的趨勢(shì)。整體視電阻率等值線呈層狀分布，呈單斜特征，與實(shí)際地層反映一致，傾角約13°。標(biāo)高450 m以淺，視電阻率小于50 Ω·m，呈現(xiàn)相對(duì)低阻，為第四系的電性反映。在標(biāo)高450 m以深，視電阻率呈現(xiàn)相對(duì)高阻，為中生界下白堊統(tǒng)九佛堂組及阜新組的電性反映。在樁號(hào)71 800～72 140 m，標(biāo)高200～350 m，視電阻率等值線梯度變化較大變化，推斷為煤層采空區(qū)的影響范圍，但該區(qū)域從上至下賦存4煤層、52煤層、5煤層，無(wú)法具體反映各煤層采空區(qū)的具體位置。

5.3?分析對(duì)比及成果驗(yàn)證

為驗(yàn)證采空區(qū)的真實(shí)情況，在6煤層掘進(jìn)巷道（物探施工區(qū)間，該位置巷道還未開(kāi)拓）沿上山樁號(hào)71 820～72 100 m向山上布置鉆孔7個(gè)ZK01—ZK07，分別位于樁號(hào)71 820 m、樁號(hào)71 860 m、樁號(hào)719 00 m、樁號(hào)71 940 m、樁號(hào)72 040 m、72 080 m、樁號(hào)72 120 m鉆孔間距不一，所有鉆孔均鉆至4煤層。其中ZK01只揭露了5煤層采空區(qū)，ZK02、ZK03、ZK04鉆孔揭露了5煤層、52煤層、4煤層采空區(qū)，ZK05、ZK06、ZK07鉆孔揭露了52煤層、4煤層采空區(qū)。

揭露采空區(qū)在地震波頻率諧振呈現(xiàn)波阻抗低值閉合圈，在CSAMT法斷面圖中呈現(xiàn)相對(duì)低阻。在CSAMT法斷面圖中，5煤層在樁號(hào)71 970～72 140 m賦存完好，視電阻率與采空區(qū)無(wú)明顯差異，煤層間距較小，CSAMT法在垂向上無(wú)法分辨采空區(qū)；在地震波頻率諧振中，5煤層在樁號(hào)71 970～72 140 m與采空區(qū)位置差異明顯。

7個(gè)鉆孔揭露的采空區(qū)位置及深度，與地震波頻率諧振勘探法解譯相對(duì)吻合較好（圖2c）。而可控源音頻大地電磁法只反映該區(qū)域整體呈現(xiàn)相對(duì)低阻，控制了采空區(qū)的整體邊界，無(wú)法具體確定異常的垂向深度。

6?結(jié)論

（1）CSAMT法具有施工便利，生產(chǎn)效率高，勘探深度最大可達(dá)2 000 m，水平方向分辨率高等特點(diǎn)，可以大致對(duì)地層進(jìn)行層位劃分、圈定采空區(qū)及構(gòu)造的位置。但是由于采空區(qū)厚度小，體積效應(yīng)影響，對(duì)于多煤層采空區(qū)勘查，能夠確定采空區(qū)邊界范圍，垂向分辨能力達(dá)不到劃分多層采空區(qū)，無(wú)法確定煤層具體采掘情況，需結(jié)合勘查區(qū)內(nèi)其他資料或則其他方法補(bǔ)充勘查才能進(jìn)行解釋判斷。

（2）地震波頻率諧振勘探技術(shù)（SFRT法）具有施工靈活，可以單點(diǎn)施工，不受地形限制，抗干擾能力強(qiáng)，勘探深度大，水平方向和垂向都相對(duì)靈敏，尤其對(duì)煤層及煤層采空具有異于圍巖低波阻抗響應(yīng)，通過(guò)地震波頻率諧振勘探可精準(zhǔn)確定煤層及采空埋深、導(dǎo)水裂隙帶發(fā)育高度以及采空水平邊界。

相比可控源音頻大地電磁測(cè)深法，地震波頻率諧振勘探技術(shù)對(duì)多煤礦采空區(qū)勘探分辨率更高。地震波頻率諧振不但能夠解譯多層煤層賦存位置，而且能夠反映多煤層中每一層煤層采空區(qū)的邊界及其影響范圍，但可控源音頻大地電磁測(cè)深法受體積效應(yīng)影響無(wú)法查明多煤層采空區(qū)具體情況。

參考文獻(xiàn)：

［1］?朱聰聰，李海.地震頻率諧振技術(shù)在深埋煤層采空區(qū)探測(cè)中的應(yīng)用［J］.工程地球物理學(xué)報(bào)，2021，18（5）：774779.

［2］?張晗，盧瑋，黃烜，等.綜合地球物理方法在范縣地?zé)峥辈橹袑?duì)比試驗(yàn)研究［J］.地質(zhì)裝備，2021，22（4）：1523.

［3］?張健橋，尹維民，張興洲.綜合物探方法在厚覆蓋煤礦采空區(qū)探測(cè)中的應(yīng)用［J］.地質(zhì)論評(píng)，2021，67（增刊1）：5153.

［4］?張文君.諧波地震在沁水盆地深部煤層氣勘探開(kāi)發(fā)中的初步應(yīng)用［J］.山西煤炭，2021，41（4）：7076.

［5］?徐吉祥，張曉亮，李瀟，等.地表淺部地震勘探方法在城市隱伏活動(dòng)斷裂調(diào)查中的應(yīng)用［J］.城市地質(zhì)，2022，17（1）：7984.

［6］?李紅光，侯錚.基于諧振原理的地震勘探方法在資源勘查中的應(yīng)用［J］.世界有色金屬，2020（9）：159160.

［7］?王亞璐，齊志龍，許國(guó)偉.地震波頻率諧振法在精河縣斷裂構(gòu)造勘查中應(yīng)用及地?zé)崆熬胺治觯跩］.新疆地質(zhì)，2023，41（2）：300304.

［8］?薛愛(ài)民，李達(dá)，宋紅星，等.應(yīng)用振動(dòng)噪聲的頻率諧振效應(yīng)對(duì)地球進(jìn)行成像［J］.地質(zhì)論評(píng)，2021，65（增刊1）：4748.

［9］?李華，楊劍，王橋，等.結(jié)合混合源面波與三分量頻率諧振的淺震勘探在城市地質(zhì)結(jié)構(gòu)精細(xì)劃分中的應(yīng)用［J］.地球物理學(xué)進(jìn)展，2020，35（3）：11491155.

［10］?劉長(zhǎng)會(huì).綜合物探技術(shù)在淺埋煤層采空區(qū)勘查中的應(yīng)用［J］.工程地球物理學(xué)報(bào)，2011，8（1）：5154.

［11］?周治民.煙臺(tái)市棲霞中橋地區(qū)巖溶塌陷地質(zhì)環(huán)境調(diào)查中物探方法應(yīng)用淺析［J］.科技經(jīng)濟(jì)與資源環(huán)境，2019，27（16）：112114.

［12］?郭強(qiáng)，張健橋，耿千頃，等.綜合物探方法在不穩(wěn)定斜坡地質(zhì)災(zāi)害勘察中的應(yīng)用研究［J］.山東國(guó)土資源，2023，39（4）：4549.

［13］?胥溢.臨夏州供水保障工程CSAMT數(shù)據(jù)解譯有效性分析［J］.四川有色金屬，2023（2）：4648.

［14］?杜賢軍，李海波，韓冬.可控源音頻大地電磁測(cè)深法在大傾斜煤層采空區(qū)勘查中的應(yīng)用［J］.山東國(guó)土資源，2023，39（5）：5458.

［15］?韋乖強(qiáng).CSAMT對(duì)含煤地層巖漿巖侵入體探測(cè)的有效性分析［J］.煤炭技術(shù)，2017，36（4）：121123.

［16］?白錦琳.可控源音頻大地電磁法在采空區(qū)探測(cè)中的應(yīng)用［J］.山東國(guó)土資源，2018，34（9）：8286.

［17］?雷達(dá)，孟小紅.復(fù)雜地形條件下的可控源音頻大地電磁測(cè)深數(shù)據(jù)二維反演技術(shù)及應(yīng)用效果［J］.物探與化探，2004，28（4）：323326.

［18］?郭國(guó)強(qiáng)，郭朋，張超，等.可控源法在海陽(yáng)市郭城地區(qū)鑫泰金礦床成礦規(guī)律中的應(yīng)用［J］.山東國(guó)土資源，2014，30（6）：7780.

［19］?馬曉東，劉鵬，陳玉菡.CSAMT法在莘縣馬西林場(chǎng)地?zé)崞詹橹械膽?yīng)用研究［J］.山東國(guó)土資源，2013，29（12）：4446.

［20］?杜賢軍.CSAMT法在費(fèi)縣地?zé)峥辈橹械膽?yīng)用［J］．山東國(guó)土資源，2017，33（10）：7377.

Comprehensive Geophysical Exploration and Comparative Study of Multi-coal Strata Goaf Area

DU Xianjun， NIU Hao， LI Mingyang， HAN Dong

（No.5 Exploration Brigade of Shandong Coalfield Geology Bureau， ?Shandong Mining Disaster Prevention and Emergency Rescue Survey Center， Shandong Ji'nan 250100， China）

Abstract：In order to explore the goaf area and situation of multiple coal seams， the seismic wave frequency resonance exploration method and the controllable source audio geodetic electromagnetic bathymetry method are used for comprehensive exploration. In this paper， the method principle of seismic wave frequency resonance exploration method， controllable source audio geodetic electromagnetic bathymetry method and physical characteristics of the goaf have been explained. Combining with actual cases， the abnormal characteristics of the goaf have been analyzed from the wave impedance section and apparent resistivity section， and the goaf area has been delineated. The inferred results are consistent with actual situation after drilling verification， and the advantages and disadvantages of the two exploration methods have been analyzed. In this study， the application of seismic wave frequency resonance exploration technology to the detection of multi-coal seam goaf has achieved good results. It has important references and guiding significance for the detection of multi-coal seam goaf.

Key words：Goaf; wave impedance; apparent resistivity; seismic wave frequency resonance exploration; controllable source audio geodetic electromagnetic bathymetry

收稿日期：20230628；修訂日期：20231207；編輯：王敏基金項(xiàng)目：山東省煤田地質(zhì)局“地震波頻率諧振在煤田地質(zhì)不良地質(zhì)現(xiàn)象的研究與識(shí)別”，魯煤地科字〔2023〕8號(hào)作者簡(jiǎn)介：杜賢軍（1983—），男，江西撫州人，高級(jí)工程師，主要從事電磁法勘查工作；Email：187315076@QQ.com