馮月麗,劉月田,丁祖鵬,陳健,毛小龍,周月波

(1.中國石油大學(xué)(北京)油氣資源與探測國家重點實驗室,北京102249;2.中海油研究總院,北京 100028;3.中國石化中原油田分公司石油工程研究院,河南 濮陽 457001)

在研究裂縫性滲流介質(zhì)的變形特征方面,國內(nèi)外很多專家學(xué)者已經(jīng)做了大量的工作[1-12],其中,大部分?jǐn)?shù)學(xué)模型中,考慮裂縫性介質(zhì)變形時,使用的都是通過實驗數(shù)據(jù)擬合得到的經(jīng)驗公式[1-5],或者引入滲透率模數(shù)[6-8].這些處理方法存在的主要問題,是介質(zhì)的變形系數(shù)不易確定,因此,實際應(yīng)用中存在一定的局限性.本文以容易測量的基本彈性參數(shù)為基礎(chǔ),建立了一組平行裂縫的變形計算模型,并在此基礎(chǔ)上給出了裂縫滲透率和介質(zhì)整體滲透率張量的計算方法.該模型可以解決變形系數(shù)不易確定的問題,同時能夠綜合考慮多因素對裂縫變形的影響,對預(yù)測復(fù)雜裂縫網(wǎng)絡(luò)的變形、計算裂縫性滲流介質(zhì)的滲透率動態(tài)變化具有重要參考價值.

1 動態(tài)滲透率模型

1.1 特征單元體裂縫總條數(shù)計算

取滲流單元內(nèi)包含一組平行裂縫的典型微元(特征單元體)作為研究對象,建立參照直角坐標(biāo)系xoy,裂縫走向與x軸的夾角為β,平面上裂縫(藍(lán)線標(biāo)識)等間距d分布,縱向上裂縫垂直切穿整個單元體.特征單元體的寬、長分別為L1,L2,高為h,單元體四周均為固定邊界(見圖1).圖中紅色虛線將裂縫分為2部分.其中:虛線左側(cè)沿AB方向的裂縫間距為d1;虛線右側(cè)沿CB方向的裂縫間距為d2.

圖1 一組平行裂縫理論模型示意

AB邊上的裂縫條數(shù)n1(取正整數(shù))為

CB邊上的裂縫條數(shù)n2(取正整數(shù))為

假設(shè)L1=L2=L,則特征單元體裂縫總條數(shù)n為

1.2 裂縫動態(tài)滲透率計算模型

為建立裂縫變形計算模型,提出2個假設(shè):

1)基質(zhì)部分為各向同性線彈性體,同時假設(shè)基質(zhì)系統(tǒng)彈性模量為Em,總應(yīng)力變化量為Δσt,平均流體壓力變化量為Δpm,應(yīng)變?yōu)棣舖,有效應(yīng)力系數(shù)為αm,則由有效應(yīng)力原理[12]可得:

2)裂縫的應(yīng)力-位移滿足線性關(guān)系.設(shè)裂縫在初始壓力p0時開度為b0,裂縫系統(tǒng)彈性模量為Ef,裂縫面所受法向應(yīng)力變化量為Δσnf,法向位移為Δunf,則有:

設(shè)Δuym為基質(zhì)沿y方向位移,Δuyf為裂縫沿y方向的總位移,則整個特征單元體在y方向的總位移Δuy=Δuym+Δuyf.

單元體四周邊界均為固定邊界,因此沿x,y方向的總位移均為0,即:

式中:Δux,Δuxm,Δuxf分別為單元體、基質(zhì)、裂縫沿 x 方向的總位移.

裂縫系統(tǒng)沿y方向的總位移為

基質(zhì)系統(tǒng)沿y方向的總位移為

式中:L0為模型初始寬度,mm;Δσyt為模型沿y方向的總應(yīng)力變化量,MPa;Δpym為基質(zhì)系統(tǒng)沿y方向上承受的流體壓力變化量,MPa.

聯(lián)立式(3)、式(5)、式(7),可得n條裂縫沿y方向的總位移:

式中:by0為裂縫沿y方向的初始開度,mm;αf為裂縫系統(tǒng)有效應(yīng)力系數(shù);Δpyf為裂縫系統(tǒng)沿y方向上承受的流體壓力變化量,MPa;d0為裂縫在壓力p0時的初始裂縫間距,mm.

