首页 / 专利库 / 物理 / 机械波 / 地震波 / 一种利用斜井直达地震波走时计算地层地震波速度的方法

一种利用斜井直达地震波走时计算地层地震波速度的方法

阅读:626发布:2020-05-12

专利汇可以提供一种利用斜井直达地震波走时计算地层地震波速度的方法专利检索,专利查询,专利分析的服务。并且本 发明 公开了一种利用斜井直达 地震 波 走时计算 地层 地震波 速度的方法,该方法假设激发地震波的 位置 为炮点S,其二维坐标为(xs,zs),沿钻井井轨迹有N个地震波接收点,这些接收点从上到下按从小到大的顺序依次记为R1,R2,R3,…,RN,相应的接收点坐标记为(x1,z1),(x2,z2),(x3,z3),…,(xN,zN),各接收点所接收到直达地震波的实际走时依次记为t1,t2,t3,…,tN。本发明适用于斜井、并可估算出光滑的地层地震 波速 度。,下面是一种利用斜井直达地震波走时计算地层地震波速度的方法专利的具体信息内容。

1.一种利用斜井直达地震波走时计算地层地震波速度的方法,其特征在于,包括以下步骤:
1)设激发地震波的位置为炮点S,其二维坐标为(xs,zs),沿钻井井轨迹有N个地震波接收点,这些接收点从上到下按从小到大的顺序依次记为R1,R2,R3,...,RN,相应的接收点坐标记为(x1,z1),(x2,z2),(x3,z3),...,(xN,zN),各接收点所接收到直达地震波的实际走时依次记为t1,t2,t3,...,tN;
设地下有N个平地层,从上到下按从小到大的顺序从1开始依次编号,各地层底界面纵坐标从上到下依次对应接收点的纵坐标z1,z2,z3,...,zN,各地层速度从上到下依次为v1,v2,v3,...,vN;
2)按地震波以直线传播的假设从上到下逐层计算地层速度,计算公式为:
其中,Lj表示从炮点S到接收点Rj的地震波传播路径在第j个地层里的长度,其计算公式为:
3)对以上计算得到的层速度vi,采用2M+1点滑动窗口平均进行光滑处理,其中i=1,2,
3,...,N,计算公式如下:
其中,M的大小反映计算结果的光滑程度,M的值在M<(N-1)/2条件下任意选取;令迭代次数iter的值为1;
4)根据地震波Snell定律计算在速度vi下的理论直达地震波走时Ti,并计算理论走时Ti和实际走时ti的均方根误差rms_error,其中i=1,2,3,...N;理论直达地震波走时Ti的计算方法采用射线追踪法;
5)设定一个正数eps以及一个最大迭代次数Nmax;若计算的理论走时Ti和实际走时ti的均方差rms_error>eps,并且迭代次数iter然后令迭代次数iter的值增1并返回步骤4);否则,输出vi(i=1,2,3,...,N)作为最终计算结果,计算结束。
2.根据权利要求1所述的利用斜井直达地震波走时计算地层地震波速度的方法,其特征在于,步骤4)中,追踪过程采用密集射线打靶法:首先从炮点S开始,在沿着以炮点为圆心、以炮点S到接收点Ri方向为中心的90°扇形范围内,以小度间隔发射密集的地震波射线,射线穿过各地层界面后的出射角θi按以下Snell定律计算
其中,θi-1表示射线从第i-1层进入第i层的入射角;射线在各地层界面的交点坐标(Xi,Zi)采用以下两式计算:
Zi=zi, i=1,2,3,...,N     (6)
然后选取这些射线与垂线x=xi的交点到Ri最近的射线作为成功射线,并根据下式计算地震波走时
其中,Lj由下式计算
由上可计算得到各接收点直达地震波走时以及相应的射线路径在各地层界面上的交点坐标;理论走时Ti和实际走时ti的均方根误差rms_error用下式计算:

说明书全文

一种利用斜井直达地震波走时计算地层地震波速度的方法

技术领域

[0001] 本发明属于地震勘探技术领域,涉及一种计算地层地震波速度的方法,尤其是一种利用斜井直达地震波走时计算地层地震波速度的方法。

背景技术

[0002] 地震勘探领域,井中地震方法在钻井井壁放置检波器接收人工激发的地震波,根据所接收的地震波传播时间(也称为走时)可估计地层的地震波传播速度或用于反射波数据处理分析,进而可研究地层的孔隙度、含油气、泊松比等物理性质,应用于油气勘探和开发。目前利用斜井直达地震波走时估计地层速度的方法一般基于垂直井假设。

