Impulse Response Method for Pipeline Systems Equipped with Water Hammer Protection Devices

Surge protection devices, such as surge tanks and air chambers, have been modeled with the impulse response method for transient analysis of water distribution systems. The lumped inertia model and continuity equation are used to represent nonpipe hydraulic elements. Results of pressure or discharge variations obtained by using the impulse response method and the method of characteristics are in good agreement. The impulse response method provides total pressure and discharge along any pipeline segment by direct integration of the ratio of complex head or complex discharge to a complex downstream discharge, respectively. A modification is proposed so that transition between turbulent and laminar flows can be considered. The representation of hydraulic devices has been incorporated into the impedance matrix method, which was developed for heterogeneous and multilooped pipe network systems. The potential advantages of the proposed method over other conventional approaches were investigated by applying the proposed method to hypothetical pipe network systems.     

Dynamic Memory Computation of Impedance Matrix Method

The computational efficiency of the impedance matrix method has been greatly improved for large pipe networks with various dimensions and complexity. Several numerical methods for solving linear system were modified to deal with the complex domain operation and used into impedance evaluation. Two different memory reduction schemes were developed based on one-dimensional storage and implemented with the biconjugate gradient method and the Gaussian elimination scheme, respectively. A new implementation of the impedance matrix method, namely, the dynamic memory allocation scheme, was introduced to efficiently model hydraulic transients in pipeline systems that have large topological structures. Three hypothetical pipe networks, the multiseries system, the multilooped system, and the multiblock system, were used to test the performance of the developed schemes. The impact of randomizing pipeline parameters, i.e., friction factor, length, and wave speed, on computation efficiency was evaluated and compared. The dynamic memory allocation scheme not only reduces costs substantially in CPU execution time and memory space compared to other schemes but also shows significant potential as a real-time unsteady flow predictor for large pipe networks.     

Transient Modeling of Arbitrary Pipe Networks by a Laplace-Domain Admittance Matrix

An alternative to the modeling of the transient behavior of pipeline systems in the time-domain is to model these systems in the frequency-domain using Laplace transform techniques. Despite the ability of current methods to deal with many different hydraulic element types, a limitation with almost all frequency-domain methods for pipeline networks is that they are only able to deal with systems of a certain class of configuration, namely, networks not containing second-order loops. This paper addresses this limitation by utilizing graph theoretic concepts to derive a Laplace-domain network admittance matrix relating the nodal variables of pressure and demand for a network comprised of pipes, junctions, and reservoirs. The adopted framework allows complete flexibility with regard to the topological structure of a network and, as such, it provides an extremely useful general basis for modeling the frequency-domain behavior of pipe networks. Numerical examples are given for a 7- and 51-pipe network, demonstrating the utility of the method.     

Pipe Index Vector: A Method to Improve Genetic-Algorithm-Based Pipe Optimization

This paper proposes a methodology for the optimal design of water distribution systems based on genetic algorithms. The objective of the optimization is to minimize the capital cost, subject to ensuring adequate pressures at all nodes during peak demands. The proposed method is novel in that it involves the use of a pipe index vector to control the genetic algorithm search. The pipe index vector is a measure of the relative importance of pipes in a network in terms of their impact on the hydraulic performance of the network. By using the pipe index vector it is possible to exclude regions of the search space where impractical and infeasible solutions exist. By reducing the search space it is possible to generate feasible solutions more quickly and hence process much healthier populations than would be the case in a standard genetic algorithm. This results in optimal solutions being found in a fewer number of generations resulting in a substantial saving in terms of computational time. The method has been tested on several networks, including networks used for benchmark testing least cost design algorithms, and has been shown to be efficient and robust.     

Improved Pressurized Pipe Network Hydraulic Solver for Applications in Irrigation Systems

