CN110043254B - 一种基于电缆地层测试资料地层有效渗透率的获取方法 - Google Patents

一种基于电缆地层测试资料地层有效渗透率的获取方法 Download PDF

Info

Publication number
CN110043254B
CN110043254B CN201910391039.3A CN201910391039A CN110043254B CN 110043254 B CN110043254 B CN 110043254B CN 201910391039 A CN201910391039 A CN 201910391039A CN 110043254 B CN110043254 B CN 110043254B
Authority
CN
China
Prior art keywords
fluidity
reservoir
radial
cable
curve
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
Application number
CN201910391039.3A
Other languages
English (en)
Other versions
CN110043254A (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.)
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Zhanjiang Branch
Original Assignee
China National Offshore Oil Corp CNOOC
CNOOC China Ltd Zhanjiang Branch
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 China National Offshore Oil Corp CNOOC, CNOOC China Ltd Zhanjiang Branch filed Critical China National Offshore Oil Corp CNOOC
Priority to CN201910391039.3A priority Critical patent/CN110043254B/zh
Publication of CN110043254A publication Critical patent/CN110043254A/zh
Application granted granted Critical
Publication of CN110043254B publication Critical patent/CN110043254B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B47/00Survey of boreholes or wells
    • EFIXED CONSTRUCTIONS
    • E21EARTH DRILLING; MINING
    • E21BEARTH DRILLING, e.g. DEEP DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B49/00Testing the nature of borehole walls; Formation testing; Methods or apparatus for obtaining samples of soil or well fluids, specially adapted to earth drilling or wells

Abstract

本发明公开了一种基于电缆地层测试资料地层有效渗透率的获取方法,包括如下步骤:S1、分析电缆地层测试压力恢复资料,计算得到径向流流度MR;S2、通过电缆地层测试资料得到的径向流流度MR与常规测井成果孔隙度曲线Φ、泥质含量VCL、储层厚度H之间建立关系模型,利用测井曲线计算静态径向流流度MR‑LOG;S3、将静态流度曲线MR‑LOG与电缆地层测试的流度进行刻度,得到有效的校正后的流度曲线MR‑LOGC;S4、将动态流度曲线MR‑LOGC与DST测试的流度进行尺度转换,从而获得储层的有效渗透率及地层系数,便于更好地分析储层的产能。

Description

