YDT88無線電波透視儀在地質異常體探測中的應用

楊高峰，衛金善，楊新亮，竇文武

(山西晉煤集團技術研究院有限責任公司，山西 晉城 048006)

0 引言

受多期次構造運動疊加影響，華北石炭二疊系煤田地質條件復雜。煤礦回采工作面常因未能提前查清陷落柱、斷層等地質異常體而導致減產、停產、甚至重新布置工作面，嚴重制約礦井安全高效生產[1-2]。當前，用于探測回采工作面地質異常體的物探技術主要有地面三維地震法、井下無線電波透視法和槽波地震法等，其中尤以無線電波透視法最為廣泛使用[3-5]。前人通過無線電波透視技術對回采工作面內陷落柱、隱伏斷層、沖刷帶、煤層破碎帶等地質異常體進行了探測，有效指導了安全回采[6-8]。隨著生產工藝改進，回采工作面設計寬度逐漸增大，如此則導致無線電波透視施工時電磁波不能有效穿透煤層，致使物探結果不能準確反映地質異常體發育特征。低頻率、大透距無線電波透視儀的研發與應用成為解決這一問題的關鍵。有學者通過優化探測工藝、降低發射頻率、改變接收天線大小形狀等方式對大面寬回采工作面進行探測[9-11]，一定程度上提高了探測精度。

以山西晉煤集團某礦4308回采工作面為研究對象，通過YDT88無線電波透視儀對其內部地質異常體發育情況進行探測，并結合實際回采驗證情況，以期明確該儀器在大面寬回采工作面地質異常體探測中的可行性，并為晉城礦區其他回采工作面地質異常體探測提供一定借鑒。

1 回采工作面探測

1.1 工作面概況

4308回采工作面由43081巷、43082巷及切眼圈定，走向長1 520 m，傾向寬300 m，回采煤層為二疊系山西組3號無煙煤，煤層結構簡單，煤厚4.8～5.6 m。工作面整體為一傾向北西的單斜構造，受東部區域構造晉獲斷裂帶控制，該工作面地質異常體發育。圈定過程中兩順槽曾揭露斷距1 m以下小型斷層5條、斷距6 m以上斷層1條、43081巷曾揭露1個陷落柱，強行割矸距離為55 m。由此可見，為保證該工作面安全高效回采，須對其內部地質異常體發育情況進行精細化探測與分析。

1.2 探測方案

無線電波透視技術原理：電磁波在地下巖層中傳播時，由于巖層電性不同，它們對電磁波能量吸收有一定差異，電阻率低的巖層具有較大的吸收作用。另外地質異常體的界面會對電磁波產生折射、反射等作用，致使電磁波能量衰減和損耗。因此，如果在電磁波穿越煤層的途徑中，遇到與煤層電性不同的陷落柱、斷層等地質異常體，電磁波能量就會被其吸收或完全屏蔽，信號顯著減弱，形成透視異常。通過對透視異常進行分析，即可獲得地質異常體發育特征[12]。

探測頻率：電磁波在煤層中的透射距離與其頻率有關，頻率高，透射距離小，頻率低，則透視距離大。YDT88無線電波透視儀內設88 K、158 K、365 K、965 K這4種工作頻率，透距0～350 m。為了突出“透視異常”，考慮工作面規模和產狀，采用158 kHz 作為電磁波激勵頻率。

測點布置：如圖1所示，本次探測采用一發一收式的定點法。首先于43081巷布置發射天線，43082巷布置接收天線，由切眼向巷口方向依次發射和接收。通常發射測點間距為50 m，接收測點間距為10 m，為實現對地質異常體的精細探測，本次設計發射測點間距為40 m，接收測點間距為5 m。單條順槽發射測點37個，各發射測點均對應17個接收測點，單次發射接收合計獲取629個場強數據。43081巷完成發射工作后，調換發射天線和接收天線，在43082巷發射電磁波，43081巷進行數據采集，兩順槽共獲取1 258個實測場強值。

圖1 4308工作面無線電波透視射線示意

干擾排除：電力干擾直接影響電磁波的傳播，進而影響實測場強數據的真實性。因此，根據發射接收測點數量和順槽長度，結合以往施工經驗，施工時需確保工作面周邊停電至少180 min。此外，施工選在掘進配套設備全部退出而回采設備未安裝期間進行，盡量減少鐵器對電磁波的干擾。

2 數據分析及驗證

2.1 數據分析

工作面場強吸收系數計算：層析成像技術是醫學CT的延伸，在無線電波透視數據處理中應用廣泛。它基于工作面不同介質對電磁波的吸收差異，通過數學算法反演各介質單元吸收系數，最終以圖像形式直觀展現各種地質體的發育特征[13]。

根據電磁波理論，耗散介質中的電磁波發射的場強H0與接收到的場強HP存在下述關系

(1)

