一种微地震震相识别初至拾取方法、设备和存储介质
技术领域
本发明涉及震相识别技术领域,更具体地说,涉及一种基于到时曲线拟合的微地震震相识别初至拾取方法、一种计算机设备和一种存储有计算机可读指令的存储介质。
背景技术
非常规油气藏储层的渗透率极低,在开采过程中需要采用水平井和分段压裂进行储层改造。微地震监测技术可以实时监测水力压裂过程,利用监测数据定位诱发微地震事件,根据破裂事件点的分布特征描述压裂裂缝发育的几何形态和空间分布,对评价压裂改造效果并优化后期压裂设计和井位布局有重要意义。微地震监测数据处理过程中,微地震事件识别及初至拾取是微地震震源定位的基础,而且初至拾取精度直接影响的事件定位、震源机制分析以及裂缝成像效果。
微地震事件识别以有效信号和环境噪声的差异识别为基础,常用于信号识别的属性特征包括能量、偏振、频谱、统计和波形相似等方面。基于能量属性的长短时窗能量比(STA/LTA)方法原理简单,容易实现,并且能够满足实时处理的要求是微地震监测数据处理常用的方法,但该方法的缺点是基于单道能量特征在低信噪比事件的识别往往无法取得令人满意的效果。基于波形相似特征的微地震事件识别方法包括利用事件内道间波形相似的方法以及模板匹配识别的方法,能够提高低信噪比事件信号的识别能力,是目前微地震事件识别的重要手段。模板匹配的微地震事件识别方法依赖波形模板的选取,而基于事件内道间波形相似的方法相比之下更具普适性的特点,但其存在计算效率低和规则干扰误拾率高的问题。
微地震初至拾取方法可以分为基于单道记录的方法和基于多道记录的方法。基于单道记录的初至拾取方法可以通过计算和分析信号的时频分布、瞬时属性、能量属性、偏振属性等来实现,应用在低信噪比实际数据时所得初至结果往往会出现个别道的误差较大。基于多道记录的互相关方法充分利用了道间记录的相似特征,能够提高初至拾取结果的质量。然而,当微地震记录信噪比低或者存在多震相信号时,互相关函数的峰值不明显或者连续出现多个大小近于相似的峰值,此时最大值对应的时间并不能可靠地代表真正的目标震相位置,校正后的初至结果并不理想。
发明内容
本发明要解决的技术问题在于,针对现有技术的上述目标震相位置不可靠的缺陷,提供一种基于到时曲线拟合的微地震震相识别初至拾取方法、一种计算机设备和一种存储有计算机可读指令的存储介质。
本发明解决其技术问题所采用的技术方案是:构造一种基于到时曲线拟合的微地震震相识别初至拾取方法,包括:
步骤1,利用射孔记录的震相到时建立到时曲线拟合公式;
步骤2,沿着拟合曲线选取滑动窗口计算相似系数及叠加道能量,以能量加权相似系数值作为微地震信号存在与否的依据;
步骤3,基于最优拟合曲线所得叠加道计算时差校正后各道记录的剩余时差,设定约束时窗,判断各道剩余时差校正量是否小于给定时窗大小;
步骤4,叠加满足条件的剩余时差校正后的波形记录,在叠加道初至到时的基础上得到各道的准确初至到时,不满足约束条件的各道记录到时通过插值得到。
在本发明所述的基于到时曲线拟合的微地震震相识别初至拾取方法中,在利用射孔记录的震相到时建立到时曲线拟合公式的步骤中,包括分析射孔记录的震相到时规律,选用合适的拟合公式并给定拟合参数范围的步骤。
在本发明所述的基于到时曲线拟合的微地震震相识别初至拾取方法中,初至曲线拟合公式为
t=t0+K×dT×m
其中t0为相对记录时刻,m为检波器号,K为拟合曲线拟合曲线参数,dT为采样时间间隔,K×dT可以视为道间时差。圆圈和三角形分别为P波S波初至信息,实线和虚线为对应的线性拟合结果。
在本发明所述的基于到时曲线拟合的微地震震相识别初至拾取方法中,沿着拟合曲线选取滑动窗口计算相似系数及叠加道能量,以能量加权相似系数值作为微地震信号存在与否的依据的步骤中,包括:
利用拟合曲线沿着时间轴在连续的地震记录上逐个采样点进行移动,同时通过拟合参数不断地修改曲线;
在移动过程中计算出以不同拟合曲线为起点的一定时窗范围内的多道相似系数,多道相似系数的计算公式为:
式中M为道数,P为拟合参数,
分别为各道时窗的起点与终点,
为时窗长度,A
x,A
y,A
z分别为时窗内地震波三个分量的振幅,S
x,S
y,S
z分别为三个分量记录的多道相似系数;此时定义总多道相似系数为
当拟合曲线接近真实初至时,有效微地震信号沿着拟合曲线所得叠加道波形能量增强,随机噪声能量压制,所得叠加道表示为
其中tm为拟合曲线初至时刻;
预先给定一个能量加权相似系数的阈值,如果所得能量加权相似系数大于阈值,则在微地震记录的相应时刻处附近存在一个有效的微地震信号且拟合曲线反映的是真实初至曲线趋势。
在本发明所述的基于到时曲线拟合的微地震震相识别初至拾取方法中,为了避免非微地震信号的误拾取,定义能量加权相似系数I计算公式为I=Eq·S
其中指数q取值为0~1,E为沿着拟合曲线所得叠加道的能量,计算公式表示为
在本发明所述的基于到时曲线拟合的微地震震相识别初至拾取方法中,基于最优拟合曲线所得叠加道计算时差校正后各道记录的剩余时差,设定约束时窗,判断各道剩余时差校正量是否小于给定时窗大小的步骤中,包括:
通过震相识别得到震相时刻以及初至曲线趋势T(m);
基于已得震相时刻以及初至曲线趋势叠加时差校正后的记录形成叠加道,求取叠加道与时差校正后的记录相对时差;
设定约束时窗δT,判断各道时差校正量|ΔT(m)|是否小于δT。
在本发明所述的基于到时曲线拟合的微地震震相识别初至拾取方法中,叠加满足条件的剩余时差校正后的波形记录,在叠加道初至到时的基础上得到各道的准确初至到时,不满足约束条件的各道记录到时通过插值得到的步骤中,包括:
叠加满足时窗约束条件的剩余时差校正后的波形记录,利用长短时窗能量比初至拾取方法可以得到最终叠加道的相对初至到时;
满足约束条件的各道记录的准确到时可以公式T′(m)=T(m)+ΔT(m)+T0计算;不满足约束条件的各道记录到时通过对T′插值得到。
在本发明所述的基于到时曲线拟合的微地震震相识别初至拾取方法中,在得到各道的准确初至到时的步骤之后,还包括对拾取结果进行验证的步骤,采用STA/LTA方法对实际数据进行初至拾取处理,并进行对比。
本发明解决其技术问题所采用的技术方案是:构造一种计算机设备,包括输入输出单元、存储器和处理器,所述存储器中存储有计算机可读指令,所述计算机可读指令被所述处理器执行时,使得所述处理器执行如前述技术方案所述的微地震震相识别初至拾取方法中的步骤。
本发明解决其技术问题所采用的技术方案是:构造一种存储有计算机可读指令的存储介质,所述计算机可读指令被一个或多个处理器执行时,使得一个或多个处理器执行如前述技术方案所述的微地震震相识别初至拾取方法中的步骤。
本发明中的基于到时曲线拟合的微地震震相识别及初至拾取方法,沿着拟合曲线计算能量加权相似系数进行信号识别,能够快速有效识别出微地震信号(P波S波震相),已知的震相初至趋势的基础上,计算目标震相的剩余时差校正量,设定约束时窗可以有效避免低信噪比的道记录初至拾取出现较大的误差情况,提高初至拾取的效率和准确度,这对于改进压裂施工效率,提高油气藏改造的成功率具有重要的意义。
附图说明
下面将结合附图及实施例对本发明作进一步说明,附图中:
图1是本发明提供的一种基于到时曲线拟合的微地震震相识别初至拾取方法的流程示意图;
图2为本发明的一具体实例中射孔记录及其震相到时拟合示意图;
图3为本发明的一具体实例中基于到时曲线拟合的微地震信号识别示意图;
图4为本发明的一具体实例中基于到时曲线趋势的先验信息的初至优化示意图;
图5为本发明的一具体实例中本发明方法所得初至结果与STA/LTA所得结果对比图。
具体实施方式
为了对本发明的技术特征、目的和效果有更加清楚的理解,现对照附图详细说明本发明的具体实施方式。
如图1所示,图1为本发明的基于多道扫描叠加的井下微地震事件自动识别方法的流程图。
参照图1,在步骤110,利用射孔记录或者强能量事件信号的震相到时建立到时曲线拟合公式。
对于均匀速度模型以及垂直排列的观测系统,检波器接收地震波的初至曲线可以表示为抛物线,当监测距离r远大于检波器间距Δz且检波器组布设在压裂段一侧时,此时初至曲线可近似为一次曲线
t=t0+K×dT×m
其中k,d为一次曲线拟合参数,d为相对时间。通过采用一次曲线拟合的方法得到曲线同样能够近似反映真实到时曲线趋势。
井中微地震监测仪器空间位置相邻,检波器接收到微地震震相在记录上表现出相似的走时特征。射孔初至信息一定程度上能够表示该压裂段的微地震事件震相曲线趋势。如图2所示,图中左右两侧图分别为一实例射孔记录及初至拟合结果,该地区监测区域地层起伏变化相对较小,监测距离远大于检波器间距,记录中同相轴道间时差变化不大。根据射孔记录初至的走时规律可以通过一次曲线拟合的方法得到近似的初至曲线。拟合公式为
t=t0+K×dT×m
其中t0为相对记录时刻,m为检波器号,K为拟合曲线拟合曲线参数,dT为采样时间间隔,K×dT可以视为道间时差。圆圈和三角形分别为P波S波初至信息,实线和虚线为对应的线性拟合结果。
在步骤120,沿着拟合曲线选取滑动窗口计算相似系数及叠加道能量的步骤可包括如下:选用一个滑动时窗的长度;利用拟合曲线沿着时间轴在连续的地震记录上逐个采样点进行移动,同时通过拟合参数不断地修改曲线。
利用拟合曲线沿着时间轴在连续的地震记录上逐个采样点进行移动,同时通过拟合参数不断地修改曲线。在移动过程中计算出以不同拟合曲线为起点的一定时窗范围内的多道相似系数,多道相似系数的计算公式为:
式中M为道数,P为拟合参数,
分别为各道时窗的起点与终点,
为时窗长度,A
x,A
y,A
z分别为时窗内地震波三个分量的振幅,S
x,S
y,S
z分别为三个分量记录的多道相似系数。此时定义总多道相似系数为
当拟合曲线接近真实初至时,有效微地震信号沿着拟合曲线所得叠加道波形能量增强,随机噪声能量压制,所得叠加道可以表示为
其中tm为拟合曲线初至时刻。多道相似系数未考虑到信号的能量特征,所以对于某些能量弱的相关环境噪声也较为敏感。为了避免这些非微地震信号的误拾取,定义能量加权相似系数I计算公式为
I=Eq·S
其中指数q取值为0~1,E为沿着拟合曲线所得叠加道的能量,计算公式可以表示为
预先给定一个能量加权相似系数的阈值,如果所得能量加权相似系数大于该阈值,则可认为在微地震记录的相应时刻处附近存在一个有效的微地震信号且拟合曲线反映的是真实初至曲线趋势。如图2所示,从图2中可以看出当时窗内存在有效信号,由于有效信号具有一定的道间相似性,则该时窗内的多道相似系数较大,且有效微地震信号沿着拟合曲线所得叠加道波形能量增强,随机噪声能量压制。从图2中能量加权相似系数剖面中可以清楚看出存在着1个局部极大值区域,可识别出微地震震相时刻t及拟合参数对应的位置Kmax。
在步骤130,由于识别所得震相的时刻和同相轴趋势曲线只能大概的反映出地震波的初至位置,因此仍需要对震相初至进行进一步优化。在到时曲线的先验信息基础上以一定时窗约束初至拾取可以避免低信噪比道初至的误拾。步骤包括:基于最优拟合曲线得到时差校正记录及叠加道;计算叠加道与时差校正后各道记录的剩余时差;设定约束时窗,判断各道剩余时差校正量是否小于给定时窗大小。将微地震信号识别时刻以及对应的拟合参数带入公式可以计算得到微地震同相轴的到时趋势曲线,为了得到更为准确的初至信息,需要利用叠加道求取各道初至到时的剩余时差校正量,如图4所示。初步时差校正后的记录如图4a,图4a中所示为P波震相Z分量记录的时差校正记录,然后利用叠加道计算各道相对剩余时差得到时差校正后记录和最终叠加道波形。设定约束时窗δT,判断各道时差校正量|ΔT(m)|是否小于δT。如果是,则保留该道i=m;否则j=m。
在步骤140,叠加满足条件的剩余时差校正后的波形记录(即i=m),通过长短能量比方法可以得到叠加道的相对初至到时T0,在叠加道初至到时的基础上得到各道的准确初至到时
T(m)=T(m)+ΔT(m)+T0
=t+Kmax·dT·m+ΔT(m)+T0
不满足约束条件的各道记录到时通过插值得到。
为了与这一方法的应用效果进行对比,同时也采用了STA/LTA方法对该实际数据进行了初至拾取处理。图5a为采用本发明方法所得初至结果,以初至结果时差校正后的记录如图5c所示,图5b为采用STA/LTA方法所得初至结果,以初至结果时差校正后的记录如图5d所示。时差校正后的记录中波形对齐可以说明初至拾取的效果,通过对比两种方法得到的初至结果可以看出,利用能量比法求出的初至时间在低信噪比道存在一定的误差,本发明提出的方法利用了到时曲线规律的先验信息利用互相关求取剩余时差能够有效地避免低信噪比情况下初至误差的出现,且可在较为合理的时窗内得到更为准确的初至到时。
此外本发明还提供了一种计算机设备,包括输入输出单元、存储器和处理器,所述存储器中存储有计算机可读指令,所述计算机可读指令被所述处理器执行时,使得所述处理器执行如前述技术方案所述的微地震震相识别初至拾取方法中的步骤。
本发明提供了一种存储有计算机可读指令的存储介质,所述计算机可读指令被一个或多个处理器执行时,使得一个或多个处理器执行如前述技术方案所述的微地震震相识别初至拾取方法中的步骤。
本发明中的基于到时曲线拟合的微地震震相识别及初至拾取方法,沿着拟合曲线计算能量加权相似系数进行信号识别,能够快速有效识别出微地震信号(P波S波震相),已知的震相初至趋势的基础上,计算目标震相的剩余时差校正量,设定约束时窗可以有效避免低信噪比的道记录初至拾取出现较大的误差情况,提高初至拾取的效率和准确度,这对于改进压裂施工效率,提高油气藏改造的成功率具有重要的意义。
上面结合附图对本发明的实施例进行了描述,但是本发明并不局限于上述的具体实施方式,上述的具体实施方式仅仅是示意性的,而不是限制性的,本领域的普通技术人员在本发明的启示下,在不脱离本发明宗旨和权利要求所保护的范围情况下,还可做出很多形式,这些均属于本发明的保护之内。