CN107942385A - 一种精细旅行时计算的动校正方法 - Google Patents

一种精细旅行时计算的动校正方法 Download PDF

Info

Publication number
CN107942385A
CN107942385A CN201711130848.6A CN201711130848A CN107942385A CN 107942385 A CN107942385 A CN 107942385A CN 201711130848 A CN201711130848 A CN 201711130848A CN 107942385 A CN107942385 A CN 107942385A
Authority
CN
China
Prior art keywords
msubsup
root mean
bearing calibration
shot point
geophone station
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.)
Granted
Application number
CN201711130848.6A
Other languages
English (en)
Other versions
CN107942385B (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.)
China National Petroleum Corp
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 CN201711130848.6A priority Critical patent/CN107942385B/zh
Publication of CN107942385A publication Critical patent/CN107942385A/zh
Application granted granted Critical
Publication of CN107942385B publication Critical patent/CN107942385B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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
    • G01V1/30—Analysis
    • G01V1/308—Time lapse or 4D effects, e.g. production related effects to the formation
    • 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
    • G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362—Effecting static or dynamic corrections; Stacking
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/50—Corrections or adjustments related to wave propagation
    • G01V2210/52—Move-out correction
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/60—Analysis
    • G01V2210/61—Analysis by combining or comparing a seismic data set with other data
    • G01V2210/612—Previously recorded data, e.g. time-lapse or 4D
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00—Details of seismic processing or analysis
    • G01V2210/70—Other details related to processing
    • G01V2210/74—Visualisation of seismic data

Landscapes

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

Abstract

本发明提供了一种精细旅行时计算的动校正方法。所述动校正方法包括：A、对于地震数据中任一采样点，获取炮点均方根速度vS和检波点均方根速度vR；B、根据式1计算动校正的旅行时t，并进行动校正，式1为：其中，t0为炮点到检波点的自激自收时间，xS为炮点的偏移距，xR为检波点的偏移距；C、重复所述步骤A和B，直到地震数据中所有采样点处理完毕。本发明适用于地下情况复杂的地震数据处理，能够获得更准确均方根速度和动校正量，实现同相叠加，最终提升成像质量。

Description

