基于水化學及D、18O的柳江盆地東宮河流域地下水循環(huán)特征解析

宋 洋, 遲寶明, 谷洪彪, 秦文婧, 杜利偉

(1.中國地震局工程力學研究所, 哈爾濱 150080; 2.防災科技學院, 河北 三河 065201)

?

基于水化學及D、18O的柳江盆地東宮河流域地下水循環(huán)特征解析

宋 洋1,2, 遲寶明1,2, 谷洪彪2, 秦文婧2, 杜利偉2

(1.中國地震局工程力學研究所, 哈爾濱 150080; 2.防災科技學院, 河北 三河 065201)

運用水化學成分和環(huán)境同位素作為水循環(huán)研究的示蹤劑，揭示柳江盆地東宮河流域地下水循環(huán)特征。通過現(xiàn)場調查，設計和布置水樣的采樣點，采集和測試大氣降水、地表水及地下水樣品，對比不同時期水樣的水化學特征，結合氫氧同位素組成，最終確定不同類型地下水的補給、徑流及排泄條件。研究表明：最具供水意義的地下水主要來自東宮河流域厚層的府君山組灰?guī)r裂隙含水層，巖溶水自北向南流，與柳江盆地區(qū)域巖溶水流動方向一致。接受北部山區(qū)大氣降水入滲補給，豐水期經(jīng)由第四系弱透水層在東部落村以泉水形式排泄，在盆地北緣至東部落以西河段以地下徑流的形式向東宮河排泄?；鶐r裂隙水次之，其分布具有強烈的非均質各向異性的特點，與盆地西北部山區(qū)的饅毛組泥頁巖相比，盆地東北部花崗巖山區(qū)的裂隙水，水質好，補給水頭最高，補給源穩(wěn)定、補給路徑較短，在李莊及其附近存在自流區(qū)。研究成果為當?shù)氐乃Y源開發(fā)與利用提供了科學依據(jù)。

地下水循環(huán); 東宮河流域; 水化學成分; 氫氧同位素

地下水是賦存于復雜地質構造單元中具有動態(tài)特征的水體，地下水循環(huán)是水文循環(huán)的重要組成部分[1]，查明地下水循環(huán)的補給、徑流、排泄條件及水質演化規(guī)律是地下水合理開發(fā)利用的前提。秦皇島市地處我國嚴重缺水的華北地區(qū)，人均水資源擁有量僅為全國平均水平的1/4[2]，于1987年6月投產(chǎn)運行的秦皇島市第二供水水源地——柳江水源地的運行[3]，為緩解日益突出的水資源供需矛盾提供了有力保障。柳江水源地位于柳江盆地核心區(qū)，距秦皇島市以北約15 km的燕山腳下，盆地內(nèi)薈萃了從新太古至新生代數(shù)億年中國華北地區(qū)的地殼運動、巖漿活動、沉積環(huán)境變化及生物進化等地質現(xiàn)象的精華，被公認為“天然地質博物館”。良好的地質構造環(huán)境也為地下水的賦存提供了條件，尤其是寒武—奧陶系碳酸鹽巖含水層，水量充沛、水質優(yōu)良，各種水文地質現(xiàn)象齊全。然而，在近30 a強烈的人為干預下，20世紀80年代末期柳江盆地就因地下水的超采出現(xiàn)了農(nóng)用井干涸、地面塌陷及地裂縫的現(xiàn)象[4]。截止到2000年，盆地內(nèi)水泥廠遍布，大氣污染嚴重，同時地下水水質趨于惡化。柳江盆地從地質到水文地質條件均發(fā)生了巨大變化，地下水系統(tǒng)的循環(huán)條件及循環(huán)途徑今非昔比。目前國內(nèi)關于柳江盆地地質及水文地質條件的描述及研究成果絕大部分來源于20世紀80年代長春地質學院編撰的實習指導書[5]及水源地建立初期的水文地質勘察資料，到目前為止未見有利用同位素等先進的方法論證該地區(qū)地下水循環(huán)規(guī)律的相關報道。為滿足水文地質專業(yè)教學的迫切需要，有必要對盆地內(nèi)地下水與地表水補排關系、地下水運移規(guī)律等方面進行重新論證，同時希望本次研究成果能為當?shù)氐乃Y源開發(fā)與利用提供科學依據(jù)。

