Applicability of six potential evapotranspiration estimation methods in Liaoning Province
CAO Yongqiang, LIU Mingyang, LI Yuanfei, XIAO ChunliuSchool of Geographical Sciences, Liaoning Normal University, Dalian 116029, China收稿日期:2019-03-15修回日期:2019-07-16网络出版日期:2019-10-29
基金资助: |
Received:2019-03-15Revised:2019-07-16Online:2019-10-29
作者简介 About authors
曹永强,男,内蒙古乌盟丰镇市人,教授,主要从事水文水资源研究E-mail:
摘要
关键词:
Abstract
Keywords:
PDF (14191KB)元数据多维度评价相关文章导出EndNote|Ris|Bibtex收藏本文
本文引用格式
曹永强, 刘明阳, 李元菲, 肖春柳. 不同潜在蒸散发估算方法在辽宁省的适用性分析. 资源科学[J], 2019, 41(10): 1780-1790 doi:10.18402/resci.2019.10.02
CAO Yongqiang.
1 引言
蒸发表示液态水从地面转化为水蒸气,最终进入大气的全过程的总称,主要包括各类水体、土壤和植物等水分的蒸发[1]。由于实际蒸散发(AET)观测数据的缺乏,人们通常利用潜在蒸散发(PET)来核算实际蒸散发。潜在蒸散发是指在水分供给充足前提下,某一下垫面达到的最大可能蒸散量[2]。因此在全球气候变化的大背景下,能够精准测量和估算潜在蒸散发,不仅有利于水资源管理和流域水循环规律等方面的研究,而且对于水库设计、灌溉用水调配和生态环境治理等方面均有十分重要的现实意义[3]。目前,Penman-Monteith(PM)模型是世界粮农组织(FAO)推荐的使用精度较高的PET估算方法[4,5]。一般关于PET时空分布的研究均是采用以PM模型计算结果为基础,再结合气象、地理位置等影响因素,来探究PET时空差异的原因。例如韦振峰等[6]使用PM公式对1981—2010年陕西潜在蒸散量进行研究,发现纬度跨度大和地形复杂是陕西潜在蒸散量空间差异的主要原因;杨志刚等[7]利用PM公式分析了雅鲁藏布江流域近50年潜在蒸散发变化,发现流域平均年、干季和湿季PET呈微弱的增加趋势,增幅均在3.00 mm/10a以下。然而,赵玲玲等[8]认为在众多潜在蒸散发量估算方法中PM法虽然具有很高的物理意义,但是需要众多的气象数据,不可避免地出现数据难以获取、确定参数困难等问题。同样国外****一直在寻找数据量需求少、方法简单、适用度高且满足精度的PET估算方法来代替传统的PM法,如基于辐射类的Mak法和PT法等,基于温度的Linacer法和Kharrufa法等,但这些PET估算方法往往具有特定的气候条件和地区性限制,直接从国外引进可能产生较大的误差[9]。若推广应用则需对经验系数或计算结果进行重新修正,如赵捷等[10]选择6种基于辐射的PET估算方法对黑河流域各站点潜在蒸散量进行估算,发现经过参数修正后,Mak法估算精度最高。
辽宁省是中国农业大省,作物种植面积约为3227.2 khm2,受极端气候变化影响,经常出现干旱、洪涝等自然灾害,导致农作物受灾减产[11]。潜在蒸散发异常容易影响该地区的水热平衡和干湿状况,可能导致区域内水资源短缺和灌溉调配计划紊乱,严重威胁到作物的正常生长发育。针对以上问题,选择合适的PET估算方法对认识辽宁省区域水文循环规律、指导农业生产和保护生态环境就显得格外重要。本文选择辽宁省历史气象资料,以PM法计算出1968年以来辽宁省潜在蒸散发为标准,根据误差评价指标和相关性指标,分析Rohwer(Roh)法、Hamon(Ham)法、Hargreaves(Har)法、Makkink(Mak)法和Priestley-Taylor(PT)法等5种潜在蒸散发估算方法在辽宁省的应用效果;同时对5种潜在蒸散发估算方法的相对误差空间变化特征进行系统分析,进而探讨空间异质性,以期为辽宁省潜在蒸散发量估算、水文循环规律和农作物灌溉管理等方面提供理论依据。
2 研究区域概况、数据来源及研究方法
2.1 研究区域概况
辽宁省地处中国东北地区最南端,位于118°53′E—125°46′E、38°43′N—43°26′N,属于温带大陆性气候,四季分明,平均气温在2~10 ℃之间波动,多 年平均降水量为500~1000 mm,平均太阳辐射总量在100~200 cal/cm2之间,平均日照时数2500 h左 右[12,13,14]。由于受到地理位置、地形和气候的影响,辽宁省降水分布不均匀,下垫面差异明显,导致省内旱涝灾害多发,潜在蒸散发作为气候变化和水汽循环中缺一不可的组成环节就显得格外重要。因此本文关于辽宁省潜在蒸散发的研究对实现本地区水资源可持续发展和生态文明建设有着现实的意义。2.2 数据来源及处理
选取1968—2017年辽宁省及周边27个气象站点(省内23个气象站点)的逐日气温、日照时数、平均风速、相对湿度及平均水汽压等气象资料。以上气象资料均来源于国家气象数据平台(http://data.cma.cn/),气象站点分布如图1所示。图1
新窗口打开|下载原图ZIP|生成PPT图1辽宁省气象站点分布图
Figure 1Spatial distribution of weather stations in Liaoning Province
2.3 研究方法
2.3.1 潜在蒸散发量估算方法介绍本文采用国际上较为公认的6种潜在蒸散发估算方法,分别为PM法、Roh法、Ham法、Har法、Mak法和PT法,每种方法都有特殊的形式和意义。PM法兼顾作物生理特征与空气动力学系数变化,适用于不同区域,但因其参数复杂、计算繁琐、基础数据要求较多具有一定局限性[15]。Roh法是一种估算自由水面蒸散发的方法,主要以空气动力学为原理,考虑风速和水汽压的影响,主要适用于干旱半干旱地区[16]。Ham法认为气温是蒸散发的主要驱动因子[17]。Har法是一种应用气温和太阳辐射资料的能量平衡的潜在蒸散发估算方法,主要适用于干旱半干旱草地[18]。Mak法将入射短波太阳辐射代替了净辐射同时忽略空气动力项,计算简单且应用较广,主要适用于湿润草地[19]。PT法同样忽略空气动力项,主要适用于湿润地区[20]。表1列出了各潜在蒸散法估算方法的计算公式,并对公式中各参数的物理意义进行了说明。
Table 1
表1
表1潜在蒸散发计算方法公式表
Table 1
方法名称 | 计算公式 | 说明 |
---|---|---|
PM法 | Δ为饱和水汽压与空气温度关系曲线的斜率,kPa/℃;Rn为到达地面的净辐射量,MJ/m2d;G为土壤热通量密度,MJ/m2d;γ为湿度计常数,kPa/℃;U2为地面以上2 m高处风速,m/s;es为空气饱和水汽压,kPa;ea为空气实际水汽压,kPa。各因子计算公式不再详述,具体参考文献[12]中的计算步骤 | |
Roh法 | α为经验系数,取0.44[16];U2、es、ea含义同PM法 | |
Ham法 | α为经验系数,取0.84[17];N为日照时数,h;T为日平均气温,℃;Pt对应的饱和蒸汽密度,g/m3 | |
Har法 | α为经验系数,取0.0135[18];λ为水汽化潜热,MJ/Kg;T为日平均气温,℃;RS为太阳辐射,MJ/(m2·d) | |
Mak法 | α为经验系数,取0.7[19];Δ、γ含义同PM法;RS、λ同Har法 | |
PT法 | α为经验系数,取1.26[20];Δ、γ、Rn含义同PM法 |
新窗口打开|下载CSV
2.3.2 精度评价标准
由于AET实际测量需要布设昂贵的大型称重式蒸渗仪,FAO专家建议将PM法计算结果看作参考潜在蒸散发,要求其他PET估算方法在新地区使用前应与其进行有效验证[21],并且中国许多****像夏军、赵玲玲、赵福年、廖晓芳等人[22,23,24]均证明PM法计算的潜在蒸散发与蒸发皿蒸发量相关性较高。因此本文采用平均绝对误差(MAE)、平均相对误差(MRE)以及确定性系数(R2)用于评价5种潜在蒸散发估算方法与PM法之间的精度,MAE、MRE越接近0,R2越接近1,精度越高。具体公式如下:
式中:n为样本总数;AE、E分别为其他方法计算的潜在蒸散发和PM法计算的潜在蒸散发;
3 结果与分析
3.1 辽宁省多年平均PET空间分布
利用ArcGIS对PM法计算出的辽宁省多年平均PET进行空间插值(图2)。从图中可以发现,辽宁省近50年平均PET为462.11~1109.58 mm,呈现出由西向东方向逐渐减少的趋势。从局部来看,位于西部的朝阳站、叶伯寿站以及彰武站一带的PET值较高,均值在1005.51 mm左右,这主要是由于辽西北地区紧邻内蒙古草原沙化带,常年日照强烈、风沙大、降雨稀少[25]。低值区出现在辽东的宽甸站、桓仁站以及清原站一带,均值在493.02 mm左右,这主要是由于这一带受植被覆盖面积大和黄渤海间接影响,湿度较大,降水较多,因此PET较低。图2
新窗口打开|下载原图ZIP|生成PPT图2辽宁省近50年平均PET空间分布图
Figure 2Average potential evapotranspiration (PET) spatial distribution map of Liaoning Province in the past 50 years
3.2 未修正经验系数时估算结果分析
3.2.1 辽宁省潜在蒸散发年值估算由于辽宁省冬季气温较低,使得某些估算结果出现负值,因此本文参考前人研究,在计算结果累加时将其剔除[10]。根据6种不同的方法,估算辽宁省近50年间PET随时间变化的曲线(图3),从图中可以看出,PM法计算结果总体呈现下降的趋势,PET值在595.07~804.3 mm之间波动,其中最小值出现在2015年,最大值则出现在1986年,这与吴霞等人[26]的研究成果基本一致。Mak法和PT法的估算结果与PM法的结果相差较小,且年际变化趋势大致相同,基本上同时达到波峰和波谷,表明这两种方法可以较好地反映出辽宁省PET的年际变化,而其他方法差异较明显,这与这些公式使用的经验系数选用固定值,没有考虑时间、下垫面等因素的变化有一定关系[27]。
图3
新窗口打开|下载原图ZIP|生成PPT图3未修正经验系数时6种估算方法年PET
Figure 3Results of six methods of estimation for annual PET when the empirical coefficient is not corrected
图4
新窗口打开|下载原图ZIP|生成PPT图4未修正经验系数时6种估算方法月PET
Figure 4Results of six methods of estimation for monthly potential evapotranspiration (PET) when the empirical coefficient is not corrected
根据精度评价方法,将5种估算方法计算出的PET值进行精度量化(表2)。从表中可以发现,PM法计算出的近50年辽宁省PET多年平均值为714.07 mm,其余5种计算结果在465.55~1019.59 mm之间,计算结果参差不齐。从精度评价指标来看,PT法和Mak法误差较小,为20.59%、23.27%,Ham法误差最大为42.77%。从R2指标来看,Mak法是5种估算方法中R2最高的,为0.87,其次为PT法。因此结合上文,在未修正经验系数时,Mak法、PT法与PM法计算结果较为相近,估算辽宁省PET年值的效果较好、精度较高。造成三者计算结果相近可能是由于辽宁省近50年PET对太阳辐射的变化较为敏感[28],这与Mak法和PT法的计算原理大致相同,故计算结果相近。
Table 2
表2
表2未修正经验系数时6种估算方法年PET效果
Table 2
方法名称 | PET年平均值/mm | 平均绝对误差/mm | 平均相对误差/% | 确定性系数 |
---|---|---|---|---|
PM法 | 714.07 | * | * | * |
Roh法 | 456.55 | 257.25 | 30.53 | 0.65 |
Ham法 | 1019.59 | 305.52 | 42.77 | 0.59 |
Har法 | 456.68 | 257.39 | 36.02 | 0.73 |
Mak法 | 880.26 | 166.19 | 23.27 | 0.87 |
PT法 | 852.86 | 138.79 | 20.59 | 0.78 |
新窗口打开|下载CSV
3.2.2 辽宁省潜在蒸散发月值估算
为了更好地精准分析辽宁省5种估算方法的应用效果和变化规律,本文对辽宁省PET月值进行研究(图4、表3)。从图4中可以发现,6种估算方法所得到的PET月值曲线变化趋势基本一致,均呈现出“两边低,中间高”的特征。Mak法、PT法与PM法的曲线最为接近,但夏季月值略高。Ham法曲线峰值出现时间(7月)比其他估算方法(6月)滞后一个月,这可能与其将温度作为唯一考虑的气象因子有关。从表3来看,除Ham法以外,其余估算方法的确定性系数均大于0.9。在5种估算方法中Mak法和PT法的相对误差均在20%以内,说明其误差较小。因此,在未修正经验系数时估算辽宁省PET月值,PT法与PM法最为接近,其次是Mak法。
Table 3
表3
表3未修正经验系数时6种估算方法月PET估算效果
Table 3
方法名称 | PET月平均值/mm | 平均绝对误差/mm | 平均相对误差/% | 确定性系数 |
---|---|---|---|---|
PM法 | 59.53 | * | * | * |
Roh法 | 38.02 | 19.05 | 38.83 | 0.90 |
Ham法 | 85.72 | 18.45 | 37.28 | 0.88 |
Har法 | 37.54 | 24.04 | 63.27 | 0.91 |
Mak法 | 74.43 | 9.79 | 18.95 | 0.97 |
PT法 | 69.10 | 6.19 | 13.96 | 0.99 |
新窗口打开|下载CSV
为了全面分析5种潜在蒸散发估算方法在辽宁省地区的适用性,仅仅从年尺度、月尺度研究还不够充分,为此本文进一步对各估算方法在空间尺度的适用性进行分析。以PM法计算的PET值空间分布为标准,对其余估算方法较PM法的相对误差(MRE)进行空间特征插值,进而分析出5种估算方法在空间尺度上的使用效果(图5)。
图5
新窗口打开|下载原图ZIP|生成PPT图5未修正经验系数时5种估算方法的MRE空间特征
Figure 5Mean relative error (MRE) spatial characteristics of five estimation methods when the empirical coefficients are not corrected
从图2和图5可以看出,未修正经验系数时5种估算方法在整个辽宁省适用情况和MRE的空间变化差异性较明显,对于开原、章党、沈阳和清原等辽中部偏北地区误差较大。5种方法中,PT法在整个区域适用性最好,平均MRE为20.59%,辽宁省南部尤其是沿海地区误差最小,辽西北地区误差较大。而Ham法平均MRE最大,空间上呈现出由北向南MRE越来越小的趋势,低值区出现在丹东、庄河辽东沿海一带和岫岩。结合各估算方法MRE空间变化特征和已有的辽宁省各气象因子敏感性变化规律分析[28],辽宁省各估算方法MRE与各自气象因子敏感性空间变化大致一致,说明各估算方法的适用范围差异与其主要气象因子空间分布和敏感性变化有关。未修正经验系数时,5种估算方法在辽东北地区估算精度较低,PT法估算精度优于其他估算方法。
3.3 辽宁省潜在蒸散发估算方法参数修正
辽宁省潜在蒸散发估算方法较PM法在计算结果、变化趋势和适用地区都存在较大差异性,因此需要对5种估算方法的经验系数进行修正,以期提高潜在蒸散发估算方法在辽宁省的适用性和可靠性。对比上文统计结果发现,5种估算方法计算的月PET值与PM法计算的回归效果更为显著,因此本文选择以5种估算方法计算的月PET值为自变量,PM法计算的月PET值为因变量,建立27个站点的线性回归方程,斜率即为经验系数的修正倍数,然后对经验系数进行回归修正,进而提出符合各站点特征的估算公式,最终达到提高精度的效果,具体修正结果见表4。Table 4
Table 4Correction of empirical coefficients of estimation methods for potential evapotranspiration (PET) in Liaoning Province
方法名称 | 原始经验系数 | 修正后平均经验系数 |
---|---|---|
Roh法 | 0.44 | 1.27 |
Ham法 | 0.84 | 0.56 |
Har法 | 0.0135 | 0.0142 |
Mak法 | 0.70 | 0.71 |
PT法 | 1.26 | 1.19 |
新窗口打开|下载CSV
3.4 修正经验系数后估算结果分析
3.4.1 辽宁省潜在蒸散发年值估算修正经验系数后,利用6种方法对辽宁省近50年间PET进行重新估算,制作年际变化曲线和精度量化表(图6和表5)。从图6中可以发现,经过修正后5种估算方法计算结果与PM法结果的年际变化趋势大致相同,基本上同时出现波峰和波谷。Roh法和Har法计算结果仍然低估PET值,而其他3种方法则高估PET值。从表5中可以发现,修正经验系数后,5种估算方法计算结果不同程度地接近PM法,精度明显提高,Mak法和PT法相对误差较修正前减小至10%以内,Ham法也缩小至17.46%。从R2指标来看,5种估算方法相关性均较修正前均有明显提高,全部在0.7以上,其中Har法、Mak法和PT法则在0.8以上,说明整体相关性较好。综上,在修正经验系数后,Mak法和PT法用于估算辽宁省年PET效果最好、精度最高。修正前后Mak法和PT法均对辽宁省年PET估算效果较好,再一次证明辽宁省PET对太阳辐射的变化较为敏感。
图6
新窗口打开|下载原图ZIP|生成PPT图6修正经验系数后6种估算方法年PET变化(1968-2017)
Figure 6Results of six methods of estimation after the correction of the empirical coefficients, 1968-2017
Table 5
表5
表5修正经验系数后6种估算方法年PET效果
Table 5
方法名称 | PET年平均值/mm | 平均绝对误差/mm | 平均相对误差/% | 确定性系数 |
---|---|---|---|---|
PM法 | 714.07 | * | * | * |
Roh法 | 567.95 | 146.12 | 21.00 | 0.74 |
Ham法 | 837.72 | 123.65 | 17.46 | 0.70 |
Har法 | 626.01 | 88.06 | 12.33 | 0.82 |
Mak法 | 776.26 | 62.19 | 9.72 | 0.88 |
PT法 | 774.87 | 60.08 | 8.17 | 0.86 |
新窗口打开|下载CSV
3.4.2 辽宁省潜在蒸散发月值估算
为检查经验系数的修正效果,对辽宁省PET月值进行重新估算,绘制出PET月变化趋势和精度量化表(图7和表6)。从图7中看出,修正经验系数后,6种方法所得到的PET月份变化趋势基本一致,均呈现出“低-高-低”的特点。除Roh法以外,其余4种估算方法得到的夏季PET值均略高于PM法。从表6中可以发现,修正经验系数后,5种方法估算的月PET值精度和相关性较修正前都有不同程度的提高。除Roh法和Har法外,其余3种方法估算PET月值略高于PM法。Har法估算的PET月平均值虽然与PM法较接近,但平均绝对误差和平均相对误差却是5种估算方法中最高的,这主要是由于Har法估算的各月PET值与PM法相差较大造成的。修正后,PT法的相对误差在5种估算方法中最小,为8.01%。因此,修正经验系数后,若考虑对辽宁省潜在蒸散发PET月值估算,5种估算方法中推荐使用PT法。
图7
新窗口打开|下载原图ZIP|生成PPT图7修正经验系数后6种方法月PET值
Figure 7Results of monthly potential evapotranspiration (PET) value of the six estimation methods after the correction of the empirical coefficients
Table 6
表6
表6修正经验系数后6种估算方法月PET值效果
Table 6
方法名称 | PET月平均值/mm | 平均绝对误差/mm | 平均相对误差/% | 确定性系数 |
---|---|---|---|---|
PM法 | 59.53 | * | * | * |
Roh法 | 47.05 | 12.84 | 24.68 | 0.90 |
Ham法 | 70.72 | 9.69 | 14.52 | 0.88 |
Har法 | 52.83 | 18.09 | 53.11 | 0.91 |
Mak法 | 65.24 | 5.70 | 12.94 | 0.97 |
PT法 | 64.63 | 4.77 | 8.01 | 0.99 |
新窗口打开|下载CSV
3.4.3 辽宁省潜在蒸散发方法空间分布
对修正后各站点PET估算方法的相对误差(MRE)进行空间特征插值,从而得出5种估算方法修正经验系数后在空间尺度上的应用效果(图8),对比图5和图8可以看出,经过修正经验系数后,5种估算方法在整个辽宁省区域的PET估算精度均有提高,但仍然存在对辽东北地区估算精度不高的问题,这可能是由于辽东北地区纬度较高,多为山林地,气候影响因素多样,下垫面复杂,导致本修正方法修正效果不明显。其中PT法在整个区域适用性最好,平均MRE为7.98%,而Har法精度较修正前提高最多。Roh法修正效果虽然有明显的提高,但是其适用性仍然是5种估算方法中最低的,这可能与其只考虑了空气动力项有关。对比修正前后,各估算方法MRE高值区均有所变化,这说明修正经验系数在辽宁省不同地区的修正效果并不一致,具有一定地域特点。因此综上,修正经验系数后,5种估算方法适用性均有不同程度的提高,其中PT法估算误差小于其他估算方法。
图8
新窗口打开|下载原图ZIP|生成PPT图8修正经验系数后5种估算方法的MRE空间特征
Figure 8Mean relative error (MRE) spatial characteristics of five estimation methods after the correction of the empirical coefficients
4 讨论
不同潜在蒸散发的估算方法所得到的结果差异不尽相同,这不仅受各站点地理位置、气候条件和下垫面情况影响,还与方法、参数的选择有关。Mak法和PT法无论修正前后相对于另外几种方法均有更好的适用效果,主要原因是PM法是由能量平衡学和空气动力学组成的计算公式,包括辐射能量项和空气动力项,辽宁省位于东北三省最南端,紧邻黄渤海,饱和水汽压差相对较小,空气动力项对潜在蒸散发贡献率较小,则辐射能量项为主导部分。而Mak法和PT法恰恰是以辐射项为主的估算方法,故适用效果较好。这些年辽宁省暖湿现象严重,导致潜在蒸散发变化巨大,因此能够准确快速估算潜在蒸散发至关重要。但现阶段水文模型的水文循环模拟中常出现其他简化的估算公式模拟结果与PM法模拟结果相当甚至效果更好的现象,陈桂华等[29]****认为PM法需要数据众多,且其中许多参数获取难度较大,对于许多传统概念型水文模型来讲,PM法对于数据输入过于敏感,可能出现不兼容的现象。因此,在一定条件下,修正简化后的潜在蒸散发估算公式不仅可以发挥其所需数据量少、计算方式简单等特点,同时还可以达到一定的精度要求,这为今后辽宁省区域水资源模拟与配置、农业规划管理提论了必不可少的理论支撑。
5 结论
本文以PM法计算的1968年以来辽宁省PET值为标准,根据误差评价指标和相关性指标,对Roh法、Ham法、Har法、Mak法和PT法等5种潜在蒸散发估算方法在辽宁省的适用性进行对比分析,可得出以下几点结论:(1)采用初始经验系数时,Mak法和PT法适用于估算辽宁省PET年值。对于PET月值的估算,PT法和Mak法与PM法估算结果的误差较小且相关性较高。空间上,5种估算方法在辽宁省的适用性存在明显的空间差异性,其中PT法估算精度优于其他估算方法,Ham法误差较大。
(2)修正经验系数后,5种估算方法各项数据均接近PM法的计算结果且年际变化趋势大致相同,其中Mak法和PT法适用性最好、精度最高。对于PET月值的估算,PT法是5种估算方法中最精确的方法,相对误差为8.01%。修正经验系数后,5种估算方法在空间尺度上适用性均有不同程度的提高,其中PT法估算误差小于其他估算方法,但仍然存在对开原、章党等辽东北地区精度不高的问题。
总的来说,PT法在辽宁省潜在蒸散发估算的应用中适用效果最好,其次是Mak法。本文仅以PM法月值为标准,对5种估算方法经验系数进行修正,但无论哪种估算方法在修正后仍会在某些区域出现误差较大的问题,因此选用合适的修正方法是下一阶段的主要研究重点。
参考文献 原文顺序
文献年度倒序
文中引用次数倒序
被引期刊影响因子
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 2]
[J].
[本文引用: 2]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 2]
[J].
[本文引用: 2]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 2]
[J].
[本文引用: 2]
[J]. ,
[本文引用: 2]
[J].
[本文引用: 2]
[J]. ,
[本文引用: 2]
[J]. ,
[本文引用: 2]
[J]. ,
[本文引用: 2]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 1]
[J].
[本文引用: 1]
[J]. ,
[本文引用: 2]
[J].
[本文引用: 2]
[本文引用: 1]
[本文引用: 1]