GESTAR is an advanced computational hydraulic software tool specially adapted for the design, planning, and management of pressurized irrigation networks. A summary is given of the most significant characteristics of GESTAR. The hydraulic solver for quasi-steady scenarios uses specific strategies and incorporates several new features that improve the algorithms for pipe network computation, overcoming some of the problems that arise when attempting to apply drinking water software, using the gradient method, to irrigation systems. It is shown that the gradient method is a nodal method variant, where flow rates are relaxed using head loss formula exponents. Although relaxation produces a damping effect on instabilities, it is still unable to solve some of the numerical problems common to the nodal methods. In this contribution the results of the research on computational strategies capable of dealing with low resistance elements, hydrant modelling, multiple regulation valves, numerous emitters, and pumps with complex curves are presented, obtaining accurate results even in conditions where other software fails to converge. GESTAR incorporates all these computational techniques, achieving a high convergence rate and robustness. Furthermore, GESTAR’s solver algorithm was easily adapted to incorporate inverse analysis options for optimum network control and parameter calibration. Illustrative examples are provided, documenting the improved numerical techniques and examples of GESTAR’s performance in comparison with EPANET2, a widely used gradient method-based hydraulic solver.     

Optimal Design of Pressurized Irrigation Subunit

A linear programming (LP) model is presented for optimal design of the pressurized irrigation system subunit. The objective function of the LP is to minimize the equivalent annual fixed cost of pipe network of the irrigation system and its annual operating energy cost. The hydraulic characteristics in the irrigation subunit are ensured by using the length, energy conservation, and pressure head constraints. The input data are the system layout, segment-wise cost and hydraulic gradients in all the alternative pipe diameters, and energy cost per unit head of pumping water through the pipeline network. The output data are: segment-wise lengths of different diameters, operating inlet pressure head, and equivalent annual cost of the pipeline network. The explicit optimal design is demonstrated with design examples on lateral and submain or manifold of pressurized irrigation systems. The effect of the equations for friction head loss calculation on optimization procedure is investigated through the design example for microirrigation manifold. The performance evaluation of the proposed model in comparison with the analytical methods, graphical methods, numerical solutions, and dynamic programming optimization model reveals the good performance of the proposed model. The verification of operating inlet pressure head obtained by the proposed model with accurate numerical step-by-step method suggested that it is mostly accurate.     

Simultaneous Zonation and Calibration of Pipe Network Parameters

A procedure for simultaneously zoning and calibrating pipe network parameters is proposed and applied to the determination of pipe resistance coefficients. The methodology is aimed at grouping the parameters of all the pipes into a small number of zone parameters, constrained to keep the difference between the computed and the measured water heads below a given tolerance. It is shown that, in the case of nonlooped networks, the methodology leads to a linear minimization problem where the objective function is a measure of the heterogeneity of the estimated parameters. In the case of looped networks an iterative procedure, where the linear problem is coupled with a nonlinear problem having a restricted number of decision variables, is proposed and demonstrated. Application of the procedure to a hypothetical example is shown. In a lab experiment, the model of a pipe network made by two different types of links has been calibrated using two measurement points and three different measured data sets, each of which was obtained by adjusting the valves located in the network to modify the pressure field. Comparison of the measured and the estimated resistance coefficients shows good correlation.     

Numerical Sensitivity Study of Unsteady Friction in Simple Systems with External Flows

This paper investigates the importance of unsteady friction effects when performing water hammer analyses for pipe systems with external fluxes due to demands, leaks, and other system elements. The transient energy equation for a system containing an orifice-type external flow is derived from the two-dimensional, axial momentum equation. A quasi-two-dimensional flow model is used to evaluate the relative energy contribution of total friction, unsteady friction, and the external flow, in a 1,500?m pipeline, with orifice flows ranging from steady-state flows of 2–70% of the mean pipe flow, and a Reynolds number of 600,000. It is found that for initial lateral flows larger than around 30% of the mean flow, unsteady friction effects can probably be neglected, whereas for external flows smaller than this, unsteady friction should generally be considered. Overall, the relative role of unsteady friction is found to diminish as the external flux increases, implying that unsteady friction is not critical for systems with large external flows. These results imply that unsteady friction may have a significant impact on the validity of transient leak detection techniques that have been derived assuming quasi-steady friction. To demonstrate this point, an existing transient leak detection method, originally derived under quasi-steady conditions, is tested with unsteady friction included.     

