一种微地震事件检测方法及系统技术方案

本发明专利技术公开一种微地震事件检测方法及系统。该方法包括：获取按时间顺序采集的待检测的微地震数据序列；建立第一时间窗和紧随第一时间窗的第二时间窗，第一时间窗和第二时间窗的容量均为奇数；对微地震数据序列的左侧和右侧进行数据填充；从第一时间窗的中心对准微地震数据序列的第一个数据滑动第一时间窗和第二时间窗，直至第一时间窗的中心对准微地震数据序列的最后一个数据，每滑动一个数据的步长，对第一时间窗内的数据与第二时间窗内的数据进行卡方检验，从而确定微地震数据序列中的各个数据处是否发生微地震波动。本发明专利技术的方法及系统，能够有效避免噪声对检测精度的影响，提高检测准确度。

【技术实现步骤摘要】
一种微地震事件检测方法及系统
本专利技术涉及微地震事件
，特别是涉及一种微地震事件检测方法及系统。
技术介绍
水力压裂微震监测技术是近年来在低渗透率储层压裂、油藏驱动和水驱前缘等领域发展起来的一项重要新技术，也是页岩气开发的重要支撑技术。该项技术在邻井中布置多级三分量检波器排列，监测压裂井目的层段在水力压裂过程中所产生的微震事件，反演微震事件求取震源位置等参数，从而描述水力压裂过程中裂缝生长的几何形状及空间分布，实时提供水力压裂产生裂缝的长度、高度、宽度及方位，实现页岩气的工业化开发。水力压裂微震检测是当前页岩气开发领域科学研究的热点和难点。微震监测系统中重要的一项工作是微震事件的定位。定位精度是影响微震监测系统应用效果的最为重要的因素，而微震事件定位的准确程度则主要依赖于波动初至(又可称为初至)读取的准确性等有关因素。但问题是，初至拾取并不如想象中的那般简单。受地面仪器采动以及地质构造的影响，岩石破裂形式十分复杂，继而产生各种形式和能量的微震波动，其形式可多达几十甚至上百种，不仅主频、延时和能量等方面有差异，而且在初至位置附近的波形形态差异巨大，这种波形特征的不统一为初至拾取带来了很大困难。除了初至点特征复杂外，初至拾取还面临着另外一个更大的挑战：微震记录是海量数据。并且，这些数据中有很大一部分都是人类或者机械活动所造成的噪声和干扰，与微震无关。除此之外，微震信号本身也并不纯粹，例如我国学者窦林名教授等认为微震信号包括多种信号。目前生产上多采取人工方法分析微震记录从而判断微震事件，然而人工方法无法准确避免噪声和干扰，精度与可靠性差。
技术实现思路
本专利技术的目的是提供一种微地震事件检测方法及系统，有效避免噪声对检测精度的影响，提高检测准确度。为实现上述目的，本专利技术提供了如下方案：一种微地震事件检测方法，包括：获取按时间顺序采集的待检测的微地震数据序列；建立第一时间窗和第二时间窗，所述第一时间窗和所述第二时间窗的容量均为L个数据，其中L为奇数；所述第二时间窗紧随所述第一时间窗；对所述微地震数据序列的左侧和右侧进行数据填充，使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述微地震数据序列的第一个数据滑动至所述第一时间窗的中心对准所述微地震数据序列的最后一个数据时，所述第一时间窗和所述第二时间窗内均填充满L个数据；从所述第一时间窗的中心对准所述微地震数据序列的第一个数据滑动所述第一时间窗和所述第二时间窗，直至所述第一时间窗的中心对准所述微地震数据序列的最后一个数据，每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述微地震数据序列中的各个数据处是否发生微地震波动；记录所有发生微地震波动的数据在所述微地震数据序列中的序号。可选的，所述对所述微地震数据序列的左侧和右侧进行数据填充，具体包括：将所述微地震数据序列的左侧补充至少个数据；左侧补充的至少个数据均与所述微地震数据序列最左侧的数据相同；将所述微地震数据序列的右侧补充至少个数据；右侧补充的至少个数据均与所述微地震数据序列最右侧的数据相同。可选的，所述每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述微地震数据序列中的各个数据处是否发生微地震波动，具体包括：求取所述第一时间窗内所有数据的最大值和最小值，得到取值区间；将所述取值区间划分为K段，得到K段小区间；计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n)P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n)Q(k)；其中n表示位于所述第一时间窗中心的数据在所述微地震数据序列中的序号；k为小区间序号，k＝1,2,…,K；计算数量差值比率：根据参数K和置信度α，通过查询卡方表获取卡方值χ2(α,K-1)；判断所述数量差值比率是否大于或等于所述卡方值，得到第一判断结果；若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值，则确定在所述微地震数据序列的第n个数据处发生微地震波动；若所述第一判断结果表示所述数量差值比率小于所述卡方值，则确定在所述微地震数据序列的第n个数据处未发生微地震波动。本专利技术还公开一种微地震事件检测系统，包括：序列获取模块，用于获取按时间顺序采集的微地震数据序列；时间窗建立模块，用于建立第一时间窗和第二时间窗，所述第一时间窗和所述第二时间窗的容量均为L个数据，其中L为奇数；所述第二时间窗紧随所述第一时间窗；数据填充模块，用于对所述微地震数据序列的左侧和右侧进行数据填充，使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述微地震数据序列的第一个数据滑动至所述第一时间窗的中心对准所述微地震数据序列的最后一个数据时，所述第一时间窗和所述第二时间窗内均填充满L个数据；微震事件检测模块，用于从所述第一时间窗的中心对准所述微地震数据序列的第一个数据滑动所述第一时间窗和所述第二时间窗，直至所述第一时间窗的中心对准所述微地震数据序列的最后一个数据，每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述微地震数据序列中的各个数据处是否发生微地震波动；记录模块，用于记录所有发生微地震波动的数据在所述微地震数据序列中的序号。可选的，所述数据填充模块，具体包括：序列左填充单元，用于将所述微地震数据序列的左侧补充至少个数据；左侧补充的至少个数据均与所述微地震数据序列最左侧的数据相同；序列右填充单元，用于将所述微地震数据序列的右侧补充至少个数据；右侧补充的至少个数据均与所述微地震数据序列最右侧的数据相同。可选的，所述微震事件检测模块，具体包括：区间计算单元，用于求取所述第一时间窗内所有数据的最大值和最小值，得到取值区间；区间划分单元，用于将所述取值区间划分为K段，得到K段小区间；区间数据量计算单元，用于计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n)P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n)Q(k)；其中n表示位于所述第一时间窗中心的数据在所述微地震数据序列中的序号；k为小区间序号，k＝1,2,…,K；比率计算单元，用于计算数量差值比率：卡方值查询单元，用于根据参数K和置信度α，通过查询卡方表获取卡方值χ2(α,K-1)；判断单元，用于判断所述数量差值比率是否大于或等于所述卡方值，得到第一判断结果；微地震波动确定单元，用于若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值，则确定在所述微地震数据序列的第n个数据处发生微地震波动；若所述第一判断结果表示所述数量差值比率小于所述卡方值，则确定在所述微地震数据序列的第n个数据处未发生微地震波动。根据本专利技术提供的具体实施例，本专利技术公开了以下技术效果：本专利技术的微地震事件检测方法及系统，利用卡方检验的方法识别微地震波动，关注微地震数据落入某个区间的数量，能够避免依赖于微震数据本身的大小，从而有效避免噪声的影响，提高检测准确度。附图说明为了更清楚地说明本专利技术实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本专利技术的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动性的前提下，还可以根据这本文档来自技高网...

【技术保护点】
1.一种微地震事件检测方法，其特征在于，包括：获取按时间顺序采集的待检测的微地震数据序列；建立第一时间窗和第二时间窗，所述第一时间窗和所述第二时间窗的容量均为L个数据，其中L为奇数；所述第二时间窗紧随所述第一时间窗；对所述微地震数据序列的左侧和右侧进行数据填充，使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述微地震数据序列的第一个数据滑动至所述第一时间窗的中心对准所述微地震数据序列的最后一个数据时，所述第一时间窗和所述第二时间窗内均填充满L个数据；从所述第一时间窗的中心对准所述微地震数据序列的第一个数据滑动所述第一时间窗和所述第二时间窗，直至所述第一时间窗的中心对准所述微地震数据序列的最后一个数据，每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述微地震数据序列中的各个数据处是否发生微地震波动；记录所有发生微地震波动的数据在所述微地震数据序列中的序号。

【技术特征摘要】
1.一种微地震事件检测方法，其特征在于，包括：获取按时间顺序采集的待检测的微地震数据序列；建立第一时间窗和第二时间窗，所述第一时间窗和所述第二时间窗的容量均为L个数据，其中L为奇数；所述第二时间窗紧随所述第一时间窗；对所述微地震数据序列的左侧和右侧进行数据填充，使所述第一时间窗和所述第二时间窗从所述第一时间窗的中心对准所述微地震数据序列的第一个数据滑动至所述第一时间窗的中心对准所述微地震数据序列的最后一个数据时，所述第一时间窗和所述第二时间窗内均填充满L个数据；从所述第一时间窗的中心对准所述微地震数据序列的第一个数据滑动所述第一时间窗和所述第二时间窗，直至所述第一时间窗的中心对准所述微地震数据序列的最后一个数据，每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述微地震数据序列中的各个数据处是否发生微地震波动；记录所有发生微地震波动的数据在所述微地震数据序列中的序号。2.根据权利要求1所述的一种微地震事件检测方法，其特征在于，所述对所述微地震数据序列的左侧和右侧进行数据填充，具体包括：将所述微地震数据序列的左侧补充至少个数据；左侧补充的至少个数据均与所述微地震数据序列最左侧的数据相同；将所述微地震数据序列的右侧补充至少个数据；右侧补充的至少个数据均与所述微地震数据序列最右侧的数据相同。3.根据权利要求1所述的一种微地震事件检测方法，其特征在于，所述每滑动一个数据的步长，对所述第一时间窗内的数据与所述第二时间窗内的数据进行卡方检验，从而确定所述微地震数据序列中的各个数据处是否发生微地震波动，具体包括：求取所述第一时间窗内所有数据的最大值和最小值，得到取值区间；将所述取值区间划分为K段，得到K段小区间；计算所述第一时间窗内的数据介于每个所述小区间内的数量F(n)P(k)和所述第二时间窗内的数据介于每个所述小区间内的数量F(n)Q(k)；其中n表示位于所述第一时间窗中心的数据在所述微地震数据序列中的序号；k为小区间序号，k＝1,2,…,K；计算数量差值比率：根据参数K和置信度α，通过查询卡方表获取卡方值χ2(α,K-1)；判断所述数量差值比率是否大于或等于所述卡方值，得到第一判断结果；若所述第一判断结果表示所述数量差值比率大于或等于所述卡方值，则确定在所述微地震数据序列的第n个数据处发生微地震波动；若所述第一判断结果表示所述数量差值比率小于所述卡方值，则确定在所述微地震数据序列的第n个数据处未发生微地震波动。4.一种微地震事件检测系统，其特征在于，包括：序列获取...

【专利技术属性】
技术研发人员：翟明岳，
申请(专利权)人：广东石油化工学院，
类型：发明
国别省市：广东,44