洞庭湖入江水道斷面調(diào)整模式研究

余 蕾，李凌云，盧金友，王加虎,鄭 珊

(1.長(zhǎng)江科學(xué)院 水利部江湖治理與防洪重點(diǎn)實(shí)驗(yàn)室，武漢 430010； 2.河海大學(xué) 水文水資源學(xué)院，南京 210098；3.武漢大學(xué) 水資源與水電工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，武漢 430072)

洞庭湖入江水道斷面調(diào)整模式研究

余 蕾1,2，李凌云1，盧金友1，王加虎2,鄭 珊3

(1.長(zhǎng)江科學(xué)院 水利部江湖治理與防洪重點(diǎn)實(shí)驗(yàn)室，武漢 430010； 2.河海大學(xué) 水文水資源學(xué)院，南京 210098；3.武漢大學(xué) 水資源與水電工程科學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室，武漢 430072)

為研究洞庭湖入江水道斷面的調(diào)整規(guī)律，選擇城陵磯(七里山)水文斷面，根據(jù)1962—2010年的實(shí)測(cè)水沙數(shù)據(jù),運(yùn)用小波分析方法對(duì)流量、含沙量序列進(jìn)行多時(shí)間尺度分析，探明了洞庭湖入江水道水沙條件的周期變化規(guī)律；并選取城陵磯(七里山)枯水河槽斷面面積為特征變量，結(jié)合流量小波系數(shù)對(duì)水沙條件序列進(jìn)行時(shí)段劃分；基于滯后響應(yīng)原理，建立了適用于城陵磯河段斷面面積對(duì)水沙條件的分時(shí)段變化調(diào)整模型。結(jié)果表明，建立的模型可以較好地模擬洞庭湖入江水道斷面的調(diào)整規(guī)律，計(jì)算值與實(shí)測(cè)值擬合系數(shù)R2為0.77。

河床演變；小波分析；多時(shí)間尺度；滯后響應(yīng)；斷面面積

1 研究背景

河床演變泛指河流從河源至河口所流經(jīng)河谷的各個(gè)部分的形成和發(fā)展的整個(gè)歷史進(jìn)程[1]。在有沖有淤的演變過(guò)程中，河床會(huì)進(jìn)行自動(dòng)調(diào)整，建立一種與水沙條件或河床邊界條件相適應(yīng)的反饋機(jī)制。水沙條件是沖積河流河床形態(tài)調(diào)整的主要驅(qū)動(dòng)因子，水沙擾動(dòng)將會(huì)促使河床形態(tài)發(fā)生調(diào)整，不平衡輸沙是河床演變的實(shí)質(zhì)[2]。

大量研究表明，沖積河流河床演變受當(dāng)年和前期水沙條件的影響。錢意穎等[3]認(rèn)為河床形態(tài)是水沙多年累積作用的結(jié)果，將平灘流量與往年的來(lái)水來(lái)沙條件聯(lián)系起來(lái)。林秀芝等[4]通過(guò)分析渭河下游華縣站平灘流量與多年(1～5 a)平均水量滑動(dòng)值的關(guān)系，發(fā)現(xiàn)華縣站平灘流量Qb與當(dāng)年來(lái)水量Wi的0.7倍和往年來(lái)水量Wi-1的0.3倍加權(quán)平均水量相關(guān)最好，同時(shí)發(fā)現(xiàn)華縣站平灘流量Qb與當(dāng)年最大3日流量平均值Q3b及往年平灘流量Qi-1的加權(quán)組合值(權(quán)重分別為0.8,0.2)的相關(guān)性較好。Surian等[5]根據(jù)大量的實(shí)測(cè)資料總結(jié)出大部分河床再造過(guò)程是非線性的指數(shù)衰減函數(shù)的調(diào)整模式，在水沙發(fā)生擾動(dòng)時(shí)河床會(huì)立刻做出調(diào)整，迅速向新的平衡狀態(tài)靠近，河床再造的速度快。隨著時(shí)間的推移，其靠近的速度越來(lái)越慢?；谶@一原理，吳保生[6-7]提出了沖積河流河床演變的調(diào)整模式，即

dy/dt=β(ye-y) 。

(1)

式中：y為特征變量；ye為特征變量的相對(duì)平衡值；t為時(shí)間；β為特征變量的調(diào)整速率，可根據(jù)實(shí)測(cè)資料率定。

