CN107688199B - 一种适用于二维低信噪比资料的叠加速度场构建方法 - Google Patents

一种适用于二维低信噪比资料的叠加速度场构建方法 Download PDF

Info

Publication number
CN107688199B
CN107688199B CN201710747551.8A CN201710747551A CN107688199B CN 107688199 B CN107688199 B CN 107688199B CN 201710747551 A CN201710747551 A CN 201710747551A CN 107688199 B CN107688199 B CN 107688199B
Authority
CN
China
Prior art keywords
speed
velocity
velocity field
low snr
dimensional
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
CN201710747551.8A
Other languages
English (en)
Other versions
CN107688199A (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 CN201710747551.8A priority Critical patent/CN107688199B/zh
Publication of CN107688199A publication Critical patent/CN107688199A/zh
Application granted granted Critical
Publication of CN107688199B publication Critical patent/CN107688199B/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/303—Analysis for determining velocity profiles or travel times
    • 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/62—Physical property of subsurface
    • G01V2210/622—Velocity, density or impedance
    • G01V2210/6222—Velocity; travel time

Landscapes

  • Engineering & Computer Science (AREA)
  • Remote Sensing (AREA)
  • Physics & Mathematics (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)
  • Magnetic Resonance Imaging Apparatus (AREA)
  • Image Generation (AREA)

Abstract

本发明提供了一种适用于二维低信噪比资料的叠加速度场构建方法。所述方法包括：对二维测线的共中心点道集进行常速扫描，按照从浅到深的顺序选择速度趋势，得到第一速度体Vij；对所述第一速度体Vij进行插值处理，得到第二速度体intVij；随后变比例生成第三速度体fmVij；对所述第三速度体fmVij中的每个分别进行叠加处理，根据叠加品质优选速度，然后插值构件叠加速度场。本发明能够实现通过二维低信噪比资料构建叠加速度场，能够有效解决以往依靠能量谱、超道集、叠加小剖面无法选取低信噪比资料速度的问题。

Description

一种适用于二维低信噪比资料的叠加速度场构建方法
技术领域
本发明涉及油气藏地震勘探资料处理领域，具体来讲，涉及一种适用于二维低信噪比资料的叠加速度场构建方法。
背景技术
近年来，油气勘探领域逐渐向复杂地区延伸，这些地区地震地表条件复杂，高、低速面波及散射发育，地震资料信噪比低，能量散失快，成像效果不理想。特别是在山体主体部位，由于出露灰岩，激发接收条件差，地腹波场(例如，中、深层波场)复杂，资料甚至难以成像。
资料叠加成像，是速度逐渐逼近的过程。常规的均方根速度拾取方法，是在能量谱上选取能量集聚区，应用该速度确保超道集拉平并在叠加小剖面上成像。叠加小剖面通常以某一叠加速度为准，按比例生成不同的速度，将控制点周围的地震道做叠加处理，形成不同的小叠加剖面。
然而，发明人发现：对于低信噪比资料，上面描述的常规方法不可行，最根本原因是由于信噪比过低，能量谱上没有集中的能量团供选择，超道集也没有明显的有效信号供拉平。
因此，本发明了一种适用于二维低信噪比资料的叠加速度场构建方法，以改善剖面叠加成像效果。
发明内容
本发明的目的在于解决现有技术存在的上述不足中的至少一项。
例如，本发明的目的之一在于解决二维低信噪比资料较难构建叠加速度场的问题。
为了实现上述目的，本发明提供了一种适用于二维低信噪比资料的叠加速度场构建方法，所述方法包括步骤：对二维测线的共中心点道集进行常速扫描，按照从浅到深的顺序选择速度趋势，得到第一速度体Vij；对所述第一速度体Vij进行插值处理，得到第二速度体intVij；按照式1变比例生成第三速度体fmVij，式1为fmVij＝(1-fm)×intVij，其中，
f0在-0.6～-0.4的范围内确定，m＝0,1,2,3…k，k为大于3的自然数；对所述第三速度体fmVij中的每个分别进行叠加处理，根据叠加品质优选速度，然后插值构件叠加速度场。
与现有技术相比，本发明的有益效果包括：能够提供一种适用于二维低信噪比资料的叠加速度场构建方法，能够有效解决以往依靠能量谱、超道集、叠加小剖面无法选取低信噪比资料速度的问题。
具体实施方式
在下文中，将结合示例性实施例来详细说明本发明的适用于二维低信噪比资料的叠加速度场构建方法。
在本发明的一个示例性实施例中，适用于二维低信噪比资料的叠加速度场构建方法可通过以下步骤来实现：
首先，进行等间隔常速扫描，确定速度趋势。具体来讲，对二维测线的共中心点道集进行常速扫描，按照从浅到深的顺序选择速度趋势，得到第一速度体Vij。这里，按照从浅到深的顺序选择速度趋势的步骤可包括：对于极浅层，参照静校正环节的大炮初至层析模型选择速度；对于中深层，根据叠加剖面的成像品质优劣选择。其中，极浅层深度即静校正环节的大炮初至层析反演得到的深度。此外，二维测线的共中心点道集优选为校正到浮动基准面，经去噪、振幅补偿、反褶积的预处理后的二维测线的共中心点道集。
然后，对所述第一速度体Vij进行插值处理，得到第二速度体intVij。
接下来，采用变比例扫描，优选速度。具体来讲，按照下面式1变比例生成第三速度体fmVij。
fmVij＝(1-fm)×intVij (1)
其中，
f0可以在-0.6～-0.4的范围(包含端点值)内选择一确定值。m＝0,1,2,3…k，k为自然数。k可以为16至30。例如，k可以在20～26的范围(包含端点值)内选择一确定值。
随后，对第三速度体fmVij中的每个(即，m的每个具体取值所对应的第三速度体)分别进行叠加处理，根据叠加品质优选速度，然后插值处理(例如，线性插值)以构件叠加速度场。根据叠加品质优选速度的过程可由本领域普通技术人员来完成。
在本发明的另一个示例性实施例中，适用于二维低信噪比资料的叠加速度场构建方法通过在常速扫描基础上，变比例扫描，改善叠加成像效果。
其处理对象可以为：校正到浮动基准面，经去噪、振幅补偿、反褶积等一系列预处理后的二维测线的共中心点道集。例如，常速扫面的过程可以为：根据经验给一个初始速度1stVel，按速度间隔velInt扫描得到每一轮的叠加速度vel，考虑复杂地区表层速度极低，本发明1stVel为1000m/s，velInt取200m/s，n为扫描个数，取30个，以保证扫描的最大速度为7000m/s。
vel＝1stVel+n*velInt (2)
应用n个常速度分别进行叠加，将叠加剖面按速度从小到大进行排序，按照从浅到深的顺序选择速度趋势。针对极浅层，参照静校正环节的大炮初至层析模型选择速度。中深层，根据成像品质优劣选择，得到一个初步的速度体Vij，这里i代表速度控制点索引，j是在该速度控制点从上到下的时间索引。
之后，在常速扫描基础上进行变比例扫描，优选速度。因为Vij在时间、空间上都是不等间距的，所以要对其进行插值(例如，线性插值)，获得新的速度体intVij，在该基础上按上面的式(1)变比例生成新的速度体fmVij。在本示例性实施例中，f0取-0.5，k取26，则相应的fm值如下表所示。
f<sub>0</sub> f<sub>1</sub> f<sub>2</sub> f<sub>3</sub> f<sub>4</sub> f<sub>5</sub> f<sub>6</sub> f<sub>7</sub> f<sub>8</sub> f<sub>9</sub> f<sub>10</sub> f<sub>11</sub>
-0.5 -0.4 -0.3 -0.25 -0.2 -0.15 -0.1 -0.08 -0.06 -0.04 -0.02 0
f<sub>26</sub> f<sub>25</sub> f<sub>24</sub> f<sub>23</sub> f<sub>22</sub> f<sub>21</sub> f<sub>19</sub> f<sub>17</sub> f<sub>15</sub> f<sub>14</sub> f<sub>13</sub> f<sub>12</sub>
0.52 0.42 0.32 0.27 0.22 0.17 0.12 0.1 0.08 0.06 0.04 0.02
在本发明的方法中，变比例(式1)扫描有两个优点，一是能确保扫描较宽的范围，可对上一次的扫描结果有充分的补充，二是在基于速度趋势，在接近上一次常速扫描结果时，缩小扫描步长，更为精准。m为变比例因子扫描次数。
应用m个速度体分别进行叠加处理，导入交互软件，按变比例因子m从小到大排序，根据叠加品质优选速度。优选的方式是应用鼠标左键在成像处单击，即时保存单击处的共中心点索引、时间索引和速度值，供后续输出成文件，然后插值构建叠加速度场。
综上所述，本发明的方法尤其适用于二维低信噪比资料的叠加速度场构建，其通过常速扫描和变比例扫描相结合，确保低信噪比资料的叠加速度从无到有，从有到优，可控性更强，更适合山地资料速度场变化剧烈的特性，可改善叠加成像效果。
尽管上面已经结合示例性实施例及描述了本发明，但是本领域普通技术人员应该清楚，在不脱离权利要求的精神和范围的情况下，可以对上述实施例进行各种修改。

Claims (4)

1.一种适用于二维低信噪比资料的叠加速度场构建方法，所述方法对二维测线的共中心点道集进行常速扫描，按照从浅到深的顺序选择速度趋势，得到第一速度体Vij，其特征在于，所述方法还包括步骤：
对所述第一速度体Vij进行插值处理，得到第二速度体intVij；
按照式1变比例生成第三速度体fmVij，式1为fmVij＝(1-fm)×intVij，其中，
f0在-0.6～-0.4的范围内确定，m＝0,1,2,3…k，k为大于3的自然数；
对所述第三速度体fmVij中的每个分别进行叠加处理，根据叠加品质优选速度，然后插值构件叠加速度场。
2.根据权利要求1所述的适用于二维低信噪比资料的叠加速度场构建方法，其特征在于，所述按照从浅到深的顺序选择速度趋势的步骤包括：对于极浅层，参照静校正环节的大炮初至层析模型选择速度；对于中深层，根据叠加剖面的成像品质优劣选择。
3.根据权利要求1所述的适用于二维低信噪比资料的叠加速度场构建方法，其特征在于，k在20～26的范围内确定。
4.根据权利要求1所述的适用于二维低信噪比资料的叠加速度场构建方法，其特征在于，所述二维测线的共中心点道集为校正到浮动基准面，经去噪、振幅补偿、反褶积的预处理后的二维测线的共中心点道集。
CN201710747551.8A 2017-08-28 2017-08-28 一种适用于二维低信噪比资料的叠加速度场构建方法 Active CN107688199B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710747551.8A CN107688199B (zh) 2017-08-28 2017-08-28 一种适用于二维低信噪比资料的叠加速度场构建方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710747551.8A CN107688199B (zh) 2017-08-28 2017-08-28 一种适用于二维低信噪比资料的叠加速度场构建方法

Publications (2)

Publication Number Publication Date
CN107688199A CN107688199A (zh) 2018-02-13
CN107688199B true CN107688199B (zh) 2019-03-29

Family

ID=61155520

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710747551.8A Active CN107688199B (zh) 2017-08-28 2017-08-28 一种适用于二维低信噪比资料的叠加速度场构建方法

Country Status (1)

Country Link
CN (1) CN107688199B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113514880B (zh) * 2020-04-10 2024-02-20 中国石油化工股份有限公司 一种二维测线交点处叠加速度构建方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0866343A2 (en) * 1997-03-21 1998-09-23 Exxon Production Research Company Method for determining seismic velocities
WO2005010797A2 (en) * 2003-07-23 2005-02-03 Lee Wook B Improved 3d veloctiy modeling, with calibration and trend fitting using geostatistical techniques, particularly advantageous for curved-ray prestack time migration and for such migration followed by prestack depth migration
CN102096103A (zh) * 2010-12-03 2011-06-15 中国石油天然气集团公司 用于低信噪比地震资料的速度分析方法
CN103226210A (zh) * 2013-04-01 2013-07-31 中国石油天然气股份有限公司 一种圈闭落实方法
CN104678440A (zh) * 2015-02-15 2015-06-03 山东科技大学 井约束二维地震变速度场非线性误差校正方法
CN104932015A (zh) * 2015-06-19 2015-09-23 中国石油集团川庆钻探工程有限公司地球物理勘探公司 构建地震数据的速度模型的方法
CN105891887A (zh) * 2016-04-11 2016-08-24 北京科胜伟达石油科技股份有限公司 基于叠加数据的速度纵横向高密度分析方法
CN106249290A (zh) * 2016-08-25 2016-12-21 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种利用多级数据融合建立表层速度结构模型的方法
CN106324679A (zh) * 2016-08-04 2017-01-11 北京金双狐油气技术有限公司 量版法求取层速度并建立平均速度场的方法
CN106772596A (zh) * 2016-12-08 2017-05-31 中国石油天然气集团公司 一种确定叠前时间偏移速度场的方法及装置

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0866343A2 (en) * 1997-03-21 1998-09-23 Exxon Production Research Company Method for determining seismic velocities
WO2005010797A2 (en) * 2003-07-23 2005-02-03 Lee Wook B Improved 3d veloctiy modeling, with calibration and trend fitting using geostatistical techniques, particularly advantageous for curved-ray prestack time migration and for such migration followed by prestack depth migration
CN102096103A (zh) * 2010-12-03 2011-06-15 中国石油天然气集团公司 用于低信噪比地震资料的速度分析方法
CN103226210A (zh) * 2013-04-01 2013-07-31 中国石油天然气股份有限公司 一种圈闭落实方法
CN104678440A (zh) * 2015-02-15 2015-06-03 山东科技大学 井约束二维地震变速度场非线性误差校正方法
CN104932015A (zh) * 2015-06-19 2015-09-23 中国石油集团川庆钻探工程有限公司地球物理勘探公司 构建地震数据的速度模型的方法
CN105891887A (zh) * 2016-04-11 2016-08-24 北京科胜伟达石油科技股份有限公司 基于叠加数据的速度纵横向高密度分析方法
CN106324679A (zh) * 2016-08-04 2017-01-11 北京金双狐油气技术有限公司 量版法求取层速度并建立平均速度场的方法
CN106249290A (zh) * 2016-08-25 2016-12-21 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种利用多级数据融合建立表层速度结构模型的方法
CN106772596A (zh) * 2016-12-08 2017-05-31 中国石油天然气集团公司 一种确定叠前时间偏移速度场的方法及装置

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
"低信噪比条件下地震叠加速度的提取技术";李列 等;《海洋地质前沿》;20121031;第28卷(第10期);第53-59页

Also Published As

Publication number Publication date
CN107688199A (zh) 2018-02-13

Similar Documents

Publication Publication Date Title
WO2006054181A1 (en) Method for processing at least two sets of seismic data
CN107843922A (zh) 一种基于地震初至波和反射波走时联合的层析成像方法
Trad et al. Fast and robust deblending using apex shifted Radon transform
CN104090301B (zh) 一种求取三维高频静校正量的方法
CN109856676B (zh) 一种实现地震共反射面叠加参数优化的方法
NO180026B (no) Fremgangsmåte for korreksjon av amplitudeverdier ved seismiske målinger
CN107229071B (zh) 一种地下构造反演成像方法
WO2005026776A1 (en) Wide-offset-range pre-stack depth migration method for seismic exploration
CN107688199B (zh) 一种适用于二维低信噪比资料的叠加速度场构建方法
CN110389377B (zh) 基于波形互相关系数相乘的微震偏移成像定位方法
CN104316961B (zh) 获取风化层的地质参数的方法
CN106054252B (zh) 一种叠前时间偏移的方法及装置
Abedi et al. Three-parameter normal moveout correction in layered anisotropic media: A stretch-free approach
CN106199704A (zh) 一种三维三分量海底电缆地震资料速度建模方法
AU2004208729A1 (en) Method for seismic migration using explicit depth extrapolation operators with dynamically variable operator length
CN110716231B (zh) 基于共聚焦域稀疏反演的海上多震源波场分离方法及系统
US5067113A (en) Efficient generation of traveltime tables for two dimensional and three dimensional prestack depth migration
AU2003295905B2 (en) Method for processing seismic traces to provide seismic data with enhanced signal-to-noise ratio
CN106896408B (zh) 一种角度域叠前时间偏移方法
CN111766626A (zh) 一种实现地震数据三维空间拟合的算法
CN106646614B (zh) 面向时差校正的叠合波分离方法
CA2533011C (en) Method for bispectral picking of anelliptical nmo correction parameters
CN105425300A (zh) 一种剩余静校正方法
CN106569278B (zh) 一种多道相似相干速度谱计算方法
CN107942385B (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
CB02 Change of applicant information

Address after: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan

Applicant after: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

Address before: 610213 No. 1, No. 1, No. 1, Huayang Avenue, Huayang Town, Shuangliu County, Chengdu, Sichuan

Applicant before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

CB02 Change of applicant information
TA01 Transfer of patent application right

Effective date of registration: 20180402

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

Applicant after: BGP INC., CHINA NATIONAL PETROLEUM Corp.

Address before: No. 216, No. 216, Huayang Avenue, Tianfu New District, Sichuan, Sichuan

Applicant before: GEOPHYSICAL EXPLORATION COMPANY OF CNPC CHUANQING DRILLING ENGINEERING Co.,Ltd.

TA01 Transfer of patent application right
GR01 Patent grant
GR01 Patent grant
TR01 Transfer of patent right

Effective date of registration: 20201116

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