A Hydraulic-Thermal Coupling Prediction and Simulation Method and System for Central Heating Systems

Abstract

This invention discloses a hydraulic-thermal coupling prediction simulation method and system for centralized heating systems, relating to the field of urban heating pipe network technology. The method includes setting the heat source supply water temperature, the target load of the heat exchange station, and the outdoor air temperature as time-varying boundary conditions, generating a boundary dataset using an initial simulation condition set; employing a heat exchange station heat exchange state solution method, calculating the primary flow rate and primary return water temperature of the heat exchange station based on the boundary dataset to obtain the dynamic parameters of the heat exchange station; using a steady-state hydraulic analysis algorithm, solving the pipe network hydraulics based on the dynamic parameters of the heat exchange station to generate the pipe segment flow distribution; and based on the pipe segment flow distribution, the primary return water temperature of the heat exchange station, and the heat source supply water temperature, employing a dynamic thermodynamic calculation method to solve for pipe temperature changes, updating the temperature distribution of all network nodes, and iteratively advancing the simulation timeline to generate the pipe network flow distribution, temperature distribution, and hydraulic operating parameters. This invention improves the completeness of hydraulic-thermal coupling prediction and the practicality of simulation analysis.

Description

Technical Field

[0001] This invention relates to the field of urban heating network technology, and in particular to a hydraulic-thermal coupling prediction and simulation method and system for centralized heating systems. Background Technology

[0002] Currently, many cities in my country have built large-scale multi-heat-source ring-shaped heating networks and are developing smart operation platforms based on large-scale heating network simulation technology.

[0003] Current simulation methods mostly only consider the simulation calculation of the heating network itself, lacking consideration of the heat exchange capacity of heat exchange stations. In the application of simulation strategies, most methods can only manually set the boundary flow rates of heat exchange stations, making it difficult to match them with actual flow rate and return water temperature changes at heat exchange stations. Therefore, they can only be used for designing pipeline network operation schemes, lacking the ability to predict real-time dynamics of the pipeline network. Summary of the Invention

[0004] The predictive simulation described in this invention refers to a dynamic simulation process that, based on the current simulation state set and combined with time-varying boundary conditions such as the heat source water supply temperature, the target load of the heat exchange station, and the outdoor air temperature, continuously and recursively solves the heat exchange state of the heat exchange station, the hydraulic state of the pipeline network, and the thermal state of the pipeline network according to the simulation time interval, thereby obtaining the pipeline network flow distribution, temperature distribution, and hydraulic operating parameters at each subsequent simulation time.

[0005] In view of the aforementioned existing problems, the present invention is proposed.

[0006] Therefore, this invention provides a hydraulic-thermal coupling prediction simulation method for centralized heating systems to solve the problem of insufficient hydraulic-thermal coupling prediction capability under continuous operation at multiple times.

[0007] To solve the above-mentioned technical problems, the present invention provides the following technical solution:

[0008] In a first aspect, the present invention provides a hydraulic-thermal coupling prediction and simulation method for centralized heating systems, comprising: collecting pipeline parameter data, heat source operation data, heat exchange station operation data, and outdoor air temperature data; generating an initial simulation condition set by initializing pipe segment flow rate and node temperature and setting an initial simulation time; setting the heat source supply water temperature, heat exchange station target load, and outdoor air temperature as time-varying boundary conditions, and generating a boundary dataset using the initial simulation condition set; employing a heat exchange station heat exchange state solution method to calculate the primary side flow rate and primary side return water temperature of the heat exchange station based on the boundary dataset to obtain the dynamic parameters of the heat exchange station; using a steady-state hydraulic analysis algorithm to solve the pipeline hydraulic parameters based on the heat exchange station dynamic parameters to generate the pipe segment flow distribution; and using a dynamic thermodynamic calculation method based on the pipe segment flow distribution, the heat exchange station primary side return water temperature, and the heat source supply water temperature to solve the pipeline temperature change, update the temperature distribution of all network nodes, and iterate the simulation time cycle to generate the pipeline flow distribution, temperature distribution, and hydraulic condition parameters.

[0009] As a preferred embodiment of the hydraulic-thermal coupling prediction and simulation method for centralized heating systems described in this invention, the pipeline parameter data includes: pipeline connection relationship, number of nodes in each pipeline segment, pipeline impedance, pipeline length, pipeline diameter, pipeline height difference, pipeline thermal resistance, pump parameters, initial temperature, heat exchange station node number, heat source node number, heat source constant pressure point pressure, simulation duration, and simulation time interval.

[0010] The heat source operation data includes heat source supply water temperature, heat source return water temperature, heat source inlet pressure, heat source outlet pressure, heat source circulation flow rate, heat source heating power, and heat source regulation setpoint.

[0011] The heat exchange station operation data includes the heat exchange station target load, heat exchange station actual load, primary side supply water temperature, heat exchange station primary side return water temperature, primary side flow rate, secondary side supply water temperature, and secondary side return water temperature.

[0012] The meteorological temperature data includes current outdoor temperature, historical outdoor temperature, and predicted outdoor temperature.

[0013] As a preferred embodiment of the hydraulic-thermal coupling prediction and simulation method for centralized heating systems described in this invention, the specific steps for generating the initial simulation condition set are as follows:

[0014] Based on the pipeline network parameter data, a two-dimensional structure of the initial pipeline segment flow matrix is ​​established, where the first dimension represents the simulation time and the second dimension represents the pipeline segment number. The initial flow value of each pipeline segment is set to zero to generate the initial pipeline segment flow matrix.

[0015] Based on the pipeline network parameter data, a three-dimensional structure of the initial moment node temperature matrix is ​​established, where the first dimension represents the simulation moment, the second dimension represents the pipe segment number, and the third dimension represents the node number inside the pipe segment. The temperature of each node is set to the initial temperature to generate the initial moment node temperature matrix.

[0016] The initial simulation time is set to zero, and the initial pipe segment flow matrix and the initial time node temperature matrix are used as initialization inputs to generate the initial simulation condition set.

[0017] As a preferred embodiment of the hydraulic-thermal coupling prediction and simulation method for centralized heating systems described in this invention, the specific steps for generating the boundary dataset are as follows:

[0018] Extract pipe segment flow rate information and node temperature information at the current simulation moment from the initial simulation condition set to establish the current simulation state set;

[0019] Based on the heat source operation data, heat exchange station operation data and meteorological temperature data, time series matching and processing are performed to form the heat source water supply temperature time series, the heat exchange station target load time series and the outdoor temperature time series.

[0020] Based on the current simulation state set, extract the heat source water supply temperature data at the current simulation moment from the heat source water supply temperature time series.

[0021] Based on the current simulation state set, extract the target load data of the heat exchange station at the current simulation moment from the target load time series of the heat exchange station;

[0022] Based on the current simulation state set, extract the outdoor temperature data at the current simulation moment from the outdoor temperature time series;

[0023] The heat source water supply temperature data, heat exchange station target load data, and outdoor air temperature data are set as time-varying boundary conditions, and a boundary dataset is generated by combining them with the current simulation state set.

[0024] As a preferred embodiment of the hydraulic-thermal coupling prediction and simulation method for centralized heating systems described in this invention, the specific steps for obtaining the dynamic parameters of the heat exchange station are as follows:

[0025] The target load data of each heat exchange station at the current simulation moment is obtained from the boundary dataset, and the primary side water supply temperature data of each heat exchange station at the current simulation moment is obtained from the current simulation state set. At the same time, the secondary side water supply temperature set value and the target load data of the heat exchange station at the previous simulation moment are obtained from the heat exchange station operation data.

[0026] The heat exchange station state solution method is adopted. Based on the primary side water supply temperature data, the secondary side water supply temperature setpoint, the target load data of the heat exchange station at the current simulation time, and the target load data of the heat exchange station at the previous simulation time, the state equation of the heat exchange station is solved, and the primary side flow rate and the primary side return water temperature of each heat exchange station are calculated to obtain the dynamic parameters of the heat exchange station.

