阵列声波测井识别和评价岩体弱面的方法、系统、设备及介质

摘要

本发明提供了一种阵列声波测井识别和评价岩体弱面的方法、系统、设备及介质。所述方法包括：利用自动增益恢复系数对采集到的阵列声波测井波形数据进行增益恢复；利用变异函数对增益恢复后的阵列声波测井波形数据进行统计变异处理；对统计变异处理后的波形数据以变密度的方式进行灰度成图，绘制得到经统计变异处理后的第一波形变密度图，第一波形变密度图上显示有岩体弱面的位置及规模；以及计算岩体弱面表征指数，并基于岩体弱面表征指数来评价岩体弱面。根据本发明，可基于阵列声波测井资料来快速、直观地识别和评价井下深部地层岩体弱面。

权利要求

1.一种阵列声波测井识别和评价岩体弱面的方法，其特征在于，所述方法包括：利用自动增益恢复系数对采集到的阵列声波测井波形数据进行增益恢复；

利用变异函数对增益恢复后的阵列声波测井波形数据进行统计变异处理；

对统计变异处理后的波形数据以变密度的方式进行灰度成图，绘制得到经统计变异处理后的第一波形变密度图，第一波形变密度图上显示有岩体弱面的位置及规模；以及计算岩体弱面表征指数，并基于岩体弱面表征指数来评价岩体弱面，其中，所述岩体弱面表征指数越大，则岩体弱面越容易出现滑移，所述计算岩体弱面表征指数具体包括：对统计变异处理后的波形数据开时窗进行声波均方根幅度计算，得到第一均方根幅度；

对相同时窗标准硬地层经统计变异处理后的波形数据开时窗进行声波均方根幅度计算，得到第二均方根幅度；和

以第一均方根幅度与第二均方根幅度的比值作为所述岩体弱面表征指数。

2.根据权利要求1所述的阵列声波测井识别和评价岩体弱面的方法，其特征在于，所述变异函数为：

式中：

Vsize为波形阵列数据的时间采样点数，不同测井系列，Vsize值不同；

为 深度点第M个时间采样点经统计变异处理技术处理后的反映岩体弱面特征波形数据；

N为计算点的深度；

N+1为计算点的下一个深度；

M为当前计算点的波形数据阵列的时间采样点，不同测井系列，M值不同；

SMOTHR为深度域内进行比较的采样点数，所述采样点数为奇数，取值为1,3,5,7,9…；

ABS表示对括弧内的WAVE(N,M)‑WAVE(N+1,M)取绝对值；

WAVE(N,M)为N深度点第M个时间采样点自动增益恢复后的阵列声波测井波形数据；

WAVE(N+1,M)为N+1深度点第M个时间采样点自动增益恢复后的阵列声波测井波形数据。

3.根据权利要求1所述的阵列声波测井识别和评价岩体弱面的方法，其特征在于，所述方法还包括：将经统计变异处理技术处理后的阵列声波测井波形数据与自动增益恢复后的阵列声波测井波形数据相加，得到总体变异波形数据，以变密度的方式进行灰度成图，绘制得到经总体统计变异处理后的第二波形变密度图。

4.根据权利要求3所述的阵列声波测井识别和评价岩体弱面的方法，其特征在于，所述方法还包括：基于所述岩体弱面表征指数绘制岩体弱面表征指数曲线，给定分类阈值或分类标准，对岩体弱面分布层段进行自动识别与岩体弱面分级评价。

5.根据权利要求4所述的阵列声波测井识别和评价岩体弱面的方法，其特征在于，所述方法还包括：

绘制岩体弱面的解释评价成果图，所述解释评价成果图上包括第一曲线道、第二曲线道、第三曲线道、第四曲线道、第五曲线道、第六曲线道、第七曲线道以及岩体弱面分级评价道中的一个或多个，

其中，第一曲线道上绘制有井径、无铀伽玛和自然伽玛曲线；

第二曲线道上绘制有纵波时差、补偿中子和补偿密度曲线；

第三曲线道上绘制有深侧向和浅侧向曲线；

第四曲线道上绘制有经过所述增益恢复后的第三波形变密度图；

第五曲线道上绘制有所述第一波形变密度图；

第六曲线道上绘制有所述第二波形变密度图；

第七曲线道上绘制有所述岩体弱面表征指数曲线；

岩体弱面分级评价道上按不同图例绘制了对应深度段岩体弱面的分级评价结果。

6.根据权利要求1所述的阵列声波测井识别和评价岩体弱面的方法，其特征在于，所述方法还包括：

对阵列声波测井波形数据进行预处理，所述预处理包括数据格式解编与数据加载、测井深度校正、阵列声波测井波形数据质量检查，其中，利用阵列声波波形均方根振幅曲线进行阵列声波测井波形数据质量检查，以甄别包括设备故障造成的错误记录波形数据，选择期望的波形数据进行波形增益恢复和后续处理。

7.一种阵列声波测井识别和评价岩体弱面的系统，其特征在于，所述系统包括：增益恢复单元，被配置为利用自动增益恢复系数对采集到的阵列声波测井波形数据进行增益恢复；

统计变异处理单元，被配置为利用变异函数对增益恢复后的阵列声波测井波形数据进行统计变异处理；

灰度成图单元，被配置为对统计变异处理后的波形数据以变密度的方式进行灰度成图，绘制得到经统计变异处理后的第一波形变密度图，以显示岩体弱面的位置及规模；以及岩体弱面表征指数计算单元，被配置为：对统计变异处理后的波形数据开时窗进行声波均方根幅度计算，得到第一均方根幅度；对相同时窗标准硬地层统计变异处理后的波形数据开时窗进行声波均方根幅度计算，得到第二均方根幅度；以第一均方根幅度与第二均方根幅度的比值作为所述岩体弱面表征指数，其中，所述岩体弱面表征指数用于评价岩体弱面，所述岩体弱面表征指数越大，则岩体弱面越容易出现滑移。

8.根据权利要求7所述的阵列声波测井识别和评价岩体弱面的系统，其特征在于，所述变异函数为：

式中：

Vsize为波形阵列数据的时间采样点数，不同测井系列，Vsize值不同；

为 深度点第M个时间采样点经统计变异处理技术处理后的反映岩体弱面特征波形数据；

N为计算点的深度；

N+1为计算点的下一个深度；

M为当前计算点的波形数据阵列的时间采样点，不同测井系列，M值不同；

SMOTHR为深度域内进行比较的采样点数，所述采样点数为奇数，取值为1,3,5,7,9…；

