套損井與取心井相似井段識別及其巖石力學(xué)參數(shù)確定方法

邢岳堃， 張廣清， 李世遠(yuǎn)， 王元元， 楊 瀟

(中國石油大學(xué)(北京)石油工程學(xué)院，北京 102249)

?鉆井完井?

套損井與取心井相似井段識別及其巖石力學(xué)參數(shù)確定方法

邢岳堃， 張廣清， 李世遠(yuǎn)， 王元元， 楊 瀟

(中國石油大學(xué)(北京)石油工程學(xué)院，北京 102249)

為了確定老油田巖心樣本嚴(yán)重匱乏的高套損率區(qū)塊巖石靜態(tài)彈性力學(xué)參數(shù)(靜態(tài)彈性模量與泊松比)，研究了套損井與取心井相似井段識別及其巖石靜態(tài)彈性力學(xué)參數(shù)確定方法。對取心井段的巖心進行巖石力學(xué)試驗，通過分析該井段測井?dāng)?shù)據(jù)與巖石靜態(tài)彈性力學(xué)參數(shù)的相關(guān)性，提出了基于組合模量確定靜態(tài)彈性力學(xué)參數(shù)的模型；通過分析取心井與套損井測井?dāng)?shù)據(jù)的的相關(guān)性，提出了識別取心井與套損井巖石力學(xué)參數(shù)相似井段的判別準(zhǔn)則；建立了通過識別套損井與取心井相似井段來確定套損井段巖石靜態(tài)彈性力學(xué)參數(shù)的方法(SIIM)。由試驗與計算結(jié)果可知：與傳統(tǒng)計算模型中動、靜態(tài)彈性力學(xué)參數(shù)的相關(guān)性相比，考慮有效水平地應(yīng)力、動態(tài)彈性力學(xué)參數(shù)的組合模量與靜態(tài)彈性力學(xué)參數(shù)的相關(guān)系數(shù)分別提高了10%與20%；取心井與鄰近套損井相似井段的識別標(biāo)準(zhǔn)為兩井段的聲波時差與自然伽馬相關(guān)系數(shù)最大和差值數(shù)組標(biāo)準(zhǔn)差最小。研究結(jié)果表明，SIIM方法在高套損率區(qū)塊多井、多井段靜態(tài)彈性力學(xué)參數(shù)的確定中具有較高的適用性與準(zhǔn)確性，在現(xiàn)場應(yīng)用中可取得較好的效果。

彈性模量；泊松比；測井；取心井；套損井；相似井段

造成井下套管損傷的原因十分復(fù)雜，國內(nèi)外相關(guān)研究表明[1-5],注采過程中圍巖靜態(tài)彈性力學(xué)參數(shù)對套管擠壓作用的影響很大，因此獲得套損井段巖石靜態(tài)彈性力學(xué)參數(shù)對制定針對性的套損防治措施意義重大。求取巖石靜態(tài)彈性力學(xué)參數(shù)的方法主要有2種，一是在實驗室對巖樣進行實測，二是用地球物理測井資料求取巖石的動態(tài)彈性力學(xué)參數(shù)，然后應(yīng)用相關(guān)的轉(zhuǎn)換模型轉(zhuǎn)換為靜態(tài)彈性力學(xué)參數(shù)[6]。在目標(biāo)井段取心進行巖石力學(xué)試驗，是獲得巖石靜態(tài)彈性力學(xué)參數(shù)直接且準(zhǔn)確的方法，但由于老油田高套損率區(qū)塊套損井?dāng)?shù)量多且?guī)r心樣本嚴(yán)重匱乏，該方法難以實施。相比于試驗獲取巖石力學(xué)參數(shù)的局限性，應(yīng)用地球物理測井資料計算巖石靜態(tài)彈性力學(xué)參數(shù)的方法[7-9]具有測井資料獲取容易及能夠連續(xù)表征地層信息的優(yōu)點。然而，該方法的準(zhǔn)確性很大程度上取決于動靜態(tài)彈性力學(xué)參數(shù)轉(zhuǎn)換模型[10-11]在目標(biāo)井段的適用性，且轉(zhuǎn)換模型需要通過巖石力學(xué)試驗進行標(biāo)定，對于巖心嚴(yán)重匱乏的套損井，在鄰井取心進行巖石力學(xué)試驗標(biāo)定獲得的動靜態(tài)彈性力學(xué)參數(shù)轉(zhuǎn)換模型是否適用于套損井段，尚需進一步研究論證。綜合以上2種方法的優(yōu)點，筆者提出在取心井相似井段(與鄰井多口套損井段相似)取心進行巖石力學(xué)參數(shù)測試，將巖石力學(xué)參數(shù)與測井?dāng)?shù)據(jù)的轉(zhuǎn)換模型應(yīng)用于套損井段的方法。其中，包含2個關(guān)鍵問題：1)如何確定適用于高套損率區(qū)塊多井、多井段的巖石力學(xué)參數(shù)(尤其是靜態(tài)彈性力學(xué)參數(shù))與測井?dāng)?shù)據(jù)相關(guān)關(guān)系的計算模型；2)如何建立準(zhǔn)確有效的相似井段識別方法。

