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, Engelung-Hansen Formula, Ackers-White 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 [14]. 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 half-tide 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, Engelund-Hansen Formula, Ackers-White 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. Engelund-Hansen Formula, Ackers-White 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 single-wide 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/m3), is the index, is the average velocity (m/s), is the gravity acceleration (m/s2), 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/cm3), is the bulk density of muddy water (g/cm3), 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 first-order or second-order 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. [3638], 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/m3); 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 km2. 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], Einstein-Chien [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), Meyer-Peter-Muller [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.

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).