楊心超,朱海波,崔樹(shù)果,王 瑜,龐 銳,李華昌

(1.中國(guó)石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103；2.中國(guó)石油化工股份有限公司西南油氣分公司,四川成都610094)

P波初動(dòng)震源機(jī)制解在水力壓裂微地震監(jiān)測(cè)中的應(yīng)用

楊心超1,朱海波1,崔樹(shù)果1,王 瑜1,龐 銳1,李華昌2

(1.中國(guó)石油化工股份有限公司石油物探技術(shù)研究院,江蘇南京211103；2.中國(guó)石油化工股份有限公司西南油氣分公司,四川成都610094)

微地震監(jiān)測(cè)技術(shù)是水力壓裂過(guò)程監(jiān)測(cè)及壓裂效果評(píng)價(jià)的重要手段。P波初動(dòng)震源機(jī)制解是研究震源參數(shù)的有效工具?；诟顸c(diǎn)搜索算法求解水力壓裂裂縫的P波初動(dòng)震源機(jī)制解,進(jìn)而得到其方位參數(shù)——方位角、傾角和滑動(dòng)角,為壓裂裂縫的解釋和壓裂效果評(píng)價(jià)提供有效依據(jù)。簡(jiǎn)單介紹了利用P波初動(dòng)極性信息反演水力壓裂裂縫參數(shù)的方法原理,并通過(guò)水力壓裂裂縫模型對(duì)該方法的應(yīng)用效果進(jìn)行了測(cè)試,最后將該方法應(yīng)用于四川某頁(yè)巖氣壓裂井的實(shí)際地面微地震監(jiān)測(cè)資料,得到了合理的單裂縫解釋結(jié)果。

水力壓裂；微地震監(jiān)測(cè)；震源機(jī)制；P波初動(dòng)極性

在致密氣、頁(yè)巖氣以及煤層氣等非常規(guī)油氣資源的開(kāi)發(fā)過(guò)程中,微地震監(jiān)測(cè)技術(shù)是壓裂過(guò)程監(jiān)測(cè)及壓裂效果評(píng)估的重要手段。壓裂微地震監(jiān)測(cè)技術(shù)的研究主要包含兩部分內(nèi)容:一個(gè)是震源定位技術(shù),另一個(gè)是震源機(jī)制(focal mechanism)描述技術(shù)。震源機(jī)制解對(duì)于了解壓裂區(qū)的油氣藏特征,例如天然裂縫發(fā)育特點(diǎn)、應(yīng)力狀態(tài)以及壓裂裂縫破裂機(jī)制等有著重要的指導(dǎo)作用,同時(shí)也是建立離散裂縫網(wǎng)格及估算有效壓裂體積的重要參數(shù)[1]。

對(duì)于地面微地震監(jiān)測(cè)的情況,普遍選擇利用P波信息進(jìn)行震源機(jī)制反演。盡管也可以綜合利用P波和S波信息,但是地面微地震監(jiān)測(cè)的高覆蓋次數(shù)、寬方位角以及大偏移距等特點(diǎn),使得僅利用P波信息就可以得到穩(wěn)定的震源機(jī)制解,同時(shí)避免了S波波形提取困難以及S波速度誤差對(duì)反演結(jié)果的影響[2]。P波信息又分為P波初動(dòng)極性信息和P波振幅信息,其中,P波初動(dòng)極性的物理圖像明確,是穩(wěn)定的地震波信息。在20世紀(jì)的一二十年代,人們就已經(jīng)在同一次天然地震不同臺(tái)站的記錄中發(fā)現(xiàn)了P波初動(dòng)極性的四象限分布特性[3],隨后將其與雙力偶(double-couple)力系聯(lián)系起來(lái),并在30年代發(fā)展了最初的P波初動(dòng)震源機(jī)制解求解方法。由于其簡(jiǎn)便快捷且結(jié)果較為可靠,利用P波初動(dòng)極性來(lái)求解震源機(jī)制解的方法被廣泛應(yīng)用于地震事件分析、震源參數(shù)求取以及應(yīng)力場(chǎng)研究等方面[4-6]。

