CN-117906066-B - Calculation method and device for hydraulic parameters of natural gas station

CN117906066BCN 117906066 BCN117906066 BCN 117906066BCN-117906066-B

Abstract

The application provides a method and a device for calculating hydraulic parameters of a natural gas station. The method comprises the steps of obtaining structural parameters of a target natural gas station and structural parameters of a target pipeline, constructing a target natural gas station model according to the structural parameters of the target natural gas station, constructing a target pipeline model according to the structural parameters of the target pipeline, determining first target pressure corresponding to an outlet node of a compressor, determining pipeline boundary conditions corresponding to pipeline boundary nodes, calculating gas density, boundary flow and boundary pressure corresponding to the boundary nodes of the station according to the pipeline boundary conditions, determining first target flow corresponding to the compressor and an adjusting valve and second target pressure corresponding to the station according to the gas density, the boundary flow and the boundary pressure, calculating power of the compressor according to the first target flow corresponding to the compressor, and correcting the first target flow corresponding to the compressor by utilizing the first target pressure corresponding to the outlet node of the compressor to obtain second target flow corresponding to the compressor.

Inventors

LI WEIJIA
GONG JING
FAN DI
LIU PENG
HUANG WEIHE

Assignees

中国石油大学(北京)

Dates

Publication Date: 20260508
Application Date: 20240110

Claims (10)

