劉 犟，鄭榮才，茍幸福，王昌勇

（1.油氣藏地質(zhì)及開發(fā)工程國家重點實驗室（成都理工大學），成都610059；2.中國石油天然氣集團公司 長慶油田超低滲透油藏第四項目部，甘肅 慶陽745100）

研究區(qū)位于鄂爾多斯盆地西南部，靠近吳旗和白豹生油洼陷中心，區(qū)域構造主要以天環(huán)拗陷為主［1］（圖1）。雖然湖泊水體性質(zhì)的研究屬于基礎地質(zhì)范疇，但因長8油層組埋藏深度相對較大，導致研究難度加大，相關研究較少。本文利用地球化學公式計算鎮(zhèn)原地區(qū)長8油層組沉積時的古鹽度，并且結合沉積相和古地理特征，恢復研究區(qū)的古岸線、古水道和古鹽度，闡釋古鹽度的分帶性，重現(xiàn)湖泊鹽度特征和變化規(guī)律，對長8油層組儲層特征研究及油藏勘探開發(fā)具有重要意義。

1 古鹽度計算

在計算古鹽度時，通常有3種方法：用于研究碎屑巖古鹽度環(huán)境的微量元素法、適用區(qū)別淡水環(huán)境與咸水環(huán)境的同位素法、用于辨別半咸水環(huán)境的沉積磷酸鹽法。其中微量元素法［2］因硼元素容易確定，且在反映鹽度方面敏感性強，經(jīng)常被當作確定鹽度的指標元素；同位素法面向的是區(qū)別淡水環(huán)境與咸水環(huán)境，因而更多適用于研究海相碳酸鹽巖［3，4］的鹽度；沉積磷酸鹽法容易受黏土或成巖作用影響［5］，因而在確定鹽度時精確度較低。因此，為更準確地計算古鹽度，本文采用微量元素法。

1.1 計算方法

在地球化學分析研究中，因硼的含量容易被確定，并且對沉積環(huán)境的古鹽度敏感，因而在選取反映古鹽度的指標［6］時，硼元素通常不可少。在自然水體中，硼含量與水體的總體鹽度呈線性函數(shù)關系。但因自然水體中含有吸附硼的黏土，所以水體鹽度最終與硼含量呈雙對數(shù)關系，用如下方程［7］表示

式中：wB表示硼被吸附的質(zhì)量分數(shù)（wB/10－6）；S表示鹽度（鹽度即為NaCl的質(zhì)量分數(shù)：wNaCl/‰）；C1、C2為常數(shù)。在該方程表示的古鹽度理論基礎上，還可延伸出Adamas和Couch兩種計算古鹽度的公式。

圖1 鎮(zhèn)原地區(qū)構造位置及長8油層組沉積相圖Fig.1 The sedimentary facies of Member 8oil layer of Yanchang Formation in the Zhengyuan area and the location of study area

（1）Adamas公式［8］為

其中：S指古鹽度；wB指“相當硼”的現(xiàn)實質(zhì)量分數(shù)（wB/%），在實際計算中，通常要把其單位轉換為wB/10－6。需注意的是，不同區(qū)域的“相當硼”計算公式不同。如在某泥巖樣品中，伊利石含量居多，其所含有的“相當硼”則需用另一種公式，即Walker校正公式［9］來計算

其中：8.5作為常數(shù)存在，是指純伊利石內(nèi)所含有的K2O理論濃度；而wB－測和wK2O－測分別指實際測試的B含量、K2O含量。

本文主要采用了17個泥巖樣本做微量元素分析和X射線衍射分析。分析結果顯示，在鎮(zhèn)原地區(qū)長8油層組區(qū)域的泥巖中，黏土礦物含量占據(jù)最大比例，平均質(zhì)量分數(shù)達到47.76%；而石英、鉀長石等含量相對較低（表1）。此外，微量元素分析結果表明，鉀長石內(nèi)部K的質(zhì)量分數(shù)通常在3%左右，不會對古鹽度計算產(chǎn)生明顯影響，因而可不納入計算；但鉀長石內(nèi)的K元素對黏土礦物的影響卻較大，可直接作為黏土礦物中的K元素參與計算。同時，利用微量元素分析法對B元素與K元素的含量進行測算，結果表明，這2種元素的含量均可作為古鹽度定量計算的依據(jù)。所以，計算鎮(zhèn)原地區(qū)長8油層組整體的古鹽度時最終可得出如表1所示的古鹽度范圍，其平均值達到41.57‰。