发明内容

[0003] 本发明的目的在于克服上述现有技术的缺点,提供一种利用斜井直达地震波走时计算地层地震波速度的方法。
[0004] 本发明的目的是通过以下技术方案来实现的:
[0005] 这种利用斜井直达地震波走时计算地层地震波速度的方法,包括以下步骤:
[0006] 1)设激发地震波的位置为炮点S,其二维坐标为(xs,zs),沿钻井井轨迹有N个地震波接收点,这些接收点从上到下按从小到大的顺序依次记为R1,R2,R3,...,RN,相应的接收点坐标记为(x1,z1),(x2,z2),(x3,z3),...,(xN,zN),各接收点所接收到直达地震波的实际走时依次记为t1,t2,t3,...,tN;
[0007] 设地下有N个平地层,从上到下按从小到大的顺序从1开始依次编号,各地层底界面纵坐标从上到下依次对应接收点的纵坐标z1,z2,z3,...,zN,各地层速度从上到下依次为v1,v2,v3,...,vN;
[0008] 2)按地震波以直线传播的假设从上到下逐层计算地层速度,计算公式为:
[0009]
[0010] 其中,Lj表示从炮点S到接收点Rj的地震波传播路径(也称为地震射线)在第j个地层里的长度,其计算公式为:
[0011]
[0012] 3)对以上计算得到的层速度vi(i=1,2,3,...,N)采用(2M+1)点滑动窗口平均进行光滑处理,计算公式如下:
[0013]
[0014] 其中,M的大小反映计算结果的光滑程度,M的值在M<(N-1)/2条件下任意选取,计算时根据实际情况确定;令迭代次数iter的值为1;
[0015] 4)根据地震波Snell定律计算在速度vi(i=1,2,3,...,N)下的理论直达地震波走时Ti(i=1,2,3,...,N),并计算理论走时Ti(i=1,2,3,...,N)和实际走时ti(i=1,2,3,...,N)的均方根误差rms_error;理论直达地震波走时的计算方法采用射线追踪法,追踪过程采用密集射线打靶法;首先从炮点S开始,在沿着以炮点为圆心、以炮点S到接收点Ri方向为中心的90°扇形范围内,以小度(如0.01°)间隔发射密集的地震波射线,射线穿过各地层界面后的出射角θi按以下Snell定律计算
[0016]
[0017] 其中,θi-1表示射线从第i-1层进入第i层的入射角;射线在各地层界面的交点坐标(Xi,Zi)采用以下两式计算:
[0018]
[0019] Zi=zi,i=1,2,3,...,N  (6)
[0020] 然后选取这些射线与垂线x=xi的交点到Ri最近的射线作为成功射线,并根据下式计算地震波走时
[0021]
[0022] 其中,Lj由下式计算
[0023]
[0024] 由上可计算得到各接收点直达地震波走时以及相应的射线路径在各地层界面上的交点坐标。理论走时Ti(i=1,2,3,...,N)和实际走时ti(i=1,2,3,...,N)的均方根误差rms_error用下式计算
[0025]
[0026] 5)设定一个正数eps以及一个最大迭代次数Nmax,例如eps=0.01及Nmax=100。若计算的理论走时Ti(i=1,2,3,...,N)和实际走时ti(i=1,2,3,...,N)的均方差rms_error>eps,并且迭代次数iter
[0027]
[0028] 然后令迭代次数iter的值增1并返回步骤4);否则,输出vi(i=1,2,3,...,N)作为最终计算结果,计算结束。
[0029] 本发明具有以下有益效果:
[0030] 本发明公开一种利用斜井直达地震波走时、适用于斜井、并可估算出光滑的地层地震波速度的方法。附图说明
[0031] 图1为井中地震波观测示意图;
[0032] 图2为井中实测地震波记录图;
[0033] 图3(a)为拾取的实际直达地震波走时,图3(b)为计算的地层速度曲线
[0034] 图4为利用所计算地层速度处理得到的反射地震波成像剖面图。

具体实施方式

