陳建生，何會祥，王　濤

(1. 河海大學(xué)巖土力學(xué)與堤壩工程教育部重點實驗室，江蘇 南京　210098；2. 河海大學(xué)巖土工程科學(xué)研究所，江蘇 南京　210098)

?

基于熵權(quán)-可變模糊集模型的堤壩滲漏探測

陳建生1,2，何會祥1,2，王濤1,2

(1. 河海大學(xué)巖土力學(xué)與堤壩工程教育部重點實驗室，江蘇 南京210098；2. 河海大學(xué)巖土工程科學(xué)研究所，江蘇 南京210098)

針對堤壩中集中滲漏通道探測問題，基于信息熵理論和可變模糊集理論，建立用于堤壩滲漏探測的熵權(quán)-可變模糊集模型。以南方某水庫為例，首先確定各探測指標(biāo)，采用客觀賦值法中的熵權(quán)法確定不同探測指標(biāo)的權(quán)重；然后計算各鉆孔水樣相對于降水、滲漏水和地下水的綜合相對隸屬度以及各個鉆孔的H值；最后依據(jù)相應(yīng)的判斷標(biāo)準(zhǔn)將鉆孔水樣進(jìn)行分類，確定鉆孔水的來源。結(jié)果表明，通過該模型確定此大壩存在4條集中滲漏通道，與前人結(jié)果一致。

信息熵；熵權(quán)法；可變模糊集；堤壩滲漏；滲漏通道探測

堤壩是防御洪水的重要工程設(shè)施，而滲漏是危害堤壩安全的一個重大隱患。對堤壩滲漏通道，尤其是集中滲漏通道位置的準(zhǔn)確探測在堤壩滲漏分析中顯得尤為重要。達(dá)西滲流理論應(yīng)用的關(guān)鍵在于滲透系數(shù)的準(zhǔn)確確定。對于均質(zhì)場，滲透系數(shù)的確定較為準(zhǔn)確，所得結(jié)果符合精度要求；但是對于非均質(zhì)、高度離散、不連續(xù)的裂隙巖體，所得結(jié)果在多數(shù)情況下與實際情況出入較大[1]。董海洲等[2-3]、王新建等[4-5]將熱力學(xué)第一定律中的內(nèi)能應(yīng)用于描述滲流場。流體內(nèi)能的傳遞不僅有對流方式，還有熱傳導(dǎo)，內(nèi)能是靠介質(zhì)傳遞的，介質(zhì)的幾何形狀對滲透性有很大影響，但是對于熱傳導(dǎo)系數(shù)的影響相對較小[1]。因此在滲流量計算以及集中滲漏通道的準(zhǔn)確探測方面，基于熱傳導(dǎo)理論建立的熱源法滲流模型比傳統(tǒng)的達(dá)西滲流理論具有得天獨厚的優(yōu)勢。

由于示蹤方法具有直觀的特點，因此其被廣泛應(yīng)用于堤壩滲漏探測中[6]。陳建生等[7]利用溫度示蹤方法，成功探測出陡河水庫左壩肩繞壩滲漏通道，并計算出繞壩滲漏量；段珣[8]利用地層中溫度的異常變化，探測出集中滲漏水從壩基處水平排水條帶匯集到壩后量水堰；陳建生等[9]、張令等[10]以環(huán)境同位素和水化學(xué)作為天然示蹤劑研究堤壩滲漏。然而目前的示蹤方法往往僅限于某單一方法，或幾種方法簡單結(jié)合，缺乏系統(tǒng)性。筆者擬基于水化學(xué)和環(huán)境同位素示蹤方法，綜合考慮分屬于水化學(xué)和環(huán)境同位素的各探測指標(biāo)，采用熵權(quán)法確定指標(biāo)權(quán)重，再運用可變模糊集理論計算各鉆孔水樣的綜合相對隸屬度，依據(jù)最大隸屬度原則對各鉆孔水樣進(jìn)行分類，并縱向分析級別特征值在壩軸線上的變化趨勢，以期綜合分析堤壩滲漏問題。

1　熵權(quán)-可變模糊集模型

1.1熵權(quán)法確定指標(biāo)權(quán)重

