邊界層對致密油藏逆向滲吸的數(shù)值模擬

徐中一，方思冬

(中國石化 石油勘探開發(fā)研究院，北京 100083)

目前致密油藏開發(fā)過程中的滲吸傳質(zhì)作用十分復(fù)雜，其復(fù)雜性主要集中在兩個(gè)方面，第一由于縫網(wǎng)的存在大幅度增加了基質(zhì)-水的接觸面積，導(dǎo)致滲吸的范圍增加；第二因?yàn)橹旅軆拥暮淼兰?xì)小達(dá)到了微納米級別[1-4]，已經(jīng)達(dá)到了流體分子尺度的級別，因此微納米喉道內(nèi)的微尺度效應(yīng)將十分明顯，在微尺度效應(yīng)的影響下普通油藏適用的達(dá)西定律在致密油藏中可能不再適用。

為了準(zhǔn)確模擬致密儲層滲吸流動，本文首先基于耗散粒子動力學(xué)(DPD)方法在微納米喉道中模擬了油水兩相的滲吸流動過程并證明了微納米喉道中存在邊界層效應(yīng)，通過調(diào)研現(xiàn)有的文獻(xiàn)以及DPD模擬結(jié)果綜合表征了邊界層厚度的變化規(guī)律，在此基礎(chǔ)上利用孔隙網(wǎng)絡(luò)模型計(jì)算致密儲層的有效滲透率、相滲曲線及毛管力曲線；同時(shí)，利用已有的滲吸速度計(jì)算模型結(jié)合得到的關(guān)鍵滲流參數(shù)，得到關(guān)于致密儲層的滲吸速度的變化規(guī)律，將實(shí)驗(yàn)結(jié)果與計(jì)算結(jié)果對比驗(yàn)證本文提出的關(guān)鍵滲流參數(shù)的計(jì)算方法是否正確，最后，利用數(shù)值模擬方法闡明致密儲層中考慮邊界層的滲吸對開發(fā)效果的影響。

1 致密儲層的微觀流動

1.1 滲吸過程中邊界層現(xiàn)象兩相的情況

采用耗散粒子動力學(xué)(DPD)方法對納米喉道內(nèi)的滲吸現(xiàn)象進(jìn)行了模擬，整個(gè)模型由3種珠子構(gòu)成，油珠子、水珠子和二氧化硅珠子，其中二氧化硅珠子用作構(gòu)建石英微圓管，二氧化硅分子被分為兩類，第一類SiO2珠子用于構(gòu)建石英微圓管內(nèi)部，第二類SiO2珠子用于構(gòu)建石英微圓管的管壁，兩種SiO2粒子的受力不同，第一類石英微圓管內(nèi)部的SiO2珠子受力主要是SiO2珠子間的結(jié)構(gòu)力，第二類石英管表面的SiO2除了與第一類SiO2珠子間存在結(jié)構(gòu)力外還有與流體間存在相互作用力，微圓管內(nèi)壁表面的SiO2作為非極性分子與水中的氫離子結(jié)合形成羥基，帶有羥基的SiO2珠子與H2O珠子之間會形成氫鍵，導(dǎo)致了壁面對水的吸引[5]。將二氧化硅珠子均勻排列構(gòu)建微圓管并作為流體的流動通道，油珠子和水珠子分別構(gòu)成了滲吸過程中的油水兩相，粒子的基本參數(shù)如表1所示，在表1中的amu為原子質(zhì)量單位，它的定義為碳12元素原子質(zhì)量的1/12;?為半徑的長度單位埃，?比納米小一個(gè)數(shù)量級，1埃等于0.1nm。

