DB61∕T 1480-2021 薄油层测井识别与计算规范(陕西省).pdf

ICS75.010CCS E11DB61陕西省地方标准DB 61/T 14802021薄油层测井识别与计算规范Logging identification and evaluation specification for thin oil reservoirs2021-08-12 发布2021-09-12 实施陕西省市场监督管理局发 布DB61/T 14802021I目次前言.II1范围.12规范性引用文件.13术语和定义.14符号与缩略语.15测井识别.46薄油层参数测井计算.77应用.11附录 A（资料性）测井响应围岩影响校正曲线图.12附录 B（资料性）常见岩石的骨架参数及孔隙流体参数.13附录 C（资料性）薄油层测井识别与计算应用实例.14DB61/T 14802021II前言本文件按照GB/T 1.12020标准化工作导则 第1部分：标准话文件的结构和起草规则的规定起草。本文件由陕西延长石油（集团）有限责任公司提出。本文件由陕西省能源局归口。本文件起草单位：陕西延长石油（集团）有限责任公司研究院、西北大学、西安煤航遥感信息有限公司。本文件主要起草人：陈义国、贺永红、马芳侠、王超、方晓君、李广涛、李晓路、许璟、王辉、段亮。本文件由陕西延长石油（集团）有限责任公司研究院负责解释。本文件首次发布。联系信息如下：单位：陕西延长石油（集团）有限责任公司研究院电话：[phone]地址：陕西省西安市雁塔区唐延路61号邮编：710065本文件的发布机构提请注意，声明符合本文件时，可能涉及到5.1条与一种三角洲前缘薄油藏识别方法相关的专利的使用。本文件的发布机构对该专机构有效性和范围无任何立场。本文件的发布机构对有效性和该专利持有人已向本文件的发布机构保证，他愿意同任何申请人在合理且无歧视的条款和条件下，就专利授权许可进行谈判。该专利持有人的声明已在本文件的发布机构备案，相关信息可以通过以下联系方式获得：专利持有人姓名：陈义国地址：陕西省西安市雁塔区唐延路61号延长石油科研中心请注意除上述专利外，本文件的某些内容仍可能涉及专利，本文件的发布机构不承担识别这些专利的责任。DB61/T 148020211薄油层测井识别与计算规范1范围本文件规定了鄂尔多斯盆地中生界的薄油层测井识别与储层参数测井计算方法。本文件适用于鄂尔多斯盆地中生界的薄油层测井识别与计算。2规范性引用文件下列文件中的内容通过文中的规范性引用而构成本文件必不可少的条款。其中，注日期的引用文件，仅该日期对应的版本适用于本文件。不注日期的引用文件，其最新版本（包括所有的修改单）适用于本文件。DZ/T 0217-2020石油天然气储量估算规范3术语和定义下列术语和定义适用于本文件3.1薄油层thin reservoirs储层有效厚度不大于1.5 m，具有工业开采价值的油层。4符号与缩略语4.1 测井资料处理计算公式中符号名称和计量单位0a、1a、2a 回归模型系数；lC 弱纵向分辨率测井曲线测量值；liC 弱分辨率测井曲线采样深度点i对应的地层“真值”；eChfC对lC的回归曲线；hfC 拟弱分辨率曲线；eChC对lC的回归曲线；hC 强纵向分辨率测井曲线测量值；DB61/T 148020212hiC 强分辨率测井曲线采样深度点i对应的地层“真值”；G 强分辨率曲线响应函数的系数矩阵；ig 仪器响应函数系数；J 弱分辨率曲线响应函数的系数矩阵；ij 仪器响应函数系数；R 相关系数；t 滤波系数；12,.,Tnjjj 弱分辨率曲线响应函数的系数向量；12(21),.,mmng gg 强分辨率曲线响应函数的系数向量；12,.,Tka aa 滤波器的系数向量。4.2 薄油层相关参数测井计算公式中符号名称和计量单位AC 薄层声波时差测井值，单位为微秒每米（s/m）；a 比例系数；B 黏土交换阳离子当量电导率，单位为西门子每米（S/m）；b 岩电系数；CEC 阳离子交换能力，单位为毫摩尔每立方厘米（mmol/cm3）；c 泥质指数；pC 压实校正系数；d 钙质指数；e 声波时差泥质校正系数与体积密度泥质校正系数之差；*F 泥质砂岩地层因素；GR 薄层自然伽马曲线的测井值，单位为API；GCUR 希尔奇指数；minGR 自然伽马曲线在纯砂岩处的测井值，单位为API；maxGR 自然伽马曲线在纯泥岩处的测井值，单位为API；GRI 自然伽马相对值；SPI 自然电位相对值；K 渗透率，单位为毫达西（mD）；DB61/T 148020213m 胶结指数；n 饱和度指数；*n 岩石不含黏土的饱和度指数；vQ 阳离子交换容量，单位为毫摩尔每立方厘米（mmol/cm3）；caR 钙质电阻率，单位为欧姆米（m）；clR 黏土电阻率，单位为欧姆米（m）；mfR 泥浆滤液电阻率，单位为欧姆米（m）；shR 泥质电阻率，单位为欧姆米（m）；tR 薄层电阻率测井值，单位为欧姆米（m）；wR 地层水电阻率，单位为欧姆米（m）；wfR 地层可动水电阻率，单位为欧姆米（m）；wbR 束缚水电阻率，单位为欧姆米（m）；xoR 冲洗带电阻率，单位为欧姆米（m）；SP 薄层自然电位测井值，单位为毫伏（mv）；minSP 自然电位在纯砂岩处的测井值，单位为毫伏（mv）；maxSP 自然电位在纯泥岩处的测井值，单位为毫伏（mv）；moS 可动油饱和度；oS 储层含油饱和度；orS 残余油饱和度；wS 含水饱和度；wiS 束缚水饱和度；wtS 地层总含水饱和度；shV泥质含量；clV 黏土体积，单位为立方厘米（cm3）；caV 钙质体积，单位为立方厘米（cm3）；N 薄层含氢指数；Nma 矿物骨架含氢指数；Nf 储层流体含氢指数；薄层孔隙度；t 总孔隙度；e 有效孔隙度；b 薄层体积密度测井值，单位为克每立方厘米（g/cm3）；ma 骨架密度值，单位为克每立方厘米（g/cm3）；DB61/T 148020214f 流体密度值，单位为克每立方厘米（g/cm3）；G 平均颗粒密度，单位为克每立方厘米（g/cm3）；t 声波时差测井值，单位为微秒每米（s/m）；mat 岩石骨架声波时差值，单位为微秒每米（s/m）；ft 流体声波时差值，单位为微秒每米（s/m）。5测井识别5.1测井资料处理5.1.1测井识别前应按以下步骤对声波时差、自然伽马、电阻率数据进行处理。5.1.1.1测井曲线离散化纵向分辨率强于0.5 m的强纵向分辨率测井曲线hC按式（1）离散化，纵向分辨率弱于1.5 m的弱纵向分辨率测井曲线lC按式（2）离散化：mihiihCgC1.(1)niliilCjC1.(2)5.1.1.2匹配滤波按式（3）对hC滤波处理，得拟弱分辨率曲线hfC：hfhCfC.(3)1fFJ G.(4)5.1.1.3滤波器系数向量计算使hfC与lC匹配相同的纵向分辨率，滤波器的系数向量按式（5）矩阵求取：DB61/T 14802021512121122121.mmnn mmggggjajaggggjag.(5)5.1.1.4强分辨率曲线滤波按式（6）对强分辨率曲线hC滤波处理，得到与lC具有相同纵向分辨率的滤波曲线hfC：thfhCC.(6)5.1.1.5回归建模对测井曲线上每个深度采样点数值，对hfC与lC进行相关分析，根据相关系数R最大化原则，确定一个包含处理点在内的最佳相关区间，建立最佳函数回归模型，以建立二次函数模型（7）为例：2210hfhfeCaCaaC.(7)5.1.1.6强分辨率曲线计算强分辨率曲线回归估算模型为式（8）：2210hheCaCaaC.(8)5.1.1.7误差校正a)对于数据建模回归计算后，得式（9）：leelCCCC.(9)b)按式（10）对分辨率匹配结果校正，以此反演得到薄层测井响应真值：221llCR CRC.(10)c)当hfC与lC相关程度高时，R趋近于 1，C接近于lC；当hfC与lC相关程度较差时，R趋近于 0，C接近于lC。5.1.1.8测井响应围岩影响校正曲线图参见附录 A。5.2薄油层单井识别DB61/T 1480202165.2.1采用声波时差与深感应电阻率数据绘制交会图（见图 1），确定声波时差与深感应电阻率下限值。5.2.2采用声波时差与含油饱和度数据绘制交会图（见图 2），确定含油饱和度下限值。5.2.3采用孔隙度与日产油数据绘制交会图（见图 3），确定孔隙度下限值。5.2.4采用渗透率与日产油数据绘制交会图（见图 4），确定渗透率下限值。图 1声波时差与深感应电阻率交会图图 2声波时差与含油饱和度交会图DB61/T 148020217图 3有效厚度孔隙度下限图图 4有效厚度渗透率下限图5.2.5图版符合率应符合 DZ/T 02172020 第 6.2.2 条要求。6薄油层参数测井计算6.1薄层泥质含量计算6.1.1常规砂泥岩地层的薄层泥质含量计算对常规砂泥岩地层，按式（11）、（12）计算薄层泥质含量：minmaxminGRGRGRGRIGR.(11)DB61/T 1480202181212GCURIGCURshGRV.(12)6.1.2高放射性砂泥岩地层的薄层泥质含量计算对高伽马砂泥岩地层，储层自然伽马普遍高于80API，利用自然电位测井值，按式（13）、（14）计算薄层泥质含量：minmaxminSPSPSPSPISP.(13)1212shGCURIGCURSPV.(14)6.1.3岩心数据回归建模计算薄层泥质含量利用典型薄层岩心资料刻度测井资料，通过多元统计分析，建立地区高精度回归模型式（15）计算：RTACSPGRfV,sh.(15)6.2薄层孔隙度计算6.2.1对不发育裂缝的薄层利用声波时差测井值，按式（16）计算薄层孔隙度：mafmaptttt1C.(16)常见岩石的骨架参数及孔隙流体参数参见附录B。6.2.2对发育裂缝且含天然气的薄层利用密度测井值，按式（17）计算薄层孔隙度：faa-mbm.(17)常见岩石的骨架参数及孔隙流体参数参见附录B。6.2.3对发育裂缝且不含天然气的薄层利用密度测井值，按式（18）计算薄层孔隙度：faa-mbm.(18)常见岩石的骨架参数及孔隙流体参数参见附录B。6.2.4对发育裂缝且不含天然气的薄层DB61/T 148020219利用中子测井值，按式（19）计算薄层孔隙度：NmaaNNNfm.(19)常见岩石的骨架参数及孔隙流体参数参见附录B。6.2.5岩心数据回归建模计算薄层孔隙度利用薄层岩心资料刻度测井资料，通过多元数据相关分析，建立地区高精度回归模型式（20）计算：=fAC Vsh.(20)6.3薄层渗透率计算6.3.1Timur 经验公式对薄层岩心化验分析资料较少时，按式（21）估算薄层渗透率：24.4316.0wiSK.(21)6.3.2岩心数据回归建模计算薄层渗透率利用典型薄层岩心资料刻度测井资料，通过多元统计分析，建立地区高精度回归模型式（22）计算：shVfK,.(22)6.4薄层含水饱和度计算6.4.1泥质含量低于 15%的砂岩薄层含水饱和度计算对于泥质含量低于15%的砂岩薄层，按式（23）计算含水饱和度：nwwmta b RSR.(23)6.4.2泥质含量高于 15%的泥质砂岩薄层含水饱和度计算6.4.2.1泥质均匀分布的泥质砂岩薄层含水饱和度计算黏土均匀分布的泥质砂岩薄层，按式（24）计算含水饱和度：21wwmwclcltSRaSRVR.(24)6.4.2.2泥质层状分布的泥质砂岩薄层含水饱和度计算黏土呈层状分布的泥质砂岩薄层，按式（25）计算含水饱和度：DB61/T 14802021102122wclclwtRe Ve VSR.(25)6.4.2.3富含单一导电性黏土的泥质砂岩薄层含水饱和度计算富含蒙皂石等单一导电性黏土的泥质砂岩薄层，按式（26）、（27）计算含水饱和度：wwwt)(-nw1SSRQBRFRV.(26)tGtCECQ1v.(27)6.4.2.4富含多种导电性黏土的泥质砂岩薄层含水饱和度计算富含蒙皂石、绿泥石、伊利石等多种导电性黏土的泥质砂岩薄层，泥质阳离子交换能力难以准确确定，认为泥质砂岩导电性是黏土束缚水与孔隙自由水并联的结果，按式（28）、（29）计算含水饱和度：2122yyttwfwtRRS.(28)()2wiwbwfwbSRRyR.(29)6.4.2.5复杂岩性薄层含水饱和度计算含多种黏土及钙质的复杂岩性薄层，按式（30）、（31）、（32）计算含水饱和度：twecacashcshwRRaRVRVS2dn-1.(30)21cshV.(31)21dcaV.(32)6.5薄层含油饱和度计算已知含水饱和度时，按式（33）计算薄层含油饱和度：woSS1.(33)6.6薄层残余油饱和度计算按式（34）计算薄层残余油饱和度：DB61/T 1480202111xommfnorRRbaS1.(34)6.7薄层可动油饱和度计算已知含油饱和度及残余油饱和度时，按式（35）计算薄层可动油饱和度：oromoSSS.(35)7应用薄油层测井识别与计算应用实例见附录C。DB61/T 1480202112附录A（资料性）测井响应围岩影响校正曲线图A.1测井响应围岩影响校正曲线见图A.1图 A.1CH46 井测井响应围岩影响校正曲线图DB61/T 1480202113AA附录B（资料性）常见岩石的骨架参数及孔隙流体参数B.1常见岩石的骨架参数及孔隙流体参数见表B.1表 B.1常见岩石的骨架参数及孔隙流体参数类别声波时差体积密度g/cm3补偿中子s/ms/ft砂岩骨架182.055.52.65-0.05淡水钻井液620.0189.01.001.00盐水钻井液608.0185.01.101.00DB61/T 1480202114BB附录C（资料性）薄油层测井识别与计算应用实例C.1薄油层测井识别与计算应用实例见图C.1图 C.1薄油层测井识别与计算应用实例_