CN111077589A - 干岩石纵横波速度比的确定方法和装置 - Google Patents

干岩石纵横波速度比的确定方法和装置 Download PDF

Info

Publication number
CN111077589A
CN111077589A CN201811224519.2A CN201811224519A CN111077589A CN 111077589 A CN111077589 A CN 111077589A CN 201811224519 A CN201811224519 A CN 201811224519A CN 111077589 A CN111077589 A CN 111077589A
Authority
CN
China
Prior art keywords
target
data
logging
longitudinal
interval
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
CN201811224519.2A
Other languages
English (en)
Other versions
CN111077589B (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201811224519.2A priority Critical patent/CN111077589B/zh
Publication of CN111077589A publication Critical patent/CN111077589A/zh
Application granted granted Critical
Publication of CN111077589B publication Critical patent/CN111077589B/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
    • G01V11/00Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00

Landscapes

  • Physics & Mathematics (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • General Physics & Mathematics (AREA)
  • Geophysics (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

本申请公开了一种干岩石纵横波速度比的确定方法和装置。该方法包括:确定与待确定的干岩石纵横波速度比相对应的目标层段,其中,目标层段位于目标测量井中;获取目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合,其中,目标测井数据经由声波测井和密度测井获取;依据目标测井数据集合中的目标测井数据,确定目标层段对应的干岩石纵横波速度比。通过本申请,解决了相关技术中无法便捷、准确地确定干岩石纵横波速度比的问题。

Description

干岩石纵横波速度比的确定方法和装置
技术领域
本申请涉及地质勘探领域,具体而言,涉及一种干岩石纵横波速度比的确定方法和装置。
背景技术
早期用地震振幅异常检测烃类的方法是一种定性技术,如地震上的“亮点”、“暗点”和“平点”技术等。随着岩性油气藏勘探的深入,这种简单的定性方法不再适用。而基于叠前地震资料直接进行流体预测方法,提供了一种定量的方法来区分不同的孔隙类型流体。Russell(2003)提出的流体属性因子可较好的区分储层中的孔隙流体,对流体识别具有重要意义。但如果流体属性因子中的比例系数,干岩石纵横波速比选择的不合适,流体属性因子就不仅和孔隙流体有关,还与岩石骨架岩性有关,这会给流体识别带来多解性。因此,干岩石纵横波速度比是流体属性因子能有效预测储层流体的关键因素。
但是,传统获得干岩石纵横波速度比的方法是采用经验值
Figure BDA0001835609700000011
或在实验室测量该值,其中,经验值一般偏离实际工区的实际值,基于此计算得到的流体属性因子指示孔隙流体效果差,水藏描述不清楚,且流体属性因子还包含岩性的信息;而实验室测量,不仅经济成本高、耗时长,而且难以克服测量过程中的不确定性和测量误差。
针对相关技术中无法便捷、准确地确定干岩石纵横波速度比的问题,目前尚未提出有效的解决方案。
发明内容
本申请提供一种干岩石纵横波速度比的确定方法和装置,以解决相关技术中无法便捷、准确地确定干岩石纵横波速度比的问题。
根据本申请的一个方面,提供了一种干岩石纵横波速度比的确定方法。该方法包括:确定与待确定的干岩石纵横波速度比相对应的目标层段,其中,所述目标层段位于目标测量井中;获取所述目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合,其中,所述目标测井数据经由声波测井和密度测井获取;依据所述目标测井数据集合中的目标测井数据,确定所述目标层段对应的干岩石纵横波速度比。
可选的,获取所述目标层段中多个采样点分别对应的目标测井数据包括:获取所述目标测量井中目标层段对应的测井总数据,其中,所述测井总数据由所述目标测量井进行声波测井处理和密度测井处理而获取的;判断所述目标层段对应的测井总数据中是否存在失真数据;在所述目标层段对应的测井总数据中不存在失真数据的情况下,则依据所述目标测量井对应的测井总数据,确定所述目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合;在所述目标层段对应的测井总数据中存在失真数据的情况下,则对所述测井总数据中的未失真数据进行多元线性拟合,确定所述未失真数据中多种数据之间的数学关系;以依据所述未失真数据中多种数据之间的数学关系,重新确定失真层段对应的测井计算数据,依据所述测井计算数据更替所述测井总数据中的失真数据,进而确定被更替所述失真数据后的测井总数据中不存在失真数据,其中,失真数据为所述失真层段对应的测井数据。
可选的,在所述目标层段对应的测井总数据包含:横波时差曲线、纵波时差曲线、以及密度曲线的情况下,判断所述目标层段对应的测井总数据中是否存在失真数据包括:获取所述目标层段对应的测井曲线,其中,所述测井曲线包含以下至少之一:自然电位曲线、井径曲线、自然伽马曲线、深电阻率曲线、中电阻率曲线、浅电阻率曲线、声波曲线、密度曲线、中子曲线;依据所述目标层段对应的测井曲线对所述目标层段对应的测井总数据进行分析,确定中所述目标层段对应的测井总数据是否存在失真数据。
可选的,在每个所述采样点对应的目标测井数据包括:所述采样点对应的纵波速度值、所述采样点对应的横波速度值、所述采样点对应的密度数值的情况下,依据所述目标测井数据集合中的目标测井数据,确定所述目标层段对应的干岩石纵横波速度比包括:依据所述目标测井数据集合中的目标测井数据,确定所述干岩石纵横波速度比对应的待计算值,其中,所述待计算值中包含:每个所述采样点对应的密度与纵波速度平方的乘积、每个所述采样点对应的密度与横波速度平方的乘积、所述目标层段包含的多个采样点对应的密度与纵波速度平方的乘积的平均值,所述目标层段包含的多个采样点对应的密度与横波速度平方的乘积的平均值、所述采样点的数量;依据所述干岩石纵横波速度比对应的待计算值确定所述目标层段对应的干岩石纵横波速度比。
可选的,在目标地区包含至少一个所述目标测量井的情况下,依据每个所述目标测量井对应的多个采样点分别对应的目标测井数据,确定所述目标地区对应的干岩石纵横波速度比包括:获取每个所述目标测量井对应多个采样点分别对应的目标测井数据与所述目标地区对应的干岩石纵横波速度比之间的关系公式:
Figure BDA0001835609700000031
其中,m为目标测量井的数目,n为目标层段中的采样点的数量,
Figure BDA0001835609700000032
为目标层段中第j采样点的密度与纵波速度平方的乘积,
Figure BDA0001835609700000033
为目标层段中第j采样点的密度与横波速度平方的乘积,
Figure BDA0001835609700000034
为包含的多个采样点对应的密度与纵波速度平方的乘积的平均值,
Figure BDA0001835609700000035
为所述目标层段包含的多个采样点对应的密度与横波速度平方的乘积的平均值,
Figure BDA0001835609700000036
为所述目标地区对应的干岩石纵横波速度比;依据所述关系公式和所述多个采样点分别对应的目标测井数据,确定所述目标地区对应的干岩石纵横波速度比。
可选的,在所述目标地区对应的干岩石纵横波速度比之后,所述方法还包括:依据所述目标地区对应的干岩石纵横波速度比和多个所述目标测量井中每个所述采样点对应的目标测井数据,确定多个所述目标测量井对应的流体属性因子;获取所述目标地区包含的每个目标测量井对应的测井解释,并依据多个所述目标测量井对应的流体属性因子和多个所述目标测量井对应的测井解释建立所述目标地区对应的识别模型,其中,所述识别模型用于确定所述目标地区的储层流体的位置。
根据本申请的另一方面,提供了一种干岩石纵横波速度比的确定装置。该装置包括:第一确定单元,用于确定与待确定的干岩石纵横波速度比相对应的目标层段,其中,所述目标层段位于目标测量井中;获取单元,用于获取所述目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合,其中,所述目标测井数据经由声波测井和密度测井获取;第二确定单元,用于依据所述目标测井数据集合中的目标测井数据,确定所述目标层段对应的干岩石纵横波速度比。
可选的,所述获取单元包括:第一获取模块,用于获取所述目标测量井中目标层段对应的测井总数据,其中,所述测井总数据由所述目标测量井进行声波测井处理和密度测井处理而获取的;判断模块,用于判断所述目标层段对应的测井总数据中是否存在失真数据;第一确定模块,用于在所述目标层段对应的测井总数据中不存在失真数据的情况下,则依据所述目标测量井对应的测井总数据,确定所述目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合;重新确定模块,用于在所述目标层段对应的测井总数据中存在失真数据的情况下,则对所述测井总数据中的未失真数据进行多元线性拟合,确定所述未失真数据中多种数据之间的数学关系;以依据所述未失真数据中多种数据之间的数学关系,重新确定失真层段对应的测井计算数据,依据所述测井计算数据更替所述测井总数据中的失真数据,进而确定被更替所述失真数据后的测井总数据中不存在失真数据,其中,失真数据为所述失真层段对应的测井数据。
根据本申请的另一方面,提供了一种存储介质,所述存储介质包括存储的程序,其中,所述程序执行上述任意一项所述的干岩石纵横波速度比的确定方法。
根据本申请的另一方面,提供了一种处理器,所述处理器用于运行程序,其中,所述程序运行时执行上述任意一项所述的干岩石纵横波速度比的确定方法。
通过本申请,采用以下步骤:确定与待确定的干岩石纵横波速度比相对应的目标层段,其中,所述目标层段位于目标测量井中;获取所述目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合,其中,所述目标测井数据经由声波测井和密度测井获取;依据所述目标测井数据集合中的目标测井数据,确定所述目标层段对应的干岩石纵横波速度比,解决了相关技术中无法便捷、准确地确定干岩石纵横波速度比的问题。
需要说明的是:在常规测井工程中往往包含着对目标测量井进行声波测井处理和密度测井处理,基于此,本申请实施例通过获取目标测量井中目标层段对应的声波测井数据和密度测井数据,并基于该声波测井数据和密度测井数据,即可确定出目标测量井中目标层段对应的干岩石横纵波速度比,达到了快速且便捷的确定干岩石横纵波速度比的技术效果。
此外,相对于现有技术中的采用经验值确定干岩石横纵波速度比的技术方案,本申请实施例通过依据目标测量井中目标层段对应的具体的目标测量数据集合(目标层段中多个采样点分别对应的声波测井数据和密度测井数据),来确定该目标层段对应的干岩石横纵波速度比,达到了准确的确定干岩石横纵波速度比的技术效果。
附图说明
构成本申请的一部分的附图用来提供对本申请的进一步理解,本申请的示意性实施例及其说明用于解释本申请,并不构成对本申请的不当限定。在附图中:
图1是根据本申请实施例提供的干岩石纵横波速度比的确定方法的流程图;
图2是根据本申请实施例提供的一种可选的Marmousi模型的示意图;
图3是根据本申请实施例提供的一种可选的纵横波速度的剖面示意图;
图4是根据本申请实施例提供的一种可选的密度的剖面示意图;
图5是根据本申请实施例提供的一种可选的横波速度的剖面示意图;
图6是根据本申请实施例提供的一种可选的剪切模量的剖面示意图;
图7是根据本申请实施例提供的一种可选的横纵波速度比的剖面示意图;
图8是根据本申请实施例提供的一种可选的流体属性因子的剖面示意图一;
图9是根据本申请实施例提供的一种可选的流体属性因子的剖面示意图二;以及
图10是根据本申请实施例提供的干岩石纵横波速度比的确定装置的示意图。
具体实施方式
需要说明的是,在不冲突的情况下,本申请中的实施例及实施例中的特征可以相互组合。下面将参考附图并结合实施例来详细说明本申请。
为了使本技术领域的人员更好地理解本申请方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分的实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
需要说明的是,本申请的说明书和权利要求书及上述附图中的术语“第一”、“第二”等是用于区别类似的对象,而不必用于描述特定的顺序或先后次序。应该理解这样使用的数据在适当情况下可以互换,以便这里描述的本申请的实施例。此外,术语“包括”和“具有”以及他们的任何变形,意图在于覆盖不排他的包含,例如,包含了一系列步骤或单元的过程、方法、系统、产品或设备不必限于清楚地列出的那些步骤或单元,而是可包括没有清楚地列出的或对于这些过程、方法、产品或设备固有的其它步骤或单元。
根据本申请的实施例,提供了一种干岩石纵横波速度比的确定方法。
图1是根据本申请实施例的干岩石纵横波速度比的确定方法的流程图。如图1所示,该方法包括以下步骤:
步骤S102,确定与待确定的干岩石纵横波速度比相对应的目标层段,其中,目标层段位于目标测量井中。
步骤S104,获取目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合,其中,目标测井数据经由声波测井和密度测井获取。
步骤S106,依据目标测井数据集合中的目标测井数据,确定目标层段对应的干岩石纵横波速度比。
本申请实施例提供的干岩石纵横波速度比的确定方法,通过确定与待确定的干岩石纵横波速度比相对应的目标层段,其中,目标层段位于目标测量井中;获取目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合,其中,目标测井数据经由声波测井和密度测井获取;依据目标测井数据集合中的目标测井数据,确定目标层段对应的干岩石纵横波速度比,解决了相关技术中无法便捷、准确地确定干岩石纵横波速度比的问题。
需要说明的是:在常规测井工程中往往包含着对目标测量井进行声波测井处理和密度测井处理,基于此,本申请实施例通过获取目标测量井中目标层段对应的声波测井数据和密度测井数据,并基于该声波测井数据和密度测井数据,即可确定出目标测量井中目标层段对应的干岩石横纵波速度比,达到了快速且便捷的确定干岩石横纵波速度比的技术效果。
此外,相对于现有技术中的采用经验值确定干岩石横纵波速度比的技术方案,本申请实施例通过依据目标测量井中目标层段对应的具体的目标测量数据集合(目标层段中多个采样点分别对应的声波测井数据和密度测井数据),来确定该目标层段对应的干岩石横纵波速度比,达到了准确的确定干岩石横纵波速度比的技术效果。
在一种可选的示例中,目标层段中多个采样点分别对应的目标测井数据是由测井总数据确定的,其中,测井总数据是目标测量井中的目标层段进行声波测井处理和密度测井处理而获取的。
而为了保证测井总数据的能够真实表达目标测量井的属性参数,避免测井总数据中的失真数据影响干岩石横纵波速度比的确定结果,可选的,在本申请实施例提供的干岩石纵横波速度比的确定方法中,获取目标层段中多个采样点分别对应的目标测井数据包括:获取目标测量井中目标层段对应的测井总数据,其中,测井总数据由目标测量井进行声波测井处理和密度测井处理而获取的;判断目标层段对应的测井总数据中是否存在失真数据;在目标层段对应的测井总数据中不存在失真数据的情况下,则依据目标测量井对应的测井总数据,确定目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合;在目标层段对应的测井总数据中存在失真数据的情况下,则对测井总数据中的未失真数据进行多元线性拟合,确定未失真数据中多种数据之间的数学关系;以依据未失真数据中多种数据之间的数学关系,重新确定失真层段对应的测井计算数据,依据测井计算数据更替测井总数据中的失真数据,进而确定被更替失真数据后的测井总数据中不存在失真数据,其中,失真数据为失真层段对应的测井数据。
也即,考虑到在真实测井工程中,对目标测量井中的目标层段进行声波测井处理和密度测井处理时,获取的测井总数据中容易包含因各种因素导致的失真数据,因此,在获取到目标测量井中目标层段对应的测井总数据之后,对该测井总数据进行数据质量评估,也即,评估测井总数据中是否存在失真数据,若不存在失真数据则直接依据测井总数据确定目标层段中多个采样点分别对应的目标测井数据,若存在失真数据则依据未失真数据重新确定失真层段的测井计算数据,进而依据测井计算数据和未失真数据来确定目标层段中多个采样点分别对应的目标测井数据。
举例说明,若目标层段对应的测井总数据包含:横波时差曲线、纵波时差曲线、以及密度曲线时,则判断目标层段对应的测井总数据中是否存在失真数据包括:获取目标层段对应的测井曲线,其中,测井曲线包含以下至少之一:自然电位曲线、井径曲线、自然伽马曲线、深电阻率曲线、中电阻率曲线、浅电阻率曲线、声波曲线|(包含横波时差曲线、纵波时差曲线)、密度曲线、中子曲线;依据目标层段对应的测井曲线对目标层段对应的测井总数据进行分析,确定中目标层段对应的测井总数据是否存在失真数据。
也即,通过测井曲线综合分析的方式,对测井总数据进行质量控制,评估目标层段对应的测井总数据中是否包含失真数据。
举例说明,在实际测井过程中,测井总数据容易发生失真的情况,若继续使用失真的测井总数据计算目标层段中多个采样点分别对应的目标测井数据,则会降低最终的干岩石纵横波速度比的真实性。为了避免上述情况发生,在测井总数据中包含失真的情况下,依据失真数据确定目标层段中的失真层段和未失真层段,其中,失真层段对应的测井数据包含失真数据,而未失真层段对应的测井数据则不包含失真数据,基于未失真层段对应的测井数据确定测井数据中多种数据之间的数学关系,进而依据该数学关系和失真层段对应的未失真数据,重新确定失真层段对应的测井计算数据,进而使用测井计算数据更替测井总数据中的失真数据,达到令测井总数据中不存在失真数据的技术效果,其中,该测井计算数据与失真数据的数据种类相同,且测井计算数据与失真数据对应的目标测量井中的层段也相同。
也即,以测井总数据中包含:横波时差数据、纵波时差数据、以及密度数据为例,在测井总数据中包含一段失真的密度数据,此时,则确定该失真的密度数据对应的失真层段,以及不包含该失真的密度曲线数据的未失真层段,依据未失真层段的横波时差数据、纵波时差数据、以及密度数据,确定横波时差数据、纵波时差数据和密度数据之间的数学关系,此时依据该数学关系和失真层段对应的横波时差数据、纵波时差数据,确定该失真层段对应的密度数据,此时,使用重新确定的密度数据替换测井总数据中失真的密度数据,进而达到了令测井总数据中不包含失真数据的技术效果。
举例说明,若目标层段对应的测井总数据包含:横波时差曲线、纵波时差曲线、以及密度曲线,且目标层段对应的测井总数据不存在失真数据时,则依据测井总数据确定目标层段中多个采样点分别对应的目标测井数据包括:基于测井总数据包含的横波时差曲线和纵波时差曲线,通过取时差倒数的方式,确定目标层段中多个采样点分别对应的纵波速度和横波速度;以及基于测井总数据中包含的密度曲线,确定目标层段中多个采样点分别对应的密度值。
此外,也可以基于测井总数据包含的横波时差曲线和纵波时差曲线,通过取时差倒数的方式,获取目标层段对应的横波速度曲线和纵波速度曲线,进而基于该横波速度曲线和纵波速度曲线确定目标层段中多个采样点分别对应的横波速度和纵波速度。
进一步地,在依据目标测井数据集合中的目标测井数据,确定目标层段对应的干岩石纵横波速度比之前,还需要确定出目标测井数据集合中的目标测井数据与目标层段对应的干岩石纵横波速度比之间的数值关系,也即,确定多个采样点分别对应的目标测井数据与目标层段对应的干岩石纵横波之间的关系公式。
在目标测井数据包含纵波速度、横波速度和密度时,该关系公式的推导步骤具体如下:
基于Biot-Gassmann理论,在低频条件下,饱和岩石的拉梅系数λsat和体积模量Ksat可以表示为流体项和干岩石项的和:
λsat=λdry2M (1)
Ksat=Kdry2M (2)
μsat=μdry=μ (3)
其中,λdry和Kdry分别是干岩石的拉梅系数和体积模量,μsat和μdry分别是饱和岩石和干岩石的剪切模量,β是Biot系数,M表示能使水进入地层而不改变岩层体积形态所需的压力。
根据Biot(1941)和Gassmann(1951)方程,β和M可以表示为:
Figure BDA0001835609700000081
Figure BDA0001835609700000091
其中,Ksolid是固体颗粒的体积模量,φ是岩石的孔隙度,Kf是有效孔隙流体的体积模量。
根据弹性孔隙介质理论,地震波在岩石中传播时,纵横波速度可写为:
Figure BDA0001835609700000092
Figure BDA0001835609700000093
其中,ρsat是饱和岩石的密度。将公式(2)(3)代入(6)中可得:
Figure BDA0001835609700000094
Russell等(2003)重写了饱和岩石的纵横波速度:
Figure BDA0001835609700000095
Figure BDA0001835609700000096
其中,s是干骨架项,可表示为λdry+2μ或Kdry+(4/3)μ;f是混合流体/孔隙项,可表示为β2M。
联立公式(9)(10),流体属性因子f可表示为:
Figure BDA0001835609700000097
其中,
Figure BDA0001835609700000098
只与孔隙流体有关,而与岩石岩性无关。因此,可建立如下等式:
Figure BDA0001835609700000101
其中,r(f,μ)表示流体属性因子和剪切模量的互相关系数。f(j)和μ(j)分别为目标层段测井数据中第j点的流体属性因子和剪切模量值,
Figure BDA0001835609700000102
Figure BDA0001835609700000103
分别为流体属性因子和剪切模量的平均值,n为目标层测井数据的采样点数。
通过上述等式,可以推导出干岩石的纵横波速度比
Figure BDA0001835609700000104
Figure BDA0001835609700000105
其中,
Figure BDA0001835609700000106
Figure BDA0001835609700000107
分别表示目标层段测井数据中第j点的密度与纵波速度平方的乘积和密度与横波速度的乘积,
Figure BDA0001835609700000108
Figure BDA0001835609700000109
分别为目标层段测井数据密度与纵波速度平方的乘积和密度与横波速度平方的乘积的平均,n为目标层测井数据的采样点数。
基于上述推导出的公式(13)可知,在每个采样点对应的目标测井数据包括:采样点对应的纵波速度值、采样点对应的横波速度值、采样点对应的密度数值的情况下,依据目标测井数据集合中的目标测井数据,确定目标层段对应的干岩石纵横波速度比包括:依据目标测井数据集合中的目标测井数据,确定干岩石纵横波速度比对应的待计算值,其中,待计算值中包含:每个采样点对应的密度与纵波速度平方的乘积、每个采样点对应的密度与横波速度平方的乘积、目标层段包含的多个采样点对应的密度与纵波速度平方的乘积的平均值,目标层段包含的多个采样点对应的密度与横波速度平方的乘积的平均值、采样点的数量;依据干岩石纵横波速度比对应的待计算值确定目标层段对应的干岩石纵横波速度比。
进一步地,依据干岩石纵横波速度比对应的待计算值确定目标层段对应的干岩石纵横波速度比具体为:依据干岩石纵横波速度比对应的待计算值和目标公式,确定目标层段对应的干岩石纵横波速度比,其中,目标公式为:
Figure BDA0001835609700000111
为了确定整个目标地区对应的干岩石纵横波速度比时,可以通过在该目标地区选定多个目标测量井,并获取该多个目标测量井分别对应的干岩石纵横波速度比,进而通过求平均确定整个目标地区对应的干岩石纵横波速度比。
举例说明,在目标地区包含至少一个目标测量井的情况下,依据每个目标测量井对应的多个采样点分别对应的目标测井数据,确定目标地区对应的干岩石纵横波速度比包括:获取每个目标测量井对应多个采样点分别对应的目标测井数据与目标地区对应的干岩石纵横波速度比之间的关系公式:
Figure BDA0001835609700000112
其中,m为目标测量井的数目,n为目标层段中的采样点的数量,
Figure BDA0001835609700000113
为目标层段中第j采样点的密度与纵波速度平方的乘积,
Figure BDA0001835609700000114
为目标层段中第j采样点的密度与横波速度平方的乘积,
Figure BDA0001835609700000115
为包含的多个采样点对应的密度与纵波速度平方的乘积的平均值,
Figure BDA0001835609700000116
为目标层段包含的多个采样点对应的密度与横波速度平方的乘积的平均值,
Figure BDA0001835609700000117
为目标地区对应的干岩石纵横波速度比;依据关系公式和多个采样点分别对应的目标测井数据,确定目标地区对应的干岩石纵横波速度比,其中,每个采样点对应的目标测井数据包括:采样点对应的纵波速度值、采样点对应的横波速度值、采样点对应的密度数值。
最后,在知晓目标测量井中目标层段的干岩石纵横波速度比的情况下,要想对目标测量井中目标层段的储层流体进行预测,则可以依据目标层段对应的干岩石纵横波速度比和目标层段中每个采样点对应的目标测井数据,确定目标层段对应的流体属性因子,以及获取目标层段所处的目标测量井对应的测井解释,进而依据流体属性因子和测井解释建立目标测量井对应的识别模型,其中,识别模型用于确定目标测量井中目标层段的储层流体的位置。
同理,在知晓目标地区对应的干岩石横纵波速度比的情况下,想要对目标地区的储层流体进行预测,则可以依据目标地区对应的干岩石纵横波速度比和多个目标测量井中每个采样点对应的目标测井数据,确定多个目标测量井对应的流体属性因子,以及获取目标地区包含的每个目标测量井对应的测井解释,进而依据多个目标测量井对应的流体属性因子和多个目标测量井对应的测井解释建立目标地区对应的识别模型,其中,识别模型用于确定目标地区的储层流体的位置。
也即,利用干岩石纵横波速比和纵横波速度、密度计算出流体属性因子f。对计算的流体属性因子f剖面和过该剖面井的测井解释结果进行标定处理,即依据测井解释结果,把油、气、水层与流体属性因子进行对应标定,此时,可以依据油、气、水层在该流体属性因子剖面上的显示特征,可对目标地区的储层进行流体检测,从而确定有利钻井区。
综上,本申请通过依据流体属性因子与剪切模量之间的岩石物理关系,推导出基于常规测井数据(纵横波速度和密度)即可计算出干岩石纵横波速度比的关系公式,达到了避免实验室测量干岩石纵横波速度比带来的高经济成本和测量误差,以及避免了采用岩石纵横波速度比经验值时导致的客观误差。此外,在依据本申请提供的确定方法,确定出干岩石纵横波速度比之后,基于该干岩石纵横波速度比计算出的流体属性因子识别孔隙流体效果最佳,进而达到了提高储层流体预测的效率和准确性,为油气检测提供了一种可靠方法的技术效果。
最后,针对本申请提供的干岩石纵横波速度比的确定方法进行测验验证说明,具体的,获取某预设地区的Marmousi模型示意图(图2),基于该某预设地区的Marmousi模型,分别确定该某预设地区对应的多个剖面示意图。
图3为某预设地区对应的纵横波速度的剖面示意图,该示意图较难区分油层和水层。
图4为某预设地区对应的密度的剖面示意图,该示意图中的密度变化区间较小,导致该示意图的真实性易受到密度数据的误差影响。
图5为某预设地区对应的横波速度的剖面示意图,该示意图只能识别砂泥岩。
图6为某预设地区对应的剪切模量的剖面示意图,该示意图也只能识别砂泥岩。
图7为某预设地区对应的横纵波速度比的剖面示意图,该示意图能对油、气、水进行较好的区分,但是仍然难以区分泥岩和油层。
图8为某预设地区对应的流体属性因子的剖面示意图,该示意图的水藏描述存在不清楚的问题,同时流体属性因子还包含了岩性信息,其中,该流体属性因子由干岩石纵横波速度比的经验值
Figure BDA0001835609700000131
计算所得。
图9为某预设地区对应的流体属性因子的剖面示意图,该示意图不仅仅能只识别油、气、水藏,同时还不再包含岩性信息,达到了对储层中的流体清楚识别的技术效果,其中,该流体属性因子由干岩石纵横波速度比的计算值
Figure BDA0001835609700000132
计算所得,该干岩石纵横波速度比的计算值由本申请实施例提供的干岩石纵横波速度比的确定方法确定。
显而易见,基于本申请实施例提供的干岩石纵横波速度比的确定方法得到的流体属性因子可以有效的指示不同类型的流体,解决了流体属性因子的不易精准且便捷获取的技术问题,提高了储层流体预测的效率和准确性。
需要说明的是,在附图的流程图示出的步骤可以在诸如一组计算机可执行指令的计算机系统中执行,并且,虽然在流程图中示出了逻辑顺序,但是在某些情况下,可以以不同于此处的顺序执行所示出或描述的步骤。
本申请实施例还提供了一种干岩石纵横波速度比的确定装置,需要说明的是,本申请实施例的干岩石纵横波速度比的确定装置可以用于执行本申请实施例所提供的用于干岩石纵横波速度比的确定方法。以下对本申请实施例提供的干岩石纵横波速度比的确定装置进行介绍。
图10是根据本申请实施例的干岩石纵横波速度比的确定装置的示意图。如图10所示,该装置包括:第一确定单元11、获取单元13和第二确定单元15。
第一确定单元11,用于确定与待确定的干岩石纵横波速度比相对应的目标层段,其中,目标层段位于目标测量井中。
获取单元13,用于获取目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合,其中,目标测井数据经由声波测井和密度测井获取。
第二确定单元15,用于依据目标测井数据集合中的目标测井数据,确定目标层段对应的干岩石纵横波速度比。
可选地,在本申请实施例提供的干岩石纵横波速度比的确定装置中,获取单元13包括:第一获取模块,用于获取目标测量井中目标层段对应的测井总数据,其中,测井总数据由目标测量井进行声波测井处理和密度测井处理而获取的;判断模块,用于判断目标层段对应的测井总数据中是否存在失真数据;第一确定模块,用于在目标层段对应的测井总数据中不存在失真数据的情况下,则依据目标测量井对应的测井总数据,确定目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合;重新确定模块,用于在目标层段对应的测井总数据中存在失真数据的情况下,则对测井总数据中的未失真数据进行多元线性拟合,确定未失真数据中多种数据之间的数学关系;以依据未失真数据中多种数据之间的数学关系,重新确定失真层段对应的测井计算数据,依据测井计算数据更替测井总数据中的失真数据,进而确定被更替失真数据后的测井总数据中不存在失真数据,其中,失真数据为失真层段对应的测井数据。
可选地,在本申请实施例提供的干岩石纵横波速度比的确定装置中,判断模块包括:获取子模块,用于在目标层段对应的测井总数据包含:横波时差曲线、纵波时差曲线、以及密度曲线的情况下,获取目标层段对应的测井曲线,其中,测井曲线包含以下至少之一:自然电位曲线、井径曲线、自然伽马曲线、深电阻率曲线、中电阻率曲线、浅电阻率曲线、声波曲线、密度曲线、中子曲线;确定子模块,用于依据目标层段对应的测井曲线对目标层段对应的测井总数据进行分析,确定中目标层段对应的测井总数据是否存在失真数据。
可选地,在本申请实施例提供的干岩石纵横波速度比的确定装置中,第二确定单元15包括:第二确定模块,用于在每个采样点对应的目标测井数据包括:采样点对应的纵波速度值、采样点对应的横波速度值、采样点对应的密度数值的情况下,依据目标测井数据集合中的目标测井数据,确定干岩石纵横波速度比对应的待计算值,其中,待计算值中包含:每个采样点对应的密度与纵波速度平方的乘积、每个采样点对应的密度与横波速度平方的乘积、目标层段包含的多个采样点对应的密度与纵波速度平方的乘积的平均值,目标层段包含的多个采样点对应的密度与横波速度平方的乘积的平均值、采样点的数量;第三确定模块,用于依据干岩石纵横波速度比对应的待计算值确定目标层段对应的干岩石纵横波速度比。
可选地,在本申请实施例提供的干岩石纵横波速度比的确定装置中,装置还包括:第三确定单元,用于在目标地区包含至少一个目标测量井的情况下,依据每个目标测量井对应的多个采样点分别对应的目标测井数据,确定目标地区对应的干岩石纵横波速度比,其中,第三确定单元包括:第二获取模块,用于获取每个目标测量井对应多个采样点分别对应的目标测井数据与目标地区对应的干岩石纵横波速度比之间的关系公式:
Figure BDA0001835609700000141
其中,m为目标测量井的数目,n为目标层段中的采样点的数量,
Figure BDA0001835609700000142
为目标层段中第j采样点的密度与纵波速度平方的乘积,
Figure BDA0001835609700000151
为目标层段中第j采样点的密度与横波速度平方的乘积,
Figure BDA0001835609700000152
为包含的多个采样点对应的密度与纵波速度平方的乘积的平均值,
Figure BDA0001835609700000153
为目标层段包含的多个采样点对应的密度与横波速度平方的乘积的平均值,
Figure BDA0001835609700000154
为目标地区对应的干岩石纵横波速度比;第四确定模块,用于依据关系公式和多个采样点分别对应的目标测井数据,确定目标地区对应的干岩石纵横波速度比。
可选地,在本申请实施例提供的干岩石纵横波速度比的确定装置中,方法还包括:第五确定单元,用于在目标地区对应的干岩石纵横波速度比之后,依据目标地区对应的干岩石纵横波速度比和多个目标测量井中每个采样点对应的目标测井数据,确定多个目标测量井对应的流体属性因子;第二建立单元,用于获取目标地区包含的每个目标测量井对应的测井解释,并依据多个目标测量井对应的流体属性因子和多个目标测量井对应的测井解释建立目标地区对应的识别模型,其中,识别模型用于确定目标地区的储层流体的位置。
本申请实施例提供的干岩石纵横波速度比的确定装置,通过第一确定单元11确定与待确定的干岩石纵横波速度比相对应的目标层段,其中,目标层段位于目标测量井中;获取单元13获取目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合,其中,目标测井数据经由声波测井和密度测井获取;第二确定单元15依据目标测井数据集合中的目标测井数据,确定目标层段对应的干岩石纵横波速度比,解决了相关技术中无法便捷、准确地确定干岩石纵横波速度比的问题。
需要说明的是:在常规测井工程中往往包含着对目标测量井进行声波测井处理和密度测井处理,基于此,本申请实施例通过获取目标测量井中目标层段对应的声波测井数据和密度测井数据,并基于该声波测井数据和密度测井数据,即可确定出目标测量井中目标层段对应的干岩石横纵波速度比,达到了快速且便捷的确定干岩石横纵波速度比的技术效果。
此外,相对于现有技术中的采用经验值确定干岩石横纵波速度比的技术方案,本申请实施例通过依据目标测量井中目标层段对应的具体的目标测量数据集合(目标层段中多个采样点分别对应的声波测井数据和密度测井数据),来确定该目标层段对应的干岩石横纵波速度比,达到了准确的确定干岩石横纵波速度比的技术效果。
干岩石纵横波速度比的确定装置包括处理器和存储器,上述第一确定单元11、获取单元13和第二确定单元15等均作为程序单元存储在存储器中,由处理器执行存储在存储器中的上述程序单元来实现相应的功能。
处理器中包含内核,由内核去存储器中调取相应的程序单元。内核可以设置一个或以上,通过调整内核参数来快速、便捷且精准的确定干岩石横纵波速度比。
存储器可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM),存储器包括至少一个存储芯片。
本发明实施例提供了一种存储介质,其上存储有程序,该程序被处理器执行时实现干岩石纵横波速度比的确定方法。
本发明实施例提供了一种处理器,处理器用于运行程序,其中,程序运行时执行干岩石纵横波速度比的确定方法。
本发明实施例提供了一种设备,设备包括处理器、存储器及存储在存储器上并可在处理器上运行的程序,处理器执行程序时实现以下步骤:确定与待确定的干岩石纵横波速度比相对应的目标层段,其中,目标层段位于目标测量井中;获取目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合,其中,目标测井数据经由声波测井和密度测井获取;依据目标测井数据集合中的目标测井数据,确定目标层段对应的干岩石纵横波速度比。
可选的,获取目标层段中多个采样点分别对应的目标测井数据包括:获取目标测量井中目标层段对应的测井总数据,其中,测井总数据由目标测量井进行声波测井处理和密度测井处理而获取的;判断目标层段对应的测井总数据中是否存在失真数据;在目标层段对应的测井总数据中不存在失真数据的情况下,则依据目标测量井对应的测井总数据,确定目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合;在目标层段对应的测井总数据中存在失真数据的情况下,则对测井总数据中的未失真数据进行多元线性拟合,确定未失真数据中多种数据之间的数学关系;以依据未失真数据中多种数据之间的数学关系,重新确定失真层段对应的测井计算数据,依据测井计算数据更替测井总数据中的失真数据,进而确定被更替失真数据后的测井总数据中不存在失真数据,其中,失真数据为失真层段对应的测井数据。
可选的,在目标层段对应的测井总数据包含:横波时差曲线、纵波时差曲线、以及密度曲线的情况下,判断目标层段对应的测井总数据中是否存在失真数据包括:获取目标层段对应的测井曲线,其中,测井曲线包含以下至少之一:自然电位曲线、井径曲线、自然伽马曲线、深电阻率曲线、中电阻率曲线、浅电阻率曲线、声波曲线、密度曲线、中子曲线;依据目标层段对应的测井曲线对目标层段对应的测井总数据进行分析,确定中目标层段对应的测井总数据是否存在失真数据。
可选的,在每个采样点对应的目标测井数据包括:采样点对应的纵波速度值、采样点对应的横波速度值、采样点对应的密度数值的情况下,依据目标测井数据集合中的目标测井数据,确定目标层段对应的干岩石纵横波速度比包括:依据目标测井数据集合中的目标测井数据,确定干岩石纵横波速度比对应的待计算值,其中,待计算值中包含:每个采样点对应的密度与纵波速度平方的乘积、每个采样点对应的密度与横波速度平方的乘积、目标层段包含的多个采样点对应的密度与纵波速度平方的乘积的平均值,目标层段包含的多个采样点对应的密度与横波速度平方的乘积的平均值、采样点的数量;依据干岩石纵横波速度比对应的待计算值确定目标层段对应的干岩石纵横波速度比。
可选的,在目标地区包含至少一个目标测量井的情况下,依据每个目标测量井对应的多个采样点分别对应的目标测井数据,确定目标地区对应的干岩石纵横波速度比包括:获取每个目标测量井对应多个采样点分别对应的目标测井数据与目标地区对应的干岩石纵横波速度比之间的关系公式:
Figure BDA0001835609700000171
其中,m为目标测量井的数目,n为目标层段中的采样点的数量,
Figure BDA0001835609700000172
为目标层段中第j采样点的密度与纵波速度平方的乘积,
Figure BDA0001835609700000173
为目标层段中第j采样点的密度与横波速度平方的乘积,
Figure BDA0001835609700000174
为包含的多个采样点对应的密度与纵波速度平方的乘积的平均值,
Figure BDA0001835609700000175
为目标层段包含的多个采样点对应的密度与横波速度平方的乘积的平均值,
Figure BDA0001835609700000176
为目标地区对应的干岩石纵横波速度比;依据关系公式和多个采样点分别对应的目标测井数据,确定目标地区对应的干岩石纵横波速度比。
可选的,在目标地区对应的干岩石纵横波速度比之后,方法还包括:依据目标地区对应的干岩石纵横波速度比和多个目标测量井中每个采样点对应的目标测井数据,确定多个目标测量井对应的流体属性因子;获取目标地区包含的每个目标测量井对应的测井解释,并依据多个目标测量井对应的流体属性因子和多个目标测量井对应的测井解释建立目标地区对应的识别模型,其中,识别模型用于确定目标地区的储层流体的位置。本文中的设备可以是服务器、PC、PAD、手机等。
本申请还提供了一种计算机程序产品,当在数据处理设备上执行时,适于执行初始化有如下方法步骤的程序:确定与待确定的干岩石纵横波速度比相对应的目标层段,其中,目标层段位于目标测量井中;获取目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合,其中,目标测井数据经由声波测井和密度测井获取;依据目标测井数据集合中的目标测井数据,确定目标层段对应的干岩石纵横波速度比。
可选的,获取目标层段中多个采样点分别对应的目标测井数据包括:获取目标测量井中目标层段对应的测井总数据,其中,测井总数据由目标测量井进行声波测井处理和密度测井处理而获取的;判断目标层段对应的测井总数据中是否存在失真数据;在目标层段对应的测井总数据中不存在失真数据的情况下,则依据目标测量井对应的测井总数据,确定目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合;在目标层段对应的测井总数据中存在失真数据的情况下,则对测井总数据中的未失真数据进行多元线性拟合,确定未失真数据中多种数据之间的数学关系;以依据未失真数据中多种数据之间的数学关系,重新确定失真层段对应的测井计算数据,依据测井计算数据更替测井总数据中的失真数据,进而确定被更替失真数据后的测井总数据中不存在失真数据,其中,失真数据为失真层段对应的测井数据。
可选的,在目标层段对应的测井总数据包含:横波时差曲线、纵波时差曲线、以及密度曲线的情况下,判断目标层段对应的测井总数据中是否存在失真数据包括:获取目标层段对应的测井曲线,其中,测井曲线包含以下至少之一:自然电位曲线、井径曲线、自然伽马曲线、深电阻率曲线、中电阻率曲线、浅电阻率曲线、声波曲线、密度曲线、中子曲线;依据目标层段对应的测井曲线对目标层段对应的测井总数据进行分析,确定中目标层段对应的测井总数据是否存在失真数据。
可选的,在每个采样点对应的目标测井数据包括:采样点对应的纵波速度值、采样点对应的横波速度值、采样点对应的密度数值的情况下,依据目标测井数据集合中的目标测井数据,确定目标层段对应的干岩石纵横波速度比包括:依据目标测井数据集合中的目标测井数据,确定干岩石纵横波速度比对应的待计算值,其中,待计算值中包含:每个采样点对应的密度与纵波速度平方的乘积、每个采样点对应的密度与横波速度平方的乘积、目标层段包含的多个采样点对应的密度与纵波速度平方的乘积的平均值,目标层段包含的多个采样点对应的密度与横波速度平方的乘积的平均值、采样点的数量;依据干岩石纵横波速度比对应的待计算值确定目标层段对应的干岩石纵横波速度比。
可选的,在目标地区包含至少一个目标测量井的情况下,依据每个目标测量井对应的多个采样点分别对应的目标测井数据,确定目标地区对应的干岩石纵横波速度比包括:获取每个目标测量井对应多个采样点分别对应的目标测井数据与目标地区对应的干岩石纵横波速度比之间的关系公式:
Figure BDA0001835609700000191
其中,m为目标测量井的数目,n为目标层段中的采样点的数量,
Figure BDA0001835609700000192
为目标层段中第j采样点的密度与纵波速度平方的乘积,
Figure BDA0001835609700000193
为目标层段中第j采样点的密度与横波速度平方的乘积,
Figure BDA0001835609700000194
为包含的多个采样点对应的密度与纵波速度平方的乘积的平均值,
Figure BDA0001835609700000195
为目标层段包含的多个采样点对应的密度与横波速度平方的乘积的平均值,
Figure BDA0001835609700000196
为目标地区对应的干岩石纵横波速度比;依据关系公式和多个采样点分别对应的目标测井数据,确定目标地区对应的干岩石纵横波速度比。
可选的,在目标地区对应的干岩石纵横波速度比之后,方法还包括:依据目标地区对应的干岩石纵横波速度比和多个目标测量井中每个采样点对应的目标测井数据,确定多个目标测量井对应的流体属性因子;获取目标地区包含的每个目标测量井对应的测井解释,并依据多个目标测量井对应的流体属性因子和多个目标测量井对应的测井解释建立目标地区对应的识别模型,其中,识别模型用于确定目标地区的储层流体的位置。
本领域内的技术人员应明白,本申请的实施例可提供为方法、系统、或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器,使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中,使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品,该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上,使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理,从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中,计算设备包括一个或多个处理器(CPU)、输入/输出接口、网络接口和内存。
存储器可能包括计算机可读介质中的非永久性存储器,随机存取存储器(RAM)和/或非易失性内存等形式,如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括,但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器(CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带,磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质,可用于存储可以被计算设备访问的信息。按照本文中的界定,计算机可读介质不包括暂存电脑可读媒体(transitory media),如调制的数据信号和载波。
还需要说明的是,术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下,由语句“包括一个……”限定的要素,并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
本领域技术人员应明白,本申请的实施例可提供为方法、系统或计算机程序产品。因此,本申请可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且,本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
以上仅为本申请的实施例而已,并不用于限制本申请。对于本领域技术人员来说,本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等,均应包含在本申请的权利要求范围之内。

Claims (10)

1.一种干岩石纵横波速度比的确定方法,其特征在于,包括:
确定与待确定的干岩石纵横波速度比相对应的目标层段,其中,所述目标层段位于目标测量井中;
获取所述目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合,其中,所述目标测井数据经由声波测井和密度测井获取;
依据所述目标测井数据集合中的目标测井数据,确定所述目标层段对应的干岩石纵横波速度比。
2.根据权利要求1所述的方法,其特征在于,获取所述目标层段中多个采样点分别对应的目标测井数据包括:
获取所述目标测量井中目标层段对应的测井总数据,其中,所述测井总数据由所述目标测量井进行声波测井处理和密度测井处理而获取的;
判断所述目标层段对应的测井总数据中是否存在失真数据;
在所述目标层段对应的测井总数据中不存在失真数据的情况下,则依据所述目标测量井对应的测井总数据,确定所述目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合;
在所述目标层段对应的测井总数据中存在失真数据的情况下,则对所述测井总数据中的未失真数据进行多元线性拟合,确定所述未失真数据中多种数据之间的数学关系;以依据所述未失真数据中多种数据之间的数学关系,重新确定失真层段对应的测井计算数据,依据所述测井计算数据更替所述测井总数据中的失真数据,并确定被更替所述失真数据后的测井总数据中不存在失真数据,其中,失真数据为所述失真层段对应的测井数据。
3.根据权利要求2所述的方法,其特征在于,在所述目标层段对应的测井总数据包含:横波时差曲线、纵波时差曲线、以及密度曲线的情况下,判断所述目标层段对应的测井总数据中是否存在失真数据包括:
获取所述目标层段对应的测井曲线,其中,所述测井曲线包含以下至少之一:自然电位曲线、井径曲线、自然伽马曲线、深电阻率曲线、中电阻率曲线、浅电阻率曲线、声波曲线、密度曲线、中子曲线;
依据所述目标层段对应的测井曲线对所述目标层段对应的测井总数据进行分析,确定中所述目标层段对应的测井总数据是否存在失真数据。
4.根据权利要求1所述的方法,其特征在于,在每个所述采样点对应的目标测井数据包括:所述采样点对应的纵波速度值、所述采样点对应的横波速度值、所述采样点对应的密度数值的情况下,依据所述目标测井数据集合中的目标测井数据,确定所述目标层段对应的干岩石纵横波速度比包括:
依据所述目标测井数据集合中的目标测井数据,确定所述干岩石纵横波速度比对应的待计算值,其中,所述待计算值中包含:每个所述采样点对应的密度与纵波速度平方的乘积、每个所述采样点对应的密度与横波速度平方的乘积、所述目标层段包含的多个采样点对应的密度与纵波速度平方的乘积的平均值,所述目标层段包含的多个采样点对应的密度与横波速度平方的乘积的平均值、所述采样点的数量;
依据所述干岩石纵横波速度比对应的待计算值确定所述目标层段对应的干岩石纵横波速度比。
5.根据权利要求1所述的方法,其特征在于,在目标地区包含至少一个所述目标测量井的情况下,依据每个所述目标测量井对应的多个采样点分别对应的目标测井数据,确定所述目标地区对应的干岩石纵横波速度比包括:
获取每个所述目标测量井对应多个采样点分别对应的目标测井数据与所述目标地区对应的干岩石纵横波速度比之间的关系公式:
Figure FDA0001835609690000021
其中,m为目标测量井的数目,n为目标层段中的采样点的数量,
Figure FDA0001835609690000022
为目标层段中第j采样点的密度与纵波速度平方的乘积,
Figure FDA0001835609690000023
为目标层段中第j采样点的密度与横波速度平方的乘积,
Figure FDA0001835609690000024
为包含的多个采样点对应的密度与纵波速度平方的乘积的平均值,
Figure FDA0001835609690000025
为所述目标层段包含的多个采样点对应的密度与横波速度平方的乘积的平均值,
Figure FDA0001835609690000026
为所述目标地区对应的干岩石纵横波速度比;
依据所述关系公式和所述多个采样点分别对应的目标测井数据,确定所述目标地区对应的干岩石纵横波速度比。
6.根据权利要求5所述的方法,其特征在于,在所述目标地区对应的干岩石纵横波速度比之后,所述方法还包括:
依据所述目标地区对应的干岩石纵横波速度比和多个所述目标测量井中每个所述采样点对应的目标测井数据,确定多个所述目标测量井对应的流体属性因子;
获取所述目标地区包含的每个目标测量井对应的测井解释,并依据多个所述目标测量井对应的流体属性因子和多个所述目标测量井对应的测井解释建立所述目标地区对应的识别模型,其中,所述识别模型用于确定所述目标地区的储层流体的位置。
7.一种干岩石纵横波速度比的确定装置,其特征在于,包括:
第一确定单元,用于确定与待确定的干岩石纵横波速度比相对应的目标层段,其中,所述目标层段位于目标测量井中;
获取单元,用于获取所述目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合,其中,所述目标测井数据经由声波测井和密度测井获取;
第二确定单元,用于依据所述目标测井数据集合中的目标测井数据,确定所述目标层段对应的干岩石纵横波速度比。
8.根据权利要求7所述的装置,其特征在于,所述获取单元包括:
第一获取模块,用于获取所述目标测量井中目标层段对应的测井总数据,其中,所述测井总数据由所述目标测量井进行声波测井处理和密度测井处理而获取的;
判断模块,用于判断所述目标层段对应的测井总数据中是否存在失真数据;
第一确定模块,用于在所述目标层段对应的测井总数据中不存在失真数据的情况下,则依据所述目标测量井对应的测井总数据,确定所述目标层段中多个采样点分别对应的目标测井数据,得到目标测井数据集合;
重新确定模块,用于在所述目标层段对应的测井总数据中存在失真数据的情况下,则对所述测井总数据中的未失真数据进行多元线性拟合,确定所述未失真数据中多种数据之间的数学关系;以依据所述未失真数据中多种数据之间的数学关系,重新确定失真层段对应的测井计算数据,依据所述测井计算数据更替所述测井总数据中的失真数据,进而确定被更替所述失真数据后的测井总数据中不存在失真数据,其中,失真数据为所述失真层段对应的测井数据。
9.一种存储介质,其特征在于,所述存储介质包括存储的程序,其中,所述程序执行权利要求1至6中任意一项所述的干岩石纵横波速度比的确定方法。
10.一种处理器,其特征在于,所述处理器用于运行程序,其中,所述程序运行时执行权利要求1至6中任意一项所述的干岩石纵横波速度比的确定方法。
CN201811224519.2A 2018-10-19 2018-10-19 干岩石纵横波速度比的确定方法和装置 Active CN111077589B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811224519.2A CN111077589B (zh) 2018-10-19 2018-10-19 干岩石纵横波速度比的确定方法和装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811224519.2A CN111077589B (zh) 2018-10-19 2018-10-19 干岩石纵横波速度比的确定方法和装置

Publications (2)

Publication Number Publication Date
CN111077589A true CN111077589A (zh) 2020-04-28
CN111077589B CN111077589B (zh) 2022-06-03

Family

ID=70309565

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811224519.2A Active CN111077589B (zh) 2018-10-19 2018-10-19 干岩石纵横波速度比的确定方法和装置

Country Status (1)

Country Link
CN (1) CN111077589B (zh)

Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2463631C1 (ru) * 2011-04-08 2012-10-10 Федеральное государственное бюджетное научное учреждение "Научно-исследовательский институт аэрокосмического мониторинга "АЭРОКОСМОС" Способ обнаружения очагов землетрясений сетью сейсмостанций
CN103630939A (zh) * 2013-11-06 2014-03-12 中国石油天然气集团公司 一种气层识别评价方法
CN103675911A (zh) * 2014-01-06 2014-03-26 中国石油大学(华东) 一种基于纵波和转换横波联合反演截距和梯度的方法
CN104181585A (zh) * 2014-06-27 2014-12-03 中国石油化工股份有限公司 一种地球物理勘探横波估算方法和系统
CN104181610A (zh) * 2014-08-07 2014-12-03 中国石油天然气集团公司 一种射线路径弹性反演方法以及系统
CN104950331A (zh) * 2015-06-30 2015-09-30 中国石油天然气集团公司 一种砂泥岩储层的孔隙度与泥质含量的地震预测方法
CN105988136A (zh) * 2015-01-30 2016-10-05 中石化石油工程技术服务有限公司 一种利用纵横波速度信息进行含气性分析的方法
CN107367762A (zh) * 2017-06-27 2017-11-21 中国石油天然气股份有限公司 一种确定储层参数的方法及装置
CN107422384A (zh) * 2017-03-20 2017-12-01 中石化石油工程技术服务有限公司 一种利用测井资料仿垂直地震剖面资料的方法

Patent Citations (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2463631C1 (ru) * 2011-04-08 2012-10-10 Федеральное государственное бюджетное научное учреждение "Научно-исследовательский институт аэрокосмического мониторинга "АЭРОКОСМОС" Способ обнаружения очагов землетрясений сетью сейсмостанций
CN103630939A (zh) * 2013-11-06 2014-03-12 中国石油天然气集团公司 一种气层识别评价方法
CN103675911A (zh) * 2014-01-06 2014-03-26 中国石油大学(华东) 一种基于纵波和转换横波联合反演截距和梯度的方法
CN104181585A (zh) * 2014-06-27 2014-12-03 中国石油化工股份有限公司 一种地球物理勘探横波估算方法和系统
CN104181610A (zh) * 2014-08-07 2014-12-03 中国石油天然气集团公司 一种射线路径弹性反演方法以及系统
CN105988136A (zh) * 2015-01-30 2016-10-05 中石化石油工程技术服务有限公司 一种利用纵横波速度信息进行含气性分析的方法
CN104950331A (zh) * 2015-06-30 2015-09-30 中国石油天然气集团公司 一种砂泥岩储层的孔隙度与泥质含量的地震预测方法
CN107422384A (zh) * 2017-03-20 2017-12-01 中石化石油工程技术服务有限公司 一种利用测井资料仿垂直地震剖面资料的方法
CN107367762A (zh) * 2017-06-27 2017-11-21 中国石油天然气股份有限公司 一种确定储层参数的方法及装置

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
杜启振: "基于横纵波速度比值迭代的纵横波联合反演方法", 《石油地球物理勘探》 *
霍国栋 等: "纵向和横向同时约束AVO反演", 《地球物理学报》 *

Also Published As

Publication number Publication date
CN111077589B (zh) 2022-06-03

Similar Documents

Publication Publication Date Title
CA3070868C (en) Resource density screening tool
US11105943B2 (en) System and method for mapping hydrocarbon source rock using seismic attributes
US10007015B2 (en) Methods, systems and devices for predicting reservoir properties
US20200095858A1 (en) Modeling reservoir permeability through estimating natural fracture distribution and properties
US20160161635A1 (en) A simulation-to-seismic workflow construed from core based rock typing and enhanced by rock replacement modeling
Teixeira et al. Quantitative seismic interpretation integrated with well-test analysis in turbidite and presalt reservoirs
Germán Rubino et al. Seismic characterization of thin beds containing patchy carbon dioxide-brine distributions: A study based on numerical simulations
CN108873065B (zh) 砂岩优质储层预测方法及装置
Zaree et al. Estimating fracture intensity in hydrocarbon reservoir: an approach using DSI data analysis
CN111381292B (zh) 一种预测砂岩含烃储层的测井解释方法与装置
RU2630852C1 (ru) Способ прогноза эффективной емкости коллекторов на основе получаемых поляризационных параметров и проводимости для выбранного типа среды
CN111812716B (zh) 页岩气储层总有机碳含量叠前定量预测方法、装置及设备
CN112185469B (zh) 一种预测海域天然气水合物有利聚集区的方法
Orellana et al. Influence of variograms in 3D reservoir-modeling outcomes: An example
US11346968B2 (en) Estimation of reservoir flow properties from seismic data
CN111077589B (zh) 干岩石纵横波速度比的确定方法和装置
Guo et al. Frequency-dependent P-wave anisotropy due to scattering in rocks with aligned fractures
Sauermilch et al. A comparative analysis of time–depth relationships derived from scientific ocean drilling expeditions
Gutierrez et al. Rock physics workflows for exploration in frontier basins
US20190129048A1 (en) System and method for assessing the presence of hydrocarbons in a subterranean reservoir based on seismic inversions
CN113052356B (zh) 油井单井产能预测方法、装置、电子设备和存储介质
Yemets et al. Evaluation of the Void Space Structure and Reservoir Properties of Rocks in the Visevian and Turnaisian Formations of the Berezivske Field of the Dnipro-Donetsk Depression
CN111856571B (zh) 测井频段速度向地震频段速度转换方法、装置及设备
RU2566424C2 (ru) Способ и технологический комплекс для анализа нелинейных свойств среды с целью расширения спектра регистрируемого волнового сигнала
CN111751900B (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