地下水中的水化學成分和環(huán)境同位素組分蘊含著地下水的補給來源、運移途徑等重要信息，是當前國內(nèi)外用于研究水循環(huán)中不同水體補給來源[6-8]、確定地表水和地下水轉化關系的較為有效且先進的技術手段[9-14]。本文選取柳江盆地的二級水文地質單元東宮河流域地下水系統(tǒng)作為研究對象，運用水化學成分和氫氧穩(wěn)定同位素作為地下水循環(huán)研究的示蹤劑，揭示柳江盆地東宮河流域地下水水化學特征及地下水與地表水的補排關系，為進一步探討柳江盆地的區(qū)域地下水循環(huán)特征提供技術支持與參考。

1 研究區(qū)概況

東宮河流域位于秦皇島市撫寧縣石門寨鎮(zhèn)東北部，四面環(huán)山，為柳江向斜東翼東部低山丘陵區(qū)一南北向延伸的小型山間盆地，南北長約9 km，東西寬約6 km，面積約54 km2。海拔100～500 m，最高峰位于盆地東北部,海拔約480 m,最低點位于盆地南部東宮河出口處,海拔約105 m。盆地東北部及西、南面整體地形坡度較大,海拔多在280 m以上,地形切割深度大于200 m,北、東部坡度較緩,地形切割深度約150 m。盆地內(nèi)部地勢較為平坦,沿河流干流流向主要分布有蛤蟆泉、李莊、程莊、溫莊、黃土營以及東部落村等自然村落。

研究區(qū)地處于中緯度地帶,東經(jīng)109°37′13.3″—109°41′,北緯40°07′11.1″—40°11′11.9″,屬暖溫帶半濕潤大陸型季風氣候區(qū)。四季分明,年氣溫差較大,夏季平均氣溫高于20℃,冬季平均氣溫低于5℃,年平均氣溫25.7℃,無霜期年均174 d。年降水量多寡變化顯著,一般為400～1 000 mm,年平均降雨量為695.5 mm,70%的降水集中于7—8月份。年均蒸發(fā)量為1 646.8 mm,最大蒸發(fā)量出現(xiàn)在5月份。

東宮河為研究區(qū)內(nèi)唯一的地表水系,主要由三條支流匯合而成，分別發(fā)源于王家峪、蛤蟆泉、娃娃峪的低山丘陵區(qū)，在黃土營處匯流，沿東部落村西側向南流，于下花野村東南部匯入大石河，并隨之最終匯入石河水庫。屬季節(jié)性山區(qū)河流，干流大體流向為自北向南流。受大氣降雨影響，夏季流量較大，暴雨后洪水立刻上漲，且暴漲暴落，分別于2012年7月14日、2013年7月13日黃土營橋下實測東宮河的流量分別為0.748，0.870 m3/s。冬季則因降水減少徑流量變小，東部落村西至黃土營以北河段形成干谷，冬夏兩季徑流量差異較大。

盆地位于柳江向斜東翼，巖層整體走向近南北，傾向西，傾角約16°，在東北角略微翹起，自東向西由老到新地層展布情況及巖性描述如下：

新元古界青白口系下馬嶺組(qx)砂頁巖和景兒峪組(qj)砂巖、泥灰?guī)r，裂隙不發(fā)育，透水性差，厚度合計約110 m，在盆地東部及東北部山區(qū)廣泛分布。古生界寒武系下統(tǒng)府君山組(∈1f)豹皮狀含瀝青質白云質灰?guī)r，位于東部落村的標準剖面厚度為146 m，與下伏地層呈平行不整合接觸關系，在東部丘陵區(qū)與盆地交界地帶廣泛分布，垂向節(jié)理裂隙發(fā)育，地表可見各種溶蝕現(xiàn)象，縱張裂隙特別發(fā)育，呈高角度近乎直立，張開程度好，且填充物少而粗或者未充填[15]。在東部落村村南發(fā)育一規(guī)模較大的落水洞，洞口近圓形，直徑約3 m，深度較大，據(jù)當?shù)卮迕穹从吃撀渌磧?nèi)水量較為充足，抽水時水中伴有砂礫石等河流相物質。古生界寒武系下統(tǒng)饅頭組(∈1m)泥頁巖和毛莊組(∈1mo)頁巖，與下伏的府君山組呈平行不整合接觸關系，由于西側山前的正斷層(70°∠65°)的存在，使得盆地中部第四紀松散堆積物以下至西側半山腰處均有分布，總厚度為183 m，巖層富水性和透水性均較差。上覆的寒武系中統(tǒng)徐莊組(∈2x)粉砂巖和張夏組(∈2z)頁巖，裂隙不發(fā)育，透水性和富水性也較差。盆地最上部覆蓋有厚度不等的第四系(Q)松散堆積物，西部東宮河河谷區(qū)厚度最大，為5～10 m厚的沖、洪積物，向東逐漸減小，至東部落村厚度減至3～5 m，成分以沉積、殘積物為主。除此之外，在盆地周邊出露有大規(guī)模的多期侵入巖體，盆地南部為中生界的閃長玢巖、東北角為中元古界肉紅色混合巖化鉀長花崗巖、西北部邊緣地帶為花崗斑巖。