聲波測井?dāng)?shù)據(jù)可以表征巖石的動態(tài)彈性力學(xué)參數(shù)，而靜態(tài)彈性力學(xué)參數(shù)能更準(zhǔn)確地描述套損井段巖石的力學(xué)特性。W.A.Zisman[12]早在1933年就指出巖石動靜態(tài)彈性力學(xué)參數(shù)之間存在差異，因此確定套損井段巖石動、靜態(tài)彈性力學(xué)參數(shù)的相關(guān)關(guān)系，是建立測井?dāng)?shù)據(jù)與巖石靜態(tài)彈性力學(xué)參數(shù)轉(zhuǎn)換模型的關(guān)鍵；N.R.Warpinski等人[13]對比了聲波測井解釋值與靜態(tài)測試值，但沒有提出包含地應(yīng)力的動靜態(tài)參數(shù)轉(zhuǎn)換模型；國內(nèi)外許多學(xué)者[14-15]研究了巖石動、靜態(tài)彈性力學(xué)參數(shù)之間的關(guān)系，發(fā)現(xiàn)巖石動、靜態(tài)彈性模量之間具有較好的相關(guān)性，但動、靜態(tài)泊松比之間的關(guān)系不明顯；樓一珊等人[16]提出了巖石靜態(tài)彈性模量與泊松比在主應(yīng)力差確定的情況下，分別與動態(tài)彈性模量及泊松比呈線性關(guān)系的轉(zhuǎn)換模型，并廣泛用于石油工程領(lǐng)域，但地層真實的主應(yīng)力差難以確定。以上研究存在共同的制約因素，即通過試驗標(biāo)定的動靜態(tài)彈性力學(xué)參數(shù)轉(zhuǎn)換模型，若計算井段過長會有很大的計算誤差，且不適用于多井、多井段的計算。從油藏地質(zhì)特征角度看，相似井段的識別與地層的劃分、對比有一定的相似性，其中基于旋回對比法則的高分辨率層序地層分析[17-21]、模式識別[22]等是地層劃分的有效方法；在高分辨率層序格架內(nèi)開展小層對比[23]是獲取等時對比層的常用方法。上述方法主要用于描述儲層，且小層對比獲得相似地層的方法受人為因素的影響，因而該方法不能解決高套損率區(qū)塊套損井段數(shù)量巨大與出現(xiàn)位置隨機的問題，更難以實現(xiàn)識別套損井與取心井的相似井段、獲取套損井段靜態(tài)彈性力學(xué)參數(shù)的最終目的。

針對以上2個關(guān)鍵問題，筆者基于巖石力學(xué)試驗和測井?dāng)?shù)據(jù)分析，提出了基于組合模量(考慮有效水平地應(yīng)力與動態(tài)彈性力學(xué)參數(shù))計算靜態(tài)彈性力學(xué)參數(shù)的模型(下文簡稱為組合模量模型)；進而分析取心井與套損井聲波時差、自然伽馬測井曲線的相關(guān)性，建立了在取心井查找與套損井段巖石力學(xué)參數(shù)相似井段的識別準(zhǔn)則；最終將組合模量模型應(yīng)用于套損井段，較好地解決了巖心匱乏高套損率區(qū)塊獲得多井、多井段巖石力學(xué)參數(shù)的問題。

1 套損井、取心井及測井資料的選取

選取西部某油田某區(qū)塊巖心資料匱乏的3口套損井，分別記為S1井、S2井和S3井。套損井信息見表1。

表1 套損井信息Table 1 Data of casing failure wells

取心井為Q1井，S1井、S2井、S3井與Q1井的井距分別為190.00，245.00和275.00 m。