ABS表示对括弧内的WAVE(N,M)‑WAVE(N+1,M)取绝对值；

WAVE(N,M)为N深度点第M个时间采样点自动增益恢复后的阵列声波测井波形数据；

WAVE(N+1,M)为N+1深度点第M个时间采样点自动增益恢复后的阵列声波测井波形数据。

9.根据权利要求7所述的阵列声波测井识别和评价岩体弱面的系统，其特征在于，所述系统还包括：总体变异波形变密度生成单元，被配置为将经统计变异处理技术处理后的阵列声波测井波形数据与自动增益恢复后的阵列声波测井波形数据相加，得到总体变异波形数据，以变密度的方式进行灰度成图，绘制得到经总体统计变异处理后的第二波形变密度图。

10.根据权利要求9所述的阵列声波测井识别和评价岩体弱面的系统，其特征在于，所述系统还包括：岩体弱面表征指数分级评价单元，被配置为基于所述岩体弱面表征指数绘制岩体弱面表征指数曲线，给定分类阈值或分类标准，对岩体弱面分布层段进行自动识别与岩体弱面分级评价。

11.根据权利要求10所述的阵列声波测井识别和评价岩体弱面的系统，其特征在于，所述系统还包括：

岩体弱面的解释评价成果图形成单元，用于绘制岩体弱面的解释评价成果图，所述解释评价成果图上包括第一曲线道、第二曲线道、第三曲线道、第四曲线道、第五曲线道、第六曲线道、第七曲线道以及岩体弱面分级评价道中的一个或多个，其中，第一曲线道上绘制有井径、无铀伽玛和自然伽玛曲线；

第二曲线道上绘制有纵波时差、补偿中子和补偿密度曲线；

第三曲线道上绘制有深侧向和浅侧向曲线；

第四曲线道上绘制有经过所述增益恢复后的第三波形变密度图；

第五曲线道上绘制有所述第一波形变密度图；

第六曲线道上绘制有所述第二波形变密度图；

第七曲线道上绘制有所述岩体弱面表征指数曲线；

岩体弱面分级评价道上按不同图例绘制了对应深度段岩体弱面的分级评价结果。

12.根据权利要求7所述的阵列声波测井识别和评价岩体弱面的系统，其特征在于，所述系统还包括：预处理模块，被配置为对阵列声波测井波形数据进行预处理，所述预处理包括数据格式解编与数据加载、测井深度校正、阵列声波测井波形数据质量检查，其中，利用阵列声波波形均方根振幅曲线进行阵列声波测井波形数据质量检查，以甄别包含设备故障造成的错误记录波形数据，选择期望的波形数据进行波形增益恢复和后续处理。

13.一种电子设备，其特征在于，所述电子设备包括：处理器；和

存储器，存储有计算机程序，当所述计算机程序被处理器执行时，实现如权利要求1‑6中任一项所述的方法。

14.一种存储有计算机程序的计算机可读存储介质，其特征在于，当所述计算机程序在被处理器执行时实现如权利要求1‑6中任一项所述的方法。

说明书

阵列声波测井识别和评价岩体弱面的方法、系统、设备及介质

技术领域

[0001] 本发明涉及阵列声波测井资料处理技术领域，特别地，涉及一种阵列声波测井识别和评价岩体弱面的方法、阵列声波测井识别和评价岩体弱面系统、电子设备以及计算机
可读存储介质。

背景技术

[0002] 工程岩石力学中，岩体工程的安危受岩体中存在的弱面（亦称结构面）控制或者受岩块（亦称结构体）和弱面的共同控制。弱面是岩体内存在的各种地质界面，包括：不连续面
和物质分异面，如：裂隙、节理、层理、软弱夹层、断层、断裂破碎带等。岩石中弱面的识别与
评价是工程岩石力学的重要内容。

[0003] 传统上人们对地面岩石或取心岩石弱面的研究，通常采用罗盘、尺子等工具到现场进行测量，然后对数据进行统计分析。这种方法既费力气又费时间，而且对于一些难以进
入的场地，获取数据是非常困难的，甚至是危险的。目前，采用先进的ShapeMetrix 3D数字
摄影测量系统或地面三维激光扫描(TLS)技术等，可对地面或现场岩体表面大量结构面进
行测量，通过图像处理等技术得到结构面详细信息，从而为节理、裂隙信息的快速获取开辟
新的途径，解决了以人工现场接触测量为主、劳动强度大、效率低下的难题。

[0004] 在石油行业，深部岩体弱面（特别是裂缝）对油气的勘探开发起着越来越重要的作用，人们对它的研究已有近百年历史，研究的方法有物探、测井、数理统计、非线性和力学等
多种手段。但是，对于深部岩体弱面的预测与评价迄今仍然是一个公认的世界级难题。近年
来，页岩气等非常规油气的水平井开发中，大液量、大砂量、高排量、高泵压的非常规储层体
积压裂改造作业方式，在提高油气井产量的同时，也诱发了井下套管变形问题，且套管变形
比例呈明显上升的趋势，严重影响页岩气等非常规油气的单井产量和压裂施工作业时效。
不但造成经济的巨大损失，而且也给油气井安全生产、环境保护等带来较大隐患。因此，对
岩体弱面的研究和评价就显得尤为重要。

发明内容

[0005] 本发明的目的在于解决现有技术存在的上述不足中的至少一项。例如，本发明的目的之一在于提供一种基于阵列声波测井资料来识别和评价井下深部地层岩体弱面的方
法及系统。

[0006] 为了实现上述目的，本发明提供了一种阵列声波测井识别和评价岩体弱面的方法，所述方法包括：利用自动增益恢复系数对采集到的阵列声波测井波形数据进行增益恢
复；利用变异函数对增益恢复后的阵列声波测井波形数据进行统计变异处理；对统计变异
处理后的波形数据以变密度的方式进行灰度成图，绘制得到经统计变异处理后的第一波形
变密度图，第一波形变密度图上显示有岩体弱面的位置及规模；以及计算岩体弱面表征指
数，并基于岩体弱面表征指数来评价岩体弱面，其中，所述岩体弱面表征指数越大，则岩体
弱面越容易出现滑移，所述计算岩体弱面表征指数具体包括：对统计变异处理后的波形数
据开时窗进行声波均方根幅度计算，得到第一均方根幅度；对相同时窗标准硬地层统计变
异处理后的波形数据开时窗进行声波均方根幅度计算，得到第二均方根幅度；和以第一声
波均方根幅度与第二均方根幅度的比值作为所述岩体弱面表征指数。

[0007] 在根据本发明的阵列声波测井识别和评价岩体弱面的方法的一个示例性实施例中，所述变异函数可以为：

[0008]

[0009] 式中：

[0010] Vsize为波形阵列数据的时间采样点数，不同测井系列，Vsize值不同；

[0011] ΔWAVE为N深度点经统计变异处理技术处理后的反映岩体弱面特征波形数据；

[0012] N为计算点的深度；

[0013] N+1为计算点的下一个深度；

[0014] M为当前计算点的波形数据阵列的时间采样点，不同测井系列，M值不同；

[0015] SMOTHR为深度域内进行比较的采样点数，所述采样点数为奇数，取值为1,3,5,7,9…；

[0016] WAVE为自动增益恢复后的阵列声波测井波形数据。

[0017] 在根据本发明的阵列声波测井识别和评价岩体弱面的方法的一个示例性实施例中，所述方法还可以包括：将经统计变异处理技术处理后的阵列声波测井波形数据与自动
增益恢复后的阵列声波测井波形数据相加，得到总体变异波形数据，以变密度的方式进行
灰度成图，绘制得到经总体统计变异处理后的第二波形变密度图，所述第二波形变密度图
上显示有岩体弱面层段与均质岩体层段之间的差异。

[0018] 在根据本发明的阵列声波测井识别和评价岩体弱面的方法的一个示例性实施例中，所述方法还可以包括：基于所述岩体弱面表征指数绘制岩体弱面表征指数曲线，给定分
类阈值或分类标准，对岩体弱面分布层段进行自动识别与岩体弱面分级评价。

[0019] 在根据本发明的阵列声波测井识别和评价岩体弱面的方法的一个示例性实施例中，所述方法还可以包括：绘制岩体弱面的解释评价成果图，所述解释评价成果图上包括第
一曲线道、第二曲线道、第三曲线道、第四曲线道、第五曲线道、第六曲线道、第七曲线道以
及岩体弱面分级评价道中的一个或多个。其中，第一曲线道上绘制有井径、无铀伽玛和自然
伽玛曲线；第二曲线道上绘制有纵波时差、补偿中子和补偿密度曲线；第三曲线道上绘制有
深侧向和浅侧向曲线；第四曲线道上绘制有经过所述增益恢复后的第三波形变密度图；第
五曲线道上绘制有所述第一波形变密度图；第六曲线道上绘制有所述第二波形变密度图；
第七曲线道上绘制有所述岩体弱面表征指数曲线；岩体弱面分级评价道上按不同图例绘制
了对应深度段岩体弱面的分级评价结果。

[0020] 在根据本发明的阵列声波测井识别和评价岩体弱面的方法的一个示例性实施例中，所述方法还可以包括：对阵列声波测井波形数据进行预处理，所述预处理包括数据格式
解编与数据加载、测井深度校正、阵列声波测井波形数据质量检查，其中，利用阵列声波波
形均方根振幅曲线进行阵列声波测井波形数据质量检查，甄别包括因设备故障等造成的错
误记录波形数据，选择期望的波形数据进行波形增益恢复和后续处理。

[0021] 本发明还提供了一种阵列声波测井识别和评价岩体弱面的系统，所述系统包括：增益恢复单元，被配置为利用自动增益恢复系数对采集到的阵列声波测井波形数据进行增
益恢复；统计变异处理单元，被配置为利用变异函数对增益恢复后的阵列声波测井波形数
据进行统计变异处理；灰度成图单元，被配置为对统计变异处理后的波形数据以变密度的
方式进行灰度成图，以显示岩体弱面的位置及规模；以及岩体弱面表征指数计算单元，被配
置为：对统计变异处理后的波形数据开时窗进行声波均方根幅度计算，得到第一均方根幅
度；对相同时窗标准硬地层统计变异处理后的波形数据开时窗进行声波均方根幅度计算，
得到第二均方根幅度；以第一均方根幅度与第二均方根幅度的比值作为所述岩体弱面表征
指数，其中，所述岩体弱面表征指数用于评价岩体弱面，所述岩体弱面表征指数越大，则岩
体弱面越容易出现滑移。

[0022] 在根据本发明的阵列声波测井识别和评价岩体弱面的系统的一个示例性实施例中，所述变异函数可以为：

[0023]

[0024] 式中：

[0025] Vsize为波形阵列数据的时间采样点数，不同测井系列，Vsize值不同；

[0026] ΔWAVE为N深度点经统计变异处理技术处理后的反映岩体弱面特征波形数据；

[0027] N为计算点的深度；

[0028] N+1为计算点的下一个深度；

[0029] M为当前计算点的波形数据阵列的时间采样点，不同测井系列，M值不同；

[0030] SMOTHR为深度域内进行比较的采样点数，所述采样点数为奇数，取值为1,3,5,7,9…；

[0031] WAVE为自动增益恢复后的阵列声波测井波形数据。

[0032] 在根据本发明的阵列声波测井识别和评价岩体弱面的系统的一个示例性实施例中，所述系统还可以包括：总体变异波形变密度生成单元，被配置为将经统计变异处理技术
处理后的阵列声波测井波形数据与自动增益恢复后的阵列声波测井波形数据相加，得到总
体变异波形数据，以变密度的方式进行灰度成图，绘制得到经总体统计变异处理后的第二
波形变密度图，所述第二波形变密度图上显示有岩体弱面层段与均质岩体层段之间的差
异。

[0033] 在根据本发明的阵列声波测井识别和评价岩体弱面的系统的一个示例性实施例中，所述系统还可以包括：岩体弱面表征指数分级评价单元，被配置为基于所述岩体弱面表
征指数绘制岩体弱面表征指数曲线，给定分类阈值或分类标准，对岩体弱面分布层段进行
自动识别与岩体弱面分级评价。

[0034] 在根据本发明的阵列声波测井识别和评价岩体弱面的系统的一个示例性实施例中，所述系统还可以包括：岩体弱面的解释评价成果图形成单元，用于绘制岩体弱面的解释
评价成果图，所述解释评价成果图上包括第一曲线道、第二曲线道、第三曲线道、第四曲线
道、第五曲线道、第六曲线道、第七曲线道以及岩体弱面分级评价道中的一个或多个。其中，
第一曲线道上绘制有井径、无铀伽玛和自然伽玛曲线；第二曲线道上绘制有纵波时差、补偿
中子和补偿密度曲线；第三曲线道上绘制有深侧向和浅侧向曲线；第四曲线道上绘制有经
过所述增益恢复后的第三波形变密度图；第五曲线道上绘制有所述第一波形变密度图；第
六曲线道上绘制有所述第二波形变密度图；第七曲线道上绘制有所述岩体弱面表征指数曲
线；岩体弱面分级评价道上按不同图例绘制了对应深度段岩体弱面的分级评价结果。

