唐 杰 溫 雷 王 浩 張文征

(中國石油大學(華東)地球科學與技術學院，山東青島 266580)

1 引言

本文基于獲得的干燥和飽水正交各向異性介質(zhì)彈性參數(shù)，分析了正交各向異性介質(zhì)中的剪張源地震矩張量，研究了各向異性對雙力偶分量(DC)、各向同性分量(ISO)和補償線性偶極子分量(CLVD)的影響。

2 正交各向異性介質(zhì)中剪張源模型

2.1 含流體垂直正交裂隙介質(zhì)分析

圖1為包含垂直正交排列裂隙的正交各向異性介質(zhì)(圖1a)及其裂縫區(qū)剪切和擴張示意圖(圖1b)。研究表明，已存在裂縫介質(zhì)可能產(chǎn)生微地震活動，通常與裂縫擴張和可能發(fā)生的孔隙壓力擴散導致的應力變化有關。當流體進入已存在裂縫時可能誘發(fā)剪切破裂，若已存在裂縫貫穿水力壓裂裂縫，所導致的裂縫形變較為復雜，可能包含張開分量，會使裂縫具有復雜性，并構成水力壓裂裂縫網(wǎng)絡的組成部分[12]。

圖1 垂直正交各向異性OA介質(zhì)(a)及其裂縫區(qū)剪切和擴張(b)的示意圖

Schoenberg[13]把裂縫當成一個特殊界面處理，應力通過該界面時是連續(xù)的，而位移存在不連續(xù)現(xiàn)象，并通過引入韌度矩陣表示其相互關系。由于裂縫的存在，需在原有背景介質(zhì)的韌度矩陣上加上裂縫韌度修正項。Schoenberg模型一般用四個參數(shù)表示介質(zhì)的彈性性質(zhì)，即兩個拉梅參數(shù)(λ與μ)和裂縫法向韌度ZN及剪切韌度ZT，這大大簡化了描述裂縫介質(zhì)的參數(shù)個數(shù)。當存在多組裂縫時，可對每一組進行一定的坐標變換，然后進行線性相加(在裂縫隙密度較低時)，最后得到介質(zhì)的等效彈性參數(shù)。裂隙對有效彈性性質(zhì)的貢獻可用韌度或剛度表征，兩種不同的表達式在裂隙密度低時得到的結果一致，但在裂隙密度較高時，用韌度常數(shù)比剛度常數(shù)更準確。從物理角度看，這是因為這種情況是將裂隙看作是介質(zhì)產(chǎn)生應變的來源(裂隙比巖石基質(zhì)更易屈服)，韌度常數(shù)是各個裂隙位移的總和，可表示為

(1)

Krief等[14]認為干燥巖石(其孔隙度為γ)的體積模量K、剪切模量μ和固體基質(zhì)的體積模量Km、剪切模量μm滿足下列關系

K=Km(1-γ)3/(1-γ)

(2)

μ=μm(1-γ)3/(1-γ)

(3)

上述兩式是經(jīng)驗速度—孔隙度公式，通過它們能模擬不同孔隙度的巖石彈性特征，即背景各向同性介質(zhì)的彈性模量參數(shù)。

對于包含垂直正交排列裂隙的干燥正交對稱介質(zhì)，其剛度矩陣滿足[15]

(4)

式中：Q=λ+2μ是縱波模量，其中λ和μ分別是背景各向同性材料的拉梅參數(shù)和剪切模量，滿足λ=K-2μ/3；ΔN和ΔT分別是法向和切向裂隙弱度，滿足ΔN=ZNQ/(1+ZNQ)，ΔT=μZT/(1+μZT)，下標1和2表示兩組裂隙的參數(shù)，ZN是附加的裂隙法向韌度，ZT是附加的裂隙切向韌度； 參數(shù)r=λ/(λ+2μ),g=μ/(λ+2μ);i,j=1,2,…,6。

裂紋的切向弱度不隨流體類型變化，干燥巖石的法向裂紋弱度和切向裂紋弱度滿足[16]

(5)

式中σ是背景巖石的泊松比，這里認為切向裂紋弱度僅與泊松比相關。對于包含垂直正交排列裂隙的飽和正交對稱介質(zhì)，結合各向異性流體替換方程能夠得到飽和巖石的剛度參數(shù)[17]。

