唐軍，章成廣，張碧星，師芳芳（1.油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室（長(zhǎng)江大學(xué)）；2.長(zhǎng)江大學(xué)地球物理與石油資源學(xué)院；.中國(guó)科學(xué)院聲學(xué)研究所）

?

基于聲波-變密度測(cè)井的固井質(zhì)量評(píng)價(jià)方法

唐軍1， 2，章成廣1， 2，張碧星3，師芳芳3
（1.油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室（長(zhǎng)江大學(xué)）；2.長(zhǎng)江大學(xué)地球物理與石油資源學(xué)院；3.中國(guó)科學(xué)院聲學(xué)研究所）

摘要：采用數(shù)值模擬與刻度井群測(cè)井相結(jié)合的方法進(jìn)行套管井聲場(chǎng)研究，分析套管井竄槽角度、Ⅱ界面（水泥與地層膠結(jié)界面）膠結(jié)質(zhì)量對(duì)聲波-變密度測(cè)井的影響，在此基礎(chǔ)上對(duì)利用CBL/VDL測(cè)井資料評(píng)價(jià)固井質(zhì)量的方法進(jìn)行改進(jìn)，提出新的評(píng)價(jià)標(biāo)準(zhǔn)。分別采用實(shí)軸積分法和2.5維有限差分算法對(duì)軸對(duì)稱聲場(chǎng)與非軸對(duì)稱聲場(chǎng)進(jìn)行數(shù)值模擬，并與刻度井測(cè)井?dāng)?shù)據(jù)對(duì)比后發(fā)現(xiàn)：數(shù)值模擬方法可靠；隨著竄槽角度減小，套管首波幅度減小，且水泥漿密度越大，套管波衰減越快；聲波測(cè)井能辨別的最小竄槽角度為45°；通過波形圖在時(shí)域內(nèi)很難提取地層波信息，可利用主頻識(shí)別地層波，且在灰?guī)r地層中當(dāng)Ⅱ界面環(huán)隙寬度較小時(shí)地層波才可作為首波接收。依據(jù)套管井聲場(chǎng)研究結(jié)果，建立了新的基于聲波-變密度測(cè)井的固井質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)。圖8表3參20

關(guān)鍵詞：固井質(zhì)量；聲波-變密度測(cè)井；套管井聲場(chǎng)；竄槽角度；界面膠結(jié)指數(shù)

0 引言

為保證油氣井的安全生產(chǎn)及延長(zhǎng)套管使用壽命，套管與地層間的空隙一般用水泥充填，套管、水泥以及地層三者要緊密膠結(jié)。由于井深一般大于1 000 m，固井質(zhì)量的有效快速檢測(cè)顯得尤為重要。聲波測(cè)井是固井質(zhì)量檢測(cè)的主要手段，一般包括Ⅰ界面（套管與水泥膠結(jié)界面）和Ⅱ界面（水泥與地層膠結(jié)界面）膠結(jié)質(zhì)量評(píng)價(jià)，研究認(rèn)為套管波信息可用于Ⅰ界面評(píng)價(jià)，地層波信息可用于Ⅱ界面評(píng)價(jià)［1-2］。

井孔聲場(chǎng)研究是運(yùn)用聲波測(cè)井檢測(cè)固井質(zhì)量的基礎(chǔ)。國(guó)內(nèi)外研究者開展了相關(guān)研究［3-11］，但均是基于嚴(yán)格的波動(dòng)方程的理論研究，而在實(shí)際套管井中，地層、井內(nèi)流體的非均質(zhì)性和聲波檢測(cè)儀器等都會(huì)影響井內(nèi)聲場(chǎng)。雖然一些學(xué)者利用數(shù)值模擬技術(shù)研究了測(cè)井儀器、套管類型、裂縫、地層傾角等因素對(duì)套管井聲場(chǎng)的影響［12-17］，但受實(shí)驗(yàn)條件所限，研究采用的物理模型往往只有幾米，遠(yuǎn)小于實(shí)際套管井的尺寸，因此研究結(jié)果存在局限性。此外，聲波測(cè)井儀器種類多，其中聲波-變密度（CBL/VDL）測(cè)井儀器因其檢測(cè)成本低廉、數(shù)據(jù)采集快捷，在國(guó)內(nèi)外油田中被大量使用，但是關(guān)于竄槽角度和界面膠結(jié)質(zhì)量評(píng)價(jià)的研究不多。

本文根據(jù)實(shí)際建造的刻度井群，采用2.5維有限差分法計(jì)算非軸對(duì)稱聲場(chǎng)，并與刻度井中利用聲波-變密度測(cè)井儀器檢測(cè)的套管首波進(jìn)行比對(duì)，驗(yàn)證數(shù)值模擬方法的可靠性，分析Ⅰ界面竄槽角度對(duì)聲波-變密度測(cè)井的影響；采用實(shí)軸積分法計(jì)算軸對(duì)稱聲場(chǎng)，分析Ⅱ界面環(huán)隙寬度對(duì)套管波、地層波的影響。在分析水泥漿密度、竄槽角度等對(duì)聲波-變密度測(cè)井影響的基礎(chǔ)上，對(duì)利用CBL/VDL測(cè)井資料評(píng)價(jià)固井質(zhì)量的方法進(jìn)行改進(jìn)。

1 套管井聲場(chǎng)研究方法

本文采用實(shí)際刻度井群聲波測(cè)量與數(shù)值模擬結(jié)合的方法開展套管井聲場(chǎng)研究。一方面，利用刻度井群聲波測(cè)井結(jié)果驗(yàn)證數(shù)值模擬的準(zhǔn)確性、可靠性；另一方面，通過數(shù)值模擬降低研究成本，彌補(bǔ)刻度井群模型數(shù)量上的不足。

1.1 刻度井群概況

新疆某測(cè)井公司施工建造了功能完善的水泥膠結(jié)刻度井群，主要用于下井儀器的刻度和標(biāo)定。該井群共由8口井組成，標(biāo)記為1—8號(hào)，水泥漿密度在1.20～2.25 g/cm3，地層巖性包括泥巖、砂巖、石灰?guī)r、花崗巖4種。其中，3號(hào)、4號(hào)兩口井模擬的是扇形竄槽的情況，套管直徑均為139.70 mm、壁厚7.72 mm，水泥環(huán)厚度38 mm，地層內(nèi)徑216 mm，外徑1 216 mm；從下到上分別模擬了套管偏心和Ⅰ界面22.5°、45°、90°、180°、360°竄槽的情況（見圖1）；3號(hào)井固井水泥漿密度為1.89 g/cm3，4號(hào)井固井水泥漿密度為1.20 g/cm3。由圖1可知，隨著竄槽角度變小，套管首波幅度逐漸變小，竄槽角度與套管波幅度變化趨勢(shì)一致，說明聲波測(cè)井可用于竄槽角度的評(píng)價(jià)。

圖1　3號(hào)刻度井模型和CBL-VDL測(cè)井波形

1.2 套管井聲場(chǎng)數(shù)值模擬

水泥在套管與地層間的固結(jié)情況可以分為2類，一類是環(huán)形未膠結(jié)，即套管與水泥或者水泥與地層間存在環(huán)狀空隙，另一類是水泥在套管與地層間部分膠結(jié)。前者為軸對(duì)稱聲場(chǎng)，后者為非軸對(duì)稱聲場(chǎng)。

1.2.1 軸對(duì)稱套管井聲場(chǎng)

本文采用實(shí)軸積分法對(duì)軸對(duì)稱聲場(chǎng)進(jìn)行數(shù)值模擬。套管、水泥、地層均沿井軸對(duì)稱，可視為徑向可分的多層介質(zhì)圓柱形開放聲波導(dǎo)，有解析解，井內(nèi)聲場(chǎng)為：

1.2.2 非軸對(duì)稱套管井聲場(chǎng)

圖2為非軸對(duì)稱套管井模型，對(duì)應(yīng)油田實(shí)際生產(chǎn)情況為膠結(jié)水泥在某層段出現(xiàn)缺失，并在套管和地層間形成扇型空隙，被稱為竄槽。為了節(jié)省計(jì)算機(jī)存儲(chǔ)空間及提高計(jì)算速度，本文采用2.5維有限差分算法對(duì)非軸對(duì)稱聲場(chǎng)進(jìn)行數(shù)值模擬。

