清河水库鲢鳙鱼类资源声学评估——回波计数与回波积分法的比较

王靖, 张超, 王丹, 蔺丹清, 李娜娜, 汤勇

王靖, 张超, 王丹, 蔺丹清, 李娜娜, 汤勇. 清河水库鲢鳙鱼类资源声学评估——回波计数与回波积分法的比较[J]. 南方水产科学, 2010, 6(5): 50-55. DOI: 10.3969/j.issn.1673-2227.2010.05.008
引用本文: 王靖, 张超, 王丹, 蔺丹清, 李娜娜, 汤勇. 清河水库鲢鳙鱼类资源声学评估——回波计数与回波积分法的比较[J]. 南方水产科学, 2010, 6(5): 50-55. DOI: 10.3969/j.issn.1673-2227.2010.05.008
WANG Jing, ZHANG Chao, WANG Dan, LIN Danqing, LI Nana, TANG Yong. Acoustic assessment of silver carp and bighead carp in Qinghe Reservoir: comparing echo-counting and echo-integrating methods[J]. South China Fisheries Science, 2010, 6(5): 50-55. DOI: 10.3969/j.issn.1673-2227.2010.05.008
Citation: WANG Jing, ZHANG Chao, WANG Dan, LIN Danqing, LI Nana, TANG Yong. Acoustic assessment of silver carp and bighead carp in Qinghe Reservoir: comparing echo-counting and echo-integrating methods[J]. South China Fisheries Science, 2010, 6(5): 50-55. DOI: 10.3969/j.issn.1673-2227.2010.05.008

清河水库鲢鳙鱼类资源声学评估——回波计数与回波积分法的比较

基金项目: 

国家自然科学基金项目 10774021

教育部留学基金项目“近岸岩礁鱼类单体三维目标强度的研究” [2008]890号

辽宁省教育厅计划项目 2009A170

详细信息
    作者简介:

    王靖(1981-),女,硕士研究生,从事渔业水声学研究。E-mail: wj3010@163.com

    通讯作者:

    汤勇, E-mail: tang@dlfu.edu.cn

  • 中图分类号: S932.4;S932.9+15