[0035] 在根据本发明的阵列声波测井识别和评价岩体弱面的系统的一个示例性实施例中，所述系统还可以包括：预处理模块，被配置为对阵列声波测井波形数据进行预处理，所
述预处理包括数据格式解编与数据加载、测井深度校正、阵列声波测井波形数据质量检查，
其中，利用阵列声波波形均方根振幅曲线进行阵列声波测井波形数据质量检查，甄别包含
因设备故障等造成的错误记录波形数据，选择期望的波形数据进行波形增益恢复和后续处
理。

[0036] 本发明还提供了一种设备，所述设备包括：处理器；存储器，存储有计算机程序，当所述计算机程序被处理器执行时，实现如上任一项所述的方法。

[0037] 本发明还提供了一种存储有计算机程序的计算机可读存储介质，当所述计算机程序在被处理器执行时实现如上任一项所述的方法。

[0038] 与现有技术相比，本发明的有益效果包括：可识别和评价井下深部地层岩体弱面；能突出早期被忽略的有关地层应力弱面结构（如：裂隙、节理、层理、软弱夹层、断层、断裂破
碎带等）的测井信息，从而达到地层应力弱面结构精细评价的目的；所需测井数据易得，原
理可靠，方法可行，设计合理，效果明显，且无需额外增加测井采集费用。

附图说明

[0039] 图1为根据本发明示例性实施例的阵列声波测井识别和评价岩体弱面的方法的流程图。

[0040] 图2为本发明的阵列声波测井波形质量控制与优选最佳波形进行统计变异处理的甄别示例。

[0041] 图3为本发明的阵列声波波形增益恢复前后对比示例。

[0042] 图4示出了根据本发明的示例性实施例的方法对阵列声波波形数据进行统计变异处理数据后的变密度成图。

[0043] 图5示出了根据本发明示例性实施例的方法对岩体弱面阵列声波波形数据统计变异处理后的解释评价成果图。

[0044] 图6示出了根据本发明示例性实施例的阵列声波测井识别和评价岩体弱面的系统的框图。

[0045] 附图标记说明：

[0046] 100‑增益恢复单元，200‑统计变异处理单元，300‑灰度成图单元，400‑岩体弱面表征指数计算单元，500‑总体变异波形变密度生成单元，600‑岩体弱面表征指数分级评价单
元，700‑解释评价成果图形成单元，800‑预处理模块。

具体实施方式

[0047] 在下文中，将结合附图和示例性实施例来详细说明本发明的阵列声波测井识别和评价岩体弱面的方法、系统、设备及介质。

[0048] 术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含
地包括一个或者更多个该特征。在本发明的描述中，除非另有说明，“多个”的含义是两个或
两个以上。

[0049] 随着页岩气等非常规油气勘探开发认识的提高，岩体弱面（包括：裂隙、节理、层理、软弱夹层、断层、断裂破碎带等）研究在页岩气地质甜点评价、页岩气水平井储层压裂改
造、以及套管变形研究与预防中的作用更加凸显，识别和评价岩体弱面的工作越来越受到
地质工程方面重视。

[0050] 目前，考虑到成本等问题，更多地利用阵列声波测井波形变密度、阵列声波处理时差、以及阵列声波幅度衰减等方式来识别和评价地层深部岩体弱面，但仍存在识别和评价
岩体弱面不直观、数据处理难度大、解释评价耗时长等限制，尤其对深部页岩气水平井岩体
弱面的快速直观识别和评价方面，问题更加突出。在本发明中，所提及“深部页岩气水平井
岩体弱面”中的“深部”是指非地面位置，例如，地表以下作业人员无法抵达的位置。

[0051] 由此，本发明提供了一种阵列声波测井快速直观识别和评价地层深部岩体弱面（亦称结构面）的方法及系统（装置）。

[0052] 本发明的主要技术构思是：对于均质岩体，在一定三维空间范围内，其声学特征应基本相同。具体在阵列声波测井上，应表现为随测量深度的变化，其波形的幅度、频率等应
相同。若出现差异，排除测量仪器故障原因造成的，应该是岩体弱面特征（如：裂隙、节理、层
理、软弱夹层、断层、断裂破碎带等）发生变化造成的。从统计学上来说，这种变化在数量标
志上则表现为变异。这种变异，则主要反映了岩体弱面的改变，指示岩体弱面的存在。

[0053] 为了识别和评价井下深部地层岩体弱面，为地质和工程应用提供岩体弱面（如：裂隙、节理、层理、软弱夹层、断层、断裂破碎带等）特征和信息，本发明提供了一种阵列声波测
井快速直观识别和评价岩体弱面的新方法，通过构建一个统计学意义上的变异函数，对阵
列声波测井波形数据进行快速处理后，利用变密度灰度成图方式，直观表征岩体弱面存在
的位置，进行岩体弱面分级评价。

[0054] 根据本发明的阵列声波测井识别和评价岩体弱面的方法可包括：利用自动增益恢复系数对采集到的阵列声波测井波形数据进行增益恢复；利用变异函数对增益恢复后的阵
列声波测井波形数据进行统计变异处理；对统计变异处理后的波形数据以变密度的方式进
行灰度成图，绘制得到经统计变异处理后的第一波形变密度图，第一波形变密度图上显示
有岩体弱面的位置及规模；以及计算岩体弱面表征指数，并基于岩体弱面表征指数来评价
岩体弱面。其中，所述岩体弱面表征指数越大，则岩体弱面越容易出现滑移，所述计算岩体
弱面表征指数具体包括：对统计变异处理后的波形数据开时窗进行声波均方根幅度计算，
得到第一均方根幅度；对相同时窗标准硬地层统计变异处理后的波形数据开时窗进行声波
均方根幅度计算，得到第二均方根幅度；以第一均方根幅度与第二均方根幅度的比值作为
所述岩体弱面表征指数。

[0055] 根据本发明提供了一种阵列声波测井快速直观识别和评价岩体弱面的新方法，构建一个统计学意义上的变异函数，对阵列声波测井波形数据进行快速处理后，通过变密度
灰度成图方式，可直观表征岩体弱面存在的位置，还可对岩体弱面进行分级评价。其对岩体
弱面进行分级评价的基础是：裂缝、断层、岩性、井眼扩径等在统计变异处理技术处理后的
波形成果图上的测井响应特征差异，原理可靠，方法可行，设计合理。

