在线免费观看成年人视频-在线免费观看国产-在线免费观看国产精品-在线免费观看黄网站-在线免费观看精品

產品分類

當前位置: 首頁 > 工業電氣產品 > 端子與連接器 > 線路板連接器 > FFC連接器

類型分類:
科普知識
數據分類:
FFC連接器

最小化運行費用的時變重構全局優化算法

發布日期:2022-04-17 點擊率:75

<script var cpro_id = "u1457042";

<iframe id="iframeu1457042_0" fcam?rdid=1457042&dc=2&di=u1457042&dri=0&dis=0&dai=3&ps=345x362&dcb=BAIDU_SSP_define&dtm=BAIDU_DUP_SETJSONADSLOT&dvi=0.0&dci=-1&dpt=none&tsr=0&tpr=1459711229288&ti=%E6%9C%80%E5%B0%8F%E5%8C%96%E8%BF%90%E8%A1%8C%E8%B4%B9%E7%94%A8%E7%9A%84%E6%97%B6%E5%8F%98%E9%87%8D%E6%9E%84%E5%85%A8%E5%B1%80%E4%BC%98%E5%8C%96%E7%AE%97%E6%B3%95_%E7%94%B5%E6%B0%94%E8%87%AA%E5%8A%A8%E5%8C%96%E6%8A%80%E6%9C%AF%E7%BD%91&ari=1&dbv=0&drs=1&pcs=645x335&pss=970x346&cfv=0&cpl=22&chi=50&cce=true&cec=gbk&tlm=1402381900&ltu=http%3A%2F%2Fwww.dqjsw.com.cn%2Fdiangongdianzi%2Fdianlidiangong%2F3074.html&ecd=1&psr=1366x768&par=1366x728&pis=-1x-1&ccd=24&cja=false&cmi=34&col=zh-CN&cdo=-1&tcn=1459711229&qn=9dc1042a8b99bde6&tt=1459711229266.143.215.216" vspace="0" hspace="0" marginwidth="0" marginheight="0" scrolling="no" style="border:0; vertical-align:bottom;margin:0;" allowtransparency="true" align="center,center" width="200" height="200" frameborder="0">

1  引言
   
配電網絡重構通過改變網絡拓撲結構來獲得運行效益最優的運行方式,其研究內容可分為兩大類:
       第一類是針對恒定負荷的靜態重構,主要方法有啟發式方法,如支路交換法[1]、最優流模式法等;傳統數學優化方法,如線性整數規劃法、分支定界法等;智能計算方法,如模擬退火法、神經網絡法、遺傳算法等。第二類是時變重構(也稱動態重構),是在負荷不斷變化的情況下求解一段時間內的重構方案。主要方法有啟發式方法[2, 3],啟發式方法與智能計算方法相結合的混合算法[4],等等。
   
實際配電系統負荷變化頻繁,靜態重構只給出開關設備一個時刻的動作方案,各時刻間的優化方案缺乏聯系,所以缺乏實用性。靜態重構屬于大規模、多目標、非線性、混合整型、組合優化問題,屬于NP難問題,顯然,針對給定時間區間的時變重構問題時間復雜性和空間復雜性更高。目前見諸文獻的方法根據啟發式規則只提供了局優解。
   
本文以中壓配電網絡為研究對象,提出了以運行費用最小為目標的時變重構全局優化算法。
2  數學模型
   
本文算法的目標函數可表示為

式中  FNBN分別為給定時間區間內的運行費用(包括開關操作費用和電能損耗費用)、可操作的開關總數和所劃分的時間段總數;xk,i為開關k在第i時段的狀態,取值為0或1,分別為開關處于斷開位置和閉合位置;xk,0為開關k的初始狀態;PLi為第i時段系統的有功功率損耗,kW;ti為第i時段的長度, h;CBkCLi分別為開關k操作一次的費用,元和第i時段的網損電價,元/kWh。
   
重構后的網絡應滿足潮流方程約束、各種安全約束和放射狀運行約束。此外,時變重構問題還必須滿足開關操作次數約束,即在所研究的時段內滿足式(2)所描述的系統中所有可操作開關總的操作次數限制和式(3)所描述的單臺開關的操作次數限制。

       式中  NANSk分別為給定時間區間內最大允許的總的開關操作次數和開關k最大允許的操作次數。
3  時變重構全局優化算法
3.1  算法步驟

   
為了在給定時間區間內實現中壓配電網的運行費用最小,本文算法進行如下操作:
  
