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

    不同生態(tài)基流在水庫優(yōu)化調度中的合理性分析

    2018-09-10 20:05:34吳承君吳丹暉黃顯峰
    人民黃河 2018年11期
    關鍵詞:粒子群算法發(fā)電量水電站

    吳承君 吳丹暉 黃顯峰

    摘要:針對不同來水期生態(tài)需水的差異,提出改進最小月平均徑流法計算不同來水期的生態(tài)基流。通過建立基于能值分析的不同生態(tài)基流約束下水電站經(jīng)濟運行模型,以發(fā)電量最大為目標,根據(jù)水庫下游需水是否遭受破壞,提出了需水破壞系數(shù)的概念,基于水電站調度結果及需水破壞系數(shù),利用基于能值分析的不同生態(tài)基流約束下水電站運行評價函數(shù),從定性和定量兩方面分析了不同生態(tài)基流的合理性及可靠性。陸渾水庫實例計算表明:利用改進最小月平均徑流法計算的生態(tài)基流需水破壞系數(shù)為0.19,能值評價函數(shù)值為4.42×1020sej ,貨幣價值評價函數(shù)值為4.65億元;改進最小月平均徑流法相比常用水文學生態(tài)基流計算方法具有更好的適用性,以及針對年內不同來水期的適應性。

    關鍵詞:水電站;發(fā)電量;生態(tài)基流;能值分析;粒子群算法

    中圖分類號:TV213.4 文獻標志碼:A

    隨著對生態(tài)環(huán)境保護重視程度的不斷提高,河道生態(tài)基流成為研究的熱點。目前,關于河道生態(tài)基流的確定方法有水文學法、水力學法、水文一生物分析法等,其中水文學法計算方便、所需資料易于獲取,是目前最常用的生態(tài)基流計算方法,但其結果通常為單一流量值,對于年內來水不均的河流適用性不強。針對生態(tài)基流計算方法的多樣性,諸多學者[1-5]對不同生態(tài)基流確定方法的計算結果進行了對比分析,評價各方法的合理性和適用性,但主要是對生態(tài)基流計算結果的對比,缺乏對不同生態(tài)基流產生經(jīng)濟效益的定量評價。

    能值分析理論由Odum于1987年提出,其以太陽能為單位,將不同形式的物質或能量轉化為統(tǒng)一的衡量標準,可直觀反映不同要素之間的關系[6]。呂翠美[7]首次將能值分析理論應用到水資源生態(tài)經(jīng)濟價值計算中;李懷恩等[8]利用能值分析理論,探討了河道生態(tài)基流價值在年內的變化情況;賀成龍[9-10]研究了中國水力發(fā)電的能值轉換率,并利用能值足跡模型對三峽工程建設期和運行期的生態(tài)效益進行了定量評價;李友輝等[11]利用能值理論分析了水資源的投入產出;田桂桂等[12]利用能值理論對生態(tài)水系的生態(tài)環(huán)境效益進行了定量分析和統(tǒng)一度量。根據(jù)已有研究成果,針對水文學法生態(tài)基流計算過程中存在的問題,本文提出了改進的最小月平均徑流法,為評估利用改進的最小月平均徑流法計算的生態(tài)基流結果的合理性和適用性,建立基于能值分析的不同生態(tài)基流約束下水電站經(jīng)濟運行模型,利用能值進行統(tǒng)一度量,結合水電站調度結果,從定性及定量兩方面分析不同生態(tài)基流計算結果的合理性和可靠性。

    1 基于能值分析的不同生態(tài)基流約束下水電站經(jīng)濟運行模型

    電力市場環(huán)境下,滿足合約電量要求、提高現(xiàn)貨交易電量是水電站經(jīng)濟運行的主要途徑[13]。以不同生態(tài)基流約束下水電站發(fā)電量最大為目標函數(shù),建立基于能值分析的不同生態(tài)基流約束下水電站經(jīng)濟運行模型。

    1.1 目標函數(shù)及評價函數(shù)

    以水電站發(fā)電量最大為目標,目標函數(shù)為

    基于能值分析的不同生態(tài)基流約束下水電站運行的評價函數(shù)為式中:max E為調度期內最大發(fā)電量;Nt為t時段的出力;max En為第n種生態(tài)基流約束下的最大發(fā)電量;ηE為水電站發(fā)電能值轉換率;P為電價;ηB為價格能值轉換率;k為水庫所在地能值貨幣比率;Qni,為第n種生態(tài)基流在i月的流量;e為水電站所在地單方水能值;△t為常數(shù),其值為2.6×106;ε為水庫下游需水破壞系數(shù),表示水庫下游用水戶受影響遭破壞程度,(n'為水庫調度期間下游用水戶遭遇破壞的年數(shù);N'為水庫調度總年數(shù)),ε∈(0~1)。

    由于來水的不確定性,當年內來水較少或年內來水嚴重不均無法或部分時段難以滿足下游需水要求時,水電站經(jīng)濟運行將受到嚴重限制,因此提出需水破壞系數(shù)的概念,需水破壞系數(shù)反映水庫在多年運行調度中供水無法滿足水庫下游用水需求出現(xiàn)的頻率。

    1.2 約束條件

    模型約束條件主要為水庫調度約束。水量平衡約束為

    Vt+1=Vt+(It-Ot)Δt(3)式中:Vt、Vt+1分別為t時段初和t時段末的蓄水量;It、Ot分別為t時段的入庫流量和出庫流量。

    水庫蓄水位約束為

    Ztmin≤Zt≤Ztmax(4)式中:Ztmin、Ztmax分別為t時段庫水位最低、最高約束。出庫流量約束為式中:Qtmin、Qtmax分別為t時段水庫最小、最大下泄流量,其中Qtmin為水庫下游河道生態(tài)基流Qt,st、灌溉需水流量Qt,gg等各種需水流量總和;Qt為t時段水庫下泄流量。

    電站出力約束為

    Nmin≤Nt≤Nmax(6)式中:Nmin、Nmax分別為電站允許出力下限、上限。

    1.3 能值分析

    能值分析是對可更新資源、不可更新資源等的系統(tǒng)分析,產品產出包括正效應產出和負效應產出兩大類。其中,正效應產出是對社會、生態(tài)或經(jīng)濟發(fā)展具有積極影響的產出,其他則為負效應產出。能值轉換率和能值貨幣比率是能值價值分析的基礎,根據(jù)國民經(jīng)濟和社會發(fā)展情況,利用相應的能值轉換率分析計算國民經(jīng)濟及社會發(fā)展總能值及能值貨幣比率。能值轉換率及能值貨幣比率為式中:ηi為第i種物質的能值轉換率;Ji為形成第a種物質需使用的太陽能;Mi為第i種物質的質量或能量;k為能值貨幣比率;PG為國家(或區(qū)域)國民生產總值。

    根據(jù)式(7)和式(8),利用能值轉換率和能值貨幣比率可以反推不同物質具有的能值及貨幣價值。

    2 改進最小月平均徑流法

    針對水文學法計算生態(tài)基流時未考慮或未全面考慮不同來水期生態(tài)需水不同的問題,提出改進最小月平均徑流法,利用以發(fā)電量最大為目標的水電站經(jīng)濟運行模型及能值分析理論,分析評價不同生態(tài)基流的合理性和適用性。根據(jù)長系列入庫徑流資料,通過劃分不同類型(豐、平、枯)的來水年,利用改進的最小月平均徑流法確定不同來水期河道的生態(tài)基流。

    (1)來水年型劃分。目前,常用的來水年型劃分方法包括保證率劃分法和距平百分率法[14]。由于距平百分率法不受系列個數(shù)的影響,劃分更為客觀[15],因此采用距平百分率法劃分來水年型,其中豐水年距平百分率μ>20%、平水年-20%<μ≤20%、枯水年μ≤-20%。距平百分率μ為

    μ=(Qi-Qa)/Qa×100%(9)式中:Qi為第i年平均徑流量,m3/s;Qa為多年平均徑流量,m3/s。

    (2)改進最小月平均徑流法。Tennant法和最小月平均徑流法作為確定生態(tài)基流最常用的兩種方法,其結果通常為單一流量值,不適用于季節(jié)性河流。通常情況下,河流分為豐水期、平水期和枯水期,豐水期生態(tài)需水量較大,枯水期生態(tài)需水量相對較小。本文根據(jù)最小月平均徑流法原理,綜合考慮不同水期生態(tài)基流存在差異的特點,根據(jù)劃分的豐水年、平水年和枯水年,分別計算相應來水期下河道的生態(tài)基流Qj:式中:Qj(j=1,2,3)分別為豐水期、平水期、枯水期下的生態(tài)基流量,m3/s;Qij表示j年型下第i年的最小月平均流量,m3/s;nj表示j年型下的年數(shù)。

    3 模型求解算法

    粒子群算法具有全局尋優(yōu)性、魯棒性、并行性及編程易實現(xiàn)性等特點,是目前水庫調度中應用比較成熟的一種智能優(yōu)化算法[16-17],因此采用粒子群算法求解基于能值分析的不同生態(tài)基流約束下的水電站經(jīng)濟運行模型。由于粒子群算法在生成初始解時未考慮水庫調度的多約束性,因此絕大多數(shù)初始解質量較低[18]。筆者采用罰函數(shù)法解決水電站經(jīng)濟運行模型中的強約束問題,詳見文獻[19]。為提高粒子群算法前期的全局尋優(yōu)能力及后期局部尋優(yōu)能力,采用隨迭代次數(shù)的增加,慣性權重線性遞減的處理方法。式中:i為當前迭代次數(shù);m為最大迭代次數(shù);ω為慣性權重;ωmax為慣性權重系數(shù)的最大值;ωmin為慣性權重系數(shù)的最小值。

    4 實例分析

    以河南省洛陽市篙縣陸渾水庫為研究對象進行分析。陸渾水庫位于黃河二級支流伊河上,在防洪、發(fā)電、灌溉和水產養(yǎng)殖等方面發(fā)揮了重要作用,受溫帶季風氣候影響,水庫所在流域降水量隨季節(jié)變化較大,汛期為6-9月。

    4.1 生態(tài)基流計算

    根據(jù)陸渾水庫1980-2010年共31 a的長系列來水資料,通過劃分來水年型,分別采用Tennant法、最小月平均徑流法和改進最小月平均徑流法計算生態(tài)基流。Tennant法計算結果見表1。最小月平均徑流法計算得出生態(tài)基流為4.12m3/s,年徑流量為1.300億m3。改進最小月平均徑流法計算得到豐水期、平水期、枯水期生態(tài)基流分別為5.83、3.33、3.25m3/s,年徑流量為1.370億m3??梢钥闯?,改進最小月平均徑流法、最小月平均徑流法及Tennant法中“一般”標準下的生態(tài)基流在年內不同時期雖有差異,但年徑流量總體相差不大,表明改進最小月平均徑流法計算得到的生態(tài)基流是可靠的。Tennant法中“最小”標準下的生態(tài)基流明顯低于最小月平均徑流法和改進最小月平均徑流法計算的生態(tài)基流,“良好”“好”及“很好”標準下的年徑流量明顯較其他方法計算的年徑流量大。

    4.2 能值計算

    (1)篙縣能值貨幣比率。由于未搜集到篙縣2016年國民經(jīng)濟和社會發(fā)展公報,根據(jù)式(2),考慮到能值貨幣比率的改變對研究結論無影響,因此將洛陽市2016年能值貨幣比率近似看作篙縣2016年能值貨幣比率。洛陽市2016年能值分析見表2,其中,“可更新資源”各形式能量原始數(shù)據(jù)的確定方法見文獻[20-21],其余原始數(shù)據(jù)主要由洛陽市2016年國民經(jīng)濟和社會發(fā)展統(tǒng)計公報獲得。

    根據(jù)式(8),由表2計算得到洛陽市2016年能值貨幣比率為9.50×1011sej/元,即6.28×1012sej/美元(按2016年6月人民幣對美元匯率6.61計),其值低于我國(19.12×1012sej/美元)平均值,高于北京(0.51×1012sej/美元)和發(fā)達國家如美國(2.55×1012sej/美元)平均值。

    (2)篙縣水體服務功能能值價值計算。水體服務功能能值價值計算主要包括休閑娛樂能值價值、維持生物多樣性能值價值、氣候調節(jié)能值價值、水質凈化能值價值及科學研究能值價值。根據(jù)文獻[8],水體服務功能能值價值主要體現(xiàn)在前兩種能值價值方面,后三種占比相對較小,因此僅對水體服務功能的休閑娛樂能值價值及維持生物多樣性能值價值進行分析。

    休閑娛樂能值價值計算主要為旅游收入的計算,根據(jù)篙縣富水旅游的特點,假定篙縣旅游收入60%(根據(jù)篙縣水景觀或水上旅游景點與篙縣總景點的比例估算)吸引力來自水體。根據(jù)篙縣和洛陽市2015年國民經(jīng)濟和社會發(fā)展公報,其2015年旅游收入分別為33.4億元、780.0億元,假設嵩縣和洛陽市的旅游收入同比增長,根據(jù)洛陽市2016年旅游收入,估算得篙縣2016年旅游收入為38.75億元,利用洛陽市2016年能值貨幣比率,根據(jù)式(8)計算得篙縣水體休閑娛樂能值價值:38.75×108×60%×9.50×1011=2.21×1021sej。

    維持生物多樣性能值價值主要針對區(qū)域內水生生物進行計算,計算公式為

    φ=ργF0/Fe(12)式中:φ為維持生物多樣性能值價值;P為平均每個物種的太陽能值(1.26×1026sej/種[8]);γ為區(qū)域內水生生物物種數(shù);F0為水生生物活動面積;Fe為地球表面積(5.21×1014m2)。

    在維持生物多樣性能值價值分析計算時,考慮區(qū)域物種的相似性,根據(jù)伊河、伊河支流沙河生物群落研究成果,流域內鑒定出魚類50種、浮游植物105種、底棲動物95種,共有水生生物物種250種。篙縣水面面積約為3.33×107m2,根據(jù)式(12)計算得到維持生物多樣性能值價值為2.01×1021sej。

    根據(jù)篙縣水資源總量(11億m3)計算得篙縣單方水能值為(2.21×1021+2.02×1021)/(11×108)=3.84×1012sej,從而單方水貨幣價值為(3.84×1012)/(9.50×1011)=4.04元。根據(jù)單方水能值及單方水能值價值計算不同生態(tài)基流能值及其價值,見表3。由表3可以看出,“Tennant法很好”生態(tài)基流能值及其貨幣價值最大,“Tennant法最小”生態(tài)基流能值及貨幣價值最小,生態(tài)基流能值及其貨幣價值與生態(tài)基流年徑流量之間存在明顯正相關關系。

    4.3 不同生態(tài)基流約束下水電站運行調度

    陸渾水庫水電站出力系數(shù)設為8.0,初始化粒子群個數(shù)為100,迭代200次,粒子移動最大速度為1.5,權重系數(shù)最大值為1.2,最小值為0.8。針對不同生態(tài)基流約束條件下的水電站優(yōu)化調度程序分別運行20次,取最大值作為不同生態(tài)基流約束下的水電站最大發(fā)電量。根據(jù)來水年份發(fā)電量距平百分率由大到小排序,進行各年不同來水、不同生態(tài)基流約束下的水電站調度,結果見圖1和表4。

    由圖1可以看出:不同生態(tài)基流約束下,其水電站發(fā)電量調度結果均隨年內來水量的減小整體呈下降趨勢;不同生態(tài)基流約束對水電站豐水年型發(fā)電量的影響無明顯差異,對水電站平水年和枯水年型發(fā)電量的影響較大,整體表現(xiàn)為隨著生態(tài)基流量的增大,發(fā)電量呈減小趨勢。由表4可以看出,隨著生態(tài)基流量的增大,水庫年內來水不能滿足水庫下游需水要求的情況增多,需水破壞系數(shù)增大,即生態(tài)流量并非越大越好,在一定程度上,隨著生態(tài)流量的增大,水庫供水無法滿足下游需水的概率增大,多年平均發(fā)電量整體表現(xiàn)為降低趨勢。從中長期角度分析,生態(tài)基流與水電站發(fā)電之間存在明顯的反比關系,即生態(tài)基流能值(或貨幣價值)越大,水電站發(fā)電量越小,水電站的實際經(jīng)濟運行效益越低。值得指出的是由于1986年來水較枯,因此無論生態(tài)基流大小如何,水庫下泄流量均不能滿足下游需水要求。

    兼顧生態(tài)環(huán)境保護及水電站經(jīng)濟運行效益,生態(tài)基流的確定既不應過大也不能太小,從定性分析角度考慮,利用“Tennant法一般”“最小月平均徑流法”和“改進最小月平均徑流法”計算的生態(tài)基流更滿足決策者偏好。

    4.4 生態(tài)基流綜合評價

    根據(jù)發(fā)電量調度結果,計算發(fā)電量能值價值。根據(jù)文獻[10],發(fā)電量能值轉換率為5.69×1011sej/(kW·h)??紤]河南省不同用戶及階梯電價標準,取電價為0.6元/(kW·h)。發(fā)電量能值價值見表5。結合生態(tài)基流能值價值及發(fā)電量能值價值,考慮需水破壞系數(shù)影響,根據(jù)式(2)計算不同生態(tài)基流約束下的總能值貨幣價值,見表6。

    由表6可以看出,考慮需水破壞系數(shù)后,“Tennant法良好”能值最大,貨幣價值最大,其次為“改進最小月平均徑流法”,“Tennant法很好”能值及貨幣價值最小。從需水破壞系數(shù)角度考慮,“Tennant法良好”“Tennant法好”和“Tennant法很好”三種生態(tài)基流約束條件下,需水破壞系數(shù)較大,水庫來水通常無法滿足該生態(tài)基流條件下的下游需水要求。結合定性分析及定量計算結果,本文提出的改進最小月平均徑流法更具合理性和可靠性。

    5 結論

    根據(jù)年內來水不均的特點,提出改進最小月平均徑流法計算生態(tài)基流,建立了基于能值分析的不同生態(tài)基流約束下水電站經(jīng)濟運行模型,利用水電站發(fā)電量及能值分析理論從定性和定量兩方面分析不同生態(tài)基流的合理性和可靠性。實例研究表明,改進最小月平均徑流法相比常用水文學生態(tài)基流計算方法具有更好的適用性,以及針對年內不同來水期的適應性,體現(xiàn)了改進最小月平均徑流法計算生態(tài)基流的合理性和可靠性。

    參考文獻:

    [1]許亞雄,趙志懷.基于改進Tennant法的河流生態(tài)基流估算[J].水電能源科學,2017,35(2):38-40.

    [2]徐宗學,武瑋,于松延.生態(tài)基流研究:進展與挑戰(zhàn)[J].水力發(fā)電學報,2016,35(4):1-11.

    [3]商玲,李宗禮,孫偉,等.基于HIMS模型的西營河流域河道內生態(tài)基流估算[J].水土保持研究,2014,21(1):100-103.

    [4]肖玉成,董飛,張新華,等.基于SWAT分布式水文模型的河道內生態(tài)基流[J].四川大學學報(工程科學版),2013,45(1):85-90.

    [5]張帆,徐建新,郭文獻.河流生態(tài)環(huán)境需水量計算及優(yōu)化分析模型應用[J].人民黃河,2011,33(5):72-73,129.

    [6]ODUM H T,ODUM E C,BLISSETY M.Ecology and Econ-omy:Energy Analysis and Public Policy in Texas Results ofPolicy Research Project[R].Austin Texas:LBJ School ofPublic Affairs,1987:421.

    [7]呂翠美.區(qū)域水資源生態(tài)經(jīng)濟價值的能值研究[D].鄭州:鄭州大學,2009:36.

    [8]李懷恩,徐梅梅,趙宇.基于能值法的河道生態(tài)基流價值及其年內時間變化研究[J].生態(tài)經(jīng)濟,2017,33(3):160-164.

    [9]賀成龍.三峽工程的能值足跡與生態(tài)承載力[J].自然資源學報,2017,32(2):329-341.

    [10]賀成龍.水力發(fā)電的能值轉換率計算方法[J].自然資源學報,2016,31(11):1958-1968.

    [11]李友輝,段曹斌,孔瓊菊.基于能值理論的水資源投入產出分析[J].人民黃河,2013,35(12):57-60.

    [12]田桂桂,吳澤寧,郭溪.生態(tài)水系生態(tài)環(huán)境效益能值評估方法及其應用[J].人民黃河,2014,36(8):76-78,82.

    [13]原文林,吳澤寧,黃強.電力市場環(huán)境下梯級水庫發(fā)電優(yōu)化調度的協(xié)進化粒子群算法應用研究[J].水力發(fā)電學報,2011,30(3):65-70.

    [14]水利部水利信息中心.水文情報預報規(guī)范:GB/T22482-2008[S].北京:中國標準出版社,2008:10.

    [15]吳喜軍,李懷恩,董穎,等.基于基流比例法的渭河生態(tài)基流計算[J].農業(yè)工程學報,2011,27(10):154-159.

    [16]鄒強,王學敏,李安強,等.基于并行混沌量子粒子群算法的梯級水庫群防洪優(yōu)化調度研究[J].水利學報,2016,47(8):967-976.

    [17]周建中,李英海,肖舸,等.基于混合粒子群算法的梯級水電站多目標優(yōu)化調度[J].水利學報,2010,41(10):1212-1219.

    [18]鄧顯羽,彭勇,葉碎高,等.粒子群算法在水庫(群)優(yōu)化調度研究中的應用綜述[J].水利水電科技進展,2010,30(5):90-94.

    [19]馬細霞,儲冬冬.粒子群優(yōu)化算法在水庫調度中的應用分析[J].鄭州大學學報(工學版),2006,27(4):121-124.

    [20]齊雪艷.東居延海濕地生態(tài)經(jīng)濟系統(tǒng)價值的能值分析[D].鄭州:鄭州大學,2013:37-39.

    [21]田桂桂.基于物質循環(huán)的生態(tài)用水價值能值評估方法研究[D].鄭州:鄭州大學,2016:32-37.

    猜你喜歡
    粒子群算法發(fā)電量水電站
    里底水電站工程
    西北水電(2022年1期)2022-04-16 05:20:06
    9月份全國發(fā)電量同比下降0.4%
    白鶴灘水電站有多“?!?
    金橋(2021年8期)2021-08-23 01:06:54
    2019年全國發(fā)電量同比增長3.5%
    中國煤炭(2020年2期)2020-01-21 02:49:28
    全國上半年水力發(fā)電量同比增長2.9%
    水電站設計(2018年3期)2018-03-26 03:52:32
    小型水電站技術改造探討
    電力市場交易背景下水電站優(yōu)化調度研究
    基于粒子群算法的產業(yè)技術創(chuàng)新生態(tài)系統(tǒng)運行穩(wěn)定性組合評價研究
    預測(2016年5期)2016-12-26 10:04:59
    肯尼亞火電發(fā)電量下降
    電站輔機(2016年4期)2016-05-17 03:52:38
    交通堵塞擾動下多車場車輛路徑優(yōu)化
    商(2016年5期)2016-03-28 18:10:26
    国产免费福利视频在线观看| 国产深夜福利视频在线观看| 日韩一卡2卡3卡4卡2021年| 国产有黄有色有爽视频| 少妇的逼水好多| 91成人精品电影| 国产老妇伦熟女老妇高清| 国产成人免费无遮挡视频| 国产一区亚洲一区在线观看| 丝袜人妻中文字幕| 久久女婷五月综合色啪小说| 国产不卡av网站在线观看| 国产精品蜜桃在线观看| 在现免费观看毛片| 日日爽夜夜爽网站| 欧美少妇被猛烈插入视频| 日产精品乱码卡一卡2卡三| 女人久久www免费人成看片| 国产精品香港三级国产av潘金莲 | 美女中出高潮动态图| 韩国高清视频一区二区三区| 日韩成人av中文字幕在线观看| 国产片内射在线| 中文字幕最新亚洲高清| www.熟女人妻精品国产| 有码 亚洲区| 国产精品人妻久久久影院| 亚洲av男天堂| 亚洲一级一片aⅴ在线观看| 亚洲综合色惰| 日韩一卡2卡3卡4卡2021年| 精品久久久久久电影网| 日本午夜av视频| 日本色播在线视频| 欧美日韩精品成人综合77777| 日韩伦理黄色片| 制服诱惑二区| 母亲3免费完整高清在线观看 | 如日韩欧美国产精品一区二区三区| 亚洲av成人精品一二三区| 看非洲黑人一级黄片| 亚洲精品自拍成人| 欧美激情高清一区二区三区 | 国产在线免费精品| 国产又爽黄色视频| 一个人免费看片子| tube8黄色片| 午夜福利乱码中文字幕| 国产野战对白在线观看| 香蕉精品网在线| 国产乱来视频区| 亚洲成人一二三区av| 春色校园在线视频观看| 中文字幕av电影在线播放| 日韩中文字幕视频在线看片| 少妇被粗大猛烈的视频| 婷婷色麻豆天堂久久| 久久女婷五月综合色啪小说| 国产成人精品福利久久| 久久亚洲国产成人精品v| 午夜福利乱码中文字幕| 亚洲三级黄色毛片| 一个人免费看片子| 日日啪夜夜爽| 极品少妇高潮喷水抽搐| 熟妇人妻不卡中文字幕| 丝瓜视频免费看黄片| 十八禁高潮呻吟视频| 午夜福利一区二区在线看| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | videos熟女内射| 巨乳人妻的诱惑在线观看| 涩涩av久久男人的天堂| 精品人妻在线不人妻| 免费在线观看完整版高清| 欧美人与善性xxx| 国产探花极品一区二区| 中文字幕人妻丝袜一区二区 | 国产色婷婷99| 国产精品麻豆人妻色哟哟久久| 黑人巨大精品欧美一区二区蜜桃| 香蕉精品网在线| 如何舔出高潮| 免费看av在线观看网站| 黄网站色视频无遮挡免费观看| 水蜜桃什么品种好| 热99国产精品久久久久久7| 99久久综合免费| 久久 成人 亚洲| 精品少妇久久久久久888优播| 日本av手机在线免费观看| 国精品久久久久久国模美| 99久国产av精品国产电影| 免费看av在线观看网站| 国产免费又黄又爽又色| 国产在线视频一区二区| 90打野战视频偷拍视频| 亚洲中文av在线| 免费大片黄手机在线观看| www.自偷自拍.com| av电影中文网址| 天堂8中文在线网| av网站在线播放免费| 可以免费在线观看a视频的电影网站 | 亚洲欧美一区二区三区久久| 咕卡用的链子| 亚洲av免费高清在线观看| 欧美精品亚洲一区二区| 国语对白做爰xxxⅹ性视频网站| 人人澡人人妻人| 91在线精品国自产拍蜜月| 国产1区2区3区精品| av卡一久久| 男女国产视频网站| 成人18禁高潮啪啪吃奶动态图| 久久热在线av| 国产精品熟女久久久久浪| 亚洲精品成人av观看孕妇| 人妻系列 视频| 国产片特级美女逼逼视频| 777久久人妻少妇嫩草av网站| 下体分泌物呈黄色| 亚洲欧美一区二区三区久久| 亚洲欧美一区二区三区久久| 亚洲成人手机| 精品少妇内射三级| 久久国产亚洲av麻豆专区| 久久 成人 亚洲| 亚洲成人手机| 色播在线永久视频| 美国免费a级毛片| 女性被躁到高潮视频| 中文字幕制服av| 天天影视国产精品| 亚洲图色成人| av在线app专区| 高清在线视频一区二区三区| 伊人久久国产一区二区| 三上悠亚av全集在线观看| 三上悠亚av全集在线观看| 在线亚洲精品国产二区图片欧美| av不卡在线播放| 国产在视频线精品| av女优亚洲男人天堂| 曰老女人黄片| 男人爽女人下面视频在线观看| 亚洲精品成人av观看孕妇| 久久久久精品性色| 少妇人妻精品综合一区二区| 日韩中文字幕欧美一区二区 | 大陆偷拍与自拍| 久久久久久久久久人人人人人人| 丝袜美腿诱惑在线| 欧美xxⅹ黑人| 寂寞人妻少妇视频99o| 在线看a的网站| 母亲3免费完整高清在线观看 | 大码成人一级视频| 91久久精品国产一区二区三区| 午夜久久久在线观看| 亚洲一级一片aⅴ在线观看| 免费不卡的大黄色大毛片视频在线观看| 日本vs欧美在线观看视频| 亚洲一级一片aⅴ在线观看| 黑人猛操日本美女一级片| 国产片内射在线| 男男h啪啪无遮挡| 国产av国产精品国产| 国产1区2区3区精品| 美国免费a级毛片| 国产片特级美女逼逼视频| 2021少妇久久久久久久久久久| 五月天丁香电影| 韩国av在线不卡| 欧美 日韩 精品 国产| 男女午夜视频在线观看| av在线老鸭窝| 尾随美女入室| 欧美+日韩+精品| 波野结衣二区三区在线| 免费观看在线日韩| 大话2 男鬼变身卡| 亚洲精品自拍成人| 波多野结衣一区麻豆| 极品人妻少妇av视频| 免费不卡的大黄色大毛片视频在线观看| 免费久久久久久久精品成人欧美视频| 久久ye,这里只有精品| 9色porny在线观看| 国产免费视频播放在线视频| 欧美日韩国产mv在线观看视频| 欧美人与性动交α欧美软件| 人妻少妇偷人精品九色| 性色av一级| 亚洲美女搞黄在线观看| 国产亚洲精品第一综合不卡| 日韩视频在线欧美| 九色亚洲精品在线播放| 婷婷成人精品国产| 国产乱人偷精品视频| 高清视频免费观看一区二区| 亚洲精品美女久久久久99蜜臀 | av天堂久久9| 久久精品久久精品一区二区三区| 亚洲综合精品二区| 国产精品免费视频内射| 女的被弄到高潮叫床怎么办| 丁香六月天网| 久久精品aⅴ一区二区三区四区 | 男女啪啪激烈高潮av片| 亚洲av国产av综合av卡| 少妇被粗大猛烈的视频| 少妇的丰满在线观看| 91久久精品国产一区二区三区| 久久久久久久精品精品| 久久午夜福利片| 成年人午夜在线观看视频| 国产日韩欧美视频二区| 黑人巨大精品欧美一区二区蜜桃| 国产亚洲av片在线观看秒播厂| 久久毛片免费看一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 国产在线视频一区二区| 尾随美女入室| 天天躁夜夜躁狠狠久久av| 亚洲国产av新网站| 久久久久久久久免费视频了| 夜夜骑夜夜射夜夜干| 久热久热在线精品观看| 亚洲精品久久成人aⅴ小说| 在线天堂中文资源库| 免费高清在线观看日韩| 亚洲av国产av综合av卡| 久久精品国产自在天天线| 国产男人的电影天堂91| 天美传媒精品一区二区| 亚洲欧美清纯卡通| 免费看av在线观看网站| 丁香六月天网| av有码第一页| 精品福利永久在线观看| 久久久久视频综合| av在线播放精品| 一边摸一边做爽爽视频免费| 天天影视国产精品| 国产欧美亚洲国产| av网站免费在线观看视频| 亚洲av男天堂| 两个人看的免费小视频| 国产精品久久久久成人av| 国产精品久久久久久精品古装| 青春草视频在线免费观看| kizo精华| www.av在线官网国产| 久久99一区二区三区| 亚洲精品久久久久久婷婷小说| 精品久久久久久电影网| 青春草亚洲视频在线观看| 美女福利国产在线| 亚洲五月色婷婷综合| av在线app专区| 精品少妇黑人巨大在线播放| 国产成人精品久久久久久| 制服诱惑二区| 成年女人毛片免费观看观看9 | 一级片免费观看大全| 99热国产这里只有精品6| 国产亚洲最大av| 午夜老司机福利剧场| 国产精品一二三区在线看| 少妇精品久久久久久久| 免费大片黄手机在线观看| 成人黄色视频免费在线看| 亚洲精品国产av成人精品| 国产乱人偷精品视频| 国产不卡av网站在线观看| 美女高潮到喷水免费观看| 国产女主播在线喷水免费视频网站| 最近最新中文字幕免费大全7| 久久精品久久精品一区二区三区| 一二三四在线观看免费中文在| 老女人水多毛片| 五月开心婷婷网| 日本av免费视频播放| 丰满乱子伦码专区| 国产精品久久久久成人av| 18禁动态无遮挡网站| 久久亚洲国产成人精品v| 欧美xxⅹ黑人| 免费在线观看黄色视频的| av天堂久久9| 曰老女人黄片| 在线观看人妻少妇| 国产在线免费精品| 亚洲人成77777在线视频| 26uuu在线亚洲综合色| 99香蕉大伊视频| 欧美精品av麻豆av| 美女国产高潮福利片在线看| 国产精品麻豆人妻色哟哟久久| 女人久久www免费人成看片| av免费观看日本| 日日爽夜夜爽网站| 99久久人妻综合| 亚洲经典国产精华液单| 日本wwww免费看| 免费观看性生交大片5| 免费看av在线观看网站| 91精品三级在线观看| 久久99蜜桃精品久久| 成人18禁高潮啪啪吃奶动态图| 啦啦啦啦在线视频资源| 在线精品无人区一区二区三| 如日韩欧美国产精品一区二区三区| 黄色怎么调成土黄色| 纯流量卡能插随身wifi吗| 在线精品无人区一区二区三| 国产精品无大码| 黄片播放在线免费| 黄网站色视频无遮挡免费观看| 一区福利在线观看| 国产福利在线免费观看视频| 少妇的丰满在线观看| 午夜福利影视在线免费观看| 精品亚洲成a人片在线观看| 女性生殖器流出的白浆| 精品国产乱码久久久久久小说| 我要看黄色一级片免费的| 亚洲第一av免费看| av在线观看视频网站免费| 亚洲精品久久午夜乱码| 国产精品99久久99久久久不卡 | 香蕉精品网在线| 深夜精品福利| 2022亚洲国产成人精品| xxxhd国产人妻xxx| 水蜜桃什么品种好| 精品国产一区二区三区四区第35| 国产成人精品在线电影| 99久久中文字幕三级久久日本| 天美传媒精品一区二区| 久久人人97超碰香蕉20202| 亚洲成av片中文字幕在线观看 | 亚洲国产欧美在线一区| 2021少妇久久久久久久久久久| a级片在线免费高清观看视频| 中文天堂在线官网| 一级毛片电影观看| 中文字幕亚洲精品专区| 精品国产乱码久久久久久男人| 国产97色在线日韩免费| 天美传媒精品一区二区| 蜜桃国产av成人99| 看免费成人av毛片| 国产欧美日韩综合在线一区二区| 黄色怎么调成土黄色| 晚上一个人看的免费电影| 99re6热这里在线精品视频| 午夜av观看不卡| 黄色怎么调成土黄色| 晚上一个人看的免费电影| 岛国毛片在线播放| 又大又黄又爽视频免费| 久久 成人 亚洲| 久久久久久人妻| 亚洲成色77777| av片东京热男人的天堂| 亚洲国产欧美网| 精品人妻一区二区三区麻豆| 大陆偷拍与自拍| 桃花免费在线播放| 一级毛片我不卡| 亚洲少妇的诱惑av| 大陆偷拍与自拍| 亚洲一级一片aⅴ在线观看| 亚洲精品国产av成人精品| 国产一区有黄有色的免费视频| 久久精品国产亚洲av高清一级| 亚洲av福利一区| 日韩 亚洲 欧美在线| 久久婷婷青草| 1024视频免费在线观看| 9色porny在线观看| 飞空精品影院首页| 国产成人精品福利久久| 亚洲国产av新网站| 一区二区三区乱码不卡18| 午夜福利视频在线观看免费| 成年美女黄网站色视频大全免费| 少妇的丰满在线观看| 曰老女人黄片| 欧美激情 高清一区二区三区| 欧美日韩视频精品一区| 国产av国产精品国产| 老熟女久久久| 狠狠婷婷综合久久久久久88av| 国产成人精品婷婷| 97精品久久久久久久久久精品| 亚洲,欧美精品.| 观看av在线不卡| 少妇猛男粗大的猛烈进出视频| 久久精品国产鲁丝片午夜精品| av.在线天堂| 91久久精品国产一区二区三区| 久久精品久久久久久久性| 老司机影院毛片| 亚洲精品国产一区二区精华液| 亚洲欧美日韩另类电影网站| 五月开心婷婷网| 有码 亚洲区| 亚洲精品中文字幕在线视频| 女人久久www免费人成看片| 国产亚洲欧美精品永久| 久久精品夜色国产| 亚洲,欧美精品.| 欧美日韩精品成人综合77777| 高清不卡的av网站| 少妇精品久久久久久久| 久久99一区二区三区| av天堂久久9| 精品亚洲成国产av| av网站免费在线观看视频| 日本爱情动作片www.在线观看| 国产野战对白在线观看| 18在线观看网站| 久久韩国三级中文字幕| 久久精品熟女亚洲av麻豆精品| 只有这里有精品99| 亚洲 欧美一区二区三区| 国产在线一区二区三区精| 精品人妻偷拍中文字幕| 久久 成人 亚洲| 亚洲三级黄色毛片| 日本色播在线视频| 老司机影院成人| 久久鲁丝午夜福利片| 一本色道久久久久久精品综合| 亚洲欧美日韩另类电影网站| 色94色欧美一区二区| 欧美少妇被猛烈插入视频| 丝瓜视频免费看黄片| 国产精品不卡视频一区二区| 亚洲欧美成人综合另类久久久| 水蜜桃什么品种好| 老司机影院毛片| 青春草亚洲视频在线观看| 97在线人人人人妻| 母亲3免费完整高清在线观看 | 亚洲欧美清纯卡通| 深夜精品福利| 各种免费的搞黄视频| 麻豆精品久久久久久蜜桃| 99精国产麻豆久久婷婷| 黄色怎么调成土黄色| 一区福利在线观看| av片东京热男人的天堂| 免费在线观看完整版高清| 人人妻人人爽人人添夜夜欢视频| 亚洲欧美成人综合另类久久久| 亚洲精品第二区| 亚洲精品美女久久久久99蜜臀 | 亚洲欧美成人精品一区二区| 久久久精品94久久精品| 日本欧美视频一区| 欧美bdsm另类| 青春草国产在线视频| 国产福利在线免费观看视频| videossex国产| 欧美老熟妇乱子伦牲交| av片东京热男人的天堂| 久久精品久久久久久久性| 你懂的网址亚洲精品在线观看| 两个人看的免费小视频| 国产色婷婷99| 视频区图区小说| 秋霞在线观看毛片| 在线天堂最新版资源| videossex国产| 国产免费现黄频在线看| 天美传媒精品一区二区| 国产精品不卡视频一区二区| 国产激情久久老熟女| 日本欧美视频一区| 国产国语露脸激情在线看| 国产乱来视频区| 亚洲av日韩在线播放| 国产免费视频播放在线视频| 亚洲精华国产精华液的使用体验| av网站在线播放免费| 国产成人精品一,二区| 国产成人91sexporn| 在线亚洲精品国产二区图片欧美| 母亲3免费完整高清在线观看 | 亚洲伊人色综图| 五月伊人婷婷丁香| 婷婷成人精品国产| 制服诱惑二区| 伊人久久大香线蕉亚洲五| 五月伊人婷婷丁香| 超碰97精品在线观看| 老鸭窝网址在线观看| 中文乱码字字幕精品一区二区三区| 人人妻人人澡人人爽人人夜夜| 久久99热这里只频精品6学生| 日韩在线高清观看一区二区三区| 国产视频首页在线观看| 一区二区日韩欧美中文字幕| 亚洲av福利一区| 久久久国产欧美日韩av| 亚洲av成人精品一二三区| 午夜精品国产一区二区电影| 日本色播在线视频| 成人二区视频| 久久精品亚洲av国产电影网| 久久久久网色| 色94色欧美一区二区| 国产成人av激情在线播放| 80岁老熟妇乱子伦牲交| 国产精品人妻久久久影院| 成年av动漫网址| 亚洲美女黄色视频免费看| 另类精品久久| 国产黄色免费在线视频| 亚洲av免费高清在线观看| 久久国产精品大桥未久av| 久久久久久久精品精品| 久久毛片免费看一区二区三区| 亚洲精品国产一区二区精华液| 欧美精品高潮呻吟av久久| 五月伊人婷婷丁香| 如何舔出高潮| 哪个播放器可以免费观看大片| 婷婷色av中文字幕| 两性夫妻黄色片| 丝袜脚勾引网站| 国产av码专区亚洲av| 久久 成人 亚洲| 亚洲精品美女久久久久99蜜臀 | 国产精品麻豆人妻色哟哟久久| 日韩精品有码人妻一区| 激情视频va一区二区三区| 丝袜人妻中文字幕| 九色亚洲精品在线播放| 一级毛片黄色毛片免费观看视频| 久久国产亚洲av麻豆专区| 国产在线视频一区二区| 黑人欧美特级aaaaaa片| 亚洲国产精品999| 久久午夜综合久久蜜桃| 亚洲视频免费观看视频| 久久久精品94久久精品| 久久婷婷青草| 啦啦啦在线观看免费高清www| 高清在线视频一区二区三区| 欧美bdsm另类| 老汉色∧v一级毛片| 麻豆av在线久日| 国产片特级美女逼逼视频| 欧美日韩一区二区视频在线观看视频在线| 99久久综合免费| 亚洲中文av在线| 黄色毛片三级朝国网站| 欧美激情极品国产一区二区三区| 黄片播放在线免费| 亚洲精品一二三| 一级毛片电影观看| 久久这里有精品视频免费| 熟女少妇亚洲综合色aaa.| 国产又爽黄色视频| 精品午夜福利在线看| 丝袜喷水一区| 成人国产av品久久久| 欧美日韩综合久久久久久| 卡戴珊不雅视频在线播放| 在线天堂中文资源库| 久久精品久久精品一区二区三区| 一本一本久久a久久精品综合妖精 国产伦在线观看视频一区 | 9191精品国产免费久久| 亚洲人成77777在线视频| av国产久精品久网站免费入址| 性高湖久久久久久久久免费观看| 啦啦啦视频在线资源免费观看| 不卡视频在线观看欧美| 国产成人av激情在线播放| 亚洲男人天堂网一区| 久久人人爽av亚洲精品天堂| 日韩在线高清观看一区二区三区| 欧美国产精品一级二级三级| 热99久久久久精品小说推荐| 激情五月婷婷亚洲| 精品久久蜜臀av无| 日韩欧美一区视频在线观看| 伦理电影免费视频| 欧美+日韩+精品| 观看美女的网站| 免费大片黄手机在线观看| 9191精品国产免费久久| 国产av码专区亚洲av| 不卡av一区二区三区| 电影成人av| 亚洲三区欧美一区| 免费在线观看黄色视频的| 97在线人人人人妻| 国产成人免费观看mmmm| 国产成人精品一,二区| 欧美精品国产亚洲| 日韩电影二区| 青春草亚洲视频在线观看| 欧美日本中文国产一区发布| 嫩草影院入口| av.在线天堂| 国产 精品1| 国产极品天堂在线| 久久久久久免费高清国产稀缺| 街头女战士在线观看网站|