CN112147686A - 多期发育火成岩叠前深度偏移成像速度建模方法及系统 - Google Patents
多期发育火成岩叠前深度偏移成像速度建模方法及系统 Download PDFInfo
- Publication number
- CN112147686A CN112147686A CN201910577201.0A CN201910577201A CN112147686A CN 112147686 A CN112147686 A CN 112147686A CN 201910577201 A CN201910577201 A CN 201910577201A CN 112147686 A CN112147686 A CN 112147686A
- Authority
- CN
- China
- Prior art keywords
- initial
- igneous rock
- layer
- velocity
- bottom interface
- 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
Links
- 239000011435 rock Substances 0.000 title claims abstract description 165
- 238000003384 imaging method Methods 0.000 title claims abstract description 85
- 238000013508 migration Methods 0.000 title claims abstract description 59
- 230000005012 migration Effects 0.000 title claims abstract description 59
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000011161 development Methods 0.000 title claims abstract description 15
- 238000012937 correction Methods 0.000 claims abstract description 45
- 239000010410 layer Substances 0.000 claims description 102
- 239000011229 interlayer Substances 0.000 claims description 13
- 238000012545 processing Methods 0.000 claims description 11
- 230000005856 abnormality Effects 0.000 abstract description 4
- 230000002159 abnormal effect Effects 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 230000008901 benefit Effects 0.000 description 4
- 238000004587 chromatography analysis Methods 0.000 description 4
- 230000015572 biosynthetic process Effects 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 208000035126 Facies Diseases 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000000704 physical effect Effects 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/303—Analysis for determining velocity profiles or travel times
-
- 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
-
- 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/51—Migration
- G01V2210/512—Pre-stack
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
公开了一种多期发育火成岩叠前深度偏移成像速度建模方法及系统。该方法可以包括:步骤1:针对地震数据,获得初始层速度模型;步骤2:针对初始层速度模型进行校正,获得校正底界面与校正层速度;步骤3:根据校正底界面与校正层速度,获得校正层速度模型;步骤4:以校正层速度模型替换初始层速度模型,重复步骤1‑4,直至火成岩底界地层的反射地震同相轴拉平,获得最终的校正层速度模型;步骤5:根据最终的校正层速度模型,获得火成岩叠前深度偏移成像速度模型。本发明通过建立准确的深度域火成岩层速度模型,进行准确的地震叠前深度偏移成像,消除由于火成岩速度不准在地震剖面上引起的假断裂以及地震同相轴上提或下拉的异常。
Description
技术领域
本发明涉及油气藏地震勘探领域,更具体地,涉及一种多期发育火成岩叠前深度偏移成像速度建模方法及系统。
背景技术
火成岩油气藏勘探是目前油气勘探的热点之一,在我国东部和中西部油田存在大量的火成岩油气藏。但是相比沉积岩,火成岩的形成非常复杂,很多地区的火成岩往往多期发育,而且每次发育的厚度不同,火成岩多期发育区域受到断层通道的控制,火成岩的空间分布变化剧烈。在对火成岩油气藏进行地震勘探的过程中,由于火成岩的纵波速度的在纵向和横向上都变化剧烈,在地震数据剖面上产生了大量的局部异常,导致下方地层及目的层的地震反射同相轴信噪比降低,在地震勘探时间剖面中有可能解释出假的构造,影响了进一步的解释工作,增加了勘探的风险。为了解决相应的时间地震剖面上构造假象,需要在深度域进行精确的纵波地震速度建模,并采用精确的叠前深度偏移成像方法对地震勘探数据进行成像。
对深度偏移成像来说,成像精度很大程度上取决于速度模型的精度。但是,当前的地震处理流程都是从时间域的共中心点道集上求取叠加速度,然后进行时间偏移迭代生成共反射点道集来拾取均方根速度,最后从均方根速度来求取地层的层速度。如果地下介质为层状而且速度横向变化不大,则上面的方法可以得到较为准确的层速度。但是对于火成岩地层,岩相复杂多变,其形态难以刻画,很难建立准确的火成岩深度域速度模型。火成岩岩性和岩相的复杂多变造成地层速度横向变化剧烈,反射层析中的射线路径非常复杂,造成速度建模的多解性。另外火成岩深度成像的共成像点道集在下方地层会产生杂乱反射,因此共成像点道集的剩余时差能量谱上会产生大量假象。由于常规的网格层析方法都是利用剩余时差能量普拾取剩余时差,然后利用剩余时差来更新速度模型,错误的剩余时差必然会在速度模型中产生更大的误差。因此,常规的深度域速度建模方法对于火成岩多期发育地区是失效的。因此,有必要开发一种多期发育火成岩叠前深度偏移成像速度建模方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种多期发育火成岩叠前深度偏移成像速度建模方法及系统,其能够通过建立准确的深度域火成岩层速度模型,对于火成岩多期发育地区的地震数据能够进行准确的地震叠前深度偏移成像,消除之前由于火成岩速度不准在地震剖面上引起的假断裂以及地震同相轴上提或下拉的异常。
根据本发明的一方面,提出了一种多期发育火成岩叠前深度偏移成像速度建模方法。所述方法可以包括:步骤1:针对地震数据,获得初始层速度模型;步骤2:针对所述初始层速度模型进行校正,获得校正底界面与校正层速度;步骤3:根据所述校正底界面与所述校正层速度,获得校正层速度模型;步骤4:以所述校正层速度模型替换所述初始层速度模型,重复步骤1-4,直至火成岩底界地层的反射地震同相轴拉平,获得最终的校正层速度模型;步骤5:根据所述最终的校正层速度模型,获得火成岩叠前深度偏移成像速度模型。
优选地,所述步骤1包括:针对地震数据进行时间域处理,计算地下地层的均方根速度,获得时间域偏移叠加剖面;根据所述时间域偏移叠加剖面,确定火成岩的顶界和底界,进而计算层间双程旅行时;将所述均方根速度转换为初始层速度,获得初始层速度模型。
优选地,将火成岩的底界和定界相减,获得火成岩的层间双程旅行时。
优选地,所述步骤2包括:针对所述初始层速度模型进行目标线偏移,优化火成岩之上的地层层速度,输出优化后的共成像点道集;根据所述优化后的共成像点道集,拾取火成岩深度域的初始顶界面和初始底界面,在所述初始底界面上拾取深度误差,进而获得校正底界面;根据所述校正底界面与所述初始顶界面,获得校正层速度。
优选地,所述根据所述优化后的共成像点道集,拾取火成岩深度域的初始顶界面和初始底界面,在所述初始底界面上拾取深度误差,进而获得校正底界面包括:将所述优化后的共成像点道集进行近偏移距和远偏移距叠加,获得近偏移距叠加剖面与远偏移距叠加剖面;根据所述优化后的共成像点道集与全偏移距叠加剖面,拾取火成岩深度域的初始顶界面和初始底界面;比较所述近偏移距叠加剖面与所述远偏移距叠加剖面,在所述近偏移距叠加剖面上沿着所述初始底界面上拾取深度误差,进而获得校正底界面。
根据本发明的另一方面,提出了一种多期发育火成岩叠前深度偏移成像速度建模系统,其特征在于,该系统包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:步骤1:针对地震数据,获得初始层速度模型;步骤2:针对所述初始层速度模型进行校正,获得校正底界面与校正层速度;步骤3:根据所述校正底界面与所述校正层速度,获得校正层速度模型;步骤4:以所述校正层速度模型替换所述初始层速度模型,重复步骤1-4,直至火成岩底界地层的反射地震同相轴拉平,获得最终的校正层速度模型;步骤5:根据所述最终的校正层速度模型,获得火成岩叠前深度偏移成像速度模型。
优选地,所述步骤1包括:针对地震数据进行时间域处理,计算地下地层的均方根速度,获得时间域偏移叠加剖面;根据所述时间域偏移叠加剖面,确定火成岩的顶界和底界,进而计算层间双程旅行时;将所述均方根速度转换为初始层速度,获得初始层速度模型。
优选地,将火成岩的底界和定界相减,获得火成岩的层间双程旅行时。
优选地,所述步骤2包括:针对所述初始层速度模型进行目标线偏移,优化火成岩之上的地层层速度,输出优化后的共成像点道集;根据所述优化后的共成像点道集,拾取火成岩深度域的初始顶界面和初始底界面;在所述初始底界面上拾取深度误差,进而获得校正底界面;根据所述校正底界面与所述初始顶界面,获得校正层速度。
优选地,所述根据所述优化后的共成像点道集,拾取火成岩深度域的初始顶界面和初始底界面,在所述初始底界面上拾取深度误差,进而获得校正底界面包括:将所述优化后的共成像点道集进行近偏移距和远偏移距叠加,获得近偏移距叠加剖面与远偏移距叠加剖面;根据所述优化后的共成像点道集与全偏移距叠加剖面,拾取火成岩深度域的初始顶界面和初始底界面;比较所述近偏移距叠加剖面与所述远偏移距叠加剖面,在所述近偏移距叠加剖面上沿着所述初始底界面上拾取深度误差,进而获得校正底界面。
本发明的方法和装置具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的多期发育火成岩叠前深度偏移成像速度建模方法的步骤的流程图。
图2a和图2b分别示出了根据本发明的一个实施例的应用本方法前后的叠前深度偏移成像结果的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的多期发育火成岩叠前深度偏移成像速度建模方法的步骤的流程图。
在该实施例中,根据本发明的多期发育火成岩叠前深度偏移成像速度建模方法可以包括:步骤1:针对地震数据,获得初始层速度模型;步骤2:针对初始层速度模型进行校正,获得校正底界面与校正层速度;步骤3:根据校正底界面与校正层速度,获得校正层速度模型;步骤4:以校正层速度模型替换初始层速度模型,重复步骤1-4,直至火成岩底界地层的反射地震同相轴拉平,获得最终的校正层速度模型;步骤5:根据最终的校正层速度模型,获得火成岩叠前深度偏移成像速度模型。
在一个示例中,步骤1包括:针对地震数据进行时间域处理,计算地下地层的均方根速度,获得时间域偏移叠加剖面;根据时间域偏移叠加剖面,确定火成岩的顶界和底界,进而计算层间双程旅行时;将均方根速度转换为初始层速度,获得初始层速度模型。
在一个示例中,将火成岩的底界和定界相减,获得火成岩的层间双程旅行时。
在一个示例中,步骤2包括:针对初始层速度模型进行目标线偏移,优化火成岩之上的地层层速度,输出优化后的共成像点道集;根据优化后的共成像点道集,拾取火成岩深度域的初始顶界面和初始底界面;在初始底界面上拾取深度误差,进而获得校正底界面;根据校正底界面与初始顶界面,获得校正层速度。
在一个示例中,根据优化后的共成像点道集,拾取火成岩深度域的初始顶界面和初始底界面,在初始底界面上拾取深度误差,进而获得校正底界面包括:将优化后的共成像点道集进行近偏移距和远偏移距叠加,获得近偏移距叠加剖面与远偏移距叠加剖面;根据优化后的共成像点道集与全偏移距叠加剖面,拾取火成岩深度域的初始顶界面和初始底界面;比较近偏移距叠加剖面与远偏移距叠加剖面,在近偏移距叠加剖面上沿着初始底界面上拾取深度误差,进而获得校正底界面。
具体地,火成岩是在高温地质条件下由岩浆或熔岩流冷凝结晶而成的岩石,通常分为喷发岩(火山岩)和侵入岩两大类。侵入岩和喷出岩又可划分为多种亚类,每种亚类由于成份和岩相的不同,其地震成像速度变化很大。火成岩一般致密坚硬,与沉积岩在岩石物理性质方面相差很大,具有密度大、电阻率高、磁化率高等特点。其地震波速度普遍高于围岩,与围岩之间常存在明显的密度差,因此在其顶面可形成一个较强的波阻抗反射界面,常在常规的地震剖面上形成低频中强反射,对下部岩层常产生屏蔽作用。火成岩对下伏地层构造的影响程度与多种因素有关:如火成岩的形态、厚度、速度及其与围岩速度的差值等。火成岩体厚度越大、速度越高及其与围岩的速度差值越大,下伏地层反射同相轴的上移越大;高速火成岩体厚度的突变可导致下伏地层反射同相轴呈凸凹不平的反射特征,对于本身比较弱的同相轴,其连续性变差,有假断层产生,形成杂乱反射。要想对火成岩及下部地层精确成像,需要在地震剖面上对火成岩进行识别,然后进行地质形态刻画,同时确定相应的速度,然后开展深度域速度建模和优化。
根据本发明的多期发育火成岩叠前深度偏移成像速度建模方法可以包括:
步骤1:针对地震数据进行时间域处理,计算地下地层的均方根速度,获得时间域偏移叠加剖面,均方根速度是基于地下地层为水平层状的假设,不考虑地层倾角的影响,主要用来求取火成岩周围沉积地层的层速度,时间域的地震偏移叠加剖面用来作为刻画火成岩形态的初始参考信息;
收集工区内或附近的测井数据,包括声波测井曲线以及井分层等信息,以此为基础,有经验的地质解释人员可以对时间域偏移叠加剖面进行精确的层位解释,精确确定火成岩的顶界和底界,层位的解释精度与火成岩建模的精度密切相关,进而将火成岩的底界和定界相减,获得火成岩的层间双程旅行时;
将均方根速度通过迪克斯(Dix)公式转换为初始层速度,获得初始层速度模型,作为整个深度偏移建模的基础。
步骤2:针对初始层速度模型进行目标线偏移,优化火成岩之上的地层层速度,对火成岩和火成岩之下的层位进行约束,保持其层速度不变,输出优化后的共成像点道集,优化方法可以选用比较成熟的网格层析方法;
优化后的共成像点道集是拉平的,则火成岩之上的地层速度最佳优化,将优化后的共成像点道集进行近偏移距和远偏移距叠加,获得近偏移距叠加剖面与远偏移距叠加剖面,若最大偏移距为6000米,那么3000米前为近偏移距,之后为远偏移距;根据优化后的共成像点道集与全偏移距叠加剖面,拾取火成岩深度域的初始顶界面和初始底界面;比较近偏移距叠加剖面与远偏移距叠加剖面,在近偏移距叠加剖面上沿着初始底界面上拾取深度误差,进而获得校正底界面
根据优化后的共成像点道集与全偏移距叠加剖面,拾取火成岩深度域的初始顶界面和初始底界面;
地震波的传播可以近似为射线路径传播,由于时间域处理过程中假设地层为水平层状,因此局部速度异常区域在不同偏移距的叠加剖面上会有不同的反映,通过比较近偏移距叠加剖面和远偏移距叠加剖面,如果存在速度异常,速度异常反映在深度剖面上则火成岩的底界存在深度误差,则在近偏移距叠加剖面上沿着初始底界面拾取深度误差,并将拾取的深度误差加到初始底界面上来校正火成岩底界面,进而获得校正底界面;
根据校正底界面与初始顶界面,从初始顶界面的各个网格点出发,求取与校正底界面的垂直射线,作为地震波的传播路径,与层间双程旅行时做运算可以求取地震波在火成岩内的传播速度,获得校正层速度。
步骤3:根据校正底界面与校正层速度,获得校正层速度模型。
步骤4:以校正层速度模型替换初始层速度模型,重复步骤1-4,直至火成岩底界地层的反射地震同相轴拉平,获得最终的校正层速度模型。
步骤5:获得最终的校正层速度模型后,可以继续进行下面地层的速度建模工作,可以采用常规的网格层析速度建模方法继续进行速度建模,迭代过程中不能再更新火成岩及其以上地层的层速度,最终完成整个工区的深度域速度建模工作,获得火成岩叠前深度偏移成像速度模型。
本方法通过建立准确的深度域火成岩层速度模型,对于火成岩多期发育地区的地震数据能够进行准确的地震叠前深度偏移成像,消除之前由于火成岩速度不准在地震剖面上引起的假断裂以及地震同相轴上提或下拉的异常。
应用示例
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
根据本发明的多期发育火成岩叠前深度偏移成像速度建模方法可以包括:
步骤1:针对地震数据进行时间域处理,计算地下地层的均方根速度,获得时间域偏移叠加剖面;根据时间域偏移叠加剖面,确定火成岩的顶界和底界,进而将火成岩的底界和定界相减,获得火成岩的层间双程旅行时;将均方根速度转换为初始层速度,获得初始层速度模型;
步骤2:针对初始层速度模型进行目标线偏移,优化火成岩之上的地层层速度,输出优化后的共成像点道集;将优化后的共成像点道集进行近偏移距和远偏移距叠加,获得近偏移距叠加剖面;在近偏移距叠加剖面拾取火成岩深度域的初始顶界面和初始底界面;在初始底界面上拾取深度误差,进而获得校正底界面;根据校正底界面与初始顶界面,获得校正层速度;
步骤3:根据校正底界面与校正层速度,获得校正层速度模型;
步骤4:以校正层速度模型替换初始层速度模型,重复步骤1-4,直至火成岩底界地层的反射地震同相轴拉平,获得最终的校正层速度模型;
步骤5:根据最终的校正层速度模型,获得火成岩叠前深度偏移成像速度模型。
对某火成岩多期发育地区的实际三维地震数据进行利用本方法进行了深度速度建模和成像。该地区在在二叠系广泛发育火成岩,在火成岩的喷发和溢流过程中形成了各种岩性的火成岩,包括玄武岩、英安岩、凝灰岩以等。在对该地区采集的地震数据的时间域处理过程中,发现这些岩石混在一起产生了大量的横向速度异常,严重影响了下部的时间域地震成像,导致下部地层时间域地震信号信噪比降低,同相轴也发生了扭曲变形,最严重的地方出现了垂直的假断层。为了消除这些异常显示,首先对该地区的地震数据进行了时间域处理得到均方根速度模型,在此基础上利用迪克斯公式转化为深度域初始层速度模型。
图2a和图2b分别示出了根据本发明的一个实施例的应用本方法前后的叠前深度偏移成像结果的示意图。
图2a是没有使用本方法对火成岩建模之前的偏移结果,由于火成岩内部速度异常的影响因此底界有明显的波浪状起伏,而且底界的下伏地层反射杂乱,看不出有效的反射。通过本方法对火成岩的层速度进行精确建模之后,图2b消除了底界的波浪状起伏,同时下部的恢复了真实的层状反射构造,消除了之前的反射假象,方框下部的反射同相轴更为平直和聚焦。
综上所述,本发明通过建立准确的深度域火成岩层速度模型,对于火成岩多期发育地区的地震数据能够进行准确的地震叠前深度偏移成像,消除之前由于火成岩速度不准在地震剖面上引起的假断裂以及地震同相轴上提或下拉的异常。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
根据本发明的实施例,提供了一种多期发育火成岩叠前深度偏移成像速度建模系统,其特征在于,该系统包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:步骤1:针对地震数据,获得初始层速度模型;步骤2:针对初始层速度模型进行校正,获得校正底界面与校正层速度;步骤3:根据校正底界面与校正层速度,获得校正层速度模型;步骤4:以校正层速度模型替换初始层速度模型,重复步骤1-4,直至火成岩底界地层的反射地震同相轴拉平,获得最终的校正层速度模型;步骤5:根据最终的校正层速度模型,获得火成岩叠前深度偏移成像速度模型。
在一个示例中,步骤1包括:针对地震数据进行时间域处理,计算地下地层的均方根速度,获得时间域偏移叠加剖面;根据时间域偏移叠加剖面,确定火成岩的顶界和底界,进而计算层间双程旅行时;将均方根速度转换为初始层速度,获得初始层速度模型。
在一个示例中,将火成岩的底界和定界相减,获得火成岩的层间双程旅行时。
在一个示例中,步骤2包括:针对初始层速度模型进行目标线偏移,优化火成岩之上的地层层速度,输出优化后的共成像点道集;根据优化后的共成像点道集,拾取火成岩深度域的初始顶界面和初始底界面;在初始底界面上拾取深度误差,进而获得校正底界面;根据校正底界面与初始顶界面,获得校正层速度。
在一个示例中,根据优化后的共成像点道集,拾取火成岩深度域的初始顶界面和初始底界面,在初始底界面上拾取深度误差,进而获得校正底界面包括:将优化后的共成像点道集进行近偏移距和远偏移距叠加,获得近偏移距叠加剖面与远偏移距叠加剖面;根据优化后的共成像点道集与全偏移距叠加剖面,拾取火成岩深度域的初始顶界面和初始底界面;比较近偏移距叠加剖面与远偏移距叠加剖面,在近偏移距叠加剖面上沿着初始底界面上拾取深度误差,进而获得校正底界面。
本系统通过建立准确的深度域火成岩层速度模型,对于火成岩多期发育地区的地震数据能够进行准确的地震叠前深度偏移成像,消除之前由于火成岩速度不准在地震剖面上引起的假断裂以及地震同相轴上提或下拉的异常。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种多期发育火成岩叠前深度偏移成像速度建模方法,其特征在于,包括:
步骤1:针对地震数据,获得初始层速度模型;
步骤2:针对所述初始层速度模型进行校正,获得校正底界面与校正层速度;
步骤3:根据所述校正底界面与所述校正层速度,获得校正层速度模型;
步骤4:以所述校正层速度模型替换所述初始层速度模型,重复步骤1-4,直至火成岩底界地层的反射地震同相轴拉平,获得最终的校正层速度模型;
步骤5:根据所述最终的校正层速度模型,获得火成岩叠前深度偏移成像速度模型。
2.根据权利要求1所述的多期发育火成岩叠前深度偏移成像速度建模方法,其中,所述步骤1包括:
针对地震数据进行时间域处理,计算地下地层的均方根速度,获得时间域偏移叠加剖面;
根据所述时间域偏移叠加剖面,确定火成岩的顶界和底界,进而计算层间双程旅行时;
将所述均方根速度转换为初始层速度,获得初始层速度模型。
3.根据权利要求2所述的多期发育火成岩叠前深度偏移成像速度建模方法,其中,将火成岩的底界和定界相减,获得火成岩的层间双程旅行时。
4.根据权利要求1所述的多期发育火成岩叠前深度偏移成像速度建模方法,其中,所述步骤2包括:
针对所述初始层速度模型进行目标线偏移,优化火成岩之上的地层层速度,输出优化后的共成像点道集;
根据所述优化后的共成像点道集,拾取火成岩深度域的初始顶界面和初始底界面,在所述初始底界面上拾取深度误差,进而获得校正底界面;
根据所述校正底界面与所述初始顶界面,获得校正层速度。
5.根据权利要求4所述的多期发育火成岩叠前深度偏移成像速度建模方法,其中,所述根据所述优化后的共成像点道集,拾取火成岩深度域的初始顶界面和初始底界面,在所述初始底界面上拾取深度误差,进而获得校正底界面包括:
将所述优化后的共成像点道集进行近偏移距和远偏移距叠加,获得近偏移距叠加剖面与远偏移距叠加剖面;
根据所述优化后的共成像点道集与全偏移距叠加剖面,拾取火成岩深度域的初始顶界面和初始底界面;
比较所述近偏移距叠加剖面与所述远偏移距叠加剖面,在所述近偏移距叠加剖面上沿着所述初始底界面上拾取深度误差,进而获得校正底界面。
6.一种多期发育火成岩叠前深度偏移成像速度建模系统,其特征在于,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
步骤1:针对地震数据,获得初始层速度模型;
步骤2:针对所述初始层速度模型进行校正,获得校正底界面与校正层速度;
步骤3:根据所述校正底界面与所述校正层速度,获得校正层速度模型;
步骤4:以所述校正层速度模型替换所述初始层速度模型,重复步骤1-4,直至火成岩底界地层的反射地震同相轴拉平,获得最终的校正层速度模型;
步骤5:根据所述最终的校正层速度模型,获得火成岩叠前深度偏移成像速度模型。
7.根据权利要求6所述的多期发育火成岩叠前深度偏移成像速度建模系统,其中,所述步骤1包括:
针对地震数据进行时间域处理,计算地下地层的均方根速度,获得时间域偏移叠加剖面;
根据所述时间域偏移叠加剖面,确定火成岩的顶界和底界,进而计算层间双程旅行时;
将所述均方根速度转换为初始层速度,获得初始层速度模型。
8.根据权利要求7所述的多期发育火成岩叠前深度偏移成像速度建模系统,其中,将火成岩的底界和定界相减,获得火成岩的层间双程旅行时。
9.根据权利要求6所述的多期发育火成岩叠前深度偏移成像速度建模系统,其中,所述步骤2包括:
针对所述初始层速度模型进行目标线偏移,优化火成岩之上的地层层速度,输出优化后的共成像点道集;
根据所述优化后的共成像点道集,拾取火成岩深度域的初始顶界面和初始底界面,在所述初始底界面上拾取深度误差,进而获得校正底界面;
根据所述校正底界面与所述初始顶界面,获得校正层速度。
10.根据权利要求9所述的多期发育火成岩叠前深度偏移成像速度建模系统,其中,所述根据所述优化后的共成像点道集,拾取火成岩深度域的初始顶界面和初始底界面,在所述初始底界面上拾取深度误差,进而获得校正底界面包括:
将所述优化后的共成像点道集进行近偏移距和远偏移距叠加,获得近偏移距叠加剖面与远偏移距叠加剖面;
根据所述优化后的共成像点道集与全偏移距叠加剖面,拾取火成岩深度域的初始顶界面和初始底界面;
比较所述近偏移距叠加剖面与所述远偏移距叠加剖面,在所述近偏移距叠加剖面上沿着所述初始底界面上拾取深度误差,进而获得校正底界面。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910577201.0A CN112147686B (zh) | 2019-06-28 | 2019-06-28 | 多期发育火成岩叠前深度偏移成像速度建模方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910577201.0A CN112147686B (zh) | 2019-06-28 | 2019-06-28 | 多期发育火成岩叠前深度偏移成像速度建模方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112147686A true CN112147686A (zh) | 2020-12-29 |
CN112147686B CN112147686B (zh) | 2022-06-21 |
Family
ID=73869542
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910577201.0A Active CN112147686B (zh) | 2019-06-28 | 2019-06-28 | 多期发育火成岩叠前深度偏移成像速度建模方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112147686B (zh) |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050237857A1 (en) * | 2004-04-22 | 2005-10-27 | Junru Jiao | Horizon-based residual depth migration velocity analysis |
US20100074053A1 (en) * | 2008-07-18 | 2010-03-25 | William Marsh Rice University | Methods for concurrent generation of velocity models and depth images from seismic data |
CN102279416A (zh) * | 2010-12-30 | 2011-12-14 | 北京诺克斯达石油科技有限公司 | 远近偏移距叠加剖面相位频率的匹配方法 |
CN104375172A (zh) * | 2013-08-15 | 2015-02-25 | 中国石油集团东方地球物理勘探有限责任公司 | 火山岩下伏地层的构造形态正确成像的方法 |
CN104950330A (zh) * | 2014-03-25 | 2015-09-30 | 中国石油化工股份有限公司 | 气顶油藏深度域成像的速度建模方法及系统 |
CN108072892A (zh) * | 2016-11-09 | 2018-05-25 | 中国石油化工股份有限公司 | 一种自动化的地质构造约束层析反演方法 |
CN109738945A (zh) * | 2018-11-08 | 2019-05-10 | 成都捷科思石油天然气技术发展有限公司 | 一种利用叠前深度偏移成果直接生成构造图的方法 |
CN109839660A (zh) * | 2018-11-08 | 2019-06-04 | 成都捷科思石油天然气技术发展有限公司 | 一种利用叠前道集数据建立速度深度模型的方法 |
CN109884700A (zh) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | 多信息融合地震速度建模方法 |
-
2019
- 2019-06-28 CN CN201910577201.0A patent/CN112147686B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20050237857A1 (en) * | 2004-04-22 | 2005-10-27 | Junru Jiao | Horizon-based residual depth migration velocity analysis |
US20100074053A1 (en) * | 2008-07-18 | 2010-03-25 | William Marsh Rice University | Methods for concurrent generation of velocity models and depth images from seismic data |
CN102279416A (zh) * | 2010-12-30 | 2011-12-14 | 北京诺克斯达石油科技有限公司 | 远近偏移距叠加剖面相位频率的匹配方法 |
CN104375172A (zh) * | 2013-08-15 | 2015-02-25 | 中国石油集团东方地球物理勘探有限责任公司 | 火山岩下伏地层的构造形态正确成像的方法 |
CN104950330A (zh) * | 2014-03-25 | 2015-09-30 | 中国石油化工股份有限公司 | 气顶油藏深度域成像的速度建模方法及系统 |
CN108072892A (zh) * | 2016-11-09 | 2018-05-25 | 中国石油化工股份有限公司 | 一种自动化的地质构造约束层析反演方法 |
CN109738945A (zh) * | 2018-11-08 | 2019-05-10 | 成都捷科思石油天然气技术发展有限公司 | 一种利用叠前深度偏移成果直接生成构造图的方法 |
CN109839660A (zh) * | 2018-11-08 | 2019-06-04 | 成都捷科思石油天然气技术发展有限公司 | 一种利用叠前道集数据建立速度深度模型的方法 |
CN109884700A (zh) * | 2019-03-20 | 2019-06-14 | 中国石油化工股份有限公司 | 多信息融合地震速度建模方法 |
Non-Patent Citations (6)
Title |
---|
曾庆才,等: "盐下高陡构造成像技术-以塔里木盆地库车坳陷克深地区为例", 《石油勘探与开发》 * |
曾庆才,等: "盐下高陡构造成像技术-以塔里木盆地库车坳陷克深地区为例", 《石油勘探与开发》, 31 December 2017 (2017-12-31), pages 871 - 879 * |
温铁民,等: "基于地质模型的火山岩体速度建模方法", 《CPS/SEG北京2018国际地球物理会议暨展览电子论文集》 * |
温铁民,等: "基于地质模型的火山岩体速度建模方法", 《CPS/SEG北京2018国际地球物理会议暨展览电子论文集》, 31 December 2018 (2018-12-31), pages 294 - 296 * |
黄棱,等: "速度优化与建模技术在德惠断陷火山岩成像中的应用", 《石油地球物理勘探》 * |
黄棱,等: "速度优化与建模技术在德惠断陷火山岩成像中的应用", 《石油地球物理勘探》, 30 November 2017 (2017-11-30), pages 34 - 40 * |
Also Published As
Publication number | Publication date |
---|---|
CN112147686B (zh) | 2022-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3011366B1 (en) | Seismic data processing | |
EP3069171B1 (en) | Seismic data processing | |
US20100118653A1 (en) | Vertical seismic profiling velocity estimation method | |
CN102147481B (zh) | 用于三维表面相关多次波预测中的数据重建的基于倾角的校正 | |
WO2017035104A1 (en) | Velocity model seismic static correction | |
CN104991268B (zh) | 一种真振幅偏移成像方法 | |
CN111722284B (zh) | 一种基于道集数据建立速度深度模型的方法 | |
CN111123359B (zh) | 随钻测井与地层格架约束的井周地震成像探测方法及装置 | |
CN102879816A (zh) | 一种地震多次波偏移方法 | |
Wang et al. | Multiparameter TTI tomography of P-wave reflection and VSP data | |
Witten et al. | Microseismic image-domain velocity inversion: Marcellus Shale case study | |
WO2015058015A1 (en) | Iterative stacking of seismic image partitions | |
Gras et al. | Full-waveform inversion of short-offset, band-limited seismic data in the Alboran Basin (SE Iberia) | |
CN113805237A (zh) | 偏移陆地交叉排列地震的使用压缩感测模型的方法和系统 | |
CN104536041B (zh) | 一种地震观测系统参数的优化方法 | |
CA2926179C (en) | System and method for seismic adaptive optics | |
CN102565852B (zh) | 针对储层含油气性检测的角度域叠前偏移数据处理方法 | |
AU2011279350B2 (en) | Method for accentuating specular and non-specular seismic events from within shallow subsurface rock formations | |
WO2012054124A1 (en) | System and method for characterization with non-unique solutions of anisotropic velocities | |
CN105425300B (zh) | 一种剩余静校正方法 | |
CN112147686B (zh) | 多期发育火成岩叠前深度偏移成像速度建模方法及系统 | |
Nosjean et al. | Statics: from imaging to interpretation pitfalls and an efficient way to overcome them | |
CN111142159B (zh) | 中深层异常体地震速度建模方法及系统 | |
CN104237939B (zh) | 一种三维地震观测系统多频率偏移噪声获取方法及装置 | |
CN110941029B (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 |