縫洞型油藏溶洞內(nèi)水驅(qū)流態(tài)研究

劉承婷 李江 宋洋

摘要：為研究縫洞型油藏溶洞中注水速度和黏度對(duì)水驅(qū)波及面積和油水界面的影響，利用Fluent數(shù)值模擬軟件從流體力學(xué)角度對(duì)單縫洞油藏溶洞內(nèi)流體流動(dòng)進(jìn)行了數(shù)值模擬，采用VOF模型對(duì)油水兩項(xiàng)界面進(jìn)行跟蹤。研究發(fā)現(xiàn)，入口水速一定時(shí)，波及面積隨著黏度的增大逐漸增大，油水界面呈對(duì)稱均勻推進(jìn)，黏度較小的原油水相會(huì)趨向中心線向洞口流動(dòng)并產(chǎn)生漩渦，隨著黏度逐漸增大，漩渦逐漸消失，驅(qū)替方式逐漸轉(zhuǎn)變成活塞式驅(qū)替；原油黏度一定時(shí)，波及面積隨著入口速度的增加而減小，油水界面呈對(duì)稱趨勢(shì)均勻推進(jìn)，水速較大時(shí)水相會(huì)趨向中心線向洞口流動(dòng)并有漩渦產(chǎn)生，隨著入口水速增大，趨勢(shì)逐漸明顯。研究結(jié)果可為進(jìn)一步完善縫洞型油藏溶洞內(nèi)水驅(qū)過(guò)程中油水兩相的流態(tài)分析提供參考。

關(guān)鍵詞：多相流體力學(xué)；單縫洞油藏；流體力學(xué)；波及面積；油水界面；Fluent數(shù)值模擬

中圖分類(lèi)號(hào)：TE341文獻(xiàn)標(biāo)志碼：Adoi： 10.7535/hbgykj.2018yx03004

縫洞型碳酸鹽巖油藏主要由細(xì)小的裂縫與溶洞連接而成，縫洞作為主要儲(chǔ)集體，裂縫為流動(dòng)通道的強(qiáng)非均質(zhì)油藏[1-3]。裂縫低孔高滲的特性和溶洞高孔高滲的特性符合單孔介質(zhì)的特征，使縫洞型油藏不同于常規(guī)的雙孔介質(zhì)油藏，開(kāi)采機(jī)理與普通雙孔介質(zhì)油藏有所不同[4-6]。邢翠巧等[7]針對(duì)縫洞型碳酸鹽巖油藏經(jīng)過(guò)構(gòu)造及巖溶作用之后的情況，建立了溶洞及裂縫發(fā)育的雙孔模型，研究了縫洞型碳酸鹽巖油藏實(shí)際儲(chǔ)層特征和流體的流動(dòng)機(jī)理。KANG等[8]和WU等[9]應(yīng)用連續(xù)介質(zhì)滲流理論建立了縫洞型油藏多相流問(wèn)題的數(shù)值模擬方法。鄭小敏等[10]利用全直徑縫洞型巖芯模型，通過(guò)物理實(shí)驗(yàn)?zāi)M水驅(qū)油實(shí)驗(yàn)研究了縫洞型碳酸鹽巖油藏的滲流特點(diǎn)。康志宏[11]針對(duì)塔河油田4區(qū)奧陶系縫洞型碳酸鹽巖油藏開(kāi)展?jié)B流機(jī)理模擬試驗(yàn)研究，得到碳酸鹽巖油藏中裂縫（或大的溶洞）是主要滲流通道的結(jié)論。