已有研究結(jié)果表明,水力壓裂裂縫主要發(fā)育在天然裂縫或斷層的附近[7-8],受地下局部應(yīng)力場(chǎng)及地質(zhì)條件等因素的影響,壓裂裂縫的方向、大小及破裂速度多樣[9]。壓裂裂縫走向基本垂直或近似垂直于最小水平主應(yīng)力方向[10]；壓裂裂縫的大小和破裂速度則反映了裂縫破裂過(guò)程中應(yīng)力降的變化,體現(xiàn)在微地震信號(hào)上就是信號(hào)能量的強(qiáng)弱。受地面微地震監(jiān)測(cè)資料信噪比低的影響,利用弱微地震信號(hào)進(jìn)行震源機(jī)制研究有一定的困難,所以選取信噪比相對(duì)較高的強(qiáng)微地震事件進(jìn)行震源機(jī)制反演是目前廣泛采用的方法[11-12]。Baig等[13]指出,弱微地震事件(震級(jí)M<0)一般由地層中天然裂縫的剪切-拉伸活動(dòng)造成,而強(qiáng)微地震事件(震級(jí)M>0)則對(duì)應(yīng)于已有斷層的剪切運(yùn)動(dòng)。因此,在研究強(qiáng)微地震事件時(shí),對(duì)震源作純剪切(雙力偶)源假設(shè),并根據(jù)P波初動(dòng)極性的四象限分布現(xiàn)象來(lái)求解震源機(jī)制解是合理的。Tibi等[9]在雙力偶震源的假設(shè)下,利用P波初動(dòng)極性信息反演得到了單個(gè)微地震事件對(duì)應(yīng)壓裂裂縫的方位角、傾角、滑動(dòng)角以及其矩張量,并將一組微地震事件的矩張量求和得到了復(fù)合矩張量,用來(lái)評(píng)價(jià)壓裂裂縫網(wǎng)絡(luò)的復(fù)雜性。利用P波初動(dòng)極性信息并結(jié)合SH波極性信息,Rutledge等[8]對(duì)Cotton Valley氣田某壓裂井的壓裂裂縫進(jìn)行了震源機(jī)制分析,并將這些壓裂裂縫的成因解釋為天然裂縫沿近似最大水平主應(yīng)力方向滑動(dòng)的結(jié)果。

本文首先簡(jiǎn)單介紹利用P波初動(dòng)極性信息求解水力壓裂裂縫P波初動(dòng)震源機(jī)制解的方法原理；然后通過(guò)理論模型記錄對(duì)方法的應(yīng)用效果進(jìn)行測(cè)試；最后將該方法應(yīng)用于四川某頁(yè)巖氣壓裂井的實(shí)際地面微地震監(jiān)測(cè)資料,并依據(jù)P波初動(dòng)解對(duì)該壓裂井的壓裂裂縫進(jìn)行震源機(jī)制分析。

1 方法原理

利用P波初動(dòng)極性信息反演水力壓裂裂縫P波初動(dòng)震源機(jī)制解的理論基礎(chǔ)是[3]:斷層在地下發(fā)生純剪切錯(cuò)動(dòng)時(shí),在地面不同地點(diǎn)接收到的直達(dá)P波的初動(dòng)極性不同(圖1),其在震源球上的投影具有四象限分布的現(xiàn)象(圖2)。因此,在觀測(cè)資料充分的條件下,可以利用從微地震事件中提取的P波初動(dòng)極性在震源球上的投影來(lái)求取斷層節(jié)面解。

圖1 剪切錯(cuò)動(dòng)斷層輻射P波初動(dòng)極性圖示

圖2 P波初動(dòng)極性與斷層節(jié)面在震源球上的投影

利用P波初動(dòng)極性求解震源機(jī)制解的方法有多種,這些方法大致可分為迭代法和網(wǎng)格搜索法兩類。由于迭代算法對(duì)資料品質(zhì)的依賴性更大,所以網(wǎng)格搜索法是目前主流的求解方法[5]。本文采用了許忠淮等[4]提出的求解P波初動(dòng)震源機(jī)制解的格點(diǎn)搜索算法,對(duì)斷層節(jié)面對(duì)應(yīng)的震源坐標(biāo)架在地平坐標(biāo)架中的位置進(jìn)行全空間掃描,求取與觀測(cè)P波初動(dòng)極性符號(hào)擬合最佳而矛盾符號(hào)比ψ最小的一對(duì)正交節(jié)面(矛盾符號(hào)比定義為ψ=矛盾符號(hào)數(shù)/符號(hào)總數(shù),0≤ψ≤1)。

