CN105891887A - 基于叠加数据的速度纵横向高密度分析方法 - Google Patents

基于叠加数据的速度纵横向高密度分析方法 Download PDF

Info

Publication number
CN105891887A
CN105891887A CN201610220911.4A CN201610220911A CN105891887A CN 105891887 A CN105891887 A CN 105891887A CN 201610220911 A CN201610220911 A CN 201610220911A CN 105891887 A CN105891887 A CN 105891887A
Authority
CN
China
Prior art keywords
velocity
different
imaging
high density
data
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
CN201610220911.4A
Other languages
English (en)
Other versions
CN105891887B (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.)
Beijing Co-Sail Petroleum Technology CoLtd
Original Assignee
Beijing Co-Sail Petroleum Technology CoLtd
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 Beijing Co-Sail Petroleum Technology CoLtd filed Critical Beijing Co-Sail Petroleum Technology CoLtd
Priority to CN201610220911.4A priority Critical patent/CN105891887B/zh
Publication of CN105891887A publication Critical patent/CN105891887A/zh
Application granted granted Critical
Publication of CN105891887B publication Critical patent/CN105891887B/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

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)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明公开了一种基于叠加数据速度纵横向高密度分析的方法，包括以下步骤：建立原始速度模型；按照不同比例形成新的速度模型；按照纵线速度间隔和横线速度间隔，得到不同速度的纵线成像叠加数据体和横线成像叠加数据体；分别把不同速度模型、对应不同速度的纵线成像结果、对应不同速度的横线成像结果输入到显示窗口中；把对应的测井声波曲线输入到分析窗口中；合成记录制作；结合合成记录根据地质认识标定主要反射层系；对比不同速度下纵波成像结果，确定对应位置速度值；完成纵线、横线速度值拾取；平滑新的全区速度模型得到新的成像速度模型；迭代速度分析，得到最终的成像速度模型。可以提高叠加速度或者偏移速度的准确率和效率。

Description

