亚洲狠狠干,亚洲国产福利精品一区二区,国产八区,激情文学亚洲色图

一種重力異常邊界增強(qiáng)的方法及裝置的制造方法

文檔序號(hào):10533443閱讀:471來(lái)源:國(guó)知局
一種重力異常邊界增強(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)行布格校正,若在布格校 正之后仍然存在的異常稱為布格重力異常??梢詫⑿U蟮闹亓τ^測(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ū)S玫挠?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é)任公司
網(wǎng)友詢問(wèn)留言 已有0條留言
  • 還沒(méi)有人留言評(píng)論。精彩留言會(huì)獲得點(diǎn)贊!
1