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

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

Info

Publication number
CN107688199A
CN107688199A CN201710747551.8A CN201710747551A CN107688199A CN 107688199 A CN107688199 A CN 107688199A CN 201710747551 A CN201710747551 A CN 201710747551A CN 107688199 A CN107688199 A CN 107688199A
Authority
CN
China
Prior art keywords
mrow
speed
msub
mtd
dimentional
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
CN201710747551.8A
Other languages
English (en)
Other versions
CN107688199B (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
Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
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 Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd filed Critical Geophysical Prospecting Co of CNPC Chuanqing Drilling Engineering Co Ltd
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)
  • Image Generation (AREA)
  • Magnetic Resonance Imaging Apparatus (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值如下表所示。
f0 f1 f2 f3 f4 f5 f6 f7 f8 f9 f10 f11
-0.5 -0.4 -0.3 -0.25 -0.2 -0.15 -0.1 -0.08 -0.06 -0.04 -0.02 0
f26 f25 f24 f23 f22 f21 f19 f17 f15 f14 f13 f12
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，其中，
<mrow> <msub> <mi>f</mi> <mrow> <mi>m</mi> <mo>+</mo> <mn>1</mn> </mrow> </msub> <mo>=</mo> <mfenced open = "{" close = ""> <mtable> <mtr> <mtd> <mrow> <msub> <mi>f</mi> <mi>m</mi> </msub> <mo>+</mo> <mn>0.1</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mn>0.3</mn> <mo>&lt;</mo> <mo>|</mo> <msub> <mi>f</mi> <mi>m</mi> </msub> <mo>|</mo> <mo>&amp;le;</mo> <mn>0.6</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>f</mi> <mi>m</mi> </msub> <mo>+</mo> <mn>0.05</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mn>0.1</mn> <mo>&lt;</mo> <mo>|</mo> <msub> <mi>f</mi> <mi>m</mi> </msub> <mo>|</mo> <mo>&amp;le;</mo> <mn>0.3</mn> </mrow> </mtd> </mtr> <mtr> <mtd> <mrow> <msub> <mi>f</mi> <mi>m</mi> </msub> <mo>+</mo> <mn>0.02</mn> <mo>,</mo> </mrow> </mtd> <mtd> <mrow> <mo>|</mo> <msub> <mi>f</mi> <mi>m</mi> </msub> <mo>|</mo> <mo>&amp;le;</mo> <mn>0.1</mn> </mrow> </mtd> </mtr> </mtable> </mfenced> <mo>,</mo> </mrow>
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 true CN107688199A (zh) 2018-02-13
CN107688199B 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)

Cited By (1)

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

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
李列 等: ""低信噪比条件下地震叠加速度的提取技术"", 《海洋地质前沿》 *

Cited By (2)

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

Also Published As

Publication number Publication date
CN107688199B (zh) 2019-03-29

Similar Documents

Publication Publication Date Title
CN108196305B (zh) 一种山地静校正方法
CN101598808B (zh) 一种提高地震资料成像质量的方法
CN109782339A (zh) 一种基于3D-DnCNN网络的叠后三维地震资料随机噪声压制方法
CN109696906B (zh) 基于小波修正贝叶斯卷积能量的水下机器人推进器故障诊断方法
CN104570076A (zh) 一种基于二分法的地震波初至自动拾取方法
CN103927730B (zh) 基于Primal Sketch修正及矩阵填充的图像降噪方法
CN107390266A (zh) 基于角道集的速度更新方法及叠前深度偏移速度建模方法
CN107688199B (zh) 一种适用于二维低信噪比资料的叠加速度场构建方法
CN104155694B (zh) 一种反射转换横波共检波点叠加剖面的剩余静校正方法
Zhang et al. Improving the quality of prestack inversion by prestack data conditioning
CN109884693B (zh) 自适应走向速度谱求取方法和系统
CN104656140A (zh) 一种基于垂直地震衰减规律的中值滤波方法以及系统
CN104155688A (zh) 一种高精度加权叠加方法
CN106226815B (zh) 转换波等效速度比谱的获取方法及转换波速度分析方法
CN107942389A (zh) 用于压制邻炮干扰的方法、系统和计算机可读介质
CN106371140B (zh) 一种提高中深层地震资料分辨率的方法
CN104122583B (zh) 一种拓宽地震数据频谱的方法和装置
CN106896408A (zh) 一种角度域叠前时间偏移方法
CN112213784B (zh) 复杂地表地震数据一次处理快速静校正方法
CN106569278A (zh) 一种多道相似相干速度谱计算方法
CN107942385B (zh) 一种精细旅行时计算的动校正方法
CN105487111B (zh) 一种自动方位角道集部分叠加方法
CN106338762A (zh) 基于数据域导向的模型道重建反聚焦成像方法
CN109884701A (zh) 地质体散射角度导向深度成像方法
CN106249286A (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: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company Ltd.

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

Applicant before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company 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: Dongfang Geophysical Exploration Co., Ltd., China Petrochemical Corp.

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

Applicant before: China National Petroleum Corporation Chuanqing Drilling Engineering Geophysical Exploration Company 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