CN111665563A - 基于聚焦分析的叠前偏移垂向分辨率评价方法 - Google Patents

基于聚焦分析的叠前偏移垂向分辨率评价方法 Download PDF

Info

Publication number
CN111665563A
CN111665563A CN201910166749.6A CN201910166749A CN111665563A CN 111665563 A CN111665563 A CN 111665563A CN 201910166749 A CN201910166749 A CN 201910166749A CN 111665563 A CN111665563 A CN 111665563A
Authority
CN
China
Prior art keywords
resolution
seismic
imaging
wave
evaluating
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
CN201910166749.6A
Other languages
English (en)
Other versions
CN111665563B (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.)
Sinopec Oilfield Service Corp
Sinopec Petroleum Engineering Geophysics Co Ltd
Sinopec Petroleum Engineering Geophysics Co Ltd Shengli Branch
Original Assignee
Sinopec Oilfield Service Corp
Sinopec Petroleum Engineering Geophysics Co Ltd
Sinopec Petroleum Engineering Geophysics Co Ltd Shengli Branch
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 Sinopec Oilfield Service Corp, Sinopec Petroleum Engineering Geophysics Co Ltd, Sinopec Petroleum Engineering Geophysics Co Ltd Shengli Branch filed Critical Sinopec Oilfield Service Corp
Priority to CN201910166749.6A priority Critical patent/CN111665563B/zh
Publication of CN111665563A publication Critical patent/CN111665563A/zh
Application granted granted Critical
Publication of CN111665563B publication Critical patent/CN111665563B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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. analysis, for interpretation, for correction
    • G01V1/34Displaying seismic recordings or visualisation of seismic data or attributes
    • G01V1/345Visualisation of seismic data or attributes, e.g. in 3D cubes
    • 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. analysis, for interpretation, for correction
    • G01V1/30Analysis

Abstract

本发明提供一种基于聚焦分析的叠前偏移垂向分辨率评价方法,包括:步骤1,建立复杂构造区的速度模型;步骤2,从目标点激发地震波,使得波场上行至地表并接收;步骤3,当上行波场达到地表后,应用观测系统获得相应的地震记录;步骤4,由采集的地震记录,通过单程波外推技术,进行反向外推;步骤5,采用叠前深度偏移的零时刻成像原理对波场进行成像;步骤6,获得整个空间上的聚焦分辨率图像;步骤7,获得震源、检波点对目标点成像质量影响的分析结果;步骤8,进行观测系统排列的调整及优化。该基于聚焦分析的叠前偏移垂向分辨率评价方法能够精确地对地震观测系统的垂向成像效果进行评价分析。

Description

