APP下载

含衰减地层微地震震源机制反演及其反演分辨率

2016-09-29翟鸿宇常旭王一博姚振兴

地球物理学报 2016年8期
关键词:张量震源反演

翟鸿宇,常旭,王一博,姚振兴

1 中国科学院地质与地球物理研究所 中国科学院页岩气与地质工程重点实验室,北京 100029 2 中国科学院大学,北京 100049



含衰减地层微地震震源机制反演及其反演分辨率

翟鸿宇1,2,常旭1*,王一博1,姚振兴1

1 中国科学院地质与地球物理研究所 中国科学院页岩气与地质工程重点实验室,北京100029 2 中国科学院大学,北京100049

微地震震源机制的反演对于非常规油气开发具有至关重要的作用.微地震信号主频高、能量小,容易受地层吸收衰减作用的影响使其波形发生畸变,本文提出了一种考虑地层吸收衰减作用的微震源机制反演方法,并利用费雷谢偏导矩阵的SVD分解(特征值分解)方法,分析研究了地层的吸收衰减因子的变化对于微地震震源机制反演分辨率的影响,根据理论计算给出了不同地震数据对各种震源机制反演的适用条件.理论计算证明,采用直达P波和S波数据联合反演震源的T值,单独利用直达P波反演震源的k值,可以有效降低地层吸收衰减作用对反演结果的影响.

地层吸收衰减;微地震;震源机制反演;反演分辨率

1 引言

众多研究者已经对微地震监测开展了广泛的研究,通过分析处理水力压裂过程中产生的微震信号,我们可以反演得到微地震发生的震源位置、震级大小、震源机制等信息,从而得到储层中裂缝的发育情况.这些信息为最终的储层改造评价提供了有力证据(Eisner et al.,2010,2011;Kendall et al.,2011;Maxwell,2011).

尽管前人提出了不同的震源机制反演方法并对反演过程中存在的一些问题进行了分析,但地层的吸收衰减作用对震源机制反演的影响并没有得到充分的研究.并且微震信号具有主频高,能量小的特性,地层的吸收衰减必然会严重影响微震信号的波形、频谱和信噪比,从而降低震源机制的反演精度.因此,本文首先根据地层的吸收衰减特性建立了含衰减地层的微地震震源机制频率域反演方法.然后,通过数值计算,并利用反演分辨率矩阵分析了地层吸收衰减的计算误差对不同震源机制反演结果的影响.

2 震矩张量与微地震波场关系

根据Backus和Mulcahy(1976a,1976b)最早对震源机制的定义,任何震源都是由于地层中的过量应力σ造成的,通过对应力过量变化率进行时间与空间上的积分,可以用地震矩对震源进行定义:

(1)

并且,这一严格定义同样可以用震矩张量M的形式进行表示:

(2)

其中,M0表示地震矩幅度,每一个元素Mij表示沿着±i方向,大小相等,并且在j方向相距无限小的一对力所形成的力偶(i,j=1,2,3).由于任何地震的震矩张量必须满足净应力为零和净应力矩也为零的基本条件(Udías et al.,2014).因此,震矩张量M中的元素必须满足Mij=Mji这一条件,即,M是一个对称矩阵,只有6个独立的元素(M11、M12、M13、M22、M23、M33).

利用震矩张量,我们可以对岩石破裂的几种基本形式进行定量的描述.例如:

补偿线性偶极子震源(CLVD)

以上定义中,只是部分震源机制的表示形式.对于同一类型的震源,由于所受的力偶的不唯一性,其震矩张量的表示形式也有多样.如双偶极子震源还可以表示为:

但是,当实际地层中的岩石发生破裂时,其诱发的微地震震源是十分复杂的.利用正则化分解,可以将任意的震源机制分解为三种基本震源机制(ISO、CLVD、DC震源)的组合(Julian et al.,1998;Chapman and Leaney,2012),即

(3)

