CN102779199A - Method for realizing dynamic simulation process of natural gas pipe network conveying system - Google Patents

Method for realizing dynamic simulation process of natural gas pipe network conveying system Download PDF

Info

Publication number
CN102779199A
CN102779199A CN2011101246518A CN201110124651A CN102779199A CN 102779199 A CN102779199 A CN 102779199A CN 2011101246518 A CN2011101246518 A CN 2011101246518A CN 201110124651 A CN201110124651 A CN 201110124651A CN 102779199 A CN102779199 A CN 102779199A
Authority
CN
China
Prior art keywords
hydraulic power
therrmodynamic
solution
dynamic simulation
power system
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
CN2011101246518A
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.)
Petrochina Co Ltd
Original Assignee
Petrochina Co Ltd
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 Petrochina Co Ltd filed Critical Petrochina Co Ltd
Priority to CN2011101246518A priority Critical patent/CN102779199A/en
Publication of CN102779199A publication Critical patent/CN102779199A/en
Pending legal-status Critical Current

Links

Images

Landscapes

  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method for realizing a dynamic simulation process of a natural gas pipe network conveying system. Relates to the technical field of pipeline systems. The method comprises the following steps: the frog leap strategy separates a hydraulic system and a thermodynamic system; dealing with thermodynamic system instability in a 'windward' format; solving a linear equation system; and obtaining a solving result of the dynamic simulation numerical model of the natural gas pipe network conveying system. The coupling solution of the hydraulic system and the thermodynamic system is stable, short in calculation time and high in solution efficiency.

Description