[0027] As a preferred embodiment of the hydraulic-thermal coupling prediction and simulation method for centralized heating systems described in this invention, the specific steps for generating the flow distribution of the pipe segment are as follows:

[0028] Based on the pipeline network parameter data, the connection relationship between each pipe segment and the topological nodes at both ends of the pipe segment is analyzed to form the pipeline network topology information, and the closed loop relationship formed by each pipe segment is identified to establish the basic correlation matrix and basic loop matrix of the pipeline network.

[0029] Extract the pipe segment connection relationship and flow direction relationship of each node from the basic correlation matrix of the pipeline network, and establish the flow conservation relationship of each node based on the heat source circulation flow and the primary side flow of the heat exchange station node, and construct the node flow balance equation.

[0030] Extract the pipe segment connection relationship and flow direction relationship of each closed loop from the basic loop matrix of the pipe network, and combine the pressure drop of each pipe segment and the pump head to form a pressure balance relationship, and construct the loop pressure balance equation;

[0031] Based on the dynamic parameters of the heat exchange station and the flow balance equation of the nodes, the flow of the remaining branches in the pipe section is initially set, and the branch flow is calculated to generate the hydraulic iteration initial value matrix at the current moment.

[0032] A steady-state hydraulic analysis algorithm is used, with the initial value matrix of the hydraulic iteration at the current moment as the initial value, to iteratively solve the nodal flow balance equation and the loop pressure balance equation to generate the flow distribution of the pipe segment.

[0033] As a preferred embodiment of the hydraulic-thermal coupling prediction and simulation method for centralized heating systems described in this invention, the specific steps for generating the pipe network flow distribution, temperature distribution, and hydraulic operating parameters are as follows:

[0034] The boundary dataset is updated based on the current simulation state set, and the heat exchange state of the heat exchange station is predicted based on the updated boundary dataset to obtain the heat exchange response result set of the station.

[0035] Steady-state hydraulic calculations and dynamic thermal calculations are performed on the heat exchange response result set of the site to generate the flow distribution, temperature distribution and hydraulic operating parameters of the pipeline network.

[0036] The simulation time is advanced according to the preset simulation time interval, and the simulation is iterated in cycles based on the flow distribution of pipe sections, the temperature distribution of nodes in the whole network, and the hydraulic operating parameters at each simulation time until the preset total simulation time is reached.

[0037] As a preferred embodiment of the hydraulic-thermal coupling prediction and simulation method for centralized heating systems described in this invention, the specific steps for advancing the current simulation time according to a preset simulation time interval are as follows:

[0038] The current simulation time is incremented according to the simulation time interval, and it is determined whether the incremented simulation time reaches the preset total simulation duration.

[0039] If the incremented simulation time does not reach the preset total simulation time, the updated pipeline flow distribution and temperature distribution will be used as the current simulation state set for the next round of simulation to continue the loop calculation.

[0040] The loop iteration ends when the incremented simulation time reaches the preset total simulation duration.

[0041] As a preferred embodiment of the hydraulic-thermal coupling prediction and simulation method for centralized heating systems described in this invention, the specific steps for determining the primary side water supply temperature data and the secondary side water supply temperature setpoint of the heat exchange station are as follows:

[0042] Extract the node temperature information at the current simulation moment from the node temperature information in the current simulation state set;

[0043] According to the node number of the heat exchange station, the water supply temperature of each heat exchange station node at the current simulation time is extracted from the node temperature information at the current simulation time to form the primary side water supply temperature data of the heat exchange station.

[0044] Extract the secondary water supply temperature setting information of the heat exchange station at the current simulation moment from the heat exchange station operation data;

[0045] Based on the heat exchange station node number, extract the secondary side water supply temperature setting value of each heat exchange station at the current simulation time from the heat exchange station secondary side water supply temperature setting information at the current simulation time.

[0046] Secondly, this invention provides a hydraulic-thermal coupling prediction and simulation system for a centralized heating system, comprising: an initial setting module for collecting pipeline parameter data, heat source operation data, heat exchange station operation data, and meteorological temperature data; initializing pipe segment flow rate and node temperature based on the pipeline parameter data and setting the initial simulation time; and generating an initial simulation condition set; a boundary setting module for setting the heat source supply water temperature, heat exchange station target load, and outdoor air temperature as time-varying boundary conditions; and generating a boundary dataset using the initial simulation condition set; and a heat exchange prediction module for predicting the heat exchange status of the heat exchange station. The system employs a steady-state solution method, which calculates the primary flow rate and primary return water temperature of the heat exchange station based on the boundary dataset to obtain the dynamic parameters of the heat exchange station. A hydraulic solution module utilizes a steady-state hydraulic analysis algorithm to solve the hydraulic problems of the pipe network based on the dynamic parameters of the heat exchange station, generating the flow distribution of pipe segments. A thermal iteration module uses a dynamic thermal calculation method based on the flow distribution of pipe segments, the primary return water temperature of the heat exchange station, and the supply water temperature of the heat source to solve the pipe temperature changes, update the temperature distribution of all nodes in the network, and iterates continuously during the simulation to generate the flow distribution, temperature distribution, and hydraulic operating parameters of the pipe network.

[0047] The beneficial effects of this invention are as follows: by using dynamic thermodynamic calculation methods to solve the pipeline temperature change based on the flow distribution of the pipe section, the primary return water temperature of the heat exchange station, and the supply water temperature of the heat source, the temperature distribution of all network nodes is updated and the simulation time is cyclically iterated, realizing the continuous simulation of the operation process of the centralized heating system. It can be used to output the flow distribution, temperature distribution and hydraulic operating parameters of the pipeline network, thereby improving the completeness of hydraulic-thermal coupling prediction and the practicality of simulation analysis. Attached Figure Description

[0048] To more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings used in the following description of the embodiments will be briefly introduced. Obviously, the drawings described below are only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort.

[0049] Figure 1 This is a flowchart of a hydraulic-thermal coupling prediction and simulation method for centralized heating systems.

[0050] Figure 2 This is a schematic diagram of a hydraulic-thermal coupling prediction simulation system for a centralized heating system.

[0051] Figure 3 This is a flowchart for calculating the dynamic parameters of a heat exchange station.

[0052] Figure 4 This is a flowchart for dynamic thermal analysis of the pipeline network. Detailed Implementation

[0053] To make the above-mentioned objects, features and advantages of the present invention more apparent and understandable, the specific embodiments of the present invention will be described in detail below with reference to the accompanying drawings.

[0054] Many specific details are set forth in the following description in order to provide a full understanding of the invention. However, the invention may also be practiced in other ways different from those described herein, and those skilled in the art can make similar extensions without departing from the spirit of the invention. Therefore, the invention is not limited to the specific embodiments disclosed below.

[0055] Secondly, the term "one embodiment" or "embodiment" as used herein refers to a specific feature, structure, or characteristic that may be included in at least one implementation of the present invention. The phrase "in one embodiment" appearing in different places in this specification does not necessarily refer to the same embodiment, nor is it a single or selective embodiment that is mutually exclusive with other embodiments.

[0056] Reference Figures 1-4 As an embodiment of the present invention, this embodiment provides a hydraulic-thermal coupling prediction and simulation method for a centralized heating system, comprising the following steps:

[0057] S1. Collect pipeline parameter data, heat source operation data, heat exchange station operation data and outdoor temperature data. Generate the initial simulation condition set by initializing the pipe section flow rate and node temperature and setting the initial simulation time.

[0058] S1.1 Pipeline network parameter data includes pipe segment connection relationships, number of nodes in each pipe segment, pipe segment impedance, pipe segment length, pipe segment diameter, pipe segment height difference, pipe segment thermal resistance, water pump parameters, initial temperature, heat exchange station node number, heat source node number, heat source constant pressure point pressure, simulation duration and simulation time interval.

