劉凱+郭高軒+謝振華+李志萍+劉穎超+寇文杰

摘要：北京市平原區(qū)地下水系統(tǒng)劃分為3個相對獨立的區(qū)域地下水系統(tǒng)，選取其中典型的永定河及潮白河水系進行分析，通過17個取樣點不同深度的83組樣品的礦物成份、粘土礦物成份、常量元素及微量元素分析其沉積環(huán)境及沉積特征。通過聚類分析、礦物組合和風化程度計算，確定了兩大水系沉積特征既有相同點，又有差異性。首先用聚類分析半定量分析出兩大系統(tǒng)的相同點和差異性，通過風化程度計算出潮白河系統(tǒng)松散沉積物風化程度大于永定河系統(tǒng)。其次通過礦物成分衍射實驗證明兩條剖面主要沉積物礦物成分及其中粘土礦物類型均相同，兩條剖面均處于河流沉積環(huán)境且水分條件優(yōu)越，其中永定河剖面方解石及白云石百分含量和高于潮白河，潮白河剖面角閃石百分含量高于永定河。再次通過巖石化學成分分析兩條剖面常量元素、微量元素主要類型相同，根據(jù)不同元素組合定量分析永定河及潮白河淺層松散沉積物差異性沉積特征。最后使用半定量、定量方法綜合分析了Y6、Y14取樣點沉積環(huán)境的典型性特征。這些方法的應用為河流沉積時序、古河道變遷提供了最直接的量化證明，為地下水系統(tǒng)劃分、地下水系統(tǒng)邊界確定提供了最有利最直接的證據(jù)，也將為未來水化學場研究、水資源合理開發(fā)利用等提出新的思路和方法。

關鍵詞：永定河；潮白河；沉積特征； 聚類分析；礦物成分

中圖分類號：P641文獻標志碼：A文章編號：

1672-1683（2015）001-0188-04

The characteristics analysis of shallow sediments

LIU Kai，GUO Gao-xuan，XIE Zhen-hua，LI Zhi-ping，LIU Ying-chao，KOU Wen-jie

（Beijing Institute of Hydrogeology and Engineering Geology，Beijing 100195，China）

Abstract：The beijing plain groundwater system is divided into three relatively independent groundwater subsystem，this study analyzes the sedimentary environment and sedimentary characteristics of the typical profile in the Yongding river and the Chaobai river by the mineral composition，clay mineral composition，major elements and the trace elements of the 83 different depth samples in 17 sampling points.In this paper，the similarities and difference of two systems was determined by clustering analysis，typical mineral assemblage method and weathering degree method.Firstly semi-quantitative analysis the similarities and difference of the typical systems，the weathering degree of Chaobai river was calculated much more bigger than the Yongding river.Then it was found by XRD that the the mineral composition and the clay mineral component are same in the two systems，also they are blong to the fluvial deposit environment and they have a sufficient hydrodynamic condition.The sum of calcspar and dolomite in Yongding river is more than Chaobai river，but the hornblende is less than the Chaobai river.The two systems have the same major elements and trace elements，then quantitative analyzes the sedimentary characteristics differences between the systems by the different element associations.Finally determines the typical sedimentary characteristics in sampling points Y6 and Y14.The result of study is significant for the study of evolution，vicissitudes and identification of ancient river source，and the application of these methods provide the most direct quantitative proof on division and boundary confirmationof groundwater systems，at the same time this article puts forward a new method of hydrogeochemical environment analysis and the reasonable exploitation of water resource.endprint

Key words：Yongding river；Chaobai river；sedimentary characteristics；cluster analysis；the mineral composition

北京地處華北平原北部，地形西北高、東南低，西部和北部山脈屬太行山脈和燕山山脈，市區(qū)位于向東南傾斜的沖洪積平原上。根據(jù)北京市平原區(qū)水系及其河流沖洪積扇的分布特征，整體上將北京市平原區(qū)地下水系統(tǒng)劃分為3個相對獨立的區(qū)域地下水系統(tǒng)[1]（見圖1），本文選取永定河水系和潮白河水系兩大典型水系分析淺層沉積特征分析，闡述兩大水系沉積環(huán)境，為地下水化學特征及地下水系統(tǒng)邊界驗證、古河道變遷、水資源合理開發(fā)利用等提供了直接的科學依據(jù)。

永定河及潮白河兩條典型剖面上由于河流流經(jīng)母巖區(qū)性

圖1北京市平原區(qū)地下水系統(tǒng)及取樣點分布