由于巖石中聲波速度的大小取決于彈性模量的大小[24]，巖石總的伽馬射線強度可用于研究井剖面的地層性質(zhì)[25]。因此，筆者選擇分析聲波時差和自然伽馬2種測井?dāng)?shù)據(jù)，以獲得巖石的物性與巖性信息，用于相似井段的識別。

2 巖石靜態(tài)彈性力學(xué)參數(shù)與測井?dāng)?shù)據(jù)的相關(guān)性

考慮套損井段圍巖所處地層的井深及上覆巖層壓力的差異，巖石的靜態(tài)彈性力學(xué)參數(shù)不僅與由測井?dāng)?shù)據(jù)獲得的動態(tài)彈性力學(xué)參數(shù)相關(guān)，還與所處的應(yīng)力狀態(tài)相關(guān)[16,24]，提出了組合模量(綜合水平地應(yīng)力與動態(tài)彈性力學(xué)參數(shù))與靜態(tài)彈性力學(xué)參數(shù)的相關(guān)計算模型，并對組合模量的表達形式及與靜態(tài)彈性力學(xué)參數(shù)的相關(guān)性進行了研究。

2.1 巖石力學(xué)試驗

Q1井巖心存放時間較長，由于長期取樣，全直徑巖心破壞嚴(yán)重。因此，以距離相似井段最近、取心數(shù)量最多為原則，篩選出9組標(biāo)準(zhǔn)巖心試樣(長為50.0 mm，直徑為25.0 mm)，用三軸伺服巖石壓縮試驗機先以3.0 MPa/min的速率施加圍壓，然后以0.02 mm/min的軸向變形進行三軸壓縮試驗，結(jié)果見表2。

表2 巖石力學(xué)試驗結(jié)果Table 2 Test results for rock mechanics

2.2 考慮水平地應(yīng)力的靜態(tài)彈性力學(xué)參數(shù)計算模型

巖石動靜態(tài)彈性力學(xué)參數(shù)的相關(guān)性分析,是通過測井?dāng)?shù)據(jù)預(yù)測巖石靜態(tài)彈性力學(xué)參數(shù)的一項重要內(nèi)容。石油工程領(lǐng)域廣泛使用的一種動靜態(tài)彈性力學(xué)參數(shù)轉(zhuǎn)化關(guān)系(后文簡稱該轉(zhuǎn)化關(guān)系為傳統(tǒng)模型)為：巖石靜態(tài)彈性模量與泊松比在主應(yīng)力差確定的情況下，分別與動態(tài)彈性模量及泊松比呈線性關(guān)系[16]，其表達式為：

(1)

式中：Ed和Es分別為動、靜態(tài)彈性模量，GPa;νd和νs分別為動、靜態(tài)泊松比；a1與k1為巖石動、靜態(tài)彈性模量線性關(guān)系式的截距與斜率;a2與k2為巖石動、靜態(tài)泊松比線性關(guān)系式的截距與斜率。

然而，傳統(tǒng)模型在解決套損井段巖石力學(xué)參數(shù)方面存在2個主要問題：1)套損井段主應(yīng)力差難以確定；2)因巖樣的差異性，線性函數(shù)的截距和斜率往往與主應(yīng)力差對數(shù)值的相關(guān)性很差。

若將地層視為橫觀各向同性巖層，有效圍壓為水平有效地應(yīng)力分量，在巖石力學(xué)三軸試驗中圍壓與巖石抗壓強度、彈性模量呈正相關(guān)關(guān)系[24]。筆者在現(xiàn)有研究的基礎(chǔ)上，提出了用組合模量計算靜態(tài)彈性力學(xué)參數(shù)的模型，其中組合模量α和β分別與Es和νs呈線性關(guān)系，其表達式為：

(2)

(3)

式中：α為與動態(tài)彈性模量、有效圍壓相關(guān)的組合模量，GPa；β為與動態(tài)泊松比、有效圍壓相關(guān)的組合模量，GPa；C1和M1為巖石組合模量α與靜態(tài)彈性模量線性關(guān)系式的截距和斜率；C2和M2為巖石組合模量β與靜態(tài)泊松比線性關(guān)系式的截距和斜率;pc為三軸試驗中的圍壓(對應(yīng)于橫觀各向同性巖層的水平有效地應(yīng)力分量)，GPa。

選用3組試驗數(shù)據(jù)分析靜態(tài)彈性力學(xué)參數(shù)與組合模量的相關(guān)性(A組與B組為國內(nèi)學(xué)者[15,26]研究得到的砂巖三軸試驗數(shù)據(jù)，Q組為表2中的試驗結(jié)果)，結(jié)果如圖1所示。由圖1可知，3組試驗所得結(jié)果均顯示，靜態(tài)彈性力學(xué)參數(shù)與組合模量具有很高的相關(guān)性。

將組合模量模型的相關(guān)系數(shù)與式(1)表示的傳統(tǒng)模型進行比較(見表3)，組合模量模型中Es與α以及νs與β的R2分別比傳統(tǒng)模型高10%與20%。傳統(tǒng)模型得到的動、靜態(tài)泊松比出現(xiàn)了正負(fù)相關(guān)關(guān)系不確定的情況，而組合模量與靜態(tài)彈性力學(xué)參數(shù)始終為正相關(guān)關(guān)系，解決了動靜態(tài)泊松比相關(guān)關(guān)系難以確定的問題。

圖1 靜態(tài)彈性力學(xué)參數(shù)與組合模量的相關(guān)性Fig.1 Correlation among static elastic parameters and combined modulus

Table 3 Comparison between combined modulus model and conventional model

模型類型相關(guān)系數(shù)A組B組Q組平均值組合模量模型R2(Es，α)090092083088R2(νs，β)094073087085傳統(tǒng)模型R2(Es，Ed)088081071080R2(νs，νd)050079?085071

注：*為負(fù)相關(guān)關(guān)系，A組、B組、C組中其余皆為正相關(guān)關(guān)系。

組合模量考慮了地應(yīng)力對巖石的擠壓作用,因此與巖石彈性力學(xué)參數(shù)具有更好的相關(guān)性。如果水平地應(yīng)力難以確定，在求解組合模量模型系數(shù)時可用有效上覆巖層壓力代替pc。綜上所述，組合模量模型具有很好的適用性。

3 巖石力學(xué)參數(shù)相似井段的識別

3.1 井段相似識別標(biāo)準(zhǔn)及計算方法

測井?dāng)?shù)據(jù)與巖石靜態(tài)彈性力學(xué)參數(shù)的相關(guān)關(guān)系體現(xiàn)在組合模量模型中，因此巖性、物性測井?dāng)?shù)據(jù)(如聲波時差和自然伽馬測井)的變化與巖石力學(xué)參數(shù)的變化一一對應(yīng)?；跀?shù)理統(tǒng)計，筆者提出對于取心井與鄰近套損井的相似井段，兩者測井?dāng)?shù)據(jù)的(聲波時差與自然伽馬)相關(guān)系數(shù)最大和差值數(shù)組標(biāo)準(zhǔn)差最小是巖石力學(xué)參數(shù)相似的識別標(biāo)準(zhǔn)，相似井段計算流程如圖2所示。

圖2 套損井與取心井相似井段計算流程Fig.2 Process flow in calculation for similar intervals in wells with casing failure and coring wells

套損井段聲波時差(ac)和自然伽馬(gc)數(shù)據(jù)點的總數(shù)均為n1，取心井全井段聲波時差(a)和自然伽馬(g)數(shù)據(jù)點的總數(shù)均為n2，且n2>n1。將ac，gc分別與a，g按深度自淺至深分別進行n2-n1次差值計算，則有：

(4)

(5)

式中：i=1,2,3,…,n1；j=1,2,3,…,n2-n1+1。

每次計算完畢，得到n2-n1個長度為n1的差值數(shù)組ad與gd。考慮ad與gd的數(shù)量級差異，引入削弱系數(shù)am和gm，降低因數(shù)量級差異對計算結(jié)果的影響，則有：

(6)

(7)

式中:上標(biāo)“-”表示取平均值。下同。

經(jīng)過n2-n1次計算獲得每個測井?dāng)?shù)據(jù)差值數(shù)組的標(biāo)準(zhǔn)差與測井?dāng)?shù)據(jù)的相關(guān)系數(shù)，并分別形成數(shù)組σ與r，可表示為：

(9)

式中：i=1,2,3,…,n1；j=1,2,3,…,n2-n1+1。

若標(biāo)準(zhǔn)差數(shù)組σ的最小值min[σ]與相關(guān)系數(shù)數(shù)組r的最大值max[r]對應(yīng)的取心井段相重合，可初步判定兩個井段的巖石力學(xué)參數(shù)相似。

3.2 計算結(jié)果與分析

計算可得取心井各深度對應(yīng)的標(biāo)準(zhǔn)差以及相關(guān)系數(shù)，min[σ]與max[r]對應(yīng)井段上下15.00 m的計算結(jié)果如圖3所示，其中每個數(shù)據(jù)點對應(yīng)的深度為每次參與循環(huán)計算數(shù)組的起始深度。圖3中，用黑色方框標(biāo)注max[r]與min[σ]的數(shù)據(jù)點。

