CN101963674A - 共转换点抽道集方法 - Google Patents
共转换点抽道集方法 Download PDFInfo
- Publication number
- CN101963674A CN101963674A CN 201010229078 CN201010229078A CN101963674A CN 101963674 A CN101963674 A CN 101963674A CN 201010229078 CN201010229078 CN 201010229078 CN 201010229078 A CN201010229078 A CN 201010229078A CN 101963674 A CN101963674 A CN 101963674A
- Authority
- CN
- China
- Prior art keywords
- ccp
- point
- velocity
- wave
- transfer point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种共转换点抽道集方法,该方法包括:根据炮点与检波点连线与CCP面元的交点计算转换点深度;根据转换点深度的深度域速度谱获得转换点的速度比;根据转换点的速度比计算输入道转换波的反射时间段,将反射时间段内接收的地震数据映射到CCP面元中,其中,在映射过程中将非双曲时距关系校正成双曲时距关系;生成双曲精确速度比CCP道集。通过本发明,能够提高处理效率和精度。
Description
技术领域
本发明涉及多分量转换波地震数据的成像领域,尤其涉及一种共转换点(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为输入道转换波的反射时间,t0i为自激自收时间,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,yr),它们的中心点为CMP点,转换波的转换点为CP点,其地面投影为C点,炮检距为2h,炮点与转换点距离为Xc,检波点与转换点距离为Xr,CMP点与C点距离为D。根据 求取转换点渐近线与炮点、检波点连线的交点C距炮点S的距离其中: 为炮检距,r=Vs/Vp为横波、纵波速度比。根据S(xs,ys)、R(xr,yr)、2h、Xc求取 C点坐标为:C(xc,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 %Inline方向起始位置
xie=xmax-xmax%Δx+Δx
若xmax%Δx=0,xie=xie-Δx %Inline方向终止位置
Crossline方向的范围:
ycs=ymin-ymin%Δy %Crossline方向起始位置
yce=ymax-ymax%Δy+Δy
若ymax%Δy=0,yce=yce-Δy %Crossline方向终止位置
其中Δx、Δy分别为Inline线和Crossline线的间距。
(2)计算C、S连线占据的Inline线、Crossline线的交点坐标,
与Inline线的交点:
其中,m=j-1用于统计交点数。
与Crossline线的交点:
验证交点是否与Inline交点重合,即判断xi%Δx是否为0,若不等时m++。
(3)交点重新排列,确定所处面元位置。
炮检连线与Inline线、Crossline线的交点共计m个。加上C、S点,共计m+2个点,按Inline方向排序。
x:x0,x1,x2,x3,…,xm,xm+1
x0=xmin,xm+1=xmax
在xis-xie范围内检索,若xi≥x&&xi+1≤x+Δx,0≤i≤m则第i,i+1两点连线位于Inline方向x,x+Δx范围内。
取yh=max(yi,yi+1),yl=min(yi,yi+1)若yl≥y&&yh≤y+Δy,0≤i≤m则第i,i+1两点连线位于Crossline方向y,y+Δy范围内。
(xi,yi)、(xi+1,yi+1)点连线位于x,x+Δx、y,y+Δy确定的面元内。
步骤3:根据炮检连线与成像面元的交点计算转换点深度。
根据坐标位置可以求取各点与炮点的距离记为Xci,记对应的转换点的深度为Zci,则有:
求取自激自收时间t0i。
接着,将t0i-t0(i+1)的接收的地震信号置于(xi,yi)、(xi+1,yi+1)所在的面元炮检距为2h的位置处。
步骤4:初始双曲CCP道集映射。
在上述的步骤中,获得了以下参数:转换点与炮点的距离Xc、炮检距2h、转换点深度Zc,另外输入深度域的纵波速度谱获得转换点处的纵波均方根速度Vp与常速度比r,则输入道的转换波的反射时间为
然后计算出所有面元交点的反射时间,将同一面元的反射时间段的地震数据映射到该反射面元的CCP道集中,在映射过程中,通过下式将非双曲的时间校正为双曲时间:
需要说明,在计算输入道的转换波反射时间时还可以采用转换波的高阶时距方程,该高阶时距方程是双平方根时距方程的高阶近似。
步骤5:用常规纵波数据处理软件(如Foucs、CGG等)分析等效C波均方根速度。
步骤6:深度域速度比谱换算。
将纵波、等效C波均方根速度谱换算至深度域,两者相除,得到深度域均方根纵、横波速度比谱。
步骤7:双曲精确速度比CCP道集生成。
迭代步骤1-6,但是在步骤4时采用精确深度域速度比谱替代固定速度比,完成全部双曲精确速度比CCP道集的生成。
综上所示,根据本发明实施例,通过根据面元交点计算转换点深度,无需迭代计算转换点的位置,提高了处理效率和精度;并且,校正转换波的“非双曲”以后获得的道集可用常规纵波的速度分析方法进行速度分析,比速度比谱分析方法更加精确与快速。
以上所述仅为本发明的实施例而已,并不用于限制本发明,对于本领域 的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (5)
1.一种共转换点CCP抽道集方法,其特征在于,包括:
根据炮点与检波点连线与CCP面元的交点计算转换点深度;
根据所述转换点深度的深度域速度谱获得所述转换点的速度比;
根据所述转换点的速度比计算输入道转换波的反射时间段,将所述反射时间段内接收的地震数据映射到CCP面元中,其中,在映射过程中将非双曲时距关系校正成双曲时距关系;
生成双曲精确速度比CCP道集。
2.根据权利要求1所述的方法,其特征在于,在根据炮点与检波点连线与CCP面元的交点计算转换点深度之前,所述方法进一步包括:
计算炮点与检波点连线与CCP面元的交点。
3.根据权利要求2所述的方法,其特征在于,在映射过程中通过以下公式将非双曲时距关系校正成双曲时距关系:
其中,tc为输入道转换波的反射时间,t0i为自激自收时间,2h为炮检距、Vc为转换点的纵波速度。
4.根据权利要求1所述的方法,其特征在于,在生成双曲精确速度比CCP道集之前,所述方法进一步包括:
使用纵波速度分析方法分析转换波的等效均方根速度。
5.根据权利要求4所述的方法,其特征在于,在使用纵波速度分析方法分析转换波的等效均方根速度之后,所述方法进一步包括:
获取深度域均方根纵波、横波速度比谱。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010229078 CN101963674A (zh) | 2010-07-09 | 2010-07-09 | 共转换点抽道集方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN 201010229078 CN101963674A (zh) | 2010-07-09 | 2010-07-09 | 共转换点抽道集方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN101963674A true CN101963674A (zh) | 2011-02-02 |
Family
ID=43516625
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN 201010229078 Pending CN101963674A (zh) | 2010-07-09 | 2010-07-09 | 共转换点抽道集方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN101963674A (zh) |
Cited By (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103645498A (zh) * | 2013-11-08 | 2014-03-19 | 中国石油大学(北京) | 一种直接生成共转换点道集的方法 |
CN105891881A (zh) * | 2014-05-15 | 2016-08-24 | 蔡晓刚 | 一种水平界面vti介质转换点ccp道集抽取方法 |
CN106019373A (zh) * | 2016-06-24 | 2016-10-12 | 中国石油天然气集团公司 | 一种共转换点道集抽取方法及装置 |
CN104267432B (zh) * | 2014-10-16 | 2017-02-08 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于规则化的转换波共转换点道集高精度抽取方法 |
CN107422372A (zh) * | 2016-05-24 | 2017-12-01 | 中国石油化工股份有限公司 | 转换波面元映射快速ccp抽道集方法 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6208587B1 (en) * | 1997-07-10 | 2001-03-27 | Pgs Tensor, Inc. | Method of detecting seismic events and for detecting and correcting geometry and statics error in seismic data |
CN101251604A (zh) * | 2008-04-12 | 2008-08-27 | 中国石油集团西北地质研究所 | 二参数转换波速度分析及动校正方法 |
-
2010
- 2010-07-09 CN CN 201010229078 patent/CN101963674A/zh active Pending
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6208587B1 (en) * | 1997-07-10 | 2001-03-27 | Pgs Tensor, Inc. | Method of detecting seismic events and for detecting and correcting geometry and statics error in seismic data |
CN101251604A (zh) * | 2008-04-12 | 2008-08-27 | 中国石油集团西北地质研究所 | 二参数转换波速度分析及动校正方法 |
Non-Patent Citations (1)
Title |
---|
《CREWES Research Report》 19941231 Shaowu Wang et al. Fast 3-D converted-wave depth-variant common conversion point stacking 25-1~25-14 1-5 第6卷, * |
Cited By (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103645498A (zh) * | 2013-11-08 | 2014-03-19 | 中国石油大学(北京) | 一种直接生成共转换点道集的方法 |
CN103645498B (zh) * | 2013-11-08 | 2015-08-12 | 中国石油大学(北京) | 一种直接生成共转换点道集的方法 |
CN105891881A (zh) * | 2014-05-15 | 2016-08-24 | 蔡晓刚 | 一种水平界面vti介质转换点ccp道集抽取方法 |
CN104267432B (zh) * | 2014-10-16 | 2017-02-08 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于规则化的转换波共转换点道集高精度抽取方法 |
CN107422372A (zh) * | 2016-05-24 | 2017-12-01 | 中国石油化工股份有限公司 | 转换波面元映射快速ccp抽道集方法 |
CN106019373A (zh) * | 2016-06-24 | 2016-10-12 | 中国石油天然气集团公司 | 一种共转换点道集抽取方法及装置 |
CN106019373B (zh) * | 2016-06-24 | 2018-06-01 | 中国石油天然气集团公司 | 一种共转换点道集抽取方法及装置 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102967882B (zh) | 地层的层速度模型的建模方法 | |
CN101630016B (zh) | 一种提高垂直地震剖面成像质量的方法 | |
CN102841376A (zh) | 一种基于起伏地表的层析速度反演方法 | |
US7630865B2 (en) | Complex analysis of kinematics for non-hyperbolic moveout corrections | |
CN103149586B (zh) | 一种倾斜层状粘弹性介质中波场正演模拟方法 | |
CN104570102B (zh) | 一种近地表速度模型与中深层速度模型的融合方法 | |
CN101915938B (zh) | 一种转换波的偏移成像方法及装置 | |
CN105549081A (zh) | 各向异性介质共炮域高斯束偏移成像方法 | |
CN101598808B (zh) | 一种提高地震资料成像质量的方法 | |
CN102914796B (zh) | 一种基于高斯束的获取纵横波偏移速度的控制方法 | |
CN105589100A (zh) | 一种微地震震源位置和速度模型同时反演方法 | |
CN101963674A (zh) | 共转换点抽道集方法 | |
CN102866426B (zh) | 一种利用avo大角度道集分析岩体油气信息的方法 | |
CN104570073B (zh) | 一种适用于复杂高陡构造的双反射地震波成像方法 | |
CN104570124A (zh) | 一种适合井间地震大角度反射条件的延拓成像方法 | |
CN102053260B (zh) | 获得地震纵波的方位速度的方法及处理地震数据的方法 | |
CN102901984A (zh) | 真地表地震数据倾角道集构建方法 | |
CN106338766B (zh) | 基于分步傅里叶算法的叠前时间偏移方法 | |
CN107515420A (zh) | 一种用于局部相关同相轴的走时与梯度精确拾取方法 | |
CN106199704A (zh) | 一种三维三分量海底电缆地震资料速度建模方法 | |
CN103823241B (zh) | 一种在炮检距域中求取折射剩余静校正量的方法 | |
CN101609167A (zh) | 基于起伏地表的井间地震波动方程叠前深度偏移成像方法 | |
Gong et al. | Combined migration velocity model-building and its application in tunnel seismic prediction | |
Biloti et al. | Multiparametric traveltime inversion | |
CN108693560A (zh) | 一种基于互相关道的散射波成像方法及*** |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C12 | Rejection of a patent application after its publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20110202 |