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

    基于磷蝦群算法的汽輪機組最優(yōu)初壓研究

    2017-09-03 09:26:39牛培峰劉阿玲馬云鵬李國強
    動力工程學報 2017年8期
    關鍵詞:初壓熱耗率磷蝦

    牛培峰, 陳 科, 劉阿玲, 馬云鵬, 趙 振, 李國強

    (1. 燕山大學 電氣工程學院, 河北秦皇島 066004; 2. 山東省萊西市院上鎮(zhèn)中心中學, 山東青島 266600)

    基于磷蝦群算法的汽輪機組最優(yōu)初壓研究

    牛培峰1, 陳 科1, 劉阿玲2, 馬云鵬1, 趙 振1, 李國強1

    (1. 燕山大學 電氣工程學院, 河北秦皇島 066004; 2. 山東省萊西市院上鎮(zhèn)中心中學, 山東青島 266600)

    為找到汽輪機變工況運行時的最優(yōu)初壓,利用改進的磷蝦群算法(A-KH)和快速學習網(wǎng)(FLN)建立熱耗率預測模型,然后利用A-KH算法的全局搜索能力,在可行的壓力區(qū)間內(nèi)對所建模型熱耗率最低時對應的主蒸汽壓力進行尋優(yōu),并將優(yōu)化后的最優(yōu)初壓曲線與廠家設計壓力曲線進行對比.結果表明:優(yōu)化后的最優(yōu)初壓曲線能有效降低汽輪機組的熱耗率,對汽輪機的安全經(jīng)濟運行更具有指導意義.

    汽輪機; 熱耗率; 最優(yōu)初壓; 磷蝦群算法; 快速學習網(wǎng)

    我國用電結構的改變導致電網(wǎng)的峰谷差增大,也就意味著大型機組不得不進行調(diào)峰運行,汽輪機組經(jīng)常在變負荷狀態(tài)下運行,為了使機組能夠在變負荷狀態(tài)下保持較高的熱效率,常采用復合滑壓運行方式,即定-滑-定運行方式.因此,確定定-滑-定運行分界點和尋找最優(yōu)初壓對提高熱電廠的熱經(jīng)濟性具有重要作用.

    目前,熱電廠汽輪機的運行初壓一般都是按照廠家提供的滑壓運行曲線來設定的,并沒有考慮汽輪機組運行的實際情況,無法正確描述機組在變負荷狀態(tài)下的運行狀態(tài).因此,要確保汽輪機組在變負荷運行時依然保持較好的運行情況,就必須對汽輪機的運行初壓進行優(yōu)化.在一定程度上,對汽輪機的運行初壓進行優(yōu)化可以歸結為求解熱耗率最小時對應的主蒸汽壓力,即最優(yōu)初壓[1],其確定方法有2種:理論尋優(yōu)法和實驗尋優(yōu)法[2-3].理論尋優(yōu)法利用汽輪機各級的幾何參數(shù)進行熱力學計算,從而得出最優(yōu)初壓,但這種方法計算量很大且所需數(shù)據(jù)較難采集;實驗尋優(yōu)法通過選定幾組不同的主蒸汽壓力進行實驗,并計算出對應主蒸汽壓力的熱耗率,對比實驗結果找出最優(yōu)初壓,但這種方法選取的初壓點有限,可能導致所尋壓力并非是最優(yōu)初壓.

    磷蝦群算法(KH)[4]是一種新的智能優(yōu)化算法,由于其操作簡單、搜索多樣性強和調(diào)整的參數(shù)較少而應用廣泛.然而在對復雜問題進行優(yōu)化時,KH算法容易陷入局部最優(yōu).為解決上述問題,筆者提出了一種改進的磷蝦群算法(A-KH),實驗表明,A-KH算法具有較高的收斂精度,且跳出局部最優(yōu)的能力也得到加強.

    以600 MW超臨界汽輪機組為研究對象,筆者采用A-KH算法和快速學習網(wǎng)(FLN)[5]建立了汽輪機組熱耗率預測模型,該模型具有較高的預測精度和較強的泛化能力,能夠很好地對汽輪機組的熱耗率進行預測.然后,利用A-KH算法去搜索相應負荷下熱耗率最低時的主蒸汽壓力,最終給出了汽輪機最優(yōu)初壓曲線.

    1 磷蝦群算法和快速學習網(wǎng)

    1.1 磷蝦群算法

    KH算法的靈感來自于磷蝦覓取食物和相互之間的交流[6].在KH算法中,每只磷蝦位置代表一個可行解.在海洋生活中,磷蝦個體的位置會受到以下3個因素影響[7].

    (1)種群位置遷移引起的個體游動.

    Ni,new=Nmaxαi+ωnNi,old

    (1)

    式中:αi為移動方向,αi=αi,local+αi,target,αi,local為局部鄰近個體的感應方向向量,αi,target為種群最優(yōu)個體的方向向量;Nmax為最大感應速度,取0.01 m/s;ωn為慣性權值,其值范圍為[0,1];Ni,old為上次個體游動的位置變化;Ni,new為下一次個體游動的位置變化.

    (2) 覓食行為.

    Fi=Vfβi+ωfFi,old

    (2)

    式中:βi為磷蝦個體覓食的方向向量,βi=βi,food+βi,best,βi,food為食物源方向向量,βi,best為當前最優(yōu)磷蝦位置的方向向量;Vf為覓食的速度,取0.02 m/s;ωf為慣性權值,其值范圍為[0,1];Fi,old為第i個磷蝦個體上次覓食運動的位置變化;Fi為第i個磷蝦個體覓食運動的位置變化.

    (3) 磷蝦個體的隨機擴散.

    Di=Dmaxδ

    (3)

    式中:Dmax為磷蝦個體最大的擴散速度;δ為一個隨機的擴散方向向量,其值范圍為[-1,1];Di為第i個磷蝦個體隨機擴散的位置變化.

    針對不同的磷蝦個體,越靠近食物源位置,擴散現(xiàn)象越不明顯.上述3個因素共同決定了磷蝦個體的游動方向,即朝著適應度值最小的方向改變.其中種群位置遷移引起的個體游動和覓食行為都包含一個全局搜索策略和一個局部搜索策略,2種策略并行進行,使得該算法成為一種穩(wěn)定有效的優(yōu)化算法[8].

    綜上所述,從t到t+Δt時間內(nèi)的位置矢量為:

    (4)

    式中:ΔT為速度矢量的步長調(diào)節(jié)因子.

    ΔT的定義如下:

    (5)

    式中:Uj和Lj分別為第j個變量的上、下邊界;N為變量總數(shù);Ct的取值范圍為[0,2].

    1.2 快速學習網(wǎng)

    FLN是一種新型雙并聯(lián)前饋神經(jīng)網(wǎng)絡,其最大的特點是在單隱藏層神經(jīng)網(wǎng)絡的基礎上增加了輸入層與輸出層的直接聯(lián)系[9].因此,F(xiàn)LN可以看成是一種隱藏層到輸出層的非線性與輸入層到輸出層的線性組合模型,其應用前景廣泛.

    對于任意N個隨機樣本{(xi,yi),i=1,2,…,N},其中,xi=[xi1,xi2,…,xin]T∈Rn表示第i個樣本的n維特征向量,yi=[yi1,yi2,…,yil]∈Rl表示第i個樣本的l維輸出向量.令隱藏層神經(jīng)元數(shù)為m個,隱藏層激活函數(shù)為g(x),則FLN的輸出如下:

    (6)

    式中:j=1,2,…,N;Woi為輸入層與輸出層之間的連接權值;Wk,in為輸入層到第k個隱藏層神經(jīng)元之間的連接權值;Wk,oh為第k個隱藏層神經(jīng)元到輸出層之間的連接權值;bk為第k個隱藏層神經(jīng)元的閾值.

    式(6)用矩陣形式可表達為:

    (7)

    (8)

    (9)

    式中:W為輸出權值矩陣;G為隱藏層輸出矩陣;Y為期望輸出;l為輸出層節(jié)點個數(shù);Woi和Woh分別為輸入層到輸出層的權值矩陣和隱藏層到輸出層的權值矩陣.

    式(7)用最小二乘范數(shù)解求法可得

    (10)

    (11)

    FLN算法步驟如下:(1) 隨機生成輸入權值矩陣Win和隱藏層閾值矩陣b;(2) 通過式(9)計算隱藏層輸出矩陣G;(3) 通過式(10)計算輸出權值矩陣W;(4) 通過式(11)將輸出權值矩陣W分割為Woi和Woh.

    2 改進的磷蝦群算法

    與其他優(yōu)化算法相比,KH算法的優(yōu)勢在于其搜索多樣性強、操作簡單和調(diào)節(jié)的參數(shù)少[10].但KH算法中的搜索完全依賴隨機性,導致該算法收斂精度低,且不易跳出局部最優(yōu)[11].針對以上KH算法的缺陷,筆者提出了A-KH算法.改進后的算法具有較強的全局搜索能力,同時收斂速度和收斂精度也有較大幅度的提高,主要有以下4個改進點.

    (1) 采用反向?qū)W習算法[12]進行種群位置初始化,可以有效提高初始化種群位置的質(zhì)量.

    在A-KH算法中,基于反向?qū)W習算法的種群位置初始化過程如下(其中NP為可行解的維數(shù),NK為種群數(shù)量,G為當前迭代次數(shù)):

    ① 初始化種群位置P(G=0)={xij},i=1,2,…,NP,j=1,2,…,NK.

    (12)

    式中:xi,min和xi,max分別為第i維元素的最小值和最大值.

    ③ 從組合種群位置{P(G=0)∪P′(G=0)}中選擇NK個適應度值較小的位置作為初始種群位置.

    (2) 從式(5)可以看出,Ct在KH算法中扮演著重要的角色,為了更好地平衡算法前期的探索能力和后期的開采能力,引入線性遞減的Ct:

    (13)

    式中:Ct∈[0,2];a=0.45,b=0.1;Gmax為最大迭代次數(shù).

    (3) 從式(1)和式(2)可以看出,慣性權值(wn,wf)決定了KH算法的收斂速度和位置更新方式.在原始的KH算法中,(wn,wf)隨著迭代次數(shù)的增加呈線性遞減,這有可能會影響該算法的收斂精度和收斂速度.為了解決該問題,引入了sine混沌圖譜[13]來更新(wn,wf)的值.

    (14)

    式中:?=0.8;xk∈(0,1);c∈(0,4];θ=0.1;k=1,2,…,Gmax.

    Xi(t+Δt)=Xi(t)wij+

    (15)

    wij=φ=1/[1+exp(-fj/h)G]

    (16)

    (17)

    式中:h為在第一次迭代中種群位置最優(yōu)的適應度值;fj為第j個磷蝦的適應度值.

    A-KH算法的流程如下:(1) 初始化參數(shù)設定,即種群數(shù)量NK、最大迭代次數(shù)Gmax、最大感應速度Nmax、覓食速度Vf和磷蝦個體的最大擴散速度Dmax等;(2) 在可行域空間內(nèi)利用式(12)反向?qū)W習算法初始化磷蝦種群的位置;(3) 計算每一個磷蝦相應的適應度值,并選擇NK個適應度值較小的位置作為初始種群位置;(4) 采用式(14)sine混沌圖譜更新慣性權值(wn,wf);(5) 計算由磷蝦種群位置遷移引起的個體游動、覓食行為和隨機擴散對磷蝦個體位置的變化量;(6) 采用式(13)計算Ct;(7) 利用式(15)產(chǎn)生下一代種群的位置;(8) 計算新產(chǎn)生位置矢量的適應度值,重復步驟(5)~步驟(8),直到達到最大迭代次數(shù),算法結束;(9) 輸出最優(yōu)個體,即找到的最優(yōu)解.

    為了更好地體現(xiàn)A-KH算法搜索的高效性,采用6個典型的數(shù)值測試函數(shù)對KH算法和A-KH算法進行仿真實驗,6個測試函數(shù)見表1,其中前3個為單峰測試函數(shù),后3個為多峰測試函數(shù).A-KH算法與KH算法的參數(shù)設置相同:種群數(shù)量為40,最大迭代次數(shù)為500,可行解的維數(shù)設置為30,運行次數(shù)為20,最大感應速度Nmax=0.01 m/s,覓食速度Vf=0.02 m/s,磷蝦個體最大擴散速度Dmax=0.005 m/s.對20次尋優(yōu)結果的平均值和均方差進行記錄,結果見表2.

    從表2可以看出,在設定的參數(shù)條件下, A-KH算法具有更優(yōu)的尋優(yōu)效果和更高的收斂精度.在對f1、f2、f4、f5和f6函數(shù)進行尋優(yōu)時,A-KH算法能夠找到理論最優(yōu)值或非常接近,可以近似看成理論最優(yōu)值.對于函數(shù)f3,A-KH算法搜索的精度雖然不高,但也比KH算法提高了1個數(shù)量級.與KH算法相比,A-KH算法不僅搜索的精度更高,而且跳出局部最優(yōu)的能力也更強,可以應用于模型參數(shù)的尋優(yōu)和實際工程問題的優(yōu)化.

    表1 測試函數(shù)

    表2 2種算法測試效果的對比

    3 汽輪機組熱耗率建模

    汽輪機組的熱耗率是指汽輪發(fā)電機組每發(fā)1 kW·h電量所消耗的熱量,是反映機組能量轉(zhuǎn)換過程中一項重要的經(jīng)濟指標[14].選擇熱耗率預測模型參數(shù)時,通過比較輸入?yún)?shù)和輸出參數(shù)對熱力性能實驗結果的影響,并結合相關資料[8,15],選擇發(fā)電負荷、主蒸汽壓力、主蒸汽溫度、再熱器出口蒸汽壓力、再熱器出口蒸汽溫度、再熱器入口蒸汽壓力、再熱器入口蒸汽溫度、再熱減溫水流量、過熱減溫水流量、汽輪機背壓、循環(huán)水進水口溫度和給水流量共12個參數(shù)作為熱耗率預測模型的輸入,熱耗率作為輸出.其中,主蒸汽壓力對熱耗率值有直接的影響,熱耗率的計算詳見文獻[16].

    3.1 模型建立及參數(shù)優(yōu)化

    以某熱電廠600 MW超臨界汽輪機組為研究對象,從集散控制系統(tǒng)(DCS)中采集10 d機組正常運行數(shù)據(jù),其中春、夏、秋、冬四季各采集的天數(shù)為2、2、3、3,每隔2 h采集一次,每天12組,每一組均包含上述12個自變量參數(shù),總共采集了120組多工況運行數(shù)據(jù),發(fā)電負荷的范圍在298.455~563.555 MW內(nèi),其他11個參數(shù)的變化區(qū)間均在額定范圍內(nèi).其中,前8 d共96組運行數(shù)據(jù)作為訓練樣本,后2 d共24組運行數(shù)據(jù)作為測試樣本,用于檢測所建模型的泛化能力和預測精度.

    采用FLN模型建立熱耗率預測模型,由于FLN的輸入權值和隱藏層閾值是隨機產(chǎn)生的,很難確保所建的熱耗率預測模型具有較強的泛化能力和較高的預測精度.因此,采用A-KH算法和FLN模型進行綜合建模,A-KH算法用來對FLN模型的輸入權值和隱藏層閾值進行尋優(yōu),目標函數(shù)以適應度值最小為原則,以最大迭代次數(shù)為循環(huán)終止條件,當尋優(yōu)結束后,將最優(yōu)尋優(yōu)參數(shù)作為FLN模型的輸入權值和隱藏層閾值,即完成A-KH-FLN熱耗率預測模型(以下簡稱A-KH-FLN模型)的建立.A-KH算法的參數(shù)設置與前面相同.FLN模型的參數(shù)設置如下:3層網(wǎng)絡結構為12-20-1,隱藏層激活函數(shù)為sigmoid,輸入權值和隱藏層閾值的尋優(yōu)范圍均為[-1,1],通過多次實驗,確定循環(huán)迭代次數(shù)為200.具體流程圖見圖1.

    在模型結構參數(shù)優(yōu)化中,目標函數(shù)為:

    i=1,2,…,N

    (18)

    圖1 采用A-KH算法優(yōu)化FLN模型結構參數(shù)的流程圖

    3.2 模型性能分析

    圖2給出了A-KH-FLN模型對訓練樣本的熱耗率預測情況.由圖2可知,A-KH-FLN模型的預測值與實際值的跟蹤情況良好,表明A-KH-FLN模型對訓練樣本具有很高的擬合精度和較強的跟蹤能力.

    圖2 訓練樣本熱耗率預測值與實際值的對比

    Fig.2 Comparison of training data between predicted results and actual measurements

    圖3給出了A-KH-FLN模型對測試樣本的熱耗率預測情況.由圖3可知,A-KH-FLN模型具有較高的預測精度和較強的泛化能力.

    圖3 測試樣本熱耗率預測值與實際值的對比

    Fig.3 Comparison of testing data between predicted results and actual measurements

    為了更好地體現(xiàn)A-KH-FLN模型的預測效果,表3給出了訓練樣本和測試樣本的均方根誤差(RMSE)、平均絕對百分比誤差(MAPE)和平均絕對誤差(MAE)3個性能評價指標,其中KH-FLN模型是由KH算法和FLN模型建立的熱耗率預測模型.由表3可知,A-KH-FLN模型訓練樣本的RMSE為6.106 1,MAPE為6.042 3×10-9,MAE為4.615 8;測試樣本的RMSE為11.413 3,MAPE為2.091 2×10-5,MAE為9.322 3,各項指標均優(yōu)于其他2種模型,說明A-KH-FLN模型具有良好的預測性能和泛化能力,能夠用于工程應用.

    表3 各模型對熱耗率的預測精度

    圖4為3種模型對測試樣本的預測誤差曲線.

    圖4 各模型對測試樣本的預測誤差曲線

    由圖4可知,相比于其他2種模型,A-KH-FLN模型的誤差曲線波動平穩(wěn),誤差較小,說明A-KH-FLN模型的泛化能力較穩(wěn)定.

    4 汽輪機組最優(yōu)初壓研究

    4.1 優(yōu)化問題描述

    汽輪機在變工況條件下運行時,有且只有主蒸汽壓力和主蒸汽流量是可以調(diào)節(jié)的.根據(jù)線性規(guī)劃條件式可知,給定主蒸汽壓力,主蒸汽流量也就可以確定了.因此,在某一工況下,主蒸汽壓力值與汽輪機組熱耗率值是一一對應的.當熱耗率最低時對應的主蒸汽壓力為最優(yōu)初壓.相關的數(shù)學模型描述如下:

    (19)

    式中:HR為汽輪機組的熱耗率,kJ/(kW·h);NgE為汽輪機組運行負荷,MW;po為當前主蒸汽壓力,MPa;X為影響熱耗率的其他因素;Ngmin和Ngmax分別為汽輪機組允許的最小和最大負荷,MW;pod為機組的額定主蒸汽壓力,MPa;Ngd為汽輪機組額定負荷,MW.

    4.2 采用A-KH算法優(yōu)化汽輪機初壓

    A-KH算法參數(shù)設置如下:種群數(shù)量為10,最大迭代次數(shù)為50,其余參數(shù)設置與第2節(jié)相同.在所建熱耗率預測模型的基礎上,以最大迭代次數(shù)為尋優(yōu)終止條件,采用A-KH算法對主蒸汽壓力進行尋優(yōu),并找出該負荷下對應的最優(yōu)初壓.

    在全部工況中,挑選5種典型工況:300 MW、360 MW、420 MW、480 MW和540 MW.5種工況優(yōu)化前后的實驗結果見表 4.由表4可知,420 MW下優(yōu)化后熱耗率值下降最多,為56.19 kJ/(kW·h),各負荷下對應的熱耗率降低值平均為28.35 kJ/(kW·h).在電廠中,通過A-KH算法優(yōu)化得到最優(yōu)初壓,將其作為該負荷下的主蒸汽壓力設定值,可以有效降低汽輪機組的熱耗率,長期如此運行能夠提高電廠的熱效率,同時可以帶來可觀的經(jīng)濟效益.

    圖5給出了汽輪機組最優(yōu)初壓的變化曲線.分析圖5可知,汽輪機變工況運行時,最優(yōu)初壓是一條曲率隨負荷增加逐漸降低的平滑曲線,而不是廠家提供的一條斜直線.結果表明,經(jīng)過A-KH算法優(yōu)化后的汽輪機最優(yōu)初壓曲線能夠更好地指導汽輪機組的安全運行.

    表4 各負荷下優(yōu)化前后的結果對比

    圖5 最優(yōu)初壓曲線與設計壓力曲線的對比

    5 結 論

    (1) 針對KH算法易陷入局部最優(yōu)且收斂精度低等缺點,提出了A-KH算法,經(jīng)測試函數(shù)驗證,A-KH算法具有更高的搜索精度和更強的跳出局部最優(yōu)的能力.

    (2) A-KH-FLN模型能夠準確反映輸入與輸出之間的非線性映射關系,相比于其他模型,該模型具有更高的預測精度和更強的泛化能力.

    (3) 在熱耗率預測模型基礎上,利用A-KH算法的全局搜索能力對汽輪機組熱耗率最低時的主蒸汽壓力進行尋優(yōu),相比于廠家設計壓力曲線,得到的最優(yōu)初壓曲線能有效降低汽輪機組的熱耗率,對汽輪機的安全經(jīng)濟運行更具有指導意義.

    [1] 劉偉, 葉亞蘭, 司風琪, 等. 基于BP神經(jīng)網(wǎng)絡和SA-BBO算法的汽輪機組最優(yōu)運行初壓的確定[J]. 熱能動力工程, 2013, 28(1): 18-22.

    LIU Wei, YE Yalan, SI Fengqi, et al. Determination of the optimum initial operation pressure of a steam turbine unit based on a BP (back propagation) neural network and SA-BBO (simulated annealing biogeography-based optimization) algorithm[J]. Journal of Engineering for Thermal Energy and Power, 2013, 28(1): 18-22.

    [2] 郭江龍, 常澍平, 姚力強, 等. 大型汽輪機復合滑壓運行參數(shù)尋優(yōu)方法研究[J]. 汽輪機技術, 2010, 52(6): 467-469.

    GUO Jianglong, CHANG Shuping, YAO Liqiang, et al. Method research on parameter optimization in composite sliding pressure operation of large scale steam turbine[J]. Turbine Technology, 2010, 52(6): 467-469.

    [3] 周志平, 范鑫, 李明, 等. 超臨界600 MW機組滑壓運行參數(shù)優(yōu)化分析與試驗[J]. 熱力發(fā)電, 2011, 40(10): 50-54.

    ZHOU Zhiping, FAN Xin, LI Ming, et al. Analysis and test of parameters' optimization for sliding pressure operation of supercritical 600 MW unit[J]. Thermal Power Generation, 2011, 40(10): 50-54.

    [4] GANDOMI A H, ALAVI A H. Krill herd: a new bio-inspired optimization algorithm[J]. Communications in Nonlinear Science and Numerical Simulation, 2012, 17(12): 4831-4845.

    [5] LI Guoqiang, NIU Peifeng, DUAN Xiaolong, et al. Fast learning network: a novel artificial neural network with a fast learning speed[J]. Neural Computing and Applications, 2014, 24(7/8): 1683-1695.

    [6] WANG Gaige, GANDOMI A H, ALAVI A H. Stud krill herd algorithm[J]. Neurocomputing, 2014, 128: 363-370.

    [7] WANG Gaige, GUO Lihong, WANG Heqi, et al. Incorporating mutation scheme into krill herd algorithm for global numerical optimization[J]. Neural Computing and Applications, 2014, 24(3/4): 853-871.

    [8] 牛培峰, 楊瀟, 馬云鵬, 等. 基于改進的磷蝦群優(yōu)化算法的汽輪機初壓優(yōu)化研究[J]. 動力工程學報, 2015, 35(9): 709-714.

    NIU Peifeng, YANG Xiao, MA Yunpeng, et al. Optimization on initial pressure of a steam turbine based on improved krill herd algorithm[J]. Journal of Chinese Society of Power Engineering, 2015, 35(9): 709-714.

    [9] 李國強. 新型人工智能技術研究及其在鍋爐燃燒優(yōu)化中的應用[D]. 秦皇島: 燕山大學, 2013.

    [10] LI Junpeng, TANG Yinggan, HUA Changchun, et al. An improved krill herd algorithm: krill herd with linear decreasing step[J]. Applied Mathematics and Computation, 2014, 234: 356-367.

    [11] WANG Gaige, GANDOMI A H, ALAVI A H. An effective krill herd algorithm with migration operator in biogeography-based optimization[J]. Applied Mathematical Modelling, 2014, 38(9/10): 2454-2462.

    [12] RAHNAMAYAN S, TIZHOOSH H R, SALAMA M M A. Opposition-based differential evolution (ODE) with variable jumping rate[C]//Proceedings of IEEE Symposium on Foundations of Computational Intelligence. Honolulu, USA: IEEE, 2007: 81-88.

    [13] WANG Gaige, GUO Lihong, GANDOMI A H, et al. Chaotic krill herd algorithm[J]. Information Sciences, 2014, 274: 17-34.

    [14] 劉超, 牛培峰, 游霞. 反向建模方法在汽輪機熱耗率建模中的應用[J]. 動力工程學報, 2014, 34(11): 867-872, 902.

    LIU Chao, NIU Peifeng, YOU Xia. Application of reversed modeling method in prediction of steam turbine heat rate[J]. Journal of Chinese Society of Power Engineering, 2014, 34(11): 867-872, 902.

    [15] 云曦, 閻維平. 火電廠汽輪機組影響熱耗率計算的因素[J]. 東北電力技術, 2007, 28(3): 15-18.

    YUN Xi, YAN Weiping. Factors effecting heat consumption calculation for steamed turbine of fossil-fired power plant[J]. Northeast Electric Power Technology, 2007, 28(3): 15-18.

    [16] 盛德仁, 李蔚, 陳堅紅, 等. 汽輪機組熱耗率分析及實時計算[J]. 熱力發(fā)電, 2003, 32(5): 16-18.

    SHENG Deren, LI Wei, CHEN Jianhong, et al. Analysis and real-time calculation of specific heat consumption for steam turbine unit[J]. Thermal Power Generation, 2003, 32(5): 16-18.

    Study on the Optimal Initial Pressure of a Steam Turbine Unit Based on Krill Herd Algorithm

    NIUPeifeng1,CHENKe1,LIUAling2,MAYunpeng1,ZHAOZhen1,LIGuoqiang1

    (1. School of Electrical Engineering, Yanshan University, Qinhuangdao 066004, Hebei Province, China; 2. Yuanshang Town Central Middle School, Qingdao 266600, Shandong Province, China)

    To obtain the optimal initial pressure of a steam turbine under variable working conditions, a heat rate prediction model was established using operating data of the turbine by improved krill herd algorithm (A-KH) and fast learning network (FLN), with which the optimal initial steam pressure corresponding to the lowest heat rate was searched in the range of premitted pressure by taking use of the global searching ability of A-KH algorithm, and the optimal initial pressure curve was subsequently compared with the design curve. Results show that via the optimal initial pressure curve, the heat rate could be lowered effectively, which may serve as a reference for safety and economic operation of steam turbines.

    steam turbine; heat rate; optimal initial pressure; krill herd algorithm; fast learning network

    1674-7607(2017)08-0615-07

    TK267

    A

    470.30

    2016-08-30

    2016-10-08

    國家自然科學基金資助項目(61573306,61403331);河北省自然科學基金資助項目(F2016203427);中國博士后科學基金資助項目(2015M571280)

    牛培峰(1958-),男,吉林舒蘭人,教授,博士生導師,研究方向為復雜工業(yè)系統(tǒng)的智能建模與智能控制和流程工業(yè)綜合自動化. 陳 科(通信作者),男,碩士研究生,電話(Tel.):15603372641;E-mail:m15603372641@163.com.

    猜你喜歡
    初壓熱耗率磷蝦
    磷蝦真是“蝦無敵”
    基于低碳理念的重載瀝青路面輪胎初壓工藝
    價值工程(2022年6期)2022-02-21 03:47:54
    瀝青路面碾壓施工膠輪初壓工藝
    價值工程(2022年5期)2022-02-21 02:25:30
    南極磷蝦粉在水產(chǎn)飼料中的應用
    湖南飼料(2021年4期)2021-10-13 07:32:46
    320MW抽凝供熱機組蝶閥開度調(diào)整試驗與仿真研究
    機電信息(2021年5期)2021-02-22 06:47:26
    “美味”的磷蝦
    降低燃氣發(fā)電機組煤耗方法探討
    魅力中國(2016年52期)2017-09-01 12:02:09
    基于相關向量機的汽輪機最優(yōu)運行初壓的確定
    化工學報(2016年9期)2016-10-14 00:03:26
    西門子二次再熱機組限壓/初壓切換過程典型故障分析
    廣西電力(2016年5期)2016-07-10 09:16:44
    “美味”的磷蝦
    最近2019中文字幕mv第一页| 一本色道久久久久久精品综合| 亚洲电影在线观看av| 大陆偷拍与自拍| 夜夜看夜夜爽夜夜摸| 天天躁夜夜躁狠狠久久av| 91午夜精品亚洲一区二区三区| 久久久午夜欧美精品| 99热国产这里只有精品6| 成人黄色视频免费在线看| 亚洲美女视频黄频| 国产黄a三级三级三级人| 亚洲最大成人av| 一二三四中文在线观看免费高清| 色婷婷久久久亚洲欧美| 亚洲精品第二区| 亚洲精品国产av蜜桃| 新久久久久国产一级毛片| 精品一区二区三卡| 内射极品少妇av片p| 亚洲av二区三区四区| 国产午夜福利久久久久久| 久久久亚洲精品成人影院| 日本午夜av视频| 国产精品一及| 成人亚洲欧美一区二区av| 午夜免费观看性视频| 久久精品人妻少妇| 国产成人精品一,二区| 国产淫片久久久久久久久| 国产成人免费观看mmmm| 1000部很黄的大片| 美女脱内裤让男人舔精品视频| 亚洲丝袜综合中文字幕| 亚洲欧美中文字幕日韩二区| 日日摸夜夜添夜夜爱| 亚洲内射少妇av| 我的女老师完整版在线观看| 嫩草影院新地址| 亚洲最大成人手机在线| 午夜精品国产一区二区电影 | 国产精品人妻久久久影院| 大片电影免费在线观看免费| 成人无遮挡网站| 少妇熟女欧美另类| 免费人成在线观看视频色| 免费观看av网站的网址| 免费播放大片免费观看视频在线观看| av在线播放精品| 国产亚洲av嫩草精品影院| 美女高潮的动态| 伊人久久国产一区二区| 最近中文字幕2019免费版| 日产精品乱码卡一卡2卡三| 少妇人妻久久综合中文| 精品久久久精品久久久| 蜜桃亚洲精品一区二区三区| 99久久九九国产精品国产免费| 欧美精品一区二区大全| 亚洲精品色激情综合| 免费播放大片免费观看视频在线观看| 91精品伊人久久大香线蕉| 国产美女午夜福利| tube8黄色片| 成人亚洲精品一区在线观看 | 久久精品国产亚洲av天美| 成人欧美大片| 美女脱内裤让男人舔精品视频| 国产精品无大码| 大码成人一级视频| 久久精品久久精品一区二区三区| 丝瓜视频免费看黄片| 另类亚洲欧美激情| 纵有疾风起免费观看全集完整版| 最近最新中文字幕免费大全7| 国产亚洲午夜精品一区二区久久 | 欧美激情国产日韩精品一区| 国产爱豆传媒在线观看| 大码成人一级视频| 日本一二三区视频观看| 直男gayav资源| 在线播放无遮挡| 国产成人精品久久久久久| 美女主播在线视频| 欧美性感艳星| 最近2019中文字幕mv第一页| 国产欧美日韩精品一区二区| 亚洲欧美一区二区三区黑人 | 国产毛片在线视频| 亚洲成人久久爱视频| 18禁裸乳无遮挡免费网站照片| 国产精品精品国产色婷婷| 亚洲综合色惰| 国产精品偷伦视频观看了| 亚洲激情五月婷婷啪啪| 日本黄大片高清| 国产成人精品福利久久| 国产精品一区二区三区四区免费观看| 国产有黄有色有爽视频| 视频区图区小说| 亚洲精品国产色婷婷电影| 丰满少妇做爰视频| 内地一区二区视频在线| 在线精品无人区一区二区三 | 色5月婷婷丁香| 寂寞人妻少妇视频99o| 又大又黄又爽视频免费| 久久久久久久午夜电影| 日本欧美国产在线视频| 中文字幕人妻熟人妻熟丝袜美| 国产成人一区二区在线| 日韩av不卡免费在线播放| 七月丁香在线播放| 中文精品一卡2卡3卡4更新| 99视频精品全部免费 在线| 日韩欧美精品免费久久| 国产一区二区亚洲精品在线观看| 人妻夜夜爽99麻豆av| 国产精品不卡视频一区二区| 国产色爽女视频免费观看| 狠狠精品人妻久久久久久综合| 三级国产精品片| 午夜福利视频1000在线观看| 日韩免费高清中文字幕av| 69av精品久久久久久| 国产成人免费无遮挡视频| 久久人人爽av亚洲精品天堂 | 国产午夜精品久久久久久一区二区三区| 日产精品乱码卡一卡2卡三| eeuss影院久久| 亚洲经典国产精华液单| 欧美日本视频| 黄色配什么色好看| 蜜臀久久99精品久久宅男| 2018国产大陆天天弄谢| 下体分泌物呈黄色| 老女人水多毛片| 午夜免费男女啪啪视频观看| av网站免费在线观看视频| 97超视频在线观看视频| 2021少妇久久久久久久久久久| 久久精品国产鲁丝片午夜精品| 国产成人freesex在线| 80岁老熟妇乱子伦牲交| 国产免费一级a男人的天堂| 国产精品一及| 久久国内精品自在自线图片| 国产日韩欧美在线精品| 国产精品一区二区性色av| 久久久久久久久大av| 国产精品久久久久久av不卡| 国产91av在线免费观看| 少妇高潮的动态图| 最新中文字幕久久久久| 欧美成人一区二区免费高清观看| 乱码一卡2卡4卡精品| 亚洲精品乱码久久久v下载方式| 中文资源天堂在线| 亚洲精品视频女| 亚洲精品乱久久久久久| 欧美性感艳星| 如何舔出高潮| 一级爰片在线观看| 亚洲经典国产精华液单| 国内揄拍国产精品人妻在线| 国产免费一区二区三区四区乱码| 国产成人a区在线观看| 一个人看的www免费观看视频| xxx大片免费视频| 99久久精品国产国产毛片| 青春草国产在线视频| 伊人久久精品亚洲午夜| 国产精品无大码| 精品午夜福利在线看| 久久久久久久精品精品| 久久精品国产亚洲av天美| 另类亚洲欧美激情| 我的老师免费观看完整版| 国产精品成人在线| a级毛片免费高清观看在线播放| 久久综合国产亚洲精品| 国产精品99久久99久久久不卡 | 欧美高清性xxxxhd video| 黄色配什么色好看| 亚洲av成人精品一二三区| av又黄又爽大尺度在线免费看| 丰满少妇做爰视频| 日韩电影二区| 男女边摸边吃奶| 久久ye,这里只有精品| 久久久久性生活片| 亚洲无线观看免费| 日本午夜av视频| 亚洲欧美日韩卡通动漫| 内地一区二区视频在线| 永久网站在线| 丰满乱子伦码专区| 成人无遮挡网站| 日日啪夜夜爽| 高清午夜精品一区二区三区| 国产精品国产三级国产av玫瑰| 一级毛片我不卡| 国产大屁股一区二区在线视频| 大陆偷拍与自拍| 亚洲av男天堂| 中文乱码字字幕精品一区二区三区| 国产成人精品一,二区| 亚洲欧美精品自产自拍| a级毛片免费高清观看在线播放| 国产精品99久久99久久久不卡 | 亚洲人成网站高清观看| 九九在线视频观看精品| 最近最新中文字幕免费大全7| 一本久久精品| 午夜免费观看性视频| 夜夜看夜夜爽夜夜摸| 免费电影在线观看免费观看| 亚洲国产最新在线播放| 久久久久久久亚洲中文字幕| 国产爱豆传媒在线观看| 午夜亚洲福利在线播放| av线在线观看网站| 成人亚洲欧美一区二区av| 国产精品一及| 男人和女人高潮做爰伦理| 在线观看一区二区三区激情| 麻豆久久精品国产亚洲av| 777米奇影视久久| 最近手机中文字幕大全| 亚洲精华国产精华液的使用体验| 亚洲欧美一区二区三区黑人 | 成人亚洲欧美一区二区av| 国产有黄有色有爽视频| 久久久亚洲精品成人影院| 综合色av麻豆| 亚洲国产精品成人久久小说| 男男h啪啪无遮挡| 99热这里只有是精品在线观看| 久久久久久久国产电影| 国产日韩欧美亚洲二区| 国产精品麻豆人妻色哟哟久久| 成人欧美大片| 麻豆成人av视频| 女人久久www免费人成看片| 可以在线观看毛片的网站| 在线观看一区二区三区激情| 十八禁网站网址无遮挡 | 久久久午夜欧美精品| 日本黄大片高清| 精品国产三级普通话版| 深爱激情五月婷婷| 99久久精品国产国产毛片| 久久久久久久久久久丰满| 欧美xxⅹ黑人| av在线天堂中文字幕| 91精品国产九色| 欧美老熟妇乱子伦牲交| 国内精品美女久久久久久| 男女边吃奶边做爰视频| 欧美成人精品欧美一级黄| 女人久久www免费人成看片| 视频区图区小说| 禁无遮挡网站| 亚洲三级黄色毛片| 日韩av在线免费看完整版不卡| 噜噜噜噜噜久久久久久91| 亚洲精品国产av成人精品| av女优亚洲男人天堂| 久久久欧美国产精品| 亚洲精品,欧美精品| 在线观看免费高清a一片| 亚洲精品一区蜜桃| 日韩欧美 国产精品| 99九九线精品视频在线观看视频| 久热久热在线精品观看| 亚洲av二区三区四区| 亚洲欧洲国产日韩| 好男人在线观看高清免费视频| 国产午夜精品一二区理论片| 日韩精品有码人妻一区| 麻豆成人午夜福利视频| 天天一区二区日本电影三级| 国产一区亚洲一区在线观看| 美女高潮的动态| 亚洲最大成人中文| 国产成人福利小说| 成年av动漫网址| 成人亚洲欧美一区二区av| 一边亲一边摸免费视频| 亚洲国产成人一精品久久久| 我的老师免费观看完整版| 超碰av人人做人人爽久久| 国产乱人偷精品视频| 纵有疾风起免费观看全集完整版| 日日摸夜夜添夜夜爱| 青春草视频在线免费观看| 1000部很黄的大片| 日韩av在线免费看完整版不卡| 国产精品人妻久久久影院| 能在线免费看毛片的网站| 欧美变态另类bdsm刘玥| 国产精品.久久久| 男女边吃奶边做爰视频| av在线app专区| 肉色欧美久久久久久久蜜桃 | 秋霞伦理黄片| 五月伊人婷婷丁香| 乱系列少妇在线播放| 欧美激情国产日韩精品一区| 久久精品人妻少妇| 久热这里只有精品99| 亚洲欧美成人精品一区二区| 精品熟女少妇av免费看| 免费观看a级毛片全部| av免费在线看不卡| 综合色av麻豆| 亚洲欧美清纯卡通| 午夜精品国产一区二区电影 | 搡老乐熟女国产| 热99国产精品久久久久久7| 丝瓜视频免费看黄片| 看黄色毛片网站| 在线a可以看的网站| 亚洲精品日韩av片在线观看| 在线观看一区二区三区| 大片电影免费在线观看免费| 能在线免费看毛片的网站| 男人狂女人下面高潮的视频| 一级毛片aaaaaa免费看小| 三级经典国产精品| 久久精品人妻少妇| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 一边亲一边摸免费视频| 日韩成人av中文字幕在线观看| 欧美zozozo另类| 在线亚洲精品国产二区图片欧美 | 成人一区二区视频在线观看| 亚洲精品日韩在线中文字幕| 久久97久久精品| 国产色婷婷99| 免费高清在线观看视频在线观看| 美女内射精品一级片tv| 久久精品熟女亚洲av麻豆精品| 2021天堂中文幕一二区在线观| 欧美zozozo另类| 青春草视频在线免费观看| 日本色播在线视频| 久久久久网色| 最近最新中文字幕大全电影3| 日日摸夜夜添夜夜添av毛片| 少妇人妻精品综合一区二区| 国产亚洲一区二区精品| 男女下面进入的视频免费午夜| 一级a做视频免费观看| 大片免费播放器 马上看| 日日啪夜夜爽| 国产一区二区三区综合在线观看 | 熟女人妻精品中文字幕| 狠狠精品人妻久久久久久综合| 联通29元200g的流量卡| 国产精品不卡视频一区二区| 日韩欧美 国产精品| 国产真实伦视频高清在线观看| av在线天堂中文字幕| 国产 精品1| av在线老鸭窝| 国产熟女欧美一区二区| 中文欧美无线码| 有码 亚洲区| 简卡轻食公司| 久久女婷五月综合色啪小说 | 免费少妇av软件| av又黄又爽大尺度在线免费看| 久久久久久久亚洲中文字幕| 国产av不卡久久| 精品视频人人做人人爽| 精品一区在线观看国产| 午夜免费男女啪啪视频观看| 国产精品女同一区二区软件| 看免费成人av毛片| 女人十人毛片免费观看3o分钟| 九草在线视频观看| 久久韩国三级中文字幕| .国产精品久久| 欧美日韩国产mv在线观看视频 | 高清日韩中文字幕在线| 国产精品久久久久久久久免| 中文天堂在线官网| 99久久精品国产国产毛片| 国产探花极品一区二区| 国内揄拍国产精品人妻在线| 亚洲精品成人av观看孕妇| 国产精品一区二区在线观看99| 18禁在线播放成人免费| 亚洲人成网站在线观看播放| 黑人高潮一二区| 亚洲av不卡在线观看| 神马国产精品三级电影在线观看| 亚洲精品国产av成人精品| 岛国毛片在线播放| 综合色丁香网| 成人午夜精彩视频在线观看| 久久人人爽人人爽人人片va| 男女那种视频在线观看| 国内少妇人妻偷人精品xxx网站| 深爱激情五月婷婷| 美女被艹到高潮喷水动态| 国产精品av视频在线免费观看| 国产精品国产三级国产专区5o| 亚洲成色77777| 久久午夜福利片| 99精国产麻豆久久婷婷| 国产免费一区二区三区四区乱码| 亚洲在久久综合| 亚洲欧美日韩卡通动漫| 欧美+日韩+精品| 麻豆国产97在线/欧美| 狠狠精品人妻久久久久久综合| 久久久亚洲精品成人影院| 校园人妻丝袜中文字幕| 国产伦理片在线播放av一区| 91久久精品国产一区二区三区| 黄色欧美视频在线观看| 五月玫瑰六月丁香| av黄色大香蕉| 街头女战士在线观看网站| 日本黄色片子视频| 日韩欧美 国产精品| 日本一二三区视频观看| 成人午夜精彩视频在线观看| 国产亚洲最大av| 人人妻人人看人人澡| 久久久国产一区二区| 国产 一区 欧美 日韩| 亚洲美女视频黄频| 91精品国产九色| 免费看a级黄色片| 国产av国产精品国产| 赤兔流量卡办理| 少妇人妻 视频| 亚洲高清免费不卡视频| 日韩三级伦理在线观看| 深爱激情五月婷婷| 欧美3d第一页| 国产熟女欧美一区二区| 国产精品精品国产色婷婷| 涩涩av久久男人的天堂| 亚洲av成人精品一二三区| 久久99热6这里只有精品| 欧美精品人与动牲交sv欧美| 午夜老司机福利剧场| www.av在线官网国产| 色视频在线一区二区三区| 午夜免费男女啪啪视频观看| 91久久精品国产一区二区三区| 乱码一卡2卡4卡精品| 久久精品久久久久久噜噜老黄| 久久久久性生活片| 久久久午夜欧美精品| 久久国内精品自在自线图片| 人妻制服诱惑在线中文字幕| 观看美女的网站| 国产午夜精品一二区理论片| 丰满乱子伦码专区| 天堂网av新在线| 99久久中文字幕三级久久日本| 亚洲av不卡在线观看| 亚洲精品,欧美精品| 99热6这里只有精品| 亚洲av国产av综合av卡| 最近2019中文字幕mv第一页| 我的老师免费观看完整版| 2021少妇久久久久久久久久久| av天堂中文字幕网| 一二三四中文在线观看免费高清| 亚洲最大成人中文| 青青草视频在线视频观看| 91精品一卡2卡3卡4卡| 一级毛片电影观看| 久久这里有精品视频免费| 3wmmmm亚洲av在线观看| 国产精品一二三区在线看| 91久久精品电影网| 超碰av人人做人人爽久久| 久热这里只有精品99| 最近最新中文字幕免费大全7| 国产黄色视频一区二区在线观看| 日本-黄色视频高清免费观看| 好男人视频免费观看在线| 久久久久久久久久久免费av| 色播亚洲综合网| 亚洲在久久综合| 有码 亚洲区| 如何舔出高潮| 国产视频内射| 亚洲美女搞黄在线观看| 国产探花极品一区二区| 少妇被粗大猛烈的视频| 大片免费播放器 马上看| 国产免费又黄又爽又色| 日韩欧美 国产精品| 久久久久久久午夜电影| 国产成人aa在线观看| 久久久久国产网址| 18禁裸乳无遮挡动漫免费视频 | 亚洲av免费高清在线观看| 亚洲精品一二三| 亚洲av国产av综合av卡| 久久久久久久午夜电影| 大香蕉97超碰在线| 国产探花极品一区二区| 亚洲国产欧美在线一区| 免费看日本二区| 日韩成人伦理影院| 国产精品一及| 欧美3d第一页| 成年女人看的毛片在线观看| 青青草视频在线视频观看| 国产日韩欧美亚洲二区| 一个人看视频在线观看www免费| 日韩三级伦理在线观看| 欧美+日韩+精品| 男插女下体视频免费在线播放| 亚洲精品日韩在线中文字幕| 亚洲久久久久久中文字幕| 国产黄a三级三级三级人| 亚洲性久久影院| 最近手机中文字幕大全| 国产午夜精品久久久久久一区二区三区| 亚洲天堂国产精品一区在线| 狂野欧美激情性xxxx在线观看| 久久女婷五月综合色啪小说 | 成人亚洲精品一区在线观看 | 成人二区视频| 中国三级夫妇交换| 久久久久久久久久成人| 春色校园在线视频观看| 免费看a级黄色片| 国产 一区 欧美 日韩| 偷拍熟女少妇极品色| 久久久久久久久久人人人人人人| 久久久久性生活片| 中国三级夫妇交换| 亚洲性久久影院| 毛片女人毛片| 久久99精品国语久久久| 国产精品蜜桃在线观看| 亚洲国产精品专区欧美| 成人美女网站在线观看视频| 国产亚洲91精品色在线| 国产又色又爽无遮挡免| 日本爱情动作片www.在线观看| 亚洲国产日韩一区二区| 欧美成人一区二区免费高清观看| 亚洲欧美一区二区三区黑人 | 三级国产精品欧美在线观看| 亚洲精品日本国产第一区| 欧美精品人与动牲交sv欧美| 色视频www国产| 最新中文字幕久久久久| 日日啪夜夜爽| 中文精品一卡2卡3卡4更新| 赤兔流量卡办理| 综合色丁香网| 成年版毛片免费区| 一区二区三区乱码不卡18| 精品熟女少妇av免费看| 狂野欧美白嫩少妇大欣赏| 一区二区三区精品91| 国产欧美日韩一区二区三区在线 | 成人无遮挡网站| 免费观看av网站的网址| 国产黄色视频一区二区在线观看| 亚洲欧美精品专区久久| 国产91av在线免费观看| 亚洲av.av天堂| 欧美一级a爱片免费观看看| 成人特级av手机在线观看| 亚洲国产精品999| 啦啦啦中文免费视频观看日本| 亚洲欧美日韩东京热| 国产亚洲最大av| 精品视频人人做人人爽| 在线观看av片永久免费下载| 亚洲成色77777| 中国三级夫妇交换| 国产爽快片一区二区三区| 日本-黄色视频高清免费观看| 九色成人免费人妻av| 国产日韩欧美亚洲二区| 亚洲av免费在线观看| 久久久精品免费免费高清| 一级毛片久久久久久久久女| 免费大片黄手机在线观看| 亚洲va在线va天堂va国产| 两个人的视频大全免费| 免费黄色在线免费观看| 亚洲精品日韩av片在线观看| 久久精品国产鲁丝片午夜精品| 中国三级夫妇交换| 秋霞伦理黄片| 国产精品女同一区二区软件| 亚洲最大成人中文| 性色avwww在线观看| 国内少妇人妻偷人精品xxx网站| 国产乱来视频区| eeuss影院久久| 在线观看人妻少妇| 精品久久久久久久久av| 久久久久久久久久成人| 久久午夜福利片| 亚洲欧美一区二区三区黑人 | av专区在线播放| 能在线免费看毛片的网站| 免费观看a级毛片全部|