一种基于电缆地层测试资料地层有效渗透率的获取方法
技术领域
本发明属于石油天然气勘探开发电缆地层测试技术领域,特别涉及一种基于电缆地层测试资料地层有效渗透率的获取方法。
背景技术
电缆地层测试技术作为唯一测量储层动态(有效)渗透率的测井技术在油气田勘探开发中得到了越来越多的应用。电缆地层测试器通过单井多点采集压力测试数据,可以求取与储层物性和产能相关的渗透率数值,同时用于确定储层的有效厚度、油气水界面和储层间的封堵情况。
而非电缆地层测试的测井技术手段所获得随深度剖面连续的储层参数,主要反映储层的静态特征,却极少直接反映储层的动态特征。实际上,储层有效渗透率的求取方法是通过试油、试采或油藏模拟等方法得到的,利用常规测井资料主要以建立经验关系的方法来实现,其精度、适用性都较差。
早期普遍采用半对数曲线分析法(Horner、MDH)进行试井解释,这就是常规的试井解释方法。当测不到半对数直线段时,或者半对数曲线从向开始难于判断时常规试井解释的应用受到局限。
发明内容
本发明的目的是提供一种基于电缆地层测试压力恢复资料获取地层有效渗透率的获取方法,通过电缆地层测试资料求取径向流流度MR,利用反映储层垂向渗透率差异的常规测井曲线,利用敏感的孔隙度、泥质含量、储层厚度曲线,经刻度计算储层静态径向流流度MR-LOG,通过DST测试流度,实现动态流度MR-LOGC与DST流度之间的尺度转换,从而得到评价储层的有效渗透率的目的。
为了解决上述问题,本发明提供一种基于电缆地层测试资料地层有效渗透率的获取方法,所述方法包括如下步骤:
S1、分析电缆地层测试压力恢复资料,计算得到径向流流度MR
S2、考察储层的各向异性或非均质性,在电缆地层测试流度与DST流度之间需要采用反映储层垂向非均质性的指示曲线,通过电缆地层测试资料得到的径向流流度MR与常规测井曲线建立关系模型,利用测井曲线计算静态径向流流度MR-LOG
S3、将静态径向流流度MR-LOG与电缆地层测试的流度进行刻度,利用电缆地层测试压力恢复资料计算测压点的径向流流度MR,分区块、分油/气模型的公式,得到了经径向流流度MR校正后的动态径向流流度MR-LOGC;其中,建立测压点的流度与对应深度的测井孔隙度、泥质含量的关系。当电缆地层测试器作业时,由于地层钻开后的测井时间、井况、泥浆性能、侵入深度、储层矿物、黏土的敏感性原因,使得电缆地层测试器所测量的结果更接近于储层污染后的有效流度,而测井曲线所反映的信息更接近于储层物性反映的信息,因此刻度时,通过刻度井的测井资料和电缆地层测试资料,分析后进行刻度;
S4、将动态径向流流度MR-LOGC与DST流度进行尺度转换,从而获得储层的有效渗透率及地层系数,所述DST试井解释是根据试井中所测得的资料,包括压力和产量,结合其它资料来判断油气藏类型、测试井类型和井底完善程度,实现确定测试井的特性参数,渗透率、储量、地层温度。
进一步地,所述步骤S2中常规测井成果包括有孔隙度曲线Φ、泥质含量VCL、储层厚度H。
进一步地,所述步骤S2中计算储层静态流度的流度,其计算公式为:
Figure GDA0003838622160000031
式中,MR_LOG为常规测井曲线模型计算得到的静态径向流流度,单位为×10-3μm2/mPa·s;
MR为电缆地层测试资料计算得到的径向流流度,单位为×10-3μm2/mPa·s;Φ为测井解释孔隙度,单位为v/v;
VCL为泥质含量,单位为v/v;
H为测试点所在储层厚度,单位为m;
进一步地,所述步骤S3中所述刻度定义为:
使用本井模型进行刻度:
使用本井模型进行刻度:
当储层段的电缆地层测试样品点达到20个以上、规律明显时,利用本井电缆地层测试压力恢复资料计算的径向流流度MR对由储层孔隙度、泥质含量及储层厚度计算的静态径向流流度MR_LOG进行刻度;
当储层段的电缆地层测试样品点较少,电缆地层测试压力恢复资料计算的径向流流度与孔隙度曲线、泥质含量曲线有明显的相关性时,参照该井其它井段的样品点进行刻度;
当储层段的电缆地层测试样品点计算的流度异常,同该井其它井段的样品点相比,相同物性条件下,流度值偏低或偏高较多时,参照符合该井物性规律的样品点进行本井刻度;
经过以上刻度方法得到反映储层动态渗流能力的动态径向流流度MR_LOGC
使用地区经验模型进行刻度:
当储层段样品点小于20、规律不明显时,使用地区经验模型计算出的静态径向流流度MR-LOG
当储层段样品点丰富,径向流流度MR与利用测井资料结合地区经验模型计算出的静态径向流流度MR-LOG匹配时,使用地区经验模型计算得到的静态径向流流度MR-LOG可直接作为动态径向流流度MR-LOGC使用。
进一步地,所述步骤S4中动态径向流流度MR-LOGC与DST流度的尺度转换方法为:
计算压降曲线或压力恢复曲线得到储层流动系数、地层系数、有效渗透率,
流动系数:
Figure GDA0003838622160000041
地层系数:
Figure GDA0003838622160000042
有效渗透率:
Figure GDA0003838622160000043
式中,K为地层渗透率,单位为×10-3μm2
h为地层厚度,单位为m;
μ为流体粘度,单位为mPa·s;
q为井的地面产量,单位为m3/d;
B为原油体积系数;
MDST为DST测试流度,单位为×10-3μm2/mPa·s,
得到经电缆地层测试资料解释结果刻度的动态径向流流度,计算的动态径向流流度与储层厚度的累计面积即是储层流动系数,其公式为:
Mh=∑MR-LOGC×Δh
式中:Mh为储层流动系数,单位为×10-3μm2·m/mPa·s;
MR-LOGC为动态径向流流度,单位为×10-3μm2/mPa·s;
△h为储层厚度,单位为m,
对DST测试流度,即地层系数与粘度、储层厚度乘积的比值与动态流度进行统计回归分析,得到关系模型:
Figure GDA0003838622160000051
式中,A、B、C为模型系数,通过回归分析得到流体类型及研究区的模型系数。
与现有技术相比,本发明依据电缆地层测试器压力分析方法,求取径向流流度MR,利用反映储层垂向渗透率差异的常规测井曲线,经刻度计算储层静态径向流流度MR-LOG,通过DST测试流度,实现动态径向流流度MR-LOGC与DST流度之间的尺度转换,从而得到评价储层的有效渗透率的目的。
本发明附加技术特征所具有的有益效果将在本说明书具体实施方式部分进行说明。
附图说明
图1是本发明一种基于电缆地层测试资料地层有效渗透率的获取方法的流程图;
图2是本发明实施例的气层MR与模型计算的MR-LOG结果分析图;
图3是本发明实施例的WA区块油层MR与模型计算的MR-LOG结果分析图;
图4是本发明实施例的WB区块油层MR与模型计算的MR-LOG结果分析图;
图5是本发明实施例动态径向流流度曲线刻度图;
图6是本发明实施例DST流度与电缆地层测试流度尺度转化的模型图;
图7是本发明实施例测试层DST流度与测井计算的流度对比图。
具体实施方式
下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而不是全部的实施例。基于本发明中的实施例,本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例,都属于本发明保护的范围。
结合图1所示,一种基于电缆地层测试压力恢复资料地层有效渗透率的获取方法,该方法包括如下步骤:
S1、在流型分析的基础上,根据电缆地层测试资料解释理论及方法计算得到电缆地层测试径向流流度MR
S2、利用测井曲线计算静态径向流流度MR-LOG,通过电缆地层测试资料得到的径向流流度MR与常规测井成果孔隙度曲线Φ、泥质含量VCL、储层厚度H之间建立关系模型,计算储层静态径向流流度。
通过电缆地层测试资料分析处理结果得到的径向流流度MR与常规测井成果孔隙度曲线Φ、泥质含量VCL、储层厚度H之间建立关系模型,如图2是气层径向流流度MR与常规测井曲线所建模型计算数值的关系图,该模型为:
Figure GDA0003838622160000071
式中,MR_LOG为常规测井曲线模型计算得到的静态径向流流度,单位为×10-3μm2/mPa·s;
MR为电缆地层测试资料计算得到的径向流流度,单位为×10-3μm2/mPa·s;
Φ为测井解释孔隙度,单位为v/v;
VCL为泥质含量,单位为v/v;
H为测试点所在储层厚度,单位为m。
结合图3、图4,分别是WA区块、WB区块的油层径向流流度MR与常规测井曲线所建模型计算数值的关系图,模型表达式分别为:
Figure GDA0003838622160000072
Figure GDA0003838622160000073
MR-LOG与储层的孔隙度(Φ)、泥质含量(VCL)、储层厚度(H)建立关系时,将电缆地层测试资料解释的径向流流度MR作为因变量,而孔隙度(Φ)、泥质含量(VCL)、储层厚度(H)为自变量。
S3、动静态流度刻度,利用电缆地层测试资料计算测压点的径向流流度,计算测压点的流度与对应深度的测井孔隙度、泥质含量建立了分区块、分油/气模型的公式,采用建立的地区经验模型:
结合图5所示动态径向流流度曲线刻度方法及刻度模板,图中,MR-LOG是利用地区经验公式计算的静态径向流流度,MR是电缆地层测试压力恢复资料计算的径向流的流度,MR-LOGC是经刻度的动态径向流流度。可见,在未刻度前,电缆地层测试资料计算的径向流流度与测井反映的流度曲线数值差异较大,经储层段电缆地层测试资料计算的径向流流度刻度后测井资料计算的流度精度显著提升。
S4、有效(动态)流度MR-LOGC与DST流度的尺度转换:
DST试井解释就是根据试井中所测得的资料,包括压力和产量,判断油气藏类型、测试井类型和井底完善程度,并确定测试井的特性参数,包括渗透率、储量、地层温度。通过压降曲线或压力恢复曲线储层参数得到储层流动系数、地层系数、有效渗透率。
流动系数:
Figure GDA0003838622160000081
地层系数:
Figure GDA0003838622160000082
有效渗透率:
Figure GDA0003838622160000083
式中,K为地层渗透率,单位为×10-3μm2
h为地层厚度,单位为m;
μ为流体粘度,单位为mPa·s;
q为井的地面产量,单位为m3/d;
B为原油体积系数;
MDST为DST测试流度,单位为×10-3μm2/mPa·s,
得到经电缆地层测试资料解释结果刻度的动态径向流流度,计算的动态径向流流度与储层厚度的累计面积即是储层流动系数,即:
Mh=∑MR-LOGC×Δh
式中:Mh为储层流动系数,单位为×10-3μm2·m/mPa·s;
MR-LOGC为动态径向流流度,单位为×10-3μm2/mPa·s;
△h为储层厚度,单位为m,
研究区域有DST测试层26层/25口井,其中,测试出气的有8层/8口井,测试出油的有18层/17口井。对DST测试流度与动态径向流流度MR_LOGC进行统计回归分析,得到关系模型:
Figure GDA0003838622160000091
式中,A、B、C为模型系数,不同流体类型及研究区的各模型系数及相关系数如表1所示。
表1模型系数及相关系数
Figure GDA0003838622160000092

