余洋洋 梁春濤 杜瑤 王宇璽 莊園旭

1)四川省地震局，成都 610041 2)成都理工大學(xué)，地球探測(cè)與信息技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室，成都 610059

0 引言

在頁(yè)巖氣開(kāi)發(fā)中，通常用水力壓裂的方式提高致密頁(yè)巖的滲透率，對(duì)水力壓裂效果的評(píng)價(jià)一般通過(guò)微地震監(jiān)測(cè)實(shí)現(xiàn)(撒利明等，2012；劉振武等，2013)。結(jié)合微地震事件的位置分布和震源機(jī)制信息，可以有效解釋壓裂施工區(qū)域裂縫的發(fā)育和應(yīng)力的分布，從而達(dá)到評(píng)價(jià)壓裂效果、指導(dǎo)壓裂設(shè)計(jì)和施工的目的(尹陳，2017；Herrmann et al，2011)。微地震監(jiān)測(cè)可以分為井中監(jiān)測(cè)和地面監(jiān)測(cè)，與井中監(jiān)測(cè)相比，地面監(jiān)測(cè)雖然有信號(hào)衰減嚴(yán)重、信噪比低等缺點(diǎn)，但在某些缺少鉆井或者鉆井困難的區(qū)域，地面監(jiān)測(cè)得到越來(lái)越廣泛的應(yīng)用(秦俐等，2013；余洋洋等，2017)。此外，地面監(jiān)測(cè)系統(tǒng)滿足了震源機(jī)制反演需要較好的方位角覆蓋的要求。

Kao等(2004)提出了震源掃描算法(Source Scanning Algorithm，簡(jiǎn)稱(chēng)SSA)，該方法被廣泛應(yīng)用于地面監(jiān)測(cè)定位微地震事件，但其無(wú)法反演震源機(jī)制，且對(duì)于信號(hào)較弱的剪切震源的識(shí)別能力較弱。Liang等(2016)在震源掃描算法的基礎(chǔ)上發(fā)展了微地震定位與震源機(jī)制聯(lián)合反演算法(Joint Source Scanning Algorithm，簡(jiǎn)稱(chēng)JSSA)，該方法可以同時(shí)反演微地震事件的位置和震源機(jī)制信息，且能夠較好地識(shí)別剪切震源。與其他方法相比，JSSA算法不需要清晰的波形以及準(zhǔn)確的到時(shí)信息，操作簡(jiǎn)便。作為較新的方法，JSSA算法已經(jīng)通過(guò)一系列理論合成數(shù)據(jù)和實(shí)際壓裂數(shù)據(jù)的測(cè)試，驗(yàn)證了其準(zhǔn)確性和穩(wěn)定性。

本研究采用微地震定位與震源機(jī)制聯(lián)合反演方法處理威H3-1井的實(shí)際壓裂數(shù)據(jù)。威H3-1井位于四川省內(nèi)江市資中縣銀山鎮(zhèn)，是長(zhǎng)寧-威遠(yuǎn)頁(yè)巖氣示范區(qū)的一口水平井，開(kāi)采施工時(shí)間自2014年10月29日至11月10日。該井與威203井屬同場(chǎng)井，在威203井開(kāi)采完成之后對(duì)其進(jìn)行鉆探施工，周邊還有威202井、威204井、威205井以及楠木寺斷層和板栗埡斷層(周太郎，2017)。根據(jù)中國(guó)地震臺(tái)網(wǎng)目錄(1)http://www.ceic.ac.cn/history.查詢統(tǒng)計(jì)2014年10月20日至11月20日威H3-1井周邊30km范圍內(nèi)的地震事件，如表1 所示。事件集中在壓裂施工區(qū)域西北側(cè)楠木寺斷裂附近，如圖1 所示。在壓裂施工期間及施工結(jié)束后10天內(nèi)，共有10個(gè)天然地震事件，最大震級(jí)為1.5級(jí)，最小為0.6級(jí)，平均震級(jí)1.0級(jí)，震源深度2～16km不等，平均深度7.1km，且事件均位于威H3-1井西北側(cè)楠木寺斷層附近。