第3期劉承婷，等：縫洞型油藏溶洞內(nèi)水驅(qū)流態(tài)研究河北工業(yè)科技第35卷 縫洞型油藏由于儲(chǔ)集體類(lèi)型多樣，微孔隙、微裂縫、裂縫和中型溶洞以及大型溶洞、大型裂縫的分布，在尺度和空間結(jié)構(gòu)上都具有很大的差異，流體流動(dòng)類(lèi)型復(fù)雜，滲流、裂縫流和洞穴流都有存在[12-15]。趙歡等[16]利用油藏?cái)?shù)值模擬方法，研究了不同聚合物參數(shù)對(duì)“二三結(jié)合”開(kāi)發(fā)模式原射孔層和補(bǔ)孔層的壓力、含水飽和度等滲流特征的影響。劉金玉[17]利用COMSOL Multiphysics多物理場(chǎng)軟件對(duì)單縫概念模型進(jìn)行了水驅(qū)油數(shù)值模擬，分析了油水流態(tài)、波及面積、水驅(qū)前緣及剩余油分布的基本特征。目前對(duì)縫洞型碳酸鹽巖油藏水驅(qū)的研究多是從滲流角度出發(fā)，但是由于縫洞型油藏內(nèi)部結(jié)構(gòu)復(fù)雜，裂縫和溶洞的存在情況無(wú)法準(zhǔn)確掌握，特別是針對(duì)一些大裂縫和大溶洞，內(nèi)部屬于未填充狀態(tài)，多孔介質(zhì)沒(méi)有在內(nèi)部分布，所以內(nèi)部流體流動(dòng)不屬于滲流范疇，它的流動(dòng)規(guī)律符合Navier-Stokes方程，因此需采用流體力學(xué)對(duì)其進(jìn)行研究。

本文從流體力學(xué)角度出發(fā)，針對(duì)未填充縫洞型碳酸鹽巖油藏進(jìn)行了油水兩相流動(dòng)模擬，研究單裂縫-溶洞中流體流動(dòng)時(shí)油藏內(nèi)原油黏度、注水速度對(duì)油水界面、波及面積的影響。

1數(shù)學(xué)模型

1.1連續(xù)方程

連續(xù)性方程是流體運(yùn)動(dòng)學(xué)的基本方程，是質(zhì)量守恒原理的流體力學(xué)表達(dá)式。

連續(xù)方程微分形式見(jiàn)式（1）。

ρt+Δ·（ρu）=0，（1）

式中：ρ表示流體密度，kg/s；t表示時(shí)間，s；u表示流速，m/s。

模型設(shè)定的油水兩相流動(dòng)為穩(wěn)態(tài)流動(dòng)，即：

Δ·（ρu）=0 。 （2）

油水兩相都是不可壓縮流體，即：

ρt=0 。 （3）

將以上條件代入連續(xù)方程，得到油水兩相流的連續(xù)方程表達(dá)式：

Δ·u=0 。 （4）

式（4）即用來(lái)描述縫洞型介質(zhì)兩相流動(dòng)機(jī)理的連續(xù)方程，可以通過(guò)連續(xù)方程得到流體的流速，代入納維-斯托克斯方程。

1.2納維-斯托克斯方程

納維-斯托克斯方程簡(jiǎn)稱N-S方程，用來(lái)描述黏性不可壓縮流體動(dòng)量守恒的運(yùn)動(dòng)方程[18]。

溶洞型油藏中油水兩相流動(dòng)屬于黏性流體的不可壓縮流動(dòng)，由于需要研究油水界面在不同黏度原油和不同入口速度下的推進(jìn)情況，所以采用張量形式的不可壓縮黏性N-S方程，即：

ρduidt=ρFipijxj+ρν2uix2j，（5）

式中：ρduidt表示單位體積上的慣性力；ρFipijxj表示單位體積上的質(zhì)量力；ρν2uix2j表示單位體積上的黏滯力。

由于油水兩項(xiàng)在大溶洞和大裂縫中的運(yùn)動(dòng)屬于自由流動(dòng)，遵循N-S定律，同時(shí)為了結(jié)合VOF方法追蹤溶洞中油水界面問(wèn)題，所以采用張量形式的不可壓縮黏性N-S方程，用于研究未填充溶洞中的油水兩相流動(dòng)問(wèn)題。

1.3VOF方法及原理

VOF方法基本原理是通過(guò)對(duì)網(wǎng)格單元中流體之間的體積比函數(shù)F來(lái)確定自有界面，追蹤界面的推進(jìn)情況。由Hirt和Nichols首次提出后經(jīng)NOH等改進(jìn)，在保證質(zhì)量守恒的前提下，PUCKETT等采用方向分裂算法對(duì)輸送方程進(jìn)行優(yōu)化，利用流體體積函數(shù)方程張量表示方式如下：

