Disclosure of Invention
The invention aims to solve the technical problem of providing a pipeline compressor energy consumption measuring method, a system and a computer storage medium, which can accurately calculate compression work in the whole changeable process of the pipeline compressor, and further can provide reliable basis for the optimized operation of the pipeline compressor.
The technical scheme for solving the technical problems is as follows: a pipeline compressor energy consumption measurement method comprises the following steps,
s1, dividing the whole variable compression process of a pipeline compressor into a plurality of variable compression sub-processes according to inlet and outlet pressures of the pipeline compressor;
s2, sequentially solving the variable compression sub-work of each variable compression sub-process based on a variable power solving model;
and S3, superposing the variable compression sub-work of all the variable compression sub-processes to obtain the variable compression work of the whole variable compression process of the pipeline compressor.
On the basis of the technical scheme, the invention can be improved as follows.
Further, the S2 is specifically that,
s21, based on the multiple transformation solving model and the set condition M, according to the known inlet condition K of the pipeline compressor 0 Calculating a variable compression sub-process G 1 Outlet condition K of (2) 1 And variable compression sub-work W pol(1) ;
S22, variable compression sub-process G x Outlet condition K of (2) x As a variable compression sub-process G x+1 Inlet condition F of (2) x+1 The method comprises the steps of carrying out a first treatment on the surface of the In the variable compression sub-process G x+1 Based on the variable power solving model and the set condition M, according to the variable compression subprocess G x+1 Inlet condition F of (2) x+1 Calculating a variable compression sub-process G x+1 Outlet condition K of (2) x+1 And variable compression sub-work W pol(x+1) The method comprises the steps of carrying out a first treatment on the surface of the Where x=1, 2,3. N-1, n is a pipeline total number of variable compression sub-processes of the compressor.
Further, the known inlet condition K 0 Including the inlet pressure P1 and inlet temperature T1 of the pipeline compressor; wherein the known inlet condition K 0 Is a variable compression sub-process G 1 Inlet conditions of (2);
the setting condition M comprises a preset pressure ratio Alfa, a preset iteration termination criterion epsilon and a preset temperature increment alpha;
variable compression sub-process G x Outlet condition K of (2) x Including the outlet pressure P (x)j Outlet temperature T (x)j Specific enthalpy of outlet h (x)j And the outlet specific entropy s (x)j ;
Variable compression sub-process G x+1 Inlet condition F of (2) x+1 Including the inlet pressure P (x+1)i Inlet temperature T (x+1)i Specific enthalpy of inlet h (x+1)i And inlet specific entropy s (x+1)i The method comprises the steps of carrying out a first treatment on the surface of the And P is (x+1)i =P (x)j ,T (x+1)i =T (x)j ,h (x+1)i =h (x)j ,s (x+1)i =s (x)j ;
The calculation termination condition of the S2 is a variable compression subprocess G n Is out of (a)Oral Condition K n Outlet pressure P in (a) (n)j Greater than or equal to the known outlet condition of the pipeline compressor, wherein the known outlet condition of the pipeline compressor is the outlet pressure P2 of the pipeline compressor.
Further, the variable power solving model comprises a Mallin model and a variable process compression power definition model;
the mullen model is specifically described as,
wherein W is
p ′
ol(x) Is a variable compression sub-process G
x Calculated value of variable compression sub-work, h
(x)i 、s
(x)i And T
(x)i Respectively the multiple compression subprocess G
x Inlet specific enthalpy, inlet specific entropy and inlet temperature, and h
(x)i =h
(x-1)j ,s
(x)i =s
(x-1)j ,T
(x)i =T
(x-1)j ;
The variable process compression work definition model is specifically W p ″ ol(x) =η pol (h (x)j -h (x)i ) Wherein W is p ″ ol(x) Is a variable compression sub-process G x Is defined by the variable compression sub-work values, eta pol Is the variable compression efficiency of the pipeline compressor in the whole variable compression process.
Further, let the variable compression sub-process G 1 Inlet pressure P of (2) (1)i Let variable compression sub-process G =p1 1 Is at an inlet temperature T (1)i =T1;
The step S21 is specifically that,
s211 at inlet pressure P (1)i And inlet temperature T (1)i Under the condition of calling NIST database, calculating variable compression subprocess G 1 Is the inlet specific enthalpy h of (2) (1)i And inlet specific entropy s (1)i ;
S212, according to inlet pressure P (1)i And a preset pressure ratio Alfa, by the formula P (1)j =P (1)i * Alfa, calculating to obtain a variable compression sub-process G 1 Outlet pressure P of (2) (1)j ;
S213, assume a polytropic compression sub-process G 1 Is the outlet specific entropy s (1)j Specific entropy s with inlet (1)i Equal, i.e. let s (1)j =s (1)i The method comprises the steps of carrying out a first treatment on the surface of the At outlet pressure P (1)j And the outlet specific entropy s (1)j Under the condition of calling NIST database, calculating variable compression subprocess G 1 Is set to the outlet temperature T s(1)j The method comprises the steps of carrying out a first treatment on the surface of the And let the variable compression sub-process G 1 Outlet temperature T of (2) (1)j =T s(1)j ;
S214, at outlet temperature T (1)j And outlet pressure P (1)j Under the condition of calling NIST database, calculating variable compression subprocess G 1 Is the outlet specific enthalpy h of (2) (1)j And the outlet specific entropy s (1)j ;
S215, according to the inlet specific enthalpy h in S211 (1)i Inlet specific entropy s (1)i And inlet temperature T (1)i And the outlet specific enthalpy h in the S214 (1)j Specific entropy s of outlet (1)j And outlet temperature T (1)j Calculating a variable compression sub-process G by using a Mallin model 1 Calculated value W 'of variable compression sub-work' pol(1) And calculating a variable compression sub-process G by using a variable process compression work definition model 1 Variable compressor work definition value W pol(1) ;
S216, judging the variable compression sub-process G 1 Calculated value W 'of variable compression sub-work' pol(1) And a variable compression sub-work definition value W pol(1) If the absolute value of the difference is smaller than or equal to the preset iteration termination criterion epsilon, if not, the outlet temperature T is calculated (1)j Adding a preset temperature increment alpha, namely, let T (1)j =T (1)j +α, update outlet temperature T (1)j Repeating the steps S214-S215 until the variable compression sub-process G 1 Calculated value W 'of variable compression sub-work' pol(1) And a variable compression sub-work definition value W pol(1) The absolute value of the difference is smaller than or equal to the iteration termination criterion epsilon;
s217, when the variable compression sub-process G 1 Calculated value W 'of variable compression sub-work' pol(1) Variable compressor work settingSense value W pol(1) When the absolute value of the difference is smaller than or equal to the iteration termination criterion epsilon, the variable compression subprocess G 1 Calculated value W 'of variable compression sub-work' pol(1) And a variable compression sub-work definition value W pol(1) As the average value of the variable compression sub-process G 1 Variable compression sub-work W pol(1) 。
Further, variable compression sub-process G 1 Final outlet temperature T (1)j For the final updated outlet temperature T in S216 (1)j Variable compression sub-process G 1 Final outlet specific enthalpy h (1)j And the outlet specific entropy s (1)j For, the outlet temperature T in the step S214 (1)j Updated to the final updated outlet temperature T in S216 (1)j Then, the calculated outlet specific enthalpy h (1)j And the outlet specific entropy s (1)j 。
Further, let the variable compression sub-process G x+1 Inlet pressure P of (2) (x+1)i =P (x)j Make the variable compression sub-process G x+1 Is at an inlet temperature T (x+1)i =T (x)j The method comprises the steps of carrying out a first treatment on the surface of the Wherein, x=1, 2, 3....n.. N. -1, n is pipeline compression the total number of variable compression sub-processes of the machine;
the step S22 is specifically that,
s221, at inlet pressure P (x+1)i And inlet temperature T (x+1)i Under the condition of calling NIST database, calculating variable compression subprocess G x+1 Is the inlet specific enthalpy h of (2) (x+1)i And inlet specific entropy s (x+1)i ;
S222, according to inlet pressure P (x+1)i And a preset pressure ratio Alfa, by the formula P (x+1)j =P (x+1)i * Alfa, calculating to obtain a variable compression sub-process G x+1 Outlet pressure P of (2) (x+1)j ;
S223, assume a polytropic compression sub-process G x+1 Is the outlet specific entropy s (x+1)j Specific entropy s with inlet (x+1)i Equal, i.e. let s (x+1)j =s (x+1)i The method comprises the steps of carrying out a first treatment on the surface of the At outlet pressure P (x+1)j And the outlet specific entropy s (x+1)j Under the condition of calling NIST database, calculating variable compression subprocess G x+1 Is set to the outlet temperature T s(x+1)j The method comprises the steps of carrying out a first treatment on the surface of the And let the variable compression sub-process G x+1 Outlet temperature T of (2) (x+1)j =T s(x+1)j ;
S224, at outlet temperature T (x+1)j And outlet pressure P (x+1)j Under the condition of calling NIST database, calculating variable compression subprocess G x+1 Is the outlet specific enthalpy h of (2) (x+1)j And the outlet specific entropy s (x+1)j ;
S225, according to the inlet specific enthalpy h in S221 (x+1)i Inlet specific entropy s (x+1)i And inlet temperature T (x+1)i And the outlet specific enthalpy h in the S224 (x+1)j Specific entropy s of outlet (x+1)j And outlet temperature T (x+1)j Calculating a variable compression sub-process G by using a Mallin model x+1 Calculated value W 'of variable compression sub-work' pol(x+1) And calculating a variable compression sub-process G by using a variable process compression work definition model x+1 Variable compressor work definition value W pol(x+1) ;
S226, judging the variable compression sub-process G x+1 Calculated value W 'of variable compression sub-work' pol(x+1) And a variable compression sub-work definition value W pol(x+1) If the absolute value of the difference is smaller than or equal to the preset iteration termination criterion epsilon, if not, the outlet temperature T is calculated (x+1)j Adding a preset temperature increment alpha, namely, let T (x+1)j =T (x+1)j +α, update outlet temperature T (x+1)j Repeating the steps S224-S225 until the variable compression sub-process G x+1 Calculated value W 'of variable compression sub-work' pol(x+1) And a variable compression sub-work definition value W pol(x+1) The absolute value of the difference is smaller than or equal to the iteration termination criterion epsilon;
s227, when the variable compression sub-process G x+1 Calculated value W 'of variable compression sub-work' pol(x+1) And a variable compression sub-work definition value W pol(x+1) When the absolute value of the difference is smaller than or equal to the iteration termination criterion epsilon, the variable compression subprocess G x+1 Calculated value W 'of variable compression sub-work' pol(x+1) And a variable compression sub-work definition value W pol(x+1) Average value of (2)As a variable compression sub-process G x+1 Variable compression sub-work W pol(x+1) 。
Further, variable compression sub-process G x+1 Final outlet temperature T (x+1)j For the final updated outlet temperature T in said S226 (x+1)j Variable compression sub-process G x+1 Final outlet specific enthalpy h (x+1)j And the outlet specific entropy s (x+1)j For, the outlet temperature T in the S224 (x+1)j Updated to the final updated outlet temperature T in S226 (x+1)j Then, the calculated outlet specific enthalpy h (x+1)j And the outlet specific entropy s (x+1)j 。
Based on the energy consumption measuring method of the pipeline compressor, the invention further provides an energy consumption measuring system of the pipeline compressor.
A pipeline compressor energy consumption measurement system comprises the following modules,
the variable compression process decomposition module is used for dividing the whole variable compression process of the pipeline compressor into a plurality of variable compression subprocesses according to the inlet pressure and the outlet pressure of the pipeline compressor;
the variable compression sub-work calculation module is used for sequentially solving the variable compression sub-work of each variable compression sub-process based on a variable power solving model;
and the variable compression sub-work superposition summation module is used for superposing variable compression sub-work of all variable compression sub-processes to obtain variable compression work of the whole variable compression process of the pipeline compressor.
Based on the pipeline compressor energy consumption measurement method, the invention further provides a computer storage medium.
A computer storage medium comprising a memory, and a computer program stored on the memory, which when executed by a processor implements a pipeline compressor energy measurement method as described above.
The beneficial effects of the invention are as follows: the invention calculates the compression work of the multi-transformation process by adopting the direct integral model, and the direct integral model divides the whole multi-transformation compression process into a plurality of subprocesses and solves each subprocess respectively, so that the invention can better simulate the change of natural gas in the compression process, thereby the application range of the direct integral model is wider and the invention can be suitable for the compression process with phase transformation. In addition, the invention adopts a more accurate variable power solving model, namely a Malen model, to solve the variable compression sub-power of the variable compression sub-process, and solves the variable compression sub-power of the whole variable compression process from the first variable compression sub-process to the last variable compression sub-process, and then adds the variable compression sub-power of all the variable compression sub-processes to obtain the variable compression power of the whole variable compression process, so that the variable compression power of the whole variable compression process can be calculated more accurately; meanwhile, the calculation method avoids calculation errors caused by solving the polytropic process index m.
Detailed Description
The principles and features of the present invention are described below with reference to the drawings, the examples are illustrated for the purpose of illustrating the invention and are not to be construed as limiting the scope of the invention.
As shown in fig. 1, a pipeline compressor energy consumption measurement method includes the steps of,
s1, dividing the whole variable compression process of a pipeline compressor into a plurality of variable compression sub-processes according to inlet and outlet pressures of the pipeline compressor;
s2, sequentially solving the variable compression sub-work of each variable compression sub-process based on a variable power solving model;
and S3, superposing the variable compression sub-work of all the variable compression sub-processes to obtain the variable compression work of the whole variable compression process of the pipeline compressor.
In this particular embodiment:
in S1, the entire variable compression process of the pipeline compressor may be equally divided into a plurality of variable compression sub-processes.
The specific step S2 is that,
s21, based on the multiple transformation solving model and the set condition M, according to the known inlet condition K of the pipeline compressor 0 Calculating a variable compression sub-process G 1 Outlet condition K of (2) 1 And variable compression sub-work W pol(1) ;
S22, variable compression sub-process G x Outlet condition K of (2) x As a variable compression sub-process G x+1 Inlet condition F of (2) x+1 The method comprises the steps of carrying out a first treatment on the surface of the In the variable compression sub-process G x+1 Based on the variable power solving model and the set condition M, according to the variable compression subprocess G x+1 Inlet condition F of (2) x+1 Calculating a variable compression sub-process G x+1 Outlet condition K of (2) x+1 And variable compression sub-work W pol(x+1) The method comprises the steps of carrying out a first treatment on the surface of the Where x=1, 2,3. N-1, n is a pipeline total number of variable compression sub-processes of the compressor.
Still further, the known inlet condition K 0 Including the inlet pressure P1 and inlet temperature T1 of the pipeline compressor; wherein the known inlet condition K 0 Is a variable compression sub-process G 1 Inlet conditions of (2);
the setting condition M comprises a preset pressure ratio Alfa, a preset iteration termination criterion epsilon and a preset temperature increment alpha;
variable compression sub-process G x Outlet condition K of (2) x Including the outlet pressure P (x)j Outlet temperature T (x)j Specific enthalpy of outlet h (x)j And the outlet specific entropy s (x)j ;
Variable compression sub-process G x+1 Inlet condition F of (2) x+1 Including the inlet pressure P (x+1)i Inlet temperature T (x+1)i Specific enthalpy of inlet h (x+1)i And inlet specific entropy s (x+1)i The method comprises the steps of carrying out a first treatment on the surface of the And P is (x+1)i =P (x)j ,T (x+1)i =T (x)j ,h (x+1)i =h (x)j ,s (x+1)i =s (x)j ;
The calculation termination condition of the S2 is a variable compression subprocess G n Outlet of (2)Condition K n Outlet pressure P in (a) (n)j Greater than or equal to the known outlet condition of the pipeline compressor, wherein the known outlet condition of the pipeline compressor is the outlet pressure P2 of the pipeline compressor.
Further, the variable power solving model comprises a Mallin model and a variable process compression power definition model;
the mullen model is specifically described as,
wherein W'
pol(x) Is a variable compression sub-process G
x Calculated value of variable compression sub-work, h
(x)i 、s
(x)i And T
(x)i Respectively the multiple compression subprocess G
x Inlet specific enthalpy, inlet specific entropy and inlet temperature, and h
(x)i =h
(x-1)j ,s
(x)i =s
(x-1)j ,T
(x)i =T
(x-1)j ;
The variable process compression work definition model is specifically W pol(x) =η pol (h (x)j -h (x)i ) Wherein W' pol(x) Is a variable compression sub-process G x Is defined by the variable compression sub-work values, eta pol Is the variable compression efficiency of the pipeline compressor in the whole variable compression process.
Natural gas physical parameters can be obtained by calling the National Institute of Standards and Technology (NIST) database. For example, v, s, and h can be solved under the condition of P, T; v, h, T can be solved under the known P, s condition; v, s, T can be solved under the known P, h condition; wherein P is pressure, T is temperature, v is specific volume, s is specific entropy, and h is specific enthalpy.
Specifically, let the variable compression sub-process G 1 Inlet pressure P of (2) (1)i Let variable compression sub-process G =p1 1 Is at an inlet temperature T (1)i =T1;
The step S21 is specifically that,
s211 at inlet pressure P (1)i And inlet temperature T (1)i Under the condition of calling NIST database and calculating variable compressionsub-Process G 1 Is of inlet specific volume v (1)i Specific enthalpy of inlet h (1)i And inlet specific entropy s (1)i ;
S212, according to inlet pressure P (1)i And a preset pressure ratio Alfa, by the formula P (1)j =P (1)i * Alfa, calculating to obtain a variable compression sub-process G 1 Outlet pressure P of (2) (1)j ;
S213, assume a polytropic compression sub-process G 1 Is the outlet specific entropy s (1)j Specific entropy s with inlet (1)i Equal, i.e. let s (1)j =s (1)i The method comprises the steps of carrying out a first treatment on the surface of the At outlet pressure P (1)j And the outlet specific entropy s (1)j Under the condition of calling NIST database, calculating variable compression subprocess G 1 Is set to the outlet temperature T s(1)j The method comprises the steps of carrying out a first treatment on the surface of the And let the variable compression sub-process G 1 Outlet temperature T of (2) (1)j =T s(1)j ;
S214, at outlet temperature T (1)j And outlet pressure P (1)j Under the condition of calling NIST database, calculating variable compression subprocess G 1 Specific volume v of outlet of (2) (1)j Specific enthalpy of outlet h (1)j And the outlet specific entropy s (1)j ;
S215, according to the inlet specific enthalpy h in S211 (1)i Inlet specific entropy s (1)i And inlet temperature T (1)i And the outlet specific enthalpy h in the S214 (1)j Specific entropy s of outlet (1)j And outlet temperature T (1)j Calculating a variable compression sub-process G by using a Mallin model 1 Calculated value W 'of variable compression sub-work' pol(1) And calculating a variable compression sub-process G by using a variable process compression work definition model 1 Variable compressor work definition value W pol(1) ;
S216, judging the variable compression sub-process G 1 Calculated value W 'of variable compression sub-work' pol(1) And a variable compression sub-work definition value W pol(1) If the absolute value of the difference is smaller than or equal to the preset iteration termination criterion epsilon, if not, the outlet temperature T is calculated (1)j Adding a preset temperature increment alpha, namely, let T (1)j =T (1)j +α, update outlet temperature T (1)j Repeating the steps S214-S215 until the variable compression sub-process G 1 Calculated value W 'of variable compression sub-work' pol(1) And a variable compression sub-work definition value W pol(1) The absolute value of the difference is smaller than or equal to the iteration termination criterion epsilon;
s217, when the variable compression sub-process G 1 Calculated value W 'of variable compression sub-work' pol(1) And a variable compression sub-work definition value W pol(1) When the absolute value of the difference is smaller than or equal to the iteration termination criterion epsilon, the variable compression subprocess G 1 Calculated value W 'of variable compression sub-work' pol(1) And a variable compression sub-work definition value W pol(1) As the average value of the variable compression sub-process G 1 Variable compression sub-work W pol(1) 。
Further, variable compression sub-process G 1 Final outlet temperature T (1)j For the final updated outlet temperature T in S216 (1)j Variable compression sub-process G 1 Final outlet specific enthalpy h (1)j And the outlet specific entropy s (1)j For, the outlet temperature T in the step S214 (1)j Updated to the final updated outlet temperature T in S216 (1)j Then, the calculated outlet specific enthalpy h (1)j And the outlet specific entropy s (1)j 。
Specifically, let the variable compression sub-process G x+1 Inlet pressure P of (2) (x+1)i =P (x)j Make the variable compression sub-process G x+1 Is at an inlet temperature T (x+1)i =T (x)j The method comprises the steps of carrying out a first treatment on the surface of the Wherein, x=1, 2, 3....n.. N. -1, n is pipeline compression the total number of variable compression sub-processes of the machine;
the step S22 is specifically that,
s221, at inlet pressure P (x+1)i And inlet temperature T (x+1)i Under the condition of calling NIST database, calculating variable compression subprocess G x+1 Is of inlet specific volume v (x+1)i Specific enthalpy of inlet h (x+1)i And inlet specific entropy s (x+1)i ;
S222, according to inlet pressure P (x+1)i And a preset pressure ratio Alfa, by the formula P (x+1)j =P (x+1)i * Alfa, calculating to obtain a variable compression sub-process G x+1 Outlet pressure P of (2) (x+1)j ;
S223, assume a polytropic compression sub-process G x+1 Is the outlet specific entropy s (x+1)j Specific entropy s with inlet (x+1)i Equal, i.e. let s (x+1)j =s (x+1)i The method comprises the steps of carrying out a first treatment on the surface of the At outlet pressure P (x+1)j And the outlet specific entropy s (x+1)j Under the condition of calling NIST database, calculating variable compression subprocess G x+1 Is set to the outlet temperature T s(x+1)j The method comprises the steps of carrying out a first treatment on the surface of the And let the variable compression sub-process G x+1 Outlet temperature T of (2) (x+1)j =T s(x+1)j ;
S224, at outlet temperature T (x+1)j And outlet pressure P (x+1)j Under the condition of calling NIST database, calculating variable compression subprocess G x+1 Specific volume v of outlet of (2) (x+1)j Specific enthalpy of outlet h (x+1)j And the outlet specific entropy s (x+1)j ;
S225, according to the inlet specific enthalpy h in S221 (x+1)i Inlet specific entropy s (x+1)i And inlet temperature T (x+1)i And the outlet specific enthalpy h in the S224 (x+1)j Specific entropy s of outlet (x+1)j And outlet temperature T (x+1)j Calculating a variable compression sub-process G by using a Mallin model x+1 Calculated value W 'of variable compression sub-work' pol(x+1) And calculating a variable compression sub-process G by using a variable process compression work definition model x+1 Variable compressor work definition value W pol(x+1) ;
S226, judging the variable compression sub-process G x+1 Calculated value W 'of variable compression sub-work' pol(x+1) And a variable compression sub-work definition value W pol(x+1) If the absolute value of the difference is smaller than or equal to the preset iteration termination criterion epsilon, if not, the outlet temperature T is calculated (x+1)j Adding a preset temperature increment alpha, namely, let T (x+1)j =T (x+1)j +α, update outlet temperature T (x+1)j Repeating the steps S224-S225 until the variable compression sub-process G x+1 Calculated value W 'of variable compression sub-work' pol(x+1) And a variable compression sub-work definition value W pol(x+1) The absolute value of the difference is smaller than or equal to the iteration termination criterion epsilon;
s227, when the variable compression sub-process G x+1 Calculated value W 'of variable compression sub-work' pol(x+1) And a variable compression sub-work definition value W pol(x+1) When the absolute value of the difference is smaller than or equal to the iteration termination criterion epsilon, the variable compression subprocess G x+1 Calculated value W 'of variable compression sub-work' pol(x+1) And a variable compression sub-work definition value W pol(x+1) As the average value of the variable compression sub-process G x+1 Variable compression sub-work W pol(x+1) 。
Further, variable compression sub-process G x+1 Final outlet temperature T (x+1)j For the final updated outlet temperature T in said S226 (x+1)j Variable compression sub-process G x+1 Final outlet specific enthalpy h (x+1)j And the outlet specific entropy s (x+1)j For, the outlet temperature T in the S224 (x+1)j Updated to the final updated outlet temperature T in S226 (x+1)j Then, the calculated outlet specific enthalpy h (x+1)j And the outlet specific entropy s (x+1)j 。
In summary, the variable compression work of the whole variable compression process is:
in this specific embodiment, the S21 and the S22 are integrated together to obtain the integrated solution shown in fig. 2. As shown in fig. 2:
the conditions are known: inlet temperature T of pipeline compressor 1 Inlet pressure p 1 Outlet pressure p 2 Variable compression efficiency eta in the whole variable compression process pol 。
Setting conditions: the iteration termination criterion epsilon, the small temperature increment alpha and the pressure ratio Alfa of the variable compression sub-process.
The comprehensive scheme is as follows:
the whole variable compression process is first divided into a number of variable compression sub-processes. First multipleInlet temperature T of variable compression sub-process i =T 1 Inlet pressure p i =p 1 The following variable compression sub-work W of any one variable compression sub-process poli The calculation steps of (a) are explained.
(1) Invoking NIST database to calculate variable compression subprocess inlet temperature T i Inlet pressure p i Specific volume v of natural gas under the condition i Specific enthalpy h i Sum of specific entropy s i ;
(2) Variable compression sub-process outlet pressure p j =p i * Alfa, judgment of p j And p is as follows 2 If p j >p 2 Let p j =p 2 ;
(3) Variable compression sub-process outlet pressure p j Assuming that the natural gas specific entropy of the outlet of the variable compression sub-process is equal to that of the inlet, and s is the same as that of the inlet i Invoking NIST database to calculate specific enthalpy h of natural gas under the condition sj Temperature T sj And specific volume v sj . Let the outlet temperature T of the variable compression sub-process j =T sj ;
(4) The outlet temperature of the variable compression sub-process is T j The outlet pressure is p j Calling NIST database to calculate specific volume v of natural gas under the condition j Specific enthalpy h j Sum of specific entropy s j ;
(5) Step (1) obtaining an inlet parameter, step (4) obtaining an outlet parameter, and substituting the calculation results into the step (4) to calculate so as to obtain variable power W poli1 Substituting into (5) to calculate to obtain variable power W poli2 ;
(6) Definition W poli1 And W is poli2 The absolute value of the difference is abs (W poli1 -W poli2 ). If abs (W) poli1 -W poli2 )>ε, give T j Adding alpha and updating T j Repeating steps (4) and (5) until abs (W poli1 -W poli2 )≤ε;
(7) Let W poli =(W poli1 +W poli2 ) And/2, obtaining the variable work W of the variable compression sub-process poli Is changeable in value of (a)Variable power W in compression sub-process poli Is completed.
The outlet condition of the last variable compression sub-process is the inlet condition of the next variable compression sub-process, i.e. after calculating the variable power of one variable compression sub-process, let p i =p j ,T i =T j The computation of the multiple work of the next multiple compression sub-process can be entered, and the end condition of the computation of the whole multiple compression process is the outlet pressure p of a certain multiple compression sub-process j ≥p 2 . And adding the variable work of all the variable compression sub-processes to obtain the variable compression work of the whole variable compression process.
In this embodiment, the inlet and outlet parameters of the pipeline compressor for a natural gas are as follows: the inlet temperature was 307.55K, the inlet pressure was 6927kPa, the outlet temperature was 338.95K, the outlet pressure was 6927kPa, and the variable compression efficiency η pol 87.4%. The variable power of the compression process of the natural gas pipeline compressor is calculated by the method of the invention and the method of the prior art (the calculation method of the formula (1) in the background art), and the result shown in the following table 1 is obtained. As can be seen from Table 1, the calculation accuracy of the method of the present invention is greatly improved compared with that of the conventional method.
Table 1: different methods for calculating the results of variable compression work of pipeline compressors
Based on the energy consumption measuring method of the pipeline compressor, the invention further provides an energy consumption measuring system of the pipeline compressor.
As shown in fig. 3, a pipeline compressor energy consumption measurement system, comprising the following modules,
the variable compression process decomposition module is used for dividing the whole variable compression process of the pipeline compressor into a plurality of variable compression subprocesses according to the inlet pressure and the outlet pressure of the pipeline compressor;
the variable compression sub-work calculation module is used for sequentially solving the variable compression sub-work of each variable compression sub-process based on a variable power solving model;
and the variable compression sub-work superposition summation module is used for superposing variable compression sub-work of all variable compression sub-processes to obtain variable compression work of the whole variable compression process of the pipeline compressor.
Based on the pipeline compressor energy consumption measurement method, the invention further provides a computer storage medium.
A computer storage medium comprising a memory, and a computer program stored on the memory, which when executed by a processor implements a pipeline compressor energy measurement method as described above.
The invention calculates the compression work of the multi-transformation process by adopting the direct integral model, and the direct integral model divides the whole multi-transformation compression process into a plurality of subprocesses and solves each subprocess respectively, so that the invention can better simulate the change of natural gas in the compression process, thereby the application range of the direct integral model is wider and the invention can be suitable for the compression process with phase transformation. In addition, the invention adopts a more accurate variable power solving model, namely a Malen model, to solve the variable compression sub-power of the variable compression sub-process, and solves the variable compression sub-power of the whole variable compression process from the first variable compression sub-process to the last variable compression sub-process, and then adds the variable compression sub-power of all the variable compression sub-processes to obtain the variable compression power of the whole variable compression process, so that the variable compression power of the whole variable compression process can be calculated more accurately; meanwhile, the calculation method avoids calculation errors caused by solving the polytropic process index m.
The foregoing description of the preferred embodiments of the invention is not intended to limit the invention to the precise form disclosed, and any such modifications, equivalents, and alternatives falling within the spirit and scope of the invention are intended to be included within the scope of the invention.