CN109991663B - 工区地震速度体校正方法和装置 - Google Patents
工区地震速度体校正方法和装置 Download PDFInfo
- Publication number
- CN109991663B CN109991663B CN201711478060.4A CN201711478060A CN109991663B CN 109991663 B CN109991663 B CN 109991663B CN 201711478060 A CN201711478060 A CN 201711478060A CN 109991663 B CN109991663 B CN 109991663B
- Authority
- CN
- China
- Prior art keywords
- velocity
- seismic
- well point
- layer
- velocity gradient
- 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
Links
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/28—Processing seismic data, e.g. analysis, for interpretation, for correction
- G01V1/36—Effecting static or dynamic corrections on records, e.g. correcting spread; Correlating seismic signals; Eliminating effects of unwanted energy
- G01V1/362—Effecting static or dynamic corrections; Stacking
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V2210/00—Details of seismic processing or analysis
- G01V2210/50—Corrections or adjustments related to wave propagation
- G01V2210/52—Move-out correction
Abstract
本发明提供一种工区地震速度体校正方法和装置,通过获取与目的层和多个控制层中每层对应的第一地震层速度,并根据第一地震层速度计算得到目的层和各控制层的第一地震层速度梯度;然后获取井点位置的井点速度梯度系数比,对第一地震层速度梯度进行校正,获得校正后的第二地震层速度梯度;接着,以第二地震层速度梯度对第一地震层速度进行校正,获得校正后的第二地震层速度;最后,根据第二地震层速度与目的层和各个控制层的时间层位信息,得到校正后的工区地震速度体,结合了目的层和各控制层的井点参考速度,确定地震速度数据体平面梯度分布规律,有效地校正地震层速度平面趋势,提高台盆区地震工区速度体精度。
Description
技术领域
本发明涉及石油地震勘探技术,尤其涉及一种工区地震速度体校正方法和装置。
背景技术
在石油地震勘探过程中,台盆区的低幅度构造落实具有幅度低、面积小、埋藏深等特点,随着勘探开发的深入,构造落实难度加大,精度要求越来越高。
目前对速度体的各速度平面误差校正的方法,主要是利用各种网格算法对离散不规则分布的井点地震速度误差点进行网格化,得到误差平面图进行地震速度平面的校正。
由于井间及无井区地区的校正量随着不同网格算法变化而变化,在非井点地区的地震速度规律往往存在异常,导致非井点地区速度梯度的校正量准确不高。
发明内容
本发明提供一种工区地震速度体校正方法和装置,在保证井点处速度梯度准确的基础上提高了非井点处的速度梯度的准确性,从而提高了工区地震速度体的校正准确性。
根据本发明的第一方面,提供一种工区地震速度体校正方法,包括:
获取与目的层和多个控制层中每层对应的第一地震层速度,并根据所述第一地震层速度计算得到所述目的层和各控制层的第一地震层速度梯度;
获取井点位置的井点速度梯度系数比,并以所述井点速度梯度系数比,对所述第一地震层速度梯度进行校正,获得校正后的第二地震层速度梯度,其中,所述井点速度梯度系数比为所述第一地震层速度梯度与井点参考速度梯度之比,所述井点参考速度梯度为根据各井点在所述目的层和各控制层的井点参考速度确定的速度梯度;
以所述第二地震层速度梯度对所述第一地震层速度进行校正,获得校正后的第二地震层速度;
根据所述第二地震层速度与所述目的层和各个控制层的时间层位信息,得到校正后的工区地震速度体。
作为一种实现方式,所述根据第一地震层速度计算得到所述目的层和各控制层的第一地震层速度梯度,包括:
根据所述目的层和各控制层的第一地震层速度,得到同一个层速度平面各点位置处的地震速度值;
根据以下公式一得到目的层和各控制层的第一地震层速度梯度;
▽v地震=dv地震/dl公式一
其中,▽v地震为一个层速度平面的第一地震层速度梯度,dv地震为所述层速度平面各点之间的地震速度值变化量,dl为所述层速度平面各点的位置距离。
作为一种实现方式,所述获取井点位置的井点速度梯度系数比,包括:
获取各井点在所述目的层和各控制层的井点参考速度;
根据以下公式二确定井点参考速度梯度;
▽v井=dv井/dl公式二
其中,▽v井为一个层速度平面的井点参考速度梯度,dv井为所述层速度平面各井点之间的所述井点参考速度的变化量,dl为所述层速度平面各井点的位置距离;
根据以下公式三得到井点位置的井点速度梯度系数比;
r=(▽v地震/▽v井)公式三
其中,r为井点位置的井点速度梯度系数比,▽v地震为所述井点位置的第一地震层速度梯度,▽v井为所述井点位置的井点参考速度梯度。
作为一种实现方式,所述以所述井点速度梯度系数比,对所述第一地震层速度梯度进行校正,获得校正后的第二地震层速度梯度,包括:
对各井点的所述井点速度梯度系数比在所述目的层和各控制层上进行插值处理,得到与所述目的层和各控制层对应的速度梯度系数比平面,所述梯度系数比平面中的值为所述目的层和各控制层的对应点的速度梯度系数比;
用所述目的层和各控制层对应的第一地震层速度梯度,除以所述目的层和各控制层对应的速度梯度系数比平面中的值,得到所述目的层和各控制层对应的速度梯度校正平面;
将在井点处的所述第一地震层速度梯度,确定为所述井点处的第二地震层速度梯度,并将在非井点处的所述速度梯度校正平面对应的速度梯度,确定为所述非井点处的第二地震层速度梯度;
以所述第二地震层速度梯度对所述第一地震层速度进行校正,获得校正后的第二地震层速度。
作为一种实现方式,所述以所述第二地震层速度梯度对所述第一地震层速度进行校正,获得校正后的第二地震层速度,包括:
根据所述第二地震层速度梯度和各井点处的所述第一地震层速度,获得所述目的层和各控制层对应的速度校正平面;
将在井点处的所述第一地震层速度,确定为所述井点处的第二地震层速度;
将在非井点处的所述速度校正平面的地震层速度,确定为所述非井点处的第二地震层速度。
作为一种实现方式,在所述获取井点位置的井点速度梯度系数比之前,还包括:
获取所述第一地震层速度在各井点处的地震速度值;
获取各井点在所述目的层和各控制层的井点参考速度
判断所述地震速度值与所述井点参考速度在所述各井点处的变化趋势是否一致,若不一致,则执行所述获取井点位置的井点速度梯度系数比。
作为一种实现方式,所述井点参考速度包括测井速度值,或垂直地震剖面VSP速度值,或钻井分层速度值;
所述获取各井点在所述目的层和各控制层的井点参考速度,包括:
以声波测井法获得各井点在所述目的层和各控制层的测井速度值;或者
从VSP资料中获得各井点在所述目的层和各控制层的VSP速度值;或者
从钻井分层资料获得钻井分层深度值Δh,获取所述目的层和各控制层的的双程旅行时T0,获得各井点在所述目的层和各控制层的钻井分层速度值v=Δh/Δt0*2000。
作为一种实现方式,所述获取与目的层和多个控制层中每层对应的第一地震层速度,包括:
在工区内确定目的层和位于所述目的层之上的多个控制层;
获取与目的层和多个控制层中每层对应的深度域的原始速度体;
对所述深度域的原始速度体进行深时转换,得到时间域的原始速度体以及用于指示所述目的层和各个所述控制层的时间层位信息的等T0图;
从所述时间域的原始速度体中,得到与目的层和多个控制层中每层对应的第一地震层速度。
根据本发明的第二方面,提供一种工区地震速度体校正装置,包括:
第一地震层速度梯度获得模块,用于获取与目的层和多个控制层中每层对应的第一地震层速度,并根据所述第一地震层速度计算得到所述目的层和各控制层的第一地震层速度梯度;
第一校正模块,用于获取井点位置的井点速度梯度系数比,并以所述井点速度梯度系数比,对所述第一地震层速度梯度进行校正,获得校正后的第二地震层速度梯度,其中,所述井点速度梯度系数比为所述第一地震层速度梯度与井点参考速度梯度之比,所述井点参考速度梯度为根据各井点在所述目的层和各控制层的井点参考速度确定的速度梯度;
第二校正模块,用于以所述第二地震层速度梯度对所述第一地震层速度进行校正,获得校正后的第二地震层速度;
工区地震速度体建立模块,用于根据所述第二地震层速度与所述目的层和各个控制层的时间层位信息,得到校正后的工区地震速度体。
作为一种实现方式,所述第一地震层速度梯度获得模块具体用于:
根据所述目的层和各控制层的第一地震层速度,得到同一个层速度平面各点位置处的地震速度值;
根据以下公式一得到目的层和各控制层的第一地震层速度梯度;
▽v地震=dv地震/dl公式一
其中,▽v地震为一个层速度平面的第一地震层速度梯度,dv地震为所述层速度平面各点之间的地震速度值变化量,dl为所述层速度平面各点的位置距离。
本发明提供一种工区地震速度体校正方法和装置,通过获取与目的层和多个控制层中每层对应的第一地震层速度,并根据第一地震层速度计算得到目的层和各控制层的第一地震层速度梯度,从而得到一个待校正的对象;然后,获取井点位置的井点速度梯度系数比,并以井点速度梯度系数比,对第一地震层速度梯度进行校正,提高了校正的准确性,获得校正后的第二地震层速度梯度,其中,井点速度梯度系数比为第一地震层速度梯度与井点参考速度梯度之比,井点参考速度梯度为根据各井点在目的层和各控制层的井点参考速度确定的速度梯度;接着,以第二地震层速度梯度对第一地震层速度进行校正,获得校正后的第二地震层速度,能够在井点处以及非井点处都获得较准确的第二地震层速度;最后,根据第二地震层速度与目的层和各个控制层的时间层位信息,得到校正后的工区地震速度体,结合了目的层和各控制层的井点参考速度,确定地震速度数据体平面梯度分布规律,有效地校正地震层速度平面趋势,提高台盆区地震工区速度体精度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图做简单地介绍,显而易见地,下面描述中的附图是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以根据这些附图获得其他的附图。
图1为本发明实施例提供的一种工区地震速度体校正方法流程示意图;
图2为本发明实施例提供的一种目的层的第一地震层速度示意图;
图3为本发明实施例提供的一种控制层的第一地震层速度梯度示意图;
图4为本发明实施例提供的一种控制层的第二地震层速度梯度示意图;
图5为本发明实施例提供的一种目的层的第二地震层速度示意图;
图6是本发明实施例提供的一种地震速度值与井点参考速度在多个井点处的趋势进行比较示例;
图7是本发明实施例提供的一种工区地震速度体校正装置结构示意图。
具体实施方式
为使本发明实施例的目的、技术方案和优点更加清楚,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
本发明的说明书和权利要求书及上述附图中的术语“第一”、“第二”、“第三”“第四”等(如果存在)是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本发明的实施例能够以除了在这里图示或描述的那些以外的顺序实施。
应当理解,在本申请的各种实施例中,各过程的序号的大小并不意味着执行顺序的先后,各过程的执行顺序应以其功能和内在逻辑确定,而不应对本申请实施例的实施过程构成任何限定。
应当理解,在本申请中,“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
应当理解,在本申请中,“多个”是指两个或两个以上。“和/或”仅仅是一种描述关联对象的关联关系,表示可以存在三种关系,例如,和/或B,可以表示:单独存在A,同时存在A和B,单独存在B这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。
应当理解,在本申请中,“与A相应的B”表示B与A相关联,根据A可以确定B。根据A确定B并不意味着仅仅根据A确定B,还可以根据A和/或其他信息确定B。A与B的匹配,是A与B的相似度大于或等于预设的阈值。
取决于语境,如在此所使用的“若”可以被解释成为“在……时”或“当……时”或“响应于确定”或“响应于检测”。
下面以具体地实施例对本发明的技术方案进行详细说明。下面这几个具体的实施例可以相互结合,对于相同或相似的概念或过程可能在某些实施例不再赘述。
图1为本发明实施例提供的一种工区地震速度体校正方法流程示意图。图1所示的方法流程可以是由工区地震速度体校正装置实现的,该装置可以是专门的校正设备,也可以是包含于工区管理系统中的内置模块。其获取参数的方式可以是从预存储的资料中获取,也可以是通过相关的探测仪器、采集模块进行测量获取的。图1所示方法可以包括以下步骤:
S101,获取与目的层和多个控制层中每层对应的第一地震层速度,并根据所述第一地震层速度计算得到目的层和各控制层的第一地震层速度梯度。
地震层速度是在层状地层中,地震波在每一个划分层中传播的速度。不同结构的地层具有不同的地震层速度,它直接反映地层的岩性,还能通过地震层速度来划分地层。例如本实施例中的目的层和多个控制层。相应地,与目的层和多个控制层中每层对应的第一地震层速度,具体可以理解为分别是目的层中地震波的传播速度,和多个控制层中地震波的传播速度。
图2为本发明实施例提供的一种目的层的第一地震层速度示意图。图2中圆圈均表示井点,并以字符A、B、C、D、E和F进行区别。其中,在A和B井点左侧测得的第一地震层速度为3358m/s,在E井点左侧测得的第一地震层速度为3301m/s。图2中的曲线用于示意速度等值线。图2中由多条速度等值线构成了目的层的层速度平面。黑色箭头指示了第一地震层速度的减小方向。
目的层可以理解为需要进行油气开采的地层,多个控制层可以理解为位于目的层之上的比较重要的地层。为了给油气开采做准备,不仅仅需要对目的层进行分析,还需要考虑控制层的结构情况。
第一地震层速度的一种获取途径,可以是从预先存储的地层资料中获取,例如,预先通过地震测井或声波测井,测量得到包含地层速度谱的测量结果,并根据地震层速度将地层解释为目的层、多个控制层以及不重要的其他层,每层对应有第一地震层速度。预先测量的信息都存储于资料库中,在进行工区地震速度体校正时,直接从资料库中预存储的资料中获取与目的层和多个控制层中每层对应的第一地震层速度。
第一地震层速度的另一种获取途径,可以是实地测量。例如向测量子模块发送测量指令,以使得测量子模块开始对工区地层进行地震测井或声波测井测得与目的层和多个控制层中每层对应的第一地震层速度。
图3为本发明实施例提供的一种控制层的第一地震层速度梯度示意图。
在获得第一地震层速度后,通过预设的公式可以计算得到目的层和各控制层的第一地震层速度梯度。第一地震层速度梯度反映了在各层中不同位置点处第一地震层速度的变化快慢。例如图3示出的一个控制层的第一地震层速度梯度示意图中,圆圈均表示井点,并以字符A、B、C、D、E和F进行区别,曲线用于示意速度梯度等值线。
S102,获取井点位置的井点速度梯度系数比,并以井点速度梯度系数比,对第一地震层速度梯度进行校正,获得校正后的第二地震层速度梯度。
其中,井点速度梯度系数比为第一地震层速度梯度与井点参考速度梯度之比,井点参考速度梯度为根据各井点在目的层和各控制层的井点参考速度确定的速度梯度。
具体地,先根据各井点在目的层和各控制层的井点参考速度,确定井点参考速度梯度,然后用第一地震层速度梯度与井点参考速度梯度之比,作为井点速度梯度系数比。井点速度梯度系数比仍然是由井点数据得到的速度梯度比。由于第一地震层速度梯度是连续的梯度平面数据,不是离散的点值。因此,在以井点速度梯度系数比,对第一地震层速度梯度进行校正之前,还需要对井点速度梯度系数比使用最小二乘法等网格算法进行距离加权网格,进行插值和平滑,从离散的数据变为连续的数据。得到与第一地震层速度梯度类似分布的梯度平面后,再对第一地震层速度梯度进行校正。
井点参考速度可以理解为是测井速度值,或垂直地震剖面(VerticalSeismicProfiling,简称VSP)速度值,或者钻井分层速度值。
测井速度值是在测井过程中获取的各井点在目的层和各控制层的井点参考速度。石油钻井时,在钻到设计井深深度后都必须进行测井,以获得各种石油地质及工程技术资料,作为完井和开发油田的原始资料。
在一种实现方式中,井点参考速度为测井速度值,获取各井点在所述目的层和各控制层的井点参考速度可以理解为:以声波测井法获得各井点在所述目的层和各控制层的测井速度值。
在另一种实现方式中,井点参考速度为钻井分层速度值,获取各井点在所述目的层和各控制层的井点参考速度可以理解为:从钻井分层资料中获得各井点在所述目的层和各控制层之间的钻井分层深度值Δh,与从地震资料解释过程中得到的目的层和各控制层的双程旅行时值Δt0,获得钻井分层速度值v=Δh/Δt0*2000。
在再一种实现方式中,井点参考速度为VSP速度值,获取各井点在所述目的层和各控制层的井点参考速度可以理解为:从VSP资料中获得各井点在所述目的层和各控制层的VSP速度值。VSP速度值是通过VSP获得的各井点在目的层和各控制层的井点参考速度。VSP与通常地面观测的地震剖面相对应,是在地表附近的一些点上激发地震波,在沿井孔不同深度布置的一些多级多分量的检波点上进行观测。
图4为本发明实施例提供的一种控制层的第二地震层速度梯度示意图。图4中,圆圈均表示井点,并以字符A、B、C、D、E和F进行区别,曲线用于示意速度梯度等值线。图4是图3所示梯度示意图校正后的示例。
S103,以第二地震层速度梯度对第一地震层速度进行校正,获得校正后的第二地震层速度。
第二地震层速度梯度是经过校正的梯度,以校正后的梯度结合井点处的第一地震层速度对非井点处的第一地震层速度进行校正,可以得到校正后的速度。具体地,以每个井点处的第一地震层速度作为定值,以第二地震层速度梯度作为速度变化率,可以得到非井点处的第二地震层速度。井点处的第二地震层速度与其第一地震层速度一致。
图5为本发明实施例提供的一种目的层的第二地震层速度示意图。图5中,圆圈均表示井点,并以字符A、B、C、D、E和F进行区别,曲线用于示意速度梯度等值线。图5是图2所示速度示意图校正后的示例。
S104,根据第二地震层速度与目的层和各个控制层的时间层位信息,得到校正后的工区地震速度体。
时间层位信息可以理解为等T0图。等T0图可由时间剖面的数据直接绘制出,可以反映构造的基本形态,其偏移也小。等T0图也是由野外地震采集及室内地震数据处理后的深度数据进行深时转换得到的。
第二地震层速度是与目的层和各个控制层对应的速度平面,而结合目的层和各个控制层的时间层位信息,可以得到工区地震速度体。
本发明实施例提供一种工区地震速度体校正方法,通过获取与目的层和多个控制层中每层对应的第一地震层速度,并根据第一地震层速度计算得到目的层和各控制层的第一地震层速度梯度,从而得到一个待校正的对象;然后,获取井点位置的井点速度梯度系数比,并以井点速度梯度系数比,对第一地震层速度梯度进行校正,提高了校正的准确性,获得校正后的第二地震层速度梯度,其中,井点速度梯度系数比为第一地震层速度梯度与井点参考速度梯度之比,井点参考速度梯度为根据各井点在目的层和各控制层的井点参考速度确定的速度梯度;接着,以第二地震层速度梯度对第一地震层速度进行校正,获得校正后的第二地震层速度,能够在井点处以及非井点处都获得较准确的第二地震层速度;最后,根据第二地震层速度与目的层和各个控制层的时间层位信息,得到校正后的工区地震速度体,结合了目的层和各控制层的井点参考速度,确定地震速度数据体平面梯度分布规律,有效地校正地震层速度平面趋势,提高台盆区地震工区速度体精度。
图6是本发明实施例提供的一种地震速度值与井点参考速度在多个井点处的趋势进行比较示例。以下结合图6,对地震速度分布是否精确,是否需要进行校正提供了一种评价和判断的手段。在图1所示实施例的基础上,在获取井点位置的井点速度梯度系数比之前,还可以先执行是否需要进行校正的判断步骤,具体地:首先获取第一地震层速度在各井点处的地震速度值,并获取各井点在目的层和各控制层的井点参考速度,由此得到的地震速度值和井点参考速度都是在井点处的点值,以使得相同维度的数值可以进行比较。判断地震速度值与井点参考速度在各井点处的变化趋势是否一致,例如以图5所示的直方图对地震速度值与井点参考速度在多个井点处的趋势进行比较。每一组柱体从左至右依次表示在一处井点的地震速度值、测井速度值和VSP速度值。其中,在第四组合第五组未加入VSP速度值。两条虚线表示变化趋势,Q1为测井速度值的变化趋势,Q2为地震速度值的变化趋势。从图上可以看出,在第一组至第三组之间,Q1呈现先下降后上升的变化趋势,而Q2仅呈现上升的变化趋势。可见,图5所示地震速度值与井点参考速度在各井点处的变化趋势不一致。
若地震速度值与井点参考速度在各井点处的变化趋势不一致,表明工区地震速度体的误差较大,则执行获取井点位置的井点速度梯度系数比;若地震速度值与井点参考速度在各井点处的变化趋势一致,表明工区地震速度体的误差可以忽略不计,则无需对工区地震速度体进行校正。
在上述实施例的基础上,根据第一地震层速度计算得到所述目的层和各控制层的第一地震层速度梯度的过程,可以是:
首先,根据目的层和各控制层的第一地震层速度,得到同一个层速度平面各点位置处的地震速度值。
然后,根据以下公式一得到目的层和各控制层的第一地震层速度梯度;
▽v地震=dv地震/dl 公式一
其中,▽v地震为一个层速度平面的第一地震层速度梯度,dv地震为层速度平面各点之间的地震速度值变化量,dl为层速度平面各点的位置距离。
由于一个工区内地下地质情况从整体上变化不是十分剧烈,各地层岩性变化的平缓程度决定着地层层速度的变化梯度,结合构造地质因素,在已钻井点分布合理和断层断裂强度较小的情况下,本实施例用已知井点的速度梯度规律对地层的层速度梯度变化进行约束。
又因为造成目前井控处理使用的地震速度数据体分布异常的原因,在一定程度上体现在目的层以上的控制层的层速度在无井处的地震层速度梯度分布异常导致,本发明实施例采用体现井上速度梯度分布规律的井点速度梯度系数比,对地震速度梯度分布规律进行校正,得到更为精确的三维地震速度场。
在上述实施例的基础上,在得到第一地震层速度梯度之后,获取井点位置的井点速度梯度系数比的过程,可以是:
第一步,获取各井点在目的层和各控制层的井点参考速度。
第二步,根据以下公式二确定井点参考速度梯度。
▽v井=dv井/dl 公式二
其中,▽v井为一个层速度平面的井点参考速度梯度,dv井为层速度平面各井点之间的井点参考速度的变化量,dl为层速度平面各井点的位置距离。
第三步,根据以下公式三得到井点位置的井点速度梯度系数比。
r=(▽v地震'/▽v井') 公式三
其中,r为井点位置的井点速度梯度系数比,▽v地震'为井点位置的第一地震层速度梯度,▽v井'为井点位置的井点参考速度梯度。
在上述实施例的基础上,以井点速度梯度系数比,对第一地震层速度梯度进行校正,获得校正后的第二地震层速度梯度的过程,可以是:
首先,对各井点的井点速度梯度系数比在目的层和各控制层上进行插值处理,得到与目的层和各控制层对应的速度梯度系数比平面,梯度系数比平面中的值为目的层和各控制层的对应点的速度梯度系数比。其中,插值处理可以是在规则网格处理,或非规则网格处理的基础上,用距离加权法、方位法、趋势面法、叠加法、加权最小二乘法和克里金法等进行网格数值化,最后生成等值线形成速度梯度系数比平面。非规则网格处理例如可以是三角网格处理、任意四边形网格处理等。插值处理例如可以是最小二乘法等网格算法进行距离加权网格计算,从而对各井点的井点速度梯度系数比平滑处理得到速度梯度系数比平面。各井点的井点速度梯度系数比是在各井点处的离散的数据点,而梯度系数比平面是在目的层和各控制层上的连续的数据平面。
然后,用目的层和各控制层对应的第一地震层速度梯度,除以目的层和各控制层对应的速度梯度系数比平面中的值,得到目的层和各控制层对应的速度梯度校正平面。具体地,目的层和各控制层对应的第一地震层速度梯度和目的层和各控制层对应的速度梯度系数比平面都是由连续数据构成的平面,相处之后得到目的层和各控制层对应的速度梯度校正平面。
再然后,将在井点处的第一地震层速度梯度,确定为井点处的第二地震层速度梯度,并将在非井点处的速度梯度校正平面对应的速度梯度,确定为非井点处的第二地震层速度梯度。
最后,以第二地震层速度梯度对第一地震层速度进行校正,获得校正后的第二地震层速度。具体地,可以是根据第二地震层速度梯度和各井点处的所述第一地震层速度,获得目的层和各控制层对应的速度校正平面。例如以各井点处的所述第一地震层速度为定点,并以第二地震层速度梯度作为速度变化率,则可以从定点向外推算出非井点处的地震速度。将在井点处的第一地震层速度,确定为井点处的第二地震层速度,并将在非井点处的速度校正平面的地震层速度,确定为非井点处的第二地震层速度。由此,获得了校正后的第二地震层速度。
在上述实施例的基础上个,获取与目的层和多个控制层中每层对应的第一地震层速度的过程,具体可以是:经过野外地震采集及室内地震数据处理后,再经过地震解释人员根据实际生产情况解释,首先,在工区内确定目的层和位于目的层之上的多个控制层。然后,获取与目的层和多个控制层中每层对应的深度域的原始速度体,通常最初的原始数据都是深度域数据。再然后,对深度域的原始速度体进行深时转换,得到时间域的原始速度体以及用于指示目的层和各个控制层的时间层位信息的等T0图。最后,从时间域的原始速度体中,得到与目的层和多个控制层中每层对应的第一地震层速度。
深时转换,是将地震数据从深度域转换到时间域的过程,依靠速度信息展示地层的构造情况。
本发明实施例分析了基于井控处理的速度场速度梯度分布异常情况,结合地下地层构造因素,利用已知井点的测井、VSP、分层资料得到的速度梯度对地震速度梯度进行校正,从而科学有效地校正原始地震速度梯度分布,切实提高构造落实精度。本发明实施例利用了井点处的测井、VSP、及钻井分层数据的速度梯度信息对地震速度进行校正,在校正中不仅利用了井点处的速度信息,还在井间及远离井的区域获得了控制手段,提高了速度体的精度。
图7是本发明实施例提供的一种工区地震速度体校正装置结构示意图。如图7所示的装置主要包括:
第一地震层速度梯度获得模块11,用于获取与目的层和多个控制层中每层对应的第一地震层速度,并根据所述第一地震层速度计算得到所述目的层和各控制层的第一地震层速度梯度.
第一校正模块12,用于获取井点位置的井点速度梯度系数比,并以所述井点速度梯度系数比,对所述第一地震层速度梯度进行校正,获得校正后的第二地震层速度梯度,其中,所述井点速度梯度系数比为所述第一地震层速度梯度与井点参考速度梯度之比,所述井点参考速度梯度为根据各井点在所述目的层和各控制层的井点参考速度确定的速度梯度。
第二校正模块13,用于以所述第二地震层速度梯度对所述第一地震层速度进行校正,获得校正后的第二地震层速度.
工区地震速度体建立模块14,用于根据所述第二地震层速度与所述目的层和各个控制层的时间层位信息,得到校正后的工区地震速度体。
图7所示实施例的工区地震速度体校正装置对应地可用于执行图1所示方法实施例中的步骤,其实现原理和技术效果类似,此处不再赘述。
在上述实施例的基础上,第一地震层速度梯度获得模块具体用于:
根据所述目的层和各控制层的第一地震层速度,得到同一个层速度平面各点位置处的地震速度值;
根据以下公式一得到目的层和各控制层的第一地震层速度梯度;
▽v地震=dv地震/dl 公式一
其中,▽v地震为一个层速度平面的第一地震层速度梯度,dv地震为所述层速度平面各点之间的地震速度值变化量,dl为所述层速度平面各点的位置距离。
最后应说明的是:以上各实施例仅用以说明本发明的技术方案,而非对其限制;尽管参照前述各实施例对本发明进行了详细的说明,本领域的普通技术人员应当理解:其依然可以对前述各实施例所记载的技术方案进行修改,或者对其中部分或者全部技术特征进行等同替换;而这些修改或者替换,并不使相应技术方案的本质脱离本发明各实施例技术方案的范围。
Claims (9)
1.一种工区地震速度体校正方法,其特征在于,包括:
获取与目的层和多个控制层中每层对应的第一地震层速度,并根据所述第一地震层速度计算得到所述目的层和各控制层的第一地震层速度梯度;
获取井点位置的井点速度梯度系数比,并以所述井点速度梯度系数比,对所述第一地震层速度梯度进行校正,获得校正后的第二地震层速度梯度,其中,井点位置的井点速度梯度系数比为井点位置的第一地震层速度梯度与井点位置的井点参考速度梯度之比,所述井点参考速度梯度为根据各井点在所述目的层和各控制层的井点参考速度确定的速度梯度;
以所述第二地震层速度梯度对所述第一地震层速度进行校正,获得校正后的第二地震层速度;
根据所述第二地震层速度与所述目的层和各个控制层的时间层位信息,得到校正后的工区地震速度体;
所述以所述井点速度梯度系数比,对所述第一地震层速度梯度进行校正,获得校正后的第二地震层速度梯度,包括:
对各井点的所述井点速度梯度系数比在所述目的层和各控制层上进行插值处理,得到与所述目的层和各控制层对应的速度梯度系数比平面,所述梯度系数比平面中的值为所述目的层和各控制层的对应点的速度梯度系数比;
用所述目的层和各控制层对应的第一地震层速度梯度,除以所述目的层和各控制层对应的速度梯度系数比平面中的值,得到所述目的层和各控制层对应的速度梯度校正平面;
将在井点处的所述第一地震层速度梯度,确定为所述井点处的第二地震层速度梯度,并将在非井点处的所述速度梯度校正平面对应的速度梯度,确定为所述非井点处的第二地震层速度梯度。
4.根据权利要求1所述的方法,其特征在于,所述以所述第二地震层速度梯度对所述第一地震层速度进行校正,获得校正后的第二地震层速度,包括:
根据所述第二地震层速度梯度和各井点处的所述第一地震层速度,获得所述目的层和各控制层对应的速度校正平面;
将在井点处的所述第一地震层速度,确定为所述井点处的第二地震层速度;
将在非井点处的所述速度校正平面的地震层速度,确定为所述非井点处的第二地震层速度。
5.根据权利要求1所述的方法,其特征在于,在所述获取井点位置的井点速度梯度系数比之前,还包括:
获取所述第一地震层速度在各井点处的地震速度值;
获取各井点在所述目的层和各控制层的井点参考速度
判断所述地震速度值与所述井点参考速度在所述各井点处的变化趋势是否一致,若不一致,则执行所述获取井点位置的井点速度梯度系数比。
6.根据权利要求3或5所述的方法,其特征在于,所述井点参考速度包括测井速度值,或垂直地震剖面VSP速度值,或钻井分层速度值;
所述获取各井点在所述目的层和各控制层的井点参考速度,包括:
以声波测井法获得各井点在所述目的层和各控制层的测井速度值;或者
从VSP资料中获得各井点在所述目的层和各控制层的VSP速度值;或者
从钻井分层资料获得钻井分层深度值Δh,获取所述目的层和各控制层的的双程旅行时值Δt0,获得各井点在所述目的层和各控制层的钻井分层速度值v=Δh/Δt0*2000。
7.根据权利要求1所述的方法,其特征在于,所述获取与目的层和多个控制层中每层对应的第一地震层速度,包括:
在工区内确定目的层和位于所述目的层之上的多个控制层;
获取与目的层和多个控制层中每层对应的深度域的原始速度体;
对所述深度域的原始速度体进行深时转换,得到时间域的原始速度体以及用于指示所述目的层和各个所述控制层的时间层位信息的等T0图;
从所述时间域的原始速度体中,得到与目的层和多个控制层中每层对应的第一地震层速度。
8.一种工区地震速度体校正装置,其特征在于,包括:
第一地震层速度梯度获得模块,用于获取与目的层和多个控制层中每层对应的第一地震层速度,并根据所述第一地震层速度计算得到所述目的层和各控制层的第一地震层速度梯度;
第一校正模块,用于获取井点位置的井点速度梯度系数比,并以所述井点速度梯度系数比,对所述第一地震层速度梯度进行校正,获得校正后的第二地震层速度梯度,其中,井点位置的井点速度梯度系数比为井点位置的第一地震层速度梯度与井点位置的井点参考速度梯度之比,所述井点参考速度梯度为根据各井点在所述目的层和各控制层的井点参考速度确定的速度梯度;
第二校正模块,用于以所述第二地震层速度梯度对所述第一地震层速度进行校正,获得校正后的第二地震层速度;
工区地震速度体建立模块,用于根据所述第二地震层速度与所述目的层和各个控制层的时间层位信息,得到校正后的工区地震速度体;
所述以所述井点速度梯度系数比,对所述第一地震层速度梯度进行校正,获得校正后的第二地震层速度梯度,包括:
对各井点的所述井点速度梯度系数比在所述目的层和各控制层上进行插值处理,得到与所述目的层和各控制层对应的速度梯度系数比平面,所述梯度系数比平面中的值为所述目的层和各控制层的对应点的速度梯度系数比;
用所述目的层和各控制层对应的第一地震层速度梯度,除以所述目的层和各控制层对应的速度梯度系数比平面中的值,得到所述目的层和各控制层对应的速度梯度校正平面;
将在井点处的所述第一地震层速度梯度,确定为所述井点处的第二地震层速度梯度,并将在非井点处的所述速度梯度校正平面对应的速度梯度,确定为所述非井点处的第二地震层速度梯度。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711478060.4A CN109991663B (zh) | 2017-12-29 | 2017-12-29 | 工区地震速度体校正方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201711478060.4A CN109991663B (zh) | 2017-12-29 | 2017-12-29 | 工区地震速度体校正方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN109991663A CN109991663A (zh) | 2019-07-09 |
CN109991663B true CN109991663B (zh) | 2021-11-02 |
Family
ID=67109793
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201711478060.4A Active CN109991663B (zh) | 2017-12-29 | 2017-12-29 | 工区地震速度体校正方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109991663B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113064215A (zh) * | 2020-01-02 | 2021-07-02 | 中国石油天然气股份有限公司 | 构造图井校的方法、装置及系统 |
CN111624649B (zh) * | 2020-06-05 | 2022-05-20 | 中油奥博(成都)科技有限公司 | 一种零偏移距vsp建立横向变速层速度模型的方法和装置 |
Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106569260A (zh) * | 2015-10-10 | 2017-04-19 | 中国石油化工股份有限公司 | 针对高陡构造地区二维地震变速空校成图方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US4849887A (en) * | 1987-08-28 | 1989-07-18 | Amoco Corporation | Horizon velocity analysis |
US5285422A (en) * | 1992-10-23 | 1994-02-08 | Western Atlas International, Inc. | Method for compensating 3D DMO for the effects of an inhomogeneous earth |
CN102129084B (zh) * | 2010-12-17 | 2012-10-24 | 中国石油天然气股份有限公司 | 一种井控获取地震薄储层速度的方法及装置 |
CN102176053B (zh) * | 2011-01-27 | 2012-10-17 | 中国科学院地质与地球物理研究所 | 提升波动方程叠前深度偏移成像效果的方法 |
CN103543466B (zh) * | 2012-07-17 | 2016-05-04 | 中国石油化工股份有限公司 | 一种时间域地震层速度反演方法 |
CN102967882B (zh) * | 2012-11-16 | 2015-08-05 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 地层的层速度模型的建模方法 |
CN105182409B (zh) * | 2015-08-31 | 2018-09-04 | 中国石油天然气集团公司 | 一种速度模型建立方法 |
CN105842736B (zh) * | 2016-05-27 | 2018-07-10 | 中国石油集团东方地球物理勘探有限责任公司 | 一种构建地层速度模型的方法 |
-
2017
- 2017-12-29 CN CN201711478060.4A patent/CN109991663B/zh active Active
Patent Citations (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106569260A (zh) * | 2015-10-10 | 2017-04-19 | 中国石油化工股份有限公司 | 针对高陡构造地区二维地震变速空校成图方法 |
Also Published As
Publication number | Publication date |
---|---|
CN109991663A (zh) | 2019-07-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109425896B (zh) | 白云岩油气储层分布预测方法及装置 | |
US11209561B2 (en) | Generation of fault displacement vector and/or fault damage zone in subsurface formation using stratigraphic function | |
CN104047597B (zh) | 含油气泥页岩地层测井曲线标准化方法 | |
US5835882A (en) | Method for determining barriers to reservoir flow | |
US20160124116A1 (en) | Generation of structural elements for subsurface formation using stratigraphic implicit function | |
CN107783187B (zh) | 一种将测井速度和地震速度结合建立三维速度场的方法 | |
CA2823710C (en) | Methods and systems regarding models of underground formations | |
EA020635B1 (ru) | Способ сейсмической разведки для поиска углеводородов с применением поля средних скоростей сейсмической волны | |
WO1998034190A9 (en) | Method for determining distribution of reservoir permeability, porosity and pseudo relative permeability | |
EP0956533A1 (en) | Method for determining distribution of reservoir permeability, porosity and pseudo relative permeability | |
US10732310B2 (en) | Seismic attributes derived from the relative geological age property of a volume-based model | |
CN107329171A (zh) | 深度域储层地震反演方法及装置 | |
CN105259581B (zh) | 一种地震资料时深转换方法 | |
Kianoush et al. | Compressional and shear interval velocity modeling to determine formation pressures in an oilfield of SW Iran | |
US10451769B2 (en) | Method for determining petrophysical properties from logging measurements | |
CN112114380A (zh) | 一种用于砂砾岩油藏储层的预测方法 | |
CN109991663B (zh) | 工区地震速度体校正方法和装置 | |
CN112505754B (zh) | 基于高精度层序格架模型的井震协同划分沉积微相的方法 | |
CN108072901B (zh) | 一种获得准确的静水压力和上覆地层压力的方法及系统 | |
CN107589457A (zh) | 测井资料尺度粗化的方法和装置 | |
CN113267810B (zh) | 地震勘探全深度速度建模方法及装置 | |
CN113514904B (zh) | 地层参数模型建立方法及装置 | |
CN110954951B (zh) | 一种利用地震切片进行四级层序地层划分与识别的方法 | |
CN110927819B (zh) | 一种裂缝发育程度表征方法 | |
RU2201606C1 (ru) | Способ типизации и корреляции нефтегазопродуктивных горных пород по скважинным спектрально-временным параметрам |
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 |