質(zhì)不同，其沉積環(huán)境各有差異，其各自沉積環(huán)境造就典型剖面的風化及組合特征、沉積物的搬運特征等各有差異[2-4]。永定河主要流經(jīng)西山，西山屬太行山脈，地表組成物質(zhì)以石灰?guī)r類為主，其次是火山碎屑巖類；石灰?guī)r風化主要是化學溶解作用，風化后土層較薄，稍黏、鈣質(zhì)豐富，抗酸力強。潮白河主要流經(jīng)北山燕山山脈，地表組成物質(zhì)以花崗巖、片麻巖類為主，其次為石灰?guī)r和砂礫巖類。

本文物源分析樣品取自北京市平原區(qū)17個取樣點不同層位（見圖1），選取其中83個典型樣品分別進行巖礦實驗、礦物實驗和黏土礦物實驗。

礦物之間具有嚴格的共生關系，所以礦物組合是沉積環(huán)境變化極為敏感的指示劑，由于不同地區(qū)沉積特點不同而其具體組合形式、分析方法不同[3]。同一沉積平原中，同時期沉積物的碎屑組分一致，不同時期沉積物所含的碎屑物質(zhì)不同，利用不同時期水平方向上礦物種類和含量變化，可分析其沉積特征。本文在聚類分析[5-6]、相關分析等方法的基礎上確定研究區(qū)的礦物組合形式及特征。

1聚類分析

聚類分析是根據(jù)事物本身的特征對被研究對象進行分類的數(shù)理統(tǒng)計方法，是進行元素相關分析和地球化學分區(qū)的一種有效手段。其基本思想是先定義樣品間距離及類與類之間的距離；開始時各樣品看成一類，將距離最近的兩類合并；重新計算新類與其它類的距離，再將距離最近的兩類合并；再計算新類與其它類的距離……這樣一步步的計算，每一步減少一類，直至所有樣品都合并成一類為止。

在計算過程中常分為三個步驟：數(shù)據(jù)標準化，計算分類統(tǒng)計量，建立分類系統(tǒng)；其結果以分類譜系圖形表示。本文數(shù)據(jù)標準化采用總和標準化，即分別求出各聚類要素所對應的數(shù)據(jù)總和，以各要素數(shù)據(jù)除以該要素數(shù)據(jù)總和。距離是事物之間差異性的測度，差異性越大，相似性越小，故距離是系統(tǒng)聚類分析的依據(jù)和基礎；系統(tǒng)聚類方法里常用的距離計算有歐氏距離、切比雪夫距離、馬氏距離、蘭氏距離等。常用的聚類方法有最短距離聚類法、直接聚類法、最長距離聚類法、中間距離法、離差平方和聚類法等。

本文結合樣品實際情況，以8種常規(guī)元素（SiO2、Al2O3、K2O、MgO、CaO、Na2O、Fe2O3及FeO）進行聚類[7]，聚類方法采用最短距離法，距離采用歐氏距離平方法，聚類的距離比較短，聚類圖譜分類明顯（見圖2）。

圖2歐式距離/最短距離聚類圖

在圖1地下水系統(tǒng)劃分基礎上結合圖2可以分析所取樣品的半定量關系，對所取樣品進行歸類。當聚類距離為235時，所有樣品同屬于北京平原區(qū)松散沉積物沉積環(huán)境中。當聚類距離為195時，樣品Y1-Y10（Y6除外）屬于永定河沖洪積扇沉積環(huán)境，而樣品Y11-Y17屬于潮白河沖洪積扇沉積環(huán)境，這也定性說明兩個系統(tǒng)的沉積環(huán)境各有其獨特的沉積特征。當聚類距離在179～189之間時，樣品Y6、Y7及Y14具有獨特的沉積特征，其沉積環(huán)境既不完全屬于潮白河系統(tǒng)，也不完全屬于永定河系統(tǒng)，結合北京市平原區(qū)系統(tǒng)劃分，可以看出這幾個樣品均處于永定河古河道及潮白河兩個沖洪積扇交匯沉積地帶。

2風化程度

Nesbitt 和Young（1982，1984，1989）[8]引進了化學變異指數(shù)（CIA），此指數(shù)提供了一種定量化硅酸鹽礦物風化（化學變異） 度的方法。CIA 的計算依據(jù)以下公式：

CIA=[（Al2O3）/（CaO+Na2O+K2O+Al2O3）]×100

此處氧化物是以克分子含量形式表示的，CaO只是硅酸鹽礦物中的CaO，而磷酸鹽和碳酸鹽礦物中的CaO 除外。CIA 值也可用A-CN-K三角圖來表示，三個端元分別為Al2O3，CaO +Na2O 和K2O，其含量為克分子含量[9]。

根據(jù)樣品分析結果，計算得83個永定河系統(tǒng)樣品化學變異指數(shù)值分布在46～56之間，潮白河樣品化學變異指數(shù)分布在60～70之間，兩個典型剖面化學變異指數(shù)值存在明顯差異，潮白河系統(tǒng)物質(zhì)風化程度大于永定河系統(tǒng)，也定量計算出兩條典型剖面的風化沉積程度區(qū)別。

