#### Abstract

Water hammer analysis is a fundamental work of pipeline systems design process for water distribution networks. The main characteristics for mine drainage system are the limited space and high cost of equipment and pipeline changing. In order to solve the protection problem of valve-closing water hammer for mine drainage system, a water hammer protection method for mine drainage system based on velocity adjustment of HCV (Hydraulic Control Valve) is proposed in this paper. The mathematic model of water hammer fluctuations is established based on the characteristic line method. Then, boundary conditions of water hammer controlling for mine drainage system are determined and its simplex model is established. The optimization adjustment strategy is solved from the mathematic model of multistage valve-closing. Taking a mine drainage system as an example, compared results between simulations and experiments show that the proposed method and the optimized valve-closing strategy are effective.

#### 1. Introduction

Mine drainage system is an important part in the safety production of coal mine [1]. Due to the limited space and the high cost of equipment changing, water hammer is a common phenomenon in mine drainage system and its harm is inestimable [2, 3]. The minor injuries of water hammer can cause severe shock or even pipes breaking, and the major injuries can cause equipment damaging, pumping station flooding, or even injuries to underground staff. There are several traditional water hammer protection methods such as installing vacuum valves, exhaust valve, and pressure tank. However, when these water hammer protection methods are used to mine drainage system, the original piping arrangement has to be changed because additional equipment is needed. It is unenforceable to do this in this limited mine space and the cost is too high. In this situation, controlling the time and velocity of valve-closing is an effective means to protect water hammer in mine drainage system. A water hammer is easily to be formed in pipeline if the valve-closing is fast. On the contrary, the capability of pumps is inefficient if the valve-closing is too slow because pumps in mine drainage system are centrifugal. This inefficient operation has damage for pumps, because power provided by pumps is converted into heat. Therefore, it is important for us to research the critical valve-closing velocity for Hydraulic Control Valve (HCV). Researches of theoretical system and engineering applications for water hammer protection of pipeline fluid delivery areas are increasingly improved [4–6]. But it is still a difficult problem in mine drainage system because of the limited mine space [7, 8].

At present, the theories and methods of water hammer and its protection means become more and more mature. In summary, researches about this problem are mainly in three aspects: hydraulic transient modeling, the calculation, and protection methods of water hammer.

*(1**) Hydraulic Transient Modeling*. For the research on hydraulic transient, from the mathematical derivation of the 18th century to the graphical analysis of the mid-20th century and to the current computer digital simulation, scholars have already made a lot of research results. The major achievements are getting the relationship between multiphase and multicomponent transient flow equations, water hammer equations, and the control equations, such as Joukowsky equation [9]. Based on the transient flow simulation theory, Colombo et al. [10] proposed an aqueducts fault detection technology, Lee et al. [11] proposed the pipe network leak and deterioration over time detection technology by the time domain reflectometry (TDR), Arbon et al. [12] proposed pipeline corrosion and blockage detection technology, Gong et al. [13] proposed a detection technology for pipe friction, wall thickness, velocity, position, and the length of the pipes, and Ferrante et al. [14] presented a leak detection method with coupling wavelet analysis and a Lagrangian model techniques. Meniconi et al. [15] presented a pipe system diagnosis method with the small amplitude sharp pressure waves.

*(2**) The Calculation Methods of Hydraulic Transient for Water Hammer*. The calculation methods of hydraulic transient for water hammer include arithmetic method, graphic method, and numerical method.

*(i) Arithmetic Method*. Before the 1930s, the hydraulic transient calculation of water hammer used Allievi equations mostly [16]. Allievi equations can be called Arithmetic method which is used to solve the problems of water hammer that with simple boundary conditions and its workload is very large.

*(ii) Graphic Method*. Graphic method is developed in 1930s to 1960s. Bergeron, Parmakian, and so forth [17] are committed to develop this method. Boundary conditions and the process of water hammer fluctuation are expressed through coordinate graphics of and according to this method. Due to the graphics, it is simple and intuitive for the hydraulic transient calculation of water hammer. However, the accuracy is not high because this method is restricted by calculating means and assumptions.

