蔡計光， 房好青， 蘇軍偉， 王樂

(1.中國石化西北油田分公司石油工程技術(shù)研究院， 烏魯木齊 830011； 2.中國石化縫洞型油藏提高采收率重點實驗室， 烏魯木齊 830011； 3.西安交通大學(xué)人居環(huán)境與建筑工程學(xué)院， 西安 710049; 4.西安石油大學(xué)機械工程學(xué)院， 西安 710065)

碳酸鹽巖縫洞型儲層受原始巖性、構(gòu)造和巖溶的綜合影響，具有空間孔隙結(jié)構(gòu)復(fù)雜、非均質(zhì)性強等特點[1]。酸化工藝是碳酸鹽巖縫洞型儲層地質(zhì)改造的主要手段[2-3]，對于酸化過程的研究，能夠揭示縫洞對碳酸鹽巖儲層酸化過程的影響機制，進而為優(yōu)化改良酸化工藝、提高酸化效果提供借鑒和指導(dǎo)。

目前，不同學(xué)者分別基于經(jīng)驗?zāi)Ｐ蚚4]、分形模型[5]、網(wǎng)絡(luò)模型[6-7]和雙尺度連續(xù)模型[8-9]，對碳酸鹽巖儲層酸化進行研究。紀國法等[10]基于徑向蚓孔擴展模型和等效滲透率模型建立了碳酸鹽巖基質(zhì)酸化產(chǎn)能預(yù)測新模型，為產(chǎn)能預(yù)測提供了理論依據(jù)。Tansey等[11]提出了基于有限元的非均相反應(yīng)遷移和基質(zhì)溶解的孔隙網(wǎng)絡(luò)模型，極大地改善了基質(zhì)溶解的準確預(yù)測。相對于孔隙網(wǎng)絡(luò)模型和經(jīng)驗?zāi)Ｐ?，Panga等[12]建立的雙尺度連續(xù)模型由于考慮了達西尺度和孔隙尺度特征，能夠再現(xiàn)酸和固體的質(zhì)量變化，有著較低的計算成本等優(yōu)勢，已經(jīng)逐漸發(fā)展為碳酸鹽巖儲層酸化研究的主流數(shù)值模型[13-16]。Luo等[8]研究了裂縫型雙礦物碳酸鹽巖中活性酸運移機理，發(fā)現(xiàn)白云石體積分數(shù)、白云石表面反應(yīng)常數(shù)和擴散系數(shù)的變化會改變?nèi)芙庑螒B(tài)和溶解性能。廖毅等[17]建立了理想天然裂縫模型，考慮了裂縫長度數(shù)量和方向?qū)λ嵛g蚓孔形態(tài)特征與酸液注入的影響。柳明等[18]基于雙尺度數(shù)值模型研究了碳酸鹽巖油藏酸化酸蝕蚓孔的分形型，發(fā)現(xiàn)面溶蝕、蚓孔和均一溶蝕的分形維數(shù)分別為1.46、1.50和1.44。齊寧等[19]利用離散裂縫網(wǎng)絡(luò)模型，求解了徑向流裂縫性地層蚓孔擴展問題，發(fā)現(xiàn)裂縫周圍存在局部控制域，控制域尺寸與裂縫在酸液注入方向上的投影有關(guān)，且蚓孔的生長軌跡由控制域覆蓋蚓孔尖端的裂縫決定，與其他裂縫無關(guān)。理論上來說，碳酸鹽巖縫洞型儲層中縫洞分布具有一定的隨機性，以往的研究一方面假設(shè)了裂縫的理想化分布，一方面缺少對于縫洞這一組合形態(tài)的考量。而這些對于蚓孔方向和形態(tài)會發(fā)生重要改變，且耦合注入酸質(zhì)量分數(shù)的變化，蚓孔的物理性質(zhì)變化將更加復(fù)雜。

