CN101520516B - 一种三维地震记录的叠前似二维变换的方法 - Google Patents

一种三维地震记录的叠前似二维变换的方法 Download PDF

Info

Publication number
CN101520516B
CN101520516B CN2008100072964A CN200810007296A CN101520516B CN 101520516 B CN101520516 B CN 101520516B CN 2008100072964 A CN2008100072964 A CN 2008100072964A CN 200810007296 A CN200810007296 A CN 200810007296A CN 101520516 B CN101520516 B CN 101520516B
Authority
CN
China
Prior art keywords
geophone
dimensional
record
shot
transform
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.)
Active
Application number
CN2008100072964A
Other languages
English (en)
Other versions
CN101520516A (zh
Inventor
詹毅
赵波
周熙襄
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
BGP Inc
Original Assignee
BGP Inc
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by BGP Inc filed Critical BGP Inc
Priority to CN2008100072964A priority Critical patent/CN101520516B/zh
Priority to PCT/CN2009/000199 priority patent/WO2009105964A1/zh
Priority to US12/919,410 priority patent/US20110010098A1/en
Publication of CN101520516A publication Critical patent/CN101520516A/zh
Application granted granted Critical
Publication of CN101520516B publication Critical patent/CN101520516B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28—Processing seismic data, e.g. for interpretation or for event detection
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
    • G06T11/00—2D [Two Dimensional] image generation
    • G06T11/20—Drawing from basic elements, e.g. lines or circles
    • G06T11/206—Drawing of charts or graphs

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Remote Sensing (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geophysics (AREA)
  • Theoretical Computer Science (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明涉及地球物理勘探三维地震资料叠前处理，是三维地震记录的叠前似二维变换的方法。步骤是：采集三维地震数据按炮集排列，从第一炮的第一条测线开始，算出所有检波点的炮检距，以炮点为圆心，以各检波点的炮检距为半径画圆，将各检波点旋转到直线上，完成该炮所有测线的似二维变换得到高精度的三维地震数据图，使用常规的方法进行三维地震记录的去噪和静校正处理。本发明可将二维折射波静校正技术应用到三维资料处理，初至波和线性干扰都能很好地恢复了其线性特征，有助于后续处理的进行，可很好地压制线性干扰，去噪效果良好。

Description

一种三维地震记录的叠前似二维变换的方法
技术领域
本发明涉及地球物理勘探技术，属于三维地震资料叠前处理，是一种将成熟的二维地震资料叠前处理方法直接用于三维地震资料处理的三维地震记录的叠前似二维变换的方法。
背景技术
在地球物理勘探技术中，对采集而来的数据需要进行各种处理，以增强地震反演的准确性。目前，对于三维地震记录的处理，由于三维地震数据的特殊性，在二维记录中的线性特征呈非线性化，从而导致在数据处理时无法直接将成熟的二维叠前处理运用在三维数据上。为了在叠前压制线性干扰时保持线性干扰的“线性性”，只能先将三维记录抽成共偏移距道集，或者进行速度校正，然后再使用二维压制线性干扰的方法进行去噪处理，应用不便。另外，在二维折射波静校正方法中，常常采用线性拟合初至波的方法实现静校正量的求取，而在三维情况下，初至波除了炮点在某一条测线之上的记录呈线性状态外，其它测线的初至都呈非线性，而且随着炮点到测线的距离不同，初至曲线的曲率也在变化，这样就给初至拟合带来了困难，增加了静校正的难度，也影响了静校正处理的效率和精度。
发明内容
本发明提供一种将三维地震记录进行变换，在此基础上可以方便的进行三维地震记录的去噪和静校正处理，可以提高处理的效率和精度的三维地震记录的叠前似二维变换的方法。
本发明采用以下技术方案实现：
1)用常规的方法野外采集三维地震数据，进行解编后，按照炮集记录方式排列；
2)从第一炮的第一条测线开始，根据炮点(S)的坐标(X，Y)和所有检波点的坐标(Xi，Yi)计算出该炮所有检波点的炮检距(offseti)；
所述的炮检距计算是： offsetk = ( X k - X ) 2 + ( Y k - Y ) 2 ) ;
其中：k为第k个检波点的炮检距；X，Y为炮点的坐标，Xk，Yk为第k个检波点的坐标；
3)以穿过炮点(S)和最小炮检距的检波点(R1)的直线(L2’)作为变换后的坐标轴；
4)以炮点(S)为圆心，以各检波点的炮检距(offseti)为半径画圆，将各检波点旋转到直线(L2’)轴上；
所述的炮检距的检波点(R1’)与R重合，其余检波点按照R2’，…Rk’…Rn’形成一种新的排列方式，完成该条测线的似二维变换；
5)换该炮的另一条测线，重复上述2～4的步骤，完成该炮所有测线的似二维变换；
6)其余各炮记录重复2～5的步骤依次进行，完成地震资料叠前三维记录的似二维变换；
7)得到高精度的三维地震数据图，可以使用常规的方法进行三维地震记录的去噪和静校正处理。
本发明可将成熟的二维折射波静校正技术应用到三维资料处理之中，经过本发明的处理，初至波和线性干扰都能很好地恢复了其线性特征，有助于后续处理的进行，可很好地压制线性干扰，去噪效果良好。
附图说明
图1是本发明三维记录的似二维变换示意图；
图2是实际三维单炮记录中某三条测线按照常规排列方式示意图；
图3是图2的记录经过本发明似二维变换后示意图。
具体实施方案
本发明针对地震资料叠前三维记录，采用坐标变换和投影的方法将其变换为似二维记录，
本发明的实现过程如图1所示，图中S代表炮点，R1，R2，…Rk…Rn代表1～n个检波点位置，L2代表这n个检波点所在的测线坐标，offsetk和offsetn分别代表第k个和第n个检波点的真实炮检距；L2’是变换后新的排列方式坐标，R1’…Rk’…Rn’代表变换后n个检波点在新的坐标上的位置。
具体步骤如下：
1)用常规的方法野外采集三维地震数据，进行解编后，按照炮集记录方式排列；
2)从第一炮的第一条测线开始，根据炮点(S)的坐标(X，Y)和所有检波点的坐标(Xi，Yi)计算出该炮所有检波点的炮检距；第k个检波点的炮检距是 offsetk = ( X k - X ) 2 + ( Y k - Y ) 2 ) ;
3)以穿过炮点(S)和最小炮检距的检波点R1的直线L2’作为变换后的坐标轴(如图1所示)；
4)以炮点S为圆心，以各检波点的炮检距offseti为半径画圆，将各检波点旋转到L2’轴上，形成一种新的排列方式，完成该条测线的似二维变换；
5)换该炮的另一条测线，重复上述2～4的步骤，完成该炮所有测线的似二维变换；
6)其余各炮记录重复2～5的步骤依次进行，即可完成地震资料叠前三维记录的似二维变换；
7)得到高精度的三维地震数据图，可以使用常规的方法进行三维地震记录的去噪和静校正处理。
图2是本发明实际三维叠前单炮记录中抽取出来的某一炮的三条测线，按照常规排列方式的示意图，其中纵坐标是时间，横坐标是道数。从图中可以观察到，随着炮点与测线距离的不同，初至波和记录中的线性干扰呈现的状态也不相同，距离越小，则线性性越好，距离越大，非线性特征越强。在这种情况下，就压制线性干扰而言，成熟的二维去噪方法，例如倾角扫描压制线性干扰等方法很难取得好的去噪效果；就折射波静校正而言，由于各条测线初至波的曲率都不相同，极大地增加了初至波的拟合难度，也无法将成熟的二维折射波静校正技术应用到三维资料处理之中。
图3是图2的记录，按照上述方法进行似二维变换后的示意图，其中横坐标是炮检距，纵坐标为时间。对比图2和图3可以发现，经过变换后，记录中的初至波和线性干扰都很好地恢复了其线性特征，有助于后续处理的进行。

