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

    產(chǎn)水氣井井筒溫度壓力計(jì)算方法

    2014-06-09 12:33:45蘇海洋位云生
    計(jì)算物理 2014年5期
    關(guān)鍵詞:氣液井口井筒

    李 波, 甯 波, 蘇海洋, 劉 虹, 位云生

    (中國石油勘探開發(fā)研究院,北京 100083)

    產(chǎn)水氣井井筒溫度壓力計(jì)算方法

    李 波, 甯 波, 蘇海洋, 劉 虹, 位云生

    (中國石油勘探開發(fā)研究院,北京 100083)

    應(yīng)用考慮氣液相間滑脫與流體沿井筒截面非均勻分布的漂移模型,結(jié)合動(dòng)量守恒、能量守恒和井筒傳熱學(xué),建立考慮井斜變化的氣井井筒溫度、壓力耦合預(yù)測(cè)模型,通過15口實(shí)測(cè)井?dāng)?shù)據(jù)對(duì)模型進(jìn)行驗(yàn)證.結(jié)果表明:所建模型可以準(zhǔn)確計(jì)算產(chǎn)水氣井井筒的溫度和壓力,計(jì)算結(jié)果平均絕對(duì)誤差3.60%,滿足工程計(jì)算要求.

    多相流;漂移模型;溫度;壓力;致密砂巖氣井

    0 引言

    地層流體(油、氣、水)從井底向井口流動(dòng)過程中往往伴隨著許多復(fù)雜的物理現(xiàn)象,流動(dòng)型態(tài)一般包括泡狀流、段塞流、過渡流和環(huán)霧流[1-7].井筒多相流動(dòng)常用的處理方法有三種:經(jīng)驗(yàn)?zāi)P?,均質(zhì)模型和機(jī)理模型.經(jīng)驗(yàn)?zāi)P透鶕?jù)多相流體運(yùn)動(dòng)的物理概念和物性參數(shù),利用因次分析法或根據(jù)基本運(yùn)動(dòng)方程得到描述某種特殊流動(dòng)的無因次參數(shù)建立多相流動(dòng)的經(jīng)驗(yàn)關(guān)系式.截止目前,國內(nèi)外已有許多經(jīng)驗(yàn)?zāi)P停?-7],但在實(shí)際運(yùn)用中,經(jīng)驗(yàn)?zāi)P托枰罅康挠?jì)算過程來判斷流型.此外,在某些流型之間的過渡區(qū)可能出現(xiàn)壓降和持液率的不連續(xù)現(xiàn)象,這樣就使得在求解井筒壓力和持液率等參數(shù)時(shí)出現(xiàn)不收斂或不連續(xù).均質(zhì)模型采用單相流的處理方法,認(rèn)為氣體和液體的流速相等,但這與實(shí)際情況不符,特別是在氣體流動(dòng)含量較低、流動(dòng)速度較慢的情況下[8-10].機(jī)理模型通常被認(rèn)為是最精確的模型,因?yàn)樗窃谠敿?xì)描述不同流型的物理現(xiàn)象的基礎(chǔ)上提出的,但是其方程系統(tǒng)復(fù)雜,求解困難,以及在某些流型的過渡區(qū)可能出現(xiàn)壓降和持液率不連續(xù)從而使計(jì)算不收斂,因此沒有被廣泛運(yùn)用.

    由于漂移模型具有相對(duì)簡(jiǎn)單、連續(xù)可微的優(yōu)點(diǎn),適用于井筒多相流動(dòng)建模.為此,對(duì)漂移模型開展了大量的研究[8-16],本文對(duì)漂移模型的重要參數(shù)(分布系數(shù)和漂移速度)進(jìn)行比較優(yōu)選,建立適合井筒多相流動(dòng)的漂移模型.通過漂移模型計(jì)算井筒任意深度的截面含氣率及流體物性參數(shù),再根據(jù)井筒傳熱學(xué)、動(dòng)量守恒和能量守恒,建立井筒溫度和壓力耦合模型,通過分段迭代法對(duì)模型進(jìn)行求解,并運(yùn)用15口產(chǎn)水氣井實(shí)測(cè)數(shù)據(jù)對(duì)模型進(jìn)行驗(yàn)證,計(jì)算一口實(shí)例水平井井筒壓力、溫度分布.

    1 漂移模型

    在井筒氣液兩相流動(dòng)過程中,由于氣相和液相的密度不一樣,存在氣相超越液相流動(dòng)的現(xiàn)象,即氣液兩相間存在滑脫,因此采用簡(jiǎn)單的均相模型來計(jì)算井筒的溫度和壓力分布會(huì)有一定的偏差.1965年,Zuber&Findlay[8]首次提出考慮了氣液相間的滑脫現(xiàn)象,以及氣相沿井筒徑向的變密度分布的漂移模型.而后,Bankoff[9],Wallis[10],Ishii[15-16],Hasan&Kabir[19-20],Ansari[5]等對(duì)漂移模型進(jìn)行了改進(jìn),使其廣泛用于核工業(yè)和石油工業(yè)中的多相流動(dòng)建模.由于漂移模型具有連續(xù)、可微和相對(duì)簡(jiǎn)單的特點(diǎn),Holmes等[13]于1997年將漂移模型引入數(shù)值模擬復(fù)雜結(jié)構(gòu)井井筒等溫多相流動(dòng)建模,但是由于之前的模型參數(shù)大多適合小直徑管流,不能直接用于井筒或地面設(shè)備中的多相流動(dòng)建模[14-15].為此,Oddie&Shi等[14-16]開展了大量的實(shí)驗(yàn),對(duì)模型參數(shù)進(jìn)行了回歸分析,并建立了適合大尺寸井筒多相流動(dòng)的漂移模型.本研究的目的是通過研究前人所建立的漂移模型,建立適合大尺寸井筒非等溫多相流動(dòng)的漂移模型,為準(zhǔn)確計(jì)算井筒溫度和壓力分布提供前提.

    1.1 漂移模型的基本關(guān)系式

    氣液兩相流動(dòng)情況下,氣相與液相間的滑脫現(xiàn)象是由于以下兩種機(jī)理結(jié)合產(chǎn)生的.一是氣體管橫截面上速度剖面和相分布剖面不一致,氣相在管中央處分布最高,管中央處混合物的流動(dòng)速度也最快,如圖1所示.因此,整個(gè)截面上氣相的平均速度要比液相的平均速度大.另外一種機(jī)理就是氣體受浮力的作用,傾向于垂直向上運(yùn)動(dòng).

    考慮以上兩種機(jī)理的模型,氣體的截面平均流速可表示為

    式中,vg為氣體的截面平均速度,m·s-1;Co為分布系數(shù)(描述兩相的分布特性,即流動(dòng)型態(tài)特性),無因次;vm為混合物的截面平均流速(無滑脫),m·s-1;vd為氣體的漂移速度(描述浮力效應(yīng)),m·s-1.

    混合物的截面平均速度是氣相和液相折算速度之和,即

    式中,vsg為氣體折算速度,m·s-1;vsl為液體折算速度,m·s-1;αg為截面平均含氣率;vl為液相的流動(dòng)速度,m·s-1.

    圖1 漂移模型示意圖Fig.1 Drift-flux model

    1.2 漂移模型的關(guān)鍵參數(shù)

    1.2.1 分布系數(shù)Co

    Zuber&Findlay給出Co的取值范圍是1~1.5[8].對(duì)于泡流和段塞流,Aziz、Ansari和Hasan&Kabir給出Co取值為1.2;而對(duì)于環(huán)霧流,Co取值接近于1.0[4-5].此外,在αg為1.0(高氣液比)時(shí),Co的值也應(yīng)為1.0.與此同時(shí),αgCo的值不能超過1.0.因此,Co表達(dá)式必須使其值在泡流和段塞流情況下有恒定不變的值,并且在αg為1.0時(shí)或混合物速度增大時(shí),Co取值減小為1.0,即

    為了使Co在αg和vm較大時(shí)取值減小到1.0,γ采用如下表達(dá)式

    在以上條件不能滿足時(shí),可能會(huì)出現(xiàn)不穩(wěn)定的情況.在B<(2-A)/A時(shí),以上兩個(gè)條件都成立.

    1.2.2 漂移速度vd

    漂移速度是氣液混合物流動(dòng)過程中氣液相對(duì)運(yùn)動(dòng)的截面平均速度值.在截面含氣率αg較小時(shí),氣體漂

    式中,vsgf為淹沒速度,即氣相折算流速到達(dá)足以使沿管壁液膜不會(huì)掉落下來時(shí)的速度值,流體的流動(dòng)將出現(xiàn)環(huán)狀流.參數(shù)A,B和Fv需要通過實(shí)驗(yàn)進(jìn)行調(diào)整,根據(jù)Oddie&Shi等[14-16]人的實(shí)驗(yàn)及模型參數(shù)優(yōu)化結(jié)果,氣液兩相流情況下,A取1.2,B取0.3,F(xiàn)v取1.0.

    淹沒速度用Wallis等[3]給出的關(guān)系式進(jìn)行求解

    式中,ku是臨界Kutateladze數(shù),具體算法請(qǐng)參考文獻(xiàn)[15];vc為特征流速,由式(8)進(jìn)行計(jì)算.

    顯然,氣相折算速度(vsg=αgvm)隨αg和vm的增加而增加.因此,Co必須滿足以下關(guān)系式移速度與氣體從靜止液體中向上移動(dòng)速度有關(guān),Harmathy[11]得到其值為1.53 vc,vc表達(dá)式為

    式中,σgl為氣-液相的表面張力,N·m-1;g為重力加速度,m·s-2;ρl為液相密度,kg·m-3;ρg為氣相密度,kg·m-3.

    在αg取較大值時(shí),可能出現(xiàn)氣液相沿不同方向流動(dòng)的逆流情況[14].在逆流的限制條件(vsg>0,vsl<0)的情況下,漂移速度是αg的函數(shù),即有

    對(duì)于處在以上兩種極端情況之間的過渡區(qū),則通過選用兩個(gè)臨界的截面含氣率a1和a2,然后插值進(jìn)行確定漂移速度.根據(jù)Oddie&Shi等[14-16]人的實(shí)驗(yàn)結(jié)果,a1取0.06,a2取0.21.

    綜合考慮以上情況,漂移速度的總表達(dá)式為

    式中,當(dāng)αg≤a1時(shí),K(αg)=1.53/Co;當(dāng)αg≥a2時(shí),K(αg)=Ku();在a1<αg<a2時(shí),通過線性插值得到其值.

    對(duì)于傾斜管流,鑒于Beggs&Brill對(duì)傾斜管流的處理方法,基于垂直管中的流動(dòng)情況,Hasan&Kabir[4]提出采用下式進(jìn)行修正

    其中,θ為傾斜角,(°);m(θ)為傾斜影響因子;vdθ是傾角為θ時(shí)的漂移速度,m·s-1;vd0為垂直流下的漂移速度,m·s-1.

    對(duì)于氣液兩相流動(dòng)系統(tǒng),Hasan-Kabir[4]提出m(θ)采用下面式子進(jìn)行計(jì)算

    根據(jù)Oddie&Shi等[14]在大管徑傾斜管流實(shí)驗(yàn),傾斜影響因子采用下式計(jì)算結(jié)果與實(shí)驗(yàn)結(jié)果吻合更好

    由以上研究可知,通過確定出漂移模型的兩個(gè)重要參數(shù)(分布系數(shù)和漂移速度),便可以得到井筒的截面處含氣率,進(jìn)而確定該井筒處流體的平均物性參數(shù).但從式(4)、(3)、(10)可知:分布系數(shù)、漂移速度和截面含氣率之間互為隱函數(shù),因此需要進(jìn)行迭代求解.

    2 井筒溫度壓力計(jì)算模型

    地層流體從井底流向井口過程中,不斷向地層中散發(fā)熱量,以及需要克服重力和摩擦力作用,井筒溫度、

    壓力不斷降低,而流體物性又與溫度、壓力之間相互影響,因此需要建立溫度、壓力和物性耦合模型,采用迭代法進(jìn)行求解.

    2.1 井筒壓力計(jì)算模型

    油氣井井筒壓力分布計(jì)算是油氣井生產(chǎn)設(shè)計(jì)、動(dòng)態(tài)分析和預(yù)測(cè)的核心問題.在此,假設(shè)流體在井筒中的流動(dòng)為一維穩(wěn)定流動(dòng),從管流系統(tǒng)任取一控制單元,建立如圖2所示的坐標(biāo)系.由動(dòng)量守恒定理有

    式中,ρm為混合物密度,kg·m-3;A為井筒截面,m2;D為井筒直徑,m;θ為井筒傾斜角,(°);vm為混合物流速,m·s-1;w為管壁摩擦力,N;dp為壓力變化量,Pa;dz為井筒長(zhǎng)度,m.

    整理得到壓力梯度

    圖2 井筒一維流動(dòng)示意圖Fig.2 One-dimensional flow in wellbores

    通過引入摩擦阻力系數(shù)fm,壓力梯度方程可表示為

    式中,ε為管壁粗糙度,mm;D為管徑,mm;Λ為無因次參數(shù),Λ=(ε/D)1.109 8/2.825 7+(7.149/Rem)0.898 1;Rem為無因次雷諾數(shù),Rem=Dυmρm/μm.

    混合物的密度ρm和混合物的粘度μm分別采用體積加權(quán)和質(zhì)量加權(quán):

    式中φ為截面體積含氣率(通過漂移模型計(jì)算),小數(shù);x為截面質(zhì)量含氣率,小數(shù);μg、μl為氣體、液體粘度,mPa·s.

    2.2 井筒溫度計(jì)算模型

    地層中的流體通過井筒從井底流到井口的過程中,地層溫度與井筒內(nèi)流體之間存在溫度差,井筒中流體與地層之間存在熱量交換,此外,流體自身的各種能量也會(huì)相互轉(zhuǎn)化.根據(jù)Ramey、Willhite、Hasan&Kabir、毛偉和廖新維等[17-23]人研究,井筒傳熱可作以下簡(jiǎn)化假設(shè):①流體在井筒中的流動(dòng)為一維穩(wěn)定流動(dòng);②氣液間不存在質(zhì)量交換;③井筒與地層之間只進(jìn)行熱量徑向傳遞,熱量徑向傳遞包括兩個(gè)過程:井筒到水泥外沿之間的穩(wěn)態(tài)傳熱,以及水泥環(huán)外向地層深處的非穩(wěn)態(tài)傳熱.

    由能量守恒定律有

    摩擦阻力系數(shù)fm可以采用Chen方法進(jìn)行計(jì)算,

    式中,h為流體的比焓,J;q為流體的徑向熱容量,J·s-1.

    由熱力學(xué)第一定律,得到流體比焓梯度

    式中,cP為流體的定壓比熱容,J·(kg·℃)-1;Tf為油管內(nèi)流體流動(dòng)溫度,℃;αJ為焦耳—湯姆遜系數(shù),℃·Pa-1.

    在單位長(zhǎng)度井段上,流體從油管到井壁的熱流量梯度為

    式中,Tf為井筒流體溫度,℃;Th為井壁溫度,℃;Uto為總傳熱系數(shù),W·(m·℃)-1;Wm為質(zhì)量流量,kg·s-1.

    根據(jù)Hasan&Kabir[17]提出的無因次時(shí)間函數(shù)f(tD),引入松弛因子A,由式(20)~(22)可得到井筒溫度梯度方程為

    式中,Te是地層溫度,℃;ke是地層導(dǎo)熱系數(shù),W·(m·℃)-1;rto是套管外徑,m.考慮焦耳—湯姆遜效應(yīng),穩(wěn)定流動(dòng)情況下(23)式的解析解為

    式中,Teout、Tein是計(jì)算段的出口和入口地層溫度,℃;zout、zin為計(jì)算段出口和入口深度,m;gT為計(jì)算井段處的地溫梯度,℃·m-1.

    2.3 井筒溫度、壓力計(jì)算

    由于井筒中流體物性、溫度和壓力之間關(guān)系緊密,相互影響,因此在計(jì)算井筒溫度和壓力時(shí),需要進(jìn)行耦合迭代求解.迭代求解面臨一個(gè)重要的問題是如何給定初值,因?yàn)槌跏贾禃?huì)影響迭代的收斂速度,本文采用(25)式計(jì)算井筒流出初始溫度值.

    溫度、壓力耦合迭代計(jì)算詳細(xì)步驟如下:

    1)給定井底條件:Tfin=TeBH,Tein=TeBH,pin=pBH,并選取合適的ΔL將井筒分為n段.

    2)根據(jù)(25)式計(jì)算粗略出口溫度Tfout.

    3)估計(jì)ΔL段壓力變化量Δp估計(jì),利用漂移模型計(jì)算各相流體的流動(dòng)速度和持率.

    5)根據(jù)(16)式計(jì)算ΔL段壓降梯度,然后計(jì)算出ΔL段的壓力變化量Δp計(jì)算.

    6)比較壓力變化估計(jì)量Δp估計(jì)和計(jì)算值Δp計(jì)算,若二者之差不滿足條件,則將值Δp計(jì)算作為初值,重復(fù)3)~6)步,直到滿足條件為止.

    7)計(jì)算井筒總傳熱系數(shù)和焦耳—湯姆遜系數(shù),利用(24)式精確計(jì)算段出口段的溫度值.

    8)比較第3)步與第7)步的溫度變化值,若二者之差滿足要求,則進(jìn)行下一段溫度壓力計(jì)算;否則,將第7)步計(jì)算的溫度變化值作為溫度初始值,跳到第2)步進(jìn)行,重新計(jì)算,直到滿足要求為止.

    9)計(jì)算該計(jì)算段對(duì)應(yīng)的深度Li、壓力pi和出口溫度Ti

    10)將Li處的壓力pi和Tf作為下一計(jì)算段的入口參數(shù),重復(fù)2)~9)步計(jì)算下一深度Li+1的壓力pi+1和溫度Ti+1,直到算完整個(gè)井筒溫度和壓力.

    3 模型驗(yàn)證與實(shí)例計(jì)算

    3.1 模型驗(yàn)證

    A氣田為典型致密砂巖氣藏,氣藏平均壓力系數(shù)為0.87,溫度110℃,天然氣相對(duì)密度0.598.氣藏油水關(guān)系復(fù)雜,氣井普遍產(chǎn)水,產(chǎn)出水包括凝析水和層間水.取A氣田15口典型氣井測(cè)試資料對(duì)所建模型進(jìn)行評(píng)價(jià),氣井產(chǎn)氣量為(3.98~14.59)×104m3·d-1,產(chǎn)水量為(1.11~18.34)m3·d-1,平均氣水比1.14× 104m3·m-3,井底流壓(21.54~28.62)MPa,測(cè)試井口壓(21.54~28.62)MPa,測(cè)試井口壓力(11.43~22.43)MPa,井的基本參數(shù)與實(shí)測(cè)如表1所示.油管直徑62 mm,管壁相對(duì)粗糙度0.004 5.

    由模型計(jì)算井口壓力與實(shí)測(cè)井口壓力對(duì)比可知(表2、表3、圖3),模型計(jì)算井口壓力與實(shí)測(cè)壓力誤差均小于10%,平均誤差-0.38%,平均誤差標(biāo)準(zhǔn)偏差為4.21%,平均絕對(duì)誤差為3.60%,平均絕對(duì)誤差標(biāo)準(zhǔn)差為2.22%,是參與評(píng)價(jià)模型中誤差最小的,其次是Hagedorn-Brown計(jì)算結(jié)果較好.而在油井管流計(jì)算中常用的Beggs-Brill和Orkiszewsiki模型,計(jì)算結(jié)果誤差較大,原因是這兩種模型都是通過實(shí)驗(yàn)建立的經(jīng)驗(yàn)關(guān)系式來計(jì)算井筒截面處的流體物性,受實(shí)驗(yàn)參數(shù)范圍(如:流量、氣液比等)和實(shí)驗(yàn)條件(如:管徑、管長(zhǎng)等)的限制,這些經(jīng)驗(yàn)關(guān)系式對(duì)低流速、較低氣液比的油井多相流計(jì)算具有較高精度,而對(duì)高流量、高氣液比氣井則適應(yīng)性較差,計(jì)算結(jié)果誤差較大.

    表1 氣井基本參數(shù)表Table 1 Basic data of gas wells

    表2 模型計(jì)算井口壓力與實(shí)測(cè)井口壓力對(duì)比Table 2 Comparison between calculated wellhead pressures and measured wellhead pressures

    通過所建模型計(jì)算還得到了井口溫度,但由于缺少現(xiàn)場(chǎng)井口溫度測(cè)試數(shù)據(jù),不能進(jìn)行溫度誤差分析.根據(jù)文獻(xiàn)[21]的觀點(diǎn),如果溫度和壓力是通過耦合計(jì)算得到的,只要壓力計(jì)算結(jié)果準(zhǔn)確,可以認(rèn)為溫度計(jì)算結(jié)果也可靠.此外,利用本文所建模型對(duì)文獻(xiàn)[23]中的實(shí)例氣井進(jìn)行計(jì)算,計(jì)算溫度、壓力值與實(shí)測(cè)溫度、壓力值誤差分別為1.54%和2.34%,說明本文計(jì)算方法是準(zhǔn)確的.因此,可以認(rèn)為本文建立的溫度壓力預(yù)測(cè)模型可靠,滿足工程計(jì)算要求.

    表3 計(jì)算結(jié)果誤差分析Table 3 Error analysis of results

    3.2 實(shí)例井計(jì)算

    長(zhǎng)深D平6井位于長(zhǎng)嶺斷陷中部隆起帶哈爾金構(gòu)造上,完鉆井深4 207 m,水平段長(zhǎng)1 007 m,水平段中深3 864 m,采用水平段裸眼封隔器+滑套完井.為滿足產(chǎn)氣與臨界攜液要求,采用73 mm油管進(jìn)行生產(chǎn).該井壓裂投產(chǎn)后日產(chǎn)氣19.53×104m3,水氣比1.77 m3·(104m3)-1,天然氣相對(duì)密度為0.59,井底流壓33.50 MPa,井底溫度133℃.

    由井筒溫度與壓力分布計(jì)算結(jié)果(圖4)可知,地層流體從井底向地面流動(dòng)過程中,由于克服重力、摩擦力和速度變化,井筒壓力不斷降低;井筒溫度不斷降低,但始終大于地層溫度,并且差值越來越大,這主要是由于流體從地層深部攜帶的熱量由于溫差需要不斷向地層中散發(fā)熱量造成的.通過計(jì)算:流體達(dá)到井口后溫度為65.68℃,油壓為20.61 MPa,計(jì)算油壓與井口實(shí)測(cè)油壓20.42 MPa,誤差僅為1.89%,模型計(jì)算精度較高.此外,在水平井的水平段中點(diǎn)至造斜段(3 120 m~3 864 m),井筒和溫度下降較慢;從井底到井口的總壓降為12.892 0 MPa,其中加速度壓降為0.001 9 MPa,摩擦壓降3.801 3 MPa,重力壓降9.088 8 MPa,可以看出井筒中的壓降主要是由摩擦力和重力引起的,而加速壓降則相對(duì)較小.

    圖3 模型計(jì)算壓力與實(shí)測(cè)井口壓力對(duì)比Fig.3 Comparison of calculated wellhead pressure with measured values

    4 結(jié)論

    1)漂移模型考慮了氣液相間的滑脫和流體沿井筒截面的非均勻分布,且具有相對(duì)簡(jiǎn)單、連續(xù)可微的特點(diǎn),可以準(zhǔn)確計(jì)算任意傾角情況下井筒截面處的流體分布,適合用于井筒多相流動(dòng)建模.

    2)基于漂移模型的井筒溫度、壓力耦合預(yù)測(cè)模型可以計(jì)算產(chǎn)水氣井井筒的溫度和壓力分布,以及流體的物性參數(shù),計(jì)算過程相對(duì)簡(jiǎn)單.對(duì)高氣液比產(chǎn)水氣井,所建模型計(jì)算精度高于常用的Hagedorn-Brown、Beggs-Brill和Orkiszewsiki模型,平均絕對(duì)誤差僅為3.60%,滿足工程計(jì)算要求.

    圖4 井筒溫度和壓力分布Fig.4 Temperature and pressure distribution in wellbores

    [1]Hagedorn A R,Brown K E.Experimental study of pressure gradients occurring during continuous two-phase flow in small diameter vertical conduits[J].JPT,1965,(4):475-484.

    [2]Orkiszewski J.Predicting two-phase pressure drops in vertical pipes[J].JPT,1967,23(6):829-838.

    [3]Beggs H D,Brill J P.A study of tow-phase flow in inclined pipes[J].JPT,1973,25(5):607-617.

    [4]Hasan A R,Kabir C S.Predicting multiphase flow behavior in a deviated well[J].SPE Prod Eng,1988,474(11):475-482.

    [5]Ansari A M,Sylvester N D,Sarica,et al.A comprehensive mechanistic model for upward two-phase flow in wellbores[J]. SPE Prod&Faci,1994,9(2):143-151.

    [6]Gomez L E,Ovadia S,Zelimir S,et al.Unified mechanistic model for steady-state two-phase flow:Horizontal to vertical upward flow[J].SPE Journal,2000,5(3):339-350.

    [7]Kaya A S,Sarica C,Bill J P.Mechanistic modeling of two phase flow in deviated wells[J].SPE Prod&Faci,2001,(9):156-166.

    [8]Zuber N,F(xiàn)indlay J.Average volumetric concentration in two-phase flow systems[J].Heat Transfer Trans ASME,1965,87:453-468.

    [9]Bankoff S G.A variable density,single fluid model for two-phase flow with particular reference to steam-water flow[J].Heat Transfer Trans ASME,1960,82:265-272.

    [10]Wallis G B,Makkenchery S.The hanging film phenomenon in vertical annular two-phase flow[J].Trans ASME,1974:186-197.

    [11]Harmathy T Z.Velocity of large drops and bubbles in media of infinite and restricted extent[J].AIChE J,1960,6(2):281-288.

    [12]Hibiki T,Ishii M.One-dimensional drift-flux model for two-phase flow in a large diameter pipe[J].International Journal of Heat and Mass Transfer,2003,46:1773-1790.

    [13]Holmes J A,Barkve T,Lund O.Application of a multi-segment well model to simulate flow in advanced wells[R].SPE 50646,1998.

    [14]Oddie G,Shi H,Durlofsky,et al.Experimental study of two and three phase flows in large diameter inclined pipes[J]. Multiphase Flow,2003,29:527-558.

    [15]Shi H,Holmes J A,Durlofsky,et al.Drift-flux modeling of two-phase flow in wellbores[J].SPEJ,2005,10(1):24-33. [16]Shi H,Holmes J A,Diaz L R,et al.Drift-flux parameters for three-phase steady-state flow in wellbores[J].SPEJ,2005,10 (2):130-137.

    [17]Ramey H J.Wellbore heat transmission[J].Trans,AIME,1962,14(4):427-435.

    [18]Willhite G P.Over-all heat transfer coefficients in stream and hot water injection wells[J].JPT,1967:607-615.

    [19]Hasan A R,Kabir C S.Aspects of heat transfer during two-phase flow in wellbores[J].SPE Prod&Faci,1994,9(3):211 -216.

    [20]Hasan A R,Kabir C S,Wang X.A robust steady-state model for flowing-fluid temperature model in complex wells[J].SPE Prod&Oper,2009,24(2),269-276.

    [21]Mao Wei,Liang Zheng.New method of calculating wellbore temperature in gas wells[J].Journal of Southwest Petroleum Institute,1999,21(1):56-66.

    [22]Mao Wei,Liang Zheng.Coupling analysis of the pressure and temperature in gas well borehole[J].Natural Gas Industry,1999,19(6):66-68.

    [23]Liao Xinwei,F(xiàn)eng Jilei.Pressure-temperature coupling calculation of transient wellbore heat transfer in deep geopressured gas reservoir[J].Petroleum Exploration and Development,2005,32(1):67-69.

    [24]Wu Xiaodong,Wu Han,Han Guoqing,et al.A new model for calculating wellbore temperature and pressure distribution of a high-H2S gas well considering the influence of the sulfur release in wellbores[J].Natural Gas Industry,2011,31(9):69-72.

    A Calculation Method for Temperature and Pressure Distribution in Water Production Gas Wellbores

    LI Bo,NING Bo,SU Haiyang,LIU Hong,WEI Yunsheng
    (Research Institute of Petroleum Exploration&Development,Beijing 100083,China)

    Based on drift model which copes with gas-liquid phase slippage and non-uniform distribution along section of wellbore,a coupled temperature-pressure prediction model was established for varied angle gas wellbores considering momentum conservation,conservation of energy and wellbore heat transmission.The model was validated with 15 tested gas wells.The calculation showed satisfying results with average absolute error of 3.6%which meets engineering requirements.

    multi-phase flow;drift model;temperature;pressure;tight sandstone gas wells

    date:2013-08-21;Revised date:2014-03-04

    TE319

    A

    2013-08-21;

    2014-03-04

    國家科技重大專項(xiàng)(2011ZX05015)資助項(xiàng)目

    李波(1985-),男,貴州思南,博士生,從事氣藏工程與多相管流研究,E-mail:libosonova@163.com

    1001-246X(2014)05-0573-08

    猜你喜歡
    氣液井口井筒
    140MPa井口壓裂四通管道沖蝕分析
    微重力下兩相控溫型儲(chǔ)液器內(nèi)氣液界面仿真分析
    氣液分離罐液位計(jì)接管泄漏分析
    CO2 驅(qū)低液量高氣液比井下氣錨模擬與優(yōu)化
    CO2驅(qū)油系統(tǒng)的井口預(yù)警系統(tǒng)設(shè)計(jì)
    礦井井筒煤柱開采技術(shù)措施
    煤峪口礦西三井筒提升中心的測(cè)定
    復(fù)雜地段副斜井井筒施工方法的選擇
    人間(2015年21期)2015-03-11 15:24:48
    氣田開發(fā)試采井口流程探討
    井底之蛙
    国产亚洲一区二区精品| 日本午夜av视频| 搡老乐熟女国产| 亚洲精品国产av成人精品| 啦啦啦啦在线视频资源| 一本一本久久a久久精品综合妖精| 成年美女黄网站色视频大全免费| 国产精品人妻久久久影院| av国产久精品久网站免费入址| 啦啦啦啦在线视频资源| 欧美激情极品国产一区二区三区| 久久精品久久久久久噜噜老黄| 青青草视频在线视频观看| 九色亚洲精品在线播放| 亚洲视频免费观看视频| 欧美成狂野欧美在线观看| 亚洲精品国产一区二区精华液| 最近手机中文字幕大全| 日韩视频在线欧美| 电影成人av| 在线av久久热| 中文乱码字字幕精品一区二区三区| 日本91视频免费播放| 国产三级黄色录像| 欧美人与性动交α欧美精品济南到| 亚洲国产精品一区三区| 免费在线观看影片大全网站 | 在线av久久热| 国产片特级美女逼逼视频| 久久久国产欧美日韩av| 男女免费视频国产| 久久久久视频综合| 亚洲欧美中文字幕日韩二区| 久久毛片免费看一区二区三区| 欧美在线一区亚洲| 亚洲精品国产色婷婷电影| 水蜜桃什么品种好| 亚洲国产精品国产精品| 久久精品aⅴ一区二区三区四区| 国产一卡二卡三卡精品| 美女高潮到喷水免费观看| 肉色欧美久久久久久久蜜桃| a级毛片在线看网站| 国产精品一国产av| 中文字幕最新亚洲高清| 超碰97精品在线观看| 国产午夜精品一二区理论片| 成人免费观看视频高清| 久久99一区二区三区| 大香蕉久久成人网| 日日爽夜夜爽网站| 亚洲精品一区蜜桃| 亚洲欧美一区二区三区黑人| 成年女人毛片免费观看观看9 | 日本黄色日本黄色录像| 久久这里只有精品19| 国产成人免费无遮挡视频| 老汉色∧v一级毛片| 在线精品无人区一区二区三| 一区二区av电影网| 老司机深夜福利视频在线观看 | 成年美女黄网站色视频大全免费| 免费久久久久久久精品成人欧美视频| 亚洲中文字幕日韩| 岛国毛片在线播放| 首页视频小说图片口味搜索 | 国产成人影院久久av| 国产片内射在线| 亚洲精品第二区| 亚洲人成网站在线观看播放| 国产精品成人在线| 99国产精品99久久久久| 丝袜喷水一区| 1024视频免费在线观看| 午夜免费鲁丝| 老熟女久久久| 欧美日韩亚洲综合一区二区三区_| 免费不卡黄色视频| 无限看片的www在线观看| tube8黄色片| 波野结衣二区三区在线| 一本色道久久久久久精品综合| 国产在视频线精品| 搡老岳熟女国产| 欧美日本中文国产一区发布| 欧美精品一区二区大全| av网站在线播放免费| 手机成人av网站| av欧美777| 99热全是精品| 国产在线视频一区二区| 性少妇av在线| 下体分泌物呈黄色| 欧美黄色片欧美黄色片| 精品国产一区二区久久| 欧美精品啪啪一区二区三区 | 老司机在亚洲福利影院| 精品人妻一区二区三区麻豆| 天天操日日干夜夜撸| 国产精品三级大全| 黄色 视频免费看| www.av在线官网国产| 黄色片一级片一级黄色片| 亚洲av片天天在线观看| 午夜影院在线不卡| 免费高清在线观看日韩| 亚洲一卡2卡3卡4卡5卡精品中文| 国产精品 欧美亚洲| 中文字幕高清在线视频| 国产视频一区二区在线看| 精品人妻一区二区三区麻豆| 亚洲免费av在线视频| 久久 成人 亚洲| 精品亚洲乱码少妇综合久久| 日韩一本色道免费dvd| 男女国产视频网站| 成人国产一区最新在线观看 | 欧美 亚洲 国产 日韩一| 乱人伦中国视频| 久久av网站| 男的添女的下面高潮视频| 国产人伦9x9x在线观看| 亚洲专区中文字幕在线| 国产成人精品无人区| 国产片内射在线| 亚洲精品中文字幕在线视频| 国产有黄有色有爽视频| 国产精品一区二区在线不卡| 十八禁高潮呻吟视频| 天堂中文最新版在线下载| 欧美国产精品一级二级三级| 91字幕亚洲| 欧美激情高清一区二区三区| av天堂久久9| 老司机靠b影院| 久久女婷五月综合色啪小说| 男女高潮啪啪啪动态图| 老熟女久久久| 视频在线观看一区二区三区| 久久精品久久精品一区二区三区| 国产日韩欧美在线精品| 国产一区二区三区av在线| 丁香六月欧美| 女人被躁到高潮嗷嗷叫费观| 1024香蕉在线观看| 亚洲av日韩在线播放| 亚洲欧洲国产日韩| 中文字幕精品免费在线观看视频| 乱人伦中国视频| 黄色视频不卡| 久久免费观看电影| 亚洲av男天堂| 午夜免费观看性视频| 性高湖久久久久久久久免费观看| 亚洲精品国产一区二区精华液| 亚洲成人免费电影在线观看 | 最近手机中文字幕大全| 亚洲,欧美,日韩| 亚洲国产av新网站| 侵犯人妻中文字幕一二三四区| www.999成人在线观看| 精品国产一区二区久久| 亚洲成人免费电影在线观看 | a 毛片基地| 欧美黑人欧美精品刺激| 成人国语在线视频| 国产一级毛片在线| 欧美在线一区亚洲| av在线播放精品| 国产在线观看jvid| 别揉我奶头~嗯~啊~动态视频 | 无遮挡黄片免费观看| 母亲3免费完整高清在线观看| 久久综合国产亚洲精品| 亚洲国产毛片av蜜桃av| 欧美日韩一级在线毛片| 亚洲精品一卡2卡三卡4卡5卡 | 黄频高清免费视频| 两个人看的免费小视频| 日韩视频在线欧美| 99精品久久久久人妻精品| 高清欧美精品videossex| kizo精华| 精品久久久久久电影网| 亚洲五月色婷婷综合| www.av在线官网国产| 婷婷丁香在线五月| 伊人久久大香线蕉亚洲五| 男女下面插进去视频免费观看| 青青草视频在线视频观看| 亚洲少妇的诱惑av| 麻豆av在线久日| 真人做人爱边吃奶动态| 国产不卡av网站在线观看| 亚洲av成人精品一二三区| 黄色视频在线播放观看不卡| 国产精品一区二区在线观看99| 亚洲国产中文字幕在线视频| 精品亚洲成a人片在线观看| 婷婷色麻豆天堂久久| 午夜日韩欧美国产| 欧美激情 高清一区二区三区| 性少妇av在线| 晚上一个人看的免费电影| 亚洲国产精品一区二区三区在线| 国产免费视频播放在线视频| 赤兔流量卡办理| 午夜老司机福利片| 在线观看www视频免费| 大陆偷拍与自拍| 成人18禁高潮啪啪吃奶动态图| 精品亚洲成a人片在线观看| 亚洲人成电影观看| 成年人午夜在线观看视频| 十八禁网站网址无遮挡| 亚洲欧洲日产国产| av在线老鸭窝| 少妇人妻 视频| 黄色视频在线播放观看不卡| 久久精品人人爽人人爽视色| 免费观看人在逋| 国产精品三级大全| 波多野结衣av一区二区av| 只有这里有精品99| 国产99久久九九免费精品| 久久久国产欧美日韩av| 国产国语露脸激情在线看| 丝袜脚勾引网站| 人人妻人人爽人人添夜夜欢视频| 大码成人一级视频| 亚洲伊人色综图| 中国美女看黄片| 亚洲国产欧美网| 久久亚洲国产成人精品v| 欧美黄色淫秽网站| 999久久久国产精品视频| 亚洲黑人精品在线| 精品少妇内射三级| 亚洲第一av免费看| 久久久亚洲精品成人影院| www.熟女人妻精品国产| 日韩中文字幕欧美一区二区 | 欧美大码av| av天堂在线播放| 嫁个100分男人电影在线观看 | 十八禁人妻一区二区| 欧美乱码精品一区二区三区| 黄色一级大片看看| 国产高清国产精品国产三级| 国产一区有黄有色的免费视频| 日本午夜av视频| 国产精品一区二区精品视频观看| 免费观看人在逋| 一级毛片女人18水好多 | 看免费av毛片| 另类亚洲欧美激情| 亚洲国产欧美在线一区| 亚洲少妇的诱惑av| 国产成人av教育| 一级黄色大片毛片| 伊人亚洲综合成人网| 亚洲国产欧美网| 亚洲av美国av| 久久久久精品国产欧美久久久 | 久久国产精品男人的天堂亚洲| 中国美女看黄片| 日日爽夜夜爽网站| av在线app专区| 男女午夜视频在线观看| 亚洲欧美精品自产自拍| 制服人妻中文乱码| 晚上一个人看的免费电影| 国产av精品麻豆| 国产精品亚洲av一区麻豆| 欧美乱码精品一区二区三区| 国产精品久久久av美女十八| 永久免费av网站大全| 丰满人妻熟妇乱又伦精品不卡| 欧美日韩黄片免| 男女下面插进去视频免费观看| 午夜日韩欧美国产| 视频区欧美日本亚洲| 久久人人爽av亚洲精品天堂| 亚洲视频免费观看视频| 色综合欧美亚洲国产小说| 亚洲av在线观看美女高潮| 国产一区二区三区av在线| 国产黄色视频一区二区在线观看| 精品福利永久在线观看| 亚洲精品日本国产第一区| 国产精品一区二区在线观看99| 久久久久久久久免费视频了| 宅男免费午夜| 自拍欧美九色日韩亚洲蝌蚪91| 男人舔女人的私密视频| 亚洲人成电影免费在线| 少妇粗大呻吟视频| 久9热在线精品视频| 亚洲欧美日韩另类电影网站| 久久久久久久久久久久大奶| 色婷婷久久久亚洲欧美| 免费女性裸体啪啪无遮挡网站| 精品高清国产在线一区| 性高湖久久久久久久久免费观看| 高清不卡的av网站| 在线天堂中文资源库| 两个人看的免费小视频| 国产成人一区二区在线| 久久久久精品国产欧美久久久 | 国产精品国产av在线观看| 欧美黄色淫秽网站| 精品亚洲成国产av| 99久久精品国产亚洲精品| 狠狠精品人妻久久久久久综合| 国产精品国产三级专区第一集| 男男h啪啪无遮挡| 亚洲精品国产av蜜桃| 午夜精品国产一区二区电影| 一二三四社区在线视频社区8| 亚洲激情五月婷婷啪啪| 七月丁香在线播放| 国产精品国产三级国产专区5o| 91精品三级在线观看| 香蕉国产在线看| 久久精品国产综合久久久| 超碰成人久久| 赤兔流量卡办理| 免费在线观看黄色视频的| 最近手机中文字幕大全| 国产日韩欧美亚洲二区| 久久久久久久久久久久大奶| 少妇的丰满在线观看| 不卡av一区二区三区| 国产欧美日韩一区二区三 | 在线观看免费视频网站a站| 麻豆国产av国片精品| 不卡av一区二区三区| www.av在线官网国产| 国产男人的电影天堂91| 欧美日韩av久久| 老汉色av国产亚洲站长工具| 高清黄色对白视频在线免费看| 桃花免费在线播放| 国产成人免费无遮挡视频| 欧美黑人精品巨大| 免费在线观看黄色视频的| 久久人人爽人人片av| 亚洲成人国产一区在线观看 | 婷婷成人精品国产| 国产一区二区激情短视频 | 欧美乱码精品一区二区三区| 另类亚洲欧美激情| 亚洲一区二区三区欧美精品| 亚洲欧美一区二区三区久久| 国产日韩欧美亚洲二区| 如日韩欧美国产精品一区二区三区| 两个人看的免费小视频| 亚洲精品国产区一区二| 国产成人精品久久二区二区91| 男男h啪啪无遮挡| 日韩 欧美 亚洲 中文字幕| 亚洲欧美成人综合另类久久久| 高清不卡的av网站| 国产欧美日韩一区二区三 | 欧美精品人与动牲交sv欧美| 十分钟在线观看高清视频www| 国产日韩欧美视频二区| 人体艺术视频欧美日本| 亚洲情色 制服丝袜| 每晚都被弄得嗷嗷叫到高潮| 黑丝袜美女国产一区| 欧美日韩视频高清一区二区三区二| av在线app专区| 婷婷色麻豆天堂久久| 一本综合久久免费| 欧美变态另类bdsm刘玥| 青草久久国产| 老司机在亚洲福利影院| 日韩,欧美,国产一区二区三区| av网站在线播放免费| 国产精品一区二区精品视频观看| 欧美日韩综合久久久久久| 成人18禁高潮啪啪吃奶动态图| 国产精品麻豆人妻色哟哟久久| 最新的欧美精品一区二区| 久久精品国产综合久久久| 国产在线观看jvid| 国产黄色免费在线视频| 黄色a级毛片大全视频| 久热这里只有精品99| 亚洲成国产人片在线观看| 亚洲精品国产av蜜桃| 久久人妻福利社区极品人妻图片 | 熟女av电影| 久9热在线精品视频| 男人爽女人下面视频在线观看| 在线精品无人区一区二区三| 久久久精品国产亚洲av高清涩受| 人妻人人澡人人爽人人| 亚洲第一青青草原| 大片电影免费在线观看免费| 国产精品久久久久久精品电影小说| 深夜精品福利| 久久久久国产一级毛片高清牌| 女人被躁到高潮嗷嗷叫费观| 婷婷色麻豆天堂久久| 国产免费一区二区三区四区乱码| www.自偷自拍.com| 久久久国产精品麻豆| 亚洲免费av在线视频| 侵犯人妻中文字幕一二三四区| 免费看不卡的av| 一区二区三区精品91| 久久精品久久精品一区二区三区| 一个人免费看片子| 高清黄色对白视频在线免费看| 91精品伊人久久大香线蕉| 考比视频在线观看| 日韩大码丰满熟妇| 高清欧美精品videossex| 国产欧美日韩精品亚洲av| 菩萨蛮人人尽说江南好唐韦庄| 91九色精品人成在线观看| kizo精华| 国产一区二区三区av在线| 久久国产亚洲av麻豆专区| 国产精品久久久av美女十八| 精品国产乱码久久久久久小说| 美国免费a级毛片| 亚洲精品美女久久久久99蜜臀 | 99久久人妻综合| 91国产中文字幕| 热99国产精品久久久久久7| 欧美激情 高清一区二区三区| 国产一区亚洲一区在线观看| 国产精品免费大片| 97人妻天天添夜夜摸| 国产高清不卡午夜福利| 欧美久久黑人一区二区| 男人添女人高潮全过程视频| 青草久久国产| 午夜免费鲁丝| 午夜福利乱码中文字幕| 日本vs欧美在线观看视频| 一级黄色大片毛片| 波野结衣二区三区在线| 91麻豆av在线| 伊人久久大香线蕉亚洲五| 只有这里有精品99| 一边摸一边抽搐一进一出视频| 免费在线观看完整版高清| 自拍欧美九色日韩亚洲蝌蚪91| 天堂俺去俺来也www色官网| 国产伦人伦偷精品视频| 亚洲av美国av| 极品人妻少妇av视频| 五月天丁香电影| 国产主播在线观看一区二区 | 一边摸一边做爽爽视频免费| 精品少妇黑人巨大在线播放| 国产成人免费观看mmmm| 自线自在国产av| 少妇裸体淫交视频免费看高清 | 爱豆传媒免费全集在线观看| 婷婷成人精品国产| 亚洲欧美一区二区三区黑人| 中国国产av一级| 精品亚洲乱码少妇综合久久| 国产精品久久久久久人妻精品电影 | 日韩中文字幕欧美一区二区 | 热re99久久国产66热| 高清欧美精品videossex| 亚洲人成网站在线观看播放| 精品久久蜜臀av无| 悠悠久久av| 亚洲一区二区三区欧美精品| 一本久久精品| 蜜桃在线观看..| 狂野欧美激情性xxxx| 一边摸一边做爽爽视频免费| 狠狠婷婷综合久久久久久88av| 看免费av毛片| 18禁国产床啪视频网站| 国产精品免费大片| 欧美av亚洲av综合av国产av| 国产欧美日韩综合在线一区二区| 又紧又爽又黄一区二区| 免费观看av网站的网址| 一级a爱视频在线免费观看| 一区在线观看完整版| 日本91视频免费播放| 亚洲精品国产色婷婷电影| 久久鲁丝午夜福利片| 人妻人人澡人人爽人人| 亚洲一区二区三区欧美精品| 自线自在国产av| 欧美日韩亚洲国产一区二区在线观看 | 国产精品二区激情视频| 日韩免费高清中文字幕av| 日韩av在线免费看完整版不卡| 欧美成人精品欧美一级黄| 久久久久视频综合| 亚洲男人天堂网一区| 亚洲国产最新在线播放| 亚洲一区二区三区欧美精品| av线在线观看网站| svipshipincom国产片| 在线天堂中文资源库| 热99国产精品久久久久久7| 国产精品成人在线| 无遮挡黄片免费观看| 成人18禁高潮啪啪吃奶动态图| 波野结衣二区三区在线| 又紧又爽又黄一区二区| 日韩制服骚丝袜av| 热99久久久久精品小说推荐| 王馨瑶露胸无遮挡在线观看| 久久精品成人免费网站| 成人国产一区最新在线观看 | 欧美激情极品国产一区二区三区| 我要看黄色一级片免费的| 日韩大码丰满熟妇| 极品少妇高潮喷水抽搐| 在线观看国产h片| 成在线人永久免费视频| 亚洲专区国产一区二区| 久久这里只有精品19| av视频免费观看在线观看| 久久国产亚洲av麻豆专区| 欧美精品一区二区大全| 天天躁夜夜躁狠狠躁躁| av网站在线播放免费| 亚洲精品美女久久av网站| 欧美精品啪啪一区二区三区 | 少妇的丰满在线观看| 亚洲男人天堂网一区| 又紧又爽又黄一区二区| 老司机午夜十八禁免费视频| 午夜福利一区二区在线看| 国产成人a∨麻豆精品| 男女之事视频高清在线观看 | 天天躁夜夜躁狠狠久久av| 中文字幕制服av| 亚洲精品国产一区二区精华液| 大香蕉久久网| 欧美精品亚洲一区二区| 日韩制服骚丝袜av| 欧美老熟妇乱子伦牲交| 最近手机中文字幕大全| 后天国语完整版免费观看| 少妇精品久久久久久久| 亚洲欧美日韩高清在线视频 | 波野结衣二区三区在线| 秋霞在线观看毛片| 成人国产av品久久久| 晚上一个人看的免费电影| 又粗又硬又长又爽又黄的视频| 亚洲精品美女久久av网站| 亚洲国产精品成人久久小说| 制服诱惑二区| 巨乳人妻的诱惑在线观看| 亚洲精品久久午夜乱码| 交换朋友夫妻互换小说| 久久人人97超碰香蕉20202| √禁漫天堂资源中文www| 波野结衣二区三区在线| 国产精品一二三区在线看| 成人国语在线视频| 国产在线视频一区二区| 欧美精品亚洲一区二区| 99热网站在线观看| 欧美人与善性xxx| 成年人午夜在线观看视频| 精品国产乱码久久久久久小说| av视频免费观看在线观看| 亚洲国产最新在线播放| 如日韩欧美国产精品一区二区三区| 精品国产一区二区三区四区第35| 国产成人免费无遮挡视频| 视频在线观看一区二区三区| 免费黄频网站在线观看国产| 精品少妇久久久久久888优播| 纵有疾风起免费观看全集完整版| 天天躁夜夜躁狠狠躁躁| 国产高清国产精品国产三级| 精品免费久久久久久久清纯 | 国产男女内射视频| 午夜福利免费观看在线| 韩国高清视频一区二区三区| 中文精品一卡2卡3卡4更新| 性高湖久久久久久久久免费观看| 国产日韩欧美视频二区| 亚洲国产欧美在线一区| 久久精品人人爽人人爽视色| 国产亚洲一区二区精品| 2018国产大陆天天弄谢| 夫妻性生交免费视频一级片| 另类精品久久| 日本猛色少妇xxxxx猛交久久| 大陆偷拍与自拍| 啦啦啦啦在线视频资源| 日本猛色少妇xxxxx猛交久久| 色网站视频免费| 美女视频免费永久观看网站| 国产成人啪精品午夜网站| 成人亚洲欧美一区二区av| 性色av一级| 国产成人啪精品午夜网站| 丝袜美腿诱惑在线| 国产精品国产三级国产专区5o| 午夜激情久久久久久久| 精品久久久久久久毛片微露脸 | 免费黄频网站在线观看国产| 啦啦啦在线免费观看视频4| 最新的欧美精品一区二区|