陳飛宇,張廣智,周 游,劉俊州,韓 磊

(1.中國(guó)石油大學(xué)(華東)深層油氣重點(diǎn)實(shí)驗(yàn)室,山東青島266580;2.中國(guó)石油大學(xué)(華東),山東青島266580;3.中國(guó)石油化工股份有限公司石油勘探開發(fā)研究院,北京100083)

彈性模量的測(cè)量可分為動(dòng)力法和靜力法[1]。在靜力法中,靜態(tài)楊氏模量通常利用排水條件下的單軸和三軸試驗(yàn)得到,表征不受孔隙流體影響的情況下巖石的力學(xué)性質(zhì);在動(dòng)力法中,可以在室內(nèi)或野外測(cè)量縱、橫波速度以及密度,據(jù)此計(jì)算出動(dòng)態(tài)模量。儲(chǔ)層巖石的動(dòng)靜態(tài)模量(楊氏模量或泊松比)之間的關(guān)系對(duì)油氣儲(chǔ)集層的表征和鉆井工程有較為重要的意義[2]。

一般來說,靜態(tài)模量測(cè)量需從所研究的地層中取出巖心樣品,這種方法對(duì)巖石具有破壞性、耗時(shí)較多、價(jià)格較貴、需要巖石樣本的數(shù)量大[3],由于實(shí)際條件限制,常常不能獲得井中每段地層的巖石樣本,因此得到的沿井筒方向信息不連續(xù)[4]。KING[5]和NAJIBI等[3]通過大量的巖石物理實(shí)驗(yàn),提出了通過擬合動(dòng)靜態(tài)模量之間線性或非線性的經(jīng)驗(yàn)公式計(jì)算無巖心段巖石的靜態(tài)模量;ASEF等[6]考慮了孔隙度對(duì)動(dòng)靜態(tài)模量之間的影響,建立了一種不同孔隙度下動(dòng)靜態(tài)模量轉(zhuǎn)換的半經(jīng)驗(yàn)公式。王飛等[7]研究了溫度和壓力對(duì)頁巖動(dòng)靜態(tài)模量的影響,并建立了一定溫度、壓力條件下的動(dòng)靜態(tài)模量的轉(zhuǎn)換模型。邊會(huì)媛等[8]在實(shí)驗(yàn)室條件下研究了致密砂巖動(dòng)靜態(tài)模量之間的轉(zhuǎn)換關(guān)系與溫度和壓力之間的關(guān)系。盛英帥等[9]通過模擬致密砂巖儲(chǔ)層環(huán)境,研究了不同孔隙度下動(dòng)靜態(tài)模量的經(jīng)驗(yàn)公式。葛洪魁等[10]認(rèn)為動(dòng)靜態(tài)模量轉(zhuǎn)換的機(jī)制是由多種因素控制,如:頻散效應(yīng)、排水條件、應(yīng)變振幅等,但并沒有給出這些影響因素的數(shù)學(xué)表達(dá)。從研究現(xiàn)狀可以看出,靜態(tài)模量的預(yù)測(cè)大多是通過實(shí)驗(yàn)測(cè)量或經(jīng)驗(yàn)公式擬合的方法進(jìn)行的。這些經(jīng)驗(yàn)公式雖然可以用于預(yù)測(cè)靜態(tài)模量,但缺乏明確的物理含義,難以解釋動(dòng)靜態(tài)模量之間的復(fù)雜關(guān)系,應(yīng)用范圍有限。

在測(cè)井中,由于慢地層的限制和井筒不穩(wěn)定的影響,很難獲得橫波速度。通過巖石物理模型可以獲得儲(chǔ)層巖石的縱、橫波速度,計(jì)算出動(dòng)態(tài)模量[11]。XU等[12]最早提出了通過V-R-H模型、K-T模型和微分等效模型構(gòu)建砂巖巖石物理模型的方法。劉倩等[13]針對(duì)致密砂巖低孔、低滲的特點(diǎn),提出了含有孤立孔隙的巖石物理建模方法。張益明等[14]針對(duì)致密砂巖復(fù)雜的孔隙結(jié)構(gòu)提出了雙重孔隙模型來計(jì)算巖石骨架的彈性模量的方法。印興耀等[15-16]指出,在巖石物理建模中裂縫形狀的表征與復(fù)雜裂縫的描述是影響模型精確性的因素之一。

