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

    超聲成像檢測研究進展

    2022-12-02 06:55:22張碧星閻守國
    關鍵詞:信號檢測方法

    張碧星,張 萍,閻守國,黃 娟

    (1 中國科學院 聲學研究所 聲場聲信息國家重點實驗室,北京 100190;2 中國科學院大學, 北京 100190)

    超聲成像檢測是一種利用超聲波在介質(zhì)內(nèi)的傳播特性并以反射或透射聲波作為信息載體的可視化圖像檢測方法。近年來,超聲成像檢測技術取得了長足的發(fā)展,與其他成像方法相比,例如:電子計算機斷層掃描(computed tomography,CT)、磁共振成像(magnetic resonance imaging,MRI)、紅外熱成像(infrared thermal imaging,ITI)等,超聲成像檢測具有安全性高、速度快、成本低等優(yōu)勢,因此被廣泛應用于工業(yè)無損檢測(nondestructive testing,NDT)[1-2]和醫(yī)學診斷[3]等領域。

    在工業(yè)無損檢測中,對于各種設備材料,其內(nèi)部經(jīng)常存在多種缺陷,有的缺陷在使用前就已經(jīng)存在,有的缺陷是在材料服役過程中由于周圍環(huán)境的影響和作用(如腐蝕和疲勞)而產(chǎn)生。這些缺陷對設備材料的安全使用至關重要,不同類型的缺陷對材料的危害程度也有所不同[4],故對設備材料及結(jié)構進行無損檢測和安全評估是一項非常重要的工作[5]。然而,由于缺陷往往存在于材料內(nèi)部,無法通過視覺方式直接觀察,必須利用特殊方法進行檢測[6],超聲成像檢測是一種重要的無損檢測方法。隨著科學技術的不斷發(fā)展,人們對超聲成像檢測的要求越來越高,希望觀察到更加細微的缺陷及其結(jié)構。本文針對工業(yè)無損檢測領域,對幾種重要的超聲成像檢測方法及其發(fā)展現(xiàn)狀進行回顧和總結(jié),在此基礎上,對超聲多波成像檢測方法的研究和發(fā)展趨勢進行分析和展望。

    1 基于自適應和動態(tài)聚焦的成像檢測

    本節(jié)主要分析聲波時間反轉(zhuǎn)、超聲相控陣、合成孔徑和全聚焦等成像檢測方法的發(fā)展現(xiàn)狀。聲波時間反轉(zhuǎn)法是一種自適應聚焦的成像方法,超聲相控陣、合成孔徑和全聚焦方法則屬于動態(tài)聚焦的成像檢測方法。

    1.1 時間反轉(zhuǎn)成像

    聲波時間反轉(zhuǎn)(time reverse, TR)是一種自適應聚焦技術,它不需要介質(zhì)和換能器陣列性質(zhì)和結(jié)構的先驗知識就可以實現(xiàn)聲波自適應聚焦。該方法是法國科學家Fink最早將光學中連續(xù)波的相位共軛法引入到聲學領域并推廣到脈沖波的時間反轉(zhuǎn)法,實現(xiàn)了聲波的自適應聚焦與成像檢測,在超聲無損檢測等領域得到了快速發(fā)展和應用。

    在沒有耗散的聲傳播介質(zhì)中,波動方程具有時間反演不變性,如果聲波位移場u(r,t)是波動方程的解,那么u(r,-t)也是波動方程的解,u(r, -t)就是時間反轉(zhuǎn)聲場。在實際情況中,要構建完整的時間反轉(zhuǎn)聲場是很困難的,也是不可能的,大多都是采用時間反轉(zhuǎn)鏡(time reversal mirror, TRM)進行近似處理。假如一個聲源S向外發(fā)射聲波并在介質(zhì)中傳播,在不同位置用換能器或陣列進行接收,不同陣元接收到的波形具有不同的到時和波形結(jié)構。將換能器各陣元記錄到的信號按時間先后進行反序處理,將反序后的信號重新加載到對應的各陣元上進行激勵,這時各陣元發(fā)出的聲波經(jīng)過介質(zhì)后將自動地在原聲源位置處實現(xiàn)聚焦。這時,每一個換能器陣元就是一個時間反轉(zhuǎn)鏡,類似于光學中的平面鏡,使聲波按原路返回。但由于換能器陣元的個數(shù)有限,不能在空間每一點上進行采樣,且每一陣元記錄聲波的時間點數(shù)也有限,因而時間反轉(zhuǎn)聲場在原聲源處存在空間和時間上的旁瓣。陣元數(shù)越多,各陣元記錄的時間點越多,則時間反轉(zhuǎn)聲場的空間和時間旁瓣就越小。

    Fink研究組在聲波時間反轉(zhuǎn)研究上做出了杰出的工作,取得了大量研究成果,不僅在原理和方法上對時間反轉(zhuǎn)自適應聚焦進行了研究和論證[7-11],而且還針對固體介質(zhì)[12]、分層介質(zhì)[13]、波導介質(zhì)[14]、非均勻介質(zhì)[15]等開展了大量的時間反轉(zhuǎn)理論和實驗研究,并將時間反轉(zhuǎn)法應用到航空航天用鈦合金材料檢測[16]和醫(yī)學上腎結(jié)石及膽結(jié)石粉碎[17]研究中。這些研究表明,時間反轉(zhuǎn)法在各種復雜介質(zhì)和結(jié)構中都能實現(xiàn)自適應聚焦,在不規(guī)則界面和強噪聲介質(zhì)中具有獨特的技術優(yōu)勢。之后,F(xiàn)ink研究組還發(fā)展了TR循環(huán)迭代方法[18],可選擇具體目標信號進行TR循環(huán)迭代處理,使目標信號不斷聚焦和加強,而其他信號相對減弱,該方法在弱信號檢測中特別有效。并在此基礎上提出了DORT方法[19-20],通過TR算子本征值和本征矢量,得到了與每個本征值對應的目標信號,將DORT方法和TR循環(huán)迭代相結(jié)合,就可對所有目標信號進行識別和檢測,實現(xiàn)了多目標的檢測與區(qū)分。

    TR是自適應的聚焦方法,能使來自不同路徑的波[21-22]以及不同類型的波[14, 23]實現(xiàn)同時聚焦。聲波時間反轉(zhuǎn)法是聲學互易原理的結(jié)果,其自適應聚焦特性具有重要的科學意義。時間反轉(zhuǎn)法雖然能自適應地實現(xiàn)聲波聚焦,但時間反轉(zhuǎn)法本身不能對缺陷實現(xiàn)定位,必須借助其他方法或者介質(zhì)的先驗知識等才能對缺陷實現(xiàn)定位。早期的時間反轉(zhuǎn)法不能用于定位,后來采用虛擬TR過程并結(jié)合介質(zhì)先驗知識形成了虛擬時間反轉(zhuǎn)法[24],實現(xiàn)了目標的檢測和定位。

    近年來,人們對TR法進行了不斷改進,如TR-MUSIC法[25],該方法是基于多重信號分類(multiple signal classification, MUSIC)的TR算法,可得到散射體的超分辨率圖像,還可應用于對點散射體密度和相關參數(shù)的有效估計[26]。研究表明,TR-MUSIC方法能分辨出距離瑞利極限更近的橫向目標,實現(xiàn)超分辨成像,還可在強噪聲情況下,抑制偽影并得到穩(wěn)定的圖像[27]。另一個典型的改進是ACU-TR法[28],該方法是結(jié)合了空氣耦合超聲(air-coupled ultrasound)的TR算法,它可以補償衍射效應,從而進行復合材料板材的無損檢測,極大提高了空氣耦合超聲無損檢測圖像的橫向分辨率。另外,將TR與逆時偏移相結(jié)合可區(qū)分界面和界面附近的缺陷目標[29]。雖然近年關于TR研究的文獻較多,但大多都是將TR應用于具體介質(zhì)和目標對象的自適應聚焦與成像研究。

    TR聚焦方法不適用于非線性和耗散介質(zhì),但經(jīng)過適當改進后也能在非線性和耗散介質(zhì)中較好地實現(xiàn)聚焦和成像檢測[30-31]。例如:Dauson等[32]采用時間反轉(zhuǎn)非線性彈性波譜法可以將聲能聚焦到材料中所需位置,可以檢測到實驗室井筒的早期損壞跡象。

    1.2 超聲相控陣成像

    超聲相控陣技術是源于電磁波相控陣雷達技術的原理,超聲相控陣初期主要在醫(yī)學診斷和醫(yī)學治療[33]上使用,現(xiàn)在被廣泛應用于工業(yè)無損檢測中。

    超聲相控陣由多個獨立的陣元組成換能器陣列,通過電子系統(tǒng)控制每個陣元發(fā)射信號的延時和幅度形成相控陣聚焦聲束,通過不斷控制和改變聚焦聲束的形狀和方向,在不需要移動或少移動相控陣換能器探頭位置情況下,實現(xiàn)較大范圍的掃描與成像檢測。隨著聲束偏轉(zhuǎn)方向和焦點位置的不斷變化,聲束對介質(zhì)區(qū)域?qū)崿F(xiàn)全方位的掃查,這是聚焦發(fā)射過程。隨后對換能器各陣元接收到的信號根據(jù)焦點位置不同施加相應的時間延遲進行虛擬聚焦,以實現(xiàn)對相應聲束掃描區(qū)域的成像,這是聚焦接收過程。相對于傳統(tǒng)的單探頭換能器成像方法,超聲相控陣利用各陣元的靈活性對信號的聚焦發(fā)射和聚焦接收兩個過程均進行動態(tài)聚焦,從而可提高回波信號的信噪比和成像分辨率,并達到實時成像的效果。

    隨著計算機技術的快速發(fā)展和換能器陣列制作工藝的不斷提高,超聲相控陣檢測作為一種新型超聲成像檢測技術,最近二三十年來在工業(yè)無損檢測領域得到了飛速發(fā)展。1992 年,美國通用電氣公司成功研制了數(shù)字式超聲相控陣實時成像系統(tǒng);之后,Hatfield等[34]提出了高集成度超聲相控陣系統(tǒng),實現(xiàn)了手持式操作,加拿大R/D Tech公司較早地推出了便攜式相控陣成像設備,以及相關的相控陣產(chǎn)品[35],在超聲相控陣產(chǎn)業(yè)化應用方面占據(jù)領先地位。隨后,國內(nèi)外許多研究機構和公司在超聲相控陣研究上也取得了很大進步,相繼推出了越來越成熟的商業(yè)儀器和設備[36-37]。

    對于外形復雜、具有不規(guī)則界面的被檢對象,傳統(tǒng)的超聲檢測非常困難,常常在遇到復雜界面時,需要改變探頭的位置和方向,檢測效果差。而超聲相控陣檢測技術具有聲束靈活可控、覆蓋面積大、檢測精度高等優(yōu)勢,因而在各種復雜情況下都取得了成功的應用,例如:在汽輪機葉片和渦輪圓盤檢測[38-40]、管道焊縫檢測[41]、火車輪軸檢測[42]、核電站[43]和航空材料檢測[44]、大型鍛造鋼件檢測[45]、復合材料檢測[46]、井壁成像檢測[47-48]等應用領域,超聲相控陣技術發(fā)揮了極為重要的作用。

    相對超聲相控陣的成功應用,超聲相控陣換能器輻射聲場的研究與分析相對較晚。陳啟敏等[49]對相控陣超聲換能器聲場進行了實驗分析;Wooh等[50]較為系統(tǒng)地研究了一維線性相控陣聲場的聚焦特性,詳細分析了相控陣換能器參數(shù)對旁瓣和柵瓣的影響。對于環(huán)形相控陣聲場,Dupenloup等[51]和Zhang等[52]研究了超聲聚焦特性及其與陣列參數(shù)的關系。對于柱面凹形相控陣聲場,F(xiàn)leury等[53]和Zhang等[54]分析了超聲聚焦過程并對套管井壁進行了成像研究,董晗[55]將柱面凹形相控陣探頭應用于井壁成像檢測系統(tǒng),克服了常規(guī)技術中探頭旋轉(zhuǎn)帶來的困難。對于柱面凹形相控陣,王文龍等[56]、張碧星等[57]和Wang等[58]深入研究了聚焦聲場特性及柵瓣控制方法。對于二維相控陣列,Mckee等[59]對水浸狀態(tài)下具有雙曲面輪廓的試件進行三維成像研究,龍絨蓉等[60]發(fā)現(xiàn)二維圓形陣列比二維矩形陣列具有更窄的主瓣寬度和更低的第一級旁瓣。對于柔性相控陣,Jocelyn等[61]和Casula等[62]研究了在不規(guī)則表面檢測時的聲聚焦過程,這些研究為復雜條件下超聲相控陣列探頭的設計提供了可靠的技術參數(shù)。

    此外,在超聲相控陣聲場仿真分析中,高斯聲束法[63-64]可以快速有效地模擬相控陣輻射聲場,相比于波動方程的精確求解能降低幾個數(shù)量級的計算量。Newberry等[65]將高斯函數(shù)展開法應用于圓形活塞聲源輻射聲場的計算;Spies等[66-67]利用多元高斯聲束模型計算了非均勻多層介質(zhì)內(nèi)的超聲傳播;Huang等[68-69]將多元高斯聲束模型擴展到各向異性奧氏體不銹鋼介質(zhì)中;趙新玉等[70]利用多高斯聲束模型進一步模擬計算了超聲相控陣聲場的聚焦特性;Ye等[71]提出了一種基于相控陣的線性相位多高斯波束模型,成功預測了異種金屬焊縫中的相控陣聚焦聲場。高斯聲束法為相控陣聲場仿真提供了較好的計算基礎。

    常規(guī)超聲相控陣成像檢測技術基于信號延時和幅度來控制聲束的聚焦和偏轉(zhuǎn),可以實現(xiàn)不同角度的實時掃描成像。但是為了確保成像實時性,超聲相控陣列發(fā)射聲束的數(shù)量和聚焦點數(shù)相對受限,導致檢測精度和分辨率相對一般,這使得成像精度更高的超聲相控陣后處理成像研究得到了越來越多的重視。

    1.3 合成孔徑成像

    20世紀60年代,F(xiàn)laherty等[72]在超聲成像領域提出了合成孔徑的概念。此后,合成孔徑聚焦技術(synthetic aperture focusing technique, SAFT)廣泛地應用于超聲檢測領域。相比于傳統(tǒng)聚焦技術,SAFT作為一種超聲后處理方法,能夠?qū)⑿】讖綋Q能器陣列合成為一個大孔徑陣列,以提高縱向檢測深度和圖像分辨率,在使用同樣換能器陣列探頭的情況下,可得到更高分辨率的重建圖像,為缺陷的成像和定性分析提供更好的技術基礎。

    SAFT成像方法采用孔徑較小的換能器陣列(單陣元或陣元數(shù)較少的陣列)在某位置發(fā)射并接收來自反射體的回波信號,之后該陣列移動到另一位置再發(fā)射并接收回波信號,以后以相同的掃描步距依次在不同位置處發(fā)射并接收反射回波信號[73];最后,將不同位置處該陣列接收到的回波信號按照特定的時間延時和幅度規(guī)則進行疊加形成波束形成算法[74],從而得到目標點的聚焦和成像信號。此過程相當于把小孔徑在不同位置的信號看為一個位置固定不動的大孔徑中不同陣元接收到的信號,是將一個小孔徑合成為大孔徑的過程。

    合成孔徑效果主要是由波束合成和成像方法所決定,國內(nèi)外工作者開展了大量的研究。近年來,國外學者主要對SAFT成像算法進行了多種改進,Chang等[75]提出用平面波代替球面波的合成孔徑聚焦方法以減小超聲衍射擴展效應的影響,其橫向輻射波束寬度小于傳統(tǒng)SAFT的波束寬度;Skjelvareid等[76]提出了圓柱形掃描的SAFT算法,對于寬波束換能器可產(chǎn)生更低的旁瓣和更高的分辨率;針對超聲成像的具體對象,將SAFT和其他方法相結(jié)合會取得更好的效果。Skjelvareid等[77]將合成孔徑聚焦與虛擬源相結(jié)合擴大了換能器的掃描范圍,可應用于管道表面的腐蝕檢測并獲得較高的橫向分辨率。針對工業(yè)大型鍛件中小缺陷的檢測問題,F(xiàn)endt等[78]將SAFT和迭代反演方法相結(jié)合,可對多種小缺陷進行檢測和定位;Castano等[79]研制出適用于SAFT檢測的水射流耦合系統(tǒng),可應用于鋼軌的在線SAFT成像檢測;Lin等[80]進一步將SAFT應用到多層混凝土的檢測與成像,提高了混凝土成像檢測的準確性。

    國內(nèi)對超聲合成孔徑的研究相對落后。孫寶申等[81-83]較早研究了合成孔徑超聲成像及相關算法。近年來,國內(nèi)針對SAFT的研究也逐漸增多,但大多是跟蹤國外的研究并稍做拓展。朱新杰等[84]研制出合成孔徑超聲SH導波成像檢測系統(tǒng),并應用于鋼板的合成孔徑SH導波成像檢測。王東亞等[85]將FPGA應用于合成孔徑均勻采樣的波束合成算法,可實現(xiàn)SAFT的實時成像。杜英華等[86]發(fā)現(xiàn)多陣元非聚焦的聲場更適合于合成孔徑聚焦超聲成像,并根據(jù)這個原則進行了成像實驗,結(jié)果表明多陣元合成孔徑聚焦超聲成像可以獲得比單陣元合成孔徑聚焦成像更高的成像質(zhì)量和分辨能力。謝雪等[87]將SAFT應用于超聲衍射時差(time of fight diffraction, TOFD)檢測,有效降低了干擾衍射信號的影響。周正干等[88]將頻域合成孔徑聚焦技術和半波高法結(jié)合進行缺陷定量評價,該方法可提高缺陷檢測精度且受檢測深度的影響不明顯。富志凱等[89]提出基于幅度補償?shù)某昐AFT成像方法,與傳統(tǒng)SAFT相比其具有更高的信噪比和成像分辨率。為實現(xiàn)雙層介質(zhì)的快速超聲成像,陳堯等[90]等將虛擬源與頻域SAFT相結(jié)合得到了合成孔徑檢測圖像,可直觀顯示非規(guī)則界面形狀及內(nèi)部缺陷位置。張棣等[91]采用超聲合成孔徑技術,將輪廓線轉(zhuǎn)換為云數(shù)據(jù),并根據(jù)探頭移動方向與輪廓線的位置關系可以對物體實現(xiàn)三維成像。

    1.4 全聚焦成像

    2005年,Holmes等[92]首次提出了超聲全聚焦算法(total focusing method, TFM)的概念,該方法采用超聲相控陣換能器上的所有陣元對成像區(qū)域內(nèi)的每個擬定焦點進行聚焦[93]。首先,利用超聲換能器陣列進行全矩陣數(shù)據(jù)采集,具體方法是第一個陣元發(fā)射激勵信號所有陣元接收反射回波;然后,第二個陣元發(fā)射激勵信號所有陣元接收回波;以此類推,所有陣元依次發(fā)射激勵信號然后所有陣元接收回波信號。對于特定的焦點,計算換能器所有陣元關于該點的時間延遲,再將全矩陣回波數(shù)據(jù)按照這個延時規(guī)則進行疊加得到關于該焦點的聚焦及成像信號,當這個特定焦點遍布掃描區(qū)域時,就得到了整個掃描區(qū)域的圖像。

    與常規(guī)相控陣超聲成像方法相比,TFM可以使各陣元聲束在檢測范圍內(nèi)的每一個虛擬焦點處聚焦,具有很高的檢測精度和信噪比。但是,TFM要求每個陣元依次發(fā)射所有陣元接收,存在巨大的數(shù)據(jù)處理和運算,對實時成像提出了巨大挑戰(zhàn)。近年來,國內(nèi)外研究者在硬件和軟件上開展了大量研究工作,使TFM成像越來越成熟,越來越實用。Wilcox和Holmes等[94-95]在TFM基礎上提出了向量全聚焦成像算法(vector total focusing method, VTFM),實現(xiàn)了對缺陷類型和幾何特征的區(qū)分。考慮到固體中的波型轉(zhuǎn)換,Zhang等[96]和Budyn等[97]提出了對同一區(qū)域進行多方位多視圖的全聚焦方法(multi-view total focusing method, MTFM),MTFM可對缺陷進行分類,能區(qū)分波長尺度的體積型缺陷和裂縫型缺陷。Saini等[98]研究和測試了縱波傳播模式下TFM對表面裂紋的表征能力,通過數(shù)值計算確定了該方法測量裂紋的尺寸范圍,優(yōu)化了成像算法的可控參數(shù)。Villaverde等[99-100]將TFM與時間反轉(zhuǎn)及空間編碼相結(jié)合,提高了TFM在復合材料和粗晶等強噪聲材料中的檢測效果,較大程度地提高了信噪比。對于焊縫檢測,Sumana等[101]提出了可增加特定方向傳輸能量的角波束虛擬全聚焦方法(angle beam virtual source full matrix capture-total focusing method, ABVSFMC-TFM),該方法在鎳基合金鍛件檢測中取得了良好的應用效果。

    國內(nèi)研究者也對TFM成像方法展開了深入研究。周正干等[102]對TFM在楔塊中的能量衰減進行了校準,改善了TFM圖像中的能量均勻性,降低了缺陷的漏檢率和檢測誤差。胡宏偉等[103]研究了兩層介質(zhì)全聚焦成像時的稀疏矩陣算法,提高了傳統(tǒng)TFM的計算效率。李文濤等[104]提出了基于環(huán)形超聲陣列換能器的全聚焦方法,實現(xiàn)了鈦合金試樣中平底孔和橫孔缺陷的高精度檢測。張昊等[105]研究了斷層掃描的TFM三維成像檢測方法,得到了缺陷的三維圖像及缺陷分布情況。焦敬品等[106]采用相位特征參數(shù)對TFM成像進行幅值加權處理,發(fā)現(xiàn)相位加權的TFM成像效果明顯優(yōu)于常規(guī)TFM,在裂紋方向識別和長度測量上更加準確。楊敬等[107]提出了一種各向異性焊縫檢測的TFM,可明顯降低奧氏體鋼焊縫區(qū)域內(nèi)缺陷定位的誤差。吳斌等[108]將多模式TFM成功應用于奧氏體鋼小徑管的焊縫檢測,提高了信噪比,實現(xiàn)了小徑管的多模式TFM檢測。

    超聲相控陣、合成孔徑和全聚焦都是關于超聲陣列的成像方法,能實現(xiàn)聲束可調(diào),將各陣元信號按照一定的延時和幅度規(guī)則進行疊加處理。近年來,人們又將超聲相控陣成像與合成孔徑及全聚焦方法結(jié)合起來,在提高圖像分辨率等方面取得長足進展。Sinclair等[109]利用合成孔徑算法改善了超聲相控陣超聲檢測圖像分辨率;施克仁等[110]在提高超聲相控陣三維成像分辨率方面提出了相位畸變修正方法;Quaegebeur等[111]使用一種基于相關性的全聚焦方法,實現(xiàn)了超聲相控陣發(fā)射接收檢測數(shù)據(jù)的成像。這些研究對超聲成像技術的發(fā)展起到了推進作用,隨著電子硬件和軟件技術的發(fā)展,超聲陣列成像將朝著三維實時成像方向發(fā)展。

    2 基于頻散導波和多波聚焦的成像檢測

    本節(jié)主要分析寬帶頻散導波和多波的聚焦與成像。寬帶導波具有頻散效應,而多波是指介質(zhì)中存在縱波、橫波和導波等多種波型的傳播。寬帶頻散導波成像和多波成像具有獨特的技術優(yōu)勢,目前的研究還處在理論和實驗研究階段,有待進一步發(fā)展。

    2.1 寬帶導波的聚焦與成像

    寬帶導波聚焦是一個富有挑戰(zhàn)性的研究課題。目前,超聲導波的聚焦與成像在超聲檢測中越來越重要,已有不少的理論與應用研究報道,但基本都局限于常速度或窄帶導波的聚焦與檢測,將導波的傳播速度作為常數(shù)來處理。然而,導波的頻散效應在實際檢測中不可避免,研究寬帶頻散導波的聚焦與成像具有重要的意義。

    導波具有多模和頻散特性。針對導波多模的復雜性,通常策略是有選擇性地激發(fā)單一模式,然而不同模式適用的目標不同,因此單一模式的選擇與缺陷的識別存在矛盾。單一模式導波的選擇與控制,通常需要采用準單頻激發(fā)信號,而不是脈沖信號,還需要考慮分辨率與所需模式的特征等[112-118]。

    為了避免導波的頻散效應,一般選擇具有常速度或準常速度的模式來進行檢測[119],已有大量工作將超聲相控陣技術應用于常速度的導波進行檢測[112, 115-116,118,120]。然而將導波控制在常速度范圍內(nèi)并不容易,為克服導波的頻散特性,可通過先驗頻散曲線進行補償來消除導波的頻散效應[121],實現(xiàn)頻散導波的聚焦與成像[122-123]。

    近年來,本課題組一直致力于寬帶導波的聚焦與成像檢測研究。張碧星等[124]提出了寬帶導波相控陣的聚焦方法,實現(xiàn)了寬帶導波的聚焦與成像。首先,假定在介質(zhì)參數(shù)已知的前提下,對導波的激勵信號進行預處理,達到導波在預設焦點處聚焦的目的。在直角坐標系下,假設導波沿x方向傳播,聲源位于坐標系原點,則

    (1)

    其中:u(x,t)和u(x,ω)分別為導波在時間域和頻率域的位移分量;F(ω)是激勵聲源的頻譜函數(shù);ω是角頻率。如果導波在所考慮的頻段范圍內(nèi)具有N個模式,則

    (2)

    其中:Aj(ω)和kj分別為第j個模式的幅度和波數(shù),kj=ω/cj,cj為第j個模式的相速度,是頻率的函數(shù)。

    對于給定的介質(zhì)參數(shù),不難得到各導波模式的頻散曲線及相速度。由式(1)可知,對于不同的聲源激勵函數(shù),可得到不同的導波位移;并且式(2)表明不同模式不同頻率分量的導波具有不同的相位,從而造成散焦。對于給定的介質(zhì)參數(shù)和目標位置x,可以設計出一種聲源激勵函數(shù),使得不同模式不同頻率分量的導波都同相到達目標點實現(xiàn)聚焦。例如:將原點處的聲源頻譜函數(shù)由F(ω)變?yōu)?/p>

    (3)

    則目標點x處的導波位移分量相應地由u(x,t)變?yōu)?/p>

    (4)

    由于對u(x,ω)取絕對值,不同模式不同頻率分量的導波都將同相到達目標點,從而實現(xiàn)聚焦。

    顯然,還有多種聲源頻譜函數(shù)都能實現(xiàn)聚焦,下面給出了兩種不同頻譜函數(shù)和目標點導波位移形式:

    (5)

    (6)

    公式(5)給出的就是時間反轉(zhuǎn)法的結(jié)果,可見在寬帶導波聚焦方法中,時間反轉(zhuǎn)法是有效方法之一。事實上,時間反轉(zhuǎn)法是一種自適應聚焦法,對于非導波情況都能實現(xiàn)聚焦。

    以上方法表明,通過設計新的聲源激勵函數(shù)就能實現(xiàn)寬帶導波聚焦,不同聲源激勵函數(shù)聚焦的效果稍有不同。文獻[125]對這些方法進行了詳細研究和分析,結(jié)果表明由公式(3)和(4)給出的聚焦效果最好;時間反轉(zhuǎn)法雖然是一種很好的自適應聚焦方法,但對寬帶導波聚焦,在時間上的壓縮效果不如公式(3)。本課題組[126]還針對含低速夾層的3層半空間介質(zhì)模型,采用16陣元的超聲換能器陣列,得到半空間自由表面上的聲場分布,如圖1所示。在焦點處,各種頻率的導波同時到達,在時間上高度壓縮實現(xiàn)聚焦;而其他位置處,不同頻率的導波到達時刻不同,波形在時序上延展較長。

    圖1 寬帶導波在分層半空間表面上聚焦聲場的瞬態(tài)分布

    聚焦的導波信號遇到缺陷后會發(fā)生散射,導波散射后仍然會因為頻散效應造成散焦,因此在寬帶導波的聚焦接收與成像過程中需要對頻散效應進行補償。對于多模式的寬帶導波,本課題組提出了一種基于頻散補償?shù)娜毕荻ㄎ怀上袼惴╗127-128],通過對蘭姆波回波信號進行聚焦接收處理實現(xiàn)了板中缺陷的檢測與成像。

    頻散導波的聚焦與成像需要針對每個焦點計算每一陣元的發(fā)射信號,與普通超聲相控陣相比,計算量稍大;與聲波時間反轉(zhuǎn)成像的計算時間基本相同。

    2.2 超聲多波成像檢測

    近年來,融合多種成像方式和多波型成像是國內(nèi)外聲學成像領域的重要發(fā)展趨勢,具有代表性的成像方法有多模式全聚焦成像[93]、多波聚焦與成像[129]、多波多分量勘探[130]等。Zhang等[93]提出了多模式全聚焦方法(multi-mode total focusing method, MTFM),利用多種模式轉(zhuǎn)換的聲波對目標散射系數(shù)矩陣進行分析與預測,和常規(guī)全聚焦方法相比,可以顯著提高成像分辨率。Zhang等[129]提出了多波聚焦與成像方法,將固體介質(zhì)中的縱波、橫波和導波等多種波型同時進行聚焦與成像,和單一波型成像方法相比,不僅提高了回波信噪比,而且還增大了聲波的掃描與成像范圍。多波多分量勘探法主要應用于地球物理勘探中,是一種聚焦接收處理方法。

    本課題組針對固體中的縱波和橫波特點提出了多波聚焦的雙脈沖法[131]。目前的超聲相控陣只考慮了一種波型,因而換能器各陣元只發(fā)射一個脈沖信號,通過控制各陣元信號的時延和幅度來實現(xiàn)聚焦。在固體介質(zhì)中,由于存在縱波和橫波兩種波型,如果換能器各陣元發(fā)射2個脈沖信號,使第1個脈沖產(chǎn)生的橫波和第2個脈沖產(chǎn)生的縱波同時到達目標,那么就使縱波和橫波在焦點處實現(xiàn)了同時聚焦,即雙脈沖法的思路。當然第1個脈沖產(chǎn)生的橫波和第2個脈沖產(chǎn)生的縱波將會在焦點周圍形成旁瓣。

    雙脈沖法可以控制焦點處聲場的偏振方向。第1個脈沖產(chǎn)生的橫波S和第2個脈沖產(chǎn)生的縱波P同時到達焦點處(圖2),它們的振動方向相互垂直,合成新的質(zhì)點偏振方向M。焦點處的橫波幅度受第1個脈沖信號幅度影響較大,而縱波幅度受第2個脈沖影響較大。因而,可以通過控制雙脈沖法中兩個脈沖信號的幅度比使焦點處的質(zhì)點偏振方向沿任意方向。

    圖2 聲場偏振方向示意圖

    在接收聚焦中,也可實現(xiàn)多波聚焦。針對雙脈沖法,對換能器各陣元接收到的信號先按縱波時延進行疊加得到信號1,然后再按橫波時延進行疊加得到信號2,最后將這兩個信號1和2扣除縱波和橫波的時延差進行相加,使信號1和2中的目標信號同相相加,得到多波聚焦的成像信號。在實際操作中,不斷改變焦點位置并進行多波的聚焦發(fā)射和聚焦接收,直到焦點覆蓋整個掃描區(qū)域,就完成了整個區(qū)域的多波聚焦與成像檢測。

    采用表1所示的換能器參數(shù)以及將含有缺陷的試塊(圖3)作為被檢對象,試塊內(nèi)均勻地分布著13個直徑為2 mm的圓形通孔缺陷,分別采用縱波聚焦、橫波聚焦和多波聚焦得到了該試塊的成像結(jié)果[132],見圖4。

    表1 實驗所用相控陣換能器參數(shù)

    圖3 標準試塊結(jié)構示意圖

    圖4 不同聚焦方法的扇掃成像圖

    通過對比3種方法的成像結(jié)果可知,相比于單縱波扇掃成像多波扇掃聚焦成像結(jié)合了縱波和橫波各自的激發(fā)與傳播特性,有效提高了成像分辨率及掃描范圍。

    雙脈沖法適用于介質(zhì)中只存在縱波和橫波的情況,在含有縱波、橫波和導波的情況下,可采用時間反轉(zhuǎn)法來實現(xiàn)多波聚焦。對于給定的焦點P,假設在焦點處放置一個給定的假想聲源激勵聲場,并計算得到換能器第i個陣元接收到的信號Ai(t),將這個計算得到的信號時間反轉(zhuǎn)后(時間序列上反序)并加載到換能器對應陣元上進行激勵,按照時間反轉(zhuǎn)法的原理,這時聲場將自動聚焦在焦點位置處,即縱波、橫波和導波都在焦點處聚焦,實現(xiàn)了多波聚焦。假想聲源是為了使換能器激發(fā)的聲場能在焦點處聚焦而假設的聲源,通過假想聲源可以得到換能器各陣元的激勵信號Ai(t),在上述過程中,焦點處的假想聲源并不存在。

    在此基礎上,本課題組還提出了一種多波全聚焦的成像方法[133-134],該方法只需各陣元進行一次激勵,換能器各陣元接收到含有縱波和橫波的缺陷回波信號。將換能器陣元i發(fā)射時,陣元j接收到的信號記為Sij(t),首先將所有的Sij(t)按照縱波時延疊加得到時間信號SP(t),然后再按照橫波時延疊加得到時間信號SS(t),最后將SP(t)和SS(t)按照時間差(Δt)進行延時疊加,得到總信號,即是焦點P處的成像信號,可以表示為

    (7)

    利用公式(7)可得到多波全聚焦的成像結(jié)果。采用圖3的定制試塊進行成像檢測實驗,該試塊的材料為20#碳鋼,實驗所用的換能器參數(shù)如表1所示。對實驗數(shù)據(jù)分別進行全聚焦算法成像和多波全聚焦算法成像,結(jié)果如圖5所示。

    圖5 兩種算法的缺陷檢測成像結(jié)果

    多波全聚焦過程中計算兩次時間信號的疊加,成像時間約為單波全聚焦成像時間的2倍。但是,通過對比圖5a和圖5b可以看出,多波全聚焦成像檢測到的缺陷明顯優(yōu)于全聚焦算法成像。當偏轉(zhuǎn)角度較小時,縱波檢測起主導作用;隨著偏轉(zhuǎn)角度的增大,橫波能量逐漸增加,多波全聚焦法得到的結(jié)果信噪比增加。多波全聚焦沒有發(fā)射聚焦過程,效果沒有多波聚焦效果好,但比單個波型的全聚焦效果好。由于全聚焦的技術優(yōu)勢,多波全聚焦將會得到進一步的應用。

    以上研究和分析表明,融合多種成像方法和結(jié)合多種類型信息成像已經(jīng)成為該領域的重要發(fā)展趨勢[135],這也為超聲相控陣成像檢測提供了新的研究思路。

    3 結(jié)語

    超聲成像檢測技術作為一種常規(guī)的檢測方法,在無損檢測等眾多領域中發(fā)揮著至關重要的作用。本文對超聲成像檢測方法在無損檢測領域內(nèi)的發(fā)展歷程進行了綜述,歸納了6種超聲成像檢測方法的研究進展,對比分析了這些成像檢測方法的優(yōu)缺點及應用場景,討論了超聲成像檢測技術的發(fā)展趨勢,并指出融合多種成像方法和多信息成像將可能是今后超聲成像領域的重要發(fā)展方向。

    猜你喜歡
    信號檢測方法
    “不等式”檢測題
    “一元一次不等式”檢測題
    “一元一次不等式組”檢測題
    信號
    鴨綠江(2021年35期)2021-04-19 12:24:18
    完形填空二則
    基于FPGA的多功能信號發(fā)生器的設計
    電子制作(2018年11期)2018-08-04 03:25:42
    可能是方法不對
    小波變換在PCB缺陷檢測中的應用
    用對方法才能瘦
    Coco薇(2016年2期)2016-03-22 02:42:52
    基于LabVIEW的力加載信號采集與PID控制
    看非洲黑人一级黄片| 成人亚洲精品一区在线观看| 久久天堂一区二区三区四区| 精品国产一区二区三区四区第35| 免费少妇av软件| 欧美在线黄色| 欧美日韩福利视频一区二区| 久久人妻熟女aⅴ| 国产精品.久久久| 亚洲精品美女久久av网站| 久久精品aⅴ一区二区三区四区| 无遮挡黄片免费观看| 精品一区二区三卡| 一本一本久久a久久精品综合妖精| 亚洲国产精品999| 亚洲情色 制服丝袜| 黑人欧美特级aaaaaa片| xxxhd国产人妻xxx| 女人爽到高潮嗷嗷叫在线视频| a级毛片黄视频| 老司机靠b影院| 如何舔出高潮| 久久人人97超碰香蕉20202| 国产日韩欧美亚洲二区| 日本欧美视频一区| 人妻 亚洲 视频| 亚洲男人天堂网一区| 亚洲精品视频女| 欧美日韩亚洲国产一区二区在线观看 | 国产在线免费精品| 国产一区二区在线观看av| 日日摸夜夜添夜夜爱| 国产一级毛片在线| 日韩,欧美,国产一区二区三区| 只有这里有精品99| 男女高潮啪啪啪动态图| 久热爱精品视频在线9| 日韩一区二区视频免费看| 亚洲第一区二区三区不卡| 国产成人免费观看mmmm| 热99国产精品久久久久久7| 一个人免费看片子| av网站免费在线观看视频| 色综合欧美亚洲国产小说| 又黄又粗又硬又大视频| 久久婷婷青草| av网站免费在线观看视频| 亚洲精品国产一区二区精华液| av不卡在线播放| 在线观看人妻少妇| 满18在线观看网站| 亚洲精品国产色婷婷电影| 黄色视频在线播放观看不卡| 熟女av电影| 国产精品成人在线| 一二三四中文在线观看免费高清| videos熟女内射| 在线天堂最新版资源| 色婷婷av一区二区三区视频| 国产精品熟女久久久久浪| 日本欧美视频一区| 又大又爽又粗| 啦啦啦在线免费观看视频4| 晚上一个人看的免费电影| 在线天堂中文资源库| 在线观看免费午夜福利视频| av卡一久久| 老熟女久久久| 国产免费福利视频在线观看| 国产片内射在线| 久久久久国产精品人妻一区二区| 久久鲁丝午夜福利片| 美女主播在线视频| 99久久99久久久精品蜜桃| 亚洲国产精品国产精品| 天堂中文最新版在线下载| 777米奇影视久久| 中文乱码字字幕精品一区二区三区| 国产xxxxx性猛交| 精品福利永久在线观看| 一级毛片电影观看| 国产精品久久久久久久久免| 国产免费现黄频在线看| 中文精品一卡2卡3卡4更新| 国产老妇伦熟女老妇高清| 青春草视频在线免费观看| 久久久亚洲精品成人影院| 久久久精品94久久精品| 欧美日韩亚洲高清精品| 香蕉丝袜av| 亚洲精品乱久久久久久| 免费黄色在线免费观看| 19禁男女啪啪无遮挡网站| 亚洲欧美成人综合另类久久久| xxx大片免费视频| 夫妻午夜视频| 国产黄色免费在线视频| 精品一区二区三区四区五区乱码 | 精品人妻在线不人妻| 1024香蕉在线观看| 天天操日日干夜夜撸| 菩萨蛮人人尽说江南好唐韦庄| 欧美日韩视频高清一区二区三区二| 亚洲 欧美一区二区三区| 亚洲精品国产一区二区精华液| 久久影院123| 高清欧美精品videossex| 亚洲国产成人一精品久久久| 热99久久久久精品小说推荐| 精品国产乱码久久久久久男人| 久久性视频一级片| 好男人视频免费观看在线| av在线app专区| 亚洲七黄色美女视频| 亚洲精品国产一区二区精华液| 一级黄片播放器| 在线精品无人区一区二区三| 狂野欧美激情性xxxx| av女优亚洲男人天堂| 考比视频在线观看| 国产极品粉嫩免费观看在线| 成年人午夜在线观看视频| 亚洲人成网站在线观看播放| 一区在线观看完整版| 精品少妇内射三级| h视频一区二区三区| 欧美亚洲 丝袜 人妻 在线| 黑人欧美特级aaaaaa片| 在线观看人妻少妇| 青草久久国产| 热re99久久国产66热| 无遮挡黄片免费观看| 黄频高清免费视频| 一区在线观看完整版| 久久人妻熟女aⅴ| 热99久久久久精品小说推荐| 欧美老熟妇乱子伦牲交| av有码第一页| 午夜老司机福利片| 啦啦啦中文免费视频观看日本| 午夜91福利影院| 性色av一级| 高清在线视频一区二区三区| 亚洲第一青青草原| 国产女主播在线喷水免费视频网站| 一区二区av电影网| 国产高清国产精品国产三级| a级毛片黄视频| 欧美老熟妇乱子伦牲交| 久久精品久久久久久噜噜老黄| 国产99久久九九免费精品| 一边摸一边抽搐一进一出视频| 精品国产一区二区三区四区第35| 久久精品国产亚洲av高清一级| 日日撸夜夜添| 天堂俺去俺来也www色官网| 美女高潮到喷水免费观看| e午夜精品久久久久久久| 曰老女人黄片| 国产一区二区激情短视频 | 黄频高清免费视频| 久久人人爽av亚洲精品天堂| 国产一区有黄有色的免费视频| 久久久亚洲精品成人影院| 久久99热这里只频精品6学生| 国产在线一区二区三区精| 亚洲一码二码三码区别大吗| 久久青草综合色| 精品人妻一区二区三区麻豆| 最新的欧美精品一区二区| 日日爽夜夜爽网站| 一区二区日韩欧美中文字幕| 99热国产这里只有精品6| 国产精品熟女久久久久浪| 午夜福利,免费看| 精品国产一区二区久久| 亚洲精品成人av观看孕妇| 19禁男女啪啪无遮挡网站| 精品午夜福利在线看| 精品国产超薄肉色丝袜足j| 黑丝袜美女国产一区| 男女边吃奶边做爰视频| 国产在线一区二区三区精| 老鸭窝网址在线观看| xxxhd国产人妻xxx| av视频免费观看在线观看| 国产极品粉嫩免费观看在线| 秋霞在线观看毛片| 久久人妻熟女aⅴ| 黑人欧美特级aaaaaa片| 国产亚洲av高清不卡| 婷婷成人精品国产| 天堂8中文在线网| 国产免费视频播放在线视频| 美国免费a级毛片| 18禁观看日本| 欧美日本中文国产一区发布| 十八禁人妻一区二区| 卡戴珊不雅视频在线播放| 亚洲第一av免费看| 国产国语露脸激情在线看| a级毛片黄视频| 19禁男女啪啪无遮挡网站| xxx大片免费视频| 亚洲av综合色区一区| 在线天堂最新版资源| 久久热在线av| 国产日韩欧美在线精品| 亚洲欧美成人综合另类久久久| 日韩 亚洲 欧美在线| 久久久久久久久久久免费av| 久久午夜综合久久蜜桃| 日韩电影二区| 一级毛片 在线播放| 国产亚洲av高清不卡| 欧美日韩综合久久久久久| 国产片内射在线| 精品人妻一区二区三区麻豆| 精品国产乱码久久久久久小说| 欧美日韩av久久| 看十八女毛片水多多多| 天堂8中文在线网| 亚洲第一av免费看| 99热网站在线观看| 欧美日韩亚洲高清精品| 国产精品偷伦视频观看了| 超碰成人久久| 久久女婷五月综合色啪小说| 在线观看www视频免费| 亚洲中文av在线| 亚洲国产中文字幕在线视频| 精品亚洲成国产av| 视频区图区小说| 最近的中文字幕免费完整| 亚洲成人av在线免费| 久久午夜综合久久蜜桃| 国产成人系列免费观看| 国产精品无大码| 老司机影院成人| www.av在线官网国产| svipshipincom国产片| 伊人久久国产一区二区| 69精品国产乱码久久久| 97精品久久久久久久久久精品| 黑丝袜美女国产一区| 国产精品女同一区二区软件| 精品人妻在线不人妻| 日韩一卡2卡3卡4卡2021年| 青春草国产在线视频| 狠狠婷婷综合久久久久久88av| 欧美日韩视频高清一区二区三区二| 日韩视频在线欧美| 最新在线观看一区二区三区 | 精品少妇一区二区三区视频日本电影 | av免费观看日本| 一边摸一边抽搐一进一出视频| 久久久久久久久久久免费av| 亚洲 欧美一区二区三区| 80岁老熟妇乱子伦牲交| 在线观看三级黄色| 亚洲精品一区蜜桃| 国产国语露脸激情在线看| 麻豆精品久久久久久蜜桃| 男女之事视频高清在线观看 | 国产老妇伦熟女老妇高清| 久久久久久久久久久免费av| 啦啦啦在线免费观看视频4| kizo精华| 丰满乱子伦码专区| 男女无遮挡免费网站观看| 久久久久精品久久久久真实原创| 国产精品久久久久久久久免| 日本av免费视频播放| 亚洲欧美一区二区三区国产| 满18在线观看网站| 热99国产精品久久久久久7| 成人亚洲精品一区在线观看| 中文字幕av电影在线播放| a级片在线免费高清观看视频| 18禁裸乳无遮挡动漫免费视频| 老汉色∧v一级毛片| 麻豆av在线久日| av不卡在线播放| 午夜福利在线免费观看网站| 国产免费视频播放在线视频| 岛国毛片在线播放| 成人漫画全彩无遮挡| 男女无遮挡免费网站观看| 国产成人免费观看mmmm| 亚洲国产欧美日韩在线播放| 成人黄色视频免费在线看| 国产精品一区二区精品视频观看| 交换朋友夫妻互换小说| 曰老女人黄片| 亚洲av欧美aⅴ国产| 午夜福利影视在线免费观看| 99精国产麻豆久久婷婷| 亚洲国产最新在线播放| 女人久久www免费人成看片| 黄片播放在线免费| 91成人精品电影| 99国产综合亚洲精品| 欧美黑人精品巨大| 99国产综合亚洲精品| 波多野结衣一区麻豆| 久久99热这里只频精品6学生| 国产成人系列免费观看| 国产黄色免费在线视频| 久久影院123| 国产成人欧美| 叶爱在线成人免费视频播放| 2018国产大陆天天弄谢| 亚洲精品在线美女| 国产精品一区二区在线观看99| 在线观看三级黄色| 丝袜在线中文字幕| 丝袜美腿诱惑在线| 久久久久精品国产欧美久久久 | 九草在线视频观看| 亚洲成国产人片在线观看| 亚洲五月色婷婷综合| 岛国毛片在线播放| 色综合欧美亚洲国产小说| 蜜桃国产av成人99| 最近中文字幕2019免费版| 成人影院久久| 日日啪夜夜爽| 18在线观看网站| a 毛片基地| 精品免费久久久久久久清纯 | 一级爰片在线观看| 十分钟在线观看高清视频www| 亚洲图色成人| 桃花免费在线播放| 亚洲国产毛片av蜜桃av| 亚洲国产欧美在线一区| 精品视频人人做人人爽| 2018国产大陆天天弄谢| 成人午夜精彩视频在线观看| 国产高清国产精品国产三级| 国产成人一区二区在线| 国产成人av激情在线播放| 亚洲欧美成人精品一区二区| 国产高清国产精品国产三级| 午夜福利,免费看| 中国国产av一级| 国产成人午夜福利电影在线观看| 久久影院123| 亚洲精品一二三| 久久ye,这里只有精品| 精品人妻在线不人妻| www.精华液| 精品一品国产午夜福利视频| svipshipincom国产片| 一区二区日韩欧美中文字幕| 成人影院久久| 亚洲成av片中文字幕在线观看| 两个人看的免费小视频| 深夜精品福利| 久久精品久久精品一区二区三区| 国产一区二区三区综合在线观看| 久久 成人 亚洲| 黄色一级大片看看| 精品少妇内射三级| 在线天堂最新版资源| 国产在视频线精品| 国产伦人伦偷精品视频| 欧美国产精品va在线观看不卡| 亚洲欧美精品自产自拍| 水蜜桃什么品种好| 欧美日韩成人在线一区二区| 婷婷色av中文字幕| 别揉我奶头~嗯~啊~动态视频 | 国产精品一二三区在线看| 免费人妻精品一区二区三区视频| 蜜桃国产av成人99| 亚洲欧洲国产日韩| 老熟女久久久| 成年女人毛片免费观看观看9 | 精品少妇一区二区三区视频日本电影 | 午夜免费观看性视频| 久久精品久久久久久噜噜老黄| a 毛片基地| 啦啦啦 在线观看视频| 精品久久蜜臀av无| 老司机靠b影院| 香蕉丝袜av| 中文字幕最新亚洲高清| 欧美日韩福利视频一区二区| 最近手机中文字幕大全| 久久久精品免费免费高清| 一区在线观看完整版| 卡戴珊不雅视频在线播放| 精品人妻在线不人妻| 精品一区二区三区四区五区乱码 | 看非洲黑人一级黄片| 午夜91福利影院| av视频免费观看在线观看| 国产一区二区在线观看av| 久久久久久免费高清国产稀缺| 制服丝袜香蕉在线| 亚洲色图综合在线观看| 一区二区日韩欧美中文字幕| 精品一区二区三区四区五区乱码 | 精品亚洲成国产av| 欧美日韩亚洲国产一区二区在线观看 | 建设人人有责人人尽责人人享有的| 色精品久久人妻99蜜桃| 国产精品国产三级国产专区5o| 只有这里有精品99| 亚洲精品国产av蜜桃| 日韩制服骚丝袜av| 巨乳人妻的诱惑在线观看| www.精华液| 1024香蕉在线观看| 日日摸夜夜添夜夜爱| 成人亚洲欧美一区二区av| 夜夜骑夜夜射夜夜干| 成人毛片60女人毛片免费| 国产精品一区二区精品视频观看| 91aial.com中文字幕在线观看| 中文字幕制服av| 极品人妻少妇av视频| 亚洲成人免费av在线播放| 狂野欧美激情性xxxx| 日韩一区二区视频免费看| 亚洲av电影在线观看一区二区三区| 国产成人一区二区在线| 午夜日韩欧美国产| 国产成人精品无人区| 亚洲精品成人av观看孕妇| 亚洲熟女精品中文字幕| 999精品在线视频| 婷婷色麻豆天堂久久| 欧美在线黄色| 可以免费在线观看a视频的电影网站 | 国产精品久久久久久精品电影小说| 一区二区三区激情视频| 国产日韩欧美在线精品| 午夜福利免费观看在线| 国产精品成人在线| 99久久精品国产亚洲精品| 777米奇影视久久| 亚洲国产欧美网| 黑人巨大精品欧美一区二区蜜桃| 深夜精品福利| 亚洲激情五月婷婷啪啪| 午夜免费鲁丝| 欧美亚洲日本最大视频资源| 国产成人午夜福利电影在线观看| 亚洲精品aⅴ在线观看| xxx大片免费视频| 无限看片的www在线观看| 久久人人97超碰香蕉20202| 欧美日本中文国产一区发布| 日本av免费视频播放| 免费少妇av软件| 亚洲精品久久成人aⅴ小说| 男女之事视频高清在线观看 | 亚洲av在线观看美女高潮| 精品国产露脸久久av麻豆| 精品一区二区三区四区五区乱码 | 精品国产一区二区三区久久久樱花| 国产一区亚洲一区在线观看| 9热在线视频观看99| 久久婷婷青草| 亚洲一区二区三区欧美精品| av一本久久久久| 老司机靠b影院| 亚洲av日韩精品久久久久久密 | 欧美老熟妇乱子伦牲交| 69精品国产乱码久久久| 亚洲婷婷狠狠爱综合网| 国产亚洲av片在线观看秒播厂| 最近中文字幕2019免费版| 各种免费的搞黄视频| 最近的中文字幕免费完整| 看免费成人av毛片| 十分钟在线观看高清视频www| 丝袜脚勾引网站| 啦啦啦 在线观看视频| 精品一区二区三区四区五区乱码 | 如日韩欧美国产精品一区二区三区| 高清不卡的av网站| 成年动漫av网址| 亚洲熟女精品中文字幕| 国产免费视频播放在线视频| 女人高潮潮喷娇喘18禁视频| 一级,二级,三级黄色视频| av国产精品久久久久影院| 伊人久久国产一区二区| 成人毛片60女人毛片免费| 人人澡人人妻人| 97在线人人人人妻| 久久久久久久久久久久大奶| 久久久久久久精品精品| 亚洲国产毛片av蜜桃av| 久久精品亚洲av国产电影网| 久久精品国产综合久久久| 日韩视频在线欧美| 欧美日韩视频精品一区| 日韩欧美精品免费久久| 欧美变态另类bdsm刘玥| 国产成人91sexporn| 水蜜桃什么品种好| av有码第一页| 日韩av免费高清视频| 啦啦啦啦在线视频资源| 在线天堂最新版资源| www.精华液| 嫩草影视91久久| 亚洲一区二区三区欧美精品| 久久综合国产亚洲精品| 乱人伦中国视频| 午夜日本视频在线| 91精品伊人久久大香线蕉| 考比视频在线观看| 亚洲美女黄色视频免费看| av电影中文网址| 亚洲国产精品999| 狂野欧美激情性xxxx| 国产一级毛片在线| 一级片免费观看大全| 一本大道久久a久久精品| 午夜福利在线免费观看网站| 久久97久久精品| 欧美激情 高清一区二区三区| 亚洲成色77777| 国产精品久久久久久人妻精品电影 | 国产av码专区亚洲av| 免费看不卡的av| 搡老岳熟女国产| 日本91视频免费播放| av线在线观看网站| 欧美av亚洲av综合av国产av | 少妇 在线观看| 欧美最新免费一区二区三区| 久久久国产一区二区| 亚洲精品国产区一区二| 国产有黄有色有爽视频| 我的亚洲天堂| 国产成人系列免费观看| 日本色播在线视频| 久热爱精品视频在线9| 国产成人精品在线电影| 免费av中文字幕在线| 涩涩av久久男人的天堂| 国产精品无大码| 免费在线观看视频国产中文字幕亚洲 | 亚洲精品久久成人aⅴ小说| 哪个播放器可以免费观看大片| 又黄又粗又硬又大视频| 亚洲在久久综合| av网站在线播放免费| 免费黄网站久久成人精品| 满18在线观看网站| 精品一区二区三区av网在线观看 | 国产男女内射视频| 亚洲精品中文字幕在线视频| 亚洲av欧美aⅴ国产| 91精品伊人久久大香线蕉| 久久毛片免费看一区二区三区| 少妇的丰满在线观看| 久久久久久久国产电影| 欧美精品一区二区大全| 最近中文字幕高清免费大全6| 欧美精品高潮呻吟av久久| 日本色播在线视频| 亚洲欧美日韩另类电影网站| 性高湖久久久久久久久免费观看| 久久久国产一区二区| 亚洲美女搞黄在线观看| 在线观看免费日韩欧美大片| 2021少妇久久久久久久久久久| 国产精品二区激情视频| 亚洲国产中文字幕在线视频| 亚洲伊人色综图| 成人手机av| 亚洲欧美激情在线| 亚洲精品aⅴ在线观看| 国产免费一区二区三区四区乱码| 久久鲁丝午夜福利片| 在线观看免费高清a一片| 欧美日韩亚洲综合一区二区三区_| 99九九在线精品视频| 国产有黄有色有爽视频| 免费日韩欧美在线观看| 亚洲欧美一区二区三区黑人| 国产精品嫩草影院av在线观看| 人成视频在线观看免费观看| 老司机在亚洲福利影院| 日韩一本色道免费dvd| 亚洲熟女毛片儿| 亚洲天堂av无毛| 久久久久久人妻| 女人精品久久久久毛片| 日韩 欧美 亚洲 中文字幕| 丰满乱子伦码专区| 久久精品国产亚洲av涩爱| 少妇猛男粗大的猛烈进出视频| 成人亚洲欧美一区二区av| 亚洲精品日本国产第一区| 母亲3免费完整高清在线观看| 国产日韩欧美视频二区| 大香蕉久久成人网| 丝袜在线中文字幕| 亚洲欧洲国产日韩| 国产亚洲精品第一综合不卡| 久久久久国产一级毛片高清牌| 午夜福利乱码中文字幕| 久久影院123| 亚洲精品美女久久av网站| 中文字幕高清在线视频| 一级毛片黄色毛片免费观看视频| a级毛片黄视频|