CN113236237B - 一种基于常规测井曲线的致密砂岩储层有效性评价方法 - Google Patents
一种基于常规测井曲线的致密砂岩储层有效性评价方法 Download PDFInfo
- Publication number
- CN113236237B CN113236237B CN202110511151.3A CN202110511151A CN113236237B CN 113236237 B CN113236237 B CN 113236237B CN 202110511151 A CN202110511151 A CN 202110511151A CN 113236237 B CN113236237 B CN 113236237B
- Authority
- CN
- China
- Prior art keywords
- curve
- cnl
- effectiveness
- logging
- gas
- 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 37
- 239000004576 sand Substances 0.000 claims abstract description 55
- 238000011156 evaluation Methods 0.000 claims abstract description 19
- 238000010606 normalization Methods 0.000 claims abstract description 7
- 238000004364 calculation method Methods 0.000 claims description 9
- 239000011435 rock Substances 0.000 claims description 7
- 238000011160 research Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 4
- 230000000704 physical effect Effects 0.000 abstract description 7
- 238000011161 development Methods 0.000 abstract description 5
- 238000005516 engineering process Methods 0.000 abstract description 4
- 239000007789 gas Substances 0.000 description 41
- 230000015572 biosynthetic process Effects 0.000 description 10
- 230000000694 effects Effects 0.000 description 8
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 7
- 229910052739 hydrogen Inorganic materials 0.000 description 7
- 239000001257 hydrogen Substances 0.000 description 7
- 239000000243 solution Substances 0.000 description 7
- 238000009412 basement excavation Methods 0.000 description 4
- 238000004458 analytical method Methods 0.000 description 3
- 238000005553 drilling Methods 0.000 description 3
- 235000019738 Limestone Nutrition 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012512 characterization method Methods 0.000 description 2
- 239000004927 clay Substances 0.000 description 2
- 230000007547 defect Effects 0.000 description 2
- 239000006028 limestone Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000002285 radioactive effect Effects 0.000 description 2
- OKTJSMMVPCPJKN-UHFFFAOYSA-N Carbon Chemical compound [C] OKTJSMMVPCPJKN-UHFFFAOYSA-N 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 229910052799 carbon Inorganic materials 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000010586 diagram Methods 0.000 description 1
- 239000008398 formation water Substances 0.000 description 1
- 238000002347 injection Methods 0.000 description 1
- 239000007924 injection Substances 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 239000000203 mixture Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000035699 permeability Effects 0.000 description 1
- 238000003904 radioactive pollution Methods 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- E—FIXED CONSTRUCTIONS
- E21—EARTH OR ROCK DRILLING; MINING
- E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
- E21B49/00—Testing 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
Landscapes
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Geology (AREA)
- Mining & Mineral Resources (AREA)
- Physics & Mathematics (AREA)
- Environmental & Geological Engineering (AREA)
- Fluid Mechanics (AREA)
- General Life Sciences & Earth Sciences (AREA)
- Geochemistry & Mineralogy (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于常规测井曲线的致密砂岩储层有效性评价方法，涉及气藏勘探领域，属于井筒地球物理测井评价技术，解决了现有技术中窄河道砂体发育但物性和含气性较差的无效储层以及高伽马有效储层的难以识别问题。本发明所述的基于常规测井曲线的致密砂岩储层有效性评价方法包括以下步骤：声波测井曲线归一化处理；中子测井曲线归一化处理；电阻率曲线取对数值后归一化处理；判断泥岩和砂体；计算声波、中子归一化后的幅度差曲线；计算声波、电阻率归一化后的幅度差曲线；计算砂体有效性评价指标。本发明在缺乏元素测井等特殊测井资料的情况下，利用常规测井曲线，以测井曲线重叠法为基础，通过计算砂体有效性评价指标能快速直观的评价致密砂岩储层有效性，可为下一步该区勘探开发的各项工作提供有力技术支持和保障。
Description
技术领域
本发明属于井筒地球物理勘探领域，涉及测井资料地质评价技术，具体涉及一种基于常规测井曲线的致密砂岩储层有效性评价方法。
背景技术
金华—中台山区块目前的钻采数据和地质研究成果表明侏罗系沙溪庙组气藏具有良好的勘探开发前景。四川盆地中部侏罗系沙溪庙组将是未来勘探开发的重点层系。四川盆地侏罗系沙溪庙组储层具有低孔隙度、低渗透率的特点，为典型的致密砂岩储层。随着区域钻采数据的增多和研究的不断深入，发现该区部分井段存在测井曲线指示砂体发育，但试油结果均证实为无效储层。测井曲线有含气性指示但物性较差，岩心物性分析数据证实孔隙度低于区域储层下限值，为致密无效储层。另外还发现存在一些井段具有高伽马测井响应特征，但钻井和录井均有良好油气显示，岩心分析物性也较好的储层，这类储层在缺少元素测井时易被错误解释为无效储层。因此该区沙溪庙组致密砂岩储层的有效性测井评价面临巨大挑战。
地层含放射性成分、黏土类型的变化以及黏土颗粒吸附放射性有机分子、邻井注水时造成的放射性污染等都有可能造成储层显示高伽马值特征。因此，建立高伽马值储层的识别方法，进而进行含气性评价尤为重要。目前国内外针对高伽马储层测井评价已形成了以下相关技术，如张涛(2012)采用Pe-GR交会图法识别高伽马储层；程超(2008)采用元素测井资料和岩心分析资料识别高伽马储层；曲娟(2013)以常规测井为基础，首先分别计算出声波孔隙度、密度孔隙度和中子孔隙度曲线，再根据三条孔隙度的差值及比值法来评价高伽马储层。这些高伽马储层识别方法需要特殊测井资料或者需要经过比较复杂的测井处理，不利于现场工程师使用。发明专利(CN202010944199.9)公开了一种致密砂岩低阻储层挖掘效应的定量表征及含气性评价方法。该方法依据气测录井资料计算出地层中挖掘效应指数，实现了挖掘效应现象从定性描述到定量表征。发明专利(CN202010714325.1)公开了一种致密砂岩无水气藏含气性评价方法。该方法充分利用天然气在中子测井曲线上的挖掘效应以进行含气性指标参数的定量表征，通过分类评价有效的解决了现有技术致密砂岩无水气藏在缺乏地层水等资料的情况下难以准确评价储层含气性的技术问题。以上两种方法在金华—中台山区块均有良好的应用，但未对该区高伽马储层的适应性进行论证。发明专利(CN201510047950.4)公开了一种利用纵横波速度资料进行储层含气性分析的方法，该方法通过求取一系列岩石弹性力学参数来进行含气性评价，与张鹏等(2020)基于Xu-White模型提出的一种致密砂岩气层识别新方法一样，使用前提均需要获取地层纵横波资料。ΔlgR技术是Exxon和Esso公司于1979年开发的一种烃源岩测井评价方法。该方法核心是将线性刻度的声波时差曲线与对数刻度的电阻率曲线重叠，确定富含有机质层段与重叠曲线间的对应关系。Passey(1990)等在此基础上提出了考虑烃源岩热变质指数(LOM)、不同成熟度条件下烃源岩的测井评价法，即ΔlgR法，随后该方法在烃源岩评价领域有了广泛的应用(张小莉等，1998；张志伟等，2000)。李延钧等(2013)探索了用ΔlgR法开展页岩气储层有机碳含量评价的可行性，现已成为页岩储层TOC参数计算的重要方法之一。通过广泛调研，该方法目前还尚未应用到常规储层有效性的评价上。
发明内容
本发明的目的在于，提供一种基于常规测井曲线的致密砂岩储层有效性评价方法，以AC、CNL、RT三条原始测井曲线为基础，经归一化后采用两两重叠计算其幅度差，无需经过复杂计算便能快速直观的评价致密砂岩储层有效性。实践证明该方法在金华秋林区块侏罗系致密砂岩储层的有效性评价中效果良好，方便现场工程师应用，不仅能有效解决致密无效储层和高伽马有效储层的识别难题，还能进行有效储层含气性评价。
为实现上述目的，本发明采用的技术方案如下：
本发明所述的一种基于常规测井曲线的致密砂岩储层有效性评价方法，包括以下步骤：
步骤1.声波时差测井曲线(AC)归一化处理，得到AC_01曲线；
步骤2.中子测井曲线(CNL)归一化处理，得到CNL_01曲线；
步骤3.对电阻率曲线(Rt)取对数值，得到logRt曲线；
步骤4.对logRt曲线进行归一化处理，得到logRt_01曲线；
步骤5.利用CNL_01曲线判断出砂体和泥岩；
步骤6.计算CNL_01曲线和AC_01曲线的幅度差DACCNL_01；
步骤7.计算logRt_01曲线和AC_01曲线的幅度差DlogRtAC_01；
步骤8.计算储层有效性指标GAS_FLAG曲线，并利用该曲线对步骤5识别出的砂体进行储层有效性评价。
本发明的部分实施方案中，所述步骤1中AC_01曲线的计算公式为：
式中AC_01为中子测井曲线归一化后的曲线，无量纲；AC为实测中子测井曲线，单位为μs/ft；ACmin为研究井段中子曲线最小值，单位为μs/ft；ACmax为研究井段中子曲线最大值，单位为μs/ft。
本发明的部分实施方案中，所述步骤2中CNL_01曲线的计算公式为：
式中CNL_01为中子测井曲线归一化后的曲线，无量纲；CNL为实测中子测井曲线，单位为％；CNLmin为研究井段中子曲线最小值，单位为％；CNLmax为研究井段中子曲线最大值，单位为％。
本发明的部分实施方案中，所述步骤3中logRt曲线为实测电阻率曲线的对数值，步骤4中logRt_01曲线的计算公式为：
式中logRt_01为实测电阻率测井曲线取对数后归一化后的曲线，无量纲；logRt为实测电阻率测井曲线的对数，单位为Ω·m；logRtmin为研究井段中子曲线最小值，单位为Ω·m；logRtmax为研究井段中子曲线最大值，单位为Ω·m。
本发明的部分实施方案中，所述步骤5中，利用步骤2中求取的CNL_01曲线判断砂体和泥岩；其中，CNL_01曲线≥0.4时为泥岩，CNL_01曲线<0.4时为砂体。CNL为视石灰岩孔隙度曲线，表征的是地层中的含氢指数，理论上归一化后的CNL_01曲线越接近0，则岩石越致密，值越大含氢指数越高。由于泥岩的含氢指数高，因此高CNL_01值表征泥岩，低值则为砂体。纯泥岩的含氢指数理论值介于0.4～0.45，实践表明研究区砂体和泥岩的阈值取0.4应用效果良好。该方法能有效识别高自然伽马砂体，规避了在缺乏元素测井时利用自然伽马曲线无法划分高伽马储层的缺陷。
本发明的部分实施方案中，所述步骤6中的DACCNL_01曲线为步骤1计算得到的AC_01曲线与步骤2计算得到的CNL_01曲线的差值，即DACCNL_01＝AC_01-CNL_01。DACCNL_01能有效表征砂岩的有效孔隙度，DACCNL_01值越小，有效孔隙度越小；DACCNL_01值越大，有效孔隙度越大。
本发明的部分实施方案中，所述步骤7中的DlogRtAC_01曲线为步骤3计算得到的logRt_01曲线与步骤1计算得到的AC_01曲线的差值，即DlogRtAC_01＝logRt_01-AC_01。DlogRtAC_01能有效表征砂体的含气性，砂体DlogRtAC_01值越小含气性越好；砂体DlogRtAC_01值越大，含气性越差。
本发明的部分实施方案中，步骤8中的所述的储层有效性指标GAS_FLAG曲线为步骤6计算得到的DACCNL_01曲线与步骤7计算得到的DlogRtAC_01曲线幅度差，即GAS_FLAG＝DACCNL_01-DlogRtAC_01。
GAS_FLAG能较好的表征砂体的有效性，不仅能有效识别物性和含气性较差的无效储层，还可对有效储层进行含气性评价。由于DACCNL_01越大物性越好，DlogRtAC_01值越小含气性越好，所以两条曲线度在砂岩段重叠时，正幅度差越大储层有效性越好，负幅度差越大为无效储层；两条曲线近似重叠时，即DACCNL_01与DlogRtAC_01之差近似等于0时为储层有效性的分界。因此，在研究区GAS_FLAG曲线小于0值的砂体为无效储层，GAS_FLAG大于0值的砂体为有效储层，GAS_FLAG值越大，含气性越好。
本发明的部分实施方案中，GAS_FLAG值为0.0～0.15时，砂体为含气层；
GAS_FLAG值大于0.15小于等于0.25时，砂体为差气层；
GAS_FLAG值大于0.25时，砂体为气层。
与现有技术相比，本发明具有以下有益效果：
本发明构思巧妙，设计科学，在缺乏元素测井等特殊测井资料的情况下，本发明创造性地基于常规测井曲线，以AC-CNL曲线重叠、AC-RT重叠为基础，通过计算曲线幅度差得到储层有效性评价指标，该指标能快速直观的评价致密砂岩储层有效性。实际应用效果表明本发明方法具有快速直观、准确性高的特点，能为下一步该区勘探开发的各项工作提供有力技术支持和保障。
附图说明
图1是本发明技术流程图；
图2是本发明四川盆地某井沙溪庙组致密砂岩储层有效性评价效果图。
具体实施方式
为了使本发明的目的、技术方案及优点更加清楚明白，以下结合附图，对本发明进一步详细说明。显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。在本发明的描述中，需要说明的是，术语“中心”、“上”、“下”、“左”、“右”、“竖直”、“水平”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本发明和简化描述，而不是指示或暗示所指的装置或元件必须具有特定的方位或者以特定的方位构造和操作，因此其不能理解为对本发明的限制。此外，术语“第一”、“第二”、“第三”仅用于描述目的，而不能理解为指示或暗示相对重要性。在本发明的描述中，需要说明的是，除非另有明确的规定和限定，术语“安装”、“相连”、“连接”应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或一体地连接；当然的，还可以是机械连接，也可以是电连接；另外的，还可以是直接相连，也可以是通过中间媒介间接相连，或者可以是两个元件内说明书部的连通。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本发明中的具体含义。
本发明所述的一种基于常规测井曲线的致密砂岩储层有效性评价方法，包括以下步骤：
步骤1.声波时差测井曲线(AC)归一化处理，得到AC_01曲线；
步骤2.中子测井曲线(CNL)归一化处理，得到CNL_01曲线；
步骤3.对电阻率曲线(Rt)取对数值，得到logRt曲线；
步骤4.对logRt曲线进行归一化处理，得到logRt_01曲线；
步骤5.利用CNL_01曲线判断出砂体和泥岩；
步骤6.计算CNL_01曲线和AC_01曲线的幅度差DACCNL_01；
步骤7.计算logRt_01曲线和AC_01曲线的幅度差DlogRtAC_01；
步骤8.计算储层有效性指标GAS_FLAG曲线，并利用该曲线对步骤5识别出的砂体进行储层有效性评价。
所述声波时差测井曲线(AC)、中子测井曲线(CNL)、电阻率测井曲线(Rt)均指原始的测井曲线。
所述步骤1中AC_01曲线的计算公式为：
式中AC_01为中子测井曲线归一化后的曲线，无量纲；AC为实测中子测井曲线，单位为μs/ft；ACmin为研究井段中子曲线最小值，单位为μs/ft；ACmax为研究井段中子曲线最大值，单位为μs/ft。
所述步骤2中CNL_01曲线的计算公式为：
式中CNL_01为中子测井曲线归一化后的曲线，无量纲；CNL为实测中子测井曲线，单位为％；CNLmin为研究井段中子曲线最小值，单位为为％；CNLmax为研究井段中子曲线最大值，单位为％。
所述步骤3中logRt曲线为实测电阻率曲线(Rt)的对数值，步骤4中logRt_01曲线的计算公式为：
式中logRt_01为实测电阻率测井曲线取对数后归一化后的曲线，无量纲；logRt为实测电阻率测井曲线的对数，单位为Ω·m；logRtmin为研究井段中子曲线最小值，单位为Ω·m；logRtmax为研究井段中子曲线最大值，单位为Ω·m。
步骤5中，利用步骤2中求取的CNL_01曲线判断砂体和泥岩；其中，CNL_01曲线≥0.4时为泥岩，CNL_01曲线<0.4时为砂体。
所述步骤6中的DACCNL_01曲线为步骤1计算得到的AC_01曲线与步骤2计算得到的CNL_01曲线的差值，即DACCNL_01＝AC_01-CNL_01。
所述步骤7中的DlogRtAC_01曲线为步骤3计算得到的logRt_01曲线与步骤1计算得到的AC_01曲线的差值，即DlogRtAC_01＝logRt_01-AC_01。
所述步骤8中的所述的储层有效性指标GAS_FLAG曲线为步骤6计算得到的DACCNL_01曲线与步骤7计算得到的DlogRtAC_01曲线幅度差，即GAS_FLAG＝DACCNL_01-DlogRtAC_01。GAS_FLAG曲线小于0值的砂体为无效储层，GAS_FLAG大于0值的砂体为有效储层，GAS_FLAG值越大，含气性越好。GAS_FLAG值为0.0～0.15时，砂体为含气层；GAS_FLAG值大于0.15小于等于0.25时，砂体为差气层；GAS_FLAG值大于0.25时，砂体为气层。
为了使本领域技术人员能够更好地理解本发明的技术方案，下面四川盆地某井为例对本发明进行详细描述：
(1)声波测井曲线归一化处理，得到AC_01曲线：利用声波测井曲线数据，按公式(1)计算声波时差归一化曲线AC_01。
式中AC_01为中子测井曲线归一化后的曲线，无量纲；AC为实测中子测井曲线，单位为μs/ft；ACmin为研究井段中子曲线最小值，单位为μs/ft；ACmax为研究井段中子曲线最大值，单位为μs/ft。
(2)中子测井曲线归一化处理，得到CNL_01曲线：利用中子测井曲线数据，按公式(2)计算声波时差归一化曲线CNL_01。
式中CNL_01为中子测井曲线归一化后的曲线，无量纲；CNL为实测中子测井曲线，单位为％；CNLmin为研究井段中子曲线最小值，单位为为％；CNLmax为研究井段中子曲线最大值，单位为％。
(3)电阻率曲线取对数值，得到logRt曲线。
(4)对logRt曲线进行归一化处理，得到logRt_01曲线。利用电阻率取对数后的logRt曲线，按公式(3)计算电阻率曲线取对数值后的归一化曲线logRt_01。
式中logRt_01为实测电阻率测井曲线取对数后归一化后的曲线，无量纲；logRt为实测电阻率测井曲线的对数，单位为Ω·m；logRtmin为研究井段中子曲线最小值，单位为Ω·m；logRtmax为研究井段中子曲线最大值，单位为Ω·m。
(5)利用CNL_01曲线判断出砂体和泥岩。CNL为视石灰岩孔隙度曲线，表征的是地层中的含氢指数。理论上归一化后的CNL_01曲线越接近0，则岩石越致密，值越大含氢指数越高。由于泥岩的含氢指数高，因此高CNL_01值表征泥岩，低值则为砂体。研究区砂体和泥岩的阈值取0.4能很好的划分砂体和泥岩。CNL_01曲线≥0.4时为泥岩，CNL_01曲线<0.4时为砂体。该方法能有效识别高自然伽马砂体，规避了在缺乏元素测井时利用自然伽马曲线无法划分高伽马储层的缺陷。
(6)计算CNL_01曲线和AC_01曲线的幅度差DACCNL_01。DACCNL_01曲线为步骤1计算得到的AC_01曲线与步骤2计算得到的CNL_01曲线的差值，即DACCNL_01＝AC_01-CNL_01。其值能有效表征砂岩的有效孔隙度，DACCNL_01值越小，有效孔隙度越小；DACCNL_01值越大，有效孔隙度越大。
(7)计算logRt_01曲线和AC_01曲线的幅度差DlogRtAC_01。DlogRtAC_01曲线为步骤4计算得到的logRt_01曲线与步骤1计算得到的AC_01曲线的差值，即DlogRtAC_01＝logRt_01-AC_01。其值能有效表征砂体的含气性，砂体DlogRtAC_01值越小含气性越好；砂体DlogRtAC_01值越大，含气性越差。
(8)计算储层有效性指标GAS_FLAG曲线，并利用该曲线对步骤5识别出的砂体进行储层有效性评价。储层有效性指标GAS_FLAG曲线为步骤6计算得到的DACCNL_01曲线与步骤7计算得到的DlogRtAC_01曲线幅度差，即GAS_FLAG＝DACCNL_01-DlogRtAC_01。该指标曲线能较好的表征砂体的有效性，不仅能有效识别物性和含气性较差的无效储层，还可对有效储层进行含气性评价。GAS_FLAG曲线小于0值的砂体为无效储层，大于0值的砂体为有效储层，且值越大，含气性越好，GAS_FLAG值为0.0～0.15时，砂体为含气层；GAS_FLAG值大于0.15小于等于0.25时，砂体为差气层；GAS_FLAG值大于0.25时，砂体为气层。
结果如附图2所示，说明步骤8指示效果较好。
实践结果表明本发明的基于常规测井曲线的致密砂岩储层有效性评价方法可以快速准确的评价储层的含气性。
最后应说明的是：以上各实施例仅仅为本发明的较优实施例用以说明本发明的技术方案，而非对其限制，当然更不是限制本发明的专利范围；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围；也就是说，但凡在本发明的主体设计思想和精神上作出的毫无实质意义的改动或润色，其所解决的技术问题仍然与本发明一致的，均应当包含在本发明的保护范围之内；另外，将本发明的技术方案直接或间接的运用在其他相关的技术领域，均同理包括在本发明的专利保护范围。
Claims (7)
1.一种基于常规测井曲线的致密砂岩储层有效性评价方法，其特征在于，包括以下步骤：
步骤1.声波时差测井曲线AC归一化处理，得到AC_01曲线；
步骤2.中子测井曲线CNL归一化处理，得到CNL_01曲线；
步骤3.对电阻率曲线Rt取对数值，得到logRt曲线；
步骤4.对logRt曲线进行归一化处理，得到logRt_01曲线；
步骤5.利用CNL_01曲线判断出砂体和泥岩；
步骤6.计算CNL_01曲线和AC_01曲线的幅度差DACCNL_01；
步骤7.计算logRt_01曲线和AC_01曲线的幅度差DlogRtAC_01；
步骤8.计算储层有效性指标GAS_FLAG曲线，并利用该曲线对步骤5识别出的砂体进行储层有效性评价；
所述步骤6中的DACCNL_01曲线为步骤1计算得到的AC_01曲线与步骤2计算得到的CNL_01曲线的差值，即DACCNL_01＝AC_01-CNL_01；
所述步骤7中的DlogRtAC_01曲线为步骤3计算得到的logRt_01曲线与步骤1计算得到的AC_01曲线的差值，即DlogRtAC_01＝logRt_01-AC_01；
所述步骤8中的储层有效性指标GAS_FLAG曲线为步骤6计算得到的DACCNL_01曲线与步骤7计算得到的DlogRtAC_01曲线幅度差，即GAS_FLAG＝DACCNL_01-DlogRtAC_01。
2.根据权利要求1所述的一种基于常规测井曲线的致密砂岩储层有效性评价方法，其特征在于，所述步骤1中AC_01曲线的计算公式为：
式中AC_01为中子测井曲线归一化后的曲线，无量纲；AC为实测中子测井曲线，单位为μs/ft；ACmin为研究井段中子曲线最小值，单位为μs/ft；ACmax为研究井段中子曲线最大值，单位为μs/ft。
3.根据权利要求2所述的一种基于常规测井曲线的致密砂岩储层有效性评价方法，其特征在于，所述步骤2中CNL_01曲线的计算公式为：
式中CNL_01为中子测井曲线归一化后的曲线，无量纲；CNL为实测中子测井曲线，单位为％；CNLmin为研究井段中子曲线最小值，单位为％；CNLmax为研究井段中子曲线最大值，单位为％。
4.根据权利要求3所述的一种基于常规测井曲线的致密砂岩储层有效性评价方法，其特征在于，所述步骤3中logRt曲线为实测电阻率曲线的对数值，步骤4中logRt_01曲线的计算公式为：
式中logRt_01为实测电阻率测井曲线取对数后归一化后的曲线，无量纲；logRt为实测电阻率测井曲线的对数，单位为Ω·m；logRtmin为研究井段中子曲线最小值，单位为Ω·m；logRtmax为研究井段中子曲线最大值，单位为Ω·m。
5.根据权利要求4所述的一种基于常规测井曲线的致密砂岩储层有效性评价方法，其特征在于，所述步骤5中，利用步骤2中求取的CNL_01曲线判断砂体和泥岩；其中，CNL_01曲线>0.4时为泥岩，CNL_01曲线≤0.4时为砂体。
6.根据权利要求1所述的一种基于常规测井曲线的致密砂岩储层有效性评价方法，其特征在于，GAS_FLAG曲线小于等于0值的砂体为无效储层，GAS_FLAG大于0值的砂体为有效储层。
7.根据权利要求6所述的一种基于常规测井曲线的致密砂岩储层有效性评价方法，其特征在于，GAS_FLAG值为0.0～0.15时，砂体为含气层；
GAS_FLAG值大于0.15小于等于0.25时，砂体为差气层；
GAS_FLAG值大于0.25时，砂体为气层。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110511151.3A CN113236237B (zh) | 2021-05-11 | 2021-05-11 | 一种基于常规测井曲线的致密砂岩储层有效性评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110511151.3A CN113236237B (zh) | 2021-05-11 | 2021-05-11 | 一种基于常规测井曲线的致密砂岩储层有效性评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113236237A CN113236237A (zh) | 2021-08-10 |
| CN113236237B true CN113236237B (zh) | 2022-05-03 |
Family
ID=77133362
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110511151.3A Active CN113236237B (zh) | 2021-05-11 | 2021-05-11 | 一种基于常规测井曲线的致密砂岩储层有效性评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113236237B (zh) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116299733B (zh) * | 2023-04-06 | 2025-09-16 | 重庆科技学院 | 利用常规测井曲线识别地层中含硫储层的方法 |
| CN119644460A (zh) * | 2023-09-18 | 2025-03-18 | 中国石油天然气股份有限公司 | 一种测井曲线标准化方法、装置、电子设备及存储介质 |
| CN119291775B (zh) * | 2024-12-10 | 2025-02-18 | 东北石油大学三亚海洋油气研究院 | 一种窄河道含气砂体的分布范围预测方法及相关装置 |
Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103233728A (zh) * | 2013-05-15 | 2013-08-07 | 中国石油大学(华东) | 一种LogR-ΔT测井评价的泥页岩油含量预测方法 |
| CA2867583A1 (en) * | 2013-10-15 | 2015-04-15 | Korea Institute Of Geoscience And Mineral Resources (Kigam) | Fracking method for fracking intervals of a horizontal drilling zone in a sweet spot range based on measurements of resistivity and neutron logging data in the horizontal drillingzone |
| CN105370272A (zh) * | 2015-12-07 | 2016-03-02 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 页岩气储层的测井评价方法 |
| CN106021793A (zh) * | 2016-06-01 | 2016-10-12 | 中国地质大学（武汉） | 基于存储系数和渗流系数的低渗储层甜点评价方法 |
| CN107315203A (zh) * | 2017-08-04 | 2017-11-03 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种致密砂岩储层甜点识别及综合评价方法 |
| CN108150161A (zh) * | 2017-12-27 | 2018-06-12 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 页岩含气性评价方法装置 |
| CN108252709A (zh) * | 2016-12-28 | 2018-07-06 | 中国石油化工股份有限公司 | 一种致密砂岩油藏的油水性质识别方法及系统 |
| CN109856674A (zh) * | 2019-03-04 | 2019-06-07 | 西南石油大学 | 工程甜点评测方法及装置 |
| US10371857B1 (en) * | 2013-05-29 | 2019-08-06 | DataInfoCom USA, Inc. | System and method for well log analysis |
| CN110318745A (zh) * | 2019-06-10 | 2019-10-11 | 中国石油大学(华东) | 一种沉积微相约束下粒径岩性测井评价方法 |
| CN111812736A (zh) * | 2020-07-22 | 2020-10-23 | 西南石油大学 | 一种致密砂岩无水气藏含气性评价方法 |
| CN112083515A (zh) * | 2020-09-10 | 2020-12-15 | 西南石油大学 | 致密砂岩低阻储层挖掘效应的定量表征及含气性评价方法 |
-
2021
- 2021-05-11 CN CN202110511151.3A patent/CN113236237B/zh active Active
Patent Citations (12)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103233728A (zh) * | 2013-05-15 | 2013-08-07 | 中国石油大学(华东) | 一种LogR-ΔT测井评价的泥页岩油含量预测方法 |
| US10371857B1 (en) * | 2013-05-29 | 2019-08-06 | DataInfoCom USA, Inc. | System and method for well log analysis |
| CA2867583A1 (en) * | 2013-10-15 | 2015-04-15 | Korea Institute Of Geoscience And Mineral Resources (Kigam) | Fracking method for fracking intervals of a horizontal drilling zone in a sweet spot range based on measurements of resistivity and neutron logging data in the horizontal drillingzone |
| CN105370272A (zh) * | 2015-12-07 | 2016-03-02 | 中国石油集团川庆钻探工程有限公司地球物理勘探公司 | 页岩气储层的测井评价方法 |
| CN106021793A (zh) * | 2016-06-01 | 2016-10-12 | 中国地质大学（武汉） | 基于存储系数和渗流系数的低渗储层甜点评价方法 |
| CN108252709A (zh) * | 2016-12-28 | 2018-07-06 | 中国石油化工股份有限公司 | 一种致密砂岩油藏的油水性质识别方法及系统 |
| CN107315203A (zh) * | 2017-08-04 | 2017-11-03 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 一种致密砂岩储层甜点识别及综合评价方法 |
| CN108150161A (zh) * | 2017-12-27 | 2018-06-12 | 中国石油化工股份有限公司江汉油田分公司勘探开发研究院 | 页岩含气性评价方法装置 |
| CN109856674A (zh) * | 2019-03-04 | 2019-06-07 | 西南石油大学 | 工程甜点评测方法及装置 |
| CN110318745A (zh) * | 2019-06-10 | 2019-10-11 | 中国石油大学(华东) | 一种沉积微相约束下粒径岩性测井评价方法 |
| CN111812736A (zh) * | 2020-07-22 | 2020-10-23 | 西南石油大学 | 一种致密砂岩无水气藏含气性评价方法 |
| CN112083515A (zh) * | 2020-09-10 | 2020-12-15 | 西南石油大学 | 致密砂岩低阻储层挖掘效应的定量表征及含气性评价方法 |
Non-Patent Citations (3)
| Title |
|---|
| 大宁-吉县地区二叠系致密砂岩含气性评价;李小刚等;《云南化工》;20170630;第44卷(第6期);第90-92页 * |
| 常规测井评价复杂砂岩储层孔隙结构的方法与应用———以莺歌海盆地黄流组储层为例;何胜林等;《地球物理学进展》;20171231;第32卷(第6期);第2642-2649页 * |
| 杭锦旗地区二叠系含气性评价;沙文武;《石油天然气学报(江汉石油学院学报)》;20051031;第27卷(第5期);第704-706页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113236237A (zh) | 2021-08-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN113236237B (zh) | 一种基于常规测井曲线的致密砂岩储层有效性评价方法 | |
| CN110847901B (zh) | 一种变矿化度地层水下的致密砂岩储层流体识别方法 | |
| CN106950347B (zh) | 一种评价泥页岩各组分体积的方法 | |
| WO2017024530A1 (zh) | 一种计算烃源岩中有机碳含量的方法 | |
| CN116047602B (zh) | 基于生烃数值模拟的ii型水合物饱和度预测方法 | |
| CN112746835B (zh) | 一种优化的深层页岩气地质甜点测井综合评价方法 | |
| CN103470250B (zh) | 一种测定地层孔隙结构以及流体特性的方法及设备 | |
| CN103775057A (zh) | 一种致密油气藏有效储层的判识方法及装置 | |
| CN105986813A (zh) | 准致密储层快速评价方法及准致密储层多指标评价方法 | |
| CN103291287A (zh) | 一种孔洞型储层有效性级别测定方法 | |
| CN108240952A (zh) | 一种解析法计算页岩含气量的方法 | |
| CN112083515B (zh) | 致密砂岩低阻储层挖掘效应的定量表征及含气性评价方法 | |
| CN104863574A (zh) | 一种适用于致密砂岩储层的流体识别方法 | |
| CN107092036A (zh) | 一种基于储层真电阻率反演的火山岩地层流体识别方法及系统 | |
| CN112943229A (zh) | 一种储气库盖层突破压力连续预测方法 | |
| CN111381292B (zh) | 一种预测砂岩含烃储层的测井解释方法与装置 | |
| CN111812736A (zh) | 一种致密砂岩无水气藏含气性评价方法 | |
| CN112346147A (zh) | 一种基于中子密度孔隙度差的储层评价方法 | |
| CN117569799A (zh) | 利用冲洗带和地层含水孔隙度识别砂砾岩气层与干层方法 | |
| CN117345198A (zh) | 一种陆相页岩油甜点测井分析方法和装置 | |
| Zahaf et al. | Vertical permeability from in situ horizontal measurements in shaly-sand reservoirs | |
| CN113720745A (zh) | 含碳屑碎屑岩储层地球物理测井计算孔隙度的方法 | |
| CN118964792A (zh) | 一种深层致密砂岩气藏有效储层快速划分方法 | |
| CN114295674B (zh) | 一种确定沉积盆地中原位热导率的方法 | |
| CN106907140B (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 |