石端勝，袁偉杰，王宏申，劉躍忠，王錦林，張志軍，叢越男

（1.中海油能源發(fā)展股份有限公司工程技術分公司，天津300452；2.中海石油（中國）有限公司秦皇島32-6 作業(yè)公司，天津300450）

渤海B油田主要開發(fā)沙河街組， 膠質瀝青質含量較高，為5.33%～32.4%，原始地層壓力47 MPa，平均地層壓力30 MPa，油藏溫度120 ℃。隨著油田不斷的生產(chǎn)開發(fā)，地層壓力逐漸下降，在近井地帶及井筒附近膠質瀝青質容易析出，對油井產(chǎn)能造成嚴重影響［1-4］。本文利用氣-液-固三相平衡熱力學方程判別油田膠質瀝青沉積析出，由此建立單井徑向流模型，進行數(shù)值模擬研究，對膠質瀝青沉積規(guī)律進行模擬分析和預測，旨在為油田開展相關解堵工藝措施提供重要依據(jù)。

1 膠質瀝青質沉積模擬計算判別

1.1 模擬計算目的

通過建立井下油管內膠質瀝青沉積過程氣-液-固三相平衡熱力學模型，可以描述溫度、壓力和組成等熱力學條件對氣-液-固三相相態(tài)的影響，以及在氣-液-固三相相態(tài)轉換過程中和在氣-液-固三相平衡時各組分的變化規(guī)律，以確定膠質瀝青沉積物的形成原因及關鍵因素臨界點［5-7］。

1.2 模擬計算方法

（1）在求解氣-液-固三相平衡問題時，采用PR狀態(tài)方程計算氣-液平衡。

（2）對液-固平衡，將固相視為非理想的固態(tài)溶液，應用正規(guī)溶液理論來處理。

（3）用有限的高壓模擬實驗數(shù)據(jù)回歸熱力學模型中的參數(shù)。

（4）原油組分特征化方法采用前面原油族組成分析結果。

1.3 熱力學計算模型

根據(jù)熱力學相平衡原理，氣、液、固三相的逸度分別表示為：

氣相和液相逸度計算（PR狀態(tài)方程）：

式中am=a（T）ac

液相逸度計算（PR狀態(tài)方程）：將上式按壓縮因子Z的形式表示，則有

狀態(tài)方程參數(shù)am和bm可由純組分的相應參數(shù)ai和bi按混合規(guī)則求定：

混合物中i組分的逸度系數(shù)的表達式為：

膠質瀝青液相逸度：

固相膠質瀝青逸度計算：

膠質瀝青是否發(fā)生沉積的判斷條件是：ln fa≥lnfa*即膠質瀝青液相逸度大于瀝青固相逸度， 則膠質瀝青發(fā)生沉積；lnfa＜lnfa*即膠質瀝青液相逸度小于瀝青固相逸度， 則膠質瀝青不發(fā)生沉積；lnfa=lnfa*對應于膠質瀝青的沉積初始壓力。

1.4 熱力學計算結果

熱力學模型計算選取B油田原油組分組成實驗分析數(shù)據(jù)，進行原油組分擬特征化，利用熱力學模型進行求解，見表1～3。

表1 B油田原油組分組成實驗分析數(shù)據(jù)

在目前渤海B油田平均地層壓力30 MPa， 油藏溫度120 ℃左右條件下，計算得出ln fa≤lnfa*，即膠質瀝青液相逸度大于膠質瀝青固相逸度，則膠質瀝青發(fā)生沉積。

表2 原油組分擬特征化結果

表3 膠質瀝青沉積條件計算結果

2 單井徑向流模型的建立

由熱力學計算模型得知B油田膠質瀝青液相逸度大于膠質瀝青固相逸度，易發(fā)生膠質瀝青沉積。開展膠質瀝青沉積數(shù)模研究分析［8-10］，結合近井地帶壓力梯度分布， 得出油井膠質瀝青沉積規(guī)律，以及膠質瀝青在產(chǎn)層流體流動中對產(chǎn)能的影響。

2.1 數(shù)值模擬理論基礎

運用三相平衡原理，當油藏條件改變打破平衡時，便發(fā)生反應，產(chǎn)生膠質瀝青沉積。

