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

    基于改進(jìn)螢火蟲(chóng)算法的礦井水害避災(zāi)路徑規(guī)劃

    2021-06-18 07:04:52朱希安王占剛劉德民
    中國(guó)礦業(yè) 2021年6期
    關(guān)鍵詞:當(dāng)量螢火蟲(chóng)亮度

    王 鵬,朱希安,王占剛,劉德民

    (1.北京信息科技大學(xué)信息與通信工程學(xué)院,北京 100101;2.華北科技學(xué)院安全工程學(xué)院,北京 101601)

    在煤礦生產(chǎn)過(guò)程中,從業(yè)人員在巷道中隨時(shí)面臨著礦井突水等自然災(zāi)害的威脅。當(dāng)突水災(zāi)害發(fā)生時(shí),如何使受困人員盡快擺脫險(xiǎn)境,是一個(gè)值得考慮的問(wèn)題。地下井巷道交錯(cuò)復(fù)雜,從事故發(fā)生地到安全區(qū)域通常會(huì)有很多條路徑,在諸多路徑中選擇一條最優(yōu)路徑符合從業(yè)人員的安全利益[1-2]。

    有學(xué)者對(duì)礦井突水逃生路徑規(guī)劃問(wèn)題進(jìn)行了研究。苑亞南等[3]通過(guò)綜合考慮突水點(diǎn)處流體力學(xué)特性、巷道基本信息,提出了基于改進(jìn)的D-K路徑搜索算法的多路徑逃生模型;蔡明杰等[4]通過(guò)重新設(shè)計(jì)遺傳算子,引入新的交叉、變異算子,較好地提高了遺傳算法在規(guī)劃路徑時(shí)的收斂速度和可靠性;劉夢(mèng)杰等[5]通過(guò)將傳統(tǒng)A*算法單向搜索方式改為雙向搜索并改進(jìn)評(píng)估函數(shù)的定義,使得每一步搜索都根據(jù)評(píng)估函數(shù)朝著某個(gè)方向搜索。隨著計(jì)算機(jī)運(yùn)算能力的加強(qiáng),蟻群算法等元啟發(fā)式算法逐漸成為路徑規(guī)劃的主流算法。但是,上述算法都存在以下缺點(diǎn)中的一個(gè)或多個(gè):易陷入局部最優(yōu)、收斂效果不佳、無(wú)效計(jì)算過(guò)多等。鑒于上述背景,本文通過(guò)引入巷道當(dāng)量長(zhǎng)度以調(diào)整螢火蟲(chóng)的初始化方式,改變螢火蟲(chóng)的間距、最大亮度和相對(duì)亮度等要素的定義,結(jié)合王家?guī)X礦實(shí)例實(shí)現(xiàn)對(duì)避災(zāi)路徑的最優(yōu)規(guī)劃。

    1 巷道當(dāng)量長(zhǎng)度

    考慮到受困人員的身體狀況以及每條巷道的客觀干擾因素,諸如風(fēng)速風(fēng)向、巷道坡度、局部障礙物類型、有毒有害氣體濃度等都會(huì)對(duì)逃生巷道的安全性產(chǎn)生影響,所以不能簡(jiǎn)單地求解線路的最短距離并以此充當(dāng)最優(yōu)路徑。因此,本文將當(dāng)量長(zhǎng)度的概念引入避災(zāi)路徑的規(guī)劃中,通過(guò)綜合評(píng)價(jià)巷道中對(duì)逃生安全有重要影響的因素,計(jì)算巷道的當(dāng)量長(zhǎng)度[6],并將該值作為巷道網(wǎng)絡(luò)結(jié)構(gòu)的權(quán)重。

    1.1 影響系數(shù)的定義

    通過(guò)對(duì)巷道實(shí)際情況的研判,本文重點(diǎn)的干擾因素為風(fēng)速(含風(fēng)向)、巷道坡度、局部障礙物類型、巷道類型,對(duì)應(yīng)的影響系數(shù)[7]分別為δ1、δ2、δ3、δ4,干擾因素的影響系數(shù)見(jiàn)表1。 設(shè)不受某項(xiàng)因素影響時(shí)通過(guò)巷道Rij的速度為V0,受某項(xiàng)因素影響時(shí)通過(guò)巷道Rij的速度為V,則影響系數(shù)的表達(dá)式為式(1)。

    表1 巷道干擾因素及影響系數(shù)Table 1 Disturbance factors and influence coefficients of roadway

    (1)

    1.2 巷道當(dāng)量長(zhǎng)度的求解

    將各巷道影響系數(shù)加權(quán)到對(duì)應(yīng)的巷道實(shí)際長(zhǎng)度即可得巷道當(dāng)量長(zhǎng)度。設(shè)定節(jié)點(diǎn)i和節(jié)點(diǎn)j(i,j須為可直接相連的兩點(diǎn))之間的巷道實(shí)際長(zhǎng)度為lij,δ1,δ2,δ3,…,δn為該巷道干擾因素的影響系數(shù),巷道Rij的當(dāng)量長(zhǎng)度wij定義見(jiàn)式(2)。

    wij=[1+δ1+δ2+δ3+…+δn]lij

    (2)

    若有一條從突水點(diǎn)至安全點(diǎn)的路徑U,包含p段巷道,第q條巷道的當(dāng)量長(zhǎng)度用wq表示,則路徑U的當(dāng)量長(zhǎng)度WU可表示為式(3)。

    (3)

    當(dāng)量長(zhǎng)度表示通過(guò)該逃生路徑的難易程度,路徑的當(dāng)量長(zhǎng)度值越大,代表該路徑越難以通過(guò);路徑的當(dāng)量長(zhǎng)度值越小,代表該路徑越容易通過(guò)??蓪⒏倪M(jìn)后的螢火蟲(chóng)算法中的目標(biāo)函數(shù)用當(dāng)量長(zhǎng)度表示,巷道當(dāng)量長(zhǎng)度越小則目標(biāo)函數(shù)越優(yōu),反之則越差。

    2 標(biāo)準(zhǔn)螢火蟲(chóng)算法

    2.1 標(biāo)準(zhǔn)螢火蟲(chóng)算法基本原理

    螢火蟲(chóng)算法(firefly algorithm,F(xiàn)A)是一種基于螢火蟲(chóng)的發(fā)光模式和運(yùn)動(dòng)模式的元啟發(fā)式算法。其原理是螢火蟲(chóng)隨機(jī)地分布在特定的場(chǎng)域中,亮度較小的螢火蟲(chóng)會(huì)飛向亮度較大的螢火蟲(chóng),隨著移動(dòng)的過(guò)程在亮度較大的螢火蟲(chóng)處聚集。螢火蟲(chóng)算法是通過(guò)其運(yùn)動(dòng)、更新位置的過(guò)程來(lái)探尋最優(yōu)解的過(guò)程[8]。

    螢火蟲(chóng)相對(duì)亮度的表達(dá)式見(jiàn)式(4)。

    (4)

    式中:I0為螢火蟲(chóng)距離為零時(shí)的亮度,即最大亮度;γ為光吸收系數(shù);Rij為螢火蟲(chóng)i和j的間距。

    螢火蟲(chóng)吸引度表達(dá)式見(jiàn)式(5)。

    (5)

    式中:β0為兩只螢火蟲(chóng)距離為零時(shí)的吸引力大小,即最大吸引力;γ、Rij含義同式(4)。

    任意兩只螢火蟲(chóng)的間距通過(guò)笛卡爾距離定義,表達(dá)式見(jiàn)式(6)。

    (6)

    式中:xi、xj為螢火蟲(chóng)i、j在場(chǎng)域中的位置;d為解空間的總維度。

    螢火蟲(chóng)移動(dòng)后位置的表達(dá)式見(jiàn)式(7)。

    xi(t+1)=xi(t)+β(Rij)(xj(t)-xi(t))+αεi

    (7)

    式中:xi(t)、xj(t)為螢火蟲(chóng)i、j在第t次迭代所處位置;α為隨機(jī)化參數(shù),α∈[0,1];εi為搜索精度,εi∈[-0.5,0.5]。

    2.2 標(biāo)準(zhǔn)螢火蟲(chóng)算法描述

    在尋優(yōu)問(wèn)題中,目標(biāo)函數(shù)表示為螢火蟲(chóng)的亮度,亮度越大代表該螢火蟲(chóng)所處的空間解越好。螢火蟲(chóng)算法首先設(shè)定一些基本參數(shù),通過(guò)循環(huán)迭代螢火蟲(chóng)的亮度并更新螢火蟲(chóng)的亮度和位置,直至滿足終止條件[9]。標(biāo)準(zhǔn)螢火蟲(chóng)算法步驟如下所述。①設(shè)定基本參數(shù):β0、γ、α、N(螢火蟲(chóng)數(shù)量)、T(最大迭代次數(shù));②隨機(jī)初始化螢火蟲(chóng)位置,根據(jù)目標(biāo)函數(shù)計(jì)算其亮度;③計(jì)算螢火蟲(chóng)的亮度和吸引度,確定螢火蟲(chóng)的移動(dòng)方向;④確定螢火蟲(chóng)更新移動(dòng)后的位置,對(duì)其亮度進(jìn)行重新計(jì)算;⑤若達(dá)到最大迭代次數(shù)則結(jié)束算法,否則返回到第②步繼續(xù)搜索;⑥輸出結(jié)果,算法程序結(jié)束。標(biāo)準(zhǔn)螢火蟲(chóng)算法流程圖如圖1所示。

    圖1 標(biāo)準(zhǔn)螢火蟲(chóng)算法流程圖Fig.1 Flow chart of standard firefly algorithm

    標(biāo)準(zhǔn)螢火蟲(chóng)算法雖然具有參數(shù)少、操作性強(qiáng)等優(yōu)點(diǎn),但也具有易陷入局部最優(yōu)的缺點(diǎn)。螢火蟲(chóng)都是飛向比自己亮度更大的螢火蟲(chóng),但當(dāng)螢火蟲(chóng)個(gè)體全都集中于場(chǎng)域中的某一個(gè)或某幾個(gè)位置時(shí),其搜索能力顯著下降。本文通過(guò)改進(jìn)螢火蟲(chóng)基本要素的定義及搜索策略,將螢火蟲(chóng)算法應(yīng)用于礦井突水逃生路徑規(guī)劃的研究。

    3 改進(jìn)的螢火蟲(chóng)算法

    3.1 改進(jìn)的螢火蟲(chóng)算法描述

    設(shè)定在某場(chǎng)域中有N只螢火蟲(chóng),M個(gè)離散點(diǎn)(巷道節(jié)點(diǎn)),針對(duì)礦井逃生路徑,令每只螢火蟲(chóng)代表一條完整的從突水點(diǎn)至安全點(diǎn)的路徑。首先,對(duì)螢火蟲(chóng)進(jìn)行初始化可得N條逃生路徑,計(jì)算每條路徑的當(dāng)量長(zhǎng)度,保留當(dāng)量長(zhǎng)度最小的路徑,定義其為全局最優(yōu)路徑。此處,將螢火蟲(chóng)最大亮度定義為該條路徑當(dāng)量長(zhǎng)度的倒數(shù)。其次,比較螢火蟲(chóng)的最大亮度值和種群最大亮度平均值,舍棄最大亮度值低于種群平均亮度值的路徑,保留大于等于種群平均亮度值的路徑,對(duì)舍棄的路徑重新初始化。最后,計(jì)算螢火蟲(chóng)之間的相對(duì)亮度,以確定螢火蟲(chóng)的移動(dòng)方向并進(jìn)行移動(dòng)。將移動(dòng)后的路徑當(dāng)量長(zhǎng)度與當(dāng)前全局最優(yōu)路徑對(duì)比,若移動(dòng)后有路徑的當(dāng)量長(zhǎng)度小于當(dāng)前全局最優(yōu)路徑,則保留該路徑的當(dāng)量長(zhǎng)度并將其記為全局最優(yōu)路徑;反之則不變。通過(guò)數(shù)次迭代即可得到全局最優(yōu)路徑,該路徑可滿足礦井水害避災(zāi)路徑規(guī)劃的需求。

    3.2 改進(jìn)的螢火蟲(chóng)算法具體步驟

    1) 初始化螢火蟲(chóng)。如圖2所示,設(shè)定節(jié)點(diǎn)a0有n個(gè)與之直接相連的節(jié)點(diǎn)序列a={a1,a2,a3,…,an},本文通過(guò)輪盤賭結(jié)合路徑權(quán)重的方法初始化螢火蟲(chóng)。若ap為節(jié)點(diǎn)序列a中的一點(diǎn),則節(jié)點(diǎn)a0選擇ap為下一節(jié)點(diǎn)的概率計(jì)算見(jiàn)式(8)。

    圖2 螢火蟲(chóng)初始化的選擇概率Fig.2 Selection probability of firefly initialization

    (8)

    式中:k∈{1,2,…,n};w0p為節(jié)點(diǎn)a0與ap相連路徑的權(quán)重;w0k定義同理。若兩節(jié)點(diǎn)不相連,則兩節(jié)點(diǎn)互相選擇的概率均為0。

    為選擇一條從突水點(diǎn)至安全點(diǎn)的有效路徑,必須保證已選擇的點(diǎn)不會(huì)再被選擇,故建立相鄰點(diǎn)矩陣和移動(dòng)選擇矩陣。相鄰點(diǎn)矩陣通過(guò)設(shè)定任意兩點(diǎn)的關(guān)系為1或0來(lái)判斷是否直接相連,若節(jié)點(diǎn)a0和a1直接相連,則表示為Con(a0,a1)=1;若節(jié)點(diǎn)a0和b0不直接相連,則表示為Con(a0,b0)=0。移動(dòng)選擇矩陣中的值表示為任意一點(diǎn)選擇另一直接相連點(diǎn)的概率,若節(jié)點(diǎn)a0選擇到了a1為其下一節(jié)點(diǎn),則Con(a0,a1)=1,并修改a1對(duì)a0的移動(dòng)概率為P(a0,a1)=0,可保證節(jié)點(diǎn)a1不會(huì)再選擇到節(jié)點(diǎn)a0。以突水點(diǎn)為起點(diǎn)進(jìn)行初始化直至選擇到安全點(diǎn)為止,即可得到一條完整路徑。

    2) 舍棄劣質(zhì)螢火蟲(chóng)路徑。通過(guò)求得初始化后螢火蟲(chóng)路徑的最大亮度平均值,逐個(gè)將螢火蟲(chóng)的最大亮度值與最大亮度平均值進(jìn)行比較,當(dāng)某螢火蟲(chóng)亮度值大于或等于亮度平均值時(shí),則保留該條路徑;反之,則不予保留。并返回步驟1)重新初始化舍棄掉的螢火蟲(chóng)路徑。

    3) 最大亮度。由于螢火蟲(chóng)都是向比自己亮度大的螢火蟲(chóng)移動(dòng),所以將螢火蟲(chóng)的最大亮度定義為目標(biāo)函數(shù)的倒數(shù),即螢火蟲(chóng)將通過(guò)迭代移動(dòng)至目標(biāo)函數(shù)(當(dāng)量長(zhǎng)度值)比自己小的螢火蟲(chóng)個(gè)體。

    4) 相對(duì)亮度和吸引度。由于螢火蟲(chóng)序列中節(jié)點(diǎn)數(shù)可能不一致,所以對(duì)螢火蟲(chóng)的節(jié)點(diǎn)序列進(jìn)行補(bǔ)0操作。設(shè)某螢火蟲(chóng)i包含n個(gè)節(jié)點(diǎn),記為Xi={xi1,xi2,…,xin};某螢火蟲(chóng)j包含m個(gè)節(jié)點(diǎn),記為Xj={xj1,xj2,…,xjm},在螢火蟲(chóng)i序列的終點(diǎn)xin前補(bǔ)(M-n)個(gè)0,在螢火蟲(chóng)j序列的終點(diǎn)xjn前補(bǔ)(M-m)個(gè)0,使每個(gè)序列中節(jié)點(diǎn)個(gè)數(shù)均為節(jié)點(diǎn)總數(shù)。螢火蟲(chóng)i到螢火蟲(chóng)j的距離表示為式(9)。

    (9)

    式中,xik、xjk為螢火蟲(chóng)i、j的第k個(gè)節(jié)點(diǎn)。

    5) 移動(dòng)。從兩條路徑中選擇任意一個(gè)公共點(diǎn),將公共點(diǎn)后的一段路徑進(jìn)行替換以完成移動(dòng)。如螢火蟲(chóng)i的路徑序列為Xi={xi1,xi2,…,xiM};某螢火蟲(chóng)j的路徑序列為Xj={xj1,xj2,…,xjM},選擇一公共點(diǎn)(不考慮終點(diǎn))xia=xjb,若螢火蟲(chóng)i移動(dòng)至j,則將xjb后長(zhǎng)度小于等于μ的一段路徑取代xia后長(zhǎng)度小于等于μ的一段路徑,如果替換后螢火蟲(chóng)i仍為從起點(diǎn)到終點(diǎn)的完整路徑,則結(jié)束移動(dòng);反之,則隨機(jī)另找一公共點(diǎn)進(jìn)行上述操作,直至找到一條可行路線。若沒(méi)有找到可行路線,則重新移動(dòng)。移動(dòng)后對(duì)螢火蟲(chóng)的最大亮度、相對(duì)亮度、吸引度更新,若更新后的值小于全局最優(yōu)解的值,則將全局最優(yōu)解替換為更新后的值。

    6) 終止迭代。若達(dá)到最大迭代次數(shù),則終止循環(huán),輸出全局最優(yōu)解。反之,則返回步驟2)繼續(xù)循環(huán)。改進(jìn)的螢火蟲(chóng)算法流程圖如圖3所示。

    圖3 改進(jìn)的螢火蟲(chóng)算法流程圖Fig.3 Flow chart of improved firefly algorithm

    4 改進(jìn)的螢火蟲(chóng)算法在王家?guī)X礦的實(shí)現(xiàn)

    本文選取王家?guī)X部分巷道構(gòu)建了如圖4所示的無(wú)向網(wǎng)絡(luò)結(jié)構(gòu)圖,圖中的圓圈表示巷道節(jié)點(diǎn),線段表示巷道,圖中未等比例構(gòu)建。統(tǒng)計(jì)的巷道干擾因素、影響系數(shù)及巷道當(dāng)量長(zhǎng)度見(jiàn)表2。通過(guò)使用改進(jìn)的螢火蟲(chóng)算法和A*算法規(guī)劃逃生路徑,比較和分析兩個(gè)算法規(guī)劃的路徑結(jié)果。

    圖4 巷道簡(jiǎn)化架構(gòu)圖Fig.4 Simplified architecture diagram of roadway

    表2 王家?guī)X局部巷道實(shí)際長(zhǎng)度及干擾因素信息Table 2 Information on actual length and interference factors of local roadways in Wangjialing

    為驗(yàn)證算法在路徑規(guī)劃中的可行性,本文基于MATLAB環(huán)境進(jìn)行仿真。因算法采用輪盤賭的方式進(jìn)行節(jié)點(diǎn)選擇,路徑規(guī)劃結(jié)果會(huì)受到隨機(jī)因素的影響,為了減小隨機(jī)因素所帶來(lái)的差異,本文設(shè)置多組突水點(diǎn)和安全點(diǎn)以提高實(shí)驗(yàn)的可靠性。改進(jìn)的螢火蟲(chóng)算法參數(shù)取值為:巷道節(jié)點(diǎn)數(shù)M=16,螢火蟲(chóng)個(gè)數(shù)N=40,最大迭代次數(shù)T=40,光吸收系數(shù)γ=1.0,替換長(zhǎng)度數(shù)μ=600。改進(jìn)后的螢火蟲(chóng)算法和A*算法路徑規(guī)劃的測(cè)試結(jié)果見(jiàn)表3和表4。

    通過(guò)分析表3和表4的實(shí)驗(yàn)結(jié)果可知,改進(jìn)的螢火蟲(chóng)算法所規(guī)劃路徑的實(shí)際長(zhǎng)度和當(dāng)量長(zhǎng)度比A*算法求得的長(zhǎng)度都短,表明改進(jìn)的螢火蟲(chóng)算法更易找到最優(yōu)解。因?yàn)楸疚目紤]了巷道中的干擾因素,故本算法規(guī)劃的逃生路徑能更好地滿足受困人員的安全利益,在礦井事故應(yīng)急救援應(yīng)用中具有較大的價(jià)值。

    表3 優(yōu)化算法與A*算法規(guī)劃路徑線路對(duì)比Table 3 Comparison of optimization algorithm and A* algorithm for planning route routes

    表4 優(yōu)化算法與A*算法規(guī)劃路徑長(zhǎng)度對(duì)比Table 4 Comparison of optimization algorithm and A* algorithm planning path length

    5 結(jié) 論

    1) 提出了改進(jìn)的螢火蟲(chóng)算法,重新規(guī)劃螢火蟲(chóng)算法的初始化方式、螢火蟲(chóng)之間的距離等要素,改變搜索策略,改進(jìn)后的算法具有全局優(yōu)化能力強(qiáng)、穩(wěn)定性強(qiáng)、收斂效果佳等優(yōu)點(diǎn)。

    2) 考慮了影響巷道通行的諸多不利因素,并以當(dāng)量長(zhǎng)度的形式數(shù)值化地表示通過(guò)巷道的難易程度,該改進(jìn)算法可規(guī)劃出滿足受困人員安全利益的路徑。

    3) 通過(guò)實(shí)例對(duì)比,改進(jìn)的螢火蟲(chóng)算法相比于A*算法規(guī)劃出的路徑實(shí)際長(zhǎng)度和當(dāng)量長(zhǎng)度均有減少。即螢火蟲(chóng)算法可以應(yīng)用于礦井水害避災(zāi)路徑規(guī)劃的場(chǎng)景中,有效降低受困人員的傷亡率。

    猜你喜歡
    當(dāng)量螢火蟲(chóng)亮度
    亮度調(diào)色多面手
    螢火蟲(chóng)
    螢火蟲(chóng)
    亮度一樣嗎?
    基于斬波調(diào)制的LED亮度控制
    人生的亮度
    黃河之聲(2016年24期)2016-02-03 09:01:52
    抱抱就不哭了
    夏天的螢火蟲(chóng)
    超壓測(cè)試方法對(duì)炸藥TNT當(dāng)量計(jì)算結(jié)果的影響
    国产精品一区二区三区四区久久 | 每晚都被弄得嗷嗷叫到高潮| 97人妻天天添夜夜摸| 中出人妻视频一区二区| 中国美女看黄片| 亚洲熟妇熟女久久| 日韩免费av在线播放| 精品少妇一区二区三区视频日本电影| 淫妇啪啪啪对白视频| 人人妻人人澡人人看| 在线观看免费日韩欧美大片| 69av精品久久久久久| 久久99一区二区三区| 国产精品爽爽va在线观看网站 | 母亲3免费完整高清在线观看| 亚洲狠狠婷婷综合久久图片| 欧美日本中文国产一区发布| 丰满的人妻完整版| 国产亚洲精品第一综合不卡| 天天躁夜夜躁狠狠躁躁| 9191精品国产免费久久| 日日夜夜操网爽| 视频区图区小说| 国产熟女午夜一区二区三区| 亚洲五月天丁香| 亚洲七黄色美女视频| 丝袜美腿诱惑在线| 琪琪午夜伦伦电影理论片6080| av天堂在线播放| 欧美一区二区精品小视频在线| 日韩免费高清中文字幕av| 99国产精品99久久久久| 看片在线看免费视频| 久久青草综合色| 50天的宝宝边吃奶边哭怎么回事| 成人18禁在线播放| 黑人猛操日本美女一级片| 国产成人欧美在线观看| 性少妇av在线| 在线观看免费日韩欧美大片| 久久久久国内视频| 久久天躁狠狠躁夜夜2o2o| 男女之事视频高清在线观看| 久久人妻福利社区极品人妻图片| 久久香蕉国产精品| 欧美激情极品国产一区二区三区| 美女高潮喷水抽搐中文字幕| 亚洲视频免费观看视频| 在线永久观看黄色视频| 宅男免费午夜| 国产91精品成人一区二区三区| 另类亚洲欧美激情| 80岁老熟妇乱子伦牲交| 日本 av在线| av在线播放免费不卡| 国产又爽黄色视频| 久久国产精品男人的天堂亚洲| 亚洲国产毛片av蜜桃av| 无遮挡黄片免费观看| 无遮挡黄片免费观看| 午夜影院日韩av| 亚洲欧美日韩无卡精品| 成人三级黄色视频| 啦啦啦免费观看视频1| 国产成人系列免费观看| 在线视频色国产色| 亚洲国产毛片av蜜桃av| av欧美777| avwww免费| 黄片小视频在线播放| 88av欧美| 欧美激情极品国产一区二区三区| 女性生殖器流出的白浆| 亚洲av电影在线进入| 国产一区二区三区在线臀色熟女 | 淫妇啪啪啪对白视频| 欧美最黄视频在线播放免费 | 又紧又爽又黄一区二区| 在线观看66精品国产| 日韩中文字幕欧美一区二区| 亚洲中文字幕日韩| 女人被躁到高潮嗷嗷叫费观| 91字幕亚洲| 人人妻,人人澡人人爽秒播| 一级a爱片免费观看的视频| 午夜a级毛片| 18禁国产床啪视频网站| 黄片小视频在线播放| 亚洲精品美女久久久久99蜜臀| 亚洲精品美女久久久久99蜜臀| 亚洲精品粉嫩美女一区| av在线天堂中文字幕 | 精品熟女少妇八av免费久了| 日本一区二区免费在线视频| 国产麻豆69| 女人精品久久久久毛片| 亚洲第一欧美日韩一区二区三区| 久久久久国产一级毛片高清牌| 久久久久国产一级毛片高清牌| 久久伊人香网站| 黑人欧美特级aaaaaa片| 久久青草综合色| 在线国产一区二区在线| 色播在线永久视频| 欧美乱色亚洲激情| 曰老女人黄片| 在线观看免费视频日本深夜| 麻豆久久精品国产亚洲av | 国产精华一区二区三区| 成年人黄色毛片网站| 一进一出抽搐动态| 久久久久久亚洲精品国产蜜桃av| 人妻久久中文字幕网| 多毛熟女@视频| 久久人人爽av亚洲精品天堂| 国产1区2区3区精品| 亚洲欧美一区二区三区久久| 99国产精品一区二区三区| 97人妻天天添夜夜摸| 在线观看www视频免费| 女人被躁到高潮嗷嗷叫费观| 久久人人精品亚洲av| 成人三级黄色视频| av视频免费观看在线观看| 久久精品影院6| 色综合站精品国产| 久久久久久久精品吃奶| a在线观看视频网站| 99香蕉大伊视频| 亚洲欧美激情综合另类| 亚洲精品国产一区二区精华液| √禁漫天堂资源中文www| 精品第一国产精品| 91精品国产国语对白视频| 女人爽到高潮嗷嗷叫在线视频| 岛国在线观看网站| 在线观看免费视频日本深夜| 成人国产一区最新在线观看| 久久国产乱子伦精品免费另类| 精品一区二区三卡| 69av精品久久久久久| 成人黄色视频免费在线看| 日本撒尿小便嘘嘘汇集6| 欧美中文综合在线视频| 最近最新中文字幕大全电影3 | 侵犯人妻中文字幕一二三四区| 99久久久亚洲精品蜜臀av| 精品久久久久久,| 亚洲激情在线av| 成年人免费黄色播放视频| 亚洲性夜色夜夜综合| 亚洲欧美一区二区三区黑人| 免费看十八禁软件| 久久香蕉激情| 欧美中文日本在线观看视频| 成人影院久久| 国产熟女xx| 久热爱精品视频在线9| 国产99白浆流出| 欧美日本亚洲视频在线播放| 亚洲av美国av| 亚洲九九香蕉| 又大又爽又粗| 国产精品秋霞免费鲁丝片| 狂野欧美激情性xxxx| xxx96com| 91九色精品人成在线观看| 在线观看免费视频日本深夜| 欧美+亚洲+日韩+国产| 国产xxxxx性猛交| 国产一区二区三区综合在线观看| 日韩欧美免费精品| 国产精品偷伦视频观看了| 久久久久久亚洲精品国产蜜桃av| 9色porny在线观看| 9191精品国产免费久久| 久热爱精品视频在线9| 女人高潮潮喷娇喘18禁视频| 亚洲视频免费观看视频| 黑丝袜美女国产一区| 又大又爽又粗| 操美女的视频在线观看| 色尼玛亚洲综合影院| 午夜日韩欧美国产| 亚洲精品一卡2卡三卡4卡5卡| 亚洲精品在线观看二区| 午夜免费观看网址| 国内毛片毛片毛片毛片毛片| cao死你这个sao货| 国产成人av教育| 国产精品 欧美亚洲| 真人做人爱边吃奶动态| 在线观看免费视频日本深夜| 精品久久久久久,| 亚洲久久久国产精品| 日日干狠狠操夜夜爽| 午夜精品久久久久久毛片777| 一级片'在线观看视频| 中文字幕av电影在线播放| 成年版毛片免费区| 亚洲全国av大片| avwww免费| 国产一区二区三区在线臀色熟女 | 女性被躁到高潮视频| 嫩草影视91久久| 99热只有精品国产| 免费搜索国产男女视频| 久久久精品国产亚洲av高清涩受| 久久影院123| 亚洲专区字幕在线| 精品电影一区二区在线| 精品国产一区二区三区四区第35| 国产深夜福利视频在线观看| 国产高清激情床上av| 日本欧美视频一区| 亚洲av日韩精品久久久久久密| 久久久久久久精品吃奶| 免费在线观看亚洲国产| 久久热在线av| 欧美日韩av久久| 亚洲自拍偷在线| 亚洲精品国产精品久久久不卡| 亚洲一区二区三区欧美精品| 欧美黄色淫秽网站| 久久 成人 亚洲| 中文亚洲av片在线观看爽| 国产成人系列免费观看| 五月开心婷婷网| 国产黄色免费在线视频| 久久久国产一区二区| 日日爽夜夜爽网站| 动漫黄色视频在线观看| 在线观看一区二区三区| 亚洲性夜色夜夜综合| 如日韩欧美国产精品一区二区三区| 欧洲精品卡2卡3卡4卡5卡区| 美国免费a级毛片| 国产精品综合久久久久久久免费 | 一二三四在线观看免费中文在| 欧美日本亚洲视频在线播放| 黄色视频不卡| www.www免费av| 国产亚洲欧美精品永久| 精品卡一卡二卡四卡免费| 国产99白浆流出| 18禁国产床啪视频网站| 人成视频在线观看免费观看| 操美女的视频在线观看| 丁香六月欧美| 成人特级黄色片久久久久久久| 日韩欧美一区视频在线观看| 久久中文字幕人妻熟女| 在线观看免费午夜福利视频| 国产精品久久久人人做人人爽| 免费高清在线观看日韩| 视频区图区小说| 50天的宝宝边吃奶边哭怎么回事| 久热爱精品视频在线9| 国产在线观看jvid| 90打野战视频偷拍视频| 一级作爱视频免费观看| av国产精品久久久久影院| 亚洲中文字幕日韩| 日韩精品免费视频一区二区三区| av在线天堂中文字幕 | 淫妇啪啪啪对白视频| 久热爱精品视频在线9| 免费少妇av软件| 午夜亚洲福利在线播放| 亚洲一卡2卡3卡4卡5卡精品中文| 国产人伦9x9x在线观看| 一级a爱视频在线免费观看| 亚洲成人精品中文字幕电影 | 两性夫妻黄色片| 精品第一国产精品| 丰满饥渴人妻一区二区三| 国产不卡一卡二| 日本免费a在线| 在线av久久热| 午夜日韩欧美国产| 日韩免费av在线播放| 国产黄a三级三级三级人| 波多野结衣一区麻豆| www.熟女人妻精品国产| 欧美日韩精品网址| 国产主播在线观看一区二区| 交换朋友夫妻互换小说| 欧美一级毛片孕妇| 欧美性长视频在线观看| 国产野战对白在线观看| 丝袜美腿诱惑在线| 两人在一起打扑克的视频| 国产成人免费无遮挡视频| 三上悠亚av全集在线观看| 免费高清视频大片| 黄片大片在线免费观看| 多毛熟女@视频| 国产乱人伦免费视频| 亚洲国产看品久久| 黄频高清免费视频| 大码成人一级视频| 免费在线观看视频国产中文字幕亚洲| 涩涩av久久男人的天堂| 成在线人永久免费视频| 亚洲免费av在线视频| 无限看片的www在线观看| 亚洲第一av免费看| 国产成人精品无人区| 精品久久久久久成人av| 国产精品1区2区在线观看.| 两个人免费观看高清视频| 夜夜爽天天搞| 正在播放国产对白刺激| 成人特级黄色片久久久久久久| 国产一卡二卡三卡精品| 十八禁人妻一区二区| 亚洲精品国产色婷婷电影| 欧美精品亚洲一区二区| 国产精品国产av在线观看| 一个人观看的视频www高清免费观看 | 亚洲一区二区三区色噜噜 | 久久久久国内视频| 成人18禁在线播放| 麻豆成人av在线观看| 在线永久观看黄色视频| 丝袜人妻中文字幕| 国产成人av激情在线播放| 国产99白浆流出| 天堂√8在线中文| 99在线视频只有这里精品首页| 亚洲 欧美一区二区三区| 999久久久国产精品视频| 欧美日韩中文字幕国产精品一区二区三区 | 老司机在亚洲福利影院| 在线观看一区二区三区| 亚洲国产欧美一区二区综合| 午夜成年电影在线免费观看| 变态另类成人亚洲欧美熟女 | 80岁老熟妇乱子伦牲交| 免费看十八禁软件| 9色porny在线观看| 悠悠久久av| 亚洲狠狠婷婷综合久久图片| 中出人妻视频一区二区| 老司机在亚洲福利影院| 久久天堂一区二区三区四区| 久久精品亚洲av国产电影网| xxxhd国产人妻xxx| 久久人人97超碰香蕉20202| 亚洲av电影在线进入| 少妇的丰满在线观看| 又黄又爽又免费观看的视频| 欧美在线黄色| 久久精品91无色码中文字幕| 久久热在线av| 亚洲专区中文字幕在线| 女性生殖器流出的白浆| 两人在一起打扑克的视频| 国产精品乱码一区二三区的特点 | 天天躁狠狠躁夜夜躁狠狠躁| cao死你这个sao货| 国产aⅴ精品一区二区三区波| 动漫黄色视频在线观看| 欧美激情极品国产一区二区三区| 精品久久久久久电影网| 超碰成人久久| av视频免费观看在线观看| 中文欧美无线码| 麻豆一二三区av精品| 久久国产乱子伦精品免费另类| 交换朋友夫妻互换小说| 最新美女视频免费是黄的| 国产精品99久久99久久久不卡| 宅男免费午夜| 99香蕉大伊视频| 好看av亚洲va欧美ⅴa在| 亚洲国产看品久久| a在线观看视频网站| 99精品在免费线老司机午夜| 久久欧美精品欧美久久欧美| 性色av乱码一区二区三区2| 夫妻午夜视频| 超碰97精品在线观看| 久久久久久久午夜电影 | 纯流量卡能插随身wifi吗| 叶爱在线成人免费视频播放| 久久伊人香网站| 亚洲精品中文字幕在线视频| 精品久久久久久久毛片微露脸| 少妇的丰满在线观看| 在线观看66精品国产| 在线天堂中文资源库| av片东京热男人的天堂| 国产精品1区2区在线观看.| 亚洲精品美女久久久久99蜜臀| 18禁黄网站禁片午夜丰满| 757午夜福利合集在线观看| 国产精品免费一区二区三区在线| 18禁黄网站禁片午夜丰满| 日韩国内少妇激情av| 他把我摸到了高潮在线观看| 欧美日韩一级在线毛片| 无限看片的www在线观看| 欧美午夜高清在线| 国产成人一区二区三区免费视频网站| 欧美人与性动交α欧美软件| 精品一品国产午夜福利视频| 如日韩欧美国产精品一区二区三区| 精品国产国语对白av| 男女下面进入的视频免费午夜 | 日本黄色视频三级网站网址| 香蕉丝袜av| 黄色丝袜av网址大全| 欧美一区二区精品小视频在线| 国内毛片毛片毛片毛片毛片| 夫妻午夜视频| 精品电影一区二区在线| 亚洲国产精品一区二区三区在线| 精品一区二区三区av网在线观看| 神马国产精品三级电影在线观看 | 狠狠狠狠99中文字幕| 成人18禁高潮啪啪吃奶动态图| 日韩精品青青久久久久久| 免费搜索国产男女视频| 久久香蕉国产精品| 久久久久国产一级毛片高清牌| 在线观看www视频免费| 91av网站免费观看| 91大片在线观看| 18禁观看日本| 日韩人妻精品一区2区三区| 可以在线观看毛片的网站| 啦啦啦免费观看视频1| 亚洲色图 男人天堂 中文字幕| 亚洲免费av在线视频| 欧美精品亚洲一区二区| 黄色视频不卡| 一级片'在线观看视频| 12—13女人毛片做爰片一| 高清黄色对白视频在线免费看| 免费观看精品视频网站| 欧美在线黄色| 99精品久久久久人妻精品| 曰老女人黄片| 亚洲av日韩精品久久久久久密| 国产亚洲欧美精品永久| 亚洲三区欧美一区| 精品人妻1区二区| 国产亚洲欧美在线一区二区| 欧美中文日本在线观看视频| 亚洲三区欧美一区| 91av网站免费观看| 50天的宝宝边吃奶边哭怎么回事| 一级a爱片免费观看的视频| 制服诱惑二区| 一本综合久久免费| 天堂动漫精品| 精品欧美一区二区三区在线| 久久婷婷成人综合色麻豆| 男女下面插进去视频免费观看| 免费在线观看日本一区| 男人舔女人下体高潮全视频| 中文字幕人妻熟女乱码| 曰老女人黄片| 国产成人精品在线电影| 搡老乐熟女国产| 亚洲欧美日韩无卡精品| 国产极品粉嫩免费观看在线| 午夜激情av网站| 日韩视频一区二区在线观看| 国产精品一区二区在线不卡| 亚洲美女黄片视频| 久久香蕉激情| 亚洲午夜精品一区,二区,三区| 亚洲av电影在线进入| 91大片在线观看| 国产真人三级小视频在线观看| 99精品久久久久人妻精品| 一二三四在线观看免费中文在| 免费高清在线观看日韩| 最新美女视频免费是黄的| 欧美av亚洲av综合av国产av| 真人一进一出gif抽搐免费| 欧美日韩精品网址| 国产精品自产拍在线观看55亚洲| 黄片播放在线免费| av在线播放免费不卡| 日韩一卡2卡3卡4卡2021年| 精品国产亚洲在线| 亚洲av熟女| 一区二区三区国产精品乱码| 一进一出好大好爽视频| 久久久国产成人精品二区 | 日韩av在线大香蕉| 好看av亚洲va欧美ⅴa在| 国产精品二区激情视频| 黑人猛操日本美女一级片| 国产亚洲精品第一综合不卡| 亚洲人成电影免费在线| 久久国产精品影院| 69av精品久久久久久| 欧美乱码精品一区二区三区| 久久香蕉国产精品| 99久久久亚洲精品蜜臀av| 夜夜爽天天搞| 老司机亚洲免费影院| 国产一区二区三区综合在线观看| 国产人伦9x9x在线观看| 在线观看一区二区三区激情| 午夜老司机福利片| 亚洲精品在线美女| 国产精品 国内视频| 国产免费男女视频| videosex国产| 人人妻人人添人人爽欧美一区卜| 久久精品国产亚洲av香蕉五月| e午夜精品久久久久久久| 欧美一区二区精品小视频在线| 视频在线观看一区二区三区| 亚洲一区中文字幕在线| 国产精品av久久久久免费| 国产99久久九九免费精品| 日韩欧美三级三区| 国产激情久久老熟女| 亚洲欧美日韩无卡精品| 国产成人一区二区三区免费视频网站| 天堂俺去俺来也www色官网| 色精品久久人妻99蜜桃| 国产精品乱码一区二三区的特点 | 中文字幕人妻丝袜一区二区| 丰满人妻熟妇乱又伦精品不卡| av在线播放免费不卡| 免费不卡黄色视频| 天天躁夜夜躁狠狠躁躁| 麻豆国产av国片精品| 999久久久精品免费观看国产| www.www免费av| 又紧又爽又黄一区二区| 欧美成人性av电影在线观看| 老司机深夜福利视频在线观看| 久久久精品欧美日韩精品| 国产99久久九九免费精品| 十八禁网站免费在线| 黄色a级毛片大全视频| 免费在线观看完整版高清| 女人高潮潮喷娇喘18禁视频| 人人妻人人澡人人看| 午夜视频精品福利| 国产精品成人在线| 亚洲 欧美 日韩 在线 免费| 1024视频免费在线观看| 水蜜桃什么品种好| 国产成人精品无人区| 高清欧美精品videossex| 精品国产乱码久久久久久男人| 中文字幕高清在线视频| 亚洲精品美女久久av网站| 黄片播放在线免费| 午夜福利,免费看| 中文字幕最新亚洲高清| 男人舔女人下体高潮全视频| 一级片'在线观看视频| 亚洲国产看品久久| 精品午夜福利视频在线观看一区| 久久久国产欧美日韩av| 老司机午夜十八禁免费视频| 99国产极品粉嫩在线观看| 久久这里只有精品19| 欧美+亚洲+日韩+国产| 日韩有码中文字幕| 动漫黄色视频在线观看| 久久性视频一级片| 超碰97精品在线观看| 制服诱惑二区| 身体一侧抽搐| 亚洲一卡2卡3卡4卡5卡精品中文| 国产无遮挡羞羞视频在线观看| 超色免费av| 妹子高潮喷水视频| 欧美日韩视频精品一区| 不卡一级毛片| 久久午夜亚洲精品久久| 激情视频va一区二区三区| 黄片大片在线免费观看| 日韩中文字幕欧美一区二区| 国产色视频综合| 欧美亚洲日本最大视频资源| 久9热在线精品视频| 欧美一区二区精品小视频在线| 欧美日韩一级在线毛片| 成人国产一区最新在线观看| 国产成人av教育| 国产精品影院久久| 亚洲欧美激情综合另类| 日韩欧美国产一区二区入口| 老司机午夜十八禁免费视频| 亚洲美女黄片视频| 色综合欧美亚洲国产小说| 国产真人三级小视频在线观看| 夫妻午夜视频| 欧美老熟妇乱子伦牲交| 我的亚洲天堂| 亚洲中文av在线| netflix在线观看网站| 日韩欧美免费精品| www日本在线高清视频| 国产成人精品久久二区二区91| 手机成人av网站| 亚洲人成电影观看| 欧美成人午夜精品| 午夜精品久久久久久毛片777| 好男人电影高清在线观看| aaaaa片日本免费| 亚洲国产精品sss在线观看 | 丝袜美足系列|