Research Article  Open Access
Zhihui Ni, Qiang Zeng, Wu Lichun, "Determination of the Sediment Carrying Capacity Based on Perturbed Theory", The Scientific World Journal, vol. 2014, Article ID 240858, 10 pages, 2014. https://doi.org/10.1155/2014/240858
Determination of the Sediment Carrying Capacity Based on Perturbed Theory
Abstract
According to the previous studies of sediment carrying capacity, a new method of sediment carrying capacity on perturbed theory was proposed. By taking into account the average water depth, average flow velocity, settling velocity, and other influencing factors and introducing the median grain size as one main influencing factor in deriving the new formula, we established a new sediment carrying capacity formula. The coefficients were determined by the principle of dimensional analysis, multiple linear regression method, and the least square method. After that, the new formula was verified through measuring data of natural rivers and flume tests and comparing the verified results calculated by Cao Formula, Zhang Formula, Li Formula, EngelungHansen Formula, AckersWhite Formula, and Yang Formula. According to the compared results, it can be seen that the new method is of high accuracy. It could be a useful reference for the determination of sediment carrying capacity.
1. Introduction
Sediment carrying capacity reflects the account of entrainment and transportation by the flow under the certain boundary condition. It is a comprehensive index characterizing the sediment carrying capacity of flow under the conditions of equilibrium of scouring and deposition. For a long time, the sediment carrying capacity is an important issue in the field of sediment study; it involves sediment transport, pollutant diffusion, evolution of beach, sea bed erosion, and many other issues.
For half a century, scholars at home and abroad had done a lot of studies on the sediment carrying capacity and got a series of results [1â€“4]. At present, the more common sediment carrying capacity formula is obtained by Zhang [5] because of its harmony dimension, clear mechanical mechanism, and well expressed contradictory relationship between turbulence and gravity effects. Zhang Formula is derived by the measured data of Yangtze River, Yellow River, channels, and flume tests based on the hypothesis of turbulence restriction. But in practical engineering applications, it is difficult to find a unified and value. Yang [6] established a new sediment carrying capacity formula including the sandy bed load, and this formula is based on the research of unit stream power by theoretical models. Also, Yang made detailed analyses and comparisons between Velikanovâ€™s gravitational power theory [7] and his unit stream power theory from theoretical point of view using field and laboratory data for verification in â€śSediment Transport Theory and Practiceâ€ť [8]. By considering the bed load transport rate and combining it with the law of suspended sediment along the vertical distribution, Wang [9] obtained the sediment carrying capacity formula of suspended sediment. Based on the theory of sediment movement and energy balance, Xing et al. [10] derived a halftide average sediment carrying capacity formula at Yangtze river Delta that combined the researches of Dengâ€™s and Douâ€™s. There are so many sediment carrying capacity formulas like the ones above [11, 12]. From the above analysis, it is easy to find that most of the sediment carrying capacity formulas have adopted two common hydraulic factors (velocity and depth) and sediment factor (settling velocity), and the difference is the index of each factor.
By considering the impact of sediment grain size, Mai and Zhao [13] proposed a sediment carrying capacity formula based on the measured data of Yellow River. Chien and Wan [14] summarized and analyzed the Einstein Formula, Bagnold Formula, EngelundHansen Formula, AckersWhite Formula, Zhang Formula, Sha formula, and other formulas and suggested that the relative roughness is a main influencing factor in deriving the sediment carrying capacity in â€śMechanics of Sediment Transport.â€ť
The perturbed theory started from 1882. It is closely united to fluid mechanics and other mechanics disciplines. Ni and Zhou [15] obtained the law of suspended sediment concentration by momentum equation based on the perturbed theory.
Internationally, there are so many achievements on sediment carrying capacity. EngelundHansen Formula, AckersWhite Formula, and Yang Formula are commoner ones. Engelund and Hansen [16] assumed that the mean velocity of sediment is proportional to friction velocity and obtained a new sediment carrying capacity formula in the stage of dune and sand waves by introducing flow intensity parameter and sediment transport intensity parameter , , , , where is the Darcy coefficient, is sediment transport intensity parameter, is flow intensity parameter, and is singlewide transport rate of total sediment. Different to Engelund and Hansen, Ackers and White [17] assumed that the efficiency of sediment carrying by the flow is related to sediment movement intensity, effective tractive force, and weight of surface sediment particles in the water, and then they verified their conjecture by a large number of measured data. At last, they derived a new sediment carrying capacity formula based on the flume tests, , , where is sediment transport parameter, is flow intensity parameter, , , , and are fitting parameters, is friction velocity, is mean velocity of each profile, is sediment grain size, are bulk densities of sediment and water, and is water depth. Thinking of the theoretical model of unit stream power, Yang [6] derived a new bed load sediment carrying capacity formula, , where is sediment carrying capacity, is the settling velocity, is starting velocity, and , , , , , and are the fitting parameters.
Also, some famous sediment carrying capacity formulas had been derived in China. According to the hypothesis of turbulence restriction, Zhang [5] established the energy balance equations and then derived the formula () based on the principle of energy balance, where is the coefficient including dimension (kg/m^{3}), is the index, is the average velocity (m/s), is the gravity acceleration (m/s^{2}), is the water depth (m), and is the settling velocity (mm/s). Cao et al. [18] improved the Zhang Formula according to the balance of saturated sediment data, and he obtained a new sediment carrying capacity formula (), where is bulk density of sediment (g/cm^{3}), is the bulk density of muddy water (g/cm^{3}), is the hydraulic radius (m), and is the corresponding settling velocity (mm/s). Li [19] thought that there is certain linear relationship between and relative roughness and proposed a new sediment carrying capacity formula () based on it, where , , and is average velocity of water vertical line (m/s).
The paper will analyze the major factors by using perturbed theory and then derive a new sediment carrying capacity formula. After that, the measured data of home and abroad rivers will be used to verify and analyze the new formula. Simultaneously, some tank tests by famous scholars should be used for further validation.
2. Methodology Description
2.1. Introduction of Perturbed Theory
Perturbed theory gradually developed in solving nonlinear equations with small parameter. It differs from conventional ones because the main task of the evolutionary method established on the perturbed theory is to solve the series of small parameter in expansion. Then the original problem is transformed to gradually solve the coefficient of the expansion. In this way, the value will be close to the true value, and it is easier to other iterative methods and always the reason can be solved.
When using perturbed theory to solve the problem, firstly we must utilize dimensionless method to analyze the equation and then propose the perturbed variable. For example, we utilized dimensionless method to analyze the influencing factor , and at the same time, we proposed the perturbed variable ; then we adopted perturbed theory to analyze the influencing factor ; after that we have
Due to the gradualness of asymptotic series, usually taking a few can get very high precision and sometimes even reach the exact solution. At present, we usually take the firstorder or secondorder type of the analytic expression.
2.2. Derivation of the New Equation
The factors influencing the sediment carrying capacity included velocity (), water depth (), gradient ratio (), gravity action result in flow (), viscosity (), effective bulk density (), settling velocity (), median grain size (), and composition of the river bed and river width (). Therefore, the sediment carrying capacity formula can be written as a function of the following form:
Taking , , and as the basic variables, then (2) can be written as a dimensionless formula:
For the natural river, is a constant, the effects of flow viscosity and the side walls. Therefore, (3) can be converted into
According to the home and abroad researches it can be found that, in fact, the basic expression of most experienced and semiempirical formulas derived from natural rivers is (4). This paper will analyze by perturbed theory combining the contrast relation between turbulence intensity of the flow () and gravity action () with the contrast relation between turbulence intensity of the flow () and relative roughness ().
Analyzing the contrast relation between turbulence intensity of flow () and the relative roughness () on perturbed theory, we get where is the perturbation variable.
Analyzing the contrast relation between turbulence intensity of flow () and gravity action () on perturbed theory, we get where is the perturbation variable.
Equations (5) and (6) can be transformed into whereâ€‰â€‰ and are small values; therefore, and are smaller. Because of that, the values of and could be neglected. Combining the above expressions, (7) can be written as follows:
To make (8) more intuitive, we use replacing , replacing , and replacing ; then we can get
2.3. Determination of Physical Parameters
According to Ni et al. [36â€“38], the least square method combination of enumeration can be used to determine the coefficients. Concrete steps are as follows.(1)The establishment of objective function â€‰where , , , , and represent the th layer flow velocity, median grain size, water depth, settling velocity, and sediment carrying capacity. is the total number of the layers considered.(2) When the objective function reaches the minimum, the following conditions are satisfied: â€‰where is the total number of the layers for fitting. Solving (11) with the matrix method gives , , and .
Analyzing 1633 domestic and foreign measured data by multivariate linear regression analysis combination of the least square method, the values of coefficients can be obtained as follows:
Then the sediment carrying capacity formula can be written as follows: where is the sediment carrying capacity (kg/m^{3}); is average flow velocity (m/s); is the median grain size (mm); is the average water depth (m); is the average settling velocity (mm/s).
Then, (13) is the formula of the sediment carrying capacity.
3. Verification and Comparison
The measured data covers several types of rivers at home and abroad. In the paper, we will adopt 2633 measured data to fit and verify the derived formula. The parameter values of 2633 measured data in different rivers are shown in Table 1.

