CN114578445A - 基于重力资料确定断裂位置的方法及装置 - Google Patents

Abstract

本发明提供一种基于重力资料确定断裂位置的方法及装置，方法包括：获取重力资料的基础数据布格重力异常数据；由布格重力异常网格化数据获得剩余重力异常网格化数据，对剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；对伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；根据伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置。本发明能够使得利用重力资料研究断裂的水平梯度极大值条带宽度降低一半，从而提高了断裂位置的确定精度，同时提高了规模小、深度大的断裂的识别能力。

Description

基于重力资料确定断裂位置的方法及装置

技术领域

本发明涉及重力勘探资料处理技术领域，具体涉及一种基于重力资料确定断裂位置的方法及装置。

背景技术

断裂是自然界的一种基本地质现象，人类活动及能源矿产勘探开发与之密切相关，有关技术人员一直在探索提高断裂位置解释精度的方法。由于不同地层和岩性的岩石密度不同，造成断裂两侧存在重力差异，使得重力资料研究断裂具备物理基础。重力方法研究人员一直在探索不断提高重力资料确定断裂位置精度的新方法。布格重力异常是重力勘探的基础资料，断裂在布格重力异常和从布格重力异常中去除区域重力场后得到的剩余重力异常上表现为重力梯度变化带(有时简称为重力梯级带)，人们可以通过布格重力异常或剩余重力异常的梯度变化带的位置了解断裂的大致位置，但由于重力梯度变化带存在一定宽度，造成无法直接从该带上判读断裂的准确位置。为了获得较准确的断裂位置，重力研究工作者不断探索处理解释方法，归纳起来，可以分为两个步骤，第一步是对布格重力异常或剩余重力异常进行资料处理，获得可以突出反映断裂位置的重力异常，传统和代表性方法是进行水平导数计算，得到重力水平导数图(也称为重力水平梯度图)，第二步是绘制反映断裂位置的重力异常图，在图上识别断裂位置，对重力水平导数图而言，是在图上识别极大值位置，将极大值位置连线解释为断裂位置。上述方法是正确有效的，但存在断裂位置精度不理想的问题。

由于重力梯度变化带宽度较大，使得水平导数的极大值虽然明显收窄在一个相对较狭窄的范围，但从这个较狭窄的带中确定断裂的准确位置仍然存在精度不理想的问题。为更好地解决该问题，研究人员采取了多种技术措施。在上述第一步中采用的措施有几类：第一类是对布格重力异常或剩余重力异常进行小子域滤波、小波变换等处理，使重力梯度变化带宽度收窄，然后再进行水平导数计算；第二类是对布格重力异常或剩余重力异常计算重力水平总梯度、重力总梯度等，取代重力水平导数，使梯度极大值带的范围收窄；第三类是在计算水平导数后再求取垂直导数、高通滤波等，使极大值带的范围进一步收窄。这三类方法均属于对布格重力异常或剩余重力异常求取水平导数的方法，主要是在求取导数前或后进行了信息增强，一些研究人员也把上述方法中的两种以上方法联合使用以增强处理效果。在第二步中采用的措施有几类：第一类是采用影像法，即把水平梯度数据模拟为类似三维起伏地形，再以一定高度角的光照射，光线明暗变化的界限是极大值的位置，解释为断裂的位置；第二类是寻值法，通过数学运算，寻找极大值的位置，以极大值连线作为断裂位置。这两类方法均是在第一步获得的重力异常数据的基础上进行的，它可以提高断裂解释精度，但其精度受限于第一步获得的重力异常资料的精度。

如何收窄水平梯度的极值带宽度，是提高基于重力资料确定断裂位置的精度的亟待解决的技术问题。

发明内容

针对现有技术中的问题，本发明提供一种基于重力资料确定断裂位置的方法及装置，提高确定断裂位置的精度。

为解决上述技术问题，本发明提供以下技术方案：

第一方面，本发明提供一种基于重力资料确定断裂位置的方法，包括：

获取重力资料的基础数据布格重力异常数据；

由所述布格重力异常网格化数据获得剩余重力异常网格化数据，对所述剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；

对所述伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；

根据所述伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；

在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置。

其中，在所述获取重力资料的基础数据布格重力异常数据之后，还包括：