震源機(jī)制反演得到一組正交節(jié)面,其中一個(gè)為斷層面,另一個(gè)為輔助面。但是,對(duì)于純剪切源(雙力偶源)來(lái)說(shuō)，這兩個(gè)面是等價(jià)的。所以,僅由單獨(dú)的P波初動(dòng)震源機(jī)制解無(wú)法確定斷層面信息,還需要借由其它輔助信息來(lái)確定。例如,可以利用由測(cè)井資料等得到的地層應(yīng)力狀態(tài)來(lái)指導(dǎo)斷層面的選取；或者在得到大量微地震事件的震源機(jī)制解后,將其震源機(jī)制球按事件定位坐標(biāo)畫在同一幅圖上,相鄰微地震事件的裂縫產(chǎn)狀應(yīng)該是比較接近的,從而判斷用哪組節(jié)面解來(lái)表示裂縫的產(chǎn)狀。

2 理論模型測(cè)試

設(shè)定理論模型的觀測(cè)系統(tǒng)如圖3所示,8條測(cè)線以井口(坐標(biāo)(0,0,0))為中心呈放射狀排列,相鄰測(cè)線方位角相差45°;每條線50道接收,道間距為25m,最小偏移距為25m,最大偏移距為1250m。震源位于地下1200m深度處,其坐標(biāo)為(100m,100m,1200m)。這種檢波器放射狀布設(shè)、震源深度與最大偏移距之比約為1的觀測(cè)系統(tǒng)設(shè)計(jì)方案在水力壓裂地面微地震監(jiān)測(cè)中被普遍采用[2]。

圖3 理論模型觀測(cè)系統(tǒng)圖示

設(shè)定兩類震源模型:分別沿裂縫的走向方向滑動(dòng)(滑動(dòng)角為0)和沿傾向方向滑動(dòng)(滑動(dòng)角為90°)的壓裂裂縫模型。裂縫方位角的設(shè)定原則上是任意的。裂縫的傾角設(shè)定參考了Fisher等[14]根據(jù)大量實(shí)際資料得到的壓裂裂縫傾角與地層埋深的統(tǒng)計(jì)關(guān)系。Fisher等[14]的研究結(jié)果表明,在埋深超過(guò)1000m時(shí),由于上覆地層壓力的影響,壓裂裂縫基本表現(xiàn)為高角度裂縫。按照以上原則,設(shè)計(jì)了2個(gè)壓裂裂縫模型:模型1,方位角為30°,傾角為70°,滑動(dòng)角為0；模型2,方位角為30°,傾角為70°,滑動(dòng)角為90°。

采用Chen[15]提出的基于廣義反透射系數(shù)的震源模擬方法和圖4所示的水平層狀速度模型,對(duì)兩個(gè)壓裂裂縫模型進(jìn)行正演模擬,得到的理論微地震記錄如圖5所示。

圖4 水平層狀速度模型

從圖5理論微地震記錄中提取初至P波的初動(dòng)極性,利用格點(diǎn)搜索算法(搜索步長(zhǎng)為1°×1°×1°)進(jìn)行壓裂裂縫模型的參數(shù)反演,反演結(jié)果與實(shí)際模型參數(shù)的對(duì)比結(jié)果如表1所示。圖6給出了反演得到的節(jié)面解與P波初動(dòng)極性在震源球上的赤平投影。圖6中,藍(lán)色圓圈代表正極性(P波初至運(yùn)動(dòng)方向向上),紅色圓圈代表負(fù)極性(P波初至運(yùn)動(dòng)方向向下),兩條黑色曲線代表反演得到的兩個(gè)節(jié)面解。

