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

    改進(jìn)鼠群優(yōu)化算法在地下巖體破裂事件定位中的應(yīng)用研究

    2023-11-21 05:55:18白力名簡興祥
    物探化探計算技術(shù) 2023年5期
    關(guān)鍵詞:檢波器定位精度種群

    白力名, 簡興祥, 張 偉

    ( 1.成都理工大學(xué) 地球物理學(xué)院,成都 610059;2.中國地質(zhì)調(diào)查局 成都地質(zhì)調(diào)查中心,成都 610081;3.電子科技大學(xué) 資源與環(huán)境學(xué)院,成都 611731)

    0 引言

    地下巖體在開挖擾動作用下,圍巖內(nèi)的彈塑性勢能集中程度達(dá)到某一臨界值后,會誘發(fā)巖體中微裂隙的產(chǎn)生,同時伴隨著彈性勢能的釋放,這種現(xiàn)象稱為巖體破裂[1]。巖體破裂通常伴生在地下巖體中,通過實時監(jiān)測巖石破裂來對巖體演化進(jìn)行解釋和評價。因此,巖體破裂作為一種實時、被動的監(jiān)測手段被應(yīng)用到巖石力學(xué)與工程監(jiān)測領(lǐng)域中[2]。為了使實時監(jiān)測技術(shù)更好地發(fā)揮邊坡巖體穩(wěn)定性評價的作用[3],提高地下巖體破裂事件的定位精度,選擇合適的定位算法是巖體破裂事件的主要影響因素。

    最常見的定位算法,包括最小二乘法、Geiger 法[4]和單純形法,但都存在不同程度的局限性:①最小二乘法不一定能解出最優(yōu);②Geiger 法對初始值要求高;③單純形法易陷于局部最優(yōu)。為了提高定位精度,國內(nèi)、外學(xué)者做了許多研究,劉培洵等[5]采用絕對偏差最小法進(jìn)行巖石微破裂定位,減小了定位誤差;張曉平等[6]通過布谷鳥算法進(jìn)行層狀巖體破裂事件定位,提高了全局最優(yōu)能力。隨著計算機(jī)科學(xué)的不斷發(fā)展,越來越多的智能算法得到了快速發(fā)展[7],極大地提高了定位精度。最早提出的智能算法包括粒子群算法[8]、模擬退火算法[9]、遺傳算法[10]等,它們有著不同的優(yōu)缺點:①模擬退火算法可獲得全局最優(yōu)解,但搜索速度較慢; ②遺傳算法具有全局尋優(yōu)特點,但設(shè)置參數(shù)多;③粒子群算法簡單高效、全局尋優(yōu)能力好,但易陷入局部最優(yōu)。針對單個算法的局限性,許多學(xué)者進(jìn)行智能算法的改進(jìn)或多種算法相結(jié)合,達(dá)姝瑾等[11]提出了改進(jìn)粒子群算法在巖石聲發(fā)射定位中的研究,避免了局部最優(yōu)問題;LAGOS等[12]將粒子群算法與模擬退火法進(jìn)行聯(lián)合定位,提高了搜索效率,這些方法極大地提高了定位精度。

    隨著智能化的迅猛發(fā)展,許多智能優(yōu)化算法被開發(fā)出來[13]。比如鼠群優(yōu)化算法(Rat swarm optimization,Rso)[14]、灰狼算法、布谷鳥算法,其中鼠群優(yōu)化算法是Gaurav Dhiman等[15]在2020年提出的一種原理簡單、易于編程實現(xiàn)以及參數(shù)少、全局尋優(yōu)能力強(qiáng),搜索效率高等特點的元啟發(fā)式智能優(yōu)化算法。一經(jīng)問世就引起了各領(lǐng)域?qū)W者的廣泛關(guān)注,目前,鼠群優(yōu)化算法主要應(yīng)用在醫(yī)學(xué)診斷[16]、復(fù)合材料檢測[17]、電力系統(tǒng)[18]等領(lǐng)域中,但尚未應(yīng)用于地下巖體破裂事件定位中。因此,筆者主要探討針對改進(jìn)鼠群優(yōu)化算法在地下巖體破裂事件定位中的應(yīng)用研究。

    1 改進(jìn)鼠群優(yōu)化算法定位原理

    1.1 鼠群優(yōu)化算法

    鼠群優(yōu)化算法主要來自于大自然中老鼠的追逐和攻擊獵物的行為。為了避免老鼠相互碰撞,在追逐時每只老鼠的初始位置不同,老鼠在攻擊獵物時,鼠群做出逐次逼近的運(yùn)動形態(tài)。該優(yōu)化算法主要通過老鼠追逐獵物過程模擬局部搜索,鼠群攻擊過程模擬全局搜索。鼠群追逐獵物如圖1所示。

    (1)

    (2)

    t=0,1,2,…,MaxIteration

    (3)

    式中:R為控制變量A的參數(shù);MaxIteration為最大迭代次數(shù);t為當(dāng)前迭代次數(shù);參數(shù)R和C是取值在[1,5]和[0,2]之間的隨機(jī)數(shù)。

    當(dāng)老鼠與獵物進(jìn)行斗爭時,通過在搜索空間中不斷更新迭代來確定最優(yōu)位置。

    (4)

    1.2 Levy飛行策略

    Levy飛行是法國數(shù)學(xué)家萊維在上世紀(jì)30年代提出了基于一種概率分布[19]。它是一種非高斯隨機(jī)過程,通過小步長搜索和偶爾的大步長搜索相結(jié)合的搜索模式,增強(qiáng)了局部搜索能力和全局尋優(yōu)的能力,如圖2所示。

    圖2 二維Levy飛行示意圖Fig.2 2D Levy flight diagram

    Levy飛行的更新公式為式(5)。

    Yi(t+1)=Yi(t)+m⊕Levy(γ)

    (5)

    式中:Yi(t)為Yi在第t代的位置;m為i老鼠個體位置的隨機(jī)數(shù),由式(6)決定;⊕表示點對點乘法;Levy(γ)表示Levy飛行搜索路徑。

    m=random(size(iposition))

    (6)

    Levy(γ)=0.01(u*|v|γ*(Yi(t)-Yibest))

    (7)

    式中:iposition表示i老鼠個體位置;γ的取值范圍一般是1<γ<3[20],在本文中,γ取1.5;Yibest為歷史最優(yōu)i老鼠的位置;u和v服從公式(8)和(9)所示的正態(tài)分布:

    (8)

    (9)

    (10)

    (11)

    1.3 時差定位方法

    時差定位(Time difference of arrival, TDOA)是一種雙曲線定位法,該方法是通過不同檢波器接收到破裂點的時間差信息來建立方程組便可得到定位點[21]。在地下巖體破裂時,根據(jù)波的傳播特性,采用不同檢波器監(jiān)測到的P波到時差建立方程組,進(jìn)行定位求解。檢波器與巖體破裂如圖3所示。

    圖3 檢波器與巖體破裂示意圖Fig.3 Schematic diagram of geophone and rock mass fracture

    假設(shè)模擬的第i個檢波器坐標(biāo)為τ(xi,yi,zi),模擬的第j個巖體破裂事件坐標(biāo)為∈(xj,yj,zj),其中i,j=1、2、3、z,并且假設(shè)檢波器與巖體破裂事件位于同一巖層中。第i個檢波器與巖體破裂事件的P波測量時差關(guān)系為:

    (12)

    其中:Δ為第i個檢波器接收到的P波到時差;v為已知巖層速度。

    不同檢波器之間測量的P波到時差為:

    Δtn,m=tn-tm

    (13)

    通過已知波形速度V0,根據(jù)差值法(TD),計算理論模型中檢波器與破裂事件位置的時間差為:

    V0

    (14)

    式中:(xn,yn,zn)和(xm,ym,zm)分別為第n個和第m個檢波器位置;(x0,y0,z0)為破裂事件位置。

    根據(jù)檢波器之間的測量時差與計算時差的差值作為構(gòu)建的目標(biāo)函數(shù)。即公式(15)作為本文改進(jìn)的鼠群優(yōu)化算法的適應(yīng)度函數(shù),當(dāng)適應(yīng)度函數(shù)的取值越小時,則代表種群對函數(shù)的適應(yīng)性越好[15]。

    (15)

    1.4 改進(jìn)鼠群優(yōu)化算法的實現(xiàn)

    根據(jù)地下巖體破裂事件與檢波器的位置關(guān)系,利用雙曲線時差定位方法構(gòu)建目標(biāo)函數(shù)作為本文算法的適應(yīng)度函數(shù),將適應(yīng)度函數(shù)代入基于改進(jìn)的Levy飛行鼠群優(yōu)化算法中進(jìn)行優(yōu)化求解,將其最終優(yōu)化解作為定位坐標(biāo)解。算法流程圖如圖4所示,實驗步驟如下:

    圖4 算法流程圖Fig.4 Algorithm flow chart

    1)設(shè)置參數(shù),其中包括鼠群規(guī)模N、最大迭代次數(shù)T、維度、上下邊界等。

    2)隨機(jī)初始化鼠群位置。

    3)將初始化的鼠群位置進(jìn)行Levy飛行。

    4)將產(chǎn)生的每一個隨機(jī)位置代入適應(yīng)度函數(shù)計算其適應(yīng)度值。

    5)對每一個位置按鼠群優(yōu)化算法尋優(yōu)并處理邊界,選出當(dāng)前迭代次數(shù)下,適應(yīng)度值最小的位置,并與上一代比較,更新最優(yōu)位置。

    6)檢驗當(dāng)前迭代是否是最大迭代次數(shù),如未滿足最大迭代次數(shù),返回步驟5)再次進(jìn)行運(yùn)算。否則,終止算法。

    7)輸出最佳結(jié)果。

    2 模型計算

    2.1 模型設(shè)計

    筆者為了驗證基于Levy飛行的鼠群優(yōu)化算法對地下巖體破裂事件定位的有效性和可行性,設(shè)計了兩種不同地形的三維地下巖體破裂事件監(jiān)測系統(tǒng)模型。地下巖體破裂事件監(jiān)測系統(tǒng)如圖5所示。

    圖5 地下巖體破裂事件監(jiān)測系統(tǒng)Fig.5 Underground rock mass fracture event monitoring system

    針對兩種不同模型分別設(shè)置了9個檢波器,15 個破裂點。檢波器坐標(biāo)如表1所示,破裂點坐標(biāo)如表2所示。

    表1 檢波器坐標(biāo)Tab.1 Geophone coordinates

    表2 破裂點坐標(biāo)Tab.2 Rupture point coordinates

    2.2 參數(shù)驗證及實驗結(jié)果

    在模型實驗中,檢波器個數(shù)、種群數(shù)量、迭代次數(shù)以及不同地形對定位精度起著關(guān)鍵作用。針對兩種不同地形分別采用本文算法對破裂點1號(110,90,60)坐標(biāo)進(jìn)行定位誤差參數(shù)分析。在空間中最少需要四個臺站才可以判斷破裂點位置,因此筆者設(shè)置檢波器個數(shù)從4個增至9個,在不同地形、不同檢波器個數(shù)下,進(jìn)行了定位誤差對比,如圖6所示。

    圖6 檢波器對定位誤差的影響Fig.6 Influence of geophone on positioning error

    當(dāng)檢波器從4個增至6個時,針對兩種不同地形,本文算法對定位誤差都逐漸減小;當(dāng)檢波器從7個增至9個時,其定位誤差逐漸增大;故筆者采用檢波器個數(shù)為6。由實驗數(shù)據(jù)分析,當(dāng)檢波器一定時,水平地形的定位誤差要略小于傾斜地形的定位誤差。

    種群的數(shù)目對定位精度有著一定的影響,種群數(shù)量較大會導(dǎo)致算法的尋優(yōu)時間較長。而較小的種群數(shù)量會導(dǎo)致算法陷入局部最優(yōu)。因此,設(shè)定種群的數(shù)目從10等距增加到110。隨著種群個數(shù)的增加,定位精度也相對提高,而當(dāng)種群個數(shù)超過60后,定位精度的變化范圍不再明顯,基本維持不變,故筆者采用種群個數(shù)為60。種群對定位誤差的影響,如圖7所示。

    圖7 種群對定位誤差的影響Fig.7 Influence of population on positioning error

    迭代次數(shù)對定位精度也有著一定的影響,較大的迭代次數(shù)會導(dǎo)致算法的尋優(yōu)時間較長;而較小的迭代次數(shù)會導(dǎo)致算法陷入局部最優(yōu)。迭代次數(shù)從10等值增加到80,當(dāng)?shù)螖?shù)達(dá)到50次后,再增加迭代次數(shù),定位精度變化很小,從而設(shè)定最大迭代次數(shù)為50。迭代次數(shù)對定位影響,如圖8所示。

    圖8 迭代次數(shù)對定位誤差的影響Fig.8 Influence of iteration number on positioning error

    通過對改進(jìn)的鼠群優(yōu)化算法的影響參數(shù)實驗分析,筆者采用六個檢波器,分別是編號為A、C、D、F、G和I。設(shè)置種群數(shù)為60,迭代次數(shù)為50的參數(shù),針對兩種不同的地形模型對十五個破裂點進(jìn)行反演定位求解,得到破裂點1號(110,90,60)的三維定位位置圖和十五個破裂點定位結(jié)果統(tǒng)計表,分別如圖9和表3所示。

    表3 十五個破裂點定位結(jié)果Tab.3 Location results of fifteen fracture points

    圖9 破裂點1號的三維定位位置結(jié)果Fig.9 3D positioning results of rupture point

    2.3 評價與對比

    通過引入均方根誤差公式和迭代距離(Generational Distance,GD)評價指標(biāo),進(jìn)一步驗證本文算法的魯棒性和收斂性。

    均方根誤差公式(Root Mean Square Error RMSE)是評價定位精度的指標(biāo)之一,當(dāng)誤差越小,表示精度越高;同時,當(dāng)誤差趨于穩(wěn)定,表示算法的魯棒性較好。均方根誤差公式為式(16)。

    RMSE=

    (16)

    其中:(xj,yj,zj)為第j個真實解的坐標(biāo);(xbesti,ybesti,zbesti)為最優(yōu)解坐標(biāo);i為計算次數(shù)。

    分別采用狼群優(yōu)化算法(GWO)、鼠群優(yōu)化算法(RSO)以及本文改進(jìn)鼠群優(yōu)化算法分別對兩種不同的地形模型的十五個破裂點進(jìn)行均方根誤差計算,對每一個破裂點進(jìn)行10次計算,取平均值作為定位誤差結(jié)果,不同算法的定位誤差結(jié)果,如圖10所示。

    圖10 不同算法的定位誤差結(jié)果Fig.10 Positioning error results of different algorithms

    通過對圖10進(jìn)行分析,可以看出本文算法的定位誤差小于另外兩種算法的定位誤差,并且定位誤差結(jié)果都較趨于穩(wěn)定,進(jìn)一步說明本文算法的定位精度要高于其他兩種算法的定位精度,同時,表明該算法的魯棒性較好。

    迭代距離是衡量真實解與近似解集之間的間隔距離的平均值,當(dāng)GD值越小,表示該算法具有越好的收斂性,越能逼近真實解。迭代距離公式為式(17)。

    (17)

    其中:(xT,yT,zT)為真實解的坐標(biāo);(xRm,yRm,zRm)為第m個近似解的坐標(biāo);m為近似解的個數(shù)。

    同樣,采用上述三種算法,對兩種模型的十五個破裂點分別進(jìn)行迭代距離計算。不同算法的迭代距離統(tǒng)計,如圖11所示。

    圖11 不同算法的迭代距離統(tǒng)計Fig.11 Generational distance statistics of different algorithms

    通過對圖11進(jìn)行分析,可以看出本文算法的GD值小于另外兩種算法的GD值,說明本文算法的收斂性要好于另外兩種算法的收斂性,進(jìn)一步驗證了該算法性能的優(yōu)越性。

    3 結(jié)論

    通過對改進(jìn)鼠群優(yōu)化算法在地下巖體破裂事件的定位研究,得到以下幾點認(rèn)識:

    1)筆者提出了基于改進(jìn)鼠群優(yōu)化算法,增強(qiáng)了算法在進(jìn)行地下巖體破裂事件定位時的局部搜索能力和全局尋優(yōu)能力。通過模型設(shè)計,依據(jù)實驗所獲得的數(shù)據(jù)及三維定位結(jié)果分析,驗證了本文算法的有效性和可行性。

    2)采用均方根誤差公式和迭代距離評價指標(biāo),驗證了該算法性能的優(yōu)越性;在同等條件下,通過與另外兩種算法對比,本文定位算法具有較高的精度,較好的魯棒性和收斂性。

    3)由于本文方法僅是理論實驗,未在實際中應(yīng)用,還需結(jié)合實際地下層狀巖體破裂事件,考慮如何準(zhǔn)確構(gòu)建多層介質(zhì)模型和精確拾取到時等問題。因此,下一步主要研究工作是解決多層介質(zhì)模型和拾取到時等相關(guān)問題,并將理論研究成果成功應(yīng)用到實際工區(qū)中,從而實現(xiàn)快速、準(zhǔn)確地定位,以達(dá)到地下巖體破裂實時監(jiān)測的目的。

    猜你喜歡
    檢波器定位精度種群
    邢氏水蕨成功繁衍并建立種群 等
    北斗定位精度可達(dá)兩三米
    軍事文摘(2023年4期)2023-04-05 13:57:35
    山西省發(fā)現(xiàn)刺五加種群分布
    檢波器容差對地震信號接收的影響研究
    GPS定位精度研究
    智富時代(2019年4期)2019-06-01 07:35:00
    組合導(dǎo)航的AGV定位精度的改善
    一種井下檢波器測試儀的設(shè)計
    電子測試(2018年6期)2018-05-09 07:31:34
    基于高低頻聯(lián)測的常規(guī)檢波器數(shù)據(jù)低頻振幅和相位同時恢復(fù)方法
    崗更湖鯉魚的種群特征
    一種增大動圈式檢波器靜磁場的方法及一種動圈式檢波器
    97超碰精品成人国产| 亚洲欧美日韩东京热| 国产精品.久久久| 日韩一本色道免费dvd| 亚洲乱码一区二区免费版| 精品午夜福利在线看| 国产黄片视频在线免费观看| 超碰97精品在线观看| 午夜视频国产福利| 精品人妻偷拍中文字幕| 丰满乱子伦码专区| 一级黄片播放器| 大又大粗又爽又黄少妇毛片口| 高清午夜精品一区二区三区| 国产精品一及| 毛片女人毛片| av在线老鸭窝| 国产免费一级a男人的天堂| 久久精品夜夜夜夜夜久久蜜豆| 国产精品女同一区二区软件| 一级av片app| 午夜视频国产福利| 中文欧美无线码| 中文字幕熟女人妻在线| 少妇熟女欧美另类| 国产精品99久久久久久久久| av免费观看日本| 欧美性猛交╳xxx乱大交人| 精品人妻熟女av久视频| 看黄色毛片网站| 欧美精品一区二区大全| 一级毛片电影观看 | 亚洲av二区三区四区| 少妇人妻一区二区三区视频| av视频在线观看入口| 又爽又黄无遮挡网站| 日韩在线高清观看一区二区三区| 日本色播在线视频| 免费搜索国产男女视频| 中文字幕人妻熟人妻熟丝袜美| 免费无遮挡裸体视频| 日本午夜av视频| 久久综合国产亚洲精品| 特大巨黑吊av在线直播| 国产精品1区2区在线观看.| 26uuu在线亚洲综合色| 精品免费久久久久久久清纯| 成人二区视频| 久久久久久久久久久免费av| 日韩大片免费观看网站 | 精品少妇黑人巨大在线播放 | 尤物成人国产欧美一区二区三区| 岛国毛片在线播放| 亚洲精品成人久久久久久| 成人漫画全彩无遮挡| 免费在线观看成人毛片| 久久鲁丝午夜福利片| av国产免费在线观看| 一夜夜www| 又爽又黄无遮挡网站| 国产精品99久久久久久久久| 色综合亚洲欧美另类图片| 能在线免费看毛片的网站| 别揉我奶头 嗯啊视频| 波野结衣二区三区在线| 我的女老师完整版在线观看| 国产精品嫩草影院av在线观看| 男女那种视频在线观看| 91久久精品电影网| 亚洲av电影在线观看一区二区三区 | 欧美又色又爽又黄视频| 中国国产av一级| 欧美xxxx性猛交bbbb| a级一级毛片免费在线观看| 中文天堂在线官网| 超碰av人人做人人爽久久| 麻豆av噜噜一区二区三区| 中文资源天堂在线| av在线观看视频网站免费| 国产高清国产精品国产三级 | 精品久久久久久久久亚洲| 免费观看性生交大片5| 晚上一个人看的免费电影| 熟妇人妻久久中文字幕3abv| 亚洲美女视频黄频| 好男人在线观看高清免费视频| 亚洲在线观看片| 麻豆久久精品国产亚洲av| 久久精品影院6| 国产真实乱freesex| 国产在线男女| 欧美日韩精品成人综合77777| 两个人视频免费观看高清| 精品99又大又爽又粗少妇毛片| 国产成人freesex在线| www.av在线官网国产| 美女大奶头视频| 一区二区三区四区激情视频| 国产大屁股一区二区在线视频| 最近中文字幕2019免费版| 亚洲av免费在线观看| 一本一本综合久久| 欧美日韩国产亚洲二区| 99久久精品热视频| 国语对白做爰xxxⅹ性视频网站| 国产精品久久久久久精品电影小说 | 亚洲最大成人中文| 七月丁香在线播放| 亚洲欧美一区二区三区国产| 高清视频免费观看一区二区 | 国产精品久久久久久av不卡| 天堂网av新在线| 欧美一区二区精品小视频在线| 国产精品三级大全| 欧美成人免费av一区二区三区| 久久久久久久久久久免费av| 国产国拍精品亚洲av在线观看| 国产精品美女特级片免费视频播放器| 亚洲精华国产精华液的使用体验| 久久婷婷人人爽人人干人人爱| 国语对白做爰xxxⅹ性视频网站| 久久久精品大字幕| 亚洲精品日韩在线中文字幕| 欧美zozozo另类| 亚洲精品一区蜜桃| 亚洲欧美精品自产自拍| 久久久久久国产a免费观看| 韩国av在线不卡| 欧美+日韩+精品| 免费观看a级毛片全部| 亚洲欧美日韩卡通动漫| 日本一本二区三区精品| 国产精品无大码| 欧美zozozo另类| 国产大屁股一区二区在线视频| 日韩欧美 国产精品| АⅤ资源中文在线天堂| 久久精品久久久久久久性| 一个人看视频在线观看www免费| 九九久久精品国产亚洲av麻豆| 99久国产av精品国产电影| 免费一级毛片在线播放高清视频| 免费观看a级毛片全部| 欧美成人a在线观看| 国产色爽女视频免费观看| 男女啪啪激烈高潮av片| 少妇猛男粗大的猛烈进出视频 | 成人二区视频| 身体一侧抽搐| 久久精品久久精品一区二区三区| 99热这里只有是精品50| 精华霜和精华液先用哪个| 午夜免费激情av| 成人综合一区亚洲| 3wmmmm亚洲av在线观看| 日韩 亚洲 欧美在线| 91狼人影院| 热99在线观看视频| 51国产日韩欧美| 熟妇人妻久久中文字幕3abv| 五月玫瑰六月丁香| 欧美极品一区二区三区四区| 亚洲精品aⅴ在线观看| 超碰av人人做人人爽久久| 99久久成人亚洲精品观看| 我的老师免费观看完整版| 我要看日韩黄色一级片| 久久久久久久亚洲中文字幕| 免费看a级黄色片| 91精品伊人久久大香线蕉| 乱码一卡2卡4卡精品| 久久久久久九九精品二区国产| 久久99热这里只频精品6学生 | 欧美日韩精品成人综合77777| 最近手机中文字幕大全| 久久久久免费精品人妻一区二区| 亚洲怡红院男人天堂| 国产精品一区www在线观看| 国产精品乱码一区二三区的特点| 亚洲欧美成人综合另类久久久 | 国产三级在线视频| 大香蕉97超碰在线| 九草在线视频观看| 久久久精品94久久精品| 日韩欧美国产在线观看| 成人三级黄色视频| 亚洲欧美清纯卡通| 国产精品乱码一区二三区的特点| 男女边吃奶边做爰视频| 免费黄网站久久成人精品| 欧美精品一区二区大全| 丰满少妇做爰视频| 国产精品99久久久久久久久| 亚洲内射少妇av| 久久久国产成人精品二区| 精品不卡国产一区二区三区| av在线天堂中文字幕| 精品人妻一区二区三区麻豆| 91午夜精品亚洲一区二区三区| 亚洲av熟女| 亚洲在线观看片| 亚洲精品日韩av片在线观看| 人妻少妇偷人精品九色| 国产熟女欧美一区二区| 观看美女的网站| 老女人水多毛片| 国产亚洲5aaaaa淫片| 六月丁香七月| 国产成人一区二区在线| 99在线人妻在线中文字幕| 久久精品国产99精品国产亚洲性色| 九九在线视频观看精品| 中文字幕熟女人妻在线| 99热精品在线国产| 精品一区二区免费观看| 少妇的逼水好多| 亚洲欧美成人综合另类久久久 | 亚洲自偷自拍三级| 欧美极品一区二区三区四区| 日本五十路高清| 午夜免费激情av| 麻豆成人av视频| 99久久中文字幕三级久久日本| 国产亚洲一区二区精品| 国产免费福利视频在线观看| 久久午夜福利片| 97超碰精品成人国产| 免费电影在线观看免费观看| 亚洲欧美精品综合久久99| 国产精品久久久久久久电影| 欧美激情在线99| 秋霞伦理黄片| 色吧在线观看| 欧美3d第一页| 身体一侧抽搐| 男人的好看免费观看在线视频| 欧美激情国产日韩精品一区| 亚洲精品国产成人久久av| 99久久成人亚洲精品观看| 免费黄色在线免费观看| 欧美区成人在线视频| 丰满乱子伦码专区| 国产淫语在线视频| 色吧在线观看| 久久精品国产鲁丝片午夜精品| 亚洲精品亚洲一区二区| 午夜精品在线福利| 国产精品一区二区在线观看99 | 日本爱情动作片www.在线观看| 亚洲av成人精品一区久久| av免费观看日本| 久久精品夜色国产| 国产精品综合久久久久久久免费| 最近手机中文字幕大全| 三级男女做爰猛烈吃奶摸视频| 狂野欧美激情性xxxx在线观看| 国产乱来视频区| 你懂的网址亚洲精品在线观看 | 夜夜看夜夜爽夜夜摸| 久久综合国产亚洲精品| 中文字幕久久专区| 99久久中文字幕三级久久日本| 蜜桃亚洲精品一区二区三区| 久久99蜜桃精品久久| 久久久久网色| 高清午夜精品一区二区三区| av国产久精品久网站免费入址| 久久亚洲国产成人精品v| 亚洲熟妇中文字幕五十中出| 真实男女啪啪啪动态图| 亚洲精品成人久久久久久| 51国产日韩欧美| 亚洲电影在线观看av| 免费大片18禁| 99热精品在线国产| 亚洲av福利一区| av线在线观看网站| 黄片wwwwww| 乱系列少妇在线播放| 三级男女做爰猛烈吃奶摸视频| 天美传媒精品一区二区| 亚洲三级黄色毛片| 久久久久久久久久黄片| 尤物成人国产欧美一区二区三区| 国产精品av视频在线免费观看| 人妻系列 视频| 午夜福利在线观看免费完整高清在| 91精品伊人久久大香线蕉| 国产伦一二天堂av在线观看| 国产视频首页在线观看| 看黄色毛片网站| 亚洲国产欧美在线一区| 亚洲精品一区蜜桃| 最新中文字幕久久久久| 嘟嘟电影网在线观看| 日日啪夜夜撸| 久久人人爽人人片av| 国产成人精品久久久久久| 99久久无色码亚洲精品果冻| 伦理电影大哥的女人| 欧美性感艳星| 国产午夜精品论理片| АⅤ资源中文在线天堂| 亚洲三级黄色毛片| 长腿黑丝高跟| 久久久精品大字幕| 亚洲熟妇中文字幕五十中出| 真实男女啪啪啪动态图| 国产精品国产三级国产av玫瑰| 女的被弄到高潮叫床怎么办| 精品久久久久久久久亚洲| 久久精品久久久久久久性| 亚洲精品自拍成人| 国产精品美女特级片免费视频播放器| 国产爱豆传媒在线观看| 精品一区二区三区视频在线| 午夜福利网站1000一区二区三区| 国产精品熟女久久久久浪| 一级黄色大片毛片| 日韩高清综合在线| 最近中文字幕2019免费版| 国产成人a∨麻豆精品| 天天一区二区日本电影三级| 日日干狠狠操夜夜爽| 人人妻人人看人人澡| 又粗又硬又长又爽又黄的视频| 国产 一区 欧美 日韩| 国内揄拍国产精品人妻在线| 欧美极品一区二区三区四区| 久久精品夜夜夜夜夜久久蜜豆| 极品教师在线视频| 免费观看a级毛片全部| 亚洲一区高清亚洲精品| 亚洲精品日韩av片在线观看| 女的被弄到高潮叫床怎么办| 男人舔女人下体高潮全视频| 欧美成人午夜免费资源| 日韩精品有码人妻一区| 国产91av在线免费观看| 少妇丰满av| 蜜桃亚洲精品一区二区三区| 欧美bdsm另类| 久久精品国产亚洲av天美| 五月伊人婷婷丁香| 美女高潮的动态| 免费看a级黄色片| 久久人人爽人人片av| 中文字幕av成人在线电影| 好男人在线观看高清免费视频| 免费观看在线日韩| 毛片女人毛片| 边亲边吃奶的免费视频| 亚洲人与动物交配视频| 亚洲精品久久久久久婷婷小说 | 国产av码专区亚洲av| videossex国产| 国产精品1区2区在线观看.| 色播亚洲综合网| 日本熟妇午夜| 欧美激情在线99| 少妇高潮的动态图| 99久久无色码亚洲精品果冻| 色网站视频免费| 欧美成人a在线观看| 国语自产精品视频在线第100页| 国产黄片视频在线免费观看| 白带黄色成豆腐渣| 亚洲精品影视一区二区三区av| 国产精品国产高清国产av| 99国产精品一区二区蜜桃av| 国产精品.久久久| 免费搜索国产男女视频| 精品国产一区二区三区久久久樱花 | 国产久久久一区二区三区| 欧美bdsm另类| 国产乱人偷精品视频| 热99在线观看视频| 国产亚洲av片在线观看秒播厂 | 久久久久久大精品| 色5月婷婷丁香| 乱人视频在线观看| 特级一级黄色大片| av视频在线观看入口| 少妇高潮的动态图| 一区二区三区乱码不卡18| 一级黄色大片毛片| 在现免费观看毛片| 又粗又爽又猛毛片免费看| 亚洲成人久久爱视频| 男人舔奶头视频| 久久精品夜色国产| 深夜a级毛片| 亚洲av不卡在线观看| 美女大奶头视频| 女的被弄到高潮叫床怎么办| 色综合亚洲欧美另类图片| 乱码一卡2卡4卡精品| 免费不卡的大黄色大毛片视频在线观看 | 99久久精品一区二区三区| 床上黄色一级片| 久久这里有精品视频免费| 丰满乱子伦码专区| 亚洲在久久综合| 国产黄片视频在线免费观看| 免费av不卡在线播放| 午夜福利在线观看吧| 亚洲中文字幕日韩| 久久久久网色| 蜜臀久久99精品久久宅男| 国产视频内射| 五月伊人婷婷丁香| 亚洲天堂国产精品一区在线| 久久精品影院6| 波多野结衣巨乳人妻| 欧美性猛交黑人性爽| 国产又色又爽无遮挡免| 亚洲av一区综合| 日韩国内少妇激情av| 一区二区三区四区激情视频| 精品国产一区二区三区久久久樱花 | 国产亚洲av片在线观看秒播厂 | 欧美一区二区亚洲| 精品酒店卫生间| 成人无遮挡网站| 亚洲av二区三区四区| 男插女下体视频免费在线播放| 久久精品91蜜桃| 久久久色成人| 亚洲精品亚洲一区二区| 极品教师在线视频| av免费在线看不卡| 久久精品熟女亚洲av麻豆精品 | 亚洲五月天丁香| 欧美日韩一区二区视频在线观看视频在线 | 日韩av在线免费看完整版不卡| 非洲黑人性xxxx精品又粗又长| 亚洲成人精品中文字幕电影| 欧美极品一区二区三区四区| 亚洲精品乱久久久久久| 秋霞伦理黄片| 搡老妇女老女人老熟妇| 人人妻人人澡人人爽人人夜夜 | 午夜激情福利司机影院| av国产免费在线观看| 高清在线视频一区二区三区 | 高清在线视频一区二区三区 | 中文字幕人妻熟人妻熟丝袜美| 国产精品一区二区三区四区久久| 午夜视频国产福利| 一区二区三区四区激情视频| av.在线天堂| 亚洲人与动物交配视频| 日韩中字成人| 国产亚洲av片在线观看秒播厂 | 日本爱情动作片www.在线观看| h日本视频在线播放| 久久久久九九精品影院| 日本免费一区二区三区高清不卡| 嫩草影院新地址| 波多野结衣巨乳人妻| 欧美一区二区亚洲| 一级黄色大片毛片| 国产美女午夜福利| 3wmmmm亚洲av在线观看| 成人美女网站在线观看视频| 18禁在线播放成人免费| 三级经典国产精品| 日韩亚洲欧美综合| 日韩精品青青久久久久久| 晚上一个人看的免费电影| 亚洲av男天堂| 亚洲丝袜综合中文字幕| 少妇裸体淫交视频免费看高清| 免费av不卡在线播放| 2021少妇久久久久久久久久久| 22中文网久久字幕| 亚洲av不卡在线观看| 97超视频在线观看视频| 欧美日本视频| 天天躁日日操中文字幕| 麻豆久久精品国产亚洲av| 久久99热6这里只有精品| 亚洲av成人精品一二三区| 又黄又爽又刺激的免费视频.| 亚洲欧洲日产国产| 亚洲国产精品合色在线| 美女cb高潮喷水在线观看| 男人舔奶头视频| 中文字幕免费在线视频6| 精品午夜福利在线看| 久久久久免费精品人妻一区二区| 一边摸一边抽搐一进一小说| 亚洲精品亚洲一区二区| 亚洲国产精品久久男人天堂| 欧美精品国产亚洲| 国产成人freesex在线| 男人的好看免费观看在线视频| 男人舔奶头视频| 成年av动漫网址| 国产成年人精品一区二区| 亚洲国产日韩欧美精品在线观看| 免费搜索国产男女视频| ponron亚洲| 亚洲国产成人一精品久久久| 天堂√8在线中文| 搞女人的毛片| 男人的好看免费观看在线视频| 日日啪夜夜撸| 国产在视频线在精品| 淫秽高清视频在线观看| av天堂中文字幕网| 一级黄片播放器| 九九久久精品国产亚洲av麻豆| 国产伦精品一区二区三区四那| 亚洲av熟女| 国产免费一级a男人的天堂| 人妻系列 视频| 久久精品夜夜夜夜夜久久蜜豆| 免费人成在线观看视频色| 我要搜黄色片| 国产极品精品免费视频能看的| 美女高潮的动态| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产美女午夜福利| 亚洲欧美成人综合另类久久久 | 岛国在线免费视频观看| 国产极品精品免费视频能看的| 啦啦啦啦在线视频资源| 黑人高潮一二区| 久久久久久久久中文| 嫩草影院入口| 国产乱人视频| 国产在线一区二区三区精 | 午夜精品一区二区三区免费看| 在线观看美女被高潮喷水网站| 国产高清不卡午夜福利| 免费黄色在线免费观看| 亚洲国产精品专区欧美| 日韩成人伦理影院| 午夜福利视频1000在线观看| 精品人妻视频免费看| 天堂网av新在线| 免费观看性生交大片5| 亚洲成色77777| 99热6这里只有精品| 国产精品久久久久久精品电影| 久久精品熟女亚洲av麻豆精品 | 男人的好看免费观看在线视频| 亚洲精品,欧美精品| 真实男女啪啪啪动态图| 天美传媒精品一区二区| 1000部很黄的大片| 国产精品蜜桃在线观看| 国产av不卡久久| 欧美不卡视频在线免费观看| 国产黄色视频一区二区在线观看 | 日韩一区二区视频免费看| 一个人观看的视频www高清免费观看| 日韩人妻高清精品专区| 亚洲美女搞黄在线观看| 国产精品爽爽va在线观看网站| 精品不卡国产一区二区三区| 日本黄色片子视频| 国产高清视频在线观看网站| 黄色一级大片看看| 成人性生交大片免费视频hd| 男女下面进入的视频免费午夜| 亚洲av中文字字幕乱码综合| 成人无遮挡网站| 村上凉子中文字幕在线| 中文字幕久久专区| 国产精品久久久久久av不卡| 日产精品乱码卡一卡2卡三| 丝袜喷水一区| 两个人视频免费观看高清| 国内少妇人妻偷人精品xxx网站| 国产91av在线免费观看| 淫秽高清视频在线观看| 国产精品久久久久久久电影| 久热久热在线精品观看| 国产国拍精品亚洲av在线观看| 久久久久久久久久久免费av| 免费在线观看成人毛片| 国产成年人精品一区二区| 九草在线视频观看| 免费黄色在线免费观看| 日韩欧美三级三区| 黄色一级大片看看| 韩国av在线不卡| 亚洲18禁久久av| 亚洲色图av天堂| 午夜激情福利司机影院| 99久久精品热视频| 久久久久久久久久成人| 国产人妻一区二区三区在| 国产精品久久久久久精品电影| 国产成人一区二区在线| 国产老妇伦熟女老妇高清| 久久婷婷人人爽人人干人人爱| av在线亚洲专区| 亚洲欧美中文字幕日韩二区| 久久精品久久久久久噜噜老黄 | 91精品伊人久久大香线蕉| av在线观看视频网站免费| 中文字幕久久专区| 国产精品久久视频播放| 亚洲精品,欧美精品| 精品国内亚洲2022精品成人| 精品久久久久久久末码| 亚洲精品,欧美精品| 国产精品电影一区二区三区| 欧美xxxx黑人xx丫x性爽| 村上凉子中文字幕在线| av专区在线播放| 男插女下体视频免费在线播放|