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

    面向微電網(wǎng)的WSN分簇算法

    2018-03-16 06:17:09劉振興謝祥中

    施 義,劉振興,蔡 彬,謝祥中

    (1.武漢科技大學(xué) 信息科學(xué)與工程學(xué)院,湖北 武漢 430081;2.倫敦大學(xué)學(xué)院 UCL動(dòng)力工程系,倫敦 英國(guó) WC1E 7JE)

    0 引 言

    作為智能電網(wǎng)領(lǐng)域的重要組成部分,微電網(wǎng)中存在各類(lèi)分布式電源,網(wǎng)絡(luò)中使用大量的監(jiān)測(cè)設(shè)備來(lái)滿(mǎn)足其雙向、實(shí)時(shí)、高效的通信要求[1,2]。針對(duì)微電網(wǎng)就地控制層中采集節(jié)點(diǎn)多、節(jié)點(diǎn)位置分散等復(fù)雜的網(wǎng)絡(luò)情況[3,4],利用無(wú)線(xiàn)傳感器網(wǎng)絡(luò)(wireless sensor networks,WSN)技術(shù)來(lái)構(gòu)建該層的數(shù)據(jù)監(jiān)控網(wǎng)絡(luò)是保證電網(wǎng)安全穩(wěn)定運(yùn)行的重要手段。

    考慮微電網(wǎng)的實(shí)際特點(diǎn),監(jiān)控網(wǎng)絡(luò)選擇具有能量效率高、可擴(kuò)展性強(qiáng)的分簇結(jié)構(gòu)[5]。傳統(tǒng)的分簇算法如LEACH算法采用均勻分簇,由于簇間路由是單跳的方式,存在遠(yuǎn)離匯聚點(diǎn)的簇頭節(jié)點(diǎn)因通信代價(jià)高而過(guò)早死亡的問(wèn)題[6]。EEUC算法引入非均勻分簇概念,在簇頭選擇階段,其通信代價(jià)大,簇頭選擇和成簇等階段缺乏節(jié)點(diǎn)剩余能量的考慮,導(dǎo)致靠近匯聚點(diǎn)的簇頭節(jié)點(diǎn)沒(méi)有足夠多的能量完成大量數(shù)據(jù)的轉(zhuǎn)發(fā)而提前死亡[7,8],影響了網(wǎng)絡(luò)整體的生命周期。

    針對(duì)LEACH算法和EEUC算法的不足,結(jié)合微電網(wǎng)中分布式電源的分布情況,本文提出了一種基于節(jié)點(diǎn)實(shí)時(shí)能量的無(wú)線(xiàn)傳感器網(wǎng)絡(luò)非均勻分簇算法(node real-time energy based uneven clustering algorithm for WSN in micro-grid,NREUC)。NREUC算法在簇頭競(jìng)爭(zhēng)半徑中引入節(jié)點(diǎn)實(shí)時(shí)能量作為部分權(quán)重,在成簇階段,普通節(jié)點(diǎn)依據(jù)成簇判斷因子來(lái)選擇要加入的簇,達(dá)到非均勻分簇的目的。實(shí)驗(yàn)結(jié)果表明,與LEACH算法和EEUC算法相比,NREUC算法顯著延長(zhǎng)了網(wǎng)絡(luò)中第一個(gè)節(jié)點(diǎn)死亡時(shí)間,死亡節(jié)點(diǎn)的分布是均勻分布,避免了監(jiān)控盲區(qū)的出現(xiàn),提高了網(wǎng)絡(luò)的監(jiān)控質(zhì)量。

    1 微電網(wǎng)監(jiān)控網(wǎng)絡(luò)

    微電網(wǎng)的運(yùn)行控制中,分布式電源的監(jiān)控網(wǎng)絡(luò)如圖1所示。其中通信系統(tǒng)可以分為控制中心層、集中控制層和就地控制層[9]。終端電源單元將自身的狀態(tài)信息通過(guò)無(wú)線(xiàn)傳感器網(wǎng)絡(luò)發(fā)送給數(shù)據(jù)匯聚層的監(jiān)控單元,再由數(shù)據(jù)匯聚層統(tǒng)一通過(guò)光纖傳遞的方式傳遞給控制中心層進(jìn)行調(diào)控[10]。本文中,假設(shè)網(wǎng)絡(luò)里所有的監(jiān)控節(jié)點(diǎn)和區(qū)域匯聚點(diǎn)位置是固定的,每個(gè)監(jiān)控節(jié)點(diǎn)是同構(gòu)的,都可以成為簇頭、候選簇頭或普通節(jié)點(diǎn)。同時(shí),網(wǎng)絡(luò)中鏈路是對(duì)稱(chēng)的,已知對(duì)方發(fā)射功率時(shí),節(jié)點(diǎn)能夠根據(jù)接收信號(hào)強(qiáng)度計(jì)算距離[11]。

    圖1 微電網(wǎng)中分布式電源模塊監(jiān)控網(wǎng)絡(luò)

    根據(jù)圖1可得圖2監(jiān)控網(wǎng)絡(luò)的仿真場(chǎng)景,其中小圓圈是監(jiān)控節(jié)點(diǎn),五角星是匯聚點(diǎn)。微電網(wǎng)中分布式電源模塊監(jiān)控網(wǎng)絡(luò)中簇規(guī)模的整體趨勢(shì)是靠近匯聚節(jié)點(diǎn)的簇規(guī)模小數(shù)量多,遠(yuǎn)離匯聚節(jié)點(diǎn)的簇規(guī)模大數(shù)量小,這樣讓靠近匯聚節(jié)點(diǎn)的簇頭有更多的能量轉(zhuǎn)發(fā)消息,不至于過(guò)早死亡形成“熱區(qū)”。

    圖2 監(jiān)控網(wǎng)絡(luò)仿真場(chǎng)景

    2 分簇算法

    2.1 簇頭選舉

    無(wú)線(xiàn)傳感器網(wǎng)絡(luò)在進(jìn)行分簇時(shí),首先是簇頭選舉的過(guò)程。在簇頭選舉中,NREUC算法利用簇頭競(jìng)爭(zhēng)半徑來(lái)實(shí)現(xiàn)簇頭節(jié)點(diǎn)的非均勻分布,在競(jìng)爭(zhēng)半徑Ri計(jì)算公式中引入距離和能量?jī)蓚€(gè)因素

    (1)

    式中:dmax、dmin——網(wǎng)絡(luò)中節(jié)點(diǎn)到匯聚節(jié)點(diǎn)的最大與最小距離,d(si,s0)是節(jié)點(diǎn)i與匯聚節(jié)點(diǎn)的距離,R0是簇頭競(jìng)爭(zhēng)半徑的最大值,Ep是網(wǎng)絡(luò)中存活節(jié)點(diǎn)的平均能量與節(jié)點(diǎn)i的實(shí)時(shí)能量的比值,c是用于控制取值范圍的參數(shù)。

    考慮到節(jié)點(diǎn)的通信能耗,在NREUC算法的簇頭競(jìng)爭(zhēng)半徑計(jì)算中,節(jié)點(diǎn)與匯聚節(jié)點(diǎn)間距離仍是主要因素,距離匯聚節(jié)點(diǎn)越遠(yuǎn),競(jìng)爭(zhēng)半徑越大。節(jié)點(diǎn)自身實(shí)時(shí)能量是次要因素,對(duì)簇頭競(jìng)爭(zhēng)半徑起到微調(diào)的作用,能量大的節(jié)點(diǎn)的競(jìng)爭(zhēng)半徑偏大,但不會(huì)影響網(wǎng)絡(luò)中簇規(guī)模分布的整體結(jié)構(gòu)。

    2.2 成簇階段

    簇頭節(jié)點(diǎn)確定后,進(jìn)入成簇階段。在本文NREUC算法中,普通節(jié)點(diǎn)根據(jù)成簇判斷因子決定要加入的簇

    (2)

    式中:C(i,j,k)是第k輪,節(jié)點(diǎn)j到簇頭i的成簇判斷因子,Re(i,k)為簇頭i的實(shí)時(shí)能量,Ek是當(dāng)前網(wǎng)絡(luò)中節(jié)點(diǎn)的平均能量,d(si,sj)為節(jié)點(diǎn)j到簇頭i的距離。這樣在考慮就近成簇減少通信代價(jià)的同時(shí),讓剩余能量大的簇頭能夠接收更多的節(jié)點(diǎn),達(dá)到能量均衡的目的。

    2.3 NREUC算法

    圖3 NREUC算法

    初始化網(wǎng)絡(luò)時(shí),匯聚節(jié)點(diǎn)會(huì)向所有節(jié)點(diǎn)廣播消息,以便讓所有節(jié)點(diǎn)計(jì)算出與匯聚節(jié)點(diǎn)間的距離,用于后面的非均勻分簇以及節(jié)點(diǎn)向基站傳輸信息時(shí)發(fā)送功率的選擇。NREUC算法采用循環(huán)機(jī)制,每輪中包括簇頭的選舉及簇的形成、簇間路由建立和數(shù)據(jù)的傳輸。網(wǎng)絡(luò)中的節(jié)點(diǎn)有4種狀態(tài):普通節(jié)點(diǎn)狀態(tài),候選簇頭狀態(tài),簇頭狀態(tài)和死亡狀態(tài)。每個(gè)節(jié)點(diǎn)設(shè)有簇頭標(biāo)志位G,標(biāo)志位G初始態(tài)置0,在節(jié)點(diǎn)當(dāng)選簇頭后,標(biāo)志位G會(huì)置1。節(jié)點(diǎn)的狀態(tài)每輪會(huì)動(dòng)態(tài)變化,簇頭在候選簇頭中產(chǎn)生,其中候選簇頭節(jié)點(diǎn)的簇頭競(jìng)爭(zhēng)半徑為Ri,由式(1)得到。大小不同的競(jìng)爭(zhēng)半徑形成規(guī)模不等的簇。

    NREUC算法的具體實(shí)現(xiàn)步驟如下(如圖3所示):

    步驟1 標(biāo)志位G的初始化

    對(duì)所有節(jié)點(diǎn)的簇頭標(biāo)志位G置0。

    步驟2 候選簇頭選舉

    先判斷簇頭標(biāo)志位G,允許標(biāo)志位為0的節(jié)點(diǎn)進(jìn)行候選簇頭的選舉,否則節(jié)點(diǎn)休眠直至簇頭競(jìng)選結(jié)束。再按概率選出部分節(jié)點(diǎn)成為候選簇頭,其它落選節(jié)點(diǎn)休眠直至簇頭競(jìng)選結(jié)束。

    步驟3 簇頭選舉

    依據(jù)式(1)求出的簇頭競(jìng)爭(zhēng)半徑來(lái)實(shí)現(xiàn)簇頭的選舉。每輪簇頭競(jìng)選中,候選簇頭節(jié)點(diǎn)i有一張?jiān)摴?jié)點(diǎn)的鄰候選簇頭信息表,表內(nèi)的任意節(jié)點(diǎn)j都滿(mǎn)足節(jié)點(diǎn)j與節(jié)點(diǎn)i間的距離小于這兩個(gè)節(jié)點(diǎn)的競(jìng)爭(zhēng)半徑的最大值的條件。實(shí)時(shí)能量最大的候選簇頭最先被選出擔(dān)任簇頭,同時(shí)它的鄰候選簇頭信息表中的候選簇頭都退選成為普通節(jié)點(diǎn)。然后再在剩余的候選簇頭中,選出實(shí)時(shí)能量最大的節(jié)點(diǎn)成為簇頭,它的鄰候選簇頭信息表中的候選簇頭也都退選成為普通節(jié)點(diǎn),依次循環(huán)至沒(méi)有候選簇頭。在節(jié)點(diǎn)出任簇頭后,該節(jié)點(diǎn)標(biāo)志位G置1。

    步驟4 成簇階段

    網(wǎng)絡(luò)中的普通節(jié)點(diǎn)依據(jù)成簇判斷因子的式(2)來(lái)選擇簇頭,形成規(guī)模大小不等的簇。

    步驟5 信息傳輸階段

    本網(wǎng)絡(luò)中采用簇內(nèi)單跳,簇首間多跳的方式進(jìn)行信息的傳輸。

    步驟6 簇頭標(biāo)志G的再判斷

    信息傳輸完成后,計(jì)算各節(jié)點(diǎn)在本輪結(jié)束時(shí)的剩余能量,將剩余能量大于或等于平均剩余能量的簇頭節(jié)點(diǎn)的標(biāo)志位G置0,讓其可以繼續(xù)參與到后面的簇頭競(jìng)選中。

    步驟7 死亡節(jié)點(diǎn)個(gè)數(shù)判斷

    若網(wǎng)絡(luò)中的節(jié)點(diǎn)死亡數(shù)超過(guò)預(yù)設(shè)最大值,則網(wǎng)絡(luò)宣告死亡,否則進(jìn)入新一輪的網(wǎng)絡(luò)分簇、信息傳遞中。

    3 算法仿真

    3.1 仿真參數(shù)設(shè)置

    本文在MATLAB中分別編碼LEACH算法、EEUC算法和NREUC算法,并進(jìn)行仿真性能分析。實(shí)驗(yàn)中,在邊長(zhǎng)200 m的方形區(qū)域內(nèi)隨機(jī)分布400個(gè)節(jié)點(diǎn),網(wǎng)絡(luò)中消耗的能量主要由發(fā)射電路能耗和功率放大電路的能耗組成。當(dāng)節(jié)點(diǎn)發(fā)射比特?cái)?shù)為k比特的數(shù)據(jù)到相聚d距離的位置時(shí),通信能耗計(jì)算為

    (3)

    式中:Eelec——發(fā)射電路消耗的能量,εfs、εmp——兩類(lèi)信道模型中功率放大電路的最大消耗能量。相關(guān)仿真場(chǎng)景參數(shù)見(jiàn)表1。

    表1 仿真場(chǎng)景參數(shù)

    在式(1)簇頭競(jìng)爭(zhēng)半徑Ri的計(jì)算中,系數(shù)c和最大競(jìng)爭(zhēng)半徑R0的值是需要預(yù)先設(shè)定的。其中,為方便與EEUC算法比較,系數(shù)c的值沿用原EEUC的值0.5。對(duì)于候選簇頭的最大競(jìng)爭(zhēng)半徑R0,它影響最后簇的生成個(gè)數(shù),最大競(jìng)爭(zhēng)半徑大時(shí),則生成的簇個(gè)數(shù)少,從而影響最終網(wǎng)絡(luò)中簇分布情況。這里采用實(shí)驗(yàn)方法得出本仿真環(huán)境下最大競(jìng)爭(zhēng)半徑的最適值為90 m。

    3.2 仿真結(jié)果分析

    網(wǎng)絡(luò)仿真中,本文主要在網(wǎng)絡(luò)生命周期及網(wǎng)絡(luò)的死亡節(jié)點(diǎn)分布上進(jìn)行比較分析。網(wǎng)絡(luò)生命周期越長(zhǎng),特別是網(wǎng)絡(luò)中第一個(gè)節(jié)點(diǎn)的死亡時(shí)間越長(zhǎng),則表明網(wǎng)絡(luò)中能量越均衡,節(jié)點(diǎn)不會(huì)過(guò)早死亡。網(wǎng)絡(luò)的死亡節(jié)點(diǎn)分布影響著整個(gè)網(wǎng)絡(luò)的監(jiān)控質(zhì)量,死亡節(jié)點(diǎn)分布越均勻,監(jiān)控質(zhì)量越好。

    3.2.1 網(wǎng)絡(luò)生命周期

    針對(duì)采集信息精確度要求較高的網(wǎng)絡(luò),節(jié)點(diǎn)的初始死亡時(shí)間決定著網(wǎng)絡(luò)整體的生命周期。如微電網(wǎng)中,單元級(jí)別的設(shè)備狀態(tài)數(shù)據(jù)采集需要達(dá)到全覆蓋和高精度的要求,盲區(qū)的出現(xiàn)會(huì)影響整個(gè)網(wǎng)絡(luò)的運(yùn)行安全。

    圖4中比較了3種分簇算法的網(wǎng)絡(luò)生命周期,NREUC算法相比于LEACH算法和EEUC算法,網(wǎng)絡(luò)中第一個(gè)節(jié)點(diǎn)的死亡時(shí)間得到了很大的延長(zhǎng),網(wǎng)絡(luò)整體的存活時(shí)間也得到了增加。NREUC算法中,簇間多跳路由的方式解決了LEACH算法中遠(yuǎn)離匯聚點(diǎn)的簇頭節(jié)點(diǎn)因通信消耗大而過(guò)早死亡的問(wèn)題。在候選簇頭的選擇上,NREUC算法加入簇頭標(biāo)志位的判斷,讓上輪中擔(dān)任過(guò)簇頭的這類(lèi)能量消耗大的節(jié)點(diǎn)不參與之后的簇頭競(jìng)選,提高候選簇頭質(zhì)量,避免不必要的能量消耗,緩解了EEUC算法中簇頭選擇能量開(kāi)銷(xiāo)大的問(wèn)題。在NREUC算法簇頭競(jìng)爭(zhēng)半徑計(jì)算中,考慮節(jié)點(diǎn)自身能量和其與匯聚節(jié)點(diǎn)間距離兩個(gè)因素,使分簇的規(guī)模較EEUC算法更加合理。根據(jù)圖4中NREUC算法曲線(xiàn)的衰減走勢(shì)可以看出,NREUC算法在同等的能耗條件下提高了網(wǎng)絡(luò)中的能量均衡程度,在網(wǎng)絡(luò)生存周期后段,每輪中的存活節(jié)點(diǎn)占有率都要比LEACH算法和EEUC算法高。

    圖4 3種分簇算法的網(wǎng)絡(luò)生命周期

    3.2.2 網(wǎng)路死亡節(jié)點(diǎn)分布

    根據(jù)圖5和圖6中網(wǎng)絡(luò)的死亡節(jié)點(diǎn)分布可以看出,圖5中EEUC算法中死亡節(jié)點(diǎn)大多集中在靠近匯聚節(jié)點(diǎn)的區(qū)域,圖6中NREUC算法中死亡節(jié)點(diǎn)分布均勻。

    圖5 EEUC算法的死亡節(jié)點(diǎn)分布

    圖6 NREUC算法的死亡節(jié)點(diǎn)分布

    NREUC算法在每輪的簇頭選舉階段都會(huì)考慮各個(gè)節(jié)點(diǎn)的剩余能量問(wèn)題,減少低能量節(jié)點(diǎn)成為簇頭的可能性,在簇頭競(jìng)爭(zhēng)半徑計(jì)算中加入簇頭自身能量的考慮因素,利用大小不等的簇頭競(jìng)爭(zhēng)半徑達(dá)到網(wǎng)絡(luò)中簇頭節(jié)點(diǎn)非均勻分布的目的。同時(shí)在成簇時(shí),NREUC算法將簇頭能量考慮進(jìn)去,普通節(jié)點(diǎn)在選擇簇頭時(shí)不再是單純的就近選擇,而是根據(jù)節(jié)點(diǎn)與簇頭間距離及簇頭剩余能量?jī)蓚€(gè)因素進(jìn)行判斷,最終使得網(wǎng)絡(luò)中最后的死亡節(jié)點(diǎn)分布是均勻的,解決了EEUC算法中靠近匯聚點(diǎn)的簇頭節(jié)點(diǎn)因轉(zhuǎn)發(fā)過(guò)多數(shù)據(jù)而提前死亡的問(wèn)題,達(dá)到提高網(wǎng)絡(luò)監(jiān)控質(zhì)量的目的。

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

    本文提出了一種面向微電網(wǎng)的WSN非均勻分簇算法,算法的核心思想是在簇頭選擇和成簇階段均加入節(jié)點(diǎn)的實(shí)時(shí)能量因素進(jìn)行考慮,以均衡網(wǎng)絡(luò)中的能量消耗,延長(zhǎng)網(wǎng)絡(luò)的生存時(shí)間。針對(duì)微電網(wǎng)監(jiān)控網(wǎng)絡(luò)中節(jié)點(diǎn)分布廣、監(jiān)控精度高的特點(diǎn),該算法相比于傳統(tǒng)的分簇算法,網(wǎng)絡(luò)中能量均衡的程度更高,節(jié)點(diǎn)初始死亡時(shí)間更晚,網(wǎng)絡(luò)生命期更長(zhǎng)。

    雖然本文的NREUC算法對(duì)無(wú)線(xiàn)傳感器網(wǎng)絡(luò)中能量均衡起到了較好的效果,但在分簇時(shí)需要一定的消息開(kāi)銷(xiāo),可能會(huì)影響網(wǎng)絡(luò)生存時(shí)間。下一步的工作考慮減少網(wǎng)絡(luò)中每輪的能量消耗,同時(shí)驗(yàn)證該算法在實(shí)際微電網(wǎng)環(huán)境中應(yīng)用的可行性。

    [1]YANG Xinfa,SU Jian,LYU Zhipeng,et al.Overview on micro-grid technology[J].Proceedings of the Chinese Society of Electrical Engineering,2014,34(1):57-70(in Chinese).[楊新法,蘇劍,呂志鵬,等.微電網(wǎng)技術(shù)綜述[J].中國(guó)電機(jī)工程學(xué)報(bào),2014,34(1):57-70.]

    [2]Arefifar SA,Mohamed ARI,El-Fouly T.Optimized multiple microgrid-based clustering of active distribution systems consi-dering communication and control requirements[J].IEEE Transactions on Industrial Electronics,2015,62(2):711-723.[3]WU Xiong,WANG Xiuli,LIU Shimin,et al.Summary of research on microgrid energy management system[J].Electric Power Automation Equipment,2014,34(10):7-14(in Chinese).[吳雄,王秀麗,劉世民,等.微電網(wǎng)能量管理系統(tǒng)研究綜述[J].電力自動(dòng)化設(shè)備,2014,34(10):7-14.]

    [4]Angelo PM,Ted S.A survey of techniques used to control microgrid generation and storage during island operation[J].IEEE Transactions on Power Systems,2015,30(2):701-709.

    [5]Bin LI,Wang WJ,Yin QY,et al.An energy-efficient geographic routing based on cooperative transmission in wireless sensor networks[J].Science China Information Sciences,2013,56(7):1-10.

    [6]FENG Lin,RAN Xiaomin,MEI Guanlin.WSN coverage tactics based on balanced node energy consumption[J].Computer Engineering and Design,2015,36(9):2334-2339(in Chinese).[馮琳,冉曉旻,梅關(guān)林.基于節(jié)點(diǎn)能耗均衡的WSN覆蓋策略[J].計(jì)算機(jī)工程與設(shè)計(jì),2015,36(9):2334-2339.]

    [7]YUE Liying,DAI Yueming,WU Dinghui.Energy optimized uneven clustering routing protocol in wireless sensor networks[J].Computer Engineering and Applications,2015,51(15):80-85(in Chinese).[岳麗穎,戴月明,吳定會(huì).一種能量?jī)?yōu)化WSNs非均勻分簇路由協(xié)議[J].計(jì)算機(jī)工程與應(yīng)用,2015,51(15):80-85.]

    [8]WEN Peizhi,XU Chenjiao,DENG Zhenrong,et al.Cluster-routing protocol for heterogenous multi-level wireless sensor networks[J].Computer Engineering and Design,2016,37(6):1471-1477(in Chinese).[溫佩芝,許晨蛟,鄧珍榮,等.多級(jí)異構(gòu)無(wú)線(xiàn)傳感器網(wǎng)絡(luò)分簇路由協(xié)議[J].計(jì)算機(jī)工程與設(shè)計(jì),2016,37(6):1471-1477.]

    [9]ZHANG Xinchang.Solution and application of micro-grid ope-ration control[J].Power System Protection and Control,2014,42(10):141-146(in Chinese).[張新昌.微電網(wǎng)運(yùn)行控制解決方案及應(yīng)用[J].電力系統(tǒng)保護(hù)與控制,2014,42(10):141-146.]

    [10]Tan X,Li Q,Wang H.Advances and trends of energy sto-rage technology in microgrid[J].International Journal of Electrical Power & Energy Systems,2013,44(1):179-191.

    [11]JIANG Changjiang,SHI Weiren,TANG Xianlun,et al.Energy-balanced unequal clustering routing protocol for wireless sensor networks[J].Journal of Software,2012,34(5):1222-1232(in Chinese).[蔣暢江,石為人,唐賢倫,等.能量均衡的無(wú)線(xiàn)傳感器網(wǎng)絡(luò)非均勻分簇路由協(xié)議[J].軟件學(xué)報(bào),2012,34(5):1222-1232.]

    最黄视频免费看| 看免费成人av毛片| 午夜福利视频精品| av免费观看日本| 免费大片黄手机在线观看| 国产乱人偷精品视频| 国产精品久久久久成人av| 大片电影免费在线观看免费| 国产黄频视频在线观看| 不卡视频在线观看欧美| 久热久热在线精品观看| 精品国产国语对白av| 女人精品久久久久毛片| 高清黄色对白视频在线免费看| 尾随美女入室| 男人舔女人的私密视频| av电影中文网址| 午夜久久久在线观看| 黑人欧美特级aaaaaa片| 国产成人a∨麻豆精品| 人人妻人人添人人爽欧美一区卜| 多毛熟女@视频| 久久影院123| 久久人人爽av亚洲精品天堂| 18禁裸乳无遮挡动漫免费视频| 99热国产这里只有精品6| 亚洲精品乱码久久久久久按摩| 欧美精品一区二区大全| 亚洲欧美精品自产自拍| 你懂的网址亚洲精品在线观看| 欧美xxⅹ黑人| 国产伦理片在线播放av一区| 亚洲国产毛片av蜜桃av| 蜜桃国产av成人99| 国产精品国产三级国产专区5o| a 毛片基地| 热99国产精品久久久久久7| 免费观看性生交大片5| 免费大片18禁| 搡老乐熟女国产| 制服人妻中文乱码| 欧美精品国产亚洲| 日本色播在线视频| 精品福利永久在线观看| 成年美女黄网站色视频大全免费| 伊人亚洲综合成人网| 亚洲精品久久久久久婷婷小说| 性高湖久久久久久久久免费观看| 亚洲成人一二三区av| 97人妻天天添夜夜摸| 51国产日韩欧美| 欧美日韩av久久| 侵犯人妻中文字幕一二三四区| 丝袜脚勾引网站| 人妻人人澡人人爽人人| 少妇猛男粗大的猛烈进出视频| 在线观看美女被高潮喷水网站| 午夜久久久在线观看| 有码 亚洲区| 亚洲内射少妇av| 少妇的丰满在线观看| 亚洲国产av新网站| 成人国产麻豆网| av福利片在线| 看非洲黑人一级黄片| 观看av在线不卡| 国产精品一国产av| 欧美亚洲 丝袜 人妻 在线| 视频在线观看一区二区三区| 男女免费视频国产| 日韩熟女老妇一区二区性免费视频| 免费观看av网站的网址| 成人亚洲精品一区在线观看| 国产成人精品一,二区| 精品亚洲成a人片在线观看| 婷婷成人精品国产| 国产高清三级在线| 免费看不卡的av| 黑丝袜美女国产一区| 蜜臀久久99精品久久宅男| freevideosex欧美| 男女啪啪激烈高潮av片| a级毛片在线看网站| 精品熟女少妇av免费看| 午夜福利网站1000一区二区三区| 制服诱惑二区| 天天躁夜夜躁狠狠久久av| 最近最新中文字幕大全免费视频 | 日韩在线高清观看一区二区三区| 日韩,欧美,国产一区二区三区| 欧美日本中文国产一区发布| 少妇的逼水好多| 青春草视频在线免费观看| 男人添女人高潮全过程视频| 精品亚洲乱码少妇综合久久| 欧美丝袜亚洲另类| 999精品在线视频| 美女视频免费永久观看网站| 久久久久久久国产电影| 国产精品国产三级国产专区5o| 国内精品宾馆在线| 日本欧美国产在线视频| 国产免费一区二区三区四区乱码| 日日爽夜夜爽网站| 国产激情久久老熟女| 免费播放大片免费观看视频在线观看| freevideosex欧美| 91午夜精品亚洲一区二区三区| 人妻一区二区av| 亚洲丝袜综合中文字幕| 99热全是精品| 黄色视频在线播放观看不卡| 国产精品嫩草影院av在线观看| 毛片一级片免费看久久久久| 日本av手机在线免费观看| 观看美女的网站| 精品少妇内射三级| 精品久久久精品久久久| 国语对白做爰xxxⅹ性视频网站| 毛片一级片免费看久久久久| 啦啦啦视频在线资源免费观看| 一区二区三区四区激情视频| 亚洲欧美日韩卡通动漫| 天堂8中文在线网| 黑人欧美特级aaaaaa片| 丝袜美足系列| 国产精品秋霞免费鲁丝片| 最新的欧美精品一区二区| 亚洲人成77777在线视频| 色5月婷婷丁香| 国产又色又爽无遮挡免| 美女中出高潮动态图| 国产免费又黄又爽又色| 晚上一个人看的免费电影| 99热这里只有是精品在线观看| 在线观看www视频免费| 韩国av在线不卡| 国产黄色视频一区二区在线观看| av不卡在线播放| 日韩欧美精品免费久久| 大片免费播放器 马上看| 午夜老司机福利剧场| 看免费av毛片| 综合色丁香网| av天堂久久9| 91久久精品国产一区二区三区| 精品一区在线观看国产| 中文字幕av电影在线播放| 好男人视频免费观看在线| 最近手机中文字幕大全| 日韩欧美精品免费久久| 国产精品久久久久成人av| av女优亚洲男人天堂| 欧美 亚洲 国产 日韩一| 五月天丁香电影| 婷婷成人精品国产| 十八禁高潮呻吟视频| 国产精品国产av在线观看| 日韩在线高清观看一区二区三区| 亚洲丝袜综合中文字幕| 亚洲欧美成人精品一区二区| 亚洲久久久国产精品| 99re6热这里在线精品视频| 在线精品无人区一区二区三| 亚洲在久久综合| 人人妻人人澡人人爽人人夜夜| 亚洲第一区二区三区不卡| 大话2 男鬼变身卡| 国产黄频视频在线观看| 波多野结衣一区麻豆| 精品国产一区二区久久| 男女边吃奶边做爰视频| 人人妻人人添人人爽欧美一区卜| 中文字幕av电影在线播放| 日韩欧美精品免费久久| 欧美精品亚洲一区二区| 国产成人欧美| 国产精品久久久av美女十八| 熟女av电影| 90打野战视频偷拍视频| 亚洲精品久久久久久婷婷小说| 国产精品秋霞免费鲁丝片| 26uuu在线亚洲综合色| 成年美女黄网站色视频大全免费| 22中文网久久字幕| 国产极品天堂在线| 如日韩欧美国产精品一区二区三区| 久久这里只有精品19| 男女边吃奶边做爰视频| 两个人看的免费小视频| 国产在视频线精品| 国产亚洲一区二区精品| 久久久a久久爽久久v久久| 97超碰精品成人国产| 少妇高潮的动态图| 天天操日日干夜夜撸| 亚洲熟女精品中文字幕| 色吧在线观看| 一边亲一边摸免费视频| 国产1区2区3区精品| 亚洲国产色片| 九九爱精品视频在线观看| 少妇精品久久久久久久| h视频一区二区三区| 人人妻人人爽人人添夜夜欢视频| 青春草视频在线免费观看| 国产不卡av网站在线观看| 久久久久久人妻| 一区二区日韩欧美中文字幕 | 国产日韩欧美在线精品| 又黄又爽又刺激的免费视频.| 久久鲁丝午夜福利片| 国产精品久久久久久久电影| 国产精品国产三级专区第一集| 日产精品乱码卡一卡2卡三| 日韩精品有码人妻一区| 亚洲一级一片aⅴ在线观看| av又黄又爽大尺度在线免费看| 免费大片18禁| 欧美日韩av久久| 乱人伦中国视频| 高清视频免费观看一区二区| 两个人免费观看高清视频| 精品国产一区二区三区久久久樱花| 满18在线观看网站| 欧美亚洲 丝袜 人妻 在线| 久久久久视频综合| 精品卡一卡二卡四卡免费| 大香蕉97超碰在线| 女的被弄到高潮叫床怎么办| 国产深夜福利视频在线观看| 亚洲国产看品久久| 日本黄色日本黄色录像| 亚洲av欧美aⅴ国产| 99香蕉大伊视频| 激情五月婷婷亚洲| 校园人妻丝袜中文字幕| 国产在视频线精品| 性色av一级| 又黄又爽又刺激的免费视频.| 亚洲丝袜综合中文字幕| 人人妻人人澡人人爽人人夜夜| 午夜日本视频在线| 一级毛片 在线播放| 9色porny在线观看| 亚洲av免费高清在线观看| 久久久久久久大尺度免费视频| 考比视频在线观看| 天堂中文最新版在线下载| 国产一区二区三区av在线| 精品午夜福利在线看| 最近中文字幕高清免费大全6| 美女国产高潮福利片在线看| 精品99又大又爽又粗少妇毛片| 在线观看人妻少妇| 亚洲精品美女久久av网站| 国产黄频视频在线观看| 亚洲内射少妇av| 最后的刺客免费高清国语| 国产爽快片一区二区三区| 蜜桃在线观看..| 永久网站在线| 制服诱惑二区| 久久ye,这里只有精品| 少妇人妻精品综合一区二区| 欧美97在线视频| 亚洲中文av在线| av免费观看日本| 久久久精品94久久精品| 久久韩国三级中文字幕| 老司机亚洲免费影院| 飞空精品影院首页| 国产又爽黄色视频| 成人亚洲欧美一区二区av| 中文欧美无线码| 亚洲av.av天堂| 久久久久人妻精品一区果冻| 成人无遮挡网站| 中文乱码字字幕精品一区二区三区| 国产精品蜜桃在线观看| 夜夜爽夜夜爽视频| 十八禁网站网址无遮挡| 香蕉丝袜av| 亚洲国产av影院在线观看| 国产精品一区二区在线不卡| 日韩免费高清中文字幕av| 如何舔出高潮| 热99久久久久精品小说推荐| 久久人妻熟女aⅴ| 国产精品秋霞免费鲁丝片| 黄片播放在线免费| 国产黄色视频一区二区在线观看| 乱码一卡2卡4卡精品| 侵犯人妻中文字幕一二三四区| 久久精品久久久久久久性| 久久精品熟女亚洲av麻豆精品| 亚洲精品日韩在线中文字幕| xxx大片免费视频| 一二三四中文在线观看免费高清| 九色亚洲精品在线播放| 国产色爽女视频免费观看| 欧美日韩综合久久久久久| 搡女人真爽免费视频火全软件| 在线观看美女被高潮喷水网站| 人妻系列 视频| 另类亚洲欧美激情| 22中文网久久字幕| 亚洲国产欧美在线一区| 国产亚洲精品久久久com| 美女视频免费永久观看网站| 热99久久久久精品小说推荐| 99视频精品全部免费 在线| 制服诱惑二区| 久久久久久久精品精品| 国产无遮挡羞羞视频在线观看| 午夜免费男女啪啪视频观看| 国产精品国产三级国产av玫瑰| 欧美+日韩+精品| 交换朋友夫妻互换小说| 日本色播在线视频| 久久人人爽av亚洲精品天堂| 2021少妇久久久久久久久久久| 久久久a久久爽久久v久久| 久久久久久久久久人人人人人人| 欧美bdsm另类| 成人影院久久| 亚洲欧美日韩另类电影网站| 免费看不卡的av| 这个男人来自地球电影免费观看 | 亚洲欧美日韩另类电影网站| 少妇的逼水好多| 日本av免费视频播放| 黄色毛片三级朝国网站| 三上悠亚av全集在线观看| 成年人免费黄色播放视频| 男人操女人黄网站| 久久精品久久精品一区二区三区| 中文字幕亚洲精品专区| 成人免费观看视频高清| 久久精品久久久久久噜噜老黄| 亚洲 欧美一区二区三区| 热re99久久精品国产66热6| 99国产精品免费福利视频| 99香蕉大伊视频| 观看美女的网站| 国产爽快片一区二区三区| 亚洲成色77777| 一本色道久久久久久精品综合| 国产亚洲精品久久久com| av黄色大香蕉| 久久人妻熟女aⅴ| 黑人高潮一二区| 高清欧美精品videossex| 国产色爽女视频免费观看| 永久网站在线| 91久久精品国产一区二区三区| 黄网站色视频无遮挡免费观看| 国产精品国产av在线观看| av在线观看视频网站免费| www.熟女人妻精品国产 | 国产一区二区三区综合在线观看 | 国产精品久久久久成人av| 夜夜爽夜夜爽视频| 插逼视频在线观看| av视频免费观看在线观看| 午夜激情久久久久久久| 最近2019中文字幕mv第一页| 久久久久久久大尺度免费视频| 国产在线一区二区三区精| 人人妻人人爽人人添夜夜欢视频| 久久99一区二区三区| xxx大片免费视频| 亚洲欧美清纯卡通| 午夜精品国产一区二区电影| 久久精品aⅴ一区二区三区四区 | 亚洲成人av在线免费| 人人妻人人澡人人看| 人人妻人人添人人爽欧美一区卜| 国产精品久久久久久精品古装| 国产精品国产av在线观看| 日本黄大片高清| xxx大片免费视频| 久久精品夜色国产| 午夜视频国产福利| 一级片免费观看大全| 国产成人精品婷婷| 97人妻天天添夜夜摸| 成年动漫av网址| 亚洲av日韩在线播放| 有码 亚洲区| 精品一区二区三区视频在线| 日日啪夜夜爽| 日韩熟女老妇一区二区性免费视频| 99热国产这里只有精品6| 91精品国产国语对白视频| 热re99久久精品国产66热6| 伦精品一区二区三区| 如日韩欧美国产精品一区二区三区| av在线播放精品| 国产精品久久久久久精品古装| 99国产精品免费福利视频| 中文欧美无线码| 一二三四中文在线观看免费高清| 自线自在国产av| 免费av不卡在线播放| 最近最新中文字幕免费大全7| 亚洲精品国产av成人精品| 精品一品国产午夜福利视频| 中文字幕另类日韩欧美亚洲嫩草| av在线观看视频网站免费| 777米奇影视久久| 成人漫画全彩无遮挡| 飞空精品影院首页| 亚洲av电影在线进入| 在线观看三级黄色| 99热全是精品| 在线天堂中文资源库| av片东京热男人的天堂| 99久久精品国产国产毛片| 午夜福利,免费看| 日本vs欧美在线观看视频| 国产黄色视频一区二区在线观看| 亚洲国产日韩一区二区| 自拍欧美九色日韩亚洲蝌蚪91| 日韩一区二区三区影片| 国产精品.久久久| 欧美少妇被猛烈插入视频| 18禁动态无遮挡网站| 肉色欧美久久久久久久蜜桃| 国产激情久久老熟女| 下体分泌物呈黄色| 国产探花极品一区二区| 国产精品人妻久久久影院| 蜜桃国产av成人99| 成人毛片60女人毛片免费| 亚洲 欧美一区二区三区| 亚洲欧美色中文字幕在线| 一区二区三区精品91| 嫩草影院入口| 日韩在线高清观看一区二区三区| 亚洲欧美成人综合另类久久久| 国产精品熟女久久久久浪| 国产在线视频一区二区| 大码成人一级视频| 秋霞在线观看毛片| 美女脱内裤让男人舔精品视频| 丰满乱子伦码专区| 亚洲伊人色综图| 亚洲精品乱码久久久久久按摩| 一二三四中文在线观看免费高清| 亚洲av日韩在线播放| 我的女老师完整版在线观看| 麻豆乱淫一区二区| 美女脱内裤让男人舔精品视频| 在线观看国产h片| 伊人久久国产一区二区| 欧美亚洲 丝袜 人妻 在线| 一本久久精品| 99国产精品免费福利视频| 黑人高潮一二区| 啦啦啦啦在线视频资源| 日本-黄色视频高清免费观看| 最近中文字幕高清免费大全6| 成年女人在线观看亚洲视频| 亚洲av成人精品一二三区| 免费黄色在线免费观看| 日本免费在线观看一区| 国产综合精华液| 日本欧美视频一区| 视频区图区小说| 欧美丝袜亚洲另类| 高清黄色对白视频在线免费看| 久久 成人 亚洲| 18+在线观看网站| 美女国产视频在线观看| 麻豆乱淫一区二区| 狠狠婷婷综合久久久久久88av| 两个人看的免费小视频| 国产日韩欧美在线精品| 69精品国产乱码久久久| 水蜜桃什么品种好| 成人亚洲欧美一区二区av| 男女国产视频网站| 97在线人人人人妻| 亚洲熟女精品中文字幕| 国产精品99久久99久久久不卡 | 亚洲综合精品二区| 国产伦理片在线播放av一区| 亚洲精品色激情综合| 亚洲国产精品国产精品| 午夜av观看不卡| 最近手机中文字幕大全| 天天躁夜夜躁狠狠躁躁| 国产日韩欧美视频二区| 日日啪夜夜爽| 精品一区二区三卡| 久久久久久久亚洲中文字幕| 精品99又大又爽又粗少妇毛片| 美女视频免费永久观看网站| 下体分泌物呈黄色| 成年人午夜在线观看视频| 午夜福利影视在线免费观看| 欧美+日韩+精品| 激情五月婷婷亚洲| av黄色大香蕉| 看非洲黑人一级黄片| 亚洲五月色婷婷综合| 国产有黄有色有爽视频| 91精品国产国语对白视频| 日韩成人伦理影院| 各种免费的搞黄视频| 纯流量卡能插随身wifi吗| 亚洲一区二区三区欧美精品| 亚洲国产最新在线播放| 久久女婷五月综合色啪小说| 亚洲欧洲精品一区二区精品久久久 | 天美传媒精品一区二区| 男人添女人高潮全过程视频| 日韩制服骚丝袜av| 晚上一个人看的免费电影| 免费高清在线观看视频在线观看| 我要看黄色一级片免费的| 国产在线视频一区二区| 亚洲精品国产色婷婷电影| av在线老鸭窝| 草草在线视频免费看| 免费黄色在线免费观看| 美女大奶头黄色视频| 美国免费a级毛片| videosex国产| 综合色丁香网| 午夜免费男女啪啪视频观看| 亚洲综合色网址| 日本欧美视频一区| 日日撸夜夜添| 黑人巨大精品欧美一区二区蜜桃 | 9热在线视频观看99| 国产成人精品福利久久| 黄色怎么调成土黄色| 日韩成人伦理影院| 国国产精品蜜臀av免费| 免费少妇av软件| 一区二区av电影网| 三级国产精品片| 亚洲精品日韩在线中文字幕| 国产精品嫩草影院av在线观看| 九九在线视频观看精品| 91精品伊人久久大香线蕉| 日韩,欧美,国产一区二区三区| 国产日韩一区二区三区精品不卡| 午夜视频国产福利| 五月玫瑰六月丁香| 国产一区有黄有色的免费视频| 日日撸夜夜添| 国产日韩欧美亚洲二区| 亚洲欧美成人综合另类久久久| 国产一级毛片在线| 在线观看免费日韩欧美大片| 人妻系列 视频| a级毛片在线看网站| 久久久久精品性色| 在线看a的网站| 五月天丁香电影| 国产福利在线免费观看视频| 国产精品久久久久成人av| 精品国产露脸久久av麻豆| 丝袜美足系列| 免费高清在线观看日韩| 久久韩国三级中文字幕| 女性生殖器流出的白浆| 国产1区2区3区精品| 国产一区二区激情短视频 | 国产有黄有色有爽视频| 丰满少妇做爰视频| 欧美最新免费一区二区三区| 国产精品久久久久久av不卡| 国产一区二区激情短视频 | 亚洲经典国产精华液单| 午夜福利视频精品| 女人久久www免费人成看片| 一二三四中文在线观看免费高清| 91精品三级在线观看| 久久综合国产亚洲精品| 亚洲欧美成人精品一区二区| 97在线视频观看| 欧美变态另类bdsm刘玥| www.熟女人妻精品国产 | 色吧在线观看| 亚洲精品一区蜜桃| 精品一区二区免费观看| 久久精品国产a三级三级三级| 七月丁香在线播放| 韩国精品一区二区三区 | 少妇 在线观看| freevideosex欧美| 亚洲精品第二区| 亚洲精品,欧美精品| 婷婷成人精品国产| 18在线观看网站| 国产精品国产三级国产专区5o| freevideosex欧美| 99热这里只有是精品在线观看| 亚洲一区二区三区欧美精品| 国产乱来视频区| av播播在线观看一区| 欧美亚洲日本最大视频资源| 国产一区二区三区av在线| 纯流量卡能插随身wifi吗| 国产麻豆69| 丝袜喷水一区| 国产一区亚洲一区在线观看| 亚洲第一av免费看| 大片电影免费在线观看免费| 亚洲人成网站在线观看播放| 日韩中文字幕视频在线看片|