本文通過構(gòu)建孔隙裂縫型致密砂巖巖石物理模型獲得動(dòng)態(tài)模量,通過分析頻散效應(yīng)、排水效應(yīng)、應(yīng)變振幅對(duì)動(dòng)靜態(tài)模量轉(zhuǎn)換的物理影響,建立一套動(dòng)靜態(tài)模量的轉(zhuǎn)換流程。最后利用川西地區(qū)須家河組致密砂巖驗(yàn)證了該方法的可靠性,為該地區(qū)地下致密砂巖靜態(tài)模量的預(yù)測(cè)提供參考。

1 技術(shù)方法

1.1 構(gòu)建致密砂巖巖石物理模型及動(dòng)態(tài)模量計(jì)算

致密砂巖儲(chǔ)層具有礦物成分多樣、孔隙結(jié)構(gòu)復(fù)雜、孔隙中流體分布不均勻、微裂縫發(fā)育和強(qiáng)烈的非均質(zhì)性等特點(diǎn)[11]。致密砂巖主要礦物為石英、長(zhǎng)石、泥質(zhì)等。不同礦物之間主要以線性和凹凸?fàn)罱佑|,使得巖石的孔隙形狀具有明顯的不規(guī)則狀。另外,壓實(shí)作用、膠結(jié)作用等也是導(dǎo)致孔隙結(jié)構(gòu)復(fù)雜的原因[17]。致密砂巖的非均質(zhì)性主要體現(xiàn)在孔隙結(jié)構(gòu)的復(fù)雜和孔隙流體分布狀態(tài)的多樣,大多以斑塊飽和為主[18]。

由于致密砂巖埋深大,巖石較為致密,常規(guī)巖石物理模型采用V-R-H模型計(jì)算的巖石基質(zhì)只是將巖石中的礦物含量混合起來,上下限的范圍較大,模量誤差較大。因此采用更為準(zhǔn)確的加權(quán)平均邊界模型來計(jì)算巖石基質(zhì)的模量。我們?cè)赬u-White建模[12]的思路上,考慮致密砂巖孔隙、裂縫的特征,構(gòu)建孔隙裂縫型巖石物理等效模型,流程如下:

第①步:根據(jù)致密砂巖礦物成分及其含量,考慮各礦物組分的接觸關(guān)系,用H-S邊界加權(quán)平均[19]計(jì)算由基本礦物(石英、長(zhǎng)石、泥質(zhì))混合組成各向同性的巖石基質(zhì)模量。

KHS+=Λ(μmax)

(1)

KHS-=Λ(μmin)

(2)

μHS+=Γ[ζ(Kmax,μmax)]

(3)

μHS-=Γ[ζ(Kmin,μmin)]K

(4)

Km=a+KHS++a-KHS-

(5)

μm=a+μHS++a-μHS-

(6)

(7)

式中:Λ(z)=〈1/[K(r)+4/(3z)]〉-1-4/(3z);Γ(z)=〈1/[μ(r)+z]〉-1-z;ζ(K,μ)=μ(9K+8μ)/6(K+2μ);〈·〉表示對(duì)每種礦物按其體積含量進(jìn)行加權(quán)平均;a+,a-分別為HS+上限和HS-下限的權(quán)值,且a++a-=1;fi是第i組分的體積分?jǐn)?shù);Km,μm,ρm分別為巖石基質(zhì)的體積模量、剪切模量和密度。

第②步:致密砂巖孔隙較為孤立,連通性較差,孔隙之間的影響小,采用K-T模型[20]在巖石基質(zhì)中加入孤立孔隙形成干燥的巖石骨架。此時(shí)巖石的體積模量K0和剪切模量μ0為:

(8)

(9)

(10)

式中:Ki,μi分別為孔隙包含物的體積模量和剪切模量。由于加入的是孔隙空腔,則Ki,μi等效為0,此時(shí)有:

∑xi(Ki-Km)Pmi=-φPKm

(11)

∑xi(μi-μm)Qmi=-φQμm

(12)

式中:φ為孔隙度;Pmi,Qmi與為孔隙形狀有關(guān)的參數(shù),表示孔隙結(jié)構(gòu)對(duì)巖石彈性模量的影響。

(13)

(14)

(15)

式中:λ和μ為各向同性背景介質(zhì)的拉梅常數(shù);e為裂縫密度;U11和U33為與裂縫有關(guān)的彈性參數(shù):

(16)

(17)

式中:l,k為與裂縫特征有關(guān)的參數(shù)。

(18)

(19)

式中:k′和μ′分別為裂縫填充物的體積模量和剪切模量。

第④步:計(jì)算氣、水混合流體的模量。由于致密砂巖的流體大多呈現(xiàn)為斑塊飽和狀態(tài),利用Voigt平均計(jì)算孔隙中混合流體的體積模量(Kf)與密度(ρf):

Kf=SwKw+(1-Sw)Kg

(20)

ρf=Swρw+(1-Sw)ρg

(21)

式中:Kw,Kg分別為水和氣的體積模量;ρw,ρg分別為水和氣的密度;Sw為含水飽和度。

(22)

第⑥步:巖石縱橫波速度及巖石彈性模量的求取。根據(jù)剛度系數(shù)矩陣計(jì)算飽和巖石縱、橫波速度和動(dòng)態(tài)模量:

(23)

(24)

(25)

致密砂巖巖石物理建模具體流程如圖1所示。

圖1 致密砂巖巖石物理建模流程

1.2 動(dòng)態(tài)模量到靜態(tài)模量的轉(zhuǎn)換

巖石動(dòng)靜態(tài)模量的差異是由多種因素共同造成的,主要包括:頻散效應(yīng)、排水效應(yīng)、應(yīng)變振幅等,下面分別分析這些因素的影響,建立動(dòng)態(tài)到靜態(tài)模量的轉(zhuǎn)換流程。

1.2.1 頻散效應(yīng)的影響

對(duì)于動(dòng)態(tài)模量而言,一般實(shí)驗(yàn)室中測(cè)量方法使用的超聲波的頻率為105～107Hz,地震方法的頻率是1～100Hz,聲波測(cè)井方法的頻率為103～104Hz;靜態(tài)模量可看成巖石力學(xué)測(cè)試過程中應(yīng)變率在頻率為1Hz時(shí)的彈性波的應(yīng)變率[1],而頻散效應(yīng)使得高頻模量高于低頻模量,因此在動(dòng)靜態(tài)模量的轉(zhuǎn)換中首先考慮頻散效應(yīng),即將高頻的動(dòng)態(tài)模量轉(zhuǎn)變到低頻模量。

波在巖石中傳播表現(xiàn)出的速度頻散和衰減,主要由孔隙、流體、滲透率和巖石的基質(zhì)所造成[23],在飽和巖石中速度的頻散與衰減主要由波在微裂縫和粒間孔隙之間的局部流體流動(dòng)導(dǎo)致的[24-26]。致密砂巖經(jīng)過強(qiáng)烈的壓實(shí)作用,膠結(jié)作用和成巖作用,普遍表現(xiàn)為低孔、低滲、孔隙連通性較差的特性。在1～1000Hz頻帶范圍,致密砂巖的縱橫波速度的頻散與衰減較為顯著,隨著頻率的不斷升高,縱橫波速度的增長(zhǎng)逐漸降低并趨于一個(gè)定值[27]。采用一種有效介質(zhì)的模型(孔隙-裂縫的頻散模型)描述高頻極限與低頻極限下的彈性模量[28]。

飽和巖石高頻極限下的彈性模量為:

(26)

(27)

(28)

(29)

式中:Em為巖石基質(zhì)的楊氏模量;α為裂縫的高寬比。

干燥巖石高頻極限下的彈性模量計(jì)算:令(26)式和(27)式中δp和δc→∞,進(jìn)而δp/(1+δp)和δc/(1+δc)→1,干燥巖石與飽和巖石的剪切模量具有相同的表達(dá)式。

飽和巖石低頻極限下的體積模量為:

(30)

利用參數(shù)Disp表征高頻和低頻模量之間頻散與衰減的大小:

(31)

式中:M為彈性模量,可以是體積模量、剪切模量和楊氏模量。下面以楊氏模量為例探究孔隙、裂縫以及裂縫的高寬比α對(duì)頻散大小的影響,當(dāng)致密砂巖的體積模量為36GPa、剪切模量為31GPa時(shí)的頻散與孔隙、裂縫之間的關(guān)系如圖2所示。隨著孔隙、裂縫的增加,頻散也逐漸增大,其中孔隙對(duì)頻散影響相對(duì)較小,裂縫隨著裂縫高寬比的增加頻散逐漸降低。