对所述布格重力异常数据进行网格化插值处理，得到布格重力异常网格化数据；

确定布格重力异常网格化数据中存在区域背景重力场，则由所述布格重力异常网格化数据获得剩余重力异常网格化数据。

其中，所述由所述布格重力异常网格化数据获得剩余重力异常网格化数据，包括：

采用位场延拓法或波数域滤波法由所述布格重力异常网格化数据获得剩余重力异常网格化数据。

其中，所述对所述剩余重力异常网格化数据进行伪磁力异常转换，包括：

利用同一均匀地质体的磁位与引力位的泊松公式，对剩余重力异常网格化数据进行伪磁力异常转换。

第二方面，本发明提供一种基于重力资料确定断裂位置的装置，包括：

获取单元，用于获取重力资料的基础数据布格重力异常数据；

提取单元，用于由所述布格重力异常网格化数据获得剩余重力异常网格化数据，对所述剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；

计算单元，用于对所述伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；

绘图单元，用于根据所述伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；

确定单元，用于在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置。

其中，还包括：

插值单元，用于对所述布格重力异常数据进行网格化插值处理，得到布格重力异常网格化数据；

重力场单元，用于确定布格重力异常网格化数据中存在区域背景重力场，则由所述布格重力异常网格化数据获得剩余重力异常网格化数据。

其中，所述提取单元包括：

第一提取子单元，用于采用位场延拓法或波数域滤波法由所述布格重力异常网格化数据获得剩余重力异常网格化数据。

其中，所述提取单元包括：

第二提取子单元，用于利用同一均匀地质体的磁位与引力位的泊松公式，对剩余重力异常网格化数据进行伪磁力异常转换。

第三方面，本发明提供一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述程序时实现所述的基于重力资料确定断裂位置的方法的步骤。

第四方面，本发明提供一种计算机可读存储介质，其上存储有计算机程序，该计算机程序被处理器执行时实现所述的基于重力资料确定断裂位置的方法的步骤。

由上述技术方案可知，本发明提供一种基于重力资料确定断裂位置的方法及装置，通过获取重力资料的基础数据布格重力异常数据；由所述布格重力异常网格化数据获得剩余重力异常网格化数据，对所述剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；对所述伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；根据所述伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置，能够使得利用重力资料研究断裂的水平梯度极大值条带宽度降低一半，从而提高了断裂位置的确定精度，同时提高了规模小、深度大的断裂的识别能力。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明实施例中的基于重力资料确定断裂位置的方法的第一流程示意图。

图2为本发明实施例中的基于重力资料确定断裂位置的方法的第二流程示意图。

图3为本发明实施例中的基于重力资料确定断裂位置的方法中台阶模型示意图。

图4为本发明实施例中的基于重力资料确定断裂位置的方法中台阶模型重力异常曲线示意图。

图5为本发明实施例中的基于重力资料确定断裂位置的方法中台阶模型伪磁力异常曲线示意图。

图6为本发明实施例中的基于重力资料确定断裂位置的方法中台阶模型伪磁力异常水平导数曲线及半极值宽度示意。

图7为本发明实施例中的基于重力资料确定断裂位置的方法中台阶模型重力异常水平导数曲线及半极值宽度示意。

图8为本发明实施例中的基于重力资料确定断裂位置的方法中布格重力异常图示意图。

图9为本发明实施例中的基于重力资料确定断裂位置的方法中伪磁力异常图示意图。

图10为本发明实施例中的基于重力资料确定断裂位置的方法中伪磁力异常水平总梯度图示意图。

图11为本发明实施例中的基于重力资料确定断裂位置的方法中布格重力异常水平总梯度图示意图。

图12为本发明实施例中的基于重力资料确定断裂位置的结构示意图。

图13为本发明实施例中的电子设备的结构示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚，下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整的描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明提供一种基于重力资料确定断裂位置的方法的实施例，参见图1，所述基于重力资料确定断裂位置的方法具体包含有如下内容：

S101：获取重力资料的基础数据布格重力异常数据；

在本步骤中，在研究区部署重力勘探，进行重力资料采集，获得测点坐标和重力值，进行正常场改正、高程改正、中间层改正、地形改正等各项改正计算，得到布格重力异常数据。

S102：由所述布格重力异常网格化数据获得剩余重力异常网格化数据，对所述剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；

在本步骤中，采用位场延拓法或波数域滤波法由所述布格重力异常网格化数据获得剩余重力异常网格化数据。

其中，布格重力异常数据包括了布格重力异常网格化数据。

利用同一均匀地质体的磁位与引力位的泊松公式，对剩余重力异常网格化数据进行伪磁力异常转换。

本步骤中的，磁位与引力位的泊松公式为：

式中U为地质体的磁位，V为地质体的引力位，J为磁化强度，G为引力常数，δ为密度。

S103：对所述伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；

在本步骤中，水平导数可以是水平方向导数、水平总梯度或总梯度。

S104：根据所述伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；

在本步骤中，图件可以是平面等值线图或者三维影像图。

S105：在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置。

在本步骤中，在伪磁力异常水平导数图上识别极大值位置，连线各相关的极大值位置，得到断裂平面位置，实现重力资料精确确定断裂位置。

从上述描述可知，本发明提供一种基于重力资料确定断裂位置的方法，通过获取重力资料的基础数据布格重力异常数据；由所述布格重力异常网格化数据获得剩余重力异常网格化数据，对所述剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；对所述伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；根据所述伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置，能够使得利用重力资料研究断裂的水平梯度极大值条带宽度降低一半，从而提高了断裂位置的确定精度，同时提高了规模小、深度大的断裂的识别能力。

在本发明的一实施例中，参见图2，所述基于重力资料确定断裂位置的方法的步骤S101之后，具体包含有如下内容：

S106：对所述布格重力异常数据进行网格化插值处理，得到布格重力异常网格化数据；

在本步骤中网格间距取重力点距，网格化范围为数据范围。本步骤的目的是为后续的处理工作按需要的格式准备数据。

S107：确定布格重力异常网格化数据中存在区域背景重力场，则由所述布格重力异常网格化数据获得剩余重力异常网格化数据。

在本步骤中，结合已知历史地质资料和历史勘探成果，确定布格重力异常中是否存在区域背景重力场。如果存在，则对布格重力异常网格化数据进行剩余重力异常提取处理，获得剩余重力异常网格化数据；如果不存在明显区域背景场，则可不做剩余重力异常提取处理。

剩余重力异常的提取可采用位场延拓法、波数域滤波法等各种有效方法，求得的剩余重力异常应最大程度地压制区域背景场。

从上述描述可知，本实施例提供的技术方案，采用重力资料确定断裂位置的方法，使重力资料确定的断裂位置更准确。

本发明实施例提供一种应用基于重力资料确定断裂位置的方法的具体实例，具体包括如下内容：

下面分别展示利用模型重力剖面数据和实测重力平面数据确定断裂位置的实施步骤和方法效果。

首先展示模型重力剖面数据确定断裂位置的实施步骤和效果。

模型选用垂直台阶二维模型，台阶顶面埋深500米，底面埋深1000米，台阶(断裂)位置在中部0米处，两侧密度差1.0g/cm³，见图3。重力计算剖面垂直穿过台阶。通过重力正演计算，得到了上述模型的重力异常曲线，见图4，这相当于实施本发明的步骤一，获取重力异常数据。

对重力剖面研究而言，不需要进行数据网格化，因此不实施步骤二。

在设置模型时，没有叠加区域背景场，因此不实施步骤三。

实施步骤四，把模型重力异常(图4)转换为伪磁力异常，得到模型伪磁力异常曲线(图5)。在转换中，磁场强度设为50000nT，磁倾角设为90°，相当于垂直磁化，磁化率设为1000×10^-5SI。对比图5和图4可以看出，在剖面中部对应台阶(断裂)位置0米附近，伪磁力异常梯度变化带宽度(图7)比重力异常梯度变化带(图4)要窄很多，主要局限在中部较小的范围内。

实施步骤五，对伪磁力异常求取水平导数，得到伪磁力异常水平导数数据。

实施步骤六，以伪磁力异常水平导数数据绘制伪磁力异常水平导数曲线(图6)。按照重磁力研究中经常采用的半极值异常宽度的指标进行分析，模型伪磁力异常水平导数半极值异常宽度约为700米。