綜上所述，研究區(qū)為一較為完整的水文地質單元。區(qū)內(nèi)含水層劃分為：第四系孔隙潛水含水層、府君山組灰?guī)r巖溶裂隙水含水層以及中元古界花崗巖基巖裂隙水含水層。巖溶裂隙水含水層富水性最好，其隔水頂為饅毛組的泥頁巖，隔水底板為下馬嶺組和景兒峪組砂頁巖、泥灰?guī)r，盆地南部侵入巖體以及為側向阻水邊界。根據(jù)巖溶發(fā)育的規(guī)律，盆地東部及東北部山區(qū)灰?guī)r直接出露地表，水循環(huán)速度快，巖溶發(fā)育，成為巖溶裂隙含水層的側向補給邊界，相較之下，盆地西側灰?guī)r深埋區(qū)，地下水補給路徑較長，水循環(huán)速率緩慢，不滿足巖溶條件，故判定其為側向隔水邊界。盆地西側山前的正斷層，兩盤主要巖性為塑性的泥頁巖，說明其透水性、導水性均較差，雖無鉆孔資料，無法判斷其是否完全切穿饅毛組地層，但這并不影響最終結果判定，即盆地西部邊界從地表至下馬嶺組和景兒峪組砂頁巖、泥灰?guī)r，透水性、導水性均較差，共同構成一側向隔水邊界。第四系孔隙潛水含水層，厚度不均，主要沿盆地東宮河河谷區(qū)展布，透水性、導水性由西向東逐漸變差。此外，在盆地東北部出露的花崗巖，裂隙較為發(fā)育，且補給水頭較高，在李莊周邊形成自溢區(qū)。

2 研究方法

2.1 樣品采集

采集水樣包括大氣降水、地表水、第四系孔隙水、泉水、基巖裂隙水及巖溶裂隙水。大氣降水采集時間段為2013-03-12—2013-09-23，用于測試δ18O，δD值，建立當?shù)亟邓€方程[16-17]。為滿足流域內(nèi)部地下水的循環(huán)規(guī)律研究的需求，避免隨機性，其他水樣的采集分別于豐水期(2012年7月和2013年7月)、枯水期2013-10月進行，采集水樣信息詳見表1。地下水樣品直接從井采集，采樣前先抽1～3 mim的水，δ18O、δD測試樣品選用100 ml采樣瓶，常規(guī)水化學組分測試樣品選用500 ml采樣瓶進行現(xiàn)場密封，并現(xiàn)場速測pH，Eh，電導率，DO，TDS。

表1 水樣采集基本信息

2.2 樣品測試

水樣常規(guī)水化學組分和δ18O，δD由中國科學院地理科學與資源研究所測定。其中Ca2+，K+，Mg2+，Na+測定儀器為美國珀金埃爾默儀器有限公司電感耦合等離子體光譜儀(ICP-OES 5300 DV，PerkinElmer)，測定方法為電感耦合等離子體原子發(fā)射光譜分析法[18]，Cl-，SO42-，NO3-由美國賽默飛離子色譜儀(Ion chromotography system，ICS-2100，Thermo)利用離子色譜法[19]測定，HCO3-采用電位滴定法[20]測定，測定儀器為德國Buerkle自動滴定器(Automatic Titrating Burette，The Company Buerkle )。δ18O，δD則采用高分辨率的直接吸收激光光譜法[21]進行測量，利用液態(tài)水同位素分析儀(DLT-100)測定，測試精度分別為±0.1‰和±1‰。

