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

    儲罐底板鋼點蝕過程中聲發(fā)射信號的聚類分析

    2015-11-22 06:54:02畢海勝李自力胡德棟ISAACTokuGyamerah吳向陽
    關(guān)鍵詞:信號

    畢海勝,李自力,胡德棟,羅 琴,ISAAC Toku-Gyamerah,吳向陽

    (1.青島科技大學(xué)機電工程學(xué)院,山東青島266061;2.中國石油大學(xué)儲運與建筑工程學(xué)院,山東青島266580;3.中石化石油工程設(shè)計有限公司數(shù)字化工程設(shè)計所,山東東營257000;4.延長油田股份有限公司杏子川采油廠勘探開發(fā)研究所,陜西延安717400)

    儲罐底板鋼點蝕過程中聲發(fā)射信號的聚類分析

    畢海勝1,李自力2,胡德棟1,羅 琴3,ISAAC Toku-Gyamerah2,吳向陽4

    (1.青島科技大學(xué)機電工程學(xué)院,山東青島266061;2.中國石油大學(xué)儲運與建筑工程學(xué)院,山東青島266580;3.中石化石油工程設(shè)計有限公司數(shù)字化工程設(shè)計所,山東東營257000;4.延長油田股份有限公司杏子川采油廠勘探開發(fā)研究所,陜西延安717400)

    聯(lián)合采用聲發(fā)射和電化學(xué)技術(shù)研究儲罐底板鋼試樣在w(NaCl)=3.0%,pH=2.0的酸性溶液中的點蝕特征,基于K-均值聚類算法對點蝕聲發(fā)射信號特征參數(shù)進行聚類分析,從而提取各類信號的自身特征。將分類后的信號作為樣本訓(xùn)練BP人工神經(jīng)網(wǎng)絡(luò),成功對平行試驗采集的聲發(fā)射信號進行識別。研究結(jié)果表明,底板鋼在酸性條件下的點蝕過程主要產(chǎn)生氫氣泡、膜破裂和蝕坑生長這3類典型的聲發(fā)射信號,通過聚類方法可以區(qū)分這3類信號,并能用神經(jīng)網(wǎng)絡(luò)對聲發(fā)射源進行有效識別。這對現(xiàn)場常壓金屬儲罐底板腐蝕聲發(fā)射檢測結(jié)果的解釋和評價具有指導(dǎo)意義,有助于提高檢測結(jié)果可靠性,降低儲罐運行風(fēng)險,保證其運行安全。

    底板鋼;點蝕;聲發(fā)射;K-均值聚類;Gabor小波變換

    Q235低碳鋼在石油化工行業(yè)應(yīng)用廣泛,是國內(nèi)常壓金屬儲罐底板普遍使用的材料之一,常壓金屬儲罐在原油、成品油以及石化產(chǎn)品的儲運過程中起著不可替代作用[1]。然而超過80%的儲罐關(guān)閉、穿孔和泄漏事故都是由于儲罐底板腐蝕引起的,底板的腐蝕穿孔嚴重威脅到儲罐的運行安全。聲發(fā)射作為一種動態(tài)在線檢測技術(shù),相對于傳統(tǒng)離線無損檢測技術(shù),具有低成本、快速高效全覆蓋檢測,能夠?qū)顒有匀毕葸M行在線監(jiān)控和早期破壞預(yù)警的特點,因而在常壓金屬儲罐底板腐蝕檢測方面得到了推廣應(yīng)用[2]。聲發(fā)射是材料中局域源能量快速釋放而產(chǎn)生瞬態(tài)彈性波的現(xiàn)象[3],目前被廣泛應(yīng)用于材料屈服變形、疲勞斷裂、腐蝕損傷和表面摩擦等過程的在線監(jiān)測。Ferrer等[4-5]研究證明,腐蝕過程產(chǎn)生的聲發(fā)射信號與腐蝕因素和電化學(xué)機理密切相關(guān)。Fregonese等[6]對AISI 316L不銹鋼點蝕的誘導(dǎo)期和擴展期進行了研究,認為點蝕擴展階段能夠產(chǎn)生豐富的聲發(fā)射信號;Darowicki等[7]從概率統(tǒng)計角度將點蝕累積分布函數(shù)與聲發(fā)射撞擊累積分布函數(shù)進行了關(guān)聯(lián);Prateepasen等研究表明[8-9],點蝕過程氫氣泡破裂產(chǎn)生的聲發(fā)射頻率高于125 kHz,并提出了氫氣泡破裂頻率的計算公式。金屬自身的變形和斷裂屬于直接聲發(fā)射源,而在腐蝕過程中氫氣泡的活動、鈍化膜的破裂、點蝕坑的生長和擴展、腐蝕產(chǎn)物層堆積、摩擦和開裂等均屬于二次聲發(fā)射源[10]。受聲發(fā)射源自身特性、聲信號到傳感器的傳播途徑、傳感器的特性和聲發(fā)射儀器測量系統(tǒng)等多種因素的影響,接收到的聲發(fā)射信號和聲發(fā)射源信號之間存在較大差距,因此如何對聲發(fā)射信號進行分析和處理以獲取有用的聲發(fā)射源信息成為聲發(fā)射檢測技術(shù)應(yīng)用的瓶頸[11]。對于常壓金屬儲罐底板腐蝕聲發(fā)射檢測來說,通過單一聲發(fā)射參數(shù)分析得出的結(jié)果往往與現(xiàn)場實際差距較大,提取腐蝕聲發(fā)射信號特征以評定腐蝕聲發(fā)射源的特性是現(xiàn)場檢測的難題。基于此筆者通過聲發(fā)射技術(shù)并結(jié)合電化學(xué)“原位”測量技術(shù)對底板鋼點蝕進行監(jiān)測,采用Gabor小波變換提取各類聲發(fā)射源信號的時頻域局部特征,進一步利用BP神經(jīng)網(wǎng)絡(luò)對各類聲發(fā)射源進行判斷和識別。

    1 試驗材料和方法

    1.1 試驗材料

    試驗材料為Q235鋼,并采用線切割加工成正方形試樣(29 mm×29 mm,厚度(2±0.1)mm)。試樣先采用320~1 200#金相砂紙依次打磨,去離子水清洗表面雜質(zhì),再用丙酮除脂,酒精棉擦拭后放入干燥器,干燥稱重,標號備用。

    1.2 試驗裝置

    試驗裝置如圖1所示。試驗容器為250 mm× 250 mm×250 mm的正方體容器,容器頂部不完全密封,容器底部中央開有階梯方孔,分別為50 mm×50 mm×5 mm的外孔和32 mm×32 mm×5 mm的內(nèi)孔。容器全部采用有機玻璃加工而成,目的是避免非試驗部分腐蝕產(chǎn)生聲發(fā)射信號對試樣腐蝕的干擾。將封裝好的試樣嵌入到容器底部的中心階梯方孔,四周邊緣處采用硅酮密封膠793密封。聲發(fā)射傳感器通過凡士林與試樣背面進行良好的聲耦合,并通過同軸電纜與聲發(fā)射儀的采集通道相連接,實時監(jiān)測來自試樣工作面的腐蝕信號,工作面直接與腐蝕溶液接觸。容器與基座、基座與試驗臺之間均襯有5 mm厚的橡膠墊以減少外界振動干擾。為驗證聲發(fā)射測試結(jié)果,同時在容器內(nèi)構(gòu)造了標準三電極體系,并與電化學(xué)工作站相連接,同步進行電化學(xué)測試。試驗在室溫下進行。

    圖1 試驗系統(tǒng)示意圖Fig.1 Schematic diagram of experiment system

    1.3 測試儀器

    1.3.1 聲發(fā)射測試系統(tǒng)

    采用德國華倫公司(Vallen Systeme)研發(fā)的AMSY-5第二代全數(shù)字、全波形、強抗干擾聲發(fā)射儀。聲發(fā)射傳感器型號為VS150-RIC,頻率為100~450 kHz,中心頻率為150 kHz,集成前置放大34 dB,采集閾值30.2 dB(ref.1 μV),采樣頻率為2 MHz,硬件濾波器帶寬25~500 kHz,靈敏度-57.5 dB(1 V/0.1 Pa)。

    1.3.2 電化學(xué)測試系統(tǒng)

    采用美國AMETEX公司生產(chǎn)的PAR 2273電化學(xué)工作站和配套的PowerSuite軟件構(gòu)成電化學(xué)測試和分析系統(tǒng)。該系統(tǒng)的電流最小分辨率為1.2 fA,最小電位步長2.5 μV。試驗采用標準三電極體系,工作電極(WE)為封裝好的底板鋼試樣,參比電極(RE)體系由飽和甘汞電極(SCE)和直鹽橋組成,輔助電極(CE)為鉑片電極(Pt)。

    1.3.3 表面形貌觀察

    采用德國Zeiss體式顯微鏡和Zeiss Axio 3D共聚焦金相顯微鏡分別對腐蝕后的試樣進行表面形貌和三維蝕坑形貌觀察,并生成典型蝕坑的三維圖像。配置有micro-image analysis&process(Miaps)金相圖像分析軟件,可以對蝕坑數(shù)量、面積和最大蝕坑深度進行統(tǒng)計分析。

    2 試驗結(jié)果與討論

    向電解池中充入腐蝕溶液至80%液位,并靜置

    2 h待溶液穩(wěn)定。先前通過動電位極化測得該條件下底板鋼的點蝕電位在0.8~1.0 V,因此通過設(shè)定電位0.8 V,連續(xù)恒電位極化3 h來加速點蝕。試樣在酸性條件下主要發(fā)生以下電化學(xué)反應(yīng):

    加速點蝕過程中,氫氣泡的振蕩、運動和崩潰,氯離子侵蝕作用下氧化膜的剝離和破裂,以及腐蝕產(chǎn)物的堆積和摩擦都產(chǎn)生了明顯的聲發(fā)射信號。

    2.1 聲發(fā)射特征參數(shù)提取和相關(guān)分析

    2.1.1 特征參數(shù)提取

    由于直接采用波形或頻譜進行聚類分析會使計算量很大,且單一的頻譜很難描述信號之間的差異,因此必須對聲發(fā)射波形進行特征提取。典型的聲發(fā)射特征參數(shù)包括幅度A、上升時間R、能量E、持續(xù)時間D、振鈴計數(shù)C,除了這5個特征參數(shù)外,還另外引進了3個衍生的聲發(fā)射參數(shù)[12]:

    (Ⅰ)RA值:上升時間與幅度的比值,表征信號的斜率,單位為ms V-1。

    (Ⅱ)AF值:振鈴數(shù)與持續(xù)時間的比值,表征信號的平均頻率,單位為kHz。

    (Ⅲ)RD值:上升時間與持續(xù)時間之比,是一個無量綱數(shù),表征信號波形的陡峭程度。

    特征參數(shù)提取后,為了避免具有不同物理意義和量綱的輸入變量不平等使用,消除指標之間的量綱影響,對特征向量原始數(shù)據(jù)進行Z-score標準化處理。

    2.1.2 特征參數(shù)相關(guān)分析

    在神經(jīng)網(wǎng)絡(luò)模式識別分析中,要求兩兩變量之間盡可能是相關(guān)性不大的獨立變量,這樣可避免出現(xiàn)向量“維數(shù)災(zāi)難”[13]。在聚類和神經(jīng)網(wǎng)絡(luò)分析之前,須對聲發(fā)射各特征參數(shù)之間進行相關(guān)分析。判斷兩個變量x和y之間的相關(guān)程度,通常采用皮爾遜相關(guān)系數(shù)r。

    其中

    對試驗采集到的聲發(fā)射信號8個特征參數(shù)A、R、E、D、C、RA值、AF值和RD值進行皮爾遜相關(guān)分析,相關(guān)系數(shù)矩陣見表1。

    表1 皮爾遜相關(guān)系數(shù)矩陣Table 1 Pearson correlation coefficients matrix

    從相關(guān)系數(shù)矩陣可看出,RA值與R相關(guān)系數(shù)超過了0.9,近似線性相關(guān),同時使用這兩個特征參數(shù)并不能提供額外的信息,且會增加聚類分析的計算量,更會對后面的神經(jīng)網(wǎng)絡(luò)分析產(chǎn)生不利的影響。同時R受采集閾值影響較大,不及RA值包含有更多的聲發(fā)射信號信息。本文中聚類分析選用A、E、D、C、RA、AF和RD這7個特征參數(shù)。

    2.2 基于K-均值的聚類分析

    K-均值聚類法,又稱快速聚類法,是由James MacQueen于1967年提出的,它將數(shù)據(jù)看成K維空間上的點,以距離作為測度個體“親疏程度”的指標。在聚類過程中,樣本所屬的類會不斷調(diào)整,目的是使各個樣本與所在類均值的誤差平方和達到最小。假設(shè)給定一個待分類數(shù)據(jù)集,把它分屬于k個子集C1,C2,…,Ck,滿足

    當達到設(shè)定的預(yù)定誤差時,停止迭代。其中mj是類別Cj的聚類中心[14-17]。

    試驗觀察到,在恒電位極化的前250 s時段,電流密度比較大,氫離子還原產(chǎn)生了大量的氫氣泡,并且在100~230 s觀察到試樣表面生成的氧化膜開始剝離、皺縮和破裂,氫氣泡潰滅和氧化膜破裂密集產(chǎn)生了大量的聲發(fā)射信號,初步確定該時段K均值聚類數(shù),K=2,并進行聚類計算。圖2為250 s之前聲發(fā)射信號幅度隨時間的分布。從圖2可看出,聚類1中聲發(fā)射信號幅值相對較低,分布在32~45 dB范圍內(nèi),聚類2中信號主要集中在45~60 dB較高幅值區(qū)間。

    圖2 聲發(fā)射信號幅度隨時間分布(250 s之前)Fig.2 AE signal amplitude distribution with time(before 250 s)

    圖3 為聲發(fā)射信號振鈴數(shù)、能量與持續(xù)時間的關(guān)聯(lián)關(guān)系。由圖3可以看出,聚類算法將聲發(fā)射信號很好地劃分為2類。聚類1中信號振鈴數(shù)基本小于25,能量值集中在100 eu(1 eu=10-14V2s)以下,持續(xù)時間很少超過300 μs;聚類2中信號振鈴數(shù)主要分布在25~150較寬范圍內(nèi),且能量高達1 500 eu,持續(xù)時間相對較長。

    圖3 聲發(fā)射信號振鈴數(shù)、能量與持續(xù)時間的關(guān)聯(lián)關(guān)系(250 s之前)Fig.3 AE signal count-duration correlations and energy-duration correlations(before 250 s)

    在氧化膜局部剝離破裂的同時,點蝕開始進入誘導(dǎo)期。經(jīng)過一個短暫的孕育之后,大約在1 000 s之后,點蝕進入穩(wěn)定發(fā)展階段[18]。此階段聲發(fā)射強度依然很高,聲發(fā)射信號主要來自于氫氣泡活動、點蝕生長擴展和腐蝕產(chǎn)物剝離與開裂,初步確定聚類數(shù)K=3,對點蝕穩(wěn)定發(fā)展階段的聲發(fā)射信號進行聚類分析,如圖4、5所示。

    圖4為AE信號能量與持續(xù)時間分布。由圖4可以看出,聚類算法將3類信號很好地區(qū)分開來。與之前的250 s時段的AE信號相比較,出現(xiàn)了大量的新聚類3信號,其持續(xù)時間集中在100 μs以下,同時聚類1信號急劇減少。從圖5可看出,主要是聚類2和聚類3這兩類信號分布于整個采樣時間,而新出現(xiàn)的聚類3信號占絕對主導(dǎo)。

    在恒電位極化進行1000 s之后,觀察到氫氣泡的活動顯著減弱,電流密度趨近穩(wěn)定,試樣發(fā)生穩(wěn)定的點蝕,且表面腐蝕產(chǎn)物開始逐層堆積并出現(xiàn)開裂。綜上分析初步推斷,聚類1主要由氫氣泡產(chǎn)生,聚類2主要由前期氧化膜破裂和后期腐蝕產(chǎn)物活動所產(chǎn)生,而聚類3主要由后期大量點蝕擴展產(chǎn)生。

    圖4 聲發(fā)射信號能量與持續(xù)時間的關(guān)聯(lián)關(guān)系(1 000 s之后)Fig.4 AE signal energy-duration correlations(after 1 000 s)

    圖5 聲發(fā)射信號幅值隨時間分布(1000 s之后)Fig.5 AE signal amplitude distribution with time(after 1 000 s)

    為了進一步證實上述推斷,采用Gabor小波變換(Gabor wavelet transform)對這3類聲發(fā)射信號進行時頻局域化分析,分別提取其頻帶能量特征[19-21],如圖6、7、8所示,圖中顏色標尺表示信號在某個時間點某個頻率分量的能量大小,紅色為最大,粉紅為最小。

    從圖6可看出,氫氣泡AE信號近似于諧振信號[4],這與氣泡在溶液中的振蕩和空化過程密切相關(guān)。從Gabor小波變換時頻圖上看,信號的能量集中度高,中心能量對應(yīng)的頻帶范圍為130~170 kHz,峰值頻率為160 kHz。Minnaert提出了計算單氣泡聲共振諧振頻率的Minnaert公式[22-25],而Prateepasen A等[8,26-28]在考慮到單個穩(wěn)定氫氣泡的最大直徑Dmax和破裂瞬間壓力pb之后,將公式修正為

    Prateepasen等深入研究后建議pb取3.110×107Pa。試驗中利用高速攝像機捕捉到最大穩(wěn)定氫氣泡直徑約為0.9 mm,取γ=1.4,ρ0=1.0×103kg·m-3,代入式(7)可計算出氫氣泡破裂的頻率fb≥127.8 kHz,這與上述的氫氣泡信號的頻帶范圍基本吻合。

    圖6 氫氣泡AE信號波形及其小波變換時頻圖Fig.6 Waveform and time-frequency plane of AE signal of hydrogen bubble

    圖7 氧化膜破裂AE信號波形及其小波變換時頻圖Fig.7 Waveform and time-frequency plane of AE signal of corrosion film breakage

    相對于氫氣泡AE信號,膜破裂AE信號能量分布范圍相對寬一些(圖7),能量集中在90~200 kHz,膜破裂的微觀過程復(fù)雜,信號所包含的頻率成分較豐富,最大能量對應(yīng)的時刻為65 μs,相對滯后;而點蝕AE信號能量發(fā)散較小(圖8),頻率集中在95~130 kHz,最大能量出現(xiàn)的時刻更提前,在15 μs左右。這與Prateepasen等研究不銹鋼點蝕時所獲得的氣泡和點蝕擴展信號的頻帶特征基本相符,進一步驗證對3類聲發(fā)射源的推斷是正確的。

    圖8 點蝕生長AE信號波形及其小波變換時頻圖Fig.8 Waveform and time-frequency plane of AE signal of pitting growth

    2.3 BP人工神經(jīng)網(wǎng)絡(luò)識別

    為了驗證3種聚類信號對應(yīng)的3種不同的腐蝕聲發(fā)射源的準確性,又做了多組平行試驗。以A、E、D、C、RA、AF和RD值這7個特征參數(shù)作為BP神經(jīng)網(wǎng)絡(luò)的輸入向量,3類聲發(fā)射源為輸出向量,隱含層節(jié)點數(shù)為6,建立一個7-6-3單隱含層BP神經(jīng)網(wǎng)絡(luò)[29-31]。隨機抽取已分出來的3類聲發(fā)射源信號各1000組數(shù)據(jù)中2 250組數(shù)據(jù)作為訓(xùn)練樣本,剩余750組作為測試樣本,每組抽取的各類聲發(fā)射源樣本個數(shù)不相等。利用已訓(xùn)練好的神經(jīng)網(wǎng)絡(luò)來對平行試驗采集到的聲發(fā)射信號進行識別。識別出來的各類信號其能量與持續(xù)時間分布如圖9所示。

    從圖9可以看出,平行試驗采集的AE信號能量與持續(xù)時間也存在3個明顯的分區(qū),且正確對應(yīng)3類聲發(fā)射源。由于采用的是恒電位極化最后1 h的聲發(fā)射數(shù)據(jù),這期間AE信號主要來自點蝕擴展和腐蝕產(chǎn)物剝離、開裂,氫氣泡AE信號相對較少。占主導(dǎo)的點蝕AE信號集中分布在低持續(xù)時間、低能量區(qū),腐蝕產(chǎn)物開裂AE信號的持續(xù)時間較前者高,能量分布相對發(fā)散。試驗后采用3D共聚焦金相顯微鏡對試樣表面形貌和三維蝕坑形貌進行觀察,如圖10所示。試樣表面出現(xiàn)了大量隨機分布的、呈蜂窩狀、大小不一的蝕坑,經(jīng)軟件統(tǒng)計分析,蝕坑平均深度達195 μm,證實試樣表面確實發(fā)生了明顯的點蝕。

    圖9 聲發(fā)射信號能量與持續(xù)時間的關(guān)聯(lián)Fig.9 AE signal energy-duration correlations

    圖10 試樣點蝕后表面形貌Fig.10 Surface morphology of specimen after pitting corrosion

    3 結(jié) 論

    (1)底板鋼試樣在酸性NaCl溶液發(fā)生點蝕的過程中主要存在3類聲發(fā)射源:氫氣泡、膜剝離和破裂、點蝕坑生長,K-means聚類算法可將這3類典型AE信號很好地區(qū)分開來。

    (2)在點蝕進入穩(wěn)定階段后,蝕坑的生長、擴展和腐蝕產(chǎn)物的剝離、開裂是主要的AE信號,氫氣泡AE信號相對減弱,這與電化學(xué)測試結(jié)果相一致。

    (3)BP人工神經(jīng)網(wǎng)絡(luò)能較好地識別出這3類聲發(fā)射源,通過提取各類源信號特征,并與現(xiàn)場聲發(fā)射檢測結(jié)果相結(jié)合,考慮實際AE信號的衰減并作適當修正,這對常壓金屬儲罐底板腐蝕聲發(fā)射檢測結(jié)果的解釋和評價具有指導(dǎo)意義。

    [1] 閻永貴,吳建華,陳光章.原油儲罐內(nèi)底板的腐蝕防護現(xiàn)狀與展望[J].腐蝕與防護,2002,23(5):199-201. YAN Yonggui,WU Jianhua,CHEN Guangzhang.Recent development and prospect corrosion and protection of inside bottom of crude oil tanks[J].Corrosion&Protection,2002,23(5):199-201.

    [2] CHEN Jianfei,BI Haisheng,WANG Qiang.The application of acoustic emission technology in oil and gas storage and transportation equipment[J].Advanced Materials Research,2013,807:2652-2657.

    [3] ASME.Standard terminology for nondestructive examinations:E1316-06a[S].United States:ASTM,2006.

    [4] FERRER F,F(xiàn)AURE T,GOUDIAKAS J,et al.Acoustic emission study of active-passive transitions during carbon steel erosion-corrosion in concentrated sulfuric acid[J]. Corrosion Science,2002,44(7):1529-1540.

    [5] FERRER F,IDRISSI H,MAZILLE H.A study of abrasion corrosion of AISI 304L austenitic stainless steel in saline solution using acoustic emission technique[J]. NDT&E International,2000,33(6):363-371.

    [6] FREGONESE M,IDRISSI H,MAZILLE H,et al.Initiation and propagation steps in pitting corrosion of austenitic stainless steels:monitoring by acoustic emission[J].Corrosion Science,2001,43(4):627-641.

    [7] DAROWICKI K,MIRAKOWSKI A,KRAKOWIAK S. Investigation of pitting corrosion of stainless steel by means of acoustic emission and potential dynamic methods[J].Corrosion Science,2003,45(8):1747-1756.

    [8] PRATEEPASEN A,JIRARUNGSATIAN C.Implementation of acoustic emission source recognition for corrosion severity prediction[J].Corrosion,2011,67(5):1-11.

    [9] JOMDECHA C,PRATEEPASEN A,KAEWTRAKULPONG P.Study on source location using an acoustic emission system forvarious corrosion types[J].NDT&E International,2007,40(8):584-593.

    [10] 楊明緯,耿榮生.聲發(fā)射檢測[M].北京:機械工業(yè)出版社,2005.

    [11] 施科仁.無損檢測新技術(shù)[M].北京:清華大學(xué)出版社,2007.

    [12] GIUSEPPE C,LUIGI C,EDOARDO P.Evaluation of increasing damage severity in concrete structures by cluster analysis of acoustic emission signals:29th European Conference on Acoustic Emission Testing,September 8-10,2010[C].Vienna:Ewgae,c2010.

    [13] SAENKHUM N,PRATEEPASEN A,KEAWTRAKULPONG P.Classification of corrosion detected by acoustic emission:International Mechanical Engineering Congress,November 15-21,2003[C].Washington D C:ASME,c2003.

    [14] DIDIER L S,IDRISSI H,ROUE L.In-situ study of the cracking of metal hydride electrodes by acoustic emission technique[J].Journal of Power Source,2008,179(1):412-416.

    [15] ALIREZA F,SALVATORE S,EHSAN D N,et al. Damage assessment of reinforced concrete shear walls by acoustic emission:NDE/NDT for Highways and Bridges:Structural Materials Technology(SMT),August 21-24,2012[C].New York:ASNT,c2012.

    [16] ALIREZA F,EHSAN D N,SALVATORE S,et al.Monitoring crack propagation in reinforced concrete shear walls by acoustic emission[J].Journal of Structural Engineering,2012,139(12):1-10.

    [17] RIAHI M,SHAMEKH H.Health monitoring of aboveground storage tanks'floors:a new methodology based on practical experience[J].Journal of Nondestructive Testing(Russia),2006,42(8):437-543.

    [18] 曹楚南,張鑒清.電化學(xué)阻抗譜導(dǎo)論[M].北京:科學(xué)出版社,2002.

    [19] DING Y,REUBEN R L,STEEL J A.A new method for waveform analysis for estimating AE wave arrival times using wavelet decomposition[J].NDT&E International,2004,37(4):279-290.

    [20] HAMSTAD M A,O'GALLAGHER A,GARY A J.A wavelet transform applied to acoustic emission signals:part 1:source identification[J].Journal of Acoustic E-mission,2002,20:39-61.

    [21] JOCHEN,HARTMUT.Latest improvements on freeware AGU vallen wavelet European:29th European Conferen on Acoustic Emission Testing,September 8-10,2010[C].Vienna:Ewgae,c2010.

    [22] JIRARUNGSATIAN C,PRATEEPASEN A.Pitting and uniform corrosion source recognition using acoustic emission parameters[J].Corrosion Science,2010,52(1):187-197.

    [23] PAUL A,HWANG W J.Low-frequency resonant scattering of bubble clouds[J].Journal of Atmospheric and Oceanic Technology,2000,17:847-853.

    [24] LEIGHTON T G,PHELPS A D,RAMBLE D G,et al. Comparison of the abilities of eight acoustic techniques to detect and size a single bubble[J].Ultrasonics,1996,34(6):661-667.

    [25] JONATHAN W R B,JULIE V.The use of passive measurement of acoustic emissions from chemical engineering processes[J].Chemical Engineering Science,2001,56(5):1749-1767.

    [26] PRATEEPASEN A,JIRARUNGSATEAN C,TUENGSOOK P.Identification of AE source in corrosion process[J].Journal of Key Engineering Materials,2006,321:545-548.

    [27] 王偉魁,杜剛,曾周末,等.酸性NaCl溶液中304控氮不銹鋼腐蝕過程的聲發(fā)射特征[J].化工學(xué)報,2010,61(4):916-922. WANG Weikui,DU Gang,ZENG Zhou-mo,et al.A-coustic emission characteristics of corrosion process of 304 nitrogen controlled stainless steel in acidic NaCl solution[J].Journal of Chemical Industry and Engineering,2010,61(4):916-922.

    [28] 王偉魁,杜剛,曾周末,等.304控氮不銹鋼應(yīng)力腐蝕過程中聲發(fā)射信號聚類分析[J].化工學(xué)報,2011,62(4):1027-1033. WANG Weikui,DU Gang,ZENG Zhoumo,et al.Cluster analysis of acoustic emission signals during 304 NG stainless steel corrosion process[J].Journal of Chemical Industry and Engineering,2011,62(4):1027-1033.

    [29] RIAHI M,SHAMEKH H,KHOSROWZADEH B.Differentiation of leakage and corrosion signals in acoustic emission testing of aboveground storage tank floors with artificial neural networks[J].Journal of Nondestructive Testing(Russia),2008,44(6):436-441.

    [30] DAVIDE C,MARIO G,ALBERTO M.Development of an artificial neural network processing technique for the analysis of damage evolution in pultruded composites with acoustic emission[J].Composites Part B:Engineering,2014,56:948-959.

    [31] ANTONY J C,SUMANGALA K.Nondestructive evaluation of prestressed concrete beams using an artificial neural network(ANN)approach[J].Structural Health Monitoring,2012,5(4):313-323.

    (編輯 沈玉英)

    Cluster analysis of acoustic emission signals during tank bottom steel pitting corrosion process

    BI Haisheng1,LI Zili2,HU Dedong1,LUO Qin3,ISAAC Toku-Gyamerah2,WU Xiangyang4
    (1.College of Electromechanical Engineering in Qingdao University of Science&Technology,Qingdao 266061,China;2.College of Pipeline and Civil Engineering in China University of Petroleum,Qingdao 266580,China;3.Digitalization Engineering Department,SINOPEC Petroleum Engineering Corporation,Dongying 257000,China;4.Institute of Exploration and Development in Xingzichuan Oil Production Plant,Yanchang Oilfield,Yanan 717400,China)

    The pitting characteristics of tank bottom steel sample were studied by combined acoustic emission(AE)and electrochemical techniques in acidic NaCl solution(w=3.0%,pH=2.0).The AE signals characteristic parameters were classified using K-means clustering algorithm and each cluster signal characteristic was also extracted.The classified signals were trained using BP artificial neural network,and the AE signals from parallel experiments were successfully identified.The results show that the oscillation,movement and burst of hydrogen bubbles,breakage of passive film,growth and propagation of pit are the typical AE sources in pitting,which could be effectively classified using cluster analysis and identified by artificial neural network.It has guiding significance for interpreting and evaluating the AE on-site testing result of bottom corrosion of atmospheric storage tank,improving the reliability of testing result,reducing risk and ensuring the safety of tank.

    tank bottom steel;pitting corrosion;acoustic emission;K-means clustering;Gabor wavelet transform

    TQ 050.9

    A

    畢海勝,李自力,胡德棟,等.儲罐底板鋼點蝕過程中聲發(fā)射信號的聚類分析[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2015,39(6):145-152.

    BI Haisheng,LI Zili,HU Dedong,et al.Cluster analysis of acoustic emission signals during tank bottom steel pitting corrosion process[J].Journal of China University of Petroleum(Edition of Natural Science),2015,39(6):145-152.

    1673-5005(2015)06-0145-08

    10.3969/j.issn.1673-5005.2015.06.020

    2015-03-02

    國家自然科學(xué)基金項目(51301201);山東省自然科學(xué)基金項目(ZR2013EMQ014)

    畢海勝(1985-),男,講師,博士,研究方向為油氣儲運系統(tǒng)安全工程技術(shù)。E-mail:bhs888796@163.com。

    猜你喜歡
    信號
    信號
    鴨綠江(2021年35期)2021-04-19 12:24:18
    完形填空二則
    7個信號,警惕寶寶要感冒
    媽媽寶寶(2019年10期)2019-10-26 02:45:34
    孩子停止長個的信號
    《鐵道通信信號》訂閱單
    基于FPGA的多功能信號發(fā)生器的設(shè)計
    電子制作(2018年11期)2018-08-04 03:25:42
    基于Arduino的聯(lián)鎖信號控制接口研究
    《鐵道通信信號》訂閱單
    基于LabVIEW的力加載信號采集與PID控制
    Kisspeptin/GPR54信號通路促使性早熟形成的作用觀察
    91av网站免费观看| 欧美黄色淫秽网站| 亚洲狠狠婷婷综合久久图片| 日韩三级视频一区二区三区| 不卡av一区二区三区| 亚洲伊人色综图| 一级片免费观看大全| 久久久久久国产a免费观看| 看免费av毛片| 久久国产精品人妻蜜桃| 一卡2卡三卡四卡精品乱码亚洲| 免费久久久久久久精品成人欧美视频| 看免费av毛片| 999久久久国产精品视频| 精品国产一区二区三区四区第35| 久久中文字幕人妻熟女| 老司机午夜十八禁免费视频| 91精品国产国语对白视频| 国产免费男女视频| 国产精品九九99| 国产成人精品在线电影| 别揉我奶头~嗯~啊~动态视频| 黄色丝袜av网址大全| 国产精品亚洲美女久久久| 999久久久精品免费观看国产| 纯流量卡能插随身wifi吗| 欧美激情 高清一区二区三区| 一个人观看的视频www高清免费观看 | 90打野战视频偷拍视频| 国产成人一区二区三区免费视频网站| 精品高清国产在线一区| 此物有八面人人有两片| 免费高清在线观看日韩| 亚洲欧美激情综合另类| 亚洲国产日韩欧美精品在线观看 | 欧美黄色片欧美黄色片| 欧美老熟妇乱子伦牲交| 中文字幕精品免费在线观看视频| 88av欧美| 国产精品,欧美在线| 久久精品91蜜桃| 两性午夜刺激爽爽歪歪视频在线观看 | 桃红色精品国产亚洲av| 在线永久观看黄色视频| 一区二区日韩欧美中文字幕| av视频免费观看在线观看| 欧美+亚洲+日韩+国产| 丝袜在线中文字幕| 欧美成人午夜精品| 国产免费男女视频| 亚洲在线自拍视频| 国产精品电影一区二区三区| 国产成人欧美在线观看| 久久精品aⅴ一区二区三区四区| 夜夜看夜夜爽夜夜摸| 日韩精品中文字幕看吧| 看黄色毛片网站| 亚洲精品国产色婷婷电影| 校园春色视频在线观看| 大码成人一级视频| 久久草成人影院| 女人精品久久久久毛片| 性少妇av在线| 女警被强在线播放| 国产欧美日韩综合在线一区二区| 国产极品粉嫩免费观看在线| av福利片在线| 18禁裸乳无遮挡免费网站照片 | 女人被狂操c到高潮| 亚洲中文日韩欧美视频| 色综合欧美亚洲国产小说| 日本精品一区二区三区蜜桃| 久久青草综合色| 91老司机精品| 丰满人妻熟妇乱又伦精品不卡| 亚洲黑人精品在线| 制服丝袜大香蕉在线| 国产精品精品国产色婷婷| 在线天堂中文资源库| 19禁男女啪啪无遮挡网站| 精品人妻在线不人妻| 婷婷精品国产亚洲av在线| 亚洲人成网站在线播放欧美日韩| 九色亚洲精品在线播放| 精品久久久久久久人妻蜜臀av | 精品少妇一区二区三区视频日本电影| 亚洲中文字幕一区二区三区有码在线看 | 亚洲全国av大片| 亚洲国产毛片av蜜桃av| 久久中文看片网| 免费看美女性在线毛片视频| 亚洲精品久久国产高清桃花| 非洲黑人性xxxx精品又粗又长| 校园春色视频在线观看| 日日摸夜夜添夜夜添小说| 久久久久国产精品人妻aⅴ院| 曰老女人黄片| 亚洲专区国产一区二区| 一边摸一边做爽爽视频免费| 国产亚洲av高清不卡| 亚洲精品在线美女| 亚洲自拍偷在线| 宅男免费午夜| 国产一区二区三区在线臀色熟女| 一级,二级,三级黄色视频| 久99久视频精品免费| 亚洲avbb在线观看| 99久久99久久久精品蜜桃| 变态另类成人亚洲欧美熟女 | 热99re8久久精品国产| 国产高清videossex| 亚洲av五月六月丁香网| 国产一区在线观看成人免费| 国产一区二区激情短视频| 亚洲精华国产精华精| 在线播放国产精品三级| 久久久久久人人人人人| 国产伦人伦偷精品视频| 欧美另类亚洲清纯唯美| 999久久久精品免费观看国产| 在线天堂中文资源库| 国产精品秋霞免费鲁丝片| 91在线观看av| 亚洲激情在线av| 给我免费播放毛片高清在线观看| 色播在线永久视频| 日韩欧美一区视频在线观看| 如日韩欧美国产精品一区二区三区| 一级毛片女人18水好多| av天堂在线播放| 一夜夜www| 法律面前人人平等表现在哪些方面| 在线观看免费视频日本深夜| 欧美日韩中文字幕国产精品一区二区三区 | 成人18禁在线播放| 12—13女人毛片做爰片一| 一个人免费在线观看的高清视频| 美女扒开内裤让男人捅视频| 日本黄色视频三级网站网址| 大陆偷拍与自拍| 大码成人一级视频| 日本黄色视频三级网站网址| 黑丝袜美女国产一区| 亚洲成国产人片在线观看| 欧美另类亚洲清纯唯美| 国产亚洲精品综合一区在线观看 | 亚洲三区欧美一区| 国产成人精品在线电影| 在线播放国产精品三级| 国产精品影院久久| a级毛片在线看网站| 国产人伦9x9x在线观看| 乱人伦中国视频| 欧美日韩中文字幕国产精品一区二区三区 | 99国产精品99久久久久| av欧美777| 亚洲五月婷婷丁香| 国产伦一二天堂av在线观看| 午夜精品国产一区二区电影| 欧美在线黄色| 久久亚洲精品不卡| 国产精品九九99| 免费人成视频x8x8入口观看| 亚洲九九香蕉| 露出奶头的视频| a级毛片在线看网站| 热99re8久久精品国产| av天堂久久9| 欧美国产精品va在线观看不卡| 激情视频va一区二区三区| 国产精品久久久av美女十八| 琪琪午夜伦伦电影理论片6080| www.www免费av| 亚洲av电影不卡..在线观看| 国产成人精品久久二区二区免费| 精品国产亚洲在线| 两个人免费观看高清视频| 久久中文字幕人妻熟女| 久久九九热精品免费| 性欧美人与动物交配| 9191精品国产免费久久| 国产精品综合久久久久久久免费 | 啪啪无遮挡十八禁网站| 两性夫妻黄色片| 桃色一区二区三区在线观看| 最近最新中文字幕大全电影3 | 午夜久久久久精精品| 天天一区二区日本电影三级 | 精品一区二区三区四区五区乱码| 亚洲情色 制服丝袜| 成年版毛片免费区| 国产在线观看jvid| 99国产综合亚洲精品| 日本撒尿小便嘘嘘汇集6| 一级毛片高清免费大全| 久久久久久人人人人人| 亚洲精品美女久久久久99蜜臀| 大香蕉久久成人网| 国产av一区在线观看免费| 99精品在免费线老司机午夜| 在线观看www视频免费| 两个人免费观看高清视频| 久久久久久久精品吃奶| 老司机午夜福利在线观看视频| 午夜福利影视在线免费观看| 91成人精品电影| 黑丝袜美女国产一区| 亚洲人成伊人成综合网2020| 亚洲国产欧美网| 99国产极品粉嫩在线观看| 国产成人精品在线电影| aaaaa片日本免费| 黑人操中国人逼视频| 日本五十路高清| 国产精品精品国产色婷婷| www国产在线视频色| 久久精品亚洲精品国产色婷小说| 老汉色av国产亚洲站长工具| 国产人伦9x9x在线观看| 99re在线观看精品视频| 日韩一卡2卡3卡4卡2021年| 黄色成人免费大全| 亚洲狠狠婷婷综合久久图片| 精品国产乱子伦一区二区三区| 精品久久久久久,| 国产亚洲av嫩草精品影院| 国产免费男女视频| 国产成人精品在线电影| 一级毛片女人18水好多| 国产91精品成人一区二区三区| 女人精品久久久久毛片| 日韩免费av在线播放| 亚洲五月天丁香| 99精品在免费线老司机午夜| 亚洲成人国产一区在线观看| www.精华液| 嫩草影视91久久| 视频在线观看一区二区三区| 欧美日韩瑟瑟在线播放| 免费高清在线观看日韩| 可以在线观看毛片的网站| 精品一品国产午夜福利视频| 可以免费在线观看a视频的电影网站| 999久久久精品免费观看国产| 51午夜福利影视在线观看| 久久人妻av系列| 久久国产乱子伦精品免费另类| 国产色视频综合| 在线播放国产精品三级| 一区二区三区高清视频在线| 这个男人来自地球电影免费观看| 欧美性长视频在线观看| 99久久久亚洲精品蜜臀av| 激情视频va一区二区三区| 精品国产一区二区三区四区第35| 国产亚洲av嫩草精品影院| 欧美精品啪啪一区二区三区| 亚洲狠狠婷婷综合久久图片| 亚洲欧洲精品一区二区精品久久久| 亚洲精品一区av在线观看| 久久性视频一级片| 又紧又爽又黄一区二区| 高清黄色对白视频在线免费看| 久久 成人 亚洲| 亚洲激情在线av| 久久草成人影院| 精品电影一区二区在线| 国产欧美日韩一区二区精品| 精品国产超薄肉色丝袜足j| 满18在线观看网站| 国产成+人综合+亚洲专区| 99在线人妻在线中文字幕| 久久草成人影院| 亚洲一区二区三区色噜噜| 国产伦人伦偷精品视频| 亚洲av片天天在线观看| 88av欧美| 一卡2卡三卡四卡精品乱码亚洲| 此物有八面人人有两片| 女警被强在线播放| 午夜视频精品福利| 亚洲精品久久国产高清桃花| 国产精品免费视频内射| 最近最新中文字幕大全免费视频| 亚洲专区字幕在线| 国产精品九九99| 成人18禁高潮啪啪吃奶动态图| 国产不卡一卡二| 两人在一起打扑克的视频| 亚洲精品中文字幕在线视频| 亚洲五月婷婷丁香| 国产亚洲精品久久久久久毛片| 亚洲一码二码三码区别大吗| 欧美国产精品va在线观看不卡| videosex国产| 国产精品野战在线观看| 国产黄a三级三级三级人| 国产av一区在线观看免费| 亚洲第一av免费看| 国产单亲对白刺激| 欧美av亚洲av综合av国产av| 亚洲成a人片在线一区二区| 老熟妇仑乱视频hdxx| 久久久久九九精品影院| 国产精品亚洲一级av第二区| 免费看十八禁软件| 麻豆成人av在线观看| 性色av乱码一区二区三区2| 丝袜美足系列| 久久久久久久午夜电影| 亚洲 欧美一区二区三区| 久久精品国产清高在天天线| 欧美另类亚洲清纯唯美| 一a级毛片在线观看| 国产区一区二久久| 人人妻,人人澡人人爽秒播| 欧美色视频一区免费| 亚洲欧美激情在线| 少妇裸体淫交视频免费看高清 | 99精品久久久久人妻精品| 亚洲电影在线观看av| 电影成人av| 757午夜福利合集在线观看| 国产一区二区三区在线臀色熟女| 亚洲av片天天在线观看| 免费高清在线观看日韩| 亚洲国产日韩欧美精品在线观看 | 69精品国产乱码久久久| 老汉色∧v一级毛片| 高清黄色对白视频在线免费看| 91大片在线观看| 国产真人三级小视频在线观看| 欧美日韩一级在线毛片| www日本在线高清视频| 亚洲中文av在线| 国产99久久九九免费精品| 国产私拍福利视频在线观看| 午夜久久久久精精品| 色综合亚洲欧美另类图片| 午夜久久久在线观看| 国产真人三级小视频在线观看| 亚洲免费av在线视频| 午夜免费成人在线视频| 1024视频免费在线观看| 人人澡人人妻人| 午夜福利视频1000在线观看 | 亚洲国产高清在线一区二区三 | 国产精品乱码一区二三区的特点 | 精品国产乱码久久久久久男人| 人人妻人人澡人人看| 亚洲av片天天在线观看| 欧美在线黄色| 99香蕉大伊视频| 国产精品乱码一区二三区的特点 | 久久热在线av| 婷婷精品国产亚洲av在线| 一边摸一边抽搐一进一小说| 久久精品国产综合久久久| 欧美午夜高清在线| 亚洲欧美激情综合另类| 少妇 在线观看| 久久久久久国产a免费观看| 美女扒开内裤让男人捅视频| 老司机靠b影院| 男女下面插进去视频免费观看| 日韩大尺度精品在线看网址 | 亚洲成国产人片在线观看| 亚洲中文字幕日韩| 中文字幕另类日韩欧美亚洲嫩草| 中文字幕av电影在线播放| 日韩欧美三级三区| 可以免费在线观看a视频的电影网站| 亚洲欧美激情在线| 欧美日韩亚洲综合一区二区三区_| 91麻豆av在线| 亚洲国产欧美日韩在线播放| 18美女黄网站色大片免费观看| www日本在线高清视频| 一夜夜www| 搞女人的毛片| 国产成人精品无人区| www日本在线高清视频| 日本a在线网址| 搞女人的毛片| cao死你这个sao货| 在线观看66精品国产| www.www免费av| 亚洲中文字幕日韩| 久久久国产成人免费| 成人av一区二区三区在线看| 非洲黑人性xxxx精品又粗又长| 神马国产精品三级电影在线观看 | av网站免费在线观看视频| 国语自产精品视频在线第100页| 亚洲国产欧美一区二区综合| 成人特级黄色片久久久久久久| 青草久久国产| 侵犯人妻中文字幕一二三四区| 国产成人影院久久av| 日韩中文字幕欧美一区二区| 亚洲成av人片免费观看| 日本免费一区二区三区高清不卡 | 男人舔女人的私密视频| 日本撒尿小便嘘嘘汇集6| 欧美日本亚洲视频在线播放| 成人三级做爰电影| 日韩 欧美 亚洲 中文字幕| 美女扒开内裤让男人捅视频| 91麻豆精品激情在线观看国产| 久久影院123| 在线免费观看的www视频| 淫妇啪啪啪对白视频| 长腿黑丝高跟| 99久久99久久久精品蜜桃| 国产av精品麻豆| 无限看片的www在线观看| 午夜亚洲福利在线播放| 精品福利观看| 男人的好看免费观看在线视频 | 9热在线视频观看99| 免费人成视频x8x8入口观看| 亚洲欧美激情在线| 午夜精品久久久久久毛片777| 身体一侧抽搐| 久久久精品国产亚洲av高清涩受| 成人三级黄色视频| 日韩一卡2卡3卡4卡2021年| 999久久久精品免费观看国产| 亚洲欧美一区二区三区黑人| 50天的宝宝边吃奶边哭怎么回事| 国产熟女xx| 身体一侧抽搐| 人成视频在线观看免费观看| 亚洲情色 制服丝袜| 国产视频一区二区在线看| 亚洲精品美女久久av网站| 午夜久久久在线观看| 色在线成人网| 三级毛片av免费| 亚洲精品在线美女| 一级,二级,三级黄色视频| 久久精品亚洲熟妇少妇任你| 90打野战视频偷拍视频| 一级毛片女人18水好多| 亚洲成人国产一区在线观看| 九色国产91popny在线| 啦啦啦 在线观看视频| 亚洲午夜理论影院| 最近最新中文字幕大全免费视频| 亚洲人成网站在线播放欧美日韩| 欧美日韩一级在线毛片| 欧美av亚洲av综合av国产av| 亚洲精品国产色婷婷电影| 国产欧美日韩精品亚洲av| 精品卡一卡二卡四卡免费| 狠狠狠狠99中文字幕| 一级黄色大片毛片| 久99久视频精品免费| 亚洲美女黄片视频| 午夜成年电影在线免费观看| 久久精品影院6| 亚洲av成人一区二区三| 69av精品久久久久久| 日韩中文字幕欧美一区二区| 热99re8久久精品国产| 一级a爱片免费观看的视频| 在线播放国产精品三级| 国产一区二区在线av高清观看| av电影中文网址| 高清黄色对白视频在线免费看| 国产亚洲精品久久久久久毛片| www.www免费av| ponron亚洲| 久久久国产成人精品二区| 两个人免费观看高清视频| 日韩大尺度精品在线看网址 | 757午夜福利合集在线观看| 免费搜索国产男女视频| 12—13女人毛片做爰片一| 桃色一区二区三区在线观看| 日韩欧美一区二区三区在线观看| 国产精品野战在线观看| 高清黄色对白视频在线免费看| 亚洲国产高清在线一区二区三 | 可以在线观看的亚洲视频| 欧美av亚洲av综合av国产av| 最新美女视频免费是黄的| 欧美激情极品国产一区二区三区| 不卡av一区二区三区| a在线观看视频网站| 精品电影一区二区在线| 成熟少妇高潮喷水视频| 一二三四社区在线视频社区8| 女人爽到高潮嗷嗷叫在线视频| 夜夜爽天天搞| ponron亚洲| 久久这里只有精品19| 最好的美女福利视频网| www.熟女人妻精品国产| 美女高潮到喷水免费观看| 搞女人的毛片| 久久九九热精品免费| 此物有八面人人有两片| 久久伊人香网站| 长腿黑丝高跟| 少妇熟女aⅴ在线视频| 嫁个100分男人电影在线观看| 此物有八面人人有两片| 国产精品永久免费网站| 午夜福利免费观看在线| 亚洲国产精品合色在线| 十八禁人妻一区二区| 亚洲伊人色综图| 日本在线视频免费播放| 女生性感内裤真人,穿戴方法视频| 不卡av一区二区三区| 国产成人影院久久av| 成人av一区二区三区在线看| 亚洲一区高清亚洲精品| 午夜福利,免费看| 在线观看舔阴道视频| 视频在线观看一区二区三区| 国产精品久久久久久亚洲av鲁大| 午夜久久久在线观看| 亚洲精品国产精品久久久不卡| 欧美国产精品va在线观看不卡| 如日韩欧美国产精品一区二区三区| 色在线成人网| 91老司机精品| 亚洲欧美一区二区三区黑人| 黄色 视频免费看| 成年女人毛片免费观看观看9| 人人妻人人爽人人添夜夜欢视频| 又黄又粗又硬又大视频| 亚洲五月婷婷丁香| 天天躁狠狠躁夜夜躁狠狠躁| 久久久久久久久免费视频了| 亚洲九九香蕉| 精品一区二区三区视频在线观看免费| 夜夜夜夜夜久久久久| 欧美激情 高清一区二区三区| 好看av亚洲va欧美ⅴa在| 国产成人影院久久av| 最好的美女福利视频网| 两个人看的免费小视频| 久久中文字幕一级| 国产午夜福利久久久久久| 国产成人欧美在线观看| 亚洲精品美女久久久久99蜜臀| 天天一区二区日本电影三级 | 国产精品久久久av美女十八| 人妻丰满熟妇av一区二区三区| 51午夜福利影视在线观看| 麻豆久久精品国产亚洲av| 精品午夜福利视频在线观看一区| 亚洲熟妇中文字幕五十中出| 99国产精品一区二区三区| 麻豆国产av国片精品| 国产不卡一卡二| 麻豆一二三区av精品| 电影成人av| 久久久久亚洲av毛片大全| 久久午夜综合久久蜜桃| 日本免费a在线| 妹子高潮喷水视频| 大码成人一级视频| 久久香蕉国产精品| 黄色成人免费大全| 女人精品久久久久毛片| 亚洲久久久国产精品| www日本在线高清视频| 精品熟女少妇八av免费久了| 啦啦啦免费观看视频1| 中文字幕av电影在线播放| 中亚洲国语对白在线视频| 亚洲少妇的诱惑av| 一进一出抽搐动态| 精品国产美女av久久久久小说| 怎么达到女性高潮| 视频区欧美日本亚洲| 一二三四在线观看免费中文在| aaaaa片日本免费| 美女高潮喷水抽搐中文字幕| 国产精品乱码一区二三区的特点 | 中文字幕人成人乱码亚洲影| 国产1区2区3区精品| 在线观看午夜福利视频| 亚洲自拍偷在线| 脱女人内裤的视频| 久久亚洲真实| 免费看a级黄色片| 亚洲欧美日韩高清在线视频| 免费搜索国产男女视频| 女人高潮潮喷娇喘18禁视频| 夜夜躁狠狠躁天天躁| 国产精品久久久久久人妻精品电影| 国产精品影院久久| 我的亚洲天堂| 亚洲黑人精品在线| 亚洲精品一区av在线观看| 精品国产乱子伦一区二区三区| 久久久久久久久久久久大奶| 亚洲精品在线观看二区| 久久人人爽av亚洲精品天堂| 人成视频在线观看免费观看| 婷婷精品国产亚洲av在线| 黄色a级毛片大全视频| av网站免费在线观看视频| 免费看a级黄色片| 免费在线观看日本一区| 亚洲五月天丁香| 亚洲欧美日韩无卡精品| 脱女人内裤的视频| 午夜亚洲福利在线播放| 少妇 在线观看|