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

    地下目標(biāo)體重力張量數(shù)據(jù)的邊界識(shí)別與增強(qiáng)算法的研究

    2017-02-09 03:08:26王慶賓黃佳喜馮進(jìn)凱
    測(cè)繪工程 2017年4期
    關(guān)鍵詞:效果水平

    張 超,王慶賓,黃佳喜,馮進(jìn)凱

    (信息工程大學(xué) 地理空間信息學(xué)院,河南 鄭州 450052)

    地下目標(biāo)體重力張量數(shù)據(jù)的邊界識(shí)別與增強(qiáng)算法的研究

    張 超,王慶賓,黃佳喜,馮進(jìn)凱

    (信息工程大學(xué) 地理空間信息學(xué)院,河南 鄭州 450052)

    邊界增強(qiáng)與識(shí)別在重力數(shù)據(jù)處理中占據(jù)重要地位,與傳統(tǒng)重力異常數(shù)據(jù)相比,重力張量及其高階分量對(duì)于直接反映異常體的邊界具有更高的精度。當(dāng)異常數(shù)據(jù)中的所有網(wǎng)格點(diǎn)的值均較低時(shí),通過Sigmoid變換,可以實(shí)現(xiàn)高異常值網(wǎng)格數(shù)據(jù)的拉升,同時(shí)壓縮低灰度級(jí)像素,從而凸顯地質(zhì)體邊界,提高邊界增強(qiáng)后圖像的識(shí)別效果。文中利用張量及其分量構(gòu)建常用的邊界識(shí)別算法,通過組合體模型進(jìn)行多種邊界識(shí)別算法的試算,以比較分析各自的效果,并對(duì)結(jié)果進(jìn)行Sigmoid變換。結(jié)果表明:對(duì)于張量高階分量組合形式,水平梯度模、解析信號(hào)能基本反映淺異常體的邊界,gzz水平梯度模能較好反映淺異常體邊界,但三者均不能識(shí)別深異常體邊界;Tilt梯度、Theta和ITA3效果不理想;ITA2能在有效均衡不同強(qiáng)度異常信號(hào)的同時(shí),清晰地識(shí)別不同深度異常體的邊界;采用Sigmoid變換,可以明顯提高邊界識(shí)別的識(shí)別效果。

    重力張量;邊界識(shí)別;水平梯度;解析信號(hào);Sigmoid變換

    重力勘探方法屬于應(yīng)用地球物理地質(zhì)勘探與礦藏尋找技術(shù)的一個(gè)重要分支。由于礦產(chǎn)資源或地下地質(zhì)構(gòu)造的密度以及磁性一般同附近的巖體有著顯著的不同,因此,重力勘探在尋找礦產(chǎn)資源與地質(zhì)構(gòu)造研究中具有重要的應(yīng)用價(jià)值。重力勘探具有快速、經(jīng)濟(jì)、范圍大等優(yōu)點(diǎn),是進(jìn)行構(gòu)造劃分、異常圈定和發(fā)掘礦產(chǎn)資源等快速而有效的手段[1]。地質(zhì)體邊界是地質(zhì)體的斷裂帶、相異地質(zhì)體之間的分界線。在這些位置處,由于重力能在一定程度上表現(xiàn)出數(shù)值的不連續(xù)性,呈現(xiàn)出重力異常的較大變化。利用這一特點(diǎn),可以通過重力異常的較大變化率尋找對(duì)應(yīng)的地質(zhì)體邊界線,從而為地下地質(zhì)結(jié)構(gòu)的分析提供支持。

    邊界增強(qiáng)與識(shí)別在重力數(shù)據(jù)處理中占據(jù)重要地位。原始實(shí)測(cè)重力異常數(shù)據(jù)對(duì)地下異常體信息反映并不突出,直接根據(jù)此數(shù)據(jù)并不能準(zhǔn)確劃分出地質(zhì)體的邊界。因此需要對(duì)原始異常數(shù)據(jù)進(jìn)行相應(yīng)濾波處理來增強(qiáng)對(duì)地質(zhì)體邊界的識(shí)別。重力異常水平梯度的極大值和垂直梯度的零值基本反映地質(zhì)體邊界,邊界提取算法也是基于此提出的[2]。常見的位場(chǎng)數(shù)據(jù)邊界提取算法主要包括Pedersen等提出的兩個(gè)張量不變量組合、水平梯度模法(The total horizontal derivatives,THDR)[3]、解析信號(hào)法(The analytic signal amplitude,ASA)[4]、Tilt梯度法[5]、Tilt梯度的水平導(dǎo)數(shù)法[6]、Theta圖法[7]、gzz的水平梯度模法等等,然后利用圖中的極值點(diǎn)、零點(diǎn)或其他特征點(diǎn)來識(shí)別邊界。本文利用張量及其高階分量構(gòu)建了常用的邊界識(shí)別算法,通過組合體模型進(jìn)行多種邊界識(shí)別算法的試算,比較分析各自的效果,并對(duì)結(jié)果進(jìn)行Sigmoid變換。

    1 重力異常及梯度張量

    建立一個(gè)局部坐標(biāo)系,以地面某一點(diǎn)O為坐標(biāo)原點(diǎn),坐標(biāo)軸ζ軸沿重力方向垂直向下,ξ,η軸在水準(zhǔn)面上,以下所有的計(jì)算都在局部坐標(biāo)系下進(jìn)行。設(shè)異常體剩余密度為σ,則引力位

    (1)

    重力異常就是引力位沿ζ軸方向的導(dǎo)數(shù),即

    (2)

    (3)

    由于在地球外部,位場(chǎng)滿足Laplace方程Uxx+Uyy+Uzz=0 ,以及重力場(chǎng)的無旋性,即Uxy=Uyx,Uxz=Uzx,Uyz=Uzy,因此重力梯度張量中只有5個(gè)獨(dú)立分量。

    2 邊界增強(qiáng)與識(shí)別算法

    2.1 邊界識(shí)別算法

    常用的異常數(shù)據(jù)處理濾波器主要都是基于水平梯度和垂直導(dǎo)數(shù)的各種組合來構(gòu)建,以增強(qiáng)異常地質(zhì)體的邊界特征,從而實(shí)現(xiàn)邊界識(shí)別。與傳統(tǒng)重力異常數(shù)據(jù)相比,重力張量數(shù)據(jù)對(duì)于直接反映異常體的邊界具有更高的精度。通過梯度張量及其高階分量的不同組合,把多個(gè)分量的信息集中增強(qiáng),突出邊界信息,從而輔助數(shù)據(jù)處理和解釋。

    對(duì)于各張量組合形式的表達(dá)式見表1和表2。

    表1 各張量組合形式的表達(dá)式

    表2 常見張量高階組合形式的表達(dá)式

    2.2 重力張量數(shù)據(jù)對(duì)比度增強(qiáng)

    采用對(duì)比度增強(qiáng)方法對(duì)重力張量網(wǎng)格數(shù)據(jù)進(jìn)行對(duì)比度拉伸,增強(qiáng)數(shù)據(jù)的對(duì)比度,擴(kuò)展其中的強(qiáng)異常與弱異常區(qū)域的對(duì)比效果,可以增強(qiáng)邊界識(shí)別和檢測(cè)的效果。這里主要運(yùn)用Sigmoid變換對(duì)網(wǎng)格值進(jìn)行處理。

    Sigmoid變換方程式原型為

    (4)

    Sigmoid曲線呈現(xiàn)為S狀,通過將其作用于異常數(shù)據(jù),可以實(shí)現(xiàn)高異常值網(wǎng)格數(shù)據(jù)的拉伸,同時(shí)壓縮低灰度級(jí)像素[8]。這在一定意義上符合異常體邊界初步增強(qiáng)后,進(jìn)一步提高邊界對(duì)比度的要求,從而凸顯地質(zhì)體的邊界,提高邊界增強(qiáng)后圖像的視覺效果。通過引入合適的參數(shù),可以在一定程度上改變響應(yīng)曲線的形狀,從而使變換算法更加實(shí)用。引入?yún)?shù)后的模型為

    (5)

    式中:fmax為異常數(shù)據(jù)最大值;I(i,j) 為原位于網(wǎng)格(i,j)的異常值;f(i,j)為網(wǎng)格(i,j)處理結(jié)果異常值;α,β為調(diào)整變換模型的參數(shù)。其中α為用于控制形狀的參數(shù),β則用于控制相位。在實(shí)際應(yīng)用時(shí),可以根據(jù)當(dāng)前處理的異常數(shù)據(jù)的不同特性選擇適當(dāng)?shù)摩?,β值,從而達(dá)到有針對(duì)性的改善異常數(shù)據(jù)分布質(zhì)量,實(shí)現(xiàn)對(duì)比度拉升。參數(shù)中的α值決定變換函數(shù)的陡峭程度,而β值則確定變換函數(shù)在坐標(biāo)水平軸上的平移分量。分析可得:當(dāng)β取值很小時(shí),可以實(shí)現(xiàn)壓縮低異常值部分,拉伸高異常值部分。當(dāng)β取值較大時(shí),則全圖的所有網(wǎng)格點(diǎn)的值均會(huì)得以拉升,其中高異常值的像素的拉升度要強(qiáng)于低灰階的拉升度。由該函數(shù)性質(zhì)特點(diǎn)可知,當(dāng)異常數(shù)據(jù)中的所有網(wǎng)格點(diǎn)的值均較低時(shí),采用這種方式得到的增強(qiáng)效果比較明顯。

    3 模型試算

    為充分驗(yàn)證表2邊界識(shí)別算法對(duì)地下目標(biāo)體的識(shí)別效果,這里通過建立規(guī)則的地下空洞模型進(jìn)行計(jì)算分析。本文所選地面觀測(cè)數(shù)據(jù)區(qū)域?yàn)?01×401點(diǎn)距(點(diǎn)距單位為m)的平面網(wǎng)格面積,并通過規(guī)則模型正演計(jì)算得到重力異常數(shù)據(jù)。本文布設(shè)一個(gè)由球體和長(zhǎng)方體構(gòu)成的組合體模型:其中球體半徑為30 m,中心點(diǎn)坐標(biāo)為(-110,-110),上頂面埋深為30 m,與圍巖密度差為2.67 g/cm2;長(zhǎng)方體的長(zhǎng)寬高為120 m×30 m×30 m,中心點(diǎn)坐標(biāo)為(50,50),上頂面埋深為15 m,與圍巖密度差為2.67 g/cm2。

    3.1 常見張量高階組合形式的邊界識(shí)別效果分析

    分別利用水平梯度模法、解析信號(hào)法、Tilt梯度法、Theta圖法、gzz的水平梯度模法、改進(jìn)二階傾斜角等方法ITA2對(duì)重力異常進(jìn)行處理,結(jié)果見圖1。

    圖1 模型重力張量組合形式邊界識(shí)別效果

    從以上試算結(jié)果中可以看出:水平梯度模和解析信號(hào)能基本反映出淺異常體平面位置和邊界,但基本不能反映深異常體的邊界;雖然Tilt梯度和Theta圖能有效均衡反映不同深度的異常體信號(hào),但Tilt梯度不能反映異常體的邊界,而Theta圖反映的邊界模糊且發(fā)散,并且兩者均出現(xiàn)邊界效應(yīng);gzz的水平梯度模能夠較好地反映出淺異常體的邊界,但不能反映深異常體的邊界;改進(jìn)二階傾斜角ITA2能有效均衡反映不同深度的異常體信號(hào),識(shí)別出的邊界略大于真實(shí)邊界,對(duì)異常體邊界反映效果最好。

    3.2 改進(jìn)二階傾斜角ITA2和改進(jìn)三階傾斜角ITA3對(duì)比

    由3.1中試算結(jié)果可知,ITA2對(duì)不同深度的異常體邊界識(shí)別效果最好。理論上,重力張量高階導(dǎo)數(shù)對(duì)異常體的邊界識(shí)別具有更高的分辨率,能更加突出異常體邊界。這里構(gòu)建ITA3,并比較ITA2和ITA3對(duì)不同深度的異常體模型的邊界識(shí)別效果。

    (6)

    試驗(yàn)?zāi)P团c3.1中模型一樣,結(jié)果如圖2所示。

    圖2 ITA2和ITA3比較結(jié)果

    從圖2可以看出:與ITA2相比,ITA3所識(shí)別的淺異常體的邊界形狀清晰準(zhǔn)確,識(shí)別的深異常體邊界清晰,但卻明顯小于真實(shí)異常體邊界,而且產(chǎn)生較大的噪聲和邊界效應(yīng)。總的來說,識(shí)別效果較差,明顯不如ITA2識(shí)別效果,可能是由于高階導(dǎo)數(shù)的計(jì)算增大了噪聲的干擾。

    3.3 Sigmoid變換分析

    分別對(duì)前文中水平梯度模、解析信號(hào)、Tilt梯度、Theta圖法、gzz的水平梯度模、改進(jìn)二階傾斜角ITA2和ITA3進(jìn)行Sigmoid變換,結(jié)果如圖3所示。

    圖3 Sigmoid變換結(jié)果

    從圖3對(duì)比可以看出:水平梯度模、解析信號(hào)、Tilt梯度、Theta圖法、gzz的水平梯度模、改進(jìn)二階傾斜角ITA2和ITA3等異常數(shù)據(jù)的對(duì)比度得到明顯的提高,凸顯異常體的邊界,提高邊界識(shí)別圖像的識(shí)別效果;但是從(d)和(f)可以看出,Sigmoid變換在增強(qiáng)有用異常信號(hào)的同時(shí)也使得數(shù)據(jù)中的噪聲被增大。

    4 結(jié) 論

    本文通過對(duì)組合體模型試算分析,全面比較多種張量高階組合形式的邊界識(shí)別方法的應(yīng)用效果,并對(duì)各種邊界識(shí)別算法得到的張量異常數(shù)據(jù)進(jìn)行Sigmoid變換以增強(qiáng)數(shù)據(jù)對(duì)比度。分析結(jié)果表明:水平梯度模和解析信號(hào)能基本反映出淺異常體平面位置和邊界,但不能反映出深異常體的邊界;Tilt梯度和Theta圖能有效均衡反映不同深度的異常體信號(hào),但Tilt梯度不能反映異常體的邊界,而Theta圖反映的邊界模糊且發(fā)散,并且兩者均出現(xiàn)邊界效應(yīng);gzz的水平梯度模能夠較好地反映出淺異常體的邊界,但不能反映深異常體的邊界;ITA3所識(shí)別的淺異常體的邊界形狀清晰準(zhǔn)確,但由于高階導(dǎo)數(shù)的計(jì)算增大噪聲的干擾,從而產(chǎn)生較大的噪聲和邊界效應(yīng);改進(jìn)二階傾斜角ITA2能有效均衡反映不同深度的異常體信號(hào),識(shí)別出的邊界略大于真實(shí)邊界,對(duì)異常體邊界反映效果最好;當(dāng)異常數(shù)據(jù)中的所有網(wǎng)格點(diǎn)的值均較低時(shí),采用Sigmoid變換可以凸顯異常體的邊界,提高邊界圖像的識(shí)別效果。ITA2和Sigmoid變換能夠在異常體的邊界識(shí)別工作中發(fā)揮更大的作用,從而促進(jìn)重力反演解釋工作的順利進(jìn)行。

    [1] 馬國(guó)慶.位場(chǎng)(重&磁)及其梯度異常自動(dòng)解釋方法研究[D].長(zhǎng)春:吉林大學(xué),2013.

    [2] 吳亮,劉長(zhǎng)弘,王慶賓,等.重力異常對(duì)地下異常體邊界的識(shí)別[J].測(cè)繪科學(xué),2015,40(12):257-262.

    [3] 馬國(guó)慶,李麗麗,杜曉娟.磁張量數(shù)據(jù)的邊界識(shí)別和解釋方法[J].石油地球物理勘探,2012,47(5):815-821.

    [4] NABIGHIAN M N.The analytic signal of two dimensional magnetic bodies with polygonal cross section:its properties and use for automated anomaly interpretation [J].Geophysics,2012,37(3):507.

    [5] 王明,郭志宏,駱遙,等.Tilt-Euler方法在位場(chǎng)數(shù)據(jù)處理及解釋中的應(yīng)用[J].物探與化探,2012(2):126-132.

    [6] VERDUZCO B,F(xiàn)AIRHEAD J D,GREEN C M,et al.The meter reader-new insights into magnetic derivatives for structural mapping [J],Leading Edge,2012 (2):116.

    [7] WIJNS C,PEREZ C,KOWALCZYK P.Theta map:edge detection in magnetic data[J].Geophysics,2005,70(4):39-43.

    [8] 高云軍,劉召芹,萬文輝,等.DEM網(wǎng)格間距對(duì)重力近區(qū)地形改正精度影響驗(yàn)證[J].測(cè)繪工程,2016,25(8):14-18.

    [9] 賈魯,王世忠,王慶賓,等.基于Android平臺(tái)的兩種重力垂直梯度解算方法研究[J].測(cè)繪工程,2016,25(7):57-61.

    [責(zé)任編輯:張德福]

    Research on underground anomaly boundary recognition and enhancement algorithms based on gravitational tensor

    ZHANG Chao, WANG Qingbin, HUANG Jiaxi, FENG Jinkai

    (School of Surveying and Mapping, Information Engineering University, Zhengzhou 450052,China)

    Boundary enhancement and recognition occupies the important position in magnetic data processing. Compared with the traditional gravity anomaly data, gravitational tensor data has a higher accuracy in directly reflecting the boundary of the anomaly. High abnormal values can be pulled up with low gray-scale pixels compressed at the same time by Sigmoid transformation, when all the values of the grid points are lower, thus highlighting geological boundary and enhancing image recognition effect. This paper calculates many kinds of gravity anomaly processing algorithms based on tensor components, including Sigmoid transformation, and then analyzes their effects through combination model. The result shows that: the horizontal gradient,analytic signal and the horizontal gradient of,gzzcan reflect the boundary of the shallow anomaly with the horizontal gradient of,gzza better effect, but none of them can recognize the boundary of deep anomaly. The result of Tilt gradient and Theta diagram method is not ideal, so it is with the ITA3. By contrast, the ITA2can effectively balance the different intensity anomaly and at the same time clearly identify the different depth of anomaly boundary. Sigmoid transformation is proposed, which can distinctly improve the recognition effect of boundary identification.

    gravitational tensor;boundary recognition;horizontal gradient;analytical signal;sigmoid transformation

    引用著錄:張 超,王慶賓,黃佳喜,等.地下目標(biāo)體重力張量數(shù)據(jù)的邊界識(shí)別與增強(qiáng)算法的研究[J].測(cè)繪工程,2017,26(4):16-21.

    10.19349/j.cnki.issn1006-7949.2017.04.004

    2016-06-09

    973國(guó)家高科技計(jì)劃資助項(xiàng)目(6132220202)

    張 超(1990-),男,碩士研究生.

    P223

    A

    1006-7949(2017)04-0016-06

    猜你喜歡
    效果水平
    張水平作品
    按摩效果確有理論依據(jù)
    作家葛水平
    火花(2019年12期)2019-12-26 01:00:28
    加強(qiáng)上下聯(lián)動(dòng) 提升人大履職水平
    迅速制造慢門虛化效果
    抓住“瞬間性”效果
    中華詩詞(2018年11期)2018-03-26 06:41:34
    模擬百種唇妝效果
    Coco薇(2016年8期)2016-10-09 02:11:50
    老虎獻(xiàn)臀
    3D—DSA與3D—CTA成像在顱內(nèi)動(dòng)脈瘤早期診斷中的應(yīng)用效果比較
    組合練習(xí)難度大,貼近實(shí)戰(zhàn)效果佳
    亚洲精品成人久久久久久| 蜜臀久久99精品久久宅男| 国产探花极品一区二区| 亚洲无线观看免费| 亚洲国产精品国产精品| 午夜福利成人在线免费观看| 麻豆成人午夜福利视频| 精品久久久久久久末码| 亚洲婷婷狠狠爱综合网| 欧美性感艳星| 日本精品一区二区三区蜜桃| 十八禁网站免费在线| 国产精品久久视频播放| 精品人妻一区二区三区麻豆 | 村上凉子中文字幕在线| 欧美日韩乱码在线| 成年免费大片在线观看| av在线亚洲专区| 亚洲激情五月婷婷啪啪| 久久久久久国产a免费观看| 免费观看的影片在线观看| 成人亚洲精品av一区二区| 色5月婷婷丁香| 久久欧美精品欧美久久欧美| 淫秽高清视频在线观看| 久久99热6这里只有精品| 真人做人爱边吃奶动态| 蜜臀久久99精品久久宅男| 色综合站精品国产| 国产视频内射| 国产免费男女视频| 国产不卡一卡二| 综合色丁香网| 成熟少妇高潮喷水视频| 伊人久久精品亚洲午夜| 搡老岳熟女国产| 天美传媒精品一区二区| 亚洲av电影不卡..在线观看| .国产精品久久| 看非洲黑人一级黄片| 少妇裸体淫交视频免费看高清| aaaaa片日本免费| 国产欧美日韩精品一区二区| 久久精品人妻少妇| 综合色丁香网| 精品99又大又爽又粗少妇毛片| 熟妇人妻久久中文字幕3abv| 国产真实伦视频高清在线观看| 99九九线精品视频在线观看视频| 干丝袜人妻中文字幕| 中国美白少妇内射xxxbb| 男人舔奶头视频| 国产精品三级大全| 一进一出抽搐动态| 国产伦在线观看视频一区| 国产男靠女视频免费网站| a级毛片免费高清观看在线播放| 国产人妻一区二区三区在| 人妻丰满熟妇av一区二区三区| 免费观看在线日韩| 久久精品国产清高在天天线| 校园人妻丝袜中文字幕| 日本a在线网址| 变态另类丝袜制服| 成人漫画全彩无遮挡| 国产一区二区激情短视频| 色哟哟·www| 长腿黑丝高跟| 国产男人的电影天堂91| 亚洲国产精品sss在线观看| 亚洲婷婷狠狠爱综合网| 亚洲欧美日韩卡通动漫| 日本免费a在线| 一区二区三区免费毛片| 美女cb高潮喷水在线观看| 亚洲aⅴ乱码一区二区在线播放| 激情 狠狠 欧美| 亚洲精品粉嫩美女一区| 变态另类丝袜制服| 丰满人妻一区二区三区视频av| 亚洲美女视频黄频| av国产免费在线观看| 美女大奶头视频| 久久精品国产自在天天线| 色噜噜av男人的天堂激情| 久久精品久久久久久噜噜老黄 | 不卡一级毛片| 嫩草影院精品99| 麻豆精品久久久久久蜜桃| 欧美激情国产日韩精品一区| 丰满的人妻完整版| 午夜爱爱视频在线播放| 日韩成人伦理影院| 少妇人妻一区二区三区视频| 女人十人毛片免费观看3o分钟| 精品日产1卡2卡| 一个人看的www免费观看视频| 99久久精品热视频| 99在线视频只有这里精品首页| 欧美潮喷喷水| 欧美日韩在线观看h| 搡老妇女老女人老熟妇| 直男gayav资源| 国产熟女欧美一区二区| 亚洲无线在线观看| 欧美另类亚洲清纯唯美| 变态另类成人亚洲欧美熟女| 亚洲av美国av| 日本撒尿小便嘘嘘汇集6| 亚洲性久久影院| 91av网一区二区| 99国产精品一区二区蜜桃av| 丝袜喷水一区| 免费看a级黄色片| 免费看美女性在线毛片视频| 国产久久久一区二区三区| 欧美区成人在线视频| 国产精品一区二区三区四区久久| 免费在线观看成人毛片| 国产精品1区2区在线观看.| 香蕉av资源在线| 一本一本综合久久| 欧美激情国产日韩精品一区| 我的女老师完整版在线观看| 狂野欧美白嫩少妇大欣赏| 国产精品不卡视频一区二区| 日本 av在线| 一边摸一边抽搐一进一小说| a级一级毛片免费在线观看| 99国产极品粉嫩在线观看| 最近中文字幕高清免费大全6| 内地一区二区视频在线| 在线天堂最新版资源| 国产男人的电影天堂91| 极品教师在线视频| 国产视频一区二区在线看| 天堂√8在线中文| 麻豆乱淫一区二区| 黄色一级大片看看| 国产真实伦视频高清在线观看| 日韩制服骚丝袜av| 精品少妇黑人巨大在线播放 | or卡值多少钱| 日本欧美国产在线视频| 成年av动漫网址| 午夜免费男女啪啪视频观看 | 亚洲av免费高清在线观看| 日本三级黄在线观看| 国产精品99久久久久久久久| 国产乱人视频| 久久精品影院6| 午夜福利视频1000在线观看| 国产久久久一区二区三区| 亚洲最大成人中文| av在线观看视频网站免费| 久久人人爽人人爽人人片va| 一个人观看的视频www高清免费观看| 久久精品夜色国产| 久久久a久久爽久久v久久| 插逼视频在线观看| 丰满人妻一区二区三区视频av| 男女视频在线观看网站免费| 老司机午夜福利在线观看视频| 波多野结衣巨乳人妻| 身体一侧抽搐| 久久精品久久久久久噜噜老黄 | 国产精品乱码一区二三区的特点| 在线观看av片永久免费下载| 久久婷婷人人爽人人干人人爱| 乱人视频在线观看| 亚洲欧美清纯卡通| 一级av片app| 联通29元200g的流量卡| 国产综合懂色| 中文字幕人妻熟人妻熟丝袜美| 天堂网av新在线| 欧美潮喷喷水| 大香蕉久久网| 我的老师免费观看完整版| 欧美日韩在线观看h| 亚洲一级一片aⅴ在线观看| av在线天堂中文字幕| 日本免费a在线| 好男人在线观看高清免费视频| 麻豆乱淫一区二区| av在线老鸭窝| 九九在线视频观看精品| 国产亚洲91精品色在线| 精品欧美国产一区二区三| 中国美白少妇内射xxxbb| 日韩欧美精品v在线| videossex国产| 国产探花在线观看一区二区| 搡女人真爽免费视频火全软件 | 成人无遮挡网站| 91午夜精品亚洲一区二区三区| 成人特级黄色片久久久久久久| 三级国产精品欧美在线观看| 人人妻,人人澡人人爽秒播| 国产精品久久久久久精品电影| 久久精品综合一区二区三区| 国产精品精品国产色婷婷| 俺也久久电影网| 欧洲精品卡2卡3卡4卡5卡区| 精品一区二区三区av网在线观看| 伊人久久精品亚洲午夜| 插逼视频在线观看| 久久欧美精品欧美久久欧美| 色综合亚洲欧美另类图片| 嫩草影院新地址| 麻豆国产av国片精品| 欧洲精品卡2卡3卡4卡5卡区| 国产精品人妻久久久久久| 三级国产精品欧美在线观看| 国产不卡一卡二| 99久久无色码亚洲精品果冻| 亚洲无线观看免费| 有码 亚洲区| 蜜桃久久精品国产亚洲av| 精品人妻偷拍中文字幕| 男插女下体视频免费在线播放| 亚洲成av人片在线播放无| 日日撸夜夜添| 欧美丝袜亚洲另类| 人妻夜夜爽99麻豆av| 岛国在线免费视频观看| 久久久欧美国产精品| 免费av毛片视频| 秋霞在线观看毛片| 精品熟女少妇av免费看| 亚洲婷婷狠狠爱综合网| 日本熟妇午夜| 欧美性猛交╳xxx乱大交人| 久久久久精品国产欧美久久久| 神马国产精品三级电影在线观看| 国产欧美日韩一区二区精品| 精品不卡国产一区二区三区| 美女xxoo啪啪120秒动态图| 亚洲人成网站高清观看| 日韩制服骚丝袜av| 成人亚洲精品av一区二区| 午夜福利高清视频| 亚洲最大成人手机在线| 人人妻人人澡人人爽人人夜夜 | 在线看三级毛片| 全区人妻精品视频| 亚洲aⅴ乱码一区二区在线播放| 国产激情偷乱视频一区二区| 黑人高潮一二区| 欧美成人免费av一区二区三区| 午夜福利高清视频| 在现免费观看毛片| 又黄又爽又免费观看的视频| 观看美女的网站| 国产探花在线观看一区二区| 国产精品免费一区二区三区在线| 直男gayav资源| 2021天堂中文幕一二区在线观| 99久国产av精品国产电影| 少妇高潮的动态图| 在线观看美女被高潮喷水网站| 亚洲三级黄色毛片| 丰满的人妻完整版| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品日韩av片在线观看| 美女 人体艺术 gogo| 日本黄色视频三级网站网址| 美女内射精品一级片tv| 久久鲁丝午夜福利片| 国国产精品蜜臀av免费| 淫妇啪啪啪对白视频| 日韩精品有码人妻一区| 亚洲欧美精品综合久久99| 欧美一区二区精品小视频在线| 国模一区二区三区四区视频| 国产黄片美女视频| 女人十人毛片免费观看3o分钟| 少妇的逼水好多| 亚洲精品色激情综合| 亚洲美女搞黄在线观看 | 99久久成人亚洲精品观看| 国产精品无大码| 国产精品美女特级片免费视频播放器| 悠悠久久av| 成人无遮挡网站| 看免费成人av毛片| 少妇裸体淫交视频免费看高清| 美女被艹到高潮喷水动态| 久久99热这里只有精品18| 男人狂女人下面高潮的视频| 亚洲乱码一区二区免费版| 能在线免费观看的黄片| 在线天堂最新版资源| 亚洲国产精品成人综合色| 久久亚洲国产成人精品v| 天堂√8在线中文| 亚洲天堂国产精品一区在线| 国产精品av视频在线免费观看| 成人鲁丝片一二三区免费| 男女啪啪激烈高潮av片| 亚洲人与动物交配视频| 国产高清视频在线播放一区| 亚洲欧美日韩高清专用| 亚洲电影在线观看av| 亚洲精品影视一区二区三区av| 国产免费男女视频| 观看美女的网站| 日韩欧美在线乱码| 中文字幕熟女人妻在线| 波多野结衣巨乳人妻| 亚洲av成人av| 啦啦啦观看免费观看视频高清| 久久久久久伊人网av| 中文字幕免费在线视频6| 日韩欧美精品v在线| 午夜激情福利司机影院| 此物有八面人人有两片| 国产片特级美女逼逼视频| 高清午夜精品一区二区三区 | 亚洲第一电影网av| 色综合亚洲欧美另类图片| 欧美极品一区二区三区四区| 色噜噜av男人的天堂激情| 又爽又黄无遮挡网站| 亚洲精品在线观看二区| 99热这里只有是精品在线观看| 悠悠久久av| 毛片女人毛片| 日本黄大片高清| 国产成人影院久久av| 亚洲色图av天堂| 岛国在线免费视频观看| 少妇高潮的动态图| 毛片一级片免费看久久久久| 国产精品1区2区在线观看.| 久久6这里有精品| 国产成人一区二区在线| 免费大片18禁| 亚洲性久久影院| 最近视频中文字幕2019在线8| 看片在线看免费视频| 少妇的逼好多水| 午夜精品一区二区三区免费看| 国产高清激情床上av| 人妻久久中文字幕网| 日本与韩国留学比较| 麻豆精品久久久久久蜜桃| 又粗又爽又猛毛片免费看| 乱系列少妇在线播放| 亚洲一区二区三区色噜噜| 国产精品一区二区性色av| 久久国内精品自在自线图片| 色哟哟哟哟哟哟| 一本一本综合久久| 91久久精品电影网| 在线免费观看不下载黄p国产| 亚洲国产高清在线一区二区三| 亚洲高清免费不卡视频| 精品少妇黑人巨大在线播放 | 精品日产1卡2卡| 免费看光身美女| 久久国产乱子免费精品| 99久久精品国产国产毛片| 人妻丰满熟妇av一区二区三区| 神马国产精品三级电影在线观看| 免费大片18禁| 国产精品亚洲美女久久久| 中文字幕熟女人妻在线| 国产伦一二天堂av在线观看| 69人妻影院| 日本熟妇午夜| 日日干狠狠操夜夜爽| 一个人看视频在线观看www免费| 免费观看人在逋| 国产免费男女视频| 99国产精品一区二区蜜桃av| 国语自产精品视频在线第100页| 亚洲aⅴ乱码一区二区在线播放| 国产精品一区二区三区四区免费观看 | 欧美极品一区二区三区四区| 99久久精品热视频| 日本一本二区三区精品| 床上黄色一级片| 国产av麻豆久久久久久久| 变态另类丝袜制服| 亚洲精华国产精华液的使用体验 | 97超级碰碰碰精品色视频在线观看| 嫩草影视91久久| 51国产日韩欧美| 国产乱人视频| 成年免费大片在线观看| 精品久久久久久久久亚洲| 91精品国产九色| 91av网一区二区| 亚洲av电影不卡..在线观看| 99热网站在线观看| 久久人妻av系列| 在线免费观看的www视频| 久久久a久久爽久久v久久| 黑人高潮一二区| 成人综合一区亚洲| 露出奶头的视频| 91在线观看av| 男女之事视频高清在线观看| h日本视频在线播放| 欧美成人一区二区免费高清观看| 国产极品精品免费视频能看的| 1024手机看黄色片| 欧美又色又爽又黄视频| 极品教师在线视频| 草草在线视频免费看| 亚洲天堂国产精品一区在线| 国产国拍精品亚洲av在线观看| 天堂av国产一区二区熟女人妻| a级毛片免费高清观看在线播放| 亚洲精华国产精华液的使用体验 | 久久久国产成人精品二区| 国产午夜精品久久久久久一区二区三区 | 麻豆一二三区av精品| 我要看日韩黄色一级片| 黄色配什么色好看| videossex国产| av.在线天堂| 在线免费十八禁| 免费av不卡在线播放| 国产高清视频在线观看网站| 波多野结衣高清作品| 国产一区二区三区在线臀色熟女| 淫妇啪啪啪对白视频| 久久久午夜欧美精品| av女优亚洲男人天堂| 欧美日韩一区二区视频在线观看视频在线 | 成人无遮挡网站| 黄色欧美视频在线观看| 亚洲av免费在线观看| 直男gayav资源| 精品一区二区免费观看| 欧美另类亚洲清纯唯美| 亚洲精品成人久久久久久| 色视频www国产| 99久久久亚洲精品蜜臀av| 精品欧美国产一区二区三| 日本黄色视频三级网站网址| 男插女下体视频免费在线播放| 国内精品一区二区在线观看| 岛国在线免费视频观看| 久久久久久久久久成人| 禁无遮挡网站| 成人性生交大片免费视频hd| 日本成人三级电影网站| 狠狠狠狠99中文字幕| 亚洲中文日韩欧美视频| 精品久久久久久久久av| 亚洲图色成人| 美女大奶头视频| 色视频www国产| 男人舔奶头视频| 人妻少妇偷人精品九色| 小说图片视频综合网站| 日韩av在线大香蕉| 国产成人freesex在线 | 国内精品久久久久精免费| 久久亚洲精品不卡| 久久草成人影院| 亚洲精品日韩在线中文字幕 | 18禁裸乳无遮挡免费网站照片| 非洲黑人性xxxx精品又粗又长| 欧美不卡视频在线免费观看| 日韩强制内射视频| 最新在线观看一区二区三区| 老熟妇仑乱视频hdxx| 久久九九热精品免费| 精品不卡国产一区二区三区| 一个人观看的视频www高清免费观看| 伊人久久精品亚洲午夜| 舔av片在线| 中文字幕av在线有码专区| 91久久精品国产一区二区三区| 可以在线观看毛片的网站| 丰满人妻一区二区三区视频av| 在线播放国产精品三级| 搡老岳熟女国产| 国产69精品久久久久777片| 亚洲真实伦在线观看| 国产成人aa在线观看| 悠悠久久av| 1024手机看黄色片| 精品一区二区三区人妻视频| 男插女下体视频免费在线播放| 国产老妇女一区| 中国美女看黄片| 亚洲色图av天堂| 国产成人a∨麻豆精品| 亚洲专区国产一区二区| 毛片一级片免费看久久久久| 婷婷精品国产亚洲av| 日本欧美国产在线视频| 无遮挡黄片免费观看| 久久久久精品国产欧美久久久| 女生性感内裤真人,穿戴方法视频| 免费看光身美女| 国产亚洲av嫩草精品影院| 日本五十路高清| 国产久久久一区二区三区| 少妇的逼好多水| 不卡视频在线观看欧美| 美女cb高潮喷水在线观看| 亚洲真实伦在线观看| 亚洲成人久久爱视频| 少妇人妻一区二区三区视频| 亚洲在线自拍视频| 欧美性感艳星| 国产精品久久久久久久电影| 桃色一区二区三区在线观看| 国产精品,欧美在线| 97超碰精品成人国产| 六月丁香七月| 精品人妻偷拍中文字幕| 九九爱精品视频在线观看| 变态另类丝袜制服| 少妇高潮的动态图| 国产美女午夜福利| 草草在线视频免费看| 麻豆一二三区av精品| 少妇猛男粗大的猛烈进出视频 | 日韩中字成人| 日韩大尺度精品在线看网址| 91在线精品国自产拍蜜月| 亚洲欧美精品自产自拍| 亚洲精品456在线播放app| 日本精品一区二区三区蜜桃| 不卡视频在线观看欧美| 国产片特级美女逼逼视频| 亚洲精品成人久久久久久| 国产色婷婷99| 99国产精品一区二区蜜桃av| 22中文网久久字幕| 欧美高清性xxxxhd video| 97在线视频观看| 亚洲av五月六月丁香网| 一进一出抽搐动态| 日韩制服骚丝袜av| 日韩成人伦理影院| 国产精华一区二区三区| 久久人人精品亚洲av| 国产在线精品亚洲第一网站| 精品免费久久久久久久清纯| 国产高清不卡午夜福利| 麻豆精品久久久久久蜜桃| 国产精品免费一区二区三区在线| 国产伦在线观看视频一区| 国内精品宾馆在线| 黄片wwwwww| 国产伦精品一区二区三区四那| 99国产精品一区二区蜜桃av| 国产乱人视频| 亚洲国产精品国产精品| 欧美性猛交╳xxx乱大交人| 夜夜看夜夜爽夜夜摸| 国产单亲对白刺激| 亚洲av成人av| 免费无遮挡裸体视频| 在线免费观看不下载黄p国产| 亚洲综合色惰| 亚洲,欧美,日韩| 1024手机看黄色片| 日韩精品有码人妻一区| 国产探花在线观看一区二区| 日韩强制内射视频| 不卡视频在线观看欧美| 亚洲婷婷狠狠爱综合网| 久久久欧美国产精品| 婷婷精品国产亚洲av在线| 一级av片app| 国产精品一及| 国产在线男女| 日日摸夜夜添夜夜爱| 午夜福利在线观看免费完整高清在 | 国产精品一及| 国产欧美日韩精品一区二区| 久久久精品欧美日韩精品| 亚洲,欧美,日韩| 欧美成人a在线观看| 91久久精品电影网| 色视频www国产| 一区福利在线观看| 亚洲在线观看片| 黄色视频,在线免费观看| 日日摸夜夜添夜夜添av毛片| 国产精品久久视频播放| 中文字幕av成人在线电影| 国产单亲对白刺激| 亚州av有码| 国产伦一二天堂av在线观看| 在线观看午夜福利视频| 天天躁夜夜躁狠狠久久av| 中文字幕免费在线视频6| 亚洲熟妇中文字幕五十中出| av在线亚洲专区| 18禁裸乳无遮挡免费网站照片| 91午夜精品亚洲一区二区三区| 精品午夜福利视频在线观看一区| 亚洲av免费在线观看| 国产午夜精品论理片| 插阴视频在线观看视频| 人妻夜夜爽99麻豆av| 如何舔出高潮| 国产成年人精品一区二区| 看片在线看免费视频| 日韩av在线大香蕉| 成人毛片a级毛片在线播放| 国产精品伦人一区二区| 国产精品国产三级国产av玫瑰| 97超碰精品成人国产| 嫩草影院入口|