姚藝, 王一博，王偉君，陳凌

1 中國科學(xué)院地質(zhì)與地球物理研究所, 北京 100029 2 中國科學(xué)院大學(xué)，北京 100049 3 中國地震局地震預(yù)測研究所，北京 100036

0 引言

分布式光纖聲波傳感技術(shù)(Distributed Acoustic Sensing, DAS)是一種新型的地震采集技術(shù).它利用光纖既是傳輸介質(zhì)也是傳感設(shè)備這一優(yōu)勢，獲取光纖上每一點的應(yīng)變/應(yīng)變率.其主要工作原理是，外部擾動(如地震波等)引起光纖中微小不均勻體的位置發(fā)生變化，從而導(dǎo)致光纖中的散射光信號發(fā)生變化，通過解調(diào)散射光信號便可以獲得沿光纖上每個點的應(yīng)變信息(Parker et al., 2014).

相較于傳統(tǒng)的檢波器，DAS具有以下優(yōu)點：首先是高密度數(shù)據(jù)采樣，DAS可以達到米級的空間采樣率，對觀測區(qū)域的覆蓋密度更高，只要布設(shè)一條光纜，就可以獲得數(shù)千上萬個地震觀測點；其次，DAS的傳感介質(zhì)是光纖，光纖本身成本低廉，適用于長距離布設(shè)；最后，光纖還具有耐高溫和耐腐蝕性等特性，可以永久性布設(shè)，長時間采集地震數(shù)據(jù).

DAS的應(yīng)用領(lǐng)域非常多，已經(jīng)從早期的定性測量，比如周界入侵、油氣管道安全監(jiān)測等，逐漸發(fā)展為定量測量.井中監(jiān)測是DAS最早的地震學(xué)應(yīng)用，主要包含垂直地震剖面成像(Bakku et al., 2014)、微地震監(jiān)測(Lellouch et al., 2020)等.DAS在近地表結(jié)構(gòu)以及城市地下空間探測中也發(fā)揮了重要作用，如通過DAS噪聲分析可以獲得500 m以內(nèi)的淺層結(jié)構(gòu)信息(Shragge et al，2021).此外，DAS還應(yīng)用于火山(Nishimura et al., 2021)、冰川(Walter et al., 2020)等地震監(jiān)測中.

DAS與傳統(tǒng)的三分量檢波器不同，其對沿著光纖軸向方向敏感.因此DAS能夠提供單分量的地震數(shù)據(jù)，且這些數(shù)據(jù)與地震波的傳播方向和振動方向相關(guān).地震勘探已經(jīng)逐漸開始使用三分量地震數(shù)據(jù)進行成像等工作，因此亟需三分量的分布式光纖聲波傳感技術(shù).有學(xué)者提出通過設(shè)計螺旋纏繞的光纜結(jié)構(gòu)，實現(xiàn)地震波場的三分量應(yīng)變測量(Ning and Sava, 2018; Innanen et al., 2019; 曹丹平等，2022).

盡管DAS已經(jīng)在地震勘探等領(lǐng)域得到很多應(yīng)用，但還有一些理論知識有待研究.本文致力于研究DAS的應(yīng)變輻射花樣及微振動響應(yīng)特征，首先通過理論推導(dǎo)獲取應(yīng)變輻射花樣，分析其與速度輻射花樣的區(qū)別；然后通過數(shù)值模擬，展示不同震源機制情況下的應(yīng)變輻射花樣以及不同測線的DAS微振動響應(yīng).

1 震源模型

震源模型主要包含單力源模型及力偶模型，其中單力源模型中的震源是某個方向的作用力，比如錘擊實驗中的錘擊震源.力偶模型中的震源通常是一對方向相反的作用力，一般用矩張量來表示(Gilbert, 1970).震源矩張量可分解為各向同性震源(Isotropic，ISO)、雙力偶震源(Double Couple，DC)以及補償線性矢量偶極子震源(Compensated Linear Vector Dipole, CLVD)(Knopoff and Randall，1970).各向同性震源用歸一化單位矩陣表示，代表爆炸震源或者壓縮震源，在各向同性均勻介質(zhì)中只產(chǎn)生P波.雙力偶震源也即剪切震源，是較為常見的震源模型.CLVD震源一般與ISO震源共存，不涉及體積變化，其某個軸的膨脹(壓縮)會被另兩個軸的壓縮(膨脹)進行補償.

圖1 一般位錯模型示意圖

(1)

