1 Introduction
Hydrazine and its derivatives, which are highly toxic and strong reducing agents, are miscible with water and soluble in most organic solvents. Hydrazine, methylhydrazine, 1,1-dimethylhydrazine, phenylhydrazine, 2-hydroxyethylhydrazine, and tert-butylhydrazine are the most common hydrazine derivatives. They are important chemical raw materials and widely used in polymers, metallic materials, metallurgy, printing and dyes, textiles, medicine, disease treatment, etc. [1-3]. In particular, they are used as propellants for altitude control engines in missiles, satellites, and the aerospace industry because they readily combust to produce hydrogen and nitrogen [4-6]. Owing to their strong reduction properties, they have been extensively studied in the Purex process as a support reducing stabilizer for the recovery of Pu(Ⅳ) using Fe(II), U(IV), hydroxylamine, and urea. They have also been considered as organic salt-free reductants for Pu(Ⅳ) and Np(Ⅵ) during the post-treatment of nuclear fuel [7]. Accordingly, significant research effort has been targeted at their application as support reducing agents or reducing agents for spent nuclear fuel reprocessing in recent years [8].
A significant amount of hydrazine derivatives are released into the environment as pollutants [9-12]. They are strong stimulants and may induce acute poisoning via injection, inhalation, skin poisoning, and gastrointestinal absorption or cause tumor growth. Thus, hydrazine derivatives are considered to be very harmful for living organisms and the environment. Once they are released to the environment, they can physically adsorb and interact chemically with the soil, resulting in a polluted environment that is difficult to reclaim. Some of the degradation products of hydrazine derivatives are more toxic that the derivatives themselves. Therefore, efficient wastewater disposal and decomposition treatment of effluent containing hydrazine derivatives is important.
The stoichiometry, kinetics, and mechanisms of the decomposition of hydrazine and its derivatives have yet to be fully elucidated [13]. Currently, the following types of treatment methods for hydrazine compounds are used globally: (1) Physical methods, such as incineration, ion exchange, adsorption, etc., are simple to use but readily produce secondary pollution; also, the regeneration period and lifetime of the ion resin are limited. (2) Biodegradation, which includes activated sludge (microbial) and aquatic plant treatment, is an environmentally friendly method; however, it cannot be applied to non-industrial effluent because the pollutants could bioaccumulate in organisms and avoid fundamental degradation [14]. (3) Chemical oxidation methods, such as catalytic oxidation, electrolyzed oxidizing water (EOW) oxidation, and photocatalytic oxidation, etc. [15], quickly oxidize and remove the organic matter by destroying the organic molecular structure. These methods are suitable for industrial processes; however, they are restricted to wastewater with low pollutant concentrations. (4) Other new processing technologies, such as low temperature plasma technology [16, 17], catalyzed supercritical water oxidation [18, 19], Fenton oxidation, and membrane bioreactor treatment [20], have also been applied. Fenton’s method has been widely studied because of the high activity of hydroxyl radicals towards the oxidation and decomposition of a variety of organic substances that are typically difficult to oxidize.
Radiation technology has been used in the environmental field [21], which is an important aspect of the peaceful use of atomic energy in the 21st century, as proposed by the International Atomic Energy Agency. Radiation technology is a non-polluting, inexpensive, simple, and efficient degradation method. The application of beta radiation technology to remove hydrazine pollutants should generate the target outcome with a lower energy consumption and at a faster rate. Hydrazine and its derivatives have been actively researched as support reducing stabilizers and novel salt free reductants; therefore, their stability under radiation exposure during the post-treatment process is important for the optimal selection of reducing agents for spent fuel reprocessing.
The quantitative structure-activity relationships (QSAR) method [22-24] is widely used during the development of new medicines and is a well-recognized research method. However, it is rarely used in other fields, especially that of radiochemistry. In this study, five hydrazine derivatives, i.e., hydrazine hydrate, methylhydrazine, dimethylhydrazine, tert-butylhydrazine, and 2-hydroxyethylhydrazine, were irradiated using a pulsed electron beam. The data were then analyzed to establish the QSAR model. This investigation of the radiative stability of hydrazine and its derivatives under β radiolysis expands the fundamental knowledge on the radiochemistry of inorganic compounds and provides abundant theoretical evidence of the degradation of hydrazine. It also provides an important guideline for the degradation of hydrazine and its derivatives in spent nuclear fuel.
2 Source and structure of the target compound
2.1 Chemicals and equipment
All chemicals were of analytical grade and used without further purification. Hydrazine hydrate (85%) was purchased from Shanghai Hushi Laboratorial Equipment Co. Ltd. (China). tert-Butylhydrazine (98%) and 2-hydroxyethylhydrazine (95%) were purchased from Shanghai Macklin Biochemical Co. Ltd. (China). Methylhydrazine, dimethylhydrazine, and p-dimethylaminobenzaldehyde were purchased from Sinopharm Chemical Reagent Co. Ltd. (China). Ethanol, nitric acid, and methane sulfonic acid (MSA) were purchased from Beijing Chemical Works (China). A beta particle accelerator (2 MeV), ion chromatograph (CIC-100), ultraviolet spectrophotomer (Perkin Elmer Lamda 45), and automatic potentiometric titrator (Mettler-Toledo G20 compact titrator) were used during the experiment. The beta particle accelerator was rented from the Department of Nuclear Technology Applications in China Institute of Atomic Energy, and the required safety protection was provided. The ion chromatograph was purchased from Qingdao Sheng Han Chromatography Technology Co. Ltd. (China); a CS12A column was used as the cation column, and the baseline noise was stabilized at around 20 mV before detection.
2.2 Experimental method
There are a few published research papers on the stability of hydrazine and its derivatives under radiation. For example, a group at the National Institute of Atomic Energy researched their stability under γ radiation. However, γ rays are difficult to use because of their high energy and radiation dose, which makes them unsuitable for large-scale industrial applications. There is no suitable kinetic data regarding the β radiation stability of hydrazine and its derivatives. Therefore, in this paper, source data were obtained using experimental and computational methods. The initial solutions of hydrazine and its derivatives were diluted to 0.05 mol/L; then, 250 mL aliquots of the solutions were prepared and irradiated using a pulsed electron beam from a beta particle accelerator. The accelerator used a 30 mA current, 50 Hz frequency, 4 μs pulse width, and 0.25 KGy/s dose. Then, 0.1 mL samples of the irradiated solutions were extracted using a syringe after the following irradiation durations: 0, 0.5, 1, 2, 3, 4, 5, 6, 7, 9, 11, 13, 15, 18, and 21 min.
2.3 Analytic procedures
Using methane sulfonic acid as the eluent, the concentrations of the hydrazine, tert-butylhydrazine, 1,1-dimethylhydrazine, and 2-hydroxyethylhydrazine solutions were separately determined using ion chromatography. The concentration of methylhydrazine was determined using an ultraviolet visible spectrophotometer with the color rendering of p-dimethylaminobenzaldehyde. The concentrations of the irradiated solutions were determined after the following irradiation durations: 0, 0.5, 1, 2, 3, 4, 5, 6, 7, 9, 11, 13, 15, 18, and 21 min. Then, the kinetic data for the decomposition of the hydrazine derivatives under beta radiation were obtained. The half-reaction times of hydrazine and its derivatives under β-radiation are expressed by τ50, which is the time required for the compound to decompose to 50% of the original concentration: A smaller τ50 value indicates a faster decomposition rate. The molecular formulas of various hydrazine derivatives and their τ50 values are listed in Table 1.
Derivatives | τ50 (s) |
---|---|
Hydrazine, N2H4 | 8155.6 |
methyl hydrazine, CH3N2H3 | 6689.2 |
1,1-dimethylhydrazine, (CH3)2N2H2 | 791.6 |
tert-butylhydrazine, (CH3)3CHN2H3 | 426.2 |
2-hydroxyethylhydrazine, HOCH2CH2N2H3 | 204.5 |
3 Calculation methods
The initial structures of the hydrazine compounds were generated using ChemOffice software, and geometric pre-optimization was conducted using built-in molecular mechanics MM2 modules in Chem3D Ultra. The minimum energy conformation of the hydrazine derivative was obtained after optimization. Geometric optimization and energy calculations were conducted on a variety of target compounds by applying the density functional theory B3LYP method of the Gaussian 03 program package with a 6-311+(3d, 3p) basis set. The vibration frequencies calculated for all of the molecular structures did not have imaginary frequencies, which indicates that the configuration obtained from the calculation was the most stable. Using the HyperChem software package, calculations were conducted on the optimized molecules and the physicochemical parameters were obtained; these included the total energy (Etotal), highest occupied molecular orbital energy (EHOMO), lowest unoccupied molecular orbital energy (ELUMO), molecular orbital transition energy (ΔE), molecular dipole moment (μ), hydrophobic parameter (logP), molecular refractivity (R), molecular volume (V), molecular superficial area (A), molecular molar mass (M), molecular polarizability (P), and hydration energy (HE). Mathematical statistical software, SPSS, was applied to perform a correlation analysis and stepwise regression analysis on the τ50 and physicochemical parameters of the selected hydrazine derivatives. Finally, a QSAR equation with sufficient correlation was obtained. The convergence precision of the quantum chemical calculations is a default value.
4 Results and discussion
4.1 3D stable molecular structures of the hydrazine compounds
Optimized geometric construction and frequency analysis of the hydrazine derivatives were conducted by adopting the density functional theory B3LYP method of the Gaussian 03 program package with the 6-311+(3d, 3p) basis set. No imaginary frequency occurred in the calculation results, indicating that they were reasonable and reliable. This also confirms that the most stable configuration was obtained. The optimized 3D structures are shown in Fig. 1.
-201805/1001-8042-29-05-005/alternativeImage/1001-8042-29-05-005-F001.jpg)
4.2 Values of the quantum chemical parameters
The following parameter values were acquired from the quantum chemical calculations: ELUMO, EHOMO, ΔE (ΔE = ELUMO − EHOMO), μ, and Etotal. In addition, the following parameters were obtained using the molecular simulation software, HyperChem: logP, R, V, molecular superficial area (A - approx, G - grid), M, P, and HE. The degradation rates of hydrazine and its derivatives under β irradiation are expressed by τ50. The quantum parameter values are shown in Table 2.
EHOMO | ELUMO | ΔE | Etotal | μ | logP | R | V | A | G | M | P | HE | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|
hydrazine | -0.218 | 0.011 | 0.229 | -111.912 | 0.000 | -0.680 | 10.570 | 183.870 | 135.890 | 158.850 | 32.050 | 3.480 | -21.860 |
methylhydrazine | -0.243 | -0.002 | 0.241 | -151.237 | 1.740 | -0.460 | 14.280 | 239.450 | 143.590 | 190.220 | 46.070 | 5.310 | -13.590 |
1,1-dimethylhydrazine | -0.203 | -0.003 | 0.201 | -190.554 | 0.564 | -0.550 | 19.180 | 291.360 | 134.660 | 217.930 | 60.100 | 7.150 | -.640 |
tert-butylhydrazine | -0.239 | -0.004 | 0.235 | -269.218 | 1.746 | 0.760 | 27.970 | 393.630 | 142.060 | 275.410 | 88.150 | 10.820 | -4.730 |
2-hydroxyethylhydrazine | -0.247 | -0.005 | 0.241 | -265.806 | 1.183 | -0.910 | 20.570 | 319.840 | 246.110 | 239.980 | 76.100 | 7.780 | -17.450 |
4.3 Structure-performance equation for the hydrazine derivatives under β irradiation
Before a QSAR regression equation was established, the quantum chemical parameters adopted were applied for preprocessing. The Pearson correlation coefficient analysis parameters that were significantly related to the β irradiation half-reaction time were obtained to facilitate elucidation of the relationship between the parameters of each structure and the degradation rate of the target compound.
Statistical software, SPSS, was applied to conduct a correlation analysis on the QSAR parameters to obtain the Pearson correlation coefficient matrix (Table 3).
τ50 | EHOMO | ELUMO | ΔE | Etotal | μ | logP | R | V | A | G | M | P | HE | |
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|
τ50 | 1 | 0.085 | 0.782 | 0.221 | 0.909 | -0.338 | -0.292 | -0.883 | -0.866 | -0.438 | -0.891 | -0.906 | -0.863 | -0.632 |
EHOMO | 0.085 | 1 | 0.448 | -0.933 | 0.462 | -0.762 | -0.177 | -0.339 | -0.288 | -0.564 | -0.369 | -0.401 | -0.310 | -0.350 |
ELUMO | 0.782 | 0.448 | 1 | -0.095 | 0.798 | -0.770 | -0.262 | -0.782 | -0.735 | -0.425 | -0.793 | -0.795 | -0.756 | -0.572 |
ΔE | 0.221 | -0.933 | -0.095 | 1 | -0.192 | 0.538 | 0.091 | 0.062 | 0.024 | 0.456 | 0.091 | 0.126 | 0.040 | -0.620 |
Etotal | 0.909 | 0.462 | 0.798 | -0.192 | 1 | -0.569 | -0.419 | -0.944 | -0.919 | -0.569 | -0.960 | -0.984 | -0.920 | -0.403 |
μ | -0.338 | -0.762 | -0.770 | 0.538 | -0.569 | 1 | 0.536 | 0.627 | 0.593 | 0.173 | 0.629 | 0.605 | 0.621 | 0.301 |
logP | -0.292 | -0.177 | -0.262 | 0.091 | -0.419 | 0.536 | 1 | 0.665 | 0.712 | -0.425 | 0.631 | 0.562 | 0.710 | 0.535 |
R | -0.883 | -0.339 | -0.782 | 0.062 | -0.944 | 0.627 | 0.665 | 1 | 0.996 | 0.266 | 0.998 | 0.988 | 0.998 | 0.620 |
V | -0.866 | -0.288 | -0.735 | 0.024 | -0.919 | 0.593 | 0.712 | 0.996 | 1 | 0.198 | 0.991 | 0.974 | 0.999 | 0.644 |
A | -0.438 | -0.564 | -0.425 | 0.456 | -0.569 | 0.173 | -0.425 | 0.266 | 0.198 | 1 | 0.319 | 0.412 | 0.203 | -0.369 |
G | -0.891 | -0.369 | -0.793 | 0.091 | -0.960 | 0.629 | 0.631 | 0.998 | 0.991 | 0.319 | 1 | 0.995 | 0.993 | 0.586 |
M | -0.906 | -0.401 | -0.795 | 0.126 | -0.984 | 0.605 | 0.562 | 0.988 | 0.974 | 0.412 | 0.995 | 1 | 0.976 | 0.523 |
P | -0.863 | -0.310 | -0.756 | 0.040 | -0.920 | 0.621 | 0.710 | 0.998 | 0.999 | 0.203 | 0.993 | 0.976 | 1 | 0.646 |
HE | -0.632 | 0.350 | -0.572 | -0.620 | -0.733 | -0.403 | 0.301 | 0.535 | 0.620 | 0.644 | -0.369 | 0.523 | 0.646 | 1 |
It is evident from the data in Table 3 that the correlations between τ50 and Etotal and M are significant. G, R, V, ELUMO, HE, A, μ, P, ΔE, and logP also had strong correlations with τ50. However, there was no significant correlation between EHOMO and τ50.
The Pearson correlation analysis shows that there is a close relationship between the total energy of the molecules and τ50 of the hydrazine compounds radiolysis reaction under beta irradiation; their relationship diagram was established, as depicted in Fig. 2.
-201805/1001-8042-29-05-005/alternativeImage/1001-8042-29-05-005-F002.jpg)
In Fig. 2, a relationship between the Etotal and τ50 of the hydrazine derivatives hydrolysis reaction under beta irradiation is slightly evident. In general, τ50 is positively correlated with Etotal; however, the correlation cannot be defined using a mathematical relationship. Accordingly, we used a multiple linear stepwise regression method to analyze the data.
4.4 QSAR equation for the hydrazine derivatives under β irradiation
The QSAR equation for the hydrazine derivatives (i.e., hydrazine, methylhydrazine, 1,1-dimethylhydrazine, tert-butylhydrazine, and 2-hydroxyethylhydrazin) under β irradiation was established using a multiple linear stepwise regression method, and the results are shown in Tables 4, 5, 6 and 7.
Model | Variables entered | Variables removed | method |
---|---|---|---|
1 | Molecule total energy | stepwise | |
2 | Orbital transition energy | stepwise |
Model | R | R square | Adjusted R square | Standard error of the estimate |
---|---|---|---|---|
1 | 0.909a | 0.827 | 0.769 | 1848.089752 |
2 | 0.995b | 0.989 | 0.979 | 559.650556 |
Model | Sum of squares | df | Mean square | F | Significance | |
---|---|---|---|---|---|---|
1 | Regression | 4.894×107 | 1 | 4.894×107 | 14.329 | 0.032a |
Residual | 1.025×107 | 3 | 3.415×106 | |||
Total | 5.918×107 | 4 | ||||
2 | Residual | 5.856×107 | 2 | 2.928×107 | 93.481 | 0.011b |
Residual | 6.264×105 | 2 | 3.132×105 | |||
Total | 5.918×107 | 4 |
Model | Unstandardized Coefficients | Standardized coefficients | t | Significance | ||
---|---|---|---|---|---|---|
B | standard error | Beta | ||||
1 | (Constant) | 13205.114 | 2755.880 | 4.792 | 0.017 | |
total energy | 50.326 | 13.295 | 0.909 | 3.785 | 0.032 | |
2 | (Constant) | -7583.464 | 3842.799 | -1.973 | 0.187 | |
Molecule total energy | 54.687 | 4.102 | 0.988 | 13.331 | 0.006 | |
Orbital transition energy | 94333.586 | 17021.516 | 0.411 | 5.542 | 0.031 |
The stepwise regression analysis results suggest the following: 1) Significant correlation exists between Etotal and τ50 as well as Etotal and ΔE, as shown in Tables 4 and 5. The square root of the determination coefficient, determination coefficient, and adjusted determination coefficient were 0.995, 0.989, and 0.979, respectively; these values indicate that the equations could account for 97.9% of the change in the τ50 and that the correlations between independent and dependent variable are evident. (2) The square of the equation, residual sum of the squares, total sum of the squares, F-statistics value, and significance were 5.858 × 107, 0.955, 6.264 × 105, 93.481, and 0.011< 0.05, respectively, indicating that a significant linear relationship exists between the dependent and independent variables. (3) From Table 7, it is evident that the QSAR equation for the hydrazine derivatives and τ50 of the hydrazine derivatives hydrolysis reaction under beta radiation was Y = −7583.464 + 54.687X1 + 94333.586X2, where X1 and X2 represent Etotal and ΔE, respectively. The significance levels of the regression coefficients are 0.006 and 0.031; these values are less than 0.05. This model fully represents the relationship between the hydrazine derivatives and β radiolysis stability.
4.5 The relationship between the hydrazine derivative structure and half-reaction time for β irradiation
In the structural activity equation, the τ50 for hydrazine compound radiolysis relates closely with Etotal and ΔE. A diagram of their relationship is shown in Fig. 3.
-201805/1001-8042-29-05-005/alternativeImage/1001-8042-29-05-005-F003.jpg)
The analysis results indicate the following: 1) The relationship between the structure of the hydrazine derivative and reaction time can be described quantitatively by the QSAR equation even for complex hydrazine derivative structures. 2) The decomposition rate under β irradiation and structure of the hydrazine derivatives are negatively correlated and can be described quantitatively by Etotal and ΔE. 3) The relationship between the decomposition rate under β irradiation and structure of the hydrazine derivatives can be described quantitatively by the QSAR equation, which simplifies the complex problem and provides direction for the subsequent study of hydrazine and its derivatives.
5 Conclusion
In this work, the relationship between the structure of hydrazine derivatives and β radiolysis stability was studied resulting in a quantitative structural activity relationship equation. The most stable configurations of the five hydrazine derivatives was obtained by geometrical optimization and energy calculations. The relationships between the half lifetimes of radiolysis and structural parameters of the hydrazine derivatives were analyzed using the Pearson correlation coefficient, and there was significant correlation with the total energy and molecular mass. There was also strong correlation with the molecular superficial area, molecular refractivity, molecular volume, lowest unoccupied molecular orbital energy, orbital transition energy, etc. The 3D-QSAR structure equations between the structure and rate of radiolysis were also established using SPSS software. According to the model equation, the total energy of the molecule and orbital transition energy are the main factors that affect the stability of hydrazine derivatives under beta radiolysis, and the correlation is negative. In other words, the physicochemical parameters that describe the properties of hydrazine derivatives were obtained using quantitative calculations and related chemistry software. It can be concluded that the hydrazine derivatives with higher total molecular and orbital transition energies are more stable under beta irradiation, as evidenced by the slower β radiolysis rate.
Hydrazine and its derivatives
, Kirk-othmer Encyclope of Chemical Technology. 19, 562-607 (2004).Inorganic hydrazine derivatives: synthesis. properties and applications
. Inorganic Hydrazine Derivatives, 219-252 (2014).The crystal structure of hydrazine
. Acta. Cryst. 4, 10-14 (1951). doi: 10.1107/S0365110X51000027History of liquid propellant rocket engines in the united states
. J. Propul. Power. 19 (6), 978-1007 (2003). doi: 10.2514/2.6942Liquid fuels and propellants for aerospace propulsion: 1903-2003
. J. Propul. Power. 19 (6), 1089-1107 (2003). doi: 10.2514/2.6946Characterization of iridium catalyst for decomposition of hydrazine hydrate for hydrogen generation
. Catal. Lett. 109:181-186 (2006). doi: 10.1007/s10562-006-0081-3Review of recovery of Pu(Ⅳ) in organic waste
. Journal of Nuclear Radiochemistry. 38, 79-85 (2016). doi: 10.7538/hhx.2016.38.02.0079 (in Chinese)The investigation of the quantitative structure-activity relationships between the structure of hydrazine derivatives and the reduction of Np(VI)
. J. Radioanal. Nucl. Ch. 311:1565-1575 (2017). doi: 10.1007/s10967-016-5137-2Dosimeter for measuring exposure to hydrazine and hazardous hydrazine derivatives
. Atmos. Environ. 23: 1989 (1967). doi: 10.1016/0004-6981(89)90036-XMutagenicity and toxicity of carcinogenic and other hydrazine derivatives:correlation between toxic potency in animals and toxic potency in Salmonella typhimurium TA1538
. App. Environ. Micro. 44, 801-808 (1982).Spectrophotometric determination of hydrazine, methylhydrazine, and 1,1-dimethylhydrazine with preliminary derivatization by 5-nitro-2-furaldehyd
. J. Anal. Chem, 2017, 72: 171-177. doi: 10.1134/S106193481702006XEffects of hydrazine and its derivatives on the development of intestinal brush border enzymes
. Toxicol. App. Pharm. 49, 305-311 (1979). doi: 10.1016/0041-008X(79)90255-2The platinum catalyzed hydrazine decomposition in non-nitrate acidic media
. App. Catal. A. 242, 1-10 (2003). doi: 10.1016/S0926-860X(02)00368-XProgress in treatment technology of unsymmetrical dimethylhydrazine wastewater
. Chemical propellants & polymeric materials. 13, 21-25 (2015).Photocatalytic oxidation of 1,1-dimethyl hydrazine vapours on TiO2: FTIR insitu studies
. Catal. Today. 122,178-185 (2007). doi: 10.1016/j.cattod.2007.01.029Low-temperature plasma-induced degradation of aqueous 2,4-dinitrophenol
. J. Hazard. Mater. 154, 506-512 (2008). doi: 10.1016/j.jhazmat.2007.10.053Low-temperature plasma processing of materials: past, present, and future
. Plasma. Process. Polym. 2,7-15 (2005). doi: 10.1002/ppap.20040007Carbon catalyzed supercritical water oxidation of phenol
. J. Supercrit. Fluid. 22, 149-156 (2002). doi: 10.1016/S0896-8446(01)00112-7Supercritical water oxidation research and development update
. Environ. Prog. Sustain. 14 (3), 182-192 (1995). doi: 10.1002/ep.670140318Treatment of dye aqueous solutions using nanofiltration polyamide composite membranes for the dye wastewater reuse
. Dyes. Pigments. 76, 429-434 (2008). doi: 10.1016/j.dyepig.2006.09.007Radiation technology and environmental protection
, Atomic Energy. 70, 438-443 (1991). doi: 10.1007/BF01124523Use of quantitative structure-activity relationships (QSAR) in drug design (review)
. J. Pharm. Chem. 14, 678-691 (1980). doi: 10.1007/BF00765654A review of quantitative structure-activity relationship (QSAR) models
. Hi-Stat Discussion Paper. Series. 0719, 1-35 (2007).The investigation of the quantitative structure-activity relationships between the structure of hydrazine derivatives and the reduction of Pu(IV)
. Nucl. Tech. 39, 080302 (2016). (in Chinese) doi: 10.11889/j.0253-3219.2016.hjs.39.080302