[0056] 根据本发明，对裸眼井阵列声波波形测井数据进行统计变异处理，能突出早期被忽略的有关地层应力弱面结构的测井信息，在不需要花费测井资料采集费用的情况下，从
而达到地层应力弱面结构精细评价的目的，它是岩体弱面评价所需费用最少的。统计变异
处理技术的应用，可以直接计算深度域内各深度间波形数据的起伏变化情况，根据处理后
的波形变密度图像，可以更清楚地、直观地识别裂缝、断层、岩性、井眼扩径等弱面的差别，
甚至确定裂缝、断层、岩性、井眼扩径等在不同井段的弱面微小变化。

[0057] 图1为根据本发明示例性实施例的阵列声波测井识别和评价岩体弱面方法的流程图。

[0058] 在第一示例性实施例中，根据本发明的阵列声波测井识别和评价岩体弱面的方法包括：测井数据预处理、波形数据增益恢复、波形数据统计变异处理、岩体弱面表征指数计
算、处理数据快速直观成图以及岩体弱面解释评价。

[0059] 图1中给出了本发明的阵列声波测井波形处理解释评价的步骤。如图1中所示，该方法具体包括：

[0060] 步骤S1：数据准备，包括收集整理常规测井数据、阵列声波数据、以及其它相关测井与地质资料。

[0061] 根据本发明，阵列声波测井波形数据可通过常用的阵列声波测井仪器采集得到。阵列声波测井仪器包括但不限于过钻头存储式阵列声波测井仪器TBDS、偶极横波成像测井
仪DSI、SMLS‑6210泵出式存储阵列声波测井仪器、WSTT交叉偶极子阵列声波测井仪器、MPAL
多极子阵列声波测井仪器、以及ECLIPSE5700成像测井系列所属DAC单极阵列声波测井仪
器、多极阵列声波测井仪MAC、XMAC交叉偶极子阵列声波测井仪器、XMAC‑Ⅱ交叉偶极子阵列
声波测井仪器、XMAC‑F1交叉偶极子阵列声波测井仪器等仪器所采集的波形数据。

[0062] 步骤S2：阵列声波预处理。

[0063] 阵列声波预处理包括：数据格式解编与数据加载、测井深度校正以及波形数据质量检查。

[0064] 除了对阵列声波测井数据进行数据格式解编、数据加载、深度校正等常规的预处理工作外，更重要地是对波形阵列声波测井波形数据质量检查，其目的主要是甄别包含因
设备故障等原因造成的错误记录波形数据，选择期望的（例如，低噪声、正确的）波形数据进
行波形增益恢复和后续处理解释评价。

[0065] 图2为本发明的阵列声波测井波形质量控制与优选最佳波形进行统计变异处理的甄别示例。数据来源于中国四川一口页岩气井TBDS阵列声波测井实施例，从图2中的斯通利
波均方根幅度曲线（箭头符号所指曲线）以及3#接收器的波形外观可以看出，阵列声波测井
斯通利波3#接收器的波形异常，在选择波形进行统计变异处理时，3#接收器的波形不宜作
为波形统计变异处理的目标。

[0066] S3、对阵列声波测井波形数据进行增益恢复。

[0067] 该步骤可包括：确定仪器类型和确定增益恢复系数。

[0068] 阵列声波测井中，在没有超出测井数据采集系统的动态接收范围的情况下，总希望阵列声波测井采集到的数据幅度尽可能地大。然而，实际测量的时候，从一个深度点到下
一个深度点，因岩体岩块成分和弱面（如：裂隙、节理、层理、软弱夹层、断层、断裂破碎带等）
的变化，会导致阵列声波测井采集到的波形幅度变化可能很大，从而超出测井数据采集系
统的动态接收范围，若采用固定的增益控制办法，采集到的波形数据就会出现饱和，则无法
取得反映岩体岩块成分和弱面真实变化的测井采集效果。为此，阵列声波测井仪采用自动
增益控制办法，通过计算出一个合适的增益值来优化波形。当在进行波形数据处理时，为了
得到真实的波形，必须去除增益控制，达到还原原始测量波形的目的。

[0069] 对阵列声波测井波形数据进行自动增益恢复计算公式如下式1所示：

[0070] 式1：

[0071] WAVE=a*WAVE0

[0072] 在式1中：WAVE0为阵列声波测井波形数据；WAVE为自动增益恢复后的阵列声波测井波形数据，a为自动增益恢复系数。

[0073] 自动增益恢复系数a可参照对应仪器生产厂家提供的计算公式。

[0074] 例如，对于ECLIPSE5700成像测井系列阵列声波仪，其自动增益恢复系数a的计算式如下式2。

[0075] 式2：

[0076] a=1/（10*（3*Trgain/20））

[0077] 在式2中：Trgain为ECLIPSE5700成像测井系列阵列声波自动增益数据。

[0078] 又例如，对于ThruBit存储式测井系列TBDS阵列声波仪，其自动增益恢复系数a的计算式见式3。

[0079] 式3：

[0080] a=1/Trgain。

[0081] 在式3中：Trgain为ThruBit存储式测井系列阵列声波仪自动增益数据。

[0082] 图3为本发明的阵列声波波形增益恢复前后对比示例。数据来源于中国四川一口页岩气井TBDS阵列声波测井实施例，图3中绘制的波形采用刻度完全一致。从图3中可以看
出，阵列声波波形增益恢复前的WAVE0与阵列声波波形增益恢复后的WAVE差别迥异。因此，
在进行波形统计变异处理前，需要对采集的波形数据进行增益恢复，以保证处理结果的正
确性。

[0083] 步骤S4：波形数据统计变异处理。

[0084] 步骤S41、波形统计变异处理。

[0085] 具体地，可构建一个统计学意义上的变异函数，并对经增益恢复后的阵列声波测井波形数据进行统计变异处理。波形数据统计变异处理函数见下式4。

[0086] 式4：

[0087]

[0088] 在式4中：Vsize为波形阵列数据的时间采样点数，不同测井系列，Vsize值不同；ΔWAVE为N深度点经统计变异处理技术处理后的反映岩体弱面特征波形数据；N为计算点的深
度；N+1为计算点的下一个深度；M为当前计算点的波形数据阵列的时间采样点，不同测井系
列，M值不同；SMOTHR为深度域内进行比较的采样点数，一般为奇数，取值为1,3,5,7,9…。
WAVE为自动增益恢复后的阵列声波测井波形数据。

[0089] 步骤S42、波形总体变异处理。

