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

    一種相似性框架下基于非線性擴散過程的剩余壽命估計模型

    2016-04-01 03:52:59裴洪胡昌華王兆強張正新張會會
    中國測試 2016年11期
    關(guān)鍵詞:估計值相似性監(jiān)測點

    裴洪,胡昌華,王兆強,張正新,張會會

    (火箭軍工程大學302室,陜西西安710025)

    一種相似性框架下基于非線性擴散過程的剩余壽命估計模型

    裴洪,胡昌華,王兆強,張正新,張會會

    (火箭軍工程大學302室,陜西西安710025)

    剩余壽命(residual life,RL)估計是預測與健康管理(prognostics and health management,PHM)的關(guān)鍵環(huán)節(jié)。目前,傳統(tǒng)的基于相似性的RL估計模型僅利用失效設(shè)備的退化數(shù)據(jù),忽略截斷設(shè)備和運行設(shè)備的退化數(shù)據(jù),難以保證RL的估計精度。針對該問題,在相似性框架下提出一種新的基于擴散過程的RL估計模型。首先基于擴散過程對截斷設(shè)備進行退化建模和RL估計,然后通過比較參考設(shè)備(失效設(shè)備與截斷設(shè)備)和運行設(shè)備間的相似性,同時利用參考設(shè)備與運行設(shè)備的退化數(shù)據(jù)實現(xiàn)運行設(shè)備的RL估計。最后仿真實驗驗證該文模型優(yōu)于傳統(tǒng)的基于相似性的模型。研究結(jié)果表明:該文模型能夠顯著提高運行設(shè)備RL的估計精度,具有潛在的工程應用價值。

    剩余壽命;擴散過程;相似性;參考設(shè)備

    0 引言

    隨著對工程設(shè)備安全性和可靠性要求的不斷提高,預測與健康管理(PHM)技術(shù)受到了學術(shù)界和工業(yè)界的普遍關(guān)注[1]。作為一種新興技術(shù),PHM能夠提高設(shè)備的可靠性和安全性、降低失效事件發(fā)生的概率,對于航空航天、武器裝備及工業(yè)制造等領(lǐng)域至關(guān)重要[2,3]。作為PHM的重要環(huán)節(jié),剩余壽命(RL)估計可為PHM技術(shù)中決策活動(如最優(yōu)維護、備件訂購等)的有效開展提供前提和基礎(chǔ)。因而,近年來RL估計成為當前研究中的熱點問題,并涌現(xiàn)出了大量的研究成果[4-6]。

    現(xiàn)有的RL估計方法主要可以分為失效物理分析和數(shù)據(jù)驅(qū)動2類方法[1]。失效物理分析方法通過事先獲取設(shè)備的物理失效機制進而實現(xiàn)其RL估計。在實際工程中,對于大型化或復雜化設(shè)備來說,事先獲取物理失效機制的難度很大,同時對于一些價格昂貴的設(shè)備來說,物理模型的獲取成本過高[6]。隨著傳感器技術(shù)和狀態(tài)檢測技術(shù)的發(fā)展,數(shù)據(jù)驅(qū)動的方法凸顯出特有的優(yōu)越性。當前數(shù)據(jù)驅(qū)動的方法大致分為傳統(tǒng)的RL估計方法和基于退化建模的RL估計方法。傳統(tǒng)的RL估計方法通過對設(shè)備的壽命數(shù)據(jù)進行分析進而確定壽命的概率分布。然而,隨著科技的進步,實際設(shè)備的壽命和可靠性不斷提高,通常難以獲取足夠的壽命數(shù)據(jù),尤其是對于一些價格昂貴的設(shè)備來說,獲取壽命數(shù)據(jù)的經(jīng)濟成本過高甚至難以承受,這使得傳統(tǒng)的RL方法在實際中常常難以適用?;谕嘶瘮?shù)據(jù)的方法不僅能夠降低實驗成本,而且能夠提高RL估計的精度,已成為目前的主流方法[4]。事實上,不僅可以從總體中獲得失效設(shè)備或截斷設(shè)備的退化數(shù)據(jù),而且可以獲得運行設(shè)備自身的退化數(shù)據(jù)。如何有效地利用這些數(shù)據(jù)實現(xiàn)設(shè)備RL的估計是當前研究的重點[1-8]。

    由于設(shè)備內(nèi)部結(jié)構(gòu)和使命任務的相似性,同類設(shè)備間的退化過程或退化數(shù)據(jù)必然具有某種相似性。近年來,考慮同類設(shè)備間相似性的退化建模和RL估計方法,即基于相似性的方法,已逐步引起了國內(nèi)外學者的關(guān)注,并取得了較好的實際應用效果。Wang等[9]針對失效設(shè)備足夠多的情況通過構(gòu)建退化數(shù)據(jù)庫,提出了一種基于相似性的RL估計的方法。Zio等[10]基于相似性方法采用隸屬度函數(shù)衡量設(shè)備間相似度,進而確定RL估計值。文獻[11]基于在不同監(jiān)測點運行設(shè)備與參考設(shè)備間相似性所占權(quán)重具有差異性,提出一種廣義的相似性模型進行RL估計,進一步推廣了基于相似性的方法,提高了RL估計精度。文獻[12]在相似性的框架下,利用失效設(shè)備的退化數(shù)據(jù)主要研究了權(quán)重函數(shù)對估計精度的影響、RL估計的魯棒性以及RL估計的不確定性。然而,上述方法均忽略了截斷設(shè)備的退化數(shù)據(jù)和運行設(shè)備的自身退化數(shù)據(jù),在一定程度上無法保證RL估計精度。You等[13]在相似性框架下利用失效設(shè)備和截斷設(shè)備的退化數(shù)據(jù)對RL進行估計,然而該方法采用傳統(tǒng)相似性方法對截斷設(shè)備的RL進行估計,并且將截斷設(shè)備看成失效設(shè)備,具有一定的局限性。鑒于此,本文在相似性框架下提出一種基于非線性擴散過程的RL估計方法。利用參考設(shè)備的退化數(shù)據(jù)和運行設(shè)備的退化數(shù)據(jù),通過比較參考設(shè)備與運行設(shè)備間的相似性確定參考設(shè)備的權(quán)重,與文獻[13]的不同在于基于擴散過程模型確定截斷設(shè)備的RL,進而實現(xiàn)運行設(shè)備的RL估計。數(shù)值仿真驗證了本文模型能夠提高RL的估計精度。

    1 問題描述

    基于相似性的RL估計方法主要通過分析失效設(shè)備與截斷設(shè)備的RL權(quán)重,進而實現(xiàn)運行設(shè)備的RL估計?;谙嗨菩缘腞L估計方法的關(guān)鍵要素如圖1所示。其中,D表示運行設(shè)備和參考設(shè)備間相似性測量的時間范圍,s(·)為測量兩種設(shè)備之間相似性的測度函數(shù),如概率或歐式距離[9],表示在第k個狀態(tài)監(jiān)測點運行設(shè)備與第i個參考設(shè)備間的相似性。權(quán)重函數(shù)w(·)根據(jù)不同參考設(shè)備與運行設(shè)備間的相似性確定其權(quán)重,并將相似性代入權(quán)重函數(shù)中得到權(quán)重w(so?ri(k)),用wi(k)表示。PRL0(k)表示運行設(shè)備在第k個狀態(tài)監(jiān)測點的RL估計值,RRLri(Ni(k))表示第i個參考設(shè)備在第k個狀態(tài)監(jiān)測點的RL,Ni(k)為監(jiān)測點數(shù)目。

    圖1 基于相似性的RL估計方法的關(guān)鍵要素

    傳統(tǒng)的基于相似性的方法均假定參考設(shè)備為失效設(shè)備,而工程設(shè)備的失效往往會造成人員傷亡和財產(chǎn)損失,故有限的失效設(shè)備限制了傳統(tǒng)的基于相似性方法的適用性。文獻[13]同時考慮了失效設(shè)備和截斷設(shè)備的退化數(shù)據(jù),拓寬了參考設(shè)備的范疇,但在實際工程中,截斷設(shè)備即在實際運行過程中未發(fā)生失效的設(shè)備,因而無法預知截斷設(shè)備的實際RL,因而文獻[13]僅簡單地將截斷設(shè)備看成失效設(shè)備,仍采用傳統(tǒng)相似性的方法進行RL估計,難以保證估計的精度。針對該問題,本文主要考慮基于擴散過程的模型對截斷設(shè)備的RL進行估計,進而實現(xiàn)運行設(shè)備的RL估計。

    基于以上分析,本文主要研究以下問題:

    1)如何基于擴散過程模型對截斷設(shè)備的RL進行估計。

    2)如何將步驟1)得到的截斷設(shè)備的RL估計結(jié)果融入相似性的框架內(nèi),利用運行設(shè)備和參考設(shè)備的退化數(shù)據(jù)提高RL估計精度。

    2 相似性框架下基于擴散過程的RL估計

    2.1 截斷設(shè)備的RL估計

    擴散過程模型已廣泛應用于設(shè)備退化模型構(gòu)建和RL估計之中[7,14]。常見的基于擴散過程的退化模型[7]可以表示為

    式中:Y(t)——設(shè)備在t時刻的退化量;

    y0——退化量初始值;

    λ——漂移系數(shù)的尺度參數(shù);

    θ——漂移系數(shù)的形狀參數(shù);

    σB——擴散系數(shù);

    B(t)——標準布朗運動。

    對于一批設(shè)備來說,由于制造工藝及運行環(huán)境的差異性,每個設(shè)備的退化路徑總會存在差異性。因而,令λ為服從正態(tài)分布的隨機參數(shù),即,以刻畫出同類不同個體設(shè)備間的差異性[7,14]。

    設(shè)備的失效時間通常定義為公式(1)所示的退化過程首達既定失效閾值ω的時間,即首達時間[14-15]?;谑走_時間的概念,tk時刻的RL可以定義為,根據(jù)文獻[7],不考慮λ的隨機性時tk時刻RL的概率密度函數(shù)(probability density function,PDF)可以表示為

    可以得到隨機參數(shù)作用下RL的PDF,其中,f(λ)為隨機參數(shù)λ的PDF。

    假定有N組截斷設(shè)備的歷史退化數(shù)據(jù),對于截斷設(shè)備i(i=1,2,…,N),在ti1,ti2,…,tini時刻,測得了設(shè)備i的ni個退化量,記作Yi=(Yi(ti1),Yi(ti2),…,Yi(tini))T。因而,截斷設(shè)備i的退化模型可以表示為

    其中j=1,2,…,ni。令,根據(jù)文獻[7],Yi服從多變量正態(tài)分布,其均值和協(xié)方差分別為

    分別求似然函數(shù)對μλ和σλ的一階偏導數(shù),可以得到如下結(jié)果:

    參數(shù)估計完成后,將進行RL的PDF的推導,對于式(1)定義的非線性隨機退化過程,則tk時刻以w為閾值的RL的PDF可以表示為

    對于第i臺截斷設(shè)備在tR時刻的RL的期望值可表示為,那么如何利用這類數(shù)據(jù),即如何將擴散過程融入相似性框架中,將在下節(jié)討論。

    2.2 相似性框架下運行設(shè)備的RL估計

    圖2總結(jié)了相似性框架下基于擴散過程的RL估計方法。其中,Q為失效設(shè)備的數(shù)目,P為截斷設(shè)備的數(shù)目,SRL(k)為截斷設(shè)備在第k個狀態(tài)監(jiān)測點的RL。A步中,利用截斷設(shè)備的退化數(shù)據(jù),基于擴散過程方法對截斷設(shè)備的RL進行估計。B步中,利用全部退化過程,在該框架下估計運行設(shè)備第k個狀態(tài)監(jiān)測點的RL。從圖2可以看出,H為非負整數(shù)并且(H+1)·Δt=D,即圖1中的時間間隔,B步中,利用連續(xù)H+1個狀態(tài)監(jiān)測點測量參考設(shè)備和運行設(shè)備間的相似性。實際工程中,H或D由專家給定,根據(jù)H+1個監(jiān)測點可以獲取運行設(shè)備的當前狀態(tài),下節(jié)將深入研究H對RL估計的影響。

    圖2 相似性框架下基于擴散過程的RL估計流程

    令y0(k·Δt)表示運行設(shè)備第k個監(jiān)測點的退化量,其中Δt為狀態(tài)監(jiān)測間隔,k為正的監(jiān)測數(shù)。令yri(k·Δt)表示參考設(shè)備i第k個監(jiān)測點的退化量。本文選擇歐式距離函數(shù)測量運行設(shè)備和參考設(shè)備間的相似性:

    式中,s(k,H,i,m)表示運行設(shè)備在[(k-H)·Δt,k·Δt]區(qū)間的退化過程與參考設(shè)備i在[(m-H)· Δt,m·Δt]區(qū)間的退化過程間的相似性,對于運行設(shè)備有k≥H+1,對于失效設(shè)備有H+1≤m≤Mi,對于截斷設(shè)備有H+1≤m≤Gi。其中,Mi表示失效設(shè)備i的監(jiān)測點數(shù)目,Gi表示截斷設(shè)備i的監(jiān)測點數(shù)目。根據(jù)式(7),在第k個監(jiān)測點,參考設(shè)備i與運行設(shè)備間的相似性可以表示為

    其中q,p為正整數(shù)。

    在計算so?ri(k)時,Ni(k)能夠表示為

    第i個參考設(shè)備在Ni(k)監(jiān)測點的RL估計值可以表示為

    根據(jù)文獻[11],權(quán)重函數(shù)可定義為

    其中,zi(k)定義為

    式(16)中,參數(shù)α決定了分配到每臺參考設(shè)備的權(quán)重,文獻[12]研究了α對基于相似性的退化模型的影響。

    根據(jù)式(15)中定義的權(quán)重和式(14)中參考設(shè)備的RL,可以得到運行設(shè)備在第k個狀態(tài)監(jiān)測點的RL估計值:

    通過上述分析,基于圖2中的流程框圖,運行設(shè)備RL估計的算法概括如下。

    步驟1:利用截斷設(shè)備的歷史退化數(shù)據(jù)實現(xiàn)基于擴散過程模型的參數(shù)估計;

    步驟2:基于式(14)利用步驟1的參數(shù)估計結(jié)果實現(xiàn)參考設(shè)備(截斷設(shè)備或失效設(shè)備)i在監(jiān)測點Ni(k)的RL估計;

    步驟3:基于相似性定義和式(15),利用參考設(shè)備的全部退化數(shù)據(jù)與運行設(shè)備自身的退化數(shù)據(jù)確定每個參考設(shè)備相對于運行設(shè)備的權(quán)重wi(k);

    步驟4:基于式(17)能夠確定運行設(shè)備在第k個狀態(tài)監(jiān)測點的RL估計值PRL0(k),實現(xiàn)運行設(shè)備RL的實時估計。

    3 實驗研究

    本節(jié)通過數(shù)值仿真產(chǎn)生10組失效設(shè)備、9組截斷設(shè)備及1組運行設(shè)備的退化數(shù)據(jù),假定退化數(shù)據(jù)單位為mm,利用本文方法與傳統(tǒng)的基于相似性的方法[16]分別估計運行設(shè)備的RL,并比較本文方法與傳統(tǒng)的基于相似性的方法間的RL估計結(jié)果。

    3.1 數(shù)值仿真

    利用式(1)定義的模型產(chǎn)生仿真數(shù)據(jù),參數(shù)設(shè)定如下:μλ=1,σλ=0.1,σB=1,θ=0.9,且Δt=1h。根據(jù)2.1節(jié)提到的首達時間定義,一旦退化數(shù)據(jù)超出了閾值,意味著設(shè)備發(fā)生失效,本文令ω=50 mm。圖3描繪了20組帶有失效閾值的退化過程。

    根據(jù)式(1)隨機仿真得到20組超過失效閾值的退化過程,作為參考設(shè)備和運行設(shè)備的退化過程,對本文方法進行評估。由于仿真得到的數(shù)據(jù)均超出了失效閾值,故為失效設(shè)備的退化過程。為獲取截斷設(shè)備的退化過程,隨機選取9組仿真過程作進一步裁剪,裁剪后的退化過程可作為截斷設(shè)備的退化過程。隨機選取常數(shù)Tr,如Tr=60 h,并且移走Tr與失效時刻間的退化過程,即可獲取截斷設(shè)備的退化過程。根據(jù)截斷設(shè)備的數(shù)目選擇Tr,一旦Tr確定后,設(shè)備的退化過程在監(jiān)測點Tr不會超出閾值。基于圖3的仿真數(shù)據(jù),圖4給出了9組截斷設(shè)備的退化過程。

    圖3 20組帶有閾值的仿真退化軌跡

    圖4 9組截斷退化過程

    3.2 結(jié)果和討論

    首先,將3.1節(jié)中的20組仿真數(shù)據(jù)代入傳統(tǒng)的基于相似性的方法和本文模型中,運行設(shè)備的RL估計值如圖5所示。從圖5可以明顯看出,本文模型和基于傳統(tǒng)相似性方法均能顯示出運行設(shè)備的RL趨勢,但本文模型更加精確。在退化過程的初期,運行設(shè)備的退化數(shù)據(jù)較少,兩種方法的RL估計值均偏離了實際RL,但本文方法的偏離程度相對較小;隨著退化數(shù)據(jù)的不斷獲取,兩種方法的偏離程度逐漸減小,當實際的RL很小時(t=40 h后),兩種方法的PRL非常接近。因為當運行設(shè)備將要失效時,能夠獲得足夠多的運行設(shè)備的退化數(shù)據(jù),能夠準確計算運行設(shè)備與參考設(shè)備間的相似性,因而對于兩種方法,運行設(shè)備的RL估計值PRL(k)與實際RL相差不大。

    圖5 運行設(shè)備的PRL

    圖6 RL估計值的估計誤差

    為了定量比較運行設(shè)備的RL估計精度,這里引入RL估計值的估計誤差(AE)評價RL估計的精度,估計誤差可以表示為

    其中,PRL0(k)表示運行設(shè)備在第k個狀態(tài)監(jiān)測點的RL估計值,ARL0(k)表示設(shè)備在第k個狀態(tài)監(jiān)測點的實際RL。RL估計值的估計誤差如圖6所示??梢钥闯?,在運行設(shè)備的壽命周期內(nèi),本文模型的RL估計誤差均小于基于傳統(tǒng)相似性方法的RL估計誤差,且本文模型的RL估計誤差波動較小,因而本文模型優(yōu)于傳統(tǒng)的基于相似性的方法。

    根據(jù)上述研究,非負整數(shù)與RL估計有著密切關(guān)系。本文選擇非負整數(shù)H=2,5,8,10,13,18,并且利用3.1節(jié)中的20組仿真數(shù)據(jù)研究H的影響。為量化H的影響,RL估計效果G可定義為

    其中,n為運行設(shè)備的全部狀態(tài)監(jiān)測點,ARL0(k)表示設(shè)備在第k個狀態(tài)監(jiān)測點的實際RL。

    根據(jù)式(19),能夠得到RL估計效果,結(jié)果如圖7所示。

    圖7 不同非負整數(shù)H下的RL估計效果

    從圖中能夠直觀看出,H的取值與RL估計效果有著密切關(guān)系。當H值很小時(如H=2),RL估計效果G的值很大,即本文方法效果較差。在一定范圍內(nèi),隨著H的增加,RL的估計效果將明顯提高。因此,當缺乏先驗知識時,H值應該較大一點,但當H值超過一定范圍內(nèi),RL的估計效果將會降低,因為H值選擇太大時,本文方法的復雜性將大大增加。故根據(jù)圖7可得H=8時RL估計效果最佳。

    4 結(jié)束語

    本文針對非線性隨機退化設(shè)備,在相似性的框架下,提出一種基于非線性擴散過程的RL估計方法。首先基于擴散過程構(gòu)建隨機退化模型,利用截斷設(shè)備的退化數(shù)據(jù)實現(xiàn)其RL的估計;然后基于相似性原理確定參考設(shè)備相對于運行設(shè)備的權(quán)重,最后實現(xiàn)運行設(shè)備的RL估計。數(shù)值仿真驗證了本文方法優(yōu)于傳統(tǒng)的基于相似性的方法。研究結(jié)果表明,本文綜合利用失效設(shè)備和截斷設(shè)備的退化數(shù)據(jù),提高運行設(shè)備的RL估計精度,因而在工程實際中,更具有潛在應用價值。

    [1]PECHT M G.Prognostics and Health Management of Electronics[M].John Wiley:New Jersey,2008:201-211.

    [2]司小勝,胡昌華,周東華.帶測量誤差的非線性退化過程建模與剩余壽命估計[J].自動化學報,2013,39 (5):590-601.

    [3]李明福,胡昌華,周志杰,等.基于退化數(shù)據(jù)的貯存設(shè)備最優(yōu)檢測策略[J].系統(tǒng)工程與電子技術(shù),2015,37 (5):1219-1223.

    [4]SI X S,WANG W B,HU C H,et al.Remaining useful life estimation—A review on the statistical data driven approaches[J].European Journal of Operational Research,2011,213(1):1-14.

    [5]WANG Z Q,HU C H,WANG W B,et al.A simulationbased remaining useful life prediction method considering the influence of maintenance activities[C]//Proc.of the Prognostics and System Health Management(PHM) Conference,2014:284-289.

    [6]SHEN Y,TANG L C,XIE M.A model for upside-down bathtub-shaped mean residual life and its properties[J].IEEE Transactions on Reliability,2009,58(2):425-431.

    [7]SI X S,WANG W B,HU C H,et al.Remaining useful life estimation based on a nonlinear diffusion degradation process[J].IEEE Transactions on Reliability,2012,61 (1):50-67.

    [8]WANG W B,CARR M,XU W J,et al.A model for residual life prediction based on Brownian motion with an adaptive drift[J].Microelectronics Reliability,2011,51(1):285-293.

    [9]WANG T,YU J,SIEGEL D,et al.A similarity-based prognostics approach for remaining useful life estimation of engineered systems[C]//Proceedings of International Conference on Prognostics Health Manage,Denver,2008.

    [10]ZIO E,MAIO F D.A data-driven approach for predicting the remaining useful life in dynamic failure scenarios of a nuclear system[J].Reliability Engineering and System Safety,2010,95(1):49-57.

    [11]YOU M Y,MENG G.A generalized similarity measure for similarity-based residual life prediction[J].Proceedings of the Institution of Mechanical Engineers,Part E:Journal of Process Mechanical Engineering,2011,225(3):151-160.

    [12]YOU M Y,MENG G.Toward effective utilization of similarity based residual life prediction methods:Weight allocation,prediction robustness,and prediction uncertainty[J].Proceedings of the Institution of Mechanical Engineers,Part E:Journal of Process Mechanical Engineering,2013,227(1):74-84.

    [13]YOU M Y,MENG G.A framework of similarity-based residual life prediction approaches using degradation histories with failure,preventive maintenance,and suspension events[J].IEEE Transactions on Reliability,2013,62(1):127-135.

    [14]WANG Z Q,HU C H,WANG W B,et al.An Additive Wiener Process-Based Prognostic Model for Hybrid Deteriorating Systems[J].IEEE Transactions on Reliability,2014,63(1):208-222.

    [15]LEE M L T,WHITMORE G A.Threshold regression for survival analysis:Modeling event times by a stochastic process reaching a boundary[J].Statistical Science,2006,21(4):501-513.

    (編輯:李妮)

    A model for residual life estimation based on the nonlinear diffusion process under the framework of similarity

    PEI Hong,HU Changhua,WANG Zhaoqiang,ZHANG Zhengxin,ZHANG Huihui
    (302 Unit,Rocket Force University of Engineering,Xi’an 710025,China)

    Residual life(RL)estimation is a key part in prognostics and health management(PHM).In current literature,traditional similarity-based RL estimation model can only utilize the degradation data of failed devices,while ignoring the degradation data of suspended devices and operating device,which cannot guarantee the accuracy of the RL estimation.Aiming at this issue,this paper develops a new diffusion process driven RL estimation model in the framework of similarity.Firstly,the degradation model is constructed and the RL of the suspended devices is estimated.And then the observed degradation data from the reference devices(i.e.failed and suspended devices)and the operating device are utilized to estimate the RL of the operating device by comparing the similarity between the operating device and reference devices.Finally,a numerical simulation is provided to substantiate the superiority of the proposed model over the traditional similarity-based approach.The research result shows that the proposed model can remarkably improve the accuracy of the RL estimation for the operating device,which can be potentially applied in practice.

    residual life;diffusion process;similarity;reference devices

    A

    1674-5124(2016)11-0006-07

    10.11857/j.issn.1674-5124.2016.11.002

    2016-04-13;

    2016-06-20

    國家杰出青年基金(61025014);國家自然科學基金(61174030,61374120,61573365)

    裴洪(1992-),男,安徽霍邱縣人,碩士研究生,專業(yè)方向為預測與健康管理。

    猜你喜歡
    估計值相似性監(jiān)測點
    一類上三角算子矩陣的相似性與酉相似性
    天津南港LNG接收站沉降監(jiān)測點位布設(shè)
    煤氣與熱力(2022年4期)2022-05-23 12:44:56
    撫河流域綜合治理監(jiān)測布局優(yōu)化
    淺析當代中西方繪畫的相似性
    河北畫報(2020年8期)2020-10-27 02:54:20
    全站儀極坐標法監(jiān)測點穩(wěn)定性分析方法研究
    一道樣本的數(shù)字特征與頻率分布直方圖的交匯問題
    統(tǒng)計信息
    2018年4月世界粗鋼產(chǎn)量表(續(xù))萬噸
    低滲透黏土中氯離子彌散作用離心模擬相似性
    我省舉辦家畜血吸蟲病監(jiān)測點培訓班
    毛片一级片免费看久久久久| 舔av片在线| 日韩伦理黄色片| 欧美日韩精品成人综合77777| 欧美激情国产日韩精品一区| 亚洲精品一二三| 日本-黄色视频高清免费观看| 亚洲欧美日韩无卡精品| 一个人免费在线观看电影| 亚洲国产色片| 国内精品宾馆在线| 黄片wwwwww| 淫秽高清视频在线观看| 精品人妻一区二区三区麻豆| 又粗又硬又长又爽又黄的视频| 美女高潮的动态| 女人十人毛片免费观看3o分钟| 嫩草影院精品99| 男女边吃奶边做爰视频| 色吧在线观看| 色哟哟·www| 国产亚洲午夜精品一区二区久久 | 在线a可以看的网站| av在线亚洲专区| 不卡视频在线观看欧美| 免费看光身美女| 精品不卡国产一区二区三区| 中文资源天堂在线| 又爽又黄无遮挡网站| 亚洲国产色片| 久久这里只有精品中国| 女人久久www免费人成看片| 亚洲天堂国产精品一区在线| 最近最新中文字幕大全电影3| 九九爱精品视频在线观看| 大香蕉97超碰在线| 一本久久精品| 高清欧美精品videossex| 国产精品三级大全| 亚洲国产av新网站| 老司机影院毛片| 欧美三级亚洲精品| 亚洲精品视频女| 男人舔女人下体高潮全视频| 97超视频在线观看视频| 综合色av麻豆| 一二三四中文在线观看免费高清| 亚洲av电影不卡..在线观看| 国产精品一区二区在线观看99 | 色哟哟·www| 插阴视频在线观看视频| 国产成人精品一,二区| 色网站视频免费| 国内揄拍国产精品人妻在线| 性插视频无遮挡在线免费观看| 成人二区视频| 欧美xxxx黑人xx丫x性爽| 国内少妇人妻偷人精品xxx网站| 人妻夜夜爽99麻豆av| 超碰av人人做人人爽久久| 18+在线观看网站| 成人午夜精彩视频在线观看| 亚洲在久久综合| 中文天堂在线官网| 美女国产视频在线观看| 一级毛片电影观看| 校园人妻丝袜中文字幕| 国产一区二区在线观看日韩| 免费不卡的大黄色大毛片视频在线观看 | 亚洲欧美一区二区三区国产| 人妻少妇偷人精品九色| 99热6这里只有精品| 免费看av在线观看网站| 亚洲成人久久爱视频| 熟女人妻精品中文字幕| 亚洲最大成人av| 亚洲自偷自拍三级| 精品久久国产蜜桃| 久久精品久久精品一区二区三区| 日本一二三区视频观看| 国产中年淑女户外野战色| 久久久久久久大尺度免费视频| 真实男女啪啪啪动态图| 国产午夜精品论理片| 欧美日韩一区二区视频在线观看视频在线 | 大香蕉97超碰在线| 欧美日韩亚洲高清精品| 午夜久久久久精精品| 久久久久久久亚洲中文字幕| 免费看av在线观看网站| 亚洲av不卡在线观看| 69av精品久久久久久| 午夜精品一区二区三区免费看| 少妇人妻精品综合一区二区| 日本黄色片子视频| 十八禁国产超污无遮挡网站| 熟妇人妻久久中文字幕3abv| 免费观看av网站的网址| 成人毛片60女人毛片免费| 插逼视频在线观看| 夜夜看夜夜爽夜夜摸| 中文字幕av成人在线电影| 国内揄拍国产精品人妻在线| 激情五月婷婷亚洲| 97热精品久久久久久| 青春草国产在线视频| 国产高清三级在线| 亚洲在久久综合| 一区二区三区高清视频在线| 最新中文字幕久久久久| 亚洲av一区综合| 久久久成人免费电影| 青春草亚洲视频在线观看| 久99久视频精品免费| 人人妻人人看人人澡| 亚洲欧美日韩无卡精品| 2021少妇久久久久久久久久久| 日韩视频在线欧美| 亚洲精品乱码久久久久久按摩| 综合色av麻豆| 国产单亲对白刺激| 99久久九九国产精品国产免费| 亚洲怡红院男人天堂| 精品久久久精品久久久| 国产色婷婷99| 汤姆久久久久久久影院中文字幕 | 一区二区三区免费毛片| 一二三四中文在线观看免费高清| 亚洲第一区二区三区不卡| videossex国产| 日韩 亚洲 欧美在线| 亚洲国产色片| 91精品伊人久久大香线蕉| 99热网站在线观看| 成人特级av手机在线观看| 天堂√8在线中文| 免费看a级黄色片| 亚洲熟女精品中文字幕| 免费不卡的大黄色大毛片视频在线观看 | 观看美女的网站| 最近手机中文字幕大全| 亚洲精品456在线播放app| 看黄色毛片网站| 寂寞人妻少妇视频99o| 国产不卡一卡二| 偷拍熟女少妇极品色| 成年人午夜在线观看视频 | 美女脱内裤让男人舔精品视频| 亚洲成人精品中文字幕电影| 麻豆国产97在线/欧美| 一本久久精品| 丝瓜视频免费看黄片| 一个人看视频在线观看www免费| 97人妻精品一区二区三区麻豆| 欧美成人a在线观看| 欧美一级a爱片免费观看看| 亚洲精品乱码久久久久久按摩| 69人妻影院| 精品久久久久久久久久久久久| 麻豆久久精品国产亚洲av| 淫秽高清视频在线观看| 国产 一区 欧美 日韩| 在线观看av片永久免费下载| 久久6这里有精品| 天堂影院成人在线观看| 免费黄频网站在线观看国产| 黄色一级大片看看| 国产精品国产三级国产av玫瑰| 毛片女人毛片| 国产成人freesex在线| 99热网站在线观看| 免费在线观看成人毛片| 国产有黄有色有爽视频| 嫩草影院精品99| 国产91av在线免费观看| 日韩精品青青久久久久久| 欧美zozozo另类| 欧美一级a爱片免费观看看| 亚洲真实伦在线观看| 丝袜美腿在线中文| a级一级毛片免费在线观看| 亚洲熟女精品中文字幕| 在线免费十八禁| 特级一级黄色大片| 国产成人精品婷婷| 高清毛片免费看| av网站免费在线观看视频 | 大又大粗又爽又黄少妇毛片口| av在线观看视频网站免费| 日本免费a在线| 国产 一区 欧美 日韩| 日韩精品青青久久久久久| 久久久精品欧美日韩精品| 中文字幕免费在线视频6| 亚洲欧美一区二区三区黑人 | 白带黄色成豆腐渣| 久久人人爽人人爽人人片va| 亚洲精品中文字幕在线视频 | 视频中文字幕在线观看| 亚洲一区高清亚洲精品| 纵有疾风起免费观看全集完整版 | 日韩电影二区| 亚洲欧洲国产日韩| 韩国av在线不卡| 欧美激情国产日韩精品一区| 男女边吃奶边做爰视频| 国产精品伦人一区二区| 日韩欧美一区视频在线观看 | 精品一区二区三区视频在线| 白带黄色成豆腐渣| 偷拍熟女少妇极品色| 老师上课跳d突然被开到最大视频| 男女边摸边吃奶| 欧美日韩视频高清一区二区三区二| 高清在线视频一区二区三区| 天天躁夜夜躁狠狠久久av| 欧美一区二区亚洲| 亚洲丝袜综合中文字幕| 久久99热这里只有精品18| 能在线免费观看的黄片| 亚洲成色77777| 综合色av麻豆| 亚洲三级黄色毛片| 国产成人免费观看mmmm| 日韩一区二区三区影片| 国产av码专区亚洲av| 免费大片18禁| 亚洲人成网站在线观看播放| 亚洲精品色激情综合| 日韩一区二区三区影片| 国产精品综合久久久久久久免费| 插阴视频在线观看视频| 免费无遮挡裸体视频| 精品久久久久久久久av| 日日啪夜夜爽| 五月伊人婷婷丁香| 99热这里只有是精品在线观看| 黄色一级大片看看| 99九九线精品视频在线观看视频| 精品一区二区三卡| 亚洲怡红院男人天堂| 午夜激情福利司机影院| 国产精品久久久久久久电影| 国产黄色小视频在线观看| 亚洲在线自拍视频| 久久99热这里只有精品18| 亚洲精品成人av观看孕妇| 熟女电影av网| 久久综合国产亚洲精品| 免费人成在线观看视频色| 日本色播在线视频| 国产精品久久久久久精品电影小说 | 色视频www国产| 国产成人精品久久久久久| 校园人妻丝袜中文字幕| 国产精品不卡视频一区二区| 精品久久久久久久久av| 成人美女网站在线观看视频| 久久这里只有精品中国| 中文字幕制服av| 男人舔奶头视频| 99久久精品国产国产毛片| 岛国毛片在线播放| 精品一区在线观看国产| 国产精品麻豆人妻色哟哟久久 | 又黄又爽又刺激的免费视频.| 国产 一区精品| 亚洲国产色片| 日本黄大片高清| 欧美zozozo另类| 欧美高清成人免费视频www| 免费电影在线观看免费观看| 99久国产av精品| 在线观看av片永久免费下载| 国内精品宾馆在线| 高清午夜精品一区二区三区| 91在线精品国自产拍蜜月| 亚洲人成网站在线观看播放| 97超视频在线观看视频| 午夜久久久久精精品| 欧美日韩精品成人综合77777| 少妇人妻精品综合一区二区| 国内少妇人妻偷人精品xxx网站| 18禁裸乳无遮挡免费网站照片| 日韩,欧美,国产一区二区三区| 亚洲国产欧美人成| 成年女人在线观看亚洲视频 | 99久久精品国产国产毛片| 色综合色国产| 69av精品久久久久久| 中国美白少妇内射xxxbb| 国产亚洲91精品色在线| 亚洲在线观看片| 自拍偷自拍亚洲精品老妇| 成人av在线播放网站| 国产老妇女一区| 成人午夜高清在线视频| 内射极品少妇av片p| av免费在线看不卡| 精品人妻一区二区三区麻豆| 亚洲伊人久久精品综合| 亚洲精品乱码久久久v下载方式| 日韩欧美精品v在线| 国产精品久久久久久精品电影小说 | 直男gayav资源| 亚洲成人中文字幕在线播放| 国产精品久久久久久久久免| 亚洲熟女精品中文字幕| 亚洲精品中文字幕在线视频 | 国产精品久久久久久精品电影小说 | 天堂影院成人在线观看| 婷婷六月久久综合丁香| 国产精品久久久久久久久免| 国产亚洲5aaaaa淫片| 精品一区二区三区视频在线| 午夜福利网站1000一区二区三区| 国产精品久久久久久av不卡| 网址你懂的国产日韩在线| 一夜夜www| 赤兔流量卡办理| 草草在线视频免费看| 日韩一区二区三区影片| 亚洲真实伦在线观看| 51国产日韩欧美| 国产女主播在线喷水免费视频网站 | 午夜免费男女啪啪视频观看| 亚洲国产欧美人成| 波多野结衣巨乳人妻| 九色成人免费人妻av| 久久久久久久亚洲中文字幕| 日韩三级伦理在线观看| 国产大屁股一区二区在线视频| 亚洲精品日韩av片在线观看| 麻豆精品久久久久久蜜桃| 国产女主播在线喷水免费视频网站 | 亚洲在久久综合| 听说在线观看完整版免费高清| 精品久久久久久久人妻蜜臀av| 2018国产大陆天天弄谢| 99热这里只有是精品50| 日韩欧美三级三区| 男女那种视频在线观看| 亚洲av免费在线观看| 久久精品国产鲁丝片午夜精品| 久久久久精品久久久久真实原创| 成人鲁丝片一二三区免费| 神马国产精品三级电影在线观看| av女优亚洲男人天堂| 寂寞人妻少妇视频99o| 久热久热在线精品观看| 久久精品人妻少妇| 亚洲人与动物交配视频| 寂寞人妻少妇视频99o| 精品久久久精品久久久| 亚洲第一区二区三区不卡| 欧美日韩视频高清一区二区三区二| 97超碰精品成人国产| av又黄又爽大尺度在线免费看| 亚洲久久久久久中文字幕| 少妇丰满av| 91精品国产九色| 日韩人妻高清精品专区| 少妇熟女欧美另类| 中文欧美无线码| 国产69精品久久久久777片| 亚洲精品,欧美精品| 亚洲精品第二区| 中文字幕av成人在线电影| 久99久视频精品免费| av网站免费在线观看视频 | 人体艺术视频欧美日本| 国产精品久久视频播放| 肉色欧美久久久久久久蜜桃 | 老女人水多毛片| 日韩大片免费观看网站| 亚洲在线观看片| 亚洲熟妇中文字幕五十中出| 一级毛片aaaaaa免费看小| 免费av毛片视频| 国产在视频线在精品| 26uuu在线亚洲综合色| 久久久久久久国产电影| 美女国产视频在线观看| 老师上课跳d突然被开到最大视频| 久久久久久久久大av| 又爽又黄a免费视频| 我的女老师完整版在线观看| 国产成人a区在线观看| 国产亚洲91精品色在线| 中文字幕av在线有码专区| 嫩草影院精品99| 尾随美女入室| 偷拍熟女少妇极品色| 日日摸夜夜添夜夜添av毛片| 午夜亚洲福利在线播放| 网址你懂的国产日韩在线| 久久午夜福利片| 亚洲经典国产精华液单| 最后的刺客免费高清国语| 亚洲国产高清在线一区二区三| 久久久久九九精品影院| 亚洲精品aⅴ在线观看| 国产片特级美女逼逼视频| ponron亚洲| 在线免费十八禁| 男的添女的下面高潮视频| 简卡轻食公司| 国产黄频视频在线观看| 特大巨黑吊av在线直播| 日韩制服骚丝袜av| 亚洲天堂国产精品一区在线| 亚洲国产欧美在线一区| 人人妻人人看人人澡| 偷拍熟女少妇极品色| 日本熟妇午夜| 视频中文字幕在线观看| 十八禁网站网址无遮挡 | 色网站视频免费| 禁无遮挡网站| 纵有疾风起免费观看全集完整版 | 欧美高清成人免费视频www| 国产精品日韩av在线免费观看| 看黄色毛片网站| 免费大片18禁| 人人妻人人澡人人爽人人夜夜 | 国产精品不卡视频一区二区| 两个人视频免费观看高清| 99久久九九国产精品国产免费| 欧美成人一区二区免费高清观看| 在线天堂最新版资源| 成年人午夜在线观看视频 | 校园人妻丝袜中文字幕| 三级男女做爰猛烈吃奶摸视频| 一级毛片我不卡| 免费观看在线日韩| 亚洲国产精品成人久久小说| 日韩av不卡免费在线播放| 国产老妇伦熟女老妇高清| 国产精品综合久久久久久久免费| 久久精品久久久久久久性| 国产精品福利在线免费观看| a级毛色黄片| 精品久久久久久久久亚洲| 国产精品蜜桃在线观看| 亚洲精品视频女| 国产老妇伦熟女老妇高清| 天天躁夜夜躁狠狠久久av| 美女高潮的动态| 又大又黄又爽视频免费| 亚洲欧美中文字幕日韩二区| 亚洲,欧美,日韩| 一个人免费在线观看电影| 国产成人一区二区在线| 一级a做视频免费观看| 看黄色毛片网站| or卡值多少钱| 午夜视频国产福利| 淫秽高清视频在线观看| 特级一级黄色大片| 一级毛片我不卡| 男的添女的下面高潮视频| 麻豆av噜噜一区二区三区| 美女大奶头视频| av线在线观看网站| av又黄又爽大尺度在线免费看| 天堂av国产一区二区熟女人妻| 欧美日韩视频高清一区二区三区二| 网址你懂的国产日韩在线| 精品久久久精品久久久| 亚洲精品亚洲一区二区| 99热这里只有是精品在线观看| 国产一区二区在线观看日韩| 天堂影院成人在线观看| 91久久精品国产一区二区三区| 欧美激情国产日韩精品一区| 亚洲综合精品二区| 国产精品综合久久久久久久免费| 国产精品伦人一区二区| 亚洲一级一片aⅴ在线观看| 午夜免费观看性视频| 乱人视频在线观看| 成人高潮视频无遮挡免费网站| 能在线免费看毛片的网站| 国内少妇人妻偷人精品xxx网站| 亚洲国产最新在线播放| 纵有疾风起免费观看全集完整版 | 高清午夜精品一区二区三区| 亚洲精品第二区| 国产一区二区在线观看日韩| 国产免费又黄又爽又色| or卡值多少钱| 亚洲av二区三区四区| 男人和女人高潮做爰伦理| 免费高清在线观看视频在线观看| 国产久久久一区二区三区| 国产国拍精品亚洲av在线观看| 人妻一区二区av| 美女内射精品一级片tv| 白带黄色成豆腐渣| 午夜爱爱视频在线播放| 26uuu在线亚洲综合色| 亚洲人成网站高清观看| 亚洲av一区综合| 在线天堂最新版资源| 久久精品夜色国产| 国产午夜精品久久久久久一区二区三区| 伊人久久国产一区二区| 少妇的逼好多水| 成人高潮视频无遮挡免费网站| 最近手机中文字幕大全| 亚洲第一区二区三区不卡| 麻豆av噜噜一区二区三区| 精品一区二区免费观看| 久久久久久久久中文| 亚洲aⅴ乱码一区二区在线播放| 永久网站在线| 女人被狂操c到高潮| 亚洲综合色惰| 午夜免费男女啪啪视频观看| 好男人视频免费观看在线| 人体艺术视频欧美日本| 99热全是精品| 又大又黄又爽视频免费| 久久久精品免费免费高清| 国产黄片视频在线免费观看| 亚洲欧美一区二区三区黑人 | 欧美精品国产亚洲| 欧美高清成人免费视频www| 亚洲婷婷狠狠爱综合网| 亚洲精品自拍成人| 国产69精品久久久久777片| 五月伊人婷婷丁香| 精品国内亚洲2022精品成人| 亚洲内射少妇av| 黄色欧美视频在线观看| 精品久久国产蜜桃| 18+在线观看网站| 免费观看性生交大片5| 亚洲人成网站在线观看播放| 成人高潮视频无遮挡免费网站| 亚洲精品久久午夜乱码| 国产 亚洲一区二区三区 | 又爽又黄无遮挡网站| 成人国产麻豆网| 真实男女啪啪啪动态图| 九九久久精品国产亚洲av麻豆| 老师上课跳d突然被开到最大视频| 色吧在线观看| 亚洲三级黄色毛片| 精品一区在线观看国产| 十八禁网站网址无遮挡 | av在线老鸭窝| 国产成人免费观看mmmm| av在线老鸭窝| 欧美另类一区| 色5月婷婷丁香| 91精品伊人久久大香线蕉| 国产高清不卡午夜福利| 五月天丁香电影| 黄色日韩在线| 好男人在线观看高清免费视频| 免费观看在线日韩| 简卡轻食公司| 欧美日韩在线观看h| 99热这里只有精品一区| av免费在线看不卡| 精品久久久久久久人妻蜜臀av| 国国产精品蜜臀av免费| 国产乱人偷精品视频| 亚洲18禁久久av| 人人妻人人看人人澡| 亚洲av一区综合| 夫妻性生交免费视频一级片| 男插女下体视频免费在线播放| 亚洲国产高清在线一区二区三| 简卡轻食公司| 汤姆久久久久久久影院中文字幕 | 国产精品.久久久| 丰满少妇做爰视频| 国产av码专区亚洲av| av女优亚洲男人天堂| freevideosex欧美| 久热久热在线精品观看| 日本wwww免费看| 亚洲三级黄色毛片| 色综合站精品国产| 日本熟妇午夜| 国产精品久久视频播放| 我的老师免费观看完整版| 狂野欧美白嫩少妇大欣赏| 亚洲精品日韩av片在线观看| 国产乱人偷精品视频| 亚洲av成人av| 亚洲性久久影院| 国产精品久久久久久av不卡| 国内精品美女久久久久久| 99热6这里只有精品| 久久精品国产自在天天线| 精品人妻一区二区三区麻豆| 国产黄a三级三级三级人| 日韩在线高清观看一区二区三区| 两个人视频免费观看高清| 日韩一本色道免费dvd| 免费黄频网站在线观看国产| 久久6这里有精品| 亚洲精品乱码久久久久久按摩| av在线老鸭窝| 在线观看免费高清a一片| 自拍偷自拍亚洲精品老妇| 国产淫语在线视频| 大片免费播放器 马上看| 麻豆久久精品国产亚洲av| 一级av片app| 日日撸夜夜添| 看黄色毛片网站|