CN107942385B - 一种精细旅行时计算的动校正方法 - Google Patents
一种精细旅行时计算的动校正方法 Download PDFInfo
- Publication number
- CN107942385B CN107942385B CN201711130848.6A CN201711130848A CN107942385B CN 107942385 B CN107942385 B CN 107942385B CN 201711130848 A CN201711130848 A CN 201711130848A CN 107942385 B CN107942385 B CN 107942385B
- Authority
- CN
- China
- Prior art keywords
- point
- dynamic correction
- shot
- root
- travel time
- 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
Links
- 238000012937 correction Methods 0.000 title claims abstract description 39
- 238000000034 method Methods 0.000 title claims abstract description 23
- 238000012545 processing Methods 0.000 claims abstract description 17
- 238000005070 sampling Methods 0.000 claims abstract description 15
- 239000000523 sample Substances 0.000 claims abstract description 14
- 238000001514 detection method Methods 0.000 claims abstract description 10
- 238000003384 imaging method Methods 0.000 abstract description 6
- 239000003208 petroleum Substances 0.000 description 3
- 239000004576 sand Substances 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000007812 deficiency Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
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
- 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
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 (2)
1.一种精细旅行时计算的动校正方法,其特征在于,所述动校正方法包括:
A、对于地震数据中任一采样点,获取炮点均方根速度vS和检波点均方根速度vR;
B、根据式1计算动校正的旅行时t,并进行动校正,
式1为:
在式1中,t0为炮点到检波点的自激自收时间,xS为炮点的偏移距,xR为检波点的偏移距;
C、重复所述步骤A和B,直到地震数据中所有采样点处理完毕;
其中,所述步骤A包括:获取炮点到反射样点路径上的所有均方根速度场并进行中值、均值和均方根计算中的一种或多种,得到所述炮点均方根速度vS;获取反射样点到接收点路径上的所有均方根速度场并进行中值、均值和均方根计算中的一种或多种,得到所述检波点均方根速度vR。
2.根据权利要求1所述的精细旅行时计算的动校正方法,其特征在于,所述动校正方法用于复杂地区的地震数据处理。
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 CN107942385A (zh) | 2018-04-20 |
CN107942385B true 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) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111856581B (zh) * | 2020-07-27 | 2022-02-22 | 广州海洋地质调查局 | 一种obs时钟漂移校正方法及处理终端 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101984366A (zh) * | 2010-09-29 | 2011-03-09 | 北京吉星吉达科技有限公司 | 基于共孔径面的起伏地表叠前时间偏移方法及设备 |
CN102162858A (zh) * | 2010-12-06 | 2011-08-24 | 中国海洋石油总公司 | 利用非对称走时进行动校正速度分析的方法 |
CN102353992B (zh) * | 2011-06-02 | 2013-04-03 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于中心点规则化的起伏地表动校正方法 |
CN102565857B (zh) * | 2011-12-16 | 2014-09-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 自动剩余动校正方法 |
CN104122582B (zh) * | 2014-06-24 | 2017-06-23 | 夏正元 | 利用叠加速度求取高精度地震波速度的方法 |
CN104769457A (zh) * | 2014-07-25 | 2015-07-08 | 杨顺伟 | 一种基于旅行时法确定地层裂隙裂缝走向的方法及装置 |
CN104216016B (zh) * | 2014-08-12 | 2017-01-25 | 中国石油天然气集团公司 | 一种纵波约束扫描的转换波动校正方法以及系统 |
CN106324673B (zh) * | 2015-07-03 | 2018-05-11 | 中国石油化工股份有限公司 | 一种用于进行全方位纵波速度分析的方法和装置 |
CN105866833A (zh) * | 2016-06-16 | 2016-08-17 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | Vsp-cdp叠加方法及三维vsp覆盖次数计算方法 |
CN107179548B (zh) * | 2017-07-10 | 2019-03-29 | 中国石油集团东方地球物理勘探有限责任公司 | 一种基于真地表的叠前地震成像方法 |
-
2017
- 2017-11-15 CN CN201711130848.6A patent/CN107942385B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN107942385A (zh) | 2018-04-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN101609160B (zh) | 一种抗背景噪音强的地震数据初至波自动拾取方法 | |
CN104732531B (zh) | 一种高分辨率遥感图像信噪比曲线自适应获取方法 | |
CN103293554B (zh) | 一种基于地震数据的静校正方法 | |
CN105259571B (zh) | 一种地层倾角检测方法 | |
CN105204065B (zh) | 一种拾取初至波的方法和装置 | |
CN105510975B (zh) | 提高地震数据信噪比的方法及装置 | |
CN109782335A (zh) | 低畸变宽频扫描信号设计方法 | |
CN102053260B (zh) | 获得地震纵波的方位速度的方法及处理地震数据的方法 | |
CN106483563A (zh) | 基于互补集合经验模态分解的地震能量补偿方法 | |
CN107942385B (zh) | 一种精细旅行时计算的动校正方法 | |
CN113740915B (zh) | 一种球坐标系重力和接收函数联合反演地壳结构参数的方法 | |
CN106990434B (zh) | 椭圆展开转换波成像方法及系统 | |
CN106814394B (zh) | 台站式检波器与常规检波器联合低频恢复方法 | |
CN107179548B (zh) | 一种基于真地表的叠前地震成像方法 | |
CN109581494B (zh) | 叠前偏移方法及装置 | |
CN112213784B (zh) | 复杂地表地震数据一次处理快速静校正方法 | |
CN105223630B (zh) | 基于地质模型的全方位观测系统参数论证方法 | |
CN110703327A (zh) | 一种全频带成像方法 | |
CN106569278B (zh) | 一种多道相似相干速度谱计算方法 | |
CN108363739B (zh) | 一种基于稀疏采集的地震资料高低频拓展方法 | |
CN106646614B (zh) | 面向时差校正的叠合波分离方法 | |
CN105445719B (zh) | 一种三维激光扫描仪数据滤波方法 | |
El-Emam et al. | Maximizing the value of Single-Sensor measurements, Kuwait experience | |
CN102778692B (zh) | 一种利用地震数据确定品质因数的方法 | |
CN106443778B (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 | ||
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. |