為此，現(xiàn)構(gòu)建縫洞的隨機初始化算法，實現(xiàn)縫洞內(nèi)固相含率、滲透率等物理性質(zhì)的真實刻畫，并基于雙尺度連續(xù)模型，研究不同酸濃度、表面化學(xué)反應(yīng)速率以及注入速度等差異下對縫洞型儲層酸蝕過程及蚓孔發(fā)育的影響。試圖厘清縫洞型儲層酸化過程中的物理化學(xué)機理，進而為優(yōu)化改良酸化工藝提供借鑒和指導(dǎo)。

1 數(shù)值模型及算法

針對碳酸鹽巖儲層中酸蝕過程建立了以下數(shù)學(xué)模型，該模型通過Navier-Stokes-Darcy方程模擬碳酸鹽巖內(nèi)流體的流動行為[20]，耦合化學(xué)反應(yīng)方程描述鹽酸對碳酸鹽巖的溶蝕行為，并隨著溶蝕模擬的進行動態(tài)調(diào)整孔隙率等物理量。

1.1 流體動力學(xué)方程

1.1.1 質(zhì)量守恒方程

(1)

式(1)中：Uf為流體的表觀速度；ρf為流體的密度；t為時間；φf為孔隙率。

(2)

式(2)中：V為總體積；Vf為流體的體積；其中，φf=1時，為流體區(qū)域，即巖石完全溶蝕；0<φf<1的位置為多孔區(qū)，即酸液未將巖石完全溶蝕。

1.1.2 動量守恒方程

流體在碳酸鹽巖儲層中的流動遵循達西定律，在孔隙尺度下流體的動量守恒方程為

(3)

需要說明的是，式(3)是孔隙尺度Naiver-Stokes方程的進一步發(fā)展，在考慮流體相的慣性效應(yīng)和黏性耗散外，最后一項考慮了多孔結(jié)構(gòu)的阻尼效應(yīng)。該方程與達西方程相比，可以適用于酸化過程的所有區(qū)域[13]。

1.1.3 酸液輸運方程

酸液在縫洞中流動，伴有擴散和化學(xué)反應(yīng)行為，該行為可以用以下方程來描述：

(4)

式(4)中：Cf為流體區(qū)內(nèi)平均酸濃度；R(Cs)為礦物與酸反應(yīng)過程中單位時間單位體積酸的消耗量；Cs為礦物表面的酸濃度；Deff為酸的有效擴散系數(shù)。酸的濃度在方程中表示為酸的質(zhì)量分數(shù)。反應(yīng)速率可用一級動力學(xué)描述[12, 21]，即

R(Cs)=ksCs

(5)

式(5)中：ks為化學(xué)反應(yīng)速率，為0.015 m/s；Cs為質(zhì)量分數(shù)。酸必須從體相中擴散到界面上才能與固體發(fā)生反應(yīng)，其輸運通量Jacid與傳質(zhì)系數(shù)kc的關(guān)系可以表示為

Jacid=kc(Cf-Cs)

(6)

當(dāng)擴散過程和反應(yīng)過程兩者平衡時，下面關(guān)系式成立：

kcCs=kc(Cf-Cs)

(7)

通過式(7)可得

(8)

將式(8)代入式(5)可得

(9)

(10)

當(dāng)擴散主控(kc?ks)時，有效傳質(zhì)系數(shù)keff≈kc；而當(dāng)反應(yīng)主控(kc?ks)時，keff≈ks。對于傳質(zhì)系數(shù)kc通常通過Sherwood數(shù)Sh計算[22]。

(11)

式(11)中：rp為孔的半徑；Dm為體相擴散系數(shù)，為2.15×10-9m2/s。

(12)

式(12)中：φf 0為初始孔隙度；rp 0為初始孔隙半徑，取為1×10-5m；r為經(jīng)驗參數(shù)，取為3。

方程(11)中，Sh數(shù)可采用式(13)計算：

(13)