基于叠加数据的速度纵横向高密度分析方法
技术领域
本发明属于地震资料处理速度分析技术领域，具体地涉及一种基于叠加数据的速度纵横向高密度分析方法。
背景技术
叠加数据体高密度速度分析原理是基于成像结果和井信息标定相结合的一种地震资料处理速度分析流程。在地震勘探中，随着寻找大规模构造圈闭越来越少，对微幅构造和隐蔽油气藏越来越重视，寻找这类小构造需要地震资料高精度成像，在地震成像方法确定的情况下，提高地震处理的速度精度成为提高成像精度的重要方法之一。目前的速度分析方法一般是以共中心点道集或者成像道集为基础，采用动校正公式或者剩余校正公式，把道集中的同相轴校平，最后形成新的速度模型。由于以上速度分析中缺乏叠加剖面信息和测井数据信息，仅仅以同相轴是否拉平作为唯一标准，使得速度分析判别标准比较单一，速度精度往往低于3%，因此迫切需要一种新的速度分析方法，需要综合叠加成果、测井数据等信息，使得速度精度误差低于1%。
发明内容
针对上述技术问题，本发明旨在提供一种基于叠加数据的速度纵横向高密度分析方法，可以提高地震资料处理中速度分析的精度，并大大提高速度分析的效率。
为达到上述目的，本发明的技术方案是：
一种基于叠加数据的速度纵横向高密度分析方法，其特征在于，包含以下步骤：
S01：建立地质成像的初始速度模型；
S02：按照不同比例参数，在初始速度模型的基础上形成不同比例的速度模型；
S03：根据不同工区和速度分析阶段，采用不同的纵线间隔和横线间隔，利用步骤S02的不同比例的速度模型分别做成像处理；
S04：把步骤S03的成像结果按照纵线和横线进行排列，并对成像结果进行可视化显示；
S05：在地质研究基础下对测井声波曲线进行编辑；
S06：利用整理好的声波测井数据制作合成记录，并进行可视化显示；
S07：根据实际成像结果与井旁道合成地震记录进行对比，确定合成记录和成像结果的对应关系，进行相关层位标定；
S08：在全区按照纵线和横线根据标定对应关系进行速度拾取；
S09：对离散点进行速度插值形成新的速度模型；
S10：对成像结果进行综合分析，如果满足地质要求，则输出速度模型；如果不满足地质要求，则把该速度模型作为新的初始速度模型重复步骤S02，直到满足要求为止。
优选的，所述步骤S01中所述初始速度模型是叠加速度或者DMO速度，经平滑后得到的速度模型。
优选的，所述步骤S02中不同比例速度模型是针对初始速度模型或者修正后的速度模型按照不同的比例系数相乘得到的不同比例的速度模型，根据分析速度精度和处理进度，在不同阶段比例系数不同，先期比例间隔为3-5%，后期比例间隔在0.5-1%之间。
优选的，所述步骤S03在初始速度分析阶段，纵线间隔和横线间隔分析点间隔为1000米X1000米，随着速度分析迭代的推进，最终的速度分析点间隔为200米X200米。
优选的，所述步骤S05还包括删除异常值，所述异常值为大于或者小于平均值而又与地质分层不相符的数据。
优选的，所述步骤S06中的合成记录为利用测井声波测井数据，根据理论子波合成井旁道理论地震数据记录，所选子波为雷克子波，主频根据不同地震资料特性进行选择。
优选的，所述步骤S08的速度拾取为根据地质含义和标定对比，确定合适的成像速度。
与现有技术相比，本发明的有益效果是：
在本发明中，根据一组不同速度模型在纵、横线上成像不同，反映出构造形态，断层位置、断距，储层厚度等地质方面的差异，确定优选速度；结合地质认识和地震合成记录进行对比对速度进一步优化，确定速度细节；由于是结合合成记录在叠加数据体上进行速度分析和拾取，大大提高生产效率，经过估算，拾取速度的效率提高了40-60%。本发明中的基于叠加数据的速度纵横向高密度分析方法在地震资料处理中有广泛的应用前景。
附图说明
图1为本发明基于叠加数据的速度纵横向高密度分析方法的流程图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚明了，下面结合具体实施方式以及附图，对本发明进一步详细说明。应该理解，这些描述只是示例性的，而并非要限制本发明的范围。此外，在以下说明中，省略了对公知结构和技术的描述，以避免不必要地混淆本发明的概念。
实施例：
如图1所示，一种基于叠加数据的速度纵横向高密度分析方法，包含以下步骤：
在步骤101，建立初始速度模型，初始速度模型是叠加速度或者DMO速度，经平滑后得到的速度模型，该速度模型是一二进制文件，每个采样点对应一个速度值，能够大体反映地下构造形态。流程引入步骤102.
在步骤102-1、102-2、…、102-N中，计算不同比例速度模型，是针对初始速度模型或者修正后的速度模型按照不同与不同的比例系数相乘得到的不同比例的速度模型。根据分析速度精度和处理进度，在不同阶段比例系数不同，先期比例为0.9、0.95、1.0、1.05、1.10，中期为0.94、0.97、1.00、1.03、1.06，后期比例系数为0.98、0.99、1.00、1.01、1.02，根据需要还可以细化为0.99、0.995、1.000、1.005、1.015等。形成不同比例数据体后流程进入步骤103。
在步骤103-1、103-2、…、103-N中利用以上不同速度模型，采用不同的纵线间隔和横线间隔进行成像处理，在初始速度分析阶段，分析点间隔一般为1000米X1000米，随着速度分析迭代的推进，最终的速度分析点间隔为200米X200米，如果需要可在构造复杂区进行加密分析，达到100米X100米的密度，即进行所称的高密度速度分析，把成像结果输入到步骤104中。
在步骤104-1中按照纵线抽取不同速度的成像结果，在步骤104-2中按照横线抽取不同速度的成像结果，然后输入到步骤105中进行可视化显示。
在步骤106中，在地质研究基础下对测井声波曲线进行编辑，删除异常值，这些异常值主要是指明显大于或者小于平均值而又与地质分层不相符的数据，完成数据整理后，把对应的测井声波曲线输入到分析窗口中，输入到步骤107中进行合成记录制作。
在步骤107中，利用整理好的声波测井数据，根据理论子波合成井旁道理论地震数据记录，这里所选子波为雷克子波，主频根据不同地震资料特性进行选择，把生成好的合成记录输入显示窗口，进行步骤108工作。
在步骤108中结合地质研究和层位划分，把合成记录和不同速度成像结果相结合，进行井震标定，确定标志层、目的层等，一旦确定了井震关系，就可以进行速度拾取，进入步骤109。
在步骤109和110中，根据地质含义和井震标定，分别在纵线和横线上进行速度拾取，开始阶段是进行稀疏速度拾取，后续阶段进行高密度速度拾取，然后进行步骤111。
在步骤111中对新拾取的速度进行插值，每点对应一个速度值，生成二进制文件备用。
步骤112中，结合步骤109和步骤110对成像结果进行综合分析，判断是否进行新的一轮速度分析，如果成像结果中主要地质层位满足根据钻井信息所标定的实际地质层位则把步骤111中生成的速度模型输出，作为最终的成像速度，如果不满足地质要求，则把步骤111中的速度模型作为初始速度（修正后的速度）输入步骤102中进行新一轮的速度分析。
在本发明中，根据基于叠加数据的速度纵横向高密度分析方法得到的速度可广泛应用于地震资料处理的各个阶段，运行效率高，能够大大提高资料成像的精度。
应当理解的是，本发明的上述具体实施方式仅仅用于示例性说明或解释本发明的原理，而不构成对本发明的限制。因此，在不偏离本发明的精神和范围的情况下所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。此外，本发明所附权利要求旨在涵盖落入所附权利要求范围和边界、或者这种范围和边界的等同形式内的全部变化和修改例。

