CN109254330B - 火成岩储层的裂缝段的识别方法及装置 - Google Patents

火成岩储层的裂缝段的识别方法及装置 Download PDF

Info

Publication number
CN109254330B
CN109254330B CN201710566959.5A CN201710566959A CN109254330B CN 109254330 B CN109254330 B CN 109254330B CN 201710566959 A CN201710566959 A CN 201710566959A CN 109254330 B CN109254330 B CN 109254330B
Authority
CN
China
Prior art keywords
data
reservoir
sampling point
group
normalized data
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
CN201710566959.5A
Other languages
English (en)
Other versions
CN109254330A (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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN201710566959.5A priority Critical patent/CN109254330B/zh
Publication of CN109254330A publication Critical patent/CN109254330A/zh
Application granted granted Critical
Publication of CN109254330B publication Critical patent/CN109254330B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01VGEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00Prospecting 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)
  • Geophysics And Detection Of Objects (AREA)
  • Analysing Materials By The Use Of Radiation (AREA)

Abstract

本发明公开了一种火成岩储层的裂缝段的识别方法及装置,属于油气藏开发技术领域。该方法包括:获取油气井的火成岩储层中的多个采样点的测井数据,该测井数据包括密度、补偿中子和深侧向电阻率;基于多个采样点的测井数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子;对于该多个采样点中的每个采样点,如果该采样点所在的储层段的裂缝发育因子不小于第一预设值,则确定该采样点所在的储层段为火成岩储层的裂缝段。由于该多个采样点的测井数据能够简单方便地获取得到,且本发明实施例中仅基于该多个采样点的测井数据,就能完成对火成岩储层的裂缝段的识别,因此,将大大节省火成岩储层的裂缝段的识别成本,且提高识别效率。

Description