实施步骤七，对伪磁力异常水平导数极大值位置进行解释，极大值位置就是断裂的位置。从曲线上看，极大值位于0米上，指示了断裂的位置在0米附近，与台阶模型中台阶(断裂)的位置一致。从而实现了台阶模型重力数据识别断裂位置。

为了与重力异常水平导数确定断裂位置的现有方法进行效果对比，对模型重力异常(图4)直接求取了水平导数，重力异常水平导数曲线见图7。重力异常的水平导数极大值位于0米，与伪磁力异常水平导数反映的位置一致，说明本发明的方法是正确的；重力异常水平导数的半极值异常宽度约为1400米，是伪磁力异常水平导数半极值异常宽度的2倍，说明本发明的方法精度更高。可见，伪磁力异常水平导数对断裂位置的反映明显好于重力异常水平导数对断裂位置的反映，即本发明的方法优于直接利用重力异常水平导数确定断裂位置的现有方法。

然后展示实测重力平面数据确定断裂位置的实施步骤和方法效果。研究区选在某盆地的盆地边界处，断裂斜穿研究区。

实施步骤一，在研究区部署了500米×500米的重力测网，通过重力采集，得到了各测点坐标和重力值，进行正常场改正、高程改正、中间层改正、地形改正等各项改正计算，得到了布格重力异常数据。

实施步骤二，对布格重力异常数据做网格化插值处理，得到布格重力异常网格化数据。其中网格间距取重力点距，即500米，网格化范围为数据范围，图中空白处为未部署重力测点的地区，网格化方法采用现有的克里京插值方法。布格重力异常图示意图见图8。

实施步骤三，收集和分析区域地质资料知道，本区处于平稳地台区，理论上本区不存在明显区域背景场，通过收集和分析大范围重力资料知道，本区不存在明显的区域背景场，因此本实施例不需要作剩余重力异常提取，本区布格重力异常就相当于剩余重力异常，因此后续处理基于布格重力异常数据进行。

实施步骤四，利用同一均匀地质体磁位与引力位的泊松公式

对布格重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据。各转换参数按测区实际参数选取。伪磁力异常图见图9。

实施步骤五，对伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据。这里采用的水平梯度计算方法是求取水平总梯度，与后面用于对比效果的重力异常处理为同样的方法，目的是便于效果对比。

实施步骤六，采用绘图软件，以伪磁力异常水平导数网格化数据绘制伪磁力异常水平导数平面等值线图，见图10。以网格化数据绘制平面等值线图是常规方法，现有各绘图软件均可实现，这里的图形形式是最常见的平面等值线图。

实施步骤七，在伪磁力异常水平总梯度图(图10)上识别伪磁力异常水平导数的极大值位置，连线相关极大值位置就得到了断裂的平面位置，从而实现了基于重力资料确定断裂位置。这种以伪磁力异常水平总梯度极大值位置确定断裂位置的方法与以重力异常水平总梯度极大值位置确定断裂位置的做法相同，因此图中未标绘断裂位置。

为了与重力异常水平导数确定断裂位置的现有方法进行效果对比，采用常规方法对实测重力异常(图8)直接计算了水平总梯度，绘制了重力异常水平总梯度图(图11)。利用重力异常水平总梯度资料解释断裂位置的方法是从水平总梯度图上确定极大值的位置并解释为断裂的位置，因此，当极大值条带较窄时，确定的极大值位置就更准确，解释的断裂位置精度就更高，反之，当极大值条带较宽时，最终解释的断裂位置精度就会较低。对比伪磁力异常水平总梯度图(图10)和重力异常水平总梯度图(图11)可以看出，伪磁力异常水平总梯度极大值条带要明显比重力异常水平总梯度极大值条带窄，因此利用伪磁力异常水平总梯度解释的断裂位置必然比利用重力异常水平总梯度解释的断裂位置更准确。另外可以看到，在伪磁力异常水平总梯度图(图10)上，在极大值条带(反映断裂位置)中部存在较明显的断开和相对错动，而在重力异常水平总梯度图(图11)上则难以看出该现象，这进一步反映了本发明对断裂的分辨能力明显好于直接利用重力异常水平总梯度确定断裂位置的现有方法。

