CN105302979A - Modeling method and system of valve groups in two-phase fluid network model - Google Patents

Abstract

The invention relates to a modeling method and system of valve groups in a two-phase fluid network model. The modeling method comprises the steps that the fluid density and the current valve position of each valve of each valve group are obtained; the current flow capacity of each valve is calculated according to the fluid density and the current valve position of each valve of each valve group; the comprehensive flow capacity of each valve group is calculated according to the current flow capacity of each valve of each valve group; a branch grid of each valve group is established in a simulation system according to the comprehensive flow capacity, wherein the branch grid comprises two nodes and a branch connected between the two nodes and equivalent to the valve group. Due to the fact that each valve group is divided into one branch grid and the two nodes and the comprehensive flow capacity of the branch is obtained by adopting an algebraic method according to the fluid density and the current valve position of each valve of each valve group, the number of the nodes is decreased, and the calculation amount can be decreased. In addition, the number of the nodes is decreased, data points needed to be monitored on site are also decreased, and accordingly the cost is reduced.

Description

The modeling method of valve sets and system in two-p hase fluid network model

Technical field

The present invention relates to Power Plant Simulation field, particularly relate to modeling method and the system of valve sets in two-p hase fluid network model.

Background technology

Along with the high speed development of China's power industry, the unit of Large Copacity, high parameter becomes main flow, analogue system is also slowly transferred to the vital role that can be and generate and provide prioritization scheme and Instructing manufacture by original pure training system, that is current High Precision Simulation system is slowly slowly drawn close to production by backstage, produce the actual production at scene generating and more directly affect, this just has higher requirement to the modeling of replicating machine and computational accuracy.And along with in-circuit emulation technology be following development trend, this has just had higher requirement to the rationality of simulation modeling.

In the modeling of two-p hase fluid network model, the division of branch road grid is basis, sets up two-p hase fluid network model according to bypass flow.The division of existing branch road grid is using any one valve arrangement as a branch road, and equipment and equipment connection place as internal node process, then go to calculate every bar bypass flow.Due to branch road, node is more, and on the one hand, calculated amount is comparatively large, on the other hand, needs to detect to obtain node data and branch data to each node and branch road at the scene, and need to be equipped with a large amount of spot measurement device, cost is higher.

Summary of the invention

Based on this, be necessary to provide a kind of cost low and the modeling method of valve sets and system in the two-p hase fluid network model that calculated amount is little.

A modeling method for valve sets in two-p hase fluid network model, comprising:

Obtain the current valve position of each valve of fluid density and described valve sets;

The current negotiability of described each valve is calculated according to the current valve position of each valve of described fluid density and described valve sets;

The combined flow ability of described valve sets is calculated according to the current negotiability of each valve of described valve sets;

In analogue system, set up the branch road grid of described valve sets according to described combined flow ability, described branch road grid comprises two nodes and is connected to the branch road being equivalent to described valve sets between described two nodes.

Wherein in a kind of embodiment, the formula calculating the current negotiability of described each valve according to the current valve position of each valve of described fluid density and described valve sets is:

Cv＝CvMAX*Vp*Vp*ρ

Wherein, Cv is current negotiability, and Vp is the current valve position of valve, and ρ is fluid density, and CvMAX is valve is maximum flow ability.

Wherein in a kind of embodiment, the combined flow ability that the current negotiability of described each valve according to described valve sets calculates described valve sets comprises:

According to the annexation between valve described in described valve sets, the flow using the current negotiability of described each valve as described each valve place branch road, calculates the combined flow ability of described valve sets respectively by algebraic method.

Wherein in a kind of embodiment, the step of the current valve position of each valve of described acquisition fluid density and described valve sets is: the current valve position obtaining each valve of described fluid density and the described valve sets monitored, or, obtain the current valve position of the described fluid density of user's input and each valve of described valve sets.

Wherein in a kind of embodiment, described analogue system is DCOSE.

The present invention also provides the modeling of valve sets in a kind of two-p hase fluid network model, comprising:

Acquisition module, for obtaining the current valve position of each valve of fluid density and described valve sets;

First computing module, the current valve position for each valve according to described fluid density and described valve sets calculates the current negotiability of described each valve;

Second computing module, the current negotiability for each valve according to described valve sets calculates the combined flow ability of described valve sets;

MBM, for setting up the branch road grid of described valve sets in analogue system according to described combined flow ability, described branch road grid comprises two nodes and is connected to the branch road being equivalent to described valve sets between described two nodes.

Wherein in a kind of embodiment, the formula calculating the current negotiability of described each valve according to the current valve position of each valve of described fluid density and described valve sets is:

