虞紹永

（達(dá)坦能源集團，上海 200433）

目前，超低滲透儲層（頁巖/致密砂巖）普遍采用水平井分段壓裂技術(shù)開發(fā)，分段壓裂作業(yè)完成后，需要進(jìn)行壓裂液返排，返排時間一般為15～30 d。壓裂液返排制度（包括返排時間、返排期油嘴變化和返排速度）在一定程度上會影響裂縫導(dǎo)流能力，從而影響水平井產(chǎn)量。過去人們對壓裂液返排歷史數(shù)據(jù)的重視度和利用率不高，實際上，這些數(shù)據(jù)對描述水力裂縫特征、深入分析壓裂作業(yè)效果及產(chǎn)量預(yù)測具有重要作用，通過分析這些數(shù)據(jù)，可以得到初始有效裂縫（水力裂縫與天然裂縫）孔隙體積、返排過程中的裂縫體積損失等參數(shù)，進(jìn)而分析評價返排制度對壓裂效果的影響。同時，井下測量技術(shù)的發(fā)展也極大提高了返排時數(shù)據(jù)采集的密度和精度，使得利用早期返排歷史數(shù)據(jù)進(jìn)行水平井分段壓裂效果評價成為可能。迄今為止，國外學(xué)者在這方面做了許多理論研究，提出了很多返排早期（單相流生產(chǎn)階段）返排量和壓力的計算方法，以及基于壓裂液返排歷史數(shù)據(jù)的有效裂縫孔隙體積計算模型[1-8]。然而，這些模型均需要已知水力壓裂裂縫壓縮系數(shù)（以下簡稱裂縫壓縮系數(shù)）cf，而cf的計算結(jié)果存在較大的不確定性，直接影響了模型計算結(jié)果的準(zhǔn)確性。為此，筆者基于已有的計算模型，應(yīng)用產(chǎn)量不穩(wěn)定試井中的流動物質(zhì)平衡方法[9]，結(jié)合傳統(tǒng)的產(chǎn)量遞減方法[10]，提高了裂縫壓縮系數(shù)的計算精度，從而提高了有效裂縫孔隙體積的準(zhǔn)確性，為描述水力裂縫特征、評價壓裂后效果提供了參考。

1 壓裂液返排歷史的主要特征

水平井分段壓裂后返排作業(yè)時間一般都不長，雖然不同水平井的返排制度不同，但其返排歷史的主要特征基本類似。壓裂液返排期可以分為3個階段（如圖1所示）。

圖1 壓裂液返排期 3 個階段的主要特征示意Fig.1 Main characteristics of the three phases during the flowback of fracturing fluid

1）第一階段：單相壓裂液生產(chǎn)。在該階段，隨返排時間增長，套壓及井底流壓迅速降低，返排量達(dá)到峰值后緩慢下降。

2）第二階段：單相壓裂液生產(chǎn)。在該階段，套壓及井底流壓基本穩(wěn)定，后期隨著氣體開始產(chǎn)出并急劇增加，出現(xiàn)井筒卸載效應(yīng)（即儲層中的氣體進(jìn)入井筒后導(dǎo)致靜液壓力降低），套壓會有所升高。

3）第三階段：多相流（壓裂液、油/氣）生產(chǎn)。在該階段，返排量、井底流壓急劇下降，套壓持續(xù)下降，產(chǎn)氣量快速升高。

2 初始壓裂效果參數(shù)的定義

水平井在進(jìn)行水力壓裂時，井筒周圍會形成一個高壓力區(qū)，迫使壓裂液不斷向地層濾失。然而，由于致密儲層滲透性超低，濾失于基質(zhì)內(nèi)的壓裂液很少，絕大部分的壓裂液將進(jìn)入水力壓裂裂縫中，從而獲得很大的裂縫孔隙體積。當(dāng)壓裂作業(yè)結(jié)束、地面泵壓為零后，裂縫孔隙體積可能仍會增大。因此，評價壓裂效果，就是要準(zhǔn)確計算初始有效裂縫孔隙體積。由于地層存在天然裂縫，而水力壓裂裂縫只能溝通部分天然裂縫，未被溝通的天然裂縫，盡管在微地震資料中明顯可見，但對水平井產(chǎn)能沒有貢獻(xiàn)，因此不能算為有效裂縫孔隙體積。壓裂液返排過程中，部分裂縫閉合會導(dǎo)致裂縫孔隙體積減小，因此有效裂縫孔隙體積就是初始有效裂縫孔隙體積減去返排過程中損失的裂縫孔隙體積，即:

式中：Vef為有效裂縫孔隙體積，m3；Vfi為初始有效裂縫孔隙體積，計算方法見文獻(xiàn)[11]，m3；ΔVef為壓裂液返排過程中損失的裂縫孔隙體積，m3；cf為裂縫壓縮系數(shù)，MPa-1；Δpf為返排期間的壓力差，MPa；pfi為開始返排時裂縫內(nèi)的初始壓力（可用當(dāng)時的井底流壓來代替），MPa；pwf為井底流壓，MPa。

為了更好地表征壓裂液返排過程中裂縫閉合、水膨脹及地層流體（油/氣）膨脹的影響，引入了相應(yīng)代表返排驅(qū)替能量的指數(shù)因子[11]。

1）裂縫壓實驅(qū)替指數(shù)因子ICD。用該因子來反映壓裂液返排時有效裂縫閉合引起的壓裂液產(chǎn)出情況，其表達(dá)式為：

式中：ICD為裂縫壓實驅(qū)替指數(shù)因子；ct為綜合壓縮系數(shù)，MPa-1。

2）油/氣膨脹驅(qū)替指數(shù)因子IHD。用該因子來反映壓裂液返排時由油/氣膨脹引起的壓裂液產(chǎn)出情況，其表達(dá)式為：

式中：IHD為油/氣膨脹驅(qū)替指數(shù)因子；Sg為裂縫內(nèi)含氣飽和度；cg為裂縫內(nèi)氣體壓縮系數(shù)，MPa-1。

現(xiàn)場觀察發(fā)現(xiàn)，在壓裂液返排早期，氣井一般會同時產(chǎn)氣和產(chǎn)水，而油井則僅產(chǎn)水，因此Sg為0，則IHD也為 0。

3）水膨脹驅(qū)替指數(shù)因子IWD。用該因子來反映返排時水膨脹引起的壓裂液產(chǎn)出情況，其表達(dá)式為：

式中：IWD為水膨脹驅(qū)替指數(shù)因子；Sw為裂縫內(nèi)含水飽和度；cw為裂縫內(nèi)水壓縮系數(shù)。

在壓裂液返排期間，油氣井的ICD、IHD和IWD的和為1。圖2所示為北美某水平井的3個驅(qū)替指數(shù)因子隨返排時間的變化曲線[11]。由圖2可知，由于ICD的值最大，所以有效裂縫閉合是該井壓裂液返排的主要驅(qū)動機制。隨著時間的推移，ICD不斷減小,IHD不斷增大（即返排液中出現(xiàn)油氣并不斷增多），油/氣膨脹成為壓裂液返排的重要驅(qū)動機制。

圖2 北美某水平井的驅(qū)替指數(shù)因子隨返排時間的變化[11]Fig.2 Variation of displacement index with flowback time for a well in North America[11]

3 有效裂縫孔隙體積的計算方法

有效裂縫孔隙體積的準(zhǔn)確計算，關(guān)鍵在于提高裂縫壓縮系數(shù)cf的準(zhǔn)確性。為此，筆者結(jié)合產(chǎn)量不穩(wěn)定試井中流動物質(zhì)平衡方法和傳統(tǒng)的產(chǎn)量遞減方法，提高了初始有效裂縫體積和裂縫壓縮系數(shù)計算的準(zhǔn)確性，形成了計算裂縫有效孔隙體積的新方法。

3.1 初始有效裂縫孔隙體積的計算

由式（1）可知，初始有效裂縫孔隙體積Vfi及裂縫壓縮系數(shù)cf是計算有效裂縫孔隙體積的2個關(guān)鍵參數(shù)。考慮到部分壓裂液會進(jìn)入基質(zhì)中（體積非常?。蚨跏加行Я芽p孔隙體積Vfi一定小于壓裂液注入量VTinj。為此，作如下假設(shè)：1）返排壓裂液主要從有效裂縫流出，忽略從基質(zhì)流出的壓裂液；2）返排開始時，有效裂縫中僅充滿壓裂液。

B.Bai等人[12]發(fā)現(xiàn)，壓裂液返排量可以采用Arps產(chǎn)量調(diào)和遞減模型計算[3]，即：

式中：Di為壓裂液返排量的初始遞減率，1/年；qw為壓裂液產(chǎn)量（即返排速度，下同），m3/d；qwi為壓裂液初始產(chǎn)量，m3/d；Wp為壓裂液累計返排量，m3。

根據(jù)式（7），若以Wp為橫坐標(biāo)、qw的對數(shù)為縱坐標(biāo)做圖，可發(fā)現(xiàn)兩者之間基本呈線性關(guān)系，如圖3(a)所示（BT 代表油/氣開始產(chǎn)出的點，下同）。若壓裂液返排時井底壓力變化很大，需要將壓裂液產(chǎn)量進(jìn)行歸一化處理[13]，即縱坐標(biāo)變?yōu)閴毫w一化產(chǎn)量的對數(shù)，則兩者之間仍然基本呈線性關(guān)系（如圖3（b）所示），對應(yīng)的表達(dá)式為：

忽略從基質(zhì)流出的壓裂液，有以下近似關(guān)系[12-13]：

圖3 壓裂液累計返排量與產(chǎn)量的關(guān)系Fig.3 Variation of cumulative fracturing fluid flowback with production

式中：Bwi為壓裂液的體積系數(shù)，約為1.0；為最終累計返排量（即圖3中擬合直線與x軸交點的坐標(biāo)值），m3。

由式（9）可知，在缺少初始有效裂縫孔隙體積Vfi時，可以用最終累計返排量來代替。

3.2 壓裂液返排的主要流動特征

大量研究發(fā)現(xiàn)，頁巖及致密砂巖油氣藏中，水平井分段壓裂后壓裂液返排通常經(jīng)歷3個流動階段，如圖4所示[14-15]。

圖4 壓裂液返排主要流動特征示意及特征診斷圖Fig.4 Main flow characteristics and diagnostic plot of the fracturing fluid flowback

1）裂縫內(nèi)的早期線性流或徑向流階段（圖4中的FP1）。該階段中壓裂液沿著主裂縫流向井筒，主要流動特征為早期線性流或徑向流。若為徑向流，則可以利用該階段壓裂液返排歷史數(shù)據(jù)分析求得裂縫滲透率。但在實際應(yīng)用中，由于返排歷史數(shù)據(jù)的質(zhì)量問題，且早期受井筒存儲效應(yīng)的影響，很難在流態(tài)診斷圖上找到該流動階段，即很難利用該階段的返排歷史數(shù)據(jù)求解裂縫滲透率。

2）擬穩(wěn)態(tài)流階段（圖4中的FP2）。該階段為裂縫衰竭/閉合階段，即隨著壓裂液產(chǎn)出，裂縫中的壓力波不斷向外傳播，逐漸傳播到裂縫與基質(zhì)的相交處。此時，以裂縫與基質(zhì)相交處為邊界，邊界內(nèi)（裂縫內(nèi)）流體的流動原理與在推導(dǎo)物質(zhì)平衡方程時假設(shè)油氣藏為一個封閉油藏的原理類似。隨著壓裂液不斷流出，裂縫內(nèi)的壓力不斷下降，因此裂縫受壓縮閉合（類似于裂縫孔隙體積減?。Ｔ撾A段壓裂液的主要流動特征為產(chǎn)量不穩(wěn)定試井（RTA）中的擬穩(wěn)態(tài)流，可以通過流動物質(zhì)平衡方程來計算裂縫孔隙體積。

3）不穩(wěn)態(tài)線性流階段（圖4中的FP3）。隨著壓裂液不斷流向井筒，地層流體開始流入裂縫，與壓裂液形成多相流并同時流向井筒，該階段的主要流動特征一般為不穩(wěn)態(tài)線性流，可以利用該階段的返排歷史數(shù)據(jù)求得裂縫半長及基質(zhì)滲透率。

在以歸一化壓力導(dǎo)數(shù)d(Δp/q)/dlntc為縱坐標(biāo)、物質(zhì)平衡時間tc為橫坐標(biāo)的流動特征診斷圖上，3個壓裂液返排階段的流動特征分別是（如圖4（b）所示）：1）曲線斜率為0.5或0時，為不穩(wěn)態(tài)線性流或徑向流；2）曲線斜率為1.0時，為擬穩(wěn)態(tài)流；3）曲線斜率等于0.5時，為不穩(wěn)態(tài)線性流。

3.3 裂縫壓縮系數(shù)cf的計算

國內(nèi)外眾多學(xué)者開展了裂縫壓縮系數(shù)計算方法的研究[7,16]，其中R.Aguilera等人[16]通過室內(nèi)試驗得出裂縫壓縮系數(shù)與作用在裂縫上的凈壓力pn及裂縫內(nèi)“礦物質(zhì)”的含量有關(guān)。凈壓力是指巖石的最小主應(yīng)力σmin減去裂縫內(nèi)的流動壓力pf，即：

式中：pn為作用在裂縫上的凈壓力，MPa；σmin為巖石最小主應(yīng)力，MPa；pf為裂縫內(nèi)的流動壓力，MPa。

巖石的最小主應(yīng)力σmin可用裂縫關(guān)閉壓力pc（可由裂縫診斷注入測試（DFIT）數(shù)據(jù)分析獲?。﹣泶妫芽p內(nèi)的流動壓力pf可以用井底流壓pwf來代替，則式（10）變?yōu)椋?/p>

式中：pc為裂縫關(guān)閉壓力，MPa。

R.Aguilera等人[16]繪制的圖版如圖5所示，不同曲線代表裂縫內(nèi)“礦物質(zhì)”（即支撐劑）含量不同。J.D.Williams-Kovacs等人[17-18]認(rèn)為，“礦物質(zhì)”含量是指裂縫中支撐劑的體積分?jǐn)?shù)，并與裂縫孔隙度（φf）有關(guān)，或者說，“礦物質(zhì)”含量可以由裂縫孔隙度來表征。由于裂縫孔隙度的真實值較難確定，一般可以給出一個區(qū)間值。因此，“礦物質(zhì)”含量也是一個區(qū)間值，裂縫壓縮系數(shù)也相應(yīng)地是一個區(qū)間值。筆者對圖5中的曲線進(jìn)行了回歸分析，得到近似回歸方程為：

圖5 Aguilera等人的裂縫壓縮系數(shù)圖版[16]Fig. 5 Type curves for fracture compressibility coefficient by Aguilera et al.[16]

式中：a和b為回歸擬合參數(shù)。

回歸分析得到，當(dāng)“礦物質(zhì)”含量分別為0，10%，20%，30%，40%和50%時，a分別為0.023 8，0.026 1，0.026 4，0.027 4，0.045 6和0.053 0，b分別為0.640，0.663，0.683，0.710，0.809和0.893。

上述裂縫壓縮系數(shù)計算方法均存在未考慮支撐劑在裂縫中的不均勻分布、不是每口井都有DFIT數(shù)據(jù)、裂縫孔隙度難以確定等局限性。為此，H.Sabbir等人[13]提出以下假設(shè)：1）裂縫最初由壓裂液飽和；2）初期裂縫中不存在流動的油氣，忽略由基質(zhì)流入裂縫的油氣；3）在擬穩(wěn)態(tài)流階段，壓力隨時間的變化是常量，通過裂縫邊界控制流（BDF）階段范圍內(nèi)的數(shù)據(jù)在時間和空間上取平均值來描述裂縫可壓縮性?；谠摷僭O(shè)，他們提出了改進(jìn)的裂縫壓縮系數(shù)計算方法：

1）整理壓裂液返排期的流量及井底流壓數(shù)據(jù)；

2）利用體積系數(shù)將地面流量轉(zhuǎn)換成地下流量；

