CN112068203A - 提高地震数据纵向分辨率的方法及装置 - Google Patents

提高地震数据纵向分辨率的方法及装置 Download PDF

Info

Publication number
CN112068203A
CN112068203A CN202011050625.0A CN202011050625A CN112068203A CN 112068203 A CN112068203 A CN 112068203A CN 202011050625 A CN202011050625 A CN 202011050625A CN 112068203 A CN112068203 A CN 112068203A
Authority
CN
China
Prior art keywords
reflected wave
seismic data
data
frequency
seismic
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.)
Pending
Application number
CN202011050625.0A
Other languages
English (en)
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.)
Petrochina Co Ltd
Original Assignee
Petrochina 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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN202011050625.0A priority Critical patent/CN112068203A/zh
Publication of CN112068203A publication Critical patent/CN112068203A/zh
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V1/00Seismology; Seismic or acoustic prospecting or detecting
    • G01V1/28Processing seismic data, e.g. for interpretation or for event detection
    • G01V1/36Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
    • G01V1/362Effecting static or dynamic corrections; Stacking
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/50Corrections or adjustments related to wave propagation

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

本发明公开了一种提高地震数据纵向分辨率的方法及装置,其中该方法包括:对获取的地震数据进行预处理;根据预处理后的地震数据确定每个炮点和检波点的平均反射波质心频率;对平均反射波质心频率进行偏移距校正;利用校正后的平均反射波质心频率计算每个炮点和检波点的Q值;利用Q值及稳定的吸收衰减补偿算法对地震数据中的反射波数据进行吸收衰减补偿处理,输出补偿后的反射波数据。本发明可以提高地震数据纵向分辨率。

Description

