CN103410504B - 一种水平井/大斜度井确定地层真电阻率的方法及装置 - Google Patents
一种水平井/大斜度井确定地层真电阻率的方法及装置 Download PDFInfo
- Publication number
- CN103410504B CN103410504B CN201310308564.7A CN201310308564A CN103410504B CN 103410504 B CN103410504 B CN 103410504B CN 201310308564 A CN201310308564 A CN 201310308564A CN 103410504 B CN103410504 B CN 103410504B
- Authority
- CN
- China
- Prior art keywords
- well
- resistivity
- logging
- drilling
- log
- 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
- 238000000034 method Methods 0.000 title claims abstract description 29
- 238000005553 drilling Methods 0.000 claims abstract description 114
- 230000004044 response Effects 0.000 claims abstract description 94
- 230000002452 interceptive effect Effects 0.000 claims abstract description 9
- 230000015572 biosynthetic process Effects 0.000 claims description 49
- 230000008859 change Effects 0.000 claims description 13
- 239000012530 fluid Substances 0.000 claims description 8
- 230000005251 gamma ray Effects 0.000 claims description 8
- 239000004576 sand Substances 0.000 claims description 8
- 238000004364 calculation method Methods 0.000 claims description 5
- 208000035126 Facies Diseases 0.000 claims description 4
- 239000000463 material Substances 0.000 claims description 4
- 241001074085 Scophthalmus aquosus Species 0.000 claims description 3
- 238000003384 imaging method Methods 0.000 claims description 3
- 238000011156 evaluation Methods 0.000 abstract description 4
- 238000004088 simulation Methods 0.000 abstract description 4
- 238000012360 testing method Methods 0.000 abstract description 3
- 238000005259 measurement Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 3
- 230000005611 electricity Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 3
- 239000011435 rock Substances 0.000 description 2
- 239000007787 solid Substances 0.000 description 2
- 101000918657 Homo sapiens L-xylulose reductase Proteins 0.000 description 1
- 102100029137 L-xylulose reductase Human genes 0.000 description 1
- 102220350010 c.119C>A Human genes 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000013329 compounding Methods 0.000 description 1
- 230000008021 deposition Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000003208 petroleum Substances 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 102220206869 rs1057522105 Human genes 0.000 description 1
- 102220086595 rs267607538 Human genes 0.000 description 1
- 102220076630 rs796052858 Human genes 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种水平井/大斜度井确定地层真电阻率的方法及装置,该方法包括:获取随钻电阻率测井资料、辅助测井曲线资料、邻井或导眼井的测井资料以及井附近区域地震资料,确定水平井段/大斜度井段钻探的目的层,对目的层进行分层建立初始地层模型并确定初始地层模型的电阻率初值,生成对应层的随钻电阻率模拟响应曲线;判断对应层的随钻电阻率模拟响应曲线与对应层的随钻电阻率原始响应曲线是否一致:一致,将初始地层模型的电阻率初值做为地层真电阻率值输出;判断不一致,通过交互式反演修改电阻率初值。本发明大大提高油层的识别能力和含油饱和度评价的精度,为射孔试油提供了较为坚实的依据。
Description
技术领域
本发明关于石油测井技术,特别是关于确定地层真电阻率的技术,具体的讲是一种水平井/大斜度井确定地层真电阻率的方法及装置。
背景技术
水平井或大斜度井,由于受到泥浆侵入、井眼、围岩、井斜及地层电各向异性等的影响,使得测量的电阻率偏离真实电阻率,特别是在砂泥岩互层等各向异性地层,从而使得通过电阻率曲线进行油层识别的能力降低,计算的含油饱和度误差较大。因此,需要利用反演方法还原真实地层面目。
但是,传统的反演方法一直以来难以实用化,其原因在于它们都是纯数学方法,只要满足模拟出的测井曲线与实测曲线误差小于给定的ε即可,而不管在岩石物理、地质环境及工程条件等方面是否合理,因此,利用传统反演方法得到的结果不够合理。
发明内容
为了提高油层的识别能力和含油饱和度的评价精度,本发明实施例提供了一种水平井/大斜度井确定地层真电阻率的方法,该包括:
步骤1,获取待评价的井的随钻电阻率测井资料、辅助测井曲线资料、邻井或导眼井的测井资料以及待评价井所在区域的地震资料;
步骤2,根据所述随钻电阻率曲线测井资料、辅助测井曲线资料和邻井或导眼井的测井资料确定水平井段/大斜度井段钻探的目的层、随钻电阻率原始响应曲线;
步骤3,根据所述辅助测井曲线资料、邻井或导眼井的测井资料对所述目的层进行分层建立初始地层模型并确定初始地层模型的电阻率初值;
步骤4,根据所述初始地层模型和电阻率初值,应用快速正演法生成对应层的随钻电阻率模拟响应曲线;
步骤5,判断所述对应层的随钻电阻率模拟响应曲线与对应层的随钻电阻率原始响应曲线是否一致:判断一致,执行步骤6;判断不一致,执行步骤7;
步骤6将所述初始地层模型的电阻率初值做为地层真电阻率值输出;
步骤7,根据所述对应层的随钻电阻率原始响应曲线、对应层的地质特征通过交互式反演修改所述电阻率初值,重复步骤4至步骤7。
此外,本发明还公开了一种水平井/大斜度井确定地层真电阻率的装置,该装置包括:
资料获取模块,用于获取待测井的随钻电阻率测井资料、辅助测井曲线资料、邻井或导眼井的测井资料以及待测井附近区域地震资料,所述随钻电阻率测井资料包括随钻电阻率测井曲线,所述辅助测井曲线资料包括辅助测井曲线,所述邻井或导眼井的测井资料包括邻井或导眼井测井曲线;
目的层确定模块,应用根据所述随钻电阻率曲线测井资料、辅助测井曲线资料和邻井或导眼井的测井资料确定水平井段/大斜度井段钻探的目的层、随钻电阻率原始响应曲线;
模型建立模块,用于根据所述辅助测井曲线资料、邻井或导眼井的测井资料对所述目的层进行分层建立初始地层模型并确定初始地层模型的电阻率初值;
模拟响应曲线生成模块,用于根据所述初始地层模型和电阻率初值,应用快速正演法生成对应层的随钻电阻率模拟响应曲线;
判断模块,用于判断所述对应层的随钻电阻率模拟响应曲线与对应层的随钻电阻率原始响应曲线是否一致;
修正模块,用于根据所述对应层的随钻电阻率原始响应曲线、对应层的地质特征通过交互式反演修改所述电阻率初值;
输出模块,用于输出地层真电阻率值。
本发明利用随钻电阻率测井资料计算地层真电阻率的方法及装置充分考虑了地质背景、油藏条件、工程环境等因素,主动修改地层模型,以模拟出的测井曲线与实测曲线是否吻合作约束,从而使得反演结果更为合理。因此,利用本发明的方法可以最大限度的恢复地层真实面目,大大提高油层的识别能力和含油饱和度评价的精度,为射孔试油提供了较为坚实的依据。
为让本发明的上述和其他目的、特征和优点能更明显易懂,下文特举较佳实施例,并配合所附图式,作详细说明如下。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案,下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提下,还可以根据这些附图获得其他的附图。
图1为本发明一种水平井/大斜度井确定地层真电阻率的方法的流程图;
图2为对地层进行分层的步骤流程图;
图3为修改初始的电阻率的流程图;
图4为本发明水平井/大斜度井确定地层真电阻率装置的结构框图;
图5为目的层确定模块的框图;
图6为模型建立模块的框图;
图7为修正模块的框图;
图8为利用本发明确定地层真电阻率的方法的随钻电阻率测井解释成果图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
如图1所示,本发明提供了一种水平井/大斜度井确定地层真电阻率的方法,该包括:
步骤S101,获取待评价的井的随钻电阻率测井资料、辅助测井曲线资料、邻井或导眼井的测井资料以及待评价井所在区域的地震资料;
步骤S102,根据随钻电阻率曲线测井资料、辅助测井曲线资料和邻井或导眼井的测井资料确定水平井段/大斜度井段钻探的目的层;本实施例中,根据井眼轨迹将随钻电阻率测井曲线和辅助测井曲线投影到垂直深度,确定随钻电阻率测井曲线和辅助测井曲线在直井模式时的响应特征;根据随钻电阻率测井曲线和辅助测井曲线在直井模式时的响应特征与邻井或导眼井的测井曲线响应特征对比,确定水平井段/大斜度井段钻探的目的层;本发明中邻井或导眼井的测井曲线响应特征,即为原始测井仪器测量得到的,然后根据解释总结每一目的层对应曲线值的范围。本发明技术方案中的测井原始曲线,都是通过测井仪器得到的。
步骤S103,根据辅助测井曲线资料、邻井或导眼井的测井资料对目的层进行分层建立初始地层模型并确定初始地层模型的电阻率初值;
步骤S104,根据初始地层模型和电阻率初值,应用快速正演法生成对应层的随钻电阻率模拟响应曲线;
步骤S105,判断对应层的随钻电阻率模拟响应曲线与对应层的随钻电阻率原始响应曲线是否一致:判断一致,执行步骤S106;判断不一致,执行步骤S107;
步骤S106,将初始地层模型的电阻率初值做为地层真电阻率值输出;
步骤S107,根据对应层的随钻电阻率原始响应曲线、对应层的地质特征通过交互式反演修改所述电阻率初值,重复步骤S104至步骤S107。
在本实施例中,步骤S101中的辅助测井曲线资料包括自然伽马测井曲线、井径测井曲线、密度测井曲线以及自然电位测井曲线。辅助测井曲线资料主要用于在步骤S102中确定目的层段以及在步骤S103中对目的层井段进行分层。根据地层沉积规律,利用随钻电阻率测井资料结合这些辅助测井资料,将地层进行精细分层,特别是地层中的薄夹层,需要仔细予以识别。
步骤S103的对地层进行分层的步骤如图2所示,其包括:
步骤S1031,根据所述邻井或导眼井测井曲线特征依据斜率变化趋势确定层界面。此步骤中,斜率是指曲线上下两点幅度差与深度差的比值;当斜率由正变负或者由负变正时,该深度点作为特征点,两个特征点的中点即为层界面。
步骤S1032,根据所述辅助测井曲线资料和确定的层界面,分别将所述邻井或导眼井电阻率测井曲线特征不明显的泥岩夹层、致密层、垮塌层分出来。此处的辅助测井曲线资料包括自然伽马测井曲线、井径测井曲线、密度测井曲线以及自然电位测井曲线。
步骤S1033,根据所述邻井或导眼井电阻率曲线资料和确定的各地层单元,确定各对应层的电阻率初值。
本实施例中还包括,根据所述待测井附近区域地震资料、微电阻率扫描成像测井资料或邻井测井资料,确定目的层段各对应层的初始地层倾角值。在建立初始地层模型时,地层是水平的,但实际中,地层不一定是水平的,因此需要根据区域地震资料确定地层大致的倾向和倾角。如果有电阻率成像,可以大致的判断井轨迹与地层的穿层关系。
在本实施例的步骤S104中,根据对应层的初始的电阻率,生成对应层的随钻电阻率模拟响应曲线,具体包括:根据对应层的初始的电阻率,结合随钻电阻率测井仪器的仪器模型参数,通过快速正演法生成对应层的随钻电阻率模拟响应曲线;其中,随钻电阻率测井仪器的仪器模型参数包括:仪器长度、仪器半径、各单元线圈距、线圈匝数、工作频率。
在本实施例的步骤S105中,判断对应层的随钻电阻率模拟响应曲线与对应层的所述随钻电阻率原始响应曲线是否一致,包括:
设定对应层的随钻电阻率模拟响应曲线与对应层的所随钻电阻率原始响应曲线间的相关系数为S:
N-目的层段内的测量点数;
X-所述对应层的随钻电阻率原始响应曲线目的层段内的数据集;
Y-所述对应层的随钻电阻率模拟响应曲线目的层段内的数据集;
当所述相关系数S的计算结果大于0.9,即认为对应层的随钻电阻率模拟响应曲线与对应层的所述随钻电阻率原始响应曲线一致。
本实施例中,相关系数S为根据实际情况设定的一个阈值,随钻电阻率模拟响应曲线与随钻电阻率原始响应曲线进行对应,是指同线圈距同工作频率。在本实施例中将S值设定为0.9。当然,本发明不限于此,还可根据其他实际情况将S设置为其他值,例如0.8、0.88、0.91等等。
本发明的步骤S107中,根据所述随钻电阻率原始响应曲线,结合对应层的地质特征,通过交互式反演修改电阻率初值(包括水平电阻率和垂直电阻率),如图3所示,其包括:
步骤S1071,根据对应层所在井的邻井资料以及本井的录井资料、气测资料、取心资料判断对应层的流体性质以及所述流体性质下的地层电阻率范围;
步骤S1072,根据对应层的倾角测井资料或地震资料判断对应层的倾向及倾角范围;
步骤S1073,根据待评价井所在区域内地层对比结合地区沉积相特征,判断对应层砂体厚度在径向上的变化;
步骤S1074,根据对应层的随钻电阻率模拟响应曲线与对应的随钻电阻率原始响应曲线间的相关程度,并结合所述地层电阻率范围、倾向、倾角范围以及砂体径向上的变化,修改所述初始的电阻率、地层倾向、地层倾角和地层边界。
其中,本井的录井资料、气测资料、取心资料、倾角测井资料或地震资料都是实际测量得到的,在应用本专利方法解释实际井时,这些资料都要预先收集,若油田本身没有测量这些资料,则不需要这些资料,这些资料都是辅助判断手段。
本发明还公开了一种水平井/大斜度井确定地层真电阻率的装置,如图4所示,本实施例的水平井/大斜度井确定地层真电阻率的装置40包括:
资料获取模块401,用于获取待评价井的随钻电阻率测井资料、辅助测井曲线资料、邻井或导眼井的测井资料以及待评价井所在区域的地震资料;
目的层确定模块402,应用根据所述随钻电阻率曲线测井资料、辅助测井曲线资料和邻井或导眼井的测井资料确定水平井段/大斜度井段钻探的目的层;
模型建立模块403,用于根据所述辅助测井曲线资料、邻井或导眼井的测井资料对所述目的层进行分层建立初始地层模型并确定初始地层模型的电阻率初值;
模拟响应曲线生成模块404,用于根据所述初始地层模型和电阻率初值,应用快速正演法生成对应层的随钻电阻率模拟响应曲线;
判断模块405,用于判断所述对应层的随钻电阻率模拟响应曲线与对应层的随钻电阻率原始响应曲线是否一致;
修正模块406,用于根据所述对应层的随钻电阻率原始响应曲线、对应层的地质特征通过交互式反演修改所述电阻率初值;
输出模块407,用于输出地层真电阻率值。
本实施例中,如图5所示,目的层确定模块402包括:
响应特征生成单元4021,用于根据井眼轨迹将随钻电阻率测井曲线和辅助测井曲线投影到垂直深度,确定随钻电阻率测井曲线和辅助测井曲线在直井模式时的响应特征;
目的层确定单元4022,用于根据随钻电阻率测井曲线和辅助测井曲线在直井模式时的响应特征与邻井或导眼井的测井曲线响应特征对比,确定水平井段/大斜度井段钻探的目的层。
如图6所示,模型建立模块403包括:
层界面确定单元4031,用于根据邻井或导眼井测井曲线特征依据邻井或导眼井测井曲线的斜率变化确定层界面;
分层单元4032,用于根据辅助测井曲线资料和确定的层界面,将所述随钻电阻率测井曲线特征不明显的泥岩夹层、致密层、垮塌层分出来;
电阻率确定单元4033,用于根据邻井或导眼井测井资料和确定的层界面,确定初始地层模型的电阻率初值;
初始地层倾角值确定单元4034,用于根据待测井附近区域地震资料、微电阻率扫描成像测井资料或邻井测井资料、确定的层界面,确定目的层段各对应层的初始地层倾角值。
如图7所示,修正模块406包括:
地层电阻率范围确定单元4061,用于根据对应层所在井的邻井测井资料以及本井的录井资料、气测资料、取心资料判断所述对应层的流体性质以及所述流体性质下的地层电阻率范围;
倾角范围确定单元4062,用于根据所述对应层的倾角测井资料或地震资料判断对应层的倾向及倾角范围;
径向变化判断单元4063,用于根据区域内地层对比结合地区沉积相特征,判断对应层砂体厚度在径向上的变化;
修改单元4064,用于根据所述对应层的随钻电阻率模拟响应曲线与对应的随钻电阻率原始响应曲线间的相关程度,并结合地层电阻率范围、倾向、倾角范围以及砂体径向上的变化,修改所述初始的水平电阻率、垂直电阻率、地层倾向、地层倾角和地层边界。
具体实施方式:
图8为利用本发明的利用随钻电阻率测井资料计算地层真电阻率的方法和装置生成的随钻电阻率测井解释成果图。
在此具体实施例中,选用水平井中普遍应用的随钻电阻率测井仪的测井资料作处理对象。图8分四个象限,左上角为曲线道头,右上角为曲线水平显示,左下角为曲线的垂直投影,右下角为地层模型(颜色越深表示地层电阻率值越低)。在右上角中:
第一道是自然伽马曲线(GR),用曲线801表示自然伽马曲线。
第二道是深度道,表示测量井段的距离井口的横向距离,其中括号中的深度表示测量深度。
第三道是地层模型道,表示经分层后各层的水平电阻率(FormationResistivityH)和垂直电阻率(FormationResistivityV),通过该道可修改地层模型值。其中,实线表示水平电阻率,虚线表示垂直电阻率。
第四道是反演后模拟的随钻电阻率与原始随钻电阻率响应曲线的对比结果,P16H、P22H、P28H、P34H、P40H分别表示源距为16、22、28、34、40英尺时,频率在2MHz下的相位电阻率曲线,P16Hm、P22Hm、P28Hm、P34Hm、P40Hm分别表示模拟出的对应不同源距的随钻电阻率模拟响应曲线。曲线中从上到下依次排列从大到小的不同源距的曲线,其中,虚线表示模拟出的随钻电阻率模拟响应曲线,实线表示测量出的随钻电阻率原始响应曲线。
右下角表示的是地层模型,显示的是井眼轨迹与地层关系,其中颜色的深浅表示地层电阻率值的高低。
当人机交互时,可修改第三道的地层模型值,以及右下角的地层模型层边界面,使得第四道的随钻电阻率模拟响应曲线与实测的随钻电阻率原始响应曲线相一致时,表示此时的地层模型已接近真实的地层,当满足相关系数S时即可输出所求的各层的电阻率。
本发明的此种方法可以大大提高油层的识别能力和含油饱和度的评价精度,为射孔和试油方案的设计提供坚实的依据。
本发明实施例的利用随钻电阻率测井资料计算地层真电阻率的方法及装置,充分考虑了地质背景、油藏条件、工程环境等因素,主动修改地层模型,以模拟出的测井曲线与实测曲线是否吻合作约束,从而使得计算结果更为合理。因此,利用本发明的参数计算方法及装置可以最大限度地恢复地层真实面目,大大提高油层的识别能力和含油饱和度评价的精度。
本发明中应用了具体实施例对本发明的原理及实施方式进行了阐述,以上实施例的说明只是用于帮助理解本发明的方法及其核心思想;同时,对于本领域的一般技术人员,依据本发明的思想,在具体实施方式及应用范围上均会有改变之处,综上所述,本说明书内容不应理解为对本发明的限制。
Claims (16)
1.一种水平井/大斜度井确定地层真电阻率的方法,其特征在于,所述的方法包括:
步骤1,获取待评价井的随钻电阻率曲线测井资料、辅助测井曲线资料、邻井或导眼井的测井资料以及待评价井所在区域的地震资料,所述随钻电阻率曲线测井资料包括随钻电阻率测井曲线;
步骤2,根据所述随钻电阻率曲线测井资料、辅助测井曲线资料和邻井或导眼井的测井资料确定水平井段/大斜度井段钻探的目的层;
步骤3,根据所述辅助测井曲线资料、邻井或导眼井的测井资料对所述目的层进行分层建立初始地层模型并确定初始地层模型的电阻率初值;
步骤4,根据所述初始地层模型和电阻率初值,应用快速正演法生成对应层的随钻电阻率模拟响应曲线;
步骤5,判断所述对应层的随钻电阻率模拟响应曲线与对应层的随钻电阻率原始响应曲线是否一致:判断一致,执行步骤6;判断不一致,执行步骤7;
步骤6将所述初始地层模型的电阻率初值作为地层真电阻率值输出;
步骤7,根据所述对应层的随钻电阻率原始响应曲线、对应层的地质特征通过交互式反演修改所述电阻率初值,重复步骤4至步骤7;
其中,所述的根据所述对应层的随钻电阻率原始响应曲线、对应层的地质特征通过交互式反演修改所述电阻率初值包括:
根据对应层所在井的邻井测井资料以及本井的录井资料、气测资料、取心资料判断所述对应层的流体性质以及所述流体性质下的地层电阻率范围;
根据所述对应层的倾角测井资料或地震资料判断对应层的倾向及倾角范围;
根据区域内地层对比结合地区沉积相特征,判断对应层砂体厚度在径向上的变化;
根据所述对应层的随钻电阻率模拟响应曲线与对应层的随钻电阻率原始响应曲线间的相关程度,并结合所述地层电阻率范围、倾向、倾角范围以及砂体厚度在径向上的变化,修改所述电阻率初值、地层倾向、地层倾角和地层边界。
2.如权利要求1所述的水平井/大斜度井确定地层真电阻率的方法,其特征在于,
所述的辅助测井曲线资料包括:自然伽马测井曲线、井径测井曲线、密度测井曲线、自然电位测井曲线、声波测井曲线、中子测井曲线、微电阻率扫描成像测井资料、方位伽马成像资料或者方位密度成像资料;
所述的邻井或导眼井的测井资料还包括:自然伽马测井曲线、电阻率测井曲线。
3.如权利要求1所述的水平井/大斜度井确定地层真电阻率的方法,其特征在于,
所述的电阻率初值包括:水平电阻率和垂直电阻率。
4.如权利要求2所述的水平井/大斜度井确定地层真电阻率的方法,其特征在于,所述的根据所述随钻电阻率曲线测井资料、辅助测井曲线资料和邻井或导眼井的测井资料确定水平井段/大斜度井段钻探的目的层包括:
根据井眼轨迹将随钻电阻率测井曲线和辅助测井曲线投影到垂直深度,确定随钻电阻率测井曲线和辅助测井曲线在直井模式时的响应特征;
根据随钻电阻率测井曲线和辅助测井曲线在直井模式时的响应特征与邻井或导眼井的测井曲线响应特征对比,确定水平井段/大斜度井段钻探的目的层。
5.如权利要求3或4所述的水平井/大斜度井确定地层真电阻率的方法,其特征在于,所述的根据所述辅助测井曲线资料、邻井或导眼井的测井资料对所述目的层进行分层建立初始地层模型并确定初始地层模型的电阻率初值包括:
根据邻井或导眼井测井曲线特征依据邻井或导眼井的测井曲线的斜率变化确定层界面;
根据所述辅助测井曲线资料和确定的层界面,将所述随钻电阻率测井曲线特征不明显的泥岩夹层、致密层、垮塌层分出来;
根据所述邻井或导眼井测井资料和确定的层界面,确定初始地层模型的电阻率初值;
根据所述待评价井所在区域地震资料、微电阻率扫描成像测井资料或邻井测井资料,确定目的层段各对应层的初始地层倾角值。
6.如权利要求1所述的水平井/大斜度井确定地层真电阻率的方法,其特征在于,所述的根据所述初始地层模型和电阻率初值,应用快速正演法生成对应层的随钻电阻率模拟响应曲线包括:
根据所述的电阻率初值和初始地层模型,结合随钻电阻率测井仪器的仪器模型参数,通过快速正演法生成对应层的随钻电阻率模拟响应曲线;
其中,所述随钻电阻率测井仪器的仪器模型参数包括:仪器长度、仪器半径、各阵列单元线圈距、线圈匝数、工作频率。
7.如权利要求1所述的水平井/大斜度井确定地层真电阻率的方法,其特征在于,所述的判断所述对应层的随钻电阻率模拟响应曲线与对应层的随钻电阻率原始响应曲线是否一致包括:
设定对应层的随钻电阻率模拟响应曲线与对应层的所述随钻电阻率原始响应曲线间的相关系数为S:
N-目的层段内的测量点数;
X-所述对应层的随钻电阻率原始响应曲线目的层段内的数据集;
Y-所述对应层的随钻电阻率模拟响应曲线目的层段内的数据集;
当所述相关系数S的计算结果大于预设阈值,判断对应层的随钻电阻率模拟响应曲线与对应层的所述随钻电阻率原始响应曲线一致。
8.如权利要求7所述的水平井/大斜度井确定地层真电阻率的方法,其特征在于,所述的预设阈值为0.9。
9.一种水平井/大斜度井确定地层真电阻率的装置,其特征在于,所述的装置包括:
资料获取模块,用于获取待评价井的随钻电阻率曲线测井资料、辅助测井曲线资料、邻井或导眼井的测井资料以及待评价井附近区域地震资料,所述随钻电阻率曲线测井资料包括随钻电阻率测井曲线,所述辅助测井曲线资料包括辅助测井曲线,所述邻井或导眼井的测井资料包括邻井或导眼井测井曲线;
目的层确定模块,应用根据所述随钻电阻率曲线测井资料、辅助测井曲线资料和邻井或导眼井的测井资料确定水平井段/大斜度井段钻探的目的层、随钻电阻率原始响应曲线;
模型建立模块,用于根据所述辅助测井曲线资料、邻井或导眼井的测井资料对所述目的层进行分层建立初始地层模型并确定初始地层模型的电阻率初值;
模拟响应曲线生成模块,用于根据所述初始地层模型和电阻率初值,应用快速正演法生成对应层的随钻电阻率模拟响应曲线;
判断模块,用于判断所述对应层的随钻电阻率模拟响应曲线与对应层的随钻电阻率原始响应曲线是否一致;
修正模块,用于根据所述对应层的随钻电阻率原始响应曲线、对应层的地质特征通过交互式反演修改所述电阻率初值;
输出模块,用于输出地层真电阻率值;
其中,所述的修正模块包括:
地层电阻率范围确定单元,用于根据对应层所在井的邻井测井资料以及本井的录井资料、气测资料、取心资料判断所述对应层的流体性质以及所述流体性质下的地层电阻率范围;
倾角范围确定单元,用于根据所述对应层的倾角测井资料或地震资料判断对应层的倾向及倾角范围;
径向变化判断单元,用于根据区域内地层对比结合地区沉积相特征,判断对应层砂体厚度在径向上的变化;
修改单元,用于根据所述对应层的随钻电阻率模拟响应曲线与对应层的随钻电阻率原始响应曲线间的相关程度,并结合所述地层电阻率范围、倾向、倾角范围以及砂体厚度在径向上的变化,修改电阻率初值、地层倾向、地层倾角和地层边界。
10.如权利要求9所述的水平井/大斜度井确定地层真电阻率的装置,其特征在于,
所述的辅助测井曲线资料包括:自然伽马测井曲线、井径测井曲线、密度测井曲线、自然电位测井曲线、声波测井曲线、中子测井曲线、微电阻率扫描成像测井资料、方位伽马成像资料或者方位密度成像资料;
所述的邻井或导眼井测井资料还包括:自然伽马测井曲线、电阻率测井曲线。
11.如权利要求9所述的确水平井/大斜度井定地层真电阻率的装置,其特征在于,
所述的电阻率初值包括:水平电阻率和垂直电阻率。
12.如权利要求10所述的水平井/大斜度井确定地层真电阻率的装置,其特征在于,所述的目的层确定模块包括:
响应特征生成单元,用于根据井眼轨迹将随钻电阻率测井曲线和辅助测井曲线投影到垂直深度,确定随钻电阻率测井曲线和辅助测井曲线在直井模式时的响应特征;
目的层确定单元,用于根据随钻电阻率测井曲线和辅助测井曲线在直井模式时的响应特征与邻井或导眼井的测井曲线响应特征对比,确定水平井段/大斜度井段钻探的目的层。
13.如权利要求11或12所述的水平井/大斜度井确定地层真电阻率的装置,其特征在于,所述的模型建立模块包括:
层界面确定单元,用于根据邻井或导眼井测井的曲线特征依据邻井或导眼井测井曲线的斜率变化确定层界面;
分层单元,用于根据所述辅助测井曲线资料和确定的层界面,将所述随钻电阻率测井曲线特征不明显的泥岩夹层、致密层、垮塌层分出来;
电阻率确定单元,用于根据所述邻井或导眼井测井资料和确定的层界面,确定初始地层模型的电阻率初值;
初始地层倾角值确定单元,用于根据所述待评价井附近区域地震资料、微电阻率扫描成像测井资料或邻井测井资料、确定的层界面,确定目的层段各对应层的初始地层倾角值。
14.如权利要求9所述的水平井/大斜度井确定地层真电阻率的装置,其特征在于,所述的模拟响应曲线生成模块,根据所述的电阻率初值和初始地层模型,结合随钻电阻率测井仪器的仪器模型参数,通过快速正演法生成对应层的随钻电阻率模拟响应曲线;
其中,所述随钻电阻率测井仪器的仪器模型参数包括:仪器长度、仪器半径、各阵列单元线圈距、线圈匝数、工作频率。
15.如权利要求9所述的水平井/大斜度井确定地层真电阻率的装置,其特征在于,所述的判断模块判断所述对应层的随钻电阻率模拟响应曲线与对应层的随钻电阻率原始响应曲线是否一致包括:
设定对应层的随钻电阻率模拟响应曲线与对应层的所述随钻电阻率原始响应曲线间的相关系数为S:
N-目的层段内的测量点数;
X-所述对应层的随钻电阻率原始响应曲线目的层段内的数据集;
Y-所述对应层的随钻电阻率模拟响应曲线目的层段内的数据集;
当所述相关系数S的计算结果大于预设阈值,判断对应层的随钻电阻率模拟响应曲线与对应层的所述随钻电阻率原始响应曲线一致。
16.如权利要求15所述的水平井/大斜度井确定地层真电阻率的装置,其特征在于,所述的预设阈值为0.9。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310308564.7A CN103410504B (zh) | 2013-07-22 | 2013-07-22 | 一种水平井/大斜度井确定地层真电阻率的方法及装置 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201310308564.7A CN103410504B (zh) | 2013-07-22 | 2013-07-22 | 一种水平井/大斜度井确定地层真电阻率的方法及装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN103410504A CN103410504A (zh) | 2013-11-27 |
CN103410504B true CN103410504B (zh) | 2016-08-10 |
Family
ID=49603538
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201310308564.7A Active CN103410504B (zh) | 2013-07-22 | 2013-07-22 | 一种水平井/大斜度井确定地层真电阻率的方法及装置 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN103410504B (zh) |
Families Citing this family (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103883318B (zh) * | 2014-03-05 | 2017-02-01 | 中国石油天然气股份有限公司 | 一种中低矿化度地层水储层的参数反演方法及装置 |
CN104102814B (zh) * | 2014-06-11 | 2017-07-11 | 中国科学院地质与地球物理研究所 | 一种基于大地电磁数据反演电阻率和磁化率的方法及系统 |
CN105317431B (zh) * | 2014-07-30 | 2018-09-25 | 中国石油化工股份有限公司 | 一种用于解释和评价水平井测井参数的方法 |
CN105735977A (zh) * | 2014-12-12 | 2016-07-06 | 中石化胜利石油工程有限公司钻井工艺研究院 | 一种随钻测井曲线实时对比方法 |
CN104775811B (zh) * | 2015-04-13 | 2017-09-01 | 中国海洋石油总公司 | 一种地层各向异性信息的提取及校正方法和系统 |
CN104895561A (zh) * | 2015-04-28 | 2015-09-09 | 中国海洋石油总公司 | 基于并行计算技术的随钻电阻率测井联合反演方法 |
CN105114063B (zh) * | 2015-08-04 | 2017-09-29 | 中国石油大学(华东) | 一种斜井双侧向测井曲线校正方法 |
CN105507891B (zh) * | 2015-11-30 | 2018-06-01 | 中国石油天然气股份有限公司 | 获取大斜度井地层电阻率各向异性系数的方法及装置 |
CN106204725B (zh) * | 2016-07-05 | 2019-06-11 | 中国石油天然气股份有限公司 | 一种水平井轨迹校正方法及装置 |
CN107762490A (zh) * | 2016-08-17 | 2018-03-06 | 中国石油化工股份有限公司 | 一种水平井基于双侧向测井真电阻率反演方法 |
CN106869915B (zh) * | 2017-03-14 | 2020-05-08 | 中国石油天然气股份有限公司 | 一种水平井井间隔夹层预测方法及装置 |
CN110552689A (zh) * | 2018-05-15 | 2019-12-10 | 中国石油化工股份有限公司 | 一种确定随钻仪器到地层边界距离的方法 |
CN114542055A (zh) * | 2020-11-24 | 2022-05-27 | 中国石油天然气股份有限公司 | 斜井或倾斜地层的阵列感应测井数据聚焦处理方法及装置 |
CN113779812B (zh) * | 2021-09-28 | 2024-09-27 | 中海油田服务股份有限公司 | 一种大斜度井环境中利用随钻测井数据的薄夹层识别方法 |
CN116398126B (zh) * | 2023-06-07 | 2023-08-25 | 山东万洋石油科技有限公司 | 基于套管开窗水平井小直径随钻电阻率的数据处理方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102140889A (zh) * | 2011-03-09 | 2011-08-03 | 邱世军 | 一种水平井随钻测井地质导向方法 |
CN102162355A (zh) * | 2011-03-09 | 2011-08-24 | 邱世军 | 一种水平井着陆时地质导向方法 |
CN102322256A (zh) * | 2011-06-02 | 2012-01-18 | 中国石油天然气股份有限公司 | 水平井着陆导向方法 |
CN102789003A (zh) * | 2012-08-15 | 2012-11-21 | 中国石油天然气股份有限公司 | 利用阵列感应测井资料进行地层参数反演的方法及装置 |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7191850B2 (en) * | 2004-10-28 | 2007-03-20 | Williams Danny T | Formation dip geo-steering method |
US8473212B2 (en) * | 2009-06-04 | 2013-06-25 | Schlumberger Technology Corporation | Log processing in highly deviated wellbores |
-
2013
- 2013-07-22 CN CN201310308564.7A patent/CN103410504B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102140889A (zh) * | 2011-03-09 | 2011-08-03 | 邱世军 | 一种水平井随钻测井地质导向方法 |
CN102162355A (zh) * | 2011-03-09 | 2011-08-24 | 邱世军 | 一种水平井着陆时地质导向方法 |
CN102322256A (zh) * | 2011-06-02 | 2012-01-18 | 中国石油天然气股份有限公司 | 水平井着陆导向方法 |
CN102789003A (zh) * | 2012-08-15 | 2012-11-21 | 中国石油天然气股份有限公司 | 利用阵列感应测井资料进行地层参数反演的方法及装置 |
Non-Patent Citations (1)
Title |
---|
水平井测井解释技术综述;周灿灿等;《地球物理学进展》;20060331;第21卷(第1期);152-160 * |
Also Published As
Publication number | Publication date |
---|---|
CN103410504A (zh) | 2013-11-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103410504B (zh) | 一种水平井/大斜度井确定地层真电阻率的方法及装置 | |
CN106054248B (zh) | 一种基于大面积致密储层地震岩石物理反演方法 | |
CN103573250B (zh) | 一种计算水平井井眼到地层上下界面距离的方法 | |
CN104502971B (zh) | 页岩储层总有机碳含量和脆性空间展布的检测方法 | |
CN103135135B (zh) | 一种基于疏松砂岩模型进行烃类定量预测的方法和装置 | |
AU2008292169B2 (en) | Monitoring a region of interest in a subsurface formation | |
CN102636812B (zh) | 一种获得碳酸盐岩储层储集空间体积的方法 | |
CN105182424B (zh) | 一种基于斑块饱和模型定量预测储层孔隙度的方法和装置 | |
CN105653815A (zh) | 一种基于岩石物理模型理论的定量解释油藏流体分布方法 | |
CN105510993A (zh) | 前陆盆地深埋挤压型复杂膏盐岩层识别和分布预测方法 | |
CN103733089B (zh) | 用于包括不确定性估计的地下表征的系统和方法 | |
CN105629325B (zh) | 前陆盆地冲积扇精细刻画与预测方法 | |
CN106096081B (zh) | 缝洞型底水油藏的储量估算方法 | |
CN106154322A (zh) | 测井曲线校正方法和装置 | |
CN107942378A (zh) | 一种河流相低含砂率储层预测方法 | |
CN104834007B (zh) | 地震反演过程中计算碳酸盐岩缝洞型储层充填程度的方法 | |
CN107526108B (zh) | 碳酸盐岩缝洞体积校正方法 | |
CN108665545A (zh) | 录井参数三维地质模型建立方法 | |
CN105464650A (zh) | 一种随钻测井解释方法 | |
Singleton | Geophysical data processing, rock property inversion, and geomechanical model building in a Midland Basin development project, Midland/Ector counties, Texas | |
CN110119587B (zh) | 一种地层分层方法及装置 | |
CN107945271A (zh) | 基于地质块体追踪的三维压力场建模方法 | |
CN113885096A (zh) | 高频层序划分与小层对比方法、装置、电子设备及介质 | |
CN108072901A (zh) | 一种获得准确的静水压力和上覆地层压力的方法及系统 | |
SZABÓ | Pore Pressure Prediction In Pannonian Hydrocarbon Reservoir Systems Using An Integrated Interpretation Approach |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant |