水力压裂诱发的剪张型微地震震源机制矩张量反演方法

唐 杰 温 雷 李 聪 戚瑞轩

(①中国石油大学(华东)地球科学与技术学院，山东青岛 266580；②海洋国家实验室海洋矿产资源评价与探测技术功能实验室，山东青岛 266071)

0 引言

本文首先介绍利用P波初动振幅信息反演水力压裂裂缝震源机制解的方法原理，然后通过理论模型记录测试、分析方法的应用效果。

1 矩张量理论分析

1.1 源张量与矩张量

采用断面法向矢量n=(n1,n2,n3)T及位错方向矢量v=(v1,v2,v3)T表示源张量[15]

(1)

式中P为源强度。

由剪张源模型(图1)可知

(2)

(3)

走向角φ、 倾角η、 滑动角θ和张裂角α(图1)通常是空间位置和时间的函数，为了简化起见，在研究剪张错动的点源近似时可视为常量。走向角、倾角和滑动角描述了断面法向和断面位错矢量的切向，张裂角描述了位错矢量偏离断面的程度。张裂角α∈[-90°，90°]，即：当α>0°时为剪张源，对于纯张裂源，α=90°；当α<0°时，为压裂源，对于纯压裂源，α=-90°；对于纯剪切源，α=0°。

图1 剪张源模型

微地震震源机制可用矩张量的形式表示[16]

(4)

式中：M为对称矩阵，满足Mij=Mji，其中独立元素有6个，Mij表示沿着±i方向大小相等并且在j方向相距无限小的力偶(i、j=1，2，3)；M0为地震矩幅度。

对于各向同性介质而言，矩张量分量Mij和源张量分量Dij满足

Mij=λDkkδij+2μDij

(5)

式中：λ和μ为震源区介质的拉梅常数；Dkk表示下标相同的源张量分量；δ为克罗内克函数，当i=j时，δij=1，当i≠j时，δij=0。辐射花样可观测破裂面的几何分布和破裂机制类型，本质上是对沙滩球显示的拓展。图2为不同张裂角震源的纵、横波辐射花样图。

1.2 矩张量分解

对微地震震源矩张量M采用特征值分解法分解为双力偶部分MDC、补偿线性矢量偶极成分MCLVD和各向同性部分MISO[17]

M=MISO+MCLVD+MDC

(6)

其中

图2 不同张裂角震源的纵、横波辐射花样图

各分量所占的比例满足

(7)

式中Mmax为M的最大特征值。在各向同性介质中，非DC分量决定源的类型，即：当非DC分量都为正值时，源是张裂源(裂纹张开)；当非DC分量都为负值时，源是压缩源(裂纹闭合)；DC分量总是正的；当为纯剪切源时DC分量占比为100%。

1.3 矩张量图示

Hudson等[19]定义了一种震源类型的图解表示法，可实现各种破裂模式的概念化和可视化，可以研究矩张量中剪切组分和张性组分的分布情况。首先对矩张量进行正则化分解，然后计算其特征值Mj，并按大小排序

M1≥M2≥M3

(8)

(9)

(10)

通过下列参数化运算，可以计算张性组分k和剪切组分T

(11)

图3 不同张裂角(α=90°、45°、10°、0°、-10°、-45°、-90°)

横坐标表示剪切组分T值，纵坐标表示张性组分k值，两者的取值范围均为-1～1。各向同性爆炸型震源(ISO)中的膨胀震源(E)和压缩震源(I)的T-k坐标分别为(0，1)和(0，-1)；纯剪切源的坐标为(0，0)，位于图的中心位置；两类错动方向相反的补偿线性偶极子震源(CLVD+与CLVD-)的T-k坐标分别为(-1，0)和(1，0)

(12)

从而将矩张量的6个独立元素简化为Hudson图中的点(T，k)[20](图3)。

2 基于振幅信息的矩张量反演

微地震震源矩张量反演方法很多，包括采用极性、振幅以及全波形或集中考虑特定的体波相位(如P波和S波)的反演方法等。根据输入数据的质量和数量，可获得全张量或者部分张量，如果限制反演为偏量部分能够减少未知参数的数目，从而增强反演问题的稳定性[21-22]。