从上述描述可知，本发明实施例提供的基于重力资料确定断裂位置的方法，采用该项技术，使得利用重力资料研究断裂的水平梯度极大值条带宽度降低一半，从而提高了断裂位置的确定精度，同时提高了规模小、深度大的断裂的识别能力。

本发明实施例提供一种能够实现所述基于重力资料确定断裂位置的方法中全部内容的基于重力资料确定断裂位置的装置的具体实施方式，参见图12，所述基于重力资料确定断裂位置的装置具体包括如下内容：

获取单元10，用于获取重力资料的基础数据布格重力异常数据；

提取单元20，用于由所述布格重力异常网格化数据获得剩余重力异常网格化数据，对所述剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；

计算单元30，用于对所述伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；

绘图单元40，用于根据所述伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；

确定单元50，用于在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置。

其中，还包括：

插值单元，用于对所述布格重力异常数据进行网格化插值处理，得到布格重力异常网格化数据；

重力场单元，用于确定布格重力异常网格化数据中存在区域背景重力场，则由所述布格重力异常网格化数据获得剩余重力异常网格化数据。

其中，所述提取单元包括：

第一提取子单元，用于采用位场延拓法或波数域滤波法由所述布格重力异常网格化数据获得剩余重力异常网格化数据。

其中，所述提取单元包括：

第二提取子单元，用于利用同一均匀地质体的磁位与引力位的泊松公式，对剩余重力异常网格化数据进行伪磁力异常转换。

本发明提供的基于重力资料确定断裂位置的装置的实施例具体可以用于执行上述实施例中的基于重力资料确定断裂位置的方法的实施例的处理流程，其功能在此不再赘述，可以参照上述方法实施例的详细描述。

从上述描述可知，本发明实施例提供的基于重力资料确定断裂位置的装置，使得利用重力资料研究断裂的水平梯度极大值条带宽度降低一半，从而提高了断裂位置的确定精度，同时提高了规模小、深度大的断裂的识别能力。

本申请提供一种用于实现所述基于重力资料确定断裂位置方法中的全部或部分内容的电子设备的实施例所述电子设备具体包含有如下内容：

处理器(processor)、存储器(memory)、通信接口(Communications Interface)和总线；其中，所述处理器、存储器、通信接口通过所述总线完成相互间的通信；所述通信接口用于实现相关设备之间的信息传输；该电子设备可以是台式计算机、平板电脑及移动终端等，本实施例不限于此。在本实施例中，该电子设备可以参照实施例用于实现所述基于重力资料确定断裂位置方法的实施例及用于实现所述基于重力资料确定断裂位置装置的实施例进行实施，其内容被合并于此，重复之处不再赘述。

图13为本申请实施例的电子设备9600的系统构成的示意框图。如图13所示，该电子设备9600可以包括中央处理器9100和存储器9140；存储器9140耦合到中央处理器9100。值得注意的是，该图13是示例性的；还可以使用其他类型的结构，来补充或代替该结构，以实现电信功能或其他功能。

一实施例中，基于重力资料确定断裂位置功能可以被集成到中央处理器9100中。其中，中央处理器9100可以被配置为进行如下控制：

获取重力资料的基础数据布格重力异常数据；由所述布格重力异常网格化数据获得剩余重力异常网格化数据，对所述剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；对所述伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；根据所述伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置。

从上述描述可知，本申请的实施例提供的电子设备，通过获取重力资料的基础数据布格重力异常数据；由所述布格重力异常网格化数据获得剩余重力异常网格化数据，对所述剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；对所述伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；根据所述伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置，能够使得利用重力资料研究断裂的水平梯度极大值条带宽度降低一半，从而提高了断裂位置的确定精度，同时提高了规模小、深度大的断裂的识别能力。

在另一个实施方式中，基于重力资料确定断裂位置装置可以与中央处理器9100分开配置，例如可以将基于重力资料确定断裂位置配置为与中央处理器9100连接的芯片，通过中央处理器的控制来实现基于重力资料确定断裂位置功能。

如图13所示，该电子设备9600还可以包括：通信模块9110、输入单元9120、音频处理器9130、显示器9160、电源9170。值得注意的是，电子设备9600也并不是必须要包括图13中所示的所有部件；此外，电子设备9600还可以包括图13中没有示出的部件，可以参考现有技术。

