基于垂直地震剖面资料的地层速度确定方法技术

提供一种基于垂直地震剖面资料的地层速度确定方法，包括：(A)利用垂直地震剖面资料获取目标区域浅部各地层的待处理地层速度；(B)获取目标区域深部中假设存在的多个地震波接收点所对应的反射地震波同相轴的斜率；(C)利用所述反射地震波同相轴的斜率确定目标区域深部各地层的地层速度；(D)对目标区域深部各地层的地层速度进行稳定性校正处理，以获取目标区域深部各地层的待处理地层速度；(E)对目标区域浅部和深部的所有待处理地层速度进行光滑处理，以得到目标区域各地层的最终地层速度。根据所述方法，能够获取目标区域钻井中最深地震波接收点以下地层的地震波速度，可以有效地提高了垂直地震剖面反射地震波数据成像的精度。

【技术实现步骤摘要】
基于垂直地震剖面资料的地层速度确定方法
本专利技术总体来说涉及地球物理信号处理领域。更具体地讲，涉及一种基于垂直地震剖面资料的地层速度确定方法。
技术介绍
垂直地震剖面勘探方法是一种钻井中地震观测技术，与地面地震观测技术相比，垂直地震剖面资料的信噪比高，分辨率高，波的运动学和动力学特征明显。垂直地震剖面勘探方法提供了地下地层结构同地面测量参数之间最直接的对应关系，可以为地面地震资料处理解释提供精确的时深转换和速度模型，为零相位子波分析提供支持。通常，垂直地震剖面勘探方法是在地表附近激发地震波而在井下接收地震波响应，由于该方法是在钻井中观测数据，因此，检波器靠近目标层，可观测到上行波和下行波，并且，根据垂直地震剖面所接收的直达地震波走时可估算出地层速度，同时，还可利用所述地层速度进一步对地层岩石性质进行分析和判断，并作为进一步反射地震波成像的基础数据。但是，目前在利用垂直地震剖面所接收的地震波走时来估计地层速度的方法中，大多数方法只考虑地层最深地震波接收点以上的地层速度，并没有包括最深地震波接收点以下的地层速度，因此，利用垂直地震剖面勘探方法所获取的地层速度并不能直接用于最深地震波接收点以下地层部分的反射地震波数据成像。因此，需要一种可直接用于最深地震波接收点以下地层部分的反射地震波数据成像的地层速度获取方法。
技术实现思路
本专利技术示例性实施例的目的在于提供一种基于垂直地震剖面资料的地层速度确定方法，以克服现有技术中无法获取目标区域钻井中最深地震波接收点以下的地层速度的问题。根据本专利技术的示例性实施例，提供一种基于垂直地震剖面资料的地层速度确定方法，包括：(A)利用垂直地震剖面资料获取目标区域浅部各地层的待处理地层速度，其中，目标区域被划分为浅部和深部，其中，浅部指示目标区域钻井中的临界地震波接收点以上的区域；深部指示目标区域钻井中的临界地震波接收点以下的预定区域；(B)获取目标区域深部中假设存在的多个地震波接收点所对应的反射地震波同相轴的斜率；(C)利用所述反射地震波同相轴的斜率确定目标区域深部各地层的地层速度；(D)对目标区域深部各地层的地层速度进行稳定性校正处理，以获取目标区域深部各地层的待处理地层速度；(E)对目标区域浅部和深部的所有待处理地层速度进行光滑处理，以得到目标区域各地层的最终地层速度。可选地，步骤(A)可包括：利用设置于目标区域浅部钻井中的每个地震波接收点的深度信息和直达地震波走时确定与目标区域浅部钻井中的每个地震波接收点对应的地层的待处理地层速度，其中，目标区域浅部钻井中的每个地震波接收点的深度信息等于与其对应的地层的深度信息，其中，所述地震波沿直线传播。可选地，在步骤(A)中，可利用下面的等式确定与目标区域浅部钻井中的每个地震波接收点对应的地层的待处理地层速度：其中，vi表示目标区域浅部的第i地层的待处理地层速度，z1表示目标区域浅部的第1地层的深度信息，t1表示目标区域浅部钻井中的第一个地震波接收点R1所接收的直达地震波走时，zi表示目标区域浅部的第i地层的深度信息，zi-1表示目标区域浅部的第i-1地层的深度信息，ti表示目标区域浅部钻井中的第i个地震波接收点Ri所接收的直达地震波走时，ti-1表示目标区域浅部钻井中的第i-1个地震波接收点Ri-1所接收的直达地震波走时，其中，第i个地震波接收点Ri与第i地层相对应，1≤i≤N，N为正整数，RN表示临界地震波接收点。可选地，步骤(B)可包括：(B1)获取目标区域深部中各个假设存在的地震波接收点所接收的直达地震波走时；(B2)根据在步骤(B1)获取的各个假设存在的地震波接收点所接收的直达地震波走时确定在每个假设存在的地震波接收点所对应的地层上产生的反射地震波到达临界地震波接收点处的反射时间；(B3)利用在每个假设存在的地震波接收点所对应的地层上产生的反射地震波到达临界地震波接收点处的反射时间以及从垂直地震剖面资料获取的在每个假设存在的地震波接收点所对应的地层上产生的反射地震波到达目标区域浅部钻井中的第一个地震波接收点处的反射时间，确定每个假设存在的地震波接收点所对应的反射地震波同相轴的斜率。可选地，在步骤(B1)中，可利用下面的等式获取目标区域深部中各个假设存在的地震波接收点所接收的直达地震波走时：其中，tN+j表示目标区域深部假设存在的第N+j个地震波接收点RN+j所接收的直达地震波走时，tN表示临界地震波接收点RN所接收的直达地震波走时，LN+j表示目标区域深部假设存在的第N+j个地震波接收点RN+j所对应的地震反射线在第N+j地层的长度，vN+j表示目标区域深部假设存在的第N+j个地震波接收点RN+j所对应的地层速度，其中，vN+j＝vN，vN表示临界地震波接收点RN所对应的地层速度，M表示大于等于1的正整数。可选地，可利用下面的等式确定目标区域深部假设存在的第N+j个地震波接收点RN+j所对应的地震反射线在第N+j地层的长度：其中，zN+j表示目标区域深部假设存在的第N+j个地震波接收点RN+j所对应的第N+j地层的深度信息，zN+j-1表示目标区域深部假设存在的第N+j-1个地震波接收点RN+j-1所对应第N+j-1地层的深度信息，xs表示地震波激发点s在地面沿x轴方向的坐标信息，其中，zN+j＝zN+j(zN-zN-1),1≤j≤M，其中，zN表示临界地震波接收点RN所对应的地层的深度信息，zN-1表示第N-1个地震波接收点RN-1所对应的地层的深度信息。可选地，在步骤(B3)中，可利用下面的等式确定每个假设存在的地震波接收点所对应的反射地震波同相轴的斜率：其中，kN+j表示目标区域深部假设存在的第N+j个地震波接收点RN+j所对应的反射地震波同相轴的斜率，T1N+j表示在第N+j个假设存在的第N+j个地震波接收点RN+j所对应的地层上产生的反射地震波到达目标区域浅部钻井中的第一个地震波接收点R1的反射时间，TNN+j表示在第N+j个假设存在的第N+j个地震波接收点RN+j所对应的地层上产生的反射地震波波到达临界地震波接收点RN的反射时间，其中，TNN+j＝tN+2(tN+j-tN)。可选地，在步骤(C)中，可对各个假设存在的地震波接收点所对应的反射地震波同相轴的斜率取倒数，将其作为目标区域深部与各个假设存在的地震波接收点相应的各地层的地层速度。可选地，步骤(D)可包括：(D1)将当前获取的目标区域深部的各地层的地层速度作为第一组地层速度；(D2)重新获取目标区域中假设存在的多个地震波接收点所对应的反射地震波同相轴的斜率，利用所述反射地震波同相轴的斜率确定目标区域深部各地层的地层速度，并将重新获取的目标区域深部各地层的地层速度作为第二组地层速度，其中，利用下面的等式获取在假设存在的第N+j个地震波接收点RN+j所对应的地层上产生的反射地震波到达临界地震波接收点RN的反射时间：其中，TN'N+j表示重新获取的在第N+j个假设存在的第N+j个地震波接收点RN+j所对应的地层上产生的反射地震波波到达临界地震波接收点RN的反射时间；(D3)将第一组地层速度中与第二组地层速度中对应地层的地层速度的差值的绝对值中的最大值与预定阈值进行比较；(D4)当第一组地层速度中与第二组地层速度中对应地层的地层速度的差值的绝对值中的最大值本文档来自技高网...