Aki和Richards(2002)将各向同性介质中,点震源的远场波场函数与震源机制的关系定义为:

(4)

(5)

(6)

其中,向量n=(ni,nj,nk)表示震源到检波点的方向向量,δij为克罗内克函数,Mjk表示震矩张量中的元素.

3 地层的吸收衰减效应

地震波在真实地层中传播的过程中,其能量不仅会发生几何扩散,还会由于地层介质颗粒间的内摩擦作用,孔隙流体的流动,黏弹性介质的黏性张弛等作用而产生吸收衰减和地震波频散现象.而对于地层的吸收衰减效应主要使用地层的品质因子Q值进行表征,其严格定义可表示为:

(7)

其中,E0表示地震波的原始能量,ΔE表示地震波传播一个周期后所消散的能量(Tonn,1991).在常规的地震勘探中,由于地震子波的频带范围有限,我们假设地层的吸收衰减作用与地震波的频率无关,即地层的Q值在常规地震勘探频带范围内保持不变(Zhang and Ulrych,2002).因此,地层对地震波的衰减过程可以表示为:

(8)

其中,A(f)表示检波器接收到的波形频谱,S(f)表示震源的频谱特征,t表示炮点到检波点的走时.

但是,微地震信号的主频值和频带范围都远远高于常规地震勘探的地震信号特征.并且地层Q值也表现出频率依赖性,对于不同类型的介质和裂缝分布特征,地层的QP和QS在不同频带范围内具有较大的差异,且随着频率具有不同的变化特征(Pride et al.,2004;Brajanovski et al.,2006;Müller et al.,2010;Quintal et al.,2014).实际的应用中,为了分析和计算的方便,可以假设在微地震事件的有效频带范围内,地层的Q值仍保持其频率独立性(Eaton et al.,2014).因此,在反演的过程中本文仅利用微地震波形有效频带内的数据作为反演过程的输入数据.根据公式(4)—(6)和公式(8),我们可以得到,简化后的受到地层吸收衰减作用影响的波场函数表达式:

(9)

根据(9)式,本文正演出主能量频带范围受到地层吸收衰减作用影响的微地震波场,以此模拟微地震事件在地层中的传播过程.

4 微震源机制反演与解释

在均匀各向同性介质中,根据公式(9)的定义,震源机制与地震波场函数之间的关系可以用矩阵的形式表示为:

(10)

(11)

(12)

另一方面,根据波场数据与震矩张量之间的关系可知,公式(10)是一个线性不适定方程组.为了求解该方程组,我们利用MATLAB中的最小二乘算法,即,通过求取观测数据与正演数据间的最小平方差Γ(公式(13)),便可以反演得到微地震震源机制Mjk.

(13)

在实际的计算过程中,本文通过傅里叶变换将时间域数据转换到频率域进行计算,以提高反演的速度.

对于反演得到的震矩张量,我们将采用Hudson图示法的形式进行解释分析.这一图示法最早由Hudson等(1989)提出,对于任意的震矩张量进行正则化分解,首先计算其特征值并按照大小进行排序:

(14)

(15)

(16)

(17)

剪切组分T为:

(18)

通过这一参数化运算,震矩张量的6个独立元素可以简化为两个元素(T、k)进行表示(如图 1).其中,横坐标表示剪切组分T值,纵坐标表示张性组分k值,T、k的取值范围均为-1~1.通过该图我们可以清楚直观的对水力压裂过程中发生的全部微地震震源机制进行统计性的分类描述.

不同的震源机制位于图 1中不同的T-k坐标位置处.例如,各向同性爆炸型震源张量(MISO)中的膨胀震源(E)和压缩震源(I)的T-k坐标分别为(0,1)和(0,-1),而纯剪切型震源张量(MDC),即双偶极子震源的T-k坐标为(0,0),位于Hudson图示的中心位置.两类错动方向相反的补偿线性偶极子震源(M+CLVD与M-CLVD)的T-k坐标则分别为(-1,0)和(1,0).

