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

    油氣鉆井巖屑沉降阻力計算模型

    2018-06-26 02:08:24孫曉峰張克博曲晶瑀屈俊波
    天然氣工業(yè) 2018年5期
    關鍵詞:末速度雷諾數(shù)巖屑

    孫曉峰 張克博 陳 燁 曲晶瑀 湯 捷 屈俊波

    東北石油大學石油工程學院

    0 引言

    顆粒在介質(zhì)中沉降的阻力研究廣泛應用于選礦、煤炭、化工及石油工程等領域,石油天然氣工業(yè)領域典型的應用包括鉆井工程中巖屑和壓裂液中支撐劑的沉降。在鉆井工程中經(jīng)常會發(fā)生停泵工況,此時,為避免巖屑由于重力作用沉降在造斜段處下滑埋井眼等安全事故的發(fā)生,需要準確計算巖屑顆粒的沉降速度、預測井底巖屑床厚度。因此,深入研究巖屑顆粒在液體中的自由沉降規(guī)律,對于高效、安全地鉆井至關重要。

    1687年,牛頓(Newton)給出了平板在流體中運動的阻力公式[1];1851年Stokes[2]認為在顆粒雷諾數(shù)較低的情況下,可以忽略壓差阻力,并推導出了球體顆粒在液體中運動的黏性阻力公式;1867年,Rittinger[3]根據(jù)Newton提出的平板在流體中的阻力公式,推導出球體在流體中運動的阻力公式,后人稱之為Newton-Rittinger公式,最初的Newton-Rittinger公式與實驗值誤差較大,后來在不斷實踐與應用中進行了修正[4-6]。錢寧和萬兆惠[7]及Peden和Luo[8]認為Stokes公式應用于雷諾數(shù)小于0.1的情況,但是多數(shù)文獻[4-6,9-12]認為Stokes公式在雷諾數(shù)小于1時使用;Newton-Rittinger公式應用于湍流區(qū),適用雷諾數(shù)范圍[4-5]在103~105;1980年,Kao和Hwang[13]用Allen公式計算顆粒沉降速度,該公式適用的雷諾數(shù)范圍[6,14]介于30~300;對于過渡區(qū)的始段與末段選礦工程[4-6,14]中也給出了相應的經(jīng)驗公式。

    2007年Shah等[11]研究了多位學者[15-21]關于顆粒在非牛頓流體中的沉降實驗數(shù)據(jù),認為顆粒沉降阻力公式中的系數(shù)是雷諾數(shù)與流性指數(shù)n的函數(shù),其計算模型適用雷諾數(shù)介于0.001~1 000,Peden和Luo[8]也得出過同樣的結(jié)論,其模型適用的雷諾數(shù)范圍小于200;Shahi和Kuru[22-23]在2015年和2016年利用PIS(Particle Image Shadowgraph)技術分別研究了顆粒的沉降規(guī)律,得出了半經(jīng)驗公式,并對Shah等[11]的模型進行了修正,有效地提高了預測精度;2013年,Terfous等[24]利用球體顆粒的阿基米德數(shù),求解顆粒雷諾數(shù)和阻力系數(shù),該方法適用于雷諾數(shù)介于(0.1~5)×104的球體顆粒沉降;2012年和2016年,Elgaddaf i等[25-26]以聚陰離子纖維素(PAC)懸浮液、低毒礦物油以及Herschel-Bulkley流體作為基礎液,并在其中添加不同濃度的纖維,發(fā)現(xiàn)在基礎液中添加纖維能為顆粒增加額外的纖維阻力,有效降低了顆粒的沉降速度;2010年和2014年Bro?ek和Surowiak等[27-28]考慮了顆粒雷諾數(shù)大于500時的情況,建立了沉降末速度預測模型,所得沉降預測模型與Newton-Rittinger公式很接近。

    雖然選礦工程中各流態(tài)都有相應的球形顆粒沉降末速度預測模型,但是在選擇計算模型之前必須判斷雷諾數(shù),而顆粒雷諾數(shù)又與沉降末速度相關,盡管使用圖表法[4-6]可以不用判斷雷諾數(shù),但是這會造成較大的誤差且使用起來不夠方便;而近些年對于沉降的研究多以球體或天然沙粒為實驗對象,討論顆粒在特定流體中的沉降規(guī)律,且不同的模型有不同的適用范圍。而在鉆井工程中,巖屑顆粒的形態(tài)與球體及天然沙粒有較大差別,因此,筆者以上返巖屑顆粒為實驗對象,具體研究其沉降規(guī)律,并試圖提出一種適用顆粒雷諾數(shù)范圍更廣的顆粒沉降阻力預測模型。

    1 巖屑在牛頓流體中的沉降阻力模型

    巖屑顆粒在沉降時受到的阻力除了與流體密度、沉降速度、顆粒尺寸及流體黏度有關外,還與巖屑的表面粗糙度和形狀系數(shù)相關。利用π定理求巖屑顆粒在液體中沉降時所受阻力,即:

    這里,我們選取顆粒當量直徑、流體密度和沉降速度作為基本變量,則利用π定理,有則此時顆粒所受阻力的函數(shù)關系可表達為:

    即:

    由式(2)可知,顆粒沉降時所受阻力是多個變量的函數(shù),需要進行一定的簡化。由于巖屑顆粒的形態(tài)千變?nèi)f化,為方便描述它的尺寸,一般將其用當量球體表示,則此時顆粒形狀系數(shù)為常量1;同時,顆粒雷諾數(shù)可以表示為 ;此外,由于巖屑顆粒形態(tài)的差異化,很難將表面粗糙度δ作為參變量來確定其對沉降阻力的影響,但是,通過量綱分析法可以確定自由沉降速度是粗糙度的函數(shù),而自由沉降速度又與顆粒雷諾數(shù)相關,因此,式(2)可以化簡為:

    式中Ff表示顆粒沉降所受阻力,N;ρf表示沉降液密度,kg/m3;μ表示沉降液黏度,Pa·s;dp表示顆粒當量直徑,m;v表示顆粒自由沉降末速度,m/s;δ表示絕對粗糙度,m;χ表示顆粒形狀系數(shù),無量綱;表示顆粒沉降阻力系數(shù),無量綱。

    顆粒在液體中沉降時主要受到重力、浮力以及阻力,當顆粒達到沉降末速度時其受力處于平衡狀態(tài),則此時有:

    令式(3)等于式(4),得顆粒沉降阻力系數(shù)表達式:

    式中V表示巖屑顆粒的體積,m3;ρs表示巖屑顆粒密度,kg/m3;g表示重力加速度,m/s2。

    當把巖屑視為當量球體時,并且在顆粒雷諾數(shù)小于1的情況下,顆粒沉降所受到的壓差阻力可以忽略不計,主要考慮黏性阻力,基于此,1851年Stokes[2]推導得出球形顆粒沉降時的黏性阻力式:

    將式(6)帶入式(4),并用球體顆粒的體積表達式代替巖屑體積,此時可得沉降末速度(vs)表達式為:

    隨著顆粒沉降速度的增大,黏性阻力所占比重逐漸減小,壓差阻力不再忽略,當達到湍流時(一般認為是顆粒雷諾數(shù)大于1 000),黏性阻力對沉降影響較小,主要考慮壓差阻力,此時用修正后的Newton-Rittinger式為:

    當在紊流繞流區(qū)間完全忽略黏性阻力時,壓差阻力等于浮重,此時式(8)中阻力系數(shù)取 ,將式(8)帶入式(4)中,湍流區(qū)顆粒沉降末速度(vN-R)表達式為:

    在過渡區(qū)(即顆粒雷諾數(shù)大于1小于1 000時),顆粒沉降既受到黏性阻力,又有壓差阻力,任何一項都不可忽略,當受力平衡時,黏性阻力與壓差阻力之和等于有效浮重,表達式為:

    隨著流態(tài)從層流到過渡區(qū),再到紊流區(qū),顆粒所受黏性阻力占比從最大逐漸減小到可以忽略,而壓差阻力的趨勢正好相反,但是總阻力依然等于浮重,得:

    我們稱A和B分別為黏性阻力占比系數(shù)與壓差阻力占比系數(shù)。實際上,隨著顆粒雷諾數(shù)的變化,黏性阻力及壓差阻力不能簡單地用Stokes公式或Newton-Rittinger公式乘以一個常系數(shù),由式(3)可知阻力系數(shù)為雷諾數(shù)的函數(shù),對比式(10),認為阻力占比系數(shù)是顆粒雷諾數(shù)的函數(shù),即A=g(Rep),B=1 - g(Rep)。

    2 巖屑在冪律流體中的沉降阻力模型

    現(xiàn)場除少數(shù)鉆井液體系為牛頓流體外,多數(shù)鉆井液體系都符合冪律流型。對于冪律流體來說,阻力系數(shù)除了與顆粒雷諾數(shù)相關外,還與流性指數(shù)相關[8,11,20],但當顆粒雷諾數(shù)小于10時,顆粒在冪律流體沉降的阻力系數(shù)接近于在牛頓流體中的值[8],此時,Stokes定律依然適用。

    對于冪律流體,其黏度表達式[16]為:

    式中K表示稠度系數(shù),Pa·sn;n表示流性指數(shù),無量綱。

    將上式帶入到式(7)中,則斯托克斯沉降速度式變?yōu)椋?/p>

    當沉降處于紊流區(qū)時,顆粒只受壓差阻力,阻力系數(shù)為常量,與雷諾數(shù)及流性指數(shù)無關,所以冪律流體中顆粒沉降的壓差阻力式不變,沉降末速度vN-R也不變,此時顆粒沉降阻力預測模型式(10)變?yōu)椋?/p>

    由于顆粒在牛頓流體中沉降時,阻力系數(shù)為雷諾數(shù)的函數(shù),認為阻力占比系數(shù)也為雷諾數(shù)的函數(shù),同樣地,由于冪律流體中顆粒沉降阻力系數(shù)是雷諾數(shù)和流性指數(shù)的函數(shù)[8,11,20],則可以認為冪律流體中的顆粒沉降阻力占比系數(shù)A和B也為雷諾數(shù)和流性指數(shù)的函數(shù),即冪律流體中,X、Y為無量綱參數(shù),是流性指數(shù)n的函數(shù)。

    3 巖屑顆粒在飽和鹽水與冪律流體中的沉降實驗結(jié)果及分析

    3.1 沉降實驗

    筆者實驗所用巖屑分別選取于鶴煤3井500~800 m井段的碳酸鹽巖和松遼盆地X 井400~700 m沉積巖井段,共計18組,用不同孔徑的標準篩充分篩選后,分別選取粒徑在0.112~3.000 mm與0.45~1.25 mm范圍內(nèi)的巖屑顆粒作為沉降實驗研究對象。實驗用透明有機玻璃管模擬井筒,其內(nèi)徑為100 mm,長度為2 m,沉降液為常溫下的飽和鹽水和不同濃度的聚丙烯酰胺(PAM)溶液(PAM溶液的流變參數(shù)如表1所示),巖屑顆粒密度范圍介于2.1~2.4 g/cm3。

    表1 PAM溶液流變參數(shù)表

    巖屑粒徑與圓筒內(nèi)徑的比值小于0.03,盡可能地減小了管壁附加的阻尼效應。通過測量巖屑經(jīng)過玻璃管上兩個標記處的時間差來計算顆粒的沉降末速度,初始計時點距離玻璃管頂端100 cm,而終了計時點距離頂端170 cm。為減小測量誤差,采用以下方法進行實驗:每組沉降實驗先測量10次,求這10次測量的平均值,再進行第11次測量,求得前11次實驗的平均值;如果這兩個平均值之間的誤差不超過1%,則結(jié)束本組測量,如果大于1%,繼續(xù)實驗,直至最后一次求得的值與前一次所求值誤差不超過1%時一組測量結(jié)束;更換實驗對象,繼續(xù)上述步驟,至實驗結(jié)束。

    3.2 顆粒沉降實驗結(jié)果及分析

    3.2.1 實驗結(jié)果

    根據(jù)第3節(jié)的顆粒沉降實驗,巖屑在飽和鹽水溶液沉降實驗共獲得了91組數(shù)據(jù),不同樣本與粒徑范圍的巖屑顆粒沉降末速度如表2所示。

    根據(jù)表2的實驗結(jié)果,以飽和鹽水各密度范圍巖屑顆粒沉降末速度的平均值為依據(jù),繪制出巖屑顆粒沉降末速度隨粒度變化的曲線(圖1)。

    由圖1可知,顆粒沉降末速度隨著巖屑密度與巖屑粒徑的增加而增加,但這種變化并不是線性的,而是呈對數(shù)曲線變化的趨勢(底數(shù)大于1的對數(shù)函數(shù)),即巖屑沉降末速度的增速隨著粒徑的增大而減緩,雖然巖屑顆粒的外觀呈現(xiàn)出不同形態(tài),但這種變化趨勢與球體顆粒在牛頓流體中的沉降規(guī)律是一致的[24]。

    表2 飽和鹽水溶液中巖屑沉降末速度實驗值表

    圖1 不同密度的飽和鹽水中巖屑沉降末速度變化曲線圖

    根據(jù)理論分析及實驗數(shù)據(jù),可以得出顆粒雷諾數(shù)Rep與黏性阻力占比系數(shù)A的函數(shù)關系式:

    再聯(lián)立式(11),可得顆粒雷諾數(shù)Rep與壓差阻力占比系數(shù)B的函數(shù)關系式:

    則此時將式(15a)和式(15b)帶入式(10),預測模型變?yōu)椋?/p>

    巖屑在PAM溶液沉降實驗共獲得了70組數(shù)據(jù),不同樣本與粒徑范圍的巖屑顆粒沉降末速度如表3所示。

    根據(jù)表3的實驗結(jié)果,以不同粒徑范圍沉降末速度的平均值為依據(jù),作巖屑顆粒在冪律流體中沉降末速度隨粒度的變化曲線(圖2)。

    由圖2可知,顆粒沉降末速度隨著粒徑增加而增加的規(guī)律與前文巖屑在牛頓流體中的沉降規(guī)律并不相同。巖屑在牛頓流體中沉降時隨粒徑增大呈對數(shù)曲線變化的趨勢(底數(shù)大于1的對數(shù)函數(shù)),而巖屑在冪律流體中沉降時隨粒徑增大呈指數(shù)曲線變化的趨勢(底數(shù)大于1的指數(shù)函數(shù)),即巖屑沉降末速度的增速隨著粒徑的增大而增大,雖然巖屑顆粒的外觀呈現(xiàn)出不同形態(tài),但這種變化趨勢與球體顆粒在冪律流體中的沉降規(guī)律是一致的[24]。為了研究顆粒在冪律流體中沉降時阻力占比系數(shù)A和B的變化規(guī)律,結(jié)合本文實驗數(shù)據(jù)及相關研究數(shù)據(jù),得到不同流性指數(shù)下無量綱參數(shù)X與Y的值如表4所示。

    根據(jù)表4的數(shù)據(jù),回歸得到無量綱參數(shù)X和Y關于流性指數(shù)n與稠度系數(shù)K乘積的函數(shù)關系:

    此時,式(14)變?yōu)椋?/p>

    3.2.2 顆粒沉降末速度求解步驟

    步驟(1):對于冪律流體,根據(jù)已知流變參數(shù),求出無量綱參數(shù)X、Y的值。

    表3 不同濃度PAM溶液中巖屑沉降末速度實驗值表

    圖2 不同濃度PAM溶液中巖屑沉降末速度變化曲線圖

    步驟(2):根據(jù)巖屑顆粒的已知參數(shù),利用式(7)、(9)可計算出牛頓流體中巖屑顆粒的vs和vN-R,而利用公式(13)、(9)則可求解冪律流體顆粒的和vN-R。

    步驟(3):對于牛頓流體,將步驟(2)所求得的vs和vN-R值帶入式(16),得到顆粒阻力表達式,對于冪律流體,將步驟(2)所求得的vs'和 vN-R帶入

    式(18),得到顆粒阻力表達式。

    步驟(4):根據(jù)步驟(3)并聯(lián)立式(4)可以求出巖屑的顆粒雷諾數(shù)Rep。

    根據(jù)步驟(4)求得的顆粒雷諾數(shù)及其定義,求得顆粒沉降末速度v。

    3.2.3 結(jié)果分析

    當顆粒處在斯托克斯區(qū)域時,可以忽略壓差阻力,此時A=1,B=0,由式(15a)可得此時Rep=1.11,即式(14)所適用的顆粒雷諾數(shù)范圍的下限為1.11。當顆粒雷諾數(shù)不斷增大時,有:

    由式(19)可知,當時,系數(shù)A=0,B=1,此時顆粒沉降可以只考慮壓差阻力,即預測模型適用顆粒雷諾數(shù)范圍的上限是正無窮。然而由于式(8)的顆粒雷諾數(shù)適用范圍介于103~105。因此,預測模型適用的顆粒雷諾數(shù)范圍介于1.11~105。

    在巖屑沉降時,隨著顆粒雷諾數(shù)的逐漸增大,巖屑的黏性阻力占比系數(shù)以較快的變化速度在不斷地減小,而壓差阻力占比系數(shù)的變化趨勢正好相反,當顆粒雷諾數(shù)達到2.557 2時,黏性與壓差阻力占比系數(shù)相等(圖3),顆粒雷諾數(shù)達到100時,變化趨勢趨于平緩,此后,壓差阻力占比系數(shù)無限趨向于1,而黏性阻力占比系數(shù)無限趨近于0。

    表4 不同流性指數(shù)下無量綱參數(shù)X與Y值表

    圖3 阻力占比系數(shù)隨顆粒雷諾數(shù)變化曲線圖

    假設一定粒徑與密度的顆粒在靜止流體中沉降,保持液體的密度近似不變,僅適度改變沉降液的黏度,此時顆粒所受浮力不變,即阻力一定,研究顆粒受到的黏性阻力與壓差阻力隨顆粒雷諾數(shù)變化的趨勢(圖4)。雷諾數(shù)較小時顆粒沉降只受黏性阻力,當顆粒雷諾數(shù)逐漸增加時,黏性阻力以較快的速率減小,雷諾數(shù)為100時,黏性阻力占比為12%,雷諾數(shù)為500時,占比為10%,可見在顆粒雷諾數(shù)在達到100后下降趨勢減緩,而壓差阻力的變化趨勢正好與之相反。同時,由圖4可知,在顆粒雷諾數(shù)達到2.944 6時,顆粒受到的黏性阻力與壓差阻力大小相等,但此時的雷諾數(shù)值較圖3中占比系數(shù)相等時的值大,這主要是因為在計算時式(14)中的相關值取較小數(shù)π/20,即在湍流區(qū)使用Newton-Rittinger公式時并沒有完全忽略黏性阻力。

    圖4 顆粒所受的黏性阻力和壓差阻力隨顆粒雷諾數(shù)變化曲線圖

    對預測模型的阻力系數(shù)進行誤差分析,可以得出誤差主要控制介于±20%的范圍,在試驗采集的91組數(shù)據(jù)中,有86組數(shù)據(jù)的阻力系數(shù)誤差在這一范圍內(nèi)。對預測模型沉降末速度進行誤差分析,發(fā)現(xiàn)有超過81%的誤差控制在±8%以內(nèi),其中,有54組數(shù)據(jù)的絕對誤差小于5%,33組數(shù)據(jù)絕對誤差介于5%~10%,3組數(shù)據(jù)絕對誤差介于10%~20%,只有1組數(shù)據(jù)的絕對誤差大于20%,達到22.56%,同時,該數(shù)據(jù)點所對應的阻力系數(shù)絕對誤差也較大,達到33.43%,分析數(shù)據(jù)發(fā)現(xiàn),該數(shù)據(jù)點的實驗值與同性質(zhì)巖屑顆粒的沉降末速度實驗值相差較大,認為造成該點數(shù)值誤差較大的原因極有可能是實驗誤差所致。最終,預測模型所得的沉降末速度與實驗值的絕對平均誤差為5.02%。

    對于顆粒在冪律流體中的沉降,由于X和Y是流性指數(shù)n的函數(shù),需要根據(jù)流體的流變參數(shù)確定X和Y值,再令黏性阻力占比系數(shù)A=1,可以確定出不同流變參數(shù)下式(18)的適用下限。同時,對某一特定的冪律流體而言,X可以求得且大于0,Y小于0,則A、B的變化趨勢與顆粒在牛頓流體中沉降的變化趨勢相同。當Rep不變,流性指數(shù)n與稠度系數(shù)K的乘積變化時,A的變化規(guī)律如圖5所示。由圖5可知,不同顆粒雷諾數(shù)下黏性阻力占比系數(shù)A隨流性指數(shù)n與稠度系數(shù)K的乘積增大而減小。

    圖5 黏性阻力占比系數(shù)A隨n與K乘積的變化曲線圖

    因為式(17)中無量綱參數(shù)X與Y的求解是依據(jù)有限的數(shù)據(jù)擬合得到的,所以預測模型式(18)僅適用有限的流性指數(shù)范圍,根據(jù)表4的數(shù)據(jù),可以確定式(18)適用于流性指數(shù)介于0.062 3到1之間的冪律流體,適用的顆粒雷諾數(shù)下限由冪律流體的流性指數(shù)n與稠度系數(shù)K的乘積決定,適用上限不受流變參數(shù)的影響,為Newton-Rittinger公式的適用范圍,最大可達到105。

    當巖屑顆粒在牛頓流體中沉降時,使用式(16)求解出顆粒雷諾數(shù),若顆粒雷諾數(shù)介于1.11~105時,根據(jù)顆粒雷諾數(shù)的定義由求解出的雷諾數(shù)值計算沉降末速度;若由式(16)所得的雷諾數(shù)小于1.11,則此時得出的顆粒雷諾數(shù)值與顆粒沉降的真實值有較大差距,不能由顆粒雷諾數(shù)的定義直接反求顆粒沉降末速度,但此時已經(jīng)判斷出雷諾數(shù)小于1.11,因而可以用Stokes定律,即式(7)直接求解沉降末速度。因此,根據(jù)筆者提出的牛頓阻力沉降公式,可以計算顆粒雷諾數(shù)介于0~105的巖屑在牛頓流體中的沉降。

    同理,求解顆粒在冪律流體中的沉降末速度時,先由冪律流體的流變參數(shù)求得預測模型式(18)的適用下限,再由式(18)求得顆粒雷諾數(shù),若其值在預測模型的適用范圍內(nèi),可直接求解顆粒沉降末速度值;若所得值小于預測模型的適用下限,此時選用斯托克斯沉降公式的變形式(13)求解顆粒沉降末速度值。因此,根據(jù)本文提出的顆粒在冪律流體中的沉降模型,可以計算顆粒雷諾數(shù)介于0~105,流性指數(shù)介于0.062 3~1的顆粒沉降。

    由于文中將巖屑顆粒視為球體,使用體積當量直徑所得顆粒當量直徑(dp)結(jié)果較為準確,求解式[4]為:

    式中d1、d2表示相鄰篩孔尺寸,m;ψ表示篩分粒度與體積當量直徑比值,一般在1.15~1.30,無量綱。

    4 結(jié)論

    1)提出了一種顆粒沉降的阻力計算模型,當鉆井液為牛頓流體時,通過該模型可以預測顆粒雷諾數(shù)介于0~105的巖屑沉降末速度;當鉆井液為冪律流體時,可以預測顆粒雷諾數(shù)介于0~105、流性指數(shù)介于0.062 3~1的巖屑沉降末速度,這一范圍可以滿足鉆井工程中對巖屑沉降的預測需求。

    2)當顆粒在牛頓流體中沉降時,隨著顆粒雷諾數(shù)的增大,巖屑沉降所受黏性阻力逐漸減小,壓差阻力逐漸增大;分析表明:顆粒雷諾數(shù)小于2.944 6時黏性阻力大于壓差阻力,顆粒雷諾數(shù)大于2.944 6時壓差阻力大于黏性阻力;顆粒雷諾數(shù)小于1.11時,巖屑沉降主要考慮黏性阻力,介于1.11~500時,巖屑沉降受到黏性阻力與壓差阻力的共同作用,大于500時,壓差阻力占比達到90%,在巖屑沉降中占主導作用。

    3)顆粒在冪律流體中沉降時,所受到的黏性阻力和壓差阻力不僅與顆粒雷諾數(shù)相關,還與流性指數(shù)及稠度系數(shù)相關。

    [ 1 ] 伊薩克·牛頓. 自然哲學之數(shù)學原理[M]. 王克迪, 譯. 西安:陜西人民出版社, 2001.Newton I. The Principia: Mathematical principles of natural philosophy[M]. Wang Kedi trans. Xi'an: Shaanxi People's Publishing House, 2001.

    [ 2 ] Stokes GG. On the eあect of the internal friction of fl uids on the motion of pendulums[J]. Transactions of the Cambridge Philosophical Society, 1851, 9: 8-106.

    [ 3 ] Rittinger PR. Lehrbuch der aufbereitungskunde[M]. Berlin: Ernst& Kern, 1867.

    [ 4 ] 孫玉波. 重力選礦[M]. 北京: 冶金工業(yè)出版社, 1982.Sun Yubo. Gravity beneficiation[M]. Beijing: Metallurgical Industry Press, 1982.

    [ 5 ] 姚書典. 重選原理[M]. 北京: 冶金工業(yè)出版社, 1992.Yao Shudian. Gravity selection principle[M]. Beijing: Metallurgical Industry Press, 1992.

    [ 6 ] 謝廣元. 選礦學[M]. 徐州: 中國礦業(yè)大學出版社, 2001.Xie Guangyuan. Mineral processing[M]. Xuzhou: China University of Mining and Technology Press, 2001.

    [ 7 ] 錢寧, 萬兆惠. 泥沙運動力學[M]. 北京: 科學出版社, 1983.Qian Ning & Wan Zhaohui. Mechanics of sediment transport[M].Beijing: Science Press, 1983.

    [ 8 ] Peden JM & Luo Yuejin. Settling velocity of variously shaped particles in drilling and fracturing fl uids[J]. SPE Drilling Engineering, 1987, 2(4): 337-343.

    [ 9 ] Cheng Niansheng. Simplif i ed settling velocity formula for sediment particle[J]. Journal of Hydraulic Engineering, 1997, 123(2):149-152.

    [10] Guo Junke. Logarithmic matching and its applications in computational hydraulics and sediment transport[J]. Journal of Hydraulic Research, 2002, 40(5): 555-565.

    [11] Shah SN, EI Fadili Y & Chhabra RP. New model for single spherical particle settling velocity in power law (visco-inelastic) fl uids[J]. International Journal of Multiphase Flow, 2007, 33(1): 51-66.

    [12] Camenen B. Simple and general formula for the settling velocity of particles[J]. Journal of Hydraulic Engineering, 2007, 133(2):229-233.

    [13] Kao DTY & Hwang ALY. Determination of particle settling velocity in Heterogeneous suspensions and its eあects on energy loss prediction in solid-liquid freight pipelines[J]. Journal of Powder and Bulk Solids Technology, 1980, 4(1): 31-40.

    [14] 桂夏輝, 李延鋒, 劉炯天, 王永田, 曹亦俊. 液固流化床內(nèi)顆粒沉降特性試驗研究[J]. 煤炭學報, 2010, 35(8): 1374-1379.Gui Xiahui, Li Yanfeng, Liu Jiongtian, Wang Yongtian & Cao Yijun. Study on settlement characteristic of the grain in fl uidized bed[J]. Journal of China Coal Society, 2010, 35(8): 1374-1379.

    [15] Chhabra RP. Steady non-Newtonian fl ow about a rigid sphere[J].Encyclopedia of Fluid Mechanics, 1986, 1: 983-1033.

    [16]Chhabra RP. Motion of spheres in power law (viscoinelastic)fl uids at intermediate Reynolds numbers, a unif i ed approach[J].Chemical Engineering and Processing: Process Intensification,1990, 28(2): 89-94.

    [17] Dallon DS. A drag-coeきcient correlation for spheres settling in Ellis fl uids[D]. Salt Lake City, UT: University of Utah, 1967.

    [18] Lali AM, Khare AS, Joshi JB & Nigam KDP. Behaviour of solid particles in viscous non-Newtonian solutions: Settling velocity,wall eあects and bed expansion in solid–liquid fl uidized beds[J].Powder Technology, 1989, 57(1): 39-50.

    [19] Prakash S. Experimental evaluation of terminal velocity in non-Newtonian fl uids in the turbulent region[J]. Indian Chemical Engineer, 1983, 25: 1-4.

    [20] Shah SN. Proppant settling correlations for non-Newtonian fl uids under static and dynamic conditions[J]. Society of Petroleum Engineers Journal, 1982, 22(2): 164-170.

    [21] Shah SN. Proppant settling correlations for non-Newtonian fl uids[J]. SPE Production Engineering, 1986, 1(6): 446-448.

    [22] Shahi S & Kuru E. An experimental investigation of settling velocity of natural sands in water using particle image shadowgraph[J]. Powder Technology, 2015, 281: 184-192.

    [23] Shahi S & Kuru E. Experimental investigation of the settling velocity of spherical particles in power-law fl uids using particle image shadowgraph technique[J]. International Journal of Mineral Processing, 2016, 153: 60-65.

    [24] Terfous A, Hazzab A & Ghenaim A. Predicting the drag coeき-cient and settling velocity of spherical particles[J]. Powder Technology, 2013, 239: 12-20.

    [25] Elgaddafi R, Ahmed R, George M & Growcock F. Settling behavior of spherical particles in fi ber-containing drilling fl uids[J].Journal of Petroleum Science and Engineering, 2012, 84-85: 20-28.

    [26] Elgaddaf i R, Ahmed R & Growcock F. Settling behavior of particles in fi ber-containing Herschel Вulkley fl uid[J]. Powder Technology, 2016, 301: 782-793.

    [27] Bro?ek M & Surowiak A. Argument of separation at upgrading in the jig[J]. Archives of Mining Sciences, 2010, 55(1): 21-40.

    [28] Surowiak A & Bro?ek M. Methodology of calculation the terminal settling velocity distribution of spherical particles for high values of the Reynold's number[J]. Archives of Mining Sciences,2014, 59(1): 269-282.

    [29] Kelessidis VC. Terminal velocity of solid spheres falling in Newtonian and non-Newtonian liquids[J]. Tech. Chron. Sci., 2003,24(1/2): 43-54.

    [30] Kelessidis VC & Mpandelis G. Measurements and prediction of terminal velocity of solid spheres falling through stagnant pseudoplastic liquids[J]. Powder Technology, 2004, 147(1/3): 117-125.

    [31] Ford JT, Oyeneyin MB, Gao EH, Williamson RS & Peel LC.The formulation of milling fluids for efficient hole cleaning:An experimental investigation[C]//European Petroleum Conference, 25-27 October 1994, London, UK.DOI:https://dx.doi.org/10.2118/28819-MS.

    [32] Miura H, Takahashi T, Ichikawa J & Kawase Y. Bed expansion in liquid–solid two-phase fluidized beds with Newtonian and non-Newtonian fluids over the wide range of Reynolds numbers[J]. Powder Technology, 2001, 117(3): 239-246.

    [33] Reynolds PA & Jones TER. An experimental study of the settling velocities of single particles in non-Newtonian fl uids[J]. International Journal of Mineral Processing, 1989, 25(1/2): 47-77.

    猜你喜歡
    末速度雷諾數(shù)巖屑
    巖屑床破壞器在水平井斜井段的清潔效果研究
    論細碎巖屑地質(zhì)錄井
    勻變速直線運動的速度與位移的關系
    基于Transition SST模型的高雷諾數(shù)圓柱繞流數(shù)值研究
    K 近鄰分類法在巖屑數(shù)字圖像巖性分析中的應用
    錄井工程(2017年1期)2017-07-31 17:44:42
    巖屑實物錄井成果網(wǎng)上應用研究
    錄井工程(2017年1期)2017-07-31 17:44:38
    逆向思維求解勻減速到零問題中的“陷阱”
    考試周刊(2017年4期)2017-01-19 16:39:59
    失穩(wěn)初期的低雷諾數(shù)圓柱繞流POD-Galerkin 建模方法研究
    基于轉(zhuǎn)捩模型的低雷諾數(shù)翼型優(yōu)化設計研究
    勻變速直線運動的速度與位移的關系
    麻豆乱淫一区二区| 国产探花在线观看一区二区| 女的被弄到高潮叫床怎么办| 天美传媒精品一区二区| 亚洲中文字幕日韩| 联通29元200g的流量卡| 两个人视频免费观看高清| 国产男靠女视频免费网站| 狠狠狠狠99中文字幕| 久久精品国产亚洲av涩爱 | 午夜久久久久精精品| 久久久久免费精品人妻一区二区| 波多野结衣高清无吗| 美女黄网站色视频| 亚洲欧美清纯卡通| 国产又黄又爽又无遮挡在线| 熟女电影av网| 91午夜精品亚洲一区二区三区| 中文字幕av成人在线电影| 欧美不卡视频在线免费观看| 俄罗斯特黄特色一大片| 国产精品久久久久久久久免| 我要搜黄色片| 插逼视频在线观看| 在线观看一区二区三区| 卡戴珊不雅视频在线播放| 最新在线观看一区二区三区| 国产午夜精品论理片| 中文字幕av成人在线电影| 国产欧美日韩精品亚洲av| 性色avwww在线观看| 美女高潮的动态| 又黄又爽又刺激的免费视频.| 直男gayav资源| 韩国av在线不卡| 在线播放国产精品三级| 日日摸夜夜添夜夜爱| 国产成人精品久久久久久| 尾随美女入室| 日韩大尺度精品在线看网址| 国产熟女欧美一区二区| 色综合亚洲欧美另类图片| 最近中文字幕高清免费大全6| 久久人人爽人人片av| 国产综合懂色| 在线天堂最新版资源| 一a级毛片在线观看| 免费大片18禁| 亚洲av不卡在线观看| 国产午夜精品论理片| 蜜臀久久99精品久久宅男| 国产成人福利小说| 午夜激情福利司机影院| 中文资源天堂在线| 亚洲国产精品sss在线观看| 舔av片在线| 可以在线观看的亚洲视频| 丝袜喷水一区| 1000部很黄的大片| 国产男人的电影天堂91| 亚洲av一区综合| 天美传媒精品一区二区| 在线免费观看的www视频| 一边摸一边抽搐一进一小说| 亚洲熟妇中文字幕五十中出| av天堂在线播放| 我的老师免费观看完整版| 国产免费男女视频| 大又大粗又爽又黄少妇毛片口| 午夜免费男女啪啪视频观看 | 久久久久九九精品影院| 校园人妻丝袜中文字幕| 一级毛片电影观看 | 欧美日韩国产亚洲二区| 欧美国产日韩亚洲一区| 国产中年淑女户外野战色| 91精品国产九色| 亚洲图色成人| 亚洲国产高清在线一区二区三| 日日撸夜夜添| 男人的好看免费观看在线视频| 久久精品国产鲁丝片午夜精品| .国产精品久久| avwww免费| 欧美在线一区亚洲| 观看免费一级毛片| 久久热精品热| 观看美女的网站| 成年女人毛片免费观看观看9| 熟妇人妻久久中文字幕3abv| 国产一区二区在线观看日韩| 嫩草影院新地址| 久久草成人影院| 亚洲最大成人中文| 99久久中文字幕三级久久日本| 99久久精品国产国产毛片| 亚洲av成人精品一区久久| 亚洲国产高清在线一区二区三| 波多野结衣高清作品| 1024手机看黄色片| 国产伦在线观看视频一区| 午夜福利18| 欧美激情在线99| 特大巨黑吊av在线直播| 欧美高清性xxxxhd video| 国产精品av视频在线免费观看| 97超视频在线观看视频| 99国产精品一区二区蜜桃av| 亚洲国产精品成人综合色| 欧美日本视频| 免费观看的影片在线观看| 亚洲熟妇中文字幕五十中出| 亚洲综合色惰| 国产不卡一卡二| 国产精品av视频在线免费观看| 亚洲美女黄片视频| 少妇的逼水好多| 久久久久性生活片| 成人性生交大片免费视频hd| 97热精品久久久久久| 嫩草影院新地址| 国产亚洲av嫩草精品影院| 91久久精品国产一区二区成人| 国产91av在线免费观看| 成人毛片a级毛片在线播放| 九色成人免费人妻av| 国产成人影院久久av| 热99re8久久精品国产| 午夜爱爱视频在线播放| 免费看av在线观看网站| 久久久久久久亚洲中文字幕| 成人美女网站在线观看视频| 97超碰精品成人国产| 又黄又爽又免费观看的视频| 1000部很黄的大片| 国产精品久久久久久亚洲av鲁大| 国产精品av视频在线免费观看| 国产一区二区三区在线臀色熟女| 18禁在线无遮挡免费观看视频 | 少妇熟女欧美另类| 国产 一区 欧美 日韩| 在现免费观看毛片| 国产欧美日韩一区二区精品| 综合色丁香网| 亚洲国产欧洲综合997久久,| 国产一区亚洲一区在线观看| 丰满乱子伦码专区| 嫩草影院新地址| 国产精华一区二区三区| 亚洲一区高清亚洲精品| 成人美女网站在线观看视频| 国产 一区精品| а√天堂www在线а√下载| 色av中文字幕| 在线观看66精品国产| 3wmmmm亚洲av在线观看| 天美传媒精品一区二区| 啦啦啦啦在线视频资源| 亚洲无线观看免费| 一a级毛片在线观看| videossex国产| 尤物成人国产欧美一区二区三区| 日本熟妇午夜| 2021天堂中文幕一二区在线观| 男女做爰动态图高潮gif福利片| 色播亚洲综合网| 亚洲av成人精品一区久久| 国产精品久久久久久亚洲av鲁大| 亚洲精品国产av成人精品 | 国产一区二区激情短视频| 欧美丝袜亚洲另类| 亚洲av五月六月丁香网| 99热只有精品国产| 欧美色欧美亚洲另类二区| 欧美不卡视频在线免费观看| 国产欧美日韩精品亚洲av| 欧美丝袜亚洲另类| 免费在线观看影片大全网站| 欧美成人a在线观看| 99久国产av精品| av福利片在线观看| 国产蜜桃级精品一区二区三区| 欧美一区二区国产精品久久精品| 男女做爰动态图高潮gif福利片| 午夜福利在线在线| 两个人的视频大全免费| 国产探花极品一区二区| h日本视频在线播放| 午夜福利在线在线| 日本精品一区二区三区蜜桃| 俄罗斯特黄特色一大片| 99在线视频只有这里精品首页| 亚洲久久久久久中文字幕| 国内精品一区二区在线观看| 国产精品三级大全| 国产91av在线免费观看| 99国产极品粉嫩在线观看| 美女内射精品一级片tv| 亚洲三级黄色毛片| 99国产极品粉嫩在线观看| 亚洲五月天丁香| 亚洲在线自拍视频| 久久这里只有精品中国| 最新在线观看一区二区三区| 国产高清视频在线播放一区| 久久6这里有精品| 老师上课跳d突然被开到最大视频| 国内精品一区二区在线观看| 3wmmmm亚洲av在线观看| 成人特级av手机在线观看| a级毛片a级免费在线| 国产不卡一卡二| 国产精品一区二区性色av| 国产一区二区在线观看日韩| 国产高潮美女av| 久久99热这里只有精品18| 超碰av人人做人人爽久久| 嫩草影院新地址| 看非洲黑人一级黄片| 国产午夜精品久久久久久一区二区三区 | 日日干狠狠操夜夜爽| 国产黄色视频一区二区在线观看 | 大又大粗又爽又黄少妇毛片口| 一级a爱片免费观看的视频| 性插视频无遮挡在线免费观看| 久久精品国产鲁丝片午夜精品| 赤兔流量卡办理| 欧美最新免费一区二区三区| 全区人妻精品视频| 一区二区三区高清视频在线| 亚洲激情五月婷婷啪啪| 久久亚洲精品不卡| 欧美在线一区亚洲| 内地一区二区视频在线| 日本熟妇午夜| 日本与韩国留学比较| 狠狠狠狠99中文字幕| 成人一区二区视频在线观看| 免费av毛片视频| 色视频www国产| 在线天堂最新版资源| 日本免费a在线| 国产 一区 欧美 日韩| 亚洲国产精品sss在线观看| 99在线视频只有这里精品首页| 黄色一级大片看看| 舔av片在线| 国产精品永久免费网站| 三级国产精品欧美在线观看| 国产国拍精品亚洲av在线观看| 97热精品久久久久久| 桃色一区二区三区在线观看| 成熟少妇高潮喷水视频| 蜜桃久久精品国产亚洲av| 亚洲精华国产精华液的使用体验 | 一本精品99久久精品77| 内射极品少妇av片p| 婷婷色综合大香蕉| 国产黄色视频一区二区在线观看 | 国产成人一区二区在线| 床上黄色一级片| 国产精品久久电影中文字幕| 色在线成人网| 日本撒尿小便嘘嘘汇集6| 欧美三级亚洲精品| 六月丁香七月| 亚洲最大成人中文| 少妇裸体淫交视频免费看高清| 成人美女网站在线观看视频| 久久这里只有精品中国| 亚洲一区高清亚洲精品| 日本五十路高清| 国产亚洲精品久久久com| 色综合色国产| 午夜亚洲福利在线播放| 成人一区二区视频在线观看| 国产男人的电影天堂91| 99热这里只有是精品50| 免费av毛片视频| 老司机午夜福利在线观看视频| 免费在线观看影片大全网站| 一本一本综合久久| 插逼视频在线观看| 成人无遮挡网站| 日韩,欧美,国产一区二区三区 | 天堂√8在线中文| 中文字幕免费在线视频6| 又粗又爽又猛毛片免费看| 国产精华一区二区三区| 在线天堂最新版资源| 丝袜美腿在线中文| 欧美绝顶高潮抽搐喷水| 免费观看人在逋| 小蜜桃在线观看免费完整版高清| 国产亚洲91精品色在线| 精品久久久噜噜| 91av网一区二区| 一级毛片我不卡| 国产色婷婷99| 在线a可以看的网站| 国产欧美日韩精品一区二区| 免费黄网站久久成人精品| 亚洲人成网站在线观看播放| 又粗又爽又猛毛片免费看| 激情 狠狠 欧美| 亚洲欧美精品综合久久99| 18+在线观看网站| 午夜影院日韩av| 69av精品久久久久久| 在线a可以看的网站| 综合色丁香网| 内射极品少妇av片p| 亚洲最大成人av| 日韩大尺度精品在线看网址| 国产男人的电影天堂91| 成人综合一区亚洲| 日韩欧美一区二区三区在线观看| 国产黄色视频一区二区在线观看 | 日本精品一区二区三区蜜桃| 日日干狠狠操夜夜爽| 韩国av在线不卡| 一卡2卡三卡四卡精品乱码亚洲| 熟妇人妻久久中文字幕3abv| 国产午夜精品久久久久久一区二区三区 | 我的女老师完整版在线观看| 岛国在线免费视频观看| 观看免费一级毛片| 人人妻,人人澡人人爽秒播| 国内精品宾馆在线| 日本黄色视频三级网站网址| 91久久精品国产一区二区成人| 不卡视频在线观看欧美| 免费看美女性在线毛片视频| 可以在线观看毛片的网站| 色在线成人网| 99久国产av精品国产电影| 啦啦啦韩国在线观看视频| 精品一区二区三区人妻视频| 麻豆国产97在线/欧美| 男插女下体视频免费在线播放| 成人一区二区视频在线观看| 嫩草影视91久久| 午夜亚洲福利在线播放| 精品日产1卡2卡| 在线观看66精品国产| 国产av不卡久久| 啦啦啦啦在线视频资源| 亚洲av中文av极速乱| 精品一区二区三区视频在线| 国产精品一区二区性色av| 搡老妇女老女人老熟妇| 亚洲av成人精品一区久久| 中文字幕av成人在线电影| 久久6这里有精品| 久久久精品大字幕| 亚洲精品一区av在线观看| 久久久久久国产a免费观看| 免费人成视频x8x8入口观看| 国产在线精品亚洲第一网站| 三级国产精品欧美在线观看| 有码 亚洲区| 亚洲va在线va天堂va国产| 久久久久免费精品人妻一区二区| 亚洲精品成人久久久久久| 最近视频中文字幕2019在线8| 成人亚洲精品av一区二区| 悠悠久久av| 九九热线精品视视频播放| 一本精品99久久精品77| 日本精品一区二区三区蜜桃| 联通29元200g的流量卡| 九九久久精品国产亚洲av麻豆| 欧美xxxx性猛交bbbb| 日本成人三级电影网站| 激情 狠狠 欧美| 天堂影院成人在线观看| 在线天堂最新版资源| 国产精品人妻久久久影院| 日韩 亚洲 欧美在线| av专区在线播放| 亚洲精品在线观看二区| 我要搜黄色片| 成年av动漫网址| 成人特级黄色片久久久久久久| 亚洲无线在线观看| 观看免费一级毛片| 精品久久久久久久久亚洲| 国产单亲对白刺激| 久久热精品热| 夜夜夜夜夜久久久久| 久久人人精品亚洲av| 免费看日本二区| 3wmmmm亚洲av在线观看| 久久久久久久久久久丰满| 久久久久久伊人网av| 午夜精品国产一区二区电影 | 听说在线观看完整版免费高清| 欧美中文日本在线观看视频| 99热精品在线国产| 国产一区二区三区av在线 | 日日摸夜夜添夜夜爱| 精品99又大又爽又粗少妇毛片| 在线天堂最新版资源| a级毛片a级免费在线| 久久久a久久爽久久v久久| 亚洲国产精品sss在线观看| 国产视频内射| 蜜臀久久99精品久久宅男| 麻豆成人午夜福利视频| 欧美最新免费一区二区三区| 亚洲精品粉嫩美女一区| 别揉我奶头~嗯~啊~动态视频| 国产高清三级在线| 国产熟女欧美一区二区| АⅤ资源中文在线天堂| 国产精品一二三区在线看| 在线天堂最新版资源| 成人漫画全彩无遮挡| 亚洲av成人精品一区久久| 人人妻,人人澡人人爽秒播| 69人妻影院| 国产欧美日韩精品一区二区| 久久久精品94久久精品| 在线国产一区二区在线| 狠狠狠狠99中文字幕| 国语自产精品视频在线第100页| 久久久久国内视频| 一区二区三区高清视频在线| 色视频www国产| 婷婷色综合大香蕉| 老师上课跳d突然被开到最大视频| 五月玫瑰六月丁香| 免费av观看视频| 日本爱情动作片www.在线观看 | 18禁在线无遮挡免费观看视频 | 夜夜夜夜夜久久久久| 欧美日韩精品成人综合77777| 久久精品国产亚洲av天美| 免费不卡的大黄色大毛片视频在线观看 | 精品久久久久久久久亚洲| 高清毛片免费看| 亚州av有码| 免费人成在线观看视频色| 九九热线精品视视频播放| 欧美人与善性xxx| 免费搜索国产男女视频| 日本黄色视频三级网站网址| 午夜激情欧美在线| 少妇被粗大猛烈的视频| 午夜久久久久精精品| 变态另类丝袜制服| 国产精品99久久久久久久久| ponron亚洲| 亚洲欧美日韩高清专用| 2021天堂中文幕一二区在线观| 欧美成人免费av一区二区三区| 国产成人91sexporn| 中文字幕久久专区| 超碰av人人做人人爽久久| 欧美绝顶高潮抽搐喷水| 美女xxoo啪啪120秒动态图| 波多野结衣巨乳人妻| 日韩三级伦理在线观看| 亚洲国产精品国产精品| 精品久久久久久久久亚洲| 97超级碰碰碰精品色视频在线观看| 插逼视频在线观看| 日本免费a在线| 亚洲欧美成人精品一区二区| 国产精品人妻久久久影院| 成年av动漫网址| 五月伊人婷婷丁香| 久久人人爽人人片av| 精品一区二区三区人妻视频| 97热精品久久久久久| 亚洲欧美日韩东京热| 一级毛片久久久久久久久女| 国产亚洲精品av在线| 校园人妻丝袜中文字幕| 高清午夜精品一区二区三区 | 国产av一区在线观看免费| 亚洲内射少妇av| 午夜免费男女啪啪视频观看 | 国产免费一级a男人的天堂| 国产日本99.免费观看| 成人国产麻豆网| 成人综合一区亚洲| 日本免费a在线| 国产老妇女一区| 别揉我奶头~嗯~啊~动态视频| 无遮挡黄片免费观看| 久久久久久久久久久丰满| 亚洲精品粉嫩美女一区| 看免费成人av毛片| 99riav亚洲国产免费| 此物有八面人人有两片| 久久精品国产亚洲av天美| 亚洲一级一片aⅴ在线观看| 在线播放国产精品三级| 色综合亚洲欧美另类图片| 国产大屁股一区二区在线视频| 久久精品国产清高在天天线| 亚洲国产精品国产精品| 精品久久久久久久久亚洲| 欧美成人精品欧美一级黄| 国产真实伦视频高清在线观看| 麻豆av噜噜一区二区三区| 一区二区三区高清视频在线| 日本欧美国产在线视频| 内射极品少妇av片p| 一级毛片电影观看 | 老熟妇仑乱视频hdxx| 久久久色成人| .国产精品久久| 91在线精品国自产拍蜜月| 亚洲成av人片在线播放无| 亚洲精品国产成人久久av| 波多野结衣巨乳人妻| 久久久精品94久久精品| 国内精品久久久久精免费| 成年女人永久免费观看视频| 麻豆一二三区av精品| 狠狠狠狠99中文字幕| 亚洲一区二区三区色噜噜| 久久精品国产99精品国产亚洲性色| 久久久久久久久久黄片| 波多野结衣巨乳人妻| av卡一久久| 亚洲国产高清在线一区二区三| 国产精品一二三区在线看| 久久精品国产清高在天天线| 搡老熟女国产l中国老女人| 亚洲va在线va天堂va国产| 亚洲真实伦在线观看| 国内久久婷婷六月综合欲色啪| 欧美性猛交黑人性爽| 99久久九九国产精品国产免费| 精品欧美国产一区二区三| 91av网一区二区| 国内精品久久久久精免费| 91在线精品国自产拍蜜月| 悠悠久久av| 成年免费大片在线观看| 免费看a级黄色片| 精品日产1卡2卡| 亚洲欧美成人综合另类久久久 | 波多野结衣高清无吗| 人妻少妇偷人精品九色| 国产高清视频在线播放一区| 国产亚洲精品久久久com| 插逼视频在线观看| 久久九九热精品免费| 日本黄大片高清| 亚洲国产精品成人久久小说 | 天天躁夜夜躁狠狠久久av| 精品99又大又爽又粗少妇毛片| 国产一区二区激情短视频| 人人妻人人澡人人爽人人夜夜 | 国产精品一区二区性色av| 精品午夜福利视频在线观看一区| 亚洲精品粉嫩美女一区| 黄色日韩在线| 亚洲成a人片在线一区二区| 99热网站在线观看| 精品国内亚洲2022精品成人| 91久久精品国产一区二区三区| 熟妇人妻久久中文字幕3abv| 精品免费久久久久久久清纯| 看十八女毛片水多多多| 嫩草影院入口| 亚洲专区国产一区二区| 六月丁香七月| 伦理电影大哥的女人| 天美传媒精品一区二区| 两性午夜刺激爽爽歪歪视频在线观看| 我要搜黄色片| 97超碰精品成人国产| 久久久精品欧美日韩精品| 伦精品一区二区三区| 深爱激情五月婷婷| 欧美日韩国产亚洲二区| 国产高清有码在线观看视频| 一个人看的www免费观看视频| 久久人人爽人人片av| 日韩欧美三级三区| 久久久久国产精品人妻aⅴ院| 久久精品夜色国产| 国产精品免费一区二区三区在线| 特级一级黄色大片| 精品人妻熟女av久视频| 久久久久久久久大av| av在线亚洲专区| 国产淫片久久久久久久久| 久久精品国产亚洲av香蕉五月| 亚洲无线观看免费| 国产淫片久久久久久久久| 国产蜜桃级精品一区二区三区| 亚洲性久久影院| 日本欧美国产在线视频| 久久精品夜夜夜夜夜久久蜜豆| 国产亚洲精品av在线| 精品欧美国产一区二区三| 久久这里只有精品中国| 免费一级毛片在线播放高清视频| 一个人免费在线观看电影| 欧美区成人在线视频| 好男人在线观看高清免费视频| av在线观看视频网站免费| 97超级碰碰碰精品色视频在线观看| 美女免费视频网站| 蜜桃久久精品国产亚洲av| 久久久久久伊人网av| 亚洲在线观看片| 亚洲18禁久久av| 日韩在线高清观看一区二区三区|