的反射地震波同相轴的 斜率； (C) 利用所述反射地震波同相轴的斜率确定目标区域深部各地层的地层速度； (D) 对目标区域深部各地层的地层速度进行稳定性校正处理，以获取目标区域深部各 地层的待处理地层速度； (E) 对目标区域浅部和深部的所有待处理地层速度进行光滑处理，以得到目标区域各 地层的最终地层速度。2. 如权利要求1所述的方法，其中，步骤（A)包括：利用设置于目标区域浅部钻井中的 每个地震波接收点的深度信息和直达地震波走时确定与目标区域浅部钻井中的每个地震 波接收点对应的地层的待处理地层速度，其中，目标区域浅部钻井中的每个地震波接收点 的深度信息等于与其对应的地层的深度信息，其中，所述地震波沿直线传播。3. 如权利要求2所述的方法，其中，在步骤（A)中，利用下面的等式确定与目标区域浅 部钻井中的每个地震波接收点对应的地层的待处理地层速度：其中，\表示目标区域浅部的第i地层的待处理地层速度，z:表示目标区域浅部的第1地层的深度信息，h表示目标区域浅部钻井中的第一个地震波接收点Ri所接收的直达地震 波走时，21表示目标区域浅部的第i地层的深度信息，zii表示目标区域浅部的第i-Ι地层 的深度信息，^表示目标区域浅部钻井中的第i个地震波接收点1所接收的直达地震波走 时，hi表示目标区域浅部钻井中的第i-Ι个地震波接收点Li所接收的直达地震波走时， 其中，第i个地震波接收点民与第i地层相对应，1彡i彡N，N为正整数，RN表示临界地震 波接收点。4. 如权利要求1所述的方法，其中，步骤⑶包括： (B1)获取目标区域深部中各个假设存在的地震波接收点所接收的直达地震波走时； (B2)根据在步骤（B1)获取的各个假设存在的地震波接收点所接收的直达地震波走时 确定在每个假设存在的地震波接收点所对应的地层上产生的反射地震波到达临界地震波 接收点处的反射时间； (B3)利用在每个假设存在的地震波接收点所对应的地层上产生的反射地震波到达临 界地震波接收点处的反射时间以及从垂直地震剖面资料获取的在每个假设存在的地震波 接收点所对应的地层上产生的反射地震波到达目标区域浅部钻井中的第一个地震波接收 点处的反射时间，确定每个假设存在的地震波接收点所对应的反射地震波同相轴的斜率。5. 如权利要求4所述的方法，其中，在步骤（B1)中，利用下面的等式获取目标区域深部 中各个假设存在的地震波接收点所接收的直达地震波走时：其中，tN+]表示目标区域深部假设存在的第N+j个地震波接收点RN+]所接收的直达地震 波走时，tN表示临界地震波接收点RN所接收的直达地震波走时，LN+]表示目标区域深部假设 存在的第N+j个地震波接收点RN+]所对应的地震反射线在第N+j地层的长度，vN+]表示目标 区域深部假设存在的第N+j个地震波接收点RN+j所对应的地层速度，其中，vN+j=vN，^表 示临界地震波接收点RN所对应的地层速度，Μ表示大于等于1的正整数。6. 如权利要求5所述的方法，其中，利用下面的等式确定目标区域深部假设存在的第 N+j个地震波接收点RN+j所对应的地震反射线在第N+j地层的长度：其中，zN+]表示目标区域深部假设存在的第N+j个地震波接收点RN+]所对应的第N+j地 层的深度信息，i表示目标区域深部假设存在的第Ν+j-l个地震波接收点RN+]i所对应 第Ν+j-l地层的深度信息，xs表示地震波激发点s在地面沿X轴方向的坐标信息，其中，zN+] =zN+j(ZN-ZNi)，1 <j<M，其中，示临界地震波接收点R派对应的地层的深度信息， zNi表示第N-1个地震波接收点RNi所对应的地层的深度信息。7. 如权利要求6所述的方法，其中，在步骤（B3)中，利用下面的等式确定每个假设存在 的地震波接收点所对应的反射地震波同相轴的斜率：其中，kN+]表示目标区域深部假设存在的第N+j个地震波接收点RN+]所对应的反射地 震波同相轴的斜率，TlN+j表示在第N+j个假设存在的第N+j个地震波接收点RN+j所对应的 地层上产生的反射地震波到达目标区域浅部钻井中的第一个地震波接收点&的反射时间， TNN+j表示在第N+j个假设存在的第N+j个地震波接收点RN+j所对应的地层上产生的反射地 震波波到达临界地震波接收点RN的反射时间，其中，TNN+j=tN+2(tN+j-tN)。8. 如权利要求1所述的方法，其中，在步骤（C)中，对各个假设存在的地震波接收点所 对应的反射地震波同相轴的斜率取倒数，将其作为目标区域深部与各个假设存在的地震波 接收点相应的各地层的地层速度。9. 如权利要求7所述的方法，其中，步骤⑶包括： (D1)将当前获取的目标区域深部的各地层的地层速度作为第一组地层速度； (D2)重新获取目标区域中假设存在的多个地震波接收点所对应的反射地震波同相轴 的斜率，利用所述反射地震波同相轴的斜率确定目标区域深部各地层的地层速度，并将重 新获取的目标区域深部各地层的地层速度作为第二组地层速度， 其中，利用下面的等式获取在假设存在的第N+j个地震波接收点&+]所对应的地层上 产生的反射地震波到达临界地震波接收点RN的反射时间：其中，TN'N+j表示重新获取的在第N+j个假设存在的第N+j个地震波接收点RN+j所对应 的地层上产生的反射地震波波到达临界地震波接收点RN的反射时间； (D3)将第一组地层速度中与第二组地层速度中对应地层的地层速度的差值的绝对值 中的最大值与预定阈值进行比较； (D4)当第一组地层速度中与第二组地层速度中对应地层的地层速度的差值的绝对值 中的最大值小于预定阈值时，将第一组地层速度作为目标区域深部各地层的待处理地层速 度；当第一组地层速度中与第二组地层速度中对应地层的地层速度的差值的绝对值中的最 大值不小于预定阈值时，将第二组地层速度作为当前获取的目标区域深部的各地层的地层 速度，重复步骤（D1)至（D4)。
【专利摘要】提供一种基于垂直地震剖面资料的地层速度确定方法，包括：(A)利用垂直地震剖面资料获取目标区域浅部各地层的待处理地层速度；(B)获取目标区域深部中假设存在的多个地震波接收点所对应的反射地震波同相轴的斜率；(C)利用所述反射地震波同相轴的斜率确定目标区域深部各地层的地层速度；(D)对目标区域深部各地层的地层速度进行稳定性校正处理，以获取目标区域深部各地层的待处理地层速度；(E)对目标区域浅部和深部的所有待处理地层速度进行光滑处理，以得到目标区域各地层的最终地层速度。根据所述方法，能够获取目标区域钻井中最深地震波接收点以下地层的地震波速度，可以有效地提高了垂直地震剖面反射地震波数据成像的精度。
【IPC分类】G01V1/28, G01V1/42
【公开号】CN105259578
【申请号】CN201510700222
【发明人】罗坤, 秦俐, 曹立斌, 李豪, 唐建
【申请人】中国石油集团川庆钻探工程有限公司地球物理勘探公司
【公开日】2016年1月20日
【申请日】2015年10月26日