Claims (3)

1.一种三维地震记录的叠前似二维变换的方法，其特征在于采用以下步骤实现：
1)用常规的方法野外采集三维地震数据，进行解编后，按照炮集记录方式排列；
2)从第一炮的第一条测线开始，根据炮点(S)的坐标(X，Y)和所有检波点的坐标(Xi，Yi)计算出该炮所有检波点的炮检距(offseti)；
3)以穿过炮点(S)和最小炮检距的检波点(R1)的直线(L2’)作为变换后的坐标轴；
4)以炮点(S)为圆心，以各检波点的炮检距(offseti)为半径画圆，将各检波点旋转到变换后的直线(L2’)轴上；
5)换该炮的另一条测线，重复上述2)～4)的步骤，完成该炮所有测线的似二维变换；
6)其余各炮记录重复2)～5)的步骤依次进行，完成地震资料叠前三维记录的似二维变换；
7)得到高精度的三维地震数据图，使用常规的方法进行三维地震记录的去噪和静校正处理。
2.根据权利要求1所述的三维地震记录的叠前似二维变换的方法，其特征在于步骤2)所述的炮检距计算是：
Figure FSB00000627419700011
其中：offsetk为第k个检波点的炮检距；X，Y为炮点的坐标，Xk，Yk为第k个检波点的坐标。
3.根据权利要求1所述的三维地震记录的叠前似二维变换的方法，其特征在于所述的最小炮检距的检波点旋转变换后的检波点(R1’)与最小炮检距的检波点R1重合，其余旋转变换后的检波点按照R2’，…Rk’…Rn’排列，完成测线的似二维变换。
CN2008100072964A 2008-02-25 2008-02-25 一种三维地震记录的叠前似二维变换的方法 Active CN101520516B (zh)

