基于電偶極子的煤巖破裂電磁輻射研究

王小旗 侯 爭

（1.河南理工大學(xué)計算機(jī)科學(xué)與技術(shù)學(xué)院，河南 焦作454000；2.河南工業(yè)和信息化職業(yè)學(xué)院，河南 焦作454000）

0 引言

作為能源大國，煤礦安全一直是全國工業(yè)安全工作的重點(diǎn)。而煤礦瓦斯災(zāi)害，卻是造成重大傷亡事故、影響生產(chǎn)進(jìn)度、損失經(jīng)濟(jì)利益的主要原因，其中煤與瓦斯突出災(zāi)害的預(yù)防，是其中一個重要環(huán)節(jié)。研究可靠的突出預(yù)測方法，可有效保障礦井安全生產(chǎn)，提高礦井經(jīng)濟(jì)效益。因此，煤巖破裂的預(yù)測研究具有重大的實(shí)際意義及社會意義。

煤巖破裂預(yù)測作為瓦斯災(zāi)害防治綜合措施是最重要的環(huán)節(jié)。目前，煤巖破裂預(yù)測方法主要有以下幾種：聲發(fā)射監(jiān)測方法，電磁輻射預(yù)測方法，以及利用環(huán)境監(jiān)測系統(tǒng)連續(xù)監(jiān)測工作面的瓦斯涌出變化特征，分析瓦斯涌出與突出關(guān)系的預(yù)測方法。其中電磁輻射預(yù)測方法目前正在廣泛研究開展中，但是如何在此基礎(chǔ)上，開發(fā)出從軟件到硬件實(shí)現(xiàn)整合的電磁輻射預(yù)測系統(tǒng)，將使煤巖破裂防治技術(shù)大幅度提升，對改善當(dāng)前煤礦安全狀況具有極其深遠(yuǎn)的意義。

1 煤巖破裂電磁輻射分析

1.1 理論分析

煤巖體同其他固體材料一樣，都是由成千上萬的電子、原子等基本粒子組成，當(dāng)煤巖體受載變形破裂時，電子等帶電粒子變速運(yùn)動就會向外輻射電磁波，這就是電磁輻射現(xiàn)象。當(dāng)煤巖相鄰顆粒之間發(fā)生非均勻變形時，界面處的電平衡遭到破壞，在受拉界面會積累許多自由電荷，而在受壓處則積累了同樣數(shù)量的相反電荷，這就相當(dāng)于電偶極子。因此，可以將煤巖破裂時的輻射源等效為電偶極子。

依據(jù)電磁波的傳播理論，煤巖破裂時輻射的電磁波在分層媒質(zhì)中的傳播，可等效于輻射源在由“巖石層－煤層－巖石層”組成的三層媒質(zhì)模型。研究煤巖破裂時電偶極子場源電磁波的正向傳播規(guī)律可轉(zhuǎn)化為分析該輻射源輻射的電磁波在各媒質(zhì)層中電場和磁場的場量表達(dá)式。由于垂直和水平電偶極子是電磁輻射的基本模型，任何輻射源都可以分解為垂直和水平電偶極子的加和。因此，本文將針對電偶極子場源進(jìn)行研究。

1.2 電磁波輻射場強(qiáng)的計算

由于水平電偶極子與垂直電偶極子場源輻射場強(qiáng)的計算思路基本相同，文章只討論垂直電偶極子場源輻射場強(qiáng)的計算。

圖1 垂直電偶極子情形

為方便理解建立直角坐標(biāo)系。以電偶極子所在點(diǎn)為原點(diǎn)，該電偶極子的長度為dl，其電流為I·ejwt（A）。以垂直指向大地的方向為z軸正方向，如圖1所示。根據(jù)實(shí)際情況，第一、三兩層為巖石層，其電磁參數(shù)分別為μ1、ε1、σ1和μ3、ε3、σ3；第二層為煤層，其電磁參數(shù)為μ2、ε2和σ2。

垂直電偶極子情形下只有z分量，則垂直電偶極子場源情形下中間媒質(zhì)層中各部分的電場和磁場場強(qiáng)表達(dá)式如下所示：

2 系統(tǒng)軟件開發(fā)

