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

    熱固耦合作用下的套管?水泥環(huán)?地層多層組合系統(tǒng)應(yīng)力分析

    2017-05-19 03:43:56趙新波韓生超楊秀娟閆相禎張立松
    關(guān)鍵詞:內(nèi)層泊松比內(nèi)壁

    趙新波,韓生超,楊秀娟,閆相禎,3,張立松

    熱固耦合作用下的套管?水泥環(huán)?地層多層組合系統(tǒng)應(yīng)力分析

    趙新波1,2,韓生超1,楊秀娟1,閆相禎1,3,張立松1

    (1.中國石油大學(xué)(華東)儲運(yùn)與建筑工程學(xué)院,山東青島,266580;2.青島理工大學(xué)理學(xué)院,山東青島,266033;3.中國石油大學(xué)(華東)油氣CAE技術(shù)研究中心,山東青島,266580)

    考慮非均勻地應(yīng)力、內(nèi)壓和地層高溫對高溫高壓井套管?水泥環(huán)?地層多層組合系統(tǒng)的作用,應(yīng)用彈性力學(xué)理論,結(jié)合邊界條件、接觸條件和連續(xù)條件,獲得熱固耦合作用下套管?水泥環(huán)?地層多層組合系統(tǒng)應(yīng)力分布的解析解。討論水泥環(huán)厚度、彈性模量、泊松比、地層溫度、套管?水泥環(huán)層數(shù)對最內(nèi)層套管內(nèi)壁M ises應(yīng)力分布規(guī)律的影響。研究結(jié)果表明:上述因素對最內(nèi)層套管內(nèi)壁M ises應(yīng)力的分布均有影響;對比套管?水泥環(huán)?地層的單層和多層組合系統(tǒng),多層組合系統(tǒng)的最內(nèi)層套管內(nèi)壁的M ises應(yīng)力較單層系統(tǒng)明顯降低,這說明多層組合系統(tǒng)具有重要的工程價值。

    套管?水泥環(huán)?地層;熱固耦合;多層組合系統(tǒng);彈性力學(xué)

    準(zhǔn)確獲取套管?水泥環(huán)?地層多層組合系統(tǒng)的應(yīng)力分布對于防止套管損壞具有重要工程價值?,F(xiàn)階段,高溫高壓井(HPHT)越來越多,套管處于非均勻高地應(yīng)力、高內(nèi)壓和高地層溫度的聯(lián)合作用下,這進(jìn)一步加劇了套管的損壞,故開展熱固耦合作用下的套管?水泥環(huán)?地層多層組合系統(tǒng)的受力特性分析具有重要價值。國內(nèi)外學(xué)者取得了大量套管?水泥環(huán)?地層三層組合系統(tǒng)的受力特性研究成果。YIN等[1]基于彈性力學(xué)理論,研究了平面問題垂直井筒受地應(yīng)力作用時的力學(xué)行為;殷有泉等[2?5]先后用彈性理論及數(shù)值方法求解得到井筒受地應(yīng)力荷載問題的理論解,并研究了系統(tǒng)接觸表面的應(yīng)力分布規(guī)律;LAST等[6]對構(gòu)造應(yīng)力引起的套管位移進(jìn)行了研究;李軍等[7]用彈性理論及有限元方法對傾斜井組合系統(tǒng)受力進(jìn)行了研究。前人的研究大多基于套管?水泥環(huán)?地層三層組合系統(tǒng),難以將研究結(jié)果直接應(yīng)用于實際的多層套管受力分析。王宴濱等[8]研究了套管?水泥環(huán)?地層多層組合系統(tǒng)受力特性,然而未考慮地層溫度的作用,并不能直接用于高溫高壓井多層組合系統(tǒng)受力特性分析?;诖耍疚淖髡咭蕴坠?水泥環(huán)?地層多層組合系統(tǒng)為對象,應(yīng)用彈性力學(xué)方法,考慮熱固耦合作用,獲取多層組合系統(tǒng)在非均勻地應(yīng)力下的應(yīng)力分布規(guī)律,并討論水泥環(huán)厚度、水泥環(huán)彈性模量、水泥環(huán)泊松比、套管?水泥環(huán)層數(shù)對最內(nèi)層套管內(nèi)壁M ises應(yīng)力分布規(guī)律的影響。套管?水泥環(huán)?地層多層組合系統(tǒng)對于保證最內(nèi)層套管安全具有重要作用。

    1 熱固耦合作用下的多層組合系統(tǒng)力學(xué)模型

    假設(shè)高溫高壓井套管?水泥環(huán)?地層多層組合系統(tǒng)無縱向變形,即不考慮垂向主應(yīng)變的變化,此時套管?水泥環(huán)?地層多層組合系統(tǒng)可簡化為平面應(yīng)變問題。根據(jù)套管?水泥環(huán)?地層的實際服役條件,建立高溫高壓井套管?水泥環(huán)?地層多層組合系統(tǒng)熱固耦合模型如圖1所示(其中,hσ和Hσ分別為最小和最大非均勻地應(yīng)力)。

    在圖1所示的熱固耦合模型中,套管及水泥環(huán)為理想圓形、厚度均勻的各向同性材料。多層組合系統(tǒng)受到的載荷主要為非均勻地應(yīng)力Hσ和hσ、套管內(nèi)壓p0、套管內(nèi)溫度T0與地層溫度Tf的溫差載荷。整個系統(tǒng)共包括n層,其中,第n層為地層,第n?1層為最外層的水泥環(huán),第n?2層為最外層套管,以此類推。第i層的參數(shù)為內(nèi)半徑Ri?1、外半徑Ri、彈性模量Ei、泊松比μi。

    圖1 多層組合系統(tǒng)熱固耦合模型Fig.1 Thermo-structural coup ling analysismodel of multilayer composite system

    2 熱固耦合作用下的多層組合系統(tǒng)應(yīng)力分析

    2.1 邊界應(yīng)力分量的坐標(biāo)變換

    在圖1所示的熱固耦合模型中,需將直角坐標(biāo)變換為極坐標(biāo)進(jìn)行應(yīng)力求解。選取極坐標(biāo)系下的模型外邊界(r=Rn處)足夠遠(yuǎn),即不考慮外邊界溫度變化[9]。因此,由溫度變化引起的內(nèi)外邊界熱應(yīng)力分量為0。利用坐標(biāo)變化[10?11],組合系統(tǒng)應(yīng)力邊界條件可描述為

    2.2 熱固耦合分析

    對于各向同性材料,穩(wěn)態(tài)溫度場及包含熱力學(xué)耦合項的結(jié)構(gòu)場控制方程為[12]

    溫度場邊界條件為

    結(jié)構(gòu)場邊界條件為其中:δ為傳導(dǎo)系數(shù);h為對流換熱系數(shù);t1為溫度;Δt1為變溫;Γ為溫度場區(qū)域邊界;uj為位移向量;fj為體力分量;Ejklm為結(jié)構(gòu)彈性模量張量;t1*為已知溫度;Q*為邊界已知熱流量; u*為已知位移向量;f*為已

    1

    j

    j知外力矢量;εkj,σjk和αjk分別為應(yīng)變、應(yīng)力和熱膨脹系數(shù)張量;應(yīng)力σjk為溫度場t1的函數(shù),下標(biāo)j和k表示對坐標(biāo)的微分;nk為面分量。

    由式(2)中第4式可知,σjk由結(jié)構(gòu)應(yīng)力和溫度應(yīng)力線性疊加得到,因此,熱固耦合問題中的應(yīng)力分量可直接由平均應(yīng)力場和偏差應(yīng)力場引起的結(jié)構(gòu)應(yīng)力以及溫度應(yīng)力疊加得到。

    2.3 平均應(yīng)力場作用下的組合系統(tǒng)結(jié)構(gòu)應(yīng)力分析

    由彈性力學(xué)可得組合系統(tǒng)中第i層圓環(huán)受均勻壓力的拉梅解答[9],其相對于外邊界均勻壓力的無因次表達(dá)式為

    由極坐標(biāo)幾何方程和物理方程[9],并結(jié)合均勻內(nèi)外壓作用下圓環(huán)環(huán)向位移為0的條件,可得第i層圓環(huán)徑向位移的方程為

    由多層組合系統(tǒng)接觸面徑向位移相等可得n?1個方程,形式如下:

    式中:系數(shù)ai,bi和ci為包含材料參數(shù)和幾何尺寸的表達(dá)式;i=1,2,3,…,n?1。

    2.4 偏差應(yīng)力場作用下的組合系統(tǒng)結(jié)構(gòu)應(yīng)力分析

    在偏差應(yīng)力場作用下,通過應(yīng)力函數(shù)求解得到組合系統(tǒng)第i層圓環(huán)相對于外邊界非均勻應(yīng)力的無因次表達(dá)式為

    在式(8)~(10)中,共含有4n個未知量。為此,求解偏差應(yīng)力場下的組合系統(tǒng)應(yīng)力需建立4n個方程。利用偏差應(yīng)力場邊界條件及各層交界面的應(yīng)力和位移連續(xù)性條件可組成4n個方程。

    2.5 熱載荷作用下的組合系統(tǒng)溫度應(yīng)力分析

    2.5.1 溫度場分析

    工程上多數(shù)材料的熱導(dǎo)率λ是溫度的函數(shù)[13],組合系統(tǒng)第i層圓環(huán)的熱傳導(dǎo)率可表示為λi=λ0i(1+eiθ),進(jìn)而可求得第i層圓環(huán)的線熱流量qi:

    式中:θi?1為第i層圓環(huán)的內(nèi)壁溫度;θi為第i層圓環(huán)的外壁溫度;λ0i為第i層圓環(huán)溫度為0℃時材料熱導(dǎo)率;ei為第i層圓環(huán)熱導(dǎo)率隨溫度的變化梯度。

    由能量守恒原理,可得求解組合系統(tǒng)穩(wěn)定溫度場的n?1個方程,形式為qi=qi+1(i=1,2,3,…,n?1),利用試算法[14]可求得n?1個交界面的溫度θi(i=1,2,3,…, n?1),得到其溫度場分布為

    2.5.2 溫度應(yīng)力分析

    通過求解滿足平衡微分方程的位移勢函數(shù),獲得相應(yīng)于位移特解的應(yīng)力分量。由滿足相容條件的應(yīng)力函數(shù)Φ=2?1·Fr2獲得應(yīng)力補(bǔ)充解σr=σφ=F,τrφ=0(F為常數(shù)),最終求得第i層圓環(huán)的溫度應(yīng)力,其相對于內(nèi)壓p0的無因次表達(dá)式為

    其中:Fi和Hi為第i層與力學(xué)參數(shù)有關(guān)的常數(shù);iα為第i層材料的線脹系數(shù);i=1,2,3,…,n。利用極坐標(biāo)下的幾何方程和物理方程,第i層圓環(huán)徑向位移可表示為

    在式(13)和(14)中,共含有2n個未知量。為此,需建立2n個方程求解溫度應(yīng)力。利用溫度場應(yīng)力邊界條件及各層交界面的應(yīng)力和位移連續(xù)性條件可組成2n個方程。

    3 算例分析

    以新疆油田KS-XX井為例,該井地層壓力系數(shù)高、地溫梯度高、地層情況復(fù)雜,為此在2 800m段采用多層組合系統(tǒng)結(jié)構(gòu)。各層管的幾何尺寸、彈性模量、泊松比、線脹系數(shù)和熱導(dǎo)率如表1所示,其中:最大水平主應(yīng)力Hσ為70MPa、最小水平主應(yīng)力σh為62MPa;最內(nèi)層套管內(nèi)徑R0為64.95mm;內(nèi)部壓力p0為28MPa;套管溫度T0為60℃;地層溫度θf為340℃。

    表1 計算參數(shù)Table1 Calculated parameters

    3.1 高溫高壓井多層組合系統(tǒng)應(yīng)力分析

    利用本文推導(dǎo)的理論,編寫程序并采用上述參數(shù)計算組合系統(tǒng)多層組合系統(tǒng)關(guān)鍵位置處的M ises應(yīng)力,計算結(jié)果如圖2所示。

    由圖2可知:在9層的套管?水泥環(huán)?地層多層組合系統(tǒng)中,套管內(nèi)壁的M ises應(yīng)力明顯高于水泥環(huán)內(nèi)壁的M ises應(yīng)力,這與套管及水泥環(huán)的彈性模量相差較大密切相關(guān)。套管M ises應(yīng)力的最大值與最小水平主應(yīng)力hσ方向重合,最小值與最大水平主應(yīng)力σH方向重合。套管M ises應(yīng)力沿井周分布不均勻,呈現(xiàn)正弦曲線分布規(guī)律,這主要源自偏差應(yīng)力場的作用。套管M ises應(yīng)力的非均勻性在第1層和第4層套管內(nèi)壁體現(xiàn)最明顯,應(yīng)力非均勻系數(shù)[8]分別高達(dá)1.678和1.569。對于第9層的套管?水泥環(huán)?地層組合系統(tǒng),第4層套管是系統(tǒng)的最外層套管,直接受到非均勻地應(yīng)力作用,這使得該層套管M ises應(yīng)力非均勻性較強(qiáng)。

    圖2 多層組合系統(tǒng)各層內(nèi)壁的M ises應(yīng)力Fig.2 M ises stress distribution on inner facesofmultilayer com posite system model

    對于第1層套管的內(nèi)壁,直接面臨臨空面,其發(fā)生的應(yīng)力集中必然使得該層套管M ises應(yīng)力非均勻性最高。由于最內(nèi)層套管的應(yīng)力非均勻性較高,使得套管的抗外擠能力大大降低。在第2層和第3層套管內(nèi)壁M ises應(yīng)力非均勻性下降較多,應(yīng)力非均勻系數(shù)分別為1.330和1.389。對于水泥環(huán)內(nèi)壁的M ises應(yīng)力,其應(yīng)力非均勻性明顯降低,這與其彈性模量明顯低于套管的彈性模量直接相關(guān)。

    根據(jù)圖2所示結(jié)果,由于第1層、第4層套管的組合應(yīng)力高于第2層、第3層套管,因此,各層套管應(yīng)分別進(jìn)行鋼級選擇。為了充分發(fā)揮套管管材的力學(xué)性質(zhì),第1層和第4層套管的鋼級應(yīng)高于第2層和第3層套管的鋼級。以新疆油田KS-XX井為例,第1層和第4套管建議選擇N80套管(最小屈服強(qiáng)度為552 MPa),而第2層和第3層套管建議選擇K55套管(最小屈服強(qiáng)度為379MPa)。

    3.2 影響因素分析

    根據(jù)以上分析結(jié)果,在套管?水泥環(huán)?地層組合系統(tǒng)中,最內(nèi)層套管內(nèi)壁的M ises應(yīng)力最大,即最危險區(qū)域為最內(nèi)層套管內(nèi)壁。換句話說,只要最內(nèi)層套管內(nèi)壁的M ises應(yīng)力小于許用應(yīng)力,套管強(qiáng)度即滿足。因此,為了保證套管安全,應(yīng)重點關(guān)注最內(nèi)層套管內(nèi)壁的M ises應(yīng)力。考慮到套管的彈性模量與泊松比基本不變,討論水泥環(huán)厚度、彈性模量、泊松比、地層溫度及套管及水泥環(huán)層數(shù)對最內(nèi)層套管內(nèi)壁M ises應(yīng)力的影響。設(shè)套管壁厚為7.5mm,彈性模量為210 GPa,泊松比為0.25;水泥環(huán)壁厚為25mm,彈性模量為11GPa,泊松比為0.3;地層彈性模量2GPa,泊松比為0.3;地層溫度為340℃,內(nèi)壁溫度為60℃。

    3.2.1 水泥環(huán)厚度

    保持其他參數(shù)不變,改變水泥環(huán)厚度,分析水泥環(huán)厚度對最內(nèi)層套管內(nèi)壁的M ises應(yīng)力的影響。水泥環(huán)厚度變化范圍為10~150mm。

    圖3 水泥環(huán)壁厚對最內(nèi)層套管內(nèi)壁M ises應(yīng)力的影響Fig.3 Influenceof cement thicknesson Misesstress distribution of innermost casing insidewall

    由圖3可知:隨著水泥環(huán)厚度的增加,最內(nèi)層套管內(nèi)壁的M ises應(yīng)力增大,但增長趨勢變小;當(dāng)水泥環(huán)厚度從10mm增加到30mm時,最大M ises應(yīng)力從417.237MPa增加到470.991MPa,增加53.754MPa,然而,當(dāng)水泥環(huán)厚度從130mm增加到150mm時,最大M ises應(yīng)力從656.744MPa增加到677.162MPa,增加20.418MPa;另外,隨著水泥環(huán)厚度增加,最內(nèi)層套管內(nèi)壁M ises應(yīng)力的非均勻性降低。當(dāng)水泥環(huán)厚度為10mm時,非均勻系數(shù)為2.499 5,然而,當(dāng)水泥環(huán)厚度增加到150 mm時,非均勻系數(shù)降低為1.389 3。

    3.2.2 水泥環(huán)彈性模量

    保持其他參數(shù)不變,改變水泥環(huán)彈性模量,分析水泥環(huán)彈性模量對最內(nèi)層套管內(nèi)壁的M ises應(yīng)力的影響。水泥環(huán)彈性模量變化范圍為10~80GPa。

    由圖4可知:隨著水泥環(huán)彈性模量的增加,最內(nèi)層套管內(nèi)壁M ises應(yīng)力呈現(xiàn)先增加后降低的趨勢。在水泥環(huán)彈性模量達(dá)到30 GPa時,最內(nèi)層套管內(nèi)壁M ises應(yīng)力達(dá)到極大值。當(dāng)繼續(xù)增大水泥環(huán)彈性模量,最內(nèi)層套管內(nèi)壁M ises降低。上述分析說明,對于多層組合系統(tǒng),較低的水泥環(huán)彈性模量對于保護(hù)最內(nèi)層套管安全是有益的。

    圖4 水泥環(huán)彈性模量對最內(nèi)層套管內(nèi)壁Mises應(yīng)力的影響Fig.4 Influence of cementelasticity moduluson M ises stress distribution of innermost casing insidew all

    3.2.3 水泥環(huán)泊松比

    保持其他參數(shù)不變,改變水泥環(huán)泊松比,分析水泥環(huán)泊松比對最內(nèi)層套管內(nèi)壁的M ises應(yīng)力的影響。水泥環(huán)泊松比變化范圍為0.1~0.4。

    由圖5可知:隨著水泥環(huán)泊松比增大,最內(nèi)層套管內(nèi)壁M ises應(yīng)力降低,但降低幅度較??;隨著水泥環(huán)泊松比的變化,套管M ises應(yīng)力的非均勻性變化不大。這說明水泥環(huán)泊松比的變化對最內(nèi)層套管內(nèi)壁M ises應(yīng)力的影響較小。

    3.2.4 地層溫度

    保持其他參數(shù)不變,改變地層溫度,分析地層溫度對最內(nèi)層套管內(nèi)壁的M ises應(yīng)力的影響。地層溫度變化范圍為140~490℃。

    圖5 水泥環(huán)泊松比對最內(nèi)層套管內(nèi)壁M ises應(yīng)力的影響Fig.5 Influence of cement Poisson’s ratio on M ises stress distribution of innermost casing insidewall

    圖6 地層溫度對最內(nèi)層套管內(nèi)壁Mises應(yīng)力的影響Fig.6 Influence of formation temperature on M ises stress distribution of innermost casing insidewall

    由圖6可知:隨著地層溫度增加,最內(nèi)層套管內(nèi)壁的M ises應(yīng)力增大。由式(12)可知,距離套管中心相等的位置溫度載荷相等,相應(yīng)的溫度應(yīng)力相等。換句話說,對于距離套管中心相等的位置,溫度應(yīng)力是一個均勻應(yīng)力場,這不會改變最內(nèi)層套管內(nèi)壁M ises應(yīng)力沿井周的變化趨勢,即M ises應(yīng)力沿井周的變化趨勢由非均勻地應(yīng)力控制。

    3.2.5 套管及水泥環(huán)層數(shù)

    水泥環(huán)壁厚15mm,保持其他參數(shù)不變,改變套管及水泥環(huán)層數(shù),分析套管及水泥環(huán)層數(shù)對最內(nèi)層套管內(nèi)壁的M ises應(yīng)力的影響。套管層數(shù)變化范圍為1~4。

    圖7 套管層數(shù)對最內(nèi)層套管內(nèi)壁Mises應(yīng)力的影響Fig.7 Influenceof casing layer numberon Mises stress distribution of innermost casing insidewall

    由圖7可知:隨著多層組合系統(tǒng)層數(shù)增大,最內(nèi)層套管內(nèi)壁的M ises應(yīng)力急劇減小。當(dāng)套管?水泥環(huán)層數(shù)從1增加到4時,M ises應(yīng)力最大值從947.207MPa降低到437.232MPa,降幅509.975MPa;非均勻應(yīng)力系數(shù)從2.863減小到2.119。從以上分析可知,套管?水泥環(huán)?地層多層組合系統(tǒng),對于降低最內(nèi)層套管內(nèi)壁的M ises應(yīng)力及非均勻應(yīng)力系數(shù)具有重要作用,具有工程應(yīng)用的合理性。

    4 結(jié)論

    1)考慮熱固耦合作用,建立了套管?水泥環(huán)?地層多層組合系統(tǒng)的力學(xué)模型,利用彈性力學(xué)理論獲取了套管?水泥環(huán)?地層多層組合系統(tǒng)應(yīng)力分布的解析解。該解析解表明,溫度對套管?水泥環(huán)?地層多層組合系統(tǒng)的當(dāng)量應(yīng)力影響不能忽略。

    2)水泥環(huán)厚度和地層溫度的增加,最內(nèi)層套管內(nèi)壁的M ises應(yīng)力增大;水泥環(huán)彈性模量的增加,最內(nèi)層套管內(nèi)壁M ises應(yīng)力呈現(xiàn)先增大后減小的趨勢;水泥環(huán)泊松比的變化對最內(nèi)層套管內(nèi)壁M ises應(yīng)力的影響較??;多層組合系統(tǒng)層數(shù)增多,最內(nèi)層套管內(nèi)壁的M ises應(yīng)力急劇減小。

    3)在套管?水泥環(huán)?地層多層組合系統(tǒng)中,最內(nèi)層套管和最外層套管的M ises應(yīng)力比中間層的高。該系統(tǒng)對于降低最內(nèi)層套管內(nèi)壁的M ises應(yīng)力、保證服役期間的套管安全具有重要工程價值。

    [1]YIN Fei,GAO Deli.Mechanicalanalysisof casings in boreholes, under non uniform remote crustal stress fields:analytical& numerical methods[J].Computer M odeling in Engineering& Science,2012,89(1):25?37.

    [2]殷有泉,陳朝偉,李平恩.套管?水泥環(huán)?地層應(yīng)力分布的理論解[J].力學(xué)學(xué)報,2006,38(6):835?842.

    YIN Youquan,CHEN Zhaowei,LIPingen.Theoretical solutions of stress distribution in casing-cement and stratum system[J]. Chinese Journal of Theoretical and Applied Mechanics,2006, 38(6):835?842

    [3]房軍,岳伯謙,趙懷文,等.非均勻地應(yīng)力作用下套管和水泥環(huán)表面受力特性分析[J].石油大學(xué)學(xué)報(自然科學(xué)版),1997, 21(1):46?48.

    FANG Jun,YUE Boqian,ZHAO Huaiwen,et al.Analysis of surface loading on casing and cement sheath under non uniform geologic stress[J].Journal of China University of Petroleum (NaturalScience Edition),1997,21(1):46?48

    [4]房軍,趙懷文,岳伯謙,等.非均勻地應(yīng)力作用下套管與水泥環(huán)的受力分析[J].石油大學(xué)學(xué)報(自然科學(xué)版),1995,19(6): 52?57.

    FANG Jun,ZHAN Huaiwen,YUE Boqian,et al.Analysis of sheath loading property of casing and cement under non uniform geologic stress[J].Journal of China University of Petroleum (NaturalScience Edition),1995,19(6):52?57

    [5]房軍,谷玉洪,米豐珍.非均勻載荷作用下套管擠壓失效數(shù)值分析[J].石油機(jī)械,1999,27(7):34?37.

    FANG Jun,GU Yuhong,M IFengzhen.A numerical analysis of casing collapse under non uniform load[J].China Petroleum Machinery,1999,27(7):34?37.

    [6]LAST N C,MUJICA M,PATTILLO P D,et al.Evaluation, impact,and management of casing deformation caused by tectonic forces in Andean Foothills,Colombia[J].SPE Drill Complet,2006,21(2):116?124.

    [7]李軍,陳勉,張輝.定向井套管應(yīng)力隨地應(yīng)力條件的變化規(guī)律研究[J].石油學(xué)報,2005,26(1):109?112.

    LI Jun,CHEN M ian,ZHANG Hui.Study on varying rules of casing stress in directional hole w ith in-situ stress conditions[J]. Acta Petrolei Sinica,2005,26(1):109?112.

    [8]王宴濱,高德利,房軍.套管?水泥環(huán)?地層多層組合系統(tǒng)受力特性分析[J].應(yīng)用力學(xué)學(xué)報,2014,31(3):387?393.

    WANG Yanbin,GAO Deli,FANG Jun.Mechanical characteristics analysis of casing-cement ring-formation multilayer composite system[J].Chinese Journal of Applied Mechanics,2014,31(3):387?393.

    [9]徐芝綸.彈性力學(xué)(上冊)[M].北京:人民教育出版社,1978: 60?77.

    XU Zhilun.Elastic mechanics(volume1)[M].Beijing:People’s Education Press,1978:60?77.

    [10]楊桂通.彈塑性力學(xué)引論[M].北京:清華大學(xué)出版社,2004: 115?124.

    YANGGuitong.Induction to elasticity and plasticity[M].Beijing: Tsinghua University Press,2004:115?124.

    [11]李志明,殷有泉.油水井套管外擠力計算及其力學(xué)基礎(chǔ)[M].北京:石油工業(yè)出版社,2006:5?38.

    LIZhiming,YIN Youquan.Oil and water wells outside casing extrusion force calculation and mechanical foundation[M]. Beijing:Petroleum Industry Press,2006:5?38.

    [12]左孔天,錢勤,趙雨東,等.熱固耦合結(jié)構(gòu)的拓?fù)鋬?yōu)化設(shè)計研究[J].固體力學(xué)學(xué)報,2005,26(4):447?452.

    ZUO Kongtian,QIAN Qin,ZHAO Yudong,et al.Research on the topology optim ization about thermo-structural coup ling field[J].Chinese Journal of Solid M echanics,2005,26(4): 447?452.

    [13]戴鍋生.傳熱學(xué)[M].北京:高等教育出版社,1999:31?35.

    DAI Guosheng.Heat transfer[M].Beijing:Higher Education Press,1999:31?35.

    [14]陳艷華,張振迎,景寶國.試算法在多維非穩(wěn)態(tài)導(dǎo)熱問題中的應(yīng)用[J].建筑節(jié)能,2010,227(38):68?71.

    CHEN Yanhua,ZHANG Zhenying,JING Baoguo.Trial algorithm applied on multi-dimensional transient heat conduction[J].Building Energy Efficiency,2010,227(38): 68?71.

    (編輯 趙俊)

    M echanical characteristicsanalysisof casing?cement sheath?formationmultilayer com posite system w ith thermo-structural coup ling effects

    ZHAO Xinbo1,2,HAN Shengchao1,YANG Xiujuan1,YAN Xiangzhen1,3,ZHANG Lisong1

    (1.College of Pipelineand Civil Engineering,China University of Petroleum(EastChina),Qingdao 266580,China; 2.College of Science,Qingdao Technology University,Qingdao 266033,China; 3.Oiland Gas CAE Technology Research Center,China University of Petroleum(EastChina),Qingdao 266580,China)

    According to elastic mechanics,analytical solution of the equations under the action of thermo-structural coupling wasobtained by considering contact,continuity and boundary conditions in themultilayer com posite system of casing?cement sheath?formation in HPHT production well.The multilayer composite system was subjected to non-uniform in-situ stress,uniform inner pressure and tem perature pressure.The influences of cement thickness,cement elastic modulus,cement Poisson’s ratio,formation temperature and layer number of casing?cement sheath on M ises stress distribution of the insidewall of innermost layer casingwere discussed.The results show thatall the above factors have great influences on the M ises stress distribution.Comparing monolayer and multilayer casing?cement sheath?formation system of casing?cement sheath?formation,the M ises stress distribution of the inside wall of innermost layer casing inmultilayer casing?cementsheath?formation system is distinctly decreased,and so themultilayer casing?cement sheath?formation system hasan importantengineering value to ensure the safety of casing strength during service.

    casing?cementsheath?formation;thermo-structural coupling;multilayer com positesystem;elasticmechanics

    TE329

    A

    1672?7207(2017)03?0837?07

    10.11817/j.issn.1672-7207.2017.03.035

    2016?03?06;

    2016?05?30

    國家自然科學(xué)基金資助項目(51274231,51374228);國家自然科學(xué)基金?石化聯(lián)合基金資助項目(U1262208);中國石油天然氣集團(tuán)公司重點實驗室課題資助項目(2014A-4214);青島市應(yīng)用基礎(chǔ)研究計劃項目(15-9-1-71-jch)(Projects(51274231,51374228) supported by the National Natural Science Foundation of China;Project(U1262208)supported by the National Science Foundation of China and Petrochem ical Foundation;Project(2014A-4214)supported by the Key Laboratory of China National Petroleum Corporation;Project (15-9-1-71-jch)supported by the Applied Basic Research Program of Qingdao City)

    趙新波,博士,講師,從事油氣田地下工程力學(xué)研究;E-mail:zxbups@163.com

    猜你喜歡
    內(nèi)層泊松比內(nèi)壁
    ◆ 裝飾板材
    ◆ 裝飾板材
    裝飾板材
    ◆ 裝飾板材
    具有負(fù)泊松比效應(yīng)的紗線研發(fā)
    垣曲北白鵝墓地出土的青銅匽姬甗(M3:10)
    文物季刊(2021年1期)2021-03-23 08:12:58
    負(fù)泊松比功能的結(jié)構(gòu)復(fù)合紡紗技術(shù)進(jìn)展
    毛紡科技(2020年6期)2021-01-06 03:41:48
    考慮粘彈性泊松比的固體推進(jìn)劑蠕變型本構(gòu)模型①
    固體推進(jìn)劑粘彈性泊松比應(yīng)變率-溫度等效關(guān)系
    膀胱內(nèi)壁子宮內(nèi)膜異位1例
    亚洲狠狠婷婷综合久久图片| 久9热在线精品视频| 一卡2卡三卡四卡精品乱码亚洲| 免费一级毛片在线播放高清视频| 国产午夜福利久久久久久| eeuss影院久久| 国产精品女同一区二区软件 | 99热精品在线国产| 中文亚洲av片在线观看爽| 3wmmmm亚洲av在线观看| 丝袜美腿在线中文| 亚洲av不卡在线观看| 在线观看舔阴道视频| 91字幕亚洲| 欧美精品国产亚洲| av女优亚洲男人天堂| 欧美极品一区二区三区四区| 搞女人的毛片| 久久人人精品亚洲av| 成年女人永久免费观看视频| 在线免费观看不下载黄p国产 | 久久久久久九九精品二区国产| 美女黄网站色视频| 久久午夜亚洲精品久久| 国产精品电影一区二区三区| 动漫黄色视频在线观看| 亚洲18禁久久av| 日韩欧美精品免费久久 | 有码 亚洲区| 久久久色成人| 99热精品在线国产| 99国产极品粉嫩在线观看| 有码 亚洲区| 日韩有码中文字幕| 91在线观看av| 中文字幕人成人乱码亚洲影| 亚洲精华国产精华精| 亚洲av电影不卡..在线观看| av天堂中文字幕网| 国产三级中文精品| 亚洲熟妇熟女久久| 中文字幕精品亚洲无线码一区| 成人美女网站在线观看视频| 国内少妇人妻偷人精品xxx网站| www.999成人在线观看| 男人和女人高潮做爰伦理| 日本一本二区三区精品| 在线观看av片永久免费下载| 极品教师在线免费播放| 极品教师在线视频| .国产精品久久| 深夜a级毛片| 我要看日韩黄色一级片| 国产v大片淫在线免费观看| 亚洲欧美日韩东京热| 日韩中字成人| 久久亚洲真实| 日本一本二区三区精品| 欧美+亚洲+日韩+国产| 亚洲av日韩精品久久久久久密| 亚洲五月婷婷丁香| xxxwww97欧美| 天美传媒精品一区二区| 亚洲欧美日韩高清在线视频| 国产毛片a区久久久久| 久久精品人妻少妇| 欧美精品国产亚洲| 女人十人毛片免费观看3o分钟| 日韩有码中文字幕| 天堂网av新在线| 亚洲人成网站在线播放欧美日韩| 精品国产三级普通话版| 亚洲精品色激情综合| 久久天躁狠狠躁夜夜2o2o| 国产69精品久久久久777片| 国产av麻豆久久久久久久| 婷婷亚洲欧美| 亚洲欧美日韩无卡精品| 18禁在线播放成人免费| 99在线视频只有这里精品首页| 国产在视频线在精品| 亚洲,欧美,日韩| 国产精品一区二区性色av| 国产av在哪里看| 天堂√8在线中文| 嫁个100分男人电影在线观看| 亚洲一区高清亚洲精品| 午夜精品一区二区三区免费看| 两性午夜刺激爽爽歪歪视频在线观看| 国产av在哪里看| 亚洲av一区综合| 一级av片app| 国产精品乱码一区二三区的特点| 亚洲,欧美精品.| 禁无遮挡网站| 国产亚洲精品综合一区在线观看| 欧美bdsm另类| 国产亚洲精品久久久com| 久久精品国产99精品国产亚洲性色| 亚洲精品456在线播放app | 欧美精品国产亚洲| 高清毛片免费观看视频网站| 午夜激情福利司机影院| 国产av麻豆久久久久久久| 亚洲一区二区三区色噜噜| 男女下面进入的视频免费午夜| 少妇裸体淫交视频免费看高清| 色视频www国产| 国产精品不卡视频一区二区| 日韩欧美 国产精品| 日韩三级伦理在线观看| 18禁在线无遮挡免费观看视频| 免费黄色在线免费观看| 国产免费又黄又爽又色| 国产有黄有色有爽视频| 少妇丰满av| 国产精品女同一区二区软件| 又爽又黄无遮挡网站| 亚洲欧美成人精品一区二区| 中文天堂在线官网| 精品久久国产蜜桃| av卡一久久| 性色av一级| 久久久久久久久久成人| 亚洲精品乱码久久久久久按摩| 成人国产av品久久久| 18禁动态无遮挡网站| 日本一本二区三区精品| 中国美白少妇内射xxxbb| 看十八女毛片水多多多| 亚洲天堂av无毛| 一本色道久久久久久精品综合| 亚洲国产欧美人成| 亚洲国产色片| 欧美97在线视频| 大陆偷拍与自拍| 日本与韩国留学比较| 亚洲自偷自拍三级| 成人特级av手机在线观看| 性色avwww在线观看| 免费看av在线观看网站| 少妇的逼水好多| 亚洲第一区二区三区不卡| 精品国产三级普通话版| 深爱激情五月婷婷| 熟女人妻精品中文字幕| 国产 一区 欧美 日韩| 人妻制服诱惑在线中文字幕| 亚洲精品视频女| 少妇的逼好多水| 性色avwww在线观看| 亚洲久久久久久中文字幕| 看免费成人av毛片| 亚洲成人精品中文字幕电影| 国产精品一及| 五月开心婷婷网| 日本免费在线观看一区| 高清av免费在线| 毛片女人毛片| 大又大粗又爽又黄少妇毛片口| 欧美日韩综合久久久久久| 高清日韩中文字幕在线| 毛片一级片免费看久久久久| 国产一级毛片在线| 在线精品无人区一区二区三 | 久久久精品欧美日韩精品| 51国产日韩欧美| av在线老鸭窝| 女人久久www免费人成看片| 日本午夜av视频| 亚洲色图综合在线观看| 国产黄色视频一区二区在线观看| 99热这里只有是精品50| eeuss影院久久| 99久国产av精品国产电影| 午夜免费男女啪啪视频观看| 欧美bdsm另类| 一级黄片播放器| 色综合色国产| 久久精品综合一区二区三区| 久久精品国产鲁丝片午夜精品| av.在线天堂| 中国国产av一级| 亚洲人成网站在线播| 午夜福利在线观看免费完整高清在| 激情 狠狠 欧美| 国产在线一区二区三区精| 亚洲精品aⅴ在线观看| 久久久久久久精品精品| 午夜免费观看性视频| 久久精品久久久久久噜噜老黄| 亚洲最大成人手机在线| 高清视频免费观看一区二区| 五月玫瑰六月丁香| 国产久久久一区二区三区| 日韩欧美 国产精品| 国精品久久久久久国模美| 一级爰片在线观看| 3wmmmm亚洲av在线观看| 久久久久性生活片| 制服丝袜香蕉在线| 久久6这里有精品| 国产成人a区在线观看| 国产亚洲5aaaaa淫片| 美女内射精品一级片tv| 色5月婷婷丁香| 免费av毛片视频| 狂野欧美激情性bbbbbb| 天堂中文最新版在线下载 | 国产淫语在线视频| 日本wwww免费看| 乱系列少妇在线播放| 在线精品无人区一区二区三 | 久久久亚洲精品成人影院| 2018国产大陆天天弄谢| 亚洲精品一区蜜桃| 日韩欧美精品v在线| 日本一二三区视频观看| 岛国毛片在线播放| av在线天堂中文字幕| 国产日韩欧美亚洲二区| 最新中文字幕久久久久| 国产视频内射| 22中文网久久字幕| 性色avwww在线观看| 人妻一区二区av| 成人国产av品久久久| 日韩不卡一区二区三区视频在线| 久久久久性生活片| 我的老师免费观看完整版| 69人妻影院| 天堂网av新在线| 欧美日韩精品成人综合77777| 国产精品久久久久久精品电影小说 | 嘟嘟电影网在线观看| 成年女人看的毛片在线观看| 欧美成人一区二区免费高清观看| 日韩,欧美,国产一区二区三区| 一级片'在线观看视频| 亚洲国产成人一精品久久久| 国产精品国产三级专区第一集| 亚洲国产最新在线播放| 国产精品一区二区三区四区免费观看| 中国国产av一级| 一级爰片在线观看| 亚洲av.av天堂| 美女高潮的动态| 永久免费av网站大全| 久久韩国三级中文字幕| 搞女人的毛片| 听说在线观看完整版免费高清| 在线观看一区二区三区激情| 日韩中字成人| 成人二区视频| 亚洲av国产av综合av卡| 亚洲成人久久爱视频| 免费av不卡在线播放| 在线观看国产h片| 国产一区有黄有色的免费视频| 国产亚洲5aaaaa淫片| 国产老妇伦熟女老妇高清| 国产毛片a区久久久久| 亚洲最大成人av| 真实男女啪啪啪动态图| 亚洲国产av新网站| 中文字幕亚洲精品专区| 国产午夜精品久久久久久一区二区三区| 大香蕉97超碰在线| 国产毛片a区久久久久| 狂野欧美激情性xxxx在线观看| 亚洲av成人精品一区久久| 日韩中字成人| 香蕉精品网在线| 天天躁夜夜躁狠狠久久av| 国产免费一区二区三区四区乱码| 国产欧美亚洲国产| 少妇裸体淫交视频免费看高清| 日韩一区二区视频免费看| 国产成人91sexporn| 成人二区视频| 亚洲自拍偷在线| 久久久久精品性色| 插逼视频在线观看| 国产免费一区二区三区四区乱码| 国产白丝娇喘喷水9色精品| 免费大片18禁| 亚洲精华国产精华液的使用体验| 国产亚洲午夜精品一区二区久久 | 日本猛色少妇xxxxx猛交久久| 别揉我奶头 嗯啊视频| 国产 一区精品| 亚洲欧美一区二区三区国产| 两个人的视频大全免费| 2021少妇久久久久久久久久久| 黄色日韩在线| 色播亚洲综合网| 久久久久久久久久成人| 青青草视频在线视频观看| 可以在线观看毛片的网站| 99热网站在线观看| 日本三级黄在线观看| 亚洲天堂av无毛| 三级国产精品片| 日韩av不卡免费在线播放| 成人综合一区亚洲| 五月玫瑰六月丁香| 午夜精品一区二区三区免费看| 高清日韩中文字幕在线| videos熟女内射| 搡女人真爽免费视频火全软件| 亚洲人与动物交配视频| 男的添女的下面高潮视频| 2021少妇久久久久久久久久久| av天堂中文字幕网| 男男h啪啪无遮挡| 国产探花极品一区二区| 69人妻影院| 国产成人一区二区在线| 综合色丁香网| 久久精品国产亚洲av涩爱| 大片电影免费在线观看免费| 日韩视频在线欧美| 亚洲精品乱久久久久久| 最近最新中文字幕大全电影3| 亚洲怡红院男人天堂| 91精品国产九色| 国产精品国产三级专区第一集| 日韩 亚洲 欧美在线| 永久网站在线| 亚洲精品日本国产第一区| 精品熟女少妇av免费看| 好男人在线观看高清免费视频| 九九在线视频观看精品| 大片免费播放器 马上看| 亚洲欧美日韩东京热| 美女xxoo啪啪120秒动态图| 一级毛片aaaaaa免费看小| 亚州av有码| 97在线人人人人妻| 七月丁香在线播放| 久久久久久久亚洲中文字幕| 国产成年人精品一区二区| 尤物成人国产欧美一区二区三区| 99久久中文字幕三级久久日本| 久久精品国产a三级三级三级| 视频区图区小说| 欧美 日韩 精品 国产| 亚洲精品日韩av片在线观看| 在线亚洲精品国产二区图片欧美 | 联通29元200g的流量卡| 久久久国产一区二区| 午夜精品一区二区三区免费看| 亚洲精品色激情综合| 真实男女啪啪啪动态图| 中文精品一卡2卡3卡4更新| 极品少妇高潮喷水抽搐| 一区二区av电影网| 国产综合懂色| 久久人人爽人人爽人人片va| 可以在线观看毛片的网站| 国产成人福利小说| 亚洲av中文字字幕乱码综合| 久久影院123| 欧美高清成人免费视频www| 日韩强制内射视频| 久久久久久久久久成人| 69av精品久久久久久| 色综合色国产| 日韩欧美精品免费久久| 寂寞人妻少妇视频99o| 99久久精品国产国产毛片| 人体艺术视频欧美日本| 午夜福利视频1000在线观看| 制服丝袜香蕉在线| 国产美女午夜福利| freevideosex欧美| 亚洲精品中文字幕在线视频 | 久久鲁丝午夜福利片| 啦啦啦中文免费视频观看日本| 亚洲国产日韩一区二区| 91精品伊人久久大香线蕉| 女人久久www免费人成看片| 国产成人a∨麻豆精品| 成年版毛片免费区| 一区二区三区乱码不卡18| 精品国产一区二区三区久久久樱花 | 久久人人爽av亚洲精品天堂 | 亚洲成人中文字幕在线播放| 精品99又大又爽又粗少妇毛片| 少妇人妻 视频| 国产乱人偷精品视频| 亚洲伊人久久精品综合| 免费少妇av软件| 少妇的逼好多水| 成人午夜精彩视频在线观看| 中国国产av一级| 欧美人与善性xxx| 超碰av人人做人人爽久久| 国产成人freesex在线| 国产黄片视频在线免费观看| .国产精品久久| 日韩免费高清中文字幕av| 亚洲国产高清在线一区二区三| 一本一本综合久久| 国产永久视频网站| 啦啦啦啦在线视频资源| 男女下面进入的视频免费午夜| 国产爽快片一区二区三区| 亚洲美女视频黄频| 老女人水多毛片| 自拍欧美九色日韩亚洲蝌蚪91 | 国产真实伦视频高清在线观看| 日本wwww免费看| 国产精品99久久99久久久不卡 | 免费看日本二区| 人妻 亚洲 视频| 91久久精品电影网| 欧美一区二区亚洲| 中文在线观看免费www的网站| 亚洲成人中文字幕在线播放| 狂野欧美白嫩少妇大欣赏| 大又大粗又爽又黄少妇毛片口| 国产亚洲av片在线观看秒播厂| 人体艺术视频欧美日本| 国产 一区精品| 欧美成人a在线观看| 欧美潮喷喷水| 亚洲精品成人久久久久久| 简卡轻食公司| 亚洲第一区二区三区不卡| 在线观看免费高清a一片| a级一级毛片免费在线观看| 亚洲伊人久久精品综合| 热re99久久精品国产66热6| 午夜福利在线观看免费完整高清在| 在线观看三级黄色| 亚洲精品中文字幕在线视频 | 久久影院123| 久久这里有精品视频免费| 亚洲精品乱久久久久久| 蜜桃久久精品国产亚洲av| 精品人妻视频免费看| 久久人人爽人人片av| 中文天堂在线官网| 禁无遮挡网站| 色网站视频免费| 国产精品久久久久久av不卡| 亚洲国产日韩一区二区| 嫩草影院精品99| 青春草视频在线免费观看| 我的女老师完整版在线观看| 九九爱精品视频在线观看| 亚洲精品日韩av片在线观看| 天堂俺去俺来也www色官网| 久久久亚洲精品成人影院| 水蜜桃什么品种好| 日本色播在线视频| 国产亚洲av片在线观看秒播厂| 国产成人福利小说| 日韩国内少妇激情av| 97热精品久久久久久| 交换朋友夫妻互换小说| 国产69精品久久久久777片| 搡女人真爽免费视频火全软件| 乱码一卡2卡4卡精品| 高清av免费在线| 人人妻人人爽人人添夜夜欢视频 | av免费在线看不卡| av专区在线播放| 一级毛片久久久久久久久女| 国产老妇女一区| 水蜜桃什么品种好| 看非洲黑人一级黄片| 又黄又爽又刺激的免费视频.| 18禁在线无遮挡免费观看视频| 简卡轻食公司| 国产成人freesex在线| 男女边吃奶边做爰视频| 亚洲图色成人| 精品国产一区二区三区久久久樱花 | 国产毛片a区久久久久| 国产高清有码在线观看视频| 内射极品少妇av片p| 免费不卡的大黄色大毛片视频在线观看| 色视频www国产| 亚洲av日韩在线播放| 禁无遮挡网站| 日本-黄色视频高清免费观看| 久久人人爽人人片av| 视频区图区小说| 中文字幕久久专区| 午夜免费鲁丝| 99视频精品全部免费 在线| 最近2019中文字幕mv第一页| 日本午夜av视频| av在线老鸭窝| 我要看日韩黄色一级片| 网址你懂的国产日韩在线| 国产精品久久久久久精品电影小说 | 视频区图区小说| 国产淫语在线视频| 久久久久久久精品精品| 亚洲成色77777| 国产成人aa在线观看| 欧美bdsm另类| 欧美潮喷喷水| 日日摸夜夜添夜夜爱| 人妻夜夜爽99麻豆av| 亚洲精品456在线播放app| 一本色道久久久久久精品综合| 欧美变态另类bdsm刘玥| 在线观看免费高清a一片| 国产高潮美女av| 又爽又黄无遮挡网站| 一本—道久久a久久精品蜜桃钙片 精品乱码久久久久久99久播 | 亚洲美女视频黄频| 亚洲av中文字字幕乱码综合| 国产精品久久久久久av不卡| 男女啪啪激烈高潮av片| 久久精品国产a三级三级三级| 色哟哟·www| 精品酒店卫生间| 亚洲国产精品成人久久小说| 看非洲黑人一级黄片| 男女无遮挡免费网站观看| 观看美女的网站| 天堂中文最新版在线下载 | 男女啪啪激烈高潮av片| 亚洲图色成人| 搡老乐熟女国产| 国产精品一区二区在线观看99| 亚洲欧美中文字幕日韩二区| 日韩精品有码人妻一区| 成人午夜精彩视频在线观看| 18禁在线无遮挡免费观看视频| 亚洲激情五月婷婷啪啪| 久久精品国产鲁丝片午夜精品| 91久久精品国产一区二区三区| 五月开心婷婷网| 熟女人妻精品中文字幕| 国产伦精品一区二区三区视频9| 少妇的逼水好多| 久久韩国三级中文字幕| 亚洲精品亚洲一区二区| av黄色大香蕉| 搡老乐熟女国产| 在线天堂最新版资源| 女人被狂操c到高潮| 精品午夜福利在线看| 欧美人与善性xxx| 在现免费观看毛片| eeuss影院久久| 久久这里有精品视频免费| 男男h啪啪无遮挡| 男人舔奶头视频| 亚洲天堂国产精品一区在线| 亚洲av一区综合| 三级国产精品欧美在线观看| 久热久热在线精品观看| 国产成人精品婷婷| 成年版毛片免费区| 久久精品国产亚洲网站| 各种免费的搞黄视频| 国产一区二区三区综合在线观看 | 另类亚洲欧美激情| 亚洲综合精品二区| 深夜a级毛片| 男女无遮挡免费网站观看| 一个人观看的视频www高清免费观看| 国产爱豆传媒在线观看| 日本三级黄在线观看| 我要看日韩黄色一级片| 欧美xxⅹ黑人| 国产欧美日韩一区二区三区在线 | 网址你懂的国产日韩在线| 久久人人爽人人片av| 97超视频在线观看视频| 久久久久国产网址| 亚洲成人久久爱视频| 美女国产视频在线观看| 国产伦理片在线播放av一区| 在现免费观看毛片| 亚洲av成人精品一二三区| 久久亚洲国产成人精品v| 在现免费观看毛片| 成人国产麻豆网| 人人妻人人看人人澡| 高清欧美精品videossex| eeuss影院久久| 免费看不卡的av| 亚洲最大成人手机在线| 亚洲欧洲国产日韩| 国产伦理片在线播放av一区| 国产老妇伦熟女老妇高清| 日本wwww免费看| 国产伦在线观看视频一区| 在现免费观看毛片| 国产精品秋霞免费鲁丝片| 国产伦在线观看视频一区| 久久久久精品性色| 99热网站在线观看| 天天一区二区日本电影三级| 久久久久国产精品人妻一区二区| 一本色道久久久久久精品综合| 九草在线视频观看| 久久久久性生活片| 久久人人爽人人爽人人片va| av天堂中文字幕网| 美女视频免费永久观看网站| 激情五月婷婷亚洲| 国产精品久久久久久久电影| 亚洲精品国产av蜜桃| 99九九线精品视频在线观看视频| 日韩欧美 国产精品| 中文在线观看免费www的网站|