數(shù)字圖像重建技術(shù)在巖石滲流類課程教學(xué)中的應(yīng)用

宋睿+崔夢夢++鄭玲靜

摘 要：提升高校巖石滲流類課程的教學(xué)效果對保障國家建設(shè)和能源開采的人才供應(yīng)具有重要意義。針對傳統(tǒng)教學(xué)中存在的不足，本文提出了一種結(jié)合巖石微觀數(shù)字圖像及其重建模型的新型課堂教學(xué)方法。以某砂巖微CT圖像為例，本文詳細(xì)介紹了在巖石多孔介質(zhì)特性、單相-兩相流體輸運(yùn)特征教學(xué)中的課程案例，詳細(xì)論述了該方法在各關(guān)鍵知識點講解中的應(yīng)用。分析表明，本文提出的授課方法有助于改善空洞、乏味的授課方式，從而提升了該類課程的教學(xué)效果。

關(guān)鍵詞：巖石 數(shù)字圖像 重建模型

中圖分類號：TP39 文獻(xiàn)標(biāo)識碼：A 文章編號：1674-098X（2017）10（c）-0238-04

Abstract： It is of significant importance to ensure the demands of technical talents for national construction and energy extraction of China by promoting the teaching efficiency of high education courses related to seepage mechanism in rock and soil. Regarding to the existing shortage of traditional teaching method， this paper presents an improved teaching case using the digital images and reconstructed models of rock. The detailed teaching cases including the nonhomogeneous porous media characteristics of rock and single-/two- phase flow are designed using the sandstone micro-CT images. Then the key knowledge points are analysed. It is indicated that the proposed teaching methods have been proved to enrich the teaching mode and improve the teaching efficiency.

Key Words： Courses related to seepage； Digital images； Reconstructed rock models

隨著我國經(jīng)濟(jì)建設(shè)的高速發(fā)展，大規(guī)模的土木工程建設(shè)（如建筑、交通、水利水電等）、資源礦產(chǎn)開采（如煤炭、石油、天然氣及地?zé)岬荣Y源的開采）、地下儲庫建設(shè)（如地下儲氣庫，核廢料、化學(xué)廢棄物及CO2地質(zhì)封存）等一系列國計民生工程急需大批有相關(guān)專業(yè)知識背景的人才[1-2]。作為這些大型工程的共同施工對象，巖石類材料是高校相關(guān)工程專業(yè)的重點研究對象。巖石是由不同種類的固體礦物組成的骨架和由骨架分隔成大量密集成群的微小孔隙構(gòu)成，孔隙中往往存在著液相、氣相系統(tǒng)，是一種典型的非均質(zhì)多孔介質(zhì)材料。相關(guān)專業(yè)的高等教育課程開設(shè)有巖土力學(xué)（屬于固體力學(xué)）及滲流力學(xué)（屬于流體力學(xué)）兩大類課程，包括土力學(xué)、巖石力學(xué)、滲流力學(xué)、水文地質(zhì)、水力學(xué)、地下水動力學(xué)等專業(yè)基礎(chǔ)理論課程，涵蓋了工程類高等教育的教學(xué)難點和重點[3]。因此，提升巖石相關(guān)課程的教學(xué)效果，對保障我國經(jīng)濟(jì)建設(shè)和能源開采的人才供應(yīng)具有重要的意義。

然而，此類課程公式較多，推導(dǎo)過程復(fù)雜，課堂計算練習(xí)多局限于一維或規(guī)則二維模型[4]。由于巖石孔隙尺寸狹小，傳統(tǒng)教學(xué)中常通過展示示意圖或巖石薄片的形式向?qū)W生講解巖石的多孔介質(zhì)特征，這類二維圖形往往難以讓學(xué)生對巖石空間特征產(chǎn)生直觀、立體的印象，甚至產(chǎn)生理解誤區(qū)。如部分學(xué)生看到巖石切片會誤認(rèn)為巖石孔隙并未連通，是孤立的，在看到連通的二維孔隙時又誤以為孔隙如同連接的水管一樣等。此外，作為一種不透明、孔隙狹小、嚴(yán)重非均質(zhì)的多孔介質(zhì)，巖石的微觀孔隙特征及其中流體的輸運(yùn)過程往往難以觀測，限制了高校滲流類課程的教學(xué)效果。

近年來，隨著以微觀-CT為代表的巖石微成像技術(shù)的發(fā)展，研究人員可獲取微納米級別的三維巖石微觀圖像，這為相關(guān)課程的教學(xué)提供了新的手段。本文以巖石微觀CT數(shù)字圖像及其重建模型為例，合理利用新興研究成果，與實際緊密結(jié)合，避免空洞的課堂教學(xué)方式，以達(dá)到提升巖石滲流類課程教學(xué)效果的目的。

1 巖石數(shù)字圖像及模型重建技術(shù)簡介