Claims (1)

1.一种基于电缆地层测试资料地层有效渗透率的获取方法,其特征在于,所述的方法包括如下步骤:
S1、分析电缆地层测试压力恢复资料,计算得到径向流流度MR
S2、考察储层的各向异性或非均质性,在电缆地层测试与DST之间需要采用反映储层垂向非均质性的测井曲线,通过电缆地层测试压力恢复资料得到的径向流流度MR与常规测井成果建立关系模型,利用测井曲线计算静态径向流流度MR-LOG
S3、将静态径向流流度曲线MR-LOG与电缆地层测试的流度进行刻度,利用电缆地层测试压力恢复资料计算测压点的径向流流度MR,分区块、分油/气建立了径向流流度MR与静态径向流流度MR-LOG刻度模型,得到刻度校正后的动态径向流流度曲线MR-LOGC
S4、将动态径向流流度曲线MR-LOGC与DST测试的流度进行尺度转换,从而获得储层的有效渗透率及地层系数;
所述步骤S2中常规测井成果包括有孔隙度曲线Φ、泥质含量VCL、储层厚度H;
所述步骤S2中计算静态径向流流度,其计算公式为:
Figure FDA0003851253990000011
式中,MR_LOG为关系模型计算得到的静态径向流流度,单位为×10-3μm2/mPa·s;
MR为电缆地层测试压力恢复资料计算得到的径向流流度,单位为×10-3μm2/mPa·s;Φ为测井解释孔隙度,单位为v/v;
VCL为泥质含量,单位为v/v;
H为测试点所在储层厚度,单位为m;
所述步骤S3中所述刻度定义为:
使用本井模型进行刻度:
当储层段的电缆地层测试样品点达到20个以上、规律明显时,利用本井电缆地层测试压力恢复资料计算的径向流流度MR对由储层孔隙度、泥质含量及储层厚度计算的静态径向流流度MR_LOG进行刻度;
当储层段的电缆地层测试样品点较少,电缆地层测试压力恢复资料计算的径向流流度与孔隙度曲线、泥质含量曲线有明显的相关性时,参照该井其它井段的样品点进行刻度;
当储层段的电缆地层测试样品点计算的流度异常,同该井其它井段的样品点相比,相同物性条件下,流度值偏低或偏高较多时,参照符合该井物性规律的样品点进行本井刻度;
经过以上刻度方法得到反映储层动态渗流能力的动态径向流流度MR_LOGC
使用地区经验模型进行刻度:
当储层段样品点较少、规律不明显时,使用地区经验模型计算出的静态径向流流度MR-LOG
当储层段样品点大于20个以上,径向流流度MR与利用测井资料结合地区经验模型计算出的静态径向流流度MR-LOG匹配时,使用地区经验模型计算得到的静态径向流流度MR-LOG直接作为动态径向流流度MR-LOGC使用;
所述步骤S4中动态径向流流度MR-LOGC与DST流度的尺度转换方法为:
计算压降曲线或压力恢复曲线得到储层流动系数、地层系数、有效渗透率,
流动系数:
Figure FDA0003851253990000031
地层系数:
Figure FDA0003851253990000032
有效渗透率:
Figure FDA0003851253990000033
式中,K为地层渗透率,单位为×10-3μm2
h为地层厚度,单位为m;
μ为流体粘度,单位为mPa·s;
q为井的地面产量,单位为m3/d;
B为原油体积系数;
MDST为DST测试流度,单位为×10-3μm2/mPa·s;
得到经电缆地层测试压力恢复资料解释结果刻度的动态径向流流度MR_LOGC,计算的动态径向流流度MR_LOGC曲线与储层厚度的累计面积即是储层流动系数,其公式为:
Mh=∑MR-LOGC×Δh
式中:Mh为储层流动系数,单位为×10-3μm2·m/mPa·s;
MR-LOGC为动态径向流流度,单位为×10-3μm2/mPa·s;
△h为储层厚度,单位为m;
对DST测试流度与动态径向流流度进行统计回归分析,得到关系模型:
Figure FDA0003851253990000041
式中,A、B、C为模型系数,通过回归分析得到流体类型及研究区的模型系数。
CN201910391039.3A 2019-05-11 2019-05-11 一种基于电缆地层测试资料地层有效渗透率的获取方法 Active CN110043254B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201910391039.3A CN110043254B (zh) 2019-05-11 2019-05-11 一种基于电缆地层测试资料地层有效渗透率的获取方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201910391039.3A CN110043254B (zh) 2019-05-11 2019-05-11 一种基于电缆地层测试资料地层有效渗透率的获取方法

