-
光谱吸收理论表明,不同气体由于其分子结构、能级的不同,对不同的光波吸收具有一定的选择吸收性——即“指纹”特性,利于此特性,可以对气体进行质和量的探测分析。
DIAL系统在工作时,向大气发射两束波长间隔很近的激光,一束光的波长位于甲醛的光谱吸收尖峰上,能够被气体强烈地吸收,称之为
${\lambda _{\rm on}}$ 光;另一束光的波长位于甲醛吸收峰的外侧,基本不被其吸收或吸收较弱,称之为${\lambda _{\rm off}}$ 。根据激光雷达原理,两束波长分别为${\lambda _{\rm on}}$ 和${\lambda _{\rm off}}$ 、发射激光峰值功率为P0、脉冲宽度为τ的光脉冲在大气中传输时,在距离R处的后向散射信号${P}_{(\rm on,}{}_{\rm off)}\left(R\right)$ 可表示为:$$ \begin{split} {P_{\left( {\rm on,off} \right)}}\left( R \right) =& {P_0}\eta \left( {\frac{A}{{{R^2}}}} \right)\left( {\frac{{c\tau }}{2}} \right){\beta _{\left( {\rm on,off} \right)}}\left( R \right)\exp \cdot\\ &\left[ { - 2\int_0^R {[{\alpha _{\left( {\rm on,off} \right)}}\left( r \right) + N(r){\sigma _{\left( {\rm on,off} \right)}}]{\rm{d}}r} } \right] \end{split}$$ (1) 式中:c为光速;A为望远镜有效接收面积;η为接收系统效率;β(R)为后向散射系数;α(r)为由大气散射导致的消光系数,由气体吸收导致的消光系数可表示为气体浓度N(r)与吸收截面σ的乘积。由于
${\lambda _{\rm on}}$ 和${\lambda _{\rm off}}$ 十分接近,可以认为大气中其他气体以及发射和接收系统对两束激光的散射影响基本一致,其后向散射系数β(λ, R)和散射导致的消光系数α(λ, r)在${\lambda _{\rm on}}$ 和${\lambda _{\rm off}}$ 处也是相等的,由${\lambda _{\rm on}}$ 、${\lambda _{\rm off}}$ 的后向散射光信号相比,可得到气体的浓度为:$$ N(R)\text{=}-\frac{1}{2(\sigma \left({\lambda }_{\rm on}\right)-\sigma \left({\lambda }_{\rm off}\right))}\frac{d}{{\rm{d}}R}\mathrm{ln} \left(\frac{{P}_{\rm on}\left(R\right)}{{P}_{\rm off}\left(R\right)}\right) $$ (2) 上式表明:通过测量不同高度处
${\lambda _{\rm on}}$ 和${\lambda _{\rm off}}$ 光波的后向散射光强度,已知被探测气体在这两个波段上的吸收截面情况下,就可以获得被测气体的浓度廓线。需要说明的是,气体分子的差分吸收截面会对浓度探测产生一定的影响,定义吸收截面差分:$\Delta \sigma {\text{ = }}\sigma \left( {{\lambda _{\rm on}}} \right) - \sigma \left( {{\lambda _{\rm off}}} \right)$ ,当$ \Delta \sigma $ 较大时,若气体浓度较大,则吸收加剧,会导致探测距离受限;当$ \Delta \sigma $ 较小时,若气体浓度较小,则吸收减弱,会导致探测灵敏度较低。因此,选择合理的探测波长、获得合适的$ \Delta \sigma $ 对于开放光程下甲醛浓度探测十分重要。根据DIAL原理,设计的甲醛浓度探测系统结构如图1所示:主要由激光发射系统、接收系统、数据采集以及处理系统组成。其中,两个量子级联激光器依次发出红外激光脉冲,经分色镜汇合后再与可见光波段的准直激光耦合,经准直、扩束后射入大气,在传输过程中与被探测甲醛、大气等相互作用,产生的回波信号经光学系统接收,由探测器进行光电转换,经滤波、反演处理后,得到浓度廓线。对于差分吸收激光雷达探测系统,通过控制量子级联激光器的注入电流及工作温度,将其输出的探测光分别调谐到被探测甲醛的特征吸收线。
${\lambda _{\rm on}}$ 和${\lambda _{\rm off}}$ 能否精确锁定甲醛的吸收线尖峰、吸收谷将直接影响探测性能,故为其配置锁频控制模块。此外,为了监控激光器的输出光能,将其一小部分通过探测器反馈回激光控制器,实现稳定的光能输出。 -
HITRAN数据库是一个汇编了众多气体分子吸收光谱学参数的数据库,根据HITRAN2012数据库,甲醛的强吸收谱线主要位于中、远红外波段。由于气体分子的复杂性,其吸收谱线包含多个吸收谱带,在选择探测波长时,除了考虑激光器探测波长的可获得性,还要考虑干扰气体、大气窗口特性等因素。相比于远红外波段,甲醛在中红外波段有多个吸收峰谷,因此,选择2600~3100 cm−1之间波段作为探测激光的初始筛选波段。
甲醛为痕量气体,其在大气中的含量为ppt(10−12)量级,即使在重污染区,其浓度一般也仅为ppm量级(体积浓度为1 ppm的甲醛对应的质量浓度为1.34 mg/m3,以某市为例,2019年环保监测站点醛类浓度为44.7 μg/m3,相当于33 ppt,而化工重污染区其浓度为其百倍左右)[18]。而标准大气中N2、O2、CO2、CH4以及H2O的浓度相对于甲醛高出很多,因此,在选择探测波长时,还需进一步考虑干扰气体的影响。根据后续探测需求,文中先设定探测甲醛的浓度为0.1~5 ppm,探测光程为1500 m,温度为300 K,压强为1 atm(1 atm=1.013×105),根据大气中上述干扰气体的平均浓度,进一步选择甲醛的探测波长。在选择探测波长时要充分考虑气体吸收谱线的展宽特性并依据下列原则进行选择:合适的
${\lambda _{\rm on}}$ 的吸收截面,以满足激光雷达最大探测范围的要求;足够大的${\lambda _{\rm on}}$ 与${\lambda _{\rm off}}$ 差分截面,以满足距离分辨率的要求。吸光度是描述气体对不同波长光束的吸收程度。据朗伯-比尔定律,当入射光强为
$ {I_0}\left( \lambda \right) $ 、波长为$ \lambda $ 的光束通过光程为$ L $ 、浓度为$ c $ 的气室后,若出射光强为$ I\left( \lambda \right) $ ,则吸光度$ \alpha \left( \lambda \right) $ 可表示为:$$ \alpha \left( \lambda \right) = {\text{ - ln}}\frac{{I\left( \lambda \right)}}{{{I_0}\left( \lambda \right)}} = \sigma \left( \lambda \right)cL $$ (3) 式中:
$ \sigma \left( \lambda \right) $ 为气体在波长$ \lambda $ 处的吸收截面,单位为cm2/molecule,该量是一个与吸收线强、吸收线型(与温度、压力相关)相关的量。图2所示为求解得到的标准大气中不同组分气体的吸光度曲线。计算时,
$ \sigma \left( \lambda \right) $ 可查阅HITRAN2012数据库,大气组分气体的浓度参照美国标准大气的浓度参数(N2浓度:77.026%;O2浓度:20.0%等),甲醛浓度人为设定为1 ppm;光程设定为1500 m,吸收线型选择Voigt线型函数(温度为300 K,压强为1 atm)。根据波长选择原则,文中选择探测的
${\lambda _{\rm on}}$ 和${\lambda _{\rm off}}$ 分别为2 778.48 cm−1和2 777.82 cm−1。同时,从图2可以看出,在大气中H2O和CH4对于甲醛的影响较大,因此,文中在研究激光雷达系统对于甲醛的探测能力时,选择H2O和CH4作为主要干扰气体。查阅HITRAN2012数据库,可得到上述三种气体在温度为300 K、压强为1 atm时在选定波长处的吸收截面值如表1所示。表 1 甲醛、水汽和甲烷气体的吸收截面
Table 1. Absorption cross section of formaldehyde, vapor and methane
Absorption cross section/cm2·molecule−1 Formaldehyde Vapor Methane ${\sigma _{{\rm{on}}} }$ ${\text{2} }{\text{.907 190 05} } \times {10^{ - 19} }$ ${\text{1} }{\text{.401 362 6} } \times {10^{ - 26} }$ ${\text{4} }{\text{.958 401 5} } \times {10^{ - 22} }$ ${\sigma _{\rm off} }$ ${\text{1} }{\text{.551 551 58} } \times {10^{ - 19} }$ ${\text{9} }{\text{.131 776 4} } \times {10^{ - 27} }$ ${\text{7} }{\text{.142 580 7} } \times {10^{ - 22} }$ -
激光雷达回波信号是进行气体浓度反演的基础,其功率的大小决定了系统的探测距离,根据选择的探测波长,可对DIAL的探测性能进行研究。研究时, DIAL的系统参数如表2所示;此外,根据标准大气组份比例,设定H2O的浓度为1.86%,CH4的浓度为1.7 ppm。
表 2 差分吸收激光雷达系统参数
Table 2. System parameters of DIAL
Parameter Value Laser wavelength/cm−1 λon : 2 778.48 λoff : 2 777.51 Pulse energy/mJ 1 Detector bandwidth/MHz 5 Detect active area/mm2 1×1 Pulse width/ns 25 Pulse repetition frequency/Hz 250 Laser divergence angle/mrad <1 Detector type VIGO-PVI-4 TE-4 Detector sensitivity/A·W−1 ≥1 Receiving area diameter 300 Telescope type MEADE-LX600-ACF12 Occlusion ratio 16.86% 由于被探测甲醛的浓度很低,导致后向散射光由吸收引起的光强变化非常微弱,能否实现有效探测的关键是探测器的探测灵敏度,研究选择波兰VIGO公司的HgCdTe型探测器,其参数如表3所示。
表 3 VIGO-PVI-4 TE-4探测器性能参数
Table 3. Performance parameters of VIGO-PVI-4 TE-4
Parameter Detector:VIGO-PVI-4 TE-4 Material HgCdTe Optium wavelength ${\lambda _{\rm opt} }$/$\text µ{\text{m} }$ 2.2-4.2 Detectivity $ {D^*} $/$\rm cm \cdot H{z^{1/2} }\cdot W^{-1}$ $ \geqslant 6.0 \times {10^{11}} $ Current responsivity Ri/A·W−1 $ \geqslant 1.0 $ Time constant t/ns $ \leqslant 100 $ Resistance-optical area product $ R \cdot A $/$\Omega \cdot \rm c{m^{-2}}$ $ \geqslant 800 $ Optical area A/$\rm mm^2$ 1 根据光电探测理论,结合以上参数,由归一化探测率
$ {D^*} $ 可得探测率$ D $ :$$ D = \frac{{{D^*}}}{{{{({A_{\rm d}}\Delta f)}^{1/2}}}} = \frac{{6.0 \times {{10}^{11}}}}{{22.360\;6}} = 2.683 \times {10^{10}}\;{{\rm W}^{ - 1}} $$ (4) 而噪声等效功率
$\rm NEP$ 为:$$ {\rm NEP} = \frac{1}{D} = 3.727 \times {10^{ - 11}}\;{\rm W} $$ (5) NEP是描述探测器探测微弱信号的能力,大于此值,表示信号能被探测器探测到。
-
回波信号的功率主要由雷达方程的消光项和雷达硬件参数决定,消光项包括甲醛吸收引起的消光和大气气溶胶、气体分子的消光。大气气溶胶、其他气体分子的消光可以根据美国标准大气模型获得。考虑到实际应用,重化工污染区域的污染气体的分布高度主要在1 km以内,所以探测区域选择在1.5 km以下的区域进行分析;为了模拟重化工污染区甲醛的分布,设定在整个探测路径有浓度为0.5 ppm的甲醛分布,另设定在0.6~0.8 km距离处,有一个最大浓度为1 ppm、浓度变化曲线为半波正弦的高浓度甲醛分布带。
图3所示为依据上述参数模拟获得的激光雷达回波信号功率随高度的分布图。由图可以看出,由于甲醛的吸收,在整个探测区域,
${\lambda _{\rm on}}$ 的回波功率衰减比${\lambda _{\rm off}}$ 大;在0.6~0.8 km处的高浓度甲醛分部区域,${\lambda _{\rm on}}$ 出现了较明显的下降;图中,紫色的点划线为探测器的噪声等效功率,对应的${\lambda _{\rm on}}$ 距离为1.1 km左右,这说明系统可以实现1 km以内的浓度探测。系统信噪比SNR是反映系统探测性能的一项关键指标,其值越大,表示信号的质量越佳。SNR的计算有模拟法和光子计数法,文中采用的是光子计数法[7],图4所示为求解得到的SNR随高度的变化图,从图中可以看出,随着系统探测距离的增加,SNR呈现出逐步减小的趋势,但在1.20 km以下,SNR大于100,表示系统可以探测幅值变化1%的信号,说明探测系统的SNR较高,能实现微弱的甲醛浓度的探测[19]。
通常,甲醛在重化工污染区域的分布是极其不规律、不稳定的,文中设定了不同高度、不同浓度的甲醛分布,并依据表1~3中探测系统的参数进行了模拟探测,结果如图5所示。从图中可以看出,由于甲醛气体的吸收,导致回波信号的信噪比降低,但在1 km范围内,反演结果与设定值有较好的一致性。
-
重化工污染区域污染气体无组织排放的不确定性会造成目标气体浓度发生变化,为此,需要研究目标气体浓度变化对系统探测能力的影响。由图6(a)可以看出,当甲醛在整个探测光程内的平均浓度为3 ppm时,由于吸收强烈,其最大探测距离仅为0.35 km左右,而当平均浓度小于1 ppm时,其探测距离为1.1 km;图6(b)所示为不同甲醛浓度对应的最远探测距离的拟合变化曲线。考虑到激光雷达探测的盲区影响,选取0.4 km以上的区间,这时,探测区域的甲醛平均浓度小于2.02 ppm。
图 6 不同甲醛浓度对回波信号功率及探测距离的影响
Figure 6. Effects of different formaldehyde concentrations on echo signal power and detection distance
图7所示为不同甲醛浓度对应的反演误差,其中图7(a)为不同甲醛浓度对应的反演误差,从图中可以看出,在该研究范围内,随着甲醛浓度的增加,模拟反演的浓度误差逐渐减小;图7(b)为不同甲醛浓度对应的平均浓度反演误差,从图中可以看出,浓度反演误差小于10%的部分为大于0.017 ppm。
-
在开放光程下的进行甲醛浓度探测时,干扰气体会对探测结果造成一定的影响,选择探测波段内影响最大的两种气体——H2O和CH4,作为干扰气体,研究其浓度变化对甲醛探测的影响。图8所示是以大气中H2O的平均浓度1.86%为单位浓度,H2O浓度分别为标准浓度的0.25、0.5、1、2、6、8、10倍时,
${\lambda _{\rm on}}$ 的回波信号功率。图8(a)中,随着H2O浓度的升高,激光雷达回波信号功率并没有产生明显的变化,这是因为在文中选择的探测波长处,H2O的吸收截面相对于CH4气体要小7个数量级。图8(b)为不同H2O浓度对应的反演误差的拟合曲线,由图可以看出,H2O浓度在在标准H2O浓度的0.25~10倍范围内变化时,引起的最大反演误差小于1.12%。同理,以大气中CH4的平均浓度1.7 ppm为单位浓度,将CH4的浓度设定为标准浓度的0.25、0.5、1、2、6、8、10倍时,得到不同浓度的CH4对
${\lambda _{\rm on}}$ 回波信号功率的影响,如图9所示。从图9(a)可以看出,随着CH4浓度的升高,激光雷达回波信号功率变化比H2O浓度变化引起的变化要大,这是因为与H2O相比,CH4的吸收截面相对于待测气体甲醛进小3个数量级。图9(b)为不同CH4浓度对应的反演误差的拟合曲线,由图可以看出,CH4浓度在在标准CH4浓度的0.25~10倍范围内变化时,引起的最大反演误差小于5%。综合探测距离及探测极限,当前选择的系统参数能对浓度为0.017~1.5 ppm范围内、距离为0.45~1.1 km范围内的甲醛进行探测,H2O及CH4对甲醛探测引入的误差分别在1.12%和5%之内。
Wavelength selection and detection performance of DIAL for formaldehyde gas detection
-
摘要: 为了对重化工污染区甲醛气体进行探测,开展了差分吸收激光雷达甲醛气体浓度探测研究。基于差分吸收激光雷达原理,针对甲醛气体在中红外波段有较强的吸收波段,并考虑开放光程下大气干扰气体的影响,选择了系统的探测波长
${\lambda _{\rm on}}$ 、${\lambda _{\rm off}}$ ;依据差分吸收激光雷达浓度反演方法,结合系统参数,对设计的DIAL系统的探测距离、气体浓度及干扰气体的可能影响进行了仿真研究。结果表明,该系统有望对浓度为0.017~1.5 ppm (1 ppm=10−6),距离为0.4~1.1 km范围内的甲醛气体进行探测,相对误差小于5%,可满足对重化工污染区甲醛气体浓度探测需求。文中的研究可为应用于化工重污染区域甲醛气体无组织排放监测的中红外差分吸收激光雷达系统的研制提供理论依据和技术基础。Abstract: In order to detect formaldehyde gas in severe chemical pollution area, the research of concentration detection for formaldehyde by differential absorption lidar (DIAL) is carried out in this paper. Based on the principle of differential absorption lidar and the strong absorption of formaldehyde gas in the mid infrared band, the detection wavelengths of${\lambda _{\rm on}}$ and${\lambda _{\rm off}}$ for the formaldehyde are selected in consideration of the influence of atmospheric interference gas under open optical path. According to the concentration inversion method of differential absorption lidar and combined with the system parameters, the effects of the detection distance, gas concentration and interfering gas on the system perform are studied. The results show that the system can detect formaldehyde concentration in the range of 0.017-1.5 ppm (1 ppm=10−6) and within the distance of 0.4-1.1 km, and the relative error is less than 5%, which can meet its demand for formaldehyde detection in severe chemical pollution area. This paper can provide a theoretical and technical basis for the development of mid infrared differential absorption lidar system for monitoring the unorganized emission of formaldehyde in severe chemical pollution area.-
Key words:
- DIAL /
- formaldehyde gas /
- absorption cross section /
- concentration inversion
-
表 1 甲醛、水汽和甲烷气体的吸收截面
Table 1. Absorption cross section of formaldehyde, vapor and methane
Absorption cross section/cm2·molecule−1 Formaldehyde Vapor Methane ${\sigma _{{\rm{on}}} }$ ${\text{2} }{\text{.907 190 05} } \times {10^{ - 19} }$ ${\text{1} }{\text{.401 362 6} } \times {10^{ - 26} }$ ${\text{4} }{\text{.958 401 5} } \times {10^{ - 22} }$ ${\sigma _{\rm off} }$ ${\text{1} }{\text{.551 551 58} } \times {10^{ - 19} }$ ${\text{9} }{\text{.131 776 4} } \times {10^{ - 27} }$ ${\text{7} }{\text{.142 580 7} } \times {10^{ - 22} }$ 表 2 差分吸收激光雷达系统参数
Table 2. System parameters of DIAL
Parameter Value Laser wavelength/cm−1 λon : 2 778.48 λoff : 2 777.51 Pulse energy/mJ 1 Detector bandwidth/MHz 5 Detect active area/mm2 1×1 Pulse width/ns 25 Pulse repetition frequency/Hz 250 Laser divergence angle/mrad <1 Detector type VIGO-PVI-4 TE-4 Detector sensitivity/A·W−1 ≥1 Receiving area diameter 300 Telescope type MEADE-LX600-ACF12 Occlusion ratio 16.86% 表 3 VIGO-PVI-4 TE-4探测器性能参数
Table 3. Performance parameters of VIGO-PVI-4 TE-4
Parameter Detector:VIGO-PVI-4 TE-4 Material HgCdTe Optium wavelength ${\lambda _{\rm opt} }$ /$\text µ{\text{m} }$ 2.2-4.2 Detectivity $ {D^*} $ /$\rm cm \cdot H{z^{1/2} }\cdot W^{-1}$ $ \geqslant 6.0 \times {10^{11}} $ Current responsivity Ri/A·W−1 $ \geqslant 1.0 $ Time constant t/ns $ \leqslant 100 $ Resistance-optical area product $ R \cdot A $ /$\Omega \cdot \rm c{m^{-2}}$ $ \geqslant 800 $ Optical area A/ $\rm mm^2$ 1 -
[1] Golden R. Identifying an indoor air exposure limit for formaldehyde considering both irritation and cancer hazards [J]. Critical Reviews in Toxicology, 2011, 41(8): 671-721. [2] Liu Wenqing, Cheng Zhenyi, Liu Jianguo, et al. Stereoscopic monitoring technology and applications for the atmospheric environment in China [J]. Chin Sci Bull, 2016, 61(30): 3196. doi: 10.1360/N972016-00394 [3] Di Huige, Hua Dengxin. Research status and progress of Lidar for atmosphere in China [J]. Infrared and Laser Engineering, 2021, 50(3): 20210032. (in Chinese) doi: 10.3788/IRLA20210032 [4] Chen Yuqin, Tan Wenyuan, Fu Dayou, et al. Research progress of formaldehyde detection [J]. Applied Chemical Industry, 2018, 47(6): 1258-1262. [5] Xu Fei, Zhou Xiaobin, Liu Zhengbo, et al. Near-infrared optical-feedback linear cavity-enhanced absorption spectroscopy [J]. Optics and Precision Engineering, 2021, 29(5): 933-939. (in Chinese) doi: 10.37188/OPE.2020.0657 [6] Liu Qiuwu, Wang Xiaobin, Chen Yafeng, et al. Detection of atmospheric NO2 concentration by differential absorption lidar based on dye laser [J]. Acta Optica Sinica, 2017, 37(4): 0428004. (in Chinese) [7] Shao Jiangfeng, Hua Dengxin, Wang Li, et al. Development of ultraviolet dual-wavelength lidar and analysis of its signal-to-noise ratio [J]. Acta Optica Sinica, 2020, 40(12): 1201004. (in Chinese) [8] Gong Yu, Bu Lingbing, Yang Bin, et al. High repetition rate mid-infrared differential absorption lidar for atmospheric pollution detection [J]. Sensors, 2020, 20(8): 2211. doi: 10.3390/s20082211 [9] Qi Zhong, Zhang Teng, Han Ge, et al. A nonlinear merging method of analog and photon signals for CO2 detection in lower altitudes using differential absorption lidar [J]. Optics Communications, 2017, 388: 68-76. doi: 10.1016/j.optcom.2016.11.071 [10] Veerabuthiran S, Razdan A K, Jindal M K, et al. Open field testing of mid IR DIAL for remote detection of thiodiglycol vapor plumes in the topographic target configuration [J]. Sens Actuator B-Chem, 2019, 298: 7. [11] Weibring P, Edner H, Svanberg S. Versatile mobile lidar system for environmental monitoring [J]. Appl Opt, 2003, 42(18): 3583-3594. doi: 10.1364/AO.42.003583 [12] Hong Guanglie, Li Jiatang, Kong Wei, et al. 935 nm differential absorption lidar system and water vapor profiles in convective boundary layer [J]. Acta Optica Sinica, 2017, 37(2): 0201003. (in Chinese) [13] Zhang Yang, Huang Weidong, Dong Changzhe, et al. Research on the development of the detection satellite technology in oceanographic lidar [J]. Infrared and Laser Engineering, 2020, 49(11): 20201045. (in Chinese) doi: 10.3788/IRLA20201045 [14] Edner H, Fredriksson K, Sunesson A, et al. Mobile remote sensing system for atmospheric monitoring [J]. Appl Opt, 1987, 26(19): 4330. doi: 10.1364/AO.26.004330 [15] Gholamreza Shayeganrad. Single laser-based differential absorption lidar (DIAL) for remote profiling atmospheric oxygen [J]. Optics and Lasers in Engineering, 2018, 111: 80-85. doi: 10.1016/j.optlaseng.2018.07.015 [16] Ge Ye, Shu Rong, Hu Yihua, et al. System design and performance simulation of ground-based differential absorption lidar for water-vapor measurements [J]. Acta Phys Sin, 2014, 63(20): 204301. (in Chinese) [17] Xu Ling, Bu Lingbing, Cai Haoze, et al. Wavelength selection and detection capability simulation of the midinfrared DIAL for NO2 detecion [J]. Infrared and Laser Engineering, 2018, 47(10): 1030002. (in Chinese) doi: 10.3788/IRLA201847.1030002 [18] Yao Yizhou, Zhang Xianbao. Analysis of pollution characteristics and sources of volatile organic compounds in Zhenjiang [J]. The Administration and Technique of Environ-mental Monitoring, 2021, 33(5): 60-63. [19] Mei Liang. Applications of laser spectroscopy in enviroumental monitoring[D]. Hangzhou: Zhejiang University, 2013. (in Chinese)