一种精细旅行时计算的动校正方法
技术领域
本发明属于石油天然气地震勘探调查领域，更具体地讲，涉及一种石油地震勘探的地震资料处理解释方法。
背景技术
在石油地震勘探常规资料处理中，通常假设地下介质是水平层状，层内速度横向不变或变化比较缓慢。但在实际复杂地区生产应用中，地下介质速度的不均匀，纵、横向变化会引起地震波传播旅行时的不均匀，会导致动校正量计算不准确，使得同一共中心点道集无法同相叠加，最终影响成像质量。
发明内容
针对现有技术中存在的问题，本发明的目的在于解决现有技术存在的上述不足中的至少一项。例如，本发明的目的在于提供一种能够获得更准确的动校正量、实现同相叠加的动校正方法。
为了实现上述目的，本发明提供了一种精细旅行时计算的动校正方法。所述动校正方法包括：A、对于地震数据中任一采样点，获取炮点均方根速度vS和检波点均方根速度vR；B、根据式1计算动校正的旅行时t，并进行动校正；C、重复所述步骤A和B，直到地震数据中所有采样点处理完毕。其中，式1为：
在式1中，t0为炮点到检波点的自激自收时间，xS为炮点的偏移距，xR为检波点的偏移距。
在本发明的一个示例性实施例中，所述步骤A可以包括：获取炮点到反射样点路径上的所有均方根速度场并进行中值、均值和均方根计算中的一种或多种，得到所述炮点均方根速度vS；获取反射样点到接收点路径上的所有均方根速度场并进行中值、均值和均方根计算中的一种或多种，得到所述检波点均方根速度vR。
在本发明的一个示例性实施例中，所述动校正方法可以用于复杂地区的地震数据处理。
与现有技术相比，本发明的有益效果包括：适用于地下情况复杂的地震数据处理，能够获得更准确的动校正量，实现同相叠加，最终提升成像质量。
附图说明
图1示出了根据本发明的示例性实施例的精细旅行时计算的动校正方法的流程图。
具体实施方式
在下文中，将结合示例性实施例和附图来详细说明本发明的精细旅行时计算的动校正方法。
在石油地震勘探常规资料处理中，通常假设地下介质是水平层状，层内速度横向不变或变化比较缓慢，所以在同一共中心点道集内，其动校正使用的旅行时计算公式为公式1：
其中，t为偏移距x对应的旅行时，t0为自激自收时间，x为偏移距，v为均方根速度。
在实际复杂地区生产应用中，地下介质速度的不均匀，纵、横向变化会引起地震波传播旅行时的不均匀，即不严格满足公式1，继续使用固定不变的均方根速度v，会导致动校正量计算不准确，使得同一共中心点道集无法同相叠加，最终影响成像质量。
在复杂地区，由于速度场横向变化较大，同一共中心点道集的同一深度的动校正使用相同的均方根速度明显不合理。本发明提出了基于精细旅行时计算的复杂地区动校正方法，首先通过分别计算炮点和检波点旅行时的思路，为后续处理奠定基础，然后通过分别获得炮点、检波点的更精细的均方根速度场，并在动校正计算中应用，从而达到获得更准确的动校正量，实现同相叠加，最终提升成像质量。
本发明的技术思路是：(1)分别计算炮点和检波点旅行时：通过双平方根的计算公式(公式2)，来获得动校正的旅行时，以利于后面的叠加与成像不同需求的处理。(2)获得炮点、检波点的更精细的均方根速度场：如动校正计算中，只使用反射点位置的均方根速度，如果速度横向变化剧烈，且偏移距太大，这个速度明显无法近似炮点到反射点位置，反射点到接收点位置的均方根速度场。本发明通过获取炮点到反射点路径上的均方根速度场，并将这些速度进行中值、均值和均方根计算中的一种或多种计算(即任意一种、两种或三种计算均可)，获得新的炮点均方根速度vS，同理，通过获取反射点到接收点路径上的均方根速度场，并将这些速度进行中值、均值和均方根计算中的一种或多种计算，获得新的检波点均方根速度vR。将新的炮点、检波点均方根速度代入炮点和检波点旅行时计算公式(公式2)，获得更准确的动校正量。
根据本发明示例性实施例的精细旅行时计算的动校正方法包括以下步骤：
A、对于地震数据中任意一个采样点，获取炮点均方根速度vS和检波点均方根速度vR。
在本实施例中，获取炮点均方根速度vS和检波点均方根速度的方式优选为：获取炮点到反射样点路径上的所有均方根速度场并进行中值、均值和均方根计算一种或多种计算，得到炮点均方根速度vS；获取反射样点到接收点路径上的所有均方根速度场并进行中值、均值和均方根计算中的一种或多种计算，得到检波点均方根速度vR。其中，中值是取多个数据中，由小到大排序的中间的那个数；均值是多个数据中进行算术平均得到的那个数。
B、根据公式2计算动校正的旅行时t，并进行动校正。
在上式中，t为偏移距xS、xR对应的旅行时，t0为炮点到检波点的自激自收时间，xS和xR分别为炮点和检波点的偏移距，vS和vR分别为炮点和检波点的均方根速度。
C、重复所述步骤A和B，直到地震数据中所有采样点处理完毕。
下面将结合具体示例来对本发明的精细旅行时计算的动校正方法作进一步的说明。
示例
图1示出了根据本发明的示例性实施例的精细旅行时计算的动校正方法的流程图。如图1所示，动校正方法包括以下步骤：
(1)以输入均方根速度场和待处理的地震数据。
(2)根据定义的CDP进行依次处理。
(3)在同一CDP内，根据定义的采样依次处理。
(4)对于每一个地震数据样点，获取炮点到反射样点路径上的所有均方根速度场，并将这些速度进行中值、均值、均方根计算，获得新的炮点均方根速度，同理，通过获取反射样点到接收点路径上的所有均方根速度场，并将这些速度进行中值、均值和均方根计算，获得新的检波点均方根速度。
(5)将新的炮点、检波点均方根速度代入炮点和检波点旅行时计算公式2，获得更准确的动校正量，并进行动校正。
(6)判断当前CDP内的采样是否处理完毕，若“否”，则返回步骤(3)继续，若“是”，进入下一步.
(7)判断所有CDP是否处理完毕，若“否”，则返回步骤(2)继续，若“是”，进入下一步。
(8)输出动校正处理后的地震数据，计算结束。
综上所述，本发明通过获取炮点到反射点路径上的均方根速度场，并将这些速度进行中值、均值和均方根等计算，获得新的炮点均方根速度，同理，通过获取反射点到接收点路径上的均方根速度场，并将这些速度进行中值、均值和均方根等计算，获得新的检波点均方根速度。将新的炮点、检波点均方根速度代入炮点和检波点旅行时计算公式，获得更准确的动校正量。本发明基于精细旅行时计算的复杂地区动校正方法特别适用于地下情况复杂的地震数据处理，具有计算步骤简便和处理效果好等特点，在目前地震资料高精度处理中具有广阔的应用前景。
尽管上面已经结合附图和示例性实施例描述了本发明，但是本领域普通技术人员应该清楚，在不脱离权利要求的精神和范围的情况下，可以对上述实施例进行各种修改。

