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

    基于元學(xué)習(xí)的設(shè)備故障知識(shí)圖譜構(gòu)建及推理方法

    2023-12-04 11:18:40王曉茜竇潤(rùn)亮季海鵬
    關(guān)鍵詞:故障診斷故障設(shè)備

    劉 晶,唐 震,王曉茜,竇潤(rùn)亮,季海鵬

    (1.河北工業(yè)大學(xué) 人工智能與數(shù)據(jù)科學(xué)學(xué)院,天津 300400;2.河北省數(shù)據(jù)驅(qū)動(dòng)工業(yè)智能工程研究中心,天津 300400;3.天津開發(fā)區(qū)精諾瀚海數(shù)據(jù)科技有限公司,天津 300400;4.天津大學(xué) 管理與經(jīng)濟(jì)學(xué)部,天津 300072;5.河北工業(yè)大學(xué) 材料科學(xué)與工程學(xué)院,天津 300400)

    0 引言

    工業(yè)4.0背景下,設(shè)備運(yùn)行過程越來越復(fù)雜,某些關(guān)鍵部件的故障可能導(dǎo)致整個(gè)生產(chǎn)系統(tǒng)停機(jī),因此設(shè)備故障的早期診斷和預(yù)測(cè)成為研究熱點(diǎn)[1]。近年來,數(shù)據(jù)驅(qū)動(dòng)的故障診斷方法[2]解決了因機(jī)理復(fù)雜難以建立精確數(shù)學(xué)模型的問題,在理論上大幅提高了故障診斷的準(zhǔn)確率,受到廣泛關(guān)注。在實(shí)際應(yīng)用中,由于生產(chǎn)現(xiàn)場(chǎng)的不穩(wěn)定導(dǎo)致設(shè)備運(yùn)行數(shù)據(jù)具有強(qiáng)噪聲和復(fù)雜性,運(yùn)用數(shù)據(jù)驅(qū)動(dòng)的故障診斷方法的實(shí)際故障診斷正確率遠(yuǎn)低于實(shí)驗(yàn)室結(jié)果,為了進(jìn)一步提高生產(chǎn)現(xiàn)場(chǎng)的準(zhǔn)確率,需要人機(jī)交互的智能輔助決策方法。

    2012年谷歌公司提出知識(shí)圖譜技術(shù)[3],通過構(gòu)建結(jié)構(gòu)化的知識(shí)庫(kù),推理算法求解故障成因,實(shí)現(xiàn)由知識(shí)圖譜引導(dǎo)的設(shè)備故障診斷?;谥R(shí)圖譜的知識(shí)推理算法主要包括基于距離向量的TransE模型[4]及其衍生的模型,以及基于語(yǔ)義匹配的RESCAL模型[5]、ComplEx模型[6]、DistMult模型[7]等,目前已有學(xué)者將知識(shí)推理算法應(yīng)用于設(shè)備故障診斷。喬驥等[8]將電網(wǎng)調(diào)度中大量文本規(guī)程記錄組織為知識(shí)圖網(wǎng)絡(luò),構(gòu)建了面向電網(wǎng)故障處理輔助決策的知識(shí)圖譜應(yīng)用框架,實(shí)現(xiàn)了電網(wǎng)故障處置的解析判別、輔助決策和多維人機(jī)交互。為進(jìn)一步實(shí)現(xiàn)故障數(shù)據(jù)和知識(shí)輔助之間的聯(lián)動(dòng),周彬等[9]提出一種面向設(shè)備點(diǎn)檢故障根因分析的因果知識(shí)建模方法,該方法通過提取設(shè)備故障運(yùn)維因果知識(shí)圖,結(jié)合故障知識(shí)中混雜影響因素,實(shí)現(xiàn)設(shè)備故障語(yǔ)義關(guān)系發(fā)掘,并實(shí)現(xiàn)圖譜節(jié)點(diǎn)間隱含的因果性語(yǔ)義鏈路補(bǔ)全。為了實(shí)現(xiàn)知識(shí)庫(kù)的擴(kuò)展,劉華一等[10]提出一種基于語(yǔ)義網(wǎng)的機(jī)床故障診斷知識(shí)擴(kuò)展方法和一種應(yīng)用該方法的知識(shí)擴(kuò)展平臺(tái)體系結(jié)構(gòu),該方法利用語(yǔ)義網(wǎng)技術(shù)實(shí)現(xiàn)了對(duì)機(jī)床故障知識(shí)收集效率的提升,驗(yàn)證了故障診斷知識(shí)庫(kù)在實(shí)例擴(kuò)展方面的有效性。上述方法從故障文本中提取結(jié)構(gòu)化數(shù)據(jù),構(gòu)建設(shè)備故障知識(shí)圖譜,均取得了較好的效果。但上述方法需要依賴于大量的故障樣本作為訓(xùn)練數(shù)據(jù),由于設(shè)備的高可靠性,故障數(shù)據(jù)較難獲得,如何在小樣本情況下由知識(shí)圖譜實(shí)現(xiàn)設(shè)備故障的推理非常重要。

    為實(shí)現(xiàn)小樣本情況下知識(shí)圖譜的故障推理,需要知識(shí)圖譜推理算法擁有元學(xué)習(xí)的能力。元學(xué)習(xí)[11]的目標(biāo)是使算法適應(yīng)當(dāng)前任務(wù),且能夠在新任務(wù)上利用較少的訓(xùn)練樣本取得更好的泛化表現(xiàn)?;谧罹叽硇缘脑獙W(xué)習(xí)MAML(model-agnostic meta-learning)算法[12]模型無(wú)關(guān)的思想,目前已有可在小樣本情況下實(shí)現(xiàn)知識(shí)推理的模型,包括MetaKGR(meta-based multi-hop reasoning method)[13]模型、FSRL(few-shot relation learning model)[14]模型、GANA(gated and attentive neighbor aggregator)[15]模型等,而對(duì)于非高頻三元組如何實(shí)現(xiàn)有效預(yù)測(cè),CHEN等[16]通過模型無(wú)關(guān)元學(xué)習(xí)對(duì)實(shí)體嵌入過程中的參數(shù)進(jìn)行訓(xùn)練,可以在小樣本情況下取得有價(jià)值的實(shí)體信息,實(shí)現(xiàn)了在稀疏數(shù)據(jù)集中的鏈接預(yù)測(cè)。

    由上述研究可知,目前由知識(shí)圖譜引導(dǎo)的設(shè)備故障診斷算法缺乏元學(xué)習(xí)的能力,難以在在小樣本情況下推理故障成因,且在進(jìn)行知識(shí)實(shí)體的推理時(shí),未考慮知識(shí)圖譜推理算法的負(fù)樣本生成策略,導(dǎo)致隨機(jī)生成的負(fù)樣本帶來了不必要的模型訓(xùn)練開銷。通識(shí)領(lǐng)域常采用隨機(jī)替換頭尾節(jié)點(diǎn)的方法生成負(fù)樣本,但隨機(jī)替換產(chǎn)生的負(fù)樣本對(duì)于設(shè)備故障診斷并無(wú)實(shí)際意義,反而會(huì)導(dǎo)致可能相關(guān)聯(lián)的故障間因?yàn)閷?duì)負(fù)例的學(xué)習(xí)而失去聯(lián)系,或是故障診斷模型在負(fù)例中學(xué)習(xí)不到任何有價(jià)值的信息。針對(duì)上述問題,本文提出一種基于元學(xué)習(xí)的設(shè)備故障知識(shí)圖譜構(gòu)建及推理方法,該方法通過提取設(shè)備時(shí)序信號(hào)和文本規(guī)則等多源信息構(gòu)建知識(shí)圖譜,實(shí)現(xiàn)了同一故障簇可視化及相似故障查詢,提供對(duì)故障診斷決策的輔助。并通過基于元學(xué)習(xí)的故障鏈接預(yù)測(cè)算法(Meta-Fault link prediction algorithm,MetaF),實(shí)現(xiàn)在小樣本情況下故障類型的鏈接預(yù)測(cè)和設(shè)備故障診斷的推理,為實(shí)現(xiàn)更精確的模型訓(xùn)練,采用在同一故障簇鄰域的負(fù)樣本生成策略,避免了隨機(jī)替換頭尾實(shí)體導(dǎo)致預(yù)測(cè)算法無(wú)法在負(fù)例中學(xué)習(xí)到有用的知識(shí)。最后在公開數(shù)據(jù)集上驗(yàn)證了模型故障診斷能力,證明了本文故障診斷提供輔助決策的有效性。

    1 基于元學(xué)習(xí)的設(shè)備故障知識(shí)圖譜構(gòu)建及推理方法

    基于元學(xué)習(xí)的設(shè)備故障知識(shí)圖譜構(gòu)建及推理方法框架如圖1所示。該方法包括基于設(shè)備故障信息知識(shí)圖譜構(gòu)建和基于元學(xué)習(xí)的故障鏈接預(yù)測(cè)算法MetaF兩個(gè)模塊。其中,基于設(shè)備故障信息知識(shí)圖譜構(gòu)建實(shí)現(xiàn)對(duì)設(shè)備文本數(shù)據(jù)的分類提取和設(shè)備時(shí)序數(shù)據(jù)的特征提取,得到設(shè)備故障規(guī)則鏈的描述和故障特征信息,結(jié)合故障規(guī)則和信號(hào)特征后實(shí)現(xiàn)知識(shí)圖譜的構(gòu)建?;谠獙W(xué)習(xí)的故障鏈接預(yù)測(cè)算法MetaF分別接收支持集和查詢集的任務(wù)實(shí)現(xiàn)小樣本情況下的故障診斷,支持集訓(xùn)練數(shù)據(jù)通過元學(xué)習(xí)N-way、K-shot的方式劃分得到多組任務(wù),輸入MetaF訓(xùn)練算法,查詢集則驗(yàn)證模型在新鏈接預(yù)測(cè)任務(wù)上的泛化效果。在MetaF中采用同一故障簇鄰域的負(fù)樣本生成策略進(jìn)行負(fù)樣本的生成,通過TransE算法進(jìn)行故障實(shí)體的鏈接預(yù)測(cè),最終輸出對(duì)故障類型的診斷推理結(jié)果,并提供相似故障列表查詢,為故障診斷提供輔助決策。

    1.1 設(shè)備故障信息知識(shí)圖譜的構(gòu)建

    針對(duì)傳統(tǒng)故障診斷方法缺乏故障信息結(jié)構(gòu)化管理能力的問題,提出一種構(gòu)建基于設(shè)備故障信息知識(shí)圖譜方法,并將其數(shù)據(jù)進(jìn)行劃分以適應(yīng)于元學(xué)習(xí)訓(xùn)練。構(gòu)建知識(shí)圖譜的方法如圖2所示,對(duì)于設(shè)備運(yùn)行生成的故障信息,首先從設(shè)備數(shù)據(jù)的文本描述中提取知識(shí)圖譜的故障規(guī)則鏈,其次從設(shè)備運(yùn)行數(shù)據(jù)中提取信號(hào)特征信息,最后將信號(hào)特征和故障規(guī)則鏈結(jié)合,生成統(tǒng)一描述形式的節(jié)點(diǎn)信息,完成知識(shí)圖譜的構(gòu)建。

    對(duì)于設(shè)備故障信息知識(shí)圖譜的定義如式(1)所示:

    (1)

    式中:E={e1,e2,…,en}(n<|E|)為n個(gè)實(shí)體集合,R={r1,r2,…,rm}(m<|R|)為實(shí)體間關(guān)系集合。構(gòu)建知識(shí)圖譜,需要從設(shè)備故障數(shù)據(jù)的文本數(shù)據(jù)中抽取故障規(guī)則,為不同故障類型定義分類標(biāo)準(zhǔn)。設(shè)備文本數(shù)據(jù)可以總結(jié)出多種工況特征,而設(shè)備運(yùn)行數(shù)據(jù)是在多種工況特征疊加下產(chǎn)生的結(jié)果,因此定義圖譜中的兩類節(jié)點(diǎn),分別為故障類型節(jié)點(diǎn)和故障時(shí)序片段節(jié)點(diǎn)。故障類型節(jié)點(diǎn)為某一工況下的實(shí)體表示,故障時(shí)序片段節(jié)點(diǎn)為多種工況疊加下,設(shè)備運(yùn)行時(shí)序信號(hào)采樣片段的實(shí)體表示。以某一工況下設(shè)備運(yùn)行記錄文件為例,將文本抽取為故障規(guī)則的方法如圖3所示。假設(shè)取得一個(gè)由“外圈故障”“014深度損傷”“風(fēng)扇端”這3種類型工況疊加得到的混合工況記錄文件,則分別定義“故障點(diǎn)位置”“損傷深度”“檢測(cè)位置”3種故障類型節(jié)點(diǎn)與3種工況類型對(duì)應(yīng),并定義一個(gè)故障時(shí)序片段節(jié)點(diǎn)與該混合工況記錄對(duì)應(yīng),且存在<故障時(shí)序片段節(jié)點(diǎn),belong_to,故障類型節(jié)點(diǎn)>的三元組關(guān)系。

    由此得到知識(shí)圖譜中三元組的定義,如式(2)所示。對(duì)于圖譜中的所有故障規(guī)則鏈s=(h,r,t),均有s?S,r為實(shí)體節(jié)點(diǎn)間關(guān)系描述。

    S={(h,r,t)|h,t∈E,r∈R}。

    (2)

    進(jìn)一步的,為獲得圖譜節(jié)點(diǎn)的嵌入向量形式描述,需要對(duì)設(shè)備時(shí)序數(shù)據(jù)進(jìn)行信號(hào)分解并提取特征值,將其存儲(chǔ)為知識(shí)圖譜的形式化表達(dá)結(jié)構(gòu),便于知識(shí)圖譜鏈接預(yù)測(cè)算法直接處理。對(duì)選定的時(shí)序片段序列,對(duì)其進(jìn)行變分模態(tài)分解(Variational Mode Decomposition,VMD),VMD分解通過求變分模型最優(yōu)解確定信號(hào)的各本征模態(tài)函數(shù)(Intrinsic Mode Functions,IMFs)分量,最終實(shí)現(xiàn)信號(hào)頻譜劃分和各本征模態(tài)函數(shù)分量的分離。求解模態(tài)分量方法如式(3)~式(5)所示。求解約束變分問題如公式(3)所示:

    s.t.

    (3)

    式中:{uk}為分解得到的IMFs分量,{ωk}為IMFs分量的中心頻率,K為分解的IMFs數(shù)量。采用二次懲罰因子β和拉格朗日乘數(shù)λ將約束問題轉(zhuǎn)化為非約束問題,如式(4)所示:

    L({uk},{ωk},λ)=

    (4)

    最后將上述問題通過交替方向乘子法求解,uk和ωk的最小值求解方法如式(5)所示:

    (5)

    其中:n為迭代次數(shù),當(dāng)?shù)綕M足收斂條件(6)時(shí)停止;ε為給定的判定精度,

    (6)

    根據(jù)VMD分解得到K個(gè)模態(tài)分量IMFs={IMFs1,IMFs2,…,IMFsK}后,為使知識(shí)圖譜鏈接預(yù)測(cè)算法獲取時(shí)序信號(hào)特征數(shù)據(jù),需要對(duì)IMFs分量求解其波形因子得其一維信號(hào)特征。波形因子是訊號(hào)的均方根值和整流平均值的比值,常用于時(shí)序信號(hào)的故障診斷任務(wù),其計(jì)算如式(7)所示:

    (7)

    由IMFs提取的{form1,form2,…,formk}波形因子矩陣為故障特征向量矩陣,矩陣每列為L(zhǎng)維的向量,例如formk={line1,line2,…,lineL},對(duì)應(yīng)在第k維模態(tài)分解的向量上L個(gè)知識(shí)圖譜節(jié)點(diǎn)數(shù)據(jù)。

    綜上,完成設(shè)備故障信息知識(shí)圖譜的構(gòu)建,故障規(guī)則鏈和故障特征值分別存儲(chǔ)于故障圖譜規(guī)則層和故障圖譜數(shù)據(jù)層。在圖譜的兩個(gè)層級(jí)間定義模式映射<故障時(shí)序片段節(jié)點(diǎn),故障特征節(jié)點(diǎn)>,映射表示同一故障時(shí)序片段在兩個(gè)不同層級(jí)之間的不同存儲(chǔ)方式。

    根據(jù)設(shè)備故障信息知識(shí)圖譜的數(shù)據(jù)形式,需要將其三元組劃分為元學(xué)習(xí)任務(wù)的形式。因?yàn)樵O(shè)備故障信息構(gòu)建的圖譜中,實(shí)體數(shù)量遠(yuǎn)遠(yuǎn)多于實(shí)體間關(guān)系的數(shù)量,由此設(shè)置關(guān)系為任務(wù)的劃分依據(jù)。此外,故障實(shí)體間的語(yǔ)義信息是由關(guān)系進(jìn)行連接的,由此通過關(guān)系進(jìn)行元學(xué)習(xí)任務(wù)參數(shù)訓(xùn)練是合理的。根據(jù)N-way,K-shot的原則,對(duì)元學(xué)習(xí)訓(xùn)練任務(wù)的三元組劃分方法如圖4所示。對(duì)于三元組數(shù)據(jù)集D,劃分為支持集訓(xùn)練數(shù)據(jù)Dtrain和查詢集驗(yàn)證數(shù)據(jù)Dtest,并對(duì)Dtrain進(jìn)一步分為n個(gè)任務(wù)組:Task1、Task2、…、Taskn。對(duì)于每一個(gè)任務(wù)組,其內(nèi)部均選取同一種關(guān)系的三元組,分別作為該組任務(wù)的訓(xùn)練支持集Strain和訓(xùn)練查詢集Qtrain。

    如圖4中的任務(wù)劃分方式,即使某一類的故障樣本數(shù)量較稀少,也可通過在任務(wù)中重復(fù)使用該樣本以實(shí)現(xiàn)數(shù)據(jù)增強(qiáng),解決了樣本數(shù)量少和稀缺故障難提取的問題。

    1.2 基于元學(xué)習(xí)的故障鏈接預(yù)測(cè)算法

    針對(duì)實(shí)際工況下故障數(shù)據(jù)稀缺的問題提出MetaF算法,該算法基于元學(xué)習(xí)的鏈接預(yù)測(cè)算法融入了同一故障簇鄰域的負(fù)樣本生成策略。對(duì)于知識(shí)圖譜算法而言,生成負(fù)樣本的目標(biāo)是提供高質(zhì)量的負(fù)三元組樣本,隨機(jī)生成的負(fù)樣本對(duì)于故障診斷鏈接預(yù)測(cè)算法的訓(xùn)練并無(wú)實(shí)際意義,會(huì)對(duì)最終訓(xùn)練結(jié)果產(chǎn)生影響。受知識(shí)圖譜推理算法RippleNet[17]啟發(fā),同一故障簇鄰域的負(fù)樣本生成策略模擬了水波紋的傳播方式,在知識(shí)圖譜中以需要替換的頭實(shí)體或尾實(shí)體為種子,在限定范圍的區(qū)域內(nèi)以種子為中心向外逐層查詢其他節(jié)點(diǎn)與待替換節(jié)點(diǎn)的相似度,其擴(kuò)散查詢方式如圖5所示。

    (8)

    (9)

    (10)

    (11)

    (12)

    當(dāng)替換頭實(shí)體時(shí),生成的負(fù)樣本三元組S如式(13)所示,當(dāng)替換尾實(shí)體時(shí),生成的負(fù)樣本三元組S如式(14)所示。因?yàn)樘鎿Q的節(jié)點(diǎn)是與起始節(jié)點(diǎn)相似度最小的,所以對(duì)應(yīng)的三元組應(yīng)被作為故障推理模型訓(xùn)練時(shí)的負(fù)例學(xué)習(xí),以獲得較高質(zhì)量的故障關(guān)聯(lián)信息學(xué)習(xí)。

    (13)

    (14)

    MetaF算法結(jié)構(gòu)如圖6所示。設(shè)置MetaF為兩層結(jié)構(gòu),通過鄰域負(fù)樣本生成策略構(gòu)建三元組負(fù)例,分別實(shí)現(xiàn)在支持集中關(guān)系元信息的學(xué)習(xí),以及查詢集中新實(shí)例的任務(wù)學(xué)習(xí)。MetaF在每一輪任務(wù)學(xué)習(xí)中進(jìn)行關(guān)系元的更新和故障類型的鏈接預(yù)測(cè),即通過訓(xùn)練數(shù)據(jù)更新關(guān)系元和梯度元,預(yù)測(cè)與給定缺失三元組具有關(guān)系的故障實(shí)體節(jié)點(diǎn),預(yù)測(cè)結(jié)果為故障診斷的結(jié)果。

    MetaF在小樣本情況下的故障診斷主要步驟如下:

    在支持集中,對(duì)輸入的三元組(hi,ti,ri),通過關(guān)系元學(xué)習(xí)器計(jì)算其關(guān)系元信息RTr,在關(guān)系元學(xué)習(xí)器中L層的神經(jīng)網(wǎng)絡(luò)計(jì)算對(duì)應(yīng)輸入三元組的關(guān)系元信息R(hi,ti),如式(15)~式(17)所示,其中L為神經(jīng)網(wǎng)絡(luò)層數(shù)且l∈(1,2,…,L-1),Wl為第l層權(quán)重,bl為第l層偏置項(xiàng)。

    (15)

    (16)

    (17)

    對(duì)于提取到的多個(gè)關(guān)系元信息,使用求平均的方法計(jì)算該趟任務(wù)的關(guān)系元信息RTr,如式(18)所示,其中K表示該趟任務(wù)參與訓(xùn)練的關(guān)系元數(shù)量。

    (18)

    依據(jù)知識(shí)圖譜鏈接預(yù)測(cè)算法TransE,計(jì)算當(dāng)前關(guān)系元RTr對(duì)于三元組的尾節(jié)點(diǎn)S(hi,ti),此處使用的評(píng)分函數(shù)如式(19)所示,其中‖x‖表示計(jì)算向量x的L-2范數(shù),且S(hi,ti)取值越低表示三元組計(jì)算結(jié)果越準(zhǔn)確,反之則表示該三元組為負(fù)樣本。

    (19)

    (20)

    依據(jù)損失函數(shù)L(Sr)的計(jì)算,可以計(jì)算該關(guān)系元的梯度以指導(dǎo)參數(shù)變化。梯度元GTr的計(jì)算如式(21)所示:

    (21)

    (22)

    實(shí)現(xiàn)關(guān)系元的更新后,即可在查詢集中鏈接預(yù)測(cè)未知三元組(hj,tj,rj)′的對(duì)應(yīng)故障尾節(jié)點(diǎn)S(hj,tj),且在查詢集中的得分函數(shù)和損傷分別如式(23)和式(24)所示:

    (23)

    (24)

    2 實(shí)驗(yàn)及結(jié)果分析

    2.1 鏈接預(yù)測(cè)實(shí)驗(yàn)

    為驗(yàn)證MetaF在知識(shí)圖譜網(wǎng)絡(luò)中的鏈接預(yù)測(cè)能力,以及小樣本情況下算法的有效性,選擇公開數(shù)據(jù)集NELL-One作為實(shí)驗(yàn)的基準(zhǔn)數(shù)據(jù)集。NELL數(shù)據(jù)集是在互聯(lián)網(wǎng)信息中自動(dòng)抽取新知識(shí)的知識(shí)庫(kù),NELL-One是NELL知識(shí)圖譜數(shù)據(jù)集的子集,其根據(jù)關(guān)系出現(xiàn)的次數(shù),選擇小于500次且大于50次的三元組,提取并重新組合為數(shù)據(jù)集,常用于驗(yàn)證算法在小樣本情況下知識(shí)圖譜的鏈接預(yù)測(cè)能力。

    為評(píng)估模型效果,根據(jù)關(guān)系出現(xiàn)的頻率確定三元組出現(xiàn)的頻率,定義出現(xiàn)頻率較低的關(guān)系為目標(biāo)關(guān)系,用于評(píng)估小樣本情況下的效果,非低頻關(guān)系為背景關(guān)系。采用Pre-Train和In-Train兩種模型訓(xùn)練方法分別評(píng)估其效果,Pre-Train在模型訓(xùn)練前結(jié)合背景關(guān)系進(jìn)行模型預(yù)訓(xùn)練,In-Train為直接使用背景關(guān)系和目標(biāo)關(guān)系直接訓(xùn)練模型。選擇MRR、Hits@10、Hits@5、Hits@1作為實(shí)驗(yàn)評(píng)估指標(biāo),MRR為測(cè)試集數(shù)據(jù)上的平均排名,Hits@10、Hits@5、Hits@1分別為正確結(jié)果在預(yù)測(cè)結(jié)果中排名前10、前5、前1的比例。對(duì)比試驗(yàn)選擇的模型為FSRL[14]、GANA[15]、MetaR[16]算法,分別在1-shot和5-shot兩種情況下進(jìn)行實(shí)驗(yàn),實(shí)驗(yàn)結(jié)果如表1和表2所示。

    表1 1-shot實(shí)驗(yàn)結(jié)果 %

    表2 5-shot實(shí)驗(yàn)結(jié)果 %

    根據(jù)實(shí)驗(yàn)結(jié)果可看出,由MetaR優(yōu)化得到的MetaF在1-shot和5-shot兩種情況下,于In-Train的模型訓(xùn)練中均取得了優(yōu)于原始模型的效果。對(duì)于FSRL模型,MetaF在MRR上獲得了3.7%(1-shot)和2.5%(5-shot)的提升,在Hits@10上獲得了2.1%(1-shot)和3.2%(5-shot)的提升;而對(duì)于GANA模型,MetaF在MRR上獲得了2.5%(1-shot)和0.8%(5-shot)的提升,在Hits@10上獲得了1.9%(1-shot)和0.9%(5-shot)的提升,說明MetaF在鏈接預(yù)測(cè)的性能上與當(dāng)前主流方法基本一致。在Pre-Train訓(xùn)練中,MetaF與MetaR相比并未取得明顯的進(jìn)步,由于Pre-Trian過程會(huì)對(duì)實(shí)體嵌入進(jìn)行預(yù)訓(xùn)練,且根據(jù)負(fù)樣本生成策略增加了數(shù)據(jù)集中僅出現(xiàn)一次的三元組數(shù)量,因此影響了MetaF在知識(shí)圖譜中鏈接預(yù)測(cè)的表現(xiàn)。對(duì)于In-Train訓(xùn)練,MetaF較MetaR在MRR上取得了5.2%(1-shot)和2.2%(5-shot)的進(jìn)步,而對(duì)于發(fā)現(xiàn)相似故障列表最有評(píng)估價(jià)值的Hits@10,MetaF較MetaR取得了7.7%(1-shot)和2.9%(5-shot)的進(jìn)步。說明MetaF可有效在包含背景知識(shí)的訓(xùn)練過程In-Trian中進(jìn)行鏈接預(yù)測(cè)任務(wù)。

    2.2 故障診斷實(shí)驗(yàn)

    故障診斷實(shí)驗(yàn)分為設(shè)備故障信息知識(shí)圖譜的構(gòu)建和MetaF故障推理兩個(gè)部分,首先抽取故障數(shù)據(jù)集特征信息并構(gòu)建知識(shí)圖譜,將圖譜以三元組的形式存儲(chǔ)于txt文本文檔和csv文件中,并將圖譜導(dǎo)入Neo4j數(shù)據(jù)庫(kù),MetaF模型使用PyTorch包構(gòu)建模型進(jìn)行訓(xùn)練。

    2.2.1 數(shù)據(jù)集介紹

    實(shí)驗(yàn)選取設(shè)備故障診斷領(lǐng)域常用的美國(guó)凱斯西儲(chǔ)大學(xué)(Case Western Reserve University,CWRU)滾動(dòng)軸承故障數(shù)據(jù)集,已有多篇文獻(xiàn)使用該數(shù)據(jù)集驗(yàn)證算法[18-19]。CWRU數(shù)據(jù)集包含由加速度傳感器采集的軸承振動(dòng)信號(hào)數(shù)據(jù)和對(duì)應(yīng)的機(jī)理說明文檔。CWRU數(shù)據(jù)包括在軸承試驗(yàn)臺(tái)驅(qū)動(dòng)端(Drive)、風(fēng)扇端(Fan)和基座端(Base)采集的振動(dòng)信號(hào)。采用電火花加工單點(diǎn)損傷,損傷直徑有0.007英寸、0.014英寸、0.021英寸3種,損傷部位有內(nèi)圈、外圈和滾動(dòng)體3種,其中外圈損傷點(diǎn)分別在3點(diǎn)鐘、6點(diǎn)鐘、12點(diǎn)鐘的方位設(shè)置。

    2.2.2 設(shè)備故障信息知識(shí)圖譜的構(gòu)建

    實(shí)驗(yàn)選取了CWRU數(shù)據(jù)集中所有12 kHz數(shù)據(jù),根據(jù)數(shù)據(jù)集的文本規(guī)則描述,對(duì)數(shù)據(jù)進(jìn)行分類提取,總結(jié)出5種分類可以為數(shù)據(jù)進(jìn)行劃分,具體為:端點(diǎn)、故障點(diǎn)位置、損傷深度、載荷和檢測(cè)位置,如表3所示。

    表3 CWRU數(shù)據(jù)集分類規(guī)則

    根據(jù)以上分類規(guī)則,可得圖譜中節(jié)點(diǎn)間的對(duì)應(yīng)關(guān)系,如表4所示。根據(jù)每一個(gè)故障采樣片段在5個(gè)維度上的屬性值,得到故障時(shí)序片段節(jié)點(diǎn)和故障類型節(jié)點(diǎn)間的連接關(guān)系。

    表4 節(jié)點(diǎn)間對(duì)應(yīng)關(guān)系

    根據(jù)設(shè)備故障信息知識(shí)圖譜的節(jié)點(diǎn)定義和數(shù)據(jù)關(guān)系,對(duì)軸承故障信號(hào)劃分片段并進(jìn)行變分模態(tài)分解,生成的特征向量與節(jié)點(diǎn)結(jié)合構(gòu)建知識(shí)圖譜,提取到1 196個(gè)實(shí)體和6種關(guān)系,構(gòu)建得到的知識(shí)圖譜局部結(jié)構(gòu)如圖7所示,其描述了CWRU軸承故障數(shù)據(jù)在端點(diǎn)、故障點(diǎn)位置、損傷深度、載荷和檢測(cè)位置這5種維度下的關(guān)聯(lián)方式,并實(shí)現(xiàn)了不同類型故障下的故障片段間的知識(shí)鏈接。知識(shí)圖譜中任意兩個(gè)故障時(shí)序片段節(jié)點(diǎn)間的距離均為兩跳,中間由與其相關(guān)聯(lián)的故障類型節(jié)點(diǎn)連接,這樣的結(jié)構(gòu)便于在某一區(qū)域內(nèi)鏈接預(yù)測(cè)故障片段真正的故障類型。

    2.2.3 故障推理實(shí)驗(yàn)和分析

    為驗(yàn)證在小樣本條件下,MetaF可以較好地完成故障推理任務(wù),選取4種不同工況下的設(shè)備正常運(yùn)行數(shù)據(jù)為基準(zhǔn),并以9:1的比例將正常數(shù)據(jù)與故障數(shù)據(jù)混合,構(gòu)建小樣本故障數(shù)據(jù)集,數(shù)據(jù)集的劃分如表5所示。

    表5 小樣本故障數(shù)據(jù)集劃分方法

    分別測(cè)試4組小樣本數(shù)據(jù)集在MetaF數(shù)據(jù)集上的表現(xiàn)效果,以平均排名MRR和精確率Precision為評(píng)價(jià)指標(biāo),驗(yàn)證算法推理故障實(shí)體的能力。MetaF由鄰域負(fù)樣本生成模塊和故障鏈接預(yù)測(cè)模塊組成,因此設(shè)置移除鄰域負(fù)樣本生成模塊作為消融實(shí)驗(yàn)對(duì)照組,記為模型M_min,該組實(shí)驗(yàn)僅使用隨機(jī)生成的負(fù)樣本,以驗(yàn)證負(fù)樣本生成策略的有效性,故障診斷實(shí)驗(yàn)和消融實(shí)驗(yàn)如表6所示。

    表6 故障診斷實(shí)驗(yàn)和消融實(shí)驗(yàn) %

    由表6的故障診斷實(shí)驗(yàn)結(jié)果可以看出,MetaF在4組故障樣本稀疏的數(shù)據(jù)集上均表現(xiàn)出較為均衡的效果,并未因?yàn)楣收项愋偷牟煌瑢?dǎo)致模型性能的變化。對(duì)于表6中作為對(duì)照的消融實(shí)驗(yàn)結(jié)果,MetaF比作為對(duì)照的M_min,在MRR和Precision上均有提升,對(duì)于4個(gè)數(shù)據(jù)集,MetaF在MRR上的提升分別為6.8%、6.4%、7.3%和7.2%,并且MetaF在Precision上的提升分別為14.1%、14.0%、13.6%和13.2%,說明鄰域負(fù)樣本生成模塊的加入有效提高了MetaF推理故障實(shí)體的能力,同時(shí)表明MetaF在小樣本數(shù)據(jù)的故障診斷中表現(xiàn)良好。

    由于傳統(tǒng)基于數(shù)據(jù)驅(qū)動(dòng)的故障診斷方法無(wú)法結(jié)構(gòu)化管理故障特征信息,在故障數(shù)據(jù)內(nèi)容發(fā)生變化時(shí),需要重新訓(xùn)練模型。而基于知識(shí)圖譜的方法因?yàn)橐延袑?duì)故障信息的描述,在理論上具有一定程度的優(yōu)勢(shì),選取DS_01小樣本數(shù)據(jù)集為實(shí)驗(yàn)數(shù)據(jù),驗(yàn)證知識(shí)圖譜方法對(duì)比數(shù)據(jù)驅(qū)動(dòng)方法,在模型收斂速度和運(yùn)行時(shí)間上的優(yōu)勢(shì)。選取故障診斷領(lǐng)域常用的LSTM算法[20]與MetaF算法對(duì)比,模型的損失值曲線對(duì)比如圖8所示,平均運(yùn)行時(shí)間對(duì)比如圖9所示。

    根據(jù)實(shí)驗(yàn)結(jié)果顯示,MetaF獲得了比數(shù)據(jù)驅(qū)動(dòng)的LSTM算法更快的模型收斂速度和更短的平均運(yùn)行時(shí)間,說明基于知識(shí)圖譜的方法在模型的運(yùn)行速度和運(yùn)行時(shí)間上優(yōu)于基于數(shù)據(jù)的方法,MetaF依賴結(jié)構(gòu)化的故障信息存儲(chǔ),有效利用了已知信息加快故障樣本向量計(jì)算的速度,實(shí)現(xiàn)故障診斷的效率提升。

    2.2.4 相似故障查詢結(jié)果和分析

    MetaF實(shí)現(xiàn)相似故障查詢的流程如圖10所示,在輸入故障采樣片段后,MetaF通過知識(shí)圖譜中已有故障信息確定樣本故障類型并更新模型,得到故障類型后更新知識(shí)圖譜,最后在更新的節(jié)點(diǎn)鄰域內(nèi)進(jìn)行可視化的相似故障查詢。

    實(shí)驗(yàn)隨機(jī)選取CWRU數(shù)據(jù)驗(yàn)證集中的一段采樣片段作為輸入,采樣長(zhǎng)度為1 600,將該采樣片段進(jìn)行模態(tài)分解后輸入MetaF模型進(jìn)行故障類型的診斷,并根據(jù)診斷結(jié)果進(jìn)行相似故障查詢,查詢結(jié)果如圖11所示。根據(jù)輸入的故障片段,可查詢出一組目標(biāo)相似故障簇,以及3組其他相似故障簇。

    對(duì)于查詢得到的相似故障列表,根據(jù)故障向量間的歐式距離進(jìn)行故障相似度的評(píng)分,對(duì)于任一組故障向量P={p1,p2,…,pk}和Q={q1,q2,…,qk},其故障相似度Fault_Similarity計(jì)算方法如式(25)所示,式中k表示向量的維數(shù)。

    (25)

    根據(jù)輸入的故障片段,分別與求得的相似故障列表中每一項(xiàng)計(jì)算其故障相似度,得到評(píng)分結(jié)果統(tǒng)計(jì)如圖12所示。

    根據(jù)評(píng)分結(jié)果,可得片段141130020、片段141230070、片段141230040、片段141230010這4組數(shù)據(jù)與輸入故障最相似,進(jìn)一步查表可得這4組數(shù)據(jù)與輸入故障均為故障類型faultType = {point:Drive,located:B,deepth:007,load:1}。在推理的14個(gè)相似故障中發(fā)現(xiàn)4個(gè)故障樣本與輸入故障相似,則稱目標(biāo)相似故障簇的樣本檢出率為28.6%。通過該組相似故障查詢實(shí)驗(yàn),可說明本文提出的方法在故障診斷應(yīng)用上的有效性。為進(jìn)一步說明MetaF的故障推理能力較于MetaR有提升,分別對(duì)比兩種模型在相似故障查詢中目標(biāo)相似故障簇的樣本檢出率。設(shè)置實(shí)驗(yàn)分別在相似故障列表長(zhǎng)度為10、15、20的情況下進(jìn)行,每組實(shí)驗(yàn)進(jìn)行10次取其結(jié)果的平均值,實(shí)驗(yàn)結(jié)果如表7所示。

    表7 目標(biāo)相似故障簇的樣本檢出率 %

    由表7的實(shí)驗(yàn)結(jié)果可看出,MetaF在相似故障列表長(zhǎng)度分別為10、15、20時(shí),MetaF相較于MetaR提升了11.1%、6.4%、8.5%,說明同一故障簇鄰域負(fù)樣本生成策略的加入提升了模型故障診斷的能力,并通過相似故障查詢有效輔助了故障決策。

    2.2.5 實(shí)際工程應(yīng)用案例分析

    以電力變壓器故障數(shù)據(jù)診斷為例,驗(yàn)證MetaF發(fā)現(xiàn)設(shè)備故障和相似故障查詢的有效性。電力變壓器運(yùn)行數(shù)據(jù)包括正常運(yùn)行以及低能放電、低溫過熱、高能放電、高溫過熱、中低溫過熱和中溫過熱這6種故障。依據(jù)5種變壓器綜合故障特征,定義變壓器數(shù)據(jù)標(biāo)簽為fault_tag={feature01,feature02,feature03,feature04,feature05},5種特征分別對(duì)應(yīng)變壓器五種油中溶解氣體含量計(jì)算得到的特征值,將已有數(shù)據(jù)劃分為DGA_01={正常運(yùn)行,低能放電,低溫過熱},DGA_02={正常運(yùn)行,高能放電,高溫過熱},DGA_03={正常運(yùn)行,中低溫過熱、中溫過熱}3組實(shí)驗(yàn)數(shù)據(jù),將正常數(shù)據(jù)和故障數(shù)據(jù)以9:1的比例混合,分別在3個(gè)數(shù)據(jù)集上進(jìn)行知識(shí)圖譜構(gòu)建,并在圖譜基礎(chǔ)上完成故障診斷任務(wù)。變壓器故障圖譜的構(gòu)建及其故障診斷結(jié)果分別如圖13和圖14所示。

    圖14分別是MetaF在3組變壓器故障數(shù)據(jù)的精確率和相似故障的樣本檢出率,可看出相較于CWRU數(shù)據(jù)集,變壓器數(shù)據(jù)獲得了更高的相似故障樣本檢出率,這是由于變壓器數(shù)據(jù)構(gòu)建的知識(shí)圖譜結(jié)構(gòu)相對(duì)簡(jiǎn)單,因此獲得了更高的相似故障命中。根據(jù)實(shí)驗(yàn)驗(yàn)證了MetaF在完成較高準(zhǔn)確度的故障診斷時(shí),還可提供較好的相似故障查詢,有效輔助了設(shè)備故障診斷決策。

    3 結(jié)束語(yǔ)

    針對(duì)小樣本情況下難以構(gòu)建故障數(shù)據(jù)知識(shí)圖譜,從而導(dǎo)致推理算法無(wú)法準(zhǔn)確實(shí)現(xiàn)故障診斷的問題,本文提出一種基于元學(xué)習(xí)的設(shè)備故障知識(shí)圖譜構(gòu)建及推理方法。該方法包括兩部分:第一部分為根據(jù)故障數(shù)據(jù)文本描述提取數(shù)據(jù)的分類規(guī)則,構(gòu)建本體模型的故障規(guī)則鏈,并提取故障數(shù)據(jù)中的特征信息作為故障規(guī)則鏈的數(shù)據(jù)信息,實(shí)現(xiàn)設(shè)備故障信息知識(shí)圖譜的構(gòu)建;第二部分利用同一故障簇鄰域范圍內(nèi)的負(fù)樣本生成策略,通過MetaF模型進(jìn)行小樣本情況下的故障診斷和相似故障列表查詢,實(shí)現(xiàn)對(duì)設(shè)備運(yùn)行數(shù)據(jù)的故障診斷決策輔助。實(shí)驗(yàn)結(jié)果表明,MetaF有效提高了在知識(shí)圖譜上推理的準(zhǔn)確率,并能有效利用設(shè)備運(yùn)行數(shù)據(jù)構(gòu)建故障診斷領(lǐng)域知識(shí)圖譜,提供相似故障列表和可視化的相似故障簇,為故障診斷提供了有效的知識(shí)支撐。

    猜你喜歡
    故障診斷故障設(shè)備
    諧響應(yīng)分析在設(shè)備減振中的應(yīng)用
    故障一點(diǎn)通
    基于MPU6050簡(jiǎn)單控制設(shè)備
    電子制作(2018年11期)2018-08-04 03:26:08
    奔馳R320車ABS、ESP故障燈異常點(diǎn)亮
    500kV輸變電設(shè)備運(yùn)行維護(hù)探討
    因果圖定性分析法及其在故障診斷中的應(yīng)用
    故障一點(diǎn)通
    江淮車故障3例
    原來他們都是可穿戴設(shè)備
    基于LCD和排列熵的滾動(dòng)軸承故障診斷
    av免费观看日本| 22中文网久久字幕| 99久久人妻综合| 最近中文字幕高清免费大全6| 国产精品人妻久久久影院| 欧美日韩精品成人综合77777| 免费高清在线观看日韩| 99国产精品免费福利视频| 久久精品aⅴ一区二区三区四区 | 国产极品粉嫩免费观看在线| 国产av精品麻豆| 久久久欧美国产精品| 国产色爽女视频免费观看| 久久 成人 亚洲| 亚洲,一卡二卡三卡| 搡老乐熟女国产| 亚洲一级一片aⅴ在线观看| 久久久久精品久久久久真实原创| 91午夜精品亚洲一区二区三区| 亚洲第一区二区三区不卡| 日韩熟女老妇一区二区性免费视频| 欧美日本中文国产一区发布| 久久毛片免费看一区二区三区| 女的被弄到高潮叫床怎么办| 99热这里只有是精品在线观看| a级毛色黄片| 亚洲av免费高清在线观看| 我的女老师完整版在线观看| 街头女战士在线观看网站| tube8黄色片| 日韩三级伦理在线观看| 国产成人91sexporn| 街头女战士在线观看网站| 亚洲av国产av综合av卡| 国产有黄有色有爽视频| 高清视频免费观看一区二区| 青春草视频在线免费观看| 男女午夜视频在线观看 | 不卡视频在线观看欧美| 国产淫语在线视频| 精品人妻一区二区三区麻豆| 九色亚洲精品在线播放| 老司机影院成人| 99久国产av精品国产电影| 大片免费播放器 马上看| 日本av手机在线免费观看| 日韩三级伦理在线观看| 又大又黄又爽视频免费| 中文字幕人妻熟女乱码| 国产熟女午夜一区二区三区| 香蕉国产在线看| 国产极品天堂在线| 久久综合国产亚洲精品| 精品人妻熟女毛片av久久网站| a级片在线免费高清观看视频| 高清在线视频一区二区三区| 午夜福利视频精品| 只有这里有精品99| 国产免费现黄频在线看| 亚洲少妇的诱惑av| 激情视频va一区二区三区| 日韩成人av中文字幕在线观看| 成人国语在线视频| 欧美亚洲 丝袜 人妻 在线| 国产欧美亚洲国产| 三上悠亚av全集在线观看| 一边摸一边做爽爽视频免费| 婷婷色综合www| 精品少妇内射三级| 国产男人的电影天堂91| 国产精品国产三级专区第一集| 日韩欧美一区视频在线观看| 久久99热6这里只有精品| 丰满少妇做爰视频| 一级a做视频免费观看| 国产黄色免费在线视频| 宅男免费午夜| 精品一区二区三卡| 女的被弄到高潮叫床怎么办| 免费看av在线观看网站| 如何舔出高潮| 又黄又爽又刺激的免费视频.| 不卡视频在线观看欧美| 一级毛片黄色毛片免费观看视频| 丝袜人妻中文字幕| 免费在线观看黄色视频的| 国产激情久久老熟女| 久久ye,这里只有精品| 一二三四中文在线观看免费高清| 欧美人与性动交α欧美软件 | 精品一区二区三卡| 亚洲国产日韩一区二区| 久久99一区二区三区| 婷婷色综合大香蕉| 久久人人爽人人爽人人片va| 美女xxoo啪啪120秒动态图| 熟女人妻精品中文字幕| 午夜视频国产福利| 丝袜人妻中文字幕| 黑人欧美特级aaaaaa片| 亚洲精品中文字幕在线视频| 亚洲性久久影院| 成人手机av| 免费大片18禁| 黄色 视频免费看| av有码第一页| 日韩伦理黄色片| 91精品国产国语对白视频| 在线亚洲精品国产二区图片欧美| 国产一区二区激情短视频 | 综合色丁香网| 纯流量卡能插随身wifi吗| 岛国毛片在线播放| 亚洲av成人精品一二三区| 亚洲国产最新在线播放| 十八禁高潮呻吟视频| 国产av精品麻豆| 亚洲精品色激情综合| 热99国产精品久久久久久7| 国产xxxxx性猛交| 纵有疾风起免费观看全集完整版| 夫妻午夜视频| 久久久久国产精品人妻一区二区| 97超碰精品成人国产| 国产精品国产三级专区第一集| 亚洲国产精品专区欧美| 色吧在线观看| 2018国产大陆天天弄谢| 精品少妇久久久久久888优播| 色视频在线一区二区三区| 午夜福利影视在线免费观看| 三级国产精品片| 老熟女久久久| 久久这里有精品视频免费| 啦啦啦中文免费视频观看日本| 我要看黄色一级片免费的| 亚洲av综合色区一区| 十八禁高潮呻吟视频| 人妻 亚洲 视频| 不卡视频在线观看欧美| 日本免费在线观看一区| 日韩一本色道免费dvd| 亚洲天堂av无毛| 七月丁香在线播放| 日本猛色少妇xxxxx猛交久久| 最近中文字幕高清免费大全6| 国产亚洲精品久久久com| 最近2019中文字幕mv第一页| 日韩一区二区三区影片| 国产高清国产精品国产三级| 亚洲 欧美一区二区三区| 国产精品久久久久久精品电影小说| 精品人妻一区二区三区麻豆| 国产成人精品福利久久| 香蕉国产在线看| 免费观看无遮挡的男女| 日日啪夜夜爽| 狠狠婷婷综合久久久久久88av| xxx大片免费视频| 下体分泌物呈黄色| 亚洲欧美成人精品一区二区| 街头女战士在线观看网站| 日日爽夜夜爽网站| 亚洲伊人色综图| 国产熟女欧美一区二区| 啦啦啦在线观看免费高清www| 美女国产高潮福利片在线看| 两个人免费观看高清视频| 日产精品乱码卡一卡2卡三| 国产精品女同一区二区软件| 晚上一个人看的免费电影| 黄色配什么色好看| 桃花免费在线播放| 欧美国产精品va在线观看不卡| 亚洲国产看品久久| 男女啪啪激烈高潮av片| 亚洲情色 制服丝袜| 国产精品久久久久久久久免| 亚洲成av片中文字幕在线观看 | 黄色毛片三级朝国网站| 亚洲四区av| 天天躁夜夜躁狠狠躁躁| 亚洲图色成人| 高清视频免费观看一区二区| 久久女婷五月综合色啪小说| 18禁动态无遮挡网站| 十八禁高潮呻吟视频| 80岁老熟妇乱子伦牲交| 婷婷色麻豆天堂久久| 黑人欧美特级aaaaaa片| 99热全是精品| 欧美激情国产日韩精品一区| 哪个播放器可以免费观看大片| 天堂8中文在线网| 久久人人爽人人片av| 免费看不卡的av| 精品国产露脸久久av麻豆| 精品久久久久久电影网| 欧美人与性动交α欧美软件 | 母亲3免费完整高清在线观看 | 日本黄大片高清| 色婷婷av一区二区三区视频| 国产1区2区3区精品| 一边摸一边做爽爽视频免费| 在线观看免费视频网站a站| 美女脱内裤让男人舔精品视频| 69精品国产乱码久久久| 男男h啪啪无遮挡| 新久久久久国产一级毛片| 99视频精品全部免费 在线| 一级毛片 在线播放| √禁漫天堂资源中文www| 亚洲情色 制服丝袜| 国产免费一区二区三区四区乱码| 国产精品久久久久久精品电影小说| 麻豆精品久久久久久蜜桃| 久久精品国产亚洲av天美| 少妇的逼水好多| 日韩成人伦理影院| 90打野战视频偷拍视频| 边亲边吃奶的免费视频| 国产欧美日韩综合在线一区二区| 爱豆传媒免费全集在线观看| 欧美性感艳星| 最近手机中文字幕大全| 久久久久国产精品人妻一区二区| 在线观看www视频免费| videos熟女内射| 免费人成在线观看视频色| 久久亚洲国产成人精品v| 亚洲精品日韩在线中文字幕| 人人妻人人添人人爽欧美一区卜| 亚洲国产欧美日韩在线播放| 91精品三级在线观看| 精品国产乱码久久久久久小说| 啦啦啦中文免费视频观看日本| 一区二区三区乱码不卡18| 侵犯人妻中文字幕一二三四区| 青春草视频在线免费观看| 日本与韩国留学比较| 久久精品久久精品一区二区三区| 妹子高潮喷水视频| videosex国产| 中文字幕制服av| 大陆偷拍与自拍| 少妇人妻久久综合中文| av女优亚洲男人天堂| 国产成人午夜福利电影在线观看| 如日韩欧美国产精品一区二区三区| 91精品国产国语对白视频| 十八禁高潮呻吟视频| 午夜福利网站1000一区二区三区| 欧美国产精品一级二级三级| 妹子高潮喷水视频| 国产av精品麻豆| 欧美少妇被猛烈插入视频| 天堂俺去俺来也www色官网| 免费看光身美女| 日韩精品有码人妻一区| 成人无遮挡网站| 亚洲丝袜综合中文字幕| 国产成人精品无人区| 狂野欧美激情性xxxx在线观看| 国产精品久久久久久久久免| 免费高清在线观看日韩| 成人国产麻豆网| 国产免费视频播放在线视频| 曰老女人黄片| 久久综合国产亚洲精品| 麻豆乱淫一区二区| 18禁国产床啪视频网站| 日本av免费视频播放| 人人澡人人妻人| 婷婷色av中文字幕| 看免费成人av毛片| 在线观看三级黄色| 成年人午夜在线观看视频| 人妻 亚洲 视频| 日韩伦理黄色片| a级毛片黄视频| 久久久久久久久久成人| av免费观看日本| 啦啦啦视频在线资源免费观看| 国产精品久久久久久久久免| 国产成人精品在线电影| 欧美日韩视频精品一区| 一二三四中文在线观看免费高清| 精品亚洲成a人片在线观看| 永久免费av网站大全| 国产激情久久老熟女| 大香蕉97超碰在线| 欧美精品高潮呻吟av久久| 寂寞人妻少妇视频99o| 多毛熟女@视频| 18禁动态无遮挡网站| 美女脱内裤让男人舔精品视频| 国产欧美另类精品又又久久亚洲欧美| 国产片特级美女逼逼视频| 亚洲丝袜综合中文字幕| 天天操日日干夜夜撸| 国产激情久久老熟女| 欧美xxxx性猛交bbbb| 色5月婷婷丁香| 一级毛片黄色毛片免费观看视频| 看非洲黑人一级黄片| 日日爽夜夜爽网站| 女人被躁到高潮嗷嗷叫费观| 亚洲性久久影院| 久久99精品国语久久久| 91在线精品国自产拍蜜月| tube8黄色片| 汤姆久久久久久久影院中文字幕| 亚洲,一卡二卡三卡| 一个人免费看片子| 亚洲精品日韩在线中文字幕| 国精品久久久久久国模美| 九九爱精品视频在线观看| 国产成人午夜福利电影在线观看| 肉色欧美久久久久久久蜜桃| 日本wwww免费看| videossex国产| 亚洲国产毛片av蜜桃av| 狠狠婷婷综合久久久久久88av| 观看av在线不卡| 两个人看的免费小视频| 日韩 亚洲 欧美在线| 久久久久久久亚洲中文字幕| 免费av不卡在线播放| 两性夫妻黄色片 | 成人毛片a级毛片在线播放| 99久久中文字幕三级久久日本| 国产精品久久久久久av不卡| 亚洲欧美日韩另类电影网站| 精品亚洲成a人片在线观看| 少妇高潮的动态图| 桃花免费在线播放| 婷婷色综合www| 人妻 亚洲 视频| 国产熟女午夜一区二区三区| 国产精品一区二区在线不卡| 亚洲国产最新在线播放| 日日啪夜夜爽| 久久精品国产亚洲av涩爱| 亚洲精品日韩在线中文字幕| 中文精品一卡2卡3卡4更新| 美女内射精品一级片tv| 精品99又大又爽又粗少妇毛片| 18禁国产床啪视频网站| 日韩av不卡免费在线播放| 免费久久久久久久精品成人欧美视频 | 免费在线观看完整版高清| 肉色欧美久久久久久久蜜桃| 亚洲欧美成人综合另类久久久| 欧美日韩视频高清一区二区三区二| 又黄又粗又硬又大视频| 99久久人妻综合| 欧美xxⅹ黑人| 免费大片18禁| 欧美 日韩 精品 国产| videosex国产| 成人国语在线视频| 人妻少妇偷人精品九色| 少妇被粗大的猛进出69影院 | 中国国产av一级| 国产乱人偷精品视频| 22中文网久久字幕| 久久99热这里只频精品6学生| 国产精品99久久99久久久不卡 | 午夜免费鲁丝| 久久精品国产鲁丝片午夜精品| 国产在线视频一区二区| 日韩熟女老妇一区二区性免费视频| 中文天堂在线官网| 99热国产这里只有精品6| 黑人猛操日本美女一级片| 中文乱码字字幕精品一区二区三区| 涩涩av久久男人的天堂| 欧美最新免费一区二区三区| 中国三级夫妇交换| 一区二区三区四区激情视频| 女人被躁到高潮嗷嗷叫费观| 久久99一区二区三区| 成人影院久久| 建设人人有责人人尽责人人享有的| 国产精品国产av在线观看| 你懂的网址亚洲精品在线观看| 啦啦啦在线观看免费高清www| 亚洲五月色婷婷综合| 99国产精品免费福利视频| 人人妻人人澡人人看| 亚洲综合精品二区| 国产又爽黄色视频| 亚洲一码二码三码区别大吗| 少妇的逼好多水| 好男人视频免费观看在线| tube8黄色片| 亚洲精品日韩在线中文字幕| av国产精品久久久久影院| 如何舔出高潮| 90打野战视频偷拍视频| 大码成人一级视频| 王馨瑶露胸无遮挡在线观看| 国产精品国产三级国产专区5o| 久久人人97超碰香蕉20202| 国产深夜福利视频在线观看| 涩涩av久久男人的天堂| 国产精品久久久av美女十八| 成年动漫av网址| 日韩,欧美,国产一区二区三区| 精品熟女少妇av免费看| 久久久久精品久久久久真实原创| 人人澡人人妻人| 毛片一级片免费看久久久久| 亚洲婷婷狠狠爱综合网| 久久国内精品自在自线图片| 曰老女人黄片| 自线自在国产av| 人妻 亚洲 视频| 男女无遮挡免费网站观看| 蜜桃在线观看..| 欧美xxⅹ黑人| 大香蕉久久成人网| 精品一区二区三区四区五区乱码 | 视频区图区小说| 91在线精品国自产拍蜜月| 女的被弄到高潮叫床怎么办| 校园人妻丝袜中文字幕| 国产精品嫩草影院av在线观看| 夜夜爽夜夜爽视频| 亚洲色图 男人天堂 中文字幕 | 美女视频免费永久观看网站| 国产成人精品福利久久| av又黄又爽大尺度在线免费看| 婷婷色av中文字幕| 久久久久国产网址| 亚洲欧美成人综合另类久久久| 性高湖久久久久久久久免费观看| 汤姆久久久久久久影院中文字幕| 丰满迷人的少妇在线观看| 国产深夜福利视频在线观看| av网站免费在线观看视频| 国产成人免费无遮挡视频| 赤兔流量卡办理| 国产xxxxx性猛交| 国产精品成人在线| 18禁在线无遮挡免费观看视频| 亚洲精华国产精华液的使用体验| 亚洲综合色惰| 国产精品嫩草影院av在线观看| 国产日韩欧美在线精品| 成人毛片a级毛片在线播放| 国产成人精品婷婷| 亚洲欧美中文字幕日韩二区| 久久婷婷青草| 国产精品 国内视频| 多毛熟女@视频| 亚洲欧美一区二区三区国产| 日韩一区二区视频免费看| 亚洲精品自拍成人| 中文字幕最新亚洲高清| 免费看av在线观看网站| 久久久久久久久久成人| 婷婷色综合大香蕉| 亚洲高清免费不卡视频| 十八禁高潮呻吟视频| 国产男女超爽视频在线观看| 亚洲av.av天堂| 亚洲欧美中文字幕日韩二区| 久久国产精品大桥未久av| 成人18禁高潮啪啪吃奶动态图| 巨乳人妻的诱惑在线观看| 26uuu在线亚洲综合色| 99热全是精品| 最近最新中文字幕大全免费视频 | 欧美人与善性xxx| 日本vs欧美在线观看视频| 国产老妇伦熟女老妇高清| 香蕉丝袜av| 久久 成人 亚洲| 亚洲国产精品999| 各种免费的搞黄视频| 九色亚洲精品在线播放| 欧美亚洲日本最大视频资源| 中文字幕另类日韩欧美亚洲嫩草| 午夜久久久在线观看| 亚洲欧美日韩另类电影网站| 亚洲,欧美精品.| 五月伊人婷婷丁香| 美女内射精品一级片tv| 国产又爽黄色视频| 大香蕉久久成人网| 91精品伊人久久大香线蕉| av线在线观看网站| 国产av码专区亚洲av| xxx大片免费视频| 亚洲人与动物交配视频| 亚洲国产精品999| 精品一区二区免费观看| 日本欧美国产在线视频| 国产av一区二区精品久久| 日本色播在线视频| 嫩草影院入口| 国产亚洲av片在线观看秒播厂| 国产福利在线免费观看视频| 久久久精品94久久精品| 韩国高清视频一区二区三区| 黄色配什么色好看| av福利片在线| 色婷婷久久久亚洲欧美| 亚洲精华国产精华液的使用体验| 久久精品国产自在天天线| 亚洲精品日本国产第一区| 国产伦理片在线播放av一区| 观看av在线不卡| av电影中文网址| kizo精华| 亚洲美女视频黄频| 在线 av 中文字幕| 捣出白浆h1v1| 我的女老师完整版在线观看| 国产日韩欧美视频二区| 女人久久www免费人成看片| 欧美国产精品一级二级三级| 久久久久久人人人人人| 国产免费一级a男人的天堂| 青青草视频在线视频观看| 国产亚洲av片在线观看秒播厂| 这个男人来自地球电影免费观看 | 国语对白做爰xxxⅹ性视频网站| 人妻 亚洲 视频| 亚洲欧洲国产日韩| 曰老女人黄片| 90打野战视频偷拍视频| 九色成人免费人妻av| 国产精品久久久久久精品古装| 肉色欧美久久久久久久蜜桃| 国产视频首页在线观看| 国产深夜福利视频在线观看| 日韩一本色道免费dvd| 咕卡用的链子| 亚洲性久久影院| 成人漫画全彩无遮挡| h视频一区二区三区| 国产精品人妻久久久影院| 狠狠精品人妻久久久久久综合| av播播在线观看一区| 欧美3d第一页| 免费在线观看黄色视频的| 久久精品国产a三级三级三级| 亚洲美女视频黄频| 国产精品99久久99久久久不卡 | 国产亚洲精品久久久com| 成年美女黄网站色视频大全免费| 国产免费一区二区三区四区乱码| 99re6热这里在线精品视频| 亚洲欧洲精品一区二区精品久久久 | av天堂久久9| 久久这里有精品视频免费| 久久人人爽人人爽人人片va| 男女午夜视频在线观看 | 另类精品久久| 边亲边吃奶的免费视频| 丰满少妇做爰视频| 亚洲中文av在线| 久久97久久精品| 欧美激情 高清一区二区三区| 欧美少妇被猛烈插入视频| 久久国产精品男人的天堂亚洲 | 国产精品国产三级国产专区5o| 亚洲一码二码三码区别大吗| 中文字幕精品免费在线观看视频 | 国产精品欧美亚洲77777| 精品视频人人做人人爽| 一本色道久久久久久精品综合| 久热久热在线精品观看| 超色免费av| 美女国产高潮福利片在线看| 精品熟女少妇av免费看| 亚洲av免费高清在线观看| 99国产综合亚洲精品| 王馨瑶露胸无遮挡在线观看| 看十八女毛片水多多多| av国产精品久久久久影院| 18禁国产床啪视频网站| av一本久久久久| 精品国产一区二区三区四区第35| 五月玫瑰六月丁香| 国产xxxxx性猛交| 大片电影免费在线观看免费| 国产免费又黄又爽又色| 日韩中文字幕视频在线看片| 国产深夜福利视频在线观看| 美女主播在线视频| 狠狠精品人妻久久久久久综合| 国产精品无大码| 欧美亚洲日本最大视频资源| 大香蕉久久成人网| 建设人人有责人人尽责人人享有的| a级毛片在线看网站| 卡戴珊不雅视频在线播放| 91精品三级在线观看| 在线观看国产h片| 两个人免费观看高清视频| 青春草亚洲视频在线观看| 人人妻人人爽人人添夜夜欢视频| 日本av手机在线免费观看| 大香蕉久久网| av国产精品久久久久影院| 久久久久国产网址| 午夜福利视频精品| 亚洲av电影在线进入| 老司机亚洲免费影院| 超色免费av| 内地一区二区视频在线| 成人国产麻豆网|