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

    抽水蓄能機(jī)組甩負(fù)荷時(shí)無葉區(qū)壓力幅值特性

    2023-08-27 13:29:48趙毅鋒張飛陳善貴徐亞鵬
    人民黃河 2023年7期

    趙毅鋒 張飛 陳善貴 徐亞鵬

    摘 要:甩負(fù)荷是抽水蓄能機(jī)組的正常設(shè)計(jì)工況,此時(shí)無葉區(qū)經(jīng)歷過流系統(tǒng)最為嚴(yán)重的壓力脈動(dòng),對(duì)機(jī)組運(yùn)行穩(wěn)定產(chǎn)生深遠(yuǎn)影響。長(zhǎng)期以來,對(duì)無葉區(qū)壓力及其脈動(dòng)研究主要集中于數(shù)值模擬、模型觀測(cè),對(duì)真機(jī)試驗(yàn)研究少。以某抽水蓄能電站6 臺(tái)機(jī)組單機(jī)甩額定負(fù)荷時(shí)無葉區(qū)壓力信號(hào)為分析對(duì)象,基于Savitzky-Golay 濾波器對(duì)壓力信號(hào)進(jìn)行處理,分離出了趨勢(shì)壓力與脈動(dòng)壓力,闡釋了趨勢(shì)壓力與脈動(dòng)壓力的主要影響因素,并指明了亟須解決的問題。研究表明:甩負(fù)荷過程中無葉區(qū)趨勢(shì)壓力主要受轉(zhuǎn)速影響,與轉(zhuǎn)速的相關(guān)系數(shù)在0.99 左右;脈動(dòng)壓力峰峰值受動(dòng)靜干涉、旋轉(zhuǎn)失速等影響,是額定負(fù)荷穩(wěn)定運(yùn)行時(shí)的7~ 22 倍,顯著高于數(shù)值模擬結(jié)果。

    關(guān)鍵詞:抽水蓄能機(jī)組;甩負(fù)荷;無葉區(qū);趨勢(shì)壓力;脈動(dòng)壓力

    中圖分類號(hào):TH734 文獻(xiàn)標(biāo)志碼:A doi:10.3969/ j.issn.1000-1379.2023.07.025

    引用格式:趙毅鋒,張飛,陳善貴,等.抽水蓄能機(jī)組甩負(fù)荷時(shí)無葉區(qū)壓力幅值特性[J].人民黃河,2023,45(7):134-139.

    隨著我國(guó)“碳達(dá)峰、碳中和”政策的實(shí)施,抽水蓄能發(fā)展步入了快車道。截至2021 年8 月,我國(guó)建成投運(yùn)抽水蓄能電站總裝機(jī)容量3 249 萬kW,在建抽水蓄能電站總裝機(jī)容量5 373 萬kW,在運(yùn)、在建電站規(guī)模均居世界第一。根據(jù)《抽水蓄能中長(zhǎng)期發(fā)展規(guī)劃(2021—2035)》,到2025 年,抽水蓄能投產(chǎn)總規(guī)模6 200萬kW 以上;到2030 年,投產(chǎn)總規(guī)模1.2 億kW左右。光伏、風(fēng)電等間歇性可再生能源的大規(guī)模并網(wǎng),抽水蓄能機(jī)組運(yùn)行強(qiáng)度顯著提升,相關(guān)影響機(jī)組運(yùn)行安全穩(wěn)定的現(xiàn)象引起普遍關(guān)注,其中無葉區(qū)壓力脈動(dòng)導(dǎo)致的機(jī)組及廠房問題表現(xiàn)最為凸出,如張河灣抽水蓄能電站因無葉區(qū)動(dòng)靜干涉而引發(fā)的廠房局部共振現(xiàn)象[1] 、黑麋峰抽水蓄能電站機(jī)組無葉區(qū)壓力脈動(dòng)大導(dǎo)致的轉(zhuǎn)輪裂紋及廠房振動(dòng)問題[2] 等。

    無葉區(qū)壓力脈動(dòng)是引起抽水蓄能機(jī)組振動(dòng)問題的主要激振源,是水泵水輪機(jī)流道系統(tǒng)中幅值最大的壓力脈動(dòng)[3] ,其主要表現(xiàn)形式包括旋轉(zhuǎn)失速、動(dòng)靜干涉等。傳統(tǒng)上,研究人員主要關(guān)心穩(wěn)態(tài)情況下的無葉區(qū)壓力脈動(dòng),且主要采用模型試驗(yàn)和流動(dòng)計(jì)算方法進(jìn)行研究,模型試驗(yàn)方面:Hasmatuchi 等[4] 通過模型試驗(yàn)分別觀測(cè)到水泵水輪機(jī)發(fā)電工況下偏離最優(yōu)工況和泵工況下的旋轉(zhuǎn)失速現(xiàn)象,Ran 等[5] 通過模型試驗(yàn)分別觀測(cè)到了水泵水輪機(jī)在水泵工況和水輪機(jī)工況的流動(dòng)不穩(wěn)定現(xiàn)象;流動(dòng)計(jì)算方面:李德友[6] 采用剪切應(yīng)力輸運(yùn)SST k-w 模型分別實(shí)現(xiàn)了原型和模型水泵水輪機(jī)在水泵工況下無葉區(qū)壓力脈動(dòng)的流動(dòng)分析;Zhang 等[7]采用不同的湍流模型針對(duì)S 區(qū)內(nèi)的流動(dòng)進(jìn)行了分析,獲得了S 區(qū)形成的機(jī)理。近年來,隨著計(jì)算流體動(dòng)力學(xué)的發(fā)展以及計(jì)算能力的提升,采用一維特征線法結(jié)合三維流動(dòng)計(jì)算法針對(duì)甩負(fù)荷過程進(jìn)行分析的手段逐漸展開,如:Zhang 等[8] 采用該方法對(duì)模型水泵水輪機(jī)甩負(fù)荷進(jìn)行建模分析,獲得了甩負(fù)荷過程中流道內(nèi)壓力脈動(dòng)、轉(zhuǎn)輪受力情況; Yin 等[9] 對(duì)抽水蓄能電站雙機(jī)同時(shí)甩負(fù)荷工況進(jìn)行了研究,指出雙甩較單甩導(dǎo)致更高的無葉區(qū)壓力脈動(dòng);張曉曦等[10] 指出了甩負(fù)荷過程中尾水管水柱分離的機(jī)理及其嚴(yán)重影響。同時(shí),針對(duì)模型機(jī)組,研究者們[11-12] 開展了過渡過程無葉區(qū)壓力脈動(dòng)及其流場(chǎng)測(cè)試試驗(yàn)研究;針對(duì)原型機(jī)組,李立等[13] 、楊建東等[14] 研究了甩負(fù)荷過程中的無葉區(qū)壓力脈動(dòng)特征,結(jié)果表明,多數(shù)研究采用模型試驗(yàn)對(duì)數(shù)值方法進(jìn)行驗(yàn)證,并在此基礎(chǔ)上展開過渡過程工況流動(dòng)機(jī)理研究。考慮到水電機(jī)組一次開停機(jī)過程等效機(jī)組正常運(yùn)行15 h 的損傷[15] ,而甩負(fù)荷過程是較開停機(jī)過程更大風(fēng)險(xiǎn)的過渡過程工況,對(duì)機(jī)組的損傷將更嚴(yán)重;同時(shí),在甩負(fù)荷過程中,無葉區(qū)壓力是整個(gè)過流系統(tǒng)中壓力上升最大的區(qū)域。因此,對(duì)甩負(fù)荷這類過渡過程工況壓力及其脈動(dòng)進(jìn)行研究,有助于工程技術(shù)人員認(rèn)識(shí)甩負(fù)荷對(duì)機(jī)組造成的損傷,對(duì)于提升機(jī)組預(yù)期壽命有很大幫助。然而受限于原型試驗(yàn)的復(fù)雜性及風(fēng)險(xiǎn)性,針對(duì)原型機(jī)組過渡過程壓力脈動(dòng),特別是甩負(fù)荷下無葉區(qū)壓力脈動(dòng)原型機(jī)組觀測(cè)較少,通常只在機(jī)組投入商業(yè)運(yùn)行前的調(diào)試階段進(jìn)行這類試驗(yàn),且相關(guān)研究局限于單次觀測(cè)結(jié)果分析,對(duì)于數(shù)據(jù)的重復(fù)性缺乏足夠分析。

    2019 年8 月—2021 年2 月,某抽水蓄能電站完成了6 臺(tái)機(jī)組的啟動(dòng)調(diào)試任務(wù),調(diào)試過程中每臺(tái)機(jī)組均進(jìn)行了甩負(fù)荷試驗(yàn)。為掌握機(jī)組性能,每臺(tái)機(jī)組均在相同位置就近設(shè)置了無葉區(qū)壓力脈動(dòng)測(cè)點(diǎn),因而對(duì)其進(jìn)行研究,以掌握甩負(fù)荷過程中的無葉區(qū)壓力脈動(dòng)規(guī)律。此成果有助于頂蓋、轉(zhuǎn)輪等關(guān)鍵部件受力分析及疲勞評(píng)估,指導(dǎo)實(shí)現(xiàn)機(jī)組安全運(yùn)行評(píng)價(jià),對(duì)領(lǐng)域內(nèi)研究學(xué)者及工程技術(shù)人員有重要的參考價(jià)值。

    1 試驗(yàn)概況

    1.1 電站簡(jiǎn)介

    某抽水蓄能電站安裝6 臺(tái)單機(jī)容量為300 MW 的立軸單級(jí)混流式水泵水輪機(jī)-發(fā)電電動(dòng)機(jī)組,年均發(fā)電量30.15 億kW·h,年均抽水用電量40.2 億kW·h。該電站在電網(wǎng)中擔(dān)任調(diào)峰、調(diào)頻、調(diào)相、黑啟動(dòng)和事故備用等任務(wù)。電站引水系統(tǒng)采用3 洞6 機(jī)斜井式布置,引水調(diào)壓室采用阻抗+上室結(jié)構(gòu)形式,布置在引水隧洞末端;尾水系統(tǒng)采用2 機(jī)合1 洞布置形式,尾水調(diào)壓室采用阻抗+上室結(jié)構(gòu)形式,設(shè)置在尾水岔管下游20 m處。水泵水輪機(jī)主要參數(shù)見表1。

    1.2 測(cè)試簡(jiǎn)況

    每臺(tái)機(jī)組啟動(dòng)調(diào)試期間均進(jìn)行了甩負(fù)荷試驗(yàn),通常甩負(fù)荷試驗(yàn)重點(diǎn)關(guān)注蝸殼進(jìn)口壓力、尾水錐管進(jìn)口壓力和機(jī)組轉(zhuǎn)速,同時(shí)每臺(tái)機(jī)組甩負(fù)荷時(shí)亦對(duì)無葉區(qū)壓力進(jìn)行了同步測(cè)量,對(duì)有功功率、發(fā)電機(jī)出口開關(guān)位置、導(dǎo)葉關(guān)閉規(guī)律等也進(jìn)行了同步測(cè)量。測(cè)試點(diǎn)滿足國(guó)標(biāo)的相關(guān)要求[16] 。無葉區(qū)壓力測(cè)點(diǎn)位置及相應(yīng)照片見圖1。其中,HC01 為無葉區(qū)壓力測(cè)點(diǎn),HC02 為頂蓋下轉(zhuǎn)輪外緣側(cè)壓力測(cè)點(diǎn)。考慮采用長(zhǎng)測(cè)壓管路導(dǎo)致的脈動(dòng)壓力放大效應(yīng)[17] ,壓力測(cè)點(diǎn)采用就近布置方式。導(dǎo)葉關(guān)閉規(guī)律和機(jī)組轉(zhuǎn)速取自調(diào)速器電氣控制柜模擬量輸出面板, 壓力傳感器采用GE 公司的PTX5072 型傳感器,誤差為±0.2%,頻響范圍為0 ~ 5kHz(-3 dB);數(shù)據(jù)采集系統(tǒng)為HBM 公司的QuantumXMX840A-P,模數(shù)轉(zhuǎn)換24 位,采樣率為1 200 Hz。

    6 臺(tái)機(jī)組甩額定負(fù)荷試驗(yàn)含甩25%、50%、75%、100%額定負(fù)荷,其中甩100%額定負(fù)荷時(shí)無葉區(qū)壓力最大,相應(yīng)甩負(fù)荷時(shí)參數(shù)見表2,實(shí)測(cè)機(jī)組轉(zhuǎn)速與導(dǎo)葉關(guān)閉規(guī)律見圖2(為便于比對(duì),導(dǎo)葉開度和機(jī)組轉(zhuǎn)速分別以額定開度、額定轉(zhuǎn)速為基準(zhǔn)折算成相對(duì)值),實(shí)測(cè)無葉區(qū)壓力見圖3。在圖2 和圖3 中,為便于對(duì)比分析,將6 臺(tái)機(jī)組實(shí)際甩負(fù)荷的時(shí)間統(tǒng)一調(diào)整至發(fā)電機(jī)出口開關(guān)斷開時(shí)間,且斷開前保留10 s 數(shù)據(jù)、斷開后保留90 s 數(shù)據(jù),一共100 s 數(shù)據(jù)。

    由表2 和圖2、圖3 可見,2U 和3U 在水頭相當(dāng)?shù)那闆r下,2U 導(dǎo)葉開度明顯偏大,無葉區(qū)壓力脈動(dòng)也較其他機(jī)組偏大。主要原因是:2U 在帶負(fù)荷過程中控制環(huán)與其中一個(gè)活動(dòng)導(dǎo)葉之間的連板脫落,在導(dǎo)葉自關(guān)閉特性下導(dǎo)致機(jī)組帶同樣負(fù)荷時(shí)正常導(dǎo)葉需要更大的開度;同時(shí),該導(dǎo)葉關(guān)閉導(dǎo)致無葉區(qū)流態(tài)對(duì)稱性受到影響,引起無葉區(qū)壓力脈動(dòng)偏大。因此,2U 無葉區(qū)壓力脈動(dòng)數(shù)據(jù)雖不代表正常甩負(fù)荷情況下的脈動(dòng)水平,但反映了異常情況下機(jī)組脈動(dòng)情況,作為運(yùn)行異常數(shù)據(jù)具有極其難得的參考價(jià)值。

    2 時(shí)域特征分析

    甩負(fù)荷時(shí)無葉區(qū)壓力脈動(dòng)由水擊壓力疊加脈動(dòng)壓力構(gòu)成。水擊壓力的波動(dòng)幅值取決于導(dǎo)葉關(guān)閉規(guī)律,波動(dòng)周期取決于引水管路特性,屬于低頻成分;脈動(dòng)壓力幅值受非黏滯性勢(shì)流干涉和黏滯性尾跡干涉[18] ,屬于高頻成分。因此,通過合適的濾波器能夠?qū)⑺畵魤毫兔}動(dòng)壓力分離。很多方法可以用于提取信號(hào)中的低頻成分,如低通濾波器、局域均值分解、經(jīng)驗(yàn)?zāi)B(tài)分解、Savitzky-Golay 濾波器等。局域均值分解與經(jīng)驗(yàn)?zāi)B(tài)分解方法存在模態(tài)混疊問題;低通濾波器種類多且編程實(shí)現(xiàn)略復(fù)雜;Savitzky-Golay 濾波器為線性濾波器,算法較低通濾波器易于實(shí)現(xiàn),且該方法已在水力機(jī)械壓力信號(hào)處理方面得到應(yīng)用[14,19] ,因此本研究采用Savitzky-Golay 濾波器對(duì)壓力數(shù)據(jù)進(jìn)行低通濾波??紤]到無葉區(qū)壓力信號(hào)中低頻成分主要是水擊壓力,水擊壓力波動(dòng)周期通常在數(shù)秒至數(shù)十秒,低通濾波器的截止頻率設(shè)置為0.1 倍轉(zhuǎn)頻即可實(shí)現(xiàn)低頻成分的提取[20] ;同時(shí),將低頻成分從原始信號(hào)中扣除作為脈動(dòng)壓力的高頻成分。壓力信號(hào)處理流程見圖4。

    采用Savitzky-Golay 濾波器對(duì)6 臺(tái)機(jī)組無葉區(qū)壓力進(jìn)行處理,處理結(jié)果見圖5。由圖5(a)可見,6 臺(tái)機(jī)組甩負(fù)荷后無葉區(qū)壓力低頻趨勢(shì)基本一致,表現(xiàn)出了良好一致性。圖5(a)中趨勢(shì)壓力波峰數(shù)量與圖2(a)中轉(zhuǎn)速波峰數(shù)量一致,主要原因是:受限于S 區(qū)特性及初始導(dǎo)葉開度(取決于試驗(yàn)時(shí)水頭),當(dāng)水頭較低時(shí)導(dǎo)葉開度較大,導(dǎo)葉關(guān)閉時(shí)間較長(zhǎng),造成轉(zhuǎn)輪在S 區(qū)來回穿插,導(dǎo)致較高水頭時(shí)波峰數(shù)量多1 次;2U 由于一個(gè)導(dǎo)葉失控關(guān)閉導(dǎo)致開度偏大,因此較正常運(yùn)行開度偏大,亦造成關(guān)閉時(shí)間變長(zhǎng)。

    圖5(b)中機(jī)組甩負(fù)荷過程脈動(dòng)壓力屬于典型的非穩(wěn)態(tài)信號(hào),不能直接計(jì)算時(shí)域特征值。參考文獻(xiàn)[20]針對(duì)甩負(fù)荷時(shí)尾水錐管脈動(dòng)壓力處理方式,采用4 個(gè)旋轉(zhuǎn)周期窗口對(duì)應(yīng)數(shù)據(jù),步長(zhǎng)?。?個(gè)旋轉(zhuǎn)周期,對(duì)應(yīng)數(shù)據(jù)對(duì)脈動(dòng)壓力數(shù)據(jù)進(jìn)行滑移計(jì)算,以獲得整個(gè)過程中脈動(dòng)壓力峰峰值的變化規(guī)律。考慮到無葉區(qū)趨勢(shì)壓力、脈動(dòng)壓力與機(jī)組轉(zhuǎn)速相關(guān),圖6 分別給出了無葉區(qū)脈動(dòng)壓力峰峰值、趨勢(shì)壓力與轉(zhuǎn)速關(guān)系對(duì)比曲線,表3 給出了趨勢(shì)壓力、脈動(dòng)壓力與轉(zhuǎn)速的相關(guān)系數(shù),表4給出了脈動(dòng)壓力峰峰值與相對(duì)值。

    由圖6(a)可知6 臺(tái)機(jī)組脈動(dòng)壓力峰峰值變化趨勢(shì)基本一致,且與機(jī)組轉(zhuǎn)速正相關(guān)(表3 可見脈動(dòng)壓力峰峰值與機(jī)組轉(zhuǎn)速相關(guān)系數(shù)為0.67~0.76),存在兩次顯著的脈動(dòng)壓力峰峰值極大值,分別對(duì)應(yīng)前兩次機(jī)組轉(zhuǎn)速極大值點(diǎn),在機(jī)組轉(zhuǎn)速達(dá)到最大值的第一個(gè)升速過程中,脈動(dòng)壓力達(dá)到最大值。為說明脈動(dòng)壓力峰峰值之間的一致性,計(jì)算6 臺(tái)機(jī)組之間的脈動(dòng)壓力峰峰值相關(guān)系數(shù),計(jì)算結(jié)果見表5,從表中可見:最大值為0.984 3(3U 與6U),最小值為0.866 2(2U 與6U),相關(guān)系數(shù)水平處于高水平,進(jìn)一步驗(yàn)證了脈動(dòng)壓力變化趨勢(shì)的一致性。絕對(duì)值方面,機(jī)組之間存在明顯差異,2U 脈動(dòng)壓力峰峰值在甩負(fù)荷后的不同時(shí)刻均處于6 臺(tái)機(jī)組中較高水平,表明單個(gè)導(dǎo)葉拒動(dòng)明顯增大了無葉區(qū)脈動(dòng)壓力水平。由表4 可見,6 臺(tái)機(jī)組脈動(dòng)壓力最大峰峰值均處于非常高的水平,較正常帶額定負(fù)荷運(yùn)行時(shí)脈動(dòng)壓力峰峰值顯著增大,其中1U、2U、6U脈動(dòng)壓力峰峰值的絕對(duì)值均超過6 000 kPa 并超過了毛水頭,說明在甩額定負(fù)荷時(shí)無葉區(qū)存在極為強(qiáng)烈的脈動(dòng)壓力。

    無葉區(qū)大幅脈動(dòng)壓力將對(duì)轉(zhuǎn)輪產(chǎn)生不利影響,是造成轉(zhuǎn)輪失效的主要原因,相關(guān)案例已有報(bào)道[2,21] 。在甩負(fù)荷過程中,由于脈動(dòng)壓力加大,因此對(duì)轉(zhuǎn)輪造成的破壞效果將顯著增強(qiáng)。針對(duì)模型水泵水輪機(jī)甩負(fù)荷時(shí)的脈動(dòng)壓力,文獻(xiàn)[19]指出無葉區(qū)脈動(dòng)壓力在甩負(fù)荷后達(dá)到甩前穩(wěn)態(tài)工況的2 倍;文獻(xiàn)[22]指出,甩負(fù)荷過程中無葉區(qū)脈動(dòng)壓力幅值增加主要原因是轉(zhuǎn)輪進(jìn)口回流的出現(xiàn)和發(fā)展增大了無葉區(qū)內(nèi)流體的湍動(dòng)能,加強(qiáng)了無葉區(qū)內(nèi)的動(dòng)靜干涉現(xiàn)象,最大脈動(dòng)壓力幅值可達(dá)到初始值的5 倍。從實(shí)際測(cè)試結(jié)果來看,模型試驗(yàn)與數(shù)值模擬顯著低估了甩負(fù)荷過程中的脈動(dòng)壓力。由于模型試驗(yàn)、數(shù)值模擬與實(shí)測(cè)結(jié)果有差異,大部分抽水蓄能電站采用一管多機(jī)形式布置,多機(jī)同時(shí)甩負(fù)荷時(shí)無葉區(qū)脈動(dòng)壓力較單機(jī)時(shí)將顯著增大[23] ,因此如果采用模型試驗(yàn)與數(shù)值模擬結(jié)果對(duì)轉(zhuǎn)輪進(jìn)行設(shè)計(jì),將嚴(yán)重高估轉(zhuǎn)輪的預(yù)期壽命,造成轉(zhuǎn)輪過早疲勞損傷并失效。

    圖6(b)可見:無葉區(qū)趨勢(shì)壓力與機(jī)組轉(zhuǎn)速變化趨勢(shì)基本一致,趨勢(shì)壓力變化與機(jī)組轉(zhuǎn)速正相關(guān),兩者之間的相關(guān)系數(shù)均在0.99 左右,說明無葉區(qū)趨勢(shì)壓力變化主要取決于機(jī)組轉(zhuǎn)速,而由導(dǎo)葉快速關(guān)閉產(chǎn)生的水擊壓力影響可忽略,這一觀測(cè)結(jié)果與文獻(xiàn)[19]對(duì)模型機(jī)組無葉區(qū)趨勢(shì)壓力的觀測(cè)結(jié)果一致。根據(jù)轉(zhuǎn)輪入口處的速度三角形,轉(zhuǎn)輪入口處流速由圓周速度和相對(duì)速度確定,考慮到無葉區(qū)趨勢(shì)壓力與機(jī)組轉(zhuǎn)速的強(qiáng)相關(guān)性,可以判斷趨勢(shì)壓力取決于圓周速度而與相對(duì)速度無關(guān)。發(fā)生這一現(xiàn)象的本質(zhì)原因可以分為兩個(gè)時(shí)間段進(jìn)行解釋:在導(dǎo)葉未全關(guān)的時(shí)候,水泵水輪機(jī)轉(zhuǎn)輪流道狹長(zhǎng),具有明顯的S 區(qū);在全特性曲線上,正常滿負(fù)荷運(yùn)行時(shí)機(jī)組運(yùn)行工況點(diǎn)位置較靠近S 區(qū),機(jī)組一旦甩負(fù)荷后將迅速進(jìn)入S 區(qū),并往返循環(huán)多次[20] ;而在S 區(qū)內(nèi),存在流動(dòng)不穩(wěn)定導(dǎo)致的旋轉(zhuǎn)失速現(xiàn)象[4] ,轉(zhuǎn)輪流道堵塞,同時(shí)在無葉區(qū)內(nèi)靠近活動(dòng)導(dǎo)葉出口處存在明顯的旋渦與回流從而形成水環(huán)[24] ,進(jìn)一步堵塞流道,引起壓力上升。在導(dǎo)葉全關(guān)后,無葉區(qū)內(nèi)流道完全堵塞,造成壓力上升。

    3 結(jié)論

    本研究以某抽水蓄能電站6 臺(tái)機(jī)組單機(jī)甩額定負(fù)荷時(shí)無葉區(qū)脈動(dòng)壓力信號(hào)為分析對(duì)象,基于Savitzky-Golay 濾波器分離了趨勢(shì)壓力與脈動(dòng)壓力,對(duì)影響壓力上升及脈動(dòng)變化水平的因素進(jìn)行了分析。通過研究獲得了以下結(jié)論:

    1)甩負(fù)荷過程中無葉區(qū)趨勢(shì)壓力與機(jī)組轉(zhuǎn)速相關(guān)系數(shù)在0.99 左右,表明壓力升高主要受機(jī)組轉(zhuǎn)速影響,壓力上升主要由流道堵塞效應(yīng)導(dǎo)致;脈動(dòng)壓力峰峰值與機(jī)組轉(zhuǎn)速相關(guān)系數(shù)為0.67 ~0.76,表明脈動(dòng)壓力峰峰值與機(jī)組轉(zhuǎn)速正相關(guān)。

    2)甩負(fù)荷過程中無葉區(qū)存在強(qiáng)烈的由動(dòng)靜干涉、旋轉(zhuǎn)失速等導(dǎo)致的脈動(dòng)壓力,且脈動(dòng)壓力峰峰值超過了毛水頭,為正常額定負(fù)荷運(yùn)行的7~22 倍,顯著高于數(shù)值模擬結(jié)果。

    參考文獻(xiàn):

    [1] ZHANG F,LOWYS P Y,HOUDELINE J B,et al.Pump?Tur?bine Rotor?Stator Interaction Induced Vibration: ProblemResolution and Experience[C] //30th IAHR Symposium onHydraulic Machinery and Systems. IOP Conf. Series: Earthand Environmental Science,2021,774(1):012124.

    [2] 龐希斌,彭碩群,祝加勇,等.水泵水輪機(jī)轉(zhuǎn)輪裂紋成因分析及處理[J].水電站機(jī)電技術(shù),2019,42(2):51-54.

    [3] ZUO Z,LIU S. Flow?Induced Instabilities in Pump?Turbinesin China[J]. Engineering,2017, 3(4):504-511.

    [4] HASMATUCHI V, FARHAT M, ROTH S, et al. ExperimentalEvidence of Rotating Stall in a Pump?Turbine at Off?DesignConditions in GENERATING mode[J]. Journal of Fluid En?gineering,2011,133:051104.

    [5] RAN H,LUO X.Experimental Study of Instability Character?istics in Pump Turbines[J].Journal of Hydraulic Research,2018,56(6),871-876.

    [6] 李德友.水泵水輪機(jī)駝峰區(qū)流動(dòng)機(jī)理及瞬態(tài)特性研究[D].哈爾濱:哈爾濱工業(yè)大學(xué),2017:101-109.

    [7] ZHANG W,CHEN Z,ZHU B, et al.Pressure Fluctuation andFlow Instability in S?Shaped Region of a Reversible Pump?Turbine[J]. Renewable Energy,2020,154:826-840.

    [8] ZHANG X, CHENG Y, XIA L,et al.Looping Dynamic Charac?teristics of a Pump?Turbine in the S?Shaped Region DuringRunaway [ J ]. Journal of Fluids Engineering, 2016,138:091102.

    [9] YIN C, ZENG W, YANG J.Transient Simulation and Analysisof the Simultaneous Load Rejection Process in Pumped StoragePower Stations Using a 1-D- 3-D Coupling Method[J].Journal of Hydrodynamics,2021,33(5):979-991.

    [10] 張曉曦,陳秋華.水泵水輪機(jī)甩負(fù)荷過渡過程尾水管水柱分離數(shù)值模擬[J].水動(dòng)力學(xué)研究與進(jìn)展,2019,34(6):749-754.

    [11] TRIVEDI C,AGNALT E,DAHLHAUG O G.ExperimentalInvestigation of a Francis Turbine During Exigent Rampingand Transition into Total Load Rejection[J]. Journal of Hy?draulic Engineering,2018,144(6):04018027.

    [12] LIU D,XU W,ZHAO Y,et al. Experimental Study of theFlow Field of a High Head Model Pump Turbine Based onPIV Technique[J]. Journal of Hydrodynamics,2021,33(5):1045-1055.

    [13] 李立,陳源,伍志軍,等.基于小波分解的抽水蓄能機(jī)組甩負(fù)荷反演預(yù)測(cè)[J]. 人民黃河,2021,43(9):144 -149,155.

    [14] 楊建東,胡金弘,曾威,等.原型混流式水泵水輪機(jī)過渡過程中的壓力脈動(dòng)[J].水利學(xué)報(bào),2016,47(7):858-864.

    [15] NILSSON O,SJELVGREN D. Hydro Unit Start?Up Costsand Their Impact on the Short Term Scheduling Strategies ofSwedish Power Procedures [ J]. IEEE Transactions onPower Systems,1997,12(1):38-44.

    [16] 中華人民共和國(guó)國(guó)家質(zhì)量監(jiān)督檢驗(yàn)檢疫總局,中國(guó)國(guó)家標(biāo)準(zhǔn)化管理委員會(huì).水力機(jī)械(水輪機(jī)、蓄能泵和水泵水輪機(jī)) 振動(dòng)和脈動(dòng)現(xiàn)場(chǎng)測(cè)試規(guī)程:GB/ T 17189—2017[S].北京:中國(guó)標(biāo)準(zhǔn)出版社,2018:13-15.

    [17] 張飛,郭磊,宮讓勤,等.基于線性摩擦模型的水力機(jī)械測(cè)壓管路特性[J].工程熱物理學(xué)報(bào),2018,39(8):1725-1730.

    [18] DRING R P,JOSLYN H D,HARDIN L W,et al.TurbineRotor?Stator Interaction[J].Journal of Engineering for GasTurbines and Power,1982,104(4):729-742.

    [19] CHIRAG T,MICHEL C, BHUPENDRA G,et al.TransientPressure Measurements on a High Head Model Francis Tur?bine During Emergency Shutdown,Total Load Rejection,and Runaway [J]. ASME Journal of Fluids Engineering,2014,136(12):121107.

    [20] 張飛,王憲平.抽水蓄能機(jī)組甩負(fù)荷試驗(yàn)時(shí)尾水錐管壓力[J].農(nóng)業(yè)工程學(xué)報(bào),2020,36(20):93-101.[21] EDUARD E,CARME V, HUANG X ,et al.Failure Investi?gation of a Large Pump?Turbine Runner[J]. EngineeringFailure Analysis,2012,23:27-34.

    [22] 周勤,夏林生,張春澤,等.水泵水輪機(jī)甩負(fù)荷過渡過程中的壓力脈動(dòng)和轉(zhuǎn)輪受力[J].水利學(xué)報(bào),2018,49(11):1429-1438.

    [23] HU J , YANG J ,ZENG W,et al.Effect of the Speed Factoron the Amplitude of the Blade Passing Frequency in theVaneless Space of a Pump Turbine in Turbine Mode[J].Journal of Fluids Engineering,2021,143:111203.

    [24] 朱偉,肖業(yè)祥,姚洋陽(yáng),等.混流式水泵水輪機(jī)小開度S特性區(qū)內(nèi)流特性分析[J].水力發(fā)電學(xué)報(bào),2015,34(10):138-144.

    【責(zé)任編輯 簡(jiǎn) 群】

    日日撸夜夜添| 99久久九九国产精品国产免费| 欧美中文日本在线观看视频| 免费大片18禁| 精华霜和精华液先用哪个| 国语自产精品视频在线第100页| 亚洲av电影不卡..在线观看| 亚洲欧美日韩卡通动漫| 五月伊人婷婷丁香| av福利片在线观看| 免费观看精品视频网站| 天堂网av新在线| 三级毛片av免费| 午夜亚洲福利在线播放| 一区福利在线观看| 久久香蕉精品热| 观看免费一级毛片| 国产精品美女特级片免费视频播放器| 亚洲狠狠婷婷综合久久图片| 亚洲成人中文字幕在线播放| 在线观看66精品国产| 亚洲三级黄色毛片| 日韩 亚洲 欧美在线| 观看美女的网站| 久久精品国产清高在天天线| 在线观看免费视频日本深夜| 免费在线观看日本一区| 成年免费大片在线观看| 免费搜索国产男女视频| 欧美色视频一区免费| 无人区码免费观看不卡| 12—13女人毛片做爰片一| 内地一区二区视频在线| 天美传媒精品一区二区| 日韩欧美一区二区三区在线观看| 日本黄大片高清| 尾随美女入室| 亚洲图色成人| 午夜精品久久久久久毛片777| 国产精品精品国产色婷婷| 国产伦精品一区二区三区视频9| 夜夜看夜夜爽夜夜摸| 国产亚洲精品综合一区在线观看| 欧美成人免费av一区二区三区| 俺也久久电影网| 亚洲最大成人av| 热99在线观看视频| 午夜福利在线在线| 中文亚洲av片在线观看爽| 麻豆一二三区av精品| 欧洲精品卡2卡3卡4卡5卡区| 真人一进一出gif抽搐免费| 岛国在线免费视频观看| 精品乱码久久久久久99久播| 午夜福利在线观看吧| 九九久久精品国产亚洲av麻豆| 一区二区三区高清视频在线| 日韩国内少妇激情av| 久久人妻av系列| 欧美3d第一页| 久久九九热精品免费| 欧美一区二区亚洲| 99久久中文字幕三级久久日本| 91午夜精品亚洲一区二区三区 | 久久久久久九九精品二区国产| 欧美绝顶高潮抽搐喷水| 99热这里只有是精品50| 99热这里只有是精品在线观看| www.www免费av| 欧美性感艳星| 一区福利在线观看| 国产精品99久久久久久久久| 给我免费播放毛片高清在线观看| 久久久久久久午夜电影| 男女做爰动态图高潮gif福利片| 国产高清视频在线播放一区| 中文资源天堂在线| av在线老鸭窝| 男人狂女人下面高潮的视频| 精品欧美国产一区二区三| 精品免费久久久久久久清纯| 国产亚洲91精品色在线| 精品久久久久久久久亚洲 | 精品久久久久久久末码| 狂野欧美激情性xxxx在线观看| 国产免费av片在线观看野外av| 国产av在哪里看| 永久网站在线| 成年版毛片免费区| 亚洲成av人片在线播放无| 级片在线观看| 老司机午夜福利在线观看视频| 国产淫片久久久久久久久| 99久久无色码亚洲精品果冻| av视频在线观看入口| 亚洲精品色激情综合| 色噜噜av男人的天堂激情| 亚洲av中文av极速乱 | av视频在线观看入口| 淫秽高清视频在线观看| 亚洲精品色激情综合| 一区福利在线观看| 亚洲精品一区av在线观看| av在线蜜桃| 国产精品1区2区在线观看.| 在线国产一区二区在线| 男女那种视频在线观看| 人妻久久中文字幕网| 一级黄片播放器| 国产欧美日韩一区二区精品| 性插视频无遮挡在线免费观看| 黄色视频,在线免费观看| 一卡2卡三卡四卡精品乱码亚洲| 国产精品久久电影中文字幕| 麻豆av噜噜一区二区三区| 亚洲欧美日韩高清专用| 一边摸一边抽搐一进一小说| 校园人妻丝袜中文字幕| a级毛片免费高清观看在线播放| 天堂动漫精品| 日本a在线网址| 久久精品综合一区二区三区| 99热只有精品国产| 久久精品国产亚洲av香蕉五月| 免费看a级黄色片| 老熟妇仑乱视频hdxx| 一区二区三区激情视频| 国产毛片a区久久久久| 在线观看舔阴道视频| 国产精品电影一区二区三区| 观看免费一级毛片| 高清在线国产一区| 欧美日韩亚洲国产一区二区在线观看| 一级黄片播放器| 国产成人aa在线观看| 给我免费播放毛片高清在线观看| 午夜福利在线观看吧| 91久久精品国产一区二区成人| 99九九线精品视频在线观看视频| 亚洲精品粉嫩美女一区| 精品午夜福利在线看| 一进一出好大好爽视频| 国产一区二区在线观看日韩| 亚洲自拍偷在线| 他把我摸到了高潮在线观看| 人人妻,人人澡人人爽秒播| 国产女主播在线喷水免费视频网站 | 国产成人a区在线观看| 久久久久久久久久久丰满 | 中国美白少妇内射xxxbb| 99久久精品一区二区三区| 亚洲欧美日韩高清在线视频| 精品人妻一区二区三区麻豆 | 搡老妇女老女人老熟妇| 国产精品精品国产色婷婷| 狂野欧美激情性xxxx在线观看| 琪琪午夜伦伦电影理论片6080| 天堂√8在线中文| 亚洲18禁久久av| 国产精品人妻久久久影院| 亚洲综合色惰| 丰满的人妻完整版| 亚洲综合色惰| 91麻豆精品激情在线观看国产| www.色视频.com| 丰满的人妻完整版| 免费av毛片视频| 亚洲综合色惰| 色吧在线观看| 国产私拍福利视频在线观看| 香蕉av资源在线| 高清日韩中文字幕在线| 高清日韩中文字幕在线| 国产精品,欧美在线| 一进一出好大好爽视频| 免费av毛片视频| 99久久中文字幕三级久久日本| 高清日韩中文字幕在线| 免费电影在线观看免费观看| 久久精品国产清高在天天线| 色吧在线观看| 欧美黑人欧美精品刺激| 国产国拍精品亚洲av在线观看| 国产亚洲91精品色在线| 九色成人免费人妻av| 床上黄色一级片| 国国产精品蜜臀av免费| 国产精品爽爽va在线观看网站| 欧美性猛交╳xxx乱大交人| 久久久久久大精品| 国产精品精品国产色婷婷| 俄罗斯特黄特色一大片| 国内毛片毛片毛片毛片毛片| 男女之事视频高清在线观看| 岛国在线免费视频观看| 毛片一级片免费看久久久久 | 久久草成人影院| 色av中文字幕| av黄色大香蕉| 久久久久久久亚洲中文字幕| 精品国内亚洲2022精品成人| 人人妻人人看人人澡| 国产免费男女视频| 久久久久久久久久成人| 国产探花极品一区二区| 亚洲色图av天堂| 此物有八面人人有两片| 波多野结衣高清无吗| 国产午夜精品久久久久久一区二区三区 | 禁无遮挡网站| 简卡轻食公司| 韩国av在线不卡| 91麻豆av在线| 精品久久久久久久久久久久久| 亚洲人成网站高清观看| 日韩一区二区视频免费看| 在线播放无遮挡| 亚洲国产精品合色在线| 国产精品亚洲一级av第二区| 身体一侧抽搐| av中文乱码字幕在线| 国产一区二区在线观看日韩| 丰满的人妻完整版| АⅤ资源中文在线天堂| 我要看日韩黄色一级片| 一级毛片久久久久久久久女| 毛片一级片免费看久久久久 | 欧美高清性xxxxhd video| 真实男女啪啪啪动态图| 国产伦精品一区二区三区四那| av在线亚洲专区| 国产精品98久久久久久宅男小说| 精品一区二区三区视频在线观看免费| 成人国产一区最新在线观看| 国产真实伦视频高清在线观看 | 成年免费大片在线观看| 亚洲不卡免费看| 啦啦啦观看免费观看视频高清| 夜夜夜夜夜久久久久| 天堂影院成人在线观看| 欧美+日韩+精品| 亚洲人成网站在线播放欧美日韩| 日韩欧美国产一区二区入口| 99热网站在线观看| 制服丝袜大香蕉在线| 别揉我奶头~嗯~啊~动态视频| 亚洲男人的天堂狠狠| 亚洲精华国产精华液的使用体验 | 午夜精品在线福利| 久99久视频精品免费| 最近在线观看免费完整版| 国产欧美日韩一区二区精品| 久久久国产成人免费| 内地一区二区视频在线| 中文字幕av在线有码专区| 国产精品1区2区在线观看.| 午夜免费成人在线视频| 国产精品国产三级国产av玫瑰| 一级毛片久久久久久久久女| 男人和女人高潮做爰伦理| 亚洲精品456在线播放app | 国产成年人精品一区二区| 日日干狠狠操夜夜爽| 国产美女午夜福利| av黄色大香蕉| 欧美激情国产日韩精品一区| 我的老师免费观看完整版| 99久久精品热视频| 我要搜黄色片| 一进一出抽搐动态| 中文字幕人妻熟人妻熟丝袜美| 国产伦一二天堂av在线观看| 一区二区三区高清视频在线| 一进一出好大好爽视频| 国产高清有码在线观看视频| 三级毛片av免费| 亚洲av成人精品一区久久| 国产精品98久久久久久宅男小说| 国产精品无大码| 亚洲内射少妇av| 亚洲精品久久国产高清桃花| 性欧美人与动物交配| 99热这里只有是精品50| 亚洲av五月六月丁香网| 亚洲精品一区av在线观看| 五月玫瑰六月丁香| 免费搜索国产男女视频| 97超级碰碰碰精品色视频在线观看| 又黄又爽又刺激的免费视频.| 又黄又爽又免费观看的视频| 亚洲精品亚洲一区二区| 色吧在线观看| 免费人成视频x8x8入口观看| 中文字幕精品亚洲无线码一区| 亚洲天堂国产精品一区在线| 亚洲最大成人中文| 人人妻,人人澡人人爽秒播| 波多野结衣高清作品| 中亚洲国语对白在线视频| 国产单亲对白刺激| 欧美高清成人免费视频www| 麻豆久久精品国产亚洲av| 久久精品国产清高在天天线| 一区二区三区四区激情视频 | 日韩欧美精品v在线| a级毛片免费高清观看在线播放| 免费观看的影片在线观看| 精品久久久久久成人av| 亚洲va日本ⅴa欧美va伊人久久| 久久欧美精品欧美久久欧美| 床上黄色一级片| 国产精品嫩草影院av在线观看 | 亚洲经典国产精华液单| 精品久久久久久久久亚洲 | 91久久精品国产一区二区成人| 干丝袜人妻中文字幕| 国产在线男女| 老熟妇仑乱视频hdxx| 国产高清视频在线观看网站| www.www免费av| 欧美国产日韩亚洲一区| 人妻制服诱惑在线中文字幕| 亚洲乱码一区二区免费版| 亚洲av日韩精品久久久久久密| 午夜爱爱视频在线播放| 日韩大尺度精品在线看网址| 欧美区成人在线视频| 黄色女人牲交| 亚洲精品影视一区二区三区av| 亚洲av电影不卡..在线观看| 少妇熟女aⅴ在线视频| 国产精品无大码| 超碰av人人做人人爽久久| aaaaa片日本免费| 中文资源天堂在线| av在线观看视频网站免费| 中文字幕熟女人妻在线| 非洲黑人性xxxx精品又粗又长| 国产在线男女| 嫩草影视91久久| 亚洲av日韩精品久久久久久密| 可以在线观看的亚洲视频| 天堂网av新在线| 免费在线观看影片大全网站| 国产成年人精品一区二区| 春色校园在线视频观看| 久久人妻av系列| 国产午夜福利久久久久久| 色哟哟·www| 国产精品伦人一区二区| 中文亚洲av片在线观看爽| av视频在线观看入口| 真人做人爱边吃奶动态| 深夜精品福利| 男人狂女人下面高潮的视频| 动漫黄色视频在线观看| 亚洲黑人精品在线| 国产伦一二天堂av在线观看| 国产精品一区二区三区四区久久| 校园春色视频在线观看| 日韩人妻高清精品专区| 99riav亚洲国产免费| 伊人久久精品亚洲午夜| 哪里可以看免费的av片| 中文字幕人妻熟人妻熟丝袜美| 亚洲美女黄片视频| 久99久视频精品免费| 日本成人三级电影网站| 欧美性猛交╳xxx乱大交人| 国产高清三级在线| 色噜噜av男人的天堂激情| 精品人妻熟女av久视频| 精品福利观看| 国产精品爽爽va在线观看网站| 亚洲avbb在线观看| 亚洲成人免费电影在线观看| 久久久久久久久大av| 午夜精品久久久久久毛片777| 国产在线精品亚洲第一网站| 一级a爱片免费观看的视频| 麻豆国产av国片精品| 小说图片视频综合网站| а√天堂www在线а√下载| 18+在线观看网站| 国产精品,欧美在线| 久久精品国产99精品国产亚洲性色| 精品久久久久久久久久免费视频| 91午夜精品亚洲一区二区三区 | 在线国产一区二区在线| 婷婷六月久久综合丁香| 麻豆成人午夜福利视频| 亚洲中文日韩欧美视频| 亚洲国产精品sss在线观看| www.www免费av| 国产成人一区二区在线| 在线观看一区二区三区| 久久天躁狠狠躁夜夜2o2o| 国产精品99久久久久久久久| 成年女人永久免费观看视频| 亚洲人成网站在线播| 看免费成人av毛片| 亚洲18禁久久av| 蜜桃亚洲精品一区二区三区| 成人特级黄色片久久久久久久| 一夜夜www| 久久久久久久久久久丰满 | 精品午夜福利在线看| 亚洲欧美激情综合另类| 日韩欧美 国产精品| 偷拍熟女少妇极品色| 99热这里只有精品一区| av在线蜜桃| 日韩欧美精品免费久久| 人人妻人人看人人澡| 搡老岳熟女国产| 午夜精品久久久久久毛片777| 欧美人与善性xxx| 91av网一区二区| 久久午夜亚洲精品久久| 看片在线看免费视频| 国内揄拍国产精品人妻在线| 亚洲欧美日韩卡通动漫| 精品久久久久久久久久免费视频| 99久久久亚洲精品蜜臀av| 亚洲无线在线观看| 精品久久久噜噜| 亚洲综合色惰| 蜜桃亚洲精品一区二区三区| 国产极品精品免费视频能看的| 男女之事视频高清在线观看| 欧美潮喷喷水| 午夜福利欧美成人| 麻豆成人午夜福利视频| av福利片在线观看| 国产一区二区激情短视频| 精品午夜福利视频在线观看一区| 国产精品一区二区性色av| 中国美女看黄片| 午夜免费男女啪啪视频观看 | 自拍偷自拍亚洲精品老妇| 亚洲欧美日韩高清专用| 亚洲av五月六月丁香网| 美女黄网站色视频| 国产日本99.免费观看| 欧美日韩黄片免| 国产精品一区二区性色av| 一级毛片久久久久久久久女| 国产精品,欧美在线| 亚洲欧美清纯卡通| 亚洲国产日韩欧美精品在线观看| 男女做爰动态图高潮gif福利片| 国产真实乱freesex| 欧美日韩瑟瑟在线播放| 长腿黑丝高跟| 99久国产av精品| 欧美日韩黄片免| 午夜亚洲福利在线播放| 国产精品野战在线观看| 国产伦精品一区二区三区视频9| 韩国av在线不卡| 大型黄色视频在线免费观看| 99久久精品国产国产毛片| 亚洲av二区三区四区| a级毛片a级免费在线| 久久精品国产亚洲av天美| 国产乱人视频| 一个人观看的视频www高清免费观看| 99热精品在线国产| 一区福利在线观看| 99热6这里只有精品| 一区二区三区激情视频| 国产伦人伦偷精品视频| 免费在线观看日本一区| 午夜a级毛片| 一级黄片播放器| 一级黄色大片毛片| 俄罗斯特黄特色一大片| 精品一区二区三区av网在线观看| 亚洲成人中文字幕在线播放| 久久久国产成人免费| 精品福利观看| 我要看日韩黄色一级片| 国产欧美日韩精品一区二区| 88av欧美| 人妻丰满熟妇av一区二区三区| 国产私拍福利视频在线观看| 午夜精品久久久久久毛片777| 熟女电影av网| 亚洲欧美日韩卡通动漫| 搡老熟女国产l中国老女人| 国产精品一区二区三区四区久久| 人妻久久中文字幕网| 亚洲精品亚洲一区二区| 高清毛片免费观看视频网站| 国产黄a三级三级三级人| 中文资源天堂在线| 搡女人真爽免费视频火全软件 | videossex国产| 蜜桃亚洲精品一区二区三区| 国产极品精品免费视频能看的| 一级毛片久久久久久久久女| 国产在线精品亚洲第一网站| 欧美日韩综合久久久久久 | 久久精品国产亚洲av香蕉五月| 嫁个100分男人电影在线观看| 一个人看视频在线观看www免费| 深夜精品福利| 免费黄网站久久成人精品| 中文亚洲av片在线观看爽| 在线免费观看的www视频| 两人在一起打扑克的视频| 老熟妇仑乱视频hdxx| 一个人观看的视频www高清免费观看| 欧美绝顶高潮抽搐喷水| 国产高清视频在线观看网站| 啦啦啦啦在线视频资源| 免费观看在线日韩| 日韩欧美精品免费久久| 亚洲午夜理论影院| 亚洲综合色惰| 国产中年淑女户外野战色| 国产精品国产高清国产av| 久久精品综合一区二区三区| 日本三级黄在线观看| av天堂在线播放| 校园人妻丝袜中文字幕| 99在线人妻在线中文字幕| 尤物成人国产欧美一区二区三区| 国产一区二区三区av在线 | 国产三级中文精品| 别揉我奶头~嗯~啊~动态视频| 成人特级黄色片久久久久久久| 亚洲av熟女| 特级一级黄色大片| www.www免费av| 国产三级中文精品| 精品久久久久久久久久久久久| 床上黄色一级片| 亚洲最大成人中文| 好男人在线观看高清免费视频| 亚洲欧美日韩无卡精品| 九色国产91popny在线| 黄色视频,在线免费观看| 精品不卡国产一区二区三区| 亚洲人成伊人成综合网2020| 色综合亚洲欧美另类图片| 波多野结衣巨乳人妻| 亚洲aⅴ乱码一区二区在线播放| 999久久久精品免费观看国产| 性欧美人与动物交配| 免费不卡的大黄色大毛片视频在线观看 | 国产黄色小视频在线观看| 在线天堂最新版资源| 桃色一区二区三区在线观看| 久久精品国产清高在天天线| 国产精品亚洲美女久久久| 亚洲最大成人手机在线| 人人妻人人澡欧美一区二区| 午夜久久久久精精品| 黄色视频,在线免费观看| 成人特级黄色片久久久久久久| 久久久国产成人精品二区| 深夜精品福利| 国内精品一区二区在线观看| 国产真实乱freesex| 国产一区二区在线观看日韩| 国产精品电影一区二区三区| 啪啪无遮挡十八禁网站| 欧美+亚洲+日韩+国产| 国内久久婷婷六月综合欲色啪| a级毛片a级免费在线| 国产高清不卡午夜福利| 麻豆精品久久久久久蜜桃| 亚洲av美国av| 亚洲最大成人中文| 18禁黄网站禁片午夜丰满| 草草在线视频免费看| 人人妻人人看人人澡| 日韩高清综合在线| 亚洲成人精品中文字幕电影| 黄色欧美视频在线观看| 一区二区三区四区激情视频 | 啦啦啦观看免费观看视频高清| 91午夜精品亚洲一区二区三区 | 亚洲黑人精品在线| 亚洲无线观看免费| 黄色视频,在线免费观看| 老司机深夜福利视频在线观看| 亚洲七黄色美女视频| 亚洲专区国产一区二区| 午夜精品在线福利| 日日啪夜夜撸| 亚洲美女搞黄在线观看 | 日韩欧美一区二区三区在线观看| 欧美国产日韩亚洲一区| 国产蜜桃级精品一区二区三区| 日本一本二区三区精品| 亚洲熟妇中文字幕五十中出| 窝窝影院91人妻| 色综合站精品国产| 国产免费男女视频| 亚洲五月天丁香| 亚洲中文字幕一区二区三区有码在线看| 桃红色精品国产亚洲av| 国产精品一区二区性色av| 熟女人妻精品中文字幕| 亚洲美女黄片视频| 欧美极品一区二区三区四区| 色哟哟哟哟哟哟| 最近视频中文字幕2019在线8| 国产69精品久久久久777片| 天天一区二区日本电影三级| 欧美性猛交黑人性爽| 国产精品99久久久久久久久| 亚洲黑人精品在线|