從模型1的反演結(jié)果(表1)可以看出,對(duì)于沿走向方向滑動(dòng)開(kāi)裂類型的壓裂裂縫,由于其在實(shí)驗(yàn)中的觀測(cè)系統(tǒng)范圍內(nèi)表現(xiàn)出了明顯的“四象限分布”現(xiàn)象,所以使用P波初動(dòng)極性信息可以很好地反演出該壓裂裂縫模型的參數(shù)信息——方位角、傾角和滑動(dòng)角。而對(duì)于模型2沿傾向方向滑動(dòng)開(kāi)裂類型的壓裂裂縫,由于受觀測(cè)范圍的限制,利用P波初動(dòng)極性信息只能準(zhǔn)確確定一個(gè)節(jié)面的位置,而無(wú)法確定另外一個(gè)節(jié)面的位置,反映在裂縫參數(shù)上則表現(xiàn)為可以準(zhǔn)確反演壓裂裂縫的方位角和傾角,而無(wú)法準(zhǔn)確反演滑動(dòng)角。

圖7給出的是模型2對(duì)應(yīng)的理論節(jié)面解與本次實(shí)驗(yàn)中得到的P波初動(dòng)極性在震源球上的赤平投影,圖中藍(lán)色虛線代表了準(zhǔn)確反演出所有模型參數(shù)所需要的最小額外觀測(cè)范圍。而在實(shí)際生產(chǎn)工作中,受采集成本的限制,如此大的觀測(cè)范圍是無(wú)法滿足的,這也是P波初動(dòng)極性反演在應(yīng)用于水力壓裂地面微地震監(jiān)測(cè)時(shí)的局限性。可是即便如此,可靠的方位角和傾角信息仍然可以對(duì)壓裂裂縫的解釋工作發(fā)揮非常大的指導(dǎo)和幫助作用。

圖5 壓裂裂縫模型1(a)和模型2(b)的正演模擬微地震記錄

表1 模型參數(shù)及反演結(jié)果

模型震源(裂縫)類型模型參數(shù)反演結(jié)果方位角/(°)傾角/(°)滑動(dòng)角/(°)方位角/(°)傾角/(°)滑動(dòng)角/(°)矛盾符號(hào)比1走滑30．0070．00028．9071．000．950．04122傾滑30．0070．0090．0030．3070．00121．100．1138

圖6 壓裂裂縫模型1(a)和模型2(b)模擬數(shù)據(jù)反演節(jié)面解

圖7 模型2的理論節(jié)面解

3 實(shí)際資料應(yīng)用試驗(yàn)

選取了四川某頁(yè)巖氣壓裂井地面微地震監(jiān)測(cè)資料進(jìn)行本文方法的應(yīng)用試驗(yàn)。該次監(jiān)測(cè)所采用的觀測(cè)系統(tǒng)如圖8所示。圖8中,10條測(cè)線以井口(紅點(diǎn))為中心呈放射狀排列,相鄰測(cè)線方位角相差36°；每條線125道接收,道間距為25m,最小偏移距為300m,最大偏移距為3400m。兩個(gè)射孔段深度分別為3294～3326m和3351～3369m。從監(jiān)測(cè)數(shù)據(jù)中共檢測(cè)到1073個(gè)微地震事件,圖9給出了其水平方向定位結(jié)果,可以看到,壓裂裂縫網(wǎng)絡(luò)主要沿北東—南西方向擴(kuò)展發(fā)育。從所有微地震事件中選取116個(gè)高信噪比的強(qiáng)微地震事件進(jìn)行震源機(jī)制分析。

圖8 實(shí)際地面微地震監(jiān)測(cè)觀測(cè)系統(tǒng)平面示意圖解

圖9 實(shí)際監(jiān)測(cè)資料微地震事件水平方向定位結(jié)果 (圓圈大小代表事件能量相對(duì)大小)

從選取的強(qiáng)微地震事件中觀察到了兩類震源機(jī)制,分別記為FM1和FM2。圖10給出了分屬FM1類和FM2類的兩個(gè)微地震事件的原始記錄。利用記錄上提取的P波初動(dòng)極性信息,基于格點(diǎn)嘗試法分別求解了兩個(gè)事件的P波初動(dòng)震源機(jī)制解。反演時(shí)使用了圖4所示的水平層狀速度模型,該速度模型由壓裂井的聲波測(cè)井資料求得,并利用射孔資料進(jìn)行了校正。反演結(jié)果如表2所示。圖11 給出了兩個(gè)事件各自的節(jié)面解和P波初動(dòng)極性在震源球上的赤平投影。