图1 Hudson图示法Fig.1 Hudson diagram

5 数值计算分析

为了分析地层的吸收衰减作用对微震震源机制反演的影响,本文设计了一个均匀各向同性的地质模型.三口垂直观测井位于距离震源点距离为150 m的圆弧上,每口井中模拟9个三分量检波器进行观测(如图 2).在这一观测系统下,所获得的数据具有充足的观测方位角,可以避免因观测方位角不足而带来的反演误差.另一方面,该模型的各项物性参数如表 1所示,其中,纵横波速度与地层密度参考延长组地层的速度和密度特征,地层的Q值则参考前人的研究成果进行设定,在本文的所有数值计算中均采用600 Hz的雷克子波作为震源子波;考虑到实际地层Q值的复杂性,模型中对其设定了一个较大的误差范围-90%~165%.本文主要分析在不同地层Q值误差条件下,微震源机制反演结果变化情况.

根据图2的地质模型与观测系统,我们利用DC型震源进行数值计算,得到了直达P波和S波的三分量微地震记录(如图 3).由于DC型震源机制具有特定能量辐射样式,三口井中所观测到的地震记录具有明显的差异,这与常规地震勘探中利用假设爆炸震源所获得的地震记录有本质的不同.

在数值计算过程中,本文利用不同类型微震源正演得到的微地震数据作为观测数据,然后根据公式(10)反演对应的微震源机制.对于不同的观测数据,反演结果的精度也不尽相同,这是根据不适定方程(10)中的系数矩阵的条件数决定的(常旭等,1999).为分析不同的Q值误差对微震源机制反演的影响,本文采用了Grechka(2010,2015)提出的费雷谢偏导矩阵SVD分解方法,对震矩张量M每个元素的反演分辨率进行了详细分析.费雷谢偏导矩阵SVD分析方法首先求取观测数据关于震矩张量6个独立元素Mij的偏导数矩阵F,通过对矩阵F进行SVD分解:

表1 模型参数与Q值的误差范围

图2 观测系统(a)观测系统俯视图;(b)单口观测井与震源位置示意图.Fig.2 Observation system(a)Downward view of observation system;(b)Position relationship of single well and the source.

图3 DC类型微震源三分量地震记录Fig.3 Three-component microseismic record of DC source

(19)

可以得到分解后的奇异值矩阵Λ和特征向量矩阵W,对比分析这两个矩阵我们可以获得不同观测数据对于震矩张量M的6个独立元素的反演分辨率(如图 4).其中,图 4a和图 4b表示在地层Q值没有误差的情况下,分别利用P波和S波联合反演和单独利用P波反演的震源机制反演分辨率矩阵.

当获得的地层Q值没有误差,且观测方位角充足时,无论利用直达P波和S波数据联合反演,还是单独利用直达P波进行反演,均可以得到准确的反演结果.但是,不同类型的输入数据对于震矩张量的不同元素的反演敏感度是不相同的.从图 4中可以看出,无论使用何种类型数据进行反演,经过费雷谢偏导矩阵SVD分解后,其对应的奇异值Λi均远大于零,这验证了当Q值没有误差时可以通过当前观测系统反演得到震矩张量.但是,对比图 4a和图 4b中的特征矩阵W,当利用P波和S波数据联合反演时,待反演震矩张量中的M11、M12和M13所对应的特征向量具有更好的稀疏性;而只利用直达P波进行反演时,输入数据的约束效果刚好相反,即M22、M23和M33所对应的特征向量具有更好的稀疏性.

然而,在实际的微地震反演过程中,模型的各项参数总是存在不同幅度的误差.尤其是地层的Q值计算误差,对于具有高主频的微地震信号更容易产生较大的干扰.因此,本文在数值计算的反演过程中,假定地层的品质因子受到不同程度的误差干扰,而震源子波和地层的其他物性参数均不含有误差,以此重点分析地层不同幅度的Q值衰减误差对震源机制反演结果的影响程度.