3）計算產(chǎn)量歸一化壓力RNP及其導(dǎo)數(shù)dRNP，以物質(zhì)平衡時間tc為橫坐標(biāo)，產(chǎn)量歸一化壓力RNP及其導(dǎo)數(shù)dRNP為縱坐標(biāo)，在雙對數(shù)坐標(biāo)系上繪制了如圖4（b）所示的壓裂液流動特征診斷圖；

4）在壓裂液流動特征診斷圖上找出斜率為1的流動段（邊界控制的擬穩(wěn)態(tài)流）的數(shù)據(jù)，并繪制在線性坐標(biāo)系上，橫坐標(biāo)為物質(zhì)平衡時間tc，縱坐標(biāo)為產(chǎn)量歸一化壓力RNP，擬合出一條直線，直線段斜率記為mBDF；

5）在半對數(shù)坐標(biāo)系中，作壓力歸一化產(chǎn)量PNR與累計返排量Wp的關(guān)系曲線，擬合直線外推到廢棄壓裂液返排速度0.16 m3/d、廢棄壓力0.69 MPa處，廢棄點的橫坐標(biāo)值即為最終累計返排量；

6）根據(jù)式（9），將初始有效裂縫孔隙體積Vfi近似為最終累計返排量，再利用下式求出cf：

式中：Δpf為返排期擬穩(wěn)態(tài)流階段的壓力差，MPa；pwfe為擬穩(wěn)態(tài)流結(jié)束處的井底流壓，即在步驟4）中找出的最后一個數(shù)據(jù)點對應(yīng)的井底流壓，MPa。

分析認(rèn)為，該計算方法存在3個缺陷：1）因為其理論依據(jù)是假設(shè)平均裂縫孔隙體積等于初始裂縫孔隙體積與最終裂縫孔隙體積的算術(shù)平均值，但實際上裂縫孔隙體積隨著壓力、裂縫壓縮系數(shù)、平均含水飽和度等參數(shù)的變化而不斷變化，且不是一個簡單的線性變化過程；2）在用調(diào)和遞減預(yù)測模型計算最終累計返排量時未考慮基質(zhì)向裂縫的滲流，影響了初始有效裂縫孔隙體積計算結(jié)果的準(zhǔn)確性；3）未考慮支撐劑的壓實效應(yīng)對裂縫壓縮系數(shù)的影響。因此，筆者對該方法進(jìn)行了修正，采用C.R.Clarkson等人[15]提出的流動物質(zhì)平衡方法計算初始有效裂縫孔隙體積，以壓力歸一化產(chǎn)量PNR為縱坐標(biāo)、歸一化累計返排量為橫坐標(biāo)作流動物質(zhì)平衡分析圖（如圖6所示），擬合直線與x軸交點的橫坐標(biāo)即為初始有效裂縫孔隙體積。

圖6 壓裂液返排歷史的流動物質(zhì)平衡分析示意Fig.6 Flowing material balance analysis of fracturing fluid flowback history

其中，歸一化累計返排量為：

式中：W'p為歸一化累計返排量，m3/MPa；φfi為充滿支撐劑裂縫的初始孔隙度；φf為充滿支撐劑裂縫的孔隙度；為裂縫內(nèi)平均壓力，MPa。

在利用式（17）計算綜合壓縮系數(shù)ct時，需要已知裂縫壓縮系數(shù)cf。因此，計算過程將是一個迭代過程，可以將由調(diào)和遞減模型計算得到的最終累計返排量近似為初始有效裂縫孔隙體積，并作為迭代計算的初始值。由此，裂縫壓縮系數(shù)的計算流程需要增加第7步，即用第6步得到的裂縫壓縮系數(shù)cf值，再由圖6的流動物質(zhì)平衡曲線得到初始有效裂縫孔隙體積Vfi值，再由式（13）計算cf值。由此迭代，最終分別得到收斂的cf值和Vfi值。

4 計算實例

以國內(nèi)某致密油水平井A為例，該井采用分段壓裂完井，共注入壓裂液約17 600 m3，壓裂液返排約1個月后開始有油產(chǎn)出，壓裂液返排歷史如圖7所示。

圖7 水平井A的壓裂液返排歷史示意Fig.7 Fracturing fluid flowback history of horizontal well A