(1)計算與識別網絡的初始信息,包括節點電壓、支路功率、有功損耗、開關狀態,等等;
  
(2)將給定的時間區間劃分為多個時間段。在每個時間段內,根據負荷預測的結果將負荷近似為恒功率負荷。以網損最小為目標,使用核心模式遺傳最短路算法求解靜態重構子問題,得出多個候選的可行放射狀網絡結構。
 
 (3)根據上步得出的各個時段的可行放射狀網絡結構,應用動態規劃,求解運行費用最小的全局最優時變重構策略。
   
可見,本文算法包含核心模式遺傳最短路算法和動態規劃兩個核心,下面分別予以說明。
3.2  CSGSA求解靜態重構子問題
    用無向元件網絡圖G來表示所研究的中壓配電網,BS LS分別表示支路集合與負荷節點集合。大部分見諸文獻的重構方法是針對BS進行操作的,試圖找到BS的一個子集BS'?BS,使得BS中的所有支路組成一個滿足各種約束條件的放射狀運行網絡,這類方法被稱為組合支路法。CSGSA基于組合負荷的思想,其研究對象為LS而非BS,通過為LS中的每個負荷分別尋找供電路徑來逐步獲得放射狀網絡結構。負荷尋路的先后次序決定結果網絡,CSGSA就是在使用遺傳算法搜尋最優的負荷尋路次序的過程中并行得出多個可行的放射狀網絡結構。CSGSA以網損最小為目標,算法流程如圖1所示。
       圖中第(1)步為將LS中的負荷表示成負荷串的形式,其排列次序即為負荷尋路的先后次序;

    第(3)步為形成有向網絡拓撲圖
   
第(4)步為使用最短路法在GSqi中搜尋當前負荷的最佳供電路徑。所用弧權值公式定義為

    式(4)的作用是使尋路過程中網損增量相對較小,剩余容量相對較大的弧優先被選中,從而平衡結果網絡各個支路的負荷,避免由于先尋路負荷的供電路徑過于集中而導致的后尋路負荷無法找到合法供電路徑的現象;
   
第(5)步至第(7)步的理論基礎是遺傳算法的模式定理。研究發現,根據Sqi中負荷排列的先后次序可將LS分解為2個互補的子集LS¢SqiLS2Sqi,其中,LSSqi中的負荷被稱為核心模式負荷,它們對放射狀網絡結構的生成起著關鍵性的影響。LS2Sqi中的負荷被稱為非核心模式負荷,它們通過改變GSqi中的支路權值來影響后續負荷的尋路結果,因此對網絡結構的生成也有影響,但該影響是次要的。例如,對Sqi中的負荷Li尋路,將得到供電路徑P。假設P上存在還未使用最短路法尋路的負荷Lα1, ···, Lαα,根據配電網絡的放射狀約束易知,其供電路徑已由Li確定,因此有:LS2Sqi中的負荷混雜形成染色體Sqi,先尋路負荷對后尋路負荷供電路徑的搜索有著程序很難控制的復雜影響,因此該編碼方式得出的染色體各基因的關系是強耦合的,具有高階、長定義長度的特點,不適合遺傳操作。
    第(5)步實現為Sqi尋路的同時,將LS分解為LS'SqiLS2 Sqi,從而突出核心模式負荷的作用,屏蔽非核心模式負荷的影響。具體操作如下:建立空的核心模式染色體,為Sqi中的負荷依次尋路的同時,剔除相關的非核心模式負荷,同時向中以從首端至尾端的順序添加核心模式負荷,以從尾端至首端的順序添加相關的非核心模式負荷,最后用所得到的核心模式染色體替換Sqi來參加遺傳運算。所得S'qi的模式如下所示:

式中  *為無需使用最短路法尋路的非核心模式負荷。

   
的首段Lg1,…,LggLS'Sqi由核心模式負荷組成,尾段由非核心模式負荷組成,負荷排列的先后順序由其在Sqi中的相對位置決定。可見,該類染色體都具備低階、短定義長度的特點,所以,通過選擇、交叉和變異等遺傳運算,一些較差的核心模式逐步被淘汰,一些較好的核心模式逐步被遺傳和進化,最終可得到問題的最優解。這種方法不但大大改善了算法的收斂性,提高了搜尋最優解的能力,同時省去了為非核心模式負荷使用最短路法尋優的過程。通常,核心模式負荷數g遠小于非核心模式負荷數β,因此,計算速度得到大幅度提高。
   
