一种确定地震波走时的方法及装置制造方法及图纸

本申请实施例公开了一种确定地震波走时的方法及装置。所述方法包括：根据第一地层介质的传播速度和炮点到第一地层介质与第二地层介质之间的分界面的距离，确定第一地层介质对应的直达波走时，根据第二地层介质的传播速度、炮点到检波点的距离和炮点到分界面的距离，确定第二地层介质对应的直达波走时，根据第一地层介质的传播速度和地层介质的厚度，确定第一地层介质对应的垂直双程旅行时；基于第一地层介质对应的直达波走时和第二地层介质对应的直达波走时，分别确定多个走时关联系数；根据多个走时关联系数和垂直双程旅行时，确定目标地震波走时。本申请实施例提供的技术方案，可以提高在横向变速介质的情况下计算地震波走时的效率。

【技术实现步骤摘要】
一种确定地震波走时的方法及装置
本申请涉及油田的勘探、开发和开采
，特别涉及一种确定地震波走时的方法及装置。
技术介绍
随着地震勘探技术的发展，油气勘探的难度和深度也越来越大，对地震资料的信噪比和分辨率要求也越来越高。在地震数据采集时，由震源(炮点)激发，会产生一个地震子波，这个地震子波在地下介质中向下传播，当遇到一个反射界面时，地震子波会发生反射和透射，一部分地震子波的能量反射后改变方向，向上传播，到达地面接收点，并被布设在地面接收点的检波器接收。地下不同反射界面所发生的反射，依次会向地表传播，在地面的检波器会依次接收到这些地下反射的地震子波，构成了检波器记录的地震数据。地震波走时是指地震波从震源传到接收点所经过的时间。而地震波走时计算是地震数据处理技术中一个非常重要的技术，在地震数据处理方法中，许多技术都涉及到走时计算，例如，叠后时间偏移、叠后深度偏移、叠前时间偏移、叠前深度偏移、动校正、速度分析、层析反演、正演模拟等技术。其中，偏移处理的过程主要是根据已经建立的地下反射界面的深度-速度模型，计算地震子波从炮点传播到地下反射界面再返回到地面接收点的走时，然后，根据所计算的走时和地震数据记录到的地震子波，反演地下反射界面的真实位置。因此，走时计算是偏移处理最为关键的技术。目前，常规偏移处理(包括叠前、叠后、时间和深度偏移等)中的地震波走时通常采用射线追踪方法来计算。该方法的主要过程是首先对地下反射界面的深度-速度模型进行网格化处理，然后在每个网格上进行走时计算。然而，针对水平层状介质横向速度变化的情况，该方法采用的网格较小，计算费时，不能满足实际数据处理需要。
技术实现思路
本申请实施例的目的是提供一种确定地震波走时的方法及装置，以提高在横向变速介质的情况下计算地震波走时的效率。为解决上述技术问题，本申请实施例提供一种确定地震波走时的方法及装置是这样实现的：一种确定地震波走时的方法，提供有目的工区中地层反射界面的速度与深度模型；其中，所述速度与深度模型中包括沿水平方向相邻且地震波传播速度不同的第一地层介质和第二地层介质，以及地层介质的厚度；所述第一地层介质和所述第二地层介质的上方分别设置有炮点和检波点；所述方法包括：根据所述第一地层介质的传播速度和炮点到所述第一地层介质与所述第二地层介质之间的分界面的距离，确定所述第一地层介质对应的直达波走时，根据所述第二地层介质的传播速度、所述炮点到检波点的距离和所述炮点到所述分界面的距离，确定所述第二地层介质对应的直达波走时，以及根据所述第一地层介质的传播速度和所述地层介质的厚度，确定所述第一地层介质对应的垂直双程旅行时；基于所述第一地层介质对应的直达波走时和所述第二地层介质对应的直达波走时，分别确定多个走时关联系数；根据所述多个走时关联系数和所述垂直双程旅行时，确定目标地震波走时。优选方案中，分别确定多个走时关联系数，包括：基于所述第一地层介质对应的直达波走时和所述第二地层介质对应的直达波走时，分别确定多个直达波交叉走时；根据所述多个直达波交叉走时，分别所述确定多个走时关联系数。优选方案中，所述多个直达波交叉走时中包括零阶直达波交叉走时、一阶直达波交叉走时、二阶直达波交叉走时、三阶直达波交叉走时、四阶直达波交叉走时和五阶直达波交叉走时；相应的，采用下述公式分别确定所述多个直达波交叉走时：其中，τn表示所述多个直达波交叉走时中的第n阶直达波交叉走时，n＝0,1,2,3,4,5，当n＝0时，τ0表示所述零阶直达波交叉走时，当n＝1时，τ1表示所述一阶直达波交叉走时，当n＝2时，τ2表示所二阶直达波交叉走时，n＝3时，τ3表示所述三阶直达波交叉走时，n＝4时，τ4表示所述四阶直达波交叉走时，n＝5时，τ5表示所述五阶直达波交叉走时；td表示所述第一地层介质对应的直达波走时，tx表示所述第二地层介质对应的直达波走时，V1表示所述第一地层介质的传播速度，V2表示所述第二地层介质的传播速度。优选方案中，所述多个走时关联系数中包括零阶偏移走时关联系数、一阶偏移走时关联系数、二阶偏移走时关联系数、三阶偏移走时关联系数、四阶偏移走时关联系数和五阶偏移走时关联系数；相应的，分别确定多个走时关联系数，包括：根据所述零阶直达波交叉走时确定所述零阶偏移走时关联系数；根据所述零阶直达波交叉走时和所述一阶直达波交叉走时确定所述一阶偏移走时关联系数；根据所述零阶直达波交叉走时、所述一阶直达波交叉走时和所述二阶直达波交叉走时，确定所述二阶偏移走时关联系数；根据所述零阶直达波交叉走时、所述一阶直达波交叉走时、所述二阶直达波交叉走时和所述三阶直达波交叉走时，确定所述三阶偏移走时关联系数；根据所述零阶直达波交叉走时、所述一阶直达波交叉走时、所述二阶直达波交叉走时、所述三阶直达波交叉走时和所述四阶直达波交叉走时，确定所述四阶偏移走时关联系数；根据所述零阶直达波交叉走时、所述一阶直达波交叉走时、所述二阶直达波交叉走时、所述三阶直达波交叉走时、所述四阶直达波交叉走时和所述五阶偏移走时关联系数，确定所述五阶偏移走时关联系数。优选方案中，采用下述公式确定所述零阶偏移走时关联系数：其中，β0表示所述零阶偏移走时关联系数，τ0表示所述零阶直达波交叉走时；采用下述公式确定所述一阶偏移走时关联系数：其中，β1表示所述一阶偏移走时关联系数，τ1表示所述一阶直达波交叉走时；采用下述公式确定所述二阶偏移走时关联系数：其中，β2表示所述二阶偏移走时关联系数，τ2表示所述二阶直达波交叉走时；采用下述公式确定所述三阶偏移走时关联系数：其中，β3表示所述三阶偏移走时关联系数，τ3表示所述三阶直达波交叉走时；采用下述公式确定所述四阶偏移走时关联系数：其中，β4表示所述四阶偏移走时关联系数，τ4表示所述四阶直达波交叉走时；采用下述公式确定所述五阶偏移走时关联系数：其中，β5表示所述五阶偏移走时关联系数，τ5表示所述五阶直达波交叉走时。优选方案中，采用下述公式确定所述第一地层介质对应的垂直双程旅行时：其中，t0表示所述第一地层介质对应的垂直双程旅行时，V1表示所述第一地层介质的传播速度，h表示所述地层介质的厚度。优选方案中，采用下述公式确定所述第一地层介质对应的直达波走时：其中，td表示所述第一地层介质对应的直达波走时，V1表示所述第一地层介质的传播速度，d表示所述炮点到所述第一地层介质与所述第二地层介质之间的分界面的距离。优选方案中，采用下述公式确定所述第二地层介质对应的直达波走时：其中，tx表示所述第二地层介质对应的直达波走时，V2表示所述第二地层介质的传播速度，x表示所述炮点到检波点的距离，d表示所述炮点到所述第一地层介质与所述第二地层介质之间的分界面的距离。优选方案中，采用下述公式确定目标地震波走时：其中，t表示所述目标地震波走时，β0、β1、β2、β3、β4和β5分别表示所述多个走时关联系数中包括的零阶偏移走时关联系数、一阶偏移走时关联系数、二阶偏移走时关联系数、三阶偏移走时关联系数、四阶偏移走时关联系数和五阶偏移走时关联系数，t0表示所述第一地层介质对应的垂直双程旅行时。一种确定地震波走时的装置，所述装置提供目的工区中地层反射界面的速度与深度模型；其中，所述速度与深度模型中包括沿水平方向相邻且地震波传播速度不本文档来自技高网...