接下來只要準(zhǔn)確構(gòu)建粒子間相互作用力就可以準(zhǔn)確模擬微納米喉道中的滲吸流動，根據(jù)DPD理論珠子間的相互作用力主要有3種[6-7]：保守力、耗散力和隨機(jī)力，其中耗散力的設(shè)置依據(jù)田虓豐博士論文的設(shè)置，耗散力半徑為1.5(12?)時(shí)，耗散力系數(shù)為45(0.885 42 amu/fs)時(shí)，密度為3(1g/cm3)的流體粘度為1 mPa·s[7]。此時(shí)當(dāng)體系內(nèi)有油相后，粒子的密度為5(1.6 g/cm3)時(shí)[7]，將粒子的截?cái)喟霃秸{(diào)大，設(shè)置為20?時(shí)，此時(shí)流體粘度的粘度增加為3 mPa·s。而隨機(jī)力和耗散力有直接的關(guān)系，當(dāng)耗散力確定后，軟件會根據(jù)耗散力的設(shè)置自動設(shè)置隨機(jī)力[5]，保守力參數(shù)的計(jì)算根據(jù)Flory-Huggins理論,分子間的保守力計(jì)算如式(1)和式(2)[8-9]：

(1)

其中χ為Flory-Huggins相互作用參數(shù)，F(xiàn)lory-Huggins的相互作用參數(shù)一般通過式(2)得到[11]：

(2)

式中：χ為Flory-Huggins相互作用參數(shù)，J/(mol·K)；V為粒子平均摩爾體積，cm3/mol；R為氣體常數(shù)，8.314 J/(mol·K)；T為溫度,K;δA,δB為粒子的溶解度參數(shù)，(J/cm3)0.5。

利用Material Studio(MS) 軟件中的Forcite模塊計(jì)算溶解度參數(shù)，將計(jì)算獲得的溶解度參數(shù)輸入MS軟件的DPD模塊獲得最終的各珠子間的保守力系數(shù)如表2所示。

另外為了準(zhǔn)確的描述管壁親水作用，需要在普通耗散粒子動力學(xué)方法的基礎(chǔ)上添加管壁SiO2珠子和水珠子間額外的吸引力，根據(jù)實(shí)驗(yàn)室測得的邊界層厚度變化規(guī)律以及極限邊界層厚度[12]，選用MS軟件所提供的純吸引力函數(shù)LJ96X[如式(3)]，且其參數(shù)為D0=1，R0=1，E=1時(shí)，其擬合得到的邊界層厚度的變化規(guī)律與實(shí)際相符。

(3)

表1 DPD粒子基本參數(shù)Table 1 Basic parameters of DPD particles

表2 粒子保守力系數(shù)Table 2 Particle conservative force coefficients

式中：F為吸引力，N；D0為吸引力參數(shù)，kcal/mol；R為粒子間距離，?；R0為粒子間平衡距離，?；E為偏移距離，?。

將表2中計(jì)算得到的粒子間作用力參數(shù)以及LJ96X吸引力函數(shù)輸入MS軟件中的DPD模塊，進(jìn)行模擬計(jì)算，計(jì)算結(jié)果如圖1所示，圖1a顯示的是滲吸模擬過程中珠子分布圖，圖1b顯示的是水珠子的密度分布圖；從圖中可以看出，在滲吸階段水珠子在管壁珠子的吸引力作用下進(jìn)入微圓管壁，并在管壁處形成了邊界層，遠(yuǎn)離管壁處的水珠子則形成彎液面向管內(nèi)推進(jìn)。

1.2 模型的抽提及邊界層的定量表征

根據(jù)DPD研究顯示的邊界層現(xiàn)象，抽提得到致密儲層滲吸時(shí)的模型(圖2)，模型顯示當(dāng)考慮致密多孔介質(zhì)喉道為4～200 nm時(shí)，邊界層現(xiàn)象都十分明顯。根據(jù)鄒才能對致密儲集層特征的描述，致密油儲集層孔喉直徑主要分布在50～800 nm,孔隙度小于10%，滲透率為1×10-6～1×10-3μm2 [10]。根據(jù)上述研究，在致密儲層微納米喉道中，不能忽視邊界的存在；邊界層對致密儲層中滲吸過程的影響如圖2所示，微納米喉道中的邊界層會減小有效流動空間，如圖2b所示rori為致密儲層多孔介質(zhì)的原始喉道半徑，邊界層被認(rèn)為是靜止的不參與喉道內(nèi)流體的流動，原始喉道半徑減去邊界層厚度h后則為喉道內(nèi)流體的有效流動半徑reff，相比于中低滲油藏，致密儲層由于本身喉道細(xì)小加上邊界層的作用導(dǎo)致有效流動半徑的進(jìn)一步減小。

為了定量表征邊界層對流動的影響，根據(jù)本文DPD模擬結(jié)果以及李洋等人實(shí)驗(yàn)結(jié)果[12]，邊界層厚度可以定量表征為公式(4)，公式(4)具有一定的適用范圍，應(yīng)用時(shí)喉道壁面需完全潤濕，壓力梯度不得小于0.02 MPa/m，同時(shí)喉道半徑需小于0.2 μm，管內(nèi)的流體需是牛頓流體[9]。

(4)

式中：r為孔隙的喉道半徑,μm；h為孔隙中的邊界層厚度,μm；μ為流體粘度,mPa·s；Δp為壓力梯度，MPa/m。

2 致密儲層關(guān)鍵參數(shù)的計(jì)算

根據(jù)徐中一關(guān)于滲吸的文章,準(zhǔn)確計(jì)算滲吸量需要準(zhǔn)確計(jì)算滲透率、油水相滲、毛管力曲線3個(gè)關(guān)鍵滲流參數(shù)[9]，本文使用孔隙網(wǎng)絡(luò)模型，在流動模擬的過程中考慮邊界層對流動的影響，直接計(jì)算以上3個(gè)關(guān)鍵滲流參數(shù)。

圖1 耗散粒子動力學(xué)(DPD)模擬結(jié)果Fig.1 Diagram of Dissipative Particle Dynamics (DPD) simulationa.粒子圖；b.流體密度圖

圖2 基于DPD實(shí)驗(yàn)建立的考慮邊界層的兩相流動模型Fig.2 A two-phase flow model considering boundary layers based on DPD experimentsa.不考慮邊界層時(shí)的滲吸情況；b.考慮邊界層的滲吸情況

2.1 滲透率的計(jì)算

依據(jù)核磁公振實(shí)驗(yàn)對實(shí)際致密巖心在飽和水狀態(tài)下測得的T2譜經(jīng)過轉(zhuǎn)化后獲得的實(shí)際致密多孔介質(zhì)中孔隙大小的分布，如圖3a所示，將致密儲層的多孔介質(zhì)處理為三維的隨機(jī)網(wǎng)絡(luò)模型，將多孔介質(zhì)中較大的空間簡化為球型孔隙，較小的空間簡化為圓柱狀喉道，節(jié)點(diǎn)代表孔隙在三維空間的位置，節(jié)點(diǎn)間的連線代表孔隙間的喉道，其中無論是孔隙和喉道都要滿足致密儲層多孔介質(zhì)內(nèi)孔隙和喉道的分布，最后形成的孔隙網(wǎng)絡(luò)模型如圖3所示。孔隙和喉道大小都滿足截?cái)喔咚鬼f布爾分布[公式(5)][13-15]，通過調(diào)節(jié)參數(shù)rmax、rmin使得孔隙網(wǎng)絡(luò)模型中的孔隙分布范圍與圖3a實(shí)測的孔隙分布范圍相符。

