李澤亮, 張予生, 屈艷梅, 忽磊磊, 胡秀, 馬勇軍

(1.中國(guó)石油長(zhǎng)慶油田蘇里格南作業(yè)分公司, 陜西 西安 710018; 2.中國(guó)石油集團(tuán)測(cè)井有限公司吐哈事業(yè)部, 新疆 哈密 839009; 3.西安石油大學(xué)石油工程學(xué)院, 陜西 西安 710065)

0 引 言

溫西三塊三間房組油層處于辮狀河三角洲前緣水下分流河道及其邊前緣沉積微相帶中,儲(chǔ)集體由主體河道多次疊加形成,三角洲前緣總體受水流和波浪雙重作用,儲(chǔ)集砂體結(jié)構(gòu)特點(diǎn)復(fù)雜,注水開發(fā)水淹層測(cè)井響應(yīng)復(fù)雜多變。在水下分流河道及其河口壩發(fā)育部位,沉積物在平面上錯(cuò)綜分布,剖面上相互連通,致使儲(chǔ)層內(nèi)部高滲透碎屑沉積物在平面上呈狹長(zhǎng)條帶狀出現(xiàn),在橫剖面上高滲透層與非滲透層分流間泥質(zhì)沉積間隔分布。形成的主體河道相對(duì)高滲透層段在剖面上具有非均質(zhì)分布的多層段性,它們?cè)谧⑺_發(fā)中高滲透層段首先水淹,且水淹嚴(yán)重,形成滲流通道。油層在注入水的長(zhǎng)期沖刷下,水淹厚度向上、下2個(gè)方向逐漸擴(kuò)大,在其河道正向單滲砂層和疊置滲透層間差異層內(nèi)部水淹程度是多變的,在主體層段強(qiáng)水淹背景下,其中、上部可形成弱水淹或未水淹部位,或者在疊置高低滲油層交互狀況下,在高滲透油層干擾下造成相對(duì)低滲透層間差異弱水淹或未水淹層。特別在河道邊遠(yuǎn)梢席狀砂、遠(yuǎn)砂沉積模式在注水開發(fā)中受縱橫向組合高滲透河道層間干擾,造成許多未水淹和弱水淹層段,導(dǎo)致相對(duì)薄差油層未能有效開發(fā)利用。

該區(qū)目的層段儲(chǔ)集砂體不但受沉積時(shí)期水動(dòng)力條件、物源供應(yīng)能力、沉積速度變化及特征控制,特別受注水開發(fā)水淹程度變化的影響。水淹層注水沖刷孔隙半徑普遍加大,迂曲度減小,連通性變好,孔隙結(jié)構(gòu)趨好。它們控制著水淹儲(chǔ)層巖石物理性質(zhì)的變化,反映出儲(chǔ)層孔隙度、滲透率及其受淹測(cè)井曲線受損的幅度、厚度形態(tài)、接觸關(guān)系與圓滑程度變化。水淹層具有十分明顯的非均質(zhì)、非線性分布和測(cè)井響應(yīng)復(fù)雜的特點(diǎn)[1-4]。

為此,有必要利用該區(qū)目的層段水淹層各類測(cè)井、巖心和試油試采等資料,研究水淹層巖石物理相的多種信息,分析建立不同類別巖石物理相水淹儲(chǔ)層孔隙度、滲透率參數(shù)解釋模型。

1 水淹層分類及其評(píng)價(jià)

該區(qū)三間房組水淹層巖石物理相和測(cè)井響應(yīng)特征表明,巖石物理相分類集中體現(xiàn)出水淹層沉積成巖作用特征及其巖性、物性、孔隙類型、孔喉結(jié)構(gòu)和測(cè)井響應(yīng)對(duì)儲(chǔ)層巖石物理相的控制作用。特別在注水開發(fā)中,油層物理性質(zhì)將發(fā)生一系列的變化,在強(qiáng)水洗地帶,地層的含油性及油水分布、地層水礦化度和電阻率、黏土礦物微觀結(jié)構(gòu)、孔隙度、滲透率、潤(rùn)濕性、驅(qū)油效率、地層壓力和溫度、放射性均會(huì)發(fā)生明顯變化,而且這些變化同油層物理性質(zhì)、注入水型及其注入量等有關(guān)。不同注入期,這些變化也不同,因而導(dǎo)致地質(zhì)情況更加復(fù)雜多變[4-7]。為此,依據(jù)該區(qū)水淹層巖石物理相形成的地質(zhì)和注水開發(fā)條件,分析該區(qū)目的層段儲(chǔ)層沉積、成巖作用和巖性、物性、孔隙類型、孔喉結(jié)構(gòu)、水淹特征及測(cè)井響應(yīng)特征,在研究區(qū)劃分出好、中、差的Ⅰ、Ⅱ、Ⅲ類巖石物理相類型。

