999精品在线视频,手机成人午夜在线视频,久久不卡国产精品无码,中日无码在线观看,成人av手机在线观看,日韩精品亚洲一区中文字幕,亚洲av无码人妻,四虎国产在线观看 ?

催化裂化裝置四旋分離系統(tǒng)內(nèi)氣相流場的數(shù)值研究

2014-09-16 03:52:04王建軍金有海
石油煉制與化工 2014年10期

王 銳,王建軍,趙 艷,金有海

(中國石油大學(xué)(華東)化工學(xué)院,山東 青島 266555)

催化裂化裝置四旋分離系統(tǒng)內(nèi)氣相流場的數(shù)值研究

王 銳,王建軍,趙 艷,金有海

(中國石油大學(xué)(華東)化工學(xué)院,山東 青島 266555)

采用雷諾應(yīng)力模型(RSM)對兩種催化裂化裝置四旋分離系統(tǒng)內(nèi)氣相流動的三維流場進行數(shù)值模擬,分析了兩種第四級旋風(fēng)分離系統(tǒng)(簡稱四旋)內(nèi)氣相流場的特點。結(jié)果表明:在傳統(tǒng)式結(jié)構(gòu)中,儲料罐內(nèi)部分區(qū)域氣流向上運動,不利于顆粒沉降,料腿內(nèi)產(chǎn)生氣流反竄現(xiàn)象,四旋內(nèi)存在偏流,影響顆粒輸送,結(jié)構(gòu)設(shè)計不合理是導(dǎo)致顆粒堵塞的主要原因;改進式結(jié)構(gòu)中,儲料罐內(nèi)氣相速度很小,有利于顆粒沉降,料腿中周期性的二次渦影響顆粒的氣力輸送,可能對壁面產(chǎn)生磨損,四旋分離空間內(nèi)流場穩(wěn)定,有利于顆粒分離,改進式結(jié)構(gòu)可以有效避免顆粒堵塞現(xiàn)象的發(fā)生。兩種結(jié)構(gòu)內(nèi)部氣相流場的對比為合理設(shè)計四旋分離系統(tǒng)結(jié)構(gòu)提供了參考依據(jù)。

催化裂化 四旋分離系統(tǒng) 氣相流場 數(shù)值模擬

煉油廠催化裂化裝置能量回收系統(tǒng)中,第三級旋風(fēng)分離器(簡稱“三旋”)是其關(guān)鍵設(shè)備之一[1]。三旋的作用是將再生器中前兩級旋風(fēng)分離器未能從再生煙氣中分離出來的催化劑微粒分離出來,為煙氣透平機提供凈化的煙氣,確保透平機長期穩(wěn)定運轉(zhuǎn)[2]。為了保證三旋高效運行,需要設(shè)置臨界噴嘴對其卸料口進行抽氣,但生產(chǎn)中發(fā)現(xiàn),臨界噴嘴磨損比較嚴(yán)重,因而許多煉油廠在三旋與臨界噴嘴間設(shè)置第四級旋風(fēng)分離系統(tǒng)(四旋),將三旋排出煙氣中的催化劑進一步分離后再進入臨界噴嘴。現(xiàn)階段對于四旋分離系統(tǒng)的研究較少,多數(shù)集中在對旋風(fēng)分離器的研究上。

目前,對于旋風(fēng)分離器內(nèi)氣相流場的模擬,多采用k-ε模型、代數(shù)應(yīng)力模型(ASM)、雷諾應(yīng)力模型(RSM)以及大渦模擬等。旋風(fēng)分離器內(nèi)氣流強旋轉(zhuǎn)運動造成流場的各向異性,而k-ε模型基于各向同性假設(shè),ASM考慮湍流各向異性,但對各向異性特征的描述能力有限,均不符合實際流動情況;大渦模擬雖然可以直接模擬各向異性湍流中的大渦,但計算量太大,較難實現(xiàn)[3];而RSM考慮各向異性,計算量適中,而且近年來許多專家利用不同的湍流模型對旋風(fēng)分離器內(nèi)的流場進行了模擬,其中RSM模型得到的結(jié)果與實驗結(jié)果吻合較好[4-6]。本課題采用RSM湍流模型對工程應(yīng)用中無翼閥的四旋分離系統(tǒng)內(nèi)氣相流場進行研究,為進一步優(yōu)化系統(tǒng)結(jié)構(gòu)、提高分離性能奠定基礎(chǔ)。