Gas distributing system induction system dynamic simulation realization process method
Technical field
The present invention is a kind of gas distributing system induction system dynamic simulation realization process method.Relate to the piping system technical field.
Background technology
The emulation of gas distributing system is to utilize the mathematical model that rock gas flows in pipeline to describe; It is one group of nonlinear partial differential equations; Academia does not provide the analytic solution that meet the practical applications needs, in engineering, adopts numerical method to find the solution usually.Pipe network emulation comprises static Simulation and transient state emulation two aspects, and static Simulation is the strong instrument that carries out design of pipe networks, is the important evidence of choosing pipe network structure and parameter scheme, and its accuracy directly influences whether safe and reliable and economical rationality of design proposal.Transient state emulation can simulative tube network operation parameter over time; Can be directed against the different processes design proposal; For example peak regulation, pipeline emptying and fracture burst accident carry out transient state and calculate, according to result of calculation; Propose good solution, thereby ensure the efficient scheduling management and in time find and handle burst accident.Gas distributing system induction system dynamic simulation compares with static Simulation that to find the solution difficulty bigger, and range of application is wider.In 2005 " University Of Tianjin's paper " " theoretical analysis of gas distribution system, gas network security and emulation and applied research " to disclose a kind of serve as the steady dynamic simulation theoretical model of gas ductwork under the gentle non-isothermal condition such as basis foundation etc. with three big conservation equations; The result compares and analyzes through the contrast gained; Draw the applicability of its engineering, and then reference is provided for the optimization of gas ductwork.
But at present gas distributing system induction system dynamic simulation procedural implementation method exists waterpower, the unsettled difficult problem of therrmodynamic system Solving Coupled, and computing time is long, it is low to find the solution efficient.
The method of finding the solution the dynamic simulation numerical model commonly used comprises approximate analytical method (linearization) and numerical method two big classes.The processing of approximate analytical method is with the linear partial linearization of the branch of mathematical model, changes into partial differential equations, obtains analytic solution according to initial value then, and linearizing method mainly contains average, proper that interpolate value and least square method.Approximate analytical method mainly be a little that computing velocity is fast, shortcoming is that computational accuracy is relatively poor, is applicable to simulation pipe network operation situation in the long-time scope.Numerical method is the maturation of Along with computer technology and growing up gradually, mainly is method of finite difference, and the principle of this method is to change into difference equation to PDE, finds the solution difference equation with iterative numerical approach then.PDE change into method of Difference Equations and mainly contain method of characteristic, implicit expression center method and Taylor expansion method.The advantage of numerical solution is that computational accuracy is high and calculate more flexibly, and shortcoming is that to find the solution speed fast not as the approximate analysis solution.
For method of characteristic; Non-linear algebraic equation group to each node and the pipe interior net point of system are independently set up corresponding low dimension can turn to the PDE on the pipeline characteristic difference equation, is characterized in need not finding the solution huge Nonlinear System of Equations; Be easy to find the solution; Take less calculator memory, but in order to satisfy the stability of finding the solution, the time level is often obtained very little.Yet; For the gas distributing system system, because the compressibility of gas, it is much littler than the unstable degree that fluid pipeline occurs that it unstable degree occurs; Thereby the time level got for a short time and there is no need, and do not meet the requirement of practical applications to the dynamic simulation time step.
For implicit difference method, consider all nodes and pipe interior net point, set up unified high dimensional nonlinear Algebraic Equation set.Since will be with treating that unknown quantity simultaneous all on the seeking time level gets up to find the solution, the Nonlinear System of Equations that need find the solution is very huge, and finding the solution needs the long time.But this method can guarantee the stability found the solution, and time step can obtain bigger.
Summary of the invention
The objective of the invention is to invent that a kind of waterpower, therrmodynamic system Solving Coupled are stable, computing time is short, find the solution the high gas distributing system induction system dynamic simulation realization process method of efficient.
Gas distributing system induction system dynamic simulation procedural implementation method (see figure 1) of the present invention is:
The strategy that " leapfrogs " separates waterpower, therrmodynamic system;
The format analysis processing of " facining the wind " therrmodynamic system instability;
Solving Linear;
Obtain gas distributing system induction system dynamic simulation numerical model solving result.
" leapfroging " strategy of said Hydraulic Power System and therrmodynamic system is:
Many Hydraulic Power System and therrmodynamic system are coupled of the whole bag of tricks of finding the solution the dynamic simulation numerical model at present commonly used carried out simultaneous solution; In the actual emulation process; The acute variation of the outer boundary of system (pressure, temperature) can cause that serious " sawtooth " phenomenon appears in temperature curve in the solution procedure; Cause the concussion of pressure and flow distribution curve; Even the temperature computation value of physical significance appears not meeting, until result of calculation disperse, system crash.
For " sawtooth " of treatment temperature curve and to the interference problem of Hydraulic Power System to therrmodynamic system; Whole gas distributing system system is divided into Hydraulic Power System and therrmodynamic system is discrete respectively and simultaneous; Each the time step find the solution separately respectively; " extrapolated value " in step imported Hydraulic Power System into and found the solution when earlier thermal parameter being taked preceding two, and the Hydraulic Power System solving result is passed to therrmodynamic system, carries out finding the solution of therrmodynamic system.
It finds the solution the flow process (see figure 2):
Hydraulic Power System and therrmodynamic system are discrete respectively;
Hydraulic Power System and therrmodynamic system are set up original state respectively;
The thermal parameter extrapolated value of therrmodynamic system is handled;
The Hydraulic Power System original state is set up the back and is handled the back acting in conjunction in Hydraulic Power System with the thermal parameter extrapolated value;
Arrive therrmodynamic system behind the Hydraulic Power System; Therrmodynamic system arrives the thermal parameter extrapolated value again and handles;
Finish behind the Hydraulic Power System.
The strategy that " leapfrogs " is alternately found the solution Hydraulic Power System and therrmodynamic system in time; Each energy equation obtains pressure with the front and flow value is the basic calculation temperature, and the temperature value that two time steps " extrapolation " obtain before each Hydraulic Power System is basic calculation pressure and flow; Its advantage replaces a large matrix to find the solution with two less Matrix Solving of scale, has improved computing velocity significantly; Because the special construction of minor matrix can be stored with sparse matrix technology; Although be not subject matter, off-line dynamic simulation numerical model is found the solution still have crucial meaning but fully improve counting yield for the modern computing motor speeds;
Said " facining the wind " format analysis processing therrmodynamic system instability (see figure 3) is:
The instability that basic reason that " sawtooth " appear in curve is therrmodynamic system appears in temperature when Solving Coupled Hydraulic Power System and therrmodynamic system; After " leapfroging " strategy separation waterpower, therrmodynamic system; When increasing computing velocity; Also unsettled therrmodynamic system is separated and focused on, the convergence exponent number of therrmodynamic system on time shaft brought up to 2 rank from 1 rank, greatly improved the stability of therrmodynamic system.
In addition, when the pipeline inner fluid flow direction changed, temperature curve " sawtooth " can occur equally, thereby causes the collapse of Hydraulic Power System.This is because in finding the solution the dynamic simulation system process; The very important point is speed and a direction of confirming that flow field information is propagated; So that numeric format is consistent with the physical property that flows; Otherwise the result who obtains sizable concussion can occur near rough zone, even the result who obtains is unacceptable.In order to address the above problem, introduced " facining the wind " form; Upstreame scheme is to confirm the difference orientation according to the characteristic velocity direction in flow field; Meet the perturbation wave propagation law physically; This obtains very little restriction for satisfied stability requires time horizon to have overcome method of characteristic simultaneously, and adopts central difference to cause the deficiency of solving result concuss because of ignoring the flowing information direction of propagation.
In gas distributing system induction system dynamic simulation numerical model is found the solution specifically; Adopt the piece-wise constant function to approach temperature curve; In the direction of each computing node, decide the interpolation form (regressive interpolation and forward interpolation) of the derivative of temperature to guarantee it is to carry out interpolation processing all the time along the mobile direction of fluid according to flow velocity; This method has solved the stability problem at the heat equation that flows to not timing effectively.
Said Solving Linear is topological structure (see figure 4) and the characteristic according to pipe network; NLS to forming is carried out topological sorting; Make that the Jacobian matrix that forms is the block diagonal angle sparse matrix of a band edge; Adopt the Harwell-Boeing form to carry out carrying out exchanges data with the Solving Linear device, inner pretreated consuming time when finding the solution with minimizing, improve on the whole and find the solution efficient; Practice shows, for the matrix of 30,000 dimensions, and its linear solution consuming time at 10 Milliseconds.
The present invention proposes is applicable to and has overcome waterpower, the unsettled difficult problem of therrmodynamic system Solving Coupled by gas distributing system induction system dynamic simulation numerical model method for solving, shortens computing time, improves and find the solution efficient, satisfies application requirements in the actual production.
Description of drawings
Fig. 1 dynamic simulation procedural implementation method process flow diagram
Fig. 2 tactful synoptic diagram that " leapfrogs "
Fig. 3 tactful synoptic diagram that facings the wind
Fig. 4 pipe network topology diagram
Embodiment
Embodiment. embodiment of the present invention is described and the present invention is further described with this example.This example is that gas distributing system induction system dynamic simulation numerical model is found the solution the efficient instance, and its flow process is as shown in Figure 1:
The strategy that " leapfrogs " separates waterpower, therrmodynamic system;
The format analysis processing of " facining the wind " therrmodynamic system instability;
Solving Linear;
Obtain gas distributing system induction system dynamic simulation numerical model solving result.
For the counting yield of verifying the gas distributing system dynamic simulation procedural implementation method that the present invention proposes and the accuracy of finding the solution; As embodiment, carried out dynamic simulation result contrast with domestic certain gas distributing system with foreign commercial pipeline simulation software Stoner Pipeline Simulator (SPS) (gas component, ground temperature distribute, pressure is set off in first battle, terminal enters the station pressure, compressor and valve state) under same condition.
Long defeated gas distributing system comprises 3 sources of the gas, 7 compressor stations, 6 defeated stations of branch, 32 valve chambers, 45 pipeline sections, main line length 1568km among the embodiment.Be without loss of generality, provide the topological structure at one of them compressor station, 2 valve chambers, 1 defeated station of branch at this, as shown in Figure 4.
Begin to start the operation dynamic simulation from zero flow, as a child system was basicly stable when simulation time accumulative total reaches 16, and gas distributing system intake condition and outlet parameter result of calculation contrast situation are respectively shown in table 1 and table 2:
The contrast of table 1 embodiment pipe network start-up course emulation intake condition result of calculation
Figure BDA0000061108650000051
Table 2 embodiment pipe network start-up course outlet parameter simulation calculation result contrast
Figure BDA0000061108650000052
Annotate: S representes on the hurdle SPS result of calculation, and R representes that the present invention proposes implementation method result of calculation, and E representes relative error
Comparison result from table, the result of calculation goodness of fit of the result of calculation of this software and SPS is than higher, and wherein the average relative error of pressure is in 3%, and the average relative error of flow is in 3.5%, and the average relative error of temperature is in 1%.The counting yield aspect: in dominant frequency is 2.6GHZ, in save as on the PC of 2G, when computing node was 2000, the dynamic simulation single step was calculated average out to consuming time 0.210 second, was linear growth with interstitial content in the analogue system.
This example is through on probation, and it has overcome waterpower, the unsettled difficult problem of therrmodynamic system Solving Coupled, shortens computing time, improves and find the solution efficient, satisfies application requirements in the actual production.