Claims (3)

1.一种精细旅行时计算的动校正方法，其特征在于，所述动校正方法包括：
A、对于地震数据中任一采样点，获取炮点均方根速度vS和检波点均方根速度vR；
B、根据式1计算动校正的旅行时t，并进行动校正，
式1为：
<mrow> <mi>t</mi> <mo>=</mo> <msqrt> <mrow> <msubsup> <mi>t</mi> <mn>0</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>x</mi> <mi>s</mi> <mn>2</mn> </msubsup> <mo>/</mo> <msubsup> <mi>v</mi> <mi>s</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> <mo>+</mo> <msqrt> <mrow> <msubsup> <mi>t</mi> <mn>0</mn> <mn>2</mn> </msubsup> <mo>+</mo> <msubsup> <mi>x</mi> <mi>R</mi> <mn>2</mn> </msubsup> <mo>/</mo> <msubsup> <mi>v</mi> <mi>R</mi> <mn>2</mn> </msubsup> </mrow> </msqrt> </mrow>
在式1中，t0为炮点到检波点的自激自收时间，xS为炮点的偏移距，xR为检波点的偏移距；
C、重复所述步骤A和B，直到地震数据中所有采样点处理完毕。
2.根据权利要求1所述的精细旅行时计算的动校正方法，其特征在于，所述步骤A包括：
获取炮点到反射样点路径上的所有均方根速度场并进行中值、均值和均方根计算中的一种或多种，得到所述炮点均方根速度vS；
获取反射样点到接收点路径上的所有均方根速度场并进行中值、均值和均方根计算中的一种或多种，得到所述检波点均方根速度vR。
3.根据权利要求1所述的精细旅行时计算的动校正方法，其特征在于，所述动校正方法用于复杂地区的地震数据处理。
CN201711130848.6A 2017-11-15 2017-11-15 一种精细旅行时计算的动校正方法 Active CN107942385B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201711130848.6A CN107942385B (zh) 2017-11-15 2017-11-15 一种精细旅行时计算的动校正方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201711130848.6A CN107942385B (zh) 2017-11-15 2017-11-15 一种精细旅行时计算的动校正方法

Publications (2)

Publication Number Publication Date
CN107942385A true CN107942385A (zh) 2018-04-20
CN107942385B CN107942385B (zh) 2020-02-07

Family

ID=61931279

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201711130848.6A Active CN107942385B (zh) 2017-11-15 2017-11-15 一种精细旅行时计算的动校正方法

Country Status (1)