*(iii) Numerical Method*. From 1960s, some numerical methods appeared that can be aided by computers, such as Characteristic (MOC) [18], Wave Characteristic Method (WCM) [19], Implicit Method [20, 21], and Finite Element Method (FEM) [22, 23]. The WCM can solve water hammer problems of complex piping systems and boundary conditions. It is the most common method because of the high accuracy and computing. The Implicit Method divides pipeline into several segments and solves equations of the entire pipeline system simultaneously in each segment. The advantage of Implicit Method can be described in a way that a longer segment is selected and the number of calculations is reduced. However, there is more time needed for calculation in large and complex pipeline network system [24]. FEM with flexibility is used in pipe network system which have complex boundary conditions. However, it has a limitation in solving hydraulic transient problems.

*(3**) Water Hammer Protection*. Prevention and controlling means for water hammer are researched with the development of their theory and calculation. Wylie [25–27] has researched several protective devices for water hammer, such as air valves, check valves, pressure tank, and surge tank. Lee [28] and Stephenson [29] discussed the performances of air valves in water hammer protection. However, these researches have not given a quantitative calculation for the problem of water hammer protection.

This paper proposed a water hammer protection method based on velocity adjustment of HCV to deal with the problem of valve-closing water hammer in mine drainage system. The mathematic model of water hammer fluctuations is founded based on MOC according to the hydraulic transient. Then, boundary conditions of water hammer controlling for mine drainage system are determined and the simplex model is established. Finally, the optimization adjustment strategy is solved with simulation and experiment.

The remainder of this paper is organized as follows. Section 2 provides the mathematic model of the propagation and superposition model for water hammer fluctuations. Section 3 provides an optimization method to determine the adjustment strategy for HCV. Section 4 presents a case study. Concluding remarks are offered in the last section of this paper.

#### 2. The Process of Water Hammer Fluctuations

##### 2.1. Mathematic Model

*(1**) Foundation Equations*. The momentum equations of transient flow can be expressed aswhere is longitudinal mean velocity, is pressure head, is friction factor, is diameter of pipeline, is gravitational acceleration, is propagation distance along pipe, and is time.

The left of the equation represents the time-varying inertia force in a unit volume. In the right, the first item represents the pressure in a unit volume of fluid. The second item is friction loss pressure in a unit length of pipe, and the third item is pressure of transient flow.

The continuity equation is used to describe transient flow and it can be expressed as follows:where is propagation speed of water hammer wave and is the angle between the axis of pipe and the horizontal line.

*(2**) The Process of Water Hammer Fluctuations*. Pressures that are caused by friction loss and transient flow velocity are ignored. Furthermore, is considered for continuity equations. Then, the foundation equations can be simplified as

The second-order partial differential equations about and are obtained through the partial derivative and are taken for variables and in (3). Therefore, the momentum equations for water hammer wave in pipeline drainage system are expressed as

The functions and are introduced along with characteristic lines . Then, the general solution of (4) can be described aswhere is the initial pressure head, is the initial flow velocity, and and are direct and reflection fluctuation functions.

##### 2.2. Propagation and Superposition Model of Water Hammer Wave Based on Characteristic Line Method

Pumps and valves are the generating sources of water hammer wave in hydraulic transition process of pipeline fluid delivery system. These two kinds of water hammer wave start at the same time and propagation directions are superimposed. The superposition leads to strengthening of fluctuations for water hammer. In order to control the water hammer pressure effectively, the superposition effects of water hammer wave need to be weakened as much as possible.

The fundamental equations of water hammer can be rewritten as

These equations are combined linearly using an unknown multiplier . Let be the linear combination. The coefficient can be determined by

In the constraint of characteristic line equation , (6) converts to ordinary differential equation as follows:

and are negligible when and the tilt angle of pipeline is less than 25°. The average flow velocity is replaced by flow . The characteristic equations are changed intowhere is the characteristic line of forward wave in -axis, is the characteristic line of reflected wave in -axis, and is cross-sectional area of pipe.

The integral operation and differential conversion are introduced to (9). The discrete characteristic equations of water hammer are obtained as follows. Equation (10) and the characteristic line grid are shown in Figure 1:

By solving the above equations, transient variables and in the node of the pipe over the time are expressed aswhere .

Consider

The initial condition of water hammer is that of a steady flow, and in (11) are known, and the parameters on each time step can be solved.

Hydraulic characteristics of HCV can be described through loss coefficient and capacity coefficient . Consider where is cross-sectional area of the upstream of HCV and is the pressure difference of HCV in upper and lower.

Let be the relative opening displacement of HCV. The dimensionless flow through HCV can be described aswhere is pressure difference on HCV. It can be solved as . The subscript “0” represents HCV being open completely.

#### 3. Optimization Method for Velocity Adjustment of HCV

##### 3.1. Boundary Conditions

Boundary conditions for pump of outlet are shown in Figure 2. If the inlet loss of pump is ignored, the inlet of pool, pump, and HCV can be combined and dealt with as a boundary of characteristic line in foundational equations of water hammer.

The equilibrium condition for water head is established in this boundary. It can be described aswhere is the depth of pool inlet, is the work lift of the pump, is the transient resistance loss of HCV, and is the pressure head of HCV outlet.

We take the balance equations of water head which in the characteristic line and inertia equations of pump unit together. And then the boundary conditions can be expressed as follows based on the characteristic lines of pumps:where the pump-lift is . The flow of pump is . The rotational speed is . The torque is . The subscript represents rated working conditions. Then, (16) is generated as follows:where ; ; and is rotational inertia. and are rotational speed and torque in the initial state. Then, (17) can be rewritten as follows when variables and are introduced:The problem is solved as follows:in which and are solved through (18).

The iterative errors should meet (19). Consider

##### 3.2. Optimization Method of HCV Adjustment

The boundary conditions of multistage valve-closing are introduced in the characteristic line calculation of water hammer for which the best procedure for multistage valve-closing is considered. The displacement of valve-closing at any time can be determined through linear procedure in multistage valve closure.

Through a large number of practical tests, the valve closure with a uniform or uniform variable speed is found to have a small influence on the water hammer in the front 2/3 of valve-closing stroke, while there is a large influence on the later 1/3 of valve-closing stroke. Therefore, the valve closure process is divided into multistage. By doing so, the uniform and uniform variable speeds are both introduced to valve closure process. In the following section we will discuss uniform or uniform variable speed calculation method in valve closure process.

Assuming that the times of valve-closing in multi-stage are , the corresponding displacements of valve-closing are . Then, the displacements of valve-closing at any time are determined as follows:where is displacement of valve-closing at and is the total stroke of valve-closing.

For the valve-closing displacement at any time , we can calculate the dimensionless opening degree using (21) through linear interpolation according to equidistant and its input value:

In the scope of velocity adjustment for HCV, the velocity, displacement, and time of valve-closing have various combinations. This paper adopts a multistage adjustment strategy according to simulation results of water hammer.

The valve closure time for the first stage is equal or close to the time of the flow of pump that is zero. The valve closure time for the following stages is one in four of the previous stage. Consequently, the initial valve closure procedure is determined. Then, search computing is taken to determine the optimization valve closure procedure.

#### 4. Case Study and Experiment

##### 4.1. Working Conditions

The mine with normal water quality was considered. The water normal inflow was 50 m^{3}/h and the maximum inflow was 250 m^{3}/h. Water was discharged into ground by two drain pipes of which the specifications were* D*273 × 12, respectively (one was used and another was spare). The pipes in this existing and antiquated mine drainage system were seamless steel pipe. The allowed pressure was as lower at 6.4 MPa. The wellhead elevation was 997.991 m and the borehole inclination was 20°. There were three multistage centrifugal pumps and two HCVs were used in Figure 3. And the specification of pumps was MD280-65 × 6 ( m^{3}/h, m). One pump was used and others were maintenance when the mine was with a normal inflow, or two pumps were used and one was spare when the mine was with a maximum inflow.

**(a)**

**(b)**