表12014年10月20日—11月20日威H3-1周邊30km范圍內(nèi)地震統(tǒng)計(jì)

圖 1 壓裂施工區(qū)域周邊地質(zhì)概況及天然地震事件分布

威H3-1井工區(qū)地形較為平坦，海拔340～440m，水平井位于3000m深度附近，長(zhǎng)度1870m，共19個(gè)壓裂段自水平井西北端至東南端均勻分布，采用滑套多級(jí)壓裂技術(shù)。地面監(jiān)測(cè)系統(tǒng)由1193個(gè)檢波器組成，呈10個(gè)分支的星狀排列，每個(gè)分支有90～160個(gè)檢波器，相鄰檢波器間距約為30m。威H3-1井的工區(qū)狀況較為復(fù)雜，根據(jù)該井試油地質(zhì)設(shè)計(jì)方案等資料，工區(qū)周邊最大水平主應(yīng)力方向?yàn)镋W向，與水平井存在約45°夾角。

1 理論方法研究

傳統(tǒng)的震源掃描算法(SSA)需要對(duì)事件的空間坐標(biāo)x、y和z三個(gè)參數(shù)進(jìn)行掃描，聯(lián)合反演算法(JSSA)則需要對(duì)事件的空間坐標(biāo)x、y、z以及震源機(jī)制的走向、傾角、滑動(dòng)角共6個(gè)參數(shù)進(jìn)行掃描(Liang et al，2016)。

(1)

圖 2 地面監(jiān)測(cè)系統(tǒng)及虛擬震源空間示意圖

(2)

權(quán)重函數(shù)與虛擬震源i和檢波器j的空間幾何位置有關(guān)，其表達(dá)式為

(3)

JSSA算法實(shí)現(xiàn)過(guò)程的實(shí)質(zhì)是將波形資料從觀測(cè)時(shí)空域轉(zhuǎn)換到震源機(jī)制時(shí)空域，對(duì)于疊加能量bmi，在震源機(jī)制時(shí)空域中根據(jù)下式計(jì)算其每個(gè)采樣點(diǎn)的“亮度”函數(shù)

Rk=STAk/LTAk

(4)

其中，Rk為第k個(gè)采樣點(diǎn)的“亮度”函數(shù)；STAk為第k個(gè)采樣點(diǎn)之后的短時(shí)窗內(nèi)的均方根振幅，反映地震信號(hào)；LTAk為第k個(gè)采樣點(diǎn)之前的長(zhǎng)時(shí)窗內(nèi)的均方根振幅，反映背景噪聲。

在一定時(shí)空內(nèi)，若該采樣點(diǎn)的“亮度”函數(shù)R值大于預(yù)設(shè)的某一數(shù)值Rmin，那么該采樣點(diǎn)被認(rèn)為是一個(gè)事件，定義為該數(shù)值Rmin為拾取事件的閾值(李文軍等，2006；Grigoli et al，2013)?！傲炼取焙瘮?shù)R大于閾值的采樣點(diǎn)所對(duì)應(yīng)的時(shí)刻為發(fā)震時(shí)刻，對(duì)應(yīng)的虛擬震源點(diǎn)的坐標(biāo)為事件的空間位置，對(duì)應(yīng)的震源機(jī)制信息為事件的震源機(jī)制解。

由于地面監(jiān)測(cè)到的大部分微地震信號(hào)的主頻在20～30Hz，因此，在本研究中STA的時(shí)窗長(zhǎng)度固定為0.05s，接近主頻信號(hào)的一個(gè)周期(Eisner et al，2013)。基于大量測(cè)試，LTA的時(shí)窗長(zhǎng)度為1.0s，拾取事件的閾值Rmin=3.0(吳治濤等，2010；余洋洋等，2017)。由于實(shí)際壓裂數(shù)據(jù)的振幅為處理后的結(jié)果，無(wú)法確定震級(jí)，因而本研究得到的R值為事件的相對(duì)大小。