Ⅰ類中高單滲砂層巖石物理相處于淺水辮狀河三角洲前緣水下分流河道主體微相帶中,沉積中形成礫狀不等粒砂巖、中粗粒砂巖為主的長(zhǎng)石質(zhì)巖屑砂巖沉積的最粗部分,成巖中壓實(shí)、微弱膠結(jié)形成有利的強(qiáng)溶解中孔中高滲單滲砂層能量厚度儲(chǔ)層。它們?cè)谟吞镒⑺_發(fā)中基本上處于滲流通道主峰附近,孔隙結(jié)構(gòu)及其孔隙度、滲透率更趨增大,形成滲流通道附近強(qiáng)水淹層及較強(qiáng)水淹層。測(cè)井響應(yīng)相對(duì)于泥巖呈現(xiàn)“三大減三小減一高一變”特征,即自然電位、自然伽馬和密度減小大;聲波時(shí)差、中子孔隙度、井徑呈微減或微縮,變化穩(wěn)定;徑向電阻率不存在低電阻率環(huán)帶及正差異,數(shù)值高并有變化趨勢(shì)。反映出一種強(qiáng)水淹中高單滲砂層能量巖石物理相流動(dòng)單元。

Ⅱ類層間差異滲砂層巖石物理相處于河道高滲主體微相干擾下河道次主體微相帶中,沉積中在主流線之外形成的中、細(xì)砂巖大多處于河道單滲砂層中上部或河道疊置滲層間干擾差異層內(nèi),成巖中強(qiáng)壓實(shí)、較差膠結(jié)形成中溶解中孔隙度中低滲透率層間差異砂層能量厚度儲(chǔ)層。它們?cè)谟吞镒⑺_發(fā)中處于滲流通道側(cè)翼或者邊緣,孔隙結(jié)構(gòu)及其孔隙度、滲透率總體趨大,但不明顯,形成滲流通道側(cè)翼的中水淹層及邊遠(yuǎn)區(qū)的弱水淹層。測(cè)井響應(yīng)相對(duì)于泥巖呈現(xiàn)“六減小一較高”特征,即自然電位、自然伽馬、密度減小較大,聲波時(shí)差、中子孔隙度、井徑減低或微縮;徑向電阻率低阻環(huán)帶及正差異不明顯,數(shù)值較高。反映出一種中水淹及中弱水淹單滲砂能量巖石物理相流動(dòng)單元。

Ⅲ類低滲砂層巖石物理相處于河道邊遠(yuǎn)梢席狀砂、遠(yuǎn)砂微相帶中,沉積中受湖浪影響巖性細(xì),形成以中細(xì)砂巖、細(xì)砂巖及粉砂巖組成厚度薄的席狀或非席狀砂體,成巖中強(qiáng)壓實(shí)和膠結(jié)作用形成低孔隙度低滲透率席狀砂、遠(yuǎn)砂剩余能量?jī)?chǔ)層。它們?cè)谟吞镒⑺_發(fā)中處于非主力油層單滲砂層能量厚度小,注水開發(fā)中基本上不形成滲流通道,儲(chǔ)層處于弱水淹或未水淹狀態(tài)。測(cè)井響應(yīng)相對(duì)于泥巖呈現(xiàn)“五減小一升高不變”特征,即自然電位、自然伽馬、密度減小較小;聲波時(shí)差、中子孔隙度、井徑更減低或不明顯;徑向電阻率呈低電阻率環(huán)帶及正差異,數(shù)值有所降低。特別在注水驅(qū)油使儲(chǔ)層中黏土受注水浸泡發(fā)生膨脹,孔喉變窄,儲(chǔ)層孔隙結(jié)構(gòu)、孔隙度、滲透率特征及其測(cè)井響應(yīng)較差,反映出一種弱水淹及未水淹低滲砂層剩余能量巖石物理相流動(dòng)單元。