式(13)中：Sh為3.66；b為依賴于孔隙結(jié)構(gòu)的參數(shù)，取值為0.7；Sc為施密特數(shù)；雷諾數(shù)Rep求解公式為

(14)

施密特數(shù)Sc求解公式為

(15)

Deff表示酸的有效擴散系數(shù)，計算公式為

Deff=φfDm

(16)

單位體積的礦物表面積αv，公式為

(17)

式(17)中：η為經(jīng)驗參數(shù)，取為1；αv0為初始的礦物表面積，平均值為500 m-1。隨著化學(xué)反應(yīng)的進行，巖心的滲透率[23-25]為

(18)

式(18)中：δ為經(jīng)驗參數(shù)，取為6。巖石孔隙度的更新公式為

(19)

式(19)中：β100為單位質(zhì)量的純酸溶解的礦物的質(zhì)量，取為1.37；ρm為礦物密度，為2 710 kg/m3。

1.2 數(shù)值算法

基于OpenFOAM計算力學(xué)類庫，針對1.1節(jié)中的數(shù)學(xué)模型，時間項采用歐拉格式(Euler)離散，梯度項采用高斯線性差分(Gauss linear)離散，拉普拉斯項采用正交的高斯線性差分(Gauss linear orthogonal)離散。模擬過程中不同物理量的殘差設(shè)置為10-6，最大庫朗數(shù)為0.4，采用自適應(yīng)時間步長。

1.3 數(shù)學(xué)模型驗證

為了驗證本文數(shù)學(xué)模型的準確性，采用Furui等[26]的實驗數(shù)據(jù)進行了驗證。該實驗采用巖心驅(qū)替裝置研究了注入速度對酸液濾失的影響。在2.54 cm寬、15.24 cm長的高孔隙度白堊巖內(nèi)采用質(zhì)量分數(shù)為28%的鹽酸進行驅(qū)替，實驗得到蚓孔突破儲層后得到的無量綱溶蝕蚓孔體積。采用本文所述的數(shù)值模型對該酸化實驗進行模擬，圖1所示為本文的數(shù)學(xué)模型所得結(jié)果與實驗的對比，顯示在不同的注入速度下，數(shù)值模擬能夠量化無量綱溶蝕孔隙體積。

圖1 實驗值與數(shù)值模擬值對比Fig.1 Comparison between experimental data and numerical simulation values

2 模擬與算例設(shè)置

2.1 物理模型及網(wǎng)格劃分

圖2為本文數(shù)值模擬計算所采用的碳酸鹽巖儲層理想物理模型，該圓形物理模型半徑為0.15 m，中心為注入孔，平均孔隙度為0.25，平均滲透率為50 mD。采用四面體結(jié)構(gòu)化均一網(wǎng)格對該物理模型進行網(wǎng)格劃分。整個計算域生成的總網(wǎng)格數(shù)量為196 256，網(wǎng)格尺度為0.000 6 m。

圖2 物理模型Fig.2 Physical model

2.2 邊界條件及數(shù)值算法

2.2.1 邊界條件

數(shù)值模擬計算中，主要物理量及物理邊界對應(yīng)的數(shù)值邊界條件如表1所示。

表1 數(shù)值邊界條件

2.2.2 隨機孔洞生成

為了在計算域網(wǎng)格中表征縫洞不同物理量在計算域中的分布，采用如下算法：①以系統(tǒng)時間為隨機分布量，通過設(shè)置洞和裂縫個數(shù)，隨機生成不同位置的天然裂縫和洞；②在縫洞內(nèi)改變其中的固相分布、滲透率倒數(shù)、礦物表面積等參數(shù)值，滿足數(shù)值模擬過程中縫洞的物理特性。圖3給出了縫洞型碳酸鹽初始固相分布、滲透率倒數(shù)分布及碳酸鹽巖礦物表面積分布云圖。

2.3 算例設(shè)置