對(duì)式(1)進(jìn)行積分求解，即

y=(1-e-βt)ye+e-βtye0。

(2)

式中ye0為t=0時(shí)刻的平衡值。

在吳保生[6-7]研究的基礎(chǔ)上，李凌云[8]根據(jù)水沙條件對(duì)平灘流量的實(shí)際影響情況，調(diào)整了初始時(shí)刻平灘流量的平衡值Qe0權(quán)重，建立了黃河下游平灘流量的滯后響應(yīng)模型，模擬效果良好。

廖治棋[9]基于李凌云[8]的研究，改進(jìn)了滯后響應(yīng)模型的結(jié)構(gòu)和參數(shù)，并將其運(yùn)用到長(zhǎng)江中游的荊江河段上，初步揭示了長(zhǎng)江流域的滯后響應(yīng)現(xiàn)象。

總的來(lái)看，上述研究考慮了前期一定時(shí)期內(nèi)水沙作用對(duì)河床形態(tài)的累計(jì)影響，建立河床變量與水沙條件之間的響應(yīng)關(guān)系,可以很好地描述河床演變的調(diào)整規(guī)律。吳保生[6]提出的滯后響應(yīng)模型對(duì)河床演變的內(nèi)在機(jī)理進(jìn)行了更加全面的闡述，不足之處在于，適用范圍局限在我國(guó)的多沙河流,并且以往的研究均延續(xù)了逐年模擬的特點(diǎn)。但是這種方法對(duì)于河道斷面較為穩(wěn)定的河流，逐年模擬存在一定的困難。

河床演變是一個(gè)宏觀過(guò)程，本文擬在吳保生[7]建立的滯后響應(yīng)模型及李凌云[8]改進(jìn)的滯后響應(yīng)模型的基礎(chǔ)上運(yùn)用小波分析方法[10]對(duì)洞庭湖匯入長(zhǎng)江河道的斷面——城陵磯(七里山)水文斷面進(jìn)行水沙的多時(shí)間尺度研究，結(jié)合長(zhǎng)江流域“水多沙少”的特點(diǎn)，基于年平均流量序列的小波系數(shù)對(duì)水沙條件序列進(jìn)行時(shí)段劃分，建立模擬長(zhǎng)江流域河道斷面面積分時(shí)段變化的模型，研究斷面面積對(duì)水沙條件變化的響應(yīng)調(diào)整規(guī)律。

2 城陵磯站水沙的多時(shí)間尺度分析

城陵磯(七里山)水文斷面位于洞庭湖約8 km 處的出口水道上，下游4 km處是洞庭湖與長(zhǎng)江的交匯處，江湖關(guān)系十分復(fù)雜[11]。

本文擬選取城陵磯(七里山)水文站1962—2010年的流量、含沙量資料，利用小波方法分析城陵磯段水沙的多時(shí)間尺度規(guī)律。

2.1 小波分析方法

對(duì)于給定的Morlet子小波函數(shù)，即

ψa,b(t)=e-0.5t2eiw0t。

(3)

式中：w0為常數(shù);i為虛數(shù)。

若ψa,b(t)是由式(3)給出的子小波，給定的能量信號(hào)f(t)∈L2(R),其連續(xù)小波變換為

(4)

將小波系數(shù)的平方值在b域上積分，就可得到小波方差Var(a)，即

(5)

小波分析方法可以呈現(xiàn)水文時(shí)間序列的多時(shí)間尺度特征，從不同尺度上展現(xiàn)水文序列的波動(dòng)特性，反映水文過(guò)程的變化趨勢(shì)、預(yù)測(cè)未來(lái)的發(fā)展趨勢(shì)[13]。

2.2 城陵磯站水沙序列分析

將城陵磯(七里山)站1962—2010年的年平均流量、年平均含沙量資料進(jìn)行標(biāo)準(zhǔn)化處理后,把Morlet小波函數(shù)代入式(4)，分別計(jì)算出Morlet小波變換系數(shù)的實(shí)部。正負(fù)相位分別用實(shí)線和虛線繪制(如圖1)，實(shí)線代表流量處于豐水期或者含沙量處于偏高期，虛線表示流量處于枯水期或者含沙量處于偏低期。

圖1 城陵磯站流量和含沙量小波系數(shù)實(shí)部時(shí)頻分布Fig.1 Time-frequency distributions of discharge and sediment concentration of Chenglingji station