上述3類不同類別巖石物理相處于不同沉積成巖儲(chǔ)集相帶及注水開發(fā)環(huán)境中,儲(chǔ)層具有不同巖性、物性、孔隙圖像、孔隙結(jié)構(gòu)及水淹特征,反映出不同類別水淹儲(chǔ)層參數(shù)和測(cè)井響應(yīng)特征與差異。為此,利用水淹儲(chǔ)層巖石物理相分類參數(shù)和測(cè)井響應(yīng)曲線特征,采用物性標(biāo)準(zhǔn)與自然電位、深電阻率、自然伽馬、井徑、密度、聲波時(shí)差、中子孔隙度、徑向電阻率差異及儲(chǔ)層質(zhì)量指標(biāo)(見表1),利用灰色理論水淹儲(chǔ)層評(píng)價(jià)方法進(jìn)行被評(píng)價(jià)井點(diǎn)數(shù)據(jù)的綜合分析處理。采用矩陣分析、標(biāo)準(zhǔn)化、標(biāo)準(zhǔn)指標(biāo)絕對(duì)差的極值加權(quán)組合方法技術(shù),利用灰色理論集成和綜合多種信息,確定和劃分出該區(qū)目的層段Ⅰ、Ⅱ、Ⅲ類(中高單滲透砂、層間差異滲透砂、低滲透砂)巖石物理相儲(chǔ)層[7-10]。

表1 研究區(qū)三間房組水淹層巖石物理相評(píng)價(jià)劃分標(biāo)準(zhǔn)表

*非法定計(jì)量單位,1 mD=9.87×10-4μm2,下同

2 水淹層參數(shù)分類建模

不同類別巖石物理相水淹層滲透率、孔隙度主要取決于組成巖石顆粒大小、孔隙類型、組織結(jié)構(gòu)、水淹狀況及所含流體性質(zhì)等,明顯反映出不同類別巖石物理相儲(chǔ)層參數(shù)分布概率的差異。Ⅰ類中高單滲砂層巖石物理相巖性以礫狀不等粒砂巖、中粗粒砂巖為主,粒級(jí)較粗,孔隙類型主要為殘余粒間孔和較大粒間溶孔組合為主,注水開發(fā)中處于滲流通道主峰,形成中高單滲強(qiáng)水淹層,具較好的物性及孔隙結(jié)構(gòu)特征,儲(chǔ)層滲透率、孔隙度參數(shù)(含相應(yīng)測(cè)井響應(yīng)參數(shù))的分布及標(biāo)準(zhǔn)趨于相對(duì)集中的高值范圍。而Ⅲ類低滲砂層巖石物理相巖性以細(xì)砂巖及粉砂巖組成厚度薄的席狀或非席狀砂體,孔隙類型以粒間孔和微孔隙為主,注水開發(fā)中基本不形成滲流通道,儲(chǔ)層處于弱水淹或未水淹狀態(tài),物性及孔隙結(jié)構(gòu)差,儲(chǔ)層參數(shù)滲透率、孔隙度(含相應(yīng)測(cè)井響應(yīng)參數(shù))的分布及標(biāo)準(zhǔn)趨于相對(duì)集中的低值范圍。Ⅱ類層間差異層巖石物理相儲(chǔ)層的巖性、物性、水淹狀況、孔隙類型與結(jié)構(gòu)的分布及標(biāo)準(zhǔn)則趨于居中。特別Ⅰ、Ⅱ類強(qiáng)中水淹層中,由于注入水的沖刷,巖石孔壁上貼附的黏土被剝落,含油砂巖較大孔隙中的黏土被沖散、沖走,溝通孔隙的喉道半徑加大,孔隙變得干凈、暢通,孔隙半徑普遍增大,迂曲度減小、連通性變好,實(shí)際上縮短了流體滲流途徑,巖石孔隙結(jié)構(gòu)明顯趨好,因而水淹層孔隙度、滲透率均有不同程度增大變化。尤其是水淹后孔隙喉道比例增大,并在水洗中產(chǎn)生微裂縫,造成中強(qiáng)水淹層滲透率有更大幅度增高變化。因此,必須研究巖石物理相分類才能準(zhǔn)確進(jìn)行水淹層孔隙度、滲透率參數(shù)建模[10-12]。

在該區(qū)水淹層巖石物理相分類評(píng)價(jià)基礎(chǔ)上,采用分類巖心刻度測(cè)井方法更為有效地區(qū)分和建立測(cè)井水淹層參數(shù)解釋模型。

圖1至圖6為Ⅰ類巖石物理相巖心分析孔隙度、滲透率、聲波時(shí)差、中子孔隙度、密度測(cè)井響應(yīng)擬合關(guān)系圖,圖中數(shù)據(jù)點(diǎn)的均勻程度及擬合線性關(guān)系具有相對(duì)集中趨勢(shì),利用相應(yīng)統(tǒng)計(jì)擬合關(guān)系都可以較為準(zhǔn)確表達(dá)和建立孔隙度、滲透率參數(shù)解釋模型。

圖1 Ⅰ類巖石物理相巖心分析孔隙度與聲波時(shí)差關(guān)系圖

圖4 Ⅰ類巖石物理相聲波時(shí)差與巖心分析滲透率關(guān)系圖

圖5 Ⅰ類巖石物理相中子孔隙度與巖心分析滲透率關(guān)系圖

圖6 Ⅰ類巖石物理相密度與巖心分析滲透率關(guān)系圖

表2為該區(qū)水淹層不同巖石物理相類別建立的測(cè)井水淹層孔隙度、滲透率參數(shù)解釋模型。

表2中計(jì)算參數(shù)ρb、Δt、φN、α、GR1分別為密度、聲波時(shí)差、中子孔隙度測(cè)井值及自然電位、自然伽馬減小系數(shù),它們的減小系數(shù)為測(cè)井響應(yīng)值相對(duì)于泥質(zhì)巖的減小程度。表2中巖石物理相分類后以聲波時(shí)差、中子孔隙度與巖心分析建立的孔隙度模型擬合度為高,反映孔隙度參數(shù)模型趨于較好;滲透率模型擬合難度增大,以孔隙度、聲波時(shí)差、中子孔隙度和密度擬合的滲透率趨于較高[12-14]。

表2 研究區(qū)分類巖石物理相水淹層孔隙度、滲透率參數(shù)解釋模型

3 水淹層參數(shù)精細(xì)建模方法

3.1 孔隙度參數(shù)

可以利用上述聲波時(shí)差、中子孔隙度分類模型及其相應(yīng)聲波時(shí)差、中子孔隙度交會(huì)模型求取孔隙度,更為準(zhǔn)確地利用綜合信息求取分類巖石物理相孔隙度參數(shù)值[15-16]。

聲波時(shí)差測(cè)井響應(yīng)值計(jì)算孔隙度模型

φ1=A1Δt+B1

(1)

中子孔隙度測(cè)井響應(yīng)值計(jì)算孔隙度模型

φ2=A2φN+B2

(2)

分別制作聲波時(shí)差與中子孔隙度測(cè)井值交會(huì)關(guān)系圖(見圖7),得到相關(guān)公式

Δt=AφN+B

(3)

圖7 聲波時(shí)差、中子孔隙度求取孔隙度交會(huì)關(guān)系圖

圖7中,計(jì)算層點(diǎn)孔隙度對(duì)應(yīng)聲波時(shí)差、中子孔隙度分別為Δti、φNi,對(duì)應(yīng)到聲波時(shí)差與中子孔隙度關(guān)系圖中分別為F與E點(diǎn),顯然這2點(diǎn)不在Δt=AφN+B直線上,從而造成了分別采用聲波時(shí)差和中子孔隙度計(jì)算的孔隙度相差較大。采用該層點(diǎn)聲波時(shí)差、中子孔隙度對(duì)應(yīng)圖中交會(huì)直線E、F中點(diǎn)M,把該中點(diǎn)M作為綜合后層點(diǎn)計(jì)算點(diǎn),分別采用聲波時(shí)差和中子孔隙度模型計(jì)算孔隙度,求取兩者孔隙度參數(shù)綜合擬合值為其更可靠的孔隙度值[15-16]。

則該層點(diǎn)孔隙度值為

φi=φ1i+φ2i2=A14(Δti+AφNi+B)+

A24φNi+Δti-BA+B1+B22

(4)

利用不同巖石物理相聲波時(shí)差、中子孔隙度交會(huì)統(tǒng)計(jì)模式,建立聲波時(shí)差、中子孔隙度綜合擬合計(jì)算孔隙度參數(shù)精細(xì)解釋模型(見表3)。

3.2 滲透率參數(shù)

