三分量地震数据处理以及水力压裂地震解释方法

摘要

系统和方法包括用于呈现同步地震数据和压裂处理时间的解释结果的计算机实施方法。生成从三分量传感器获得的传感器读数的标准格式的地震数据集。将与这些传感器读数相对应的坐标和记录时间添加到该标准格式的地震数据集中。通过将地震记录时间与压裂处理时间进行同步，从该标准格式的地震数据集中生成同步地震数据。通过从该同步地震数据中去除死道和异常数据样本来生成质量受控的同步地震数据。对该质量受控的同步地震数据执行时频分析，包括执行短时傅里叶变换以分析傅里叶谱随时间的变化。基于该时频分析，在不同的时间样本处从每个频谱中提取共振频率。将基于的解释结果呈现给用户。

说明书

优先权声明

本申请要求2020年4月6日提交的美国专利申请号16/840,903的优先权，该专利申请的全部内容通过引用并入本文。

背景技术

本公开适用于监测压裂(fracturing)(或水力压裂(fracking))过程的技术。

水力压裂地震(Seismic While Fracking，SWF)是一种用于评估水力压裂期间水力裂缝的质量的地震技术。常规的SWF技术可以包括微震监测方法，其可能需要许多传感器。常规技术还可能需要大量的处理时间和计算资源，并且可能是劳动密集型的。

发明内容

本公开描述了可以用于基于同步地震数据和压裂处理时间生成解释结果的技术。

在一些实施方式中，一种计算机实施的方法包括以下内容。生成从三分量传感器获得的传感器读数的标准格式的地震数据集。将与这些传感器读数相对应的坐标和记录时间添加到该标准格式的地震数据集中。通过将地震记录时间与压裂处理时间进行同步，从该标准格式的地震数据集中生成同步地震数据。通过从该同步地震数据中去除死道和异常数据样本来生成质量受控的同步地震数据。对该质量受控的同步地震数据执行时频分析，包括执行短时傅里叶变换以分析傅里叶谱随时间的变化。基于该时频分析，在不同的时间样本处从每个频谱中提取共振频率。将基于的解释结果呈现给用户。

前面描述的实施方式可以使用以下各项来实施：计算机实施的方法；存储用于执行计算机实施的方法的计算机可读指令的非暂态计算机可读介质；以及计算机实施的系统，该计算机实施的系统包括与硬件处理器可互操作地耦接的计算机存储器，该硬件处理器被配置为执行计算机实施的方法/存储在非暂态计算机可读介质上的指令。

本说明书所描述的主题可以以特定实施方式实施，以实现以下优点中的一个或多个。首先，本公开中描述的技术可以提供相对于使用从裂缝散射的体波的常规技术的改进。体波往往很弱，因此可能需要宽方位角和密集的接收器阵列来重建波场。相比之下，本公开中描述的技术可以替代地使用管波，管波在井筒内以更少的能量耗散传播，因此需要更少的接收器。其次，可以确定一个或几个三分量传感器来有效地记录反射的管波。第三，本公开的技术可以减少计算、人力和运输成本。第四，本公开的技术可以使用反射的管波来提供对压裂质量的高效且可靠的评估。

本说明书的主题的一个或多个实施方式的细节在具体实施方式、附图和权利要求中进行阐述。从具体实施方式、权利要求以及附图，本主题的其他特征、方面和优点将变得显而易见。

附图说明

图1是根据本公开的一些实施方式的用于压裂地震测试(seismic whilefracturing test)的井轨迹的几何图。

图2A和图2B是示出根据本公开的一些实施方式的传感器记录的原始数据的图。

图3A和图3B是示出根据本公开的一些实施方式的对地震数据进行的时频分析的绘图。

图4A至图4C是示出了根据本公开的一些实施方式的在第15压裂阶段采集的地震数据的解释的图。

图5是根据本公开的一些实施方式的用于基于同步地震数据和压裂处理时间生成解释结果的方法的示例的流程图。

图6是展示了根据本公开的一些实施方式的用于提供与如在本公开中描述的所描述的算法、方法、功能、过程、流程和规程相关联的计算功能的示例计算机系统的框图。

不同附图中类似的附图标记和名称表示类似的元件。

具体实施方式

以下具体实施方式描述了用于基于同步地震数据和压裂处理时间来生成解释结果的技术。可以对所公开的实施方式进行各种修改、变更和置换，并且这些修改、变更和置换对于本领域的普通技术人员来说将是显而易见的，并且所定义的一般原理可以应用于其他实施方式和应用，而不脱离本公开的范围。在一些情况下，获得对所描述主题的理解不必要的细节可以省略，以免用不必要的细节混淆一个或多个所描述的实施方式，并且因为这些细节在本领域普通技术人员的技能范围内。本公开不旨在限于所描述或展示的实施方式，而是要被赋予与所描述的原理和特征一致的最广泛的范围。

