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

    農(nóng)田輸水渠(管)道組合設計研究
    ——基于Python遺傳算法

    2020-10-19 12:40:38徐彬冰孫梟沁佘冬立
    節(jié)水灌溉 2020年10期
    關鍵詞:斗渠遺傳算法效益

    徐彬冰,孫梟沁,李 麗,王 琴,佘冬立

    (1.江蘇農(nóng)墾集團,南京 210008;2.河海大學農(nóng)業(yè)科學與工程學院,南京 210098;3. 南京江地土地開發(fā)咨詢服務有限責任公司,南京 210024;4. 河海大學設計研究院有限公司,南京 210098)

    0 引 言

    在糧食大量需求背景下,我國提出了嚴守18億畝耕地紅線的政策,并將全國總用水量的65%左右用于農(nóng)業(yè)灌溉[1]。而目前,我國人均水資源約為世界水平的25%,人均土地資源約為世界水平的40%,人地資源和人水資源矛盾嚴重[2]。近年來由于建設占地、生態(tài)退耕等因素,全國耕地總面積減少,且耕地后備資源日漸不足;此外由于在輸配水技術和管理技術等方面的問題,我國灌溉水利用效率低,灌排溝渠等基礎設施占地率較高。因此,如何提高灌溉水利用效率、增加耕地面積和糧食產(chǎn)量,以滿足國家糧食安全的要求,成為自然資源領域的重要議題。

    為解決這一問題,其一就是要從改善灌溉輸水條件,減少基礎設施占地率。相關研究指出,采用防滲渠道或管道灌溉能夠提高灌溉水利用系數(shù),且占地面積少,灌溉及時、高效,從而實現(xiàn)增產(chǎn)[3,4]。但是傳統(tǒng)的灌溉渠道橫斷面和管道的設計主要采用試算法或圖解法,這些方法不僅計算量大而且精度不易控制。近年來,隨著計算機技術迅速發(fā)展,相對成熟的算法已經(jīng)應用于解決傳統(tǒng)灌溉渠(管)道設計的問題[5-7]。在這些算法中,遺傳算法是模擬生物界遺傳進化而構造出來的一種算法,以概率選擇為主要手段,不需關心問題的內(nèi)在規(guī)律,過程簡單,能夠?qū)崿F(xiàn)全局隨機搜索,能夠很好解決目前灌溉渠(管)道設計的不足。然而,前期的模型研究主要集中于渠道輸水過水性能[6,8],而針對灌溉渠(管)道系統(tǒng)投資研究較少,多級渠道組合優(yōu)化也鮮有研究?;诖耍疚囊越K省臨海農(nóng)場為研究對象,基于PYTHON遺傳算法建立以年需投資費用最低和年凈效益最大的灌溉渠(管)道的優(yōu)化模型,并確定兩級田間灌溉渠(管)道最優(yōu)組合形式,為灌溉渠(管)道的設計及組合優(yōu)化提供新的方法。

    1 模型建立

    本文選擇現(xiàn)澆混凝土梯形渠道、現(xiàn)澆混凝土U形渠道和PE暗管進行農(nóng)田輸水渠(管)道設計組合研究,3種渠(管)道的斷面形式如圖1所示。結合斷面示意圖與明渠均勻流公式,得到梯形渠道水力最佳斷面和U形渠道斷面部分參數(shù)計算方法[9],并建立以年需投資費用最低和年凈效益最大為目標的優(yōu)化模型。

    圖1 渠(管)道斷面示意圖Fig.1 Section diagram of the channel and the pipeline

    灌溉渠(管)道設計時,一方面要降低工程造價,增加效益,另一方面又要滿足設計的約束條件,這兩者之間是互相影響的。因此,如何選擇合適的設計參數(shù)使得工程造價最低,效益最大,并滿足設計的約束條件是灌溉渠(管)道優(yōu)化設計需要解決的問題。

    1.1 工程建設費用

    1.1.1 材料費用

    灌溉渠(管)道采用不同的斷面形式、不同的材料,會產(chǎn)生不同的工程量,進而影響工程造價。因此,本文選擇3種不同斷面形式和建筑材料的灌溉渠(管)道作為“變量”,計算工程造價。

    (1)現(xiàn)澆混凝土梯形渠道。材料費用主要由渠身材料費用和模板費用構成,根據(jù)圖1計算其材料費用:

    (1)

    式中:b為渠道底寬,m;h為渠道高度,m;m為渠道邊坡系數(shù);δ為渠身混凝土厚度,m;l為渠道長度,m;E1為混凝土單價,元/m2;E2為混凝土模板單價,元/m2。

    (2)現(xiàn)澆混凝土U形渠道。材料費用主要由渠身材料費用、模板費用和壓頂費用構成。由圖1計算材料費用:

    (2)

    式中:r為圓弧段的半徑,m;α為直線段的傾斜角,°;h2為圓弧段以上水深,m;h3為安全超高,一般取0.1~0.2 m,本文計算取0.1 m;E3為 U形渠道壓頂單價,元/m3。

    (3)PE管道。材料費用主要由管道材料費用構成,即:

    F0=E4l

    (3)

    式中:E4為直徑為D的管道單價,元/m。

    唐亮等[10]發(fā)現(xiàn)管道單價與管道直徑呈冪函數(shù)關系,即:

    E4=βDγ

    (4)

    式中:β、γ為擬合系數(shù);D為管道直徑,m。

    1.1.2 人工費用和機械費

    人工費和機械費等費用按材料費用的比例計算,即:

    F1=τF0

    (5)

    式中:τ為比例系數(shù)。

    1.1.3 年運行成本

    年運行成本主要包括能耗費、工程維修費、水費等費用,除了能耗費和水費外,其他的費用可按材料費用的比例計算[11]。運行成本可表示為:

    G=σF0+C+O

    (6)

    式中:σ為運行成本系數(shù);C為能耗費,萬元;O為水資源費,萬元。

    (7)

    O=E6m毛S

    (8)

    式中:E5為電費單價,元/kWh;Q泵為泵站流量,m3;T為水泵年工作時間,h;H泵為泵站揚程,m;η裝為泵站裝機效率;E6為水價,元/m3;m毛為毛灌水定額,m3/hm2;S為灌溉面積,hm2。

    1.1.4 年需投資費用

    F=α(1+τ)F0+G

    (9)

    (10)

    式中:α為均付因子;e為資金年利率,取7%;t為經(jīng)濟計算期,取20年;F為年需投資費用,萬元。

    本文采用式(9)作為年需投資費用優(yōu)化模型目標函數(shù)。

    1.2 灌溉效益

    灌溉渠(管)道效益包括新增耕地效益與作物增產(chǎn)效益。即:

    (11)

    作物產(chǎn)量的增加不僅與灌溉輸水方式有關,而且受到作物品質(zhì),管護方式、市場行情等方面影響,灌溉產(chǎn)生的效益一般采用分攤系數(shù)法計算。王靜等研究發(fā)現(xiàn)江蘇省平均灌溉效益分攤系數(shù)為0.226[12]。

    相關研究指出[13,14],采用防滲渠道和管道地下進行輸水灌溉,既能夠有效節(jié)省土地,新增耕地,又可提高灌溉保證率,增加作物產(chǎn)量。渠道輸水灌溉平均增產(chǎn)約5%,管道輸水灌溉平均增產(chǎn)約10%。對于相同的農(nóng)田,其作物種植面積、價格與原產(chǎn)量均相同,不同的農(nóng)田灌溉設施通過改變新增耕地面積與作物增產(chǎn)產(chǎn)量來影響年效益。在已知灌區(qū)現(xiàn)狀土渠的斷面數(shù)據(jù)和產(chǎn)量的情況下,可計算不同農(nóng)田灌溉設施產(chǎn)生的新增耕地面積與作物增產(chǎn)產(chǎn)量?;诖?,本文的灌區(qū)效益(新增耕地效益與作物增產(chǎn)效益)可按式(12)與式(13)計算。

    (1)作物增產(chǎn)產(chǎn)量:

    (12)

    式中:k為作物增產(chǎn)率,%。梯形渠道和U形渠道取5%[13],PE管道取10%[13]。

    (2)新增耕地面積:

    A′=l(b0-b)

    (13)

    式中:b0為采用土渠時渠道上口寬,m;b為梯形渠道和U形渠道的上口寬,m。PE管道取0 m。

    W=B-α(1+τ)F0-G

    (14)

    式中:W為灌區(qū)年凈效益,元。

    本文采用式(14)作為年凈效益優(yōu)化模型目標函數(shù)。

    1.3 約束條件

    上述的優(yōu)化模型僅考慮經(jīng)濟性要求,雖然滿足年需投資費用最低或年凈效益最大的目標,但未考慮渠(管)道的輸水能力和流速限制,也未考慮地形高差的影響,導致部分最優(yōu)解不滿足水力約束條件。因此,必須確定農(nóng)田灌溉渠(管)道的水力約束條件,才能確保優(yōu)化模型得到的結果符合實際。

    1.3.1 輸水能力約束

    渠(管)道計算流量應等于或接近渠道設計流量,誤差不超過5%[15],即:

    (15)

    (16)

    式中:C為謝才系數(shù);A為過流斷面面積,m2;i為渠道比降;R為水力半徑,m;Q為計算流量,m3/s。

    1.3.2 流速約束

    渠(管)道中流速采用謝才公式計算:

    (17)

    無鋼筋的混凝土渠道流速不宜大于2.5 m/s[15],管道設計流速宜控制在0.9~1.5 m/s。

    1.3.3 水位約束

    (1)渠道水位約束。渠道入口處應有足夠的水頭,從而滿足控制點地面高程的要求。

    H≥H0+li+∑hj+h′

    (18)

    (19)

    式中:H為渠道進水口處水頭,m;H0為灌區(qū)控制點地面高程,m;∑hj為通過渠系建筑物時的局部水頭損失,m;h′為控制點地面與附近末級固定渠道設計水位的高差,一般取0.1~0.2 m,本文取0.1 m;ξ為局部阻力系數(shù),本文取0.5[16]。

    (2)管道水位約束。管道入口處應有足夠的水頭,滿足控制點灌水器工作水頭要求。

    H≥H0+∑hf+∑hj+HΔ

    (20)

    (21)

    式中:∑hf為管道沿程水頭損失;f為摩阻系數(shù),取9.48×105;η為流量指數(shù),取1.77;λ為管徑指數(shù),取4.77;HΔ為灌水器水頭。

    2 遺傳算法

    2.1 基本概念

    遺傳算法是模擬生物界的遺傳和進化過程而建立起來的一種并行隨機優(yōu)化算法,其對目標函數(shù)、設計變量及可行域沒有特殊要求,適用于傳統(tǒng)搜索方法解決不了的復雜和非線性問題[17]。遺傳算法需要針對具體問題,尋找合適的適應度函數(shù),確定設計變量的編碼方式,并設計相應的選擇、交叉、變異等遺傳因子。遺傳算法流程圖如圖2所示。

    圖2 遺傳算法流程圖Fig.2 The calculation process of genetic algorithm

    2.2 遺傳算法應用設計

    本文采用Python編寫基于二進制編碼的遺傳算法優(yōu)化程序,其計算過程簡述如下:

    (1)變量編碼。梯形渠道、U形渠道和PE管道設計變量的選擇如表1所示??紤]到計算機計算的方便性,每個變量在其域值內(nèi)生成一個隨機二進制數(shù)構成一個染色體,多個帶有染色體的個體組合為一個種群。

    表1 設計變量選擇表Tab.1 The selection of design variable

    (2)產(chǎn)生初始種群。隨機產(chǎn)生N個個體形成的初始種群,由這N個初始種群開始進行進化計算。

    (3)計算適應度及評價。若采用年需投資費用作為目標函數(shù),其值越小,則個體的適應度越強;若采用年凈效益作為目標函數(shù),其值越大,則個體的適應度越強。

    (4)選擇操作。完成適應度計算后,采用賭輪盤選擇方法從初始種群選擇部分個體生存并保留到下一代,適應度越大被選擇的概率越高。

    (5)交叉操作。按交叉概率PC選擇PCN/2個個體作為父輩,隨機交換某一位置的變量,產(chǎn)生新的個體。

    (6)變異操作。按變異概率Pm選擇PmN個個體,重新生成某一位置的變量,產(chǎn)生新的個體。

    2.3 遺傳算法參數(shù)選擇

    遺傳算法的控制參數(shù)有初始種群數(shù)N,交叉概率PC和變異概率Pm,不同的參數(shù)組合對于遺傳算法的運行影響較大[18]。DE JONG等[19]系統(tǒng)研究不同參數(shù)組合對于遺傳算法運行影響,并提出一組標準參數(shù),即:N=50,PC=0.60,Pm=0.001。本文遺傳算法的參數(shù)選用DE JONG提出的標準參數(shù)。

    2.4 約束條件處理

    本文選用PYTHON算法中NUMPY函數(shù)對約束條件進行處理。若不滿足約束條件,將目標函數(shù)值標記為不合法,輸出結果時將不合法的目標函數(shù)值舍去,從而達到約束的目的。

    3 算例分析

    為驗證本文方法的實際應用效果,以江蘇省臨海農(nóng)場渠道規(guī)劃整治工程為例,采用本文建立的優(yōu)化模型進行渠(管)道設計。該工程規(guī)劃圖如圖3所示,該工程控制灌溉面積為2 km2,分為兩個輪灌組進行灌溉。工程采用泵站提水灌溉,泵站設計水位為1.6 m,泵站揚程為3.17 m。該工程設計資料如表2所示。

    圖3 臨海農(nóng)場項目規(guī)劃示意圖Fig.3 The project planning diagram of Linhai Farm

    表2 工程設計資料Tab.2 The data of project design

    將上述設計資料代入優(yōu)化模型,分別計算以年需投資費用和年凈效益為目標函數(shù)的優(yōu)化模型。以U形斗渠與農(nóng)管組合為例,分析優(yōu)化模型的計算過程與計算結果。模型計算過程如圖4所示,計算結果如表3所示。

    從圖4可以看出,由于賭輪盤選擇方式的使用,兩種優(yōu)化模型種群個體目標最優(yōu)值和種群平均值變化明顯,算法收斂速

    圖4 模型優(yōu)化過程示意圖Fig 4 The calculation process of model optimization

    表3 最優(yōu)結果典型值Tab.3 Typical value of the optimal result

    度較快[20],分別在90次和76次遺傳代數(shù)時得到最優(yōu)解。在計算機上運行時間為1.30 s,計算時間短[21]。

    由表3可知,當U形斗渠傾斜角取18.12°,農(nóng)管管徑取280 mm時,該工程的年需投資費用最小,為39.36 萬元;而當U形斗渠傾斜角取6.41°,農(nóng)管管徑取280 mm時,年凈效益最大,為6.29 萬元。此外,兩個模型的校核流量均小于10-15m3/m3,在U形渠道中流速為0.70 m/s,在管道中流速為1.08 m/s,兩種優(yōu)化模型結果均滿足其設計的約束條件。由上述結果可知,遺傳算法能夠在較短的時間內(nèi)得到較好的結果,將遺傳算法應用于灌溉渠(管)道設計是可行的。

    4 田間灌溉渠系組合方式研究

    由上述結果可知,遺傳算法能夠應用于灌溉輸水渠(管)道設計,下面將探究田間灌溉渠系最優(yōu)組合方式?;谀陜粜б孀畲髢?yōu)化模型,采用梯形渠道、U形渠道和管道三種不同方式的斗渠與農(nóng)渠組合,探究年凈效益最大的田間灌溉渠系組合方式,即為最優(yōu)組合方式。模型計算結果如表4所示。

    由表4可知,不同形式的斗渠顯著影響年需投資費用(P<0.05),其中斗管的年需投資費用最高,平均為56.95 萬元,遠遠高于梯形斗渠的33.90 萬元和U形斗渠的34.05 萬元。相應的,年凈效益也隨著斗渠形式的變化而顯著改變,其中梯形斗渠年凈效益和U形斗渠年凈效益遠高于斗管的年凈效益。農(nóng)渠的不同形式顯著影響襯砌渠道灌溉和管道灌溉新增的耕地面積。根據(jù)表4,采用農(nóng)管新增耕地面積最大,平均為8.79萬m2,而采用梯形農(nóng)渠和U形農(nóng)渠的新增耕地僅為6.26和7.69 萬m2。

    在9種農(nóng)田灌溉渠系組合形式中,梯形斗渠和農(nóng)管的組合與U形斗渠和農(nóng)管的組合年凈效益高于其他類型的組合方式。梯形斗渠與農(nóng)管的組合年凈效益最大為6.54 萬元,其新增耕地為8.64 萬m2;采用U形斗渠和農(nóng)管的組合較梯形斗渠和農(nóng)管組合相比,年凈效益降低3.97%,新增耕地增加1.37%;而斗管與農(nóng)管的組合方式雖然新增耕地最高為8.96 萬m2,但其年需投資費用最高為72.84 萬元,年凈效益最低為-20.71 萬元,故在實際工程中不考慮采用這種組合形式。綜上所示,采用梯形斗渠與農(nóng)管和U形斗渠與農(nóng)管的組合方式較為理想。

    表4 綜合年需投資費用最小灌溉系統(tǒng)優(yōu)化結果表Tab.4 Optimization results of minimum annual investment cost for irrigation system

    5 討 論

    本文結果表明,基于遺傳算法的優(yōu)化模型能夠在較短的時間內(nèi)處理農(nóng)田灌溉輸水渠(管)道組合設計非線性優(yōu)化問題,且計算結果較為穩(wěn)定[22]。本文構建的優(yōu)化模型在100次遺傳代數(shù)左右得到最優(yōu)解,模型收斂速度快,與李云峰等[20]建立的以供水成本最低的目標函數(shù)的渠道優(yōu)化模型相似,其建立的優(yōu)化模型在50次左右遺傳代數(shù)時得到最優(yōu)解。劉波等[21]基于遺傳算法建立的農(nóng)田有壓管道系統(tǒng)整體優(yōu)化方案計算時間為6.250 s,而本模型計算時間為1.30 s,計算速度較快。由表3可知,模型能夠在得到最優(yōu)解的同時,較好地處理約束條件。與MATLAB工具箱中的遺傳算法相比[20],本文采用的PYTHON遺傳算法處理約束條件不需要構造懲罰函數(shù),而是將不符合約束條件的解標記為不合法,輸出時去除不合法的解,從而得到更準確的優(yōu)化結果。從圖4中可以看出,兩種不同的模型均在100次遺傳代數(shù)左右得到最優(yōu)解,在后期停滯不前。這一方面可能是因為模型中的變量較少,模型收斂速度快;另一方面可能是因為本文采用的是DE JONG提出的遺傳算法參數(shù),參數(shù)不會隨著遺傳過程的改變而自動調(diào)整,從而發(fā)生了早熟收斂。李康順等[23]提出將改進遺傳算法根據(jù)群體的分散程度和個體的適應度而自動控制算法參數(shù),克服傳統(tǒng)遺傳算法易早熟的問題同時保持群體的多樣性,加快收斂速度。因此,后期可采用李康順等提出的改進遺傳算法優(yōu)化灌溉渠(管)道優(yōu)化設計與選型。

    此外,本文研究的是以年凈效益為目標函數(shù)的田間灌溉輸水渠(管)道設計與組合方式。蔣曉紅等[6]以渠道設計水深作為優(yōu)化變量,將均勻流公式計算設計流量轉化為非線性優(yōu)化問題的目標函數(shù)進行渠道設計。李云峰等[20]以灌溉用水成本為目標函數(shù),建立設計梯形渠道和U形渠道的非線性優(yōu)化模型。本文的優(yōu)化模型綜合考慮了渠(管)道設計的水力性能和經(jīng)濟效益,并確定兩級渠(管)道最優(yōu)的組合形式,能夠更好地指導實際工程的設計。但是灌溉水資源成本和灌溉效益選用的均為當?shù)囟嗄昶骄鶖?shù)據(jù),并未考慮作物實際的灌溉制度。GONZLEZ等[24]根據(jù)土壤水分平衡方程和遺傳算法建立灌溉水調(diào)度優(yōu)化模型,采用該模型優(yōu)化作物的灌溉制度后能夠節(jié)約15%的灌溉水資源成本。劉向等[25]通過對新疆棗棉田研究,結合作物水分生產(chǎn)函數(shù)模型和水量平衡理論,建立棗棉田灌溉——收益模型,得到在有限灌水量條件下的單位面積最大經(jīng)濟收益及棗棉各生長階段的最優(yōu)灌水量分配方案。綜上,下一步研究可將作物實際的灌溉制度與本文研究的灌溉渠(管)道設計組合成果相結合,采用改進遺傳算法,優(yōu)化本文構建的模型。

    6 結 論

    (1)基于遺傳算法模型優(yōu)化灌溉渠(管)道設計,均在100次遺傳代數(shù)內(nèi)得到最優(yōu)結果,計算時間為1.30 s,速度快,優(yōu)化結果好。遺傳算法可應用于灌溉渠(管)道的設計。

    (2)采用梯形斗渠和U形斗渠能顯著減少約40%年需投資費用,確保正年凈效益(P<0.05);采用農(nóng)管能顯著增加10%以上新增耕地。

    (3)采用梯形斗渠與農(nóng)管組合和U形斗渠與農(nóng)管組合年凈效益較大,均超過6 萬元/a,較其他組合新增耕地增加10%以上。實際工程中可采用此方式作為田間灌溉渠系最優(yōu)組合方式。

    猜你喜歡
    斗渠遺傳算法效益
    草粉發(fā)酵 喂羊效益高
    蓮魚混養(yǎng) 效益提高一倍
    基于灌溉渠水質(zhì)保護的都江堰灌區(qū)斗渠廊道LID生態(tài)措施研究
    冬棚養(yǎng)蝦效益顯著,看技術達人如何手到“錢”來
    果園有了“鵝幫工” 一舉多得效益好
    基于自適應遺傳算法的CSAMT一維反演
    一種基于遺傳算法的聚類分析方法在DNA序列比較中的應用
    基于遺傳算法和LS-SVM的財務危機預測
    某縣新增千億斤糧食產(chǎn)能規(guī)劃渠道工程典型設計
    治淮(2016年7期)2016-09-02 02:55:41
    基于改進的遺傳算法的模糊聚類算法
    啪啪无遮挡十八禁网站| 十八禁网站免费在线| 99re在线观看精品视频| 欧美日本中文国产一区发布| 欧美亚洲 丝袜 人妻 在线| 亚洲成a人片在线一区二区| 黄色a级毛片大全视频| 黄色成人免费大全| 免费少妇av软件| 黑丝袜美女国产一区| 久久香蕉激情| 免费观看精品视频网站| 999久久久精品免费观看国产| 精品一区二区三卡| 亚洲一码二码三码区别大吗| av网站在线播放免费| 国产成+人综合+亚洲专区| 精品电影一区二区在线| 五月开心婷婷网| 午夜精品在线福利| 日日夜夜操网爽| 精品国产超薄肉色丝袜足j| 久热这里只有精品99| 满18在线观看网站| 免费久久久久久久精品成人欧美视频| 欧美 亚洲 国产 日韩一| 婷婷成人精品国产| 午夜福利一区二区在线看| 久久中文字幕一级| 亚洲免费av在线视频| 国产精品电影一区二区三区 | 精品福利永久在线观看| 亚洲成国产人片在线观看| 色综合婷婷激情| 中文字幕制服av| 欧美日本中文国产一区发布| 久久久精品国产亚洲av高清涩受| 欧美日韩国产mv在线观看视频| 亚洲 国产 在线| 国产欧美日韩精品亚洲av| 又大又爽又粗| 91老司机精品| 999久久久精品免费观看国产| 午夜福利,免费看| 久久久久国内视频| 日韩欧美在线二视频 | 18在线观看网站| 成熟少妇高潮喷水视频| 变态另类成人亚洲欧美熟女 | 午夜精品久久久久久毛片777| 下体分泌物呈黄色| 美女福利国产在线| 欧美不卡视频在线免费观看 | 欧美乱码精品一区二区三区| 国产视频一区二区在线看| 三上悠亚av全集在线观看| 亚洲片人在线观看| 国产三级黄色录像| 欧美在线一区亚洲| 日韩有码中文字幕| 欧洲精品卡2卡3卡4卡5卡区| 窝窝影院91人妻| 9热在线视频观看99| 嫩草影视91久久| 一a级毛片在线观看| 精品一区二区三卡| 国产亚洲精品第一综合不卡| 国产高清视频在线播放一区| 又紧又爽又黄一区二区| a级毛片黄视频| 美女视频免费永久观看网站| tube8黄色片| 欧美黑人欧美精品刺激| 亚洲九九香蕉| 国产精品偷伦视频观看了| 国产97色在线日韩免费| 亚洲情色 制服丝袜| 99久久99久久久精品蜜桃| 91精品国产国语对白视频| 国产成人免费无遮挡视频| 人妻丰满熟妇av一区二区三区 | 一级毛片女人18水好多| 身体一侧抽搐| 99精国产麻豆久久婷婷| 成人国语在线视频| 亚洲va日本ⅴa欧美va伊人久久| 天天影视国产精品| 国产成人精品久久二区二区91| 大码成人一级视频| 色婷婷久久久亚洲欧美| 久久精品熟女亚洲av麻豆精品| 天天操日日干夜夜撸| svipshipincom国产片| 女同久久另类99精品国产91| 日本黄色日本黄色录像| 如日韩欧美国产精品一区二区三区| 看免费av毛片| 极品少妇高潮喷水抽搐| 日韩欧美免费精品| e午夜精品久久久久久久| 欧美乱色亚洲激情| 欧美av亚洲av综合av国产av| 狠狠婷婷综合久久久久久88av| 久久亚洲真实| 成人国产一区最新在线观看| 老熟女久久久| 亚洲av成人不卡在线观看播放网| 亚洲精品中文字幕一二三四区| 99热只有精品国产| 在线永久观看黄色视频| 日韩视频一区二区在线观看| 69av精品久久久久久| 老汉色av国产亚洲站长工具| 久久久国产精品麻豆| 免费看十八禁软件| 国产区一区二久久| 91老司机精品| 久久久久国产精品人妻aⅴ院 | 捣出白浆h1v1| 狠狠狠狠99中文字幕| av福利片在线| 日韩欧美在线二视频 | 国产男女超爽视频在线观看| 91国产中文字幕| 国产亚洲欧美精品永久| 黄片播放在线免费| 999久久久精品免费观看国产| 日韩三级视频一区二区三区| 在线av久久热| 亚洲三区欧美一区| 午夜福利,免费看| 久久午夜亚洲精品久久| 免费观看精品视频网站| 成人av一区二区三区在线看| 少妇裸体淫交视频免费看高清 | 亚洲avbb在线观看| 国产精品久久视频播放| 国产麻豆69| 精品久久久久久久久久免费视频 | 操出白浆在线播放| 69精品国产乱码久久久| 国产激情欧美一区二区| 国产精品乱码一区二三区的特点 | 久久午夜综合久久蜜桃| 久久 成人 亚洲| 国产成人欧美| 99精品欧美一区二区三区四区| 电影成人av| 国产激情久久老熟女| 男女高潮啪啪啪动态图| 无限看片的www在线观看| 久久国产精品影院| 午夜免费成人在线视频| 人妻 亚洲 视频| 国产一区二区三区在线臀色熟女 | 久久午夜亚洲精品久久| 久久久久视频综合| 精品少妇久久久久久888优播| 超碰97精品在线观看| 99精品欧美一区二区三区四区| 国产av又大| 午夜免费鲁丝| 丝瓜视频免费看黄片| 久久中文字幕一级| 色综合欧美亚洲国产小说| 9热在线视频观看99| 亚洲片人在线观看| 免费看十八禁软件| 9色porny在线观看| 少妇粗大呻吟视频| 天堂√8在线中文| 亚洲国产欧美一区二区综合| 亚洲熟女毛片儿| 欧美 亚洲 国产 日韩一| 国产一区二区三区视频了| videos熟女内射| 新久久久久国产一级毛片| 老汉色∧v一级毛片| 亚洲精品中文字幕在线视频| 欧美 日韩 精品 国产| 成人av一区二区三区在线看| 99久久人妻综合| 亚洲精品久久成人aⅴ小说| 午夜福利一区二区在线看| 国产成人精品在线电影| 国产精品九九99| 大香蕉久久网| 丝袜人妻中文字幕| 人人妻人人澡人人看| 国产亚洲欧美精品永久| 国产免费男女视频| 九色亚洲精品在线播放| 国产av一区二区精品久久| a在线观看视频网站| 人人澡人人妻人| 国产亚洲一区二区精品| 欧美老熟妇乱子伦牲交| 制服诱惑二区| 亚洲专区字幕在线| 在线观看一区二区三区激情| 亚洲五月天丁香| 亚洲欧洲精品一区二区精品久久久| 侵犯人妻中文字幕一二三四区| 91麻豆精品激情在线观看国产 | 国产亚洲精品久久久久5区| 亚洲精品乱久久久久久| a级毛片黄视频| 在线永久观看黄色视频| 中文欧美无线码| 国产激情欧美一区二区| 人人妻,人人澡人人爽秒播| 亚洲人成77777在线视频| 亚洲成a人片在线一区二区| 国产乱人伦免费视频| 亚洲视频免费观看视频| 在线观看免费午夜福利视频| 亚洲一区中文字幕在线| 国产精品免费一区二区三区在线 | 激情视频va一区二区三区| 亚洲人成77777在线视频| 国产精品欧美亚洲77777| 精品电影一区二区在线| 亚洲aⅴ乱码一区二区在线播放 | 韩国精品一区二区三区| 亚洲七黄色美女视频| 久久人人爽av亚洲精品天堂| 国产在线观看jvid| 日韩大码丰满熟妇| 久久 成人 亚洲| 日本欧美视频一区| 桃红色精品国产亚洲av| 精品卡一卡二卡四卡免费| 搡老岳熟女国产| 老鸭窝网址在线观看| 亚洲在线自拍视频| 91成人精品电影| 日韩大码丰满熟妇| 国产精品亚洲av一区麻豆| 又黄又爽又免费观看的视频| 国产男女超爽视频在线观看| 亚洲avbb在线观看| 国产精品一区二区精品视频观看| 久久国产亚洲av麻豆专区| 国产成人啪精品午夜网站| 欧美人与性动交α欧美软件| 亚洲全国av大片| 亚洲av成人av| 麻豆国产av国片精品| 99久久综合精品五月天人人| 男男h啪啪无遮挡| 亚洲va日本ⅴa欧美va伊人久久| 久久国产精品影院| 国产成人系列免费观看| 午夜福利欧美成人| 成年人免费黄色播放视频| 69精品国产乱码久久久| 欧美精品av麻豆av| 又紧又爽又黄一区二区| 亚洲,欧美精品.| 制服诱惑二区| 国产99久久九九免费精品| 欧美亚洲日本最大视频资源| 国产一区二区激情短视频| 无遮挡黄片免费观看| 日韩大码丰满熟妇| 欧美乱色亚洲激情| 美女高潮到喷水免费观看| 成人18禁在线播放| 国产亚洲精品第一综合不卡| 女警被强在线播放| 国产亚洲精品久久久久久毛片 | 国精品久久久久久国模美| 日韩欧美一区视频在线观看| 国产区一区二久久| 老汉色∧v一级毛片| 大型av网站在线播放| 在线观看免费高清a一片| 丰满的人妻完整版| 国产亚洲欧美98| 精品国产超薄肉色丝袜足j| 日本a在线网址| 色播在线永久视频| 人妻一区二区av| 99久久综合精品五月天人人| 午夜免费成人在线视频| 夜夜爽天天搞| 一本大道久久a久久精品| 99国产精品免费福利视频| 老司机午夜福利在线观看视频| 免费在线观看影片大全网站| 欧美亚洲 丝袜 人妻 在线| 咕卡用的链子| 国产精品一区二区在线观看99| 国内久久婷婷六月综合欲色啪| 成人国产一区最新在线观看| 99香蕉大伊视频| 国产精品久久久久成人av| 美女高潮喷水抽搐中文字幕| 免费观看a级毛片全部| 国产高清国产精品国产三级| 国产精品国产高清国产av | 欧美色视频一区免费| www.熟女人妻精品国产| 久热爱精品视频在线9| 免费黄频网站在线观看国产| 亚洲av日韩精品久久久久久密| 免费在线观看黄色视频的| 亚洲精品av麻豆狂野| 国产深夜福利视频在线观看| 黄色视频,在线免费观看| 国产97色在线日韩免费| 国产精品99久久99久久久不卡| 亚洲av美国av| 欧美成人免费av一区二区三区 | 久久性视频一级片| 亚洲一区中文字幕在线| 成人亚洲精品一区在线观看| 国产成人av激情在线播放| av视频免费观看在线观看| 国产97色在线日韩免费| 久久人人97超碰香蕉20202| 超色免费av| 不卡av一区二区三区| 免费高清在线观看日韩| 窝窝影院91人妻| 宅男免费午夜| av一本久久久久| tube8黄色片| 国产极品粉嫩免费观看在线| 中文字幕制服av| 国产色视频综合| 又紧又爽又黄一区二区| 免费观看人在逋| 窝窝影院91人妻| 一本一本久久a久久精品综合妖精| 午夜成年电影在线免费观看| 99re在线观看精品视频| 在线av久久热| 亚洲成人免费电影在线观看| 国产淫语在线视频| 一级a爱视频在线免费观看| 亚洲精品在线观看二区| 99精国产麻豆久久婷婷| 操美女的视频在线观看| 亚洲色图av天堂| 午夜福利,免费看| 午夜老司机福利片| 亚洲va日本ⅴa欧美va伊人久久| 热99久久久久精品小说推荐| 国产精华一区二区三区| 99国产精品99久久久久| 啦啦啦免费观看视频1| 青草久久国产| 美女福利国产在线| av在线播放免费不卡| 成年人黄色毛片网站| 在线观看舔阴道视频| 两个人免费观看高清视频| 高清毛片免费观看视频网站 | 欧美日韩国产mv在线观看视频| 久久久久久久午夜电影 | 午夜视频精品福利| 亚洲熟女精品中文字幕| 麻豆乱淫一区二区| 日韩大码丰满熟妇| 一级黄色大片毛片| 大型av网站在线播放| 亚洲va日本ⅴa欧美va伊人久久| 国产精品一区二区在线观看99| 怎么达到女性高潮| 在线免费观看的www视频| 精品无人区乱码1区二区| 色老头精品视频在线观看| 91成年电影在线观看| √禁漫天堂资源中文www| 久久久久久免费高清国产稀缺| 狠狠狠狠99中文字幕| 美国免费a级毛片| 色94色欧美一区二区| 成人精品一区二区免费| 色94色欧美一区二区| 亚洲av片天天在线观看| 国产成人精品在线电影| 天天躁日日躁夜夜躁夜夜| 亚洲av日韩精品久久久久久密| 中文字幕高清在线视频| videosex国产| 女人爽到高潮嗷嗷叫在线视频| 香蕉丝袜av| 久久精品亚洲精品国产色婷小说| 18禁美女被吸乳视频| 国产精品久久久久成人av| 狠狠狠狠99中文字幕| 欧美日韩瑟瑟在线播放| 亚洲欧美一区二区三区黑人| 亚洲精品一卡2卡三卡4卡5卡| 一级片免费观看大全| 国产精品一区二区精品视频观看| 18在线观看网站| 激情视频va一区二区三区| 成人手机av| 成人特级黄色片久久久久久久| 十八禁人妻一区二区| 精品一区二区三区视频在线观看免费 | 老熟女久久久| 老司机午夜十八禁免费视频| 久久国产乱子伦精品免费另类| 国产成人欧美在线观看 | 黄色怎么调成土黄色| 亚洲av成人av| netflix在线观看网站| 老司机影院毛片| 好看av亚洲va欧美ⅴa在| 国产精品av久久久久免费| 韩国av一区二区三区四区| 久久精品人人爽人人爽视色| 亚洲一区二区三区不卡视频| 大型av网站在线播放| 777久久人妻少妇嫩草av网站| 丝袜人妻中文字幕| 超碰97精品在线观看| 午夜精品久久久久久毛片777| 日韩免费av在线播放| 美女 人体艺术 gogo| 亚洲少妇的诱惑av| 欧美黄色淫秽网站| 捣出白浆h1v1| 亚洲av片天天在线观看| 天堂中文最新版在线下载| 色老头精品视频在线观看| 电影成人av| 国产亚洲av高清不卡| 久久精品亚洲av国产电影网| 999久久久国产精品视频| 夜夜夜夜夜久久久久| 搡老岳熟女国产| 美女扒开内裤让男人捅视频| 人妻 亚洲 视频| 搡老乐熟女国产| 精品乱码久久久久久99久播| 啦啦啦视频在线资源免费观看| 一个人免费在线观看的高清视频| 亚洲va日本ⅴa欧美va伊人久久| 亚洲aⅴ乱码一区二区在线播放 | 一本一本久久a久久精品综合妖精| 精品少妇久久久久久888优播| 人人妻人人澡人人爽人人夜夜| 黑人操中国人逼视频| 国产精品自产拍在线观看55亚洲 | 久久久久久久久久久久大奶| 侵犯人妻中文字幕一二三四区| 两性夫妻黄色片| 久久婷婷成人综合色麻豆| 一级作爱视频免费观看| 亚洲精品在线观看二区| 欧美成人免费av一区二区三区 | 久久中文字幕一级| 成人永久免费在线观看视频| 亚洲五月天丁香| 不卡一级毛片| 亚洲欧美激情综合另类| 一a级毛片在线观看| 久久人妻av系列| 欧美激情极品国产一区二区三区| 在线观看舔阴道视频| 久久99一区二区三区| 亚洲熟女毛片儿| 国产熟女午夜一区二区三区| 国产野战对白在线观看| 多毛熟女@视频| 亚洲成av片中文字幕在线观看| 99国产精品99久久久久| 亚洲七黄色美女视频| 午夜两性在线视频| 99热国产这里只有精品6| videos熟女内射| 亚洲精品国产一区二区精华液| 国产一卡二卡三卡精品| 国产精品电影一区二区三区 | 少妇 在线观看| 男人操女人黄网站| 国产野战对白在线观看| 在线播放国产精品三级| 国产一区二区三区综合在线观看| av电影中文网址| 电影成人av| 亚洲精品中文字幕在线视频| 精品久久久久久电影网| 久久精品熟女亚洲av麻豆精品| 久久热在线av| 亚洲性夜色夜夜综合| 国产精品久久久久成人av| 久久亚洲真实| 欧美日韩亚洲高清精品| 国产成人欧美在线观看 | 啦啦啦视频在线资源免费观看| 成在线人永久免费视频| 欧美不卡视频在线免费观看 | 久久久久久久久免费视频了| 国产免费现黄频在线看| 又黄又爽又免费观看的视频| 一区在线观看完整版| 国产精品影院久久| 国产精品 欧美亚洲| 亚洲成av片中文字幕在线观看| 9191精品国产免费久久| 免费人成视频x8x8入口观看| 国产亚洲精品久久久久久毛片 | 国产一区有黄有色的免费视频| 成人18禁在线播放| 99国产精品免费福利视频| 精品久久久精品久久久| 亚洲国产欧美一区二区综合| 亚洲 国产 在线| 一级片'在线观看视频| 美女高潮喷水抽搐中文字幕| 在线国产一区二区在线| 日本a在线网址| 国产精品亚洲一级av第二区| 18禁黄网站禁片午夜丰满| videos熟女内射| 夫妻午夜视频| 亚洲熟妇中文字幕五十中出 | videos熟女内射| 在线播放国产精品三级| 成人亚洲精品一区在线观看| 欧美精品亚洲一区二区| 女警被强在线播放| 美女午夜性视频免费| 久久热在线av| 久久人妻福利社区极品人妻图片| 午夜两性在线视频| 欧美黄色淫秽网站| 亚洲,欧美精品.| 国产成人精品久久二区二区免费| 精品免费久久久久久久清纯 | 热99re8久久精品国产| 国产一区在线观看成人免费| 夜夜躁狠狠躁天天躁| 欧美日韩乱码在线| 妹子高潮喷水视频| 大码成人一级视频| 深夜精品福利| 亚洲av成人一区二区三| 最近最新免费中文字幕在线| ponron亚洲| 婷婷丁香在线五月| 亚洲久久久国产精品| 亚洲中文日韩欧美视频| 99国产精品一区二区蜜桃av | 18禁裸乳无遮挡动漫免费视频| 国产91精品成人一区二区三区| 精品国产国语对白av| 中亚洲国语对白在线视频| 嫩草影视91久久| 色尼玛亚洲综合影院| 黄色视频不卡| 老熟妇乱子伦视频在线观看| 香蕉国产在线看| 91成人精品电影| 欧美大码av| 人人妻人人澡人人爽人人夜夜| 黄色成人免费大全| 免费在线观看亚洲国产| 国产一卡二卡三卡精品| 成人永久免费在线观看视频| 黄片大片在线免费观看| 国产野战对白在线观看| 欧美大码av| 大陆偷拍与自拍| 黄色成人免费大全| 日韩有码中文字幕| 久久精品国产亚洲av香蕉五月 | 久久久久久久精品吃奶| 每晚都被弄得嗷嗷叫到高潮| 欧美成人免费av一区二区三区 | 最新美女视频免费是黄的| 叶爱在线成人免费视频播放| 日本黄色视频三级网站网址 | 美女高潮喷水抽搐中文字幕| 欧美 亚洲 国产 日韩一| 一区在线观看完整版| 国产淫语在线视频| 亚洲国产欧美网| 国产黄色免费在线视频| 女人高潮潮喷娇喘18禁视频| 久久午夜综合久久蜜桃| 亚洲欧美一区二区三区黑人| 中文字幕制服av| 正在播放国产对白刺激| 日韩欧美免费精品| 免费在线观看影片大全网站| 黑人猛操日本美女一级片| 天天添夜夜摸| 国产亚洲一区二区精品| 熟女少妇亚洲综合色aaa.| 国产免费现黄频在线看| 变态另类成人亚洲欧美熟女 | 人人妻人人澡人人看| 在线观看免费视频网站a站| 在线国产一区二区在线| 久久亚洲真实| 日本五十路高清| 我的亚洲天堂| av网站免费在线观看视频| 黑人巨大精品欧美一区二区蜜桃| 久久人人爽av亚洲精品天堂| 久久99一区二区三区| 韩国精品一区二区三区| 国产成人av激情在线播放| 交换朋友夫妻互换小说| videos熟女内射| 亚洲 欧美一区二区三区| 亚洲黑人精品在线| 精品亚洲成国产av|