碳酸鹽巖孔隙性預測方法研究

魯小飛

?

碳酸鹽巖孔隙性預測方法研究

魯小飛

（中國地震局地震研究所 中國地震局地震大地測量重點實驗室，武漢 430071）

碳酸鹽巖儲層由于其結構復雜、孔隙性多樣，儲層非均質性強，常規巖石物理模型，或經驗模型難以準確描述儲層的孔隙性結構和彈性參數。用Gassmann方程的推廣分析碳酸鹽巖中不同孔隙相礦物組分對彈性參數的影響，研究碳酸鹽巖中白云巖孔隙性的對巖石物性的整體影響，用理論模型分析實際工區。研究結果表明Gassmann的推廣模型具有較高的實用性，可以為碳酸鹽巖儲層孔隙性預測提供一種新的思路。

碳酸鹽巖；孔隙性；巖石物理模型；預測方法

碳酸鹽巖通常由幾種礦物組分構成其巖石基質，巖石中存在不同大小、形狀、連通或非連通的孔隙，孔隙中往往充填著不同性質的物質（如流體或碎屑物），構成一種非均質性很強的復合介質，儲層孔隙空間及孔隙結構十分復雜具有較強非均質性。儲集空間既有孔隙也有裂縫，其中孔隙包括原生孔隙（如晶間微孔及粒間微孔隙）和次生孔隙（主要包括溶蝕孔、洞），裂縫包括微裂縫、構造縫和溶蝕縫等。由于碳酸鹽巖性脆易溶，導致在不同的地區碳酸鹽巖儲層主要儲集空間類型差異較大，有的以裂縫型為主、有的以溶孔洞型為主。因此,對碳酸鹽巖儲層評價需建立在對巖石孔隙系統的充分了解之上。

國內外許多研究已經證明碳酸鹽巖中的孔隙系統對巖石的彈性性質影響較大。目前孔隙度預測方法太多，但主要分為兩大類,一種是線性方法，主要采用經驗公式用測井曲線對孔隙度進行線性或非線性擬合，這種方法適用于簡單的地質條件。對于復雜的地質條件，許多研究者也提出了一些改進的關系式（如：Raymer等、Hamilton等、Han等、Kahraman等）。另一種是非線性方法，有人工神經網絡[1]、數據組合處理方法、隨機模擬[1]等多種方法。這類方法人工干預較小，不考慮實際巖石模型中的內在因素的影響，比如孔隙度大小，孔隙形狀。因此結果分析中或多或少存在一定的問題。該研究根據碳酸鹽巖孔隙性因素對巖石物理、彈性參數的影響分析，孔隙介質的孔隙相出發，建立了孔隙介質巖石物理模型，并分析了其在實際工區的應用。

1 理論基礎

當孔隙介質是由兩個孔隙相組合而成，每個孔隙相都可以用常規Gassman公式[2]單獨描述時，推廣的Gassman公式（Berryman和Milton,1999）描述出靜態或低頻流體填充孔隙介質的有效體積模量，這是對常規Gassmann公式的改進。常規Gassmann方程假設孔隙介質是由單一的、統計均勻的孔隙構成的，具有單一孔隙空間和單一固體礦物。像Gassmann公式一樣，推廣的Gassmann公式[2]與孔隙空間的幾何形態無關。該模型假設兩孔隙組成部分在接觸點處相連接并且充滿孔隙介質的整個空間，而且他假設頻率足夠低使得粘性和慣性作用可忽略不計[2]，而且由應力誘發的孔隙壓力在每個孔隙組分中是均勻的，盡管孔隙組分之間可能會不相同，但是當巖石的巖性由主要有一種主成分構成時，我們可以把孔隙壓力看作是均勻的。

其中，K是流體體積模量，常量有效固體的體積模量K和有效流體的體積模量K可由兩個構成孔隙介質的孔隙相的模量表示。產生這些結果的關鍵思想是，每當兩個標量場，例如P和P，可以在只有兩個組分的線性混合中獨立地變化時，存在一個遞變比的特殊值p/，該值對應于介質的總體膨脹或收縮但沒有任何相對形狀的變化，這就保證存在一套一致的關系式，使得K和K可以用巖石的骨架模量和組分模量表示。

圖1 白云巖含量與縱橫波速度比關系圖（三條曲線表示孔隙度不同）

2 理論模型分析