1 幾何與數(shù)學(xué)模型

1.1 幾何模型與網(wǎng)格劃分

傳統(tǒng)式與改進式四旋分離系統(tǒng)結(jié)構(gòu)尺寸如圖1所示。兩種結(jié)構(gòu)的不同點在于:傳統(tǒng)式結(jié)構(gòu)中三旋卸料口出氣先進入儲料罐,反轉(zhuǎn)后經(jīng)封頭上部管路進入四旋,即氣體先進儲料罐,再進四旋;改進式結(jié)構(gòu)中,三旋卸料口直接與四旋入口相連,即氣體先進四旋再進儲料罐,儲料罐頂部接平衡管,用于對四旋分離系統(tǒng)進行泄氣。儲料罐筒體段直徑為2 400 mm,錐段錐角為60°,上端接標(biāo)準(zhǔn)橢圓形封頭;旋風(fēng)分離器筒體段直徑為450 mm,進氣口尺寸為288 mm×116 mm,模型以旋風(fēng)分離器環(huán)形空間的頂蓋中心為坐標(biāo)原點,Z軸沿軸向向下為正。

四旋分離系統(tǒng)網(wǎng)格如圖2所示。為了提高網(wǎng)格質(zhì)量,對分離系統(tǒng)采取了分區(qū)生成網(wǎng)格的方法,并對不同區(qū)域分別采用結(jié)構(gòu)化與非結(jié)構(gòu)化網(wǎng)格進行劃分,網(wǎng)格總數(shù)分別為1 110 000和1 100 000。

1.2 數(shù)學(xué)模型和邊界條件

Fluent軟件的RSM模型中,雷諾應(yīng)力各分量的輸運方程為:

圖1 傳統(tǒng)式與改進式四旋分離系統(tǒng)結(jié)構(gòu)(單位:mm)

圖2 四旋分離系統(tǒng)網(wǎng)格

式中:ρ為流體密度,kgm3;t為時間,s;ui,uj,uk分別為速度在各方向(i,j,k)上的分量,ms;x=(xi,xj,xk),為笛卡爾坐標(biāo),m;Di,j為擴散項;Pi,j為應(yīng)力產(chǎn)生項;Gi,j為浮力產(chǎn)生項;Фi,j為壓力應(yīng)變再分配項;εi,j為耗散項;Fi,j為旋轉(zhuǎn)系統(tǒng)產(chǎn)生項;Suser為自定義源項。

采用RSM模型模擬旋風(fēng)分離器內(nèi)非穩(wěn)態(tài)不可壓縮湍流流動,使用有限體積法建立離散方程,求解控制方程時采用QUICK差分格式和SIMPLE算法[7-8]。

入口邊界條件:入口氣流為實際工況條件下的空氣,密度為2.922 kg/m3,黏度為3.028 2×10-5kg/(m·s);認(rèn)為來流氣體的湍流已經(jīng)充分發(fā)展,直接給出入口截面法向處的時均速度,大小為18 m/s;改進式四旋分離系統(tǒng)中,泄氣口處氣體速度為0.9 m/s,即泄氣率為5%。

出口邊界條件:按充分發(fā)展管流條件處理,除壓力外所有變量在出口截面法向上的梯度為零。為保證充分發(fā)展條件成立,在計算中將旋風(fēng)分離器排氣管路加長。

壁面邊界條件:采用無滑移邊界條件,對近壁網(wǎng)格點用壁面函數(shù)近似處理。

2 傳統(tǒng)式四旋分離系統(tǒng)內(nèi)氣相流場分布

目前,在催化裂化裝置中主要采用的是傳統(tǒng)式四旋分離系統(tǒng)結(jié)構(gòu),但在現(xiàn)場應(yīng)用中發(fā)現(xiàn),四旋下部錐段出現(xiàn)了催化劑顆粒堵塞的現(xiàn)象,高溫催化劑顆粒長時間堆積在分離器內(nèi)使其壁面呈現(xiàn)紅色。對傳統(tǒng)式四旋分離系統(tǒng)內(nèi)氣相流場分布進行數(shù)值研究,有助于找到顆粒堵塞的原因、發(fā)現(xiàn)結(jié)構(gòu)設(shè)計上的不足,為結(jié)構(gòu)優(yōu)化提供指導(dǎo)。

