用J函數法求取碳酸鹽巖儲層飽和度方法探討

司馬立強，李 清，楊 毅，2，陳 強，3

（1.西南石油大學地球科學與技術學院，成都610500；2.中海石油（中國）有限公司湛江分公司，廣東湛江524057；3.中國石油塔里木油田分公司，新疆庫爾勒841000）

用J函數法求取碳酸鹽巖儲層飽和度方法探討

司馬立強1，李 清1，楊 毅1，2，陳 強1，3

（1.西南石油大學地球科學與技術學院，成都610500；2.中海石油（中國）有限公司湛江分公司，廣東湛江524057；3.中國石油塔里木油田分公司，新疆庫爾勒841000）

準確求取非均質碳酸鹽巖儲層含水飽和度一直是測井解釋工作的一大難題（尤其是在缺乏密閉取心等直接確定巖心含水飽和度資料的情況下）。根據毛管壓力曲線、物性分析數據及RFT測試等資料，并以中東某油田為例，探討了利用毛管壓力曲線計算非均質碳酸鹽巖儲層原始含水飽和度的方法（J函數法）。利用此方法計算出的原始含水飽和度與相滲實驗分析結果及試油結果均吻合較好。此次研究為缺乏巖心含水飽和度資料的非均質碳酸鹽巖儲層含水飽和度的求取提供了一種合理、可靠的計算方法。

碳酸鹽巖儲層；毛管壓力曲線；J函數法；原始含水飽和度

0 引言

原始含水飽和度的確定是三維儲層地質建模的一大關鍵步驟，其對后續油藏動態評價、儲量計算及開發方案的編制等均具有重要意義［1-2］。準確求取非均質碳酸鹽巖儲層含水飽和度一直是測井解釋工作的一大難題（尤其是在缺乏密閉取心等直接確定巖心含水飽和度資料的情況下）。以往儲層的含水飽和度大多是根據阿爾奇公式或者其衍生的公式確定，且主要受電阻率和孔隙度的影響。由于影響電阻率的因素很多，因此利用其公式計算的含水飽和度值，在沒有考慮其他影響因素的情況下容易造成計算結果存在較大誤差［3］。研究發現，含水飽和度與巖石的微觀結構有關，主要是巖石毛管壓力函數［2］。利用毛管壓力曲線確定原始含水飽和度已成為一種常用的手段，具體方法有：含油高度法、Purcell法和J函數法等。這3種方法都是以求取J函數為基礎，其中含油高度法偏重理論，僅通過含油高度便可計算出含水飽和度，但該方法對地層的影響因素考慮不夠；Purcell法和J函數法比含油高度法多加入了孔隙度和滲透率等參數，彌補了含油高度法的缺陷，但J函數法比Purcell法的計算過程更加簡單［4］。因此，筆者僅探討采用J函數法求取研究區的含水飽和度。

筆者以中東某油田為例，探討J函數法在非均質碳酸鹽巖儲層飽和度模型建立中的應用。根據實驗提供的3口取心井的27塊巖樣的物性分析資料，對J函數曲線進行分類，擬合實驗條件下各類J函數與巖心標準化飽和度的關系式，接著通過油藏條件下的毛管壓力、潤濕角及界面張力等參數計算J函數值，進而確定儲層中任意一點的含水飽和度。根據飽和度的計算結果與試油資料、相滲資料等的對比分析，最終對其解釋結果進行評價。

1 方法原理

由于儲層具有非均質性，任何一條毛管壓力曲線都只能描述儲層中取樣點的特征，而無法用于反映整體油藏的特征［5］。因此，Leverett通過建立毛管壓力曲線與儲層巖石物理性質之間的關系，利用J函數對毛管壓力曲線上的毛管壓力值（Pc）進行轉換，使得各巖樣的資料點相對集中，進而反映儲層的各項特征［2，5-7］。利用J函數對毛管壓力數據進行處理，其定義式為

在式（1）中可看出：J函數與巖樣的Pc，σ，θ，K和φ有關，與電阻率值無關。將求得的J函數與實驗條件下的Swn擬合，得出Swn與J函數的關系式。以Swn為縱坐標、J函數為橫坐標建立兩者的關系式。常用的擬合函數有冪函數、對數函數以及指數函數等［8-11］。根據研究區毛管壓力曲線的變化規律，采用指數函數擬合Swn與J的關系，即

其中

在進行Swn與J函數擬合時，J函數用的是實驗條件下求得的巖心毛管壓力Pc，而油藏條件下毛管壓力的求取，則根據油藏的毛管壓力與油水的重力差平衡原理［7］進行計算，即

式（1）～（4）中：J為J函數，無因次量；Pc為毛管壓力，Pa；σ為界面張力，mN/m；θ為潤濕角，（°）；a和b為常數，無因次量；K為滲透率，mD；φ為孔隙度，%；Swn為巖心標準化含水飽和度，%；Sw為原始含水飽和度，%；Swi為束縛水飽和度，%；ρw和ρo分別為油、水的密度，g/cm3；H為自由水面以上的油柱高度，m。

