密閉熱流體研究報告401

1、項目驗收材料井筒密閉熱流體循環(huán)降粘工藝研究研 究 報 告2012年3月目 錄第一章 前言1第二章 井筒密閉熱流體循環(huán)工藝分析52.1井筒中循環(huán)水、產(chǎn)出液流動的溫度場計算52.2井筒中循環(huán)水、產(chǎn)出液流動的壓力場計算112.2.1井筒氣液兩相管流模型112.2.2產(chǎn)出液在油管內(nèi)的兩相流動122.2.3熱水在隔熱管、抽油桿環(huán)空內(nèi)的流動132.2.4井筒壓力場的計算132.3水的物性參數(shù)求解172.3.1水的熱力學性質(zhì)172.3.2未飽和水的熱遷移性質(zhì)20第三章 抽油工藝分析213.1抽油桿柱的載荷計算213.1.1懸點靜載荷213.1.2懸點慣性載荷223.1.3懸點振動載荷223.1.4摩擦載荷2

2、33.2懸點最大和最小載荷233.3抽油機扭矩23第四章 計算程序開發(fā)與設(shè)計254.1 軟件的研制開發(fā)目的254.2 軟件需求分析與功能說明254.3 開發(fā)環(huán)境的選擇254.4 軟件設(shè)計264.5 軟件的使用264.5.1軟件的安裝264.5.2進入系統(tǒng)264.5.3降粘計算過程264.5.4操作274.5.5退出274.6 軟件特點274.7 軟件的界面27第五章 計算實例分析315.1油井基本數(shù)據(jù)315.2循環(huán)熱水流量對降粘工藝的影響315.3循環(huán)熱水溫度對降粘工藝的影響35第六章 結(jié)論與建議39參考文獻40第一章 前言注蒸汽開采稠油是一種行之有效的方法，已成為國內(nèi)外廣泛采用的開采稠油、超

3、稠油的常規(guī)方法。世界上稠油資源極為豐富，隨著石油科技的日益發(fā)展，世界高凝油藏、稠油油藏不斷被發(fā)現(xiàn)，全世界稠油的總資源量約為已經(jīng)探明的常規(guī)原油儲量的6倍，約為15500108t。稠油油藏的分布范圍十分廣泛，世界上各產(chǎn)油國幾乎都有稠油資源，目前發(fā)現(xiàn)的稠油資源主要集中在美國、加拿大、委內(nèi)瑞拉和中國。我國的稠油資源也非常豐富，預(yù)計儲量在80108t以上，主要分布在松遼、渤海灣、準噶爾等盆地。稠油產(chǎn)量已成為我國原油產(chǎn)量構(gòu)成的重要組成部分，且逐年增長。一些稠油區(qū)塊的油質(zhì)較稠，雖然在油藏條件下原油具有一定的流動性和較好的供油能力，但原油進入井筒向地面流動的過程中，受井筒溫度、原油含水等因素的影響，粘度發(fā)生急

4、劇變化，逐漸失去流動性部分稠油井或超稠油井依靠地層的能量不能將原油舉升到地面，無法實現(xiàn)油井配產(chǎn)要求。為解決上述問題，必須加強井筒降粘工藝技術(shù)研究，目前國內(nèi)外在稠油開采過程中常用的井筒降粘工藝主要有：電加熱、熱流體循環(huán)、摻稀油及摻化學劑乳化降粘工藝等。井筒電加熱降粘工藝由于耗電量大，運行成本和初始投資高昂，降粘效果較好，但對于含水量大、礦化度高、超稠油降粘效果不理想；井筒化學降粘工藝地面流程相對復(fù)雜，運行成本較高，而且生產(chǎn)的原油分離工藝復(fù)雜，化學降粘劑需要根據(jù)不同油井原油物性進行篩選。井筒摻稀降粘工藝則利用區(qū)塊自己生產(chǎn)的稀油進行，成本低，降粘效果好，但地面生產(chǎn)工藝相對復(fù)雜，當在油井含水達到一定程

5、度，產(chǎn)出液出現(xiàn)反向乳化的情況下，降粘效果較差。熱流體循環(huán)工藝主要依靠稠油的強熱敏感性，適用于油層較淺、高凝固點、高粘度稠油油藏。熱流體循環(huán)工藝的特點：（1）工藝原理簡單，一次投資大；（2）循環(huán)熱水溫度應(yīng)保證產(chǎn)液溫度不低于粘溫曲線拐點值，防止由于溫度太低，溫度突變（可出現(xiàn)對摻入水加熱情況）；（3）井筒熱流體循環(huán)工藝是針對稠油依賴于溫度的強熱敏感性；（4）井筒熱流體循環(huán)工藝基本不受井深的限制，它除了提高產(chǎn)液的溫度外，還可以通過提高井筒中混合液（產(chǎn)液+摻入的熱流體）的含水量來降低粘度。常用的熱流體有熱水、水蒸氣、柴油等。常用的熱流體循環(huán)可分為井筒熱流體循環(huán)和空心抽油桿熱流體循環(huán)；井筒熱流體循環(huán)按其循

6、環(huán)方式又可分為開式循環(huán)、閉式循環(huán)。開式循環(huán)是指循環(huán)液從油套環(huán)空進入，在井下某一深度與產(chǎn)出液混合后一起被舉升到地面，開式循環(huán)又分為正循環(huán)和反循環(huán)，見圖1。閉式循環(huán)又稱分程式熱流體循環(huán)，見圖2。與開式熱流體循環(huán)不同的是：循環(huán)熱流體與從油層采出的原油不相摻混，因而解決了由于熱流體施加于油層壓力而導致的油井舉升困難問題。閉式熱流體循環(huán)中，熱流體所攜帶的熱量是通過油管壁傳遞給油管中的油氣混合物，使之升溫而達到防蠟、降粘的目的。采用的熱流體可以是柴油、水或其他流體，當然也可以為原油，但要求這種原油含蠟少、凝固點低、粘度不甚高等特點，以保證在關(guān)井后便于再啟動。空心抽油桿熱流體循環(huán)：該工藝的熱效率大大高于應(yīng)用

