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

    基于近紅外圖像的巖土濕度快速檢測(cè)系統(tǒng)

    2021-12-14 02:47:20錢昱沛陳向東陳一健夏文鶴
    物聯(lián)網(wǎng)技術(shù) 2021年12期
    關(guān)鍵詞:像素點(diǎn)高斯巖土

    錢昱沛,陳向東,丁 星,李 皋,陳一健,夏文鶴

    (1.西南交通大學(xué) 信息科學(xué)與技術(shù)學(xué)院,四川 成都 611756;2.西南石油大學(xué) 油氣藏地質(zhì)及開發(fā)工程國家重點(diǎn)實(shí)驗(yàn)室,四川 成都 610500;3.西南石油大學(xué) 石油與天然氣工程學(xué)院,四川 成都 610500;4.西南石油大學(xué) 電氣信息學(xué)院,四川 成都 610500)

    0 引 言

    氣體鉆井技術(shù)是油氣資源勘探開發(fā)的一種重要技術(shù)手段,在提高機(jī)械鉆井效率等方面具有獨(dú)特的優(yōu)勢(shì);但是在氣體鉆井過程中遇到地層出水后,井下容易發(fā)生環(huán)空堵塞、鉆頭卡鉆、井壁坍塌等復(fù)雜事故[1]。因此,在鉆井過程中需要快速檢測(cè)掘進(jìn)過程中返出巖屑的濕度水平來預(yù)警地層出水狀況。目前,巖土濕度的測(cè)量方式以烘干法最為準(zhǔn)確,但是無法達(dá)到快速檢測(cè)的目的;而常用的快速檢測(cè)方法中,最常見的為探針法,但無法實(shí)現(xiàn)大面積檢測(cè)。文獻(xiàn)[2-3]研究了通過近紅外光譜來分析土壤濕度的方法,能夠達(dá)到較高的檢測(cè)精度,但并未涉及有干擾物影響的復(fù)雜環(huán)境檢測(cè)。文獻(xiàn)[4]采用基于圖像處理的方法來檢測(cè)巖土濕度,通過圖像處理算法有效解決干擾物的影響,但其檢測(cè)靈敏度不高,而且受環(huán)境可見光影響較大,無法適應(yīng)野外鉆井的應(yīng)用場(chǎng)景。

    針對(duì)上述問題,本文設(shè)計(jì)了一種基于近紅外圖像的巖土濕度快速檢測(cè)系統(tǒng)。

    1 系統(tǒng)總體方案

    巖土濕度快速檢測(cè)的總體方案為:利用樹莓派控制USB攝像頭采集850 nm波段近紅外巖土圖像數(shù)據(jù),并通過WiFi實(shí)時(shí)發(fā)送至上位機(jī)。在上位機(jī)中,首先將圖像從RGB顏色空間變換到Lab顏色空間;然后通過混合高斯模型的建立分割出巖土區(qū)域,將分割后的只包含巖土區(qū)域的圖像轉(zhuǎn)換為灰度空間,采用一種基于小波變換的算法降低散斑噪聲干擾,通過DBSCAN聚類算法找到只包含巖土類的像素點(diǎn)并計(jì)算其灰度中心值,并將其作為分段回歸模型的輸入,輸出巖土的預(yù)測(cè)濕度值。圖1為所述系統(tǒng)框架。

    圖1 系統(tǒng)框架

    2 基于樹莓派的圖像采集和傳輸

    本系統(tǒng)主要利用便攜式嵌入式設(shè)備樹莓派控制攝像頭采集圖像,并將采集的圖像發(fā)送至上位機(jī)。采用基于TCP傳輸協(xié)議的Socket技術(shù)實(shí)現(xiàn)樹莓派與上位機(jī)之間的通信[5]。樹莓派作為硬件控制端,首先初始化Socket模塊,綁定一個(gè)IP端口并監(jiān)聽,等待上位機(jī)的連接請(qǐng)求。上位機(jī)為PC機(jī),同樣初始化Socket模塊,通過connect()函數(shù)給硬件控制端發(fā)送連接請(qǐng)求來建立起與遠(yuǎn)程硬件之間的通信橋梁。如果連接成功,上位機(jī)可以控制樹莓派實(shí)時(shí)拍攝圖像并接收樹莓派連續(xù)發(fā)送的圖像數(shù)據(jù)流,實(shí)現(xiàn)系統(tǒng)的快速檢測(cè)功能。上位機(jī)和下位機(jī)通信流程如圖2所示。

    圖2 上位機(jī)和下位機(jī)通信

    3 近紅外圖像處理算法

    3.1 基于Lab顏色空間的圖像分割

    3.1.1 Lab顏色空間模型的轉(zhuǎn)換

    巖土圖像容易受背景復(fù)雜性、光照多變性等環(huán)境因素干擾,在Lab顏色空間下建立混合高斯模型可以根據(jù)巖土的亮度特征和顏色特征來區(qū)分巖土中的干擾物。Lab顏色空間是一種基于人的視覺感應(yīng)并使其均勻化的顏色系統(tǒng),其中包含L、a、b三個(gè)特征值,L值表示亮度特征,a和b表示顏色的對(duì)立維度特征[6]。CMOS近紅外攝像頭采集到的近紅外圖像為RGB顏色空間下的圖像,首先需要映射到XYZ顏色空間[7],其關(guān)系為:

    然后通過XYZ顏色空間進(jìn)一步轉(zhuǎn)換為Lab顏色空間,轉(zhuǎn)換公式為:

    其中:

    式中,m=0.008 856。

    經(jīng)過轉(zhuǎn)換后,Lab顏色空間下的近紅外巖土圖像可以表示為X={x1,x2,...,xi,...,xN},其中xi=(xL,xa,xb)表示該圖像的每個(gè)像素點(diǎn),是一個(gè)3維向量,xL、xa、xb分別表示該像素點(diǎn)的L、a、b值。

    3.1.2 混合高斯模型的建立

    混合高斯模型是根據(jù)事物的正態(tài)分布特性采用多個(gè)高斯模型的加權(quán)和來表征數(shù)據(jù)的一個(gè)混合模型[8]。在巖土的近紅外圖像中其像素點(diǎn)存在聚集性分布,因此可以利用混合高斯模型對(duì)其進(jìn)行建模,從而分割出檢測(cè)區(qū)域中較大干擾物,降低復(fù)雜環(huán)境下誤檢測(cè)的概率。近紅外圖像的概率密度為:

    3.1.3 基于EM算法的參數(shù)學(xué)習(xí)

    EM(Expectation-Maximization)算法是混合高斯模型的常用參數(shù)學(xué)習(xí)算法[9],通過對(duì)參數(shù)θ設(shè)定初值使其變?yōu)橐阎P?,根?jù)該模型進(jìn)行E步的計(jì)算得到每個(gè)像素點(diǎn)的概率密度;然后再通過整幅圖像的概率分布進(jìn)行M步計(jì)算,采用極大似然估計(jì)法求取高斯模型的參數(shù)θ,更新E步的假設(shè)參數(shù)并重復(fù)上述過程,使其達(dá)到最終的最優(yōu)化收斂。

    首先假定圖像中的高斯分布模型是已知的,給定初值θ0。通過如下公式進(jìn)行E步計(jì)算:

    其中,ωik表示圖像中像素點(diǎn)xi在第k個(gè)高斯分布下的概率,且1≤k≤K,1≤m≤K,1≤i≤N。

    通過E步計(jì)算出概率后進(jìn)行M步計(jì)算,并更新該模型的參數(shù),具體表達(dá)式如下:

    重復(fù)以上E步和M步,直至參數(shù)達(dá)到最優(yōu)的收斂。

    3.1.4 通過先驗(yàn)知識(shí)確定巖土的高斯分布

    通過混合高斯模型可以將采集到的圖像分為K個(gè)高斯分布,這個(gè)過程屬于無監(jiān)督學(xué)習(xí),不需要先驗(yàn)知識(shí)。這K個(gè)高斯分布模型中就包括了本文需要的巖土區(qū)域的高斯模型。問題的關(guān)鍵在于并不能事先知道哪個(gè)高斯模型是該巖土的高斯模型,需要通過多次實(shí)驗(yàn)來確定。

    本文經(jīng)過多次實(shí)驗(yàn)后確定近紅外圖像中所含K個(gè)高斯模型,均值為μk=(Lk,ak,bk),其中,Lk為第k個(gè)高斯分布圖像像素區(qū)域的亮度值,ak、bk為其顏色值,將其作為該巖土區(qū)域高斯分布的先驗(yàn)知識(shí)。

    圖3(a)是含干擾物原始近紅外圖像,對(duì)其分別采用基于RGB顏色空間和基于Lab顏色空間的混合高斯模型進(jìn)行圖像分割處理,得到分割后的圖像分別見圖3(b)、圖3(c)。從圖3的對(duì)比效果可以看出,基于Lab顏色空間的混合高斯模型可以有效將巖土像素點(diǎn)從圖像的背景以及干擾中分割出來;同時(shí)相較于RGB顏色空間的建模能夠較好地分割出巖土區(qū)域中的陰影部分,提高檢測(cè)精度。

    圖3 圖像分割效果對(duì)比

    3.2 巖土灰度圖像去噪

    3.2.1 圖像顏色空間轉(zhuǎn)換

    單色灰度空間直觀地反映了巖土的亮度特征,能夠較好地分析不同濕度下的巖土對(duì)于光照的吸收反射程度。經(jīng)過圖像分割后的灰度圖像可表示為f(x,y),其中f(x,y)可由2個(gè)分量來表征,分別是入射分量和反射分量[10],如公式(11)所示:

    式中:i(x,y)的性質(zhì)取決于近紅外光源的光照強(qiáng)度;r(x,y)的性質(zhì)取決于待測(cè)巖土混合物成分,主要影響變量為濕度。因濕度不同的巖土對(duì)近紅外光的吸收不同,因此反射分量的變化會(huì)引起數(shù)字圖像函數(shù)值的變化。

    在數(shù)字圖像處理過程中,對(duì)攝像頭傳感器接收的電壓波形進(jìn)行了取樣和量化[11],單色圖像在任意坐標(biāo)(x0,y0)處的強(qiáng)度(灰度)表示為:

    l的取值范圍為lmin≤l≤lmax,其中區(qū)間[lmin,lmax]為灰度級(jí),范圍為[0,255]。

    3.2.2 基于小波變換的圖像去噪

    近紅外圖像中通常存在散斑噪聲,其中既包含加性噪聲的成分,也包含乘性噪聲的成分,表達(dá)式為:

    其中:f(x,y)為原始圖像信號(hào);nm為乘性噪聲信號(hào);na為加性噪聲信號(hào);(x,y)為含噪信號(hào)。因此,本文采用一種基于小波變換的近紅外圖像去噪方法來解決噪聲[12]。

    信號(hào)中的加性噪聲采用閾值去噪的方法濾除,濾除之后信號(hào)表達(dá)式(13)變?yōu)椋?/p>

    通過對(duì)公式兩邊取對(duì)數(shù)可將乘性噪聲變?yōu)榧有栽肼?,公式如下?/p>

    對(duì)log4(x,y)進(jìn)行離散小波變換,然后通過硬閾值規(guī)則,消除加性噪聲log4nm。

    首先選擇近似于巖土頻率信號(hào)的小波基函數(shù),對(duì)該信號(hào)進(jìn)行等間隔抽樣,經(jīng)過抽樣后得到信號(hào)樣本序列為CJ+1[k],然后對(duì)該序列進(jìn)行N級(jí)離散小波變換。目的是為了得到圖像中每個(gè)像素點(diǎn)的頻率特征,并將該像素點(diǎn)的噪聲信號(hào)分離出來[13]。通過以上過程得到了N級(jí)不同尺度的離散小波展開系數(shù)dJ[k],dJ-1[k],...,dJ-N+1[k]以及一級(jí)近似展開系數(shù)CJ-N+1[k]。

    對(duì)于小波變換的各級(jí)展開系數(shù),首先給定每級(jí)系數(shù)的閾值界限,將低于該閾值界限的系數(shù)設(shè)置為零,高于該閾值的系數(shù)保持不變,公式為:

    3.3 巖土區(qū)域密度聚類

    巖土是復(fù)雜的混合物,雖然混合高斯模型算法能夠從近紅外圖像中分割出較大干擾物,但復(fù)雜環(huán)境下的巖土包含的一些較小的干擾物也會(huì)影響檢測(cè)的準(zhǔn)確性,如石塊以及溝壑陰影。除此之外,圖像檢測(cè)屬于大面積的區(qū)域性檢測(cè),而巖土區(qū)域中還存在著水分分布不均的情況。

    針對(duì)以上問題,本文采用了密度聚類算法中的DBSCAN算法來實(shí)現(xiàn)聚類。該算法是一種基于密度的聚類算法,通過將高密度區(qū)域內(nèi)具有相同特征的點(diǎn)劃分為同一類實(shí)現(xiàn)聚類[14]??梢愿鶕?jù)區(qū)域的不同性質(zhì)和形狀聚類出具有同一特征的類別,并且可以很好地識(shí)別出噪聲點(diǎn)、排除干擾。

    DBSCAN算法需要2個(gè)參數(shù):以每個(gè)核心點(diǎn)為中心構(gòu)成鄰域的半徑ε、構(gòu)成一個(gè)聚類簇的區(qū)域內(nèi)至少包含的關(guān)聯(lián)點(diǎn)個(gè)數(shù)minPts[15]。通過設(shè)置以上參數(shù)進(jìn)行多次迭代直至最終收斂,該算法能夠找到巖土類別在灰度圖像中的簇。本文通過大量實(shí)驗(yàn)發(fā)現(xiàn),當(dāng)取ε=4、minPts=(w×h)/5時(shí)可以得到最優(yōu)化的聚類效果,其中w為近紅外巖土圖像區(qū)域的橫向像素點(diǎn)個(gè)數(shù),h為圖像區(qū)域的縱向像素點(diǎn)個(gè)數(shù)。使用DBSCAN算法能夠分割出類似于巖土的巖石塊以及巖土圖像中溝壑造成的陰影,找出巖土的像素點(diǎn)。前后的對(duì)比如圖4所示,噪聲點(diǎn)統(tǒng)一用白色點(diǎn)代替。聚類后的巖土區(qū)域可能存在水分分布不均的情況,因此需要計(jì)算該區(qū)域的灰度中心點(diǎn)并將其作為該區(qū)域濕度代表。

    圖4 巖土聚類前后效果對(duì)比

    4 實(shí)驗(yàn)部分

    4.1 模擬實(shí)驗(yàn)環(huán)境搭建

    本系統(tǒng)主要采集了不同產(chǎn)地下不同類型的三種巖土做對(duì)比實(shí)驗(yàn),樣本分別為碳酸鹽巖屑、天然農(nóng)田土和綠化帶黏土。對(duì)3種巖土分別隨機(jī)采集30組不同濕度的樣本,將每組樣本隨機(jī)抽樣分為訓(xùn)練集和測(cè)試集。

    測(cè)量巖土的實(shí)際濕度值采用的是傳統(tǒng)的烘干稱重法,可以準(zhǔn)確地測(cè)得巖土樣本的實(shí)際濕度值。其中所用到的烘干稱重儀器有:質(zhì)量精度為0.1 g的天平、盛土所用的不銹鋼容器、50 mL的量杯、ZK-82B型電熱真空干燥箱。

    采集近紅外圖像數(shù)據(jù)用的是800萬像素的索尼IMX317系列的415 USB攝像頭模組,鏡頭上貼了850 nm的窄帶濾光片防止攝像頭采集到其他波段的信息,近紅外主動(dòng)光源選用850 nm的40 W廣角近紅外補(bǔ)光燈。硬件控制平臺(tái)采用Raspberry Pi 3 Model B。

    采集數(shù)據(jù)參數(shù)設(shè)置包括攝像頭模組的鏡頭和近紅外補(bǔ)光燈與巖土樣品的距離,為48 cm。通過USB攝像頭的SDK控制攝像頭的白平衡、曝光等參數(shù)保持不變。

    4.2 實(shí)驗(yàn)過程

    取完全干燥條件下質(zhì)量為M1的巖土,將其放在質(zhì)量為M0的不銹鋼容器里。用量杯盛取一定質(zhì)量的清水加入該巖土中使其過飽和,然后攪拌均勻使巖土充分吸收水分,再將其放入真空干燥箱中進(jìn)行烘干。平均每烘干5 min后對(duì)盛有巖土的容器稱得一次質(zhì)量并放在實(shí)驗(yàn)箱內(nèi)的平臺(tái)上拍攝一次該濕度下的圖像數(shù)據(jù)。每次烘干后稱得的質(zhì)量為M2,重復(fù)該實(shí)驗(yàn)過程直至巖土完全干燥,對(duì)其他2種巖土重復(fù)上述實(shí)驗(yàn)過程。

    巖土的濕度公式為:

    在模擬實(shí)驗(yàn)環(huán)境下通過上述方法采集到的數(shù)據(jù)見表1所列。

    表1 巖土濕度與巖土圖像的灰度中心值

    4.3 實(shí)驗(yàn)結(jié)果與討論

    4.3.1 灰度特征曲線

    不同類別巖土的圖像灰度中心值隨其對(duì)應(yīng)的濕度變化曲線如圖5所示。通過圖5可以看出,圖像灰度中心值隨濕度的變化趨勢(shì)近似線性;但是不同類別巖土在每個(gè)濕度下對(duì)應(yīng)的灰度中心值不同,且灰度值隨濕度變化的斜率不同,反映了不同類別巖土的近紅外特性不同。產(chǎn)生這種差異的影響因素比較多,如化學(xué)成分的組成和含量不同、對(duì)水分的吸收能力不同等。

    圖5 巖土像素點(diǎn)的灰度中心值隨其濕度的變化曲線

    根據(jù)巖土濕度與灰度中心值的對(duì)應(yīng)關(guān)系,本文建立了巖土濕度與拍攝圖像中巖土的灰度中心值之間的分段線性回歸模型。除此之外,為了對(duì)低濕度進(jìn)行及時(shí)且準(zhǔn)確的檢測(cè),將巖土濕度分為2個(gè)類別來提升低濕度下的預(yù)測(cè)準(zhǔn)確率,其模型為:

    式中,t0代表劃分高濕度與低濕度臨界點(diǎn)的灰度值。在t處于不同的灰度值時(shí),Y的期望值分別為:

    該模型的損失函數(shù)定義為:

    模型的參數(shù)估計(jì)值為:

    式中,β={β0,β1,β2}。

    4.3.2 分段線性回歸

    實(shí)驗(yàn)中將采集的不同類別巖土在不同濕度下的30組數(shù)據(jù),以13%的濕度為分界線分為高濕度和低濕度2類,其中每類分為4個(gè)子集進(jìn)行4折交叉驗(yàn)證,得到如圖6所示的對(duì)比圖。

    圖6 不同類別巖土在交叉驗(yàn)證下的濕度實(shí)測(cè)值與預(yù)測(cè)值對(duì)比

    本文使用濕度的平均絕對(duì)誤差來衡量系統(tǒng)的準(zhǔn)確性,平均絕對(duì)誤差MAD表示為:

    式中:N為訓(xùn)練樣本的個(gè)數(shù);Tp(i)表示樣本的預(yù)測(cè)濕度值;Tf(i)表示樣本的實(shí)測(cè)濕度值。

    由式(25)計(jì)算出巖土濕度模型的誤差見表2所列。

    表2 巖土濕度模型的誤差

    通過對(duì)比實(shí)驗(yàn),得到了表2所列的不同類型巖土的識(shí)別結(jié)果,可以看出3種巖土在低濕度下的預(yù)測(cè)結(jié)果比高濕度下的預(yù)測(cè)結(jié)果更為準(zhǔn)確。對(duì)比碳酸鹽巖屑、天然農(nóng)田土和綠化帶黏土的實(shí)驗(yàn)結(jié)果,本系統(tǒng)對(duì)于不同產(chǎn)地不同類型的巖土同樣具有很高的檢測(cè)準(zhǔn)確率。在實(shí)際應(yīng)用過程中應(yīng)該首先針對(duì)不同的巖土進(jìn)行分類,然后通過建立對(duì)應(yīng)巖土類別的數(shù)據(jù)庫實(shí)現(xiàn)不同巖土濕度的預(yù)測(cè)。

    4.3.3 系統(tǒng)測(cè)試

    本系統(tǒng)進(jìn)行快速檢測(cè)測(cè)試時(shí)采用Intel i5-1135G7 CPU、Win10操作系統(tǒng)的PC機(jī)作為上位機(jī),其運(yùn)行效果如圖7所示。

    圖7 系統(tǒng)測(cè)試上位機(jī)運(yùn)行效果

    通過多次的實(shí)驗(yàn)驗(yàn)證,檢測(cè)得到一幅規(guī)格為800×600的圖像時(shí),各功能運(yùn)行大致時(shí)間見表3所列。由表3可知,系統(tǒng)需要花費(fèi)的總體時(shí)間約為50 s,基本滿足巖土濕度快速檢測(cè)的需求。

    表3 系統(tǒng)各部分運(yùn)行時(shí)間

    5 結(jié) 語

    本文針對(duì)巖土濕度快速檢測(cè)提出了基于近紅外圖像的處理方法,在不同光照環(huán)境和復(fù)雜背景下具有很強(qiáng)的魯棒性,針對(duì)不同類別的巖土還有很強(qiáng)的適應(yīng)性;同時(shí)搭建了模擬環(huán)境驗(yàn)證了該方法的可行性,針對(duì)巖土的濕度實(shí)現(xiàn)了高精度測(cè)量,對(duì)于氣體鉆井地層出水的快速檢測(cè)具有一定的參考價(jià)值。

    猜你喜歡
    像素點(diǎn)高斯巖土
    小高斯的大發(fā)現(xiàn)
    天才數(shù)學(xué)家——高斯
    基于canvas的前端數(shù)據(jù)加密
    基于逐像素點(diǎn)深度卷積網(wǎng)絡(luò)分割模型的上皮和間質(zhì)組織分割
    有限域上高斯正規(guī)基的一個(gè)注記
    《巖土力學(xué)》2014年第9 期被EI 收錄論文(40 篇,收錄率100 %)
    《巖土力學(xué)》2014年第7 期被EI 收錄論文(40 篇,收錄率100 %)
    《巖土力學(xué)》2014年第5 期被EI 收錄論文(41 篇,收錄率100%)
    《巖土力學(xué)》2014年第6期被EI收錄論文(43篇,收錄率100%)
    基于Node-Cell結(jié)構(gòu)的HEVC幀內(nèi)編碼
    国产毛片在线视频| 国产成人91sexporn| 欧美亚洲 丝袜 人妻 在线| 精品午夜福利在线看| 欧美精品国产亚洲| 亚洲av福利一区| 精品少妇内射三级| 国产精品一区www在线观看| 天堂俺去俺来也www色官网| 在线观看www视频免费| 中文字幕制服av| 亚洲第一av免费看| 久久国产精品大桥未久av | 嫩草影院入口| 黑人巨大精品欧美一区二区蜜桃 | 美女中出高潮动态图| 成人亚洲精品一区在线观看| 91aial.com中文字幕在线观看| 色哟哟·www| 亚洲精品国产av成人精品| 亚洲欧美中文字幕日韩二区| 三上悠亚av全集在线观看 | 精品国产国语对白av| 嘟嘟电影网在线观看| 亚洲人成网站在线观看播放| 日韩av不卡免费在线播放| 日日撸夜夜添| 能在线免费看毛片的网站| 久久久久久久大尺度免费视频| 人人妻人人看人人澡| 夜夜看夜夜爽夜夜摸| 少妇人妻 视频| 国产精品女同一区二区软件| 性色av一级| 国产av一区二区精品久久| 久久人妻熟女aⅴ| 少妇人妻久久综合中文| av国产久精品久网站免费入址| 国产精品伦人一区二区| 久久久久国产网址| 国产欧美日韩综合在线一区二区 | 国产在线男女| 成人毛片a级毛片在线播放| 中国三级夫妇交换| 亚洲av免费高清在线观看| 乱码一卡2卡4卡精品| 国产成人精品福利久久| 亚洲av成人精品一二三区| 日本av手机在线免费观看| 国产视频内射| 国产又色又爽无遮挡免| 2022亚洲国产成人精品| 美女xxoo啪啪120秒动态图| 狂野欧美白嫩少妇大欣赏| 国产色婷婷99| 国产免费一级a男人的天堂| 国产精品国产三级专区第一集| 亚洲美女视频黄频| 黄色怎么调成土黄色| 91精品国产国语对白视频| 亚洲精品中文字幕在线视频 | 欧美国产精品一级二级三级 | 一区在线观看完整版| 中国国产av一级| 亚洲人与动物交配视频| 18+在线观看网站| 日韩欧美 国产精品| 97超视频在线观看视频| 国内精品宾馆在线| 久久这里有精品视频免费| 三级国产精品欧美在线观看| 免费av中文字幕在线| 亚洲国产精品成人久久小说| 久久久欧美国产精品| 欧美精品人与动牲交sv欧美| 国产亚洲91精品色在线| 男人爽女人下面视频在线观看| 国产成人精品一,二区| 国产有黄有色有爽视频| 亚洲国产欧美在线一区| 国产黄色视频一区二区在线观看| 国产在线男女| 少妇人妻一区二区三区视频| 一二三四中文在线观看免费高清| 欧美成人精品欧美一级黄| 内射极品少妇av片p| 国产精品无大码| 日韩成人av中文字幕在线观看| 啦啦啦啦在线视频资源| 色94色欧美一区二区| 亚洲国产av新网站| 99九九在线精品视频 | 少妇人妻 视频| av卡一久久| 亚洲人成网站在线播| 大香蕉久久网| 成人国产av品久久久| 春色校园在线视频观看| 极品少妇高潮喷水抽搐| 久久久久久久久久人人人人人人| 亚洲伊人久久精品综合| 美女福利国产在线| 精品一区二区三卡| 日产精品乱码卡一卡2卡三| 搡老乐熟女国产| 国产中年淑女户外野战色| 国产精品国产av在线观看| 亚洲精品视频女| 成人18禁高潮啪啪吃奶动态图 | 国产日韩欧美亚洲二区| 久久人妻熟女aⅴ| 午夜激情久久久久久久| 天堂俺去俺来也www色官网| 欧美人与善性xxx| 亚洲美女黄色视频免费看| 中文字幕av电影在线播放| 99热这里只有精品一区| 久久99一区二区三区| 97超视频在线观看视频| 亚洲欧美一区二区三区国产| 人人妻人人爽人人添夜夜欢视频 | 国产日韩欧美视频二区| 亚洲无线观看免费| 国产熟女午夜一区二区三区 | 久久久久久久久久久丰满| 精品国产一区二区三区久久久樱花| 各种免费的搞黄视频| 国内精品宾馆在线| 亚洲性久久影院| 日本午夜av视频| 国产成人a∨麻豆精品| 日日啪夜夜爽| 男女边摸边吃奶| 亚洲国产精品一区三区| 十八禁高潮呻吟视频 | 99热这里只有是精品在线观看| 亚洲性久久影院| 国产欧美日韩综合在线一区二区 | 国产探花极品一区二区| 最近最新中文字幕免费大全7| 久久国内精品自在自线图片| 色5月婷婷丁香| 观看av在线不卡| 男男h啪啪无遮挡| av在线app专区| 肉色欧美久久久久久久蜜桃| 久久久久人妻精品一区果冻| 老熟女久久久| 亚洲精品aⅴ在线观看| 成人影院久久| 特大巨黑吊av在线直播| 国产精品一区二区在线观看99| 亚洲怡红院男人天堂| 日韩欧美精品免费久久| 日韩精品免费视频一区二区三区 | 尾随美女入室| 亚洲伊人久久精品综合| 99热6这里只有精品| 亚洲在久久综合| 26uuu在线亚洲综合色| 能在线免费看毛片的网站| 精品一区二区三卡| 一区二区三区精品91| 久久久久久久精品精品| 中文在线观看免费www的网站| 永久网站在线| 免费av中文字幕在线| 亚洲av综合色区一区| 亚洲中文av在线| 99久久中文字幕三级久久日本| 五月天丁香电影| 免费看不卡的av| 少妇人妻一区二区三区视频| 亚洲va在线va天堂va国产| 在线观看免费视频网站a站| 国产黄色免费在线视频| 欧美日韩国产mv在线观看视频| 成人综合一区亚洲| 一区在线观看完整版| 久久午夜福利片| 99热这里只有是精品50| 亚洲激情五月婷婷啪啪| 国模一区二区三区四区视频| 日日摸夜夜添夜夜添av毛片| 日韩精品有码人妻一区| 人妻少妇偷人精品九色| 国产精品.久久久| 中国三级夫妇交换| 亚洲av综合色区一区| 色网站视频免费| 国产成人精品婷婷| 内地一区二区视频在线| 日韩av在线免费看完整版不卡| 国产精品女同一区二区软件| 国产色婷婷99| 只有这里有精品99| 高清欧美精品videossex| h视频一区二区三区| 自拍偷自拍亚洲精品老妇| 一二三四中文在线观看免费高清| 成人午夜精彩视频在线观看| 91aial.com中文字幕在线观看| 美女内射精品一级片tv| 久久亚洲国产成人精品v| 少妇人妻 视频| 久久久久精品性色| www.av在线官网国产| 亚洲一级一片aⅴ在线观看| 51国产日韩欧美| 十八禁网站网址无遮挡 | 国产av一区二区精品久久| 观看av在线不卡| 欧美 亚洲 国产 日韩一| 欧美变态另类bdsm刘玥| 99热网站在线观看| 国产视频内射| 在线看a的网站| 在线观看美女被高潮喷水网站| 高清在线视频一区二区三区| 国产精品偷伦视频观看了| 大片免费播放器 马上看| 一级毛片aaaaaa免费看小| 欧美日韩在线观看h| 最近手机中文字幕大全| 亚洲在久久综合| 97超视频在线观看视频| 亚洲国产成人一精品久久久| 亚洲欧美精品专区久久| 夜夜骑夜夜射夜夜干| 日韩在线高清观看一区二区三区| 爱豆传媒免费全集在线观看| 中文字幕亚洲精品专区| 亚洲av国产av综合av卡| 美女cb高潮喷水在线观看| 国产精品国产av在线观看| 人人妻人人澡人人看| 麻豆成人午夜福利视频| av福利片在线观看| 日韩成人av中文字幕在线观看| 免费黄频网站在线观看国产| 亚洲精品亚洲一区二区| 亚洲精品日本国产第一区| 久久毛片免费看一区二区三区| 99久久精品热视频| 建设人人有责人人尽责人人享有的| 久久女婷五月综合色啪小说| 免费大片18禁| 80岁老熟妇乱子伦牲交| 免费不卡的大黄色大毛片视频在线观看| 日韩一区二区三区影片| 三级国产精品欧美在线观看| 大片免费播放器 马上看| 赤兔流量卡办理| 日韩制服骚丝袜av| 国产精品伦人一区二区| 亚洲av中文av极速乱| 六月丁香七月| 国产伦理片在线播放av一区| 在线观看三级黄色| 亚洲图色成人| 一区在线观看完整版| 水蜜桃什么品种好| 国产成人精品久久久久久| 久久99热这里只频精品6学生| 欧美97在线视频| 99热国产这里只有精品6| 久久精品夜色国产| 久久人妻熟女aⅴ| 一区二区三区四区激情视频| 国产伦理片在线播放av一区| 夫妻午夜视频| 日韩视频在线欧美| 十分钟在线观看高清视频www | 亚洲激情五月婷婷啪啪| 老司机亚洲免费影院| 亚洲国产日韩一区二区| 一级二级三级毛片免费看| 国产高清三级在线| 日本vs欧美在线观看视频 | av免费在线看不卡| 狠狠精品人妻久久久久久综合| 边亲边吃奶的免费视频| 亚洲综合精品二区| 日韩熟女老妇一区二区性免费视频| 久久久精品94久久精品| 欧美97在线视频| 一区二区av电影网| 亚洲国产欧美日韩在线播放 | 国产91av在线免费观看| 亚洲真实伦在线观看| 99久久精品国产国产毛片| 成年人免费黄色播放视频 | av在线老鸭窝| 少妇人妻 视频| 国产精品一二三区在线看| 国产成人精品久久久久久| 国精品久久久久久国模美| 久久婷婷青草| av又黄又爽大尺度在线免费看| 免费看日本二区| 精品人妻偷拍中文字幕| 人妻人人澡人人爽人人| 亚洲av二区三区四区| 成人18禁高潮啪啪吃奶动态图 | 高清视频免费观看一区二区| √禁漫天堂资源中文www| 又爽又黄a免费视频| 欧美日韩视频精品一区| 亚洲精品第二区| 在线天堂最新版资源| 亚洲色图综合在线观看| 欧美+日韩+精品| 久久人人爽人人爽人人片va| 又黄又爽又刺激的免费视频.| 国产一区二区在线观看日韩| 久久精品久久久久久久性| 边亲边吃奶的免费视频| 亚洲av免费高清在线观看| 国产白丝娇喘喷水9色精品| 一级二级三级毛片免费看| 久久狼人影院| 黑人高潮一二区| 国产爽快片一区二区三区| 只有这里有精品99| 99热全是精品| 天天躁夜夜躁狠狠久久av| 春色校园在线视频观看| 99精国产麻豆久久婷婷| 99九九在线精品视频 | 老司机亚洲免费影院| 人人妻人人澡人人爽人人夜夜| 青青草视频在线视频观看| 久久婷婷青草| 国产精品免费大片| 午夜影院在线不卡| 国产精品国产三级国产av玫瑰| 欧美日韩亚洲高清精品| 深夜a级毛片| 乱系列少妇在线播放| 在线天堂最新版资源| 国产成人免费无遮挡视频| 蜜桃在线观看..| 亚洲内射少妇av| tube8黄色片| 国产极品粉嫩免费观看在线 | 九九爱精品视频在线观看| 国产精品成人在线| 日韩一区二区三区影片| 国产亚洲5aaaaa淫片| 国产黄色视频一区二区在线观看| 亚洲av欧美aⅴ国产| 亚洲va在线va天堂va国产| 女性生殖器流出的白浆| 久久 成人 亚洲| 久久精品国产亚洲av涩爱| 一级二级三级毛片免费看| av视频免费观看在线观看| 中文精品一卡2卡3卡4更新| 新久久久久国产一级毛片| 精品国产一区二区久久| 晚上一个人看的免费电影| 亚洲欧美精品专区久久| 亚洲av综合色区一区| 2022亚洲国产成人精品| 人妻系列 视频| 91午夜精品亚洲一区二区三区| 国产老妇伦熟女老妇高清| 中文字幕人妻丝袜制服| 大香蕉97超碰在线| 国产免费福利视频在线观看| 久久久久网色| 22中文网久久字幕| 亚洲国产色片| 天天操日日干夜夜撸| 2021少妇久久久久久久久久久| 高清午夜精品一区二区三区| 春色校园在线视频观看| 日本黄大片高清| 高清欧美精品videossex| 国产av精品麻豆| 国产深夜福利视频在线观看| 亚洲av国产av综合av卡| 水蜜桃什么品种好| 男人舔奶头视频| 国产美女午夜福利| 91午夜精品亚洲一区二区三区| 人妻少妇偷人精品九色| 亚洲精品乱码久久久久久按摩| 国产淫片久久久久久久久| 制服丝袜香蕉在线| 欧美日韩av久久| 纵有疾风起免费观看全集完整版| 精品一区二区三区视频在线| 久久热精品热| 水蜜桃什么品种好| 久久久午夜欧美精品| 51国产日韩欧美| 久久久亚洲精品成人影院| 国产老妇伦熟女老妇高清| 国产精品一区二区性色av| 色哟哟·www| 男人添女人高潮全过程视频| 国产亚洲欧美精品永久| 午夜影院在线不卡| 国产国拍精品亚洲av在线观看| 美女中出高潮动态图| 又粗又硬又长又爽又黄的视频| 一级爰片在线观看| 国产男女内射视频| 少妇的逼水好多| 日韩强制内射视频| 在线精品无人区一区二区三| 一二三四中文在线观看免费高清| 久久国产精品大桥未久av | 国产精品一区二区性色av| 麻豆成人午夜福利视频| av在线观看视频网站免费| 婷婷色av中文字幕| 亚洲欧美成人综合另类久久久| 精品99又大又爽又粗少妇毛片| 精品一区二区免费观看| 亚洲三级黄色毛片| 久久午夜综合久久蜜桃| 亚洲精品一二三| 中国国产av一级| 国产在线视频一区二区| 最后的刺客免费高清国语| 在线看a的网站| 狂野欧美白嫩少妇大欣赏| 国产精品福利在线免费观看| 欧美成人午夜免费资源| xxx大片免费视频| 精品视频人人做人人爽| 亚洲情色 制服丝袜| 国产在线视频一区二区| 一区二区三区乱码不卡18| 国产免费福利视频在线观看| 久久久久网色| 街头女战士在线观看网站| 建设人人有责人人尽责人人享有的| 日本午夜av视频| 亚洲国产精品国产精品| 大片电影免费在线观看免费| 欧美老熟妇乱子伦牲交| 男男h啪啪无遮挡| 欧美三级亚洲精品| 51国产日韩欧美| 深夜a级毛片| 乱人伦中国视频| 自拍欧美九色日韩亚洲蝌蚪91 | 国产精品久久久久久av不卡| 国产免费福利视频在线观看| 三级国产精品片| 26uuu在线亚洲综合色| 精品亚洲成国产av| 激情五月婷婷亚洲| 亚洲国产日韩一区二区| 免费观看的影片在线观看| 欧美精品一区二区大全| 建设人人有责人人尽责人人享有的| 少妇被粗大猛烈的视频| 亚洲av成人精品一区久久| 多毛熟女@视频| 精品少妇黑人巨大在线播放| 精品久久久精品久久久| 不卡视频在线观看欧美| 六月丁香七月| 久久久久久久久久久久大奶| 如何舔出高潮| 十分钟在线观看高清视频www | 极品人妻少妇av视频| 一级爰片在线观看| 男女国产视频网站| av天堂久久9| 国产爽快片一区二区三区| 99久久精品热视频| 乱人伦中国视频| 一边亲一边摸免费视频| 两个人免费观看高清视频 | 亚洲丝袜综合中文字幕| 午夜激情久久久久久久| 99国产精品免费福利视频| 青春草亚洲视频在线观看| 成年人免费黄色播放视频 | 久久精品熟女亚洲av麻豆精品| 国产欧美另类精品又又久久亚洲欧美| 免费黄网站久久成人精品| www.色视频.com| 女的被弄到高潮叫床怎么办| 大码成人一级视频| 成人综合一区亚洲| 黄色怎么调成土黄色| a级片在线免费高清观看视频| 91在线精品国自产拍蜜月| 国产精品99久久久久久久久| 亚洲欧美清纯卡通| 亚洲电影在线观看av| 99热6这里只有精品| 2022亚洲国产成人精品| 久久久久久久久久久丰满| 日日撸夜夜添| 91精品一卡2卡3卡4卡| av天堂久久9| 成人亚洲精品一区在线观看| 国模一区二区三区四区视频| 亚洲精品视频女| 在线免费观看不下载黄p国产| 观看av在线不卡| 狠狠精品人妻久久久久久综合| 在线观看三级黄色| 日韩一区二区视频免费看| 天美传媒精品一区二区| 国产一区二区三区综合在线观看 | 亚洲综合色惰| a级一级毛片免费在线观看| 免费久久久久久久精品成人欧美视频 | av在线观看视频网站免费| 天堂俺去俺来也www色官网| 美女国产视频在线观看| 免费黄频网站在线观看国产| 国产日韩一区二区三区精品不卡 | 亚洲性久久影院| 九色成人免费人妻av| 人妻制服诱惑在线中文字幕| 中国国产av一级| 99精国产麻豆久久婷婷| 免费观看的影片在线观看| 只有这里有精品99| 成年人午夜在线观看视频| 国产av精品麻豆| 国产精品熟女久久久久浪| 亚洲欧美精品专区久久| 精品99又大又爽又粗少妇毛片| 九色成人免费人妻av| 国产亚洲av片在线观看秒播厂| 久久影院123| 国产精品国产av在线观看| 久久人人爽av亚洲精品天堂| 精品国产露脸久久av麻豆| 亚洲欧美清纯卡通| 欧美精品一区二区免费开放| 18禁动态无遮挡网站| 国产亚洲欧美精品永久| 蜜桃在线观看..| 一级毛片我不卡| 欧美xxⅹ黑人| 亚洲av不卡在线观看| 日韩制服骚丝袜av| 久久亚洲国产成人精品v| 色视频在线一区二区三区| 寂寞人妻少妇视频99o| videossex国产| 永久免费av网站大全| 精品99又大又爽又粗少妇毛片| 亚洲国产av新网站| 熟女av电影| 国产av国产精品国产| 久久久久久久久久久久大奶| 丰满饥渴人妻一区二区三| 十分钟在线观看高清视频www | 国产一区二区在线观看日韩| 亚洲美女搞黄在线观看| 搡老乐熟女国产| 内射极品少妇av片p| 亚洲精品日本国产第一区| 国产永久视频网站| 亚洲欧美中文字幕日韩二区| 自拍偷自拍亚洲精品老妇| 在线观看一区二区三区激情| 九色成人免费人妻av| 人人妻人人澡人人爽人人夜夜| 免费在线观看成人毛片| 日本午夜av视频| 久热这里只有精品99| 校园人妻丝袜中文字幕| 伊人亚洲综合成人网| 国产精品一二三区在线看| 亚洲自偷自拍三级| 日韩一区二区视频免费看| videossex国产| 欧美区成人在线视频| 国产精品一区www在线观看| 精品一区二区免费观看| 午夜影院在线不卡| 久久国产亚洲av麻豆专区| 成人毛片a级毛片在线播放| 99热这里只有是精品在线观看| 精品午夜福利在线看| 性高湖久久久久久久久免费观看| 99久久精品国产国产毛片| 少妇被粗大猛烈的视频| 国产又色又爽无遮挡免| 亚洲精品日本国产第一区| 国产欧美日韩综合在线一区二区 | 亚洲av.av天堂| 国产色婷婷99| 国产免费一区二区三区四区乱码| av.在线天堂| 在线亚洲精品国产二区图片欧美 | 青青草视频在线视频观看| 久久久久久久久久久免费av| 国产在视频线精品| 免费人妻精品一区二区三区视频| 欧美一级a爱片免费观看看| 精品一区在线观看国产| 一级毛片黄色毛片免费观看视频| 在线免费观看不下载黄p国产| 精品国产一区二区久久| 乱人伦中国视频| 日韩成人av中文字幕在线观看| 国内少妇人妻偷人精品xxx网站| 香蕉精品网在线| 黑人猛操日本美女一级片| 国产精品一区二区三区四区免费观看| 在线免费观看不下载黄p国产|