一種重力異常邊界增強(qiáng)的方法及裝置的制造方法
【專利摘要】本申請(qǐng)實(shí)施例提供的一種重力異常邊界增強(qiáng)的方法及裝置。所述方法包括：獲取目標(biāo)工區(qū)的重力異常數(shù)據(jù)，并對(duì)所述重力異常數(shù)據(jù)進(jìn)行網(wǎng)格化；采用空間域算子逐個(gè)提取第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重力異常邊界極值信息；基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，計(jì)算第二網(wǎng)格點(diǎn)的相干增強(qiáng)值。本申請(qǐng)?zhí)峁┑闹亓Ξ惓＿吔缭鰪?qiáng)方法對(duì)所述重力異常數(shù)據(jù)進(jìn)行空間域處理，避免傳統(tǒng)頻率域算法對(duì)數(shù)據(jù)噪聲的放大作用，可以獲取清晰、連續(xù)的重力異常邊界信息。
【專利說(shuō)明】
一種重力異常邊界増強(qiáng)的方法及裝置
技術(shù)領(lǐng)域
[0001] 本發(fā)明涉及石油和固體礦產(chǎn)勘探中的地球物理勘探技術(shù)領(lǐng)域，尤其涉及一種重力 異常邊界增強(qiáng)的方法及裝置。
【背景技術(shù)】
[0002] 在地球物理勘探技術(shù)領(lǐng)域，通過(guò)重力勘探可以了解地質(zhì)場(chǎng)源的分布，其中包括場(chǎng) 源的形態(tài)和規(guī)模等。重力異常是地下密度體不均勻分布的綜合反應(yīng)，利用重力異常資料中 的異常邊線性梯級(jí)帶、異常的扭曲、異常的分區(qū)特征等信息數(shù)據(jù)可以識(shí)別異常地質(zhì)邊界，比 如地層斷裂位置等。而與地質(zhì)邊界有關(guān)的信息通常表現(xiàn)為高頻率、低幅值，往往很難直接從 布格重力異常圖中提取出地質(zhì)邊界的水平位置。通常需要先對(duì)重力異常數(shù)據(jù)的邊界信息進(jìn) 行增強(qiáng)后再提取異常地質(zhì)邊界信息。
[0003] 現(xiàn)有技術(shù)中的異常邊界增強(qiáng)的方法主要有導(dǎo)數(shù)類算法，所述導(dǎo)數(shù)類算法主要包括 垂向?qū)?shù)法、總水平導(dǎo)數(shù)法、總梯度模法、斜導(dǎo)數(shù)法、西塔(Theta)圖法等。其中，所述總水平 導(dǎo)數(shù)法可以采用對(duì)橫軸和豎軸方向水平導(dǎo)數(shù)的平方和進(jìn)行開(kāi)方的方法求得。導(dǎo)數(shù)類算法的 共同特點(diǎn)在于利用垂向?qū)?shù)或?qū)?shù)的組合對(duì)信號(hào)中的高頻成分進(jìn)行放大，使重力異常圖中 的極大值或拐點(diǎn)近似位于地質(zhì)體邊界的正上方，后續(xù)可以通過(guò)邊緣檢測(cè)的方法提取出異常 地質(zhì)邊界位置。
[0004]導(dǎo)數(shù)類算法一般采用頻率域傅氏變換的方法進(jìn)行，因此，導(dǎo)數(shù)類算法往往對(duì)原始 數(shù)據(jù)的噪聲有放大作用，并且導(dǎo)數(shù)極值的方向性和連續(xù)性較差，對(duì)弱異常的邊界識(shí)別能力 較弱，導(dǎo)致對(duì)異常地質(zhì)邊界的提取精度較低。

【發(fā)明內(nèi)容】