7、油管循環(huán)的效率，對于平均溫度可以提高40左右。其原因是循環(huán)液始終與套管接觸，減少了熱傳遞。 （a）反循環(huán) （b）正循環(huán)圖1 開式熱流體循環(huán)工藝圖 2 常規(guī)閉式熱流體循環(huán)工藝為了進一步提高產(chǎn)出液的溫度，降低熱量損失，提出一種新型閉式熱流體循環(huán)工藝，見圖3。該工藝采用空心抽油桿內(nèi)加一隔熱管，熱流體從隔熱管內(nèi)注入，通過隔熱管和抽油桿間環(huán)空循環(huán)放熱后流出；油管和抽油桿間是產(chǎn)出液的流道，保持油套環(huán)空。該循環(huán)工藝的優(yōu)勢是：充分利用循環(huán)熱水攜帶的熱量以及油套環(huán)空的隔熱效果，提高產(chǎn)出液的出口溫度，使產(chǎn)出液在井筒中維持較高的溫度，達到井筒降粘的目的。油層油層253mm隔熱管101.6mm油 管套 管抽油泵熱水產(chǎn)

8、出液505mm空心桿圖3 新型熱水閉式循環(huán)降粘舉升工藝示意圖封割器因此，對稠油組分、油品性質(zhì)、流變性、粘溫關(guān)系等進行綜合分析的基礎(chǔ)上，對比分析前期井筒降粘工藝，根據(jù)產(chǎn)液性質(zhì)、油品性質(zhì)、含水率等優(yōu)化閉式熱水循環(huán)降粘工藝參數(shù)，從而進一步提高系統(tǒng)效率，有效指導油氣生產(chǎn)，為形成完善的配套工藝技術(shù)提供技術(shù)支撐。第二章 井筒密閉熱流體循環(huán)工藝分析整個閉式循環(huán)降粘工藝中，熱水在隔熱管內(nèi)向下流動，同時將一部分熱量傳遞給隔熱管和空心抽油桿環(huán)空中向上流動的循環(huán)排出水；同時排出水又與空心抽油桿和油管環(huán)空中向上流動的產(chǎn)出液進行熱量的交換；而產(chǎn)出液又通過油套環(huán)空、水泥環(huán)將熱量傳遞給地層，整個過程伴隨著動量、熱量的傳遞。

9、2.1井筒中循環(huán)水、產(chǎn)出液流動的溫度場計算根據(jù)合理的假設(shè)將上所述管柱結(jié)構(gòu)示意圖簡化為圖4所示的物理模型。由于注入時采用隔熱管，這樣將熱水攜帶的熱量絕大部分帶入底部，通過往上的流動與產(chǎn)出液進行熱量交換，從而確保產(chǎn)出液的溫度不至于由于散熱降低。同時由于保存了油套環(huán)空，從而使得產(chǎn)出液散熱的速度有所減緩，也有利于維持產(chǎn)出液的溫度。圖 4 新型閉式熱流體循環(huán)工藝熱水在井筒隔熱管內(nèi)自上往下流動，到空心桿底部后從空心桿與隔熱管環(huán)空中返回地面，產(chǎn)出液則在空心桿與油管環(huán)空中自下向上流動，油套管之間的環(huán)空中可能是空氣、水或油水混合物。由于水泥環(huán)外側(cè)的地層體積和熱容量比水泥環(huán)內(nèi)側(cè)的井筒大得多，水泥環(huán)與地層交界處的溫

10、度隨熱水循環(huán)時間的推移緩慢上升，且溫度與未受熱影響遠處地層溫度的差值不大。相反，循環(huán)熱水溫度與地層溫度的差值較大。因而可以假定水泥環(huán)內(nèi)側(cè)(以下簡稱“井筒內(nèi)”)熱量傳遞是穩(wěn)態(tài)的，而水泥環(huán)外側(cè)(以下簡稱“地層內(nèi)”)熱量傳遞則按非穩(wěn)態(tài)處理。井筒熱量傳遞方式比較復(fù)雜，自內(nèi)向外有：注入熱水與隔熱管內(nèi)壁之間的受迫對流換熱，隔熱管的導熱，隔熱管外側(cè)和空心桿內(nèi)側(cè)環(huán)空的對流換熱，空心桿的導熱，空心桿外側(cè)和油管內(nèi)側(cè)的對流換熱，油管的導熱，油套管環(huán)空中對流輻射同時作用的換熱以及水泥環(huán)的導熱等環(huán)節(jié)，這些熱阻徑向串并聯(lián)組成“井筒內(nèi)”的傳熱。如果把這一傳熱過程按非穩(wěn)態(tài)處理，求解是非常困難的。井筒內(nèi)穩(wěn)態(tài)傳熱的合理假設(shè)把這一

11、時間因素劃歸到井筒外側(cè)地層的導熱中。后者屬圓柱坐標系中非穩(wěn)態(tài)純導熱問題。如果把井筒長度分成足夠短的若干段，不計及每段長度方向的熱傳遞，則為徑向一維非穩(wěn)態(tài)導熱的典型命題，可用解析法或數(shù)值法求解。井筒內(nèi)穩(wěn)定傳熱的假設(shè)并非意味井筒內(nèi)各處溫度(包括)不隨時間變化。由地層非穩(wěn)態(tài)導熱得出的隨時間變化的溫度，再反饋到井筒內(nèi)穩(wěn)態(tài)傳熱公式的溫差(-)中去，間接地反映了井筒內(nèi)傳熱的非穩(wěn)態(tài)性質(zhì)。因此有人把井筒內(nèi)的傳熱稱為“擬穩(wěn)態(tài)”。熱流體在井筒流動過程中，井筒中徑向熱流量，即由隔熱管徑向流向井筒周圍的熱流量，就是井筒熱損失。2.1.1假設(shè)條件井口產(chǎn)出液的壓力、溫度保持不變，動液面在一定時間內(nèi)保持不變；油管與套管形成

12、的環(huán)形空間充滿低壓空氣；根據(jù)ramey和satter的方法，將井筒的徑向傳熱看作是由油管中心到水泥環(huán)外緣的一維穩(wěn)態(tài)傳熱和水泥環(huán)外緣到地層之間的一維不穩(wěn)態(tài)傳熱；忽略地層導熱系數(shù)沿井深方向的變化，并視為一個常數(shù)。2.1.2井筒熱傳遞數(shù)學模型根據(jù)圖4，熱水從隔熱管中注入，從隔熱管和空心抽油桿之間的環(huán)空返回；原油從空心抽油桿和油管環(huán)空流出，他們之間進行了熱量的傳遞。沿井深方向取微元長度，根據(jù)能量平衡可得： (1)式中，為隔熱管內(nèi)外流體間、油管內(nèi)外流體間、環(huán)形空間流體與地層間的傳熱系數(shù)，w/(m)；為循環(huán)熱水的質(zhì)量流量，kg/s；為產(chǎn)出液的質(zhì)量流量，kg/s；為隔熱管內(nèi)熱水的溫度，；為空心抽油桿內(nèi)循環(huán)熱