Partial Blockage Detection in Pipelines by Frequency Response Method

A new technique is presented utilizing the frequency response for the detection of partial blockages in a pipeline. In the system frequency response, a partial blockage increases the amplitude of the pressure oscillations at even harmonics. Such an increase in amplitude has an oscillatory pattern, the frequency and amplitude of which may be used to predict the location and size of a partial blockage. In this technique, the pressure transient history at only one location is sufficient, and the history of the transient in the pipe prior to blockage is not needed, which is an advantage over a number of other available techniques, in addition to being simpler to use. It is shown that the technique successfully detects the location of a blockage in a number of simple systems with blockage size as small as 10%. The technique is verified by comparing the computed results with those computed by the method of characteristics and with measurements from simple laboratory setups. A number of practical issues and limitations for field implementations are discussed.     

Explicit Power Formula for the Darcy–Weisbach Pipe Flow Equation: Application in Optimal Pipeline Design

Although the Darcy–Weisbach equation combined with the Colebrook–White semitheoretical formula for calculating the friction coefficient is a highly accurate generalized pipe-water flow resistance equation, most users prefer the use of simple, explicit power law form formulas. Because of their simplicity (despite their limitations) the purely empirical power formulas of Hazen–Williams and Manning remain the most popular pipe flow resistance equations used in routine hydraulic engineering applications. In this paper, a new simple power law form formula is derived to approximate the generalized Darcy–Weisbach combined with the Colebrook–White equation. The two main pipe flow parameters, such as the discharge (or velocity) and the diameter, appeared explicitly in the proposed formula. The suggested power-form formula compared with the Darcy–Weisbach and Coolbrook–White equation yields a maximum relative error of about ±4.5%. The power-form suggested formula is dimensionally homogeneous and its accuracy is sufficient for practical engineering applications. A correction factor is introduced for the variation of kinematic viscosity with temperature. The usefulness of the formula is demonstrated in an application concerning the optimal design of a delivery pipeline with pumping. The power form of the friction formula facilitates the formulation of the problem leading to the derivation of a simple equation from which the economic diameter is explicitly calculated.     

Calibration of Water Distribution Hydraulic Models Using a Bayesian-Type Procedure

Estimating model parameters is a difficult, yet critical step in the use of water distribution system models. Most of the optimization-based approaches developed so far concentrate primarily on efficient and effective ways of obtaining optimal calibration parameter values. At the same time, very little effort has been made to determine the uncertainties (i.e., errors) associated with those values (and related model predictions). So far, this has typically been done using the first-order second moment (FOSM) method. Even though reasonably computationally efficient, the FOSM approach relies on several restrictive assumptions and requires computationally demanding calculation of derivatives. To overcome these limitations, the recently developed shuffled complex evolution metropolis (SCEM-UA) global optimization algorithm is linked to the Epanet2 hydraulic model and used to solve a least-squares-type calibration problem. The methodology is tested and verified on the Anytown literature case study. The main advantage of the SCEM-UA algorithm over existing approaches is that both calibration parameter values and associated uncertainties can be determined in a single optimization model run. In addition, no model linearity or parameter normality assumptions have to be made nor any derivatives calculated. The main drawback of the SCEM-UA methodology is that it could, potentially, be computationally demanding, although this is not envisaged as a major problem with current computers.     

高炉干式除尘喷水（碱液）降温系统