火成岩储层的裂缝段的识别方法及装置
技术领域
本发明涉及油气藏开发技术领域,特别涉及一种火成岩储层的裂缝段的识别方法及装置。
背景技术
随着油气藏开发技术的不断发展,火成岩储层以其产层厚、产率高、储量大等特点,已经逐渐成为重要的油气藏开发目标。由于火成岩储层的裂缝段的渗透能力往往比较高,所以如果在火成岩储层的裂缝段进行油气开发,则可以大大提高油气产量。为此,在对火成岩储层进行油气开发之前,需要识别火成岩储层的裂缝段。
目前,识别火成岩储层的裂缝段时,会先对在该火成岩储层中钻开的油气井进行取芯,得到多个岩芯,再对该多个岩芯进行分析,以识别该火成岩储层的裂缝段。然而,油气井的取芯过程通常比较复杂,且后续的岩芯分析过程也比较繁琐,从而导致火成岩储层的裂缝段的识别效率较低,识别成本较高。
发明内容
为了解决相关技术中火成岩储层的裂缝段的识别效率较低且识别成本较高的问题,本发明实施例提供了一种火成岩储层的裂缝段的识别方法及装置。所述技术方案如下:
一方面,提供了一种火成岩储层的裂缝段的识别方法,所述方法包括:
获取油气井的火成岩储层中的多个采样点的测井数据,所述测井数据包括密度、补偿中子和深侧向电阻率;
基于所述多个采样点的测井数据,确定所述多个采样点中每个采样点所在的储层段的裂缝发育因子;
对于所述多个采样点中的每个采样点,如果所述采样点所在的储层段的裂缝发育因子不小于第一预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段。
可选地,所述基于所述多个采样点的测井数据,确定所述多个采样点中每个采样点所在的储层段的裂缝发育因子,包括:
将所述多个采样点的测井数据划分为第一组测井数据、第二组测井数据和第三组测井数据,所述第一组测井数据为所述多个采样点的密度,所述第二组测井数据为所述多个采样点的补偿中子,所述第三组测井数据为所述多个采样点的深侧向电阻率;
对所述第一组测井数据进行归一化处理,得到第一组归一化数据,对所述第二组测井数据进行归一化处理,得到第二组归一化数据,对所述第三组测井数据进行归一化处理,得到第三组归一化数据;
基于所述第一组归一化数据、所述第二组归一化数据和所述第三组归一化数据,确定所述多个采样点中每个采样点所在的储层段的裂缝发育因子。
可选地,所述基于所述第一组归一化数据、所述第二组归一化数据和所述第三组归一化数据,确定所述多个采样点中每个采样点所在的储层段的裂缝发育因子,包括:
分别确定所述第一组归一化数据的平均值、所述第二组归一化数据的平均值和所述第三组归一化数据的平均值;
对于所述多个采样点中的每个采样点,分别从所述第一组归一化数据、所述第二组归一化数据和所述第三组归一化数据中获取所述采样点的测井数据的归一化数据;
基于所述第一组归一化数据的平均值、所述第二组归一化数据的平均值、所述第三组归一化数据的平均值和所述采样点的测井数据的归一化数据,确定所述采样点所在的储层段的裂缝发育因子。
可选地,所述基于所述第一组归一化数据的平均值、所述第二组归一化数据的平均值、所述第三组归一化数据的平均值和所述采样点的测井数据的归一化数据,确定所述采样点所在的储层段的裂缝发育因子,包括:
基于所述第一组归一化数据的平均值、所述第二组归一化数据的平均值、所述第三组归一化数据的平均值和所述采样点的测井数据的归一化数据,通过如下第一指定公式确定所述采样点所在的储层段的裂缝发育因子;
第一指定公式:
Figure BDA0001348634010000021
其中,K为所述采样点所在的储层段的裂缝发育因子,DEN为所述采样点的密度的归一化数据,CNL是为所述采样点的补偿中子的归一化数据,
Figure BDA0001348634010000031
为所述第一组归一化数据的平均值,
Figure BDA0001348634010000032
为所述第二组归一化数据的平均值,
Figure BDA0001348634010000033
为所述第三组归一化数据的平均值。
可选地,所述如果所述采样点所在的储层段的裂缝发育因子不小于第一预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段,包括:
如果所述采样点所在的储层段的裂缝发育因子不小于所述第一预设值且不大于第二预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段中的次裂缝发育区;
如果所述采样点所在的储层段的裂缝发育因子大于所述第二预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段中的裂缝发育区。
另一方面,提供了一种火成岩储层的裂缝段的识别装置,所述装置包括:
获取模块,用于获取油气井的火成岩储层中的多个采样点的测井数据,所述测井数据包括密度、补偿中子和深侧向电阻率;
第一确定模块,用于基于所述多个采样点的测井数据,确定所述多个采样点中每个采样点所在的储层段的裂缝发育因子;
第二确定模块,用于对于所述多个采样点中的每个采样点,如果所述采样点所在的储层段的裂缝发育因子不小于第一预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段。
可选地,所述第一确定模块包括:
划分子模块,用于将所述多个采样点的测井数据划分为第一组测井数据、第二组测井数据和第三组测井数据,所述第一组测井数据为所述多个采样点的密度,所述第二组测井数据为所述多个采样点的补偿中子,所述第三组测井数据为所述多个采样点的深侧向电阻率;
处理子模块,用于对所述第一组测井数据进行归一化处理,得到第一组归一化数据,对所述第二组测井数据进行归一化处理,得到第二组归一化数据,对所述第三组测井数据进行归一化处理,得到第三组归一化数据;
第一确定子模块,用于基于所述第一组归一化数据、所述第二组归一化数据和所述第三组归一化数据,确定所述多个采样点中每个采样点所在的储层段的裂缝发育因子。
可选地,所述第一确定子模块包括:
第一确定单元,用于分别确定所述第一组归一化数据的平均值、所述第二组归一化数据的平均值和所述第三组归一化数据的平均值;
获取单元,用于对于所述多个采样点中的每个采样点,分别从所述第一组归一化数据、所述第二组归一化数据和所述第三组归一化数据中获取所述采样点的测井数据的归一化数据;
第二确定单元,用于基于所述第一组归一化数据的平均值、所述第二组归一化数据的平均值、所述第三组归一化数据的平均值和所述采样点的测井数据的归一化数据,确定所述采样点所在的储层段的裂缝发育因子。
可选地,所述第二确定单元主要用于:
基于所述第一组归一化数据的平均值、所述第二组归一化数据的平均值、所述第三组归一化数据的平均值和所述采样点的测井数据的归一化数据,通过如下第一指定公式确定所述采样点所在的储层段的裂缝发育因子;
第一指定公式:
Figure BDA0001348634010000041
其中,K为所述采样点所在的储层段的裂缝发育因子,DEN为所述采样点的密度的归一化数据,CNL是为所述采样点的补偿中子的归一化数据,
Figure BDA0001348634010000042
为所述第一组归一化数据的平均值,
Figure BDA0001348634010000043
为所述第二组归一化数据的平均值,
Figure BDA0001348634010000044
为所述第三组归一化数据的平均值。
可选地,所述第二确定模块包括:
第二确定子模块,用于如果所述采样点所在的储层段的裂缝发育因子不小于所述第一预设值且不大于第二预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段中的次裂缝发育区;
第三确定子模块,用于如果所述采样点所在的储层段的裂缝发育因子大于所述第二预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段中的裂缝发育区。
本发明实施例中,获取火成岩储层的多个采样点的测井数据后,由于测井数据可以反映岩石的总孔隙度,而岩石的总孔隙度能够反映裂缝发育程度,因此,可以基于该多个采样点中每个采样点的测井数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子。之后,对于该多个采样点中的每个采样点,如果该采样点所在的储层段的裂缝发育因子不小于第一预设值,则可以确定该储层段为火成岩储层的裂缝段。由于多个采样点的测井数据能够简单方便地获取得到,且本发明实施例中仅基于该多个采样点的测井数据,就能完成对火成岩储层的裂缝段的识别,因此,将大大节省火成岩储层的裂缝段的识别成本,且提高识别效率。
附图说明
为了更清楚地说明本发明实施例中的技术方案,下面将对实施例描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本发明的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动的前提下,还可以基于这些附图获得其他的附图。
图1是本发明实施例提供的一种火成岩储层的裂缝段的识别方法的流程图;
图2A是本发明实施例提供的另一种火成岩储层的裂缝段的识别方法的流程图;
图2B是本发明实施例提供的多个测井数据对应的三条测井曲线的示意图;
图2C是本发明实施例提供的确定储层段的裂缝发育因子的操作的流程图;
图3A是本发明实施例提供的第一种火成岩储层的裂缝段的识别装置的结构示意图;
图3B是本发明实施例提供的第二种火成岩储层的裂缝段的识别装置的结构示意图;
图3C是本发明实施例提供的第三种火成岩储层的裂缝段的识别装置的结构示意图;
图3D是本发明实施例提供的第四种火成岩储层的裂缝段的识别装置的结构示意图。
具体实施方式
为使本发明的目的、技术方案和优点更加清楚,下面将结合附图对本发明实施方式作进一步地详细描述。
图1是本发明实施例提供的一种火成岩储层的裂缝段的识别方法的流程图。参见图1,该方法主要包括如下步骤。
步骤101:获取油气井的火成岩储层中的多个采样点的测井数据,该测井数据包括密度、补偿中子和深侧向电阻率。
步骤102:基于多个采样点的测井数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子。
步骤103:对于多个采样点中的每个采样点,如果采样点所在的储层段的裂缝发育因子不小于第一预设值,则确定采样点所在的储层段为该火成岩储层的裂缝段。
本发明实施例中,获取火成岩储层的多个采样点的测井数据后,由于测井数据可以反映岩石的总孔隙度,而岩石的总孔隙度能够反映裂缝发育程度,因此,可以基于该多个采样点中每个采样点的测井数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子。之后,对于该多个采样点中的每个采样点,如果该采样点所在的储层段的裂缝发育因子不小于第一预设值,则可以确定该储层段为火成岩储层的裂缝段。由于多个采样点的测井数据能够简单方便地获取得到,且本发明实施例中仅基于该多个采样点的测井数据,就能完成对火成岩储层的裂缝段的识别的方案,因此,将大大节省火成岩储层的裂缝段的识别成本,且提高识别效率。
可选地,基于多个采样点的测井数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子,包括:
将多个采样点的测井数据划分为第一组测井数据、第二组测井数据和第三组测井数据,第一组测井数据为该多个采样点的密度,第二组测井数据为该多个采样点的补偿中子,第三组测井数据为该多个采样点的深侧向电阻率;
对第一组测井数据进行归一化处理,得到第一组归一化数据,对第二组测井数据进行归一化处理,得到第二组归一化数据,对第三组测井数据进行归一化处理,得到第三组归一化数据;
基于第一组归一化数据、第二组归一化数据和第三组归一化数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子。
可选地,基于第一组归一化数据、第二组归一化数据和第三组归一化数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子,包括:
分别确定第一组归一化数据的平均值、第二组归一化数据的平均值和第三组归一化数据的平均值;
对于该多个采样点中的每个采样点,分别从第一组归一化数据、第二组归一化数据和第三组归一化数据中获取该采样点的测井数据的归一化数据;
基于第一组归一化数据的平均值、第二组归一化数据的平均值、第三组归一化数据的平均值和采样点的测井数据的归一化数据,确定采样点所在的储层段的裂缝发育因子。
可选地,基于第一组归一化数据的平均值、第二组归一化数据的平均值、第三组归一化数据的平均值和采样点的测井数据的归一化数据,确定采样点所在的储层段的裂缝发育因子,包括:
基于第一组归一化数据的平均值、第二组归一化数据的平均值、第三组归一化数据的平均值和采样点的测井数据的归一化数据,通过如下第一指定公式确定采样点所在的储层段的裂缝发育因子;
第一指定公式:
Figure BDA0001348634010000071
其中,K为该采样点所在的储层段的裂缝发育因子,DEN为该采样点的密度的归一化数据,CNL是为采样点的补偿中子的归一化数据,
Figure BDA0001348634010000072
为第一组归一化数据的平均值,
Figure BDA0001348634010000073
为第二组归一化数据的平均值,
Figure BDA0001348634010000074
为第三组归一化数据的平均值。
可选地,如果采样点所在的储层段的裂缝发育因子不小于第一预设值,则确定该采样点所在的储层段为火成岩储层的裂缝段,包括:
如果该采样点所在的储层段的裂缝发育因子不小于第一预设值且不大于第二预设值,则确定该采样点所在的储层段为火成岩储层的裂缝段中的次裂缝发育区;
如果该采样点所在的储层段的裂缝发育因子大于第二预设值,则确定该采样点所在的储层段为火成岩储层的裂缝段中的裂缝发育区。
上述所有可选技术方案,均可按照任意结合形成本发明的可选实施例,本发明实施例对此不再一一赘述。
图2A是本发明实施例提供的一种火成岩储层的裂缝段的识别方法的流程图。参见图2A,该方法主要包括如下步骤。
步骤201:获取油气井的火成岩储层中的多个采样点的测井数据,该测井数据包括密度、补偿中子和深侧向电阻率。
具体地,可以对油气井进行测井,得到油气井的多个测井数据;基于该多个测井数据,确定油气井的火成岩储层的位置;基于火成岩储层的位置,从该多个测井数据中获取多个采样点的测井数据。
需要说明的是,在油气藏的开发过程中,将油气井钻至预设的深度后,为了了解油气井的井下情况,通常会对油气井进行测井,以得到油气井的测井资料,该测井资料中一般会包括密度、补偿中子、深侧向电阻率等。也即是,实际应用中,该多个采样点的测井数据无需通过较为复杂的操作(如取芯)得到,而仅通过常规的测井就可以简单方便的采集得到。
另外,多个采样点位于火成岩储层中,且该多个采样点可以为火成岩储层中的多个测井点,或者为火成岩储层中的多个预设深度处的位置点,该多个预设深度可以预先进行设置,本发明实施例对此不作限定。
再者,密度用于指示岩石的体积密度,主要反映岩石的总孔隙度,密度越小,岩石的裂缝发育程度越高;补偿中子用于指示岩石的总孔隙度,补偿中子指示的总孔隙度值越高,岩石的裂缝发育程度越高;深侧向电阻率用于反映岩石的总孔隙度,深侧向电阻率越大,岩石的裂缝发育程度越低,深侧向电阻率越小,岩石的裂缝发育程度越高。
其中,对油气井进行测井,得到油气井的多个测井数据的操作可以为:通过多个测井设备对油气井进行测井,得到油气井的测井资料,之后,从该测井资料中获取多个测井数据。当然,也可以通过其它方式对油气井进行测井,得到油气井的多个测井数据,本发明实施例对此不作限定。
需要说明的是,多个测井设备可以包括密度测井仪、补偿中子测井仪、深侧向测井仪等。
其中,基于该多个测井数据,确定油气井的火成岩储层的位置的操作可以为:确定该多个测井数据对应的测井曲线;基于该多个测井数据对应的测井曲线,通过指定分析法确定油气井的火成岩储层的位置。当然,也可以通过其它方式基于该多个测井数据,确定油气井的火成岩储层的位置,本发明实施例对此不作限定。
需要说明的是,指定分析法可以预先进行设置,如指定分析法可以为测井相自动分析法等。
其中,基于火成岩储层的位置,从该多个测井数据中获取多个采样点的测井数据的操作可以为:如果该多个采样点为火成岩储层中的多个测井点,则从该多个测井数据中获取在该火成岩储层的位置处测得的测井数据,并将获取的测井数据中在同一深度测得的测井数据确定为同一采样点的测井数据;如果该多个采样点为火成岩储层中的多个预设深度处的位置点,则确定该多个测井数据对应的三条测井曲线,从这三条测井曲线中分别获取多个预设深度对应的测井数据,并将获取的测井数据中对应同一预设深度的测井数据确定为同一采样点的测井数据。
例如,如果该多个采样点为火成岩储层中的多个测井点,假设火成岩储层的位置为2760米~2810米,则可以从该多个测井数据中获取在2760米~2810米处测得的测井数据,并将2760米~2810米处测得的测井数据中在同一深度测得的测井数据确定为同一采样点的测井数据。
再例如,如果该多个采样点为火成岩储层中的多个预设深度处的位置点,则可以确定该多个测井数据对应的三条测井曲线。如图2B所示,假设这三条测井曲线分别为测井曲线1、测井曲线2和测井曲线3,且测井曲线1、测井曲线2、测井曲线3与该多个测井数据中的密度、补偿中子、深侧向电阻率一一对应,且假设多个预设深度为2770米和2780米,则可以从这三条测井曲线中分别获取2770米和2780米对应的测井数据,并将2770米对应的测井数据确定为一个采样点的测井数据,将2780米对应的测井数据确定为另一个采样点的测井数据。
步骤202:基于该多个采样点的测井数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子。
需要说明的是,对于该多个采样点中的每个采样点,该采样点所在的储层段的裂缝发育因子用于指示该储层段的裂缝发育程度,也即是,如果该储层段的裂缝发育因子越大,则表明该储层段的裂缝发育程度越高,如果该储层段的裂缝发育因子越小,则表明该储层段的裂缝发育程度越低。
另外,由于采样点的密度、补偿中子和深侧向电阻率能够反映该采样点处的岩石的总孔隙度,而岩石的总孔隙度能够反映储层段的裂缝发育程度,因此,可以基于多个采样点的测井数据,来确定该多个采样点中每个采样点所在的储层段的裂缝发育因子。
由于该多个采样点的测井数据的地质因素、响应机理、量纲和采集装置等存在差异,所以基于该多个采样点的测井数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子时,需要先对该多个采样点的测井数据进行归一化处理,得到该多个采样点的测井数据的归一化数据,再基于该多个采样点的测井数据的归一化数据来确定每个采样点所在的储层段的裂缝发育因子。因此,参见图2C,步骤202可以包括如下步骤2021-步骤2023。
步骤2021:将该多个采样点的测井数据划分为第一组测井数据、第二组测井数据和第三组测井数据,第一组测井数据为该多个采样点的密度,第二组测井数据为该多个采样点的补偿中子,第三组测井数据为该多个采样点的深侧向电阻率。
由于该多个采样点的测井数据包括密度、补偿中子和深侧向电阻率,因此,可以按照密度、补偿中子和深侧向电阻率,将该多个采样点的测井数据分为三组测井数据,以便可以基于这三组测井数据在后续的步骤2022中实现对该多个采样点的测井数据的归一化处理。
例如,该多个采样点包括采样点1和采样点2,且采样点1的测井数据包括密度1、补偿中子1和深侧向电阻率1,采样点2的测井数据包括密度2、补偿中子2和深侧向电阻率2,此时可以将该多个采样点的测井数据划分为第一组测井数据、第二组测井数据和第三组测井数据,且第一组测井数据为密度1和密度2,第二组测井数据为补偿中子1和补偿中子2,第三组测井数据为深侧向电阻率1和深侧向电阻率2。
步骤2022:对第一组测井数据进行归一化处理,得到第一组归一化数据;对第二组测井数据进行归一化处理,得到第二组归一化数据;对第三组测井数据进行归一化处理,得到第三组归一化数据。
需要说明的是,第一组归一化数据、第二组归一化数据和第三组归一化数据即为该多个采样点的测井数据的归一化数据,且第一组归一化数据为该多个采样点的密度的归一化数据,第二组归一化数据为该多个采样点的补偿中子的归一化数据,第三组归一化数据为该多个采样点的深侧向电阻率的归一化数据。
其中,对第一组测井数据进行归一化处理,得到第一组归一化数据的操作可以为:对于该多个采样点的测井数据中的每个密度,基于第一组测井数据,通过如下第二指定公式确定该密度的归一化数据;将该多个采样点的测井数据中每个密度的归一化数据确定为第一组归一化数据;
第二指定公式:
Figure BDA0001348634010000101
其中,X'为该密度的归一化数据,X为该密度,Xmin为第一组测井数据中最小的测井数据,Xmax为第一组测井数据中最大的测井数据。
当然,实际应用中,也可以通过其它方式对第一组测井数据进行归一化处理,得到第一组归一化数据,本发明实施例对此不作限定。
其中,对第二组测井数据进行归一化处理,得到第二组归一化数据的操作可以为:对于该多个采样点的测井数据中的每个补偿中子,基于第二组测井数据,通过如下第三指定公式确定该补偿中子的归一化数据;将该多个采样点的测井数据中每个补偿中子的归一化数据确定为第二组归一化数据;
第三指定公式:
Figure BDA0001348634010000111
其中,Y'为该补偿中子的归一化数据,Y为该补偿中子,Ymin为第二组测井数据中最小的测井数据,Ymax为第二组测井数据中最大的测井数据。
当然,实际应用中,也可以通过其它方式对第二组测井数据进行归一化处理,得到第二组归一化数据,本发明实施例对此不作限定。
其中,对第三组测井数据进行归一化处理,得到第三组归一化数据的操作与上述步骤2022中对第一组测井数据进行归一化处理,得到第一组归一化数据的操作类似,本发明实施例对此不再赘述。
步骤2023:基于第一组归一化数据、第二组归一化数据和第三组归一化数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子。
具体地,可以分别确定第一组归一化数据的平均值、第二组归一化数据的平均值和第三组归一化数据的平均值;对于该多个采样点中的每个采样点,分别从第一组归一化数据、第二组归一化数据和第三组归一化数据中获取该采样点的测井数据的归一化数据;基于第一组归一化数据的平均值、第二组归一化数据的平均值、第三组归一化数据的平均值和该采样点的测井数据的归一化数据,确定该采样点所在的储层段的裂缝发育因子。
当然,实际应用中,也可以基于第一组归一化数据、第二组归一化数据和第三组归一化数据,通过其它方式确定该多个采样点中每个采样点所在的储层段的裂缝发育因子,本发明实施例对此不作限定。
其中,基于第一组归一化数据的平均值、第二组归一化数据的平均值、第三组归一化数据的平均值和该采样点的测井数据的归一化数据,确定该采样点所在的储层段的裂缝发育因子的操作可以为:基于第一组归一化数据的平均值、第二组归一化数据的平均值、第三组归一化数据的平均值和该采样点的测井数据的归一化数据,通过如下第一指定公式确定该采样点所在的储层段的裂缝发育因子;
第一指定公式:
Figure BDA0001348634010000121
其中,K为该采样点所在的储层段的裂缝发育因子,DEN为该采样点的密度的归一化数据,CNL是为该采样点的补偿中子的归一化数据,
Figure BDA0001348634010000122
为第一组归一化数据的平均值,
Figure BDA0001348634010000123
为第二组归一化数据的平均值,
Figure BDA0001348634010000124
为第三组归一化数据的平均值。
当然,实际应用中,也可以基于第一组归一化数据的平均值、第二组归一化数据的平均值、第三组归一化数据的平均值和该采样点的测井数据的归一化数据,通过其它方式确定该采样点所在的储层段的裂缝发育因子,本发明实施例对此不作限定。
步骤203:对于多个采样点中的每个采样点,判断该采样点所在的储层段的裂缝发育因子是否小于第一预设值,如果否,则执行步骤204,如果是,则执行步骤205。
需要说明的是,第一预设值可以预先进行设置,如第一预设值可以为0.005等。
另外,当该采样点所在的储层段的裂缝发育因子不小于第一预设值时,表明该储层段的裂缝发育因子较大,此时该储层段的裂缝发育程度较高;当该采样点所在的储层段的裂缝发育因子小于第一预设值时,表明该储层段的裂缝发育因子较小,此时该储层段的裂缝发育程度较低。
步骤204:确定该采样点所在的储层段为火成岩储层的裂缝段。
当该采样点所在的储层段的裂缝发育因子不小于第一预设值时,表明该储层段的裂缝发育程度较高,因此,可以确定该储层段为火成岩储层的裂缝段。
具体地,如果该采样点所在的储层段的裂缝发育因子不小于第一预设值且不大于第二预设值,则确定该采样点所在的储层段为火成岩储层的裂缝段中的次裂缝发育区;如果该采样点所在的储层段的裂缝发育因子大于第二预设值,则确定该采样点所在的储层段为火成岩储层的裂缝段中的裂缝发育区。
需要说明的是,第二预设值可以预先进行设置,如第二预设值可以为0.0065等。
如果该采样点所在的储层段的裂缝发育因子不小于第一预设值且不大于第二预设值,则表明该储层段的裂缝发育程度较高,但尚未达到很高的程度,因此,可以确定该储层段为火成岩储层的裂缝段中的次裂缝发育区;如果该采样点所在的储层段的裂缝发育因子大于第二预设值,则表明该储层段的裂缝发育程度很高,因此,可以确定该储层段为火成岩储层的裂缝段中的裂缝发育区。
步骤205:确定该采样点所在的储层段为火成岩储层的致密段。
如果该采样点所在的储层段的裂缝发育因子小于第一预设值,则表明该储层段的裂缝发育程度较低,则可以确定该储层段不为火成岩储层的裂缝段,而是为火成岩储层的致密段。
本发明实施例中,获取火成岩储层的多个采样点的测井数据后,由于测井数据可以反映岩石的总孔隙度,而岩石的总孔隙度能够反映裂缝发育程度,因此,可以基于该多个采样点中每个采样点的测井数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子。之后,对于该多个采样点中的每个采样点,如果该采样点所在的储层段的裂缝发育因子不小于第一预设值,则可以确定该储层段为火成岩储层的裂缝段,如果该采样点所在的储层段的裂缝发育因子小于第一预设值,则可以确定该储层段为火成岩储层的致密段。由于多个采样点的测井数据能够简单方便地获取得到,且本发明实施例中仅基于该多个采样点的测井数据,就能完成对火成岩储层的裂缝段和致密段的识别的方案,因此,将大大节省火成岩储层的裂缝段和致密段的识别成本,且提高识别效率。
图3A是本发明实施例提供的一种火成岩储层的裂缝段的识别装置的结构示意图。参见图3A,该装置包括:
获取模块301,用于获取油气井的火成岩储层中的多个采样点的测井数据,该测井数据包括密度、补偿中子和深侧向电阻率。
第一确定模块302,用于基于多个采样点的测井数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子。
第二确定模块303,用于对于该多个采样点中的每个采样点,如果该采样点所在的储层段的裂缝发育因子不小于第一预设值,则确定该采样点所在的储层段为该火成岩储层的裂缝段。
可选地,如图3B所示,第一确定模块302包括:
划分子模块3021,用于将多个采样点的测井数据划分为第一组测井数据、第二组测井数据和第三组测井数据,第一组测井数据为多个采样点的密度,第二组测井数据为多个采样点的补偿中子,第三组测井数据为多个采样点的深侧向电阻率。
处理子模块3022,用于对第一组测井数据进行归一化处理,得到第一组归一化数据,对第二组测井数据进行归一化处理,得到第二组归一化数据,对第三组测井数据进行归一化处理,得到第三组归一化数据。
第一确定子模块3023,用于基于第一组归一化数据、第二组归一化数据和第三组归一化数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子。
可选地,如图3C所示,第一确定子模块3023包括:
第一确定单元30231,用于分别确定第一组归一化数据的平均值、第二组归一化数据的平均值和第三组归一化数据的平均值。
获取单元30232,用于对于该多个采样点中的每个采样点,分别从第一组归一化数据、第二组归一化数据和第三组归一化数据中获取采样点的测井数据的归一化数据。
第二确定单元30233,用于基于第一组归一化数据的平均值、第二组归一化数据的平均值、第三组归一化数据的平均值和采样点的测井数据的归一化数据,确定该采样点所在的储层段的裂缝发育因子。
可选地,第二确定单元30233主要用于:
基于第一组归一化数据的平均值、第二组归一化数据的平均值、第三组归一化数据的平均值和采样点的测井数据的归一化数据,通过如下第一指定公式确定该采样点所在的储层段的裂缝发育因子;
第一指定公式:
Figure BDA0001348634010000141
其中,K为该采样点所在的储层段的裂缝发育因子,DEN为该采样点的密度的归一化数据,CNL是为该采样点的补偿中子的归一化数据,
Figure BDA0001348634010000142
为第一组归一化数据的平均值,
Figure BDA0001348634010000143
为第二组归一化数据的平均值,
Figure BDA0001348634010000144
为第三组归一化数据的平均值。
可选地,如图3D所示,第二确定模块303包括:
第二确定子模块3031,用于如果该采样点所在的储层段的裂缝发育因子不小于第一预设值且不大于第二预设值,则确定该采样点所在的储层段为火成岩储层的裂缝段中的次裂缝发育区。
第三确定子模块3032,用于如果该采样点所在的储层段的裂缝发育因子大于第二预设值,则确定该采样点所在的储层段为火成岩储层的裂缝段中的裂缝发育区。
本发明实施例中,获取火成岩储层的多个采样点的测井数据后,由于测井数据可以反映岩石的总孔隙度,而岩石的总孔隙度能够反映裂缝发育程度,因此,可以基于该多个采样点中每个采样点的测井数据,确定该多个采样点中每个采样点所在的储层段的裂缝发育因子。之后,对于该多个采样点中的每个采样点,如果该采样点所在的储层段的裂缝发育因子不小于第一预设值,则可以确定该储层段为火成岩储层的裂缝段,如果该采样点所在的储层段的裂缝发育因子小于第一预设值,则可以确定该储层段为火成岩储层的致密段。由于多个采样点的测井数据能够简单方便地获取得到,且本发明实施例中仅基于该多个采样点的测井数据,就能完成对火成岩储层的裂缝段和致密段的识别的方案,因此,将大大节省火成岩储层的裂缝段和致密段的识别成本,提高识别效率。
需要说明的是:上述实施例提供的火成岩储层的裂缝段的识别方法在火成岩储层的裂缝段识别时,仅以上述各功能模块的划分进行举例说明,实际应用中,可以基于需要而将上述功能分配由不同的功能模块完成,即将设备的内部结构划分成不同的功能模块,以完成以上描述的全部或者部分功能。另外,上述实施例提供的火成岩储层的裂缝段的识别装置与火成岩储层的裂缝段的识别方法实施例属于同一构思,其具体实现过程详见方法实施例,这里不再赘述。
本领域普通技术人员可以理解实现上述实施例的全部或部分步骤可以通过硬件来完成,也可以通过程序来指令相关的硬件完成,所述的程序可以存储于一种计算机可读存储介质中,上述提到的存储介质可以是只读存储器,磁盘或光盘等。
以上所述仅为本发明的较佳实施例,并不用以限制本发明,凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。