碳酸鹽巖的彈性參數主要受巖石的礦物組分、孔隙度、流體性質、溫度、壓力及其他影響地震波傳播的因數影響[3]。通常情況下，碳酸鹽巖的礦物組分主要是白云巖、灰巖、泥巖及其過渡成分的組成。在碳酸鹽巖的礦物成分比較簡單的情況下，地質背景并不復雜的情況下主要受孔隙度和礦物組分的影響。

圖2 孔隙度與縱橫波速度比關系圖（三條曲線表示不同白云巖和灰巖含量）

圖3 孔隙度與縱橫波速度比關系圖(表示孔隙縱橫比)

根據上述模型建立的彈性參數模型圖，一般情況下，碳酸鹽巖的比較致密，礦物由主要的一種或兩種礦物組成（如：白云巖、灰巖和泥質等）[4]，礦物含量表征彈性參數變化規律，然后在分析礦物孔隙的特征，進而建立孔隙性特征到彈性參數的變化規律。由圖1可以看出，隨著白云巖含量的增加，縱橫波速度比減小而縱波阻抗逐漸增大。同一白云巖含量，孔隙度大的縱波阻抗較小，但是縱橫波速度比較大。這是由于橫波速度在巖石骨架中傳播，孔隙度對橫波速度的影響要比縱波速度的影響小的多，所以隨著孔隙度的增加，縱橫波速度比逐漸增大（圖2），而縱波阻抗之間逐漸減小，這是必然的。而相同孔隙度，白云巖含量較多的巖石縱橫波速度比小，縱波阻抗較大，這和白云巖的性質有關[5]（白云巖的體積模量和剪切模量較大）。在孔隙度較小的情況下，縱橫波速比變化不是很大，在孔隙度下，縱橫波大速度比陡増。

碳酸鹽巖的孔隙系統遠比碎屑巖的復雜，碎屑巖中的孔隙主要屬粒間孔隙類型，而碳酸鹽巖卻包含多種類型孔隙[5]，如印?？紫?、孔洞、粒間孔隙和粒內孔隙，且巖石組分也存在差異。因此我們研究碳酸鹽巖的孔隙性要考慮巖石的孔隙結構（孔隙大小，形狀，分布情況等），理論上相同孔隙度和礦物組分的情況下，圓形孔隙的彈性模量要比扁形或裂隙形空隙大（如圖3）。實際情況中，可分析具體礦物中較易形成的孔隙類型，探求碳酸鹽巖的巖石物理性質。

圖4 工區測井綜合解釋圖

3 實例分析

四川盆地龍女寺地區為典型的碳酸鹽巖儲層，儲層的沉積類型為屬陸表海碳酸鹽臺地。礦物成分主要以白云巖，灰巖為主，中間夾雜有少量的泥和砂。前期解釋的儲層巖性主要以殘余生屑白云巖、泥晶白云巖為主。巖石組合也比較多樣，孔隙度較低，并且多以巖溶縫洞型孔隙為主。由上述模型的巖石物理量版圖1可知道，白云巖含量大時彈性參數縱波阻抗大，縱橫波速比較小。由圖4可以看出，以灰巖為主的地層比以白云巖為主的地層的縱橫波時差大，從而驗證了Vp和Vs和白云巖含量在一定程度上有一定的相關性。

我們可以遵照理論模型圖，把縱橫波速度比和λ*ρ的交匯作為區分巖性的參數，圖5和圖6為研究區磨溪29井的這兩參數的交會圖。我們把礦物含量方解石在75%的作為灰巖，白云石含量在75%以上的作為白云巖，剩下的作為泥巖處理。由圖5，可知白云巖和灰巖用縱橫波速度比和λ*ρ的交匯可以很好的區分開來，而泥巖的縱橫波速度及λ*ρ都很小，并且泥巖和白云巖有部分重疊。本工區主要為白云巖儲層，白云巖儲層孔隙度的變化范圍在0~10%左右，如圖6所示，白云巖層孔隙度分布在0~10%之間，而白云巖含量低的地層孔隙度大多都在4%以下。

圖5 磨溪29井λ*ρ與縱橫波速度比交會圖

圖6 孔隙度與縱橫波速度比交會圖

