网刊加载中。。。

使用Chrome浏览器效果最佳,继续浏览,你可能不会看到最佳的展示效果,

确定继续浏览么?

复制成功,请在其他浏览器进行阅读

基于综合经济效益最优的光伏温室太阳辐射分配方法  PDF

  • 王医 1,2
  • 周伟绩 3
  • 赵易 1
  • 吴晓 1
1. 石河子大学水利建筑工程学院,石河子 832000; 2. 重庆大学建筑城规学院,重庆 400045; 3. 石河子大学兵团能源发展研究院,石河子 832000

中图分类号: S625.1

最近更新:2023-10-13

DOI:10.13300/j.cnki.hnlkxb.2023.05.029

  • 全文
  • 图表
  • 参考文献
  • 作者
  • 出版信息
EN
目录contents

摘要

针对光伏温室光伏发电和农业生产光照分配不合理导致整体收益较低的问题,提出了基于综合效益最优的光伏温室光伏发电与农业生产太阳辐射分配方法。首先,建立光伏温室的光伏发电收益、农作物收益、运行维护成本模型,其中,以农作物直角双曲线修正光响应模型作为光伏温室农作物产量计算方法,将光照与农作物产量和收益紧密结合;其次,以光伏温室综合效益最优为目标,以光伏发电光照强度、温室农作物生长光照条件为约束条件,建立光伏温室整体收益数学模型。以新疆地区某光伏温室实际运行数据为基础,对比分析多场景下光伏温室收益,结果显示,选择光伏温室农作物最小日光有效光照辐射值光照分配占比,光伏温室整体收益最大;考虑光伏温室改造与运行维护花费,在不改变原有结构的情况下,农作物光照饱和点对应的光照分配比例可作为光伏温室整体收益最大点。

光伏温室(photovoltaic greenhouse,PVG)是近几年设施农业领域的研究热点之一。光伏温室的光伏组件与农业温室大棚复用土地和建筑框架,既节约了建设成本,也为农业生产提供安全可靠绿色电能。目前有关光伏温室的研究,主要集中于光伏温室有利于农作物生长的环境控制方

1-3、适合于不同地区和外环境光伏温室种植的农作物类4-7以及使用ANSYS8和Ecotect9等软件及CFD10-12等分析方法建立光伏温室内环境模型,以达到最佳生长环境要求,使光伏温室内环境更适宜农作物生长等方面。

光伏发电与温室农业生产的能量均来源于光照,并且PVG中光伏温室光伏发电(photovoltaic power generation in photovoltaic greenhouse,PVPG-PVG)和光伏温室农业生产(agricultural production in photovoltaic greenhouse,AP-PVG)共用土地的关系,导致单位面积内太阳能辐射既要供给光伏用于发电又要供给温室用于农业生产。由于太阳辐射在PVPG-PVG与AP-PVG中产生的经济效益不同,在不影响农作物正常生长的情况下,合理分配PVPG-PVG与AP-PVG光照比例对提升PVG整体经济效益非常重要。基于此,Li

13研究了我国PVG中光伏发电“全额上网”方式及其利弊,指出该方法并不适用于PVG;Cossu14以欧洲典型商用PVG为例,分析了PVG光伏覆盖率对农作物生长的影响;Marucci15研究表明,PVG不同倾斜角对内部小气候的影响不同;Bambara16研究了半透明光伏组件在温室中的使用可能性,结果表明在成本降低23%的情况下,能够替代温室敷层;Yano17研究表明,PVG中光伏发电对农作物有遮阴效益,提出可利用光伏组件间距调节对农作物的遮阴效果;Perez-alonso18提出采用棋盘布局方法来布置PVG中光伏组件位置,能够使PVG农作物获得更多光照;Ezzaeri19研究了地中海南部地区PVG冬季和夏季小气候与番茄产量之间的关系,结果表明40%光伏组件覆盖情况下,光伏组件对番茄产量影响很小;王20研究了PVG中光伏组件以直线排布、交错排布和棋盘排布的方式排列,在发电量、温室内环境方面的优缺点,结果表明直线排布更适合PVG光伏组件,但没有给出PVG采用直线排布时光伏组件面积占比。上述相关研究表明,光照辐射对PVG整体效益影响主要聚焦于光伏组件的排布方法、覆盖率及类型等方面,鲜见考虑PVPG-PVG与AP-PVG光照分配问题的研究。

