CN111624662B - 速度场校正方法、装置和设备 - Google Patents

速度场校正方法、装置和设备 Download PDF

Info

Publication number
CN111624662B
CN111624662B CN201910148879.7A CN201910148879A CN111624662B CN 111624662 B CN111624662 B CN 111624662B CN 201910148879 A CN201910148879 A CN 201910148879A CN 111624662 B CN111624662 B CN 111624662B
Authority
CN
China
Prior art keywords
velocity
target
model
layer
depth migration
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
CN201910148879.7A
Other languages
English (en)
Other versions
CN111624662A (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.)
China National Petroleum Corp
BGP Inc
Original Assignee
China National Petroleum Corp
BGP Inc
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 China National Petroleum Corp, BGP Inc filed Critical China National Petroleum Corp
Priority to CN201910148879.7A priority Critical patent/CN111624662B/zh
Publication of CN111624662A publication Critical patent/CN111624662A/zh
Application granted granted Critical
Publication of CN111624662B publication Critical patent/CN111624662B/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/40Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging
    • G01V1/44Seismology; Seismic or acoustic prospecting or detecting specially adapted for well-logging using generators and receivers in the same well
    • G01V1/48Processing data
    • G01V1/50Analysing data
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V2210/00Details of seismic processing or analysis
    • G01V2210/60Analysis
    • G01V2210/61Analysis by combining or comparing a seismic data set with other data
    • G01V2210/616Data from specific type of measurement
    • G01V2210/6169Data from specific type of measurement using well-logging

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Acoustics & Sound (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geology (AREA)
  • Remote Sensing (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本发明提供了一种速度场校正方法、装置和设备,其中,该方法包括:获取工区中目标井处的地震资料,其中,所述地震资料包括:初始深度偏移速度;利用随钻地震技术,获取目标井处的目标层位的VSP层速度;计算所述VSP层速度与所述初始深度偏移速度的比例系数;建立目标井处地质构造模型;确定比例系数三维模型;根据所述初始深度偏移速度和所述比例系数三维模型,确定校正后的深度偏移层速度场。通过上述方案可以实时确定钻井实测的VSP层速度与钻前的初始深度偏移速度之间的差异,对初始深度偏移速度进行校正,从而快速准确地调整钻井轨迹,提高储层钻遇率;通过采用随钻地震技术,在校正过程中无需暂停钻井工作,降低了勘探开发的风险。

Description

速度场校正方法、装置和设备
技术领域
本发明涉及石油及天然气开发技术领域,特别涉及一种速度场校正方法、装置和设备。
背景技术
在油井设计和钻井过程中,一般先建立钻前地下速度模型,并使用该模型对地震数据进行偏移成像;随后解释人员对地震成像数据体进行解释,进行储层识别以及井位确定等工作。如果钻前地下速度模型不够精确,钻探目标位置就会有偏移,从而导致钻井的准确度不高。
实际钻井过程中,由于获得的钻井实测数据与钻前地下速度模型存在一定差异,现有技术中,通常是将实际钻井过程中获得的实际速度数据对偏移成像进行一维的纵向拉伸,以实现深度预测的更新。在钻前地下速度模型不够精确、钻探目标横向位置有偏移的情况下,一维的纵向拉伸并不能完全将钻探目标归位到正确的空间位置。在上述情况下,需要用侧钻或酸化压裂等技术进行弥补,以找到目标储层的空间位置。然后,通过侧钻或酸化压裂等技术进行弥补的方式来更新钻前地下速度模型,需要暂停钻井工作,从而影响钻井的周期和效率,并增大了钻探风险和钻井成本。
针对上述问题,目前尚未提出有效的解决方案。
发明内容
本申请实施例提供了一种速度场校正方法、装置和设备,以解决现有技术中无法有效校正钻前地下速度模型问题。
本申请实施例提供了一种速度场校正方法,包括:获取工区中目标井处的地震资料,其中,所述地震资料包括:初始深度偏移速度;利用随钻地震技术,获取目标井处的目标层位的VSP层速度;计算所述VSP层速度与所述初始深度偏移速度的比例系数;根据所述地震资料建立目标井处地质构造模型;根据所述比例系数和所述地质构造模型,确定比例系数三维模型;根据所述初始深度偏移速度和所述比例系数三维模型,确定校正后的深度偏移层速度场。
在一个实施例中,根据所述比例系数和所述地质构造模型,确定比例系数三维模型,包括:在目标层位处为单个目标井时,将所述目标井处的比例系数沿所述地质构造模型外推,得到所述比例系数三维模型;在目标层位为多个目标井时,通过反距离加权插值法确定所述多个目标井之间的比例系数,将所述多个目标井处的比例系数以及所述多个目标井之间的比例系数沿所述地质构造模型外推,得到所述比例系数三维模型。
在一个实施例中,计算所述VSP层速度与所述初始深度偏移速度的比例系数,包括:将所述VSP层速度与所述初始深度偏移速度做相除运算;将相除运算得到的值,作为所述比例系数。
在一个实施例中,根据所述地震资料建立目标井处地质构造模型,包括:获取所述目标井处的初始深度偏移层速度;根据所述初始深度偏移层速度和所述初始深度偏移速度,确定速度控制层层位数据;对所述速度控制层层位数据进行沿层内插,得到目标井处地质构造模型。
在一个实施例中,在得到校正后的深度偏移层速度场后,还包括:获取所述目标井的各向异性参数;根据所述各向异性参数和所述地质构造模型,建立各向异性三维数据模型;根据所述各向异性三维数据模型和校正后的深度偏移层速度场,建立各向异性深度偏移层速度模型。
在一个实施例中,建立各向异性三维数据模型包括:在目标层位处为单个目标井时,将所述目标井处各向异性参数沿所述地质构造模型外推,得到所述各向异性三维数据模型;在目标层位处为多个目标井时,通过反距离加权插值法确定所述多个目标井之间的各向异性参数,将所述多个目标井处各向异性参数以及所述多个目标井之间的各向异性参数沿所述地质构造模型外推,得到所述各向异性三维数据模型。
在一个实施例中,获取目标井处的目标层位的VSP层速度,包括:确定所述初始深度偏移速度模型在深度上的起算面;根据所述初始深度偏移速度模型在深度上的起算面,作为起算面确定所述VSP层速度。
本申请实施例还提供了一种速度场校正装置,包括:第一获取模块,用于获取工区中目标井处的地震资料,其中,所述地震资料包括:初始深度偏移速度;第二获取模块,用于利用随钻地震技术,获取目标井处的目标层位的VSP层速度;计算模块,用于计算所述VSP层速度与所述初始深度偏移速度的比例系数;建立模块,用于根据所述地震资料,建立目标井处地质构造模型;第一确定模块,用于根据所述比例系数和所述地质构造模型,确定比例系数三维模型;第二确定模块,用于根据所述初始深度偏移速度和所述比例系数三维模型,确定校正后的深度偏移层速度场。
在一个实施例中,所述第一确定模块包括:第一外推单元,用于在目标层位处为单个目标井时,将所述目标井处的比例系数沿所述地质构造模型外推,得到所述比例系数三维模型;第一处理单元,用于在目标层位为多个目标井时,通过反距离加权插值法确定所述多个目标井之间的比例系数,将所述多个目标井处的比例系数以及所述多个目标井之间的比例系数沿所述地质构造模型外推,得到所述比例系数三维模型。
在一个实施例中,所述计算模块包括:运算单元,用于将所述VSP层速度与所述初始深度偏移速度做相除运算;第二处理单元,用于将相除运算得到的值,作为所述比例系数。
在一个实施例中,所述建立模块包括:第一获取单元,用于获取所述目标井处的初始深度偏移层速度;确定单元,用于根据所述初始深度偏移层速度和所述初始深度偏移速度,确定速度控制层层位数据;内插单元,用于对所述速度控制层层位数据进行沿层内插,得到目标井处地质构造模型。
在一个实施例中,还包括:第二获取单元,用于获取所述目标井的各向异性参数;
第一建立单元,用于根据所述各向异性参数和所述地质构造模型,建立各向异性三维数据模型;第二建立单元,用于根据所述各向异性三维数据模型和校正后的深度偏移层速度场,建立各向异性深度偏移层速度模型。
在一个实施例中,所述第二建立单元包括:第二外推单元,用于在目标层位处为单个目标井时,将所述目标井处各向异性参数沿所述地质构造模型外推,得到所述各向异性三维数据模型;第三处理单元,用于在目标层位处为多个目标井时,通过反距离加权插值法确定所述多个目标井之间的各向异性参数,将所述多个目标井处各向异性参数以及所述多个目标井之间的各向异性参数沿所述地质构造模型外推,得到所述各向异性三维数据模型。
本申请实施例还提供了一种速度场校正设备,包括处理器以及用于存储处理器可执行指令的存储器,所述处理器执行所述指令时实现所述速度场校正方法的步骤。
本申请实施例还提供了一种计算机可读存储介质,其上存储有计算机指令,所述指令被执行时所述速度场校正方法的步骤。
在本申请实施例中,提供了一种速度场校正方法,通过利用随钻地震技术,获取目标井处的目标层位的VSP层速度,计算VSP层速度与初始深度偏移速度的比例系数,可以实时确定钻井实测的VSP层速度与钻前的初始深度偏移速度之间的差异;通过初始深度偏移速度和比例系数三维模型,确定校正后的深度偏移层速度场,使得在钻前地下速度模型不够精确的情况下,可以对初始深度偏移速度进行校正,从而可以快速、准确地调整钻井轨迹,提高储层钻遇率;进一步的,由于采用了随钻地震技术,使得在速度场校正过程中无需暂停钻井工作,降低了勘探开发的风险。
附图说明
此处所说明的附图用来提供对本申请的进一步理解,构成本申请的一部分,并不构成对本申请的限定。在附图中:
图1是根据本申请实施例提供的速度场校正方法步骤示意图;
图2是根据本申请具体实施例提供的速度场校正方法示意图;
图3是根据本申请具体实施例提供的原始偏移速度与VSP层速度对比示意图;
图4是根据本申请具体实施例提供的井点VSP层速度与原始偏移速度的比例系数曲线示意图;
图5是根据本申请具体实施例提供的工区地质框架模型剖面示意图;
图6是根据本申请具体实施例提供的原始偏移速度剖面与VSP校正后速度剖面对比示意图;
图7是根据本申请具体实施例提供的重新处理前后深度偏移剖面对比示意图;
图8是根据本申请实施例提供的速度场校正装置示意图;
图9是根据本申请实施例提供的速度场校正电子设备示意图。
具体实施方式
下面将参考若干示例性实施方式来描述本申请的原理和精神。应当理解,给出这些实施方式仅仅是为了使本领域技术人员能够更好地理解进而实现本申请,而并非以任何方式限制本申请的范围。相反,提供这些实施方式是为了使本申请公开更加透彻和完整,并且能够将本公开的范围完整地传达给本领域的技术人员。
本领域的技术人员知道,本申请的实施方式可以实现为一种系统、装置设备、方法或计算机程序产品。因此,本申请公开可以具体实现为以下形式,即:完全的硬件、完全的软件(包括固件、驻留软件、微代码等),或者硬件和软件结合的形式。
考虑到现有的更新钻前地下速度模型的技术方案中,在钻前地下速度模型不够精确、钻探目标横向位置有偏移的情况下,将实际钻井过程中获得的实际速度数据对偏移成像进行一维的纵向拉伸,无法准确将目标归位到正确的空间位置;通过侧钻或酸化压裂等技术进行弥补的方式来更新钻前地下速度模型,需要暂停钻井工作,从而影响钻井的周期和效率,并增大了钻探风险和钻井成本。
基于上述无法有效更新钻前地下速度模型的问题,本发明实施例提供了一种速度场校正方法,如图1所示,可以包括以下步骤:
步骤S101:获取工区中目标井处的地震资料,其中,地震资料包括:初始深度偏移速度。
由于基于VSP(Vertical Seismic Profiling,垂直地震剖面)随钻的速度场校正是在原始深度偏移处理的基础上开展的,因此,在进行速度场校正之前,可以先获取工区中有VSP测井资料的目标井处的地震资料,上述地震资料可以包括但不限于:初始深度偏移速度、初始深度偏移层速度、初始深度偏移资料。
其中,VSP是一种地震观测方法,VSP具有地震波单程衰减、地震信号频率较高、速度分析精度高、可以估算各向异性参数等优点。
进一步的,钻前地下速度模型可以为前期深度偏移处理中建立的初始深度偏移层速度模型,可以将叠前时间偏移中的均方根速度通过DIX公式(Dix formula,迪克斯公式)转换得到初始深度偏移层速度,并通过沿层剩余分析及网格层析方法优化得到钻前地下速度模型。
步骤S102:利用随钻地震技术,获取目标井处的目标层位的VSP层速度。
考虑到钻井实测数据与钻前地下速度模型存在一定差异,因此,可以在钻井过程中,利用钻头钻进过程中钻头与钻遇地层之间的撞击、摩擦所产生的微弱地震信号作为信号源(震源),实时获取目标井处的目标层位的VSP层速度,从而能够快速地确定钻井实测数据与钻前地下速度模型之间的差异,并对钻前地下速度模型进行校正。其中,VSP层速度在深度上的起算面与上述初始深度偏移速度的深度起算面可以是一致的。
步骤S103:计算VSP层速度与初始深度偏移速度的比例系数。
在VSP层速度在深度上的起算面与上述初始深度偏移速度的深度起算面一致的情况下,计算目标井处各个层位VSP层速度与初始深度偏移速度的比例系数,得到比例系数曲线。其中,可以将目标井处各个层位的VSP层速度与初始深度偏移速度作相除运算,并将作相除运算的到的值,作为目标井处各个层位的比例系数。
步骤S104:根据地震资料建立目标井处地质构造模型。
考虑到需要建立三维模型,可以先获取目标井处的初始深度偏移层速度,并利用初始深度偏移层速度及初始深度偏移资料解释速度控制层,得到速度控制层层位数据,并利用速度控制层层位沿层内插方式,建立目标井处地质构造模型。
由于平面层速度场所反映的是一个粗糙的速度场,没有考虑到地层重叠、褶皱或者缺失、剥蚀等变化。因此,可以将解释得到的速度控制层层位数据通过沿层内插的方式,结合地质要素,建立地质构造模型。其中,在已观测点的区域内估算未观测点的数据的过程称为内插;在已观测点的区域外估算未观测点的数据的过程称为外推。
步骤S105:根据比例系数和地质构造模型,确定比例系数三维模型。
在工区中目标层位只有单个目标井的VSP测井资料的情况下,可以将目标井的比例系数沿上述地质构造模型外推,得到比例系数三维模型。其中,上述外推法是根据一组观测值,计算观测范围以外同一对象近似值的方法,通过外推法可以根据已知部分的地质规律来预测或推断未知部分的情况。例如,在储量计算及圈定矿体时,可根据已有的勘探工程资料,结合地质构造及矿体的变化趋势,推断未知部分的矿体可能分布界线。
在工区中目标层位有多个井的VSP测井资料的情况下,可以将多个井的比例系数沿上述地质构造模型外推,得到比例系数三维模型。其中,上述多个井之间的比例系数的值可以通过反距离加权插值法确定,将多个目标井处的比例系数结合多个目标井之间的比例系数,沿上述地质构造模型外推。其中,反距离加权插值法,是基于相近相似的原理:即两个物体离得近,它们的性质就越相似,反之,离得越远则相似性越小。以插值点与样本点间的距离为权重进行加权平均,离插值点越近的样本点赋予的权重越大。
步骤S106:根据初始深度偏移速度和比例系数三维模型,确定校正后的深度偏移层速度场。
可以将上述比例系数三维模型与初始深度偏移速度相乘,得到校正后的深度偏移层速度场,在将上述比例系数三维模型与初始深度偏移速度相乘时,将对应深度、目标层位的数据相乘。
地震偏移处理又称为偏移归位处理、地震成像处理,使反射波和绕射波空间正确归位;使绕射波、弯曲界面上的回转波、陡倾界面或断层界面上的反射波收敛到正确的空间位置上,从而得到地下真实地质构造的过程。在地震资料叠前偏移处理中,如果忽视各向异性,可能导致地质体垂向深度与横向位置的偏差,也会引起陡倾地层信息的丢失。因此,在地震资料处理中,特别是在速度分析和偏移成像两个步骤中需要考虑速度的各向异性。
可以在确定校正后的深度偏移速度场后,获取目标井上分层数据与地震分层数据,并根据上述目标井上分层数据与地震分层数据,确定目标井上各向异性参数,得到各向异性参数曲线。并利用正后的深度偏移层速度场以及各向异性参数三维模型,开展叠前深度偏移处理。
其中,在工区中目标层位只有单个目标井的VSP测井资料的情况下,可以将目标井的各向异性参数沿上述地质构造模型外推,得到各向异性参数三维模型。其中,上述外推法是根据一组观测值,计算观测范围以外同一对象近似值的方法。地质工作中常用这种方法,即根据已知部分的地质规律来预测或推断未知部分的情况。例如,在储量计算及圈定矿体时,可根据已有的勘探工程资料,结合地质构造及矿体的变化趋势,推断未知部分的矿体可能分布界线。
进一步的,在工区中目标层位有多个井的VSP测井资料的情况下,可以将多个井的各向异性参数沿上述地质构造模型外推,得到各向异性参数三维模型。其中,上述多个井之间的各向异性参数的值可以通过反距离加权插值法确定,将多个目标井处的各向异性参数结合多个目标井之间的各向异性参数,沿上述地质构造模型外推。反距离加权插值法,是基于相近相似的原理,即两个物体离得近,它们的性质就越相似,反之,离得越远则相似性越小。以插值点与样本点间的距离为权重进行加权平均,离插值点越近的样本点赋予的权重越大。
从以上的描述中,可以看出,本申请实施例实现了如下技术效果:通过利用随钻地震技术,获取目标井处的目标层位的实时VSP层速度,计算VSP层速度与初始深度偏移速度的比例系数,可以实时确定钻井实测的VSP层速度与钻前的初始深度偏移速度之间的差异;通过初始深度偏移速度和比例系数三维模型,确定校正后的深度偏移层速度场,使得在钻前地下速度模型不够精确的情况下,可以快速、准确地校正初始深度偏移速度场,从而提高深度偏移成像的精度、提高储层钻遇率;并且由于采用随钻地震技术,在速度场校正过程中无需暂停钻井工作,降低了勘探开发的风险。
下面结合一个具体实施例对上述方法进行说明,然而,值得注意的是,该具体实施例仅是为了更好地说明本申请,并不构成对本申请的不当限定。
本发明实施提供了一种速度场校正方法,以塔里木盆地FY202-1X井区实际地震工区为例,如图2所示,可以包括:
步骤S201:VSP速度资料的整理与分析。
首先建立塔里木盆地FY202-1X井区地震工区,加载工区中目标井的原始偏移速度、井上VSP层速度,计算井上VSP速度与原始偏移速度比例系数。
具体的,可以通过VSP随钻地震导向技术,获取工区中目标井的VSP层速度,并获取工区进行VSP随钻测量的原始偏移速度。进一步的,为了使得计算比例系数时可以逐层计算分析,需要确保VSP层速度曲线在深度上的起算面与加载进的原始偏移速度深度起算面一致,其中,原始偏移速度与VSP层速度对比图可以如图3所示。
提取有VSP层速度资料的目标井处的原始偏移速度资料,计算井上VSP速度与原始偏移速度的比例系数。并按照以下公式计算各目标井处目标层位的VSP层速度与原始偏移速度的比例系数:
比例系数=VSP层速度/原始偏移速度
得到如图4所示的井点VSP层速度与原始偏移速度的比例系数曲线,其中,横坐标为井点VSP层速度与原始偏移速度的比值,纵坐标为深度值。
步骤S202:利用工区速度控制层,建立工区地质框架模型。
利用原始偏移层速度资料以及原始深度偏移资料解释速度控制层,得到速度控制层层位数据,为建立工区地质框架模型做准备;通过使用速度控制层层位沿层内插方式,建立工区地质框架模型。工区地质框架模型剖面如图5所示,其中,横坐标第一行为地震道号、第二行为地震线号,表示平面位置,纵坐标为深度值。
步骤S203:生成比例系数的空间模型,得到比例系数三维体。
由于每口井VSP测量深度的不同,在不同的层位可以采用不同建模策略:在目标层位只有单井VSP测量资料时,将井点处比例系数沿地质构造模型外推,得到比例系数三维模型;在目标层位有多个井的VSP测量资料是,将比例系数沿地质构造模型外推,井间不同位置比例系数值可以通过反距离加权插值法综合周围各个井的VSP测量资料确定,得到多井及地质构造模型共同控制的比例系数三维模型。
其中,反距离加权插值法是基于相近相似的原理,即两个物体离得近,它们的性质就越相似,反之,离得越远则相似性越小。以插值点与样本点间的距离为权重进行加权平均,离插值点越近的样本点赋予的权重越大。
步骤S204:由原始偏移速度与比例系数三维模型计算得到VSP校正后速度体。
将原始偏移速度体和比例系数体逐层相乘,计算VSP校正后速度体。原始偏移速度剖面与VSP校正后速度剖面对比可以如图6所示,其中,横坐标第一行为地震道号、第二行为地震线号,表示平面位置,纵坐标为速度值。从图6中不难发现,校正后的速度剖面可以更加精确的体现地下地层信息。进一步的,还可以利用校正后的速度体及结合地质构造模型建立的各向异性参数体,开展叠前深度偏移处理。
利用VSP校正后的速度体,开展FY202-1X井区地震资料叠前深度偏移处理,VSP速度校正前后深度偏移剖面对比,如图7所示,可以发现,利用VSP速度校正后的深度偏移数据体对地下缝洞体的成像更为清晰,可以更好的指导后期钻井轨迹的调整,并有效提高了储层钻遇率。
基于同一发明构思,本申请实施例中还提供了一种速度场校正装置,如下面的实施例所述。由于速度场校正装置解决问题的原理与速度场校正方法相似,因此速度场校正装置的实施可以参见速度场校正方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。图8是本申请实施例的速度场校正装置的一种结构框图,如图8所示,包括:第一获取模块801、第二获取模块802、计算模块803、建立模块804、第一确定模块805以及第二确定模块806,下面对该结构进行说明。
第一获取模块801,用于获取工区中目标井处的地震资料,其中,地震资料包括:初始深度偏移速度。
第二获取模块802,用于利用随钻地震技术,获取目标井处的目标层位的VSP层速度。
考虑到钻井实测数据与钻前地下速度模型存在一定差异,因此,可以在钻井过程中,利用钻头钻进过程中钻头与钻遇地层之间的撞击、摩擦所产生的微弱地震信号作为信号源(震源),实时获取目标井处的目标层位的VSP层速度,从而能够快速地确定钻井实测数据与钻前地下速度模型之间的差异,并对钻前地下速度模型进行校正。其中,VSP层速度在深度上的起算面与上述初始深度偏移速度的深度起算面可以是一致的。
计算模块803,用于计算VSP层速度与初始深度偏移速度的比例系数。
在一个实施例中,计算模块可以包括:运算单元,用于将VSP层速度与初始深度偏移速度做相除运算;第二处理单元,用于将相除运算得到的值,作为比例系数。
建立模块804,用于根据地震资料,建立目标井处地质构造模型。
在一个实施例中,建立模块可以包括:第一获取单元,用于获取目标井处的初始深度偏移层速度;确定单元,用于根据初始深度偏移层速度和初始深度偏移速度,确定速度控制层层位数据;以及内插单元,用于对速度控制层层位数据进行沿层内插,得到目标井处地质构造模型。
第一确定模块805,用于根据比例系数和地质构造模型,确定比例系数三维模型。
在一个实施例中,第一确定模块可以包括:第一外推单元,可以用于在目标层位处为单个目标井时,将所述目标井处的比例系数沿所述地质构造模型外推,得到所述比例系数三维模型;以及第一处理单元,可以用于在目标层位为多个目标井时,通过反距离加权插值法确定所述多个目标井之间的比例系数,将所述多个目标井处的比例系数以及所述多个目标井之间的比例系数沿所述地质构造模型外推,得到所述比例系数三维模型。
第二确定模块806,用于根据初始深度偏移速度和比例系数三维模型,确定校正后的深度偏移层速度场。
在一个实施例中,上述速度场校正装置还可以包括:第二获取单元,用于获取目标井的各向异性参数;第一建立单元,用于根据各向异性参数和地质构造模型,建立各向异性三维数据模型;第二建立单元,用于根据各向异性三维数据模型和校正后的深度偏移层速度场,建立各向异性深度偏移层速度模型。
在一个实施例中,上述第二建立单元可以包括:第二外推单元,可以用于在目标层位处为单个目标井时,将所述目标井处各向异性参数沿所述地质构造模型外推,得到所述各向异性三维数据模型;第三处理单元,可以用于在目标层位处为多个目标井时,通过反距离加权插值法确定所述多个目标井之间的各向异性参数,将所述多个目标井处各向异性参数以及所述多个目标井之间的各向异性参数沿所述地质构造模型外推,得到所述各向异性三维数据模型。
本申请实施方式还提供了一种电子设备,具体可以参阅图9所示的基于本申请实施例提供的速度场校正方法的电子设备组成结构示意图,电子设备具体可以包括输入设备91、处理器92、存储器93。其中,输入设备91具体可以用于获取工区中目标井处的地震资料、目标井处的目标层位的VSP层速度,其中,地震资料包括:初始深度偏移速度。处理器92具体可以用于计算VSP层速度与初始深度偏移速度的比例系数;根据地震资料建立目标井处地质构造模型;根据比例系数和地质构造模型,确定比例系数三维模型;根据初始深度偏移速度和比例系数三维模型,确定校正后的深度偏移层速度场。存储器93具体可以用于存储工区中目标井处的地震资料、目标井处的目标层位的VSP层速度、校正后的深度偏移层速度场等参数。
在本实施方式中,所述输入设备具体可以是用户和计算机系统之间进行信息交换的主要装置之一。所述输入设备可以包括键盘、鼠标、摄像头、扫描仪、光笔、手写输入板、语音输入装置等;输入设备用于把原始数据和处理这些数的程序输入到计算机中。所述输入设备还可以获取接收其他模块、单元、设备传输过来的数据。所述处理器可以按任何适当的方式实现。例如,处理器可以采取例如微处理器或处理器以及存储可由该(微)处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌入微控制器的形式等等。所述存储器具体可以是现代信息技术中用于保存信息的记忆设备。所述存储器可以包括多个层次,在数字系统中,只要能保存二进制数据的都可以是存储器;在集成电路中,一个没有实物形式的具有存储功能的电路也叫存储器,如RAM、FIFO等;在系统中,具有实物形式的存储设备也叫存储器,如内存条、TF卡等。
在本实施方式中,该电子设备具体实现的功能和效果,可以与其它实施方式对照解释,在此不再赘述。
本申请实施方式中还提供了一种基于速度场校正方法的计算机存储介质,所述计算机存储介质存储有计算机程序指令,在所述计算机程序指令被执行时实现:根据初始深度偏移速度和比例系数三维模型,确定校正后的深度偏移层速度场。
在本实施方式中,上述存储介质包括但不限于随机存取存储器(Random AccessMemory,RAM)、只读存储器(Read-Only Memory,ROM)、缓存(Cache)、硬盘(Hard DiskDrive,HDD)或者存储卡(Memory Card)。所述存储器可以用于存储计算机程序指令。网络通信单元可以是依照通信协议规定的标准设置的,用于进行网络连接通信的接口。
在本实施方式中,该计算机存储介质存储的程序指令具体实现的功能和效果,可以与其它实施方式对照解释,在此不再赘述。
显然,本领域的技术人员应该明白,上述的本申请实施例的各模块或各步骤可以用通用的计算装置来实现,它们可以集中在单个的计算装置上,或者分布在多个计算装置所组成的网络上,可选地,它们可以用计算装置可执行的程序代码来实现,从而,可以将它们存储在存储装置中由计算装置来执行,并且在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤,或者将它们分别制作成各个集成电路模块,或者将它们中的多个模块或步骤制作成单个集成电路模块来实现。这样,本申请实施例不限制于任何特定的硬件和软件结合。
应该理解,以上描述是为了进行图示说明而不是为了进行限制。通过阅读上述描述,在所提供的示例之外的许多实施方式和许多应用对本领域技术人员来说都将是显而易见的。因此,本申请的范围不应该参照上述描述来确定,而是应该参照前述权利要求以及这些权利要求所拥有的等价物的全部范围来确定。
以上所述仅为本申请的优选实施例而已,并不用于限制本申请,对于本领域的技术人员来说,本申请实施例可以有各种更改和变化。凡在本申请的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本申请的保护范围之内。

Claims (15)

1.一种速度场校正方法,其特征在于,包括:
获取工区中目标井处的地震资料,其中,所述地震资料包括:初始深度偏移速度;
利用随钻地震技术,获取目标井处的目标层位的VSP层速度;
计算所述VSP层速度与所述初始深度偏移速度的比例系数;
根据所述地震资料建立目标井处地质构造模型;
根据所述比例系数和所述地质构造模型,确定比例系数三维模型;
根据所述初始深度偏移速度和所述比例系数三维模型,确定校正后的深度偏移层速度场。
2.根据权利要求1所述的方法,其特征在于,根据所述比例系数和所述地质构造模型,确定比例系数三维模型,包括:
在目标层位处为单个目标井时,将所述目标井处的比例系数沿所述地质构造模型外推,得到所述比例系数三维模型;
在目标层位为多个目标井时,通过反距离加权插值法确定所述多个目标井之间的比例系数,将所述多个目标井处的比例系数以及所述多个目标井之间的比例系数沿所述地质构造模型外推,得到所述比例系数三维模型。
3.根据权利要求1所述的方法,其特征在于,计算所述VSP层速度与所述初始深度偏移速度的比例系数,包括:
将所述VSP层速度与所述初始深度偏移速度做相除运算;
将相除运算得到的值,作为所述比例系数。
4.根据权利要求1所述的方法,其特征在于,根据所述地震资料建立目标井处地质构造模型,包括:
获取所述目标井处的初始深度偏移层速度;
根据所述初始深度偏移层速度和所述初始深度偏移速度,确定速度控制层层位数据;
对所述速度控制层层位数据进行沿层内插,得到目标井处地质构造模型。
5.根据权利要求1所述的方法,其特征在于,在得到校正后的深度偏移层速度场后,还包括:
获取所述目标井的各向异性参数;
根据所述各向异性参数和所述地质构造模型,建立各向异性三维数据模型;
根据所述各向异性三维数据模型和校正后的深度偏移层速度场,建立各向异性深度偏移层速度模型。
6.根据权利要求5所述的方法,其特征在于,建立各向异性三维数据模型包括:
在目标层位处为单个目标井时,将所述目标井处各向异性参数沿所述地质构造模型外推,得到所述各向异性三维数据模型;
在目标层位处为多个目标井时,通过反距离加权插值法确定所述多个目标井之间的各向异性参数,将所述多个目标井处各向异性参数以及所述多个目标井之间的各向异性参数沿所述地质构造模型外推,得到所述各向异性三维数据模型。
7.根据权利要求1所述的方法,其特征在于,获取目标井处的目标层位的VSP层速度,包括:
确定所述初始深度偏移速度模型在深度上的起算面;
根据所述初始深度偏移速度模型在深度上的起算面,作为起算面确定所述VSP层速度。
8.一种速度场校正装置,其特征在于,包括:
第一获取模块,用于获取工区中目标井处的地震资料,其中,所述地震资料包括:初始深度偏移速度;
第二获取模块,用于利用随钻地震技术,获取目标井处的目标层位的VSP层速度;
计算模块,用于计算所述VSP层速度与所述初始深度偏移速度的比例系数;
建立模块,用于根据所述地震资料,建立目标井处地质构造模型;
第一确定模块,用于根据所述比例系数和所述地质构造模型,确定比例系数三维模型;
第二确定模块,用于根据所述初始深度偏移速度和所述比例系数三维模型,确定校正后的深度偏移层速度场。
9.根据权利要求8所述的装置,其特征在于,所述第一确定模块包括:
第一外推单元,用于在目标层位处为单个目标井时,将所述目标井处的比例系数沿所述地质构造模型外推,得到所述比例系数三维模型;
第一处理单元,用于在目标层位为多个目标井时,通过反距离加权插值法确定所述多个目标井之间的比例系数,将所述多个目标井处的比例系数以及所述多个目标井之间的比例系数沿所述地质构造模型外推,得到所述比例系数三维模型。
10.根据权利要求8所述的装置,其特征在于,所述计算模块包括:
运算单元,用于将所述VSP层速度与所述初始深度偏移速度做相除运算;
第二处理单元,用于将相除运算得到的值,作为所述比例系数。
11.根据权利要求8所述的装置,其特征在于,所述建立模块包括:
第一获取单元,用于获取所述目标井处的初始深度偏移层速度;
确定单元,用于根据所述初始深度偏移层速度和所述初始深度偏移速度,确定速度控制层层位数据;
内插单元,用于对所述速度控制层层位数据进行沿层内插,得到目标井处地质构造模型。
12.根据权利要求8所述的装置,其特征在于,还包括:
第二获取单元,用于获取所述目标井的各向异性参数;
第一建立单元,用于根据所述各向异性参数和所述地质构造模型,建立各向异性三维数据模型;
第二建立单元,用于根据所述各向异性三维数据模型和校正后的深度偏移层速度场,建立各向异性深度偏移层速度模型。
13.根据权利要求12所述的装置,其特征在于,所述第二建立单元包括:
第二外推单元,用于在目标层位处为单个目标井时,将所述目标井处各向异性参数沿所述地质构造模型外推,得到所述各向异性三维数据模型;
第三处理单元,用于在目标层位处为多个目标井时,通过反距离加权插值法确定所述多个目标井之间的各向异性参数,将所述多个目标井处各向异性参数以及所述多个目标井之间的各向异性参数沿所述地质构造模型外推,得到所述各向异性三维数据模型。
14.一种速度场校正设备,包括处理器以及用于存储处理器可执行指令的存储器,所述处理器执行所述指令时实现权利要求1至7中任一项所述方法的步骤。
15.一种计算机可读存储介质,其上存储有计算机指令,所述指令被执行时实现权利要求1至7中任一项所述方法的步骤。
CN201910148879.7A 2019-02-28 2019-02-28 速度场校正方法、装置和设备 Active CN111624662B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910148879.7A CN111624662B (zh) 2019-02-28 2019-02-28 速度场校正方法、装置和设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910148879.7A CN111624662B (zh) 2019-02-28 2019-02-28 速度场校正方法、装置和设备

Publications (2)

Publication Number Publication Date
CN111624662A CN111624662A (zh) 2020-09-04
CN111624662B true CN111624662B (zh) 2023-02-10

Family

ID=72269805

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910148879.7A Active CN111624662B (zh) 2019-02-28 2019-02-28 速度场校正方法、装置和设备

Country Status (1)

Country Link
CN (1) CN111624662B (zh)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117849864A (zh) * 2022-09-30 2024-04-09 中国石油天然气股份有限公司 随钻vsp井驱地震成像方法及装置

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102967882B (zh) * 2012-11-16 2015-08-05 中国石油集团川庆钻探工程有限公司地球物理勘探公司 地层的层速度模型的建模方法
WO2015196779A1 (zh) * 2014-06-24 2015-12-30 夏正元 利用叠加速度求取高精度地震波速度的方法
CN105842736B (zh) * 2016-05-27 2018-07-10 中国石油集团东方地球物理勘探有限责任公司 一种构建地层速度模型的方法
CN109143375B (zh) * 2018-08-20 2020-06-09 中国石油天然气集团有限公司 一种确定钻井方位的方法及装置

Also Published As

Publication number Publication date
CN111624662A (zh) 2020-09-04

Similar Documents

Publication Publication Date Title
US7840394B2 (en) Method for generating a 3D earth model
US9081918B2 (en) Methods and systems regarding models of underground formations
US20150370934A1 (en) Completion design based on logging while drilling (lwd) data
US9182510B2 (en) Methods and systems of incorporating pseudo-surface pick locations in seismic velocity models
EP3047096A1 (en) Identifying geological formation depth structure using well log data
CA3032780C (en) Geostatistical analysis of microseismic data in fracture modeling
EP2862008B1 (en) Stratigraphic modeling using production data density profiles
AU2011374333B2 (en) Method and system of correlating a measured log to a predicted log
CN111624662B (zh) 速度场校正方法、装置和设备
CN118091786B (zh) 基于井下测井资料的裂缝线密度校正方法和装置
CN118091786A (zh) 基于井下测井资料的裂缝线密度校正方法和装置
CN115951408A (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