圖3 取心井與套損井相似井段識別標(biāo)準(zhǔn)計算結(jié)果Fig.3 Calculation results of identification criterion for similar intervals in casing failure wells and coring wells

由圖3可知:在相似地層上下15.00 m的區(qū)間內(nèi)，聲波時差與自然伽馬的相關(guān)系數(shù)曲線波動趨勢基本一致，且2條曲線峰值點與谷值點對應(yīng)于同一深度；聲波時差與自然伽馬的標(biāo)準(zhǔn)差曲線波動趨勢基本一致，且2條曲線最低點大致對應(yīng)同一深度；且max[r]與min[σ]的數(shù)據(jù)點均對應(yīng)同一取心井段。

取心井與套損井相似地層對照情況見表4。

表4 取心井與套損井相似井段對比Table 4 Data of similar intervals in casing failure wells and coring wells

由表4可知，3口套損井均可在取心井中找到與套損井段處于同一地質(zhì)年代、地質(zhì)分層且測井?dāng)?shù)據(jù)相關(guān)性很高的取心井段，并分別繪制相似地層的測井曲線(見圖4)。S1井、S2井和S3井均可在取心井Q1井中找到與套損井段測井曲線旋回趨勢基本一致的相似井段。因此，可用取心井中相似井段的巖心進行巖石力學(xué)試驗(見表2)，進而確定組合模量模型中的系數(shù)，用于套損井段靜態(tài)彈性力學(xué)參數(shù)的計算。

4 套損井段巖石靜態(tài)彈性力學(xué)參數(shù)計算

筆者應(yīng)用SIIM方法計算了S1井、S2井和S3井等3口井套損井段的靜態(tài)彈性力學(xué)參數(shù)，主要計算步驟為：

1) 根據(jù)式(4)—式(9)所述的識別準(zhǔn)則，利用取心井的聲波與自然伽馬測井?dāng)?shù)據(jù),分析取心井與S1井、S2井和S3井套損井段巖石力學(xué)參數(shù)相似的井段(見表4)。

2) 在取心井的相似井段取心進行巖石力學(xué)試驗，結(jié)果見表2。

3) 根據(jù)巖石力學(xué)試驗結(jié)果與測井?dāng)?shù)據(jù)確定組合模量模型中的系數(shù)，使該模型同時適用于取心井段與套損井的相似井段。橫波測井?dāng)?shù)據(jù)往往較難獲得，該情況下可通過巖石剪切模量與巖石縱波模量之間的關(guān)系[27-28]獲得干巖石的橫波波速，計算式為：

(10)

(11)

式中：Md為干巖石縱波模量，10-3g/(mm·s2)；vp與vs分別為縱波波速與橫波波速，mm/s；ρ為巖樣密度，g/cm3。

進而可將測井?dāng)?shù)據(jù)轉(zhuǎn)換為動態(tài)彈性力學(xué)參數(shù)，計算式為：

(12)

(13)

圖4 套損井與取心井相似井段測井曲線Fig.4 Well logs for similar intervals between casing failure wells and coring well

結(jié)合巖石力學(xué)試驗結(jié)果確定的組合模量模型中的系數(shù)，將式(3)、式(10)—式(13)代入式(2)，可得：

(14)

(15)

若實際應(yīng)用中水平地應(yīng)力難以確定，在確定組合模量模型中的系數(shù)時，pc可選用有效上覆巖層壓力，并用于最后的套損井段靜態(tài)彈性力學(xué)參數(shù)計算。

4) 套損井段巖石靜態(tài)彈性力學(xué)參數(shù)的確定。利用式(14)和式(15)將S1井、S2井和S3井等3口井套損井段的測井?dāng)?shù)據(jù)轉(zhuǎn)換為靜態(tài)彈性力學(xué)參數(shù)，并獲得如圖5所示的參數(shù)剖面。

5 現(xiàn)場應(yīng)用

SIIM方法可通過編寫計算機程序應(yīng)用于巖心樣本匱乏的高套損率區(qū)塊，以解決確定多井、多井段巖石力學(xué)參數(shù)的問題。例如,A1—A5為5個高套損率區(qū)塊，每個區(qū)塊均有大量的套損井及有限的取心井。其中，A1、A2和A3等3個區(qū)塊由于投產(chǎn)時間較早，套損尤其嚴(yán)重且取心井中可以取樣的巖心極其匱乏，長期以來套損井巖石力學(xué)參數(shù)的預(yù)測因缺乏切實有效的方法未能大規(guī)模開展，影響了針對性套損預(yù)防與治理措施的制定與實施。在A1—A5等5個高套損率區(qū)塊應(yīng)用SIIM方法選取24口取心井，對397口套損井進行了相似井段的識別與巖石力學(xué)參數(shù)計算，應(yīng)用效果見表5。

由表5可知，該方法適用于344口套損井，總體適用率為86.6%，且各區(qū)塊砂巖套損井段適用率高于泥巖套損井段。因此，SIIM方法克服了油田高套損率區(qū)塊多井、多井段巖石靜態(tài)彈性力學(xué)參數(shù)無法大規(guī)模預(yù)測的難題，有助于制定套管損傷預(yù)防及治理措施，降低開采成本。

圖5 套損井段靜態(tài)彈性力學(xué)參數(shù)剖面Fig.5 Profiles of static elastic modulus for intervals with casing damage

Table 5 Performances of SIIM in blocks with high rates of casing damage

區(qū)塊編號取心井?dāng)?shù)套損井?dāng)?shù)應(yīng)用井?dāng)?shù)泥巖砂巖泥巖砂巖A1569545749A2693488143A3650424138A45317215A52021018合計24215182181163

6 結(jié) 論

1) 與傳統(tǒng)計算模型中動、靜態(tài)彈性力學(xué)參數(shù)的相關(guān)性相比，考慮有效水平地應(yīng)力、動態(tài)彈性力學(xué)參數(shù)的組合模量與靜態(tài)彈性力學(xué)參數(shù)(靜態(tài)彈性模量與泊松比)的相關(guān)性分別提高了10%與20%，在計算靜態(tài)彈性力學(xué)參數(shù)時具有更高的準(zhǔn)確性與更好的適用性。

2) 對于取心井與鄰近套損井，巖石力學(xué)參數(shù)相似井段的識別標(biāo)準(zhǔn)為兩井段的聲波時差與自然伽馬相關(guān)系數(shù)最大和差值數(shù)組標(biāo)準(zhǔn)差最小。

3) 通過識別套損井與取心井相似井段確定套損井巖石靜態(tài)彈性力學(xué)參數(shù)的方法，在巖心樣本匱乏的條件下，能夠較為準(zhǔn)確地獲得高套損率區(qū)塊多井、多井段巖石力學(xué)靜態(tài)彈性力學(xué)參數(shù)。

4) SIIM方法通過編寫計算機程序可廣泛應(yīng)用于巖心樣本匱乏的高套損率區(qū)塊，總體適用率為86.6%，有助于制定針對性的套損防治措施，降低開采成本。

[1] GHOLAMI R,RASOULI V,AADNOY B,et al.Geomechanical and numerical studies of casing damages in a reservoir with solid production[J].Rock Mechanics & Rock Engineering,2016,49(4):1441-1460.

[2] 宋杰.砂巖油藏地應(yīng)力及巖石力學(xué)參數(shù)與套管損壞相關(guān)性[J].地質(zhì)力學(xué)學(xué)報，2014，20(3)：324-330. SONG Jie.The correlations between geological stress and casing damage and between rock mechanics parameters and casing damage in sandstone reservoir[J].Journal of Geomechanics,2014,20(3):324-330.

[3] XING Y K,ZHANG G Q,LI S Y,et al.Investigation on the integrity of cement annulus during multi-stage horizontal fracturing with FEM/DEM[C]//Proceedings of the 7th International Conference on Discrete Element Methods,December 02,2016.Singapore:Springer,2016.

[4] 張廣清，陳勉，王群立.泥頁巖地層水井套管載荷、臨界注水時間確定[J].哈爾濱工業(yè)大學(xué)學(xué)報，2006，38(8)：1273-1277. ZHANG Guangqing,CHEN Mian,WANG Qunli.Study on casing load in injections wells of shale formation and the critical injection time[J].Journal of Harbin Institute of Technology,2006,38(8):1273-1277.

[5] 張廣清，陳勉，耿宇迪，等.松軟地層套管損壞的影響因素研究[J].石油勘探與開發(fā)，2005，32(1)：112-115. ZHANG Guangqing,CHEN Mian,GENG Yudi,et al.Casing damage mechanism in sand producing formations[J].Petroleum Exploration & Development,2005,32(1):112-115.