Priority Applications (3)

Application Number Priority Date Filing Date Title
CN2008100072964A CN101520516B (zh) 2008-02-25 2008-02-25 一种三维地震记录的叠前似二维变换的方法
PCT/CN2009/000199 WO2009105964A1 (zh) 2008-02-25 2009-02-24 一种三维地震记录的叠前似二维变换的方法
US12/919,410 US20110010098A1 (en) 2008-02-25 2009-02-24 Method of pre-stack two-dimension-like transformation of three-dimensional seismic record

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2008100072964A CN101520516B (zh) 2008-02-25 2008-02-25 一种三维地震记录的叠前似二维变换的方法

Publications (2)

Publication Number Publication Date
CN101520516A CN101520516A (zh) 2009-09-02
CN101520516B true CN101520516B (zh) 2012-02-29

Family

ID=41015519

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2008100072964A Active CN101520516B (zh) 2008-02-25 2008-02-25 一种三维地震记录的叠前似二维变换的方法

Country Status (3)

Country Link
US (1) US20110010098A1 (zh)
CN (1) CN101520516B (zh)
WO (1) WO2009105964A1 (zh)

Families Citing this family (13)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102004264B (zh) * 2010-10-18 2015-09-23 中国石油化工股份有限公司 一种地震采集资料质量定量分析与评价方法
CN103076626B (zh) * 2011-10-25 2016-02-03 中国石油化工股份有限公司 一种波场净化处理方法
CN103837896B (zh) * 2014-03-10 2016-08-17 中国石油集团川庆钻探工程有限公司地球物理勘探公司 对地震资料进行叠前去噪的方法
CN104769458A (zh) * 2014-07-15 2015-07-08 杨顺伟 一种基于柯西分布的叠后波阻抗反演方法
CN104133245B (zh) * 2014-07-22 2018-01-02 中国石油天然气集团公司 一种地震资料的静校正方法及系统
CN104375181B (zh) * 2014-11-24 2017-03-08 中国石油天然气集团公司 一种快速生成三维共检波点初至道集的方法
CN105158795B (zh) * 2015-08-27 2018-03-02 中国石油集团川庆钻探工程有限公司地球物理勘探公司 利用地层叠前纹理属性值来检测缝洞的方法
CN105911585B (zh) * 2016-07-05 2018-05-15 中国石油集团东方地球物理勘探有限责任公司 一种地震记录规则干扰波的提取方法及装置
US10935682B2 (en) 2016-10-03 2021-03-02 Schlumberger Technology Corporation Downhole seismic sensing synchronization systems and methods
CN111965729B (zh) * 2019-05-20 2023-04-07 中国石油天然气集团有限公司 可控震源组合中心实时监控方法、系统及装置
CN111078953B (zh) * 2019-11-27 2023-11-07 滨州学院 基于导航文件直接重构segy数据坐标的方法
CN113945979B (zh) * 2020-07-17 2023-10-27 中国石油天然气股份有限公司 一种基于数据重构的三维锥形滤波方法及装置
US12072460B2 (en) 2021-12-15 2024-08-27 Saudi Arabian Oil Company System and method for determining a set of first breaks of a seismic dataset

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2301889A (en) * 1995-06-06 1996-12-18 Inst Francais Du Petrole Synthesising zero-offset seismic data
CN1369716A (zh) * 2002-03-13 2002-09-18 大庆油田有限责任公司 一种平点地震处理解释油气检测方法
CN1888934A (zh) * 2006-07-17 2007-01-03 石殿祥 一种三维地震偏移成像的块追赶方法
CN1896768A (zh) * 2005-07-14 2007-01-17 中国石油化工股份有限公司 一种用于显示三维地震数据反射结构的处理方法
CN101105536A (zh) * 2006-07-12 2008-01-16 中国石油集团东方地球物理勘探有限责任公司 一种复杂地表区的条带状地震采集方法

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4803669A (en) * 1986-12-19 1989-02-07 Atlantic Richfield Company Three-component three-dimensional seismic data acquisition
CA2176058C (en) * 1996-05-08 1999-03-02 William Nicholas Goodway Three-dimensional seismic acquisition
GB9813851D0 (en) * 1998-06-27 1998-08-26 Geco Prakla Uk Ltd Seismic data acquisition and processing method
CN1128373C (zh) * 2001-01-21 2003-11-19 周熙襄 用初至波对地震勘探资料进行短波长静校正处理的方法
US6691039B1 (en) * 2002-08-30 2004-02-10 John M. Robinson Removal of noise from seismic data using improved radon transformations
US6912468B2 (en) * 2003-08-14 2005-06-28 Westerngeco, L.L.C. Method and apparatus for contemporaneous utilization of a higher order probe in pre-stack and post-stack seismic domains
CN100349004C (zh) * 2004-12-29 2007-11-14 中国石油天然气集团公司 三维叠前地震数据中规则噪声的压制方法
EP1960234A1 (en) * 2005-12-16 2008-08-27 Mentor Graphics Corporation Two-dimensional orthogonal wire harness representation
CN100487489C (zh) * 2006-07-12 2009-05-13 中国石油集团东方地球物理勘探有限责任公司 一种高精度的深度域叠前地震数据反演方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
GB2301889A (en) * 1995-06-06 1996-12-18 Inst Francais Du Petrole Synthesising zero-offset seismic data
CN1369716A (zh) * 2002-03-13 2002-09-18 大庆油田有限责任公司 一种平点地震处理解释油气检测方法
CN1896768A (zh) * 2005-07-14 2007-01-17 中国石油化工股份有限公司 一种用于显示三维地震数据反射结构的处理方法
CN101105536A (zh) * 2006-07-12 2008-01-16 中国石油集团东方地球物理勘探有限责任公司 一种复杂地表区的条带状地震采集方法
CN1888934A (zh) * 2006-07-17 2007-01-03 石殿祥 一种三维地震偏移成像的块追赶方法

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
B.O.Ekren，B.Ursin，高顺莉.真振幅频率—波数恒定炮检距偏移.《国外油气勘探》.2000,第12卷(第04期),全文. *
JP特开9-90051A 1997.04.04
赵景霞, 张叔伦等.三维并行合成震源记录叠前深度偏移.《地球物理学报》.2006,第49卷(第1期),全文. *
钱忠平,赵波.三维地震资料叠前规则噪声压制.《CPS/SEG2004国际地球物理会议》.2004,全文. *