5.1随机微震源机制反演

在本小节中,我们利用423个不同类型的随机微震源作为输入的微地震震源(如图 5),并根据表 1中的参数所示,将地层的品质因子误差引入反演过程中.然后,采用直达P波和S波数据联合反演得到了受到不同Q值误差影响后的微震震源机制反演结果.

图6展示了当地层Q值误差从-90%逐渐增大到165%的过程中,随机微震源的部分反演结果.

图4 微震源机制反演分辨率分析Fig.4 Inversion resolution of microseismic focal mechanism

图5 原始输入随机震矩张量Fig.5 Initial input random moment tensors

从图中可以看出地层的品质因子误差对震源机制的反演结果产生了严重干扰.并且,当地层的品质因子具有正向误差时,随机微震源的反演结果出现了向ISO型震源汇聚的现象;而当地层的品质因子具有负向误差时,随机微震源反演结果开始向TC型震源汇聚.由此可以说明,当我们过高地估计地层的品质因子时,最终的反演结果更偏向于爆炸型震源,而当我们低估地层的品质因子时,最终的反演结果更偏向于TC型震源.并且,对比正向Q值误差与负向Q值误差对应的反演结果可以发现,地层品质因子的负向误差对微地震震源机制反演的影响更大.

但是,不同类型的微震源对于地层品质因子误差的抗干扰能力并不相同,而利用随机微震源,我们不能确定地层品质因子的误差对于各种类型微震源的震源机制反演结果的影响程度.因此,本文设计了第二个数值算例,用于分析特定类型的张性震源和剪切震源在地层吸收衰减效应影响下的反演结果误差.

5.2含随机扰动的四类基本微震源反演

在实际的水力压裂地层中,微震源并不是单一类型的震源机制组成.根据公式(3),任意的震源机制都可以利用四类基本震源机制(ISO、DC、CLVD-和CLVD+型震源)进行描述.为了模拟实际地层中复杂震源机制,本小节对四类基本震源引入一定组分含量的随机扰动震源项,即:

(20)

其中,MST表示基础微震源,MPer表示随机扰动微震源项,φ表示扰动震源项所占的组分系数(本文计算中,φ值为30%).然后,利用图2的地质模型和观测系统正演出每个微震源对应的直达P波与S波地震记录,并作为反演过程的输入数据.

通过数值计算,分析地层Q值存在误差时不同数据类型微震源反演结果的影响(如图 7).其中,图 7a、7b、7c和7d分别表示利用P波和S波数据联合反演含随机扰动ISO、DC、CLVD-和CLVD+类型震源的结果,图 7e、7f、7g和7h则表示单独利用P波数据进行反演的结果.图中不同颜色代表不同地层Q值误差情况下的反演结果,当地层Q值没有误差时,真实的反演结果位于图中的“*”号位置.根据图中反演结果的变化情况,可以得出以下结论:

(1)利用直达P波和S波联合反演的结果更容易受到地层Q值误差的影响,其反演结果主要表现在张性组分的误差变化,地层Q值减小使得所有震源的反演结果向爆炸型震源偏移,地层Q值增大则产生相反的偏移效果;

(2)当单独利用P波数据进行反演时,其反演结果整体好于直达P波和S波联合反演的结果,但地层Q值减小对震源机制反演的误差影响更大;

(3)不同类型的震源对于地层Q值误差的敏感度不同,地层的Q值减小对于反演结果的影响要大于地层Q值增大带来的干扰.

为了定量地分析反演结果与反演误差,我们计算了四种类型随机震源机制反演结果的剪切组分(T值)与张性组分(k值)误差(如图8).通过对比可以看出,利用不同数据进行震源机制反演所得到的T值与k值反演误差具有明显不同的变化规律.