Mij=cijklDkl,

(2)

(3)

由式(2)可知，震源區(qū)的矩張量不僅與所處區(qū)域的彈性參數(shù)有關(guān)，也與斷層滑移的方向有關(guān).通過設(shè)計不同的震源參數(shù)，可以得到不同的震源模型.

2 微振動輻射花樣

單力源的位移記錄表達式為(Aki and Richard，2002)：

(4)

(5)

(6)

由于應(yīng)變是位移的空間導(dǎo)數(shù)，為此通過位移對空間方向進行求導(dǎo)，可以得到應(yīng)變的數(shù)學(xué)表達式：

(7)

f(xi)=xixp,

(8)

g(xi)=r-3,

(9)

(10)

根據(jù)求導(dǎo)法則，可推導(dǎo)單力源的P波應(yīng)變?yōu)椋?/p>

(11)

在式(11)中，前半部分與1/r2成正比，為中場項，后半部分與1/r成正比，為應(yīng)變的遠場項.同樣對應(yīng)變的輻射花樣選取遠場項部分，將上述方法應(yīng)用于單力源的S波可以得到：

(12)

因此單力源的P波和S波的應(yīng)變輻射公式為：

(13)

(14)

類似地，根據(jù)Aki和Richard(2002)的位移計算可知，矩張量產(chǎn)生的位移場包含近場項、中場項以及遠場項.Song和Toks?z(2011)認為遠場項在大約1.5個波長距離后成為主導(dǎo)項.因此對于輻射圖案的分析，也都是選擇遠場項.在均質(zhì)、各向同性介質(zhì)中，由矩張量震源引起的P波位移的表達式為：

(15)

其中遠場項可表示為：

(16)

(17)

位移輻射花樣表達式為：

(18)

(19)

將位移對空間方向求導(dǎo)，可得P波應(yīng)變的表達式：

(20)

同理，S波的應(yīng)變表達式為：

(21)

同樣選取遠場項，矩張量震源的應(yīng)變輻射花樣公式為：

(22)

(23)

表1列出了本文所用的震源模型，其中后5種震源均包含DC分量.圖2—圖9展示了表1對應(yīng)震源模型的P波和S波在x1x3平面上的速度輻射花樣和應(yīng)變輻射花樣.

圖2 單力源在x1x3平面內(nèi)x1和x3分量的速度輻射花樣和應(yīng)變輻射花樣

圖3 ISO震源在x1x3平面內(nèi)x1和x3分量的速度輻射花樣和應(yīng)變輻射花樣

圖4 CLVD震源在x1x3平面內(nèi)x1和x3分量的速度輻射花樣和應(yīng)變輻射花樣

圖5 DC震源在x1x3平面內(nèi)x1和x3分量的速度輻射花樣和應(yīng)變輻射花樣

圖6 逆沖斷層震源在x1x3平面內(nèi)x1和x3分量的速度輻射花樣和應(yīng)變輻射花樣

圖7 正斷層震源在x1x3平面內(nèi)x1和x3分量的速度輻射花樣和應(yīng)變輻射花樣

圖8 一般位錯震源在x1x3平面內(nèi)x1和x3分量的速度輻射花樣和應(yīng)變輻射花樣

圖9 全矩張量震源在x1x3平面內(nèi)x1和x3分量的速度輻射花樣和應(yīng)變輻射花樣

表1 本文采用的震源模型

從輻射花樣圖中可以看出，無論是何種震源，P波的速度輻射花樣與應(yīng)變輻射花樣形狀上非常相似，但S波的速度輻射花樣與應(yīng)變輻射花樣的差別較大，尤其是以DC成分為主的震源類型，可以看到這種震源的S波速度輻射花樣分為6瓣，而在應(yīng)變輻射花樣中分為8瓣.這也證明了這些震源的DAS應(yīng)變響應(yīng)中，S波極性將會發(fā)生更多的反轉(zhuǎn).另外，無論是P波還是S波，其應(yīng)變輻射花樣的極性與速度輻射花樣的極性是對稱相反的.考慮x1分量極性變化，如果速度的輻射花樣兩端極性相反，那么應(yīng)變輻射花樣的兩端極性是一致的，反過來說如果速度的輻射花樣兩端極性相同，那么應(yīng)變輻射花樣的兩端極性是相反的.這些現(xiàn)象在后續(xù)的理論合成記錄中也能夠得到佐證.同樣的，x3分量的極性變化也是如此，在這里將不再贅述.