[6] 路保平，鮑洪志.巖石力學(xué)參數(shù)求取方法進展[J].石油鉆探技術(shù)，2005，33(5)：44-47. LU Baoping,BAO Hongzhi.Advances in calculation methods for rock mechanics parameters[J].Petroleum Drilling Techniques,2005,33(5):44-47.

[7] 尹帥，單鈺銘，周文，等.Hoek-Brown準(zhǔn)則在致密砂巖彈性力學(xué)參數(shù)測井解釋中的應(yīng)用[J].石油鉆探技術(shù)，2015，43(1)：88-95. YIN Shuai,SHAN Yuming,ZHOU Wen,et al.Application of Hoek-Brown criterion for tight sandstone elastic parameters in log interpretation[J].Petroleum Drilling Techniques,2015,43(1):88-95.[8] 歸榕，萬永平.基于常規(guī)測井?dāng)?shù)據(jù)計算儲層巖石力學(xué)參數(shù):以鄂爾多斯盆地上古生界為例[J].地質(zhì)力學(xué)學(xué)報，2012，18(4)：418-424. GUI Rong,WAN Yongping.Rock mechanics parameter calculation based on conventional logging data:a case study of Upper Paleozoic in Ordos Basin[J].Journal of Geomechanics,2012,18(4):418-424.

[9] 賀順義，師永民，謝楠，等.根據(jù)常規(guī)測井資料求取巖石力學(xué)參數(shù)的方法[J].新疆石油地質(zhì)，2008，29(5)：662-664. HE Shunyi,SHI Yongmin,XIE Nan,et al.The method for acquirement of conventional logging response-based lithomechanical parameters[J].Xinjiang Petroleum Geology,2008,29(5):662-664.

[12] ZISMAN W A.Comparison of the statically and seismologically determined elastic constants of rocks[J].Proceedings of the National Academy of Sciences of the United States of America,1933,19(7):680-686.

[13] WARPINSKI N R,PETERSON R E,BRANAGAN P T,et al.In situ stress and moduli:comparison of values derived from multiple techniques[R].SPE 49190,1998.

[14] JIZBA D,MAVKO G,NUR A.Static and dynamic moduli of tight gas sandstones[M].Tulsa： Society of Exploration Geophysicists,1990:827-829.

[15] 邊會媛，王飛，張永浩，等.儲層條件下致密砂巖動靜態(tài)彈性力學(xué)參數(shù)實驗研究[J].巖石力學(xué)與工程學(xué)報，2015，34(增刊1)：3045-3054. BIAN Huiyuan,WANG Fei,ZHANG Yonghao,et al.Experimental study of dynamic and static elastic parameters of tight sandstones under reservoir conditions[J].Chinese Journal of Rock Mechanics & Engineering,2015,34(supplement 1):3045-3054.

[16] 樓一珊，莊錦江，黃榮樽.巖石動、靜彈性力學(xué)參數(shù)相關(guān)性的研究及其應(yīng)用[J].江漢石油學(xué)院學(xué)報，1989，11(2)：62-69. LOU Yishan,ZHUANG Jinjiang,HUANG Rongzun.Study on the interrelation of dynamic and static elastic parameters of rocks and its application[J].Journal of Jianghan Petroleum Institute,1989,11(2):62-69.

[17] CROSS T A.Controls on coal distribution in transgressive-regressive caycles,upper cretaceous,western interior,U.S.A.[J].Special Publications,1988,42:371-380.

[18] CROSS T A,LESSENGER M A.Sediment volume partitioning:rationale for stratigraphic model evaluation and high-resolution stratigraphic correlation[C]//GRADSTEIN F M,SANDVIK K O,MILTON N J,et al. Sequence strati graphy concepts and application.Stavange:Norwegian Petroleum Society,1998:171-195.

[19] CROSS T A.Stratigraphic controls on reservoir attributes in continental strata[J].Earth Science Frontiers,2000,7(4):322-350.

[20] AITKEN J F,HOWELL J A.High resolution sequence stratigraphy:innovations,applicationsand future prospects[J].Geological Society London Special Publications,1996,104(1):1-9.

[21] van WAGONER J C,MITCHUM R M,CAMPION K M,et al.Siliciclastic sequence stratigraphy in well logs,cores,and outcrops:concepts for high-resolution correlation of time and facies[M].Tulsa:the American Association of Petroleum Geologists,1990:43-62.