[0059] Specifically, the pipe segment connection relationships, heat exchange station node numbers, and heat source node numbers are obtained through analysis of the centralized heating network topology diagram, design drawings, or as-built drawings; the number of nodes in each pipe segment is determined based on the pipe segment discrete calculation requirements, combined with the pipe segment length and the set discrete accuracy; the pipe segment impedance is calculated based on the pipe segment length, pipe diameter, pipe wall roughness, and local resistance conditions; the pipe segment length, pipe segment diameter, and pipe segment height difference are obtained from design parameters, survey results, or as-built data; the pipe segment thermal resistance is calculated based on the pipe material, insulation layer structure, burial method, and heat exchange conditions of the surrounding medium; the pump parameters are obtained from pump equipment data, performance curves, or operating parameters; the initial temperature is set based on the monitoring temperature at the start of the simulation or the preset operating conditions; the pressure at the heat source constant pressure point is obtained from the operating set value or real-time monitoring value on the heat source side; the simulation duration and simulation time interval are set according to the simulation task scope and time resolution requirements.

[0060] S1.2 The heat source operation data includes the heat source supply water temperature, heat source return water temperature, heat source inlet pressure, heat source outlet pressure, heat source circulation flow rate, heat source heating power, and heat source regulation setting value.

[0061] Specifically, the supply and return water temperatures of the heat source are obtained in real time by temperature sensors installed on the supply and return water pipelines. The inlet and outlet pressures of the heat source are measured by pressure sensors or pressure transmitters installed at the inlet and outlet of the heat source. The circulating flow rate of the heat source is measured by a flow meter installed on the circulating pipeline. The heating power of the heat source is calculated based on the circulating flow rate and the temperature difference between the supply and return water temperatures. The heat source regulation setpoint is determined based on the heating scheduling strategy, controller settings, or operation control commands.

[0062] S1.3 The heat exchange station operation data includes the target load of the heat exchange station, the actual load of the heat exchange station, the primary side supply water temperature, the primary side return water temperature of the heat exchange station, the primary side flow rate, the secondary side supply water temperature, and the secondary side return water temperature.

[0063] Specifically, the target load of the heat exchange station is determined based on the preset heating scheduling plan, user-side heat demand forecast results, building area parameters, historical load patterns, and outdoor temperature compensation rules. The actual load of the heat exchange station is calculated based on the primary side flow rate and the temperature difference between the primary side supply water temperature and the primary side return water temperature of the heat exchange station. The primary side supply water temperature and the primary side return water temperature of the heat exchange station are obtained in real time by temperature sensors installed on the primary side supply water pipeline and return water pipeline. The primary side flow rate is obtained by flow meter installed on the primary side pipeline. The secondary side supply water temperature and the secondary side return water temperature are obtained in real time by temperature sensors installed on the secondary side supply water pipeline and return water pipeline.

[0064] S1.4 Meteorological temperature data includes current outdoor temperature, historical outdoor temperature and forecast outdoor temperature.

[0065] Specifically, the current outdoor temperature is obtained in real time by outdoor temperature sensors installed at the heating area site; the historical outdoor temperature is obtained by retrieving historical meteorological records stored in the meteorological monitoring platform, regional meteorological database, or historical sampling data of the heating season; and the predicted outdoor temperature is obtained directly from the forecast data released by the meteorological forecast platform, or calculated by using time series forecasting methods based on historical outdoor temperature, current outdoor temperature, and meteorological change trends.

[0066] S1.5. Based on the pipeline network parameter data, establish a two-dimensional structure of the initial pipeline flow matrix, where the first dimension represents the simulation time and the second dimension represents the pipeline segment number. Set the initial flow value of each pipeline segment to zero to generate the initial pipeline flow matrix.

[0067] Specifically, based on the pipe segment connection relationships and the number of nodes in each pipe segment in the pipeline network parameter data, all pipe segment numbers in the pipeline network parameter data are counted and the column index in the two-dimensional structure is determined; based on the simulation duration and simulation time interval in the pipeline network parameter data, the simulation time sequence is divided and the row index in the two-dimensional structure is determined; a two-dimensional storage structure is constructed according to the correspondence between simulation time and pipe segment number, so that each row corresponds to one simulation time and each column corresponds to one pipe segment; after completing the construction of the two-dimensional storage structure, the position of each pipe segment at the initial simulation time is assigned to zero one by one to obtain the initial pipe segment flow value, and the initial pipe segment flow value is written into the two-dimensional structure to generate the initial pipe segment flow matrix.

[0068] S1.6. Based on the pipeline network parameter data, establish a three-dimensional structure of the node temperature matrix at the initial moment, where the first dimension represents the simulation moment, the second dimension represents the pipe segment number, and the third dimension represents the node number inside the pipe segment. Set the temperature of each node to the initial temperature to generate the node temperature matrix at the initial moment.

[0069] Specifically, based on the simulation duration and simulation time interval in the pipeline network parameter data, the simulation time sequence of the first dimension of the initial moment node temperature matrix is ​​determined; based on the pipe segment connection relationship in the pipeline network parameter data, the pipe segment number sequence of the second dimension of the initial moment node temperature matrix is ​​determined; based on the number of nodes in each pipe segment in the pipeline network parameter data, the range of internal node numbers for each pipe segment in the third dimension of the initial moment node temperature matrix is ​​determined; after determining the correspondence between the simulation time sequence, the pipe segment number sequence, and the range of internal node numbers, a three-dimensional structure is constructed according to the index relationship of the first, second, and third dimensions; after completing the construction of the three-dimensional structure, the positions of the internal nodes of each pipe segment at the initial simulation moment are assigned initial temperatures one by one, and the assigned node temperatures are written into the three-dimensional structure according to the simulation time, pipe segment number, and internal node number to generate the initial moment node temperature matrix.

[0070] S1.7 Set the initial simulation time to zero, and use the initial pipe segment flow matrix and the initial time node temperature matrix as initialization inputs to generate the initial simulation condition set.

[0071] Specifically, the initial simulation time is set to zero, and a correspondence is established between the initial simulation time and the first-dimensional index of the initial pipe segment flow matrix and the first-dimensional index of the initial time node temperature matrix. According to the correspondence of the initial simulation time, the flow position of each pipe segment in the initial pipe segment flow matrix at the initial simulation time is calibrated, and the internal node temperature position of each pipe segment in the initial time node temperature matrix at the initial simulation time is calibrated. After the initial simulation time calibration is completed, the initial pipe segment flow matrix and the initial time node temperature matrix are associated and integrated according to the initial simulation time, so that the pipe segment flow information in the initial pipe segment flow matrix and the node temperature information in the initial time node temperature matrix are stored correspondingly under the same initial simulation time, generating the initial simulation condition set.

[0072] S2. Set the heat source water supply temperature, the target load of the heat exchange station and the outdoor air temperature as time-varying boundary conditions, and generate the boundary dataset using the initial simulation condition set.

[0073] S2.1 Extract the pipe flow rate information and node temperature information at the current simulation moment from the initial simulation condition set to establish the current simulation state set.

[0074] Specifically, based on the time index position of the current simulation time in the initial simulation condition set, the initial pipe segment flow matrix and the initial node temperature matrix for the current simulation time are located. According to the time index of the current simulation time, the flow value of each pipe segment at the current simulation time is read one by one from the initial pipe segment flow matrix to form the pipe segment flow information for the current simulation time. According to the time index of the current simulation time, the temperature value of each node within each pipe segment at the current simulation time is read one by one from the initial node temperature matrix to form the node temperature information for the current simulation time. The pipe segment flow information and the node temperature information at the current simulation time are then associated and organized according to the corresponding relationship at the same time to establish the current simulation state set.

[0075] S2.2. Based on the heat source operation data, heat exchange station operation data and meteorological temperature data, time series matching and processing are performed to form the heat source water supply temperature time series, the heat exchange station target load time series and the outdoor temperature time series.