根據(jù)水平井A的壓裂液返排數(shù)據(jù)，繪制壓裂液返排量遞減分析圖，如圖8所示。從圖8可以看出，由壓裂液返排量遞減分析或壓力歸一化產(chǎn)量遞減分析得到的最終累計返排量基本一致，約為6 855.8 m3。用產(chǎn)量歸一化壓力及其導(dǎo)數(shù)、真實時間及物質(zhì)平衡時間作壓裂液流動特征診斷圖，見圖9。從圖9可以看出，流動曲線存在明顯的斜率為1的流動段，表明返排期單相壓裂液生產(chǎn)階段，流體流動出現(xiàn)了非常明顯的擬穩(wěn)態(tài)流。

選取斜率為1的流動段的數(shù)據(jù)作產(chǎn)量歸一化壓力導(dǎo)數(shù)與物質(zhì)平衡時間的關(guān)系曲線（如圖10所示），求出直線段斜率mBDF為 0.251 8。按上述裂縫壓縮系數(shù)和初始有效裂縫孔隙體積計算方法，將最終累計返排量的預(yù)測值作為初始有效裂縫孔隙體積，由式（13）計算裂縫壓縮系數(shù)，作流動物質(zhì)平衡分析圖（如圖11所示），將數(shù)據(jù)回歸直線與x軸交點的橫坐標(biāo)作為新的初始有效裂縫孔隙體積，代入式（13）再次計算裂縫壓縮系數(shù)，不斷迭代計算，最終得到初始有效裂縫孔隙體積Vfi為6 625.6 m3，裂縫壓縮系數(shù)cf為0.021 75 MPa-1。

圖8 水平井A的壓裂液返排量遞減分析結(jié)果Fig. 8 Decline analysis of fracturing fluid flowback of horizontal well A

圖9 水平井A的壓裂液返排流動特征分析診斷圖Fig. 9 Flow characteristic diagnostic plot of fracturing fluid flowback of horizontal well A

圖10 水平井A的壓裂液返排分析特種曲線Fig. 10 Special curve of the fracturing fluid flowback for horizontal well A

圖11 水平井A的壓裂液返排期流動物質(zhì)平衡分析圖Fig. 11 Flowing material balance analysis of the fracturing fluid flowback for horizontal well A

水平井A從開始返排到擬穩(wěn)態(tài)結(jié)束時的壓力差Δpf約為20.13 MPa，因此裂縫體積損失ΔVef=VficfΔpf=20.13×0.021 75×6 625.6=2 902.0 m3，則投產(chǎn)后的有效裂縫孔隙體積Vef=Vfi-ΔVef= 6 625.6 - 2 902.0 =3 723.6 m3，約為壓裂液總注入量的21%。也就是說，只 有約21%的壓裂液對形成裂縫導(dǎo)流能力有貢獻(xiàn)。

5 結(jié) 論

1）壓裂液返排的地下流體流動機理是裂縫的閉合和流體的膨脹。

2）受裂縫壓縮系數(shù)計算準(zhǔn)確性的影響，基于壓裂液返排歷史數(shù)據(jù)的有效裂縫孔隙體積常規(guī)計算方法存在計算結(jié)果不夠準(zhǔn)確的問題，筆者結(jié)合產(chǎn)量不穩(wěn)定試井中的流動物質(zhì)平衡方法和傳統(tǒng)的產(chǎn)量遞減方法，以提高初始有效裂縫體積和裂縫壓縮系數(shù)的計算準(zhǔn)確性為核心，提出了計算有效裂縫孔隙體積的新方法。

3）計算實例表明，返排期單相壓裂液生產(chǎn)階段，流體流動出現(xiàn)了非常明顯的擬穩(wěn)態(tài)流，只有部分注入的壓裂液對形成裂縫導(dǎo)流能力有貢獻(xiàn)。

4）致密油氣藏水平井壓裂作業(yè)中，地層或多或少地會出現(xiàn)漏失，因此，筆者提出的計算有效裂縫孔隙體積的新方法在對分析致密油氣藏水平井壓裂液返排歷史數(shù)據(jù)具有比較強的實用性。