本研究以PVG整体经济效益最大为目标,以农作物生长内环境要求以及光伏发电光照要求等为约束,建立基于PVG整体经济收益最优的太阳辐射分配优化数学模型,以期为光伏温室设计提供理论参考。

1 材料与方法

1.1 PVPG-PVG年收益模型

为了计算简便,本研究做了以下假设:(1)暂不考虑地域、温度、海拔等气象和地理条件和损耗等因素对光伏组件发电量的限制和影响;(2)暂不考虑弃光限电,假设电网能够无限接纳PVG发电上网。基于此,综合考虑基于安装容量的PVG中光伏发电量计算方法和基于安装面积的PVG中光伏发电量计算方法的特点,采用年太阳光辐射度代替每天太阳光辐射度。PVPG-PVG年发电量可表示为:

E(x)=K1K2P(1-x)S (1)

其中,Ex)为PVPG-PVG年发电量,kW·h;K1为PVG中光伏组件的转换率;K2为PVG中光伏逆变器的转化效率;P为PVG所在地区的年太阳辐照量,MJ/m2,由单位换算3.6 MJ=1 (kW·h);x为PVG中农业生产太阳辐射分配比例;S为PVG南向屋顶面积总和,m2

基于“自发自用,余电上网”的运行方式,PVG电价在PVG建设之初已经与国家能源局和国家电网签订了《PVG并网收购电价协议》,明确规定了PVG中光伏发电“自发自用”电价和上网电价。假设PVPG-PVG“自发自用”部分年电量为E1x),上网部分的电量为E2x),PVPG-PVG产生经济效益的电能为上网部分,而“自发自用”部分的电能减小从电网下电电量,在电网下电部分一并核算。因此,若上网部分电价为Q2,则:

E(x)=E1(x)+E2(x) (2)
Mph(x)=E2(x)Q2 (3)

式(2)~(3)中:Mphx)为PVPG-PVG年收益,万元;Q2单位为元/(kW·h)。

1.2 PVG内环境调控电能费用模型

考虑到PVG农业生产需以电能为基础的内环境调控,电网为PVG的稳定运行起支撑作用,当PVPG-PVG功率大于光伏温室内环境调节负荷需求时,接纳剩余电能;当PVPG-PVG不足时,为PVG提供内环境调节电能。根据实地调研,PVG在并网处有计量装置,能够准确计量上网和下网电量。PVG年耗电量为E3x),PVPG-PVG年供给AP-PVG电能为E1x),暂不考虑电网为PVG做备用电源而产生的附加费用,则PVG农业生产内环境调节费用模型为:

Mgrid=(E3(x)-E1(x))Q3 (4)

其中:Mgrid为PVG内环境调控费用;Q3为下网电价,元/(kW·h)。

1.3 AP-PVG年收益模型

AP-PVG年收益与年产量密切相关,需要先计算年产量,而AP-PVG中光照辐射能量与农作物产量之间的关系非常复杂,因此,目前还没有确切的数学表达式能够表征二者之间的关系。光照辐射能量为农作物生长提供了基础条件,而农作物的产量是生长累计的结果,因此,可以采用农作物光响应模型近似模拟农作物光照与产量之间的关系。基于此,本研究提出了基于农作物直角双曲线修正光响应模型的AP-PVG年产量数学模型。采用农作物直角双曲线修正光响应模型表征PVG农作物光照与干物质累积的关系,进而得到AP-PVG年产量。直角双曲线修正光响应模

21式(5)所示。

Pn(I)=αIPmaxαI+Pmax-Rd (5)
I=K3λ×106Pxt (6)

式(5)~(6)中,PnI)为农作物净光合速率,μmol/(m2·s);α为光响应系数;I为农作物光合有效辐射,μmol/(m2·s);Pmax为最大净光合速率,μmol/(m2·s);Rd为农作物呼吸速率,μmol/(m2·s);λ取0.53,为光合辐射单位μmol/(m2·s)转化为W/m2的系数;K3为转化系数;t为年平均有效时间,h。

PVG中,温室农作物生长获得太阳辐射分配为Px,W/m2,将农作物呼吸速率进行单位转化之后,PVG农作物光照-产量模型如(7)式:

f(x)=44ηPn(I)Kprodtresp10-6 (7)

