【技术实现步骤摘要】
一种利用累积相似性的微震事件检测方法和系统
[0001]本专利技术涉及石油领域,特别是涉及一种微震事件检测方法及系统。
技术介绍
[0002]水力压裂微震监测技术是近年来在低渗透率储层压裂、油藏驱动和水驱前缘等领域发展起来的一项重要新技术,也是页岩气开发的重要支撑技术。该项技术在邻井中布置多级三分量检波器排列,监测压裂井目的层段在水力压裂过程中所产生的微震事件,反演微震事件求取震源位置等参数,从而描述水力压裂过程中裂缝生长的几何形状及空间分布,实时提供水力压裂产生裂缝的长度、高度、宽度及方位,实现页岩气的工业化开发。水力压裂微震检测是当前页岩气开发领域科学研究的热点和难点。从社会和国家的需求角度考虑,开展微震监测系统方面的研究十分重要,具有重大的社会和经济价值。
[0003]微震监测系统中重要的一项工作是微震事件的定位。定位精度是影响微震监测系统应用效果的最为重要的因素,而微震事件定位的准确程度则主要依赖于波动初至(又可称为初至)读取的准确性等有关因素。但问题是,初至拾取并不如想象中的那般简单。受地面仪器采动以及地质构造的影响,岩石破裂形式十分复杂,继而产生各种形式和能量的微震波动,其形式可多达几十甚至上百种,不仅主频、延时和能量等方面有差异,而且在初至位置附近的波形形态差异巨大,这种波形特征的不统一为初至拾取到来了很大困难。进一步的研究还表明,微震震源机制也会影响初至点特征:硬岩剪切作用产生的微震波动大多能量大、主频较高、延时短以及最大峰值位置紧跟初始初至,这类波的初至点清晰、起跳延时短,拾取较为容易;但拉 ...
【技术保护点】
【技术特征摘要】
1.一种利用累积相似性的微震事件检测方法,其特征在于,包括:步骤101获取按时间顺序采集的信号序列S;步骤102对信号序列进行中值滤波处理,具体为:中值滤波后的信号序列记为S
mean
,其第n个元素记为求取方法为:其中:Mean[*]表示对序列*中的元素取中间值;*表示任一序列或任一元素;为所述信号序列S的第|n-K+1|
N
个元素;为所述信号序列S的第|n-K+2|
N
个元素;为所述信号序列S的第|n+K|
N
个元素;|*|
N
表示以N为模对*取余数;为中值滤波窗口长度;为所述信号序列S的方差;N为所述信号序列S的长度;n=1,2,
···
,N为元素序号;步骤103生成N个窗口信号序列,具体为:第o个窗口信号序列记为b
o
,其第l个元素为生成公式为:其中:为窗口元素序号;o=1,2,
···
,N为窗口序号;为所述中值滤波后的信号序列S
mean
的第|o+l|
N
个元素;步骤104求取N个候选微震事件起点和终点,具体为:第o个候选微震事件起点记为其求取方法为:所述第o个候选微震事件终点记为其求取方法为:其中:m
o
为所述第o个窗口信号序列b
o
的均值;σ
o
为所述第o个窗口信号序列b
o
的均方差;步骤105求取N个候选微震事件矢量,具体为:第o个候选微震事件矢量记为h
o
,其求取公式为:其中:
为所述第o个窗口信号序列b
o
的第元素;为所述第o个窗口信号序列b
o
的第元素;为所述第o个窗口信号序列b
o
的第元素;步骤106获取典型微震事件波形,具体为:通过实际调查,由专家研判所调查地区的典型微震事件波形,并记为c,其第z个元素记为c
z
;其中:z=1,2,
···
,N
c
为波形元素序号;N
c
为所述典型微震事件波形的长度,通过实际调查研判得到;步骤107求取N个相似矩阵,具体为:第o个相似矩阵记为D
o
,其第i行第j列元素为求取公式为其中:为所述第o个候选微震事件矢量h
o
的第j个元素;c
i
为所述所述典型微震事件波形c的第i个元素;i=1,2,
···
,N
c
为行序号;为列序号;步骤108求取N个拉伸比较路径和典型路径,具体为:第o个拉伸比较路径记为p
o
;其第1个元素为其第N
p
个元素为第o个典型路径记为r
o
,其第1个元素为其第N
p
个元素为所述第o个拉伸比较路径p
o
的第i
P
个元素记为所述第o个典型路径r
o
的第i
P
个元素记为求取公式为:s.t.1<i
P
≤N
p
其中:为路径长度;步骤109计算N个累计相似性值,具体...
还没有人留言评论。发表了对其他浏览者有用的留言会获得科技券。