CN104636592B - 通过常规测井获取油层指示标志的方法 - Google Patents
通过常规测井获取油层指示标志的方法 Download PDFInfo
- Publication number
- CN104636592B CN104636592B CN201410776712.2A CN201410776712A CN104636592B CN 104636592 B CN104636592 B CN 104636592B CN 201410776712 A CN201410776712 A CN 201410776712A CN 104636592 B CN104636592 B CN 104636592B
- Authority
- CN
- China
- Prior art keywords
- oil reservoir
- warning mark
- formation water
- conventional logging
- resistivity
- 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.)
- Expired - Fee Related
Links
Landscapes
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种通过常规测井获取油层指示标志的方法,包括如下步骤:步骤一、从数据文件中读取常规测井数据;步骤二、测试并计算本地区地层参数;步骤三、根据测井曲线计算孔隙度φ;步骤四、计算视地层水电阻率;步骤五、计算深度点的油层指示标志;步骤六、一个深度点计算完成后,回到步骤一,直到整个深度段计算完成。本发明可有效划分出可产油层段、准确性高、显著降低解释员的工作强度,可以广泛应用于石油勘探与开发领域。
Description
技术领域
本发明涉及石油勘探与开发,特别是涉及一种通过常规测井获取油层指示标志的方法。
背景技术
地层电阻率是反映储层含油或含气的参数之一,但地层水矿化度纵向或横向的变化会引起电阻率的变化,干扰对油气水层的判断:泥质含量的影响,使得油层处的电阻率降低,干扰对油气水层的判断:井眼扩径使得测井所求孔隙度比真实孔隙度高,会将水层误判为油层。
发明内容
本发明的目的是为了克服上述背景技术的不足,提供一种通过常规测井获取油层指示标志的方法,可有效划分出可产油层段、准确性高、显著降低解释员的工作强度。
本发明提供的一种通过常规测井获取油层指示标志的方法,包括如下步骤:步骤一、从数据文件中读取常规测井数据;步骤二、测试并计算本地区地层参数;步骤三、根据测井曲线计算孔隙度φ;步骤四、计算视地层水电阻率;步骤五、计算深度点的油层指示标志;步骤六、一个深度点计算完成后,回到步骤一,直到整个深度段计算完成。
在上述技术方案中,所述步骤一中,所述常规测井数据包括:深度、声波时差、井径、中子、密度、自然伽马、深探测电阻率、浅探测电阻率和自然电位。
在上述技术方案中,所述步骤二中,所述本地区地层参数包括:阿尔奇胶结指数m,阿尔奇常数a,地层水电阻率Rw,视地层水电阻率最小值Rwmn,视地层水电阻率最大值Rwmx,自然伽马测井值最大值GRmax,自然伽马测井值最小值GRmin,钻头直径BITS。
在上述技术方案中,测试并计算地层水电阻率Rw时,需要根据每个深度点进行温度校正。
在上述技术方案中,所述步骤三中,所述测井曲线包括声波时差、中子和密度。
在上述技术方案中,所述步骤四中,所述视地层水电阻率:式中:Rwa-视地层水电阻率,RD-深探测电阻率,φ-地层孔隙度,m-阿尔奇胶结指数,a-阿尔奇常数。
在上述技术方案中,所述步骤五中,所述深度点的油层指示标志:式中:OIL-油层指示标志,Rwa-视地层水电阻率,Rwmn-视地层水电阻率最小值,Rwmx-视地层水电阻率最大值,GR-自然伽马测井值,GRmin-自然伽马测井值最小值,GRmax-自然伽马测井值最大值,CAL-井径,BITS-钻头直径。
本发明通过常规测井获取油层指示标志的方法,具有以下有益效果:在精细测井解释之前,用本发明计算的油层指示标志,可有效划分出可产油层段,将解释员从繁重的劳动中解放出来。目前的试验结果表明,凡是有指示的地方,都有油产出,凡是没有指示的地方,试油要么出水,要么为干层。因此采用本发明,大幅降低对非产层试油的花费,并使解释人员专注于对有油层指示层段的分析。
附图说明
图1为本发明通过常规测井获取油层指示标志的方法的流程示意图;
图2为本发明通过常规测井获取油层指示标志的方法中实施例雅布赖地区雅探L井试油的产层段示意图。
具体实施方式
下面结合附图及实施例对本发明作进一步的详细描述,但该实施例不应理解为对本发明的限制。
参见图1,本发明通过常规测井获取油层指示标志的方法,包括如下步骤:
步骤一、从数据文件中读取常规测井数据,所述常规测井数据包括:深度(m)、声波时差(μs/m)、井径(in或cm,in即英寸)、中子(%)、密度(g/cm3)、自然伽马(API)、深探测电阻率(Ω·m)、浅探测电阻率(Ω·m)和自然电位(mV);
步骤二、测试并计算本地区地层参数,所述本地区地层参数包括:
阿尔奇胶结指数m,
阿尔奇常数a,
地层水电阻率Rw,Rw需要根据每个深度点进行温度校正,(Ω·m),
视地层水电阻率最小值Rwmn,(Ω·m),
视地层水电阻率最大值Rwmx,(Ω·m),
自然伽马测井值最大值GRmax,(API),
自然伽马测井值最小值GRmin,(API),
钻头直径BITS,(in或cm);
步骤三、通过相关软件,根据声波时差、中子和密度这些测井曲线计算孔隙度φ;
步骤四、计算视地层水电阻率:
式中:Rwa-视地层水电阻率,(Ω·m),
RD-深探测电阻率,(Ω·m),
φ-地层孔隙度,小数,
m-阿尔奇胶结指数,
a-阿尔奇常数;
步骤五、计算深度点的油层指示标志:
式中:OIL-油层指示标志,
Rwa-视地层水电阻率,(Ω·m),
Rwmn-视地层水电阻率最小值,(Ω·m),
Rwmx-视地层水电阻率最大值,(Ω·m),
GR-自然伽马测井值,(API),
GRmin-自然伽马测井值最小值,(API),
GRmax-自然伽马测井值最大值,(API),
CAL-井径,(in或cm),
BITS-钻头直径,(in或cm);
步骤六、一个深度点计算完成后,回到步骤一,直到整个深度段计算完成。
本发明通过常规测井获取油层指示标志的方法的工作原理如下:
(1)视地层水电阻率的具体含义
在地层水矿化度不变的地区,可将该地区视为标准水层,可应用包括阿尔奇公式反求地层水电阻率:
式中:Rw-地层孔隙中水的电阻率,(Ω·m),
R0-100%含水地层的电阻率,(Ω·m),
φ-地层孔隙度,小数,
m-阿尔奇胶结指数,
a-阿尔奇常数。
如果找不到标准水层,但仍按这一公式计算地层孔隙空间流体的电阻率,称之为视地层水电阻率:
式中:Rwa-视地层水电阻率,(Ω·m),
RD-深探测电阻率,(Ω·m),
φ-地层孔隙度,小数,
m-阿尔奇胶结指数,
a-阿尔奇常数;
(2)影响视地层水电阻率的因素
矿化度变化:在解释层段内,逐点计算出一条视地层水电阻率曲线,取各层平均值的最小值作为Rw。此外,在解释层段内,如果地层水矿化度是变化的,则应考虑各层矿化度的变化。
泥质:在油层处,往往Rwa较高,可以用Rwa来定性指示其含油性,当然,如果存在泥质的影响,油层的Rwa也可能并不高;
井眼扩径:Rwa还取决于含水孔隙度,在井眼条件差的地方,孔隙度测井受井径扩径或井壁不光滑的影响将使孔隙度增加,增加的多少往往与井径有一定的关系,使得在扩径的地方,本来不是油层的,也求得较高的Rwa,从而引起误判。
(3)测井估算油层指示标志
在考虑泥质、井径与地层水矿化度变化情况下的油层指示标志如下:
式中:OIL-油层指示标志,
Rwa-视地层水电阻率,(Ω·m),
Rwmn-视地层水电阻率最小值(理想情况下作为实际地层水电阻率),(Ω·m),
Rwmx-视地层水电阻率最大值(理想情况下看做实际油层视地层水电阻率最大值),(Ω·m),
GR-自然伽马测井值,(API),
GRmin-自然伽马测井值最小值,(API),
GRmax-自然伽马测井值最大值,(API),
CAL-井径,(in或cm),
BITS-钻头直径,(in或cm)。
油层指示标志实例
作为一种新的油水层识别方法,通过常规测井获取油层指示标志的方法在有显示的层段,都能与录井结果相对应,下面以雅布赖地区试油的产层段为实施例加以阐述:
图2为雅布赖地区雅探L井2340-2680m井段处理显示结果,其中,2648~2655m段油层指示标志OIL大于0,指示为油层。表1为雅探L井2345-2665m井段采样间隔为1m的数据表。
表1雅探L井2345-2665m井段数据表(采样间隔1m)
显然,本领域的技术人员可以对本发明进行各种改动和变型而不脱离本发明的精神和范围。这样,倘若本发明的这些修改和变型属于本发明权利要求及其等同技术的范围之内,则本发明也意图包含这些改动和变型在内。
本说明书中未作详细描述的内容属于本领域专业技术人员公知的现有技术。
Claims (6)
1.一种通过常规测井获取油层指示标志的方法,其特征在于:包括如下步骤:
步骤一、从数据文件中读取常规测井数据;
步骤二、测试并计算本地区地层参数;
步骤三、根据测井曲线计算孔隙度φ;
步骤四、计算视地层水电阻率;
步骤五、计算深度点的油层指示标志,所述步骤五中,所述深度点的油层指示标志:
式中:OIL-油层指示标志,
Rwa-视地层水电阻率,
Rwmn-视地层水电阻率最小值,
Rwmx-视地层水电阻率最大值,
GR-自然伽马测井值,
GRmin-自然伽马测井值最小值,
GRmax-自然伽马测井值最大值,
CAL-井径,
BITS-钻头直径;
步骤六、一个深度点计算完成后,回到步骤一,直到整个深度段计算完成。
2.根据权利要求1所述的通过常规测井获取油层指示标志的方法,其特征在于:所述步骤一中,所述常规测井数据包括:深度、声波时差、井径、中子、密度、自然伽马、深探测电阻率、浅探测电阻率和自然电位。
3.根据权利要求2所述的通过常规测井获取油层指示标志的方法,其特征在于:所述步骤二中,所述本地区地层参数包括:阿尔奇胶结指数m,阿尔奇常数a,地层水电阻率Rw,视地层水电阻率最小值Rwmn,视地层水电阻率最大值Rwmx,自然伽马测井值最大值GRmax,自然伽马测井值最小值GRmin,钻头直径BITS。
4.根据权利要求3所述的通过常规测井获取油层指示标志的方法,其特征在于:测试并计算地层水电阻率Rw时,需要根据每个深度点进行温度校正。
5.根据权利要求3或4所述的通过常规测井获取油层指示标志的方法,其特征在于:所述步骤三中,所述测井曲线包括声波时差、中子和密度。
6.根据权利要求5所述的通过常规测井获取油层指示标志的方法,其特征在于:所述步骤四中,所述视地层水电阻率:
式中:Rwa-视地层水电阻率,
RD-深探测电阻率,
φ-地层孔隙度,
m-阿尔奇胶结指数,
a-阿尔奇常数。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410776712.2A CN104636592B (zh) | 2014-12-15 | 2014-12-15 | 通过常规测井获取油层指示标志的方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201410776712.2A CN104636592B (zh) | 2014-12-15 | 2014-12-15 | 通过常规测井获取油层指示标志的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN104636592A CN104636592A (zh) | 2015-05-20 |
CN104636592B true CN104636592B (zh) | 2017-08-15 |
Family
ID=53215333
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201410776712.2A Expired - Fee Related CN104636592B (zh) | 2014-12-15 | 2014-12-15 | 通过常规测井获取油层指示标志的方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN104636592B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107701177A (zh) * | 2017-08-25 | 2018-02-16 | 中国石油天然气股份有限公司 | 一种地质勘探方法和装置 |
CN109707378B (zh) * | 2019-02-20 | 2021-08-10 | 西北大学 | 一种基于泥浆侵入特性及纵向对比的低阻油层识别方法 |
CN112668814B (zh) * | 2019-09-30 | 2024-03-08 | 中国石油天然气股份有限公司 | 胶结等级确定方法及装置 |
US20210373190A1 (en) * | 2020-06-02 | 2021-12-02 | Halliburton Energy Services, Inc. | Evaluation and visualization of azimuthal resistivity data |
Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102518428A (zh) * | 2011-11-11 | 2012-06-27 | 中国石油天然气股份有限公司 | 油水层识别方法及装置 |
CN103437760A (zh) * | 2013-08-30 | 2013-12-11 | 中国石油天然气股份有限公司 | 一种使用阵列感应数据快速评价油水层的方法 |
CN103867197A (zh) * | 2014-04-04 | 2014-06-18 | 中国石油集团川庆钻探工程有限公司 | 复杂岩性天然气层声波时差判别法 |
-
2014
- 2014-12-15 CN CN201410776712.2A patent/CN104636592B/zh not_active Expired - Fee Related
Patent Citations (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102518428A (zh) * | 2011-11-11 | 2012-06-27 | 中国石油天然气股份有限公司 | 油水层识别方法及装置 |
CN103437760A (zh) * | 2013-08-30 | 2013-12-11 | 中国石油天然气股份有限公司 | 一种使用阵列感应数据快速评价油水层的方法 |
CN103867197A (zh) * | 2014-04-04 | 2014-06-18 | 中国石油集团川庆钻探工程有限公司 | 复杂岩性天然气层声波时差判别法 |
Non-Patent Citations (2)
Title |
---|
地堑式陷盆地储层流体综合识别-以伊通盆地岔路河断陷为例;许圣传等;《吉林大学学报(地球科学版)》;20110331;第41卷(第2期);359-364 * |
视地层水电阻率正态分布分析及应用;张豆娟等;《断块油气田》;20090131;第16卷;123页 * |
Also Published As
Publication number | Publication date |
---|---|
CN104636592A (zh) | 2015-05-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US9091781B2 (en) | Method for estimating formation permeability using time lapse measurements | |
CN104636592B (zh) | 通过常规测井获取油层指示标志的方法 | |
CA2762261C (en) | Process for determining mobile water saturation in a reservoir formation | |
GB2574132A (en) | Chemically tagged drilling fluid additives | |
CN103437760B (zh) | 一种使用阵列感应数据快速评价油水层的方法 | |
AU2011302598B2 (en) | System and method for sweet zone identification in shale gas reservoirs | |
CN107832574B (zh) | 一种基于随钻测井的水平井水淹层解释方法 | |
Kutasov et al. | Determination of formation temperature from bottom-hole temperature logs—a generalized Horner method | |
CN103995297A (zh) | 一种裂缝性碳酸盐岩地层电阻率校正的方法 | |
US11008835B2 (en) | Horizontal reservoir description systems | |
CN107861162B (zh) | 基于微电极测井资料的天然裂缝识别方法及*** | |
CN109711090B (zh) | 一种环空流体综合摩阻系数确定方法及装置 | |
CN103883320B (zh) | 随钻时间推移测井方法 | |
CA2977216C (en) | Surface excitation ranging methods and systems employing a ground well and a supplemental grounding arrangement | |
GB2560840A (en) | Multiple depth of investigation nuclear magnetic resonance logging for determining the porosity and pore type of subterranean formations | |
US11249218B2 (en) | Horizontal reservoir description systems | |
EP3329305B1 (en) | Apparatus and method for determining earth fluid formation | |
Akram et al. | A model to predict wireline formation tester sample contamination | |
CN104699975B (zh) | 从声电效应井下探测器测量数据中提取参数的方法 | |
Vasvari | On the applicability of Dual Laterolog for the deter-mination of fracture parameters in hard rock aquifers | |
CN104879126A (zh) | 一种基于阵列感应测井的碎屑岩储层流体识别方法 | |
CN104018832A (zh) | 一种地层水的电阻率的测量方法 | |
US20120199345A1 (en) | Unconventional Gas Fracture Logging Method and Apparatus | |
Dubey et al. | Aquifer Parameterization in an Alluvial Area: Varanasi District, Uttar Pradesh, India-A Case Study | |
CN108021726A (zh) | 储层有效性的评价方法和装置 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20170815 Termination date: 20171215 |
|
CF01 | Termination of patent right due to non-payment of annual fee |