Ft+ujFxj=0。（6）

由于油的密度小于水的密度，兩相之間存在明顯界面，所以采用VOF模型對(duì)油水界面進(jìn)行跟蹤。根據(jù)單位網(wǎng)格內(nèi)油水體積比F得到油水界面形狀和位置。當(dāng)F=0時(shí)，說(shuō)明單元內(nèi)全部是油相；當(dāng)F=1時(shí)，說(shuō)明單元內(nèi)全部是水相，已經(jīng)將油驅(qū)出；當(dāng)0根據(jù)流體體積函數(shù)方程張量表示方式，通過(guò)流體流動(dòng)時(shí)間、速度以及流動(dòng)方向上的位置可以計(jì)算出單位網(wǎng)格內(nèi)油水體積比，利用瞬態(tài)計(jì)算方法，得到不同時(shí)刻的油水界面分布情況。

2數(shù)值計(jì)算

2.1初始條件

單裂縫是縫洞型油藏的基本流動(dòng)單元，對(duì)流體在油藏中的流動(dòng)起決定性作用。假設(shè)油和水都是不可壓縮流體，密度為常數(shù)。由于溶洞三度空間長(zhǎng)度相近，同時(shí)為了減小計(jì)算量，將溶洞簡(jiǎn)化為正方形二維模型，出入口位于左右兩側(cè)中心，左端為出口，右端為入口。入口定義為速度入口，出口定義為壓力出口。

2.2模型參數(shù)及網(wǎng)格劃分

參數(shù)確定：直徑大于100 mm屬于巨洞；直徑在10～100 mm屬于大洞；直徑在5～10 mm之間屬于中洞；直徑在2～5 mm之間屬于小洞。由于溶洞空隙三度空間相近，且都大于2×103 μm，所以選用三度尺寸均為50 mm的溶洞，縫寬為0.6 mm。實(shí)際縫洞型碳酸鹽巖油藏地層原油黏度一般為2～30 mPa·s，所以選用原油黏度為5，15和25 mPa·s。原油密度選取860 kg/m3，水的密度為996 kg/m3。

溶洞模型為邊長(zhǎng)50 mm的正方形二維模型，出入口縫寬為0.6 mm，利用gambit對(duì)模型進(jìn)行六面體網(wǎng)格劃分，見(jiàn)圖1。根據(jù)雷諾數(shù)計(jì)算，采用Fluent中的層流模型，且為了方便追蹤油水界面，采用VOF方法進(jìn)行瞬態(tài)計(jì)算，收斂精度為0.000 01。

首先，計(jì)算入口速度為0.08 m/s、黏度為5 mPa·s的溶洞內(nèi)油藏的流態(tài)，得到水驅(qū)油水分界面分布，波及面積及剩余油分布，依次改變黏度為15和25 mPa·s，得到不同黏度下的油水界面發(fā)展情況、波及面積及剩余油分布情況。改變?cè)宛ざ龋玫较嗤肟谒俣认虏煌ざ鹊挠退植记闆r，分析原油黏度和流態(tài)之間關(guān)系；再改變?nèi)肟谒俣?，依次設(shè)定入口水速為0.08，0.16，0.24 m/s，得到原油黏度為15 mPa·s時(shí)不同入口速度下的油水分布情況，分析入口水速和流態(tài)之間的關(guān)系。

3計(jì)算結(jié)果和分析

3.1一定入口速度下黏度變化對(duì)波及面積的影響

通過(guò)圖2可以看到在入口速度為0.08 m/s時(shí)完整的水驅(qū)過(guò)程波及面積變化，黏度對(duì)單縫洞型油藏中水驅(qū)波及面積的影響很大：在黏度為5 mPa·s時(shí)，水相趨向中心線，向洞口流動(dòng)現(xiàn)象明顯，油藏內(nèi)形成漩渦，波及面積最小，水相占據(jù)比例為27%；當(dāng)黏度增大到15 mPa·s時(shí)，開(kāi)始向活塞驅(qū)替轉(zhuǎn)變，沒(méi)有形成漩渦，波及面積增大，水相占據(jù)比例為52%；當(dāng)黏度增大到25 mPa·s時(shí)，呈現(xiàn)水滴狀，波及面積達(dá)到最大，水相占據(jù)比例為61%。