上述JSSA算法的實(shí)現(xiàn)過(guò)程即為微地震定位與震源機(jī)制的聯(lián)合反演過(guò)程。本研究以P波為標(biāo)準(zhǔn)進(jìn)行極性校正和動(dòng)校正處理，經(jīng)過(guò)處理后P波的疊加能量得到增強(qiáng)，S波和其他轉(zhuǎn)換波的疊加能量則會(huì)減弱，這樣便能有效地通過(guò)P波進(jìn)行事件的識(shí)別。

另外，由于JSSA算法需要對(duì)(x，y，z)以及(φ，δ，λ)6個(gè)參數(shù)進(jìn)行掃描，而在震源機(jī)制時(shí)空域的疊加能量道集一般有幾百萬(wàn)甚至上千萬(wàn)條，逐個(gè)計(jì)算道集采樣點(diǎn)的“亮度”函數(shù)，計(jì)算量巨大。因此，采用多次掃描的方式對(duì)JSSA算法進(jìn)行優(yōu)化，以提高其計(jì)算效率。在震源機(jī)制時(shí)空域，首先以大尺度的步長(zhǎng)掃描震源機(jī)制和空間位置；如果存在微地震事件，那么在該事件所在的大網(wǎng)格內(nèi)以小的空間步長(zhǎng)進(jìn)行二次掃描，確定其精確的空間位置；最后，以第二步確定的空間位置為中心，縮小空間掃描范圍，再以小尺度的步長(zhǎng)在縮小后的空間范圍內(nèi)掃描震源機(jī)制和空間位置。優(yōu)化加速后，在確保精度的前提下JSSA算法的計(jì)算效率提高了約400倍，雖然仍無(wú)法滿足實(shí)時(shí)反演的需求，但也足夠?qū)?shí)際資料進(jìn)行及時(shí)處理和信息反饋，達(dá)到指導(dǎo)壓裂的目的(Yu et al，2018)。對(duì)于JSSA算法，Liang等(2016)和Yu等(2018)通過(guò)一系列不同信噪比、不同震源機(jī)制、不同震源位置的理論合成數(shù)據(jù)測(cè)試驗(yàn)證了其準(zhǔn)確性，通過(guò)實(shí)際壓裂數(shù)據(jù)反演結(jié)果的分析討論，驗(yàn)證了其適用性和反演結(jié)果的合理性。

通過(guò)JSSA算法得到微地震事件的空間位置和震源機(jī)制之后，可以通過(guò)各個(gè)事件的走向φ、傾角δ和滑動(dòng)角λ計(jì)算事件的最大壓應(yīng)力軸P軸，如公式(5)～(7)所示(Kanamori et al，1974;Aki et al，1980;Stein et al，2003)。其中，地震事件的滑動(dòng)矢量為

(5)

地震事件破裂面的法向量為

(6)

基于破裂的法向量和滑動(dòng)矢量，可以計(jì)算最大壓應(yīng)力方向，即P軸

P=n-d

(7)

通過(guò)對(duì)地震事件P軸的統(tǒng)計(jì)分析，可以揭示發(fā)震斷層的應(yīng)力分布(易桂喜等，2012；楊宜海等，2015)，達(dá)到解釋壓裂施工區(qū)域應(yīng)力狀態(tài)，分析事件發(fā)震機(jī)理，從而評(píng)價(jià)壓裂效果的目的。

2 威H3-1井的實(shí)際應(yīng)用

2.1 空間定位結(jié)果分析

對(duì)于威H3-1井，結(jié)合測(cè)井?dāng)?shù)據(jù)和實(shí)際射孔數(shù)據(jù)，選取3500～6000m/s變化的等效速度對(duì)各個(gè)壓裂段的射孔事件進(jìn)行重定位，最終選取(5200±50)m/s的等效速度進(jìn)行聯(lián)合反演。對(duì)于每個(gè)壓裂段，等效速度根據(jù)各段射孔重定位結(jié)果選取，通過(guò)測(cè)試，以射孔為中心選用半徑為1500m的掃描范圍進(jìn)行聯(lián)合反演較為合理，掃描深度范圍為水平井±500m。

