• <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
    精品久久蜜臀av无| 亚洲精品国产精品久久久不卡| 亚洲国产欧美一区二区综合| 91麻豆av在线| 9191精品国产免费久久| 韩国精品一区二区三区| 国产有黄有色有爽视频| 欧美日韩一级在线毛片| 91大片在线观看| av电影中文网址| 精品国产超薄肉色丝袜足j| 国产一区有黄有色的免费视频| 视频区欧美日本亚洲| 亚洲精品自拍成人| 欧美国产精品一级二级三级| 777久久人妻少妇嫩草av网站| 欧美av亚洲av综合av国产av| 免费观看人在逋| 久久这里只有精品19| 亚洲精品久久午夜乱码| 欧美日韩亚洲高清精品| 国产精品.久久久| 成人黄色视频免费在线看| 丰满迷人的少妇在线观看| 热99国产精品久久久久久7| 人人妻,人人澡人人爽秒播| 国产欧美亚洲国产| 麻豆国产av国片精品| 9色porny在线观看| 丰满饥渴人妻一区二区三| 国产精品久久久人人做人人爽| 天堂中文最新版在线下载| 日本vs欧美在线观看视频| 国产免费视频播放在线视频| 国产深夜福利视频在线观看| 叶爱在线成人免费视频播放| 久久精品国产99精品国产亚洲性色 | 一本综合久久免费| 在线观看一区二区三区激情| 交换朋友夫妻互换小说| 欧美国产精品一级二级三级| 建设人人有责人人尽责人人享有的| 两个人免费观看高清视频| 99国产精品一区二区蜜桃av | 我要看黄色一级片免费的| 国产男靠女视频免费网站| 女人久久www免费人成看片| 男人操女人黄网站| 亚洲欧美日韩高清在线视频 | 亚洲人成伊人成综合网2020| 香蕉久久夜色| 精品国产一区二区三区四区第35| 夜夜爽天天搞| 亚洲国产中文字幕在线视频| 久久国产精品影院| 久久av网站| 免费人妻精品一区二区三区视频| 搡老乐熟女国产| 麻豆成人av在线观看| 久热这里只有精品99| 天天操日日干夜夜撸| 亚洲国产欧美在线一区| 一级a爱视频在线免费观看| 日韩精品免费视频一区二区三区| 这个男人来自地球电影免费观看| 欧美大码av| 欧美成人午夜精品| 99在线人妻在线中文字幕 | 高潮久久久久久久久久久不卡| 亚洲熟女毛片儿| 亚洲成人国产一区在线观看| 中文欧美无线码| 国产一区二区 视频在线| 亚洲熟女毛片儿| 午夜福利,免费看| 一区二区三区精品91| 91大片在线观看| 一区二区三区激情视频| 国产精品九九99| 99久久国产精品久久久| 99国产精品一区二区三区| 在线观看一区二区三区激情| 色精品久久人妻99蜜桃| 一级毛片女人18水好多| 人人妻,人人澡人人爽秒播| 免费高清在线观看日韩| av网站在线播放免费| 国产精品久久久久久精品电影小说| 日韩 欧美 亚洲 中文字幕| 亚洲黑人精品在线| 9色porny在线观看| 日本五十路高清| 久久人妻福利社区极品人妻图片| www.精华液| 在线观看66精品国产| 亚洲精品一二三| 又大又爽又粗| 日日摸夜夜添夜夜添小说| 亚洲 国产 在线| 免费高清在线观看日韩| 亚洲国产中文字幕在线视频| 999精品在线视频| 日韩三级视频一区二区三区| 一区二区av电影网| 成人永久免费在线观看视频 | 搡老乐熟女国产| 另类精品久久| 免费高清在线观看日韩| √禁漫天堂资源中文www| 女同久久另类99精品国产91| 亚洲精品美女久久久久99蜜臀| 麻豆国产av国片精品| 99国产精品99久久久久| 亚洲国产看品久久| 亚洲精品自拍成人| 热re99久久国产66热| 99国产极品粉嫩在线观看| 香蕉丝袜av| 亚洲avbb在线观看| 亚洲五月色婷婷综合| 亚洲av美国av| 亚洲成a人片在线一区二区| 国产精品熟女久久久久浪| 亚洲精品久久成人aⅴ小说| 一级,二级,三级黄色视频| 成人影院久久| 精品亚洲乱码少妇综合久久| 老鸭窝网址在线观看| 搡老乐熟女国产| 国产成人欧美| 老熟女久久久| 狠狠狠狠99中文字幕| 欧美日韩亚洲国产一区二区在线观看 | 在线观看免费高清a一片| 国产在线一区二区三区精| 肉色欧美久久久久久久蜜桃| 亚洲avbb在线观看| 久久久久久久大尺度免费视频| 视频区欧美日本亚洲| 国产单亲对白刺激| 久热爱精品视频在线9| 中文字幕人妻丝袜一区二区| 亚洲精品乱久久久久久| 丰满饥渴人妻一区二区三| 亚洲中文av在线| 国产成人一区二区三区免费视频网站| 97人妻天天添夜夜摸| 深夜精品福利| 国产精品亚洲av一区麻豆| 老熟妇乱子伦视频在线观看| 啦啦啦 在线观看视频| 性少妇av在线| 亚洲精品成人av观看孕妇| 国产一区二区三区综合在线观看| 婷婷丁香在线五月| 飞空精品影院首页| 国产成人精品无人区| 久久热在线av| 国产一区有黄有色的免费视频| 中国美女看黄片| 啦啦啦中文免费视频观看日本| 国产亚洲精品第一综合不卡| 久久精品亚洲av国产电影网| 两性午夜刺激爽爽歪歪视频在线观看 | 老司机深夜福利视频在线观看| av在线播放免费不卡| 国产成人影院久久av| 日日摸夜夜添夜夜添小说| 亚洲欧美色中文字幕在线| 久久久久国内视频| 久久精品国产综合久久久| 免费少妇av软件| 桃红色精品国产亚洲av| 精品高清国产在线一区| 蜜桃在线观看..| 免费在线观看黄色视频的| 超碰成人久久| 国产一区二区激情短视频| 成人国语在线视频| 国产一区二区三区视频了| 国产精品av久久久久免费| 国产精品久久电影中文字幕 | 精品亚洲成a人片在线观看| 人成视频在线观看免费观看| 久久久久久久久免费视频了| 美女午夜性视频免费| a在线观看视频网站| 中文字幕高清在线视频| 另类亚洲欧美激情| 国产区一区二久久| 国产一区二区三区视频了| 欧美日韩亚洲综合一区二区三区_| 国产精品免费大片| 乱人伦中国视频| 久久免费观看电影| 少妇猛男粗大的猛烈进出视频| 美女午夜性视频免费| 如日韩欧美国产精品一区二区三区| 9191精品国产免费久久| 91精品三级在线观看| 一本—道久久a久久精品蜜桃钙片| 亚洲五月色婷婷综合| 成人亚洲精品一区在线观看| 国产日韩一区二区三区精品不卡| 亚洲少妇的诱惑av| 久久精品91无色码中文字幕| 成人永久免费在线观看视频 | 9191精品国产免费久久| 中文字幕人妻丝袜一区二区| 国产精品欧美亚洲77777| 少妇粗大呻吟视频| 又大又爽又粗| 久久精品国产综合久久久| 欧美午夜高清在线| 极品人妻少妇av视频| 女同久久另类99精品国产91| 久久国产精品大桥未久av| 黄色丝袜av网址大全| 成人国产一区最新在线观看| 欧美精品人与动牲交sv欧美| 一区二区三区乱码不卡18| 宅男免费午夜| 老司机福利观看| 少妇粗大呻吟视频| 欧美亚洲 丝袜 人妻 在线| 午夜91福利影院| 另类亚洲欧美激情| 日本vs欧美在线观看视频| 欧美成人午夜精品| 国产精品av久久久久免费| 久久国产精品人妻蜜桃| 国产野战对白在线观看| 欧美日韩福利视频一区二区| 中文亚洲av片在线观看爽 | 中文字幕色久视频| 久久人人97超碰香蕉20202| 麻豆av在线久日| 在线观看舔阴道视频| 一进一出好大好爽视频| 成人国产av品久久久| 麻豆av在线久日| 欧美日韩国产mv在线观看视频| 国产精品国产av在线观看| 亚洲国产av新网站| 欧美日韩亚洲高清精品| 91国产中文字幕| 国产精品久久久久成人av| 99精品欧美一区二区三区四区| 国产成人欧美在线观看 | 精品人妻熟女毛片av久久网站| 少妇的丰满在线观看| 美女视频免费永久观看网站| 亚洲一区二区三区欧美精品| 一边摸一边抽搐一进一出视频| 国产欧美亚洲国产| 在线播放国产精品三级| 日韩 欧美 亚洲 中文字幕| bbb黄色大片| 国产精品亚洲一级av第二区| 人成视频在线观看免费观看| 夜夜骑夜夜射夜夜干| 91精品三级在线观看| 国产成人啪精品午夜网站| 老熟女久久久| 精品国产乱码久久久久久小说| 一二三四在线观看免费中文在| 国产精品一区二区在线观看99| 两个人看的免费小视频| 国产精品国产av在线观看| 在线av久久热| 亚洲av成人不卡在线观看播放网| 亚洲天堂av无毛| 九色亚洲精品在线播放| 久久99一区二区三区| 黄色毛片三级朝国网站| 十八禁网站免费在线| 欧美日韩亚洲国产一区二区在线观看 | 亚洲一卡2卡3卡4卡5卡精品中文| 捣出白浆h1v1| 人人妻人人澡人人看| 欧美 亚洲 国产 日韩一| 99re在线观看精品视频| 大型av网站在线播放| 成人影院久久| 日本黄色视频三级网站网址 | 久久国产精品大桥未久av| 自拍欧美九色日韩亚洲蝌蚪91| 一二三四在线观看免费中文在| 中文字幕人妻丝袜一区二区| 精品国产乱码久久久久久小说| 一区二区三区乱码不卡18| 精品视频人人做人人爽| 嫁个100分男人电影在线观看| 脱女人内裤的视频| 乱人伦中国视频| 大型av网站在线播放| 久久热在线av| 亚洲欧美色中文字幕在线| 成人黄色视频免费在线看| 国产精品久久久久成人av| 亚洲男人天堂网一区| 最新在线观看一区二区三区| 国产亚洲一区二区精品| 国产成人精品在线电影| 99热国产这里只有精品6| 亚洲专区字幕在线| 成年动漫av网址| av网站免费在线观看视频| 9色porny在线观看| 91字幕亚洲| 国产伦人伦偷精品视频| 法律面前人人平等表现在哪些方面| 在线观看www视频免费| 亚洲中文av在线| 日韩中文字幕视频在线看片| 黄色成人免费大全| 成年人免费黄色播放视频| 国产精品亚洲av一区麻豆| 一二三四在线观看免费中文在| 99精品欧美一区二区三区四区| 日韩一卡2卡3卡4卡2021年| 色婷婷久久久亚洲欧美| 亚洲,欧美精品.| 国产一区二区在线观看av| 欧美精品啪啪一区二区三区| 美女高潮喷水抽搐中文字幕| 成人影院久久| 国产91精品成人一区二区三区 | 欧美老熟妇乱子伦牲交| 中文字幕最新亚洲高清| 他把我摸到了高潮在线观看 | 十八禁网站免费在线| av天堂久久9| h视频一区二区三区| 首页视频小说图片口味搜索| 涩涩av久久男人的天堂| 在线av久久热| 国产在线观看jvid| 人妻久久中文字幕网| 亚洲成人免费av在线播放| 一区二区av电影网| 十八禁高潮呻吟视频| 日本一区二区免费在线视频| 自线自在国产av| av福利片在线| 欧美亚洲日本最大视频资源| tocl精华| 久久久久久久大尺度免费视频| 亚洲五月婷婷丁香| 一区二区三区激情视频| 中文欧美无线码| 日本黄色视频三级网站网址 | 9色porny在线观看| 亚洲,欧美精品.| 欧美大码av| 两个人免费观看高清视频| 大码成人一级视频| 99re在线观看精品视频| 人妻久久中文字幕网| 国产精品久久久av美女十八| 亚洲国产欧美在线一区| 成年动漫av网址| 高清av免费在线| 天天影视国产精品| 高清av免费在线| 精品一区二区三卡| 欧美国产精品一级二级三级| 1024香蕉在线观看| 91老司机精品| 一本一本久久a久久精品综合妖精| 国产aⅴ精品一区二区三区波| 国产高清激情床上av| 欧美日韩亚洲国产一区二区在线观看 | 亚洲国产欧美日韩在线播放| 久久午夜亚洲精品久久| 黄色 视频免费看| 高清毛片免费观看视频网站 | 亚洲情色 制服丝袜| 1024香蕉在线观看| 色视频在线一区二区三区| 日韩欧美三级三区| 精品少妇久久久久久888优播| www.999成人在线观看| 精品少妇久久久久久888优播| 一二三四社区在线视频社区8| 高潮久久久久久久久久久不卡| 亚洲一码二码三码区别大吗| 男女床上黄色一级片免费看| 窝窝影院91人妻| 国产不卡av网站在线观看| 亚洲精品国产色婷婷电影| kizo精华| 色视频在线一区二区三区| 捣出白浆h1v1| 丝袜美足系列| 99国产综合亚洲精品| 一区二区三区乱码不卡18| 91av网站免费观看| 国产精品久久久久久精品电影小说| 夜夜爽天天搞| 丰满饥渴人妻一区二区三| tocl精华| 韩国精品一区二区三区| 一本—道久久a久久精品蜜桃钙片| 亚洲国产毛片av蜜桃av| 夫妻午夜视频| 欧美黑人欧美精品刺激| 国产精品电影一区二区三区 | 中国美女看黄片| 亚洲欧美日韩另类电影网站| 亚洲av国产av综合av卡| 午夜两性在线视频| 久9热在线精品视频| 中文字幕制服av| bbb黄色大片| 久久精品国产a三级三级三级| 无人区码免费观看不卡 | 黑人猛操日本美女一级片| 精品人妻在线不人妻| svipshipincom国产片| 久久毛片免费看一区二区三区| 亚洲欧洲精品一区二区精品久久久| 国产一卡二卡三卡精品| 看免费av毛片| 老司机影院毛片| 女人爽到高潮嗷嗷叫在线视频| 叶爱在线成人免费视频播放| 亚洲一区中文字幕在线| 最新的欧美精品一区二区| 日韩中文字幕欧美一区二区| 最黄视频免费看| 天天躁夜夜躁狠狠躁躁| 久久久久久亚洲精品国产蜜桃av| 99re在线观看精品视频| av天堂久久9| 亚洲成av片中文字幕在线观看| 久久99热这里只频精品6学生| 日本黄色日本黄色录像| 亚洲国产欧美网| 欧美日韩亚洲综合一区二区三区_| 午夜福利影视在线免费观看| 精品福利永久在线观看| 一区二区三区激情视频| bbb黄色大片| 精品亚洲成a人片在线观看| 可以免费在线观看a视频的电影网站| 欧美在线黄色| 最新美女视频免费是黄的| 国产欧美日韩精品亚洲av| 久久国产亚洲av麻豆专区| 国产av精品麻豆| 热99国产精品久久久久久7| 婷婷成人精品国产| 欧美精品一区二区大全| 在线观看免费视频日本深夜| 中文字幕人妻熟女乱码| 日本五十路高清| 777米奇影视久久| 久久精品熟女亚洲av麻豆精品| 黑人猛操日本美女一级片| 一本一本久久a久久精品综合妖精| av国产精品久久久久影院| 国产男女超爽视频在线观看| 国产亚洲精品一区二区www | 精品欧美一区二区三区在线| 国产精品国产av在线观看| 色视频在线一区二区三区| 久久久国产精品麻豆| 亚洲成人手机| 亚洲七黄色美女视频| 久热爱精品视频在线9| 日本av免费视频播放| 亚洲av成人不卡在线观看播放网| 一本—道久久a久久精品蜜桃钙片| 波多野结衣av一区二区av| 熟女少妇亚洲综合色aaa.| 热99久久久久精品小说推荐| 日本av手机在线免费观看| 国产在线免费精品| 亚洲va日本ⅴa欧美va伊人久久| 最近最新免费中文字幕在线| 老司机靠b影院| 精品一区二区三区四区五区乱码| 久久国产精品男人的天堂亚洲| 黑人猛操日本美女一级片| 亚洲av美国av| 大片免费播放器 马上看| 中文字幕人妻丝袜制服| 久久国产精品影院| 麻豆国产av国片精品| 老熟女久久久| 法律面前人人平等表现在哪些方面| 一本一本久久a久久精品综合妖精| 国产免费视频播放在线视频| 男女免费视频国产| 少妇 在线观看| 欧美日韩视频精品一区| 十八禁高潮呻吟视频| 国产精品久久电影中文字幕 | 日韩视频一区二区在线观看| a级毛片黄视频| 高清毛片免费观看视频网站 | 怎么达到女性高潮| 久久中文字幕一级| 不卡一级毛片| 老司机午夜十八禁免费视频| 国产免费av片在线观看野外av| 亚洲色图综合在线观看| 满18在线观看网站| 国产成人欧美| 自拍欧美九色日韩亚洲蝌蚪91| 美女高潮喷水抽搐中文字幕| 亚洲欧美一区二区三区久久| 黄频高清免费视频| 精品国产超薄肉色丝袜足j| 亚洲成a人片在线一区二区| 日韩欧美一区二区三区在线观看 | 又黄又粗又硬又大视频| 91国产中文字幕| 精品亚洲成a人片在线观看| 亚洲精品国产色婷婷电影| 色94色欧美一区二区| 亚洲色图av天堂| 一区在线观看完整版| 天天躁狠狠躁夜夜躁狠狠躁| 天堂8中文在线网| 少妇裸体淫交视频免费看高清 | 高清av免费在线| 欧美日本中文国产一区发布| 久久久久网色| 亚洲人成电影免费在线| 免费少妇av软件| xxxhd国产人妻xxx| 国产成人精品无人区| 成人av一区二区三区在线看| av网站在线播放免费| 亚洲熟妇熟女久久| 欧美精品啪啪一区二区三区| 久久香蕉激情| 一区福利在线观看| 少妇猛男粗大的猛烈进出视频| 啪啪无遮挡十八禁网站| 天天躁夜夜躁狠狠躁躁| 国产精品 欧美亚洲| 91成年电影在线观看| 少妇 在线观看| 老司机午夜十八禁免费视频| 亚洲精品一二三| 国产亚洲精品一区二区www | 午夜免费鲁丝| 久久人妻熟女aⅴ| 欧美亚洲日本最大视频资源| 人妻久久中文字幕网| 欧美激情极品国产一区二区三区| 一边摸一边抽搐一进一出视频| 欧美日韩一级在线毛片| 日韩欧美一区视频在线观看| 丝袜美足系列| 高清毛片免费观看视频网站 | 国产成+人综合+亚洲专区| 久久性视频一级片| 国产精品麻豆人妻色哟哟久久| 高潮久久久久久久久久久不卡| 满18在线观看网站| 亚洲全国av大片| 高清在线国产一区| 叶爱在线成人免费视频播放| 妹子高潮喷水视频| 一级片'在线观看视频| av有码第一页| 久热爱精品视频在线9| 免费不卡黄色视频| 他把我摸到了高潮在线观看 | 亚洲一区二区三区欧美精品| 国产男靠女视频免费网站| 国产在线精品亚洲第一网站| 亚洲av日韩在线播放| 天天躁夜夜躁狠狠躁躁| 欧美精品av麻豆av| 中文字幕另类日韩欧美亚洲嫩草| 国产成人欧美| 午夜成年电影在线免费观看| 成在线人永久免费视频| 国产无遮挡羞羞视频在线观看| 人妻久久中文字幕网| 欧美日韩中文字幕国产精品一区二区三区 | 男人舔女人的私密视频| 亚洲精品国产色婷婷电影| 少妇裸体淫交视频免费看高清 | 99精品在免费线老司机午夜| 亚洲国产欧美一区二区综合| 免费看a级黄色片| 人人妻人人添人人爽欧美一区卜| 人成视频在线观看免费观看| 国产在线一区二区三区精| 久久久久久久久免费视频了| 国产成人精品久久二区二区91| www.自偷自拍.com| 国产精品一区二区在线不卡| 国产免费福利视频在线观看| 亚洲全国av大片| 成人18禁高潮啪啪吃奶动态图| 男女免费视频国产| 日韩制服丝袜自拍偷拍| 欧美久久黑人一区二区| 男女下面插进去视频免费观看| 久久精品国产a三级三级三级| 人人澡人人妻人| 人成视频在线观看免费观看| 日本五十路高清| xxxhd国产人妻xxx| 麻豆av在线久日| 丁香欧美五月| 精品少妇内射三级| 成年人黄色毛片网站|