權(quán)重的確定是計算樣本相對隸屬度的關(guān)鍵，其恰當(dāng)與否直接影響結(jié)果的準(zhǔn)確性。目前權(quán)重的確定方法主要有主觀賦值法和客觀賦值法，主觀賦值法是指根據(jù)以前的統(tǒng)計數(shù)據(jù)由專家主觀判斷取得權(quán)值，這種方法與當(dāng)前樣本數(shù)據(jù)沒有關(guān)聯(lián)，和實際情況也脫節(jié)較大，客觀性較差，工作量大；客觀賦值法是根據(jù)樣本數(shù)據(jù)統(tǒng)計得出權(quán)值，此方法可以有效避免人為主觀因素帶來的誤差。

熵原指的是體系的混亂程度，是熱力學(xué)中的概念，后來被運用到信息論當(dāng)中，并產(chǎn)生了信息熵原理，其理論依據(jù)是：某一指標(biāo)的離散程度越高，其熵值越大；反之，熵值越小。當(dāng)由若干個指標(biāo)對事物進(jìn)行綜合評價時，對其中一個指標(biāo)而言，若各樣本值區(qū)別不大，則在綜合評價中該指標(biāo)的作用不大；反之，若該指標(biāo)的各個樣本值存在很大波動，即離散程度很高，則表明該指標(biāo)在綜合評價中具有很大影響?；谛畔㈧卦泶_定的權(quán)重客觀性較好，可信度較高，因此，筆者采用此方法確定指標(biāo)權(quán)重，其主要計算步驟見文獻(xiàn)[11-12]。

1.2可變模糊集和各探測指標(biāo)的相對隸屬度確定

可變模糊集理論是陳守煜基于Zadeh創(chuàng)建的模糊集合理論建立的。由于Zadeh創(chuàng)建的模糊集合隸屬函數(shù)存在絕對化、唯一化亦即靜態(tài)化的理論缺陷，因此，陳守煜運用自然辯證法中中介、差異、共維、兩極的概念以及客觀事物都存在的矛盾運動規(guī)律，建立了相對隸屬函數(shù)概念與模型，并首先將其運用于模糊水文水資源學(xué)[13-14]，而后被廣泛應(yīng)用于其他實際工程領(lǐng)域[15-16]，但在堤壩滲漏探測方面還鮮有應(yīng)用。由可變模糊集理論計算各探測指標(biāo)的相對隸屬度的主要步驟如下[17-18]：

圖1點x與區(qū)間[a,b]、[c,d]的位置關(guān)系
Fig. 1Position of point x relative to sections [a,b] and [c,d]

a. 計算各探測指標(biāo)的相對差異度D

A

(u)，A為論域上的一個模糊概念(事物、現(xiàn)象)，u為論域上的任意元素。依據(jù)可變模糊集定義，[

a

,

b

]為實軸上模糊可變集合的吸引域，[

c

,

d

]為包含[

a

,

b

]的某一上、下界范圍域，如圖1所示。

圖1中M為吸引域區(qū)間[a,b]中DA(u)=1的某一點值，從物理角度分析，一般M為[a,b]的中點值。點x是區(qū)間[c,d]內(nèi)的任意點量值，x落入M點左側(cè)時的相對差異函數(shù)模型為

(1)

x落入M點右側(cè)時,其相對差異函數(shù)模型為

(2)

式中：β——非負(fù)指數(shù)，通?？扇?，即相對差異函數(shù)模型為線性函數(shù)。

由式(1)、式(2)可以算得每個待測樣本中各探測指標(biāo)相對于不同級別的相對差異度DA(uth)(t=1,2,...,n，為指標(biāo)數(shù)；h=1,2,...,c，為級別數(shù))。

b. 依據(jù)模糊集合理論的余集定義計算各指標(biāo)對應(yīng)的相對隸屬度:

(3)

1.3計算綜合相對隸屬度和級別特征值

(4)

其中

式中：wt——指標(biāo)熵權(quán)重；p——距離參數(shù)；α——優(yōu)化準(zhǔn)則參數(shù)。

各樣本的級別特征值[18]：

(5)

其中

2　工 程 實 例

2.1工程概況

