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

    水下定位聲線擾動分析與分段指數(shù)權(quán)函數(shù)設(shè)計

    2021-08-14 02:27:10王薪普薛樹強曲國慶劉以旭楊文龍
    測繪學(xué)報 2021年7期
    關(guān)鍵詞:權(quán)函數(shù)聲線聲速

    王薪普,薛樹強,曲國慶,劉以旭,楊文龍

    1.山東理工大學(xué),山東 淄博 255049;2.中國測繪科學(xué)研究院,北京 100830;3.山東科技大學(xué),山東 青島 266590

    水下高精度時空服務(wù)是海洋開發(fā)活動的重要支撐[1]。我國陸地大地基準(zhǔn)網(wǎng)建設(shè)相對完善,并且長期復(fù)測維持,具備厘米級精度[2-4]。隨著我國北斗系統(tǒng)的建成,構(gòu)建從深空到深海全域無縫的綜合PNT體系將成為未來大地測量與導(dǎo)航的重要技術(shù)方向[5-6]。文獻[7—8]較全面地描述海底大地基準(zhǔn)網(wǎng)建設(shè)的主要技術(shù)問題和海洋導(dǎo)航技術(shù)的研究現(xiàn)狀,水下高精度聲吶定位模型和算法是海底空間基準(zhǔn)構(gòu)建及其聲吶導(dǎo)航定位應(yīng)用的重要基礎(chǔ)。

    由于衛(wèi)星定位采用的無線電磁波在海水中衰減嚴(yán)重,而聲波在水中具有良好的傳播性,因此水下高精度導(dǎo)航定位多采用聲學(xué)定位系統(tǒng)[9-10]。然而,海洋環(huán)境復(fù)雜,水中介質(zhì)不均勻,聲速會隨溫度、深度、鹽度的變化而變化,聲波在海水中傳播時存在聲線彎曲和延遲誤差,是制約水下高精度定位的瓶頸問題[11-12]。一方面,水下定位通常需要通過聲速剖面觀測獲取定位所需的聲速信息。然而,而聲速剖面存在測量誤差以及時空代表性誤差等[13-15],進而影響水下高精度定位。另一方面,由于入射角越大聲線彎曲越明顯,從而影響測距精度進而影響定位精度[16-20]。

    為解決聲線彎曲誤差對高精度定位的影響,文獻[21—24]對聲線彎曲誤差問題提出了不同的解決方法。在衛(wèi)星定位方面,多位學(xué)者對隨機模型進行了大量研究,提出高度角隨機模型、信噪比隨機模型、最小范數(shù)二次無偏估計模型及其簡化形式、實時鐘差融合權(quán)函數(shù)設(shè)計等隨機模型[25-27]。文獻[28]給出彈性定位導(dǎo)航授時(positioning,navigation and timing,PNT)的基本定義,提出需要研究海洋環(huán)境自適應(yīng)彈性定位模型。構(gòu)建合理的隨機模型是水下高精度定位的關(guān)鍵技術(shù)之一,然而,目前多采用觀測量等權(quán)模型,已有高度角相關(guān)模型也多仿照GNSS提出。為提高水下定位精度,文獻[27]在分析隨機模型不完善基礎(chǔ)上,提出4種基于入射角的水下定位隨機模型,改善等權(quán)隨機模型的定位精度,其影響機理和模型的合理性還有待深入研究。

    本文考慮海洋聲速場變化對聲線的擾動,基于常梯度聲線跟蹤模型推導(dǎo)入射角、梯度以及深度擾動對聲線傳播水平位移的影響函數(shù)關(guān)系,通過試驗分析入射角出現(xiàn)擾動時對水平位移的影響,依此建立入射角相關(guān)的水下定位隨機模型,采用淺水湖水下實測數(shù)據(jù)進行試驗分析驗證,與已有的權(quán)函數(shù)模型進行對比分析。

    1 入射角、梯度及深度擾動對聲線傳播的影響

    受海水溫度、壓力、鹽度等因素影響,聲速在垂向方向存在明顯的分層變化現(xiàn)象,導(dǎo)致聲線傳播發(fā)生折射彎曲,影響定位精度。通常,在同一分層中,假設(shè)聲速梯度是一個常數(shù),其聲線折射如圖1所示。

    圖1 聲線折射Fig.1 Illustration on the sound ray refraction

    依據(jù)常梯度聲線跟蹤模型,聲線掠過某一水層時的水平位移y可表示為[13]

    y=R(cosβ-cosα)

    (1)

    式中,α為入射角;β為層內(nèi)的出射角度;R=-C0/(gsinα),為聲線傳播的曲率半徑;C0表示入射聲速。α與β滿足Snell方程,即

    sinβ=sinα(C0+gh)/C0

    (2)

    式中,g=ΔC/h為聲速梯度。顯然,若將β、R視為過程變量,可將水平位移表示為α、g、h的函數(shù),即

    y=f(α,g,h)=

    (3)

    在聲線跟蹤定位算法中,入射角是需要迭代計算的量,即可表示為聲速剖面觀測值的函數(shù),因此也是一個隨機量。由式(3)可知聲線傳播的水平位移由入射角α、梯度g和深度h共同決定。為研究上述變量對水平位移的影響,對式(3)中的α、g、h求偏導(dǎo)

    (4)

    (5)

    (6)

    經(jīng)過化簡可得式(3)全微分為

    dy=(Rsinα-(Rtanβsinβ+y)cotα)dα+

    (7)

    類似于聲速剖面簡化思想,為了節(jié)省計算效率,建議將聲速剖面簡化為3~4個分層,其中每一層的聲速梯度利用最小二乘方法進行線性擬合,可得如下常梯度聲速剖面模型

    Ci=Ci,0+gih

    (8)

    式中,Ci,0為在第i個分層內(nèi)的初始聲速;gi為層內(nèi)的梯度(采用層內(nèi)聲速剖面數(shù)據(jù)的最小二乘估計);h∈[0,li]為層內(nèi)深度變量;li為第i層的厚度。此時,水平位移總擾動則可由多層擾動累加得到,即

    dy=dy1+dy2+…

    (9)

    需要指出,上述分層簡化策略只是為了計算方便。

    為了直觀給出水平位移與入射角的響應(yīng)關(guān)系,采用60 m實測聲速剖面進行試驗。聲速剖面和擬合剖面如圖2所示,聲速剖面大致可以分為兩部分,以分段方式進行聲速剖面擬合,計算入射角擾動對水平位移的影響。測試深度分別為30 m和60 m,入射角范圍為20°~80°,分別加入中誤差為0.005 rad~0.02 rad(0.286 5°~1.145 9°)的隨機誤差,模擬次數(shù)為1000次,統(tǒng)計1000次的水平位移擾動信息,計算水平位移擾動的樣本方差,結(jié)果如圖3所示。

    圖2 分段線性擬合聲速剖面Fig.2 Piecewise linear fitting sound velocity profile

    圖3 不同深度時入射角擾動對水平位移擾動的方差Fig.3 Variance of incidence angle disturbance to horizontal displacement disturbance at different depths

    給梯度加入-0.05~-0.15誤差,計算梯度擾動對水平位移的擾動,結(jié)果如圖4所示。

    圖4 不同深度下聲速梯度擾動對水平位移擾動Fig.4 The horizontal displacement is disturbed by gradient disturbance at different depths

    從圖4可以得出以下結(jié)論:

    (1)在入射角較小時,聲速場梯度變化對水平位移的擾動較小。

    (2)在相同水深和入射角下,水平位移的擾動隨著梯度變化增大而增大。

    (3)在相同深度和梯度下,水平位移擾動總是隨著入射角的增大而增大。

    在實際測量中,聲速場存在或多或少的測量誤差,并且聲速場是隨時間變化的,在聲線跟蹤計算時只能利用某一時刻的聲速場觀測信息代表聲吶觀測時刻的聲速場信息,因此總是存在聲速梯度誤差及其時變影響。據(jù)此,可基于上述響應(yīng)關(guān)系構(gòu)建水下定位隨機模型。

    2 顧及入射角的權(quán)函數(shù)設(shè)計

    假設(shè)Ti(i=1,2,…,n)時刻由GNSS和姿態(tài)測量給出船載聲吶換能器的位置Xi=(xi,yi,zi),海底應(yīng)答器待定坐標(biāo)為X=(x,y,z),聲吶信號傳播時間為ti,海底應(yīng)答器定位觀測模型為

    Li=cti=f(Xi,X)+Δi+εi(i=1,2,…,n)

    (10)

    (11)

    隨機模型的不完善通常歸結(jié)于定權(quán)不準(zhǔn)確[27],觀測權(quán)函數(shù)設(shè)計是構(gòu)建隨機模型的關(guān)鍵。對于相同的聲速場誤差和深度,大入射角對聲線的水平位移影響越大,進而影響水下定位精度。因此,構(gòu)建以下入射角相關(guān)觀測權(quán)函數(shù)

    (12)

    觀測方差大致可分為兩部分,即入射角相關(guān)部分和入射角不相關(guān)部分,后者主要來源于聲吶儀器測量誤差。當(dāng)入射角大時,則入射角相關(guān)誤差占主要成分,由圖3可知,此時儀器測量誤差可以忽略。從圖3中可以看出,無論30 m深度還是60 m深度時,當(dāng)入射角大于50°時,入射角擾動對水平位移的影響開始有上升趨勢,在深度60 m,入射角80°時,即使入射角加入0.005 rad的隨機誤差,水平位移的擾動方差達到了84.206 3 m2??梢姡?dāng)入射角較大時,入射角發(fā)生微小變化對水平的影響是指數(shù)增長的,即方差隨入射角的變化呈指數(shù)形式變化。為此,本文使用指數(shù)函數(shù)模型擬合觀測方差

    (13)

    式中,a、b為待估計量,可根據(jù)圖3相關(guān)數(shù)據(jù)進行擬合得到。以60 m水深0.005 rad誤差為例,擬合結(jié)果如圖5所示。

    圖5 方差擬合Fig.5 Variance fitting

    入射角擾動或聲速場擾動變化的具體量級未知或很難精確獲取。試驗表明,加入不同量級的擾動誤差,對聲線的擾動絕對量級影響很大,但對式(12)建立的權(quán)函數(shù)的擾動卻很小,其原因在于觀測權(quán)是觀測值間精度的相對量。如圖6所示,給出不同量級的入射角隨機擾動誤差,擬合入射角相關(guān)方差函數(shù)式(13),并采用式(12)構(gòu)建觀測權(quán)函數(shù),最終得到的權(quán)函數(shù)設(shè)計幾乎完全相同。說明本文所建立的權(quán)函數(shù)具有很好的適應(yīng)性,而入射角的擾動量級設(shè)置較為靈活。試驗表明,依據(jù)聲速梯度擾動對聲線傳播擾動影響進行權(quán)函數(shù)設(shè)計時,與依據(jù)入射角擾動得到的權(quán)幾乎完全相同。

    從上述權(quán)函數(shù)曲線可以發(fā)現(xiàn),分段指數(shù)權(quán)函數(shù)曲線和IGG3抗差方案具有類似的權(quán)函數(shù)圖形,因此其在抑制大觀測誤差方面的功效可較好體現(xiàn)[31-32]。

    圖7顯示在入射角0°~89°時的不同權(quán)函數(shù)的對比情況,直觀體現(xiàn)出分段指數(shù)權(quán)函數(shù)模型和已有的入射角指數(shù)函數(shù)模型、分段余弦函數(shù)模型。

    從圖7可以看出,指數(shù)模型在抑制大入射角時效果不夠明顯,分段余弦函數(shù)模型在分段節(jié)點處出現(xiàn)權(quán)值不連續(xù)現(xiàn)象。

    3 實測數(shù)據(jù)分析

    為驗證本文的權(quán)函數(shù)模型的有效性,采用淺水湖實測數(shù)據(jù)進行驗證分析。以圓標(biāo)校聲線跟蹤定位結(jié)果為未知點真值,測量船軌跡和定位點位置如圖8(a)所示,圓標(biāo)校引起對稱性可很好地消除聲線系統(tǒng)誤差影響,可獲取高精度平面定位[33]。

    采用如圖8(b)所示的導(dǎo)航軌跡數(shù)據(jù)比較分段指數(shù)權(quán)函數(shù)模型與等權(quán)模型、指數(shù)函數(shù)模型、分段余弦函數(shù)模型的定位效果,其C2、C6、C8各個控制點對應(yīng)的入射角如圖9所示。

    為分析本算例入射角相關(guān)誤差的大小及其與本文權(quán)函數(shù)設(shè)計的關(guān)系,采用圓測線標(biāo)校水下控制點結(jié)果作為參考值,計算觀測距離誤差,入射角和觀測誤差絕對值的對應(yīng)關(guān)系,如圖10所示(以C2為例)。結(jié)果表明,大入射角觀測的誤差相對更大,說明當(dāng)入射角大時,對測距的影響較大,從而影響定位精度,所以對大入射角觀測值進行降權(quán)處理是有必要的。

    基于前文試驗結(jié)果,使用入射角閾值為50°,以聲線跟蹤算法計算3個定位點C2、C6、C8的圓標(biāo)校結(jié)果,比較等權(quán)模型、入射角指數(shù)模型、分段余弦模型、分段指數(shù)權(quán)函數(shù)模型與圓標(biāo)校定位結(jié)果差異,結(jié)果見表1。

    從表1可以看出,分段指數(shù)權(quán)函數(shù)模型和分段余弦權(quán)函數(shù)模型的定位結(jié)果更接近圓標(biāo)校定位結(jié)果,兩者定位結(jié)果沒有顯著差異。以C2為例,相較于圓標(biāo)校結(jié)果,分段指數(shù)權(quán)函數(shù)模型比分段余弦權(quán)函數(shù)模型在X方向提高了0.019 m,Y方向降低0.013 m,Z方向提高了0.005 m;分段指數(shù)權(quán)函數(shù)模型的定位結(jié)果相較于等權(quán)模型的定位結(jié)果在X方向提升了4.134 m,Y方向提升了0.237 m,Z方向提升了3.098 m;指數(shù)模型和等權(quán)模型結(jié)果相近。C6分段指數(shù)權(quán)函數(shù)模型的定位結(jié)果相較于等權(quán)模型的定位結(jié)果在X方向上提升8.002 m;在Z方向上提升了5.965 m。而C8的定位精度改善不論是分段指數(shù)權(quán)函數(shù)模型還是分段余弦函數(shù)模型都主要在X方向。

    表1 不同方法定位結(jié)果對比Tab.1 Comparison of positioning results by different methods

    為了更加詳細(xì)地對比本文分段指數(shù)權(quán)函數(shù)模型、指數(shù)模型和分段余弦權(quán)函數(shù)模型,表2從協(xié)因數(shù)陣對角線之和、單位權(quán)中誤差和點位誤差3個方面進行對比。

    表2 3種模型對比Tab.2 Comparison of three models

    圖6 加入不同誤差時入射角與權(quán)值關(guān)系Fig.6 Relation ship of incident angle and weight

    圖7 不同權(quán)函數(shù)模型對比Fig.7 Comparison of different weight function models

    圖8 C2、C6、C8測量船航行軌跡及應(yīng)答器平面位置Fig.8 C2,C6,C8 tracking ship and transponder plane location

    圖9 C2、C6、C8控制點對應(yīng)的入射角Fig.9 C2,C6,C8 incident angle corresponding to the control point

    圖10 入射角與觀測誤差(取絕對值)的關(guān)系 Fig.10 Relation ship between incident angle and observation error (take the absolute value)

    結(jié)合圖7和表2可以看出權(quán)值P越小單位權(quán)中誤差越小,但是單位權(quán)中誤差應(yīng)表征實際的測量精度。本文提出的分段指數(shù)權(quán)函數(shù)模型更接近聲吶距離測量標(biāo)稱精度,即約為0.15 m。此外,協(xié)因數(shù)陣的對角元素之和代表觀測信息的利用率和觀測幾何精度。在單位權(quán)方差相同時,觀測幾何精度越小越好。指數(shù)模型的幾何精度最高,本文提出的模型幾何精度次之,但是指數(shù)函數(shù)模型權(quán)分配不合理,對大入射角的抑制作用不明顯。

    4 結(jié) 論

    類似于GNSS定位觀測,水下聲吶定位觀測中存在明顯的高度角相關(guān)誤差,且高度角越小誤差越大。研究表明,入射角相關(guān)誤差的產(chǎn)生主要來源于大入射角觀測對相同聲速場擾動、聲速測量誤差具有放大效應(yīng),入射角越大,這種放大效應(yīng)越明顯。在實測數(shù)據(jù)分析中,也發(fā)現(xiàn)了這種入射角相關(guān)誤差,且其變化趨勢與理論結(jié)果具有較好的一致性,均隨入射角變大呈現(xiàn)指數(shù)增長趨勢?;谏鲜鲈?,本文構(gòu)建的水下聲吶定位隨機模型,可有效抑制大入射角觀測條件下大誤差的影響,從而提高定位的精度和可靠性。

    聲速變化及其觀測的不確定以及入射角擾動大小很難精確給出。本文研究發(fā)現(xiàn),由于觀測權(quán)反映觀測間的相對精度,上述各類不確定性的絕對量級對觀測權(quán)值分配的影響很小,因此在實際應(yīng)用中可靈活設(shè)置。試驗表明,本文提出的分段指數(shù)權(quán)函數(shù)模型相對現(xiàn)有模型具有較為明顯的優(yōu)勢,其觀測不確定性意義明確,既可控制大入射角觀測影響,也可保留足夠的幾何觀測信息。

    需要指出的是,本文將入射角相關(guān)誤差歸結(jié)到隨機模型中,即設(shè)計了入射角相關(guān)的分段指數(shù)權(quán)函數(shù)模型。同時,該誤差也可歸入函數(shù)模型加以參數(shù)化估計,有望改善定位精度。實踐中,也可直接采用高度角變量(即觀測點與定位點之間連線與水平方向的夾角)構(gòu)造權(quán)函數(shù)的方法改善定位精度。本文的分段指數(shù)權(quán)函數(shù)模型目前只采用淺水湖水深60 m實測數(shù)據(jù)進行試驗驗證,在深海定位方面還有待繼續(xù)研究。

    猜你喜歡
    權(quán)函數(shù)聲線聲速
    基于改進權(quán)函數(shù)的探地雷達和無網(wǎng)格模擬檢測混凝土結(jié)構(gòu)空洞缺陷工程中的數(shù)學(xué)問題
    水聲中非直達聲下的聲速修正方法①
    一類廣義的十次Freud-型權(quán)函數(shù)
    基于聲線法的特殊體育館模型中聲場均勻性分析
    異徑電磁流量傳感器權(quán)函數(shù)分布規(guī)律研究*
    糾纏的曲線
    優(yōu)雅(2017年3期)2017-03-09 17:02:52
    聲速是如何測定的
    三維溫度梯度場中本征聲線軌跡的求取*
    跨聲速風(fēng)洞全模顫振試驗技術(shù)
    機翼跨聲速抖振研究進展
    亚洲欧美一区二区三区黑人 | 又粗又硬又长又爽又黄的视频| 色视频在线一区二区三区| 国产淫片久久久久久久久| 直男gayav资源| 久久亚洲国产成人精品v| 久久热精品热| 亚洲欧洲日产国产| 亚洲欧美日韩另类电影网站 | 国产精品国产av在线观看| 久久久久久久久久久免费av| 在线看a的网站| 免费观看性生交大片5| 一区二区三区乱码不卡18| 人妻少妇偷人精品九色| av在线亚洲专区| 亚洲丝袜综合中文字幕| 国产精品不卡视频一区二区| 久久女婷五月综合色啪小说 | 国产一区二区在线观看日韩| 亚洲美女视频黄频| www.色视频.com| 五月玫瑰六月丁香| 国产高清国产精品国产三级 | 在线免费十八禁| 熟妇人妻不卡中文字幕| 亚洲精品一区蜜桃| 国产69精品久久久久777片| 蜜桃久久精品国产亚洲av| 亚洲精品乱码久久久v下载方式| 久久久久久久亚洲中文字幕| 午夜免费男女啪啪视频观看| 日本熟妇午夜| 美女脱内裤让男人舔精品视频| 午夜福利高清视频| 亚洲欧美日韩无卡精品| 我要看日韩黄色一级片| 久久国产乱子免费精品| 久久久久久久久久成人| 97超碰精品成人国产| 国产在线男女| 国产一区二区三区综合在线观看 | 又爽又黄无遮挡网站| 精品久久久精品久久久| 我的女老师完整版在线观看| 亚洲一级一片aⅴ在线观看| 亚洲一级一片aⅴ在线观看| 日韩强制内射视频| 噜噜噜噜噜久久久久久91| 免费观看av网站的网址| 亚洲色图av天堂| 亚洲欧美精品自产自拍| 亚洲精品视频女| 日本wwww免费看| 久久精品久久久久久噜噜老黄| 一本色道久久久久久精品综合| 嫩草影院新地址| 18禁在线播放成人免费| 国产精品嫩草影院av在线观看| 国产人妻一区二区三区在| 亚洲熟女精品中文字幕| 亚洲人成网站高清观看| 欧美激情久久久久久爽电影| 91aial.com中文字幕在线观看| 国产一区二区亚洲精品在线观看| 天美传媒精品一区二区| 丝袜美腿在线中文| 久久精品夜色国产| 亚洲不卡免费看| 国产一区有黄有色的免费视频| 80岁老熟妇乱子伦牲交| av在线亚洲专区| 熟女电影av网| 国产精品无大码| 嫩草影院精品99| 水蜜桃什么品种好| 亚洲天堂国产精品一区在线| 欧美高清性xxxxhd video| av线在线观看网站| 亚洲国产av新网站| 亚洲精品久久久久久婷婷小说| av国产免费在线观看| 七月丁香在线播放| 18禁在线播放成人免费| 婷婷色综合www| 国产成人aa在线观看| 欧美极品一区二区三区四区| 亚洲三级黄色毛片| 国产精品一区二区三区四区免费观看| 1000部很黄的大片| 又大又黄又爽视频免费| 国产精品三级大全| 日韩欧美精品免费久久| 国产精品国产三级国产专区5o| 国产午夜精品一二区理论片| 欧美日韩精品成人综合77777| 国精品久久久久久国模美| 免费看日本二区| 成人欧美大片| 超碰97精品在线观看| 国产精品嫩草影院av在线观看| 大片免费播放器 马上看| 午夜激情福利司机影院| 九九爱精品视频在线观看| 夜夜爽夜夜爽视频| 看十八女毛片水多多多| 久久热精品热| 国产精品一及| 69人妻影院| 亚洲欧美中文字幕日韩二区| 亚洲图色成人| 亚洲欧美一区二区三区国产| 中文字幕亚洲精品专区| 久久久精品免费免费高清| 国产老妇女一区| 91aial.com中文字幕在线观看| 国产成人精品婷婷| 少妇猛男粗大的猛烈进出视频 | 成人亚洲精品一区在线观看 | 欧美日韩综合久久久久久| 免费高清在线观看视频在线观看| 国产成年人精品一区二区| 久久99热这里只频精品6学生| 下体分泌物呈黄色| 久久精品国产鲁丝片午夜精品| 欧美+日韩+精品| 亚洲av国产av综合av卡| 亚洲内射少妇av| 亚洲精品国产av成人精品| 日本黄色片子视频| 国产欧美日韩精品一区二区| 老司机影院毛片| 亚洲熟女精品中文字幕| 久久精品综合一区二区三区| 日本一二三区视频观看| 欧美极品一区二区三区四区| 欧美高清成人免费视频www| 亚洲第一区二区三区不卡| 日本午夜av视频| 成人特级av手机在线观看| 毛片女人毛片| 国产91av在线免费观看| 国产精品国产av在线观看| 日本免费在线观看一区| 日本与韩国留学比较| 国产精品人妻久久久影院| 久久97久久精品| 国产精品99久久99久久久不卡 | 精品99又大又爽又粗少妇毛片| 老师上课跳d突然被开到最大视频| 在线精品无人区一区二区三 | 国产黄频视频在线观看| av又黄又爽大尺度在线免费看| 嫩草影院入口| 99热这里只有精品一区| 日韩av免费高清视频| 日韩不卡一区二区三区视频在线| 国产综合精华液| 午夜爱爱视频在线播放| 国产免费又黄又爽又色| 亚洲国产高清在线一区二区三| 一级毛片 在线播放| 免费观看的影片在线观看| 免费少妇av软件| 三级国产精品欧美在线观看| 一级片'在线观看视频| 国产人妻一区二区三区在| 精品99又大又爽又粗少妇毛片| 亚洲国产精品999| 国产一区二区三区av在线| av在线播放精品| 亚洲精品国产av蜜桃| 精品少妇黑人巨大在线播放| 日韩欧美精品v在线| 国产成年人精品一区二区| 日韩欧美精品免费久久| 精品亚洲乱码少妇综合久久| 日韩成人伦理影院| 日本免费在线观看一区| av播播在线观看一区| 男女国产视频网站| 下体分泌物呈黄色| 青春草视频在线免费观看| 在线天堂最新版资源| 美女国产视频在线观看| 毛片一级片免费看久久久久| 欧美丝袜亚洲另类| 国产又色又爽无遮挡免| 男女啪啪激烈高潮av片| 精品久久久久久久久av| 久久久色成人| 观看免费一级毛片| 欧美xxⅹ黑人| 亚洲精品色激情综合| 精品国产乱码久久久久久小说| 国产午夜精品久久久久久一区二区三区| 欧美精品国产亚洲| av国产精品久久久久影院| 中文乱码字字幕精品一区二区三区| 国产 一区 欧美 日韩| 伊人久久精品亚洲午夜| 国产极品天堂在线| 日韩av免费高清视频| 在线 av 中文字幕| 天堂网av新在线| 亚洲第一区二区三区不卡| 一级毛片 在线播放| 一区二区三区乱码不卡18| 国产大屁股一区二区在线视频| 夜夜看夜夜爽夜夜摸| 亚州av有码| 18禁裸乳无遮挡免费网站照片| 尤物成人国产欧美一区二区三区| .国产精品久久| 大码成人一级视频| 狂野欧美白嫩少妇大欣赏| 秋霞伦理黄片| 婷婷色麻豆天堂久久| 少妇人妻一区二区三区视频| 亚洲欧美成人综合另类久久久| 久久人人爽av亚洲精品天堂 | 99久国产av精品国产电影| 欧美一级a爱片免费观看看| 亚洲欧美日韩无卡精品| 精品一区在线观看国产| 欧美xxxx性猛交bbbb| 中文字幕久久专区| 在线观看av片永久免费下载| 久久人人爽人人爽人人片va| 精品一区二区三区视频在线| 黄色欧美视频在线观看| 在线看a的网站| 777米奇影视久久| 亚洲国产精品成人久久小说| 可以在线观看毛片的网站| 肉色欧美久久久久久久蜜桃 | 老女人水多毛片| av天堂中文字幕网| 一级毛片 在线播放| 男人和女人高潮做爰伦理| 青春草视频在线免费观看| 亚洲一区二区三区欧美精品 | 熟女av电影| 亚洲精品日韩av片在线观看| 一本一本综合久久| 国产探花极品一区二区| 综合色av麻豆| 秋霞伦理黄片| 国产亚洲91精品色在线| 亚洲精品乱码久久久v下载方式| 久久久午夜欧美精品| 亚洲va在线va天堂va国产| 亚洲在久久综合| 国产成人精品婷婷| 国产色爽女视频免费观看| 国产一区二区三区av在线| 简卡轻食公司| 在线看a的网站| 男男h啪啪无遮挡| 国产又色又爽无遮挡免| 亚洲欧洲日产国产| 老师上课跳d突然被开到最大视频| 一区二区三区精品91| 老司机影院成人| 一个人看的www免费观看视频| 久久久精品欧美日韩精品| 在线观看一区二区三区激情| 日本色播在线视频| 丰满少妇做爰视频| 搡老乐熟女国产| 丰满少妇做爰视频| 国产一区有黄有色的免费视频| 日本黄大片高清| 精品一区在线观看国产| 亚洲最大成人手机在线| 18禁裸乳无遮挡免费网站照片| 丝瓜视频免费看黄片| 99精国产麻豆久久婷婷| 精品午夜福利在线看| 最近2019中文字幕mv第一页| 中文资源天堂在线| 日本色播在线视频| 性插视频无遮挡在线免费观看| 亚洲人成网站高清观看| 日韩av在线免费看完整版不卡| 综合色丁香网| av黄色大香蕉| 国产爱豆传媒在线观看| 国产精品女同一区二区软件| 欧美高清性xxxxhd video| 国产午夜福利久久久久久| 久久久精品欧美日韩精品| 九九爱精品视频在线观看| 国内精品美女久久久久久| 久久精品综合一区二区三区| 欧美丝袜亚洲另类| 亚洲精华国产精华液的使用体验| 亚洲熟女精品中文字幕| 午夜亚洲福利在线播放| 国产欧美日韩一区二区三区在线 | 99精国产麻豆久久婷婷| 欧美国产精品一级二级三级 | 自拍欧美九色日韩亚洲蝌蚪91 | 七月丁香在线播放| 下体分泌物呈黄色| 欧美人与善性xxx| 久久99热这里只有精品18| 国产淫片久久久久久久久| 亚洲最大成人中文| 69人妻影院| 欧美日本视频| 亚洲精品国产色婷婷电影| 国产免费又黄又爽又色| 国产男人的电影天堂91| 国产高清三级在线| 亚洲精品色激情综合| 国产乱人偷精品视频| 久久精品久久久久久噜噜老黄| 啦啦啦啦在线视频资源| 中文欧美无线码| 午夜激情福利司机影院| freevideosex欧美| 国产乱来视频区| 在线观看av片永久免费下载| 日本一二三区视频观看| 91精品一卡2卡3卡4卡| 伊人久久精品亚洲午夜| 精品国产露脸久久av麻豆| 国产欧美另类精品又又久久亚洲欧美| 国产黄色视频一区二区在线观看| 波多野结衣巨乳人妻| 日本午夜av视频| 免费观看的影片在线观看| 欧美日本视频| 国产精品av视频在线免费观看| 亚洲精品日本国产第一区| 国产探花极品一区二区| 免费不卡的大黄色大毛片视频在线观看| 99热全是精品| 黄色日韩在线| 国产精品三级大全| 狠狠精品人妻久久久久久综合| 亚洲av日韩在线播放| 精品少妇久久久久久888优播| 又黄又爽又刺激的免费视频.| 午夜福利网站1000一区二区三区| 女人十人毛片免费观看3o分钟| 男人添女人高潮全过程视频| 美女主播在线视频| 午夜精品国产一区二区电影 | 男男h啪啪无遮挡| 大香蕉久久网| 18禁在线无遮挡免费观看视频| 国产精品.久久久| 啦啦啦中文免费视频观看日本| 亚洲精品一二三| 欧美国产精品一级二级三级 | 国产精品麻豆人妻色哟哟久久| 国产精品一及| 在线观看一区二区三区激情| 国产国拍精品亚洲av在线观看| 午夜精品一区二区三区免费看| xxx大片免费视频| 一级爰片在线观看| 2018国产大陆天天弄谢| 波多野结衣巨乳人妻| 97在线视频观看| 国产综合懂色| 国产伦精品一区二区三区视频9| 亚洲不卡免费看| 一区二区av电影网| av黄色大香蕉| xxx大片免费视频| 欧美一级a爱片免费观看看| 精品午夜福利在线看| 欧美一级a爱片免费观看看| 国产精品人妻久久久影院| 成人美女网站在线观看视频| 九九久久精品国产亚洲av麻豆| 99久久人妻综合| 少妇人妻久久综合中文| 老女人水多毛片| 爱豆传媒免费全集在线观看| 日韩视频在线欧美| a级毛片免费高清观看在线播放| 少妇的逼水好多| 噜噜噜噜噜久久久久久91| 在线播放无遮挡| 97超视频在线观看视频| 蜜桃亚洲精品一区二区三区| 欧美极品一区二区三区四区| 在线精品无人区一区二区三 | 少妇人妻一区二区三区视频| 热99国产精品久久久久久7| 伦精品一区二区三区| 亚洲av成人精品一二三区| 水蜜桃什么品种好| 久久女婷五月综合色啪小说 | 又爽又黄无遮挡网站| 午夜福利网站1000一区二区三区| 免费观看无遮挡的男女| 小蜜桃在线观看免费完整版高清| 亚洲伊人久久精品综合| 亚洲精品第二区| 九九在线视频观看精品| 亚洲精品久久午夜乱码| a级一级毛片免费在线观看| 九九久久精品国产亚洲av麻豆| 欧美3d第一页| 91午夜精品亚洲一区二区三区| 欧美xxxx黑人xx丫x性爽| 国产精品99久久久久久久久| 成人特级av手机在线观看| 国产日韩欧美在线精品| 亚洲欧美成人综合另类久久久| 久久精品久久久久久久性| 少妇人妻一区二区三区视频| 在线免费观看不下载黄p国产| 免费观看在线日韩| av卡一久久| 人妻制服诱惑在线中文字幕| 三级国产精品欧美在线观看| 精品一区二区免费观看| 亚洲精品成人av观看孕妇| 国产精品久久久久久久电影| 国产黄片美女视频| 久久人人爽人人爽人人片va| 国产91av在线免费观看| 三级国产精品欧美在线观看| 亚洲aⅴ乱码一区二区在线播放| 欧美高清成人免费视频www| 大码成人一级视频| 日本三级黄在线观看| av国产精品久久久久影院| 欧美潮喷喷水| 97人妻精品一区二区三区麻豆| 亚洲精品成人av观看孕妇| 国产女主播在线喷水免费视频网站| 亚洲欧美一区二区三区国产| 成人黄色视频免费在线看| 99热这里只有精品一区| 亚洲精品色激情综合| 人妻 亚洲 视频| 青春草亚洲视频在线观看| 亚洲国产日韩一区二区| 能在线免费看毛片的网站| 人妻夜夜爽99麻豆av| 久久人人爽人人片av| 成人黄色视频免费在线看| 别揉我奶头 嗯啊视频| av又黄又爽大尺度在线免费看| 80岁老熟妇乱子伦牲交| av国产精品久久久久影院| 大陆偷拍与自拍| 国产片特级美女逼逼视频| freevideosex欧美| 成年av动漫网址| 97超视频在线观看视频| 18+在线观看网站| 亚洲,一卡二卡三卡| 精品国产一区二区三区久久久樱花 | 久久久亚洲精品成人影院| 美女高潮的动态| 在线观看一区二区三区| 99热这里只有是精品50| 嫩草影院精品99| 国产成人精品福利久久| 日本欧美国产在线视频| 亚洲精品,欧美精品| 亚洲美女搞黄在线观看| 九色成人免费人妻av| 国产精品99久久久久久久久| 性色avwww在线观看| 中文天堂在线官网| 婷婷色av中文字幕| 丰满人妻一区二区三区视频av| 人妻 亚洲 视频| 久久精品久久精品一区二区三区| 免费观看在线日韩| 色视频www国产| 99久久精品热视频| 美女被艹到高潮喷水动态| 色视频在线一区二区三区| 欧美日韩综合久久久久久| 精品人妻熟女av久视频| 久久精品国产自在天天线| 国产一区有黄有色的免费视频| 亚洲最大成人av| 五月开心婷婷网| 日本一二三区视频观看| 一区二区三区精品91| 亚洲av在线观看美女高潮| 有码 亚洲区| 天天一区二区日本电影三级| 2021少妇久久久久久久久久久| 男女边摸边吃奶| 特大巨黑吊av在线直播| 中文精品一卡2卡3卡4更新| 欧美老熟妇乱子伦牲交| 亚洲精品影视一区二区三区av| 国产免费一级a男人的天堂| 亚洲成人一二三区av| 99久久中文字幕三级久久日本| 看黄色毛片网站| 91精品伊人久久大香线蕉| 三级国产精品欧美在线观看| 美女视频免费永久观看网站| xxx大片免费视频| 欧美人与善性xxx| 亚洲欧美清纯卡通| 国产精品福利在线免费观看| 乱码一卡2卡4卡精品| 狂野欧美激情性xxxx在线观看| 国产黄片美女视频| 最近的中文字幕免费完整| 日韩伦理黄色片| 国产日韩欧美在线精品| 又粗又硬又长又爽又黄的视频| 久久精品熟女亚洲av麻豆精品| 免费av观看视频| 亚洲国产欧美人成| 精品一区二区三卡| 国产成人福利小说| 成人一区二区视频在线观看| 久久精品久久久久久噜噜老黄| 晚上一个人看的免费电影| 2022亚洲国产成人精品| 性色av一级| 国产真实伦视频高清在线观看| 水蜜桃什么品种好| 亚洲精品乱码久久久久久按摩| 国产黄色视频一区二区在线观看| 白带黄色成豆腐渣| 丰满乱子伦码专区| 久久亚洲国产成人精品v| 卡戴珊不雅视频在线播放| 美女脱内裤让男人舔精品视频| 久久国内精品自在自线图片| 成人毛片60女人毛片免费| 久久久久国产精品人妻一区二区| 蜜臀久久99精品久久宅男| 婷婷色av中文字幕| 久久精品国产自在天天线| 欧美成人精品欧美一级黄| www.色视频.com| 精品久久久噜噜| 国产亚洲最大av| 黄色配什么色好看| 亚洲精品国产av成人精品| 国产永久视频网站| 99精国产麻豆久久婷婷| 能在线免费看毛片的网站| 在线天堂最新版资源| 男人爽女人下面视频在线观看| 2021天堂中文幕一二区在线观| 九九爱精品视频在线观看| 黑人高潮一二区| 亚洲人成网站在线观看播放| 我要看日韩黄色一级片| 日韩成人伦理影院| 视频中文字幕在线观看| 亚洲国产日韩一区二区| 深爱激情五月婷婷| 女人久久www免费人成看片| 春色校园在线视频观看| 纵有疾风起免费观看全集完整版| 国产乱来视频区| 国产精品秋霞免费鲁丝片| 夜夜爽夜夜爽视频| 国产亚洲午夜精品一区二区久久 | 久久精品久久久久久噜噜老黄| 内地一区二区视频在线| 美女国产视频在线观看| 男女那种视频在线观看| 各种免费的搞黄视频| 高清在线视频一区二区三区| 2018国产大陆天天弄谢| 麻豆成人av视频| 国产高清有码在线观看视频| 欧美最新免费一区二区三区| 国产老妇伦熟女老妇高清| 日本与韩国留学比较| 亚洲精品中文字幕在线视频 | 亚洲精品国产av成人精品| 亚洲精品亚洲一区二区| 欧美丝袜亚洲另类| 美女主播在线视频| 内地一区二区视频在线| 国产精品国产三级国产av玫瑰| 日韩电影二区| 久久影院123| 国产精品.久久久| 国内精品宾馆在线| 亚洲av中文av极速乱| 亚洲精品456在线播放app| 美女国产视频在线观看| 可以在线观看毛片的网站| 18禁在线播放成人免费| h日本视频在线播放| 国产精品熟女久久久久浪| 亚洲精品视频女| 中文字幕亚洲精品专区| 国产精品秋霞免费鲁丝片| 中文字幕久久专区| 国产综合精华液| 亚洲成人av在线免费| 乱系列少妇在线播放| 久久国产乱子免费精品| 国产综合懂色| 男女国产视频网站| 国产成人精品福利久久| 男人和女人高潮做爰伦理| 自拍偷自拍亚洲精品老妇| 日本三级黄在线观看| 一级a做视频免费观看|