基于聚焦分析的叠前偏移垂向分辨率评价方法
技术领域
本发明涉及地震采集观测系统评价领域,特别是涉及到一种基于聚焦分析的叠前偏移垂向分辨率评价方法。
背景技术
在三维勘探的观测系统设计中,除了希望获得高质量、高信噪比的数据之外,还希望获得更高的成像精度,但基于层状介质假设的设计方法无法很好地反映叠前偏移成像精度。目前观测系统评价方法所得到的覆盖次数和面元分析结果对成像精度的评价比较有限,且无垂向成像效果的评价方法。为此我们发明了一种新的基于聚焦分析的叠前偏移垂向分辨率评价方法,解决了以上技术问题。
发明内容
本发明的目的是提供一种基于聚焦分析的垂向分辨率评价方法,针对给定的观测系统,确定其叠前深度偏移的垂向分辨率,同时为了提高计算效率,提出了并行计算基于聚焦分析的叠前偏移垂向分辨率评价的方法。
本发明的目的可通过如下技术措施来实现:基于聚焦分析的叠前偏移垂向分辨率评价方法,该基于聚焦分析的叠前偏移垂向分辨率评价方法包括:步骤1,建立复杂构造区的速度模型;步骤2,从目标点激发地震波,使得波场上行至地表并接收;步骤3,当上行波场达到地表后,应用观测系统获得相应的地震记录;步骤4,由采集的地震记录,通过单程波外推技术,进行反向外推;步骤5,采用叠前深度偏移的零时刻成像原理对波场进行成像;步骤6,获得整个空间上的聚焦分辨率图像;步骤7,获得震源、检波点对目标点成像质量影响的分析结果;步骤8,进行观测系统排列的调整及优化。
本发明的目的可通过如下技术措施来实现:
在步骤1中,根据需要选择地震照明效果不佳的区域作为评价的目标点,收集工区前期地质、地震、测录井和钻井资料,结合地震解释成果和地质认识,建立工区纵波速度模型,并根据勘探要求确定地质目标。
在步骤2中,从目标点激发震源,采用高阶有限差分离散双程声波波动方程进行波场模拟。
在步骤3中,根据步骤2中模拟的地震波场传播至地表,应用地表布设的观测排列对波场进行采样处理,分别得到震源排列对应的记录和检波点排列对应的记录。
在步骤4中,由地表向目标点深度进行波场延拓,聚焦算子采用傅立叶有限差分算子,将步骤3中地表检波点所得共聚焦点道集,作为波场外推的边界,其计算过程完全模拟地震波上行波传播的逆过程,由此,将震源排列成像分辨率Rs,检波点排列成像分辨率Rd定义如式(1):
Figure BDA0001985523340000021
其中:I表示单位矩阵,S、D分别表示随着深度聚焦后的震源、检波点采样记录,Fs、Fd表示震源、检波点聚焦算子。
在步骤5中,在深度延拓的过程中,采用叠前深度偏移的零时刻成像原理对波场进行成像,根据公式(2)获取当前深度的在各个方向的成像分辨率值。
在步骤6中,根据等式(1)中震源、检波点分辨率的定义,对所有深度进行步骤5中的成像处理,获得整个空间上的聚焦分辨率图像,在聚焦过程中,针对前人无垂向分辨率评价的方法,垂向、横向分辨率图像提取方法按照下式进行:
Rxoy(x,y)=U(x,y,z=sz,t=0) (2)
其中:R表示各个平面内的聚焦分辨率图像,U表示在延拓记录,sx,sy,sz表示目标点坐标,XOZ与YOZ为垂向分辨率函数,XOY为横向分辨率函数。
在步骤7中,通过对震源和检波点排列分辨率图像进行处理,分别获得震源、检波点对目标点成像质量的影响,将震源和检波点排列聚焦图像的乘积作为当前观测系统在该目标区域的总成像分辨率图像。
在步骤7中,在滚动排列情况中,检波点排列的成像分辨率与炮点排列有关,将滚动排列情况下的检波点排列成像分辨率RD_roll定义如式(3)所示:
Figure BDA0001985523340000031
其中,s为当前炮索引,NumShot为观测系统总炮数,RD为当前炮对应的检波点聚焦束结果。
在步骤7中,在实际处理中,当炮点数目过多,会导致计算量过大,采用等式(4),每一次聚焦过程中,隐式差分计算完成相邻深度记录的更新,其中,p表示不同深度的延拓记录,a、a*是与模型、频率等固有参数的变量;因此,多炮情况下的检波点聚焦过程存在独立性,即实现p1,p2,p3和p4的聚焦计算可并行处理;
Figure BDA0001985523340000032
其中,z表示深度,k表示当前深度,1表示单位1。
在步骤8中,根据分辨率图像的结果进行观测系统排列的调整及优化,分辨率图像能量团越集中,成像分辨率越高。
本发明中的基于聚焦分析的叠前偏移垂向分辨率评价方法,是针对复杂构造区面向地质目标的叠前深度偏移垂向分辨率的评价方法,能够精确地对地震观测系统的垂向成像效果进行评价分析。本发明的主要优势在于:(1)在共聚焦点模拟的计算过程中采用声波波动方程方法,波动方程能够精确地描述地震波场在地下介质中的传播特征,在复杂构造地区具有更好的适用性。(2)该方法中垂向分辨率的评价计算基于波动方程叠前深度偏移技术,其评价结果与真实偏移结果更为接近。(3)根据该方法中获得的垂向分辨率图形,能够从空间上,直观地评价当前观测系统对目标点的成像精度和质量,克服了以往只能从横向上评价目标点成像精度方法的不足。(4)本发明针对计算效率问题,提出了并行计算算法。
附图说明
图1为本发明的一个具体实例中用于实验的工区三维模型示意图;
图2为本发明的一个具体实例中三维工区的观测系统图;
图3为本发明的一个具体实例中所计算所得纵向分辨率图;
图4为本发明的基于聚焦分析的叠前偏移垂向分辨率评价方法的一个具体实例的流程图。
具体实施方式
为使本发明的上述和其他目的、特征和优点能更明显易懂,下文特举出较佳实施例,并配合附图所示,作详细说明如下。
如图4所示,图4为本发明的基于聚焦分析的叠前偏移垂向分辨率评价方法的流程图。
步骤101,建立复杂构造区的速度模型,根据需要选择地震照明效果不佳的区域作为评价的目标点。收集工区前期地质、地震、测录井和钻井资料,结合地震解释成果和地质认识,建立如图1所示的工区纵波速度模型,并根据勘探要求确定地质目标,这里定为(X=100,Y=80,Z=40)。
步骤102,从目标点激发地震波,使得波场上行至地表并接收。
由于聚焦计算对象是震源与下行波传播算子、检波点与上行波算子,首先要生成对应的聚焦点道集。其过程为,从目标点激发震源,采用高阶有限差分离散双程声波波动方程进行波场模拟。
步骤103,当上行波场达到地表后,应用当前评价的观测系统(图2)获得相应的地震记录。
根据步骤102中模拟的地震波场传播至地表,应用地表布设的观测排列对波场进行采样处理,分别得到震源排列对应的记录和检波点排列对应的记录。
步骤104,由采集的地震记录,通过单程波外推技术,进行反向外推,即由地表向目标点深度进行波场延拓。聚焦算子采用傅立叶有限差分算子,将步骤103中地表检波点所得共聚焦点道集,作为波场外推的边界,其计算过程完全模拟地震波上行波传播的逆过程。由此,可将震源排列成像分辨率Rs,检波点排列成像分辨率Rd定义如式(1):
Figure BDA0001985523340000051
其中:I表示单位矩阵,S、D分别表示随着深度聚焦后的震源、检波点采样记录,Fs、Fd表示震源、检波点聚焦算子。
步骤105,在深度延拓的过程中,采用叠前深度偏移的零时刻成像原理对波场进行成像,根据下面的等式(2)可获取当前深度的在各个方向的成像分辨率值。
步骤106,根据等式(1)中震源、检波点分辨率的定义,对所有深度进行(5)中的成像处理,获得整个空间上的聚焦分辨率图像。在聚焦过程中,针对前人无垂向分辨率评价的方法,本发明提出垂向、横向分辨率图像提取方法按照下式进
Rxoy(x,y)=U(x,y,z=sz,t=0) (2)
其中:R表示各个平面内的聚焦分辨率图像,U表示在延拓记录,sx,sy,sz表示目标点坐标,XOZ与YOZ为垂向分辨率函数,XOY为横向分辨率函数。
步骤107,通过对震源和检波点排列分辨率图像进行上述处理,可分别获得震源、检波点对目标点成像质量的影响。将震源和检波点排列聚焦图像的乘积作为当前观测系统在该目标区域的总成像分辨率图像。如图3所示,分别为XOY,XOZ和YOZ方向上总成像分辨率图像的结果。
在滚动排列情况中,检波点排列的成像分辨率与炮点排列有关。本发明将滚动排列情况下的检波点排列成像分辨率RD_roll定义如式(3)所示。
Figure BDA0001985523340000061
其中,s为当前炮索引,NumShot为观测系统总炮数,RD为当前炮对应的检波点聚焦束结果。
在实际处理中,当炮点数目过多,会导致计算量过大,本发明创新性的给出以下并行解决方案,如等式(4)中,每一次聚焦过程中,隐式差分计算完成相邻深度记录的更新。其中,p表示不同深度的延拓记录,a、a*是与模型、频率等固有参数的变量。因此,多炮情况下的检波点聚焦过程存在独立性,即实现p1,p2,p3和p4的聚焦计算可并行处理。
Figure BDA0001985523340000062
其中,z表示深度,k表示当前深度,1表示单位1。
步骤108,根据上述分辨率图像的结果进行观测系统排列的调整及优化。一般情况下,分辨率图像能量团越集中,成像分辨率越高。
本发明的基于聚焦分析的叠前偏移垂向分辨率评价方法,是对复杂构造地区,面向特定地下目标的观测系统叠前偏移成像质量和成像精度的评价方法。该方法基于聚焦分析和叠前深度偏移原理,并结合声波波动方程模拟和混合域单程波外推技术,获得当前观测系统情况下,叠前偏移垂向分辨率的计算结果。

Claims (11)

1.基于聚焦分析的叠前偏移垂向分辨率评价方法,其特征在于,该基于聚焦分析的叠前偏移垂向分辨率评价方法包括:
步骤1,建立复杂构造区的速度模型;
步骤2,从目标点激发地震波,使得波场上行至地表并接收;
步骤3,当上行波场达到地表后,应用观测系统获得相应的地震记录;
步骤4,由采集的地震记录,通过单程波外推技术,进行反向外推;
步骤5,采用叠前深度偏移的零时刻成像原理对波场进行成像;
步骤6,获得整个空间上的聚焦分辨率图像;
步骤7,获得震源、检波点对目标点成像质量影响的分析结果;
步骤8,进行观测系统排列的调整及优化。
2.根据权利要求1所述的基于聚焦分析的叠前偏移垂向分辨率评价方法,其特征在于,在步骤1中,根据需要选择地震照明效果不佳的区域作为评价的目标点,收集工区前期地质、地震、测录井和钻井资料,结合地震解释成果和地质认识,建立工区纵波速度模型,并根据勘探要求确定地质目标。
3.根据权利要求1所述的基于聚焦分析的叠前偏移垂向分辨率评价方法,其特征在于,在步骤2中,从目标点激发震源,采用高阶有限差分离散双程声波波动方程进行波场模拟。
4.根据权利要求1所述的基于聚焦分析的叠前偏移垂向分辨率评价方法,其特征在于,在步骤3中,根据步骤2中模拟的地震波场传播至地表,应用地表布设的观测排列对波场进行采样处理,分别得到震源排列对应的记录和检波点排列对应的记录。
5.根据权利要求1所述的基于聚焦分析的叠前偏移垂向分辨率评价方法,其特征在于,在步骤4中,由地表向目标点深度进行波场延拓,聚焦算子采用傅立叶有限差分算子,将步骤3中地表检波点所得共聚焦点道集,作为波场外推的边界,其计算过程完全模拟地震波上行波传播的逆过程,由此,将震源排列成像分辨率Rs,检波点排列成像分辨率Rd定义如式(1):
Figure FDA0001985523330000021
其中:I表示单位矩阵,S、D分别表示随着深度聚焦后的震源、检波点采样记录,Fs、Fd表示震源、检波点聚焦算子。
6.根据权利要求5所述的基于聚焦分析的叠前偏移垂向分辨率评价方法,其特征在于,在步骤5中,在深度延拓的过程中,采用叠前深度偏移的零时刻成像原理对波场进行成像,获取当前深度的在各个方向的成像分辨率值。
7.根据权利要求6所述的基于聚焦分析的叠前偏移垂向分辨率评价方法,其特征在于,在步骤6中,根据等式(1)中震源、检波点分辨率的定义,对所有深度进行步骤5中的成像处理,获得整个空间上的聚焦分辨率图像,在聚焦过程中,针对前人无垂向分辨率评价的方法,垂向、横向分辨率图像提取方法按照下式进行:
Rxoy(x,y)=U(x,y,z=sz,t=0) (2)
其中:R表示各个平面内的聚焦分辨率图像,U表示在延拓记录,sx,sy,sz表示目标点坐标,XOZ与YOZ为垂向分辨率函数,XOY为横向分辨率函数。
8.根据权利要求1所述的基于聚焦分析的叠前偏移垂向分辨率评价方法,其特征在于,在步骤7中,通过对震源和检波点排列分辨率图像进行处理,分别获得震源、检波点对目标点成像质量的影响,将震源和检波点排列聚焦图像的乘积作为当前观测系统在该目标区域的总成像分辨率图像。
9.根据权利要求8所述的基于聚焦分析的叠前偏移垂向分辨率评价方法,其特征在于,在步骤7中,在滚动排列情况中,检波点排列的成像分辨率与炮点排列有关,将滚动排列情况下的检波点排列成像分辨率RD_roll定义如式(3)所示:
Figure FDA0001985523330000031
其中,s为当前炮索引,NumShot为观测系统总炮数,RD为当前炮对应的检波点聚焦束结果。
10.根据权利要求9所述的基于聚焦分析的叠前偏移垂向分辨率评价方法,其特征在于,在步骤7中,在实际处理中,当炮点数目过多,会导致计算量过大,采用等式(4),每一次聚焦过程中,隐式差分计算完成相邻深度记录的更新,其中,p表示不同深度的延拓记录,a、a*是与模型、频率等固有参数的变量;因此,多炮情况下的检波点聚焦过程存在独立性,即实现p1,p2,p3和p4的聚焦计算并行处理;
Figure FDA0001985523330000032
其中,z表示深度,k表示当前深度,1表示单位1。
11.根据权利要求1所述的基于聚焦分析的叠前偏移垂向分辨率评价方法,其特征在于,在步骤8中,根据分辨率图像的结果进行观测系统排列的调整及优化,分辨率图像能量团越集中,成像分辨率越高。
CN201910166749.6A 2019-03-05 2019-03-05 基于聚焦分析的叠前偏移垂向分辨率评价方法 Active CN111665563B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910166749.6A CN111665563B (zh) 2019-03-05 2019-03-05 基于聚焦分析的叠前偏移垂向分辨率评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910166749.6A CN111665563B (zh) 2019-03-05 2019-03-05 基于聚焦分析的叠前偏移垂向分辨率评价方法

Publications (2)

Publication Number Publication Date
CN111665563A true CN111665563A (zh) 2020-09-15
CN111665563B CN111665563B (zh) 2024-02-20

Family

ID=72381663

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910166749.6A Active CN111665563B (zh) 2019-03-05 2019-03-05 基于聚焦分析的叠前偏移垂向分辨率评价方法

Country Status (1)

Country Link
CN (1) CN111665563B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114879257A (zh) * 2022-03-31 2022-08-09 中国石油大学(华东) 一种地震成像分辨分析方法、装置及存储介质
CN115291300A (zh) * 2022-09-28 2022-11-04 山东省鲁南地质工程勘察院(山东省地质矿产勘查开发局第二地质大队) 一种基于复杂地质条件下的三维勘探分析系统及评估方法

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1673774A (zh) * 2004-03-26 2005-09-28 中国石油天然气集团公司 地震观测系统优化设计的层状介质双聚焦方法及其应用
US20090192714A1 (en) * 2008-01-29 2009-07-30 Baker Hughes Incorporated Characterization of Fracture Length and Formation Resistivity from Array Induction Data
CN101937102A (zh) * 2010-01-05 2011-01-05 中国石油化工集团公司 三维观测系统聚焦性能分析方法
CN104142518A (zh) * 2013-10-29 2014-11-12 中国石油化工股份有限公司 一种地震观测系统叠前时间偏移响应的分析方法
CN104570073A (zh) * 2013-10-16 2015-04-29 中国石油化工股份有限公司 一种适用于复杂高陡构造的双反射地震波成像方法
CN104914466A (zh) * 2015-06-26 2015-09-16 中国石油大学(华东) 一种提高地震资料分辨率的方法
US20150331123A1 (en) * 2014-05-15 2015-11-19 Acoustic Zoom, Inc. System and Method for Adaptive Classification and Filtering for Imaging Diffused Energy from the Earths Subsurface
CN106873031A (zh) * 2017-02-15 2017-06-20 中国科学院地质与地球物理研究所 一种三维地震观测系统垂向分辨率定量分析评价方法
US9726771B1 (en) * 2013-07-19 2017-08-08 Z Terra Inc. Time migration diffraction imaging systems and methods
CN108919351A (zh) * 2018-05-15 2018-11-30 中国海洋石油集团有限公司 基于逆时聚焦原理进行观测系统双向聚焦性的评价方法

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1673774A (zh) * 2004-03-26 2005-09-28 中国石油天然气集团公司 地震观测系统优化设计的层状介质双聚焦方法及其应用
US20090192714A1 (en) * 2008-01-29 2009-07-30 Baker Hughes Incorporated Characterization of Fracture Length and Formation Resistivity from Array Induction Data
CN101937102A (zh) * 2010-01-05 2011-01-05 中国石油化工集团公司 三维观测系统聚焦性能分析方法
US9726771B1 (en) * 2013-07-19 2017-08-08 Z Terra Inc. Time migration diffraction imaging systems and methods
CN104570073A (zh) * 2013-10-16 2015-04-29 中国石油化工股份有限公司 一种适用于复杂高陡构造的双反射地震波成像方法
CN104142518A (zh) * 2013-10-29 2014-11-12 中国石油化工股份有限公司 一种地震观测系统叠前时间偏移响应的分析方法
US20150331123A1 (en) * 2014-05-15 2015-11-19 Acoustic Zoom, Inc. System and Method for Adaptive Classification and Filtering for Imaging Diffused Energy from the Earths Subsurface
CN104914466A (zh) * 2015-06-26 2015-09-16 中国石油大学(华东) 一种提高地震资料分辨率的方法
CN106873031A (zh) * 2017-02-15 2017-06-20 中国科学院地质与地球物理研究所 一种三维地震观测系统垂向分辨率定量分析评价方法
CN108919351A (zh) * 2018-05-15 2018-11-30 中国海洋石油集团有限公司 基于逆时聚焦原理进行观测系统双向聚焦性的评价方法

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114879257A (zh) * 2022-03-31 2022-08-09 中国石油大学(华东) 一种地震成像分辨分析方法、装置及存储介质
CN114879257B (zh) * 2022-03-31 2023-10-20 中国石油大学(华东) 一种地震成像分辨分析方法、装置及存储介质
CN115291300A (zh) * 2022-09-28 2022-11-04 山东省鲁南地质工程勘察院(山东省地质矿产勘查开发局第二地质大队) 一种基于复杂地质条件下的三维勘探分析系统及评估方法

Also Published As

Publication number Publication date
CN111665563B (zh) 2024-02-20

Similar Documents

Publication Publication Date Title
CN108802813B (zh) 一种多分量地震资料偏移成像方法及系统
CN106405651B (zh) 一种基于测井匹配的全波形反演初始速度模型构建方法
CN106443765B (zh) 基于多分量观测系统的城市工程地震探测综合成像方法
US7830748B2 (en) Method for acoustic imaging of the earth's subsurface using a fixed position sensor array and beam steering
CN111158049B (zh) 一种基于散射积分法的地震逆时偏移成像方法
CN109557582B (zh) 一种二维多分量地震资料偏移成像方法及系统
WO2012139082A1 (en) Event selection in the image domain
CN111722284B (zh) 一种基于道集数据建立速度深度模型的方法
CN111352151B (zh) 一种井中虚源法
CN109655918B (zh) 地面浅井微地震监测观测台站位置确定方法及系统
CN111665563A (zh) 基于聚焦分析的叠前偏移垂向分辨率评价方法
CN107656308B (zh) 一种基于时间深度扫描的共散射点叠前时间偏移成像方法
CN104316961A (zh) 获取风化层的地质参数的方法
CN112230274B (zh) 面向随钻导向的声波方程频率域逆时偏移快速成像方法
CN110579799B (zh) 一种等旅行时间间隔的地震采集观测方法及系统
CN108919351A (zh) 基于逆时聚焦原理进行观测系统双向聚焦性的评价方法
CN112748463A (zh) 一种基于深度学习照明分析的局部偏移成像方法
CN108375794B (zh) 基于对称观测的vsp缝洞绕射成像技术方法
US6662111B2 (en) Method for analyzing reflection curvature in seismic data volumes
CN112305615B (zh) 一种地震资料角度域共成像点道集提取方法及系统
Kazei et al. Acquisition and near-surface impacts on VSP mini-batch FWI and RTM imaging in desert environment
CN112462428B (zh) 一种多分量地震资料偏移成像方法及系统
CN115993650B (zh) 一种基于棱柱波的地震干涉成像方法
Rollins et al. Calibrating 3D VSP image area with finite difference modeling–A case study in Gulf of Mexico
CN113835123B (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