CN112147699A - 结合过井剖面的vsp深度域层速度模型构建方法及系统 - Google Patents
结合过井剖面的vsp深度域层速度模型构建方法及系统 Download PDFInfo
- Publication number
- CN112147699A CN112147699A CN201910579729.1A CN201910579729A CN112147699A CN 112147699 A CN112147699 A CN 112147699A CN 201910579729 A CN201910579729 A CN 201910579729A CN 112147699 A CN112147699 A CN 112147699A
- Authority
- CN
- China
- Prior art keywords
- vsp
- layer
- well
- determining
- time domain
- 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
- 238000010276 construction Methods 0.000 title claims description 12
- 238000000034 method Methods 0.000 claims abstract description 30
- 238000006243 chemical reaction Methods 0.000 claims abstract description 14
- 230000015572 biosynthetic process Effects 0.000 claims description 14
- 238000010586 diagram Methods 0.000 description 12
- 238000012545 processing Methods 0.000 description 11
- 230000005012 migration Effects 0.000 description 7
- 238000013508 migration Methods 0.000 description 7
- 230000008901 benefit Effects 0.000 description 4
- 238000011160 research Methods 0.000 description 3
- 230000008859 change Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000008569 process Effects 0.000 description 1
- 238000003672 processing method Methods 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/40—Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
-
- 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
- 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
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Remote Sensing (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
公开了一种结合过井剖面的VSP深度域层速度模型构建方法及系统。该方法包括:根据零井源距VSP资料,计算VSP井位置处的地下地层的深度域层速度进行时深转换,获得VSP井位置处的地下地层的时间域层速度;根据过井剖面的地震层位,确定主要反射层;根据过井剖面的反射同相轴,确定每个主要反射层的次级反射层;根据VSP井位置处的地下地层的时间域层速度,确定对应的主要反射层与次要反射层的层速度,确定VSP时间域层速度模型;将VSP时间域层速度模型进行时深转换,获得VSP深度域层速度模型。本发明通过过井剖面构建VSP深度域层速度模型,提高VSP速度模型的准确度,消除由常规VSP速度建模方法所带来的局限性。
Description
技术领域
本发明涉及地球物理勘探领域,更具体地,涉及一种结合过井剖面的VSP深度域层速度模型构建方法及系统。
背景技术
随着非零偏VSP勘探技术在地震勘探领域的推广应用,有关非零偏VSP地震资料处理技术的研究也越来越受到业界的重视,目前,对非零偏VSP资料的处理主要还是采用常规的VSP-CDP转换以及Kirchhoff积分偏移,然而,不管是VSP-CDP转换还是Kirchhoff积分偏移,都要求事先建立速度模型,只有在已知速度模型的前提下,才有可能接下来进行VSP-CDP转换或Kirchhoff积分偏移。目前,建立VSP速度模型主要还是通过直接利用零井源距VSP初至时间求取地层层速度的方式来建立,所以,建立的VSP速度模型其实为水平层速度模型,这样的速度模型往往导致VSP-CDP转换以及Kirchhoff积分偏移存在以下局限性:1、VSP-CDP转换仅对水平层状介质或横向变速缓慢的介质才能取得比较好的处理结果;2、利用Kirchhoff积分偏移方法对VSP数据成像时,要求速度场横向变化不能太大。换句话说,按照常规的VSP速度模型构建方法,只有当横向速度变化较为缓慢的情况下,才可以得到比较好的处理结果,如果地层横向速度变化较大,地下地质产状不是水平产状,则很难得到高质量的非零偏VSP资料处理结果,因此,为了使构建的VSP速度模型更符合地下的实际构造情况,有必要寻找一种更有效的VSP速度模型的构建方法,在过去的这些年内,不少从事非零偏VSP资料处理方法研究的科研人员为此做了大量的工作,也取得了不少研究成果,例如通过利用VSP旅行时进行层析反演来得到速度模型的层速度,但是该方法由于受到VSP旅行时信息量较少的限制,所以导致反演得到的层速度不是很理想。因此,有必要开发一种结合过井剖面的VSP深度域层速度模型构建方法及系统。
公开于本发明背景技术部分的信息仅仅旨在加深对本发明的一般背景技术的理解,而不应当被视为承认或以任何形式暗示该信息构成已为本领域技术人员所公知的现有技术。
发明内容
本发明提出了一种结合过井剖面的VSP深度域层速度模型构建方法及系统,其能够通过过井剖面构建VSP深度域层速度模型,提高VSP速度模型的准确度,消除由常规VSP速度建模方法所带来的局限性。
根据本发明的一方面,提出了一种结合过井剖面的VSP深度域层速度模型构建方法。所述方法可以包括:根据零井源距VSP资料,计算VSP井位置处的地下地层的深度域层速度进而进行时深转换,获得VSP井位置处的地下地层的时间域层速度;根据过井剖面的地震层位,确定VSP时间域层速度模型的主要反射层;根据所述过井剖面的反射同相轴,确定所述VSP时间域层速度模型的每个主要反射层的次级反射层;根据所述VSP井位置处的地下地层的时间域层速度,确定对应的所述主要反射层与所述次要反射层的层速度,进而确定VSP时间域层速度模型;将所述VSP时间域层速度模型的时间域层速度转换为深度域层速度,获得VSP深度域层速度模型。
优选地,所述根据过井剖面的地震层位,确定主要反射层包括:根据所述VSP时间域层速度模型的空间范围,确定所述过井剖面对应的范围;在所述范围内的过井剖面上确定主要地震层位,以VSP井所在位置为起点,对每个主要地震层位的主同相轴进行追踪,确定所述主同相轴上每个共深度点所对应的时间;根据所述每个共深度点所对应的时间,确定所述VSP时间域层速度模型的主要反射层在每个共深度点所对应的时间位置。
优选地,所述根据所述过井剖面的反射同相轴,确定每个主要反射层的次级反射层包括:判断所述过井剖面的反射层内是否出现反射同相轴;若是,则该反射层内存在次级反射层,以所述VSP井所在位置为起点,对所述反射同相轴进行追踪,确定所述反射同相轴在所述VSP时间域层速度模型上每个共深度点所对应的时间,进而确定所述次级反射层在每个共深度点所对应的时间位置;若否,则该反射层内不存在次级反射层。
优选地,所述根据所述VSP井位置处的地下地层的时间域层速度,确定对应的所述主要反射层与所述次要反射层的层速度为:将反射层在VSP井位置处所对应的地层的时间域层速度作为该反射层的层速度。
优选地,所述根据所述VSP井位置处的地下地层的时间域层速度,确定对应的所述主要反射层与所述次要反射层的层速度为:若反射层的厚度大于VSP的层厚度,则计算所述反射层对应的多个时间域层速度的平均速度,作为所述反射层的层速度。
根据本发明的另一方面,提出了一种结合过井剖面的VSP深度域层速度模型构建系统,其特征在于,该系统包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:根据零井源距VSP资料,计算VSP井位置处的地下地层的深度域层速度进而进行时深转换,获得VSP井位置处的地下地层的时间域层速度;根据过井剖面的地震层位,确定VSP时间域层速度模型的主要反射层;根据所述过井剖面的反射同相轴,确定所述VSP时间域层速度模型的每个主要反射层的次级反射层;根据所述VSP井位置处的地下地层的时间域层速度,确定对应的所述主要反射层与所述次要反射层的层速度,进而确定VSP时间域层速度模型;将所述VSP时间域层速度模型的时间域层速度转换为深度域层速度,获得VSP深度域层速度模型。
优选地,所述根据过井剖面的地震层位,确定主要反射层包括:根据所述VSP时间域层速度模型的空间范围,确定所述过井剖面对应的范围;在所述范围内的过井剖面上确定主要地震层位,以VSP井所在位置为起点,对每个主要地震层位的主同相轴进行追踪,确定所述主同相轴上每个共深度点所对应的时间;根据所述每个共深度点所对应的时间,确定所述VSP时间域层速度模型的主要反射层在每个共深度点所对应的时间位置。
优选地,所述根据所述过井剖面的反射同相轴,确定每个主要反射层的次级反射层包括:判断所述过井剖面的反射层内是否出现反射同相轴;若是,则该反射层内存在次级反射层,以所述VSP井所在位置为起点,对所述反射同相轴进行追踪,确定所述反射同相轴在所述VSP时间域层速度模型上每个共深度点所对应的时间,进而确定所述次级反射层在每个共深度点所对应的时间位置;若否,则该反射层内不存在次级反射层。
优选地,所述根据所述VSP井位置处的地下地层的时间域层速度,确定对应的所述主要反射层与所述次要反射层的层速度为:将反射层在VSP井位置处所对应的地层的时间域层速度作为该反射层的层速度。
优选地,所述根据所述VSP井位置处的地下地层的时间域层速度,确定对应的所述主要反射层与所述次要反射层的层速度为:若反射层的厚度大于VSP的层厚度,则计算所述反射层对应的多个时间域层速度的平均速度,作为所述反射层的层速度。
本发明的方法和装置具有其它的特性和优点,这些特性和优点从并入本文中的附图和随后的具体实施方式中将是显而易见的,或者将在并入本文中的附图和随后的具体实施方式中进行详细陈述,这些附图和具体实施方式共同用于解释本发明的特定原理。
附图说明
通过结合附图对本发明示例性实施例进行更详细的描述,本发明的上述以及其它目的、特征和优势将变得更加明显,其中,在本发明示例性实施例中,相同的参考标号通常代表相同部件。
图1示出了根据本发明的结合过井剖面的VSP深度域层速度模型构建方法的步骤的流程图。
图2a和图2b分别示出了根据本发明的一个实施例的零井源距VSP资料的X分量和Z分量的示意图。
图3示出了根据本发明的一个实施例的时间域纵波速度与横波速度的示意图。
图4示出了根据本发明的一个实施例的过井剖面的示意图。
图5示出了根据本发明的一个实施例的VSP深度域层速度模型的示意图。
图6a和图6b分别示出了根据本发明的一个实施例的VSP深度域纵波、横波偏移剖面的示意图。
图7a和图7b分别示出了根据本发明的一个实施例的VSP时间域纵波、转换波偏移剖面的示意图。
图8a、图8b分别示出了根据常规VSP速度建模得到的VSP时间域纵波、转换波偏移剖面的示意图。
具体实施方式
下面将参照附图更详细地描述本发明。虽然附图中显示了本发明的优选实施例,然而应该理解,可以以各种形式实现本发明而不应被这里阐述的实施例所限制。相反,提供这些实施例是为了使本发明更加透彻和完整,并且能够将本发明的范围完整地传达给本领域的技术人员。
图1示出了根据本发明的结合过井剖面的VSP深度域层速度模型构建方法的步骤的流程图。
在该实施例中,根据本发明的结合过井剖面的VSP深度域层速度模型构建方法可以包括:步骤101,根据零井源距VSP资料,计算VSP井位置处的地下地层的深度域层速度,进而进行时深转换,获得VSP井位置处的地下地层的时间域层速度;步骤102,根据过井剖面的地震层位,确定VSP时间域层速度模型的主要反射层;步骤103,根据过井剖面的反射同相轴,确定VSP时间域层速度模型的每个主要反射层的次级反射层;步骤104,根据VSP井位置处的地下地层的时间域层速度,确定对应的主要反射层与次要反射层的层速度,进而确定VSP时间域层速度模型;步骤105,将VSP时间域层速度模型的时间域层速度转换为深度域层速度,获得VSP深度域层速度模型。
在一个示例中,根据过井剖面的地震层位,确定主要反射层包括:根据VSP时间域层速度模型的空间范围,确定过井剖面对应的范围;在范围内的过井剖面上确定主要地震层位,以VSP井所在位置为起点,对每个主要地震层位的主同相轴进行追踪,确定主同相轴上每个共深度点所对应的时间;根据每个共深度点所对应的时间,确定VSP时间域层速度模型的主要反射层在每个共深度点所对应的时间位置。
在一个示例中,根据过井剖面的反射同相轴,确定每个主要反射层的次级反射层包括:判断过井剖面的反射层内是否出现反射同相轴;若是,则该反射层内存在次级反射层,以VSP井所在位置为起点,对反射同相轴进行追踪,确定反射同相轴在VSP时间域层速度模型上每个共深度点所对应的时间,进而确定次级反射层在每个共深度点所对应的时间位置;若否,则该反射层内不存在次级反射层。
在一个示例中,根据VSP井位置处的地下地层的时间域层速度,确定对应的主要反射层与次要反射层的层速度为:将反射层在VSP井位置处所对应的地层的时间域层速度作为该反射层的层速度。
在一个示例中,根据VSP井位置处的地下地层的时间域层速度,确定对应的主要反射层与次要反射层的层速度为:若反射层的厚度大于VSP的层厚度,则计算反射层对应的多个时间域层速度的平均速度,作为反射层的层速度。
具体地,根据本发明的结合过井剖面的VSP深度域层速度模型构建方法可以包括:
根据零井源距VSP资料,计算VSP井位置处的地下地层的深度域层速度进而进行时深转换,获得VSP井位置处的地下地层的时间域层速度。
根据VSP时间域层速度模型的空间范围,确定过井剖面对应的范围;在范围内的过井剖面上确定主要地震层位,以VSP井所在位置为起点,对每个主要地震层位的主同相轴进行追踪,确定主同相轴在VSP时间域层速度模型上每个共深度点所对应的时间;根据每个共深度点所对应的时间,确定主要反射层在每个共深度点所对应的时间位置。判断过井剖面的反射层内是否出现反射同相轴;若是,则该反射层内存在次级反射层,以VSP井所在位置为起点,对反射同相轴进行追踪,确定反射同相轴在VSP时间域层速度模型上每个共深度点所对应的时间,进而确定次级反射层在每个共深度点所对应的时间位置;若否,则该反射层内不存在次级反射层。
根据VSP井位置处的地下地层的时间域层速度,确定对应的主要反射层与次要反射层的层速度,可以将反射层在VSP井位置处所对应的地层的时间域层速度作为该反射层的层速度;若反射层的厚度大于VSP的层厚度,则计算反射层对应的多个时间域层速度的平均速度,作为反射层的层速度,进而确定VSP时间域层速度模型。将VSP时间域层速度模型的时间域层速度转换为深度域层速度,获得VSP深度域层速度模型。
本方法通过过井剖面构建VSP深度域层速度模型,提高VSP速度模型的准确度,消除由常规VSP速度建模方法所带来的局限性。
应用示例
为便于理解本发明实施例的方案及其效果,以下给出一个具体应用示例。本领域技术人员应理解,该示例仅为了便于理解本发明,其任何具体细节并非意在以任何方式限制本发明。
图2a和图2b分别示出了根据本发明的一个实施例的零井源距VSP资料的X分量和Z分量的示意图。
图3示出了根据本发明的一个实施例的时间域纵波速度与横波速度的示意图。
根据如图2a、图2b所示的零井源距VSP资料的X分量和Z分量,计算VSP井位置处的地下地层的时间域层速度,包括纵波速度与横波速度,如图3所示。
图4示出了根据本发明的一个实施例的过井剖面的示意图。
根据VSP时间域层速度模型的空间范围,确定如图4所示的过井剖面对应的范围;在范围内的过井剖面上确定主要地震层位,以VSP井所在位置为起点,对每个主要地震层位的主同相轴进行追踪,确定主同相轴在非VSP时间域层速度模型上每个共深度点所对应的时间;根据每个共深度点所对应的时间,确定主要反射层在每个共深度点所对应的时间位置。判断过井剖面的反射层内是否出现反射同相轴;若是,则该反射层内存在次级反射层,以VSP井所在位置为起点,对反射同相轴进行追踪,确定反射同相轴在VSP时间域层速度模型上每个共深度点所对应的时间,进而确定次级反射层在每个共深度点所对应的时间位置;若否,则该反射层内不存在次级反射层。
图5示出了根据本发明的一个实施例的VSP深度域层速度模型的示意图。
根据VSP井位置处的地下地层的时间域层速度,确定对应的主要反射层与次要反射层的层速度,将反射层在VSP井位置处所对应的地层的时间域层速度作为该反射层的层速度;若反射层的厚度大于VSP的层厚度,则计算反射层对应的多个时间域层速度的平均速度,作为反射层的层速度,进而确定VSP时间域层速度模型。将VSP时间域层速度模型的时间域层速度转换为深度域层速度,获得VSP深度域层速度模型,如图5所示。
图6a和图6b分别示出了根据本发明的一个实施例的VSP深度域纵波、横波偏移剖面的示意图。
图7a和图7b分别示出了根据本发明的一个实施例的VSP时间域纵波、转换波偏移剖面的示意图。
利用图5所示的深度域层速度模型对非零偏VSP资料的X、Z分量的上行波场进行叠前逆时深度偏移,得到如图6a、图6b所示的VSP深度域纵、横波偏移剖面,将VSP深度域纵、横波偏移剖面转换成如图7a、图7b所示的VSP时间域纵波、转换波偏移剖面。
图8a、图8b分别示出了根据常规VSP速度建模得到的VSP时间域纵波、转换波偏移剖面的示意图。
对比图7a、7b、8a与8b发现:在非零偏VSP资料处理过程中,利用常规VSP速度建模方法构建的VSP速度模型与利用本发明提出的VSP速度建模方法相比,最后处理出的非零偏VSP剖面在处理质量上相差很大,利用常规VSP速度建模方法构建的VSP速度模型对非零偏VSP进行处理,不能得到反射层清晰的VSP处理剖面,无法在剖面上区分反射层的顶、底界面;而利用本发明提出的VSP速度建模方法得到的VSP速度模型对非零偏VSP资料进行处理,则可以得到反射层非常清晰的VSP处理剖面。
综上所述,本发明通过过井剖面构建VSP深度域层速度模型,提高VSP速度模型的准确度,消除由常规VSP速度建模方法所带来的局限性。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
根据本发明的实施例,提供了一种结合过井剖面的VSP深度域层速度模型构建系统,其特征在于,该系统包括:存储器,存储有计算机可执行指令;处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:根据零井源距VSP资料,计算VSP井位置处的地下地层的深度域层速度进而进行时深转换,获得VSP井位置处的地下地层的时间域层速度;根据过井剖面的地震层位,确定VSP时间域层速度模型的主要反射层;根据过井剖面的反射同相轴,确定VSP时间域层速度模型的每个主要反射层的次级反射层;根据VSP井位置处的地下地层的时间域层速度,确定对应的主要反射层与次要反射层的层速度,进而确定VSP时间域层速度模型;将VSP时间域层速度模型的时间域层速度转换为深度域层速度,获得VSP深度域层速度模型。
在一个示例中,根据过井剖面的地震层位,确定主要反射层包括:根据VSP时间域层速度模型的空间范围,确定过井剖面对应的范围;在范围内的过井剖面上确定主要地震层位,以VSP井所在位置为起点,对每个主要地震层位的主同相轴进行追踪,确定主同相轴上每个共深度点所对应的时间;根据每个共深度点所对应的时间,确定VSP时间域层速度模型的主要反射层在每个共深度点所对应的时间位置。
在一个示例中,根据过井剖面的反射同相轴,确定每个主要反射层的次级反射层包括:判断过井剖面的反射层内是否出现反射同相轴;若是,则该反射层内存在次级反射层,以VSP井所在位置为起点,对反射同相轴进行追踪,确定反射同相轴在VSP时间域层速度模型上每个共深度点所对应的时间,进而确定次级反射层在每个共深度点所对应的时间位置;若否,则该反射层内不存在次级反射层。
在一个示例中,根据VSP井位置处的地下地层的时间域层速度,确定对应的主要反射层与次要反射层的层速度为:将反射层在VSP井位置处所对应的地层的时间域层速度作为该反射层的层速度。
在一个示例中,根据VSP井位置处的地下地层的时间域层速度,确定对应的主要反射层与次要反射层的层速度为:若反射层的厚度大于VSP的层厚度,则计算反射层对应的多个时间域层速度的平均速度,作为反射层的层速度。
本系统通过过井剖面构建VSP深度域层速度模型,提高VSP速度模型的准确度,消除由常规VSP速度建模方法所带来的局限性。
本领域技术人员应理解,上面对本发明的实施例的描述的目的仅为了示例性地说明本发明的实施例的有益效果,并不意在将本发明的实施例限制于所给出的任何示例。
以上已经描述了本发明的各实施例,上述说明是示例性的,并非穷尽性的,并且也不限于所披露的各实施例。在不偏离所说明的各实施例的范围和精神的情况下,对于本技术领域的普通技术人员来说许多修改和变更都是显而易见的。
Claims (10)
1.一种结合过井剖面的VSP深度域层速度模型构建方法,其特征在于,包括:
根据零井源距VSP资料,计算VSP井位置处的地下地层的深度域层速度,进而进行时深转换,获得VSP井位置处的地下地层的时间域层速度;
根据过井剖面的地震层位,确定VSP时间域层速度模型的主要反射层;
根据所述过井剖面的反射同相轴,确定所述VSP时间域层速度模型的每个主要反射层的次级反射层;
根据所述VSP井位置处的地下地层的时间域层速度,确定对应的所述主要反射层与所述次要反射层的层速度,进而确定VSP时间域层速度模型;
将所述VSP时间域层速度模型的时间域层速度转换为深度域层速度,获得VSP深度域层速度模型。
2.根据权利要求1所述的结合过井剖面的VSP深度域层速度模型构建方法,其中,所述根据过井剖面的地震层位,确定主要反射层包括:
根据所述VSP时间域层速度模型的空间范围,确定所述过井剖面对应的范围;
在所述范围内的过井剖面上确定主要地震层位,以VSP井所在位置为起点,对每个主要地震层位的主同相轴进行追踪,确定所述主同相轴上每个共深度点所对应的时间;
根据所述每个共深度点所对应的时间,确定所述VSP时间域层速度模型的主要反射层在每个共深度点所对应的时间位置。
3.根据权利要求1所述的结合过井剖面的VSP深度域层速度模型构建方法,其中,所述根据所述过井剖面的反射同相轴,确定每个主要反射层的次级反射层包括:
判断所述过井剖面的反射层内是否出现反射同相轴;
若是,则该反射层内存在次级反射层,以所述VSP井所在位置为起点,对所述反射同相轴进行追踪,确定所述反射同相轴在所述VSP时间域层速度模型上每个共深度点所对应的时间,进而确定所述次级反射层在每个共深度点所对应的时间位置;
若否,则该反射层内不存在次级反射层。
4.根据权利要求1所述的结合过井剖面的VSP深度域层速度模型构建方法,其中,所述根据所述VSP井位置处的地下地层的时间域层速度,确定对应的所述主要反射层与所述次要反射层的层速度为:
将反射层在VSP井位置处所对应的地层的时间域层速度作为该反射层的层速度。
5.根据权利要求1所述的结合过井剖面的VSP深度域层速度模型构建方法,其中,所述根据所述VSP井位置处的地下地层的时间域层速度,确定对应的所述主要反射层与所述次要反射层的层速度为:
若反射层的厚度大于VSP的层厚度,则计算所述反射层对应的多个时间域层速度的平均速度,作为所述反射层的层速度。
6.一种结合过井剖面的VSP深度域层速度模型构建系统,其特征在于,该系统包括:
存储器,存储有计算机可执行指令;
处理器,所述处理器运行所述存储器中的计算机可执行指令,执行以下步骤:
根据零井源距VSP资料,计算VSP井位置处的地下地层的深度域层速度进而进行时深转换,获得VSP井位置处的地下地层的时间域层速度;
根据过井剖面的地震层位,确定VSP时间域层速度模型的主要反射层;
根据所述过井剖面的反射同相轴,确定所述VSP时间域层速度模型的每个主要反射层的次级反射层;
根据所述VSP井位置处的地下地层的时间域层速度,确定对应的所述主要反射层与所述次要反射层的层速度,进而确定VSP时间域层速度模型;
将所述VSP时间域层速度模型的时间域层速度转换为深度域层速度,获得VSP深度域层速度模型。
7.根据权利要求6所述的结合过井剖面的VSP深度域层速度模型构建系统,其中,所述根据过井剖面的地震层位,确定主要反射层包括:
根据所述VSP时间域层速度模型的空间范围,确定所述过井剖面对应的范围;
在所述范围内的过井剖面上确定主要地震层位,以VSP井所在位置为起点,对每个主要地震层位的主同相轴进行追踪,确定所述主同相轴上每个共深度点所对应的时间;
根据所述每个共深度点所对应的时间,确定所述VSP时间域层速度模型的主要反射层在每个共深度点所对应的时间位置。
8.根据权利要求6所述的结合过井剖面的VSP深度域层速度模型构建系统,其中,所述根据所述过井剖面的反射同相轴,确定每个主要反射层的次级反射层包括:
判断所述过井剖面的反射层内是否出现反射同相轴;
若是,则该反射层内存在次级反射层,以所述VSP井所在位置为起点,对所述反射同相轴进行追踪,确定所述反射同相轴在所述VSP时间域层速度模型上每个共深度点所对应的时间,进而确定所述次级反射层在每个共深度点所对应的时间位置;
若否,则该反射层内不存在次级反射层。
9.根据权利要求6所述的结合过井剖面的VSP深度域层速度模型构建系统,其中,所述根据所述VSP井位置处的地下地层的时间域层速度,确定对应的所述主要反射层与所述次要反射层的层速度为:
将反射层在VSP井位置处所对应的地层的时间域层速度作为该反射层的层速度。
10.根据权利要求6所述的结合过井剖面的VSP深度域层速度模型构建系统,其中,所述根据所述VSP井位置处的地下地层的时间域层速度,确定对应的所述主要反射层与所述次要反射层的层速度为:
若反射层的厚度大于VSP的层厚度,则计算所述反射层对应的多个时间域层速度的平均速度,作为所述反射层的层速度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910579729.1A CN112147699B (zh) | 2019-06-28 | 2019-06-28 | 结合过井剖面的vsp深度域层速度模型构建方法及系统 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910579729.1A CN112147699B (zh) | 2019-06-28 | 2019-06-28 | 结合过井剖面的vsp深度域层速度模型构建方法及系统 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112147699A true CN112147699A (zh) | 2020-12-29 |
CN112147699B CN112147699B (zh) | 2022-06-21 |
Family
ID=73891575
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910579729.1A Active CN112147699B (zh) | 2019-06-28 | 2019-06-28 | 结合过井剖面的vsp深度域层速度模型构建方法及系统 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112147699B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024030565A1 (en) * | 2022-08-05 | 2024-02-08 | Saudi Arabian Oil Company | Method to correct vsp data |
Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101071175A (zh) * | 2006-05-11 | 2007-11-14 | 中国石油集团东方地球物理勘探有限责任公司 | 零井源距垂直地震剖面纵横波数据深度域走廊叠加剖面处理方法 |
US20090213692A1 (en) * | 2008-02-22 | 2009-08-27 | Martinez Ruben D | Method for three dimensional seismic travel time tomography in transversely isotropic media |
US20100118653A1 (en) * | 2008-11-08 | 2010-05-13 | Ruiqing He | Vertical seismic profiling velocity estimation method |
CN102967882A (zh) * | 2012-11-16 | 2013-03-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 地层的层速度模型的建模方法 |
CN104360385A (zh) * | 2014-11-26 | 2015-02-18 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 结合地震地质认识综合建立初始深度层速度模型的方法 |
CN105259578A (zh) * | 2015-10-26 | 2016-01-20 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于垂直地震剖面资料的地层速度确定方法 |
CN105842736A (zh) * | 2016-05-27 | 2016-08-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种构建地层速度模型的方法 |
AU2014386285A1 (en) * | 2014-03-11 | 2016-08-11 | Landmark Graphics Corporation | Validation of depth-depth curves using time-seismic depth interval velocity |
CN109100784A (zh) * | 2018-06-08 | 2018-12-28 | 恒泰艾普(北京)能源科技研究院有限公司 | 三维vsp源检互换全波场成像方法 |
-
2019
- 2019-06-28 CN CN201910579729.1A patent/CN112147699B/zh active Active
Patent Citations (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101071175A (zh) * | 2006-05-11 | 2007-11-14 | 中国石油集团东方地球物理勘探有限责任公司 | 零井源距垂直地震剖面纵横波数据深度域走廊叠加剖面处理方法 |
US20090213692A1 (en) * | 2008-02-22 | 2009-08-27 | Martinez Ruben D | Method for three dimensional seismic travel time tomography in transversely isotropic media |
US20100118653A1 (en) * | 2008-11-08 | 2010-05-13 | Ruiqing He | Vertical seismic profiling velocity estimation method |
CN102967882A (zh) * | 2012-11-16 | 2013-03-13 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 地层的层速度模型的建模方法 |
AU2014386285A1 (en) * | 2014-03-11 | 2016-08-11 | Landmark Graphics Corporation | Validation of depth-depth curves using time-seismic depth interval velocity |
CN104360385A (zh) * | 2014-11-26 | 2015-02-18 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 结合地震地质认识综合建立初始深度层速度模型的方法 |
CN105259578A (zh) * | 2015-10-26 | 2016-01-20 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 基于垂直地震剖面资料的地层速度确定方法 |
CN105842736A (zh) * | 2016-05-27 | 2016-08-10 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 一种构建地层速度模型的方法 |
CN109100784A (zh) * | 2018-06-08 | 2018-12-28 | 恒泰艾普(北京)能源科技研究院有限公司 | 三维vsp源检互换全波场成像方法 |
Non-Patent Citations (4)
Title |
---|
L. J. HOFMAN,等: "A Shallow Seismic Velocity Model for the Groningen Area in the Netherlands", 《JOURNAL OF GEOPHYSICAL RESEARCH: SOLID EARTH》 * |
刘兵,等: "琼东南盆地坡折带斜井时深关系研究", 《天然气工业》 * |
张山,等: "复杂山前带有偏VSP数据处理方法研究", 《石油物探》 * |
李文杰,等: "利用VSP资料反演地层层速度的一种新途径", 《石油物探》 * |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024030565A1 (en) * | 2022-08-05 | 2024-02-08 | Saudi Arabian Oil Company | Method to correct vsp data |
Also Published As
Publication number | Publication date |
---|---|
CN112147699B (zh) | 2022-06-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US8576663B2 (en) | Multicomponent seismic inversion of VSP data | |
Nakanishi et al. | Detailed structural image around splay‐fault branching in the Nankai subduction seismogenic zone: Results from a high‐density ocean bottom seismic survey | |
CN107817526B (zh) | 叠前地震道集分段式振幅能量补偿方法及系统 | |
CN101738636B (zh) | 一种三维vsp高斯束法多波联合偏移成像方法 | |
CN105093292A (zh) | 一种地震成像的数据处理方法和装置 | |
CN105259581B (zh) | 一种地震资料时深转换方法 | |
Carpenter et al. | Estimating site response with recordings from deep boreholes and HVSR: Examples from the Mississippi embayment of the central United States | |
Poulin et al. | Focal-time analysis: A new method for stratigraphic depth control of microseismicity and induced seismic events | |
AU2004220765A1 (en) | Method for stable estimation of anisotropic parameters for P-wave prestack imaging | |
Pech et al. | Quartic moveout coefficient: 3D description and application to tilted TI media | |
CN104422959A (zh) | 一种检测储层边界的曲率属性方法 | |
Skuce et al. | Passive-roof duplexes under the Rocky Mountain foreland basin, Alberta | |
CN112147699B (zh) | 结合过井剖面的vsp深度域层速度模型构建方法及系统 | |
Ren et al. | Prestack reverse time imaging in tunnels based on the decoupled nonconversion elastic-wave equation | |
Roslov et al. | Deep seismic investigations in the Barents and Kara Seas | |
Yan et al. | Seismic imaging and velocity analysis for an Alberta Foothills seismic survey | |
CN111175825A (zh) | 深度域速度建模方法 | |
Harvey et al. | Structural variations of the crust in the Southwestern Cape, deduced from seismic receiver functions | |
Sorlien et al. | Late Miocene–Quaternary fault evolution and interaction in the southern California Inner Continental Borderland | |
US4847813A (en) | Method for extending the lateral subsurface coverage in VSP surveys | |
Maione | Discovery of ring faults associated with salt withdrawal basins, Early Cretaceous age, in the East Texas Basin | |
Rohilla et al. | Shear‐wave velocity structure of the Koyna–Warna Region, western India, through modeling of P‐receiver functions | |
Li et al. | SV-P extraction and imaging for far-offset vertical seismic profile data | |
Stephenson et al. | High-resolution seismic reflection surveys and modeling across an area of high damage from the 1994 Northridge earthquake, Sherman Oaks, California | |
Murty et al. | Mapping the sub-trappean Mesozoic sediments in the western part of Narmada–Tapti region of Deccan Volcanic Province, India |
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 |