吳洋，嚴家斌

(1.中南大學 地球科學與信息物理學院，湖南 長沙 410083； 2.有色資源與地質(zhì)災害探查湖南省重點實驗室，湖南 長沙 410083)

0 引言

大地極化聲子測深技術(shù)是由烏克蘭科學家尤·阿·博格達諾夫提出的一種被動源地球物理電磁勘探方法，地球內(nèi)部巖石的應力變化導致巖石產(chǎn)生變形以及破壞斷裂，從而產(chǎn)生聲發(fā)射和電磁輻射，該方法通過測量電磁波輻射強度來研究地球內(nèi)部結(jié)構(gòu)[1]。大地形變產(chǎn)生的電磁輻射具有較寬的頻率范圍，理論上，該方法可用于研究大地的地質(zhì)結(jié)構(gòu)。目前，在國內(nèi)已經(jīng)開始進行了一些試驗和研究[2]。

大地構(gòu)造運動(如地震)會導致地球內(nèi)部應力場發(fā)生巨大變化，造成巖石中斷裂和裂隙產(chǎn)生。同樣，由于頁巖儲層的低孔和極低基質(zhì)滲透率的特點，頁巖氣的開發(fā)主要采用壓裂技術(shù)，頁巖氣通常含有天然裂縫，一般會因為膠結(jié)而封閉，這些天然裂縫容易在壓裂中破裂[3]，這個過程也會造成巖石的破裂。巖石破裂是指巖石在一定的應力條件下，巖石中裂紋開始萌生或者擴大。1966年Cook通過實驗研究了巖石的應力—應變曲線，巖石的失穩(wěn)破裂并非發(fā)生在峰值強度點，而是在峰值載荷后。Cook的研究開啟了巖石破裂過程的新篇章[4]。華保欽歸納了巖石破裂的3個應力場條件：莫爾圓直徑增大；異常壓力發(fā)育時，造成一定深度下孔隙流體壓力與上覆負荷之比增大，有效應力減小；最小主應力為張力，最大主應力與最小主應力的大小差距小[5]。唐志強等使用分離式SHPB裝置進行巖石動力學實驗，發(fā)現(xiàn)火山碎屑巖以劈裂破壞為主，砂巖以壓碎破壞為主[6]。高美奔等采用統(tǒng)計方法對不同實驗條件下花崗巖的宏觀破裂特征進行分析，發(fā)現(xiàn)低溫低壓、低溫高壓和高溫高壓時巖樣以剪切破壞為主，高溫低壓時為張拉破壞。巖樣性能主要受圍壓影響，溫度對巖樣破裂方式有一定影響[7]。Glynn O.Cress等通過單軸壓縮實驗，記錄到巖石破裂時產(chǎn)生的電磁輻射[8]。龔強等利用張開位移法研究了巖石產(chǎn)生的電磁輻射頻率與巖石彈性參數(shù)之間的關(guān)系，表明電磁輻射頻率除了與試件的尺寸有關(guān)外，還與巖石的彈性參數(shù)有關(guān)[9]。

