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

    基于狀態(tài)曲線的風電機組運行工況異常檢測

    2019-08-13 08:11:08孫群麗劉長良
    熱力發(fā)電 2019年7期
    關鍵詞:距角變槳風電

    孫群麗,劉長良,周 瑛

    基于狀態(tài)曲線的風電機組運行工況異常檢測

    孫群麗1,劉長良2,周 瑛3

    (1.華北電力大學科技學院,河北 保定 071003;2.華北電力大學新能源電力系統(tǒng)國家重點實驗室,北京 102206;3.石家莊鐵道大學四方學院,河北 石家莊 051132)

    為了提高風電機組的利用率和發(fā)電量,減少風電機組維修和更換費用,需要對其運行工況進行監(jiān)測。本文首先對機組的風速-功率、風速-轉速、風速-槳距角、轉速-功率、轉速-槳距角5種狀態(tài)曲線進行理論介紹,然后結合實際運行數(shù)據(jù)對其進行了分析。結果表明:由于風速的隨機性和風電機組的慣性,前3種曲線不能很好地區(qū)分機組的正常運行狀態(tài)和故障狀態(tài),而轉速-功率、轉速-槳距角能夠對機組的異常情況進行準確的監(jiān)測;以轉速-功率、轉速-槳距角狀態(tài)曲線為基礎,分析了機組不同運行工況在狀態(tài)曲線上的分布,對各個不同工況分別建立相應的評價體系,通過故障實例分析,表明本文方法能提前感知異常情況,有效提高系統(tǒng)的狀態(tài)監(jiān)測精度。

    風電機組;狀態(tài)曲線;轉速-功率;轉速-槳距角;異常檢測;運行工況;監(jiān)控與數(shù)據(jù)采集系統(tǒng)

    風電機組的可靠性關乎著風電場的安全經(jīng)濟運行,因而對風電機組的運行狀態(tài)進行監(jiān)測就顯得尤為重要[1]。文獻[2-3]利用自適應神經(jīng)模糊推理系統(tǒng)(ANFIS)分別對機組的變槳系統(tǒng)和整機進行狀態(tài)監(jiān)視,建立了故障模型,針對故障的情況給出警告的輸出值,但他們都是利用機組已知故障來建立模型,在實際診斷中對故障數(shù)據(jù)和故障類型有一定的要求,因此應用受到一定的限制。文獻[4]從機組運行數(shù)據(jù)出發(fā),利用Relief特征參數(shù)選擇的方法,找出與變槳系統(tǒng)密切相關的運行參數(shù),在這些數(shù)據(jù)的基礎上建立變槳系統(tǒng)的支持向量機回歸模型對變槳系統(tǒng)的狀態(tài)進行監(jiān)測。文獻[5-9]用數(shù)據(jù)挖掘的方法建立了風電機組的非線性識別系統(tǒng),對機組正常運行時的狀態(tài)進行監(jiān)測,一旦發(fā)現(xiàn)異常提前發(fā)出報警信息;但是在使用機組監(jiān)控和數(shù)據(jù)采集系統(tǒng)(SCADA)數(shù)據(jù)來建立訓練模型時,由于數(shù)據(jù)是某一段時間內(nèi)的平均值,這就造成模型的精確度不高。文獻[10]利用Copula函數(shù)建立了風電機組功率曲線的概率模型,然后對機組的運行數(shù)據(jù)進行分析,結果表明可對機組的葉片退化、偏航和變槳的異常情況進行檢測。文獻[11]為了對機組的有功功率數(shù)據(jù)進行分析,采用了Hotelling T2的方法對機組的運行狀態(tài)進行識別。文獻[12]利用計算流體動力學(CFD)技術得到修正后的風速建立功率特性曲線,用此來對機組的性能進行評價;文獻[13]用INNER-DBSCAN算法和功率曲線模型的數(shù)據(jù)驅動實時檢測方法,對數(shù)據(jù)先進行清洗和聚類,然后構造正常運行時的功率曲線模型來識別風機的異常情況。文獻[14]以風力機功率曲線為基礎,建立了風力機效能轉換的指標,對風資源與風場的發(fā)電量進行了準確的描述;但是由于風速的快速隨機的變化,以及設備的惰性造成風速-功率曲線在使用中存在很大的局限性。

    結合上述的研究思想,本文首先從理論上利用機組的技術參數(shù)和運行原理分析機組幾種典型的狀態(tài)曲線,指出它們在確定機組運行狀態(tài)時的優(yōu)缺點,最后選擇轉速-功率,轉速-槳距角狀態(tài)曲線來對機組的運行工況進行監(jiān)測,并建立相應的評判標準。

    1 狀態(tài)曲線

    1.1 理論狀態(tài)曲線

    風力發(fā)電機組是利用葉輪從風中獲取能量通過傳動系統(tǒng)帶動發(fā)電機發(fā)電,從而產(chǎn)生電能。目前并網(wǎng)運行的風力發(fā)電機組均采用變速恒頻的方式,能夠最大限度地吸收風能,提高機組的利用效率,該機組在整個運行過程中包括空載、并網(wǎng)、最大風能捕獲、恒轉速運行和恒功率運行幾個階段,不同運行階段的風速、功率、槳距角、發(fā)電機轉速之間的關系如圖1所示。圖1中in為切入風速,out為切出風速,rtd為額定風速,nrtd為額定轉速下對應的風速;rtd為額定功率,in為并網(wǎng)轉速,rtd為額定轉速;1為機組正常運行過程中的最大槳距角(一般在30°左右),2為機組處于待機狀態(tài)的槳距角(一般在50°左右),off為機組處于順槳狀態(tài)的槳距角(一般在90°左右[15])。

    1.2 實際狀態(tài)曲線

    本文采用某風電場機組5天正常運行的數(shù)據(jù)來描述上述5種性能曲線,結果如圖2所示。

    圖2a)為實際風速-功率運行曲線。圖形標注的部分表示運行中偏離嚴重的部位,這些偏離一方面由設備故障導致,另一方面風速的變動也會導致狀態(tài)曲線的偏移,如虛線圓標注部分。當風速降到額定風速rtd以下,此時由于轉子的慣性輸出的功率依然是額定功率rtd;在黑色方框標注部分也存在同樣的情況,此時風速已經(jīng)小于切入風速in,本應該脫離電網(wǎng),功率小于0,但顯示的功率卻大于0;此外,在橢圓區(qū)域所包圍的散點圖右上部分,當風速在增加的過程中高于額定風速rtd時,同樣由于轉子的慣性,導致實際的功率偏低。因此,用風速-功率曲線來判斷系統(tǒng)的性能顯然不合適。

    圖2b)為風速-葉輪轉速的實際運行曲線。曲線中方塊標注部分是當風速小于切入風速時轉速還不為0,圓形標注部分是當風速降低到nrtd以下時,轉速卻仍然為額定轉速rtd,造成這種現(xiàn)象的原因,可能是因為轉子的慣性[16-17],也可能是因為機組故障。

    圖2c)為風速-槳距角的實際運行曲線。圖中方形表示當風速在切入風速in和額定風速rtd之間時槳距角大于0的部分,引起這部分偏移的原因可能是設備故障,也可能是風速變化較快,而變槳系統(tǒng)有一定滯后性所致[17];橢圓表示當風速10 min平均值小于脫網(wǎng)風速in時或者是發(fā)電機輸出功率達到一定負值時,風電機組不允許長期在電網(wǎng)中運行,必須脫網(wǎng),處于自由狀態(tài),槳距角處于順槳狀態(tài)[18],這是機組正常運行現(xiàn)象,但是當機組出現(xiàn)故障時也會有上述的運行過程。

    圖2d)和圖2e)分別為葉輪轉速-功率,葉輪轉速-槳距角實際運行曲線。對比圖1d)和圖1e)可見,實際運行曲線與理論曲線有很好的匹配度,正常運行狀態(tài)下偏離理論曲線的距離很小。圖2e)中圓形和向上的箭頭表示當風速小于切入風速時,機組處于待機狀態(tài)的運動過程,當風速達到足以吹動葉輪旋轉并滿足啟動條件時,風機開始進入啟動階段,槳距角沿著向下箭頭的方向由順槳狀態(tài)減小到達橢圓所示的部分進行加速,接著槳距角繼續(xù)減小,最終達到最小值[18]。

    2 基于狀態(tài)曲線的運行工況分析

    在分析比較上述5種狀態(tài)曲線的基礎上,選取葉輪轉速-功率、葉輪轉速-槳距角狀態(tài)曲線對機組在不同階段的工況進行分析。

    以變槳系統(tǒng)的運行狀況為例,當機組運行在最佳功率跟蹤階段時,如果出現(xiàn)變槳故障導致槳距角增大,最佳風能利用系數(shù)減小,功率和驅動力矩減小,這將導致轉速下降[19],如圖3中“.”所示的區(qū)域;當機組運行在恒轉速和恒功率階段時,變槳系統(tǒng)如果出現(xiàn)故障導致槳距角減小,將使功率增大,轉速升高[17],最終也將偏離理論運行曲線,如圖3中“o”和“+”所示區(qū)域。

    3 故障分析過程

    3.1 基于狀態(tài)曲線的故障監(jiān)測模型

    對比圖2和圖3中的葉輪轉速-功率、葉輪轉速-槳距角的狀態(tài)曲線可知,正常運行數(shù)據(jù)與理論曲線之間存在一定的距離,這些距離在正常的偏離范圍,當需要檢測的數(shù)據(jù)與理論狀態(tài)曲線的偏離超過這個范圍時可認定為故障。故障監(jiān)測流程如圖4所示。

    圖4 故障監(jiān)測流程

    首先通過機組的技術參數(shù),建立理論狀態(tài)曲線,然后計算實時在線數(shù)據(jù)到理論曲線的最短距離np、nβ,其中nplimnp)&nβlimnβ)表示若此最短距離np,nβ超過限定的閾值limnp、limnβ,則認為系統(tǒng)出現(xiàn)故障。另外,當風速在不同的范圍時機組運行在不同的狀態(tài),對應理論曲線上的不同階段,因此要有以風速為輸入量的狀態(tài)識別環(huán)節(jié)。

    監(jiān)測點到理論曲線的距離由下式計算:

    式中,為實時監(jiān)測到的點,ai為理論曲線上的點。圖5為不同運行工況下的偏離范圍。

    3.2 最小距離的計算

    由于機組運行參數(shù)的量綱不一致,在計算最小距離時存在諸多不便,因此需要對其進行歸一化,取它們相對于額定參數(shù)的比例作為計算值,表1給出了某1.5 MW機組的相關技術參數(shù)??紤]到機組的狀態(tài)曲線是分段的,實際運行數(shù)據(jù)到曲線的最小距離則可通過求帶約束的最小化問題式(2)來實現(xiàn)。

    表1 技術參數(shù)

    Tab.1 The technical parameters

    其中,函數(shù)方程和1的取值范圍由表1中的機組技術參數(shù)和理論曲線確定。

    3.3 偏離閾值的確定

    由于風速的隨機性,風電機組的運行狀態(tài)在最大功率追蹤、恒轉速控制和恒功率控制、停機和啟動幾個階段之間切換。為了加強識別的精度,應根據(jù)不同階段的特點確定對應轉速-功率、轉速-槳距角的偏離閾值。

    在最大功率追蹤階段(MPPT),cut-in+2<

    式中,limnp為運行工況點到轉速-功率曲線的偏離閾值,limnβ為運行工況點到轉速-槳距角曲線的偏離閾值。

    在恒功率階段(CP),rtd<

    在恒轉速階段(CS),rtd–2<

    機組在運行中會因為風速的變化而處于啟動或者停機過程,在啟動過程中,槳距角連續(xù)減小,而轉速連續(xù)增大,如式(9):

    在停機過程中,槳距角連續(xù)增大,而轉速連續(xù)減小,如式(10):

    當轉速和槳距角出現(xiàn)上述變化趨勢時,是正常的運行現(xiàn)象。

    4 故障實例分析

    為了驗證本文所提出的機組狀態(tài)異常辨識方法的有效性,應用某風場1.5 MW機組的運行數(shù)據(jù)進行測試,該機組在2014年1月15日09:50報“變槳安全鏈故障”。變槳安全鏈是把可能對機組造成嚴重損害的故障節(jié)點用繼電器串聯(lián)起來,一旦某個節(jié)點出現(xiàn)故障,則整條線路斷開變槳系統(tǒng)進入緊急停機過程,迅速向90°順槳。運行人員通過排查確認故障是變槳滑環(huán)接觸不良所致。變槳滑環(huán)負責主控和變槳系統(tǒng)之間的信息連接,當其接觸不良時,隨著轉速增大信號丟失概率增加,主控將失去對變槳系統(tǒng)的控制,出現(xiàn)變槳故障。選取該機組前5天的運行數(shù)據(jù),對應的葉輪轉速-功率,葉輪轉速-槳距角關系如圖6所示。從圖6可以看出:在系統(tǒng)報故障前有大量的數(shù)據(jù)偏離正常運行的范圍,這些偏離的點是從風速超過額定值的恒功率運行階段開始的,隨著風速的減小,由于變槳滑環(huán)的故障,槳距角不能有效地減小,從而使風輪獲得的風能減小,轉速偏低,功率偏小,如圖6a)中的圓形所示;當風速繼續(xù)減小,機組進入恒轉速,最佳功率跟蹤階段時,槳距角的偏大使運行數(shù)據(jù)的偏離如圖6a)中虛線橢圓和實線橢圓形所示;同時槳距角的偏離如圖6b)中橢圓區(qū)域所示,箭頭向上和向下的區(qū)域是由于風速在切入值附近變化時引起的啟動和停機過程,滿足式(9)和式(10)。

    基于上述故障檢測流程,當偏離距離np和nβ都小于限定的閾值時,系統(tǒng)處于正常運行狀態(tài),顯示為1;當某一個或者兩個大于閾值時,系統(tǒng)處于故障狀態(tài),顯示為0。從2014年1月10日開始每隔1 min取1次數(shù)據(jù),共取得8 640組數(shù)據(jù)。圖7為異常時間顯示,其中圖7a)為SCADA系統(tǒng)的異常報警時間,圖7b)為本文方法的異常報警時間。從圖7可以看出:在實際運行過程中可能會有偶爾超過閾值的情況,不用發(fā)出報警信號,如圖7b)中橢圓標注的區(qū)域;只有在2014年1月13日09:23持續(xù)出現(xiàn)“0”時系統(tǒng)才發(fā)出報警信息;本文方法在系統(tǒng)出現(xiàn)輕微偏離時即發(fā)出報警信號,比SCADA系統(tǒng)提前48 h 27 min。

    5 結 論

    1)通過對葉輪轉速-功率,葉輪轉速-槳距角狀態(tài)曲線在理論上和實際中正常運行情況下的對比發(fā)現(xiàn),二者具有較高的一致性,能夠較好地反映機組的異常工況。

    2)利用狀態(tài)曲線進行機組的異常情況監(jiān)測時能夠知道異常發(fā)生的機理和過程。

    3)對不同的運行工況分別建立評價閾值,可提高運行狀態(tài)異常識別的精度。

    4)本文方法在系統(tǒng)出現(xiàn)輕微偏離時即可發(fā)出報警信號,比SCADA系統(tǒng)提前48 h 27 min。

    [1] 陳雪峰, 李繼猛, 程航, 等. 風力發(fā)電機狀態(tài)監(jiān)測和故障診斷技術的研究與進展[J]. 機械工程學報, 2011, 47(9): 45-52.CHEN Xuefeng, LI Jimeng, CHENG Hang, et al. Research and application of condition monitoring and fault diagnosis technology in wind turbines[J]. Journal of Mechanical Engineering, 2011, 47(9): 45-52.

    [2] ZHOU Q, XIONG T T, WANG M B, et al. Diagnosis and early warning of wind turbine faults based on cluster analysis theory and modified ANFIS[J]. Energies, 2017, 10(7): 898.

    [3] CHEN B, MATTHEWS P, TAVNER P. Automated on-line fault prognosis for wind turbine pitch systems using supervisory control and data acquisition[J]. IET Renewable Power Generation, 2014, 9(5): 503-513.

    [4] 李輝, 楊超, 李學偉, 等. 風機電動變槳系統(tǒng)狀態(tài)特征參量挖掘及異常識別[J]. 中國電機工程學報, 2014, 34(12): 1922-1930. LI Hui, YANG Chao, LI Xuewei, et al. Conditions characteristic parameters mining and outlier identification for electric pitch system of wind turbine[J]. Proceedings of the CSEE, 2014, 34(12): 1922-1930.

    [5] BI R, ZHOU C, HEPBURN D M, et al. A NBM based on P-N relationship for DFIG wind turbine fault detection[C]. 2015 International Conference on Smart Grid and Clean Energy Technologies (ICSGCE). IEEE, 2016.

    [6] 董興輝, 閆慧麗, 張曉亮. 基于FDA貢獻圖的風電機組變槳系統(tǒng)故障定位[J]. 可再生能源, 2017, 35(1): 93-100. DONG Xinghui, YAN Huili, ZHANG Xiaoliang. Fault location of pitch system based on FDA contribution diagram[J]. Renewable Energy Resources. 2017, 35(1): 93-100.

    [7] 肖成, 劉作軍, 張磊. 基于SCADA系統(tǒng)的風電變槳故障預測方法研究[J]. 可再生能源, 2017, 35(2): 278-284. XIAO Cheng, LIU Zuojun, ZHANG Lei. Variable pitch fault prediction of wind power system based on SCADA system[J]. Renewable Energy Resources, 2017, 35(2): 278-284.

    [8] 巨林倉, 史貝貝, 楊清宇, 等. 基于LM算法建立風電機組神經(jīng)網(wǎng)絡故障預警診斷模型[J]. 熱力發(fā)電, 2010, 39(12): 44-49. JU Lincang, SHI Beibei, YANG Qingyu, et al. Establishment of fault warning diagnosis model for wind power units based on lm algorithm of neural network[J]. Thermal Power Generation, 2010, 39(12): 44-49.

    [9] 徐創(chuàng)學, 劉魯清, 謝云明, 等. 基于實時數(shù)據(jù)的風力發(fā)電機組性能在線評價[J]. 熱力發(fā)電, 2015, 44(4): 88-91. XU Chuangxue, LIU Luqing, XIE Yunming, et al. Real-time data based online evaluation of output performance for wind turbine units[J]. Thermal Power Generation, 2015, 44(4): 88-91.

    [10] WANG Y, INFIELD D, STEPHEN B, et al. Copula- based model for wind turbine power curve outlier rejection[J]. Wind Energy, 2014, 17(11): 1677-1688.

    [11] KUSIAK A, VERMA A. Monitoring wind farms with performance curves[J]. IEEE Transactions on Sustainable Energy, 2013, 4(1): 192-199.

    [12] 邢作霞, 張鵬飛, 李媛. 基于CFD技術的風電機組功率特性評估方法[J]. 太陽能學報, 2015, 36(12): 2882-2887. XING Zuoxia, ZHANG Pengfei, LI Yuan. Evaluation method of power characteristics of wind turbine based on CFD technology[J]. Acta Energiae Solaris Sinica, 2015, 36(12): 2882-2887.

    [13] 陳警鈺, 陳玉航. 基于INNER-DBSCAN和功率曲線模型的風機異常狀態(tài)檢測[J]. 電力科學與工程, 2017, 33(8): 27-34.CHEN Jingyu, CHEN Yuhang. Condition monitoring for wind turbines based on INNER-DBSCAN and power curvepattern[J]. Electric Power Science and Engineering, 2017, 33(8): 27-34.

    [14] 張家安, 劉寧, 吳林林, 等. 基于風力機功率特性的風資源評估指標體系[J]. 太陽能學報, 2017, 38(2): 498-502. ZHANG Jiaan, LIU Ning, WU Linlin, et al.index system for wind resource assessment based on power characteristics of wind turbine[J]. Acta Energiae Solaris Sinica, 2017, 38(2): 498-502.

    [15] 夏長亮. 雙饋風力發(fā)電系統(tǒng)設計與并網(wǎng)運行[M]. 北京: 科學出版社, 2014: 22-23.XIA Changliang. The design and grid-connected operations for double fed wind power generation system[M]. Beijing: Science Press, 2014: 22-23.

    [16] 王明軍. 風電機組功率曲線問題及爭議淺析[J]. 風能, 2016(5): 42-45. WANG Mingjun. Power curve problem and dispute analysis of wind turbines[J]. Wind Energy, 2016(5): 42-45.

    [17] 王明軍, 高原生. 風電機組實際運行功率曲線影響因素分析[J]. 風能, 2013(4): 74-79.WANG Mingjun, GAO Yuansheng. An analysis of influence factors on wind turbine’s actual power curve[J]. Wind Energy, 2013(4): 74-79.

    [18] 吳雙群. 風力發(fā)電原理[M]. 北京: 北京大學出版社, 2011: 77.WU Shuangqun. Principle of the wind driven generator[M]. Beijing: Peking University Press, 2011: 77.

    [19] 趙振宙. 風力機原理[M]. 北京: 中國水利水電出版社, 2016: 45. ZHAO Zhenzhou. Principle of wind turbine[M]. Beijing: China Water & Power Press, 2016: 45.

    Abnormal detection of operation conditions of wind turbine based on state curve

    SUN Qunli1, LIU Changliang2, ZHOU Ying3

    (1. Science and Technology College, North China Electric Power University, Baoding 071003, China; 2. State Key Laboratory of Alternate Electrical Power System with Renewable Energy Sources, North China Electric Power University, Beijing 102206, China; 3. Sifang College, Shijiazhuang Tiedao University, Shijiazhuang 051132, China)

    In order to improve the utilization rate and power generation of wind turbines, and reduce the cost of maintenance and replacement of wind turbines, it is necessary to monitor the operating conditions. In this paper, five state curves of the unit, including the wind speed-power curve, the wind speed-rotor speed curve, the wind speed-pitch angle curve, the rotor speed-power curve, and the rotor speed-pitch angle curve, are introduced theoretically. Then, they are analyzed combining with the actual operational data. The results show that, due to the randomness of wind speed and the inertia of wind turbines, the first three curves does not allow the wind turbine to distinguish normal operation status from a fault status. However, the rotor speed-power curve and rotor speed-pitch angle curve are able to accurately monitor the abnormal conditions of the wind turbine. On the basis of the last two state curves, the distribution of different operating conditions of the unit is analyzed, and corresponding evaluation systems are established for each operating condition. The fault case study indicates that, this method can sense the unit's abnormal situation in advance, and effectively improve the state monitoring accuracy of the system.

    wind turbine unit, state curve, rotor speed-power, rotor speed-pitch angle, anomaly detection, operating condition, monitoring and data acquisition system

    TK83

    A

    10.19666/j.rlfd.201901028

    孫群麗, 劉長良, 周瑛. 基于狀態(tài)曲線的風電機組運行工況異常檢測[J]. 熱力發(fā)電, 2019, 48(7): 110-116. SUN Qunli, LIU Changliang, ZHOU Ying. Abnormal detection of operation conditions of wind turbine based on state curve[J]. Thermal Power Generation, 2019, 48(7): 110-116.

    2019-01-26

    中央高校基本科研業(yè)務費專項基金資助項目(9161717007);北京市自然科學基金(4182061)

    Supported by:Fundamental Research Funds for the Central Universities (9161717007); Natural Science Foundation of Beijing (4182061)

    孫群麗(1982—),女,講師,博士研究生,主要研究方向為電廠設備的智能控制和故障診斷技術,hdsunqunli@163.com。

    (責任編輯 馬昕紅)

    猜你喜歡
    距角變槳風電
    基于數(shù)據(jù)驅動的風電機組最優(yōu)槳距角 辨識方法
    考慮風速差異的風電場減載方案與一次調(diào)頻策略
    電力建設(2022年7期)2022-07-04 02:50:50
    風力發(fā)電機組最優(yōu)槳距角自尋優(yōu)控制算法設計
    東方汽輪機(2020年3期)2020-11-04 11:07:34
    兆瓦級風電機組變槳距系統(tǒng)設計
    大電機技術(2018年4期)2018-09-26 03:11:38
    海上風電躍進隱憂
    能源(2018年6期)2018-08-01 03:42:00
    分散式風電破“局”
    能源(2018年6期)2018-08-01 03:41:56
    風電:棄風限電明顯改善 海上風電如火如荼
    能源(2018年8期)2018-01-15 19:18:24
    變速風力發(fā)電機組變槳距建模與仿真
    基于最小化多變量的獨立變槳距控制研究
    電源技術(2016年2期)2016-02-27 09:05:15
    重齒風電
    風能(2016年12期)2016-02-25 08:46:38
    99国产精品一区二区蜜桃av| 亚洲av五月六月丁香网| 久久久精品大字幕| 国内揄拍国产精品人妻在线| 老司机福利观看| 久久精品国产亚洲av香蕉五月| 欧美一级a爱片免费观看看 | 久久人妻福利社区极品人妻图片| 大型黄色视频在线免费观看| 怎么达到女性高潮| tocl精华| 老司机深夜福利视频在线观看| 亚洲,欧美精品.| 欧美另类亚洲清纯唯美| 国产91精品成人一区二区三区| 一本大道久久a久久精品| 成人午夜高清在线视频| 亚洲av美国av| 香蕉久久夜色| 国产av一区在线观看免费| 午夜精品在线福利| 国产精品久久久久久精品电影| 日本免费一区二区三区高清不卡| 亚洲av成人精品一区久久| 国产视频内射| 99国产精品一区二区蜜桃av| 亚洲乱码一区二区免费版| 欧美一级a爱片免费观看看 | 精品久久久久久久久久免费视频| 日日干狠狠操夜夜爽| 免费无遮挡裸体视频| 国产成人啪精品午夜网站| 欧美成人性av电影在线观看| 婷婷精品国产亚洲av在线| 女人爽到高潮嗷嗷叫在线视频| 一进一出抽搐动态| 国产一区二区在线av高清观看| 国产午夜精品论理片| 欧美一级毛片孕妇| 亚洲精华国产精华精| 男女床上黄色一级片免费看| 嫩草影视91久久| 国产午夜精品论理片| 少妇的丰满在线观看| 嫁个100分男人电影在线观看| 五月玫瑰六月丁香| 国内少妇人妻偷人精品xxx网站 | 蜜桃久久精品国产亚洲av| 欧美日韩福利视频一区二区| 少妇裸体淫交视频免费看高清 | 一级毛片女人18水好多| 女人被狂操c到高潮| 一本一本综合久久| 舔av片在线| 女生性感内裤真人,穿戴方法视频| 91麻豆精品激情在线观看国产| 色在线成人网| 又粗又爽又猛毛片免费看| 久久久久国产精品人妻aⅴ院| 又黄又粗又硬又大视频| 日韩精品中文字幕看吧| 久久国产精品人妻蜜桃| 久久久久国产一级毛片高清牌| 国产熟女午夜一区二区三区| 丝袜人妻中文字幕| av视频在线观看入口| 久久久久久久午夜电影| 99久久精品热视频| 国产三级中文精品| 国产99白浆流出| 欧美乱色亚洲激情| 欧美一级a爱片免费观看看 | 夜夜躁狠狠躁天天躁| 可以免费在线观看a视频的电影网站| 久久精品国产亚洲av高清一级| 床上黄色一级片| 在线观看一区二区三区| 亚洲国产精品sss在线观看| 欧美日韩黄片免| 91九色精品人成在线观看| 日韩大尺度精品在线看网址| 美女扒开内裤让男人捅视频| 欧美丝袜亚洲另类 | 国产日本99.免费观看| 免费观看人在逋| 亚洲成人国产一区在线观看| www.自偷自拍.com| 午夜福利高清视频| 一本综合久久免费| 国产aⅴ精品一区二区三区波| 男人舔奶头视频| 欧美最黄视频在线播放免费| 99riav亚洲国产免费| 欧美日韩中文字幕国产精品一区二区三区| 97碰自拍视频| 婷婷亚洲欧美| 欧美中文综合在线视频| 99在线视频只有这里精品首页| 一个人免费在线观看的高清视频| 在线免费观看的www视频| av在线播放免费不卡| 亚洲精品在线观看二区| 国产精品亚洲一级av第二区| 脱女人内裤的视频| 午夜成年电影在线免费观看| 亚洲精品国产一区二区精华液| 俄罗斯特黄特色一大片| 亚洲男人天堂网一区| 一进一出好大好爽视频| 在线观看免费视频日本深夜| 看免费av毛片| 精品国产乱子伦一区二区三区| 欧美成人免费av一区二区三区| 久久久久久亚洲精品国产蜜桃av| 一个人免费在线观看电影 | 国内揄拍国产精品人妻在线| 少妇被粗大的猛进出69影院| 亚洲欧美日韩无卡精品| 丁香六月欧美| 可以免费在线观看a视频的电影网站| 男人的好看免费观看在线视频 | cao死你这个sao货| 成人国产一区最新在线观看| 成人国产综合亚洲| 国产熟女xx| 欧美乱码精品一区二区三区| 欧美黄色淫秽网站| 最近在线观看免费完整版| cao死你这个sao货| 国产精品免费一区二区三区在线| 老熟妇乱子伦视频在线观看| 欧美中文日本在线观看视频| 1024手机看黄色片| 亚洲人成77777在线视频| 在线观看66精品国产| 窝窝影院91人妻| 成人国语在线视频| 午夜精品在线福利| 久久精品国产亚洲av高清一级| 搡老妇女老女人老熟妇| 国产av一区二区精品久久| 日本免费a在线| 在线观看日韩欧美| 国产精品一区二区三区四区免费观看 | 欧美大码av| 非洲黑人性xxxx精品又粗又长| 欧美成人性av电影在线观看| 一夜夜www| 国产黄色小视频在线观看| 一个人免费在线观看的高清视频| 俄罗斯特黄特色一大片| 亚洲av熟女| 黄片大片在线免费观看| 日日干狠狠操夜夜爽| 免费观看人在逋| 最近最新中文字幕大全免费视频| 国产一区二区在线av高清观看| 精品免费久久久久久久清纯| 草草在线视频免费看| 男人舔奶头视频| 精华霜和精华液先用哪个| 99久久精品热视频| 岛国视频午夜一区免费看| 18禁国产床啪视频网站| or卡值多少钱| 亚洲av日韩精品久久久久久密| 男人舔女人的私密视频| 久久久久国产精品人妻aⅴ院| 国产伦在线观看视频一区| 色老头精品视频在线观看| 亚洲国产精品合色在线| 又黄又粗又硬又大视频| 手机成人av网站| 999久久久精品免费观看国产| 国产精品久久久久久亚洲av鲁大| 全区人妻精品视频| e午夜精品久久久久久久| 欧洲精品卡2卡3卡4卡5卡区| 精品日产1卡2卡| 午夜精品久久久久久毛片777| 操出白浆在线播放| 久久久久久亚洲精品国产蜜桃av| 欧美日韩亚洲国产一区二区在线观看| 国产精品亚洲美女久久久| 久久久久久免费高清国产稀缺| 久久人人精品亚洲av| a级毛片a级免费在线| 精品欧美国产一区二区三| 两个人看的免费小视频| 欧美成人午夜精品| 巨乳人妻的诱惑在线观看| 国产精品自产拍在线观看55亚洲| 欧美日韩乱码在线| 国产精品日韩av在线免费观看| 国产av在哪里看| 亚洲电影在线观看av| 国产男靠女视频免费网站| 18禁黄网站禁片免费观看直播| 国产精品影院久久| 黄片小视频在线播放| 999精品在线视频| 午夜福利高清视频| 国产高清激情床上av| 中文字幕人妻丝袜一区二区| 国产一区二区三区在线臀色熟女| 变态另类丝袜制服| 19禁男女啪啪无遮挡网站| 90打野战视频偷拍视频| 国产成人系列免费观看| avwww免费| 国产免费av片在线观看野外av| 欧美黄色片欧美黄色片| 黄色 视频免费看| 国产高清视频在线播放一区| 一级黄色大片毛片| 国产伦一二天堂av在线观看| 日本一本二区三区精品| 成人一区二区视频在线观看| 国产成人精品久久二区二区免费| 日本 欧美在线| 又大又爽又粗| 午夜影院日韩av| 在线观看免费视频日本深夜| 国产免费男女视频| 日韩精品中文字幕看吧| 两性夫妻黄色片| 欧美成人性av电影在线观看| 男女下面进入的视频免费午夜| 99热只有精品国产| 久久精品综合一区二区三区| 欧美性猛交╳xxx乱大交人| 制服丝袜大香蕉在线| 午夜福利免费观看在线| 麻豆国产97在线/欧美 | 嫁个100分男人电影在线观看| 在线观看www视频免费| 国产一级毛片七仙女欲春2| 亚洲成人中文字幕在线播放| 日本精品一区二区三区蜜桃| 听说在线观看完整版免费高清| 最近最新中文字幕大全免费视频| 精品欧美国产一区二区三| 可以免费在线观看a视频的电影网站| 国产在线精品亚洲第一网站| 久久人妻av系列| 两性夫妻黄色片| 欧美人与性动交α欧美精品济南到| 在线看三级毛片| 高清在线国产一区| 女人被狂操c到高潮| 一卡2卡三卡四卡精品乱码亚洲| 欧美另类亚洲清纯唯美| 少妇裸体淫交视频免费看高清 | 中文字幕高清在线视频| 欧美成狂野欧美在线观看| 两人在一起打扑克的视频| 最新美女视频免费是黄的| 一个人免费在线观看电影 | 男女之事视频高清在线观看| 久久性视频一级片| 国产成人一区二区三区免费视频网站| 一卡2卡三卡四卡精品乱码亚洲| 欧美不卡视频在线免费观看 | 午夜精品久久久久久毛片777| 国产成人av教育| 高清在线国产一区| 精品久久久久久久久久久久久| 亚洲av成人精品一区久久| 国产成人av教育| 制服人妻中文乱码| 免费看美女性在线毛片视频| 夜夜爽天天搞| 国产av又大| 1024视频免费在线观看| 欧美色欧美亚洲另类二区| 搡老妇女老女人老熟妇| 国产精品一区二区精品视频观看| 日韩免费av在线播放| 国产一区二区在线观看日韩 | 亚洲,欧美精品.| 男女午夜视频在线观看| 欧美av亚洲av综合av国产av| 亚洲真实伦在线观看| 亚洲精品国产精品久久久不卡| 欧美性猛交黑人性爽| av在线天堂中文字幕| 中亚洲国语对白在线视频| 欧美中文日本在线观看视频| 大型av网站在线播放| 久久精品国产清高在天天线| 黄频高清免费视频| 99久久综合精品五月天人人| 中文资源天堂在线| 久久久久免费精品人妻一区二区| 999精品在线视频| 女生性感内裤真人,穿戴方法视频| 日本一区二区免费在线视频| 久久人妻av系列| 成人国产一区最新在线观看| 最近最新中文字幕大全免费视频| 亚洲第一电影网av| 精品少妇一区二区三区视频日本电影| 美女高潮喷水抽搐中文字幕| 99riav亚洲国产免费| 国产精品电影一区二区三区| 老汉色∧v一级毛片| 精品乱码久久久久久99久播| 热99re8久久精品国产| 淫妇啪啪啪对白视频| 麻豆久久精品国产亚洲av| 美女高潮喷水抽搐中文字幕| 日韩大码丰满熟妇| 怎么达到女性高潮| 波多野结衣高清无吗| 听说在线观看完整版免费高清| 色综合欧美亚洲国产小说| 老司机在亚洲福利影院| 色综合亚洲欧美另类图片| 久久性视频一级片| 啦啦啦观看免费观看视频高清| 国产爱豆传媒在线观看 | 免费看日本二区| av免费在线观看网站| 色精品久久人妻99蜜桃| 成人亚洲精品av一区二区| 国产亚洲精品久久久久久毛片| 99久久精品热视频| 国产熟女午夜一区二区三区| 久久精品国产综合久久久| 国产一区在线观看成人免费| 午夜精品在线福利| 国产伦在线观看视频一区| 国产精品久久久久久精品电影| 最近在线观看免费完整版| 在线国产一区二区在线| 男人舔女人下体高潮全视频| 日本一二三区视频观看| 国产亚洲精品av在线| 精华霜和精华液先用哪个| 这个男人来自地球电影免费观看| 国产探花在线观看一区二区| 亚洲乱码一区二区免费版| 久久精品91无色码中文字幕| 亚洲在线自拍视频| 婷婷精品国产亚洲av在线| 久久亚洲真实| 天堂av国产一区二区熟女人妻 | 91麻豆精品激情在线观看国产| 欧美乱妇无乱码| 无限看片的www在线观看| 在线永久观看黄色视频| 黄色片一级片一级黄色片| 色老头精品视频在线观看| 精品乱码久久久久久99久播| 国产精品美女特级片免费视频播放器 | 最近最新免费中文字幕在线| 黑人欧美特级aaaaaa片| ponron亚洲| 国产亚洲精品久久久久5区| 欧美大码av| 此物有八面人人有两片| 男女视频在线观看网站免费 | 变态另类丝袜制服| 美女 人体艺术 gogo| www日本在线高清视频| 国产主播在线观看一区二区| 在线观看日韩欧美| 久久久久免费精品人妻一区二区| 国产精品久久久久久精品电影| 99热6这里只有精品| 亚洲专区国产一区二区| 十八禁网站免费在线| 国产视频内射| 欧美成人午夜精品| 欧美日韩中文字幕国产精品一区二区三区| 婷婷丁香在线五月| ponron亚洲| 母亲3免费完整高清在线观看| 丁香六月欧美| 久久这里只有精品中国| 女警被强在线播放| 亚洲精品一卡2卡三卡4卡5卡| 国产午夜精品久久久久久| 狂野欧美激情性xxxx| 悠悠久久av| 久久亚洲真实| 韩国av一区二区三区四区| 国产欧美日韩一区二区三| 精品一区二区三区四区五区乱码| 亚洲国产精品合色在线| 久久精品综合一区二区三区| 级片在线观看| 日韩欧美精品v在线| 麻豆av在线久日| 在线观看舔阴道视频| 亚洲精品色激情综合| 一进一出抽搐动态| 精品福利观看| 极品教师在线免费播放| 亚洲av成人精品一区久久| 午夜福利在线在线| 亚洲一码二码三码区别大吗| 国产亚洲精品av在线| 长腿黑丝高跟| 国产精品永久免费网站| 欧美一区二区精品小视频在线| 日本a在线网址| 免费一级毛片在线播放高清视频| 少妇被粗大的猛进出69影院| 热99re8久久精品国产| 黑人巨大精品欧美一区二区mp4| 亚洲性夜色夜夜综合| 亚洲男人的天堂狠狠| 欧美日本亚洲视频在线播放| 18禁美女被吸乳视频| 亚洲精品av麻豆狂野| 免费在线观看影片大全网站| 欧美日韩亚洲综合一区二区三区_| 久久久久国内视频| 在线观看美女被高潮喷水网站 | 香蕉久久夜色| 成人av在线播放网站| 久久国产精品人妻蜜桃| 午夜免费成人在线视频| 欧美乱妇无乱码| 日本熟妇午夜| 国产精品久久久人人做人人爽| 两人在一起打扑克的视频| 悠悠久久av| 麻豆av在线久日| 两性午夜刺激爽爽歪歪视频在线观看 | 99久久无色码亚洲精品果冻| 国产午夜精品论理片| 麻豆国产97在线/欧美 | 日韩国内少妇激情av| 91成年电影在线观看| 精品午夜福利视频在线观看一区| 一本精品99久久精品77| 国产精品 欧美亚洲| 久久久久亚洲av毛片大全| 久久精品成人免费网站| 免费在线观看视频国产中文字幕亚洲| 91成年电影在线观看| 99精品在免费线老司机午夜| 在线永久观看黄色视频| 婷婷丁香在线五月| 人人妻人人澡欧美一区二区| 无人区码免费观看不卡| 女同久久另类99精品国产91| 岛国在线免费视频观看| 女人爽到高潮嗷嗷叫在线视频| 久久精品影院6| xxx96com| 欧美午夜高清在线| 亚洲男人的天堂狠狠| av中文乱码字幕在线| 成年人黄色毛片网站| 国产精品一区二区精品视频观看| 国内毛片毛片毛片毛片毛片| 啦啦啦韩国在线观看视频| 88av欧美| 午夜成年电影在线免费观看| 搡老岳熟女国产| 夜夜夜夜夜久久久久| 精品电影一区二区在线| 国产av又大| 一级黄色大片毛片| 两个人视频免费观看高清| 欧美乱码精品一区二区三区| 九九热线精品视视频播放| 怎么达到女性高潮| 免费av毛片视频| 欧美极品一区二区三区四区| 女人高潮潮喷娇喘18禁视频| 别揉我奶头~嗯~啊~动态视频| 全区人妻精品视频| 男女视频在线观看网站免费 | 精品久久久久久久毛片微露脸| 欧美一级a爱片免费观看看 | 可以免费在线观看a视频的电影网站| 91字幕亚洲| 成人av在线播放网站| 亚洲中文字幕一区二区三区有码在线看 | 黄色女人牲交| 欧美成人免费av一区二区三区| 2021天堂中文幕一二区在线观| 国内精品久久久久久久电影| 成人av一区二区三区在线看| www.www免费av| 国产黄片美女视频| 亚洲国产精品成人综合色| 可以在线观看毛片的网站| 哪里可以看免费的av片| 日本一本二区三区精品| 国产精品一及| 国产片内射在线| 男女做爰动态图高潮gif福利片| 欧美中文综合在线视频| 日韩精品免费视频一区二区三区| 首页视频小说图片口味搜索| 亚洲精品在线观看二区| 免费在线观看亚洲国产| 成人av一区二区三区在线看| 青草久久国产| 老汉色∧v一级毛片| 给我免费播放毛片高清在线观看| 久久香蕉精品热| 麻豆av在线久日| 在线a可以看的网站| 在线观看午夜福利视频| 又紧又爽又黄一区二区| 久久精品91无色码中文字幕| 国产精品一区二区精品视频观看| cao死你这个sao货| 曰老女人黄片| 国产黄色小视频在线观看| www日本在线高清视频| 一夜夜www| 婷婷亚洲欧美| 男女视频在线观看网站免费 | 国产高清有码在线观看视频 | 亚洲激情在线av| 国产aⅴ精品一区二区三区波| 国产私拍福利视频在线观看| 国产人伦9x9x在线观看| 欧美午夜高清在线| 成人18禁高潮啪啪吃奶动态图| 亚洲av成人av| 国产99久久九九免费精品| 国产高清视频在线观看网站| 男女那种视频在线观看| 免费av毛片视频| 久久热在线av| 欧美黑人欧美精品刺激| 欧美大码av| 久久中文看片网| 一级毛片女人18水好多| 国产午夜精品久久久久久| 中文字幕久久专区| 亚洲精品美女久久av网站| 国产探花在线观看一区二区| 女生性感内裤真人,穿戴方法视频| 国产激情偷乱视频一区二区| 亚洲国产欧美人成| 国产亚洲精品综合一区在线观看 | 亚洲av中文字字幕乱码综合| 琪琪午夜伦伦电影理论片6080| 99国产精品99久久久久| 一个人免费在线观看的高清视频| 国产亚洲精品av在线| 欧美日韩亚洲综合一区二区三区_| 亚洲一区二区三区不卡视频| 亚洲精品中文字幕一二三四区| 亚洲欧美日韩东京热| 亚洲国产精品成人综合色| 国产乱人伦免费视频| 最近最新中文字幕大全免费视频| 激情在线观看视频在线高清| 精品久久久久久久末码| 又紧又爽又黄一区二区| 日本一二三区视频观看| 麻豆久久精品国产亚洲av| 国产av又大| 又粗又爽又猛毛片免费看| 桃色一区二区三区在线观看| 欧美日本亚洲视频在线播放| 可以免费在线观看a视频的电影网站| 国产一区在线观看成人免费| 午夜福利成人在线免费观看| 国产高清激情床上av| 黄色丝袜av网址大全| 岛国视频午夜一区免费看| 国产99久久九九免费精品| 九九热线精品视视频播放| 一本综合久久免费| 成人午夜高清在线视频| 小说图片视频综合网站| 香蕉国产在线看| 国产蜜桃级精品一区二区三区| 天堂影院成人在线观看| 五月玫瑰六月丁香| 看片在线看免费视频| 欧美性猛交╳xxx乱大交人| 国产伦一二天堂av在线观看| 99国产精品一区二区蜜桃av| 又紧又爽又黄一区二区| tocl精华| 一区二区三区高清视频在线| 欧美色欧美亚洲另类二区| av国产免费在线观看| 法律面前人人平等表现在哪些方面| 日韩 欧美 亚洲 中文字幕| 色av中文字幕| 欧美日本亚洲视频在线播放| 国产成人aa在线观看| 亚洲乱码一区二区免费版| 成人国语在线视频| 国内少妇人妻偷人精品xxx网站 | 久久久久九九精品影院| 两个人的视频大全免费| 精品熟女少妇八av免费久了| 村上凉子中文字幕在线| 一级作爱视频免费观看| 女人高潮潮喷娇喘18禁视频| 特大巨黑吊av在线直播| 1024香蕉在线观看| 免费一级毛片在线播放高清视频| 欧美性猛交╳xxx乱大交人| 欧美精品啪啪一区二区三区| 一级片免费观看大全| 999精品在线视频| 亚洲人成伊人成综合网2020| 欧美在线黄色| 一本精品99久久精品77| 19禁男女啪啪无遮挡网站|