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

    一種適合開(kāi)采沉降實(shí)時(shí)觀測(cè)的單頻單歷元算法

    2015-05-06 01:26:52陳西斌張書畢
    金屬礦山 2015年2期
    關(guān)鍵詞:歷元正則實(shí)數(shù)

    陳西斌 張書畢 鮑 國(guó)

    (1.中國(guó)礦業(yè)大學(xué)環(huán)境與測(cè)繪學(xué)院,江蘇 徐州 221116;2.空軍勤務(wù)學(xué)院機(jī)場(chǎng)工程與保障系,江蘇 徐州 221116)

    一種適合開(kāi)采沉降實(shí)時(shí)觀測(cè)的單頻單歷元算法

    陳西斌1張書畢1鮑 國(guó)2

    (1.中國(guó)礦業(yè)大學(xué)環(huán)境與測(cè)繪學(xué)院,江蘇 徐州 221116;2.空軍勤務(wù)學(xué)院機(jī)場(chǎng)工程與保障系,江蘇 徐州 221116)

    為了彌補(bǔ)現(xiàn)有GPS單頻單歷元算法在礦區(qū)開(kāi)采沉降中的不足,更加精確地實(shí)時(shí)監(jiān)測(cè)煤炭開(kāi)采引起的地表移動(dòng),基于Tikhonov正則化方法及阻尼LAMBDA方法,提出了一種適合礦區(qū)沉降實(shí)時(shí)觀測(cè)的新方法。首先,根據(jù)單頻單歷元方程的特點(diǎn),利用Tikhonov正則化方法修正單頻單歷元法方程矩陣,由最小二乘法得到可靠的模糊度實(shí)數(shù)解;其次,依據(jù)阻尼LAMBDA方法搜索固定整周模糊度,其中阻尼因子可根據(jù)礦區(qū)首期變形監(jiān)測(cè)結(jié)果確定;最后,利用該算法處理某礦區(qū)地表沉降的實(shí)測(cè)GPS數(shù)據(jù),得到了測(cè)站北、東、天頂三方向的定位精度為毫米級(jí),且單歷元解算時(shí)間遠(yuǎn)小于采樣率。試驗(yàn)結(jié)果表明,該算法在礦區(qū)沉降中能滿足毫米級(jí)的觀測(cè),同時(shí)可滿足礦區(qū)沉降實(shí)時(shí)觀測(cè),對(duì)于實(shí)現(xiàn)對(duì)礦區(qū)沉降的實(shí)時(shí)觀測(cè)具有一定的參考價(jià)值。

    開(kāi)采沉降 GPS單頻單歷元 Tikhonov正則化 阻尼LAMBDA 阻尼因子

    為了保護(hù)工作面巷道、地面建筑物、構(gòu)筑物、水體及鐵路免受煤炭開(kāi)采掘進(jìn)的影響,必須充分研究地下開(kāi)采引起的巖層與地表移動(dòng)規(guī)律,以便最大限度減少地面財(cái)產(chǎn)的損失。由于GPS具有高精度、快速、全天候、多時(shí)段、操作方便等特點(diǎn),因而被廣泛應(yīng)用于開(kāi)采沉陷以及對(duì)開(kāi)采區(qū)內(nèi)各種建筑物、構(gòu)筑物的變形監(jiān)測(cè)中。在實(shí)時(shí)GPS觀測(cè)中,模糊度快速解算常用有OTF方法[1]、快速解算模糊度方法[2]等,該類方法均需要幾個(gè)歷元甚至幾十個(gè)歷元的初始化,導(dǎo)致解算時(shí)間過(guò)長(zhǎng),缺少實(shí)時(shí)性與有效性。GPS接收機(jī)易受到礦區(qū)復(fù)雜的觀測(cè)環(huán)境影響,導(dǎo)致周跳問(wèn)題嚴(yán)重,影響了定位解算結(jié)果精度。而利用單頻單歷元GPS相位觀測(cè)技術(shù),不僅克服了初始化時(shí)間過(guò)長(zhǎng)問(wèn)題以及降低了周跳的影響,而且單頻接收機(jī)硬件成本較雙頻接收機(jī)低很多。

    關(guān)于單頻單歷元相位解算,王德軍通過(guò)附加軌跡約束條件,并結(jié)合雙頻線性組合的寬巷觀測(cè)值,首先固定寬巷模糊度,然后固定載波L1 或L2 的模糊度,實(shí)現(xiàn)了單歷元固定模糊度的可能性[3];韓保民通過(guò)增加較強(qiáng)的坐標(biāo)三分量先驗(yàn)條件,改善單歷元法矩陣的病態(tài)性,直接得到模糊度實(shí)數(shù)解,在此基礎(chǔ)上根據(jù)沉降觀測(cè)中存在的特點(diǎn),成功地固定了模糊度,獲得了厘米量級(jí)的單歷元解算精度[4];唐衛(wèi)明提出了基線長(zhǎng)為約束的LAMBDA方法,在解算中使用兩步搜索,實(shí)現(xiàn)了單歷元定位[5]。該類算法的特點(diǎn)在于利用了先驗(yàn)信息,篩選了模糊度備選解或壓縮了模糊度搜索空間,減少了納偽時(shí)間的發(fā)生?;谏鲜鏊惴ǖ奶攸c(diǎn),結(jié)合礦區(qū)變形監(jiān)測(cè)的特點(diǎn),在分析了單歷元方程存在秩虧的基礎(chǔ)之上,采用Tikhonov正則化方法[6]有效地解決了秩虧問(wèn)題,得到了可靠的模糊度實(shí)數(shù)解;選擇阻尼LAMBDA方法[7]搜索固定模糊度。阻尼因子的選擇可根據(jù)礦區(qū)首次觀測(cè)求出的基線方差作為先驗(yàn)方差,根據(jù)單位權(quán)方差與先驗(yàn)方差求得阻尼因子,改進(jìn)了適合礦區(qū)實(shí)時(shí)變形監(jiān)測(cè)的單頻單歷元相位求解算法。

    1 單頻單歷元方程病態(tài)性的證明

    假設(shè)觀測(cè)歷元為t時(shí)刻,測(cè)站i、j共同觀測(cè)p1,p2,…,pk,…pn等n顆衛(wèi)星,以i測(cè)站為基準(zhǔn)站,j測(cè)站為流動(dòng)站,假設(shè)第pk顆衛(wèi)星的高度角最大,選其為參考衛(wèi)星。在測(cè)站間與衛(wèi)星間作雙差,得單歷元觀測(cè)方程

    (1)

    式(1)可簡(jiǎn)化成

    (2)

    式中,

    P與Q對(duì)應(yīng)雙差載波觀測(cè)值權(quán)矩陣與協(xié)因數(shù)矩陣。

    依據(jù)最小二乘原理解算式(2),可得

    (3)

    法方程為

    (4)

    令法矩陣系數(shù)陣為

    將法矩陣進(jìn)行奇異值分解,得

    F=USVT,

    (5)

    式中,U、V均為n+2階正交矩陣,由于F為對(duì)稱正定矩陣,所以U=V,S為F的奇異值矩陣,降序排列。將U、V、S按照以下形式分塊:

    (6)

    S1的3個(gè)奇異值很小,接近于0,間接致使法矩陣條件數(shù)很大,S1與S2之間有1個(gè)比較大的跳躍,導(dǎo)致了法矩陣病態(tài),這便為單頻單歷元雙差載波定位中法矩陣病態(tài)的共同點(diǎn)[8-10]。

    2 基于Tikhonov正則化方法獲得模糊度實(shí)數(shù)解

    Tikhonov正則化方法是由Tikhonov 1963 年提出來(lái)求解病態(tài)方程的常用方法,其實(shí)質(zhì)是通過(guò)附加部分參數(shù)或者全部參數(shù)的加權(quán)平方和極小條件,增加約束,補(bǔ)充先驗(yàn)信息來(lái)克服矩陣的病態(tài)性,使方程解穩(wěn)定且唯一。由式(3)構(gòu)造估計(jì)準(zhǔn)則函數(shù)[11]

    (7)

    (8)

    式中,R為正則化矩陣;R為不同對(duì)應(yīng)不同的正則化形式與結(jié)果;正則化參數(shù)α不同,解的穩(wěn)定性不同。圖1 為某算例某一歷元法矩陣正則化前后奇異值比較。

    圖1 法矩陣奇異值修正前后比較

    由圖1可知,正則化前法方程有3個(gè)奇異值遠(yuǎn)小于其余奇異值,證明了第1小節(jié)闡述的奇異值譜存在突變的情況;正則化后法方程最小的3個(gè)奇異值改變明顯,而法方程的特征譜變化不大。此外,正則化后,法矩陣條件數(shù)由4.591 0×105減少至1.228 5×103, 因此正則化使法方程的病態(tài)性得到了很大改善。

    3 改進(jìn)LAMBDA方法固定模糊度

    LAMBDA方法由Teunissen針對(duì)雙差模糊度之間相關(guān)性提出的一種降相關(guān)平差方法,是目前模糊度固定有效的方法之一。由第2節(jié)中得到的模糊度實(shí)數(shù)解,再依據(jù)LAMBDA方法搜索到模糊度固定解。LAMBDA實(shí)數(shù)模糊度與協(xié)方差矩陣[7]可表示成

    (10)

    (11)

    (12)

    (13)

    式中,Z為元素都為整數(shù)的矩陣。

    LAMBDA方法搜索到正確的模糊度的前提是必須提供精度較高的模糊度實(shí)數(shù)解及協(xié)因數(shù)陣,而單歷元方程秩虧,得到的實(shí)數(shù)解偏差較大,因此采用LAMBDA法很難搜索到正確模糊度,因此在式(10)中加入阻尼因子可得

    (14)

    (15)

    4 試驗(yàn)分析

    試驗(yàn)數(shù)據(jù)來(lái)自安徽省恒源煤電股份有限公司某礦區(qū)GPS變形監(jiān)測(cè)數(shù)據(jù),基準(zhǔn)站架設(shè)在基巖穩(wěn)定,地下及周邊無(wú)開(kāi)采的地帶,測(cè)站名為“stat”;監(jiān)測(cè)站架設(shè)在開(kāi)采影響范圍之內(nèi),預(yù)估計(jì)最大下沉速度即將來(lái)臨,測(cè)站名“dyna”?;€長(zhǎng)1 024.85 m,觀測(cè)時(shí)間2014年5月24日UTC 4時(shí)00分00秒至4時(shí)20分00秒。共觀測(cè)歷元1 201個(gè),采樣率為1 s,衛(wèi)星截止高度角為15o,每個(gè)歷元可觀測(cè)到8 顆衛(wèi)星。接收機(jī)板卡采用司南公司生產(chǎn)的K500 三系統(tǒng)單頻GNSS板卡(偽距精度為10 cm,載波相位精度為1.0 mm,采用司南導(dǎo)航公司研發(fā)的Compass Solution(簡(jiǎn)稱CS軟件)GNSS靜動(dòng)態(tài)后處理軟件,得到的每個(gè)歷元的位置信息,作為參考值。

    選取第1個(gè)歷元的觀測(cè)結(jié)果作比較,用CS軟件解算的L1雙差模糊度真值為

    為了比較文中所提出算法的有效性,采用2種計(jì)算方案:①方案1,普通最小二乘得到模糊度實(shí)數(shù)解,LAMBDA算法得到模糊度整數(shù)解;②方案2,Tikhonov正則化算法修正法方程得到模糊度的實(shí)數(shù)解,阻尼LAMBDA算法得到模糊度整數(shù)解。結(jié)果見(jiàn)表1。

    表1 固定模糊度比較

    由表2 可知,普通最小二乘得到的模糊度實(shí)數(shù)解與LAMBDA方法固定的模糊度與真值差較大;而經(jīng)過(guò)正則化修正之后的最小二乘模糊度實(shí)數(shù)解與阻尼LAMBDA方法固定的模糊度與真值差異較小。驗(yàn)證了Tikhonov正則化方法與阻尼LAMBDA方法結(jié)合可得到可靠的整周模糊度。

    采用Tikhonov正則化方法結(jié)合阻尼LAMBDA方法固定模糊度的單歷元在3個(gè)方向上解算參數(shù)見(jiàn)表2。

    表2 單歷元解算參數(shù)

    由表2可知,采用文中提出的模糊度固定方法的單歷元解算成功率為100%,且每個(gè)歷元解算平均時(shí)間遠(yuǎn)小于采樣率,可適用于1,2,10 Hz甚至更高采樣率變形監(jiān)測(cè)。北、東、天頂?shù)?個(gè)方向的定位精度在亞厘米量級(jí)甚至毫米級(jí),能夠滿足礦區(qū)的日常實(shí)時(shí)監(jiān)測(cè)要求。

    5 結(jié) 語(yǔ)

    在分析了單頻單歷元法方程特點(diǎn)的基礎(chǔ)之上,利用Tikhonov正則化方法修正了法方程矩陣,得到可靠的模糊度實(shí)數(shù)解;并結(jié)合阻尼LAMBDA搜索方法得到整周模糊度,解算結(jié)果精度維持在亞厘米級(jí)甚至毫米級(jí)。文中研究成果可以廣泛應(yīng)用在開(kāi)采沉陷中的地表移動(dòng)、變形沉降的日常實(shí)時(shí)監(jiān)測(cè)中。文中試驗(yàn)數(shù)據(jù)采用的是短基線,中長(zhǎng)基線、毫米級(jí)的單頻單歷元算法將是下一步的研究方向。

    [1] Abidin H Z.On-the-fly ambiguity resolution: formulation and results[J].Manuscripta Geodaetica,1993(18):380-405.

    [2] Xu P L,Shi C,Liu J N.Integer estimation methods for GPS ambiguity resolution: an application oriented review and improvement [J].Survey Review,2012(4):59-71.

    [3] 王德軍,熊永良,劉 寧,等.寬巷組合及軌跡約束下的列車在軌動(dòng)態(tài)單歷元定位算法[J].測(cè)繪學(xué)報(bào),2013,43(2):131-136. Wang Dejun,Xiong Yongliang,Liu Ning,et al.Train on track kinematic single epoch positioning algorithm based on wide-lane combination and track constraints[J].Acta Geodaetica et Cartographica Sinica,2013,43(2):131-136.

    [4] 韓保民,歐吉坤.一種附約束的單頻單歷元GPS雙差相位解算方法[J].測(cè)繪學(xué)報(bào),2002,31(4):300-303. Han Baomin,Ou Jikun.A GPS single epoch phase processing algorithm with constraints for single frequency receiver[J].Acta Geodaetica et Cartographica Sinica,2002,31(4):300-303.

    [5] 唐衛(wèi)明,孫紅星,劉經(jīng)南.附有基線長(zhǎng)度約束的單頻數(shù)據(jù)單歷元LAMBDA方法整周模糊度確定[J].武漢大學(xué)學(xué)報(bào):信息科學(xué)版,2005,30(5):444-446. Tang Weimin,Sun Hongxing,Liu Jingnan.Ambiguity resolution of single epoch single frequency data with baseline length constraint using LAMBDA algorithm [J].Geomatics and Information Science of Wuhan University,2005,30(5):444-446.

    [6] Tikhonov A N.Regularization of ill-posed problems[J].Dokl Akad Nauk SSSR,1963(1):49-52.

    [7] Teunissen P.J.G.The least-squares ambiguity de-correlation adjustment: a method for fast GPS integer ambiguity estimation[J].Journal of Geodesy,1995(1):65-82.

    [8] 歐吉坤.測(cè)量平差中不適定問(wèn)題解的統(tǒng)一表達(dá)與選權(quán)擬合法[J].測(cè)繪學(xué)報(bào),2004,33(4):283-288. Ou Jikun.Uniform expression of solutions of ill-posed problems in surveying adjustment and the fitting method by selection of the parameter weights[J].Acta Geodaetica et Cartographica Sinica,2004,33(4):283-288.

    [9] 王振杰,趙斌臣,黃汝杰.一種適合變形監(jiān)測(cè)的單頻單歷元算法[J].金屬礦山,2008(2):83-86. Wang Zhenjie,Zhao Binchen,Huang Rujie.A single epoch algorithm for deformation monitoring using single frequency GPS receivers[J].Metal Mine,2008(2):83-86.

    [10] 劉根友,郝曉光.GPS單歷元阻尼LAMBDA算法在水庫(kù)形變監(jiān)測(cè)中的應(yīng)用[J].大的測(cè)量與地球動(dòng)力學(xué),2004(2):64-69. Liu Genyou,Hao Xiaoguang.Damped LAMBDA algorithm of single epoch GPS positioning for reservoir-induced deformation monitoring[J].Journal of Geodesy and Geodynamics,2004(2):64-69.

    [11] Jazaeri S.Fast integer least-squares estimation for GNSS high-dimensional ambiguity resolution using lattice theory[J].Journal of Geodesy,2012,86(2):123-136.

    [12] Teunissenp J G.The LAMBDA method for the GNSS compass[J].Artificial Satellites,2006,41(3):89-103.

    [13] Chen W T,Qin H L.New method for single epoch,single frequency land vehicle attitude determination using low-end GPS receiver[J].GPS Solutions,2012(16):329-338.

    [14] Teunissenp J G.Integer least-squares theory for the GNSS compass[J].Journal of Geodesy,2010,84(1):433-447.

    [15] 劉小生,羅任秀.邊坡工程變形監(jiān)測(cè)系統(tǒng)研究[J].金屬礦山,2005(7):45-47. Liu Xiaosheng,Luo Renxiu.Research on monitoring system for slope engineering deformation[J].Metal Mine,2005(7):45-47.

    [16] Qin H L,Chen W T.Performance analysis and improved method for single epoch,single frequency attitude determination using GPS[J].Journal of PLA University of Science and Technology,2011,12(5):431-435.

    [17] 周玉新,項(xiàng)宏海,汪 斌.銅鑫礦業(yè)有限公司尾礦壩GPS位移檢測(cè)系統(tǒng)[J].金屬礦山,2013(6):100-103. Zhou Yuxin,Xiang Honghai,Wang Bin.GPS displacement monitoring system for tailings dam of Tongxin mining industry Co.Ltd.[J].Metal Mine,2013(6):100-103.

    (責(zé)任編輯 王小兵)

    A Single Frequency and Single Epoch Algorithm for Mining Settlement Observation in Real Time

    Chen Xibin1Zhang Shubi1Bao Guo2

    (1.SchoolofEnvironmentalScienceandSpatialInformatics,ChinaUniversityofMiningandTechnology,Xuzhou221116,China;2.AirForceServiceAcademy,AirportEngineeringandSecurityDepartment,Xuzhou221116,China)

    In order to make up for the shortage of the existing algorithm for mining subsidence by using single epoch and single frequency GPS observation,and realize real-time monitoring of the ground movement caused by coal mining more accurately,based on Tikhonov regularization method and damped LAMBDA method,a new algorithm for mining settlement in real time is proposed.Firstly,according to the characteristics of the single-frequency and single-epoch equation,Tikhonov regularization method is adopted to correct the matrix of the equation to obtain the reliable real-number solutions of ambiguity by the least squares method.Secondly,the integer ambiguity is fixed according to damped LAMBDA method,the damped factor is determined by the initial observation result.Finally,the actual GPS data in mining area is used to conduct experiment so as to test the performance of the algorithm proposed in this paper.The experimental results show that,three directions accuracy of the station is in millimeter-level,and calculating time of single epoch is less than the sampling.The algorithm proposed can satisfy the millimeter-level of observation and observation in real time,and provide a certain reference value for the real-time mining settlement observation.

    Mining settlement,GPS single epoch,Tikhonov regularization method,Damped LAMBDA method,Damping factor

    2014-12-22

    “十二五”國(guó)家科技支撐計(jì)劃項(xiàng)目(編號(hào):2011BAB04B10),國(guó)家自然科學(xué)基金項(xiàng)目(編號(hào):40972084),教育部博士點(diǎn)基金項(xiàng)目(編號(hào):2010512211013)。

    陳西斌(1987—),男,碩士研究生。

    TD176,P228.4

    A

    1001-1250(2015)-02-103-05

    猜你喜歡
    歷元正則實(shí)數(shù)
    “實(shí)數(shù)”實(shí)戰(zhàn)操練
    歷元間載波相位差分的GPS/BDS精密單點(diǎn)測(cè)速算法
    剩余有限Minimax可解群的4階正則自同構(gòu)
    類似于VNL環(huán)的環(huán)
    認(rèn)識(shí)實(shí)數(shù)
    1.1 實(shí)數(shù)
    Recent advances of TCM treatment of childhood atopic dermatitis
    比較實(shí)數(shù)的大小
    Clinical observation of Huatan Huoxue Formula in treating coronary heart disease with hyperlipidemia
    Mechanism of sex hormone level in biological clock disorder induced acne and analysis of TCM Pathogenesis
    国内毛片毛片毛片毛片毛片| 一区福利在线观看| 日本-黄色视频高清免费观看| 在线天堂最新版资源| 麻豆久久精品国产亚洲av| 日韩人妻高清精品专区| 深夜a级毛片| 少妇丰满av| 一进一出好大好爽视频| 韩国av一区二区三区四区| 一个人看视频在线观看www免费| 黄色配什么色好看| av天堂在线播放| 国产精品免费一区二区三区在线| 国产成人影院久久av| 九九爱精品视频在线观看| 搡老熟女国产l中国老女人| 精品人妻偷拍中文字幕| 在线免费观看的www视频| 女人十人毛片免费观看3o分钟| 一级黄片播放器| 国产aⅴ精品一区二区三区波| 看免费成人av毛片| 亚洲国产精品sss在线观看| 国产精品爽爽va在线观看网站| 国产淫片久久久久久久久| av黄色大香蕉| 在线观看66精品国产| 男女啪啪激烈高潮av片| 黄色视频,在线免费观看| 日韩 亚洲 欧美在线| 免费在线观看日本一区| 日韩精品青青久久久久久| 美女高潮喷水抽搐中文字幕| 成人av在线播放网站| 成人美女网站在线观看视频| 国产爱豆传媒在线观看| 久久久久久久午夜电影| 日本撒尿小便嘘嘘汇集6| 麻豆一二三区av精品| 91久久精品电影网| 国产精品爽爽va在线观看网站| 狂野欧美白嫩少妇大欣赏| 联通29元200g的流量卡| 亚洲精品一卡2卡三卡4卡5卡| 久久精品夜夜夜夜夜久久蜜豆| 搡老妇女老女人老熟妇| 免费黄网站久久成人精品| 婷婷亚洲欧美| 久久国产乱子免费精品| 69av精品久久久久久| 女的被弄到高潮叫床怎么办 | 丰满的人妻完整版| 高清毛片免费观看视频网站| 啦啦啦韩国在线观看视频| av在线蜜桃| 国产精品人妻久久久影院| 最后的刺客免费高清国语| 久久久久精品国产欧美久久久| 亚洲专区国产一区二区| 日本精品一区二区三区蜜桃| 在现免费观看毛片| 亚洲中文日韩欧美视频| 久久国产乱子免费精品| 欧美成人一区二区免费高清观看| 天堂网av新在线| 伦理电影大哥的女人| 999久久久精品免费观看国产| 小说图片视频综合网站| 国产一级毛片七仙女欲春2| 中文字幕av成人在线电影| 午夜a级毛片| 亚洲无线在线观看| 国产精品久久久久久av不卡| 久久国内精品自在自线图片| 无人区码免费观看不卡| 中文字幕免费在线视频6| 中文字幕高清在线视频| 欧美区成人在线视频| av中文乱码字幕在线| 香蕉av资源在线| a在线观看视频网站| 亚洲国产精品久久男人天堂| 免费在线观看影片大全网站| 极品教师在线免费播放| 亚洲国产精品久久男人天堂| 一进一出抽搐gif免费好疼| 亚洲熟妇熟女久久| 美女高潮的动态| 国产成人一区二区在线| 亚洲av.av天堂| 禁无遮挡网站| 日韩欧美在线乱码| 久久亚洲精品不卡| 日本一二三区视频观看| 午夜免费男女啪啪视频观看 | 久久人人爽人人爽人人片va| 麻豆av噜噜一区二区三区| 欧美日韩国产亚洲二区| 国产精品久久久久久精品电影| 久久久久国内视频| a级一级毛片免费在线观看| 在线观看av片永久免费下载| 嫩草影院新地址| 午夜福利成人在线免费观看| 小说图片视频综合网站| a级毛片a级免费在线| 亚洲人与动物交配视频| 少妇熟女aⅴ在线视频| avwww免费| 午夜老司机福利剧场| 日本一本二区三区精品| 日本五十路高清| 精品免费久久久久久久清纯| 国产综合懂色| 欧美一区二区精品小视频在线| 亚洲精品国产成人久久av| 国产亚洲欧美98| 国内精品一区二区在线观看| 波多野结衣高清作品| 别揉我奶头~嗯~啊~动态视频| 看片在线看免费视频| 久久欧美精品欧美久久欧美| 久久亚洲真实| 免费人成在线观看视频色| 又粗又爽又猛毛片免费看| 偷拍熟女少妇极品色| 91麻豆av在线| 色哟哟哟哟哟哟| 99在线人妻在线中文字幕| 日韩 亚洲 欧美在线| 国产伦在线观看视频一区| eeuss影院久久| 午夜福利在线观看吧| 午夜精品一区二区三区免费看| 老熟妇仑乱视频hdxx| 国产精品精品国产色婷婷| 91精品国产九色| 在线a可以看的网站| 亚洲av免费在线观看| 亚洲人与动物交配视频| 看十八女毛片水多多多| 亚洲av日韩精品久久久久久密| 精品无人区乱码1区二区| 两性午夜刺激爽爽歪歪视频在线观看| 在线免费观看的www视频| 国产美女午夜福利| 精品福利观看| 亚洲精品粉嫩美女一区| 亚洲欧美精品综合久久99| 一级黄色大片毛片| 亚洲精品国产成人久久av| 婷婷丁香在线五月| ponron亚洲| 国内毛片毛片毛片毛片毛片| 精品一区二区三区av网在线观看| 午夜影院日韩av| 亚洲欧美日韩卡通动漫| 99热精品在线国产| 热99在线观看视频| 国产精品美女特级片免费视频播放器| av天堂中文字幕网| 亚洲精品色激情综合| 美女 人体艺术 gogo| 舔av片在线| 在现免费观看毛片| 黄色配什么色好看| 丰满乱子伦码专区| 又爽又黄a免费视频| 日韩人妻高清精品专区| 欧美黑人巨大hd| 亚洲在线自拍视频| 欧美bdsm另类| 中文字幕精品亚洲无线码一区| 伦精品一区二区三区| 国产精品av视频在线免费观看| 国产真实乱freesex| 两性午夜刺激爽爽歪歪视频在线观看| 少妇的逼好多水| 桃色一区二区三区在线观看| 午夜福利在线观看免费完整高清在 | 成人永久免费在线观看视频| 91麻豆精品激情在线观看国产| 日韩 亚洲 欧美在线| 欧美日韩瑟瑟在线播放| 琪琪午夜伦伦电影理论片6080| 久久草成人影院| 欧美三级亚洲精品| 最近中文字幕高清免费大全6 | 久久人人精品亚洲av| 欧美成人a在线观看| 可以在线观看毛片的网站| 一级a爱片免费观看的视频| 久久久久国内视频| 一个人看视频在线观看www免费| 亚洲一级一片aⅴ在线观看| 无人区码免费观看不卡| 亚洲人成网站在线播放欧美日韩| 在线观看美女被高潮喷水网站| 久久精品国产鲁丝片午夜精品 | 韩国av在线不卡| 国产精品一区二区三区四区免费观看 | 成人一区二区视频在线观看| 搡老熟女国产l中国老女人| 久久精品国产亚洲av涩爱 | 干丝袜人妻中文字幕| 午夜日韩欧美国产| 男女啪啪激烈高潮av片| 久久久久久久午夜电影| 51国产日韩欧美| 国产精品电影一区二区三区| 悠悠久久av| 五月玫瑰六月丁香| 午夜精品在线福利| 99热这里只有是精品50| 一区二区三区激情视频| 欧美激情在线99| 日本撒尿小便嘘嘘汇集6| 国产精品乱码一区二三区的特点| 国国产精品蜜臀av免费| 中文在线观看免费www的网站| 亚洲av中文av极速乱 | 两个人的视频大全免费| 毛片一级片免费看久久久久 | 亚洲一级一片aⅴ在线观看| 成年免费大片在线观看| 精品人妻偷拍中文字幕| 成年人黄色毛片网站| 国产探花极品一区二区| 变态另类丝袜制服| 国产大屁股一区二区在线视频| 99国产极品粉嫩在线观看| 亚洲va在线va天堂va国产| 日日摸夜夜添夜夜添av毛片 | 亚洲国产日韩欧美精品在线观看| 日韩欧美三级三区| 蜜桃久久精品国产亚洲av| 成人国产一区最新在线观看| 男人和女人高潮做爰伦理| 熟妇人妻久久中文字幕3abv| 精品久久久久久久久av| 亚洲图色成人| 日韩精品青青久久久久久| 俺也久久电影网| 熟妇人妻久久中文字幕3abv| 欧美日韩中文字幕国产精品一区二区三区| 又黄又爽又刺激的免费视频.| 精品久久久久久,| 久久久色成人| 亚洲人成网站在线播放欧美日韩| 国产精品久久久久久av不卡| 看免费成人av毛片| 免费观看人在逋| 啦啦啦观看免费观看视频高清| 久久亚洲精品不卡| 男人和女人高潮做爰伦理| 亚洲成av人片在线播放无| 五月伊人婷婷丁香| 男插女下体视频免费在线播放| 大又大粗又爽又黄少妇毛片口| 人人妻人人澡欧美一区二区| 亚洲精品影视一区二区三区av| 久久久精品大字幕| 午夜福利18| 亚洲精品影视一区二区三区av| 亚洲av中文字字幕乱码综合| 日韩中字成人| 99国产精品一区二区蜜桃av| 啦啦啦观看免费观看视频高清| 日韩欧美精品v在线| 嫩草影院入口| 亚洲精品乱码久久久v下载方式| 变态另类成人亚洲欧美熟女| 丰满的人妻完整版| 嫁个100分男人电影在线观看| 日韩大尺度精品在线看网址| 国产三级中文精品| 我要看日韩黄色一级片| 特大巨黑吊av在线直播| 在线观看66精品国产| 88av欧美| 能在线免费观看的黄片| 欧美人与善性xxx| 亚洲美女视频黄频| 春色校园在线视频观看| 日韩大尺度精品在线看网址| 色在线成人网| 国产成人影院久久av| 中文字幕av在线有码专区| 亚洲精品在线观看二区| 国产亚洲av嫩草精品影院| 国产精品国产高清国产av| 最近视频中文字幕2019在线8| 免费在线观看日本一区| 十八禁网站免费在线| 天堂影院成人在线观看| 成人国产麻豆网| aaaaa片日本免费| 色哟哟哟哟哟哟| 国模一区二区三区四区视频| 国产毛片a区久久久久| av在线天堂中文字幕| 可以在线观看毛片的网站| 在线免费观看的www视频| 99热这里只有精品一区| 国内揄拍国产精品人妻在线| 成人鲁丝片一二三区免费| 22中文网久久字幕| 午夜爱爱视频在线播放| 悠悠久久av| 男插女下体视频免费在线播放| 欧美激情久久久久久爽电影| 一区二区三区四区激情视频 | 国产精品电影一区二区三区| 九色成人免费人妻av| 国产又黄又爽又无遮挡在线| 丰满的人妻完整版| 久久久久精品国产欧美久久久| www日本黄色视频网| 国内精品久久久久久久电影| 久久久午夜欧美精品| 别揉我奶头~嗯~啊~动态视频| 日本-黄色视频高清免费观看| 午夜福利在线观看吧| 伦理电影大哥的女人| 搡老妇女老女人老熟妇| 特级一级黄色大片| 免费观看人在逋| 国产精品电影一区二区三区| 日韩欧美精品v在线| 日本五十路高清| 美女免费视频网站| 69av精品久久久久久| 成熟少妇高潮喷水视频| 成人特级黄色片久久久久久久| 日韩强制内射视频| 午夜激情福利司机影院| 黄色配什么色好看| 亚洲熟妇中文字幕五十中出| 性欧美人与动物交配| 欧美区成人在线视频| a级毛片a级免费在线| 精品久久久久久久久久久久久| 在线观看av片永久免费下载| 丝袜美腿在线中文| 亚洲三级黄色毛片| 日韩中字成人| 久久久国产成人免费| 看免费成人av毛片| 天堂√8在线中文| 校园春色视频在线观看| ponron亚洲| 色播亚洲综合网| 赤兔流量卡办理| 久久久久久久久大av| 中文字幕精品亚洲无线码一区| 免费看光身美女| 亚洲一级一片aⅴ在线观看| 一边摸一边抽搐一进一小说| 我要搜黄色片| 干丝袜人妻中文字幕| 精品不卡国产一区二区三区| 精品久久国产蜜桃| 美女高潮的动态| 国产午夜福利久久久久久| 欧美性猛交╳xxx乱大交人| 亚洲美女搞黄在线观看 | h日本视频在线播放| 有码 亚洲区| 午夜福利成人在线免费观看| 12—13女人毛片做爰片一| 人妻制服诱惑在线中文字幕| 欧美bdsm另类| 三级毛片av免费| 亚洲中文字幕日韩| 丰满人妻一区二区三区视频av| 在线观看美女被高潮喷水网站| 搡老熟女国产l中国老女人| 五月玫瑰六月丁香| 别揉我奶头 嗯啊视频| 亚洲真实伦在线观看| 国内精品久久久久久久电影| 国产 一区精品| 欧美成人免费av一区二区三区| 3wmmmm亚洲av在线观看| 亚洲美女搞黄在线观看 | 午夜福利在线观看吧| 三级男女做爰猛烈吃奶摸视频| 亚洲精品影视一区二区三区av| 免费av不卡在线播放| 亚洲精品影视一区二区三区av| 变态另类成人亚洲欧美熟女| 亚洲av免费高清在线观看| 两个人的视频大全免费| 亚洲人成伊人成综合网2020| 无遮挡黄片免费观看| 如何舔出高潮| 一进一出好大好爽视频| 中文字幕久久专区| 狂野欧美激情性xxxx在线观看| 午夜免费男女啪啪视频观看 | 亚洲欧美精品综合久久99| 男人舔女人下体高潮全视频| 国内久久婷婷六月综合欲色啪| 色综合亚洲欧美另类图片| 久久国产精品人妻蜜桃| 色哟哟·www| 亚洲午夜理论影院| 狠狠狠狠99中文字幕| 精品福利观看| 久久九九热精品免费| 亚洲成人久久性| 午夜激情福利司机影院| 亚洲成人中文字幕在线播放| 国产精品永久免费网站| 欧美精品国产亚洲| 成人精品一区二区免费| 嫩草影院新地址| 国产黄片美女视频| 久久国内精品自在自线图片| 久久精品人妻少妇| 久久精品国产99精品国产亚洲性色| 亚洲精品成人久久久久久| 久久精品国产亚洲av天美| 国产91精品成人一区二区三区| 国产成人a区在线观看| 哪里可以看免费的av片| 91狼人影院| 日日干狠狠操夜夜爽| 在线观看66精品国产| 久久九九热精品免费| 在线观看av片永久免费下载| 一级av片app| 精品国产三级普通话版| 久久久精品大字幕| 亚洲性夜色夜夜综合| 人妻制服诱惑在线中文字幕| 一本久久中文字幕| 国产乱人视频| 如何舔出高潮| 国产精品国产三级国产av玫瑰| 欧美最黄视频在线播放免费| 深爱激情五月婷婷| 欧美最新免费一区二区三区| 在现免费观看毛片| 免费看a级黄色片| 蜜桃亚洲精品一区二区三区| 久久99热这里只有精品18| 国产精品一区二区免费欧美| 狂野欧美激情性xxxx在线观看| 两性午夜刺激爽爽歪歪视频在线观看| 天美传媒精品一区二区| 亚洲在线自拍视频| 97超视频在线观看视频| 日本三级黄在线观看| 日日撸夜夜添| 国产激情偷乱视频一区二区| 色哟哟·www| 99久久精品一区二区三区| 波多野结衣高清作品| 精品久久久久久久久久久久久| 国产精品无大码| 性色avwww在线观看| 国产毛片a区久久久久| 我要搜黄色片| 日韩欧美在线乱码| 亚洲无线在线观看| 日日夜夜操网爽| 亚洲国产欧洲综合997久久,| 天堂网av新在线| 一进一出好大好爽视频| 国内久久婷婷六月综合欲色啪| 男女做爰动态图高潮gif福利片| 欧美另类亚洲清纯唯美| 少妇人妻一区二区三区视频| 麻豆av噜噜一区二区三区| 欧美三级亚洲精品| 国产精品久久久久久av不卡| 国产探花极品一区二区| 淫秽高清视频在线观看| 又黄又爽又刺激的免费视频.| 久久99热这里只有精品18| 一本精品99久久精品77| 色吧在线观看| 亚洲精品456在线播放app | 国产亚洲精品久久久com| 两人在一起打扑克的视频| 99国产极品粉嫩在线观看| 狂野欧美白嫩少妇大欣赏| 日韩精品青青久久久久久| 亚洲成a人片在线一区二区| 91av网一区二区| 91精品国产九色| 久久精品夜夜夜夜夜久久蜜豆| 国产久久久一区二区三区| 亚洲美女视频黄频| 亚洲精品乱码久久久v下载方式| 日韩中字成人| 欧美日韩亚洲国产一区二区在线观看| 久久人人精品亚洲av| 国产精品伦人一区二区| 婷婷色综合大香蕉| 极品教师在线视频| 婷婷色综合大香蕉| 蜜桃久久精品国产亚洲av| 成人国产麻豆网| 在线播放国产精品三级| 日韩欧美一区二区三区在线观看| 天堂影院成人在线观看| 亚洲一级一片aⅴ在线观看| 精品久久久久久成人av| 搡老岳熟女国产| 老司机福利观看| 亚洲成人精品中文字幕电影| 97超级碰碰碰精品色视频在线观看| 此物有八面人人有两片| 久9热在线精品视频| 精品一区二区三区人妻视频| 免费观看的影片在线观看| av天堂中文字幕网| 国产黄a三级三级三级人| av在线老鸭窝| 日韩欧美三级三区| 日本黄色片子视频| 神马国产精品三级电影在线观看| 精品无人区乱码1区二区| 不卡一级毛片| 午夜精品在线福利| 国产女主播在线喷水免费视频网站 | 日韩欧美在线乱码| 麻豆国产97在线/欧美| 国产 一区 欧美 日韩| 九九久久精品国产亚洲av麻豆| 韩国av在线不卡| 亚洲精品一区av在线观看| 少妇人妻精品综合一区二区 | 两个人视频免费观看高清| 内射极品少妇av片p| 国产亚洲91精品色在线| 尤物成人国产欧美一区二区三区| 国产免费av片在线观看野外av| 精品午夜福利视频在线观看一区| 国产午夜精品久久久久久一区二区三区 | 久久久久性生活片| 国产精品综合久久久久久久免费| av女优亚洲男人天堂| 在线天堂最新版资源| 日韩欧美精品免费久久| 国产免费男女视频| 观看美女的网站| 久久久久免费精品人妻一区二区| 最近视频中文字幕2019在线8| 黄色欧美视频在线观看| 国产欧美日韩精品亚洲av| 成人欧美大片| 长腿黑丝高跟| 国产精品电影一区二区三区| 麻豆国产97在线/欧美| 国产亚洲精品综合一区在线观看| 狂野欧美激情性xxxx在线观看| 亚洲国产日韩欧美精品在线观看| 在线观看美女被高潮喷水网站| 国产精品乱码一区二三区的特点| 亚洲av熟女| 亚洲va日本ⅴa欧美va伊人久久| 综合色av麻豆| 亚洲三级黄色毛片| 国产亚洲91精品色在线| 久久久久久久久中文| 亚洲无线在线观看| 免费av观看视频| 岛国在线免费视频观看| 国产精品98久久久久久宅男小说| 春色校园在线视频观看| av.在线天堂| 亚洲精华国产精华精| 色综合色国产| 久久天躁狠狠躁夜夜2o2o| 国产午夜福利久久久久久| 国产精品嫩草影院av在线观看 | а√天堂www在线а√下载| av天堂中文字幕网| 欧美色视频一区免费| 久久久午夜欧美精品| 成年女人毛片免费观看观看9| 精品久久久久久久久亚洲 | 久久精品国产亚洲av香蕉五月| 亚洲美女搞黄在线观看 | 女生性感内裤真人,穿戴方法视频| 真人做人爱边吃奶动态| 国产精品女同一区二区软件 | av视频在线观看入口| 变态另类成人亚洲欧美熟女| 日韩亚洲欧美综合| 99国产极品粉嫩在线观看| 一级av片app| 黄色女人牲交| 99精品在免费线老司机午夜| 噜噜噜噜噜久久久久久91| 国内精品一区二区在线观看| 在线观看午夜福利视频| 亚洲一区二区三区色噜噜| 国内精品久久久久久久电影| 少妇丰满av| av在线老鸭窝| 国产欧美日韩精品亚洲av| 午夜免费男女啪啪视频观看 | 欧美日韩乱码在线| 欧美不卡视频在线免费观看| videossex国产| 亚洲成人久久性| 精品无人区乱码1区二区| 欧美日韩精品成人综合77777| av黄色大香蕉| 哪里可以看免费的av片|