表1中6種方法確定6個(gè)滲透率參數(shù),仍然有可能包含若干有利條件和不利因素,因而很難避免數(shù)值上有較大誤差。為此,利用Hodges-Lehmann方法,分析兼顧各個(gè)評(píng)估數(shù)值,較大限度減小偏離很大的測(cè)量結(jié)果影響[15-16]。首先計(jì)算出所有二元對(duì)

K(1)+K(1)2,K(1)+K(2)2,…,…,K(1)+K(m)2

K(2)+K(2)2,K(2)+K(3)2,…,K(2)+K(m)2

K(3)+K(3)2,…,K(3)+K(m)2

…,

K(m)+K(m)2

對(duì)6個(gè)數(shù),m=6,二元對(duì)的總數(shù)為m+(m-1)+…+1=21。其Hodges-Lehmann綜合估計(jì)值則為這21個(gè)數(shù)的中間數(shù)值。如果二元對(duì)是一組偶數(shù)的數(shù)值,Hodges-Lehmann綜合估計(jì)值為數(shù)值排列2個(gè)中間數(shù)值的算術(shù)平均值。從而采用分類分析的Hodges-Lehmann方法統(tǒng)計(jì)改善提高滲透率參數(shù)計(jì)算精度。

表3 分類巖石物理相聲波時(shí)差、中子孔隙度綜合擬合計(jì)算孔隙度參數(shù)精細(xì)解釋模型

圖8 溫西3-2343井層段水淹層巖石物理相分類及其參數(shù)計(jì)算成果圖

4 水淹層測(cè)井評(píng)價(jià)解釋實(shí)例

滲透率則采用該層段Ⅰ類巖石物理相孔隙度、聲波時(shí)差、中子孔隙度、自然電位減小系數(shù)、自然伽馬減小系數(shù)測(cè)井解釋滲透率模型,分別計(jì)算出102.83×10-3、137.59×10-3、158.34×10-3、125.89×10-3、113.29×10-3μm2等5個(gè)滲透率值,計(jì)算出所有二元對(duì)

102.83 120.21 130.60 114.36 108.06

137.59 147.98 131.74 125.44

158.34 142.13 135.83

125.89 119.59

113.29

剩余油飽和度參數(shù)采用層段深度水樣礦化度及地層溫度求取混合液電阻率Rz=0.55 Ω·m、深探測(cè)電阻率Rt=19.5 Ω·m、孔隙度φ=18.39%及a、b、m、n值,利用阿爾奇公式求得剩余油飽和度So=20.91%。

同樣計(jì)算上部2 374.5～2 380.5 m段的2號(hào)層段Ⅱ類巖石物理相中水淹層孔隙度為15.93%,滲透率為29.66×10-3μm2,剩余油飽和度為31.46%。

上述1、2號(hào)層段計(jì)算的孔隙度、滲透率、剩余油飽和度參數(shù)分別反映Ⅰ類巖石物理相強(qiáng)水淹和Ⅱ類巖石物理相中水淹層的水淹狀況(見表1),有效地改善和提高了水淹層參數(shù)計(jì)算的精度。

經(jīng)在中部2 382.5～2 387.0 m段的Ⅰ類巖石物理相1層段強(qiáng)水淹層射孔,日產(chǎn)油0,日出水29.3 m3,含水100%,證實(shí)了水淹層巖石物理相分類及其參數(shù)解釋的有效性。

5 結(jié) 論

(1) 該區(qū)塊三間房組油層處于辮狀河三角洲前緣水下分流河道及其邊前緣沉積微相帶中,受其水流和波浪雙重作用,儲(chǔ)集砂體結(jié)構(gòu)復(fù)雜,注水開發(fā)水淹層測(cè)井響應(yīng)特征多變,儲(chǔ)層具有非均質(zhì)、非線性分布特點(diǎn)。通過巖石物理相分類進(jìn)行水淹層參數(shù)精細(xì)建模,實(shí)現(xiàn)在同類注水開發(fā)條件水淹層及其巖性、物性、孔隙結(jié)構(gòu)、測(cè)井響應(yīng)特征的巖相建立同類測(cè)井參數(shù)模型,將水淹層非均質(zhì)、非線性問題轉(zhuǎn)化為相對(duì)均質(zhì)、線性問題。