由式(6)得:

再將式(8)、式(9)代入式(10),得:

假設(shè)基質(zhì)系統(tǒng)和裂縫系統(tǒng)中的流體壓力瞬間達到平衡,則有 Δpyf=Δpym=Δp,從而得到:

將式(12)代入式(9),得:

則裂縫法向位移Δuynf為

因此,裂縫變形后開度b為

裂縫內(nèi)蘊滲透率Kff[13]為

壓力為p0時的裂縫滲透率K0和壓力為p時的裂縫滲透率K與裂縫開度的關(guān)系[14]為

由式(17)可得平行裂縫方向的裂縫滲透率(沒有特殊說明的裂縫滲透率均為平行裂縫方向的值):

再沿裂縫發(fā)育方向建立裂縫直角坐標(biāo)系 x′o′y′,假設(shè)基質(zhì)滲透率為Km,則單元體的滲透率張量K′表示為

將K′轉(zhuǎn)換為xoy參照坐標(biāo)系下的滲透率張量K,其表達式[15]為

2 模型驗證

為驗證裂縫變形計算模型的正確性,開展了滲流介質(zhì)中包含一組裂縫的二維滲流物理模擬實驗研究.常規(guī)巖心實驗造縫困難,變形量小,不易于觀察實驗結(jié)果,因此,本次研究的實驗材料選用造縫可控,裂縫變形較容易的相似多孔介質(zhì)材料.該裂縫性滲流介質(zhì)的基本參數(shù)如表1所示.

圖2為實驗?zāi)Ｐ褪疽鈭D.由圖2可知,壓力梯度方向平行于x軸,邊界AB和CD分別為入口端、出口端,BC和AD均為封閉邊界.模型邊長為L,基質(zhì)部分為各向同性應(yīng)變介質(zhì),實驗流體為單相蒸餾水.實驗通過改變模型中的平均孔隙流體壓力,研究不同裂縫分布情況下,裂縫的開度變化及其對滲透率的影響.

表1 裂縫性滲流介質(zhì)主要物理參數(shù)

圖2 實驗?zāi)Ｐ褪疽?/p>

以模型中包含一組與x軸方向夾角(裂縫傾角)分別為0°,45°裂縫的實驗測試結(jié)果為例進行模型驗證(見圖3).由圖可以看出,實驗測試和理論計算的數(shù)據(jù)吻合度高,說明本文建立的裂縫滲透率計算模型是合理的,可以用來計算封閉邊界條件下的裂縫變形量及裂縫滲透率隨孔隙流體壓力的變化規(guī)律.

圖3 實驗測試折算與理論計算值對比

3 變形特性

3.1 b,Kf與特征單元體幾何、彈性參數(shù)的關(guān)系

3.1.1 裂縫初始開度和孔隙流體壓力

圖4為裂縫開度、滲透率隨流體壓力、裂縫初始開度變化的情況.由圖可知:初始裂縫開度相同時,裂縫開度和裂縫滲透率均隨孔隙流體壓力增大而增大;孔隙流體壓力保持不變時,裂縫初始開度越大的裂縫,變形后開度和相應(yīng)的裂縫滲透率也越大.

3.1.2 裂縫初始間距

圖5為裂縫開度、滲透率與裂縫間距的關(guān)系.由圖可知,裂縫間距越大,相同大小單元體內(nèi)包含的裂縫條數(shù)越少,因此,平均孔隙流體壓力保持不變時,裂縫間距越大,總開度越小,相應(yīng)的裂縫滲透率也越小.

圖4 裂縫開度、裂縫滲透率隨孔隙流體壓力、裂縫初始開度的變化情況

圖5 裂縫開度、裂縫滲透率與裂縫間距的關(guān)系

3.1.3 彈性模量

基質(zhì)系統(tǒng)和裂縫系統(tǒng)的彈性模量相互關(guān)聯(lián),其中一個不會脫離另一個而單獨保持不變.因此,對同一滲流介質(zhì)而言,需要討論二者同時發(fā)生變化時的裂縫滲透率變化規(guī)律.如圖6所示,基質(zhì)系統(tǒng)彈性模量越大,有效應(yīng)力改變量相同時,基質(zhì)系統(tǒng)的應(yīng)變越小.因為裂縫系統(tǒng)與基質(zhì)系統(tǒng)的位移大小相等,方向相反,因而,裂縫系統(tǒng)的變形量也越小,裂縫開度和裂縫滲透率隨著彈性模量的增大而減小.裂縫系統(tǒng)彈性模量增加時,變化規(guī)律相同.