[0005] 本申請(qǐng)的目的在于提供一種重力異常邊界增強(qiáng)的方法及裝置，可以降低重力異常 邊界信息的噪聲，提高重力異常邊界的分辨力。
[0006] 本發(fā)明提供了一種重力異常邊界增強(qiáng)的方法及裝置，所述方法及裝置具體是這樣 實(shí)現(xiàn)的：
[0007] 一種重力異常邊界增強(qiáng)的方法，所述方法包括：
[0008] 獲取目標(biāo)工區(qū)的重力異常數(shù)據(jù)，并對(duì)所述重力異常數(shù)據(jù)進(jìn)行網(wǎng)格化；
[0009] 采用空間域算子逐個(gè)提取第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重力異常邊界極值信息；
[0010] 基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，計(jì)算第二網(wǎng)格點(diǎn)的相干增強(qiáng)值。
[0011] 可選的，在本發(fā)明的一個(gè)實(shí)施例中，所述基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極 值信息，計(jì)算第二網(wǎng)格點(diǎn)的相干增強(qiáng)值，包括：
[0012] 基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，建立以第二網(wǎng)格點(diǎn)為中心、大小 為第一預(yù)設(shè)值的增強(qiáng)值計(jì)算窗口；
[0013] 提取與所述增強(qiáng)值計(jì)算窗口相鄰八個(gè)方向的、大小為所述第一預(yù)設(shè)值的相鄰數(shù)據(jù) 窗口；
[0014] 分別計(jì)算所述增強(qiáng)值計(jì)算窗口與所述相鄰八個(gè)方向的相鄰數(shù)據(jù)窗口的相關(guān)系數(shù) 值；
[0015] 將所述八個(gè)相關(guān)系數(shù)值的乘積的八次開(kāi)方值作為所述第二網(wǎng)格點(diǎn)的相干增強(qiáng)值。
[0016] 可選的，在本發(fā)明的一個(gè)實(shí)施例中，在所述對(duì)所述重力資料進(jìn)行網(wǎng)格化之后，還包 括：
[0017] 采用低通濾波法去除所述重力異常資料中的噪聲。
[0018] 可選的，在本發(fā)明的一個(gè)實(shí)施例中，所述方法還包括：
[0019] 判斷所述第二網(wǎng)格點(diǎn)的相干增強(qiáng)值是否小于預(yù)設(shè)閾值；
[0020] 當(dāng)判斷結(jié)果為是時(shí)，基于所述第一網(wǎng)格點(diǎn)的相干增強(qiáng)值信息，計(jì)算所述第二網(wǎng)格 點(diǎn)的二次相干增強(qiáng)值。
[0021] 可選的，在本發(fā)明的一個(gè)實(shí)施例中，所述采用空間域算子逐個(gè)提取第一網(wǎng)格點(diǎn)在 預(yù)設(shè)方向的重力異常邊界極值信息包括：
[0022] 獲取所述目標(biāo)工區(qū)在預(yù)設(shè)方向上的空間域算子，所述空間域算子的大小為第二預(yù) 設(shè)值；
[0023]提取第一網(wǎng)格點(diǎn)中以待計(jì)算網(wǎng)格點(diǎn)為中心、大小為第二預(yù)設(shè)值的極值計(jì)算窗口；
[0024] 將所述極值計(jì)算窗口與所述空間域算子的乘積作為所述待計(jì)算網(wǎng)格點(diǎn)的重力異 常邊界極值；
[0025] 遍歷所述第一網(wǎng)格點(diǎn)中的其他網(wǎng)格點(diǎn)，生成所述第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重力異 常邊界極值信息。
[0026] 一種重力異常邊界增強(qiáng)的裝置，所述裝置包括：
[0027]網(wǎng)格化單元，用于獲取目標(biāo)工區(qū)的重力異常數(shù)據(jù)，并對(duì)所述重力異常數(shù)據(jù)進(jìn)行網(wǎng) 格化；
[0028] 邊界極值計(jì)算單元，用于采用空間域算子逐個(gè)提取第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重力 異常邊界極值信息；
[0029] 增強(qiáng)值計(jì)算單元，用于基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，計(jì)算第二 網(wǎng)格點(diǎn)的相干增強(qiáng)值。
[0030] 可選的，在本發(fā)明的一個(gè)實(shí)施例中，所述增強(qiáng)值計(jì)算單元包括：
[0031 ]增強(qiáng)值窗口建立單元，用于基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，建立 以第二網(wǎng)格點(diǎn)為中心、大小為第一預(yù)設(shè)值的增強(qiáng)值計(jì)算窗口；
[0032]相鄰窗口提取單元，用于提取與所述增強(qiáng)值計(jì)算窗口相鄰八個(gè)方向的、大小為所 述第一預(yù)設(shè)值的相鄰數(shù)據(jù)窗口；
[0033]相關(guān)系數(shù)計(jì)算單元，用于分別計(jì)算所述增強(qiáng)值計(jì)算窗口與所述相鄰八個(gè)方向的相 鄰數(shù)據(jù)窗口的相關(guān)系數(shù)值；
[0034]開(kāi)方計(jì)算單元，用于將所述八個(gè)相關(guān)系數(shù)值的乘積的八次開(kāi)方值作為所述第二網(wǎng) 格點(diǎn)的相干增強(qiáng)值。
[0035] 可選的，在本發(fā)明的一個(gè)實(shí)施例中，所述裝置還包括：
[0036] 去噪單元，用于采用低通濾波法去除所述重力異常資料中的噪聲。
[0037] 可選的，在本發(fā)明的一個(gè)實(shí)施例中，所述裝置還包括：
[0038] 增強(qiáng)值判斷單元，用于判斷所述第二網(wǎng)格點(diǎn)的相干增強(qiáng)值是否小于預(yù)設(shè)閾值；
[0039] 二次增強(qiáng)單元，用于當(dāng)所述增強(qiáng)值判斷單元的判斷結(jié)果為是時(shí)，基于所述第一網(wǎng) 格點(diǎn)的相干增強(qiáng)值信息，計(jì)算所述第二網(wǎng)格點(diǎn)的二次相干增強(qiáng)值。
[0040] 可選的，在本發(fā)明的一個(gè)實(shí)施例中，所述邊界極值計(jì)算單元包括：
[0041] 算子獲取單元，用于獲取所述目標(biāo)工區(qū)在預(yù)設(shè)方向上的空間域算子，所述空間域 算子的大小為第二預(yù)設(shè)值；
[0042] 極值窗口提取單元，用于提取第一網(wǎng)格點(diǎn)中以待計(jì)算網(wǎng)格點(diǎn)為中心、大小為第二 預(yù)設(shè)值的極值計(jì)算窗口；
[0043] 極值計(jì)算單元，用于將所述極值計(jì)算窗口與所述空間域算子的乘積作為所述待計(jì) 算網(wǎng)格點(diǎn)的重力異常邊界極值；
[0044] 遍歷計(jì)算單元，用于遍歷所述第一網(wǎng)格點(diǎn)中的其他網(wǎng)格點(diǎn)，生成所述第一網(wǎng)格點(diǎn) 在預(yù)設(shè)方向的重力異常邊界極值信息。
[0045]本申請(qǐng)?zhí)峁┑闹亓Ξ惓＿吔缭鰪?qiáng)方法及裝置，可以提取網(wǎng)格化后的重力異常數(shù)據(jù) 在預(yù)設(shè)方向上的重力異常邊界極值信息，所述重力異常邊界的極值一般對(duì)應(yīng)著模型的邊界 點(diǎn)。然后基于所述重力異常邊界極值信息，可以計(jì)算網(wǎng)格點(diǎn)的相干增強(qiáng)值，進(jìn)一步獲取模型 邊界點(diǎn)的增強(qiáng)值。本申請(qǐng)?zhí)峁┑闹亓Ξ惓＿吔缭鰪?qiáng)方法對(duì)所述重力異常數(shù)據(jù)進(jìn)行空間域處 理，避免傳統(tǒng)頻率域算法對(duì)數(shù)據(jù)噪聲的放大作用，可以獲取清晰、連續(xù)的重力異常邊界信 息，對(duì)目標(biāo)工區(qū)地層斷裂數(shù)據(jù)的研究提供有力的基礎(chǔ)資料。
【附圖說(shuō)明】
[0046] 為了更清楚地說(shuō)明本申請(qǐng)實(shí)施例或現(xiàn)有技術(shù)中的技術(shù)方案，下面將對(duì)實(shí)施例或現(xiàn) 有技術(shù)描述中所需要使用的附圖作簡(jiǎn)單地介紹，顯而易見(jiàn)地，下面描述中的附圖僅僅是本 申請(qǐng)中記載的一些實(shí)施例，對(duì)于本領(lǐng)域普通技術(shù)人員來(lái)講，在不付出創(chuàng)造性勞動(dòng)性的前提 下，還可以根據(jù)這些附圖獲得其他的附圖。
[0047] 圖1是本發(fā)明提供的重力異常邊界增強(qiáng)方法的一種實(shí)施例的方法流程示意圖；
[0048] 圖2是本發(fā)明提供的計(jì)算重力異常邊界極值的一種實(shí)施例的方法流程圖；
[0049] 圖3是本發(fā)明提供的計(jì)算網(wǎng)格點(diǎn)相干增強(qiáng)值的一種實(shí)施例的方法流程圖；
[0050] 圖4是本發(fā)明提供的增強(qiáng)值計(jì)算窗口與相鄰八個(gè)方向的相鄰數(shù)據(jù)窗口的位置關(guān)系 圖；
[0051] 圖5是本發(fā)明提供的我國(guó)東北目標(biāo)工區(qū)A的實(shí)測(cè)布格重力異常圖；
[0052] 圖6是本發(fā)明提供的對(duì)工區(qū)A進(jìn)行向上延拓去噪后的布格重力異常圖；
[0053] 圖7是本發(fā)明提供的將工區(qū)A北東、北西兩個(gè)方向的重力異常邊界信息進(jìn)行融合后 的斷裂走勢(shì)圖；
[0054] 圖8是本發(fā)明提供的利用總水平導(dǎo)數(shù)法對(duì)工區(qū)A提取和增強(qiáng)的重力異常邊界信息 圖；
[0055] 圖9是本發(fā)明提供的重力異常邊界增強(qiáng)裝置的一種實(shí)施例的模塊結(jié)構(gòu)示意圖；
[0056] 圖10是本發(fā)明提供的增強(qiáng)值計(jì)算單元的一種實(shí)施例的模塊結(jié)構(gòu)示意圖；
[0057] 圖11是本發(fā)明提供的重力異常邊界增強(qiáng)裝置的另一種實(shí)施例的模塊結(jié)構(gòu)示意圖； [0058]圖12是本發(fā)明提供的邊界極值計(jì)算單元的一種實(shí)施例的模塊結(jié)構(gòu)示意圖。
【具體實(shí)施方式】
[0059] 為了使本技術(shù)領(lǐng)域的人員更好地理解本申請(qǐng)中的技術(shù)方案，下面將結(jié)合本申請(qǐng)實(shí) 施例中的附圖，對(duì)本申請(qǐng)實(shí)施例中的技術(shù)方案進(jìn)行清楚、完整地描述，顯然，所描述的實(shí)施 例僅僅是本申請(qǐng)一部分實(shí)施例，而不是全部的實(shí)施例?；诒旧暾?qǐng)中的實(shí)施例，本領(lǐng)域普通 技術(shù)人員在沒(méi)有做出創(chuàng)造性勞動(dòng)前提下所獲得的所有其他實(shí)施例，都應(yīng)當(dāng)屬于本申請(qǐng)保護(hù) 的范圍。
[0060] 下面結(jié)合附圖對(duì)本申請(qǐng)所述的重力異常邊界增強(qiáng)的方法進(jìn)行詳細(xì)的說(shuō)明。圖1是 本申請(qǐng)?zhí)峁┑闹亓Ξ惓＿吔缭鰪?qiáng)方法的一種實(shí)施例的方法流程示意圖。雖然本申請(qǐng)?zhí)峁┝?如下述實(shí)施例或附圖所示的方法操作步驟，但基于常規(guī)或者無(wú)需創(chuàng)造性的勞動(dòng)在所述方法 中可以包括更多或者更少的操作步驟。在邏輯性上不存在必要因果關(guān)系的步驟中，這些步 驟的執(zhí)行順序不限于本申請(qǐng)實(shí)施例提供的執(zhí)行順序。所述方法在實(shí)際中的重力異常邊界的 增強(qiáng)過(guò)程中或者裝置執(zhí)行時(shí)，可以按照實(shí)施例或者附圖所示的方法順序執(zhí)行或者并行執(zhí)行 (例如并行處理器或者多線程處理的環(huán)境）。
[0061] 具體的本申請(qǐng)?zhí)峁┑闹亓Ξ惓＿吔缭鰪?qiáng)方法的一種實(shí)施例如圖1所示，所述方法 可以包括：
[0062] S1:獲取目標(biāo)工區(qū)的重力異常數(shù)據(jù)，并對(duì)所述重力異常數(shù)據(jù)進(jìn)行網(wǎng)格化。
[0063]重力觀測(cè)數(shù)據(jù)一般是地質(zhì)勘測(cè)人員通過(guò)重力測(cè)量?jī)x觀測(cè)得到的目標(biāo)工區(qū)的重力 觀測(cè)數(shù)據(jù)。當(dāng)觀測(cè)點(diǎn)的重力觀測(cè)數(shù)據(jù)與該觀測(cè)點(diǎn)的正常重力值常常存在偏差時(shí)，可以確定 所述觀測(cè)點(diǎn)發(fā)生重力異常。引起觀測(cè)點(diǎn)的重力異常一般存在多種因素，例如所述觀測(cè)點(diǎn)與 大地水準(zhǔn)面存在一定高差，當(dāng)觀測(cè)點(diǎn)位置高于大地水平面時(shí)，重力值越小偏??；當(dāng)觀測(cè)點(diǎn)與 大地水準(zhǔn)面之間存在剩余物質(zhì)時(shí)，也會(huì)在觀測(cè)點(diǎn)處產(chǎn)生附加重力值。
[0064]本申請(qǐng)實(shí)施例中的重力異常數(shù)據(jù)可以為所述目標(biāo)工區(qū)的布格重力異常數(shù)據(jù)。本實(shí) 施例中，可以對(duì)所述重力觀測(cè)數(shù)據(jù)中引起重力異常的影響因素進(jìn)行布格校正，若在布格校 正之后仍然存在的異常稱為布格重力異常?？梢詫⑿Ｕ蟮闹亓τ^測(cè)數(shù)據(jù)中與對(duì)應(yīng)觀測(cè)點(diǎn) 的正常重力值的差值作為布格重力異常，差值為正稱正異常，差值為負(fù)稱負(fù)異常。對(duì)所述重 力觀測(cè)數(shù)據(jù)進(jìn)行布格處理之后，消除一些引起重力異常的客觀因素帶來(lái)的影響，可以大大 降低重力觀測(cè)數(shù)據(jù)龐大的數(shù)據(jù)量，顯著提高后續(xù)計(jì)算處理的效率。
[0065]在獲取目標(biāo)工區(qū)的重力異常數(shù)據(jù)之后，可以對(duì)所述重力異常數(shù)據(jù)進(jìn)行網(wǎng)格化。具 體地，可以將布格重力異常圖以一定的網(wǎng)格大小劃分成正方形規(guī)則網(wǎng)格，所述正方形規(guī)則 網(wǎng)格的大小和形狀可以根據(jù)研究目的和工作比例尺而定，網(wǎng)格大小一般為重力觀測(cè)網(wǎng)格間 距的幾倍到十幾倍。將正方形規(guī)則網(wǎng)格中包含的觀測(cè)節(jié)點(diǎn)重力異常值的平均值作為所述正 方形規(guī)則網(wǎng)格的代表值。對(duì)所述重力異常數(shù)據(jù)進(jìn)行網(wǎng)格化，可以將空間上不均勻分布的重 力異常數(shù)據(jù)均勻化，一方面可以抑制局部噪音，另一方面可以彌補(bǔ)空白網(wǎng)格的數(shù)值。
[0066]本申請(qǐng)實(shí)施例中，可以獲取目標(biāo)工區(qū)的重力異常數(shù)據(jù)，并對(duì)所述重力異常數(shù)據(jù)進(jìn) 行網(wǎng)格化，對(duì)所述重力異常數(shù)據(jù)進(jìn)行網(wǎng)格化，在一定程度上可以抑制局部噪音。
[0067]本申請(qǐng)實(shí)施例中，在所述對(duì)所述重力資料進(jìn)行網(wǎng)格化之后，還可以包括：
[0068]采用低通濾波法去除所述重力異常資料中的噪聲。
[0069]所述低通濾波法可以包括滑動(dòng)平均法、向上延拓法等。其中，對(duì)于重力異常的向上 延拓可以突出規(guī)模較大的異常體的異常特征，所述規(guī)模較大的異常體可以包括區(qū)域性的或 者深部較大規(guī)模的異常體，并壓制規(guī)模較小的異常體的異常特征，所述規(guī)模較小的異常體 可以包括局部的或者淺而小的異常體。所述向上延拓法可以消除近地表效應(yīng)，去除重力異 常數(shù)據(jù)的噪聲影響。
[0070] S2:采用空間域算子逐個(gè)提取第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重力異常邊界極值信息。
[0071] 一般地，重力異常具有方向性，同一區(qū)域的重力異常方向往往是一致的，不同地區(qū) 的重力異常方向一般具有較大差異。例如我國(guó)的華南、下?lián)P子地區(qū)、東北松遼平原地區(qū)的重 力異常呈北東向分布，內(nèi)蒙古、甘肅、新疆及黃土高原的重力異常以東西向排列為主。本實(shí) 施例中，可以根據(jù)目標(biāo)地區(qū)的區(qū)域特征，確定進(jìn)行重力異常邊界增強(qiáng)的預(yù)設(shè)方向。圖2是本 申請(qǐng)?zhí)峁┑挠?jì)算重力異常邊界極值的一種實(shí)施例的方法流程圖，如圖2所示，所述方法包 括：
[0072] S11:獲取所述目標(biāo)工區(qū)在預(yù)設(shè)方向上的空間域算子，所述空間域算子的大小為第 二預(yù)設(shè)值。
[0073] 本實(shí)施例中的空間域算子可以包括方向差分空間域算子，其中，如表1-2所示，東 西方向的差分算子和南北方向的差分算子為索貝爾（Sobel)算子，索貝爾算子可以對(duì)網(wǎng)格 化后的重力異常數(shù)據(jù)進(jìn)行加權(quán)平均濾波，具有一定的噪聲抑制能力。如表3-4所示，所述北 東方向差分算子和北西方向差分算子為差分算子，同樣可以對(duì)網(wǎng)格化后的重力異常數(shù)據(jù)有 噪聲抑制作用。當(dāng)然，本實(shí)施例中的空間域算子不僅僅限于方向差分空間域算子，還可以包 括Roberts算子、Canny算子、Prewitt算子、邊緣檢測(cè)算子、二階算子等。一般地，空間域算子 的大小是固定的，例如索貝爾算子的大小為3X3矩陣，可以獲取空間域算子的大小作為第 二預(yù)設(shè)值，用于后續(xù)的選取數(shù)據(jù)窗口大小。
[0074] 需要說(shuō)明的是，根據(jù)目標(biāo)工區(qū)的區(qū)域特征，一般只選取預(yù)設(shè)方向的方向差分算子 對(duì)所述網(wǎng)格化后的重力異常數(shù)據(jù)進(jìn)行極值處理。例如，位于我國(guó)華南的目標(biāo)工區(qū)，可以利用 北東方向差分算子提取所述目標(biāo)工區(qū)的重力異常邊界極值信息。
[0076] S12:提取第一網(wǎng)格點(diǎn)中以待計(jì)算網(wǎng)格點(diǎn)為中心、大小為第二預(yù)設(shè)值的極值計(jì)算窗 □ 〇
[0077]本實(shí)施例中的第一網(wǎng)格點(diǎn)可以為網(wǎng)格化后的重力異常數(shù)據(jù)中網(wǎng)格點(diǎn)的集合，所述 待計(jì)算網(wǎng)格點(diǎn)屬于第一網(wǎng)格點(diǎn)中的單個(gè)網(wǎng)格點(diǎn)。提取以所述待計(jì)算網(wǎng)格點(diǎn)為中心、大小為 第二預(yù)設(shè)值的極值計(jì)算窗口。例如，提取以待計(jì)算網(wǎng)格點(diǎn)為中心的大小為3X3的極值計(jì) 算窗口，所述極值計(jì)算窗口數(shù)據(jù)為：
[0079] S13:將所述極值計(jì)算窗口與所述空間域算子的乘積作為所述待計(jì)算網(wǎng)格點(diǎn)的重 力異常邊界極值。
[0080] 本申請(qǐng)實(shí)施例中，可以將所述極值計(jì)算窗口中的數(shù)據(jù)與所述空間域算子的乘積作 為所述待計(jì)算網(wǎng)格點(diǎn)的重力異常邊界極值。例如，基于上述表1-4的方向差分空間域算子以 及以待計(jì)算網(wǎng)格點(diǎn)ai, j為中心的極值計(jì)算窗口，可以分別計(jì)算所述待計(jì)算網(wǎng)格點(diǎn)ai, j在東 西、南北、北東、北西方向的重力異常邊界極值。具體地，對(duì)于東西方向，重力異常邊界極值 的公式為：
[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]對(duì)于南北方向，重力異常邊界極值的公式為：
[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]對(duì)于北東方向，重力異常邊界極值的公式為：
[0085] Ai,j = ai+i,j-ai,j+i+ai,j-i-ai-ij；
[0086]對(duì)于北西方向，重力異常邊界極值的公式為：
[0087] Ai, j = ai, j+i-ai-1, j+ai+i, j_ai, j-1。
[0088] S14:遍歷所述第一網(wǎng)格點(diǎn)中的其他網(wǎng)格點(diǎn)，生成所述第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重 力異常邊界極值信息。
[0089] 在計(jì)算得到所述待計(jì)算網(wǎng)格點(diǎn)的重力異常邊界極值后，可以按照同樣的方法計(jì)算 所述第一網(wǎng)格點(diǎn)中的其他網(wǎng)格點(diǎn)，直至計(jì)算得到第一網(wǎng)格點(diǎn)中所有網(wǎng)格點(diǎn)的重力異常邊界 極值，生成所述第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重力異常邊界極值信息。
[0090] 需要說(shuō)明的是，在所述網(wǎng)格化后的重力異常數(shù)據(jù)中，處于邊緣位置的網(wǎng)格點(diǎn)一般 無(wú)法生成極值計(jì)算窗口，此時(shí)，需要對(duì)所述重力異常數(shù)據(jù)進(jìn)行擴(kuò)邊。本申請(qǐng)實(shí)施例中，可以 采用余弦衰減至零的方法對(duì)所述網(wǎng)格化后的重力異常數(shù)據(jù)進(jìn)行擴(kuò)邊。設(shè)網(wǎng)格化后的重力異 常數(shù)據(jù)為2￥(1 :111，1:11)，111為線數(shù)，11為每條線上的點(diǎn)數(shù)，設(shè)擴(kuò)邊數(shù)為1例如，當(dāng)1 = 2時(shí)，擴(kuò)邊 后的數(shù)據(jù)為ZK( 1-t :m+t，1-t: n+t)。設(shè)置ZK( i，j) = ZY( i，j)，（ 1 彡 i彡m，1 彡 j彡n)，則ZK中的 數(shù)據(jù)可以按如下公式擴(kuò)邊：
[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] 本實(shí)施例中，可以將網(wǎng)格化后的重力異常數(shù)據(jù)轉(zhuǎn)化為極值信息，用于后續(xù)的相干 增強(qiáng)處理，可以避免傳統(tǒng)頻率域算法中對(duì)網(wǎng)格化后的重力異常數(shù)據(jù)噪聲的放大作用。
[0096] S3:基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，計(jì)算第二網(wǎng)格點(diǎn)的相干增強(qiáng) 值。
[0097] 本實(shí)施例中的第二網(wǎng)格點(diǎn)為所述第一網(wǎng)格點(diǎn)中的任意一個(gè)網(wǎng)格點(diǎn)，基于所述目標(biāo) 工區(qū)中第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，可以計(jì)算第二網(wǎng)格點(diǎn)的相干增強(qiáng)值。圖3是本 申請(qǐng)?zhí)峁┑挠?jì)算網(wǎng)格點(diǎn)相干增強(qiáng)值的一種實(shí)施例的方法流程圖，如圖3所示，所述方法包 括：
[0098] S101:基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，建立以第二網(wǎng)格點(diǎn)為中心、 大小為第一預(yù)設(shè)值的增強(qiáng)值計(jì)算窗口。
[0099] 本實(shí)施例中，可以基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，建立以第二網(wǎng) 格點(diǎn)為中心、大小為第一預(yù)設(shè)值的增強(qiáng)值計(jì)算窗口。所述第一預(yù)設(shè)值大小可以根據(jù)整體網(wǎng) 格數(shù)目設(shè)置，本實(shí)施例中，所述第一預(yù)設(shè)值可以為5 X 5。
[0100] S102:提取與所述增強(qiáng)值計(jì)算窗口相鄰八個(gè)方向的、大小為所述第一預(yù)設(shè)值的相 鄰數(shù)據(jù)窗口。
[0101] 本實(shí)施例中，可以計(jì)算所述增強(qiáng)值計(jì)算窗口與相鄰八個(gè)方向的窗口的相關(guān)系數(shù) 值。因此，可以提取與所述增強(qiáng)值計(jì)算窗口相鄰八個(gè)方向的、大小相等的相鄰數(shù)據(jù)窗口。圖4 是本申請(qǐng)?zhí)峁┑脑鰪?qiáng)值計(jì)算窗口與相鄰八個(gè)方向的相鄰數(shù)據(jù)窗口的位置關(guān)系圖，如圖4所 示，其中，西北相鄰數(shù)據(jù)窗口是將所述增強(qiáng)值計(jì)算窗口的位置向上然后向左平移一個(gè)網(wǎng)格 得到的窗口位置。
[0102] S103:分別計(jì)算所述增強(qiáng)值計(jì)算窗口與所述相鄰八個(gè)方向的相鄰數(shù)據(jù)窗口的相關(guān) 系數(shù)值。
[0103] 本實(shí)施例中，可以分別計(jì)算所述增強(qiáng)值計(jì)算窗口與所述相鄰八個(gè)方向的相鄰數(shù)據(jù) 窗口的相關(guān)系數(shù)值。其中，所述相關(guān)系數(shù)值的計(jì)算公式可以為：
[0105] 其中，X可以為增強(qiáng)值計(jì)算窗口數(shù)據(jù)，Y可以為任一方位相鄰數(shù)據(jù)窗口內(nèi)的數(shù)據(jù) X、歹分別為對(duì)應(yīng)窗口全部數(shù)據(jù)的均值，n為第一預(yù)設(shè)值，k為1~8。
[0106] 例如，基于所述第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重力異常邊界極值信息數(shù)據(jù)，X和Y可以 為下式所示數(shù)據(jù)，根據(jù)上述相關(guān)系數(shù)值的計(jì)算公式，可以計(jì)算得到第二網(wǎng)格點(diǎn)乂 13與八個(gè)相 鄰數(shù)據(jù)窗口的相關(guān)系數(shù)值。
[0108] S104:將所述八個(gè)相關(guān)系數(shù)值的乘積的八次開(kāi)方值作為所述第二網(wǎng)格點(diǎn)的相干增 強(qiáng)值。
[0109]假設(shè)計(jì)算得到的8個(gè)相關(guān)系數(shù)值為^，所述第二網(wǎng)格點(diǎn)的相干增強(qiáng)值的計(jì)算公式 可以為：
[0111] 后續(xù)的，還可以對(duì)相干增強(qiáng)后的重力異常數(shù)據(jù)進(jìn)行二次相干增強(qiáng)，所述二次相干 增強(qiáng)可以包括下述步驟：
[0112] SS1:判斷所述第二網(wǎng)格點(diǎn)的相干增強(qiáng)值是否小于預(yù)設(shè)閾值；
[0113] SS2:當(dāng)判斷結(jié)果為是時(shí)，基于所述第一網(wǎng)格點(diǎn)的相干增強(qiáng)值信息，計(jì)算所述第二 網(wǎng)格點(diǎn)的二次相干增強(qiáng)值。
[0114] 當(dāng)判斷所述第二網(wǎng)格點(diǎn)的相干增強(qiáng)值小于預(yù)設(shè)閾值時(shí)，可以基于所述第一網(wǎng)格點(diǎn) 的相干增強(qiáng)值信息，按照上述S101-S104的方法步驟計(jì)算所述第二網(wǎng)格點(diǎn)的二次相干增強(qiáng) 值。
[0115] 下面通過(guò)一個(gè)具體的應(yīng)用場(chǎng)景說(shuō)明上述實(shí)施例方法，圖5是本發(fā)明提供的我國(guó)東 北目標(biāo)工區(qū)A的實(shí)測(cè)布格重力異常圖，如圖5所示，圖中顏色的深淺表示重力異常值的高低， 顏色越深，表示重力異常值越高，反之，顏色越淺，表示重力異常值越低。利用向上延拓2千 米的方法去除圖5中布格重力異常圖的噪聲，得到圖6所示的對(duì)工區(qū)A進(jìn)行向上延拓去噪后 的布格重力異常圖。根據(jù)工區(qū)A已知的地質(zhì)情況，采用北東、北西方向差分算子提取北東、北 西兩個(gè)方向的重力異常邊界極值信息?；谔崛〉谋睎|、北西兩個(gè)方向重力異常邊界極值 信息，對(duì)圖6中的各個(gè)網(wǎng)格點(diǎn)進(jìn)行相干增強(qiáng)處理，得到工區(qū)A的第一次相干增強(qiáng)數(shù)據(jù)。然后采 用閾值0.985對(duì)所述第一次相干增強(qiáng)數(shù)據(jù)進(jìn)行二次相干增強(qiáng)處理，得到該區(qū)北東向和北西 向的重力異常邊界信息，圖7是將工區(qū)A北東、北西兩個(gè)方向重力異常邊界信息融合后的斷 裂走勢(shì)圖，反映了該地區(qū)斷陷斷裂的主體分布情況。圖8是本發(fā)明提供的利用總水平導(dǎo)數(shù)法 對(duì)工區(qū)A提取和增強(qiáng)的重力異常邊界信息圖，將圖7與圖8比較，可以發(fā)現(xiàn)本發(fā)明所提取和增 強(qiáng)的重力異常邊界信息更加清晰、連續(xù)。另外，如圖7所示本發(fā)明提供的實(shí)施例方法不僅可 以提取和增強(qiáng)工區(qū)A北東向的主要斷裂，而且還對(duì)正交方向（北西向）的扭曲和錯(cuò)斷斷裂反 映的比較清楚，表現(xiàn)了比傳統(tǒng)的倒數(shù)類算法具有更強(qiáng)的處理能力。
[0116]本申請(qǐng)?zhí)峁┑闹亓Ξ惓＿吔缭鰪?qiáng)方法，可以提取網(wǎng)格化后的重力異常數(shù)據(jù)在預(yù)設(shè) 方向上的重力異常邊界極值信息，所述重力異常邊界的極值一般對(duì)應(yīng)著模型的邊界點(diǎn)。然 后基于所述重力異常邊界極值信息，可以計(jì)算網(wǎng)格點(diǎn)的相干增強(qiáng)值，進(jìn)一步獲取模型邊界 點(diǎn)的增強(qiáng)值。本申請(qǐng)?zhí)峁┑闹亓Ξ惓＿吔缭鰪?qiáng)方法對(duì)所述重力異常數(shù)據(jù)進(jìn)行空間域處理， 避免傳統(tǒng)頻率域算法對(duì)數(shù)據(jù)噪聲的放大作用，可以獲取清晰、連續(xù)的重力異常邊界信息，對(duì) 目標(biāo)工區(qū)地層斷裂數(shù)據(jù)的研究提供有力的基礎(chǔ)資料。
[0117]本發(fā)明另一方面還提供一種重力異常邊界增強(qiáng)的裝置，圖9是本申請(qǐng)?zhí)峁┑闹亓?異常邊界增強(qiáng)裝置的一種實(shí)施例的模塊結(jié)構(gòu)示意圖，如圖9所示，所述裝置90可以包括： [0118]網(wǎng)格化單元91，用于獲取目標(biāo)工區(qū)的重力異常數(shù)據(jù)，并對(duì)所述重力異常數(shù)據(jù)進(jìn)行 網(wǎng)格化；
[0119] 邊界極值計(jì)算單元92,用于采用空間域算子逐個(gè)提取第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重 力異常邊界極值信息；
[0120] 增強(qiáng)值計(jì)算單元93,用于基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，計(jì)算第 二網(wǎng)格點(diǎn)的相干增強(qiáng)值。
[0121] 本申請(qǐng)?zhí)峁┑闹亓Ξ惓＿吔缭鰪?qiáng)裝置，可以提取網(wǎng)格化后的重力異常數(shù)據(jù)在預(yù)設(shè) 方向上的重力異常邊界極值信息，所述重力異常邊界的極值一般對(duì)應(yīng)著模型的邊界點(diǎn)。然 后基于所述重力異常邊界極值信息，可以計(jì)算網(wǎng)格點(diǎn)的相干增強(qiáng)值，進(jìn)一步獲取模型邊界 點(diǎn)的增強(qiáng)值。本申請(qǐng)?zhí)峁┑闹亓Ξ惓＿吔缭鰪?qiáng)方法對(duì)所述重力異常數(shù)據(jù)進(jìn)行空間域處理， 避免傳統(tǒng)頻率域算法對(duì)數(shù)據(jù)噪聲的放大作用，可以獲取清晰、連續(xù)的重力異常邊界信息，對(duì) 目標(biāo)工區(qū)地層斷裂數(shù)據(jù)的研究提供有力的基礎(chǔ)資料。
[0122] 圖10是本申請(qǐng)?zhí)峁┑脑鰪?qiáng)值計(jì)算單元的一種實(shí)施例的模塊結(jié)構(gòu)示意圖，如圖10所 示，所述增強(qiáng)值計(jì)算單元93可以包括：
[0123] 增強(qiáng)值窗口建立單元101，用于基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，建 立以第二網(wǎng)格點(diǎn)為中心、大小為第一預(yù)設(shè)值的增強(qiáng)值計(jì)算窗口；
[0124] 相鄰窗口提取單元102，用于提取與所述增強(qiáng)值計(jì)算窗口相鄰八個(gè)方向的、大小為 所述第一預(yù)設(shè)值的相鄰數(shù)據(jù)窗口；
[0125] 相關(guān)系數(shù)計(jì)算單元103，用于分別計(jì)算所述增強(qiáng)值計(jì)算窗口與所述相鄰八個(gè)方向 的相鄰數(shù)據(jù)窗口的相關(guān)系數(shù)值；
[0126] 開(kāi)方計(jì)算單元104,用于將所述八個(gè)相關(guān)系數(shù)值的乘積的八次開(kāi)方值作為所述第 二網(wǎng)格點(diǎn)的相干增強(qiáng)值。
[0127] 可選的，在本發(fā)明的一個(gè)實(shí)施例中，所述裝置還包括：
[0128] 去噪單元，用于采用低通濾波法去除所述重力異常資料中的噪聲。
[0129] 圖11是本申請(qǐng)?zhí)峁┑闹亓Ξ惓＿吔缭鰪?qiáng)裝置的另一種實(shí)施例的模塊結(jié)構(gòu)示意圖， 如圖11所示，所述裝置110還可以包括：
[0130] 增強(qiáng)值判斷單元111，用于判斷所述第二網(wǎng)格點(diǎn)的相干增強(qiáng)值是否小于預(yù)設(shè)閾值；
[0131] 二次增強(qiáng)單元112,用于當(dāng)所述增強(qiáng)值判斷單元的判斷結(jié)果為是時(shí)，基于所述第一 網(wǎng)格點(diǎn)的相干增強(qiáng)值信息，計(jì)算所述第二網(wǎng)格點(diǎn)的二次相干增強(qiáng)值。
[0132] 圖12是本申請(qǐng)?zhí)峁┑倪吔鐦O值計(jì)算單元的一種實(shí)施例的模塊結(jié)構(gòu)示意圖，如圖12 所示，所述邊界極值計(jì)算單元92可以包括：
[0133] 算子獲取單元121，用于獲取所述目標(biāo)工區(qū)在預(yù)設(shè)方向上的空間域算子，所述空間 域算子的大小為第二預(yù)設(shè)值；
[0134] 極值窗口提取單元122,用于提取第一網(wǎng)格點(diǎn)中以待計(jì)算網(wǎng)格點(diǎn)為中心、大小為第 二預(yù)設(shè)值的極值計(jì)算窗口；
[0135] 極值計(jì)算單元123,用于將所述極值計(jì)算窗口與所述空間域算子的乘積作為所述 待計(jì)算網(wǎng)格點(diǎn)的重力異常邊界極值；
[0136] 遍歷計(jì)算單元124,用于遍歷所述第一網(wǎng)格點(diǎn)中的其他網(wǎng)格點(diǎn)，生成所述第一網(wǎng)格 點(diǎn)在預(yù)設(shè)方向的重力異常邊界極值信息。
[0137] 本說(shuō)明書(shū)中的各個(gè)實(shí)施例均采用遞進(jìn)的方式描述，每個(gè)實(shí)施例重點(diǎn)說(shuō)明的都是與 其他實(shí)施例的不同之處，各個(gè)實(shí)施例之間相同相似的部分互相參見(jiàn)即可。尤其，對(duì)于系統(tǒng)實(shí) 施例而言，由于其基本相似于方法實(shí)施例，所以描述的比較簡(jiǎn)單，相關(guān)之處參見(jiàn)方法實(shí)施例 的部分說(shuō)明即可。
[0138] 雖然通過(guò)實(shí)施例描繪了本申請(qǐng)，本領(lǐng)域普通技術(shù)人員知道，本申請(qǐng)有許多變形和 變化而不脫離本申請(qǐng)的精神，希望所附的權(quán)利要求包括這些變形和變化而不脫離本申請(qǐng)的 精神。
[0139] 盡管本申請(qǐng)內(nèi)容中提到圖片裁剪、RGB通道顏色轉(zhuǎn)換、微調(diào)(fine-tuning)、卷積神 經(jīng)網(wǎng)絡(luò)模型層結(jié)構(gòu)、卷積、損失反饋等的圖片信息處理、神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)的描述，但是，本 申請(qǐng)并不局限于必須是完全標(biāo)準(zhǔn)或者所提及的方式的信息處理、神經(jīng)網(wǎng)絡(luò)模型結(jié)構(gòu)的情 況。本申請(qǐng)中各個(gè)實(shí)施例所涉及的上述描述僅是本申請(qǐng)中的一些實(shí)施例中的應(yīng)用，在某些 標(biāo)準(zhǔn)、模型、方法的基礎(chǔ)上略加修改后的實(shí)施方式也可以實(shí)行上述本申請(qǐng)各實(shí)施例的方案。 當(dāng)然，在符合本申請(qǐng)上述各實(shí)施例的中所述的處理方法步驟的其他無(wú)創(chuàng)造性的變形，仍然 可以實(shí)現(xiàn)相同的申請(qǐng)，在此不再贅述。
[0140] 雖然本申請(qǐng)?zhí)峁┝巳鐚?shí)施例或流程圖所述的方法操作步驟，但基于常規(guī)或者無(wú)創(chuàng) 造性的勞動(dòng)可以包括更多或者更少的操作步驟。實(shí)施例中列舉的步驟順序僅僅為眾多步驟 執(zhí)行順序中的一種方式，不代表唯一的執(zhí)行順序。在實(shí)際中的裝置或客戶端產(chǎn)品執(zhí)行時(shí)，可 以按照實(shí)施例或者附圖所示的方法順序執(zhí)行或者并行執(zhí)行(例如并行處理器或者多線程處 理的環(huán)境）。
[0141] 上述實(shí)施例闡明的裝置或模塊，具體可以由計(jì)算機(jī)芯片或?qū)嶓w實(shí)現(xiàn)，或者由具有 某種功能的產(chǎn)品來(lái)實(shí)現(xiàn)。為了描述的方便，描述以上裝置時(shí)以功能分為各種模塊分別描述。 在實(shí)施本申請(qǐng)時(shí)可以把各模塊的功能在同一個(gè)或多個(gè)軟件和/或硬件中實(shí)現(xiàn)。當(dāng)然，也可以 將實(shí)現(xiàn)某功能的模塊由多個(gè)子模塊或子單元組合實(shí)現(xiàn)。
[0142] 本申請(qǐng)中所述的方法、裝置或模塊可以以計(jì)算機(jī)可讀程序代碼方式實(shí)現(xiàn)控制器按 任何適當(dāng)?shù)姆绞綄?shí)現(xiàn)，例如，控制器可以采取例如微處理器或處理器以及存儲(chǔ)可由該(微） 處理器執(zhí)行的計(jì)算機(jī)可讀程序代碼(例如軟件或固件）的計(jì)算機(jī)可讀介質(zhì)、邏輯門(mén)、開(kāi)關(guān)、專 用集成電路(Application Specific Integrated Circuit，ASIC)、可編程邏輯控制器和嵌 入微控制器的形式，控制器的例子包括但不限于以下微控制器:ARC 62f5D、Atmel AT91SAM、 Microchip PIC18F26K20以及Silicone Labs C8051F320,存儲(chǔ)器控制器還可以被實(shí)現(xiàn)為存 儲(chǔ)器的控制邏輯的一部分。本領(lǐng)域技術(shù)人員也知道，除了以純計(jì)算機(jī)可讀程序代碼方式實(shí) 現(xiàn)控制器以外，完全可以通過(guò)將方法步驟進(jìn)行邏輯編程來(lái)使得控制器以邏輯門(mén)、開(kāi)關(guān)、專用 集成電路、可編程邏輯控制器和嵌入微控制器等的形式來(lái)實(shí)現(xiàn)相同功能。因此這種控制器 可以被認(rèn)為是一種硬件部件，而對(duì)其內(nèi)部包括的用于實(shí)現(xiàn)各種功能的裝置也可以視為硬件 部件內(nèi)的結(jié)構(gòu)?；蛘呱踔粒梢詫⒂糜趯?shí)現(xiàn)各種功能的裝置視為既可以是實(shí)現(xiàn)方法的軟件 模塊又可以是硬件部件內(nèi)的結(jié)構(gòu)。
[0143] 本申請(qǐng)所述裝置中的部分模塊可以在由計(jì)算機(jī)執(zhí)行的計(jì)算機(jī)可執(zhí)行指令的一般 上下文中描述，例如程序模塊。一般地，程序模塊包括執(zhí)行特定任務(wù)或?qū)崿F(xiàn)特定抽象數(shù)據(jù)類 型的例程、程序、對(duì)象、組件、數(shù)據(jù)結(jié)構(gòu)、類等等。也可以在分布式計(jì)算環(huán)境中實(shí)踐本申請(qǐng)，在 這些分布式計(jì)算環(huán)境中，由通過(guò)通信網(wǎng)絡(luò)而被連接的遠(yuǎn)程處理設(shè)備來(lái)執(zhí)行任務(wù)。在分布式 計(jì)算環(huán)境中，程序模塊可以位于包括存儲(chǔ)設(shè)備在內(nèi)的本地和遠(yuǎn)程計(jì)算機(jī)存儲(chǔ)介質(zhì)中。
[0144] 通過(guò)以上的實(shí)施方式的描述可知，本領(lǐng)域的技術(shù)人員可以清楚地了解到本申請(qǐng)可 借助軟件加必需的硬件的方式來(lái)實(shí)現(xiàn)。基于這樣的理解，本申請(qǐng)的技術(shù)方案本質(zhì)上或者說(shuō) 對(duì)現(xiàn)有技術(shù)做出貢獻(xiàn)的部分可以以軟件產(chǎn)品的形式體現(xiàn)出來(lái)，也可以通過(guò)數(shù)據(jù)迀移的實(shí)施 過(guò)程中體現(xiàn)出來(lái)。該計(jì)算機(jī)軟件產(chǎn)品可以存儲(chǔ)在存儲(chǔ)介質(zhì)中，如R0M/RAM、磁碟、光盤(pán)等，包 括若干指令用以使得一臺(tái)計(jì)算機(jī)設(shè)備(可以是個(gè)人計(jì)算機(jī)，移動(dòng)終端，服務(wù)器，或者網(wǎng)絡(luò)設(shè) 備等)執(zhí)行本申請(qǐng)各個(gè)實(shí)施例或者實(shí)施例的某些部分所述的方法。
[0145] 本說(shuō)明書(shū)中的各個(gè)實(shí)施例采用遞進(jìn)的方式描述，各個(gè)實(shí)施例之間相同或相似的部 分互相參見(jiàn)即可，每個(gè)實(shí)施例重點(diǎn)說(shuō)明的都是與其他實(shí)施例的不同之處。本申請(qǐng)的全部或 者部分可用于眾多通用或?qū)Ｓ玫挠?jì)算機(jī)系統(tǒng)環(huán)境或配置中。例如:個(gè)人計(jì)算機(jī)、服務(wù)器計(jì)算 機(jī)、手持設(shè)備或便攜式設(shè)備、平板型設(shè)備、移動(dòng)通信終端、多處理器系統(tǒng)、基于微處理器的系 統(tǒng)、可編程的電子設(shè)備、網(wǎng)絡(luò)PC、小型計(jì)算機(jī)、大型計(jì)算機(jī)、包括以上任何系統(tǒng)或設(shè)備的分布 式計(jì)算環(huán)境等等。
[0146] 雖然通過(guò)實(shí)施例描繪了本申請(qǐng)，本領(lǐng)域普通技術(shù)人員知道，本申請(qǐng)有許多變形和 變化而不脫離本申請(qǐng)的精神，希望所附的權(quán)利要求包括這些變形和變化而不脫離本申請(qǐng)的 精神。
【主權(quán)項(xiàng)】
1. 一種重力異常邊界增強(qiáng)的方法，其特征在于，所述方法包括： 獲取目標(biāo)工區(qū)的重力異常數(shù)據(jù)，并對(duì)所述重力異常數(shù)據(jù)進(jìn)行網(wǎng)格化； 采用空間域算子逐個(gè)提取第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重力異常邊界極值信息； 基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，計(jì)算第二網(wǎng)格點(diǎn)的相干增強(qiáng)值。2. 根據(jù)權(quán)利要求1所述的一種重力異常邊界增強(qiáng)的方法，其特征在于，所述基于所述第 一網(wǎng)格點(diǎn)的重力異常邊界極值信息，計(jì)算第二網(wǎng)格點(diǎn)的相干增強(qiáng)值，包括： 基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，建立以第二網(wǎng)格點(diǎn)為中心、大小為第 一預(yù)設(shè)值的增強(qiáng)值計(jì)算窗口； 提取與所述增強(qiáng)值計(jì)算窗口相鄰八個(gè)方向的、大小為所述第一預(yù)設(shè)值的相鄰數(shù)據(jù)窗 P； 分別計(jì)算所述增強(qiáng)值計(jì)算窗口與所述相鄰八個(gè)方向的相鄰數(shù)據(jù)窗口的相關(guān)系數(shù)值； 將所述八個(gè)相關(guān)系數(shù)值的乘積的八次開(kāi)方值作為所述第二網(wǎng)格點(diǎn)的相干增強(qiáng)值。3. 根據(jù)權(quán)利要求1所述的一種重力異常邊界增強(qiáng)的方法，其特征在于，在所述對(duì)所述重 力資料進(jìn)行網(wǎng)格化之后，還包括： 采用低通濾波法去除所述重力異常資料中的噪聲。4. 根據(jù)權(quán)利要求1所述的一種重力異常邊界增強(qiáng)的方法，其特征在于，所述方法還包 括： 判斷所述第二網(wǎng)格點(diǎn)的相干增強(qiáng)值是否小于預(yù)設(shè)閾值； 當(dāng)判斷結(jié)果為是時(shí)，基于所述第一網(wǎng)格點(diǎn)的相干增強(qiáng)值信息，計(jì)算所述第二網(wǎng)格點(diǎn)的 二次相干增強(qiáng)值。5. 根據(jù)權(quán)利要求1所述的一種重力異常邊界增強(qiáng)的方法，其特征在于，所述采用空間域 算子逐個(gè)提取第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重力異常邊界極值信息包括： 獲取所述目標(biāo)工區(qū)在預(yù)設(shè)方向上的空間域算子，所述空間域算子的大小為第二預(yù)設(shè) 值； 提取第一網(wǎng)格點(diǎn)中以待計(jì)算網(wǎng)格點(diǎn)為中心、大小為第二預(yù)設(shè)值的極值計(jì)算窗口； 將所述極值計(jì)算窗口與所述空間域算子的乘積作為所述待計(jì)算網(wǎng)格點(diǎn)的重力異常邊 界極值； 遍歷所述第一網(wǎng)格點(diǎn)中的其他網(wǎng)格點(diǎn)，生成所述第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重力異常邊 界極值信息。6. -種重力異常邊界增強(qiáng)的裝置，其特征在于，所述裝置包括： 網(wǎng)格化單元，用于獲取目標(biāo)工區(qū)的重力異常數(shù)據(jù)，并對(duì)所述重力異常數(shù)據(jù)進(jìn)行網(wǎng)格化； 邊界極值計(jì)算單元，用于采用空間域算子逐個(gè)提取第一網(wǎng)格點(diǎn)在預(yù)設(shè)方向的重力異常 邊界極值信息； 增強(qiáng)值計(jì)算單元，用于基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，計(jì)算第二網(wǎng)格 點(diǎn)的相干增強(qiáng)值。7. 根據(jù)權(quán)利要求6所述的一種重力異常邊界增強(qiáng)的裝置，其特征在于，所述增強(qiáng)值計(jì)算 單元包括： 增強(qiáng)值窗口建立單元，用于基于所述第一網(wǎng)格點(diǎn)的重力異常邊界極值信息，建立以第 二網(wǎng)格點(diǎn)為中心、大小為第一預(yù)設(shè)值的增強(qiáng)值計(jì)算窗口； 相鄰窗口提取單元，用于提取與所述增強(qiáng)值計(jì)算窗口相鄰八個(gè)方向的、大小為所述第 一預(yù)設(shè)值的相鄰數(shù)據(jù)窗口； 相關(guān)系數(shù)計(jì)算單元，用于分別計(jì)算所述增強(qiáng)值計(jì)算窗口與所述相鄰八個(gè)方向的相鄰數(shù) 據(jù)窗口的相關(guān)系數(shù)值； 開(kāi)方計(jì)算單元，用于將所述八個(gè)相關(guān)系數(shù)值的乘積的八次開(kāi)方值作為所述第二網(wǎng)格點(diǎn) 的相干增強(qiáng)值。8. 根據(jù)權(quán)利要求6所述的一種重力異常邊界增強(qiáng)的裝置，其特征在于，所述裝置還包 括： 去噪單元，用于采用低通濾波法去除所述重力異常資料中的噪聲。9. 根據(jù)權(quán)利要求6所述的一種重力異常邊界增強(qiáng)的裝置，其特征在于，所述裝置還包 括： 增強(qiáng)值判斷單元，用于判斷所述第二網(wǎng)格點(diǎn)的相干增強(qiáng)值是否小于預(yù)設(shè)閾值； 二次增強(qiáng)單元，用于當(dāng)所述增強(qiáng)值判斷單元的判斷結(jié)果為是時(shí)，基于所述第一網(wǎng)格點(diǎn) 的相干增強(qiáng)值信息，計(jì)算所述第二網(wǎng)格點(diǎn)的二次相干增強(qiáng)值。10. 根據(jù)權(quán)利要6所述的一種重力異常邊界增強(qiáng)的裝置，其特征在于，所述邊界極值計(jì) 算單元包括： 算子獲取單元，用于獲取所述目標(biāo)工區(qū)在預(yù)設(shè)方向上的空間域算子，所述空間域算子 的大小為第二預(yù)設(shè)值； 極值窗口提取單元，用于提取第一網(wǎng)格點(diǎn)中以待計(jì)算網(wǎng)格點(diǎn)為中心、大小為第二預(yù)設(shè) 值的極值計(jì)算窗口； 極值計(jì)算單元，用于將所述極值計(jì)算窗口與所述空間域算子的乘積作為所述待計(jì)算網(wǎng) 格點(diǎn)的重力異常邊界極值； 遍歷計(jì)算單元，用于遍歷所述第一網(wǎng)格點(diǎn)中的其他網(wǎng)格點(diǎn)，生成所述第一網(wǎng)格點(diǎn)在預(yù) 設(shè)方向的重力異常邊界極值信息。
【文檔編號(hào)】G01V1/36GK105891889SQ201610197300
【公開(kāi)日】2016年8月24日
【申請(qǐng)日】2016年3月31日
【發(fā)明人】趙文舉, 劉云祥, 趙荔
【申請(qǐng)人】中國(guó)石油天然氣集團(tuán)公司, 中國(guó)石油集團(tuán)東方地球物理勘探有限責(zé)任公司