CN103645498B - 一种直接生成共转换点道集的方法 - Google Patents
一种直接生成共转换点道集的方法 Download PDFInfo
- Publication number
- CN103645498B CN103645498B CN201310553398.7A CN201310553398A CN103645498B CN 103645498 B CN103645498 B CN 103645498B CN 201310553398 A CN201310553398 A CN 201310553398A CN 103645498 B CN103645498 B CN 103645498B
- Authority
- CN
- China
- Prior art keywords
- transfer point
- common
- point
- window
- delta
- 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.)
- Expired - Fee Related
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000006243 chemical reaction Methods 0.000 title claims abstract description 32
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 5
- 101100409194 Rattus norvegicus Ppargc1b gene Proteins 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 claims description 3
- 238000000926 separation method Methods 0.000 claims description 3
- 230000003068 static effect Effects 0.000 claims description 3
- 230000008030 elimination Effects 0.000 claims description 2
- 238000003379 elimination reaction Methods 0.000 claims description 2
- 230000006978 adaptation Effects 0.000 abstract description 3
- 238000003384 imaging method Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 3
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 239000012141 concentrate Substances 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明提供一种直接生成共转换点道集的方法,假设转换波地震数据已预处理好,已获得准确的动校正参数和计算转换点位置的速度等参数,按炮检距分组,分时窗叠加获得时间和空间变化的共转换点道集。本发明实现简便、适应大数据量处理。
Description
技术领域
本发明涉及转换波地震资料处理技术领域,具体地,涉及一种直接时空变生成共转换点道集的方法。
背景技术
在多分量转换波地震资料处理中,共转换点道集的抽取方法一直是一项关键技术。因为转换波的下行波是纵波、上行波是横波,造成传播路径的不对称性,转换点位于从中心点偏向接收点的一侧,且随深度及界面上下纵横波速度的变化而变化。这使得基于常规共中心点(CMP)道集的速度分析、动校正和水平叠加技术不能适用于转换波,转换波叠加成像必须按共转换点进行道集分选。目前,有关转换点位置的计算已有很多方法,主要包括渐近线法(ACP)、单层各向同性介质迭代解、多层介质指定深度各向同性近似解(Thomsen,1999)、多层介质指定旅行时各向同性近似解(Thomsen,1999)和多层VTI介质指定旅行时近似解(Yuan&Li,2001)等。另外,周竹生等(1993)导出了一种快速、高精度转换点的迭代算法;许士勇等(2002)给出了四次方程的近似表达式;姚陈等(2005)给出了三维倾角转换波共转换点的解析表达式;苑春方等(2005,2006,2006)推导出了单层各向同性介质与深度有关的唯一的解析解和含有旅行时间唯一解析解,倾斜界面上的转换点的精确解;Chang等(2005)利用倾斜界面上转换波射线路径的三角几何关系,推导出倾斜界面上的转换点的近似公式;毕丽飞等(2005)用水平界面推导转换点坐标迭代解类似的方法。但目前,共转换点道集的选排方法仍然是以整道选排为主,主要是渐近逼近道分类法、依赖深度的选排方法。渐近逼近道分类法是通过一个工区的固定速度比,计算转换点的渐近线,然后将整个工区的转换波数据抽取到渐近线道集中,即ACP道集,这种方法能一定程度上能够满足深部地层的成像要求,但是浅层的成像效果较差。依赖深度的选排方法只对某一深度最为正确,都是能在目标层聚焦成像,更适合目的层较浅的转换波处理,更适合更新速度分析,并且指定目的层计算提供了多种转换点位置的计算方法,但在所选反射层的上面或下面,反射界面的叠后成像可能会变得模糊。另外,分层选排方法能够使不同深度的转换波均能达到较好的聚焦效果,但需要在不同层之间进行道集的拼接,特别是在大数据量的情况下实际操作困难。
上述的共转换点道集的选排方法都有一定的应用条件,并不能很好的满足转换波对转换点在时间上和空间上变化的要求。
发明内容
本发明目的是提供一种能适应时间上和空间上变化的共转换点道集的选排方法,一种直接时空变生成共转换点道集的方法。
本发明通过以下技术方案实现:
一种直接生成共转换点道集的方法,具体步骤包括:
1)采集原始转换波地震数据,经预处理和速度分析,得到转换波动校正和共转换点位置计算需要的参数,并对预处理后转换波地震数据进行动校正,且通过切除消除动校拉伸的影响,
步骤1)所述的预处理包括加载观测系统、静校正、波场分离、叠前去噪、振幅恢复、预置共转换点CCP号;
2)步骤1)处理后的转换波地震数据按炮检距大小进行分组,每一组对应一个炮检距范围;`
3)输入一组炮检距范围内的数据;
4)按道头中的共转换点CCP号读取转换波用来计算转换点位置的参数;
5)对地震道进行分时窗处理,
步骤5)所述的分时窗处理是指将地震道分成多个时窗,时窗长度不变,时窗重叠,重叠时窗的重叠长度为时窗内样点数的百分比perc,perc取值为0-100;
6)从第一个时窗开始,计算时窗中点的共转换点CCP号,并由这个共转换点号重新获得该时窗用来计算转换点位置的参数,之后重新计算时窗中点相对准确的共转换点CCP号,直接把时窗内的样点值与窗函数相乘叠在共转换点地震道上,即转换波共反射面元的宽度上,相应共转换点地震道每个样点的覆盖次数加上窗函数的数值,
步骤6)所述的计算转换点位置的参数是预先速度分析得到的,是转换波均方根速度纵波均方根速度各向异性参数χeff、纵横波垂直速度比γ0和纵横波有效速度比reff,
步骤6)所述的共转换点CCP号,是指按工区面元划分的,与常规处理中纵波共中心点CMP号一致,
步骤6)所述的共转换点CCP号是由下述方法计算的:
xC=xM+D (1)
式中xM、xC为共中心点CMP和共转换点CCP点的水平坐标,D为转换点偏离中心点的水平距离,
式中xp为共转换点与炮点之间的水平距离,
已知共中心点CMP号计算相应的共转换点CCP号
其中:INTn表示取最接近的整数,取整隐含做面元化处理,Δxm为工区面元大小,
所述的共转换点与炮点之间的水平距离,采用单层各向同性介质迭代解、解析解、多层介质指定旅行时各向同性近似解、多层VTI介质指定旅行时近似解、倾斜界面迭代解、解析解计算得到,
步骤6)所述的直接把时窗内的样点值与窗函数相乘叠在共转换点地震道上,即转换波共反射面元的宽度上,是把这个时窗内的样点叠加在准确的共转换点CCP号上以及相邻的I个共转换点CCP号上,I=0,1,2,…,INT(2·m/Δx+0.5),m为转换波共反射面元的宽度,INT表示取整,
步骤6)所述的转换波共反射面元的宽度采用以下关系式确定:
其中:Δx为道间距,Δs为炮间距,r为纵横波速度比,INT表示取整;
7)重复步骤6)直到该道的每个时窗完成;
8)重复步骤4)至7)直到炮检距分组内所有道输入完毕,输出该组的共转换点地震道,
步骤8)所述的输出该组的共转换点地震道为共转换点地震道/共转换点地震道的覆盖次数,并在道头中置上该组炮检距的大小;
9)重复步骤3)至8)直到所有的炮检距分组输入完毕,就得到了时间和空间变化的共转换点叠前道集。
本发明的有益效果为:假设转换波地震数据已预处理好,已获得准确的动校正参数和计算转换点位置的速度等参数,按炮检距分组,分时窗叠加获得时间和空间变化的共转换点道集。本发明实现简便、适应大数据量处理。
附图说明
图1是应用该方法生成的共转换点道集图;
具体实施方式
本发明在转换波数据上采用以下步骤应用:
1)采集原始转换波地震数据,经预处理和速度分析,得到转换波各向异性动校正和共转换点位置各向异性计算需要的速度参数:转换波均方根速度纵横波垂直速度比γ0、纵横波有效速度比reff和各向异性参数χeff,并对预处理后转换波地震数据进行各向异性动校正,且通过切除消除动校拉伸的影响;
步骤1)所述的预处理包括加载观测系统,静校正,波场分离,叠前去噪,振幅恢复,预置共转换点(CCP)号。
2)步骤1)处理后的转换波地震数据按炮检距最小为-4000m,最大为4000m,增量为50m进行分组;`
3)输入一组炮检距范围内的数据;
4)按道头中的共转换点(CCP)号读取转换波用来计算转换点位置的速度等参数;
5)对地震道进行分时窗处理;
步骤5)所述的分时窗处理是指将地震道分成多个时窗,时窗长度为100ms,时窗重叠50%。
6)从第一个时窗开始,计算时窗中点的共转换点(CCP)号,并由这个共转换点号重新获得该时窗用来计算转换点位置的速度等参数,之后重新计算时窗中点相对准确的共转换点(CCP)号,直接把时窗内的样点值与窗函数相乘叠在共转换点地震道上,即转换波共反射面元的宽度上,相应共转换点地震道每个样点的覆盖次数加上窗函数的数值;
步骤6)所述的计算转换点位置的速度等参数是预先速度分析得到的转换波均方根速度纵横波垂直速度比γ0、纵横波有效速度比reff和各向异性参数χeff;
步骤6)所述的共转换点(CCP)号,是指按工区面元划分的,与常规处理中纵波共中心点(CMP)号一致。
步骤6)所述的共转换点(CCP)号是由下述方法计算的:
xC=xM+D(1)
式中xM、xC为共中心点(CMP)和共转换点(CCP)点的水平坐标,D为转换点偏离中心点的水平距离
式中xp共转换点与炮点之间的水平距离,
已知共中心点(CMP)号计算相应的共转换点(CCP)号
其中:INTn表示取最接近的整数,取整隐含做面元化处理,Δxm为工区面元大小,大小为25m。
步骤6)所述的共转换点与炮点之间的水平距离,可采用多层VTI介质指定旅行时近似解计算得到(下式)。
式中, x是炮检距,r0是垂直速度比,reff是有效速度比,t0ps是转换波垂直双程旅行时,vc是转换波动校正速度,ηeff是纵波各向异性参数,χeff是转换波各向异性参数。
步骤6)所述的转换波共反射面元的宽度采用以下关系式确定:
其中:Δx为道间距,Δs为炮间距,r为纵横波速度比,INT表示取整。
7)重复步骤6)直到该道的每个时窗完成;
8)重复步骤4)至7)直到炮检距分组内所有道输入完毕,输出该组的共转换点地震道;
步骤8)所述的输出该组的共转换点地震道为共转换点地震道/共转换点地震道的覆盖次数,并在道头中置上该组炮检距的大小。
9)重复步骤3)至8)直到所有的炮检距分组输入完毕,就得到了时间和空间变化的共转换点叠前道集。
图1是直接时空变生成的共转换点道集,效果较好,反射波同相轴能量突出。
Claims (1)
1.一种直接生成共转换点道集的方法,具体步骤包括:
1)采集原始转换波地震数据,经预处理和速度分析,得到转换波动校正和共转换点位置计算需要的参数,并对预处理后转换波地震数据进行动校正,且通过切除消除动校拉伸的影响,
步骤1)所述的预处理包括加载观测系统、静校正、波场分离、叠前去噪、振幅恢复、预置共转换点CCP号;
2)步骤1)处理后的转换波地震数据按炮检距大小进行分组,每一组对应一个炮检距范围;`
3)输入一组炮检距范围内的数据;
4)按道头中的共转换点CCP号读取转换波用来计算转换点位置的参数;
5)对地震道进行分时窗处理,
步骤5)所述的分时窗处理是指将地震道分成多个时窗,时窗长度不变,时窗重叠,重叠时窗的重叠长度为时窗内样点数的百分比perc,perc取值为0-100;
6)从第一个时窗开始,计算时窗中点的共转换点CCP号,并由这个共转换点号重新获得该时窗用来计算转换点位置的参数,之后重新计算时窗中点相对准确的共转换点CCP号,直接把时窗内的样点值与窗函数相乘叠在共转换点地震道上,即转换波共反射面元的宽度上,相应共转换点地震道每个样点的覆盖次数加上窗函数的数值,
步骤6)所述的计算转换点位置的参数是预先速度分析得到的,是转换波均方根速度纵波均方根速度各向异性参数χeff、纵横波垂直速度比γ0和纵横波有效速度比reff,
步骤6)所述的共转换点CCP号,是指按工区面元划分的,与常规处理中纵波共中心点CMP号一致,
步骤6)所述的共转换点CCP号是由下述方法计算的:
xC=xM+D
式中xM、xC为共中心点CMP和共转换点CCP点的水平坐标,D为转换点偏离中心点的水平距离,
式中xp为共转换点与炮点之间的水平距离,
已知共中心点CMP号计算相应的共转换点CCP号
其中:INTn表示取最接近的整数,取整隐含做面元化处理,Δxm为工区面元大小,
所述的共转换点与炮点之间的水平距离,采用单层各向同性介质迭代解、解析解、多层介质指定旅行时各向同性近似解、多层VTI介质指定旅行时近似解、倾斜界面迭代解、解析解计算得到,
步骤6)所述的直接把时窗内的样点值与窗函数相乘叠在共转换点地震道上,即转换波共反射面元的宽度上,是把这个时窗内的样点叠加在准确的共转换点CCP号上以及相邻的I个共转换点CCP号上,I=0,1,2,…,INT(2·m/Δx+0.5),m为转换波共反射面元的宽度,INT表示取整,
步骤6)所述的转换波共反射面元的宽度采用以下关系式确定:
其中:Δx为道间距,Δs为炮间距,r为纵横波速度比,INT表示取整;
7)重复步骤6)直到该道的每个时窗完成;
8)重复步骤4)至7)直到炮检距分组内所有道输入完毕,输出该组的共转换点地震道,
步骤8)所述的输出该组的共转换点地震道为共转换点地震道/共转换点地震道的覆盖次数,并在道头中置上该组炮检距的大小;
9)重复步骤3)至8)直到所有的炮检距分组输入完毕,就得到了时间和空间变化的共转换点叠前道集。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310553398.7A CN103645498B (zh) | 2013-11-08 | 2013-11-08 | 一种直接生成共转换点道集的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310553398.7A CN103645498B (zh) | 2013-11-08 | 2013-11-08 | 一种直接生成共转换点道集的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103645498A CN103645498A (zh) | 2014-03-19 |
CN103645498B true CN103645498B (zh) | 2015-08-12 |
Family
ID=50250746
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310553398.7A Expired - Fee Related CN103645498B (zh) | 2013-11-08 | 2013-11-08 | 一种直接生成共转换点道集的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103645498B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN104267432B (zh) * | 2014-10-16 | 2017-02-08 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种基于规则化的转换波共转换点道集高精度抽取方法 |
CN106019373B (zh) * | 2016-06-24 | 2018-06-01 | 中国石油天然气集团公司 | 一种共转换点道集抽取方法及装置 |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101598803A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种直接得到转换波叠加剖面的方法 |
CN101923177A (zh) * | 2010-07-09 | 2010-12-22 | 中国石油天然气股份有限公司 | 多波地震资料的速度分析方法和装置 |
CN101963674A (zh) * | 2010-07-09 | 2011-02-02 | 中国科学院地质与地球物理研究所 | 共转换点抽道集方法 |
-
2013
- 2013-11-08 CN CN201310553398.7A patent/CN103645498B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101598803A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种直接得到转换波叠加剖面的方法 |
CN101923177A (zh) * | 2010-07-09 | 2010-12-22 | 中国石油天然气股份有限公司 | 多波地震资料的速度分析方法和装置 |
CN101963674A (zh) * | 2010-07-09 | 2011-02-02 | 中国科学院地质与地球物理研究所 | 共转换点抽道集方法 |
Non-Patent Citations (3)
Title |
---|
倾斜CCP道集抽取方法及影响因素分析;于光明等;《煤田地质与勘探》;20070815(第04期);全文 * |
共转换点道集的迭代算法及应用;崔秀等;《世界地质》;19981226(第04期);全文 * |
快速有效的转换波共转换点叠加技术;许士勇等;《地球物理学报》;20020717(第04期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN103645498A (zh) | 2014-03-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102998704B (zh) | 一种地球物理勘探地震资料处理方法 | |
CN102053261B (zh) | 一种地震数据处理方法 | |
CN101598803B (zh) | 一种直接得到转换波叠加剖面的方法 | |
CN104730579B (zh) | 一种基于表层横波速度反演的纵横波联合静校正方法 | |
CN102221709B (zh) | 基于地层参数信息的速度分析与动校正方法 | |
CN102721977B (zh) | 高精度转换波各向异性叠加速度分析方法 | |
CN103576197B (zh) | 一种转换波角道集抽取方法 | |
CN101738636B (zh) | 一种三维vsp高斯束法多波联合偏移成像方法 | |
CN104142518A (zh) | 一种地震观测系统叠前时间偏移响应的分析方法 | |
CN102866426B (zh) | 一种利用avo大角度道集分析岩体油气信息的方法 | |
CN104237937A (zh) | 叠前地震反演方法及其系统 | |
CN102540252A (zh) | 基于互相关的高精度中值叠加方法 | |
CN104316965A (zh) | 一种裂缝方位和强度的预测方法及系统 | |
CN102053260B (zh) | 获得地震纵波的方位速度的方法及处理地震数据的方法 | |
CN102053262B (zh) | 获得地震转换波的方位速度的方法及处理地震数据的方法 | |
CN102901984A (zh) | 真地表地震数据倾角道集构建方法 | |
CN103645498B (zh) | 一种直接生成共转换点道集的方法 | |
CN102798888B (zh) | 一种利用非零井源距数据计算纵横波速度比的方法 | |
CN105445788B (zh) | 一种基于模型和全局寻优的速度谱自动解释方法 | |
CN103217707B (zh) | 一种直接抽取纵波时间域转换波角道集的方法 | |
CN103472482A (zh) | 基于基因排序体系的多波地震资料时间域匹配方法及系统 | |
CN106990434A (zh) | 椭圆展开转换波成像方法及系统 | |
CN104678439A (zh) | 一种多波连续介质速度分析方法 | |
CN103076629B (zh) | 基于纵横地震波速度比进行油气勘探的方法及设备 | |
CN106226814B (zh) | 利用转换横波地震数据反演储层横波阻抗和密度的方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20150812 |
|
CF01 | Termination of patent right due to non-payment of annual fee |