3 結果與分析

3.1 水化學特征

2) 地下水化學特征。賦存于第四紀沖積物中的地下水化學特征在區(qū)域上呈現(xiàn)明顯差別。測試結果顯示，從北部山區(qū)至平原區(qū)水化學類型由HCO3-Ca·Na型過渡到HCO3·Cl-Ca·Mg，至東部落轉化為HCO3-Ca·Mg型。各采樣點礦化度均較高，最大值位于黃土營(w16)，2013年7月測定值達1 251.69 mg/L。三個采樣時段均值分別為959 mg/L，983.31 mg/L，819.56 mg/L，除黃土營(w16)的pH為酸性外，其余點pH值7.23～7.54。年內(nèi)及年際水化學特征變化不大，間接反映含水層補給源的水質和水量變化較小。

區(qū)內(nèi)巖溶裂隙水賦存于府君山組灰?guī)r溶隙之中，蛤蟆泉一帶、溫莊和黃土營的巖溶裂隙水水化學類型均以HCO3·SO4-Ca型為主，至程莊為HCO3-Ca型，間有HCO3·SO4-Ca·Na型水，東部落則為HCO3-Ca·Mg型水。該層水礦化度在266.63～747.77 mg/L之間，均值分別為509.00，548.26，375.23 mg/L。pH在7.1～7.6之間，呈弱堿性。對比同一時段自北向南的水樣數(shù)據(jù)發(fā)現(xiàn)，蛤蟆泉(w3)-程莊(w9)-溫莊(w11)一線HCO3-、Ca2+濃度以及礦化度均呈現(xiàn)出逐漸升高的趨勢，至黃土營(w13)處出現(xiàn)拐點，向南至東部落村(w23)又緩慢上升，初步判定地下水流向為自北向南，見圖1A，B，C。

A，B，C為2012年7月、2013年7月和2013年10月地下水Durov 圖；a，b，c為A，B，C對應時段東宮河流域各水體δD—δ18O關系

圖1 不同時段地下水水化學特征及氫氧同位素組成綜合分析

位于研究區(qū)北端的李莊自流井w24，水質優(yōu)良，水位穩(wěn)定，根據(jù)打井深度及巖芯、巖性判定井水為基巖裂隙水，補給區(qū)來自盆地東北部陡峻的花崗巖出露區(qū)，裂隙水由東北向盆地內(nèi)部流動，補給水頭較高，在李莊及其周邊地區(qū)局部基巖裂隙含水系統(tǒng)水頭高于地表，形成自溢區(qū)。水化學類型為HCO3-Ca·Mg型。陰離子以HCO3-占優(yōu)勢，陽離子以Ca2+、Mg2+為主，其余離子濃度值均較低，礦化度均值為399 mg/L，枯水期礦化度略高于豐水期。經(jīng)調查發(fā)現(xiàn)，李莊村僅這一戶姜姓村民家為自流井，其鄰居打同樣深度的井卻未自流，說明該地區(qū)的基巖裂隙水具有強烈的不均勻性和各向異性的特點，其成因解釋見圖2。在盆地西北部山區(qū)的饅毛組泥頁巖中也賦存有少量裂隙水，如王家峪的w19，水化學類型為HCO3-Ca·Na型水，礦化度為550.53 mg/L，硝酸鹽氮及鉀、鈉離子濃度較高，水量貧乏，水質較差。

泉水在盆地中東部的東部落村集中出露。s1，s2為HCO3·SO4-Ca·Mg型水，s3，s4為HCO3-Ca·Mg 型水，兩類型水中各離子濃度測試結果相差不大。礦化度介于344.00～428.18 mg/L之間，三個采樣期均值分別為344.00，405.56，387.98 mg/L。pH在7.18～7.74之間，呈弱堿性。以2013年10月為例，礦化度略高于東部落村采集的巖溶裂隙水(w20，w21，w23)，遠小于第四系孔隙水w22，由此判定泉水主要來源于巖溶裂隙水，經(jīng)由薄層的第四系弱透水層出露地表。結合東部落標準剖面(圖3)，主要接受東北部山區(qū)大氣降水入滲補給的巖溶水，自北向南流，流經(jīng)東部落村附近，南側、西側受阻，頂、底板受限，水頭抬升接近地表，豐水期補給水量大，泉水集中出露，枯水期補給量少，水位降至地表以下。

