一种微地震地面监测速度模型校正算法的制作方法

本发明涉及一种微地震地面监测速度模型校正算法。

背景技术：

水力压裂技术在低渗油气田开发开采以及储层改造过程中发挥着至关重要的作用，通过水力压裂技术开发低渗透油气藏已成为一种普遍趋势。微地震定位技术是微地震监测工作的核心，该技术实施过程中，裂缝延伸导致周围岩石破裂，从而引发一系列可观测记录的微地震事件。而监测工区范围内的速度模型是影响微地震事件定位准确与否的主要因素，因此，如何获得一个有效的速度模型是微地震监测工程中一个关键性问题。

在微地震监测工程中，地震层析成像方法是寻求地下速度结构的有效手段，但在实际压裂监测过程中受条件制约，射孔信息量较少，地面检波器数量及覆盖范围严重不足，获得较为精细的工区地下速度结构几乎是不可能的。目前现有的微地震速度模型构建方法都以射孔重定位精度作为评判依据。anikiev采用尝试的办法，同时对各层测井初始模型进行同时增加或减小一定值，得到一个相对较准确的反演速度模型。pei等采用occam法进行地下速度模型反演，谭玉阳等采用levenberg-marquardt方法进行速度模型反演。pei等，jiang等采用极快速模拟退火方案进行速度模型校正，取得了一定效果。但以上方法均依赖于射孔记录的初至拾取，在低信噪比情况下该方法失效。

cn105807316a公开了一种《基于振幅叠加的地面观测微地震速度模型校正方法》，从逆时偏移振幅叠加法原理出发，并结合极快速模拟退火算法，不需要拾取震相初至信息，通过监测射孔点重定位精度，判断速度模型是否可用于后续微地震定位。能有效的克服现有算法存在的极大值和最大值不分的情况，准确的找到能量聚焦最大值e，准确定位射孔位置，有效的校正速度模型。

但上述方法仅考虑到各层的速度不确定性，而并没有考虑到各层位位置不确定性，导致射孔重定位仍存在一定误差。因此，如何有效克服现有的微地震地面速度模型校正缺陷从而进一步提高微地震事件定位精度是本领域亟需解决的问题。

技术实现要素：

为了解决上述问题，本发明在振幅叠加地面观测微地震速度模型校正方法的基础上，在极快速模拟退火方法过程中同时考虑各层速度不确定性以及层位不确定性，通过扩大解空间的方法提高射孔事件重定位精度。

实现本发明方法的步骤如下：

a、在地面布置n个检波器，并以射孔点位置作为中心点定义目标区域；

b、根据定位精度需要将目标区域划分成若干网格尺寸大小相同的体元；

c、根据测井曲线建立初始速度模型并读取n个检波器所获取的n道地震波资料；

d、采用震幅叠加方法定位射孔体元，并获得能量聚焦值最大的网格体元中心点坐标及相应能量聚焦值e；

e、以速度值和层位值作为不确定因素，采用极快速模拟退火法对目标层进行调整；

f、对射孔事件进行重定位，如果定位误差不满足定位要求，还需进行回火处理调整速度模型，直到定位误差足够小为止。

g、是否满足射孔精度要求：是；

h、结束。

步骤d中震幅叠加方法具体如下：

d1、计算各道检波器获取的地震波相对于参考道的走时差δti；参考道选择有相对清晰的初至同相轴及较高的信噪比；

式(1)表示第i道与参考道走时之差；式(1)中δlij为第i道与参考道在第j层的路径差，vj第j层的波速；

d2、各道检波器获取的地震波波形根据δti逆时偏移后进行震幅叠加，表达式为：

式(2)中，各道波形在第j时刻的振幅值为a(f(t),k)；m为检波器个数，n为时窗长度；k为第j时刻的第k个采样点；

将式(1)带入式(2)得到：

其中步骤e，极快速模拟退火法具体实施如下：

第一步，根据测井曲线获得初始速度向量v,v＝[vp1,vp2,vp3,...,vpn]及初始层位模型向量h，h＝[h1,h2,h3,...,hm-1]，其中，vpi为第i层p波速度，hj为第j个层界面深度坐标。