式(7)中,fx)为农作物产量,t;η为光合辐射有效占比,晴天取0.47,阴天取0.5;Kprod为产量累计系数,取5;tresp为响应时间,取3 600 s。

则AP-PVG年收益模型为:

Magri=f(x)Q4 (8)

式(8)中:Magri为AP-PVG年收益,万元;Q4为农业生产产品单价,取该产品年平均价格,万元/t。

1.4 PVG年收益模型

以PVG整体年收益最大为目标,建立基于太阳辐射分配的PVG整体收益数学模型,如式(9)所示。

M=Mph+Magri-Mgrid-Mcost (9)

其中:Mcost为PVG年平均成本,包括PVG的年建造成本、年运行和维护成本,元。

为了描述随着AP-PVG占比降低需要使用电能补光,将所需的成本折合为Mcost,表达式为:

Mcost=KcostS[1+(x-x0)2] (10)

其中,Kcost为基数维护费用,取0.002万元/m2S为PVG面积,m2x0为农作物光饱和点时AP-PVG占比。

1.5 约束条件

PVG年收益模型中包含光伏发电和农作物生长两方面。因此,约束条件中包含光伏发电约束和农业生产约束。

1)光伏发电约束条件。PVG的获得的总辐射P应该满足光伏发电的基本条件即:

PminP (11)

其中,Pmin表示光伏电站能够并网的最低辐射功率,W/m2

2)农业生产约束条件。PVG农业生产需要满足农作物生长的光照、温度及湿度要求。温湿度能量要求是由下网电量转换而来,光照是由太阳光照提供。则PVG需满足农作物日光有效光照辐射(photosynthetically active radiation,PAR)条件为:

PARminPxPARmax (12)

其中:PARmin和PARmax是种植农作物的最小及最大日光光照辐射值,μmol/(m2·s)。

1.6 计算流程

基于综合效益最优的光伏温室太阳辐射分配方法计算流程如图1所示。将公式(1)~(8)和式(10)带入式(9),得到基于综合效益最优的光伏温室太阳辐射分配数学模型,将算例参数带入公式进行计算,并分析结果是否满足边界条件,并再次对模型进行调整,最终得到计算结果。

图1  辐射分配方法计算流程

Fig.1  Calculation flow of radiation distribution method

1.7 试验条件

以中国新疆某多晶硅连栋光伏温室为试验对象,其主要运行参数如表1所示。光伏日光温室的具体参数:连栋光伏温室长为68 m,跨度为7 m,脊高3.1 m。采用的多晶硅光伏组件长度为1.65 m,宽0.992 m,额定功率为250 W,正南方位,安装倾斜角为38º。AP-PVG种植的农作物为芹菜(Apium graveolens L.),其PAR转化为光照约束条件,光照占比为[0.2,0.8],光饱和点时x=0.6。

表1  光伏温室主要参数
Table1  Main parameters of PVG
参数 Parameters数值 Value
K1 0.17
K2 0.98
P/ (MJ/m2) 18.93
t/h 800
E2/(kW·h) 8 432.4
E3/(kW·h) 3 434.6
Q1/[元/(kW·h)] 0.42
Q2/[元/(kW·h)] 0.78
Q3/[元/(kW·h)] 0.39
Q4/(万元/t) 0.2
PARmax 1 200
PARmin 650

2 结果与分析

2.1 AP-PVG年收益模型求解

采用数学计算工具Matlab编写计算程序,考虑约束条件,对AP-PVG年收益模型进行求解,计算结果如图2所示。由图2可知:PVG综合效益与AP-PVG光照分配占比成反比关系,随着农业生产光照分配占比增加,收益逐渐降低。造成此种趋势的主要原因是农作物价格太低,导致农业生产效益较低,在最大日光有效光照辐射时达到最低。这与目前“光伏温室废弃农业生产,仅用于发电”的现状密切吻合。相比于农业生产,当PVPG-PVG光照占比达0.8时,PVG整体效益达4.46万元,相比于AP-PVG占比0.8时整体收益(3.06万元)约高出50%。从整体收益方面来看,选择PVG农作物最小日光有效光照辐射值PARmin时,光照分配占比作为PVPG-PVG和AP-PVG光照分配比例,PVG整体收益最大。

图2  AP-PVG光照分配占比对PVG综合效益的影响