（2）Couch公式［9］為

式中：wB＊指 “校正硼”的現(xiàn)實含量（質(zhì)量分數(shù)），在實際計算中，通常要把單位轉換為wB＊/10－6。wB＊可以由Couch校正公式表示為

其中：wi、ws＊、wk分別指伊利石、蒙脫石、高嶺石的實際測試的質(zhì)量分數(shù)，各系數(shù)與硼元素的吸收強度呈正比關系。一般情況下，含復雜成分的泥巖多用該公式計算硼的含量；但需注意的是，該計算結果容易受成巖作用影響。

本文選用17個泥巖樣本，并用X射線衍射定量分析結果對B元素和相關的黏土礦物含量進行公式計算。經(jīng)過分析與計算，可確定該區(qū)域的古鹽度為7.39‰～27.06‰，平均值可達到11.89‰（表2）。

表1 鎮(zhèn)原地區(qū)長8油層組泥巖礦物成分和硼、鉀、相當硼含量及Adams公式古鹽度計算Table 1 X-ray diffraction data of mudstone，the data of B，K and the calculated data of equivalent boron and paleosalinity by Adams method in Member 8oil layer of Yanchang Formation in Zhenyuan area

表2 鎮(zhèn)原地區(qū)長8油層組B、黏土礦物、“相當硼”含量和Couch公式古鹽度計算Table 2 The data of B and clay minerals and the calculated“equivalent boron”content and paleosalinity by Couch method in Member 8oil layer of Yanchang Formation in Zhenyuan area

在計算古鹽度時，通常會將wSr/wBa與wB/wGa作為判別指標。而在實際的沉積物數(shù)據(jù)分析中，這2個比值也與古鹽度呈現(xiàn)可靠的正比關系［5，10］。鎮(zhèn)原地區(qū)的wSr/wBa為0.19～0.34，wB/wGa為2.33～27.28（表3）。利用 Adams公式和Couch公式對鎮(zhèn)原地區(qū)長8油層組做古鹽度定量分析，計算結果表明，古鹽度、wSr/wBa和wB/wGa的趨勢保持相似的一致（圖2）。這也是說，在研究區(qū)長8油層組沉積期湖泊水體鹽度測算方面，Adamas公式與Couch公式都具有很高的可信度。

圖2 不同古鹽度恢復方法結果對比Fig.2 Comparison of the results of paleosalinity recovery by different methods

1.2 計算結果討論

通過對比表1和表2可以發(fā)現(xiàn)，采用Adamas公式與Couch公式定量計算得出的鹽度范圍與平均值都存在明顯的差異。這主要是因為二者所適用的條件不同；加之黏土礦物成巖后生作用影響，導致差異進一步加劇。

鎮(zhèn)原地區(qū)長8油層組泥巖的黏土礦物中含有較多的伊/蒙混層，其質(zhì)量分數(shù)平均值達到49%（表2）。而在成巖演化中，伊/蒙混層中的蒙皂石會逐步變成伊利石，從而計算出的“校正硼B(yǎng)＊”含量低于原始沉積物中的換算值，使得到的古鹽度值比原始沉積物中的真實值要小。研究區(qū)長8油層組泥巖黏土礦物內(nèi)的伊利石校正含量（質(zhì)量分數(shù)）平均達到了62%，屬于以伊利石為主的黏土巖，更適合采用Adams公式進行古鹽度的計算。同時，由于鎮(zhèn)原地區(qū)在長8油層組沉積時期為淺水三角洲沉積體系，河流涌入湖泊促使湖水鹽度被稀釋，使鹽度具有帶狀分布特征。利用Adamas公式定量計算出的水體古鹽度橫跨了半咸水鹽度帶和咸水鹽度帶，相比Couch公式可靠性更高。

表3 鎮(zhèn)原地區(qū)長8油層組微量元素特征Table 3 The data of trace elements and calculated paleosalinity data by Couch method in Member 8oil layer of Yanchang Formation in Zhenyuan area

