[發(fā)明專利]裂縫性三維數(shù)字巖心油水或氣水分布模型的構(gòu)建方法

有效

• 申請(qǐng)?zhí)枺?/td>: 202010873878.1
• 申請(qǐng)日: 2020-08-26
• 公開（公告）號(hào): CN112131776B
• 公開（公告）日: 2023-05-05
• 發(fā)明（設(shè)計(jì)）人: 王謙;譚茂金;王思宇;張令坦
• 申請(qǐng)（專利權(quán)）人: 中國(guó)地質(zhì)大學(xué)（北京）
• 主分類號(hào): G06F30/25
• 分類號(hào): G06F30/25;G06T7/11;G06T7/136;G06T17/00;G06Q50/02;G06F111/10
• 代理公司: 北京興智翔達(dá)知識(shí)產(chǎn)權(quán)代理有限公司 11768
• 代理人: 郭衛(wèi)芹

1.一種裂縫性三維數(shù)字巖心油水或氣水分布模型的構(gòu)建方法，其特征在于，包括：

步驟1：基于巖心樣本的X-CT掃描技術(shù)獲得巖心掃描圖像，通過濾波、閾值分割手段重構(gòu)三維數(shù)字巖心；

步驟2：利用三維圖像處理軟件，構(gòu)建三維裂縫模型；

步驟3：利用數(shù)學(xué)形態(tài)學(xué)中的開運(yùn)算，選擇一個(gè)結(jié)構(gòu)元素半徑處理基質(zhì)孔隙的油水或氣水分布，得到第一油水或氣水分布模型，計(jì)算各組分體積，記錄飽和度S_w0；

步驟4：利用三維圖像處理技術(shù)將步驟2構(gòu)建的裂縫模型和步驟3處理后的第一油水或氣水分布模型進(jìn)行組合，得到第二油水或氣水分布模型，重新計(jì)算飽和度，記為S_w1；

步驟5：利用測(cè)井資料和巖石物理資料確定束縛水飽和度大小；當(dāng)所需的含水飽和度小于或等于束縛水飽和度時(shí)，認(rèn)為裂縫中飽含油或氣，水則以束縛水的狀態(tài)全部賦存于基質(zhì)孔隙中；當(dāng)所需含水飽和度大于束縛水飽和度時(shí)，按步驟6進(jìn)行處理；

步驟6：利用LBM模擬裂縫空間的油水或氣水分布，即步驟2所構(gòu)建的三維裂縫模型，得到第三油水或氣水分布模型；確定輸入含水飽和度大小，選擇邊界條件和迭代步數(shù)模擬裂縫的油水或氣水分布，飽和度記為Sw2；Sw2的值取決于油氣儲(chǔ)層基本特征、流體的分布情況以及Sw1的大小；

步驟7：利用三維圖像處理軟件將開運(yùn)算模擬的基質(zhì)孔隙油水或氣水分布結(jié)果與LBM模擬的裂縫中的油水或氣水分布結(jié)果進(jìn)行組合，構(gòu)建裂縫性三維數(shù)字巖心模型，利用公式(1)計(jì)算最終飽和度S_w，即為該飽和度下的裂縫性三維數(shù)字巖心油水或氣水分布模型；公式（1）如下：?(1)

其中，S_w是裂縫性三維數(shù)字巖心油水或氣水分布模型的含水飽和度，S_w1是第二油水或氣水分布模型的含水飽和度，S_w2是第三油水或氣水分布模型的含水飽和度，φ₁是第二油水或氣水分布模型中基質(zhì)的孔隙度，φ₂是第三油水或氣水分布模型中裂縫的孔隙度；

步驟8：改變結(jié)構(gòu)元素半徑大小，重復(fù)步驟3到步驟7，即可構(gòu)建不同含水飽和度的裂縫性三維數(shù)字巖心。

2.一種電子設(shè)備，其特征在于，包括存儲(chǔ)器和處理器；

所述存儲(chǔ)器中存儲(chǔ)有計(jì)算機(jī)程序；

所述處理器，用于執(zhí)行所述計(jì)算機(jī)程序以實(shí)現(xiàn)權(quán)利要求1所述的方法。

3.一種計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)，其特征在于，所述計(jì)算機(jī)可讀存儲(chǔ)介質(zhì)上存儲(chǔ)有計(jì)算機(jī)程序，所述計(jì)算機(jī)程序被處理器執(zhí)行時(shí)實(shí)現(xiàn)權(quán)利要求1所述的方法。