r=(rmax-rmin)(-δln[x(1-e-1/δ)+e1/δ])1/γ+rmin

(5)

式中：r為喉道半徑，μm；rmax為最大喉道半徑，μm；rmin為最小喉道半徑，μm；γ、δ為幾何分布參數(shù)，無因次；γ=-1，δ=-1使得喉道和孔隙規(guī)則分布；x為隨機(jī)變量，取值范圍0～1。

計(jì)算滲透率時(shí)，首先在孔隙網(wǎng)絡(luò)兩側(cè)，流體的入口和流體的出口分別設(shè)置恒定的壓力，使得孔隙網(wǎng)絡(luò)中的流體在恒定的壓力梯度下流動，孔隙網(wǎng)絡(luò)模型內(nèi)部的情況如圖4a所示，孔隙間喉道中的流體在流動的過程中存在著邊界層，邊界層的厚度根據(jù)式(4)計(jì)算；在計(jì)算過程中首先計(jì)算每個(gè)孔隙處的壓力，依據(jù)每個(gè)喉道相連的孔隙的壓力與喉道本身的長度計(jì)算喉道內(nèi)的壓力梯度，從而計(jì)算每個(gè)喉道內(nèi)的邊界層的厚度，因?yàn)檫吔鐚邮遣粎⑴c流動的所以邊界層的存在就減小了喉道內(nèi)流體的有效流動空間。最后在孔隙網(wǎng)絡(luò)模型兩端(如圖4a所示的入口和出口)設(shè)置不同的壓差,每設(shè)置一個(gè)壓差就計(jì)量孔隙網(wǎng)絡(luò)模型末端的流量，最后用公式(6)計(jì)算致密儲層的滲透率。

圖3 構(gòu)建的孔隙網(wǎng)絡(luò)模型Fig.3 Pore network modela.孔隙大小分布；b.依據(jù)孔喉分布建立的20*20*20的孔隙網(wǎng)絡(luò)模型