1. The method for calculating the hydraulic parameters of the natural gas station is characterized by comprising the following steps of: the method comprises the steps of obtaining structural parameters of a target natural gas station and structural parameters of a target pipeline, constructing a target natural gas station model according to the structural parameters of the target natural gas station, constructing a target pipeline model according to the structural parameters of the target pipeline, wherein the target natural gas station model comprises a plurality of station nodes, a plurality of compressors and a plurality of regulating valves, the station nodes comprise station internal nodes and station boundary nodes, the station internal nodes comprise compressor inlet nodes, compressor outlet nodes, regulating valve inlet nodes and regulating valve outlet nodes, the target pipeline model comprises a plurality of pipeline nodes, the pipeline nodes comprise pipeline boundary nodes and pipeline internal nodes, and the target natural gas station model is connected with the target pipeline model through the station boundary nodes; Determining a first target pressure corresponding to the compressor outlet node, and determining a pipeline boundary condition corresponding to the pipeline boundary node; Calculating the gas density, the boundary flow and the boundary pressure corresponding to the station boundary nodes according to the pipeline boundary conditions; determining a first target flow corresponding to the compressor, a first target flow corresponding to the regulating valve and a second target pressure corresponding to the station node according to the gas density, the boundary flow and the boundary pressure; Calculating the power of the compressor according to the first target flow corresponding to the compressor; and correcting the first target flow corresponding to the compressor by using the first target pressure corresponding to the outlet node of the compressor to obtain the second target flow corresponding to the compressor.
2. The method of claim 1, wherein calculating the gas density, boundary flow, boundary pressure corresponding to the yard boundary node according to the pipeline boundary condition comprises: constructing a natural gas pipeline flow equation corresponding to the target pipeline model; and carrying out discrete processing on the natural gas pipeline flow equation based on the pipeline boundary condition, and solving to obtain the gas density, the boundary flow and the boundary pressure corresponding to the station boundary node.
3. The method of claim 1, wherein determining a first target flow rate for the compressor, a first target flow rate for the regulator valve, and a second target pressure for the yard node based on the gas density, the boundary flow rate, and the boundary pressure comprises: According to the gas density, the boundary flow and the boundary pressure, solving a compressor flow pressure equation of the current round and a regulating valve flow pressure equation of the current round to obtain a first flow corresponding to the current round compressor and a first flow corresponding to the current round regulating valve; constructing a flow change equation of the current turn according to the first flow corresponding to the current turn compressor, the first flow corresponding to the current turn regulating valve and the second pressure corresponding to the current turn station node; Solving a flow change equation of the current turn to obtain the pressure increment of the current turn station yard node; Detecting whether the pressure increment of the current round station yard node is smaller than a preset pressure increment precision threshold value or not; Under the condition that the pressure increment of the current round station node is smaller than a preset pressure increment precision threshold value, taking the first flow corresponding to the current round compressor as the first target flow corresponding to the compressor, taking the first flow corresponding to the current round regulating valve as the first target flow corresponding to the regulating valve, and taking the second pressure corresponding to the current round station node as the second target pressure corresponding to the station node.
4. A method according to claim 3, wherein after detecting whether the pressure increase of the current turn yard node is less than a preset pressure increase accuracy threshold, the method further comprises: under the condition that the pressure increment of the current turn station node is larger than or equal to a preset pressure increment precision threshold value, calculating a second pressure corresponding to the next turn station node according to the pressure increment of the current turn station node; Calculating a first flow corresponding to a next-round compressor and a first flow corresponding to a next-round regulating valve according to a second pressure corresponding to a next-round station node; Reconstructing and solving a flow change equation according to the second pressure corresponding to the next-round station node, the first flow corresponding to the next-round compressor and the first flow corresponding to the next-round regulating valve to obtain the pressure increment of the next-round station node; detecting whether the pressure increment of the next-round station yard node is smaller than a preset pressure increment precision threshold value or not; and under the condition that the pressure increment of the next round is smaller than a preset pressure increment precision threshold value, taking the first flow corresponding to the next round of compressor as the first target flow corresponding to the compressor, taking the first flow corresponding to the next round of regulating valve as the first target flow corresponding to the regulating valve, and taking the second pressure corresponding to the next round of station node as the second target pressure corresponding to the station node.
5. The method of claim 3, wherein constructing a current turn flow rate change equation based on the first flow rate corresponding to the current turn compressor, the first flow rate corresponding to the current turn regulator valve, and the second pressure corresponding to the current turn yard node, comprises: Constructing a node flow balance equation of the current turn according to the first flow corresponding to the current turn compressor and the first flow corresponding to the current turn regulating valve; calculating a jacobian matrix corresponding to a node flow balance equation of the current turn according to the second pressure corresponding to the current turn station node; And obtaining a current round flow change equation according to the current round node flow balance equation and the jacobian matrix corresponding to the current round node flow balance equation.
6. A method according to claim 3, characterized in that the method further comprises: Under the condition that the pressure increment of the current round of station nodes is smaller than a preset pressure increment precision threshold value, extracting updated pressure difference corresponding to the station boundary nodes from the second target pressure; And calculating the boundary flow and boundary pressure corresponding to the boundary node of the next time step according to the updated pressure difference corresponding to the boundary node of the station and the boundary condition of the pipeline of the next time step.
7. The method of claim 1, wherein calculating compressor power based on the first target flow rate for the compressor comprises: The compressor power is calculated according to the following formula: Wherein, the Which is indicative of the power of the compressor, A first target flow corresponding to the compressor is represented, g represents gravitational acceleration, The multi-energy-changing head is represented, Indicating compressor efficiency.
8. The method of claim 1, wherein modifying the first target flow rate for the compressor with the first target pressure for the compressor outlet node to obtain the second target flow rate for the compressor comprises: Determining sampling point flow corresponding to a sampling point according to first target flow corresponding to the compressor; calculating a pressure calculation value corresponding to the sampling point according to the sampling point flow; Calculating a gain matrix corresponding to the sampling point according to the sampling point flow and the pressure calculation value; And obtaining a second target flow corresponding to the compressor according to the pressure estimated value, the gain matrix and the first target pressure corresponding to the compressor outlet node.
9. A computing device for hydraulic parameters of a natural gas farm, comprising: The system comprises a building module, a target natural gas station model, a pipeline model and a target pipeline model, wherein the building module is used for acquiring structural parameters of a target natural gas station and structural parameters of a target pipeline, building the target natural gas station model according to the structural parameters of the target natural gas station, building the target pipeline model according to the structural parameters of the target pipeline, the target natural gas station model comprises a plurality of station nodes, a plurality of compressors and a plurality of regulating valves, the station nodes comprise station internal nodes and station boundary nodes, the station internal nodes comprise compressor inlet nodes, compressor outlet nodes, regulating valve inlet nodes and regulating valve outlet nodes, the target pipeline model comprises a plurality of pipeline nodes, the pipeline nodes comprise pipeline boundary nodes and pipeline internal nodes, and the target natural gas station model is connected with the target pipeline model through the station boundary nodes; the determining module is used for determining a first target pressure corresponding to the outlet node of the compressor and determining a pipeline boundary condition corresponding to the pipeline boundary node; The first calculation module is used for calculating the gas density, the boundary flow and the boundary pressure corresponding to the station boundary node according to the pipeline boundary condition; the second calculation module is used for determining a first target flow corresponding to the compressor, a first target flow corresponding to the regulating valve and a second target pressure corresponding to the station node according to the gas density, the boundary flow and the boundary pressure; the third calculation module is used for calculating the power of the compressor according to the first target flow corresponding to the compressor; And the correction module is used for correcting the first target flow corresponding to the compressor by using the first target pressure corresponding to the outlet node of the compressor to obtain the second target flow corresponding to the compressor.
10. A computer readable storage medium, having stored thereon computer instructions which, when executed by a processor, implement the steps of the method of any of claims 1 to 8.