13、水回水溫度，；為油管內(nèi)內(nèi)產(chǎn)出液的溫度，；、分別為熱水比熱、產(chǎn)出液當量比熱，j/(kg)；為熱水向循環(huán)回水的傳熱量，w；為循環(huán)回水向產(chǎn)出液的散熱量，w；為產(chǎn)出液向地層的散熱量，w；為井深，m；為初始地層溫度。，k；為地表溫度，k；為地溫梯度，k/m；z為井深，m；。2.1.3油管中心至水泥環(huán)外緣的傳熱由穩(wěn)態(tài)傳熱公式得： （2）式中，為傳熱系數(shù)， ；r為總傳熱熱阻；為注入熱水溫度，k；為水泥環(huán)外緣初溫度，k；為井筒長度，m；為單位時間內(nèi)長度上的熱損失，kw；為隔熱管內(nèi)徑，m。熱阻包括以下九部分：（1）隔熱管內(nèi)側(cè)熱水流換熱熱阻： （3）式中，為熱水對流換熱系數(shù)，kw/（m2k）。（2）隔熱管的導熱熱

14、阻： （4）式中，為隔熱管導熱系數(shù)，kw/（mk）；為隔熱管外半徑，m。（3）隔熱管與空心抽油桿環(huán)空的循環(huán)回水對流換熱熱阻： （5）式中，為循環(huán)熱水回水對流換熱系數(shù)，kw/（m2k）。（4）空心抽油桿的導熱熱阻： （6）式中，為空心抽油桿導熱系數(shù)，kw/（mk）；為空心抽油桿內(nèi)半徑，m；為空心抽油桿外半徑，m。（5）空心抽油桿與油管環(huán)空產(chǎn)出液對流換熱熱阻： （7）式中，為產(chǎn)出液在桿管環(huán)空中的對流換熱系數(shù)，kw/（m2k）。（6）油管的導熱熱阻： （8）式中，為油管導熱系數(shù)，kw/（mk）；為油管內(nèi)半徑，m；為油管外半徑，m。（7）環(huán)空的自然對流和輻射換熱的熱阻： （9）式中，為環(huán)空內(nèi)自然對流換

15、熱系數(shù)，kw/(m2k)；為環(huán)空內(nèi)輻射換熱系數(shù)，kw/(m2k)。在計算井筒內(nèi)流體溫度分布、套管溫度及井筒熱損失率時，最關(guān)鍵的是如何確定在具體井筒機構(gòu)條件下總傳熱系數(shù)。而且，最困難的是如何準確計算出環(huán)空液體或氣體的熱對流、熱傳導及熱輻射都存在的條件下的環(huán)空傳熱系數(shù)。因為它與油管外表面性質(zhì)、液體的物理性質(zhì)（尤其是高溫下的粘度變化），油管外壁與套管內(nèi)壁之間的溫度與距離，套管內(nèi)壁表面性質(zhì)等等都有影響，計算很復(fù)雜。這里用一種傳統(tǒng)的計算方法。 確定輻射傳熱系數(shù)： 當油管環(huán)空或隔熱管與套管之間充有氣體時，輻射熱流量取決于注入管外壁溫度與套管內(nèi)壁溫度，按stefan-boltzmann定律，得： （10）星

16、號是絕對溫度（），是stefan-boltzmann常數(shù),。是由油管外壁表面向套管內(nèi)壁表面輻射散熱有效系數(shù)，它代表吸收輻射的能力。對于井筒條件： （11）式中，及是油管外壁及套管內(nèi)壁的輻射系數(shù)，無因次。是兩個表面間的總交換稀疏。對于井筒傳熱條件，平常取值為1.0。因此（11）式簡化為： （12） 由公式（10）的因子分解，可得出輻射傳熱系數(shù)的表達式為： （13） （14） 如已知及就可計算出。 確定自然對流傳熱系數(shù)：在油套環(huán)空間的熱傳導及自然對流引起的徑向熱流速度為： （15）式中，為熱傳導及熱對流引起的徑向熱流速度，w；為環(huán)空液體的等效導熱系數(shù)（也稱環(huán)空介質(zhì)當量導熱系數(shù)），即在環(huán)空平均溫度及

17、壓力下，包括自然對流影響的環(huán)空液體的綜合導熱系數(shù)，。當自然對流很小時，是環(huán)空液體或氣體的導熱系數(shù)。因為， （16） （17）的表達式為： （18） （19）式中，為環(huán)空介質(zhì)導熱系數(shù)，w/(m)；為瑞利準數(shù)，無因次；為環(huán)空介質(zhì)的熱膨脹系數(shù)，；為環(huán)空中介質(zhì)在平均溫度下的密度，kg/m2；為環(huán)空中介質(zhì)在平均溫度下的定壓比熱，j/(kg)。（8）套管的導熱熱阻： （20）式中，為套管的導熱系數(shù)，kw/(mk)；為套管外半徑，m；為套管內(nèi)半徑，m。（9）水泥環(huán)的導熱熱阻： （21）式中，為水泥環(huán)導熱系數(shù)，kw/（mk）；為水泥環(huán)外緣半徑，m。在這九項中， 和其值太小，可忽略不計，若以隔熱管外表面為基準面

18、積，總傳熱熱阻可寫為： （22）總傳熱系數(shù)為： （23） 各傳熱系數(shù)為： （24） （25） （26）2.1.4從水泥環(huán)外緣至地層的導熱 由于這是不穩(wěn)定的熱傳導，隨時間而變化，用公式可表示為： （27）式中，為地層導熱系數(shù)，j/（mk）；為無因次地層導熱時間函數(shù)。 我們可用hasan公式： （28）式中，為無因次變量，；這里，為（注汽或）生產(chǎn)時間, d；為地層平均熱擴散系數(shù), ，。式（27）中，為水泥環(huán)與地層交界面的溫度，；為平均原始地層溫度，；為隨時間變化的導熱傳熱函數(shù)。2.2井筒中循環(huán)水、產(chǎn)出液流動的壓力場計算井筒中注入的循環(huán)水為單相流，產(chǎn)出液則可能由于含油伴生氣呈現(xiàn)為兩相流。在沿井筒流動