圖4 二維孔隙網(wǎng)絡(luò)模型中驅(qū)替和滲吸的情況Fig.4 Schemes showing displacement and imbibition in a two-dimension pore network modela.驅(qū)替的情況；b.滲吸的情況

(6)

式中：k為滲透率，10-3μm2；q為出口端的流量，cm3/s；μ為流體的粘度，mPa·s；L為孔隙網(wǎng)絡(luò)模型的長度，cm；Δp為孔隙網(wǎng)絡(luò)兩側(cè)的壓力差，MPa。

通過以上方法計(jì)算出的滲透率如圖5所示，可以看出當(dāng)不考慮邊界層的多孔介質(zhì)而言，滲透率是一個(gè)僅與多孔介質(zhì)有關(guān)的參數(shù)，只要多孔介質(zhì)一定無論怎樣改變壓力梯度，滲透率的大小都不變；考慮邊界層后滲透率就不僅與多孔介質(zhì)本身有關(guān)而是與流體流動有關(guān)，從圖5可以看出當(dāng)壓差較小時(shí)，邊界層厚度較大，有效流動空間較小，因此有效滲透率偏離不考慮邊界層時(shí)的滲透率，而隨著壓力梯度的增大，邊界層厚度逐漸減小，直到壓力梯度很大達(dá)到104MPa/m時(shí)，邊界層厚度遠(yuǎn)小于喉道半徑，此時(shí)邊界層的可以忽略。因此隨著壓力梯度的增大，致密多孔介質(zhì)的滲透率逐漸接近不考慮邊界層的多孔介質(zhì)的滲透率。

2.2 相對滲透率與毛管力曲線的計(jì)算

滲吸的過程是一個(gè)兩相流動的過程，因此需要計(jì)算致密多孔介質(zhì)滲吸過程的油水兩相相滲曲線和毛管力曲線，本文同樣采用孔隙網(wǎng)絡(luò)模型，基于Blunt公布的流動模擬程序[16]，在其程序的基礎(chǔ)上增加了邊界層計(jì)算模塊，核心的計(jì)算步驟如圖6所示，主要需要計(jì)算每個(gè)孔隙處的壓力值，然后根據(jù)壓力值和喉道的長度計(jì)算每個(gè)喉道的壓力梯度再根據(jù)壓力梯度計(jì)算每個(gè)喉道的邊界層厚度，值得注意的是，如圖4b所示只有水濕的喉道需要計(jì)算邊界厚度，油濕的喉道是不存在邊界層的；計(jì)算完邊界層后再計(jì)算由于邊界層的變化導(dǎo)致的喉道的傳導(dǎo)系數(shù)的變化。

滲吸的過程不同于計(jì)算滲透率的驅(qū)替過程，是在入口和出口設(shè)定恒定的壓力；滲吸的邊界條件是在入口設(shè)置恒定的含水飽和度為1，出口端設(shè)置恒定的壓力，本文設(shè)置了大氣壓力0.101 MPa，初始時(shí)刻每個(gè)孔隙處的壓力也是大氣壓力，計(jì)算得到的相滲曲線和毛管力曲線如圖7和圖8所示。

計(jì)算結(jié)果表明邊界層對于致密儲層的相對滲透率曲線和毛管力曲線有較大的影響，從計(jì)算得到的相滲曲線圖7a來看，邊界層的存在主要減弱了油水兩相在多孔介質(zhì)中的流動能力，同時(shí)邊界層的存在會使得油水兩相區(qū)縮小，束縛水飽和度和殘余油飽和度都增加，這是因?yàn)檫吔鐚拥拇嬖谑沟迷絹碓蕉嗟挠捅弧版i”在基質(zhì)中無法通過滲吸作用流出多孔介質(zhì)。

圖5 考慮邊界層后滲透率的變化Fig.5 A plot showing the permeability variation before and after considering boundary layers

圖6 技術(shù)路線Fig.6 The workflow in technology

從毛管力曲線圖7b來看，考慮邊界層后其毛管力曲線大小明顯較不考慮邊界層的情況要大，這主要是邊界層的存在減小了喉道的有效流動空間，根據(jù)公式pc=2δcosθ/r可知當(dāng)有效流動半徑減小時(shí)，毛管力會增大；同時(shí)毛管力開始發(fā)生上翹的點(diǎn)發(fā)生了右移，這說明很多原先潤濕相可以進(jìn)入的喉道和孔隙現(xiàn)在因?yàn)檫吔鐚拥拇嬖诂F(xiàn)在無法進(jìn)入，因此總的來說毛管力的存在導(dǎo)致毛管力大小增加，另外可流動的喉道和孔隙數(shù)減少。