(2) 巖石物理相分類建立水淹層參數(shù)模型研究中,分類模型數(shù)據(jù)點(diǎn)分布的擬合具有相對(duì)集中分布趨勢(shì)及其較好線性關(guān)系。特別是分類模型分別利用聲波時(shí)差、中子孔隙度參數(shù)模型的綜合擬合值求取有效孔隙度參數(shù);滲透率分類模型采用孔隙度、聲波、中子、密度、自然電位、自然伽馬等多種曲線分別求取,并利用Hodges-Lehmann法綜合評(píng)估擬定。從而通過水淹層精細(xì)建模提高了水淹參數(shù)計(jì)算精度,為準(zhǔn)確計(jì)算水淹層參數(shù)提供了有效方法。

參考文獻(xiàn):

[1] 李菊花, 李相方, 高文軍, 等. 吐哈油田中高含水期油藏剩余油分布研究 [J]. 天然氣地球科學(xué), 2005, 16(3): 378-381.

[2] 高博禹, 彭仕宓, 王建波, 等. 剩余油形成與分布的研究現(xiàn)狀及發(fā)展趨勢(shì) [J]. 特種油氣藏, 2004, 11(4): 7-11.

[3] 宋子齊, 趙磊, 王瑞飛, 等. 利用常規(guī)測(cè)井方法識(shí)別劃分水淹層 [J]. 西安石油學(xué)院學(xué)報(bào), 2003, 18(6): 50-53.

[4] 張斌成, 石曉燕, 劉瑛, 等. 水淹層測(cè)井綜合解釋及水淹特征研究 [J]. 測(cè)井技術(shù), 2005, 29(6): 545-547.

[5] 高興軍, 宋子齊, 程仲平, 等. 影響砂巖油藏水驅(qū)開發(fā)效果的綜合評(píng)價(jià)方法 [J]. 石油勘探與開發(fā), 2003, 30(2): 68-71.

[6] 趙培華. 油田開發(fā)水淹層測(cè)井技術(shù) [M]. 北京: 石油工業(yè)出版社, 2003: 106-185.

[7] 朱衛(wèi)紅, 韓芳芳, 章成廣, 等. 輪南油田剩余油飽和度綜合解釋方法研究 [J]. 石油天然氣學(xué)報(bào), 2008, 30(3): 257-259.

[8] 熊琦華, 彭仕宓, 黃述旺, 等. 巖石物理相研究方法初探——以遼河冷凍-雷家地區(qū)為例 [J]. 石油學(xué)報(bào), 1994, 15(增刊): 68-75.

[9] 姚光慶, 趙彥超, 張森龍. 新民油田低滲細(xì)粒儲(chǔ)集砂巖巖石物理相研究 [J]. 地球科學(xué), 1995, 20(3): 355-360.

[10] 宋子齊. 灰色理論油氣儲(chǔ)層評(píng)價(jià) [M]. 北京: 石油工業(yè)出版社, 1995: 17-77.

[11] 張俊峰, 宋自恒, 王軍, 等. 馬廠油田油層水淹級(jí)別劃分方法與標(biāo)準(zhǔn) [J]. 石油天然氣學(xué)報(bào), 2006, 28(2): 79-81.

[12] 孫寶佃, 宋子齊, 成志剛, 等. 特低滲透儲(chǔ)層測(cè)井系列優(yōu)化評(píng)價(jià)研究 [J]. 石油地球物理勘探, 2012, 47(3): 483-490.

[13] 何勝林, 高楚橋, 陳嶸. 潿洲W油田中塊剩余油飽和度計(jì)算中關(guān)鍵參數(shù)確定方法研究 [J]. 石油天然氣學(xué)報(bào), 2011, 33(12): 88-91.

[14] 高興軍, 宋新民, 褚人杰, 等. 動(dòng)靜結(jié)合方法計(jì)算儲(chǔ)層水淹后地層混合液電阻率技術(shù)及其在剩余油飽和度解釋中的應(yīng)用 [J]. 地學(xué)前緣, 2005, 15(1): 160-167.

[15] 姜香云, 王志章, 吳勝和. 儲(chǔ)層三維建模及其在油藏描述中的應(yīng)用研究 [J]. 地球物理學(xué)進(jìn)展, 2006, 21(3): 902-908.

[16] 宋子齊. 測(cè)井多參數(shù)的地質(zhì)應(yīng)用 [M]. 西安: 西北工業(yè)大學(xué)出版社, 1993: 110-140.