[0090] 具体地，为了对比岩体弱面层段与均质岩体层段之间的差异，该方法还可包括计算了另外一条总体变异波形曲线TWAVE，具体地，将经统计变异处理技术处理后的阵列声波
测井波形数据与自动增益恢复后的阵列声波测井波形数据相加，得到总体变异波形曲线
TWAVE，其计算公式见式5。

[0091] 式5：

[0092] TWAVE=WAVE+ΔWAVE。

[0093] 在式5中：TWAVE表示总体变异波形数据，WAVE为自动增益恢复后的声波波形数据，ΔWAVE为经统计变异处理技术处理后的波形数据。

[0094] 图4示出了根据本发明的示例性实施例的方法对阵列声波波形数据进行统计变异处理数据后的变密度成图。数据来源于中国四川一口页岩气井XMAC‑F1阵列声波测井实施
例。在图4中，从第三曲线道采用灰度变密度方式显示的经过统计变异处理后的波形变密度
图上，可以清楚地看出“V”字形干涉条纹（图中箭头所指位置），而在第二曲线道采用灰度变
密度方式绘制了经过增益恢复后的波形变密度图上基本无法看到。第四曲线道上所显示的
是采用灰度变密度方式绘制的经过总体统计变异处理后的波形变密度图，其既能看出波形
的总体变化，又能看到“V”字形干涉条纹，还可以展现出岩体弱面层段与均质岩体层段之间
的差异。

[0095] 步骤S5：岩体弱面表征指数计算。

[0096] 岩体弱面表征指数计算就是对统计变异处理技术处理后的波形数据开时窗进行声波均方根幅度计算（也称为均方根能量），将其除以相同时窗标准硬地层的统计变异处理
技术处理后波形数据声波均方根幅度计算值（也称为均方根能量），用它们的比值来表征岩
体弱面的发育情况，即：岩体弱面表征指数。根据计算的岩体弱面表征指数曲线，给定一个
适宜的分类阈值或分类标准，可对岩体弱面分布层段自动识别与岩体弱面分级评价。

[0097] 具体地，包括以下步骤：

[0098] S51、对于经过统计变异处理后的第一波形数据，将其目标层段内的波形数据开时窗进行均方根幅度计算，得到第一均方根幅度。这里，目标层段包括标准硬地层。这里，目的
层段是指利用阵列声波测井资料开展岩体弱面识别与评价的测量井段。

[0099] S52、对于经过统计变异处理后的第一波形数据，将与步骤S51相同时窗内的标准硬地层的波形数据进行均方根幅度计算，得到第二均方根幅度。这里，标准硬地层是指目标
层段中声波纵波时差最小、横向分布稳定的区域标志层。

[0100] S53、以第一均方根幅度与第二均方根幅度之比作为岩体弱面表征指数。

[0101] 岩体弱面表征指数计算公式见下式6。

[0102] 式6：

[0103]

[0104] 在式6中：ZFLEX为岩体弱面表征指数，单位：无量纲；

[0105] Amp为统计变异处理技术处理后的波形数据开时窗声波均方根幅度，即，第一均方根幅度；

[0106] Amp0为统计变异处理技术处理后的波形数据开时窗标准硬地层的声波均方根幅度，即，第二均方根幅度。

[0107] 统计变异处理技术处理后的波形数据开时窗声波均方根幅度计算见下式7。

[0108] 式7：

[0109]

[0110] 在式7中：ΔWavei：经统计变异处理技术处理后的波形第i点的数据；N为时窗窗口内波形数据的点数。

[0111] 相同时窗标准硬地层的统计变异处理技术处理后波形数据声波均方根幅度计算见下式8。

[0112] 式8：

[0113]

[0114] 在式8中，ΔWaveBi表示相同时窗标准硬地层的经统计变异处理技术处理后波形数据波形第i点的数据；N为时窗窗口内波形数据的点数。

[0115] S54、基于岩体弱面表征指数绘制岩体弱面表征指数曲线。

[0116] 图5示出了根据本发明示例性的方法对岩体弱面阵列声波波形数据统计变异处理后的解释评价成果图。利用本发明开展岩体弱面表征指数计算示例。数据来源于中国四川
一口页岩气井XMAC‑F1阵列声波测井示例。

[0117] 在图5中，第五曲线道采用灰度变密度方式绘制了经过统计变异处理后的波形变密度图，从该曲线道采用灰度变密度方式显示的经过统计变异处理后的波形变密度图上，
可以清楚地看出“V”字形干涉条纹。该曲线道中的TTWINST曲线与TTWINEND曲线包络的统计
变异处理后的波形变密度区域，即为开时窗计算统计变异处理后波形均方根幅度区域，用
于计算岩体弱面表征指数曲线ZFLEX。这里，TTWINST曲线：利用第一波形数据开时窗计算第
一均方根幅度的起始时间，单位为：微秒。TTWINEND曲线：利用第一波形数据开时窗计算第
一均方根幅度的结束时间，单位为：微秒。

[0118] 在图5中，第七曲线道绘制了岩体弱面表征指数曲线ZFLEX及分级评价，岩体弱面表征指数越大，则岩体弱面越容易出现滑移。

[0119] 步骤S6：波形统计变异变密度成图。

[0120] 处理数据快速直观成图，包括统计变异波形变密度成图和总体变异波形变密度成图。具体地，可以对增益恢复后的波形数据、对统计变异处理后的波形数据和总体变异波形
数据以变密度的方式进行灰度成图，直观显示岩体弱面的位置及规模。

[0121] 图4示出了根据本发明的示例性实施例的方法对阵列声波波形数据进行统计变异处理数据后的变密度成图。数据来源于中国四川一口页岩气井XMAC‑F1阵列声波测井实施
例。在图4中，第一曲线道绘制了井径、无铀伽玛、自然伽玛曲线，第二曲线道采用灰度变密
度方式绘制了经过增益恢复后的波形变密度图，第三曲线道采用灰度变密度方式绘制了经
过统计变异处理后的波形变密度图，第四曲线道采用灰度变密度方式绘制了经过总体统计
变异处理后的波形变密度图。从第三曲线道采用灰度变密度方式显示的经过统计变异处理
后的波形变密度图上，可以清楚地看出“V”字形干涉条纹（图中箭头所指位置）。在第四曲线
道采用的灰度变密度方式绘制的经过总体统计变异处理后的波形变密度图上，既能看出波
形的总体变化，又能看到变化明显的“V”字形干涉条纹。

[0122] 步骤S7：岩体弱面解释评价。

[0123] 包括：弱面分级评价、综合评价与成图。