2 含水飽和度的計算

2.1 建立實驗條件下J函數與Swn的關系式

為了能夠更好地表征碳酸鹽巖儲層的非均質性，國內外許多學者在進行J函數與Swn擬合之前，一般會利用各巖樣值的大小對J函數曲線進行分類，然后再根據分類分別進行擬合［1-2，11-12］。由于研究區沒有做過有關確定界面張力及潤濕角的實驗，因此采用美國巖心公司提供的不同潤濕系統的潤濕性數據［7］（表1）。

表1 不同系統的界面張力與潤濕角數據Table1 Interfacial tension and wetting angle data in different systems

因毛管壓力實驗為空氣-水銀系統，取σ cosθ= 367 mN/m時，式（1）［用國際單位制（SI制）的實用單位］可表示為

應用公式（5）將研究區利用壓汞法獲得的27條毛管壓力曲線轉換為J函數曲線（圖1）。根據表2中27塊巖樣的值大小，將研究區的J函數分為2類，即值大于5的為Ⅰ類［圖2（a）］，值小于5的為Ⅱ類［圖2（b）］。

根據上述分類，選取各類曲線平直段的數據分別進行擬合。其中：擬合Ⅰ類曲線選取圖2（a）中飽和度為10%～90%的數據點；擬合Ⅱ類曲線利用圖2（b）中飽和度為15%～90%的數據點。通過擬合Swn與J的關系，得到Swn的計算公式，即

Ⅰ類為

Ⅱ類為

圖1 中東某油田J函數曲線Fig.1 J-function curves of an oilfield in the Middle East

表2 中東某油田27塊巖樣的巖心數據及其分類標準Table2 Core data and the classification criteria of 27 rock samples from an oilfield in the Middle East

續表2

圖2 不同類別的J函數曲線Fig.2 J-function curves of different classes

2.2 計算油藏條件下的J函數

在確定Swn與J函數的關系式后，只要求得油藏條件下的J函數值，就可以得到目的層內任意一點的Swn。根據公式（1），只有確定了油藏條件下的σ，θ，K，φ和Pc值，才能進一步確定J函數值。

（1）確定σ和θ值。由于研究區油藏條件為水-油接觸系統，因此取θ=30°，σ=30 mN/m，σcosθ= 26 mN/m（參見表1）。

（2）確定K和φ值。研究區是利用巖心標定孔隙度測井（密度、中子及聲波）求得φ，因聲波曲線與巖心分析孔隙度相關性最好，所以利用聲波曲線計算φ；K采用巖心分析孔隙度與巖心分析滲透率擬合的關系式確定。

（3）確定Pc值。按公式（4）計算油藏條件下Pc值。其中：油和水的密度由PVT分析資料確定；自由水界面（FWL）由RFT等地層壓力資料分析確定；油柱高度H為FWL的深度與目的層的深度差值。

確定好上述參數后，即可確定油藏條件下的J函數值，再根據計算出的K和φ，按照所屬類別選用不同的含水飽和度公式確定油藏的Swn。由J函數確定的含水飽和度（目的層的）是Swn，其值為可動水飽和度，不含束縛水飽和度。因此，需要通過公式（3）的換算才能得到目的層的總含水飽和度，其中Swi由研究區相滲實驗數據擬合公式計算得到［2］。擬合計算公式如下：

3 應用效果分析

3.1 與相滲分析結果對比

從理論上來講，試油結果中出油層段的Swi與Sw應相等。筆者選取了研究區3口井的試油結果，均為純油層，并且屬于同一測井小層（Sar-3）的層段。利用J函數計算得到的Sw與Sar-3小層相滲實驗數據擬合公式計算得到的Swi進行對比（圖3），結果顯示：A-2井在2 765.5～2 767.5 m和2 771.5～2 772.0 m等非儲層段、A-4井在2 659～2 660 m非儲層段以及A-5井在2 872～2 874 m非儲層段的Sw與Swi吻合較差，且Sw明顯偏大。經分析認為，造成2種方法計算非儲層段含水飽和度值差異較大的原因是：①由于表2中擬合J函數與Swn關系的數據點均為儲層段（孔隙度值均大于研究區孔隙度值的下限6%）的樣本點，擬合的公式可能不適用于非儲層對Swn的計算，因此造成Swn的結果不準確；②由于相滲實驗分析數據多為儲層段的點，擬合計算Swi的公式可能不適用于非儲層段，因此依據公式（8）建立的關系可知K值越小，Swi越大，通過其擬合計算的非儲層段的Swi可能偏大，最終將造成Sw偏大［參見公式（3）］。

圖3 A-2，A-4和A-5井J函數與相滲法計算Sw值隨深度變化對比圖Fig.3 Comparison of water saturation calculated by J-function method and relative permeability among A-2 well,A-2 well and A-5 well