由于碳酸鹽巖成巖過程的復雜性包括白云化、溶蝕、破裂、壓實作用[7]，還包括膠結、壓實和填充等成巖作用，導致碳酸鹽巖孔隙結構復雜，在相同孔隙度的情況下，巖石的縱橫波速度也相差較大。所以要研究儲層的孔隙性，孔隙結構的分析必不可少。碳酸鹽巖地層孔隙中常見的孔隙類型有晶間孔、粒間孔、粒內孔、印?？?、溶孔、格狀孔以及裂縫等。其中，晶間孔與粒間孔的孔隙形狀一般不規則，粒內孔常出現在單一顆粒或其他開放有機物骨架中，印?？壮Ｓ梢恍┯袡C沉積物溶解形成，溶孔則由易溶解物質形成，裂縫往往由外力差異（構造應力變化、孔洞坍塌、異?？紫秹毫ψ兓纫蛩兀┮稹＿@些孔隙類型中，印?？?、裂縫以及部分溶洞為次生孔隙。粒內孔、印?？住⑷芏匆约案駹羁资芰Σ灰鬃冃危芽p最容易變形。

圖7 孔隙縱橫比反演圖

我們用DEM模型[2]反演得到研究區29井的孔隙縱橫比，來表征孔隙結構，進一步驗證本模型對孔隙的實用性。根據圖7展示的孔隙縱橫比看出該井的值主要集中在0.2~0.3之間，粒間孔，裂縫比較發育。結合測井解釋圖84 640~4 090m（白云巖地層）的深度段的孔隙縱橫比要比4 600~4 630m（灰巖地層）高，由于白云巖和灰巖的礦物性質不同，白云巖的脆性較大，易形成裂縫等孔隙縱橫比小的孔隙。這為我們研究碳酸鹽巖中礦物孔隙提供了一定的現實依據。

4 結論與認識

通過對碳酸鹽巖儲層巖石物理特征綜合分析，我們提出了碳酸鹽巖的彈性參數計算模型，本模型的計算與常規均勻介質等效模型相比，應用此模型分析碳酸鹽巖中不同孔隙相礦物組分對彈性參數的影響，研究碳酸鹽巖中白云巖孔隙性的對巖石物性的整體影響。從巖石礦物的角度分析了碳酸鹽巖的孔隙性特征。通過實際工區分析，本文方法有一定的實用性，為以后碳酸鹽巖巖石物理分析提供了新的思路。

[1] 陳繼華, 肖亮, 毛志強. 溶孔洞型含氣碳酸鹽巖儲層孔隙度評價方法研究[J]. 石油天然氣學報, 2011, 33(6): 92-97.

[2] 馬沃可, 慕克吉, 德沃爾金, 等. 巖石物理手冊: 孔隙介質中地震分析工具[M]. 中國科學技術大學出版社, 2008.

[3] 劉開元. 一種新的基于孔隙結構的碳酸鹽巖孔隙度反演方法在川東地區的應用[J]. 科學技術與工程, 2013 (33): 9802-9807.

[4] Berryman J G, Milton G W. Exact results for generalized Gassmann's equations in composite porous media with two constituents[J]. Geophysics, 1991, 56(12): 1950-1960.

[5] Mazzullo S J. Geochemical and neomorphic alteration of dolomite: a review[J]. Carbonates and evaporites, 1992, 7(1): 21-37.

[6] Ruiz F, Dvorkin J. Predicting elasticity in nonclastic rocks with a differential effective medium model[J]. Geophysics, 2010, 75(1): E41-E53.

[7] Cai C, Xie Z, Worden R H, et al. Methane-dominated thermochemical sulphate reduction in the Triassic Feixianguan Formation East Sichuan Basin, China: towards prediction of fatal H2S concentrations[J]. Marine and Petroleum Geology, 2004, 21(10): 1265-1279.

Study of the Method of Porosity Prediction for Carbonate Rock

LU Xiao-fei

(Institute of Seismology, China Earthquake Administration, Wuhan 430071)

Conventional rock physical model or experimental model is difficult to accurately describe porosity structure and elastic parameters of the carbonate rock reservoir due to complicated texture, various porosity and strong heterogeneity. This paper has a discussion on the influence of different pore phase mineral components on the elastic parameters of carbonate rocks and the influence of the porosity of dolomite on physical property of the reservoir by the use of promotion of the Gassmann equation. The study results indicate that Gassmann's extension model has a high practical value.

carbonate rock; reservoir; porosity; rock physical model; prediction method

2018-05-21

魯小飛（1987-），河南駐馬店人，實習研究員，研究方向：地球物理學

P583

A

1006-0995（2019）01-0011-04

10.3969/j.issn.1006-0995.2019.01.003