3 分布式光纖聲波傳感微振動數(shù)值模擬方法

分布式光纖能記錄到軸向方向上的應(yīng)變，根據(jù)應(yīng)變與位移的關(guān)系式：

(24)

將位移數(shù)據(jù)投影到光纖上，計算沿光纖軸向的空間導(dǎo)數(shù)，可將位移數(shù)據(jù)投影到光纖上得到合成應(yīng)變數(shù)據(jù)，但這并不是實際DAS測量獲得的應(yīng)變.DAS獲得的是沿光纖軸向有限長度內(nèi)的平均應(yīng)變，這個有限長度也被稱為標準距離(標距，gauge length).DAS應(yīng)變可以用公式(25)來描述：

(25)

其中L是標距.標距的大小對DAS測量性能影響較大，標距越小，DAS空間分辨率越高，但信噪比較低；標距越大，DAS空間分辨率越低，但信噪比較高(Dean et al.,2017; Hartog, 2017).一般來說，標距需要在測量距離、空間分辨率和信噪比中尋找平衡.

公式(25)提供了計算DAS應(yīng)變的一種方法.此外，還可以利用光纖的瑞利散射干涉模型，將光纖上每點的應(yīng)變轉(zhuǎn)化為瑞利散射光的相位，模擬真實的瑞利散射光信號，然后通過解調(diào)后向瑞利散射光信號，獲得真實DAS應(yīng)變響應(yīng)(方欣棟，2018；馬國旗等，2020).

4 數(shù)值算例

為了研究不同震源機制以及不同布設(shè)條件的分布式光纖微振動響應(yīng)特征，本文建立了如圖10所示的數(shù)值模型，該模型對應(yīng)的觀測系統(tǒng)有7條地面光纖測線，前6條光纖的長度均為2000 m，第7條測線(方形測線)的長度為4000 m，此外還有一條井中光纖測線8，其長度為1200 m.DAS系統(tǒng)參數(shù)設(shè)計如下：空間采樣間隔為0.1 m，標距為1 m.震源位于地表以下200 m處，介質(zhì)的縱波速度為3000 m·s-1, 橫波速度為1731 m·s-1.采用表1設(shè)計的8種震源，使用50 Hz的雷克子波作為震源函數(shù)，根據(jù)上述理論模擬方法，可以獲得各種情況下的合成數(shù)據(jù).

圖10 (a)DAS觀測系統(tǒng)，地面有7條測線，井中有一條測線.(b)觀測系統(tǒng)的俯視圖，地面測線按照逆時針順序標記為Line 1—Line 6，方形測線標記為Line 7，井中測線標記為Line 8

圖11是單力源的速度記錄和DAS應(yīng)變響應(yīng).圖12是ISO震源的速度記錄和DAS應(yīng)變響應(yīng).由于ISO震源是一種對稱源，且震源處在測線的中心，因此6條地面測線的觀測記錄是一樣的.圖13是CLVD震源的速度記錄與DAS應(yīng)變響應(yīng).本文采用的CLVD震源相對于水平面來說是對稱的，因此6條地面測線獲得的記錄也是一樣的.對比圖11至圖13的速度記錄和DAS應(yīng)變響應(yīng)，可以看到速度記錄有極性反轉(zhuǎn)現(xiàn)象，而DAS應(yīng)變響應(yīng)沒有極性反轉(zhuǎn)現(xiàn)象.

圖11 單力源的速度記錄(左列)和DAS應(yīng)變響應(yīng)(右列)；從上至下分別代表測線1—6的記錄

圖12 ISO震源的速度記錄(左列)和DAS應(yīng)變響應(yīng)(右列)；從上至下分別代表測線1—6的記錄

圖14是DC震源的速度記錄與DAS應(yīng)變響應(yīng)，DC震源的DAS應(yīng)變響應(yīng)特征與測線布設(shè)方式有一定的相關(guān)性，例如測線4記錄不到信號，而測線1、2、3記錄的信號是相似的.與圖11至圖13觀察到的現(xiàn)象不同，圖14的DAS應(yīng)變響應(yīng)存在較多的極性反轉(zhuǎn)現(xiàn)象.圖15、圖16和圖17分別是逆沖斷層震源、正斷層震源與一般位錯震源的速度記錄與DAS應(yīng)變響應(yīng)，相較于DC震源的模擬結(jié)果，測線4可以記錄到信號，但其他測線DAS應(yīng)變響應(yīng)的極性反轉(zhuǎn)位置發(fā)生了變化.