圖1(a)顯示，城陵磯站年平均流量存在明顯的年代際(>10 a)和年際尺度(<10 a)的周期性變化特征，存在17～30 a,7～17 a及7 a以下3類尺度的周期變化規(guī)律，正負(fù)相位交替出現(xiàn)。年代際變化中，17～30 a的尺度最為突出，流量經(jīng)歷了豐—平—枯—平—豐—平—枯—平—豐—平—枯 6次循環(huán)交替，中心尺度在28 a左右；7～17 a的時(shí)間尺度非常突出，流量經(jīng)歷了豐—平—枯—平—豐—平—枯—平—豐—平—枯—平—豐—平—枯—平—豐—平—枯—平—豐—平—枯—平—豐13次循環(huán)交替，中心尺度在11 a左右。年際變化中，7 a以下的尺度則出現(xiàn)更多循環(huán)。

圖1(b)顯示，城陵磯站年平均含沙量存在>30 a，14～25 a，8～14 a及6 a以下4類時(shí)間尺度。其中，>30 a的時(shí)間尺度正負(fù)相位雖交替出現(xiàn)，由于受資料限制，中心尺度無(wú)法確定；14～25 a時(shí)間尺度中，周期性較顯著，含沙量經(jīng)歷了低—高—低—高—低—高—低—高8次循環(huán)變化，中心尺度在17 a左右；8～14 a時(shí)間尺度也較為顯著，含沙量經(jīng)歷了高—低—高—低—高—低—高—低—高—低—高—低—高13次循環(huán)變化，中心尺度在7 a左右；另外，6 a以下時(shí)間尺度則出現(xiàn)更多循環(huán)。

由式(5)計(jì)算小波方差，并繪制圖2，小波方差的波峰可以確定序列的主要變化周期。

圖2 城陵磯站流量和含沙量小波方差Fig.2 Wavelet variances of discharge and sediment concentration of Chenglingji station

年平均流量的小波方差存在3個(gè)明顯的峰值(圖2)，依次為28，11，7 a的時(shí)間尺度。其中，28，11 a為主要周期，28 a的周期震蕩最強(qiáng)。含沙量的小波方差存在2個(gè)明顯的峰值(圖2)，依次對(duì)應(yīng)著7，17 a，還可判斷圖中存在>32 a的時(shí)間尺度，由于數(shù)據(jù)有限，有待進(jìn)一步證實(shí)，所以17 a為主要周期。流量、含沙量的主要周期和次要周期共同控制著城陵磯河段水沙的變化特征，而水沙一旦變化，河床形態(tài)即會(huì)做出響應(yīng)和調(diào)整。

3 城陵磯(七里山)河床斷面調(diào)整過(guò)程的模擬

3.1 特征變量的選取

斷面面積[14-15]是衡量河道橫斷面形態(tài)的一個(gè)重要物理量，斷面面積的變化直接反映了河床的沖刷或淤積，體現(xiàn)了河床的自動(dòng)調(diào)整作用。對(duì)于城陵磯(七里山)水文斷面，汛期來(lái)沙量較大，出流受長(zhǎng)江回水頂托的影響大，速度減緩，容易在汛期淤積[16]；非汛期流速較小，造床作用比較穩(wěn)定，所以選取城陵磯(七里山)枯水河槽(汛后10月份—次年3月份)的斷面面積作為本文的特征變量。

從宏觀角度來(lái)講，河床演變涉及范圍廣，歷時(shí)長(zhǎng)，影響因素眾多。以往的研究多是逐年模擬河床形態(tài)的變化，相比而言，分時(shí)段模擬則體現(xiàn)了在某一個(gè)時(shí)段內(nèi)水沙的變化比較均勻平緩或者急劇快速的情況下，河床形態(tài)的調(diào)整過(guò)程。

3.2 分時(shí)段模型的建立

水沙條件是沖積河流河道形態(tài)調(diào)整的主要驅(qū)動(dòng)力，長(zhǎng)江流域含沙量相對(duì)較低，徑流量在塑造河床過(guò)程中起著決定性作用。因此本文以流量小波系數(shù)對(duì)水沙條件序列進(jìn)行大致劃分。

根據(jù)第3.1節(jié)所述，城陵磯站年平均流量小波系數(shù)存在28,11,7 a的尺度。第1、第2主周期分別為28,11 a，主要周期決定了流量的變化特征，因此將第1、第2主周期疊加[17]，如圖3。