式中，β—吸收系數；r—發射點與接收點之間的距離；f(θ)—方向性因子。β與介質對電磁波的吸收能力有關。一般來說，β值越小，介質對電磁波的吸收越弱；反之，β值越大，介質對電磁波的吸收就越強。

由(1)式可得

(2)

式中，dr—發射點與接收點之間電磁波射線上的射線微元；Q—在接收點觀測到的電磁波場強變化值。對(2)式進行離散化，得到第i個觀測點的線性方程

(3)

式中，βj—反映介質中第j個網格的電磁波波場吸收系數；rij—與第i個觀測點對應的第i條電磁波射線在第j個網格單元中的射線長度；Qi—第i個觀測點的觀測值。對于所有的發射點和接收點，可以得到一個線性方程組

RX=Q

(4)

式中，R—射線元組成的m×n維的大型稀疏矩陣；X=(β1，β2，…，βm)T，反映了介質中電磁波波場吸收程度(吸收系數)；Q=(q1，q2，…，qn)T，為接收點電磁波場強變化的觀測值；n—電磁波射線總數，m為網格單元總數。

求解大型稀疏矩陣的線性方程組(4)，即可獲得反映電磁波波場在介質中的吸收程度圖像，進而可以了解介質的物理性質。采用并行迭代的SIRT方法，即當所有投影數據都計算分配完后，才對網格單元圖像函數進行更新。這樣可以消除某些干擾因素，增強抗干擾性。其一般計算公式為

(5)

吸收系數CT成像：在獲取工作面各單元吸收系數值后，借助Surfer軟件插值功能實現工作面吸收系數CT成像。由圖2可知，4308工作面煤層完整區場強吸收系數為0.01～0.03，但同時存在4處高值異常，其吸收系數均大于0.03；初步分析認為，自切眼向巷口方向150～200 m(YC1)、600～750 m(YC2)、950～1 050 m(YC3)及1 300～1 520 m(YC4)處存在地質異常體。

圖2 4308工作面場強吸收系數CT圖

順槽實測場強曲線：YDT88無線電波透視儀在158 kHz頻率工作時，其初始場強為143.74 db，傳播過程中受煤層等地質體吸收，正常區域實測場強值在33～52 db之間。同時，與場強吸收系數CT圖中高值異常區大體相對應位置同樣出現4處低值異常，如圖3所示，其場強值均低于30 db。其中YC1表現為43081巷5號發射點對應的17個接受測點場強值整體較低，而周邊測點場強值均差異不大，近似連續變化。經核實施工現場，該處存在一部未及時撤出的皮帶機頭及若干托輥，據此分析YC1為鐵器干擾引起的假異常，該處不存在地質異常體；YC2和YC3場強值大部分低于13db，但YC2場強值整體較YC3更低，兩異常區內多個接收測點出現零值，曲線形態亦呈陷落柱典型的“V”字形[14]，43081巷在該區域曾連續揭露2個陷落柱，由此推斷2處異常為陷落柱發育響應所致，且YC2處陷落柱規模較YC3處陷落柱要大；YC4場強曲線整體呈階梯型鋸齒狀，為隱伏斷層響應特點，結合43082巷實際揭露斷層特征分析，該異常為斷層發育區域，同時不排除異常區內伴生次級斷層，斷層沿順槽方向影響范圍為1 300～1 520 m。

圖3 順槽實測場強曲線

2.2 驗證情況

4308工作面回采過程中在物探異常區內實際揭露3個陷落柱和2條斷層，如圖4所示。陷落柱X1和X2沿工作面走向集中發育，位于距切眼610～740 m處，二者對應YC2。X1平面近似圓形，直徑70 m。X1平面近似橢圓，軸向NNE，長軸180 m，短軸40 m；陷落柱X3平面呈不規則形狀，位于距切眼970～1 050 m處，對應YC3；斷層F1、F2對應YC4，其沿順槽方向影響范圍大致為1 320～1 500 m，其中斷層F1由43082巷掘進揭露，斷距6 m，F2斷層為其伴生小型斷層，斷距2.3 m；在分析解釋的假異常YC1處煤層發育完整，未發現地質異常體。由此可見，實際揭露情況與數據分析解釋結果有較好的對應關系。

圖4 4308工作面實際揭露地質異常體示意

3 結論

(1)YDT88無線電波透視儀在158 kHz頻率下工作時透射距離大，性能穩定，可通過加密收發測點的方式實現對300 m寬工作面內部地質異常體的有效探測。

(2)4308工作面煤層完整區場強吸收系數為0.01～0.03，實測場強值在33～52 db之間。可通過場強吸收系數CT圖和場強曲線對其內部地質異常體發育范圍及性質做出定性判斷，有效指導礦井安全回采。

(3)通過加密收發測點增加數據采集密度，一定程度上可以提高無線電波透視法對地質異常體的探測精度，但當地質異常體在工作面走向集中發育或傾向互為遮擋時，仍難以作出準確分辨，需采取其他手段進一步提高探測精度。