CN110531412B - 一种计算井中微地震事件相对方位角的方法 - Google Patents

一种计算井中微地震事件相对方位角的方法 Download PDF

Info

Publication number
CN110531412B
CN110531412B CN201910922076.2A CN201910922076A CN110531412B CN 110531412 B CN110531412 B CN 110531412B CN 201910922076 A CN201910922076 A CN 201910922076A CN 110531412 B CN110531412 B CN 110531412B
Authority
CN
China
Prior art keywords
event
template
angle
target
target event
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201910922076.2A
Other languages
English (en)
Other versions
CN110531412A (zh
Inventor
陈海潮
程冬冬
孟晓波
钮凤林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China University of Petroleum Beijing CUPB
Original Assignee
China University of Petroleum Beijing CUPB
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China University of Petroleum Beijing CUPB filed Critical China University of Petroleum Beijing CUPB
Priority to CN201910922076.2A priority Critical patent/CN110531412B/zh
Publication of CN110531412A publication Critical patent/CN110531412A/zh
Application granted granted Critical
Publication of CN110531412B publication Critical patent/CN110531412B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. analysis, for interpretation, for correction
    • G01V1/288Event detection in seismic signals, e.g. microseismics

Abstract

本发明公开了一种计算井中微地震事件相对方位角的方法,其特征在于,包括以下步骤:1)确定模板事件和目标事件;在观测站数据中取任意两组微地震事件,以其中一微地震事件为模板事件,另一个微地震事件为目标事件;2)沿目标事件的反方位角方向以预定角度旋转目标事件,并计算旋转后的目标事件在模板事件上的投影值;3)确定相对方向角;确定旋转后的目标事件在模板事件的最大投影值,最大投影值对应的角度α即为模板事件和目标事件的相对方位角,角度α={‑π+目标事件的旋转次数×预定角度}。

Description

