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

    基于特征參數(shù)的井眼軌跡模式識別與分段方法

    2022-10-16 02:29:10劉衛(wèi)星楊金顯
    電子科技 2022年10期
    關(guān)鍵詞:樣條特征參數(shù)變化率

    劉衛(wèi)星,楊金顯

    (河南理工大學(xué) 電氣工程與自動(dòng)化學(xué)院,河南 焦作 454000)

    隨著油氣勘探的發(fā)展,定向鉆井軌跡擬合精度受到越來越多的重視。定向鉆井的核心技術(shù)是監(jiān)測和控制,監(jiān)測的任務(wù)是根據(jù)測斜數(shù)據(jù)計(jì)算出實(shí)鉆軌跡;控制的任務(wù)是根據(jù)監(jiān)測結(jié)果控制鉆進(jìn)方向[1-3]。由于測斜時(shí)只能獲得一系列離散測點(diǎn)處的井深、井斜角、方位角等基本參數(shù),無法還原各測段內(nèi)井眼軌跡的真實(shí)形態(tài),因此需要基于假設(shè)條件和數(shù)學(xué)模型來計(jì)算井眼軌跡的空間坐標(biāo)。目前用于提高井眼軌跡擬合精度的方法主要為兩種:(1)硬件條件改善。通過增加測斜點(diǎn)的密度,提高軌跡的擬合度和精確度。但這種方法增加了鉆井測斜的成本,難以被廣泛應(yīng)用在實(shí)際工程中;(2)計(jì)算方法改善,也就是使用更加切合實(shí)際的井眼軌跡計(jì)算方法。雖然這種方法對計(jì)算能力有很高的要求,但是隨著計(jì)算機(jī)的迅速發(fā)展,其對計(jì)算能力的需求也將逐漸被滿足。

    目前國內(nèi)外針對提高井眼軌跡測量精度的計(jì)算方法取得了一定的研究進(jìn)展。文獻(xiàn)[4~8]假設(shè)井眼軌跡為空間斜平面內(nèi)的圓弧,分別提出了用于測斜計(jì)算的最小曲率法。這些方法經(jīng)過不斷地改進(jìn)和完善,進(jìn)一步發(fā)展為空間圓弧模型。文獻(xiàn)[9]假設(shè)井眼軌跡在垂直剖面圖和水平投影圖上均為圓弧,提出了曲率半徑模型。由于這樣的井眼軌跡被證明為圓柱螺線,因此常被各行業(yè)用作擬合軌跡的模型[10]。文獻(xiàn)[11]假設(shè)井眼軌跡的井斜變化率和方位變化率分別為常數(shù),提出了自然曲線模型。雖然研究人員陸續(xù)提出了多種井眼軌跡曲線的計(jì)算方法,但是上述軌跡計(jì)算方法都是假設(shè)測段為單一的井眼軌跡模型。需要指出的是,在實(shí)際鉆進(jìn)過程中,井眼軌跡不可能總按照某種單一的軌跡模型進(jìn)行鉆進(jìn),因此單一的軌跡模型會使井眼軌跡計(jì)算產(chǎn)生較大的誤差。為了提高軌跡擬合精度,文獻(xiàn)[6]對不同的井眼軌跡模型使用數(shù)值積分的方法進(jìn)行分析,并推導(dǎo)出了對應(yīng)模型的數(shù)值積分計(jì)算式。這種方法提高了井眼軌跡的計(jì)算精度,但是計(jì)算過程仍基于某一假設(shè)模型。隨后,又有研究人員提出將測段平均分為多條直線段,再對各直線段進(jìn)行特殊積分的分段式井眼軌跡計(jì)算方法[12-13]。該方法雖然對軌跡中間進(jìn)行分段,但是分段方式無據(jù)可依,存在較大的主觀性,使得軌跡的分段對實(shí)際軌跡曲線的匹配度較低。為解決這個(gè)問題,文獻(xiàn)[14]提出通過工具面角對井眼軌跡模式定量識別的方法來選擇符合實(shí)際井眼軌跡的軌跡模型,但該方法無法精確地判斷軌跡模型。

    采用上述軌跡計(jì)算方法對測段采用單一軌跡模型求取井眼軌跡,不僅精度較低,且軌跡分段的主觀性較大。針對這些問題,本文提出了一種利用多種軌跡模型的特征參數(shù)來對軌跡進(jìn)行模式識別和分段的方法。該方法首先利用三次b樣條插值法來構(gòu)建井斜角、方位角隨井深變化的函數(shù),然后分析特征參數(shù)與軌跡模型特征參數(shù)的區(qū)別,并建立模式識別與分段評價(jià)標(biāo)準(zhǔn)。

    1 井眼軌跡測點(diǎn)鄰近結(jié)構(gòu)

    井眼軌跡是一條光滑的曲線,軌跡上每個(gè)點(diǎn)的坐標(biāo)都隨井深的變化而變化,因此可以構(gòu)建井眼軌跡的向量形式

    r(L)=Ni+Ej+Hk

    (1)

    式中,L為井深;N、E、H分別代表南北坐標(biāo)、東西坐標(biāo)、垂深,單位均為m;i、j、k表示各坐標(biāo)軸的單位向量。根據(jù)式(1)可知,為了研究L上任選的一點(diǎn)r(LP)的近鄰形狀,在其附近選取另一點(diǎn)r(LP+ΔL),根據(jù)泰勒計(jì)算式可得

    (2)

    當(dāng)ΔL很小時(shí),ζ趨向于0。

    (3)

    式中,tp、np、bp分別表示在點(diǎn)r(LP)的單位切向量、單位法向量、單位副法向量;k為曲線的曲率;τ為曲線的撓率。軌跡坐標(biāo)增量與軌跡某點(diǎn)的曲率和撓率有關(guān)。

    (4)

    式中,kα為井斜變化率;kφ為方位變化率。從式(4)中可以看出,可由井斜變化率和方位變化率計(jì)算獲得曲率和撓率,因此井斜變化率與方位變化率也能夠決定軌跡曲線某點(diǎn)的鄰近形態(tài)。

    2 樣條函數(shù)基礎(chǔ)理論

    井眼軌跡是連續(xù)光滑的空間曲線,在鉆井過程中會存在彎曲與扭轉(zhuǎn)。井眼曲率、井斜角和方位角變化率是決定井眼鄰近形狀的因子,可以作為判斷井眼軌跡模型的參數(shù)。從力學(xué)角度分析,因受井壁的限制,鉆柱可簡化為一條彈性樣條。從數(shù)學(xué)角度分析,隨著井深的變化,井斜角與方位角可近似表達(dá)為隨井深變化的函數(shù)。從鉆井工程應(yīng)用角度分析,用樣條函數(shù)[15]模擬井眼軌跡具有高度的逼真。文獻(xiàn)[16]的理論分析證明了樣條函數(shù)在井眼軌跡分析中的可用性,說明了插值方法[17]是逼近這種函數(shù)的有效方法。

    b樣條函數(shù)是一種特殊的函數(shù),是形態(tài)數(shù)學(xué)描述的主流方法之一。在鉆井過程中,軌跡井斜角與方位角連續(xù)變化,且變化與高度相關(guān)。在井眼軌跡測段內(nèi),根據(jù)樣條函數(shù)定義,建立井斜角與方位角樣條函數(shù)方程組,如下所示

    (5)

    ΔLi=Li+1-Li。

    式中,ΔL為井深變化量,單位為m。

    對于整個(gè)井眼軌跡的n個(gè)測點(diǎn),可以得到n個(gè)線性方程組。利用測斜數(shù)據(jù)可以確定θ(L)與Φ(L)的樣條函數(shù),進(jìn)而可以求得每個(gè)測點(diǎn)井眼曲率κ、井斜變化率κθ、井眼撓率τ、方位變化率κΦ隨井深的變化函數(shù)

    (6)

    (7)

    (8)

    3 模式識別及分段

    3.1 空間圓弧模型

    空間圓弧模型表示的是某一空間斜面上的圓弧,其曲率k為常值。圖1為空間圓弧的軌跡模型示意圖。

    圖1 空間圓弧軌跡模型示意圖

    (9)

    式中,ki,i-1、εi,i-1分別為測點(diǎn)i與i-1之間測段的曲率和狗腿角;Li代表測點(diǎn)i的井深。由于撓率τ為偏離該平面的程度因子,所以該模型的特征參數(shù)為

    (10)

    式中,C為常數(shù)。

    3.2 圓柱螺旋模型

    圓柱螺線軌跡模型的特征是軌跡在垂直剖面圖和水平投影圖上的曲率為常數(shù)。在螺旋線上運(yùn)動(dòng)的點(diǎn),由直線運(yùn)動(dòng)和圓周運(yùn)動(dòng)合成,如圖2所示。

    圖2 圓柱螺旋跡模型示意圖

    當(dāng)空間點(diǎn)從p0出發(fā),在水平面以圓周的形式運(yùn)動(dòng),同時(shí)又以等速度沿Z軸正向運(yùn)動(dòng),其參數(shù)方程為

    (11)

    式中,r為圓周半徑;v為沿Z軸的速度。

    螺旋上的任意一點(diǎn)的位移矢量可表示為

    R(s)=x(s)i+y(s)j+z(s)k

    (12)

    式中,s為井深;i、j、k分別為X、Y、Z軸單位向量,則螺旋線表示為

    (13)

    式中,L為螺距;l為一個(gè)螺距內(nèi)的弧長;r為水平圓的半徑,曲率和撓率如式(14)所示。

    (14)

    從上式可知,圓柱螺旋法的曲率和撓率為常數(shù)。曲率和撓率分別為

    (15)

    式中,C1、C2為常數(shù)。

    3.3 自然參數(shù)法

    自然參數(shù)模型的特征是軌跡的井斜角變化率和方位角變化率都為常數(shù)。將特征參數(shù)分別記為kα、kφ,其表達(dá)式為

    (16)

    式中,αi、φi、Li分別代表第i個(gè)測點(diǎn)的井斜角、方位角、井深。因此其特征參數(shù)為

    (17)

    其中,C1和C2為常數(shù)。

    3.4 正切模型

    正切法軌跡模型是假設(shè)測段為直線軌跡。因此井斜與方位變化率都為0。

    3.5 模式識別標(biāo)準(zhǔn)

    井眼軌跡中測點(diǎn)鄰近點(diǎn)的形態(tài)可以根據(jù)測點(diǎn)特征參數(shù)進(jìn)行描述。通過樣條函數(shù)計(jì)算各測點(diǎn)的特征參數(shù),對比測斜模型來識別軌跡模式。假設(shè)該測段有n個(gè)測點(diǎn),根據(jù)下列算式計(jì)算出這n個(gè)測點(diǎn)的特征參數(shù)的平均值以及標(biāo)準(zhǔn)差。

    井斜變化率的平均值及標(biāo)準(zhǔn)差為

    (18)

    方位變化率的平均值及標(biāo)準(zhǔn)差為

    (19)

    曲率的平均值及標(biāo)準(zhǔn)差為

    (20)

    撓率的平均值及標(biāo)準(zhǔn)差如式(21)所示。

    (21)

    在所述識別區(qū)間內(nèi),考慮到井眼軌跡的不同特征參數(shù)對模式識別的影響,本文采用特征參數(shù)的標(biāo)準(zhǔn)差來作為模式識別的評價(jià)標(biāo)準(zhǔn)。標(biāo)準(zhǔn)差越小,則認(rèn)為該模型與實(shí)鉆軌跡越接近。因?yàn)榭臻g圓弧法與圓柱螺旋法的特征參數(shù)相同,故將空間圓弧模型與圓柱螺旋模型假設(shè)為一個(gè)模型,再通過撓率的平均值進(jìn)行判斷。如果其平均值接近0,則認(rèn)為是空間圓弧模型;反之,為圓柱螺旋模型。

    由于每種模式都有兩個(gè)特征參數(shù),為了方便敘述,分別用μ和ν來表示??臻g圓弧與圓柱螺旋線模型中,μ和ν為

    (22)

    自然參數(shù)模型中,μ和ν為

    (23)

    式中,i為測段標(biāo)記。

    由于軌跡模型的特征參數(shù)都為常值,所以可根據(jù)標(biāo)準(zhǔn)差進(jìn)行判斷,標(biāo)準(zhǔn)差越小,該模型與實(shí)鉆軌跡越相似。然而,每種軌跡模型都有兩個(gè)評價(jià)參數(shù),且不同參數(shù)間相互獨(dú)立,有時(shí)會出現(xiàn)一個(gè)特征參數(shù)標(biāo)準(zhǔn)差較大,而另一個(gè)參數(shù)標(biāo)準(zhǔn)差較小的情況??紤]到各特征參數(shù)的權(quán)重,重新構(gòu)造井眼軌跡模式的評價(jià)指標(biāo),即

    σ=γ×μ+(1-γ)×v

    (24)

    式中,σ為評價(jià)指標(biāo);γ為權(quán)重系數(shù)。

    空間圓弧與圓柱螺旋線的識別方法為:首先通過σ來識別空間圓弧模型與圓柱螺線模型,再根據(jù)其撓率的平均值進(jìn)行篩選,從而篩選出最符合實(shí)鉆軌跡的測斜模型。

    3.6 數(shù)據(jù)點(diǎn)的前尋分段

    上述過程只能判斷某段的測斜方法,但在一條完整的井眼軌跡中,軌跡的形態(tài)是多變的,并不知道軌跡在何處發(fā)生變化,因此需要遍歷所有的點(diǎn)來對全段井眼軌跡進(jìn)行分段,從而避免將某段軌跡測斜方法當(dāng)作全段軌跡的測斜方法。本文使用逐點(diǎn)前尋的方法來精確地分段。對于空間曲線,需最少用4個(gè)點(diǎn)來描述其形態(tài),因此本文選用前4個(gè)測點(diǎn)進(jìn)行模式識別。向前尋找下一個(gè)數(shù)據(jù)點(diǎn)帶入樣條函數(shù),得出其特征參數(shù),并判斷這一點(diǎn)的特征參數(shù)誤差是否在標(biāo)準(zhǔn)差內(nèi)。若滿足式(25),則繼續(xù)測試下個(gè)測點(diǎn);若不滿足條件,則判斷上一點(diǎn)為分界點(diǎn),并以上一點(diǎn)為初始點(diǎn)重新選擇測斜模式。該方法能夠最大限度地進(jìn)行分段,并求取某段的測斜模式。

    假設(shè)當(dāng)前測段所計(jì)算的特征接近某種模型,其特征參數(shù)為kα和kβ,通過樣條函數(shù)求出下一點(diǎn)的kα,p和kβ,p,并根據(jù)以下計(jì)算式進(jìn)行判斷。

    (25)

    如果滿足上述條件,則認(rèn)為此點(diǎn)也在該測段上,如果不滿足,則需重新分段。

    擬合算法的具體步驟如下:

    步驟1輸入測點(diǎn)的測量數(shù)據(jù),將其帶入樣條函數(shù),求出各特征參數(shù)隨井深的變化函數(shù);

    步驟2選取前4個(gè)點(diǎn)來判斷軌跡模式;

    步驟3向前讀取下一測點(diǎn)的特征參數(shù),判斷其是否滿足誤差。若滿足,則繼續(xù)讀取下一點(diǎn);若不滿足,則以上一點(diǎn)為起始點(diǎn),返回步驟2重新判斷軌跡模式;

    步驟4直到最后一點(diǎn)讀取完成,結(jié)束循環(huán)。

    4 模擬實(shí)驗(yàn)與分析

    通過隨鉆測量樣機(jī)設(shè)備進(jìn)行模擬實(shí)鉆實(shí)驗(yàn)。樣機(jī)與實(shí)際隨鉆測量設(shè)備功能相同,僅在尺寸和材質(zhì)上有所縮減。為了能夠清楚地描述實(shí)驗(yàn)結(jié)果,按照樣機(jī)與實(shí)際鉆井設(shè)備的比例,將采集數(shù)據(jù)中井深與坐標(biāo)增量進(jìn)行1∶10的比例尺換算(以換算為實(shí)鉆數(shù)據(jù)進(jìn)行分析)。為了模擬實(shí)鉆,在實(shí)驗(yàn)時(shí),樣機(jī)設(shè)備按照設(shè)計(jì)的軌跡方向鉆進(jìn),在一定的井深做出一些增斜或降斜操作。樣機(jī)設(shè)備每隔約100 s停鉆,在此時(shí)采集靜態(tài)測量數(shù)據(jù),用以模擬實(shí)際鉆井過程中的停鉆狀態(tài)。模擬實(shí)驗(yàn)期間,共進(jìn)行30次停鉆操作。選取鉆進(jìn)過程中100~400 m的測段進(jìn)行算法分析。樣機(jī)設(shè)備采樣頻率設(shè)為50 Hz,在鉆進(jìn)的過程中采集加速度及地磁數(shù)據(jù),并換算成井斜角、方位角及井深。

    首先將測段中井斜角、方位角和井深帶入式(5)以求取軌跡基本參數(shù),計(jì)算測點(diǎn)參數(shù)的平均值以及標(biāo)準(zhǔn)差,如表1所示。然后,根據(jù)各參數(shù)求出井眼軌跡的各基本參數(shù)的平均值以及標(biāo)準(zhǔn)差,如表2所示。

    表1 軌跡基本參數(shù)

    表2 參數(shù)的平均值以及標(biāo)準(zhǔn)差

    根據(jù)分段點(diǎn)選取原則,選擇前4個(gè)點(diǎn)的特征參數(shù)帶入模式識別評價(jià)方程,得到每個(gè)模式的評價(jià)標(biāo)準(zhǔn),如表3所示。

    表3 模型評價(jià)參數(shù)1

    由于評價(jià)標(biāo)準(zhǔn)越小,井眼軌跡越接近該標(biāo)準(zhǔn)下的軌跡模型,所以根據(jù)前4個(gè)測點(diǎn)可判斷該測段為自然參數(shù)模型。

    由于在270 m時(shí)計(jì)算獲得的井斜變化率以及方位變化率的值大于評價(jià)指標(biāo),因此可以判斷上個(gè)測點(diǎn)為分界點(diǎn)。如果不對井眼軌跡進(jìn)行分段,仍繼續(xù)使用自然參數(shù)模型,并將該點(diǎn)帶入計(jì)算,將發(fā)現(xiàn)其標(biāo)準(zhǔn)差變化很大,因此本文所提出的分界方法是有效的。

    對分界點(diǎn)后的軌跡進(jìn)行上述分析,計(jì)算各模型的評價(jià)指標(biāo),如表4所示。根據(jù)識別標(biāo)準(zhǔn)判斷剩下的軌跡與空間圓弧法或圓柱螺線法之間的關(guān)系。由于撓率的平均為0.45,基本接近0,所以判斷為空間圓弧法。

    表4 模型評價(jià)參數(shù)2

    表5以絕對誤差為判斷依據(jù),表明通過對軌跡進(jìn)行分段與模型識別能夠提高軌跡計(jì)算的精度。

    表5 軌跡誤差

    根據(jù)單一軌跡與本文所用方法描繪井眼軌跡圖,如圖3所示。從圖3可以看出,相比單一井眼軌跡,采用本文所提方法繪制的軌跡更接近參考軌跡。

    圖3 井眼軌跡的三維圖

    5 結(jié)束語

    測段采用單一軌跡模型求取井眼軌跡時(shí)精度較低,且軌跡分段的主觀性較大。針對該問題,本文從理論角度分析了三次b樣條的適用性,通過該方法計(jì)算出井眼軌跡中各參數(shù)隨井深的變化函數(shù),然后構(gòu)建模式識別與分段方法,通過識別方法判斷該測段最適合的模型,再利用數(shù)據(jù)點(diǎn)前尋法遍歷所有軌跡測點(diǎn)對其進(jìn)行分段。本文設(shè)計(jì)了隨鉆測量樣機(jī)設(shè)備的模擬鉆進(jìn)實(shí)驗(yàn)。實(shí)例誤差分析和軌跡圖分析證明了本文提出的基于特征參數(shù)的井眼模式識別和分段方法的有效性。

    猜你喜歡
    樣條特征參數(shù)變化率
    一元五次B樣條擬插值研究
    故障診斷中信號特征參數(shù)擇取方法
    基于電流變化率的交流濾波器失諧元件在線辨識方法
    湖南電力(2021年4期)2021-11-05 06:44:42
    基于特征參數(shù)化的木工CAD/CAM系統(tǒng)
    例談中考題中的變化率問題
    三次參數(shù)樣條在機(jī)床高速高精加工中的應(yīng)用
    基于PSO-VMD的齒輪特征參數(shù)提取方法研究
    三次樣條和二次刪除相輔助的WASD神經(jīng)網(wǎng)絡(luò)與日本人口預(yù)測
    軟件(2017年6期)2017-09-23 20:56:27
    基于樣條函數(shù)的高精度電子秤設(shè)計(jì)
    利用基波相量變化率的快速選相方法
    欧美日本亚洲视频在线播放| 两个人视频免费观看高清| 国产午夜精品久久久久久一区二区三区 | 毛片女人毛片| 亚洲av熟女| 国国产精品蜜臀av免费| 男女下面进入的视频免费午夜| 人人妻,人人澡人人爽秒播| 日本-黄色视频高清免费观看| 国产免费男女视频| 看黄色毛片网站| 尤物成人国产欧美一区二区三区| 中出人妻视频一区二区| 少妇熟女欧美另类| 最近在线观看免费完整版| 嫩草影院精品99| 综合色av麻豆| 国产精品精品国产色婷婷| 久久久精品欧美日韩精品| 欧美精品国产亚洲| 亚洲av.av天堂| 91麻豆精品激情在线观看国产| 亚洲欧美成人综合另类久久久 | 丝袜喷水一区| 久久人妻av系列| 国产色爽女视频免费观看| 中文在线观看免费www的网站| 又黄又爽又刺激的免费视频.| 蜜桃久久精品国产亚洲av| 人妻久久中文字幕网| 久久草成人影院| 国产乱人偷精品视频| 久99久视频精品免费| 久久久色成人| 午夜精品在线福利| 国产亚洲av嫩草精品影院| 亚洲av免费高清在线观看| 精品欧美国产一区二区三| 毛片女人毛片| 国产真实乱freesex| 国产麻豆成人av免费视频| 狂野欧美白嫩少妇大欣赏| 午夜激情福利司机影院| 精品无人区乱码1区二区| 久久久久性生活片| a级毛色黄片| 日本黄色片子视频| 国产91av在线免费观看| 国产久久久一区二区三区| 国产精品永久免费网站| 在线播放国产精品三级| 国产乱人偷精品视频| 亚洲久久久久久中文字幕| 免费黄网站久久成人精品| 成年av动漫网址| 中文字幕人妻熟人妻熟丝袜美| 天堂网av新在线| 国产精品乱码一区二三区的特点| 欧美最新免费一区二区三区| 97碰自拍视频| 亚洲18禁久久av| 波多野结衣高清作品| 男女边吃奶边做爰视频| or卡值多少钱| 91在线精品国自产拍蜜月| 久久午夜福利片| 给我免费播放毛片高清在线观看| 免费观看的影片在线观看| 99久国产av精品国产电影| 国产片特级美女逼逼视频| 亚洲成人久久爱视频| 美女高潮的动态| 免费看av在线观看网站| 欧美丝袜亚洲另类| 国产精品伦人一区二区| 欧美高清性xxxxhd video| 男女那种视频在线观看| 一级毛片aaaaaa免费看小| 搡老熟女国产l中国老女人| 国产精品乱码一区二三区的特点| 黄色日韩在线| 国产爱豆传媒在线观看| 色在线成人网| av国产免费在线观看| 亚洲美女视频黄频| 国产亚洲91精品色在线| 久99久视频精品免费| 欧美3d第一页| 直男gayav资源| 日韩大尺度精品在线看网址| 一级毛片aaaaaa免费看小| 成人永久免费在线观看视频| 不卡视频在线观看欧美| 国产高清不卡午夜福利| 老熟妇仑乱视频hdxx| 久久久午夜欧美精品| 色噜噜av男人的天堂激情| 亚州av有码| 我要看日韩黄色一级片| 免费人成视频x8x8入口观看| 久久久久九九精品影院| 久久人人爽人人片av| 国产v大片淫在线免费观看| 国产aⅴ精品一区二区三区波| 亚洲三级黄色毛片| 国产成人一区二区在线| 亚洲国产精品合色在线| 精品不卡国产一区二区三区| 一进一出抽搐gif免费好疼| 亚洲国产精品sss在线观看| 在线免费观看的www视频| 成年女人永久免费观看视频| 免费观看人在逋| 夜夜爽天天搞| 国内精品宾馆在线| 99热全是精品| 直男gayav资源| 国产精品野战在线观看| 搡老熟女国产l中国老女人| 国产熟女欧美一区二区| 在线观看av片永久免费下载| 久久精品久久久久久噜噜老黄 | 日韩大尺度精品在线看网址| 成人精品一区二区免费| 三级经典国产精品| 欧美一区二区精品小视频在线| 久久国产乱子免费精品| 国产精品,欧美在线| 免费观看人在逋| 亚洲国产精品成人综合色| 久久国内精品自在自线图片| 久久精品夜色国产| 亚洲国产精品久久男人天堂| 亚洲人成网站在线观看播放| 白带黄色成豆腐渣| 亚洲美女黄片视频| 久久久久久伊人网av| av免费在线看不卡| 国产av不卡久久| а√天堂www在线а√下载| 午夜精品在线福利| 一卡2卡三卡四卡精品乱码亚洲| 久久久精品大字幕| 黄色一级大片看看| 又粗又爽又猛毛片免费看| 少妇人妻一区二区三区视频| 日日撸夜夜添| 我要看日韩黄色一级片| 久久久久国内视频| 热99re8久久精品国产| 日本在线视频免费播放| 免费av观看视频| 免费av不卡在线播放| 在线观看av片永久免费下载| 色综合亚洲欧美另类图片| 国产一区二区三区在线臀色熟女| 日韩欧美精品免费久久| 99热6这里只有精品| 一本精品99久久精品77| 搡老熟女国产l中国老女人| 日本黄大片高清| 精品国产三级普通话版| 色综合站精品国产| 亚洲精品影视一区二区三区av| 久久久色成人| 国产v大片淫在线免费观看| 大又大粗又爽又黄少妇毛片口| 国产精品亚洲美女久久久| 亚洲第一区二区三区不卡| 精品国产三级普通话版| 美女xxoo啪啪120秒动态图| 给我免费播放毛片高清在线观看| 色视频www国产| 国产精品电影一区二区三区| 欧美最新免费一区二区三区| 欧美日本亚洲视频在线播放| av免费在线看不卡| 天美传媒精品一区二区| 国产亚洲av嫩草精品影院| 亚洲精品456在线播放app| 99精品在免费线老司机午夜| 国产白丝娇喘喷水9色精品| 插逼视频在线观看| 男女下面进入的视频免费午夜| a级毛片a级免费在线| 日韩欧美精品免费久久| 女人被狂操c到高潮| 国产伦精品一区二区三区四那| 三级男女做爰猛烈吃奶摸视频| 国产精品无大码| 国产欧美日韩精品亚洲av| 美女内射精品一级片tv| 国产日本99.免费观看| 一本精品99久久精品77| 亚洲无线观看免费| 国产高潮美女av| 国产精品久久久久久精品电影| 日韩成人伦理影院| 亚洲性久久影院| 嫩草影视91久久| 99热网站在线观看| 校园人妻丝袜中文字幕| 麻豆av噜噜一区二区三区| 啦啦啦观看免费观看视频高清| 伦理电影大哥的女人| 欧美日本亚洲视频在线播放| 色综合站精品国产| 国产美女午夜福利| 国产综合懂色| h日本视频在线播放| 国产三级在线视频| 男女边吃奶边做爰视频| 1024手机看黄色片| 最近最新中文字幕大全电影3| 日韩国内少妇激情av| 亚洲国产高清在线一区二区三| 亚洲乱码一区二区免费版| 九九热线精品视视频播放| 日本欧美国产在线视频| 乱码一卡2卡4卡精品| 波野结衣二区三区在线| 丰满的人妻完整版| 欧美zozozo另类| 国产69精品久久久久777片| 国产精品永久免费网站| 欧美激情久久久久久爽电影| 大型黄色视频在线免费观看| 天堂动漫精品| 亚洲精品久久国产高清桃花| 亚洲av电影不卡..在线观看| 97超碰精品成人国产| 亚洲欧美精品自产自拍| 搡女人真爽免费视频火全软件 | 国产麻豆成人av免费视频| 国产成人a区在线观看| 免费av不卡在线播放| av国产免费在线观看| 日韩成人av中文字幕在线观看 | 国产精品综合久久久久久久免费| 99热只有精品国产| 久久天躁狠狠躁夜夜2o2o| 欧美3d第一页| 亚洲精华国产精华液的使用体验 | 悠悠久久av| 亚洲熟妇中文字幕五十中出| 亚洲精品国产av成人精品 | 久久精品久久久久久噜噜老黄 | 国产精品国产高清国产av| 久久国产乱子免费精品| av在线播放精品| 久久人人爽人人片av| 国产探花在线观看一区二区| 三级国产精品欧美在线观看| 色吧在线观看| 午夜福利在线观看吧| 别揉我奶头 嗯啊视频| 天天一区二区日本电影三级| 欧美日本亚洲视频在线播放| 亚洲性夜色夜夜综合| 最近2019中文字幕mv第一页| 黄色欧美视频在线观看| 国产极品精品免费视频能看的| 高清日韩中文字幕在线| 亚洲婷婷狠狠爱综合网| 成人亚洲精品av一区二区| 国产毛片a区久久久久| 搡女人真爽免费视频火全软件 | 欧美丝袜亚洲另类| 99热精品在线国产| 永久网站在线| 亚洲自偷自拍三级| 男人的好看免费观看在线视频| 精品久久久久久成人av| 大香蕉久久网| 亚洲成人中文字幕在线播放| 日本免费a在线| 干丝袜人妻中文字幕| av在线亚洲专区| 最好的美女福利视频网| 国产精品精品国产色婷婷| 亚洲成a人片在线一区二区| 亚洲av中文字字幕乱码综合| 精品人妻偷拍中文字幕| 在线免费观看的www视频| 午夜福利视频1000在线观看| 国产高清三级在线| 日本-黄色视频高清免费观看| 亚洲av美国av| 日本一二三区视频观看| av在线老鸭窝| 寂寞人妻少妇视频99o| 免费不卡的大黄色大毛片视频在线观看 | 亚洲欧美精品自产自拍| 一进一出好大好爽视频| 最后的刺客免费高清国语| 久久久久久久亚洲中文字幕| 夜夜看夜夜爽夜夜摸| 在线免费观看的www视频| 免费观看的影片在线观看| 亚洲天堂国产精品一区在线| 色哟哟哟哟哟哟| 国产av不卡久久| 亚洲欧美中文字幕日韩二区| 99精品在免费线老司机午夜| 男人舔奶头视频| 乱人视频在线观看| 欧美成人精品欧美一级黄| 99久久九九国产精品国产免费| 久久这里只有精品中国| 日韩欧美精品免费久久| 精品久久久久久久久久久久久| 国产午夜福利久久久久久| 床上黄色一级片| 久久精品国产亚洲网站| 两性午夜刺激爽爽歪歪视频在线观看| 国产三级在线视频| 国产精品日韩av在线免费观看| 国产中年淑女户外野战色| 久久久久久久久久成人| 18禁黄网站禁片免费观看直播| 搡老妇女老女人老熟妇| 欧洲精品卡2卡3卡4卡5卡区| 国产老妇女一区| 亚洲内射少妇av| 国产亚洲欧美98| 看十八女毛片水多多多| 丰满人妻一区二区三区视频av| 99在线视频只有这里精品首页| 免费一级毛片在线播放高清视频| 午夜影院日韩av| 91久久精品国产一区二区成人| av黄色大香蕉| 夜夜看夜夜爽夜夜摸| 成人特级av手机在线观看| 日本撒尿小便嘘嘘汇集6| 久久精品国产亚洲av涩爱 | 亚洲成人中文字幕在线播放| 毛片一级片免费看久久久久| 天堂网av新在线| 婷婷精品国产亚洲av在线| 久久6这里有精品| 国产大屁股一区二区在线视频| 国产精品综合久久久久久久免费| 综合色av麻豆| 国产中年淑女户外野战色| a级一级毛片免费在线观看| 国产一区二区三区av在线 | 性色avwww在线观看| 国产欧美日韩精品一区二区| 欧美日本亚洲视频在线播放| 一区二区三区四区激情视频 | 在现免费观看毛片| 天堂影院成人在线观看| ponron亚洲| 一夜夜www| 真实男女啪啪啪动态图| 国内少妇人妻偷人精品xxx网站| 俄罗斯特黄特色一大片| 欧洲精品卡2卡3卡4卡5卡区| 麻豆久久精品国产亚洲av| 中国国产av一级| 亚洲国产精品久久男人天堂| 国产免费一级a男人的天堂| 欧美一区二区亚洲| 99热这里只有精品一区| 亚洲成人av在线免费| 国产精品一区二区三区四区久久| 啦啦啦啦在线视频资源| 婷婷亚洲欧美| 九色成人免费人妻av| 美女高潮的动态| 99视频精品全部免费 在线| 一级毛片电影观看 | av在线蜜桃| 男女之事视频高清在线观看| 老师上课跳d突然被开到最大视频| 久久久久久九九精品二区国产| 毛片一级片免费看久久久久| av在线老鸭窝| 久久欧美精品欧美久久欧美| 嫩草影院入口| 99久久精品一区二区三区| 18禁黄网站禁片免费观看直播| 欧美xxxx黑人xx丫x性爽| 成人综合一区亚洲| 最好的美女福利视频网| 天天躁日日操中文字幕| 欧美最黄视频在线播放免费| 亚洲自拍偷在线| 免费观看的影片在线观看| 亚洲国产日韩欧美精品在线观看| 在线观看一区二区三区| 国产欧美日韩一区二区精品| 欧美激情国产日韩精品一区| 免费电影在线观看免费观看| 老司机福利观看| 亚洲精品亚洲一区二区| 禁无遮挡网站| 亚洲经典国产精华液单| 男女边吃奶边做爰视频| 高清毛片免费看| 亚洲精品日韩在线中文字幕 | 午夜免费激情av| 91狼人影院| 深夜精品福利| 男人狂女人下面高潮的视频| 成人特级黄色片久久久久久久| 91久久精品国产一区二区三区| 精品久久久久久久久久久久久| 人人妻,人人澡人人爽秒播| 97热精品久久久久久| 国产伦在线观看视频一区| 免费电影在线观看免费观看| 国产69精品久久久久777片| 国产黄片美女视频| 丰满乱子伦码专区| 免费看av在线观看网站| 全区人妻精品视频| 日本免费a在线| 亚洲国产精品国产精品| 国产男靠女视频免费网站| 97热精品久久久久久| 国产亚洲精品av在线| 亚洲av.av天堂| 国产69精品久久久久777片| 看十八女毛片水多多多| 成人特级av手机在线观看| 麻豆精品久久久久久蜜桃| 99视频精品全部免费 在线| 少妇裸体淫交视频免费看高清| 久久人妻av系列| 国产一级毛片七仙女欲春2| 免费人成视频x8x8入口观看| 国产激情偷乱视频一区二区| avwww免费| 别揉我奶头~嗯~啊~动态视频| 亚洲av五月六月丁香网| 一进一出好大好爽视频| 国产精品久久视频播放| 国产精品一二三区在线看| 国产 一区 欧美 日韩| 欧美日韩国产亚洲二区| 三级国产精品欧美在线观看| 97碰自拍视频| 一级a爱片免费观看的视频| 久久久久久久午夜电影| 色视频www国产| 在线观看免费视频日本深夜| 欧美一区二区国产精品久久精品| 欧美一区二区亚洲| 亚洲国产欧洲综合997久久,| 男人狂女人下面高潮的视频| 国产免费一级a男人的天堂| 日韩av在线大香蕉| 在线国产一区二区在线| 97超视频在线观看视频| 欧美极品一区二区三区四区| 国产爱豆传媒在线观看| 一边摸一边抽搐一进一小说| 亚洲国产日韩欧美精品在线观看| 国产精品一区二区三区四区久久| 欧洲精品卡2卡3卡4卡5卡区| 亚洲精品一区av在线观看| 国产欧美日韩一区二区精品| 婷婷亚洲欧美| 一级毛片久久久久久久久女| av专区在线播放| 男人舔女人下体高潮全视频| 一级黄片播放器| 亚洲第一区二区三区不卡| 99国产精品一区二区蜜桃av| 日韩欧美在线乱码| 精品乱码久久久久久99久播| 人妻久久中文字幕网| 久久久午夜欧美精品| 亚洲欧美中文字幕日韩二区| 中文亚洲av片在线观看爽| 免费看美女性在线毛片视频| 久久久精品大字幕| 最好的美女福利视频网| 人妻少妇偷人精品九色| 99视频精品全部免费 在线| 最近中文字幕高清免费大全6| 国产在线男女| 久久中文看片网| 高清毛片免费看| 亚洲精品成人久久久久久| 午夜激情福利司机影院| 亚洲中文字幕日韩| 欧美性猛交╳xxx乱大交人| 久久久精品大字幕| 亚洲美女视频黄频| 欧美成人精品欧美一级黄| 日日啪夜夜撸| 国产人妻一区二区三区在| 国产精品一区二区三区四区免费观看 | 超碰av人人做人人爽久久| 午夜免费男女啪啪视频观看 | 人妻丰满熟妇av一区二区三区| 91麻豆精品激情在线观看国产| 蜜桃亚洲精品一区二区三区| 日本与韩国留学比较| 久久久久性生活片| 大型黄色视频在线免费观看| 国产淫片久久久久久久久| 亚洲欧美精品自产自拍| 午夜老司机福利剧场| 伦理电影大哥的女人| 亚洲七黄色美女视频| 日本熟妇午夜| 午夜亚洲福利在线播放| 国产成人一区二区在线| 免费电影在线观看免费观看| 高清毛片免费看| 亚洲精品成人久久久久久| 夜夜夜夜夜久久久久| 久久久午夜欧美精品| 三级毛片av免费| 午夜福利在线在线| aaaaa片日本免费| av在线观看视频网站免费| 一级黄片播放器| 神马国产精品三级电影在线观看| 别揉我奶头~嗯~啊~动态视频| 国产高潮美女av| 麻豆成人午夜福利视频| 国内少妇人妻偷人精品xxx网站| 国产精品一区www在线观看| 午夜福利在线观看免费完整高清在 | 九九久久精品国产亚洲av麻豆| 成年av动漫网址| 国产高清激情床上av| 免费黄网站久久成人精品| 久久国产乱子免费精品| 一级黄片播放器| 自拍偷自拍亚洲精品老妇| 日本撒尿小便嘘嘘汇集6| av中文乱码字幕在线| 国产片特级美女逼逼视频| 国内精品美女久久久久久| 国产欧美日韩一区二区精品| 身体一侧抽搐| 精品人妻熟女av久视频| 日本免费一区二区三区高清不卡| 久久午夜亚洲精品久久| 亚洲中文字幕一区二区三区有码在线看| 日本一二三区视频观看| 免费观看在线日韩| 亚洲四区av| 国产免费一级a男人的天堂| 国产片特级美女逼逼视频| 国产高清视频在线播放一区| 插阴视频在线观看视频| h日本视频在线播放| 一个人观看的视频www高清免费观看| 日韩国内少妇激情av| 欧美zozozo另类| 成年版毛片免费区| 一边摸一边抽搐一进一小说| 亚洲国产高清在线一区二区三| 男插女下体视频免费在线播放| 国产亚洲精品久久久久久毛片| 最近手机中文字幕大全| 国产一级毛片七仙女欲春2| 久久久久九九精品影院| 国产熟女欧美一区二区| 亚洲最大成人手机在线| 12—13女人毛片做爰片一| 狂野欧美白嫩少妇大欣赏| 国产午夜精品久久久久久一区二区三区 | 别揉我奶头~嗯~啊~动态视频| 日本免费一区二区三区高清不卡| 三级毛片av免费| 国产一区二区三区av在线 | 亚洲不卡免费看| 色噜噜av男人的天堂激情| 又爽又黄a免费视频| 亚洲综合色惰| 淫妇啪啪啪对白视频| 亚洲国产欧美人成| 日韩在线高清观看一区二区三区| 观看美女的网站| 色综合色国产| 一边摸一边抽搐一进一小说| 身体一侧抽搐| 1024手机看黄色片| 此物有八面人人有两片| 亚洲成a人片在线一区二区| 日韩欧美精品v在线| 亚洲欧美精品综合久久99| 成人国产麻豆网| 国产 一区 欧美 日韩| 国产亚洲91精品色在线| 久久久欧美国产精品| 综合色丁香网| 免费无遮挡裸体视频| 18禁裸乳无遮挡免费网站照片| 精品一区二区三区视频在线观看免费| 国产亚洲av嫩草精品影院| 久久久久久九九精品二区国产| 99热全是精品| 五月伊人婷婷丁香| 国内久久婷婷六月综合欲色啪| 国产私拍福利视频在线观看| 成人欧美大片| 国产精品女同一区二区软件| 熟妇人妻久久中文字幕3abv| 嫩草影院入口| 日本黄大片高清| 免费在线观看成人毛片| 少妇被粗大猛烈的视频| 国产大屁股一区二区在线视频| 少妇被粗大猛烈的视频| 成熟少妇高潮喷水视频| 亚洲av.av天堂|