Claims (4)

1. gas distributing system induction system dynamic simulation realization process method is characterized in that this method is:
The strategy that " leapfrogs " separates waterpower, therrmodynamic system;
The format analysis processing of " facining the wind " therrmodynamic system instability;
Solving Linear;
Obtain gas distributing system induction system dynamic simulation numerical model solving result.
2. gas distributing system induction system dynamic simulation realization process method according to claim 1 is characterized in that " leapfroging " strategy of said Hydraulic Power System and therrmodynamic system is:
Whole gas distributing system system is divided into Hydraulic Power System and therrmodynamic system is discrete respectively and simultaneous; Each the time step find the solution separately respectively; " extrapolated value " in step imported Hydraulic Power System into and found the solution when earlier thermal parameter being taked preceding two; The Hydraulic Power System solving result is passed to therrmodynamic system, carry out finding the solution of therrmodynamic system;
It finds the solution flow process:
Hydraulic Power System and therrmodynamic system are discrete respectively;
Hydraulic Power System and therrmodynamic system are set up original state respectively;
The thermal parameter extrapolated value of therrmodynamic system is handled;
The Hydraulic Power System original state is set up the back and is handled the back acting in conjunction in Hydraulic Power System with the thermal parameter extrapolated value;
Arrive therrmodynamic system behind the Hydraulic Power System; Therrmodynamic system arrives the thermal parameter extrapolated value again and handles;
Finish behind the Hydraulic Power System.
The strategy that " leapfrogs " is alternately found the solution Hydraulic Power System and therrmodynamic system in time; Each energy equation obtains pressure with the front and flow value is the basic calculation temperature, and the temperature value that two time steps " extrapolation " obtain before each Hydraulic Power System is basic calculation pressure and flow; Store with sparse matrix technology.
3. gas distributing system induction system dynamic simulation realization process method according to claim 1 is characterized in that said " facining the wind " format analysis processing therrmodynamic system instability is:
Confirm the difference orientation according to the characteristic velocity direction in flow field; In gas distributing system induction system dynamic simulation numerical model is found the solution; Adopt the piece-wise constant function to approach temperature curve; In the direction of each computing node according to flow velocity, decide the interpolation form of temperature derivative, the direction that flows along fluid is all the time carried out interpolation processing.
4. gas distributing system induction system dynamic simulation realization process method according to claim 1; It is characterized in that said Solving Linear is: according to the topological structure and the characteristic of pipe network; NLS to forming is carried out topological sorting; Make that the Jacobian matrix that forms is the block diagonal angle sparse matrix of a band edge, adopt the Harwell-Boeing form to carry out carrying out exchanges data with the Solving Linear device.
CN2011101246518A 2011-05-13 2011-05-13 Method for realizing dynamic simulation process of natural gas pipe network conveying system Pending CN102779199A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN2011101246518A CN102779199A (en) 2011-05-13 2011-05-13 Method for realizing dynamic simulation process of natural gas pipe network conveying system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN2011101246518A CN102779199A (en) 2011-05-13 2011-05-13 Method for realizing dynamic simulation process of natural gas pipe network conveying system