圖2　非軸對(duì)稱套管井模型

在井孔聲場(chǎng)問題中，井孔和地層參數(shù)只在二維平面（x-y平面）上變化，沿井軸（z軸）保持不變，這就允許將三維差分方程沿z軸作空間傅立葉變換，從而只需在二維平面網(wǎng)格上進(jìn)行一系列二維有限差分計(jì)算（對(duì)應(yīng)于一系列不同的z軸方向波數(shù)分量kz），再進(jìn)行反變換就能得到三維的結(jié)果。

對(duì)于一個(gè)非均勻的各向同性彈性介質(zhì)，在笛卡爾坐標(biāo)下，可定義其密度為ρ（x， y， z），一階和二階拉梅常數(shù)分別為λ（x， y， z）和μ（x， y， z）。該類問題屬于二階雙曲型偏微分方程求解問題，可以分解為關(guān)于速度和應(yīng)力的一階偏微分方程組，而后對(duì)該方程組沿z軸作空間傅立葉變換，可得：

引入單極子聲源：

針對(duì)單極子聲源，在x-y平面上對(duì)不同的kz進(jìn)行一系列二維有限差分計(jì)算，得到一系列二維聲場(chǎng)，再將其作傅立葉反變換，就能得到帶竄槽套管井的三維聲場(chǎng)。為保證有限差分的穩(wěn)定性，時(shí)間步長(zhǎng)的選擇需滿足：

為保證計(jì)算過程不發(fā)散，空間步長(zhǎng)的選擇需滿足：

為保證計(jì)算精確，最大的軸向波數(shù)需滿足：

2 套管井聲場(chǎng)研究結(jié)果及討論

2.1 數(shù)值模擬與刻度井測(cè)量數(shù)據(jù)對(duì)比

圖3為0°、90°、180°及360°竄槽情況下x-y平面上速度譜，可以看出，x-y平面上的速度譜能夠反映實(shí)際竄槽情況。

以刻度井群的實(shí)際尺寸及水泥、地層聲學(xué)特征為依據(jù)進(jìn)行數(shù)值模擬（參數(shù)見表1）。圖4為石灰?guī)r地層完全膠結(jié)情況下的2.5維有限差分?jǐn)?shù)模模擬與CBL/VDL測(cè)井的聲波全波列波形對(duì)比圖，可以看出，數(shù)值模擬與井中實(shí)測(cè)套管波首波到達(dá)時(shí)間一致，前3個(gè)波峰吻合較好，說明數(shù)值模擬結(jié)果可靠。

表1　數(shù)值模擬參數(shù)

2.2 扇形竄槽對(duì)聲波的影響

采用數(shù)值模擬方法獲得了22.5°、45°、90°、180°及360°竄槽的全波列波形。圖5為竄槽角度與套管波幅度的關(guān)系，其中歸一化套管波幅度為當(dāng)前套管波幅度與自由套管波幅度的比值。由圖5可知，水泥漿密度越大，隨著竄槽角度的減小套管波幅度衰減越快；竄槽角度在90°～270°時(shí)，不同水泥漿密度下的套管波幅度衰減程度差異明顯；不論水泥漿密度如何變化，竄槽角度小于45°時(shí)，套管波幅度變化很小，即聲波測(cè)井能辨別的最小竄槽角度為45°；3號(hào)刻度井（水泥漿密度1.89 g/cm3）的實(shí)測(cè)數(shù)據(jù)與數(shù)值模擬結(jié)果基本一致，進(jìn)一步驗(yàn)證了數(shù)值模擬的準(zhǔn)確性。

圖3　各種竄槽角度下的二維速度譜

圖4　數(shù)值模擬與CBL/VDL測(cè)井波形對(duì)比

圖5　不同水泥漿密度下竄槽角度與套管波幅度關(guān)系數(shù)值模擬結(jié)果

2.3 Ⅱ界面膠結(jié)質(zhì)量對(duì)聲波的影響

Ⅱ界面固井質(zhì)量受水泥環(huán)、地層等多因素的影響，目前還沒有統(tǒng)一的聲波檢測(cè)方法及標(biāo)準(zhǔn)。本文利用實(shí)軸積分法對(duì)砂巖、石灰?guī)r地層Ⅱ界面聲場(chǎng)進(jìn)行數(shù)值模擬。根據(jù)胡文祥等［8］和劉繼生等［18］的研究，石灰?guī)r地層波速度比套管波速度高，當(dāng)套管、水泥、地層三者膠結(jié)較好且源距（聲波發(fā)射探頭與接收探頭的距離）足夠大時(shí)，地層波先于套管波到達(dá)接收器，有利于開展地層波信息與Ⅱ界面膠結(jié)質(zhì)量的理論分析。

圖6、圖7分別為石灰?guī)r地層、源距1 524 mm （5 ft）、水泥漿密度1.5 g/cm3條件下Ⅱ界面環(huán)隙寬度變化時(shí)的波形和頻譜。

圖6　不同Ⅱ界面環(huán)隙寬度下全波列波形

圖7　不同Ⅱ界面環(huán)隙寬度下全波列頻譜

由圖6可知，Ⅱ界面環(huán)隙寬度為0時(shí)地層波初至先于自由套管時(shí)的套管波初至到達(dá)，Ⅱ界面環(huán)隙寬度在0.5～20.0 mm時(shí)地層初至波幅度減小，而地層后續(xù)波又與其他波疊加，在時(shí)域內(nèi)很難提取地層波信息。圖7則顯示利用主頻可以很好地區(qū)分套管波和地層波：當(dāng)環(huán)隙寬度為零時(shí)，幾乎沒有套管波，整個(gè)波列信息以地層波為主，此時(shí)波列主頻即為地層波主頻（15～17 kHz）；當(dāng)環(huán)隙寬度逐漸增大但小于5 mm時(shí)，主頻在15～17 kHz的地層波雖然存在，但會(huì)出現(xiàn)頻率大于17 kHz的波列成分；當(dāng)環(huán)隙寬度大于5 mm時(shí)，地層波基本消失，此時(shí)波列主頻即為套管波主頻（17～18 kHz）。此外，根據(jù)數(shù)值模擬結(jié)果，砂巖地層條件下，通過波形圖在時(shí)域內(nèi)也很難提取地層波信息。因此，應(yīng)利用全波列主頻識(shí)別地層波，為研究地層波信息與Ⅱ界面膠結(jié)質(zhì)量關(guān)系提供數(shù)據(jù)支撐。

3 基于聲波-變密度測(cè)井的固井質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)

固井質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)一般通過2類方法建立，一類根據(jù)水泥強(qiáng)度、膠結(jié)強(qiáng)度與測(cè)井信息，從橫向上對(duì)套管-水泥-地層三者的膠結(jié)程度進(jìn)行評(píng)價(jià)；另一類根據(jù)水泥可阻止流體通過的性能，從縱向上對(duì)水泥的封隔能力進(jìn)行評(píng)價(jià)。國(guó)內(nèi)各油田根據(jù)地層、鉆井液、水泥漿等實(shí)際情況，制定了各具特色的固井質(zhì)量評(píng)價(jià)方法和標(biāo)準(zhǔn)，但總體都來源于利用聲波測(cè)井評(píng)價(jià)固井質(zhì)量的行業(yè)標(biāo)準(zhǔn)［19］，但該標(biāo)準(zhǔn)僅考慮了聲波幅度，無法對(duì)Ⅱ界面膠結(jié)質(zhì)量進(jìn)行分析。本文根據(jù)套管井聲場(chǎng)研究結(jié)果，提出新的基于CBL/VDL測(cè)井的固井質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)。

3.1 評(píng)價(jià)參數(shù)

Ⅰ界面膠結(jié)指數(shù)定義如下：

水泥環(huán)膠結(jié)比定義如下：

依據(jù)套管井聲場(chǎng)研究結(jié)果及（15）式、（16）式，得出了4種水泥漿密度下Ⅰ界面膠結(jié)指數(shù)與水泥環(huán)膠結(jié)比關(guān)系圖版（見圖8a），可以看出，Ⅰ界面膠結(jié)指數(shù)與水泥環(huán)膠結(jié)比近似呈線性關(guān)系。一般認(rèn)為，水泥環(huán)膠結(jié)比大于0.8（即竄槽角小于72°）為膠結(jié)好，則圖8a中水泥漿密度為1.20 g/cm3時(shí)，Ⅰ界面膠結(jié)指數(shù)大于0.75為膠結(jié)好；水泥漿密度為1.89 g/cm3時(shí)，Ⅰ界面膠結(jié)指數(shù)大于0.80為膠結(jié)好；水泥漿密度為2.25 g/cm3時(shí)，Ⅰ界面膠結(jié)指數(shù)大于0.82為膠結(jié)好。因此，低密度水泥漿固井質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)應(yīng)適當(dāng)降低，而高密度水泥漿條件下應(yīng)適當(dāng)提高評(píng)價(jià)標(biāo)準(zhǔn)。

在Ⅰ界面膠結(jié)好的情況下，地層波信息與Ⅱ界面膠結(jié)質(zhì)量關(guān)系密切。由于套管波、面波的干擾，直接提取地層波比較困難，地層波幅度采用文獻(xiàn)［20］的方法計(jì)算，并經(jīng)巖性校正。同時(shí)，定義Ⅱ界面膠結(jié)指數(shù)如下：

圖8　界面膠結(jié)指數(shù)與水泥環(huán)膠結(jié)比關(guān)系

根據(jù)套管井聲場(chǎng)研究結(jié)果，Ⅱ界面竄槽角度對(duì)地層波幅度有較大影響，由此把地層波幅度轉(zhuǎn)化為Ⅱ界面膠結(jié)指數(shù)，把Ⅱ界面竄槽角度轉(zhuǎn)化為水泥環(huán)膠結(jié)比后建立兩者關(guān)系（見圖8b），可以看出，兩者近似呈線性關(guān)系。若認(rèn)為水泥環(huán)膠結(jié)比大于0.80為膠結(jié)好，則Ⅱ界面膠結(jié)指數(shù)要大于0.78；若水泥環(huán)膠結(jié)比小于0.50為膠結(jié)差，則Ⅱ界面膠結(jié)指數(shù)要小于0.45。

3.2 評(píng)價(jià)標(biāo)準(zhǔn)

根據(jù)刻度井群CBL/VDL測(cè)井?dāng)?shù)據(jù)分析與數(shù)值模擬結(jié)果，利用聲波幅度、Ⅰ界面膠結(jié)指數(shù)和Ⅱ界面膠結(jié)指數(shù)建立了新的固井質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)（見表2、表3）。

表2　常規(guī)密度（1.5～1.8 g/cm3）水泥漿固井質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)

4 結(jié)論

采用實(shí)際刻度井群聲波-變密度測(cè)井與數(shù)值模擬結(jié)合的方法開展套管井聲場(chǎng)研究。套管井聲場(chǎng)可以歸納為軸對(duì)稱聲場(chǎng)與非軸對(duì)稱聲場(chǎng)兩類問題，分別采用實(shí)軸積分法和2.5維有限差分算法對(duì)兩者進(jìn)行數(shù)值模擬。采用2.5維有限差分算法計(jì)算得到的套管首波幅度隨竄槽角度的變化趨勢(shì)與對(duì)應(yīng)刻度井實(shí)測(cè)結(jié)果一致，驗(yàn)證了數(shù)值模擬方法的可靠性。

井中實(shí)測(cè)波形及數(shù)值模擬結(jié)果均表明，在現(xiàn)有聲波-變密度測(cè)井儀器探測(cè)精度下，只有當(dāng)竄槽角度大于45°時(shí)，才能利用聲波進(jìn)行檢測(cè)。同時(shí)，水泥漿密度不同，套管波幅度隨竄槽角度變化程度存在差別，在進(jìn)行聲波檢測(cè)時(shí)應(yīng)考慮水泥漿密度的影響。

在井眼條件及聲波-變密度測(cè)井儀器適當(dāng)時(shí)，利用套管波信息可以進(jìn)行Ⅰ界面膠結(jié)質(zhì)量檢測(cè)，利用地層波信息可進(jìn)行Ⅱ界面膠結(jié)質(zhì)量檢測(cè)。本文中，當(dāng)Ⅱ界面環(huán)隙寬度小于5 mm時(shí)，石灰?guī)r地層中地層波可作為首波接收，在頻譜上容易被識(shí)別出來，且地層波與Ⅱ界面環(huán)隙寬度關(guān)系密切；當(dāng)Ⅱ界面環(huán)隙寬度大于5 mm時(shí)，石灰?guī)r地層條件下無法識(shí)別地層波，也無法實(shí)現(xiàn)Ⅱ界面固井質(zhì)量的有效評(píng)價(jià)。

根據(jù)本文刻度井實(shí)測(cè)及數(shù)值模擬結(jié)果，對(duì)現(xiàn)有利用聲波-變密度測(cè)井信息建立的固井質(zhì)量評(píng)價(jià)標(biāo)準(zhǔn)進(jìn)行改進(jìn)，建立了包含Ⅰ、Ⅱ界面膠結(jié)指數(shù)的2類水泥漿密度條件下的評(píng)價(jià)標(biāo)準(zhǔn)。

致謝：本研究得到了中國(guó)石油西部鉆探工程有限公司測(cè)井公司在刻度井群實(shí)驗(yàn)數(shù)據(jù)采集上提供的支持與幫助，在此表示感謝。

符號(hào)注釋：

A——目的層地層波振幅譜的極值面積，m2；Amax——完全膠結(jié)井段地層波振幅譜的極值面積，m2；Amin——Ⅱ界面未膠結(jié)井段地層波振幅譜的極值面積，缺失時(shí)取1，m2；A1——反射系數(shù)；BI1，BI2——Ⅰ界面、Ⅱ界面膠結(jié)指數(shù)，f；C——套管波幅度；Cmax——自由套管井段套管波幅度；Cmin——完全膠結(jié)井段套管波幅度；d——接收源距，m；fmax——聲源最高有效頻率，Hz；F（t）——聲源的時(shí)域函數(shù)；gjk——單極子聲源函數(shù)矩陣；Gxx，Gyy，Gzz——聲源函數(shù)關(guān)于z軸的傅立葉變換，均為x，y，t，Kz的函數(shù)，Pa/s；I0——第一類0階虛宗量貝塞爾函數(shù)；Ic——水泥環(huán)膠結(jié)比，f；k——波數(shù)，m-1；kr——波數(shù)的徑向分量，m-1；kz——z軸方向波數(shù)分量，m-1；kzmax——最大軸向波數(shù)，m-1；Kz——空間坐標(biāo)z的傅立葉變換；p——壓力響應(yīng)，Pa；r——與井軸的距離，m；r——與井軸的距離矢量，m；t——時(shí)間，s；Txx，Tyy，Tzz，Txy，Txz，Tyz——應(yīng)力分量沿z軸的傅立葉變換，均為x，y，t，Kz的函數(shù)，Pa；vf——聲波速度，m/s；vmax，vmin——介質(zhì)中最大、最小聲速，m/s；Vx，Vy，Vz——x，y，z方向的速度沿z軸的傅立葉變換，m/s；x，y，z——直角坐標(biāo)系，m；X（ω）——聲源頻譜函數(shù)；δ——單位階躍函數(shù)；δjk——單位階躍函數(shù)矩陣；Δt——時(shí)間步長(zhǎng)，s；Δx，Δy，Δz——空間步長(zhǎng)，m；θ——竄槽角度，（°）；λ——一階拉梅常數(shù)，Pa；μ——二階拉梅常數(shù)，Pa；ρ——密度，kg/m3；ω——聲頻率，Hz。

參考文獻(xiàn)：

［1］章成廣， 江萬哲， 潘和平.聲波測(cè)井原理及應(yīng)用［M］.北京： 石油工業(yè)出版社， 2009.ZHANG Chengguang， JIANG Wanzhe， PAN Heping.Acoustic logging and application［M］.Beijing： Petroleum Industry Press， 2009.

［2］李濤.高溫高壓套損井膨脹管修復(fù)技術(shù)［J］.石油勘探與開發(fā)，2015， 42（3）： 374-378.LI Tao.Solid expandable tubular patching technique for high-temperature and high-pressure casing damaged wells［J］.Petroleum Exploration and Development， 2015， 42（3）： 374-378.

［3］BIOT M A.Propagation of elastic waves in a cylindrical bore containing a fluid［J］.Journal of Applied Physics， 1952， 23（9）： 997-1005.