對(duì)于19個(gè)壓裂段，閾值Rmin為3.0時(shí)得到1507個(gè)事件，各壓裂段產(chǎn)生的事件如表2 所示，大量微地震事件產(chǎn)生在壓裂施工開(kāi)始的S01、S03和S05壓裂段。結(jié)合微地震事件的時(shí)空分布(圖3)，開(kāi)始?jí)毫央A段(S01～S06)整個(gè)掃描范圍均有微地震事件分布，且R值大于8.0的能量較大的事件主要集中在水平井附近以及水平井的西側(cè)；中期壓裂階段(S07～S13)事件數(shù)量較少，主要分布在水平井附近，深度集中在3000m左右，能量較大的事件數(shù)量較少；后期壓裂階段(S14～S19)事件數(shù)量和分布范圍均有一定程度增加，主要分布在水平井附近及水平井西南側(cè)，事件整體有向南發(fā)展趨勢(shì)，可能與掃描范圍的變化有關(guān)，深度主要在3000m以下，且能量較大的事件也主要集中在水平井西南側(cè)。

表2JSSA反演得到的事件數(shù)量

圖 3 Rmin為3.0時(shí)JSSA反演結(jié)果的x-y(a)和x-z平面圖(b)以壓裂井井口位置為原點(diǎn)；綠色線為水平井投影；灰色“+”連成的線條為地面檢波器；實(shí)心圓點(diǎn)的大小代表事件的能量大?。患t色“*”代表R值大于8.0的能量較大的事件

閾值Rmin取為3.0時(shí)，微地震事件分布較為密集，規(guī)律性不明顯，且對(duì)于R值較小、能量較弱的地震事件，其震源機(jī)制反演的誤差相對(duì)較大。因此，提高閾值Rmin為5.0，得到392個(gè)微地震事件，如圖4 所示。結(jié)合圖1 可見(jiàn)，閾值Rmin為5.0時(shí)微地震事件整體呈SN走向的條帶狀分布，與楠木寺斷層的走向及表1 中天然地震事件的走向相一致。水平井附近事件的深度集中在3000m左右，而水平井西側(cè)的大量事件則主要在3000m以下，甚至有向下延伸的趨勢(shì)，尤其是一些R值大于8.0的能量較大的事件。結(jié)合圖3(b)，R值在3.0～5.0之間能量較弱的事件大量分布在水平井附近的層位，說(shuō)明該部分事件的產(chǎn)生可能與水力壓裂活動(dòng)有關(guān)，而R值大于5.0的事件可能受深層斷裂活動(dòng)影響。

圖 4 Rmin為5.0時(shí)JSSA反演結(jié)果的x-y(a)和x-z平面圖(b)以壓裂井井口位置為原點(diǎn)； 綠色線為水平井投影； 灰色“+”線條為地面檢波器； 沙灘球的大小代表反演事件的能量大小

由聯(lián)合反演得到的微地震事件分布可以看出，水平井西側(cè)可能存在深于水平井層位的天然斷裂，呈SN走向，其可能與楠木寺斷裂有關(guān)。

2.2 震源機(jī)制結(jié)果分析

對(duì)于威H3-1井聯(lián)合反演產(chǎn)生微地震事件，選取閾值Rmin為5.0的392個(gè)事件的震源機(jī)制信息進(jìn)行分析。通過(guò)事件的震源機(jī)制信息計(jì)算出事件的P軸，統(tǒng)計(jì)P軸的走向和傾角，如圖5 所示。

圖 5 Rmin為5.0時(shí)JSSA反演結(jié)果的震源機(jī)制信息統(tǒng)計(jì)及P軸分布(a)事件震源機(jī)制信息統(tǒng)計(jì)； (b)事件的P軸分布，黑色線為事件P軸在x-y平面投影，其長(zhǎng)度代表了事件的相對(duì)大小； (c)P軸的走向和傾角統(tǒng)計(jì)

