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

    基于任務(wù)單元的運(yùn)糧車響應(yīng)調(diào)度優(yōu)化

    2020-11-16 07:26:12曹光喬張慶凱
    關(guān)鍵詞:運(yùn)糧非生產(chǎn)性等待時(shí)間

    曹光喬 南 風(fēng),2 陳 聰 張慶凱

    (1.農(nóng)業(yè)農(nóng)村部 南京農(nóng)業(yè)機(jī)械化研究所,南京 210014;2.中國農(nóng)業(yè)科學(xué)院 研究生院,北京 463000;3.西北農(nóng)林科技大學(xué) 機(jī)械與電子工程學(xué)院,陜西 楊凌 712100)

    目前我國新興經(jīng)營主體和適度規(guī)模的田間收獲作業(yè)服務(wù)需求不斷增長(zhǎng)[1],對(duì)農(nóng)機(jī)配置與路徑規(guī)劃的要求也在不斷提高,除了收割機(jī)單獨(dú)的調(diào)度模式以外,服務(wù)組織也迫切需要精確合理的運(yùn)糧車響應(yīng)方案。研究多約束環(huán)境的復(fù)雜調(diào)度問題,擺脫依靠傳統(tǒng)的人為經(jīng)驗(yàn),獲得科學(xué)合理的運(yùn)糧車響應(yīng)路徑優(yōu)化方案,對(duì)節(jié)約農(nóng)機(jī)資源以及提升糧食收獲和運(yùn)輸整體效率具有重要意義[2-3]。

    國內(nèi)外大型農(nóng)場(chǎng)的集中收獲作業(yè),收割機(jī)與運(yùn)糧車調(diào)度有3種響應(yīng)方式:運(yùn)糧車田間跟隨收割機(jī)進(jìn)行作業(yè);收割機(jī)專門配備固定數(shù)量的運(yùn)糧車進(jìn)行卸糧作業(yè);收割機(jī)選擇定點(diǎn)卸糧的位置,運(yùn)糧車及時(shí)到達(dá)[4-5]。這3種方式均不需要考慮多臺(tái)收割機(jī)與運(yùn)糧車的復(fù)雜響應(yīng),很少針對(duì)運(yùn)糧車制定精確路徑規(guī)劃。對(duì)于適度規(guī)模的新興經(jīng)營主體,服務(wù)組織無法為1臺(tái)收割機(jī)專門配備 2 輛以上的運(yùn)糧車,尤其在農(nóng)忙時(shí)節(jié),收割機(jī)與運(yùn)糧車數(shù)量不匹配情況普遍存在。因此,類似大型農(nóng)場(chǎng)的冗余運(yùn)糧車與收割機(jī)響應(yīng)調(diào)度模式和方法并不可取,服務(wù)組織缺乏運(yùn)輸設(shè)備不足環(huán)境下的調(diào)度規(guī)劃方案。

    隨著我國農(nóng)村勞動(dòng)力大量轉(zhuǎn)移,土地流轉(zhuǎn)和作業(yè)服務(wù)托管加速,適度規(guī)模經(jīng)營成為現(xiàn)代農(nóng)業(yè)發(fā)展的必然趨勢(shì)[6-9]。農(nóng)作物收割和運(yùn)輸作業(yè),具有很強(qiáng)的時(shí)效性,收割機(jī)和運(yùn)糧車作業(yè)調(diào)度缺乏詳細(xì)規(guī)劃和技術(shù)手段,其效率很大程度依賴于人工操作的經(jīng)驗(yàn)性[10]。家庭農(nóng)場(chǎng)、綜合農(nóng)事服務(wù)中心等模式興盛發(fā)展,其擁有的農(nóng)田數(shù)量和經(jīng)營規(guī)模不斷增加,但仍然具有小規(guī)模和零散化的特點(diǎn)[11],主要依靠農(nóng)機(jī)服務(wù)組織提供生產(chǎn)性農(nóng)事服務(wù)[12-14]。農(nóng)業(yè)環(huán)境小規(guī)模、零散化的特點(diǎn)也決定了收運(yùn)響應(yīng)需要更加靈活、科學(xué)的匹配機(jī)制,以及對(duì)運(yùn)糧車詳細(xì)、精確的調(diào)度規(guī)劃。

    國內(nèi)外針對(duì)農(nóng)業(yè)機(jī)械調(diào)度問題多集中于田間農(nóng)機(jī)資源的調(diào)配和路徑優(yōu)化問題,圍繞聯(lián)合收割機(jī)和大型農(nóng)田任務(wù)等2個(gè)對(duì)象進(jìn)行資源調(diào)度,重點(diǎn)研究聯(lián)合收割機(jī)的路徑轉(zhuǎn)移。圍繞不同類型的農(nóng)機(jī)調(diào)度問題,張璠等[15]提出了兩種基于優(yōu)先策略的多機(jī)多任務(wù)緊急調(diào)配算法,建立了以調(diào)配成本和損失最小化為目標(biāo)的緊急調(diào)配模型;吳才聰?shù)萚16-17]建立了帶時(shí)間窗約束的農(nóng)機(jī)資源調(diào)度數(shù)學(xué)模型,并以動(dòng)態(tài)規(guī)劃的思想進(jìn)行模型求解。圍繞農(nóng)機(jī)調(diào)度算法問題,許多研究認(rèn)為農(nóng)機(jī)調(diào)度即為車輛調(diào)度,提出了一些相關(guān)算法,取得了良好的成效。例如 Bochtis D D等[18]針對(duì)農(nóng)業(yè)領(lǐng)域的車輛調(diào)度問題,使用 VRP 的方法解決農(nóng)業(yè)車隊(duì)管理問題,增強(qiáng)傳統(tǒng)的農(nóng)機(jī)管理系統(tǒng);Guan S等[19]提出一種兩階段的農(nóng)機(jī)調(diào)度方案,首先對(duì)農(nóng)機(jī)資源進(jìn)行分配,然后利用遺傳算法與模擬退火算法相結(jié)合對(duì)調(diào)度模型進(jìn)行求解,同時(shí)他們也將農(nóng)業(yè)生產(chǎn)過程劃分為2個(gè)階段,采用優(yōu)先級(jí)原則對(duì)農(nóng)機(jī)資源進(jìn)行合理分配,利用遺傳算法實(shí)現(xiàn)農(nóng)機(jī)合理調(diào)度,提高資源利用率[20]。圍繞收割機(jī)與運(yùn)糧車協(xié)同調(diào)度問題,He P F等[21]將收割機(jī)收獲和運(yùn)糧車運(yùn)輸問題考慮為兩級(jí)多行程問題,根據(jù)運(yùn)糧車需要響應(yīng)的具體位置動(dòng)態(tài)變化,從而提出混合整數(shù)新型規(guī)劃模型,并提出一種啟發(fā)式算法進(jìn)行求解,屬于概念性的調(diào)度響應(yīng)問題[22];Ali O等[4]將收割機(jī)和運(yùn)糧車的響應(yīng)問題轉(zhuǎn)化為內(nèi)場(chǎng)物流的車輛路徑問題,模型適用于大面積農(nóng)田的場(chǎng)內(nèi)路徑規(guī)劃,無法解決運(yùn)糧車的場(chǎng)外路徑規(guī)劃。圍繞農(nóng)機(jī)具收獲和運(yùn)輸協(xié)同調(diào)度的研究[23-27],國內(nèi)外更多的是集中于作業(yè)理論的系統(tǒng)分析和農(nóng)機(jī)具數(shù)字化監(jiān)控系統(tǒng)和信息平臺(tái)的設(shè)計(jì)開發(fā),但是不能為用戶提供相關(guān)的調(diào)度決策。收割機(jī)和運(yùn)糧車響應(yīng)問題屬于比較復(fù)雜的多目標(biāo)組合優(yōu)化問題[28-29],在已有的少數(shù)響應(yīng)調(diào)度研究中,多將其轉(zhuǎn)化為帶時(shí)間窗的旅行商問題、帶時(shí)間窗的VPR問題或者運(yùn)籌學(xué)的運(yùn)輸問題,僅僅考慮作業(yè)隊(duì)行走路線,未結(jié)合時(shí)間因素和空間因素。

    結(jié)合農(nóng)業(yè)環(huán)境小規(guī)模零散化的特點(diǎn),收割機(jī)與運(yùn)糧車響應(yīng)模式研究與科學(xué)合理的運(yùn)糧車調(diào)度規(guī)劃研究較為薄弱,探究收割機(jī)與運(yùn)糧車的協(xié)同響應(yīng)成為啟動(dòng)本項(xiàng)研究的依據(jù)。

    本研究擬將任務(wù)單元作為收割機(jī)作業(yè)的界限,分析適度規(guī)模場(chǎng)景下的收、運(yùn)響應(yīng)作業(yè)模式,運(yùn)用動(dòng)態(tài)規(guī)劃思想去解決靜態(tài)規(guī)劃問題,以期為運(yùn)糧車科學(xué)調(diào)度提供優(yōu)化方案與決策依據(jù),提高運(yùn)輸裝備的利用效率。

    1 運(yùn)糧車與收割機(jī)響應(yīng)模型的建立

    1.1 任務(wù)單元

    本研究將包含收割機(jī)的農(nóng)作物收獲區(qū)域定義為任務(wù)單元,單臺(tái)收割機(jī)在一定規(guī)模的該區(qū)域獨(dú)立完成收獲作業(yè)。任務(wù)單元彼此間相互獨(dú)立,為同一主體所有。同一主體是指經(jīng)營農(nóng)田的種糧大戶或者家庭農(nóng)場(chǎng),為了方便收獲安排,將該主體經(jīng)營的所有土地劃分為若干任務(wù)單元。同一主體能夠滿足運(yùn)糧車進(jìn)行往復(fù)循環(huán)使用的要求,承擔(dān)多個(gè)任務(wù)單元的糧食載運(yùn)。任務(wù)單元作為收獲任務(wù)的分解,可以是零散分布的農(nóng)田或者是大面積農(nóng)田的某個(gè)劃分區(qū)域。

    運(yùn)糧車的響應(yīng)對(duì)象是收割機(jī),而收割機(jī)的純收割時(shí)間及運(yùn)糧車響應(yīng)次數(shù)又與農(nóng)田產(chǎn)量息息相關(guān),為了簡(jiǎn)化收運(yùn)響應(yīng)機(jī)制,減少多種對(duì)象的復(fù)雜約束,本研究將收割機(jī)與農(nóng)田區(qū)域逐一對(duì)應(yīng),即以任務(wù)單元作為收割機(jī)作業(yè)的界限,1個(gè)任務(wù)單元只由1臺(tái)收割機(jī)進(jìn)行作業(yè)。經(jīng)過實(shí)際調(diào)研,江蘇省 3.0~13.5 hm2(50~200畝)的家庭農(nóng)場(chǎng)和種糧大戶等主體,80% 的單田面積為0.33~1.0 hm2,因此選取適度規(guī)模經(jīng)營0.33~1.0 hm2為任務(wù)單元面積。

    1.2 響應(yīng)模式

    收割機(jī)與運(yùn)糧車響應(yīng)調(diào)度描述為:農(nóng)機(jī)服務(wù)組織為同一經(jīng)營主體所屬的任務(wù)單元服務(wù),如何實(shí)現(xiàn)較少運(yùn)糧車匹配收割機(jī)從而完成收獲、運(yùn)輸任務(wù)的調(diào)度規(guī)劃問題。

    本研究選取當(dāng)前時(shí)刻收割機(jī)在任務(wù)單元中的作業(yè)狀態(tài)(收割作業(yè)狀態(tài)、收割完成狀態(tài)、等待卸糧狀態(tài)),進(jìn)行模型建立和求解。由于任務(wù)單元面積較小,運(yùn)糧車在單元內(nèi)抵達(dá)響應(yīng)位置時(shí)所行駛的路徑要遠(yuǎn)小于運(yùn)糧車在單元外的轉(zhuǎn)移路徑,故忽略收割機(jī)和運(yùn)糧車在單元內(nèi)的行走。若同時(shí)考慮收割機(jī)、運(yùn)糧車兩者的路徑規(guī)劃,模式動(dòng)態(tài)且復(fù)雜,同時(shí)需要大量實(shí)時(shí)性數(shù)據(jù),故忽略收割機(jī)路徑規(guī)劃,簡(jiǎn)化收割機(jī)與運(yùn)糧車的響應(yīng)復(fù)雜程度,實(shí)現(xiàn)運(yùn)糧車整體的調(diào)度規(guī)劃。

    對(duì)收割機(jī)已知收割效率、收割周期等作業(yè)參數(shù)信息進(jìn)行新狀態(tài)的判斷和預(yù)測(cè),以及利用傳感器等實(shí)時(shí)性傳輸設(shè)備的反饋和糾正獲取收割機(jī)已收割量、作業(yè)時(shí)間等實(shí)時(shí)性信息。每個(gè)任務(wù)單元對(duì)應(yīng)1臺(tái)收割機(jī),糧食滿倉后,進(jìn)行完全卸糧,完成響應(yīng)后再進(jìn)行下一輪作業(yè)。所有運(yùn)糧車先前往最近的任務(wù)單元進(jìn)行響應(yīng),響應(yīng)結(jié)束后再抵達(dá)下個(gè)任務(wù)單元,最大化減少收割機(jī)非生產(chǎn)性作業(yè)時(shí)間[4](收割機(jī)未收割情況下在田間的等待時(shí)間)。每個(gè)任務(wù)單元可以由多個(gè)運(yùn)輛車同時(shí)到達(dá),以保證收割機(jī)收滿后進(jìn)行完全卸糧作業(yè)。

    本研究基于任務(wù)單元的收、運(yùn)響應(yīng)調(diào)度模式,對(duì)模型做出如下假設(shè):所有任務(wù)單元的機(jī)械設(shè)備作業(yè)參數(shù)相同,且作業(yè)過程中無故障發(fā)生;收割機(jī)連續(xù)作業(yè)至滿倉,完全卸糧后再進(jìn)行下次作業(yè);運(yùn)糧車達(dá)到滿載,才能返回卸糧站,再進(jìn)行下次分配;運(yùn)糧車運(yùn)輸?shù)诌_(dá)同一終點(diǎn)位置;忽略運(yùn)糧車和收割機(jī)折舊費(fèi)用、人工費(fèi)用等生產(chǎn)性成本。

    1.3 參數(shù)及變量說明

    1)任務(wù)單元用集合U={Vi,ai} 表示,其中Vi={V1,V2,…,Vm}、ai={a1,a2,…,am} 分別表示任務(wù)單元位置、任務(wù)單元面積,i={1,2,…,m} 為任務(wù)單元編號(hào),m為任務(wù)單元數(shù)量。

    2)收割機(jī)與運(yùn)糧車分別用集合Hi={H1,H2,…,Hm}、Tj={T1,T2,…,Tn} 表示,其中i={1,2,…,m} 為任務(wù)單元編號(hào);n為運(yùn)糧車數(shù)量,j={1,2,…,n} 為運(yùn)糧車編號(hào);令CH表示收割機(jī)容量,CT表示運(yùn)糧車容量。

    3)定義集合P={Vp,dgh} 表示路網(wǎng)中各路徑節(jié)點(diǎn)路徑信息。其中運(yùn)糧車分配節(jié)點(diǎn)的位置Vp=V0∪Vi,V0表示卸糧站,Vi表示任務(wù)單元位置;dgh表示節(jié)點(diǎn)g到h之間的距離,g∈Vp,h∈Vp。

    1.4 時(shí)間槽定義

    本研究基于每個(gè)任務(wù)單元的整體起止時(shí)間,進(jìn)行時(shí)間槽劃分。時(shí)間槽是指單個(gè)任務(wù)單元收割作業(yè)周期,該周期包含收割機(jī)滿倉前的收割狀態(tài)和滿倉后的卸糧響應(yīng)狀態(tài),是2個(gè)狀態(tài)時(shí)間段的總合。劃分依據(jù)是任務(wù)單元中收割機(jī)的作業(yè)周期,數(shù)目與收割機(jī)的收割總次數(shù)Ki相同,計(jì)算公式為:

    (1)

    式中:ai表示i任務(wù)單元的收獲面積;λ是作物面積轉(zhuǎn)化為產(chǎn)量的單位系數(shù);λai表示該任務(wù)單元總產(chǎn)量。i任務(wù)單元含有Ki個(gè)時(shí)間槽表示該任務(wù)單元必須被收割λai/CH(向上取整)次,才能完全收獲完畢。任務(wù)單元最后一次收割時(shí)需完全卸糧,保證離開該任務(wù)單元的收割機(jī)糧倉處于全空狀態(tài),能夠滿足任務(wù)單元作業(yè)結(jié)束后自由轉(zhuǎn)移的需要。

    (2)

    分別表示收割機(jī)的開始收割時(shí)間、滿倉時(shí)間、卸糧結(jié)束時(shí)間,分別表示運(yùn)糧車到達(dá)時(shí)間、開始時(shí)間、結(jié)束時(shí)間;為收割機(jī)收割作業(yè)時(shí)間,為運(yùn)糧車卸糧響應(yīng)時(shí)間;k,i,j分別為時(shí)間槽、收割機(jī)、運(yùn)糧車序號(hào)。 is time for unloading of grain by grain truck; k, i and j are the serial numbers of time slot, harvester and grain truck, respectively.圖1 運(yùn)糧車與收割機(jī)響應(yīng)時(shí)間狀態(tài)Fig.1 Response time status of grain trucks and harvester

    (3)

    1.5 響應(yīng)模型的建立

    1.5.1主要目標(biāo)

    基于任務(wù)單元的響應(yīng)模式及相關(guān)影響因素分析,收運(yùn)響應(yīng)調(diào)度的非生產(chǎn)性總成本包含收割機(jī)等待時(shí)間成本、運(yùn)糧車轉(zhuǎn)移距離成本2部分。因此將模型分為2種具體目標(biāo),分別為收割機(jī)非生產(chǎn)性作業(yè)時(shí)間(等待時(shí)間)成本最小化和運(yùn)糧車的轉(zhuǎn)移距離成本最小化。

    1)收割機(jī)非生產(chǎn)性作業(yè)時(shí)間成本最小化。模型目標(biāo)數(shù)學(xué)表達(dá)式為:

    (4)

    2)運(yùn)糧車轉(zhuǎn)移距離成本最小化。模型目標(biāo)數(shù)學(xué)表達(dá)式為:

    (5)

    式中:f2代表運(yùn)糧車轉(zhuǎn)移成本最小值;c2為運(yùn)糧車單位距離成本。轉(zhuǎn)移成本等于所有運(yùn)糧車轉(zhuǎn)移總距離與單位距離成本之積。

    1.5.2主要約束條件

    對(duì)運(yùn)糧車收割機(jī)響應(yīng)調(diào)度過程分析,工作量約束條件如下:

    (6)

    時(shí)間狀態(tài)約束條件為:

    (7)

    2 運(yùn)糧車與收割機(jī)響應(yīng)模型的解算

    2.1 解算基本原理

    本研究首先考慮時(shí)間約束,在時(shí)間約束下完成有限次的調(diào)度決策。在作業(yè)期間內(nèi),將單個(gè)任務(wù)單元的整體起始時(shí)間切割成有限時(shí)間槽,每個(gè)時(shí)間槽根據(jù)當(dāng)前狀態(tài)完成1次調(diào)度決策。所有任務(wù)單元并行計(jì)算,最終完成所有單元的調(diào)度安排,獲得完整的運(yùn)糧車調(diào)度方案。

    按照時(shí)間順序依次完成各個(gè)時(shí)間槽的決策,最終完成整體解算,從而順利解決時(shí)間和工作量約束的問題。每個(gè)時(shí)間槽在單獨(dú)決策時(shí),事先確定目標(biāo)之間的相對(duì)重要程度,以此為依據(jù)將重要目標(biāo)進(jìn)行優(yōu)先決策,次要目標(biāo)隨后決策,即多目標(biāo)問題轉(zhuǎn)化為2個(gè)單目標(biāo)問題進(jìn)行順序求解。本研究將基于收割機(jī)非生產(chǎn)性時(shí)間為首要目標(biāo)的動(dòng)態(tài)時(shí)間槽算法記為 A,將基于運(yùn)糧車運(yùn)輸距離為首要目標(biāo)的動(dòng)態(tài)時(shí)間槽算法記為 B。

    2.2 算法構(gòu)建

    動(dòng)態(tài)時(shí)間槽算法 A 中運(yùn)糧車運(yùn)輸距離和收割機(jī)等待時(shí)間同時(shí)作為優(yōu)化目標(biāo),決策中收割機(jī)等待時(shí)間處于更高優(yōu)先級(jí)。具體分為針對(duì)單個(gè)任務(wù)單元的決策算法 A1和并行計(jì)算多個(gè)任務(wù)單元的算法 A2,簡(jiǎn)化步驟如下:

    1)算法 A1。

    步驟1:所有任務(wù)單元同時(shí)開始計(jì)算,每次迭代需更新時(shí)間集合T。每個(gè)任務(wù)單元都從第一個(gè)時(shí)間槽開始計(jì)算(k=1)。

    步驟2:判斷運(yùn)糧車狀態(tài)。將所有運(yùn)糧車劃分為等待分配和已分配狀態(tài)。只有等待分配的運(yùn)糧車才能參與各個(gè)任務(wù)單元的匹配。

    步驟3:判斷收割機(jī)狀態(tài)。若該任務(wù)單元第k個(gè)時(shí)間槽下收割機(jī)處于卸糧未完成狀態(tài),直接執(zhí)行算法A2;若收割機(jī)對(duì)應(yīng)收割狀態(tài),進(jìn)入下一步。

    步驟4:優(yōu)化分配。計(jì)算當(dāng)前任務(wù)單元第k個(gè)時(shí)間槽剩余的收割作業(yè)時(shí)間,計(jì)算所有未分配運(yùn)糧車到該任務(wù)單元的距離。執(zhí)行算法A2。

    步驟5:按照算法A2指定結(jié)果分配。若未能分配,則轉(zhuǎn)入步驟3。

    步驟6:判斷該任務(wù)單元作業(yè)是否全部完成(k≥Ki)。若未完成,令k增加1后轉(zhuǎn)入步驟1,再次循環(huán)。若完成,記錄并結(jié)束該任務(wù)單元。

    步驟7:若所有任務(wù)單元都被記錄,結(jié)束算法A1。

    2)算法A2。

    步驟1:優(yōu)先按時(shí)間分配。匯總當(dāng)前所有單元收割機(jī)剩余收割時(shí)間,按照從小到大的順序依次分配即將結(jié)束收割的任務(wù)單元。

    步驟2:尋找運(yùn)糧車。預(yù)算所有單元當(dāng)前第k個(gè)時(shí)間槽需要的卸糧量,并與所有運(yùn)糧車剩余載量比較,找出滿足條件的運(yùn)糧車。

    步驟3:再次按距離分配。對(duì)步驟2滿足載量條件的運(yùn)糧車,取距離該任務(wù)單元位置最近的1個(gè)運(yùn)糧車,直接執(zhí)行步驟5。

    步驟4:組合判斷。步驟2中若無滿足載重條件的運(yùn)糧車,則將未參與的運(yùn)糧車多個(gè)組合,直到滿足;若全部不滿足,則所有運(yùn)糧車共同派往,并將該任務(wù)單元標(biāo)記為分配優(yōu)先級(jí)最高。

    步驟5:結(jié)束A2,輸出分配結(jié)果。

    此外,本研究在算法A的基礎(chǔ)上提出算法B,運(yùn)糧車運(yùn)輸距離和收割機(jī)等待時(shí)間同時(shí)作為優(yōu)化目標(biāo),決策中運(yùn)糧車轉(zhuǎn)移距離處于更高優(yōu)先級(jí)。將算法A2中步驟1、3的優(yōu)先分配規(guī)則互換,即可得算法 B。

    3 運(yùn)算分析

    國家統(tǒng)計(jì)局2018年主要農(nóng)作物產(chǎn)量統(tǒng)計(jì)年鑒數(shù)據(jù)顯示,江蘇省鹽城市谷物小麥每1 hm2產(chǎn)糧5.827 t,即λ=5.827 t/hm2。根據(jù)江蘇省常用品牌的收割機(jī)和運(yùn)糧車參數(shù),收割機(jī)糧容積為1.05~4.8 m3,本研究選用雷沃谷神GE80(4LZ-8E2)輪式谷物聯(lián)合收割機(jī),糧箱容積2.0 m3,容量1.6 t,即CH=1.6,收割效率0.533 hm2/h,數(shù)量為6輛。運(yùn)糧車一般為拖拉機(jī)自帶5~7 t容量的拖斗,本研究取值5 t,即CT=5,轉(zhuǎn)移速率為30 km/h,數(shù)量為4輛。設(shè)置收割機(jī)單位等待時(shí)間成本為70 元/h(c1=70);運(yùn)糧車單位轉(zhuǎn)移成本為4 元/h(c2=4);響應(yīng)卸糧速率為38 t/h。本次研究選取鹽城市大有縣6個(gè)農(nóng)田作業(yè)點(diǎn)實(shí)際數(shù)據(jù)進(jìn)行模擬,具體信息見表1。

    表1 任務(wù)單元實(shí)際作業(yè)信息與時(shí)間槽數(shù)目Table 1 Actual operating data and number of time slots for work units

    3.1 目標(biāo)優(yōu)化結(jié)果

    3.1.1實(shí)際算例方案

    算法A與算法B分別計(jì)算鹽城大有縣6個(gè)農(nóng)田作業(yè)點(diǎn)的實(shí)際算例,得到2種不同結(jié)果的優(yōu)化調(diào)度方案。方案呈現(xiàn)的收割機(jī)等待時(shí)間、運(yùn)糧車轉(zhuǎn)移距離及整體非生產(chǎn)性作業(yè)成本等目標(biāo)優(yōu)化結(jié)果見表2。

    表2 算法A、B對(duì)實(shí)際案例的目標(biāo)優(yōu)化結(jié)果Table 2 Algorithm A and algorithm B optimize the objective results of the actual case

    算法A得到此案例的方案在收割機(jī)等待時(shí)間上較少,算法 B得到運(yùn)糧車轉(zhuǎn)移距離較短,優(yōu)化結(jié)果符合不同策略算法A、B的預(yù)期?;诙嗄繕?biāo)動(dòng)態(tài)分割求解算法 A 產(chǎn)生方案中非生產(chǎn)性作業(yè)總成本較低,對(duì)該調(diào)度模型求解結(jié)果較好。因此本研究選擇算法 A 作為該實(shí)例的求解算法,獲得優(yōu)化方案的調(diào)度響應(yīng)甘特圖(圖2)與運(yùn)糧車到訪任務(wù)單元的路徑順序(表3)。

    紅色線框?yàn)槭崭顧C(jī)與運(yùn)糧車實(shí)際卸糧響應(yīng)時(shí)間,右側(cè)數(shù)值表示到訪的運(yùn)糧車編號(hào)。The red line represents the unloading status. The value represents the number of the grain truck responding.圖2 任務(wù)單元調(diào)度響應(yīng)甘特圖Fig.2 The scheduling response Gantt diagram for work units

    表3 運(yùn)糧車最佳路徑順序Table 3 Optimal path scheme for grain trucks to visit work units

    3.1.2仿真算例運(yùn)算分析

    為比較不同優(yōu)先決策順序的算法A、B性能差異,本研究隨機(jī)生成8 組仿真算例,對(duì)非生產(chǎn)性作業(yè)總成本、運(yùn)糧車轉(zhuǎn)移距離和收割機(jī)等待時(shí)間進(jìn)行對(duì)比分析。每組算例包含6 個(gè)任務(wù)單元,每個(gè)任務(wù)單元面積范圍為0.33~1 hm2,單元間節(jié)點(diǎn)距離范圍為0.3~2 km。算法A、B計(jì)算每組算例得到的目標(biāo)優(yōu)化結(jié)果見表4??梢? 算法 A 求解的目標(biāo)優(yōu)化結(jié)果中非生產(chǎn)性作業(yè)成本比B低2%~18%,收割機(jī)等待時(shí)間比B低10%~40%;算法A、B中運(yùn)糧車轉(zhuǎn)移距離值的差異范圍在10%以內(nèi)。算法B求解的某些轉(zhuǎn)移距離大于A的原因是:算法B優(yōu)先考慮運(yùn)糧車轉(zhuǎn)移距離,會(huì)使運(yùn)糧車集中于較近單元,最后完成距離較遠(yuǎn)單元任務(wù)時(shí),容量已無法滿足卸糧需要,需多輛共同前往,導(dǎo)致運(yùn)糧車轉(zhuǎn)移總距離反而增加。

    表4 算法A、B對(duì)仿真算例的目標(biāo)優(yōu)化結(jié)果Table 4 Comparison between two algorithms with 8 sets of work units

    因此,決策優(yōu)先考慮等待時(shí)間的算法A總體上優(yōu)于優(yōu)先考慮轉(zhuǎn)移距離的算法B。

    3.2 影響因素分析

    以算法 A 為運(yùn)糧車和收割機(jī)響應(yīng)調(diào)度模型的優(yōu)化算法,進(jìn)行影響因素仿真試驗(yàn)。

    除了任務(wù)單元和運(yùn)糧車數(shù)量以外,響應(yīng)調(diào)度模型還包括收割效率、卸糧速率以及收割機(jī)糧倉容量和運(yùn)糧車載量等作業(yè)參數(shù)。收割機(jī)卸糧較快,卸糧時(shí)間相對(duì)于收割作業(yè)周期可以忽略;收割機(jī)糧倉容積與收割效率均通過影響收割機(jī)收割作業(yè)周期長(zhǎng)短進(jìn)而影響優(yōu)化結(jié)果,作用效果相似。因此本研究選擇收割機(jī)、運(yùn)糧車數(shù)量配比以及收割效率作業(yè)參數(shù),探究其對(duì)模型性能的影響。

    3.2.1運(yùn)糧車與收割機(jī)數(shù)量配比的影響

    本研究固定6臺(tái)收割機(jī)(任務(wù)單元),選擇運(yùn)糧車數(shù)量為1~6輛(收運(yùn)數(shù)量配比1∶1~6∶1)。對(duì)每組任務(wù)單元算例(隨機(jī)生成6 個(gè)任務(wù)單元數(shù)據(jù),單元面積范圍為0.33~1 hm2,單元間節(jié)點(diǎn)距離范圍為0.3~2 km)進(jìn)行求解,運(yùn)行結(jié)果見圖3。

    由圖3可知,在收割機(jī)數(shù)量固定且多于運(yùn)糧車的情況下,隨著運(yùn)糧車數(shù)量不斷增加,收割機(jī)非生產(chǎn)性作業(yè)平均成本、平均等待時(shí)間以及單個(gè)運(yùn)糧車的平均轉(zhuǎn)移路程均不斷下降,總成本、總等待時(shí)間都隨著運(yùn)糧車增加而不斷下降。當(dāng)運(yùn)糧車數(shù)量為4~6輛(收運(yùn)配比范圍在1∶1~2∶1)時(shí),每臺(tái)收割機(jī)非生產(chǎn)性平均等待時(shí)間均在0.05 h以內(nèi),優(yōu)化效果良好(圖3(b))。固定6個(gè)任務(wù)單元情況下,服務(wù)組織配置6輛運(yùn)糧車與配置1輛運(yùn)糧車相比,非生產(chǎn)性總成本、平均等待時(shí)間和轉(zhuǎn)移距離減少較為顯著,增益較大。服務(wù)組織配置6輛運(yùn)糧車與配置4輛運(yùn)糧車相比,減少的平均等待時(shí)間與轉(zhuǎn)移距離減少幅度不明顯,運(yùn)糧車數(shù)量增加對(duì)非生產(chǎn)性總成本的降低不明顯,增大了對(duì)運(yùn)輸設(shè)備資源的占用,即收、運(yùn)數(shù)量相近時(shí)提高運(yùn)糧車數(shù)量帶來增益減小。

    圖3 運(yùn)糧車數(shù)量變化對(duì)目標(biāo)優(yōu)化結(jié)果的影響Fig.3 Influence of the quantity change of grain cart on the optimization results

    因此在實(shí)際情況下服務(wù)組織選擇4輛運(yùn)糧車能取得較好的優(yōu)化結(jié)果,并節(jié)約運(yùn)糧車數(shù)量。尤其是對(duì)比當(dāng)前大型農(nóng)場(chǎng)運(yùn)糧車冗余配置的情況,本研究?jī)?nèi)容對(duì)節(jié)約運(yùn)糧車資源具有一定的指導(dǎo)意義。

    同時(shí),固定運(yùn)糧車4輛,任務(wù)單元(或收割機(jī))數(shù)量變化(4~20),得到如下結(jié)論:在運(yùn)糧車數(shù)輛固定且少于收割機(jī)的情況下,隨著收割機(jī)數(shù)量增加,平均非生產(chǎn)性作業(yè)成本和平均等待時(shí)間呈上升趨勢(shì),總路程、等待時(shí)間和非生產(chǎn)性總成本也隨之增加,趨勢(shì)明顯。即當(dāng)運(yùn)糧車數(shù)量不足時(shí),盲目增加收割機(jī)數(shù)量,會(huì)導(dǎo)致單個(gè)收割機(jī)的平均等待時(shí)間、非生產(chǎn)性作業(yè)成本和運(yùn)糧車平均轉(zhuǎn)移路程增加,且收割機(jī)數(shù)量越多,各指標(biāo)增加趨勢(shì)更加顯著。

    3.2.2收割效率對(duì)模型的影響

    研究模型以收割機(jī)等待時(shí)間、運(yùn)糧車轉(zhuǎn)移距離為目標(biāo)。為證明該目標(biāo)模型的有效性,設(shè)計(jì)1種常見的單目標(biāo)(僅考慮最優(yōu)等待時(shí)間)模型,產(chǎn)生方案與多目標(biāo)模型方案優(yōu)化結(jié)果進(jìn)行對(duì)比。選擇收割效率0.13~0.80 hm2/h,固定其他參數(shù),探究收運(yùn)配比為3∶2(收割機(jī)6輛,運(yùn)糧車4輛)條件下收割效率對(duì)2種調(diào)度模型的影響,結(jié)果見圖4。

    由圖4可知,當(dāng)收割效率小于0.53 hm2/h,2種模式下數(shù)據(jù)重合,即方案規(guī)劃相同,等待時(shí)間、轉(zhuǎn)移路徑相同。當(dāng)效率大于0.53 hm2/h,2種模式呈現(xiàn)差異:同時(shí)考慮等待時(shí)間與轉(zhuǎn)移距離的多目標(biāo)模型產(chǎn)生方案中,隨著收割效率增加,運(yùn)糧車轉(zhuǎn)移路徑不再變化(圖4(b));單目標(biāo)傳統(tǒng)模型產(chǎn)生方案中,等待時(shí)間、轉(zhuǎn)移距離、非生產(chǎn)性成本均隨收割效率增加而增大。可以判斷,多目標(biāo)模型在收割效率大于0.53 hm2/h后優(yōu)于傳統(tǒng)單目標(biāo)模型,本研究建立多目標(biāo)模型的有效性得以證實(shí)。

    收割效率小于0.53 hm2/h時(shí),等待時(shí)間最短和轉(zhuǎn)移路程最小(圖4(a)、(b));當(dāng)收割效率小于0.40 hm2/h 時(shí),2種模式下收割機(jī)等待時(shí)間與運(yùn)糧車轉(zhuǎn)移距離均可以達(dá)到最優(yōu)(圖4(c)),使得非生產(chǎn)性成本最低。收割效率超過0.53 hm2/h時(shí),最優(yōu)等待時(shí)間逐漸增加,原因是隨著收割效率增加,收割作業(yè)完成較快,留給運(yùn)糧車自由支配的轉(zhuǎn)移時(shí)間相對(duì)減少,因此總體上收割機(jī)等待時(shí)間增加。運(yùn)糧車總轉(zhuǎn)移距離在0.40 hm2/h效率時(shí)開始增加,達(dá)到0.53 hm2/h以后,最優(yōu)轉(zhuǎn)移距離保持不變,原因是收割效率增加,但模型受運(yùn)糧車數(shù)量制約,為降低收割機(jī)少量等待時(shí)間而造成運(yùn)糧車轉(zhuǎn)移距離顯著增大,不符合最優(yōu)非生產(chǎn)性成本的要求,即收割效率超過0.53 hm2/h,非生產(chǎn)性作業(yè)成本最優(yōu)結(jié)果取決于轉(zhuǎn)移距離是否最優(yōu),不再是等待時(shí)間。

    圖4 不同模型收割效率變化對(duì)目標(biāo)優(yōu)化結(jié)果的影響Fig.4 Influence of harvesting efficiency change on the target optimization results under different models

    實(shí)際收獲作業(yè)中,將收割效率控制在0.53 hm2/h以內(nèi),使得該配置下的收運(yùn)響應(yīng)效率最高,此范圍內(nèi)收割機(jī)等待時(shí)間最短,運(yùn)糧車轉(zhuǎn)移距離產(chǎn)生的影響最小。若考慮完工時(shí)間,選擇0.53 hm2/h收割效率,能達(dá)到較佳作業(yè)狀態(tài);若不考慮完工時(shí)間,選擇0.40 hm2/h收割效率,達(dá)到最優(yōu)作業(yè)狀態(tài)。

    4 討論與結(jié)論

    實(shí)際收獲作業(yè)過程中收割機(jī)與運(yùn)糧車響應(yīng)方式較多,如模式一:收、運(yùn)響應(yīng)過程中,收割機(jī)先裝滿糧倉,再完全卸出收取的糧食。模式二:收、運(yùn)響應(yīng)過程中,收割機(jī)裝滿糧倉,隨后不用全部卸出,攜帶部分糧食進(jìn)行下次作業(yè)。模式三:收、運(yùn)響應(yīng)過程中,收割機(jī)不用裝滿糧倉,但必須完全卸出所收糧食。模式四:收、運(yùn)響應(yīng)過程中,收割機(jī)不必裝滿糧倉,隨后也不必完全卸出所收糧食。模式二適合于運(yùn)糧車較少的場(chǎng)景;模式三適用于運(yùn)糧車量冗余的場(chǎng)景;模式四中收割機(jī)能夠自由根據(jù)運(yùn)糧車剩余容量收割或卸載糧食,在運(yùn)糧車抵達(dá)時(shí)直接卸糧,無論剩余直接返回作業(yè),使運(yùn)糧車和收割機(jī)減少等待時(shí)間、運(yùn)糧車載量充分,達(dá)到了收割機(jī)與運(yùn)糧車最佳狀態(tài)。但是實(shí)現(xiàn)這種復(fù)雜的動(dòng)態(tài)響應(yīng),全靠人工經(jīng)驗(yàn),難度較大,目前沒有合適模型和算法進(jìn)行實(shí)現(xiàn)。這3種模式共同點(diǎn)是收割機(jī)可以自由根據(jù)實(shí)際約束改變每次作業(yè)的收割量,在已知產(chǎn)量的作業(yè)區(qū)域中收割次數(shù)動(dòng)態(tài)變化,增大了模型構(gòu)建和解算難度。相比之下模式一更適合于規(guī)模數(shù)量稍大的收獲任務(wù),在1個(gè)任務(wù)單元內(nèi)收割機(jī)收割次數(shù)固定,為模型和算法構(gòu)建提供了便利,同時(shí)也比較符合實(shí)際生產(chǎn)場(chǎng)景。因此本研究以模式一構(gòu)建收、運(yùn)響應(yīng)模型,定義任務(wù)單元,憑借單元產(chǎn)量獲得收割次數(shù),并根據(jù)收割頻率劃分時(shí)間槽,設(shè)計(jì)算法進(jìn)行求解。

    本研究提出基于任務(wù)單元的多目標(biāo)規(guī)劃模型,共同考慮收割機(jī)田間等待時(shí)間及運(yùn)糧車距離遠(yuǎn)近的目標(biāo),以動(dòng)態(tài)規(guī)劃的思想劃分時(shí)間槽,根據(jù)目標(biāo)不同優(yōu)先級(jí)順序設(shè)計(jì)了算法 A 和算法 B,選擇非生產(chǎn)性作業(yè)成本、運(yùn)糧車調(diào)配路程和收割機(jī)非生產(chǎn)性等待時(shí)間等目標(biāo)優(yōu)化結(jié)果進(jìn)行對(duì)比發(fā)現(xiàn),優(yōu)先考慮等待時(shí)間的算法 A 在非生產(chǎn)性總成本和等待時(shí)間等方面,均優(yōu)于優(yōu)先考慮運(yùn)糧車距離的調(diào)度算法 B,因此采用算法 A 作為該調(diào)度模型的求解算法。該算法為解決運(yùn)糧車不足和提高資源配置效率提供了一種解決方案。同時(shí),具體探究了特定型號(hào)(收、運(yùn)容量分別為1.6和5.0 t;收、運(yùn)數(shù)量配比為3∶2;收割效率為0.53 hm2/h)收割機(jī)與運(yùn)糧車響應(yīng)調(diào)度,獲得該場(chǎng)景下任務(wù)單元響應(yīng)甘特圖與運(yùn)糧車精確路徑順序(圖2,表3)。

    運(yùn)糧車資源不足的情況下,運(yùn)糧車數(shù)量與收割機(jī)數(shù)量相差顯著時(shí),運(yùn)糧車轉(zhuǎn)移距離和收割機(jī)等待時(shí)間顯著增加。兩者數(shù)量相近時(shí),運(yùn)糧車數(shù)量變化產(chǎn)生的總成本基本不變,增益減小,與配置6 輛運(yùn)糧車相比,服務(wù)組織配置4 輛運(yùn)糧車也能取得較好響應(yīng)結(jié)果,同時(shí)節(jié)約運(yùn)糧車設(shè)備資源。進(jìn)一步研究發(fā)現(xiàn),收、運(yùn)數(shù)量配比為3∶2情況下,應(yīng)保持0.40 hm2/h,獲得該配置下最優(yōu)響應(yīng)作業(yè)效率,使等待時(shí)間最短、轉(zhuǎn)移距離最小、非生產(chǎn)性成本最低。對(duì)于其他型號(hào)的收割機(jī)和運(yùn)糧車,若提供收、運(yùn)的容量以及數(shù)量配比等參數(shù),能得出對(duì)應(yīng)條件下的具體結(jié)論,獲得最優(yōu)響應(yīng)方案和作業(yè)效率。因此該模型和算法能夠?yàn)檫\(yùn)糧車不足問題提供一定的指導(dǎo)意義。

    本研究模型創(chuàng)新之處在于引入了收、運(yùn)響應(yīng)機(jī)制,更符合田間實(shí)際作業(yè)情況。但由于運(yùn)糧車到達(dá)時(shí)間、載量不斷變化,設(shè)計(jì)算法在處理更大規(guī)模數(shù)據(jù)時(shí)無法收斂,適用算法需進(jìn)一步研究;同時(shí)模型僅考慮運(yùn)糧車,未涉及收割機(jī)精確調(diào)度規(guī)劃,針對(duì)收割機(jī)與運(yùn)糧車的兩級(jí)調(diào)度研究和增加對(duì)任務(wù)單元改變的動(dòng)態(tài)規(guī)劃研究是未來智能農(nóng)機(jī)調(diào)度領(lǐng)域的重要研究方向。

    猜你喜歡
    運(yùn)糧非生產(chǎn)性等待時(shí)間
    嘉興市2016—2021 年農(nóng)藥中毒病例分析
    給學(xué)生適宜的等待時(shí)間
    ——國外課堂互動(dòng)等待時(shí)間研究的現(xiàn)狀與啟示
    本期導(dǎo)讀
    烏克蘭稱土耳其扣留俄運(yùn)糧船
    非生產(chǎn)性支出損害還是促進(jìn)了企業(yè)創(chuàng)新?
    ——直接效應(yīng)與溢出效應(yīng)的檢驗(yàn)
    非生產(chǎn)性投入改進(jìn)了新創(chuàng)企業(yè)的經(jīng)營效率嗎?
    《人民日?qǐng)?bào)》猛批“濮陽運(yùn)糧事件”
    黨史縱覽(2017年8期)2017-09-13 19:01:26
    意大利:反腐敗沒有等待時(shí)間
    公民與法治(2016年2期)2016-05-17 04:08:28
    顧客等待心理的十條原則
    視野(2015年14期)2015-07-28 00:01:44
    顧客等待心理的十條原則
    讀者(2015年12期)2015-06-19 16:09:14
    一区二区三区免费毛片| 亚洲av免费在线观看| 国产精品免费一区二区三区在线| 观看美女的网站| 国产美女午夜福利| 波多野结衣巨乳人妻| 亚洲电影在线观看av| 国产一区二区在线观看日韩| 久久6这里有精品| 精品无人区乱码1区二区| 国产色婷婷99| 夜夜看夜夜爽夜夜摸| 久久这里只有精品中国| av在线天堂中文字幕| 国产精品一区www在线观看| 亚洲图色成人| 国产老妇女一区| 男女下面进入的视频免费午夜| av在线天堂中文字幕| 久久精品夜色国产| 美女被艹到高潮喷水动态| 永久网站在线| 国产av不卡久久| 在线观看午夜福利视频| 老女人水多毛片| 亚洲精品一卡2卡三卡4卡5卡| 97人妻精品一区二区三区麻豆| 国产成人福利小说| 别揉我奶头~嗯~啊~动态视频| 观看免费一级毛片| 午夜福利在线观看免费完整高清在 | 青春草视频在线免费观看| 亚洲人成网站在线播放欧美日韩| .国产精品久久| 又粗又爽又猛毛片免费看| 亚洲自偷自拍三级| 久久午夜亚洲精品久久| 国产精品野战在线观看| 亚洲精品在线观看二区| 免费观看人在逋| 国产视频内射| 欧美绝顶高潮抽搐喷水| avwww免费| 狠狠狠狠99中文字幕| 亚洲一区高清亚洲精品| av在线亚洲专区| 欧美又色又爽又黄视频| 欧美日韩乱码在线| 校园春色视频在线观看| 在线a可以看的网站| 一区福利在线观看| 亚洲高清免费不卡视频| 尤物成人国产欧美一区二区三区| h日本视频在线播放| 亚洲精品成人久久久久久| 亚洲av一区综合| 性欧美人与动物交配| 久久久久免费精品人妻一区二区| 欧美日韩综合久久久久久| 一区二区三区高清视频在线| 国产一区二区三区av在线 | 欧美日韩综合久久久久久| 日本黄大片高清| 亚洲精品一区av在线观看| 男女下面进入的视频免费午夜| 99在线人妻在线中文字幕| 一个人免费在线观看电影| 亚洲高清免费不卡视频| 午夜爱爱视频在线播放| 午夜福利视频1000在线观看| 久久久久性生活片| av黄色大香蕉| 久久久欧美国产精品| 亚洲av二区三区四区| 久久精品人妻少妇| 久久这里只有精品中国| 久久久久久大精品| 国产精品久久久久久av不卡| 国产黄色视频一区二区在线观看 | 中文字幕熟女人妻在线| 性插视频无遮挡在线免费观看| 亚洲欧美精品综合久久99| 两性午夜刺激爽爽歪歪视频在线观看| 日本 av在线| 少妇猛男粗大的猛烈进出视频 | 女同久久另类99精品国产91| 国产高清视频在线观看网站| 秋霞在线观看毛片| 看黄色毛片网站| 中文字幕精品亚洲无线码一区| 亚洲人成网站在线播| 黄色配什么色好看| 我的女老师完整版在线观看| 99riav亚洲国产免费| 干丝袜人妻中文字幕| 成人精品一区二区免费| 亚洲欧美日韩无卡精品| 日本免费a在线| 又爽又黄无遮挡网站| 国产免费男女视频| 香蕉av资源在线| 一本一本综合久久| 伦精品一区二区三区| 深夜a级毛片| 最近中文字幕高清免费大全6| 亚洲久久久久久中文字幕| 网址你懂的国产日韩在线| 久久精品国产清高在天天线| 可以在线观看的亚洲视频| 亚洲中文日韩欧美视频| 婷婷色综合大香蕉| 插阴视频在线观看视频| 亚洲丝袜综合中文字幕| 国产男靠女视频免费网站| 在线观看美女被高潮喷水网站| 午夜精品国产一区二区电影 | 你懂的网址亚洲精品在线观看 | 我的女老师完整版在线观看| 亚洲va在线va天堂va国产| 十八禁国产超污无遮挡网站| 麻豆久久精品国产亚洲av| 精品久久国产蜜桃| 性欧美人与动物交配| 看十八女毛片水多多多| 99热网站在线观看| 国产成人aa在线观看| 日韩欧美在线乱码| www.色视频.com| 好男人在线观看高清免费视频| 国产精品永久免费网站| 12—13女人毛片做爰片一| 欧美性猛交黑人性爽| 成人二区视频| 国产一级毛片七仙女欲春2| 亚洲熟妇熟女久久| 免费看av在线观看网站| 国产熟女欧美一区二区| 激情 狠狠 欧美| 日本黄色片子视频| 在线a可以看的网站| 精品午夜福利在线看| 热99re8久久精品国产| 尾随美女入室| 免费av毛片视频| 一级毛片我不卡| 啦啦啦观看免费观看视频高清| 悠悠久久av| 乱人视频在线观看| 日韩欧美国产在线观看| 人人妻,人人澡人人爽秒播| 一区二区三区四区激情视频 | 精品一区二区免费观看| 国产精品爽爽va在线观看网站| 极品教师在线视频| 亚洲中文字幕一区二区三区有码在线看| a级毛片免费高清观看在线播放| 99久国产av精品| 女同久久另类99精品国产91| 亚洲国产欧美人成| 成年av动漫网址| 国产极品精品免费视频能看的| 长腿黑丝高跟| 精品日产1卡2卡| 人人妻人人澡人人爽人人夜夜 | 国产精品久久电影中文字幕| 亚洲中文字幕日韩| 俄罗斯特黄特色一大片| 一卡2卡三卡四卡精品乱码亚洲| 免费av不卡在线播放| 亚洲18禁久久av| 51国产日韩欧美| 欧美激情国产日韩精品一区| 天堂动漫精品| 久久精品国产亚洲网站| 国产老妇女一区| 久久久精品94久久精品| 亚洲av电影不卡..在线观看| 悠悠久久av| 免费看美女性在线毛片视频| 久久久a久久爽久久v久久| 国产精品综合久久久久久久免费| 亚洲av第一区精品v没综合| 色尼玛亚洲综合影院| 永久网站在线| 女人十人毛片免费观看3o分钟| 亚洲不卡免费看| 永久网站在线| 99久久无色码亚洲精品果冻| av专区在线播放| 日韩欧美精品v在线| 国产精品久久视频播放| 国产伦一二天堂av在线观看| 亚洲av美国av| 久久久国产成人精品二区| 亚洲经典国产精华液单| 九九热线精品视视频播放| 免费看日本二区| 亚洲专区国产一区二区| 国产精品女同一区二区软件| 最近最新中文字幕大全电影3| 嫩草影视91久久| 毛片一级片免费看久久久久| 中文资源天堂在线| 亚洲一级一片aⅴ在线观看| 非洲黑人性xxxx精品又粗又长| 人妻制服诱惑在线中文字幕| 床上黄色一级片| 欧美3d第一页| 三级国产精品欧美在线观看| 亚洲高清免费不卡视频| 精品少妇黑人巨大在线播放 | 亚洲怡红院男人天堂| 乱系列少妇在线播放| av有码第一页| 多毛熟女@视频| 在线亚洲精品国产二区图片欧美 | 色哟哟·www| 国产精品久久久久久久久免| 22中文网久久字幕| 欧美 亚洲 国产 日韩一| 亚洲精华国产精华液的使用体验| 亚州av有码| 天堂俺去俺来也www色官网| 成人亚洲欧美一区二区av| 国产亚洲午夜精品一区二区久久| av线在线观看网站| 久久99精品国语久久久| 人妻少妇偷人精品九色| 插逼视频在线观看| 亚洲色图综合在线观看| 成人午夜精彩视频在线观看| 国产精品国产三级国产av玫瑰| 青春草国产在线视频| 欧美区成人在线视频| 伦理电影免费视频| 精品久久久噜噜| 亚洲熟女精品中文字幕| 国产黄色视频一区二区在线观看| 国产成人精品久久久久久| 国产日韩一区二区三区精品不卡 | 亚洲丝袜综合中文字幕| 日日撸夜夜添| 中文字幕亚洲精品专区| 一区二区三区精品91| 一区二区三区免费毛片| 国产一区二区三区综合在线观看 | 日韩大片免费观看网站| 纯流量卡能插随身wifi吗| 三级国产精品欧美在线观看| 一区二区三区免费毛片| 岛国毛片在线播放| 国产男女内射视频| 一级爰片在线观看| 亚洲怡红院男人天堂| 建设人人有责人人尽责人人享有的| 在线天堂最新版资源| 丰满迷人的少妇在线观看| 国产高清三级在线| 91aial.com中文字幕在线观看| 最近手机中文字幕大全| 伊人久久精品亚洲午夜| 久久午夜综合久久蜜桃| 久久狼人影院| 日韩不卡一区二区三区视频在线| 另类精品久久| 蜜臀久久99精品久久宅男| 一区二区三区免费毛片| 九草在线视频观看| 国产高清国产精品国产三级| 精品一区二区三卡| 99久久中文字幕三级久久日本| 国产乱来视频区| 欧美三级亚洲精品| 国模一区二区三区四区视频| 国产伦在线观看视频一区| 五月玫瑰六月丁香| 国产高清三级在线| 免费播放大片免费观看视频在线观看| 成人二区视频| 欧美亚洲 丝袜 人妻 在线| 99久国产av精品国产电影| 久久鲁丝午夜福利片| 国产亚洲午夜精品一区二区久久| 精品一区二区三卡| 在线观看三级黄色| 成年av动漫网址| 免费高清在线观看视频在线观看| 中文乱码字字幕精品一区二区三区| 女性生殖器流出的白浆| 一级毛片久久久久久久久女| 久久毛片免费看一区二区三区| 久久久久视频综合| 国产色爽女视频免费观看| 国产极品粉嫩免费观看在线 | 新久久久久国产一级毛片| 亚洲熟女精品中文字幕| 色视频www国产| 永久网站在线| 免费观看的影片在线观看| 亚洲欧美日韩东京热| 有码 亚洲区| 伊人久久国产一区二区| 亚洲av成人精品一二三区| 一区在线观看完整版| 国产美女午夜福利| 欧美三级亚洲精品| 亚洲va在线va天堂va国产| 岛国毛片在线播放| 国精品久久久久久国模美| 国产成人午夜福利电影在线观看| 丰满乱子伦码专区| 成年女人在线观看亚洲视频| 国产男人的电影天堂91| 老女人水多毛片| 搡女人真爽免费视频火全软件| 欧美3d第一页| 国产精品无大码| 人妻少妇偷人精品九色| 男的添女的下面高潮视频| 亚洲va在线va天堂va国产| 国产精品.久久久| 成人漫画全彩无遮挡| 日韩欧美一区视频在线观看 | 成人美女网站在线观看视频| 我的女老师完整版在线观看| 高清毛片免费看| 在线观看免费日韩欧美大片 | 国产在线男女| 久久久国产一区二区| 丁香六月天网| 亚洲av.av天堂| 欧美成人午夜免费资源| 99久久精品热视频| 亚洲国产精品国产精品| 少妇人妻久久综合中文| 九九爱精品视频在线观看| 欧美日韩在线观看h| 51国产日韩欧美| av网站免费在线观看视频| 午夜福利影视在线免费观看| 熟女电影av网| 最后的刺客免费高清国语| av有码第一页| 卡戴珊不雅视频在线播放| 国产在线视频一区二区| 少妇裸体淫交视频免费看高清| 18+在线观看网站| 一级a做视频免费观看| 亚洲三级黄色毛片| 日韩av免费高清视频| 老女人水多毛片| 精品一区二区免费观看| 国产亚洲av片在线观看秒播厂| 欧美bdsm另类| 狠狠精品人妻久久久久久综合| 黄片无遮挡物在线观看| 99热这里只有是精品在线观看| 一级,二级,三级黄色视频| av女优亚洲男人天堂| 国产av码专区亚洲av| 国精品久久久久久国模美| 欧美激情国产日韩精品一区| 91在线精品国自产拍蜜月| 久久狼人影院| 精华霜和精华液先用哪个| 欧美日韩国产mv在线观看视频| 亚洲精品一二三| 五月玫瑰六月丁香| 妹子高潮喷水视频| 国产成人精品无人区| 性色avwww在线观看| 亚洲国产日韩一区二区| 亚洲欧洲精品一区二区精品久久久 | 亚洲经典国产精华液单| 在线观看免费日韩欧美大片 | av不卡在线播放| 少妇的逼水好多| 国产精品秋霞免费鲁丝片| 人妻人人澡人人爽人人| 色5月婷婷丁香| 国产欧美另类精品又又久久亚洲欧美| 亚洲精品第二区| 久久久久久久久久久免费av| 777米奇影视久久| 欧美丝袜亚洲另类| 国产白丝娇喘喷水9色精品| 一级毛片电影观看| 一个人免费看片子| 精品久久久久久电影网| 国产精品一二三区在线看| 边亲边吃奶的免费视频| 婷婷色综合www| 亚洲精品一二三| 国产国拍精品亚洲av在线观看| 亚洲精品一区蜜桃| 久久99蜜桃精品久久| 黄色视频在线播放观看不卡| 国产精品.久久久| 国产日韩一区二区三区精品不卡 | 亚洲av综合色区一区| 亚洲第一区二区三区不卡| 美女xxoo啪啪120秒动态图| 欧美精品高潮呻吟av久久| 国产真实伦视频高清在线观看| 精品一区二区三区视频在线| 免费看光身美女| 免费久久久久久久精品成人欧美视频 | 香蕉精品网在线| 天天躁夜夜躁狠狠久久av| 国产精品久久久久久av不卡| 国产日韩欧美在线精品| 深夜a级毛片| 亚洲av免费高清在线观看| 精品少妇黑人巨大在线播放| xxx大片免费视频| 亚洲人成网站在线播| freevideosex欧美| 国产精品久久久久成人av| 亚洲,欧美,日韩| 色哟哟·www| 精品人妻一区二区三区麻豆| 只有这里有精品99| 国产精品久久久久久精品电影小说| 亚洲av二区三区四区| 啦啦啦中文免费视频观看日本| 久久久久久久久久成人| 亚洲欧洲日产国产| 精品人妻一区二区三区麻豆| 成人黄色视频免费在线看| 最近2019中文字幕mv第一页| 成年美女黄网站色视频大全免费 | 国产av国产精品国产| 97超碰精品成人国产| 久久久久精品久久久久真实原创| 中文乱码字字幕精品一区二区三区| 亚洲av二区三区四区| 18禁在线播放成人免费| 青春草视频在线免费观看| 久久精品国产自在天天线| 国产精品国产三级国产专区5o| 高清午夜精品一区二区三区| 亚洲情色 制服丝袜| 青春草视频在线免费观看| 国产成人免费观看mmmm| 九九爱精品视频在线观看| 日本与韩国留学比较| 国语对白做爰xxxⅹ性视频网站| 赤兔流量卡办理| 亚洲国产精品专区欧美| 亚洲精品456在线播放app| a级毛片在线看网站| 大香蕉97超碰在线| 国产在线视频一区二区| a级毛片免费高清观看在线播放| 伊人久久国产一区二区| 黄片无遮挡物在线观看| 欧美日韩av久久| 91精品国产九色| 亚洲av日韩在线播放| 久久国产亚洲av麻豆专区| 久久av网站| 狂野欧美激情性xxxx在线观看| av在线老鸭窝| 国内揄拍国产精品人妻在线| 国产视频内射| 男的添女的下面高潮视频| 国产老妇伦熟女老妇高清| 黄色日韩在线| 日本av手机在线免费观看| 80岁老熟妇乱子伦牲交| 国产男女内射视频| 在线观看一区二区三区激情| av播播在线观看一区| 久久精品夜色国产| 亚洲国产精品一区二区三区在线| 亚洲精品久久午夜乱码| 午夜福利,免费看| 精品久久国产蜜桃| 欧美一级a爱片免费观看看| 免费观看a级毛片全部| 一级黄片播放器| 两个人的视频大全免费| 欧美 亚洲 国产 日韩一| av天堂中文字幕网| 免费观看性生交大片5| 麻豆成人av视频| 在线观看国产h片| 日韩一区二区视频免费看| 免费播放大片免费观看视频在线观看| 嘟嘟电影网在线观看| 亚洲美女黄色视频免费看| a级毛色黄片| a级毛片在线看网站| 亚洲国产精品成人久久小说| 永久免费av网站大全| 日韩免费高清中文字幕av| 看免费成人av毛片| 久久 成人 亚洲| 不卡视频在线观看欧美| 在线观看国产h片| 国产成人午夜福利电影在线观看| 精品久久国产蜜桃| 久久狼人影院| 九色成人免费人妻av| 男人添女人高潮全过程视频| 亚洲精品色激情综合| 亚洲av成人精品一区久久| 九九久久精品国产亚洲av麻豆| 中文字幕免费在线视频6| 日本午夜av视频| 妹子高潮喷水视频| freevideosex欧美| 亚洲精品成人av观看孕妇| 久久久亚洲精品成人影院| 熟女人妻精品中文字幕| 高清黄色对白视频在线免费看 | 精品少妇黑人巨大在线播放| 一级毛片久久久久久久久女| 亚洲国产精品999| 成人漫画全彩无遮挡| 另类亚洲欧美激情| 精品一区在线观看国产| 国产伦精品一区二区三区四那| 又粗又硬又长又爽又黄的视频| 亚洲三级黄色毛片| 自拍欧美九色日韩亚洲蝌蚪91 | 成人国产av品久久久| 亚洲精品,欧美精品| 99热全是精品| 国产精品成人在线| 亚洲精品中文字幕在线视频 | 国产av国产精品国产| www.av在线官网国产| 啦啦啦视频在线资源免费观看| 少妇高潮的动态图| 嫩草影院入口| 黑丝袜美女国产一区| av不卡在线播放| 永久免费av网站大全| 国产成人aa在线观看| 亚洲精品久久久久久婷婷小说| 热re99久久精品国产66热6| av免费观看日本| 精品99又大又爽又粗少妇毛片| 九九久久精品国产亚洲av麻豆| 97在线人人人人妻| 桃花免费在线播放| av黄色大香蕉| 各种免费的搞黄视频| 成人综合一区亚洲| 国产男女内射视频| 午夜视频国产福利| 在线观看国产h片| 久久久久久久久久成人| 人人澡人人妻人| 欧美激情国产日韩精品一区| 免费观看性生交大片5| 成人免费观看视频高清| 美女内射精品一级片tv| 久久久久久久久久久免费av| 最后的刺客免费高清国语| 在线看a的网站| 成年人午夜在线观看视频| 纵有疾风起免费观看全集完整版| 国产淫片久久久久久久久| 精品人妻偷拍中文字幕| 妹子高潮喷水视频| 看非洲黑人一级黄片| 亚洲国产精品国产精品| av免费在线看不卡| 在线观看一区二区三区激情| 精品久久国产蜜桃| 久久久欧美国产精品| 亚洲情色 制服丝袜| 精品人妻熟女毛片av久久网站| 91在线精品国自产拍蜜月| 看十八女毛片水多多多| 王馨瑶露胸无遮挡在线观看| av.在线天堂| 亚洲av免费高清在线观看| 美女大奶头黄色视频| 三级国产精品片| 精品国产乱码久久久久久小说| 亚洲人成网站在线播| 成人国产av品久久久| 美女大奶头黄色视频| 亚洲欧美日韩另类电影网站| 国产精品无大码| 自线自在国产av| 一本色道久久久久久精品综合| 日韩电影二区| 成人18禁高潮啪啪吃奶动态图 | 乱人伦中国视频| 97超视频在线观看视频| 国产淫语在线视频| av一本久久久久| 免费看光身美女| 日韩欧美精品免费久久| 美女脱内裤让男人舔精品视频| 日本欧美国产在线视频| 日本黄大片高清| 国产美女午夜福利| 久久99精品国语久久久| 国产精品蜜桃在线观看| 亚洲国产av新网站| 观看免费一级毛片| 国产国拍精品亚洲av在线观看| 18+在线观看网站| av在线播放精品| 欧美日韩综合久久久久久| 久久久久网色| av在线观看视频网站免费| 成人二区视频| 久久久久久久久大av| 黄色欧美视频在线观看|