2 古鹽度平面分布規(guī)律

利用Adamas公式和Couch公式對鎮(zhèn)原地區(qū)長8油層組的古鹽度進行計算，其最終的平面分布與變化規(guī)律呈現(xiàn)明顯的相似性。根據(jù)對鎮(zhèn)原長8油層組古湖泊水體的古鹽度定量計算發(fā)現(xiàn)：研究區(qū)北部殷家城－合道一帶、南部鎮(zhèn)原－太平一帶及東北部太白梁－賀旗一帶水體鹽度較高，局部如鎮(zhèn)原－太平一帶鹽度明顯大于其周邊地區(qū)；郭塬－方山鄉(xiāng)以及研究區(qū)東南一隅鹽度值較低，其中郭塬－方山鄉(xiāng)最低（圖3）。湖泊水體鹽度的這一分布特征與當時河流主要入注方向及位置有密切關系：郭塬－方山鄉(xiāng)一帶為主河道活動區(qū)，湖泊水體鹽度稀釋程度最高，因此水體鹽度值最低，應該可以代表河流作用較強的方向，即主物源方向；而其余地區(qū)隨著河流作用的減弱，湖泊水體鹽度逐漸升高。

圖3 鎮(zhèn)原地區(qū)長8油層組古鹽度平面分布與變化圖Fig.3 The paleosalinity zone of Member 8oil layer of Yanchang Formation in Zhengyuan area

鎮(zhèn)原地區(qū)長8油層組古鹽度恢復后，其結果顯示，高值鹽度區(qū)域與此前觀察巖心時所描繪的湖灣區(qū)分布位置幾乎一致。在此次研究中，選取巖心的鉆井共32口，只在18口井的巖心觀察中找到錯落有致的浪成交錯層理（圖4）。浪成交錯層理的出現(xiàn)，主要是河口壩或者遠砂壩受到河流長期波浪作用影響的水下沉積環(huán)境中形成的。而且，它的浪控帶位置和Adamas公式所恢復的≥40‰的古鹽度區(qū)域吻合。這也是說，當河流作用快速減弱時，半咸水—咸水的可能性會在湖灣區(qū)明顯加大。

圖4 鎮(zhèn)原地區(qū)長8油層組浪控帶與等鹽度關系Fig.4 The relationship between paleosalinity zone and wave active zone of Member 8oil layer of Yanchang Formation in Zhengyuan area

3 研究意義

利用Adamas公式和Couch公式恢復鎮(zhèn)原地區(qū)長8油層組所具有的古鹽度，并且分析其平面分布和變化規(guī)律，可得出以下與巖相古地理和成巖期的流體性質(zhì)相關的幾點認識：①湖岸線位置的確定必須要將古鹽度定量分析作為計算的前提。②不論是巖相古地理還是成巖期，在恢復古鹽度時都必須先明確主次物源方向，精確掌握湖盆中心。③因鎮(zhèn)原地區(qū)長8油層組所在的湖盆古鹽度呈現(xiàn)明顯的半咸水—咸水特性，所以成巖期的湖泊水或其他水體都具備鹽水特質(zhì)。而這種鹽水特質(zhì)，使得砂巖早成巖階段可以獲得較為豐富的富鈉孔隙水，輔助較早時期的綠泥石環(huán)邊及濁沸石能更好地實現(xiàn)膠結作用，從而保證原生粒間孔隙能夠長期被有效保存［11］。同時，在中成巖時期，濁沸石被有機酸溶蝕后會產(chǎn)生較多的、能提高儲層質(zhì)量的次生孔隙［12，13］。所以，在尋找優(yōu)質(zhì)儲層時，高鹽度沉積環(huán)境通常是首選區(qū)域。

4 結論

a.鎮(zhèn)原地區(qū)長8油層組古鹽度恢復結果顯示，該區(qū)域湖泊水質(zhì)呈現(xiàn)半咸水—咸水的特點，Adams公式計算結果比Couch公式計算結果更為可靠。

b.不同方位湖泊水體鹽度不同，平面分布有所差異，變化規(guī)律也不相同，鹽度主要由北東、北西向中央與東南方向增大，與河流注入方向保持一致。