1.2.2 排水效應(yīng)的影響

對(duì)于飽和巖石,靜態(tài)模量的測(cè)量過程中,在三軸壓縮機(jī)的作用下多孔巖石受到擠壓,巖石發(fā)生破碎,裂縫趨于閉合,孔隙中的流體排出與外界發(fā)生交換,進(jìn)而保持孔隙壓力處于穩(wěn)定狀態(tài)?？紫秹毫υ谝欢ǔ潭壬先趸藝鷫旱挠绊?使巖石的彈性模量降低。在動(dòng)態(tài)模量的測(cè)量中,不會(huì)對(duì)巖石造成破壞,一般認(rèn)為彈性波無法引起巖石孔隙流體的交換,是不排水的。巖石孔隙中的流體在高頻聲波的作用下,在較短的時(shí)間內(nèi)發(fā)生變形作用且不與外界發(fā)生流體交換,導(dǎo)致巖石中孔隙壓力增大,從而會(huì)抵抗巖石的壓縮,使巖石更加堅(jiān)硬。對(duì)于干燥巖石排水影響相對(duì)較小。各向同性巖石排水的影響可以用Biot-Gassmann來描述[29]:

(32)

(33)

式中:Kdra,μdra為巖石排水的體積模量、剪切模量各向異性巖石排水的影響可以用Brown和Korringa公式(公式(22))描述。

1.2.3 應(yīng)變振幅的影響

動(dòng)態(tài)模量的應(yīng)變振幅是由超聲波的機(jī)械振動(dòng)引起的,范圍為10-7～10-6,靜態(tài)模量的應(yīng)變振幅是巖石受到應(yīng)力擠壓造成的,范圍為10-5～10-3[30],兩者存在數(shù)量級(jí)的差異,這是動(dòng)靜態(tài)模量差異的主要原因。動(dòng)態(tài)測(cè)量中產(chǎn)生的應(yīng)變振幅很難使裂縫閉合。靜態(tài)測(cè)量中巖石應(yīng)力應(yīng)變的第一階段是非彈性變形,在這一階段裂縫相較于孔隙更容易閉合,孔隙閉合所需的壓力遠(yuǎn)遠(yuǎn)大于常規(guī)三軸壓縮試驗(yàn)的最大壓力,所以只需要考慮裂縫對(duì)巖石變形的影響。對(duì)于均勻完全彈性的材料例如鋼、鐵等,動(dòng)靜態(tài)模量完全相等。裂縫對(duì)巖石非彈性形變的影響,采用WALSH[31]提出的公式模擬水平應(yīng)力條件下有效彈性模量(Eeff,Keff):

(34)

(35)

式中:vdra和Edra可以利用彈性模量之間的轉(zhuǎn)換公式進(jìn)行計(jì)算。

(36)

(37)

綜合考慮3種因素,靜態(tài)模量可以由動(dòng)態(tài)模量估算得到。首先計(jì)算頻散效應(yīng),從測(cè)井頻率轉(zhuǎn)換到低頻極限的準(zhǔn)靜態(tài)模量Elf,然后在準(zhǔn)靜態(tài)模量的基礎(chǔ)上計(jì)算排水后的模量Edra,最后考慮應(yīng)變振幅效應(yīng)獲得靜態(tài)模量Esta,流程如圖3所示。

圖3 動(dòng)靜態(tài)模量轉(zhuǎn)換流程

2 實(shí)際算例

選取川西地區(qū)須家河組致密砂巖,地層深度為4650～4730m,以中砂、細(xì)砂巖為主,組成礦物為長(zhǎng)石、泥質(zhì)、石英和巖屑等,孔隙度在1%～9%。由于地層的壓實(shí)作用導(dǎo)致孔隙度和滲透率低,巖石較為致密,孔隙裂縫結(jié)構(gòu)復(fù)雜,大多處于不連通的狀態(tài),如圖4 所示。

這一深度下的致密砂巖裂縫發(fā)育,在地層壓力的作用下,水平方向的裂縫趨于閉合,垂直方向的裂縫趨于張開,導(dǎo)致速度的各向異性。

