一種基于ct圖像的滲透率預(yù)測方法
【技術(shù)領(lǐng)域】
[0001] 本發(fā)明屬于石油、地質(zhì)科學(xué)技術(shù)領(lǐng)域，涉及到一種基于CT圖像的滲透率預(yù)測方法。
【背景技術(shù)】
[0002] 多孔介質(zhì)流動一直是地質(zhì)科學(xué)、油藏工程及地下水污染等相關(guān)領(lǐng)域研究的主要問 題。其中，巖石滲透率作為影響多孔介質(zhì)流動最重要的因素之一，是進(jìn)行巖石評價(jià)的主要物 性參數(shù)。巖石滲透率包括絕對滲透率和相對滲透率，它們直接決定了巖石的儲運(yùn)能力、流體 流動規(guī)律、流體分布及流動穩(wěn)定性等。除此之外滲透率也會對巖石的其它重要物理參數(shù)(如 毛細(xì)管力分布)產(chǎn)生影響。因此，準(zhǔn)確預(yù)測巖石滲透率對相關(guān)領(lǐng)域的發(fā)展具有重要意義。
[0003] 傳統(tǒng)巖石滲透率主要由地震、巖芯分析、測井和試井等方法得到，但它們有以下缺 點(diǎn):其一，測量結(jié)果與測量儀器精度有很大關(guān)系，不同的測量儀器和測量方法會產(chǎn)生較大差 異的實(shí)驗(yàn)結(jié)果;其二，許多測試方法也只能針對一個樣品測量，測量過程耗時。其三，測井、 試井等大尺度測量方法雖然測量結(jié)果準(zhǔn)確，但卻不易操作，耗時耗力，也無法評價(jià)巖石微觀 物性對滲透率的影響;其四，對于具有各向異性的巖石的滲透率測量困難。

【發(fā)明內(nèi)容】