提高地震数据纵向分辨率的方法及装置
技术领域
本发明涉及石油天然气勘探开发技术领域,尤其涉及一种提高地震数据纵向分辨率的方法及装置。
背景技术
本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。
地层的吸收衰减效应是影响地震资料纵向分辨率的主要因素之一,随着油气勘探对地震分辨率的要求越来越高,如何有效的消除这种影响尤为重要。吸收衰减补偿作为一种理论先进且物理意义明确的提高分辨率的方法,受到越来越多的关注。
吸收衰减补偿的关键在于如何准确建立地层品质因子Q模型,目前的研究表明,Q模型的建立方法大致可以分为四类,一是利用表层吸收调查数据直接获得实际Q值,但实际生产中,这种数据很少或者没有;二是利用垂直地震剖面(Vertical Seismic Profilling,VSP)等测井数据求取近地表Q值,但VSP井的数量有限,不足以反映全工区的近地表变化情况;三是速度经验公式法,该方法适应性差,难以推广应用;四是利用地震记录初至波统计信号振幅或频率信息,利用谱比法或频移法估算地层Q值,该方法是目前普遍应用的主流方法。但是由于受近地表条件影响,往往初至波很难准确识别。尤其对于可控震源采集资料,环境噪音干扰严重,利用初至波统计的振幅和频率结果往往不准确甚至出现错误,导致估算的Q值误差大,补偿效果差,提高纵向分辨率效果不明显。
发明内容
本发明实施例提供一种提高地震数据纵向分辨率的方法,用以提高地震数据纵向分辨率,该方法包括:
对获取的地震数据进行预处理;
根据预处理后的地震数据确定每个炮点和检波点的平均反射波质心频率;
对平均反射波质心频率进行偏移距校正;
利用校正后的平均反射波质心频率计算每个炮点和检波点的Q值;
利用Q值及稳定的吸收衰减补偿算法对地震数据中的反射波数据进行吸收衰减补偿处理,输出补偿后的反射波数据。
本发明实施例还提供一种提高地震数据纵向分辨率的装置,用以提高地震数据纵向分辨率,该装置包括:
预处理模块,用于对获取的地震数据进行预处理;
确定模块,用于根据预处理后的地震数据确定每个炮点和检波点的平均反射波质心频率;
校正模块,用于对平均反射波质心频率进行偏移距校正;
计算模块,用于利用校正后的平均反射波质心频率计算每个炮点和检波点的Q值;
补偿模块,用于利用Q值及稳定的吸收衰减补偿算法对地震数据中的反射波数据进行吸收衰减补偿处理,输出补偿后的反射波数据。
本发明实施例还提供一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,所述处理器执行所述计算机程序时实现上述提高地震数据纵向分辨率方法。
本发明实施例还提供一种计算机可读存储介质,所述计算机可读存储介质存储有执行上述提高地震数据纵向分辨率方法的计算机程序。
本发明实施例中,基于统计的地震数据中反射波的质心频率计算Q值,克服了其他Q值求取方法对特定资料的需求,实用性更广泛,易于推广;尤其对于可控震源采集地震数据而言,利用反射波的质心频率确定的Q值更准确合理;并且,利用稳定的吸收衰减补偿算法,进行地层反射波数据补偿,可以同时补偿振幅和相位,补偿效果好,不改变地震波相对关系,保真度高,同时也实现了提高地震数据纵向分辨率的目的。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。在附图中:
图1为本发明实施例中一种提高地震数据纵向分辨率的方法的流程图;
图2为本发明实施例中经过预处理后的地震反射波单炮数据的示意图;
图3为本发明实施例中对图2所示的地震反射波单炮数据进行频率信息统计的频率属性平面图;
图4为本发明所提供的反射波平均质心频率随偏移距变化的示意图;
图5为本发明所提供的提高地震数据纵向分辨率方法的品质因子Q值属性平面图;
图6为本发明实施例中提高纵向分辨率后地震反射波单炮数据;
图7为本发明实施例中补偿后的Q数值的属性平面图;
图8为本发明实施例中补偿后的反射波质心频率的属性平面图;
图9为本发明实施例中补偿前的叠加剖面图;
图10为本发明实施例中补偿后的叠加剖面图;
图11为本发明实施例中补偿前后叠加剖面对应的频谱图;
图12为本发明实施例中一种提高地震数据纵向分辨率的装置的结构示意图;
图13为本发明实施例中一种计算机设备的结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚明白,下面结合附图对本发明实施例做进一步详细说明。在此,本发明的示意性实施例及其说明用于解释本发明,但并不作为对本发明的限定。
本发明实施例提供了一种提高地震数据纵向分辨率的方法,如图1所示,该方法包括步骤101至步骤105:
步骤101、对获取的地震数据进行预处理。
其中,预处理方法选择的原则是不改变频率和振幅特征,主要采用的方法包括球面扩散补偿处理和去噪处理。示例性的,经过预处理后的地震反射波单炮数据如图2所示。
预处理之后得到的反射波数据的频率特征比较稳定,统计的频率信息可靠性高,提高了品质因子Q值计算结果的准确性,进而为合理提高地震数据纵向分辨率奠定基础。
步骤102、根据预处理后的地震数据确定每个炮点和检波点的平均反射波质心频率。
具体的,步骤102根据预处理后的地震数据确定每个炮点和检波点的平均反射波质心频率可以执行为如下步骤1021至步骤1023:
步骤1021、应用快速傅里叶变换将预处理后的地震数据从时间域转换至频率域。
应用快速傅里叶变换(FFT)实现时间域至频率域的转换是已经成熟的现有技术,对于具体实现过程,在此不做赘述。
步骤1022、根据转换至频率域的地震数据确定每个反射波地震道的质心频率。
其中,对转换至频率域的地震数据进行逐道频率信息统计,得到每个反射波地震道的质心频率。示例性的,如图3所示,为对如图2所示的地震反射波单炮数据进行频率信息统计的频率属性平面图。
步骤1023、根据每个反射波地震道的质心频率计算每个炮点和检波点的平均质心频率。
具体的,根据如下公式计算每个炮点和检波点的平均质心频率:
Figure BDA0002709426890000041
上式中,A(fr)表示反射波振幅谱,fr表示反射波频率,frc表示每个反射波地震道的质心频率,i表示炮点或检波点对应的道数,N表示炮点或检波点对应的总道数。
步骤103、对平均反射波质心频率进行偏移距校正。
进行偏移距校正时,选取质心频率随偏移距变化比较稳定的范围,进一步提高Q值计算结果的可靠性。
示例性的,如图4所示,为反射波平均质心频率随偏移距变化的示意图。
步骤104、利用校正后的平均反射波质心频率计算每个炮点和检波点的Q值。
需要说明的是,计算得到的Q值,也即炮点和检波点的Q模型。
示例性的,图5提供了一种利用本发明实施例中的方法计算得到的Q值的属性平面图。
具体的,根据如下公式计算每个炮点和检波点的Q值:
Figure BDA0002709426890000051
上式中,frm表示期望的反射波频率值,该值由用户设定,但需要注意的是frm必须大于frc;t表示近地表旅行时,可以从地震数据中获取到。
步骤105、利用Q值及稳定的吸收衰减补偿算法对地震数据中的反射波数据进行吸收衰减补偿处理,输出补偿后的反射波数据。
在本步骤中,本发明实施例采用带增益控制的稳健Q补偿方法对地震数据中的反射波数据进行补偿,从而实现提高地震数据纵向分辨率的效果。如图6所示,为本发明实施例中提高纵向分辨率后的地震反射波单炮数据。
补偿的结果如图7至11所示,图7与图5比较可见,补偿后的Q数值较补偿前显著提高;图8和图3比较可见,补偿后的反射波质心频率数值明显高于补偿前的质心频率统计结果;从图9、图10和图11的叠加剖面及频谱对比可见,补偿后的剖面纵向分辨率明显高于补偿前,即达到了有效提高地震数据纵向分辨率的目的。
本发明实施例中,为了验证步骤101至步骤105方法的有效性,还可以对补偿后的反射波数据进行质量控制,确保达到有效提高地震数据纵向分辨率的效果。其中,进行质量控制的方法包括地震道集数据比较法、平面属性分析法、叠加剖面对比法和频谱对比分析法中的任意一种。
本发明实施例中,基于统计的地震数据中反射波的质心频率计算Q值,克服了其他Q值求取方法对特定资料的需求,实用性更广泛,易于推广;尤其对于可控震源采集地震数据而言,利用反射波的质心频率确定的Q值更准确合理;并且,利用稳定的吸收衰减补偿算法,进行地层反射波数据补偿,可以同时补偿振幅和相位,补偿效果好,不改变地震波相对关系,保真度高,同时也实现了提高地震数据纵向分辨率的目的。
本发明实施例中还提供了一种提高地震数据纵向分辨率的装置,如下面的实施例所述。由于该装置解决问题的原理与提高地震数据纵向分辨率的方法相似,因此该装置的实施可以参见提高地震数据纵向分辨率的方法的实施,重复之处不再赘述。
如图12所示,装置1200包括预处理模块1201、确定模块1202、校正模块1203、计算模块1204和补偿模块1205。
其中,预处理模块1201,用于对获取的地震数据进行预处理;
确定模块1202,用于根据预处理后的地震数据确定每个炮点和检波点的平均反射波质心频率;
校正模块1203,用于对平均反射波质心频率进行偏移距校正;
计算模块1204,用于利用校正后的平均反射波质心频率计算每个炮点和检波点的Q值;
补偿模块1205,用于利用Q值及稳定的吸收衰减补偿算法对地震数据中的反射波数据进行吸收衰减补偿处理,输出补偿后的反射波数据。
在本发明实施例的一种实现方式中,预处理模块1201,用于:
对获取的地震数据进行球面扩散补偿处理和去噪处理。
在本发明实施例的一种实现方式中,确定模块1202,用于:
应用快速傅里叶变换将预处理后的地震数据从时间域转换至频率域;
根据转换至频率域的地震数据确定每个反射波地震道的质心频率;
根据每个反射波地震道的质心频率计算每个炮点和检波点的平均质心频率。
在本发明实施例的一种实现方式中,装置1200还包括:
质量控制模块1206,用于对补偿后的反射波数据进行质量控制,其中,进行质量控制的方法包括地震道集数据比较法、平面属性分析法、叠加剖面对比法和频谱对比分析法中的任意一种。
本发明实施例中,基于统计的地震数据中反射波的质心频率计算Q值,克服了其他Q值求取方法对特定资料的需求,实用性更广泛,易于推广;尤其对于可控震源采集地震数据而言,利用反射波的质心频率确定的Q值更准确合理;并且,利用稳定的吸收衰减补偿算法,进行地层反射波数据补偿,可以同时补偿振幅和相位,补偿效果好,不改变地震波相对关系,保真度高,同时也实现了提高地震数据纵向分辨率的目的。
本发明实施例还提供一种计算机设备,图13为本发明实施例中计算机设备的示意图,该计算机设备能够实现上述实施例中的提高地震数据纵向分辨率的方法中全部步骤,该计算机设备具体包括如下内容:
处理器(processor)1301、存储器(memory)1302、通信接口(CommunicationsInterface)1303和通信总线1304;
其中,所述处理器1301、存储器1302、通信接口1303通过所述通信总线1304完成相互间的通信;所述通信接口1303用于实现相关设备之间的信息传输;
所述处理器1301用于调用所述存储器1302中的计算机程序,所述处理器执行所述计算机程序时实现上述实施例中的提高地震数据纵向分辨率的方法。
本领域内的技术人员应明白,本发明的实施例可提供为方法、系统、或计算机程序产品。因此,本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,所应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限定本发明的保护范围,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (10)