2.1 儲料罐內(nèi)流場分布

儲料罐內(nèi)流體速度分布云圖如圖3所示。由圖3可知,切向速度分布中,橢圓形封頭及錐形擋板內(nèi)局部區(qū)域速度達到2.0 m/s,使催化劑顆粒在該空間內(nèi)旋轉(zhuǎn)產(chǎn)生離心力,大粒徑顆粒被甩到邊壁后沿錐形擋板滑落到儲料罐底部,對顆粒起到預(yù)分離作用,而儲料罐底部灰斗空間內(nèi)切向速度幾乎為零,這樣可以防止被分離顆粒卷起并再次進入四旋而增加其負(fù)荷。軸向速度分布中,氣流由三旋卸料口以18.0 m/s的速度進入到儲料罐內(nèi),由于流通截面積變大,導(dǎo)致氣體進入罐內(nèi)后很快擴散,軸向速度急劇減小,到錐形擋板下部減小到4.2 m/s,攜帶催化劑顆粒的氣體向下運動至底部灰斗空間時速度降至0.6 m/s,有利于顆粒在此沉降。由圖3還可以看出,錐形擋板上部局部區(qū)域軸向速度出現(xiàn)負(fù)值,即氣相攜帶顆粒以1.2 m/s的軸向速度向上運動,對顆粒的沉降是不利的。

圖3 傳統(tǒng)式四旋分離系統(tǒng)儲料罐內(nèi)速度分布云圖

2.2 料腿內(nèi)流場分布

儲料罐內(nèi)靜壓分布如圖4所示。由圖4可知,灰斗空間內(nèi)壓力為1 050 Pa,四旋進氣管內(nèi)壓力為450 Pa,料腿內(nèi)部壓力為-600 Pa,因此,儲料罐內(nèi)壓力最高,與四旋進氣管和料腿間的壓差分別為600 Pa和1 650 Pa,儲料罐入口來流在壓差的作用下會分別進入四旋進氣管和料腿中。模擬結(jié)果顯示,進入料腿中的氣體質(zhì)量流量為0.486 kg/s,占入口總流量(1.757 kg/s)的27.7%,如此大的氣流量進入料腿中勢必會對料腿上部分離空間內(nèi)流場產(chǎn)生影響,四旋分離出的催化劑顆粒也無法排出至儲料罐內(nèi)。

圖4 傳統(tǒng)式四旋分離系統(tǒng)儲料罐內(nèi)靜壓分布云圖

料腿內(nèi)流線圖如圖5所示。由圖5可知,流線均向上,表明料腿內(nèi)產(chǎn)生氣流反竄現(xiàn)象,如果這部分氣體攜帶了已經(jīng)被分離出的顆粒,便會造成顆粒返混,即使四旋分離出的催化劑顆粒也無法由料腿排出。

圖5 傳統(tǒng)式四旋分離系統(tǒng)料腿內(nèi)流線圖

2.3 四旋分離空間內(nèi)流場分布

四旋分離空間內(nèi)速度沿徑向的分布如圖6所示。由圖6可知,切向速度分布呈現(xiàn)對稱性,在Z(距四旋頂部的距離)為350 mm截面處切向速度沿半徑呈現(xiàn)先增大后減小的趨勢,最大切向速度約為25 m/s,表現(xiàn)出“駝峰”狀分布特點,即準(zhǔn)強制渦和準(zhǔn)自由渦相結(jié)合的蘭金組合渦特征;隨著Z的增大,準(zhǔn)自由渦逐漸消失,在Z為650 mm截面處只表現(xiàn)出準(zhǔn)強制渦的特征,這是由于料腿中氣流上竄至分離空間,對其下部流場產(chǎn)生了干擾。軸向速度分布中,不同區(qū)域軸向速度分布分別表現(xiàn)出上行流與下行流特點,靠近一側(cè)邊壁附近為下行流,其它區(qū)域為上行流,即出現(xiàn)嚴(yán)重的偏流現(xiàn)象。上行流一側(cè)邊壁處的氣流無法攜帶催化劑顆粒向下運動,使顆粒無法排出。

圖6 傳統(tǒng)式四旋分離器分離空間內(nèi)速度沿徑向的分布Z值: ■—350 mm; ▲—500 mm; ●—650 mm