19、過程中，由于熱量不斷向地層散失，介質(zhì)的壓力、溫度不斷變化，由此將影響到介質(zhì)的物性參數(shù)的變化，而物性參數(shù)的變化又會反過來影響壓力和溫度的變化。2.2.1井筒氣液兩相管流模型（一）按外形劃分兩相流模型由于氣液兩相流中氣相密度小于液相密度，因此它們在水平管內(nèi)流動、垂直向上流動、垂直向下流動、傾斜向上流動、傾斜向下流動等不同情況下，外形互不相同。在實驗室目測或高速攝影得到的圖象說明了這一現(xiàn)象。不同學者給出了多種流型的名稱，如泡狀流、塞狀流、分層流、波狀流、彈狀流、環(huán)狀流、團塊流、細束環(huán)狀流、霧狀流、下降膜流等等。這樣給計算兩相流壓降造成了困難。這說明按外形劃分兩相流流型不是一種好方法。（二）按流體力學

20、特征劃分兩相流流型綜觀上面按外形劃分的眾多流型，按它們流體力學的特征可以歸納為三種流型：1分散流（distributed flow）氣液兩相分布比較均勻2間歇流（intermittent flow）氣液兩相分布不規(guī)則、帶有波動3分離流（segregated flow）氣液兩相有較清晰較平穩(wěn)的分界面上面10幾種流型都可以歸并到這三種流型中，設(shè)法找出劃分這三種流型的一些判據(jù)，并且分別確定這三種流型的流動特性（如摩阻系數(shù)），兩相流壓降問題就能解決。最后一個問題通常是用實驗手段解決的。這類實驗一直在進行中，力求實驗越來越接近實際的兩相流過程。到目前為止，還沒有完全解決兩相流問題。（三）兩相流流型的幾種

21、主要計算模型有些工藝過程中遇到的兩相流只呈現(xiàn)分散流、間歇流或分離流中的某一種流型；另外一些工藝過程則可能在不同工況下遇到不同的流型，因此出現(xiàn)以下幾種不同的計算模型：1均相（流）模型（homogeneous flow model）把兩相流看成一種均勻相，此均勻相參數(shù)充分反映了液相和氣相的特征。2分相（流）模型（separated flow model）把兩相流看成氣液各自分開的流動。3流型模型（flow pattern model）先判斷不同流型，不同流型有不同的規(guī)律。由于產(chǎn)出液在輸汽管線、注汽井和油層中流動會出現(xiàn)不同的流型，因此用流型模型計算產(chǎn)出液的壓降比較合適。由于產(chǎn)出液舉升過程是氣液兩相流

22、動，同時在流動過程中熱量不斷向外界散失，產(chǎn)液溫度不斷下降，混合密度、粘度、表面張力等不斷變化，是流動多變而復(fù)雜，計算繁瑣。目前國內(nèi)外由許多相關(guān)的研究，認為beggs-brill方法（簡稱b-b法）比較精確。2.2.2產(chǎn)出液在油管內(nèi)的兩相流動井筒中產(chǎn)出液兩相流動時，溫度和壓力的變化會使氣液界面間的平衡打破，通過能量的交換后重新建立新的平衡。同時油、氣、水各相的密度、粘度等參數(shù)也跟隨變化?？梢?，產(chǎn)出液在井筒中的流動是多變而復(fù)雜的，計算繁瑣。兩相流體在井筒中流動時，取垂直向下為正方向，由能量方程可得： (29)式中，為井筒中產(chǎn)出液混合物壓力，pa；為井筒深度，m；為產(chǎn)出液的密度，kg/m3；為產(chǎn)出液

23、氣液混合物的速度，m/s；為能量損失；為所做的軸功，=0。將上式變形得， (30)式中，為摩阻壓降，定義為，pa/m；為重力壓降，定義為，pa/m；為加速壓降，定義為，一般較小，可忽略，pa/m；為兩相流摩擦系數(shù)，無量綱；為氣液兩相混合流速，m/s；、分別為氣相、液相的表觀流速，m/s；為井筒傾角，（）。因此，井筒內(nèi)兩相管流的壓力降是摩擦能量損失、重位壓降損失和加速壓降損失之和，壓力降公式可表示為=+ (31) 忽略加速壓降，上式變?yōu)椋?(32) 式中，為油管當量內(nèi)直徑，m。2.2.3熱水在隔熱管、抽油桿環(huán)空內(nèi)的流動熱水在隔熱管以及隔熱管和空心抽油桿環(huán)空中的流動為單相流，根據(jù)動量方程，忽略加速

24、壓降。由式（30）可得到熱水流動的壓力損失： (33) 式中，為熱水的密度，kg/m3；為熱水的流速，m/s；為隔熱管的內(nèi)直徑（環(huán)空時為當量直徑），m；為單相流的摩擦系數(shù)，可根據(jù)雷諾數(shù)判斷流態(tài)得到。2.2.4井筒壓力場的計算早期人們直接用層流湍流概念研究氣液兩相流動，著名的lockhart-martinnelli壓降計算方法便是以層流湍流組合為基礎(chǔ)。但是隨著實驗技術(shù)進展和研究深入，人們認識到不同的交界面構(gòu)形反映了不同的水力特性，流型變化意味著相交接界面形狀的變化，因而意味著相之間的動量傳遞模式和熱量傳遞模式變化。也即改變了控制這類傳遞諸效應(yīng)之間的相對權(quán)重。描述流型方法很多，最為簡單的方法是形態(tài)

25、學法，它按兩種相的相對形態(tài)確定流型的類別以及各流型相互轉(zhuǎn)化的過渡條件。一般說來，蒸發(fā)管內(nèi)出現(xiàn)的各種流型與具有同一當?shù)亓鲃訔l件下的絕熱流到內(nèi)的流型基本相同。若蒸法管受不太高均勻熱流密度加熱，只要流道足夠長，入口為近飽和水，會依次發(fā)生如圖所示的各種流型。即單相流體，泡狀流，彈狀流，環(huán)狀流，彌散流，單相氣體流動等。井筒內(nèi)油氣水兩相管流的壓力降是摩擦損失、勢能變化和動能變化的綜合結(jié)果，即垂直多相管流的壓力梯度，是三個壓降之和：靜水壓力梯度、耗于摩阻的壓力梯度以及耗于加速度的壓力梯度。由于注入井筒中的水蒸氣含有氣液兩相，同時由于熱量不斷向地層散失，水蒸氣的溫度不斷下降，其混合密度、粘度、表面張力等不斷變