[0124] 岩体弱面解释评价是指结合统计变异处理技术处理后的波形数据的变密度灰度成图，根据计算的岩体弱面表征指数曲线，给定一个适宜的分类阈值或分类标准，对岩体弱
面分布层段进行自动识别与岩体弱面分级评价。

[0125] 图5示出了根据本发明示例性的方法对岩体弱面阵列声波波形数据统计变异处理后的解释评价成果图。数据来源于中国四川一口页岩气井XMAC‑F1阵列声波测井实施例。

[0126] 在图5中，第一曲线道绘制了井径、无铀伽玛、自然伽玛曲线，第二曲线道绘制了纵波时差、补偿中子、补偿密度曲线，第三曲线道绘制了深浅双侧向曲线，第四曲线道采用灰
度变密度方式绘制了经过增益恢复后的波形变密度图。

[0127] 在图5中，第五曲线道采用灰度变密度方式绘制了经过统计变异处理后的波形变密度图，从该曲线道采用灰度变密度方式显示的经过统计变异处理后的波形变密度图上，
可以清楚地看出“V”字形干涉条纹。

[0128] 在图5中，第六曲线道采用灰度变密度方式绘制了经过总体统计变异处理后的波形变密度图，从该曲线道采用的灰度变密度方式绘制的经过总体统计变异处理后的波形变
密度图上，既能看出波形的总体变化，又能看到变化明显的“V”字形干涉条纹。

[0129] 在图5中，第七曲线道绘制了岩体弱面表征指数曲线及分级评价，岩体弱面表征指数越大，则岩体弱面越容易出现滑移。

[0130] 在岩体弱面分级评价结论道可以按不同图例（例如灰度填充块、点填充块等）绘制对应深度段岩体弱面分级评价结论。本实施例中，对岩体弱面表征指数按0 0.2，0.2 0.4，
~ ~
0.4 0.6，0.6 0.8，0.8 1对应分为五个等级（Ⅰ级、Ⅱ级、Ⅲ级、Ⅳ级、Ⅴ级），以对岩体弱面进
~ ~ ~
行分级评价。但本发明不限于此，对岩体弱面表征指数进行分级的级数以及范围可根据需
调整。

[0131] 表1给出了某地区岩块弱面阵列声波波形统计变异处理解释评价成果表，其给出了利用本发明的方法进行的岩体弱面解释评价成果，包括不同级别弱面包围的岩体井段、
岩体厚度、以及弱面级别（级别数字越大，岩体弱面越容易出现滑移）。

[0132] 表1

[0133] 序号深度（m） 厚度（m） 岩体弱面评价结论　序号深度（m） 厚度（m）岩体弱面评价结论1 2675.5~2683.2 7.7 Ⅲ级 　29 2850.8~2852.2 1.4 Ⅲ级
2 2683.2~2683.5 0.3 Ⅴ级 　30 2852.2~2860.1 7.9 Ⅰ级
3 2683.5~2685 1.5 Ⅲ级 　31 2860.1~2861.3 1.2 Ⅲ级
4 2685~2685.9 0.9 Ⅱ级 　32 2861.3~2882.3 21 Ⅰ级
5 2685.9~2689.1 3.2 Ⅲ级 　33 2882.3~2883.7 1.4 Ⅲ级
6 2689.1~2693.1 4 Ⅰ级 　34 2883.7~2888.3 4.6 Ⅰ级
7 2693.1~2698.1 5 Ⅲ级 　35 2888.3~2890.9 2.6 Ⅲ级
8 2698.1~2701.4 3.3 Ⅰ级 　36 2890.9~2895.1 4.2 Ⅰ级
9 2701.4~2704.9 3.5 Ⅲ级 　37 2895.1~2895.4 0.3 Ⅲ级
10 2704.9~2706.8 1.9 Ⅰ级 　38 2895.4~2932.9 37.5 Ⅰ级
11 2706.8~2710.3 3.5 Ⅲ级 　39 2932.9~2933.9 1 Ⅲ级
12 2710.3~2712.7 2.4 Ⅰ级 　40 2933.9~2935.2 1.3 Ⅰ级
13 2712.7~2713.6 0.9 Ⅲ级 　41 2935.2~2935.8 0.6 Ⅴ级
14 2713.6~2714.2 0.6 Ⅱ级 　42 2935.8~2937.2 1.4 Ⅲ级
15 2714.2~2715 0.8 Ⅳ级 　43 2937.2~2939 1.8 Ⅰ级
16 2715~2715.3 0.3 Ⅰ级 　44 2939~2939.3 0.3 Ⅲ级
17 2715.3~2715.8 0.5 Ⅱ级 　45 2939.3~2939.9 0.6 Ⅴ级
18 2715.8~2777 61.2 Ⅰ级 　46 2939.9~2940.6 0.7 Ⅳ级
19 2777~2777.8 0.8 Ⅲ级 　47 2940.6~2940.9 0.3 Ⅰ级
20 2777.8~2780.8 3 Ⅰ级 　48 2940.9~2945.7 4.8 Ⅲ级
21 2780.8~2784 3.2 Ⅲ级 　49 2945.7~2947.9 2.2 Ⅰ级
22 2784~2787.4 3.4 Ⅰ级 　50 2947.9~2949.9 2 Ⅲ级
23 2787.4~2788.5 1.1 Ⅲ级 　51 2949.9~2950.9 1 Ⅱ级
24 2788.5~2801.6 13.1 Ⅰ级 　52 2950.9~2963.1 12.2 Ⅲ级
25 2801.6~2801.9 0.3 Ⅲ级 　53 2963.1~2994.8 31.7 Ⅰ级
26 2801.9~2807.2 5.3 Ⅰ级 　54 2994.8~2995.3 0.5 Ⅱ级
27 2807.2~2808.1 0.9 Ⅲ级 　55 2995.3~2996 0.7 Ⅳ级
28 2808.1~2850.8 42.7 Ⅰ级 　56 2996~2999.1 3.1 Ⅰ级

[0134] 以上仅示出根据本发明构思的阵列声波测井岩体弱面识别和评价的方法的示例性实施例，然而本发明构思不限于此，本发明的方法也可不按上述步骤顺序地执行，或者可
根据工程需要省略或添加一个或多个步骤及其子步骤。例如，在一个实施例中，可省略步骤
S42。

[0135] 图6示出了根据本发明示例性实施例的阵列声波测井识别和评价岩体弱面的系统的框图。

[0136] 在第二示例性实施例中，提供了一种阵列声波测井识别和评价岩体弱面的系统，如图6中所示，该系统可以包括：增益恢复单元100、统计变异处理单元200、灰度成图单元
300、岩体弱面表征指数计算单元400。

