Regular Article
Application of Particle Swarm Optimization (PSO) algorithm for Black Powder (BP) source identification in gas pipeline network based on 1D model
Khalifa University, Sas AlNakhl Campus, Abu Dhabi, United Arab Emirates
^{*} Corresponding author: jing.shi5@student.unsw.edu.au
Received:
20
August
2018
Accepted:
18
March
2019
Black Powder (BP) is a worldwide challenge that spans all stages of the natural gas industry from the producing wells to the consuming points. It can endanger the pipeline operations, damage instruments and contaminate customer supplies. The formation of BP inside natural gas pipeline mainly results from the corrosion of internal walls of the pipeline, which is a complex chemical reaction. This work aims to develop a novel algorithm for BP source identification within gas pipelines network based on a 1D model of BP transport and deposition. The optimization algorithm for BP source identification is developed based on the wellknown Particle Swarm Optimization (PSO) algorithm, which can solve constrained optimization problems. By applying this optimization algorithm on the gas transmission pipeline network, the BP source at different junctions could be identified and quantified simultaneously. Extensive simulation studies are conducted to validate the effectivity of the optimization algorithm.
© J. Shi et al., published by IFP Energies nouvelles, 2019
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
With the wide use of natural gas, it is challenging to maintain operational efficiency and safety for gas transmission pipeline network, which is a complex system with pipeline length varying from hundreds to thousands of kilometers (Banda et al., 2006; Tobin and Shambaugh, 2006). Black Powder (BP), occurring in both liquid and gas pipelines, is the name given to the black particulates and sediment found in pipelines. It is mainly generated due to the chemical reaction of H_{2}S, water and iron, resulting in a mixture of fine particle corrosion product and other solids, such as sands, clays, metal or construction debris, and liquid hydrocarbons chemically incorporated with any quantity of iron sulfide, iron carbonate and iron oxide contamination (Khan and Alshehhi, 2015; Khan et al., 2015; Sherik et al., 2008; Sherik, 2008; Trifilieff and Wines, 2009). Although there exist various compositions of BP, they possess some common characteristics, for example, adsorption, high specific gravity and difficult to clean. In terms of BP generation, they could come from the gas source in the gas field and also result from movement of the upstream point with gas flow, which originates from the corrosion of internal wall of any pipeline. Specifically, corrosion could only occur in the pipeline with the flow of “wet” gas (Beavers and Thompson, 2006; Baldwin, 1998), which is natural gas mixed with a certain amount of water, providing a necessary condition for corrosion. The presence of BP in gas pipeline network threatens the safety of gas industry and also has lots of undesirable influence on the operating companies all over the world (literature reviews of Khan and Alshehhi, 2015). First of all, the gas quality will be greatly affected by intolerable concentration of solid particles. Secondly, gas transmission cost will gradually increase, resulting from raised pressure drop due to changes of internal wall roughness and diameter over a long time of operation. Thirdly, the contaminations of BP could cause compressor failures, erosion of control valve and instrument clogging, etc. (Baldwin, 1998). Finally, the cleaning/removal of BP and frequent replacement of customers’ cartridge filter elements could increase the expenses each year.
To reduce the influence of BP on gas industry, many companies attempted to manage and control the BP in gas transmission pipeline network. In general, the existing methods could be classified into two different aspects (AlQabandi et al., 2015; Cattanach et al., 2011; Khan et al., 2015; Trifilieff and Wines, 2009; Tsochatzidis and Maroulis, 2007): (1) BP removal and (2) prevention methods. The approach of removing BP from gas transmission pipeline is successful in several respects, including increasing operational safety, altering operational parameters, increasing operational efficiency and facilitating effective corrosion inspection. If BP problem is not serious (small quantities), mechanical cleaning is widely used for most of the pipelines. Normally, mechanical scrapers are used to clean pipeline wall, such as mechanical pigs are deployed into a pipeline to scrape debris from pipeline wall and remove BP. Although mechanical scrap is efficient to keep the pipeline in fairly clean condition, frequent cleaning could damage the pipeline wall by exposing fresh steel surface under “wet” gas, resulting in excessive corrosion and BP generation. In addition, mechanical cleaning is not effective with complex BP formation, which calls for a combination of chemical cleaning. It has been practiced that chemical cleaning, acting as liquid soap, could significantly improve the efficiency of BP removal combined with mechanical cleaning, such as pigging operations, which acts as brush (Trifilieff and Wines, 2009). Separators and cyclofilters are also widely installed to reduce the BP concentration, where these devices could physically knock out the BP particles in the gas stream. Then, the BP particles are collected at the bottom in a collection hub. However, this method only applies to gas stream with a high concentration of solid particles and also relatively large particle size (more than 10 microns). The removal approach is practiced as a good option to protect the downstream operations, but it has several disadvantages (Trifilieff and Wines, 2009). These methods are not single trial, but frequent application, which account for a large amount of expense. Also, these solutions cannot address the location of BP formation, that is to say, they are afterward remedy. In addition, chemical cleaning needs subsequent handling procedures, which could be costly and challenging if the chemical disposal is toxic or harmful to the environment. To compensate the above drawbacks of removal approach, gas operators have an alternative solution to prevent the occurring of BP generation, which is considered as a consequence of corrosions at the internal wall of gas transmission pipeline (Cattanach et al., 2011; Trifilieff and Wines, 2009). To prevent the occurring of corrosions, the inner wall of gas transmission pipeline is normally coated with high solids solventbased epoxy polyamine films, which is used to protect the inner surface of pipelines. This approach is practiced to be cost effective; however, it is difficult to be applied to buried pipelines. Reducing water contamination is another approach to prevent the occurring of pipeline corrosions, which is based on the philosophy of internal corrosions, are largely related to the case of “wet” gas. Namely, this approach is moisture control.
Although several removal and prevention methods have been put into practice, the location of BP generation is still unpredictable and there is no research on the identification problem of BP source until now. The main reason is the unknown information about the BP generation and where it could be generated. The developed dynamic models for BP transport and deposition are limited for gas transmission pipeline network, which in turn restricts the study of modelbased methods for BP identification. It is well known that full threedimension (3D) Computational Fluid Dynamics (CFD) simulation software is specialized for multiphase flow modelling; however, it seems impossible to simulate the 3D dynamics of large pipeline networks (e.g. 100 km), which is under studied in these works (Kharoua et al., 2015, 2017). Filali et al. (2016) proposed a 1D approach for modelling transport and deposition of BP particles in gas transmission network, where the authors discussed behaviors of particles with different diameters. Moreover, two different deposition models (Fan and Ahmadi, 1993; Wood, 1981) were compared to calculate the bed height, which was validated by Discrete Phase Model (DPM) based on CFD software.
The main contribution of this work is to develop a tool to identify the BP source in gas pipeline network, which is modeled as treeshaped gas transmission network with BP dynamics of motion and distribution along each pipe. BP source identification is formulated as a constrained optimization problem, which is solved by the Particle Swarm Optimization (PSO) techniques (Delice et al., 2017; Eberhart and Shi, 1998; Kennedy and Eberhart, 1995; Min et al., 2017). PSO is a popular stochastic optimization technique with some features and advantages compared to other optimization algorithms such as Ant Colony Optimization (ATO) and Genetic Algorithm (GA) (Wiak et al., 2008; Saravanan, 2006). These features and advantages include (1) taking real numbers as particles; (2) few parameters need to be tuned and (3) simple implementation and effective global search capability, hence it could be a good candidate. In fact, PSO technique has been widely used in oil and gas industry, for example, Wu et al. (2014) optimized the operation of trunk natural gas pipelines via PSO based algorithm, and Madoliat et al. (2017) also successfully applied PSO to the transient analysis of natural gas pipeline. The PSO algorithm is used to solve a formulated optimization problem for similar applications. In particular, the basic PSO algorithm has been improved by incorporating the inertiaadaption technique to solve a constrained nonlinear optimization problem in Wu et al. (2014). Then, it was used to solve a more complex problem in Madoliat et al. (2017), where the solution of nonlinear PDE flow equations could be obtained simultaneously.
The paper is structured as follows. Section 2 will briefly discuss the 1D approach for modelling BP particles in gas networks, which is detailed in Filali et al. (2016). A general structure of treeshaped model of gas transmission pipeline network is given. In Section 3, PSObased optimization algorithm for BP identification is presented. In Section 4, extensive simulation results and discussions of the optimization algorithm for BP source identification are presented. Finally, conclusions are given in Section 5.
2 Modelling for dynamics of BP particles in gas transmission pipeline
In this section, a schematic of treeshaped gas transmission pipeline network will be built based on a set of pipeline connection rules. In addition, the methodology of a 1D approach for modelling the dynamics of transport and deposition of BP particles in the gas pipeline network will be explained (Filali et al., 2016).
2.1 A treeshaped model of gas network
In gas industry, the natural gas pipeline network is a highly integrated transmission and distribution grid that could transport natural gas from its origin to any position of high gas consumption demand. In many cases, the natural gas produced from a well has to be transported on a very long distance to the point of use. In order to maintain effective distribution of natural gas, the gas transmission system is extensive, consisting of complex pipeline topology. However, the entire gas transmission network could be separated into several subnetwork with single gas well. This simplification also contributes to the investigation of BP identification, which will be included in following sections of this paper. In this case, the gas transmission pipeline network is properly generalized as treeshaped model (Babonneau et al., 2012; Shiono and Suzuki, 2016), which could be built based on following rules for pipe connection:

