CN113866827B - 一种解释性速度建模地震成像方法、系统、介质和设备 - Google Patents
一种解释性速度建模地震成像方法、系统、介质和设备 Download PDFInfo
- Publication number
- CN113866827B CN113866827B CN202111149794.4A CN202111149794A CN113866827B CN 113866827 B CN113866827 B CN 113866827B CN 202111149794 A CN202111149794 A CN 202111149794A CN 113866827 B CN113866827 B CN 113866827B
- Authority
- CN
- China
- Prior art keywords
- imaging
- velocity
- seismic
- wave impedance
- model
- 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
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 133
- 238000000034 method Methods 0.000 claims abstract description 25
- 238000001914 filtration Methods 0.000 claims abstract description 8
- 238000005070 sampling Methods 0.000 claims description 3
- 238000004458 analytical method Methods 0.000 claims description 2
- 238000004422 calculation algorithm Methods 0.000 claims description 2
- 238000010276 construction Methods 0.000 abstract description 2
- 238000013508 migration Methods 0.000 description 16
- 230000005012 migration Effects 0.000 description 16
- 238000004590 computer program Methods 0.000 description 7
- 238000010586 diagram Methods 0.000 description 7
- 238000012545 processing Methods 0.000 description 5
- 230000006870 function Effects 0.000 description 4
- 238000004364 calculation method Methods 0.000 description 3
- 230000000694 effects Effects 0.000 description 3
- 238000006467 substitution reaction Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000008569 process 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/30—Analysis
- G01V1/301—Analysis for determining seismic cross-sections or geostructures
-
- 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/282—Application of seismic models, synthetic seismograms
-
- 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/306—Analysis for determining physical properties of the subsurface, e.g. impedance, porosity or attenuation profiles
-
- 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/50—Corrections or adjustments related to wave propagation
- G01V2210/51—Migration
-
- 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/61—Analysis by combining or comparing a seismic data set with other data
- G01V2210/614—Synthetically generated data
-
- 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
-
- 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
-
- 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/67—Wave propagation modeling
- G01V2210/679—Reverse-time modeling or coalescence modelling, i.e. starting from receivers
-
- 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
本发明属于地震勘探成像技术领域,涉及一种解释性速度建模地震成像方法、系统、介质和设备,包括:S1对给定初始速度模型进行第一次成像,获得第一次成像结果;S2对第一次成像结果进行相对波阻抗反演,得到相对波阻抗剖面;S3对相对阻抗剖面进行Curvelet滤波,得到第一次解释方案;S4将第一次解释方案和初始速度模型进行叠加从而得到新的偏移速度场;S5对新的偏移速度场进行第二次成像,获得第二次成像结果;S6对得到的第二次成像结果重复步骤S2‑S4,直到获得最终的地震成像结果。使成像剖面包含更多构造细节,从而使得地震解释人员在成像剖面上得到更加丰富的地质特征。
Description
技术领域
本发明涉及一种解释性速度建模地震成像方法、系统、介质和设备,属于地震勘探技术领域,特别属于地震勘探成像技术领域。
背景技术
速度建模、地震偏移和地震反演是勘探地震学领域的三种传统技术领域。常规速度建模包括动校正叠加、叠前时间、时深转换、叠前深度等偏移速度分析,由于常规速度建模只用到走时信息,平均效应很强,建立的速度场非常平滑,不适用于高频地震数据的偏移成像。采用目前最先进的最小二乘逆时深度偏移,基于如此平滑速度场地震偏移得到的成像剖面包含较少的结构细节。
发明内容
针对上述问题,本发明的目的是提供一种解释性速度建模地震成像方法、系统、介质和设备,其从初始速度模型出发,通过高精度偏移成像及其后续突出结构的数据处理,开展地质结构解释速度建模,再偏移成像,使得成像剖面包含更多构造细节,从而使得地震解释人员在成像剖面上得到更加丰富的地质特征。
为实现上述目的,本发明采取以下技术方案:一种解释性速度建模地震成像方法,包括:S1对给定初始速度模型进行第一次成像,获得第一次成像结果;S2对第一次成像结果进行相对波阻抗反演,得到相对波阻抗剖面;S3对相对阻抗剖面进行Curvelet滤波,得到第一次解释方案;S4将第一次解释方案和初始速度模型进行叠加从而得到新的偏移速度场;S5对新的偏移速度场进行第二次成像,获得第二次成像结果;S6对得到的第二次成像结果重复步骤S2-S4,直到获得最终的地震成像结果。
进一步,第一次成像和第二次成像通过将给定初始速度模型或新的偏移速度场输入最小二乘逆时偏移算法获得。
进一步,相对波阻抗反演以反褶积法为基础进行直接反演。
进一步,相对波阻抗反演的方法是在相对波阻抗剖面上计算相对波阻抗,以偏移速度分析得到标准阻抗,对标准阻抗进行归一化校准,得到相对速度剖面。
进一步,将每一次获得的成像结果分别与真实断溶体模型的成像结果进行相似度的计算,以对各次解释性速度建模的成像结果进行验证。
进一步,给定初始速度模型通过常规速度建模获得,并给出初始速度模型的横向和纵向网格大小及网格数,并通过有限差分方法,得到基于初始速度模型的炮集记录。
进一步,初始速度模型包括以下模型参数:横向和纵向网格大小、横向和纵向网格间距、子波时长和主频、时间采样间隔、时间总长度、震源的个数、震源之间的间距和震源的横纵坐标起始位置。
本发明还包括一种解释性速度建模地震成像系统,包括:一次成像模块,用于对给定初始速度模型进行第一次成像,获得第一次成像结果;相对波阻抗反演模块,用于对第一次成像结果进行相对波阻抗反演,得到相对波阻抗剖面;解释模块,用于对相对阻抗剖面进行Curvelet滤波,得到第一次解释方案;叠加模块,用于将第一次解释方案和初始速度模型进行叠加从而得到新的偏移速度场;第二成像模块,用于对新的偏移速度场进行第二次成像,获得第二次成像结果;循环模块,用于对得到的第二次成像结果输入相对波阻抗反演模块和解释模块进行循环,直到获得最终的地震成像结果。
本发明还包括一种存储一个或多个程序的计算机可读存储介质,一个或多个程序包括指令,指令当由计算设备执行时,使得计算设备执行根据上述任一项的解释性速度建模地震成像方法。
本发明还包括一种计算设备,包括:一个或多个处理器、存储器以及一个或多个程序,其中一个或多个程序存储在存储器中并被配置为由一个或多个处理器执行,一个或多个程序包括用于执行根据上述任一项的解释性速度建模地震成像方法。
本发明由于采取以上技术方案,其具有以下优点:
1、本发明通过更新速度模型,进行成像,得到的结果更加接近于真实速度模型的地震成像结果,增加了第一轮成像剖面未表现出的细节构造,其构造也更加收敛,有效的指导后续的地震解释工作,对于油气勘探开发具有很重要的意义。
2、速度建模、地震偏移和地震反演是勘探地震学的三种传统的技术领域。一般而言,速度建模和地震偏移具有直接关系,速度建模为地震偏移提供初始的速度场;速度建模和地震反演具有直接关系,速度建模为地震反演提供低频信息,实际应用中,地震偏移和地震反演并没有直接的联系,本发明建立了地震偏移和地震反演的联系。本发明可以广泛应用于地震偏移成像领域。
附图说明
图1是本发明一实施例中为解释性速度建模地震成像方法的流程图;
图2是本发明一实施例中真实断溶体速度模型;
图3是本发明一实施例中真实断溶体速度模型的最小二乘逆时偏移成像剖面;
图4是本发明一实施例中真实断溶体速度模型的偏移速度场;
图5是本发明一实施例中真实断溶体速度模型的偏移速度场的第一次成像结果图;
图6是本发明一实施例中真实断溶体速度模型的第一次成像结果图的相对波阻抗剖面;
图7是本发明一实施例中真实断溶体速度模型的基于相对波阻抗剖面的解释方案图;
图8是本发明一实施例中解释方案和偏移速度场的叠加结果,更新后的偏移速度场的图像;
图9是本发明一实施例中第二次成像结果图。
具体实施方式
为了使本领域技术人员更好的理解本发明的技术方向,通过具体实施例对本发明进行详细的描绘。然而应当理解,具体实施方式的提供仅为了更好地理解本发明,它们不应该理解成对本发明的限制。在本发明的描述中,需要理解的是,所用到的术语仅仅是用于描述的目的,而不能理解为指示或暗示相对重要性。
实施例一
本实施例公开了一种解释性速度建模地震成像方法,如图1所示,包括:
本实施例以断溶体速度模型为例进行说明,断溶体速度模型的真实模型图像如图2所示。基于该断溶体速度模型生成合成地震炮集数据,采用最小二乘逆时偏移成像(RTM--Reverse Time Migration),得到图3所示的真实速度模型偏移成像剖面。
S1对给定初始速度模型进行第一次成像,获得第一次成像结果。
给定初始速度模型,如图4所示,其通过常规速度建模获得,并给出初始速度模型的横向和纵向网格大小及网格数,并通过有限差分方法,得到基于初始速度模型的炮集记录。
初始速度模型包括以下模型参数:横向和纵向网格大小、横向和纵向网格间距、子波时长和主频、时间采样间隔、时间总长度、震源的个数、震源之间的间距和震源的横纵坐标起始位置。
S2对第一次成像结果进行相对波阻抗反演,得到相对波阻抗剖面。
第一次成像结果基于偏移速度模型的参数和炮集数据,采用最小二乘逆时偏移方法,对初始速度模型进行成像获得,本实施例中第一次成像结果,如图5所示,其中模型参数与初始速度模型相同
相对波阻抗反演以反褶积法为基础进行直接反演,该方法具有计算简单、不受地质资料和井资料约束的特点。其具体是在相对波阻抗剖面上计算相对波阻抗,以偏移速度分析得到标准阻抗,对标准阻抗进行归一化校准,得到相对速度剖面,该图具体如图6所示,得到相对阻抗。
S3对相对阻抗剖面进行Curvelet滤波,突出地质结构,得到第一次解释方案,其结果如图7所示。
S4将第一次解释方案和初始速度模型进行叠加从而得到新的偏移速度场,如图8所示。
S5对新的偏移速度场进行第二次成像,获得第二次成像结果,如图9所示。
第二次成像结果基于偏移速度模型的参数和炮集数据,采用最小二乘逆时偏移方法,对新的偏移速度场进行成像获得。
S6对得到的第二次成像结果重复步骤S2-S4,直到获得最终的地震成像结果。通过如此迭代,使得地震成像剖面逐渐逼近真实地质结构。
将每一次获得的成像结果分别与真实断溶体模型的成像结果进行相似度的计算,以对各次解释性速度建模的成像结果进行验证,以判断本方法的有效性。将第一次成像和第二次成像结果分别与真实断溶体模型(即基于图2中图像得到的真实速度成像剖面,即图3)所成像结果进行相似度的计算,相似度越大,说明解释性速度建模地震成像方法的效果越好。第一轮成像剖面与真实速度成像剖面相似度为71.46%,第二轮成像剖面与真实速度成像剖面的相似度为81.28%。经过二轮迭代,成像效果改善近10%。
实施例二
基于相同的发明构思,本实施例包括一种解释性速度建模地震成像系统,包括:
一次成像模块,用于对给定初始速度模型进行第一次成像,获得第一次成像结果;
相对波阻抗反演模块,用于对第一次成像结果进行相对波阻抗反演,得到相对波阻抗剖面;
解释模块,用于对相对阻抗剖面进行Curvelet滤波,得到第一次解释方案;
叠加模块,用于将第一次解释方案和初始速度模型进行叠加从而得到新的偏移速度场;
第二成像模块,用于对新的偏移速度场进行第二次成像,获得第二次成像结果;
循环模块,用于对得到的第二次成像结果输入相对波阻抗反演模块和解释模块进行循环,直到获得最终的地震成像结果。
实施例三
基于相同的发明构思,本实施例包括一种存储一个或多个程序的计算机可读存储介质,一个或多个程序包括指令,指令当由计算设备执行时,使得计算设备执行根据上述任一项的解释性速度建模地震成像方法。
实施例四
基于相同的发明构思,本实施例包括一种计算设备,包括:一个或多个处理器、存储器以及一个或多个程序,其中一个或多个程序存储在存储器中并被配置为由一个或多个处理器执行,一个或多个程序包括用于执行根据上述任一项的解释性速度建模地震成像方法。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
最后应当说明的是:以上实施例仅用以说明本发明的技术方案而非对其限制,尽管参照上述实施例对本发明进行了详细的说明,所属领域的普通技术人员应当理解:依然可以对本发明的具体实施方式进行修改或者等同替换,而未脱离本发明精神和范围的任何修改或者等同替换,其均应涵盖在本发明的权利要求保护范围之内。上述内容仅为本申请的具体实施方式,但本申请的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本申请揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本申请的保护范围之内。因此,本申请的保护范围应以权利要求的保护范围。
Claims (10)
1.一种解释性速度建模地震成像方法,其特征在于,包括:
对给定初始速度模型进行第一次成像,获得第一次成像结果;
对所述第一次成像结果进行相对波阻抗反演,得到相对波阻抗剖面;
对相对阻抗剖面进行Curvelet滤波,得到第一次解释方案;
将第一次解释方案和所述初始速度模型进行叠加从而得到新的偏移速度场;
对所述新的偏移速度场进行第二次成像,获得第二次成像结果;
对得到的第二次成像结果重复相对波阻抗反演和Curvelet滤波步骤,直到获得最终的地震成像结果。
2.如权利要求1所述的解释性速度建模地震成像方法,其特征在于,所述第一次成像和第二次成像通过将给定初始速度模型或新的偏移速度场输入最小二乘逆时偏移算法获得。
3.如权利要求1所述的解释性速度建模地震成像方法,其特征在于,所述相对波阻抗反演以反褶积法为基础进行直接反演。
4.如权利要求3所述的解释性速度建模地震成像方法,其特征在于,所述相对波阻抗反演的方法是在相对波阻抗剖面上计算相对波阻抗,以偏移速度分析得到标准阻抗,对所述标准阻抗进行归一化校准,得到相对速度剖面。
5.如权利要求1所述的解释性速度建模地震成像方法,其特征在于,将每一次获得的成像结果分别与真实断溶体模型的成像结果进行相似度的计算,以对各次解释性速度建模的成像结果进行验证。
6.如权利要求1所述的解释性速度建模地震成像方法,其特征在于,所述给定初始速度模型通过常规速度建模获得,并给出所述初始速度模型的横向和纵向网格大小及网格数,并通过有限差分方法,得到基于所述初始速度模型的炮集记录。
7.如权利要求6所述的解释性速度建模地震成像方法,其特征在于,所述初始速度模型包括以下模型参数:横向和纵向网格大小、横向和纵向网格间距、子波时长和主频、时间采样间隔、时间总长度、震源的个数、震源之间的间距和震源的横纵坐标起始位置。
8.一种解释性速度建模地震成像系统,其特征在于,包括:
一次成像模块,用于对给定初始速度模型进行第一次成像,获得第一次成像结果;
相对波阻抗反演模块,用于对所述第一次成像结果进行相对波阻抗反演,得到相对波阻抗剖面;
解释模块,用于对相对阻抗剖面进行Curvelet滤波,得到第一次解释方案;
叠加模块,用于将第一次解释方案和所述初始速度模型进行叠加从而得到新的偏移速度场;
第二成像模块,用于对所述新的偏移速度场进行第二次成像,获得第二次成像结果;
循环模块,用于对得到的第二次成像结果输入相对波阻抗反演模块和解释模块进行循环,直到获得最终的地震成像结果。
9.一种存储一个或多个程序的计算机可读存储介质,其特征在于,所述一个或多个程序包括指令,所述指令当由计算设备执行时,使得所述计算设备执行根据权利要求1至7任一项所述的解释性速度建模地震成像方法。
10.一种计算设备,其特征在于,包括:一个或多个处理器、存储器以及一个或多个程序,其中一个或多个程序存储在所述存储器中并被配置为由所述一个或多个处理器执行,所述一个或多个程序包括用于执行根据权利要求1至7任一项所述的解释性速度建模地震成像方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111149794.4A CN113866827B (zh) | 2021-09-29 | 2021-09-29 | 一种解释性速度建模地震成像方法、系统、介质和设备 |
US17/942,110 US11686870B2 (en) | 2021-09-29 | 2022-09-09 | Interpretive-guided velocity modeling seismic imaging method and system, medium and device |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111149794.4A CN113866827B (zh) | 2021-09-29 | 2021-09-29 | 一种解释性速度建模地震成像方法、系统、介质和设备 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113866827A CN113866827A (zh) | 2021-12-31 |
CN113866827B true CN113866827B (zh) | 2023-05-09 |
Family
ID=78992462
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111149794.4A Active CN113866827B (zh) | 2021-09-29 | 2021-09-29 | 一种解释性速度建模地震成像方法、系统、介质和设备 |
Country Status (2)
Country | Link |
---|---|
US (1) | US11686870B2 (zh) |
CN (1) | CN113866827B (zh) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117094233B (zh) * | 2023-10-19 | 2024-01-30 | 中国科学院地质与地球物理研究所 | 深地油气精准导航三维速度场高精度建模方法与系统 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6044039A (en) * | 1996-07-19 | 2000-03-28 | Elf Exploration Production | Geophysical attribute migration method for a medium |
CN105095634A (zh) * | 2014-05-23 | 2015-11-25 | 中国石油化工股份有限公司 | 一种基于地质块体的偏移速度模型建立方法 |
CN105093281A (zh) * | 2014-05-16 | 2015-11-25 | 中国石油化工股份有限公司 | 一种反演框架下的地震多波建模方法 |
CN109738945A (zh) * | 2018-11-08 | 2019-05-10 | 成都捷科思石油天然气技术发展有限公司 | 一种利用叠前深度偏移成果直接生成构造图的方法 |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9910174B2 (en) * | 2014-07-25 | 2018-03-06 | Seoul National University R&Db Foundation | Seismic imaging apparatus and method for performing iterative application of direct waveform inversion |
US9470811B2 (en) * | 2014-11-12 | 2016-10-18 | Chevron U.S.A. Inc. | Creating a high resolution velocity model using seismic tomography and impedance inversion |
KR102020759B1 (ko) * | 2015-10-02 | 2019-09-11 | 엑손모빌 업스트림 리서치 캄파니 | Q-보상된 전 파동장 반전 |
US10698126B2 (en) * | 2016-08-12 | 2020-06-30 | Exxonmobil Upstream Research Company | Tomographically enhanced full wavefield inversion |
BR112021011246A2 (pt) * | 2018-12-11 | 2021-08-24 | Exxonmobil Upstream Research Company | Inversão guiada por interpretação sísmica automatizada |
US10908308B1 (en) * | 2019-07-25 | 2021-02-02 | Chevron U.S.A. Inc. | System and method for building reservoir property models |
CN111505713B (zh) | 2020-01-21 | 2021-05-07 | 长江大学 | 基于多点地质统计的叠前地震反演方法 |
CN111538080B (zh) | 2020-05-28 | 2022-08-26 | 中国石油天然气集团有限公司 | 地震成像的方法 |
CN112230283B (zh) | 2020-10-12 | 2021-08-10 | 北京中恒利华石油技术研究所 | 一种基于测井曲线支持向量机建模的地震孔隙度预测方法 |
CN114428324B (zh) | 2022-04-06 | 2022-06-28 | 中国石油大学(华东) | 叠前高角度快速傅里叶变换地震成像方法、系统、设备 |
-
2021
- 2021-09-29 CN CN202111149794.4A patent/CN113866827B/zh active Active
-
2022
- 2022-09-09 US US17/942,110 patent/US11686870B2/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US6044039A (en) * | 1996-07-19 | 2000-03-28 | Elf Exploration Production | Geophysical attribute migration method for a medium |
CN105093281A (zh) * | 2014-05-16 | 2015-11-25 | 中国石油化工股份有限公司 | 一种反演框架下的地震多波建模方法 |
CN105095634A (zh) * | 2014-05-23 | 2015-11-25 | 中国石油化工股份有限公司 | 一种基于地质块体的偏移速度模型建立方法 |
CN109738945A (zh) * | 2018-11-08 | 2019-05-10 | 成都捷科思石油天然气技术发展有限公司 | 一种利用叠前深度偏移成果直接生成构造图的方法 |
Non-Patent Citations (1)
Title |
---|
起伏地表条件下的混叠数据最小二乘逆时偏移;周东红 等;《地球物理学进展》;第35卷(第4期);第1528-1538页 * |
Also Published As
Publication number | Publication date |
---|---|
US20230095632A1 (en) | 2023-03-30 |
US11686870B2 (en) | 2023-06-27 |
CN113866827A (zh) | 2021-12-31 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105467444B (zh) | 一种弹性波全波形反演方法及装置 | |
CN105277978B (zh) | 一种确定近地表速度模型的方法及装置 | |
Biswas et al. | 2D full-waveform inversion and uncertainty estimation using the reversible jump Hamiltonian Monte Carlo | |
CN103149586B (zh) | 一种倾斜层状粘弹性介质中波场正演模拟方法 | |
CN107894618B (zh) | 一种基于模型平滑算法的全波形反演梯度预处理方法 | |
CN104977607B (zh) | 利用变步长网格声波波场模拟的时间域全波形反演方法 | |
CN107765308B (zh) | 基于褶积思想与精确震源的重构低频数据频域全波形反演方法 | |
CN106932819A (zh) | 基于各向异性马尔科夫随机域的叠前地震参数反演方法 | |
KR20170118185A (ko) | 다중반사파 없는 데이터 세트를 생성하는 다단식 전 파동장 역산 프로세스 | |
CN113866827B (zh) | 一种解释性速度建模地震成像方法、系统、介质和设备 | |
WO2016176235A1 (en) | Time lapse seismic data processing | |
CN111580163B (zh) | 一种基于非单调搜索技术的全波形反演方法及系统 | |
CN108680968B (zh) | 复杂构造区地震勘探数据采集观测系统评价方法及装置 | |
CN109738944B (zh) | 基于广角反射的地震采集参数确定方法及装置 | |
CN111435174B (zh) | 强反射地区地震资料振幅补偿方法及装置 | |
CN112649887A (zh) | 基于钻井资料定量划分层序地层的方法及装置 | |
CN111208568B (zh) | 一种时间域多尺度全波形反演方法及系统 | |
CN113900141B (zh) | 油气分布预测方法及装置 | |
CN115659773A (zh) | 基于深度网络的全波形反演加速方法及相关装置 | |
CN110376642B (zh) | 一种基于锥面波的三维地震速度反演方法 | |
CN112415601A (zh) | 表层品质因子q值的确定方法及装置 | |
CN113031072A (zh) | 虚同相轴层间的多次波压制方法、装置及设备 | |
CN110161565A (zh) | 一种地震数据重建方法 | |
CN112558150A (zh) | 一种高效的频率空间域粘声介质高斯束正演系统 | |
CN115421195A (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 |