Publications (2)

Publication Number Publication Date
CN110043254A CN110043254A (zh) 2019-07-23
CN110043254B true CN110043254B (zh) 2022-11-08

Family

ID=67281563

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201910391039.3A Active CN110043254B (zh) 2019-05-11 2019-05-11 一种基于电缆地层测试资料地层有效渗透率的获取方法

Country Status (1)

Country Link
CN (1) CN110043254B (zh)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112147051B (zh) * 2019-12-25 2023-04-25 中国海洋石油集团有限公司 一种基于渗透率分布形态的测压流度标准化方法
CN112145165B (zh) * 2019-12-26 2024-04-23 中国海洋石油集团有限公司 一种微裂缝-孔隙型储层动静态渗透率转换方法
CN112412445B (zh) * 2020-12-18 2023-11-28 陕西延长石油(集团)有限责任公司 一种确定致密气藏混合井型的数字化表征方法
CN115822562B (zh) * 2022-12-28 2023-07-11 中海石油(中国)有限公司海南分公司 一种考虑层内窜流的纵向非均质气藏产能评价方法

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102073772A (zh) * 2011-01-18 2011-05-25 大庆油田有限责任公司 一种mdt压力数值测试模拟方法
CN106600439B (zh) * 2016-12-01 2020-01-10 中国石油大学(北京) 一种低流度油藏物理约束试井解释方法及装置
PE20190137A1 (es) * 2017-06-16 2019-01-22 Oil And Gas Tech Entpr C V Mejorador de la mojabilidad de un aceite pesado
CN109138975B (zh) * 2018-07-10 2022-08-02 中国海洋石油集团有限公司 一种基于时移测井数据的求解相渗特征曲线的新方法

