趙光華

(長城鉆探工程技術(shù)研究院，遼寧 盤錦 124001)

二氧化碳(CO2)吞吐可以用于多種條件的油藏[1]，如粘度不太高的稠油油藏，水驅(qū)油藏(稀油)高含水時期，油層連通性差、非均質(zhì)嚴重、無法形成注采系統(tǒng)的小斷塊，塊狀、具有頂部剩余油的傾斜油層等。中原油田胡12塊屬于非均質(zhì)嚴重的中高滲透油藏[2]，注水效果差，現(xiàn)場試注N2驅(qū)油效果不好，目前已進入高含水階段，急待推廣應用CO2吞吐采油技術(shù)。大量研究表明，分子擴散是CO2提高原油采收率的一個重要機理[3-4]。因此，在高壓和不同溫度下，CO2在原油中的擴散系數(shù)是決定CO2在原油傳質(zhì)速率的一個很重要的參數(shù)。由于地層基巖很復雜 ，注入的CO2很難與油藏中原油完全混合好[5]。多數(shù)情況下，CO2通過分子擴散溶于原油，當水相將油相與CO2氣相隔開時，水相阻礙了CO2分子的擴散，并且完全抑制了輕質(zhì)烴類化合物釋放到CO2中。然而，當CO2氣體連續(xù)注入時，CO2的分子擴散能夠通過水層將殘余油驅(qū)替，從而提高原油采收率。

氣體在原油中擴散系數(shù)的實驗測定方法可大體歸納為傳統(tǒng)方法和非傳統(tǒng)方法。傳統(tǒng)方法主要是測定氣體在不同時間的組成(氣相色譜[6]或核磁共振[7])，得到擴散系數(shù)。非傳統(tǒng)方法建立在測定由分子擴散造成的氣-液體系某一性質(zhì)的改變，包括溶質(zhì)的揮發(fā)速率[8]、氣體體積[9]、折射率[10]、毛細管內(nèi)原油-氣體的界面位置[11]、氣體包圍下懸油滴的形狀[12]及體積[13]等。本文采用動態(tài)界面張力法測定CO2在原油中的擴散系數(shù)，相對于其他實驗方法，該方法具有操作時間短、恒壓恒溫下通過測定的動態(tài)界面張力值就能同時計算出擴散系數(shù)和界面?zhèn)髻|(zhì)系數(shù)等優(yōu)點。本文的研究對象為中原油田胡12塊輕質(zhì)原油。

1 模型基礎

圖1(a)描述了在某一溫度、壓力下一個形狀完整的懸油滴。圖1(b)為周圍分布CO2的軸對稱懸油滴示意。注射探針的內(nèi)徑用rn表示，注射探針壁厚用εn表示，注射探針高度用hn表示，懸油滴包括探針內(nèi)部的油相用ω表示，注射探針頂部的切平面和探針及原油形成的界面用φn表示。

圖1 (a)一個完整懸油滴的數(shù)碼圖像；(b)圓柱坐標體系(r，z)下軸對稱懸油滴在CO2相中示意

本次研究中，懸油滴形狀和體積的改變對擴散過程帶來的影響可以忽略。假設擴散過程中懸油滴是靜止的，根據(jù)Ghez的氣體擴散理論[14]，用于描述CO2在懸油滴內(nèi)部(包括注射探針內(nèi)部油相)的濃度轉(zhuǎn)移擴散公式可以表示成:

(1)

式中，D為CO2在原油中的擴散系數(shù)，m2/s，假設它與CO2的濃度無關(guān)。

在傳質(zhì)過程的開始階段，懸油滴內(nèi)部CO2的濃度為0，因此，公式(1)的初始條件為[14]：

c(r,z,t)|t=0=0,(r,z)∈ω

(2)

注射針壁及切平面對于CO2來說是不可滲透的，相應的非滲透邊界條件為：

(3)

式中，nr和nz為方向余弦，即注射探針和原油以及切平面形成界面的單位矢量在r軸和z軸形成的分量。

本次研究考慮CO2通過原油-CO2界面的界面阻力，因此，非平衡邊界條件應用于界面得到[14]：