26、化，是流動多變而復(fù)雜，計算繁瑣。目前國內(nèi)外由許多相關(guān)的研究，認為beggs-brill方法（簡稱b-b法）比較精確。b-b方法是流型模型中的一種計算兩相流特性的方法，是beggs博士論文的核心內(nèi)容，發(fā)表在1973年。beggs在一套長達27.4m的透明的聚丙烯管中進行的兩相流實驗，介質(zhì)為空氣和水，管徑有2.54cm()和3.81cm()兩種，管子可以傾斜，角度變化為，其中為氣液垂直向下，為氣液垂直向上，為氣液水平流動，總共作了584次實驗，空氣和水的流量可以分別變化，因此可觀察到各種流型。用b-b方法計算壓力分布的步驟： 從已知壓力出發(fā)，假定一個值，這兩個壓力點之間的距離為，算出摩阻系數(shù)，進而

27、可計算出，可得到另一個，將兩個值進行比較，如果不相等，將賦予，如此循環(huán)下去，直到相等為止，計算步驟如下：（1）算出在這二個壓力點之間平均壓力； 其中，井筒深度，m；為在處的壓力，mpa。（2）確定在平均深度處的水蒸氣平均溫度，該值可由溫度場得到；（3）算出在平均壓力和平均溫度條件下液相和氣相的密度、流體含氣率；（4）算出水蒸氣的混合密度（流量密度）、就地的氣體和液體的流量； 式中，為流量密度，kg/m3；為水蒸氣注入量，m3/s；為液相流量，m3/s；為氣相流量，m3/s。（5）算出就地的表觀氣、液流速； 式中，為液相表觀流速，m/s；為氣相表觀流速，m/s；為油管截面積，m2；為混合物表觀流

28、速，m/s。（6）算出單位截面積上的氣、液，以及總的質(zhì)量流量： 式中，為單位截面積上的氣相質(zhì)量流量， kg/m2s；為單位截面積上的液相質(zhì)量流量， kg/m2s；為單位截面積上總的質(zhì)量流量， kg/m2s。（7）算出入口的液體含量（體積含液率） 式中 ：體積含液率。（8）算出弗洛德準數(shù)、液相粘度、氣相粘度、混合粘度，以及液相的表面張力； 式中，為弗洛德準數(shù)，無因次；為液相粘度，mpas；為汽相粘度，mpas；為混合粘度，mpas。（9）算出無滑脫雷諾數(shù)及液相流速準數(shù)； ， （10）為確定流動為水平時，可能存在的流態(tài)譜式，要先算出各相關(guān)參數(shù)，和； （11）用下列范圍確定流態(tài)譜式； 分異型：0.0

29、1及nfr 或 0.01及nfr； 過渡型：0.01及nfr； 間隔型：0.010.4及nfr 或 0.4及nfr； 分散型： 0.4及nfr 或 0.4及nfr。（12）算出水平滯留量（水平截面持液率） 式中的，分別按下表每種流型而定流態(tài)譜式分異型0.980.48460.0868間隔型0.8450.53510.0173分散型1.0650.58240.0609（13）算出斜度校正系數(shù) 式中的，和分別按下表不同流動狀態(tài)選值流態(tài)譜式分異型，上坡0.0113.7683.5891.614間隔型，上坡2.960.3050.44730.0978分散型，上坡不必校正 各種流型的下坡：=4.70，=-0.36

30、92，=0.1244，=-0.5056。（14）算出液相滯留量的斜度校正系數(shù)，對垂向井來說： ；（15）算出液相滯留量及兩相密度 (16)算出摩阻系數(shù)比： 式中，當值處于1 1.2區(qū)間時變?yōu)椴欢ǎ?所以當處于這一區(qū)間時，函數(shù)改用下公式：。（17）算出無滑脫的摩阻系數(shù) 或 ；（18）算出兩相摩阻系數(shù) （19）算出 （20）算出 （21）將與相比較，如果，則停止運行計算；否則，將賦予，重新計算，直到滿足上式條件為止。2.3水的物性參數(shù)求解2.3.1水的熱力學性質(zhì)（1）飽和溫度和飽和壓力水由液態(tài)水轉(zhuǎn)變?yōu)闅鈶B(tài)水蒸氣的過程稱為汽化。汽化過程可以通過水表面的蒸發(fā)，也可以通過水的內(nèi)部產(chǎn)生汽泡的沸騰來形成。水

31、表面附近動能較大的水分子克服表面張力逸出到上面的空間稱為蒸發(fā)，蒸發(fā)可以在常溫下發(fā)生。同時也有空間的水蒸汽分子碰撞回到水面，稱之為凝結(jié)或液化。當水面上水蒸氣分子數(shù)目較少時，蒸發(fā)出來的分子較多，即蒸發(fā)速度大；當水面上水蒸氣分子數(shù)目逐漸增多，蒸汽壓力增大，返回水面的蒸汽分子數(shù)目增多，即凝結(jié)速度大。當逸出分子數(shù)目等于返回分子數(shù)目，即兩者達到動態(tài)平衡時，稱之為飽和狀態(tài)。水面上的水蒸氣稱之為飽和蒸汽，液態(tài)水稱之為飽和水。此時汽水溫度相同，稱為飽和溫度ts，水面上的蒸汽壓力稱為飽和壓力ps。例如一燒杯水在30時，燒杯水面上的水蒸氣壓力達到0.004257mpa時，燒杯內(nèi)的水不會再減少，它們分別為此狀態(tài)下的t

32、s和ps。在30溫度下，水蒸氣壓力與密度（分子數(shù)目）達到了最大值，如果再有蒸汽分子逸出，必有相同數(shù)目的蒸汽分子返回，故稱為“飽和”。工程上所用的大量蒸汽都是在鍋爐中加熱產(chǎn)生的。對水加熱，不但使表面發(fā)生汽化，而且水的內(nèi)部也產(chǎn)生汽泡，形成強烈沸騰。在此過程中，水內(nèi)部汽泡內(nèi)的蒸汽壓力等于或稍大于汽泡外壁所受的壓力，汽泡升至水面而破裂，隨之水蒸汽進入空間。因飽和蒸汽壓力取決于溫度，故汽泡的形成也只能發(fā)生在與給定壓力相對應(yīng)的飽和溫度，即該壓力下水的沸點。還以燒杯中的水為例。燒杯中的水內(nèi)部要產(chǎn)生汽泡，必須讓燒杯加熱，燒杯敞口，水受到1bar大氣壓力，1bar壓力相對應(yīng)的飽和溫度為=99.634，因此水必須