圖3 城陵磯站主要周期下流量小波系數(shù)Fig.3 Wavelet coefficients of discharge of Chenglingji station under different time scales

根據(jù)流量小波系數(shù)疊加圖，將水沙條件序列依據(jù)疊加的小波系數(shù)過(guò)程線大致分為10段： 1969—1972年、1972—1976年、1976—1979年、1979—1983年、1983—1987年、1987—1990年、1990—1994年、1994—1998年、1998—2003年、2003—2009年。1/4個(gè)周期與每一時(shí)段的長(zhǎng)度大致相當(dāng)，斷面面積隨水沙條件發(fā)生改變。從疊加的關(guān)系中還顯示了城陵磯站流量小波系數(shù)在未來(lái)2～3 a依舊為負(fù)值，可以預(yù)測(cè)年平均流量處于枯水期。

根據(jù)對(duì)水沙條件序列劃分的時(shí)段中，每一時(shí)段包括n年，n年之中水沙作用是不同的。水沙的“記憶效應(yīng)”[18]對(duì)河道形態(tài)的影響較大，根據(jù)時(shí)段末水沙作用強(qiáng)、時(shí)段初作用弱的特點(diǎn)，調(diào)整水沙作用權(quán)重，其時(shí)段初(i=0)水沙權(quán)重為1/(1+2+…+n)，時(shí)段末(i=n-1)水沙權(quán)重為n/(1+2+…+n)，第k時(shí)段水沙權(quán)重為(k+1)/(1+2+…+n)。并根據(jù)長(zhǎng)江流域“水多”的特點(diǎn)，引入月均流量最大值Qm。

特征變量平衡值的影響因子是汛期平均流量Qf和汛期來(lái)沙系數(shù)ξf，與之建立指數(shù)函數(shù)關(guān)系，來(lái)沙系數(shù)是指單位流量中含沙量的多少[6]。因此，依據(jù)上述分段方法和調(diào)整的水沙權(quán)重，城陵磯河段斷面面積平衡值的計(jì)算方法如式(6)。

(6)

式中：Aej為第j(j=1,2,…,10)個(gè)時(shí)段末面積的平衡值(m2)；Qfi為第j個(gè)時(shí)段中第i年汛期平均流量(m3/s)；ξfi為第j個(gè)時(shí)段中第i年汛期平均來(lái)沙系數(shù)(kg·s/m6)；Qmi為第j個(gè)時(shí)段中第i年月均流量最大值(m3/s);K為系數(shù);a,b,d均為指數(shù)，根據(jù)實(shí)測(cè)資料率定。

由此，城陵磯河段的斷面面積對(duì)水沙的響應(yīng)調(diào)整模式為

Abj=(1-e-βt)Aej+(1-e-βt)e-βtAbj-1。

(7)

式中：Abj為第j時(shí)段經(jīng)過(guò)n年水沙累計(jì)作用的斷面面積(m2)；當(dāng)j=1時(shí)，Ab0為初始年份汛中斷面面積(m2)。

在第1時(shí)段的水沙作用(n=4)下,斷面面積由A0調(diào)整至Ab1，除第1時(shí)段的初始面積的平衡值由初始年份1969年的汛中面積來(lái)確定，余下的每一時(shí)段的初始面積平衡值均為上一時(shí)段末的斷面面積Abj-1。此時(shí)，若水沙條件繼續(xù)發(fā)生變化，斷面面積Ab1開始調(diào)整，經(jīng)過(guò)多個(gè)時(shí)段(j=1,2,…,10)調(diào)整之后，斷面面積調(diào)整至Abj。

3.3 模型的應(yīng)用

根據(jù)建立的響應(yīng)調(diào)整模型，利用城陵磯河段1969—1985年、1990—2009年實(shí)測(cè)大斷面資料和水沙資料，擬合相應(yīng)的系數(shù)和指數(shù)分別為：K=968，a=0.16，b=-0.049 3，β=0.502，d=0.009 8，代入式(6)和式(7)可得適用于城陵磯站的斷面面積計(jì)算方法為

(8)

Abj=(1-e-0.502t)Aej+(1-e-0.502t)e-0.502tAbj-1。

(9)