一些常规的水力阻抗测试或类似技术可能会使用水锤信号。然而，这些技术可能仅限于以不超过1赫兹(Hz)的频率采样的井口压力时间序列。在井眼声学背景下，其他技术可以用于分析来自裂缝的波反射。由井眼测井工具产生的高频信号(例如高于1千赫兹(kHz)的信号)可以用于对裂缝和多孔层进行成像。相比之下，本公开中描述的技术可以使用小于大约10Hz的低频信号。用于推断水力裂缝导流能力的常规技术可以依赖于在地面上安装单分量传感器。传感器可以安装在井口内的流体附近、流体处或与流体接触以采集管波。在一些实施方式中，可以安装与井口外部接触的三分量传感器。由于管波是沿井筒传播的导波，因此波动包括x分量、y分量和z分量。这些分量可以被一个或多个三分量传感器完全捕获。因此，与使用单分量传感器的常规采集技术相比，本公开中描述的技术可以利用更全面的数据来更准确和稳健地推断地下裂缝信息。

在一些实施方式中，多种技术可以用于使用在压裂处理期间从水力裂缝反射的管波来推断裂缝连通性。例如，这些技术只需要处理和解释由与井口外部接触的三分量传感器记录的几个地震数据通道。可以示出解释的结果以提供与压裂处理曲线的良好匹配。这些技术可以用于水力压裂的实时监测。术语“实时”可以对应于在指定时间段内(例如，在几小时、几分钟或几秒内)发生的事件。

水力压裂地震(SWF)系统的主要目标是在水力压裂过程期间提供实时水力裂缝信息以帮助做出决策。在常规的水力压裂作业期间，工程师可能主要依赖先前的水力压裂信息、正在进行的水力压裂测量值、以及在不同水力压裂阶段的间隔期间获得的测井记录。这种类型的信息不能实时获得，这种情况会使工程师无法做出及时的决策。典型的微震监测方法可能需要许多传感器，这会大量增加处理和解释数据所需的成本、时间和劳动力。本公开提出了用于处理和解释例如通过与井口外部接触的单个三分量传感器收集的地震数据的具有成本效益的技术。这些技术的输出可以包括水力裂缝的实时导流能力参数。

数据采集

图1是根据本公开的一些实施方式的用于压裂地震测试(seismic whilefracturing test)的井轨迹的几何图100。为了获得几何形状100，可以将单个三分量传感器放置在地面上，与井口102的外部接触。传感器可以与数据存储设备和电池连接。由于传感器只有三个通道，因此即使数据记录可能持续数天，数据量也可能很小。以这种方式，数据处理和解释可以作为水力裂缝实时监测的一部分。图1中描绘的几何形状100是相对于例如由东方向104、北方向106和深度方向108给出的x分量、y分量和z分量。几何形状100可以包括具有几何中点112的26个部分110。在一个示例中，几何中点112可以在这26个部分110的第13段与第14段之间。

三分量传感器可以设计为对宽带敏感。入射管波可以通过井筒中的压力脉冲、微震事件或重大的快速流体流速变化(比如水锤)来激发。例如，水锤可能由大量流体的快速注入或释放、或者突然流动引起。所产生的受激管波可以作为导波(例如，斯通莱波)在井筒中传播。当这些波遇到井筒的边界和裂缝时，这些导波可以转换成克劳克利斯波。在裂缝内，克劳克利斯波可以沿着主裂缝连接处来回反射，从而形成低衰减的驻波。克劳克利斯波可以在井筒内转换回导波。最后，导波可以反射回地面并由地面传感器记录。

图2A和图2B是示出根据本公开的一些实施方式的传感器记录的原始数据200的图。图2A示出了原始数据的三分量时间序列，包括分别来自x分量200a、y分量200b和z分量200c的道。分量200a-200c示出了图2A中相对于幅度202和时间204绘制的信号，其中，信号对应于井筒的连续部分206(例如，对应于部分110)。

图2B示出了在三个通道处分别针对x分量、y分量和z分量200a-200c的传感器记录的对应傅里叶谱。这些谱包括以150Hz的频率208排齐的等距谱的特征。叠加压裂阶段的时间窗之后，数据可以在不同阶段示出一致的模式。x分量、y分量和z分量200a-200c的绘图相对于幅度202、时间204(小时)、频率208(Hz)和幅度210绘制。

数据处理工作流程

为了从反射的管波估计裂缝导流能力，可以对原始数据进行分析和处理。在一些实施方式中，用于估计裂缝导流能力的数据处理工作流程可以包括以下内容。首先，在现场数据解码步骤中，可以从传感器下载原始数据并将其转换为标准格式，比如勘探地球物理学家协会(SEG)Y(SEGY)格式。其次，在SEGY道头编辑步骤中，可以将包括坐标和记录时间的SEGY道头添加到SEGY文件中。第三，在时间同步步骤中，地震记录时间可以与压裂处理时间同步。第四，在数据质量控制步骤中，可以去除死道和异常数据样本。第五，在时频分析步骤中，可以使用短时傅里叶变换来分析傅里叶谱随时间的变化。第六，可以进行共振频率提取步骤。例如，在每个时间样本处，可以通过堆叠不同频率间隔的频谱来从频谱中提取共振频率。共振频率是产生最强堆叠功率的频率区间。