Cv＝CvMAX*Vp*Vp*ρ

Wherein, Cv is current negotiability, and Vp is the current valve position of valve, and ρ is fluid density, and CvMAX is valve is maximum flow ability.

Wherein in a kind of embodiment, described second computing module, specifically according to the annexation between valve described in described valve sets, flow using the current negotiability of described each valve as described each valve place branch road, calculates the combined flow ability of described valve sets respectively by algebraic method.

Wherein in a kind of embodiment, described acquisition module, specifically for obtaining the current valve position of each valve of described fluid density and the described valve sets monitored, or, the current valve position of the described fluid density inputted specifically for acquisition user and each valve of described valve sets.

Wherein in a kind of embodiment, described analogue system is DCOSE.

The modeling method of valve sets in this two-p hase fluid network model, the current negotiability of each valve is calculated by the current valve position of each valve of the fluid density that obtains and valve sets, the combined flow ability of valve sets is calculated again according to the current negotiability of each valve of valve sets, to set up the branch road grid of valve sets in analogue system according to combined flow ability, each valve sets is divided into a branch road grid.The modeling method of this two-p hase fluid network model, because each valve sets is divided into a branch road grid, two nodes, and branch road combined flow ability adopts algebraic approach to calculate according to the current valve position of fluid density and valve, therefore, number of nodes reduces thus calculated amount is reduced, and on the other hand, number of nodes reduces, need the data point of monitoring also to reduce at the scene, thus reduce cost.

Accompanying drawing explanation

Fig. 1 is the process flow diagram of the modeling method of valve sets in a kind of two-p hase fluid network model;

Fig. 2 is a kind of schematic diagram of typical valve sets;

Fig. 3 is the grid schematic diagram adopting classic method partitioning site-branch road;

Fig. 4 is the grid schematic diagram of the method partitioning site-branch road adopting the application;

Fig. 5 is the module diagram of the modeling of valve sets in a kind of two-p hase fluid network model.

Embodiment

As shown in Figure 1, the modeling method of valve sets in a kind of two-p hase fluid network model, comprising:

S10: the valve position obtaining each valve of fluid density and valve sets.

Concrete, the valve position of fluid density and each valve can be inputted by host computer by user, or this analogue system is connected with on-the-spot pick-up unit, in real time the detection data at scene is sent to analogue system.

S30: the current negotiability calculating each valve according to the current valve position of each valve of fluid density and valve sets.

The current negotiability of valve is relevant with the current valve position of fluid density and valve, calculates the current negotiability of valve according to the current valve position of fluid density and valve.

S50: the combined flow ability calculating valve sets according to the current negotiability of each valve of valve sets.

S70: the branch road grid setting up valve sets according to combined flow ability in analogue system, branch road grid comprises two nodes and is connected to the branch road being equivalent to valve sets between two nodes.

In analogue system, set up the branch road grid of valve sets according to combined flow, each valve sets is divided into a branch road grid.

The modeling method of valve sets in this two-p hase fluid network model, the current negotiability of each valve is calculated by the current valve position of each valve of the fluid density that obtains and valve sets, the combined flow ability of valve sets is calculated again according to the current negotiability of each valve of valve sets, to set up the branch road grid of valve sets in analogue system according to combined flow ability, each valve sets is divided into a branch road grid.The modeling method of valve sets in this two-p hase fluid network model, because each valve sets is divided into a branch road grid, two nodes, and branch road combined flow ability adopts algebraic approach to calculate according to the current valve position of fluid density and valve, therefore, number of nodes reduces thus calculated amount is reduced, and on the other hand, number of nodes reduces, need the data point of monitoring also to reduce at the scene, thus reduce cost.

Concrete, in step S30, the formula calculating the current negotiability of each valve according to the current valve position of each valve of fluid density and valve sets is:

Cv＝CvMAX*Vp*Vp*ρ

Wherein, Cv is current negotiability, and Vp is the current valve position of valve, and ρ is fluid density, and CvMAX is valve is maximum flow ability.

Wherein in a kind of embodiment, step S50 comprises: according to the annexation in valve sets between valve, the flow using the current negotiability of each valve as each valve place branch road, is calculated the combined flow ability of valve sets by algebraic method respectively.

In a particular embodiment, this analogue system is DCOSE (distributed emulation back-up environment), and it comprises the foundation of real time database service, in real time computing engines, Computer Aided Modeling system, operator terminal and runs four major parts.Wherein Computer Aided Modeling system comprises again 5 solutions: the graphic configuration of fluid network automatic modeling, control loop configuration, electrical network automatic modeling, ladder diagram modeling and general fashion.