33、加熱到99.634以上，才可能在水的內(nèi)部產(chǎn)生汽泡而沸騰。水的飽和溫度與飽和壓力之間為非線性的單調(diào)變化關(guān)系，即飽和壓力增大，對應(yīng)的飽和溫度也升高。已編制成圖和表，可由查，也可由查。也有許多繁簡程度不同的函數(shù)關(guān)系式，即=f()或=f()?？梢愿鶕?jù)不同精度的要求，加以選擇使用。（2）水的臨界點當水的壓力達到22.064mpa，把水加熱到374，水立即全部汽化，沒有飽和水和飽和蒸汽的區(qū)別。該狀態(tài)稱為臨界點，臨界點的參數(shù)為：臨界壓力pc=22.064mpa，臨界溫度tc=374。（3）水的定壓加熱（蒸汽定壓發(fā)生）過程水進入鍋爐以后被加熱，在鍋爐出口以蒸汽形式流出，忽略過程中的壓力損失，可以認為是定壓加熱

34、過程，經(jīng)歷了預(yù)熱、汽化和過熱三個階段，前后共有五種狀態(tài)，見圖5。水從未飽和狀態(tài)（又稱過冷水）加熱到該壓力所對應(yīng)的飽和溫度，成為飽和液體（常稱飽和水）。繼續(xù)吸熱，飽和水開始沸騰，出現(xiàn)越來越多的汽泡，此時定壓又定溫，形成飽和水和飽和蒸汽的兩相混合物，稱為濕飽和蒸汽，簡稱濕蒸汽。繼續(xù)吸熱，直到液體水全部汽化成水蒸氣，此時已不含液態(tài)水，稱為干飽和蒸汽。對干飽和蒸汽繼續(xù)加熱，蒸汽溫度將從飽和溫度開始不斷升高，超過該壓力下的飽和溫度，稱為過熱蒸汽。圖5 水蒸汽定壓發(fā)生過程示意圖將水的定壓加熱和蒸汽定壓發(fā)生過程表示在壓力比容圖（-圖）和溫熵圖（-圖）上，如圖6所示。蒸汽的-圖和-圖可總結(jié)為一點、兩線、三區(qū)、

35、五態(tài)。一點是臨界點；兩線是飽和水線（下界限線）和飽和蒸汽線（上界限線）；三區(qū)是未飽和水區(qū)（過冷水區(qū)）下界限線左方，濕飽和汽區(qū)在上、下界限線之間，過熱蒸汽區(qū)在上界限線右方；五態(tài)（五種狀態(tài)）指的是：未飽和水狀態(tài)、飽和水狀態(tài)、濕飽和蒸汽狀態(tài)、干飽和蒸汽狀態(tài)和過熱蒸汽狀態(tài)。圖6 蒸汽的-圖和-圖（4）汽化潛熱在某一壓力下，1kg飽和水轉(zhuǎn)變?yōu)?kg干飽和蒸汽需要吸收的熱量稱為汽化潛熱，單位為kj/kg。同樣，1kg干飽和蒸汽凝結(jié)成1kg飽和水放出的熱量也與汽化潛熱相等。汽化潛熱的數(shù)值與壓力有關(guān)（見圖7），壓力越低，汽化潛熱越大。壓力達到臨界壓力=22.064mpa，汽化潛熱=0 kj/kg，即飽和水一加

36、熱就立刻變?yōu)楦娠柡驼羝?mpa)22.064rkj/kg2500圖7 汽化潛熱與壓力的變化曲線不同壓力下的汽化潛熱r可查表或函數(shù)關(guān)系式r=f(t)計算得到。（5）水和水蒸氣的主要熱力學性質(zhì) 本課題采用的是高壓下的熱水，但熱水的溫度為7595，沒有達到飽和溫度，所以應(yīng)該屬于未飽和水狀態(tài)（過冷水）。未飽和水的熱力學性質(zhì)有壓力p、溫度t、比容v、比熱焓h和比熵s等。其中壓力p和溫度t往往為已知變量，經(jīng)常使用的是比熱焓h，其單位為kj/kg。未飽和水的焓h由壓力p和溫度t確定，即。未飽和水的焓與壓力、溫度的關(guān)系比較復(fù)雜。現(xiàn)已制成圖和表，可由p，t數(shù)據(jù)直接查未飽和水和過熱蒸汽的性質(zhì)表。也有一些繁簡程度

37、不同的函數(shù)計算式，直接從這些計算式得到它們的焓值。2.3.2未飽和水的熱遷移性質(zhì)在熱采工藝的設(shè)計和運行管理中，人們最關(guān)心的是熱流體注入油層到底有多少熱量，在經(jīng)地面管線和井筒的過程中損失了多少熱量，以及如何減少這些熱量損失。為此，除了必須了解熱流體的熱力學性質(zhì)外，還須了解其熱遷移性質(zhì)，其中主要是粘度及導熱系數(shù)等。未飽和水的熱遷移性質(zhì)可以從傳熱學等專著查得。也可以從一些函數(shù)計算式中計算得到。由于熱流體循環(huán)過程中，熱水的壓力、溫度沿程變化，計算過程中需要多次反復(fù)進行。查表的方法不能適應(yīng)計算要求，因此比較精確的未飽和水的熱力學性質(zhì)、熱遷移性質(zhì)的函數(shù)計算式顯得很有用。原任職于美國德士古石油公司的美籍華人

38、錢思復(fù)的論文被廣泛引用來計算水和水蒸氣各種熱物性參數(shù)，可根據(jù)精確度的要求加以選擇使用。第三章 抽油工藝分析在稠油開采區(qū)，由于稠油粘度高，在有桿泵抽汲過程中出現(xiàn)較大的摩擦力，導致抽油機驢頭懸點載荷的最大值大大增加，而懸點最小載荷顯著地減小。抽油機驢頭懸點最大載荷不僅是設(shè)計或選擇抽油機和確定電動機動率的依據(jù)，也是設(shè)計和選擇抽油桿的重要依據(jù)。一般而言，懸點載荷的最大值表征了可供選用的抽油桿和抽油泵組合時可能的最大下泵深度。但是對于稠油并，較大的摩擦力使懸點最大載荷劇增，于是抽油機的有效提升能力降低，下泵深度就減小。因此，實測并研究稠油開采中抽油機懸點載荷的有關(guān)參數(shù)，計算出摩擦力的大小以及對懸點載荷的