针对高炉煤气全干法除尘系统酸性气体含量高,容易造成煤气管网及其附属设备严重腐蚀的问题,首钢京唐公司采用了喷水（碱液）中和、降温方法,并加强管道保温和管道设备内的防腐等措施,对减少腐蚀取得了良好效果。     

基于IPSO-GRU深度学习算法的海底管道缺陷尺寸磁记忆定量反演模型

针对海底管道缺陷磁记忆定量反演的难题,提出一种基于改进粒子群优化的门控循环神经网络模型,即IPSO-GRU模型。以两端焊有盲板的X52管道作为实验材料,其上预制有不同直径、深度的缺陷,采用TSC-5M-32磁记忆检测仪,外接11-6W非接触探头,进行水下磁记忆检测试验,提取不同缺陷尺寸的磁记忆信号特征值。考虑到磁记忆信号特征值随缺陷尺寸呈复杂的非线性变化,引入门控循环神经网络,利用其双门结构能够记忆缺陷处的信号特征,非线性回归拟合能力强的特点,构建海底管道缺陷定量反演模型,进一步考虑到模型超参数选择的随机性,采用改进粒子群算法进行超参数寻优。验证结果表明:该模型对缺陷深度反演平均精度达96%;对缺陷直径反演平均精度达93%,为海底管道缺陷的磁记忆定量化识别与反演提供了新的思路和方法。      

Considering Actual Pipe Connections in Water Distribution Network Analysis

The classical assumption of representing total demand along a pipe as two lumped withdrawals at its terminal nodes is hitherto common. It is a simplification of the network topology which is useful in order to drastically reduce the number of nodes during network simulation. Conversely, this simplification does not preserve energy balance equation of pipes and, for this reason, it is an approximation that could generate significant head loss errors. This paper presents a modification of the global gradient algorithm (GGA) which entails an enhancing of GGA (EGGA) permitting the effective introduction of the lumped nodal demands, without forfeiting correctness of energy balance, by means of a pipe hydraulic resistance correction. The robustness and convergence properties of the algorithm are compared with those of the classical GGA. Furthermore, the effectiveness of EGGA is demonstrated by computing the network pressure status under different configurations of the connections along the pipes of a test network.     

Modeling of Pipe–Soil Interaction and Its Application in Numerical Simulation

This paper presents three plasticity models that can be applied to numerically simulate pipe–soil interaction. They can be applied individually to evaluate the force–displacement response of a small plane-strain pipe section or in combination to simulate a long pipeline system. In the latter, numerous pipe–soil elements are attached to structural finite elements, each simulating localized foundation restraint along the pipeline. The three models are increasing in sophistication, mainly due to the manner in which they account for the behavior within an allowable combined loading surface. The first is based on traditional strain-hardening plasticity theory and therefore assumes a purely elastic response inside a single expandable yield surface. The second allows some plasticity due to the use of a bounding surface, and the third accounts for kinematic hardening through the introduction of a second smaller surface. The models are detailed in this paper, allowing for simple numerical implementation. Importantly, they are incorporated within the structural analysis of a pipeline and their potential to investigate generic pipeline system behavior is demonstrated. The applicability of the three models is interpreted theoretically and their differences shown through application for (1) a one pipe–soil interaction element and along (2) a 100?m segment of pipeline. The latter shows the practical application of these models to offshore pipeline engineering examples, with the influence of a free span behavior investigated. The ability to model complex cyclic loading is also shown.     

高压水除鳞系统管路预充水方式的选择

热轧生产线配套的高压水除鳞系统设备在除鳞作业间歇时对除鳞阀后管路进行预充水,能有效防止管路失水变成空管而集气,消除除鳞作业时引起管路冲击、振动,防止管路设备的损坏,其预填充水方式有低压充水和高压充水两种。非节能型除鳞泵和蓄能器混合供水的除鳞系统,以及除鳞泵直接供水的除鳞系统宜选择高压充水方式,系统增加的能量很少,安全可靠;节能型除鳞泵和蓄能器混合供水的除鳞系统宜选择低压充水方式,系统能耗低,但系统安全性相对较差。     