3礦物組合

3.1礦物成分

從礦物衍射報告可以看出，17個樣品的主要沉積成分為石英、鉀長石、斜長石、方解石、白云石、角閃石及黏土礦物（見圖3）。永定河及潮白河剖面由于流經(jīng)地層不同，具有不同的沉積特征，其中永定河剖面方解石+白云石百分含量在3%～25%，石英含量主要在40%以下，角閃石含量在06%～22%；潮白河剖面方解石+白云石百分含量在1%～2%之間，石英含量在40%～45%之間，角閃石含量在1%～4%之間。

圖3礦物成分百分含量柱狀圖

大量的黏土礦物是在表生風化作用中形成的[10-11]，黏土礦物實驗結果顯示，17個取樣點83個樣品的黏土礦物類型均為伊利石+伊利石/蒙脫石混層礦物+高嶺石+綠泥石（見圖4）。趙永勝（1993） 在研究云南星云湖時發(fā)現(xiàn)河流環(huán)境（包括分流河道和泛濫平原）不利于綠泥石的保存，而是出現(xiàn)伊利石/蒙脫石-高嶺石組合。當水分條件優(yōu)越時，伊利石/蒙脫石混層數(shù)量會增多，伊利石減少；當水份條件惡劣時，伊利石/蒙脫石混層會部分乃至絕大部分轉變成伊利石。分析結果顯示試驗樣品中伊利石/蒙脫石混層、伊利石平均相對含量分別為6435%、1730%，由此也看以看出永定河及潮白河剖面均處于河流沉積環(huán)境且水分條件優(yōu)越。endprint

圖4黏土礦物相對百分含量柱狀圖

3.2巖化分析

從17個取樣點83組巖石化學成分分析結果可以看出（見圖5），兩個典型剖面的常量元素有SiO2、Al2O3、K2O、MgO、CaO、Na2O、Fe2O3及FeO，8種常量元素百分含量高達887%～951%。

圖5常量元素百分含量柱狀圖

圖6常規(guī)元素組合三角圖

圖7常規(guī)元素構造環(huán)境判別圖

根據(jù)常規(guī)元素三角圖及構造環(huán)境判別圖可以看出潮白河及永定河明顯沉積特征的差異，使用化學變異指數(shù)值和礦物組合能更加明顯展示兩條剖面的沉積環(huán)境區(qū)別及其各自的沉積特征。此外，還能展示出在交互沉積環(huán)境中Y6、Y7、Y14的特殊沉積特征。

3.3微量元素

永定河及潮白河主要微量元素有P、Mn、F、Sr 、Cr、Cu、Zn、Pb、B、Ni、Co、As 、N、Sn、Sb、Mo、Se等，17種微量元素的百分含量僅為015%～037%。

由于不同巖層中Ca、Sr、Mg含量不同，礦物的抗風化能力不同。地下水在與圍巖接觸的過程中，巖石中的鍶元素以Sr2+ 形式進行離子交換進入地下水，成為地下水的組分之一。由于Sr、Sr/Ca、Sr/Mg具有較好相關性（見圖8），故本文選取微量元素Sr/Ca及CIA值（見圖9）進行綜合分析。

圖8Sr與Sr/Ca、Sr/Mg相關關系

從圖中可以看出，永定河中Sr/Ca值只要分布在0002 3～

圖9Sr/Ca構造環(huán)境判別圖

0011 4之間，而潮白河Sr/Ca值主要分布在0008～0027，從微量元素也可以看出兩個剖面沉積環(huán)境的不同。

4綜合應用

根據(jù)聚類分析，可以半定量分析樣品Y6、Y7、Y14的沉積是多期河流交疊、縱橫沉積形成的。對Y6、Y7、Y14取樣點不同層位分別取5個、7個、6個樣品，根據(jù)不同樣品的化學變異度、常規(guī)元素含量、礦物組合特征及構造環(huán)境判別圖進行綜合判斷。

取樣點Y6和Y7均處于永定河地下水系統(tǒng)中（圖1），但樣品Y6-5、Y6-7、Y6-8、Y7-2沉積特征均顯示出潮白河剖面特征； Y14分布在潮白河區(qū)域中，但Y14-3、Y14-4、Y14-5均分布在潮白河剖面的最下部，接近永定河剖面。結合取樣深度、顆分定名、地層剖面，Y6樣品18～94 m的沉積特征體現(xiàn)出潮白河系統(tǒng)特征，而Y6樣品18 m以上和94 m以下的沉積特征為永定河系統(tǒng)。Y7樣品08 m以上黏質(zhì)粉土的沉積特征顯示出潮白河剖面特征，而Y7取樣點其它樣品的沉積特征屬于永定河系統(tǒng)。Y14樣品0～40 m沉積特征體現(xiàn)出潮白河系統(tǒng)特征；40～116 m沉積特征為潮白河區(qū)域，但是并非完全屬于潮白河區(qū)域，可能是潮白河與永定河共同作用形成的；而Y14樣品116 m以下的沉積特征屬于永定河系統(tǒng)，詳見圖10。