出現(xiàn)這種現(xiàn)象的原因：在水驅(qū)入口速度為0.08 m/s、原油黏度從5 mPa·s增大到25 mPa·s時(shí)，流動(dòng)性變差，對(duì)注入水的擴(kuò)散阻力增大，相間磨損增加，水相推進(jìn)速度變小，油水界面到達(dá)出口的時(shí)間變長(zhǎng)，擴(kuò)散作用變大；同時(shí)，由于原油黏度增大，油水界面向前推進(jìn)的速度減慢，水相沿中心線縱向擴(kuò)散趨勢(shì)逐漸增強(qiáng)。因此，類(lèi)射流現(xiàn)象隨著原油黏度增加逐漸轉(zhuǎn)換成活塞驅(qū)替，漩渦逐漸消失，波及面積逐漸增大。

3.2一定入口速度下黏度變化對(duì)油水界面的影響

為研究水驅(qū)油過(guò)程中油水界面的發(fā)展形態(tài)，從上述3種黏度下分別取6幅不同時(shí)刻的流態(tài)圖，見(jiàn)圖3。

由圖3可知，未到達(dá)出口前，3種黏度下的油水界面都是按照均勻?qū)ΨQ的方式向前推進(jìn)，到達(dá)出口后，由于壓力變化，油水界面會(huì)產(chǎn)生“尖嘴狀”。當(dāng)黏度為5 mPa·s時(shí)，水流受原油阻力較小，水相優(yōu)先于中線位置通過(guò)，由于水相擴(kuò)散作用產(chǎn)生漩渦，部分原油受到漩渦波及，在漩渦作用下向右移動(dòng)，水相到達(dá)出口時(shí)間最短，耗時(shí)14 s；黏度增大到15 mPa·s時(shí)，水相受到的阻力增大，水相沿垂直中線方向擴(kuò)散趨勢(shì)增強(qiáng)，產(chǎn)生漩渦趨勢(shì)減弱，油水界面均勻向前推進(jìn)，水相到達(dá)出口的時(shí)間變長(zhǎng)，耗時(shí)17 s；黏度增大到25 mPa·s時(shí)，在到達(dá)出口前，油水界面前緣以圓形均勻前進(jìn)，靠近出口時(shí)，前緣呈現(xiàn)“尖嘴狀”，水相到達(dá)出口的時(shí)間最長(zhǎng)，耗時(shí)19 s。

3.3一定黏度下入口速度變化對(duì)波及面積的影響

由圖4可知，入口速度對(duì)波及面積影響顯著，油水界面到達(dá)右側(cè)縫洞出口時(shí)，波及面積隨著入口速度增加而減小。通過(guò)圖4中完整的水驅(qū)過(guò)程波及面積圖可知，入口速度對(duì)單縫洞型油藏中水驅(qū)波及面積的影響很大：在速度為0.08 m/s時(shí)，波及面積呈現(xiàn)水滴狀，油藏內(nèi)無(wú)漩渦產(chǎn)生，波及面積最大，水相占據(jù)比例為59%；當(dāng)入口速度為016 m/s時(shí)，水相流動(dòng)有向中心線靠近的趨勢(shì)，漩渦開(kāi)始形成，波及面積變小，水相占據(jù)比例為43%；當(dāng)入口速度為024 m/s時(shí)，射流現(xiàn)象明顯，產(chǎn)生漩渦，波及面積最小，油水界面推進(jìn)速度明顯下降。

出現(xiàn)這種現(xiàn)象的原因：在原油黏度為15 mPa·s時(shí)，水驅(qū)入口速度從0.08 m/s增加到0.24 m/s，由于注入速度與壓力梯度成正比，隨著入口速度增加，壓力梯度增大，油水界面向前推進(jìn)的速度隨之增加，水相推進(jìn)速度變大，油水界面到達(dá)出口時(shí)間變短，從18 s減少到8 s再減少到4 s，相間擴(kuò)散作用變?。煌瑫r(shí)，由于入口速度增大，水相更容易沿孔洞中間區(qū)域到達(dá)右側(cè)出口端。因此，波及面積隨著注入速度增加活塞驅(qū)替向射流發(fā)展，漩渦逐漸形成，波及面積逐漸下降。