1.一种提高地震数据纵向分辨率的方法,其特征在于,所述方法包括:
对获取的地震数据进行预处理;
根据预处理后的地震数据确定每个炮点和检波点的平均反射波质心频率;
对平均反射波质心频率进行偏移距校正;
利用校正后的平均反射波质心频率计算每个炮点和检波点的Q值;
利用Q值及稳定的吸收衰减补偿算法对地震数据中的反射波数据进行吸收衰减补偿处理,输出补偿后的反射波数据。
2.根据权利要求1所述的方法,其特征在于,对获取的地震数据进行预处理,包括:
对获取的地震数据进行球面扩散补偿处理和去噪处理。
3.根据权利要求1所述的方法,其特征在于,根据预处理后的地震数据确定每个炮点和检波点的平均反射波质心频率,包括:
应用快速傅里叶变换将预处理后的地震数据从时间域转换至频率域;
根据转换至频率域的地震数据确定每个反射波地震道的质心频率;
根据每个反射波地震道的质心频率计算每个炮点和检波点的平均质心频率。
4.根据权利要求1至3任一项所述的方法,其特征在于,在输出补偿后的反射波数据之后,所述方法还包括:
对补偿后的反射波数据进行质量控制,其中,进行质量控制的方法包括地震道集数据比较法、平面属性分析法、叠加剖面对比法和频谱对比分析法中的任意一种。
5.一种提高地震数据纵向分辨率的装置,其特征在于,所述装置包括:
预处理模块,用于对获取的地震数据进行预处理;
确定模块,用于根据预处理后的地震数据确定每个炮点和检波点的平均反射波质心频率;
校正模块,用于对平均反射波质心频率进行偏移距校正;
计算模块,用于利用校正后的平均反射波质心频率计算每个炮点和检波点的Q值;
补偿模块,用于利用Q值及稳定的吸收衰减补偿算法对地震数据中的反射波数据进行吸收衰减补偿处理,输出补偿后的反射波数据。
6.根据权利要求5所述的装置,其特征在于,预处理模块,用于:
对获取的地震数据进行球面扩散补偿处理和去噪处理。
7.根据权利要求5所述的装置,其特征在于,确定模块,用于:
应用快速傅里叶变换将预处理后的地震数据从时间域转换至频率域;
根据转换至频率域的地震数据确定每个反射波地震道的质心频率;
根据每个反射波地震道的质心频率计算每个炮点和检波点的平均质心频率。
8.根据权利要求5至7任一项所述的装置,其特征在于,所述装置还包括:
质量控制模块,用于对补偿后的反射波数据进行质量控制,其中,进行质量控制的方法包括地震道集数据比较法、平面属性分析法、叠加剖面对比法和频谱对比分析法中的任意一种。
9.一种计算机设备,包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序,其特征在于,所述处理器执行所述计算机程序时实现权利要求1至4任一所述方法。
10.一种计算机可读存储介质,其特征在于,所述计算机可读存储介质存储有执行权利要求1至4任一所述方法的计算机程序。
CN202011050625.0A 2020-09-29 2020-09-29 提高地震数据纵向分辨率的方法及装置 Pending CN112068203A (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011050625.0A CN112068203A (zh) 2020-09-29 2020-09-29 提高地震数据纵向分辨率的方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011050625.0A CN112068203A (zh) 2020-09-29 2020-09-29 提高地震数据纵向分辨率的方法及装置

Publications (1)

Publication Number Publication Date
CN112068203A true CN112068203A (zh) 2020-12-11

Family

ID=73684174

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011050625.0A Pending CN112068203A (zh) 2020-09-29 2020-09-29 提高地震数据纵向分辨率的方法及装置

Country Status (1)

Country Link
CN (1) CN112068203A (zh)

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103412324A (zh) * 2013-07-17 2013-11-27 西安交通大学 一种估计介质品质因子的epifvo 方法
CN105388518A (zh) * 2014-09-04 2016-03-09 中国石油化工股份有限公司 一种质心频率与频谱比联合的井中地震品质因子反演方法
US20160178772A1 (en) * 2013-05-27 2016-06-23 Statoil Petroleum As High Resolution Estimation of Attenuation from Vertical Seismic Profiles
CN109765615A (zh) * 2019-01-10 2019-05-17 中国石油天然气股份有限公司 一种地层品质因子反演方法及装置

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20160178772A1 (en) * 2013-05-27 2016-06-23 Statoil Petroleum As High Resolution Estimation of Attenuation from Vertical Seismic Profiles
CN103412324A (zh) * 2013-07-17 2013-11-27 西安交通大学 一种估计介质品质因子的epifvo 方法
CN105388518A (zh) * 2014-09-04 2016-03-09 中国石油化工股份有限公司 一种质心频率与频谱比联合的井中地震品质因子反演方法
CN109765615A (zh) * 2019-01-10 2019-05-17 中国石油天然气股份有限公司 一种地层品质因子反演方法及装置

Non-Patent Citations (9)

* Cited by examiner, † Cited by third party
Title
K. XIN ET AL.: "Adaptive Centroid Frequency Shift Q Tomography", 《76TH EAGE CONFERENCE & EXHIBITION 2014》 *
WASIU RAJI ET AL.: "Attenuation (1/Q) estimation in reflection seismic records", 《J. GEOPHYS. ENG》 *
何兵红等: "基于叠前地震道集的FVO分析方法", 《石油地球物理勘探》 *
余连勇等: "时间域质心频移法估算品质因子 Q", 《西南石油大学学报(自然科学版)》 *
余青露等: "基于峰值频率法的组合法估计叠前Q值", 《地球物理学进展》 *
崔琴等: "基于加权指数型子波的质心频移法估算品质因子", 《物探与化探》 *
蒋立等: "地表过渡带近地表Q 补偿与地表一致性反褶积处理效果对比研究", 《石油物探》 *
阎世信等: "《山地地球物理勘探技术》", 30 June 2000, 北京:石油工业出版社 *
韩文功等: "《高密度三维地震技术》", 30 November 2017, 北京:地质出版社 *

Similar Documents

Publication Publication Date Title
CN106405651B (zh) 一种基于测井匹配的全波形反演初始速度模型构建方法
CN109669212B (zh) 地震数据处理方法、地层品质因子估算方法与装置
CN107065013B (zh) 一种地震尺度下的层速度确定方法及装置
CN103424777B (zh) 一种提高地震成像分辨率的方法
CN109188520B (zh) 薄储层厚度预测方法及装置
CN112068203A (zh) 提高地震数据纵向分辨率的方法及装置
CN112415601B (zh) 表层品质因子q值的确定方法及装置
CN111435174A (zh) 强反射地区地震资料振幅补偿方法及装置
CN113031070B (zh) 一种深度域合成地震记录的制作方法
CN116413787A (zh) 利用das井中地震数据反演地层波阻抗的方法及装置
CN104991273B (zh) 提取叠前校正道集地震子波的方法
CN115079257A (zh) 基于融合网络的q值估计及地震衰减补偿方法
CN112379429A (zh) 地震数据的振幅补偿方法及装置
CN114114421A (zh) 基于深度学习的导向自学习地震数据去噪方法及装置
CN114563816A (zh) 油气藏评价阶段建立地震解释速度模型的方法及装置
CN112649876A (zh) 建立地震偏移速度模型的方法和装置
CN111880218A (zh) 基于品质因子的反演子波字典构建方法
CN111856557A (zh) 深度域合成地震记录制作方法及装置
CN112526611A (zh) 表层地震波品质因子的提取方法及装置
CN114002741B (zh) 叠前深度偏移方法及装置、计算机可读存储介质
Yao et al. Microseismic signal denoising using simple bandpass filtering based on normal time–frequency transform
CN114755720B (zh) 一种页岩储层地震预测方法、装置、设备及介质
CN112394393B (zh) 一种crp道集数据体重构的方法
CN108345035B (zh) 地球物理勘探密度信息提取的方法和装置
CN112835096B (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
RJ01 Rejection of invention patent application after publication

Application publication date: 20201211

RJ01 Rejection of invention patent application after publication