本文研究中設(shè)置了如表2所示的8個算例，分別考慮了有無縫洞(算例1和2)、注入酸濃度(算例2～4)、注入速度(算例3、5和6)、及表面化學(xué)反應(yīng)速率(算例3、7和8)4種情況，探討上述因素對縫洞型碳酸鹽巖酸蝕過程的影響。

表2 算例設(shè)置

圖3 縫洞型碳酸鹽巖計算域初始化Fig.3 Initialization of fracture vuggy carbonate calculation domain

3 結(jié)果與分析

3.1 有無縫洞對酸蝕的影響

3.1.1 溶蝕形態(tài)

圖4 碳酸鹽巖儲層有無縫洞下100 s時溶蝕形態(tài)Fig.4 Dissolution form of carbonate reservoir at 100 s with or without fractured-vuggy

圖4為有無縫洞存在下計算時間100 s時的溶蝕形態(tài)。圖4(a)中，當(dāng)碳酸鹽巖儲層不存在縫洞時，在儲層不同方向均發(fā)育了多條主蚓孔，且在主蚓孔上發(fā)育多個長度較短的次生蚓孔。這主要是由于相比于次生蚓孔來說，主蚓孔形成的高滲通道使得酸液流動的阻力降低，從而酸會優(yōu)先進入高滲帶溶蝕碳酸鹽巖，主蚓孔內(nèi)的酸蝕行為形成正反饋并促進蚓孔的進一步發(fā)育。而如圖4(b)所示，溶蝕主要集中在特定的天然縫及洞內(nèi)，形成的蚓孔在不斷發(fā)育的過程中與天然裂縫貫通，并沿著天然裂縫的走向進一步溶蝕發(fā)育。這一過程將導(dǎo)致：①改變主蚓孔的發(fā)育形態(tài)及方向；②改變酸突破碳酸鹽巖儲層的時間。此外，由于酸溶蝕縫洞較快地形成了突破碳酸鹽巖儲層的蚓孔，造成酸流動擴散主要聚集在這一通路上，而在入口周邊區(qū)域，酸與碳酸鹽巖的反應(yīng)不充分，形成的次級蚓孔長度較短。

為了進一步厘清縫洞中的酸蝕過程及蚓孔的形成機制。圖5給出了有縫洞時儲層(算例3)不同時刻溶蝕形態(tài)隨時間變化。圖5(a)中，20 s時刻酸在入口區(qū)域及相鄰的天然裂縫中與碳酸鹽巖儲層發(fā)生化學(xué)反應(yīng)，且溶蝕過程中受天然裂縫的影響，沿著天然裂縫展向上進一步溶蝕。如圖5(b)所示，隨著時間的不斷增大，溶蝕進一步發(fā)展，沿著天然裂縫和洞發(fā)育了貫通儲層的主蚓孔，即形成了特定的溶蝕通路，這也意味著高滲通道的建立。圖5(c)和圖5(d)中，天然裂縫的寬度方向被不斷溶蝕，溶蝕的蚓孔形態(tài)未發(fā)生明顯變化。值得一提的是，蚓孔在發(fā)育的過程中儲層下方的多條天然裂縫雖然對酸存在著一定的導(dǎo)流能力，但上方的通路突破儲層時間較短，造成上方區(qū)域壓力較小、滲透性較高，易于酸的流動，進而不斷被溶蝕。

圖5 碳酸鹽巖縫洞型不同時刻溶蝕形態(tài)Fig.5 Dissolution form of carbonate reservoir at differenttime with fractured-vuggy

3.1.2 溶蝕量