由于壓裂施工區(qū)域西北側(cè)存在SN走向的楠木寺天然走滑斷層，考慮受斷層影響的事件P軸主要呈SN和EW兩個(gè)方向，因此抽取P軸走向?yàn)镾N向和EW向分別±10°的102個(gè)事件進(jìn)行統(tǒng)計(jì)，如圖6(a)所示，由圖可見(jiàn)事件的范圍比較廣，主要集中在西側(cè)，由北至南均有分布，事件產(chǎn)生的時(shí)間主要集中在壓裂開(kāi)始階段和壓裂后期階段?？紤]水力壓裂產(chǎn)生的事件P軸主要垂直于水平井以及次生裂縫事件P軸平行于水平井，抽取平行于水平井和垂直于水平井方向±10°的87個(gè)事件統(tǒng)計(jì)分析，如圖6(b)所示，其中事件主要集中在水平井周邊及水平井頂端西北側(cè)，水平井周邊較近范圍內(nèi)事件的P軸均垂直于水平井，較遠(yuǎn)的范圍才有平行于水平井的P軸分布，證明了計(jì)算結(jié)果的合理性。較多的事件產(chǎn)生于壓裂開(kāi)始階段，主要分布在水平井西北側(cè)，而水平井周邊的事件在壓裂各個(gè)時(shí)期均存在。

圖 6 Rmin為5.0時(shí)與楠木寺斷裂和水平井作用相關(guān)的事件分布及P軸統(tǒng)計(jì)(a)P軸EW、SN向分別±10°的事件分布及P軸走向與傾角； (b)P軸垂直于水平井及平行于水平井走向的事件分布及P軸走向與傾角

對(duì)比圖6(a)、6(b)可以看出，天然裂縫的影響覆蓋整個(gè)壓裂施工區(qū)域，主要集中在水平井西側(cè)，而水力壓裂的影響則集中在水平井周邊區(qū)域，在水平井頂端西北側(cè)二者均有較多的事件分布，故該區(qū)域可能受二者共同作用的影響。通過(guò)對(duì)聯(lián)合反演事件P軸分布狀況的分析，表明了事件P軸計(jì)算結(jié)果的合理性，也驗(yàn)證了聯(lián)合反演得到的事件位置和震源機(jī)制信息的準(zhǔn)確性。

2.3 綜合分析解釋

結(jié)合圖1，對(duì)比壓裂施工區(qū)域周邊的10個(gè)天然地震事件，將坐標(biāo)調(diào)整為以威H3-1井口為原點(diǎn)，如圖7 所示，其中有3個(gè)天然地震發(fā)生在壓裂施工之后，最大震級(jí)為1.3級(jí)，另有7個(gè)天然地震發(fā)生在壓裂施工過(guò)程中，最大震級(jí)為1.5級(jí)。目前沒(méi)有充足的依據(jù)來(lái)判斷天然地震事件與壓裂事件的關(guān)系，但兩者的整體趨勢(shì)基本相同，均呈SN向條帶狀分布，結(jié)合地質(zhì)背景中楠木寺斷層的走向和位置，認(rèn)為這些天然地震事件和較大、較深的壓裂事件可能與斷層活動(dòng)相關(guān)。

圖 7 天然地震事件與聯(lián)合反演事件對(duì)比

綜合壓裂施工區(qū)域閾值Rmin為5.0的微地震事件分布、區(qū)域應(yīng)力分布以及天然地震事件分布，初步總結(jié)出壓裂施工區(qū)域存在的裂縫發(fā)育情況，如圖8 所示。其中，F(xiàn)1和F3為水力壓裂誘發(fā)的裂縫，F(xiàn)1為數(shù)條垂直于水平井的裂縫，F(xiàn)3為近SN走向，尺度較??；F2為天然裂縫，可能聯(lián)通了楠木寺斷層，由北至南發(fā)育，近SN走向，裂縫規(guī)模較大；M為水力壓裂和天然裂縫共同作用的區(qū)域，有大量微地震事件，應(yīng)力狀態(tài)復(fù)雜。

圖 8 威H3-1井裂縫F1、F2和F3的分布形態(tài)