Acoustic assessment of silver carp and bighead carp in Qinghe Reservoir: comparing echo-counting and echo-integrating methods

  • 摘要:

    2009年7月下旬使用挪威Simrad公司EY60型(70 kHz)分裂式波束科学鱼探仪对清河水库进行了鲢(Hypophthalmichthys molitrix)和鳙(Aristichthys nobilis)的资源声学调查。使用声学处理软件Echoview 4.0进行声学回波数据的后处理,研究库区内鲢和鳙昼夜分布特性,尝试使用回波计数法与回波积分法进行资源密度的评估,并进行了对比分析,讨论适合于水库鲢和鳙的声学调查方法。结果表明,鲢和鳙分布于水库10 m附近水层,昼间较为集中,单体回波较少,夜间较为分散,可以识别的单体回波较多。深水区域个体较大且密度较低,较浅的区域个体较小且密度较大。声学计数评估更适合夜间进行。根据夜间数据,2种评估方法的结果基本一致,回波积分法获得鱼类密度为0.063 ind · m-3,回波计数法为0.049 ind · m-3,可以确定2种方法均适合鲢和鳙的调查。

    Abstract:

    We conducted an acoustic assessment of silver carp (Hypophthalmichthys molitrix) and bighead carp (Aristichthys nobilis) resources in Qinghe Reservoir in late July of 2009 by using the Simrad EY60 split-beam echosounder. With the software Echoview 4.0, we post-processed the data to analyze the distribution of the species during day and night. Besides, we assessed the resource density by using echo-counting and echo-integrating methods, which were then compared to find out the suitable acoustic investigation method for the two species. The results show that both species gather around 10 m in the reservoir and gather more usually during the day with less single echo but disperse at night with more single echo. Large fish gather in deep water area with relative lower density but small fish in shallow water area with high density. The echo-counting method is more suitable for night assessment, and according to the night data, the estimated density is 0.063 ind · m-3 and 0.049 ind · m-3 obtained by echo-integrating and echo-counting methods, respectively.

  • 硬骨鱼类的生长主要受生长轴,即生长激素(GH)/胰岛素样生长因子(IGF)的调控,GH为脑垂体分泌的单链肽类激素,通过生长激素受体(GHR)刺激肝脏等组织分泌IGF,进而发生生物学效应。GH居于生长轴的中心位置,在鱼类的个体生长中起关键性作用[1]。研究表明鱼类垂体GH基因的转录和转译水平受性别[2-5]、季节和生理状态[4-7]及外界理化因子[8-10]等因素的影响,且研究结果显示在一定条件下垂体GH基因转录或转译水平高的鱼类其生长速度也快。因此,GH成为研究鱼类生长性能的重要指标。

    高密度养殖是提高单位面积产量、实现更大经济效益的有效手段,但不可避免地影响养殖对象的生长速度和健康。而革胡子鲶(Clarias gariepinus)比其他养殖鱼类对高密度养殖的适应性强[11-16]。天津市德仁农业发展有限公司设施化养殖的革胡子鲶,在高密度下仍然能够快速生长(最高产量可达400 kg · m-3水体),推测与其内在的生长调控因子有关。文章通过分析不同养殖密度下革胡子鲶幼鱼的日增质量和垂体GH基因mRNA的表达水平,初步探讨革胡子鲶幼鱼在不同养殖密度下生长差异和生长差异的内在原因,以期为深入探究革胡子鲶的生长调节以及高密度养殖快速生长的机理奠定基础。

    试验用鱼为天津市德仁农业发展有限公司从福建漳州购进约1.5 g的革胡子鲶饲养而成,选取大小规格相近的幼鱼为试验对象,平均体质量为10~12 g。

    试验分组及饲养管理同DAI等[17]。即试验用革胡子鲶幼鱼取回后在天津农学院天津市水产生态及养殖重点实验室于水族箱(70 cm×50 cm×50 cm)中暂养2周后,选取发育良好、体质健康的个体[平均体质量为(30.71±0.89)g]用于该试验。试验鱼随机分为35 kg · m-3、65 kg · m-3、95 kg ·m-3和125 kg · m-3 4个密度组,每组3个重复,饲养于水族箱(70 cm×50 cm×50 cm)中,实际水体积为120 L。模拟天津市德仁农业发展有限公司水泥池饲养革胡子鲶幼鱼的实际养殖条件,每日10: 00和18: 00各换水一次,每次换水量为水体积的3/4;养殖水温保持在(27±1)℃。按体质量的2%每天8: 00和16: 00各投喂一次,所用饲料为天津市天祥水产有限公司生产的鲶鱼膨化配合饲料(粗蛋白35%、粗脂肪5%、粗纤维10%、灰分10%)。养殖试验为期60 d,试验期间水体积保持不变,每30 d取样一次。每天观察并记录鱼的摄食及死亡情况。

    分别于第30和第60天各分析一次。分析时停饲24 h,并以200 mg · L-1的MS-222麻醉,以水族箱为单位称量鱼体的总质量,同时计数鱼体数量,根据下列公式计算鱼的均质量和日增质量:

    鱼的均质量=试验结束时每箱鱼总质量(g)/鱼的数量(尾)

    日增质量=[试验结束时鱼的均质量(g · 尾-1)-试验开始时鱼的均质量(g · 尾-1)]/养殖天数

    养殖第30和第60天时,每个养殖密度各取经MS-222麻醉的鱼18尾(每个水族箱取鱼6尾),解剖取出垂体后立即放入液氮中,之后于-80 ℃保存备用。以β-肌动蛋白基因作为内参基因分析垂体GH基因mRNA的相对表达量。试验中以个体为单位,根据王晓梅等[18]报道的方法进行垂体总RNA的提取、cDNA第一条链的合成、GH基因ORF区段和β-肌动蛋白基因片段的扩增以及GH基因mRNA相对表达量分析。

    试验数据以“平均值±标准差”(X+SD)表示,数据采用SPSS 11.5统计软件进行单因素方差分析(One-Way ANOVA),LSD法多重比较各组间数值的差异显著性(P<0.05)。

    根据革胡子鲶在不同密度养殖30 d和60 d后体质量(表 1)计算鱼体在0~30 d以及30~60 d 2个养殖期的日增质量。结果表明,试验期间随着养殖密度的增大出现日增质量降低的趋势(图 1),统计分析结果显示,在0~30 d的养殖期间,35 kg ·m-3和65 kg · m-3密度组的日增质量差异不显著(P>0.05),95 kg · m-3和125 kg · m-3密度组的日增质量差异也不显著(P>0.05);但35 kg · m-3和65 kg · m-3密度组的日增质量显著高于95 kg · m-3和125 kg · m-3密度组(P<0.05);养殖30~60 d期间,35 kg · m-3和65 kg · m-3密度组的日增质量显著高于125 kg · m-3密度组(P<0.05),而其余各密度组间的日增质量差异均不显著(P>0.05)。

    表  1  养殖密度对革胡子鲶体质量的影响
    Table  1.  Effect of stocking density on body weight of juvenile C.gariepinus
    指标index 养殖密度/ kg·m-3 stocking density
    35 65 95 125
    初始体质量/g initial body weight 31.01±0.47 29.77±0.32 30.82±1.17 31.22±0.90
    第30天体质量/g body weight on 30th day 64.75±1.21a 60.79±0.43ab 57.40±3.38bc 55.38±2.83c
    第60天体质重/g body weight on 60th day 100.94±3.42a 96.64±3.84a 90.70±1.19b 82.87±2.71c
    死亡率/% mortality 3.69±0.76a 3.68±1.13a 5.13±0.36b 5.40±0.43b
    注:部分数据引自文献[17]。同一行肩标字母相同或无肩标字母者表示差异不显著(P>0.05),肩标字母不同者表示差异显著(P<0.05)
    Note:Some data in this table are taken from Reference [17].Values with the same letters or no letter within the same row are not significantly different (P>0.05);whereas those with different letters are significantly different (P<0.05).
    下载: 导出CSV 
    | 显示表格
    图  1  不同养殖密度革胡子鲶日增质量
    同一线上肩标字母不同者为差异显著(P<0.05),肩标字母相同者为差异不显著(P>0.05)
    Figure  1.  Daily weight gain of C.gariepinus reared at different stocking densities
    Values with different superscripts on the same lines are significantly different (P < 0.05), while those with the same superscripts are not significantly different (P > 0.05).

    经1%琼脂糖凝胶电泳检测表明提取的革胡子鲶脑垂体总RNA完整性较好(图 2显示了4个RNA样本的电泳结果)。用核酸蛋白测定仪(德国Eppendorf公司出品,BioPhotometer型)测定出每个RNA样本光密度比(OD260/OD280)为1.80~2.00,表明总RNA纯度较好,RNA样本可以用于后续的试验。

    图  2  革胡子鲶4个个体脑垂体总RNA电泳结果
    Figure  2.  Electrophoretic patterns of total RNA extracts from pituitary of four individuals of C.gariepinus

    图 3-a图 3-b分别为养殖第30和第60天时革胡子鲶部分个体GH和β-actin基因片段扩增后的电泳检测结果。

    图  3  养殖第30天(a)和第60天时(b)革胡子鲶个体GH和β-actin基因片段PCR产物电泳图
    泳道1~3、4~6、7~8和9~12分别为养殖密度35 kg · m-3、65 kg · m-3、95 kg · m-3和125 kg · m-3的各3个体GH和β-actin基因片段RCR产物电泳图,泳道M为DNA标准分子量(BBI公司)
    Figure  3.  Electrophoretic patterns of PCR products of GH and β-actin gene fragments of individuals on 30th (a) and 60th day nursing(b)
    Lanes 1~3, 4~6, 7~8 and 9~12 are electrophoretic patterns of PCR products of GH and β-actin gene fragments separately on three individuals at stocking densities of 35 kg · m-3, 65 kg · m-3, 95 kg · m-3 and 125 kg · m-3, respectively; Lane M is DNA Marker (BBI Company)

    利用GeneTools软件分析革胡子鲶个体GH(目标基因)和β-actin(内参基因)基因片段扩增后的电泳图谱,GH与β-actin基因片段扩增产物的信号强度比即为GH基因的相对表达量。

    图 4为养殖第30和第60天时不同养殖密度下革胡子鲶GH基因mRNA平均相对表达量的统计分析结果。养殖第30天时,各养殖密度间革胡子鲶GH基因mRNA相对表达量差异均不显著(P>0.05);养殖第60天时,革胡子鲶GH基因mRNA相对表达量在35 kg · m-3和65 kg · m-3密度组间差异不显著(P>0.05),95 kg · m-3和125 kg · m-3密度组间差异也不显著(P>0.05),但35 kg · m-3和65 kg · m-3密度组GH基因mRNA的相对表达量显著高于95 kg · m-3和125 kg · m-3密度组(P<0.05)。

    图  4  不同养殖密度革胡子鲶GH基因mRNA的相对表达量
    标有不同字母者为差异显著(P<0.05),无字母或字母相同者为差异不显著(P>0.05)
    Figure  4.  Relative expression of GH mRNA in pituitary of C.gariepinus reared at different densities
    Different letters above the bar indicate significant difference (P < 0.05);while the same letters or no letter above the bar indicate no significant difference (P > 0.05)

    鱼体的生长状况是水产养殖中备受关注的重要经济性状,而GH与鱼类的生长密切相关,在鱼类的个体生长及生长调节中起关键性作用。马细兰等[2]指出垂体GH基因在转录水平上反映GH的合成速度,是反映动物生长速度的有效指标,生长快速的鱼其垂体GH mRNA表达量高;并且应用半定量RT-PCR法分析得出尼罗罗非鱼(Oreochromis niloticus)雄鱼垂体GH mRNA表达水平明显高于雌鱼,认为这可能是尼罗罗非鱼雄鱼生长快于雌鱼的主要内在原因之一。性成熟的金钱鱼(Scatophagus argus)雌鱼生长快于雄鱼,其生长差异与垂体GH mRNA的水平也呈显著的正相关[3]。同样,生长快速的雌性黄鲈(Perca flavescens)垂体GH mRNA水平显著高于雄鱼,并且在生长快速的春季显著高于秋季[4]。性成熟的鲤(Cyprinus carpio)垂体GH mRNA表达也呈现出季节性的周期变化,生长较快的季节(水温18~20 ℃)垂体GH mRNA表达比生长相对较慢的季节(水温8~10 ℃)高3倍[6]。鲇形目鳠科的大鳍鳠(Mystus macropterus)脑垂体GH的含量分别在3月和8月出现峰值,与该鱼春季和秋季的快速生长一致[7]。该目鲇科的野生鲇(Silurus asotus)垂体GH的含量在3月和7月各出现一个高峰值,研究者认为这可能与野生鲇春季和产卵后生长加快有关;并且雌性鲇鱼垂体GH的周年含量显著高于雄鱼,与雌鱼生长速度显著快于雄鱼的现象相符合[5]。SHIN等[8]研究发现光的波长可影响观赏性鱼类——克氏双带锯齿盖鱼(Amphiprion clarkii)的生长速度和垂体GH mRNA水平、生长速率与垂体GH mRNA的水平呈正相关。李云等[9]研究显示饲料中添加一定量的半胱胺盐酸盐(1 mg · g-1饲料)可显著提高斜带石斑鱼(Epinephelus coioides)的体质量增长率和垂体GH的转录水平。马细兰等[10]发现促黄体素释放激素类似物(LHRH-A)可显著上调尼罗罗非鱼生长轴相关基因(GH、GHR和IGF)的表达,从而促进鱼类的生长。上述研究说明鱼体的生长速度与垂体GH基因的转录水平呈正相关。

    由于脑垂体合成GH后经过血液运送至靶器官,再通过GHR和IGF的介导进而达到促进鱼体生长的作用,因此鱼类的生长速度并不能完全通过GH水平来体现,还与GHR和IGF的水平有关。如银大麻哈鱼(Oncorhynchus kisutch)在禁食3周后,比正常投喂组肝脏GHR mRNA水平低,提示长时间禁食可引起肝脏GHR基因表达的下降[19]。SAERA-VILA等[20]报道金鲷(Sparus aurata)肝脏和脂肪组织GHRⅠ基因表达最丰富,但禁食后可引起GHRⅠmRNA水平的下降。斑点叉尾(Ictalurus punctatus)在投喂外源皮质醇或者禁食4周后,与对照组相比,肝脏GHR mRNA的丰度以及血清中IGF的浓度均显著下降,且肝脏GHR mRNA的丰度和血清中IGF的浓度呈正相关;投喂外源皮质醇后鱼的生长显著下降,禁食组鱼的体质量减少[21]。华益民和林浩然[22]报道鲤饥饿32 d导致生长受阻,肝组织IGF-ⅠmRNA表达水平下降;再投喂后,其生长和肝组织IGF-ⅠmRNA丰度都逐渐恢复,表现出鲤的生长状况与肝组织IGF-ⅠmRNA水平的同步变化。尼罗罗非鱼饥饿28 d内血清GH含量反倒显著升高,而血清IGF含量和肝胰脏IGF mRNA表达丰度显著下降,饥饿后再投喂,第14天时血清GH含量显著升高,继续投喂7 d血清GH含量又表现出显著下降,但仍与饥饿第7天时的GH水平相当,而血清IGF含量和肝胰脏IGF mRNA表达丰度在第21天恢复投喂过程中显著升高,并恢复到饥饿前水平[23]。鲮(Cirrhinus molitorella)在投喂不同蛋白水平的饲料时,在8周的养殖期内其特定生长率与肝组织IGF mRNA的表达量呈极显著的正相关[24]。但强俊等[25]在分析温度和饲料蛋白水平对吉富品系尼罗罗非鱼幼鱼生长的影响时则发现血清GH水平与鱼的特定生长率相关性较低。黄国强等[26]研究也发现褐牙鲆(Paralichthys olivaceus)幼鱼在低温胁迫10 d导致的生长减缓能在恢复至最适生长温度的30 d内获得完全补偿生长,但血浆GH和IGF水平与生长率不存在明显的相关性。

    该研究通过分析不同养殖密度下革胡子鲶幼鱼的生长以及垂体GH基因的表达,探讨养殖密度对二者的影响以及二者的相关性。该研究的数据统计分析结果显示,0~30 d养殖期内,35 kg · m-3和65 kg · m-3密度组间的日增质量以及95 kg · m-3和125 kg · m-3密度组间的日增质量差异均不显著,但前2个密度组的日增质量显著高于后2个密度组(图 1)。但养殖第30天时各密度组间GH基因mRNA的相对表达量并未受养殖密度的影响,即各密度组间的GH基因的相对表达量差异不显著(图 4)。30~60 d养殖期内,35 kg · m-3和65 kg · m-3密度组间的日增质量与95 kg · m-3密度组差异不显著,但显著高于125 kg · m-3密度组(图 1)。而养殖第60天时,GH mRNA表达量的分析结果显示35 kg · m-3和65 kg · m-3密度组显著高于95 kg · m-3和125 kg · m-3密度组(图 4)。上述结果表明在该研究中革胡子鲶幼鱼的日增质量与垂体GH基因mRNA相对表达量不存在紧密的相关性。结合前人的研究结果,在该研究中笔者推测,仅就生长轴GH-GHR-IGF而言,革胡子鲶在较短的养殖时间内(30 d),养殖密度增大时,由于空间不足、摄食不均等原因可能影响鱼体生存状态和营养水平,而导致GH下游的生长调控因子,如GHR和IGF的合成与分泌减少,最终使鱼体的日增质量降低。但是随着养殖时间的延长(60 d),养殖密度增大到一定程度时则会影响处于生长轴中心位置的因子——GH mRNA表达的下降,进而GH的合成与分泌减少,这可能成为影响鱼体生长速度的主要内在因素。对于GH-GHR-IGF这一生长轴在革胡子鲶生长中的调节作用以及三者的相互补偿关系仍需进一步深入探究。

    致谢: 该研究得到了铁岭市清河水库管理局给予的大力支持和帮助,在此谨向铁岭市银州区李永贵、清河水库管理局薛云刚和李可义致谢。同时感谢日本北海道大学水产学部饭田浩二博士提供使用 Echoview数据处理软件方面的支持。
  • 图  1   清河水库声学调查航线,虚线表示昼间航线,实线表示夜间航线

    Figure  1.   Survey routes in Qinghe Reservoir; dotted line indicates day route and solid line indicates night route

    图  2   使用Echoview软件表示的清河水库水坝附近昼间(a)、夜间(b)SV回波图

    Figure  2.   Echogram of SV around Qinghe Reservoir dam using Echoview software. day(a), night(b)

    图  3   夜间各调查航线单体计数值分布(a)及对应的平均SV分布(b)

    Figure  3.   Distributions of individual fish number(a)and corresponding average SV(b)of survey routes at night

    图  4   单体回波检测的TS最大值分布图

    Figure  4.   Distribution diagram of maximum TS value of single echo detection

    表  1   EY60声学测量的设定参数

    Table  1   Parameter setting for EY60 acoustic detection

    参数parameter 设定值setting value
    换能器频率frequency 70 kHz
    换能器指向角directivity angle
    脉冲宽度pulse length 0.256 ms
    发射功率transmitting power 500 W
    声速sound speed 1 499.26 m·s-1
    下载: 导出CSV
  • [1] 赵宪勇, 陈毓桢. 狭鳕(Theragra chalcogramma Palas) 目标强度的现场测定[J]. 中国水产科学, 1996, 2 (4) : 19-27. https://xueshu.baidu.com/usercenter/paper/show?paperid=1n6s0e30rj330m80aw3t0mg0sr157335&site=xueshu_se&hitarticle=1
    [2] 唐启升, 王为祥, 陈毓桢, 等. 北太平洋狭鳕资源声学评估调查研究[J]. 水产学报, 1995, 19 (1) : 8-20.
    [3] 李永振, 陈国宝, 赵宪勇, 等. 南海北部海域小型非经济鱼类资源声学评估[J]. 中国海洋大学学报, 2005, 35(2): 206-212. doi: 10.3969/j.issn.1672-5174.2005.02.007
    [4] 陈国宝, 李永振, 赵宪勇, 等. 南海5类重要经济鱼类资源声学评估[J]. 海洋学报, 2006, 28 (2) : 128-134. doi: 10.3321/j.issn:0253-4193.2006.02.016
    [5] 谭细畅, 史建全, 张宏, 等. EY60回声探测仪在青海湖鱼类资源量评估中的应用[J]. 湖泊科学, 2009, 21 (6) : 865-872.
    [6] 张慧杰, 杨德国, 危起伟, 等. 葛洲坝至古老背江段鱼类的水声学调查[J]. 长江流域资源与环境, 2007, 16 (1) : 86-91. https://xueshu.baidu.com/usercenter/paper/show?paperid=676157be3348d8ba97d5ab60f58b7a7e&site=xueshu_se
    [7] 陶江平, 陈永柏, 乔晔, 等. 三峡水库成库期间鱼类空间分布的水声学研究[J]. 水生态学杂志, 2008, 1 (1) : 25-33. https://xueshu.baidu.com/usercenter/paper/show?paperid=e4a400bdedd1f8ec9d0332bdab20c9e1&site=xueshu_se
    [8] 谭细畅, 夏立启, 立川贤一, 等. 东湖放养鱼类时空分布的水声学研究[J]. 水生生物学报, 2002, 26 (6) : 585-590. doi: 10.3321/j.issn:1000-3207.2002.06.002
    [9] 谭细畅, 李新辉, 赖子尼, 等. 青皮塘产卵场广东鲂繁殖群体的水声学探测研究[J]. 生态学杂志, 2008, 27 (5) : 785-790.
    [10]

    FOOTE K G, MacLENNAN D N. Comparison of copper and tungsten carbide calibration spheres[J]. J Acoust Soc Am, 1984, 75 (2): 612-616. doi: 10.1121/1.390489

    [11] 中华人民共和国国家质量监督检验检疫总局, 中国国家标准化管理委员会. 中华人民共和国国家标准《海洋调查规范, 第6部分: 海洋生物调查GB/T 12766.6-2007. 》[S]. 北京: 中国标准出版社, 2008: 89-97. https://xueshu.baidu.com/usercenter/paper/show?paperid=106g06q0b11h00g0577w04q0qg095323&site=xueshu_se
    [12]

    SIMMONDS J, MacLENNAN D. Fisheries acoustics: theory and practice[M]. 2nd ed. Oxford: Blackwell Science, 2005: 51-53. doi: 10.1002/9780470995303

    [13]

    HORNE J K, WALLINE P D, JECH J M. Comparing acoustic model predictions to in suit backscatter measurements of fish with dual-chambered swimbladders[J]. J Fish Biol, 2000, 57(5): 1105-1121. doi: 10.1111/j.1095-8649.2000.tb00474.x

图(4)  /  表(1)
计量
  • 文章访问数:  4672
  • HTML全文浏览量:  162
  • PDF下载量:  2718
  • 被引次数: 0
出版历程
  • 收稿日期:  2010-05-02
  • 修回日期:  2010-05-12
  • 刊出日期:  2010-10-04

目录

/

返回文章
返回