由以上分析可知,傳統(tǒng)式四旋分離系統(tǒng)之所以出現(xiàn)堵塞現(xiàn)象,是因為其結(jié)構(gòu)設(shè)計不合理。由于四旋入口進氣管和料腿均深入到儲料罐內(nèi),入口來流會分別流入四旋進氣管和料腿中,在料腿中產(chǎn)生氣流反竄現(xiàn)象,對四旋分離空間流場產(chǎn)生干擾;同時,四旋進氣管與料腿出口“相通”,其壓差太小,不足以推動四旋中被分離出的催化劑順利排出,在兩種因素的共同作用下便產(chǎn)生了顆粒的堵塞現(xiàn)象。顆粒的堵塞使四旋分離器失效,影響三旋的泄氣,造成排料不暢,進而影響其分離性能,因此需要對四旋分離系統(tǒng)結(jié)構(gòu)進行改造。

3 改進式四旋分離系統(tǒng)內(nèi)氣相流場分布

為了改善傳統(tǒng)式四旋分離系統(tǒng)結(jié)構(gòu)上的不足,提出了以下改進式結(jié)構(gòu):三旋卸料口來流直接進入四旋進行凈化,凈化后的氣體由排氣管排出,被分離顆粒由料腿流入儲料罐內(nèi),儲料罐頂部出口與排氣管均連接到臨界噴嘴。

3.1 儲料罐內(nèi)流場分布

儲料罐內(nèi)速度分布云圖如圖7所示。由圖7可知,切向速度分布中,相比于傳統(tǒng)式結(jié)構(gòu),改進式結(jié)構(gòu)速度小得多,且分布均勻,不存在局部速度較大的區(qū)域,對顆粒的沉降是有利的。軸向速度分布中,改進式結(jié)構(gòu)速度更小,儲料罐內(nèi)大部分區(qū)域速度在0左右,只是在灰斗內(nèi)靠近料腿處速度較大,約為0.5 m/s,有利于料腿中催化劑顆粒的氣力輸送。儲料罐頂部出口附近軸向速度出現(xiàn)負(fù)值,這是由于出口與臨界噴嘴相連,料腿中排出的氣體在壓差作用下要經(jīng)此處流出,而流通截面積的減小導(dǎo)致此處軸向速度增大。

圖7 改進式四旋分離系統(tǒng)儲料罐內(nèi)速度分布云圖

3.2 料腿內(nèi)流場分布

儲料罐內(nèi)靜壓分布如圖8所示。由圖8可知,相比于傳統(tǒng)式結(jié)構(gòu),改進式結(jié)構(gòu)中料腿內(nèi)外壓差小得多,在200 Pa以內(nèi),四旋分離器中被分離出的催化劑顆粒進入料腿后,能夠跟隨氣流在慣性作用下流出料腿并進入儲料罐內(nèi)進行沉降。

圖8 改進式四旋分離系統(tǒng)儲料罐內(nèi)靜壓分布云圖

料腿內(nèi)流線圖如圖9所示。由圖9可知,改進式結(jié)構(gòu)料腿內(nèi)部沒有產(chǎn)生氣流反竄現(xiàn)象,但出現(xiàn)了周期性的二次渦,使得四旋分離出的顆粒隨氣體在此處盤旋無法排出,顆粒的長時間停留還會對料腿壁面產(chǎn)生磨損。

圖9 改進式四旋分離系統(tǒng)料腿內(nèi)流線圖

3.3 四旋分離空間內(nèi)流場分布

四旋分離空間內(nèi)速度沿徑向分布如圖10所示。由圖10可知,與傳統(tǒng)式結(jié)構(gòu)一樣,切向速度分布呈現(xiàn)對稱性,沿軸向切向速度稍有減小,表明氣流旋轉(zhuǎn)強度變?nèi)酰啾葌鹘y(tǒng)式結(jié)構(gòu)切向速度更大,最大值為40 m/s左右,同樣出現(xiàn)在Z為350 mm截面,同時偏心現(xiàn)象減弱。因此,改進式結(jié)構(gòu)中氣流旋轉(zhuǎn)更加強烈,使催化劑顆粒獲得更大的離心力,有利于分離過程的進行。軸向速度分布中,同樣出現(xiàn)明顯的上行流和下行流。中心附近由于受到排氣管的影響而存在滯留區(qū),速度較小,隨半徑的增大,軸向速度先增大后減小;沿軸向排氣管的影響減弱,中心附近速度增大,遠(yuǎn)離中心位置速度減小。Z為350 mm截面處一側(cè)邊壁附近出現(xiàn)小范圍上行流是因為此處流場存在二次渦流。因此,改進式結(jié)構(gòu)中氣流攜帶催化劑顆粒邊旋轉(zhuǎn)邊向下運動,到分離器底部時反轉(zhuǎn)向上,顆粒在慣性的作用下繼續(xù)向下流至灰斗及料腿中,最終沉降到儲料罐底部。