首先分析裂隙巖石的波速特征，設定基質(zhì)體積模量滿足Km=37GPa，基質(zhì)剪切模量μm=44GPa，基質(zhì)密度為2.65g/cm3，孔隙度為0.25。然后將裂隙參數(shù)加入各向同性背景中從而形成干燥OA介質(zhì)，選擇的流體為水，其密度ρfl=1.02g/cm3， 體積模量Kfl=2.7GPa，在孔隙和裂縫中加入流體形成飽和OA介質(zhì)。圖2給出了孔隙度γ=0.25、ΔT1=0.1、ΔT2=0.5時干燥和飽水正交各向異性介質(zhì)中P波、快S波和慢S波的速度變化，可見相對于干燥介質(zhì)，飽和OA介質(zhì)中各個方向的P波和S波速度表現(xiàn)出一些差異性。

圖2 干燥(上)、飽水(下)OA介質(zhì)中的P波(左)、快S波(中)及慢S波(右)速度γ=0.25， ΔT1=0.1， ΔT2=0.5; X1、X2、X3表示三個方向

圖3 干燥(a)、飽水(b)OA介質(zhì)中的剪切波分裂γ=0.25， ΔT1=0.1， ΔT2=0.5

在各向異性介質(zhì)中，當波的極化方向與各向異性方向不一致時，就會產(chǎn)生S波分裂現(xiàn)象。S波通過各向異性介質(zhì)時，沿著每一條射線路徑可分裂成兩種偏振波。它們具有不同的傳播速度和不同的偏振方向，而且在后續(xù)的各向同性介質(zhì)中傳播時可保留這種特性。圖3給出了干燥、飽水OA介質(zhì)中的剪切波分裂，可見含水后介質(zhì)中的剪切波分裂程度與干燥介質(zhì)中的相比也表現(xiàn)出一些差異性。

2.2 剪張源的源張量與矩張量

震源機制是指震源區(qū)在地震發(fā)生時的力學過程，一個剪張源是相對運動發(fā)生在斷裂面的任意滑動方向的地震源(斷層一邊的相對滑動朝向另一邊)，是水力壓裂微地震中較為合理的源模型，是微地震研究的重要機制類別。

假定斷面法向n=(n1,n2,n3)T，位錯方向v=(v1,v2,v3)T，f為v在斷層面上的投影方向，采用如圖4所示的剪張源模型，根據(jù)走向角φ、傾向角δ、滑動角θ和張裂角α的定義及相互關系，則斷面法向和位錯方向滿足[18]

(6)

(7)

圖4 剪張源模型圖示φ為走向角, δ為傾向角, θ為滑動角, α為張裂角

走向角φ、傾向角δ、滑動角θ和張裂角α通常是空間位置和時間的函數(shù)，為了簡化忽略其時空變化，研究剪張錯動的點源近似。走向角φ和傾向角δ可用于描述斷面法向方向，滑動角θ和張裂角α可用于描述位錯矢量沿著斷層面和偏離于斷層面的分量。當張裂角α>0°時，代表剪張源，對于純張裂源α=90°；α<0°，代表壓裂源，對于純壓裂源α=-90°，對于純剪切源α=0°。

剪張源的源張量滿足[19]

(8)

式中：u是位錯矢量u的大??；S為斷面面積。如果所觀測的地震波波長遠遠大于震源所涉及的尺度，則該震源可看成地震點源。

對于各向異性介質(zhì)，其矩張量和源張量之間存在如下關系[20]

Mij=cijklDkl

(9)

式中：cijkl(i，j，k，l=1，2，3)為震源區(qū)各向異性剛度參數(shù)，四階彈性張量cijkl可與二階彈性張量Cpq(p，q=1，2，…,6)相互轉化，滿足11→1，22→2，33→3，23→4，13→5及12→6。由式(9)可得干燥和飽和的OA介質(zhì)中的矩張量。

2.3 矩張量分解

微地震震源矩張量采用特征值分解法可分解為雙力偶部分(DC)、補償線性矢量偶極子成分(CLVD)以及各向同性部分(ISO)[21]

M=MISO+MCLVD+MDC

(10)

式中

各分量所占的比例滿足

(11)

式中：HISO>0代表體積膨脹，相反HISO<0代表體積壓縮；HCLVD的正負由ε控制；HDC為正。

3 正交各向異性介質(zhì)中的震源特征

下面研究各向異性參數(shù)對ISO、DC和CLVD在地震矩張量中的百分含量及P波輻射花樣影響。

3.1 非DC分量變化特征分析

矩張量分解結果能夠采用Hudson圖表示，Hudson圖(圖5)可清楚直觀地對微地震震源機制進行統(tǒng)計性分類描述[22]：首先對矩張量進行正則化分解；然后計算其特征值并按照大小進行排序，定義各向同性組分及其對應的偏量組分； 通過參數(shù)化運算，可將矩張量的6個獨立元素簡化為兩個元素(T，k)。Hudson圖中橫坐標表示剪切組分T值，縱坐標表示張性組分k值，T、k的取值范圍均為-1～1[3]。