圖10取樣點沉積環(huán)境判別圖

Y6、Y14附近的沉積特征是潮白河及永定河古河道變遷及河道變遷時序的最直接證明?？梢源_定在Y6、Y14附近先進行永定河古河道沉積，而后潮白河進行沉積，最后Y6、Y14分別接受永定河、潮白河沉積，其中Y14取樣點附近又接受多期交錯沉積。

Y6、Y14附近的沉積特征也為永定河地下子系統(tǒng)及潮白河地下水子系統(tǒng)邊界確定提供了最有利最直接的證據(jù)。

5結論

（1）通過聚類分析、風化程度計算、礦物成分分析、巖石化學常量元素、微量元素等分析，證明北京市永定河及潮白河淺層松散沉積物沉積特征既具有相同點，又具有差異性。

（2）兩條剖面主要沉積物礦物成分均為石英、鉀長石、斜長石、方解石、白云石、角閃石及粘土礦物，其中粘土礦物類型均為伊利石+伊利石/蒙脫石混層礦物+高嶺石+綠泥石，證明永定河及潮白河兩條剖面均處于河流沉積環(huán)境且水分條件優(yōu)越。此外，永定河剖面方解石及白云石百分含量和高于潮白河，潮白河剖面角閃石百分含量高于永定河。

（3）兩個典型剖面的常量元素均為SiO2、Al2O3、K2O、MgO、CaO、Na2O、Fe2O3及FeO，8種常量元素百分含量高達887%～951%；主要微量元素有P、Mn、F、Sr、Cr、Cu、Zn、Pb、B、Ni、Co、As、N、Sn、Sb、Mo、Se等，17種微量元素的百分含量僅為015%～037%。通過相關分析，不同常量元素及微量元素組合驗證分析出永定河及潮白河淺層松散沉積物具有其各自沉積特征。

（4）永定河系統(tǒng)樣品化學變異指數(shù)值分布在46～56之間，潮白河樣品化學變異指數(shù)分布在60～70之間，潮白河系統(tǒng)沉積物風化程度大于永定河系統(tǒng)。

（5）通過聚類分析半定量確定Y6、Y14取樣點沉積環(huán)境的典型性，通過風化程度計算、礦物成分分析、巖石化學常量元素組合分析等方法定量綜合分析了Y6、Y14取樣點的沉積特征。這些方法的應用為河流沉積時序、古河道變遷提供了最直接的量化證明，為地下水系統(tǒng)劃分、地下水系統(tǒng)邊界確定提供了最有利最直接的證據(jù)。由于沉積特征在很大程度上決定了地下水系統(tǒng)的水化學特征，這也將對未來水化學場研究提出新的思路和方法。

參考文獻：

[1]張安京，葉超，李宇，等.北京地下水第一版[M].北京：中國大地出版社，2008：36-39.

[2]徐亞軍，杜遠生，楊江海.沉積物物源分析研究進展[J].科技情報，2007，26（3）：26-31.

[3]Dickinson W R，Valloni R.Plate set tings and provenance of sands in modern oceans[J].Geology，1980，82-86.

[4]Dickinson W R.Provenance of Nort h American Phanerozoicsandstones in relation to tectonic set ting[J] .B ul l.Geol.S oc.A m.，1983，94：222-235.

[5]韓曉飛，朱長生，張玉芬，等.關于聚類分析和診斷模型在物源分析中的應用[J].廣西輕工業(yè)，2009，8（129）：106-108.

[6]王珊珊，劉凱，劉飛，等.聚類分析在土地工程能力評價中的應用[J].中國地質(zhì)災害與防治學報，2009，20（3）：112-117

[7]魏亮，郭華明，謝振華，等.北京平原區(qū)沉積物稀土元素地球化學特征及物源意義[J].地學前緣，2010，17（6）：72-80.

[8]Nesbitt H W，Young G M.Formation and diagenesis of weat hering profiles[J].J.Geol，1989，97（2）：129-147.

[9]謝振華，劉凱，李志萍，等.基于沉積物物質(zhì)來源的地下水化學特征分析—以北京市平原區(qū)為例[J].地學前緣，2010，17（6）：81-87.

[10]Galy A，F(xiàn)rance-Lanord C.Higher erosion rates in the Himalaya：geochemical constrains on riverine fluxes.Geology，2001.29：23-26.

[11]杜德文，孟憲偉，王永吉，等.沉積物物源組成的定量判識方法及其在沖繩海槽的應用[J].海洋與湖沼，1999，30（5）：532-539.endprint