3.4一定黏度下入口速度變化對(duì)油水界面的影響

為了研究水驅(qū)油過(guò)程中油水界面的發(fā)展形態(tài)，從上述3種速度下分別取6幅不同時(shí)刻的流態(tài)圖，見(jiàn)圖5。

由圖5可知，未到達(dá)出口前，3種速度下的油水界面都是按照均勻?qū)ΨQ的方式向前推進(jìn)，到達(dá)出口附近時(shí)，由于壓力變化，油水界面會(huì)產(chǎn)生“尖嘴狀”。當(dāng)速度為0.08 m/s時(shí)，水相優(yōu)先于中線位置通過(guò)現(xiàn)象并不明顯，油水界面到達(dá)出口時(shí)間為17 s；當(dāng)速度為0.16 m/s時(shí)，水相沿中線向右側(cè)出口推進(jìn)的趨勢(shì)增強(qiáng)，漩渦開(kāi)始產(chǎn)生，油水界面均勻向前推進(jìn)，到達(dá)出口時(shí)間為8.6 s；當(dāng)速度為0.24 m/s時(shí)，油水界面以圓形均勻前進(jìn)，到達(dá)出口時(shí)間為4 s，明顯產(chǎn)生漩渦，部分原油被卷入漩渦中。

3.5準(zhǔn)確性論證

模型采用兩套網(wǎng)格對(duì)比，對(duì)網(wǎng)格無(wú)關(guān)性進(jìn)行分析，第1套網(wǎng)格數(shù)量為16 283、節(jié)點(diǎn)數(shù)量為16 024，采用四面體結(jié)構(gòu)網(wǎng)格。第2套網(wǎng)格數(shù)量為52 717、節(jié)點(diǎn)數(shù)量為54 238，采用三角形網(wǎng)格。將兩套網(wǎng)格條件下模型的垂直中心軸線的鉛垂線上的油相體積分布進(jìn)行對(duì)比和模擬驗(yàn)證，見(jiàn)圖6。由圖6可知，平均相符程度為99.33%，可以認(rèn)為兩套網(wǎng)格的模擬結(jié)果相同，能做到網(wǎng)格無(wú)關(guān)性?；谶\(yùn)算量考慮，第1套網(wǎng)格更小，運(yùn)算更快，因此后續(xù)使用第1套結(jié)構(gòu)網(wǎng)格進(jìn)行數(shù)值模擬分析。

4結(jié)論

在單縫洞型油藏水驅(qū)開(kāi)采過(guò)程中，原油黏度和注水速度的改變會(huì)在很大程度上影響波及面積和油水界面，從而影響水驅(qū)效率。

1）相同注水速度下，黏度越大，油水界面推進(jìn)速率越低，原油黏度較低時(shí)會(huì)在出口產(chǎn)生回流，浪費(fèi)能量達(dá)不到水驅(qū)預(yù)期效果，波及面積隨著黏度增大逐漸增大，推進(jìn)速度降低。

2）原油黏度一定時(shí)，波及面積隨著入口速度的增加而減小，油水界面呈對(duì)稱趨勢(shì)均勻推進(jìn)，水速較大時(shí)會(huì)發(fā)生類(lèi)射流現(xiàn)象并有漩渦產(chǎn)生，隨著入口水速增大，類(lèi)射流現(xiàn)象和漩渦現(xiàn)象逐漸明顯。

3）由于水驅(qū)速度和原油黏度的共同影響，產(chǎn)生漩渦和類(lèi)射流現(xiàn)象，在水驅(qū)過(guò)程中造成能量浪費(fèi)，對(duì)于不同黏度原油需要選擇合理的注入水速度，達(dá)到預(yù)期水驅(qū)效果。

