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

    熱油管道油溫波動(dòng)隨機(jī)數(shù)值模擬及影響因素敏感性分析

    2011-09-28 06:03:56張文軻張勁軍
    關(guān)鍵詞:油溫進(jìn)站標(biāo)準(zhǔn)差

    張文軻,張勁軍,宇 波

    (中國(guó)石油大學(xué)城市油氣輸配技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室,北京102249)

    熱油管道油溫波動(dòng)隨機(jī)數(shù)值模擬及影響因素敏感性分析

    張文軻,張勁軍,宇 波

    (中國(guó)石油大學(xué)城市油氣輸配技術(shù)北京市重點(diǎn)實(shí)驗(yàn)室,北京102249)

    綜合采用有限容積法、有限差分法、Monte Carlo算法和POD算法建立埋地?zé)嵊凸艿姥鼐€油溫的隨機(jī)數(shù)值模擬算法,使用Sobol全局敏感性指標(biāo)進(jìn)行敏感性分析,綜合評(píng)價(jià)出站油溫、流量、壓力、埋深、埋深處自然地溫、土壤導(dǎo)熱系數(shù)、油品黏度和密度的隨機(jī)波動(dòng)對(duì)管道沿線油溫波動(dòng)的影響。計(jì)算結(jié)果表明:進(jìn)站油溫模擬結(jié)果與現(xiàn)場(chǎng)油溫均值偏差在0.1℃以內(nèi),標(biāo)準(zhǔn)差的偏差為0.006~0.023℃;出站油溫、流量、埋深處地溫和油品黏度4個(gè)參數(shù)的不確定性對(duì)四堡進(jìn)站油溫波動(dòng)的敏感性指標(biāo)之和為77.44%,河西站的該指標(biāo)為80.86%,進(jìn)站溫度的隨機(jī)數(shù)值模擬中主要考慮這4個(gè)參數(shù)的不確定性即可。

    熱油管道;油溫波動(dòng);數(shù)值模擬;敏感性分析;蒙特卡羅算法;POD算法

    原油管道運(yùn)行規(guī)程規(guī)定熱油管道的原油最低進(jìn)站溫度宜高于所輸原油凝點(diǎn)3℃[1],但管道實(shí)際運(yùn)行中,油溫與原油凝點(diǎn)都是波動(dòng)的。使用可靠性方法評(píng)價(jià)進(jìn)站油溫安全性時(shí),需要首先模擬管道沿線油溫的波動(dòng)。由于影響油溫的參數(shù)較多,若數(shù)值模擬中考慮所有參數(shù)的不確定性,工作量與計(jì)算量都太大,實(shí)際上亦無(wú)必要。在基于可靠性的評(píng)價(jià)中,重點(diǎn)是考慮不確定性強(qiáng)且影響大的參數(shù)的不確定性,其他參數(shù)可按確定性參數(shù)處理[2]。影響因素的敏感性分析,就是確定具有不確定性的各參數(shù)對(duì)結(jié)果的影響程度,為后續(xù)的隨機(jī)數(shù)值模擬提供依據(jù)。目前熱油管道水力、熱力特性的數(shù)值模擬研究均采用確定性算法,不考慮有關(guān)影響因素的隨機(jī)特性[3-5]。筆者引入Monte Carlo算法建立隨機(jī)數(shù)值模擬算法,并通過(guò)POD(最佳正交分解)算法解決數(shù)值模擬與Monte Carlo模擬結(jié)合(即隨機(jī)數(shù)值模擬)所產(chǎn)生的計(jì)算量大、耗時(shí)長(zhǎng)的問(wèn)題,在分析出站流量、油溫等參數(shù)變化時(shí)沿線油溫隨時(shí)間的變化基礎(chǔ)上,結(jié)合算例分析進(jìn)站油溫波動(dòng)對(duì)各參數(shù)的敏感性。

    1 埋地?zé)嵊凸艿罒崃δP?/h2>

    埋地?zé)嵊凸艿罒崃ο到y(tǒng)包括管內(nèi)油流之間及油流與周圍環(huán)境之間的熱傳導(dǎo)過(guò)程,較為復(fù)雜,作如下假設(shè)[5]:管內(nèi)同一截面上的油溫均勻分布;熱油管道熱力影響區(qū)簡(jiǎn)化為矩形;忽略原油的軸向傳熱,將三維不穩(wěn)定傳熱問(wèn)題轉(zhuǎn)化為二維不穩(wěn)定傳熱問(wèn)題。在上述假設(shè)下,埋地?zé)嵊凸艿罒崃δP腿鐖D1所示。

    圖1 埋地?zé)嵊凸艿朗疽鈭DFig.1 Sketch map of buried hot crude pipeline

    管內(nèi)油流的連續(xù)性方程為

    動(dòng)量方程為

    能量方程為

    式中,τ為時(shí)間,s;ρ為密度,kg/m3;A為管流橫截面面積,m2;v為截面平均流速,m/s;θ為管道傾角,rad;p為管道橫截面壓力,Pa;λ為達(dá)西水力摩阻系數(shù);D為管道直徑,m;u為原油比內(nèi)能,J/kg;s為高程,m;h為原油比焓,J/kg;q為單位管壁面積上的散熱量,J/(m2·s)。

    2 參數(shù)不確定性的分析及處理

    對(duì)埋地?zé)嵊凸艿姥鼐€油溫波動(dòng)產(chǎn)生影響的、具

    3 隨機(jī)數(shù)值模擬算法

    分別采用極坐標(biāo)網(wǎng)格和三角形網(wǎng)格對(duì)管內(nèi)及土壤區(qū)域進(jìn)行離散,離散后得到的計(jì)算區(qū)域網(wǎng)格見(jiàn)圖2。分別采用有限容積法和有限差分法對(duì)相關(guān)水力、熱力方程進(jìn)行離散求解。

    當(dāng)運(yùn)行參數(shù)、環(huán)境參數(shù)、油品物性具有不確定性時(shí),為得到管道沿線油溫的不確定性結(jié)果,可以使用Monte Carlo算法進(jìn)行隨機(jī)模擬,其本質(zhì)是一種數(shù)值有較明顯不確定性的參數(shù)包括運(yùn)行參數(shù)(出站油溫、流量和壓力)、油品物性(油品密度、油品黏度等)和環(huán)境參數(shù)(埋深、埋深處自然地溫、土壤導(dǎo)熱系數(shù)等)。

    (1)參數(shù)間相關(guān)性影響。管道運(yùn)行參數(shù)間存在一定相關(guān)性(如流量與壓力),ISO16708[2]指出當(dāng)參數(shù)間相關(guān)性系數(shù)絕對(duì)值小于0.5時(shí),各參數(shù)可做獨(dú)立變量處理。在本研究中,可能具有最密切相關(guān)性的參數(shù)是流量與壓力。為了區(qū)分這兩個(gè)參數(shù)本身的隨機(jī)波動(dòng)與它們之間的相互影響,在進(jìn)行參數(shù)不確定性分析時(shí),在相同開(kāi)泵條件下選取壓力較為平穩(wěn)的時(shí)間段對(duì)流量進(jìn)行不確定性分析[6],以及在流量平穩(wěn)的情況下對(duì)壓力的波動(dòng)進(jìn)行分析。

    (2)管道運(yùn)行參數(shù)。首先對(duì)SCADA系統(tǒng)在線監(jiān)測(cè)數(shù)據(jù)進(jìn)行統(tǒng)計(jì)以確定各變量的分布,在此基礎(chǔ)上使用Monte Carlo算法生成模擬樣本。由于出站壓力變化后會(huì)產(chǎn)生瞬變壓力波,其沿線擾動(dòng)影響造成區(qū)內(nèi)壓力分布和流量的瞬時(shí)變化,進(jìn)而產(chǎn)生油溫的波動(dòng),因此計(jì)算中考慮了出站壓力瞬時(shí)波動(dòng)的影響。

    (3)管內(nèi)油流物性。由于組成的變化,管內(nèi)原油的物性可具有較明顯的不確定性。研究中將管內(nèi)油流劃分為1 km一段,認(rèn)為每一段原油的物性一致。油品密度可采用統(tǒng)計(jì)參數(shù)直接生成,同一溫度下黏度不確定性可由正態(tài)分布描述[7]。生成模擬樣本時(shí),油品黏度均值由黏溫關(guān)系得到,標(biāo)準(zhǔn)差根據(jù)統(tǒng)計(jì)結(jié)果取值。

    (4)環(huán)境參數(shù)。管道埋深處自然地溫隨季節(jié)變化,但這種變化不屬于可靠性分析中“不確定性”的范疇。本研究所考慮的“地溫不確定性”,主要是由于測(cè)量或統(tǒng)計(jì)模型產(chǎn)生的不確定性,屬“認(rèn)知不確定性”范疇,同時(shí)將管道不同位置處因土壤物性等的變化而出現(xiàn)的地溫波動(dòng)亦作為不確定性考慮。研究中將管道分段,為每一段生成符合一定統(tǒng)計(jì)模型的地溫和土壤導(dǎo)熱系數(shù)模擬樣本。試驗(yàn)方法。若隨機(jī)事件f(x)的概率為p(f),對(duì)隨機(jī)變量x抽樣m次,若f(x)發(fā)生的次數(shù)為M,則對(duì)任意給定誤差ε,抽樣次數(shù)足夠多時(shí)下式成立:

    圖2 計(jì)算區(qū)域網(wǎng)格Fig.2 Grid of calculation region

    本研究分別對(duì)每一個(gè)隨機(jī)變量抽樣20 000次,這使得數(shù)值求解時(shí)間很長(zhǎng),同時(shí)不同時(shí)刻沿線各點(diǎn)處土壤溫度場(chǎng)和沿線油溫的存儲(chǔ)量很大,因此引入POD(最佳正交分解)算法。其基本原理[8]為,設(shè)f(x)為屬于空間Ω上的某已知物理場(chǎng),若求得一組最優(yōu)正交基{φi},使得函數(shù)f(x)可以展開(kāi)為級(jí)數(shù)形式,稱{φi}為該物理場(chǎng)的一組基函數(shù),{ai}為對(duì)應(yīng)于基函數(shù)的一組譜系數(shù),即

    在求解過(guò)程中,在將部分時(shí)刻管道沿線油溫和土壤溫度場(chǎng)數(shù)據(jù)進(jìn)行POD分解后,僅需存儲(chǔ)前4階基函數(shù)和對(duì)應(yīng)譜系數(shù)即可快速重構(gòu)得到任意時(shí)刻油溫和土壤溫度場(chǎng)數(shù)據(jù),從而達(dá)到減少計(jì)算量,提高計(jì)算效率和存儲(chǔ)效率的目的。研究[9]表明,以計(jì)算1000組不同參數(shù)組合為例,在建立基函數(shù)庫(kù)后,計(jì)算效率提高20000倍,數(shù)據(jù)存儲(chǔ)量減少80%以上。

    4 影響因素的敏感性分析方法

    埋地?zé)嵊凸艿朗且粋€(gè)涉及眾多影響因素的復(fù)雜熱力系統(tǒng)。常用于復(fù)雜系統(tǒng)的影響因素敏感性分析方法有傅里葉幅值敏感性測(cè)試[10](FAST)和Sobol全局敏感性指標(biāo)法[11]。FAST方法需要對(duì)每一個(gè)參數(shù)構(gòu)建具有幅值的搜索函數(shù),從而通過(guò)比較模型輸出結(jié)果在特定頻率下的傅里葉頻譜得到一階敏感性指標(biāo)。Sobol指標(biāo)法可由Monte Carlo算法中隨機(jī)變量與對(duì)應(yīng)函數(shù)值的統(tǒng)計(jì)參數(shù)得到,無(wú)須專門(mén)構(gòu)建搜索函數(shù),因此Sobol方法對(duì)復(fù)雜系統(tǒng)適用性更強(qiáng)。

    Sobol全局敏感性指標(biāo)法的原理[12]如下:設(shè)f(x)是定義在n維空間In上的可積函數(shù),則f(x)的ANOVA表達(dá)式定義為

    各子函數(shù)滿足正交條件,且可表示為積分形式。

    f(x)及其子函數(shù)的方差可表示為

    全局敏感性指標(biāo)S的下標(biāo)稱為指標(biāo)的階數(shù)或維數(shù),所有指標(biāo)和為1。

    在管道油溫波動(dòng)影響因素的分析中,由式(6)~(8)推得各參數(shù)一階和二階敏感性指標(biāo)為

    式中,S1i為第i個(gè)參數(shù)的一階敏感性指標(biāo);S2ij為第i個(gè)和第j個(gè)參數(shù)共同作用的二階敏感性指標(biāo);σi和σj分別為第i和j個(gè)參數(shù)單獨(dú)作用下管內(nèi)油溫隨時(shí)間變化的標(biāo)準(zhǔn)偏差,℃;σij為第i和第j個(gè)參數(shù)共同作用下管內(nèi)油溫隨時(shí)間變化的標(biāo)準(zhǔn)偏差,℃;σ為各參數(shù)單獨(dú)作用下管內(nèi)油溫隨時(shí)間變化的標(biāo)準(zhǔn)偏差之和,℃。

    可使用一階敏感性指標(biāo)描述油溫波動(dòng)對(duì)各參數(shù)波動(dòng)幅度(即標(biāo)準(zhǔn)差)的敏感性,一階敏感性指標(biāo)越小,說(shuō)明油溫波動(dòng)對(duì)該參數(shù)波動(dòng)不敏感,即該參數(shù)按確定性參數(shù)處理(標(biāo)準(zhǔn)差取0)對(duì)油溫波動(dòng)的分析結(jié)果無(wú)顯著影響。

    5 沿線油溫隨機(jī)波動(dòng)的模擬結(jié)果分析

    管道運(yùn)行過(guò)程中,各隨機(jī)參數(shù)的變化都會(huì)帶來(lái)進(jìn)站油溫的波動(dòng),但不同參數(shù)對(duì)油溫波動(dòng)的影響不同。一般來(lái)講,出站油溫和流量的波動(dòng)是造成進(jìn)站油溫波動(dòng)的最主要因素,因此使用前述隨機(jī)數(shù)值模擬算法分析出站油溫和流量單獨(dú)作用下的管道沿線油溫的波動(dòng),并討論站間距對(duì)沿線油溫隨機(jī)波動(dòng)的影響。

    對(duì)某100 km站間進(jìn)行隨機(jī)模擬,出站油溫均值tout分別取30和50℃,標(biāo)準(zhǔn)偏差Δt分別取0.5、1和2℃,其他參數(shù)均取確定值。出站和進(jìn)站油溫模擬結(jié)果見(jiàn)圖3。

    沿線油溫隨時(shí)間變化的標(biāo)準(zhǔn)差見(jiàn)圖4??梢?jiàn)管內(nèi)油溫的波動(dòng)幅度隨里程的增加逐漸衰減。這主要是因?yàn)橥寥缹?duì)油流的“溫度調(diào)節(jié)”作用:同一位置處的土壤從溫度較高的油流吸熱,同時(shí)向溫度較低的油流放熱。溫度調(diào)節(jié)作用的強(qiáng)弱是多種因素共同作用的結(jié)果。站間距越長(zhǎng),進(jìn)站油溫標(biāo)準(zhǔn)差越小。另外,出站油溫較高時(shí),油溫與土壤溫差較大,也有利于油溫波動(dòng)的衰減。

    取出站流量均值Q為1 200、1 800m3/h,標(biāo)準(zhǔn)偏差分別取80、120、160m3/h,其他參數(shù)均取確定值,沿線油溫標(biāo)準(zhǔn)差見(jiàn)圖5??梢?jiàn)管內(nèi)油溫的標(biāo)準(zhǔn)差在出站后約20km即達(dá)到穩(wěn)定,此后隨里程增加變化不明顯,僅與流量的波動(dòng)幅度有關(guān)。

    圖5 不同輸量下沿線油溫標(biāo)準(zhǔn)差Fig.5 Standard deviation of oil temperatures at different flow rates

    6 實(shí)際算例及驗(yàn)證

    分別選取西部原油管道站間距較長(zhǎng)的鄯善-四堡站間和站間距較短的翠嶺-河西站間,對(duì)輸送哈國(guó)油時(shí)各參數(shù)的波動(dòng)對(duì)進(jìn)站油溫波動(dòng)的影響進(jìn)行模擬。鄯善-四堡站間距為237.5 km,翠嶺-河西站間距為64.7 km,埋深處地溫為5.4℃。管道外徑813 mm,壁厚11 mm,熱力影響區(qū)取10 m,地表對(duì)流換熱系數(shù)取20.6 W/(m2·℃)。

    運(yùn)行參數(shù)均值和標(biāo)準(zhǔn)差根據(jù)西部管道SCADA系統(tǒng)數(shù)據(jù)統(tǒng)計(jì)得到,根據(jù)研究[6]分布模型可選正態(tài)分布。油品物性取值基于本課題組現(xiàn)場(chǎng)測(cè)試及其統(tǒng)計(jì)結(jié)果,分布模型取正態(tài)分布。土壤導(dǎo)熱系數(shù)缺乏現(xiàn)場(chǎng)實(shí)測(cè)數(shù)據(jù),UsowiczB等[13]對(duì)1 440組土壤樣本的測(cè)試結(jié)果表明,導(dǎo)熱系數(shù)數(shù)據(jù)服從正態(tài)分布,標(biāo)準(zhǔn)差取0.1~0.15。計(jì)算中各參數(shù)取值見(jiàn)表1(各參數(shù)的分布模型均為正態(tài)分布模型),油品黏度均值由下式給出:

    式中,t為管內(nèi)油溫,℃;μ為油品黏度,Pa·s;γ為剪切率,s-1。

    表1 隨機(jī)變量分布模型及分布參數(shù)取值Table 1 Distribution models and model parameters of stochastic variables

    分別對(duì)上述兩站間的進(jìn)站油溫進(jìn)行隨機(jī)數(shù)值模擬,模擬時(shí)間為680 h,對(duì)比油溫模擬結(jié)果與現(xiàn)場(chǎng)實(shí)測(cè)油溫統(tǒng)計(jì)值,其中四堡站進(jìn)站油溫均值為23.8℃,實(shí)測(cè)值為23.9℃,河西站進(jìn)站油溫均值為21.9℃,實(shí)測(cè)值為21.8℃,偏差平均為0.1℃。四堡站進(jìn)站油溫標(biāo)準(zhǔn)差為0.159℃,實(shí)測(cè)值為0.153℃,偏差為0.006℃;河西站進(jìn)站油溫標(biāo)準(zhǔn)差為0.250℃,實(shí)測(cè)值為0.227℃,偏差為0.023℃。可見(jiàn)隨機(jī)數(shù)值模擬結(jié)果與現(xiàn)場(chǎng)數(shù)據(jù)統(tǒng)計(jì)結(jié)果吻合較好。

    分別計(jì)算鄯善-四堡和翠嶺-河西站間影響進(jìn)站油溫波動(dòng)的各參數(shù)的一階全局敏感性指標(biāo),結(jié)果見(jiàn)表2??梢?jiàn)出站油溫和流量對(duì)進(jìn)站溫度波動(dòng)的影響最為顯著,其敏感性指標(biāo)之和分別為51.83%和59.76%。其中鄯善-四堡站間距較長(zhǎng),進(jìn)站油溫波動(dòng)對(duì)流量敏感性最強(qiáng),而翠嶺-河西站間較短,故出站油溫影響最大。其余隨機(jī)參數(shù)對(duì)進(jìn)站油溫波動(dòng)的影響依次為埋深處自然地溫、油品黏度、出站壓力、土壤導(dǎo)熱系數(shù)、埋深和油品密度。土壤導(dǎo)熱系數(shù)和管道埋深的影響較弱,是因?yàn)榉治鲋兄饕紤]其在管道沿線變化的不確定性,因此油流流經(jīng)不同管段時(shí),土壤導(dǎo)熱系數(shù)和埋深的漲落對(duì)進(jìn)站油溫波動(dòng)的影響可以在相當(dāng)程度上相互抵消。

    綜上所述,在表1給出的取值條件下,出站油溫、流量、地溫、原油黏度4個(gè)參數(shù)的敏感性指標(biāo)之和分別達(dá)到77.44%和80.86%,因此對(duì)進(jìn)站油溫隨機(jī)模擬時(shí),重點(diǎn)考慮這4個(gè)參數(shù)的不確定性即可。

    表2 進(jìn)站油溫波動(dòng)敏感性分析結(jié)果Table 2 Sensitivity analysis results of inlet oil temperature fluctuations

    7 結(jié) 論

    (1)采用基于有限差分法、有限容積法、Monte Carlo算法和POD分解算法建立的隨機(jī)數(shù)值模擬算法,可準(zhǔn)確模擬出站油溫、流量、壓力、管道埋深、埋深處自然地溫、土壤導(dǎo)熱系數(shù)、油品黏度和密度波動(dòng)時(shí)管道沿線油溫的變化。進(jìn)站油溫均值與現(xiàn)場(chǎng)統(tǒng)計(jì)數(shù)據(jù)偏差為0.1℃,標(biāo)準(zhǔn)差與現(xiàn)場(chǎng)統(tǒng)計(jì)值偏差為0.006~0.023℃。

    (2)當(dāng)出站油溫波動(dòng)時(shí),沿線油溫的波動(dòng)幅度隨輸送距離的延長(zhǎng)逐漸衰減,站間距越長(zhǎng)、出站油溫越高,衰減幅度越大。當(dāng)流量隨機(jī)波動(dòng)時(shí),沿線油溫的波動(dòng)幅度增大,但在較短距離內(nèi)達(dá)到一個(gè)相對(duì)穩(wěn)定的水平。

    (3)西部管道鄯善-四堡站、翠嶺-河西站出站油溫、流量、埋深處自然地溫、油品黏度4個(gè)參數(shù)對(duì)進(jìn)站油溫波動(dòng)的貢獻(xiàn)率分別達(dá)到77.44%和80.86%(表1取值條件下),在管道進(jìn)站油溫波動(dòng)的數(shù)值模擬中,重點(diǎn)考慮這4個(gè)參數(shù)的不確定性即可。

    [1] 中華人民共和國(guó)發(fā)展和改革委員會(huì).SY-T 5536-2004,原油管道運(yùn)行規(guī)程[S].北京:石油工業(yè)出版社,2004.

    [2] ISO.ISO16708-2006 Petroleum and natural gas industries-pipeline transportation systems-reliability-based limit state methods[S].Switzerland:The International Organization for Standardization,2006.

    [3] 李偉,宇波,張勁軍,等.熱油管道防腐層大修期間熱力參數(shù)的數(shù)值模擬[J].中國(guó)石油大學(xué)學(xué)報(bào):自然科學(xué)版,2009,33(1):103-106.L IWei,YU Bo,ZHANG Jin-jun,et al.Numerical simulation of thermal parameter of hot oil pipeline during coating overhauling period[J].Journal of China University of Petroleum(Edition of Natural Science),2009,33(1):103-106.

    [4] 吳玉國(guó),陳保東.冷熱原油順序輸送彎管內(nèi)混油特性數(shù)值模擬[J].中國(guó)石油大學(xué)學(xué)報(bào):自然科學(xué)版,2010,34(3):125-128.WU Yu-guo,CHEN Bao-dong.Numerical simulation on contaminated oil characteristics in bend pipe of batch pipelining of cold and hot crude oils[J].Journal of China University of Petroleum(Edition of Natural Science),2010,34(3):125-128.

    [5] 王凱,張勁軍,宇波.冷熱原油交替順序輸送中加熱時(shí)機(jī)的經(jīng)濟(jì)比選[J].中國(guó)石油大學(xué)學(xué)報(bào):自然科學(xué)版,2008,32(5):102-105.WANG Kai,ZHANG Jin-jun,YU Bo.Optimal heating ratio of batch pipelining of cold and hot crude oils[J].Journal of China University of Petroleum(Edition of Natural Science),2008,32(5):102-105.

    [6] 張文軻,張勁軍,丁建林.中洛輸油管道輸量和出站溫度概率分布[J].油氣儲(chǔ)運(yùn),2010,29(1):32-35.ZHANGWen-ke,ZHANG Jin-jun,DING Jian-lin.Study on probability distribution of flowrate and outlet temperature over Zhongluo crude oil pipeline[J].OGST,2010,29(1):32-35.

    [7] 許康,張勁軍.管輸含蠟原油黏度的不確定性及概率分布研究[J].西安石油大學(xué)學(xué)報(bào),2004,19(5):35-38.XU Kang,ZHANG Jin-jun.Uncertainty and probability distribution of the viscosity of pipeline waxy crude oil[J].Journal of Xi′an Shiyou University,2004,19(5):35-38.[8] BERKOOZ G,HOLMES P,LUMLEY J L.The proper orthogonal decomposition in the analysis of turbulent flows[J].Annu Rev Fluid Mech,1993,7(25):539-575.

    [9] 張勁軍,張文軻,宇波.基于POD的埋地?zé)嵊凸艿劳]斣賳?dòng)模擬高效算法[J].石油學(xué)報(bào),2011,32(1):167-170.ZHANG Jin-jun,ZHANGWen-ke,YU bo.An efficient algorithm of the proper orthogonal decomposition for the restart computation of buried heated crude pipelines[J],Acta Petrolei Sinica,2011,32(1):167-170.

    [10] LU Yi-chi,SITAKANTA.Sensitivity analysis of a complex proposed geologic waste disposal system using the Fourier Amplitude Sensitivity Testmethod[J].Reliability Engineering and System Safety,2001,72(3):275-291.

    [11] XU Chong-gang,GERTNER G.A general first-order global sensitivity analysis method[J].Reliability Engineering and System Safety,2008,93(1):1060-1071.[12] SOBOL L,TARANTOLA S.Estimating the approximation error when fixing unessential factors in global sensitivity analysis[J].Reliability Engineering and System Safety,2007,92(1):957-960.

    [13] USOW I CZ B,LIPIEC J,FERRERO A.Prediction of soil thermal conductivity based on penetration resistance and water content or air-filled porosity[J].Int J Heat Mass Transfer,2006,49:5010-5018.

    (編輯 沈玉英)

    Stochastic numerical simulation on oil temperature fluctuations of hot crude pipelines and sensitivity analysis to related factors

    ZHANG Wen-ke,ZHANG Jin-jun,YU Bo

    (Beijing Key Laboratory of Urban Oil and Gas Distribution Technology,China University of Petroleum,Beijing102249,China)

    An efficient numerical algorithm,which employed the finite volume,finite difference,Monte Carlo and proper orthogonal decomposition(POD)method,was proposed to simulate the stochastic fluctuation of oil temperatures.The effects of outlet oil temperature,flow rate,pressure,buried depth,soil temperature,soil thermal conductivity,oil viscosity and oil density on the stochastic fluctuation of oil temperatures were investigated through Sobol global sensitivity induces.The results show that the deviations of the mean oil temperatures between numerical simulations and field measurements were lower than 0.1℃,and the standard deviations vary from 0.006℃to 0.023℃.When considered the following 4 factors:outlet oil temperature,throughput,soil temperature and oil viscosity,case studies on the China West crude pipeline show that the sum of first order global sensitivity achieved as 77.44%on Sipu station and 80.86%on Hexi station.Therefore,it is acceptable to consider these 4 factors when simulating the fluctuation of inlet oil temperatures.

    hot oilpipeline;oil temperature fluctuation;numerical simulation;sensitivity analysis;Monte Carlo algorithm;POD algorithm

    TE 973.1

    A

    10.3969/j.issn.1673-5005.2011.02.025

    2010-10-22

    國(guó)家自然科學(xué)基金項(xiàng)目(50944030);北京市優(yōu)秀博士學(xué)位論文指導(dǎo)教師科技項(xiàng)目(YB20081141401)

    張文軻(1982-),男(漢族),山東淄博人,博士研究生,主要從事油氣儲(chǔ)運(yùn)系統(tǒng)安全工程研究。

    1673-5005(2011)02-0141-06

    猜你喜歡
    油溫進(jìn)站標(biāo)準(zhǔn)差
    用Pro-Kin Line平衡反饋訓(xùn)練儀對(duì)早期帕金森病患者進(jìn)行治療對(duì)其動(dòng)態(tài)平衡功能的影響
    進(jìn)站口上下行載頻切換時(shí)引起ATP制動(dòng)問(wèn)題分析
    春運(yùn)期間北京西站共有154.8萬(wàn)人次刷臉進(jìn)站
    祖國(guó)(2018年6期)2018-06-27 10:27:26
    閱讀(科學(xué)探秘)(2018年8期)2018-05-14 10:06:29
    一根筷子辨別油溫
    液壓系統(tǒng)油溫過(guò)高現(xiàn)象原因分析及改進(jìn)
    如何判斷油溫?
    食品與健康(2015年1期)2015-09-10 07:22:44
    對(duì)于平均差與標(biāo)準(zhǔn)差的數(shù)學(xué)關(guān)系和應(yīng)用價(jià)值比較研究
    重慶軌道交通三號(hào)線列車進(jìn)站警示功能接口電路的分析
    醫(yī)學(xué)科技論文中有效數(shù)字的確定
    三级国产精品欧美在线观看| 我的老师免费观看完整版| 欧美激情在线99| 最近在线观看免费完整版| 亚洲久久久久久中文字幕| 欧美一级a爱片免费观看看| 国产伦一二天堂av在线观看| 白带黄色成豆腐渣| 久久国产乱子伦精品免费另类| 午夜福利在线在线| 色综合站精品国产| 美女 人体艺术 gogo| 国产又黄又爽又无遮挡在线| 女同久久另类99精品国产91| 高清日韩中文字幕在线| a在线观看视频网站| 91麻豆av在线| 老司机深夜福利视频在线观看| 免费在线观看影片大全网站| 国产国拍精品亚洲av在线观看 | 亚洲av电影不卡..在线观看| 久久久久国产精品人妻aⅴ院| 嫩草影视91久久| www.999成人在线观看| 久久精品91无色码中文字幕| 尤物成人国产欧美一区二区三区| av福利片在线观看| www.色视频.com| 91在线精品国自产拍蜜月 | av欧美777| 国产一区二区三区在线臀色熟女| 国内精品久久久久久久电影| 亚洲成av人片在线播放无| 精品一区二区三区人妻视频| 波多野结衣高清作品| 午夜福利在线观看吧| 成熟少妇高潮喷水视频| 亚洲av五月六月丁香网| 男女做爰动态图高潮gif福利片| 观看美女的网站| 午夜激情福利司机影院| 久久精品国产综合久久久| 日韩欧美 国产精品| 18禁国产床啪视频网站| 美女cb高潮喷水在线观看| 国产精品永久免费网站| 亚洲人成网站在线播放欧美日韩| 美女黄网站色视频| 国产成人系列免费观看| 少妇裸体淫交视频免费看高清| 午夜激情欧美在线| 中文资源天堂在线| 亚洲国产精品sss在线观看| 国产一区二区三区在线臀色熟女| 欧美黄色淫秽网站| 成人高潮视频无遮挡免费网站| 亚洲成人久久性| 国产精品久久久久久精品电影| 成人午夜高清在线视频| 亚洲国产精品久久男人天堂| 亚洲专区国产一区二区| 久久久久久久久久黄片| 好男人在线观看高清免费视频| 国产在视频线在精品| 国产精品98久久久久久宅男小说| 国产精品久久久人人做人人爽| 女人十人毛片免费观看3o分钟| 首页视频小说图片口味搜索| 婷婷亚洲欧美| avwww免费| 制服丝袜大香蕉在线| 亚洲欧美激情综合另类| 欧美精品啪啪一区二区三区| 久久婷婷人人爽人人干人人爱| 亚洲aⅴ乱码一区二区在线播放| 女人被狂操c到高潮| 搡老熟女国产l中国老女人| 伊人久久大香线蕉亚洲五| 久久香蕉国产精品| 无人区码免费观看不卡| 久久精品亚洲精品国产色婷小说| 色噜噜av男人的天堂激情| 51午夜福利影视在线观看| 欧美3d第一页| 亚洲色图av天堂| 一a级毛片在线观看| 日日干狠狠操夜夜爽| 首页视频小说图片口味搜索| 亚洲av二区三区四区| 男人舔女人下体高潮全视频| 午夜福利在线观看吧| 丰满人妻熟妇乱又伦精品不卡| 黄色成人免费大全| 亚洲精品美女久久久久99蜜臀| 亚洲成av人片免费观看| 少妇的逼水好多| 97碰自拍视频| 性色av乱码一区二区三区2| 欧美丝袜亚洲另类 | 国产av一区在线观看免费| 婷婷亚洲欧美| 中文字幕人妻丝袜一区二区| 美女 人体艺术 gogo| 18+在线观看网站| 亚洲精品亚洲一区二区| 制服丝袜大香蕉在线| 蜜桃亚洲精品一区二区三区| av天堂中文字幕网| 亚洲av成人av| 国产97色在线日韩免费| 老熟妇乱子伦视频在线观看| 久久久久久久久久黄片| 一级黄片播放器| 亚洲午夜理论影院| avwww免费| 国产精品久久视频播放| 看免费av毛片| 日本与韩国留学比较| 免费看美女性在线毛片视频| 亚洲国产精品久久男人天堂| 国产午夜福利久久久久久| 亚洲精品美女久久久久99蜜臀| 一级作爱视频免费观看| 国产熟女xx| 最好的美女福利视频网| 亚洲五月天丁香| 国产av一区在线观看免费| 亚洲 欧美 日韩 在线 免费| 一个人免费在线观看的高清视频| 成熟少妇高潮喷水视频| 精品国产超薄肉色丝袜足j| 99在线人妻在线中文字幕| 国产精品美女特级片免费视频播放器| 国产激情欧美一区二区| 在线视频色国产色| 午夜福利在线观看吧| 男女那种视频在线观看| 91字幕亚洲| 亚洲av电影在线进入| 女人被狂操c到高潮| 国产毛片a区久久久久| 九九久久精品国产亚洲av麻豆| a在线观看视频网站| 亚洲av成人精品一区久久| 身体一侧抽搐| 亚洲狠狠婷婷综合久久图片| 成人亚洲精品av一区二区| 黄色成人免费大全| 色综合站精品国产| 老司机福利观看| av福利片在线观看| 欧美一区二区国产精品久久精品| 色av中文字幕| 精品国产超薄肉色丝袜足j| 波多野结衣高清作品| 51国产日韩欧美| 国产精品亚洲美女久久久| 久久国产乱子伦精品免费另类| 亚洲国产精品成人综合色| 久久久久精品国产欧美久久久| 三级男女做爰猛烈吃奶摸视频| 黄片小视频在线播放| 可以在线观看的亚洲视频| 亚洲一区二区三区不卡视频| 欧美av亚洲av综合av国产av| 男人的好看免费观看在线视频| 亚洲在线自拍视频| ponron亚洲| 欧美性感艳星| 制服人妻中文乱码| 亚洲精品一区av在线观看| av视频在线观看入口| 亚洲欧美精品综合久久99| 日本成人三级电影网站| 久久久国产精品麻豆| 亚洲精品456在线播放app | 中文资源天堂在线| 丰满人妻一区二区三区视频av | 九九在线视频观看精品| 亚洲性夜色夜夜综合| 成熟少妇高潮喷水视频| 日本黄色视频三级网站网址| 亚洲片人在线观看| 亚洲精品在线美女| 国产精品 国内视频| 成人av在线播放网站| 法律面前人人平等表现在哪些方面| 在线播放无遮挡| 成年女人看的毛片在线观看| 日韩 欧美 亚洲 中文字幕| 国产黄a三级三级三级人| 最近最新免费中文字幕在线| 非洲黑人性xxxx精品又粗又长| 免费在线观看日本一区| 色哟哟哟哟哟哟| 天堂动漫精品| 亚洲国产色片| 日韩大尺度精品在线看网址| 精品一区二区三区人妻视频| 午夜免费观看网址| 成人午夜高清在线视频| 久久久久精品国产欧美久久久| 久久精品91蜜桃| 国产成人福利小说| 日本撒尿小便嘘嘘汇集6| 国产精品99久久久久久久久| 国产真人三级小视频在线观看| 国产黄色小视频在线观看| 精品日产1卡2卡| 国产精品乱码一区二三区的特点| 婷婷精品国产亚洲av在线| 亚洲无线在线观看| 一级黄片播放器| 欧美高清成人免费视频www| 91字幕亚洲| 少妇的逼好多水| www.www免费av| 色av中文字幕| 757午夜福利合集在线观看| 3wmmmm亚洲av在线观看| 国产三级黄色录像| 可以在线观看的亚洲视频| 精品久久久久久久人妻蜜臀av| 12—13女人毛片做爰片一| 麻豆国产97在线/欧美| 国产精品久久久久久人妻精品电影| 亚洲成a人片在线一区二区| 99精品在免费线老司机午夜| 成人av一区二区三区在线看| 亚洲国产欧美人成| 国产成人av教育| 香蕉久久夜色| 99在线人妻在线中文字幕| 熟女少妇亚洲综合色aaa.| 亚洲 国产 在线| 老鸭窝网址在线观看| 少妇的逼水好多| 国产精品一及| 蜜桃久久精品国产亚洲av| 成人国产一区最新在线观看| 久久午夜亚洲精品久久| 久久久久国产精品人妻aⅴ院| 久久香蕉国产精品| 夜夜躁狠狠躁天天躁| 精品人妻1区二区| 免费人成视频x8x8入口观看| 一本综合久久免费| a级毛片a级免费在线| 熟女少妇亚洲综合色aaa.| 欧美+亚洲+日韩+国产| 俺也久久电影网| 亚洲va日本ⅴa欧美va伊人久久| 亚洲内射少妇av| 亚洲国产中文字幕在线视频| 国产真实伦视频高清在线观看 | 久久精品影院6| 欧洲精品卡2卡3卡4卡5卡区| 国产一区在线观看成人免费| 日本黄色片子视频| 每晚都被弄得嗷嗷叫到高潮| 美女cb高潮喷水在线观看| 欧美日韩瑟瑟在线播放| 成人特级av手机在线观看| 国产高清videossex| 在线视频色国产色| 两人在一起打扑克的视频| 身体一侧抽搐| 国产亚洲精品久久久com| 国产精品一区二区三区四区免费观看 | 狠狠狠狠99中文字幕| 午夜精品一区二区三区免费看| 3wmmmm亚洲av在线观看| 桃色一区二区三区在线观看| 亚洲av成人不卡在线观看播放网| 国产精品嫩草影院av在线观看 | 天堂影院成人在线观看| 热99在线观看视频| 色av中文字幕| 国产一区在线观看成人免费| 国产精品久久久久久久久免 | 久久精品夜夜夜夜夜久久蜜豆| 欧美色欧美亚洲另类二区| 国产一级毛片七仙女欲春2| 母亲3免费完整高清在线观看| 中文资源天堂在线| 中文字幕熟女人妻在线| 少妇的逼水好多| 十八禁网站免费在线| 成人18禁在线播放| 男女做爰动态图高潮gif福利片| 国产97色在线日韩免费| 免费看a级黄色片| 欧美日韩福利视频一区二区| 不卡一级毛片| 给我免费播放毛片高清在线观看| 九九在线视频观看精品| 日韩高清综合在线| 一个人看的www免费观看视频| 听说在线观看完整版免费高清| 一级毛片女人18水好多| 亚洲av成人av| 人妻丰满熟妇av一区二区三区| 久久久久久久久大av| 一卡2卡三卡四卡精品乱码亚洲| 国产av在哪里看| 亚洲成人久久爱视频| 精品一区二区三区视频在线观看免费| a在线观看视频网站| 免费无遮挡裸体视频| 日本免费a在线| 午夜日韩欧美国产| 久久久久久国产a免费观看| 两个人看的免费小视频| 日韩欧美三级三区| 九色国产91popny在线| 九九久久精品国产亚洲av麻豆| 亚洲第一欧美日韩一区二区三区| 在线视频色国产色| 国产精品久久久久久久久免 | 久久99热这里只有精品18| 久久草成人影院| 伊人久久大香线蕉亚洲五| 久久久国产成人免费| 全区人妻精品视频| 天天躁日日操中文字幕| 少妇熟女aⅴ在线视频| 欧美性感艳星| 久久久久免费精品人妻一区二区| 免费看美女性在线毛片视频| 特级一级黄色大片| av在线天堂中文字幕| 成人一区二区视频在线观看| 亚洲av美国av| 亚洲天堂国产精品一区在线| 亚洲,欧美精品.| 非洲黑人性xxxx精品又粗又长| 变态另类成人亚洲欧美熟女| 床上黄色一级片| 免费看日本二区| 亚洲精品一卡2卡三卡4卡5卡| 丰满人妻一区二区三区视频av | 欧美高清成人免费视频www| 成人三级黄色视频| 国产亚洲精品久久久com| 少妇的逼水好多| av欧美777| 女警被强在线播放| 欧美日韩瑟瑟在线播放| 久久精品影院6| 校园春色视频在线观看| 99久久久亚洲精品蜜臀av| 亚洲av熟女| 日本一本二区三区精品| 国产三级中文精品| 母亲3免费完整高清在线观看| 精品久久久久久久久久免费视频| 亚洲国产精品999在线| 99国产精品一区二区蜜桃av| 国产成人a区在线观看| 日本与韩国留学比较| 高清毛片免费观看视频网站| 中文资源天堂在线| 九九久久精品国产亚洲av麻豆| 亚洲人成电影免费在线| 日韩精品中文字幕看吧| 国产一区二区激情短视频| 国产精品久久电影中文字幕| 亚洲 国产 在线| 一本久久中文字幕| 特级一级黄色大片| 床上黄色一级片| 欧美中文综合在线视频| 成人国产一区最新在线观看| 免费看十八禁软件| 亚洲一区二区三区不卡视频| 十八禁网站免费在线| 成熟少妇高潮喷水视频| 又粗又爽又猛毛片免费看| 欧美成人一区二区免费高清观看| 日韩欧美在线乱码| 亚洲,欧美精品.| 99久国产av精品| 成人国产一区最新在线观看| 国产亚洲欧美98| 国产真实伦视频高清在线观看 | 国产欧美日韩一区二区精品| 好男人电影高清在线观看| 99热6这里只有精品| 日本五十路高清| 亚洲一区高清亚洲精品| 午夜福利欧美成人| 国产在视频线在精品| 久久中文看片网| 最近最新中文字幕大全电影3| 黄片小视频在线播放| 精品一区二区三区av网在线观看| 久久香蕉国产精品| 香蕉久久夜色| 99久久无色码亚洲精品果冻| 亚洲中文字幕日韩| 在线观看日韩欧美| 午夜影院日韩av| 老司机午夜福利在线观看视频| 啦啦啦免费观看视频1| 女人被狂操c到高潮| www日本黄色视频网| 99热这里只有精品一区| 日韩欧美三级三区| 成人欧美大片| 日韩欧美免费精品| 国产精品99久久99久久久不卡| 成人无遮挡网站| 国产精品久久久久久久电影 | 日韩欧美精品v在线| 国产v大片淫在线免费观看| 精品国产三级普通话版| 午夜久久久久精精品| 美女高潮喷水抽搐中文字幕| 亚洲av中文字字幕乱码综合| 我要搜黄色片| 波野结衣二区三区在线 | 婷婷亚洲欧美| 桃红色精品国产亚洲av| 中文字幕av成人在线电影| av片东京热男人的天堂| 久久精品国产自在天天线| 黄色丝袜av网址大全| 国产欧美日韩一区二区三| 夜夜躁狠狠躁天天躁| 亚洲av日韩精品久久久久久密| 91九色精品人成在线观看| 欧美成人性av电影在线观看| 三级男女做爰猛烈吃奶摸视频| 免费观看的影片在线观看| 国产高清视频在线播放一区| aaaaa片日本免费| av国产免费在线观看| 丰满人妻熟妇乱又伦精品不卡| 欧美黄色片欧美黄色片| 中文字幕av在线有码专区| 一个人免费在线观看电影| 岛国在线观看网站| 免费观看的影片在线观看| 精品人妻偷拍中文字幕| 国产精品亚洲av一区麻豆| 久久亚洲精品不卡| 人妻丰满熟妇av一区二区三区| 亚洲性夜色夜夜综合| 欧美+日韩+精品| 丝袜美腿在线中文| 欧美日韩黄片免| 国产淫片久久久久久久久 | 国产免费av片在线观看野外av| www.熟女人妻精品国产| 哪里可以看免费的av片| 国产免费男女视频| 国模一区二区三区四区视频| 日韩免费av在线播放| 国产69精品久久久久777片| 精品国产三级普通话版| 国内精品久久久久精免费| 噜噜噜噜噜久久久久久91| 欧美日韩瑟瑟在线播放| 毛片女人毛片| 国产av在哪里看| 亚洲久久久久久中文字幕| 99精品欧美一区二区三区四区| 看黄色毛片网站| 两个人视频免费观看高清| 亚洲人成网站在线播放欧美日韩| 国产成人av教育| 亚洲精品成人久久久久久| 老师上课跳d突然被开到最大视频 久久午夜综合久久蜜桃 | 亚洲国产欧美人成| 欧美日韩瑟瑟在线播放| 亚洲av免费高清在线观看| 久久久久亚洲av毛片大全| 一卡2卡三卡四卡精品乱码亚洲| 日韩欧美三级三区| 韩国av一区二区三区四区| 亚洲国产精品成人综合色| 国产午夜福利久久久久久| 免费在线观看亚洲国产| 国产伦精品一区二区三区视频9 | 少妇丰满av| 久久草成人影院| 少妇人妻一区二区三区视频| 久久人妻av系列| 一二三四社区在线视频社区8| 日本一本二区三区精品| 天堂动漫精品| 久久精品国产亚洲av涩爱 | 欧美大码av| 国产私拍福利视频在线观看| 国产色爽女视频免费观看| 深夜精品福利| 制服人妻中文乱码| 中文字幕人妻丝袜一区二区| 免费在线观看日本一区| 亚洲精品美女久久久久99蜜臀| 免费人成在线观看视频色| 长腿黑丝高跟| 国产伦精品一区二区三区四那| 国产精品日韩av在线免费观看| 免费观看精品视频网站| 一进一出好大好爽视频| 脱女人内裤的视频| 亚洲欧美一区二区三区黑人| 国产一区二区三区视频了| 国产视频内射| 高清日韩中文字幕在线| 亚洲国产日韩欧美精品在线观看 | 欧美日韩黄片免| 色视频www国产| 高清日韩中文字幕在线| 国产熟女xx| 18禁裸乳无遮挡免费网站照片| 欧美+亚洲+日韩+国产| 亚洲男人的天堂狠狠| 国产欧美日韩精品亚洲av| 久久精品国产亚洲av涩爱 | 国产av一区在线观看免费| 变态另类丝袜制服| 欧美不卡视频在线免费观看| 日韩有码中文字幕| 国产精品香港三级国产av潘金莲| av专区在线播放| 免费人成视频x8x8入口观看| 亚洲专区中文字幕在线| 亚洲人成网站在线播| 老汉色av国产亚洲站长工具| 又黄又爽又免费观看的视频| 欧美日韩一级在线毛片| 亚洲片人在线观看| 日韩大尺度精品在线看网址| 精品免费久久久久久久清纯| 身体一侧抽搐| 成熟少妇高潮喷水视频| 国产三级黄色录像| 婷婷精品国产亚洲av| 中文字幕精品亚洲无线码一区| 99久久精品热视频| 在线播放无遮挡| 一本久久中文字幕| 亚洲欧美日韩无卡精品| 中文字幕久久专区| 欧美zozozo另类| 日本免费一区二区三区高清不卡| 亚洲国产中文字幕在线视频| 久久久久久大精品| 亚洲人成网站在线播放欧美日韩| 色视频www国产| 亚洲精品色激情综合| 天天一区二区日本电影三级| 久9热在线精品视频| 国产一区在线观看成人免费| 国产激情偷乱视频一区二区| 亚洲五月天丁香| 午夜免费成人在线视频| 两个人看的免费小视频| 亚洲七黄色美女视频| 12—13女人毛片做爰片一| 精品久久久久久成人av| a级毛片a级免费在线| 九九在线视频观看精品| 少妇的丰满在线观看| 亚洲熟妇熟女久久| 亚洲精品一卡2卡三卡4卡5卡| 久久久久性生活片| 亚洲av美国av| 欧美日韩瑟瑟在线播放| 国产成人av教育| 一进一出好大好爽视频| 又黄又粗又硬又大视频| 一级作爱视频免费观看| 国产精品久久久人人做人人爽| 91九色精品人成在线观看| 成年人黄色毛片网站| 色噜噜av男人的天堂激情| 亚洲中文日韩欧美视频| 黄片大片在线免费观看| eeuss影院久久| 久久久久久久久中文| 日本一二三区视频观看| 国产一区二区在线观看日韩 | 欧美黑人欧美精品刺激| 成人国产综合亚洲| 99热6这里只有精品| 美女高潮喷水抽搐中文字幕| 欧美激情久久久久久爽电影| 非洲黑人性xxxx精品又粗又长| 国产精品免费一区二区三区在线| 少妇的逼水好多| 岛国在线免费视频观看| 国产精品99久久99久久久不卡| 日韩亚洲欧美综合| av天堂在线播放| 精品99又大又爽又粗少妇毛片 | 99久久精品一区二区三区| 极品教师在线免费播放| 老司机在亚洲福利影院| 1024手机看黄色片| 国产av在哪里看| 亚洲avbb在线观看| 黄片小视频在线播放| 在线国产一区二区在线| 天美传媒精品一区二区| 国产亚洲精品久久久久久毛片| 亚洲成人久久爱视频| 精品电影一区二区在线| 特大巨黑吊av在线直播| 亚洲,欧美精品.| 国产精品亚洲美女久久久| 国产成人影院久久av| 大型黄色视频在线免费观看|