One of the key factors in the operation of a natural gas pipeline network is the linepack in the network. The desired operation of the network as derived from estimated receipts and deliveries is expressed in terms of the desired linepack profile that must be maintained. The compressor stations in the pipeline network are then operated in a manner that generates this linepack profile. Generally, the operating points selected for the units in the compressor stations are based on experience and experimentation and are therefore not optimal. In this paper, we present a systematic approach for operating the units of a compressor station to meet a specified linepack profile. The first step in developing this approach is the derivation of a numerical method for analyzing the flow through the pipeline under transient nonisothermal conditions. We have developed and verified a fully implicit finite difference formulation that provides this analysis capability. Next, the optimization of the compressor stations is formulated as a standard nonlinear programing problem in the following form: Find the values in the design variable vector denoted by $b=[b1,b2,…,bn]T$, to minimize a given objective function $F(b)$, subject to the constraints $gj(b)⩽0$, $j=1,…,m$. Here, $n$ is the number of operational parameters whose optimal value is to be determined, while $m$ is the number of operational constraints that must be enforced. In our formulation, the design variables are chosen to be the operating speeds of the units in the compressor stations, while the objective function is taken to be the average fuel consumption rate over the interval of interest, summed over all units. The constraint functions $gj(b)$ are formulated suitably to ensure that operational limits are met at the final solution that is obtained. The optimization problem is then solved using a sequential unconstrained minimization technique (SUMT), in conjunction with a directed grid search method for solving the unconstrained subproblems that are encountered in the SUMT formulation. The evaluation of the objective function and constraint functions at each step of the optimization is done by using the fully implicit analysis method mentioned above. A representative numerical example has been solved by the proposed approach. The results obtained indicate that the method is very effective in finding operating points that are optimal with respect to fuel consumption. The optimization can be done at the level of a single unit, a single compressor station, a set of compressor stations, or an entire network. It should also be noted that the proposed solution approach is fully automated and requires no user involvement in the solution process.

1.
Marques
,
D.
, and
Morari
,
M.
, 1988, “
On-Line Optimization of Gas Pipeline Networks
,”
Automatica
0005-1098,
24
(
4
), pp.
455
469
.
2.
,
A. J.
, 1994, “
Dynamic Optimization of High Pressure Gas Networks Using Hierarchical Systems Theory
,”
26th Annual Meeting Pipeline Simulation Interest Group (PSIG)
,
San Diego, CA
, Oct. 13–14.
3.
Furey
,
B. P.
, 1993, “
A Sequential Quadratic Programming-Based Algorithm for Optimization of Gas Networks
,”
Automatica
0005-1098,
29
(
6
), pp.
1439
1450
.
4.
Vostry
,
Z.
,
Zaworka
,
J.
, and
Jenicek
,
T.
, 1994, “
Optimum Control of Gas Transport
,”
26th Annual Meeting Pipeline Simulation Interest Group (PSIG)
,
San Diego, CA
, Oct. 13–14.
5.
Carter
,
R. G.
, 1998, “
Pipeline Optimization: Dynamic Programming After 30Years
,”
30th Annual Meeting Pipeline Simulation Interest Group (PSIG)
,
Denver, CO
, Oct. 28–30.
6.
Pietsch
,
U.
,
Rachford
,
H. H.
, and
Carter
,
R. G.
, 2001, “
Investigating Real-World Applications of Transient Optimization
,”
33rd Annual Meeting Pipeline Simulation Interest Group (PSIG)
,
Salt Lake City, UT
, Oct. 17–19.
7.
Rachford
,
H. H.
, and
Carter
,
R. G.
, 2000, “
Optimizing Pipeline Control in Transient Gas Flow
,”
32nd Annual Meeting Pipeline Simulation Interest Group (PSIG)
,
Savannah, GA
, Oct. 28–30.
8.
Carter
,
R. G.
, and
Rachford
,
H. H.
, 2003, “
Optimizing Line-Pack Management to Hedge Against Future Load Uncertainty
,”
35th Annual Meeting Pipeline Simulation Interest Group (PSIG)
,
Bern, Switzerland
, Oct. 15–17.
9.
Abbaspour
,
M.
,
Chapman
,
K. S.
, and
Keshavarz
,
A.
, 2004, “
Dynamic Modeling of Non-Isothermal Gas Pipeline Systems
,”
ASME International Pipeline Conference
, Oct. 4–8, Paper No. IPC04–0081.
10.
Abbaspour
,
M.
,
Krishnaswami
,
P.
, and
Chapman
,
K. S.
, 2005, “
Non-Isothermal Compressor Station Optimization
,”
ASME J. Energy Resour. Technol.
0195-0738,
127
, pp.
131
142
.
11.
Krishnaswami
,
P.
,
Chapman
,
K. S.
, and
Abbaspour
,
M.
, 2004, “
Compressor Station Optimization for Linepack Maintenance
,”
36th Annual Meeting Pipeline Simulation Interest Group (PSIG)
,
Palm Springs, CA
, Oct. 20–22.
12.
Abbaspour
,
M.
, and
Chapman
,
K. S.
, 2003, “
Non-Isothermal Compressor Station Transient Modeling
,”
35th Annual Meeting Pipeline Simulation Interest Group (PSIG)
,
Bern, Switzerland
, Oct. 15–17.
13.
Dranchuck
,
P. M.
,
Purvis
,
R. A.
, and
Robinson
,
D. B.
, 1974, “
Computer Calculations of Natural Gas Compressibility Factors Using the Standing and Katz Correlation
,”
Institute of Petroleum Technical Series
,
IP74-008
, pp.
1
13
.
14.
Odom
,
F. M.
, 1990, “
Tutorial on Modeling of Gas Turbine Driven Centrifugal Compressors
,”
22nd Annual Meeting Pipeline Simulation Interest Group (PSIG)
,
Baltimore, MA
, Oct. 18–19.
You do not currently have access to this content.