式中：Ci為反應i的濃度因子；?為孔隙度，%；ρj為j相密度，g/cm3；Sj為j相飽和度，%。 膠質瀝青質析出反應方程式平衡常數(shù)見表4。

式中：Xji為組分i在j相中的摩爾分數(shù)；Xji平衡為相平衡常數(shù)K-value計算的平衡摩爾分數(shù)。

表4 膠質瀝青質析出反應方程式平衡常數(shù)

膠質瀝青沉積造成滲透率降低的理論：生產(chǎn)過程中，隨著流體的注入/采出，儲層內外壓差發(fā)生變化，孔隙受到壓縮/膨脹影響而體積改變，特別是有瀝青沉積的油藏，孔隙度變化更大。 而滲透率與孔隙度有一定的函數(shù)關系，孔隙度的變化也會使?jié)B透率隨之變化。

式中：k為變化后的滲透率，10-3μm2；k0為初始滲透率，10-3μm2；φ為變化后的孔隙度，%；φ0為初始孔隙度，%；kmul為倍乘系數(shù)。

2.2 數(shù)值模型建立

采用單井徑向流模型模擬近井地帶隨著生產(chǎn)開發(fā)進行膠質瀝青的沉積規(guī)律。

單井徑向流模型：

R×J×K：50×1×10

半徑：300 m

近井帶網(wǎng)格半徑：0.2 m

滲透率：13×10-3μm2

孔隙度：25%

深度：3 500 m

溫度：120 ℃

原始平均地層壓力：47 MPa

3 模型結果

3.1 壓力對膠質瀝青沉積量的影響

隨著油田開發(fā)生產(chǎn)時間的延長， 近井帶壓力逐步下降，膠質瀝青在近井帶析出并開始沉積，之后沉積量逐漸增大，同時滲透率大幅減小，見圖1。

圖1 生產(chǎn)中近井帶壓力、滲透率、膠質瀝青沉積量變化關系

當壓力下降到36～40 MPa時瀝青開始析出，同時近井帶滲透率大幅下降；壓力繼續(xù)下降膠質瀝青析出速度增加，滲透率下降幅度也增大；當壓力下降到28～33 MPa時， 膠質瀝青大量沉積造成地層堵塞，壓力下降變慢，滲透率下降幅度變小，見圖2。

圖2 不同地層壓力條件下滲透率、膠質瀝青沉積量變化關系

3.2 膠質瀝青沉積對產(chǎn)量的影響

分別對15，50，75 m3/d不同產(chǎn)量條件下進行膠質瀝青沉積預測，可知產(chǎn)量越低，膠質瀝青沉積量越大，說明產(chǎn)量降低導致地下滲流降低，增大了沉積時間，使得膠質瀝青沉積量增大，見圖3。

圖3 不同產(chǎn)量條件下膠質瀝青沉積情況

此外，當有膠質瀝青析出時，穩(wěn)產(chǎn)時間變短，產(chǎn)油量大幅下降，下降幅度高達84.64%，見圖4。

圖4 膠質瀝青沉積對產(chǎn)量影響規(guī)律

3.3 膠質瀝青質沉積對儲層傷害模擬

對比油井生產(chǎn)200 d和生產(chǎn)400 d近井地帶地層壓力、沉積量分布和滲透率變化情況，不難發(fā)現(xiàn)隨著油井生產(chǎn)時間的增加，近井帶壓力下降，膠質瀝青在近井帶開始沉積，同時近井帶滲透率也開始降低，對儲層滲透率傷害增大，見圖5。

圖5 近井地帶滲透率分布

4 結論

（1）通過熱力學模型計算，B油田膠質瀝青液相逸度大于膠質瀝青固相逸度，容易發(fā)生膠質瀝青質沉積。

（2）構建單井徑向流模型模擬膠質瀝青質沉積規(guī)律，研究結果表明： 膠質瀝青質的析出和沉積受地層壓力影響明顯。 膠質瀝青質開始析出的地層壓力范圍為36～40 MPa，膠質瀝青沉積速度加快的地層壓力范圍為28～33 MPa。

（3）隨著油井生產(chǎn)時間增加，近井帶壓力下降，膠質瀝青在近井帶開始沉積，同時近井帶滲透率也開始降低，對儲層滲透率傷害增大，對油井產(chǎn)能影響明顯，油井產(chǎn)能下降幅度可達80%以上。