【技术保护点】
一种基于垂直地震剖面资料的地层速度确定方法，包括：(A)利用垂直地震剖面资料获取目标区域浅部各地层的待处理地层速度，其中，目标区域被划分为浅部和深部，其中，浅部指示目标区域钻井中的临界地震波接收点以上的区域；深部指示目标区域钻井中的临界地震波接收点以下的预定区域；(B)获取目标区域深部中假设存在的多个地震波接收点所对应的反射地震波同相轴的斜率；(C)利用所述反射地震波同相轴的斜率确定目标区域深部各地层的地层速度；(D)对目标区域深部各地层的地层速度进行稳定性校正处理，以获取目标区域深部各地层的待处理地层速度；(E)对目标区域浅部和深部的所有待处理地层速度进行光滑处理，以得到目标区域各地层的最终地层速度。

【技术特征摘要】
1.一种基于垂直地震剖面资料的地层速度确定方法，包括：(A)利用垂直地震剖面资料获取目标区域浅部各地层的待处理地层速度，其中，目标区域被划分为浅部和深部，其中，浅部指示目标区域钻井中的临界地震波接收点以上的区域；深部指示目标区域钻井中的临界地震波接收点以下的预定区域；(B)获取目标区域深部中假设存在的多个地震波接收点所对应的反射地震波同相轴的斜率；(C)利用所述反射地震波同相轴的斜率确定目标区域深部各地层的地层速度；(D)对目标区域深部各地层的地层速度进行稳定性校正处理，以获取目标区域深部各地层的待处理地层速度；(E)对目标区域浅部和深部的所有待处理地层速度进行光滑处理，以得到目标区域各地层的最终地层速度，其中，步骤(B)包括：(B1)获取目标区域深部中各个假设存在的地震波接收点所接收的直达地震波走时；(B2)根据在步骤(B1)获取的各个假设存在的地震波接收点所接收的直达地震波走时确定在每个假设存在的地震波接收点所对应的地层上产生的反射地震波到达临界地震波接收点处的反射时间；(B3)利用在每个假设存在的地震波接收点所对应的地层上产生的反射地震波到达临界地震波接收点处的反射时间以及从垂直地震剖面资料获取的在每个假设存在的地震波接收点所对应的地层上产生的反射地震波到达目标区域浅部钻井中的第一个地震波接收点处的反射时间，确定每个假设存在的地震波接收点所对应的反射地震波同相轴的斜率。2.如权利要求1所述的方法，其中，步骤(A)包括：利用设置于目标区域浅部钻井中的每个地震波接收点的深度信息和直达地震波走时确定与目标区域浅部钻井中的每个地震波接收点对应的地层的待处理地层速度，其中，目标区域浅部钻井中的每个地震波接收点的深度信息等于与其对应的地层的深度信息，其中，所述地震波沿直线传播。3.如权利要求2所述的方法，其中，在步骤(A)中，利用下面的等式确定与目标区域浅部钻井中的每个地震波接收点对应的地层的待处理地层速度：其中，vi表示目标区域浅部的第i地层的待处理地层速度，z1表示目标区域浅部的第1地层的深度信息，t1表示目标区域浅部钻井中的第一个地震波接收点R1所接收的直达地震波走时，zi表示目标区域浅部的第i地层的深度信息，zi-1表示目标区域浅部的第i-1地层的深度信息，ti表示目标区域浅部钻井中的第i个地震波接收点Ri所接收的直达地震波走时，ti-1表示目标区域浅部钻井中的第i-1个地震波接收点Ri-1所接收的直达地震波走时，其中，第i个地震波接收点Ri与第i地层相对应，1≤i≤N，N为正整数，RN表示临界地震波接收点。4.如权利要求1所述的方法，其中，在步骤(B1)中，利用下面的等式获取目标区域深部中各个假设存在的地震波接收点所接收的直达地震波走时：其中，tN+j表示目标区域深部假设存在的第N+j个地震波接收点RN+j所接收的直达地震波走时，tN表示临界地震波接收点RN所接收的直达地震波走时，LN+j表示目标区域深部假设存在的第N+j个地震波接收点RN+j所对应的地震反射线在第N+j地层的长度，vN+j表示目标区域深部假设存在的第N+j个地震波接收点RN+j所对应的地层速...

【专利技术属性】
技术研发人员：罗坤，秦俐，曹立斌，李豪，唐建，
申请(专利权)人：中国石油集团川庆钻探工程有限公司地球物理勘探公司，
类型：发明
国别省市：四川;51