• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    利用加強(qiáng)水平方向總水平導(dǎo)數(shù)對位場全張量數(shù)據(jù)進(jìn)行邊界識別

    2015-03-16 11:07:14袁園黃大年余青露
    地球物理學(xué)報 2015年7期
    關(guān)鍵詞:方向水平

    袁園, 黃大年, 余青露

    1 吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院, 長春 130026 2 國家海洋局第二海洋研究所, 杭州 310012 3 國家海洋局海底科學(xué)重點實驗室, 杭州 310012 4 中國石化石油物探技術(shù)研究院, 南京 211103

    ?

    利用加強(qiáng)水平方向總水平導(dǎo)數(shù)對位場全張量數(shù)據(jù)進(jìn)行邊界識別

    袁園1,2,3, 黃大年1*, 余青露4

    1 吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院, 長春 130026 2 國家海洋局第二海洋研究所, 杭州 310012 3 國家海洋局海底科學(xué)重點實驗室, 杭州 310012 4 中國石化石油物探技術(shù)研究院, 南京 211103

    位場全張量梯度數(shù)據(jù)以其信息量大、含有更高頻的信號成分,能更好地描述小的異常特征等優(yōu)點在地球物理領(lǐng)域中得到廣泛應(yīng)用.邊界檢測是位場解釋中不可缺少的任務(wù),需要新的邊界探測器來處理位場梯度張量數(shù)據(jù).為了充分利用位場梯度張量數(shù)據(jù)的多信息成分,本文定義了方向總水平導(dǎo)數(shù)和加強(qiáng)方向總水平導(dǎo)數(shù),并利用其定義新的邊界檢測器.為了能同時顯示不同振幅大小異常的邊界,本文對其進(jìn)行了歸一化處理.通過模型試驗,證明了歸一化方法能更加清晰準(zhǔn)確地顯示淺部和深部的地質(zhì)體邊界信息.最后將該邊界檢測方法用于加拿大圣喬治灣實際測得全張量重力梯度數(shù)據(jù)和中國朱日和地區(qū)的磁異常數(shù)據(jù)中,并得到了較好的邊界檢測結(jié)果.

    邊界檢測; 位場梯度張量; 加強(qiáng)方向總水平導(dǎo)數(shù); 歸一化

    1 引言

    邊界檢測是位場數(shù)據(jù)解釋中的重要組成部分.目前,有許多方法用于邊界檢測和邊界加強(qiáng).水平導(dǎo)數(shù)和垂向?qū)?shù)經(jīng)常用于加強(qiáng)這些潛在特征.垂向?qū)?shù)通常被用于描述密度體或磁性體的邊界(Evjen, 1936).后來人們發(fā)現(xiàn)總水平導(dǎo)數(shù)和解析信號振幅的最大值直接描述密度和磁化強(qiáng)度邊界的檢測(Cordell and Grauch, 1985; Cordell, 1979; Nabighian, 1984).Hsu等(1996)提出加強(qiáng)解析信號方法進(jìn)行到高階導(dǎo)數(shù)以增加該方法的求解能力.Cooper和Cowan (2004)使用不同階的垂向?qū)?shù)進(jìn)行增強(qiáng)地球物理異常的細(xì)節(jié).

    Miller和Singh(1994)使用垂向?qū)?shù)與總水平導(dǎo)數(shù)的比來進(jìn)行邊界檢測.它是第一個均衡濾波器,能均衡不同振幅強(qiáng)度的異常信息.Rajagopalan和Milligan (1995)使用自動增益控制技術(shù)均衡不同大小的振幅.然而,這種方法與窗口大小的選擇相關(guān).Verduzco等(2004)用傾斜角的總水平導(dǎo)數(shù)(THDR)來進(jìn)行邊界檢測.Wijns等(2005)提出了Theta圖濾波器進(jìn)行邊界探測,它是用解析信號振幅來對總水平導(dǎo)數(shù)進(jìn)行歸一化處理.這些均衡邊界探測器雖然能同時顯示不同振幅大小的異常的邊界信息,但識別出來的邊界分辨率較低.馬國慶等(2012)利用高階導(dǎo)數(shù)來定義均衡濾波器,提高了對地質(zhì)體邊界的分辨率.但是,當(dāng)測量的數(shù)據(jù)中同時包含正負(fù)異常體的情況下,這些均衡探測器可能會引入錯誤的邊界信息.

    近十幾年,位場梯度張量數(shù)據(jù)被廣泛使用.位場梯度張量數(shù)據(jù)可以通過測量或者數(shù)值計算得到,它提供了9個張量成分,包含了比傳統(tǒng)位場數(shù)據(jù)更高頻的信號成分,因此,對于梯度張量數(shù)據(jù)的解釋能得到更高分辨率和更加詳細(xì)的地質(zhì)構(gòu)造.這要求新的方法來解釋梯度張量數(shù)據(jù),尤其是對小地質(zhì)構(gòu)造的邊界檢測.Cooper 和Cowan (2006)以及Oru?和Keskinsezer(2008)提出了方向傾斜角來探測邊界,但是只有垂直方向的傾斜角,即傳統(tǒng)的傾斜角方法,能得到較好的檢測結(jié)果.Beiki(2010)和Mikhailov等(2007)定義了方向解析信號來處理重力梯度張量數(shù)據(jù).Beiki(2010)利用方向解析信號進(jìn)行邊界檢測,比傳統(tǒng)的解析信號方法能得到更好的結(jié)果,但是該方法也不能同時描述大振幅和小振幅異常體.近年來,一些學(xué)者利用重力梯度張量數(shù)據(jù)的特征值來進(jìn)行邊界識別(Oru? et al., 2013; Sertcelik and Kafadar, 2012; Zhou et al., 2013).

    Marson和Klingele(1993)指出總水平導(dǎo)數(shù)比解析信號有更高的分辨率.因此,本文定義了方向總水平導(dǎo)數(shù)和加強(qiáng)方向總水平導(dǎo)數(shù),并根據(jù)其定義不同的邊界探測器.為了同時顯示不同振幅大小異常的邊界,給出了歸一化方法.

    2 方向總水平導(dǎo)數(shù)

    位場梯度張量數(shù)據(jù)為位場U在空間上的二階偏導(dǎo)數(shù),張量矩陣的形式:

    (1)

    Beiki(2010)和Mikhailov等(2007)定義了方向解析信號,其表達(dá)式為

    (2)

    (3)

    (4)

    其中,下標(biāo)x,y和z表示方向.

    這里,我們在x,y和z方向上定義方向總水平導(dǎo)數(shù),其表達(dá)式分別為

    (5)

    (6)

    (7)

    其中,下標(biāo)x,y和z表示方向.THDx的最大值描繪地質(zhì)體的N-S向邊界;THDy的最大值描繪地質(zhì)體的E-W向邊界;垂直方向上的總水平導(dǎo)數(shù)THDz與傳統(tǒng)使用的總水平導(dǎo)數(shù)相同,振幅的最大值描繪地質(zhì)體的邊界.

    圖1比較分析了方向解析信號振幅和方向總水平導(dǎo)數(shù)的邊界識別能力.圖1a為一垂直棱柱體產(chǎn)生的重力異常.圖1b—1d為x,y和z方向上的方向解析信號振幅;圖1e—1g為x,y和z方向上的方向總水平導(dǎo)數(shù).通過比較可知,方向總水平導(dǎo)數(shù)相對于方向解析信號振幅對于邊界識別有著更高的分辨率.因此,根據(jù)水平方向的總水平導(dǎo)數(shù)來定義邊界探測器有著更高的分辨率.利用THDx和THDy定義地質(zhì)體邊界探測器:

    (8)

    EDT的最大值定義地質(zhì)體的邊界.

    3 加強(qiáng)方向總水平導(dǎo)數(shù)

    Hsu等(1996)應(yīng)用高階導(dǎo)數(shù)定義加強(qiáng)解析信號振幅來提高邊界的分辨率.為了提高邊界檢測的分辨率,我們使用高階導(dǎo)數(shù)來定義方向總水平導(dǎo)數(shù),叫做加強(qiáng)方向總水平導(dǎo)數(shù).

    x,y和z方向上的加強(qiáng)方向總水平導(dǎo)數(shù)分別為

    (9)

    (10)

    (11)

    其中,下標(biāo)x,y和z表示方向.上標(biāo)n表示階數(shù),大于等于1.當(dāng)n=1時,與方向總水平導(dǎo)數(shù)相同.

    同樣,加強(qiáng)方向總水平導(dǎo)數(shù)ETHDx的最大值描繪地質(zhì)體的N-S向邊界;ETHDy的最大值描繪地質(zhì)體的E-W向邊界.因此,我們可以根據(jù)加強(qiáng)方向總水平導(dǎo)數(shù)來定義地質(zhì)體邊界探測器:

    (12)

    EEDTn的最大值定義地質(zhì)體的邊界.通常情況下,我們只取n=1,2.因為位場數(shù)據(jù)高階導(dǎo)數(shù)的計算將放大噪聲的影響.為了識別出來的邊界更加準(zhǔn)確,我們應(yīng)避免高階導(dǎo)數(shù)的計算.圖1h—1j為x,y和z方向上的二階加強(qiáng)方向總水平導(dǎo)數(shù).可以看出,相對于方向總水平導(dǎo)數(shù)有著更高的分辨率.這里,我們用EEDT來表示由二階加強(qiáng)方向總水平導(dǎo)數(shù)定義的邊界探測器EEDT2.在EEDT中,需要計算二階加強(qiáng)方向總水平的導(dǎo)數(shù),然而,垂向?qū)?shù)的直接計算會增加噪聲的影響.我們用如下表達(dá)式來計算二階ETHDx和ETHDy導(dǎo)數(shù):

    (13)

    (14)

    (15)

    4 歸一化方法

    邊界探測器EDT和EEDT都為非均衡濾波器,不能很好地同時顯示不同振幅大小的異常的邊界.為了避免這一缺點,對其采用了歸一化方法.

    本文采用位場數(shù)據(jù)的一階垂向?qū)?shù)來對EDT進(jìn)行歸一化處理,表達(dá)式為

    (16)

    這里引入常數(shù)p,一般大小在0.001到0.1.常數(shù)p的引入是為了避免產(chǎn)生額外錯誤的邊界.大的p值將降低均衡效果,小的p值將加強(qiáng)弱振幅異常的邊界.

    邊界探測器EEDT的歸一化表達(dá)式為NEEDT=

    (17)

    同樣,參數(shù)p與(13)式相同.這里,用如下表達(dá)式來計算Uzzz:

    (18)

    為了驗證本文方法的可行性,我們選用兩種常用的方法進(jìn)行對比.它們是傾角總水平導(dǎo)數(shù)(THDR)和Theta圖(Verduzcoetal., 2004;Wijnsetal., 2005).

    5 模型試驗

    建立如下重力異常模型一:模型中包含4個垂直棱柱體,各個棱柱體的上頂埋深分別為1km,1km,2km,2km.圖2顯示了模型體的平面圖和3維圖.各個棱柱體的厚度都為1km,剩余密度為200kg·m-3.圖3為合成的重力異常模型.重力異常模型的梯度張量數(shù)據(jù)可以通過數(shù)值計算和實際測量得到.這里,合成模型的重力梯度張量值通過理論公式計算得到(Forsberg, 1984)(圖4a—4c).圖4d—4i分別為EDT,EEDT,NEDT,NEEDT,THDR和Theta圖的邊界結(jié)果.可以看出,EDT和EEDT為非均衡濾波器,對于深部異常體的邊界檢測能力較弱,顯示結(jié)果不清晰.NEDT和NEEDT相對于兩種常用的邊界探測器THDR和Theta圖法更加清晰和準(zhǔn)確.

    為了進(jìn)一步驗證方法的有效性,建立更加復(fù)雜的重力異常模型二,其中包含正負(fù)異常,更加符合實際地質(zhì)情況.模型體的大小和位置與圖2相同.唯一不同的是棱柱體1和3的剩余密度為200kg·m-3;棱柱體2和4的剩余密度為-200kg·m-3.合成的重力異常見圖5.圖6為計算的梯度張量和相應(yīng)的邊界檢測結(jié)果.可以看出,當(dāng)同時包含正負(fù)異常時,傳統(tǒng)的邊界探測器THDR和Theta圖雖然能很好地識別地質(zhì)體的邊界,但是引入了一些額外的錯誤邊界.然而,本文提出的均衡探測器能很好地識別異常體的邊界,沒有引入錯誤的邊界信息,能更好地完成邊界識別工作.其中,NEDT中的p值為0.1;NEEDT中的p值為0.08.

    圖2 模型體平面圖(a)和三維圖(b)Fig.2 Plan view (a) and 3D view (b) of synthetic model

    圖3 合成重力異常圖,棱柱體的剩余密度為200 kg·m-3Fig.3 Synthetic gravity anomaly model. The contrasted densities of prisms are 200 kg·m-3

    為了驗證方法的穩(wěn)定性,對圖6中模型二的各個重力梯度值增加5%的高斯白噪聲,圖7a—7c為加入噪聲后的重力梯度值.圖7d—7i為噪聲數(shù)據(jù)的邊界檢測結(jié)果.其中,NEDT中的p值為0.1,NEEDT中的p值為0.0.001.可以看出,EEDT相對于EDT受噪聲的影響較大,主要是由于EEDT利用了位場數(shù)據(jù)的三階導(dǎo)數(shù),放大了噪聲的影響;由于噪聲的干擾,THDR已經(jīng)無法給出地質(zhì)體的邊界,而NEEDT和Theta圖法受到噪聲的影響較大.相對而言,NEDT抗噪能力強(qiáng),更加穩(wěn)定.

    圖5 合成的重力異常,棱柱體1和3的剩余密度為 200 kg·m-3,棱柱體2和4的剩余密度為-200 kg·m-3Fig.5 Synthetic gravity anomaly model. The contrasted densities of prisms 1 and 3 are 200 kg·m-3, the contrasted densities of prisms 2 and 4 are -200 kg·m-3

    圖6 同圖4,但為模型二的重力梯度和不同方法的邊界識別結(jié)果Fig.6 Same as Fig.4, but for gravity gradients and the edges detected by different methods of model two

    圖7 同圖4,但為加入噪聲后的模型二的重力梯度和不同方法的邊界識別結(jié)果Fig.7 Same as Fig.4, but for gravity gradients and the edges detected by different methods of model two add 5% Gaussian noise

    圖8 同圖4,但為加拿大圣喬治灣地區(qū)全張量重力梯度值邊界檢測結(jié)果Fig.8 Same as Fig. 4, but for edges of full tensor gravity gradients anomalies in St. Georges Bay, Canada

    圖9 朱日和地區(qū)磁異常邊界檢測結(jié)果 (a)化極后的磁異常數(shù)據(jù);(b)異常的EDT結(jié)果;(c)異常的EEDT結(jié)果;(d)異常的NEDT結(jié)果; (e)異常的NEEDT結(jié)果;(f)異常的THDR結(jié)果;(g)異常的Theta圖結(jié)果.Fig.9 Edge results of magnetic anomalies in the Zhurihe area (a) Reduction to pole of magnetic data; (b) EDT of the data in (a); (c) EEDT of the data in (a); (d) NEDT of the data in (a); (e) NEEDT of the data in (a); (f) THDR of the data in (a); (g) Theta map of the data in (a).

    利用本文提出的方法對磁異常數(shù)據(jù)進(jìn)行處理前,需要對磁異常數(shù)據(jù)進(jìn)行化極處理,因為磁數(shù)據(jù)及其導(dǎo)數(shù)均受磁化方向的干擾.因此采用化極后的數(shù)據(jù)進(jìn)行解釋所獲得的結(jié)果將更加準(zhǔn)確(Li, 2006).

    6 應(yīng)用到實際數(shù)據(jù)

    為了驗證方法在實際中的應(yīng)用效果,分別對加拿大圣喬治灣的全張量重力梯度張量數(shù)據(jù)和中國朱日和地區(qū)的磁異常數(shù)據(jù)進(jìn)行邊界檢測.

    圖8a—8c為Air-FTG全張量重力梯度儀在圣喬治灣測得的重力梯度張量數(shù)據(jù),飛行線距為500 m,聯(lián)絡(luò)線距為5000 m,飛行的高度為100 m;圖8d—8i為識別的邊界結(jié)果.可以看出,非均衡邊界探測器EDT和EEDT的識別能力較差,識別的邊界結(jié)果不夠清晰,且不能同時顯示不同振幅大小的異常的邊界.常用的均衡探測器THDR的顯示能力較差,不能顯示地質(zhì)體邊界.Theta圖法可以很好地顯示地質(zhì)體的邊界,但是識別出來的邊界并不十分清晰.而NEDT和NEEDT識別出來的邊界相對于Theta圖的結(jié)果更加清晰、準(zhǔn)確,尤其是NEEDT識別出來的邊界更清晰準(zhǔn)確.其中,NEDT中的p值為0.005,NEEDT中的p值為0.001.

    圖9a為朱日和地區(qū)化極后的磁異常數(shù)據(jù),勘探區(qū)域面積為73 km×117 km,采樣間隔為20 m,飛行高度為1500 m,飛行線距為500 m.勘探區(qū)域為新元古代超層序,主要由大陸沉積物組成,除一些富含鐵的沙堰堤壩外基本不含任何鐵磁性物質(zhì).磁異常數(shù)據(jù)主要由SE-NW向的堤壩產(chǎn)生的近線性異常和鐵異常構(gòu)成.圖中的黑色點線為已知的構(gòu)造.相應(yīng)的梯度值通過數(shù)值計算得到(Minkus and Hinojosa, 2001).圖9b—9g為相應(yīng)的邊界檢測結(jié)果.可以看出,NEDT和NEEDT的邊界識別效果明顯好于其他方法的,能更加清晰地顯示邊界信息.其中,NEDT和NEEDT中的p值為0.005.

    7 結(jié)論

    本文針對位場梯度張量數(shù)據(jù)包含多個信息成分的優(yōu)勢,提出方向總水平導(dǎo)數(shù)和加強(qiáng)方向總水平導(dǎo)數(shù).根據(jù)方向總水平導(dǎo)數(shù)和加強(qiáng)方向總水平導(dǎo)數(shù)建立了兩個新的邊界探測器,并對其進(jìn)行歸一化處理,以便能同時檢測不同振幅大小異常的邊界.當(dāng)異常數(shù)據(jù)同時包含正負(fù)異常時,本文提出的新探測器能很好地避免產(chǎn)生錯誤的邊界信息,而傳統(tǒng)的THDR和Theta圖將產(chǎn)生額外的錯誤邊界信息.通過模型試驗和實際數(shù)據(jù),證明了基于方向總水平導(dǎo)數(shù)和加強(qiáng)方向總水平導(dǎo)數(shù)定義的邊界探測器能很好地識別地質(zhì)體的邊界,且識別出來的邊界更加清晰、準(zhǔn)確,避免了錯誤邊界信息的引入.

    致謝 感謝加拿大自然資源部(Natural Resources Canada)提供圣喬治灣(St. Georges Bay)的Air-FTG全張量重力梯度測量數(shù)據(jù),并允許發(fā)表數(shù)據(jù)處理結(jié)果.感謝為評審本文所付出努力的專家.

    Beiki M. 2010. Analytic signals of gravity gradient tensor and their application to estimate source location.Geophysics, 75(6): I59-I74.

    Cooper G R J, Cowan D R. 2004. Filtering using variable order vertical derivatives.Computer&Geosciences, 30(5): 455-459.

    Cooper G R J, Cowan D R. 2006. Enhancing potential field data using filters based on the local phase.Computer&Geosciences, 32(10): 1585-1591.

    Cordell L. 1979. Gravimetric expression of Graben faulting in Santa Fe country and the Espanola Basin, New Mexico. New Mexico Geol. Soc. Guidebook, 30thField Conf., 1979: 59-64.

    Cordell L, Grauch V J S. 1985. Mapping basement magnetization zones from aeromagnetic data in the San Juan basin, New Mexico. ∥Hinzc W J ed. The Utility of Regional Gravity and Magnetic Anomaly. Society of Exploration Geophysics, 181-197.

    Evjen H M. 1936. The place of the vertical gradient in gravitational interpretations.Geophysics, 1(1): 127-136.

    Forsberg R. 1984. A study of terrain reductions, density anomalies and geophysical inversion methods in gravity field modelling. Dept. of Geodetic Science and Surveying, Report 355, OhioState University.

    Hsu S H, Sibuet J C, Shyu C T. 1996. High-resolution detection of geologic boundaries from potential-field anomalies: An enhanced analytic signal technique.Geophysics, 61(2): 373-386. Li X. 2006. Understanding 3D analytic signal amplitude.Geophysics, 71(2): L13-L16. Ma G Q, Huang D N, Yu P, et al. 2012. Application of improved balancing filters to edge identification of potential field data.ChineseJ.Geophys. (in Chinese), 55(12): 4288-4295, doi: 10.6038/j.issn.0001-5733.2012.12.040.

    Marson I, Klingele E E. 1993. Advantages of using the vertical gradient of gravity for 3-D interpretation.Geophysics, 58(11): 1588-1595. Mikhailov V, Pajot G, Diament M, et al. 2007. Tensor deconvolution: A method to locate equivalent sources from full tensor gravity data.Geophysics, 72(5): 161-169.

    Miller H G, Singh V. 1994. Potential field tilt—a new concept for location of potential field sources.JournalofAppliedGeophysics, 32(2-3): 213-217. Minkus K L, Hinojosa J H. 2001. The complete gravity gradient tensor derived from the vertical component of gravity: a Fourier transform technique.JournalofAppliedGeophysics, 46(3): 159-174.

    Nabighian M N. 1984. Toward a three-dimensional automatic interpretation of potential field data via generalized Hilbert transforms: Fundamental relations.Geophysics, 49(6): 780-786.

    Oru? B, Keskinsezer A. 2008. Structural setting of the northeastern

    Biga Peninsula (Turkey) from tilt derivatives of gravity gradient tensors and magnitude of horizontal gravity components.PureandAppliedGeophysics, 165(9-10): 1913-1927. Oru? B, Sert?elik I, Kafadar ?, et al. 2013. Structural interpretation of the Erzurum Basin, eastern Turkey, using curvature gravity gradient tensor and gravity inversion of basement relief.JournalofAppliedGeophysics, 88: 105-113. Rajagopalan S, Milligan P. 1994. Image enhancement of aeromagnetic data using automatic gain control.ExplorationGeophysics, 25(4): 173-178.

    Sertcelik I, Kafadar O. 2012. Application of edge detection to potential field data using eigenvalue analysis of structure tensor.JournalofAppliedGeophysics, 84: 86-94.

    Verduzco B, Fairhead J D, Green C M, et al. 2004. New insights into magnetic derivatives for structural mapping.TheLeadingEdge, 23(2): 116-119.

    Wijns C, Perez C, Kowalczyk P. 2005. Theta map: edge detection in magnetic data.Geophysics, 70(4): L39-L43.

    Zhou W N, Du X J, Li J Y. 2013. The limitation of curvature gravity gradient tensor for edge detection and a method for overcoming it.JournalofAppliedGeophysics, 98: 237-242.

    附中文參考文獻(xiàn)

    馬國慶, 黃大年, 于平等. 2012. 改進(jìn)的均衡濾波器在位場數(shù)據(jù)邊界識別中的應(yīng)用. 地球物理學(xué)報, 55(12): 4288-4295, doi: 10.6038/j.issn.0001-5733.2012.12.040.

    (本文編輯 胡素芳)

    Using enhanced directional total horizontal derivatives to detect the edges of potential-field full tensor data

    YUAN Yuan1,2,3, HUANG Da-Nian1*, YU Qing-Lu4

    1CollegeofGeoExplorationScienceandTechnology,JilinUniversity,Changchun130026,China2SecondInstituteofOceanography,StateOceanicAdministration,Hangzhou310012,China3KeyLaboratoryofSubmarineGeoscience,StateOceanicAdministration,Hangzhou310012,China4SinopecGeophysicalResearchInstitute,Nanjing211103,China

    With the ongoing development of the full potential field gradient techniques, more and more gradient tensor data have been widely used in geophysical exploration for its large amount of information and containing higher frequency signals than potential field data. The full potential field gradient technique can simultaneously measure six gradient components. Each component has its own geophysical meaning. The high frequency gradient tensor data can be used to delineate small scale anomalies.Edge detection is required for interpretation of potential field data, and has been used in exploration technology for discovery of mineral resources. The main geological edges are fault lines and the borders of geological or rock bodies of different density, magnetic nature, etc. On account of the higher-frequency signals, their interpretation enables high-resolution and detailed investigation of geological structures. Development of new methods is required to enable interpretation of these data, especially for edge detection for small geological structures.In order to make use of multiple component information, we define directional total horizontal derivatives and enhanced directional total horizontal derivatives and use them to define new edge detectors. However, the new defined edge detectors can balance the detected edge signals of different amplitude anomalies. In order to balance the edge signal amplitude, we present a normalization method, where we divide the new edge detectors by the vertical potential field gradient data. Besides, we introduce a constant parameter in the denominator of the normalization method, which can effectively avoid bringing some additional false edges when real geological bodies contain positive and negative anomalies simultaneously.These methods have been tested with synthetic data to verify that the new methods can delineate the edges of different amplitude anomalies clearly. The results show that new defined directional total horizontal derivatives and enhanced directional total horizontal derivatives have higher resolution than directional analytic signal. The normalized edge detectors can display the edges of large and small amplitude anomalies simultaneously, and avoid introducing additional false edges. To further test the stability, we demonstrate the new edge detectors with the model data corrupted with 5% Gaussian noise. Finally, we apply these methods to real full gravity gradient tensor data in St. Georges Bay, Canada and magnetic anomalies in Zhurihe area, China, and get good edge results.

    Edge detection; Potential field gradient tensor; Enhanced directional total horizontal derivative; Normalization

    10.6038/cjg20150730.

    國家高技術(shù)研究發(fā)展計劃(863計劃)課題(2014AA06A613),國家重點基礎(chǔ)研究發(fā)展計劃項目(973項目)西南印度洋洋中脊熱液成礦過程與硫化物礦區(qū)預(yù)測的第五課題“硫化物礦區(qū)特征和找礦標(biāo)志”(2012CB417305),大洋十二五項目西南印度洋脊合同區(qū)多金屬硫化物資源評價聯(lián)合資助.

    袁園,男,1988年生,助理研究員,主要從事位場數(shù)據(jù)處理及解釋方面的研究.E-mail:yuanyuan12@mails.jlu.edu.cn

    *通訊作者 黃大年,男,1958年生,教授,博士生導(dǎo)師,主要從事移動平臺探測數(shù)據(jù)處理與解釋及一體化軟件平臺開發(fā). E-mail:dnhuang@jlu.edu.cn

    10.6038/cjg20150730

    P631

    2015-02-22,2015-06-16收修定稿

    袁園, 黃大年, 余青露. 2015. 利用加強(qiáng)水平方向總水平導(dǎo)數(shù)對位場全張量數(shù)據(jù)進(jìn)行邊界識別.地球物理學(xué)報,58(7):2556-2565,

    Yuan Y, Huang D N, Yu Q L. 2015. Using enhanced directional total horizontal derivatives to detect the edges of potential-field full tensor data.ChineseJ.Geophys. (in Chinese),58(7):2556-2565,doi:10.6038/cjg20150730.

    猜你喜歡
    方向水平
    張水平作品
    2022年組稿方向
    2022年組稿方向
    2021年組稿方向
    2021年組稿方向
    2021年組稿方向
    作家葛水平
    火花(2019年12期)2019-12-26 01:00:28
    加強(qiáng)上下聯(lián)動 提升人大履職水平
    老虎獻(xiàn)臀
    位置與方向
    哪个播放器可以免费观看大片| 亚洲精品美女久久久久99蜜臀 | 女性生殖器流出的白浆| 五月伊人婷婷丁香| 乱人伦中国视频| 晚上一个人看的免费电影| 桃花免费在线播放| 边亲边吃奶的免费视频| 深夜精品福利| a 毛片基地| 人人妻人人澡人人看| 人人妻人人澡人人爽人人夜夜| av在线老鸭窝| 999久久久国产精品视频| 热re99久久国产66热| 免费av中文字幕在线| 有码 亚洲区| 一本色道久久久久久精品综合| 日韩人妻精品一区2区三区| 久久久久久久久免费视频了| 欧美精品国产亚洲| √禁漫天堂资源中文www| 男女高潮啪啪啪动态图| 青春草亚洲视频在线观看| 国产人伦9x9x在线观看 | 成年人午夜在线观看视频| 色婷婷久久久亚洲欧美| 久久国内精品自在自线图片| 国产片特级美女逼逼视频| 下体分泌物呈黄色| 久久人人爽人人片av| 建设人人有责人人尽责人人享有的| 蜜桃国产av成人99| 三上悠亚av全集在线观看| 纯流量卡能插随身wifi吗| 午夜日韩欧美国产| xxx大片免费视频| 免费观看a级毛片全部| 国产成人精品久久二区二区91 | 好男人视频免费观看在线| 啦啦啦在线观看免费高清www| 香蕉国产在线看| 久久99一区二区三区| 免费大片黄手机在线观看| 一级爰片在线观看| 国产无遮挡羞羞视频在线观看| 日韩成人av中文字幕在线观看| 日韩一卡2卡3卡4卡2021年| 婷婷色麻豆天堂久久| 亚洲成人av在线免费| 久久久精品免费免费高清| 国产欧美日韩一区二区三区在线| av女优亚洲男人天堂| 国精品久久久久久国模美| 一级毛片 在线播放| 免费高清在线观看视频在线观看| 久久久久久久久久久免费av| 中文字幕亚洲精品专区| 久久韩国三级中文字幕| 国产激情久久老熟女| av福利片在线| 母亲3免费完整高清在线观看 | 老鸭窝网址在线观看| 人人妻人人澡人人爽人人夜夜| 一区福利在线观看| 在线亚洲精品国产二区图片欧美| 亚洲欧美精品综合一区二区三区 | 女性被躁到高潮视频| av国产精品久久久久影院| 国产精品欧美亚洲77777| 亚洲国产av新网站| 欧美日韩精品网址| 亚洲精品国产一区二区精华液| 中文字幕人妻熟女乱码| 久久精品久久久久久噜噜老黄| 大香蕉久久成人网| www.av在线官网国产| 欧美激情 高清一区二区三区| 如何舔出高潮| 成人二区视频| 精品亚洲成国产av| 久久精品夜色国产| 久久免费观看电影| 亚洲激情五月婷婷啪啪| 如日韩欧美国产精品一区二区三区| 国产 一区精品| 成人手机av| 大码成人一级视频| www.熟女人妻精品国产| 国产乱人偷精品视频| 如何舔出高潮| 激情视频va一区二区三区| 中文字幕精品免费在线观看视频| 久久97久久精品| 欧美精品av麻豆av| 色吧在线观看| 69精品国产乱码久久久| 亚洲精品视频女| av.在线天堂| 午夜福利乱码中文字幕| 一级爰片在线观看| 欧美日韩亚洲高清精品| 一级黄片播放器| 麻豆乱淫一区二区| 超碰成人久久| 我要看黄色一级片免费的| 丰满乱子伦码专区| 七月丁香在线播放| 哪个播放器可以免费观看大片| 大话2 男鬼变身卡| 女性生殖器流出的白浆| 久久久欧美国产精品| 美女中出高潮动态图| 乱人伦中国视频| 纯流量卡能插随身wifi吗| 另类精品久久| 精品国产超薄肉色丝袜足j| 美女国产高潮福利片在线看| 又黄又粗又硬又大视频| 性色avwww在线观看| 少妇猛男粗大的猛烈进出视频| 午夜福利视频精品| 欧美人与性动交α欧美软件| 亚洲一级一片aⅴ在线观看| 国产 一区精品| 欧美bdsm另类| 18在线观看网站| 又大又黄又爽视频免费| 两性夫妻黄色片| 亚洲色图 男人天堂 中文字幕| 国产免费现黄频在线看| 巨乳人妻的诱惑在线观看| 18禁观看日本| 久久久a久久爽久久v久久| 九九爱精品视频在线观看| 国产xxxxx性猛交| 伦理电影大哥的女人| 欧美精品一区二区大全| 成人国产麻豆网| 亚洲欧洲精品一区二区精品久久久 | 日韩欧美精品免费久久| 日韩一区二区三区影片| 国产乱人偷精品视频| 国产亚洲精品第一综合不卡| 青春草视频在线免费观看| 国产日韩欧美在线精品| 国产精品香港三级国产av潘金莲 | 视频区图区小说| av福利片在线| 一区在线观看完整版| 男人操女人黄网站| 国产精品国产三级专区第一集| 美女大奶头黄色视频| 欧美精品一区二区大全| 精品福利永久在线观看| 国产福利在线免费观看视频| 高清黄色对白视频在线免费看| 欧美日韩一区二区视频在线观看视频在线| 女性生殖器流出的白浆| 欧美激情极品国产一区二区三区| 亚洲国产日韩一区二区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 蜜桃在线观看..| 69精品国产乱码久久久| 视频在线观看一区二区三区| 国产一区二区三区av在线| 欧美日韩视频精品一区| 国产熟女午夜一区二区三区| 亚洲熟女精品中文字幕| 看非洲黑人一级黄片| 国产精品无大码| 中文欧美无线码| 五月伊人婷婷丁香| 性少妇av在线| 女人精品久久久久毛片| 成人国产av品久久久| 欧美 亚洲 国产 日韩一| av电影中文网址| 日韩免费高清中文字幕av| 免费日韩欧美在线观看| 777久久人妻少妇嫩草av网站| 下体分泌物呈黄色| 久久久久国产精品人妻一区二区| 日韩 亚洲 欧美在线| 欧美成人午夜精品| 日韩中文字幕欧美一区二区 | 欧美精品人与动牲交sv欧美| 久久久亚洲精品成人影院| 日产精品乱码卡一卡2卡三| 最新的欧美精品一区二区| 亚洲人成电影观看| 两性夫妻黄色片| 婷婷色综合大香蕉| 日韩伦理黄色片| 日韩制服骚丝袜av| 国产一区亚洲一区在线观看| 久久午夜综合久久蜜桃| 国产精品欧美亚洲77777| 国产男女内射视频| 国产亚洲午夜精品一区二区久久| 黑人猛操日本美女一级片| 亚洲国产精品一区三区| 国产av国产精品国产| 亚洲国产欧美网| 亚洲精品日本国产第一区| 亚洲精品第二区| 欧美成人精品欧美一级黄| 国产日韩欧美亚洲二区| 美女中出高潮动态图| 18禁国产床啪视频网站| 美女高潮到喷水免费观看| 欧美xxⅹ黑人| 超色免费av| 日韩伦理黄色片| 亚洲成国产人片在线观看| 一边亲一边摸免费视频| 国产精品一区二区在线不卡| 性色avwww在线观看| 亚洲精品成人av观看孕妇| 伊人久久大香线蕉亚洲五| tube8黄色片| 人妻 亚洲 视频| 欧美日本中文国产一区发布| 男女下面插进去视频免费观看| 免费在线观看黄色视频的| 在线观看国产h片| 卡戴珊不雅视频在线播放| av国产久精品久网站免费入址| 一级a爱视频在线免费观看| 国产精品国产三级国产专区5o| 精品亚洲乱码少妇综合久久| 人妻 亚洲 视频| 国产一区二区激情短视频 | 免费在线观看完整版高清| 亚洲人成77777在线视频| 亚洲美女黄色视频免费看| 国产成人免费观看mmmm| 天堂8中文在线网| 久久精品国产综合久久久| 日韩中文字幕欧美一区二区 | 日韩不卡一区二区三区视频在线| 欧美日韩亚洲国产一区二区在线观看 | 国产成人欧美| 新久久久久国产一级毛片| 久久久久久久亚洲中文字幕| 亚洲,欧美,日韩| 最近中文字幕2019免费版| 91午夜精品亚洲一区二区三区| 日韩伦理黄色片| 18禁裸乳无遮挡动漫免费视频| 有码 亚洲区| 91精品国产国语对白视频| av卡一久久| 一二三四在线观看免费中文在| 久热久热在线精品观看| 免费少妇av软件| 成年动漫av网址| 国产1区2区3区精品| 激情五月婷婷亚洲| 丰满少妇做爰视频| 亚洲精品成人av观看孕妇| 久久精品国产鲁丝片午夜精品| 蜜桃在线观看..| 国产av精品麻豆| 日本色播在线视频| 搡老乐熟女国产| av国产久精品久网站免费入址| 亚洲国产成人一精品久久久| 久久亚洲国产成人精品v| 欧美精品人与动牲交sv欧美| 精品一区二区免费观看| 9色porny在线观看| 97人妻天天添夜夜摸| av线在线观看网站| 80岁老熟妇乱子伦牲交| 成人毛片a级毛片在线播放| 成人国产麻豆网| 欧美中文综合在线视频| √禁漫天堂资源中文www| 好男人视频免费观看在线| 天堂俺去俺来也www色官网| 欧美成人精品欧美一级黄| 成年动漫av网址| 国产免费福利视频在线观看| 亚洲综合色惰| 美女xxoo啪啪120秒动态图| 一边摸一边做爽爽视频免费| 久久久精品区二区三区| 老女人水多毛片| 国产av精品麻豆| 日本wwww免费看| 午夜福利视频精品| 九九爱精品视频在线观看| 亚洲成色77777| 免费播放大片免费观看视频在线观看| 亚洲,一卡二卡三卡| 日本91视频免费播放| 日韩一本色道免费dvd| 久久久久久久精品精品| 女人被躁到高潮嗷嗷叫费观| 在线精品无人区一区二区三| 观看av在线不卡| 色播在线永久视频| 国产日韩欧美亚洲二区| 欧美xxⅹ黑人| 精品午夜福利在线看| 国产高清国产精品国产三级| 97精品久久久久久久久久精品| 成人午夜精彩视频在线观看| 久久精品人人爽人人爽视色| 在线观看免费日韩欧美大片| 2018国产大陆天天弄谢| 春色校园在线视频观看| 狠狠精品人妻久久久久久综合| 国产精品国产av在线观看| 男女免费视频国产| 亚洲国产av影院在线观看| 天美传媒精品一区二区| 妹子高潮喷水视频| 国产精品久久久av美女十八| 成人亚洲精品一区在线观看| 男人操女人黄网站| 最近中文字幕高清免费大全6| 超碰97精品在线观看| 男人添女人高潮全过程视频| 夫妻性生交免费视频一级片| 最新的欧美精品一区二区| 爱豆传媒免费全集在线观看| 一本久久精品| 国产爽快片一区二区三区| 成人午夜精彩视频在线观看| 母亲3免费完整高清在线观看 | 日日撸夜夜添| 国产 精品1| 高清视频免费观看一区二区| 亚洲欧洲国产日韩| 欧美 亚洲 国产 日韩一| 久久久精品免费免费高清| 青草久久国产| 午夜日韩欧美国产| 欧美97在线视频| 国产一区二区三区av在线| 啦啦啦中文免费视频观看日本| 亚洲第一av免费看| 男女免费视频国产| 成人漫画全彩无遮挡| 国产精品久久久久久精品古装| 人妻 亚洲 视频| 男女免费视频国产| 一级片'在线观看视频| 亚洲四区av| 国产色婷婷99| 精品国产乱码久久久久久男人| 美女国产视频在线观看| 咕卡用的链子| 97在线人人人人妻| 国产成人精品在线电影| 久久这里只有精品19| 亚洲色图综合在线观看| 激情视频va一区二区三区| 久久亚洲国产成人精品v| 亚洲一区二区三区欧美精品| 亚洲欧洲日产国产| 又黄又粗又硬又大视频| 亚洲少妇的诱惑av| 99精国产麻豆久久婷婷| 极品人妻少妇av视频| 毛片一级片免费看久久久久| 久久久久久人妻| 大香蕉久久网| 日韩中字成人| 高清不卡的av网站| 男人舔女人的私密视频| 五月开心婷婷网| 搡女人真爽免费视频火全软件| 啦啦啦在线免费观看视频4| 婷婷成人精品国产| 在线精品无人区一区二区三| 欧美人与性动交α欧美精品济南到 | 国产 一区精品| 人人妻人人添人人爽欧美一区卜| 亚洲精品中文字幕在线视频| 一二三四在线观看免费中文在| 在线观看www视频免费| 欧美人与善性xxx| 国产 一区精品| 亚洲色图综合在线观看| 天天躁夜夜躁狠狠躁躁| 成人亚洲精品一区在线观看| 亚洲欧洲国产日韩| 国语对白做爰xxxⅹ性视频网站| 欧美成人午夜精品| 免费少妇av软件| 日本免费在线观看一区| 亚洲精品日韩在线中文字幕| 成人漫画全彩无遮挡| 999精品在线视频| 黄色配什么色好看| 亚洲av.av天堂| 亚洲一区中文字幕在线| 赤兔流量卡办理| 丝袜人妻中文字幕| 亚洲情色 制服丝袜| 老汉色∧v一级毛片| 黑丝袜美女国产一区| 中国三级夫妇交换| 下体分泌物呈黄色| 免费在线观看完整版高清| 最黄视频免费看| 春色校园在线视频观看| 黄网站色视频无遮挡免费观看| tube8黄色片| 美女大奶头黄色视频| 成人18禁高潮啪啪吃奶动态图| 亚洲伊人色综图| 久久久久久久久久久久大奶| 国产精品 欧美亚洲| 韩国av在线不卡| www.av在线官网国产| 精品久久久久久电影网| 中文乱码字字幕精品一区二区三区| 久久婷婷青草| 五月伊人婷婷丁香| 人妻人人澡人人爽人人| h视频一区二区三区| 国产精品一国产av| 咕卡用的链子| 日韩欧美一区视频在线观看| 伊人亚洲综合成人网| 日本av手机在线免费观看| 亚洲第一av免费看| 视频区图区小说| 一区二区三区激情视频| 久久久国产精品麻豆| 熟女少妇亚洲综合色aaa.| 久久精品久久久久久久性| 黑丝袜美女国产一区| 国产视频首页在线观看| 如何舔出高潮| 高清黄色对白视频在线免费看| 久久久久久久久免费视频了| 99久久人妻综合| 久久免费观看电影| 久久 成人 亚洲| 亚洲三级黄色毛片| 最黄视频免费看| 国产精品99久久99久久久不卡 | 欧美精品高潮呻吟av久久| 一级片免费观看大全| xxxhd国产人妻xxx| 一个人免费看片子| 2022亚洲国产成人精品| 欧美 亚洲 国产 日韩一| 午夜福利一区二区在线看| 中文字幕人妻熟女乱码| 免费观看无遮挡的男女| 久久久国产欧美日韩av| 人成视频在线观看免费观看| 久久精品久久精品一区二区三区| 男人添女人高潮全过程视频| 侵犯人妻中文字幕一二三四区| av有码第一页| 人人澡人人妻人| 久久久精品免费免费高清| 国产一区亚洲一区在线观看| 国产一级毛片在线| 国产女主播在线喷水免费视频网站| 日本av手机在线免费观看| 天天躁狠狠躁夜夜躁狠狠躁| 一区二区三区精品91| 日韩免费高清中文字幕av| 午夜福利,免费看| 亚洲一级一片aⅴ在线观看| 黄网站色视频无遮挡免费观看| 国产高清不卡午夜福利| 精品一区在线观看国产| av线在线观看网站| 这个男人来自地球电影免费观看 | 国产亚洲欧美精品永久| 如何舔出高潮| 国产成人精品久久二区二区91 | 这个男人来自地球电影免费观看 | av在线播放精品| 成年女人在线观看亚洲视频| 99香蕉大伊视频| 久久久久国产精品人妻一区二区| 丝袜美腿诱惑在线| 亚洲av电影在线观看一区二区三区| 亚洲av在线观看美女高潮| 两性夫妻黄色片| 丝袜脚勾引网站| 国产高清不卡午夜福利| 久久 成人 亚洲| 尾随美女入室| 精品久久蜜臀av无| 中文字幕制服av| 边亲边吃奶的免费视频| 中文字幕av电影在线播放| 精品人妻在线不人妻| 精品国产一区二区久久| 麻豆乱淫一区二区| 欧美xxⅹ黑人| 毛片一级片免费看久久久久| 熟女av电影| 又粗又硬又长又爽又黄的视频| 熟妇人妻不卡中文字幕| 成人手机av| 九色亚洲精品在线播放| 欧美激情高清一区二区三区 | 久久ye,这里只有精品| 日本午夜av视频| 王馨瑶露胸无遮挡在线观看| 亚洲男人天堂网一区| 亚洲精品国产一区二区精华液| 亚洲国产精品国产精品| 欧美bdsm另类| 亚洲三区欧美一区| 亚洲第一区二区三区不卡| 最近中文字幕高清免费大全6| 久热这里只有精品99| 中文字幕色久视频| 涩涩av久久男人的天堂| 精品一区二区三区四区五区乱码 | 国产精品.久久久| av国产精品久久久久影院| 国产成人一区二区在线| 搡老乐熟女国产| 国产又色又爽无遮挡免| 国产亚洲午夜精品一区二区久久| 国产免费又黄又爽又色| 国产精品亚洲av一区麻豆 | 精品人妻一区二区三区麻豆| 日韩制服骚丝袜av| 精品人妻一区二区三区麻豆| 亚洲av男天堂| 国产日韩欧美亚洲二区| 日本wwww免费看| 亚洲精品国产av蜜桃| 青草久久国产| 一二三四在线观看免费中文在| 欧美xxⅹ黑人| 国产免费一区二区三区四区乱码| 亚洲,欧美,日韩| 男女国产视频网站| 国产色婷婷99| 日日撸夜夜添| 男人舔女人的私密视频| 中文欧美无线码| 狂野欧美激情性bbbbbb| 亚洲第一区二区三区不卡| 精品卡一卡二卡四卡免费| 久久鲁丝午夜福利片| 亚洲精品国产色婷婷电影| 美女中出高潮动态图| 亚洲一级一片aⅴ在线观看| 欧美成人午夜免费资源| 欧美日韩视频精品一区| 人妻系列 视频| 国产成人精品久久二区二区91 | 观看av在线不卡| 麻豆精品久久久久久蜜桃| 91午夜精品亚洲一区二区三区| 国产精品嫩草影院av在线观看| 美女脱内裤让男人舔精品视频| 91精品国产国语对白视频| 91精品三级在线观看| 亚洲av电影在线观看一区二区三区| 久久这里有精品视频免费| 欧美精品国产亚洲| 亚洲国产毛片av蜜桃av| 99热网站在线观看| 国产日韩欧美亚洲二区| xxx大片免费视频| 成年女人在线观看亚洲视频| av一本久久久久| 午夜福利乱码中文字幕| 高清av免费在线| 日韩视频在线欧美| 亚洲av福利一区| 国产免费又黄又爽又色| 中文字幕色久视频| 亚洲欧美色中文字幕在线| 一本—道久久a久久精品蜜桃钙片| 久久精品国产鲁丝片午夜精品| 国产精品秋霞免费鲁丝片| 熟妇人妻不卡中文字幕| 免费在线观看黄色视频的| 亚洲国产成人一精品久久久| 婷婷色综合大香蕉| 热re99久久精品国产66热6| 日本色播在线视频| 精品99又大又爽又粗少妇毛片| 日韩伦理黄色片| 18禁观看日本| 国产精品秋霞免费鲁丝片| 久久人人97超碰香蕉20202| 国产极品粉嫩免费观看在线| 亚洲精品中文字幕在线视频| 欧美日韩视频精品一区| 成人二区视频| 中文字幕人妻熟女乱码| 中文字幕亚洲精品专区| 美女主播在线视频| 国产欧美日韩一区二区三区在线| 婷婷色综合大香蕉| 伊人久久大香线蕉亚洲五| 亚洲国产av新网站| 精品人妻一区二区三区麻豆| 日日撸夜夜添| 亚洲国产毛片av蜜桃av| 精品久久蜜臀av无| 国产免费又黄又爽又色| 欧美精品国产亚洲| 一区二区日韩欧美中文字幕| 日本猛色少妇xxxxx猛交久久|