通過計算得出的候選放射狀網絡結構成為動態規劃求解全局最優時變重構策略的前提與基礎。
3.3  以動態規劃求解全局最優時變重構策略
3.3.1 動態規劃數學模型

   
動態規劃是運籌學的一個分支,是求解多階段決策過程最優化的數學方法。本文把給定的時間區間劃分為若干階段,將時變重構問題轉化為多階段決策過程,然后使用動態規劃求解,其步驟為:
  
(1)以時間劃分階段。將所研究的時間區間劃分為N個時段,各時段長度可以相等,也可不等;
  
(2)選擇可行放射狀網絡結構作為動態規劃中的狀態(由網絡中可操作開關的開合位置來表示)。第i時段的可達狀態集合記為Si
  
(3)時段間狀態轉移決策。Ui為第i時段狀態處于si的決策變量,Di為第i時段的允許決策集合, ui(si)∈Di(si),Ti為狀態轉移函數。時變重構問題的初始狀態已知,故根據siui來確定si-1,有

   (4)指標函數和最優值函數。
   
選擇開關操作費用和電能損耗費用之和為指標函數,其表達式為

邊界條件為f0(s0)=0,即設定s0為網絡的初始狀態,作為第0階段,其費用為0。第1時段至第N時段為所研究的時間區間。
   
從第0階段開始遞推求解,記錄每個si的最優值函數及相應的si-1,最后從第N時段選出最優的sN并回推得出全局最優的重構策略。
3.3.2 開關操作
次數約束的處理
   
為了限制開關操作次數,防止由于開關操作次數約束導致的無解或采納局優解,應使某些相鄰時段可以存在相同的網絡結構,因為如果動態規劃在相鄰的時段選擇了相同的網絡結構,則表示時段間無需任何開關操作。動態規劃各個階段的狀態由CSGSA計算獲得,都是網損相對較小的網絡結構,所以本文采取了狀態鎖定措施,目的是提供網損相對較大但所需開關操作次數較少的網絡結構,也就是保證相鄰的若干時間段中存在相同的可行放射狀網絡結構。狀態鎖定措施的具體步驟是:將當前時段的狀態加入相鄰的下一時段中,判斷該狀態是否已經存在,并計算狀態是否可行,若該狀態已經存在或是狀態可行,則繼續向后面相鄰的時段添加該狀態,若狀態為不可行或是已經到達末時段,則終止當前狀態的添加過程。狀態鎖定豐富了各個時段的狀態數量,為限制開關操作次數提供了可能。
   
在動態規劃尋優過程中,如果某階段狀態給定,則此前各階段狀態決定了到目前為止已經操作的開關次數,若此時考慮開關操作次數約束,則意味著此前各階段狀態將影響該階段以后狀態的選擇,這將導致狀態變量不滿足構造動態規劃模型所必需的無后效性(馬爾科夫性)條件。雖然采取適當增多狀態變量的方法(如增加已操作開關數作為狀態變量),總能把過程變為無后效的,但這會大大增加算法的復雜性。考慮到本文算法的目標是最小化運行費用,其中包括了開關操作的費用,這就等價于隱性地限制了開關操作次數。加之狀態鎖定措施為限制開關操作次數提供了可能,因此,動態規劃遞推尋優過程中不考慮開關操作次數是否越限。當尋優過程結束后,選擇運行費用最小的重構策略,這時再判斷是否滿足開關操作次數約束,若滿足約束,則完成計算,得出全局最優時變重構策略。
   
動態規劃計算得出了多個可行重構策略,經過遞推尋優過程之后,這些可行策略的運行費用是非常接近的,所以當費用最小解違反開關約束時,可以直接從所得的眾多策略中選出滿足開關約束的費用最小策略作為最終結果。開關操作次數越限往往只對應于給定的開關可操作次數非常少的極端情況。大量計算表明,本文算法所得費用最小解所需開關操作次數很少,一般情況下能夠滿足開關操作次數的約束條件。因為一方面,CSGSA通過全局尋優和局部尋優,能夠快速得到每個時段的最優網絡結構和大量的次優網絡結構,狀態鎖定措施進一步豐富了各個時段可行網絡結構的數量,使得每個時段都包含了足夠的可行網絡結構以供選擇;另一方面,動態規劃本身具備了全局搜索能力。所以,本文算法能夠得到最優的和大量次優的時變重構策略,保證了算法的全局最優性。
4  算法復雜性分析

    CSGSA內嵌Dijkstra最短路法和遺傳算法,其計算復雜度是O(p3),其中p為參與計算的節點數量。核心模式負荷概念的引入使得p減小,因此有效地提高了算法的求解速度。
   