39、影響，無論對采油機械的設(shè)計還是稠油開采都具有現(xiàn)實意義。由于采用新閉式熱流體循環(huán)降粘工藝，所以應(yīng)該對抽油桿柱載荷進行分析，一是通過載荷計算，看工藝參數(shù)是否滿足電機功率的要求；二是便于與以往降粘工藝的比較。3.1抽油桿柱的載荷計算有桿泵采油在運行時，抽油機爐頭懸點上作用的載荷有三類：（1）靜載荷，包括抽油桿自重、柱塞上部液柱形成的靜液柱載荷。（2）動載荷，包括抽油桿柱和油管內(nèi)的流體做不等速運動而產(chǎn)生的抽油桿和液柱的動載荷，即慣性載荷。（3）摩阻載荷，包括光桿和密封盒的摩擦力，抽油桿與液體之間的摩擦力，抽油桿和油管之間的摩擦力，液體在桿管環(huán)形空間的流動阻力，液體通過泵閥和柱塞內(nèi)孔的局部水力阻力，柱塞

40、和泵筒之間半干摩擦阻力。3.1.1懸點靜載荷懸點靜載荷計算用下式表示：上行程 下行程 （34）式中，為上行程懸點靜載荷，n；為下行程懸點靜載荷，n；為油管內(nèi)流體密度，kg/m3；為每米抽油桿在液體中重力，組合桿分段計算，kgf/m；為泵深，m；為動液面，m；為柱塞面積，m2；為油管壓力，mpa；為套管壓力，mpa。3.1.2懸點慣性載荷懸點慣性載荷計算用下式表示：上行程 下行程 （35）其中，式中，為上行程懸點最大慣性載荷，n；為下行程懸點最大慣性載荷，n；為沖程，m；為沖數(shù)，min-1；為每米抽油桿在空氣中重力，kgf/m；為考慮油管過流斷面變化引起液柱加速度變化的系數(shù)；為油管的流通斷面面積

41、，m2；為平均抽油桿截面積，m2。3.1.3懸點振動載荷實際的抽油桿柱和液柱，由于它們長度很長，具有相當?shù)膹椥院涂蓧嚎s性，而抽油桿柱做周期性地上、下運動和液柱載荷周期性地作用于下端，使抽油桿產(chǎn)生彈性震動，同時液柱下端周期性地被柱塞推動而使液柱也產(chǎn)生振動，如果油管下部未錨定，在液柱載荷周期性地作用下，管柱也要產(chǎn)生振動。這三組彈性體的振動互相影響，加上阻尼作用，使得整個系統(tǒng)的振動作用相當復(fù)雜，因此，要準確地計算彈性震動載荷是很困難的，這里僅介紹一種忽略強迫振動項的簡單計算方法。 （36）其中，式中，為在抽油桿頂端產(chǎn)生的振動載荷，n；為常數(shù)，其值取決于，當此數(shù)值為01時=0，為13時=1，為35時=

42、2，為57時=3；為鋼材彈性模量，kn/m2；為抽油桿變形，m；為靜變形結(jié)束時的懸點速度，為上沖程，為下沖程，m/s；為前臂長，m；為抽油桿內(nèi)聲波傳播速度，m/s；為曲柄轉(zhuǎn)角，；為靜變形結(jié)束時的曲柄轉(zhuǎn)角，。3.1.4摩擦載荷摩擦載荷一般由5部分組成：抽油桿和油管之間的摩擦力，根據(jù)現(xiàn)場經(jīng)驗，在直井內(nèi)通常不超過抽油桿重力的1.5%；柱塞和泵筒之間半干摩擦阻力，根據(jù)現(xiàn)場經(jīng)驗，當泵徑不超過70mm時，半干摩擦力不超過1717n；抽油桿與液體之間的摩擦力，可用式（37）計算；液柱與油管之間的摩擦力，根據(jù)高粘度油井的現(xiàn)場資料統(tǒng)計，約等于0.77；液體通過游動閥的摩擦阻力，一般可忽略不計。 （37）式中，m

43、為油管內(nèi)徑與抽油桿直徑之比；為液體動力粘度，pa.s。這樣，摩擦載荷可表示為上行程 下行程 （38）式中，為上沖程摩擦載荷，n；為下沖程摩擦載荷，n。3.2懸點最大和最小載荷 （39）式中，為懸點最大載荷，n；為懸點最小載荷，n。3.3抽油機扭矩各種抽油機在運行過程中，曲柄軸上承受的扭矩，包括由懸點載荷所形成的負荷扭矩和由平衡重所形成的平衡扭矩，這兩項疊加形成凈扭矩，即曲柄軸上實際承受的扭矩。凈扭矩可以用下式計算。 （40）式中，為凈扭矩，kn.m；為前臂長，m；為后臂長，m；為游梁平衡重心至游梁支點o的距離，m；p為懸點載荷，kn；為游梁平衡重，kn；為從水平位置算起的游梁擺角，；為懸點加速

44、度，m/s2；g為重力加速度，m/s2；r為曲柄旋轉(zhuǎn)半徑，m；為曲柄與連桿的夾角，；為連桿與后臂的夾角，；為曲柄平衡重，kn。以上計算方法是確定扭矩最準確的方法，但計算比較繁瑣，在選擇抽油設(shè)備和預(yù)測新抽汲參數(shù)時需要知道最大扭矩，一般采用直接計算方法，可用下式進行。 （41）式中，為最大扭矩，kn.m。但是，式（41）計算結(jié)果偏低11%50%不等。因此根據(jù)國內(nèi)礦場經(jīng)驗，統(tǒng)計最大扭矩經(jīng)驗公式，結(jié)果比較切合實際。 （42）第四章 計算程序開發(fā)與設(shè)計4.1 軟件的研制開發(fā)目的井筒密閉熱流體循環(huán)降粘工藝設(shè)計軟件目的在于以密閉熱流體循環(huán)降粘工藝原理為基礎(chǔ)，優(yōu)選出經(jīng)濟的井筒熱流體循環(huán)降粘工藝參數(shù)，為現(xiàn)場的施

