3.1. Oscillatory Process of the Pressure Function in Space and Time
Neglecting the influence of gravitational and coriolis forces, we reduce system (3) to the following equation:
where
2a is the linearization coefficient:
This equation describes the oscillatory process of the pressure function in space and time, and is known as telegraphic equation in mathematical physics. We note that pressure fluctuations in the gas flow may have different frequencies and amplitudes depending on their cause. Given the above, pressure fluctuations are conventionally divided into high-frequency, medium-frequency, and low-frequency ones. High-frequency fluctuations are characterized by frequencies in the range of 0.4–4.0 Hz, and usually result from a jump-like change in the parameter (pressure, flow) in a certain section of the gas pipeline. The amplitude of such fluctuations can reach 1 MPa. Fluctuations propagate along the gas pipeline with the speed of sound, while the amplitude and frequency decrease. The medium-frequency range is 0.5–1.0 Hz. Such fluctuations cause smooth changes in the flow parameters over time. They propagate along the pipeline with a significantly lower decrement of damping. Low-frequency fluctuations are caused by daily unevenness of gas consumption and are in the frequency range of 10
−5–0.5 Hz. The amplitude of the pressure variation depends on the nature of the disturbance factor and may be unrestricted (for example, under conditions of filling the gas pipeline with gas). In conditions of high frequency fluctuations, the inertial forces and hydraulic resistance forces in the flow of gas play a decisive role in the formation of the process. Forces of the hydraulic resistance of the pipeline are the main source of the medium- and low-frequency fluctuations. As regards the reliable operation of the gas transmission system, a crucial role is played by high-frequency fluctuations of pressure, given the unpredictable nature of the process. Since the frequency and amplitude of pressure fluctuations caused by disturbance of the gas flow parameters are characteristics of the unsteady process, there must be a link between the amplitude-frequency characteristics and the criterion of unsteadiness [
12,
13,
14].
Thus, there is an optimization problem, which consists in determining the rational pressures of the steady process in the gas pipeline, under which minimum energy costs on transport would be achieved, on the one hand, and safe operation of pipeline systems would be ensured, on the other hand [
15]. As noted, the maximum possible pressure in the gas pipeline allows minimizing the hydraulic losses during the transportation of gas, that is, achieving the minimum energy costs. However, pressure fluctuations in unsteady processes caused by a jump-like change in parameters (most often gas flow under insufficient loading) may go beyond the limits of permissible loads. Therefore, it is necessary to choose the maximum possible pressure of the steady process in the gas pipelines, at which the superposition of the pressure amplitude in the unsteady process would not force the load of pipe walls beyond the range of permissible values [
16].
The above statement of the problem requires the solution of Equation (4) under the following initial and boundary conditions chosen for the following considerations. Prior to the beginning of an unsteady process caused by disturbance of gas flow, the gas pipeline worked in steady process conditions with the lengthwise distribution of pressure according to a parabolic law.
where
is the pressure at distance
from the beginning of the gas pipeline with length
;
is the pressure at the beginning and at the end of the gas pipeline, respectively.
At given pressures
, a certain mass productivity
of the gas pipeline is provided, which can be increased or decreased to a certain magnitude at any time under insufficient loading. We assume that starting from the moment
, the supply of gas to the gas pipeline has not changed, whereas the offtake of gas at the end of the route has changed by
. Then the boundary conditions for the realization of Equation (2) will take the following form:
where
. Align correctly
Using the first equation of system (3) and neglecting all types of power consumption except hydraulic resistance, we obtain:
where
is the cross-sectional area of the pipeline.
Equation (4) under the initial (5) and boundary (7) conditions is solved by the Fourier method:
3.2. Modeling of Unsteady Pressure Fluctuations in the Gas Pipeline with Enabled En-Route Offtake of Gas
The amplitudes of pressure fluctuations were estimated in time and along the route of the gas pipeline, when disturbances in the form of jump-like changes in productivity occurred at the beginning or at the end of the gas pipeline. We considered pressure fluctuations in the initial section
P(0,
t), where the greatest values of the absolute pressure in superposition with amplitude fluctuations can exceed the permissible load. Pressure fluctuations depend on a jump-like change in the gas flow as a disturbance factor, as well as on the coordinates of disturbance, absolute values of pressure and temperature, physical properties of gas [
12,
13,
14]. They were calculated from Equation (8) with different values of the above parameters as model parameters.
It is known that when temperature and other basic physical properties of gas change in the ranges that correspond to the real operating conditions of gas pipelines, their influence on the amplitude of pressure fluctuations is insignificant. The main parameters that determine the amplitude and frequency of pressure fluctuations in the unsteady process caused by a jump-like change in the gas flow are the working pressure, the value of the gas flow, and the linear coordinate of the gas offtake. Unsteady processes at different values of working pressure and linear coordinate of the gas offtake are simulated in
Figure 1.
Let us denote the physical properties of gas by the gas constant R, and the thermal regime by the absolute temperature T. In this case, given the results of the model application, it can be stated that the product RT is a characteristic parameter of the dynamics of the unsteady process, which determines (although not to a large extent) the amplitude and frequency of pressure fluctuations. In other words, when changing the gas constant and the temperature of gas in such a way that their product remains unchanged, the nature of pressure fluctuations in the unsteady process will be the same. On the other hand, the RT product characterizes the propagation rate of small disturbances of gas, that is, the speed of sound . Hence, in order to characterize the unsteady process, it is necessary to choose the propagation rate of sound in the gas flow as the independent parameter.
The analysis of results presented in
Figure 1 shows that during a jump-like change in the gas flow, the pressure wave propagates with the speed of sound to the initial cross section of the gas pipeline where the pressure is maximal, and causes the oscillatory process, the amplitude and frequency of which increase gradually and reach the maximum after (4–6) periods of oscillation. After that the amplitude and frequency of the process begin to decrease. The total duration of the oscillation process with high values of amplitude and frequency is 50–60 seconds, which is considered as a short-term overload of the pipeline. Subsequently, the amplitude and frequency of pressure fluctuations decrease significantly, and within 0.5–0.583 hours the oscillation process fades away completely, transferring the operation of the gas pipeline to a new steady regime.
A number of calculations made in accordance with the proposed mathematical model allowed establishing a number of regularities in the oscillation process caused by a jump-like change in the gas flow [
16]. In particular, it was established that with a decrease in the working pressure in the gas pipeline, the total duration of the high-frequency oscillation process and the unsteady process in general increases, and the frequency and amplitude of oscillations decrease. Thus, when the working pressure is reduced from 7.5 MPa to 7.0 MPa (6.7%), the maximum amplitude of pressure fluctuations decreases from 0.199 MPa to 0.18 MPa, that is, by 9.5%, and when the working pressure is reduced to 6.5 MPa (by 13.3%), the amplitude of pressure fluctuations is 34.7%. The maximum frequency of the oscillation process at a working pressure of 7.5 MPa is 0.44 Hz, and with a decrease in pressure to 7.0 MPa it decreases by 21.8%, and with a further decrease in pressure to 6.5 MPa it decreases by 39.4%. From the physical point of view, this is due to a decrease in the elasticity of the medium, in which oscillating waves propagate, leading to a decrease in the propagation velocity of disturbances and, consequently, to an increase in the duration of the unsteady process and its high-frequency band [
16,
17].
The linear coordinate of the localized gas offtake has a significant impact on the nature of the unsteady process, in particular, on the amplitude and frequency of pressure fluctuations. When the localized offtake is close to the initial cross-section of the linear region, where the constant pressure is maximal, the amplitude of pressure fluctuations increases, and the frequency decreases. When the localized offtake is in the middle of the linear section of the gas pipeline, the maximum amplitude of the pressure fluctuation is 18.3% lower than that of the localized offtake located in the initial section of the gas pipeline. If the localized offtake is transferred to the final part of the linear section, the amplitude decreases by 48.5% under other identical conditions, and the frequency of oscillations varies in a smaller range. So, in the first case (when the localized offtake is in the middle of the linear section), the frequency increases by 7.2%, while in the second case (when the localized offtake is at the end of the linear section) – by 11.4%.
A change in the speed of sound has a lesser influence on the nature of pressure fluctuations in the unsteady process caused by a jump-like change in the gas flow. In addition, an increase in the speed of sound leads to an increase in the amplitude and frequency of oscillations. With an increase in the speed of sound from 400 m/s to 440 m/s (10%), the amplitude increases by 5.7%, and the frequency is 3.1%. With an increase in the speed of sound to 480 m/s (by 20%), the amplitude increases by 8.5%, and frequencies – by 4.9%.
The value of the localized gas offtake has the greatest influence on the amplitude and frequency of pressure fluctuations at the beginning of the linear section of the gas pipeline in the unsteady process caused by a jump-like change in gas flow. If a jump-like gas flow caused by the localized offtake is 10% of the total gas flow in the gas pipeline under steady conditions, the maximum amplitude of pressure fluctuations in the unsteady process caused by a sudden leakage will be 0.154 MPa. With an increase in the gas flow caused by the localized offtake to 20%, the amplitude of pressure fluctuations increases to 0.287 MPa, that is, by 2.45 times. If the gas flow due to the localized offtake is 30% of the total gas flow in the gas pipeline under steady conditions, the amplitude of pressure fluctuations will be 0.517 MPa, that is, it will increase by 3.55 times, while at a jump-like increase in the gas flow due to the localized offtake to 50% of the gas flow in the gas pipeline under steady conditions, the amplitude of pressure fluctuations will be 1.14 MPa, which may endanger the safe operation of the gas pipeline due to a short-term overload.
The maximum frequency of pressure fluctuations in an unsteady process decreases with an increase in the gas flow caused by the localized offtake. With an increase in the gas flow caused by the localized offtake from 10% to 20%, the maximum frequency of pressure fluctuations decreases by 5.8%, and with a further increase in the gas flow caused by the localized offtake to 30%, the maximum frequency of pressure fluctuations decreases by 12.3%.
3.3. Ability to Regulate Pressure by Disconnecting Individual Compressor Stations
The analysis of the results of simulation of an unsteady process in a gas pipeline caused by a jump-like change in the gas flow under conditions of insufficient loading shows that with the maximum allowable constant pressure at the beginning of the linear section of the gas pipeline, the amplitude of the pressure fluctuation can lead to a short-term overload of the pipe walls, that is, to an unsafe operation of the gas transmission system. Therefore, when controlling operating conditions, it is necessary to make decisions on providing the permissible pressure at the beginning of the linear section (at the outlet of the compressor station) in cases of a jump-like change in the gas flow. One of the options is disconnecting individual compressor stations. As shown in [
18,
19,
20], depending on the quantity and serial numbers of the working stations, it is possible to achieve the required performance of the gas transmission system. At the same time, given the low efficiency of gas-pumping units with gas turbine drives, disconnecting individual CSs can be the most effective method of controlling performance from the energetic point of view. It is obvious that this method of control can be used for seasonal regulation of performance, and it should be borne in mind that shutdown and re-start of the CS will require additional energy costs.
However, from the technological point of view, shutdown and re-start of the compressor station may lead to an unsteady process, the duration of which should be predicted in order to provide consumers with gas. It is necessary to take into account the influence of unsteady processes in long-distance gas transmission systems with a large number of compressor stations [
12,
13,
14]. Prediction and analysis of the said conditions in the gas transmission system, as well as the estimation of energy losses, are possible using the mathematical model of the unsteady gas flow in pipes taking into account an increase in pressure at compressor stations and discontinuity of flow:
where
P(
x,
t) is the pressure in the gas pipeline as a function of the linear coordinate
x and time
t;
is an increase in pressure at the compressor station with coordinate
;
) is the Dirac source function, which simulates an increase in pressure at the compressor station,
is the density of gas;
w is the linear velocity of gas;
d is the internal diameter of the gas pipeline;
is the coefficient of hydraulic resistance.
We note that in order to simulate the fading vibration processes in the gas pipeline, the equation of gas flow should include inertial hydraulic losses and frictional losses. The given system of differential Equations (9) is reduced to the equation:
where
) is the linear derivative of the Dirac function from the linear coordinate;
is the speed of sound in gas;
We assume that a gas transmission system with the length
L has
m intermediate compressor stations, which start to work simultaneously at the point of time
t = 0, and that station
k is disconnected at the point of time
. In this case, Equation (10) will take the following form:
where
is the single Heaviside function. We assume that at the initial moment of time, the gas pipeline was stopped, and constant pressure
Po was maintained throughout its length. Then the initial conditions will be:
Starting from the moment , constant pressure is maintained at the beginning of the gas pipeline, and constant pressure is maintained at its end.
For the solution of the mathematical model, integral transformations were used, in particular: Fourier sine-transformation and Laplace transformation [
10,
11]
The application of inverse transformations of Laplace and Fourier after simple transformations allows us to obtain the dependence of the pressure variation on the length and time of the unsteady process in the following form:
The first four components of the solution of Equation (12) characterize steady operating conditions of the gas transmission system. The fifth component describes the unsteady process caused by the simultaneous connection of all the compressor stations at the moment
t = 0. The last component modulates the unsteady process caused by the disconnection of the
k-th compressor station from the moment
. If
, that is, the process that occurs within the pipeline is considered after a significant period of time from the moment of connecting all CSs, the initial unsteadiness will not have a significant impact due to a higher order of smallness of the
multiplier, and the solution to the problem of disconnecting the
k-th compressor station can be represented as:
The solution of Equation (13) describes an unsteady process caused by disconnecting the
k-th compressor station and does not take into account the unsteadiness of the initial process of connecting all the CSs. Therefore, the countdown can start from the moment of disconnecting the
k-th compressor station. In this case, we obtain:
Equation (14) allows predicting the nature of the unsteady process in long-distance gas transmission systems with a large number of compressor stations caused by the shutdown and re-start of one of the stations.
To estimate the duration of the unsteady process, it is necessary to construct the dependence of fluctuations over the period of the mass flow of gas as the most inertial characteristic in the initial or final section of the gas pipeline [
16,
19,
20].
To this end, we use the equation of the gas flow from system (9). Obviously, for the initial (
x = 0) or finite cross section (
x = L), the delta-function of Dirac is
, therefore:
To simplify the computational process, we neglect inertial losses in the initial and final cross-sections, that is, we accept that
= 0 [
3]. This, of course, is associated with a certain error in the calculation of the mass flow of gas, however, in the forecast calculations, it is not important to determine the absolute value of gas consumption, but the dynamics of its change over time. In addition, by using the linearization of the equation of the gas flow, we obtain:
,
Using Equation (14) after differentiation we obtain:
The obtained dependencies,
Figure 2, allow predicting the nature of the mass flow fluctuation over time at the beginning and end of a large-distance gas transmission system, which employs
m intermediate compressor stations, caused by shutdown or re-start of the
k-th compressor station (
k = 1,2, …,
m).
It is known that oscillatory processes in a gas pipeline create loading conditions different from normal operating condition, which may cause dynamic non-equilibrium processes in the pipe material [
21], unpredictable accumulation of local micro-strains, and changes in the crack resistance of pipe steels. Therefore, minimizing such negative phenomena is not only a matter of optimizing the gas pumping technology, but also ensuring the integrity of the main gas pipeline.