[0076] Specifically, the heat source water supply temperature corresponding to each recorded moment is extracted from the heat source operation data, the heat exchange station target load corresponding to each recorded moment is extracted from the heat exchange station operation data, and the outdoor temperature corresponding to each recorded moment is extracted from the meteorological temperature data. The recorded moments corresponding to the heat source water supply temperature, heat exchange station target load, and outdoor temperature are uniformly sorted and matched moment by moment. The matched heat source water supply temperature, heat exchange station target load, and outdoor temperature are then arranged in the same time order to fill in missing moments and form a time series of heat source water supply temperature, heat exchange station target load, and outdoor temperature.

[0077] S2.3. Based on the current simulation state set, extract the heat source water supply temperature data at the current simulation moment from the heat source water supply temperature time series.

[0078] Specifically, based on the current simulation time information in the current simulation state set, the corresponding time position of the current simulation time in the heat source water supply temperature time series is determined; according to the time correspondence between the current simulation time and the heat source water supply temperature time series, the heat source water supply temperature time series is matched to locate the temperature record position of the current simulation time; the heat source water supply temperature value is read from the temperature record position of the current simulation time, and the read heat source water supply temperature value is correlated with the current simulation time in the current simulation state set to obtain the heat source water supply temperature data of the current simulation time.

[0079] S2.4. Based on the current simulation state set, extract the target load of the heat exchange station at the current simulation moment from the target load time series of the heat exchange station.

[0080] Specifically, based on the current simulation time information in the current simulation state set, the corresponding time position of the current simulation time in the target load time series of the heat exchange station is determined; according to the time correspondence between the current simulation time and the target load time series of the heat exchange station, the target load time series of the heat exchange station is matched to locate the load record position of the current simulation time; the target load value of each heat exchange station is read one by one from the load record position of the current simulation time, and the read target load value of each heat exchange station is correlated with the current simulation time in the current simulation state set to obtain the target load data of the heat exchange station at the current simulation time.

[0081] S2.5. Based on the current simulation state set, extract the outdoor temperature data at the current simulation moment from the outdoor temperature time series.

[0082] Specifically, based on the current simulation time information in the current simulation state set, the corresponding time position of the current simulation time in the outdoor temperature time series is determined; according to the time correspondence between the current simulation time and the outdoor temperature time series, the outdoor temperature time series is time-matched to locate the temperature recording position of the current simulation time; the outdoor temperature value is read from the temperature recording position of the current simulation time, and the read outdoor temperature value is correlated with the current simulation time in the current simulation state set to obtain the outdoor temperature data of the current simulation time.

[0083] S2.6 Set the heat source water supply temperature data, heat exchange station target load data and outdoor air temperature data as time-varying boundary conditions, and generate a boundary dataset by combining them with the current simulation state set.

[0084] Specifically, the heat source water supply temperature data, heat exchange station target load data, and outdoor temperature data are aligned according to the current simulation time, and a time-to-time correspondence is established among these three data points. Based on this time-to-time correspondence, the heat source water supply temperature data, heat exchange station target load data, and outdoor temperature data are uniformly calibrated as time-varying boundary conditions for the current simulation time. After calibrating the time-varying boundary conditions, they are associated and integrated with the pipe flow information and node temperature information in the current simulation state set according to the current simulation time. This ensures that the heat source water supply temperature data, heat exchange station target load data, outdoor temperature data, pipe flow information, and node temperature information are stored correspondingly at the same current simulation time, generating a boundary dataset.

[0085] S3. Using the heat exchange state solution method of the heat exchange station, calculate the primary flow rate and primary return water temperature of the heat exchange station based on the boundary dataset to obtain the dynamic parameters of the heat exchange station.

[0086] S3.1 Obtain the target load data of each heat exchange station at the current simulation time from the boundary dataset, and obtain the primary side water supply temperature data of each heat exchange station at the current simulation time based on the current simulation state set. At the same time, obtain the secondary side water supply temperature set value of the heat exchange station and the target load data of the heat exchange station at the previous simulation time based on the heat exchange station operation data.

[0087] Specifically, based on the current simulation time index in the boundary dataset, the target load record position of each heat exchanger station at the current simulation time is located one by one, and the target load value of each heat exchanger station at the current simulation time is read to form the target load data of each heat exchanger station at the current simulation time; based on the node temperature information and heat exchanger station node number in the current simulation state set, the inlet node position of each heat exchanger station at the current simulation time is matched one by one, and the supply water temperature value of the inlet node of each heat exchanger station is read from the node temperature information at the current simulation time to form the target load data of each heat exchanger station at the current simulation time. The primary side water supply temperature data of the heat exchange station; after reading the target load data and primary side water supply temperature data of each heat exchange station at the current simulation time, the secondary side water supply temperature setpoint of each heat exchange station is retrieved one by one according to the correspondence of the heat exchange station node numbers. Based on the time correspondence between the current simulation time and the previous simulation time, the target load value of each heat exchange station at the previous simulation time is read one by one from the boundary dataset. At the same time, the secondary side water supply temperature setpoint and the target load data of the heat exchange station at the previous simulation time are obtained based on the heat exchange station operation data.

[0088] S3.2 Extract the node temperature information at the current simulation moment from the node temperature information in the current simulation state set.

[0089] Specifically, based on the current simulation time information in the current simulation state set, the corresponding time position of the current simulation time in the node temperature information is determined; according to the time correspondence between the current simulation time and the node temperature information, the node temperature information is time-matched to locate the temperature recording range of the current simulation time; after completing the location of the temperature recording range of the current simulation time, the internal node temperature values ​​of each pipe segment at the current simulation time are read one by one, and the internal node temperature values ​​of each pipe segment are organized according to the pipe segment number and the internal node number of the pipe segment; after completing the organization of the internal node temperature values ​​of each pipe segment, the internal node temperature values ​​of each pipe segment are correlated with the current simulation time to obtain the node temperature information of the current simulation time.

[0090] S3.3. According to the node number of the heat exchange station, extract the water supply temperature of each heat exchange station node at the current simulation time from the node temperature information at the current simulation time to form the primary side water supply temperature data of the heat exchange station.

[0091] Specifically, according to the correspondence between the node number of the heat exchange station and the node temperature information at the current simulation moment, the temperature record position of each heat exchange station node at the current simulation moment is located one by one; based on the temperature record position of each heat exchange station node at the current simulation moment, the water supply temperature value of each heat exchange station node at the current simulation moment is read; the water supply temperature values ​​of each heat exchange station node at the current simulation moment are organized according to the heat exchange station node number, and the organized water supply temperature values ​​of each heat exchange station node are written to the data location to form the primary side water supply temperature data of the heat exchange station.

[0092] S3.4 Extract the secondary water supply temperature setting information of the heat exchange station at the current simulation moment from the heat exchange station operation data.

[0093] Specifically, based on the time correspondence in the heat exchange station operation data at the current simulation time, the location of the secondary side water supply temperature setting record of each heat exchange station at the current simulation time is determined, and the secondary side water supply temperature setting value of each heat exchange station at the current simulation time is read one by one to form the secondary side water supply temperature setting information of the heat exchange station at the current simulation time.

[0094] S3.5. According to the heat exchange station node number, extract the secondary side water supply temperature setting value of each heat exchange station at the current simulation time from the heat exchange station secondary side water supply temperature setting information at the current simulation time.

[0095] Specifically, according to the correspondence between the heat exchange station node number and the secondary side water supply temperature setting information of the heat exchange station at the current simulation time, the setting value record position of each heat exchange station at the current simulation time is located one by one; based on the setting value record position of each heat exchange station at the current simulation time, the secondary side water supply temperature setting value of each heat exchange station at the current simulation time is read; the secondary side water supply temperature setting values ​​of each heat exchange station at the current simulation time are organized according to the heat exchange station node number to form the secondary side water supply temperature setting value of the heat exchange station.