当地层Q值的正向误差逐渐增大,即地层Q值逐渐增大时,利用直达P波和S波联合反演的震源机制T值结果受到地层Q值误差的影响很小,而其k值反演结果对于地层Q值误差的增大很敏感;然而,利用直达P波单独反演时,不同震源类型的反演误差变化很大;其中,ISO和DC型随机震源的k值反演误差较小,T值反演误差较大,而CLVD-和CLVD+型随机震源的T值误差和k值误差均随地层Q值误差的增大而增大,两者对比可知,CLVD-型随机震源的T-k值误差的增大速率更低.

另一方面,当地层Q值的负向误差逐渐增大,即地层Q值逐渐减小时,ISO型随机震源的T-k值的反演误差出现了不同的变化规律,即单独利用直达P波反演得到的误差明显小于直达P波和S波数据联合反演的误差.DC、CLVD-和CLVD+型随机震源的T-k反演误差随Q值误差增大的速率明显大于单独利用P波反演的结果;并且,对于T值,利用直达P波和S波反演效果较好,而对于k值,单独利用直达P波反演效果较好.

图6 随机微震源反演结果对比Fig.6 Comparison of inversion results for 423 random micro earthquake sources

图7 不同Q值误差条件下微震源机制反演结果对比Fig.7 Comparison of microseismic focal mechanism inversion results under different errors of Q factor

图8 不同Q值误差条件下震源机制反演误差对比Fig.8 Comparison of microseismic focal mechanism inversion errors under different errors of Q factor

图9 含Q值误差的震源机制反演分辨率对比Fig.9 Comparison of microseismic focal mechanism inversion resolutions with Q factor errors

综合上述所有类型震源机制的反演结果和其对应的反演误差,结合各类震源机制对于地层Q值误差的敏感度,本文认为在震源机制反演的过程中,可以利用直达P波和S波数据联合反演震源的T值,单独利用直达P波反演震源的k值,从而提高反演算法的稳定性和抗干扰性.

5.3震源机制反演分辨率分析

为了深入地探索地层Q值误差对震源机制反演的影响,本小节同样采用费雷谢偏导矩阵SVD分解方法对不同Q值条件下的震源机制反演分辨率进行了对比分析(如图 9).其中,图 9a、9c为利用直达P波和S波联合反演时,地层分别具有-60%和150%误差时的反演分辨率;图 9b、9d为利用直达P波单独反演时,地层分别具有-60%和150%误差时的反演分辨率.将图9与图4对比分析可知,在反演过程中不同观测数据对地层Q值误差的敏感度是不同的.

当利用直达P波和S波数据联合反演时,地层Q值负向误差逐渐增大(如图 9a),即地层衰减效果增强使得奇异值Λi的分布稀疏性降低,这会影响最终解的稳定性.另外,奇异值Λi所对应的特征向量Wi也发生了改变,这使得观测数据对M11、M12和M13的反演约束性降低,从而导致反演结果的误差增大.另一方面,随着地层Q值正向误差逐渐增大(如图 9c),即地层衰减效果减弱时,费雷谢偏导矩阵对应的奇异值Λi稀疏性变化不大,而特征向量W5和W6的独立性明显下降,这使得从观测数据中无法求解出准确的M23和M33.

当利用直达P波单独反演时,无论地层Q值的误差正向增大还是负向增大,费雷谢偏导矩阵对应的特征向量矩阵Wi没有发生明显变化,这说明观测数据对于震源机制的约束性没有改变.但是,当地层Q值的误差负向增大时,奇异值Λi的稀疏性降低,最终导致利用直达P波反演时的误差小于直达P波和S波联合反演的误差,且对于地层Q值的正向误差不敏感.由此可知,当观测系统的方位角足够时,应只利用P波数据进行反演工作,以此提高反演结果的抗干扰性.

6 结论