[0004] 為了解決傳統(tǒng)巖石滲透率測量方法的局限與不足，本發(fā)明了提出一種基于CT圖像 的滲透率預(yù)測方法，目的是為了快速準(zhǔn)確的得到巖石的絕對滲透率值和相對滲透率值。測 量方法主要包括3個過程:CT掃描，巖石表征單元體積確定和滲流模擬。
[0005] 本發(fā)明的技術(shù)方案是:一種基于CT圖像的滲透率預(yù)測方法，包括以下步驟：
[0006] 1.CT 掃描
[0007] 利用工業(yè)X射線CT對實(shí)際巖石進(jìn)行掃描，得到數(shù)字巖芯圖像并對圖像進(jìn)行裁剪;將 數(shù)字巖芯圖像二值化，得到數(shù)字巖芯圖像的骨架和孔隙；
[0008] 2.確定數(shù)字巖芯圖像的表征單元體積
[0009] (1)從1 %巖芯體積開始對數(shù)字巖芯圖像隨機(jī)選取至少10個1 %巖芯體積的子體 積，計(jì)算每個子體積的孔隙度并記錄;然后以2%巖芯體積的速率逐步同步擴(kuò)大每個子體積 并記錄新得到子體積的孔隙度，直到所有子體積的孔隙度相對大小大于95%為止，記錄此 時子體積的體積VI;其中，子體積的孔隙度計(jì)算公式為：
：Φ代表孔隙 度，sunw#表子體積孔隙像素點(diǎn)總數(shù)，sumgrain代表子體積骨架像素點(diǎn)總數(shù)）。
[0010] (2)從1 %巖芯體積開始對數(shù)字巖芯圖像隨機(jī)選取至少10個1 %巖芯體積的子體 積，計(jì)算每個子體積的分形維數(shù)并記錄。然后以2%巖芯體積的速率逐步同步擴(kuò)大每個子體 積并記錄新得到子體積的分形維數(shù)，直到所有子體積的分形維數(shù)相對大小大于95%為止， 記錄此時子體積的體積V2。其中，子體積的分形維數(shù)可以利用計(jì)盒算法思想，使用基于3D圖 像的分形盒維數(shù)算法計(jì)算得出。
[0011 ] (3)從1 %巖芯體積開始對數(shù)字巖芯圖像隨機(jī)選取至少10個1 %巖芯體積的子體 積，計(jì)算每個子體積的滲透率張量并記錄。然后以2%巖芯體積的速率逐步同步擴(kuò)大每個子 體積并記錄新得到子體積的滲透率張量，直到所有子體積的滲透率張量其矩陣對應(yīng)值相對 大小都大于90%為止，記錄此時子體積的體積V3。其中，子體積的滲透率張量可以基于廣義 達(dá)西定律和質(zhì)量守恒方程，利用有限元思想計(jì)算得到。
[0012] ⑷比較VI，V2，V3,其中的最大體積就是表征單元體積。
[0013] 3.滲流模擬:根據(jù)步驟(4)得到的表征單元體積從數(shù)字巖芯圖像中隨機(jī)選取表征 單元體積，利用滲透率計(jì)算模型對表征單元體積進(jìn)行滲流過程數(shù)值模擬(如孔隙網(wǎng)絡(luò)模型， CFD模型等），得到表征單元體積的絕對滲透率和相對滲透率，此時的絕對滲透率值和相對 滲透率值便可作為實(shí)際巖石滲透率的預(yù)測值。
[0014] 本發(fā)明的有益效果是：
[0015] 1.過程簡單，很容易實(shí)現(xiàn)不同地質(zhì)條件下巖石滲透率的多次測量。
[0016] 2.測量方法可對具有各向異性的巖石的滲透率進(jìn)行準(zhǔn)確預(yù)測。
[0017] 3.測量結(jié)果可以直接作為數(shù)值模擬的輸入?yún)?shù)，測量方法適用于任何氣-液-固或 液-液-固三相或多相體系的巖石滲透率預(yù)測。
[0018] 4.本發(fā)明也可推廣到巖石其它物性的預(yù)測。
【附圖說明】
[0019] 圖1是工業(yè)X射線CT掃描巖石獲取數(shù)字巖芯圖像的示意圖。
[0020] 圖2是表征單元體積獲取的流程示意圖。
[0021 ]圖中：1工業(yè)X射線CT; 2實(shí)際巖石樣品；3數(shù)字巖芯圖像;4子體積；
[0022] 5計(jì)算機(jī);6表征單元體積。
【具體實(shí)施方式】
[0023]以下結(jié)合技術(shù)方案和附圖詳細(xì)敘述本發(fā)明的【具體實(shí)施方式】。
[0024] 實(shí)施例
[0025] CT 掃描
[0026] 利用工業(yè)X射線CT1對實(shí)際巖石樣品2進(jìn)行掃描，得到數(shù)字巖芯圖像3并對圖像進(jìn)行 裁剪;將數(shù)字巖芯圖像3二值化，得到數(shù)字巖芯圖像3的骨架和孔隙。
[0027]表征單元體積確定
[0028] (1)從1%巖芯體積開始對數(shù)字巖芯圖像3隨機(jī)選取至少10個子體積4,利用計(jì)算機(jī) 5計(jì)算每個子體積4的孔隙度并記錄。然后以2%巖芯體積的速率逐步同步擴(kuò)大每個子體積4 并記錄新得到子體積4的孔隙度，直到所有子體積4的孔隙度相對大小大于95 %為止，記錄 此時子體積4的體積VI;其中，子體積4的孔隙度計(jì)算公式為：
：Φ代表 孔隙度，smwA表子體積孔隙像素點(diǎn)總數(shù)，sumgrain代表子體積骨架像素點(diǎn)總數(shù)）。
[0029] (2)從1%巖芯體積開始對數(shù)字巖芯圖像3隨機(jī)選取至少10個1 %巖芯體積的子體 積4,利用計(jì)算機(jī)5計(jì)算每個子體積4的分形維數(shù)并記錄。然后以2%巖芯體積的速率逐步同 步擴(kuò)大每個子體積4并記錄新得到子體積4的分形維數(shù)，直到所有子體積4的分形維數(shù)相對 大小大于95%為止，記錄此時子體積4的體積V2。其中，子體積4的分形維數(shù)可以利用計(jì)盒算 法思想，使用基于3D圖像的分形盒維數(shù)算法計(jì)算得出。
[0030] (3)從1 %巖芯體積開始對數(shù)字巖芯圖像3隨機(jī)選取至少10個1 %巖芯體積的子體 積4，利用計(jì)算機(jī)5計(jì)算每個子體積4的滲透率張量并記錄。然后以2 %巖芯體積的速率逐步 同步擴(kuò)大每個子體積4并記錄新得到子體積4的滲透率張量，直到所有子體積4的滲透率張 量其矩陣對應(yīng)值相對大小都大于90 %為止，記錄此時子體積4的體積V3。其中，子體積4的滲 透率張量可以基于廣義達(dá)西定律和質(zhì)量守恒方程，利用有限元思想計(jì)算得到。
[0031] (4)比較VI，V2，V3的大小，其中的最大體積就是表征單元體積6的大小。
[0032]滲流模擬
[0033]根據(jù)步驟(4)得到的表征單元體積6的大小從數(shù)字巖芯圖像中隨機(jī)選取表征單元 體積6,利用滲透率計(jì)算模型對表征單元體積6進(jìn)行滲流過程數(shù)值模擬（如孔隙網(wǎng)絡(luò)模型， CFD模型等），得到表征單元體積6的絕對滲透率和相對滲透率，此時的絕對滲透率值和相對 滲透率值便可作為實(shí)際巖石樣品2滲透率的預(yù)測值。
【主權(quán)項(xiàng)】
1. 一種基于CT圖像的滲透率預(yù)測方法，其特征在于，步驟如下： (1) CT掃描 利用工業(yè)X射線CT對實(shí)際巖石樣品進(jìn)行掃描，得到實(shí)際巖石樣品的數(shù)字巖忍圖像并對 圖像進(jìn)行裁剪;將數(shù)字巖忍圖像二值化，得到數(shù)字巖忍圖像的骨架和孔隙； (2) 確定數(shù)字巖忍圖像的表征單元體積 1) 從1%巖忍體積開始對數(shù)字巖忍圖像隨機(jī)選取至少10個1%巖忍體積的子體積，計(jì)算 每個子體積的孔隙度并記錄;然后W2%巖忍體積的速率逐步同步擴(kuò)大每個子體積并記錄 新得到子體積的孔隙度，直到所有子體積的孔隙度相對大小都大于95%為止，記錄此時子 體積的體積VI;其中，子體積的孔隙度計(jì)算公式為：4為孔隙度， sumpore為子體積孔隙像素點(diǎn)總數(shù)，sumgrain為子體積骨架像素點(diǎn)總數(shù)； 2) 從1%巖忍體積開始對數(shù)字巖忍圖像隨機(jī)選取至少10個1%巖忍體積的子體積，計(jì)算 每個子體積的分形維數(shù)并記錄;然后W2%巖忍體積的速率逐步同步擴(kuò)大每個子體積并記 錄新得到子體積的分形維數(shù)，直到所有子體積的分形維數(shù)相對大小都大于95%為止，記錄 此時子體積的體積V2;其中，子體積的分形維數(shù)利用計(jì)盒算法，使用基于3D圖像的分形盒維 數(shù)算法得出； 3) 從1%巖忍體積開始對數(shù)字巖忍圖像隨機(jī)選取至少10個1%巖忍體積的子體積，計(jì)算 每個子體積的滲透率張量并記錄;然后W2%巖忍體積的速率逐步同步擴(kuò)大每個子體積并 記錄新得到子體積的滲透率張量，直到所有子體積的滲透率張量其矩陣對應(yīng)值相對大小都 大于90%為止，記錄此時子體積的體積V3;其中，子體積的滲透率張量基于廣義達(dá)西定律和 質(zhì)量守恒方程，利用有限元計(jì)算得到； 4) 比較Vl，V2和V3，其中的最大體積就是表征單元體積； (3) 滲流模擬 根據(jù)步驟4)確定的表征單元體積從數(shù)字巖忍圖像中隨機(jī)選取表征單元體積，利用滲透 率計(jì)算模型對表征單元體積進(jìn)行滲流過程數(shù)值模擬，得到表征單元體積的絕對滲透率和相 對滲透率，絕對滲透率值和相對滲透率值即為實(shí)際巖石滲透率的預(yù)測值。2. 根據(jù)權(quán)利要求1所述的滲透率預(yù)測方法，其特征在于，所述的滲流過程數(shù)值模擬方法 為孔隙網(wǎng)絡(luò)模型或CFD模型。
【專利摘要】本發(fā)明屬于石油、地質(zhì)科學(xué)技術(shù)領(lǐng)域，提供了一種基于CT圖像的滲透率預(yù)測方法。CT掃描：利用工業(yè)X射線CT對實(shí)際巖石樣品進(jìn)行掃描，得到實(shí)際巖石樣品的數(shù)字巖芯圖像并對圖像進(jìn)行裁剪；將數(shù)字巖芯圖像二值化，得到數(shù)字巖芯圖像的骨架和孔隙；確定數(shù)字巖芯圖像的表征單元體積；滲流模擬：絕對滲透率值和相對滲透率值即為實(shí)際巖石滲透率的預(yù)測值。本發(fā)明的滲透率預(yù)測方法過程簡單，很容易實(shí)現(xiàn)不同地質(zhì)條件下巖石滲透率的多次測量。測量方法可對具有各向異性的巖石的滲透率進(jìn)行準(zhǔn)確預(yù)測。測量結(jié)果可以直接作為數(shù)值模擬的輸入?yún)?shù)，測量方法適用于任何氣‐液‐固或液‐液‐固三相或多相體系的巖石滲透率預(yù)測。本發(fā)明也可推廣到巖石其它物性的預(yù)測。
【IPC分類】G01N15/08
【公開號】CN105510204
【申請?zhí)枴緾N201510880934
【發(fā)明人】宋永臣, 呂鵬飛, 劉瑜, 蔣蘭蘭, 滕瑩, 武博浩, 陸國歡, 王夢婷
【申請人】大連理工大學(xué)
【公開日】2016年4月20日
【申請日】2015年12月3日