Claims (4)

1.一种火成岩储层的裂缝段的识别方法,其特征在于,所述方法包括:
获取油气井的火成岩储层中的多个采样点的测井数据,所述测井数据包括密度、补偿中子和深侧向电阻率;
将所述多个采样点的测井数据划分为第一组测井数据、第二组测井数据和第三组测井数据,所述第一组测井数据为所述多个采样点的密度,所述第二组测井数据为所述多个采样点的补偿中子,所述第三组测井数据为所述多个采样点的深侧向电阻率;
对所述第一组测井数据进行归一化处理,得到第一组归一化数据;对所述第二组测井数据进行归一化处理,得到第二组归一化数据;对所述第三组测井数据进行归一化处理,得到第三组归一化数据;
分别确定所述第一组归一化数据的平均值、所述第二组归一化数据的平均值和所述第三组归一化数据的平均值;
对于所述多个采样点中的每个采样点,分别从所述第一组归一化数据、所述第二组归一化数据和所述第三组归一化数据中获取所述采样点的测井数据的归一化数据;
基于所述第一组归一化数据的平均值、所述第二组归一化数据的平均值、所述第三组归一化数据的平均值和所述采样点的测井数据的归一化数据,通过如下第一指定公式确定所述采样点所在的储层段的裂缝发育因子,所述裂缝发育因子用于指示裂缝发育程度;
第一指定公式:
Figure FDA0002598281200000011
其中,K为所述采样点所在的储层段的裂缝发育因子,DEN为所述采样点的密度的归一化数据,CNL为所述采样点的补偿中子的归一化数据,
Figure FDA0002598281200000012
为所述第一组归一化数据的平均值,
Figure FDA0002598281200000013
为所述第二组归一化数据的平均值,
Figure FDA0002598281200000014
为所述第三组归一化数据的平均值;
对于所述多个采样点中的每个采样点,如果所述采样点所在的储层段的裂缝发育因子不小于第一预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段,所述裂缝段包括裂缝发育区和次裂缝发育区。
2.如权利要求1所述的方法,其特征在于,所述如果所述采样点所在的储层段的裂缝发育因子不小于第一预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段,包括:
如果所述采样点所在的储层段的裂缝发育因子不小于所述第一预设值且不大于第二预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段中的次裂缝发育区;
如果所述采样点所在的储层段的裂缝发育因子大于所述第二预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段中的裂缝发育区。
3.一种火成岩储层的裂缝段的识别装置,其特征在于,所述装置包括:
获取模块,用于获取油气井的火成岩储层中的多个采样点的测井数据,所述测井数据包括密度、补偿中子和深侧向电阻率;
第一确定模块,用于基于所述多个采样点的测井数据,确定所述多个采样点中每个采样点所在的储层段的裂缝发育因子,所述裂缝发育因子用于指示裂缝发育程度;
第二确定模块,用于对于所述多个采样点中的每个采样点,如果所述采样点所在的储层段的裂缝发育因子不小于第一预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段,所述裂缝段包括裂缝发育区和次裂缝发育区;
其中,所述第一确定模块包括:
划分子模块,用于将所述多个采样点的测井数据划分为第一组测井数据、第二组测井数据和第三组测井数据,所述第一组测井数据为所述多个采样点的密度,所述第二组测井数据为所述多个采样点的补偿中子,所述第三组测井数据为所述多个采样点的深侧向电阻率;
处理子模块,用于对所述第一组测井数据进行归一化处理,得到第一组归一化数据,对所述第二组测井数据进行归一化处理,得到第二组归一化数据,对所述第三组测井数据进行归一化处理,得到第三组归一化数据;
第一确定子模块,用于基于所述第一组归一化数据、所述第二组归一化数据和所述第三组归一化数据,确定所述多个采样点中每个采样点所在的储层段的裂缝发育因子;
其中,所述第一确定子模块包括:
第一确定单元,用于分别确定所述第一组归一化数据的平均值、所述第二组归一化数据的平均值和所述第三组归一化数据的平均值;
获取单元,用于对于所述多个采样点中的每个采样点,分别从所述第一组归一化数据、所述第二组归一化数据和所述第三组归一化数据中获取所述采样点的测井数据的归一化数据;
第二确定单元,用于基于所述第一组归一化数据的平均值、所述第二组归一化数据的平均值、所述第三组归一化数据的平均值和所述采样点的测井数据的归一化数据,确定所述采样点所在的储层段的裂缝发育因子;
其中,所述第二确定单元主要用于:
基于所述第一组归一化数据的平均值、所述第二组归一化数据的平均值、所述第三组归一化数据的平均值和所述采样点的测井数据的归一化数据,通过如下第一指定公式确定所述采样点所在的储层段的裂缝发育因子;
第一指定公式:
Figure FDA0002598281200000031
其中,K为所述采样点所在的储层段的裂缝发育因子,DEN为所述采样点的密度的归一化数据,CNL为所述采样点的补偿中子的归一化数据,
Figure FDA0002598281200000032
为所述第一组归一化数据的平均值,
Figure FDA0002598281200000033
为所述第二组归一化数据的平均值,
Figure FDA0002598281200000034
为所述第三组归一化数据的平均值。
4.如权利要求3所述的装置,其特征在于,所述第二确定模块包括:
第二确定子模块,用于如果所述采样点所在的储层段的裂缝发育因子不小于所述第一预设值且不大于第二预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段中的次裂缝发育区;
第三确定子模块,用于如果所述采样点所在的储层段的裂缝发育因子大于所述第二预设值,则确定所述采样点所在的储层段为所述火成岩储层的裂缝段中的裂缝发育区。
CN201710566959.5A 2017-07-12 2017-07-12 火成岩储层的裂缝段的识别方法及装置 Active CN109254330B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710566959.5A CN109254330B (zh) 2017-07-12 2017-07-12 火成岩储层的裂缝段的识别方法及装置

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710566959.5A CN109254330B (zh) 2017-07-12 2017-07-12 火成岩储层的裂缝段的识别方法及装置