第二步，计算目标函数e(v),具体计算方法如下：

在测井数据中，选择任意一道作为参考道m，采用射线追踪的方法求取其它各道相对于该参考道的理论走时差：

δt^cal＝[t1-tm,t2-tm,...,tn-tm](4)

根据该理论走时差对各道检波器获得的射孔记录波形进行逆时偏移叠加，目标函数(叠加能量)的数学表达式如下：

其中a为第i个体元在第j时刻的震幅大小，第i个体元的中心点坐标为(xi,yi,zi)，m为检波器个数，n为时窗长度；

第三步，计算模拟退火初始温度t0，获取初始温度的解决方法如下：

先给初始温度一个很小的正值，然后不断乘以一个恒大于1的数，直到满足对任何模型的接受概率接近于1为止；

其中，v0为初始速度模型，v1为第一次迭代模型结果。

第四步，进行模拟退火降温计算，极快速模拟退火降温公式如下：

tk＝t0exp(-ck^1/2n)(7)

其中k为迭代次数，t0为初始退火温度，c为给定常数用来调整算法退火温度，这里c＝0.5；n为需要调整速度的层数；

在进行模拟退火计算过程中，还需调整速度向量v和层位模型向量h，具体公式为：

其中，与为各层速度的最大最小边界值，其中与为各层位最大最小边界值。其中sgn为符号函数，x为随机变量，取值范围在[-1,1]之间。μ∈[0,1]；产生x表达式为：

每次迭代计算模拟退火接受概率计算如下：当e(v′)≥e(v)时，v′替代v作为当前最优解，当e(v′)＜e(v)时,以概率

替代当前最优解，其中tk为第k次迭代时的温度值。最后以e(v)>e作为本次模拟退火运算叠代终止条件。

本发明的有益效果：

现有技术中的极快速模拟退火算法过程中仅考虑了各层的速度向量v＝(vp1,vp2,vp3,...,vpn)^t的不确定性。而并没有考虑到工区地下层位测量不准而导致的第i道与参考道在第j层的路径差δlij的不确定性。这就导致速度模型解空间较小，虽经模拟退火运算，但仍无法完全消除射孔重定位误差。本专利提出在原算法模拟退火计算过程中加入各层位不确定因素，通过增大解空间的方式提高射孔点重定位精度，能够很好的解决射孔重定位不准的问题，并有效提高射孔点周围微地震定位可信度。

附图说明

图1基于振幅叠加的地面观测微地震速度模型校正方法流程图；

图2二维层状起伏结构地层模型图；

图3检波器各道正演模拟波形图；

图4(a)常规振幅叠加算法射孔重定位结果；

图4(b)算法改进后的射孔重定位结果；

图5震源定位分布结果；

图6震源定位误差结果图；

图7(a),(b)为经过本文改进后射孔重定位结果；

图7(c),(d)为采用常规振幅叠加速度模型校正算法射孔重定位结果。

具体实施方式

实施例1

下面结合附图和模拟地层模型对本发明技术方案以具体实施例的方式进行清楚、完整的描述。具体实施步骤如图1所示，具体如下：

a、如图2所示，选择二维层状起伏结构地层模型进行验证本文的改进算法，实际信号采用有限差分波动方程进行模拟，根据所模拟的地层模型建立初始平层状地层模型，后续的正演采用射线追踪方法，反演方法采用上述振幅叠加速度模型校正方法，并加入本文所提出的层位不确定因素，本次实验假定中间各层位上下浮动不超过50m。检波器排布与震源位置如图2所示，其中首个检波器坐标为(200,0)，每个检波器间距为50m，震源位置为(50,1400)。表1为根据测井曲线建立的初始速度模型参数。图3为各道检波器正演模拟波形。

b、将目标区域划分为网格尺寸为10m的体元，每一个网格中心可以被看作是微地震事件发生的潜在位置；