3.2 與試油結果對比

圖4為A-4井測井解釋圖，其最后一道為J函數計算的飽和度值。該井通過RFT資料確定的自由水界面（FWL）為2 810 m。在深度為2 663～2 780 m的2個試油層段其試油結果均顯示產油、不產水，該層段計算的含水飽和度值均低于50%；在深度為2 780～2 810 m的試油層段，隨深度的加深電阻率值逐漸降低，且現場生產顯示該段油水同產，其解釋的含水飽和度由50%逐漸增大到100%。綜上可知，解釋的含水飽和度值與實際流體分布情況較吻合。

圖4 A-4井測井解釋結果Fig.4 Log interpretation results of A-4 well

4 結論

（1）J函數法能夠較為準確地確定非均質碳酸鹽巖儲層的原始含水飽和度。

（2）利用J函數確定原始含水飽和度的過程中，應用了測井解釋的孔隙度、滲透率以及RFT確定的FWL等數據，這些數據的準確與否為Sw的計算精度奠定了基礎，但同時也使得J函數確定Sw具有一定的局限性。

［1］Obeida T A，Al-Mehairi Y S，Suryanarayana K.Calculation of fluid saturations from log-derived J-functions in giant complex middle-East carbonate reservoir［C］.SPE Annual Conference and Exhibition，Dallas，2005.

［2］胡勇，于興河，陳恭洋，等.平均毛管壓力函數分類及其在流體飽和度計算中的應用［J］.石油勘探與開發，2012，39（6）：733-738.

［3］陳康，趙新軍.確定儲層含油飽和度的測井——毛管資料綜合法［J］.江漢石油學院學報，1995，17（2）：51-55.

［4］管耀.利用毛管壓力資料求原始含油飽和度方法探討——以綏中36-2油田23井為例［J］.勘探地球物理進展，2009，32（5）：365-369.

［5］何更生.油層物理［M］.北京：石油工業出版社，1994：210-213.

［6］Elgaghah S A，Tiab D，Osisanya S O.A new approach for obtaining J-function in clean and shaly reservoir using in situ measurements［J］.Journal of Canadian Petroleum Technology，2001，40（7）：30-37.

［7］歐陽健.石油測井解釋與儲層描述［M］.北京：石油工業出版社，1994：221-229.

［8］廖敬，彭彩珍，呂文均，等.毛管壓力曲線平均化及J函數處理［J］.特種油氣藏，2008，15（6）：73-75.

［9］王培春，呂宏志.低阻油層含水飽和度計算方法研究——以渤海A油田新近系低阻油層為例［J］.中國海上油氣，2010，22（2）：104-107.

［10］羅厚義，湯達禎，湯永梅.利用毛管壓力預測碳酸鹽巖油藏油水界面［J］.油氣地質與采收率，2013，20（2）：71-73.

［11］王慶勇.利用毛細管壓力計算低滲儲層可動水飽和度［J］.斷塊油氣田，2013，20（1）：85-88.

［12］鹿克峰，徐振中，馮景林.一種定量表征油水過渡區飽和度分布的實用方法［J］.中國海上油氣，2011，23（6）：387-390.

（本文編輯：楊琦）

Using J-function method to calculate saturation of carbonate reservoirs

SIMA Liqiang1，LI Qing1，YANG Yi1，2，CHEN Qiang1，3
（1.School of Geoscience and Technology，Southwest Petroleum University，Chengdu 610500，China；2.Zhanjiang Branch，CNOOC，Zhanjiang 524057，Guangdong，China；3.PetroChina Tarim Oilfield Company，Kolar 841000，Xinjiang，China）

Calculating the initial water saturation of heterogeneous carbonate reservoirs is always a challenge to log interpretation,especially under the condition of lack of data such as sealed core which can directly determine the water saturation.Taking an oilfield in the Middle East as an example,according to the data of capillary pressure curve, physical properties analysis and RFT test,this paper discussed J-function method that uses capillary pressure curve to calculate initial water saturation of heterogeneous carbonate reservoirs.The water saturation calculated by this method agrees quite well with that interpreted by relative permeability experiment data and oil test results.This study provides a reasonable and reliable method to ensure water saturation of heterogeneous carbonate with no water saturation of cores. Key words：carbonate reservoir；capillarypressure curve；J-function method；initial water saturation

TE122.23

A

1673-8926（2014）06-0106-05

2014-05-09；

2014-07-20

國家重大科技專項“孔隙型碳酸鹽巖儲層測井評價技術研究”（編號：2011ZX05031-003-006HZ）資助

司馬立強（（1961-），男，教授，博士生導師，主要從事油氣測井方法、解釋及地質應用方面的科研與教學工作。地址：（610500）四川省成都市新都區西南石油大學地球科學與技術學院。E-mail：smlq2000@126.com。