謝丹艷, 張宏兵, 孫樹(shù)林

(河海大學(xué)地球科學(xué)與工程學(xué)院, 江蘇 南京 210098)

0 引 言

水平井工藝在中國(guó)有了一定發(fā)展[1-4],但是缺乏水平井相應(yīng)的固井質(zhì)量評(píng)價(jià)方法。水平井固井中存在有別于垂直井的多種問(wèn)題,周燦燦等[5]在水平井綜合解釋中提出的井壁裂縫、泥漿侵入、儀器偏心等因素對(duì)水平井聲波測(cè)井結(jié)果的影響。針對(duì)水平井可能存在水泥環(huán)上層弓形缺失、水平井弓形缺失部分存在混漿等多種復(fù)雜情況,傳統(tǒng)的解析方法無(wú)法實(shí)現(xiàn)水平井井孔聲場(chǎng)模擬,需要采用通用性能好、幾乎適合分析所有聲場(chǎng)的三維有限差分方法進(jìn)行井孔聲場(chǎng)的數(shù)值模擬。劉繼生等[6]、宋若龍等[7]利用三維有限差分方法模擬水泥環(huán)扇區(qū)缺失的非軸對(duì)稱井孔聲場(chǎng)。林偉軍等[8]、陳德華等[9]利用同樣方法分別計(jì)算了傾斜、裂隙等復(fù)雜地層的井孔聲場(chǎng);叢健生等[10]同樣利用有限差分方法模擬了了水平井地層界面聲波測(cè)井的響應(yīng);李瑞豐等[11]利用三維有限差分方法在水平井段存在套管偏心和巖屑床條件下正演模擬了井中聲場(chǎng)分布。采用三維有限差分方法模擬水平井的復(fù)雜井孔聲場(chǎng)可行。

在三維有限差分模擬中,針對(duì)計(jì)算空間有限引入的人為邊界所產(chǎn)生的虛假反射,Wang等[12]提出了非裂化的完全匹配方法(NPML)弱化這種邊界反射;Song等[13]將NPML應(yīng)用于模擬孔隙介質(zhì)中彈性波傳播問(wèn)題,認(rèn)為NPML吸收邊界,不僅從效果上等同于SPML,而且在算法實(shí)現(xiàn)上更簡(jiǎn)易,且不占更多的計(jì)算內(nèi)存。在水平井井孔中存在多種復(fù)雜的情況,本文對(duì)其中可能存在的水泥環(huán)上層弓形缺失的多種情況進(jìn)行分析。

1 水平井井孔聲場(chǎng)模擬方法

1.1 水平井井孔各種模型

圖1是水泥環(huán)的上層存在一個(gè)弓形的水泥缺失。o為聲源,井孔中介質(zhì)從里到外標(biāo)號(hào)為1、2、3、4依次代表孔內(nèi)流體、鋼套管、水泥環(huán)、地層等,h代表弓高,不同弓高會(huì)產(chǎn)生2種不同的缺失形狀,其中h=30 mm(即水泥環(huán)厚度)為分界點(diǎn)。圖1(a)為h<30 mm,僅第Ⅱ界面部分區(qū)域耦合不好;圖1(b)為h≥30 mm,在第Ⅰ、第 Ⅱ界面均存在耦合不好;圖1(c)為井壁四周接收器布置方位。

(a)第Ⅱ界面耦合不好的弓形缺失 (b)第Ⅰ、第Ⅱ界面均耦合不好的弓形缺失 (c)井壁四周接收器布置方位

1.2 高階交錯(cuò)網(wǎng)格的三維有限差分方法

采用速度和應(yīng)力建立井孔聲場(chǎng)的一階偏微分方程,使用三維有限差分的方法分解聲場(chǎng)的一階偏微分方程。采用交錯(cuò)網(wǎng)格的差分格式提高運(yùn)算效率,函數(shù)一階微分的8階交錯(cuò)網(wǎng)格差分格式近似為

(1)

由式(1)可知,函數(shù)微分的采樣點(diǎn)位置和函數(shù)的采用點(diǎn)位置要相差半個(gè)步長(zhǎng)。如果時(shí)間上采用二階差分,空間上采用8階差分,則a0=1225/1025,a1=-245/3072,a2=49/5120,a3=-5/7168?？梢垣@得速度和應(yīng)力9個(gè)變量的差分公式[14]

(2)

只需要將這9個(gè)微分方程中的微分算子換成相應(yīng)的有限差分格式,就能得到相應(yīng)的有限差分方程。有限差分模擬結(jié)果需要滿足收斂性和穩(wěn)定性[13],穩(wěn)定性的時(shí)間步長(zhǎng)為

(3)

確保收斂性的空間間隔又必須滿足

(4)

式中,vmax和vmin分別指介質(zhì)的最大和最小聲速;fmax是指聲源的最高有效頻率。

1.3 參數(shù)的選取

聲源采用高斯包絡(luò)調(diào)制的正弦波,中心頻率為18 kHz,帶寬6 kHz。選用套管半徑為0.062 85 mm,套管壁厚度0.007 mm。實(shí)際計(jì)算的區(qū)域?yàn)閤、y軸正負(fù)方向各0.12 m,原點(diǎn)位于中心;z軸負(fù)向取0.25 m,正向取1.625 m;收斂邊界取0.07 m。根據(jù)三維有限差分的收斂穩(wěn)定條件,取時(shí)間步長(zhǎng)為1.5×10-7s,取x、y方向的步長(zhǎng)為0.002 m,z方向的步長(zhǎng)為0.025 m,全空間劃分為189×189×84個(gè)網(wǎng)格。

2 三維有限差分方式的驗(yàn)證

為了驗(yàn)證三維有限差分法的可靠性,把有限差分的計(jì)算結(jié)果與實(shí)軸積分法的結(jié)果進(jìn)行對(duì)比,膠結(jié)完好的套管井聲全波波形對(duì)比結(jié)果如圖2所示。

圖2 膠結(jié)完好的套管井模型的三維有限差分法和實(shí)軸積分法計(jì)算結(jié)果對(duì)比圖

圖2中黑色曲線為實(shí)軸積分法的計(jì)算結(jié)果,紅色曲線為三維有限差分法的計(jì)算結(jié)果,對(duì)比2種方法的計(jì)算結(jié)果可以發(fā)現(xiàn),2種方法計(jì)算得到的裸眼井聲全波曲線,首波都是地層縱波,其次是地層橫波和斯通利波,地層縱、橫波的波幅和相位幾乎完全重合,斯通利波也基本吻合,說(shuō)明三維有限差分方法的計(jì)算結(jié)果可靠。

3 水平井弓形缺失聲場(chǎng)波形分析

對(duì)水平井可能存在的水泥環(huán)上層弓形缺失、弓形缺失部分存在混漿的情況進(jìn)行數(shù)值模擬,介質(zhì)參數(shù)詳見(jiàn)表1。

表1 模型參數(shù)

3.1 水泥環(huán)上層弓形缺失

在水泥灌漿過(guò)程中,由于重力作用,在水泥環(huán)固結(jié)后可能會(huì)在上層形成弓形的缺失。為研究缺失高度對(duì)聲場(chǎng)的影響,分別模擬弓高為5、10、15、20、25、30、40、50、60 mm等9種情況,提取源距為1.0 m處的波形(見(jiàn)圖3)。

圖3 不同弓高的聲場(chǎng)波形比較圖(源距1.0 m)

圖4 實(shí)軸積分計(jì)算的不同竄槽類型波形圖(源距1.0 m)

圖3是源距1.0 m處接收到的不同缺失高度的聲場(chǎng)波形。提取首波進(jìn)行研究,從圖3可知,缺失高度在5～25 mm的4組首波(即實(shí)線)旅行時(shí)間要比缺失高度在30 mm以上的(即虛線)來(lái)得晚,其中缺失高度在5～25 mm的首波旅行時(shí)間為0.385 ms左右,缺失高度大于30 mm的首波旅行時(shí)間約為0.245 ms。圖4為實(shí)軸積分法計(jì)算出的垂直井第Ⅰ界面、第Ⅱ界面竄槽水層厚度為10 mm的模型和空套管模型在1.0 m處的聲場(chǎng),其中第Ⅰ界面竄槽的首波旅行時(shí)間約為0.225 ms,空套管的首波旅行時(shí)間約為0.19 ms,兩者都比較接近第Ⅰ界面套管波的旅行時(shí)間;而第Ⅱ界面竄槽時(shí)首波即第Ⅱ界面特征波的旅行時(shí)間約為0.33 ms。則初步認(rèn)為圖3中的0.245 ms為第Ⅰ界面套管波旅行時(shí)間,0.385 ms為第Ⅱ界面特征波旅行時(shí)間。之所以比實(shí)軸積分計(jì)算的時(shí)間都略大,可能是弓形缺失為局部缺失,其信號(hào)強(qiáng)度較周向缺失小,可觀察的首波信號(hào)出現(xiàn)得較晚。結(jié)合圖1(a)的模型,當(dāng)缺失在30 mm之內(nèi)在時(shí)第Ⅱ界面部分區(qū)域存在弓形水層,耦合不好,有可能產(chǎn)生第Ⅱ界面特征波;同樣,結(jié)合圖1(b)的模型,當(dāng)缺失高度超過(guò)30 mm,相當(dāng)于在井孔部分區(qū)域形成了空套管,而空套管的首波即為第Ⅰ界面套管波。另外,在弓形缺失模型中,只有局部區(qū)域產(chǎn)生竄槽,其他區(qū)域則膠結(jié)完好,由于膠結(jié)完好的首波為地層波,其速度較套管波和特征波都小,所以弓形缺失首波以體現(xiàn)套管波和特征波為主。最終可以認(rèn)為缺失高度在30～60 mm的波形其首波為套管波,缺失高度在5～25 mm首波為特征波,拐點(diǎn)30 mm正是水泥環(huán)的厚度。

同時(shí),缺失高度在5～25 mm的4組首波幅度,隨著缺失高度增加,首波幅度隨之增大;同樣,缺失高度在30 mm以上的5組首波幅度也是隨著缺失高度增大而增大。

3.2 混漿

在水平井固井過(guò)程中,有可能出現(xiàn)的另一種情況是由于重力作用,當(dāng)水平井水泥固結(jié)后,在上層形成弓形的水泥和空氣或者水的混合物,稱之為混漿。選擇3種不同的混漿(Mixing2、Mixing3和Mixing4)進(jìn)行數(shù)值模擬(參數(shù)見(jiàn)表1),同時(shí)與弓形區(qū)域?yàn)樗那闆r(Mixing1)進(jìn)行對(duì)比,提取源距為1.0 m處的水平井混漿聲場(chǎng)首波的比較圖(見(jiàn)圖5)。

圖5 水平井上層混漿聲場(chǎng)首波比較圖(源距1.0 m)

圖6 首波幅度和不同混漿物關(guān)系圖(源距1.0 m)

圖5為源距1.0 m位置接收到的水平井上層水泥弓形缺失聲場(chǎng)首波圖,而缺失部分為水或者是泥水混合物等多種情況。由于選取的弓高為40 mm,根據(jù)3.1的結(jié)論認(rèn)為首波是套管波,同時(shí)計(jì)算得到1.0 m處的首波旅行時(shí)間約為0.215 ms,符合套管波旅行時(shí)間。圖6為首波幅度和不同混漿物關(guān)系圖,幅度進(jìn)行了同比放大。由圖5、圖6可知,不同混漿參數(shù)對(duì)聲場(chǎng)首波的幅度和旅行時(shí)間影響都不是很明顯,但是與內(nèi)部物質(zhì)為水的情況相比,首波的幅度明顯減小,這個(gè)可以解釋為套管外部的物質(zhì)密度小,套管外壁的反射就越強(qiáng)烈,使得套管波的幅度就越大。

4 弓形缺失方位確定

對(duì)于弓形缺失的水平井聲場(chǎng),可以在井壁四周布置一系列接收器,通過(guò)不同方位接收到的聲場(chǎng)判斷缺失的具體方位。井壁四周的接收器位置如圖1(c)所示,聲源依舊在井軸位置,在井壁四周布置8個(gè)接收器,分別標(biāo)注為1～8號(hào)接收器,接收器之間間隔為45 °。對(duì)水平井上層弓形缺失高度為40 mm的情況進(jìn)行數(shù)值模擬。在垂直源距為1.0 m處接收信號(hào),采集每個(gè)接收器首波的幅度和旅行時(shí)間,得到四周接收的首波幅度和旅行時(shí)間的比較圖(見(jiàn)圖7)。

圖7(a)為源距1.0 m處井壁四周8個(gè)接收器接收到的首波幅度比較圖(幅度進(jìn)行了同比放大)。從圖7(a)可知,0 °、45 °、90 °和135 °等4個(gè)接收器接收到的首波波谷的幅度明顯比其余4個(gè)位置的幅度小,也就是靠近缺失位置接收到的首波幅度明顯減小,而遠(yuǎn)離缺失位置的接收器接收到的首波幅度就較大;圖7(b)為源距1.0 m處井壁四周8個(gè)接收器接收到的首波旅行時(shí)間比較圖,和首波幅度的變化規(guī)律相似,0 °、45 °、90 °和135 °等4個(gè)接收器接收到的首波旅行時(shí)間明顯比其余4個(gè)位置的時(shí)間小,也就是靠近缺失位置接收到的首波旅行時(shí)間明顯減小,而遠(yuǎn)離缺失位置的接收器接收到的首波旅行時(shí)間就較大。

在井壁四周布置接收器,通過(guò)分析不同方位首波的幅度和旅行時(shí)間的變化,可以比較準(zhǔn)確地判斷出弓形缺失的缺失方位。

圖7 井壁四周接收的首波變化圖

5 結(jié) 論

(1) 水平井弓形缺失的弓高對(duì)于聲場(chǎng)首波的幅度和旅行時(shí)間影響很明顯,其中水泥環(huán)的厚度30 mm是一個(gè)重要拐點(diǎn),弓高低于30 mm首波為第Ⅱ界面特征波,高于30 mm首波為套管波。

(2) 弓形缺失部分的物質(zhì)不同也影響著聲場(chǎng),不同混漿物對(duì)于聲場(chǎng)首波的幅度和旅行時(shí)間影響不是特別明顯,但是與內(nèi)部為水的情況相比,首波幅度就明顯減小。通過(guò)在井壁四周布置接收器分析首波幅度和旅行時(shí)間的變化,可以確定弓形缺失的方位。

(3) 在水平井井孔聲場(chǎng)模擬中,水泥環(huán)上層弓形缺失的影響不可忽略,需要對(duì)可能出現(xiàn)的情況作出有效判斷,才能正確評(píng)估水平井的固井質(zhì)量。

參考文獻(xiàn):

[1] 張國(guó)仿. 提高水平井固井質(zhì)量工藝技術(shù) [J]. 江漢石油職工大學(xué)學(xué)報(bào), 2002, 15(3): 20-22.

[2] 孫莉, 黃曉川, 向興華. 國(guó)內(nèi)水平井固井技術(shù)及發(fā)展 [J]. 鉆采工藝, 2005, 28(5): 23-30.

[3] 馮京海, 郝新朝, 白亮清, 等. 冀東油田水平井固井技術(shù) [J]. 石油鉆采工藝, 2007, 29(S1): 32-39.

[4] 胡黎明, 趙留陽(yáng), 郭振斌. 提高水平井固井質(zhì)量的措施 [J]. 中國(guó)石油和化工標(biāo)準(zhǔn)與質(zhì)量, 2012, 32(4): 217-219.

[5] 周燦燦, 王昌學(xué). 水平井測(cè)井解釋技術(shù)綜述 [J]. 地球物理學(xué)進(jìn)展, 2006, 21(1): 152-160.

[6] 劉繼生, 馬俊, 王克協(xié), 等. 非軸對(duì)稱井孔聲場(chǎng)三維應(yīng)力—速度有限差分?jǐn)?shù)值模擬 [C]∥2001年中國(guó)地球物理學(xué)會(huì)年刊——中國(guó)地球物理學(xué)會(huì)第十七屆年會(huì)論文集, 北京, 2001.

[7] 宋若龍. 非軸對(duì)稱套管井聲場(chǎng)并行計(jì)算及聲波固井質(zhì)量評(píng)估理論與方法研究 [D]. 長(zhǎng)春: 吉林大學(xué), 2008.

[8] 林偉軍, 王秀明, 張海瀾. 傾斜地層中的井孔聲場(chǎng) [J]. 地球物理學(xué)報(bào), 2006, 49(1): 284-294.

[9] 陳德華, 叢建生, 徐德龍, 等. 裂縫性地層中的井孔聲場(chǎng)模擬 [J]. 大慶石油學(xué)院學(xué)報(bào), 2004, 28(3): 4-6.

[10] 叢健生, 喬文孝. 水平井地層界面聲波測(cè)井響應(yīng)模擬分析 [J]. 測(cè)井技術(shù), 2008, 32(1): 29-32.

[11] 李瑞豐, 馮逾, 李東旭. 水平井兩種復(fù)雜井周狀態(tài)下聲波測(cè)井響應(yīng)的模擬分析 [J]. 石油儀器, 2013, 27(4): 22-25.

[12] Wang Tsili, Tang Xiaoming. Finite-difference Modeling of Elastic Wave Propagation: Anonsplitting Perfectly Matched Layer Approach [J]. Geophysics, 2004, 68: 1749-1755.

[13] Song Ruolong, Ma Jun, Wang Kexie. The Application of the Nonsplitting Perfectly Matched Layer in Numerical Modeling of Wave Propagation Poroelastic Media [J]. Applied Geophysics, 2005, 2(4): 216-222.

[14] 張海瀾, 王秀明, 張碧星. 井孔的聲場(chǎng)和波 [M]. 北京: 科學(xué)出版社, 2003: 148-154.