［1］楊俊杰.鄂爾多斯盆地構造演化與油氣分布規(guī)律［M］.北京：石油工業(yè)出版社，2002.Yang J J.Basin Tectonic Evolution and Oil and Gas Distribution in Erdos［M］.Beijing：Petroleum Industry Press，2002.（In Chinese）

［2］鄭榮才，柳梅青.鄂爾多斯盆地長6油層組古鹽度研究［J］.石油與天然氣地質(zhì)，1999，20（1）：20－25.Zheng R C，Liu M Q.Study on palaeosalinity of Chang-6oil reservoir set in Ordos Basin［J］.Oil ＆Gas Geology，1999，20（1）：20－25.（In Chinese）

［3］劉寶珺，曾允孚.巖相古地理基礎和工作方法［M］.北京：地質(zhì)出版社，1985：321－323.Liu B J，Zeng Y F.Geographical Basis and Methods of Lithofacies［M］.Beijing：Geological Publishing House，1985：321－323.（In Chinese）

［4］Matthew C P，Eduard G R，Henry P S，etal.Highresolution paleosalinity reconstruction from Laguna dela Leche，north coastal Cuba，using Sr，O，and C isotopes［J］.Science Direct，2007，245：535－550.

［5］Muller G.Sedimentary phosphate method for estimating paleosalinities：Limited applicability［J］.Sci-ence，1969，163：812.

［6］周仰康，何錦文，王子玉.硼作為古鹽度指標的應用［C］//沉積學和有機地球化學學術會議論文選集.北京：科學出版社，1984：55－57.Zhou Y K，He J W，Wang Z Y.Application of B as an index for palaeosalinity［C］//Collected Papers of Sedimentology and Organic Geochemistry Academic Conference.Beijng：Science Press，1984：55－57.（In Chinese）

［7］Walker C T，Price N B.Departure curves for computing paleosalinity from boron in illites and shales［J］.AAPG Bulletin，1963，47（5）：833－841.

［8］Adams T D，Haynes J R，Walker C T.Boron in Holocene illites of the dovey estuary，wales，and its relationship to paleosalinity in cyclothems［J］.Sedimentology，1965，4：189－195.

［9］Walker C T.Evaluation of boron as a paleosalinity indicator and its application to offshore prospects［J］.AAPG Bulletin，1968，52：751－766.

［10］殷科華，葉德書，沈大興，等.息烽—遵義早石炭世大塘期巖相古地理特征［J］.沉積學報，2009，27（4）：606－613.Yin K H，Ye D S，Shen D X，etal.Characteristics of lithofacies paleogeography during Datangian Age（Early Carboniferous）in Xifeng and Zunyi，Guizhou［J］.Acta Sedimentologica Sinica，2009，27（4）：606－613.（In Chinese）

［11］王昌勇，王成玉，梁曉偉，等.鄂爾多斯盆地姬塬地區(qū)上三疊統(tǒng)延長組長8油層組成巖相［J］.石油學報，2011，32（4）：597－604.Wang C Y，Wang C Y，Liang X W，etal.Diagenetic facies of the Chang 8oil-bearing layer of the Upper Triassic Yanchang Formation in the Jiyuan area，Ordos Basin［J］.Acta Petrolei Sinica，2011，32（4）：597－604.（In Chinese）

［12］楊曉萍，裘懌楠.鄂爾多斯盆地上三疊統(tǒng)延長組濁沸石的形成機理、分布規(guī)律與油氣關系［J］.沉積學報，2002，20（4）：628－632.Yang X P，Qiu Y N.Formation process and distribution of lau-montite cements in Yanchang Formation（Upper Triassic）of Ordos Basin［J］.Acta Sedimentologica Sinica，2002，20（4）：628－632.（In Chinese）

［13］梁浩，羅權生，孔宏偉，等.三塘湖盆地火山巖中沸石的成因及其儲層意義［J］.沉積學報，2011，29（3）：537－543.Liang H，Luo Q S，Kong H W，etal.Formation and distribution of zeolite in volcanic rock and its effect on reserviors in Santanhu Basin［J］.Acta Sedimentologica Sinica，2011，29（3）：537－543.（In Chinese）