Optimal Design of Level 1 Redundant Water Distribution Networks Considering Nodal Storage

A water distribution network (WDN) is designed to meet time-varying demands with sufficient pressure, taking into consideration an appropriate demand during peak hours. Therefore, a network has inherent redundancy in the sense that under abnormal conditions such as those arising due to pipe breaks or pump failures, deficiency in supply during peak hours can be met through additional supply during off-peak periods. However, this necessitates a storage facility at the consumer end of the network, which is normally available in the form of a sump or an overhead tank in developing countries. Such a storage enables the consumer to store water during the off-peak period and then use it during the peak period. Reliability of a WDN is assessed herein considering nodal storage, and an iterative method is proposed for the optimal design of Level 1 redundant WDNs, i.e., networks that can sustain a single pipe failure without affecting consumer services either in part or in full. The method is illustrated through an example and the designs of a network with and without storage are compared. Provision of a nodal storage is found to reduce the total cost of the network.     

Calibration of On-Demand Irrigation Network Models

In this study, a new procedure for calibrating on-demand irrigation network models was developed. This procedure used a new objective function called maximum data with a reasonable error (MDRE) for calibrating the network. It was compared with the two more commonly used objective functions in calibration procedures that are the simple least squares (SLS) and the maximum likelihood estimator for the heteroscedastic error case (HMLE). In order to carry out the calibration, a quasi-Newton optimization method was used having as variable the Hazen-Williams head losses coefficient (C). This procedure was applied to an on-demand irrigation network located in Tarazona de La Mancha (Albacete, Spain) where flow and pressure at hydrant level was measured. The calibration procedure using the MDRE objective function was applied considering all the pressure control points simultaneously and the obtained results were compared with the results of considering the pressure control points independently. Therefore, the effect of the location of the pressure control point was studied. Results showed that, when the proposed objective function was used, the root mean squared error (RMSE) comparing the measured and simulated data after calibration was lower than when the SLS or HMLE objective functions were used. The location of the pressure control points throughout the irrigation network could affect the results; therefore, it was more accurate to use all the control points simultaneously than independently in the calibration process.     

Frequency-Domain Modeling of Transients in Pipe Networks with Compound Nodes Using a Laplace-Domain Admittance Matrix

An alternative to modeling the transient behavior of pipeline systems in the time domain is to model these systems in the frequency domain using Laplace transform techniques. A limitation with traditional frequency-domain pipeline models is that they are only able to deal with systems of a limited class of configuration. Despite the development of a number of recent Laplace-domain network models for arbitrarily configured systems, the current formulations are designed for systems comprised only of pipes and simple node types such as reservoirs and junctions. This paper presents a significant generalization of existing network models by proposing a framework that allows not only complete flexibility with regard to the topological structure of a network, but also, encompasses nodes with dynamic components of a more general class (such as air vessels, valves, and capacitance elements). This generalization is achieved through a novel decomposition of the nodal dynamics for inclusion into a Laplace-domain network admittance matrix. A symbolic example is given demonstrating the development of the network admittance matrix and numerical examples are given comparing the proposed method to the method of characteristics for 11-pipe and 51-pipe networks.     

基于API579准则对弯管腐蚀数值模拟的分析

为定量描述高压回注管道弯管处由于流体速度与方向的改变而引起的物理冲刷腐蚀,本文以川西北气矿的高压回注水管道为例,利用Fluent软件建立了垂直弯管冲刷腐蚀模型,并利用API579准则对相应弯管处的剩余强度进行评价和研究。结果表明,腐蚀速率最大区域发生在弯管段外壁。同时,根据数值模拟计算对各控制参数根据遗传算法进行优化处理。针对管道内壁结垢无法测量准确数据的情况,建立了一套相对完整的评价优化分析系统。