對于巖石中電磁輻射產(chǎn)生的機制，許多學者提出了基于實驗的觀點，認為巖石產(chǎn)生的電磁輻射與聲發(fā)射有關(guān)[10-12]。聲發(fā)射是指巖石變形和破裂過程中以彈性波形式釋放的應變能現(xiàn)象，聲發(fā)射是巖石變形和破裂的直接體現(xiàn)。20世紀30年代，Obert L等發(fā)現(xiàn)巖石破裂的聲發(fā)射現(xiàn)象[13]。 Kaiser觀察到金屬鋅、銅、鋁及鉛都有聲發(fā)射現(xiàn)象，并發(fā)現(xiàn)了聲發(fā)射的不可逆效應——Kaiser效應。唐春安通過模擬實驗得到了巖石聲發(fā)射的3種基本模式：主震型、前震-主震-余震型、群震型[14]。張國凱等采用巖石加載系統(tǒng)獲得脆性花崗巖單軸壓縮全應力—應變曲線，同步測試聲發(fā)射和超聲波變化，研究了聲發(fā)射隨巖石破裂的變化規(guī)律[15]。趙伏軍等進行了刀具靜力侵入花崗巖破碎實驗，發(fā)現(xiàn)刀具靜力侵入巖石破碎的過程中，聲發(fā)射和電磁輻射有著很好的正相關(guān)性[16]。曹惠馨等通過實驗發(fā)現(xiàn)電磁輻射與聲發(fā)射的最大值在時間上對應的約占70%，認為巖石中的電磁輻射信號主要是由巖石破裂引起，除此之外可能還與壓電效應等有關(guān)[17]。Yamada等進行了巖石破裂電磁輻射與聲發(fā)射的實驗，發(fā)現(xiàn)電磁輻射的起始時間與聲發(fā)射的起始時間一致，并且拉伸裂紋比剪切裂紋產(chǎn)生電磁輻射的效率更高，認為裂紋是產(chǎn)生電磁輻射的必要條件[18]。郭自強等通過實驗發(fā)現(xiàn)聲發(fā)射與電磁輻射的關(guān)系具有多樣性，既有伴隨聲發(fā)射的電磁輻射，也有不伴隨聲發(fā)射的電磁輻射，認為巖石破裂產(chǎn)生的電磁輻射應來源于多種機制[19]。徐為民等發(fā)現(xiàn)巖石破裂時，電磁脈沖與聲發(fā)射事件并不能完全一一對應，認為電磁輻射直接由巖石破碎產(chǎn)生[20]。王恩元等進行了煤體受載變形破裂的實驗研究，結(jié)果表明煤體變形破裂時能夠產(chǎn)生電磁輻射。但電磁輻射信號較聲發(fā)射信號豐富，認為電磁輻射與煤體的變形破裂過程相關(guān)[21]。

研究者認為震電效應可能是引起巖石中產(chǎn)生電磁輻射的原因之一[22]。震電效應是指地震波和電磁波能量的耦合與轉(zhuǎn)換，即機械力與電磁力的耦合，壓電效應是其中的一種。Gokhberg認為巖石的壓電效應等力電現(xiàn)象是巖石在受載情況下，由分離電荷張弛產(chǎn)生的電磁輻射效應[23]。劉煜洲等認為巖石破裂的低頻電磁輻射是由壓電效應產(chǎn)生的[24]。Huang從壓電效應與斷層位錯理論模型出發(fā)，解釋了同震電磁輻射信號的產(chǎn)生機制[25]。劉志祥研究了沖擊載荷作用下花崗巖產(chǎn)生的電磁輻射，分析了花崗巖產(chǎn)生電磁輻射的壓電效應和斷裂動力學機理[26]。李夕兵等認為壓電效應是巖石產(chǎn)生電磁輻射的原因，并從壓電效應出發(fā)，闡明了應力波和電磁波的耦合機制，解釋了地震和巖石破碎產(chǎn)生的聲、光、電現(xiàn)象[27]。雖然目前學界對于巖石破碎產(chǎn)生電磁輻射的機制還不能完全了解，但都不否認壓電效應對于巖石產(chǎn)生電磁輻射的貢獻。

盡管博格達諾夫提出了大地極化聲子的概念，國內(nèi)的研究學者們也進行了實驗研究，但其極化機理并不清楚，一般認為這種大地極化聲子可能是由大地構(gòu)造活動(如地震)或人類活動(如頁巖氣壓裂)產(chǎn)生的壓電效應所致。本文嘗試從大地巖石破碎時的電磁發(fā)射過程出發(fā)，以壓電效應為基礎，開展對大地極化聲子測深技術(shù)的模擬研究，通過巖石破裂的壓電應變矩陣和壓電本構(gòu)方程，導出巖石中電磁輻射與外力作用下的變換關(guān)系，模擬多種裂隙形態(tài)下的電磁輻射效應，分析電磁輻射特征。

1 壓電電磁輻射原理

大地極化聲子測深技術(shù)的原理是機械運動使巖石的應力狀態(tài)發(fā)生變化，巖石因受到擠壓產(chǎn)生變形直至破碎，產(chǎn)生極化聲子，從而發(fā)出電磁信號。這種極化聲子可能是由大地構(gòu)造活動或人類活動產(chǎn)生的壓電效應所致。