溶蝕量可以通過碳酸鹽巖含量進行表征。圖6所示為有無縫洞存在下的碳酸鹽巖平均固含率隨時間變化。當(dāng)無縫洞時，在儲層各個方向上均發(fā)生溶蝕，固含率下降較快；且在44 s時，由于主蚓孔溶蝕突破儲層，溶蝕過程其后主要發(fā)生在主蚓孔內(nèi)酸流動方向的展向上，這限制了固含率的下降幅度。縫洞的存在會加速溶蝕突破儲層這一過程，表現(xiàn)為在0～20 s時刻突破儲層前固相體積分數(shù)下降較快，而隨后下降較慢。這是由于：①縫洞的存在縮短了溶蝕貫通儲層的時間；②溶蝕的范圍圍繞酸在裂縫或洞流動的展向區(qū)域。此外，對比有無縫洞時儲層0 s和100 s的固相體積分數(shù)，顯示縫洞存在下儲層溶蝕量下降幅度為6.3%低于無縫洞時的9.5%。

圖6 有無縫洞碳酸鹽巖儲層內(nèi)固相體積分數(shù)隨時間變化Fig.6 Variation of solid volume fraction in acid salt reservoir with or without fractured-vuggy

3.1.3 壓力差

圖7 有無縫洞出入口壓差隨時間變化Fig.7 Variation of pressure in acid salt reservoir with or without fractured-vuggy

出入口壓力差對于判斷溶蝕動態(tài)過程具有重要意義。如圖7所示為有無縫洞時監(jiān)測得到的出入口壓差隨時間變化。顯示無縫洞時，注入酸后出入口壓差峰值明顯大于有縫洞時；且隨著時間的增大，酸不斷溶蝕碳酸鹽巖降低了儲層整體等效滲流阻力，導(dǎo)致壓力不斷下降至穩(wěn)定值(稍大于0 Pa)。而存在縫洞時，在溶蝕過程中連接縫和洞時會造成出入口壓差明顯的波動，且蚓孔突破儲層后出入口壓差達到穩(wěn)定值的時間明顯較短。

3.2 不同因素對縫洞型碳酸鹽巖儲層酸化過程影響

3.2.1 不同注入酸濃度對酸化過程影響

(1)溶蝕形態(tài)。如圖8所示，為100 s時刻不同注入酸濃度下縫洞型碳酸鹽巖的溶蝕形態(tài)。由圖8可知，隨著注入酸質(zhì)量分數(shù)的增大，被溶蝕的天然裂縫寬度明顯增大，但形成的蚓孔走向和形態(tài)相同。其原因是：①酸質(zhì)量分數(shù)的增高加大了參與化學(xué)反應(yīng)的氫離子濃度，進而蚓孔的寬度不斷增大,但受酸傳質(zhì)效應(yīng)的主導(dǎo)，化學(xué)反應(yīng)速率提高不明顯，這也是蚓孔寬度增大幅度不大的主要原因；②由于天然裂縫和洞對酸液的導(dǎo)流效應(yīng)，酸質(zhì)量分數(shù)的增大并未改變蚓孔的走向。

圖8 不同酸濃度下100 s時溶蝕形態(tài)Fig.8 Dissolution form of carbonate reservoir at 100 s with fractured-vuggy under different acid concentrations

(2)溶蝕量。圖9給出了不同酸注入濃度下固相體積分數(shù)隨時間變化?？p洞型儲層內(nèi)，注入酸質(zhì)量分數(shù)越大，不同時間下，平均固相體積分數(shù)越小。隨著酸濃度的不斷增大，固相體積分數(shù)的下降幅度不斷增大，但曲線仍然呈現(xiàn)先快速降低后緩慢降低的特點。且對比0 s以及100 s時刻的固相溶蝕量發(fā)現(xiàn)，隨著酸質(zhì)量濃度的不斷增大溶蝕量不斷增大，分別為5.7%、6.3%和7.0%。

圖9 不同酸濃度下固相體積分數(shù)隨時間變化Fig.9 Variation of solid volume fraction with time at different acid concentrations