微地震震源机制的准确反演对于水力压裂的成功实施具有重要的指导意义,通过Hudson图示法可以统计整个压裂过程中微地震震源机制的主要类型与性质,从而有效地评估整个压裂过程.由于微地震信号具有主频高,能量小的特征,即使采用井下观测系统,检波器接收到的微地震信号仍然会受到地层吸收衰减效应的严重影响.因此,本文重点分析了在均匀各向同性的吸收衰减地层条件下,微地震震源机制反演问题.考虑到实际地层Q值计算结果不稳定和误差较大的特性,本文利用数值算例分析不同地层Q值误差条件下不同观测数据对于微地震震源机制反演的分辨率变化.同时,数值计算的结果表明不同类型的微地震数据在反演过程中对于地层Q值的误差具有不同的敏感度.在观测系统的方位角足够的情况下,震源机制反演应优先考虑单独使用P波数据,或者采用直达P波和S波数据联合反演震源的T值,单独利用直达P波反演震源的k值,以此提高反演结果的稳定性和准确性.

通过理论分析和数值计算,本文也为提高微地震监测技术提供了理论基础,所获得的结论也可用于实际微地震监测的评价.然而,在实际的水力压裂过程中,目标层的物性参数产生了巨大的改变,这不仅会改变微地震震源机制的特征和微地震波场的传播规律,也会对整个反演过程产生严重的影响.这也是微地震监测技术当今所要面对的重大挑战,更深入的理论研究与详尽的岩石物理实验测试仍需开展.

References

Aki K,Richards P G.2002.Quantitative Seismology.2nd ed.US:University Science Books.

Backus G,Mulcahy M.1976a.Moment tensors and other phenomenological descriptions of seismic sources—I.Continuous displacements.Geophysical Journal International,46(2):341-361.Backus G,Mulcahy M.1976b.Moment tensors and other phenomenological descriptions of seismic sources—II.Discontinuous displacements.Geophysical Journal International,47(2):301-329.Baig A,Urbancic T.2010.Microseismic moment tensors:A path to understanding frac growth.The Leading Edge,29(3):320-324.

Brajanovski M,Müller T M,Gurevich B.2006.Characteristic frequencies of seismic attenuation due to wave-induced fluid flow in fractured porous media.Geophysical Journal International,166(2):574-578,doi:10.1111/j.1365-246X.2006.03068.x.

Chang X,Lu M X,Liu Y K.1999.Error analysis and appraisals for three general solutions in seismic tomography.Chinese Journal of Geophysics (in Chinese),42(5):695-701.

Chapman C H,Leaney W S.2012.A new moment-tensor decomposition for seismic events in anisotropic media.Geophysical Journal International,188(1):343-370,doi:10.1111/j.1365-246X.2011.05265.x.

Dando B D E,Chambers K,Velasco R.2014.A robust method for determining moment tensors from surface microseismic data.//SEG Technical Program Expanded Abstracts.SEG,2261-2266,doi:10.1190/segam2014-1481.1.

Das I,Zoback M D.2013.Long-period long-duration seismic events during hydraulic stimulation of shale and tight-gas reservoirs-Part 2:Location and mechanisms.Geophysics,78(6):KS109-KS117,doi:10.1190/Geo2013-0165.1.

Du J,Warpinski N R.2011.Uncertainty in FPSs from moment-tensor inversion.Geophysics,76(6):WC65-WC75,doi:10.1190/Geo2011-0024.1.

Eaton D W.2009.Resolution of microseismic moment tensors:A synthetic modeling study.//2009 SEG Annual Meeting,Society of Exploration Geophysicists.SEG,3569-3573.Eaton D W,Forouhideh F.2010.Microseismic moment tensors:The good,the bad and the ugly.CSEG Recorder,35(9):45-49.

Eaton D W,Forouhideh F.2011.Solid angles and the impact of receiver-array geometry on microseismic moment-tensor inversion.Geophysics,76(6):WC77-WC85,doi:10.1190/Geo2011-0077.1.Eaton D W,van der Baan M,Birkelo B,et al.2014.Scaling relations and spectral characteristics of tensile microseisms:evidence for opening/closing cracks during hydraulic fracturing.Geophysical Journal International,196(3):1844-1857,doi:10.1093/gji/ggt498.

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