Also Published As

Publication number Publication date
CN101520516A (zh) 2009-09-02
US20110010098A1 (en) 2011-01-13
WO2009105964A1 (zh) 2009-09-03

Similar Documents

Publication Publication Date Title
CN101520516B (zh) 一种三维地震记录的叠前似二维变换的方法
CN102012521B (zh) 一种地震储层预测中叠前裂缝的检测方法
CN102213769A (zh) 一种利用三维垂直地震剖面资料确定各向异性参数的方法
CN104459794B (zh) 共反射点道集时变时间差值的校正方法及装置
CN102590862B (zh) 补偿吸收衰减的叠前时间偏移方法
CN101598811B (zh) 一种二维垂直地震剖面数据计算炮点静校正的方法
CN102073067B (zh) 一种提高地震数据初至波自动拾取效率的方法
CN104570125A (zh) 一种利用井数据提高成像速度模型精度的方法
CN102692648A (zh) 一种基于大炮初至的折射波剩余静校正方法
CN105259570A (zh) 一种海上时移地震接收点位移校正方法
CN104570122B (zh) 基于反射波的地震数据静校正方法及装置
CN103091709B (zh) 获得纵波、转换波地震数据时间匹配关系的方法和装置
CN109765616B (zh) 一种保幅波场延拓校正方法及系统
CN102636809B (zh) 一种传播角度域共成像点道集的生成方法
CN102116869A (zh) 高精度叠前域最小二乘偏移地震成像技术
CN107561589A (zh) 一种近地表横波层速度模型建立方法
CN101545986A (zh) 基于最大能量旅行时计算的三维积分叠前深度偏移方法
CN104570076A (zh) 一种基于二分法的地震波初至自动拾取方法
CN104570102A (zh) 一种近地表速度模型与中深层速度模型的融合方法
CN105093319A (zh) 基于三维地震数据的地面微地震静校正方法
CN103576200A (zh) 低信噪比区浅层波阻抗界面静校正方法
CN105204065A (zh) 一种拾取初至波的方法和装置
CN109239773A (zh) 一种高阶模式瑞雷波的重建方法
CN102053260B (zh) 获得地震纵波的方位速度的方法及处理地震数据的方法
CN102692651A (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
C14 Grant of patent or utility model
GR01 Patent grant
EE01 Entry into force of recordation of patent licensing contract

Application publication date: 20090902

Assignee: NATIONAL ENGINEERING RESEARCH CENTER OF OIL AND GAS EXPLORATION SOFTWARE CO., LTD.

Assignor: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp.

Contract record no.: 2013990000780

Denomination of invention: Method for pre-stack like-two-dimensional transform of three-dimensional seismic record

Granted publication date: 20120229

License type: Exclusive License

Record date: 20131121

LICC Enforcement, change and cancellation of record of contracts on the licence for exploitation of a patent or utility model