巖石數(shù)字圖像技術(shù)是指通過一系列先進(jìn)成像手段，例如光學(xué)顯微鏡-攝像機(jī)成像、電子掃描顯微鏡、原子力顯微鏡等二維成像手段，或共聚焦離子顯微鏡、微觀CT等三維成像手段，獲得微-納米級別的巖石微觀圖像。通過對巖石微觀圖像的觀察，可以直觀獲得巖石固體顆粒尺寸、形態(tài)，礦物類型及分布特征，和孔隙尺寸、形態(tài)、連結(jié)性等。特別是三維巖石數(shù)字圖像，更能加深初學(xué)者對巖石多孔介質(zhì)特性和非均質(zhì)性的感性認(rèn)識和空間想象力。同時，利用對巖石骨架和孔隙面積/體積的統(tǒng)計計算，可以進(jìn)一步獲取巖石孔隙度、粒徑分布和孔徑分布數(shù)據(jù)。

巖石模型重建是以巖石的數(shù)字化圖像為基礎(chǔ)，通過降噪、閾值分割等圖像處理技術(shù)，結(jié)合模型重建算法將數(shù)字圖像轉(zhuǎn)換為拓?fù)浣Y(jié)構(gòu)相似或相同的、可用于數(shù)值模擬研究的計算模型，常見的有網(wǎng)格模型和等效孔隙網(wǎng)絡(luò)模型兩類。在滲流力學(xué)課程方面，可以開展單相、兩相及三相滲流數(shù)值模擬，除了獲取傳統(tǒng)實驗得到的壓力-出口流量曲線、相對滲透率曲線，還可清晰觀測到流體在孔隙內(nèi)的主要流動通道、兩/三相驅(qū)替過程、殘余液體（如殘余油、污染物等）分布位置及形成過程，加深學(xué)生對這些抽象機(jī)理現(xiàn)象的認(rèn)知程度，提升教學(xué)效果。endprint

2 課程設(shè)計

本節(jié)將以教學(xué)設(shè)計實例具體分析巖石數(shù)字圖像和重建模型在相關(guān)課程中的應(yīng)用，及其可行性、優(yōu)越性。

2.1 巖石孔隙結(jié)構(gòu)及礦物分布三維觀測

在巖石類課程中，巖石所具有的多孔介質(zhì)特性是其他相關(guān)課程的基礎(chǔ)。結(jié)合巖石三維數(shù)字圖像（圖1，分辨率2.51um/像素，980×1006×960像素），可以更為有效地開展教學(xué)。通過展示圖1（a）和圖1（b），能直觀看到巖石礦物類型的非均質(zhì)性及空間分布，在圖1（b）中表現(xiàn)為圖像灰度值的差異；通過展示圖1（c），可以加強(qiáng)學(xué)生對巖石骨架微結(jié)構(gòu)特征的認(rèn)識，為后續(xù)巖石力學(xué)和破裂力學(xué)的學(xué)習(xí)奠定基礎(chǔ)；通過展示圖1（d），可以開展巖石孔隙特征教學(xué)，加深學(xué)生對巖石孔隙形狀、連通性、空間分布的認(rèn)識；通過局部放大圖像1（d），可展示巖石有效孔隙、盲端孔隙、孤立孔隙的分布及結(jié)構(gòu)特征；通過對比巖石骨架像素體積和孔隙像素體積，可開展巖石孔隙度求解，強(qiáng)化學(xué)生對孔隙度的理解和認(rèn)識。

2.2 巖石單相滲流三維可視化數(shù)值實驗

巖石孔隙流體的輸運(yùn)特征是滲流類課程的重難點，本節(jié)給出了巖石單相滲流三維可視化數(shù)值實驗的課程設(shè)計，具體流程為：以砂巖S1為例，在巖石三維圖像的基礎(chǔ)上，利用MIMICS軟件對重建得到的巖石孔隙幾何模型進(jìn)行網(wǎng)格劃分，得到如圖2（a）所示的巖樣網(wǎng)格模型；將網(wǎng)格導(dǎo)入CFD模擬軟件Fluent，流體輸運(yùn)方程采用NS方程組，沿圖中所示Z方向施加進(jìn)出口壓力，以模擬巖石孔隙中的水在壓力梯度驅(qū)動下的滲流過程。模擬得到水在孔隙中的速度場及壓力場分布如圖2b所示，從速度場和壓力場云圖中能夠讓學(xué)生直觀地觀察到由于孔隙結(jié)構(gòu)的非均質(zhì)性，流體主要沿著尺寸較大、連通性較好的孔隙向前推進(jìn)（以減小流動過程中的阻力），從而使學(xué)生加深孔隙中流體的主要流動通道這一概念；同時利用模擬軟件得到的模型出口水流量，結(jié)合達(dá)西滲流定律可以計算出模型的滲透率，從而加強(qiáng)學(xué)生對達(dá)西定律的理解與認(rèn)識；通過改變進(jìn)出口壓力梯度，再次重復(fù)計算模型的滲透率，來驗證滲透率是巖石固有物性參數(shù)而并不會隨出口壓力不同而改變，從而強(qiáng)化該知識點的教學(xué)效果；通過改變進(jìn)出口邊界條件，如將進(jìn)出口壓力轉(zhuǎn)換為速度邊界條件，強(qiáng)化學(xué)生對流體靜壓和動壓的概念及其相互轉(zhuǎn)換。