圖10 改進式四旋分離器分離空間內(nèi)速度沿徑向的分布Z值: ■—350 mm; ▲—500 mm; ●—650 mm

由以上分析可知:在改進式結(jié)構(gòu)中,由于來流直接進入四旋分離器,且儲料罐頂部出口考慮了泄氣問題,料腿中沒有出現(xiàn)氣流反竄現(xiàn)象;四旋分離器內(nèi)流場穩(wěn)定,有利于顆粒分離與排出。因此,改進式結(jié)構(gòu)中不會出現(xiàn)顆粒堵塞現(xiàn)象,但料腿內(nèi)存在周期性二次渦,催化劑顆粒在此長時間停留會對壁面產(chǎn)生磨損。

4 結(jié) 論

(1) 傳統(tǒng)式四旋分離系統(tǒng)儲料罐內(nèi)氣相切向速度分布對催化劑顆粒有預(yù)分離作用,部分區(qū)域氣流向上運動,對顆粒沉降不利;料腿內(nèi)外壓差較大,產(chǎn)生氣流反竄現(xiàn)象;四旋內(nèi)存在偏流,影響顆粒輸運。結(jié)構(gòu)設(shè)計不合理是導(dǎo)致顆粒堵塞的主要原因。

(2) 改進式結(jié)構(gòu)中儲料罐內(nèi)氣相速度很小,有利于顆粒沉降;料腿中周期性的二次渦,影響催化劑顆粒的氣力輸送,并且可能對壁面產(chǎn)生磨損;四旋分離空間內(nèi)流場穩(wěn)定,有利于顆粒分離。改進式結(jié)構(gòu)可以有效避免顆粒堵塞現(xiàn)象的發(fā)生。

(3) Fluent軟件中的RSM模型可以很好地模擬出無翼閥的四旋分離系統(tǒng)內(nèi)部氣相流場分布,解釋傳統(tǒng)式結(jié)構(gòu)中顆粒堵塞的原因,可為結(jié)構(gòu)優(yōu)化提供指導(dǎo)。

[1] 黃榮臻,閆濤,房家貴.催化裂化第三級旋風(fēng)分離器的現(xiàn)狀和發(fā)展方向[J].石油化工設(shè)備技術(shù),2005,26(1):29-31

[2] 陳秉林.三級旋風(fēng)分離器[J].通風(fēng)除塵,1984(1):7-8

[3] 王樂勤,郝宗睿,王循明,等.筒體長度對旋風(fēng)分離器內(nèi)流場影響的數(shù)值模擬[J].工程熱物理學(xué)報,2009,30(2):223-226

[4] Hu Liyuan,Zhou Lixing,Zhang Jian,et al.Studies on strongly swirling flows in the full space of a volute cyclone separator[J].AIChE Journal,2005,51(3):740-749

[5] 王宏偉,王建軍,金有海.導(dǎo)葉式旋風(fēng)管入口環(huán)形空間內(nèi)氣相流場數(shù)值模擬[J].石油化工設(shè)備,2006,35(2):33-36

[6] 魏新利,張海紅,王定標(biāo),等.旋風(fēng)分離器內(nèi)顆粒軌跡的數(shù)值模擬[J].鄭州大學(xué)學(xué)報(工學(xué)版),2004,25(3):14-17

[7] 金有海,姬廣勤,曹晴云,等.旋風(fēng)分離器排氣管內(nèi)氣相流場的數(shù)值模擬[J].中國石油大學(xué)學(xué)報(自然科學(xué)版),2008,32(6):109-112