The modeling method of this two-p hase fluid network model, is applicable to valve sets, and such as, in power plant system water pitch valve sets, boiler water filling pitch valve sets etc. on oxygen-eliminating device, valve sets comprises at least two branch roads in parallel, each branch road comprises at least one valve.

For Fig. 2 typical case valve sets, be described.

In traditional fluidic network theory, be using any one equipment as a branch road, equipment and equipment connection place are as internal node process.If by the typical valve sets partitioning site-branch road grid in Fig. 2,6 nodes and 7 branch roads can be divided into, specifically as shown in Figure 3.

Adopt the modeling method of the application, each valve sets is divided into a branch road grid and two nodes, as shown in Figure 4.

Suppose that the fluid of each node is incompressible fluid, then valve 1, valve 2 are identical with the flow of valve 3, and valve 4, valve 5 are identical with the flow of valve 6.

According to the formula of the current negotiability of each valve, the current negotiability of each valve can be calculated:

Cv ₁＝CvMAX ₁*Vp ₁*Vp ₁*ρ；

Cv ₂＝CvMAX ₂*Vp ₂*Vp ₂*ρ；

Cv ₃＝CvMAX ₃*Vp ₃*Vp ₃*ρ；

Cv ₄＝CvMAX ₄*Vp ₄*Vp ₄*ρ；

Cv ₅＝CvMAX ₅*Vp ₅*Vp ₅*ρ；

Cv ₆＝CvMAX ₆*Vp ₆*Vp ₆*ρ；

Cv ₇＝CvMAX ₇*Vp ₇*Vp ₇*ρ；

Wherein, lower target 1,2,3,4,5,6,7 is number corresponding with valve respectively.

By the Meshing Method in Fig. 3, then have 6 nodes, be respectively node 1 ~ node 6, its pressure is P1, P2, P3, P4, P5, P6 respectively, has 7 branch roads, is respectively branch road 1 ~ branch road 7.

For branch road 1, wherein containing equipment valve 1, its upper node is No. 1 node lower node is No. 3 nodes, then the flow flow through in valve 1 can be expressed as:

F ₁ ²＝CvMAX ₁*Vp ₁*Vp ₁*ρ*(P3-P1)(1)

Accordingly, the flow that other valve flows through can be expressed as:

F ₂ ²＝CvMAX ₂*Vp ₂*Vp ₂*ρ*(P4-P3)(2)

F ₃ ²＝CvMAX ₃*Vp ₃*Vp ₃*ρ*(P2-P4)(3)

F ₄ ²＝CvMAX ₄*Vp ₄*Vp ₄*ρ*(P5-P1)(4)

F ₅ ²＝CvMAX ₅*Vp ₅*Vp ₅*ρ*(P6-P5)(5)

F ₆ ²＝CvMAX ₆*Vp ₆*Vp ₆*ρ*(P2-P6)(6)

F ₇ ²＝CvMAX ₇*Vp ₇*Vp ₇*ρ*(P2-P1)(7)

Wherein: subscript 1,2,3,4,5,6,7 represents corresponding every bar branch road; F is bypass flow; CvMAX is this branch road maximum current capacity; Vp is this valve position; ρ is the fluid density of circulation.

For in the mathematical model of Fig. 3,

Then by formula (1), (2), (3) are cumulative, abbreviation arranges and can obtain:

$\frac{{F_1}^2}{{\mathrm{CV}}_1}+\frac{{F_2}^2}{{\mathrm{CV}}_2}+\frac{{F_3}^2}{{\mathrm{CV}}_3}=(P3-P1)+(P4-P3)+(P2-P4)---\left(8\right)$

For in the mathematical model of Fig. 3, suppose that the fluid of each node is incompressible fluid, then have:

F ₁＝F ₂＝F ₃

Therefore: $F_1^2=(P2-P1)/(\frac{1}{{\mathrm{CV}}_1}+\frac{1}{{\mathrm{CV}}_2}+\frac{1}{{\mathrm{CV}}_3})---\left(9\right)$

In like manner, the 4th, 5,6 bypass flow can in like manner be calculated as

$F_4^2=(P2-P1)/(\frac{1}{{\mathrm{CV}}_4}+\frac{1}{{\mathrm{CV}}_5}+\frac{1}{{\mathrm{CV}}_6})$

Article 7, bypass flow is

$F_7^2=(P2-P1)*{\mathrm{CV}}_7---\left(10\right)$

Total flow F between node 1 and node 2 is

F＝F ₁+F ₄+F ₇

