2. 重庆大学 输配电装备及系统安全与新技术国家重点实验室, 重庆 400044;
3. 中国建筑西南设计研究院有限公司, 成都 610041
2. State Key Laboratory of Power Transmission Equipment&System Security and New Technology, Chongqing University, Chongqing 400044, P. R. China;
3. China Southwest Architectual Design and Resarch Instiute Co. Ltd., Chengdu 640041, P. R. China
微电网是近年来高速发展的综合性电力系统,可促进电网对可再生能源的消纳,将提供给网侧的能源多类型化。依据不同的运行状态,微电网存在并网和孤岛这两种不同的模式[1-2],当微网的并网条件不匹配时,将由并网过渡到孤岛状态[3]。
储能技术在微电网消纳可再生能源过程中有着重要作用,可破解目前电力“即发即用”的同步性问题,同时也可以对新能源的输出功率波动进行补偿[4-5]。
蓄电池超级电容器混合储能既能够充分应用超级电容功率型储能器件的物理特性、又可优化蓄电池能量型储能器件的充放电过程,延长蓄电池使用寿命,得到了业界的普遍关注[6]。
混合储能系统控制策略的核心在于平抑母线的功率波动,并且在功率型储能器件和能量型储能器件间进行合理的功率分配。文献[7]采用自适应小波包分解和模糊优化控制对母线的电压波动进行二次功率分配,优化储能系统的荷电状态变化过程。文献[8]提出一种基于移动平均滤波算法的自适应能量控制策略,有效解决了脉动负荷功率在蓄电池和超级电容间的功率分配问题。文献[9]利用高通、带通和低通滤波器等3个隐形滤波器分离不同频的波动功率,达到控制效果。上述文献的控制策略并未考虑超级电容的容量大小,在平抑功率波动时若超级电容的剩余容量超过了阈值,会导致系统不能达到控制目标。在此基础上,文献[10]提出了一种主从结构式自适应控制策略,并利用功率前馈解除了母线电压与扰动输入间的耦合关系,在延长蓄电池使用寿命的同时也将超级电容器的端电压有效维持在一个范围。但是该控制策略并不能抑制耦合扰动输入对超级电容端电压的影响,使得控制效果大打折扣。基于此,我们分析了互补PWM(pulse width modulation)导通控制方法在混合储能中的优势,并对其进行小信号建模,提出了一种双重解耦控制方案,通过仿真研究验证了所提控制策略的正确性。
1 系统模型建立 1.1 系统状态结构如图 1[11]所示,微电网孤岛状态下包含电源、储能和负载3个部分,它们分别通过对应的控制器与母线相连。
双向DC/DC变换器在储能器件和母线间起着连接桥梁的作用,可降低母线的电压波动,控制储能器件的充放电过程,同时也可降低储能器件的容量配置[12]。在进行分析时,通常用可控电流源im等效替代储能系统的外部扰动输入[13],如图 2所示。
双向DC/DC变换器具有成本低和工作效率较高的特点[14],可降低系统的重量和体积。互补导通模式可使变换器的开关损耗降至零,达到软开关的效果,且响应速度更快[15]。在微电网中,瞬时冲击较多,频繁的功率流向切换可导致较大的开关损耗,因此采用互补PWM导通模式更为适合。
对图 2所示的电路进行小信号建模,可以得到下列状态空间表达式[11]:
$ \begin{array}{*{20}{c}} {\left[ {\begin{array}{*{20}{c}} {{{\dot {\hat i}}_{{\rm{bi}}}}}\\ {{{\dot {\hat i}}_{{\rm{sci}}}}}\\ {{{\dot {\hat v}}_{{\rm{bus}}}}} \end{array}} \right] = \left( {\begin{array}{*{20}{c}} { - \frac{{{R_{{\rm{lb}}}} + {R_{\rm{b}}}}}{{{L_{\rm{b}}}}}}&0&{ - \frac{{\left( {1 - {D_1}} \right)}}{{{L_{\rm{b}}}}}}\\ 0&{ - \frac{{{R_{{\rm{lsc}}}} + {R_{{\rm{sc}}}}}}{{{L_{{\rm{sc}}}}}}}&{ - \frac{{\left( {1 - {D_2}} \right)}}{{{L_{{\rm{sc}}}}}}}\\ {\frac{{1 - {D_1}}}{{{C_{\rm{o}}}}}}&{\frac{{1 - {D_2}}}{{{C_{\rm{o}}}}}}&0 \end{array}} \right)\left( {\begin{array}{*{20}{c}} {{{\hat i}_{{\rm{bi}}}}}\\ {{{\hat i}_{{\rm{sci}}}}}\\ {{{\hat v}_{{\rm{bus}}}}} \end{array}} \right) + \left( {\begin{array}{*{20}{c}} {\frac{{{V_{{\rm{bus }}}}}}{{{L_{\rm{b}}}}}}&0\\ 0&{\frac{{{V_{{\rm{bus }}}}}}{{{L_{{\rm{sc}}}}}}}\\ { - \frac{{{I_{\rm{b}}}}}{{{C_{\rm{o}}}}}}&{ - \frac{{{I_{{\rm{sc}}}}}}{{{C_{\rm{o}}}}}} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {{{\hat d}_1}}\\ {{{\hat d}_2}} \end{array}} \right) + }\\ {\left( {\begin{array}{*{20}{c}} {\frac{1}{{{L_{\rm{b}}}}}}&0&0\\ 0&{\frac{1}{{{L_{{\rm{sc}}}}}}}&0\\ 0&0&{ - \frac{1}{{{C_{\rm{o}}}}}} \end{array}} \right)\left( {\begin{array}{*{20}{c}} {{{\hat v}_{\rm{b}}}}\\ {{{\hat v}_{{\rm{sc}}}}}\\ {{{\hat i}_{\rm{m}}}} \end{array}} \right)。} \end{array} $ | (1) |
式中:Lb、Lsc分别为蓄电池和超级电容端变换器的电感;Rlb、Rlsc分别为Lb、Lsc的内阻;vbus为母线电压;ibi、isci分别为储能系统的输出电流;vb、vsc分别为2个储能端变换器的端电压;d1、d2分别为开关S2、S4的占空比;D1、D2分别为开关S2、S4的稳态占空比;
根据上述模型,可分别求出isci对d2和ibi对d1的传递函数为[11]:
$ {G_{{\rm{iscd}}2}}\left( s \right) = \frac{{{{\hat i}_{{\rm{sc}}}}}}{{{{\hat d}_2}}} = \frac{{\left( {s + \frac{{{R_{{\rm{lb}}}} + {R_{\rm{b}}}}}{{{L_{\rm{b}}}}}} \right)\left[ {\left( {\frac{{1 - {D_2}}}{{{L_{{\rm{sc}}}}}}} \right)\left( { - \frac{{{I_{{\rm{sc}}}}}}{{{C_{\rm{o}}}}}} \right) - \frac{{{V_{{\rm{bus}}}}}}{{{L_{{\rm{sc}}}}}}s} \right] + \frac{{{V_{{\rm{bus}}}}}}{{{L_{{\rm{sc}}}}}}\left( {\frac{{1 - {D_1}}}{{{L_{\rm{b}}}}}} \right)\left( { - \frac{{1 - {D_1}}}{{{C_{\rm{o}}}}}} \right)}}{{ - \left( {s + \frac{{{R_{{\rm{lb}}}} + {R_{\rm{b}}}}}{{{L_{\rm{b}}}}}} \right)\left[ {s\left( {s + \frac{{{R_{{\rm{lsc}}}} + {R_{{\rm{sc}}}}}}{{{L_{{\rm{sc}}}}}}} \right) + \frac{{{{\left( {1 - {D_2}} \right)}^2}}}{{{C_{\rm{o}}}{L_{{\rm{sc}}}}}}} \right] - \frac{{{{\left( {1 - {D_1}} \right)}^2}}}{{{L_{\rm{b}}}{C_{\rm{o}}}}}\left( {s + \frac{{{R_{{\rm{lsc}}}} + {R_{{\rm{sc}}}}}}{{{L_{{\rm{sc}}}}}}} \right)}}, $ | (2) |
$ {G_{{\rm{ibd}}1}}\left( s \right) = \frac{{{{\hat i}_{\rm{b}}}}}{{{{\hat d}_1}}} = \frac{{\left( {s + \frac{{{R_{{\rm{lsc}}}} + {R_{{\rm{sc}}}}}}{{{L_{{\rm{sc}}}}}}} \right)\left[ {\left( {\frac{{1 - {D_1}}}{{{L_{\rm{b}}}}}} \right)\left( { - \frac{{{I_{\rm{b}}}}}{{{C_{\rm{o}}}}}} \right) - \frac{{{V_{{\rm{bus}}}}}}{{{L_{\rm{b}}}}}s} \right] + \frac{{{V_{{\rm{bus}}}}}}{{{L_{\rm{b}}}}}\left( {\frac{{1 - {D_2}}}{{{L_{{\rm{sc}}}}}}} \right)\left( { - \frac{{1 - {D_2}}}{{{C_{\rm{o}}}}}} \right)}}{{ - \left( {s + \frac{{{R_{{\rm{lb}}}} + {R_{\rm{b}}}}}{{{L_{\rm{b}}}}}} \right)\left[ {s\left( {s + \frac{{{R_{{\rm{lsc}}}} + {R_{{\rm{sc}}}}}}{{{L_{{\rm{sc}}}}}}} \right) + \frac{{{{\left( {1 - {D_2}} \right)}^2}}}{{{C_{\rm{o}}}{L_{{\rm{sc}}}}}}} \right] - \frac{{{{\left( {1 - {D_1}} \right)}^2}}}{{{L_{\rm{b}}}{C_{\rm{o}}}}}\left( {s + \frac{{{R_{{\rm{lsc}}}} + {R_{\rm{s}}}}}{{{L_{{\rm{sc}}}}}}} \right)}}。$ | (3) |
式中:Gibd1、Giscd2分别为蓄电池电流对开关S2的传递函数、超级电容电流对开关S4的传递函数;Vbus为母线电压稳态值。
2 混合储能双重解耦控制策略在孤岛状态下母线电压稳定、蓄电池充放电电流平滑及超级电容的能量密度小是该系统控制策略所要达到的目标。
2.1 主从结构双重解耦控制策略鉴于蓄电池与超级电容的性能互补特性,我们提出一种双重解耦控制策略,如图 3所示。
超级电容的控制环路维持母线电压的稳定,在电流内环中引入扰动前馈量,将扰动变化对系统的影响考虑在控制内并予以反馈。考虑到超级电容的容量,所提控制策略利用蓄电池的控制环路将超级电容的端电压维持在一定的设定值附近,使超级电容不会过充过放,同时也将扰动变化对其影响考虑在控制内环并用前馈控制予以反馈。
2.2 前馈项解耦系数推导双重解耦控制策略中前馈环解除了固有输入与扰动输入间的耦合关系,提高两个控制环路的独立性,同时也消除了扰动对控制目标的影响,使得系统对扰动干扰的抵抗力进一步提高。
对所提策略进行动态分析,由图 4(a)可以推导出母线电压与参考电压输入和扰动电流输入的表达式:
$ {U_{{\rm{bus}}}} = {T_{{\rm{ref}}}}{U_{{\rm{bus}} - {\rm{ref}}}} + {Z_{\rm{o}}}\left( { - {i_{\rm{m}}} + {i_{{\rm{bo}}}}} \right), $ | (4) |
其中
$ {T_{{\rm{ref}}}} = \frac{{{G_{{\rm{c}} - {\rm{usc}}}}{G_{{\rm{c - isc}}}}{K_{{\rm{pwm}}}}{G_{{\rm{isc}} - {\rm{d}}2}}}}{\Delta }\frac{{{I_{{\rm{sco}}}}}}{{{I_{{\rm{sci}}}}}}\frac{1}{{s{C_{\rm{o}}}}}, $ | (5) |
$ {Z_{\rm{o}}} = \frac{{1 + {G_{{\rm{c}} - {\rm{isc}}}}{K_{{\rm{pwm}}}}{G_{{\rm{isc}} - {\rm{d}}2}}{H_1} + {G_{{\rm{c}} - {\rm{isc}}}}{K_{{\rm{pwm}}}}{G_{{\rm{isc}} - {\rm{d}}2}}\frac{{{I_{{\rm{sco}}}}}}{{{I_{{\rm{sci}}}}}}{G_{\rm{x}}}}}{\Delta }\frac{1}{{s{C_{\rm{o}}}}}, $ | (6) |
$ \Delta = 1 + {G_{{\rm{c}} - {\rm{isc}}}}{K_{{\rm{pwm}}}}{G_{{\rm{isc}} - {\rm{d}}2}}{H_1} + {G_{{\rm{c}} - {\rm{usc}}}}{G_{{\rm{c}} - {\rm{isc}}}}{K_{{\rm{pwm}}}}{G_{{\rm{isc}} - {\rm{d}}2}}\frac{{{I_{{\rm{sco}}}}}}{{{I_{{\rm{sci}}}}}}{H_2}。$ | (7) |
式中:Tref为电压系数,Zo为扰动阻抗,Ubus-ref为母线电压参考值,Gc-usc为超级电容端电压外环的控制器传递函数,Gc-isc为超级电容端电流内环的控制器传递函数,Kpwm为比例系数,Isco为超级电容端稳态输出电流,Isci为超级电容端稳态输入电流,H1为超级电容端控制内环反馈通道系数,H2为超级电容端控制外环反馈通道系数。
由式(4)可知,由于存在扰动阻抗Zo,扰动电流将导致母线电压波动。若存在前馈项系数Gx使得Zo在理论上达到0,可使扰动电流对于母线电压的影响消失。令Zo=0,可以得到:
$ {G_{\rm{x}}} = - \frac{{1 + {G_{{\rm{c}} - {\rm{isc}}}}{K_{{\rm{pwm}}}}{G_{{\rm{isc}} - {\rm{d}}2}}{H_1}}}{{{G_{{\rm{c}} - {\rm{isc}}}}{K_{{\rm{pwm}}}}{G_{{\rm{isc}} - {\rm{d}}2}}}}\frac{{{I_{{\rm{sci}}}}}}{{{I_{{\rm{sco}}}}}}, $ | (8) |
$ {U_{{\rm{bus}}}} = {T_{{\rm{ref}}}}{U_{{\rm{bus}} - {\rm{ref}}}}, $ | (9) |
即直流母线电压理论上完全不再受扰动电流突变及耦合电流的影响,暂态时母线电压仍保持稳定。
同理可以推导出蓄电池端控制环路的前馈项系数Gxb的表达式为:
$ {G_{{\rm{xb}}}} = \frac{{1 + {G_{{\rm{c}} - {\rm{ib}}}}{K_{{\rm{pwm}}}}{G_{{\rm{ib - d}}1}}{H_3}}}{{{G_{{\rm{c}} - {\rm{ib}}}}{K_{{\rm{pwm}}}}{G_{{\rm{ib - d}}1}}}}\frac{{{I_{{\rm{bi}}}}}}{{{I_{{\rm{bo}}}}}}, $ | (10) |
式中:Gc-ub为蓄电池端电压外环的控制器传递函数,Gc-ib为蓄电池端电流内环的控制器传递函数,Ibo为蓄电池端稳态输出电流,Ibi为蓄电池端稳态输入电流,H3为蓄电池端控制内环反馈通道系数,H4为蓄电池端控制外环反馈通道系数。
超级电容电压理论上完全不再受扰动电流突变及耦合电流的影响,暂态时超级电容电压仍保持稳定。
2.3 混合储能系统功率分配系统截止频率越高,响应速度就会越快,而开环和闭环截止频率存在正相关的联系,因此系统的反应速度可用开环截止频率来表征[16]。根据2.2节所述,在我们提出的双重解耦控制结构中存在4个主回路闭环,若区别设置4个环路的开环频率,可使得波动功率的高低频分量自动隔离。
将蓄电池和超级电容的外环控制器开环频率区别设置,二者对于波动功率的响应速度也会明显不同。当外部环境变化导致储能系统功率缺额出现波动时,由超级电容器优先响应以弥补缺额,在此过程中其端电压亦会发生改变。根据蓄电池控制环路的目标,在超级电容端电压波动后,蓄电池会缓慢反应,利用自身蓄能为超级电容进行充放电操作,将超级电容的剩余电量在动态平衡保持稳定。在这个动态过程中,超级电容响应快,补偿高频功率波动;而蓄电池响应缓慢,补偿低频功率波动,降低蓄电池端电流变化率,延长其使用寿命。在整个平抑功率波动的过程中,并不影响母线电压的稳定。
3 仿真研究采用Matlab/Simulink对所提控制策略进行仿真验证。系统框图如图 4所示,电路参数见表 1,与文献[10]保持一致,开关频率50 kHz,超级电容电流内环、电压外环的开环截止频率设置为5 kHz和500 Hz,蓄电池电流内环、电压外环的开环截止频率分别设置为5 kHz和5 Hz。负载扰动用阶跃突变的极端情况来模拟。
如图 5所示,负载在1.000 s发生突变,储能系统的充电电流阶跃性变大,当负载电流从0.4 A跳变为0.8 A时,蓄电池和超级电容的充电电流都处于增大状态,但后者的电流变化率明显比前者高,系统在0.025 s内即恢复到稳态,证明超级电容承担了系统的瞬时冲击功率,而蓄电池对长时功率波动缺额进行了弥补。双重解耦控制策略使得高频和低频波动功率在两种储能装置内得到有效分配并补偿,同时蓄电池的充电电流平稳增大,变化率较小。
图 6为负载在1.000 s发生突变时直流母线电压瞬态响应的波形图,当负载电流变化时,母线电压依然保持稳定,基本没有波动,电压波动率小于0.1%,双重解耦控制策略可有效抑制母线的功率波动。
图 7为负载在1.30 s和1.60 s相继发生突变时,超级电容和蓄电池电流瞬态响应的波形图。
在1.30 s时刻,负载电流从0.8 A跳变为-3.2 A时,蓄电池和超级电容的充电电流都有所减小,但超级电容承担了系统的瞬时冲击功率,电流变化率较为明显,而蓄电池电流变化较为缓慢。在1.60 s时刻,负载电流从-3.2 A跳变为-11.6 A时,储能系统由充电状态转变为放电状态,蓄电池和超级电容的电流也都由正转负,此时超级电容的电流变化率依然较高,蓄电池的充电电流缓慢变化,系统在0.05 s内即恢复到稳态。证明超级电容承担了系统的瞬时冲击功率,而蓄电池对长时功率波动缺额进行了弥补。
图 8为负载在1.30 s和1.60 s相继发生突变时直流母线电压瞬态响应的波形图。在负载变化时,系统由充电状态转为放电状态,母线电压始终维持在100 V,电压波动率很小,证明所提控制策略对母线功率波动存在很强的抑制效果。
针对微电网孤岛运行状态下对混合储能系统的要求,提出一种主从结构双重解耦控制策略,利用功率前馈解除了母线电压与扰动输入间的耦合关系,也抑制了耦合扰动输入对超级电容端电压的影响,将其端电压有效维持在一个范围内,解决了传统控制策略下超级电容的过充过放问题,使蓄电池的充放电电流变化平滑,降低其变化率,延长了使用寿命。最后,仿真研究验证了所提控制策略的正确性和有效性。
[1] |
邱麟, 许烈, 郑泽东, 等. 微电网运行模式平滑切换的控制策略[J]. 电工技术学报, 2014, 29(2): 171-176. QIU Lin, XU Lie, ZHENG Zedong, et al. Control method of microgrid seamless switching[J]. Transactions of China Electrotechnical Society, 2014, 29(2): 171-176. (in Chinese) DOI:10.3969/j.issn.1000-6753.2014.02.019 |
[2] |
毕晓辉, 孙韬, 毛瑞, 等. 光储微电网混合储能系统的控制策略及开关优化[J]. 重庆大学学报, 2016, 39(6): 11-18. BI Xiaohui, SUN Tao, MAO Rui, et al. Control strategy and protection switch optimization of hybrid energy storage system for PV-storage microgrid system[J]. Journal of Chongqing University, 2016, 39(6): 11-18. (in Chinese) |
[3] |
Lu M S, Chang C L, Lee W J, et al. Combining the wind power generation system with energy storage equipment[J]. IEEE Transactions on Industry Applications, 2009, 45(6): 2109-2115. DOI:10.1109/TIA.2009.2031937 |
[4] |
Cipcigan L M, Taylor P C. Investigation of the reverse power flow requirements of high penetrations of small-scale embedded generation[J]. IET Renewable Power Generation, 2007, 1(3): 160. DOI:10.1049/iet-rpg:20070011 |
[5] |
Yao L Z, Yang B, Cui H F, et al. Challenges and progresses of energy storage technology and its application in power systems[J]. Journal of Modern Power Systems and Clean Energy, 2016, 4(4): 519-528. DOI:10.1007/s40565-016-0248-x |
[6] |
Zhou H H, Bhattacharya T, Tran D, et al. Composite energy storage system involving battery and ultracapacitor with dynamic energy management in microgrid applications[J]. IEEE Transactions on Power Electronics, 2011, 26(3): 923-930. |
[7] |
吴杰, 丁明. 采用自适应小波包分解的混合储能平抑风电波动控制策略[J]. 电力系统自动化, 2017, 41(3): 7-12. WU Jie, DING Ming. Wind power fluctuation smoothing strategy of hybrid energy storage system using self-adaptive wavelet packet decomposition[J]. Automation of Electric Power Systems, 2017, 41(3): 7-12. (in Chinese) |
[8] |
田明杰, 吴俊勇, 郝亮亮, 等. 基于多端口DC/DC变换器的混合储能系统自适应能量控制策略[J]. 电网技术, 2015, 39(12): 3378-3385. TIAN Mingjie, WU Junyong, HAO Liangliang, et al. Adaptive energy control strategy for hybrid energy storage system based on multiport DC/DC converters[J]. Power System Technology, 2015, 39(12): 3378-3385. (in Chinese) |
[9] |
李武华, 徐驰, 禹红斌, 等. 直流微网系统中混合储能分频协调控制策略[J]. 电工技术学报, 2016, 31(14): 84-92. LI Wuhua, XU Chi, YU Hongbin, et al. Frequency dividing coordinated control strategy for hybrid energy storage system of DC micro-grid[J]. Transactions of China Electrotechnical Society, 2016, 31(14): 84-92. (in Chinese) DOI:10.3969/j.issn.1000-6753.2016.14.010 |
[10] |
侯世英, 余海威, 李琦, 等. 微电网孤岛运行混合储能自适应控制策略[J]. 电力系统自动化, 2017, 41(17): 15-21. HOU Shiying, YU Haiwei, LI Qi, et al. Adaptive control strategy of hybrid energy storage in microgrid islanded operation state[J]. Automation of Electric Power Systems, 2017, 41(17): 15-21. (in Chinese) DOI:10.7500/AEPS20161213003 |
[11] |
余海威.微电网孤岛运行混合储能系统控制策略研究[D].重庆: 重庆大学, 2018. YU Haiwei. Research on control strategy of hybrid energy storage system in microgrid islanded operation state[D]. Chongqing: Chongqing University, 2018. (in Chinese) |
[12] |
桑丙玉, 陶以彬, 郑高, 等. 超级电容_蓄电池混合储能拓扑结构和控制策略研究[J]. 电力系统保护与控制, 2014, 42(2): 1-6. SANG Bingyu, TAO Yibin, ZHENG Gao, et al. Research on topology and control strategy of the super-capacitor and battery hybrid energy storage[J]. Power System Protection and Control, 2014, 42(2): 1-6. (in Chinese) |
[13] |
Kuperman A, Aharon I, Malki S, et al. Design of a semiactive battery-ultracapacitor hybrid energy source[J]. IEEE Transactions on Power Electronics, 2013, 28(2): 806-815. |
[14] |
张方华, 朱成花, 严仰光, 等. 双向DC-DC变换器的控制模型[J]. 中国电机工程学报, 2005, 25(11): 46-49. ZHANG Fanghua, ZHU Chenghua, YAN Yangguang, et al. The controlled model of bi-directional DC-DC converter[J]. Proceedings of the CSEE, 2005, 25(11): 46-49. (in Chinese) DOI:10.3321/j.issn:0258-8013.2005.11.009 |
[15] |
张国驹, 唐西胜, 周龙, 等. 基于互补PWM控制的Buck/Boost双向变换器在超级电容器储能中的应用[J]. 中国电机工程学报, 2011, 31(6): 15-21. ZHANG Guoju, TANG Xisheng, ZHOU Long, et al. Research on complementary PWM controlled buck/boost bi-directional converter in supercapacitor energy storage[J]. Proceedings of the CSEE, 2011, 31(6): 15-21. (in Chinese) |
[16] |
胡寿松. 自动控制原理[M]. 5版. 北京: 科学出版社, 2007: 226. HU Shousong. Principles of automatic control[M]. 5th ed. Beijing: Science Press, 2007: 226. (in Chinese) |