圖2 李莊自流井成因分析示意圖

圖3 東部落泉成因示意圖

3.2 氫氧同位素組成

研究區(qū)降水中的D和18O組成的變化主要受季風氣候影響，得到當?shù)亟邓€方程為δD=8.069δ18O+9.141‰，(R2=0.598，n=35)。其中δD值介于-10.23‰～-2.78‰，平均值為-5.90‰，δ18O介于-75.10‰～-0.14‰，平均值為-38.49‰，兩者隨時間變化趨勢一致。

地表水中D和18O的組成取決于地表水補給源中的相應組成，并受到水面蒸發(fā)過程產(chǎn)生的分餾作用的影響。當?shù)睾哟◤搅鱽碓粗饕写髿饨邓⑸嫌蝸硭偷叵滤韧緩健?/p>

總的來看，研究區(qū)地表水中的D和18O組成呈以下特征：三個時期δD均值分別為-62.8‰，-57.63‰，-56.6‰，δ18O均值分別為-8.8‰，-8.0‰，-7.65‰，地表水豐水期較枯水期δD和δ18O貧化。2012年7月采樣期正值雨季，地表水測點基本落在大氣降水線上，直接接受大氣降水補給，且沒有經(jīng)歷較強的蒸發(fā)作用，符合當?shù)叵募緷駸岬臍夂颉Ｒ源藶榛鶞?，其余兩時段測試濃度略偏離大氣降水線，采樣點數(shù)量僅為3個，統(tǒng)計規(guī)律不明顯。同時除大氣降水補給外，在北部山區(qū)至程莊段，2013年7月w8實測地下水位高于東宮河水位，地下水補給地表水，且由北向南補給量逐漸減少，至東部落村西河段，根據(jù)公式(1)，(2) 估算地下水排泄到河流中的份額為f=0.092%，說明此處地下水向東宮河排泄量已微乎其微。具體估算步驟如下：

按照18O在水體中質量守恒的原理，其質量守恒方程[22]為：

CsQs=CgQg+Cb(Qs-Qg)

(1)

由公式(1) 推導出地下水的排泄量占河水流量的百分比公式[22]為：

(2)

式中：Cs——取樣點河水中的δ18O值；Cg——取樣點地下水中的δ18O值；Cb——河流上游來水中δ18O值；Qs——取樣點河水流量；Qg——地下水排泄量；f——地下水的排泄量占河水流量的百分比。

以2012年7月采樣數(shù)據(jù)為例，東部落村西側河流采樣點r6的Cs值為-8.79‰，其上游河水樣r4的Cb值為-8.85‰，地下水取樣點w22的Cg為-8.20‰，代入公式(2) 即得結果。

如圖1a所示，第四系孔隙水三個時段絕大多數(shù)采樣點接近當?shù)仄骄髿饨邓€，直接接受大氣降水入滲補給。另外，參考礦化度均值，第四系孔隙水是地下水中礦化度最高的，間接驗證了地下水的運移方向，即深層地下水(如巖溶水)向淺層地下水排泄的事實。采樣點中僅2013年7月黃土營的w16點偏離較大，與同在黃土營且井深相同的w15對比發(fā)現(xiàn)，兩者水化學類型相同，但各離子濃度w16均高于w15，總的礦化度甚至相差約一倍，由此推斷w16的來自深層地下水補給量較大，地下水位埋深較淺，蒸發(fā)濃縮作用強烈。

巖溶裂隙水在豐水期，采集水樣對應的δD，δ18O關系趨于一直線，并與大氣降水線呈角度相交，見圖1a，b。2012年7月的水樣δD—δ18O關系趨勢線與大氣降水線交角較小，由北向南氫氧同位素濃度逐漸增大，水化學類型由重碳酸型逐漸過渡到重碳酸—硫酸型，水化學作用以溶濾作用為主。而2013年7月的水樣與大氣降水線交角較大，蛤蟆泉—程莊—溫莊一線δD基本保持不變，而δ18O則逐漸增大?？菟冢瑯悠份^為分散地分布于大氣降水線下方，蛤蟆泉w3—程莊w9—溫莊w11—黃土營w13—東部落w20一線，與大氣降水線的偏移距離逐漸增大。整體上，無論豐水期還是枯水期自北向南水樣與大氣降水線偏差均呈現(xiàn)逐漸增大的趨勢，說明北部山區(qū)巖溶裂隙出露區(qū)為巖溶水主要補給區(qū)，隨著巖溶水向南部運移，補給路徑增大，大氣降水直接入滲補給強度減弱。