如图13所示，中央处理器9100有时也称为控制器或操作控件，可以包括微处理器或其他处理器装置和/或逻辑装置，该中央处理器9100接收输入并控制电子设备9600的各个部件的操作。

其中，存储器9140，例如可以是缓存器、闪存、硬驱、可移动介质、易失性存储器、非易失性存储器或其它合适装置中的一种或更多种。可储存上述与失败有关的信息，此外还可存储执行有关信息的程序。并且中央处理器9100可执行该存储器9140存储的该程序，以实现信息存储或处理等。

输入单元9120向中央处理器9100提供输入。该输入单元9120例如为按键或触摸输入装置。电源9170用于向电子设备9600提供电力。显示器9160用于进行图像和文字等显示对象的显示。该显示器例如可为LCD显示器，但并不限于此。

该存储器9140可以是固态存储器，例如，只读存储器(ROM)、随机存取存储器(RAM)、SIM卡等。还可以是这样的存储器，其即使在断电时也保存信息，可被选择性地擦除且设有更多数据，该存储器的示例有时被称为EPROM等。存储器9140还可以是某种其它类型的装置。存储器9140包括缓冲存储器9141(有时被称为缓冲器)。存储器9140可以包括应用/功能存储部9142，该应用/功能存储部9142用于存储应用程序和功能程序或用于通过中央处理器9100执行电子设备9600的操作的流程。

存储器9140还可以包括数据存储部9143，该数据存储部9143用于存储数据，例如联系人、数字数据、图片、声音和/或任何其他由电子设备使用的数据。存储器9140的驱动程序存储部9144可以包括电子设备的用于通信功能和/或用于执行电子设备的其他功能(如消息传送应用、通讯录应用等)的各种驱动程序。

通信模块9110即为经由天线9111发送和接收信号的发送机/接收机9110。通信模块(发送机/接收机)9110耦合到中央处理器9100，以提供输入信号和接收输出信号，这可以和常规移动通信终端的情况相同。

基于不同的通信技术，在同一电子设备中，可以设置有多个通信模块9110，如蜂窝网络模块、蓝牙模块和/或无线局域网模块等。通信模块(发送机/接收机)9110还经由音频处理器9130耦合到扬声器9131和麦克风9132，以经由扬声器9131提供音频输出，并接收来自麦克风9132的音频输入，从而实现通常的电信功能。音频处理器9130可以包括任何合适的缓冲器、解码器、放大器等。另外，音频处理器9130还耦合到中央处理器9100，从而使得可以通过麦克风9132能够在本机上录音，且使得可以通过扬声器9131来播放本机上存储的声音。

本发明的实施例还提供能够实现上述实施例中的基于重力资料确定断裂位置的方法中全部步骤的一种计算机可读存储介质，所述计算机可读存储介质上存储有计算机程序，该计算机程序被处理器执行时实现上述实施例中的基于重力资料确定断裂位置的方法的全部步骤，例如，所述处理器执行所述计算机程序时实现下述步骤：

获取重力资料的基础数据布格重力异常数据；由所述布格重力异常网格化数据获得剩余重力异常网格化数据，对所述剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；对所述伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；根据所述伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置。

从上述描述可知，本发明实施例提供的计算机可读存储介质，通过获取重力资料的基础数据布格重力异常数据；由所述布格重力异常网格化数据获得剩余重力异常网格化数据，对所述剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；对所述伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；根据所述伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置，能够使得利用重力资料研究断裂的水平梯度极大值条带宽度降低一半，从而提高了断裂位置的确定精度，同时提高了规模小、深度大的断裂的识别能力。

虽然本发明提供了如实施例或流程图所述的方法操作步骤，但基于常规或者无创造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤执行顺序中的一种方式，不代表唯一的执行顺序。在实际中的装置或客户端产品执行时，可以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处理的环境)。

本领域技术人员应明白，本说明书的实施例可提供为方法、装置(系统)或计算机程序产品。因此，本说明书实施例可采用完全硬件实施例、完全软件实施例或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、装置(系统)和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