Description

Calculation method and device for hydraulic parameters of natural gas station Technical Field The specification belongs to the technical field of natural gas pipeline system monitoring and control, and particularly relates to a calculation method and a calculation device for hydraulic parameters of a natural gas station. Background The natural gas station is a junction of a natural gas pipe network, and various devices for pressurizing, depressurizing, heating and cooling, such as a compressor, a pressure regulating valve, a filtering device and the like, are included in the natural gas station, and are connected with the pipeline in a complex series connection and a complex parallel connection. The natural gas station has the advantages of large number of internal devices, various connection modes and complex process flow. But the station equipment is seriously aged, and the problem of energy consumption increase is generated. On one hand, the determination of the equipment of the backward capacity station can be judged only by operation and maintenance personnel through manual experience and gradually updated or eliminated according to an annual maintenance and transformation plan, and the working accuracy is lacking, and on the other hand, the operation and maintenance personnel cannot scientifically and effectively determine the equipment output according to the energy consumption characteristics of the equipment, so that the station equipment operates in a non-optimal energy efficiency ratio state, and the energy utilization rate is insufficient. In order to solve the problems, the monitoring of energy consumption of the equipment in the operation process needs to be enhanced, a source of high energy consumption is found, energy-saving management and equipment replacement are implemented, and then the emission of carbon dioxide is reduced. In order to know the flowing condition of gas in the natural gas station, the operation efficiency of equipment is estimated, the energy consumption of the equipment is monitored, a simulation model of the natural gas station is required to be established, the simulation model of the natural gas station is subjected to simulation operation, flow data and pressure data in the natural gas station are calculated, and the flow data and the pressure data are used for calculating the energy consumption of the equipment. In the prior art, when the pipeline and the station are jointly calculated, the natural gas station cannot be subjected to refined simulation modeling due to the limitation of calculation speed and calculation accuracy, and the structure and equipment in the natural gas station are generally simplified, so that the simulation operation is difficult to accurately perform, the flow data and the pressure data in the natural gas station are difficult to calculate, and the accurate energy consumption data are difficult to calculate. In view of the above technical problems, no effective solution has been proposed at present. Disclosure of Invention The application provides a calculation method and a calculation device for hydraulic parameters of a natural gas station, which can determine the hydraulic parameters (including flow data of a compressor and a regulating valve, pressure data of station nodes and compressor power) of the natural gas station with any topological structure. The embodiment of the application aims to provide a calculation method of hydraulic parameters of a natural gas station, which comprises the following steps: the method comprises the steps of obtaining structural parameters of a target natural gas station and structural parameters of a target pipeline, constructing a target natural gas station model according to the structural parameters of the target natural gas station, constructing a target pipeline model according to the structural parameters of the target pipeline, wherein the target natural gas station model comprises a plurality of station nodes, a plurality of compressors and a plurality of regulating valves, the station nodes comprise station internal nodes and station boundary nodes, the station internal nodes comprise compressor inlet nodes, compressor outlet nodes, regulating valve inlet nodes and regulating valve outlet nodes, the target pipeline model comprises a plurality of pipeline nodes, the pipeline nodes comprise pipeline boundary nodes and pipeline internal nodes, and the target natural gas station model is connected with the target pipeline model through the station boundary nodes; Determining a first target pressure corresponding to the compressor outlet node, and determining a pipeline boundary condition corresponding to the pipeline boundary node; Calculating the gas density, the boundary flow and the boundary pressure corresponding to the station boundary nodes according to the pipeline boundary conditions; determining a first target flow corresponding to the compressor, a first target flow corresponding to th