Fig.2  Relationship between the economic benefits of PVG and the proportion light distribution of AP-PVG

2.2 不同农作物对PVG整体收益的影响

实际生产中需要考虑PVG农作物最小日光有效光照辐射时的补光花费以及补光相关设施费用投入。图3为6种(分别编号为A、B、C、D、E和F)不同作物的收益与农业生产光照分配占比的关系图,其中已经考虑补光对其整体收入的影响农作物价格按照2 000元/t计算。由图3可知:AP-PVG种植不同农作物时,在光照占比趋于两端时,整体价格差别较大,越趋于中间,整体价格基本相等。造成这种现象的主要原因为:农作物光照分配比例在光饱和点附近时,PVG整体能耗最小,运行维护费用最小。而且,农作物生产光照分配占比越趋近于最小日光有效光照辐射值,单位收益增幅越是缓慢,表明在趋近最小日光有效光照辐射值,PVG投资越大,虽然整体收益仍呈现增长趋势,但增幅下降较大。考虑实际光照分配的可操作性,以及“双碳”目标下节能减排的重要地位,PVG整体收益最优点选择农作物光饱和点对应的农业生产光照分配比例最为合适。

图3  不同农作物对整体收益的影响

Fig.3  Impact of different crops on overall economic benefits

2.3 农作物价格对PVG整体收益的影响

农作物价格对PVG整体收益的影响如图4所示。由图4可知,在生产和管理维护成本基本不变的情况下,农作物价格上涨,PVG整体收益呈现上涨趋势,且增幅几乎相等。随着农作物价格增长,PVG整体收益由光照分配占比在最小日光有效光照辐射值比光饱和点的27.5%缩减至14.2%,表明光照分配比例选择光饱和点时,农作物价格越高,PVG整体效益有显著提升。

图4  农作物价格对PVG整体收益的影响

Fig.4  Impact of crop prices on overall economic benefits of PVG

2.4 农作物混种对PVG整体收益的影响

目前,超大型温室PVG农作物混种依然存在。通过分析PVG综合经济最优模型可知,农作物混种仅影响了农作物光照-产量模型和能量流动模型,约束条件及计算方法均与原模型一致。农作物混种PVG第i种农作物光照-产量模型如式(13)所示:

fi(x)=44ηPn,i(I)Kprodtresp10-6 (13)

农作物混种PVG农业生产年收益模型如式(14)所示:

Magri'=i=1nfi(x)Qi (14)

Magri'为农作物混种PVG农业生产年收益,万元;Qi为农作物混种PVG中各种农业生产产品单价,万元/t,取该产品年平均价格。由分析可得:光伏发电价格波动及电网下电价格波动对PVG经济收益影响与种植单一农作物相同,而农作物价格波动对经济收益影响较大。图5以2种农作物混种为例,分析农作物价格波动对农作物混种PVG经济收益影响。由图5可知,种植单一农作物,PVG整体经济效益会随着农作物价格增加而增加;相比于种植单一农作物,农作物混种整体经济收益与各农作物价格密切相关;相比于单一种植,混合种植对农作物价格波动的抗风险能力有所提升;农作物混种更能够匹配其饱和光照点,降低建造和维护费用,增加整体收益。需要注意的是,混种的农作物需要光照饱和点匹配,以便PVG整体光照分配协调。

图5  农作物混种对PVG整体收益的影响

Fig.5  Impact of mixed planting of crops on overall economic benefits of PVG

3 讨论

目前,光伏温室设计时光伏组件铺设面积基于经验值,而未经过科学论证,可能导致光伏温室内农作物光照不够充分,严重影响其产量和质量,亦可能导致光伏组件面积铺设较少,严重影响光伏发电效益。光伏温室中PVPG-PVG与AP-PVG复用土地,共同分享太阳辐射,而太阳辐射量直接影响PVPG-PVG发电量,进而影响收益。对于AP-PVG,太阳辐射影响其农作物产量,假设农作物价格不变,则太阳辐射对收益影响较大。PVPG-PVG与AP-PVG收益主要决定因素为太阳辐射,因此,合理分配太阳光照是提高PVG整体经济效益的关键。“自发自用,余电上网”PVG运行模式中,光伏发电量可应用于农业生产,增加农作物产值,从而增加农业生产收益。针对上述问题,本研究提出了基于综合效益最优的光伏温室光伏发电与光伏温室农业生产太阳辐射分配方法。