3 考慮邊界層情況下傳質(zhì)計(jì)算及驗(yàn)證

為了驗(yàn)證以上計(jì)算結(jié)果的正確性，本文挑選了致密巖心進(jìn)行了自發(fā)逆向滲吸實(shí)驗(yàn)，并將實(shí)驗(yàn)測得的滲吸采收率曲線與計(jì)算得到的滲吸采收率進(jìn)行對比，實(shí)驗(yàn)首先將巖心在真空泵中抽真空后飽和配置的重水，然后放入核磁共振儀中測量T2譜，并利用油驅(qū)水實(shí)驗(yàn)制造束縛水，最后取出巖心并除去巖心表面的浮油，最后把巖心放入盛滿水的燒杯中進(jìn)行滲吸實(shí)驗(yàn)(圖8)，實(shí)驗(yàn)過程采用稱重法對滲吸出的油量進(jìn)行稱重，滲吸實(shí)驗(yàn)共持續(xù)6d左右達(dá)到最終采收率，實(shí)驗(yàn)的巖心的基本參數(shù)列于表3中。

圖7 滲流參數(shù)計(jì)算結(jié)果Fig.7 Calculated results of seepage parametersa.相對滲透率曲線;b.毛管力曲線

本文第二節(jié)的孔隙網(wǎng)絡(luò)模型的孔喉分布是基于本次核磁公振測得的T2譜轉(zhuǎn)化獲得的，另外第二節(jié)計(jì)算的滲透率、相對滲透率、毛管力曲線就是基于本次實(shí)驗(yàn)的巖心測得的孔喉分布建立的孔隙網(wǎng)絡(luò)模型計(jì)算的，將計(jì)算得到的滲透率曲線、相對滲透率曲線、毛管力曲線綜合起來利用RIBEM方法[17]計(jì)算滲吸量，結(jié)果如圖9所示，考慮邊界層后滲吸速度和滲吸的最終采收率都大幅降低，當(dāng)不考慮邊界層時(shí)計(jì)算得到的采收率為18%，達(dá)到最終采收率所需的時(shí)間為40 h，而考慮邊界層時(shí)最終采收率僅為10%，達(dá)到最終采收率所需的時(shí)間為6 d左右，說明致密巖心的滲吸速度遠(yuǎn)小于不考慮邊界層的情況。

4 油藏尺度數(shù)值模擬計(jì)算

最后為了研究在油藏尺度下滲吸對致密油藏開發(fā)效果的影響，本文采用非結(jié)構(gòu)PEBI網(wǎng)格貼合人工水力壓裂縫對致密儲層中水平井周圍的縫網(wǎng)進(jìn)行網(wǎng)格剖分如圖10所示，圖中人工裂縫顯式描述，誘導(dǎo)縫和天然裂縫采用雙重介質(zhì)模型描述[18-23],將本文第二節(jié)計(jì)算得到的滲流參數(shù)應(yīng)用到當(dāng)前的油水兩相數(shù)值模擬計(jì)算中，通過對比不同時(shí)刻油藏的含水飽和度和產(chǎn)量，評價(jià)邊界層對致密儲層開發(fā)效果的影響。

圖8 滲吸實(shí)驗(yàn)巖心示意圖Fig.8 The schematic diagram showing the core imbibition experiments

表3 巖心基本參數(shù)Table 3 Basic parameters of core

圖9 計(jì)算滲吸采收率與實(shí)驗(yàn)結(jié)果對比Fig.9 The comparison between the calculated and measured imbibition recovery

圖10 考慮縫網(wǎng)的PEBI網(wǎng)格模型Fig.10 A PEBI fracture network model

本部分模擬了水平井的注水吞吐開采，吞吐一共持續(xù)五個(gè)輪次，每個(gè)輪次持續(xù)50 d，每個(gè)輪次內(nèi)前10 d注水，燜井20 d，后20 d開井生產(chǎn)。與模擬有關(guān)的參數(shù)如表4所示，第一輪次、第三輪次及第五輪次的模擬結(jié)果如圖11及圖12所示，圖中主要展示了不同輪次燜井末期的含水飽和度分布。其中圖11展示了不考慮邊界層情況下的含水飽和度分布，圖12展示了考慮邊界層情況下的含水飽和度分布，油藏的基本參數(shù)如表4所示。