各向異性介質(zhì)中由于剪張破裂產(chǎn)生的非DC分量依賴于破裂類型、各向異性的類型、裂縫參數(shù)以及斷面的方向等，對于OA介質(zhì)，裂隙介質(zhì)中由于剪張破裂產(chǎn)生的各分量如圖5所示，斷層面的幾何參數(shù)固定。由該圖可知對于各向同性介質(zhì)當張裂角為0°時，為純剪切源，只包含DC分量，而介質(zhì)為各向異性時，由于震源區(qū)介質(zhì)參數(shù)的影響此時矩張量中仍然包含非DC分量。研究表明較小角度的張裂角，也會產(chǎn)生較大的非雙力偶分量。此外飽和介質(zhì)中的ISO分量和DC分量相對干燥介質(zhì)而言有所增加；當走向角不沿著原有裂紋方向時會對各個分量產(chǎn)生更大的影響。

圖6a給出了背景介質(zhì)孔隙度變化時的剪張源震源矩張量分解結果，其中斷層面走向30°，傾向角為90°，滑動角為0°，張裂角為10°，可見背景介質(zhì)孔隙度變化會對剪張源矩張量產(chǎn)生影響，各向同性介質(zhì)中各分量比例不隨孔隙度的變化而變化，而各向異性介質(zhì)由于背景孔隙度的變化會影響其各向異性參數(shù)，進而會對震源矩張量產(chǎn)生影響。圖6b給出了裂隙切向弱度參數(shù)變化時的剪張源震源矩張量分解結果，可看出裂隙切向弱度參數(shù)變化會對剪張源矩張量產(chǎn)生影響。震源區(qū)各向異性介質(zhì)由于切向弱度參數(shù)的變化會影響其各向異性參數(shù)，進而會對震源矩張量產(chǎn)生影響。圖6c給出了孔隙流體為油水混合物，隨著水飽和度變化時的剪張源震源矩張量分解結果，可見介質(zhì)流體參數(shù)變化會對剪張源矩張量產(chǎn)生影響。圖6d給出了走向角φ從0°到90°時的矩張量分解結果，可見當介質(zhì)為各向同性時，走向角變化時剪張源矩張量的各個分量不會發(fā)生變化，而當震源區(qū)介質(zhì)存在各向異性時，走向角變化時剪張源矩張量的各個分量也會發(fā)生變化。

圖5 剪張源的Hudson圖(上)及矩張量分解結果(下)從左到右走向角為0°、90°和30°，傾向角為90°，滑動角為0°，張裂角為-90°～+90°；各向同性介質(zhì)(黑色)，干燥OA介質(zhì)(紅色)，飽和OA介質(zhì)(綠色)

3.2 各向異性介質(zhì)中剪張源震源機制沙灘圖

圖7為震源機制沙灘球示意圖(其中P和T軸分別位于膨脹波和壓縮波象限)，可見不同介質(zhì)中相同剪張源模型對應的震源機制在沙灘球圖中的極性分布區(qū)域表現(xiàn)出差異性，含流體介質(zhì)中的震源機制沙灘球正極性分布范圍相對較大。不同介質(zhì)中相同剪張源模型對應的震源機制在Hudson圖中的位置不同，即是不同類型介質(zhì)中震源機制所包含的各個分量存在差異。研究表明干燥與飽和OA介質(zhì)震源機制中的CLVD和ISO的比是不同的： 干燥裂紋的ISO分量較高，水飽和裂紋產(chǎn)生較高的CLVD分量，而ISO較小。此外，斷面走向也會產(chǎn)生不同的震源機制結果。以上研究表明可利用震源機制中的DC和非DC分量特征研究介質(zhì)的各向異性參數(shù)。

圖6 矩張量分解結果(a)背景介質(zhì)孔隙度為0～0.3； (b)裂隙弱度參數(shù)為0～0.3； (c)含水飽和度為0～100%；(d)走向角為0°～90°； 各向同性介質(zhì)(黑色)，干燥OA介質(zhì)(紅色)，飽和OA介質(zhì)(綠色)

圖7 震源機制沙灘球示意圖從左到右依次為各向同性、干燥和飽水正交各向異性介質(zhì)； 從上到下斷面走向依次為0°、90°和30°，傾向角為90°，滑動角為0°，張裂角為10°

3.3 各向異性介質(zhì)中剪張源的輻射花樣