矩张量反演利用点源在均匀弹性介质中P波和S波的质点运动方程计算格林函数。P波和S波在x点、t时刻的振幅分别满足

(13)

(14)

式中：i为分量号；r为源到接收器的距离；ρ为密度；VP为P波速度；VS为S波速度；γi、γj和γk为从源到接收器的方位余弦；Mjk为矩张量分量；δij为克罗内克函数。

忽略震源时间函数时求解式(13)、式(14)的问题简化为

d=Gm

式中：d为包含N个振幅的矢量；G为格林矩阵，表征震源和接收器之间的振幅变化；m为地震矩矢量，包含矩张量的6个独立参数。上式的矩阵形式为

(15)

通过

m=G-1d

计算矩张量。上式为超定方程，求解需要更多的观测数据。为此，采用G-g代替G-1进行反演

G-g=[GTG]-1GT

(16)

由于实际地震数据中含有噪声，因此需要6个以上的振幅值准确求解m。在获得m后，需要稳定性分析以证明解的可靠性。反演的稳定性依赖于很多方面，如接收器的分布以及信号的信噪比等。如果要计算剪张源的角度参数，可以通过非线性反演方法获得，或者先求出Mij，再计算各个角度参数。针对震源机制的反演问题，Li等[23]和Tan等[24]研究了高频波形匹配和非线性反演方法 。

实际矩张量和反演矩张量的误差为

(17)

3 微地震矩张量反演测试分析

矩张量反演的置信度取决于多种因素。信噪比的局限性造成的输入信号振幅的不确定性、信号初动估计不准与定位不准导致的不确定性以及对穿过界面的反射(透射)损失考虑不周全、射线路径聚焦效应考虑不当、地震衰减与各向异性解释不当等都会影响矩张量反演结果。辐射花样的震源球采样对矩张量反演结果影响很大，测井资料由于井的数据观测范围有限和不规则方向孔径使矩张量反演结果不稳定。在没有先验约束条件时，利用单观测井数据通常不足以唯一地确定矩张量；当有3口观测井数据时单独使用P波或P波与S波联合反演可以获得全部矩张量；当有2口观测井数据时则须采用P波与S波联合反演才可获得全部矩张量[25]。地面监测能够在各个方向提供良好的震源球覆盖，单分量地面接收器可获得纵波振幅和初动，若没有横波信息约束，也会影响矩张量反演的稳定性。

3.1 数据误差对反演结果的影响

设定速度模型(图4a)的观测系统台站为星形排列(图4b)，震源位于地下2600m处，选取相应震源机制(图4c)，获得的MDC、MISO和MCLVD的百分比分别为61.34%、21.48%和17.18%(图4d)。图5为理论模型矩张量反演结果。由图可见：基于振幅信息的偏量解在源中没有体积变化，因此不包含MISO，仅包含MCLVD和MDC(图5c)；基于振幅信息的MDC解在源中没有体积变化和MCLVD成分，因此仅包含MDC成分(图5d)。

图4 理论模型的Hudson图

图5 理论模型矩张量反演结果

本文的反演结果由振幅信息获得，可根据反演结果绘制沙滩球，沙滩球使用纵波初动表征的震源机制图。震源机制反演可得到一组节面，其中一个为断面，另一个为辅助面。对于纯剪切源(DC)来说，这两个面是等价的，还需要借助其他辅助信息(如应力或地质信息等)确定具体节面类型。对于地面微地震监测数据处理，可以采用P波振幅信息反演震源机制。地面微地震监测具有较高的覆盖次数、较宽的方位角以及较大炮检距等特点，因此利用P波振幅信息可以得到稳定的震源机制解，文中采用全张量反演、偏量解和DC解反演方法。为了测试不同参数误差的影响，采用Monte-Carlo测试方法，对于不同误差均进行100次随机试验，采用MISO分量和MCLVD分量百分比分布图以及矩张量Hudson图与钻石图展示反演结果。根据输入数据的质量和数量，可获得全张量或者部分张量，如果限制反演为偏量部分可减少未知参数的数目、增强反演的稳定性。