The automatic drainage system was chosen in this mine. Check valves were installed in pump outlet. There are HCVs in pipeline shown in Figure 3(a) where the stroke was 220 mm with the structure of HCV shown in Figure 4. During normal pump stopping, HCVs were closed first, and then pumps were stopped. The water hammer phenomenon was significant in drains pipeline during this process. Therefore, we should take some measures to prevent the water hammer phenomenon.

##### 4.2. Simulations

An algorithm was developed using commercial simulation software Fluent6.3 based on the water hammer theory described in Section 2. The acceleration, amplitude, pipeline pressure, and valve closure displacement were obtained in the cases of uniform and uniform variable speed through transient simulation for water hammer. The correctness of simulation and the proposed theory were then verified against experimental results. We simulated the valve-closing in single-stage with a constant velocity and in multistage with a variable velocity based on the proposed theory.

*(1**) The Simulation of Valve-Closing in Single-Stage with a Constant Velocity*. The velocity of valve-closing in single-stage was determined through repeated simulation. Then, times of equipment in this situation were obtained using (20) to (21) as shown in Table 1.

The simulation results of valve-closing parameters in single-stage with are shown in Figures 5–12.

For Figures 5–12, the static pressure of pipeline was at 3.06 MPa for pump stopping. After the vacuum was completed, the pump was started and HCV was open with a constant velocity. At this moment, pipeline pressure was increased along with the valve-opening displacement of HCV. Then, pressure fluctuations appeared. The valve-opening displacement of HCV was 42.6 mm and the maximum pressure was 3.84 MPa when s. After this moment, the pressure was steady. The valve-opening displacement of HCV was 70.2 mm and the maximum pressure was 3.25 MPa when s. The drainage system pressure reached a steady state along with the increasing of valve-opening displacement. The pressure of pipeline began to drop when s and valve-closing displacement was 53.4 mm. Water hammer phenomenon appeared when s. When s, the maximum boost of water hammer was m and the pressure of pipeline at this moment was 4.77 Mpa. It is 1.47 times the normal discharge pressure (3.25 MPa) and the pressure increased by 50% suddenly. The acceleration and amplitude of pipeline in directions of , , and were changed significantly when water hammer occurred. Maximum acceleration in directions of , , and was 4.83 m/s^{2}, 6.46 m/s^{2}, and 5.45 m/s^{2} and the maximum amplitude was 15.8 mm, 25.8 mm, and 15.3 mm, respectively. Additionally, the total stroke of HCV was assumed to be . Valve-closing with rapid and constant velocity were taken in the former and a variable velocity was taken in the latter according to the simulation results of valve-closing water hammer. Due to the velocity of valve-closing that has a large impact on water hammer parameters, it is an effective approach to protect water hammer by choosing a reasonable velocity for valve-closing.

The allowed pressure of pipe is lower and its value was 6.4 MPa. However, the pressure of pipeline reached 4.8 MPa in the case of 22 s valve-closing from the above analysis. This value is near to the allowed pressure of pipe with the pressure being not high, and the water hammer phenomenon was not significant compared with other mine drainage systems. In order to reduce the pressure of pipeline to protection pipeline, it is necessary to develop water hammer protection.

*(2**) The Simulation of Valve-Closing in Multistage with a Variable Velocity*. The optimization strategy of valve-closing adjustment for HCV is generated using boundary conditions according to optimization method of valve-closing adjustment that is proposed in this paper. It is a multistage valve-closing with a variable velocity. In the front 2/3 of valve-closing stroke, the velocity is constant and its value is . In the later 1/3 of valve-closing stroke, the velocity is a variable (the deceleration is constant) and the initial value is . The total time of valve-closing is 43.823 s. Times of equipment in multistage valve-closing are shown in Table 2.

The simulation results of valve-closing parameters in multistage are shown in Figures 13–20.

The pipeline pressuring, acceleration, and amplitude in the directions of , , and are significantly reduced in multistage valve-closing with a variable velocity. Valve-closing with a constant velocity is started at 46.007 s, and valve-closing with a variable velocity is started at 60.6 s. The HCVs are closed and pumps are stopped at 89.83 s. In this process, the maximum pressure of water hammer is 3.49 MPa. No effect of water hammer occurred under this situation.