Gas transmission pipeline network without loop;

Junction connected with three pipes;

Oneway flow of gas stream.
Based on these connection rules, a schematic of treeshaped gas pipeline transmission network is given in Figure 1. There is only one source for gas supply, and each junction has three pipe connection, including two main pipes and one branch pipe. The arrow on the pipe shows the flow direction of gas stream. The customer is located at the end of each branch pipe. The sources, main pipes, branch pipes and customers are denoted as S, P, B, and M respectively.
Fig. 1 Schematic of gas transmission pipeline network. 
2.2 OneD approach of modelling BP particle dynamics
In this section, the 1D approach for modelling dynamics of BP particles in gas pipeline (Filali et al., 2016) is described. It is well known that the CFD software is popular to simulate the dynamics of multiphase flow, which could give a solution to the simulation of BP particles in the gas stream. However, the computation load of CFD is extremely heavy for gas pipeline network, which could extend to hundreds of kilometers. In this case, Filali et al. (2016) proposed a simplified 1D approach of modelling, where the flow of gas is continuously mixed with BP particles, and the behavior of particle movement is modelled based on the dusty gas assumption and the usage of analytical solutions of steady scalar advection/reaction equation.
The governing equation for gas stream mixed with solid particles is given by the following onedimensional advection–diffusionreaction equation:(1)where C is the particle concentration in the gas flow, U is the average velocity of gas flow inside the pipeline, D _{diff} is the diffusion coefficient, S is a term to describe the deposition, pickup or generation of solid particles. For fully developed turbulent flow, the axial diffusion term is negligible compared to the dominant advection term and the following expression is valid,(2)
For steady state solution, equation (1) could be simplified as:(3)where, β _{dep} is the deposition rate and it could be calculated based on two methods, for example, Wood (1981), Fan and Ahmadi (1993), and ξ _{gen} is the generation rate of the BP inside the pipeline and it is considered as an unknown parameter.
3 PSObased optimization algorithm for BP identification
In this section, the objective of BP identification is presented as a constrained optimization problem, which will be solved by a PSObased optimization algorithm. The PSO algorithm is a populationbased evolutionary search algorithm inspired by social behavior of animals such as bird flocking. Basically, each particle in the swarm has a position and velocity, with its position representing a candidate solution in the multidimensional solution space and velocity indicating moves from one position to another. A fitness function is defined to evaluate each particle until certain convergence criteria is satisfied. During the searching process, the particle with fitness value will be selected as local/global best particle. The updated equations of position and velocity are given in the literature (Marini and Walczak, 2015) as follows:(4) (5)where, x _{ij}(t) denotes the position of particle i in dimension j at time t, with j = 1, 2, 3, …, n _{ x }, and n _{ x } is the dimension of the solution space. The updating of position is calculated by adding the velocity v _{ij}(t) to current position. y _{ij}(t) is the personal best solution in dimension j for particle i, while represents the best solution found so far. c _{1} and c _{2} are positive acceleration factors used to scale the contribution of cognitive and social component. are random variables with normal distribution.
The structure of this section is outlined as follows: firstly, the objective of BP identification is formulated, which is developed based on the 1D model of BP particles transport and deposition. Secondly, the BP source could be identified and quantified by sole running of the PSObased optimization algorithm. Finally, a proof will be provided to support the uniqueness of the formulated optimization problem.
3.1 Formulation of BP identification problem
It is noted that additional unknown sources of BP generation are considered in this study. Along with that, the objective of this paper is to identify the BP sources and quantify the generated BP concentration using several measurements of BP concentration at the end node of each branch pipe, which is a point of natural gas use (client). As most of the gas transmission pipeline network is buried underground, therefore the sensor used for BP concentration measurement is only installed at each client (AbouKhousa et al., 2015). The presence of an optimization algorithm for BP source identification is developed based on a 1D model, which has several assumptions (Filali et al., 2016), such as, oneway flow of gas stream, constant velocity of gas flow at each pipe, average mixture of solid particles and gas stream, fixed size/diameter of solid particle, etc., but some additional assumptions should be also taken into consideration. The accuracy of BP source identification is largely restricted because of limited measurement of BP concentration along the pipe, which mainly results from two aspects: (1) lack of reliable measuring device that can be installed underground; (2) the “additional” BP generation from the internal walls is unpredictable, that is, the occurrence of chemical corrosion is highly random and the location is difficult to be identified. Therefore, the gas supply station is assumed in this paper as the main BP generation source, while the junctions stand for additional source, which essentially represent the amount of BP generated along its downstream. Herein, two assumptions are given for the problem formulation of BP identification:

BP source is located in junction only.

Additional BP particle generation occurs in the main pipe only.
According to Figure 1, there is one main BP particle source and several unknown additional sources located at each junction, which are the connection points for every two main pipes and one branch pipe. Let us consider the gas transmission pipeline network as the system. Measurements of BP concentration at each client are the system inputs and estimations of BP concentration at each source are the system outputs. Therefore, the following constrained optimization problem (Aguirre et al., 2007; Cagnina et al., 2008; Hu and Eberhart, 2002; Liu, 2008; Parsopoulos and Vrahatis, 2002) can be formulated with weighted sum of errors (e.g. absolute error) as its cost function.(6) where, S _{1}, S _{2}, S _{3}, …, S _{ i } represent the estimation of BP concentration at sources (S _{1} is the main source, and S _{2}, S _{3}, …, S _{ i } are the additional source). m is the number of clients and BP measurement points. M _{ k } is the measurement of BP concentration at the client k, which is measured by BP concentration sensor located at the end of each branch pipe and it is simulated using the 1D model presented in the previous section. is the estimated value of BP concentration, λ _{ k } is the weighting parameters for each system input, which could be selected as inverse of each measurement of BP concentration at corresponding client. The optimization algorithm for BP source identification is presented as a flowchart in Figure 2.
Fig. 2 (a) PSO based BP source optimization algorithm. (b) Flowchart of OneD model running. 
Fig. 3 Generalized treeshaped gas transmission pipeline network with additional BP source. 
3.2 Discussion on unique solution of BP optimization problem
Proposition 3.1 Based on the given assumptions in Section 3.1 , the solution of the formulated optimization problem is unique if and only if the number of measurement is equal to the number of unknown BP source.
The proof is referred to the Appendix A.
4 Simulation studies
In this section, the optimization algorithm for BP source identification is applied to a real gas transmission pipeline network, including 15 pipes, 16 nodes, and 8 junctions. The geometric parameters of this gas pipeline network are practical data from an existing network. The topology of this network is given in Figure 1, as well as its geometric parameters in Table 1. The remaining parts are extensive simulation studies, including BP identification algorithm applied on gas network with perfect measurements, and sensitivity studies on this algorithm. The sensitivity studies are divided into two aspects, including investigations on disturbed measurements and model parameter mismatch. Finally, some discussions are presented.
Geometric parameters and flow conditions for gas network.
4.1 Ideal situations
In this section, the optimization algorithm for BP identification is validated on the gas network given in Table 1. It is noted that, initially, the measurements of BP concentration at the end node of each branch pipe (client) are simulated based a sequence of BP sources given as follows.
By sole running of the 1D model, the measurements of BP concentration at each client could be generated according to the BP concentration at each source given in Table 2. Then these measurements shown in Table 3, will be saved and assumed as known parameters, i.e. BP concentration measurements. This is the first step for following simulation studies. The second step is the estimation of BP concentration at each source using the measurements of BP concentration. It is noticeable that the aim of BP optimization problem is the reverse process, where the BP concentration at each source will be estimated by application of the optimization algorithm.
Real value of BP concentration at each source.
Measurements of BP concentration at each client.
Simulation results of identification for BP concentration at sources are given in Figure 4. In this simulation study, the parameters of PSO are selected as: , n = 300, m = 50, c _{1} = 2, c _{2} = 2, where, λ _{ k } is the weighting parameter in equation (6), n is the quantity of particles, m is the maximum number of iteration (termination condition), and c _{1}, c _{2} are the acceleration constants. It is noticeable that the PSO parameters are tuned by trial and error, which can ensure satisfactory performance with acceptable computational load. PSO is a stochastic optimization technique, therefore the algorithm is repeated three times (blue, red and green curves in each figure) to show/compare the stability of the optimization results. The dotted line in each figure represents the real value of BP concentration at each source, which is corresponding to Table 3. In addition, the statistic results of each test are summarized in Table 4. The estimation value and absolute error are given for each test and source, and an overall average with the Standard Derivation (SD) are given for all sources. It could be concluded that the optimization algorithm for BP source identification is satisfactory and able to achieve small estimation errors with perfect measurements of BP concentration at clients.
Fig. 4 Estimations of BP concentration at each source, figures (a)–(h) represent sources 1–8. 
Simulation results of optimization algorithm under ideal situations.
4.2 Sensitivity studies
As stated in Section 4.1, the BP concentration at sources can be well estimated/reproduced by the PSO algorithm when the model is completely known and there are no uncertainties in measurements. However, the efficacy of this approach is still questionable in practice, because all the models are at best only approximations for reality due to various model errors or parametric uncertainties. This problem can be formulated as follows:(7)with θ the vector of geometric and flow conditions parameters, ϵ the measurement error, M _{ k } the real measurements, the 1D model approximation and S* the vector of BP concentration at sources which are unknown. The aim is to effectively recover S* via PSO algorithm relying on M _{ k }, regardless of the presence of uncertainties in θ and ϵ. A good approach to solve this problem is associated to robust model calibration and inverse problem theory (Tarantola, 1987; Kaipio, 2008).
As a preliminary study, three case studies are simulated in this section. Firstly, the measurement error on M _{3} and M _{5} are simulated. Secondly, the disturbances in flow measurements are considered. Thirdly, parametric uncertainties are simulated for pipeline roughness. It is noted that each case study is repeated three times with same parameters to prove the algorithm’s stability and reliability.
4.2.1 Black Powder measurement uncertainty
In this section, the optimization algorithm for BP source identification is tested under the situations of measurement uncertainty, where M _{3} and M _{5} are increased by 5% respectively. Basically, ϵ is regarded as an additive nonzero noise and no variations in θ in this scenario. The tuned parameters for the algorithm are the same with ideal situations. Some figures and statistic data are presented as preliminary results. It is noted that the problem formulation in equation (7) is a deterministic case, which will be generalized in the next stage of experimental validations and tackled by the likelihood approaches (Tarantola, 1987), where measurement errors can be modelled with a probabilistic distribution and confidence intervals will be computed for estimated BP concentration at sources within Bayesian framework.
4.2.1.1 Uncertainty in M _{3} measurement
In this case, the BP concentration measurement M _{3} is increased by 5%. The performance of identification is shown in Figure 5, and also the corresponding statistic data is presented in Table 5. It can be seen from Figures 5c and 5d, 10% and 5% estimation error for sources 3 and 4 respectively, which are expected and acceptable. Although the inaccurate measurement occurring at M _{3} affects the BP concentration estimation at sources 3 and 4, and it has negligible influence on the BP concentration estimation at sources 5, 6, and the downward sources.
Fig. 5 Estimations of BP concentration at sources for 5% increase in M _{3} BP measurements; figures (a)–(h) represent sources 1–8. 
Simulation results of optimization algorithm under situation for 5% increase in M _{3}.
4.2.1.2 Uncertainty in M _{5} measurement
In this case study, M _{5} is increased by 5% and the performance of the optimization algorithm is presented in Table 6. It is noted that the estimation errors of sources 5 and 6 are 15% and 10% respectively, which are high relatively compared with their downward sources. The estimation error of sources 7 and 8 could be negligible.
Simulation results of optimization algorithm under situation for 5% increase in M _{5}.
4.2.2 Mass flow rate measurement uncertainty
Similarly, variations in θ are simulated in this section, where the mass flow rate in pipe 3 (Q _{3}) and pipe 5 (Q _{5}) are increased by 10% and 5% respectively. The tuned parameters for PSO are same with the ideal situations. One objective of this paper is to study the robustness of the model fidelity to measurement data given parametric uncertainties. Although these simulation tests are not representative in practice where noise/uncertainties always exist in a stochastic way, they can provide some preliminary results at this stage.
4.2.2.1 Uncertainty in Q _{3} measurement
In this case study, the mass flow rate measurement of Q _{3} is increased by 10%. The simulation results are shown in Figure 6, along with its statistic data in Table 7. It can be seen in Figure 6, the variation of flow rate Q _{3} has a very small influence on the optimization, which proves the robustness of the optimization algorithm under uncertainty.
Fig. 6 Estimations of BP concentration at sources for 10% increase in Q _{3}, figures (a)–(h) represent sources 1–8. 
Simulation results of optimization algorithm under situation for 10% increase in Q _{3}.
4.2.2.2 Uncertainty in Q _{5} measurement
In this case study, the mass flow rate measurement of Q _{5} is increased by 5%. The simulation results are shown in Table 8. The estimation error is less than 2% for each source, which is acceptance and negligible for this work.
Simulation results of optimization algorithm under situation for 5% increase in Q _{5}.
4.2.3 Deposition rate uncertainty
In this case study, the pipeline roughness is increased by 5%. The parameter for the algorithm is same with the ideal case and some simulation results are shown above:
It can be seen in Figure 7 that the estimation of BP concentration for each client is always higher than the expected value, and the estimation error is shown in Table 9. This is reasonable because the increase of roughness will result in more BP deposition, in other words, more BP should be generated in order to achieve the same measured BP concentration as expected. However, the average estimation error is nearly 2%, which is still acceptable in the practical situation.
Fig. 7 Estimations of BP concentration at sources for 5% increase in roughness, figures (a)–(h) represent sources 1–8. 
Simulation results of optimization algorithm under situation for 5% increase in roughness.
5 Conclusion
In this paper, the problem of BP source identification is studied by applying a PSObased optimization algorithm, which is developed upon a 1D model of BP transport and deposition. The 1D model is a simplified approach for modelling the dynamics of BP particles in gas transmission pipeline network. A schematic of treeshaped gas network is proposed, which is generalized with a set of connection rules. Through some preliminary simulation studies, the PSObased algorithm is validated to be a useful technique for BP source identification. At this stage, the model errors and uncertainties are assumed deterministic though they are not representative in practice. The main purpose of this paper is to discuss the application of PSO techniques on BP source identification. More work will be done together with the experimental studies in the next stage, including (1) analysis of likelihood approach which has specified probabilistic distribution to model error; (2) improvement of the 1D model by considering pickup case and various particle sizes.
Acknowledgments
The authors acknowledge the Gas Processing & Materials Science Research Centre (GRC) at Khalifa University – The Petroleum Institute, AbuDhabi, in addition to Abu Dhabi Gas Industries Limited (GASCO) for funding and supporting this project.
References
 AbouKhousa M., AlDurra A., AlWahedi K. (2015) Microwave sensing system for realtime monitoring of solid contaminants in gas flows, IEEE Sens. J. 15, 9, 5296–5302. [Google Scholar]
 Aguirre A.H., Zavala A.M., Diharce E.V., Rionda S.B. (2007) COPSO: constrained Optimization via PSO algorithm, Center for Research in Mathematics (CIMAT), Guanajuato, Mexico, 30 p. Technical Report No. I0704/22022007. [Google Scholar]
 AlQabandi S., Khuraibut Y., AlInzi M. (2015) Integrated Solution to Mitigate and Prevent Black Powder Phenomena in West Kuwait Operation Facilities, Proc. Conference on SPE Kuwait Oil and Gas Show, 11–14 October, Mishref, Kuwait. SPE175291MS. [Google Scholar]
 Babonneau F., Nesterov Y., Vial J.P. (2012) Design and operations of gas transmission networks, Oper. Res. 60, 1, 34–47. [Google Scholar]
 Baldwin R.M. (1998) Technical assessment: “Black Powder” in the Gas Industry – sources, characteristics, and treatment, Gas Machinery Research Council Report TA, pp. 97–104. [Google Scholar]
 Banda M.K., Herty M., Klar A. (2006) Gas flow in pipeline networks, NHM 1, 1, 41–56. [CrossRef] [Google Scholar]
 Beavers J.A., Thompson N.G. (2006) External corrosion of oil and natural gas pipelines, ASM Handbook 13, 1015–1025. [Google Scholar]
 Cagnina L.C., Esquivel S.C., Coello C.A.C. (2008) Solving engineering optimization problems with the simple constrained particle swarm optimizer, Informatica 32, 3. [Google Scholar]
 Cattanach K., Jovancicevic V., Ramachandran S. (2011) Development of new corrosion inhibitor to prevent black powder formation using quartz crystal microbalance technique, Proc. Conference on NACE international CORROSION. No. 11093. [Google Scholar]
 Delice Y., Aydoğan E.K., Özcan U., İlkay M.S. (2017) A modified particle swarm optimization algorithm to mixedmodel twosided assembly line balancing, J. Intell. Manuf. 28, 1, 23–36. [Google Scholar]
 Eberhart R.C., Shi Y. (1998) Comparison between genetic algorithms and particle swarm optimization, Evolutionary programming VII: Proc. Conference on 7th Annual Evolutionary Programming, pp. 611–616. [Google Scholar]
 Fan F.G., Ahmadi G. (1993) A sublayer model for turbulent deposition of particles in vertical ducts with smooth and rough surfaces, J. Aerosol Sci. 24, 45–64. [Google Scholar]
 Filali A., Khezzar L., Alshehhi M., Kharoua N. (2016) A oneD approach for modeling transport and deposition of Black Powder particles in gas network, J. Nat. Gas Sci. Eng. 28, 241–253. [Google Scholar]
 Hu X., Eberhart R. (2002) Solving constrained nonlinear optimization problems with particle swarm optimization, Proc. Conference on the Sixth World Multiconference on Systemics, Cybernetics and Informatics, 14–18 July, Orlando, FL, USA , Vol. 5, pp. 203–206. [Google Scholar]
 Kaipio J. (2008) Modeling of uncertainties in statistical inverse problems, J. Phys. Confer. Ser. 135, 1, 107–118. [CrossRef] [Google Scholar]
 Kennedy J., Eberhart R. (1995) Particle Swarm Optimization, Proc. on IEEE International Conference on Neural Networks. IV, 27 November–1 December, 1995, Perth, WA, Australia, Australia, pp. 1942–1948. [Google Scholar]
 Khan T.S., Alshehhi M.S. (2015) Review of black powder in gas pipelines–An industrial perspective, J. Nat. Gas Sci. Eng. 25, 66–76. [Google Scholar]
 Khan T.S., Alshehhi M.S., Stephen S., Khezzar L. (2015) Characterization and preliminary root cause identification of black powder content in a gas transmission network – a case study, J. Nat. Gas Sci. Eng. 27, 769–775. [Google Scholar]
 Kharoua N., Alshehhi M., Khezzar L. (2015) Prediction of Black Powder distribution in junctions using the Discrete Phase Model, Powder Technol. 286, 202–211. [Google Scholar]
 Kharoua N., Alshehhi M., Khezzar L., Filali A. (2017) CFD prediction of Black Powder particles’ deposition in vertical and horizontal gas pipelines, J. Petrol. Sci. Eng. 149, 822–833. [CrossRef] [Google Scholar]
 Liu C.A. (2008) New dynamic constrained optimization PSO algorithm, Natural Computation, Proc. Conference on the 4th ICNC’08, 18–20 October 2008, Jinan, China, Vol. 7, pp. 650–653. [Google Scholar]
 Madoliat R., Khanmirza E., Pourfard A. (2017) Application of PSO and cultural algorithms for transient analysis of natural gas pipeline, J. Petrol. Sci. Eng. 149, 504–514. [CrossRef] [Google Scholar]
 Marini F., Walczak B. (2015) Particle Swarm Optimization (PSO): A tutorial, Chemometr. Intell. Lab. Syst. 149, 153–165. [Google Scholar]
 Min Y., Jiayue Z., Damin Z. (2017) Immunization strategy based on discrete particle swarm optimization algorithm in BBV network, Proc. Conference IEEE Intelligent Systems and Control, 5–6 January 2017, Coimbatore, India, pp. 208–211. [Google Scholar]
 Parsopoulos K.E., Vrahatis M.N. (2002) Particle swarm optimization method for constrained optimization problems, Intell. Technol. Theory Appl. New Trends Intell. Technol. 76, 1, 214–220. [Google Scholar]
 Saravanan R. (2006) Manufacturing optimization through intelligent techniques, CRC Press, FL, USA. [Google Scholar]
 Sherik A.M., Zaidi S.R., Tuzan E.V. (2008) Black powder in gas transmission systems, Proc. Conference on CORROSION, 16–20 March, New Orleans, LA, pp. 16–20. [Google Scholar]
 Sherik A.M. (2008) Black powder: study examines sources, makeup in dry gas system, Oil Gas J. 106, 30, 54–59. [Google Scholar]
 Shiono N., Suzuki H. (2016) Optimal pipesizing problem of treeshaped gas distribution networks, Eur. J. Oper. Res. 252, 2, 550–560. [Google Scholar]
 Tarantola A. (1987) Method for data fitting and model parameter estimation. Inversion Problem Theory, Elsevier Science, New York, NY. [Google Scholar]
 Tobin J., Shambaugh P. (2006) The crucial link between natural gas production and its transportation to market, Proc. Conference on Stages in the production of pipelinequality natural gas and NGLs. Energy Information Administration, Natural Gas Annual. pp. 11. [Google Scholar]
 Trifilieff O., Wines T.H. (2009) Black powder removal from transmission pipelines: Diagnostics and solutions, Proc. Conference on Pipeline Rehabilitation and Maintenance, January 19–21, 2009, Bahrain, pp. 19–21. [Google Scholar]
 Tsochatzidis N.A., Maroulis K.E. (2007) Methods help remove black powder from gas pipelines, Oil Gas J. 105, 10, 52. [Google Scholar]
 Wiak S., Krawczyk A., Dolezel I., (eds), (2008) Intelligent computer techniques in applied electromagnetics, Springer, Berlin, Heidelberg. [CrossRef] [Google Scholar]
 Wood N.B. (1981) A simple method for the calculation of turbulent deposition to smooth and rough surfaces, J. Aerosol Sci. 12, 3, 275–290. [Google Scholar]
 Wu X., Li C., Jia W., He Y. (2014) Optimal operation of trunk natural gas pipelines via an inertiaadaptive particle swarm optimization algorithm, J. Nat. Gas Sci. Eng. 21, 10–18. [Google Scholar]
