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

    井壁穩(wěn)定問(wèn)題的不確定性分析方法探討

    2017-11-28 08:24:36勝亞楠管志川許玉強(qiáng)王慶張波
    斷塊油氣田 2017年6期
    關(guān)鍵詞:蒙特卡洛概率分布標(biāo)準(zhǔn)差

    勝亞楠,管志川,許玉強(qiáng),王慶,張波

    (中國(guó)石油大學(xué)(華東)石油工程學(xué)院,山東 青島 266580)

    井壁穩(wěn)定問(wèn)題的不確定性分析方法探討

    勝亞楠,管志川,許玉強(qiáng),王慶,張波

    (中國(guó)石油大學(xué)(華東)石油工程學(xué)院,山東 青島 266580)

    井壁穩(wěn)定性失效是鉆井工程中主要的復(fù)雜情況之一,給鉆井作業(yè)帶來(lái)了很大的困難和安全隱患。針對(duì)井壁穩(wěn)定性問(wèn)題,文中探討了4種不確定性分析方法:蒙特卡洛(Monte-Carlo)模擬法、統(tǒng)計(jì)誤差分析法、一次二階矩法及一次可靠度法。分析發(fā)現(xiàn),基于蒙特卡洛模擬法得到的結(jié)果更為保守,更有利于保持井壁穩(wěn)定性。同時(shí),對(duì)影響井壁問(wèn)題的各個(gè)因素進(jìn)行了敏感性分析,獲取了井壁穩(wěn)定性計(jì)算模型中最為敏感的因素:對(duì)于坍塌壓力最為敏感的因素是最大水平地應(yīng)力,對(duì)于破裂壓力最為敏感的因素是最小水平地應(yīng)力。研究結(jié)果有助于確定影響井壁穩(wěn)定可靠度最為敏感的因素,并運(yùn)用一定方法提高敏感因素的準(zhǔn)確程度,最終減少坍塌及破裂壓力預(yù)測(cè)結(jié)果的不確定性,從而提高井壁穩(wěn)定的可靠度。

    井壁穩(wěn)定性;不確定性分析;蒙特卡洛模擬法;統(tǒng)計(jì)誤差分析法;一次二階矩法;一次可靠度分析法;敏感性分析

    0 引言

    油氣井工程是一項(xiàng)隱蔽的地下工程,存在著大量的隨機(jī)性、模糊性和不確定性,也是一項(xiàng)高風(fēng)險(xiǎn)作業(yè),其中,井壁穩(wěn)定性失效是鉆井工程中主要的復(fù)雜情況與事故之一[1-3]。地層坍塌及破裂壓力是保持井壁穩(wěn)定的安全鉆井液密度窗口的上、下限,地質(zhì)環(huán)境的復(fù)雜性、測(cè)井及地震解釋資料的不完備性以及分析模型的區(qū)域適用性等,都會(huì)導(dǎo)致地層坍塌及破裂壓力預(yù)測(cè)結(jié)果存在著不確定性[4-6]?;诓淮_定性的預(yù)測(cè)結(jié)果進(jìn)行鉆井設(shè)計(jì),會(huì)導(dǎo)致不可預(yù)知的風(fēng)險(xiǎn),如井壁失穩(wěn),影響鉆井安全高效施工。因此,有必要對(duì)于井壁穩(wěn)定問(wèn)題進(jìn)行不確定性分析。本文探討了4種常用的不確定性分析方法:蒙特卡洛模擬法(Monte-Carlo Method)、統(tǒng)計(jì)誤差分析法(S.E.A.M)、一次二階矩法(F.O.S.M)及一次可靠度分析方法(F.O.R.M)。分析發(fā)現(xiàn),蒙特卡洛模擬法更適合用于分析井壁穩(wěn)定問(wèn)題。同時(shí),對(duì)于影響井壁穩(wěn)定的各因素進(jìn)行了敏感性分析。

    1 不確定性分析方法

    1.1 蒙特卡洛模擬法

    該方法是求解工程技術(shù)問(wèn)題近似解的一種數(shù)值計(jì)算方法。它以統(tǒng)計(jì)抽樣理論為基礎(chǔ),利用隨機(jī)數(shù)進(jìn)行抽樣實(shí)驗(yàn)或隨機(jī)模擬,通過(guò)統(tǒng)計(jì)有關(guān)隨機(jī)變量,求得統(tǒng)計(jì)特征值,作為待解問(wèn)題的數(shù)值解[7-8]?;谠摲椒ǖ牡貙悠屏鸭疤鷫毫Σ淮_定性量化表征主要分為3個(gè)步驟:

    1)通過(guò)統(tǒng)計(jì)分析已鉆井巖石地質(zhì)力學(xué)參數(shù)的測(cè)井解釋結(jié)果,由概率統(tǒng)計(jì)原理得到其概率分布擬合函數(shù)。

    2)根據(jù)地質(zhì)力學(xué)參數(shù)的概率分布,運(yùn)用蒙特卡洛模擬方法,將模擬次數(shù)設(shè)定為N,生成N個(gè)符合地質(zhì)力學(xué)參數(shù)概率分布的隨機(jī)數(shù),并代入地層坍塌壓力和地層破裂壓力計(jì)算模型中,得到N個(gè)計(jì)算結(jié)果。

    3)統(tǒng)計(jì)分析計(jì)算結(jié)果,基于正態(tài)信息擴(kuò)散原理得到地層坍塌壓力和地層破裂壓力的概率分布,獲取其概率分布特征值:均值、方差及變異系數(shù)。

    井漏與井壁坍塌風(fēng)險(xiǎn)概率計(jì)算公式如表1所示。表中:PRL為井漏風(fēng)險(xiǎn)概率;PRC為井壁坍塌風(fēng)險(xiǎn)概率;ρL為防井漏鉆井液密度上限;ρC為防井壁坍塌鉆井液密度下限;ρd為鉆井液密度;F(ρd>ρL)為鉆井液密度大于防井漏鉆井液密度上限的概率;F(ρd<ρC)為鉆井液密小于防井壁坍塌鉆井液密度下限的概率。

    表1 井漏和井壁坍塌風(fēng)險(xiǎn)模型

    1.2 統(tǒng)計(jì)誤差分析法

    設(shè)pw為地層坍塌或地層破裂壓力,其不確定性的計(jì)算結(jié)果取決于n個(gè)不確定性的地質(zhì)力學(xué)參數(shù)X1,X2,X3,…,Xn,可以表示為

    當(dāng)pw發(fā)生Δpw的變化時(shí),可以用n個(gè)變量的變異函數(shù)表示為

    假定各個(gè)變量之間相互獨(dú)立,并且均符合正態(tài)分布概率,則式(2)可寫(xiě)成:

    式中:σ(pw),σ(Xi)分別為pw,Xi的標(biāo)準(zhǔn)差。

    式(3)還可以寫(xiě)成用變異函數(shù)來(lái)表示的形式:

    式中:Var(pw),Var(Xi)分別為pw,Xi的變異系數(shù)。

    由式(3)和式(4)可知,地層坍塌壓力或破裂壓力計(jì)算的不確定性可以用標(biāo)準(zhǔn)差來(lái)表示,二者的均值和標(biāo)準(zhǔn)差是根據(jù)n個(gè)輸入?yún)?shù)的均值和標(biāo)準(zhǔn)差計(jì)算獲得的。定義為當(dāng)輸入?yún)?shù)取值為均值加上一個(gè)標(biāo)準(zhǔn)差得到的地層坍塌壓力或破裂壓力計(jì)算結(jié)果,為當(dāng)輸入?yún)?shù)取值為均值減去一個(gè)標(biāo)準(zhǔn)差得到的地層坍塌壓力或破裂壓力計(jì)算結(jié)果,為輸入?yún)?shù)為均值時(shí)得到的地層坍塌壓力或破裂壓力計(jì)算結(jié)果。則地層坍塌壓力或破裂壓力的變異系數(shù)可以表示為

    當(dāng)井底壓力為pa時(shí),井壁穩(wěn)定可靠度系數(shù)Z為

    1.3 一次二階矩法

    該方法[9-10]是求解非時(shí)變荷載作用下結(jié)構(gòu)可靠度問(wèn)題的行之有效的近似方法,既有較高的精度,又有較高的計(jì)算效率。它是在基本變量 Xi(i=1,2,3,…,n)的概率分布尚不清楚時(shí),采用只有均值(又稱為一階原點(diǎn)矩)和標(biāo)準(zhǔn)差(又稱為二階中心矩)的數(shù)學(xué)模型去求解結(jié)構(gòu)可靠度的方法。當(dāng)以一次二階矩法估算工程結(jié)構(gòu)可靠度時(shí),可靠指標(biāo)直接和基本變量的平均值和標(biāo)準(zhǔn)差有關(guān),故其基本概括了各有關(guān)變量的統(tǒng)計(jì)特性,比較全面地反映了各種影響因素的變異性,這是傳統(tǒng)的用安全系數(shù)來(lái)評(píng)價(jià)工程結(jié)構(gòu)安全度的方法所不能做到的。同時(shí),可靠指標(biāo)是從結(jié)構(gòu)功能函數(shù)求解的,綜合考慮了結(jié)構(gòu)上的荷載和結(jié)構(gòu)本身抗力的變異性對(duì)結(jié)構(gòu)可靠度的影響。

    定義地層坍塌壓力和破裂壓力對(duì)應(yīng)的安全系數(shù):

    式中:FSCR,F(xiàn)SFR分別為地層坍塌壓力、破裂壓力安全系數(shù);pCR,pFR分別為地層坍塌壓力、破裂壓力,MPa。

    定義一次二階矩功能函數(shù)G(X)為

    式中:FS為地層坍塌或破裂壓力安全系數(shù)。

    可靠性指標(biāo) β,定義為功能函數(shù)均值E(G (X))與標(biāo)準(zhǔn)差 σ(G (X))的比值:

    式(9)也可用變異系數(shù)表示:

    其中,安全系數(shù)的變異系數(shù)Var F()S為

    定義輸入?yún)?shù)的改變量為±δ(Xi),等于每個(gè)輸入?yún)?shù)標(biāo)準(zhǔn)差的5%。定義為當(dāng)輸入?yún)?shù)取值為均值加上δ(Xi)時(shí),得到的地層坍塌壓力或破裂壓力對(duì)應(yīng)的安全系數(shù):為當(dāng)輸入?yún)?shù)取值為均值減去δ(Xi)時(shí),得到的地層坍塌壓力或破裂壓力對(duì)應(yīng)的安全系數(shù);為輸入?yún)?shù)為均值時(shí)得到的地層坍塌壓力或破裂壓力對(duì)應(yīng)的安全系數(shù)。則地層坍塌壓力或破裂壓力對(duì)應(yīng)的安全系數(shù)的變異系數(shù)可以表示為

    井壁失效概率PR定義為

    式中:Φ(β)為標(biāo)準(zhǔn)正態(tài)分布函數(shù)。

    1.4 一次可靠度分析方法

    該方法的基本思路是[11-12]:首先將一次可靠度分析功能函數(shù) Z′=g(X1,X2,X3,…,Xn)展開(kāi)成Taylor級(jí)數(shù),忽略高階項(xiàng),僅保留線性項(xiàng),再根據(jù)隨機(jī)變量X=(X1,X2,X3,…,Xn)的一階矩求取功能函數(shù)Z′的均值與標(biāo)準(zhǔn)差,從而確定結(jié)構(gòu)的可靠度指標(biāo)。根據(jù)功能函數(shù)線性化點(diǎn)的取法不同以及是否考慮基本隨機(jī)變量的分布類型。

    將功能函數(shù)在 Xi(i=1,2,3,…,n )的均值點(diǎn) μXi展開(kāi)成Taylor級(jí)數(shù),僅保留線性項(xiàng),有:

    因此,Z′的均值 μZ′、方差 σZ′為

    井壁穩(wěn)定可靠性指標(biāo)為

    由式(13)計(jì)算可得井壁失穩(wěn)概率。

    2 敏感性分析

    敏感性分析[13-15]是不確定性分析中重要的一個(gè)環(huán)節(jié)。敏感性分析的目的是找出對(duì)于結(jié)果影響最為敏感的隨機(jī)變量,從而有針對(duì)性地提高敏感因素的準(zhǔn)確性,對(duì)于那些不敏感因素,為了提高計(jì)算效率,對(duì)于其準(zhǔn)確性的要求降低,最終目的是快速高效地得到較為準(zhǔn)確的計(jì)算結(jié)果。本文敏感性分析的思路是,對(duì)于n個(gè)隨機(jī)變量,每一次只把1個(gè)變量看作隨機(jī)值,其他變量保持定值;然后,通過(guò)分析每一個(gè)隨機(jī)變量與結(jié)果之間的相互關(guān)系,最終得到最為敏感的因素。

    定義:R+為當(dāng)輸入?yún)?shù)在均值的基礎(chǔ)上增加時(shí)得到的輸出結(jié)果;R-為當(dāng)輸入?yún)?shù)在均值的基礎(chǔ)上減小時(shí)得到的輸出結(jié)果;Ro為輸入?yún)?shù)為均值時(shí)得到的輸出結(jié)果;Δ (R ,Xi)為Xi對(duì)輸出結(jié)果的影響大小。其關(guān)系表達(dá)式為

    3 實(shí)例分析與結(jié)果討論

    本文分析各個(gè)地質(zhì)力學(xué)參數(shù)的均值、標(biāo)準(zhǔn)差如表2所示,分布形式均為正態(tài)分布N (μ,σ2)。根據(jù)表2各個(gè)參數(shù)設(shè)定的概率分布,基于蒙特卡洛模擬,隨機(jī)生成4000個(gè)符合概率分布特征的隨機(jī)數(shù),代入井壁穩(wěn)定性計(jì)算公式,得到地層坍塌及破裂壓力計(jì)算結(jié)果。統(tǒng)計(jì)分析計(jì)算結(jié)果,并基于正態(tài)信息擴(kuò)散原理得到坍塌壓力和破裂壓力的概率分布、均值、方差及變異系數(shù)。

    由表1定義的風(fēng)險(xiǎn)模型,得到井漏和井壁坍塌風(fēng)險(xiǎn)。根據(jù)本文建立的其他3種不確定性分析模型的計(jì)算公式,結(jié)合各個(gè)參數(shù)的概率分布特征參數(shù),得到不同不確定性分析方法下的井漏和井壁坍塌風(fēng)險(xiǎn),結(jié)果如圖1所示。

    表2 模型中各參數(shù)概率分布形式及特征參數(shù)

    圖1 不同不確定性分析法得到的井壁坍塌和井漏風(fēng)險(xiǎn)概率

    通過(guò)對(duì)比圖1中各種不確定性分析方法可知,基于蒙特卡洛模擬方法得到的井漏及井壁坍塌風(fēng)險(xiǎn)相比于其他3種方法得到的結(jié)果都偏大。也就是說(shuō),基于蒙特卡洛得到的結(jié)果更為保守,更有利于保持井壁穩(wěn)定性。因此,推薦選用蒙特卡洛模擬方法對(duì)井壁問(wèn)題進(jìn)行不確定性分析。

    通過(guò)蒙特卡洛模擬,分別得到不同地質(zhì)力學(xué)參數(shù)概率分布下的坍塌壓力和破裂壓力的概率分布,得到其均值、標(biāo)準(zhǔn)差及變異系數(shù)。定量計(jì)算結(jié)果表明,垂向應(yīng)力為最大主應(yīng)力,坍塌壓力、破裂壓力當(dāng)量鉆井液密度分別為1.0717,2.4210 g/m3。

    對(duì)基于蒙特卡洛法的井壁穩(wěn)定可靠度隨機(jī)變量敏感性進(jìn)行分析,結(jié)果如圖2所示。敏感性分析結(jié)果如下:地層破裂壓力影響因素敏感性由大到小依次為最小水平地應(yīng)力、最大水平地應(yīng)力、孔隙壓力、巖石抗拉強(qiáng)度;地層坍塌壓力影響因素敏感性由大到小依次為最大水平地應(yīng)力、孔隙壓力、內(nèi)聚力、內(nèi)摩擦角、最小水平地應(yīng)力。

    圖2 不同地質(zhì)力學(xué)參數(shù)對(duì)應(yīng)的地層坍塌、破裂壓力變異系數(shù)

    4 結(jié)論

    1)通過(guò)優(yōu)選4種不確定性分析方法,最終確定蒙特卡洛模擬方法更加適合于分析井壁穩(wěn)定問(wèn)題。

    2)地層破裂壓力影響因素敏感性由大到小依次為最小水平地應(yīng)力、最大水平地應(yīng)力、孔隙壓力、巖石抗拉強(qiáng)度;地層坍塌壓力影響因素敏感性由大到小依次為最大水平地應(yīng)力、孔隙壓力、內(nèi)聚力、內(nèi)摩擦角、最小水平地應(yīng)力。

    3)獲取了井壁穩(wěn)定性計(jì)算模型中最為敏感的因素,從而有助于通過(guò)一定的方法和手段提高敏感因素的準(zhǔn)確程度。而對(duì)于那些不敏感因素,則可以根據(jù)井眼或者已鉆井的統(tǒng)計(jì)數(shù)據(jù)估計(jì)其范圍,最終減少坍塌及破裂壓力預(yù)測(cè)結(jié)果的不確定性,提高井壁穩(wěn)定可靠度。

    [1]陳庭根,管志川.鉆井工程理論與技術(shù)[M].東營(yíng):石油大學(xué)出版社,2000:251-254.

    [2]蔣希文.鉆井事故與復(fù)雜問(wèn)題[M].第2版.北京:石油工業(yè)出版社,2006:35-40.

    [3]OTTESEN S,ZHENG R H,MCCANNR C.Borehole stability assessment using quantitative risk analysis[R].SPE 52864,1999.

    [4]羅黎敏,勝亞楠,劉曉坡,等.含可信度地層破裂壓力的鉆前預(yù)測(cè)方法[J].斷塊油氣田,2017,24(1):112-115.

    [5]MOOS D,PESKA P,F(xiàn)INKBEINER T,et al.Comprehensive wellbore stability analysis utilizing quantitative risk assessment[J].Journal of Petroleum Scienceamp;Engineering,2003,38(3):97-109.

    [6]MOSTAFAVI V,ADNOY B S,HARELAND G.Model-based uncertainty assessment of wellbore stability analyses and downhole pressure estimations[J].Journal of American Rock Mechanics,2011,12(5):60-65.

    [7]宮鳳強(qiáng),侯尚騫,巖小明.基于正態(tài)信息擴(kuò)散原理的Mohr-Coulomb強(qiáng)度準(zhǔn)則參數(shù)概率模型推斷方法[J].巖石力學(xué)與工程學(xué)報(bào),2013,32(11):2225-2234.

    [8]茆詩(shī)松,程依明,濮曉龍.概率論與數(shù)理統(tǒng)計(jì)教程[M].北京:高等教育出版社,2004:77-80.

    [9]HARR M E.Reliability-based design in civil engineering[M].New York:McGraw-Hill,1988:103-106.

    [10]UDEGBUNAM J E ,AADNOY B S,F(xiàn)JELDE K K.Uncertainty evaluation of wellbore stability model predictions[J].Journal of Petroleum Scienceamp;Engineering,2014,124(6):254-263.

    [11]ANG H S.Probability concepts in engineering planning and design:decision,risk,and reliability[M].New Jersey :Wiley,1975:91-94.

    [12]SPIEGEL M R ,SCHILLER J J,SRINIVASAN A.Probability and statistics:based on Schaum′s outline of probability and statistics[M].New York:McGraw-hill,2001:63-69.

    [13]SHENG Y,REDDISH D,LU Z.Assessment of uncertainties in wellbore stability analysis[J].Modern Trends in Geomechanics,2006,106(6):541-557.

    [14]PAYTEN W,LAW M.Estimating the plastic collapse of pressure vessels using plasticity contours[J].International Journal of Pressure Vesselsamp;Piping,1998,75(7):529-536.

    [15]SANTOS H,PLACIDO JCR,WOLTER C.Consequences and relevance of drilling vibration on wellbore stability[R].SPE 52820,1999.

    (編輯 王淑玉)

    Discussion on uncertainty analysis method for wellbore stability

    SHENG Yanan,GUAN Zhichuan,XU Yuqiang,WANG Qing,ZHANG Bo
    (College of Petroleum Engineering,China University of Petroleum,Qingdao 266580,China)

    Failure of wellbore stability is one of the main complex situations in drilling engineering,which brings great difficulties and security risks to the drilling operation.In this paper,the wellbore stability reliability was analyzed based on the Monte-Carlo simulation,S.E.A.M,F.O.S.M,and F.O.R.M methods.The outcomes of this study show that wellbore instability risks based on Monte-Carlo simulation were more conservative.Therefore,the Monte Carlo simulation is more conducive to the wellbore stability reliability.In addition,sensitivity of random variables of wellbore stability reliability was analyzed based on the Monte-Carlo simulation method.The results show that the maximum horizontal stress is the most sensitive parameter for collapse pressure;the minimum horizontal stress has the largest effects on formation fracture pressure.The outcomes of this study made it possible for risk reducing analysis by focusing on the most sensitive parameters which are critical for the wellbore stability.The accuracy of sensitive parameters is improved to reduce the uncertainties of the collapse and fracture pressure prediction,and the reliability of wellbore stability is improved.

    wellbore stability;uncertainty analysis;Monte-Carlo simulation;S.E.A.M;F.O.S.M;F.O.R.M;sensitivity analysis

    國(guó)家自然科學(xué)基金面上項(xiàng)目“隔水管內(nèi)氣液兩相流聲傳播特性與深水鉆井氣侵監(jiān)測(cè)方法”(51574275);長(zhǎng)江學(xué)者及創(chuàng)新團(tuán)隊(duì)發(fā)展計(jì)劃資助項(xiàng)目(IRT_14R58)

    TE21

    A

    10.6056/dkyqt201706025

    2017-05-16;改回日期:2017-09-27。

    勝亞楠,男,1989年生,在讀博士研究生,主要從事鉆井工程風(fēng)險(xiǎn)評(píng)價(jià)、井身結(jié)構(gòu)優(yōu)化設(shè)計(jì)和油氣井井下信息控制等方面的研究工作。E-mail:shengyanan_upc@163.com。

    管志川,1959年生,教授,博士生導(dǎo)師,主要從事油氣井管柱力學(xué)、井下信息傳輸技術(shù)、深井超深井井身結(jié)構(gòu)設(shè)計(jì)及提速工具的研發(fā)和深井風(fēng)險(xiǎn)評(píng)價(jià)等方面的研究。電話:0532-86981764;E-mail:guanzhch@163.com。

    勝亞楠,管志川,許玉強(qiáng),等.井壁穩(wěn)定問(wèn)題的不確定性分析方法探討[J].斷塊油氣田,2017,24(6):847-850,866.

    SHENG Yanan,GUAN Zhichuan,XU Yuqiang,et al.Discussion on uncertainty analysis method for wellbore stability[J].Fault-Block Oilamp;Gas Field,2017,24(6):847-850,866.

    猜你喜歡
    蒙特卡洛概率分布標(biāo)準(zhǔn)差
    用Pro-Kin Line平衡反饋訓(xùn)練儀對(duì)早期帕金森病患者進(jìn)行治療對(duì)其動(dòng)態(tài)平衡功能的影響
    離散型概率分布的ORB圖像特征點(diǎn)誤匹配剔除算法
    征服蒙特卡洛賽道
    關(guān)于概率分布函數(shù)定義的辨析
    科技視界(2016年19期)2017-05-18 10:18:46
    利用控制變量方法縮減蒙特卡洛方差
    基于概率分布的PPP項(xiàng)目風(fēng)險(xiǎn)承擔(dān)支出測(cè)算
    蒙特卡洛模擬法計(jì)算電動(dòng)汽車(chē)充電負(fù)荷
    基于蒙特卡洛的非線性約束條件下的優(yōu)化算法研究
    對(duì)于平均差與標(biāo)準(zhǔn)差的數(shù)學(xué)關(guān)系和應(yīng)用價(jià)值比較研究
    一種相依極小P值統(tǒng)計(jì)量概率分布的近似計(jì)算方法
    自拍欧美九色日韩亚洲蝌蚪91 | 麻豆成人av视频| 亚洲成人一二三区av| 中文乱码字字幕精品一区二区三区| 丝袜美腿在线中文| 中文精品一卡2卡3卡4更新| 菩萨蛮人人尽说江南好唐韦庄| av网站免费在线观看视频| 国产有黄有色有爽视频| 伦精品一区二区三区| 国产精品人妻久久久久久| 91狼人影院| 国产成人aa在线观看| 最后的刺客免费高清国语| 天天躁日日操中文字幕| 成人国产av品久久久| 尾随美女入室| 99热这里只有是精品50| 极品少妇高潮喷水抽搐| 18禁动态无遮挡网站| 久久久久性生活片| 插阴视频在线观看视频| 爱豆传媒免费全集在线观看| 国产精品av视频在线免费观看| 国产精品成人在线| 2022亚洲国产成人精品| 国产成人精品久久久久久| 真实男女啪啪啪动态图| 欧美成人午夜免费资源| 99久久人妻综合| 极品少妇高潮喷水抽搐| 国产精品麻豆人妻色哟哟久久| 老师上课跳d突然被开到最大视频| 亚洲精品一区蜜桃| 免费黄网站久久成人精品| 亚洲精品一区蜜桃| tube8黄色片| 又黄又爽又刺激的免费视频.| 久久精品国产a三级三级三级| 狠狠精品人妻久久久久久综合| 女的被弄到高潮叫床怎么办| 一级av片app| 午夜精品国产一区二区电影 | 人人妻人人看人人澡| 国产精品不卡视频一区二区| 亚洲人成网站在线观看播放| 日韩欧美一区视频在线观看 | 我的女老师完整版在线观看| 免费播放大片免费观看视频在线观看| 一区二区三区精品91| 哪个播放器可以免费观看大片| 蜜臀久久99精品久久宅男| 免费看a级黄色片| 国产一区二区三区av在线| 亚洲成人中文字幕在线播放| 日本午夜av视频| 老女人水多毛片| 免费在线观看成人毛片| 亚洲精品第二区| 美女高潮的动态| 特大巨黑吊av在线直播| 涩涩av久久男人的天堂| av在线老鸭窝| 亚洲av中文字字幕乱码综合| 18+在线观看网站| 91狼人影院| 搞女人的毛片| 波野结衣二区三区在线| 亚洲av日韩在线播放| 欧美+日韩+精品| 国产毛片在线视频| 91狼人影院| 免费黄频网站在线观看国产| 在线亚洲精品国产二区图片欧美 | 日本与韩国留学比较| 国产亚洲5aaaaa淫片| 深夜a级毛片| eeuss影院久久| 国内揄拍国产精品人妻在线| 国产成人免费观看mmmm| 中文精品一卡2卡3卡4更新| 精品久久国产蜜桃| 一级二级三级毛片免费看| 久久6这里有精品| videossex国产| 亚洲精品影视一区二区三区av| 亚洲av免费高清在线观看| 欧美xxxx性猛交bbbb| 99热这里只有是精品在线观看| 蜜桃亚洲精品一区二区三区| 亚洲美女搞黄在线观看| 国产日韩欧美亚洲二区| 一级毛片aaaaaa免费看小| 亚洲国产精品999| 亚洲国产日韩一区二区| 欧美日韩综合久久久久久| 亚洲第一区二区三区不卡| 国产免费一级a男人的天堂| 国产精品av视频在线免费观看| 亚洲综合色惰| 日本猛色少妇xxxxx猛交久久| 嫩草影院精品99| 视频中文字幕在线观看| 91久久精品国产一区二区成人| 18禁裸乳无遮挡免费网站照片| 亚洲欧美中文字幕日韩二区| 久久精品国产亚洲网站| 成人二区视频| 又爽又黄a免费视频| 午夜免费观看性视频| 亚洲精品一二三| 久久鲁丝午夜福利片| 国产高清不卡午夜福利| 国产爱豆传媒在线观看| 街头女战士在线观看网站| 九九久久精品国产亚洲av麻豆| 亚州av有码| 狂野欧美激情性xxxx在线观看| 啦啦啦中文免费视频观看日本| 欧美zozozo另类| 日本欧美国产在线视频| 啦啦啦啦在线视频资源| 亚洲不卡免费看| 国产高清三级在线| 亚洲无线观看免费| 最近中文字幕2019免费版| 永久网站在线| 亚洲精品成人久久久久久| 丰满人妻一区二区三区视频av| 久久精品综合一区二区三区| 国产综合精华液| 欧美成人一区二区免费高清观看| 亚洲欧美一区二区三区国产| 啦啦啦啦在线视频资源| 天堂网av新在线| 久久精品国产亚洲av涩爱| 成人亚洲精品一区在线观看 | 丝瓜视频免费看黄片| 成年女人在线观看亚洲视频 | 久久人人爽人人爽人人片va| 少妇裸体淫交视频免费看高清| 国产在视频线精品| 秋霞在线观看毛片| 伊人久久国产一区二区| 97超碰精品成人国产| av福利片在线观看| 久久精品人妻少妇| 看非洲黑人一级黄片| 免费av毛片视频| 久久久久久国产a免费观看| 欧美少妇被猛烈插入视频| 综合色丁香网| 一区二区三区免费毛片| 狂野欧美激情性bbbbbb| 秋霞在线观看毛片| 亚洲欧美中文字幕日韩二区| 免费黄网站久久成人精品| 国产伦精品一区二区三区四那| 亚洲最大成人中文| 日韩av不卡免费在线播放| 黄片wwwwww| 一级a做视频免费观看| 美女内射精品一级片tv| 在线免费观看不下载黄p国产| 免费黄网站久久成人精品| 欧美97在线视频| 最近中文字幕2019免费版| 少妇人妻久久综合中文| 国模一区二区三区四区视频| 亚洲人成网站在线播| 一级爰片在线观看| 人人妻人人爽人人添夜夜欢视频 | 国产免费一级a男人的天堂| eeuss影院久久| 夫妻午夜视频| 亚洲最大成人手机在线| 在线亚洲精品国产二区图片欧美 | 国产高清不卡午夜福利| 纵有疾风起免费观看全集完整版| 亚洲国产高清在线一区二区三| 亚洲天堂av无毛| 有码 亚洲区| 亚洲综合精品二区| 女人被狂操c到高潮| 国产日韩欧美在线精品| 国产毛片a区久久久久| av在线观看视频网站免费| 丰满人妻一区二区三区视频av| 深爱激情五月婷婷| 国产精品久久久久久精品电影小说 | 深夜a级毛片| 男女边吃奶边做爰视频| 亚洲欧美清纯卡通| 18禁动态无遮挡网站| 欧美高清成人免费视频www| 午夜福利在线观看免费完整高清在| 日本午夜av视频| 舔av片在线| 水蜜桃什么品种好| 免费人成在线观看视频色| av天堂中文字幕网| 日韩欧美精品v在线| 综合色av麻豆| 久久久久国产精品人妻一区二区| 免费黄网站久久成人精品| 看黄色毛片网站| 高清午夜精品一区二区三区| 3wmmmm亚洲av在线观看| 亚洲精品国产色婷婷电影| 在线 av 中文字幕| 亚洲av不卡在线观看| 赤兔流量卡办理| 纵有疾风起免费观看全集完整版| 亚洲av男天堂| 国产成人免费无遮挡视频| 亚洲成人一二三区av| 免费看不卡的av| 亚洲精品自拍成人| 国产国拍精品亚洲av在线观看| 日本一二三区视频观看| 直男gayav资源| 亚洲欧美日韩东京热| 国产一区二区三区av在线| 国产成人福利小说| 日本一二三区视频观看| 久久精品国产a三级三级三级| 国产成人freesex在线| freevideosex欧美| 97超视频在线观看视频| 欧美zozozo另类| 少妇 在线观看| 国产探花在线观看一区二区| 国产精品.久久久| 午夜福利在线在线| 美女国产视频在线观看| 女人十人毛片免费观看3o分钟| 水蜜桃什么品种好| 亚洲人成网站高清观看| 夫妻性生交免费视频一级片| 看非洲黑人一级黄片| 久久久久久久精品精品| 秋霞在线观看毛片| 日韩欧美一区视频在线观看 | 超碰av人人做人人爽久久| 男插女下体视频免费在线播放| 成年免费大片在线观看| 亚洲精品一二三| 日本-黄色视频高清免费观看| 久久久久久伊人网av| 亚洲aⅴ乱码一区二区在线播放| 国产69精品久久久久777片| 国产欧美日韩精品一区二区| 熟女电影av网| 欧美日韩精品成人综合77777| 午夜免费鲁丝| 少妇 在线观看| kizo精华| 肉色欧美久久久久久久蜜桃 | 少妇熟女欧美另类| 成人一区二区视频在线观看| 国产精品久久久久久久电影| 少妇人妻 视频| 一个人观看的视频www高清免费观看| 国产精品伦人一区二区| 国产免费又黄又爽又色| 青春草视频在线免费观看| 人妻一区二区av| 国产男女超爽视频在线观看| 精品人妻偷拍中文字幕| 国产男女内射视频| 下体分泌物呈黄色| 国产成人a区在线观看| 日本黄大片高清| 亚洲国产av新网站| 在线观看免费高清a一片| 欧美成人精品欧美一级黄| 亚洲欧美日韩东京热| 看非洲黑人一级黄片| 青青草视频在线视频观看| 欧美成人精品欧美一级黄| 欧美高清性xxxxhd video| 老师上课跳d突然被开到最大视频| 成年免费大片在线观看| 国产男女内射视频| 国产成人一区二区在线| 国产亚洲av嫩草精品影院| 亚洲精品日韩av片在线观看| 欧美日韩在线观看h| 日韩强制内射视频| 少妇熟女欧美另类| 欧美日本视频| 国产亚洲5aaaaa淫片| av一本久久久久| 久久精品国产亚洲网站| 亚洲av不卡在线观看| 精品酒店卫生间| 九九久久精品国产亚洲av麻豆| 一级毛片我不卡| 国产 一区 欧美 日韩| 肉色欧美久久久久久久蜜桃 | 亚洲人成网站在线播| 成人亚洲欧美一区二区av| 国产一区亚洲一区在线观看| 男人和女人高潮做爰伦理| 小蜜桃在线观看免费完整版高清| 亚洲精品日韩在线中文字幕| 欧美国产精品一级二级三级 | 九九久久精品国产亚洲av麻豆| 最近中文字幕高清免费大全6| 高清毛片免费看| av天堂中文字幕网| 七月丁香在线播放| 国产老妇伦熟女老妇高清| 久久亚洲国产成人精品v| 亚洲天堂av无毛| 色播亚洲综合网| 国产综合懂色| 精品人妻偷拍中文字幕| 精品一区在线观看国产| 久久精品久久久久久久性| 久久精品国产亚洲av涩爱| 九九久久精品国产亚洲av麻豆| 免费观看的影片在线观看| 亚洲三级黄色毛片| 国产爱豆传媒在线观看| a级毛片免费高清观看在线播放| 亚洲国产色片| 高清视频免费观看一区二区| 国产一区二区三区av在线| 王馨瑶露胸无遮挡在线观看| 97在线人人人人妻| 欧美精品一区二区大全| 久久影院123| 如何舔出高潮| 99精国产麻豆久久婷婷| 一边亲一边摸免费视频| 一级黄片播放器| 99re6热这里在线精品视频| 久久久久精品久久久久真实原创| 最后的刺客免费高清国语| 最近最新中文字幕大全电影3| 亚洲一区二区三区欧美精品 | 小蜜桃在线观看免费完整版高清| 少妇人妻一区二区三区视频| 国产精品久久久久久av不卡| 天堂网av新在线| 欧美精品人与动牲交sv欧美| 校园人妻丝袜中文字幕| 尤物成人国产欧美一区二区三区| 寂寞人妻少妇视频99o| 夫妻午夜视频| 啦啦啦啦在线视频资源| 夫妻性生交免费视频一级片| 久久人人爽人人爽人人片va| 青春草国产在线视频| 国产毛片a区久久久久| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 成年免费大片在线观看| 一个人看的www免费观看视频| 18禁裸乳无遮挡动漫免费视频 | 久久精品人妻少妇| 久久久久久久精品精品| av线在线观看网站| 亚洲精品aⅴ在线观看| 久久人人爽av亚洲精品天堂 | 熟妇人妻不卡中文字幕| 寂寞人妻少妇视频99o| 综合色av麻豆| 国产欧美日韩一区二区三区在线 | 3wmmmm亚洲av在线观看| 亚洲va在线va天堂va国产| 国产精品人妻久久久久久| 日本午夜av视频| 国产亚洲av嫩草精品影院| 国产又色又爽无遮挡免| 看非洲黑人一级黄片| 亚洲国产色片| 97超碰精品成人国产| 日韩中字成人| 男人和女人高潮做爰伦理| 国产毛片在线视频| 九九久久精品国产亚洲av麻豆| 国产片特级美女逼逼视频| 国产成人午夜福利电影在线观看| 亚洲欧美成人综合另类久久久| 免费黄色在线免费观看| .国产精品久久| 男女边吃奶边做爰视频| 九九爱精品视频在线观看| 国产高潮美女av| 亚洲国产精品成人综合色| 精品视频人人做人人爽| 欧美日韩在线观看h| 亚洲精品国产av成人精品| 日韩在线高清观看一区二区三区| 亚洲国产精品专区欧美| 97精品久久久久久久久久精品| 欧美人与善性xxx| 一本色道久久久久久精品综合| 国产乱来视频区| 免费av不卡在线播放| 久久精品综合一区二区三区| 久热久热在线精品观看| 国产免费一级a男人的天堂| 国产黄片视频在线免费观看| 99热这里只有是精品50| 2021天堂中文幕一二区在线观| 一区二区三区乱码不卡18| 尾随美女入室| av又黄又爽大尺度在线免费看| 久久精品久久久久久久性| 国产成人精品久久久久久| 国产成人a∨麻豆精品| 80岁老熟妇乱子伦牲交| 七月丁香在线播放| 99热全是精品| 精品熟女少妇av免费看| 国产精品精品国产色婷婷| 久久久a久久爽久久v久久| 街头女战士在线观看网站| 国内精品宾馆在线| 亚洲,欧美,日韩| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 国产一区有黄有色的免费视频| 日产精品乱码卡一卡2卡三| 18禁裸乳无遮挡免费网站照片| 欧美老熟妇乱子伦牲交| 中国三级夫妇交换| 男人狂女人下面高潮的视频| 寂寞人妻少妇视频99o| 少妇裸体淫交视频免费看高清| 国产色婷婷99| 亚洲国产av新网站| 97热精品久久久久久| 国产成人午夜福利电影在线观看| 日韩大片免费观看网站| 国产一级毛片在线| 国产高潮美女av| 18禁在线无遮挡免费观看视频| 日韩成人伦理影院| 国产伦精品一区二区三区视频9| 欧美xxxx性猛交bbbb| 免费观看在线日韩| 日本爱情动作片www.在线观看| 男人添女人高潮全过程视频| 视频区图区小说| a级毛色黄片| 午夜福利网站1000一区二区三区| 国产高清不卡午夜福利| 国产有黄有色有爽视频| 色视频在线一区二区三区| 成人欧美大片| 亚洲av男天堂| 国产v大片淫在线免费观看| 亚洲精品亚洲一区二区| 色播亚洲综合网| 国产精品国产av在线观看| 亚洲av免费在线观看| 熟妇人妻不卡中文字幕| 国产精品福利在线免费观看| 精品一区二区三区视频在线| 夫妻性生交免费视频一级片| 久久久久久久精品精品| 欧美高清成人免费视频www| 亚洲国产欧美在线一区| 亚洲不卡免费看| av播播在线观看一区| 狂野欧美激情性bbbbbb| 亚洲欧美一区二区三区黑人 | 久热久热在线精品观看| 中文字幕av成人在线电影| 一级毛片我不卡| 人人妻人人澡人人爽人人夜夜| 超碰av人人做人人爽久久| 国产伦精品一区二区三区四那| 各种免费的搞黄视频| 国产91av在线免费观看| 午夜福利在线在线| 国产成人aa在线观看| 午夜福利在线观看免费完整高清在| 18禁在线无遮挡免费观看视频| 国产一区二区三区综合在线观看 | 天堂中文最新版在线下载 | 中文字幕人妻熟人妻熟丝袜美| 97人妻精品一区二区三区麻豆| 欧美xxxx性猛交bbbb| 国产午夜福利久久久久久| 欧美成人a在线观看| 国产白丝娇喘喷水9色精品| 视频中文字幕在线观看| 国产午夜精品一二区理论片| 在线观看三级黄色| videos熟女内射| 亚洲美女视频黄频| 久热久热在线精品观看| 亚洲真实伦在线观看| 亚洲欧美中文字幕日韩二区| 高清午夜精品一区二区三区| 性插视频无遮挡在线免费观看| 久久鲁丝午夜福利片| 久久ye,这里只有精品| 日韩电影二区| 最近最新中文字幕免费大全7| 欧美成人午夜免费资源| 看非洲黑人一级黄片| 婷婷色av中文字幕| 中文欧美无线码| 麻豆精品久久久久久蜜桃| 国产乱人视频| 最近2019中文字幕mv第一页| 国产女主播在线喷水免费视频网站| 性色avwww在线观看| 黄色日韩在线| 久久97久久精品| 欧美人与善性xxx| 久久久精品欧美日韩精品| 国产又色又爽无遮挡免| 久久久国产一区二区| 国产日韩欧美在线精品| 久久99热这里只有精品18| 亚洲久久久久久中文字幕| 久久精品久久久久久噜噜老黄| 久久久久久久大尺度免费视频| 国产精品福利在线免费观看| a级毛色黄片| 蜜臀久久99精品久久宅男| 国产美女午夜福利| av网站免费在线观看视频| 日韩亚洲欧美综合| 久久女婷五月综合色啪小说 | 精品99又大又爽又粗少妇毛片| 日本午夜av视频| 嘟嘟电影网在线观看| 神马国产精品三级电影在线观看| 日韩亚洲欧美综合| av天堂中文字幕网| 欧美成人一区二区免费高清观看| av福利片在线观看| 一区二区三区免费毛片| 久久久久久久久久人人人人人人| 久久午夜福利片| 黄色日韩在线| 少妇猛男粗大的猛烈进出视频 | 日本av手机在线免费观看| 高清视频免费观看一区二区| 亚洲av一区综合| 中文资源天堂在线| 久久久久久国产a免费观看| 午夜精品一区二区三区免费看| 成人国产av品久久久| 寂寞人妻少妇视频99o| 禁无遮挡网站| 国产成人精品久久久久久| 18禁裸乳无遮挡免费网站照片| 高清午夜精品一区二区三区| 亚洲综合精品二区| 日本-黄色视频高清免费观看| 亚洲四区av| 边亲边吃奶的免费视频| 人妻少妇偷人精品九色| 国产永久视频网站| 亚洲欧美中文字幕日韩二区| 日韩av在线免费看完整版不卡| 99热全是精品| 成人黄色视频免费在线看| 国产精品蜜桃在线观看| 亚洲精品日韩av片在线观看| 国产v大片淫在线免费观看| 男女下面进入的视频免费午夜| 久久精品国产a三级三级三级| 国产成人91sexporn| 在线观看一区二区三区激情| 又粗又硬又长又爽又黄的视频| 久久这里有精品视频免费| 激情 狠狠 欧美| 特级一级黄色大片| 色播亚洲综合网| 国产老妇女一区| 国产老妇伦熟女老妇高清| 69av精品久久久久久| 国产精品一区二区三区四区免费观看| 91精品伊人久久大香线蕉| 国产有黄有色有爽视频| 在现免费观看毛片| 精品国产一区二区三区久久久樱花 | 久久精品国产亚洲av天美| 高清日韩中文字幕在线| 午夜精品一区二区三区免费看| 日本猛色少妇xxxxx猛交久久| 91精品一卡2卡3卡4卡| 久久人人爽av亚洲精品天堂 | 亚洲成色77777| 日韩免费高清中文字幕av| 边亲边吃奶的免费视频| 少妇熟女欧美另类| 国产真实伦视频高清在线观看| 亚洲人成网站在线播| 日韩大片免费观看网站| 中文字幕久久专区| av一本久久久久| 日韩国内少妇激情av| 国产精品久久久久久av不卡| 永久网站在线| 久久久色成人| 午夜日本视频在线| 国产精品久久久久久久久免| 九九久久精品国产亚洲av麻豆| 亚洲在线观看片| 91久久精品国产一区二区三区| 天天一区二区日本电影三级| 好男人视频免费观看在线| 亚洲欧美精品自产自拍| 国产成人aa在线观看| 亚洲欧美日韩无卡精品| 极品教师在线视频|