(4)

式中，k為界面?zhèn)髻|(zhì)系數(shù)，m/s，ceq為CO2在指定壓力和溫度下在原油中的平衡濃度，mol/m3。

公式(1)～(4)可以通過引進4個無量綱變量和1個無量綱參數(shù)將其無量綱化：

(5)

式中，C為CO2的無量綱濃度；R和Z為無量綱的徑向和軸向坐標；rn為注射探針的內(nèi)徑；τ為無量綱時間；kD為傳質(zhì)Biot數(shù)；實際上，傳質(zhì)Biot數(shù)表示本體阻力rn/D與界面阻力1/k的比值[15]。

本次研究應用半離散Galerkin有限元素法來解決上述無量綱公式,具體的數(shù)學處理過程可以參照文獻[16]。

懸滴在t=0的第一個數(shù)碼圖像用來定義每個壓力下的傳質(zhì)模型的物理范圍ω，如圖1(b)所示。然后，擴散公式(1)連同公式(2)～(5)是通過采用半離散Galerkin有限元素法得到CO2在計算域Ω下的無量綱濃度轉(zhuǎn)移C(R,Z,τ)和CO2在界面處的平衡濃度Cint(τ)。CO2的平均無量綱界面濃度能夠轉(zhuǎn)換成真實時間t下濃度Cint(t)。本次研究中，假設動態(tài)界面濃度和動態(tài)界面張力的關(guān)系遵循平衡界面濃度和平衡界面張力的關(guān)系[17]，CO2在原油中的平衡濃度(溶解度)可通過Simon-Graue關(guān)系式[18]計算出。那么，動態(tài)界面張力就能夠通過計算的動態(tài)界面濃度Cint(τ)計算出來。

用目標函數(shù)表示數(shù)學計算的動態(tài)界面張力和實驗測得數(shù)據(jù)的偏差。令ti和γm(ti)，i=1,2,3,…,N是一組時間和實驗測得的動態(tài)界面張力，tj和γc(tj)，j=1,2,3,…,K是另一組時間和理論計算的動態(tài)界面張力。這里，目標函數(shù)被定義為：

(6)

從公式(6)可以看出目標函數(shù)是計算和測得的動態(tài)界面張力的均方根相對誤差。計算的動態(tài)界面張力隨著計算的CO2界面濃度Cint(t)的變化而改變，而CO2在界面的濃度是由傳質(zhì)Biot數(shù)kD和擴散系數(shù)D共同決定的。因此，目標函數(shù)值是由傳質(zhì)Biot數(shù)kD和擴散系數(shù)D共同決定的。

E=E(kD,D)

(7)

一旦全局目標函數(shù)能夠最小化，kD和D就作為真實的傳質(zhì)Biot數(shù)和真實的擴散系數(shù)。

本次研究采用一個參數(shù)同時查找方法，計算傳質(zhì)Biot數(shù)kD和D擴散系數(shù)的數(shù)學過程如下：

(1)周圍分布CO2的懸油滴的連續(xù)數(shù)碼圖像能夠得到，動態(tài)界面張力γm(ti)能夠測得。

(2)懸油滴在t=0的第一個數(shù)碼圖像定義為傳質(zhì)模型的物理范圍ω，如圖1(b)所示。

(3)在一個給定的溫度壓力下，傳質(zhì)Biot數(shù)的11個測試點選為kD,i，i=0,1,2,…10。

(4)對于一個給定的kD,i，擴散公式(1)連同公式(2)-(5)是通過半離散Galerkin有限元素法得到CO2在計算域Ω中的無量綱濃度轉(zhuǎn)移C(R,Z,τ)。

(5)對于一個給定的kD,i，擴散系數(shù)D有上下限[D0,D10]。

(6)D的初始不確定時間間隔被11個節(jié)點Dj,j=0,1,2,…10劃分為10個等同的子區(qū)間。

(7)對于一個給定的Dj，無量綱時間τ能夠轉(zhuǎn)換為真實時間t。應用測得的平衡界面張力和計算的CO2平衡濃度的關(guān)系曲線，任何時間t下的動態(tài)界面張力就能夠計算出來。利用測得和計算的動態(tài)界面張力，通過公式(6)得到相應的目標函數(shù)E(kD,i,Dj),j=0,1,2,…10。

