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

    一種自動判定保留的奇異值個數(shù)的地震隨機(jī)噪聲壓制算法

    2022-06-11 01:25:00朱躍飛曹靜杰殷晗鈞
    石油地球物理勘探 2022年3期
    關(guān)鍵詞:同相軸個數(shù)信噪比

    朱躍飛 曹靜杰* 殷晗鈞

    (①自然資源部京津冀城市群地下空間智能探測與裝備重點(diǎn)實(shí)驗(yàn)室,河北石家莊 050031; ②河北省戰(zhàn)略性關(guān)鍵礦產(chǎn)資源重點(diǎn)實(shí)驗(yàn)室,河北石家莊 050031)

    0 引言

    受人類活動、儀器、環(huán)境、天氣等多種因素影響,野外采集的地震數(shù)據(jù)往往含有各種噪聲,嚴(yán)重影響速度分析和靜校正、速度建模及偏移成像等處理的效果。因此,消除噪聲以獲取高信噪比地震數(shù)據(jù)一直是地震勘探面臨的難題[1]。地震數(shù)據(jù)去噪一般依賴信號和噪聲在頻率、統(tǒng)計(jì)規(guī)律、振幅等方面的差異分離信號和噪聲。地震噪聲分為隨機(jī)噪聲和相干噪聲。地震隨機(jī)噪聲的消除方法眾多,大體上分為濾波類方法、基于變換的方法、降秩方法和深度學(xué)習(xí)方法等。

    濾波類方法基于地震數(shù)據(jù)時間域分布特點(diǎn)構(gòu)建濾波函數(shù)去除噪聲,主要方法有中值濾波[2-3]、各向異性擴(kuò)散濾波[4]等。基于變換的方法假設(shè)地震數(shù)據(jù)經(jīng)過某個變換后的系數(shù)具有稀疏特征,選取較大的系數(shù),通過閾值運(yùn)算去掉小的系數(shù),最后反變換到時間域?qū)崿F(xiàn)去噪[5],常用的變換有傅里葉變換[6]、Radon變換[7]、Wavelet變換[8-10]、S變換[11]、曲波變換[12]等。

    深度學(xué)習(xí)去噪方法是目前的研究熱點(diǎn),基本原理是利用大量的樣本數(shù)據(jù)的特征,通過多層卷積的方式提取數(shù)據(jù)的時域特征,然后采用深度學(xué)習(xí)的非線性逼近能力調(diào)整網(wǎng)絡(luò)參數(shù),從而建立一個復(fù)雜的去噪模型實(shí)現(xiàn)去噪。目前卷積神經(jīng)網(wǎng)絡(luò)[13]、殘差學(xué)習(xí)[14-15]、生成對抗網(wǎng)絡(luò)[16]、降噪自編碼[17]等深度學(xué)習(xí)網(wǎng)絡(luò)被用于地震數(shù)據(jù)去噪。深度學(xué)習(xí)方法需對不同的數(shù)據(jù)大量訓(xùn)練,因此計(jì)算量大。

    多道奇異譜分析(MSSA)是一種基于奇異值分解的降秩去噪方法,通過奇異值分解將原始數(shù)據(jù)分解為信號子空間和噪聲子空間,然后將噪聲子空間的能量置為零(截?cái)?,再通過反變換去噪[18]。MSSA用于多道時間序列分析,是單道奇異譜分析(SSA)的推廣[19-20]。Read[21]率先將SSA拓展到多變量MSSA方法研究,基于線性同相軸的假設(shè),利用相鄰地震道的頻譜相似性與可預(yù)測性組成低秩的Hankel矩陣[22],噪聲破壞了數(shù)據(jù)頻率切片Hankel矩陣的低秩結(jié)構(gòu)[23],常用截?cái)嗥娈愔捣纸夥椒ń鉀Q低秩近似問題。在地震信號處理領(lǐng)域,MSSA和Cadzow濾波是等效的但卻來自不同的領(lǐng)域[24],即Cadzow濾波源于信號和圖像去噪,MSSA則源于分析由動力系統(tǒng)引起的時間序列,本文采用MSSA的名稱表示這類方法。Oropeza等[25]利用MSSA同時對疊前三維數(shù)據(jù)去噪和重建,數(shù)值實(shí)驗(yàn)表明無法完全消除隨機(jī)噪聲,其去噪效果有很大的提升空間。Huang等[26]將阻尼算子引入傳統(tǒng)MSSA中,提出了阻尼多道奇異譜分析(DMSSA)算法。通過融合軟閾值移動平均算子和阻尼算子的優(yōu)點(diǎn),Oboue等[27]利用魯棒阻尼降秩方法提高地震數(shù)據(jù)的信噪比。阻尼降秩方法已成為一種有效的去噪方法,可以從含噪和不完備的觀測數(shù)據(jù)中恢復(fù)有效信號。

    對于海量地震數(shù)據(jù)來說,基于降秩的方法需要將地震數(shù)據(jù)分成不同的塊,然而每個塊對應(yīng)的奇異值個數(shù)不同,目前需要人工估計(jì)每個塊的有效奇異值個數(shù),計(jì)算效率低,無法實(shí)現(xiàn)產(chǎn)業(yè)化。為此,本文利用Akaike信息準(zhǔn)則自動地確定地震信號的奇異值個數(shù),然后基于DMSSA方法去噪。首先介紹了MSSA方法的去噪原理,然后給出確定有效奇異值個數(shù)的Akaike信息準(zhǔn)則和經(jīng)驗(yàn)方法,經(jīng)驗(yàn)方法可以驗(yàn)證Akaike信息準(zhǔn)則的有效性。模擬和實(shí)際數(shù)據(jù)實(shí)驗(yàn)表明,Akaike信息準(zhǔn)則能夠自動確定有效奇異值個數(shù),避免了人工操作,有利于實(shí)現(xiàn)產(chǎn)業(yè)化。

    1 算法原理

    1.1 DMSSA方法[18]

    假設(shè)一個含噪三維地震數(shù)據(jù)為D(x,y,t),其中x=(x1,…,xm)、y=(y1,…,yn)表示空間坐標(biāo),t=(t1,…,ts)表示時間坐標(biāo),m、n為道數(shù),s為采樣點(diǎn)數(shù)。根據(jù)DMSSA理論,使用以下步驟去噪。

    (1)通過離散傅里葉變換將D(x,y,t)從時間域變換為頻率域數(shù)據(jù)F(x,y,ω),其中ω=(ω1,…,ωj)為離散的頻率序列,j為頻率切片個數(shù)。

    (2)在給定的頻率范圍內(nèi)將不同頻率切片數(shù)據(jù)排列成塊Hankel矩陣。當(dāng)頻率為ωi(i=1,…,j)時,有

    (1)

    首先,將F(x,y,ωi)的每一行構(gòu)造成Hankel矩陣

    (2)

    Rk表示由F(x,y,ωi)的第k行構(gòu)造的Hankel矩陣,大小為v×h,v=

    n/2

    +1,h=n-v+1,

    ·

    表示向下取整,其中2

    (3)

    H為(v×l)×(h×f)階塊Hankel矩陣,l=

    m/2

    +1,f=m-l+1,2

    (3)對H進(jìn)行奇異值分解,并且選擇和截?cái)嗥娈愔?,是MSSA類方法的關(guān)鍵。如果有效信號對應(yīng)的奇異值個數(shù)為N,則奇異值對角矩陣僅保留前N個奇異值,而其他所有奇異值均設(shè)置為零。對H進(jìn)行奇異值分解,得到

    (4)

    其中

    (5)

    式中:U為H的左奇異值向量組成的(v×l)×(v×l)階正交矩陣;VT為H的右奇異值向量組成的 (h×f)×(h×f)階正交矩陣;Σ為按奇異值遞減順序σ1≥σ2≥…≥σd組成的對角矩陣,非零奇異值的個數(shù)d等于H的秩。

    (4)基于截?cái)嗟钠娈愔涤?jì)算去噪結(jié)果。通過將Hankel矩陣反變換到頻率域,再通過離散傅里葉逆變換得到時間—空間域去噪地震數(shù)據(jù),即

    Σ1=diag(σ1,σ2,…,σN) 0≤N

    (6)

    (7)

    對所有的頻率域數(shù)據(jù)進(jìn)行上述操作即可得到去噪地震數(shù)據(jù)。

    若采集的地震數(shù)據(jù)中不含噪聲,則Σ僅包含與有效信號相關(guān)的非零σ。若采集的地震數(shù)據(jù)中含有噪聲,所有σ都會發(fā)生改變,非零σ的個數(shù)也將增加。原始MSSA方法僅保留了N個σ,對σ的大小并沒有影響,因此去噪結(jié)果有很大的改進(jìn)空間。Huang等[18]提出的DMSSA方法可以減小σ,因此去噪效果更好。

    DMSSA方法通過添加阻尼因子減弱由噪聲引起的σ增量,即

    (8)

    (9)

    式中:T為阻尼算子;I為單位矩陣;D為阻尼因子,其值越小,阻尼效果越強(qiáng),反之亦然。DMSSA去噪的本質(zhì)就是利用D對第N+1個σ放大或縮小,然后使用前N個σ與其求差,并對第N+1之后的σ置零,以達(dá)到壓制噪聲的目的。

    確定N是DMSSA去噪最關(guān)鍵的一步,將影響噪聲抑制效果和有效信號的保護(hù)程度。式(4)和式(9)是假定N已知的情況得到的,如果選擇N太小,將損壞有效信號; 如果選擇N太大,將降低噪聲壓制效果。對于MSSA類方法,自動確定N是關(guān)鍵。面對復(fù)雜多變的實(shí)際地震數(shù)據(jù),在沒有充足的地質(zhì)資料時,確定數(shù)據(jù)塊中需要保留的N是一個值得研究的問題。實(shí)際數(shù)據(jù)去噪時要對數(shù)據(jù)分塊,需要人工估計(jì)每個塊的有效N,計(jì)算效率低,無法實(shí)現(xiàn)產(chǎn)業(yè)化。為此,首先引入一種確定有效N的經(jīng)驗(yàn)估計(jì)方法,然后給出了一種自動確定N的方法,該方法基于Akaike信息準(zhǔn)則自動確定有效N,有利于MSSA產(chǎn)業(yè)化。

    1.2 經(jīng)驗(yàn)公式法

    圖1為模擬地震數(shù)據(jù)A三維視圖。模擬數(shù)據(jù)使用主頻為40Hz的雷克子波作為震源,信噪比(SNR)定義為

    圖1 模擬地震數(shù)據(jù)A三維視圖(a)無噪數(shù)據(jù); (b)加入10%隨機(jī)噪聲數(shù)據(jù)(信噪比為-1.322dB)數(shù)據(jù)含有5個不同傾角的地震同相軸,時間方向有300個采樣點(diǎn),采樣率為2ms,主測線和聯(lián)絡(luò)測線方向均為60個采樣點(diǎn)

    (10)

    式中:d為不含噪數(shù)據(jù);r為去噪后數(shù)據(jù)。

    圖2為圖1數(shù)據(jù)在ω12處的σ曲線。由圖可見,地震數(shù)據(jù)在無噪聲時僅出現(xiàn)少數(shù)較大的σ且個數(shù)為N,其他σ均較小,因此σ曲線出現(xiàn)明顯的彎折現(xiàn)象(圖2b紅線)。從理論上來說,無噪時N以后的σ應(yīng)該全部為零。數(shù)據(jù)加入噪聲后,σ發(fā)生改變的同時也出現(xiàn)大量非零σ,因此加噪后的σ曲線下降相對平緩,但在第N和第N+1項(xiàng)之間,依然存在巨大的落差。

    通過區(qū)分?jǐn)?shù)據(jù)含噪和不含噪情況的σ,從而確定數(shù)據(jù)塊中要保留的N。基于錯位相除的思路,提出一種確定數(shù)據(jù)塊中有效N的經(jīng)驗(yàn)估計(jì)方法。首先將Σ的σ排列成新的向量

    Q=(σ1,σ2,…,σd)

    (11)

    對Q錯位相除,得

    (12)

    定義新的向量

    Q′=(q1,q2,…,qd-1)

    (13)

    式中Q′是基于Q的元素錯位相除構(gòu)造的新向量,其第i個元素等于Q的第i+1個元素除以第i個元素。若σi和σi+1的值較接近,則qi的值趨于1,否則將出現(xiàn)一個較小的值。對Q′取極小,并獲取極小值對應(yīng)的索引,有

    N=arg minQ′=(q1,q2,…,qd-1)

    (14)

    在實(shí)際處理中,只需截取Q′的前若干項(xiàng)即可。圖3為錯位相除向量曲線。由圖可見,Q′的第5個元素值明顯小于第4個元素和第6個元素,說明在該頻率下應(yīng)該保留5個奇異值(圖3b)。

    為了排除偶然性,使用上述方法對數(shù)據(jù)塊的主頻數(shù)據(jù)進(jìn)行相同的計(jì)算,并得到新向量

    flag1=(arg minQ′ωL,…,arg minQ′ωH)

    (15)

    其中ωL~ωH為信號的主要頻率范圍。通過統(tǒng)計(jì)flag1中不同的N出現(xiàn)的占比,根據(jù)經(jīng)驗(yàn)選擇合適的N。圖4為由經(jīng)驗(yàn)公式法確定的圖1a的σ分布。由圖可見:不論是否含有噪聲,在主頻范圍內(nèi)N均為5(圖4a); 含噪聲數(shù)據(jù)和不含噪聲數(shù)據(jù)N=5的占比分別為0.5、0.35,故確定N為5(圖4b)。因此,無論數(shù)據(jù)中是否包含噪聲,利用經(jīng)驗(yàn)公式法均可準(zhǔn)確估計(jì)N。經(jīng)驗(yàn)公式法所用的錯位相除策略和一階差分具有異曲同工的效果,可以驗(yàn)證其他方法的效果,為確定N提供了有效工具。

    圖2 圖1數(shù)據(jù)在ω12處的σ曲線(a)全部σ分布; (b)圖a前30個σ分布局部放大

    圖3 錯位相除向量曲線(a)完整曲線; (b)圖a前10項(xiàng)元素的局部放大

    圖4 由經(jīng)驗(yàn)公式法確定的圖1a的σ分布(a)全頻率掃描結(jié)果; (b)統(tǒng)計(jì)結(jié)果

    1.3 ADMSSA算法

    處理海量地震數(shù)據(jù)時需要分塊地震數(shù)據(jù),然后對每個數(shù)據(jù)塊去噪。但是每個數(shù)據(jù)塊對應(yīng)的N并不相等,人工預(yù)估方法不利于算法實(shí)施,因此需要研究自適應(yīng)算法。

    對σ序列Q可依據(jù)第N和第N+1個σ之間的巨大落差并伴隨嚴(yán)重的彎折現(xiàn)象確定有效N。事實(shí)上,N值的選擇就是檢測σ序列中的拐點(diǎn)位置。本文利用Akaike信息準(zhǔn)則[28-29]自動判定保留的N。首先對Q作如下變換

    fωi(σμ)=σμ+1-2σμ+σμ-11≤μ≤d

    (16)

    式(16)實(shí)際上是求σ序列曲線的二階導(dǎo)數(shù),fωi(σμ)描述σ曲線斜率的變化率。確定頻率為ωi、第R點(diǎn)的有效N值的Akaike信息準(zhǔn)則為

    AICωi(R)=Rlg[var(fωi[σi,σR])]+

    (d-R-1)lg[var(fωi[σR+1,σd])]

    (17)

    式中:var表示數(shù)據(jù)序列的方差; AICωi(R)是長度為d的序列,其全局最小值對應(yīng)的位置即為拐點(diǎn),按

    flag2=(arg minAICωL,…,arg min AICωH)

    (18)

    求出所有頻率中的最小值。式(18)中元素最小值即是整個數(shù)據(jù)塊中需要保留的N。圖5為由ADMSSA算法確定的圖1a的σ分布。由圖可見:①在N=8時AICωi(R)曲線取極小值(圖5a)。由于主頻范圍內(nèi)數(shù)值結(jié)果較穩(wěn)定,其他范圍則經(jīng)常出現(xiàn)異常值,為了提高精度將頻率控制在10~90Hz內(nèi)。②ADMSSA算法確定的N為8(圖5b),經(jīng)驗(yàn)公式法確定的N為5(圖4b),證明利用基于Akaike信息準(zhǔn)則的方法去噪,可以自動地估計(jì)出與真實(shí)值較接近的N。

    在使用上述方法確定N后,去噪過程采用DMSSA方法的框架,在ADMSSA算法中僅需要確定信號的主頻范圍,就可以自動地去噪。

    圖5 由ADMSSA算法確定的圖1a的σ分布(a)信息準(zhǔn)則曲線; (b)信息準(zhǔn)則確定的結(jié)果

    2 數(shù)值實(shí)驗(yàn)

    2.1 模擬數(shù)據(jù)實(shí)驗(yàn)

    分別使用ADMSSA、DMSSA方法對圖1b去噪,結(jié)果(圖6)表明:由ADMSSA、DMSSA方法確定的N分別為8、5,令阻尼因子D=3,去噪結(jié)果的信噪比分別為22.110(圖6a)、22.438dB(圖6c),兩者的去噪效果較接近,信噪比較高,同相軸清晰連貫,局部細(xì)節(jié)得以保留。

    圖6 圖1b的去噪效果對比(a)ADMSSA方法去噪結(jié)果; (b)ADMSSA方法去除的噪聲; (c) DMSSA方法去噪結(jié)果; (d) DMSSA方法去除的噪聲

    圖7 模擬地震數(shù)據(jù)B三維視圖(a)無噪數(shù)據(jù); (b)加入15%隨機(jī)噪聲數(shù)據(jù)(信噪比為-4.659dB)數(shù)據(jù)含有4個不同傾角的地震同相軸,時間方向有300個采樣點(diǎn),主測線和聯(lián)絡(luò)測線方向均為60個采樣點(diǎn),采用主頻為40Hz的雷克子波正演

    能量較強(qiáng)的噪聲經(jīng)常使地震信號發(fā)生嚴(yán)重畸變,導(dǎo)致塊Hankel矩陣σ變化復(fù)雜。為了驗(yàn)證ADMSSA方法的有效性和對噪聲的敏感程度,對含噪地震數(shù)據(jù)(圖7b)去噪,結(jié)果表明,DMSSA方法確定的N為3(圖8紅色實(shí)線),ADMSSA方法確定的N為7(圖8藍(lán)色實(shí)線)。

    圖9為圖7b的去噪效果對比。由圖可見,令阻尼因子D=3,ADMSSA、DMSSA方法去噪結(jié)果的信噪比分別為21.263(圖9a)、21.778dB(圖9c),即前者的信噪比略低,但去噪結(jié)果的同相軸清晰,噪聲殘留較少(圖9b),說明ADMSSA方法高效、精確。

    2.2 實(shí)際地震數(shù)據(jù)實(shí)例

    為了證明ADMSSA算法對實(shí)際地震數(shù)據(jù)的去噪效果,分別使用二維和三維疊后地震數(shù)據(jù)(圖10)驗(yàn)證。二維地震數(shù)據(jù)(圖10a)信噪比低,地震同相軸連續(xù)性差,隨機(jī)噪聲能量強(qiáng),剖面中間部分以及下部存在斷層,尤其是下部存在多處斷裂構(gòu)造。三維地震數(shù)據(jù)(圖10b)信噪比低,噪聲能量較強(qiáng),有效信號被噪聲嚴(yán)重污染,中間部分同相軸出現(xiàn)彎曲、斷裂現(xiàn)象。

    圖8 兩種方法對圖7b數(shù)據(jù)確定的N

    圖9 圖7b的去噪效果對比(a)ADMSSA方法去噪結(jié)果; (b)ADMSSA方法去除的噪聲; (c)DMSSA方法去噪結(jié)果; (d)DMSSA方法去除的噪聲

    圖11為圖10a的σ分布。由圖可見:ADMSSA方法確定的N為8(圖11a藍(lán)色實(shí)線); DMSSA方法統(tǒng)計(jì)N出現(xiàn)的百分比(圖11b)確定的N為10。圖12為二維地震數(shù)據(jù)去噪效果對比。由圖可見,令阻尼因子D=5,地震同相軸邊緣刻畫清晰,噪聲去除徹底,對構(gòu)造細(xì)節(jié)保護(hù)較好,去噪效果均較好。圖13為圖10b的σ分布。由圖可見:ADMSSA方法確定的N為5(圖13a藍(lán)色實(shí)線); DMSSA方法統(tǒng)計(jì)N出現(xiàn)的百分比確定的N為4(圖13b)。

    圖14為三維地震數(shù)據(jù)去噪效果對比。由圖可見,兩種方法去噪結(jié)果基本相同,無論是在地下結(jié)構(gòu)較穩(wěn)定區(qū)域還是在斷點(diǎn)附近,ADMSSA方法的去噪結(jié)果很好地保護(hù)了構(gòu)造細(xì)節(jié)(圖14a),同相軸的輪廓清晰,去噪效果很好。

    圖10 二維(a)、三維(b)地震數(shù)據(jù)

    圖11 圖10a的σ分布(a)DMSSA、ADMSSA方法確定的N; (b)DMSSA方法統(tǒng)計(jì)結(jié)果

    圖12 二維地震數(shù)據(jù)去噪效果對比(a)ADMSSA方法去噪結(jié)果; (b)ADMSSA方法去除的噪聲; (c)DMSSA方法去噪結(jié)果; (d)DMSSA方法去除的噪聲

    圖13 圖10b的σ分布(a) DMSSA、ADMSSA方法確定的N; (b)DMSSA方法統(tǒng)計(jì)結(jié)果

    圖14 三維地震數(shù)據(jù)去噪效果對比(a)ADMSSA方法去噪結(jié)果; (b)ADMSSA方法去除的噪聲; (c)DMSSA方法去噪結(jié)果; (d)DMSSA方法去除的噪聲

    3 討論

    對于MSSA類方法來說,劃分的數(shù)據(jù)塊的尺度對計(jì)算時間和去噪效果均有影響。ADMSSA方法同樣受到劃分?jǐn)?shù)據(jù)塊尺度的影響,因此在劃分?jǐn)?shù)據(jù)塊時應(yīng)先做試驗(yàn),再確定劃分尺度。ADMSSA方法需要選擇一個合適的頻率范圍,一般取信號有效頻率范圍即可取得較好去噪結(jié)果。DMSSA方法受阻尼因子D的影響,當(dāng)噪聲能量較強(qiáng)時,選擇D約為3,當(dāng)噪聲信號較弱時,應(yīng)選擇較大的D值。ADMSSA方法對強(qiáng)脈沖噪聲具有一定壓制作用,但是去噪結(jié)果中仍殘存一些強(qiáng)脈沖噪聲。MSSA方法對奇異值個數(shù)N依賴很強(qiáng),當(dāng)N較小時,會影響斷層的識別。對于曲率較大的彎曲同相軸,需要合理劃分?jǐn)?shù)據(jù)塊尺度。MSSA方法的計(jì)算量主要為奇異值分解,為了提升算法的計(jì)算效率,可以采用隨機(jī)奇異值分解等方法。

    DMSSA方法和ADMSSA方法確定的N存在差別,這是由于前者依靠經(jīng)驗(yàn)公式法估計(jì)N,存在一定誤差,而ADMSSA方法依靠數(shù)據(jù)信息通過概率分析提取信號特征。

    降秩類去噪算法在消除隨機(jī)噪聲的同時會對有效信號造成一定損傷,其原因是該類方法基于線性同相軸假設(shè),實(shí)際數(shù)據(jù)很難滿足條件。因此,盡管降秩類去噪算法獲得了較好效果,但是去噪效果尚有很大的提升空間。

    4 結(jié)論

    對于MSSA類去噪方法來說,確定有效的奇異值個數(shù)是關(guān)鍵。目前都是依靠人工經(jīng)驗(yàn)估計(jì)奇異值個數(shù),不利于該類方法的產(chǎn)業(yè)化。本文提出了一種確定有效奇異值個數(shù)的方法,該方法基于Akaike信息準(zhǔn)則自動區(qū)分有效信號對應(yīng)的奇異值與噪聲相關(guān)的奇異值,克服了人工選擇奇異值個數(shù)的問題,有利于海量地震數(shù)據(jù)去噪。此外,本文還提出利用經(jīng)驗(yàn)公式法驗(yàn)證ADMSSA方法的可靠性。數(shù)值實(shí)驗(yàn)證明,ADMSSA方法能夠自動地確定可靠的奇異值個數(shù),并且獲得高信噪比的去噪結(jié)果,該算法在工業(yè)化應(yīng)用中具有巨大潛力。

    猜你喜歡
    同相軸個數(shù)信噪比
    怎樣數(shù)出小正方體的個數(shù)
    虛同相軸方法及其在陸上地震層間多次波壓制中的應(yīng)用
    等腰三角形個數(shù)探索
    基于深度學(xué)習(xí)的無人機(jī)數(shù)據(jù)鏈信噪比估計(jì)算法
    怎樣數(shù)出小木塊的個數(shù)
    怎樣數(shù)出小正方體的個數(shù)
    低信噪比下LFMCW信號調(diào)頻參數(shù)估計(jì)
    電子測試(2018年11期)2018-06-26 05:56:02
    低信噪比下基于Hough變換的前視陣列SAR稀疏三維成像
    一種改進(jìn)的相關(guān)法自動拾取同相軸
    一種反射同相軸自動拾取算法
    日韩成人在线观看一区二区三区| a在线观看视频网站| 亚洲五月婷婷丁香| 国产真实乱freesex| 免费在线观看成人毛片| 欧美成人免费av一区二区三区| 在线看三级毛片| 日本在线视频免费播放| 久久热在线av| 不卡av一区二区三区| 日韩中文字幕欧美一区二区| 我的老师免费观看完整版| 天天躁狠狠躁夜夜躁狠狠躁| 精华霜和精华液先用哪个| 久久精品国产清高在天天线| 午夜免费观看网址| 亚洲国产欧美人成| 亚洲国产精品成人综合色| 成人av在线播放网站| 狂野欧美激情性xxxx| 制服人妻中文乱码| 波多野结衣巨乳人妻| 国产亚洲精品综合一区在线观看 | 国产主播在线观看一区二区| 亚洲国产精品999在线| 国产三级在线视频| 国产高清激情床上av| 国产午夜精品久久久久久| 成在线人永久免费视频| 国产精品久久久久久人妻精品电影| av超薄肉色丝袜交足视频| 母亲3免费完整高清在线观看| 午夜福利在线观看吧| 亚洲成av人片免费观看| 日本一区二区免费在线视频| 精品欧美一区二区三区在线| www.熟女人妻精品国产| 亚洲一区中文字幕在线| 男女做爰动态图高潮gif福利片| 久久久精品欧美日韩精品| 久久久国产成人免费| 亚洲在线自拍视频| 欧美最黄视频在线播放免费| 久久精品亚洲精品国产色婷小说| 亚洲精品一区av在线观看| 琪琪午夜伦伦电影理论片6080| tocl精华| 一级毛片精品| 国产三级中文精品| 女生性感内裤真人,穿戴方法视频| 少妇裸体淫交视频免费看高清 | 成人国产综合亚洲| 久久久久国产精品人妻aⅴ院| 国产成+人综合+亚洲专区| 国语自产精品视频在线第100页| 午夜两性在线视频| 嫩草影视91久久| 亚洲美女黄片视频| 一级a爱片免费观看的视频| 一区福利在线观看| www日本黄色视频网| 丰满人妻一区二区三区视频av | 熟女少妇亚洲综合色aaa.| 脱女人内裤的视频| 12—13女人毛片做爰片一| 99精品欧美一区二区三区四区| 18禁裸乳无遮挡免费网站照片| www.熟女人妻精品国产| 老汉色∧v一级毛片| 88av欧美| 欧美成狂野欧美在线观看| 观看免费一级毛片| 999久久久国产精品视频| 天堂√8在线中文| 国产av又大| 50天的宝宝边吃奶边哭怎么回事| 此物有八面人人有两片| 日本a在线网址| 欧美3d第一页| 亚洲专区国产一区二区| 免费av毛片视频| 久久精品成人免费网站| 国产99久久九九免费精品| 国产午夜福利久久久久久| 亚洲精品色激情综合| 亚洲五月婷婷丁香| 真人做人爱边吃奶动态| 女人高潮潮喷娇喘18禁视频| 国产成人精品无人区| 91麻豆av在线| 午夜两性在线视频| 国产精品一区二区三区四区久久| 一进一出抽搐动态| 欧美精品啪啪一区二区三区| 精品国产美女av久久久久小说| 色精品久久人妻99蜜桃| 亚洲在线自拍视频| 亚洲av成人不卡在线观看播放网| 午夜福利免费观看在线| 亚洲av电影不卡..在线观看| 黄色成人免费大全| 天堂√8在线中文| 女警被强在线播放| 天天一区二区日本电影三级| 国产私拍福利视频在线观看| 一进一出抽搐动态| 日韩精品中文字幕看吧| 亚洲男人天堂网一区| 日韩高清综合在线| 国产99久久九九免费精品| 91字幕亚洲| 欧美日韩中文字幕国产精品一区二区三区| www国产在线视频色| 婷婷亚洲欧美| 久久99热这里只有精品18| 丝袜人妻中文字幕| 中文字幕av在线有码专区| 伊人久久大香线蕉亚洲五| 精品久久久久久久毛片微露脸| 欧美色欧美亚洲另类二区| 一本大道久久a久久精品| 精华霜和精华液先用哪个| 母亲3免费完整高清在线观看| 级片在线观看| 神马国产精品三级电影在线观看 | 国产精品1区2区在线观看.| 国产av一区二区精品久久| 一本久久中文字幕| 国产亚洲精品综合一区在线观看 | 一本精品99久久精品77| 大型黄色视频在线免费观看| 亚洲av熟女| svipshipincom国产片| 日韩中文字幕欧美一区二区| 最近视频中文字幕2019在线8| 久久久久久亚洲精品国产蜜桃av| 小说图片视频综合网站| videosex国产| √禁漫天堂资源中文www| 欧美三级亚洲精品| 91av网站免费观看| 91老司机精品| 久久久国产欧美日韩av| 久久久久九九精品影院| 中文资源天堂在线| 宅男免费午夜| 欧美丝袜亚洲另类 | 国产视频内射| 国产精品久久电影中文字幕| 波多野结衣巨乳人妻| 亚洲 欧美 日韩 在线 免费| 亚洲成人中文字幕在线播放| 一个人免费在线观看的高清视频| 最近在线观看免费完整版| 美女黄网站色视频| 男女午夜视频在线观看| 男人舔女人下体高潮全视频| 91av网站免费观看| 后天国语完整版免费观看| 久久精品国产99精品国产亚洲性色| 欧美一级a爱片免费观看看 | 人人妻,人人澡人人爽秒播| 午夜免费成人在线视频| 18禁国产床啪视频网站| 国产精品自产拍在线观看55亚洲| 两性夫妻黄色片| 国产精品亚洲av一区麻豆| 久久久久久久久免费视频了| 色老头精品视频在线观看| 无遮挡黄片免费观看| 搡老妇女老女人老熟妇| 淫妇啪啪啪对白视频| 男人舔奶头视频| 51午夜福利影视在线观看| 欧美国产日韩亚洲一区| 男人舔女人的私密视频| 老鸭窝网址在线观看| 91老司机精品| 日韩 欧美 亚洲 中文字幕| 中国美女看黄片| 人妻丰满熟妇av一区二区三区| 99国产极品粉嫩在线观看| 日本a在线网址| 国产人伦9x9x在线观看| 免费看a级黄色片| 久久天躁狠狠躁夜夜2o2o| 国产激情欧美一区二区| 亚洲成人中文字幕在线播放| 一本一本综合久久| 在线十欧美十亚洲十日本专区| 久久精品aⅴ一区二区三区四区| 成年人黄色毛片网站| 亚洲国产日韩欧美精品在线观看 | 亚洲专区中文字幕在线| 亚洲国产欧美一区二区综合| 91在线观看av| 亚洲精品美女久久久久99蜜臀| 热99re8久久精品国产| 老司机深夜福利视频在线观看| 久久这里只有精品中国| 可以在线观看的亚洲视频| www日本在线高清视频| 国产私拍福利视频在线观看| 久久精品国产99精品国产亚洲性色| 久久香蕉精品热| av有码第一页| 中文字幕最新亚洲高清| 久久欧美精品欧美久久欧美| 一个人观看的视频www高清免费观看 | 亚洲一码二码三码区别大吗| 精品久久久久久久人妻蜜臀av| 一a级毛片在线观看| 精品欧美国产一区二区三| 女人爽到高潮嗷嗷叫在线视频| 国产精品久久久人人做人人爽| 亚洲人成电影免费在线| 黑人巨大精品欧美一区二区mp4| 久久精品综合一区二区三区| 国产av一区在线观看免费| 18禁观看日本| av超薄肉色丝袜交足视频| 免费一级毛片在线播放高清视频| 国产精品,欧美在线| aaaaa片日本免费| 国产精品99久久99久久久不卡| 亚洲国产欧美人成| xxx96com| 婷婷六月久久综合丁香| 国产97色在线日韩免费| 12—13女人毛片做爰片一| 国产成人精品无人区| 日日爽夜夜爽网站| 99热只有精品国产| www.精华液| 91字幕亚洲| 久久精品91无色码中文字幕| 99国产精品99久久久久| 国产精品久久久久久精品电影| 99热6这里只有精品| 国产精品98久久久久久宅男小说| 国产亚洲精品久久久久5区| 91老司机精品| 亚洲精品久久成人aⅴ小说| 18禁黄网站禁片午夜丰满| 婷婷精品国产亚洲av| 欧美一区二区国产精品久久精品 | 欧美av亚洲av综合av国产av| 亚洲人成电影免费在线| 动漫黄色视频在线观看| 久久精品91无色码中文字幕| 成人永久免费在线观看视频| av有码第一页| 欧美大码av| 亚洲真实伦在线观看| 丝袜人妻中文字幕| 亚洲av五月六月丁香网| 国产黄色小视频在线观看| 人妻丰满熟妇av一区二区三区| 国产又色又爽无遮挡免费看| 国产69精品久久久久777片 | 久久久水蜜桃国产精品网| 蜜桃久久精品国产亚洲av| 亚洲 国产 在线| 国产成人影院久久av| 亚洲专区字幕在线| 久久久久久国产a免费观看| 久久伊人香网站| 日本五十路高清| 91国产中文字幕| 日本熟妇午夜| 19禁男女啪啪无遮挡网站| 国产精品永久免费网站| 日韩精品免费视频一区二区三区| 国产精品亚洲av一区麻豆| 日韩欧美 国产精品| 男女视频在线观看网站免费 | 亚洲黑人精品在线| 一本大道久久a久久精品| 变态另类丝袜制服| 久久久精品大字幕| 亚洲,欧美精品.| 欧美黄色淫秽网站| 亚洲自偷自拍图片 自拍| 超碰成人久久| 久久 成人 亚洲| 国产99白浆流出| 美女高潮喷水抽搐中文字幕| 制服诱惑二区| 可以在线观看毛片的网站| 免费看a级黄色片| 日日干狠狠操夜夜爽| www.精华液| 精品一区二区三区视频在线观看免费| 19禁男女啪啪无遮挡网站| av片东京热男人的天堂| 国产免费男女视频| 一区福利在线观看| 青草久久国产| 免费一级毛片在线播放高清视频| 午夜福利18| av超薄肉色丝袜交足视频| 亚洲av成人精品一区久久| 伊人久久大香线蕉亚洲五| 亚洲欧美一区二区三区黑人| 亚洲男人的天堂狠狠| 麻豆久久精品国产亚洲av| 久久亚洲真实| 99国产精品一区二区三区| x7x7x7水蜜桃| 高清毛片免费观看视频网站| 中文字幕av在线有码专区| 欧美性长视频在线观看| 免费在线观看完整版高清| 国产亚洲欧美98| 手机成人av网站| avwww免费| 久久久国产成人免费| 不卡av一区二区三区| 九色成人免费人妻av| 成人亚洲精品av一区二区| 夜夜看夜夜爽夜夜摸| 精品久久久久久久毛片微露脸| 久久久国产成人精品二区| 老熟妇乱子伦视频在线观看| 精品乱码久久久久久99久播| 后天国语完整版免费观看| 波多野结衣巨乳人妻| 精华霜和精华液先用哪个| 91国产中文字幕| 国产成人精品久久二区二区免费| 欧美av亚洲av综合av国产av| 亚洲精品久久国产高清桃花| 久久久久久久久中文| 午夜福利欧美成人| 这个男人来自地球电影免费观看| www.精华液| 国产v大片淫在线免费观看| av视频在线观看入口| 成人一区二区视频在线观看| 搡老岳熟女国产| 黄片大片在线免费观看| 国产成人精品久久二区二区免费| 巨乳人妻的诱惑在线观看| 亚洲色图 男人天堂 中文字幕| 欧美性长视频在线观看| 亚洲,欧美精品.| 久99久视频精品免费| 国产激情欧美一区二区| 久久国产乱子伦精品免费另类| 丝袜美腿诱惑在线| 久久午夜综合久久蜜桃| 亚洲国产精品sss在线观看| 成人手机av| 亚洲欧美日韩东京热| 五月玫瑰六月丁香| 色综合欧美亚洲国产小说| 国产一区二区三区视频了| 69av精品久久久久久| 一级毛片高清免费大全| 国产高清激情床上av| aaaaa片日本免费| 国产精品一区二区精品视频观看| 久久久精品欧美日韩精品| 精品久久久久久久久久免费视频| 日本免费一区二区三区高清不卡| 变态另类丝袜制服| 国产单亲对白刺激| 超碰成人久久| 亚洲男人天堂网一区| x7x7x7水蜜桃| 亚洲成人免费电影在线观看| 性色av乱码一区二区三区2| 午夜久久久久精精品| 国产精品影院久久| 成年人黄色毛片网站| 国产精品1区2区在线观看.| 一夜夜www| 中文字幕高清在线视频| 久久中文看片网| 欧美极品一区二区三区四区| 此物有八面人人有两片| 嫁个100分男人电影在线观看| 欧美日韩瑟瑟在线播放| 国产午夜精品久久久久久| 久久久久国产一级毛片高清牌| 精品国产超薄肉色丝袜足j| 亚洲人成77777在线视频| 欧美中文综合在线视频| 大型黄色视频在线免费观看| 老鸭窝网址在线观看| 久久 成人 亚洲| 搡老岳熟女国产| 男女视频在线观看网站免费 | 可以在线观看毛片的网站| 在线免费观看的www视频| 桃色一区二区三区在线观看| 精品少妇一区二区三区视频日本电影| 级片在线观看| 男人舔奶头视频| 桃红色精品国产亚洲av| 国产熟女xx| 国产一区二区在线av高清观看| 欧美日韩亚洲国产一区二区在线观看| 身体一侧抽搐| 嫁个100分男人电影在线观看| 久热爱精品视频在线9| 五月玫瑰六月丁香| 中出人妻视频一区二区| 99国产极品粉嫩在线观看| 青草久久国产| 免费在线观看视频国产中文字幕亚洲| 国产私拍福利视频在线观看| 国产日本99.免费观看| 人人妻,人人澡人人爽秒播| 在线观看免费午夜福利视频| www日本在线高清视频| 国产av一区二区精品久久| 久久久久久亚洲精品国产蜜桃av| 亚洲精品在线美女| 免费人成视频x8x8入口观看| 看片在线看免费视频| 国产在线观看jvid| 非洲黑人性xxxx精品又粗又长| 久久香蕉精品热| 亚洲国产精品成人综合色| 一级毛片高清免费大全| 丰满人妻一区二区三区视频av | 亚洲中文字幕日韩| 99国产精品一区二区蜜桃av| 成人精品一区二区免费| АⅤ资源中文在线天堂| 国产亚洲欧美98| 亚洲 国产 在线| 欧美人与性动交α欧美精品济南到| 欧美激情久久久久久爽电影| 少妇裸体淫交视频免费看高清 | 91字幕亚洲| 黄片大片在线免费观看| 亚洲国产精品成人综合色| 亚洲国产中文字幕在线视频| 男人舔奶头视频| 1024手机看黄色片| 午夜激情av网站| 一区福利在线观看| 三级毛片av免费| 国产精品亚洲一级av第二区| 51午夜福利影视在线观看| 日韩欧美在线乱码| 少妇被粗大的猛进出69影院| 99久久精品国产亚洲精品| 两个人的视频大全免费| 国产精品久久久久久人妻精品电影| 国产免费男女视频| 国产精品自产拍在线观看55亚洲| 亚洲av片天天在线观看| 精品久久久久久久毛片微露脸| 欧美日本视频| 久久中文字幕一级| 亚洲av中文字字幕乱码综合| 99精品欧美一区二区三区四区| 观看免费一级毛片| 国产精品,欧美在线| 高清在线国产一区| 男人舔女人下体高潮全视频| 18禁黄网站禁片午夜丰满| 久久久久久人人人人人| 中亚洲国语对白在线视频| 欧美不卡视频在线免费观看 | 久久午夜综合久久蜜桃| 日本五十路高清| 后天国语完整版免费观看| 三级男女做爰猛烈吃奶摸视频| 搡老妇女老女人老熟妇| 国产亚洲精品综合一区在线观看 | 国产精品久久电影中文字幕| 女人高潮潮喷娇喘18禁视频| 久久久精品国产亚洲av高清涩受| 日韩免费av在线播放| 最近在线观看免费完整版| 色哟哟哟哟哟哟| 成在线人永久免费视频| 一本精品99久久精品77| 两个人免费观看高清视频| 中出人妻视频一区二区| www日本黄色视频网| 可以免费在线观看a视频的电影网站| 久久久久久人人人人人| 禁无遮挡网站| 亚洲熟女毛片儿| 亚洲精华国产精华精| 777久久人妻少妇嫩草av网站| 我要搜黄色片| 精品一区二区三区av网在线观看| 久久这里只有精品中国| 亚洲国产中文字幕在线视频| 丁香欧美五月| 免费看日本二区| 国产亚洲精品第一综合不卡| 三级毛片av免费| 99久久久亚洲精品蜜臀av| tocl精华| 精品国产超薄肉色丝袜足j| a级毛片在线看网站| 亚洲av五月六月丁香网| 欧美日韩亚洲综合一区二区三区_| 日日摸夜夜添夜夜添小说| 黄色女人牲交| 久久伊人香网站| 久久精品国产亚洲av香蕉五月| www.精华液| 两个人免费观看高清视频| 国产久久久一区二区三区| 香蕉av资源在线| 亚洲欧美日韩高清在线视频| 99热只有精品国产| 国产在线精品亚洲第一网站| 91麻豆精品激情在线观看国产| 久久午夜综合久久蜜桃| 亚洲,欧美精品.| 国产人伦9x9x在线观看| 草草在线视频免费看| 丰满人妻熟妇乱又伦精品不卡| 日韩欧美一区二区三区在线观看| 欧美中文综合在线视频| 午夜a级毛片| 成人国语在线视频| 久久久国产欧美日韩av| 又爽又黄无遮挡网站| avwww免费| 日本三级黄在线观看| 亚洲,欧美精品.| 两个人看的免费小视频| 日韩三级视频一区二区三区| 窝窝影院91人妻| 日韩三级视频一区二区三区| 丁香六月欧美| 欧美日韩亚洲综合一区二区三区_| 亚洲最大成人中文| 久久天躁狠狠躁夜夜2o2o| 亚洲人成网站在线播放欧美日韩| www.www免费av| 老汉色∧v一级毛片| 亚洲人成电影免费在线| 男人舔奶头视频| 欧美成人免费av一区二区三区| 日韩欧美一区二区三区在线观看| 黄色a级毛片大全视频| a级毛片a级免费在线| 欧美绝顶高潮抽搐喷水| 听说在线观看完整版免费高清| 在线国产一区二区在线| 亚洲熟妇熟女久久| 国产精品电影一区二区三区| 中文字幕熟女人妻在线| 国产精品乱码一区二三区的特点| 1024视频免费在线观看| 动漫黄色视频在线观看| 后天国语完整版免费观看| 99久久久亚洲精品蜜臀av| 国产成人一区二区三区免费视频网站| 久久久久久人人人人人| tocl精华| 特级一级黄色大片| 日韩三级视频一区二区三区| 校园春色视频在线观看| 美女免费视频网站| 夜夜躁狠狠躁天天躁| 亚洲熟妇熟女久久| netflix在线观看网站| 听说在线观看完整版免费高清| 国产精品一区二区三区四区免费观看 | 在线观看一区二区三区| 国产久久久一区二区三区| 88av欧美| 夜夜躁狠狠躁天天躁| 亚洲人成77777在线视频| 日韩欧美在线乱码| 在线观看美女被高潮喷水网站 | 成人国产综合亚洲| 欧美成人一区二区免费高清观看 | 黄色成人免费大全| 亚洲成av人片免费观看| 成人国产一区最新在线观看| 成人国语在线视频| 两个人免费观看高清视频| 午夜精品一区二区三区免费看| 天天添夜夜摸| 国产精品免费视频内射| 一本一本综合久久| 久久人人精品亚洲av| 精品久久久久久成人av| 国产精品一区二区三区四区久久| 香蕉久久夜色| 男女视频在线观看网站免费 | 成人手机av| 久久久久亚洲av毛片大全| 久久草成人影院| 黄片小视频在线播放| 国产成人av激情在线播放| 黄片小视频在线播放| 精品第一国产精品| a级毛片a级免费在线| 免费无遮挡裸体视频| 欧美色视频一区免费| 日韩欧美精品v在线| 精品久久久久久久久久久久久| 国产三级黄色录像| www日本在线高清视频| 亚洲成av人片免费观看| 免费看a级黄色片| 精品久久久久久久久久久久久| 琪琪午夜伦伦电影理论片6080| 久久久久国产精品人妻aⅴ院| 久久九九热精品免费| 成人一区二区视频在线观看|