[8] 曹晴云,姬廣勤,金有海,等.不同結(jié)構(gòu)參數(shù)下旋風(fēng)分離器氣相流場的數(shù)值研究[J].流體機械,2008,36(6):34-38

NUMERICAL SIMULATION OF GAS-PHASE FLOW FIELD IN FOURTH STAGE SEPARATION SYSTEM IN FCCU

Wang Rui, Wang Jianjun, Zhao Yan, Jin Youhai

(CollegeofChemicalEngineeringofChinaUniversityofPetroleum(Huadong),Qingdao,Shandong266555)

Based on the Reynolds Stress Model, the three-dimensional gas-phase flow fields of two kinds of the fourth stage separation systems in FCCU were simulated. The characteristics of gas flow fields in the two systems were especially analyzed. The results show that in the material storage tank of the traditional structure, the gas flows upward at some areas, which is not conducive to particle sedimentation. Reverse flow exists in the dipleg. Bias flow happens in the cyclone separator, by which particle transportation is affected. The unreasonable structure design is the main reason for particle plugging. However, in the improved one, the air flow is relatively slow, favorable to particle sedimentation. In the dipleg, the pneumatic transportation is influenced by the periodic secondary vortex that can wear the wall. The flow field in the cyclone separator is stable enough for particle separation. Particle plugging can be avoided effectively in the improved structure. The comparison of the gas flow field in the two structures provides valuable guidance for designing efficient fourth stage separation system.

FCCU; the fourth stage separation system; gas flow field; numerical simulation

2014-03-10; 修改稿收到日期: 2014-06-17。

王銳,碩士研究生,從事多相流分離方面的研究工作。

王建軍,wangjianjun0992@126.com。

主站蜘蛛池模板: 四虎精品黑人视频| 国产对白刺激真实精品91| 国产精品一区二区久久精品无码| 激情影院内射美女| 成人福利一区二区视频在线| 欧美成人手机在线观看网址| 中文字幕丝袜一区二区| 在线国产91| 99久久精品免费观看国产| 91午夜福利在线观看| 任我操在线视频| 美女被狂躁www在线观看| 日韩成人高清无码| 91无码网站| 在线网站18禁| 五月天久久综合| 国产福利在线免费| 亚洲av片在线免费观看| 国产呦精品一区二区三区下载 | 亚洲精品爱草草视频在线| 日韩福利在线视频| 无码 在线 在线| 伊人成人在线| 97视频在线观看免费视频| 992tv国产人成在线观看| 国产精品亚洲欧美日韩久久| 久久福利网| 丁香亚洲综合五月天婷婷| 久久黄色免费电影| 一级看片免费视频| 在线国产资源| 婷婷久久综合九色综合88| 无码免费试看| 国产波多野结衣中文在线播放| 国产又粗又爽视频| 午夜限制老子影院888| 国产精品成人第一区| 国产午夜福利亚洲第一| 亚洲精品无码av中文字幕| 麻豆AV网站免费进入| 亚洲成人动漫在线| 色哟哟国产精品一区二区| 少妇极品熟妇人妻专区视频| 国产免费精彩视频| 自慰高潮喷白浆在线观看| 国产高清又黄又嫩的免费视频网站| 国产三区二区| 精品久久香蕉国产线看观看gif| 国产喷水视频| 毛片久久网站小视频| 中文一级毛片| 综合成人国产| 国产二级毛片| 五月天久久综合| 欧美精品亚洲二区| 99偷拍视频精品一区二区| 久久国产精品娇妻素人| 色成人亚洲| 91精品日韩人妻无码久久| 在线观看视频99| 三上悠亚精品二区在线观看| 亚洲视频一区在线| 久草青青在线视频| 日本高清成本人视频一区| 成人国产一区二区三区| 1024你懂的国产精品| 亚洲国产av无码综合原创国产| 综合色天天| 亚洲欧美日韩中文字幕一区二区三区| 色悠久久久久久久综合网伊人| 久热这里只有精品6| 乱人伦中文视频在线观看免费| 91口爆吞精国产对白第三集| 国产乱子伦精品视频| 久久伊人操| 日韩av电影一区二区三区四区| 国产免费羞羞视频| 免费一级毛片完整版在线看| 美女一区二区在线观看| 亚洲精品欧美日韩在线| 日韩精品毛片| 亚洲中久无码永久在线观看软件|