[22] 耿耀輝，許少華，劉洪林.模式識別在測井曲線地層對比中的應(yīng)用[J].大慶石油學(xué)院學(xué)報，1995，19(4)：34-37. GENG Yaohui,XU Shaohua,LIU Honglin.Application of pattern recognition to stratigraphic correlation of well-logging curve[J].Journal of Daqing Petroleum Institute,1995,19(4):34-37.

[23] 趙翰卿.高分辨率層序地層對比與我國的小層對比[J].大慶石油地質(zhì)與開發(fā)，2005，24(1)：5-9，12. ZHAO Hanqing.High-resolution sequential stratigraphy correlation and Chinese subzone correlation[J].Petroleum Geology & Oilfield Development in Daqing,2005,24(1):5-9,12.

[24] FJAER E,HOLT R M,HORSRUD P,et al.Petroleum related rock mechanics[M].2nd ed.Amsterdam:Elsevier,2008:57，175.

[25] 洪有密.測井原理與綜合解釋[M].東營：石油大學(xué)出版社，1993：195. HONG Youmi.Principle and comprehensive interpretation of well logging[M].Dongying:Petroleum University Press,1993:195.

[26] 王艷梅.巖石力學(xué)動、靜態(tài)參數(shù)關(guān)系研究[D].成都：西南石油大學(xué)，2012. WANG Yanmei.Study on static and dynamic parameters of rock mechanics[D].Chengdu:Southwest Petroleum University,2012.

[27] MAVKO G,CHAN C,MUKERJI T.Fluid substitution:estimating changes in VP without knowing VS[J].Geophysics,1995,60(6):1750-1755.

[28] 張元中，周開金，趙建斌，等.砂泥巖地層橫波測井曲線預(yù)測方法研究[J].石油物探，2012，51(5)：508-514. ZHANG Yuanzhong,ZHOU Kaijin,ZHAO Jianbin,et al.Shear-wave logging curve prediction method for shaly sand formation[J].Geophysical Prospecting for Petroleum,2012,51(5):508-514.

[編輯 令文學(xué)]

Identification of Similar Intervals of Wells with Casing Failure and Coring Wells and the Determination of their Geomechanical Properties

XING Yuekun, ZHANG Guangqing, LI Shiyuan, WANG Yuanyuan, YANG Xiao

(CollegeofPetroleumEngineering,ChinaUniversityofPetroleum(Beijing),Beijing, 102249,China)

To determine the static elastic parameters (SEP) (static elastic modulus and Poisson′s ratios) of rock in blocks with high casing damage rates and with limited core samples, an innovative technique has been proposed to identify similar intervals of wells with casing failure and coring wells and to determine the SEP of these intervals. First, the mechanical parameters of cores from coring intervals were tested, then the correlation between logging data and rock mechanical parameters between the similar intervals of the coring well and the casing failure well was analyzed, a model from combined modulus was established to determine static elastic parameters. Then a similar SEP identification criterion was proposed. Mechanical parameters of coring intervals were determined to establish correlations among the logging data and SEPs of rocks in these intervals. Eventually, models for the SEPs based on combined modulus were proposed. By studying correlation among the logging data acquired in the wells with casing failures and coring wells, criteria for similar intervals in wells with casing failure and cored wells were highlighted. By identifying intervals with similar properties in casing failure wells and cored wells, it is possible to determine SEP in intervals of casing failure. Testing and calculation results showed that the combined modulus and SEP with consideration of effective horizontal in-situ stresses may have correlations 10% and 20% higher than those from conventional calculation models. In general, identification criteria for similar intervals in casing failure wells and coring wells can be summarized as the highest correlation between acoustic travel time and GR, and minimum deviation among difference arrays in these two intervals. Research results showed that these techniques have relatively high applicability and accuracy in determination of SEP in multiple wells and intervals in areas with severe casing damage. Satisfactory performance has been observed in field applications of these techniques.

modulus of elasticity;Poissons ratio;woll logging;coring well;casing failure well;similar intervals

2016-12-19；改回日期：2017-06-25。

邢岳堃(1991—)，男，山東泰安人，2014年畢業(yè)于中國海洋大學(xué)船舶與海洋工程專業(yè)，力學(xué)專業(yè)在讀博士生，研究方向為石油工程巖石力學(xué)和斷裂力學(xué)。E-mail：oucxyk@163.com。

張廣清，zhang0@263.net。

國家自然科學(xué)基金優(yōu)秀青年科學(xué)基金項目“石油工程巖石力學(xué)”(編號：51322404)資助。

10.11911/syztjs.201704006

TE21

A

1001-0890(2017)04-0033-08