CN115630479A - Gas tracking calculation method and device in multi-gas-source pipe network, server and storage medium - Google Patents

Gas tracking calculation method and device in multi-gas-source pipe network, server and storage medium Download PDF

Info

Publication number
CN115630479A
CN115630479A CN202211152701.8A CN202211152701A CN115630479A CN 115630479 A CN115630479 A CN 115630479A CN 202211152701 A CN202211152701 A CN 202211152701A CN 115630479 A CN115630479 A CN 115630479A
Authority
CN
China
Prior art keywords
gas
batch
target
time
pipe network
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202211152701.8A
Other languages
Chinese (zh)
Inventor
邓传忠
李伟
高达
夏如君
李力
崔嵘
何伟
曾桃
李跃喜
曾庆军
韩玉林
王继鹏
曾小明
曾楠诺
涂少勇
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
CNOOC China Ltd Hainan Branch
Original Assignee
CNOOC China Ltd Hainan Branch
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by CNOOC China Ltd Hainan Branch filed Critical CNOOC China Ltd Hainan Branch
Priority to CN202211152701.8A priority Critical patent/CN115630479A/en
Publication of CN115630479A publication Critical patent/CN115630479A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10Complex mathematical operations

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Data Mining & Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Geometry (AREA)
  • Computational Mathematics (AREA)
  • Evolutionary Computation (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Computer Hardware Design (AREA)
  • Pure & Applied Mathematics (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention is suitable for the technical field of fluid pipeline network transportation, and provides a method, a device, a server and a storage medium for gas tracking calculation in a multi-gas-source pipeline network, wherein the method comprises the following steps: s1: acquiring real-time data of all boundary nodes of the gas pipe network; the boundary nodes of the gas pipe network comprise a gas source point and a user point; s2: according to the real-time data in the step S1, dividing the gas output from the gas source point in batches to determine a plurality of target gas batches; s3: according to the target gas batches determined in the step S2, tracking calculation is carried out on each target gas batch at a mixing point, and a calculation result of each target gas batch is obtained; s4: acquiring real-time gas information in the gas pipe network at the time point to be inquired according to the calculation result of the target gas batch in the step S3; based on the pipe network online simulation technology, the output gas at the gas source point is divided into batches, and the gas groups after batch division can be more fit with the real-time state of the pipe network gas.

Description

Gas tracking calculation method and device in multi-gas-source pipe network, server and storage medium
Technical Field
The invention belongs to the technical field of fluid pipeline network conveying, and particularly relates to a method and a device for tracking and calculating gas in a multi-gas-source pipeline network, a server and a storage medium.
Background
Different gas fields (gas sources) on the sea are influenced by the condition factors of oil and gas reservoirs, the component heat values of produced natural gas are different, and the positions of each gas field and land and shore terminal users are dispersed, so that the seabed multi-gas-source combined gathering and transportation pipe network is provided with a plurality of mixing points/branch transportation points, and monitoring instruments are not arranged at the mixing points/branch transportation points, so that the component heat values of mixed gas downstream of each seabed gas mixing point and the distribution and the migration of the mixed gas in a pipeline are difficult to determine. The end user has definite requirements on the quality of the mixed gas, the production side needs to know whether the quality of the mixed gas in the submarine pipeline meets the gas supply requirement or not in time, if the quality of the mixed gas exceeds or is close to the quality of the gas which exceeds the requirements of the user, or needs to coordinate the gas source of the next mixing point in time to adjust the external gas output quantity or the quality of the gas, and the quality of the mixed gas is qualified through mixing again. Therefore, a method for controlling a mixed gas real-time tracking technology in a pipe network, controlling the mixed gas real-time pressure, flow rate, gas quality distribution and migration conditions of any node in the pipe network, coordinating the gas supply amount of each mixing point and enabling a gas field to input qualified gas quality to a user side is needed.
At present, tracking of the heat value of a mixed gas component of an offshore multi-gas-source pipe network is usually carried out based on the experience of production personnel; the method comprises the steps of judging mixed components and heat values according to instantaneous flow and components of a gas source, judging gas quality arrival time according to pipe stock and end user gas consumption cumulant, wherein an actual pipe network system is a dynamic change process, the instantaneous flow of the gas source is not equal to the instantaneous flow of a mixing point, and when gas well production is unstable or distribution changes, the components of the gas source fluctuate, so that the component condition of the mixed gas source cannot be accurately judged.
The Chinese patent discloses a fluid pipe network air source flow tracking calculation method and a system, which respectively acquire real-time state data of fluid states of all boundary nodes in the fluid pipe network and physical parameter data for expressing pipe section structures and fluid pipe physical properties; acquiring a first physical property parameter for representing the fluid state of the boundary node according to the real-time state data; acquiring a hydraulic calculation result according to the physical parameter data and the real-time state data, wherein the hydraulic calculation result at least comprises the flow and the pressure drop of each pipe section in the fluid pipe network; the flow tracking is carried out on the flow distribution from each gas source point to a user point mainly by means of the distribution proportion of the pipe sections, but in the gas tracking process, the accuracy in judging the component condition of gas after mixing of multiple gas sources, the tracking of the arrival time of fluid and the tracking of the heat value is not particularly ideal, the tracking prediction of the gas is not accurate enough, and the tracking timeliness is low.
Disclosure of Invention
The invention aims to provide a method, a device, a server and a storage medium for tracking and calculating gas in a multi-gas-source pipe network, and aims to solve the problems that the existing gas tracking method is not particularly ideal in the gas tracking process, the accuracy in judging the component condition of gas after mixing of multiple gas sources, the accuracy in tracking the arrival time of fluid and the accuracy in tracking the heat value are not high, the tracking and prediction of the gas are not accurate enough, and the tracking timeliness is low.
In order to solve the technical problem, the application provides a gas tracking calculation method in a multi-gas-source pipe network, which comprises the following steps:
s1: acquiring real-time data of all boundary nodes of the gas pipe network; the boundary nodes of the gas pipe network comprise gas source points and user points;
s2: according to the real-time data in the step S1, dividing the gas output from the gas source point in batches to determine a plurality of target gas batches;
s3: according to the target gas batches determined in the step S2, tracking calculation is carried out on each target gas batch at a mixing point, and a calculation result of each target gas batch is obtained;
s4: acquiring real-time gas information in the gas pipe network at the time point to be inquired according to the calculation result of the target gas batch in the step S3;
wherein the real-time gas information at least comprises one or more of target gas components, target gas heat value and target gas position change information.
According to the gas tracking calculation method in the multi-gas-source pipe network, based on the pipe network online simulation technology, batch division is performed on output gas at gas source points, so that gas source tracking is not limited by component changes of multiple gas sources and gas consumption changes of users, and gas groups subjected to batch division can be more fit with the real-time state of pipe network gas; the tracking calculation is carried out on the batch obtained by division at the mixing point, so that the tracking calculation result of the gas in the pipe network is more accurate, the component condition of the gas after the mixing of multiple gas sources, the tracking of the arrival time of the fluid and the tracking of the heat value are more accurate, meanwhile, the calculation amount is reduced by carrying out the tracking calculation after the batch division, the tracking calculation time is further shortened, and the timeliness of the gas tracking calculation is improved.
Preferably, the real-time data at least comprises gas data output by a gas source point, hydraulic parameters and parameter data used for representing the physical properties of a pipe network, wherein:
the gas data includes a gas mole percentage at a gas source point;
the hydraulic parameters at least comprise one or a combination of a plurality of pressures, flow rates and temperatures of an air source point and a user point;
the parameter data at least comprises one or a combination of several of the height above sea level of the pipeline, the length of the pipeline and the diameter of the pipeline.
Preferably, in step S2, the gas output from the gas source point is divided into batches, and the step of determining a plurality of target gas batches includes:
s201: the method comprises the steps of performing batch division on gas groups output by each gas source point according to a pipe network static simulation result, and determining a plurality of initial gas batches;
s202: acquiring current mixed data of gas batches actually generated in the pipe network according to the data corresponding to each initial gas batch;
s203: the current mixing data of the gas batch actually generated in step S202 is determined, and the target gas batch is determined according to the determination result.
Preferably, in step S202, the current mixing data of the actually generated gas batch at least includes one or more of a minimum change rate, a fluid update time and a batch volume flow rate, wherein:
minimum rate of change Δ: defined as the sum of the absolute values of the differences between the source point mole percent of each gas component in a gas batch and the source point mole percent of each gas component in a previous gas batch, expressed by the formula:
Figure BDA0003857576530000031
in which n denotes the number of components i in the gas batch, C i Represents the gas source point mole percent of component i in the gas batch, c i Represents the gas source point mole percentage of component i in the previous batch;
maximum fluid update time T: defining the time step length of the latest gas batch actually generated in the pipeline at the gas source and the time step length when the last batch is generated;
minimum batch volume V: the volume of the working condition of the new batch of gas actually generated in the pipeline at the gas source is defined, and the expression formula is as follows:
V=ν·t·A
wherein ν is the gas flow velocity at the judgment position, t is the fluid update time step length at the judgment position, and A is the pipeline cross-sectional area at the judgment position.
Preferably, in step S203, the step of determining the current mixing data of each gas batch acquired in step S202 and determining the target gas batch according to the determination result includes:
s2031: respectively giving a set value for minimum component change rate, minimum batch volume and maximum fluid update time: Δ set, V set, and T set;
s2032: comparing the calculated values of the minimum component change rate Δ, the minimum batch volume V and the maximum fluid update time T of the gas batch obtained in step S202 with the set values determined in step S2031 one by one;
s2033: judging whether the comparison result of the step S2032 meets the conditions that the delta is not less than the delta setting in real time, the V is not less than the V setting in real time, and the T is not less than the T setting in real time;
if so, determining the gas batch as a target gas batch, otherwise determining the previous gas batch as the target gas batch.
Preferably, in step S3, the obtaining of the calculation result of each target gas batch at the mixing point specifically includes:
s301, calculating a first product of the gas source point mole percentage of each component gas in a target gas batch and the volume flow of the target gas batch at a mixing point, and calculating a first sum of the first products corresponding to each component gas in the target gas batch;
s302: calculating a second sum of the volumetric flow rates of each component gas in the target gas batch at the mixing point;
s303: respectively carrying out quotient on the first sum value and the second sum value corresponding to each component gas to obtain the mole percentage of each component gas of the target gas batch at the mixing point;
s304: and determining the calculation result of the target gas batch according to the gas mole percentage of the target gas batch.
Preferably, in step S304, the calculation result at least includes one or more of a heat value of the target batch of gas, position change information, and a target batch arrival time, wherein:
determining the calculation formula of the target batch gas heat value according to the gas mole percentage in the target gas batch as follows:
Figure BDA0003857576530000041
wherein Q is the calorific value of the gas of the target batch, Q j Is the calorific value, C, of component j in the target batch of gas j Is the mole percentage of the gas component j in the target gas batch, and n is the number of the components in the target gas batch;
according to the gas mole percentage in the target gas batch, the calculation formula for determining the position change information of the target batch gas is as follows:
Figure BDA0003857576530000042
in the formula (I), the compound is shown in the specification,
Figure BDA0003857576530000043
indicating the position where the target lot moves at the next time step,
Figure BDA0003857576530000044
indicating the position of the initial time of the target batch, v i The flow velocity of the gas in the unit cell of the pipe section is shown, and delta t represents a time step;
determining the target batch arrival time by the calculation formula according to the gas mole percentage in the target gas batch:
Figure BDA0003857576530000045
where t represents the arrival time of the batch, t 0 Representing the generation time of a certain target batch, and deltax representing the space step of the pipeline; v. of i Indicating the flow rate of the gas in the unit cell of the pipe section.
In order to solve the above technical problem, the present application further provides a gas tracking calculation apparatus, including:
the data acquisition module is used for acquiring the real-time data of the current gas pipe network and transmitting the acquired real-time data to the batch division module;
the batch division module is used for carrying out batch division on the gas groups at the gas source points according to the step S2 and uploading the obtained target gas to the tracking calculation module in batches;
the tracking calculation module is used for performing tracking calculation on the target gas batch according to the step S3 and transmitting the obtained calculation result to the information display module;
the information display module is in signal connection with the data acquisition module, the batch division module and the tracking calculation module, and records and displays the acquired or generated data;
an information input module: the tracking calculation module is arranged on the information display module, is connected with the batch division module and the tracking calculation module, and is used for inputting data and inquiring real-time information of gas in the time node pipe network.
According to the gas tracking calculation device, the data acquisition module is arranged, so that real-time data of offshore multi-gas-source pipes can be conveniently acquired, the batch division module performs batch division on gas groups by using the real-time data as input and control conditions of an online simulation model of a pipeline system, tracking calculation is performed on hydraulic simulation calculation and component tracking calculation of natural gas pipe network flow according to target gas batches, and component data, heat value information and batch position information of different batches of natural gas are output; real-time gas information in the pipe network is displayed through the information display module, and the information input module inputs time nodes needing to be inquired, so that the real-time gas information in the pipe network can be inquired in real time; therefore, the accuracy of fluid arrival time tracking and heat value tracking is effectively improved, the gas tracking quality can be obviously improved, and the timeliness of gas tracking is improved.
In order to solve the above technical problem, the present application further provides a server, including:
a memory for storing a computer program;
and the processor is used for realizing the steps of any one of the gas tracking calculation methods in the multi-gas source pipe network when executing the computer program.
In order to solve the above technical problem, the present application further provides a storage medium, where a computer program is stored, and when the computer program is executed by a processor, the steps of any one of the above methods for calculating gas tracking in a multi-gas-source pipe network are implemented.
Compared with the prior art, the invention has the beneficial effects that: based on a pipe network online simulation technology, the output gas at the gas source point is divided into batches, so that the gas source tracking is not limited by the component change of multiple gas sources and the gas consumption change of users any more, and the gas groups after the batch division can be more fit with the real-time state of the pipe network gas; the tracking calculation is carried out on the batch obtained by division at the mixing point, so that the tracking calculation result of the gas in the pipe network is more accurate, the component condition of the gas mixed by multiple gas sources, the tracking of the arrival time of the fluid and the tracking of the heat value are more accurate, meanwhile, the calculation amount is reduced by carrying out the tracking calculation after the batch division, the tracking calculation time is further shortened, and the timeliness of the gas tracking calculation is improved.
Drawings
FIG. 1 is a flow chart of a gas tracking calculation method in a multi-gas source pipe network according to the present invention;
FIG. 2 is a schematic diagram of a dynamic tracking model of a gas tracking calculation method in a multi-gas-source pipe network according to the present invention;
fig. 3 is a schematic diagram of a gas pipe network of a gas tracking calculation method in a multi-gas-source pipe network according to embodiment 2 of the present invention;
fig. 4 is a schematic diagram of gas source components of a gas tracking calculation method in a multi-gas-source pipe network according to embodiment 2 of the present invention;
fig. 5 is a schematic diagram of batch division of a gas tracking calculation method in a multi-gas-source pipe network according to embodiment 2 of the present invention;
fig. 6 is a statistical chart of batch heat values of a gas tracking calculation method in a multi-gas source pipe network according to embodiment 2 of the present invention;
fig. 7 is a schematic diagram of the batch arrival time of the gas tracking calculation method in the multiple gas source pipe network according to embodiment 2 of the present invention;
fig. 8 is a comparison graph of the heat value tracking result of the gas tracking calculation method in the multi-gas source pipe network according to embodiment 2 of the present invention;
fig. 9 is a schematic structural diagram of a server according to embodiment 4 of the present invention.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is described in further detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are merely illustrative of the invention and are not intended to limit the invention.
Specific implementations of the present invention are described in detail below with reference to specific embodiments.
Example 1
The present embodiment provides a method for calculating gas tracking in a multiple gas source pipe network, as shown in fig. 1, which is a flowchart of the method for calculating gas tracking in a multiple gas source pipe network of the present embodiment, and the method includes the following steps:
s1: acquiring real-time data of all boundary nodes of the gas pipe network; the boundary nodes of the gas pipe network comprise a gas source point and a user point;
s2: according to the real-time data in the step S1, dividing the gas output from the gas source point in batches to determine a plurality of target gas batches;
s3: according to the target gas batches determined in the step S2, tracking calculation is carried out on each target gas batch at a mixing point, and a calculation result of each target gas batch is obtained;
s4: acquiring real-time gas information in the gas pipe network at the time point to be inquired according to the calculation result of the target gas batch in the step S3;
wherein the real-time gas information at least comprises one or more of target gas components, target gas heat value and target gas position change information.
In practical application, based on a pipe network online simulation technology, batch division is performed on output gas at a gas source point, so that gas source tracking is not limited by component changes of multiple gas sources and gas consumption changes of users, and gas groups subjected to batch division can be more suitable for the real-time state of pipe network gas; the tracking calculation is carried out on the batch obtained by division at the mixing point, so that the tracking calculation result of the gas in the pipe network is more accurate, the component condition of the gas after the mixing of multiple gas sources, the tracking of the arrival time of the fluid and the tracking of the heat value are more accurate, meanwhile, the calculation amount is reduced by carrying out the tracking calculation after the batch division, the tracking calculation time is further shortened, and the timeliness of the gas tracking calculation is improved.
Specifically, in step S1: the real-time data at least comprises gas data output by a gas source point, hydraulic parameters and parameter data used for representing the physical properties of a pipe network, wherein:
the gas data includes a gas mole percentage at a gas source point;
the hydraulic parameters at least comprise one or a combination of a plurality of pressures, flow rates and temperatures of an air source point and a user point;
the parameter data at least comprises one or a combination of several of the height above sea level of the pipeline, the length of the pipeline and the diameter of the pipeline.
It should be noted that the gas data is data of gas output from each gas source in the pipe network, and is mainly affected by the production time, gas reservoir conditions, seabed pressure, temperature and other environmental factors, so that the gas data changes in real time at different gas sources and different times, and the gas data can be more suitable for the real-time state of the gas in the pipe network by using the output gas batch division at the gas sources.
In the embodiment, during actual operation, before batch division, the on-site SCADA system acquires data such as online chromatographic components, flow and pressure of each platform and each terminal, and transmits the data to the real-time database through the industrial personal computer, and the pipe network simulation software extracts real-time boundary data from the real-time database as input and control conditions of the online simulation model of the pipeline system to drive the simulation model to perform online simulation of the pipeline system.
In one case of this embodiment, in step S2, the gas output from the gas source point is divided into batches, and the step of determining a plurality of target gas batches includes:
s201: the method comprises the steps of performing batch division on gas groups output by each gas source point according to a pipe network static simulation result, and determining a plurality of initial gas batches;
s202: acquiring current mixed data of gas batches actually generated in the pipe network according to the data corresponding to each initial gas batch;
in actual operation, the current mixing data for the actually generated gas batch includes at least one or more of a minimum rate of change, a fluid update time, and a batch volume flow rate, wherein:
minimum rate of change Δ: defined as the sum of the absolute values of the differences between the source point mole percent of each gas component in a gas batch and the source point mole percent of each gas component in a previous gas batch, expressed by the formula:
Figure BDA0003857576530000071
in which n denotes the number of components i in the gas batch, C i Represents the gas source point mole percentage of component i in the gas batch, c i Represents the gas source point mole percentage of component i in the previous batch;
maximum fluid update time T: defining the time step length of the latest gas batch actually generated in the pipeline at the gas source and the time step length of the last batch;
minimum batch volume V: the volume of the working condition of a new batch of actually generated gas in the pipeline at the gas source is defined, and the expression formula is as follows:
V=ν·t·A
wherein ν is the gas flow velocity at the judgment position, t is the fluid update time step length at the judgment position, and A is the pipeline cross-sectional area at the judgment position.
S203: the current mixing data of the gas batch actually generated in step S202 is determined, and the target gas batch is determined according to the determination result.
It should be noted that, because the initial gas batch is set according to the pipe network static simulation result in the simulation calculation process, the main purpose is to set the initial gas batch for determining the condition of actually generating the gas batch in the pipe network, and the batch actually generated by the fluid in the pipe network at the gas source is according to the actual gas flow, so that in the production of the actual gas batch, the initial gas batch set at each gas source point is needed to obtain three determination parameters of the gas batch, the actual generation batch of the fluid in the pipe network needs to be satisfied, and the three parameters reach the preset range, so as to serve as the condition for determining the gas batch as the target gas batch.
Specifically, in step S203, the step of determining the current mixing data of each gas batch acquired in step S202 and determining the target gas batch according to the determination result includes:
s2031: respectively giving a set value for the minimum component change rate, the minimum batch volume and the maximum fluid update time: Δ set, V set, and T set;
s2032: comparing the calculated values of the minimum component change rate Δ, the minimum batch volume V, and the maximum fluid update time T of the gas batch obtained in step S202 with the set values determined in step S2031 one by one;
s2033: judging whether the comparison result of the step S2032 meets the conditions that delta is greater than or equal to delta setting in real time, V is greater than or equal to V setting in real time, and T is greater than or equal to T setting in real time;
if so, determining the gas batch as a target gas batch, otherwise determining the previous gas batch as the target gas batch.
In the embodiment, in actual operation, three preset values need to be given to the three determination conditions respectively, then three parameters obtained by actually calculating the gas batches are compared with the three set values one by one, if the determination conditions are met, it represents that the output gas group at the gas source generates a target gas batch, and the volume flow of the gas batch and the corresponding components of different gases are greatly changed, so that the gas batch is determined as the target gas batch, and the target gas batch determination of the pipe network gas is more in line with the actual fluid delivery condition.
Example 2
In this embodiment, an improvement is made on the basis of the gas tracking method provided in embodiment 1, and in step S3, the obtaining of the calculation result of each target gas batch at the mixing point specifically includes:
s301, calculating a first product of the gas source point mole percentage of each component gas in the target gas batch and the volume flow of the target gas batch at the mixing point, and calculating a first sum of the first products corresponding to each component gas in the target gas batch;
s302: calculating a second sum of the volumetric flow rates of each component gas in the target gas batch at the mixing point;
s303: respectively taking the quotient of the first sum value and the second sum value corresponding to each component gas to obtain the mole percentage of each component gas at the mixing point of the target gas batch;
s304: and determining the calculation result of the target gas batch according to the gas mole percentage of the target gas batch.
In practical application, the obtained target gas batch is tracked and calculated at the mixing point, and the mol percentages of different component gases of the target gas batch at the mixing point can be calculated by using the gas source point mol percentage of each component gas in the target gas batch and the volume flow at the mixing point, so as to calculate the real-time information of the target gas batch in the pipeline.
In actual operation, the gas component content at the mixing point can also be solved by a similar method, because at least two branch lines are crossed and converged at the mixing point, the nearest target gas batch in the branch line in front of the mixing point is utilized, and the target gas batch in the branch line is equal to a plurality of gas sources, so that the mole percentage of the gas components at the mixing point is calculated, the calculation result is more accurate, the mixed natural gas component content is calculated according to the real-time components and the volume flow of each gas source and the flow proportion of each gas source, and the mixing calculation is as follows:
Figure BDA0003857576530000091
in the formula, C j Is the mole percentage of component j at the mixing point; c. C i,j Is the mole percent of the ith gas source component j; qi is the volumetric flow rate of the ith gas source component at the mixing point; l is the number of gas sources mixed at the mixing point.
It should be noted that, in simulation application, fluid media of different gas sources are different, so that components in input pipelines of the gas sources are different, a pipe section is filled with batches of corresponding gas sources before a mixing point in an initial stage of a pipe network, and the pipe section is filled with the mixed batches after the mixing point until a next gas mixing point; the composition of the gas in the same batch is consistent, and only after the composition of the target gas batch is calculated, the position change information of the target gas batch at the next time step or the predicted arrival time of the batch in a certain pipe section can be calculated.
Specifically, after the mole percentage of the target gas batch is obtained, the mole percentage can be used to calculate other data of the gas in the pipe network, so as to obtain a calculation result of the gas in the pipe network, where the calculation result at least includes one or more of a calorific value, position change information, and arrival time of the target batch of the gas, where:
determining the calculation formula of the target batch gas heat value according to the gas mole percentage in the target gas batch as follows:
Figure BDA0003857576530000092
wherein Q is the calorific value of the gas of the target batch, Q j Is the calorific value, C, of component j in the target batch of gas j Is the mole percentage of the gas component j in the target gas batch, and n is the number of the components in the target gas batch;
according to the gas mole percentage in the target gas batch, the calculation formula for determining the position change information of the target batch gas is as follows:
Figure BDA0003857576530000093
in the formula (I), the compound is shown in the specification,
Figure BDA0003857576530000094
indicating the position where the target lot moves at the next time step,
Figure BDA0003857576530000095
position, v, representing the initial time of the target batch i The flow velocity of the gas in the unit cell of the pipe section is shown, and delta t represents a time step;
determining the target batch arrival time by the calculation formula according to the gas mole percentage in the target gas batch:
Figure BDA0003857576530000101
where t represents the arrival time of the batch, t 0 Representing the generation time of a certain target batch, and deltax representing the space step of the pipeline; v. of i Representing the flow rate of the gas in the section unit cell;
as shown in fig. 2, in step S4: during actual operation, a dynamic tracking model can be established according to the gas components, the gas heat values and the gas position change information of the target gas batch in the step S3, and then the target gas components, the target gas heat values and the target gas position change information corresponding to the time point to be inquired are determined through the dynamic tracking model;
in actual operation, as shown in fig. 2, the abscissa represents the upstream and downstream positions of the pipeline, the ordinate represents the time, and the diagonally divided regions represent the natural gas batches in the pipeline; and dynamically simulating the positions of different batches in the pipeline at each time step from the starting time, and advancing the batches downstream at the gas flow rate calculated by water power so as to obtain the real-time information of the gas in the pipeline network.
Based on the above embodiments, the following experiments are combined to verify the effects of the present invention:
as shown in fig. 3, the multi-gas source pipe network includes 1 trunk line, 2 branch lines, which are connected to 3 gas sources and 1 user terminal, the components of the output gas of each gas source are different, and the requirements of the user terminal on the gas components are also different;
(1) The data acquisition process comprises the following steps: the pipe network is provided with a pipe network online simulation system, the field SCADA system collects online chromatographic components, flow, pressure and the like of each platform/terminal and transmits the data to a real-time database through an industrial personal computer, and pipe network simulation software extracts real-time boundary data from the database to serve as input and control conditions of a pipeline system online simulation model and drives the simulation model to realize the online simulation of the pipeline system;
since the gas source component data is real-time data of SCADA, 3 gas source components at the initial time are shown in FIG. 4;
(2) The batch dividing process comprises the following steps: first, determination values of three conditions for batch determination are set: minimum change rate Δ =0.006, minimum batch volume V =50 cubic meters, maximum fluid update time t =120s, so the actual target gas batch is determined under the following conditions: when the change rate of the components of the natural gas in 120s is more than 0.006 and the volume of the natural gas continuously flowing through is more than 50 cubic meters, judging that a new batch is generated;
as shown in fig. 5, from 14 to 21 months 8, the flow of the gas source and the pressure of the user are from the SCADA system real-time data; calculating to obtain that the pipeline newly generates 7 batches in a time period, and can sequentially obtain the gas components of each batch;
(3) And (3) tracking calculation process:
as shown in fig. 6, heat value calculation is performed on 7 batches of natural gas based on the component conditions of each batch, and the heat value of each batch is obtained through calculation;
as shown in fig. 7, based on the composition of each batch, calculating the position of the next time step of the gas batch according to the gas flow rate and the time step, determining the flow state of the target gas batch in the pipeline, and calculating the time of the target gas batch to reach each node, wherein the node includes a gas source point and a mixing point, so as to obtain the specific batch generation time and the specific batch arrival time;
(4) Acquiring real-time information of gas of a point to be inquired:
determining the components, the heat value and the position change information of the target gas corresponding to the time point to be inquired through a dynamic tracking model; and calculating the position of the next time step gas batch according to the gas flow rate and the time step.
As shown in fig. 8, in one embodiment, end user simulated and measured high/low heating value arrival times and heating value ranges are compared; the difference between the simulated high-order heat value arrival time and the actually-measured arrival time of the embodiment is 10min, and the time tracking precision can reach more than 98% by considering that the time from the generation of the high-order heat value from the mixing point to the transportation to the user is 9.16 h; the difference between the simulated heat value and the actual measured heat value is 5Btu/scf at most, and the tracking precision of the heat value is more than 99 percent by considering that the average actual measured heat value is 960 Btu/scf; the accurate tracking of the heat values of the components of the offshore multi-gas-source pipe network can be realized, and the requirements of field technicians on pipe network system scheduling are met.
Example 3
The embodiment provides a gas tracking calculation device, which is applied to the gas tracking calculation method in a multi-gas-source pipe network provided in any one of embodiments 1 and 2, and the method includes:
the data acquisition module is used for acquiring the real-time data of the current gas pipe network and transmitting the acquired real-time data to the batch division module;
the batch division module is used for carrying out batch division on the gas groups at the gas source points according to the step S2 and uploading the obtained target gas batches to the tracking calculation module;
the tracking calculation module is used for performing tracking calculation on the target gas batch according to the step S3 and transmitting the obtained calculation result to the information display module;
the information display module is in signal connection with the data acquisition module, the batch division module and the tracking calculation module, and records and displays the acquired or generated data;
an information input module: the tracking calculation module is arranged on the information display module, is connected with the batch division module and the tracking calculation module, and is used for inputting data and inquiring real-time information of gas in the time node pipe network.
In practical application, the data acquisition module is arranged to conveniently acquire real-time data of offshore multi-gas-source pipes, the batch division module performs batch division of gas groups by using the real-time data as input and control conditions of an online simulation model of a pipeline system, tracking calculation is performed on hydraulic simulation calculation and component tracking calculation of natural gas pipe network flow according to target gas batches, and component data, heat value information and batch position information of different batches of natural gas are output; real-time gas information in the pipe network is displayed through the information display module, and the information input module inputs time nodes needing to be inquired, so that the real-time gas information in the pipe network can be inquired in real time; therefore, the accuracy of fluid arrival time tracking and heat value tracking is effectively improved, the gas tracking quality can be obviously improved, and the timeliness of gas tracking is improved.
Example 4
The embodiment is an embodiment of a server provided by the present invention, and includes:
a memory for storing a computer program;
and the processor is used for realizing the steps of any one of the gas tracking calculation methods in the multi-gas source pipe network when executing the computer program.
Referring to fig. 9, in practical operation, the server includes: the system comprises a processor, a memory, a bus and a communication interface, wherein the processor, the communication interface and the memory are connected through the bus; a processor for executing executable modules, such as computer programs, stored in memory; in particular, the memory is used for storing a program, and the processor executes the program after receiving an execution instruction, and the method performed by the apparatus defined in the foregoing embodiments may be applied to or implemented by the processor.
Illustratively, the memory may comprise a high-speed Random Access Memory (RAM), and may also comprise a non-volatile memory, such as at least one disk memory, wherein the communication connection between the system network element and at least one other network element is implemented via at least one communication interface (which may be wired or wireless), and may use the internet, a wide area network, a local area network, a metropolitan area network, etc.
Further, the bus may be an ISA bus, a PCI bus, an EISA bus, etc., and the bus may be divided into an address bus, a data bus, a control bus, etc., which are indicated by a single double-headed arrow in fig. 9 for convenience of illustration, but do not indicate only one bus or one type of bus.
Further, the processor may be an integrated circuit chip having signal processing capabilities; in implementation, the steps of the above method may be performed by integrated logic circuits of hardware in a processor or by instructions in the form of software.
Specifically, the processor may be a general-purpose processor, including a Central Processing Unit (CPU), a Network Processor (NP), and the like, and may also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), an off-the-shelf programmable gate array (FPGA) or other programmable logic device, a discrete gate or transistor logic device, or a discrete hardware component, and may implement or execute each method, step, and logic block disclosed in the embodiments of the present invention.
Example 5
This embodiment is an embodiment of a storage medium provided by the present invention, where the storage medium stores a computer program, and when the computer program is executed by a processor, the method implements the steps of any one of the above-mentioned methods for calculating gas tracking in a multi-gas source pipe network.
The computer program product of the storage medium provided by the embodiment of the present invention includes a computer readable storage medium storing program code, the program code including instructions for executing the method described in the previous method embodiment.
The steps of the method disclosed by the embodiment of the invention can be directly embodied as the execution of a hardware decoding processor, or the combination of hardware and software modules in the decoding processor; the software module can be located in a storage medium mature in the field such as a random access memory, a flash memory, a read only memory, a programmable read only memory or an electrically erasable programmable memory, a register, etc., the storage medium is located in a memory, and a processor reads information in the memory and combines hardware thereof to complete the steps of the method.
If the method of the embodiment of the invention is realized in the form of a software functional module and is sold or used as an independent product, the method can be stored in a computer readable storage medium; based on such understanding, the technical solutions of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the methods according to the embodiments of the present invention.
Illustratively, the storage medium package may be various media that can store program codes, such as a U-disk, a removable hard disk, a Read Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk.
The method, the device, the server and the storage medium for calculating the gas tracking in the multi-gas-source pipe network are described in detail above; the principle and the implementation of the present application are explained herein by applying specific examples, and the above description of the embodiments is only used to help understand the method and the core idea of the present application; it should be noted that, for those skilled in the art, without departing from the principle of the present application, the present application can also make several improvements and modifications, and those improvements and modifications also fall into the protection scope of the claims of the present application.

Claims (10)

1. A gas tracking calculation method in a multi-gas-source pipe network is characterized by comprising the following steps:
s1: acquiring real-time data of all boundary nodes of the gas pipe network; the boundary nodes of the gas pipe network comprise a gas source point and a user point;
s2: according to the real-time data in the step S1, dividing the gas output from the gas source point in batches to determine a plurality of target gas batches;
s3: according to the target gas batches determined in the step S2, tracking calculation is carried out on each target gas batch at a mixing point, and a calculation result of each target gas batch is obtained;
s4: acquiring real-time gas information in the gas pipe network at the time point to be inquired according to the calculation result of the target gas batch in the step S3;
wherein the real-time gas information at least comprises one or more of target gas components, target gas heat value and target gas position change information.
2. The method for tracking and calculating the gas in the multi-gas-source pipe network according to claim 1, wherein in the step S1: the real-time data at least comprises gas data output by a gas source point, hydraulic parameters and parameter data used for representing the physical properties of a pipe network, wherein:
the gas data includes a gas mole percentage at a gas source point;
the hydraulic parameters at least comprise one or a combination of a plurality of pressures, flow rates and temperatures of an air source point and a user point;
the parameter data at least comprises one or a combination of several of the height above sea level of the pipeline, the length of the pipeline and the diameter of the pipeline.
3. The method according to claim 2, wherein in step S2, the gas output from the gas source point is divided into batches, and the step of determining a plurality of target gas batches comprises:
s201: dividing the gas groups output from each gas source point in batches according to the static simulation result of the pipe network, and determining a plurality of initial gas batches;
s202: acquiring current mixed data of gas batches actually generated in the pipe network according to the data corresponding to each initial gas batch;
s203: the current mixing data of the gas batch actually generated in step S202 is determined, and the target gas batch is determined according to the determination result.
4. The method according to claim 3, wherein in step S202, the actually generated current mixture data of the gas batch at least includes one or more of a minimum change rate, a fluid update time and a batch volume flow rate, wherein:
minimum rate of change Δ: defined as the sum of the absolute values of the differences in the gas source point mole percent of each gas component in the gas batch and the gas source point mole percent of each gas component in the previous gas batch, expressed by the formula:
Figure FDA0003857576520000021
in which n denotes the number of components i in the gas batch, C i Represents the gas source point mole percent of component i in the gas batch, c i Represents the gas source point mole percentage of component i in the previous batch;
maximum fluid update time T: defining the time step length of the latest gas batch actually generated in the pipeline at the gas source and the time step length when the last batch is generated;
minimum batch volume V: the volume of the working condition of the new batch of gas actually generated in the pipeline at the gas source is defined, and the expression formula is as follows:
V=ν·t·A
wherein ν is the gas flow velocity at the judgment position, t is the fluid update time step length at the judgment position, and A is the cross-sectional area of the pipeline at the judgment position.
5. The method according to claim 4, wherein in step S203, the step of determining the current mixed data of each gas batch obtained in step S202 and determining the target gas batch according to the determination result comprises:
s2031: respectively giving a set value for minimum component change rate, minimum batch volume and maximum fluid update time: Δ set, V set, and T set;
s2032: comparing the calculated values of the minimum component change rate Δ, the minimum batch volume V, and the maximum fluid update time T of the gas batch obtained in step S202 with the set values determined in step S2031 one by one;
s2033: judging whether the comparison result of the step S2032 meets the conditions that delta is greater than or equal to delta setting in real time, V is greater than or equal to V setting in real time, and T is greater than or equal to T setting in real time;
if so, determining the gas batch as a target gas batch, otherwise determining the previous gas batch as the target gas batch.
6. The method according to claim 1, wherein in step S3, the obtaining of the calculation result for each target gas batch at the mixing point specifically comprises:
s301, calculating a first product of the gas source point mole percentage of each component gas in the target gas batch and the volume flow of the target gas batch at the mixing point, and calculating a first sum of the first products corresponding to each component gas in the target gas batch;
s302: calculating a second sum of the volumetric flow rates of each component gas in the target gas batch at the mixing point;
s303: respectively carrying out quotient on the first sum value and the second sum value corresponding to each component gas to obtain the mole percentage of each component gas of the target gas batch at the mixing point;
s304: determining the calculation result of the target gas batch according to the gas mole percentage of the target gas batch.
7. The method as claimed in claim 6, wherein in step S304, the calculation result at least includes one or more of a calorific value, a position change information and an arrival time of the target batch of gas, wherein:
determining the calculation formula of the target batch gas heat value according to the gas mole percentage in the target gas batch as follows:
Figure FDA0003857576520000031
wherein Q is the calorific value of the gas of the target batch, Q j Is the calorific value, C, of component j in the target batch of gas j Is the mole percentage of the gas component j in the target gas batch, and n is the number of the components in the target gas batch;
according to the gas mole percentage in the target gas batch, the calculation formula for determining the position change information of the target batch gas is as follows:
Figure FDA0003857576520000032
in the formula (I), the compound is shown in the specification,
Figure FDA0003857576520000033
indicating the position where the target lot moves at the next time step,
Figure FDA0003857576520000034
indicating the position of the initial time of the target batch, v i The flow velocity of the gas in the unit cell of the pipe section is shown, and delta t represents a time step;
determining the target batch arrival time by the calculation formula according to the gas mole percentage in the target gas batch:
Figure FDA0003857576520000035
where t represents the arrival time of the batch, t 0 Representing the generation time of a certain target batch, and deltax representing the space step of the pipeline; v. of i Indicating the flow rate of the gas in the unit cell of the pipe section.
8. A gas tracking calculation apparatus applied to the gas tracking calculation method according to any one of claims 1 to 7, comprising:
the data acquisition module is used for acquiring the real-time data of the current gas pipe network and transmitting the acquired real-time data to the batch division module;
the batch division module is used for carrying out batch division on the gas groups at the gas source points according to the step S2 and uploading the obtained target gas batches to the tracking calculation module;
the tracking calculation module is used for performing tracking calculation on the target gas batch according to the step S3 and transmitting the obtained calculation result to the information display module;
the information display module is in signal connection with the data acquisition module, the batch division module and the tracking calculation module, and records and displays the acquired or generated data;
an information input module: the tracking calculation module is arranged on the information display module, is connected with the batch division module and the tracking calculation module, and is used for inputting data and inquiring real-time information of gas in the time node pipe network.
9. A server, comprising:
a memory for storing a computer program;
a processor for implementing the gas tracking calculation method of any one of claims 1 to 7 when executing the computer program.
10. A storage medium having stored thereon a computer program which, when executed by a processor, implements a gas tracking calculation method as claimed in any one of claims 1 to 7.
CN202211152701.8A 2022-09-21 2022-09-21 Gas tracking calculation method and device in multi-gas-source pipe network, server and storage medium Pending CN115630479A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202211152701.8A CN115630479A (en) 2022-09-21 2022-09-21 Gas tracking calculation method and device in multi-gas-source pipe network, server and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211152701.8A CN115630479A (en) 2022-09-21 2022-09-21 Gas tracking calculation method and device in multi-gas-source pipe network, server and storage medium

Publications (1)

Publication Number Publication Date
CN115630479A true CN115630479A (en) 2023-01-20

Family

ID=84902424

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211152701.8A Pending CN115630479A (en) 2022-09-21 2022-09-21 Gas tracking calculation method and device in multi-gas-source pipe network, server and storage medium

Country Status (1)

Country Link
CN (1) CN115630479A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116608420A (en) * 2023-07-21 2023-08-18 新仟意能源科技(成都)集团有限责任公司 Dynamic tracking and monitoring method and system for natural gas components

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN116608420A (en) * 2023-07-21 2023-08-18 新仟意能源科技(成都)集团有限责任公司 Dynamic tracking and monitoring method and system for natural gas components
CN116608420B (en) * 2023-07-21 2023-09-19 新仟意能源科技(成都)集团有限责任公司 Dynamic tracking and monitoring method and system for natural gas components

Similar Documents

Publication Publication Date Title
CN110108328B (en) Method for acquiring water leakage amount of leakage area of water supply pipe network
CN109344436B (en) Online simulation method for large complex natural gas pipe network system
CN106682369A (en) Heating pipe network hydraulic simulation model identification correction method and system, method of operation
CN115630479A (en) Gas tracking calculation method and device in multi-gas-source pipe network, server and storage medium
CN107014451A (en) The method of ultrasonic flow sensor coefficient is speculated based on generalized regression nerve networks
CN114491864B (en) Nuclear power pipe network model preprocessing method with parameterization and reconfigurable characteristics
CN111125931B (en) Automatic optimization type searching method and device for cooling tower
JP5020338B2 (en) Initial value generating apparatus and initial value generating method
Wang et al. An Improved Grey Prediction Model Based on Matrix Representations of the Optimized Initial Value.
CN113239645A (en) System and method for dynamically tracking and calculating natural gas components in natural gas pipe network in real time
CN108958029A (en) For parameter tuning method and system
CN103324989A (en) Urban water supply pipe network optimizing method
CN116305696A (en) Urban gas pipe network non-isothermal steady-state hydraulic calculation method, equipment, storage medium and program product based on Broyden iteration method
CN111022937A (en) Water pipe network leakage positioning system and positioning method
CN112182905B (en) Heat supply pipe network simulation method and device for comprehensive energy system
Mehta et al. Simulation of existing water distribution network at Punagam Area of Surat City using WaterGEMS software
CN117933658A (en) Gas system planning method and device, electronic equipment and storage medium
CN110795810A (en) Online hydraulic model generation method
Wang et al. Simulation analysis of the effects of friction on the performance of pressure regulator for drip tape
CN104123442B (en) A kind of online Equivalent Model modeling method of Wind turbines active power
CN114718662B (en) Control method and device for regulating valve of zero-output system of low-pressure cylinder of thermal power generating unit
CN111693423B (en) Goaf permeability coefficient inversion method based on genetic algorithm
CN102261933B (en) A kind of large user's water-use model defining method
Matiko et al. Computer-aided system for designing the thermal energy metering devices
CN117874977A (en) Offline transient simulation calculation method based on autonomous switching calculation type function

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination