CN106873031B - 一种三维地震观测系统垂向分辨率定量分析评价方法 - Google Patents

一种三维地震观测系统垂向分辨率定量分析评价方法 Download PDF

Info

Publication number
CN106873031B
CN106873031B CN201710081114.7A CN201710081114A CN106873031B CN 106873031 B CN106873031 B CN 106873031B CN 201710081114 A CN201710081114 A CN 201710081114A CN 106873031 B CN106873031 B CN 106873031B
Authority
CN
China
Prior art keywords
observation system
seismic
main lobe
seismic observation
geophone station
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
Application number
CN201710081114.7A
Other languages
English (en)
Other versions
CN106873031A (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.)
Institute of Geology and Geophysics of CAS
Original Assignee
Institute of Geology and Geophysics of CAS
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 Institute of Geology and Geophysics of CAS filed Critical Institute of Geology and Geophysics of CAS
Priority to CN201710081114.7A priority Critical patent/CN106873031B/zh
Publication of CN106873031A publication Critical patent/CN106873031A/zh
Application granted granted Critical
Publication of CN106873031B publication Critical patent/CN106873031B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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/30Analysis

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)根据沿垂向的共聚焦结果计算主瓣宽度;5)根据得到的主瓣宽度计算主瓣占幅比;6)根据得到的主瓣宽度以及主瓣占幅比,对该地震观测系统的地震偏移成像进行分析。本发明可以广泛应用于地震观测系统的设计分析中。

Description

一种三维地震观测系统垂向分辨率定量分析评价方法
技术领域
本发明涉及一种地震勘探的三维地震观测系统的设计与评价方法,尤其涉及一种三维地震观测系统垂向分辨率定量分析评价方法。
背景技术
随着地震勘探技术的发展,地震勘探逐步从高密度向高精度勘探转变,而地震观测系统是影响地震成像的关键性因素之一,对于地震观测系统的设计要求也越来越高。地震成像的质量一般用地震分辨率来描述。地震分辨率指分辨两个相邻同相轴的能力,其在时间域和空间域上又分为时间分辨率和空间分辨率。其中,空间域的分辨率又根据方向的不同分为水平分辨率与垂向分辨率。常规的垂向分辨率分析是基于几何地震学和无限频带宽度的假设,其计算方法一般是分析地震子波形态,以雷克子波为例,通过求取雷克子波相邻两拐点之间的距离,来定义时间域上该子波所能分辨的最小地层厚度。而实际地震勘探是有限频带宽度并且受空间离散采样的影响,所以导致基于几何地震学原理的常规观测系统设计,无法适用于实际的复杂地质情况,也无法直观反映最终地震偏移成像效果的优劣。
共聚焦(Common Focus Point,简称CFP)方法是一种直接将地震偏移方法应用于地震观测系统设计评价的新方法。该方法首先通过对检波点与震源点分别进行波场延拓以及聚焦运算,得到检波点聚焦矩阵与震源点聚焦矩阵。然后在频率空间域将检波点聚焦矩阵与震源点聚焦矩阵进行乘积运算,得到三维观测系统的分辨率矩阵。基于常规垂向分辨率分析的地震观测系统设计方法既没有考虑地下介质的影响,也没有考虑观测系统对地震成像的影响,所以并不能直接反映实际的地震成像情况。通过基于地震波场延拓的共聚焦分析技术,可以综合考虑地下介质、观测系统、频带宽度和主频等因素对地震成像分辨率的影响。
但是,目前共聚焦分析计算地震分辨率主要是运用于水平分辨率的计算,而水平分辨率和垂向分辨率的计算所依据的物理原理并不一致。水平分辨率主要是依据菲涅尔带对地震成像分辨率的控制,而垂向分辨率主要是依据能区分相邻地震子波的最小间隔的物理原理。所以使用共聚焦分析来计算水平分辨率的方法并不能直接用于对垂向分辨率的计算。
发明内容
针对上述问题,本发明的目的是提供一种三维地震观测系统垂向分辨率定量分析评价方法,在给定观测系统参数的情况下,可以计算出地下任意目标点处的地震垂向分辨率,从而实现了对观测系统成像特性的高精度定量分析。
为实现上述目的,本发明采取以下技术方案:一种三维地震观测系统垂向分辨率定量分析评价方法,其特征在于包括以下步骤:1)给定复杂速度模型、待评价地震观测系统和地下目标点位置信息,在复杂地质条件下对地震观测系统进行波场延拓,得到地下目标点传播到地面各检波点的频率域波场信息;2)采用共聚焦快速计算方法分别对地面各检波点的频率域波场信息进行垂向共聚焦分析,得到各检波点在地下目标点处的共聚焦结果;3)将步骤2)中得到的各检波点的共聚焦结果进行综合,也即相乘再相加,得到整个地震观测系统的共聚焦结果;4)对整个观测系统的共聚焦结果进行定量化,即根据沿垂向的共聚焦结果计算主瓣宽度,用于表征地震偏移成像的垂向分辨率,得到复杂介质条件下待评价地震观测系统对目标点成像的垂向分辨率;5)根据得到的主瓣宽度计算主瓣占幅比,用于表征地震偏移成像的垂向清晰度,得到复杂介质条件下待评价地震观测系统对目标点成像的垂向清晰度;6)根据得到的该地震观测系统在地下目标点处的垂向共聚焦结果,也即主瓣宽度以及主瓣占幅比,对该地震观测系统的地震偏移成像进行分析。
所述步骤2)中,对各检波点进行垂向共聚焦分析,得到各检波点在地下目标点处的共聚焦结果为:
b(r0,ω)=∫f(r,ω)a(r0,ω)dr,
其中,f(r,ω)是聚集算子,a(r0,ω)是检波点位置处记录到的频率域波场,r是地下目标点的位置信息,r0是地面检波点的位置信息,ω是频率信息。
所述步骤4)中,主瓣宽度的计算方法包括以下步骤:
①从整个地震观测系统的共聚焦结果中提取过目标点的沿垂向的数值切片;
②对该数值切片进行二阶求导,找到目标点两侧的两拐点,即二阶导数为零的两个坐标点;
③根据两个拐点坐标计算得到主瓣宽度,也就是垂向分辨率,计算公式为:
d=|xa-xb|,
式中,xa、xb为两拐点的横坐标,也即b″(xa)=0、b″(xb)=0。
所述步骤5)中,主瓣占幅比是指主瓣宽度内所有点的振幅值的平方和比上整个数值切片内所有点的振幅值的平方和。
所述主瓣占幅比的计算公式为:
式中,A(xi)为数值切片内各点的振幅值,xi是数值切片的横坐标,i=1,2,...,n,n为该数值切片内的点的总数。
所述步骤6)中,对该地震观测系统的地震偏移成像进行分析的方法为:主瓣宽度数值越小,说明该地震观测系统的地震成像能分辨的最小间距越小,其地震偏移成像效果越好;否则,效果越不好;主瓣占幅比反映的是成像清晰度,主瓣占幅比数值越大,说明该地震观测系统的地震成像清晰度越高,效果越好;否则,效果越不好。
本发明由于采取以上技术方案,其具有以下优点:1、本发明由于通过定义主瓣宽度和主瓣占幅比,将共聚焦结果直接用于对地震垂向分辨率的描述,这样得到的结果不仅更加接近于实际的地震成像精度,而且还能够定量分析地震观测系统对地震成像的影响,从而能够指导地震采集设计的优化,对于油气勘探开发具有很重要的意义。2、本发明实现了基于波场延拓的共聚焦快速计算方法,对常规共聚焦算法进行改进,常规算法是对地震观测系统中的检波点波场和目标点波场分别进行延拓和聚焦分析,而本发明共聚焦快速计算方法只对检波点波场进行延拓和聚焦分析,计算速度快,有利的促进了共聚焦技术在实际地震勘探中的应用,对指导地震观测系统设计,提高油气勘探效率有极大的帮助。本发明可以广泛应用于地震观测系统的垂向分辨率定量分析领域。
附图说明
图1为本发明的方法流程图;
图2为均匀介质下某观测系统的共聚焦结果垂向切片,其中实线为共聚焦结果,虚线为它的二次导数值,其两个相邻零点之间的宽度为主瓣宽度,即图中的虚线轴上的黑色实线;
图3为某复杂模型下,本发明得到的共聚焦结果(a)和它的垂向切片及二次导数值(b)。
具体实施方式
下面结合附图和实施例对本发明进行详细的描述。
如图1所示,本发明三维地震观测系统垂向分辨率定量分析评价方法,包括以下步骤:
1)给定复杂速度模型、待评价地震观测系统和地下目标点位置信息,在复杂地质条件下对待评价地震观测系统进行波场延拓,得到地下目标点传播到地面各检波点的频率域波场信息。
其中,复杂速度模型和待评价地震观测系统均为已有技术,本发明不再赘述。
2)采用共聚焦快速计算方法分别对地面各检波点的频率域波场信息进行垂向共聚焦分析,得到各检波点在地下目标点处的共聚焦结果。
与水平共聚焦分析不同,垂向共聚焦结果能够直接反映偏移后目标点处垂直切片上的真实地震子波。在理想情况下,垂向共聚焦结果仅含1个非零值,但实际情况下,由于垂向共聚焦结果受空间采样,地震子波及地下介质的影响,可能会得到非常复杂的结果。所以,通过基于波动方程的延拓方法,模拟得到复杂介质下的频率域波场,由于聚焦为波场逆传播的过程,所以通过求取传播算子的共轭得到聚集算子,将聚集算子运用于记录到的波场即可以得到复杂介质下垂向共聚焦分析的结果。
根据上述分析,得到对各检波点的频率域波场进行垂向共聚焦分析,得到各检波点在地下目标点处的共聚焦结果为:
b(r0,ω)=∫f(r,ω)a(r0,ω)dr (1)
其中,f(r,ω)是聚集算子,a(r0,ω)是检波点位置处记录到的频率域波场,r是地下目标点的位置信息,r0是地面检波点的位置信息,ω是频率信息。
3)将步骤2)中得到的各检波点的共聚焦结果进行综合,也即相乘再相加,得到整个地震观测系统的共聚焦结果。
4)对整个观测系统的共聚焦结果进行定量化,即根据沿垂向的共聚焦结果计算主瓣宽度,用于表征地震偏移成像的垂向分辨率,得到复杂介质条件下待评价地震观测系统对目标点成像的垂向分辨率。
如图2所示,主瓣宽度的计算方法包括以下步骤:
①从整个地震观测系统的共聚焦结果中提取过目标点的沿垂向的数值切片。
②对该数值切片进行二阶求导,找到目标点两侧的两拐点,即二阶导数为零的两个坐标点。
③根据两个拐点坐标计算得到主瓣宽度,也就是垂向分辨率,计算公式为:
d=|xa-xb| (3)
式中,xa、xb为两拐点的横坐标,也即b″(xa)=0、b″(xb)=0。
5)根据得到的主瓣宽度计算主瓣占幅比,用于表征地震偏移成像的垂向清晰度,得到复杂介质条件下待评价地震观测系统对目标点成像的垂向清晰度。
主瓣宽度内所有点的振幅值的平方和比上整个数值切片内所有点的振幅值的平方和就是主瓣占幅比,它能够反映地震成像的垂向清晰度,计算公式为:
式中,A(xi)为数值切片内各点的振幅值,xi是数值切片的横坐标,i=1,2,...,n,n为该数值切片内的点的总数。
如图3所示,为某三维观测系统在复杂介质情况下的共聚焦分辨率计算结果,从图中可以看出,应用本发明可以实现复杂介质下的垂向分辨率的准确计算,进而描述地震观测系统的成像特性,对采集设计有着重要的指导意义。
6)根据得到的该地震观测系统在地下目标点处的垂向共聚焦结果,也即主瓣宽度以及主瓣占幅比,对该地震观测系统的地震偏移成像进行分析。
根据得到的主瓣宽度和主瓣占幅比对该地震观测系统的地震偏移成像进行分析。其中,主瓣宽度反映的是地震观测系统在目标点处的垂向成像分辨率,分辨率数值越小,说明地震成像能分辨的最小间距越小,其效果越好;否则,效果越不好。主瓣占幅比反映的是成像清晰度,主瓣占幅比数值越大,说明地震成像清晰度越高,效果越好;否则,效果越不好。也就是说,通过本发明,可以定量的得到地震观测系统在复杂介质情况下,目标点处的成像特性,并以此来分析和指导采集系统的设计。
上述各实施例仅用于说明本发明,其中各部件的结构、连接方式和制作工艺等都是可以有所变化的,凡是在本发明技术方案的基础上进行的等同变换和改进,均不应排除在本发明的保护范围之外。

Claims (6)

1.一种三维地震观测系统垂向分辨率定量分析评价方法,其特征在于包括以下步骤:
1)给定复杂速度模型、待评价地震观测系统和地下目标点位置信息,在复杂地质条件下对地震观测系统进行波场延拓,得到地下目标点传播到地面各检波点的频率域波场信息;
2)采用共聚焦快速计算方法分别对地面各检波点的频率域波场信息进行垂向共聚焦分析,得到各检波点在地下目标点处的共聚焦结果;
3)将步骤2)中得到的各检波点的共聚焦结果进行综合,也即相乘再相加,得到整个地震观测系统的共聚焦结果;
4)对整个观测系统的共聚焦结果进行定量化,即根据沿垂向的共聚焦结果计算主瓣宽度,用于表征地震偏移成像的垂向分辨率,得到复杂介质条件下待评价地震观测系统对目标点成像的垂向分辨率;
5)根据得到的主瓣宽度计算主瓣占幅比,用于表征地震偏移成像的垂向清晰度,得到复杂介质条件下待评价地震观测系统对目标点成像的垂向清晰度;
6)根据得到的该地震观测系统在地下目标点处的垂向共聚焦结果,也即主瓣宽度以及主瓣占幅比,对该地震观测系统的地震偏移成像进行分析。
2.如权利要求1所述的一种三维地震观测系统垂向分辨率定量分析评价方法,其特征在于:所述步骤2)中,对各检波点进行垂向共聚焦分析,得到各检波点在地下目标点处的共聚焦结果为:
b(r0,ω)=∫f(r,ω)a(r0,ω)dr,
其中,f(r,ω)是聚集算子,a(r0,ω)是检波点位置处记录到的频率域波场,r是地下目标点的位置信息,r0是地面检波点的位置信息,ω是频率信息。
3.如权利要求1所述的一种三维地震观测系统垂向分辨率定量分析评价方法,其特征在于:所述步骤4)中,主瓣宽度的计算方法包括以下步骤:
①从整个地震观测系统的共聚焦结果中提取过目标点的沿垂向的数值切片;
②对该数值切片进行二阶求导,找到目标点两侧的两拐点,即二阶导数为零的两个坐标点;
③根据两个拐点坐标计算得到主瓣宽度,也就是垂向分辨率,计算公式为:
d=|xa-xb|,
式中,xa、xb为两拐点的横坐标,也即b”(xa)=0、b”(xb)=0,其中,b”()是指对数值切片中两点xa、xb求二阶导数。
4.如权利要求1所述的一种三维地震观测系统垂向分辨率定量分析评价方法,其特征在于:所述步骤5)中,主瓣占幅比是指主瓣宽度内所有点的振幅值的平方和比上整个数值切片内所有点的振幅值的平方和。
5.如权利要求4所述的一种三维地震观测系统垂向分辨率定量分析评价方法,其特征在于:所述主瓣占幅比的计算公式为:
式中,A(xi)为数值切片内各点的振幅值,xi是数值切片的横坐标,i=1,2,...,n,n为该数值切片内的点的总数。
6.如权利要求1所述的一种三维地震观测系统垂向分辨率定量分析评价方法,其特征在于:所述步骤6)中,对该地震观测系统的地震偏移成像进行分析的方法为:
主瓣宽度数值越小,说明该地震观测系统的地震成像能分辨的最小间距越小,其地震偏移成像效果越好;否则,效果越不好;
主瓣占幅比反映的是成像清晰度,主瓣占幅比数值越大,说明该地震观测系统的地震成像清晰度越高,效果越好;否则,效果越不好。
CN201710081114.7A 2017-02-15 2017-02-15 一种三维地震观测系统垂向分辨率定量分析评价方法 Active CN106873031B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710081114.7A CN106873031B (zh) 2017-02-15 2017-02-15 一种三维地震观测系统垂向分辨率定量分析评价方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710081114.7A CN106873031B (zh) 2017-02-15 2017-02-15 一种三维地震观测系统垂向分辨率定量分析评价方法

Publications (2)

Publication Number Publication Date
CN106873031A CN106873031A (zh) 2017-06-20
CN106873031B true CN106873031B (zh) 2019-01-15

Family

ID=59167388

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710081114.7A Active CN106873031B (zh) 2017-02-15 2017-02-15 一种三维地震观测系统垂向分辨率定量分析评价方法

Country Status (1)

Country Link
CN (1) CN106873031B (zh)

Families Citing this family (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN107728209B (zh) * 2017-10-13 2019-11-05 中国地质大学(北京) 一种确定人工地震记录垂向分辨率的方法
CN111665563B (zh) * 2019-03-05 2024-02-20 中石化石油工程技术服务有限公司 基于聚焦分析的叠前偏移垂向分辨率评价方法
CN111007571B (zh) * 2019-11-28 2022-05-13 吉林大学 一种基于三维构造张量的航磁数据地质体边界识别方法
CN114442145A (zh) * 2022-01-18 2022-05-06 中国石油大学(华东) 一种目标导向地震采集方法、装置、介质及设备
CN114879257B (zh) * 2022-03-31 2023-10-20 中国石油大学(华东) 一种地震成像分辨分析方法、装置及存储介质

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0610123A1 (fr) * 1993-02-01 1994-08-10 Compagnie Generale De Geophysique Procédé pour améliorer l'estimation, par analyse de focalisation, des vitesses de propagation des ondes sismiques
CN101105537A (zh) * 2006-07-12 2008-01-16 中国石油集团东方地球物理勘探有限责任公司 一种高精度的深度域叠前地震数据反演方法
CN102062872A (zh) * 2009-11-11 2011-05-18 中国石油天然气集团公司 可控震源多因子扫描方法
CN102565852A (zh) * 2010-12-08 2012-07-11 同济大学 针对储层含油气性检测的角度域叠前偏移数据处理方法
CN105487118A (zh) * 2014-09-28 2016-04-13 中国石油化工股份有限公司 一种处理叠前宽线地震数据的方法
CN106154354A (zh) * 2015-04-22 2016-11-23 中国石油化工股份有限公司 一种成像道集的应用方法

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP0610123A1 (fr) * 1993-02-01 1994-08-10 Compagnie Generale De Geophysique Procédé pour améliorer l'estimation, par analyse de focalisation, des vitesses de propagation des ondes sismiques
CN101105537A (zh) * 2006-07-12 2008-01-16 中国石油集团东方地球物理勘探有限责任公司 一种高精度的深度域叠前地震数据反演方法
CN102062872A (zh) * 2009-11-11 2011-05-18 中国石油天然气集团公司 可控震源多因子扫描方法
CN102565852A (zh) * 2010-12-08 2012-07-11 同济大学 针对储层含油气性检测的角度域叠前偏移数据处理方法
CN105487118A (zh) * 2014-09-28 2016-04-13 中国石油化工股份有限公司 一种处理叠前宽线地震数据的方法
CN106154354A (zh) * 2015-04-22 2016-11-23 中国石油化工股份有限公司 一种成像道集的应用方法

Also Published As

Publication number Publication date
CN106873031A (zh) 2017-06-20

Similar Documents

Publication Publication Date Title
CN106873031B (zh) 一种三维地震观测系统垂向分辨率定量分析评价方法
CN108037526B (zh) 基于全波波场vsp/rvsp地震资料的逆时偏移方法
CN108957549B (zh) 一种辫状河沉积非均质致密砂岩气藏地质建模方法
CN104280777B (zh) 一种压制陆上地震资料多次波干扰的方法
CN102645670B (zh) 一种基于叠加响应分析的观测系统优化设计方法
CN110261904B (zh) 基于广义s变换的近地表q值反演及分类评价方法
AU2017366758A1 (en) Diving wave illumination using migration gathers
CN101937102A (zh) 三维观测系统聚焦性能分析方法
Yuan et al. 6D phase-difference attributes for wide-azimuth seismic data interpretation
CN113433547A (zh) 一种探地雷达隐伏裂缝偏移成像方法、系统、终端及介质
CN112255690A (zh) 基于地震相位分解的自适应围岩强反射分离方法
CN112230274B (zh) 面向随钻导向的声波方程频率域逆时偏移快速成像方法
CN105093300B (zh) 一种地质体边界识别方法及装置
CN110673210B (zh) 一种地震原始数据信噪比定量分析评价方法
CN110954958A (zh) 裂缝与断层的预测方法及系统
CN108919351A (zh) 基于逆时聚焦原理进行观测系统双向聚焦性的评价方法
CN108562936A (zh) 一种裂缝预测方法、系统、存储介质及终端
CN110244383B (zh) 基于近地表数据的地质岩性综合模型创建方法
CN114428343A (zh) 基于归一化互相关的Marchenko成像方法及系统
CN111665563B (zh) 基于聚焦分析的叠前偏移垂向分辨率评价方法
CN112946742B (zh) 一种拾取精确叠加速度谱的方法
CN116224324A (zh) 基于深度学习的超分辨率3d-gpr图像的频率-波数分析方法
CN102830424A (zh) 一种检波器组合参数计算方法
CN112444859A (zh) 协同异变蚂蚁体的页岩储层断裂识别方法及系统
CN108375794A (zh) 基于对称观测的vsp缝洞绕射成像技术方法

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