如前所述,僅根據(jù)單獨(dú)微地震事件的P波初動(dòng)解無(wú)法確定壓裂裂縫的產(chǎn)狀信息,故求取了所有116個(gè)強(qiáng)微地震事件的P波初動(dòng)解,并將反演結(jié)果的沙灘球顯示與所有微地震事件的水平方向定位結(jié)果聯(lián)合顯示(圖12)。從圖12可以看出,近似北向和北偏東方向的節(jié)面普遍存在于各個(gè)強(qiáng)微地震事件的P波初動(dòng)節(jié)面解中。所以,按照相鄰壓裂裂縫的產(chǎn)狀相近的原則,可以選取表2中“節(jié)面1”的角度數(shù)據(jù)分別作為這兩個(gè)壓裂裂縫的產(chǎn)狀信息。這也說(shuō)明,作為在這次地面微地震監(jiān)測(cè)中觀測(cè)到的兩類主要的裂縫類型,FM1和FM2分別代表了沿近似北向和北偏東方向發(fā)育的剪切型裂縫,這一點(diǎn)與通過(guò)定位結(jié)果得到的壓裂裂縫網(wǎng)絡(luò)的擴(kuò)展方向相符。另外,裂縫的傾角特點(diǎn)——高角度,也與Fisher等[14]的研究結(jié)論一致。

圖10 實(shí)際地面監(jiān)測(cè)資料中分屬FM1類(a)和FM2類(b)微地震事件的原始記錄

表2 實(shí)際地面微地震監(jiān)測(cè)資料兩類震源機(jī)制反演結(jié)果

震源機(jī)制類型節(jié)面1節(jié)面2方位角/(°)傾角/(°)滑動(dòng)角/(°)方位角/(°)傾角/(°)滑動(dòng)角/(°)矛盾符號(hào)比FM13．8067．3099．76159．7724．6167．930．1028FM216．9466．90129．72132．2144．9633．710．0788

圖11 實(shí)際地面微地震監(jiān)測(cè)資料FM1(a)和FM2(b)兩類震源機(jī)制類型的節(jié)面解

圖12 實(shí)際地面微地震監(jiān)測(cè)資料震源機(jī)制反演結(jié)果與所有微地震事件的水平方向定位結(jié)果(紅點(diǎn))聯(lián)合顯示

4 結(jié)束語(yǔ)

我們介紹了P波初動(dòng)震源機(jī)制解用于水力壓裂裂縫參數(shù)反演的方法原理,并使用典型水力壓裂裂縫模型的正演模擬記錄對(duì)方法進(jìn)行了測(cè)試。測(cè)試結(jié)果表明,通過(guò)P波初動(dòng)解能夠得到壓裂裂縫方位角和傾角的準(zhǔn)確信息,而對(duì)滑動(dòng)角的反演精度則會(huì)受到觀測(cè)系統(tǒng)的制約。最后將該方法應(yīng)用于四川某頁(yè)巖氣壓裂井的地面微地震監(jiān)測(cè)數(shù)據(jù),得到了合理的單裂縫解釋結(jié)果。

[1] Eisner L,Williams-Stroud S,Hill A,et al.Beyond the dots in the box:microseismicity-constrained fracture models for reservoir simulation[J].The Leading Edge,2010,29(3):326-333

[2] Duncan P M,Eisner L.Reservoir characterization using surface microseismic monitoring[J].Geophysics,2010,75(5):139-146

[3] 安藝敬一,理查茲.定量地震學(xué):理論和方法[M].李欽祖,鄒其嘉,譯.北京:地震出版社,1980:85-92 Aki K,Richards P G.Quantitative seismology:theory and methods[M].Li Q Z,Zhou Q J,translators.Beijing:Seismological Press,1980:85-92

[4] 許忠淮,閻明,趙仲和.由多個(gè)小地震推斷的華北地區(qū)構(gòu)造應(yīng)力場(chǎng)的方向[J].地震學(xué)報(bào),1983,5(3):268-279 Xu Z H,Yan M,Zhao Z H.Evaluation of the direction of tectonic stress in north China from recorded data of a large number of small earthquakes[J].Acta Seismologica Sinica,1983,5(3):268-279

[5] 俞春泉,陶開(kāi),崔效鋒,等.用格點(diǎn)嘗試法求解P波初動(dòng)震源機(jī)制解及解的質(zhì)量評(píng)價(jià)[J].地球物理學(xué)報(bào),2009,52(5):1402-1411 Yu C Q,Tao K,Cui X F,et al.P-wave first-motion focal mechanism solutions and their quality evaluation[J].Chinese Journal of Geophysics,2009,52(5):1402-1411