图3A和图3B是示出根据本公开的一些实施方式的对地震数据的时频分析的绘图300a和300b。例如，针对在第13压裂阶段采集的地震数据示出了数据处理结果。图3A示出了相对于幅度轴302和时间轴304绘制的三分量数据，包括在第13水力压裂阶段的x分量、y分量和z分量的数据。例如，可以通过对三分量数据应用短时傅里叶变换并堆叠它们的傅里叶谱来获得绘图300a。绘图300包括x分量306、y分量308和z分量310的幅度。图3B示出了相对于频率轴312和时间轴304绘制的最终时频谱。信息可以以高速记录，但分析可以集中在小于10Hz的信号上。

数据解释

图4A至图4C是示出了根据本公开的一些实施方式的在第15压裂阶段采集的地震数据的解释的图。图4A示出了时频谱400，包括示出数据解释结果。相对于时间轴402绘制的共振频率404(例如，以Hz为单位)的行为可以与压力曲线一致。由于新生成的裂缝释放了岩层中的应力，因此可以通过压力曲线来指示裂缝的产生，其中随着浆液的稳定注入，压力迅速下降。从时频谱可以确定，新的充满流体的裂纹可以产生新的共振频率。其原因是，共振源自于克劳克利斯波在由射孔点和裂缝尖端困住的空间内的来回反射。

图4B示出了相对于时间轴408和裂缝尺度410的解释结果406。例如，可以从时间谱中提取共振频率，并且可以计算裂缝导流能力的直方图，如图4B所示。

图4C示出了压裂处理曲线412。包括在压裂处理曲线412中的直方图的趋势与压力曲线一致。为了详细阐明该解释结果，第15压裂阶段可分为三个时期。

在射孔阶段414期间(从时间A 416到时间B 418)，可以将具有低砂浓度的浆液泵送到连续油管中。砂可以通过喷嘴注入到目标储层。在该时期，地层中以天然裂缝为主，裂缝导流能力较低。

在压裂时期420期间(从时间B 418到时间C 422)，压裂处理可以在时间B 418开始，并且压力可以随着浆液注入而增加。当产生一些裂缝时，压力开始下降。为了产生更多裂缝，通常需要多轮注浆。在该示例中，有五次注入。解释结果证实裂缝导流能力随着注入规程逐渐增加。

在压裂时期424(从时间C 422到时间D 426)期间，排水量和连接管压力降低，而砂浓度量增加。解释结果表明，尽管逐渐注入携砂液，但随着压裂压力的降低，地层中的水力裂缝逐渐减小。最后，随着越来越多的支撑剂注入地层，裂缝导流能力趋于稳定。

