一种地层缝洞类型的识别方法及装置制造方法及图纸

本发明专利技术提供一种地层缝洞类型的识别方法及装置，属于地球物理探测技术领域。该方法包括：对研究区内待测井的测井曲线进行小波多尺度变换，得到待测井的二维能量谱；将待测井的二维能量谱中的能量幅值大小以及能量幅值所分布的尺度范围，与标准数据库中的数据进行比对，根据比对结果，得到待测井的缝洞类型；所述标准数据库中的数据为研究区内各缝洞类型对应的能量幅值大小以及能量幅值所分布的尺度范围，这些数据通过先将研究区内已知缝洞类型的标准井的测井曲线进行小波多尺度变换，得到标准井的二维能量谱，然后统计标准井的二维能量谱中的能量幅值大小以及能量幅值所分布的尺度范围得到。本发明专利技术能够快速、低成本的精确识别地层缝洞类型。

【技术实现步骤摘要】
一种地层缝洞类型的识别方法及装置
本专利技术涉及一种地层缝洞类型的识别方法及装置，属于地球物理探测

技术介绍
在石油、天然气、地热等资源勘探中，人们一般以测井资料来进行地层岩性、电阻率孔隙度、渗透率、含水饱和度以及裂缝和孔洞的评价。缝洞在碳酸盐岩油气和地热储层中是一类重要的能源存储空间，目前国内外涉及碳酸盐岩缝洞测井识别的方法主要有以下几种，其原理和特点如下：(1)由于常规测井对于地层裂缝和孔洞的识别精度有限，常需要借助于电成像、声成像以及偶极子声波等特殊测井方法来完成。例如公布号为CN106370817A的专利技术专利申请文件中，公开了一种基于岩心分析和电成像测井的缝洞定量表征方法，虽然电成像对于过井眼的缝洞有着精准的识别效果，但是由于较高的测井成本或者复杂的井眼状况而常常放弃测量。(2)以电阻率、声波、中子、密度、井径等常规测井曲线进行综合识别。该类方法主要是综合利用各种常规测井方法对孔洞的响应特征，建立一种加权识别函数进而寻求对缝洞的有效识别，但是由于各种常规测井对过井眼的孔洞和钻井过程中形成的井壁扩径有相似的响应特征而常常识别失败，无法区分是由于工程原因造成的井眼扩径还是天然的缝洞造成的井眼扩径。(3)从碳酸盐岩储层导电机理入手，建立基质、裂缝以及孔洞三种孔隙度类型的数学导电方程，依据对其孔隙度的求解从而识别缝洞。但是这种方法也需要利用常规测井中的深、浅电阻率数据来进行计算，是对缝洞的一种近似估算。因此，上述缝洞识别方法要么利用各种成像技术实现，但是这种方式成本高；要么通过建立数学方程实现，但是这种方式处理的数据多，且识别结果不准确；虽然也有利用常规测井曲线进行缝洞识别的方法，但由于一些因素的影响对过井眼的缝洞无法精确识别。
技术实现思路
本专利技术的目的是提供一种地层缝洞类型的识别方法及装置，用以解决目前利用常规测井方法识别地层缝洞类型带来的识别精度低的问题。为实现上述目的，本专利技术提供了一种地层缝洞类型的识别方法，该方法包括以下步骤：对研究区内待测井的测井曲线进行小波多尺度变换，得到待测井的二维能量谱，该能量谱的横坐标为小波变换的尺度、纵坐标为测井曲线的能量幅值；将待测井的二维能量谱中的能量幅值大小以及能量幅值所分布的尺度范围，与标准数据库中的数据进行比对，根据比对结果，得到待测井的缝洞类型；所述标准数据库中的数据为研究区内各缝洞类型对应的能量幅值大小以及能量幅值所分布的尺度范围，这些数据通过先将研究区内已知缝洞类型的标准井的测井曲线进行小波多尺度变换，得到标准井的二维能量谱，然后统计标准井的二维能量谱中的能量幅值大小以及能量幅值所分布的尺度范围得到。本专利技术还提供了一种地层缝洞类型的识别装置，该装置包括处理器和存储器，所述处理器执行由所述存储器存储的计算机程序，以实现上述地层缝洞类型的识别方法。本专利技术的地层缝洞类型的识别方法及装置的有益效果是：本专利技术依据裂缝、孔洞和扩径三者各自拥有不同的二维能量分布区间，基于标准井已知的缝洞类型信息和标准井的二维能量谱，建立了包含各缝洞类型对应的能量幅值大小以及能量幅值所分布的尺度范围的标准数据库，然后，通过将待测井的二维能量谱中的能量幅值大小以及能量幅值所分布的尺度范围，与标准数据库中的数据进行比对，进行待测井的缝洞类型识别。本专利技术基于常规测井曲线即可实现地层缝洞类型识别，可以在没有电成像等特殊项目测井资料的条件下，快速、低成本的精确识别地层裂缝、孔洞以及井眼扩径情况；而且能区分出是天然缝洞导致的井眼扩径还是由于钻井工程原因导致的井眼扩径，即能够精确识别出过井眼的缝洞。进一步地，在上述方法及装置中，所述缝洞类型包括裂缝、孔洞和扩径。为了更好地实现小波多尺度变换，进一步地，在上述方法及装置中，所述小波多尺度变换的小波函数为symN小波函数或dbN小波函数。为了避免小波多尺度变换后的数据冗余，进一步地，在上述方法及装置中，所述能量幅值为进行小波多尺度变换后提取的正半周幅值。进一步地，在上述方法及装置中，所述二维能量谱用能量矩阵的方式表示，所述能量矩阵为：其中，ⅹm,n为小波多尺度变换后某点的正半周幅值；n为分解层数，即尺度变量；m为深度变量。进一步地，在上述方法及装置中，所述测井曲线为井径测井曲线。附图说明图1是方法实施例1中地层缝洞类型的识别方法流程图；图2是方法实施例1中裂缝的能量分布标准模式图；图3是方法实施例1中孔洞的能量分布标准模式图；图4是方法实施例1中扩径的能量分布标准模式图；图5是方法实施例1中E井的二维能量谱及缝洞类型识别结论图；图6是装置实施例1的结构图。具体实施方式本专利技术提供了一种基于常规测井曲线的地层缝洞类型的识别方法及装置，可以在没有电成像等特殊项目测井资料的条件下，快速、低成本的精确识别地层裂缝、孔洞以及井眼扩径情况，降低了对电成像等特殊项目测井的依赖，降低了识别成本；而且能区分出是天然缝洞导致的井眼扩径还是由于钻井工程原因导致的井眼扩径，即能够精确识别出过井眼的缝洞；另外，还可以通过建立纵向深度上连续的二维能量谱，实现在多个深度范围内识别缝洞类型，能够为油气以及地热能源勘探过程中地层缝洞类型的快速、有效识别提供技术支撑。方法实施例1如图1所示，本实施例的地层缝洞类型的识别方法，包括以下步骤：1)对研究区内待测井的测井曲线进行小波多尺度变换，得到待测井的二维能量谱。其中，二维能量谱的横坐标为小波变换的尺度、纵坐标为测井曲线的能量幅值。本实施例中，利用symN小波函数进行小波多尺度变换，过程如下：式中，f(t)为井径测井曲线信号，f(t)∈L2(R)，L2(R)表示平方可积的实数空间，a为尺度因子，b为平移因子，wf(a,b)为f(t)进行小波变换后的结果。设a＝2-k，将连续小波wf(a,b)离散化：式中，为小波离散化后的结果。设ⅹm,n为小波离散化后某点的正半周幅值；n为分解层数，即尺度变量；m为深度变量。将二维能量谱用能量矩阵的方式表示，则能量矩阵表示为：本实施例中，选用的sym小波函数具有近似对称性，变换后为了避免数据冗余，仅提取正半周幅值，舍去负半周幅值。作为其他实施方式，二维能量谱还可以用能量矩阵之外的其他方式表示，例如用平面图的形式表示。2)将待测井的二维能量谱中的能量幅值大小以及能量幅值所分布的尺度范围，与标准数据库中的数据进行比对，根据比对结果，得到待测井的缝洞类型，缝洞类型包括裂缝、孔洞和扩径。其中，标准数据库中的数据为研究区内各缝洞类型对应的能量幅值大小以及能量幅值所分布的尺度范围，这些数据通过先将研究区内已知缝洞类型的标准井的测井曲线进行小波多尺度变换，得到标准井的二维能量谱，然后统计标准井的二维能量谱中的能量幅值大小以及能量幅值所分布的尺度范围得到。下面对某研究本文档来自技高网...