由于縫洞型油藏儲(chǔ)層地質(zhì)特征復(fù)雜，溶洞內(nèi)流體流動(dòng)形態(tài)受到多種因素影響。本文從流體力學(xué)角度出發(fā)研究了原油黏度和注水速度對(duì)波及面積和油水界面的影響，對(duì)內(nèi)部流場(chǎng)壓力和速度等具體分布情況缺少研究，縫洞比、儲(chǔ)層溫度、溶洞所處地層深度等因素對(duì)溶洞內(nèi)流體流態(tài)分布都有待進(jìn)一步研究。

參考文獻(xiàn)/References：

[1]任文博，陳小凡.縫洞型碳酸鹽巖油藏非對(duì)稱不穩(wěn)定注水研究[J].科學(xué)技術(shù)與工程，2013，13（27）：8120-8125.

REN Wenbo， CHEN Xiaofan. Study on asymmetric instability water filling in fractured-hole carbonate reservoirs [J]. Science and Technology， 2013，13 （27）： 8120-8125.

[2]胡蓉蓉，姚軍，孫致學(xué)，等.塔河油田縫洞型碳酸鹽巖油藏注氣驅(qū)油提高采收率機(jī)理研究[J].西安石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2015，30（2）：49-53.

HU Rongrong， YAO Jun， SUN Zhixue， et al. Research on mechanism of enhanced oil recovery by gas-flooding in fractured-cavity carbonate reservoirs in Tahe oilfield [J]. Journal of Xian Shiyou University （Natural Science Edition）， 2015， 30 （2）： 49-53.

[3]李?lèi)?ài)芬，張東，高成海.封閉定容型縫洞單元注水替油開(kāi)采規(guī)律[J].油氣地質(zhì)與采收率，2012，19（3）：94-97.

LI Aifen， ZHANG dong， GAO Chenghai.The law of water injection and oil extraction in closed volume cells[J]. Petroleum Geology and Recovery， 2012，19 （3）： 94-97.

[4]修乃嶺. 縫洞型碳酸鹽巖油藏流動(dòng)機(jī)理研究[D].廊坊：中國(guó)科學(xué)院研究生院（滲流流體力學(xué)研究所），2008.

XIU Nailing. Study on Flow Mechanism of Fractured-hole Carbonate Reservoir [D]. Langfang：Graduate School of the Chinese Academy of Sciences （Institute of Percolation and Hydrodynamics）， 2008.

[5]修乃嶺，熊偉，高樹(shù)生，等.縫洞型碳酸鹽巖油藏流動(dòng)機(jī)理初探[J].鉆采工藝，2008，31（1）：63-65.

XIU Nailing， XIONG Wei， GAO Shusheng， et al.Study on the flow mechanism of fractured-cavity carbonatite reservoir [J]. Drilling & Production Technology， 2008，31（1）： 63-65.

[6]孫玉平，修乃嶺，熊偉，等.縫洞型碳酸鹽巖油藏流動(dòng)數(shù)學(xué)模型初探[J].石油鉆探技術(shù)，2008，36（1）：65-68.

SUN Yuping， XIU Nailing， XIONG Wei， et al.Study on the flow mathematical model of fractured-cavity carbonate reservoir [J]. Petroleum Drilling Technology， 2008，36 （1）： 65-68.

[7]邢翠巧，尹洪軍，李興科，等.縫洞型碳酸鹽巖油藏試井分析模型研究[J].河北工業(yè)科技，2018，35（1）：31-36.

XING Cuiqiao，YIN Hongjun，LI Xingke，et al.Study on well test analysis model of fracture-cavity carbonate reservoir[J].Hebei Journal of Industrial Science and Technology，2018，35（1）：31-36.

[8]KANG Z J，WU Y S.Modeling multiphase flow in naturally fractured vuggy petroleum reservoirs[C]/ /The 2006 SPE Annual Technical Conference and Exhibition. Texas： SPE，2006：1-10.

[9]WU Y S，DI Y，KANG Z J，et al.A multiple-continuum model for simu-lating single-phase and multiphase flow in naturally fractured vuggy reservoirs[J].Journal of Petroleum Science and Engineering，2011，78（1）：13-22.