PVG年产出收益包括:PVPG-PVG年收益、PVG并网年收益(关口计量为正,表明PVG上网电量大于下网电量,收益为正;反之,PVG上网电量小于下网电量,收益为负)、AP-PVG年收益及年投资分摊成本。本研究以PVG年产出收益最大为目标,建立PVG太阳辐射分配模型,优化PVPG-PVG与AP-PVG太阳光辐射占比,并基于直角双曲线修正光响应模型得到的农作物产量计算方法,为光伏温室收益计算提供了基础条件,多场景计算结果表明,农作物最小日光有效光照辐射值PARmin时的光照分配占比为PVG整体收益最大值,但是考虑光伏温室改造与运行维护花费,在不改变原有结构的情况下,农作物光照饱和点对应的光照分配比例可作为PVG整体收益最大点。

由于本研究在模型建立过程中忽略了PVG的部分参数影响,如产量模型未考虑光伏温室内环境可控的特点,导致计算结果可能有部分偏差,下一步将优化农作物产量模型,提高农作物产量预估精度,考察PVG光照最优分配比例,基于光伏最优分配比例,以光伏温室整体收益为目标,优化光伏组件排列形状,可作为光伏温室设计的研究方向之一。

参考文献References

1

陈旭根,陈颖,房新月,等.基于光伏发电的温室滴灌系统设计[J].吉首大学学报(自然科学版),2019,40(4):52-55.CHEN X G,CHEN Y,FANG X Y,et al.Design of greenhouse drip irrigation system based on photovoltaic power generation[J].Journal of Jishou University (natural sciences edition),2019,40(4):52-55(in Chinese with English abstract). [百度学术] 

2

张旭.太阳能供电的连栋温室环境监测装置的设计[D].哈尔滨:东北农业大学,2019.ZHANG X.Design of environmental monitoring device for multi-span greenhouse based on solar energy supply[D].Harbin:Northeast Agricultural University,2019(in Chinese with English abstract). [百度学术] 

3

单天禹,梁春英.光伏幕布日光温室大棚测控系统设计[J].黑龙江八一农垦大学学报,2019,31(1):84-89.SHAN T Y,LIANG C Y.Design of measurement and control system for solar greenhouse with photovoltaic curtain[J].Journal of Heilongjiang Bayi Agricultural University,2019,31(1):84-89(in Chinese with English abstract). [百度学术] 

4

叶林,李建设,张光弟,等.宁夏不同光伏材料连栋温室冬季室内环境变化规律[J].北方园艺,2017(17):92-98.YE L,LI J S,ZHANG G D,et al.Change of winter indoor environment of greenhouse of different photovoltaic materials in Ningxia[J].Northern horticulture,2017(17):92-98(in Chinese with English abstract). [百度学术] 

5

吴艳红,马文舒,胡永生,等.沙漠多晶硅光伏连栋温室适宜栽培的瓜类作物筛选研究[J].山东农业科学,2017,49(3):73-77.WU Y H,MA W S,HU Y S,et al.Study on screening of melon crops suitable for cultivation in desert polycrystalline silicon photovolatic multispan greenhouse[J].Shandong agricultural sciences,2017,49(3):73-77(in Chinese with English abstract). [百度学术] 

6

祁娟霞,曹丽华,李建设,等.宁夏不同光伏温室和大棚冬季内环境比较研究[J].浙江农业学报,2017,29(3):414-420.QI J X,CAO L H,LI J S,et al.Comparative study on the environment of different photovoltaic greenhouses and greenhouse in winter in Ningxia[J].Acta agriculturae Zhejiangensis,2017,29(3):414-420(in Chinese with English abstract). [百度学术] 

7

胡永生,王荣,吴艳红,等.沙漠多晶硅光伏连栋温室茄果类作物适宜栽培的位置筛选[J].黑龙江农业科学,2016(11):82-87.HU Y S,WANG R,WU Y H,et al.Screening the solanaceous vegetable crops suitable for cultivation position in desert polycrystalline silicon photovoltaic greenhouse[J].Heilongjiang agricultural sciences,2016(11):82-87(in Chinese with English abstract). [百度学术] 

8