压裂处理曲线412包括连接管压力428、环空压力430、浆液速率432、砂浓度434和所添加液体436的线图。例如，压力(比如连接管压力428和环空压力430)可以以兆帕(MPa)为单位测量。例如，浆液速率(比如浆液速率432)可以以立方米每分钟(m

解释结果可以表明，从共振频率提取的裂缝导流能力与压裂处理曲线具有高度相关性，这可以实时提供有关压裂质量的宝贵信息。此外，13部分的压裂增压时间可以在很短的时间内进行。在整个压裂处理期间，从共振频率计算得到的导流因子可以示出稳定的裂缝分布。还可以比较不同阶段的导流能力直方图，以评估整个水力裂缝的处理期间的压裂质量。

在一些实施方式中，可以修改处理以考虑污染问题。例如，在数据采集期间，反射的管波可能会受到环境噪声(比如交通噪声、人类活动和泵噪声)的污染。这些污染可以在计算和所产生的绘图中考虑。

图5是根据本公开的一些实施方式的用于基于同步地震数据和压裂处理时间生成解释结果的方法500的示例的流程图。为了表述清晰，下面的描述总体上在本描述中的其他附图的上下文中描述了方法500。然而，应当理解，在适当时，方法500可以例如通过任何合适的系统、环境、软件和硬件、或者系统、环境、软件和硬件的组合来执行。在一些实施方式中，方法500的各个步骤可以并行、组合、循环或以任何顺序进行。

在502处，生成从三分量传感器获得的传感器读数的标准格式的地震数据集。例如，从三分量传感器获得的传感器读数包括可以小于10赫兹的信号。三分量传感器可以位于井口(比如井口102)的外部。传感器读数可以在发生在井中的压裂期间根据压裂地震测试的井轨迹的几何形状100获得。生成标准格式的地震数据集可以包括例如接收与由三分量传感器获得的传感器读数相对应的原始数据，以及通过将与传感器读数相对应的原始数据转换成标准格式来生成标准格式的地震数据集。例如，标准格式可以是勘探地球物理学家协会(SEG)Y(SEGY)格式。方法500从502前进到504。

在504处，将与传感器读数相对应的坐标和记录时间添加到标准格式的地震数据集中。作为示例，添加坐标和记录时间可以包括用坐标和记录时间来更新标准格式的地震数据集的SEGY道头。方法500从504前进到506。

在506处，通过将地震记录时间与压裂处理时间进行同步，从标准格式的地震数据集中生成同步地震数据。由于传感器放置在地面，因此同步过程很简单并且可以通过将记录时间和压裂时间都设置为GPS(全球定位系统)时间来执行。方法500从506前进到508。

在508处，通过从同步地震数据中去除死道和异常数据样本来生成质量受控的同步地震数据。异常数据可以根据其统计属性(包括但不限于均值、方差、偏斜度和峰度)自动检测。方法500从508前进到510。

在510处，对质量受控的同步地震数据执行时频分析，包括执行短时傅里叶变换以分析傅里叶谱随时间的变化。计算短时傅里叶变换的规程是将较长的时间信号分成等长的较短段，然后在每个较短段上单独计算傅里叶变换。可以使用其他时频分析工具，比如Gabor变换和S变换。方法500从510前进到512。

在512处，基于时频分析，在每个时间样本处从时频谱中提取共振频率。对于每个时间样本，该规程包括对频谱应用自动增益控制(AGC)，将具有不同间隔的频谱堆叠，最后挑选堆叠能量最高的频率间隔作为估计的共振频率，其中，应用AGC是为了使用滑动频率窗口来增加弱频率分量的幅度。对所有时间样本执行该规程。方法500从512前进到514。

在514处，将基于提取的解释结果呈现给用户。例如，呈现基于提取的解释结果可以包括生成裂缝导流能力的直方图。该信息可以在用户界面中显示给用户(例如，与水力压裂作业相关联的工程师)。在514之后，方法500可以停止。

图6是根据本公开的一些实施方式的用于提供与本公开中描述的所描述的算法、方法、功能、过程、流程和规程相关联的计算功能的示例计算机系统600的框图。所示计算机602旨在涵盖任何计算设备，比如服务器、台式计算机、膝上型/笔记本计算机、无线数据端口、智能手机、个人数据助理(PDA)、平板计算设备、或这些设备中的一个或多个处理器，包括物理实例、虚拟实例或两者兼有。计算机602可以包括可以接受用户信息的输入设备，比如小键盘、键盘和触摸屏。而且，计算机602可以包括输出设备，这些输出设备可以传送与计算机602的操作相关联的信息。该信息可以包括数字数据、视觉数据、音频信息或信息的组合。信息可以在图形用户界面(UI)(或GUI)中呈现。

计算机602可以充当客户端、网络部件、服务器、数据库、持久性或用于执行本公开中描述的主题的计算机系统的部件的角色。所示计算机602与网络630通信地耦接。在一些实施方式中，计算机602的一个或多个部件可以被配置为在不同环境内操作，包括基于云计算的环境、本地环境、全局环境以及环境的组合。

在最高层面上，计算机602是可操作以接收、发射、处理、存储和管理与所描述的主题相关联的数据和信息的电子计算设备。根据一些实施方式，计算机602还可以包括应用服务器、电子邮件服务器、web服务器、高速缓存服务器、流式数据服务器或服务器的组合，或与它们通信地耦接。

计算机602可以通过网络630从(例如，在另一个计算机602上执行的)客户端应用程序接收请求。计算机602可以通过使用软件应用程序处理接收到的请求来响应接收到的请求。请求也可以从内部用户(例如，从命令控制台)、外部(或第三)方、自动化应用程序、实体、个人、系统和计算机发送到计算机602。

计算机602的每个部件可以使用系统总线603进行通信。在一些实施方式中，计算机602的任何或所有部件(包括硬件或软件部件)可以通过系统总线603彼此或与接口604(或两者的组合)进行接口连接。接口可以使用应用编程接口(API)612、服务层613、或API612和服务层613的组合。API 612可以包括例程、数据结构和对象类的规范。API 612可以是独立于计算机语言的，也可以是依赖于计算机语言的。API 612可以指完整接口、单个功能或一组API。

服务层613可以向计算机602和可通信地耦接到计算机602的其他部件(无论是否示出)提供软件服务。使用该服务层的所有服务消费者都可以访问计算机602的功能。诸如由服务层613提供的那些软件服务可以通过定义的接口提供可重用的定义的功能。例如，接口可以是用JAVA、C++或以可扩展标记语言(XML)格式提供数据的语言编写的软件。虽然被展示为计算机602的集成部件，但在替代实施方式中，API 612或服务层613可以是相对于计算机602的其他部件和可通信地耦接到计算机602的其他部件的独立部件。此外，API 612或服务层613的任何或所有部分都可以被实施为另一个软件模块、企业应用程序或硬件模块的子模块或次模块，而不脱离本公开的范围。

计算机602包括接口604。尽管在图6中被展示为单个接口604，可以根据计算机602的特定需要、期望或特定实施方式以及所描述的功能来使用两个或更多个接口604。接口604可以由计算机602用于在分布式环境中与连接到网络630的其他系统(无论是否示出)进行通信。通常，接口604可以包括编码在可操作以与网络630通信的软件或硬件(或软件和硬件的组合)中的逻辑，或者使用该逻辑来实施。更具体地，接口604可以包括支持与通信相关联的一种或多种通信协议的软件。因此，网络630或接口的硬件可以可操作以在所示计算机602内部和外部传送物理信号。

计算机602包括处理器605。尽管在图6中被展示为单个处理器605，可以根据计算机602的特定需要、期望或特定实施方式以及所描述的功能来使用两个或更多个处理器605。通常，处理器605可以执行指令并且可以操纵数据以执行计算机602的操作，包括使用如本公开中描述的算法、方法、功能、过程、流程和规程的操作。

计算机602还包括数据库606，该数据库可以保存用于计算机602和连接到网络630的其他部件(无论是否示出)的数据。例如，数据库606可以是内存内的、常规的、或存储与本公开一致的数据的数据库。在一些实施方式中，根据计算机602的特定需要、期望或特定实施方式以及所描述的功能，数据库606可以是两种或更多种不同数据库类型的组合(例如，混合内存内数据库和常规数据库)。尽管在图6中被展示为单个数据库606，可以根据计算机602的特定需要、期望或特定实施方式以及所描述的功能来使用两个或更多个数据库(相同类型、不同类型、或类型的组合)。虽然数据库606被展示为计算机602的内部部件，但在替代实施方式中，数据库606可以在计算机602的外部。

计算机602还包括存储器607，该存储器可以保存用于计算机602或连接到网络630的部件组合(无论是否示出)的数据。存储器607可以存储与本公开一致的任何数据。在一些实施方式中，根据计算机602的特定需要、期望或特定实施方式以及所描述的功能，存储器607可以是两种或更多种不同存储器类型的组合(例如，半导体存储和磁存储的组合)。尽管在图6中被展示为单个存储器607，可以根据计算机602的特定需要、期望或特定实施方式以及所描述的功能来使用两个或更多个存储器607(相同类型、不同类型、或类型的组合)。虽然存储器607被展示为计算机602的内部部件，但在替代实施方式中，存储器607可以在计算机602的外部。

应用程序608可以是根据计算机602的特定需要、期望或特定实施方式和所描述的功能提供功能的算法软件引擎。例如，应用程序608可以用作一个或多个部件、模块或应用程序。进一步地，虽然被展示为单个应用程序608，但应用程序608可以被实施为计算机602上的多个应用程序608。另外，虽然被展示为在计算机602内部，但在替代实施方式中，应用程序608可以在计算机602外部。

计算机602还可以包括电源614。电源614可以包括可充电或不可充电电池，其可以被配置为用户可更换或非用户可更换。在一些实施方式中，电源614可以包括功率转换和管理电路，包括再充电、待机和功率管理功能。在一些实施方式中，电源614可以包括电源插头以允许将计算机602插入壁式插座或电源中，以例如为计算机602供电或对可充电电池充电。

可以有任意数量的计算机602与包含计算机602的计算机系统相关联或在其外部，其中每个计算机602通过网络630进行通信。进一步地，在不脱离本公开的范围的情况下，术语“客户端”、“用户”和其他适当的术语可以在适当时互换使用。此外，本公开设想许多用户可以使用一个计算机602并且一个用户可以使用多个计算机602。

所描述的主题实施方式可以包括单独或组合的一个或多个特征。

例如，在第一实施方式中，一种计算机实施的方法包括以下内容。生成从三分量传感器获得的传感器读数的标准格式的地震数据集。将与这些传感器读数相对应的坐标和记录时间添加到该标准格式的地震数据集中。通过将地震记录时间与压裂处理时间进行同步，从该标准格式的地震数据集中生成同步地震数据。通过从该同步地震数据中去除死道和异常数据样本来生成质量受控的同步地震数据。对该质量受控的同步地震数据执行时频分析，包括执行短时傅里叶变换以分析傅里叶谱随时间的变化。基于该时频分析，在不同的时间样本处从每个频谱中提取共振频率。将基于的解释结果呈现给用户。

前述和其他描述的实施方式可以各自可选地包括以下特征中的一个或多个：

第一特征，可与以下任何特征相结合，其中三分量传感器位于井口外部。

第二特征，可与前面或以下任何特征相结合，其中生成标准格式的地震数据集包括：接收与由该三分量传感器获得的传感器读数相对应的原始数据；以及通过将与这些传感器读数相对应的原始数据转换为标准格式来生成该标准格式的地震数据集。

第三特征，可与前面或以下任何特征相结合，其中标准格式是勘探地球物理学家协会(SEG)Y(SEGY)格式。

第四特征，可与前面或以下任何特征相结合，其中添加这些坐标和这些记录时间包括用这些坐标和记录时间来更新该标准格式的地震数据集的SEGY道头。

第五特征，可与前面或以下任何特征相结合，其中呈现基于该提取的解释结果包括生成裂缝导流能力直方图。

第六特征，可与前面或以下任何特征相结合，其中从该三分量传感器获得的传感器读数包括小于10赫兹的信号。

在第二实施方式中，一种非暂态计算机可读介质存储可由计算机系统执行以执行包括以下各项的操作的一个或多个指令。生成从三分量传感器获得的传感器读数的标准格式的地震数据集。将与这些传感器读数相对应的坐标和记录时间添加到该标准格式的地震数据集中。通过将地震记录时间与压裂处理时间进行同步，从该标准格式的地震数据集中生成同步地震数据。通过从该同步地震数据中去除死道和异常数据样本来生成质量受控的同步地震数据。对该质量受控的同步地震数据执行时频分析，包括执行短时傅里叶变换以分析傅里叶谱随时间的变化。基于该时频分析，在不同的时间样本处从每个频谱中提取共振频率。将基于的解释结果呈现给用户。

前述和其他描述的实施方式可以各自可选地包括以下特征中的一个或多个：

第一特征，可与以下任何特征相结合，其中三分量传感器位于井口外部。

第二特征，可与前面或以下任何特征相结合，其中生成标准格式的地震数据集包括：接收与由该三分量传感器获得的传感器读数相对应的原始数据；以及通过将与这些传感器读数相对应的原始数据转换为标准格式来生成该标准格式的地震数据集。

第三特征，可与前面或以下任何特征相结合，其中标准格式是勘探地球物理学家协会(SEG)Y(SEGY)格式。

第四特征，可与前面或以下任何特征相结合，其中添加这些坐标和这些记录时间包括用这些坐标和记录时间来更新该标准格式的地震数据集的SEGY道头。

第五特征，可与前面或以下任何特征相结合，其中呈现基于该提取的解释结果包括生成裂缝导流能力直方图。

第六特征，可与前面或以下任何特征相结合，其中从该三分量传感器获得的传感器读数包括小于10赫兹的信号。

在第三实施方式中，一种计算机实施的系统包括一个或多个处理器和非暂态计算机可读存储介质，该非暂态计算机可读存储介质耦接到该一个或多个处理器并存储用于由该一个或多个处理器执行的编程指令。这些编程指令指示一个或多个处理器执行包括以下各项的操作。生成从三分量传感器获得的传感器读数的标准格式的地震数据集。将与这些传感器读数相对应的坐标和记录时间添加到该标准格式的地震数据集中。通过将地震记录时间与压裂处理时间进行同步，从该标准格式的地震数据集中生成同步地震数据。通过从该同步地震数据中去除死道和异常数据样本来生成质量受控的同步地震数据。对该质量受控的同步地震数据执行时频分析，包括执行短时傅里叶变换以分析傅里叶谱随时间的变化。基于该时频分析，在不同的时间样本处从每个频谱中提取共振频率。将基于的解释结果呈现给用户。

前述和其他描述的实施方式可以各自可选地包括以下特征中的一个或多个：

第一特征，可与以下任何特征相结合，其中三分量传感器位于井口外部。

第二特征，可与前面或以下任何特征相结合，其中生成标准格式的地震数据集包括：接收与由该三分量传感器获得的传感器读数相对应的原始数据；以及通过将与这些传感器读数相对应的原始数据转换为标准格式来生成该标准格式的地震数据集。

第三特征，可与前面或以下任何特征相结合，其中标准格式是勘探地球物理学家协会(SEG)Y(SEGY)格式。

第四特征，可与前面或以下任何特征相结合，其中添加这些坐标和这些记录时间包括用这些坐标和记录时间来更新该标准格式的地震数据集的SEGY道头。

第五特征，可与前面或以下任何特征相结合，其中呈现基于该提取的解释结果包括生成裂缝导流能力直方图。

本说明书中描述的主题和功能操作的实施方式可以在数字电子电路系统中、在有形体现的计算机软件或固件中、在计算机硬件(包括在本说明书中公开的结构及其结构等同物)中、或者在其一个或多个组合中实施。所描述主题的软件实施方式可以实施为一个或多个计算机程序。每个计算机程序可以包括计算机程序指令的一个或多个模块，这些计算机程序指令编码在有形非暂态计算机可读的计算机存储介质上，以用于由数据处理装置执行或控制数据处理装置的操作。可替代地或附加地，这些程序指令可以编码在人工生成的传播信号中/上。例如，该信号可以是机器生成的电、光或电磁信号，该信号被生成以对信息进行编码以传输到合适的接收器装置以供数据处理装置执行。计算机存储介质可以是机器可读存储设备、机器可读存储基板、随机或串行存取存储器设备、或计算机存储介质的组合。

术语“数据处理装置”、“计算机”和“电子计算机设备”(或本领域普通技术人员所理解的等同物)是指数据处理硬件。例如，数据处理装置可以涵盖用于处理数据的所有种类的装置、设备和机器，包括例如可编程处理器、计算机、或多个处理器或计算机。该装置还可以包括专用逻辑电路系统，包括例如中央处理单元(CPU)、现场可编程门阵列(FPGA)或专用集成电路(ASIC)。在一些实施方式中，数据处理装置或专用逻辑电路系统(或者数据处理装置或专用逻辑电路系统的组合)可以是基于硬件或基于软件的(或者是基于硬件和基于软件的组合)。该装置可以可选地包括为计算机程序创建执行环境的代码，例如，构成处理器固件、协议栈、数据库管理系统、操作系统、或执行环境的组合的代码。本公开设想使用具有或不具有诸如LINUX、UNIX、WINDOWS、MAC OS、ANDROID或IOS等常规操作系统的数据处理装置。

计算机程序(也可以被称为或描述为程序、软件、软件应用程序、模块、软件模块、脚本或代码)可以以任何形式的编程语言来编写。编程语言可以包括例如编译语言、解释型语言、声明式语言或过程化语言。程序可以以任何形式部署，包括作为独立程序、模块、部件、子例程或单元在计算环境中使用。计算机程序可以但不必对应于文件系统中的文件。程序可以存储在持有其他程序或数据的文件的一部分(例如，在标记语言文档中存储的一个或多个脚本中)，存储在专用于所讨论程序的单个文件中，或存储在存储一个或多个模块、子程序或部分代码的多个协调文件中。可以部署计算机程序以在一个计算机上或位于例如一个站点或分布在多个站点上通过通信网络互连的多个计算机上执行。虽然各个图中所示的程序的部分可以示出为通过各种对象、方法或过程实施各种特征和功能的单独模块，但程序可以替代地包括多个子模块、第三方服务、部件和库。相反，在适当时，各种部件的特性和功能可以组合成单个部件。用于进行计算确定的阈值可以是静态地、动态地、或静态地和动态地确定的。

本说明书中描述的方法、过程或逻辑流程可以由一个或多个可编程计算机来执行，该一个或多个可编程计算机执行一个或多个计算机程序以通过对输入数据进行操作并生成输出来执行功能。方法、过程或逻辑流程也可以由专用逻辑电路系统执行，并且装置也可以实施为专用逻辑电路系统，例如CPU、FPGA或ASIC。

适合执行计算机程序的计算机可以基于一个或多个通用和专用微处理器以及其他种类的CPU。计算机的元件是用于实施或执行指令的CPU以及用于存储指令和数据的一个或多个存储器设备。通常，CPU可以从存储器接收指令和数据(并将数据写入)。计算机还可以包括或可操作地耦接到用于存储数据的一个或多个大容量存储设备。在一些实施方式中，计算机可以从包括例如磁盘、磁光盘或光盘等大容量存储设备接收数据并将数据转移到大容量存储设备。此外，计算机可以嵌入到另一个设备中，例如移动电话、个人数字助理(PDA)、移动音频或视频播放器、游戏控制台、全球定位系统(GPS)接收器或便携式存储设备，比如通用串行总线(USB)闪存驱动器。

适用于存储计算机程序指令和数据的计算机可读介质(暂态或非暂态，在适当时)可以包括所有形式的永久/非永久和易失性/非易失性存储器、介质和存储器设备。计算机可读介质可以包括例如半导体存储器设备，比如随机存取存储器(RAM)、只读存储器(ROM)、相变存储器(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、可擦除可编程只读存储器(EPROM)、电可擦除可编程只读存储器(EEPROM)和闪速存储器设备。计算机可读介质还可以包括例如磁性设备，比如磁带、盒式磁带、磁带盒和内部/可移除磁盘。计算机可读介质还可以包括磁光盘和光学存储器设备和技术，包括例如数字视频盘(DVD)、CD-ROM、DVD+/-R、DVD-RAM、DVD-ROM、HD-DVD和蓝光。存储器可以存储各种对象或数据，包括高速缓存、类、框架、应用程序、模块、备份数据、作业、网页、网页模板、数据结构、数据库表、储存库和动态信息。存储在存储器中的对象和数据的类型可以包括参数、变量、算法、指令、规则、约束和引用。此外，存储器可以包括日志、策略、安全或访问数据以及报告文件。处理器和存储器可以由专用逻辑电路系统补充或结合到专用逻辑电路系统中。

本公开中描述的主题的实施方式可以在具有显示设备的计算机上实施，该显示设备用于提供与用户的交互，包括向用户显示信息(以及从用户接收输入)。显示设备的类型可以包括例如阴极射线管(CRT)、液晶显示器(LCD)、发光二极管(LED)和等离子监视器。显示设备可以包括键盘和定点设备，包括例如鼠标、轨迹球或触控板。也可以通过使用触摸屏向计算机提供用户输入，比如具有压力敏感性的平板计算机表面或使用电容或电感测的多点触摸屏。其他种类的设备可以用于提供与用户的交互，包括接收用户反馈，包括例如感官反馈，包括视觉反馈、听觉反馈或触觉反馈。可以以声音、语音或触觉输入的形式接收来自用户的输入。另外，计算机可以通过向用户使用的设备发送文档和从用户使用的设备接收文档来与用户交互。例如，计算机可以响应从用户客户端设备上的网络浏览器接收到的请求，将网页发送到该网络浏览器。

可以以单数或复数形式使用术语“图形用户界面”或“GUI”来描述一个或多个图形用户界面以及特定图形用户界面的每个显示。因此，GUI可以表示处理信息并高效地将信息结果呈现给用户的任何图形用户界面，包括但不限于网络浏览器、触摸屏或命令行界面(CLI)。一般而言，GUI可以包括多个用户界面(UI)元素，其中一些或全部与网络浏览器相关联，比如交互式字段、下拉列表和按钮。这些和其他UI元素可以与网络浏览器的功能相关或表示其功能。

本说明书中描述的主题的实施方式可以在包括后端部件(例如，作为数据服务器)或包括中间件部件(例如，应用服务器)的计算系统中实施。此外，计算系统可以包括前端部件，例如，具有图形用户界面或网络浏览器之一或两者(用户可以通过其与计算机交互)的客户端计算机。系统的部件可以通过通信网络中的任何形式或媒介的有线或无线数字数据通信(或数据通信的组合)来互连。通信网络的示例包括局域网(LAN)、无线接入网(RAN)、城域网(MAN)、广域网(WAN)、全球微波接入互操作性(WIMAX)、无线局域网(WLAN)(例如，使用802.11a/b/g/n或802.20或协议组合)、互联网的全部或部分，或一个或多个位置处的任何其他一个或多个通信系统(或通信网络的组合)。网络可以与例如互联网协议(IP)分组、帧中继帧、异步传输模式(ATM)信元、语音、视频、数据或网络地址之间的通信类型的组合进行通信。

计算系统可以包括客户端和服务器。客户端和服务器通常可以彼此远离，并且通常可以通过通信网络进行交互。客户端和服务器的关系可以通过在各自计算机上运行并具有客户端-服务器关系的计算机程序而产生。

集群文件系统可以是可从多个服务器访问以进行读取和更新的任何文件系统类型。锁定或一致性跟踪可能不是必需的，因为交换文件系统的锁定可以在应用层完成。此外，统一码(Unicode)数据文件可以不同于非统一码数据文件。

虽然本说明书包含许多特定实施方式细节，但这些不应被解释为对可能要求保护的事物的范围的限制，而是被解释为对可能特定于特定实施方式的特征的描述。在单独的实施方式的背景下在本说明书中所描述的某些特征还可以组合地在单个实施方式中实施。与此相反，在单个实施方式的背景下描述的不同特征也可以单独地或以任何适合的子组合形式在多个实施方式中实施。此外，尽管先前描述的特征可以被描述为在某些组合中起作用并且甚至最初也是如此要求保护的，但是在一些情况下，可以从组合中除去来自所要求保护的组合的一个或多个特征，并且所要求保护的组合可以涉及子组合或子组合的变化。

已经描述了主题的特定实施方式。所描述的实施方式的其他实施方式、变更和排列在所附权利要求的范围内，这对于本领域的技术人员将是显而易见的。虽然附图或权利要求中以具体顺序描绘了操作，但这不应被理解成要求这种操作以所示的具体顺序或以有序顺序执行，或者要求所有展示的操作(一些操作可以被认为是可选的)都被执行，才能实现期望的结果。在某些情况下，多任务处理或并行处理(或多任务处理和并行处理的组合)可能是有利的，并在被认为适当的情况下执行。

此外，前述实施方式中的各种系统模块和部件的分离或集成不应被理解为在所有实施方式中都需要这种分离或集成。应当理解，所描述的程序部件和系统通常可以一起集成在单个软件产品中或封装成多个软件产品。

因此，先前描述的示例实施方式不定义或限制本公开。在不脱离本公开的精神和范围的情况下，其他改变、替代和变更也是可能的。

此外，任何要求保护的实施方式都被认为适用于至少一种计算机实施的方法；存储用于执行计算机实施的方法的计算机可读指令的非暂态计算机可读介质；以及计算机系统，该计算机系统包括与硬件处理器可互操作地耦接的计算机存储器，该硬件处理器被配置为执行计算机实施的方法或存储在非暂态计算机可读介质上的指令。