Claims (7)

1.一种基于叠加数据的速度纵横向高密度分析方法，其特征在于，包含以下步骤：
S01：建立地质成像的初始速度模型；
S02：按照不同比例参数，在初始速度模型的基础上形成不同比例的速度模型；
S03：根据不同工区和速度分析阶段，采用不同的纵线间隔和横线间隔，利用步骤S02的不同比例的速度模型分别做成像处理；
S04：把步骤S03的成像结果按照纵线和横线进行排列，并对成像结果进行可视化显示；
S05：在地质研究基础下对测井声波曲线进行编辑；
S06：利用整理好的声波测井数据制作合成记录，并进行可视化显示；
S07：根据实际成像结果与井旁道合成地震记录进行对比，确定合成记录和成像结果的对应关系，进行相关层位标定；
S08：在全区按照纵线和横线根据标定对应关系进行速度拾取；
S09：对离散点进行速度插值形成新的速度模型；
S10：对成像结果进行综合分析，如果满足根据钻井信息所标定的实际地质层位，则输出速度模型；如果不满足根据钻井信息所标定的实际地质层位，则把该速度模型作为新的初始速度模型重复步骤S02，直到满足要求为止。
2.根据权利要求1所述的基于叠加数据的速度纵横向高密度分析方法，其特征在于，所述步骤S01中所述初始速度模型是叠加速度或者DMO速度，经平滑后得到的速度模型。
3.根据权利要求1所述的基于叠加数据的速度纵横向高密度分析方法，其特征在于，所述步骤S02中不同比例速度模型是针对初始速度模型或者修正后的速度模型按照不同的比例系数相乘得到的不同比例的速度模型，根据分析速度精度和处理进度，在不同阶段比例系数不同，先期比例间隔为3-5%，后期比例间隔为0.5-1%。
4.根据权利要求1所述的基于叠加数据的速度纵横向高密度分析方法，其特征在于，所述步骤S03在初始速度分析阶段，纵线间隔和横线间隔分析点间隔为1000米X1000米，随着速度分析迭代的推进，最终的速度分析点间隔为200米X200米。
5.根据权利要求1所述的基于叠加数据的速度纵横向高密度分析方法，其特征在于，所述步骤S05还包括删除异常值，所述异常值为大于或者小于平均值而又与地质分层不相符的数据。
6.根据权利要求1所述的基于叠加数据的速度纵横向高密度分析方法，其特征在于，所述步骤S06中的合成记录为利用测井声波测井数据，根据理论子波合成井旁道理论地震数据记录，所选子波为雷克子波，主频根据不同地震资料特性进行选择。
7.根据权利要求1所述的基于叠加数据的速度纵横向高密度分析方法，其特征在于，所述步骤S08的速度拾取为根据地质含义和标定对比，确定合适的成像速度。
CN201610220911.4A 2016-04-11 2016-04-11 基于叠加数据的速度纵横向高密度分析方法 Active CN105891887B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201610220911.4A CN105891887B (zh) 2016-04-11 2016-04-11 基于叠加数据的速度纵横向高密度分析方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201610220911.4A CN105891887B (zh) 2016-04-11 2016-04-11 基于叠加数据的速度纵横向高密度分析方法

Publications (2)

Publication Number Publication Date
CN105891887A true CN105891887A (zh) 2016-08-24
CN105891887B CN105891887B (zh) 2018-06-01

Family

ID=57012325

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201610220911.4A Active CN105891887B (zh) 2016-04-11 2016-04-11 基于叠加数据的速度纵横向高密度分析方法

Country Status (1)