These natural rivers include the Niobrara River (see [20]), the Atchafalaya River (see [21]), the Rio Grande River (see [22]), the Missouri River (see [23]), the Sacramento River (see [24]) in USA, the North Saskatchewan River and the Elbow River in Canada (see [25]), Columbian River (see [26]), and some rivers in India and Pakistan [27, 28, 39].
In 1955, Cobly and Hembree observed the Niobrara River near Lusk and won the measured data of flow velocity, median grain size, water depth, sediment concentration, and so on. Niobrara River is 720â€‰km long and the basin area is 31,080â€‰km^{2}. In 1968, Toffaleti observed the Atchafalaya River, a tributary of the Red River, and the Mississippi River and got the measured data. From the observed results it can be found that the Atchafalaya River is 360â€‰km long, the range of flow velocity is 0.21 to 3.56â€‰m/s, and the range of water depth is 4.163 to 13.28â€‰m; this observed data had higher reliability. In 1972, Culbertson et al. observed the Rio Grande River, one of the largest rivers in North America; it is about 3,033â€‰km long and obtained some valuable measured data. In 1978, Shen et al. observed the Missouri River and got some measured data. From the observed results it can be found that this river is over 4300â€‰km long, the range of flow velocity is 1.389 to 1.745â€‰m/s, and the range of water depth is 2.714 to 4.238â€‰m. In 1990, Nakato observed the Sacramento River, the longest river in California, and obtained some useful measured data. From the observed results it can be seen that the river is 614.8â€‰km long and the ranges of flow velocity, water depth, and median grain size are wide; they are suitable for fitting formula.
To verify the reliability of the derived formula, the paper will use the remaining 1000 measured data to verify and analyze the derived formula. The results of verification of all formulas can be seen in Figures 1, 2, 3, 4, 5, 6, and 7.
In the figures, axis is the measured value and axis is the calculated value.
According to the verified results it can be seen that the fitting degree of derived formula is higher than other formulas. Cao Formula and Yang Formula are better. Also, combined with correlation coefficients and theories of all formulas it can be found that the formula taking into account the median grain size in calculation has high precision. Therefore, we need to take further verification of the derived formula whether the median grain size is a main factor influencing the sediment carrying capacity.
In order to further verify the derived formula, the paper will use the results of flume tests of Gilbert [30], Brownlie [40], Einstein [41], EinsteinChien [31], Brooks [32], Guy et al. [33], Song et al. [34], and Wang et al. [35]. The parameter values of measured data in the flume tests by the scholars are shown in Table 2.