動態規劃的計算復雜度與時間段數和每個時段的候選網絡結構數密切相關。用W表示N個時段的候選網絡結構總數,則動態規劃算法時間復雜性為O((N-1)W2),空間復雜性為O(NW)。本文只選擇了一個動態規劃狀態變量,狀態總數和時間段總數不多,所以計算速度非常快。如果在動態規劃尋優的過程中考慮開關操作次數約束,那么只有增加動態規劃狀態變量的個數才能維持解的最優性,最簡單的方法是再增加開關操作次數作為狀態變量。算法的時間復雜性和空間復雜性分別變為O(W2N-1)和O(WN),都是關于N的指數函數。對比可見,本文算法在最壞情況下的時空復雜性遠遠小于在尋優過程中直接考慮開關操作約束的情況,所以可有效地求解時變重構問題。
5  算例分析
       采用如圖2所示的12.66kV、33節點、5網孔系統[1],圖中給出了每條線段的編號。將系統中的負荷歸為3大類:大型工商業負荷,中小工商業負荷和居民負荷。每類負荷的變化遵循相應的負荷曲線,文[1]所給的負荷數據作為所研究時間區間內各節點的峰荷。圖3給出的負荷曲線取自實際系統1997年1月1日~10日(星期三~下星期五)的數據。設CBk為7元/次,CLi為0.4元/kWh,ti全部取1h。

    大量計算表明,針對不同的初始結構和時間區間,算法總能找到從初始結構過渡到最優結構的費用最小重構方案。當網絡負載率不高時,網絡調整到最優結構后,如果負荷的模式如類型、峰荷大小、負荷曲線等沒有太大的變化,則不需再進行任何重構操作,網絡結構將保持最優,這為季節性重構提供了參考。例如,在上述網絡和負荷模式下,最優網絡結構選擇線路7,9,14,28,36作為聯絡線。
       設初始網絡結構選擇33~37作為聯絡線,使用本文算法進行時變重構優化計算,第1種情況為將24h作為給定時間區間進行計算,優化結果作為下次計算的初始網絡結構,連續計算3天;第2種情況為將3天作為給定時間區間進行優化計算。兩種情況的計算結果見表1。可見,給定的時間區間不同,可能導致優化結果的不同。如果初始網絡結構與最優結構相差較大,由式(1)可知,時間區間長,則開關操作費用在總費用中的比例小,優化結果趨于最優網絡結構;時間區間短,則算法將逐步向最優結構過渡,以保證當前時間區間內的費用最小。表1即說明了2種情況到達最優結構的路徑不同,第1種情況的前24h屬于過渡階段,后48h重構為最優結構。第2種情況只在0點重構一次便使得72h內運行費用最小。圖4為前24h的網損比較。圖中還給出了不限制開關操作次數的網損最小結果,網損為4610.71kW,72h里需要進行26次開關操作,費用為2026.28元。可見,本文算法所得結果的網損與網損最小值非常接近,但費用卻顯著減少。以上雖然對不同時間區間的優化結果進行了比較,但實際應用中,若網絡負載率不高且當前負荷模式維持足夠長時間時,將選擇一步重構到當前負荷模式下的最優網絡結構,而不必經過若干過渡。如果網絡的負載率很高或負荷變化劇烈,那么網絡結構應進行相應調整以保證滿足電壓或容量等約束條件。表2給出了當每個節點的負荷都增長19%時的重構策略,同時給出相同情況下的網損最小結果。

    如果開關可操作次數很少,無法重構到最優的網絡結構,則算法將得出滿足開關操作次數約束的運行費用次優解。如表2所示的72h優化結果需要8次開關操作,如果限定每臺開關只能操作2次,總的開關操作不超過4次,那么算法所得網絡結構的聯絡線是11,28,33,34,36,從初始狀態只需4次開關操作即可實現,網損為4698.25kW,費用為1907.30元。此外,若使開關操作費用為零,則算法轉化為電能損失最小重構問題;若只考慮單一時段,則算法轉化為靜態重構問題;若只考慮單一時間點,同時開關操作費用為零,則算法轉化為傳統的靜態網損最小重構問題。實際應用中可能遇到最優解不唯一或者有其它一些約束需要考慮的情況,這時,算法將動態規劃尋優所得的所有可行重構策略由優至劣排序,然后全部輸出,由用戶來選擇最終的方案。可見本文算法非常實用,應用范圍非常廣泛。