Then, in Fig. 3 mathematical model, node 1 with the total flow of node 2 is:

$F=(\sqrt{\frac{1}{\frac{1}{{\mathrm{CV}}_1}+\frac{1}{{\mathrm{CV}}_2}+\frac{1}{{\mathrm{CV}}_3}}}+\sqrt{\frac{1}{\frac{1}{{\mathrm{CV}}_4}+\frac{1}{{\mathrm{CV}}_5}+\frac{1}{{\mathrm{CV}}_6}}}+\sqrt{{\mathrm{CV}}_7})*\sqrt{(P2-P1)}---\left(11\right)$

For in the mathematical model of Fig. 3, the flow F between node 1 and node 2 is

F ²＝Cv*ρ*(P2-P1)(12)

Wherein: F is bypass flow, CV is the comprehensive through-current capability of each branch road, and ρ is fluid density.

And for Fig. 3 and Fig. 4, its physical object is all the same, so its total flow is also identical, so arrange formula 12 and formula 11, the branch road combined flow ability of Fig. 4 can be obtained:

$CV=(\sqrt{\frac{1}{\frac{1}{{\mathrm{CV}}_1}+\frac{1}{{\mathrm{CV}}_2}+\frac{1}{{\mathrm{CV}}_3}}}+\sqrt{\frac{1}{\frac{1}{{\mathrm{CV}}_4}+\frac{1}{{\mathrm{CV}}_5}+\frac{1}{{\mathrm{CV}}_6}}}+\sqrt{{\mathrm{CV}}_7})---\left(13\right)$

Wherein, the current negotiability of each valve equals the combined flow ability of every bar branch road of Fig. 3, as can be seen here, the branch road combined flow ability of the application is solved by algebraic method, instead of solve bypass flow and node pressure by building math matrix, thus reduce the treatment capacity of data, improve modeling efficiency.

The present invention also provides the modeling of valve sets in a kind of two-p hase fluid network model, as shown in Figure 5, comprising:

Acquisition module 10, for obtaining the current valve position of each valve of fluid density and described valve sets.

Concrete, the valve position of fluid density and each valve can be inputted by host computer by user, or this analogue system is connected with on-the-spot pick-up unit, in real time the detection data at scene is sent to analogue system.Accordingly, acquisition module, the current valve position of the described fluid density inputted specifically for acquisition user and each valve of described valve sets, or, specifically for obtaining the current valve position of each valve of described fluid density and the described valve sets monitored.

First computing module 30, the current valve position for each valve according to described fluid density and described valve sets calculates the current negotiability of each valve.

The current negotiability of valve is relevant with the current valve position of fluid density and valve, calculates the current negotiability of valve according to the current valve position of fluid density and valve.

Second computing module 50, the current negotiability for each valve according to described valve sets calculates the combined flow ability of valve sets.

MBM 70, for setting up the branch road grid of valve sets in analogue system according to combined flow ability, branch road grid comprises two nodes and is connected to the branch road being equivalent to valve sets between two nodes.

In analogue system, set up the branch road grid of valve sets according to combined flow, each valve sets is divided into a branch road grid.

The modeling of valve sets in this two-p hase fluid network model, the current negotiability of each valve is calculated by the current valve position of each valve of the fluid density that obtains and valve sets, the combined flow ability of valve sets is calculated again according to the current negotiability of each valve of valve sets, to set up the branch road grid of valve sets in analogue system according to combined flow ability, each valve sets is divided into a branch road grid.The modeling of valve sets in this two-p hase fluid network model, because each valve sets is divided into a branch road grid, two nodes, and branch road combined flow ability adopts algebraic approach to calculate according to the current valve position of fluid density and valve, therefore, number of nodes reduces thus calculated amount is reduced, and on the other hand, number of nodes reduces, need the data point of monitoring also to reduce at the scene, thus reduce cost.

Concrete, the formula calculating the current negotiability of each valve according to the current valve position of each valve of fluid density and valve sets is:

Cv＝CvMAX*Vp*Vp*ρ

Wherein, Cv is current negotiability, and Vp is the current valve position of valve, and ρ is fluid density, and CvMAX is valve is maximum flow ability.

Wherein in a kind of embodiment, second computing module 50, specifically for according to the annexation in valve sets between valve, the flow using the current negotiability of each valve as each valve place branch road, calculates the combined flow ability of valve sets respectively by algebraic method.

