含水飽和度對致密砂巖聲波特性的影響研究

黃裕萌，熊 健，黃開樺，王 佳

（西南石油大學油氣藏地質及開發工程國家重點實驗室，四川成都 610500）

伴隨著世界油氣資源需求量的增長、常規油氣資源開采量的下降，具有較大儲量、分布廣的非常規油氣開采獲得了越來越多國家的關注。非常規油氣的開采難度、極低的產能，由于致密砂巖孔隙度極低、滲透率小，其開采難度大[1]。因此，致密砂巖的各向異性、聲波特性等方面的研究對致密砂巖氣藏的勘探開發具有重要的意義。

國內學者利用巖石的聲波特性對儲層的巖石進行了較為系統的實驗研究。張明明[2]，高可攀[3]，陳旭等[4]利用了超聲波時域信號分析的方法對不同巖性和含水飽和度的巖石的聲波特性、孔隙壓力等因素進行了研究；宋連騰等[1]研究了飽和水與干燥條件下，不同圍壓對致密砂巖彈性各向異性的影響，缺少了控制巖石在不同含水飽和度的情況；朱合華[5]，周治國[6]通過致密板巖、花崗片麻巖、黑云母片麻巖的巖石聲波的主頻頻率和波速進行研究，得到飽和水后，三種巖石的縱波波譜都發生變化，縱波主頻明顯降低；趙宇等[7]，成林等[8]考慮了煤巖的吸水率及各向異性對縱、橫波速度的影響，但是卻忽略了巖石聲波的頻譜特征里面也含有大量信息，并未對巖石的各向異性進行闡述。

本文以四川盆地須家河的致密砂巖為研究對象，研究了含水飽和度對致密砂巖聲波特性的影響，并進一步解釋了聲波傳播特性隨含水飽和度的內在機理。在此基礎上，也討論了含水飽和度對致密砂巖儲層聲波各向異性的影響。

1 實驗樣品和實驗方法

1.1 實驗樣品

以四川盆地須家河組致密砂巖為研究對象，從直徑的全直徑巖心上鉆取巖心，先在該巖心上沿圓周360°測量巖心的波速，以10°一組波速，得到波速與圓周角度的關系，從中獲得最大波速和最小波速的角度，以此鉆取巖心，并在兩個角度之間呈45°方向鉆取巖心，總共在全直徑3 個方向鉆取巖心，并假定波速最大方向為X 方向，波速最小方向為Z 方向，夾角呈45°方向為Y 方向。

測量巖樣長度、直徑、干重等參數，利用SCAR-Ⅱ型高溫高壓測量儀對巖樣孔隙度和滲透率進行測量，具體參數（見表1）。

表1 巖石的基本參數

1.2 含水飽和度

采用超聲波透射法對干燥的巖心樣品進行超聲波測試，用激發頻率為100 kHz 的超聲波探頭發出超聲波脈沖，采集接收端探頭聲波的波形數據。根據樣品的長度和超聲波通過樣品的時間，按下式計算巖樣樣品的波速：

將巖樣抽真空后在地層水中充分飽和，并靜置24 h，建立巖樣地層水飽和點。采用自然風干法分別獲得含水80 %，含水60 %及束縛水狀態下的巖樣。含水飽和度計算公式如下：

式中：Sw-含水飽和度；m烘干-巖心烘干后質量；m飽和-巖樣浸泡24 h 后的質量；m-巖心自然風干過程中的質量。

1.3 頻譜分析方法[9]

假設有一個周期信號x（t），其周期為2L，傅里葉變換的離散形式可以表示為：

其中：

式中：a0、ak、bk-傅里葉系數；N-變換區間長度；m-取樣點數。

根據上述原理將聲波探頭采樣到的數據轉化為頻域數據，在頻域上對聲波信號進行分析。利用上述理論，編制快速傅里葉變換程序，對采集的聲波時域離散信號進行分解，得到信號的頻域特征，探頭對接時的時域波形（見圖1），快速傅里葉變換后的頻域波形（見圖2）。從波形圖可以看出，探頭對接時域波形呈漏斗狀，有拖尾，頻譜圖的主頻率突出，為100 kHz，副頻率集中在25 kHz、210 kHz、400 kHz 附近。

圖1 探頭對接時的時域波形

圖2 快速傅里葉變換后的頻域波形

2 實驗結果

2.1 飽和水巖心頻域圖變化

用發射頻率為100 kHz 的超聲波脈沖通過干燥條件下和飽和水條件下的巖心，收集一個完整周期內的時域信號。將接收到的時域信號經快速傅里葉變換轉換為頻域信號，并繪出相應的頻譜圖，進而從時域和頻域2 個方面分別考察巖樣濾波作用的差異（見圖3～圖5）。

從圖3～圖5 中可以看出：在相同的發射脈沖條件下，干燥巖樣的波形頻域不光滑，主頻率突出，信號能量集中在100 kHz 附近，還存在多個低峰值；對巖石進行飽和水處理后，主頻率較干燥條件下呈現不同程度的衰減，高頻和低頻區域的能量衰減程度較大，不同方向的高頻區域和低頻區域的能量衰減程度不一樣。

圖3 X 方向上的聲波頻譜圖

圖4 Y 方向上的聲波頻譜圖

圖5 Z 方向上的聲波頻譜圖

利用快速傅里葉變換對原始信號進行處理，得到的頻譜圖能夠清楚地確定出原始信號所包含頻率值的大小。但是對于確定頻域突變點的位置，傅里葉變換顯得不足。為了確定突變點的時間信息，進一步了解巖塊內部的裂隙信息，采用Daubechies 6 小波在尺度1～35進行連續小波變換，得到的相應系數絕對值的圖像（見圖6～圖8）。

圖6 X 方向上Db6 連續小波變換系數圖

圖7 Y 方向上Db6 連續小波變換系數圖

圖8 Z 方向上Db6 連續小波變換系數圖

從圖6～圖8 可以看到，小波分解后的系數圖像明顯地分為明亮和陰暗兩種不同的色度。明亮的部分代表該巖石聲波的突變點比較少，甚至不存在；陰暗的部分代表該巖石聲波具有明顯的異常點，并且隨著分解的層數增大，該突變點會被放大的越明顯。反映到巖石本身，明亮的部分越多，代表巖石內部更加完整，裂隙發育比較少；而對應的陰暗部分則代表該巖石內部存在裂隙和孔洞，陰暗的部位越多，說明裂隙發育良好，對比表1 的孔隙率，發現符合的很好。

圖9 X 方向的不同飽和度的頻譜圖

2.2 飽和度變化的影響

圖10 Y 方向的不同飽和度的頻譜圖

圖11 Z 方向的不同飽和度的頻譜圖

表2 通過不同飽和度巖心的聲波速度

表3 不同飽和度對應的衰減系數

從圖9～圖11 的頻譜圖可以看出，隨著巖樣的含水飽和度的不斷增大，所有方向上巖心的聲波主頻率的能量都隨著減少，不同方向的高頻區域和低頻區域的能量衰減程度不一。究其原因，可能是因為不同方向的孔隙分布不一，水在孔隙中分布存在差異，造成聲波在其中傳播特征存在差異，對聲波能量的衰減影響作用不一樣，這說明了含水飽和度對不同方向的能量衰減的影響程度存在差異。

為了更加直觀地分析致密砂巖在不同含水飽和度下聲波的各向異性的變化情況（見表2、表3），將表2中的通過不同飽和度巖心的聲波速度和表3 中的不同飽和度對應的衰減系數進行了匯總。

圖12 不同飽和度對應的聲波速度

圖13 不同飽和度對應的衰減系數

從圖12、圖13 可以看出，隨著含水飽和度增加，不同方向的波速和衰減系數增大，且不同方向的增加幅度不一樣。這可能是因為在隨著巖石內部的孔隙逐漸被液體充填，聲波在通過水和巖石之間形成的波阻抗界面時將發生不同程度的折射和反射，造成波的能量在傳播的過程中發生衰減。隨著飽和度的增加，這種現象發生的更明顯，在不同方向上的衰減程度不同也體現了巖石的各向異性，說明了含水飽和度增加對不同方向的波速和衰減系數增加影響程度不一。

3 結論

（1）隨著巖石含水飽和度的增加，巖石中水的分布區域增多，聲波在巖石傳播過程中產生反射、折射現象增多，造成聲波衰減系數增大。

（2）相比于干燥情況下的巖樣，飽和水后巖樣的波速增大，而主頻率減小，聲波信號高頻部分能量嚴重衰減，低頻部分能量增加。

（3）不同方向的致密砂巖樣品的聲波速度和衰減系數隨著含水飽和度增加而增加，但增加幅度不一致，說明了含水飽和度對致密砂巖巖樣的各向異性存在較大的影響。