根據(jù)震級(jí)與累計(jì)事件數(shù)目統(tǒng)計(jì)出震級(jí)與數(shù)量的線性相關(guān)(G-R關(guān)系)，進(jìn)行直線擬合，擬合直線斜率的絕對(duì)值即為b值，其能夠反映區(qū)域發(fā)生大地震的概率，b值越小則發(fā)生大地震的概率越大(Gutenberg et al，1944；胡先明等，2010)。在頁(yè)巖氣開(kāi)發(fā)中，由水力壓裂誘發(fā)的微地震事件震級(jí)較小，對(duì)應(yīng)的b值較大，而天然斷裂形成的事件震級(jí)相對(duì)較大，對(duì)應(yīng)的b值較小。在本研究中，使用的壓裂資料為歸一化后的數(shù)據(jù)，且沒(méi)有已知震級(jí)和位置的能量較大的地震事件做為模板事件，因而無(wú)法計(jì)算聯(lián)合反演事件的震級(jí)，但事件的“亮度”R值也能夠反映事件的相對(duì)大小。因此，可以在雙對(duì)數(shù)坐標(biāo)軸中統(tǒng)計(jì)“亮度”R值與累計(jì)事件數(shù)目的關(guān)系，從中得到的k值具有與b值類(lèi)似的特性，如圖9 所示。

圖 9 Rmin為3.0時(shí)威H3-1井不同區(qū)域得到事件數(shù)目的k值統(tǒng)計(jì)(a)壓裂施工區(qū)域整體k值； (b)裂縫F1和F3周?chē)?00m事件的k值； (c)裂縫F2周?chē)?00m事件的k值

對(duì)于威H3-1井采用閾值Rmin為3.0聯(lián)合反演得到1507個(gè)事件，壓裂施工區(qū)域整體的k值絕對(duì)值為3.8237，F(xiàn)1和F3周邊551個(gè)微地震事件的k值絕對(duì)值為4.0162，F(xiàn)2周邊895個(gè)微地震事件的k值絕對(duì)值為3.7534，F(xiàn)1和F3周邊事件的k值大于工區(qū)整體k值，而F2周邊事件的k值小于工區(qū)整體k值，說(shuō)明F1和F3周邊事件可能以水力壓裂誘發(fā)為主導(dǎo)因素，F(xiàn)2周邊事件可能以天然裂縫影響為主導(dǎo)因素，且F2規(guī)模較大，可能與楠木寺斷裂存在聯(lián)系，應(yīng)加強(qiáng)監(jiān)控。

3 結(jié)論

對(duì)于威H3-1井的實(shí)際數(shù)據(jù)，應(yīng)用微地震定位與震源機(jī)制聯(lián)合反演算法，閾值Rmin為3.0時(shí)得到1507個(gè)微地震事件。選用閾值Rmin為5.0的392個(gè)較大的微地震事件，統(tǒng)計(jì)其時(shí)空分布特征和震源機(jī)制信息，通過(guò)事件P軸的分布分析壓裂施工區(qū)域的應(yīng)力狀況，驗(yàn)證了聯(lián)合反演得到的事件位置和震源機(jī)制信息的準(zhǔn)確性。結(jié)合天然地震事件的統(tǒng)計(jì)結(jié)果以及各區(qū)域事件簇的k值和區(qū)域應(yīng)力反演的綜合解釋?zhuān)治龀鰤毫咽┕^(qū)域存在天然裂縫主導(dǎo)的F2裂縫、水力壓裂主導(dǎo)形成的F1和F3多條裂縫，以及應(yīng)力狀態(tài)較為復(fù)雜的M區(qū)域。

綜上所述，聯(lián)合反演算法對(duì)于威H3-1井實(shí)際資料的反演結(jié)果是準(zhǔn)確、合理的。通過(guò)聯(lián)合反演得到的微地震事件信息并綜合其他資料，能夠?qū)毫咽┕み^(guò)程中微地震事件的誘發(fā)機(jī)制進(jìn)行詳細(xì)解釋?zhuān)砻鲏毫咽┕^(qū)域水平井西側(cè)存在著規(guī)模較大的SN向天然裂縫。威H3-1井及其周邊區(qū)域構(gòu)造活動(dòng)劇烈，易于水力壓裂形成比較發(fā)育的縫網(wǎng)，在壓裂施工時(shí)應(yīng)密切關(guān)注水平井西側(cè)SN向天然裂縫的活動(dòng)情況，避免發(fā)生安全事故。