基巖裂隙水水樣同位素組成穩(wěn)定，集中于當?shù)仄骄髿饨邓€下方，且距大氣降水線位移量基本相等，反映基巖裂隙水補給源較為穩(wěn)定且補給路徑較短，受外界氣候條件干擾較小，與水化學特征論證結果一致。

泉水除s4(2013年7月樣品)外，圖1a，b，c中泉水采樣點基本落在巖溶水所在區(qū)域范圍內(nèi)，偏離大氣降水線的距離基本相同，進一步驗證了泉水的主要來源為補給穩(wěn)定的巖溶裂隙水。

4 結 論

東宮河流域厚層的府君山組灰?guī)r裂隙含水層為研究區(qū)主要含水層，巖溶水自北向南流，與柳江盆地區(qū)域巖溶水流動方向一致。盆地北部山區(qū)灰?guī)r出露區(qū)是該含水層的主要補給區(qū)，接受大氣降水入滲補給，豐水期主要在東部落村以泉水形式排泄，同時，野外實際測量和計算結果，豐水期從盆地北緣至東部落以西河段巖溶裂隙水以地下徑流的形式向地表水排泄，判定東部落泉類型為侵蝕上升泉。

基巖裂隙水的分布具有強烈的非均質各向異性的特點。盆地東北部陡峻的花崗巖山區(qū)的裂隙水，水質好，補給水頭最高，補給源穩(wěn)定、補給路徑較短，在李莊及其附近存在自流區(qū)。盆地西北部山區(qū)的饅毛組泥頁巖中裂隙水水質水量均較差。

第四系松散堆積物厚度不均，自西向東巖性由沖、洪積物過渡為沉積、殘積物，富水性逐漸變差，直接接受大氣降水入滲補給，同時接受下伏深層巖溶水補給，水流方向與巖溶水基本一致。

[1] 華北平原地下水環(huán)境演化[M].北京：地質出版社,2000.

[2] 趙玉華.秦皇島市水資源持續(xù)利用研究[J].農(nóng)業(yè)環(huán)境與發(fā)展,1998,15(4):1-3.

[3] 田級生.柳江水源地巖溶地面塌陷[J].水文地質工程地質,1994,21(4):52-54.

[4] 張風霞.超采地下水造成嚴重后果[J].地下水,1988,(3):190.

[5] 楊丙中,李良芳,徐開志,等.石門寨地質概況及教學實習指南[M].北京：地震出版社,1992.

[6] 馬金珠,黃天明,丁貞玉,等.同位素指示的巴丹吉林沙漠南緣地下水補給來源[J].地球科學進展,2007,22(9):922-930.

[7] McDonnell J J, Stewart M K, Owens I F. Effect of catchment-scale subsurface mixing on stream isotopic response[J]. Water Resources Research,1991,27(12):3065-3073.

[8] Hinton M J, Schiff S L, English M C. Examining the contributions of glacial till water to storm runoff using two-and three-component hydrograph separations[J]. Water Resources Research,1994,30(4):983-993.

[9] 蘇小四,林學鈺.銀川平原地下水循環(huán)及其可更新能力評價的同位素證據(jù)[J].資源科學,2004,26(2):29-35.

[10] Fórizs I, Berecz T, Molnár Z, et al. Origin of shallow groundwater of Csepel Island (south of Budapest, Hungary, River Danube): isotopic and chemical approach[J]. Hydrological Processes,2005,19(17):3299-3312.

[11] Guay B E, Eastoe C J, Bassett R, et al. Identifying sources of groundwater in the lower Colorado River valley, USA, with δ18O, δD, and 3H: implications for river water accounting[J]. Hydrogeology Journal, 2006,14(1/2):146-158.

[12] Weyhenmeyer C E, Burns S J, Waber H N, et al. Isotope study of moisture sources, recharge areas, and groundwater flow paths within the eastern Batinah coastal plain, Sultanate of Oman[J]. Water Resources Research, 2002,38(10): 2-22.