[0096] S3.6. Using the heat exchange station heat exchange state solution method, based on the primary side water supply temperature data, the secondary side water supply temperature setpoint, the target load data of the heat exchange station at the current simulation time, and the target load data of the heat exchange station at the previous simulation time, the heat exchange station state equation is solved, the primary side flow rate and the primary side return water temperature of each heat exchange station are calculated, and the dynamic parameters of the heat exchange station are obtained.

[0097] Specifically, a heat exchanger state solution method is adopted. The primary side water supply temperature data, secondary side water supply temperature setpoints, target load data at the current simulation moment, and target load data from the previous simulation moment are aligned and matched according to the corresponding relationships between heat exchangers, forming the input terms of the heat exchanger state equation for each heat exchanger. The heat exchanger is divided into multiple control volumes along the flow direction, and the parameters of each control volume are organized to form a heat exchange calculation parameter set. The input terms of the heat exchanger state equation for each heat exchanger are substituted one by one into the heat exchanger state equation, and the heat exchange process at each heat exchanger at the current simulation moment is iteratively solved. The expression of the heat exchanger state equation is:

[0098] ;

[0099] in, Indicates the high temperature side Individual control body temperature Indicates the low temperature side Individual control body temperature This indicates the primary mass flow rate on the high-temperature side. This indicates the secondary side mass flow rate on the low-temperature side. Indicates the width of the flow channel between plates. Indicates the spacing between boards. Indicates the length of the flow channel. Indicates the overall heat transfer coefficient. This represents the volumetric heat capacity of water. Indicates heat exchange efficiency. Indicates the high-temperature side. Indicates the low-temperature side. It represents water;

[0100] Based on the heat exchange calculation results of each heat exchange station, the primary flow rate and primary return water temperature of each heat exchange station are calculated one by one. The primary flow rate and primary return water temperature of each heat exchange station are then correlated and organized according to the corresponding relationship of the heat exchange stations to obtain the dynamic parameters of the heat exchange stations.

[0101] It should be noted that the heat exchange station is not treated as a static boundary point. Instead, the target load of the heat exchange station at the current simulation moment, the load at the previous moment, the set values ​​of the primary side supply water temperature and the secondary side supply water temperature of the heat exchange station are all substituted into the state equation of the heat exchange station. After solving within a time period, the primary side flow rate and the primary side return water temperature of the heat exchange station are directly output and then fed into the subsequent steady-state hydraulic calculation and dynamic thermodynamic analysis, forming a dynamic coupling interface between the heat exchange process in the station and the hydraulic and thermodynamic calculation of the pipeline network.

[0102] S4. Using a steady-state hydraulic analysis algorithm, solve the hydraulics of the pipe network based on the dynamic parameters of the heat exchange station to generate the flow distribution of the pipe section.

[0103] S4.1 Based on the pipeline network parameter data, analyze the connection relationship between each pipe segment and the topological nodes at both ends of the pipe segment to form the pipeline network topology information, identify the closed loop relationship formed by each pipe segment, and establish the basic correlation matrix and basic loop matrix of the pipeline network.

[0104] Specifically, based on the pipe segment connection relationships in the pipeline network parameter data, the starting and ending topological nodes of each pipe segment are determined one by one. The correspondence between each pipe segment and the topological nodes at both ends of the segment is then organized according to the pipe segment number to form the pipeline network topology information. Based on the pipeline network topology information, each pipe segment is traversed according to the connectivity between the topological nodes at both ends of the segment, and the pipe segment numbers passed through in sequence are recorded. When the ending topological node in the path coincides with the starting topological node, all pipe segments in the corresponding path are identified as a closed loop. The combinations of pipe segments in all closed loops are organized to form the closed loop relationships formed by each pipe segment. After completing the organization of the pipeline network topology information and closed loop relationships, the correspondence between the pipe segment number and the topological nodes at both ends of the segment is written line by line into the matrix position, and the closed loop relationships formed by each pipe segment are written loop by loop into the matrix position, establishing the basic pipeline network association matrix and the basic pipeline network loop matrix.

[0105] S4.2 Extract the pipe segment connection relationship and flow direction relationship of each node from the basic correlation matrix of the pipeline network, and establish the flow conservation relationship of each node based on the heat source circulation flow and the primary side flow of the heat exchange station node, and construct the node flow balance equation.

[0106] Specifically, the non-zero matrix elements of each node are read row by row from the basic correlation matrix of the pipeline network. The pipe segment numbers and connection directions connected to each node are identified based on the column positions of the non-zero matrix elements, forming the pipe segment information connected to each node. Based on this information, the flow rates flowing into and out of each node are statistically analyzed, and these flow rates are correlated with the corresponding node boundary flows according to the same node number. After organizing the node correspondences, a balance expression relationship is established between the sum of the flow rates flowing into and out of each node and the node boundary flows for each node, forming a flow conservation relationship. Based on the pipe segment connection relationships and flow directions corresponding to each node, and combined with the heat source circulation flow rate and the primary side flow rate of the heat exchange station, flow conservation equations are established for each node, and these equations are organized according to the node number to form the node flow balance equations.

[0107] S4.3 Extract the pipe segment connection relationship and flow direction relationship of each closed loop from the basic loop matrix of the pipe network, and combine the pressure drop of each pipe segment and the pump head to form a pressure balance relationship, and construct the loop pressure balance equation.

[0108] Specifically, the non-zero matrix elements of each closed loop are read row by row from the basic loop matrix of the pipeline network. Based on the position of the non-zero matrix element in each column, the pipe segment number and the direction of each pipe segment within each closed loop are identified, forming the pipe segment information for each closed loop. Based on the pipe segment information of each closed loop, the pressure drop of each pipe segment within each closed loop is retrieved one by one, and the pressure drop of each pipe segment is aligned according to its direction of flow within each closed loop. After aligning the pressure drop of each pipe segment, the pump head of each closed loop is retrieved, and the pressure drop and pump head of each pipe segment are correlated according to the same closed loop. After aligning the relationships of each closed loop, a balance expression relationship between the pressure drop of each pipe segment and the pump head is established for each closed loop, forming a pressure balance relationship. Based on the pipe segment connection relationships and flow direction relationships corresponding to each closed loop, and combined with the pressure drop and pump head of each pipe segment, pressure balance equations for each closed loop are established, and these equations are arranged according to the closed loop number to form the loop pressure balance equations.

[0109] S4.4. Based on the dynamic parameters of the heat exchange station and the node flow balance equation, the initial setting of the residual branch flow of the pipe section is performed, and the branch flow is calculated to generate the hydraulic iteration initial value matrix at the current moment.

[0110] Specifically, based on the primary flow rate of each heat exchange station in the dynamic parameters of the heat exchange station, and according to the correspondence between the heat exchange station node number and the pipeline topology information, the primary flow rate of each heat exchange station is allocated to the pipe segments connected to each heat exchange station node. The primary flow rate of each heat exchange station is then written into the node flow balance equation of the corresponding node. Combining the inflow and outflow directions of the pipe segments connected to each node, flow conservation constraints are established for each node. After distinguishing between pipe segment branches and branches based on the basic correlation matrix and the basic loop matrix of the pipeline network, the flow rates of the pipe segment branches are assigned one by one according to the closed loop correspondence, forming the initial values ​​of the pipe segment branch flow rates. After assigning the initial values ​​of the residual branch flow rates in the pipe segment, the back-substitution method of the node flow balance equation is adopted. The initial values ​​of the residual branch flow rates in the pipe segment are substituted into the node flow balance equation as known quantities. Starting from the node containing only one unknown branch flow rate, the corresponding branch flow rates are solved node by node according to the flow conservation relationship of each node until all branch flows corresponding to the residual branch flow rates in the pipe segment are obtained. After the branch flow rates are obtained, the initial values ​​of the residual branch flow rates in the pipe segment and the branch flow rates are written into the corresponding positions of the matrix according to the pipe segment number order, and stored uniformly according to the time index of the initial simulation time to generate the hydraulic iteration initial value matrix at the current time.