[0137] 其中，增益恢复单元100被配置为利用自动增益恢复系数对采集到的阵列声波测井波形数据进行增益恢复。对阵列声波测井波形数据进行自动增益恢复计算公式可以为前
文式2所示，在此不再赘述。自动增益恢复系数a可参照对应仪器生产厂家提供的计算公式。

[0138] 统计变异处理单元200被配置为利用变异函数对增益恢复后的阵列声波测井波形数据进行统计变异处理。所述变异函数可以为前文式4，在此不再赘述。

[0139] 灰度成图单元300被配置为对统计变异处理后的波形数据以变密度的方式进行灰度成图，以显示岩体弱面的位置及规模。

[0140] 岩体弱面表征指数计算单元400被配置为：对统计变异处理后的波形数据开时窗进行声波均方根幅度计算，得到第一均方根幅度；对相同时窗标准硬地层统计变异处理后
的波形数据开时窗进行声波均方根幅度计算，得到第二均方根幅度；以第一均方根幅度与
第二均方根幅度的比值作为所述岩体弱面表征指数。岩体弱面表征指数如前文式6中所示，
在此不再赘述。其中，所述岩体弱面表征指数用于评价岩体弱面，所述岩体弱面表征指数越
大，则岩体弱面越容易出现滑移。

[0141] 可选择地，所述系统还可包括总体变异波形变密度生成单元500。总体变异波形变密度生成单元500可以被配置为将经统计变异处理技术处理后的阵列声波测井波形数据与
自动增益恢复后的阵列声波测井波形数据相加，得到总体变异波形数据，所涉及计算式见
前文式5，在此不再赘述。将总体变异波形数据以变密度的方式进行灰度成图，绘制得到经
总体统计变异处理后的第二波形变密度图，所述第二波形变密度图上显示有岩体弱面的位
置及规模。

[0142] 可选择地，所述系统还可包括：岩体弱面表征指数分级评价单元600，其被配置为基于所述岩体弱面表征指数绘制岩体弱面表征指数曲线，给定分类阈值或分类标准，对岩
体弱面分布层段进行自动识别与岩体弱面分级评价。

[0143] 可选择地，所述系统还可包括：岩体弱面的解释评价成果图形成单元700，用于绘制岩体弱面的解释评价成果图，所述解释评价成果图上包括第一曲线道、第二曲线道、第三
曲线道、第四曲线道、第五曲线道、第六曲线道、第七曲线道以及岩体弱面分级评价道中的
一个或多个。例如，图5中示出了八个道。

[0144] 在图5中，第一曲线道上绘制有岩性曲线，即，含井径、无铀伽玛和自然伽玛的曲线；第二曲线道上绘制有孔隙度曲线，即，含纵波时差、补偿中子和补偿密度的曲线；第三曲
线道上绘制有电阻率曲线，即，含深侧向和浅侧向的双侧向的曲线；第四曲线道上绘制有经
过所述增益恢复后的第三波形变密度图；第五曲线道上绘制有所述第一波形变密度图；第
六曲线道上绘制有所述第二波形变密度图；第七曲线道上绘制有所述岩体弱面表征指数曲
线；岩体弱面分级评价道上按不同图例绘制了对应深度段岩体弱面的分级评价结果。

[0145] 可选择地，所述系统还可包括预处理模块800。预处理模块800被配置为对阵列声波测井波形数据进行预处理，所述预处理可包括数据格式解编与数据加载、测井深度校正、
阵列声波测井波形数据质量检查。

[0146] 可选择地，可利用阵列声波波形均方根振幅曲线进行阵列声波测井波形数据质量检查，甄别设备故障等造成的错误记录波形数据，选择低噪声、正确的波形数据进行波形增
益恢复和后续处理。

[0147] 上述增益恢复单元100、统计变异处理单元200、灰度成图单元300、岩体弱面表征指数计算单元400、总体变异波形变密度生成单元500、岩体弱面表征指数分级评价单元
600，岩体弱面的解释评价成果图形成单元700、预处理模块800执行操作的具体方式在已经
在有关该方法的实施例中（例如上述步骤S1至步骤S7中）进行了详细描述，因此在此不再对
各个单元的具体内容进行详细阐述说明。

[0148] 根据本发明的方法可以被编程为计算机程序并且相应的程序代码或指令可以被存储在计算机可读存储介质中，当程序代码或指令被处理器执行时使得处理器执行上述阵
列声波测井识别和评价岩体弱面方法，上述处理器和存储器可以被包括在计算机设备中。

[0149] 在第三示例性实施例中，还提供一种存储有计算机程序的计算机可读存储介质。该计算机可读存储介质存储有当被处理器执行时使得处理器执行根据本发明的阵列声波
测井识别和评价岩体弱面方法的计算机程序。该计算机可读记录介质是可存储由计算机系
统读出的数据的任意数据存储装置。计算机可读记录介质的示例包括：只读存储器、随机存
取存储器、只读光盘、磁带、软盘、光数据存储装置和载波（诸如经有线或无线传输路径通过
互联网的数据传输）。

[0150] 在第四示例性实施例中，还提供一种计算机设备。该计算机设备包括处理器和存储器。存储器用于存储计算机程序。计算机程序被处理器执行使得处理器执行根据本发明
的阵列声波测井识别和评价岩体弱面方法的计算机程序。

[0151] 根据本发明，首先进行测井数据预处理，接着对波形数据进行增益恢复，再对波形增益恢复数据进行统计变异处理和岩体弱面表征指数计算，最后对处理数据快速直观成图
和岩体弱面解释评价。对裸眼井阵列声波波形测井数据进行统计变异处理，能突出早期被
忽略的有关地层应力弱面结构的测井信息，在不需要花费测井资料采集费用的情况下，从
而达到地层应力弱面结构精细评价的目的，它是岩体弱面评价所需费用最少的。

[0152] 根据本发明的方法，统计变异处理技术的应用，可以直接计算深度域内各深度间波形数据的起伏变化情况，根据处理后的波形变密度图像，可以更清楚地、直观地识别裂
缝、断层、岩性、井眼扩径等弱面的差别，甚至确定裂缝、断层、岩性、井眼扩径等在不同井段
的弱面微小变化。所需测井数据易得，效果明显，且无需额外增加测井采集费用。

[0153] 尽管上面已经通过结合示例性实施例描述了本发明，但是本领域技术人员应该清楚，在不脱离权利要求所限定的精神和范围的情况下，可对本发明的示例性实施例进行各
种修改和改变。