王璐,楊勝來*,劉義成,王云鵬,孟展,韓偉,錢坤
1 中國(guó)石油大學(xué)(北京)油氣資源與工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 102249
2 中國(guó)石油西南油氣田分公司勘探開發(fā)研究院,成都 610041
*通信作者, yangsl@cup.edu.cn
石油工程
縫洞型碳酸鹽巖儲(chǔ)層氣水兩相微觀滲流機(jī)理可視化實(shí)驗(yàn)研究
王璐1,楊勝來1*,劉義成2,王云鵬1,孟展1,韓偉1,錢坤1
1 中國(guó)石油大學(xué)(北京)油氣資源與工程國(guó)家重點(diǎn)實(shí)驗(yàn)室,北京 102249
2 中國(guó)石油西南油氣田分公司勘探開發(fā)研究院,成都 610041
*通信作者, yangsl@cup.edu.cn
四川盆地縫洞型碳酸鹽巖氣藏多屬于有水氣藏,氣水兩相滲流機(jī)理復(fù)雜。目前氣水兩相微觀滲流可視化模型多是基于理想孔隙結(jié)構(gòu)或鑄體薄片圖像制作,無法還原儲(chǔ)層中真實(shí)孔、縫、洞分布,且機(jī)理研究多集中于孔隙型和裂縫型儲(chǔ)層,缺少對(duì)孔洞型和縫洞型儲(chǔ)層的認(rèn)識(shí)。通過將巖心CT掃描與激光刻蝕技術(shù)結(jié)合,以四川盆地震旦系儲(chǔ)層中裂縫型、孔洞型和縫洞型碳酸鹽巖巖心CT掃描結(jié)果為模板,設(shè)計(jì)并研制了3類巖心的可視化模型,據(jù)此研究了氣水兩相微觀滲流機(jī)理及封閉氣、殘余水形成機(jī)理,通過ImageJ灰度分析法實(shí)現(xiàn)了氣水分布的定量表征。研究結(jié)果表明:3類模型水驅(qū)氣和氣驅(qū)水過程中的滲流規(guī)律各不相同;繞流、卡斷、盲端和角隅處形成的封閉氣普遍存在于3類模型中,此外還在“H型”孔道處、“啞鈴型”通道處和微裂縫縫網(wǎng)處形成特殊封閉氣;孔道、裂縫壁面上的束縛水膜,溶洞中部的圓潤(rùn)水團(tuán),狹窄喉道處的卡斷水柱和狹長(zhǎng)孔道處的滯留水柱是殘余水的主要形式;裂縫型模型水竄最嚴(yán)重,無水采收期最短,采出程度最低,而孔洞型模型水驅(qū)前緣推進(jìn)均勻,無水采收期最長(zhǎng),采出程度也最高。該研究實(shí)現(xiàn)了對(duì)縫洞型碳酸鹽巖不同類型儲(chǔ)層氣水兩相微觀滲流規(guī)律的精準(zhǔn)刻畫,為類似氣藏的高效開發(fā)提供了理論依據(jù)。
縫洞型碳酸鹽巖;氣水兩相滲流;可視化;封閉氣;殘余水;灰度分析
縫洞型碳酸鹽巖氣藏是經(jīng)過多期構(gòu)造運(yùn)動(dòng)與古巖溶共同作用形成的一種特殊類型氣藏,其儲(chǔ)集介質(zhì)由溶洞、溶孔和裂縫組成,具有構(gòu)造復(fù)雜、儲(chǔ)層非均質(zhì)性強(qiáng)、孔洞縫宏觀發(fā)育及氣水兩相滲流規(guī)律復(fù)雜等特征,是當(dāng)前最復(fù)雜的特殊氣藏之一[1-4]。2013年四川盆地震旦系-寒武系特大型氣田被發(fā)現(xiàn),其中發(fā)育著震旦系燈影組碳酸鹽巖縫洞型、寒武系龍王廟組白云巖孔隙型2套主要含氣儲(chǔ)集層,該類復(fù)雜氣藏開始引起廣大科研工作者的關(guān)注[5]。
長(zhǎng)期以來,國(guó)內(nèi)外氣藏的研究主要集中在開采工藝技術(shù)和數(shù)值模擬方面,對(duì)于氣水兩相微觀滲流機(jī)理研究較少。常規(guī)的氣水兩相滲流數(shù)學(xué)模型并不能真實(shí)反映巖石的孔隙結(jié)構(gòu)分布,也不能反映巖石中的氣水兩相真實(shí)滲流特征。2002年周克明等[6]首次提出以巖心樣品的鑄體薄片所代表的孔隙結(jié)構(gòu)為背景,借助激光刻蝕技術(shù),研制了氣水兩相可視化物理模型,并通過該模型研究了水驅(qū)氣機(jī)理及封閉氣的形成方式,實(shí)現(xiàn)了氣水兩相滲流過程的可視化。隨后,可視化模型被廣泛應(yīng)用到氣水兩相滲流機(jī)理的研究上[7-14]。但這些研究仍然存在許多問題:(1)可視化模型都是基于理想孔隙結(jié)構(gòu)或鑄體薄片制作的,無法還原巖心樣品中的真實(shí)孔、縫、洞分布,得到的實(shí)驗(yàn)結(jié)果與實(shí)際情況存在差異。(2)研究主要集中在孔隙型和裂縫型氣藏,而對(duì)于更為復(fù)雜的孔洞型和縫洞型氣藏尚未進(jìn)行研究。(3)氣水兩相可視化圖像結(jié)果取自于水驅(qū)氣實(shí)驗(yàn)完成后,缺少氣水兩相滲流過程中的可視化圖像與機(jī)理分析。(4)氣水兩相滲流實(shí)驗(yàn)以模擬氣藏水侵過程的水驅(qū)氣為主,而對(duì)于氣驅(qū)水過程的研究與分析較少。
針對(duì)以上問題,筆者首次將巖心CT掃描技術(shù)與激光刻蝕技術(shù)相結(jié)合,以四川盆地震旦系儲(chǔ)層中的裂縫型、孔洞型和縫洞型碳酸鹽巖巖心CT掃描圖像為模板,設(shè)計(jì)并研制了3類巖心的激光刻蝕透明仿真模型,實(shí)現(xiàn)了對(duì)實(shí)際儲(chǔ)層中孔、縫、洞分布的精準(zhǔn)刻畫。通過開展氣水兩相微觀滲流實(shí)驗(yàn),研究了3類模型的水驅(qū)氣、氣驅(qū)水微觀機(jī)理,封閉氣、殘余水形成機(jī)理以及采出方法,并利用ImageJ灰度分析法實(shí)現(xiàn)了對(duì)氣水兩相微觀分布的定量表征,形成了一套較為完善的縫洞型碳酸鹽巖儲(chǔ)層氣水兩相微觀滲流理論,為該類氣藏的高效開發(fā)提供了有效依據(jù)。
油氣儲(chǔ)層的微觀孔隙結(jié)構(gòu)特征與分布規(guī)律決定了儲(chǔ)層中流體的微觀分布關(guān)系、滲流機(jī)理與滲流規(guī)律[15-16]。為了完全模擬實(shí)際儲(chǔ)層中孔、縫、洞的結(jié)構(gòu)特征和分布規(guī)律,利用微米CT對(duì)震旦系儲(chǔ)層中具有代表性的裂縫型、孔洞型和縫洞型3塊巖心樣品進(jìn)行掃描,并對(duì)6 000余張CT掃描圖像進(jìn)行篩選和提取后,借助現(xiàn)代激光刻蝕技術(shù),制作了儲(chǔ)層結(jié)構(gòu)特征明顯的可視化微觀模型(見圖1)。在進(jìn)行CT掃描圖像抽提時(shí)發(fā)現(xiàn),大部分圖像的孔、縫、洞分布規(guī)律并不明顯,在單張圖像中往往只存在一種典型結(jié)構(gòu)特征,并不能完全代表該類型儲(chǔ)層,而參照這些圖像制作的可視化模型結(jié)構(gòu)會(huì)過于簡(jiǎn)單,既無法得到氣水兩相滲流的全部流動(dòng)特征,也無法對(duì)微觀滲流機(jī)理、封閉氣和殘余水的形成機(jī)理進(jìn)行系統(tǒng)全面地分析。為此,通過對(duì)比最終選取了結(jié)構(gòu)特征明顯,孔、縫、洞分布規(guī)律復(fù)雜的3類圖像,并以這3類圖像作為母版,在不破壞原始結(jié)構(gòu)的基礎(chǔ)上增加了一些母版缺少的結(jié)構(gòu),以期通過這些母版制作的可視化模型能夠較為全面地代表儲(chǔ)層的各種典型結(jié)構(gòu),得到更加系統(tǒng)全面地機(jī)理分析與流動(dòng)特征。其中,孔洞型模型與縫洞型模型在結(jié)構(gòu)特征與研究目的上存在本質(zhì)差別。在結(jié)構(gòu)特征上,孔洞型模型以孔隙和溶洞這2類空間為主,孔隙作為滲流通道,設(shè)計(jì)尺寸較小,毛管力較大且分布廣泛;而縫洞型模型以裂縫和溶洞為主,裂縫作為滲流通道,設(shè)計(jì)尺寸較大,毛管力較小且分布離散。在研究目的上,孔洞型模型主要研究孔隙與溶洞之間的兩相流動(dòng)規(guī)律,縫洞型模型的重點(diǎn)則在裂縫與溶洞之間。由于裂縫的導(dǎo)流能力和通道尺寸遠(yuǎn)遠(yuǎn)大于孔隙,所以在流動(dòng)特征與封閉氣、殘余水形成機(jī)理上存在較大差異。
設(shè)計(jì)模型尺寸為8.50 cm×8.50 cm,有效尺寸為6.00 cm×6.00 cm,孔隙直徑為0.099 mm~0.181 mm,裂縫直徑為0.325 mm~0.492 mm,溶洞直徑為1.161 mm~1.657 mm,孔、縫、洞尺寸比例均按照CT掃描結(jié)果設(shè)計(jì)。通過接觸角儀對(duì)模型材料與水之間的接觸角進(jìn)行了測(cè)量,結(jié)果表明接觸角為36.4°,原始潤(rùn)濕性為水濕。微觀可視化模型制作流程如下:①依據(jù)CT掃描微觀結(jié)構(gòu)圖像,利用Auto CAD軟件繪制出可視化模型藍(lán)圖。②根據(jù)模型藍(lán)圖通過數(shù)控銑床加工出掩膜版。③利用紫外光通過掩膜版照射到附有一層光刻膠薄膜的基片表面,使曝光區(qū)域的光刻膠發(fā)生化學(xué)反應(yīng)。④通過顯影技術(shù)溶解去除曝光區(qū)域或未曝光區(qū)域的光刻膠,使掩膜版上的圖形復(fù)制到光刻膠薄膜上。⑤利用刻蝕技術(shù)將圖形轉(zhuǎn)移到基片上。⑥在基片上覆蓋一層玻璃,放入高溫爐中燒結(jié),制得模型。通過測(cè)試,可視化模型耐壓超過8 MPa,彌補(bǔ)了微觀實(shí)驗(yàn)?zāi)P统叨刃?4.00 cm×4.00 cm)和實(shí)驗(yàn)壓力較低(小于0.2 MPa)的缺陷。
實(shí)驗(yàn)采用99.99%的高純氮?dú)庾鳛闅庠茨M儲(chǔ)層中的天然氣,實(shí)驗(yàn)用水是根據(jù)磨溪204井地層水分析資料配制的等礦化度標(biāo)準(zhǔn)鹽水,并用甲基藍(lán)染成藍(lán)色,這樣可以與圖像中的無色氮?dú)膺M(jìn)行區(qū)分。實(shí)驗(yàn)設(shè)備中微量泵的流量精度為0.001 mL/min,壓力精度為0.01 MPa,實(shí)驗(yàn)驅(qū)替壓差為0.2 MPa~0.5 MPa,實(shí)驗(yàn)溫度50 ℃,實(shí)驗(yàn)壓力8 MPa。
圖1 不同類型巖心CT掃描結(jié)果與可視化物理模型Fig. 1 The CT scan results and visualized physical models of different types cores
縫洞型碳酸鹽巖儲(chǔ)層氣水兩相微觀滲流機(jī)理可視化實(shí)驗(yàn)系統(tǒng)主要由可視化玻璃刻蝕模型,ISCO微量注入泵,回壓泵,中間容器,壓力傳感器,模型夾持器,氣體收集容器,光源,光學(xué)顯微鏡和數(shù)據(jù)、圖像收集裝置等組成(圖2)。實(shí)驗(yàn)步驟如下:
(1)將可視化模型安放在模型夾持器上,按實(shí)驗(yàn)流程連接好實(shí)驗(yàn)設(shè)備并檢查管線是否完好。
(2)對(duì)模型抽真空40 min后,調(diào)整回壓至6 MPa,圍壓至8 MPa,將模型加熱至50 ℃。
(3)將模型飽和地層水,然后用高純氮?dú)怛?qū)替至可視化模型不出水為止,期間用顯微照相和錄像設(shè)備記錄氣驅(qū)水動(dòng)態(tài),并計(jì)量出水量和出水速度。
(4)對(duì)模型再次抽真空,將殘留的液體抽干凈,用配置的地層水驅(qū)替至可視化模型不出氣為止,期間用顯微照相和錄像設(shè)備記錄水驅(qū)氣動(dòng)態(tài),并計(jì)量出氣量和出氣速度。
(5)實(shí)驗(yàn)完成后處理錄像和圖像,分析微觀驅(qū)替過程。
2.1.1 裂縫型模型水驅(qū)氣微觀滲流機(jī)理
在裂縫型模型中,裂縫是主要的滲流通道,與孔隙和溶洞相比具有更高的滲流能力。在水驅(qū)氣初期,氣水分布及流動(dòng)方式主要為“水包氣”。由于模型具有親水性,進(jìn)入模型中的水首先沿著裂縫壁面形成水膜,隨著水膜的逐漸增厚,裂縫中開始形成水流并以連續(xù)相的形式沿裂縫壁流動(dòng),而氣體以不連續(xù)的氣泡或氣柱在裂縫中間流動(dòng)。最終,水柱將裂縫中的大部分氣體驅(qū)出,同時(shí)在裂縫褶皺和縮頸部位發(fā)生卡斷現(xiàn)象,滯留了部分氣體(圖3)。在水驅(qū)氣中期,水竄現(xiàn)象已經(jīng)發(fā)生。隨著水的不斷注入,裂縫型模型中殘余氣柱或氣泡的能量不斷得到補(bǔ)充,當(dāng)受到的動(dòng)力大于阻力時(shí),氣柱或氣泡會(huì)聚能突破前方的水柱繼續(xù)流動(dòng)(圖4)。當(dāng)水驅(qū)氣進(jìn)入后期時(shí),殘余的少量氣泡會(huì)由于能量不足無法克服賈敏效應(yīng)而被束縛在微裂縫交叉部位,當(dāng)后續(xù)氣泡流經(jīng)該位置時(shí)會(huì)發(fā)生合并與能量傳遞,當(dāng)能量聚集到一定程度后,大氣泡的前端部分會(huì)分離出來形成小泡先通過,其余的部分只能繼續(xù)滯留,等待聚集更多的能力后才能通過(圖5)。
圖2 氣水兩相微觀滲流可視化實(shí)驗(yàn)流程Fig. 2 The visualization experiment process for gas-water two phase micro-seepage
圖3 裂縫型模型水驅(qū)氣初期“水包氣”形成過程Fig. 3 The formation process of “gas in water” in the early stage of water displacing gas
2.1.2 孔洞型模型水驅(qū)氣微觀滲流機(jī)理
對(duì)于孔洞型模型,既存在尺寸較小的孔隙,又存在尺寸較大的溶洞,其滲流能力主要受喉道的大小和分布控制。與裂縫型模型相比,毛管阻力相對(duì)較大,驅(qū)替壓差也相對(duì)較高。水驅(qū)氣過程中,由于模型親水,水首先在孔道壁上形成水膜,逐漸變厚形成水柱后與氣柱相互交替前進(jìn)(圖6(a))。同時(shí)由于部分孔隙內(nèi)前沿水膜在孔道交叉處匯聚形成水流,將中間段氣柱封鎖,形成封閉氣(圖6(b))。當(dāng)水進(jìn)入溶洞時(shí),首先將水膜延伸至溶洞壁面四圍,隨著水流的不斷注入,逐漸充填溶洞空間,并沿氣體滲流通道兩側(cè)壓縮氣體,最終將滲流通道壓死(圖6(c))。此時(shí)溶洞與孔道連接出口處會(huì)由于縮頸形成封閉氣。最終,溶洞大部分空間被水占據(jù),孔隙內(nèi)形成各種類型的封閉氣(圖6(d))。
2.1.3 縫洞型模型水驅(qū)氣微觀滲流機(jī)理
對(duì)于縫洞型微觀模型,在水驅(qū)氣初期,水先沿裂縫壁形成水膜(圖7(a)),隨后再聚集形成水柱,同時(shí)氣體在裂縫中央形成氣體心子或氣體段塞并隨水柱流動(dòng)(圖7(b))。由于裂縫的導(dǎo)流能力和通道尺寸遠(yuǎn)遠(yuǎn)大于孔洞模型中的孔道,所以水不僅以水膜的形式進(jìn)入溶洞壁面,還會(huì)直接以水柱形式進(jìn)入溶洞中部,并以近似活塞式的方式將溶洞內(nèi)的氣體快速驅(qū)出(圖7(c))。當(dāng)氣水界面到達(dá)出口處裂縫時(shí),會(huì)先沿主流線上的高滲裂縫突破,隨后再進(jìn)入與溶洞相連的微裂縫或低滲孔道(圖 7(d))。
圖4 裂縫型模型水驅(qū)氣中期氣柱聚能突破過程Fig. 4 The breakthrough process of gas column in the medium stage of water displacing gas
圖5 裂縫型模型水驅(qū)氣后期氣泡克服賈敏效應(yīng)突破過程Fig. 5 The breakthrough process of overcoming Jamin effect in the later stage of water displacing gas
圖6 孔洞型模型水驅(qū)氣過程Fig. 6 The water displacing gas process for the hole type model
圖7 縫洞型模型水驅(qū)氣過程Fig. 7 The water displacing gas process for the fracture-cave type model
2.2.1 盲端、角隅處形成的封閉氣
無論是孔隙、溶洞還是裂縫,其盲端或角隅處總會(huì)形成一定數(shù)量的封閉氣[17](表1)。盡管模型具有親水性,但是由于沒有形成有效的滲流通道,很難將其中的氣體驅(qū)出,特別是當(dāng)流動(dòng)通道上的壓力高于孔、洞、縫內(nèi)氣體壓力時(shí)。實(shí)驗(yàn)結(jié)果表明,只有降低驅(qū)替壓差才能將該封閉氣部分采出。這是因?yàn)楫?dāng)壓差降低時(shí),封閉氣發(fā)生膨脹后重新占據(jù)優(yōu)勢(shì)滲流通道,在后續(xù)水動(dòng)力的作用下被部分驅(qū)出,同時(shí)氣體的能量被逐漸消耗,壓力降低,當(dāng)與流動(dòng)通道的壓力達(dá)到平衡時(shí),氣體又一次被封閉。因此,在氣藏開發(fā)時(shí)只有降低氣藏壓力或在產(chǎn)能衰減時(shí)才能將這部分封閉氣采出,同時(shí)應(yīng)避免關(guān)井復(fù)壓操作,以免這部分封閉氣徹底被壓死。
2.2.2 繞流形成的封閉氣
繞流形成的封閉氣受毛管力和水動(dòng)力共同作用的影響,并與驅(qū)替壓差密切相關(guān)。3類模型中都存在該種形式的封閉氣,但是機(jī)理并不相同(表1)。對(duì)于孔隙型和孔洞型模型,當(dāng)驅(qū)替壓差較低時(shí),毛管力為氣水流動(dòng)的主要?jiǎng)恿?,此時(shí)封閉氣的形成主要來自繞流現(xiàn)象。當(dāng)水進(jìn)入多個(gè)孔徑不同的孔道后,在毛管力的作用下,水以較快的速度進(jìn)入較小的孔道[18],由于小孔道中水的滲流速度較快,氣體體積較小,水在模型出口處先于其他大孔道發(fā)生突破,隨后將大孔道中還未來得及驅(qū)出的氣體封閉起來形成封閉氣。當(dāng)驅(qū)替壓差較大時(shí),水動(dòng)力起主要作用時(shí),氣水兩相滲流機(jī)理與毛管力作主要?jiǎng)恿r(shí)正好相反。由于大孔道滲流阻力較小,水在水動(dòng)力的作用下優(yōu)先進(jìn)入大孔道,并且先于小孔道在模型出口處突破,隨后將小孔道中的氣體封閉起來,形成封閉氣。對(duì)于裂縫型和縫洞型模型,由于裂縫具有很高的滲流能力,毛管阻力很小,無論驅(qū)替壓差大小,注入水都會(huì)優(yōu)先進(jìn)入較大的裂縫并以較快的速度發(fā)生水竄,將部分孔隙和微裂縫中的氣體封閉起來,降低了主裂縫的補(bǔ)給能力和氣相滲透率,使實(shí)際氣藏的采氣速度和采出程度降低。
2.2.3 卡斷形成的封閉氣
實(shí)驗(yàn)結(jié)果表明,卡斷形成的封閉氣以不連續(xù)氣泡或氣柱的形式分布在孔道或裂縫中央(表1),造成該種形式封閉氣的原因是賈敏效應(yīng)。賈敏效應(yīng)又稱氣阻效應(yīng),由于地層孔隙結(jié)構(gòu)復(fù)雜,滲流通道尺寸存在較大差異,當(dāng)氣泡或者油滴通過細(xì)小孔隙喉道或裂縫褶皺變形部位時(shí),由于通道的前后半徑差使得氣泡或油滴兩端的弧面毛管力表現(xiàn)為阻力,若要通過半徑較小的通道必須拉長(zhǎng)并改變自身形狀,這種變形將消耗部分能量,從而減緩和限制氣泡或油滴的運(yùn)動(dòng),增加額外的阻力,這種阻力實(shí)質(zhì)是一種微毛管力效應(yīng)。對(duì)于孔隙型和孔洞型模型,當(dāng)氣水兩相流經(jīng)狹窄喉道時(shí),賈敏效應(yīng)的存在會(huì)產(chǎn)生附加阻力;同時(shí)由于模型親水,水相在喉道處使水膜增厚,易產(chǎn)生水鎖現(xiàn)象,使喉道直徑進(jìn)一步縮小,加劇了賈敏效應(yīng),增大了氣相流動(dòng)阻力,連續(xù)流動(dòng)的氣相必須收縮變形才能通過喉道。但是氣泡的收縮變形需要消耗自身能量,而原有能量只能使氣泡前端分離出來,形成小泡先行通過,其余部分只能滯留在喉道處,等待與后續(xù)氣泡碰撞聚能后,才能再次通過。而對(duì)于裂縫型和縫洞型模型,在比較粗糙的裂縫表面或者裂縫褶皺變形部位,同樣會(huì)因?yàn)橘Z敏效應(yīng)使連續(xù)流動(dòng)的氣體發(fā)生卡斷形成封閉氣柱或封閉氣泡。實(shí)驗(yàn)結(jié)果表明,可以通過提高驅(qū)替壓差來增大水動(dòng)力的方法采出卡斷封閉氣,也可以通過降低模型出口壓力的方法使卡斷封閉氣發(fā)生膨脹和聚集,利用自身的膨脹能力將其采出。
2.2.4 “H型”孔道形成的封閉氣
在孔隙型或孔洞型模型中還會(huì)形成“H型”孔道
封閉氣(表1)。該種類型封閉氣形成的主要機(jī)理有兩方面:一方面是由于水的毛管指進(jìn)優(yōu)先通過孔隙的兩條“邊路”向前突破,從而繞過了連接兩條“邊路”的“橋”;另一方面是由于當(dāng)水突破后,水會(huì)進(jìn)一步依靠模型的親水性進(jìn)入“H型”孔道的“橋”,壓縮“橋”上的氣體形成封閉氣。此類封閉氣也可以通過兩種方法將其采出:一是通過降低氣藏壓力使“橋”上的氣體膨脹到“邊路”上,在水動(dòng)力的作用下被帶出;二是通過增加驅(qū)替壓差來打破兩條“邊路”的壓力平衡,從而使“橋”上的氣體采出?!癏型”孔道形成的封閉氣是繞流的一種特殊形式。
表1 水驅(qū)氣過程封閉氣形成方式分類Table 1 Classi fi cation of con fi ned gas in the process of water displacing gas
2.2.5 “啞鈴型”通道形成的封閉氣
在縫洞型模型中會(huì)形成一種獨(dú)特的“啞鈴型”通道封閉氣(表1)。該種類型封閉氣的形成主要是由于注入水進(jìn)入溶洞時(shí),會(huì)優(yōu)先充填溶洞空間,同時(shí)將氣體沿出口端裂縫驅(qū)出,當(dāng)僅有一條裂縫同時(shí)作為2個(gè)溶洞的出口通道時(shí),兩部分氣體會(huì)同時(shí)向裂縫內(nèi)壓縮形成封閉氣。由于裂縫兩端溶洞內(nèi)的水體能量較大,很難通過改變驅(qū)替壓差的途徑將這部分封閉氣采出。
2.2.6 微裂縫縫網(wǎng)形成的封閉氣
對(duì)于裂縫型和縫洞型模型,會(huì)在微裂縫縫網(wǎng)中封閉一部分氣柱(表1)。這是由于大裂縫具有更高的導(dǎo)流能力和極低的毛管阻力,注入水會(huì)優(yōu)先進(jìn)入大裂縫并以較快的速度突破,這使得被大裂縫切割的微裂縫被水包圍,堵塞了氣流通道,使原本統(tǒng)一的壓力系統(tǒng)被分割成多個(gè),在微裂縫處形成“氣死區(qū)”[19]。該部分封閉氣只能通過降低模型出口壓力,利用自身的膨脹能力進(jìn)入主裂縫后被采出。6種類型封閉氣的形成機(jī)理與主要影響因素如表2所示。
水驅(qū)氣實(shí)驗(yàn)?zāi)M的是有水氣藏在開發(fā)過程中水侵對(duì)氣井生產(chǎn)規(guī)律的影響,對(duì)實(shí)驗(yàn)過程中不同模型的氣水分布進(jìn)行定量表征對(duì)該類氣藏的有效開發(fā)具有指導(dǎo)意義。由于氣水兩相微觀滲流可視化實(shí)驗(yàn)持續(xù)時(shí)間很短,期間很難對(duì)不同時(shí)刻的出氣量和出水量進(jìn)行計(jì)量,也很難利用常規(guī)方法對(duì)不同時(shí)刻下模型內(nèi)的氣水分布進(jìn)行定量表征。為了解決這一問題,決定采用ImageJ灰度分析法,首先利用視頻處理軟件截取水驅(qū)氣視頻中不同時(shí)刻下的可視化圖像,并利用Photoshop軟件進(jìn)行預(yù)處理,主要是對(duì)圖像的亮度進(jìn)行均勻調(diào)整。在此基礎(chǔ)上利用ImageJ軟件的圖像識(shí)別功能對(duì)圖像的灰度值進(jìn)行區(qū)分,通過調(diào)整閾值先識(shí)別出模型初始條件下的儲(chǔ)集空間和玻璃顆粒,再統(tǒng)計(jì)像素點(diǎn)求出儲(chǔ)集空間所占的總面積,進(jìn)而計(jì)算出孔隙度。之后利用Photoshop軟件區(qū)分不同時(shí)刻下的氣水分布,分離出注入水,最后再利用ImageJ軟件計(jì)算含水飽和度和含氣飽和度。
通過ImageJ灰度分析法得到的裂縫型、孔洞型和縫洞型模型在不同時(shí)刻下的含氣飽和度如圖8所示。裂縫型模型由于存在裂縫這一高滲通道,水驅(qū)氣初期產(chǎn)氣速度很快,含氣飽和度快速下降,此時(shí)驅(qū)替的都是裂縫和大孔道中的氣體,同時(shí)水竄現(xiàn)象嚴(yán)重,無水采收期較短,無水采出程度較低,只有51.3%。模型見水后被裂縫和大孔道封鎖的低滲區(qū)域內(nèi)氣體很難被采出,封閉氣大部分分布在盲端、角隅處以及微裂縫縫網(wǎng)內(nèi),最終采出程度只有63.5%。為了提高裂縫型儲(chǔ)層的采出程度,在進(jìn)行氣井配產(chǎn)時(shí)需要嚴(yán)格控制生產(chǎn)壓差,不僅能減緩水侵速度,延長(zhǎng)無水采收期,還可以提高波及效率。而在氣井見水后,可以采取逐級(jí)降低生產(chǎn)壓差的方式,使盲端、角隅和微裂縫縫網(wǎng)內(nèi)的封閉氣發(fā)生膨脹后重新占據(jù)優(yōu)勢(shì)滲流通道,在后續(xù)水動(dòng)力的作用下被驅(qū)出;同時(shí)也要避免關(guān)井復(fù)壓操作,以免這部分封閉氣徹底被壓死。
表2 封閉氣主要影響因素及形成機(jī)理Table 2 In fl uencing factors and forming mechanisms of trapped gas
圖8 水驅(qū)氣實(shí)驗(yàn)氣水飽和度分布曲線Fig. 8 The saturation distribution curves of water and gas in water displacing gas experiment
孔洞型模型由于滲流通道較小,滲流阻力較大,孤立的溶洞主要在局部影響水侵的流動(dòng)規(guī)律,而整體水侵前緣推進(jìn)和孔隙型儲(chǔ)層類似,近似均勻推進(jìn),因此溶洞的存在主要為氣藏提供了儲(chǔ)集空間,對(duì)整體水侵前緣的推進(jìn)影響不大,這使得初期產(chǎn)氣速度相對(duì)較慢,無水采收期較長(zhǎng),無水采收程度較高,達(dá)到65%,模型見水后主要形成繞流、卡段形式的封閉氣,采出程度最終為68.4%。為了提高孔洞型儲(chǔ)層的采出程度,建議采用逐級(jí)加壓的方式進(jìn)行開發(fā),既可以在低壓差時(shí)利用毛管力采出小孔道中的氣體,又可以在高壓差時(shí)利用水動(dòng)力采出因繞流封閉在大孔道中的氣體,還可以使卡斷形成的封閉氣泡聚能克服賈敏效應(yīng),使得采出程度進(jìn)一步提高。
對(duì)于縫洞型模型,氣水飽和度變化規(guī)律與裂縫型模型相似,因此封閉氣類型與開發(fā)方式也與裂縫型模型相同。水體會(huì)沿著裂縫快速向前推進(jìn),通過溶洞時(shí)由于其較大的儲(chǔ)氣空間減緩了水侵的進(jìn)度,這使得無水采出程度和最終采出程度較裂縫型高一些,無水采收期也稍長(zhǎng)一些。
氣驅(qū)水實(shí)驗(yàn)主要用來模擬氣藏的形成過程和地層水的流動(dòng)過程,對(duì)于本文研究的3類微觀模型,由于表面張力和毛細(xì)管力的共同作用,在飽和水過程中水會(huì)優(yōu)先占據(jù)大通道壁面和細(xì)小孔喉,最終完全充滿裂縫、溶洞、孔隙與喉道,但在盲端和角隅處只能部分飽和(圖9(a)、圖10(a)與圖11(a))。
3.1.1 裂縫型模型氣驅(qū)水微觀滲流機(jī)理
在氣驅(qū)水實(shí)驗(yàn)中,氣體的性質(zhì)決定了其能夠進(jìn)入極小空間進(jìn)行驅(qū)替。對(duì)于裂縫型模型,氣體很快便將裂縫中間部位的水驅(qū)出(圖9(b)),隨后將裂縫褶皺、角隅、盲道等位置的水驅(qū)出(圖9(c))。在氣驅(qū)水的后期,隨著氣體不斷注入,帶走裂縫壁面上的殘余水,使水膜逐漸變薄,并在裂縫交叉處匯聚成水珠被驅(qū)出,最后僅在裂縫壁上留下一層薄薄的水膜(圖9(d)),可以通過增大驅(qū)替壓差的方式使該部分殘余水膜進(jìn)一步變薄。
3.1.2 孔洞型模型氣驅(qū)水微觀滲流機(jī)理
對(duì)于孔洞型模型,氣體進(jìn)入后首先沿優(yōu)勢(shì)通道快速推進(jìn),隨著注氣量的增多,氣體逐漸占據(jù)大部分孔道空間,并將其中大部分水驅(qū)出(圖10(b))。當(dāng)氣體流經(jīng)溶洞時(shí),會(huì)先在溶洞水體外圍形成一圈細(xì)小的氣流通道,并隨著氣體的進(jìn)入逐漸擴(kuò)大通道,將邊部水驅(qū)出,最終在溶洞中央形成一個(gè)相對(duì)圓潤(rùn)的水體團(tuán),當(dāng)水體團(tuán)被驅(qū)替到一定程度后,形狀不再發(fā)生變化,此時(shí)氣體無法再將水體團(tuán)驅(qū)出(圖10(c))。氣驅(qū)水后期,注入的氣體將孔隙壁上的水膜聚集,當(dāng)達(dá)到一定厚度時(shí)被驅(qū)替而出,最終在細(xì)長(zhǎng)孔道處、狹窄喉道處和溶洞中部形成殘余水(圖10(d))。
3.1.3 縫洞型模型氣驅(qū)水微觀滲流機(jī)理
在縫洞型模型氣驅(qū)水開始后,氣體迅速將裂縫中部的水驅(qū)出,形成氣流通道(圖11(b))。當(dāng)氣體流經(jīng)溶洞時(shí),會(huì)優(yōu)先驅(qū)替溶洞外圍的水,形成氣體的滲流通道,并把中部的水體包圍起來(圖11(c))。隨著氣體的不斷注入,不斷剝蝕中部水體,并在出口裂縫處匯聚成水珠后被驅(qū)出,使溶洞內(nèi)水體逐漸縮小至一定程度后不再發(fā)生變化,形成殘余水團(tuán),也有一部分溶洞由于被多條裂縫溝通,水被全部驅(qū)出,僅在孔洞壁上形成束縛水膜(圖11(d))。
3.2.1 孔道、裂縫壁面上的殘余水膜
對(duì)于裂縫型和縫洞型模型,裂縫是主要的滲流通道。氣驅(qū)水過程中,裂縫中的水首先被驅(qū)出,由于毛管阻力很小,殘余水主要以“薄水膜”形式賦存在連通性較差的裂縫壁面(圖12(a)),對(duì)氣相滲流影響較小。而對(duì)于孔洞型模型,孔道中的殘余水在高毛管力的作用下,以“厚水膜”形式賦存,殘余水飽和度較高,這會(huì)使得氣相滲流通道變窄并增大流動(dòng)阻力,要將此種類型殘余水驅(qū)替出需要較大的驅(qū)替壓差和較長(zhǎng)的驅(qū)替時(shí)間。
3.2.2 溶洞中部的殘余水團(tuán)
溶洞中部的殘余水是孔洞型和縫洞型模型中殘余水的主要形式之一。這是因?yàn)楫?dāng)氣體流經(jīng)溶洞時(shí),會(huì)優(yōu)先沿溶洞水體外圍突破,形成并逐漸拓寬滲流通道,剝蝕中部水體并在出口通道處匯聚成水珠后被驅(qū)出。當(dāng)水體團(tuán)被驅(qū)替到一定程度后,優(yōu)勢(shì)滲流通道已經(jīng)完全形成,此時(shí)形狀不再發(fā)生變化,滯留在溶洞中部形成殘余水(圖12(b)),該部分殘余水可以通過增大驅(qū)替壓差的方式將其驅(qū)出。
3.2.3 狹窄喉道處、狹長(zhǎng)孔道處的殘余水柱
由于殘余水在喉道處受到的毛管阻力非常大,在氣驅(qū)水過程中容易發(fā)生卡斷,以“水柱”的形式滯留在整個(gè)喉道處(圖12(c)),堵塞滲流通道,這也是儲(chǔ)層發(fā)生水侵后氣相滲透率大幅下降的主要原因。在多條孔道同時(shí)存在的情況下,狹長(zhǎng)孔道的毛管阻力比較寬、較短的孔道要大得多,需要耗費(fèi)較多的能量才能將水驅(qū)出,所以在狹長(zhǎng)孔道處也容易形成殘余水。3種類型殘余水的形成機(jī)理與影響因素如表3所示。
圖9 裂縫型模型氣驅(qū)水過程Fig. 9 The gas displacing water process for the fracture type model
圖10 孔洞型模型氣驅(qū)水過程Fig. 10 The gas displacing water process for the hole type model
圖11 縫洞型模型氣驅(qū)水過程Fig. 11 The gas displacing water process for the fracture-cave type model
圖12 氣驅(qū)水實(shí)驗(yàn)殘余水形成模式Fig. 12 The formation mode of residual water in gas displacing water experiment
表3 殘余水主要影響因素及形成機(jī)理Table 3 In fl uencing factors and forming mechanisms of irreducible water
(1)通過將巖心CT掃描技術(shù)與激光刻蝕技術(shù)相結(jié)合,以四川盆地震旦系儲(chǔ)層中的裂縫型、孔洞型和縫洞型3類碳酸鹽巖巖心CT掃描圖像為模板,設(shè)計(jì)并研制了3類巖心的激光刻蝕透明仿真模型,實(shí)現(xiàn)了對(duì)實(shí)際儲(chǔ)層中孔、縫、洞分布的精準(zhǔn)刻畫。通過開展氣水兩相微觀滲流實(shí)驗(yàn),得到了3類模型水驅(qū)氣、氣驅(qū)水過程中的微觀可視化圖像,并通過進(jìn)一步分析得到了縫洞型碳酸鹽巖儲(chǔ)層氣水兩相滲流過程中的機(jī)理與特征規(guī)律。
(2)通過水驅(qū)氣可視化實(shí)驗(yàn)發(fā)現(xiàn),繞流、卡斷、盲端和角隅處形成的封閉氣普遍存在于裂縫型、孔洞型和縫洞型3類模型中?!癏型”孔道形成的封閉氣只存在于孔洞型模型中,“啞鈴型”通道形成的封閉氣只存在于縫洞型模型中,而微裂縫縫網(wǎng)形成的封閉氣存在于裂縫型和縫洞型2類模型中。上述形式的封閉氣通過改變驅(qū)替壓差或降低出口壓力的方法只能部分采出。
(3)通過氣驅(qū)水可視化實(shí)驗(yàn)發(fā)現(xiàn),孔道、裂縫壁面上的束縛水膜,溶洞中部的圓潤(rùn)水團(tuán),狹窄喉道處的卡斷水柱和狹長(zhǎng)孔道處的滯留水柱是3類模型殘余水的主要形式,這些形式的殘余水均可通過增大驅(qū)替壓差的方法部分采出。
(4)通過ImageJ灰度分析法對(duì)氣驅(qū)水實(shí)驗(yàn)不同時(shí)刻下的氣水分布進(jìn)行了定量表征。研究結(jié)果表明,裂縫型模型水竄最嚴(yán)重,無水采收期最短,無水采出程度和最終采出程度最低。縫洞型模型由于存在溶洞減緩了水侵的進(jìn)度,使得采出程度稍高于裂縫型模型??锥葱湍P陀捎跐B流阻力最大,水驅(qū)前緣推進(jìn)均勻,無水采收期最長(zhǎng),無水采出程度和最終采出程度也最高。
[1] 彭松,郭平.縫洞型碳酸鹽巖凝析氣藏注水開發(fā)物理模擬研究[J].石油實(shí)驗(yàn)地質(zhì),2014,36(5):645-649.[PENG S, GUO P. Physical simulation of exploiting fractured-vuggy carbonate gas condensate reservoirs by water injection[J]. Petroleum Geology & Experiment,2014, 36(5): 645-649.]
[2] 李陽(yáng).塔河油田碳酸鹽巖縫洞型油藏開發(fā)理論及方法[J].石油學(xué)報(bào),2013, 34(1):115-121.[LI Y. The theory and method for development of carbonate fractured-cavity reservoirs in Tahe oil fi eld[J]. Acta Petrolei Sinica, 2013, 34(1): 115-121.]
[3] 高樹生,劉華勛,任東,等.縫洞型碳酸鹽巖儲(chǔ)層產(chǎn)能方程及其影響因素分析[J].天然氣工業(yè),2015,35(9):48-54.[GAO S S, LIU H X, REN D, et al. Deliverability equation of fracture-cave carbonate reservoirs and its in fl uential factors[J]. Natural Gas Industry, 2015,35(09): 48-54.]
[4] 賈愛林,閆海軍.不同類型典型碳酸鹽巖氣藏開發(fā)面臨問題與對(duì)策[J].石油學(xué)報(bào),2014,35(3):519-527.[JIA A L, YAN H J. Problems and countermeasures for various types of typical carbonate gas reservoirs development[J]. Acta Petrolei Sinica, 2014, 35(3): 519-527.]
[5] 鄒才能,杜金虎,徐春春,等.四川盆地震旦系—寒武系特大型氣田形成分布、資源潛力及勘探發(fā)現(xiàn)[J].石油勘探與開發(fā),2014,41(3):278-293.[ZHOU C N, DU J H, XU C C, et al. Formation, distribution, resource potential and discovery of the Sinian-Cambrian giant gas fi eld, Sichuan Basin, SW China[J]. Petroleum Exploration & Development, 2014, 41(3): 278-293.]
[6] 周克明,李寧,張清秀,等.氣水兩相滲流及封閉氣的形成機(jī)理實(shí)驗(yàn)研究[J].天然氣工業(yè),2002,22(增刊1):122-125.[ZHOU K M,LI N, ZHANG Q X, et al. Experimental research on gas-water two phase fl ow and con fi ned gas formation mechanism[J]. Natural Gas Industry, 2002, 22(Supplement 1): 122-125.]
[7] 吳建發(fā),郭建春,趙金洲.裂縫性地層氣水兩相滲流機(jī)理研究[J].天然氣工業(yè),2004,24(11):85-87.[WU J F, GUO J C, ZHAO J Z.Study on gas-water two phase percolation mechanism for fractured formation[J]. Natural Gas Industry, 2004, 24(11): 85-87.]
[8] 朱華銀,周娟,萬玉金,等.多孔介質(zhì)中氣水滲流的微觀機(jī)理研究[J].石油實(shí)驗(yàn)地質(zhì),2004,26(6):571-573.[ZHU H Y, ZHOU J, WAN Y J, et al. Microscopic mechanism study of gas-water fl ow in porous media[J]. Petroleum Geology & Experiment, 2004, 26(6): 571-573.]
[9] 李登偉,張烈輝,周克明,等.可視化微觀孔隙模型中氣水兩相滲流機(jī)理[J].中國(guó)石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2008,32(3):80-83.[LI D W, ZHANG L H, ZHOU K M, et al. Gas-water two-phase fl ow mechanism in visual microscopic poremode[J]. Journal of China University of Petroleum, 2008, 32(3): 80-83.]
[10] 樊懷才,鐘兵,李曉平,等.裂縫型產(chǎn)水氣藏水侵機(jī)理研究[J].天然氣地球科學(xué),2012,23(6):1 179-1 184.[FAN H C, ZHONG B, LI X P, et al. Studies on water invasion mechanism of fractured-watered gas reservoir[J]. Natural Gas Geoscience, 2012, 23(6): 1 179-1 184.]
[11] 胡勇,邵陽(yáng),陸永亮,等.低滲氣藏儲(chǔ)層孔隙中水的賦存模式及對(duì)氣藏開發(fā)的影響[J].天然氣地球科學(xué),2011,22(1):176-181.[HU Y,SHAO Y, LU Y L, et al. Experimental study on occurrence models of water in pores and the in fl uencing to the development of tight gas reservoir[J]. Natural Gas Geoscience, 2011, 22(1): 176-181.]
[12] 鄢友軍,陳俊宇,郭靜姝,等.龍崗地區(qū)儲(chǔ)層微觀鮞粒模型氣水兩相滲流可視化實(shí)驗(yàn)及分析[J].天然氣工業(yè),2012,32(1):64-66.[Yan Y J, CHEN J Y, GUO J S, et al. A visualized experiment on gas-water two-phase seepage through oolitic reservoirs in the Longgang Gas Field, Sichuan Basin[J]. Natural Gas Industry, 2012, 32(1): 64-66.]
[13] 陳朝暉,謝一婷,鄧勇.澀北氣田疏松砂巖氣藏微觀氣水驅(qū)替實(shí)驗(yàn)[J].西南石油大學(xué)學(xué)報(bào)(自然科學(xué)版),2013,35(4):139-144.[CHEN Z H, XIE Y T, DENG Y. Experiment of microscopic displacement of gas and water in loose sandstone gas reservoir of Sebei[J]. Journal of Southwest Petroleum University (Science & Technology Edition), 2013, 35(4): 139-144.]
[14] 方飛飛,李熙喆,高樹生,等.邊、底水氣藏水侵規(guī)律可視化實(shí)驗(yàn)研究[J].天然氣地球科學(xué),2016,27(12):1-7.[Fang F F, LI X Z, GAO S S, et al. Visual simulation experimental study on water invasion rules of gas reservoir with edge and bottom water[J]. Natural Gas Geoscience, 2016, 27(12): 2 246-2 252.]
[15] 高樹生,胡志明,劉華勛,等.不同巖性儲(chǔ)層的微觀孔隙特征[J].石油學(xué)報(bào),2016,37(2):248-256.[GAO S S, HU Z M, LIU H X, et al.Microscopic pore characteristics of different lithological reservoirs[J]. Acta Petrolei Sinica, 2016, 37(2): 248-256.]
[16] 王家祿,高建,劉莉.應(yīng)用CT技術(shù)研究巖石孔隙變化特征[J].石油學(xué)報(bào),2009,30(6):887-893.[WANG J L, GAO J, LIU L. Porosity characteristics of sandstone by X-ray CT scanning system[J]. Acta Petrolei Sinica, 2009, 30(6): 887-893.]
[17] RENG D M, ZHANG L H, ZU S Q. Study on transportation numerical simulation of coalbed methane reservoir[J]. Journal of Hydro-dunamics (Series B), 2003, 15(6): 63-67.
[18] ZHANG L H, FENG G Q, LI X P, et al. Water breakthrough simulation in naturally fractured gas reservoirs with water drive[J]. Journal of Hydrodunamics (Series B), 2005, 17(4): 466-472.
[19] CIESLINSKI J T, MOSDORF R. Gas bubble dynamics-experiment and fractal analysis [J]. International Journal of Heat & Mass Transfer, 2005, 48(9): 1 808-1 818.
AbstractMost fracture-cavity carbonate reservoirs in the Sichuan Basin are water-bearing gas reservoirs and the gas-water two phase seepage mechanism is complex. However, the existing visual micromodels of gas-water two phase micro seepage are based on an ideal porous structure or casting sheet image, which cannot describe the distribution of pores, fractures and cavities in actual reservoirs. In addition, the research into mechanism is focused on the pore type and fracture type reservoirs, and there is a lack of the understanding of the cavity type and fracture-cavity type reservoirs. By combining the CT scanning technology and laser etching techniques, the CT scan images of fracture type, cavity type and fracture-cavity type carbonate cores in the Sinian system of the Sichuan Basin are used as the template to design and develop the three types of visual micromodels. Then the mechanism of gas-water two phase micro seepage and the formation mechanisms of trapped gas and irreducible water were studied. The quantitative characterization of gas-water distribution was also obtained by the gray scale analysis. The results show that the seepage mechanisms of three models are different in process of gas fl ooding and water fl ooding. The trapped gas caused by circum fl uence, cut off phenomenon and blind corners is prevalent in three types of models. In addition, the special trapped gas is also remains in “H” channels, “dumbbell-shaped” channels and micro-fracture channels. The irreducible water fi lms on the channel walls, the round water group in the middle of the cavities, the water columns into the narrow throats are the main forms of residual water. The fracture type model has the most serious water channeling phenomenon, the shortest water-free production period and the lowest recovery ef fi ciency, while the cavity type model has a uniform water drive front, the longest water-free production period and the highest recovery ef fi ciency. This study provides an accurate description of gas-water two phase micro-seepage mechanisms of different formations in fracture-cavity carbonate reservoirs. This provides a theoretical basis for the ef fi cient development of similar gas reservoirs.
Keywordsfracture-cavity carbonate reservoirs; gas-water two-phase fl ow; visual micromodels; trapped gas; irreducible water;gray intensity analysis
(編輯 馬桂霞)
Visual experimental investigation of gas-water two phase micro seepage mechanisms in fracture-cavity carbonate reservoirs
WANG Lu1, YANG Shenglai1, LIU Yicheng2, WANG Yunpeng1, MENG Zhan1, HAN Wei1, QIAN Kun1
1 State Key Laboratory of Petroleum Resource and Engineering, China University of Petroleum-Beijing, Beijing 102249, China
2 Exploration and Development Research Institute of Southwest Oil & Gas Field Company, PetroChina, Chengdu 610041, China
2016-12-22
國(guó)家科技重大專項(xiàng)“深層碳酸鹽巖氣藏高效開發(fā)技術(shù)” (2016ZX05015-003)和國(guó)家重點(diǎn)基礎(chǔ)研究發(fā)展計(jì)劃(“973”計(jì)劃)(2015CB250900)聯(lián)合資助
王璐, 楊勝來, 劉義成, 王云鵬, 孟展, 韓偉, 錢坤.縫洞型碳酸鹽巖儲(chǔ)層氣水兩相微觀滲流機(jī)理可視化實(shí)驗(yàn)研究. 石油科學(xué)通報(bào),2017, 03: 364-376
WANG Lu, YANG Shenglai, LIU Yicheng, WANG Yunpeng, MENG Zhan, HAN Wei, QIAN Kun. Visual experimental investigation of gas-water two phase micro seepage mechanisms in fracture-cavity carbonate reservoirs. Petroleum Science Bulletin, 2017, 03: 364-376.
10.3969/j.issn.2096-1693.2017.03.034
doi:10.3969/j.issn.2096-1693.2017.03.034