Country Link
CN (1) CN105891887B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107688199A (zh) * 2017-08-28 2018-02-13 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种适用于二维低信噪比资料的叠加速度场构建方法
CN111797508A (zh) * 2020-06-12 2020-10-20 中冶建筑研究总院有限公司 一种基于监测技术的钢屋架安全性实时评定方法

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030030565A1 (en) * 2001-08-07 2003-02-13 Nsk Ltd. Wireless sensor, rolling bearing with sensor, management apparatus and monitoring system
US20090290451A1 (en) * 2006-06-02 2009-11-26 Banik Niranjan C Subsalt velocity model building
CN102565853A (zh) * 2011-12-20 2012-07-11 中国石油集团川庆钻探工程有限公司地球物理勘探公司 对地质构造的速度模型进行建模的方法
CN104932015A (zh) * 2015-06-19 2015-09-23 中国石油集团川庆钻探工程有限公司地球物理勘探公司 构建地震数据的速度模型的方法

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20030030565A1 (en) * 2001-08-07 2003-02-13 Nsk Ltd. Wireless sensor, rolling bearing with sensor, management apparatus and monitoring system
US20090290451A1 (en) * 2006-06-02 2009-11-26 Banik Niranjan C Subsalt velocity model building
CN102565853A (zh) * 2011-12-20 2012-07-11 中国石油集团川庆钻探工程有限公司地球物理勘探公司 对地质构造的速度模型进行建模的方法
CN104932015A (zh) * 2015-06-19 2015-09-23 中国石油集团川庆钻探工程有限公司地球物理勘探公司 构建地震数据的速度模型的方法

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘和年 等: "地层倾角约束自适应孔径叠前时间偏移", 《石油地球物理勘探》 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107688199A (zh) * 2017-08-28 2018-02-13 中国石油集团川庆钻探工程有限公司地球物理勘探公司 一种适用于二维低信噪比资料的叠加速度场构建方法
CN107688199B (zh) * 2017-08-28 2019-03-29 中国石油集团东方地球物理勘探有限责任公司 一种适用于二维低信噪比资料的叠加速度场构建方法
CN111797508A (zh) * 2020-06-12 2020-10-20 中冶建筑研究总院有限公司 一种基于监测技术的钢屋架安全性实时评定方法
CN111797508B (zh) * 2020-06-12 2024-03-01 中冶建筑研究总院有限公司 一种基于监测技术的钢屋架安全性实时评定方法

Also Published As

Publication number Publication date
CN105891887B (zh) 2018-06-01

Similar Documents

Publication Publication Date Title
CN111239802B (zh) 基于地震反射波形和速度谱的深度学习速度建模方法
CN105842736B (zh) 一种构建地层速度模型的方法
CN108333628B (zh) 基于正则化约束的弹性波最小二乘逆时偏移方法
CA2920008C (en) Method and device for the generation and application of anisotropic elastic parameters in horizontal transverse isotropic (hti) media
US7970546B1 (en) Diplet-based imaging of seismic data in shot or receiver records
EP3507626B1 (en) Attenuation of multiple reflections
CN107831542B (zh) Ddw高精度深度域井震匹配方法
CN107817526B (zh) 叠前地震道集分段式振幅能量补偿方法及系统
CN102841375A (zh) 一种复杂条件下基于角度域共成像点道集的层析速度反演方法
JP2011509412A (ja) 地震探査データのスペクトルシェーピングインバージョン法及びマイグレーション法
WO2017035104A1 (en) Velocity model seismic static correction
US8139441B2 (en) Method for bispectral picking of anelliptical NMO correction parameters
Padhi et al. 2-D ocean temperature and salinity images from pre-stack seismic waveform inversion methods: An example from the South China Sea
CN102937720A (zh) 井控提高地震资料分辨率的方法
CN111983683B (zh) 一种少井条件下盐下湖相灰岩储层预测方法和系统
CN109507726A (zh) 时间域弹性波多参数全波形的反演方法及系统
Talukdar et al. Sub-basalt imaging of hydrocarbon-bearing Mesozoic sediments using ray-trace inversion of first-arrival seismic data and elastic finite-difference full-wave modeling along Sinor–Valod profile of Deccan Syneclise, India
CN105891887A (zh) 基于叠加数据的速度纵横向高密度分析方法
CN113031069A (zh) 一种岩溶地区的多元信息约束智能层析静校正方法
CN114861515A (zh) 层速度数据体的计算方法、装置、设备及介质
CN107561581B (zh) 基于相关系数法建立多井模型的方法
CN107526102B (zh) 纵波与转换波联合偏移速度建模方法和装置
CN105527648A (zh) 用于各向异性参数反演的敏感度矩阵的计算方法及系统
CN111175822B (zh) 改进直接包络反演与扰动分解的强散射介质反演方法
Eppinger et al. 2d near‐surface full‐waveform tomography reveals bedrock controls on critical zone architecture

Legal Events

Date Code Title Description
C06 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