Also Published As

Publication number Publication date
CN110043254A (zh) 2019-07-23

Similar Documents

Publication Publication Date Title
CN110043254B (zh) 一种基于电缆地层测试资料地层有效渗透率的获取方法
CN102953726B (zh) 一种水驱油田优势通道识别方法及装置
CN110321595B (zh) 一种测井提取静态品质系数的断层封闭性评价方法
CN108713089B (zh) 基于钻孔流体和钻探录井估计地层性质
RU2362875C2 (ru) Способ определения давления в подземных пластах
US10480315B2 (en) Average/initial reservoir pressure and wellbore efficiency analysis from rates and downhole pressures
CN108694264B (zh) 一种确定页岩气储层渗透率的方法
CN104899411B (zh) 一种储层产能预测模型建立方法和系统
US20090114009A1 (en) Method for analysis of pressure response in underground formations
CN110162851B (zh) 一种电缆地层测试泵抽数值模拟及其过程的数值校正方法
CN111255442B (zh) 一种利用干扰试井理论评价压裂裂缝方法
CN109138975B (zh) 一种基于时移测井数据的求解相渗特征曲线的新方法
CN112145165B (zh) 一种微裂缝-孔隙型储层动静态渗透率转换方法
CN111963149B (zh) 一种考虑滞地液量增压的压裂后地层压力求取方法
CN110017136B (zh) 一种基于视水层电阻率的水淹层识别与产水率预测方法
US3550445A (en) Method for testing wells for the existence of permeability damage
CN109975189B (zh) 孔隙型砂岩储层产能预测方法与装置
CN110792425B (zh) 一种测量地层流体含水率的方法
Jackson et al. Pressure Measurement and Pressure Gradient Analysis: How Reliable for Determining Fluid Density and Compositional Gradients?
CN109994161B (zh) 趋势基线法结合动态联动法计算地层有机碳含量的方法
CN108756867B (zh) 基于声波测井曲线和电阻率测井曲线进行压裂选层的方法
Proett et al. Low Permeability Interpretation Using a New Wireline Formation Tester" Tight Zone" Pressure Transient Analysis
CN113719271B (zh) 一种试井设计参数矫正方法
CN117077419B (zh) 一种新型缝洞型油藏的地层压力分析方法
CN114492235B (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