采用文獻(xiàn)[19]中的工程實例和數(shù)據(jù)對模型的正確性和實用性進(jìn)行驗證。水庫為一小型水庫，大壩為均質(zhì)土壩，壩頂長度為500 m、寬度為7.5 m，最大壩高為58.5 m，流域面積3.42 km2，水庫正常高水位81.8 m，相應(yīng)庫容為520×104m3。庫區(qū)內(nèi)廣布第四系松散堆積物，以殘坡積、沖洪積物為主，其次為坡洪積物。由于庫區(qū)位于排牙山背斜的南東翼，且屬較古老巖體，受多期構(gòu)造運動影響，構(gòu)造形跡較復(fù)雜，斷裂構(gòu)造發(fā)育。壩基及溢洪道開挖時發(fā)現(xiàn)多條較大斷層，對大壩的滲漏造成不可忽視的影響。大壩建設(shè)過程及建成后一直存在較大滲漏，總滲漏量達(dá)83.8 m3/d。

2.2指標(biāo)選取與分類標(biāo)準(zhǔn)

對于堤壩滲漏探測問題，指標(biāo)選取是否恰當(dāng)會直接影響探測結(jié)果的準(zhǔn)確性。探測指標(biāo)的合理選取可以依據(jù)示蹤法相關(guān)理論來進(jìn)行。示蹤法[20]按示蹤劑所屬類別不同分為水化學(xué)方法和環(huán)境同位素方法。由于地層中長期存在水巖相互作用、水土相互作用，因此地層中的水必然包含了能夠揭示地下水歷史以及所經(jīng)過的地層等方面的信息，水化學(xué)方法就是基于此來分析堤壩滲漏的部位、強度等。環(huán)境同位素方法的依據(jù)是大氣降水中稀有同位素存在的高程效應(yīng)、緯度效應(yīng)、陸地效應(yīng)、季節(jié)效應(yīng)，來源不同的水體會呈現(xiàn)不同的同位素特征，因此通過測定鉆孔水樣中的δD和δ18O可以分辨鉆孔水的來源[21]。綜上所述，選取pH、電導(dǎo)率、ρ(Cl-)、δD、δ18O共5項作為探測指標(biāo)。

由于實例中只給出了水庫溪水(近似為降水)和量水堰水(近似為地下水)的各探測指標(biāo)實測值，沒有給出降水、滲漏水和地下水各指標(biāo)的變化范圍?，F(xiàn)按照以下方法確定其變化范圍[18,20]：以量水堰水和水庫溪水的探測指標(biāo)測量值的平均值作為滲漏水的參考實測值，以水庫溪水和滲漏水探測指標(biāo)值的平均值作為降水探測指標(biāo)的上限和滲漏水探測指標(biāo)的下限，以滲漏水和量水堰水探測指標(biāo)數(shù)值的平均值作為滲漏水探測指標(biāo)的上限和地下水探測指標(biāo)的下限。由此可確定滲漏水探測指標(biāo)的變化范圍、降水探測指標(biāo)的上限、地下水探測指標(biāo)的下限。降水下限和地下水上限可按照滲漏水探測指標(biāo)的區(qū)間長度，從降水上限和地下水下限分別進(jìn)行外延得到。降水、滲漏水和地下水各探測指標(biāo)的變化范圍見表1。

表1　鉆孔水樣分類標(biāo)準(zhǔn)Table 1　Classification standard of borehole water samples

2.3計算過程

依據(jù)表1、文獻(xiàn)[17]可以確定可變模糊集合差異函數(shù)計算的各項參數(shù)[a,b]、[c,d]取值矩陣Iab、Icd：

(6)

(7)

M取[a,b]中的最優(yōu)點，即左端點a。計算各探測指標(biāo)的相對差異度、相對隸屬度；基于熵權(quán)法，得pH、電導(dǎo)率、ρ(Cl-)、δD、δ18O這5項探測指標(biāo)的權(quán)重依次為0.195 3、0.199 7、0.155 1、0.209 0和0.240 8；最后應(yīng)用文獻(xiàn)[18]中的模糊可變識別模型取優(yōu)化準(zhǔn)則(α=1，p=1)計算16個鉆孔水樣對各級別的綜合相對隸屬度并作歸一化處理。歸一化后的隸屬度ui(i=1,2,3分別對應(yīng)降水、滲漏水和地下水)和級別特征值H見表2。

表2　各鉆孔水樣的歸一化隸屬度和級別特征值Table 2　Membership degrees of borehole water samples after normalization and characteristic values for different grades of water samples

2.4結(jié)果分析