【技术保护点】
1.一种地层缝洞类型的识别方法，其特征在于，该方法包括以下步骤：/n对研究区内待测井的测井曲线进行小波多尺度变换，得到待测井的二维能量谱，该能量谱的横坐标为小波变换的尺度、纵坐标为测井曲线的能量幅值；/n将待测井的二维能量谱中的能量幅值大小以及能量幅值所分布的尺度范围，与标准数据库中的数据进行比对，根据比对结果，得到待测井的缝洞类型；/n所述标准数据库中的数据为研究区内各缝洞类型对应的能量幅值大小以及能量幅值所分布的尺度范围，这些数据通过先将研究区内已知缝洞类型的标准井的测井曲线进行小波多尺度变换，得到标准井的二维能量谱，然后统计标准井的二维能量谱中的能量幅值大小以及能量幅值所分布的尺度范围得到。/n

【技术特征摘要】
1.一种地层缝洞类型的识别方法，其特征在于，该方法包括以下步骤：
对研究区内待测井的测井曲线进行小波多尺度变换，得到待测井的二维能量谱，该能量谱的横坐标为小波变换的尺度、纵坐标为测井曲线的能量幅值；
将待测井的二维能量谱中的能量幅值大小以及能量幅值所分布的尺度范围，与标准数据库中的数据进行比对，根据比对结果，得到待测井的缝洞类型；
所述标准数据库中的数据为研究区内各缝洞类型对应的能量幅值大小以及能量幅值所分布的尺度范围，这些数据通过先将研究区内已知缝洞类型的标准井的测井曲线进行小波多尺度变换，得到标准井的二维能量谱，然后统计标准井的二维能量谱中的能量幅值大小以及能量幅值所分布的尺度范围得到。


2.根据权利要求1所述的地层缝洞类型的识别方法，其特征在于，所述缝洞类型包括裂缝、孔洞和扩径。


3.根据权利要求1或2所述的地层缝洞...

【专利技术属性】
技术研发人员：李健伟，冉利民，赵永刚，田飞，王磊，周贤斌，李功强，齐真真，刘亮，刘珺，黄潮，白彬艳，于路鹏，
申请(专利权)人：中石化石油工程技术服务有限公司，中石化华北石油工程有限公司，
类型：发明
国别省市：北京;11