[13] Scholl M A, Ingebritsen S E, Janik C J, et al. Use of precipitation and groundwater isotopes to interpret regional hydrology on a tropical volcanic island: Kilauea volcano area, Hawaii[J]. Water Resources Research, 1996,32(12):3525-3537.

[14] Guglielmi Y, Mudry J. Estimation of spatial and temporal variability of recharge fluxes to an alluvial aquifer in a fore land area by water chemistry and isotopes[J]. Groundwater, 1996,34(6):1017-1023.

[15] 田級生.柳江盆地巖溶水富水規(guī)律剖析[J].石家莊經(jīng)濟學院學報,1985,8(4)：21-24.

[16] 李暉,蔣忠誠,王月,等.新疆地區(qū)大氣降水中穩(wěn)定同位素的變化特征[J].水土保持研究,2009,16(5):157-161.

[17] 郭曉軍,蘇鳳洹,洪勇,等.蔣家溝流域雨季降水中氫氧同位素特征分析[J].水土保持研究,2012,19(2):82-85.

[18] 阮桂色.電感耦合等離子體原子發(fā)射光譜(ICP—AES)技術的應用進展[J].中國無機分析化學,2012,1(4):15-18.

[19] 杜濤,王瑩,高超,等.離子色譜法測定飲用水中無機陰離子[J].沈陽師范大學學報:自然科學版,2011,29(2):260-263.

[20] 國家環(huán)境保護總局.水和廢水監(jiān)測分析方法[M].北京：中國環(huán)境科學出版社,2002.

[21] 黃一民.洞庭湖流域水循環(huán)中穩(wěn)定同位素的變化特征及其影響機制研究[D].長沙：湖南師范大學,2013.

[22] 張應華,仵彥卿,丁建強,等.運用氧穩(wěn)定同位素研究黑河中游盆地地下水與河水轉化[J].冰川凍土,2005,27(1):106-110.

Study on Groundwater Cycle in Donggong River Basin Using Data of Hydrochemistry, D and18O

SONG Yang1,2, CHI Baoming1,2, GU Hongbiao2, QIN Wenjing2, DU Liwei2

(1.InstituteofEngineeringMechanics,CEA,Harbin150080,China; 2.InstituteofDisasterPrevention,Sanhe,Hebei065201,China)

By using environmental isotopes and hydrochemical components as tracers to survey water cycling, this paper revealed the groundwater cycling characteristics in Liujiang Donggong River Basin. Based on the field survey, designing and assigning the location of sampling points, collecting and testing atmospheric precipitation, surface water and groundwater samples, comparing different period chemical characteristics of water samples, combining hydrogen and oxygen isotopic composition, the conditions of recharge, runoff and discharge in different kinds of groundwater were determined. Results showed that groundwater originating from Fujunshan group limestone aquifer with thick layer in the Cambrian system was the most significance of water supply in Donggong River Basin; karst water flow direction was from north to south, which was in line with the regional karst water flow of Liujiang Basin. It received precipitation infiltration recharge from the northern mountains. In the wet period the discharge mainly included two forms： one was the spring discharge through Quaternary aquitard in Dongbuluo village, the other was in the form of groundwater runoff discharge to Donggong River in the northern margin of the basin to the west of Dongbuluo village. The bedrock fissure water followed, and its distribution had the strong inhomogeneous anisotropic characteristics. Compared to the Man-Mao group shale in the northwestern mountains of basin, the water quality in granite fissure from the northeastern mountains of basin was good, and had the highest recharge head, stable recharge source, and shorter recharge route. Meanwhile, artesian zone existed in Lizhuang and its vicinity. The research results can provide the scientific basis for the development and utilization of local water resources.

groundwater cycle; Donggong River basin; hydrochemistry; hydrogen and oxygen isotope

2014-04-23

2014-11-14

中國地震局教師科研基金(20130103);河北省教育廳高等學校科學研究計劃項目(科技類青年基金)(QN20132004)

宋洋(1982—),女,吉林敦化人,講師,博士,主要從事水文地質、防災減災工程與防護工程研究。E-mail:songy0433@126.com

遲寶明(1957—),男,遼寧瓦房店人,教授,主要從事水文地質、防災減災工程與防護工程研究。E-mail:chibaoming@126.com

P641.3

1005-3409(2015)02-0090-06