CN105510967A - 一种实现三维地震数据偏移归位的处理方法、装置 - Google Patents
一种实现三维地震数据偏移归位的处理方法、装置 Download PDFInfo
- Publication number
- CN105510967A CN105510967A CN201511022795.7A CN201511022795A CN105510967A CN 105510967 A CN105510967 A CN 105510967A CN 201511022795 A CN201511022795 A CN 201511022795A CN 105510967 A CN105510967 A CN 105510967A
- Authority
- CN
- China
- Prior art keywords
- data
- migration
- time
- seismic data
- 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.)
- Granted
Links
- 238000013508 migration Methods 0.000 title claims abstract description 105
- 230000005012 migration Effects 0.000 title claims abstract description 105
- 238000003672 processing method Methods 0.000 title abstract 2
- 238000000034 method Methods 0.000 claims abstract description 70
- 230000008569 process Effects 0.000 claims abstract description 32
- 238000012545 processing Methods 0.000 claims abstract description 32
- 238000012937 correction Methods 0.000 claims abstract description 11
- 230000005284 excitation Effects 0.000 claims description 50
- 230000015572 biosynthetic process Effects 0.000 claims description 25
- 238000003786 synthesis reaction Methods 0.000 claims description 24
- 229940054107 isochron Drugs 0.000 claims description 21
- MOYKHGMNXAOIAT-JGWLITMVSA-N isosorbide dinitrate Chemical compound [O-][N+](=O)O[C@H]1CO[C@@H]2[C@H](O[N+](=O)[O-])CO[C@@H]21 MOYKHGMNXAOIAT-JGWLITMVSA-N 0.000 claims description 21
- 238000001914 filtration Methods 0.000 claims description 18
- 238000012546 transfer Methods 0.000 claims description 10
- 230000000052 comparative effect Effects 0.000 claims description 9
- 230000008859 change Effects 0.000 abstract description 18
- 238000006243 chemical reaction Methods 0.000 abstract description 2
- 239000011229 interlayer Substances 0.000 abstract 1
- 230000002194 synthesizing effect Effects 0.000 abstract 1
- 238000010586 diagram Methods 0.000 description 18
- 238000013507 mapping Methods 0.000 description 8
- 230000006870 function Effects 0.000 description 7
- 238000003325 tomography Methods 0.000 description 6
- 239000003245 coal Substances 0.000 description 4
- 238000005070 sampling Methods 0.000 description 4
- 229910052500 inorganic mineral Inorganic materials 0.000 description 3
- 239000011707 mineral Substances 0.000 description 3
- 238000003384 imaging method Methods 0.000 description 2
- 238000004519 manufacturing process Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- 230000002123 temporal effect Effects 0.000 description 2
- 241001269238 Data Species 0.000 description 1
- 235000014676 Phragmites communis Nutrition 0.000 description 1
- 230000002159 abnormal effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 239000007795 chemical reaction product Substances 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 238000009826 distribution Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000008676 import Effects 0.000 description 1
- 238000005065 mining Methods 0.000 description 1
- 239000000047 product Substances 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
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/301—Analysis for determining seismic cross-sections or geostructures
- G01V1/302—Analysis for determining seismic cross-sections or geostructures in 3D data cubes
-
- 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/34—Displaying seismic recordings or visualisation of seismic data or attributes
- G01V1/345—Visualisation of seismic data or attributes, e.g. in 3D cubes
-
- 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/64—Geostructures, e.g. in 3D data cubes
-
- 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
- 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
技术领域
本发明属于地震勘探数据处理技术领域,尤其涉及一种实现三维地震数据偏移归位的处理方法、装置。
背景技术
众所周至,本领域人员在地震等时平面图偏移归位处理中,长期以来使用的主要是图偏移方法,常用的方法有t0梯度法及直接归位法等。这些方法均是直接对图纸进行偏移归位,解释人员利用这些方法工作时的劳动强度大、工作效率低,成图精度也不高。
一般的,所述的偏移归位通常是把水平叠加剖面上偏移的反射层进行反偏移,使地层的真实位置形态得到恢复,保障地震数据资料准确解释、分析。当地层水平时,上述方法显示出的地层构造通常与实际情况比较相符。而对于多层目的层倾角、倾向变化大的情况下的地震数据处理进行偏移归位时,界面法线平面与铅垂面并不正交,地层角越大,两者的差别也越大,往往容易产生较大的误差或者错误,结果导致绘制出的构造图误差较大,成图精度低或出现异常,迫使作业人员重新核查计算,劳动强度更大、工作效率。复杂断层、倾角变化大等会对地震资料解释的产生较大大的影响,使得时间剖面上记录点位置与反射点的位置不相符合,记录点的显示位置总数相对于反射点向界面的下倾方向移动,不利于地震资料的地震解释。
现有技术中对地震数据进行偏移归位的处理工作时的劳动强度大,成图精度也不高。尤其是对于多层目的层倾角和倾向变化大的复杂地震数据会还会存在着进行资料解释的工作效率和解释精度均较低的问题。
发明内容
本发明目的在于提供一种实现三维地震数据偏移归位的处理方法、装置,可以提高工作效率和成图精度,同时实现大倾角区的二维地震资料的三维偏移归位,克服复杂断层、倾角变化大对地震资料解释的影响,解决多目的层间倾角、倾向变化大条件下高精度归位的问题。
本发明提供的一种实现三维地震数据偏移归位的处理方法、装置是这样实现的:
一种实现三维地震数据偏移归位处理方法,所述方法包括:
将采集获取的二维地震数据合成三维地震时间数据体;
将所述三维地震时间数据体处理为自激自收三维地震资料数据体;
利用对二维地震资料处理中的倾角时差校正速度进行插值得到三维地震速度对所述自激自收三维地震资料数据体进行叠后时间偏移,得到三维叠后时间偏移数据体;
对三维叠后时间偏移数据体进行时深转换,生成偏移归位后三维地震数据。
可选的实施例中,所述将采集获取的二维地震数据合成三维地震时间数据体包括:
对采集获取的二维地震数据进行处理得到地震叠加时间剖面,对所述地震叠加时间剖面上的目的层进行对比解释,获取二维地震解释t0等时图;
将各个目的层的二维地震解释t0等时图进行网格化匹配计算得到各层双程旅行时脉冲反射系数的三维地震时间数据体。
可选的实施例中,所述将所述三维地震时间数据体处理为自激自收三维地震资料数据体包括:
以频带为带限的雷克子波与三维地震时间数据体进行褶积计算得到自激自收三维地震资料数据体,所述雷克子波的频带宽度可以设置为10Hz至120Hz。
可选的实施例中,所述将所述三维地震时间数据体处理为自激自收三维地震资料数据体包括:
对所述三维三维地震时间数据体进行带通滤波得到自激自收三维地震资料数据体,所述带通滤波的频带宽度可以设置为10Hz至120Hz。
可选的实施例中,所述方法还包括:
利用所述偏移归位后三维地震数据生成用于地震资料解释的构造图。
一种实现三维地震数据偏移归位的处理装置,所述装置包括:
三维数据合成模块,用于将采集获取的二维地震数据合成三维地震时间数据体;
自激自收处理模块,用于将所述三维地震时间数据体处理为自激自收三维地震资料数据体;
偏移处理模块,用于利用对二维地震资料处理中的倾角时差校正速度进行插值得到三维地震速度对所述自激自收三维地震资料数据体进行叠后时间偏移,得到三维叠后时间偏移数据体;
时深转换模块,用于对三维叠后时间偏移数据体进行时深转换,生成偏移归位后三维地震数据。
可选的实施例中,所述三维数据合成模块以包括:
叠加时间剖面模块,用于对采集获取的二维地震数据进行处理得到地震叠加时间剖面;
等时图处理模块,用于对所述地震叠加时间剖面上的目的层进行对比解释处理,获取二维地震解释t0等时图;
网格化匹配计算模块,用于将各个目的层的二维地震解释t0等时图进行网格化匹配计算得到各层双程旅行时脉冲反射系数的三维地震时间数据体。
可选的实施例中,所述自激自收处理模块包括:
子波计算模块,用于生成预设频带宽度的雷克子波,所述预设频带宽度设置为10Hz至120Hz;
褶积计算模块,用于以所述生成的雷克子波与三维地震时间数据体进行褶积计算得到自激自收三维地震资料数据体。
可选的实施例中,所述自激自收处理模块包括:
带通滤波处理模块,用于对所述三维地震时间数据体进行带通滤波得到自激自收三维地震资料数据体,所述带通滤波的频带宽度设置为10Hz至120Hz。
可选的实施例中,所述装置还包括:
构造图生成模块,用于利用所述偏移归位后三维地震数据生成用于地震资料解释的构造图。
本发明提供的一种实现三维地震数据偏移归位的处理方法、装置,可以利用二维地震资料进行三维偏移归位的处理方法,能完成二维地震资料处理为t0等时图,并完成二维地震解释t0等时图的三维偏移归位。采用本发明所述的方法可以实现大倾角区的二维地震资料的三维偏移归位,大大减小复杂断层、倾角变化大对地震资料解释的影响,提高地震解释构造图的成图精度,进而减小解释人员的劳动强度,提高工作效率。本发明提供的实现三维地震数据偏移归位处理方法可以有效解决多目的层间倾角、倾向变化大条件下高精度归位的问题,并且适用范围广,不受断层复杂、倾角倾向变化大的限制,能提高评价矿藏的准确性。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是本发明一种实现三维地震数据偏移归位的处理方法一种实施例的方法流程示意图;
图2是经过处理得到的三个层位的三维叠后时间偏移数据体的立体显示示意图;
图3是本发明一种实现三维地震数据偏移归位的处理方法另一种实施例的方法流程示意图
图4是利用本发明提供的方法完成的二维地震解释构造示意图;
图5是本发明提供的实现三维地震数据偏移归位的处理装置一种实施例的模块结构示意图;
图6是本发明提供的三维数据合成模块一种实施例的模块结构示意图;
图7是本发明提供的自激自收处理模块一种实施例的模块结构示意图;
图8是本发明所述装置另一种实施例的模块结构示意图。
具体实施方式
为了使本技术领域的人员更好地理解本发明中的技术方案,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本发明保护的范围。
下面结合附图对本发明所述的一种实现三维地震数据偏移归位的处理方法进行详细的说明。图1是本发明所述一种实现三维地震数据偏移归位的处理方法一种实施例的方法流程图。虽然本发明提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或终端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
具体的如图1所示,本发明提供的一种实现三维地震数据偏移归位的处理方法的一种实施例中,所述方法可以包括:
S1:将采集获取的二维地震数据合成三维地震时间数据体T(x,y,t)。
本发明提供的方法可以利用采集获取的石油、煤炭地震勘探或其他进行地震资料解释的人工二维地震数据处理成t0等时图,完成二维地震解释t0等时图的三维偏移归位处理。
本发明中所述的合成的三维地震时间数据体通常可以包基于获取的二维地震数据资料经过在地震系统中如GeoFrame和GeovecteurPlus解释系统处理后得到的一个或多个目的层的三维地震时间数据体T(x,y,t)。具体的利用所述地震系统将二维地震数据合成三维地震时间数据体的方式可以根据数据的处理需求设定具体的处理方式,一般的,可以主要包括完成二维地震数据资料的处理及目的层的对比追踪解释,可以包括断层、褶皱等地质构造的组合。然后可以将二维地震解释的t0等时图进行网格化与网格匹配计算得到各层双程旅行时脉冲反射系统的三维地震时间数据体。
本发明可以以宁夏回族自治区苇州矿区苇四井田为应用场景来解释说明本发明实现利用二维地震数据资料进行三维偏移归位的方法。在本应用场景中的勘探区域为煤炭详查,设置的主侧线距为500m,联络线距为1000m。全区共布置主测线20条。本工区构造复杂,为了解决高精度成图,探明地下煤层分布的问题而实际应用的实施例。在发明中,合成三维地震时间数据体T(x,y,t)的处理过程中,三维地震时间数据体T(x,y,t)的时间采样率与记录长度可以设置为与二维地震处理叠加剖面一致。本实施例中的时间采样率为1ms,记录长度为2s。本应用场景中工区为古生界地层,速度为3500m/s,地层倾角20°,频带宽度为(10,120Hz),经过计算不产生空间假频的CDP尺寸小于21.33m,因此本实施例的空间采样率选择20m*20m。当然根据不同的处理数据,所述的空间采样率可以选取10m*10m或20m*20m,具体的可以可以后获取的时间采集数据和处理需求进行合理设置。
本发明提供一种具体的利用GeoFrame和GeovecteurPlus解释系统来合成所述三维地震时间数据体的处理方法。具体的,本发明的一种实施例中,所述将采集获取的二维地震数据合成三维地震时间数据体可以包括:
S101:对采集获取的野外二维地震数据进行室内资料处理得到时间域地震叠加剖面,对所述地震叠加时间剖面上的目的层进行对比解释,获取二维地震解释t0等时图;
本实施例中所述的室内资料处理通常是指对野外采集的地震资料进行去噪、静校正、速度分析、动校正、切除、叠加等数据处理手段,得到时间域叠加地震剖面。所述的T0等时图通常是是指在处理剖面中某地层界面的自激自收时间等值线图。所述的对比解释通常是指地震记录上利用反射波的动力学特点和运动学特点来识别和追踪同一界面的波的工作。
S102:将各个目的层的二维地震解释t0等时图进行网格化匹配计算得到各层双程旅行时脉冲反射系数的三维地震时间数据体。
本实施例中所述的网格匹配计算通常是指把二维t0时间,按照指定的空间网格大小进行三维插值,以得到(x,y)平面内各网格节点处的t0时间T(x,y,t)。
在实际生产工作中,用于资料解释的通常是由多个地震道依次排序起来形成的地震剖面或地震数据体,不同类型和传播特点的波的同相轴会在地震剖面上表现出不同的特点,这些特点通常是作业人员进行地震资料解释时在地震剖面或数据体上识别各种波的主要依据。在本发明方法中,首先,可以将实际采集得到的二维地震数据(如数字地震磁带记录)在地震处理系统上进行处理,得到地震叠加时间剖面。本实施例中可以在GeoFrame地震解释系统上完成地震叠加时间剖面上目的层(4煤和12煤层)的对比追踪,具体的可以利用测井资料(人工伽玛或声波时差)制作合成地震记录。然后可以对各个层位进行标定,如本应用场景中4煤对应的地震波可以标定为T4波,12煤对应的地震波可以标定为T12波,新生界底界面对应的地震波可以标定为T0波。本实施例中可以根据钻孔标定这三层反射波,对断层、褶积等地震构造进行了组合。
进一步的,可以将T0、T4、T12等三层的解释层位资料和断层资料传入CPS-3地质绘图系统,利用CPS-3地质绘图软件将T0、T4、T12等三层层位的二维地震解释的t0进行网格化与网格匹配,并将他传入GeovecteurPlus二维地震资料系统中。在GeovecteurPlus二维地震资料系统中计算得到T0、T4、T12的三维地震时间数据体(通常可以包括,首先得到一个平界面的三维地震时间数据体,然后利用二维地震解释的t0进行静校正得到任意倾斜界面的三维地震时间数据体)。反射界面处以脉冲反射系数表示(以一个常数表示),其它区域为零。最后可以将T0、T4、T12的单个三维地震时间数据体加在一起形成一个三层的三维地震时间数据体T(x,y,t)。
本发明中可以将采集获取的二维地震数据合成三维地震时间数据体T(x,y,t)。
S2:将所述三维地震时间数据体处理为自激自收三维地震资料数据体Astk(x,y,t)。
本发明方法在获取合成的三维地震时间数据体T(x,y,t)后可以将其处理成自激自收三维地震资料数据体Astk(x,y,t)。一般的,自激自收时间通常是指是指在地表激发以垂直角度入射在反射界面,再从反射界面垂直反射到激发点的双程旅行时间,它可以直接反应地层深度的变化趋势。这里所指的自激自收三维地震数据体是整个反射界面对应的自激自收时间在三维空间的展布数据体。
本发明提供的一种实施方式中,可以利用以频宽为带限的雷克子波或设定频带的带通滤波将所述三维地震时间数据体处理为自激自收三维地震资料数据体Astk(x,y,t)。具体的,本发明的一种实施例中,所述将所述三维地震时间数据体处理为自激自收三维地震资料数据体Astk(x,y,t)可以包括:
S201:以频带为带限的雷克子波与三维地震时间数据体T(x,y,t)进行褶积计算得到自激自收三维地震资料数据体Astk(x,y,t),所述雷克子波的频带宽度可以设置为10Hz至120Hz。
在其他的时实施方式中还可以直接进行带通滤波得到。因此,本发明的另一种实施例中,所述将所述三维地震时间数据体处理为自激自收三维地震资料数据体Astk(x,y,t)可以包括:
S202:对所述三维三维地震时间数据体T(x,y,t)进行带通滤波得到自激自收三维地震资料数据体,所述带通滤波的频带宽度可以设置为10Hz至120Hz。
当然,本发明所述的方法也包括采用其他方式将所述三维地震时间数据体处理为自激自收三维地震资料数据体的实施方式。
S3:利用对二维地震资料处理中的倾角时差校正速度进行插值得到三维地震速度对所述自激自收三维地震资料数据体进行叠后时间偏移,得到三维叠后时间偏移数据体Amig(x,y,t)。
本发明实施例中可以将在GeovecteurPlus二维地震资料系统中的倾角时差校正(DMO)速度进行插值得到三维地震速度场。利用该三维地震速度场中的速度对所述自激自收三维地震资料数据体Astk(x,y,t)进行叠后时间偏移即得到三维叠后时间偏移数据体Amig(x,y,t)。偏移到位情况可以根据井资料和其他已知资料进行偏移速度的修改,然后可以将偏移归位的三维叠后时间偏移数据体Amig(x,y,t)传入GeoFrame地震解释系统。图2是经过处理得到的三个层位的三维叠后时间偏移数据体的立体显示示意图。
S4:对三维叠后时间偏移数据体Amig(x,y,t)进行层位时深转换,生成偏移归位后三维地震数据。
进一步在GeoFrame地震解释系统上对三维叠后时间偏移数据体Amig(x,y,t)进行层位的对比解释、追踪与时深转换,归位即完成。
当然,进一步的,利用二维地震资料完成三维偏移归位后,可以基于本发明方法偏移归位处理后的三维地震数据生成二维或三维用于地震解释的构造图。因此,本发明所述方法还可以包括:
S5:利用所述偏移归位后三维地震数据生成用于地震资料解释的构造图。
图3是本发明所述方法另一种实施例的方法流程示意图。这样,由于使用本发明方法偏移处理后的三维地震数据克服了复杂断层、倾角变化大对地震资料解释的影响,解决多目的层间倾角、倾向变化大条件下高精度归位的问题,因此生成的地震构造解释图成图精度更高,可以更加准确合理的进行资料解释,为地质勘探提供更加有力的依据。图4是利用本发明提供的方法完成的二维地震解释构造示意图。
本发明提供的利用二维地震资料进行三维偏移归位的处理方法,能完成二维地震资料处理为t0等时图,并完成二维地震解释t0等时图的三维偏移归位。采用本发明所述的方法可以实现大倾角区的二维地震资料的三维偏移归位,避免了二维偏移在陡倾角和侧面干扰难以正确成像的问题,大大减小复杂断层、倾角变化大对地震资料解释的影响,提高地震解释构造图的成图精度,进而减小解释人员的劳动强度,提高工作效率。本发明提供的实现三维地震数据偏移归位处理方法可以有效解决多目的层间倾角、倾向变化大条件下高精度归位的问题,并且适用范围广,不受断层复杂、倾角倾向变化大的限制,能提高评价矿藏的准确性。
本发明提供的实现三维地震数据偏移归位的处理方法可以在多种地震资料解释系统或地震数据处理系统上实现,进而提高系统的三维地震数据偏移归位处理的精度,提高地震解释构造图的成图质量。因此,基于本发明所述的实现三维地震数据偏移归位的处理方法,本发明还提供一种实现三维地震数据偏移归位的处理装置。图5是本发明所述装置一种实施例的模块结构示意图,如图5所示,所述装置可以包括:
三维数据合成模块101,可以用于将采集获取的二维地震数据合成三维地震时间数据体;
自激自收处理模块102,可以用于将所述三维地震时间数据体处理为自激自收三维地震资料数据体;
偏移处理模块103,可以用于利用对二维地震资料处理中的倾角时差校正速度进行插值得到三维地震速度对所述自激自收三维地震资料数据体进行叠后时间偏移,得到三维叠后时间偏移数据体;
时深转换模块104,可以用于对三维叠后时间偏移数据体进行时深转换,生成偏移归位后三维地震数据。
本发明的一种实施方式中,利用二维地震解释t0图合成三维地震时间数据体。图6是本发明所述三维数据合成模块101一种实施例的模块结构示意图,如图5所示,所述三维数据合成模块101可以包括:
叠加时间剖面模块1011,可以用于对采集获取的二维地震数据进行处理得到地震叠加时间剖面;
等时图处理模块1012,可以用于对所述地震叠加时间剖面上的目的层进行对比解释处理,获取二维地震解释t0等时图;
网格化匹配计算模块1013,可以用于将各个目的层的二维地震解释t0等时图进行网格化匹配计算得到各层双程旅行时脉冲反射系数的三维地震时间数据体。
图7是本发明提供的所述自激自收处理模块102一种实施例的模块结构示意图,如图7所示,本发明的另一种实施例中,所述自激自收处理模块102可以包括:
子波计算模块1021,可以用于生成预设频带宽度的雷克子波,所述预设频带宽度可以设置为10Hz至120Hz;
褶积计算模块1022,可以用于以所述生成的雷克子波与三维地震时间数据体进行褶积计算得到自激自收三维地震资料数据体。
当然,如前所述,在其他的实施例中也可以直接对与三维地震时间数据体进行带通滤波得到自激自收三维地震资料数据体。因此,如图7所示,本发明所述装置的另一种实施例中,所述自激自收处理模块102可以包括:
带通滤波处理模块1023,可以用于对所述三维地震时间数据体进行带通滤波得到自激自收三维地震资料数据体,所述带通滤波的频带宽度可以设置为10Hz至120Hz。
进一步的,利用本发明装置进行二维地震资料的三维偏移归位后,可以基于所示偏移归位处理后的三维地震数据生成二维或三维用于地震解释的构造图。图8是本发明所述装置另一种实施例的模块结构示意图,如图8所示,所述装置还可以包括:
构造图生成模块105,可以用于利用所述偏移归位后三维地震数据生成用于地震资料解释的构造图。
采用本发明所提供的实现三维地震数据偏移归位的处理装置可以实现大倾角区的二维地震资料的三维偏移归位,大大减小复杂断层、倾角变化大对地震资料解释的影响,提高地震解释构造图的成图精度,进而减小解释人员的劳动强度,提高工作效率。本发明提供的实现三维地震数据偏移归位处理装置可以有效解决多目的层间倾角、倾向变化大条件下高精度归位的问题,并且适用范围广,不受断层复杂、倾角倾向变化大的限制,能提高评价矿藏的准确性。
尽管本发明内容中提到如GeoFrame和GeovecteurPlus解释系统、数据合成、对比解释、偏移处理等的数据采集、转换的描述,但是,本发明并不局限于必须是完全标准或者所提及的方式的处理的情况。本发明中各个实施例所涉及的上述描述仅是本发明中的一些实施例中的应用,在某些标准、方法的基础上略加修改后的处理方法也可以实行上述本发明各实施例的方案。当然,在符合本发明上述各实施例的中所述的处理方法步骤的其他无创造性的变形,仍然可以实现相同的申请,在此不再赘述。
虽然本发明提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。
上述实施例阐明的装置或模块,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本发明时可以把各模块的功能在同一个或多个软件和/或硬件中实现。当然,也可以将实现某功能的模块由多个子模块或子单元组合实现。
本发明中所述的方法、单元、模块、或装置可以以可读程序嵌入控制器方式实现,所述控制器按具体要求实现。例如,控制器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(ApplicationSpecificIntegratedCircuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC625D、AtmelAT91SAM、MicrochipPIC18F26K20以及SiliconeLabsC8051F320,存储器控制器还可以被实现为存储器的控制逻辑的一部分。本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本发明所述装置中的模块可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本发明,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本发明可借助软件加必需的硬件的方式来实现。基于这样的理解,本发明的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,也可以通过数据迁移的实施过程中体现出来。该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,终端仪表设备、服务器,或者网络设备等)执行本发明各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述,各个实施例之间相同或相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。本发明的全部或者部分可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、手持设备或便携式设备、基于微处理器的系统、可编程的电子设备、包括以上任何系统或设备的分布式计算环境等等。
虽然通过实施例描绘了本发明,本领域普通技术人员知道,本发明有许多变形和变化而不脱离本发明的精神,希望所附的权利要求包括这些变形和变化而不脱离本发明的精神。
Claims (10)
1.一种实现三维地震数据偏移归位处理方法,其特征在于,所述方法包括:
将采集获取的二维地震数据合成三维地震时间数据体;
将所述三维地震时间数据体处理为自激自收三维地震资料数据体;
利用对二维地震资料处理中的倾角时差校正速度进行插值得到三维地震速度对所述自激自收三维地震资料数据体进行叠后时间偏移,得到三维叠后时间偏移数据体;
对三维叠后时间偏移数据体进行时深转换,生成偏移归位后三维地震数据。
2.如权利要求1所述的一种实现三维地震数据偏移归位的处理方法,其特征在于,所述将采集获取的二维地震数据合成三维地震时间数据体包括:
对采集获取的二维地震数据进行处理得到地震叠加时间剖面,对所述地震叠加时间剖面上的目的层进行对比解释,获取二维地震解释t0等时图;
将各个目的层的二维地震解释t0等时图进行网格化匹配计算得到各层双程旅行时脉冲反射系数的三维地震时间数据体。
3.如权利要求1所述的一种实现三维地震数据偏移归位的处理方法,其特征在于,所述将所述三维地震时间数据体处理为自激自收三维地震资料数据体包括:
以频带为带限的雷克子波与三维地震时间数据体进行褶积计算得到自激自收三维地震资料数据体,所述雷克子波的频带宽度可以设置为10Hz至120Hz。
4.如权利要求1所述的一种实现三维地震数据偏移归位的处理方法,其特征在于,所述将所述三维地震时间数据体处理为自激自收三维地震资料数据体包括:
对所述三维三维地震时间数据体进行带通滤波得到自激自收三维地震资料数据体,所述带通滤波的频带宽度可以设置为10Hz至120Hz。
5.如权利要求1至4中任意一项所述的一种实现三维地震数据偏移归位的处理方法,其特征在于,所述方法还包括:
利用所述偏移归位后三维地震数据生成用于地震资料解释的构造图。
6.一种实现三维地震数据偏移归位的处理装置,其特征在于,所述装置包括:
三维数据合成模块,用于将采集获取的二维地震数据合成三维地震时间数据体;
自激自收处理模块,用于将所述三维地震时间数据体处理为自激自收三维地震资料数据体;
偏移处理模块,用于利用对二维地震资料处理中的倾角时差校正速度进行插值得到三维地震速度对所述自激自收三维地震资料数据体进行叠后时间偏移,得到三维叠后时间偏移数据体;
时深转换模块,用于对三维叠后时间偏移数据体进行时深转换,生成偏移归位后三维地震数据。
7.如权利要求6所述的一种实现三维地震数据偏移归位的处理装置,其特征在于,所述三维数据合成模块以包括:
叠加时间剖面模块,用于对采集获取的二维地震数据进行处理得到地震叠加时间剖面;
等时图处理模块,用于对所述地震叠加时间剖面上的目的层进行对比解释处理,获取二维地震解释t0等时图;
网格化匹配计算模块,用于将各个目的层的二维地震解释t0等时图进行网格化匹配计算得到各层双程旅行时脉冲反射系数的三维地震时间数据体。
8.如权利要求6所述的一种实现三维地震数据偏移归位的处理装置,其特征在于,所述自激自收处理模块包括:
子波计算模块,用于生成预设频带宽度的雷克子波,所述预设频带宽度设置为10Hz至120Hz;
褶积计算模块,用于以所述生成的雷克子波与三维地震时间数据体进行褶积计算得到自激自收三维地震资料数据体。
9.如权利要求6所述的一种实现三维地震数据偏移归位的处理装置,其特征在于,所述自激自收处理模块包括:
带通滤波处理模块,用于对所述三维地震时间数据体进行带通滤波得到自激自收三维地震资料数据体,所述带通滤波的频带宽度设置为10Hz至120Hz。
10.如权利要求6至9中任意一项所述的一种实现三维地震数据偏移归位的处理装置,其特征在于,所述装置还包括:
构造图生成模块,用于利用所述偏移归位后三维地震数据生成用于地震资料解释的构造图。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511022795.7A CN105510967B (zh) | 2015-12-30 | 2015-12-30 | 一种实现三维地震数据偏移归位的处理方法、装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201511022795.7A CN105510967B (zh) | 2015-12-30 | 2015-12-30 | 一种实现三维地震数据偏移归位的处理方法、装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105510967A true CN105510967A (zh) | 2016-04-20 |
CN105510967B CN105510967B (zh) | 2018-01-05 |
Family
ID=55719075
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201511022795.7A Active CN105510967B (zh) | 2015-12-30 | 2015-12-30 | 一种实现三维地震数据偏移归位的处理方法、装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105510967B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108107473A (zh) * | 2017-12-12 | 2018-06-01 | 中国石油天然气集团有限公司 | 地震偏移时间剖面构造畸变现象的校正方法 |
CN112305585A (zh) * | 2019-07-26 | 2021-02-02 | 中国石油天然气集团有限公司 | 三维数据体获取方法及装置 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201408A (zh) * | 2006-12-12 | 2008-06-18 | 中国石油天然气集团公司 | 一种复杂构造二维地震资料平点控制成图方法 |
CN101329406A (zh) * | 2007-06-20 | 2008-12-24 | 中国石油天然气集团公司 | 一种提高构造图准确度的断层精确归位方法 |
CN101598806A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种提高构造图准确度的偏差消除方法 |
CN104597496A (zh) * | 2015-01-30 | 2015-05-06 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种二维地震资料中速度的三维空间归位方法 |
AU2015201747A1 (en) * | 2014-04-14 | 2015-10-29 | Pgs Geophysical As | Seismic Data Acquisition |
-
2015
- 2015-12-30 CN CN201511022795.7A patent/CN105510967B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101201408A (zh) * | 2006-12-12 | 2008-06-18 | 中国石油天然气集团公司 | 一种复杂构造二维地震资料平点控制成图方法 |
CN101329406A (zh) * | 2007-06-20 | 2008-12-24 | 中国石油天然气集团公司 | 一种提高构造图准确度的断层精确归位方法 |
CN101598806A (zh) * | 2008-06-04 | 2009-12-09 | 中国石油天然气集团公司 | 一种提高构造图准确度的偏差消除方法 |
AU2015201747A1 (en) * | 2014-04-14 | 2015-10-29 | Pgs Geophysical As | Seismic Data Acquisition |
CN104597496A (zh) * | 2015-01-30 | 2015-05-06 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种二维地震资料中速度的三维空间归位方法 |
Non-Patent Citations (6)
Title |
---|
刘俊友: "二维地震资料进行三维归位时求取速度的方法", 《中国煤田地质》 * |
张宪哲: "利用PC―1500计算机实现地震二维勘探三维归位", 《中国煤田地质》 * |
张纪: "2D/3D工区构造成图技术流程探讨", 《海洋石油》 * |
徐佩芬,朱华荣: "二维地震资料三维归位的综合法及其应用效果", 《中国矿业大学学报》 * |
李强等: "利用三维速度场控制实现二维叠前深度偏移", 《石油地球物理勘探》 * |
陆光辉等: "二维偏移剖面闭合方法研究及应用", 《石油物探》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108107473A (zh) * | 2017-12-12 | 2018-06-01 | 中国石油天然气集团有限公司 | 地震偏移时间剖面构造畸变现象的校正方法 |
CN108107473B (zh) * | 2017-12-12 | 2020-01-10 | 中国石油天然气集团有限公司 | 地震偏移时间剖面构造畸变现象的校正方法 |
CN112305585A (zh) * | 2019-07-26 | 2021-02-02 | 中国石油天然气集团有限公司 | 三维数据体获取方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
CN105510967B (zh) | 2018-01-05 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102395902B (zh) | 使用快速面向目标照明计算的地震成像系统及方法 | |
US9759826B2 (en) | System and method for generating an implicit model of geological horizons | |
US9632192B2 (en) | Method of processing seismic data by providing surface offset common image gathers | |
AU710968B2 (en) | Method and apparatus for seismic signal processing and exploration | |
US6574564B2 (en) | 3D prestack seismic data migration method | |
CN102282481B (zh) | 基于地震能见度分析的数据采集和叠前偏移 | |
CA2818790C (en) | Seismic trace attribute | |
CA2828425A1 (en) | Sensitivity kernel-based migration velocity analysis in 3d anisotropic media | |
CN104375188A (zh) | 一种地震波透射衰减补偿方法及装置 | |
CN109001813A (zh) | 一种压制多次波的方法、装置及系统 | |
CN102636809B (zh) | 一种传播角度域共成像点道集的生成方法 | |
CN113015926A (zh) | 无源地震成像 | |
US11397273B2 (en) | Full waveform inversion in the midpoint-offset domain | |
US20150234065A1 (en) | Correction of sea surface state | |
Li et al. | Finite difference seismic forward modeling method for fluid–solid coupled media with irregular seabed interface | |
Poliannikov et al. | Interferometric hydrofracture microseism localization using neighboring fracture | |
CN104570116A (zh) | 基于地质标志层的时差分析校正方法 | |
CN105093296A (zh) | 一种优化观测系统的方法及装置 | |
CN105510967A (zh) | 一种实现三维地震数据偏移归位的处理方法、装置 | |
CN108957545A (zh) | 气枪阵列子波方向性反褶积方法及系统 | |
AU2019243987B2 (en) | Wavefield propagator for tilted orthorhombic media | |
CN107422375B (zh) | 海底节点的共反射点道集的确定方法和装置 | |
CN107229074A (zh) | 低勘探程度海相探区的变速成图方法和装置 | |
US20160047925A1 (en) | Method of Determining Seismic Acquisition Aperture | |
CN110858000B (zh) | 地震数据重构方法、装置、计算机设备及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |