一種用于微觀仿真模型的圖像處理方法

劉曉玲

摘 要：微觀滲流模擬技術是一種從孔隙層次上觀察油水運動特征、分析油水滲流規律和研究微觀驅油機理的重要方法。文章通過建立微觀驅油圖像處理方法，將驅油劑的波及系數與采收率量化，為化學驅油效果評價與機理研究起到一定作用。

關鍵詞：微觀；波及系數；采收率

微觀滲流模擬技術利用微觀仿真模型進行驅油實驗，通過圖像采集和處理技術實現驅油過程的定性和定量分析，可以詳細了解到各驅油方式的滲流特征、剩余油特征和驅替效果[1-2]。與NMRI相比，它可以實時監測驅油細節，觀察油水界面作用現象，是三采微觀驅油機理強有力的研究手段。目前微觀滲流模擬技術已應用于研究各種復雜滲流，發現和證實了多種微觀機理。例如在聚合物驅機理研究方面，原有理論認為聚合物驅只能提高波及系數，而不能提高驅油效率[3]。1990年滲流所通過微觀模擬發現，聚合物的粘彈性可使這兩者都有一定程度的提高[4]。由此可見，微觀模擬可以為化學驅提高原油采收率提供有力的理論支持。

本文建立了一種可視化仿真模型微觀驅油圖像的處理方法，將波及系數與原油采收率量化，為微觀化學驅驅油效果評價與機理研究提供一定支持。

1 原油采收率計算

通過分析微觀驅油初始、結束時刻的圖像，獲取原始含油飽和度與剩余油飽和度信息，計算原油采收率。

（1）計算模型孔隙度φ

取飽和水后的模型圖（見圖1-a），將其轉化為二值圖（圖1-b），取黑色（模型骨架）和白色（水/孔隙）的像素數目，模型孔隙度計算公式為：

式中，N1w和N1b分別為飽和水后模型反相二值圖中白色（油/孔隙）、黑色（模型骨架）的像素數目。

（a）原圖 （b）反相二值圖

圖1 飽和水后的模型圖

（2）計算含油飽和度So

取任意驅替時刻的模型圖（見圖2-a），由其反二值圖（見圖2-b）來計算此時的含油飽和度，公式為：

式中，N2w和N2b分別為此刻模型反相二值圖中白色（油/孔隙）、黑色（模型骨架）的像素數目。

（a）原圖 （b）反相二值圖

圖2 飽和油后的模型圖

由此，可以用類似方法計算原始含油飽和度Soi與殘余油飽和度Sor。然后可按下式計算采收率R：R=（Soi-Sor）/Soi。

2 波及系數計算

微觀可視化試驗的最大特點是可以直觀地獲得波及系數信息。但驅替過程中剩余油的形態比較復雜，既有油膜、油滴、油塊，又有棕褐色乳狀液，而且在驅替中有的孔隙中已有連續水流通過，但隨著驅替過程的進行，后續油流可能會在孔隙中滯留，表現為孤立油滴或油膜狀、油滴狀殘余油，粘附于孔隙壁上。

這些剩余油以不同的面積與驅油劑接觸，使得判斷某些微觀區域是否被驅油劑波及變得不容易。為定量計算波及系數，本文設定了以下標準逐一判斷某些區域是否被驅油劑波及，具體步驟如下：（1）連續水流通過的區域為驅油劑波及區域。（2）對于殘留有油膜的孔道，如果油膜厚度小于孔道半徑四分之一，則定義孔道整體為波及區域；否則，定義除油膜外的孔道區域為波及區域（圖4-Ⅰ）。（3）對于有油滴或油塊的孔道，如果油滴或油塊有3/4區域與連續水流接觸，則定義孔道整體為波及區域（圖4-Ⅱ）；否則，定義除油滴或油塊外的孔道區域為波及區域（圖4-Ⅲ）。（4）驅油劑滲入產生乳狀液區域：若滲入區域顏色明顯變淡，即為波及區域。

（a）原圖 （b）填充后

（c）分離出原油 （d）反相

圖4 無連續水流通過區域中油膜與孤立油滴波及面積說明

以圖4為例，將波及到的孔隙進行填充，利用圖像分割技術分離出原油（圖4-c），將圖像反相后（圖4-d）計算原油所占像素數目No，并計算模型骨架與未波及到的孔隙所占像素數目Nb與孔隙度？漬，波及系數計算公式為：

3 結論

（1）通過圖像處理方法計算原始含油飽和度、殘余油飽和度，實現原油采收率的量化；

（2）設定了波及區域的判斷步驟，并通過圖像處理將波及系數量化，為微觀化學驅驅油效果評價與機理研究提供一定支持。

參考文獻

[1]黃延章，于大森.微觀滲流實驗力學及其應用[M].北京：石油工業出版社，2001：61-68.

[2]黃延章，郭尚平，周娟，等.孔隙介質中堿水驅油滲流的微觀機理[M].北京：科學出版社，1990.

[3]F.H.波特曼.提高原油采收率技術[M].北京：石油工業出版社，1985.

[4]黃延章，于大森，張桂芳.聚合物驅油微觀機理研究[J].油田化學，1990，7（1）：57-60.