對(duì)式(8)和式(9)計(jì)算值與斷面面積實(shí)測(cè)值進(jìn)行比較，可得相應(yīng)的模型計(jì)算精度指標(biāo)。圖4給出了小波分析分時(shí)段后，利用分時(shí)段的模型計(jì)算城陵磯站的斷面面積，并與實(shí)測(cè)值進(jìn)行對(duì)比。

圖4 斷面面積計(jì)算值與實(shí)測(cè)值時(shí)程曲線Fig.4 Time-history curves of section area between calculated and measured values

可以看出，分時(shí)段模擬基本能夠描述河床形態(tài)的調(diào)整過(guò)程，斷面面積在1970—1975年、1990—1996年變化較大。本文所建立的模型較好地模擬了斷面面積的總體變化趨勢(shì)。整體來(lái)看，斷面面積呈現(xiàn)增—減—增—減—增—減的趨勢(shì)，河床處于沖—淤—沖—淤—沖—淤的狀態(tài)。2000年以后，斷面面積變化相對(duì)比較平穩(wěn)。

圖5 斷面面積計(jì)算值與實(shí)測(cè)值擬合曲線Fig.5 Fitted curve of section area between calculated and measured values

圖5中，對(duì)斷面面積計(jì)算值與相應(yīng)的實(shí)測(cè)值進(jìn)行了比較?？梢钥吹?，當(dāng)每一時(shí)段考慮前期n年內(nèi)的水沙累計(jì)作用(其中n值根據(jù)小波分析來(lái)確定)，二者擬合系數(shù)R2=0.77，計(jì)算值與實(shí)測(cè)值相關(guān)關(guān)系良好，計(jì)算精度較高。

4 結(jié)論與展望

(1) 城陵磯河段流量、含沙量具有多時(shí)間尺度的特征。城陵磯站的流量序列存在7，11，28 a的變化周期，其中28 a為第一主周期；含沙量序列主要存在7，17 a變化周期，其中17 a為主要周期；主、次要周期共同決定著城陵磯站流量、含沙量的豐枯、高低的多時(shí)間尺度規(guī)律。水沙條件的改變將會(huì)導(dǎo)致河床的自動(dòng)調(diào)整。

(2) 長(zhǎng)江流域“水多沙少”，通過(guò)水沙的搭配關(guān)系、河道本身的特點(diǎn)，建立了斷面面積分時(shí)段變化的調(diào)整模型。本文對(duì)分時(shí)段模型的基本原理和計(jì)算方法進(jìn)行了闡述，并應(yīng)用在洞庭湖入江水道斷面—城陵磯站(七里山)，計(jì)算值與實(shí)測(cè)值的相關(guān)關(guān)系擬合系數(shù)R2達(dá)到0.77，效果較為顯著。

(3) 河床演變是一個(gè)宏觀過(guò)程，分時(shí)段模擬不失為一種新思路，可以從較大時(shí)間尺度上把握沖積河流河床演變的規(guī)律。上游來(lái)水來(lái)沙條件是河床形態(tài)調(diào)整的主要驅(qū)動(dòng)力，但是本文對(duì)于特征變量平衡值A(chǔ)e的計(jì)算只考慮了河段的水沙條件，后續(xù)研究可就河床邊界條件及人類活動(dòng)等因素的影響予以考慮，提高模型計(jì)算的精度。

[1] 謝鑒衡.河床演變及整治[M].北京:水利水電出版社,1990.

[2] 錢 寧，張 仁，周志德.河床演變學(xué)[M].北京:科學(xué)出版社, 1989.

[3] 錢意穎，吳 知，米粹俠,等．關(guān)于在多沙河流上修建水庫(kù)保持有效庫(kù)容的初步分析[C]∥黃河水利委員會(huì).水庫(kù)泥沙報(bào)告匯編.鄭州：黃河水利委員會(huì),1972.

[4] 林秀芝, 田 勇, 伊?xí)匝? 等. 渭河下游平灘流量變化對(duì)來(lái)水來(lái)沙的響應(yīng)[J].泥沙研究,2005, (5): 1-4.

[5] SURIAN N，RINALDI M.Morphological Response to River Engineering and Management in Alluvial Channels in Italy[J].Geomorphology,2003, 50(4):307-326.

[6] 吳保生.沖河流沖積河流平灘流量的滯后響應(yīng)模型[J].水利學(xué)報(bào), 2008, 39(6):680-687.