Publications (1)

Publication Number Publication Date
CN102779199A true CN102779199A (en) 2012-11-14

Family

ID=47124111

Family Applications (1)

Application Number Title Priority Date Filing Date
CN2011101246518A Pending CN102779199A (en) 2011-05-13 2011-05-13 Method for realizing dynamic simulation process of natural gas pipe network conveying system

Country Status (1)

Country Link
CN (1) CN102779199A (en)

Cited By (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103955186A (en) * 2014-04-22 2014-07-30 中国石油大学(北京) Natural gas pipe network pipe flow state parameter determining method and device
CN103970029A (en) * 2014-05-14 2014-08-06 中国石油大学(华东) Gas pipeline dynamic simulation method based on pulse response sequences
CN106707795A (en) * 2016-12-30 2017-05-24 贵州大学 Real-time simulation model and system of natural gas pipeline network and power grid coupling system
CN107977513A (en) * 2017-11-30 2018-05-01 北京石油化工学院 A kind of temperature predicting method when natural gas dynamic flows in gas distributing system based on route searching
CN109754109A (en) * 2017-11-08 2019-05-14 中国石油天然气股份有限公司 Method and device for determining natural gas pipe network operation scheme
CN111177892A (en) * 2019-12-11 2020-05-19 中电普信(北京)科技发展有限公司 Distributed simulation system

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4559610A (en) * 1983-05-04 1985-12-17 Southwest Research Corporation Gas pumping system analog
CN1184931A (en) * 1996-12-11 1998-06-17 唐秀家 Method and apparatus for detecting and positioning leakage of fluid transferring pipeline
CN101625071A (en) * 2009-08-07 2010-01-13 天津大学 Method for measuring and locating leakage of gas pipelines
CN201763373U (en) * 2010-08-05 2011-03-16 太原理工大学 Test simulator researching goaf gas migration under different vitalization conditions

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4559610A (en) * 1983-05-04 1985-12-17 Southwest Research Corporation Gas pumping system analog
CN1184931A (en) * 1996-12-11 1998-06-17 唐秀家 Method and apparatus for detecting and positioning leakage of fluid transferring pipeline
CN101625071A (en) * 2009-08-07 2010-01-13 天津大学 Method for measuring and locating leakage of gas pipelines
CN201763373U (en) * 2010-08-05 2011-03-16 太原理工大学 Test simulator researching goaf gas migration under different vitalization conditions

Non-Patent Citations (3)

* Cited by examiner, † Cited by third party
Title
PHILIP T.KEENAN: "COLLOCATION AND UPWINDING FOR THERMAL FLOW IN", 《INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS》 *
安家荣: "长距离输油管道设计计算软件的开发", 《油气储运》 *
柏毅: "川西地区天然气管网优化及仿真技术研究", 《中国博士学位论文全文数据库(电子期刊) 工程科技I辑》 *

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103955186A (en) * 2014-04-22 2014-07-30 中国石油大学(北京) Natural gas pipe network pipe flow state parameter determining method and device
CN103955186B (en) * 2014-04-22 2016-08-24 中国石油大学(北京) Gas distributing system pipe flow condition parameter determination method and device
CN103970029A (en) * 2014-05-14 2014-08-06 中国石油大学(华东) Gas pipeline dynamic simulation method based on pulse response sequences
CN103970029B (en) * 2014-05-14 2017-03-29 中国石油大学(华东) A kind of Gas Pipeline Dynamic Simulation method based on Least square estimation
CN106707795A (en) * 2016-12-30 2017-05-24 贵州大学 Real-time simulation model and system of natural gas pipeline network and power grid coupling system
CN106707795B (en) * 2016-12-30 2023-06-06 贵州大学 Real-time simulation model and system of natural gas pipe network and power grid coupling system
CN109754109A (en) * 2017-11-08 2019-05-14 中国石油天然气股份有限公司 Method and device for determining natural gas pipe network operation scheme
CN107977513A (en) * 2017-11-30 2018-05-01 北京石油化工学院 A kind of temperature predicting method when natural gas dynamic flows in gas distributing system based on route searching
CN107977513B (en) * 2017-11-30 2021-04-16 北京石油化工学院 Natural gas dynamic flow temperature prediction method based on path search
CN111177892A (en) * 2019-12-11 2020-05-19 中电普信(北京)科技发展有限公司 Distributed simulation system
CN111177892B (en) * 2019-12-11 2023-05-02 中电普信(北京)科技发展有限公司 Distributed simulation system

Similar Documents

Publication Publication Date Title
CN102779199A (en) Method for realizing dynamic simulation process of natural gas pipe network conveying system
CN101599104A (en) A kind of analogy method of blade flutter boundary of aviation turbine engine
CN108153998B (en) Full three-dimensional pneumatic optimization design method for centrifugal blower impeller
CN103955186A (en) Natural gas pipe network pipe flow state parameter determining method and device
CN102779198B (en) Method for establishing dynamic simulation system of natural gas pipe network conveying system
CN106227967A (en) Industrial steam turbine low-pressure stage group vane type line optimization method
CN103366067B (en) The method of simulation and prediction fluid torque-converter performance
CN102819651A (en) Simulation-based parameter optimizing method for precise casting process of single crystal turbine blade
CN110083968B (en) Compressor characteristic prediction method based on correction of gas seal leakage influence numerical model
CN105631125A (en) Aerodynamic-thermal-structural coupling analysis method based on reduced-order model
CN103455725A (en) Non-steady flow simulating method for pipe network system
CN105574288A (en) Method for designing water inlet conduit three-dimensional body flow surface of high-performance large-flow pump station
CN102663241A (en) Analog calculation method for transient power of water turbine under elastic water attack
CN101739718B (en) Parameter template-based corn leaf virtual simulation modeling method
CN106777731A (en) A kind of gas ductwork hydraulic calculation method
CN110110457B (en) Noise distribution prediction method for wind power plant with complex terrain
CN111551178A (en) Shortest path-based segmented track time planning method
CN112330493B (en) Energy system modeling and comprehensive analysis method, device and storage medium
CN114117819B (en) Steady-state simulation method for hot steam network
CN203703468U (en) Reducing mitre elbow
CN102521360B (en) Raster data transmission method and system
CN105160062B (en) A kind of same journey hydraulic pipeline check method
CN103441837A (en) Four-dimensional chaotic system with constant lyapunov exponent
CN104731761B (en) Gas distributing system emulation mode and device
CN103473422B (en) Based on the axial wheel Airfoil Design method of singular point distribution

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C12 Rejection of a patent application after its publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20121114