[0111] It should also be noted that the residual branch flow rate refers to the flow rate of the pipe segment outside the spanning tree used to form the closed loop after the pipe network forms a closed loop. Since the residual branch flow rate directly determines the flow distribution relationship in each closed loop, in steady-state hydraulic calculations, the residual branch flow rate is usually initially set first, and then the flow rate of other branches is deduced by combining the node flow balance equation, so as to obtain the overall network flow distribution that satisfies the node flow conservation relationship and the loop pressure balance relationship.

[0112] S4.5. Using a steady-state hydraulic analysis algorithm, with the initial value matrix of the hydraulic iteration at the current moment as the initial value, the flow balance equation of the nodes and the pressure balance equation of the loop are solved iteratively to generate the flow distribution of the pipe section.

[0113] Specifically, a steady-state hydraulic analysis algorithm is adopted. The initial flow rate of each pipe segment is read from the hydraulic iteration initial value matrix at the current moment as the initial value of the iteration. The nodal flow balance equations corresponding to each node are combined with the loop pressure balance equations corresponding to each closed loop. The nonlinear relationship between the pressure drop and the flow rate of each pipe segment is written into the loop pressure balance equation, forming a steady-state hydraulic nonlinear equation set. The expression of the steady-state hydraulic nonlinear equation set is as follows:

[0114] ;

[0115] in, Represents the basic correlation matrix of the pipeline network. This represents the column vector of flow rates in the pipe section. Represents the column vector of node boundary flows. Represents the basic loop matrix of the pipeline network. This represents the diagonal matrix of pipe segment impedance. The column vector represents the pump head; the first equation describes the flow conservation relationship at the node, and the second equation describes the energy balance relationship in the closed loop of the pipe network;

[0116] The flow balance results of each node and the pressure balance results of each closed loop are obtained based on the steady-state hydraulic nonlinear equations. Then, based on the flow balance results of each node and the pressure balance results of each closed loop, the flow rate of each pipe segment in the hydraulic iteration initial value matrix at the current moment is gradually corrected, and the corrected flow rate of each pipe segment is rewritten into the hydraulic iteration initial value matrix at the current moment. The rewritten flow rate of each pipe segment is substituted back into the node flow balance equation and the loop pressure balance equation to repeatedly check the flow balance results of each node and the pressure balance results of each closed loop, and the flow rate of each pipe segment is iteratively updated. When the flow balance results of each node and the pressure balance results of each closed loop meet the preset balance conditions, the flow rate values ​​of each pipe segment after iterative calculation are extracted to generate the pipe segment flow distribution.

[0117] It should also be noted that the preset balance conditions are usually determined based on the convergence accuracy requirements of the steady-state hydraulic analysis algorithm, the allowable error range of the pipeline network simulation, and the accuracy requirements of the engineering calculations for flow balance and pressure balance. For example, the preset balance conditions can be set to ensure that the absolute value of the flow residual at each node is not greater than 0.001 and the absolute value of the pressure residual at each closed loop is not greater than 0.01. When the flow balance results at each node and the pressure balance results at each closed loop are satisfied at the same time, it is determined that the preset balance conditions have been met.

[0118] S5. Based on the flow distribution of the pipe section, the primary return water temperature of the heat exchange station and the supply water temperature of the heat source, the dynamic thermodynamic calculation method is used to solve the temperature change of the pipeline, update the temperature distribution of all nodes in the network, and promote the simulation to iterate in cycles to generate the flow distribution, temperature distribution and hydraulic parameters of the pipeline network.

[0119] S5.1 Update the boundary dataset according to the current simulation state set, and predict the heat exchange state of the heat exchange station based on the updated boundary dataset to obtain the heat exchange response result set of the station.

[0120] Specifically, based on the current simulation time information in the current simulation state set, the corresponding records of heat source water supply temperature, heat exchange station target load, and outdoor air temperature in the boundary dataset are matched in time. Combined with the node temperature information and pipe flow information in the current simulation state set, the boundary values ​​in the boundary dataset are updated and organized to form an updated boundary dataset. After completing the organization of the updated boundary dataset, the heat exchange calculation input items of each heat exchange station at the current simulation time are extracted one by one according to the heat exchange station node number. The heat exchange calculation input items of each heat exchange station are then substituted one by one into the heat exchange state prediction process of the heat exchange station based on the heat exchange balance relationship and heat transfer relationship. An iterative solution method is used to update the primary side flow rate and the primary side return water temperature of the heat exchange station, obtaining the primary side flow rate and primary side return water temperature of each heat exchange station at the current simulation time. Finally, the primary side flow rate and primary side return water temperature of each heat exchange station at the current simulation time are correlated and unified according to the heat exchange station node number to obtain the station heat exchange response result set.

[0121] S5.2 Perform steady-state hydraulic calculations and dynamic thermal calculations on the heat exchange response result set of the site to generate the flow distribution, temperature distribution and hydraulic operating parameters of the pipeline network.

[0122] Specifically, the primary flow rate and primary return water temperature of each heat exchange station in the heat exchange response result set are matched according to the heat exchange station node number. Combined with the pipe segment connection relationships and the number of nodes in each pipe segment from the pipeline network parameter data, the primary flow rate of each heat exchange station is allocated to the corresponding pipe segment and node location. Based on the allocated primary flow rate of each heat exchange station, the pipe segment connection relationships, and the number of nodes in each pipe segment, node flow balance equations and loop pressure balance equations are established. Then, considering the nonlinear relationship between the pressure drop and flow rate of each pipe segment, a steady-state hydraulic analysis algorithm is used for iterative solution. When the flow residual at each node and the pressure residual in each closed loop meet the preset balance conditions, the flow rate value of each pipe segment, the pressure value of each node, and the pressure drop value of each pipe segment are obtained, forming the pipeline network flow distribution and hydraulic... Operating parameters: The pipeline flow distribution, the primary return water temperature of each heat exchange station, and the heat source supply water temperature are aligned to the current simulation time. Dynamic thermodynamic calculations are performed on each node within each pipe segment along the pipe connection relationships. Each pipe segment is divided into multiple control volumes. The discrete length and cross-sectional area of ​​the pipe segment are determined based on the pipeline network parameter data. The simulation time interval and current time are determined based on the simulation time settings. The flow velocity and temperature of each control volume at the current time are determined based on the flow distribution of the pipe segment and the temperature distribution of all network nodes. The density and specific heat capacity of the medium are determined based on preset medium physical property parameters. The total heat transfer thermal resistance is determined based on the pipe structure parameters and insulation parameters. The ambient temperature is determined based on meteorological temperature data, forming dynamic calculation parameters for pipeline temperature. The dynamic thermodynamic calculation formula is as follows:

[0123] ;

[0124] in, Indicates the current time. Individual control body temperature Indicates the previous moment. Individual control body temperature This indicates the temperature of the upstream adjacent control volume at the current moment. Indicates the simulation time interval. Let A represent the discrete length of the pipe segment and A represent the cross-sectional area of ​​the pipe. Indicates the density of the medium. Indicates the specific heat capacity of the medium. Indicates the total heat transfer thermal resistance. Indicates ambient temperature. Indicates the current flow rate. Indicates the current moment;

[0125] The temperature values ​​of the internal nodes of each pipe section are obtained based on the dynamic thermodynamic calculation formula, forming a temperature distribution. The temperature distribution is then generated by associating and organizing the data according to the current simulation time, resulting in the pipe network flow distribution, temperature distribution, and hydraulic parameters.

[0126] S5.3. Advance the current simulation time according to the preset simulation time interval, and iterate cyclically according to the flow distribution of pipe sections, the temperature distribution of nodes in the whole network, and the hydraulic operating parameters at each simulation time until the preset total simulation time is reached.

[0127] S5.3.1 Increment the current simulation time according to the simulation time interval, and determine whether the incremented simulation time reaches the preset total simulation duration.

