[發明專利]一種基于Matlab分析碳酸鹽巖結構的方法

有效

• 申請號: 201810571924.5
• 申請日: 2018-06-06
• 公開（公告）號: CN108734714B
• 公開（公告）日: 2022-11-25
• 發明（設計）人: 久博;黃文輝;何明倩;劉凱;賈瑜
• 申請（專利權）人: 中國地質大學（北京）
• 主分類號: G06T7/11
• 分類號: G06T7/11;G06T7/60
• 代理公司: 北京惠科金知識產權代理有限公司 11981
• 代理人: 袁曉哲

本發明公開了一種基于Matlab分析碳酸鹽巖結構的方法，包括如下步驟：首先讀取碳酸鹽巖的鑄體薄片；然后選擇色彩模式，將碳酸鹽巖的鑄體薄片圖像二值化；再識別二值圖像的黑色像素點，計算孔隙度，同時識別二值圖像的白色像素點，分辨顆粒大小，計算粒度；然后求取碳酸鹽巖的粒度與孔徑大小頻率分布圖，最后求取碳酸鹽巖粒度與孔徑的平均直徑及方差，并輸出二值圖像和求取的參數，通過Matlab定量計算不同碳酸鹽孔隙結構，平均誤差率為+4.661％，其中粒間孔平均誤差率為?4.370％，鑄模孔平均誤差率為+6.734％，生屑孔誤差率為+6.144％，微裂縫平均誤差率為?0.314％，鑄模孔與生屑孔的孔隙結構最為適用。

技術領域

本發明涉及致密砂巖儲層孔隙表征分析技術領域，具體為一種基于Matlab分析碳酸鹽巖結構的方法。

背景技術

致密砂巖儲層通常為儲層滲透率低的砂巖儲層，致密低滲透儲層是一個相對概念，世界上并無一固定的標準和界限，它是由不同國家不同時期的資源和技術經濟條件和世界油氣價格而決定的。

儲集層是指具有連通孔隙、允許油氣在其中儲存和滲濾的巖層。世界上已發現的油氣儲量大多數來自沉積巖層，其中以砂巖和碳酸鹽巖儲集層最為重要，裂縫性泥巖和煤層也可作為儲集層；火成巖和變質巖儲集層中也有工業性油氣發現。

儲層(Reservoir)的儲集能力是由儲集層的巖石物理性質決定的，通常包括其孔隙性、滲透性；孔隙性決定了儲層儲存能力的大小，滲透性決定了儲集物的滲流能力。

碳酸鹽巖是主要由方解石和白云石等碳酸鹽礦物組成的沉積巖，其巖石結構包括顆粒、膠結物、泥晶、孔隙幾部分，顆粒可分為內碎屑、鮞粒、球粒、藻粒、團塊這幾部分，巖石微觀孔隙結構是指巖石中孔隙和吼道的幾何形狀、大小、分布及其相互連通關系。國內外多名學者對碳酸鹽巖孔隙結構、表征等進行研究，多種手段被應用于研究致密砂巖孔隙，以核磁共振NMR測井技術、流體注入實驗法和圖像觀測法這三大類為主。基于巖石電阻率參數的核磁共振測井技術，以測定儲層孔徑分布為主。以壓汞法為代表的流體注入實驗法能夠間接地通過壓汞曲線獲得儲集空間分布、孔徑大小等參數，但并不能像以掃描電鏡SEM及納米CT技術等圖像觀測法直接地觀測到儲集空間分布、孔喉特征等特征。然而以掃描電鏡及納米CT為實驗指導，研究結果大多為對孔喉的幾何形狀及聯通關系進行定性描述，對孔隙表征的定量研究仍有不足。

發明內容

為了克服現有技術方案的不足,本發明提供一種基于Matlab分析碳酸鹽巖結構的方法，能有效的解決背景技術提出的問題。

本發明解決其技術問題所采用的技術方案是：

一種基于Matlab分析碳酸鹽巖結構的方法，包括如下步驟：

步驟100、讀取碳酸鹽巖的鑄體薄片；

步驟200、選擇色彩模式，將碳酸鹽巖的鑄體薄片圖像二值化；

步驟300、識別二值圖像的黑色像素點，計算孔隙度，同時識別二值圖像的白色像素點，分辨顆粒大小，計算粒度；

步驟400、求取碳酸鹽巖的粒度與孔徑大小頻率分布圖；

步驟500、求取碳酸鹽巖粒度與孔徑的平均直徑及方差，并輸出二值圖像和求取的參數。

作為本發明一種優選的技術方案，所述步驟200中，鑄體薄片圖像二值化的方法為：將鑄體薄片圖像導入Matlab的內置Color Threshod應用中通過HSV成像模式進行圖像二值化轉換。

作為本發明一種優選的技術方案，不同的成像模式適用于不同圖像色彩空間，最終目的均為最大限度地識別碳酸鹽巖結構并將其轉為二值圖像，所述HSV模式因其可以調控圖像色調H、飽和度S以及明度V。