CN113075748B - 基于成像测井和声波远探测测井资料的裂缝有效性评价方法 - Google Patents
基于成像测井和声波远探测测井资料的裂缝有效性评价方法 Download PDFInfo
- Publication number
- CN113075748B CN113075748B CN202110261157.XA CN202110261157A CN113075748B CN 113075748 B CN113075748 B CN 113075748B CN 202110261157 A CN202110261157 A CN 202110261157A CN 113075748 B CN113075748 B CN 113075748B
- Authority
- CN
- China
- Prior art keywords
- fracture
- logging
- parameters
- effectiveness
- imaging
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 82
- 238000011156 evaluation Methods 0.000 title claims abstract description 53
- 238000001514 detection method Methods 0.000 title claims abstract description 31
- 238000000034 method Methods 0.000 claims abstract description 41
- 238000012512 characterization method Methods 0.000 claims abstract description 34
- 239000011435 rock Substances 0.000 claims description 15
- 238000005070 sampling Methods 0.000 claims description 11
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 10
- 238000004364 calculation method Methods 0.000 claims description 9
- 230000035699 permeability Effects 0.000 claims description 8
- UEUIKXVPXLWUDU-UHFFFAOYSA-N 4-diazoniobenzenesulfonate Chemical compound [O-]S(=O)(=O)C1=CC=C([N+]#N)C=C1 UEUIKXVPXLWUDU-UHFFFAOYSA-N 0.000 claims description 6
- 230000002159 abnormal effect Effects 0.000 claims description 6
- 238000011160 research Methods 0.000 claims description 5
- 238000012937 correction Methods 0.000 claims description 4
- 238000007670 refining Methods 0.000 claims description 3
- 230000015572 biosynthetic process Effects 0.000 claims description 2
- 238000011157 data evaluation Methods 0.000 abstract 1
- 206010017076 Fracture Diseases 0.000 description 219
- 208000010392 Bone Fractures Diseases 0.000 description 211
- 238000011161 development Methods 0.000 description 11
- 239000007789 gas Substances 0.000 description 11
- 238000010586 diagram Methods 0.000 description 9
- 230000000694 effects Effects 0.000 description 6
- 238000012545 processing Methods 0.000 description 6
- 238000012360 testing method Methods 0.000 description 6
- 238000005553 drilling Methods 0.000 description 4
- 208000006670 Multiple fractures Diseases 0.000 description 3
- 238000004519 manufacturing process Methods 0.000 description 3
- 239000011148 porous material Substances 0.000 description 3
- 238000011158 quantitative evaluation Methods 0.000 description 3
- 239000012530 fluid Substances 0.000 description 2
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 150000003839 salts Chemical class 0.000 description 2
- 230000035945 sensitivity Effects 0.000 description 2
- 238000012163 sequencing technique Methods 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 235000019994 cava Nutrition 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000005755 formation reaction Methods 0.000 description 1
- 229930195733 hydrocarbon Natural products 0.000 description 1
- 150000002430 hydrocarbons Chemical class 0.000 description 1
- 230000009545 invasion Effects 0.000 description 1
- 230000033001 locomotion Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000003345 natural gas Substances 0.000 description 1
- 239000003209 petroleum derivative Substances 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 239000003079 shale oil Substances 0.000 description 1
- 239000002002 slurry Substances 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开了一种基于成像测井和声波远探测测井资料的裂缝有效性评价方法，通过利用已知井的成像测井资料和声波远探测测井资料评价，得到多种裂缝属性参数曲线和/或离散样值然后分层段统计特征值，并收集对应各层段的裂缝有效性表征参数；将多种测井裂缝属性参数特征值分别与裂缝有效性表征参数绘制交会图，确定二者之间的定量表征关系及它们之间的相关性，提炼对裂缝有效性敏感的测井裂缝属性参数，并确定有效裂缝和不同等级裂缝对应的测井裂缝属性参数下限值；制定有效裂缝和裂缝等级测井综合评价标准表，用上述综合评价标准表对新井对应的测井裂缝属性参数特征值进行评价，确定裂缝有效性及等级。本发明方法简单、准确性和可靠性好。
Description
技术领域
本发明涉及应用地球物理测井领域，具体的说是一种基于成像测井和声波远探测测井资料的裂缝有效性评价方法。
背景技术
随着石油天然气工业深入发展，大型整装高孔高渗砂岩油气藏越来越难以找到，许多大公司不得不把注意力转移到裂缝－复杂岩性油气藏的勘探和开发上来。致密裂缝型油气藏是21世纪石油增储上产的重要领域之一，在我国，裂缝型低渗透储层和数量比例更为突出，裂缝型低渗储层油气产量占整个石油天然气产量的一半以上，占未来准备投产的石油天然气储量的三分之二以上。
对于低渗致密裂缝性储层，由于其基质的低孔低渗特性，裂缝则作为其主要渗流途径，在孔道间起到连接作用，提高了储层的渗透率，为储层产能的提高提供了基础。故天然裂缝的识别和特征参数及有效性的评价是这类储层评价的一个十分重要的方面。另外，对于页岩气、致密气等非常规储层，往往需要压裂改造，而压裂施工后人造裂缝的发育情况评价同样十分重要。
利用测井资料识别并精细评价裂缝是储层裂缝评价最主要的手段，国内外学者已做了大量相关的研究工作。裂缝测井评价方法主要包括常规测井评价法、成像测井评价法、阵列声波测井评价法和反射声波成像测井评价法。常规测井评价法主要利用声波、密度、中子和深浅电阻率资料识别并评价裂缝的发育情况，侧重于裂缝的定性评价。成像测井评价法主要是利用微电阻率成像和超声成像测井资料对穿过井轴的裂缝进行评价，可以获得裂缝密度(单位井段长度内的裂缝条数)、裂缝宽度(也称裂缝开度，通常指单位长度井段中各种裂缝轨迹的宽度平均值)、裂缝倾角、裂缝长度(通常指每平方米井壁上所有裂缝长度总和)、裂缝面孔率(通常指单位井段上裂缝在井壁上所占面积与成像测井覆盖井壁的面积之比)等裂缝定量参数，被认为是目前可信度最高的一种裂缝测井评价技术。在水基泥浆钻井液条件下，微电阻率成像测井已广泛应用于裂缝的识别和裂缝定量参数的计算，并取得了较好的应用效果。近年来，随着超深致密油气藏、页岩油气藏等非常规油气藏的广泛勘探开发，为了克服井眼垮塌，保护储层等工程问题，提高钻井效率，减少因泥岩膨胀、盐岩蠕变等引起的卡钻事故，大量的井采用了油基钻井液；而油基泥浆的导电性差且侵入特征与水基泥浆不同，导致裂缝处的电阻率与无裂缝处岩石骨架的电阻率差异不大，进而导致基于电学性质的裂缝识别与评价方法应用效果明显变差。超声成像测井不受泥浆电阻率的影响，可以由扫描测量记录到的回波波形处理得到的声阻抗和回波时间参数提供井壁360°方位的成像图，已广泛应用于油基泥浆井裂缝的识别和评价，且一般情况下效果优于微电阻率成像。但微电阻率成像测井和超声波成像测井由于径向探测深度很浅只能反映井壁的情况，无法评价裂缝向井外延伸的情况，对裂缝有效性评价效果有待提高。
声波远探测能够探测井外10米甚至更远的裂缝等声阻抗不连续界面，可以更好地评价裂缝(包括不过井轴的井旁裂缝)向井外延伸情况(裂缝倾角、延伸长度等)，但声波远探测裂缝评价法的分辨率明显低于成像测井裂缝评价法。
为了更好地评价有效裂缝的发育情况，需要一种更准确、客观、有效的对裂缝属性参数和有效性进行综合评价的方法。
发明内容
本发明的目的是为了解决上述技术问题，提供一种方法简单、准确性和可靠性好、可以更好地评价和描述裂缝的有效性、划分有效裂缝等级、指导储层评价的基于成像测井和声波远探测测井资料的裂缝有效性评价方法。
所述评价方法包括以下步骤：
一,利用已知井的成像测井资料定量评价裂缝参数，包括裂缝宽度、密度、长度、倾角、倾向、走向和面孔率，同时利用已知井的声波远探测测井资料评价裂缝向井外延伸长度、倾向、倾角和走向，得到多种测井裂缝属性参数曲线和/或离散样值；
二,对上述多种裂缝属性参数曲线和/或离散样值分层段统计，获得各层段的多种测井裂缝属性参数特征值，并收集该井各对应层段的裂缝有效性表征参数，所述裂缝有效性表征参数为试井渗透率数据或者产能数据；
三,将步骤二中分层统计的多种测井裂缝属性参数特征值分别与裂缝有效性表征参数绘制二者之间的交会图，并通过拟合的方式确定二者之间的定量表征关系及它们之间的相关性，提炼对裂缝有效性敏感的敏感测井裂缝属性参数，并确定有效裂缝和不同等级裂缝对应的敏感测井裂缝属性参数下限值；
四,根据获得的有效裂缝和不同等级裂缝对应的敏感测井裂缝属性参数下限值，制定有效裂缝和裂缝等级测井综合评价标准表；
五,采用步骤一和步骤二的方法获得新井分层段统计的多种测井裂缝属性参数特征值；用步骤四得到的有效裂缝和裂缝等级测井综合评价标准表对上述新井对应的测井裂缝属性参数特征值进行评价，确定裂缝有效性及等级。
所述步骤一中，采用岩心裂缝参数对成像测井资料定量评价的裂缝参数进行刻度校正，具体为：
A1通过岩心观察和描述获取已知井的岩心裂缝参数，并利用岩心地面自然伽马值把岩心深度归位到常规GR曲线统一的深度尺度上；A2将已知井的成像测井资料定量评价裂缝参数的曲线样值也归位到常规GR曲线统一的深度尺度上；
A3对照已知井的岩心照片裂缝参数和成像测井资料定量评价裂缝参数，分析两者存在的关系，得到两者之间的刻度系数，达到岩心照片裂缝参数刻度成像测井资料定量评价裂缝参数的目的。
所述步骤二中，所述测井裂缝属性参数特征值为对应各层段测井裂缝属性参数的最大值、最小值、中值、均方根值、加权平均值或算术平均值等。
所述步骤二中，当多种测井裂缝属性参数特征值为对应各层段的算术平均值时，采用下述方法来统计：
B1统计任一种测井裂缝属性参数曲线上对应层段内的全部样值；
B2对上述全部样值进行排序处理，找出样点中值fpM(排序后序号位于中间的样值，若样点总数为偶数，取排序后序号位于中间两个样值的算术平均值)；
B3计算上述全部样值与样点中值之差的绝对值fDABS，计算公式如式(4)所示；
B4从上述全部样值中剔除fDABS最大的10％对应的异常样点值，并按式(5)计算剩余样值的算术平均值作为该段测井裂缝属性参数的特征值；
B5重复B1-B4，统计其它的测井测井裂缝属性参数特征值；
fDABSi＝|fpi-fpM| (4)
其中，fDABSi为层段内测井裂缝属性参数曲线第i个样值fpi与上述样点中值fpM之差的绝对值，N为层段内测井裂缝属性参数曲线样点剔除异常样点值后剩余的样点总数，fpj为剩余样点中第j个样点值，为层段内测井裂缝属性参数特征值。
所述步骤二中，所述产能数据为无阻流量或米产能指数。
所述步骤三中，当二者的相关性系数R2≥0.4时，说明二者具有较好的相关性，则提炼该测井裂缝属性参数特征值为对裂缝有效性敏感的敏感测井裂缝属性参数。
所述步骤三中，确定有效裂缝和不同等级裂缝对应的敏感测井裂缝属性参数下限值，具体方法如下：
具体方法如下：
C1统计研究区不同等级储层对应的裂缝有效性表征参数下限值；
C 2由上述确定的测井裂缝属性参数特征值与裂缝有效性表征参数之间的定量关系确定上述裂缝有效性表征参数下限值对应的测井裂缝属性参数特征值，并将它们作为有效裂缝和不同等级裂缝对应的敏感测井裂缝属性参数下限值。
所述步骤一中，当评价对象为水基泥浆井，则采用微电阻率成像测井资料定量评价裂缝参数；当评价对象为油基泥浆井，则采用超声成像测井数据资料定量评价裂缝参数。
所述步骤一中，利用已知井的声波远探测测井资料评价裂缝向井外延伸长度、倾向、倾角和走向，其中，在直接获得裂缝走向和倾向参数曲线的基础上形成井旁构造成像图，然后在该图中拾取裂缝向井外延伸长度和倾角的参数离散样值。
有益效果：
1)利用成像测井和声波远探测测井资料同步评价裂缝属性参数，不同方法评价结果可以相互验证，提高裂缝属性参数评价的准确性和可靠性；
2)本发明方法充分利用两种不同尺度测井方法在裂缝评价中的优势，二者相互补充，由此获得的有效裂缝和裂缝等级测井综合评价标准表可以更为准确和可靠的地评价和描述裂缝的有效性，并可进一步划分有效裂缝等级；
3)本发明方法可以用于水基泥浆井和油基泥浆井的裂缝定量参数和有效性评价，并用于指导储层评价，提高复杂非常规储层测井评价的准确性和可靠性，进而为制定合理高效的开发方案提供有力依据。
附图说明
图1为本发明方法流程图；
图2为本实施例中微电阻率成像测井资料计算裂缝参数成果图；
图3为利用岩心裂缝参数刻度成像裂缝参数流程图；
图4为X地区水基泥浆井岩心裂缝宽度与电成像裂缝宽度交会图；
图5为X地区水基泥浆井岩心裂缝面孔率与电成像面孔率交会图；
图6为X井声波远探测裂缝成像图；
图7为远探测裂缝倾角和延伸长度计算方法示意图；
图8为无阻流量与裂缝宽度交会图；
图9为不同等级裂缝对应的裂缝面孔率下限值确定方法示意图。
具体实施方式
参见图1，以对X地区水基泥浆井的评价为例对本发明方法作进一步解释说明：
一，多种测井裂缝属性参数特征值曲线的获得：
1)利用已知井的成像测井资料定量评价裂缝参数，包括裂缝宽度、密度、长度、倾角、倾向、走向和面孔率：
成像测井裂缝参数计算(即定量评价)即主要通过测井资料处理软件实现的，包括但不限于裂缝密度、裂缝宽度、裂缝长度、裂缝面孔率、裂缝倾角、裂缝倾向、裂缝走向等参数。需要注意是，对于水基泥浆井，优先采用微电阻率成像测井数据计算上述裂缝属性参数；对于油基泥浆井，优先采用超声成像测井数据计算上述裂缝属性参数，以提高裂缝参数计算结果的准确性和可靠性。图2为微电阻率成像测井资料计算裂缝参数成果图。所述已知井可以为一口或多口井。
采用岩心裂缝参数对成像测井资料定量评价的裂缝参数进行刻度校正，岩心观察与描述能提供关于裂缝参数和发育情况、力学性质、充填特征和含油气性等的第一手资料，是最直接、最有效和最可靠的裂缝评价方式，故评价的裂缝参数精度较高。为了进一步提高成像裂缝参数评价的精度，利用岩心裂缝参数对成像裂缝参数进行刻度校正，成像裂缝参数的刻度实质上是通过对比分析成像资料处理得到的裂缝参数和岩心观察与描述得到的裂缝参数，建立二者之间的刻度系数，以用于校正成像资料处理得到的裂缝参数，最终达到精细评价储层裂缝发育情况的目的，参见图3，具体为：
A1通过岩心观察和描述获取已知井的岩心裂缝参数，并利用岩心地面自然伽马值把岩心深度归位到常规GR曲线统一的深度尺度上；A2将已知井的成像测井资料定量评价得到的裂缝参数的曲线样值也归位到常规GR曲线统一的深度尺度上；
A3对照已知井的岩心照片裂缝参数和成像测井资料定量评价裂缝参数，分析两者存在的关系，得到两者之间的刻度系数，达到岩心照片裂缝参数刻度成像测井资料定量评价裂缝参数的目的。
采用上述方法统计X地区7口水基泥浆井的岩心裂缝参数和微电阻率成像测井资料(简称电成像)处理得到的裂缝参数，并分别绘制岩心裂缝参数和成像裂缝参数的交会图(如图4和图5所示)，图中直线为线性拟合趋势线。由图4可以看出，整体上电成像裂缝宽度约为岩心裂缝宽度的10.926倍；由图5可以看出，整体上电成像裂缝面孔率约为岩心裂缝面孔率的6.46倍；这两个系数就是根据岩心裂缝参数和电成像裂缝参数之间的关系得到的电成像裂缝参数刻度系数。实际电成像资料处理时利用这个系数对处理得到的裂缝参数进行刻度校正后可以更好地反映裂缝的真实状态和发育情况，具体电成像裂缝参数刻度公式为：
其中，FVAH为刻度校正后的裂缝宽度，FVAHFMI为电成像处理得到的裂缝宽度；FVPA为刻度校正后的裂缝面孔率，FVPAFM为电成像处理得到的裂缝面孔率。
2)利用已知井的声波远探测测井资料评价裂缝向井外延伸长度、倾向、倾角和走向：
声波远探测测井(也称为反射声波成像测井)以辐射到井外地层中的声场能量作为入射波，探测从井旁地层界面、溶洞、裂缝或小构造等引起的声阻抗不连续界面反射回来的声场；通过分析接收到的全波列阵列波形中的反射波信号，可以对井旁存在的小地质构造进行成像。该方法是一种前沿的先进的声波测井方法，可对井周几米到十几米甚至几十米范围内的地层界面、断层、裂缝、尖灭、溶洞或盐丘等地质构造或地质体进行探测，且其分辨率和径向探测深度恰好介于地震勘探和常规声波测井之间，从而填补了油气勘探领域的一项空白。
声波远探测测井主要用于评价井旁裂缝(包括过井轴的和未过井轴的裂缝)、裂缝向井外延伸长度、裂缝倾向、倾角和走向等参数，一般在成熟的测井资料处理分析平台上即可对声波远探测测井资料进行处理，得到裂缝走向和倾向的曲线，如图6所示为X井声波远探测裂缝成像图，图中椭圆虚线框标示的为裂缝带成像；根据图6的成像结果可进一步通过人工拾取裂缝并计算裂缝倾角和延伸长度参数的离散样值，如图7所示，根据裂缝在成像图中的响应轨迹可以通过人工的方式拾取出如图中倾斜虚线所示的裂缝轨迹，并由裂缝轨迹末端作竖直虚线，由该虚线与横轴的交点可以得到裂缝末端离井轴的距离L1＝13.6m，由裂缝轨迹末端作水平虚线与深度轴相交，该交点到裂缝与井轴的交点距离为L2＝13.0m；根据参数L1和L2可进一步计算出裂缝向井外延伸长度L和裂缝与井轴的夹角α1，具体计算公式如下：，
α1＝arctan(L1/L2) (3)
如果井孔为竖直的，则裂缝与井轴的夹角α1即为裂缝的倾角α；如果井孔为倾斜的，则可通过井孔倾斜方位、倾斜角及裂缝倾向、裂缝与井轴的夹角计算裂缝倾角。例如，设井孔倾斜角为α2，若井孔倾斜方位与裂缝倾向相同，则裂缝倾角α＝α1+α2，若井孔倾斜方位与裂缝倾向相反，则裂缝倾角α＝α1-α2。
所述X井为竖直井，则按照上述方法可以计算出该裂缝带向井外延伸长度约为18.81m，裂缝倾角为46.3°，说明该裂缝带发育的裂缝有效性好。
上述拾取和计算过程可以采用人工或通过分析软件自动完成。
二,对上述多种裂缝属性参数曲线和/或离散样值分层段统计，获得各层段的多种测井裂缝属性参数特征值，收集已知井对应各层段的裂缝有效性表征参数；
1)上述多种裂缝属性参数曲线和/或离散样值分层段统计，获得各层段的多种测井裂缝属性参数特征值。
测井资料计算的裂缝属性参数曲线采样间隔一般为测井时仪器深度移动间隔，如声波远探测测井一般深度采样间隔为0.125m；为了便于分析，往往需要分层段统计裂缝属性参数的特征值，如最大值、最小值、平均值或技术人员认为其它的合适取值等。分层段的方式可以按固定的深度段分层，如每2m划分为一层；也可以按实际试井渗透率测试或油气开采测试层段来分层。
如以算术平均值作为层段裂缝属性参数的特征值，则可按如下方式统计：
B1统计任一种裂缝属性参数曲线上对应层段内的全部样值；
B2对上述全部样值进行排序处理，找出样点中值fpM；
B3计算上述全部样值与样点中值之差的绝对值fDABS，计算公式如式(4)所示；
B4从上述全部样值中剔除fDABS最大的10％对应的异常样点值，并按(5)式计算剩余样值的算术平均值作为该段测井裂缝属性参数的特征值；
B5重复B1-B4，统计其它的测井裂缝属性参数特征值；
fDABSi＝|fpi-fpM| (4)
其中，fDABSi为层段内测井裂缝属性参数曲线第i个样值fpi与上述样点中值fpM之差的绝对值，N为层段内测井裂缝属性参数曲线样点剔除异常样点值后剩余的样点总数，fpj为剩余样点中第j个样点值，为层段内测井裂缝属性参数特征值。
对于离散样值参数，则可直接按照基本的算术平均值算法计算层段裂缝属性参数特征值。
2)分层段统计裂缝有效性表征参数：
一般情况下，有效裂缝的发育能够大大改善储层的渗透性，提高产能，裂缝有效性等级越高，这种效应越明显。故选择试井渗透率数据或者产能数据(所述产能数据可以为无阻流量或米产能指数等)作为裂缝有效性和等级的表征参数。分层段统计裂缝有效性和等级表征参数，以便后续研究分析其与裂缝属性参数之间的关系，进而提炼对裂缝有效性敏感的测井裂缝属性参数。需要说明的是裂缝有效性表征参数统计层段划分情况应与测井裂缝属性参数特征值统计层段划分情况一致。
三，将步骤二中分层统计的多种测井裂缝属性参数特征值分别与裂缝有效性表征参数绘制二者之间的交会图，并通过拟合的方式确定二者之间的定量表征关系及它们之间的相关性，提炼对裂缝有效性敏感的测井裂缝属性参数。
图8所示为统计分析得到的无阻流量与裂缝面孔率交会图，可以看出，无阻流量随着裂缝面孔率的增大而增大，通过拟合的方式得到无阻流量与裂缝面孔率之间的定量关系如式(6)所示，二者之间的相关系数R2＝0.7616，说明二者相关性很好，故裂缝面孔率为裂缝有效性敏感测井参数。
依此类推，可以分析其它测井裂缝属性参数特征值与裂缝有效性表征参数之间的相关性，进而提炼裂缝有效性敏感测井参数。本发明中，测井裂缝参数与裂缝有效性表征参数之间的相关性系数R2≥0.4时，说明二者具有较好的相关性，则该测井裂缝属性参数特征值为裂缝有效性敏感测井参数。
确定有效裂缝和不同等级裂缝对应的测井裂缝属性参数下限值，具体为具体方法如下：
C1统计研究区不同等级储层对应的裂缝有效性表征参数下限值；
C2由上述确定的测井裂缝属性参数特征值与裂缝有效性表征参数之间的定量关系，确定所述裂缝有效性表征参数下限值对应的测井裂缝属性参数特征值，并将它们作为有效裂缝和不同等级裂缝对应的敏感测井裂缝属性参数下限值。
如某油田将无阻流量小于10万方/天的储层划分为低产储层，将无阻流量介于(10,50)万方/天的储层划分为中产储层，将无阻流量大于50万方/天的储层划分为高产储层，则说明10万方/天为中产储层无阻流量下限值，50万方/天为高产储层无阻流量下限值。
四,根据获得的有效裂缝和不同等级裂缝对应的测井裂缝属性参数下限值，制定有效裂缝和裂缝等级测井综合评价标准表。
根据上述无阻流量与裂缝面孔率交会图确定的二者之间定量表征关系曲线可由上述两个无阻流量下限值确定对应的裂缝面孔率下限值，如图9所示确定的两个裂缝宽度下限值分别为90μm和150μm，这两个值可分别定为有效裂缝和I类有效裂缝对应的裂缝宽度下限值，由此可以将裂缝宽度小于90μm的裂缝划分为无效裂缝，将裂缝宽度介于(90,150)μm之间的裂缝划分为II类有效裂缝，将裂缝宽度大于150μm的裂缝划分为I类有效裂缝。
依此类推可以确定不同等级裂缝对应的其它测井裂缝参数下限值，进而建立如表1所示的裂缝有效性及等级综合评价标准。
表1裂缝有效性及等级综合评价标准表
五,采用步骤一和步骤二的方法获得新井分层段统计的多种测井裂缝属性参数特征值；用步骤四得到的有效裂缝和裂缝等级测井综合评价标准表对上述新井对应的测井裂缝属性参数特征值进行评价，确定裂缝有效性及等级，本实施例综合评价标准表中对应评价的裂缝有效性敏感测井参数分别为裂缝面孔率、裂缝宽度、最大主应力方向与裂缝走向夹角、裂缝向井外延伸长度和裂缝倾角这五个参数。
需要说明的是，若某层段计算的裂缝属性参数没有全部位于同一类裂缝对应的裂缝属性参数标准范围内，则可按“少数服从多数”的原则或加权平均的方法确定最终的裂缝有效性和等级。本实施例中，新井的裂缝有效性敏感参数及裂缝有效性和等级评价结果如下表2所示：表2新井的裂缝有效性敏感参数及裂缝有效性和等级评价结果
上述方法在实际储层裂缝评价中的应用结果表明，该方法评价的裂缝有效性和等级与测试得到的裂缝有效性和等级表征参数指示的裂缝等级结果一致。
Claims (8)
1.一种基于成像测井和声波远探测测井资料的裂缝有效性评价方法，其特征在于，包括以下步骤：
一,利用已知井的成像测井资料定量评价裂缝参数，包括裂缝宽度、密度、长度、倾角、倾向、走向和面孔率，同时利用已知井的声波远探测测井资料评价裂缝向井外延伸长度、倾向、倾角和走向，得到多种测井裂缝属性参数曲线和/或离散样值；
其中，采用岩心裂缝参数对成像测井资料定量评价的裂缝参数进行刻度校正，具体为：
A1通过岩心观察和描述获取已知井的岩心裂缝参数，并利用岩心地面自然伽马值把岩心深度归位到常规GR曲线统一的深度尺度上；
A2将已知井的成像测井资料定量评价裂缝参数的曲线样值也归位到常规GR曲线统一的深度尺度上；
A3对照已知井的岩心照片裂缝参数和成像测井资料定量评价裂缝参数，分析两者存在的关系，得到两者之间的刻度系数，达到岩心照片裂缝参数刻度成像测井资料定量评价裂缝参数的目的；
二,对上述多种测井裂缝属性参数曲线和/或离散样值分层段统计，获得各层段的多种测井裂缝属性参数特征值，并收集该井各对应层段的裂缝有效性表征参数，所述裂缝有效性表征参数为试井渗透率数据或者产能数据；
三,将步骤二中分层统计的多种测井裂缝属性参数特征值分别与裂缝有效性表征参数绘制二者之间的交会图，并通过拟合的方式确定二者之间的定量表征关系及它们之间的相关性，提炼对裂缝有效性敏感的敏感测井裂缝属性参数，并确定有效裂缝和不同等级裂缝对应的敏感测井裂缝属性参数下限值；
四,根据获得的有效裂缝和不同等级裂缝对应的敏感测井裂缝属性参数下限值，制定有效裂缝和裂缝等级测井综合评价标准表；
五,采用步骤一和步骤二的方法获得新井分层段统计的多种测井裂缝属性参数特征值；用步骤四得到的有效裂缝和裂缝等级测井综合评价标准表对上述新井分层段统计的多种测井裂缝属性参数特征值进行评价，确定裂缝有效性及等级。
2.如权利要求1所述的基于成像测井和声波远探测测井资料的裂缝有效性评价方法，其特征在于，
所述步骤二中，所述测井裂缝属性参数特征值为对应各层段测井裂缝属性参数的最大值、最小值、中值、均方根值、加权平均值或算术平均值。
3.如权利要求2所述的基于成像测井和声波远探测测井资料的裂缝有效性评价方法，其特征在于，
所述步骤二中，当多种测井裂缝属性参数特征值为对应各层段的算术平均值时，采用下述方法来统计：
B1统计任一种测井裂缝属性参数曲线上对应层段内的全部样值；B2对上述全部样值进行排序处理，找出样点中值fpM；
B3计算上述全部样值与样点中值之差的绝对值fDABS，计算公式如式(4)所示；
B4从上述全部样值中剔除fDABS最大的10％对应的异常样点值，并按(5)式计算剩余样值的算术平均值作为该段测井裂缝属性参数的特征值；
B5重复B1-B4，统计其它的测井裂缝属性参数特征值；
fDABSi＝|fpi-fpM| (4)
其中，fDABSi为层段内测井裂缝属性参数曲线第i个样值fpi与上述样点中值fpM之差的绝对值，N为层段内测井裂缝属性参数曲线样点剔除异常样点值后剩余的样点总数，fpj为剩余样点中第j个样点值，为层段内测井裂缝属性参数特征值。
4.如权利要求1或2或3所述的基于成像测井和声波远探测测井资料的裂缝有效性评价方法，其特征在于，所述步骤二中，所述产能数据为无阻流量或米产能指数。
5.如权利要求1所述的基于成像测井和声波远探测测井资料的裂缝有效性评价方法，其特征在于，
所述步骤三中，当二者的相关性系数R2≥0.4时，说明二者具有较好的相关性，则提炼该测井裂缝属性参数特征值为对裂缝有效性敏感的敏感测井裂缝属性参数。
6.如权利要求1或5所述的基于成像测井和声波远探测测井资料的裂缝有效性评价方法，其特征在于，
所述步骤三中，确定有效裂缝和不同等级裂缝对应的敏感测井裂缝属性参数下限值，具体方法如下：
C1统计研究区不同等级储层对应的裂缝有效性表征参数下限值；
C2由上述确定的测井裂缝属性参数特征值与裂缝有效性表征参数之间的定量关系确定上述裂缝有效性表征参数下限值对应的测井裂缝属性参数特征值，并将它们作为有效裂缝和不同等级裂缝对应的敏感测井裂缝属性参数下限值。
7.如权利要求1所述的基于成像测井和声波远探测测井资料的裂缝有效性评价方法，其特征在于，所述步骤一中，当评价对象为水基泥浆井，则采用微电阻率成像测井资料定量评价裂缝参数；当评价对象为油基泥浆井，则采用超声成像测井数据资料定量评价裂缝参数。
8.如权利要求1所述的基于成像测井和声波远探测测井资料的裂缝有效性评价方法，其特征在于，所述步骤一中，利用已知井的声波远探测测井资料评价裂缝向井外延伸长度、倾向、倾角和走向，其中，在直接获得裂缝走向和倾向参数曲线的基础上形成井旁构造成像图，然后在该图中拾取裂缝向井外延伸长度和倾角参数的离散样值。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110261157.XA CN113075748B (zh) | 2021-03-10 | 2021-03-10 | 基于成像测井和声波远探测测井资料的裂缝有效性评价方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110261157.XA CN113075748B (zh) | 2021-03-10 | 2021-03-10 | 基于成像测井和声波远探测测井资料的裂缝有效性评价方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113075748A CN113075748A (zh) | 2021-07-06 |
| CN113075748B true CN113075748B (zh) | 2024-04-16 |
Family
ID=76612330
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110261157.XA Active CN113075748B (zh) | 2021-03-10 | 2021-03-10 | 基于成像测井和声波远探测测井资料的裂缝有效性评价方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113075748B (zh) |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118011509A (zh) * | 2024-04-09 | 2024-05-10 | 中国石油大学(华东) | 地层裂缝有效性评价方法及系统 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103376469A (zh) * | 2012-04-26 | 2013-10-30 | 中国石油集团长城钻探工程有限公司 | 一种基于超声成像测井的裂缝定量评价方法 |
| CN109425912A (zh) * | 2017-09-05 | 2019-03-05 | 中国石油化工股份有限公司 | 一种碳酸盐岩油藏储层有效裂缝的识别方法 |
| CN110208859A (zh) * | 2019-05-07 | 2019-09-06 | 长江大学 | 基于超声成像的油基泥浆井裂缝定量参数智能计算方法 |
| CN111781638A (zh) * | 2020-06-02 | 2020-10-16 | 长江大学 | 一种有效裂缝带预测方法及装置 |
| CN111880234A (zh) * | 2020-06-24 | 2020-11-03 | 西南石油大学 | 基于常规测井对致密储层裂缝的精细识别方法 |
Family Cites Families (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| WO2017127058A1 (en) * | 2016-01-20 | 2017-07-27 | Halliburton Energy Services, Inc. | Fracture interpretation with resistivity and sonic logs in biaxial anisotropic formations |
-
2021
- 2021-03-10 CN CN202110261157.XA patent/CN113075748B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN103376469A (zh) * | 2012-04-26 | 2013-10-30 | 中国石油集团长城钻探工程有限公司 | 一种基于超声成像测井的裂缝定量评价方法 |
| CN109425912A (zh) * | 2017-09-05 | 2019-03-05 | 中国石油化工股份有限公司 | 一种碳酸盐岩油藏储层有效裂缝的识别方法 |
| CN110208859A (zh) * | 2019-05-07 | 2019-09-06 | 长江大学 | 基于超声成像的油基泥浆井裂缝定量参数智能计算方法 |
| CN111781638A (zh) * | 2020-06-02 | 2020-10-16 | 长江大学 | 一种有效裂缝带预测方法及装置 |
| CN111880234A (zh) * | 2020-06-24 | 2020-11-03 | 西南石油大学 | 基于常规测井对致密储层裂缝的精细识别方法 |
Non-Patent Citations (1)
| Title |
|---|
| 油基泥浆下超声成像测井仪裂缝识别效果探析;雷均;《中国石油和化工标准与质量》;第103、104页 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113075748A (zh) | 2021-07-06 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN103135131B (zh) | 一种针对裂缝性储层预测的解释装置 | |
| CN103821495B (zh) | 测井方法 | |
| CN104514552A (zh) | 一种煤层气储层识别与丰度预测的方法 | |
| CN111596364A (zh) | 基于高精度层序地层格架的地震沉积微相组合分析方法 | |
| KR101853334B1 (ko) | 순차층서 분석을 통한 비전통자원 개발 구간 예측 방법 | |
| CN112835124B (zh) | 基于成像测井和阵列声波测井资料的裂缝有效性评价方法 | |
| CN103678778B (zh) | 一种放射性物化探信息综合方法 | |
| CN111399056A (zh) | 一种基于划分方位角滤波预测裂缝强度的方法 | |
| CN107304674B (zh) | 一种利用钻井液核磁共振评价储层含油性的录井方法 | |
| CN107218033A (zh) | 一种识别地层产状的方法和地层参数的反演方法 | |
| US20130292111A1 (en) | Method of constructing a well log of a quantitative property from sample measurements and log data | |
| CN113075748B (zh) | 基于成像测井和声波远探测测井资料的裂缝有效性评价方法 | |
| CN113216945B (zh) | 一种致密砂岩储层渗透性定量评价方法 | |
| US20220120933A1 (en) | Method of detection of hydrocarbon horizontal slippage passages | |
| CN113050168B (zh) | 基于阵列声波测井和声波远探测测井资料的裂缝有效性评价方法 | |
| CN110847887B (zh) | 一种细粒沉积陆相页岩裂缝识别评价方法 | |
| CN115857047B (zh) | 一种地震储层综合预测方法 | |
| CN115965273A (zh) | 一种页岩油水平井钻探过程中甜点评价方法 | |
| CN114167515B (zh) | 一种岩性圈闭有效性识别的方法 | |
| CN110454155B (zh) | 一种应用磁化率测井方法确定第四纪地层年代的方法 | |
| CN1538166A (zh) | 固井质量综合检测方法 | |
| CN109598049B (zh) | 钻孔岩石裂隙发育程度及区域岩石裂隙发育规律的方法 | |
| CN111963162A (zh) | 基于岩性与水性的双参数流体性质识别方法 | |
| CN115469366B (zh) | 一种油基泥浆微电阻率扫描成像测井参数确定方法 | |
| Gomez et al. | Quantifying fracture intensity: An example from the Piceance Basin |
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 | ||
| TA01 | Transfer of patent application right | ||
| TA01 | Transfer of patent application right |
Effective date of registration: 20211130 Address after: 100028 Chaoyang District, Beijing Hui Xin Street 6, Twelfth level. Applicant after: SINOPEC OILFIELD SERVICE Corp. Applicant after: Yangtze University Applicant after: SINOPEC OILFIELD SERVICE JIANGHAN Corp. Applicant after: Sinopec Jingwei Co.,Ltd. Applicant after: Jianghan logging branch of Sinopec Jingwei Co.,Ltd. Address before: 100028 Chaoyang District, Beijing Hui Xin Street 6, Twelfth level. Applicant before: SINOPEC OILFIELD SERVICE Corp. Applicant before: Yangtze University Applicant before: SINOPEC OILFIELD SERVICE JIANGHAN Corp. |
|
| GR01 | Patent grant | ||
| GR01 | Patent grant |