In a particular embodiment, this analogue system is DCOSE (distributed emulation back-up environment), and it comprises the foundation of real time database service, in real time computing engines, Computer Aided Modeling system, operator terminal and runs four major parts.Wherein Computer Aided Modeling system comprises again 5 solutions: the graphic configuration of fluid network automatic modeling, control loop configuration, electrical network automatic modeling, ladder diagram modeling and general fashion.

The modeling of this two-p hase fluid network model, is applicable to valve sets, and such as, in power plant system water pitch valve sets, boiler water filling pitch valve sets etc. on oxygen-eliminating device, valve sets comprises at least two branch roads in parallel, each branch road comprises at least one valve.

Each technical characteristic of the above embodiment can combine arbitrarily, for making description succinct, the all possible combination of each technical characteristic in above-described embodiment is not all described, but, as long as the combination of these technical characteristics does not exist contradiction, be all considered to be the scope that this instructions is recorded.

The above embodiment only have expressed several embodiment of the present invention, and it describes comparatively concrete and detailed, but can not therefore be construed as limiting the scope of the patent.It should be pointed out that for the person of ordinary skill of the art, without departing from the inventive concept of the premise, can also make some distortion and improvement, these all belong to protection scope of the present invention.Therefore, the protection domain of patent of the present invention should be as the criterion with claims.

Claims (10)

1. the modeling method of valve sets in two-p hase fluid network model, is characterized in that, comprising:

Obtain the current valve position of each valve of fluid density and described valve sets;

The current negotiability of described each valve is calculated according to the current valve position of each valve of described fluid density and described valve sets;

The combined flow ability of described valve sets is calculated according to the current negotiability of each valve of described valve sets;

In analogue system, set up the branch road grid of described valve sets according to described combined flow ability, described branch road grid comprises two nodes and is connected to the branch road being equivalent to described valve sets between described two nodes.

2. the modeling method of valve sets in two-p hase fluid network model according to claim 1, is characterized in that, the formula calculating the current negotiability of described each valve according to the current valve position of each valve of described fluid density and described valve sets is:

Cυ＝CvMAX*Vp*Vp*ρ

Wherein, C υ is current negotiability, and Vp is the current valve position of valve, and ρ is fluid density, and CvMAX is valve is maximum flow ability.

3. the modeling method of valve sets in two-p hase fluid network model according to claim 1, is characterized in that, the combined flow ability that the current negotiability of described each valve according to described valve sets calculates described valve sets comprises:

According to the annexation between valve described in described valve sets, the flow using the current negotiability of described each valve as described each valve place branch road, calculates the combined flow ability of described valve sets respectively by algebraic method.

4. the modeling method of valve sets in two-p hase fluid network model according to claim 1, it is characterized in that, the step of the current valve position of each valve of described acquisition fluid density and described valve sets is: the current valve position obtaining each valve of described fluid density and the described valve sets monitored, or, obtain the current valve position of the described fluid density of user's input and each valve of described valve sets.

5. the modeling method of valve sets in two-p hase fluid network model according to claim 1, it is characterized in that, described analogue system is DCOSE.

6. the modeling of valve sets in two-p hase fluid network model, is characterized in that, comprising:

Acquisition module, for obtaining the current valve position of each valve of fluid density and described valve sets;

First computing module, the current valve position for each valve according to described fluid density and described valve sets calculates the current negotiability of described each valve;

Second computing module, the current negotiability for each valve according to described valve sets calculates the combined flow ability of described valve sets;

MBM, for setting up the branch road grid of described valve sets in analogue system according to described combined flow ability, described branch road grid comprises two nodes and is connected to the branch road being equivalent to described valve sets between described two nodes.

7. the modeling of valve sets in two-p hase fluid network model according to claim 6, is characterized in that, the formula calculating the current negotiability of described each valve according to the current valve position of each valve of described fluid density and described valve sets is:

Cυ＝CvMAX*Vp*Vp*ρ

Wherein, C υ is current negotiability, and Vp is the current valve position of valve, and ρ is fluid density, and CvMAX is valve is maximum flow ability.

8. the modeling of valve sets in two-p hase fluid network model according to claim 1, it is characterized in that, described second computing module, specifically according to the annexation between valve described in described valve sets, flow using the current negotiability of described each valve as described each valve place branch road, calculates the combined flow ability of described valve sets respectively by algebraic method.

9. the modeling of valve sets in two-p hase fluid network model according to claim 1, it is characterized in that, described acquisition module, specifically for obtaining the current valve position of each valve of described fluid density and the described valve sets monitored, or, the current valve position of the described fluid density inputted specifically for acquisition user and each valve of described valve sets.

10. the modeling of valve sets in two-p hase fluid network model according to claim 1, it is characterized in that, described analogue system is DCOSE.