(12)發(fā)明專利道8號審查員梁策(72)發(fā)明人張啟輝李海濤王文星王科陸寬馬欣高素娟聶松馬寒松代晶晶一種受限制游走模擬以計算多孔介質內流本發(fā)明公開一種受限制游走模擬以計算多度值表征的微觀尺度連續(xù)相粒子在不同區(qū)域不分辨率的連續(xù)相算子運移方向以及數(shù)字巖心鏡法獲取的迂曲度數(shù)據(jù)，相較于現(xiàn)有的koponen、21.一種受限制游走模擬以計算多孔介質內流動迂曲度的方法，其特點在于，包括以下Step2:對處理后的巖樣三維灰度Step4:通過不斷去除巖石空隙相邊界體素，結合現(xiàn)場實驗將Step6:通過隨機游走步數(shù)與距離的關系表征樣品儲層迂曲度，迂曲度通過平均運移距2.根據(jù)權利要求1所述的一種受限制游走模擬以計算多孔介質內流動迂曲度的方法，3.根據(jù)權利要求2所述的一種受限制游走模擬以計算多孔介質內流動迂曲度的方法，Step202:基于修正后的三維數(shù)字巖心數(shù)據(jù)，通過設定灰度或密度閾值劃分骨架與孔隙，遍歷劃分后的二值化孔隙，對矩陣內標記為孔隙的坐標進行bwlabel算法計算其獨立整閾值后統(tǒng)計孔隙度φ2以及不連通孔隙φ,通過不斷調整閾值使得φ?=92-9,從而確4.根據(jù)權利要求3所述的一種受限制游走模擬以計算多孔介質內流動迂曲度的方法，灰度值進行高斯平滑操作，抹除極小值后調整分水嶺算法中孔隙灰度閾值，根據(jù)數(shù)35.根據(jù)權利要求4所述的一種受限制游走模擬以計算多孔介質內流動迂曲度的方法，當算子在隨機游走過程中運移到數(shù)字巖心矩陣中骨架位置時，則算子位置保持不動，Step502:經過Step501過程獲取數(shù)字巖心矩陣中獲取的不同方向概率，調用編程軟件Step503:使用編程軟件重復迭代Step501-Step置。6.根據(jù)權利要求5所述的一種受限制游走模擬以計算多孔介質內流動迂曲度的方法，其特點在于，所述Step503過程中根據(jù)不同需要選取不同步長以獲取不同算子的具體位置4一種受限制游走模擬以計算多孔介質內流動迂曲度的方法技術領域[0001]本發(fā)明涉及氣藏開發(fā)技術領域，特別涉及非常規(guī)油氣開發(fā)領域，具體是一種受限制游走模擬以計算多孔介質內流動迂曲度的方法。背景技術[0002]作為描述滲流通道的一個重要參數(shù)，孔道迂曲度定義為滲流過程中指定運移的實際長度與滲流通道的宏觀長度的比值。隨著近年來對非常規(guī)油氣藏的開發(fā)，油氣滲流過程愈加復雜，常規(guī)的通過孔隙結構進行簡單計算獲取的迂曲度方法所獲得的迂曲度用以解釋滲流過程存在較大誤差。[0003]近年來專家及學者們趨向于使用高精度CT反演制成數(shù)字巖心以模擬復雜非常規(guī)油氣藏中的運移，但對迂曲度的解釋則僅僅停留在孔介質礦物的不同平均粒徑量化上，忽略了真實滲流過程中粒子間作用力以及儲層巖石的各向異性影響。發(fā)明內容[0004]針對上述問題，本發(fā)明旨在提供一種技術以精確的表征復雜多孔介質的迂曲度方法，使其結果上更接近于真實滲流情況。[0005]本發(fā)明的技術方案如下：[0006]一種受限制游走模擬以計算多孔介質內流動迂曲度的方法，其特點在于，包括以[0007]Step1:對巖樣進行360度掃描，基于巖樣多角度圖像進行三維修正處理；[0008]Step2:對處理后的巖樣三維灰度數(shù)據(jù)通過圖像閾值分割分離巖樣孔隙與巖樣骨架。對比矩陣統(tǒng)計孔隙與實際測試孔隙度調整閾值數(shù)據(jù)，獲取表征真實巖樣骨架與孔隙的數(shù)據(jù)；[0009]Step3:結合現(xiàn)場實驗，使用分水嶺算法分割空隙相歐式距離矩陣，將復雜矩陣轉化為較為簡單的對應樣品空間連接點集合，設定該矩陣區(qū)域為隨機游走中無碰撞區(qū)域。[0010]Step4:通過不斷去除巖石空隙相邊界體素，結合現(xiàn)場實驗將對應獲取對應巖樣在對應掃描分辨率下單個體素的碰撞與反彈概率。[0011]Step5:選取樣品數(shù)字巖心矩陣中孔隙內任意一點，進行隨機游走模擬，游走過程中記錄對應步數(shù)對應游走距離，對游走過程中矩陣邊界進行鏡像映射。隨機游走模擬坐標在接觸數(shù)字巖心骨架采取不同措施限制步數(shù)或游走距離。[0012]Step6:通過隨機游走步數(shù)與距離的關系表征樣品儲層迂曲度。[0014]使用計算機遍歷重構后的三維數(shù)字巖心矩陣，以灰度數(shù)據(jù)量為標準，將灰度數(shù)據(jù)劃分為兩個頻段，對頻段劃分后的數(shù)據(jù)，在較難區(qū)分的部分基于對應灰度值與密度的重構5[0016]Step201:通過液體飽和排液法或氦氣法等巖樣孔隙度實驗測定方法測定樣品孔整閾值后統(tǒng)計孔隙度φ?以及不連通孔隙φ。通過不斷調整閾值使得φ?=9?-φ,從而確[0019]對數(shù)字巖心中根據(jù)Step2過程中劃分為孔隙的區(qū)域中灰度值進行高斯平滑操作，對應[x,y,z]位置矩陣中。及Step3過程中分水嶺算法分隔的獨立孔隙，則巖樣流動過程中各體素分布概率可以通過6[0037]進一步地，所述Step503過程中根據(jù)不同需要選取不同步長以獲取不同附圖說明[0040]圖1為本發(fā)明一種受限制游走模擬以計算多孔介質內流動迂曲度的方法流程示意[0042]圖3為實施例1計算得到的基于分水嶺算法計算處理后的矩陣數(shù)據(jù)(Avizo可視化);7具體實施方式[0044]下面結合附圖和實施例對本發(fā)明進一步說明。需要說明的是，在不沖突的情況下，本申請中的實施例及實施例中的技術特征可以相互結合。需要指出的是，除非另有指明，本申請使用的所有技術和科學術語具有與本申請所屬技術領域的普通技術人員通常理解的者物件涵蓋出現(xiàn)在該詞后面列舉的元件或者物件及其等同，而不排除其他元件或者物件。[0046]如附圖1所示，一種基于隨機游走精確測定多孔介質內滲流動態(tài)迂曲度的方法，包括以下步驟：[0047]Stepl:對巖樣進行360度掃描，基于巖樣多角度圖像進行三維修正處理；轉化后的三維數(shù)字巖心矩陣如表1所示，經Avizo可視化讀取的數(shù)據(jù)如說明書附圖2所示。[0048]表1轉化后的三維數(shù)字巖心數(shù)據(jù)矩陣12345789……1234123456985/5頁5/5頁[0051]Step2:對處理后的巖樣三維灰度數(shù)據(jù)通過圖像閾值分割分離巖樣孔隙與巖樣骨架。對比矩陣統(tǒng)計孔隙與實際測試孔隙度調整閾值數(shù)據(jù)，獲取表征真實巖樣骨架與孔隙的[0052]經過調整閾值，以閾值設定4817獲得的孔隙度0.143與試驗結果最相近，則以該閾值劃分孔隙，圖2為處理后的矩陣數(shù)據(jù)經Avizo可視化后的結果。[0053]Step3:結合現(xiàn)場實驗，使用分水嶺算法分割空隙相歐式距離矩陣，將復雜矩陣轉化為較為簡單的對應樣品空間連接點集合，設定該矩陣區(qū)域為隨機游走中無碰撞區(qū)域。[0054]經分水嶺算法及修正后的隨機游走無碰撞區(qū)域如圖3三種藍色位置所示。[0055]Step4:通過不斷去除巖石空隙相邊界體素，結合現(xiàn)場實驗將對應獲取對應巖樣在對應掃描分辨率下單個體素的碰撞與反彈概率。[0056]以數(shù)字巖心矩陣坐標[77,77,77]為例，基于其灰度值計算其在9個方向的運移概率為163/814,5/906,16/103,49/370,5/56,25/222,24/299,50/841,35/213。[0057]Step5:選取樣品數(shù)字巖心矩陣中孔隙內任意一點，進行隨機游走模擬，游走過程中記錄對應步數(shù)對應游走距離，對游走過程中矩陣邊界進行鏡像映射。隨機游走模擬坐標在接觸數(shù)字巖心骨架采取不同措施限制步數(shù)或游走距離。[0058]數(shù)字巖心樣品中存在不同粒徑，將數(shù)字巖心樣品分為大粒徑，小粒徑總體以及按照無制約空間分別進行10000次*100000步游走，記錄不同情況下運移距離。[0059]Step6:通過隨機游走步數(shù)與距離的關系表征樣品儲層迂曲度。[0060]如圖4所示，游走步數(shù)高于70000次時，四種不同情況下的擴散系數(shù)均趨于定值，因此，該樣品中大粒徑區(qū)域，小粒徑區(qū)域以及整體迂曲度分別為，1.7,2.1以及1.9.將此值帶入泊肅葉方程進行流動模擬，模擬結果與真實結果相差8%,相比孔隙度計算的迂曲度值，模擬誤差降低了11%。[0061]以上所述，僅是本發(fā)明的較佳實施例而已，并非對本發(fā)明作任何形式上的限制，雖然本發(fā)明已以較佳實施例揭露如上，然而并非用以限定本發(fā)明，任何熟悉本專業(yè)的技術人員，在不脫離本發(fā)明技術方案范圍內，當可利用上述揭示的技術內容作出些許更動或修飾為等同變化的等效實施例，但凡是未脫離本發(fā)明技術方案的內容，依據(jù)本發(fā)明的技術實質對以上實施例所作的任何簡單修改、等同變化與修飾，均仍屬于本發(fā)明技術方案的范圍內。9樣品CT測試轉化CT圖像為三維數(shù)字巖心并修正讀取數(shù)字巖心數(shù)據(jù)，通過分水嶺算法劃分矩陣內區(qū)域計算矩陣內