专利名称:共转换点抽道集方法
技术领域:
本发明涉及多分量转换波地震数据的成像领域,尤其涉及一种共转换点(CCP, Common Converted Point profile)抽道集方法。
背景技术:
在多分量转换波地震数据的成像技术领域,主要研究的技术方法有以下三大方 面转换波的CCP叠加成像、转换波叠前(后)时间、转换波叠前(后)深度偏移。其中,转换波的CCP叠加成像主要难题是转换波的上、下行波的不对称性,使得 常规的基于共中心点(CMP)道集的速度分析、动校正和水平叠加技术不能适用于转换波, 转换波成像必须研究新的技术。目前,在这方面技术方法的研究已有一定程度的发展,主要 体现在以下几点1、转换波转换点的计算精度提高。由于转换波上下行波传播路径的非对称性,转换点的位置偏向于接收点一侧, 而且随着深度、偏移距以及速度变化而变化。基于共转换点(CCP)道集成像的关键是求 取转换点坐标。当地下介质为水平层状时,一种现有技术(Tessmer G, Behle A. Common reflection point data-stacking technique for converted wave.Geophysical prospecting, 1988, 36 (7) :671 688)提出一个等效的单一均勻水平层模型导出的一个四 次解析表达式很好地近似转换点的坐标。另一现有技术(周竹生,王卫华.一种快速、高 精度的共转换点轨迹计算方法.石油地球物理勘探,1993,28(1) :37 45)提出了计算转 换点坐标的迭代法。再一现有技术(许士勇,马在田.快速有效的转换波共转换点叠加技 术.地球物理学报,2002,45 (4) :557 568)给出了四次方程的近似表达式,使其复杂程度 和计算量都大大减小。再一现有技术(姚陈,于光明,蔡明刚.倾角CDP和倾角CCP叠加.第 21届中国地球物理年会,2005)给出了三维倾角转换波共转换点的解析表达式,使得抽取 共转换点道集的精度提高了一个台阶。2、ACP道集到CCP道集的进步。传统的转换波抽道集的方法是通过一个工区的固定速度比,计算转换点的渐近 线,然后将整个工区的转换波数据抽取到渐近线道集中,即ACP道集,这种方法能一定程度 上能够满足深部地层的成像要求,但是浅层的成像效果较差。目前国际上发展了以精确速 度比计算得到的CCP道集,使得转换波的共转换点道集(CCP道集)从浅到深都能聚焦成像。3、速度分析精度的提高。传统的速度分析方法是必须先获得纵波速度谱,然后利用转换波的CCP道集获得 速度比谱,然后利用新的速度比迭代一次,重抽CCP道集,在进行速度分析来提高精度,这 种方法的缺点是效率地下。现在常用的方法是通过转换波CCP道集直接分析等效C波(转 换波)速度,然后利用新的速度谱重抽道集,这种方法将抽道集与速度分析的步骤分离开。对于三维三分量地震勘探,常规转换波CCP抽道集技术都要首先求解转换点的位置,然后根据转换点的位置和纵横波速度比来确定接收点的位置,再根据双平方根方程将 炮集转换波反射时间段的地震数据映射到CCP道集中,最后基于纵波速度谱,采用速度比 谱扫描的方式确定精确速度比范围,用精确速度比重新抽道集。但是,这种方式导致CCP叠 加成像方法速度慢、且精度差。
发明内容
有鉴于现有技术中的CCP叠加成像方法速度慢、且精度差的问题而做出本发明, 为此本发明的主要目的在于提供一种CCP抽道集方法,其中根据本发明实施例的共转换点CCP抽道集方法包括根据炮点与检波点连线与 CCP面元的交点计算转换点深度;根据转换点深度的深度域速度谱获得转换点的速度比; 根据转换点的速度比计算输入道转换波的反射时间段,将反射时间段内接收的地震数据映 射到CCP面元中,其中,在映射过程中将非双曲时距关系校正成双曲时距关系;生成双曲 精确速度比CCP道集。其中,在根据炮点与检波点连线与CCP面元的交点计算转换点深度之前,该方法 进一步包括计算炮点与检波点连线与CCP面元的交点。优选地,在映射过程中通过以下公式将非双曲时距关系校正成双曲时距关系tc= tl+
VvC ‘其中,t。为输入道转换波的反射时间,tM为自激自收时间,2h为炮检距、Vc为转换 点的纵波速度。其中,在生成双曲精确速度比CCP道集之前,该方法进一步包括使用纵波速度分 析方法分析转换波的等效均方根速度。并且,在使用纵波速度分析方法分析转换波的等效均方根速度之后,该方法进一 步包括获取深度域均方根纵波、横波速度比谱。与现有技术相比,根据本发明实施例的CCP抽道集方法,通过根据面元交点计算 转换点深度,无需迭代计算转换点的位置,提高了处理效率和精度。
此处所说明的附图用来提供对本发明的进一步理解,构成本申请的一部分,本发 明的示意性实施例及其说明用于解释本发明,并不构成对本发明的不当限定。在附图中图1是根据本发明实施例的CCP抽道集方法的流程图;图2是根据本发明实施例的观测系统示意图;图3是根据本发明实施例的三维系统的CCP点示意具体实施例方式本发明的核心思想在于,首先计算根据炮检点连线与CCP面元的交点,并根据该 交点映射反射点深度,并根据深度域速度谱获得深度点的速度比,从而用双平方根方程获 得输入道的反射时间段,直接将一段数据映射到CCP面元中,映射过程中将非双曲时距关 系校正成双曲时距关系,则后续的速度分析可用纵波速度分析方法直接分析转换波的等效均方根速度;利用深度域的纵波与转换波的均方根速度换算成均方根速度比,再用精确速 度比抽CCP道集,校正非双曲,最终精确成像。为使本发明的目的、技术方案和优点更加清楚,以下结合附图及具体实施例,对本 发明作进一步地详细说明。根据本发明的实施例,提供了一种CCP抽道集方法。图1是根据本发明实施例的CCP抽道集方法的流程图,如图1所示,该方法包括步骤102,根据炮点与检波点连线与CCP面元的交点计算转换点深度;步骤104,根据转换点深度的深度域速度谱获得转换点的速度比;步骤106,根据转换点的速度比计算输入道转换波的反射时间段,将反射时间段内 接收的地震数据映射到CCP面元中,其中,在映射过程中将非双曲时距关系校正成双曲时 距关系;
步骤108,生成双曲精确速度比CCP道集。
下面结合图2、图3详细描述上述处理。上述的处理可以具体包括以下7个步骤,
具体地
步骤1 输入一道数据。
如图2所示,炮点(source)、检波点(receiver)的坐标分别为S(xs,ys)、
R(xr, L),它们的中心点为CMP点,转换波的转换点为CP点,其地面投影为C点,炮检距
为2h,炮点与转换点距离为Xe,检波点与转换点距离为&,CMP点与C点距离为D。根 2h
据又=A求取转换点渐近线与炮点、检波点连线的交点C距炮点S的距离其中
2h = ^(xr-xs)2+(yr-ys)2 为炮检距,r = Vs/Vp为横波、纵波速度比。根据S(xs,ys) R(xr, yr)、2h、X。求取 C 点坐标为C(x。,yc)
步骤2 求取炮检连线与成像面元的交点。
(1)如图3所示,求取炮检连线占据的Inline、Crossline矩形区域; Xmax = max (xc, xr),Xmin = min (xc, xr),ymax = max (yc, yr),ymin = min (yc, yr) Inline方向的范围
Xis — Xmin_Xmin% 厶 X Xie = Xmax-Xmax^ Δ X+Δ X 若 xmax% Δχ = 0, Xie = Xie-Ax Crossline方向的范围
Yes ^rnin ^rnin ^ 八 ^
y。e ^max ^max ^ Y+ ^ ^
若 ymax% Δγ = 0, yce = yce- Δ y 其中Δχ、Ay分别为Inline线和Crossline线的间距。 (2)计算C、S连线占据的Inline线、Crossline线的交点坐标,
% Inline方向起始位置
% Inline方向终止位置 % Crossline方向起始位置 % Crossline方向终止位置
5
点,按Inline方向排序。χ :x0, X1, x2, x3, ...,xm, xm+1
xO — Xmin' Xm+1 — Xmax在Xis-Xie范围内检索,若Xi ^ x&&xi+1 ^ χ+Δχ,Ο ^ i彡m则第i,i+1两点连线 位于Inline方向χ,χ+Δχ范围内。取yh = max(yi yi+1),yx = min(yi yi+1)若 Y1 彡 y&&yh ^ y+Δγ,Ο ^ i ^ m 则第 i,i+1两点连线位于Crossline方向y,y+Ay范围内。(Xi, Yi), (xi+1, yi+1)点连线位于 X,χ+Δ χ、y,y+Δ y 确定的面元内。步骤3 根据炮检连线与成像面元的交点计算转换点深度。根据坐标位置可以求取各点与炮点的距离记为Xcd,记对应的转换点的深度为Z。i; 则有
Zci = Xci [(2/2- Xd P -—-Γ ,若求得Zcd深度以上的地层纵波平均速度^ ’可由
ZZZZi \ΛΖ/。, 二 f + f 二 # + # =〔1 + ? j#求取自激自收时间 t0i。接着,将的接收的地震信号置于(Xi,Yi)、(xi+1, yi+1)所在的面元炮检距 为2h的位置处。
6
与Inline线的交点
/ = 5^l与Inline线共有j_l个交点,交点坐标为 Ui =Xis +/Δχ,1</< j-l 1/,( )
Xr ~Xc
其中,m = j-l用于统计交点数。 与Crossline线的交点
ν — y
k= C\y CS % Crossline线共有k_l个交点,交点坐标为 [Xi=Xiy,)
yr-yc
验证交点是否与Inline交点重合,即判断Δ χ是否为0,若不等时m++。 (3)交点重新排列,确定所处面元位置。
炮检连线与Inline线、Crossline线的交点共计m个。加上C、S点,共计m+2个
权利要求
一种共转换点CCP抽道集方法,其特征在于,包括根据炮点与检波点连线与CCP面元的交点计算转换点深度;根据所述转换点深度的深度域速度谱获得所述转换点的速度比;根据所述转换点的速度比计算输入道转换波的反射时间段,将所述反射时间段内接收的地震数据映射到CCP面元中,其中,在映射过程中将非双曲时距关系校正成双曲时距关系;生成双曲精确速度比CCP道集。
2.根据权利要求1所述的方法,其特征在于,在根据炮点与检波点连线与CCP面元的交 点计算转换点深度之前,所述方法进一步包括计算炮点与检波点连线与CCP面元的交点。
3.根据权利要求2所述的方法,其特征在于,在映射过程中通过以下公式将非双曲时 距关系校正成双曲时距关系其中,t。为输入道转换波的反射时间,‘为自激自收时间,2h为炮检距、Vc为转换点的 纵波速度。
4.根据权利要求1所述的方法,其特征在于,在生成双曲精确速度比CCP道集之前,所 述方法进一步包括使用纵波速度分析方法分析转换波的等效均方根速度。
5.根据权利要求4所述的方法,其特征在于,在使用纵波速度分析方法分析转换波的 等效均方根速度之后,所述方法进一步包括获取深度域均方根纵波、横波速度比谱。
全文摘要
本发明公开了一种共转换点抽道集方法,该方法包括根据炮点与检波点连线与CCP面元的交点计算转换点深度;根据转换点深度的深度域速度谱获得转换点的速度比;根据转换点的速度比计算输入道转换波的反射时间段,将反射时间段内接收的地震数据映射到CCP面元中,其中,在映射过程中将非双曲时距关系校正成双曲时距关系;生成双曲精确速度比CCP道集。通过本发明,能够提高处理效率和精度。
文档编号G01V1/36GK101963674SQ201010229078
公开日2011年2月2日 申请日期2010年7月9日 优先权日2010年7月9日
发明者王赟, 芦俊 申请人:中国科学院地质与地球物理研究所