##### 4.3. Experiments

The experimental devices include multifunction data acquisition instrument, vibration sensor, pressure sensors, and cable displacement encoder. Multifunction data acquisition instrument has 16 channels and its acquisition frequency is 102.4 kHz. Data of acceleration, amplitude, flow, pressure, and displacement of valve-closing were acquired by this multifunction data acquisition instrument. Vibration sensor was used to measure piping vibration under the influence of water hammer. It was disposed along with the , , and direction of the pipe. Pressure sensors were used to measure pipeline pressure. Its measurement accuracy and range are 0–10 MPa and ±0.5 FS. Cable displacement encoder was used to posit the displacement of valve-closing. The arrangement diagram of sensors and flow chart of pumps are shown in Figures 21 and 22.

HCVs control system was designed in this experiment which is shown in Figure 23. Sensors were used to collect signals of pressure and acceleration in pipeline. PLC unit was used to analyze pipeline pressure and vibration for water hammer on the basis of pressure and acceleration signals. An expected running speed was given to HCVs. And then, voltage, current, power, and fault protection information of pumps, pressure, flow, temperature, vacuum, and HCVs status information were monitored at real time. These multiparameter data were integrated using of fuzzy control theory with calculus and control strategy. The HCVs control system can control starting and stopping of pumps automatically and can monitor the signals in system running process timely.

For the total stroke of valve , the experimental conditions are as follows. In the front 2/3 of the stroke, the velocity was constant and its value was . In the later 1/3, the velocity is variable (the deceleration is constant) and the initial value was . The full valve-closing times were 44 s. Measured values of water hammer parameters in multistage are shown in Table 3.

#### 5. Results and Discussion

The reliability and effectiveness are illustrated through the comparison between the results of simulation and experiment.

*(1**) Results of Valve-Closing in Single-Stage with a Constant Velocity and in Multistage with a Variable Velocity Are Compared*. Pipeline pressures and displacements of valve-opening under these two conditions are shown in Figures 24 and 25, respectively. The red curves represent the variations in single-stage with a constant velocity while the black curves represent the variations for multistage with a variable velocity.

Figures 24 and 25 show a significant water hammer effect in single-stage valve-closing with a constant velocity, and the valve-closing time is shorter (22 s). On the other hand, water hammer phenomenon was found to be considerably improved for multistage valve-closing condition with a variable velocity and the valve-closing time was extended (44 s). It shows that valve-closing in multistage with a variable velocity can control the water hammer phenomenon effective in mine drainage system.

*(2**) Results of the Simulation and Experiment Are Also Compared*. Results for valve-closing with a variable velocity are shown in Table 4.

It can be seen from Table 4 that the simulation for valve-closing with a variable velocity shows good agreement with the experiment. The errors between simulation and experiment are less than 5%. Therefore, the valve-closing in multistage with a variable velocity can reduce the impact from water hammer in mine drainage system.

#### 6. Conclusions and Future Work

Limited space and high cost of equipment changing are the main characteristics for mine drainage system. For the problem of valve-closing water hammer protection in mine drainage system, this paper proposed a protection method based on velocity adjustment of HCVs. The mathematic model of valve-closing water hammer is established based on fundamental equations. The strategy of valve-closing in multistage with a variable velocity is obtained by velocity adjustment strategy of HCVs, and it is effective in protecting the valve-closing water hammer without change piping arrangement or addition equipment. The results of simulations and experiments in different adjustment velocity are compared. We obtained that the water hammer phenomenon is not obvious with valve-closing velocity in the front 2/3 of its stroke, and it is very seriously within the later 1/3. Therefore, we take a two-stage valve-closing strategy. In the first stage, the velocity is constant and its value is . In the second stage, the velocity is variable (the deceleration is constant) and the initial value is .

Extending errors in the theoretical wave speed have an important influence on prediction of the closure rate before a water hammer occurs in the actual pipeline. A very challenging mathematical task is variable. This is beyond the scope of this paper and may be considered as a future work.

#### Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

#### Acknowledgment

This work was supported by the National 12.5 Technology Support Program of China under Grant no. SQ2012BAJY3504.