[0128] Specifically, based on the time correspondence between the current simulation time and the simulation time interval, the current simulation time is accumulated to obtain the incremented simulation time. After obtaining the incremented simulation time, the incremented simulation time is compared with the preset total simulation duration to determine the time position of the incremented simulation time within the preset total simulation duration. Based on the comparison result between the incremented simulation time and the preset total simulation duration, it is determined whether the incremented simulation time has reached the preset total simulation duration.

[0129] It should also be noted that the preset total simulation duration is usually based on the operating time range to be covered by the simulation task, so as to ensure that the hydraulic and thermal changes of the pipeline network within the target operating cycle can be fully captured and can reflect typical operating conditions such as all-day, peak and valley loads or seasonal changes. The value range is determined by the simulation target and time resolution, and can generally be set from single hourly operation to multi-day continuous operation. It needs to be coordinated with the simulation time interval so that the incremented simulation time can fully cover the entire preset operating cycle, while ensuring that the calculation data of each simulation time is continuous and can be used for the next hydraulic and thermal iteration.

[0130] S5.3.2 If the incremented simulation time does not reach the preset total simulation time, the updated pipeline flow distribution and temperature distribution will be used as the current simulation state set for the next round of simulation to continue the loop calculation.

[0131] Specifically, when the incremented simulation time has not reached the preset total simulation duration, the incremented simulation time is first associated with the updated pipeline flow distribution, temperature distribution, and hydraulic parameters at the same time. The flow values ​​of each pipe segment in the updated pipeline flow distribution, the temperature values ​​of each node in the temperature distribution, and the pressure values ​​of each node and pressure drop of each pipe segment in the hydraulic parameters are then calibrated. After calibration, the updated pipeline flow distribution, temperature distribution, and hydraulic parameters are written into the state position of the next simulation cycle according to the incremented simulation time, forming a continuous correspondence under the incremented simulation time. After completing the continuous correspondence, the updated pipeline flow distribution, temperature distribution, and hydraulic parameters at the incremented simulation time are used as the current simulation state set for the next simulation cycle to continue the iterative calculation.

[0132] S5.3.3 When the incremented simulation time reaches the preset total simulation duration, the loop iteration ends.

[0133] Specifically, when the incremented simulation time reaches the preset total simulation time, the comparison result between the incremented simulation time and the preset total simulation time is first confirmed, and the pipeline flow distribution, temperature distribution, and hydraulic parameters of the incremented simulation time are finally calibrated. After the final time calibration is completed, the pipeline flow distribution, temperature distribution, and hydraulic parameters of the incremented simulation time are kept at the current iteration result position, and the current simulation state set update for the next round of simulation is no longer executed. After stopping the current simulation state set update for the next round of simulation, the subsequent loop calculation process of the incremented simulation time is terminated, and the loop iteration ends.

[0134] This embodiment also provides a hydraulic-thermal coupling prediction and simulation system for a centralized heating system, including: an initial setting module, used to collect pipeline parameter data, heat source operation data, heat exchange station operation data, and meteorological temperature data; initialize pipe segment flow rate and node temperature based on pipeline parameter data and set the initial simulation time; and generate an initial simulation condition set; a boundary setting module, used to set the heat source supply water temperature, heat exchange station target load, and outdoor air temperature as time-varying boundary conditions; and generate a boundary dataset using the initial simulation condition set; and a heat exchange prediction module, used to calculate the heat exchange prediction using the heat exchange station's heat exchange state. The solution method calculates the primary flow rate and primary return water temperature of the heat exchange station based on the boundary dataset to obtain the dynamic parameters of the heat exchange station. The hydraulic solution module uses a steady-state hydraulic analysis algorithm to solve the hydraulics of the pipeline network based on the dynamic parameters of the heat exchange station and generates the flow distribution of the pipe segment. The thermal iteration module uses a dynamic thermal calculation method to solve the pipeline temperature change based on the flow distribution of the pipe segment, the primary return water temperature of the heat exchange station and the supply water temperature of the heat source, updates the temperature distribution of all nodes in the network, and promotes the simulation cycle iteration to generate the flow distribution, temperature distribution and hydraulic operating parameters of the pipeline network.

[0135] In summary, this invention achieves continuous simulation of the operation process of a centralized heating system by using a dynamic thermodynamic calculation method to solve for pipeline temperature changes based on the flow distribution of pipe sections, the primary return water temperature of the heat exchange station, and the supply water temperature of the heat source. This method updates the temperature distribution of all nodes in the network and promotes iterative simulation at each simulation moment. It can be used to output the flow distribution, temperature distribution, and hydraulic parameters of the pipeline network, thereby improving the completeness of hydraulic-thermal coupling prediction and the practicality of simulation analysis.

[0136] It should be noted that the above embodiments are only used to illustrate the technical solutions of the present invention and are not intended to limit it. Although the present invention has been described in detail with reference to preferred embodiments, those skilled in the art should understand that modifications or equivalent substitutions can be made to the technical solutions of the present invention without departing from the spirit and scope of the technical solutions of the present invention, and all such modifications or substitutions should be covered within the scope of the claims of the present invention.

Claims

1. A hydraulic-thermal coupling prediction and simulation method for centralized heating systems, characterized in that: include, Collect pipeline network parameter data, heat source operation data, heat exchange station operation data, and meteorological temperature data. Initialize pipe segment flow and node temperature based on pipeline network parameter data, set the initial simulation time, and generate the initial simulation condition set. The heat source water supply temperature, the target load of the heat exchange station, and the outdoor air temperature are set as time-varying boundary conditions, and the boundary dataset is generated using the initial simulation condition set. The heat exchanger state solution method is adopted, and the primary side flow rate and primary side return water temperature of the heat exchanger are calculated based on the boundary dataset to obtain the dynamic parameters of the heat exchanger. The specific steps are as follows. The target load data of each heat exchange station at the current simulation moment is obtained from the boundary dataset, and the primary side water supply temperature data of each heat exchange station at the current simulation moment is obtained from the current simulation state set. At the same time, the secondary side water supply temperature set value and the target load data of the heat exchange station at the previous simulation moment are obtained from the heat exchange station operation data. The heat exchange station state solution method is adopted. Based on the primary side water supply temperature data, the secondary side water supply temperature setpoint, the target load data of the heat exchange station at the current simulation time and the target load data of the heat exchange station at the previous simulation time, the state equation of the heat exchange station is solved, the primary side flow rate and the primary side return water temperature of each heat exchange station are calculated, and the dynamic parameters of the heat exchange station are obtained. The specific steps for determining the primary side water supply temperature data and the secondary side water supply temperature setpoint of the heat exchange station are as follows. Extract the node temperature information at the current simulation moment from the node temperature information in the current simulation state set; According to the node number of the heat exchange station, the water supply temperature of each heat exchange station node at the current simulation time is extracted from the node temperature information at the current simulation time to form the primary side water supply temperature data of the heat exchange station. Extract the secondary water supply temperature setting information of the heat exchange station at the current simulation moment from the heat exchange station operation data; According to the heat exchange station node number, extract the heat exchange station secondary side water supply temperature setting value of each heat exchange station at the current simulation time from the heat exchange station secondary side water supply temperature setting information at the current simulation time. Using a steady-state hydraulic analysis algorithm, the hydraulics of the pipe network are solved based on the dynamic parameters of the heat exchange station to generate the flow distribution of the pipe segments. The specific steps are as follows. Based on the pipeline network topology, a basic correlation matrix and a basic loop matrix are constructed, and a set of hydraulic equations are constructed based on the conservation of nodal flow and the pressure balance of loops. Based on the dynamic parameters of the heat exchange station, an iterative initial value matrix is ​​generated, and a steady-state hydraulic analysis algorithm is used to iteratively solve the problem to obtain the flow rate and nodal pressure distribution of the pipe section. Based on the flow distribution of pipe sections, the primary return water temperature of the heat exchange station, and the supply water temperature of the heat source, a dynamic thermodynamic calculation method is used to solve for the pipe temperature change, update the temperature distribution of all nodes in the network, and advance the simulation iteratively to generate the network flow distribution, temperature distribution, and hydraulic parameters. The specific steps are as follows. Based on the flow distribution of the pipe section, the temperature distribution of all nodes in the network, the density of the medium, the specific heat capacity of the medium, the total heat transfer resistance, and the ambient temperature, dynamic calculation parameters for the pipe temperature are generated, and dynamic thermodynamic calculations are performed. The flow balance equations at the nodes and the pressure balance equations in the loops are combined, and the flow of each pipe segment is iteratively updated based on the flow balance results of each node and the pressure balance results of each closed loop to generate the flow distribution of the pipe segment. The updated pipeline flow and temperature distributions are used as the current simulation state set for the next round of simulation, and are executed cyclically according to the preset simulation time interval until the preset total simulation time is reached.