本说明书中的各个实施例均采用递进的方式描述，各个实施例之间相同相似的部分互相参见即可，每个实施例重点说明的都是与其他实施例的不同之处。尤其，对于系统实施例而言，由于其基本相似于方法实施例，所以描述的比较简单，相关之处参见方法实施例的部分说明即可。在本文中，诸如第一和第二等之类的关系术语仅仅用来将一个实体或者操作与另一个实体或操作区分开来，而不一定要求或者暗示这些实体或操作之间存在任何这种实际的关系或者顺序。需要说明的是，在不冲突的情况下，本发明中的实施例及实施例中的特征可以相互组合。本发明并不局限于任何单一的方面，也不局限于任何单一的实施例，也不局限于这些方面和/或实施例的任意组合和/或置换。而且，可以单独使用本发明的每个方面和/或实施例或者与一个或更多其他方面和/或其实施例结合使用。

最后应说明的是：以上各实施例仅用以说明本发明的技术方案，而非对其限制；尽管参照前述各实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分或者全部技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的范围，其均应涵盖在本发明的权利要求和说明书的范围当中。

Claims (10)

1.一种基于重力资料确定断裂位置的方法，其特征在于，包括：

获取重力资料的基础数据布格重力异常数据；

由所述布格重力异常网格化数据获得剩余重力异常网格化数据，对所述剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；

对所述伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；

根据所述伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；

在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置。

2.根据权利要求1所述的基于重力资料确定断裂位置的方法，其特征在于，在所述获取重力资料的基础数据布格重力异常数据之后，还包括：

对所述布格重力异常数据进行网格化插值处理，得到布格重力异常网格化数据；

确定布格重力异常网格化数据中存在区域背景重力场，则由所述布格重力异常网格化数据获得剩余重力异常网格化数据。

3.根据权利要求1所述的基于重力资料确定断裂位置的方法，其特征在于，所述由所述布格重力异常网格化数据获得剩余重力异常网格化数据，包括：

采用位场延拓法或波数域滤波法由所述布格重力异常网格化数据获得剩余重力异常网格化数据。

4.根据权利要求1所述的基于重力资料确定断裂位置的方法，其特征在于，所述对所述剩余重力异常网格化数据进行伪磁力异常转换，包括：

利用同一均匀地质体的磁位与引力位的泊松公式，对剩余重力异常网格化数据进行伪磁力异常转换。

5.一种基于重力资料确定断裂位置的装置，其特征在于，包括：

获取单元，用于获取重力资料的基础数据布格重力异常数据；

提取单元，用于由所述布格重力异常网格化数据获得剩余重力异常网格化数据，对所述剩余重力异常网格化数据进行伪磁力异常转换，得到伪磁力异常网格化数据；

计算单元，用于对所述伪磁力异常网格化数据进行水平导数计算，得到伪磁力异常水平导数网格化数据；

绘图单元，用于根据所述伪磁力异常水平导数网格化数据来绘制伪磁力异常水平导数图；

确定单元，用于在伪磁力异常水平导数图上识别极大值位置，连线各个极大值位置，得到断裂位置。

6.根据权利要求5所述的基于重力资料确定断裂位置的装置，其特征在于，还包括：

插值单元，用于对所述布格重力异常数据进行网格化插值处理，得到布格重力异常网格化数据；

重力场单元，用于确定布格重力异常网格化数据中存在区域背景重力场，则由所述布格重力异常网格化数据获得剩余重力异常网格化数据。

7.根据权利要求5所述的基于重力资料确定断裂位置的装置，其特征在于，所述提取单元包括：

第一提取子单元，用于采用位场延拓法或波数域滤波法由所述布格重力异常网格化数据获得剩余重力异常网格化数据。

8.根据权利要求5所述的基于重力资料确定断裂位置的装置，其特征在于，所述提取单元包括：

第二提取子单元，用于利用同一均匀地质体的磁位与引力位的泊松公式，对剩余重力异常网格化数据进行伪磁力异常转换。

9.一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，其特征在于，所述处理器执行所述程序时实现权利要求1至4任一项所述的基于重力资料确定断裂位置的方法的步骤。

10.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，该计算机程序被处理器执行时实现权利要求1至4任一项所述的基于重力资料确定断裂位置的方法的步骤。

Images