選擇坐標軸x、y、z與晶軸X、Y、Z重合，壓電介質(zhì)的壓電應變矩陣[28]為：

(1)

式中：d表示壓電應變常數(shù)，第一個下標表示所產(chǎn)生電位移的方向(1、2、3分別表示x，y，z方向)，第二個下標表示作用力方向(1、2、3、4、5、6分別表示沿x、y、z三個軸方向的應力和垂直于x、y、z三個軸平面(XZ、YZ、XY平面)作用的剪切力)。

在地殼上部，石英占66.4%，并且是許多巖石的成巖礦物，這為壓電效應提供了物質(zhì)基礎，故用石英的壓電應變矩陣來代替地殼中含有石英等壓電介質(zhì)的巖石的壓電應變矩陣。石英晶體的壓電應變常數(shù)與其內(nèi)部結(jié)構(gòu)有關(guān)，屬于三角晶系32點群。有3個晶軸：x軸、y軸和z軸，當晶體沿x軸或者y軸方向受到壓力時，晶格變形，產(chǎn)生壓電效應；沿z軸方向受到力時，由于晶格的變形不會引起正負電荷中心的分離，不會產(chǎn)生壓電效應。因此，含有石英等壓電介質(zhì)的巖體的壓電應變矩陣[29]為：

(2)

設巖體中有一個沿x方向傳播的縱波：

T1=T0e-η0xei(ω0t-k0x)，

(3)

式中：ω為應力波的角頻率，T0為應力波的初始振幅值，k0是應力波的波數(shù)，η0為應力波的衰減系數(shù)。巖石受應力波作用而變形，產(chǎn)生壓電效應。根據(jù)壓電本構(gòu)方程[30]：

(4)

式中：D為電位移矢量，T為應力波，E為外加電場強度，dT為d的轉(zhuǎn)置矩陣，εT表示應力恒定時的介電常數(shù)，sE表示場強恒定時的彈性柔順常數(shù)。

由式(2)、(3)、(4)可以得到應力場對電位移矢量D壓電產(chǎn)生壓電貢獻(不考慮外加電場，即E=0)，由于應力波沿x方向?qū)α严妒┘幼饔昧?，在垂直x軸的平面上產(chǎn)生電荷，因此只需考慮式(2)中的d11參數(shù)：

D壓電=d11T1，

(5)

D壓電是諧變的，服從麥克斯韋方程組：

(6)

以及電磁媒質(zhì)的本構(gòu)關(guān)系式：

(7)

式中：B、H、D和E分別為電磁場的磁感應強度、磁場強度、電位移矢量和電場強度；μ為磁導率；Jc為傳導電流密度，在巖體中設Jc=0。

由式(6)、(7)得到：

(8)

式中：k是電磁波的復波數(shù)，σ為電導率，ω為電磁波的角頻率。

將式(5)代入式(8)，可得應力波作用下巖石產(chǎn)生電磁輻射的方程：

(9)

(10)

由于電磁波在巖石介質(zhì)中具有色散特性，將電磁波用復數(shù)形式表示，設：

E=Eeei(ωt)，

(11)

將式(11)代入式(10)可得：

(12)

利用邊界條件及電磁波的極限條件，解方程得：

(13)

應力波作用下，巖石產(chǎn)生的電磁輻射方程為：

(14)

電磁輻射信號在巖石介質(zhì)中傳播時會發(fā)生衰減及相變，引入e-ikr。由式(6)得到電磁波的復波數(shù)、衰減系數(shù)和相位系數(shù)[31]：

(15)

(16)

(17)

式中：k為復波數(shù)，α是電磁波的衰減系數(shù)，β為電磁波的相位系數(shù)，σ為大地介質(zhì)的電導率，ω為電磁波的角頻率，r表示電磁波傳播的距離。應力波作用下，巖石產(chǎn)生電磁輻射的傳播方程為：

(18)

2 發(fā)射頻率