Publications (2)

Publication Number Publication Date
CN109254330A CN109254330A (zh) 2019-01-22
CN109254330B true CN109254330B (zh) 2020-11-03

Family

ID=65051484

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710566959.5A Active CN109254330B (zh) 2017-07-12 2017-07-12 火成岩储层的裂缝段的识别方法及装置

Country Status (1)

Country Link
CN (1) CN109254330B (zh)

Family Cites Families (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103149589B (zh) * 2013-02-22 2016-01-06 中国石油天然气股份有限公司 火成岩油气勘探方法及装置
CN104635269B (zh) * 2013-11-13 2017-06-16 中国石油化工股份有限公司 一种基于叠前方位角道集预测火成岩裂缝型储层的方法
CN104747163B (zh) * 2013-12-31 2017-10-17 中国石油天然气股份有限公司 一种在致密砂岩中识别储层裂缝的方法及装置
CN106772570A (zh) * 2015-11-24 2017-05-31 中国石油化工股份有限公司 火成岩开启缝地震预测方法
CN106872669B (zh) * 2017-04-24 2018-10-16 中国海洋石油集团有限公司 一种花岗岩潜山裂缝储层电性分类方法

Also Published As

Publication number Publication date
CN109254330A (zh) 2019-01-22

Similar Documents

Publication Publication Date Title
CN106323836B (zh) 一种井壁渗透率计算方法及装置
CN107402176B (zh) 裂缝孔隙度的确定方法和装置
CN109298464B (zh) 致密砂岩储层成岩相测井识别方法及装置
CN105114067A (zh) 一种岩性测井相的方法
CN108830140A (zh) 一种基于电成像测井分形维数的火山岩岩性识别方法
CN107688037A (zh) 一种利用核磁测井t2分布确定井下岩石粒度曲线的方法
AU2013392090B2 (en) Systems and methods for optimizing existing wells and designing new wells based on the distribution of average effective fracture lengths
CN108035710B (zh) 基于数据挖掘划分深层岩石地质相的方法
CN104463686A (zh) 利用判别分析法随钻识别页岩气储层的方法
CN116201535B (zh) 油气藏目标井标志地层自动划分方法、装置及设备
CN116821786A (zh) 一种低渗储层分类方法、装置及介质
CN109254330B (zh) 火成岩储层的裂缝段的识别方法及装置
CN108518220A (zh) 基于综合流体识别指数识别储层流体的方法与装置
CN109738955B (zh) 一种基于成分-结构分类下的变质岩岩性综合判别方法
US10664635B2 (en) Determining non-linear petrofacies using cross-plot partitioning
CN117388945A (zh) 一种断缝体油藏储层类型的识别方法
CN111236934B (zh) 水淹级别确定方法和装置
CN114153002A (zh) 储层天然裂缝三维地质建模方法、装置、电子设备及介质
CN111352155B (zh) 流体识别方法及装置
CN112989257A (zh) 海相页岩油气储层的产气量测定方法
CN112990320A (zh) 一种岩性的分类方法、装置、电子设备及存储介质
CN114428365A (zh) 一种基于含气特征的致密砂岩气层测井识别方法及装置
CN110633904A (zh) 一种变质岩储层分类方法及装置
CN105528732B (zh) 一种试气井产能预测方法
CN117328852B (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