Country Link
CN (1) CN107942385B (zh)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111856581A (zh) * 2020-07-27 2020-10-30 广州海洋地质调查局 一种obs时钟漂移校正方法及处理终端

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101984366A (zh) * 2010-09-29 2011-03-09 北京吉星吉达科技有限公司 基于共孔径面的起伏地表叠前时间偏移方法及设备
CN102162858A (zh) * 2010-12-06 2011-08-24 中国海洋石油总公司 利用非对称走时进行动校正速度分析的方法
CN102353992A (zh) * 2011-06-02 2012-02-15 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于中心点规则化的起伏地表动校正方法
CN102565857A (zh) * 2011-12-16 2012-07-11 中国石油集团川庆钻探工程有限公司地球物理勘探公司 自动剩余动校正方法
CN104122582A (zh) * 2014-06-24 2014-10-29 夏正元 利用叠加速度求取高精度地震波速度的方法
CN104216016A (zh) * 2014-08-12 2014-12-17 中国石油天然气集团公司 一种纵波约束扫描的转换波动校正方法以及系统
WO2016011670A1 (zh) * 2014-07-25 2016-01-28 杨顺伟 一种基于旅行时法确定地层裂隙裂缝走向的方法及装置
CN105866833A (zh) * 2016-06-16 2016-08-17 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Vsp-cdp叠加方法及三维vsp覆盖次数计算方法
CN106324673A (zh) * 2015-07-03 2017-01-11 中国石油化工股份有限公司 一种用于进行全方位纵波速度分析的方法和装置
CN107179548A (zh) * 2017-07-10 2017-09-19 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种基于真地表的叠前地震成像方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101984366A (zh) * 2010-09-29 2011-03-09 北京吉星吉达科技有限公司 基于共孔径面的起伏地表叠前时间偏移方法及设备
CN102162858A (zh) * 2010-12-06 2011-08-24 中国海洋石油总公司 利用非对称走时进行动校正速度分析的方法
CN102353992A (zh) * 2011-06-02 2012-02-15 中国石油集团川庆钻探工程有限公司地球物理勘探公司 基于中心点规则化的起伏地表动校正方法
CN102565857A (zh) * 2011-12-16 2012-07-11 中国石油集团川庆钻探工程有限公司地球物理勘探公司 自动剩余动校正方法
CN104122582A (zh) * 2014-06-24 2014-10-29 夏正元 利用叠加速度求取高精度地震波速度的方法
WO2016011670A1 (zh) * 2014-07-25 2016-01-28 杨顺伟 一种基于旅行时法确定地层裂隙裂缝走向的方法及装置
CN104216016A (zh) * 2014-08-12 2014-12-17 中国石油天然气集团公司 一种纵波约束扫描的转换波动校正方法以及系统
CN106324673A (zh) * 2015-07-03 2017-01-11 中国石油化工股份有限公司 一种用于进行全方位纵波速度分析的方法和装置
CN105866833A (zh) * 2016-06-16 2016-08-17 中国石油集团川庆钻探工程有限公司地球物理勘探公司 Vsp-cdp叠加方法及三维vsp覆盖次数计算方法
CN107179548A (zh) * 2017-07-10 2017-09-19 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种基于真地表的叠前地震成像方法

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN111856581A (zh) * 2020-07-27 2020-10-30 广州海洋地质调查局 一种obs时钟漂移校正方法及处理终端
CN111856581B (zh) * 2020-07-27 2022-02-22 广州海洋地质调查局 一种obs时钟漂移校正方法及处理终端

Also Published As

Publication number Publication date
CN107942385B (zh) 2020-02-07

Similar Documents

Publication Publication Date Title
CN104091064B (zh) 基于优化解空间搜索法的PS‑DInSAR地表形变测量参数估计方法
WO2015196779A1 (zh) 利用叠加速度求取高精度地震波速度的方法
CN107957594A (zh) 地震数据的椭圆校正方法、动校正方法及速度谱计算方法
CN104422963B (zh) 一种变偏移距vsp资料时差校正方法
CN102540252B (zh) 基于互相关的高精度中值叠加方法
CN105093319B (zh) 基于三维地震数据的地面微地震静校正方法
CN104483705A (zh) 一种三维剩余静校正方法
CN102692648A (zh) 一种基于大炮初至的折射波剩余静校正方法
CN106930770A (zh) 基于卷积神经网络的盾构机护盾间隙估计方法
CN101650443A (zh) 视电阻率的反向传播网络计算方法
CN104155694B (zh) 一种反射转换横波共检波点叠加剖面的剩余静校正方法
CN105204065A (zh) 一种拾取初至波的方法和装置
CN104749623B (zh) 一种地震资料成像处理方法
CN102043167B (zh) 一种曲面地表地震资料速度分析的方法
CN106443791B (zh) 求取倾斜地层或各向异性地层横波剩余静校正量的方法
CN107942385A (zh) 一种精细旅行时计算的动校正方法
CN106338766A (zh) 基于分步傅里叶算法的叠前时间偏移方法
CN105093280A (zh) 表层模型对地震数据影响的低频与高频成分的分解方法
CN109085642B (zh) 一种各向异性介质微地震事件定位方法
CN106226815B (zh) 转换波等效速度比谱的获取方法及转换波速度分析方法
CN107807356A (zh) 一种gpr绕射波速度分析方法
CN107976714B (zh) 一种复杂地表分级空间距离加权的超道集计算方法
CN107942373A (zh) 基于裂缝性油气储层断裂系统检测的相干算法
CN108508479B (zh) 一种空地井立体重磁数据协同目标位置反演方法
CN107179548B (zh) 一种基于真地表的叠前地震成像方法

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20201110

Address after: 100007 Beijing, Dongzhimen, North Street, No. 9, No.

Patentee after: CHINA NATIONAL PETROLEUM Corp.

Patentee after: BGP Inc., China National Petroleum Corp.

Address before: No. 189, fan Yangxi Road, Zhuozhou City, Baoding, Hebei

Patentee before: BGP Inc., China National Petroleum Corp.

TR01 Transfer of patent right