Rabinovitch等通過大量實驗發(fā)現(xiàn)巖石產(chǎn)生的電磁輻射頻率受裂紋寬度的限制，頻率與裂紋寬度成反比，并給出了電磁輻射頻率與裂紋寬度之間的關(guān)系[32]：

(19)

式中：ω為角頻率，vR為瑞利波速，b為裂紋寬度。

根據(jù)波動理論以及瑞利波速公式[33]：

(20)

(21)

式中：f為頻率，μb為泊松比，ET為彈性模量，ρ為密度?？梢缘玫綆r石產(chǎn)生的電磁輻射頻率與裂紋寬度的關(guān)系：

(22)

顯然電磁輻射頻率與裂紋寬度及巖石的物理性質(zhì)有關(guān)。含石英類巖石的泊松比范圍為0.10～0.40，彈性模量一般為10～100 GPa，密度為2.6～2.9 kg/m3，故取泊松比μb=0.25，密度ρ=2.65×103kg/m3，彈性模量ET=50 GPa。圖1a為根據(jù)式(22)得到的電磁輻射頻率隨裂紋寬度的變化曲線，可以看到巖石破裂產(chǎn)生的電磁輻射頻率較高，一般大于105Hz，同時電磁輻射頻率隨裂紋寬度的增大呈指數(shù)減小。

取裂隙寬度b=0.005 m，密度與彈性模量與上述相同，得到電磁輻射頻率隨泊松比的變化如圖1b所示，可以看到電磁輻射頻率隨泊松比的增大呈線性衰減。類似地可以得到電磁輻射頻率隨彈性模量與密度的變化如圖1c、d所示，可以看到電磁輻射頻率隨彈性模量的增大而增大，隨密度的增大而減小。比較圖1可以看到電磁輻射頻率受泊松比和彈性模量影響較小，受裂隙寬度和密度影響較大。

圖1 電磁輻射頻率隨巖石參數(shù)的變化Fig.1 Variation curve of electromagnetic radiation frequency with rock parameters

3 數(shù)值模擬與分析

通過對裂隙單模型和組合模型的電磁輻射響應特征的數(shù)值模擬，討論基于壓電效應的大地極化聲子測深技術(shù)的有效性。

3.1 單模型

設地下電阻率為1 000 Ω·m的均勻半空間中存在1個三角形裂隙(圖2)，傾角θ分別為10°、30°、60°和90°；裂隙埋深H=500 m，裂隙長度L0=50 m，最大寬度B0=0.01 m。應力波作用時間t=1 s，應力波角頻率ω0=(2π×10) rad/s，應力波波數(shù)k0=0.01 m-1，衰減系數(shù)η0=0.05。磁導率為真空磁導率μ=4π×10-7H/m。根據(jù)含石英類巖石的介電常數(shù)為(4～9)×8.85×10-12F/m和壓電應變常數(shù)范圍為(0.5～5)×10-13C/N，故取介電常數(shù)為ε=5×8.85×10-12F/m，壓電應變常數(shù)d11=10-13C/N。泊松比μb=0.25，密度ρ=2.65×103kg/m3，彈性模量ET=50 GPa，電導率σ=0.001 S/m。設有一個水平方向的應力波，其作用力T0=107Pa，且電磁輻射由應力垂直于裂隙方向的分力產(chǎn)生。

圖2 傾斜裂隙模型示意Fig.2 Schematic diagram of inclined fracture model

圖3是巖石裂隙電磁發(fā)射的地表電場幅值。由圖3a可以看到：電場幅值曲線的變化趨勢為中間高兩邊低，電磁發(fā)射頻率越低，電場幅值就越大；隨著頻率的增加，幅值快速下降；電場幅值的極值點位置和大小隨頻率變化，頻率越低，極值點越大，越偏離中心(裂隙頂點在地表的投影)。

將不同產(chǎn)狀的裂隙電場幅值圖進行對比，可以發(fā)現(xiàn)傾斜程度越大(即傾角越小)的裂隙，電場(從低頻到高頻)的幅值和極值點越小，極值點偏離坐標軸中心也越遠。電場幅值的位置以及變化趨勢能夠反映出裂隙的位置和產(chǎn)狀。極值點的連線方向與裂隙的傾斜方向相反；連線方向與地表鉛垂方向的夾角反映了裂隙的傾斜度的大小。夾角越大，裂隙的傾斜度越大。通過曲線的頻率能夠定性分析裂隙的寬度變化。