c、如图2所示，根据测井信息(图2绿线所示)获得各个地层的速度模型，初始速度向量v＝[2200,2750,3400,4800](单位m/s)，初始层界面向量模型h＝[424,729,1230](单位：m)，其中1～4层的深度范围依次为0m～424m，424m～779m，779m～1230m，1230m～1600m，与此对应的各层速度最大边界值为v^max＝[2500,3200,4000,5500](单位m/s)，各层速度最小边界值为v^min＝[1800,2400,3000,4500](单位m/s)；并读取17道检波器模拟资料。

d、采用振幅叠加对射孔定位并获得能量聚焦值最大的网格中心点坐标及相应能量聚焦值e＝12.7；震幅叠加方法具体如下：

d1、计算各道检波器获取的地震波相对于参考道的走时差δti；参考道选择有相对清晰的初至同相轴及较高的信噪比；

式(1)表示第i道与参考道走时之差；式(1)中δlij为第i道与参考道在第j层的路径差，vj第j层的波速；

d2、各道检波器获取的地震波波形根据δti逆时偏移后进行震幅叠加，表达式为：

式(2)中，各道波形在第j时刻的振幅值为a(f(t),k)；m为检波器个数，n为时窗长度；k为第j时刻的第k个采样点；

将式(1)带入式(2)得到：

e、将第一道检波器(200,0)作为参考道，根据测井信息建立的初始速度模型，利用射线追踪技术计算模拟射孔点到各参考道的理论时差δt^cal，根据理论时差对各道进行偏移叠加，获取目标函数e(v)。

δt^cal＝[t1-tm,t2-tm,...,tn-tm](4)

根据该理论走时差对各道检波器获得的射孔记录波形进行逆时偏移叠加，目标函数(叠加能量)的数学表达式如下：

其中a为第i个体元在第j时刻的震幅大小，第i个体元的中心点坐标为(xi,yi,zi)，m为检波器个数，n为时窗长度；

f、采用极快速模拟退火算法对各层界面位置及各层速度值进行调整，首先按照上述模拟退火步骤e中第三步求取模拟退火算法初始温度t0，

计算模拟退火初始温度t0，获取初始温度的解决方法如下：

先给初始温度一个很小的正值，然后不断乘以一个恒大于1的数，直到满足对任何模型的接受概率接近于1为止；

其中，v0为初始速度模型，v1为第一次迭代模型结果。得到初始温度的近似值为100℃。

然后进行模拟退火降温计算，所采用的极快速模拟退火降温公式如下：

tk＝t0exp(-ck^1/2n)(7)

其中k为迭代次数，t0为初始退火温度，c为给定常数用来调整算法退火温度，这里c＝0.5；n为需要调整速度的层数；

在进行模拟退火计算过程中，还需调整速度模型向量v和层位模型向量h，具体公式为：

其中，与为各层速度的最大最小边界值，其中与为各层位最大最小边界值。其中sgn为符号函数，x为随机变量，取值范围在[-1,1]之间。μ∈[0,1]；产生x表达式为：

每次迭代计算模拟退火接受概率计算如下：当e(v′)≥e(v)时，v′替代v作为当前最优解，当e(v′)＜e(v)时,以概率

替代当前最优解，其中tk为第k次迭代时的温度值。最后以e(v)>e作为本次模拟退火运算叠代终止条件。

g、将经过调整后的速度模型对射孔事件进行重定位，若定位结果满足定位精度(误差在5m以内)要求，即与真实射孔位置相差甚微，；若定位结果不满足定位精度要求，即与真实射孔位置相差较大，则需对进行回火处理继续调整速度模型，直至满足射孔精度要求后结束。

如图7所示，利用常规振幅叠加速度模型构建算法得到的射孔重定位结果误差为23m，而采用本文改进后的算法重定位误差小于2m，可见本文改进后的算法比常规算法更优。

为验证该方法对震源定位可信度的提高，本专利在射孔周围设置了24个震源(图5所示)，分别采用初始速度模型，原震幅叠加微地震速度模型校正算法，以及改进后速度模型校正方法构建速度模型，并对24个震源进行重定位处理，定位分布结果如图5所示，震源定位分布误差如图6所示。从以上结果可以看出，采用本专利方法进行速度模型校正后，能够有效提高射孔周边微地震定位可信度。