By taking into account the change in difference range of median grain size, Gilbert (1914) had done several flume tests in different size of flumes; the range of median grain size is 0.305~7.01, and the range of different size flume is 0.2012~0.5974. Similar to Gilbert (1914), MeyerPeterMuller [42] had done three sets of tests, respectively: when the range of median grain size is 0.52 to 5.2â€‰mm, the corresponding flume width is 0.3539â€‰m; when the median grain size is 3.3â€‰mm, the corresponding flume width is 0.4999â€‰m; when the range of median grain size is 0.38 to 28.65â€‰mm, the corresponding flume width is 1.9998â€‰m. The flume width in Einstein (1950) flume tests is 0.2667â€‰m, and the range of corresponding median grain size is 0.11 to 0.91â€‰mm. With Chien, Einstein had done other flume tests in 1953; the flume width is 0.3048â€‰m, and the range of corresponding sediment grain size is 0.095 to 0.385â€‰mm. The flume width in Brooks [32] flume tests is 0.2667â€‰m, and the range of corresponding sediment grain size is 0.088~0.145â€‰mm. Guy et al. [33] had done two sets of flume tests in different flume: when the range of median grain size is 0.32 to 0.59â€‰mm, the corresponding flume width is 0.6096â€‰m; when the median grain size is 0.19~1.20â€‰mm, the corresponding flume width is 2.4384â€‰m. The flume width in Song et al. (1998) flume tests is 0.5998â€‰m, and the corresponding median grain size is 12.3â€‰mm. The flume width in Wang et al. (1998) flume tests is 0.5â€‰m, and the corresponding median grain size is 0.18 to 13â€‰mm.
The results of verification of all formulas can be seen in Figures 8, 9, 10, 11, 12, 13, and 14.
In the figures, axis is the measured value and axis is the calculated value. According to the verified results of the flume tests it can be found that the result of derived formula is the closest to the measured data. The further verification illustrated the rationality of the derived formula, and it can be a reference in engineering.
4. Error Analysis
According to the verified results, the errors by the formulas are shown in Tables 3 and 4.
 