圖4分別是不同傾角的裂隙的電場矢量圖(z軸向上為正方向，x軸向右為正方向)。以圖4d為例，電場的矢量方向沿一個方向變化(順時針方向)，矢量幅值變化趨勢為增大—降低—增大—降低，在正負方向上各擁有一個極值;在z軸正方向上極值與負方向上的極值相近。

將不同產(chǎn)狀的裂隙電場矢量圖進行對比，可以發(fā)現(xiàn)豎直裂隙的z軸正方向極值與負方向極值相同。傾斜程度越大(即傾角越小)的裂隙，其z軸負方向上的極值相比正方向也就越小，直到30°的裂隙的負值消失。根據(jù)電場矢量圖所顯示的地表電場的幅值、方向及變化趨勢，能夠大致判斷出產(chǎn)生電磁輻射的電場源的產(chǎn)狀和位置以及深度。

圖3 裂隙電場幅值隨傾角變化Fig.3 The electric field amplitude of the fissure varies with the inclination angle

圖4 裂隙電場矢量隨傾角變化Fig.4 Crack electric field vector changes with inclination angle

3.2 組合模型

設地下電阻率為1 000 Ω·m的均勻半空間，存在2個分開的豎直的完全一致的三角形裂隙(圖5)，裂隙埋深H=500 m，長度L0=50 m,其他參數(shù)同上。2個裂隙都受到水平方向的應力波T0=107Pa的作用。

圖5 組合裂隙模型示意Fig.5 Schematic diagram of combined fracture model

從圖6可以看到：電場幅值由雙峰組成，頻率越低，曲線的幅值越大，且所有頻率電磁輻射幅值的極值點都位于x=-200 m或x=200 m線上。從圖7可以看到，電場的矢量并不是沿一個方向變化，出現(xiàn)了重復的矢量方向，矢量在z軸正方向上有2個極值，左邊極值大于右邊；在負方向上也有2個極值，但是右邊極值大于左邊。通過圖6和圖7，能夠大致判斷出組合裂隙的位置。

通過對比單裂隙模型可以發(fā)現(xiàn)組合裂隙與單裂隙的差別：在電場幅值圖中單異常體是單峰，而組合裂隙是多峰，并有波谷存在；在電場矢量圖中雖然單裂隙電場方向是沿順時針變化的，而組合裂隙的電場方向變化趨勢會出現(xiàn)轉(zhuǎn)變，單裂隙不會出現(xiàn)重復的矢量方向，且只有一個極大值和一個極小值，但組合裂隙有多個極大值和多個極小值。運用該方法能夠分辨組合裂隙與單裂隙，同時能夠大致判斷組合裂隙的位置。

圖6 組合裂隙電場幅值Fig.6 The electric field amplitude of the combined fracture

圖7 組合裂隙電場矢量Fig.7 Electric field vector of combined fracture

4 結(jié)論

本文從巖石的壓電效應出發(fā)，推導出在均勻半空間情況下巖石的電磁輻射的產(chǎn)生和傳播方程。巖石破裂產(chǎn)生的電磁輻射頻率與裂隙的尺度及波的彈性參數(shù)有關(guān)，裂紋寬度越大，產(chǎn)生的電磁輻射頻率越低，但一般高于105Hz；同時，電磁輻射頻率還會受到巖石的泊松比、彈性模量與密度的影響。

對不同產(chǎn)狀及形態(tài)的裂隙模擬表明，巖石破裂產(chǎn)生的電磁輻射響應特征明顯，在裂隙的上方有極值存在，極值大小隨裂隙的尺度及頻率變化，并且可以通過多個頻率的極值變化特征推斷裂隙的產(chǎn)狀，利用其電磁輻射響應及矢量圖可以識別裂隙的數(shù)量。模擬研究表明基于巖石壓電效應為基礎的大地極化聲子測深技術(shù)具有潛在的應用價值。