发明内容
本发明示例性实施例的目的在于提供一种基于垂直地震剖面资料的地层速度确定方法,以克服现有技术中无法获取目标区域钻井中最深地震波接收点以下的地层速度的问题。
根据本发明的示例性实施例,提供一种基于垂直地震剖面资料的地层速度确定方法,包括:(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)当第一组地层速度中与第二组地层速度中对应地层的地层速度的差值的绝对值中的最大值小于预定阈值时,将第一组地层速度作为目标区域深部各地层的待处理地层速度;当第一组地层速度中与第二组地层速度中对应地层的地层速度的差值的绝对值中的最大值不小于预定阈值时,将第二组地层速度作为当前获取的目标区域深部的各地层的地层速度,重复步骤(D1)至(D4)。
在根据本发明示例性实施例的基于垂直地震剖面资料的地层速度确定方法中,能够获取目标区域钻井中最深地震波接收点以下地层的地震波速度,利用所述地震波速度可以有效地提高了垂直地震剖面反射地震波数据成像的精度。
具体实施方式
现将详细参照本发明的示例性实施例,所述实施例的示例在附图中示出,其中,相同的标号始终指的是相同的部件。以下将通过参照附图来说明所述实施例,以便解释本发明。
图1示出根据本发明示例性实施例的基于垂直地震剖面资料的地层速度确定方法的流程图。这里,作为示例,所述方法可由基于垂直地震剖面资料确定地层速度的设备来实现,也可完全通过计算机程序来实现。
如图1所示,在步骤S100,利用垂直地震剖面资料获取目标区域浅部各地层的待处理地层速度,其中,目标区域被划分为浅部和深部,其中,浅部指示目标区域钻井中的临界地震波接收点以上的区域;深部指示目标区域钻井中的临界地震波接收点以下的预定区域。
这里,垂直地震剖面(Verticalseismicprofile)是相对于地面地震剖面而言的,该方法是在钻井中观测地震波场,通过将检波器设置于钻井中的不同深度来记录地面震源所产生的地震信号,因此,上述地震波接收点可指示设置于钻井中的可接收地震信号的检波器,此外,作为示例,所述临界地震波接收点可指示设置在钻井最深处的可接收地震信号的检波器。
作为示例,为了获取目标区域浅部各地层的待处理地层速度,可利用设置于目标区域浅部钻井中的每个地震波接收点的深度信息和直达地震波走时确定与目标区域浅部钻井中的每个地震波接收点对应的地层的待处理地层速度,其中,所述地震波沿直线传播。
具体说来,图2示出根据本发明示例性实施例的目标区域地层结构的示意图。如图2所示,设激发地震波的位置为炮点s,其二维坐标为(xs,zs),沿钻井井轨迹有N个地震波接收点,这些地震波接收点从上到下按从小到大的顺序依次记为R1、R2...RN,相应的地震波接收点坐标记为(0,z1),(0,z2),(0,z3),...,(0,zN),各地震波接收点所接收到直达地震波的实际走时依次记为t1,、t2、...tN,其中,目标区域浅部钻井中的每个地震波接收点的深度信息和直达地震波走时可从垂直地震剖面资料中获取。这里,由于目标区域浅部钻井中的每个地震波接收点的深度信息等于与其对应的地层的深度信息,因此,可利用每个地震波接收点的深度信息和直达地震波走时来获取目标区域浅部钻井中的每个地震波接收点对应的地层的待处理地层速度。
例如,可利用下面的等式1确定与目标区域浅部钻井中的每个地震波接收点对应的地层的待处理地层速度:
等式1
其中,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表示临界地震波接收点。
此外,如图2所示,还可假设在目标区域深部存在多个地震波接收点,例如,设地下有N+M个地层,在每个地层上依然设置有地震波接收点,并且,每个地层的深度信息分别与在地层深部假设存在的相应的地震波接收点的深部信息对应。
再次参照图1,在步骤S200,获取目标区域深部中假设存在的多个地震波接收点所对应的反射地震波同相轴的斜率。以下将结合图3来具体描述根据本发明示例性实施例的确定目标区域深部假设存在的多个地震波接收点所对应的反射地震波同相轴的斜率的步骤。
如图3所示,在步骤S210,获取目标区域深部中各个假设存在的地震波接收点所接收的直达地震波走时。
例如,可利用下面的等式2获取目标区域深部中各个假设存在的地震波接收点所接收的直达地震波走时:
等式2
其中,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的正整数。
这里,可利用下面的等式3确定目标区域深部第N+j个地震波接收点RN+j所对应的地震反射线在第N+j地层的长度:
等式3
其中,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所对应的地层的深度信息。
在步骤S220,根据在步骤S210获取的各个假设存在的地震波接收点所接收的直达地震波走时确定在每个假设存在的地震波接收点所对应的地层上产生的反射地震波到达临界地震波接收点处的反射时间。
例如,可利用下面的等式4来确定在每个假设存在的地震波接收点所对应的地层上产生的反射地震波到达临界地震波接收点处的反射时间:
TNN+j=tN+2(tN+j-tN),等式4
其中,TNN+j表示在第N+j个假设存在的地震波接收点RN+j所对应的地层上产生的反射地震波到达临界地震波接收点RN的反射时间。
在步骤S230,利用在每个假设存在的地震波接收点所对应的地层上产生的反射地震波到达临界地震波接收点处的反射时间以及从垂直地震剖面资料上获取的在每个假设存在的地震波接收点所对应的地层上产生的反射地震波到达目标区域浅部钻井中的第一个地震波接收点的反射时间确定每个假设存在的地震波接收点所对应的反射波同相轴的斜率。
例如,可利用下面的等式5确定所述每个假设存在的地震波接收点所对应的反射地震波同相轴的斜率:
等式5
其中,kN+j表示目标区域深部假设存在的第N+j个地震波接收点RN+j所对应的反射地震波同相轴的斜率,T1N+j表示在第N+j个假设存在的第N+j个地震波接收点RN+j所对应的地层上产生的反射地震波到达目标区域浅部钻井中的第一个地震波接收点R1的反射时间,这里,T1N+j可从垂直地震剖面资料中获取。
再次参照图1,在步骤S300,利用所述反射地震波同相轴的斜率确定目标区域深部各地层的地层速度。
这里,作为示例,可对各个假设存在的地震波接收点所对应的反射地震波同相轴的斜率取倒数,将其作为目标区域深部与各个假设存在的地震波接收点相应的各地层的地层速度。
在步骤S400,对目标区域深部各地层的地层速度进行稳定性校正处理,以获取目标区域深部各地层的待处理地层速度。以下将结合图4来具体描述根据本发明示例性实施例的获取目标区域深部各地层的待处理地层速度的步骤。
如图4所示,在步骤S410,将当前获取的目标区域深部的各地层的地层速度作为第一组地层速度。
在步骤S420,重新获取目标区域中假设存在的多个地震波接收点所对应的反射地震波同相轴的斜率,利用所述反射地震波同相轴的斜率确定目标区域深部各地层的地层速度,并将重新获取的目标区域深部各地层的地层速度作为第二组地层速度,
其中,可利用下面的等式6代替上述等式4来重新获取的在第N+j个假设存在的第N+j个地震波接收点RN+j所对应的地层上产生的反射地震波到达临界地震波接收点RN的反射时间:
等式6
其中,TN'N+j表示利用等式6获取的在第N+j个假设存在的第N+j个地震波接收点RN+j所对应的地层上产生的反射地震波波到达临界地震波接收点RN的反射时间。
在步骤S430,将第一组地层速度中与第二组地层速度中对应地层的地层速度的差值的绝对值中的最大值与预定阈值进行比较。
当第一组地层速度中与第二组地层速度中对应地层的地层速度的差值的绝对值中的最大值小于预定阈值时,在步骤S440,将第一组地层速度作为目标区域深部各地层的待处理地层速度;当第一组地层速度中与第二组地层速度中对应地层的地层速度的差值的绝对值中的最大值不小于预定阈值时,在步骤S450,将第二组地层速度作为当前获取的目标区域深部的各地层的地层速度,返回执行步骤S410。
再次参照图1,在步骤S500,对目标区域浅部和深部的所有待处理地层速度进行光滑处理,以得到目标区域各地层的最终地层速度。例如,可采用(2P+1)点滑动窗口平均进行光滑处理,其中,P为整数,这里,P的取值越大,所得到的最终地层速度越光滑。
图5示出根据本发明示例性实施例的目标区域的地层速度沿地层深度变化的曲线示意图。如图5所示,目标区域的临界地震波接收点的深度信息为3082.5m,该深度以下的地层速度可采用本发明的方法确定。
图6示出将采用本发明的方法得到地层速度用于垂直地震剖面反射地震波成像处理的效果图,其中,图6中的(a)示出将采用本发明得到的地层速度用于垂直地震剖面反射地震波成像处理的剖面图,图6中的(b)示出同一目标区域的地面三维地震数据成像的剖面图。由图6中的(a)和(b)可见,目标区域的垂直地震剖面的反射地震波成像结果的主要层位的展布形态与地面三维地震数据剖面相接近,多数层位在深度上具有一致性,这说明通过本发明可以得到垂直地震剖面的临界地震波接收点以下地层的地震波速度,利用所述地震波速度可以有效地提高垂直地震剖面反射地震波数据成像的精度。
综上所述,在根据本发明示例性实施例的基于垂直地震剖面资料的地层速度确定方法中,能够获取目标区域钻井中最深地震波接收点以下地层的地震波速度,利用所述地震波速度可以有效地提高了垂直地震剖面反射地震波数据成像的精度。
尽管已经参照其示例性实施例具体显示和描述了本发明,但是本领域的技术人员应该理解,在不脱离权利要求所限定的本发明的精神和范围的情况下,可以对其进行形式和细节上的各种改变。