represents the percentage of scatters in which the error is less than 10% of measured data in the natural rivers; represents the percentage of scatters in which the error is less than 20% of measured data in the natural rivers. 
 
represents the percentage of scatters in which the error is less than 10% of measured data in the natural rivers; represents the percentage of scatters in which the error is less than 20% of measured data in the natural rivers. 
From the tables it can be found that the error by derived formula is the least compared to the other six formulas. When verifying with the measured data of natural rivers, the average error is 4.05%; when verifying with the measured data of flume tests, the average error is 4.57%. If calculated by other formulas, the average errors are in the range of 7.16%~37.66% and 6.88%~13.66%, respectively. Taking into account the scale effect of flume tests, it is understandable that the error by flume tests is bigger than that by natural rivers data. From the error analysis it can be seen that it is feasible to take the median grain size as one main influencing factor.
Meanwhile, according to the error analysis data it can be also found that the percentage of scatters in which the error is less than 10% and less than 20% verified by derived formula are less than those by the others formulas, respectively, are 88%, 96% and 85%, 92%.
In conclusion, the derived formula is simple in principle, practical in implementation, and reasonable in results, and it further illustrates that the view of considering the median grain size in determining the sediment carrying capacity is correct. This view is consistent with Chien.
5. Conclusions
According to the present analysis and study, some conclusions can be drawn as follows.(1)By taking into account the influence of median grain size on the sediment carrying, a new sediment carrying capacity formula combined with the turbulence intensity of the flow (), gravity action (), and relative roughness () is established.(2)From the verified results it can be seen that the derived formula has high accuracy by using the method of dimensional analysis and combined with multivariate analysis and linear least squares method, it is a strong theoretical and practical method.(3)The results further verified by flume tests showed that the precision by introducing the relative roughness in deriving formula is higher than other formulas. It means that taking into account the median grain size in sediment deposition is reasonable and feasible.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper.
Acknowledgments
The study was supported by the National Natural Science Foundation of China (Grant no. 41306078), the Natural Science Foundation Project of CQ CSTC (Grant no. cstc2012jjA30002), the Science and Technology Project of Chongqing Education Committee (Grants no. KJ110409), the National Engineering Research Center for Inland Waterway Regulation Program (Grant no. SLK2012A02), and the National Key Technology R&D Program (Grant no. 2012BAB05B03).
References
 M. H. Yu, G. L. Yang, G. F. Liu et al., â€śA preliminary study on the formula of nonuniform sediment carrying capacity,â€ť Journal of Sediment Research, no. 3, pp. 25â€“29, 2001 (Chinese). View at: Google Scholar
 R. T. Milhous, â€śClimate change and changes in sediment transport capacity in the Colorado Plateau, USA,â€ť Sediment Budgets, vol. 2, no. 292, pp. 271â€“278, 2005. View at: Google Scholar
 S. Q. Yang, S. C. Koh, I. S. Kim, and Y. Song, â€śSediment transport capacity: an improved Bagnold formula,â€ť International Journal of Sediment Research, vol. 22, no. 1, pp. 27â€“38, 2007. View at: Google Scholar
 G. Q. Wang, â€śAdvances in River Sediment Research,â€ť Journal of Sediment Research, vol. 2, pp. 64â€“81, 2007 (Chinese). View at: Google Scholar
 R. J. Zhang, River Mechanics, China Industry Press, Beijing, China, 1961.
 C. T. Yang, â€śIncipient Motion an d Sediment Transport,â€ť Journal of Hydraulic Engineering, vol. 99, no. 10, pp. 1679â€“1704, 1973. View at: Google Scholar
 M. A. Velikanov, â€śGravitational theory of sediment transport,â€ť Journal of Science of Soviet Union Geophysics, vol. 4, pp. 349â€“359, 1954. View at: Google Scholar
 C. T. Yang, Sediment Transport Theory and Practice, Krieger Publishing Company, 2003.
 S. Q. Wang, â€śThe movement of sand waves and the exchanges of bed material,â€ť Journal of Sediment Research, no. 4, pp. 14â€“23, 1992 (Chinese). View at: Google Scholar
 Y. Xing, Z. Y. Song, J. Kong et al., â€śA preliminary study on sediment carrying capacity formula of Yangtze River Estuary,â€ť Journal of China Hydrology, vol. 28, no. 1, pp. 64â€“66, 2008 (Chinese). View at: Google Scholar
 G. R. Dou, F. W. Dong, and X. P. Dou, â€śThe sediment carrying capacity of waves and tidal,â€ť Chinese Science Bulletin, vol. 40, no. 5, pp. 443â€“446, 1995. View at: Google Scholar
 J. Guo, â€śLogarithmic matching and its applications in computational hydraulics and sediment transport,â€ť Journal of Hydraulic Research, vol. 40, no. 5, pp. 555â€“565, 2002. View at: Publisher Site  Google Scholar
 Q. W. Mai and S. L. Zhao, â€śPreliminary study on the effect of finest sediment content on the sediment carrying capacity of the Yellow River,â€ť International Journal of Sediment Research, no. 2, pp. 1â€“39, 1981 (Chinese). View at: Google Scholar
 N. Chien and Z. H. Wan, Mechanics of Sediment Transport, American Society of Civil Engineers, 1999.
 J. R. Ni and D. H. Zhou, â€śInterpretation of vertical sediment distribution in terms of perturbation theory,â€ť Journal of Hydraulic Engineering, vol. 5, pp. 1â€“5, 1999. View at: Google Scholar
 F. Engelund and E. Hansen, â€śComparison between similarity theory and regime formulae,â€ť Basic ResearchProgress Report 13, Hydraulic Laboratory, Technical University of Denmark, 1967. View at: Google Scholar
 P. Ackers and W. R. White, â€śSediment transport: new approach and analysis,â€ť Journal of the Hydraulics Division, vol. 99, no. 11, pp. 2041â€“2060, 1973. View at: Google Scholar
 R. X. Cao, W. Chen, S. Q. Qian et al., â€śPreliminary approach to river bed blowing up scouring by hyperconcentration flood,â€ť Yellow River, no. 2, pp. 1â€“6, 1997 (Chinese). View at: Google Scholar
 Y. T. Li, â€śA preliminary study on plane deformation calculation of alluvial channel,â€ť Journal of Sediment Research, vol. 1, pp. 34â€“44, 1988 (Chinese). View at: Google Scholar
 B. R. Cobly and C. H. Hembree, Computation of Total Sediment Discharge, Niobrara River near Cody, Nebraska, U.S. Geological Survey, 1955, Water Supply Paper 1357.
 F. B. Toffaleti, â€śA procedure for computation of the total river sand discharge and detailed distribution, bed to surface,â€ť Tech. Rep. 5, Committee of Channel Stabilization, Corps of Engineering, U.S. Army, 1968. View at: Google Scholar
 J. K. Culbertson, C. H. Scott, and J. P. Bennett, Summary of AlluvialChannel Data From Rio Grande Conveyance Channel, New Mexico, 1965–69, U.S. Geological Survey Professional Paper, Washington, DC, USA, 1972.
 H. W. Shen, W. J. Mellems, and A. S. Barrison, â€śTemperature and Missouri river stages near Omaha,â€ť Journal of the Hydraulics Division, vol. 104, no. 1, pp. 1â€“20, 1978. View at: Google Scholar
 T. Nakato, â€śTests of selected sedimenttransport formulas,â€ť Journal of Hydraulic Engineering, vol. 116, no. 3, pp. 362â€“379, 1990. View at: Publisher Site  Google Scholar
 G. W. Samide, Sediment Transport Measurements [M.S. thesis], University of Alberta, 1971.
 NEDCO, â€śRio Magdalena and Canal del Dique project, mission tecnica ColomboHolandesa,â€ť Tech. Rep., NEDECO, The Hague, The Netherlands, 1973. View at: Google Scholar
 S. V. Chitales, Hydraulics of Stable Channels, Government of India, Ministy of Irrigation and Power, Central Water and Power Commission, 1966.
 H. M. Chaudhry, K. V. H. Smith, and H. Vigil, â€śComputation of sediment transporty in irrigation canals,â€ť Proceeding of the Institution of Civil Engineering, vol. 45, pp. 79â€“101, 1970. View at: Google Scholar
 K. Mahmood, M. I. Haque, A. M. Choudri, and M. A. Malik, ACOP Canals Equilibrium Data, Volume 10: Summary of 1974–1980 Data, Civil, Mechanical and Enviromental Engineering Department, George Washington University, 1984.
 G. K. Gilbert, â€śThe transportation of debris by running water,â€ť U.S. Geological Survey, Professional Paper 86, 1914. View at: Google Scholar
 H. A. Einstein and N. Chien, â€śTransport of sediment mixtures with large ranges of grain sizes,â€ť MRD 2, University of California, 1953. View at: Google Scholar
 N. H. Brooks, â€śMechanics of streams with movable beds of fine sand,â€ť Transactions of the American Society of Civil Engineers, vol. 23, pp. 526â€“549, 1958. View at: Google Scholar
 H. P. Guy, D. B. Simons, and E. V. Richardson, Summary of Alluvial Channel Data from Flume Experiments, 1956–1961, U.S. Geological Survey, Professional Paper 462I, 1966.
 T. Song, Y. M. Chiew, and C. O. Chin, â€śEffect of bedload movement on flow friction factor,â€ť Journal of Hydraulic Engineering, vol. 124, no. 2, pp. 165â€“175, 1998. View at: Publisher Site  Google Scholar
 Z. Y. Wang, Y. N. Xu, and X. B. Su, â€śA study on channel scour rate of sedimentladen flow and river bed inertia,â€ť Journal of Sediment Research, no. 2, pp. 1â€“9, 1998 (Chinese). View at: Google Scholar
 Z. H. Ni, Z. Y. Song, and L. C. Wu, â€śA new model for velocity profiles in estuarine waters,â€ť China Ocean Engineering, vol. 24, no. 2, pp. 255â€“264, 2010. View at: Google Scholar
 Z. H. Ni, L. C. Wu, and Z. Y. Song, â€śA method of wind velocity distribution above the sea surface in nearshore waters,â€ť Acta Aerodynamica Sinica, vol. 28, no. 5, pp. 547â€“552, 2010. View at: Google Scholar
 N. Zhihui, W. Lichun, Z. Xujin, Z. Qiang, W. Minghui, and Y. Jing, â€śA new model for the regulation width of waterway based on hydraulic geometry relation,â€ť Mathematical Problems in Engineering, vol. 2013, Article ID 182656, 8 pages, 2013. View at: Publisher Site  Google Scholar
 D. B. Simons, Theory and design of stable channels in alluvial materials [Ph.D. thesis], Colorado State University, Fort Collins, Colo, USA, 1957.
 W. R. Brownlie, â€śCompilation of alluvial channel data: laboratory and field, W. M. Keck laboratory of hydraulics and water resources,â€ť Report No. KHR43B, California Institute of Technology, Pasadena, Calif, USA, 1981. View at: Google Scholar
 International Research and Training Centre on Erosion and Sedimentation, â€śSediment transport data in laboratorial flumes,â€ť IRTCES Circular, vol. 2, pp. 1â€“102, 1987. View at: Google Scholar
 E. MeyerPeter and R. Muller, â€śFormulas for bedload transport,â€ť in Proceedings of the 2nd Meeting of the International Association for Hydraulic Structures Research, pp. 39â€“64, Stockholm, Sweden, 1948. View at: Google Scholar
Copyright
Copyright © 2014 Zhihui Ni et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.