(3)壓力差。圖10給出了不同酸濃度下出入口壓差隨時間變化。縫洞型碳酸鹽巖內(nèi)，中高濃度酸對碳酸鹽巖儲層溶蝕導(dǎo)致的貫通時間較為接近，低濃度酸對碳酸鹽巖儲層溶蝕導(dǎo)致的貫通時間較長，但均明顯小于無縫洞時碳酸鹽巖的44 s。這意味著酸濃度的提高對于蚓孔在長度上的發(fā)育影響較小。值得一提的是，出入口壓差較大的波動往往表現(xiàn)為天然縫洞間的融通，不同酸濃度作用下，出入口壓差的形態(tài)大致相近，這一定程度反映出蚓孔的形成路徑大致相同，僅在形成時間上存在一定差異。

圖10 不同酸濃度下出入口壓差隨時間變化Fig.10 Variation of inlet and outlet pressure difference with time under different acid concentration

3.2.2 不同注入速度對酸化過程影響

(1)溶蝕形態(tài)。圖11所示為不同注入速度下無量綱時間為100 s時的溶蝕形態(tài)。對注入時間采用公式進行了無量綱化，保證無量綱時間下注入量相同。

(20)

式(20)中：T為無量綱時間；u為注入速度；t為物理時間；u0和t0分別為0.003 m/s和1 s。

隨著注入速度的減小，次生蚓孔的數(shù)量明顯降低，意味著低注入速度下，酸巖反應(yīng)更加充分。圖11(a)中，注入速度最高，酸鹽反應(yīng)在次生蚓孔周圍區(qū)域還未充分反應(yīng)就進入儲層深部，存在著一定的未充分反應(yīng)完全區(qū)域。圖11(b)和圖11(c)中，當(dāng)注入速度減小為0.000 6、0.000 12 m/s時，雖然碳酸鹽巖儲層上部的天然裂縫和洞被溶通；但下部的儲層溶蝕形態(tài)有著一定區(qū)別，這可能是由于低速下溶蝕會在注入端進行面溶蝕，這也與廖毅等[17]的認識相一致。

圖11 不同注入速度無量綱時間100 s時溶蝕形態(tài)Fig.11 Dissolution form of carbonate reservoir at 100 s dimensionless time with fractured-vuggy under different injection rates

(2)溶蝕量。圖12所示為不同注入速度下固相體積分數(shù)隨時間變化。在注入速度較高時(0.003 m/s和 0.000 6 m/s)，固相體積分數(shù)隨時間變化開始較大，后逐漸減小，且兩者數(shù)值較為接近。而注入速度最低時(0.000 12 m/s)，固相體積分數(shù)隨時間變化更加接近線性，且固相體積分數(shù)在無量綱時間為0～100 s時下降幅度為5.7%，明顯高于算例3的6.3%和算例5的6.4%。這也與圖11(c)中無量綱時間100 s時的溶蝕形態(tài)相呼應(yīng)。

圖12 不同注入速度下固相體積分數(shù)隨無量綱時間變化Fig.12 Variation of solid volume fraction with dimensionless time under different injection rates

圖13 不同注入速度下出入口壓差隨無量綱時間變化Fig.13 Variation of inlet and outlet pressure difference with dimensionless time under different injection speeds

(3)壓力差。圖13給出了不同注入速度下出入口壓差隨無量綱時間變化。顯示注入速度的改變能夠明顯地降低整個碳酸鹽巖儲層內(nèi)的壓力水平。當(dāng)注入速度減小時，出入口壓差的峰值降低；出入口壓差達到平衡值(約0 Pa)的時間增大，意味著酸蝕過程突破儲層時間增大。此外，不同注入速度下，酸蝕過程形成的蚓孔方向和尺寸不同，出入口壓差隨時間變化呈現(xiàn)的波動存在著差異。

3.2.3 不同化學(xué)反應(yīng)速率對酸化過程影響