［4］BIOT M A.Theory of propagation of elastic waves in a fluid-saturated porous solid（II）： Higher frequency range［J］.Journal of the Acoustical Society of America， 1956， 28（2）： 179-191.

［5］WHITE J E， ZECHMAN R E.Computed response of an acoustic logging to0l［J］.Geophysics， 1968， 33（2）： 302-310.

［6］CHENG C H， TOKSOZ M N.Elastic wave propagation in a fluid filled borehole and synthetic acoustic logs［J］.Geophysics， 1981，46（7）： 1042-1053.

［7］ZHANG Hailan， LI Mingxuan， HAILAN Z.Numerical study of acoustic field of borehole in slow formation［J］.CJ Geophysics， 1993，36（1）： 137-144.

［8］胡文祥， 錢夢(mèng)騄.套管井聲場(chǎng)合成波形及其時(shí)頻特征分析［J］.聲學(xué)學(xué)報(bào)， 2002， 27（3）： 223-228.HU Wenxiang， QIAN Menglu.Synthetic waveforms of the acoustic field in cased boreholes and their time-frequency distribution［J］.ACTA Acoustic， 2002， 27（3）： 223-228.

［9］林偉軍， 張澄宇， 張海瀾， 等.帶扇形竄槽的套管井聲場(chǎng)［J］.聲學(xué)學(xué)報(bào)， 2005， 30（1）： 9-14.LIN Weijun， ZHANG Chengyu， ZHANG Hailan， et al.Acoustic field in a cased well with a sectorial crossing channel［J］.ACTA Acoustic，2005， 30（1）： 9-14.

［10］SONG R L， LIU J X， YAO G J.Parallel finite difference modeling of acoustic fields in nonaxisymmetric cased hole［J］.Chinese Journal of Geophysics， 2010， 53（6）： 2767-2775.

［11］SONG R L， LIU J X， HOU C H.Numerical simulation of sector bond log and improved cement bond image［J］.Geophysics， 2012， 77（4）：95-104.

［12］陳德華， 王秀明， 張海瀾， 等.水泥密度和套管尺寸對(duì)油井套管波的影響［J］.聲學(xué)學(xué)報(bào)， 2008， 33（1）： 15-20.CHEN Dehua， WANG Xiuming， ZHANG Hailan， et al.The effects of cement density and casing dimension on casing waves in oil wells［J］.ACTA Acoustic， 2008， 33（1）： 15-20.

［13］SHI Fangfang， WU Xianmei， ZHANG Bixing.Application of cylindrical linear phased array in casing borehole［J］.Chinese Journal of Acoustics， 2010， 29（1）： 65-72.

［14］韓煒， 毛捷， 金士杰.套管-水泥界面微間厚度的低頻超聲反演方法［J］.聲學(xué)學(xué)報(bào)， 2014， 39（4）： 467-472.HAN Wei， MAO Jie， JIN Shijie.A thickness determination method of thin layer between casing and cement using low frequency ultrasound［J］.ACTA Acoustic， 2014， 39（4）： 467-472.

［15］刁順， 喬文孝， 杜光升.聲波測(cè)井對(duì)微環(huán)探測(cè)能力的研究［J］.石油地球物理勘探， 2003， 38（5）： 540-542.DIAO Shun， QIAO Wenxiao， DU Guangsheng.Influence of microannulus thickness on capacity of acoustic wave-long tool［J］.Oil Geophysical Prospecting， 2003， 38（5）： 540-542.

［16］徐琰鋒， 胡文祥.縱向帶狀裂隙形貌的逆時(shí)偏移超聲成像［J］.物理學(xué)報(bào)， 2014， 63（15）： 154302.XU Yanfeng， HU Wenxiang.Ultrasonic imaging for appearance of vertical slot by reverse time migration［J］.Acta Physica Sinica， 2014，63（15）： 154302.

［17］初緯， 沈吉云， 楊云飛， 等.連續(xù)變化內(nèi)壓下套管-水泥環(huán)-圍巖組合體微環(huán)隙計(jì)算［J］.石油勘探與開發(fā)， 2015， 42（3）： 379-385.CHU Wei， SHEN Jiyun， YANG Yunfei， et al.Calculation of micro-annulus size in casing-cement sheath-formation system under continuous internal casing pressure change［J］.Petroleum Exploration and Development， 2015， 42（3）： 379-385.