Appendix A
Proof for the Proposition 3.1
1) Sufficiency: To prove the uniqueness of the optimal solution, let us consider the general structure of treeshaped gas pipeline network, as shown in Figure 3. The network consists of N pipes, N + 1 nodes and (N − 1)/2 junctions.
According to the 1D model, the concentration at the end node of each branch pipe is given as a function of BP sources, where we have unknown BP source to be identified. The following equations can be easily found if we split the influence of each source on each measurement.(A1.a) (A1.b) (A1.c)
The compact form of above equations is(A2)
It is noted that the functions involved must be invertible, which are determined by the physical properties of the gas transmission pipelines. Although there is not explicit mathematical function describing the mapping between S and M, this function must be invertible at a particular time instant because the physical properties of pipeline do not change. Let us start with the first source S _{1} in equation (A1.a), the concentration could be determined by the first measurement M _{1} only. Based on equation (A1.b), the second measurement M _{2} is represented as a function of (S _{1}) and (S _{2}). Then, f _{21} (S _{1}) can be obtained according to S _{1}, which has been calculated from equation (A1.a). Therefore, , the remaining term of M _{2}, could be used to determine the concentration at source S _{2}. Consequently, the remaining sources can be determined step by step.
2) Necessity: If the number of client node measurements is less than the number of unknown BP source, for example, () is not available, i.e. is not measured due to fault.
Let us consider equation M _{ j+1}:(A3)where, represents the influence of source S _{1} to S _{ j−1} on M _{ j+1}. f _{ j+1, j }(S _{ j }) and f _{ j+1, j+1}(S _{ j+1}) represent the influence of S _{ j } and S _{ j+1} on M _{ j+1} respectively. Therefore, S _{ j } is a preliminary condition to identify S _{ j+1}. However, S _{ j } cannot be identified because M _{ j } is not available. In other word, unique solution of S _{ j } and S _{ j+1} cannot be achieved in this situation.
All Tables
Simulation results of optimization algorithm under situation for 5% increase in M _{3}.
Simulation results of optimization algorithm under situation for 5% increase in M _{5}.
Simulation results of optimization algorithm under situation for 10% increase in Q _{3}.
Simulation results of optimization algorithm under situation for 5% increase in Q _{5}.
Simulation results of optimization algorithm under situation for 5% increase in roughness.
All Figures
Fig. 1 Schematic of gas transmission pipeline network. 

In the text 
Fig. 2 (a) PSO based BP source optimization algorithm. (b) Flowchart of OneD model running. 

In the text 
Fig. 3 Generalized treeshaped gas transmission pipeline network with additional BP source. 

In the text 
Fig. 4 Estimations of BP concentration at each source, figures (a)–(h) represent sources 1–8. 

In the text 
Fig. 5 Estimations of BP concentration at sources for 5% increase in M _{3} BP measurements; figures (a)–(h) represent sources 1–8. 

In the text 
Fig. 6 Estimations of BP concentration at sources for 10% increase in Q _{3}, figures (a)–(h) represent sources 1–8. 

In the text 
Fig. 7 Estimations of BP concentration at sources for 5% increase in roughness, figures (a)–(h) represent sources 1–8. 

In the text 