關(guān)于三層煤質(zhì)中的場強(qiáng)公式已經(jīng)明確，為了將理論成果更好的應(yīng)用到實(shí)際工程中，可將理論推演通過計算機(jī)編程來實(shí)現(xiàn)。在軟件的選擇上，鑒于Visual c++有很好的穩(wěn)定性和兼容性，我們選用Visual C++6.0開發(fā)平臺作為程序開發(fā)環(huán)境。

無論是電場場強(qiáng)E，還是磁場場強(qiáng)H，都是一個很復(fù)雜的積分算式，含有諸多參數(shù)，但最終都可轉(zhuǎn)化為電導(dǎo)率σ1、σ2、σ3，磁導(dǎo)率μ1、μ2、μ3，媒質(zhì)的介電常數(shù)角頻率ω，電偶極子內(nèi)的電流I，以及相應(yīng)的坐標(biāo)值x、y、z等。在所建三層媒質(zhì)模型中，由于大地媒質(zhì)在通常情況下為非鐵磁性媒質(zhì)，因而此時可設(shè)定典型分層媒質(zhì)模型中各層媒質(zhì)的導(dǎo)磁率相等，并且均等于空氣導(dǎo)磁率，即μ1=μ2=μ3=0。第一、三兩層為巖石層，所以σ1=σ3，σ1=σ3。這樣一、三層媒質(zhì)參數(shù)設(shè)為相等，便于編程實(shí)現(xiàn)。當(dāng)場強(qiáng)計算公式必需的參數(shù)都一定的情況下，調(diào)用復(fù)數(shù)基本乘、除法函數(shù)，讓程序中用到的復(fù)數(shù)運(yùn)算都以實(shí)部和虛部的參數(shù)形式出現(xiàn)，避免了復(fù)雜的復(fù)數(shù)運(yùn)算，再調(diào)用貝塞爾函數(shù)等，得出被積函數(shù)，準(zhǔn)備進(jìn)行積分。此時調(diào)用“積分定限函數(shù)”，確定積分限的情況下再進(jìn)行龍貝格積分，將積分結(jié)果求模值，得出運(yùn)算結(jié)果。具體的程序框圖如圖2所示。在場源參數(shù)都已知的前提下，用戶根據(jù)需要變換坐標(biāo)值的輸入范圍，計算出對應(yīng)的場強(qiáng)數(shù)據(jù)，并可將對應(yīng)場源參數(shù)、坐標(biāo)值、場強(qiáng)數(shù)據(jù)存入相應(yīng)數(shù)據(jù)庫中，作為下一步場源定位查詢的依據(jù)。

圖2 整體程序框圖

3 結(jié)論

針對目前電磁預(yù)測法中存在無法對破裂源實(shí)現(xiàn)準(zhǔn)確定位的問題，本文結(jié)合煤與瓦斯突出時的實(shí)際模型，推演出煤巖破裂時電偶極子場源電磁波的正向傳播規(guī)律，將理論推演進(jìn)行計算機(jī)編程來開發(fā)相應(yīng)系統(tǒng)軟件。系統(tǒng)包括場強(qiáng)計算子系統(tǒng)數(shù)據(jù)庫管理子系統(tǒng)等，為根據(jù)接收信號隊破裂源進(jìn)行定位提供了理論依據(jù)。

［1］王恩元,何學(xué)秋,聶百勝,等.電磁輻射法預(yù)測煤與瓦斯突出原理研究[D].中國礦業(yè)大學(xué),2000.

［2］張清毅.透地通信信道的研究—一種關(guān)于分層媒質(zhì)場強(qiáng)計算的新方法及其應(yīng)用[D].北京:中國礦業(yè)大學(xué)〈北京〉,1995.

［3］鄭鈞.D.電磁場與波[M].上海交通大學(xué)出版社,1984.

［4］王凱,俞啟香.煤與瓦斯突出的非線性特征及預(yù)測模型[M].中國礦業(yè)大學(xué)出版社,2005.

［5］S.H.Ward.Electromagnetic theory for geophysical applications，Mining Geophysics[Z].1978.

［6］楊永國,主編.VISUAL C++6.0實(shí)用教程[M].清華大學(xué)出版社,2007.