根據(jù)本文的巖石物理建模流程,計(jì)算動(dòng)態(tài)楊氏模量,圖5為基于致密砂巖巖石物理等效模型的橫波速度和彈性參數(shù)預(yù)測(cè)結(jié)果,其中圖5a中,藍(lán)線為實(shí)際測(cè)井所得的縱波速度,紅線為致密砂巖巖石物理模型預(yù)測(cè)的縱波速度,兩者匹配較高,誤差在5%以內(nèi),驗(yàn)證了模型的準(zhǔn)確性,圖5b為預(yù)測(cè)的橫波速度和動(dòng)態(tài)楊氏模量。

圖5 實(shí)際縱波速度與巖石物理建模結(jié)果及誤差(a)和預(yù)測(cè)的橫波速度與動(dòng)態(tài)楊氏模量(b)曲線

利用致密砂巖巖石物理模型計(jì)算出動(dòng)態(tài)模量,分別分析頻散效應(yīng)、排水效應(yīng)、應(yīng)變振幅對(duì)動(dòng)態(tài)模量的影響如圖6、圖7和圖8所示,其中頻散效應(yīng)和排水效應(yīng)在動(dòng)靜態(tài)轉(zhuǎn)換過程中影響權(quán)重較小,應(yīng)變振幅在轉(zhuǎn)換過程中起主導(dǎo)作用。

圖6 頻散效應(yīng)的影響

圖7 排水效應(yīng)的影響

圖8 應(yīng)變振幅效應(yīng)的影響

綜上分析,采用如圖3所示的動(dòng)靜態(tài)模量轉(zhuǎn)換流程進(jìn)行轉(zhuǎn)換,結(jié)果如圖9所示。

圖9 動(dòng)靜態(tài)模量的轉(zhuǎn)換結(jié)果

圖10為須家河組致密砂巖動(dòng)靜態(tài)模量預(yù)測(cè)結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比,可見兩者具有相同的趨勢(shì)。圖11顯示了預(yù)測(cè)值與實(shí)驗(yàn)值對(duì)比,其中圖11a為動(dòng)態(tài)模量的對(duì)比,兩者相對(duì)誤差約為4%,圖11b為靜態(tài)模量的對(duì)比,兩者相對(duì)誤差為8%,證明本文方法有較好的準(zhǔn)確性。

圖10 預(yù)測(cè)結(jié)果與實(shí)驗(yàn)結(jié)果的對(duì)比

圖11 動(dòng)態(tài)模量預(yù)測(cè)值和實(shí)驗(yàn)值對(duì)比(a)和靜態(tài)模量預(yù)測(cè)值和實(shí)驗(yàn)值對(duì)比(b)

3 結(jié)論

本文針對(duì)致密砂巖,構(gòu)建了適用的巖石物理模型及動(dòng)靜態(tài)模量的轉(zhuǎn)換流程,獲得以下結(jié)論。

1) 巖石動(dòng)靜態(tài)模量的轉(zhuǎn)換機(jī)制是復(fù)雜的,需要考慮頻散效應(yīng)、排水效應(yīng)、應(yīng)變振幅的影響。通過構(gòu)建致密砂巖巖石物理模型,計(jì)算測(cè)井尺度下的動(dòng)態(tài)模量,再考慮頻散效應(yīng)和排水效應(yīng)得到低頻排水條件下的模量,最后根據(jù)動(dòng)靜態(tài)模量之間應(yīng)變振幅的差異,實(shí)現(xiàn)靜態(tài)模量的預(yù)測(cè)。

2) 分析發(fā)現(xiàn)應(yīng)變振幅的影響是造成動(dòng)靜態(tài)模量差異的主要原因,頻散效應(yīng)與排水效應(yīng)的影響相對(duì)較小,其本質(zhì)與巖石中裂縫和微裂縫的發(fā)育有關(guān)。

3) 通過與實(shí)驗(yàn)數(shù)據(jù)對(duì)比,預(yù)測(cè)的靜態(tài)模量具有相同的變化趨勢(shì),誤差在合理的范圍內(nèi)。相較于傳統(tǒng)的動(dòng)靜態(tài)模量轉(zhuǎn)換方法,考慮了動(dòng)靜態(tài)模量轉(zhuǎn)換過程中物理因素的方法,其結(jié)果更為合理。

4) 本文方法只考慮了巖石壓縮過程中裂縫閉合對(duì)巖石彈性性質(zhì)的影響,忽略了孔隙變形的影響。