依照前面地震矩張量的解析表達，然后利用Aki等[23]給出的表達式計算遠場P波輻射花樣，為了簡化起見，參考蔡曉剛等[4]做法，假定傳播介質(zhì)是全空間各向同性而震源區(qū)是各向異性的。圖8給出了各向同性和各向異性介質(zhì)條件下震源的P波三維遠場輻射圖案，可見各向同性與OA介質(zhì)中的P波輻射花樣的極性和大小表現(xiàn)出差異性，干燥OA介質(zhì)中的振幅強于各向同性中的振幅，不同介質(zhì)中的輻射花樣的分布特征總體上是一致的，此外，震源斷面的走向會產(chǎn)生不同的P波輻射花樣。以上研究表明可以利用震源機制中的P波輻射花樣特征研究介質(zhì)的各向異性參數(shù)。

圖8 P波輻射花樣從左到右依次為各向同性、干燥和飽和正交各向異性介質(zhì)； 從上到下斷層面走向依次為0°、90°和30°，傾向角為90°，滑動角為30°，張裂角為10°

4 正交各向異性介質(zhì)中微地震波場特征

利用有限差分法可研究OA介質(zhì)中不同震源機制的微地震信號傳播[24，25]，為了減小數(shù)值頻散、滿足穩(wěn)定條件和提高計算效率，網(wǎng)格和時間參數(shù)都要隨頻率改變，建立的模型網(wǎng)格間距為Δh，時間增量為Δt，滿足Δt≤0.606Δh/Vmax，其中Vmax是模型中的最大速度。計算一個三維的正交各向異性介質(zhì)模型的微地震響應，介質(zhì)參數(shù)采用前面所采用的干燥和飽和的正交裂隙介質(zhì)的參數(shù)，模型的網(wǎng)格數(shù)為100×100×100，空間采樣間隔為5m×5m×5m，震源函數(shù)為Ricker子波，主頻為30Hz，時間采樣間隔為0.5ms，震源位于模型中心(250m，250m，250m)。

圖9為干燥正交各向異性介質(zhì)中不同類型震源的波場快照，其深度為200m，上圖為水平速度分量vx，下圖為垂直速度分量vz，從左到右依次為爆炸源、純剪切源、純張裂源和剪張源(張裂角為45°)，后三種源取斷面走向為0°，可見同種介質(zhì)中不同震源類型波場快照表現(xiàn)出明顯差異，爆炸源和純張裂源的P波能量較強，而純剪切源和剪張源有較強的S波能量。

圖10是干燥和飽水正交各向異性介質(zhì)中的剪張源波場快照，其深度為200m，剪張源張裂角為10°，上圖為水平速度分量vx，下圖為垂直速度分量vz，左邊兩列為干燥正交各向異性介質(zhì)，右邊兩列為飽水各向異性介質(zhì)，圖10a、圖10c、圖10e和圖10g的斷面走向為0°，圖10b、圖10d、圖10f和圖10h的斷面走向為90°，明顯可見震源特征和介質(zhì)各向異性會對地震波極性(正負)和振幅(大小)等產(chǎn)生影響。

圖9 干燥正交各向異性介質(zhì)中的不同類型震源波場快照(a)爆炸源，vx； (b)純剪切源，vx； (c)純張裂源，vx； (d)剪張源(張裂角45°)，vx； (e)爆炸源，vz；(f)純剪切源，vz； (g)純張裂源，vz； (h)剪張源(張裂角45°)，vz； (中、右三種源)斷面走向為0°

圖10 干燥和飽水正交各向異性介質(zhì)中的剪張源波場快照(a)干燥，φ=0°，vx； (b)干燥，φ=90°，vx； (c)飽水，φ=0°，vx； (d)飽水，φ=90°，vx；(e)干燥，φ=0°，vz；(f)干燥，φ=90°，vz； (g)飽水，φ=0°，vz； (h)飽水，φ=90°，vz

5 結論

本文分析了正交各向異性介質(zhì)中的剪張源地震矩張量，研究了各向異性對雙力偶分量(DC)、各向同性分量(ISO)和補償線性偶極子分量(CLVD)的影響和各向異性介質(zhì)中P波輻射花樣和微地震信號傳播特征，得到如下主要認識和結論：

(1)分析了干燥和飽水垂直正交裂縫介質(zhì)的地震波響應特性，裂縫本身參數(shù)及其所含流體參數(shù)都會影響地震波響應特性，各向異性介質(zhì)矩張量可以通過各向異性介質(zhì)參數(shù)和源張量獲得。

(2)各向異性介質(zhì)中，由于剪張破裂產(chǎn)生的非DC分量依賴于破裂方式、各向異性介質(zhì)的類型、裂縫參數(shù)以及斷面的方向等，對于干燥裂紋，ISO分量較高，水飽和裂紋產(chǎn)生較高的CLVD分量，而ISO較小。

(3)各向同性和各向異性介質(zhì)中剪張源的遠場P波輻射花樣存在差異，震源特征和介質(zhì)各向異性對地震波走時、極性和振幅產(chǎn)生明顯影響。