3.1.4 有效應(yīng)力系數(shù)

基質(zhì)與裂縫系統(tǒng)的有效應(yīng)力系數(shù)之間也存在一定的聯(lián)系[11-12],一個發(fā)生改變時另一個也發(fā)生變化.因此,對于同一個特征單元體而言,需要討論它們同時發(fā)生變化時,裂縫滲透率的變化規(guī)律.圖7為裂縫開度、裂縫滲透率與有效應(yīng)力系數(shù)的關(guān)系曲線.由圖可以看出,有效應(yīng)力系數(shù)越大,裂縫開度和滲透率越小.這是因為,有效應(yīng)力系數(shù)越大,當(dāng)總應(yīng)力和孔隙流體變化量不變時,有效應(yīng)力越小,基質(zhì)和裂縫系統(tǒng)的變形量越小.

3.1.5 裂縫傾角

圖8為裂縫開度、裂縫滲透率與裂縫傾角的關(guān)系曲線.由圖可知,隨著裂縫傾角增大,裂縫開度和滲透率均呈現(xiàn)先增后減的趨勢,裂縫傾角對裂縫滲透率的影響較復(fù)雜.當(dāng)β=45°時,裂縫滲透率達到最大值.

圖6 裂縫開度、裂縫滲透率與彈性模量的關(guān)系

圖7 裂縫開度、裂縫滲透率與有效應(yīng)力系數(shù)的關(guān)系

圖8 裂縫開度、裂縫滲透率與裂縫傾角的關(guān)系

3.2 特征單元體滲透率張量變化

3.2.1 裂縫性滲流介質(zhì)滲透率主值方向變化

不發(fā)育裂縫的各向同性滲流介質(zhì)中,滲透率處處相等.如圖9所示,當(dāng)介質(zhì)內(nèi)發(fā)育一組裂縫時,介質(zhì)的總體滲透率將呈現(xiàn)各向異性特征,且滲透率最大主值方向總是沿著裂縫發(fā)育方向,而最小主值方向總是垂直于裂縫發(fā)育方向.

3.2.2 裂縫變形對滲流介質(zhì)滲透率張量的影響

圖10為發(fā)育一組傾角為0°,45°裂縫時滲透率張量隨孔隙流體壓力的變化圖.由圖可知,隨流體壓力增大,介質(zhì)滲透率最大主值(橢圓長半軸的長度)增大,但主值方向不變,反之亦然.說明基質(zhì)為各向同性時,裂縫變形只對介質(zhì)的滲透率主值大小有影響.

圖9 各向異性滲透率張量隨裂縫傾角的變化情況

圖10 各向異性滲透率張量隨流體壓力的變化

4 結(jié)論

1)孔隙流體壓力增大,裂縫變形量和裂縫滲透率均增大;裂縫間距越大,裂縫滲透率越小;彈性模量越大,裂縫變形量越小;有效應(yīng)力系數(shù)越大,裂縫變形量越小;隨著裂縫傾角增大,裂縫滲透率呈先增大后減小的趨勢.

2)發(fā)育一組平行裂縫的滲流介質(zhì),其滲透率最大主值方向始終沿著裂縫發(fā)育方向,裂縫發(fā)育方向改變時,介質(zhì)滲透率最大主值方向也會隨之改變.基質(zhì)為各向同性時,裂縫隨孔隙流體壓力發(fā)生的變形只對介質(zhì)的滲透率主值大小有影響.

3)相似多孔介質(zhì)微觀物理模擬實驗測試值和理論計算結(jié)果吻合度較高,驗證了一組平行裂縫變形計算模型的準(zhǔn)確性和實用性.模型能較全面考慮裂縫系統(tǒng)和基質(zhì)系統(tǒng)的彈性參數(shù),相較于目前常用的計算裂縫介質(zhì)變形的經(jīng)驗公式,基礎(chǔ)參數(shù)更容易確定,可進一步應(yīng)用于裂縫網(wǎng)絡(luò)的變形和滲透率計算.