(8)通過對比不同的E(kD,i,Dj),j=0,1,2,…10，局部最小目標函數(shù)就得到了Emin(kD,i,Dj)，1≤j≤9。然后，D的不確定間隔變?yōu)閇Dj-1，Dj+1]。

(9)如果新的不確定間隔[Dj-1，Dj+1]足夠小，能夠滿足擴散系數(shù)的代數(shù)精度要求，那么，Dj就作為給定kD,i的局部最佳值。否則，新的不確定區(qū)間[Dj-1，Dj+1]將作為D的初始區(qū)間重復步驟(6)～(9)，直到滿足終止條件。

(10)對于所有的kD,i，i=0,1,2,…10，步驟(4)～(9)將重復進行。通過對比不同的E(kD,i,Dj),j=0,1,2, … 10，可得到局部最小函數(shù)Emin(kD,m,Dn)，1≤m≤9，kD,m和Dn就作為傳質(zhì)Biot數(shù)和擴散系數(shù)的最佳值。

(11)如果m=0或m=10，則kD,i，i=0,1,2,…10將重新被選取。對于一個新值kD,i，步驟(4)之后的數(shù)學最小化過程將重復進行直到滿足1≤m≤9。

對于原油-CO2體系kD和D的查找范圍為0.3～8.8和(0.01～10)×10-9m2/s。D的終止條件是︱Dn+1-Dn-1︱≤0.005×10-9m2/s。

2 實驗部分

實驗采用加拿大ROBINSON D B公司生產(chǎn)的JEFFRI高壓界面張力儀，該實驗裝置如文獻[19]所述，實驗用CO2由北京氦普北分氣體有限公司生產(chǎn)，純度為99.95%。油樣采用中原油田胡12塊輕質(zhì)原油，原油在300.15 K、常壓下的密度和黏度分別為0.837 g/cm3和13.0 mPa·s。該原油的組成分析結(jié)果見表1。

表1 中原油田胡12塊輕質(zhì)原油組成分析

C1～C6:7.90,C7+:92.10

本實驗采用懸滴法進行界面張力測定。應用ANDREAS等提出的選面法[20]計算界面張力，實驗方法在文獻[19]中已作詳細闡述。記錄300.15 K，0.5～4.1 MPa下CO2在原油中不同時間下的懸滴圖像并計算動態(tài)界面張力。

3 結(jié)果與討論

3.1 動態(tài)界面張力與平衡界面張力

在300.15 K下，不同壓力下動態(tài)界面張力隨時間的變化如圖2所示，動態(tài)界面張力隨著時間的增大逐漸趨于一個穩(wěn)定值，這個穩(wěn)定值就定義為平衡界面張力。從圖2可以看出動態(tài)界面張力隨著壓力的升高而降低，這是由于在高壓下CO2在原油中的平衡濃度更高，能夠有效降低原油的界面張力。

圖2 300.15 K，不同壓力下測得的動態(tài)界面張力(點)與計算的動態(tài)界面張力(線)隨時間的變化

動態(tài)界面張力在高壓下能夠更快地達到平衡界面張力，5個不同壓力下的平衡界面張力如圖3所示。從圖3可以看出平衡界面張力隨著壓力的升高而降低。通過Simon-Graue關(guān)系式[11]得出CO2在原油中的平衡濃度與壓力的關(guān)系，見圖3。

圖3 300.15 K下，CO2在原油中的平衡界面張力、平衡濃度與壓力的關(guān)系曲線

從圖3可以看出CO2的平衡濃度隨著壓力的升高而增大。在300.15 K下，CO2在原油中的平衡界面張力與平衡濃度的關(guān)系見圖4。圖4就作為300.15 K下測得的平衡界面張力與CO2在原油中平衡濃度的關(guān)系曲線。

圖4 300.15 K下，CO2在原油中的平衡界面張力與平衡濃度關(guān)系曲線

3.2 擴散系數(shù)和傳質(zhì)Biot數(shù)

