• <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)測點培訓班
    久99久视频精品免费| 亚洲熟妇熟女久久| 在线观看66精品国产| 少妇 在线观看| 老司机靠b影院| 午夜福利高清视频| 中文亚洲av片在线观看爽| 亚洲欧美精品综合一区二区三区| 久久国产精品影院| 亚洲人成网站在线播放欧美日韩| 侵犯人妻中文字幕一二三四区| 欧美精品亚洲一区二区| 妹子高潮喷水视频| 精品乱码久久久久久99久播| 真人一进一出gif抽搐免费| 午夜久久久久精精品| 精品福利观看| 国产欧美日韩一区二区三| 欧美在线一区亚洲| 亚洲avbb在线观看| 黄色成人免费大全| 99热只有精品国产| 午夜精品久久久久久毛片777| 国产人伦9x9x在线观看| 久久久久久免费高清国产稀缺| 日韩欧美一区二区三区在线观看| 18禁裸乳无遮挡免费网站照片 | 亚洲国产欧美日韩在线播放| 美女高潮到喷水免费观看| 制服诱惑二区| 丁香六月欧美| 免费看a级黄色片| 亚洲人成77777在线视频| 中文字幕人妻熟女乱码| 真人做人爱边吃奶动态| 国产一区二区三区视频了| 黄色女人牲交| 亚洲精品久久国产高清桃花| 夜夜躁狠狠躁天天躁| 亚洲色图av天堂| 男人舔女人下体高潮全视频| 1024视频免费在线观看| 黄色a级毛片大全视频| 亚洲欧美精品综合一区二区三区| 免费在线观看完整版高清| 男女视频在线观看网站免费 | 熟妇人妻久久中文字幕3abv| 美女午夜性视频免费| а√天堂www在线а√下载| 欧美中文日本在线观看视频| 精品福利观看| 中文字幕人妻丝袜一区二区| 麻豆国产av国片精品| 亚洲人成网站高清观看| 精品电影一区二区在线| 极品教师在线免费播放| 夜夜夜夜夜久久久久| 亚洲午夜理论影院| 国产精品影院久久| 亚洲欧美一区二区三区黑人| 成人18禁在线播放| 精品国产亚洲在线| 国产精品九九99| 欧美日韩福利视频一区二区| 91九色精品人成在线观看| 日本五十路高清| 成在线人永久免费视频| 亚洲一区高清亚洲精品| 美女国产高潮福利片在线看| 亚洲第一青青草原| 亚洲成人久久爱视频| 久久国产亚洲av麻豆专区| ponron亚洲| 琪琪午夜伦伦电影理论片6080| 免费在线观看黄色视频的| 丁香六月欧美| 欧美又色又爽又黄视频| 嫩草影院精品99| 欧美色欧美亚洲另类二区| 国产亚洲精品av在线| 国产三级在线视频| 免费看美女性在线毛片视频| www国产在线视频色| 欧美国产精品va在线观看不卡| 一级a爱片免费观看的视频| 日本a在线网址| 午夜福利成人在线免费观看| 亚洲色图av天堂| 色综合亚洲欧美另类图片| 一边摸一边抽搐一进一小说| 制服人妻中文乱码| 亚洲国产高清在线一区二区三 | 欧美中文综合在线视频| 精品福利观看| 国产欧美日韩一区二区三| 国产精华一区二区三区| 国内精品久久久久精免费| 国产爱豆传媒在线观看 | 日韩一卡2卡3卡4卡2021年| 真人做人爱边吃奶动态| 国产av在哪里看| 国产黄色小视频在线观看| 色av中文字幕| 黄色成人免费大全| 美女高潮喷水抽搐中文字幕| 麻豆国产av国片精品| 成人精品一区二区免费| 女人高潮潮喷娇喘18禁视频| 国产伦人伦偷精品视频| 熟女少妇亚洲综合色aaa.| 777久久人妻少妇嫩草av网站| 香蕉av资源在线| 久久性视频一级片| 国产激情欧美一区二区| 日本免费一区二区三区高清不卡| 国产av不卡久久| 欧美激情 高清一区二区三区| 国产一区二区激情短视频| 久久这里只有精品19| 成人精品一区二区免费| 村上凉子中文字幕在线| 午夜亚洲福利在线播放| 亚洲avbb在线观看| 成年版毛片免费区| 午夜福利欧美成人| 国产精品,欧美在线| 人人妻人人澡欧美一区二区| 此物有八面人人有两片| 国产三级在线视频| 天天躁夜夜躁狠狠躁躁| 又大又爽又粗| 欧美zozozo另类| 午夜福利在线观看吧| 亚洲人成伊人成综合网2020| 夜夜躁狠狠躁天天躁| 岛国在线观看网站| 18禁观看日本| 国产免费av片在线观看野外av| 俄罗斯特黄特色一大片| 欧美色视频一区免费| 国产精品永久免费网站| 日韩大尺度精品在线看网址| 亚洲五月色婷婷综合| 麻豆成人av在线观看| 国内揄拍国产精品人妻在线 | 精品久久久久久成人av| 亚洲成人免费电影在线观看| 美女高潮到喷水免费观看| 首页视频小说图片口味搜索| 免费无遮挡裸体视频| 1024香蕉在线观看| 欧美成狂野欧美在线观看| 9191精品国产免费久久| 亚洲第一av免费看| 嫁个100分男人电影在线观看| 国产精品野战在线观看| 久久国产乱子伦精品免费另类| 一个人观看的视频www高清免费观看 | 男人舔奶头视频| 精品第一国产精品| 亚洲av五月六月丁香网| 午夜福利欧美成人| 十八禁人妻一区二区| 在线播放国产精品三级| 亚洲成人久久性| 免费看美女性在线毛片视频| 久久久久精品国产欧美久久久| 不卡一级毛片| 美女高潮到喷水免费观看| 成人手机av| 看片在线看免费视频| 亚洲成人精品中文字幕电影| 午夜福利在线观看吧| 国产亚洲欧美精品永久| 国产一级毛片七仙女欲春2 | 麻豆成人av在线观看| 男女床上黄色一级片免费看| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲自偷自拍图片 自拍| 亚洲激情在线av| 在线观看免费午夜福利视频| 日本熟妇午夜| 免费电影在线观看免费观看| 99久久国产精品久久久| 亚洲午夜精品一区,二区,三区| 天堂动漫精品| 国产单亲对白刺激| 日韩中文字幕欧美一区二区| 亚洲av电影在线进入| svipshipincom国产片| ponron亚洲| 成人午夜高清在线视频 | 91av网站免费观看| 亚洲第一电影网av| ponron亚洲| 亚洲国产精品成人综合色| 欧美一级a爱片免费观看看 | 久久久久久久精品吃奶| 久久青草综合色| 村上凉子中文字幕在线| 制服丝袜大香蕉在线| xxx96com| 91成人精品电影| 日韩高清综合在线| 韩国av一区二区三区四区| 欧美午夜高清在线| x7x7x7水蜜桃| 久热爱精品视频在线9| 成人免费观看视频高清| 国产精品 欧美亚洲| 男人舔奶头视频| 1024视频免费在线观看| 国产亚洲精品av在线| 午夜免费观看网址| 日本撒尿小便嘘嘘汇集6| 久久香蕉精品热| 久久人妻福利社区极品人妻图片| 成人18禁高潮啪啪吃奶动态图| 国产精品久久久久久亚洲av鲁大| 视频在线观看一区二区三区| 9191精品国产免费久久| 宅男免费午夜| 久久久久久亚洲精品国产蜜桃av| 精品一区二区三区视频在线观看免费| 精品国产乱码久久久久久男人| 久久热在线av| 亚洲国产精品sss在线观看| 久久性视频一级片| 国产精品乱码一区二三区的特点| 黄片小视频在线播放| АⅤ资源中文在线天堂| 俄罗斯特黄特色一大片| 首页视频小说图片口味搜索| 中文字幕精品免费在线观看视频| 中文字幕人成人乱码亚洲影| av视频在线观看入口| av片东京热男人的天堂| 亚洲av片天天在线观看| 悠悠久久av| 在线免费观看的www视频| 久久久久九九精品影院| 丰满人妻熟妇乱又伦精品不卡| 亚洲精品一卡2卡三卡4卡5卡| 久热这里只有精品99| 大型av网站在线播放| 亚洲一区高清亚洲精品| 国产精品野战在线观看| 嫩草影视91久久| 国产亚洲欧美在线一区二区| 精品无人区乱码1区二区| 亚洲午夜理论影院| www日本黄色视频网| 一级片免费观看大全| 一本精品99久久精品77| 不卡一级毛片| 中文字幕人妻丝袜一区二区| 男人操女人黄网站| 最近最新中文字幕大全电影3 | 欧美日韩中文字幕国产精品一区二区三区| 国产极品粉嫩免费观看在线| 91九色精品人成在线观看| 最好的美女福利视频网| 51午夜福利影视在线观看| 一进一出好大好爽视频| 精品久久久久久久末码| 一区二区三区精品91| 亚洲精品国产一区二区精华液| 老司机深夜福利视频在线观看| 久久精品人妻少妇| АⅤ资源中文在线天堂| 欧美人与性动交α欧美精品济南到| 久久青草综合色| 精品人妻1区二区| 丝袜人妻中文字幕| 9191精品国产免费久久| 丁香欧美五月| 亚洲五月天丁香| 一本精品99久久精品77| 悠悠久久av| 亚洲av电影在线进入| 亚洲免费av在线视频| 搞女人的毛片| 一区二区日韩欧美中文字幕| a在线观看视频网站| 久久人妻av系列| 久久精品人妻少妇| 亚洲七黄色美女视频| 精品久久久久久久久久久久久 | 久久人人精品亚洲av| 男人的好看免费观看在线视频 | 婷婷六月久久综合丁香| 国产视频内射| 欧美一区二区精品小视频在线| 18禁裸乳无遮挡免费网站照片 | 中文字幕另类日韩欧美亚洲嫩草| 国产99久久九九免费精品| 国产蜜桃级精品一区二区三区| 国产av一区二区精品久久| 首页视频小说图片口味搜索| 韩国精品一区二区三区| 欧美黑人巨大hd| 国产aⅴ精品一区二区三区波| 亚洲精品在线美女| 国产私拍福利视频在线观看| 桃红色精品国产亚洲av| 中文字幕av电影在线播放| 2021天堂中文幕一二区在线观 | 成年版毛片免费区| 免费电影在线观看免费观看| 一进一出抽搐动态| 国产伦在线观看视频一区| 日韩欧美一区视频在线观看| 久久九九热精品免费| netflix在线观看网站| 一a级毛片在线观看| 好男人在线观看高清免费视频 | 香蕉丝袜av| 国产又黄又爽又无遮挡在线| 国产av在哪里看| 久99久视频精品免费| 999久久久精品免费观看国产| 在线观看午夜福利视频| 亚洲av五月六月丁香网| 久久精品国产亚洲av高清一级| 国产亚洲精品一区二区www| 琪琪午夜伦伦电影理论片6080| 夜夜躁狠狠躁天天躁| 免费电影在线观看免费观看| 高清在线国产一区| 一区二区日韩欧美中文字幕| 一卡2卡三卡四卡精品乱码亚洲| 免费看a级黄色片| 中文字幕精品亚洲无线码一区 | 91国产中文字幕| 精品不卡国产一区二区三区| 一个人观看的视频www高清免费观看 | 大型av网站在线播放| 免费在线观看视频国产中文字幕亚洲| or卡值多少钱| 中文字幕高清在线视频| 久久久久久久午夜电影| 一级毛片精品| 国产黄a三级三级三级人| 欧美一级a爱片免费观看看 | 久9热在线精品视频| 亚洲成人免费电影在线观看| 又大又爽又粗| 久热爱精品视频在线9| 久久亚洲真实| 手机成人av网站| 国产成人精品久久二区二区91| 亚洲va日本ⅴa欧美va伊人久久| 欧美激情高清一区二区三区| 搡老熟女国产l中国老女人| 国产一区二区三区在线臀色熟女| 精品久久久久久久人妻蜜臀av| 天堂√8在线中文| 精品一区二区三区av网在线观看| 99热只有精品国产| 国产在线精品亚洲第一网站| 国产一卡二卡三卡精品| 成人18禁高潮啪啪吃奶动态图| 1024香蕉在线观看| 国产又黄又爽又无遮挡在线| 国产精品,欧美在线| 欧美大码av| 宅男免费午夜| 精品第一国产精品| 久久久久久久久久黄片| 黄片播放在线免费| 亚洲国产欧洲综合997久久, | 91av网站免费观看| 搞女人的毛片| 欧美成人性av电影在线观看| 日韩精品中文字幕看吧| 亚洲七黄色美女视频| 欧美精品啪啪一区二区三区| а√天堂www在线а√下载| 在线十欧美十亚洲十日本专区| 香蕉丝袜av| 亚洲自偷自拍图片 自拍| 一夜夜www| www日本在线高清视频| 男女之事视频高清在线观看| 欧美精品啪啪一区二区三区| 亚洲国产精品成人综合色| 精品福利观看| 在线观看午夜福利视频| 黄网站色视频无遮挡免费观看| 国产亚洲av嫩草精品影院| 亚洲精品国产区一区二| 麻豆成人午夜福利视频| 亚洲专区中文字幕在线| 午夜激情福利司机影院| 国产野战对白在线观看| 中出人妻视频一区二区| 亚洲国产看品久久| 日本在线视频免费播放| 18禁国产床啪视频网站| 久久久久久久久免费视频了| 亚洲熟妇熟女久久| 人人妻人人看人人澡| 俄罗斯特黄特色一大片| 日本熟妇午夜| 久久亚洲真实| 中文亚洲av片在线观看爽| 香蕉国产在线看| 久久久久久人人人人人| 99久久国产精品久久久| 国产野战对白在线观看| 色综合欧美亚洲国产小说| 亚洲国产看品久久| 精品久久久久久久久久久久久 | 国产久久久一区二区三区| 欧美午夜高清在线| 高潮久久久久久久久久久不卡| 亚洲精品美女久久av网站| 窝窝影院91人妻| 宅男免费午夜| 欧洲精品卡2卡3卡4卡5卡区| 久久久久亚洲av毛片大全| 18美女黄网站色大片免费观看| 可以在线观看毛片的网站| 村上凉子中文字幕在线| 久久久久国内视频| 91九色精品人成在线观看| 黄片小视频在线播放| 成熟少妇高潮喷水视频| 18禁观看日本| 午夜福利视频1000在线观看| 亚洲av五月六月丁香网| 免费看美女性在线毛片视频| 啦啦啦韩国在线观看视频| 午夜a级毛片| 久久久国产成人精品二区| 日韩免费av在线播放| 国产亚洲av嫩草精品影院| 亚洲午夜精品一区,二区,三区| 久久中文字幕一级| 人妻丰满熟妇av一区二区三区| 久久99热这里只有精品18| 可以免费在线观看a视频的电影网站| 男人舔女人的私密视频| 午夜免费激情av| 欧美三级亚洲精品| 国产极品粉嫩免费观看在线| 波多野结衣av一区二区av| 巨乳人妻的诱惑在线观看| 欧美丝袜亚洲另类 | 午夜两性在线视频| 亚洲aⅴ乱码一区二区在线播放 | 熟女电影av网| 叶爱在线成人免费视频播放| 一区二区三区国产精品乱码| 免费高清视频大片| 国产av不卡久久| 亚洲精品国产一区二区精华液| 国产成人系列免费观看| 1024手机看黄色片| 午夜激情av网站| 精品久久久久久久毛片微露脸| 国产成人av激情在线播放| 久久99热这里只有精品18| 成人三级做爰电影| 亚洲国产精品sss在线观看| 日韩一卡2卡3卡4卡2021年| 日韩免费av在线播放| 后天国语完整版免费观看| 日韩中文字幕欧美一区二区| 精品久久久久久久末码| 欧美不卡视频在线免费观看 | 久久精品成人免费网站| 欧美不卡视频在线免费观看 | 99riav亚洲国产免费| 亚洲av五月六月丁香网| 久久精品国产亚洲av高清一级| 亚洲一区二区三区不卡视频| 日韩视频一区二区在线观看| 中文字幕高清在线视频| 女生性感内裤真人,穿戴方法视频| 看免费av毛片| АⅤ资源中文在线天堂| 日韩大尺度精品在线看网址| 国产精品,欧美在线| 观看免费一级毛片| 1024香蕉在线观看| 国产又黄又爽又无遮挡在线| 午夜福利一区二区在线看| www国产在线视频色| 2021天堂中文幕一二区在线观 | 91字幕亚洲| 精品熟女少妇八av免费久了| 丝袜在线中文字幕| 午夜精品在线福利| 成人特级黄色片久久久久久久| 午夜福利在线观看吧| 欧美精品亚洲一区二区| 大型av网站在线播放| 欧美日韩瑟瑟在线播放| 日本a在线网址| 757午夜福利合集在线观看| 搡老妇女老女人老熟妇| 欧美性长视频在线观看| 18禁美女被吸乳视频| 两性午夜刺激爽爽歪歪视频在线观看 | 午夜福利欧美成人| 亚洲av电影不卡..在线观看| 国产精品香港三级国产av潘金莲| 亚洲第一电影网av| 久久精品国产清高在天天线| 熟女少妇亚洲综合色aaa.| avwww免费| 国产精品久久电影中文字幕| 久久99热这里只有精品18| xxxwww97欧美| 精品国内亚洲2022精品成人| 亚洲激情在线av| 人妻丰满熟妇av一区二区三区| 国产男靠女视频免费网站| 精品欧美一区二区三区在线| 久久狼人影院| 国内久久婷婷六月综合欲色啪| 宅男免费午夜| 最近最新中文字幕大全电影3 | 天天添夜夜摸| 亚洲自偷自拍图片 自拍| 日日干狠狠操夜夜爽| 在线十欧美十亚洲十日本专区| 两个人视频免费观看高清| 亚洲av五月六月丁香网| 久久精品91无色码中文字幕| 一边摸一边做爽爽视频免费| 在线观看免费日韩欧美大片| 一区二区三区激情视频| 99riav亚洲国产免费| 精品少妇一区二区三区视频日本电影| 精品久久久久久久末码| 精品卡一卡二卡四卡免费| 嫩草影院精品99| 又黄又粗又硬又大视频| 最近最新中文字幕大全免费视频| 国产乱人伦免费视频| 亚洲五月天丁香| 一级a爱视频在线免费观看| 日韩高清综合在线| 身体一侧抽搐| 国产av又大| 好男人电影高清在线观看| 免费在线观看成人毛片| 久久久久亚洲av毛片大全| 国产单亲对白刺激| 精品国产一区二区三区四区第35| 午夜久久久在线观看| 男女做爰动态图高潮gif福利片| 亚洲天堂国产精品一区在线| 岛国视频午夜一区免费看| 免费观看精品视频网站| videosex国产| 欧美绝顶高潮抽搐喷水| 黄频高清免费视频| 欧美激情久久久久久爽电影| a在线观看视频网站| 亚洲国产欧美日韩在线播放| 三级毛片av免费| 91成人精品电影| 亚洲av电影不卡..在线观看| 欧美午夜高清在线| 99精品在免费线老司机午夜| 听说在线观看完整版免费高清| www日本在线高清视频| 久久精品成人免费网站| 老司机靠b影院| 三级毛片av免费| 美女大奶头视频| 亚洲第一电影网av| 一二三四社区在线视频社区8| 看片在线看免费视频| 少妇裸体淫交视频免费看高清 | 亚洲中文日韩欧美视频| 国内毛片毛片毛片毛片毛片| 老熟妇乱子伦视频在线观看| 高清毛片免费观看视频网站| 午夜福利高清视频| 99在线视频只有这里精品首页| 男女那种视频在线观看| 久久久国产精品麻豆| 国产乱人伦免费视频| 99精品久久久久人妻精品| 啦啦啦免费观看视频1| 国产人伦9x9x在线观看| 亚洲成人国产一区在线观看| 高清在线国产一区| 深夜精品福利| 黑人操中国人逼视频| 色精品久久人妻99蜜桃| 欧美日韩黄片免| 久久久久久亚洲精品国产蜜桃av| 久久99热这里只有精品18| 久久久久久久久免费视频了| 国产一区在线观看成人免费| 免费观看精品视频网站| 满18在线观看网站| 亚洲一区中文字幕在线| 中亚洲国语对白在线视频| 亚洲精品中文字幕在线视频| 免费在线观看亚洲国产| 搡老岳熟女国产| 正在播放国产对白刺激| 国产欧美日韩精品亚洲av| 搡老熟女国产l中国老女人| 国产精品 国内视频| 久久婷婷成人综合色麻豆| 女警被强在线播放| 久久精品aⅴ一区二区三区四区| 国产片内射在线|