[10]鄭小敏，孫雷，王雷，等.縫洞型碳酸鹽巖油藏水驅(qū)油機(jī)理物理模擬研究[J].西南石油大學(xué)學(xué)報(bào)（自然科學(xué)版），2010，32（2）：89-92.

ZHENG Xiaomin， SUN Lei， WANG Lei， et al.Physical simulation of water displacing oil mechanism for vuggy fractured carbonate rock reservoir[J]. Journal of Southwest Petroleum University（Science & Technology Edition）， 2010，32 （2）： 89-92.

[11]康志宏.縫洞型碳酸鹽巖油藏水驅(qū)油機(jī)理模擬試驗(yàn)研究[J].中國(guó)西部油氣地質(zhì)，2006，2（1）：87-90.

KANG Zhihong.Simulation test study on mechanism of water flooding in fractured carbonate reservoir [J]. Petroleum Geology， West China， 2006，2（1）： 87-90.

[12]楊堅(jiān)，程倩，李江龍.塔里木盆地塔河4區(qū)縫洞型油藏井間連通程度[J].石油與天然氣地質(zhì)，2012，33（3）：484-489.

YANG Jian，CHENG Qian，LI Jianglong.Interwell-connectivity analysis for fracture-vug reservoirs in the block-4 of Tahe oilfield，Tarim Basin[J].Oil & Gas Geology，2012，33（3）：484-489.

[13]康志江，李江龍，張冬麗.塔河縫洞型碳酸鹽巖油藏滲流特征[J].石油與天然氣地質(zhì)，2005，26（5）：634-640.

KANG Zhijiang，LI Jianglong，ZHANG Dongli.Percolation characteristics of fractured-vuggy carbonate reservoir in Tahe oilfield[J].Oil & Gas Geology，2005，26（5）：634-640.

[14]康志江，趙艷艷，張冬麗.縫洞型碳酸鹽巖油藏?cái)?shù)值模擬基礎(chǔ)理論與應(yīng)用[J].大慶石油地質(zhì)與開(kāi)發(fā)，2011，30 （5）：68-71.

KANG Zhijiang，ZHAO Yanyan，ZHANG Dongli.Basic theory and application of numerical simulation in fractured-vuggy carbonate reservoir[J].Petroleum Geology and Exploitation in Daqing，2011，30（5）：68-71.

[15]李宗宇.塔河奧陶系縫洞型碳酸鹽巖油藏開(kāi)發(fā)對(duì)策探討[J].石油與天然氣地質(zhì)，2007，28（6）：856-862.

LI Zongyu.A discussion on development strategies for the Ordovician fractured-vuggy carbonate rock reservoirs in Tahe oilfield[J].Oil & Gas Geology，2007，28（6）：856-862.

[16]趙歡，尹洪軍，劉紅，等.注聚參數(shù)對(duì)L油田B區(qū)塊“二三結(jié)合”滲流特征的影響[J].河北科技大學(xué)學(xué)報(bào)，2015，36（3）：324-329.

ZHAO Huan， YIN Hongjun， LIU Hong， et al. Effects of poly condensation parameters on the seepage characteristics of “two-three combination” in block B of L oilfield[J].Journal of Hebei University of Science and Technology，2015，36（3）：324-329.

[17]劉金玉. 縫洞型介質(zhì)兩相流動(dòng)機(jī)理的數(shù)值模擬研究[D].北京：中國(guó)石油大學(xué)，2009.

LIU Jinyu.Numerical Simulation of Two-Phase Flow Mechanism in Fractured-Cavity Media [D]. Beijing：China University of Petroleum， 2009.

[18]張少峰，曹會(huì)敏，劉燕，等.彎管中液固兩相流及壁面碰撞磨損的數(shù)值模擬[J].河北工業(yè)大學(xué)學(xué)報(bào)， 2008，37（3）：48-54.

ZHANG Shaofeng， CAO Huimin， LIU Yan， et al.Numerical simulation of collision and wear between liquid-solid two-phase flow and wall in bend [J].Journal of Hebei University of Technology， 2008，37（3）： 48-54.