圖14 DC震源的速度記錄(左列)和DAS應(yīng)變響應(yīng)(右列)；從上至下分別代表測線1—6的記錄

圖15 逆沖斷層震源的速度記錄(左列)和DAS應(yīng)變響應(yīng)(右列)；從上至下分別代表測線1—6的記錄

圖16 正斷層震源的速度記錄(左列)和DAS應(yīng)變響應(yīng)(右列)；從上至下分別代表測線1—6的記錄

圖17 一般位錯震源的速度記錄(左列)和DAS應(yīng)變響應(yīng)(右列)；從上至下分別代表測線1—6的記錄

圖18展示的是全矩張量震源的速度記錄與DAS應(yīng)變響應(yīng)，綜合圖14至圖18可以發(fā)現(xiàn)，含DC成分震源模型DAS應(yīng)變響應(yīng)具有更多的極性反轉(zhuǎn)現(xiàn)象，例如圖18中S波在測線1至測線4上產(chǎn)生的DAS應(yīng)變響應(yīng)具有三次極性反轉(zhuǎn)現(xiàn)象，而在測線5和測線6上具有兩次極性反轉(zhuǎn)現(xiàn)象.

圖18 全矩張量震源的速度記錄(左列)和DAS應(yīng)變響應(yīng)(右列)；從上至下分別代表測線1—6的記錄

圖19是測線7接收到的8種震源的速度記錄與DAS應(yīng)變響應(yīng).由于該測線是方形的，屬于二維測線，圖19中的橫坐標不再是偏移距而是光纖長度，可以看到每張圖中都會有4個同相軸.圖20是測線8接收到的8種震源的速度記錄與DAS應(yīng)變響應(yīng).為了模擬水平井光纖數(shù)據(jù)采集情況，將光纖的布設(shè)方式設(shè)計為L型，記錄到的信號分為垂直段和水平段.對比圖11至圖20，可以發(fā)現(xiàn)相同震源的速度記錄與應(yīng)變響應(yīng)是相似的.

圖19 測線7對應(yīng)的8種震源的速度記錄和DAS應(yīng)變響應(yīng)；按照由上而下，從左到右的順序，8種震源分別是：單力源(a, b)、ISO震源(c, d)、CLVD震源(e, f)、DC震源(g, h)、逆斷層震源(i, j)、正斷層震源(k, l)、一般位錯震源(m, n)、全矩張量震源(o, p)

圖20 測線8對應(yīng)的8種震源模型的速度記錄和DAS應(yīng)變響應(yīng)；按照由上而下，從左到右的順序，8種震源分別是：單力源(a, b)、ISO震源(c, d)、CLVD震源(e, f)、DC震源(g, h)、逆斷層震源(i, j)、正斷層震源(k, l)、一般位錯震源(m, n)、全矩張量震源(o, p)

5 結(jié)論

本文研究了不同震源機制、不同布設(shè)方式的DAS應(yīng)變響應(yīng)特征，主要結(jié)論有以下幾點：

(1)應(yīng)變輻射花樣與速度輻射花樣有區(qū)別.從形態(tài)上來說，各種震源機制下P波的速度輻射花樣和應(yīng)變輻射花樣形狀較為相近，S波在含DC成分震源機制條件下兩者的形態(tài)上差別明顯.以震源所在位置的垂向方向或水平向方向作為對稱軸，如果速度輻射花樣兩端的極性是相同或相反，那么相對應(yīng)的應(yīng)變輻射花樣兩端極性就是相反或相同.

(2)理論合成記錄結(jié)果表明含DC成分震源的S波DAS響應(yīng)相較于檢波器記錄，存在更多的極性反轉(zhuǎn)現(xiàn)象，且多出的極性反轉(zhuǎn)一般發(fā)生在零偏移距處.

(3)ISO、CLVD震源由于震源自身的對稱性，其應(yīng)變響應(yīng)與光纖的布設(shè)方式?jīng)]有關(guān)系；而含DC成分震源在不同布設(shè)測線上的DAS響應(yīng)變化不同，比如純DC震源在某些測線方向可能接收不到任何響應(yīng)，體現(xiàn)了DAS對方向的敏感性.

本文研究的是各向同性均勻介質(zhì)情況下的DAS響應(yīng)特征，但真實地下介質(zhì)情況非常復(fù)雜，需要進一步研究各向異性非均勻介質(zhì)中的DAS響應(yīng)特征.