6  結語
   
實際配電系統負荷變化頻繁,重構方案也應相應變化,本文給出了時變重構的全局優化算法,可在所研究的時間區間內實現運行費用最小化。算法主要由CSGSA和動態規劃組成,并對開關操作次數約束進行了處理,使之不影響動態規劃計算結果的全局最優性。大量計算表明:本文算法計算速度快,解的質量高,可用于配電網絡的實時調度,可在降低運行費用的同時提高供電質量與供電可靠性,也可將之用于規劃和離線分析等方面。

下一篇: PLC、DCS、FCS三大控

上一篇: 索爾維全系列Solef?PV

推薦產品

更多
国产精品欧美一区二区三区| 国产精品无码无卡在线观看久| 肥老熟妇伦子伦456视频| 国产乱码精品一品二品| 久久精品久久久久久久精品| 妺妺窝人体色聚窝窝www毛片| 日本XX爽21护士| 亚洲AV无码AV制服丝袜在线| 在线 | 一区二区三区| 成人无码影片精品久久久| 国产在线观看无码免费视频| 两个奶头被吃到高潮什么感觉| 日本免费人成视频在线观看| 亚洲AV成人无码精品网站老司机| 野花香高清视频在线观看免费| JAPAN强要VIDEOD警妞| 国产乱人伦中文无无码视频试看| 久久天天躁夜夜躁狠狠躁| 人与牲动ZZZXXXⅩ0000| 亚洲AV无码国产精品色午| AV天堂久久天堂AV色综合| 国产精品无卡毛片视频| 免费看漫画在线成人漫画| 玩弄chinese丰满人妻| 一边亲着一面膜下奶韩剧| 成人无码区免费A∨电影| 精品亚洲成在人线AV无码| 人妻人人澡人人添人人爽人人玩 | 绯色精品人妻av一区二区 | 亚洲精品国产成人AV蜜臀| 啊轻点灬大巴太粗太长视频| 好爽又高潮了毛片| 日本XXXX裸体XXXX在线观| 亚洲色偷偷偷综合网 | 他的粗大把她捣出白沫| 一区二区三区鲁丝不卡麻豆| AAAAA级少妇高潮大片免费看| 国产成人a在线看片色欲AV| 皇上御花园HLH| 朋友的人妻的滋味BD高清中文| 小莫骚麦歌曲播放MP3| 88国产精品欧美一区二区三区| 国产成人无码免费视频79| 久久精品这里热有精品| 人人爽人人爽人人片AV| 亚洲精品国产成人99久久6| 菠萝视频高清视频在线7| 精品国内自产拍在线观看| 人人妻人人做人人爽| 亚洲精品无码永久在线观看你懂的| 锕锕锕锕锕锕~好湿WWW| 激情综合一区二区迷情校园| 日韩精品无码免费专区网站| 亚洲色偷偷综合亚洲AV| 疯狂做受ⅩⅩⅩⅩ高潮高清视频| 久久综合久久鬼色| 无码熟妇人妻AV在线影片免费 | 极品粉嫩小泬50P| 欧美性狂猛BBBBBBXXXX| 亚洲丰满性熟妇ⅩXXOOO太阳| 被男人吃奶很爽的毛片 | 欧美顶级METART裸体全部自| 亚洲AV无码一区二区三区系列| 超碰97人人做人人爱综合| 精品国产乱子伦一区二区三区| 日本最新免费二区三区| 一个吃我奶头两个舔我下面| 国产精品自在在线午夜出白浆| 欧美日韩一区二区综合 | AAA级久久久精品无码片| 狠狠躁夜夜躁人人爽天天天天97 | 日韩无码AV一区二区| 真人无码作爱免费视频网站| 国产人成无码视频在线观看| 日本高清视频色WWWWWW色| 制服丝袜中文字幕在线| 好儿子妈妈今天就是你的女人| 色婷婷综合久久久久中文| 97久久精品人人澡人人爽| 精品无码国产污污污免费网站| 特级毛片A级毛片免费观看网站| AV 无码 在线 观看| 久久婷婷五月综合97色直播 | 国产精品Ⅴ无码大片在线看| 欧美性猛交XXXXXⅩXX| 一对一大尺度聊天APP安装| 国产在线拍偷自揄拍无码| 熟妇人妻AV中文字幕老熟妇| 啊轻点灬大JI巴太粗太长了网站| 麻花传媒MD0076沈芯语在线| 亚洲精品美女久久7777777 | 好了AV第四综合无码久久| 图片区小说区另类春色| 成人AV鲁丝片一区二区免费| 欧美成人精品视频在线不卡 | 精品一区二区成人精品| 天堂…在线最新版在线| 顶级大但人文艺术巫| 人妻少妇精品视中文字幕国语| 中文字幕AV无码一二三区电影| 精品无码国产污污污免费| 亚洲AV无码有乱码在线观看| 国产精品一区二区久久| 日本不卡一区二区三区| JEAⅠOUSVUE成熟少归A| 男人扒开女人内裤强吻桶进去 | 精JAVAPARSER乱偷| 亚洲AV无码国产综合专区 | 久久精品国产亚洲ΑV忘忧草| 小少妇ASS浓PICS| 国产精品污WWW一区二区三区| 色综合热无码热国产| ASS年轻少妇BBWPIC精品| 久久夜色撩人精品国产| 亚洲中文字幕AV不卡无码| 精品人妻无码区二区三区密桃| 亚洲AV永久无码一区二区三区| 国产免费永久精品无码| 无码一区在线视频| 国产精品人人做人人爽| 无码人妻久久久一区二区三区免费| 国产A级三级三级三级| 熟妇阿 HD中文电影| 成人午夜免费无码福利片| 色综合天天综合欧美综合| 夫妻二人体验交换夫妻的后果| 日韩一中文字无码不卡| 俄罗斯女人与马Z00Z视频| 日韩视频中文字幕精品偷拍| 俄罗斯ZOOM与人性ZOOM| 熟妇人妻AV中文字幕老熟妇 | AV无码AV高潮AV喷吹免费| 欧美精品一区二区三区人妻久久久| √天堂中文WWW官网在线| 女人张开腿扒开内裤让男生桶| 999国内精品永久免费视频| 欧美牲交黑粗硬大在线视频| 爆乳一区二区三区无码| 全免费A级毛片免费看无码| 宝贝我不想带小雨伞了| 三上悠亚公侵犯344在线观看| 国产99视频精品免费视频36| 我把英语课代表按在桌子上抄| 国产丝袜无码一区二区视频| 亚洲高清乱码午夜电影网| 久久久久久精品免费免费4K| 中文无码乱人伦中文视频在线| 女性裸体无遮挡啪啪网站| 被老外添嫩苞添高潮NP视频 | 无码人妻精品一区二| 黑人双人RAPPER的特点| 亚洲愉拍自拍欧美精品| 男男GAYGAY无套GAY无套| АⅤ天堂中文在线网官网| 色综合久久久久综合一本到桃花网| 国产高清在线A视频大全| 亚洲AV永久无码一区| 极品白嫩的小少妇| 尤物蜜芽国产成人精品区| 欧美粗大猛烈老熟妇| 大杳蕉在线影院在线播放| 无码精品一区二区三区视频色欲网 | 干了快生了的孕妇| 亚洲AV无码成人精品区瑜伽裤| 久久精品国产亚洲AV麻豆网站| 亚洲综合激情七月婷婷| 欧美黑人猛XXxXX内射| 暗交小拗女一区二区三区| 天天摸天天碰天天添中文无码| 皇帝在御花园进入贵妃的小说 | 国产精品女人呻吟在线观看| 亚洲日韩精品一区二区三区| 免费看美女脱精光的网站| 成人国产亚洲精品A区| 无码专区一VA亚洲V天堂| 久久精品国产亚洲AV麻豆AⅤ| HD女人奶水授乳MILK| 无码国产偷倩在线播放| 久久国产精品无码一区二区三区| 97国产精华最好的产品价格| 人妻系列无码专区AV在线| 国产精品久久777777| 野花日本大全免费观看2019| 人成午夜免费大片| 国产亚洲无日韩乱码| 中文字幕在线观看| 双腿吊起揉捏花蒂| 久久777国产线看观看精品| GAYⅩXX小奶受GV浪小辉| 特级毛片内射WWW无码| 老湿机香蕉久久久久久| 肥水不流外人田小说| 亚洲人成无码网站| 人妻少妇精品无码专区漫画| 国产艳妇AV在线出轨| 99久久精品日本一区二区免费| 无码熟妇人妻在线视频| 久久久久精品国产99久久综合| 丰满少妇高潮惨叫久久久| 亚洲熟女乱综合一区二区|