表4 數(shù)值模擬基本參數(shù)Table 4 Basic parameters of numerical simulation

圖11 不同輪次致密油藏考慮邊界層情況下含水飽和度分布對比Fig.11 The comparison of water saturation variation in tight reservoirs considering boundary layers at different water-injection roundsa.第一輪次；b.第三輪次； c.第五輪次

圖12 不同輪次致密油藏不考慮邊界層情況下含水飽和度分布對比Fig.12 The comparison of water saturation distribution in tight reservoirs without considering boundary layers at different water-injection roundsa.第一輪次； b.第三輪次； c.第五輪次

圖13 生產(chǎn)情況對比Fig.13 The comparison of oil production with and without considering boundary layersa.采油速度對比;b.累計(jì)采油量對比

從數(shù)值模擬橫向?qū)Ρ鹊慕Y(jié)果來看，將第五輪次與第一輪次和第三輪次進(jìn)行對比可以看出滲吸的吸水排油作用對壓裂水平井周圍含水飽和度具有一定影響，當(dāng)考慮邊界層時(shí)，由于其滲吸能力較弱(吸水排油的能力弱)，吞吐周期內(nèi)注入的水更多的滯留在裂縫介質(zhì)內(nèi)，導(dǎo)致壓裂水平井周圍的含水飽和度較高。

從數(shù)值模擬橫向?qū)Ρ鹊慕Y(jié)果來看，隨著輪次的增加，注入水波及的范圍逐漸擴(kuò)大，這意味著原先沒有被注入水波及的不發(fā)生滲吸的區(qū)域，逐漸被注入水波及開始發(fā)生滲吸，也就是說，隨著注入水波及的區(qū)域越來越大，發(fā)生滲吸的區(qū)域也越來越大；邊界層的存在會影響注入水向外波及的速度，考慮邊界層時(shí)，滲吸能力弱，更多的水滯留在裂縫中，有利于向外波及，因此當(dāng)考慮邊界層時(shí)，輪次增加更多的水能夠向外波及從而導(dǎo)致能夠發(fā)生滲吸的區(qū)域比不考慮邊界層時(shí)的區(qū)域大。

壓裂水平井的采油速度和累積采油量分別如圖13a和圖13b所示，從圖13a的采油速度可以看出，不考慮邊界層的致密儲層，其采油速度要高于考慮邊界層的儲層，相應(yīng)的考慮邊界層情況下的累積采油量大約是考慮邊界層情況累積采油量的3/8。造成這種現(xiàn)象的主要原因是考慮邊界層后，儲層的滲吸能力減弱，吸水排油的效果減弱，滯留在有效流動介質(zhì)中的水多油少，從而導(dǎo)致最終開采過程中采油速度的減小，另外隨著輪次的增加，地層能量不斷釋放，導(dǎo)致每一輪次的采油能力都在衰減，如圖13a所示，采油速度在逐輪次的遞減。

5 結(jié)論

1) 使用耗散粒子動力學(xué)方法探索發(fā)現(xiàn)了微納米喉道中油水兩相滲吸過程中的邊界層現(xiàn)象，并且得到了邊界層變化的規(guī)律。

2) 使用孔隙網(wǎng)絡(luò)模型計(jì)算了存在邊界層情況下的關(guān)鍵滲流參數(shù)，從計(jì)算結(jié)果來看，對致密儲層來講邊界層的存在會減小儲層的滲透率，邊界層對滲透率的減小程度與壓力梯度的大小有關(guān)，壓力梯度越大邊界層厚度越小，邊界層對滲透率的傷害越小，當(dāng)壓力梯度達(dá)到104MPa/m時(shí)，邊界層的存在幾乎不會對儲層造成傷害。

3) 從油藏尺度來看，致密儲層中存在的邊界層會減弱儲層的流動能力以及基質(zhì)塊的滲吸能力，如果忽略邊界層會過高估計(jì)儲層的生產(chǎn)能力，實(shí)例計(jì)算表明，考慮邊界層后累計(jì)采油量約為不考慮邊界層情況的3/8，因此真實(shí)致密儲層中采用注水吞吐以期利用滲吸作用提高采收率效果應(yīng)是十分有限的。