據(jù)表2、最大隸屬度原則可知：(a)K0-025～K0+175鉆孔中的水主要是滲漏水，為第1條滲漏通道位置；由于在K0-025、K0+025、K0+075中u1大于u3，因此降水占主要成分，即鉆孔中水主要來源于淺層庫水或地表水， K0+125、K0+175中u3大于u1，深層庫水占其主要成分；一條集中滲漏通道中滲漏水的主要來源卻不相同，因此可以斷定在K0-025～K0+175鉆孔所處位置應(yīng)該存在2條滲漏通道，分別在K0-025～K0+075附近和K0+125～K0+175附近。(b)K0+275～K0+325以及K0+425～K0+475這2段鉆孔水樣亦為滲漏水，并且各鉆孔u3均大于u1，深層庫水占其滲漏水的主要來源，因此K0+275～K0+325和K0+425～K0+475附近為該大壩第3條、第4條集中滲漏通道。

由表2的級別特征值H、級別特征值判斷準(zhǔn)則[22]可知這16個鉆孔水全為庫水滲漏，顯然這不符合實際情況?，F(xiàn)繪制H沿大壩軸線方向上的變化曲線(圖2)，縱向分析H值的變化：K0-025～K0+175、K0+275～K0+300、K0+425～K0+475這些鉆孔都近似處于曲線的極小值附近，因此可以判斷該大壩主要存在3條集中滲漏通道，分別位于K0-025～K0+175、K0+275～K0+300、K0+425～K0+475附近。綜合上述分析可知熵權(quán)-可變模糊集模型對于該大壩滲漏探測結(jié)果與前人結(jié)果具有一致性。庫區(qū)平面滲漏通道位置見圖3。大部分鉆孔水樣來源與前人分析相一致，只在K0+025、K0+075、K0+275、K0+475處的分析略有不同，這4個鉆孔的u1與u3相差非常小，可見滲漏水來源于淺層庫水或地表水和來源于深層庫水的量近似相等。

圖2　H變化曲線Fig. 2　Change curve of H value

圖3　庫區(qū)滲漏示意圖Fig. 3　Sketch of reservoir leakage

3　結(jié)　　論

a. 運用基于吸引與排斥思想的可變模糊集理論，結(jié)合信息熵，提出了適用于堤壩滲漏探測的熵權(quán)-可變模糊集模型，克服以往權(quán)值確定的主觀性，并將其用于工程實例，依據(jù)最大隸屬度原則和級別特征值2種判別方法對鉆孔水樣進(jìn)行分類，2種方法相互印證，使得結(jié)果更具有可靠性。

b. 探測結(jié)果顯示：該大壩在沿壩軸線方向上共存在4條集中滲漏通道，分別在鉆孔K0-025～K0+125、K0+175、K0+275～K0+325和K0+425～K0+475附近。這與前人研究成果大致相同，也證實了該模型的正確性和實用性。

c. 本文所提出的模型也存在一定的缺陷，比如模型中參數(shù)的取值以及吸引域和界限范圍域的偏差都會對結(jié)果產(chǎn)生某種程度上的影響，接下來會就參數(shù)變化對該模型的影響程度作進(jìn)一步的研究。

