一种重力异常边界增强的方法及装置的制造方法
【专利摘要】本申请实施例提供的一种重力异常边界增强的方法及装置。所述方法包括:获取目标工区的重力异常数据,并对所述重力异常数据进行网格化;采用空间域算子逐个提取第一网格点在预设方向的重力异常边界极值信息;基于所述第一网格点的重力异常边界极值信息,计算第二网格点的相干增强值。本申请提供的重力异常边界增强方法对所述重力异常数据进行空间域处理,避免传统频率域算法对数据噪声的放大作用,可以获取清晰、连续的重力异常边界信息。
【专利说明】
一种重力异常边界増强的方法及装置
技术领域
[0001] 本发明涉及石油和固体矿产勘探中的地球物理勘探技术领域,尤其涉及一种重力 异常边界增强的方法及装置。
【背景技术】
[0002] 在地球物理勘探技术领域,通过重力勘探可以了解地质场源的分布,其中包括场 源的形态和规模等。重力异常是地下密度体不均匀分布的综合反应,利用重力异常资料中 的异常边线性梯级带、异常的扭曲、异常的分区特征等信息数据可以识别异常地质边界,比 如地层断裂位置等。而与地质边界有关的信息通常表现为高频率、低幅值,往往很难直接从 布格重力异常图中提取出地质边界的水平位置。通常需要先对重力异常数据的边界信息进 行增强后再提取异常地质边界信息。
[0003] 现有技术中的异常边界增强的方法主要有导数类算法,所述导数类算法主要包括 垂向导数法、总水平导数法、总梯度模法、斜导数法、西塔(Theta)图法等。其中,所述总水平 导数法可以采用对横轴和竖轴方向水平导数的平方和进行开方的方法求得。导数类算法的 共同特点在于利用垂向导数或导数的组合对信号中的高频成分进行放大,使重力异常图中 的极大值或拐点近似位于地质体边界的正上方,后续可以通过边缘检测的方法提取出异常 地质边界位置。
[0004]导数类算法一般采用频率域傅氏变换的方法进行,因此,导数类算法往往对原始 数据的噪声有放大作用,并且导数极值的方向性和连续性较差,对弱异常的边界识别能力 较弱,导致对异常地质边界的提取精度较低。
【发明内容】
[0005] 本申请的目的在于提供一种重力异常边界增强的方法及装置,可以降低重力异常 边界信息的噪声,提高重力异常边界的分辨力。
[0006] 本发明提供了一种重力异常边界增强的方法及装置,所述方法及装置具体是这样 实现的:
[0007] 一种重力异常边界增强的方法,所述方法包括:
[0008] 获取目标工区的重力异常数据,并对所述重力异常数据进行网格化;
[0009] 采用空间域算子逐个提取第一网格点在预设方向的重力异常边界极值信息;
[0010] 基于所述第一网格点的重力异常边界极值信息,计算第二网格点的相干增强值。
[0011] 可选的,在本发明的一个实施例中,所述基于所述第一网格点的重力异常边界极 值信息,计算第二网格点的相干增强值,包括:
[0012] 基于所述第一网格点的重力异常边界极值信息,建立以第二网格点为中心、大小 为第一预设值的增强值计算窗口;
[0013] 提取与所述增强值计算窗口相邻八个方向的、大小为所述第一预设值的相邻数据 窗口;
[0014] 分别计算所述增强值计算窗口与所述相邻八个方向的相邻数据窗口的相关系数 值;
[0015] 将所述八个相关系数值的乘积的八次开方值作为所述第二网格点的相干增强值。
[0016] 可选的,在本发明的一个实施例中,在所述对所述重力资料进行网格化之后,还包 括:
[0017] 采用低通滤波法去除所述重力异常资料中的噪声。
[0018] 可选的,在本发明的一个实施例中,所述方法还包括:
[0019] 判断所述第二网格点的相干增强值是否小于预设阈值;
[0020] 当判断结果为是时,基于所述第一网格点的相干增强值信息,计算所述第二网格 点的二次相干增强值。
[0021] 可选的,在本发明的一个实施例中,所述采用空间域算子逐个提取第一网格点在 预设方向的重力异常边界极值信息包括:
[0022] 获取所述目标工区在预设方向上的空间域算子,所述空间域算子的大小为第二预 设值;
[0023]提取第一网格点中以待计算网格点为中心、大小为第二预设值的极值计算窗口;
[0024] 将所述极值计算窗口与所述空间域算子的乘积作为所述待计算网格点的重力异 常边界极值;
[0025] 遍历所述第一网格点中的其他网格点,生成所述第一网格点在预设方向的重力异 常边界极值信息。
[0026] 一种重力异常边界增强的装置,所述装置包括:
[0027]网格化单元,用于获取目标工区的重力异常数据,并对所述重力异常数据进行网 格化;
[0028] 边界极值计算单元,用于采用空间域算子逐个提取第一网格点在预设方向的重力 异常边界极值信息;
[0029] 增强值计算单元,用于基于所述第一网格点的重力异常边界极值信息,计算第二 网格点的相干增强值。
[0030] 可选的,在本发明的一个实施例中,所述增强值计算单元包括:
[0031 ]增强值窗口建立单元,用于基于所述第一网格点的重力异常边界极值信息,建立 以第二网格点为中心、大小为第一预设值的增强值计算窗口;
[0032]相邻窗口提取单元,用于提取与所述增强值计算窗口相邻八个方向的、大小为所 述第一预设值的相邻数据窗口;
[0033]相关系数计算单元,用于分别计算所述增强值计算窗口与所述相邻八个方向的相 邻数据窗口的相关系数值;
[0034]开方计算单元,用于将所述八个相关系数值的乘积的八次开方值作为所述第二网 格点的相干增强值。
[0035] 可选的,在本发明的一个实施例中,所述装置还包括:
[0036] 去噪单元,用于采用低通滤波法去除所述重力异常资料中的噪声。
[0037] 可选的,在本发明的一个实施例中,所述装置还包括:
[0038] 增强值判断单元,用于判断所述第二网格点的相干增强值是否小于预设阈值;
[0039] 二次增强单元,用于当所述增强值判断单元的判断结果为是时,基于所述第一网 格点的相干增强值信息,计算所述第二网格点的二次相干增强值。
[0040] 可选的,在本发明的一个实施例中,所述边界极值计算单元包括:
[0041] 算子获取单元,用于获取所述目标工区在预设方向上的空间域算子,所述空间域 算子的大小为第二预设值;
[0042] 极值窗口提取单元,用于提取第一网格点中以待计算网格点为中心、大小为第二 预设值的极值计算窗口;
[0043] 极值计算单元,用于将所述极值计算窗口与所述空间域算子的乘积作为所述待计 算网格点的重力异常边界极值;
[0044] 遍历计算单元,用于遍历所述第一网格点中的其他网格点,生成所述第一网格点 在预设方向的重力异常边界极值信息。
[0045]本申请提供的重力异常边界增强方法及装置,可以提取网格化后的重力异常数据 在预设方向上的重力异常边界极值信息,所述重力异常边界的极值一般对应着模型的边界 点。然后基于所述重力异常边界极值信息,可以计算网格点的相干增强值,进一步获取模型 边界点的增强值。本申请提供的重力异常边界增强方法对所述重力异常数据进行空间域处 理,避免传统频率域算法对数据噪声的放大作用,可以获取清晰、连续的重力异常边界信 息,对目标工区地层断裂数据的研究提供有力的基础资料。
【附图说明】
[0046] 为了更清楚地说明本申请实施例或现有技术中的技术方案,下面将对实施例或现 有技术描述中所需要使用的附图作简单地介绍,显而易见地,下面描述中的附图仅仅是本 申请中记载的一些实施例,对于本领域普通技术人员来讲,在不付出创造性劳动性的前提 下,还可以根据这些附图获得其他的附图。
[0047] 图1是本发明提供的重力异常边界增强方法的一种实施例的方法流程示意图;
[0048] 图2是本发明提供的计算重力异常边界极值的一种实施例的方法流程图;
[0049] 图3是本发明提供的计算网格点相干增强值的一种实施例的方法流程图;
[0050] 图4是本发明提供的增强值计算窗口与相邻八个方向的相邻数据窗口的位置关系 图;
[0051] 图5是本发明提供的我国东北目标工区A的实测布格重力异常图;
[0052] 图6是本发明提供的对工区A进行向上延拓去噪后的布格重力异常图;
[0053] 图7是本发明提供的将工区A北东、北西两个方向的重力异常边界信息进行融合后 的断裂走势图;
[0054] 图8是本发明提供的利用总水平导数法对工区A提取和增强的重力异常边界信息 图;
[0055] 图9是本发明提供的重力异常边界增强装置的一种实施例的模块结构示意图;
[0056] 图10是本发明提供的增强值计算单元的一种实施例的模块结构示意图;
[0057] 图11是本发明提供的重力异常边界增强装置的另一种实施例的模块结构示意图; [0058]图12是本发明提供的边界极值计算单元的一种实施例的模块结构示意图。
【具体实施方式】
[0059] 为了使本技术领域的人员更好地理解本申请中的技术方案,下面将结合本申请实 施例中的附图,对本申请实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施 例仅仅是本申请一部分实施例,而不是全部的实施例。基于本申请中的实施例,本领域普通 技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都应当属于本申请保护 的范围。
[0060] 下面结合附图对本申请所述的重力异常边界增强的方法进行详细的说明。图1是 本申请提供的重力异常边界增强方法的一种实施例的方法流程示意图。虽然本申请提供了 如下述实施例或附图所示的方法操作步骤,但基于常规或者无需创造性的劳动在所述方法 中可以包括更多或者更少的操作步骤。在逻辑性上不存在必要因果关系的步骤中,这些步 骤的执行顺序不限于本申请实施例提供的执行顺序。所述方法在实际中的重力异常边界的 增强过程中或者装置执行时,可以按照实施例或者附图所示的方法顺序执行或者并行执行 (例如并行处理器或者多线程处理的环境)。
[0061] 具体的本申请提供的重力异常边界增强方法的一种实施例如图1所示,所述方法 可以包括:
[0062] S1:获取目标工区的重力异常数据,并对所述重力异常数据进行网格化。
[0063]重力观测数据一般是地质勘测人员通过重力测量仪观测得到的目标工区的重力 观测数据。当观测点的重力观测数据与该观测点的正常重力值常常存在偏差时,可以确定 所述观测点发生重力异常。引起观测点的重力异常一般存在多种因素,例如所述观测点与 大地水准面存在一定高差,当观测点位置高于大地水平面时,重力值越小偏小;当观测点与 大地水准面之间存在剩余物质时,也会在观测点处产生附加重力值。
[0064]本申请实施例中的重力异常数据可以为所述目标工区的布格重力异常数据。本实 施例中,可以对所述重力观测数据中引起重力异常的影响因素进行布格校正,若在布格校 正之后仍然存在的异常称为布格重力异常。可以将校正后的重力观测数据中与对应观测点 的正常重力值的差值作为布格重力异常,差值为正称正异常,差值为负称负异常。对所述重 力观测数据进行布格处理之后,消除一些引起重力异常的客观因素带来的影响,可以大大 降低重力观测数据庞大的数据量,显著提高后续计算处理的效率。
[0065]在获取目标工区的重力异常数据之后,可以对所述重力异常数据进行网格化。具 体地,可以将布格重力异常图以一定的网格大小划分成正方形规则网格,所述正方形规则 网格的大小和形状可以根据研究目的和工作比例尺而定,网格大小一般为重力观测网格间 距的几倍到十几倍。将正方形规则网格中包含的观测节点重力异常值的平均值作为所述正 方形规则网格的代表值。对所述重力异常数据进行网格化,可以将空间上不均匀分布的重 力异常数据均匀化,一方面可以抑制局部噪音,另一方面可以弥补空白网格的数值。
[0066]本申请实施例中,可以获取目标工区的重力异常数据,并对所述重力异常数据进 行网格化,对所述重力异常数据进行网格化,在一定程度上可以抑制局部噪音。
[0067]本申请实施例中,在所述对所述重力资料进行网格化之后,还可以包括:
[0068]采用低通滤波法去除所述重力异常资料中的噪声。
[0069]所述低通滤波法可以包括滑动平均法、向上延拓法等。其中,对于重力异常的向上 延拓可以突出规模较大的异常体的异常特征,所述规模较大的异常体可以包括区域性的或 者深部较大规模的异常体,并压制规模较小的异常体的异常特征,所述规模较小的异常体 可以包括局部的或者浅而小的异常体。所述向上延拓法可以消除近地表效应,去除重力异 常数据的噪声影响。
[0070] S2:采用空间域算子逐个提取第一网格点在预设方向的重力异常边界极值信息。
[0071] 一般地,重力异常具有方向性,同一区域的重力异常方向往往是一致的,不同地区 的重力异常方向一般具有较大差异。例如我国的华南、下扬子地区、东北松辽平原地区的重 力异常呈北东向分布,内蒙古、甘肃、新疆及黄土高原的重力异常以东西向排列为主。本实 施例中,可以根据目标地区的区域特征,确定进行重力异常边界增强的预设方向。图2是本 申请提供的计算重力异常边界极值的一种实施例的方法流程图,如图2所示,所述方法包 括:
[0072] S11:获取所述目标工区在预设方向上的空间域算子,所述空间域算子的大小为第 二预设值。
[0073] 本实施例中的空间域算子可以包括方向差分空间域算子,其中,如表1-2所示,东 西方向的差分算子和南北方向的差分算子为索贝尔(Sobel)算子,索贝尔算子可以对网格 化后的重力异常数据进行加权平均滤波,具有一定的噪声抑制能力。如表3-4所示,所述北 东方向差分算子和北西方向差分算子为差分算子,同样可以对网格化后的重力异常数据有 噪声抑制作用。当然,本实施例中的空间域算子不仅仅限于方向差分空间域算子,还可以包 括Roberts算子、Canny算子、Prewitt算子、边缘检测算子、二阶算子等。一般地,空间域算子 的大小是固定的,例如索贝尔算子的大小为3X3矩阵,可以获取空间域算子的大小作为第 二预设值,用于后续的选取数据窗口大小。
[0074] 需要说明的是,根据目标工区的区域特征,一般只选取预设方向的方向差分算子 对所述网格化后的重力异常数据进行极值处理。例如,位于我国华南的目标工区,可以利用 北东方向差分算子提取所述目标工区的重力异常边界极值信息。
[0076] S12:提取第一网格点中以待计算网格点为中心、大小为第二预设值的极值计算窗 □ 〇
[0077]本实施例中的第一网格点可以为网格化后的重力异常数据中网格点的集合,所述 待计算网格点属于第一网格点中的单个网格点。提取以所述待计算网格点为中心、大小为 第二预设值的极值计算窗口。例如,提取以待计算网格点为中心的大小为3X3的极值计 算窗口,所述极值计算窗口数据为:
[0079] S13:将所述极值计算窗口与所述空间域算子的乘积作为所述待计算网格点的重 力异常边界极值。
[0080] 本申请实施例中,可以将所述极值计算窗口中的数据与所述空间域算子的乘积作 为所述待计算网格点的重力异常边界极值。例如,基于上述表1-4的方向差分空间域算子以 及以待计算网格点ai, j为中心的极值计算窗口,可以分别计算所述待计算网格点ai, j在东 西、南北、北东、北西方向的重力异常边界极值。具体地,对于东西方向,重力异常边界极值 的公式为:
[0081 ] Ai, j = ai-i, j+i-ai-i, j-i+2(ai, j+i-ai, j-i)+ai+i, j+i-ai+i, j-i ;
[0082]对于南北方向,重力异常边界极值的公式为:
[0083] Ai, j = ai+i, j+i-ai-i, j+i+2(ai+i, j-ai-i, j)+ai+i, j-i-ai-i, j-i ;
[0084]对于北东方向,重力异常边界极值的公式为:
[0085] Ai,j = ai+i,j-ai,j+i+ai,j-i-ai-ij;
[0086]对于北西方向,重力异常边界极值的公式为:
[0087] Ai, j = ai, j+i-ai-1, j+ai+i, j_ai, j-1。
[0088] S14:遍历所述第一网格点中的其他网格点,生成所述第一网格点在预设方向的重 力异常边界极值信息。
[0089] 在计算得到所述待计算网格点的重力异常边界极值后,可以按照同样的方法计算 所述第一网格点中的其他网格点,直至计算得到第一网格点中所有网格点的重力异常边界 极值,生成所述第一网格点在预设方向的重力异常边界极值信息。
[0090] 需要说明的是,在所述网格化后的重力异常数据中,处于边缘位置的网格点一般 无法生成极值计算窗口,此时,需要对所述重力异常数据进行扩边。本申请实施例中,可以 采用余弦衰减至零的方法对所述网格化后的重力异常数据进行扩边。设网格化后的重力异 常数据为2¥(1 :111,1:11),111为线数,11为每条线上的点数,设扩边数为1例如,当1 = 2时,扩边 后的数据为ZK( 1-t :m+t,1-t: n+t)。设置ZK( i,j) = ZY( i,j),( 1 彡 i彡m,1 彡 j彡n),则ZK中的 数据可以按如下公式扩边:
[0091] ZK(i, j)=ZY(i,l)Xcos((3iX(l-j))/2t),(l^i^m,l-t^j^0);
[0092] ZK(i,j)=ZY(i,n)X cos( (it X (j-n) )/2t), (Ki<m,n+K j<n+t);
[0093] ZK(i, j)=ZY(l, j) Xcos((3iX(l-i))/2t),(l-t^i^0,l-t^j^n+t);
[0094] ZK( i,j) = ZY(m, j) X cos( (it X(i-m))/2t), (m+KiSm+t,l_t< j<n+t)。
[0095] 本实施例中,可以将网格化后的重力异常数据转化为极值信息,用于后续的相干 增强处理,可以避免传统频率域算法中对网格化后的重力异常数据噪声的放大作用。
[0096] S3:基于所述第一网格点的重力异常边界极值信息,计算第二网格点的相干增强 值。
[0097] 本实施例中的第二网格点为所述第一网格点中的任意一个网格点,基于所述目标 工区中第一网格点的重力异常边界极值信息,可以计算第二网格点的相干增强值。图3是本 申请提供的计算网格点相干增强值的一种实施例的方法流程图,如图3所示,所述方法包 括:
[0098] S101:基于所述第一网格点的重力异常边界极值信息,建立以第二网格点为中心、 大小为第一预设值的增强值计算窗口。
[0099] 本实施例中,可以基于所述第一网格点的重力异常边界极值信息,建立以第二网 格点为中心、大小为第一预设值的增强值计算窗口。所述第一预设值大小可以根据整体网 格数目设置,本实施例中,所述第一预设值可以为5 X 5。
[0100] S102:提取与所述增强值计算窗口相邻八个方向的、大小为所述第一预设值的相 邻数据窗口。
[0101] 本实施例中,可以计算所述增强值计算窗口与相邻八个方向的窗口的相关系数 值。因此,可以提取与所述增强值计算窗口相邻八个方向的、大小相等的相邻数据窗口。图4 是本申请提供的增强值计算窗口与相邻八个方向的相邻数据窗口的位置关系图,如图4所 示,其中,西北相邻数据窗口是将所述增强值计算窗口的位置向上然后向左平移一个网格 得到的窗口位置。
[0102] S103:分别计算所述增强值计算窗口与所述相邻八个方向的相邻数据窗口的相关 系数值。
[0103] 本实施例中,可以分别计算所述增强值计算窗口与所述相邻八个方向的相邻数据 窗口的相关系数值。其中,所述相关系数值的计算公式可以为:
[0105] 其中,X可以为增强值计算窗口数据,Y可以为任一方位相邻数据窗口内的数据 X、歹分别为对应窗口全部数据的均值,n为第一预设值,k为1~8。
[0106] 例如,基于所述第一网格点在预设方向的重力异常边界极值信息数据,X和Y可以 为下式所示数据,根据上述相关系数值的计算公式,可以计算得到第二网格点乂 13与八个相 邻数据窗口的相关系数值。
[0108] S104:将所述八个相关系数值的乘积的八次开方值作为所述第二网格点的相干增 强值。
[0109]假设计算得到的8个相关系数值为^,所述第二网格点的相干增强值的计算公式 可以为:
[0111] 后续的,还可以对相干增强后的重力异常数据进行二次相干增强,所述二次相干 增强可以包括下述步骤:
[0112] SS1:判断所述第二网格点的相干增强值是否小于预设阈值;
[0113] SS2:当判断结果为是时,基于所述第一网格点的相干增强值信息,计算所述第二 网格点的二次相干增强值。
[0114] 当判断所述第二网格点的相干增强值小于预设阈值时,可以基于所述第一网格点 的相干增强值信息,按照上述S101-S104的方法步骤计算所述第二网格点的二次相干增强 值。
[0115] 下面通过一个具体的应用场景说明上述实施例方法,图5是本发明提供的我国东 北目标工区A的实测布格重力异常图,如图5所示,图中颜色的深浅表示重力异常值的高低, 颜色越深,表示重力异常值越高,反之,颜色越浅,表示重力异常值越低。利用向上延拓2千 米的方法去除图5中布格重力异常图的噪声,得到图6所示的对工区A进行向上延拓去噪后 的布格重力异常图。根据工区A已知的地质情况,采用北东、北西方向差分算子提取北东、北 西两个方向的重力异常边界极值信息。基于提取的北东、北西两个方向重力异常边界极值 信息,对图6中的各个网格点进行相干增强处理,得到工区A的第一次相干增强数据。然后采 用阈值0.985对所述第一次相干增强数据进行二次相干增强处理,得到该区北东向和北西 向的重力异常边界信息,图7是将工区A北东、北西两个方向重力异常边界信息融合后的断 裂走势图,反映了该地区断陷断裂的主体分布情况。图8是本发明提供的利用总水平导数法 对工区A提取和增强的重力异常边界信息图,将图7与图8比较,可以发现本发明所提取和增 强的重力异常边界信息更加清晰、连续。另外,如图7所示本发明提供的实施例方法不仅可 以提取和增强工区A北东向的主要断裂,而且还对正交方向(北西向)的扭曲和错断断裂反 映的比较清楚,表现了比传统的倒数类算法具有更强的处理能力。
[0116]本申请提供的重力异常边界增强方法,可以提取网格化后的重力异常数据在预设 方向上的重力异常边界极值信息,所述重力异常边界的极值一般对应着模型的边界点。然 后基于所述重力异常边界极值信息,可以计算网格点的相干增强值,进一步获取模型边界 点的增强值。本申请提供的重力异常边界增强方法对所述重力异常数据进行空间域处理, 避免传统频率域算法对数据噪声的放大作用,可以获取清晰、连续的重力异常边界信息,对 目标工区地层断裂数据的研究提供有力的基础资料。
[0117]本发明另一方面还提供一种重力异常边界增强的装置,图9是本申请提供的重力 异常边界增强装置的一种实施例的模块结构示意图,如图9所示,所述装置90可以包括: [0118]网格化单元91,用于获取目标工区的重力异常数据,并对所述重力异常数据进行 网格化;
[0119] 边界极值计算单元92,用于采用空间域算子逐个提取第一网格点在预设方向的重 力异常边界极值信息;
[0120] 增强值计算单元93,用于基于所述第一网格点的重力异常边界极值信息,计算第 二网格点的相干增强值。
[0121] 本申请提供的重力异常边界增强装置,可以提取网格化后的重力异常数据在预设 方向上的重力异常边界极值信息,所述重力异常边界的极值一般对应着模型的边界点。然 后基于所述重力异常边界极值信息,可以计算网格点的相干增强值,进一步获取模型边界 点的增强值。本申请提供的重力异常边界增强方法对所述重力异常数据进行空间域处理, 避免传统频率域算法对数据噪声的放大作用,可以获取清晰、连续的重力异常边界信息,对 目标工区地层断裂数据的研究提供有力的基础资料。
[0122] 图10是本申请提供的增强值计算单元的一种实施例的模块结构示意图,如图10所 示,所述增强值计算单元93可以包括:
[0123] 增强值窗口建立单元101,用于基于所述第一网格点的重力异常边界极值信息,建 立以第二网格点为中心、大小为第一预设值的增强值计算窗口;
[0124] 相邻窗口提取单元102,用于提取与所述增强值计算窗口相邻八个方向的、大小为 所述第一预设值的相邻数据窗口;
[0125] 相关系数计算单元103,用于分别计算所述增强值计算窗口与所述相邻八个方向 的相邻数据窗口的相关系数值;
[0126] 开方计算单元104,用于将所述八个相关系数值的乘积的八次开方值作为所述第 二网格点的相干增强值。
[0127] 可选的,在本发明的一个实施例中,所述装置还包括:
[0128] 去噪单元,用于采用低通滤波法去除所述重力异常资料中的噪声。
[0129] 图11是本申请提供的重力异常边界增强装置的另一种实施例的模块结构示意图, 如图11所示,所述装置110还可以包括:
[0130] 增强值判断单元111,用于判断所述第二网格点的相干增强值是否小于预设阈值;
[0131] 二次增强单元112,用于当所述增强值判断单元的判断结果为是时,基于所述第一 网格点的相干增强值信息,计算所述第二网格点的二次相干增强值。
[0132] 图12是本申请提供的边界极值计算单元的一种实施例的模块结构示意图,如图12 所示,所述边界极值计算单元92可以包括:
[0133] 算子获取单元121,用于获取所述目标工区在预设方向上的空间域算子,所述空间 域算子的大小为第二预设值;
[0134] 极值窗口提取单元122,用于提取第一网格点中以待计算网格点为中心、大小为第 二预设值的极值计算窗口;
[0135] 极值计算单元123,用于将所述极值计算窗口与所述空间域算子的乘积作为所述 待计算网格点的重力异常边界极值;
[0136] 遍历计算单元124,用于遍历所述第一网格点中的其他网格点,生成所述第一网格 点在预设方向的重力异常边界极值信息。
[0137] 本说明书中的各个实施例均采用递进的方式描述,每个实施例重点说明的都是与 其他实施例的不同之处,各个实施例之间相同相似的部分互相参见即可。尤其,对于系统实 施例而言,由于其基本相似于方法实施例,所以描述的比较简单,相关之处参见方法实施例 的部分说明即可。
[0138] 虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和 变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的 精神。
[0139] 尽管本申请内容中提到图片裁剪、RGB通道颜色转换、微调(fine-tuning)、卷积神 经网络模型层结构、卷积、损失反馈等的图片信息处理、神经网络模型结构的描述,但是,本 申请并不局限于必须是完全标准或者所提及的方式的信息处理、神经网络模型结构的情 况。本申请中各个实施例所涉及的上述描述仅是本申请中的一些实施例中的应用,在某些 标准、模型、方法的基础上略加修改后的实施方式也可以实行上述本申请各实施例的方案。 当然,在符合本申请上述各实施例的中所述的处理方法步骤的其他无创造性的变形,仍然 可以实现相同的申请,在此不再赘述。
[0140] 虽然本申请提供了如实施例或流程图所述的方法操作步骤,但基于常规或者无创 造性的劳动可以包括更多或者更少的操作步骤。实施例中列举的步骤顺序仅仅为众多步骤 执行顺序中的一种方式,不代表唯一的执行顺序。在实际中的装置或客户端产品执行时,可 以按照实施例或者附图所示的方法顺序执行或者并行执行(例如并行处理器或者多线程处 理的环境)。
[0141] 上述实施例阐明的装置或模块,具体可以由计算机芯片或实体实现,或者由具有 某种功能的产品来实现。为了描述的方便,描述以上装置时以功能分为各种模块分别描述。 在实施本申请时可以把各模块的功能在同一个或多个软件和/或硬件中实现。当然,也可以 将实现某功能的模块由多个子模块或子单元组合实现。
[0142] 本申请中所述的方法、装置或模块可以以计算机可读程序代码方式实现控制器按 任何适当的方式实现,例如,控制器可以采取例如微处理器或处理器以及存储可由该(微) 处理器执行的计算机可读程序代码(例如软件或固件)的计算机可读介质、逻辑门、开关、专 用集成电路(Application Specific Integrated Circuit,ASIC)、可编程逻辑控制器和嵌 入微控制器的形式,控制器的例子包括但不限于以下微控制器:ARC 62f5D、Atmel AT91SAM、 Microchip PIC18F26K20以及Silicone Labs C8051F320,存储器控制器还可以被实现为存 储器的控制逻辑的一部分。本领域技术人员也知道,除了以纯计算机可读程序代码方式实 现控制器以外,完全可以通过将方法步骤进行逻辑编程来使得控制器以逻辑门、开关、专用 集成电路、可编程逻辑控制器和嵌入微控制器等的形式来实现相同功能。因此这种控制器 可以被认为是一种硬件部件,而对其内部包括的用于实现各种功能的装置也可以视为硬件 部件内的结构。或者甚至,可以将用于实现各种功能的装置视为既可以是实现方法的软件 模块又可以是硬件部件内的结构。
[0143] 本申请所述装置中的部分模块可以在由计算机执行的计算机可执行指令的一般 上下文中描述,例如程序模块。一般地,程序模块包括执行特定任务或实现特定抽象数据类 型的例程、程序、对象、组件、数据结构、类等等。也可以在分布式计算环境中实践本申请,在 这些分布式计算环境中,由通过通信网络而被连接的远程处理设备来执行任务。在分布式 计算环境中,程序模块可以位于包括存储设备在内的本地和远程计算机存储介质中。
[0144] 通过以上的实施方式的描述可知,本领域的技术人员可以清楚地了解到本申请可 借助软件加必需的硬件的方式来实现。基于这样的理解,本申请的技术方案本质上或者说 对现有技术做出贡献的部分可以以软件产品的形式体现出来,也可以通过数据迀移的实施 过程中体现出来。该计算机软件产品可以存储在存储介质中,如R0M/RAM、磁碟、光盘等,包 括若干指令用以使得一台计算机设备(可以是个人计算机,移动终端,服务器,或者网络设 备等)执行本申请各个实施例或者实施例的某些部分所述的方法。
[0145] 本说明书中的各个实施例采用递进的方式描述,各个实施例之间相同或相似的部 分互相参见即可,每个实施例重点说明的都是与其他实施例的不同之处。本申请的全部或 者部分可用于众多通用或专用的计算机系统环境或配置中。例如:个人计算机、服务器计算 机、手持设备或便携式设备、平板型设备、移动通信终端、多处理器系统、基于微处理器的系 统、可编程的电子设备、网络PC、小型计算机、大型计算机、包括以上任何系统或设备的分布 式计算环境等等。
[0146] 虽然通过实施例描绘了本申请,本领域普通技术人员知道,本申请有许多变形和 变化而不脱离本申请的精神,希望所附的权利要求包括这些变形和变化而不脱离本申请的 精神。
【主权项】
1. 一种重力异常边界增强的方法,其特征在于,所述方法包括: 获取目标工区的重力异常数据,并对所述重力异常数据进行网格化; 采用空间域算子逐个提取第一网格点在预设方向的重力异常边界极值信息; 基于所述第一网格点的重力异常边界极值信息,计算第二网格点的相干增强值。2. 根据权利要求1所述的一种重力异常边界增强的方法,其特征在于,所述基于所述第 一网格点的重力异常边界极值信息,计算第二网格点的相干增强值,包括: 基于所述第一网格点的重力异常边界极值信息,建立以第二网格点为中心、大小为第 一预设值的增强值计算窗口; 提取与所述增强值计算窗口相邻八个方向的、大小为所述第一预设值的相邻数据窗 P; 分别计算所述增强值计算窗口与所述相邻八个方向的相邻数据窗口的相关系数值; 将所述八个相关系数值的乘积的八次开方值作为所述第二网格点的相干增强值。3. 根据权利要求1所述的一种重力异常边界增强的方法,其特征在于,在所述对所述重 力资料进行网格化之后,还包括: 采用低通滤波法去除所述重力异常资料中的噪声。4. 根据权利要求1所述的一种重力异常边界增强的方法,其特征在于,所述方法还包 括: 判断所述第二网格点的相干增强值是否小于预设阈值; 当判断结果为是时,基于所述第一网格点的相干增强值信息,计算所述第二网格点的 二次相干增强值。5. 根据权利要求1所述的一种重力异常边界增强的方法,其特征在于,所述采用空间域 算子逐个提取第一网格点在预设方向的重力异常边界极值信息包括: 获取所述目标工区在预设方向上的空间域算子,所述空间域算子的大小为第二预设 值; 提取第一网格点中以待计算网格点为中心、大小为第二预设值的极值计算窗口; 将所述极值计算窗口与所述空间域算子的乘积作为所述待计算网格点的重力异常边 界极值; 遍历所述第一网格点中的其他网格点,生成所述第一网格点在预设方向的重力异常边 界极值信息。6. -种重力异常边界增强的装置,其特征在于,所述装置包括: 网格化单元,用于获取目标工区的重力异常数据,并对所述重力异常数据进行网格化; 边界极值计算单元,用于采用空间域算子逐个提取第一网格点在预设方向的重力异常 边界极值信息; 增强值计算单元,用于基于所述第一网格点的重力异常边界极值信息,计算第二网格 点的相干增强值。7. 根据权利要求6所述的一种重力异常边界增强的装置,其特征在于,所述增强值计算 单元包括: 增强值窗口建立单元,用于基于所述第一网格点的重力异常边界极值信息,建立以第 二网格点为中心、大小为第一预设值的增强值计算窗口; 相邻窗口提取单元,用于提取与所述增强值计算窗口相邻八个方向的、大小为所述第 一预设值的相邻数据窗口; 相关系数计算单元,用于分别计算所述增强值计算窗口与所述相邻八个方向的相邻数 据窗口的相关系数值; 开方计算单元,用于将所述八个相关系数值的乘积的八次开方值作为所述第二网格点 的相干增强值。8. 根据权利要求6所述的一种重力异常边界增强的装置,其特征在于,所述装置还包 括: 去噪单元,用于采用低通滤波法去除所述重力异常资料中的噪声。9. 根据权利要求6所述的一种重力异常边界增强的装置,其特征在于,所述装置还包 括: 增强值判断单元,用于判断所述第二网格点的相干增强值是否小于预设阈值; 二次增强单元,用于当所述增强值判断单元的判断结果为是时,基于所述第一网格点 的相干增强值信息,计算所述第二网格点的二次相干增强值。10. 根据权利要6所述的一种重力异常边界增强的装置,其特征在于,所述边界极值计 算单元包括: 算子获取单元,用于获取所述目标工区在预设方向上的空间域算子,所述空间域算子 的大小为第二预设值; 极值窗口提取单元,用于提取第一网格点中以待计算网格点为中心、大小为第二预设 值的极值计算窗口; 极值计算单元,用于将所述极值计算窗口与所述空间域算子的乘积作为所述待计算网 格点的重力异常边界极值; 遍历计算单元,用于遍历所述第一网格点中的其他网格点,生成所述第一网格点在预 设方向的重力异常边界极值信息。
【文档编号】G01V1/36GK105891889SQ201610197300
【公开日】2016年8月24日
【申请日】2016年3月31日
【发明人】赵文举, 刘云祥, 赵荔
【申请人】中国石油天然气集团公司, 中国石油集团东方地球物理勘探有限责任公司