[0035] 下面结合附图对本发明做进一步详细描述:
[0036] 如图1所示,假设激发地震波的位置为炮点S,其二维坐标为(xs,zs),沿钻井井轨迹有N个地震波接收点,这些接收点从上到下按从小到大的顺序依次记为R1,R2,R3,...,RN,相应的接收点坐标记为(x1,z1),(x2,z2),(x3,z3),...,(xN,zN),各接收点所接收到直达地震波的实际走时依次记为t1,t2,t3,...,tN。
[0037] 1.假设地下有N个水平地层,从上到下按从小到大的顺序从1开始依次编号,各地层底界面纵坐标从上到下依次对应接收点的纵坐标z1,z2,z3,...,zN,各地层速度从上到下依次为v1,v2,v3,...,vN。
[0038] 2.按地震波以直线传播的假设从上到下逐层计算地层速度,计算公式为[0039]
[0040] 其中,Lj表示从炮点S到接收点Rj的地震波传播路径(也称为地震射线)在第j个地层里的长度,其计算公式为
[0041]
[0042] 3.对以上计算得到的层速度vi(i=1,2,3,...,N)采用(2M+1)点滑动窗口平均进行光滑处理,计算公式如下
[0043]
[0044] 其中,M的大小反映计算结果的光滑程度,其值可在M<(N-1)/2条件下任意选取,计算时根据实际情况确定。令迭代次数iter的值为1。
[0045] 4.根据地震波Snell定律计算在速度vi(i=1,2,3,...,N)下的理论直达地震波走时Ti(i=1,2,3,...,N),并计算理论走时Ti(i=1,2,3,...,N)和实际走时ti(i=1,2,3,...,N)的均方根误差rms_error。理论直达地震波走时的计算方法采用射线追踪法,追踪过程采用密集射线打靶法。首先从炮点S开始,在沿着以炮点为圆心、以炮点S到接收点Ri方向为中心的90°扇形范围内,以小角度(如0.01°)间隔发射密集的地震波射线,射线穿过各地层界面后的出射角θi按以下Snell定律计算
[0046]
[0047] 其中,θi-1表示射线从第i-1层进入第i层的入射角,入射角与出射角θi的图示见图1。射线在各地层界面的交点坐标(Xi,Zi)采用以下两式计算
[0048]
[0049] Zi=zi,i=1,2,3,...,N  (6)
[0050] 然后选取这些射线与垂线x=xi的交点到Ri最近的射线作为成功射线,并根据下式计算地震波走时
[0051]
[0052] 其中,Lj由下式计算
[0053]
[0054] 由上可计算得到各接收点直达地震波走时以及相应的射线路径在各地层界面上的交点坐标。理论走时Ti(i=1,2,3,...,N)和实际走时ti(i=1,2,3,...,N)的均方根误差rms_error用下式计算
[0055]
[0056] 5.设定一个正数eps以及一个最大迭代次数Nmax,例如eps=0.01及Nmax=100。若计算的理论走时Ti(i=1,2,3,...,N)和实际走时ti(i=1,2,3,...,N)的均方差rms_error>eps,并且迭代次数iter
[0057]
[0058] 然后令迭代次数iter的值增1并返回第4步骤。否则,输出vi(i=1,2,3,...,N)作为最终计算结果,计算结束。
[0059] 如图2是在斜井Hb001中接收到的地震记录。该记录相应的炮点位于地表,距离Hb001井井口136m,共有206个接收点,接收点沿斜井井壁以20m等间距分布于斜深100~4200m之间。据该地震记录可拾取到各接收点的实际直达地震波走时(见图3a)。由拾取的实际走时按本发明计算得到井下各地层地震波速度(见图3b)。该速度曲线较为光滑,可直接用于反射地震数据成像处理,所得成像结果见图4。据实钻情况,该成像结果成功预测了石炭系岩含气层的横向变化,显示了本发明的应用效果。
高效检索全球专利

专利汇是专利免费检索,专利查询,专利分析-国家发明专利查询检索分析平台,是提供专利分析,专利查询,专利检索等数据服务功能的知识产权数据服务商。

我们的产品包含105个国家的1.26亿组数据,免费查、免费专利分析。

申请试用

分析报告

专利汇分析报告产品可以对行业情报数据进行梳理分析,涉及维度包括行业专利基本状况分析、地域分析、技术分析、发明人分析、申请人分析、专利权人分析、失效分析、核心专利分析、法律分析、研发重点分析、企业专利处境分析、技术处境分析、专利寿命分析、企业定位分析、引证分析等超过60个分析角度,系统通过AI智能系统对图表进行解读,只需1分钟,一键生成行业专利分析报告。

申请试用

QQ群二维码
意见反馈