• <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
    基于改進的遺傳算法的模糊聚類算法
    97人妻天天添夜夜摸| 男人舔女人的私密视频| 80岁老熟妇乱子伦牲交| 纵有疾风起免费观看全集完整版| 亚洲国产av影院在线观看| 人人澡人人妻人| 国产精品1区2区在线观看. | 日日夜夜操网爽| av一本久久久久| 啦啦啦在线免费观看视频4| netflix在线观看网站| 99久久99久久久精品蜜桃| 午夜免费鲁丝| 丝袜脚勾引网站| 看免费av毛片| www.精华液| 一本综合久久免费| 日日爽夜夜爽网站| 我的亚洲天堂| 亚洲情色 制服丝袜| 成人三级做爰电影| 一区在线观看完整版| 欧美在线一区亚洲| 大陆偷拍与自拍| 性色av乱码一区二区三区2| 永久免费av网站大全| 午夜福利免费观看在线| 久久精品人人爽人人爽视色| 成年人免费黄色播放视频| 捣出白浆h1v1| 捣出白浆h1v1| 久久久久久久大尺度免费视频| 我要看黄色一级片免费的| 国产精品一区二区在线不卡| 人妻久久中文字幕网| 国产极品粉嫩免费观看在线| 电影成人av| 91成年电影在线观看| 美女中出高潮动态图| 老熟女久久久| 日韩中文字幕欧美一区二区| 五月天丁香电影| av国产精品久久久久影院| svipshipincom国产片| 少妇人妻久久综合中文| 桃红色精品国产亚洲av| av超薄肉色丝袜交足视频| 男女免费视频国产| 久久女婷五月综合色啪小说| 久久99一区二区三区| 成在线人永久免费视频| 国产色视频综合| 国产成人欧美在线观看 | 日本vs欧美在线观看视频| 波多野结衣一区麻豆| 少妇的丰满在线观看| 精品第一国产精品| 制服诱惑二区| 久久久国产一区二区| 啦啦啦在线免费观看视频4| 夜夜骑夜夜射夜夜干| 久久久精品区二区三区| 日韩 欧美 亚洲 中文字幕| 人妻一区二区av| 国产亚洲精品久久久久5区| 国产一区二区三区av在线| 久久久久视频综合| 9191精品国产免费久久| 一边摸一边做爽爽视频免费| 老熟妇仑乱视频hdxx| 国产真人三级小视频在线观看| 国产伦理片在线播放av一区| 国产精品 国内视频| 亚洲国产精品一区二区三区在线| 欧美黑人精品巨大| 丝袜美腿诱惑在线| 欧美成狂野欧美在线观看| 各种免费的搞黄视频| 国产成人精品无人区| 免费高清在线观看视频在线观看| 欧美另类亚洲清纯唯美| 亚洲成人手机| 丝袜美足系列| 婷婷成人精品国产| 美女中出高潮动态图| 黄片小视频在线播放| 欧美激情高清一区二区三区| 欧美人与性动交α欧美精品济南到| 婷婷丁香在线五月| 国产亚洲欧美精品永久| 久久久久久久久免费视频了| 性色av一级| 国产欧美日韩一区二区三区在线| 久久av网站| 久久精品国产亚洲av高清一级| 国产一区二区三区av在线| 日本猛色少妇xxxxx猛交久久| 少妇猛男粗大的猛烈进出视频| 一二三四在线观看免费中文在| av天堂在线播放| 人妻 亚洲 视频| 黑丝袜美女国产一区| 后天国语完整版免费观看| 久久久国产一区二区| 亚洲国产欧美网| 久久国产精品影院| 波多野结衣av一区二区av| 午夜福利,免费看| 国产伦人伦偷精品视频| 欧美日韩一级在线毛片| 亚洲精品国产一区二区精华液| 人妻久久中文字幕网| 日日夜夜操网爽| 一边摸一边做爽爽视频免费| 欧美黑人精品巨大| 人成视频在线观看免费观看| 亚洲av男天堂| 91精品国产国语对白视频| 亚洲成人免费电影在线观看| 欧美变态另类bdsm刘玥| 免费女性裸体啪啪无遮挡网站| 99久久国产精品久久久| 中文字幕高清在线视频| 超碰97精品在线观看| av欧美777| 俄罗斯特黄特色一大片| 视频区图区小说| 欧美久久黑人一区二区| 国产精品一区二区精品视频观看| 男女无遮挡免费网站观看| 日韩熟女老妇一区二区性免费视频| a在线观看视频网站| 国产激情久久老熟女| 我要看黄色一级片免费的| 久久人人97超碰香蕉20202| 亚洲精品一二三| 国产伦人伦偷精品视频| 日本vs欧美在线观看视频| 日本av手机在线免费观看| 一二三四在线观看免费中文在| 国产精品国产三级国产专区5o| 国产精品麻豆人妻色哟哟久久| 满18在线观看网站| 国产一区二区 视频在线| 亚洲美女黄色视频免费看| 午夜福利在线观看吧| 丰满迷人的少妇在线观看| 狠狠精品人妻久久久久久综合| 啦啦啦 在线观看视频| 亚洲av男天堂| 国产亚洲一区二区精品| 国内毛片毛片毛片毛片毛片| 亚洲精品国产色婷婷电影| 国产无遮挡羞羞视频在线观看| 成年人免费黄色播放视频| 欧美黄色淫秽网站| 国产日韩欧美视频二区| 91成年电影在线观看| 色综合欧美亚洲国产小说| 精品少妇一区二区三区视频日本电影| videosex国产| 法律面前人人平等表现在哪些方面 | 老熟妇仑乱视频hdxx| 亚洲九九香蕉| 欧美久久黑人一区二区| 久久久国产欧美日韩av| 国产黄频视频在线观看| 午夜福利影视在线免费观看| 亚洲第一青青草原| 一进一出抽搐动态| 精品福利观看| 欧美日韩中文字幕国产精品一区二区三区 | cao死你这个sao货| 国产1区2区3区精品| 亚洲成人免费av在线播放| 欧美在线一区亚洲| 男男h啪啪无遮挡| 搡老岳熟女国产| 亚洲av国产av综合av卡| 人妻人人澡人人爽人人| 欧美激情 高清一区二区三区| 欧美人与性动交α欧美精品济南到| 亚洲专区字幕在线| 久久中文看片网| 极品人妻少妇av视频| 欧美日韩av久久| 亚洲国产日韩一区二区| 亚洲精品国产一区二区精华液| 不卡一级毛片| 免费观看人在逋| www.熟女人妻精品国产| 午夜免费成人在线视频| 汤姆久久久久久久影院中文字幕| 99九九在线精品视频| 国产深夜福利视频在线观看| a级毛片在线看网站| 脱女人内裤的视频| 日韩精品免费视频一区二区三区| 国产深夜福利视频在线观看| 亚洲人成电影观看| 一本大道久久a久久精品| 日韩欧美一区二区三区在线观看 | 日本91视频免费播放| 国产成人影院久久av| 精品一区二区三区av网在线观看 | 97在线人人人人妻| 制服诱惑二区| 一二三四在线观看免费中文在| 久久久水蜜桃国产精品网| 欧美黄色淫秽网站| 12—13女人毛片做爰片一| 91av网站免费观看| 久久 成人 亚洲| 国产精品香港三级国产av潘金莲| 久久九九热精品免费| 欧美人与性动交α欧美软件| 日韩熟女老妇一区二区性免费视频| 国产精品免费大片| 婷婷色av中文字幕| 满18在线观看网站| 久久精品国产综合久久久| 国产日韩欧美亚洲二区| 亚洲一卡2卡3卡4卡5卡精品中文| 亚洲欧美精品综合一区二区三区| 久久久国产精品麻豆| 国产精品免费大片| 免费黄频网站在线观看国产| 久久精品久久久久久噜噜老黄| 777久久人妻少妇嫩草av网站| 女人久久www免费人成看片| 2018国产大陆天天弄谢| 嫩草影视91久久| 纵有疾风起免费观看全集完整版| 精品久久久精品久久久| 日本精品一区二区三区蜜桃| 欧美中文综合在线视频| 成年av动漫网址| 老鸭窝网址在线观看| 亚洲 欧美一区二区三区| 日韩电影二区| 国产成人免费观看mmmm| 午夜视频精品福利| 亚洲精品久久久久久婷婷小说| 精品高清国产在线一区| 宅男免费午夜| 国产男女内射视频| 国产精品影院久久| 久久久久国内视频| 岛国在线观看网站| www.av在线官网国产| 欧美另类一区| 天天影视国产精品| 在线观看免费午夜福利视频| 777久久人妻少妇嫩草av网站| 亚洲人成77777在线视频| 伦理电影免费视频| 少妇精品久久久久久久| 99久久国产精品久久久| 成人av一区二区三区在线看 | 啦啦啦啦在线视频资源| 亚洲黑人精品在线| 王馨瑶露胸无遮挡在线观看| 黄片小视频在线播放| 老司机福利观看| 国产深夜福利视频在线观看| 欧美日韩成人在线一区二区| 国产精品熟女久久久久浪| 成人免费观看视频高清| 精品第一国产精品| 另类亚洲欧美激情| 国产亚洲午夜精品一区二区久久| 欧美精品高潮呻吟av久久| 一区在线观看完整版| 国产精品免费视频内射| 亚洲欧美精品自产自拍| 亚洲av片天天在线观看| 亚洲国产日韩一区二区| 精品国产乱子伦一区二区三区 | 一区二区三区四区激情视频| 久久久精品免费免费高清| 视频在线观看一区二区三区| 成年人午夜在线观看视频| 国产精品.久久久| 成人影院久久| 另类亚洲欧美激情| 黑人猛操日本美女一级片| 男女午夜视频在线观看| 一进一出抽搐动态| 啦啦啦啦在线视频资源| 桃红色精品国产亚洲av| 一边摸一边做爽爽视频免费| 精品国产乱码久久久久久男人| 亚洲av电影在线观看一区二区三区| 午夜影院在线不卡| 国产一区二区在线观看av| 国产一级毛片在线| 欧美日韩成人在线一区二区| 丝袜脚勾引网站| 欧美激情高清一区二区三区| 天天躁日日躁夜夜躁夜夜| 亚洲久久久国产精品| 午夜久久久在线观看| 极品少妇高潮喷水抽搐| 黄色片一级片一级黄色片| 建设人人有责人人尽责人人享有的| 最近最新中文字幕大全免费视频| 久久精品亚洲熟妇少妇任你| 亚洲欧美精品自产自拍| 在线十欧美十亚洲十日本专区| 欧美日韩中文字幕国产精品一区二区三区 | 热99re8久久精品国产| 999精品在线视频| 午夜福利在线观看吧| 岛国在线观看网站| 国产成人精品无人区| 美国免费a级毛片| 亚洲国产精品999| 国产精品免费大片| h视频一区二区三区| 久久久久久亚洲精品国产蜜桃av| 亚洲伊人色综图| 免费观看a级毛片全部| 亚洲精品成人av观看孕妇| 欧美xxⅹ黑人| 老司机深夜福利视频在线观看 | 不卡av一区二区三区| 成人亚洲精品一区在线观看| 久久 成人 亚洲| 精品人妻熟女毛片av久久网站| 国产男女内射视频| 欧美黑人精品巨大| 欧美另类一区| 国产不卡av网站在线观看| 亚洲欧洲日产国产| av欧美777| 自拍欧美九色日韩亚洲蝌蚪91| 亚洲精品一区蜜桃| 亚洲伊人久久精品综合| videos熟女内射| 视频区图区小说| 亚洲国产成人一精品久久久| 免费观看a级毛片全部| 亚洲精品在线美女| 免费高清在线观看视频在线观看| 老司机在亚洲福利影院| 伊人亚洲综合成人网| 女性被躁到高潮视频| 香蕉丝袜av| 亚洲av片天天在线观看| 黄色视频,在线免费观看| 大片免费播放器 马上看| 久久久水蜜桃国产精品网| 日本av手机在线免费观看| 1024香蕉在线观看| 国产精品一区二区在线不卡| 美女高潮喷水抽搐中文字幕| 国产精品一区二区免费欧美 | 99国产精品99久久久久| 性高湖久久久久久久久免费观看| 午夜日韩欧美国产| 老司机影院毛片| 国精品久久久久久国模美| 亚洲成av片中文字幕在线观看| 亚洲情色 制服丝袜| 久久中文看片网| 亚洲精品国产色婷婷电影| 捣出白浆h1v1| 精品国内亚洲2022精品成人 | 国产日韩欧美视频二区| 欧美精品一区二区免费开放| 一区二区三区激情视频| 超色免费av| 侵犯人妻中文字幕一二三四区| 69av精品久久久久久 | 久久久久久久国产电影| 精品一区二区三区四区五区乱码| 丰满饥渴人妻一区二区三| 亚洲精品日韩在线中文字幕| 精品福利观看| 亚洲av国产av综合av卡| 日韩一区二区三区影片| 亚洲熟女精品中文字幕| 欧美+亚洲+日韩+国产| 日本wwww免费看| 老熟妇仑乱视频hdxx| 在线观看免费午夜福利视频| 精品亚洲乱码少妇综合久久| 亚洲精品乱久久久久久| avwww免费| 人妻久久中文字幕网| 美女中出高潮动态图| 欧美中文综合在线视频| 久久中文看片网| 久久综合国产亚洲精品| 国产一区二区激情短视频 | 久久99一区二区三区| 色老头精品视频在线观看| 国产免费现黄频在线看| 又黄又粗又硬又大视频| 亚洲欧美精品自产自拍| 欧美激情久久久久久爽电影 | 久久这里只有精品19| 女人高潮潮喷娇喘18禁视频| 精品久久蜜臀av无| 亚洲第一av免费看| 国产精品欧美亚洲77777| 国产精品香港三级国产av潘金莲| 成人影院久久| 欧美日本中文国产一区发布| 国产精品 欧美亚洲| 考比视频在线观看| 多毛熟女@视频| 国产精品久久久av美女十八| 黑人猛操日本美女一级片| a级片在线免费高清观看视频| av一本久久久久| 黄片小视频在线播放| 黄色怎么调成土黄色| 亚洲精品美女久久久久99蜜臀| 50天的宝宝边吃奶边哭怎么回事| 啦啦啦 在线观看视频| 国产成人a∨麻豆精品| 成人亚洲精品一区在线观看| 国产精品久久久久成人av| 亚洲欧美一区二区三区黑人| 欧美日韩精品网址| a在线观看视频网站| 亚洲精品中文字幕在线视频| 精品一区二区三区四区五区乱码| 免费少妇av软件| 中文精品一卡2卡3卡4更新| 国产在线观看jvid| 男女国产视频网站| 日韩,欧美,国产一区二区三区| 精品少妇内射三级| 中文字幕另类日韩欧美亚洲嫩草| 国产精品99久久99久久久不卡| 欧美亚洲 丝袜 人妻 在线| 热re99久久精品国产66热6| 亚洲色图 男人天堂 中文字幕| 少妇粗大呻吟视频| 99久久精品国产亚洲精品| 丝袜喷水一区| 精品一区在线观看国产| 日韩三级视频一区二区三区| 王馨瑶露胸无遮挡在线观看| 国产成人精品久久二区二区91| 男女高潮啪啪啪动态图| 精品一区二区三卡| 国产人伦9x9x在线观看| 亚洲国产日韩一区二区| 国产成+人综合+亚洲专区| 狠狠狠狠99中文字幕| 美女国产高潮福利片在线看| 免费高清在线观看视频在线观看| 亚洲专区中文字幕在线| 精品欧美一区二区三区在线| 国产91精品成人一区二区三区 | 女人被躁到高潮嗷嗷叫费观| 欧美黄色淫秽网站| 在线永久观看黄色视频| 青春草亚洲视频在线观看| 国产一区二区三区综合在线观看| 人妻人人澡人人爽人人| 亚洲激情五月婷婷啪啪| 在线亚洲精品国产二区图片欧美| 爱豆传媒免费全集在线观看| 亚洲视频免费观看视频| 美女福利国产在线| 69av精品久久久久久 | 久久av网站| 在线观看免费高清a一片| 午夜福利在线免费观看网站| 一边摸一边抽搐一进一出视频| 免费看十八禁软件| 18禁观看日本| 中文字幕色久视频| 在线永久观看黄色视频| 国产成人免费无遮挡视频| 叶爱在线成人免费视频播放| 亚洲专区国产一区二区| 一区在线观看完整版| 国产日韩欧美亚洲二区| 男女国产视频网站| 女性生殖器流出的白浆| 免费日韩欧美在线观看| 亚洲美女黄色视频免费看| 人人妻,人人澡人人爽秒播| 美女扒开内裤让男人捅视频| 久久久久精品人妻al黑| 我要看黄色一级片免费的| 这个男人来自地球电影免费观看| 中文字幕制服av| 精品一区二区三区四区五区乱码| 午夜福利乱码中文字幕| 精品国产乱子伦一区二区三区 | 水蜜桃什么品种好| 精品熟女少妇八av免费久了| 国产成人系列免费观看| 久久久久久久久免费视频了| 丝袜美腿诱惑在线| 午夜两性在线视频| 久久久精品国产亚洲av高清涩受| 日韩视频在线欧美| 久久综合国产亚洲精品| 久久中文字幕一级| 亚洲精品一卡2卡三卡4卡5卡 | 丝袜脚勾引网站| 亚洲天堂av无毛| av在线app专区| 一级片'在线观看视频| 国产伦理片在线播放av一区| 一区二区三区乱码不卡18| 免费看十八禁软件| 人人妻,人人澡人人爽秒播| 亚洲欧洲日产国产| 久久亚洲国产成人精品v| 久久热在线av| 91成年电影在线观看| 男人爽女人下面视频在线观看| av超薄肉色丝袜交足视频| 亚洲精品国产区一区二| av又黄又爽大尺度在线免费看| 一本一本久久a久久精品综合妖精| 黑人操中国人逼视频| 性色av乱码一区二区三区2| 五月天丁香电影| 麻豆国产av国片精品| 免费少妇av软件| 高潮久久久久久久久久久不卡| 黄色视频,在线免费观看| 亚洲一区中文字幕在线| 久久九九热精品免费| 热re99久久精品国产66热6| 首页视频小说图片口味搜索| 人妻一区二区av| 国产亚洲精品一区二区www | 亚洲七黄色美女视频| 国产精品久久久久久精品古装| 爱豆传媒免费全集在线观看| 在线精品无人区一区二区三| netflix在线观看网站| 18在线观看网站| 国产在线一区二区三区精| www.精华液| 日韩熟女老妇一区二区性免费视频| 精品福利观看| 久久久精品国产亚洲av高清涩受| 久久毛片免费看一区二区三区| 亚洲av成人一区二区三| a级毛片在线看网站| 国产深夜福利视频在线观看| 亚洲精品国产色婷婷电影| 国产精品自产拍在线观看55亚洲 | 亚洲精品国产av成人精品| 又大又爽又粗| 亚洲精品一区蜜桃| 成人18禁高潮啪啪吃奶动态图| 91成年电影在线观看| 国产精品1区2区在线观看. | 亚洲av电影在线进入| 侵犯人妻中文字幕一二三四区| 精品久久久久久电影网| 亚洲,欧美精品.| 亚洲自偷自拍图片 自拍| 91国产中文字幕| 国产精品久久久av美女十八| 亚洲 欧美一区二区三区| 国产精品久久久久久人妻精品电影 | 精品亚洲乱码少妇综合久久| 中国国产av一级| 国产av精品麻豆| 麻豆乱淫一区二区| 18禁黄网站禁片午夜丰满| 日韩欧美一区二区三区在线观看 | 老司机午夜十八禁免费视频| 久久精品aⅴ一区二区三区四区| 亚洲色图综合在线观看| 交换朋友夫妻互换小说| 国产深夜福利视频在线观看| 亚洲av电影在线观看一区二区三区| 肉色欧美久久久久久久蜜桃| 12—13女人毛片做爰片一| 亚洲欧洲日产国产| 亚洲性夜色夜夜综合| 国产无遮挡羞羞视频在线观看| 成人三级做爰电影| 午夜福利免费观看在线| 一级黄色大片毛片| 精品国产乱码久久久久久小说| 国产免费av片在线观看野外av| 国产亚洲精品一区二区www | 国产成人精品久久二区二区91| 天天躁夜夜躁狠狠躁躁| 中文字幕最新亚洲高清| 亚洲九九香蕉| 成人av一区二区三区在线看 | 亚洲国产日韩一区二区| 亚洲精品一卡2卡三卡4卡5卡 | 老汉色∧v一级毛片| 最新在线观看一区二区三区| 国产欧美日韩一区二区三 | 久久av网站| 色老头精品视频在线观看| 国产麻豆69| 黄片小视频在线播放| 欧美激情久久久久久爽电影 | 亚洲情色 制服丝袜| 大型av网站在线播放| 纯流量卡能插随身wifi吗| 中文字幕最新亚洲高清| 肉色欧美久久久久久久蜜桃| 亚洲精品国产av蜜桃| 久久国产精品人妻蜜桃| 俄罗斯特黄特色一大片|