CN115221693A - Heating system observable state judgment method considering quasi-dynamic pipeline temperature - Google Patents

Heating system observable state judgment method considering quasi-dynamic pipeline temperature Download PDF

Info

Publication number
CN115221693A
CN115221693A CN202210733776.9A CN202210733776A CN115221693A CN 115221693 A CN115221693 A CN 115221693A CN 202210733776 A CN202210733776 A CN 202210733776A CN 115221693 A CN115221693 A CN 115221693A
Authority
CN
China
Prior art keywords
state
observable
measurement
equation
variable
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
CN202210733776.9A
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.)
South China University of Technology SCUT
Original Assignee
South China University of Technology SCUT
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 South China University of Technology SCUT filed Critical South China University of Technology SCUT
Priority to CN202210733776.9A priority Critical patent/CN115221693A/en
Publication of CN115221693A publication Critical patent/CN115221693A/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
    • G06F17/11Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2111/00Details relating to CAD techniques
    • G06F2111/04Constraint-based CAD
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2119/00Details relating to the type or aim of the analysis or the optimisation
    • G06F2119/08Thermal analysis or thermal optimisation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Mathematical Physics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • General Engineering & Computer Science (AREA)
  • Data Mining & Analysis (AREA)
  • Software Systems (AREA)
  • Databases & Information Systems (AREA)
  • Algebra (AREA)
  • Operations Research (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • Testing And Monitoring For Control Systems (AREA)

Abstract

The invention discloses a method for judging an observable state of a heat supply system by considering quasi-dynamic pipeline temperature, which comprises the following steps: s1, determining input parameters of an observable state judgment model of a heating system; s2, establishing a heat supply system observable state judgment model according to the input parameters; and S3, solving the established heat supply system observable state judgment model to obtain the observable state of the system. When the quasi-dynamic state of the pipeline temperature is considered, the propagation speed of the temperature depends on the water flow, and the time-lag characteristic exists in the propagation process of the temperature. At the moment, the thermodynamic equation of the heating system and the estimation result of the flow influence each other, which causes that the traditional observable state judgment method is difficult to be applied to the heating system. Aiming at the problem, the invention provides a method for judging the observable state based on multiple time-lag scenes, which can ensure that the acquired observable state is effective under the scenes.

Description

Heating system observable state judgment method considering quasi-dynamic pipeline temperature
Technical Field
The invention relates to the technical field of observable analysis of a heating system, in particular to a method for judging an observable state of the heating system by considering the quasi-dynamic state of pipeline temperature.
Background
In practical applications, it is necessary to perform state estimation on the heating system to monitor the operation state. Before performing state estimation, the observability of the system needs to be analyzed to obtain input measurements for the state estimation. The observable state decision problem, which is the first problem to be addressed in observability analysis, determines the set of states that the system can observe under the collected measurements. The existing observable state judgment method mainly adopts a numerical method, needs to determine the specific form of a measurement equation, and cannot directly solve the problem of observable state judgment of a heating system. In heating systems, the temperature is carried in pipes by water. Temperature propagation in a pipe has a time-lag characteristic due to the time required for water to flow from the head end of the pipe to the tail end of the pipe. When the quasi-dynamic characteristic of the temperature of the pipeline of the heat supply system is considered, the thermal dynamic equation of the heat supply system is related to a time-lag scene, and the time-lag scene is determined by the estimation result of the water flow, so that the thermal dynamic equation and the estimation result are mutually influenced, the exact form of the thermal dynamic equation cannot be determined when the observable state is judged, and the conventional numerical method cannot be directly applied to the heat supply system.
Disclosure of Invention
The invention aims to overcome the defects and shortcomings of the prior art and provides a method for judging the observable state of a heating system by considering the quasi-dynamic state of the pipeline temperature. The method solves the problem that the thermodynamic equation of the heating system and the estimation result influence each other, and can provide a reliable observable analysis result for the heating system.
In order to achieve the purpose, the technical scheme provided by the invention is as follows: a heat supply system observable state judgment method considering pipeline temperature quasi-dynamic state is disclosed, wherein the heat supply system is composed of a heat source, a heat load, a water supply pipeline and a water return pipeline, and comprises the following steps:
s1, determining input parameters of a heat supply system observable state judgment model, including a collected measurement equation set
Figure BDA0003714901110000011
State set of heating system
Figure BDA0003714901110000012
Set of considered time-lapse scenarios
Figure BDA0003714901110000014
Measuring equation j and correlation coefficient a of system state i under time-lag scene (gamma, phi) ij,(γ,φ) Measurement equation j is
Figure BDA0003714901110000013
Individual scenarios involving a union of states
Figure BDA0003714901110000021
Number n of heating system nodes Node Set of considered time periods
Figure BDA0003714901110000022
Set of time period t heating system flow continuity equations
Figure BDA0003714901110000023
S2, establishing a target function and a constraint condition of the observable state judgment model of the heating system according to the input parameters, and performing linearization processing on a nonlinear term in the constraint condition;
and S3, solving the processed observable state judgment model of the heat supply system by using a solver to obtain the observable state of the system.
Further, the step S1 includes:
determining a set of collected metrology equations
Figure BDA0003714901110000024
It is derived from real-time measurement acquired by a measuring instrument and a state equation of a system; the real-time measurement comprises pipeline flow measurement, node temperature measurement, heat source heat supply measurement and heat load heat consumption measurement; the system's equation of state includes flow continuity equation, loop pressure drop equation and pipeline temperature snapAn equation of state;
determining a set of states of a heating system
Figure BDA0003714901110000025
The temperature of all nodes of the heating system and the flow of all pipelines are contained;
determining a set of considered time-lapse scenarios
Figure BDA0003714901110000026
The flow estimation method can be obtained by calculating the flow estimation results of pseudo measurement and historical time periods of the pipeline flow;
determining a correlation coefficient a of a measurement equation j and a system state i under a time-lag scene (gamma, phi) ij,(γ,φ) Measurement equation j is
Figure BDA0003714901110000027
Individual scenarios involving a union of states
Figure BDA0003714901110000028
By using
Figure BDA0003714901110000029
Representing a system state set related to a measurement equation j under a time-lag scene (gamma, phi); if the system status i belongs to
Figure BDA00037149011100000210
The correlation coefficient is assigned a ij,(γ,φ) =1; if the system status i does not belong to
Figure BDA00037149011100000211
The correlation coefficient is assigned a ij,(γ,φ) =0;
Figure BDA00037149011100000212
Can pass through
Figure BDA00037149011100000213
And calculating to obtain:
Figure BDA00037149011100000214
determining the number n of heating system nodes Node Set of considered time periods
Figure BDA00037149011100000215
Set of flow continuity equations of heat supply system for time period t
Figure BDA00037149011100000216
They can be obtained by analyzing the system structure.
Further, the step S2 includes the steps of:
s201, establishing an objective function of the heat supply system observable state judgment model according to the input parameters, wherein the number of observable states is maximized:
Figure BDA00037149011100000217
in the formula, n i Is a variable from 0 to 1, if state i is observable, the variable is 1, otherwise the variable is 0;
s202, establishing a constraint condition of a heat supply system observable state judgment model according to the input parameters, wherein the constraint condition is ensured to be in
Figure BDA0003714901110000031
In each scene, a group of basic measurement equations consisting of related measurement equations can be found, wherein the related measurement equations refer to
Figure BDA0003714901110000032
The basic measurement equation refers to a group of independent equations capable of solving the observable state; the established constraints specifically include:
the dependent metrology equations involve observability constraints on states that ensure that the states involved in the dependent metrology equations are
Figure BDA0003714901110000033
In each sceneAre all observable:
Figure BDA0003714901110000034
in the formula, w j Is a 0-1 variable, if the measurement equation j is relevant, the variable is 1, otherwise the variable is 0;
the source constraint of the basic measurement equations ensures that all basic measurement equations are related measurement equations:
Figure BDA0003714901110000035
in the formula, v j,(γ,φ) Is a variable of 0-1, if the measurement equation j is a basic measurement equation under a time-lag scene (gamma, phi), the variable is 1, otherwise the variable is 0;
the independence constraint of the basic metrology equations, which ensures that all basic metrology equations involving multiple states are independent of each other:
Figure BDA0003714901110000036
the mapping constraint between the observable state and the basic measurement equation ensures that the observable state and the basic measurement equation have a one-to-one mapping relationship, which is a condition that the basic measurement equation can solve the observable state and needs to meet:
Figure BDA0003714901110000037
Figure BDA0003714901110000038
in the formula, y ij,(γ,φ) Is a variable of 0 to 1, if the state i and the measurement equation j have a one-to-one correspondence relationship in a time-lag scene (gamma, phi), the variable is 1, otherwise the variable is 0;
s203, carrying out linearization processing on the nonlinear terms in the constraint conditions; in the constraints (2), (6) there are non-linear terms multiplied by a plurality of decision variables, using
Figure BDA0003714901110000039
Represents the k-th nonlinear term:
Figure BDA00037149011100000310
in the formula (I), the compound is shown in the specification,
Figure BDA0003714901110000041
a set of decision variables referred to for the k-th nonlinear term; x is a radical of a fluorine atom l Is the l decision variable;
the linear treatment is carried out on the formula (7) to obtain:
Figure BDA0003714901110000042
Figure BDA0003714901110000043
in the formula (I), the compound is shown in the specification,
Figure BDA0003714901110000044
the number of decision variables involved for the k-th nonlinear term.
Further, the step S3 includes the steps of:
s301, solving the processed observable state judgment model of the heat supply system by using a solver;
s302, acquiring an observable state of the system; the solving result of the heat supply system observable state judgment model contains a variable n i Represents the observability of state i, if n i If =1, it is indicated that state i is observable.
Compared with the prior art, the invention has the following advantages and beneficial effects:
1. the method provided by the invention considers the quasi-dynamic characteristic of the pipeline temperature in the heating system, provides the observable state judgment method based on the multi-time-lag scene, and ensures that the acquired observable state is effective in the scenes. The method can solve the problem that the thermal dynamic equation and the estimation result in the heat supply system are mutually influenced, and provides a reliable observable state result while considering the dynamic characteristic of the heat supply system.
2. The traditional observable state judgment method needs to clearly determine the specific form of a measurement equation, and can only obtain an observable state result in a time-lag scene each time. Compared with a method for judging observable states one by one for each time lag scene, the method provided by the invention considers observable judgment problems under a plurality of time lag scenes, and has the advantages of less decision variables and constraint conditions for establishing the model, smaller calculation scale and higher solving efficiency.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Fig. 2 is a classification diagram of measurement equations in the observable state determination problem.
FIG. 3 is a schematic diagram of real-time measurements collected at various time intervals in case one.
FIG. 4 is a schematic diagram of the real-time measurements collected at each time interval in case two.
Fig. 5 is a schematic diagram of real-time measurements collected at each time interval in case three.
Detailed Description
The present invention will be described in further detail with reference to examples and drawings, but the present invention is not limited thereto.
As shown in fig. 1, this embodiment provides a method for determining an observable state of a heat supply system considering quasi-dynamic pipeline temperature, where the heat supply system is composed of a heat source, a heat load, a water supply pipeline, and a water return pipeline, and the method for determining an observable state specifically includes the following steps:
s1, determining input parameters of a heat supply system observable state judgment model, including a collected measurement equation set
Figure BDA0003714901110000051
State set of heating system
Figure BDA0003714901110000052
Set of considered time-lapse scenarios
Figure BDA0003714901110000053
Measuring equation j and correlation coefficient a of system state i under time-lag scene (gamma, phi) ij,(γ,φ) Measurement equation j is
Figure BDA0003714901110000054
Individual scenarios involving a union of states
Figure BDA0003714901110000055
Number n of heating system nodes Node Set of considered time periods
Figure BDA0003714901110000056
Set of flow continuity equations of heat supply system for time period t
Figure BDA0003714901110000057
The method comprises the following specific steps:
determining a set of collected metrology equations
Figure BDA0003714901110000058
The system is derived from real-time measurement acquired by a measuring instrument and a state equation of the system; the real-time measurement comprises pipeline flow measurement, node temperature measurement, heat source heat supply measurement and heat load heat consumption measurement; the state equation of the system comprises a flow continuity equation, a loop pressure drop equation and a pipeline temperature quasi-dynamic equation;
determining a set of states of a heating system
Figure BDA0003714901110000059
The temperature of all nodes of the heating system and the flow of all pipelines are contained;
determining a set of considered time-lapse scenarios
Figure BDA00037149011100000510
The flow estimation method can be obtained by calculating the flow estimation results of pseudo measurement and historical time periods of the pipeline flow;
determining a correlation coefficient a of a measurement equation j and a system state i under a time-lag scene (gamma, phi) ij,(γ,φ) Measurement equation j is in
Figure BDA00037149011100000511
Individual scenarios involving a union of states
Figure BDA00037149011100000512
By using
Figure BDA00037149011100000513
Representing a system state set related to a measurement equation j under a time-lag scene (gamma, phi); if the system status i belongs to
Figure BDA00037149011100000514
The correlation coefficient is assigned a ij,(γ,φ) =1; if the system status i does not belong to
Figure BDA00037149011100000515
The correlation coefficient is assigned a ij,(γ,φ) =0;
Figure BDA00037149011100000516
Can pass through
Figure BDA00037149011100000517
And calculating to obtain:
Figure BDA00037149011100000518
determining the number n of heating system nodes Node Set of considered time periods
Figure BDA00037149011100000519
Set of flow continuity equations of heat supply system for time period t
Figure BDA00037149011100000520
They can be obtained by analyzing the system structure.
S2, establishing a target function and a constraint condition of the observable state judgment model of the heating system according to the input parameters, and carrying out linearization treatment on a nonlinear term in the constraint condition, wherein the concrete steps are as follows:
s201, establishing an objective function of the heat supply system observable state judgment model, wherein the number of observable states is maximized:
Figure BDA0003714901110000061
in the formula, n i Is a variable from 0 to 1, if state i is observable, the variable is 1, otherwise the variable is 0;
s202, establishing constraint conditions of an observable state judgment model of the heating system, wherein the constraint conditions are ensured to be in
Figure BDA0003714901110000062
In each scene, a group of basic measurement equations consisting of related measurement equations can be found, wherein the related measurement equations refer to
Figure BDA0003714901110000063
The basic measurement equation refers to a set of independent equations (as shown in fig. 2) capable of solving the observable state; the constraints established specifically include:
the dependent metrology equation involves an observability constraint of the state that ensures that the state involved in the dependent metrology equation is at
Figure BDA0003714901110000064
Observable in each scene:
Figure BDA0003714901110000065
in the formula, w j Is a variable of 0-1, if the measurement equation j is relevant, the variable is 1, otherwise the variable is 0;
the source constraints of the basic measurement equations ensure that all the basic measurement equations are related measurement equations:
Figure BDA0003714901110000066
in the formula, v j,(γ,φ) The variable is 0-1, if the measurement equation j is a basic measurement equation under a time-lag scene (gamma, phi), the variable is 1, otherwise, the variable is 0;
the independence constraint of the basic metrology equations, which ensures that all basic metrology equations involving multiple states are independent of each other:
Figure BDA0003714901110000067
the mapping constraint between the observable state and the basic measurement equation ensures that the observable state and the basic measurement equation have a one-to-one mapping relationship, which is a condition that the basic measurement equation can solve the observable state and needs to satisfy:
Figure BDA0003714901110000068
Figure BDA0003714901110000069
in the formula, y ij,(γ,φ) Is a variable of 0 to 1, if the state i and the measurement equation j have a one-to-one correspondence relationship in a time-lag scene (gamma, phi), the variable is 1, otherwise the variable is 0;
s203, carrying out linearization processing on the nonlinear terms in the constraint conditions; in the constraints (2), (6) there are non-linear terms multiplied by a plurality of decision variables, using
Figure BDA0003714901110000071
Represents the k-th nonlinear term:
Figure BDA0003714901110000072
in the formula (I), the compound is shown in the specification,
Figure BDA0003714901110000073
a set of decision variables referred to for the k-th nonlinear term; x is a radical of a fluorine atom l Is the l decision variable;
the linearization processing of the formula (7) can obtain:
Figure BDA0003714901110000074
Figure BDA0003714901110000075
wherein the content of the first and second substances,
Figure BDA0003714901110000076
the number of decision variables involved for the k-th nonlinear term.
S3, solving the processed heat supply system observable state judgment model by using a solver to obtain the observable state of the system, wherein the method comprises the following specific steps:
s301, solving an observable state judgment model of the heating system by using a solver;
s302, acquiring an observable state of the system; the solving result of the heat supply system observable state judgment model contains a variable n i Represents the observability of state i, if n i =1, state i is observable.
For a better understanding of the effectiveness of the invention and the advantages thereof over the prior art, reference should be made to the following examples which are set forth to illustrate, but are not to be construed as limiting the scope of the invention.
And setting the performance of the observable state judgment model of the heat supply system by using different cases. All observable state judgment models are solved by a Gurobi solver.
1. Validity of observable state judgment result
The validity of the invention is verified by the observable state judgment result of case one. In case of a 4-node heating system, taking into account 3 time intervals, the real-time measurements collected by the system at each time interval are shown in fig. 3. In case one, the time lag coefficients of all pipelines are the same, and the time lag coefficient set is considered
Figure BDA0003714901110000077
Are { (0, 1), (1, 1), (0, 2), (1, 2), (2, 2) }. The inventors compared the present invention with
Figure BDA0003714901110000078
The results of the unobservable states in each scene illustrate the effectiveness of the invention. The detailed comparison results are shown in Table 1, where the last action is the result of the invention, T 3,t-2 Representing the temperature state of node 3, T, over time period T-2 3,t-1 Representing the temperature state of node 3 for time period t-1 and 'x' representing a non-observable state.
TABLE 1 results of unobservable states under different time-lag scenarios
Time-lag scene (gamma, phi) T 3,t-2 T 3,t-1
(0,1) x
(1,1) x
(0,2) x x
(1,2) x x
(2,2) x
Consider all of the above scenarios (invention) x x
As can be seen from Table 1, the unobservable state determined by the present invention is
Figure BDA0003714901110000081
The union of the unobservable states in each scene. It follows that the observable state found by the present invention is
Figure BDA0003714901110000082
The intersection of the observable states in each scene. Thus, the observable state found by the present invention is
Figure BDA0003714901110000083
Each scene is observable. In addition, for an unobservable state, which is obtained by the present invention, it is always possible to
Figure BDA0003714901110000084
To find a scene in such a way that the state is not observable. Therefore, the number of unobservable states found by the present invention is the smallest and the number of observable states found is the largest. These results demonstrate the effectiveness of the present invention.
2. Comparison with existing methods
The performance of the present invention is illustrated by the results of the observable state determinations in case two and case three. Case two and case three are 12-node heating systems. In case two, the number of time segments considered is 3, and the real-time measurements collected for each time segment are shown in fig. 4. By using
Figure BDA0003714901110000085
Representing the set of time lag coefficients considered by pipe b during time period t. Case two assumes that the symmetrical water supply pipeline and the water return pipeline have the same time lag coefficient, so the considered time lag coefficient set
Figure BDA0003714901110000086
Time lag coefficient of 5 groups related to time period t and 5 water supply pipelines
Figure BDA0003714901110000087
Figure BDA0003714901110000088
The Cartesian product of (A), which are all taken as { (1, 1), (1, 2), (2, 2) }.
Figure BDA0003714901110000089
A total of 243 time-lapse scenarios are included. In case three, the number of considered time periods is 6. The heating system is complete in measurement configuration from time t-5 to time t-3, and the system state can be observed; the real-time measurements collected by the system from time t-2 to time t are shown in FIG. 5. Case three assumes that the time lag coefficients of the symmetrical water supply pipeline and the water return pipeline are the same, so the considered time lag coefficient set
Figure BDA00037149011100000810
Time lag coefficients of 15 groups of water supply pipelines related to time period t-2 to time period t and 5 water supply pipelines
Figure BDA00037149011100000811
Figure BDA00037149011100000812
The cartesian products of (a) and their arrangement are shown in table 2.
Figure BDA00037149011100000813
A total of 216 time-lag scenarios are included.
Table 2 three time lag coefficient sets considered by water supply pipeline
Figure BDA0003714901110000091
The existing observable state judgment method needs to clearly determine the specific form of a measurement equation, so that an observable state result under a time-lag scene can be obtained each time. Under the solution strategy based on a multi-time-lag scene, the existing method can only analyze one by one
Figure BDA0003714901110000092
The observable state under each scene judges the problem and will
Figure BDA0003714901110000093
And taking the intersection of the observable states obtained in each scene as a final result. The inventors illustrate the advantages of the present invention by comparing the results of the present invention with the prior methods in case two and case three. Specific comparison results are shown in Table 3, T 8,t-2 Representing the temperature state of node 8, T, over time period T-2 9,t-2 Indicating the temperature state of node 9, T, over time period T-2 8,t-1 Representing the temperature state of node 8, T, over time period T-1 9,t-1 Indicating the temperature state of node 9, T, over time period T-1 8,t Indicating the temperature state of node 8 for time period t.
TABLE 3 comparison of results of different solving methods
Figure BDA0003714901110000094
As can be seen from Table 3, the method can obtain the same observable state judgment result as the existing method, and has shorter solving time. The existing method only can analyze the observable state judgment problem under different time-delay scenes one by one, but the method considers the observable state judgment problem under a plurality of time-delay scenes simultaneously, has fewer decision variables and constraint conditions for establishing the model, and has smaller scale for establishing the model, so the method has higher calculation efficiency than the existing method.
The above embodiments are preferred embodiments of the present invention, but the present invention is not limited to the above embodiments, and any other changes, modifications, substitutions, combinations, and simplifications which do not depart from the spirit and principle of the present invention should be construed as equivalents thereof, and all such changes, modifications, substitutions, combinations, and simplifications are intended to be included in the scope of the present invention.

Claims (4)

1. A heat supply system observable state judgment method considering pipeline temperature quasi-dynamic state is characterized by comprising a heat source, a heat load, a water supply pipeline and a water return pipeline, and comprises the following steps:
s1, determining input parameters of a heat supply system observable state judgment model, including a collected measurement equation set
Figure FDA0003714901100000011
State set of heating system
Figure FDA0003714901100000012
Set of considered time-lapse scenarios
Figure FDA0003714901100000013
Measuring equation j and correlation coefficient a of system state i under time-lag scene (gamma, phi) ij,(γ,φ) Measurement equation j is in
Figure FDA0003714901100000014
Each scenario involves a union of states
Figure FDA0003714901100000015
Number n of heating system nodes Node Set of considered time periods
Figure FDA0003714901100000016
Set of flow continuity equations of heat supply system for time period t
Figure FDA0003714901100000017
S2, establishing a target function and a constraint condition of the observable state judgment model of the heating system according to the input parameters, and performing linearization processing on a nonlinear term in the constraint condition;
and S3, solving the processed heat supply system observable state judgment model by using a solver to obtain the observable state of the system.
2. A method for determining an observable state of a heating system considering quasi-dynamic state of pipe temperature according to claim 1, wherein the step S1 includes:
determining a set of collected metrology equations
Figure FDA0003714901100000018
It is derived from real-time measurement acquired by a measuring instrument and a state equation of a system; the real-time measurement comprises pipeline flow measurement, node temperature measurement, heat source heat supply measurement and heat load heat consumption measurement; the state equation of the system comprises a flow continuity equation, a loop pressure drop equation and a pipeline temperature quasi-dynamic equation;
determining a set of states of a heating system
Figure FDA0003714901100000019
The temperature of all nodes of the heating system and the flow of all pipelines are contained;
determination testSet of filtered time-lag scenarios
Figure FDA00037149011000000110
The flow estimation method can be obtained by calculating the flow estimation results of pseudo measurement and historical time periods of the pipeline flow;
determining a correlation coefficient a of a measurement equation j and a system state i under a time-lag scene (gamma, phi) ij,(γ,φ) Measurement equation j is
Figure FDA00037149011000000111
Each scenario involves a union of states
Figure FDA00037149011000000112
By using
Figure FDA00037149011000000113
Representing a system state set related to a measurement equation j under a time-lag scene (gamma, phi); if the system status i belongs to
Figure FDA00037149011000000114
The correlation coefficient is assigned a ij,(γ,φ) =1; if the system status i does not belong to
Figure FDA00037149011000000115
The correlation coefficient is assigned a ij,(γ,φ) =0;
Figure FDA00037149011000000116
Can pass through
Figure FDA00037149011000000117
And calculating to obtain:
Figure FDA00037149011000000118
determining the number n of heating system nodes Node Set of considered time periods
Figure FDA00037149011000000119
Set of time period t heating system flow continuity equations
Figure FDA0003714901100000021
They can be obtained by analyzing the system structure.
3. A method for determining observable states of a heating system in consideration of quasi-dynamic pipeline temperature according to claim 1, wherein the step S2 comprises the steps of:
s201, establishing an objective function of the heat supply system observable state judgment model according to the input parameters, wherein the number of observable states is maximized:
Figure FDA0003714901100000022
in the formula, n i Is a variable from 0 to 1, if state i is observable, the variable is 1, otherwise the variable is 0;
s202, establishing a constraint condition of a heat supply system observable state judgment model according to the input parameters, wherein the constraint condition is ensured to be in
Figure FDA0003714901100000023
In each scene, a group of basic measurement equations consisting of related measurement equations can be found, wherein the related measurement equations refer to
Figure FDA0003714901100000024
The basic measurement equation refers to a group of independent equations capable of solving the observable state; the established constraints specifically include:
the dependent metrology equation involves an observability constraint of the state that ensures that the state involved in the dependent metrology equation is at
Figure FDA0003714901100000029
Observable in each scene:
Figure FDA0003714901100000025
in the formula, w j Is a variable of 0-1, if the measurement equation j is relevant, the variable is 1, otherwise the variable is 0;
the source constraints of the basic measurement equations ensure that all the basic measurement equations are related measurement equations:
Figure FDA0003714901100000026
in the formula, v j,(γ,φ) Is a variable of 0-1, if the measurement equation j is a basic measurement equation under a time-lag scene (gamma, phi), the variable is 1, otherwise the variable is 0;
the independence constraint of the basic metrology equations, which ensures that all basic metrology equations involving multiple states are independent of each other:
Figure FDA0003714901100000027
the mapping constraint between the observable state and the basic measurement equation ensures that the observable state and the basic measurement equation have a one-to-one mapping relationship, which is a condition that the basic measurement equation can solve the observable state and needs to meet:
Figure FDA0003714901100000028
Figure FDA0003714901100000031
in the formula, y ij,(γ,φ) Is a 0-1 variable, if state i is in a time-lag scenario (γ, φ)If the measurement equation j has a one-to-one correspondence relationship, the variable is 1, otherwise the variable is 0;
s203, carrying out linearization processing on the nonlinear terms in the constraint conditions; in the constraints (2), (6) there are non-linear terms multiplied by a plurality of decision variables, using
Figure FDA0003714901100000032
Represents the k-th nonlinear term:
Figure FDA0003714901100000033
in the formula (I), the compound is shown in the specification,
Figure FDA0003714901100000034
a set of decision variables referred to for the k-th nonlinear term; x is a radical of a fluorine atom l Is the l decision variable;
the linear treatment is carried out on the formula (7) to obtain:
Figure FDA0003714901100000035
Figure FDA0003714901100000036
in the formula (I), the compound is shown in the specification,
Figure FDA0003714901100000037
the number of decision variables involved for the k-th nonlinear term.
4. A method for determining observable states of a heating system in consideration of quasi-dynamic pipeline temperature according to claim 1, wherein the step S3 comprises the steps of:
s301, solving the processed observable state judgment model of the heat supply system by using a solver;
s302, obtainingTaking an observable state of the system; the solving result of the heat supply system observable state judgment model comprises a variable n i Represents the observability of state i if n i =1, state i is observable.
CN202210733776.9A 2022-06-27 2022-06-27 Heating system observable state judgment method considering quasi-dynamic pipeline temperature Pending CN115221693A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202210733776.9A CN115221693A (en) 2022-06-27 2022-06-27 Heating system observable state judgment method considering quasi-dynamic pipeline temperature

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202210733776.9A CN115221693A (en) 2022-06-27 2022-06-27 Heating system observable state judgment method considering quasi-dynamic pipeline temperature

Publications (1)

Publication Number Publication Date
CN115221693A true CN115221693A (en) 2022-10-21

Family

ID=83609802

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202210733776.9A Pending CN115221693A (en) 2022-06-27 2022-06-27 Heating system observable state judgment method considering quasi-dynamic pipeline temperature

Country Status (1)

Country Link
CN (1) CN115221693A (en)

Similar Documents

Publication Publication Date Title
CN110108328B (en) Method for acquiring water leakage amount of leakage area of water supply pipe network
CN108897286B (en) Fault detection method based on distributed nonlinear dynamic relation model
CN117195007B (en) Heat exchanger performance prediction method and system
CN116595327B (en) Sluice deformation monitoring data preprocessing system and method
CN111721834A (en) Cable partial discharge online monitoring defect identification method
CN117196353A (en) Environmental pollution assessment and monitoring method and system based on big data
CN114826543A (en) AI 0T-based steam jet pump parameter transmission system and method
CN115221693A (en) Heating system observable state judgment method considering quasi-dynamic pipeline temperature
WO2020119012A1 (en) Dynamic system static gain estimation method based on historical data ramp response
CN117421951A (en) River pollution tracing method
CN111412391B (en) Pipe network leakage detection method and system
EP4163587A1 (en) Estimation device, estimation method, and estimation computer program for estimating a precipitate thickness
CN113268921B (en) Condenser cleaning coefficient estimation method and system, electronic device and readable storage medium
CN111159875A (en) Power station condenser shell side and tube side dynamic coupling mathematical model and modeling method
CN108345214B (en) Industrial process nonlinear detection method based on alternative data method
CN111104734A (en) Inversion prediction method for pumped storage power station unit load shedding test
CN115186456A (en) Heating system observability recovery method considering pipeline temperature quasi-dynamic state
CN109635404B (en) Steam pipeline pressure drop estimation method, device and system
CN108171408A (en) A kind of sewage water and water yield modeling method
TW201833518A (en) Method of real-time prognosis of flooding phenomenon in packed columns
CN112818495A (en) Novel dynamic correction method for pipeline pressure drop measurement and calculation algorithm parameters
CN103776652B (en) A kind of high-pressure heater method for testing performance and system
CN114819743B (en) Energy consumption diagnosis and analysis method for chemical enterprises
CN113901628A (en) Method for simulating hot oil pipeline
CN117824772B (en) Natural gas flow metering self-adaptive compensation method, system, terminal and medium

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