受环境干扰影响，地面微地震监测资料信噪比低，考虑随机噪声对振幅的影响，首先分析极性误差和振幅误差的影响。 图6是振幅误差为20%时的全矩张量反演结果，图7是极性误差为5%时的全矩张量反演结果。对比两图可见，振幅误差导致反演结果偏差较小(图6)，极性误差导致反演结果偏差较大(图7)。采用P波振幅反演结果的MISO反演误差相对较小，而全矩张量反演的MCLVD误差相对较大，观测误差主要影响MCLVD的结果。原因为：矩张量分解为MISO和MCLVD时，MCLVD由计算最小特征值/最大特征值的比值获得，MCLVD对反演矩张量的精度更敏感；MISO和矩张量特征值之间满足线性关系，因此相对于MCLVD，MISO对振幅误差不敏感。

3.2 模型、定位精度及观测仪器对反演结果的影响

研究表明，模型参数的正确性、震源—接收器的几何分布以及定位精度等影响反演结果的可靠性，接收器的分布和选择影响反演结果的稳定性。图8为接收器20%重采样时的全矩张量反演结果。由图可见，将数据随机舍去20%，实现20%重采样，观测系统对非DC分量的反演结果非常重要，合适的观测系统可增加地震矩张量DC和非DC分量分析的可信性。

图6 振幅误差为20%时的全矩张量反演结果

图9为定位误差小于30m时的全矩张量反演结果。由图可见，当存在定位误差时，理论计算结果和反演结果的均方根误差较小，非DC分量也受影响，采用P波时误差较大，非DC分量易被虚假的DC分量所掩盖。

图10为速度参数误差小于10%时的全矩张量反演结果。由图可见，当速度参数误差较大时，采用P波振幅时非DC分量的误差较大。

图7 极性误差为5%时的全矩张量反演结果

图8 接收器20%重采样时的全矩张量反演结果

图9 定位误差小于30m时的全矩张量反演结果

图10 速度参数误差小于10%时的全矩张量反演结果

3.3 不同震源类型的反演结果分析

不同震源类型的辐射花样间存在差异，因此有限的震源球覆盖导致反演多解性。图11为不同震源类型在振幅误差为20%及极性误差为5%时的全矩张量反演结果。由图可知，张裂角α=90°时受到极性和振幅误差的影响较大，节面分布较杂乱(图11b、图11e)，与前文的结果类似，相对于MCLVD的结果而言，MISO的误差更小。

图11 不同震源类型在振幅误差为20%及极性误差为5%时的全矩张量反演结果

4 结论

本文研究了利用P波初动振幅信息反演水力压裂裂缝震源机制解的方法原理，并通过理论模型记录测试、分析了方法的应用效果。基于P波辐射花样进行微地震震源机制反演，避免了传统P波初动极性反演方法的不确定性，可以得到更可靠的压裂裂缝解释结果。获得以下认识：

(1)矩张量反演的置信度取决于多种因素，振幅误差导致反演结果偏差较小，极性误差导致反演结果偏差较大。相对于CLVD分量，采用全矩张量反演时ISO分量对误差不敏感，这是由于在矩张量分解时，ISO分量是矩张量特征值的线性函数，而CLVD分量是其非线性函数。

(2)观测系统对于非DC分量的反演结果非常重要，合适的观测系统可增加地震矩张量DC和非DC分量分析的可信度。当模型误差较大时，采用P波振幅反演时非DC分量的误差较大。

(3)辐射花样能够观测破裂面的几何分布和破裂机制类型，不同震源类型的辐射花样间存在差异，因此有限的震源球覆盖导致反演多解性。

本文利用P波初至振幅信息反演震源机制全矩张量，由于地面资料信噪比低、数据量大、观测范围大，只利用P波资料反演可以得到相对稳定的震源机制解，若利用P波、S波振幅信息联合反演可进一步提高反演的稳定性，同时在一定程度上能消除传播路径、近地表效应和仪器响应等因素对反演结果的影响。