一种计算井中微地震事件相对方位角的方法
技术领域
本发明涉及一种计算井中微地震事件相对方位角的方法,属于井中微地震三分量信号处理领域。
背景技术
井中微地震监测技术是水力压裂裂缝评价的重要新技术之一。该项技术通过布置在邻井中的三分量检波器,记录伴随裂缝剪切滑移产生的微地震信号,对水力压裂裂缝的空间展布和储层改造体积进行定量评估。在单井监测条件下,常规处理方法是采用极化分析和初至拾取获取相对震源的方位角和纵横波时差,在此基础上进一步反演获得微地震的震源位置。因此,准确计算方位角是实现微地震事件精确定位的一项关键技术。
目前国内外最常用的偏振分析方法是矢端曲线法,其原理是:利用微地震事件在水平分量上的P波信号构造协方差矩阵,计算协方差矩阵的特征值,并得到相应的特征向量,最大特征值对应的特征向量即是P波的极化方向。该方法假定直达P波偏振方向和波传播方向一致。然而,在低信噪比条件下,P波信号的质点轨迹不一定呈现严格线性,矢端曲线法获得的偏振方位角存在较大不确定性和误差。
发明内容
针对上述问题,本发明的目的是提供一种计算井中微地震事件相对方位角的方法,该方法能够用于复杂信噪比条件下井中微地震三分量数据资料处理,能够显著提高低信噪比条件下的微地震事件的定位精度。
为实现上述目的,本发明采用以下技术方案,一种计算井中微地震事件相对方位角的方法,其特征在于,包括以下步骤:
1)确定模板事件和目标事件;
在观测站数据中取两组微地震事件,以其中一微地震事件为模板事件,另一个微地震事件为目标事件;
2)沿目标事件的反方位角方向以预定角度旋转目标事件,并计算旋转后的目标事件在模板事件上的投影值;
3)确定相对方向角;
确定旋转后的目标事件在模板事件的最大投影值,最大投影值对应的角度α即为模板事件和目标事件的相对方位角,角度α={-π+目标事件的旋转次数×预定角度}。
进一步地,在上述步骤2)中,旋转后的目标事件在模板事件上的投影值的计算公式如下:
式中,C(α)是目标事件在角度为α时在模板事件上的投影值;是目标事件向量;为模板事件向量的倒置;θ是模板事件的方位角;θ0是目标事件的方位角;α为-π与旋转角度之和,旋转角度为目标事件的旋转次数×预定角度;α的范围为-π到π。
进一步地,旋转后的目标事件在模板事件上的投影值的计算公式推导过程如下:
事件投影值的计算公式如下:
式中,是模板事件向量;为目标事件旋转后的向量;ET是模板事件在东方向分量的转置;NT是模板事件在北方向分量的转置;是旋转后的目标事件在东方向的分量;是旋转后的目标事件在北方向的分量;
参考事件向量公式表示如下:
式中,E0是目标事件在东方向的分量;N0是目标事件在北方向的分量;e0是目标事件中东方向的噪音分量;n0是目标事件中北方向的噪音分量;
其中,
旋转后的目标事件向量公式表示如下:
其中,
目标事件的转换公式如下:
式中,Γ(α)是坐标旋转变换公式;如下:
目标事件的向量公式如下:
其中,
式(9)、(10)中,E为模板事件在东方向的分量,N为模板事件在北方向的分量;e为模板事件中东方向的噪声分量;n为模板事件中北方向的噪声分量;
将式(2)至(10)结合起来,得到
在观测站实际采集的数据中,模板事件与目标事件中均包含噪音信号,假设微地震事件中的噪音是白高斯噪音,并且目标事件与模板事件中的噪音不存在相关性;则式(11)中的第二项和第三项均为常数零,得到事件投影值的计算公式为:
进一步地,还包括步骤4)计算目标事件的方位角;采用矢端曲线法计算出模板事件的方位角,目标事件的方位角为(模板事件的方位角)+(相对方位角α)。
进一步地,在上述步骤1)中,观测站是指设置在井中的多个三分量检波器,观测站数据是指从设置在井中的多个三分量检波器中获取的数据。
进一步地,在上述步骤2)中,预定角度为固定角度,固定角度值取0.01°。
进一步地,在上述步骤2)中,方位角是指微地震事件的发生位置与一个三分量检波器的设置位置的连线在水平面的投影线与正北方向的夹角,并且所述夹角沿着顺时针的旋转方向进行计算;反方位角与方位角的方向相反。
本发明采用以上技术方案,其具有如下优点:1、本发明通过反方位角旋转目标事件,并计算旋转后的目标事件在模板事件上的投影值,获得的最大投影值对应的角度即为模板事件和目标事件的相对方位角,本发明的方法能够准确计算低信噪比条件下的微地震事件相对方位角,能够用于复杂信噪比条件下井中微地震三分量数据资料处理,能够显著提高低信噪比条件下的微地震事件的定位精度。
2、本发明能够确定目标事件和模板事件的相对方位角,仅用矢端曲线法计算模板事件的方位角,进而可以确定目标事件的方位角,本发明的方法能够较准确计算低信噪比条件下的微地震事件方位角,能够用于复杂信噪比条件下井中微地震三分量数据资料处理,能够显著提高低信噪比条件下的微地震事件的定位精度。
附图说明
图1是利用本发明的方法求取相对方位角的示意图;
图2是利用矢端曲线法(Particle Motion,简写PM)求取方位角的示意图;
图3是在高信噪比(S/N=10)下应用本发明的方法与矢端曲线法求取微地震事件方位角的测试结果;其中,图a是展示高信噪比地震记录中地震信号的东方向和北方向的分量的原始波形记录;图b是投影值和方位角之间的抛物线状的点线图;图c左是平均信噪比(MeanS/N=10.5)下利用本发明的方法计算方向角的误差统计分布柱状图;图c右是平均信噪比(MeanS/N=10.5)下利用矢端曲线法计算方向角的误差统计分布柱状图;
图4是在低信噪比(S/N=1)下应用本发明的方法与矢端曲线法求取微地震事件方位角的测试结果;其中,图d是展示地震记录中低信噪比地震记录中地震信号的东方向和北方向的分量的原始波形记录;图e是投影值和方位角之间的抛物线状的点线图;图f左是平均信噪比(MeanS/N=1.5)下利用本发明的方法计算方向角的误差统计分布柱状图;图f右是平均信噪比(MeanS/N=1.5)下利用矢端曲线法计算方向角的误差统计分布柱状图。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。然而应当理解,附图的提供仅为了更好地理解本发明,它们不应该理解成对本发明的限制。
如图1所示,本发明提供了一种计算井中微地震事件相对方位角的方法,又称为最大投影值法,(Maximum Projection,简写MP),其包括以下步骤:
1)确定模板事件和目标事件;
在观测站数据中取两组微地震事件,以其中一微地震事件为模板事件,另一个微地震事件为目标事件;
观测站是指设置在井中的多个三分量检波器,观测站数据是指从设置在井中的多个三分量检波器中获取的数据。
2)沿参考事件的反方位角方向(Backazimuth)以预定角度旋转目标事件,并计算旋转后的目标事件在模板事件上的投影值;
其中,方位角(azimuth)是指微地震事件的发生位置与一个三分量检波器的设置位置的连线在水平面的投影线与正北方向的夹角,并且所述夹角沿着顺时针的旋转方向进行计算;反方位角与方位角的方向相反;
预定角度为固定角度,选取的固定角度越小,精确度越高,得到的结果越精确,但在实际中需要考虑到计算机的运行成本问题,因此需要取适当大小值,一般取0.01°。
3)确定相对方向角;
确定旋转后的目标事件在模板事件的最大投影值,最大投影值对应的角度α即为模板事件和目标事件的相对方位角,角度α={-π+目标事件的旋转次数×预定角度}。
进一步地,在步骤2)中,旋转后的目标事件在模板事件上的投影值的计算公式如下:
式中,C(α)是目标事件旋转后在模板事件上的投影值;是目标事件向量;为模板事件向量的倒置;θ是模板事件的方位角;θ0是目标事件的方位角;α为-π与旋转角度之和,旋转角度为目标事件的旋转次数×预定角度;α范围为-π到π。
进一步地,旋转后的目标事件在模板事件上的投影值的计算公式的推导过程如下:
事件投影值的计算公式如下:
式中,是模板事件向量;为目标事件旋转后的向量;ET是模板事件在东方向分量的转置;NT是模板事件在北方向分量的转置;是旋转后的目标事件在东方向的分量;是旋转后的目标事件在北方向的分量;
目标事件向量公式表示如下:
式中,E0是目标事件在东方向的分量;N0是目标事件在北方向的分量;e0是目标事件中东方向的噪音分量;n0是目标事件中北方向的噪音分量;
其中,
旋转后的目标事件向量公式表示如下:
其中,
目标事件的转换公式如下:
式中,Γ(α)是坐标旋转变换公式;如下:
模板事件的向量公式如下:
其中,
式(9)、(10)中,E为模板事件在东方向的分量,N为模板事件在北方向的分量;e为模板事件中东方向的噪声分量;n为模板事件中北方向的噪声分量;
将式(2)至(10)结合起来,得到
在观测站实际采集的数据中,模板事件与目标事件中均包含噪音信号,假设微地震事件中的噪音是白高斯噪音,并且目标事件与模板事件中的噪音不存在相关性;则式(11)中的第二项和第三项均为常数零,故得到事件投影值的计算公式为:
进一步地,本发明还包括步骤4)计算目标事件的方位角;采用矢端曲线法计算出模板事件的方位角,目标事件的方位角即为(模板事件的方位角)+(相对方位角α);这样,仅在模板事件的方位角计算中使用矢端曲线法,在其他任意微地震事件的方位角计算时,利用本发明的最大投影值法计算其与模板事件的相对方位角,进而可以获得任意微地震事件的方位角,避免采用纯矢端曲线法计算方位角的较大不确定性和误差。
利用本发明的方法和矢端曲线法分别计算高信噪比(S/N=10)和低信噪比(S/N=1)下计算微地震事件的方位角,测试结果如图3、图4所示。
结合图b、e,其中,θ`是指微地震事件的方位角;θMP为采用本发明的方法获得微地震事件的方向角,θPM为矢端曲线法获得的微地震事件的方向角,可以看出,与矢端曲线法相比,利用本发明的方法计算在高信噪比(S/N=10)或低信噪比(S/N=1)下的方位角更接近实际数值,计算更精确;结合图c、f可以看出,与矢端曲线法相比,无论对平均信噪比为10.5或平均信噪比为1.5的数据,本发明的方法计算方位角产生的误差都更小。
从上述对比可以看出,本发明的方法能够较准确计算低信噪比条件下的微地震事件方位角,能够用于复杂信噪比条件下井中微地震三分量数据资料处理,能够显著提高低信噪比条件下的微地震事件的定位精度。
本发明仅以上述实施例进行说明,各部件的结构、设置位置及其连接都是可以有所变化的。在本发明技术方案的基础上,凡根据本发明原理对个别部件进行的改进或等同变换,均不应排除在本发明的保护范围之外。

Claims (6)

1.一种计算井中微地震事件相对方位角的方法,其特征在于,包括以下步骤:
1)确定模板事件和目标事件;
在观测站数据中取两组微地震事件,以其中一微地震事件为模板事件,另一个微地震事件为目标事件;
2)沿目标事件的反方位角方向以预定角度旋转目标事件,并计算旋转后的目标事件在模板事件上的投影值;
3)确定相对方向角;
确定旋转后的目标事件在模板事件的最大投影值,最大投影值对应的角度α即为模板事件和目标事件的相对方位角,角度α={-π+目标事件的旋转次数×预定角度};
在上述步骤2)中,旋转后的目标事件在模板事件上的投影值的计算公式如下:
式中,C(α)是目标事件旋转后在模板事件上的投影值;是目标事件向量;为模板事件向量的倒置;θ是模板事件的方位角;θ0是目标事件的方位角;α为-π与旋转角度之和,旋转角度为目标事件的旋转次数×预定角度;α范围为-π到π。
2.如权利要求1所述的一种计算井中微地震事件相对方位角的方法,其特征在于,旋转后的目标事件在模板事件上的投影值的计算公式推导过程如下:
事件投影值的计算公式如下:
式中,是模板事件向量;为目标事件旋转后的向量;ET是模板事件在东方向分量的转置;NT是模板事件在北方向分量的转置;是旋转后的目标事件在东方向的分量;是旋转后的目标事件在北方向的分量;
目标事件向量公式表示如下:
式中,E0是目标事件在东方向的分量;N0是目标事件在北方向的分量;e0是目标事件中东方向的噪音分量;n0是目标事件中北方向的噪音分量;
其中,
旋转后的目标事件向量公式表示如下:
其中,
目标事件的转换公式如下:
式中,Γ(α)是坐标旋转变换公式;如下:
模板事件向量公式如下:
其中,
式(9)、(10)中,E为模板事件在东方向的分量,N为模板事件在北方向的分量;e为模板事件中东方向的噪声分量;n为模板事件中北方向的噪声分量;
将式(2)至(10)结合起来,得到
在观测站实际采集的数据中,模板事件与目标事件中均包含噪音信号,假设微地震事件中的噪音是白高斯噪音,并且目标事件与模板事件中的噪音不存在相关性;则式(11)中的第二项和第三项均为常数零,得到事件投影值的计算公式为:
3.如权利要求1所述的一种计算井中微地震事件相对方位角的方法,其特征在于:还包括步骤4)计算目标事件的方位角;采用矢端曲线法计算出模板事件的方位角,目标事件的方位角为(模板事件的方位角)+(相对方位角)。
4.如权利要求1所述的一种计算井中微地震事件相对方位角的方法,其特征在于:在上述步骤1)中,观测站是指设置在井中的多个三分量检波器,观测站数据是指从设置在井中的多个三分量检波器中获取的数据。
5.如权利要求1所述的一种计算井中微地震事件相对方位角的方法,其特征在于:在上述步骤2)中,预定角度为固定角度,固定角度值取0.01°。
6.如权利要求1所述的一种计算井中微地震事件相对方位角的方法,其特征在于:在上述步骤2)中,方位角是指微地震事件的发生位置与一个三分量检波器的设置位置的连线在水平面的投影线与正北方向的夹角,并且所述夹角沿着顺时针的旋转方向进行计算;反方位角与方位角的方向相反。
CN201910922076.2A 2019-09-27 2019-09-27 一种计算井中微地震事件相对方位角的方法 Active CN110531412B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910922076.2A CN110531412B (zh) 2019-09-27 2019-09-27 一种计算井中微地震事件相对方位角的方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910922076.2A CN110531412B (zh) 2019-09-27 2019-09-27 一种计算井中微地震事件相对方位角的方法

Publications (2)

Publication Number Publication Date
CN110531412A CN110531412A (zh) 2019-12-03
CN110531412B true CN110531412B (zh) 2020-09-29

Family

ID=68670487

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910922076.2A Active CN110531412B (zh) 2019-09-27 2019-09-27 一种计算井中微地震事件相对方位角的方法

Country Status (1)

Country Link
CN (1) CN110531412B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111257940B (zh) * 2020-02-14 2020-09-04 中国科学院地质与地球物理研究所 一种组合式海底地震仪姿态自动校正装置及方法

Family Cites Families (15)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2444954B (en) * 2006-12-20 2009-05-20 Westerngeco Seismic Holdings Method of monitoring microseismic events
US7647183B2 (en) * 2007-08-14 2010-01-12 Schlumberger Technology Corporation Method for monitoring seismic events
CN103105624B (zh) * 2011-11-11 2015-07-22 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于数据库技术的纵横波时差定位方法
CN102841373B (zh) * 2012-08-23 2015-02-04 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于方位角约束的微地震事件定位方法
CN104216005B (zh) * 2013-06-03 2017-05-03 中国石油化工股份有限公司 一种井中微地震三分量旋转方法
CN104216008B (zh) * 2013-06-05 2017-02-08 中国石油天然气集团公司 一种井中压裂微地震事件识别方法
CN104280772A (zh) * 2013-07-12 2015-01-14 中国石油化工股份有限公司 一种井中微地震震相识别方法
CN103389489B (zh) * 2013-07-31 2015-03-04 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于大斜度井的微地震监测定位方法
CN104459797B (zh) * 2013-09-22 2017-05-03 中国石油化工股份有限公司 一种井中微地震事件识别和拾取的方法
CN104182651B (zh) * 2014-09-12 2017-03-15 中国石油集团川庆钻探工程有限公司地球物理勘探公司 用于三分量检波器接收的微地震事件方位角自动质控方法
CN104502976B (zh) * 2014-12-29 2017-02-22 中国石油天然气集团公司 一种微地震监测中的井中检波器水平分量定向方法
CN106168675B (zh) * 2015-05-18 2018-05-08 中国石油化工股份有限公司 井中微地震p/s波事件识别方法和装置
CN107346349B (zh) * 2016-05-06 2020-12-01 中国石油化工股份有限公司 基于多孔多级的井中微地震方位计算方法和装置
CN107290722B (zh) * 2017-06-29 2019-11-26 中国石油大学(北京) 微震源的定位方法和装置
CN107607990B (zh) * 2017-08-07 2019-09-10 中国石油天然气集团公司 三分量检波器水平分量的方向检测方法和装置

Also Published As

Publication number Publication date
CN110531412A (zh) 2019-12-03

Similar Documents

Publication Publication Date Title
WO2017167191A1 (zh) 地震数据处理方法和装置
CN107770859A (zh) 一种考虑基站位置误差的tdoa‑aoa定位方法
CN110531412B (zh) 一种计算井中微地震事件相对方位角的方法
CN106154332A (zh) 一种井中微地震纵横波事件初至识别方法
WO2016115763A1 (zh) 共姿态道集方位角分析和校正方法及装置
CN101609160A (zh) 一种抗背景噪音强的地震数据初至波自动拾取方法
WO2009105964A1 (zh) 一种三维地震记录的叠前似二维变换的方法
CN107356947A (zh) 基于单频导航卫星数据确定卫星差分伪距偏差的方法
CN103675916A (zh) 一种三分量检波器埋置方向高精度校正的方法
CN103852785A (zh) 地层各向异性的评价方法
Xu et al. Systemic comparison of seismometer horizontal orientations based on teleseismic earthquakes and ambient‐noise data
CN109085642B (zh) 一种各向异性介质微地震事件定位方法
Zheng et al. Study on the calibration method of USBL system based on ray tracing
CN107843919B (zh) 一种微地震监测数据振幅均衡方法
CN104216005B (zh) 一种井中微地震三分量旋转方法
CN106522924A (zh) 一种随钻测量中方位角的获取方法
US5532978A (en) Method of correction of seismic data for the effects of azimuthal anisotropy
Shen et al. Site-specific real-time GPS multipath mitigation based on coordinate time series window matching
CN110398777B (zh) 一种井中检波器定位方法
CN107504862B (zh) 一种全方位高精度激光定位方法
CN110736963A (zh) 一种基于CSI的室内Wi-Fi定位方法、装置及存储介质
Li et al. Methods for Amplitude Calibration and Orientation Discrepancy Measurement: Comparing Co‐Located Sensors of Different Types in the Southern California Seismic Network
Bolarinwa et al. Calibrating the 2016 IRIS Wavefields Experiment nodal sensors for amplitude statics and orientation errors
Jarvis et al. Reduction of ambiguities in HF radar results through a revised receiving antenna array and sounding pattern
US9213100B1 (en) Bearing-only tracking for horizontal linear arrays with rapid, accurate initiation and a robust track accuracy threshold

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant