CN107807409B - 地层密度和电阻率关系的确定方法和装置 - Google Patents
地层密度和电阻率关系的确定方法和装置 Download PDFInfo
- Publication number
- CN107807409B CN107807409B CN201710811319.6A CN201710811319A CN107807409B CN 107807409 B CN107807409 B CN 107807409B CN 201710811319 A CN201710811319 A CN 201710811319A CN 107807409 B CN107807409 B CN 107807409B
- Authority
- CN
- China
- Prior art keywords
- density
- data
- relation
- resistivity
- stratum
- 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
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Abstract
本申请实施方式提供了一种地层密度和电阻率关系的确定方法和装置,该方法包括:获取目标区域的测井数据和地质分层数据;按照第一间隔,从测井数据中,获取第一采样点的电阻率数据和密度数据;根据地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到地层的第一数据体;构建关于密度和电阻率关系的多个关系函数;从多个关系函数中,确定出符合预设要求的关系函数,作为地层密度和电阻率关系。由于该方案通过综合利用测井数据和地质分层数据,构建多个关系函数,进而从多个关系函数中确定出较为准确的关系函数作为地层密度和电阻率关系,因而解决了现有方法中存在的确定地层密度和电阻率关系误差较大的技术问题。
Description
技术领域
本申请涉及地球物理勘探技术领域,特别涉及一种地层密度和电阻率关系的确定方法和装置。
背景技术
在进行地球物理勘探的过程,常常需要综合利用重力和电磁方法进行联合反演,进而可以根据反演的结果,对目标区域进行具体的地质解释和钻探。为了能较准确地进行重力和电磁的联合反演,通常需要先确定密度和电阻率关系,以上述关系作为物性基础,以便进行较为准确的联合反演。
目前,为了确定地层密度和电阻率关系,通常是先求解密度和速度关系、速度和电阻率关系,再根据密度和速度关系、速度和电阻率关系,间接确定密度和电阻率关系。由于是通过间接的方式确定密度和电阻率关系,因此,上述方法具体实施时,往往会存在所确定的密度和电阻率关系的误差较大、准确度较低的技术问题。
针对上述问题,目前尚未提出有效的解决方案。
发明内容
本申请实施方式提供了一种地层密度和电阻率关系的确定方法和装置,以解决现有方法中存在的所确定的密度和电阻率关系误差较大、准确度较低的技术问题。
本申请实施方式提供了一种地层密度和电阻率关系的确定方法,包括:
获取目标区域的测井数据和地质分层数据;
按照第一间隔,从所述测井数据中获取多个第一采样点的电阻率数据和密度数据;
根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体;
根据各个地层的第一数据体,构建关于地层密度和电阻率关系的多个关系函数;
从所述多个关系函数中,确定出符合预设要求的关系函数,作为所述目标区域的地层密度和电阻率关系。
在一个实施方式中,所述第一间隔为5米。
在一个实施方式中,从所述测井数据中,获取多个第一采样点的密度数据,包括:
确定所述测井数据中是否包括密度测井数据;
在所述测井数据中包括密度测井数据的情况下,根据所述密度测井数据,确定所述多个第一采样点的密度数据;
在所述测井数据中不包括密度测井数据的情况下,根据声波测井数据,确定所述多个第一采样点的密度数据。
在一个实施方式中,根据所述声波测井数据,确定所述多个第一采样点的密度数据,包括:
按照以下公式确定所述多个第一采样点的密度数据:
den=a(b/AC)1/4
上式中,den为密度数据,AC为声波测井数据,a为第一系数,b为第二系数。
在一个实施方式中,根据各个地层的第一数据体,构建关于地层密度和电阻率关系的多个关系函数,包括:
根据各个地层的第一数据体,按照以下公式,构建关于地层密度和电阻率关系的多个关系函数,其中,所述多个关系函数包括以下至少之一:第一关系函数、第二关系函数、第三关系函数:
y2=a0+a1x
y3=a0+a1x+a2x2
上式中,y1为多个关系函数中第一关系函数的拟合密度,y2为多个关系函数中第二关系函数的拟合密度,y3为多个关系函数中第三关系函数的拟合密度,x为电阻率的对数,a0为第一拟合系数,a1为第二拟合系数,a3为第三拟合系数。
在一个实施方式中,从所述多个关系函数中,确定出符合预设要求的关系函数,包括:
根据所述多个关系函数,确定所述多个关系函数中各个关系函数的评价参数;
根据所述各个关系函数的评价参数,确定评价参数中数值最大的关系函数,作为所述符合预设要求的关系函数。
在一个实施方式中,根据所述多个关系函数,确定所述多个关系函数中各个关系函数的评价参数,包括:
按照以下公式,确定关系函数的评价参数:
上式中,r2为关系函数的评价参数,yobs为实测密度,y为关系函数对应的拟合密度,为关系函数对应的拟合密度的平均值。
在一个实施方式中,在获取多个第一采样点的电阻率数据和密度数据后,所述方法还包括:
对所述多个第一采样点的电阻数据和密度数据进行低通滤波处理,得到处理后的第一采样点的电阻率数据和密度数据;
相应的,根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体,包括:
根据所述地质分层数据,合并位于同一地层的处理后的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体。
在一个实施方式中,根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体后,所述方法还包括:
按照第二间隔,从各个地层的第一数据体中,抽取多个第二采样点的电阻率数据和密度数据作为各个地层的第二数据体;
相应的,根据各个地层的第一数据体,构建关于地层密度和电阻率关系的多个关系函数,包括:
根据各个地层的第二数据体,构建关于地层密度和电阻率关系的多个关系函数。
在一个实施方式中,在确定所述目标区域的地层密度和电阻率关系后,所述方法还包括:
根据所述各个地层的第一数据体,建立密度电阻率交汇图;
根据所述密度电阻率交汇图,对所述目标区域的地层密度和电阻率关系进行校验。
在一个实施方式中,在确定所述目标区域的地层密度和电阻率关系后,所述方法还包括:
根据所述目标区域的电阻率和密度关系,对目标区域进行电磁和重力的联合反演,得到反演结果;
根据所述反演结果,对目标区域进行地质解释和钻探。
本申请实施方式还提供了一种地层密度和电阻率关系的确定装置,包括:
第一获取模块,用于获取目标区域的测井数据和地质分层数据;
第二获取模块,用于按照第一间隔,从测井数据中获取多个第一采样点的电阻率数据和密度数据;
合并模块,用于根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体;
构建模块,用于根据各个地层的第一数据体,构建关于地层密度和电阻率关系的多个关系函数;
确定模块,用于从所述多个关系函数中,确定出符合预设要求的关系函数,作为所述目标区域的地层密度和电阻率关系。
在一个实施方式中,所述装置包括施工模块,用于根据所述目标区域的电阻率和密度关系,对目标区域进行电磁和重力的联合反演,得到反演结果;根据所述反演结果,对目标区域进行地质解释和钻探。
在本申请实施方式中,通过考虑了地层密度数据、电阻率数据的相应特点,综合利用测井数据和地质分层数据,构建多个关系函数,进而从多个关系函数中确定出较为准确的关系函数作为地层密度和电阻率关系,从而解决了现有方法中存在的确定的地层密度和电阻率关系误差较大、准确度较低的技术问题,达到准确、快速地确定地层密度和电阻率关系的技术效果。
附图说明
为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1是根据本申请实施方式提供的地层密度和电阻率关系的确定方法的处理流程图;
图2是根据本申请实施方式提供的地层密度和电阻率关系的确定装置的组成结构图;
图3是在一个场景示例中应用本申请实施方式提供地层密度和电阻率关系的确定方法/装置的流程示意图。
具体实施方式
为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护的范围。
考虑到现有方法,在确定地层密度和电阻率关系时,往往没有具体的考虑到密度数据、电阻率数据自身的特点,常常是先确定密度和速度关系、速度和电阻率关系,再根据密度和速度关系、速度和电阻率关系,间接确定密度和电阻率关系。导致具体实施时,往往存在所确定的地层密度和电阻率关系误差较大、准确度较低的技术问题。针对产生上述技术问题的根本原因,本申请考虑可以充分地考虑密度数据、电阻率数据的特点,综合利用测井数据和地质分层数据,构建多个关系函数,进而从多个关系函数中确定出较为准确的关系函数作为地层密度和电阻率关系,从而,解决现有方法中存在的所确定的地层密度和电阻率关系的误差较大、准确度较低的技术问题,达到准确、快速地确定地层密度和电阻率关系的技术效果。
其中,上述地层密度和电阻率关系具体可以是一种函数关系,用以表征目标区域中地层密度和电阻率的对应关系。由于重力法、电磁法分别属于不同类型的方法,在综合利用重力法和电磁法进行重力和电磁的联合反演时,往往需要先确定用于反演的物性参数具有共同的物性界面,或者不同物性参数之间具有关联关系。上述地层密度和电阻率关系可以用于表征地层中密度数据和电阻率数据的关联关系。进而,可以利用上述关系作为物性基础,进行较为准确的重力和电磁的联合反演。
基于上述思考思路,本申请实施方式提供了一种地层密度和电阻率关系的确定方法。请参阅图1所示的根据本申请实施方式提供的地层密度和电阻率关系的确定方法的处理流程图。本申请实施方式提供的地层密度和电阻率关系的确定方法,具体可以包括以下步骤。
S11:获取目标区域的测井数据和地质分层数据。
在本实施方式中,所述目标区域具体可以包括多个测井。相应的,上述测井数据具体可以是上述多个测井中各个测井的测井数据,上述地质分层数据具体可以是上述多个测井中各个测井的地质分层数据。
在一个实施方式中,上述测井数据具体可以包括电阻率测井数据和密度测井数据;也可以包括电阻率测井数据和声波时差测井数据。
在一个实施方式中,上述地质分层数据具体可以是包含有对应测井的各个地层层位信息的数据。具体的,例如可以包括地层层位深度信息、地层层位划分信息等等。具体实施时,为了获取测井的地质分层数据,可以根据对应测井的电阻率测井数据,确定测井中各个位置的电阻率值;以测井中各个位置的电阻率值作为依据,对测井进行地层层位的划分,从而得到上述地质分层数据。
S12:按照第一间隔,从所述测井数据中获取多个第一采样点的电阻率数据和密度数据。
在一个实施方式中,为了快速、准确地从大量测井数据中获取具有较好表征效果的第一采样点的电阻率数据和密度数据,可以按照第一间隔,从所述测井数据中获取多个第一采样点的电阻率数据和密度数据。具体实施时,可以按照以下内容,分别从各个测井中获取多个第一采样点的电阻率数据和密度数据:按照上述第一间隔,从各个测井中分别抽取多个第一采样点;根据对应测井的测井数据,获取上述多个采样点中各个采样点的电阻率数据和密度数据。
在一个实施方式中,上述第一间隔具体可以设置为5米。相应的,具体实施时,可以每隔5米从对应测井中抽取一个测点作为一个第一采样点,再根据该测井的测井数据,获取对应的第一采样点的电阻率数据和密度数据。在本实施方式中,需要说明的是,将第一间隔设置为5米,可以避免采集大量的采样点数据,减少数据处理量,提高数据处理效率;同时按照上述第一间隔抽取多个第一采样点的电阻率数据和密度数据,具有较好的代表性,可以较为清楚、稳定地反映出对应测井中地层密度和电阻率的变化规律。当然,需要说明的是,上述所列举的5米是为了更好地说明本申请实施方式,具体实施时,也可以根据具体情况和施工要求设置其他的数值作为上述第一间隔。
在一个实施方式中,上述从所述测井数据中,获取多个第一采样点的密度数据,具体可以包括以下内容:
S1:确定所述测井数据中是否包括密度测井数据;
S2:在所述测井数据中包括密度测井数据的情况下,根据所述密度测井数据,确定所述多个第一采样点的密度数据;
S3:在所述测井数据中不包括密度测井数据的情况下,根据声波测井数据,确定所述多个第一采样点的密度数据。
在一个实施方式中,在所述测井数据不包括密度测井数据、包括声波测井数据的情况下,为了获得较为准确的第一采样点的密度数据,具体实施时,可以按照以下公式,根据测井数据中的声波测井数据,分别确定上述多个第一采样点中各个第一采样点的密度数据:
den=a(b/AC)1/4
上式中,den可以为密度数据,AC可以为声波测井数据,a可以为第一系数,b可以为第二系数。其中,上述声波测井数据具体可以是一种声波时差测井数据。
在一个实施方式中,上述第一系数、第二系数具体可以按照以下方式取值:
a=0.31,b=0.3048×106
如此,可以根据声波时差测井数据较为准确地确定出各个第一采样点的密度数据。
在一个实施方式中,具体实施时,可以根据测井数据中的电阻率测井数据,分别确定上述各个第一采样点的电阻率数据。
S13:根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体。
在一个实施方式中,具体实施时,上述根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体,具体可以包括以下内容:
S13-1:根据地质分层数据中划分的各个地层层位的深度信息、各个第一采样点的深度信息,确定出位于同一地层中的第一采样点。其中,上述位于同一地层中的第一采样点包括了不同测井中位于同一地层的多个第一采样点,上述第一采样点的深度信息,具体可以根据第一采样点所在测井的测井数据确定。
S13-2:将位于同一地层的第一采样点的电阻数据和密度数据进行合并,得到对应地层的第一数据体,进而可以获得目标区域中各个地层的第一数据体。以备后续处理。
S14:根据各个地层的第一数据体,构建关于地层密度和电阻率关系的多个关系函数。
在本实施方式中,为了较为准确地确定出目标区域的地层密度和电阻率关系,具体实施时可以先根据各个地层的第一数据体,通过图线拟合,确定多个可以用于表征地层密度和电阻率关系的多个关系函数;再从多个关系函数中选出表征效果相对较好的关系函数作为目标区域的地层密度和电阻率关系。
在一个实施方式中,上述根据各个地层的第一数据体,构建关于地层密度和电阻率关系的多个关系函数,具体可以包括:
根据各个地层的第一数据体,按照以下公式,构建关于地层密度和电阻率关系的多个关系函数,其中,所述多个关系函数包括以下至少之一:第一关系函数、第二关系函数、第三关系函数;
y2=a0+a1x
y3=a0+a1x+a2x2
上式中,y1可以为多个关系函数中第一关系函数的拟合密度,y2可以为多个关系函数中第二关系函数的拟合密度,y3可以为多个关系函数中第三关系函数的拟合密度,x可以为电阻率的对数,a0可以为第一拟合系数,a1可以为第二拟合系数,a3可以为第三拟合系数。
在本实施方式中,需要说明的是,上述所列举的三种类型的关系函数只是为了更好地说明本申请实施方式,具体实施时,也可以根据具体情况和施工要求,引入其他类型的多项式作为关系函数。例如,也可以引入四次项多项式作为一个关系函数。
S15:从所述多个关系函数中,确定出符合预设要求的关系函数,作为所述目标区域的地层密度和电阻率关系。
在一个实施方式中,上述从所述多个关系函数中,确定出符合预设要求的关系函数,具体可以包括以下内容:
S15-1:根据所述多个关系函数,确定所述多个关系函数中各个关系函数的评价参数。
S15-2:根据所述各个关系函数的评价参数,确定评价参数中数值最大的关系函数,作为所述符合预设要求的关系函数。
在一个实施方式中,为了从多个关系函数中确定出较为准确的地层密度和电阻率关系,具体实施时,可以先分别确定多个关系函数中各个关系函数的评价参数。再从多个评价参数中确定数值最大的评价参数所对应的关系函数作为符合预设要求的关系函数。其中,上述符合预设要求的关系函数在多个关系函数中,表征地层密度和电阻率关系的准确度相对较高。将上述符合预设要求的关系函数,确定为所述目标区域的地层密度和电阻率关系。
在一个实施方式中,具体实施时,可以按照以下公式,分别确定多个关系函数中各个关系函数的评价参数:
上式中,r2可以为关系函数的评价参数,yobs可以为实测密度,y可以为关系函数对应的拟合密度,y可以为关系函数对应的拟合密度的平均值。
在本申请实施例中,相较于现有技术,本申请具体考虑了地层密度数据、电阻率数据的相应特点,通过综合利用测井数据和地质分层数据,构建多个关系函数,进而从多个关系函数中确定出较为准确的关系函数作为地层密度和电阻率关系,从而,解决了现有方法中存在的所确定的地层密度和电阻率关系的误差较大、准确度较低的技术问题,达到准确、快速地确定地层密度和电阻率关系的技术效果。
在一个实施方式中,为了减少所获取的多个第一采样点的电阻率数据和密度数据误差,例如,去除第一采样点的电阻率数据和密度数据的毛刺,在获取多个第一采样点的电阻率数据和密度数据后,所述方法还可以包括以下内容:
对所述多个第一采样点的电阻数据和密度数据进行低通滤波处理,得到处理后的第一采样点的电阻率数据和密度数据。
在本实施方式中,上述低通滤波处理具体可以包括以下至少之一:中值滤波、空间滤波、滑动平均滤波、高斯滤波等。当然需要说明的是,上述所列举的多种低通滤波处理只是为了更好地说明本申请实施方式,具体实施时,也可以根据具体情况和施工要求,选择使用其他的滤波处理方式进行相应的低通滤波处理。
相应的,根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体,具体可以包括:
根据所述地质分层数据,合并位于同一地层的处理后的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体。
在一个实施方式中,为了进一步提高处理速度,在上述根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体后,所述方法还可以包括以下内容。
按照第二间隔,从各个地层的第一数据体中,抽取多个第二采样点的电阻率数据和密度数据作为各个地层的第二数据体。
在本实施方式中,上述第二间隔具体可以为目标区域中测井的总数。具体的,例如,目标区域中有10口测井,则第二间隔可以设置为10;进而可以从各个地层的第一数据体中每隔10个点采集一个点作为第二采样点,并根据测井数据,获取第二采样点的电阻率数据和密度数据。如此,可以进一步减少数据处理量,提高实施效率。
相应的,根据各个地层的第一数据体,构建关于地层密度和电阻率关系的多个关系函数,包括根据各个地层的第二数据体,构建关于地层密度和电阻率关系的多个关系函数。
在一个实施方式中,为了进一步减少误差,提高所确定的目标地层密度和电阻率关系的准确度,在确定所述目标区域的地层密度和电阻率关系后,所述方法具体还可以包括以下内容。
S1:根据所述各个地层的第一数据体,建立密度电阻率交汇图。
S2:根据所述密度电阻率交汇图,对所述目标区域的地层密度和电阻率关系进行校验。
在本实施方式中,具体实施时,上述根据所述密度电阻率交汇图,对所述目标区域的地层密度和电阻率关系进行校验,可以包括:将所确定目标区域的地层密度和电阻率关系与所述密度电阻率交汇图进行比较,当两者较为一致时,则校验成功,认为所确定的目标区域的地层密度和电阻率关系较为准确;当两只不一致时,则校验失败,认为所确定的目标区域的地层密度和电阻率关系不准确,不符合要求,需要重新确定目标区域的地层密度和电阻率关系,直到所确定的目标区域的地层密度和电阻率关系的准确度较高满足施工要求为止。
在一个实施方式中,为了对目标区域进行更为准确的地质解释,在确定所述目标区域的地层密度和电阻率关系后,所述方法具体还可以包括以下内容。
S1:根据所述目标区域的电阻率和密度关系,对目标区域进行电磁和重力的联合反演,得到反演结果。
S2:根据所述反演结果,对目标区域进行地质解释和钻探。
在本实施方式中,具体实施时,可以将上述目标区域的电阻率和密度关系作为物性基础,对目标区域进行电磁和重力的联合反演;再根据反演结果对目标区域进行具体的地质解释和钻探。如此,可以使得所获得的目标区域的地层分辨率更高,提高地球物理勘探的精度。
本申请还提供一种可以较为准确地确定目标区域的地层密度和电阻率关系的方法,具体实施过程可以包括以下内容。
S1:获取目标区域的测井数据和地质分层数据;
S2:按照第一间隔,从所述测井数据中获取多个第一采样点的电阻率数据和密度数据;
S3:对所述多个第一采样点的电阻数据和密度数据进行低通滤波处理,得到处理后的第一采样点的电阻率数据和密度数据;
S4:根据所述地质分层数据,合并位于同一地层的处理后的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体;
S5:按照第二间隔,从各个地层的第一数据体中,抽取多个第二采样点的电阻率数据和密度数据作为各个地层的第二数据体;
S6:根据各个地层的第二数据体,构建关于地层密度和电阻率关系的多个关系函数;
S7:从所述多个关系函数中,确定出符合预设要求的关系函数,作为所述目标区域的地层密度和电阻率关系;
S8:根据所述各个地层的第一数据体,建立密度电阻率交汇图;
S9:根据所述密度电阻率交汇图,对所述目标区域的地层密度和电阻率关系进行校验。
从以上的描述中,可以看出,本申请实施方式提供的地层密度和电阻率关系的确定方法,通过具体考虑了地层密度数据、电阻率数据的相应特点,综合利用测井数据和地质分层数据,构建多个关系函数,进而从多个关系函数中确定出较为准确的关系函数作为地层密度和电阻率关系,从而,解决了现有方法中存在的所确定的地层密度和电阻率关系的误差较大、准确度较低的技术问题,达到准确、快速地确定地层密度和电阻率关系的技术效果;又通过分别计算各个关系函数的评价参数,根据评价参数确定出表征效果较好的关系函数作为目标区域的地层密度和电阻率关系,提高了所确定的地层密度和电阻率关系的准确度;还通过制作密度电阻率交汇图,利用密度电阻率交汇图对所确定的地层密度和电阻率关系进行校验,进一步提高了所确定的地层密度和电阻率关系的准确度。
基于同一发明构思,本发明实施方式中还提供了一种地层密度和电阻率关系的确定装置,如下面的实施方式所述。由于装置解决问题的原理与地层密度和电阻率关系的确定方法相似,因此地层密度和电阻率关系的确定装置的实施可以参见地层密度和电阻率关系的确定方法的实施,重复之处不再赘述。以下所使用的,术语“单元”或者“模块”可以实现预定功能的软件和/或硬件的组合。尽管以下实施例所描述的装置较佳地以软件来实现,但是硬件,或者软件和硬件的组合的实现也是可能并被构想的。请参阅图2,是本申请实施方式的地层密度和电阻率关系的确定装置的一种组成结构图,该装置可以包括:第一获取模块21、第二获取模块22、合并模块23、构建模块24、确定模块25。下面对该结构进行具体说明。
第一获取模块21,具体可以用于获取目标区域的测井数据和地质分层数据。
第二获取模块22,具体可以用于按照第一间隔,从测井数据中获取多个第一采样点的电阻率数据和密度数据。
合并模块23,具体可以用于根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体。
构建模块24,具体可以用于根据各个地层的第一数据体,构建关于地层密度和电阻率关系的多个关系函数。
确定模块25,具体可以用于从所述多个关系函数中,确定出符合预设要求的关系函数,作为所述目标区域的地层密度和电阻率关系。
在一个实施方式中,为了减少误差,所述装置还包括校验模块,其中,上述校验模块具体可以包括以下结构:
建立单元,具体可以用于根据所述各个地层的第一数据体,建立密度电阻率交汇图;
校验单元,具体可以用于根据所述密度电阻率交汇图,对所述目标区域的地层密度和电阻率关系进行校验。
在一个实施方式中,为了提高目标区域地层的分辨率,以便更好地对目标区域进行地质解释,所述装置具体还可以包括施工模块,具体可以用于根据所述目标区域的电阻率和密度关系,对目标区域进行电磁和重力的联合反演,得到反演结果;根据所述反演结果,对目标区域进行地质解释和钻探。
本说明书中的各个实施例均采用递进的方式描述,各个实施例之间相同相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。尤其,对于系统实施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
需要说明的是,上述实施方式阐明的系统、装置、模块或单元,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。为了描述的方便,在本说明书中,描述以上装置时以功能分为各种单元分别描述。当然,在实施本申请时可以把各单元的功能在同一个或多个软件和/或硬件中实现。
此外,在本说明书中,诸如第一和第二这样的形容词仅可以用于将一个元素或动作与另一元素或动作进行区分,而不必要求或暗示任何实际的这种关系或顺序。在环境允许的情况下,参照元素或部件或步骤(等)不应解释为局限于仅元素、部件、或步骤中的一个,而可以是元素、部件、或步骤中的一个或多个等。
从以上的描述中,可以看出,本申请实施方式提供的地层密度和电阻率关系的确定装置,具体考虑了地层密度数据、电阻率数据的相应特点,通过构建模块综合利用测井数据和地质分层数据,构建多个关系函数,进而通过确定模块从多个关系函数中确定出较为准确的关系函数作为地层密度和电阻率关系,从而,解决了现有方法中存在的所确定的地层密度和电阻率关系的误差较大、准确度较低的技术问题,达到准确、快速地确定地层密度和电阻率关系的技术效果;又通过确定模块分别计算各个关系函数的评价参数,根据评价参数确定出表征效果较好的关系函数作为目标区域的地层密度和电阻率关系,提高了所确定的地层密度和电阻率关系的准确度;还通过制作密度电阻率交汇图,利用密度电阻率交汇图对所确定的地层密度和电阻率关系进行校验,进一步提高了所确定的地层密度和电阻率关系的准确度。
在一个具体实施场景示例中,应用本申请提供的地层密度和电阻率关系的确定方法/装置,对某区域的地层密度和电阻率关系进行具体的确定。具体实施过程,可以参阅图3所示的在一个场景示例中应用本申请实施方式提供地层密度和电阻率关系的确定方法/装置的流程示意图,按照以下内容执行。
S1:收集研究区N口测井的电阻率、密度或者声波时差测井数据(即各个测井的测井数据),收集每口测井的地质分层数据,并对所获得的数据进行如下的数据处理与统计。
1)如果测井数据有密度测井数据,则可以根据测井数据中的密度测井数据确定对应密度数据,并记den,单位为g/cm3;如果测井数据没有密度测井数据,但有声波时差测井数据,则可以按照公式(1)将声波时差测井数据转换成密度数据:
den=a(b/AC)1/4 (1)
上式中,AC可以为声波时差测井数据,单位为μs/ft;den可以为声波转换的密度,单位为g/cm3;a可以为第一系数,b可以为第二系数,其中,a=0.31,b=0.3048×106。
2)根据测井数据的采样间隔(第一间隔),抽稀每口测井的电阻率、密度数据。其中采样间隔设置在5m左右(即按照第一间隔,从所述测井数据中获取多个第一采样点的电阻率数据和密度数据)。
3)取测井电阻率值的自然对数,记为lgR,对抽稀的对数电阻率和密度数据进行低通滤波。其中,上述低通滤波方法具体可以包含:中值滤波、空间滤波、滑动平均滤波、高斯滤波等。
4)按照测井分层的深度数据(根据地质分层数据确定)将N口井的相同地层对数电阻率和密度合并(得到各个地层的第一数据体),并按照N的间隔(即第二间隔)再次将每个地层的对数电阻率和密度抽稀,抽稀后的对数电阻率和密度分别记为lgRc和Denc(即第二数据体)。
S2:利用公式(2)求取各个层位密度-对数电阻率的指数拟合、线性、二次多项式等关系式(即确定多个关系函数):
上式中,y1为多个关系函数中第一关系函数的拟合密度,y2为多个关系函数中第二关系函数的拟合密度,y3为多个关系函数中第三关系函数的拟合密度,x为电阻率的对数,a0为第一拟合系数,a1为第二拟合系数,a2为第三拟合系数。
S3:根据公式(3)求取各个层位三种关系式的拟合优度(即各个关系函数的评价参数),挑选出拟合优度最大的关系式,即为当前层位最佳的密度-对数电阻率拟合(即目标区域的地层密度和电阻率关系)。其中,确定拟合优度的公式可以表示为:
上式中,r2可以为关系函数的评价参数,yobs可以为实测密度,y可以为关系函数对应的拟合密度,可以为关系函数对应的拟合密度的平均值。
S4:根据上述确定的最佳的密度-对数电阻率拟合,对该区域进行电磁和重力的联合反演,得到反演结果;根据所述反演结果,对该区域进行具体的地质解释和钻探。
通过上述的场景示例,验证了本申请实施方式提供的地层密度和电阻率关系的确定方法,通过考虑地层密度数据、电阻率数据的相应特点,综合利用测井数据和地质分层数据,构建多个关系函数,再从多个关系函数中确定出较为准确的关系函数作为地层密度和电阻率关系,确实可以解决现有方法中存在的所确定的地层密度和电阻率关系的误差较大、准确度较低的技术问题。
尽管本申请内容中提到不同的具体实施方式,但是,本申请并不局限于必须是行业标准或实施例所描述的情况等,某些行业标准或者使用自定义方式或实施例描述的实施基础上略加修改后的实施方案也可以实现上述实施例相同、等同或相近、或变形后可预料的实施效果。应用这些修改或变形后的数据获取、处理、输出、判断方式等的实施例,仍然可以属于本申请的可选实施方案范围之内。
虽然本申请提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创造性的手段可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境,甚至为分布式数据处理环境)。术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含,从而使得包括一系列要素的过程、方法、产品或者设备不仅包括那些要素,而且还包括没有明确列出的其他要素,或者是还包括为这种过程、方法、产品或者设备所固有的要素。在没有更多限制的情况下,并不排除在包括所述要素的过程、方法、产品或者设备中还存在另外的相同或等同要素。
上述实施例阐明的装置或模块等,具体可以由计算机芯片或实体实现,或者由具有某种功能的产品来实现。为了描述的方便,描述以上装置时以功能分为各种模块分别描述。当然,在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现,也可以将实现同一功能的模块由多个子模块的组合实现等。以上所描述的装置实施例仅仅是示意性的,例如,所述模块的划分,仅仅为一种逻辑功能划分,实际实现时可以有另外的划分方式,例如多个模块或组件可以结合或者可以集成到另一个系统,或一些特征可以忽略,或不执行。
本领域技术人员也知道,除了以纯计算机可读程序代码方式实现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件模块又可以是硬件部件内的结构。
本申请可以在由计算机执行的计算机可执行指令的一般上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请,在这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可借助软件加必需的通用硬件平台的方式来实现。基于这样的理解,本申请的技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来,该计算机软件产品可以存储在存储介质中,如ROM/RAM、磁碟、光盘等,包括若干指令用以使得一台计算机设备(可以是个人计算机,移动终端,服务器,或者网络设备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
本说明书中的各个实施例采用递进的方式描述,各个实施例之间相同或相似的部分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。本申请可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算机、手持设备或便携式设备、平板型设备、多处理器系统、基于微处理器的系统、置顶盒、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布式计算环境等等。
虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和变化而不脱离本申请的精神,希望所附的实施方式包括这些变形和变化而不脱离本申请。
Claims (12)
1.一种地层密度和电阻率关系的确定方法,其特征在于,包括:
获取目标区域的测井数据和地质分层数据;
按照第一间隔,从所述测井数据中获取多个第一采样点的电阻率数据和密度数据;
根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体;
根据各个地层的第一数据体,构建关于地层密度和电阻率关系的多个关系函数;
从所述多个关系函数中,确定出符合预设要求的关系函数,作为所述目标区域的地层密度和电阻率关系;
其中,根据各个地层的第一数据体,构建关于地层密度和电阻率关系的多个关系函数,包括:
根据各个地层的第一数据体,按照以下公式,构建关于地层密度和电阻率关系的多个关系函数,其中,所述多个关系函数包括以下至少之一:第一关系函数、第二关系函数、第三关系函数:
y2=a0+a1x
y3=a0+a1x+a2x2
上式中,y1为多个关系函数中第一关系函数的拟合密度,y2为多个关系函数中第二关系函数的拟合密度,y3为多个关系函数中第三关系函数的拟合密度,x为电阻率的对数,a0为第一拟合系数,a1为第二拟合系数,a3为第三拟合系数。
2.根据权利要求1所述的方法,其特征在于,所述第一间隔为5米。
3.根据权利要求1所述的方法,其特征在于,从所述测井数据中,获取多个第一采样点的密度数据,包括:
确定所述测井数据中是否包括密度测井数据;
在所述测井数据中包括密度测井数据的情况下,根据所述密度测井数据,确定所述多个第一采样点的密度数据;
在所述测井数据中不包括密度测井数据的情况下,根据声波测井数据,确定所述多个第一采样点的密度数据。
4.根据权利要求3所述的方法,其特征在于,根据声波测井数据,确定所述多个第一采样点的密度数据,包括:
按照以下公式确定所述多个第一采样点的密度数据:
den=a(b/AC)1/4
上式中,den为密度数据,AC为声波测井数据,a为第一系数,b为第二系数。
5.根据权利要求1所述的方法,其特征在于,从所述多个关系函数中,确定出符合预设要求的关系函数,包括:
根据所述多个关系函数,确定所述多个关系函数中各个关系函数的评价参数;
根据所述各个关系函数的评价参数,确定评价参数中数值最大的关系函数,作为所述符合预设要求的关系函数。
6.根据权利要求5所述的方法,其特征在于,根据所述多个关系函数,确定所述多个关系函数中各个关系函数的评价参数,包括:
按照以下公式,确定关系函数的评价参数:
上式中,r2为关系函数的评价参数,yobs为实测密度,y为关系函数对应的拟合密度,为关系函数对应的拟合密度的平均值。
7.根据权利要求1所述的方法,其特征在于,在获取多个第一采样点的电阻率数据和密度数据后,所述方法还包括:
对所述多个第一采样点的电阻率数据和密度数据进行低通滤波处理,得到处理后的第一采样点的电阻率数据和密度数据;
相应的,根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体,包括:
根据所述地质分层数据,合并位于同一地层的处理后的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体。
8.根据权利要求1所述的方法,其特征在于,根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体后,所述方法还包括:
按照第二间隔,从各个地层的第一数据体中,抽取多个第二采样点的电阻率数据和密度数据作为各个地层的第二数据体;
相应的,根据各个地层的第二数据体,构建关于地层密度和电阻率关系的多个关系函数。
9.根据权利要求1所述的方法,其特征在于,在确定所述目标区域的地层密度和电阻率关系后,所述方法还包括:
根据所述各个地层的第一数据体,建立密度电阻率交汇图;
根据所述密度电阻率交汇图,对所述目标区域的地层密度和电阻率关系进行校验。
10.根据权利要求1所述的方法,其特征在于,在确定所述目标区域的地层密度和电阻率关系后,所述方法还包括:
根据所述目标区域的电阻率和密度关系,对目标区域进行电磁和重力的联合反演,得到反演结果;
根据所述反演结果,对目标区域进行地质解释和钻探。
11.一种地层密度和电阻率关系的确定装置,其特征在于,包括:
第一获取模块,用于获取目标区域的测井数据和地质分层数据;
第二获取模块,用于按照第一间隔,从测井数据中获取多个第一采样点的电阻率数据和密度数据;
合并模块,用于根据所述地质分层数据,合并位于同一地层的第一采样点的电阻率数据和密度数据,得到所述目标区域中各个地层的第一数据体;
构建模块,用于根据各个地层的第一数据体,构建关于地层密度和电阻率关系的多个关系函数;
确定模块,用于从所述多个关系函数中,确定出符合预设要求的关系函数,作为所述目标区域的地层密度和电阻率关系;
其中,所述构建模块具体用于根据各个地层的第一数据体,按照以下公式,构建关于地层密度和电阻率关系的多个关系函数,其中,所述多个关系函数包括以下至少之一:第一关系函数、第二关系函数、第三关系函数:
y2=a0+a1x
y3=a0+a1x+a2x2
上式中,y1为多个关系函数中第一关系函数的拟合密度,y2为多个关系函数中第二关系函数的拟合密度,y3为多个关系函数中第三关系函数的拟合密度,x为电阻率的对数,a0为第一拟合系数,a1为第二拟合系数,a3为第三拟合系数。
12.根据权利要求11所述的装置,其特征在于,所述装置包括施工模块,用于根据所述目标区域的电阻率和密度关系,对目标区域进行电磁和重力的联合反演,得到反演结果;根据所述反演结果,对目标区域进行地质解释和钻探。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710811319.6A CN107807409B (zh) | 2017-09-11 | 2017-09-11 | 地层密度和电阻率关系的确定方法和装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710811319.6A CN107807409B (zh) | 2017-09-11 | 2017-09-11 | 地层密度和电阻率关系的确定方法和装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107807409A CN107807409A (zh) | 2018-03-16 |
CN107807409B true CN107807409B (zh) | 2019-08-06 |
Family
ID=61569861
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710811319.6A Active CN107807409B (zh) | 2017-09-11 | 2017-09-11 | 地层密度和电阻率关系的确定方法和装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107807409B (zh) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108828681B (zh) * | 2018-08-01 | 2020-01-07 | 中国石油天然气集团有限公司 | 地层电阻率和极化率的确定方法和装置 |
CN113970787B (zh) * | 2020-07-22 | 2023-05-16 | 中国石油化工股份有限公司 | 物性参数反演方法、装置、计算机设备和存储介质 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0922239A1 (fr) * | 1997-06-13 | 1999-06-16 | Elf Exploration Production | Methode de caracterisation de la coherence de mesures de caracteristiques d'un milieu |
CN103809216A (zh) * | 2012-11-12 | 2014-05-21 | 中国石油天然气集团公司 | 一种电阻率数据与地震数据联合速度建场方法 |
CN105242317A (zh) * | 2015-10-14 | 2016-01-13 | 中国石油天然气集团公司 | 一种纵波速度的确定方法及装置 |
CN106522932A (zh) * | 2016-09-28 | 2017-03-22 | 中国石油天然气股份有限公司 | 获取水平井各向异性地层电阻率的方法及装置 |
-
2017
- 2017-09-11 CN CN201710811319.6A patent/CN107807409B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0922239A1 (fr) * | 1997-06-13 | 1999-06-16 | Elf Exploration Production | Methode de caracterisation de la coherence de mesures de caracteristiques d'un milieu |
CN103809216A (zh) * | 2012-11-12 | 2014-05-21 | 中国石油天然气集团公司 | 一种电阻率数据与地震数据联合速度建场方法 |
CN105242317A (zh) * | 2015-10-14 | 2016-01-13 | 中国石油天然气集团公司 | 一种纵波速度的确定方法及装置 |
CN106522932A (zh) * | 2016-09-28 | 2017-03-22 | 中国石油天然气股份有限公司 | 获取水平井各向异性地层电阻率的方法及装置 |
Non-Patent Citations (1)
Title |
---|
MT-重力模拟退火联合反演研究;陈华根 等;《地球物理学报》;20120227;第55卷(第2期);第663-670页 |
Also Published As
Publication number | Publication date |
---|---|
CN107807409A (zh) | 2018-03-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107817535B (zh) | 薄夹层的确定方法和装置 | |
CN105675635B (zh) | 致密岩石组分相对含量及脆性指数确定方法和装置 | |
CN104375188B (zh) | 一种地震波透射衰减补偿方法及装置 | |
US9523781B2 (en) | Normalization seismic attribute | |
CN109100796A (zh) | 一种天然气水合物地震数据处理方法及装置 | |
CN104122581B (zh) | 一种叠后声波阻抗反演方法 | |
CN108181649B (zh) | 一种介质分界面振幅补偿方法及装置 | |
CN107450108B (zh) | 甜点区的确定方法和装置 | |
CN105044775B (zh) | 一种地震流体反演识别方法和装置 | |
CN107366534A (zh) | 粗化渗透率的确定方法和装置 | |
CN103439740A (zh) | 基于偶极地震子波多重积分的相对阻抗预测的方法及装置 | |
CN108318937A (zh) | 地质解释方法和装置 | |
CN104502996A (zh) | 一种密度曲线校正方法及系统 | |
CN107807409B (zh) | 地层密度和电阻率关系的确定方法和装置 | |
CN107576985B (zh) | 一种地震反演的方法和装置 | |
CN108089238A (zh) | 一种利用综合地球物理资料确定深层裂谷的方法及装置 | |
Zhou et al. | Data driven modeling and prediction for reservoir characterization using seismic attribute analyses and big data analytics | |
CN107229075B (zh) | 深度域地震子波的确定方法和装置 | |
CN105301638B (zh) | 一种提取风化层底界面的方法和装置 | |
CN106501872B (zh) | 一种裂缝储层地应力特征的计算方法及装置 | |
CN108415073A (zh) | 角度域逆散射偏移成像方法及装置 | |
CN107870361B (zh) | 一种地震回折波层析成像方法、装置及终端设备 | |
CN107991716B (zh) | 薄层的沉积相图和沉积体厚度的确定方法和装置 | |
CN106842316B (zh) | 裂缝确定方法和装置 | |
Shan et al. | Spectral decomposition and a waveform cluster to characterize strongly heterogeneous paleokarst reservoirs in the Tarim basin, China |
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 |