45、工作業(yè)提供技術(shù)支持。同時將這一過程計算機化，利用計算機的速度快、精度高、能夠?qū)崿F(xiàn)一些較為復(fù)雜的數(shù)學運算與邏輯判斷等特性，實現(xiàn)對密閉熱流體循環(huán)降粘工藝的分析計算。4.2 軟件需求分析與功能說明軟件需求：1、 接受現(xiàn)場的稠油井參數(shù)輸入到軟件的輸入界面，根據(jù)密閉熱流體循環(huán)降粘工藝理論進行分析計算，優(yōu)化出合理的施工參數(shù)并輸出；2、 軟件要便于操作和使用。鑒于此，我們對軟件的功能進行了如下設(shè)計：1、 能夠接受稠油井的數(shù)據(jù)輸入，根據(jù)密閉熱流體循環(huán)降粘理論進行分析計算；2、 計算的結(jié)果能方便輸出；3、 每項功能都提供菜單，簡化軟件的使用和操作。4.3 開發(fā)環(huán)境的選擇軟件采用windows窗體界面與目前流型軟

46、件接軌，既實用又具備一定的先進性為原則，本軟件開發(fā)環(huán)境如下：1、 操作系統(tǒng)，microsoft windows2000及以上操作系統(tǒng)。微軟的操作系統(tǒng)具有極大的用戶群，而且操作界面具有同一性，極易操作。windows系統(tǒng)具有較好的運行性能，而且目前大多數(shù)計算機上（包括許多品牌機在內(nèi)）安裝的依然是該版本的操作系統(tǒng)。2、 編程環(huán)境，visual basic 6.0。該產(chǎn)品是 microsoft公司推出的visual studio 6.0開發(fā)組件中的一種可視化的編程語言。3、 編程語言，visual basic 6.0。具有先進的可視化編程工具和面向?qū)ο蟮木幊坦ぞ?，便于本軟件的各項功能的實現(xiàn)。4.4

47、軟件設(shè)計通過對軟件需求、功能要求和界面設(shè)計等多方面內(nèi)容的全面考慮，利用面向?qū)ο蟮木幊谭椒▽υ撥浖M行了設(shè)計，其軟件功能圖如圖8所示：井筒降粘優(yōu)化數(shù)據(jù)輸出模塊主應(yīng)用程序模塊界面顯示模塊數(shù)據(jù)讀入模塊幫助模塊管參數(shù)繪制圖8 軟件功能結(jié)構(gòu)圖4.5 軟件的使用4.5.1軟件的安裝 將安裝盤放入驅(qū)動器，運行其中的setup.exe文件，軟件即會自動安裝在您的計算機中。并在“開始/程序”菜單建立快捷方式。4.5.2進入系統(tǒng) 在windows桌面上，點擊“開始”欄，選擇“程序”，找到【井筒密閉熱流體循環(huán)降粘工藝設(shè)計軟件】，打開本程序。首先顯示一個進入窗體，然后進入軟件的工作平臺。4.5.3降粘計算過程 (1)

48、 基礎(chǔ)數(shù)據(jù)輸入l 選擇“數(shù)據(jù)文件”主菜單“新建文件|打開”菜單，彈出一個窗體，該窗體用來確定設(shè)計有關(guān)基礎(chǔ)數(shù)據(jù)，包括生產(chǎn)參數(shù)，管柱參數(shù)，機采參數(shù)，熱水循環(huán)參數(shù)，地層參數(shù)等。點擊“確定”按鈕，即可完成參數(shù)輸入。l 也可通過該窗體的“打開”按鈕，將以前輸入的數(shù)據(jù)文件打開進行參數(shù)輸入。 (2) 井筒降粘設(shè)計l 定管柱參數(shù)：選擇“降粘計算”主菜單“管柱參數(shù)計算”菜單，下拉子菜單“計算開始”，進行優(yōu)化設(shè)計，經(jīng)過一定時間后，計算結(jié)束。可以在“管柱參數(shù)計算”菜單下看到其它兩個子菜單“曲線繪制”、“結(jié)果保存”。分別顯示在管柱參數(shù)確定條件下計算得到的計算結(jié)果，繪制的井筒壓力、溫度分布圖。 (3) 管柱圖繪制l 選

49、擇“生成管柱圖”主菜單，彈出一個窗體，該窗體用來根據(jù)管柱參數(shù)繪制管柱圖。首先在窗體右側(cè)輸入管柱設(shè)計的相關(guān)參數(shù)，點“畫圖”按鈕即可完成管柱圖繪制，點擊“保存”按鈕即可保存圖片文件為bmp和jpg格式。4.5.4操作l 數(shù)據(jù)輸入：當需要進行分析時，首先需要打開基礎(chǔ)數(shù)據(jù)輸入窗體，輸入與分析相關(guān)的參數(shù)。輸入數(shù)據(jù)有三種不同的方式：根據(jù)需要輸入數(shù)據(jù)；根據(jù)需要修改默認值；打開以往數(shù)據(jù)。l 數(shù)據(jù)輸出：點擊各主菜單下對應(yīng)子菜單瀏覽；打開相應(yīng)的曲線圖進行查閱。4.5.5退出 進入“數(shù)據(jù)文件”菜單后，選擇“退出系統(tǒng)”或點擊主窗體右上角的“”，退出系統(tǒng)。4.6 軟件特點l 語言：vb語言，該語言是一種程序設(shè)計的高級語言，vb語言具有很多優(yōu)點：可以編譯、數(shù)據(jù)類型豐富、支持結(jié)構(gòu)化程序設(shè)計和面向?qū)ο蟪绦蛟O(shè)計。此外，vb語言易于閱讀，編譯迅速，并且可以使用多個單元文件來進行模塊化程序設(shè)計。l 標準：windows操作系統(tǒng)標準。軟件的操作與windows系統(tǒng)一致，可參閱windows系統(tǒng)操作手冊。l 操作方便：采用了下拉式菜單、右鍵彈出菜單，同時設(shè)計了快捷鍵、快捷按鈕等，并有相關(guān)的提示。l 界面美觀。4.7 軟件的界面（1）啟動窗體：（2）軟件主平臺（3）數(shù)據(jù)輸入界面（4）計算結(jié)果繪圖界面（5）計算結(jié)果保存（6