[1] 陳建生,董海洲,吳慶林,等. 虛擬熱源法研究壩基裂隙巖體滲漏通道[J].巖石力學(xué)與工程學(xué)報, 2005, 24(22): 4019-4024. (CHEN Jiansheng, DONG Haizhou, WU Qinglin, et al. Detection of leakage passage in fissure rock with assumptive heat source method[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(22): 4019-4024. (in Chinese))

[2] 董海洲,陳建生.利用孔中溫度場分布確定堤壩滲透流速的熱源法模型研究[J].水文地質(zhì)工程地質(zhì), 2003, 30(5): 40-43. (DONG Haizhou, CHEN Jiansheng. Model research of heat source method by using water temperature distribution in borehole to determine seepage velocity of dyke[J]. Hydro-geology and Engineering Geology, 2003, 30(5): 40-43. (in Chinese)) [3] 董海洲,陳建生.利用溫度示蹤方法探測基坑滲漏[J].巖石力學(xué)與工程學(xué)報,2004,23(12):2085-2090.(DONG Haizhou, CHEN Jiansheng. Study on groundwater leakage of foundation pit with temperature tracer method[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(12): 2085-2090. (in Chinese))

[4] 王新建,陳建生.堤壩集中滲漏溫度場探測模型及數(shù)值試驗[J].巖石力學(xué)與工程學(xué)報, 2006, 25(增刊2): 3794-3801. (WANG Xinjian, CHEN Jiansheng. Model for detecting of concentrated leakage in dam and dyke and its numerical experiment[J]. Chinese Journal of Rock Mechanics and Engineering, 2006, 25(Sup2): 3794-3801. (in Chinese))

[5] 王新建,陳建生.溫度探測土壩圓柱狀集中滲漏模型研究[J].水文地質(zhì)工程地質(zhì), 2006, 33(4): 31-36. (WANG Xinjian, CHEN Jiansheng. A study of temperature field tracer method used to detect tubular leakage in dam[J]. Hydrogeology and Engineering Geology, 2006, 33(4): 31-36. (in Chinese))

[6] 陳建生,董海洲.堤壩滲漏探測示蹤新理論與技術(shù)研究[M].北京：科學(xué)出版社, 2007.

[7] 陳建生,樊哲超,董海洲,等. 熱源法探測陡河水庫左壩肩繞壩滲漏通道[J].水文地質(zhì)工程地質(zhì), 2006, 33(5): 61-65. (CHEN Jiansheng, FAN Zhechao, DONG Haizhou, et al. Detecting the dam-detouring leakage passage in the left dam abutment of Douhe Reservoir with the method of heat source[J]. Hydrogeology and Engineering Geology, 2006, 33(5): 61-65. (in Chinese)) [8] 段珣.基于溫度示蹤法的堤壩滲漏探測研究[J].科學(xué)技術(shù)與工程，2014，14(9):141-143. (DUAN Xun. Study on detection dam leakage with temperature tracer method[J]. Science Technology and Engineering, 2014,14(9) :141-143.(in Chinese) )

[9] 陳建生,楊松堂,劉建剛,等. 環(huán)境同位素和水化學(xué)在堤壩滲漏研究中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報, 2004, 23(12): 2091-2095. (CHEN Jiansheng, YANG Songtang, LIU Jiangang, et al. Application of environmental isotope and hydrochemistry to leakage study[J]. Chinese Journal of Rock Mechanics and Engineering, 2004, 23(12): 2091-2095.(in Chinese))

[10] 張令,張一,黨雪梅,等. 隧洞涌水來源環(huán)境同位素及水化學(xué)判別研究[J]. 工程勘察, 2013(5):39-42. (ZHANG Ling, ZHANG Yi, DANG Xuemei, et al. Study on the sources of tunnel gushing water by using environment isotope and hydrochemical analysis[J]. Geotechnical Investigation & Surveying, 2013(5):39-42. (in Chinese))

[11] 吳景龍,孫亞麗,馬毅,等. 基于熵權(quán)的可變模糊集方法在水質(zhì)污染評價中的應(yīng)用[J].數(shù)學(xué)的實踐與認(rèn)識，2012，42(6):66-73. (WU Jinglong, SUN Yali, Ma Yi, et al. Application of variable fuzzy set theory method based on entropy weight to evaluation water pollution[J]. Mathematics in Practice and Theory, 2012, 42(6):66-73. (in Chinese))

[12] 孟憲萌,胡和平.基于熵權(quán)的集對分析模型在水質(zhì)綜合評價中的應(yīng)用[J].水利學(xué)報,2009,40(3):257-260. (MENG Xianmeng, HU Heping. Application of set pair analysis model based on entropy weight to comprehensive evaluation of water quality[J]. Journal of Hydraulic Engineering, 2009, 40(3):257-262. (in Chinese))

[13] CHEN Shouyu, FU Guangtao. A DRASTIC-based fuzzy pattern recognition methodology for groundwater vulnerability evaluation[J]. Int Assoc of Hydrolog Sci, 2003, 48(2):211-220.

[14] CHEN S Y, HOU Z C. Multi-criterion decision making for flood control operations: theory and application[J]. J of Amer Water Resour Assoc, 2004, 40(1):67-76.

[15] 陳守煜，韓曉軍.圍巖穩(wěn)定性評價的模糊可變集合工程方法[J].巖石力學(xué)與工程學(xué)報，2006,25(9)：1857-1861. (CHEN Shouyu, HAN Xiaojun. Engineering method of variable fuzzy set for assessment of surrounding rock stability[J].Chinese Journal of Rock Mechanics and Engineering,2006,25(9):1857-1861.(in Chinese))

[16] CHEN Shouyu, FU Guangtao. A fuzzy approach to the lectotype optimization of offshore platform[J]. Ocean Eng, 2003, 30(7):877-891.

[17] 陳守煜.工程可變模糊集理論與模型:模糊水文水資源學(xué)數(shù)學(xué)基礎(chǔ)[J].大連理工大學(xué)學(xué)報， 2005，45(2):308-312. (CHEN Shouyu. Theory and model of engineering variable fuzzy set:mathematical basis for fuzzy hydrology and water resources[J]. Journal of Dalian University of Technology, 2005, 45(2):308-312. (in Chinese) )

[18] 陳守煜,李敏.基于可變模糊集理論的水資源可再生能力評價模型[J].水利學(xué)報,2006,37(4):431-435.(CHEN Shouyu, LI Min. Assessment model of water resources reproducible ability based on variable fuzzy set theory[J]. Journal of Hydraulic Engineering, 2006, 37(4):431-435. (in Chinese))

[19] WANG Tao, CHEN Jiansheng, WANG Ting, et al. Entropy weight-set pair analysis based on tracer techniques for dam leakage investigation[J]. Nat Hazards, 2015,76(2):747-767.

[20] 王濤,陳建生,王婷.熵權(quán)-集對分析模型探測堤壩滲漏[J].巖土工程學(xué)報,2014,36(11):2137-2143. (WANG Tao, CHEN Jiansheng, WANG Ting. Entropy weight-set pair analysis(SPA) for the dam leakage detection[J]. Chinese Journal of Geotechnical Engineering, 2014, 36(11):2137-2143. (in Chinese))

[21] 巴金福,詹滬成,張博.判斷堤壩滲漏來源的環(huán)境同位素法[J].水利水電科技進(jìn)展,2014,34(增刊2):1-3.(BA Jinfu, ZHAN Hucheng, ZHANG Bo. Environmental isotope method to judge the source of leakage of dam[J]. Advances in Science and Technology of Water Resources,2014,34(Sup2):1-3.(in Chinese))

[22] 陳守煜,郭瑜.模糊可變集合與模型及在巖爆分級預(yù)測中的應(yīng)用[J].巖石力學(xué)與工程學(xué)報，2005,24(8):4604-4609. (CHEN Shouyu, GUO Yu. Model of variable fuzzy sets and its application to classification prediction of rockburst[J]. Chinese Journal of Rock Mechanics and Engineering, 2005, 24(8):4604-4609.)

Dam leakage detection based on entropy weight-variable fuzzy set model

CHEN Jiansheng1, 2, HE Huixiang1, 2, WANG Tao1, 2

(1.KeyLaboratoryofMinistryofEducationforGeomechanicsandEmbankmentEngineering,HohaiUniversity,Nanjing210098,China;2.GeotechnicalResearchInstitute,HohaiUniversity,Nanjing210098,China)

Aiming at detecting the concentrated leakage passage in a dam, an entropy weight-variable fuzzy set model for dam leakage detection was developed based on information entropy theory and variable fuzzy set theory. Using a dam in southern China as an example, the indices to be detected were determined first, followed by the weight of each index via the entropy weight method in the objective weight assignment method. Then, the comprehensive membership degrees of different borehole water samples to precipitation, leakage water, and groundwater, as well as the H value of each water sample, were obtained. Based on these values, the borehole water samples were classified and their origins were determined. The results show that there are four concentrated leakage passages in the dam, which accords with previous findings.

information entropy; entropy weight method; variable fuzzy set; dam leakage; leakage passage detection

1000-1980(2016)04-0358-06

10.3876/j.issn.1000-1980.2016.04.0013

2015-11-17

國家自然科學(xué)基金(51578212)；國家重點基礎(chǔ)研究發(fā)展計劃(973計劃)(2012CB417005) ；中央高?；究蒲袠I(yè)務(wù)費專項(2015B24614)

陳建生(1955—)，男，江蘇泗陽人，教授，博士，主要從事同位素水文、滲流理論科學(xué)應(yīng)用研究。E-mail: jschen@hhu.edu.cn

TV62+1

A