2. The hydraulic-thermal coupling prediction and simulation method for centralized heating systems as described in claim 1, characterized in that: The pipeline network parameter data includes the connection relationship of pipe segments, the number of nodes in each pipe segment, the impedance of pipe segments, the length of pipe segments, the diameter of pipe segments, the height difference of pipe segments, the thermal resistance of pipe segments, the pump parameters, the initial temperature, the node number of the heat exchange station, the node number of the heat source, the pressure of the heat source constant pressure point, the simulation duration and the simulation time interval. The heat source operation data includes heat source supply water temperature, heat source return water temperature, heat source inlet pressure, heat source outlet pressure, heat source circulation flow rate, heat source heating power, and heat source regulation setpoint. The heat exchange station operation data includes the heat exchange station target load, heat exchange station actual load, primary side supply water temperature, heat exchange station primary side return water temperature, primary side flow rate, secondary side supply water temperature, and secondary side return water temperature. The meteorological temperature data includes current outdoor temperature, historical outdoor temperature, and predicted outdoor temperature.

3. The hydraulic-thermal coupling prediction and simulation method for centralized heating systems as described in claim 2, characterized in that: The specific steps for generating the initial simulation condition set are as follows: Based on the pipeline network parameter data, a two-dimensional structure of the initial pipeline segment flow matrix is ​​established, where the first dimension represents the simulation time and the second dimension represents the pipeline segment number. The initial flow value of each pipeline segment is set to zero to generate the initial pipeline segment flow matrix. Based on the pipeline network parameter data, a three-dimensional structure of the initial moment node temperature matrix is ​​established, where the first dimension represents the simulation moment, the second dimension represents the pipe segment number, and the third dimension represents the node number inside the pipe segment. The temperature of each node is set to the initial temperature to generate the initial moment node temperature matrix. The initial simulation time is set to zero, and the initial pipe segment flow matrix and the initial time node temperature matrix are used as initialization inputs to generate the initial simulation condition set.

4. The hydraulic-thermal coupling prediction and simulation method for centralized heating systems as described in claim 3, characterized in that: The specific steps for generating the boundary dataset are as follows: Extract pipe segment flow rate information and node temperature information at the current simulation moment from the initial simulation condition set to establish the current simulation state set; Based on the heat source operation data, heat exchange station operation data and meteorological temperature data, time series matching and processing are performed to form the heat source water supply temperature time series, the heat exchange station target load time series and the outdoor temperature time series. Based on the current simulation state set, extract the heat source water supply temperature data at the current simulation moment from the heat source water supply temperature time series. Based on the current simulation state set, extract the target load data of the heat exchange station at the current simulation moment from the target load time series of the heat exchange station; Based on the current simulation state set, extract the outdoor temperature data at the current simulation moment from the outdoor temperature time series; The heat source water supply temperature data, heat exchange station target load data, and outdoor air temperature data are set as time-varying boundary conditions, and a boundary dataset is generated by combining them with the current simulation state set.

5. The hydraulic-thermal coupling prediction and simulation method for centralized heating systems as described in claim 4, characterized in that: The specific steps for generating the flow distribution in the pipe section are as follows. Based on the pipeline network parameter data, the connection relationship between each pipe segment and the topological nodes at both ends of the pipe segment is analyzed to form the pipeline network topology information, and the closed loop relationship formed by each pipe segment is identified to establish the basic correlation matrix and basic loop matrix of the pipeline network. Extract the pipe segment connection relationship and flow direction relationship of each node from the basic correlation matrix of the pipeline network, and establish the flow conservation relationship of each node based on the heat source circulation flow and the primary side flow of the heat exchange station node, and construct the node flow balance equation. Extract the pipe segment connection relationship and flow direction relationship of each closed loop from the basic loop matrix of the pipe network, and combine the pressure drop of each pipe segment and the pump head to form a pressure balance relationship, and construct the loop pressure balance equation; Based on the dynamic parameters of the heat exchange station and the flow balance equation of the nodes, the flow of the remaining branches in the pipe section is initially set, and the branch flow is calculated to generate the hydraulic iteration initial value matrix at the current moment. A steady-state hydraulic analysis algorithm is used, with the initial value matrix of the hydraulic iteration at the current moment as the initial value, to iteratively solve the nodal flow balance equation and the loop pressure balance equation to generate the flow distribution of the pipe segment.

6. The hydraulic-thermal coupling prediction and simulation method for centralized heating systems as described in claim 4, characterized in that: The specific steps for generating the pipeline network flow distribution, temperature distribution, and hydraulic operating parameters are as follows. The boundary dataset is updated based on the current simulation state set, and the heat exchange state of the heat exchange station is predicted based on the updated boundary dataset to obtain the heat exchange response result set of the station. Steady-state hydraulic calculations and dynamic thermal calculations are performed on the heat exchange response result set of the site to generate the flow distribution, temperature distribution and hydraulic operating parameters of the pipeline network. The simulation time is advanced according to the preset simulation time interval, and the simulation is iterated in cycles based on the flow distribution of pipe sections, the temperature distribution of nodes in the whole network, and the hydraulic operating parameters at each simulation time until the preset total simulation time is reached.

7. The hydraulic-thermal coupling prediction and simulation method for centralized heating systems as described in claim 6, characterized in that: The specific steps for advancing the current simulation time according to a preset simulation time interval are as follows. The current simulation time is incremented according to the simulation time interval, and it is determined whether the incremented simulation time reaches the preset total simulation duration. If the incremented simulation time does not reach the preset total simulation time, the updated pipeline flow distribution and temperature distribution will be used as the current simulation state set for the next round of simulation to continue the loop calculation. The loop iteration ends when the incremented simulation time reaches the preset total simulation duration.

8. A hydraulic-thermal coupling prediction and simulation system for a centralized heating system, based on the hydraulic-thermal coupling prediction and simulation method for a centralized heating system according to any one of claims 1 to 7, characterized in that: This includes an initial setup module, which is used to collect pipeline parameter data, heat source operation data, heat exchange station operation data and meteorological temperature data, initialize pipe segment flow and node temperature based on pipeline parameter data, set the initial simulation time, and generate an initial simulation condition set; The boundary setting module is used to set the heat source water supply temperature, the target load of the heat exchange station and the outdoor air temperature as time-varying boundary conditions, and generate a boundary dataset using the initial simulation condition set; The heat exchange prediction module is used to calculate the primary flow rate and primary return water temperature of the heat exchange station based on the boundary dataset using the heat exchange station heat exchange state solution method, and obtain the dynamic parameters of the heat exchange station. The hydraulic solution module is used to solve the hydraulic system of the pipe network based on the dynamic parameters of the heat exchange station using a steady-state hydraulic analysis algorithm, and generate the flow distribution of the pipe section. The thermal iteration module is used to solve the pipeline temperature change based on the pipeline flow distribution, the primary side return water temperature of the heat exchange station and the heat source supply water temperature, using a dynamic thermal calculation method, update the temperature distribution of all network nodes, and promote the simulation to iterate in cycles to generate the pipeline flow distribution, temperature distribution and hydraulic operating parameters.

Images