周莹,王双喜,刘中华,等.基于ANSYS的日光温室复合相变保温墙体的模拟研究[J].太阳能学报,2020,41(4):113-122.ZHOU Y,WANG S X,LIU Z H,et al.Simulation study on composite phase change thermal insulation walls in solar greenhouse based on ANSYS[J].Acta energiae solaris sinica,2020,41(4):113-122(in Chinese with English abstract). [百度学术] 

9

王彪,陈教料,应建阳.基于Ecotect的光伏温室内光环境分析[J].江苏农业科学,2018,46(8):221-224.WANG B,CHEN J L,YING J Y.Analysis of luminous environment of photovoltaic greenhouse using Ecotect models[J].Jiangsu agricultural sciences,2018,46(8):221-224. [百度学术] 

10

郭腾腾,涂洁磊,盛钰清,等.基于CFD稳态模拟的光伏玻璃温室温度分布模拟计算[J].云南师范大学学报(自然科学版),2018,38(2):11-15.GUO T T,TU J L,SHENG Y Q,et al.Simulation calculation of photovoltaic glass greenhouse temperature distribution based on CFD steady state simulation[J].Journal of Yunnan Normal University(natural sciences edition),2018,38(2):11-15(in Chinese with English abstract). [百度学术] 

11

成珂,马晓瑶,孙琦琦.光伏温室大棚组件布置CFD模拟研究[J].太阳能学报,2021,42(8):159-165.CHENG K,MA X Y,SUN Q Q.CFD simulation study on module layout of photovoltaic greenhouse[J].Acta energiae solaris sinica,2021,42(8):159-165(in Chinese with English abstract). [百度学术] 

12

蔡彦文.基于CFD光热流场耦合分析的光伏温室研究与设计[D].杭州:浙江工业大学,2015.CAI Y W.Study and design of photovoltaic greenhouse based on the analysis of optic-thermal coupling flow field using CFD technique[D].Hangzhou:Zhejiang University of Technology,2015(in Chinese with English abstract). [百度学术] 

13

LI C S,WANG H Y,MIAO H,et al.The economic and social performance of integrated photovoltaic and agricultural greenhouses systems:case study in China[J].Applied energy,2017,190:204-212. [百度学术] 

14

COSSU M,COSSU A,DELIGIOS P A,et al.Assessment and comparison of the solar radiation distribution inside the main commercial photovoltaic greenhouse types in Europe[J].Renewable and sustainable energy reviews,2018,94:822-834. [百度学术] 

15

MARUCCI A,CAPPUCCINI A.Dynamic photovoltaic greenhouse:energy efficiency in clear sky conditions[J].Applied energy,2016,170:362-376. [百度学术] 

16

BAMBARA J,ATHIENITIS A K.Energy and economic analysis for the design of greenhouses with semi-transparent photovoltaic cladding[J].Renewable energy,2019,131:1274-1287. [百度学术] 

17

YANO A,COSSU M.Energy sustainable greenhouse crop cultivation using photovoltaic technologies[J].Renewable and sustainable energy reviews,2019,109:116-137. [百度学术] 

18

PÉREZ-ALONSO J,PÉREZ-GARCÍA M,PASAMONTES-ROMERA M,et al.Performance analysis and neural modelling of a greenhouse integrated photovoltaic system[J].Renewable and sustainable energy reviews,2012,16(7):4675-4685. [百度学术] 

19

EZZAERI K,FATNASSI H,WIFAYA A,et al.Performance of photovoltaic Canarian greenhouse:a comparison study between summer and winter seasons[J].Solar energy,2020,198:275-282. [百度学术] 

20

王彪.光伏温室光热环境数值建模及性能分析[D].杭州:浙江工业大学,2017.WANG B.Study of solar thermal numerical model and performance analysis for the photovoltaic greenhouse[D].Hangzhou:Zhejiang University of Technology,2017(in Chinese with English abstract). [百度学术] 

21

叶子飘,李进省.光合作用对光响应的直角双曲线修正模型和非直角双曲线模型的对比研究[J].井冈山大学学报(自然科学版),2010,31(3):38-44.YE Z P,LI J S.Comparative investigation light response of photosynthesis on non-rectangular hyperbola model and modified model of rectangular hyperbola[J].Journal of Jinggangshan University (natural sciences edition),2010,31(3):38-44(in Chinese with English abstract). [百度学术]