Eisner L,Thornton M,Griffin J.2011.Challenges for microseismic monitoring.//2011 SEG Annual Meeting,Society of Exploration Geophysicists.SEG,1519-1523.Gharti H N,Oye V,Kühn D,et al.2011.Simultaneous microearthquake location and moment-tensor estimation using time-reversal imaging.//SEG Technical Program Expanded Abstracts 2011.SEG,1632-1637,doi:10.1190/1.3627516.

Grechka V.2010.Data-acquisition design for microseismic monitoring.The Leading Edge,29(3):278-282,doi:10.1190/1.3353723.

Grechka V,Blias E.2013.Analytic solutions to the joint estimation of microseismic event locations and effective velocity model.//SEG Technical Program Expanded Abstracts 2013.Houston,Texas:SEG,1972-1976,doi:10.1190/segam2013-0193.1.

Grechka V.2015.On the feasibility of inversion of single-well microseismic data for full moment tensor.Geophysics,80(4):KS41-KS49,doi:10.1190/geo2014-0471.1.

Hudson J A,Pearce R G,Rogers R M.1989.Source type plot for inversion of the moment tensor.Journal of Geophysical Research-Solid Earth,94(B1):765-774,doi:10.1029/JB094iB01p00765.

Julian B R,Miller A D,Foulger G R.1998.Non-double-couple earthquakes 1.Theory.Reviews of Geophysics,36(4):525-549,doi:10.1029/98rg00716.

Kendall M,Maxwell S,Foulger G,et al.2011.Microseismicity:Beyond dots in a box-Introduction.Geophysics,76(6):WC1-WC3,doi:10.1190/Geo-2011-1114-Spsein.1.

Leaney S,Chapman C,Ulrych T.2011.Microseismic source inversion in anisotropic media.CSEG Expanded Abstracts.

Maxwell S C.2011.What does microseismic tell us about hydraulic fracture deformation.CSEG Recorder,36(8):31-45.

Müller T M,Gurevich B,Lebedev M.2010.Seismic wave attenuation and dispersion resulting from wave-induced flow in porous rocks-A review.Geophysics,75(5):75A147-75A164,doi:10.1190/1.3463417.

Pride S R,Berryman J G,Harris J M.2004.Seismic attenuation due to wave-induced flow.Journal of Geophysical Research:Solid Earth (1978—2012),109(B1):B01201.

Quintal B,Jänicke R,Rubino J G,et al.2014.Sensitivity of S-wave attenuation to the connectivity of fractures in fluid-saturated rocks.Geophysics,79(5):WB15-WB24,doi:10.1190/Geo2013-0409.1.

Song F X,Toksöz M N.2011.Full-waveform based complete moment tensor inversion and source parameter estimation from downhole microseismic data for hydrofracture monitoring.Geophysics,76(6):WC103-WC116,doi:10.1190/Geo2011-0027.1.

Song F X,Warpinski N R,Toksöz M N.2014.Full-waveform based microseismic source mechanism studies in the Barnett Shale:Linking microseismicity to reservoir geomechanics.Geophysics,79(2):KS13-KS30,doi:10.1190/Geo2013-0094.1.Tonn R.1991.The determination of the seismic quality factor Q from Vspdata:a comparison of different computational methods.Geophysical Prospecting,39(1):1-27.Udías A,Madariaga R,Buforn E.2014.Source Mechanisms of Earthquakes:Theory and Practice.Cambridge:Cambridge University Press.

Urbancic T,Baig A.2013.How microseismic data from over a decade of cyclic steaming stimulations can yield insight into a heavy oil reservoir.//75th EAGE Conference &Exhibition-Workshops.EAGE.