(1)溶蝕形態(tài)?？偡磻?yīng)速率是由酸通過對流或擴散方式運動到礦物表面的傳質(zhì)速率和碳酸鹽巖礦物表面的化學(xué)反應(yīng)速率兩者共同決定的。因此碳酸鹽巖與酸的反應(yīng)類型可以分為：傳質(zhì)控制反應(yīng)型、表面反應(yīng)速度控制反應(yīng)型或兩者共同控制類型。圖14給出了不同酸鹽化學(xué)反應(yīng)速率下100 s時刻碳酸鹽巖的溶蝕形態(tài)。當(dāng)化學(xué)反應(yīng)速率高中低變化時，固相體積分數(shù)的分布較為近似，僅在次生蚓孔的發(fā)育上存在著一定差別。這可能是由于碳酸鹽巖儲層內(nèi)的酸鹽反應(yīng)為傳質(zhì)控制反應(yīng)類型，此時化學(xué)反應(yīng)速率的差異對于溶蝕的影響較小。

圖14 不同化學(xué)反應(yīng)速率下100 s時溶蝕形態(tài)Fig.14 Dissolution form of carbonate reservoir at 100 s with fractured-vuggy under different chemical reaction rates

(2)溶蝕量。圖15給出了不同化學(xué)反應(yīng)速率下碳酸鹽巖固含量隨時間變化。在中高化學(xué)反應(yīng)速率下，固相體積分數(shù)的變化較為一致，且在碳酸鹽巖儲層突破后下降幅度變緩。隨著化學(xué)反應(yīng)速率的增大，在0～100 s時刻內(nèi)，碳酸鹽巖總含量下降幅度分別為7%、6.3%、6.5%。中高化學(xué)反應(yīng)速率下的固體含量較為接近，這可能是由于傳質(zhì)效應(yīng)主導(dǎo)總化學(xué)反應(yīng)速率造成的。這也與一般施工條件下酸巖反應(yīng)動力學(xué)幾乎總是受擴散傳質(zhì)控制的認識一致[27]。而在低化學(xué)反應(yīng)速率下，固相體積分數(shù)下降幅度反而大，其原因可能是由于傳質(zhì)效應(yīng)與化學(xué)反應(yīng)效應(yīng)共同控制總反應(yīng)速率，導(dǎo)致總反應(yīng)速率相對較高。

(3)壓力差。圖16給出了不同化學(xué)反應(yīng)速率下出入口壓差隨時間變化。顯示隨著化學(xué)反應(yīng)速率的增大，儲層突破時間越大，但突破時間差異不顯著，均在20 s左右；且出入口壓差隨時間的變化也未呈現(xiàn)出明顯差異，變化趨勢和波動趨勢近似。這也意味著，在特定速度下，酸鹽反應(yīng)速率的大小對于形成蚓孔進而貫通儲層沒有明顯的影響。

圖15 不同化學(xué)反應(yīng)速率下固相體積分數(shù)隨時間變化Fig.15 Variation of solid volume fraction with time at different chemical reaction rates

圖16 不同化學(xué)反應(yīng)速率下出入口壓差隨時間變化Fig.16 Variation of inlet and outlet pressure difference with time under different chemical reaction rates

4 結(jié)論

(1)縫洞型油藏的酸蝕過程蚓孔發(fā)育受縫和洞共同影響，不斷經(jīng)歷沿著縫隙溶蝕后受壓力的牽引溶蝕新縫隙或洞的過程。采用普通油藏進行酸化數(shù)值模擬的研究會高估對縫洞型油藏酸化過程中的溶蝕量?？p洞的存在會導(dǎo)致溶蝕過程中壓力的波動。

(2)隨著酸濃度的升高，固相溶蝕量由5.7%增大為7.0%，儲層的貫通時間不斷減小；在酸流經(jīng)縫洞的展向上蚓孔的面積增大，但對于裂縫型油藏的蚓孔形態(tài)不會形成明顯的改變。

(3)注入速度越大，無量綱時間下酸蝕影響的區(qū)域越大；低注入速度下，會減弱次生蚓孔的發(fā)育，貫通裂縫的數(shù)量減少。在特定速度下，酸鹽反應(yīng)速率的大小對于溶蝕形態(tài)及蚓孔與儲層的貫通時間沒有明顯的影響。