【技术保护点】
1.一种确定地震波走时的方法，其特征在于，提供有目的工区中地层反射界面的速度与深度模型；其中，所述速度与深度模型中包括沿水平方向相邻且地震波传播速度不同的第一地层介质和第二地层介质，以及地层介质的厚度；所述第一地层介质和所述第二地层介质的上方分别设置有炮点和检波点；所述方法包括：根据所述第一地层介质的传播速度和炮点到所述第一地层介质与所述第二地层介质之间的分界面的距离，确定所述第一地层介质对应的直达波走时，根据所述第二地层介质的传播速度、所述炮点到检波点的距离和所述炮点到所述分界面的距离，确定所述第二地层介质对应的直达波走时，以及根据所述第一地层介质的传播速度和所述地层介质的厚度，确定所述第一地层介质对应的垂直双程旅行时；基于所述第一地层介质对应的直达波走时和所述第二地层介质对应的直达波走时，分别确定多个走时关联系数；根据所述多个走时关联系数和所述垂直双程旅行时，确定目标地震波走时。

【技术特征摘要】
1.一种确定地震波走时的方法，其特征在于，提供有目的工区中地层反射界面的速度与深度模型；其中，所述速度与深度模型中包括沿水平方向相邻且地震波传播速度不同的第一地层介质和第二地层介质，以及地层介质的厚度；所述第一地层介质和所述第二地层介质的上方分别设置有炮点和检波点；所述方法包括：根据所述第一地层介质的传播速度和炮点到所述第一地层介质与所述第二地层介质之间的分界面的距离，确定所述第一地层介质对应的直达波走时，根据所述第二地层介质的传播速度、所述炮点到检波点的距离和所述炮点到所述分界面的距离，确定所述第二地层介质对应的直达波走时，以及根据所述第一地层介质的传播速度和所述地层介质的厚度，确定所述第一地层介质对应的垂直双程旅行时；基于所述第一地层介质对应的直达波走时和所述第二地层介质对应的直达波走时，分别确定多个走时关联系数；根据所述多个走时关联系数和所述垂直双程旅行时，确定目标地震波走时。2.根据权利要求1所述的方法，其特征在于，分别确定多个走时关联系数，包括：基于所述第一地层介质对应的直达波走时和所述第二地层介质对应的直达波走时，分别确定多个直达波交叉走时；根据所述多个直达波交叉走时，分别所述确定多个走时关联系数。3.根据权利要求2所述的方法，其特征在于，所述多个直达波交叉走时中包括零阶直达波交叉走时、一阶直达波交叉走时、二阶直达波交叉走时、三阶直达波交叉走时、四阶直达波交叉走时和五阶直达波交叉走时；相应的，采用下述公式分别确定所述多个直达波交叉走时：其中，τn表示所述多个直达波交叉走时中的第n阶直达波交叉走时，n＝0,1,2,3,4,5，当n＝0时，τ0表示所述零阶直达波交叉走时，当n＝1时，τ1表示所述一阶直达波交叉走时，当n＝2时，τ2表示所二阶直达波交叉走时，n＝3时，τ3表示所述三阶直达波交叉走时，n＝4时，τ4表示所述四阶直达波交叉走时，n＝5时，τ5表示所述五阶直达波交叉走时；td表示所述第一地层介质对应的直达波走时，tx表示所述第二地层介质对应的直达波走时，V1表示所述第一地层介质的传播速度，V2表示所述第二地层介质的传播速度。4.根据权利要求3所述的方法，其特征在于，所述多个走时关联系数中包括零阶偏移走时关联系数、一阶偏移走时关联系数、二阶偏移走时关联系数、三阶偏移走时关联系数、四阶偏移走时关联系数和五阶偏移走时关联系数；相应的，分别确定多个走时关联系数，包括：根据所述零阶直达波交叉走时确定所述零阶偏移走时关联系数；根据所述零阶直达波交叉走时和所述一阶直达波交叉走时确定所述一阶偏移走时关联系数；根据所述零阶直达波交叉走时、所述一阶直达波交叉走时和所述二阶直达波交叉走时，确定所述二阶偏移走时关联系数；根据所述零阶直达波交叉走时、所述一阶直达波交叉走时、所述二阶直达波交叉走时和所述三阶直达波交叉走时，确定所述三阶偏移走时关联系数；根据所述零阶直达波交叉走时、所述一阶直达波交叉走时、所述二阶直达波交叉走时、所述三阶直达波交叉走时和所述四阶直达波交叉走时，确定所述四阶偏移走时关联系数；根据所述零阶直达波交叉走时、所述一阶直达波交叉走时、所述二阶直达波交叉走时、所述三阶直达波交叉走时、所述四阶直达波交叉走时和所述五阶偏移走时关联系数，确定所述五阶偏移走时关联系数。5.根据权利要求4所述的方...

【专利技术属性】
技术研发人员：高少武，王成祥，李建峰，张红英，王狮虎，张巍毅，
申请(专利权)人：中国石油天然气集团有限公司，中国石油集团东方地球物理勘探有限责任公司，
类型：发明
国别省市：北京,11