Yang X C,Zhu H B,Cui S G,et al.2015.Application of P-wave first-motion focal mechanism solutions in microseismic monitoring for hydraulic fracturing.Geophysical Prospecting for Petroleum (in Chinese),54(1):43-50.Zhang C J,Ulrych T J.2002.Estimation of quality factors from CMP records.Geophysics,67(5):1542-1547,doi:10.1190/1.1512799.

附中文参考文献

常旭,卢孟夏,刘伊克.1999.地震层析成像反演中3种广义解的误差分析与评价.地球物理学报,42(5):695-701.

杨心超,朱海波,崔树果等.2015.P波初动震源机制解在水力压裂微地震监测中的应用.石油物探,54(1):43-50.

(本文编辑何燕)

Inversion for microseismic focal mechanisms in attenuated strata and its resolution

ZHAI Hong-Yu1,2,CHANG Xu1*,WANG Yi-Bo1,YAO Zhen-Xing1

1 Key Laboratory of Shale Gas and Geoengineering,Institute of Geology and Geophysics, Chinese Academy of Sciences,Beijing 100029,China 2 University of Chinese Academy of Sciences,Beijing 100049,China

Microseismic focal mechanism inversion is significant for the unconventional oil and gas development.Microseismic signal is characterized by high frequency and small energy,so it is easy to be influenced by formation attenuation effect,and to cause serious waveform distortion and unreliable inversion results.After analyzing the previous microseismic moment tensor inversion methods and some of their deficiency,this paper proposes a new method for microseismic focal mechanism inversion considering the stratum attenuation effect.On the other hand,the Q factor of actual strata is hard to estimate precisely and the microseismic events are very complex.So we design a simple observation model with three vertical wells and four basic moment tensors which are added 30% random perturbation to simulate the actual field situation.In the process of numerical computation,we assume the sediments′ Q factor has different levels of error.Based on this assumption,we can calculate the inversion error of the new moment tensor inversion method and analyze the inversion resolution by using the SVD of Frechet derivative matrix due to Q factor error.According to the theoretical calculation,the application conditions of the various focal mechanism inversions can be summarized for different types of microseismic data.Furthermore,the theoretical calculation proves that only using direct P waves can obtain a better inversion result than using direct P and S waves simultaneously.After analyzing the inversion error and the inversion resolution,we conclude that the joint inversion by using direct P and S wave data can estimate a better T value and a better k value can be inversed by using direct P wave separately.Moreover the direct P wave has a less sensitivity about the sediment Q factor.Based on this conclusion we can effectively reduce the negative influence on the focal mechanism inversion result by utilizing different types of wave data to determine the different parts of microseismic focal mechanisms.

Stratum attenuation;Microseismicity;Focal mechanism inversion;Inversion resolution

翟鸿宇,常旭,王一博等.2016.含衰减地层微地震震源机制反演及其反演分辨率.地球物理学报,59(8):3025-3036,

10.6038/cjg20160825.

Zhai H Y,Chang X,Wang Y B,et al.2016.Inversion for microseismic focal mechanisms in attenuated strata and its resolution.Chinese J.Geophys.(in Chinese),59(8):3025-3036,doi:10.6038/cjg20160825.

国家自然科学基金(41230317)以及中国科学院B类战略先导专项(XDB10030500)资助.

翟鸿宇,男,在读博士研究生,主要从事微地震震源机制反演与地层衰减机制研究.E-mail:zhaihongyu@hotmail.com

常旭,女,研究员.E-mail:changxu@mail.iggcas.ac.cn

10.6038/cjg20160825

P631

2015-08-01,2016-06-24收修定稿

猜你喜欢

张量震源反演
反演对称变换在解决平面几何问题中的应用
基于ADS-B的风场反演与异常值影响研究
偶数阶张量core逆的性质和应用
四元数张量方程A*NX=B 的通解
一类麦比乌斯反演问题及其应用
Pusher端震源管理系统在超高效混叠采集模式下的应用*
一类结构张量方程解集的非空紧性
震源的高返利起步
拉普拉斯变换反演方法探讨
M-张量的更多性质