2.3 巖石兩相滲流三維可視化數(shù)值實驗

兩相滲流主要以油水兩相滲流過程為例，具體流程為：以S1為例，將孔隙的網(wǎng)格模型導(dǎo)入Fluent軟件，選擇軟件中多相流模型中的VOF（volume of fraction）模型，模型初始假定為飽和油狀態(tài)，沿圖中所示Z方向施加進(jìn)出口壓力，以模擬巖石孔隙中水驅(qū)替油的過程。模擬得到不同計算時間步模型油水飽和度分布云圖如圖3所示，在課程教學(xué)中將多幅圖像組合為動圖的形式在教學(xué)中展示，可清晰觀察到兩相驅(qū)替過程中的流動特征：水沿阻力最小的流動通道驅(qū)替油相，呈現(xiàn)出不均勻推進(jìn)的油水交界面，該現(xiàn)象稱為指進(jìn)現(xiàn)象。在本模型中表現(xiàn)為注入水從進(jìn)口沿某幾條主要孔喉鏈率先突破至出口處；在實際地層中具體表現(xiàn)為注入水沿注水井進(jìn)入流動阻力最小的通道到達(dá)采油井，是開發(fā)過程中影響采收率的主要因素。結(jié)合該云圖還可以對殘余油的主要形成與分布規(guī)律進(jìn)行分析，例如：因注入水未波及而滯留于模型邊角處的殘余油，指進(jìn)現(xiàn)象導(dǎo)致的水率先突破而造成的卡斷殘余油，殘留于孔隙死角或盲端的殘余油。通過軟件輸出的不同時間步長的出口油水兩相流量，結(jié)合達(dá)西定律和相對滲透率公式，可繪制出該過程的油水相對滲透率曲線，從而強(qiáng)化學(xué)生對油水相對滲透率曲線的物理意義和計算過程的理解和記憶。同時，通過改變模型進(jìn)出口邊界、物性參數(shù)、潤濕性等可模擬不同條件下的開發(fā)效果。例如，改變進(jìn)出口壓力以模擬不同的注水壓力或速度對注水開發(fā)效果的影響，改變注入水物性參數(shù)以模擬不同毛細(xì)管數(shù)條件下的開發(fā)效果，改變模型的潤濕性以模擬親水、親油、中間潤濕性和混合潤濕性巖石中的油水流動規(guī)律。以類比的方法，用動圖的形式強(qiáng)化學(xué)生對相關(guān)概念、公式和物理參數(shù)的理解和認(rèn)識，從而達(dá)到提升教學(xué)效果的目的。

3 結(jié)語

本文利用巖石微觀CT成像這一新興技術(shù)，結(jié)合巖石數(shù)字圖像及其重建模型，詳細(xì)論證了該課程資源在講解巖石非均質(zhì)多孔介質(zhì)特性的優(yōu)越性；并結(jié)合CFD數(shù)值模擬軟件，論述了重建模型在巖石單相滲流和兩相滲流教學(xué)中的應(yīng)用。結(jié)果表明，巖石數(shù)字圖像及其重建模型能夠加深學(xué)生對巖石孔隙度、孔徑分布、礦物組分分布、固體顆粒和孔隙形狀、尺寸等材料參數(shù)的認(rèn)識，強(qiáng)化其對巖石孔隙的空間分布及連通性、單相-兩相流體輸運(yùn)特征的理解，彌補(bǔ)了傳統(tǒng)教學(xué)過程中空洞說教的不足，從而提升了課程的教學(xué)效果。同時，該課程案例對其他多孔介質(zhì)類教學(xué)也具有較強(qiáng)的借鑒意義。

參考文獻(xiàn)

[1] 姚仰平，羅汀，朱建明，等.創(chuàng)建巖土工程學(xué)科特色課程體系培養(yǎng)創(chuàng)新人才[J].北京航空航天大學(xué)學(xué)報：社會科學(xué)版，2009，22（s1）：62-65.

[2] 林彤，唐輝明，程祖依.巖土工程專業(yè)人才產(chǎn)學(xué)合作教育研究[J].工程地質(zhì)學(xué)報，2000（S1）：628-630.

[3] 王驍.數(shù)字巖心分析資料數(shù)據(jù)庫的應(yīng)用研究[J].科技創(chuàng)新導(dǎo)報，2014（2）：23.

[4] 張婭，張亞軍，王維民，等.理論力學(xué)教學(xué)質(zhì)量提高措施探索[J].科技創(chuàng)新導(dǎo)報，2017，14（1）：167-168.endprint