[6] 胡幸平,俞春泉,陶開(kāi),等.利用P波初動(dòng)資料求解汶川地震及其強(qiáng)余震震源機(jī)制解[J].地球物理學(xué)報(bào),2008,51(6):1711-1718 Hu X P,Yu C Q,Tao K,et al.Focal mechanism solutions of Wenchuan earthquake and its strong aftershocks obtained from initial P wave polarity analysis[J].Chinese Journal of Geophysics,2008,51(6):1711-1718

[7] Iannacchione A T,Batchler T,Marshall T E.Mapping hazards with microseismic technology to anticipate roof falls—a case study[J].23rdInternational Conference on Ground Control in Mining,2004,327-333

[8] Rutledge J T,Phillips W S.Hydraulic stimulation of natural fractures as revealed by induced microearthquakes,Carthage Cotton Valley gas field,east Texas[J].Geophysics,2003,68(2):441-452

[9] Tibi R,Vermilye J,Lacazette A,et al.Assesment of hydraulic complexity and stress field variability in an unconventional reservoir from composite moment tensor of double-couple microseismic events[J].Expanded Abstracts of 83rdAnnual Internat SEG Mtg,2013,2183-2187

[10] Refunjol X E,Marfurt K J,Le Calvez J H.Inversion and attribute-assisted hydraulically induced microseismic fracture characterization in the North Texas Barnett Shales[J].The Leading Edge,2011,30(3):292-299

[11] Zhou R,Huang G D,Snelling P,et al.Magnitude calibration for microseismic events from hydraulic fracture monitoring[J].Expanded Abstracts of 83rdAnnual Internat SEG Mtg,2013,2145-2149

[12] Staněk F,Eisner L.New model explaining inverted source mechanisms of microseismic events induced by hydraulic fracturing[J].Expanded Abstracts of 83rdAnnual Internat SEG Mtg,2013,2201-2205

[13] Baig A,Urbancic T,Viegas G,et al.Can Small events (M<0) observed during hydraulic fracture stimulations initiate large events (M>0)[J].The Leading Edge,2012,31(12):1470-1474

[14] Fisher K,Warpinski N.Hydraulic-fracture-height growth:real data[J].Expanded Abstracts of 86thSPE Annual Conference,2011,SPE145949

[15] Chen X F.Seismogram synthesis in multi-layered half space,part 1:theoretical formulation[J].Earthquake Research in China,1999,13(2):149-175

(編輯：顧石慶)

Application of P-wave first-motion focal mechanism solutions in microseismic monitoring for hydraulic fracturing

Yang Xinchao1,Zhu Haibo1,Cui Shuguo1,Wang Yu1,Pang Rui1,Li Huachang2

(1.SinopecGeophysicalResearchInstitute,Nanjing211103,China;2.SinopecSouthwestOilCompany,Chengdu610094,China)

Microseismic monitoring technology is very important for hydraulic fracturing monitoring and hydraulic fracturing response evaluation.P-wave first-motion focal mechanism solution is an efficient tool for focal parameters study.P-wave first-motion focal mechanism solutions for hydraulic fracture are derived with grid search method.The fracture azimuthal parameters that are strike angle,dip angle and slip angle are calculated to serve the hydraulic fracture interpretation and hydraulic fracturing response evaluation.After a brief introduction to the principles of hydraulic fracturing parameters inversion based on P-wave first-motion polarity information,the reliability of the method is demonstrated by numerical simulation data of hydraulic fracturing models.The method is applied on surface microseismic data acquired during the hydraulic fracturing of some shale gas reservoir in Sichuan province,and reasonable hydraulic fracture interpretation result is achieved based on the P-wave first-motion solutions.

hydraulic fracturing,microseismic monitor,focal mechanism,P-wave first-motion polarity

2014-05-09;改回日期：2014-12-26。

楊心超(1986—),男,助理工程師,主要從事水力壓裂微地震監(jiān)測(cè)相關(guān)技術(shù)研究。

P631

A

1000-1441(2015)01-0043-08

10.3969/j.issn.1000-1441.2015.01.006