［18］劉繼生， 王克協(xié).用頻率-波數(shù)分析研究聲波測(cè)井全波列的各相波［J］.測(cè)井技術(shù)， 2000， 24（3）： 198-202.LIU Jisheng， WANG Kexie.Studying each mode of ascoutic full-wavetrains using frequency wavenumber analysis［J］.Well Logging Technology， 2000， 24（3）： 198-202.

［19］中華人民共和國(guó)國(guó)家發(fā)展和改革委員會(huì).固井質(zhì)量評(píng)價(jià)方法：SY/T 6592-2004［S］.北京： 石油工業(yè)出版社， 2004.National Development and Reform Commission.Procedure for cement evaluation： SY/T 6592-2004 ［S］.Beijing： Petroleum Industry Press， 2004.

［20］唐軍， 章成廣.MAK-Ⅱ聲波測(cè)井資料在定量評(píng)價(jià)固井質(zhì)量中的應(yīng)用［J］.測(cè)井技術(shù)， 2011， 35（3）： 266-269.TANG Jun， ZHANG Chengguang.Application study on cement bond quantitative evaluation by MAK-Ⅱwaveform data［J］.Well Logging Technology， 2011， 35（3）： 266-269.

（編輯 胡葦瑋）

Cement bond quality evaluation based on acoustic variable density logging

TANG Jun1， 2， ZHANG Chengguang1， 2， ZHANG Bixing3， SHI Fangfang3
（1.Key Laboratory of Exploration Technologies for Oil and Gas Resources， Ministry of Education， Yangtze University，Wuhan 430100， China； 2.School of Geophysics and Oil resources， Yangtze University， Wuhan 430100， China；
3.Institute of Acoustics， Chinese Academy of Sciences， Beijing 100190， China）

Abstract：A new method of cement bond quality evaluation was proposed by combining numerical simulation and calibrated cased hole acoustic logging data.The effects of the cement channel angle and the quality of the second bond interface （the interface of cement with formation） on acoustic variable density logging data were analyzed.Based on the analysis result， a new cement bond evaluation standard was presented after revising the traditional CBL/VDL method.The axisymmetric acoustic field was simulated by real axis integral method， while the non-axisymmetric acoustic field was simulated by 2.5-D finite differential method.After comparing with the calibrated cased hole acoustic logging data， the research has the below results： the numerical simulation result matches with the calibrated well logging data very well and the new method is reliable； the amplitude of the first acoustic arrival in the case hole decreases as the angle of cement channel decreases， and the denser the cement is， the faster the amplitude of cased hole acoustic waveform decays； the lower limit of cement channel angle is around 45 degrees which can be detected by acoustic logging； the formation acoustic waveform is not easy to be detected in time domain， however it is easy to be detected in frequency domain， especially in limestone formation， the first arrival only can be detected when the annulus width of the second bond interface is small.According to the research result of the numerical simulation of cased hole acoustic field and acoustic variable density logging data， new evaluation criteria of cement bound quality were presented.

Key words：cement bond quality； acoustic variable density logging； acoustic field in cased well； cement channel angel； interface bond index

中圖分類號(hào)：P631.8

文獻(xiàn)標(biāo)識(shí)碼：A

文章編號(hào)：1000-0747（2016）03-0469-07

DOI：10.11698/PED.2016.03.19

基金項(xiàng)目：國(guó)家自然科學(xué)基金（41274141）；油氣資源與勘探技術(shù)教育部重點(diǎn)實(shí)驗(yàn)室（長(zhǎng)江大學(xué)）開放基金（K2015-04）

第一作者簡(jiǎn)介：唐軍（1979-），男，湖北荊門人，博士，長(zhǎng)江大學(xué)地球物理與石油資源學(xué)院講師，主要從事地球物理測(cè)井方面的教學(xué)與研究工作。地址：湖北省武漢市蔡甸區(qū)大學(xué)路111號(hào)，長(zhǎng)江大學(xué)地球物理與石油資源學(xué)院，郵政編碼：430100。E-mail：tangjun@yangtzeu.edu.cn

收稿日期：2015-03-26 修回日期：2016-02-25