• <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) 群】

    国产精品伦人一区二区| 又大又黄又爽视频免费| 亚洲精品久久午夜乱码| 国产视频内射| 国产成人午夜福利电影在线观看| 这个男人来自地球电影免费观看 | 亚洲精品日韩av片在线观看| 在线观看美女被高潮喷水网站| 精品一区二区三区视频在线| 国产淫语在线视频| 男女免费视频国产| 777米奇影视久久| 99久久综合免费| 18禁动态无遮挡网站| 中文资源天堂在线| 只有这里有精品99| 精品一区在线观看国产| 一级黄片播放器| 国产成人精品一,二区| av天堂中文字幕网| 高清日韩中文字幕在线| 午夜福利影视在线免费观看| 国产日韩欧美亚洲二区| 日本-黄色视频高清免费观看| 亚洲av成人精品一二三区| 嫩草影院新地址| 国产精品一二三区在线看| 久热这里只有精品99| 国产免费一级a男人的天堂| 最近最新中文字幕免费大全7| 国产成人a∨麻豆精品| 欧美日韩亚洲高清精品| h日本视频在线播放| 日本色播在线视频| 免费观看无遮挡的男女| 高清毛片免费看| 亚洲国产日韩一区二区| 最新中文字幕久久久久| 亚洲成人手机| 精品亚洲成国产av| 大话2 男鬼变身卡| 一级毛片电影观看| 国产黄色视频一区二区在线观看| 久热这里只有精品99| 22中文网久久字幕| 99九九线精品视频在线观看视频| 超碰97精品在线观看| 中文资源天堂在线| 欧美亚洲 丝袜 人妻 在线| 一级a做视频免费观看| 美女cb高潮喷水在线观看| 新久久久久国产一级毛片| 久久鲁丝午夜福利片| 日本午夜av视频| 欧美一区二区亚洲| 欧美3d第一页| 岛国毛片在线播放| 九九久久精品国产亚洲av麻豆| 校园人妻丝袜中文字幕| 亚洲精品第二区| 成年女人在线观看亚洲视频| 夜夜看夜夜爽夜夜摸| 精品一区二区免费观看| 久久久色成人| 国产又色又爽无遮挡免| 日日啪夜夜撸| 午夜免费鲁丝| 国产深夜福利视频在线观看| 国产成人精品久久久久久| 爱豆传媒免费全集在线观看| 美女cb高潮喷水在线观看| 国产成人a区在线观看| 我的老师免费观看完整版| 午夜激情久久久久久久| 五月天丁香电影| 久久国产精品男人的天堂亚洲 | 国产精品国产三级国产专区5o| 国产又色又爽无遮挡免| 久久99精品国语久久久| 久久6这里有精品| 亚洲精品国产av蜜桃| 国产黄色免费在线视频| 人人妻人人添人人爽欧美一区卜 | 亚洲精品亚洲一区二区| 国产高清国产精品国产三级 | 国产精品av视频在线免费观看| 18禁在线无遮挡免费观看视频| 欧美成人午夜免费资源| 免费观看在线日韩| 97在线视频观看| a 毛片基地| 精品一品国产午夜福利视频| 亚洲av国产av综合av卡| 免费大片18禁| 在线亚洲精品国产二区图片欧美 | 97超碰精品成人国产| 男人舔奶头视频| 男女边吃奶边做爰视频| 久久精品久久久久久久性| 亚洲,欧美,日韩| 国产 一区精品| 久久青草综合色| 嘟嘟电影网在线观看| 99re6热这里在线精品视频| 亚洲美女黄色视频免费看| 少妇裸体淫交视频免费看高清| 亚洲精品456在线播放app| 黄色配什么色好看| 一区二区av电影网| 亚洲精品国产av成人精品| 日韩视频在线欧美| 亚洲一区二区三区欧美精品| 久久99热这里只频精品6学生| 香蕉精品网在线| 黄色一级大片看看| 我的女老师完整版在线观看| 国产人妻一区二区三区在| 久久这里有精品视频免费| 日韩av不卡免费在线播放| 多毛熟女@视频| 亚洲欧美日韩另类电影网站 | 亚洲国产欧美在线一区| 最近2019中文字幕mv第一页| 久久精品国产亚洲av涩爱| 日本午夜av视频| 大片免费播放器 马上看| 精品亚洲成国产av| 日韩人妻高清精品专区| 欧美一区二区亚洲| videos熟女内射| 寂寞人妻少妇视频99o| 精品久久国产蜜桃| 国产精品伦人一区二区| 久久鲁丝午夜福利片| 亚洲综合色惰| 丝袜脚勾引网站| 亚洲精品一二三| 日韩一本色道免费dvd| 成人漫画全彩无遮挡| 特大巨黑吊av在线直播| 满18在线观看网站| 久久精品久久精品一区二区三区| 亚洲av在线观看美女高潮| 中文字幕另类日韩欧美亚洲嫩草| 制服人妻中文乱码| 乱人伦中国视频| 性色av乱码一区二区三区2| 欧美日韩亚洲国产一区二区在线观看 | 丰满饥渴人妻一区二区三| 伊人久久大香线蕉亚洲五| 国产伦理片在线播放av一区| 亚洲黑人精品在线| av天堂久久9| 天天躁夜夜躁狠狠躁躁| 日韩人妻精品一区2区三区| 91老司机精品| 精品亚洲成a人片在线观看| 人成视频在线观看免费观看| 国产av国产精品国产| 亚洲成人免费av在线播放| 永久免费av网站大全| 制服人妻中文乱码| 91成人精品电影| 国产精品 国内视频| 脱女人内裤的视频| 国产精品久久久av美女十八| 婷婷色麻豆天堂久久| 啦啦啦在线免费观看视频4| 亚洲少妇的诱惑av| 久久性视频一级片| 一区二区三区乱码不卡18| 一区福利在线观看| 欧美激情高清一区二区三区| 黄色a级毛片大全视频| 婷婷色av中文字幕| 纵有疾风起免费观看全集完整版| 国产av国产精品国产| 欧美精品高潮呻吟av久久| 18在线观看网站| 美女福利国产在线| 亚洲人成网站在线观看播放| av片东京热男人的天堂| 亚洲精品在线美女| 国产女主播在线喷水免费视频网站| 精品第一国产精品| 国产一区有黄有色的免费视频| 嫁个100分男人电影在线观看 | 大话2 男鬼变身卡| 免费一级毛片在线播放高清视频 | 晚上一个人看的免费电影| 中文字幕人妻丝袜制服| 久久精品久久精品一区二区三区| 亚洲第一av免费看| 69精品国产乱码久久久| 国产精品一区二区精品视频观看| 成人免费观看视频高清| www.熟女人妻精品国产| 女人高潮潮喷娇喘18禁视频| 亚洲人成电影观看| 精品卡一卡二卡四卡免费| 美女午夜性视频免费| 亚洲美女黄色视频免费看| 丰满少妇做爰视频| 午夜福利一区二区在线看| 国产熟女午夜一区二区三区| 亚洲自偷自拍图片 自拍| 国产91精品成人一区二区三区 | 狠狠婷婷综合久久久久久88av| 在线观看国产h片| 99久久人妻综合| 男女下面插进去视频免费观看| 50天的宝宝边吃奶边哭怎么回事| 欧美在线黄色| 久9热在线精品视频| 亚洲欧美清纯卡通| 欧美变态另类bdsm刘玥| 日本vs欧美在线观看视频| 国产亚洲av高清不卡| 国产成人精品在线电影| 国产精品三级大全| 国产欧美日韩一区二区三区在线| 一级片免费观看大全| 亚洲,一卡二卡三卡| 精品视频人人做人人爽| 少妇人妻久久综合中文| 天堂俺去俺来也www色官网| 中文字幕人妻丝袜制服| 黄色怎么调成土黄色| 巨乳人妻的诱惑在线观看| 一区在线观看完整版| 欧美国产精品va在线观看不卡| 人人妻人人爽人人添夜夜欢视频| 国产男女内射视频| 精品人妻一区二区三区麻豆| 操出白浆在线播放| 亚洲成人免费电影在线观看 | 秋霞在线观看毛片| 满18在线观看网站| 国产成人a∨麻豆精品| 精品视频人人做人人爽| 十分钟在线观看高清视频www| 久久久久久免费高清国产稀缺| 国产成人精品久久久久久| 日日摸夜夜添夜夜爱| 赤兔流量卡办理| 国产97色在线日韩免费| 丝瓜视频免费看黄片| 久久精品成人免费网站| 女性被躁到高潮视频| 国产三级黄色录像| 校园人妻丝袜中文字幕| 国产一卡二卡三卡精品| 一级片'在线观看视频| 三上悠亚av全集在线观看| 久久免费观看电影| www日本在线高清视频| 九草在线视频观看| 9色porny在线观看| 国产午夜精品一二区理论片| 中文精品一卡2卡3卡4更新| 天天躁夜夜躁狠狠躁躁| 免费观看a级毛片全部| 亚洲国产精品成人久久小说| 黄频高清免费视频| 久久 成人 亚洲| 热99久久久久精品小说推荐| 国产亚洲欧美精品永久| 国产在线一区二区三区精| 成年人黄色毛片网站| 亚洲一区中文字幕在线| 免费观看人在逋| 看免费av毛片| 久久av网站| 国产激情久久老熟女| 男女午夜视频在线观看| 黄色片一级片一级黄色片| 午夜两性在线视频| 这个男人来自地球电影免费观看| 欧美精品亚洲一区二区| 满18在线观看网站| 国产在视频线精品| 啦啦啦 在线观看视频| 久久午夜综合久久蜜桃| 免费看不卡的av| 成年女人毛片免费观看观看9 | 久久性视频一级片| 亚洲三区欧美一区| a级毛片在线看网站| 亚洲欧美清纯卡通| 香蕉丝袜av| 国产av国产精品国产| 欧美老熟妇乱子伦牲交| 亚洲精品国产色婷婷电影| 国产精品麻豆人妻色哟哟久久| 免费在线观看黄色视频的| 秋霞在线观看毛片| 18在线观看网站| 9热在线视频观看99| 99精国产麻豆久久婷婷| 亚洲成人手机| 久久久国产精品麻豆| 黑人欧美特级aaaaaa片| 亚洲av美国av| 丰满迷人的少妇在线观看| 亚洲国产看品久久| 国产精品人妻久久久影院| 777久久人妻少妇嫩草av网站| 国产欧美日韩精品亚洲av| 国产精品 欧美亚洲| 精品国产超薄肉色丝袜足j| 国产男女超爽视频在线观看| 日日摸夜夜添夜夜爱| 亚洲成av片中文字幕在线观看| 国产精品三级大全| 久久久精品免费免费高清| 成年女人毛片免费观看观看9 | 欧美日韩福利视频一区二区| 美女视频免费永久观看网站| 亚洲熟女精品中文字幕| 欧美 日韩 精品 国产| 一区二区三区精品91| 国产精品99久久99久久久不卡| 国产有黄有色有爽视频| 女人爽到高潮嗷嗷叫在线视频| 九草在线视频观看| 99国产精品99久久久久| 国产精品一区二区在线观看99| 日日夜夜操网爽| 91精品三级在线观看| 在线观看国产h片| 欧美变态另类bdsm刘玥| 午夜久久久在线观看| 亚洲国产毛片av蜜桃av| av又黄又爽大尺度在线免费看| 精品人妻一区二区三区麻豆| 日韩大片免费观看网站| 无遮挡黄片免费观看| 操美女的视频在线观看| 波多野结衣av一区二区av| 亚洲国产av影院在线观看| 成年av动漫网址| 免费观看a级毛片全部| 久久精品久久久久久久性| 视频在线观看一区二区三区| 日韩一区二区三区影片| 欧美精品一区二区免费开放| 久久久精品区二区三区| 国产精品一国产av| 女人被躁到高潮嗷嗷叫费观| 国产国语露脸激情在线看| 午夜视频精品福利| 午夜老司机福利片| 人人妻,人人澡人人爽秒播 | 国产爽快片一区二区三区| 欧美人与性动交α欧美精品济南到| 建设人人有责人人尽责人人享有的| 亚洲激情五月婷婷啪啪| 久久人人爽人人片av| 日韩一本色道免费dvd| 国产精品99久久99久久久不卡| 男的添女的下面高潮视频| 在现免费观看毛片| 亚洲国产欧美一区二区综合| 亚洲欧美日韩另类电影网站| 青草久久国产| 黄色 视频免费看| 一区在线观看完整版| 亚洲国产欧美日韩在线播放| 亚洲天堂av无毛| 欧美日韩黄片免| 亚洲熟女精品中文字幕| 91字幕亚洲| 久久久精品国产亚洲av高清涩受| 精品人妻一区二区三区麻豆| 国产精品人妻久久久影院| 巨乳人妻的诱惑在线观看| 成人三级做爰电影| 久久国产精品男人的天堂亚洲| 国产精品人妻久久久影院| 巨乳人妻的诱惑在线观看| 国产亚洲欧美精品永久| 午夜免费男女啪啪视频观看| 成人18禁高潮啪啪吃奶动态图| 青春草视频在线免费观看| 又粗又硬又长又爽又黄的视频| 亚洲精品一区蜜桃| 桃花免费在线播放| 久久狼人影院| www.熟女人妻精品国产| 国产精品二区激情视频| av国产精品久久久久影院| 五月开心婷婷网| 精品少妇内射三级| 国产亚洲精品久久久久5区| 欧美性长视频在线观看| 侵犯人妻中文字幕一二三四区| 日韩一本色道免费dvd| 一级黄片播放器| 久久性视频一级片| 男人舔女人的私密视频| 精品一区二区三区四区五区乱码 | 亚洲av欧美aⅴ国产| 老司机在亚洲福利影院| 欧美在线一区亚洲| 日韩一区二区三区影片| 别揉我奶头~嗯~啊~动态视频 | 久久国产精品大桥未久av| 99九九在线精品视频| 丝袜美腿诱惑在线| 欧美日韩精品网址| 两性夫妻黄色片| 亚洲国产欧美日韩在线播放| 亚洲伊人久久精品综合| 日本午夜av视频| 中国国产av一级| 久久久久久久国产电影| 一级,二级,三级黄色视频| 欧美日韩视频精品一区| 欧美亚洲日本最大视频资源| 18在线观看网站| 97人妻天天添夜夜摸| 成人午夜精彩视频在线观看| 丰满迷人的少妇在线观看| 人妻 亚洲 视频| 欧美国产精品一级二级三级| 日韩精品免费视频一区二区三区| 亚洲欧洲日产国产| 国产伦理片在线播放av一区| 一区二区三区四区激情视频| 成人免费观看视频高清| 日韩免费高清中文字幕av| videosex国产| 一区二区三区乱码不卡18| 黄频高清免费视频| 免费日韩欧美在线观看| 美女大奶头黄色视频| 精品少妇一区二区三区视频日本电影| 欧美精品一区二区免费开放| 999久久久国产精品视频| 日本wwww免费看| 乱人伦中国视频| 首页视频小说图片口味搜索 | 日韩一区二区三区影片| 久久性视频一级片| 欧美日韩黄片免| 99热国产这里只有精品6| 亚洲少妇的诱惑av| 国产亚洲欧美精品永久| 欧美变态另类bdsm刘玥| 50天的宝宝边吃奶边哭怎么回事| 日韩熟女老妇一区二区性免费视频| 9191精品国产免费久久| 日本a在线网址| 97人妻天天添夜夜摸| 国产色视频综合| 蜜桃国产av成人99| 99久久精品国产亚洲精品| 在线观看人妻少妇| 多毛熟女@视频| 免费在线观看完整版高清| 国产99久久九九免费精品| 久久精品人人爽人人爽视色| 99国产精品一区二区三区| 欧美在线黄色| av片东京热男人的天堂| h视频一区二区三区| 中国国产av一级| 国产在线免费精品| 亚洲精品国产av成人精品| 欧美精品一区二区大全| 久久久精品区二区三区| 亚洲欧洲国产日韩| 国产在线观看jvid| 成人18禁高潮啪啪吃奶动态图| 老汉色av国产亚洲站长工具| 一本综合久久免费| 亚洲一区二区三区欧美精品| 国产精品久久久久久精品古装| 亚洲欧洲国产日韩| 国产淫语在线视频| 黄片播放在线免费| 亚洲av电影在线进入| 国产精品熟女久久久久浪| 亚洲激情五月婷婷啪啪| 欧美 日韩 精品 国产| 99热网站在线观看| 男女无遮挡免费网站观看| 91老司机精品| 90打野战视频偷拍视频| 无限看片的www在线观看| 日韩人妻精品一区2区三区| 国产精品一二三区在线看| 国产欧美日韩综合在线一区二区| tube8黄色片| 男女国产视频网站| 一本一本久久a久久精品综合妖精| 啦啦啦视频在线资源免费观看| 国产有黄有色有爽视频| 久久人人爽av亚洲精品天堂| 精品亚洲成国产av| 大话2 男鬼变身卡| 波多野结衣一区麻豆| av电影中文网址| 亚洲av美国av| 少妇的丰满在线观看| 国产爽快片一区二区三区| 1024视频免费在线观看| 久久精品久久精品一区二区三区| 亚洲成人免费av在线播放| 麻豆av在线久日| 亚洲国产精品成人久久小说| 国产黄频视频在线观看| 亚洲专区中文字幕在线| 日韩av免费高清视频| 人人妻人人爽人人添夜夜欢视频| 久久精品国产综合久久久| 国产精品久久久久久人妻精品电影 | 日韩视频在线欧美| 欧美精品av麻豆av| 亚洲av日韩精品久久久久久密 | 看免费成人av毛片| 不卡av一区二区三区| 国产精品国产三级专区第一集| 好男人电影高清在线观看| 大码成人一级视频| 熟女av电影| 热re99久久国产66热| cao死你这个sao货| 国产高清不卡午夜福利| 欧美日韩精品网址| 国产成人av教育| 精品视频人人做人人爽| 丁香六月欧美| 中文字幕另类日韩欧美亚洲嫩草| 黄色视频在线播放观看不卡| 亚洲欧洲国产日韩| 另类亚洲欧美激情| 丝袜美足系列| 亚洲精品久久久久久婷婷小说| 午夜免费成人在线视频| 欧美少妇被猛烈插入视频| 日韩熟女老妇一区二区性免费视频| 女人爽到高潮嗷嗷叫在线视频| 一级,二级,三级黄色视频| 狠狠精品人妻久久久久久综合| 亚洲综合色网址| 91精品三级在线观看| 亚洲欧美清纯卡通| 久久精品亚洲av国产电影网| 国产亚洲av片在线观看秒播厂| 成人影院久久| 精品国产乱码久久久久久男人| 午夜两性在线视频| 国产主播在线观看一区二区 | 美女脱内裤让男人舔精品视频| 老司机影院成人| 后天国语完整版免费观看| 欧美精品高潮呻吟av久久| 黄色怎么调成土黄色| 成年美女黄网站色视频大全免费| 另类精品久久| av又黄又爽大尺度在线免费看| 男女床上黄色一级片免费看| 免费在线观看视频国产中文字幕亚洲 | 在线观看免费日韩欧美大片| 国产亚洲精品第一综合不卡| 黄色a级毛片大全视频| 国产亚洲精品久久久久5区| 999精品在线视频| 亚洲国产欧美在线一区| 国产在线视频一区二区| 国产老妇伦熟女老妇高清| 99香蕉大伊视频| 国产精品香港三级国产av潘金莲 | 亚洲av日韩在线播放| 老汉色∧v一级毛片| 亚洲综合色网址| 久久精品人人爽人人爽视色| 国产精品国产三级国产专区5o| 18禁裸乳无遮挡动漫免费视频| 国产欧美日韩精品亚洲av| 国产老妇伦熟女老妇高清| 国产成人av教育| 亚洲人成网站在线观看播放| 国产在视频线精品| 后天国语完整版免费观看| 在线观看免费视频网站a站| 亚洲中文日韩欧美视频| 18在线观看网站| 国产精品久久久av美女十八| 免费在线观看日本一区| 亚洲黑人精品在线| 欧美+亚洲+日韩+国产| 激情五月婷婷亚洲| 中文字幕人妻丝袜制服| 国产成人欧美| 日本猛色少妇xxxxx猛交久久| 久久精品人人爽人人爽视色| 韩国高清视频一区二区三区| 国产精品人妻久久久影院| 看免费成人av毛片| 一级黄片播放器| 在线av久久热| 日本黄色日本黄色录像| 黄色视频在线播放观看不卡| e午夜精品久久久久久久| 人人妻人人爽人人添夜夜欢视频| 国产精品久久久久久精品电影小说| 亚洲欧美成人综合另类久久久| 久久天堂一区二区三区四区| 亚洲精品国产av成人精品| 99热国产这里只有精品6| 宅男免费午夜| 免费日韩欧美在线观看| 亚洲国产av新网站|