楊同玉，吳 超，吳繼偉

(1.中國石化華北分公司，河南 鄭州 450006；2.中國石化石油工程技術(shù)研究院，北京 100101；3.中國石油新疆油田分公司，新疆 克拉瑪依 834000)

0 引 言

鉆井地質(zhì)力學(xué)參數(shù)通常包括巖石變形參數(shù)、抗壓強(qiáng)度、抗剪強(qiáng)度、抗拉強(qiáng)度、可鉆性、構(gòu)造應(yīng)力、孔隙壓力、坍塌壓力、破裂壓力等。精細(xì)描述鉆井地質(zhì)力學(xué)參數(shù)的分布規(guī)律將有助于全面認(rèn)識鉆井破巖能力和井身質(zhì)量控制難點(diǎn)，為優(yōu)化井身結(jié)構(gòu)、井眼軌跡、鉆頭型號及機(jī)械水力參數(shù)、鉆具組合、井控方案、鉆井液性能等工藝措施提供依據(jù)，為實現(xiàn)安全、優(yōu)質(zhì)、高效、低成本鉆井打下堅實基礎(chǔ)[1]。

數(shù)學(xué)地質(zhì)理論認(rèn)為，地質(zhì)參數(shù)在地下三維巖層空間中的分布具有統(tǒng)計相關(guān)性與隨機(jī)性的特征[2-3]?；诘厍蛭锢硇畔⒖山⑷S地層空間的鉆井地質(zhì)力學(xué)參數(shù)模型，從而掌握區(qū)域力學(xué)特征，以合理優(yōu)化目標(biāo)井(尤其是定向井、水平井)鉆井工程技術(shù)方案。近年來，國內(nèi)外學(xué)者基于非線性優(yōu)化理論提出了多種鉆井地質(zhì)力學(xué)參數(shù)預(yù)測方法[4-10]，但多針對井壁圍巖，將其推廣至三維地層空間存在難度。有學(xué)者嘗試將地震層析與儲層建模算法用于力學(xué)參數(shù)三維建模，但現(xiàn)場應(yīng)用表明，其計算精度和分辨率均存在一定限制[11-14]。因此，有必要研發(fā)一套精度高、適用面廣的鉆井地質(zhì)力學(xué)參數(shù)空間建模技術(shù)，使其更加有利于指導(dǎo)鉆井提速提效工作。為此，提出一種基于地震屬性的鉆井地質(zhì)力學(xué)參數(shù)空間建模方法。針對目標(biāo)工區(qū)，利用完鉆井資料建立層速度與地震屬性之間的分層神經(jīng)網(wǎng)絡(luò)模型，利用其預(yù)測三維地層空間的層速度及各類鉆井地質(zhì)力學(xué)參數(shù)，進(jìn)一步確定區(qū)域地層力學(xué)分布特性，并將其用于鉆井工程優(yōu)化設(shè)計。

1 理論基礎(chǔ)

按照地球物理原理，地下界面的反射系數(shù)和地震子波經(jīng)褶積運(yùn)算可得到疊后地震記錄數(shù)據(jù)，即：

(1)

式中：D(i)、R(i)、W(i)分別為地震記錄、反射系數(shù)和地震子波；i為地層深度點(diǎn)序號；*表示褶積運(yùn)算。

根據(jù)勘探地震學(xué)理論，反射系數(shù)的數(shù)值取決于地震波反射界面兩側(cè)地層的波速和密度：

(2)

式中：ρ為地層密度，g/cm3；vp為縱波速度，m/s。

巖石物理實驗研究表明，巖石的縱波速度和密度之間存在較好的經(jīng)驗關(guān)系：

ρ=C0vpC1

(3)

式中：C0、C1為經(jīng)驗參數(shù)。

綜合分析式(1)、(2)、(3)可知，疊后地震記錄和縱波速度之間存在直接的定量關(guān)系。由于縱波速度是計算各類鉆井地質(zhì)力學(xué)參數(shù)的最重要、最基礎(chǔ)的數(shù)據(jù)，因此，可考慮使用地震記錄反演求取地層縱波速度，從而進(jìn)一步預(yù)測地層空間三維區(qū)域內(nèi)的各類力學(xué)參數(shù)。

受復(fù)雜地層條件影響，地震記錄與波速之間存在非線性映射關(guān)系，對于同一工區(qū)的相同地質(zhì)層段，所體現(xiàn)出的上述映射關(guān)系是相近的，因此，可以借鑒人工智能方法，根據(jù)完鉆井?dāng)?shù)據(jù)信息，利用神經(jīng)網(wǎng)絡(luò)算法識別這類映射關(guān)系，并將其推廣至工區(qū)中的未鉆探空間。

2 地震層速度的預(yù)測

地震屬性是從地震記錄中通過特定數(shù)學(xué)方法提取的多種參數(shù)，可從多角度反映地層的構(gòu)造、巖性、物性以及流體成分等特征，相比于原始地震記錄，其更適宜用來建立層速度與地震反射信息之間的非線性映射關(guān)系[15-18]。此次研究提取并使用的地震屬性為：①反映界面反射特征的振幅屬性，包括均方根振幅、平均絕對振幅、相鄰時窗振幅比等；②描述地震信號瞬時變化規(guī)律的瞬時屬性，包括瞬時振幅、瞬時相位、瞬時頻率等；③體現(xiàn)一定時窗內(nèi)地層巖性、構(gòu)造與物性變化的傅立葉譜屬性，包括振幅譜主頻、振幅譜極大值、平均中心頻率等；④反映地震隨機(jī)特征的功率譜屬性，包括加權(quán)平均頻率、功率譜極大頻率等。

人工智能領(lǐng)域已研發(fā)了多種神經(jīng)網(wǎng)絡(luò)方法，其中，小波神經(jīng)網(wǎng)絡(luò)模式識別速度快、容錯能力強(qiáng)[19-20]，因此，此次研究中運(yùn)用小波神經(jīng)網(wǎng)絡(luò)進(jìn)行地震層速度預(yù)測，具體步驟如下。

(1) 設(shè)定小波函數(shù)，其神經(jīng)網(wǎng)絡(luò)輸出yl′為：

(4)

式中：ak、bk為小波函數(shù)因子；ωk1、ωkj2為神經(jīng)網(wǎng)絡(luò)權(quán)值(上標(biāo)表示權(quán)值類別)；xlj表示地震屬性；l=1,2,…，m,m為學(xué)習(xí)樣本對的數(shù)量；j=1,2,…，n,n為所使用地震屬性的種類數(shù)；k=1,2,…，t,t為基函數(shù)個數(shù)。

(2) 將ak、bk、ωk1、ωkj2隨機(jī)設(shè)定初始值。對于目標(biāo)工區(qū)內(nèi)鉆遇的各套地層，將其地震屬性及相應(yīng)的測井聲波速度作為學(xué)習(xí)樣本對，導(dǎo)入神經(jīng)網(wǎng)絡(luò)進(jìn)行映射關(guān)系識別。

(5)

式中：μ為網(wǎng)絡(luò)參數(shù)；p為迭代次數(shù)；γ為學(xué)習(xí)率；ξ為動量項。

按上述計算步驟迭代計算，直至誤差滿足要求，即可完成網(wǎng)絡(luò)模式識別。在此基礎(chǔ)上，將工區(qū)未鉆探空間內(nèi)的地震屬性輸入網(wǎng)絡(luò)，得到各地震道上的層速度，從而建立工區(qū)三維地震層速度數(shù)據(jù)體。

3 鉆井地質(zhì)力學(xué)參數(shù)空間建模方法

在利用地震屬性預(yù)測得到地下空間的三維地層縱波速度數(shù)據(jù)的基礎(chǔ)上，可運(yùn)用巖石力學(xué)方法進(jìn)一步預(yù)測各類鉆井地質(zhì)力學(xué)參數(shù)[21-23]。對于鄂爾多斯南部探區(qū)，通過實測數(shù)據(jù)統(tǒng)計擬合，建立起各類巖石物理及力學(xué)參數(shù)計算模型。利用縱波速度預(yù)測得到橫波速度、密度、孔隙度、泥質(zhì)含量等巖石物理參數(shù)，進(jìn)而預(yù)測彈性模量、泊松比、黏聚力、內(nèi)摩擦角、抗拉強(qiáng)度、上覆壓力、水平地應(yīng)力等力學(xué)參數(shù)。通過統(tǒng)計擬合建立了地層孔隙壓力模型為：

(6)

式中：pe為垂直應(yīng)力，MPa；pp為孔隙壓力，MPa；Vsh為泥質(zhì)含量；φ為孔隙度；σv為上覆壓力，MPa；w1、w2、w3、w4、w5為經(jīng)驗系數(shù)，對于鄂南工區(qū)，w1=4.71,w2=-4.086,w3=-0.046 3,w4=0.018 8,w5=-0.021 5 。

井壁穩(wěn)定分析是鉆井地質(zhì)力學(xué)分析的重點(diǎn)，通常根據(jù)剪切破壞準(zhǔn)則和拉伸破壞準(zhǔn)則建立井壁失穩(wěn)力學(xué)模型，最終可得到坍塌壓力和破裂壓力的計算公式：

(7)

pf=3σh-σH-αpp+St

(8)

式中：pc為坍塌壓力，MPa；pf為破裂壓力，MPa；σH和σh分別為最大和最小水平地應(yīng)力，MPa；C為巖石的黏聚力，MPa；α為有效應(yīng)力系數(shù)；K=ctan(45-φ/2) ，φ為內(nèi)摩擦角，°；η為應(yīng)力修正系數(shù)；St為抗拉強(qiáng)度，MPa。

巖石可鉆性是表征地層抗鉆能力的重要參數(shù)，其與層速度關(guān)系密切。根據(jù)實驗室測試數(shù)據(jù)分析，建立了適用于目標(biāo)工區(qū)的可鉆性級值計算模型：

(9)

式中：Kdy1和Kdpdc分別為牙輪鉆頭和PDC鉆頭可鉆性級值；k0、k1、k2、k3為由統(tǒng)計擬合得到的經(jīng)驗系數(shù)，對于鄂南工區(qū)，k0=0.000 36,k1=1.153 00,k2=0.000 257,k3=1.175 00 。

按照上述方法，可基于地震層速度導(dǎo)出各類鉆井地質(zhì)力學(xué)參數(shù)三維數(shù)據(jù)體，建立相應(yīng)的空間數(shù)據(jù)模型。運(yùn)用三維地質(zhì)體描述軟件，可對相關(guān)力學(xué)參數(shù)進(jìn)行三維空間成像，并進(jìn)行空間剖面和切片分析，以獲得各類力學(xué)參數(shù)的空間分布特征，為鉆井優(yōu)化設(shè)計提供依據(jù)。

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

鎮(zhèn)涇探區(qū)位于鄂爾多斯盆地西南緣，區(qū)內(nèi)鉆井以水平井為主，由淺至深依次鉆遇第四系、羅漢洞組、環(huán)河組、華池組、洛河-宜君組、安定組、直羅組、延安組及延長組，完鉆井深為3 000 m左右。前期6口完鉆井由于對地層可鉆性認(rèn)識不全面，導(dǎo)致鉆頭選型針對性差，破巖效果不理想，機(jī)械鉆速有待提高。

將鎮(zhèn)涇探區(qū)內(nèi)東部井區(qū)待鉆的YS101、YS301、YS5、YS6、YS9等5口井作為試驗井，利用上述鉆井地質(zhì)力學(xué)參數(shù)空間建模方法，分層段依次由淺至深建立了該地區(qū)的層速度、彈性模量、抗壓強(qiáng)度、可鉆性級值、地應(yīng)力、孔隙壓力、坍塌壓力、破裂壓力等多種參數(shù)的空間數(shù)據(jù)體。限于篇幅，只顯示目標(biāo)井區(qū)重點(diǎn)層段的可鉆性級值及孔隙壓力、坍塌壓力、破裂壓力(簡稱地層三壓力)空間數(shù)據(jù)體三維成像效果(圖1—4)。由圖1—4可知，力學(xué)參數(shù)建模結(jié)果具有較高的分辨率。從上述數(shù)據(jù)體中提取5口試驗井的數(shù)據(jù)剖面(以YS5井為例，圖5)，進(jìn)行試驗井地層力學(xué)參數(shù)特征分析，在此基礎(chǔ)上對其進(jìn)行必要的鉆井工藝措施優(yōu)化，在試驗井完鉆后收集了相關(guān)測井、鉆井?dāng)?shù)據(jù)以驗證建模預(yù)測效果。由圖5可知，預(yù)測的坍塌壓力、鉆井液密度與井徑擴(kuò)大率等曲線的匹配程度較高。將試驗井力學(xué)參數(shù)預(yù)測值與實測結(jié)果進(jìn)行了對比(表1)，各測試點(diǎn)的預(yù)測誤差均低于10%。以上均證明鉆井地質(zhì)力學(xué)參數(shù)建模技術(shù)的可靠性。

圖1鎮(zhèn)涇東部井區(qū)洛河-宜君組孔隙壓力數(shù)據(jù)體成像效果

圖3鎮(zhèn)涇東部井區(qū)環(huán)河組破裂壓力數(shù)據(jù)體成像效果

圖2鎮(zhèn)涇東部井區(qū)延安組坍塌壓力數(shù)據(jù)體成像效果

圖4鎮(zhèn)涇東部井區(qū)安定組可鉆性級值數(shù)據(jù)體成像效果

圖5從空間數(shù)據(jù)體中提取的YS5井的層速度與力學(xué)參數(shù)預(yù)測剖面

通過分析鎮(zhèn)涇東部井區(qū)地層三壓力空間建模結(jié)果，確定了該區(qū)三壓力縱向和橫向分布特征，在重點(diǎn)考慮環(huán)河組、延安組、延長組等壓力分布異常地層的基礎(chǔ)上，結(jié)合工程需求，提出了目標(biāo)井區(qū)水平井井身結(jié)構(gòu)優(yōu)化方案。該地區(qū)前期水平井采用三開井身結(jié)構(gòu)，針對提速提效需要，提出了2種井身結(jié)構(gòu)優(yōu)化方案。一是對于該區(qū)井深超過2 800 m的水平井，仍然推薦采用三開井身結(jié)構(gòu)，為提高斜井段鉆井速度及控制阻卡，建議二開采用復(fù)合井眼(表2)。對于井深低于2 800 m的水平井，充分考慮提高鉆速、降低成本、提高固井質(zhì)量的需要，推薦二開井身結(jié)構(gòu)(表3)?；诰斫Y(jié)構(gòu)優(yōu)化方案，進(jìn)一步優(yōu)化了目標(biāo)井區(qū)的鉆井液密度(表4)。

表1 鎮(zhèn)涇東部井區(qū)試驗井測試點(diǎn)上的實測值和預(yù)測值對比

表2 鎮(zhèn)涇東部井區(qū)水平井三開井身結(jié)構(gòu)優(yōu)化方案

表3 鎮(zhèn)涇東部井區(qū)水平井二開井身結(jié)構(gòu)優(yōu)化方案

表4 鎮(zhèn)涇東部井區(qū)水平井鉆井液密度優(yōu)化方案

通過分析鎮(zhèn)涇東部井區(qū)地層可鉆性空間建模結(jié)果，確定了該區(qū)巖石可鉆性縱向和橫向分布特征，結(jié)合實際施工條件，對該區(qū)的鉆頭型號及相關(guān)工程參數(shù)進(jìn)行了優(yōu)化(表5)。

表5 鎮(zhèn)涇東部井區(qū)水平井鉆頭及相關(guān)工程參數(shù)優(yōu)化方案

5口試驗井均采用了上述工程優(yōu)化方案。實鉆表明，5口井平均機(jī)械鉆速為10.80 m/h，對比該井區(qū)前期完鉆井7.49 m/h的平均鉆速，提高了44.19%，提速效果明顯。其中，YS9井采用推薦的二開井身結(jié)構(gòu)，鉆井周期為23.67 d，較之前期完鉆井縮短了51.70%，節(jié)約成本近210×104元。試驗井的平均井徑擴(kuò)大率和漏失量相比前期完鉆井均減少40%以上。上述現(xiàn)場實際情況均表明，文中提出的鉆井地質(zhì)力學(xué)參數(shù)建模方法可有效滿足工程優(yōu)化的需求。

5 結(jié) 論

(1) 地層聲波傳播速度是定量計算各類巖石物理與力學(xué)參數(shù)的基礎(chǔ)數(shù)據(jù)，基于地震資料預(yù)測的速度信息可用于進(jìn)行鉆井地質(zhì)力學(xué)參數(shù)三維空間建模。

(2) 地震屬性和聲波速度之間存在映射關(guān)系，據(jù)此可運(yùn)用智能算法由地震信息反演得到層速度，進(jìn)而建立鉆井地質(zhì)力學(xué)參數(shù)的空間三維數(shù)據(jù)模型。

(3) 現(xiàn)場應(yīng)用情況表明，文中提出的鉆井地質(zhì)力學(xué)參數(shù)空間建模方法能達(dá)到較高的參數(shù)預(yù)測精度和成像分辨率，以此為基礎(chǔ)進(jìn)行的鉆井工藝措施優(yōu)化具有較強(qiáng)的針對性和適用性。

[1] 吳超，陳小鋒，張東清，等.巖石力學(xué)特征分析技術(shù)及其在元壩氣田的應(yīng)用[J].測井技術(shù)，2015，39(2)：247-252.

[2] 吳勝和.儲層表征與建模[M].北京：石油工業(yè)出版社，2010：290-355.

[3] 樂友喜，曾勉，問雪，等.利用序貫高斯隨機(jī)模擬分析構(gòu)造圖的不確定性[J].石油地球物理勘探，2017，52(2)：333-339.

[4] 吳超，臧艷彬，張東清，等.定向井隨鉆井壁穩(wěn)定預(yù)測方法及其應(yīng)用[J].石油學(xué)報，2015，36(10)：1290-1298.

[5] WU C，CHEN M，JIN Y.A prediction method of borehole stability based on seismic attribute technology[J]. Journal of Petroleum Science and Engineering，2009，65(3/4)：208-216.

[6] WU C，LIU J，ZHANG D，et al. A prediction of borehole stability while drilling preliminary prospecting wells based on seismic impedance[J].Petroleum Exploration and Development.2015，42(1)：390-395.

[7] 金衍，陳勉.利用地震記錄鉆前預(yù)測井壁穩(wěn)定性研究[J].石油學(xué)報，2004，25(1)：89-92.

[8] 金衍，陳勉，楊小奇.利用層速度鉆前預(yù)測安全泥漿密度窗口[J].巖石力學(xué)與工程學(xué)報，2004，23(14)： 2430-2433.

[9] WU C，CHEN M，JIN Y.Real-time prediction method of borehole stability[J].Petroleum Exploration and Development，2008，35(1)：80-84.

[10] 金衍，陳勉，張旭東.鉆前井壁穩(wěn)定預(yù)測方法的研究[J].石油學(xué)報，2001，22(3)：96-99.

[11] AZEDOVE L，NUNES R，CORREIA P，et al. Multidimensional scaling for the evaluation of a geostatistical seismic elastic inversion methodology[J].Geophysics，2014，79(1)：M1-M10.

[12] 吳超，劉建華，張東清，等.基于非線性反演的井壁穩(wěn)定隨鉆預(yù)測方法[J].巖土工程學(xué)報，2015，37(10)： 1810-1817.

[13] BOSCH M，MUKERJI T，GONZALEZ E F.Seismic inversion for reservoir properties combining statistical rock physics and geostatistics：A review[J].Geophysics，2010，75(5)：A165-A176.

[14] ESMERSOY C，RAMIREZ A，TEEBENNY S.A new fully integrated method for seismic geohazard prediction ahead of the bit while drilling[J].Leading Edge，2013，32(10)：1222-1233.

[15] 盧寶坤，史謌.測井資料與地震屬性關(guān)系研究綜述[J].北京大學(xué)學(xué)報(自然科學(xué)版)，2005，41(1)：154-160.

[16] 張陽，邱隆偉，李際，等.基于模糊C均值地震屬性聚類的沉積相分析[J].中國石油大學(xué)學(xué)報(自然科學(xué)版)， 2015，39(4)：53-61.

[17] NEVES F A，ZAHRANI M S，BREMKAMP S W，et al. Detection of potential fractures and small faults using seismic attributes[J]. Leading Edge，2004，23(9)：903-906.

[18] NATH S K，DEWANGAN P.Detection of seismic reflections from seismic attributes through fractal analysis[J]. Geophysical Prospecting，2002，50(3)：341-360.

[19] 彭琦，張茹，謝和平，等.基于AE時間序列的巖爆預(yù)測模型[J].巖土力學(xué)，2009，30(5)：1436-1440.

[20] 曹愛武，顧圣平，何露，等.基于小波神經(jīng)網(wǎng)絡(luò)的水庫泥沙預(yù)測[J].江南大學(xué)學(xué)報(自然科學(xué)版)，2015，14(3)： 338-343.

[21] ZHANG J.Borehole stability analysis accounting for anisotropies in drilling to weak bedding planes[J].International Journal of Rock Mechanics and Mining Sciences，2013，60(3)：160-170.

[22] 盧運(yùn)虎，陳勉，袁建波，等.各向異性地層中斜井井壁失穩(wěn)機(jī)理[J].石油學(xué)報，2013，34(3)：563-568.

[23] 樓一珊，金業(yè)權(quán).巖石力學(xué)與石油工程[M].北京：石油工業(yè)出版社，2006：99-137.