[7] 吳保生.沖積河流河床演變的滯后響應(yīng)模型——I模型建立[J].泥沙研究, 2008,(6):1-7.

[8] 李凌云. 黃河平灘流量的計(jì)算方法及應(yīng)用研究[D].北京:清華大學(xué), 2010.

[9] 廖治棋. 荊江河段平灘面積對(duì)水沙條件變化的滯后響應(yīng)研究[D]. 武漢:長(zhǎng)江科學(xué)院, 2014.

[10]王文圣，丁 晶，向紅蓮.小波分析在水文學(xué)中的應(yīng)用研究及展望[J].水科學(xué)進(jìn)展, 2002,13(4):515-520.

[11]盧金友，羅恒凱.長(zhǎng)江與洞庭湖關(guān)系變化研究[J].人民長(zhǎng)江, 1999, 30(4):24-27.

[12]郭文獻(xiàn)，夏自強(qiáng)，王鴻翔,等.近50年來(lái)長(zhǎng)江宜昌站水溫變化的多尺度分析[J].水利學(xué)報(bào), 2008,39(11):1197-1203.

[13]王文圣，丁 晶，李躍清.水文序列小波分析[M].北京:化學(xué)工業(yè)出版社, 2005.

[14]李凌云,盧金友,范北林,等.沖積性河流河床橫斷面形態(tài)研究與進(jìn)展 [J].長(zhǎng)江科學(xué)院院報(bào), 2014, 31(5):1-6.

[15]章運(yùn)超,李凌云,范北林,等.長(zhǎng)江沙市段斷面面積對(duì)水沙變異的滯后響應(yīng)研究 [J].長(zhǎng)江科學(xué)院院報(bào), 2016, 33(7):1-5.

[16]金升高，陳建湘，彭育元.城陵磯水文站水文斷面的代表性分析[J].水利水電快報(bào),1999, 20(18):20-23.

[17]張艷艷，鐘德鈺，吳保生.黃河平灘流量的多時(shí)間尺度現(xiàn)象[J].水科學(xué)進(jìn)展,2012, 23(3):302-309.

[18]梁志勇,楊麗豐,馮 普.黃河下游平灘河槽形態(tài)與水沙搭配之關(guān)系[J]. 水力發(fā)電學(xué)報(bào), 2005, 24(6):68-71.

(編輯：占學(xué)軍)

Model of the Adjustment of Inlet Section fromDongting Lake to Yangtze River

YU Lei1,2，LI Ling-yun1,LU Jin-you1,WANG Jia-hu2,ZHENG Shan3

(1.Key Laboratory of River Regulation and Flood Control of Ministry of Water Resources,Yangtze River Scientific Research Institute, Wuhan 430010, China; 2.College of Hydrology and Water Resources, Hohai University, Nanjing 210098,China; 3.State Key Laboratory of Water Resources and Hydropower Engineering Science, Wuhan University, Wuhan 430072,China )

In order to study the adjustment of discharge and sediment concentration of inlet section from Dongting Lake to Yangtze River, we take Chenglingji hydrological section as an example. According to measured data from 1962 to 2010, we carry out multi-time scale analysis of flow discharge and sediment concentration of Chenglingji station in association with wavelet analysis method to obtain the periodic changes of water-sediment load at the section. Then, we select the section area of low water level as characteristic variables, and use wavelet coefficient of discharge to divide discharge and sediment concentration sequences into different periods. On the basis of delayed response principle, we establish an adjustment model with different time intervals, which is suitable for the response of the water and sediment concentration in the area of the Chenglingji station of Qilishan. Results show that the model could well simulate the regularities of inlet section’s adjustment of Dongting Lake, and the fitted coefficientR2between calculated and measured data is 0.77, which demonstrates that the model has adequate accuracy.

riverbed evolution; wavelet analysis; multi-time scale; delayed response; section area

2015-09-22；

2015-10-23

國(guó)家自然科學(xué)基金項(xiàng)目(51339001,51209015)；湖北省自然科學(xué)基金項(xiàng)目(2014CFB329)

余 蕾(1990-)，女，湖北黃岡人，碩士研究生，主要從事水利工程方面研究，(電話)18061797809(電子信箱)1032602833@qq.com。

10.11988/ckyyb.20150808

2016,33(12):1-5

TV147

A

1001-5485(2016)12-0001-05