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

    新能源運行容量可信度的精細化評估方法

    2023-02-11 08:35:44王榆楗張嘯虎劉盾盾程浩忠
    可再生能源 2023年1期
    關鍵詞:系統(tǒng)可靠性風光出力

    王榆楗,吳 瓊,柳 璐,張嘯虎,桑 妲,劉盾盾,程浩忠

    (1.上海交通大學 電氣工程系,上海 200240;2.中國長江電力股份有限公司,北京 100038;3.國家電網有限公司華東分部,上海 200002;4.珠海澳大科技研究院,廣東 珠海 519085;5.貴州電網有限責任公司電力調度控制中心,貴州 貴陽 550000)

    0 引言

    以風電和光伏為代表的可再生能源裝機容量將持續(xù)快速增長,然而,風光可再生能源的間歇性與波動性導致其難以準確計入電力平衡[1]~[3]。當前電網以風電5%和光伏0%作為參與電力平衡的容量系數,無法滿足未來以可再生能源為主體的新型電力系統(tǒng)發(fā)展需求??陀^準確評估風光可再生能源可信度對于電力系統(tǒng)規(guī)劃、運行具有十分重要的意義。

    新能源機組能替代常規(guī)機組的容量即為置信容 量 或 可 信 容 量(Capacity Credit,CC)[4],約 有10余種指標用于定量評估[5],現有大量文獻分別從不同角度分別對風電、光伏容量可信度的計算方法及影響因素進行研究。文獻[6]研究風電穿透率水平、風資源優(yōu)劣程度對容量可信度的影響。文獻[7]基于人工神經網絡提出經驗模型,用于評估光伏容量可信度。文獻[8]在光伏容量可信度評估過程中考慮了天氣不確定性。文獻[9]評估了儲能對光伏容量可信度的影響。文獻[10]在可靠性計算中考慮需求響應的影響,并應用于可信度評估??紤]相關性方面,基于歷史數據的評估結果表明,風光新能源出力特性與系統(tǒng)負荷的相關性對于容量可信度評估結果具有顯著影響[11]。文獻[11]研究表明,風電出力與系統(tǒng)負荷之間的相關性越強,風電容量可信度就越大。文獻[12]通過Spearman秩相關系數描述光伏出力與負荷相關性,并用于光伏容量可信度評估。文獻[13]通過建立離散聯合概率分布對光伏出力與負荷相關性進行建模并進行可信度評估;考慮“光-荷”相關性的評估結果均表明光伏出力和負荷之間的相關性對于光伏容量可信度具有積極影響。文獻[14]通過單一Copula函數建立多風電場出力的聯合概率分布,提出了考慮多風電場出力相關性的容量可信度評估方法。文獻[15]在多光伏電站出力建模中考慮了空間相關性,并且分析了天氣類型對于光伏出力的影響。文獻[16]通過2維Archimedean Copula函數對風電、光伏出力相關性進行建模,并應用于可靠性評估。

    現 有 研 究 分 別 從“風-光”、“風-荷”、“光-荷”角度進行考慮。然而,同一時刻、不同位置的風速、光照強度和負荷分布特性各異,并且相鄰風電廠、光伏電站之間既不完全獨立,也不完全相關,具有非線性相關性,同時“風-光-荷”三者之間亦具有非線性相關性[17]。因此在可信度評估中對風電出力之間、光伏出力之間以及風光荷之間的相關性進行準確建模,其難點在于單一類型Copula函數對“風-光”或者“風-荷”相關性能夠較為準確的建模,但對于“風-光-荷”相關性建模,則準確性較低;此外采用參數分布也難以準確描述風光荷的分布特性。風光容量可信度評估均為長時間尺度(通常為1a),因此計算結果表征全年風、光新能源對該區(qū)域內電力系統(tǒng)充裕度的綜合貢獻[18],未能體現風、光新能源的時空分布特性。且僅有少量文獻同時考慮了風電、光伏聯合發(fā)電系統(tǒng)的綜合容量可信度評估[19],[20]。因此,現有可信度評估存在計算方法對于風-光-荷多維時空相關性的考慮有所欠缺,使得計算結果不夠精確、評價指標偏籠統(tǒng)、未計及系統(tǒng)運行備用,從而對于實際運行場景指導意義相對有限,有必要考慮更為精細化、具有實用性指導價值的擴展性指標與計算方法。

    本文采用D-vine Pair Copula函數對多維風光荷相關性進行建模,并在系統(tǒng)運行可靠性不變的前提下,定義各時刻運行容量可信度指標,對特定運行場景進行精細化考慮,采用我國某實際區(qū)域電網2018-2020年風光荷歷史數據,生成蘊含相關性的計算樣本,計算風光運行可信容量,以期為系統(tǒng)實際運行中的電力平衡分析提供參考。

    1 風光荷相關性數據建模及采樣

    1.1 Pair Copula風光荷相關性建模

    對 于 多 元 隨 機 變 量x=(x1,x2,…,xd),xk的 邊緣累積分布函數為Fk(xk),Sklar定理給出了多維聯合分布H(x)與Copula函數C(u)之間的關系:

    式 中:uk=Fk(xk),且uk服 從[0,1]均 勻 分 布。

    對式(1)求導,得到多維聯合概率密度函數:

    式 中:c(u)為Copula密 度 函 數;fk(xk)為xk的 邊 緣概率密度函數。

    由式(2)可見,通過估計隨機變量的邊緣密度函 數fk(xk)(k=1,2…,d),并 通 過 估 計d維Copula密度函數c(u)表示相關性結構,可求得聯合概率密 度 分 布f(x1,x2,…,xd)?,F 有Copula函 數 能 夠 對二維相關性建模較準確,但對于多維變量的相關性 建 模 則 不 夠 準 確[21],[22],而Pair Copula理 論[23]將多維Copula函數的密度函數c(u)拆分為一組2維Copula密度函數和隨機變量邊緣概率密度函數的乘積,從而能較好描地述復雜的多維相關性。

    本文采用D-vine Pair Copula結構作為c(u)的分解方式,d維D-vine結構如圖1所示。

    圖1 d維D-vine Pair Copula結構圖Fig.1 Structure diagram of d-dimensional D-vine Pair Copula

    d維D-vine Pair Copula結構有d-1層,第Li層有d-i個節(jié)點。

    結 構 中c(u)表 達 式 為

    通過Pair Copula將隨機變量兩兩合并,并采用相應的二元Copula函數描述,從而可以提高模型擬合精度,更加準確地描述多維“風-光-荷”之間不同的相關性結構。

    本文采用歐氏距離法[22],作為評價Copula模型對于原始樣本數據擬合效果的評價方法。

    記(x1i,x2i,…,xdi)為d維 隨 機 變 量x的 歷 史 數據 樣 本,其 中i=1,2,…,N,N為 歷 史 數 據 樣 本 總數。擬合得到的歐氏距離為

    式中:Cem(·)為從歷史數據中獲得的在樣本點處的實際概率值。

    采用非參數核密度估計法(Kernel Density Estimation,KDE)分別估計多維風光荷數據的概率 密 度 函 數fk(xk)、累 積 分 布 函 數Fk(xk)及 其 反 函數Fk-1(xk),再 利 用 累 積 分 布 函 數Fk(·)將 數 據xk變換為均勻分布的概率。

    隨機變量xk的核密度為

    式 中:h為 帶 寬;K(·)為 核 函 數;Xki為 隨 機 變 量xk的第i個歷史樣本。

    當樣本數據給定,fk(xk)便可由核函數以及帶寬確定。KDE中通常采用高斯核函數為

    對式(7)積分可得隨機變量xk的邊緣累積分布函數為

    1.2 建模及采樣流程

    基于D-vine Pair Copula結構以及非參數核密度估計方法,對d維風光荷相關性進行建模并采樣生成數據。步驟如下。

    ①讀 入d維 風 光 荷 歷 史 數 據(X1i,X2i,…,Xdi)i=1,2,…,N,利用KDE得到各隨機變量相應的邊緣 概 率 密 度 函 數fk(xk)和 累 積 分 布 函 數Fk(xk)。

    ②令uk=Fk(xk),將 歷 史 數 據 相 應 帶 入,得 到均 勻 分 布U=[U1,U2,…,Ud]的 樣 本 點。

    ③根據樣本數據U,對第L1層各Copula密度函數進行擬合,采用對數極大似然估計得到相應參數。依據式(6)對擬合度進行檢驗,選出第L1層各最優(yōu)2維Copula密度函數。

    ④將式(4)所示的條件分布作為新的隨機變量,并利用U的樣本和式(4)得到其樣本。采用式(4)中條件分布得到的樣本,依據式(6)對擬合度進行檢驗,選出第L2層各最優(yōu)Copula密度函數。

    ⑤通過前一層條件分布的樣本以及式(5)生成本層條件分布的樣本,與上述步驟類似,直至得到第Ld-1層的Copula密度函數,并得出式(3)。

    ⑥按數據進行采樣,生成蘊含d維風-光-荷相 關 性 的 數 據(X1,X2,…,Xd),采 樣 函 數 為

    式 中:W1,W2,…,Wd為 通 過 蒙 特 卡 洛 生 成 在[0,1]上服從均勻分布的隨機數;F-1為式(3)中各累積分布函數的逆函數。

    2 風光運行容量可信度評估方法

    2.1 風-光運行容量可信度指標

    現有容量可信度評估指標以等效可靠容量(Equivalent Firm Capacity,EFC)為 代 表,EFC為 實際系統(tǒng)除去風光機組并加入特定容量,且100%可靠的虛擬常規(guī)機組后,系統(tǒng)可靠性與實際系統(tǒng)可靠性相同,此時虛擬常規(guī)機組容量占風光新能源裝機容量的比例為風光容量可信度[18]。

    在實際運行中,為保證系統(tǒng)各時刻運行可靠性,需配置一定的備用容量,在風光新能源裝機較多的地區(qū),還需結合風光出力特性和參與電力平衡的比例,額外設置一定的負荷備用容量[24]。按照EFC定義,對于運行中風光容量可信度評估存在適應性問題。即目前運行中是以一定負荷比例設置備用,從而確定開機容量,此時存在原系統(tǒng)中的風光應該以多少容量計入系統(tǒng)開機容量的問題,開機容量改變,系統(tǒng)可靠性相應改變,從而直接影響到最終的計算結果。

    本文在系統(tǒng)運行可靠性不變的前提下,考慮系統(tǒng)運行備用,定義風光新能源各時刻的運行可信 容 量(Operational Credit Capacity,OCC),其 表達式為

    式中:Rt為系統(tǒng)t時刻的可靠性,選取失負荷概率(Loss of Load Probability,LOLP)作 為 可 靠 性 的 判斷標準,等式左側為不含風光系統(tǒng);C0為根據負荷設置備用容量后,系統(tǒng)常規(guī)機組開機容量;dt為t時刻系統(tǒng)負荷,等式右側為風光接入后系統(tǒng);C1為風光接入后系統(tǒng)常規(guī)機組開機容量;Cc-t為風光新能源在t時刻的運行可信容量,Cc-t與風光新能源裝機的比值,即為風光新能源t時刻的運行容量可信度。

    在OCC定義下,首先根據負荷設置備用,確定系統(tǒng)常規(guī)機組開機容量,并評估該時刻的系統(tǒng)可靠性;之后將風光新能源加入系統(tǒng),在等可靠性前提下,常規(guī)機組開機容量將相應減少,此時,常規(guī)機組開機容量的減少量 ΔC=C0-C1即為t時刻風光新能源運行可信容量Cc-t。

    2.2 考慮風-光-荷相關性的發(fā)電系統(tǒng)可靠性計算方法

    通過前文方法生成蘊含多維風-光-荷相關性的計算樣本后,采用基于重要抽樣法的蒙特卡洛法對發(fā)電系統(tǒng)可靠性進行計算。常規(guī)機組可靠性模型為

    式中:P(g)為系統(tǒng)狀態(tài)分布;Yi為系統(tǒng)狀態(tài)向量第i個元素;yi為常規(guī)機組i的狀態(tài);qi為常規(guī)機組i的強迫停運率。

    對于受端系統(tǒng)可靠性評估,本文將直流輸電系統(tǒng)的注入功率等效為多狀態(tài)電源,等效方式參見文獻[25]。此外,還需根據歷史運行方式考慮直流線路不同月份的輸送容量。例如對于水電送端直流,汛期(6-9月)輸送容量按額定容量計,其余月份則按60%額定容量作為雙極完全運行時的輸送容量。

    運行容量可信度評估過程中,涉及多次迭代計算系統(tǒng)可靠性,蒙特卡洛法需花費大量計算時間。為減小計算耗時,本文引入重要抽樣法,在維持樣本期望值恒定的前提下,改變現有樣本空間的概率分布,縮小計算指標的方差,以提高運算速率[26]。以 系 統(tǒng) 狀 態(tài) 函 數S(Y)的 期 望E(S)作 為 可靠性評估指標,采用方差系數作為收斂判據[26],其表達式為

    式中:Ω為狀態(tài)空間;S*(Y)為新狀態(tài)分布下的系統(tǒng)狀態(tài)函數。

    若新狀態(tài)分布P*(Y)能夠突出“重要事件”(即引起系統(tǒng)故障停運的事件),則能夠證明S*(Y)方 差 將 小 于S(Y)方 差。新 分 布P*(Y)又 稱為重要分布函數。定義系統(tǒng)重要分布函數為

    式中:k為重要乘子。

    2.3 風光運行容量可信度評估流程

    根據OCC定義,考慮相關性的風光新能源運行容量可信度評估流程,如圖2所示。

    圖2 考慮相關性的風光新能源運行容量可信度評估流程Fig.2 OCC evaluation flow chart of new energy considering relevance

    首先讀入原始系統(tǒng)數據:常規(guī)機組裝機容量、強迫停運率數據、直流受電功率、歷史數據集內每日i時刻的風電、光伏及負荷數據。根據負荷情況確定原始系統(tǒng)中常規(guī)機組開機容量C0,評估原始系統(tǒng)可靠性R0。開機容量為峰值負荷加上備用容量,系統(tǒng)備用需求源于電網運行中的不確定性,一方面是不可預計的電網事故(事故備用),另一方面是供需預測的不準確(負荷備用),通常為峰值負荷的2%~5%。然后,生成蘊含多維風-光-荷相關性的計算樣本,計算風光新能源接入后系統(tǒng)可靠性R1,通過弦截法[18]調節(jié)常規(guī)機組開機容量C1使風光新能源接入前后系統(tǒng)可靠性指標相差在一定誤差范圍內,此時,常規(guī)機組開機容量的減少量ΔC=C0-C1,即為t時刻風光新能源運行可信容量。

    3 算例分析

    我國某實際區(qū)域電網總裝機38317萬kW,電源結構如圖3所示。各子區(qū)域相對位置如圖4所示。饋入區(qū)域電網的跨區(qū)直流總額定容量為6976萬kW,網內最大負荷為31871萬kW;采用2018-2020年迎峰度夏期間(6-8月)15min間隔的風光荷歷史數據擬合邊緣分布,其中風電、光伏出力數據為出力占裝機的百分比,由5個子區(qū)域分別采集,各子區(qū)域風光裝機占比如表1所示,負荷則按最大功率標幺化。

    圖3 區(qū)域電網電源結構Fig.3 Power source structure of regional power grid

    圖4 各子區(qū)域相對位置Fig.4 The relative position of each sub-region

    表1 子區(qū)域風光新能源裝機結構Table1 Sub-regional new energy installed structure

    3.1 風-光-荷相關性模型驗證

    對歷史數據集內每日12:00-13:00的風電、光伏及負荷數據,采用KDE擬合邊緣概率密度以及累積概率,擬合結果如圖5所示。

    圖5 歷史數據集內12:00-13:00風電、光伏出力和負荷的概率密度和累積概率Fig.5 Probability density and cumulative probability of wind power,photovoltaic power and load

    由圖5可以看出,僅采用Weibull分布無法準確描述各風電出力的概率密度和累積概率。同樣,各光伏出力也難以用Beta分布進行描述。而負荷的概率密度曲線在形態(tài)上也與正態(tài)分布明顯不同。與風-光-荷服從特定參數分布相比,采用KDE有利于提高Pair Copula模型準確度。采用Wi,Pi(i=1,…,5)及L分 別 代 表 各 風 電、光 伏 及 負荷功率,對多維風光荷之間分別進行2維最優(yōu)Copula建模,最優(yōu)Copula函數如表2所示。分別考察高維normal Copula、高維t Copula,以及Dvine Pair Copula對于多維風光荷擬合效果,歐氏距離如表3所示。

    表3 不同Copula函數對于高維風光荷建模的歐式距離Table3 The Euclidean distance of different Copula functions for modeling high-dimensional wind power,photovoltaic and load

    由表2,3可以看出,大量2維變量的最優(yōu)Copula并非normal Copula或t Copula,相較于高維normal Copula與 高 維t Coupula,D-vine Pair Copula由于引入多種Copula函數,能夠更好的對多維風光荷間復雜的整體與局部相關性進行建模,從而為風光運行容量可信度計算的準確性提供保障。

    表2 各2維最優(yōu)Copula函數類型Table2 Optimal Copula function types of2dimension

    3.2 迎峰度夏期間峰值負荷日風-光運行可信容量評估

    以迎峰度夏期間出現峰值負荷(31871萬kW)時的日負荷情況為例,確定系統(tǒng)開機容量,并評估系統(tǒng)運行可靠性,在系統(tǒng)運行可靠性不變的前提下,基于本文所述方法進行該日各時刻(24 h)風、光運行可信容量評估,結果如圖6所示。由于夏季光伏僅在6:00-19:00出力,因此其余時刻僅考慮多維風-荷相關性。

    圖6 迎峰度夏期間峰值負荷日風光運行可信容量計算結果Fig.6 Calculation results of OCC of peak load during summer

    根據計算結果,日間風光運行容量可信度均大于10%,最高達到26.31%;夜間僅考慮風電,運行容量可信度最低為8.11%,其余時刻接近10%。由于風光新能源的間歇性、波動性和隨機性,目前我國調度系統(tǒng)往往在日前計劃中對風光出力考慮的比較保守,因此在未來安排系統(tǒng)日前發(fā)電計劃時,可根據負荷預測結果按本文所述方法計算風光運行可信容量,并結合風電、光伏出力預測,從而更加合理的安排常規(guī)機組的發(fā)電計劃。

    EFC定義下,僅采用6-8月風光荷歷史數據進行評估,容量可信度為19.43%,評估結果表征風光新能源對系統(tǒng)電力充裕度的綜合貢獻,忽略了風光新能源時空分布特性,對于實際運行中的電力平衡指導意義相對有限。針對不同季節(jié)的風光新能源時空分布特性,同樣可采用本文所提指標與方法進行相應計算。相較于僅以一個定值表示風光可信容量的傳統(tǒng)計算方法,本文對可信度評估指標體系及評估方法進行了拓展,且更加精細化的考慮了多維風光荷時空分布特性,以期為風光新能源參與電力平衡提供參考。

    3.3 風、光裝機容量對運行可信容量的影響

    “雙碳”背景下,風電、光伏裝機容量將進一步快速增加。為探究風光裝機容量對于運行可信容量的影響,本節(jié)在保持現有負荷不變的情況下依據現有各子區(qū)域風光裝機容量的占比關系對風光總裝機容量進行線性外推,計算午間12點不同風光總裝機容量下的風光運行可信容量,結果如圖7所示。

    圖7 不同風光裝機容量下運行可信度評估結果Fig.7 OCC results under different installed capacity

    由圖7可以看出,隨著風光裝機容量增加,可信度從1000萬kW裝機時的29.74%逐漸下降至12000萬kW裝機時的19.97%。以1000萬kW為梯度,增加同樣風光裝機的邊際效益下降,即風光運行可信容量增長幅度降低,存在飽和效應。這是由于風光裝機容量較少時,系統(tǒng)可靠性主要取決于常規(guī)機組狀態(tài),而風光裝機容量增加后,根據指標定義,常規(guī)機組開機容量減少,因此風光出力的波動性對于系統(tǒng)可靠性的影響會進一步增大,為維持系統(tǒng)等可靠性,常規(guī)機組開機容量的減少幅度也會相應減小。

    4 結論

    本文針對可信度指標尚未深入應用于運行階段的問題,定義了運行容量可信度指標,提出考慮多維風光荷相關性的運行容量可信度評估方法,該方法引入D-vine Pair Copula,采用歷年風-光-荷數據對多維風光荷相關性進行建模,利用基于重要抽樣法的蒙特卡洛法對系統(tǒng)可靠性進行評估,并采用弦截法迭代,實現對給定運行場景的風光運行可信容量評估。本文通過算例分析得到結論如下:①多維風光荷具有復雜非線性相關性,Pair Copula函數對于各隨機變量之間的相關性均采用最優(yōu)2維Copula函數建模,使得生成的風光荷數據包含更多的相關性特征,此外采用KDE能夠有效提高建模準確性。因此相較于傳統(tǒng)容量可信度計算中對于風-光-荷的考慮方法,本文采用方法考慮因素更加全面,能獲得更準確的計算結果;②本文所提指標與方法更精細化考慮了多維風光荷相關性,以某實際區(qū)域電網迎峰度夏期間峰值負荷日的負荷情況為例,計算各時刻(24 h)風光運行容量可信度,結果為8.11%~26.31%。運行中僅以風電5%和光伏0%作為參與電力平衡的容量系數會導致過多的開機,進而影響風光消納。安排系統(tǒng)日前發(fā)電計劃時,可根據計算結果并結合風光出力預測,為調度安排合理開機計劃提供參考。探究了風光裝機容量對于運行可信容量的影響,計算結果表明,隨著風光裝機容量增加,可信容量存在飽和效應。

    儲能可以彌補新能源在隨機波動性方面的先天缺陷,隨著儲能系統(tǒng)的大規(guī)模引入,風-光-儲聯合系統(tǒng)有望得到廣泛應用,在后續(xù)研究中可考慮采用序貫蒙特卡洛模擬法進行風-光-儲聯合系統(tǒng)的可靠性評估,以及考慮風光荷時間相關性的建模與采樣,進而精細化評估風-光-儲聯合系統(tǒng)的運行容量可信度。

    猜你喜歡
    系統(tǒng)可靠性風光出力
    風光新580
    汽車觀察(2021年11期)2021-04-24 20:47:38
    風光如畫
    海峽姐妹(2019年12期)2020-01-14 03:25:02
    試析提高配網系統(tǒng)可靠性的技術措施
    電子制作(2019年20期)2019-12-04 03:51:54
    電氣化鐵路牽引系統(tǒng)可靠性分析
    風光ix5:當轎跑邂逅SUV
    汽車觀察(2018年12期)2018-12-26 01:05:36
    各國首都風光
    基于故障樹模型的光伏跟蹤系統(tǒng)可靠性分析
    電測與儀表(2016年3期)2016-04-12 00:27:30
    風電場有功出力的EEMD特性分析
    要爭做出力出彩的黨員干部
    河南電力(2016年5期)2016-02-06 02:11:35
    風電場群出力的匯聚效應分析
    電測與儀表(2015年5期)2015-04-09 11:31:12
    欧美性猛交黑人性爽| 久久久久久大精品| 国产高清视频在线观看网站| 国产亚洲91精品色在线| 亚洲真实伦在线观看| 久久精品久久久久久噜噜老黄 | 99热6这里只有精品| av女优亚洲男人天堂| 日本一二三区视频观看| 国产精品久久久久久精品电影小说 | 精品一区二区三区视频在线| 少妇熟女欧美另类| 午夜激情欧美在线| 久久精品国产亚洲av天美| 久久久久性生活片| 日韩亚洲欧美综合| 国产黄色视频一区二区在线观看 | 午夜福利在线观看免费完整高清在| 黄片无遮挡物在线观看| 国产一区亚洲一区在线观看| 一级二级三级毛片免费看| 偷拍熟女少妇极品色| 国产亚洲精品久久久com| 啦啦啦啦在线视频资源| 国产一区有黄有色的免费视频 | 午夜福利视频1000在线观看| 性插视频无遮挡在线免费观看| 六月丁香七月| 午夜久久久久精精品| 国产白丝娇喘喷水9色精品| 成人三级黄色视频| 色播亚洲综合网| 亚洲av不卡在线观看| 麻豆成人av视频| 国产乱人视频| 国产精品精品国产色婷婷| 国产精品美女特级片免费视频播放器| 国产成人午夜福利电影在线观看| 国产精品福利在线免费观看| 一夜夜www| 韩国高清视频一区二区三区| 久久久久久久久久成人| 在线观看一区二区三区| 精品久久久久久久末码| 丝袜喷水一区| 91久久精品电影网| 欧美性感艳星| 中文字幕亚洲精品专区| 欧美激情久久久久久爽电影| 日本wwww免费看| 嫩草影院入口| 久久久精品大字幕| 边亲边吃奶的免费视频| 搡女人真爽免费视频火全软件| 一级毛片久久久久久久久女| 噜噜噜噜噜久久久久久91| 我的老师免费观看完整版| 全区人妻精品视频| 亚洲精品aⅴ在线观看| 国产成人freesex在线| 蜜桃久久精品国产亚洲av| 蜜桃久久精品国产亚洲av| 精品国内亚洲2022精品成人| kizo精华| 黄色配什么色好看| 成年版毛片免费区| 日韩三级伦理在线观看| 熟女人妻精品中文字幕| 高清视频免费观看一区二区 | 午夜a级毛片| 嫩草影院新地址| 如何舔出高潮| 国产熟女欧美一区二区| 97在线视频观看| 少妇人妻精品综合一区二区| 国产精品人妻久久久久久| 小说图片视频综合网站| av免费观看日本| 97在线视频观看| 日韩欧美精品免费久久| 蜜桃亚洲精品一区二区三区| 精品久久久噜噜| 天堂av国产一区二区熟女人妻| 欧美高清成人免费视频www| 岛国毛片在线播放| 身体一侧抽搐| 亚洲美女视频黄频| 99在线视频只有这里精品首页| 夜夜爽夜夜爽视频| 久久精品人妻少妇| 欧美zozozo另类| 日本三级黄在线观看| .国产精品久久| 亚洲av男天堂| 国产在视频线精品| 国产成人午夜福利电影在线观看| 久久精品久久久久久噜噜老黄 | 欧美精品一区二区大全| 亚洲伊人久久精品综合 | 午夜精品国产一区二区电影 | ponron亚洲| 欧美日韩精品成人综合77777| 久久久久久久午夜电影| 国产精品久久久久久av不卡| 久久国内精品自在自线图片| 人妻系列 视频| 成人亚洲精品av一区二区| 成人欧美大片| 色网站视频免费| 亚州av有码| 成人性生交大片免费视频hd| 日韩欧美在线乱码| 久久久久精品久久久久真实原创| 日本欧美国产在线视频| 成人午夜精彩视频在线观看| 亚洲精品成人久久久久久| 亚洲国产精品成人综合色| 中文字幕人妻熟人妻熟丝袜美| 亚洲第一区二区三区不卡| 又粗又硬又长又爽又黄的视频| 国产精品一区www在线观看| 国产又色又爽无遮挡免| 26uuu在线亚洲综合色| 一边摸一边抽搐一进一小说| 国产在线男女| 深夜a级毛片| av卡一久久| 日本一本二区三区精品| av在线蜜桃| 丰满乱子伦码专区| 国产精品久久电影中文字幕| 成人毛片60女人毛片免费| 亚洲精品国产成人久久av| 欧美3d第一页| 亚洲欧美精品自产自拍| 我要搜黄色片| 一级爰片在线观看| 精品人妻视频免费看| or卡值多少钱| 中文字幕制服av| 精品国产一区二区三区久久久樱花 | 成人漫画全彩无遮挡| 亚洲精品456在线播放app| 国产精品.久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 又黄又爽又刺激的免费视频.| 纵有疾风起免费观看全集完整版 | 97人妻精品一区二区三区麻豆| 日韩强制内射视频| 简卡轻食公司| 青春草视频在线免费观看| 国产免费福利视频在线观看| 日韩大片免费观看网站 | 国产精品久久久久久久久免| 看片在线看免费视频| 亚洲成色77777| 男女下面进入的视频免费午夜| 日本免费a在线| 岛国在线免费视频观看| 大香蕉97超碰在线| 高清在线视频一区二区三区 | 国产 一区精品| 欧美一区二区国产精品久久精品| 乱系列少妇在线播放| videos熟女内射| 青春草视频在线免费观看| 毛片女人毛片| 三级毛片av免费| 联通29元200g的流量卡| 精品国产露脸久久av麻豆 | 观看免费一级毛片| 日韩,欧美,国产一区二区三区 | 成人鲁丝片一二三区免费| 成人综合一区亚洲| 免费av不卡在线播放| www日本黄色视频网| 免费看美女性在线毛片视频| 日韩 亚洲 欧美在线| 国产av一区在线观看免费| 久久精品国产亚洲av涩爱| 美女大奶头视频| 中文字幕熟女人妻在线| 亚洲国产色片| 中文字幕人妻熟人妻熟丝袜美| 夫妻性生交免费视频一级片| 国产私拍福利视频在线观看| 亚洲va在线va天堂va国产| 免费看日本二区| 亚洲精华国产精华液的使用体验| 老司机影院毛片| 日韩精品有码人妻一区| or卡值多少钱| 99久国产av精品国产电影| 极品教师在线视频| 伦精品一区二区三区| av在线观看视频网站免费| 男女下面进入的视频免费午夜| 亚洲av电影在线观看一区二区三区 | 国产91av在线免费观看| 小蜜桃在线观看免费完整版高清| 大话2 男鬼变身卡| 国产亚洲91精品色在线| 一区二区三区乱码不卡18| 免费看av在线观看网站| 欧美潮喷喷水| 成年女人看的毛片在线观看| 色综合亚洲欧美另类图片| 久久韩国三级中文字幕| 欧美日本视频| 国产伦精品一区二区三区视频9| 午夜福利在线观看吧| 成年免费大片在线观看| 夜夜爽夜夜爽视频| 国产av码专区亚洲av| 久久精品国产自在天天线| 最近视频中文字幕2019在线8| 亚洲最大成人手机在线| 国产av一区在线观看免费| 欧美最新免费一区二区三区| 天堂√8在线中文| 麻豆av噜噜一区二区三区| 白带黄色成豆腐渣| 亚洲国产欧美人成| 国产精品蜜桃在线观看| 精品熟女少妇av免费看| 国产黄片视频在线免费观看| 日韩视频在线欧美| av免费观看日本| 国内精品宾馆在线| www.av在线官网国产| 99久久精品热视频| 在线免费观看不下载黄p国产| 成人高潮视频无遮挡免费网站| 黑人高潮一二区| 桃色一区二区三区在线观看| 中文资源天堂在线| 欧美又色又爽又黄视频| 成年女人永久免费观看视频| 插阴视频在线观看视频| 国产精品蜜桃在线观看| 亚洲精华国产精华液的使用体验| 永久免费av网站大全| 亚洲精品乱久久久久久| 久久精品久久精品一区二区三区| 大又大粗又爽又黄少妇毛片口| www日本黄色视频网| 99在线视频只有这里精品首页| 欧美+日韩+精品| 亚洲精品色激情综合| 国产单亲对白刺激| 男女啪啪激烈高潮av片| 91精品伊人久久大香线蕉| 人人妻人人看人人澡| 久久婷婷人人爽人人干人人爱| 乱系列少妇在线播放| 在线观看66精品国产| 狂野欧美激情性xxxx在线观看| 精品久久久久久成人av| 日本免费在线观看一区| 午夜亚洲福利在线播放| 亚洲欧美成人精品一区二区| 亚洲自偷自拍三级| 免费看日本二区| 边亲边吃奶的免费视频| 免费观看人在逋| 高清毛片免费看| 天天躁夜夜躁狠狠久久av| 性插视频无遮挡在线免费观看| 久久久久久久久中文| 亚洲av.av天堂| 久久久色成人| 欧美日韩综合久久久久久| 亚洲欧洲日产国产| 国产成人精品一,二区| 亚洲成人精品中文字幕电影| 熟女人妻精品中文字幕| 亚洲中文字幕一区二区三区有码在线看| 一夜夜www| 大又大粗又爽又黄少妇毛片口| 日本wwww免费看| 韩国高清视频一区二区三区| 久久久午夜欧美精品| 一级毛片我不卡| 国产成人一区二区在线| 亚洲国产欧美在线一区| 午夜久久久久精精品| 哪个播放器可以免费观看大片| 午夜精品在线福利| 国产精品一区二区在线观看99 | 国国产精品蜜臀av免费| 亚洲人成网站高清观看| 久久久久久伊人网av| 国产精品蜜桃在线观看| 亚洲乱码一区二区免费版| 亚洲欧美日韩高清专用| .国产精品久久| 午夜亚洲福利在线播放| 国产av一区在线观看免费| 午夜福利在线观看免费完整高清在| 国产免费视频播放在线视频 | 看黄色毛片网站| 黄片wwwwww| 一级黄片播放器| 夜夜看夜夜爽夜夜摸| 嫩草影院新地址| 中文字幕av在线有码专区| 国产欧美另类精品又又久久亚洲欧美| 亚洲天堂国产精品一区在线| 99热网站在线观看| 色尼玛亚洲综合影院| 国产一区二区三区av在线| 熟女人妻精品中文字幕| 久久热精品热| 最近手机中文字幕大全| 在线天堂最新版资源| 99久久无色码亚洲精品果冻| 欧美不卡视频在线免费观看| 欧美zozozo另类| 麻豆成人av视频| 九九热线精品视视频播放| 亚洲天堂国产精品一区在线| 日日撸夜夜添| 最近最新中文字幕大全电影3| 久久精品国产自在天天线| 色视频www国产| 亚洲精品乱码久久久久久按摩| 国产精品久久久久久精品电影小说 | 国产一区二区在线av高清观看| videos熟女内射| 热99在线观看视频| 亚洲国产精品久久男人天堂| www日本黄色视频网| 国产69精品久久久久777片| 日韩人妻高清精品专区| 国产精品美女特级片免费视频播放器| 精品一区二区三区人妻视频| 日韩亚洲欧美综合| 成人毛片60女人毛片免费| 一区二区三区乱码不卡18| 一个人看的www免费观看视频| 免费观看性生交大片5| 亚洲,欧美,日韩| 国产视频首页在线观看| 日韩欧美在线乱码| 乱人视频在线观看| 亚洲成av人片在线播放无| 成人特级av手机在线观看| 大又大粗又爽又黄少妇毛片口| 国产精品,欧美在线| 亚洲精品一区蜜桃| www.色视频.com| 水蜜桃什么品种好| 午夜福利网站1000一区二区三区| 国产精品一区www在线观看| 国产精品一区二区性色av| 国产亚洲av嫩草精品影院| 精品国产一区二区三区久久久樱花 | 三级国产精品片| 精品人妻偷拍中文字幕| 男人和女人高潮做爰伦理| 永久网站在线| 99在线视频只有这里精品首页| 国产一级毛片七仙女欲春2| 亚洲精品国产av成人精品| 成年版毛片免费区| 欧美成人a在线观看| 国产精品1区2区在线观看.| 成人午夜高清在线视频| 精品酒店卫生间| 日日摸夜夜添夜夜爱| 天天一区二区日本电影三级| 一个人免费在线观看电影| 久久久久免费精品人妻一区二区| 免费观看人在逋| 特级一级黄色大片| 女人十人毛片免费观看3o分钟| 国产免费男女视频| 99国产精品一区二区蜜桃av| 欧美最新免费一区二区三区| 日韩高清综合在线| 中文亚洲av片在线观看爽| 国语对白做爰xxxⅹ性视频网站| 色网站视频免费| 97热精品久久久久久| 黄色一级大片看看| 一个人观看的视频www高清免费观看| 亚洲色图av天堂| 毛片一级片免费看久久久久| 九九在线视频观看精品| 99久久无色码亚洲精品果冻| 日日摸夜夜添夜夜添av毛片| 欧美精品国产亚洲| 亚洲美女搞黄在线观看| 国产熟女欧美一区二区| 国产爱豆传媒在线观看| 久久久精品欧美日韩精品| 啦啦啦韩国在线观看视频| 久久婷婷人人爽人人干人人爱| 免费看av在线观看网站| 亚洲精品亚洲一区二区| 午夜福利视频1000在线观看| 三级经典国产精品| 美女xxoo啪啪120秒动态图| 青春草视频在线免费观看| 久久婷婷人人爽人人干人人爱| 免费看av在线观看网站| 国产爱豆传媒在线观看| 在现免费观看毛片| 网址你懂的国产日韩在线| 国产精品女同一区二区软件| 99热精品在线国产| 亚洲国产精品专区欧美| 国产亚洲一区二区精品| 男女啪啪激烈高潮av片| 午夜精品在线福利| 老司机影院毛片| 91在线精品国自产拍蜜月| 精品人妻一区二区三区麻豆| 国产伦理片在线播放av一区| 亚洲精品久久久久久婷婷小说 | 国产精品女同一区二区软件| 看片在线看免费视频| 久久人人爽人人爽人人片va| 精品99又大又爽又粗少妇毛片| av播播在线观看一区| 日韩一区二区三区影片| 亚洲av日韩在线播放| 两个人视频免费观看高清| 日韩一区二区视频免费看| 亚洲欧美成人综合另类久久久 | 国产爱豆传媒在线观看| 欧美xxxx黑人xx丫x性爽| 成人av在线播放网站| 特级一级黄色大片| 国产男人的电影天堂91| 国产色爽女视频免费观看| 成人国产麻豆网| 一卡2卡三卡四卡精品乱码亚洲| 国产精品综合久久久久久久免费| 男插女下体视频免费在线播放| 欧美激情久久久久久爽电影| 久久久精品欧美日韩精品| 欧美97在线视频| 美女黄网站色视频| 老司机福利观看| 日韩强制内射视频| 久久久久性生活片| 午夜福利在线观看免费完整高清在| 国产精品女同一区二区软件| 久久久a久久爽久久v久久| 久久99热这里只有精品18| 国产精品久久久久久久久免| 久久热精品热| 亚洲四区av| 国产免费男女视频| 国产亚洲午夜精品一区二区久久 | 少妇的逼好多水| 国产探花在线观看一区二区| 美女黄网站色视频| 精品久久久久久成人av| 全区人妻精品视频| 性插视频无遮挡在线免费观看| 高清视频免费观看一区二区 | 久久久久久国产a免费观看| 久久久久久久久久久丰满| 亚洲精品,欧美精品| 国产精品国产三级国产专区5o | 人妻系列 视频| 亚洲精品日韩在线中文字幕| 少妇熟女欧美另类| 男的添女的下面高潮视频| 中文乱码字字幕精品一区二区三区 | 欧美97在线视频| 爱豆传媒免费全集在线观看| 久久久久久久久久久丰满| 搡女人真爽免费视频火全软件| 内地一区二区视频在线| 黄片wwwwww| 精品久久久久久成人av| 婷婷色综合大香蕉| 亚洲欧美精品综合久久99| 精品久久久久久久久亚洲| 超碰97精品在线观看| 婷婷色综合大香蕉| 国产91av在线免费观看| 草草在线视频免费看| 高清毛片免费看| 久久精品国产99精品国产亚洲性色| 国产高清国产精品国产三级 | 亚洲丝袜综合中文字幕| 国产高清视频在线观看网站| 老女人水多毛片| 久久精品久久久久久噜噜老黄 | 久久久久久九九精品二区国产| 国产精品久久视频播放| 伊人久久精品亚洲午夜| 中文字幕熟女人妻在线| 可以在线观看毛片的网站| 国产精品国产三级国产av玫瑰| 成人二区视频| 国产亚洲5aaaaa淫片| 久久久久性生活片| 亚洲美女视频黄频| 亚洲丝袜综合中文字幕| 日日干狠狠操夜夜爽| 九九在线视频观看精品| 久久人人爽人人片av| 亚洲美女视频黄频| 成人一区二区视频在线观看| 校园人妻丝袜中文字幕| 高清午夜精品一区二区三区| 麻豆精品久久久久久蜜桃| 久久这里只有精品中国| 国产淫语在线视频| 亚洲欧美日韩卡通动漫| 1000部很黄的大片| 如何舔出高潮| 亚洲熟妇中文字幕五十中出| 高清视频免费观看一区二区 | 国产av一区在线观看免费| 男女边吃奶边做爰视频| 少妇熟女aⅴ在线视频| 亚洲四区av| 村上凉子中文字幕在线| 久久人妻av系列| 最近视频中文字幕2019在线8| av天堂中文字幕网| 午夜激情欧美在线| 亚洲怡红院男人天堂| 国产伦精品一区二区三区四那| 99国产精品一区二区蜜桃av| 欧美精品一区二区大全| 中文字幕亚洲精品专区| 久久精品影院6| 少妇熟女aⅴ在线视频| 好男人在线观看高清免费视频| 亚洲成人精品中文字幕电影| 又黄又爽又刺激的免费视频.| 免费大片18禁| 日本免费在线观看一区| 婷婷色av中文字幕| 国产欧美另类精品又又久久亚洲欧美| 欧美+日韩+精品| 亚洲四区av| 久久国产乱子免费精品| 啦啦啦观看免费观看视频高清| 一夜夜www| 我的老师免费观看完整版| 国产v大片淫在线免费观看| 国产高清视频在线观看网站| 天天一区二区日本电影三级| 亚洲国产色片| 国产一区二区在线观看日韩| 天堂网av新在线| 免费观看性生交大片5| 纵有疾风起免费观看全集完整版 | 级片在线观看| 噜噜噜噜噜久久久久久91| 欧美日韩综合久久久久久| 精品无人区乱码1区二区| 婷婷六月久久综合丁香| 国产在线男女| 丰满少妇做爰视频| 日本与韩国留学比较| 国产又黄又爽又无遮挡在线| 免费看美女性在线毛片视频| 亚洲中文字幕一区二区三区有码在线看| 成人鲁丝片一二三区免费| 在线免费观看的www视频| 日日摸夜夜添夜夜添av毛片| 禁无遮挡网站| 美女大奶头视频| 欧美97在线视频| 久久精品久久精品一区二区三区| 中文天堂在线官网| 欧美成人午夜免费资源| 成人特级av手机在线观看| 亚洲经典国产精华液单| 一级黄片播放器| 99久国产av精品| 国产一区二区在线av高清观看| h日本视频在线播放| 午夜激情福利司机影院| 色吧在线观看| 精品少妇黑人巨大在线播放 | 欧美日本亚洲视频在线播放| 中文在线观看免费www的网站| 欧美性猛交╳xxx乱大交人| 国产午夜福利久久久久久| 如何舔出高潮| 精品久久国产蜜桃| 日韩国内少妇激情av| 成人性生交大片免费视频hd| 日本wwww免费看| 欧美日韩国产亚洲二区| 老司机影院成人| 婷婷色综合大香蕉| 麻豆一二三区av精品| 97热精品久久久久久| 人人妻人人澡欧美一区二区| 最近最新中文字幕大全电影3| 一夜夜www| 亚洲激情五月婷婷啪啪| 久久精品国产亚洲av天美| 免费电影在线观看免费观看| 亚洲av中文字字幕乱码综合| 搡老妇女老女人老熟妇| 欧美日韩综合久久久久久| 午夜免费激情av| 寂寞人妻少妇视频99o| 欧美xxxx性猛交bbbb| 日韩欧美精品v在线| 18禁动态无遮挡网站| 久久99热这里只有精品18| 精品久久久久久久久久久久久|