通過測得的動態(tài)界面張力數(shù)值，目標函數(shù)E(D,kD)通過計算能夠最小化，這樣能夠確定CO2在原油中每個壓力下的擴散系數(shù)D和傳質(zhì)Biot數(shù)kD。在300.15 K，p=0.5 MPa下，原油-CO2體系在不同傳質(zhì)Biot數(shù)下的擴散系數(shù)D如圖5所示，從圖5可以看出，對于一個給定的傳質(zhì)Biot數(shù)，存在一個最佳擴散系數(shù)D使全局目標函數(shù)最小化。對于不同的傳質(zhì)Biot數(shù)，存在一對擴散系數(shù)D和傳質(zhì)Biot數(shù)kD最優(yōu)值，能夠使目標函數(shù)最小化。在300.15 K、p=0.5 MPa下，原油-CO2體系的最優(yōu)擴散系數(shù)D=0.48×10-9m2/s，最優(yōu)傳質(zhì)Biot數(shù)kD=4.0。

圖5 300.15 K,0.5 MPa下不同傳質(zhì)Biot數(shù)下目標函數(shù)與擴散系數(shù)的關(guān)系曲線

圖6 300.15 K,0.5 MPa下不同壓力下的全局最優(yōu)傳質(zhì)Biot數(shù)下的目標函數(shù)與擴散系數(shù)的關(guān)系曲線

圖6表示在p=0.5,1.5,2.1,3.1,4.1 MPa下，CO2的擴散系數(shù)D在全局最優(yōu)傳質(zhì)Biot數(shù)下不同的目標函數(shù)E(D,kD)。從圖6可以看出，每一條目標函數(shù)曲線都對應一個最佳擴散系數(shù)D，能夠使目標函數(shù)最小化。此時，擴散系數(shù)D和傳質(zhì)Biot數(shù)kD就作為真實值。本次研究中，300.15 K，0.5～4.1 MPa下，擴散系數(shù)D的計算結(jié)果為(0.45～4.61)×10-9m2/s，傳質(zhì)Biot數(shù)為2.5～7.1，界面?zhèn)髻|(zhì)系數(shù)為(0.85～8.33)×10-5m/s。

300.15 K，p=0.5,1.5,2.1,3.1,4.1 MPa下計算出的CO2在原油中的平衡濃度、擴散系數(shù)、傳質(zhì)Biot數(shù)、界面?zhèn)髻|(zhì)系數(shù)以及全局最小目標函數(shù)Emin(D,kD)值見表2。

從表2可以看出，理論計算與實驗測得的界面張力均方根相對誤差低于1.81%。

表2 300.15 K，不同壓力下的CO2在原油中的平衡濃度、擴散系數(shù)D、傳質(zhì)Biot數(shù)kD、界面?zhèn)髻|(zhì)系數(shù)k以及全局最小目標函數(shù)Emin(%)

從表2還可以看出,CO2在原油中的平衡濃度隨壓力的升高而增大，CO2在原油中的溶解能夠降低原油-CO2之間的界面張力，因為這樣CO2能夠更容易通過原油-CO2界面擴散進入原油相，造成高壓下有更大的擴散系數(shù)和界面?zhèn)髻|(zhì)系數(shù)，這樣，CO2-原油的動態(tài)界面張力在更短的時間內(nèi)達到平衡界面張力。

4 結(jié)論

本文采用動態(tài)界面張力法研究300.15 K，0.5～4.1 MPa下CO2在中原油田胡12塊輕質(zhì)原油中的擴散系數(shù)和界面?zhèn)髻|(zhì)系數(shù)，得到原油-CO2體系的擴散系數(shù)為(0.45～4.61)×10-9m2/s，傳質(zhì)Biot數(shù)為2.5～7.1，界面?zhèn)髻|(zhì)系數(shù)為(0.85～8.33)×10-5m/s，模型計算與實驗測得動態(tài)界面張力的均方根相對誤差低于1.81%，說明本文建立的傳質(zhì)模型具有較好的穩(wěn)定性，動態(tài)界面張力法在一定程度上能夠用于CO2在原油中擴散系數(shù)的測定。