WO2023217979A1 - Procédé et système de détermination de perte de chaleur dans une conduite de desserte d'un système d'énergie collectif - Google Patents

Procédé et système de détermination de perte de chaleur dans une conduite de desserte d'un système d'énergie collectif Download PDF

Info

Publication number
WO2023217979A1
WO2023217979A1 PCT/EP2023/062638 EP2023062638W WO2023217979A1 WO 2023217979 A1 WO2023217979 A1 WO 2023217979A1 EP 2023062638 W EP2023062638 W EP 2023062638W WO 2023217979 A1 WO2023217979 A1 WO 2023217979A1
Authority
WO
WIPO (PCT)
Prior art keywords
service
equations
lines
fluid
temperature
Prior art date
Application number
PCT/EP2023/062638
Other languages
English (en)
Inventor
Brian Kongsgaard Nielsen
Carsten Skovmose Kallesøe
Original Assignee
Grundfos Holding A/S
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 Grundfos Holding A/S filed Critical Grundfos Holding A/S
Publication of WO2023217979A1 publication Critical patent/WO2023217979A1/fr

Links

Classifications

    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01K—MEASURING TEMPERATURE; MEASURING QUANTITY OF HEAT; THERMALLY-SENSITIVE ELEMENTS NOT OTHERWISE PROVIDED FOR
    • G01K17/00—Measuring quantity of heat
    • G01K17/06—Measuring quantity of heat conveyed by flowing media, e.g. in heating systems e.g. the quantity of heat in a transporting medium, delivered to or consumed in an expenditure device

Definitions

  • the present invention relates to the determination of heat loss in a service line of a district energy system, such as a district heating system or a cooling distribution network.
  • a fluid is transported from an energy supply plant, e.g. a district heating plant, via a grid of fluid lines to a plurality of consumers.
  • the district energy system comprises a network of supply lines and a plurality of service lines.
  • the service lines branch off from the supply lines and connect the respective consumers with the supply line.
  • the grid of fluid lines for distributing the fluid to a plurality of consumers will also be referred to as district energy grid.
  • a problem for the district energy grid operator is that, contrary to the supply lines, there are currently no good ways of inspecting the conditions of the service lines.
  • the service lines typically have no such electrical detection system.
  • the service lines are typically located on private property, which creates difficulties when performing thermal inspections or other inspection types.
  • a fault in a service line can for example be a service line where the insulation has deteriorated or, even worse, a service line leaking district energy water. Therefore, it would be desirable if such an inspection can be done remotely. Moreover it would be desirable when faults in a service line can be detected early, in particular before leakage occurs.
  • EP3531368 discloses a method for determining a flow rate and a temperature of a fluid at a selected position in a district heating or cooling utility distribution network, comprising a plurality of interconnected distribution lines and smart utility meters.
  • the meters are arranged to register the energy delivered at consumer premises situated along the distribution lines, the method comprising the steps of: collecting meter data time series from the smart utility meters using an Advanced Metering Infrastructure, and calculating the temperature and the flow rate at the selected position.
  • the calculations are based on flow and temperature information derived from the collected meter data time series, the topology of the utility distribution network and heat transfer coefficients of the distribution lines.
  • the meter data time series comprises the accumulated volume of fluid delivered to the consumer, and the integrated flow-temperature product calculated by the meter.
  • a computer-implemented method for determining thermal loss in a district energy grid comprising one or more supply lines and a plurality of service lines for transporting a fluid to a plurality of consumers, each supply line feeding the fluid into a plurality of service lines and each service line connecting an consumer with one of the one or more supply lines, the method comprising:
  • the set of equations comprising a plurality of equations relating the observed fluid temperatures and the observed fluid flows with thermal loss parameters indicative of thermal losses of respective ones of said selected one or more service lines, and
  • the losses of some or all of the service lines of the district energy grid can be calculated directly from sensor data obtainable from the consumers' heat meters or other suitable logging devices.
  • Embodiments of the method and system disclosed herein do not require detailed knowledge of a network model or the modelling of temperature and flow at respective nodes of the main grid.
  • each of the selected service lines has a logging device associated with it, such that the associated logging device is configured to record sensor data pertaining to the service line with which the logging device is associated. It will be appreciated that the method may determine thermal energy losses for all or for only some service lines of the district energy grid, e.g.
  • Embodiments of the method and system described herein provide the district energy grid operator with information about service pipe losses in sub areas of the grid. Accordingly, the district grid operator can use that information to plan maintenance ahead, even when only modest knowledge of the network is available.
  • Various embodiments of the method disclosed herein allow early detection of faults in service lines, e.g. faulty insulation which might otherwise result in pipe corrosion.
  • the thermal loss parameters of the respective service lines are unknown variables of the set of equations.
  • Computing a solution of the set of equations may thus comprise numerically solving the set of equations for the thermal loss parameters.
  • the plurality of equations are linear in the service loss parameters, thus allowing for a computationally efficient solution of the equations.
  • the plurality of equations relate instant heat losses in the respective service lines to the fluid flow in the respective service pipes and a temperature difference between the observed fluid temperatures, observed for the respective service lines, and a fluid temperature in the supply line.
  • the fluid temperature of the supply line may be a further unknown of the set of equations.
  • the set of equations is an overdetermined set of equations and the solution of the set of equations is an approximate solution of the overdetermined set of equations.
  • the number of equations in the set of equations may be larger than the number of thermal loss parameters and, in particular, larger than the number of selected service lines.
  • Computing the solution of the set of equations may thus comprise computing an approximate solution, in particular a least-squares solution, to the overdetermined set of equations.
  • computing the solution of the set of equations may comprise applying a least-square method for determining the approximate solution of the overdetermined set of equations.
  • computing the solution of the set of equations comprises computing the solution as a solution to a constrained optimization problem.
  • computing the solution of the set of equations comprises minimizing or maximizing a cost function over at least one additional temperature parameter indicative of a temperature and/or a temperature loss associated with the supply line.
  • computing the solution of the set of equations comprises minimizing or maximizing a cost function over at least one ground temperature parameter that is associated with a ground temperature. Accordingly, the computation of the set of equations does not require prior knowledge of the temperature and/or a temperature loss associated with the supply line and or prior knowledge of the ground temperature.
  • the set of equations includes equations relating to different service lines, thus allowing an assessment of the state of multiple service lines and providing a computed estimate that is less sensitive to noise.
  • the one or more selected service lines include a plurality of service lines, and wherein each of the plurality of equations relates a thermal loss parameter of a respective service line with an observed fluid temperature and an observed fluid flow associated with said respective service line.
  • the selected plurality of service lines are all connected to a single supply line, in particular to a section of a single supply line that does not include any branch points other than the service lines branching off to individual consumers.
  • the plurality of equations includes, for some or all of the selected one or more service lines, multiple equations associated with respective points in time within an observation period of time, thereby providing a more reliable estimate for the thermal losses during said period of time.
  • the set of equations comprises a first equation relating an observed fluid temperature and an observed fluid flow of the first service line at a first point in time with the first thermal loss parameter.
  • the plurality of equations further comprises another equation associated with a second point in time and relating an observed fluid temperature and an observed fluid flow of the first service line at said second point in time with said first thermal loss parameter, the second point in time being different from the first point in time.
  • the one or more selected service lines include a first and a second service line, different from the first service line, and the plurality of equations comprises another equation, associated with the second service line and relating an observed fluid temperature and an observed fluid flow of said second service line with a second thermal loss parameter indicative of a thermal loss of said second service line.
  • the first equation and the another equation associated with the second service line each include a temperature parameter associated with the supply line, the temperature parameter representing an unknown variable of the first equation and of the another equation associated with the second service line, i.e. the computational model may exploit the fact that the selected service lines are all connected to the same supply line or otherwise such that the temperature in the supply line is substantially uniform across the selected service lines.
  • the present disclosure relates to different aspects including the method described above and in the following, corresponding apparatus, systems, methods, and/or products, each yielding one or more of the benefits and advantages described in connection with one or more of the other aspects, and each having one or more embodiments corresponding to the embodiments described in connection with one or more of the other aspects and/or disclosed in the appended claims.
  • embodiments of the method disclosed herein may be computer- implemented. Accordingly, disclosed herein are embodiments of a data processing system configured to perform the steps of the method described herein.
  • the data processing system may have stored thereon program code adapted to cause, when executed by the data processing system, the data processing system to perform the steps of the method described herein.
  • the data processing system may be embodied as a single computer or other data processing device, or as a distributed system including multiple computers and/or other data processing devices, e.g. a client-server system, a cloudbased system, etc.
  • the data processing system may include a data storage device for storing the computer program and sensor data.
  • the data processing system may include a communications interface for receiving sensor data from consumer's heat meters or other suitable logging devices and/or sensors associated with the individual service lines.
  • the data processing system may receive the sensor data from the heat meters and/or other logging devices and/or sensors via a suitable wired or wireless communicative connection, e.g. via a suitable communications network.
  • the data processing system may provide a user-interface for allowing a user to monitor the computed thermal loss parameters and/or other data associated with the conditions of individual service lines.
  • the data processing system may also issue warnings or alerts or other notifications responsive to detected losses or other conditions, e.g. audible or visual alerts, alerts communicated via e-mail, SMS, or other forms of notifications, and/or the like.
  • the district heating or cooling system comprises:
  • the district energy grid comprising one or more supply lines and a plurality of service lines for transporting a fluid to a plurality of consumers, each supply line feeding fluid into a plurality of service lines and each service line connecting a consumer with one of the one or more supply lines,
  • a plurality of sensors configured to output sensor data indicative of measured fluid temperatures and associated fluid flow of the transported fluid received at respective consumers via the respective service lines
  • the sensor data indicative of the measured fluid temperatures and the associated fluid flow may represent the fluid temperatures and associated fluid flow directly.
  • the fluid temperatures and the associated fluid flow may be derivable from the sensor data.
  • the sensor data indicative of the fluid flow may be sensor data directly representing the fluid flow or it may be sensor data from which the fluid flow can be derived, e.g. sensor data representing received volumes of fluid.
  • a computer program may comprise program code means adapted to cause a data processing system to perform the acts of the method disclosed above and in the following when the program code means are executed on the data processing system.
  • the computer program may be stored on a computer-readable storage medium, in particular a non-transient storage medium, or embodied as a data signal.
  • the non-transient storage medium may comprise any suitable circuitry or device for storing data, such as a RAM, a ROM, an EPROM, EEPROM, flash memory, magnetic or optical storage device, such as a CD ROM, a DVD, a hard disk, and/or the like.
  • FIG. 1 schematically shows an embodiment of a district energy system.
  • FIG. 2 schematically shows another embodiment of a district energy system.
  • FIG. 3 shows a flow diagram of an example of method for determining thermal loss in a district energy grid.
  • FIG. 4 shows a flow diagram of a method for computing estimated thermal loss parameters as a solution of an optimization problem.
  • FIGs. 5A-F illustrate examples of measured fluid temperatures and flows and examples of results computed by an embodiment of the method described herein.
  • FIG. 1 schematically shows an embodiment of a district energy system, generally designated by reference numeral 1.
  • the district energy system comprises a district energy grid 2, which comprises one or more supply lines 6 and a plurality of service lines 5 for transporting a fluid to a plurality of consumers 4.
  • Each supply line 6 feeds the fluid into a plurality of service lines 5 and each service line 5 connects a consumer 4 with one of the one or more supply lines 6.
  • the supply lines 6 and/or the service lines may be pipes or other suitable fluid conduits, preferably thermally insulated pipes.
  • the fluid may be water or another suitable medium facilitating heat transfer by advection of the medium through the district energy grid.
  • the district energy system comprises an energy supply plant 3.
  • the district energy system is a district heating system and the energy supply plant is a district heating plant which feeds heated liquid into the district energy grid 2.
  • the district energy grid may comprise a network of supply lines, also referred to as main grid, e.g. as illustrated in FIG. 2.
  • Each consumer 4 is connected to the main grid by a service line 5, e.g. by an insulated service pipe.
  • the main grid has one or more supply lines 6 and one or more return lines 7.
  • the supply lines form a supply path for feeding the fluid from the energy supply plant 3 to the individual consumers 4 and the return lines 7 define a return path for returning fluid from the consumers to the energy supply plant 3.
  • the district energy system further comprises a plurality of logging devices 8.
  • the logging devices 8 are configured to provide temperature sensor data and flow sensor data.
  • the temperature data may be indicative of the measured fluid temperature of the fluid arriving at the respective consumers, also referred to as the supply temperature
  • the service line has a receiving end and a consumer end, opposite the receiving end.
  • the service line is fluidly connected to the supply line and receives fluid from the supply line at the receiving end.
  • the consumer receives the fluid via the service line at the consumer end of the service line.
  • the supply temperature represented by the sensor data is indicative of the fluid temperature of the fluid arriving at the consumer at the consumer end of the service line.
  • the flow sensor data may be sensor data directly representing the fluid flow or it may be sensor data from which the fluid flow can be derived.
  • the flow sensor data may represent volume data indicative of the accumulated volume of fluid received by the consumer through the service line.
  • the fluid flow may then be derived as a derivative of the volume data with respect to time.
  • the logging devices are provided in the form of heat meters, that are installed at the respective consumers 4 and that are configured to log the supply temperature and the accumulated volume of fluid received by the respective consumers. It will be appreciated that, in alternative embodiments, the logging devices may include other types of sensors, e.g. they may be embodied as individual and separate sensors for measuring temperature and volume/flow, respectively.
  • the district energy system further comprises a data processing system 9, e.g. a server computer, a network of multiple computers, e.g. a client-server architecture, a distributed computing architecture, a cloud-based computing environment, and/or the like.
  • the data processing system 9 may be located at the energy supply plant 3 or at a different location.
  • the data processing system 9 receives logged sensor data from the logging devices 8.
  • the logging devices 8 at the respective consumers are communicatively coupled to the data processing system 9, e.g. directly or indirectly.
  • the logging devices 8 may communicate the logged sensor data to the data processing system via wired and/or wireless connections, e.g. via a suitable communications network, such as a cellular telecommunications network, via the internet or another suitable computer network.
  • many modern heat meters are configured for automatic communication of the logged data to a central data processing system.
  • the logging devices may communicate the sensor data continuously or intermittently, e.g. one measurement at a time or as batches of multiple, time-stamped sets of sensor data.
  • the data processing system obtains the logged sensor data from a suitable data repository where the logged sensor data has previously been stored, e.g. by the heat meters, by an energy grid maintenance or surveillance system, or the like.
  • the data processing system 9 receives the logged sensor data and is programmed to process the received sensor data to determine thermal loss in the district energy grid as described herein. Embodiments of a method implemented by the data processing system will be described in more detail below with reference to FIGs. 3 and 4.
  • the data processing system may store the received sensor data and the results of the determination of the energy loss. Alternatively or additionally, the data processing system may implement functionality for presenting the sensor data and/or the results of the energy loss determination to a user of the data processing system, e.g. in the form of graphs, tables or the like.
  • FIG. 2 schematically shows another embodiment of a district energy system 1.
  • the district energy system 1 of FIG. 2 is similar to the district energy system of FIG.
  • the logging devices 8 are preferably communicatively coupled to the data processing system 9 and configured to transmit logged sensor data to the data processing system.
  • sensor data may be read out and entered into the data processing system in a different manner, e.g. by storing the logged sensor data on a portable data carrier, and by reading the stored data by the data processing system from the portable data carrier.
  • a communication of the logged sensor data via a suitable communications, link allows a more efficient and, optionally, even (quasi) realtime transmission and analysis of the logged sensor data.
  • FIG. 2 illustrates that the energy grid 2 typically includes more than a single supply line, in particular a network of supply lines 6 that are interconnected at nodes 20.
  • the nodes 20 typically include valves, which allow the energy grid operator to direct the fluid flow along respective paths through the grid, e.g. so as to temporarily isolate certain parts of the grid from each other for maintenance, or the like.
  • open valves are designated V and closed valves are designated C.
  • the energy grid 2 includes multiple subareas 2a - 2e where all consumers in a given subarea are connected to a single supply line. In each subarea of the grid, there are multiple logging devices 8.
  • Embodiments of a process or determining thermal loss in a district energy grid will now be described with reference to FIG. 3 and with continued reference to FIGs. 1 and 2.
  • the process may be performed by the data processing system 9.
  • some embodiments of the method may be implemented by a cloud computing environment or another suitable data processing system.
  • the determined thermal loss may be calculated and represented as a suitable loss figure of merit indicative of the thermal loss.
  • the loss figure of merit will also be referred to as thermal loss parameter. It may e.g. be output in the form of a graphical indicator of that figure of merit, e.g. via an internet browser or another suitable user interface.
  • FIG. 3 shows a flow diagram of an example of a computer-implemented process for determining thermal loss in a district energy grid.
  • the process groups the service lines and the associated logging devices 8 of the district energy grid 2 into respective groups.
  • the process may set up a data structure representing the district energy grid and identifying respective selected groups of data logging devices.
  • the grouping may be performed manually, automatically or in a semi-automatic, user-assisted manner. The grouping is typically only performed once or when changes to the structure of the energy grid occur.
  • various embodiments of the process use data from a selected group of logging devices 8 within a geographically limited region of the energy grid 2 where the supply line temperature of the fluid fed from the main district energy grid into the service lines can be expected to be roughly the same for all service lines or where differences in the supply line temperature between users can easily be modelled.
  • a geographically limited region will also be referred to as a "similar-temperature area”.
  • Groups of logging devices within the energy grid may be selected in a number of ways, e.g. based on existing information about the energy grid, such as:
  • the grouping may be based on each of the above information, individually or on a combination of two or even all three types of information. It will further be appreciated that the grouping may be performed in a different manner, including based on other types of information.
  • the grouping based on street name and/or consumer address coordinates is easy to implement but is also expected to be the most inaccurate. GPS location of the logging devices is sometimes registered by the grid operator upon installation of the logging device. If available, this data may provide valuable additional information for the grouping, thus increasing the accuracy of the grouping.
  • Some district energy grid operators maintain valve closing maps in their GIS data. Accordingly, if such data is available, it may also provide a useful basis for the grouping of logging devices. Moreover, such maps may give accurate information about the flow direction in the grid, and thus help identifying similar temperature areas.
  • the logging devices in each of the sub-areas 2a - 2d may be selected as a respective group representing a similar temperature area, as all service pipes within each group are connected to a linear string of supply line without any branches.
  • several of the sub-areas may be combined into a single group.
  • a subarea may be divided into multiple groups, e.g. if the subarea includes a very long string of supply line.
  • a graph theoretical network model can be established from the district energy grid operator's GIS data. Even though the grid operator in many cases does not have detailed pipe dimensions and heat loss coefficients available they often do have a map of their pipes that may be transferred to a graph theoretical model. It is an advantage of the method disclosed herein, that it does not rely on detailed information about the energy grid, such as pipe dimensions, heat loss coefficients of the pipes etc. As can be appreciated from the above, even though access to GIS data or other grid topology data may increase the accuracy of the grouping, and thus the accuracy of the resulting heat loss computation, the grouping of logging devices may even be performed without any topology information about the grid; it may entirely be based on postal addresses of the consumers. In subsequent step S2, the process selects a group of service lines of the district energy grid 1, each selected service line having an associated data logging device 8.
  • various embodiments of the method disclosed herein perform a determination of energy loss for a thus selected group of logging devices, i.e. for a selected group of service pipes, preferably a group of service pipes connected to a single supply line or otherwise connected to the supply grid such that the supply lines of the selected group receive fluid of approximately the same or at least of similar temperature or where temperature differences between the points of the supply line where the service lines are connected can easily be modelled.
  • the selected group of service lines includes between 5 and 50 service lines with associated logging devices, a typical group size is 10 to 30 service lines. Larger groups of service lines increase the statistical accuracy of the computation, but also increase the computational complexity.
  • step S3 the process selects a suitable observation period for which thermal losses of one or more service lines of the selected group of service lines are to be computed.
  • the length of the observation period may be selected in dependence of the rate at which the logging devices record sensor data. Other factors that may influence the choice of observation period may include the size of the selected group of service lines. For example, if the group is small, the accuracy of the estimates of the heat loss may decrease. Increasing the length of the observation period may then increase the accuracy of the estimates. When the observation period is selected to be very long, the assumption of uniform supply line temperature may be less accurate.
  • step S4 the process obtains logged sensor data from the logging devices associated with the selected group of service lines.
  • the process may obtain the sensor data in a variety of ways. For example, the process may receive the sensor data directly or indirectly from the logging devices. The process may receive sensor data from an entity that transmits the sensor data to the process, or the process may request the sensor data from another entity, or retrieve the sensor data from a suitable data repository.
  • Various embodiments of the method use sensor data logged at different points in time over a selected observation period of time, e.g. a predetermined observation period, such as over an observation period of between lh and 1 week, e.g. between 6 h and 24 h, e.g. between 12 h and 24 h.
  • a predetermined observation period such as over an observation period of between lh and 1 week, e.g. between 6 h and 24 h, e.g. between 12 h and 24 h.
  • various embodiments use, for each logging device of the selected group, logged sensor data, logged at different points in time over an observation period of time.
  • the duration of the observation period may be selected long enough so that a plurality of data points logged at respective points in time can be collected, i.e. the choice of observation period may depend on the frequency at which the logging devices log measurements.
  • each logging device may log the sensor data at a rate between once every minute and once per day, such as between once every 10 minutes and once every 6h, such as hourly. It will be appreciated that some embodiments, e.g. when the sensor data is logged at a high rate, may accumulate, e.g. average, data over a sampling period and use the thus accumulated data as a logged measurement value.
  • the length of the observation period is preferably chosen short enough such that the temperature in the supply line may be regarded as substantially constant over the observation period.
  • a larger number of data points may increase the statistical accuracy of the determined thermal losses.
  • a larger number of data points also increases the computational complexity of the process.
  • the use of a large data batch from multiple logging devices and at multiple points in time improves robustness of the heat loss determination.
  • the process may obtain the sensor data in real time from the logging devices or it may retrieve one or more batches of previously logged data for a desired observation period.
  • step S5 the process defines and solves an optimization problem, the optimization problem is formulated such that its solution directly represents the thermal losses of the service lines of the selected group for the selected observation period, in particular without the need for subsequent dynamic modelling of fluid flow through the main grid and without use of measured or computed data pertaining to parts of the energy grid outside the subarea defined by the selected group of service lines.
  • the process calculates the thermal loss of a service line 5 between a supply line 6 and a logging device 8. The calculation may be based on the supply temperature and the volume and/or flow data logs of one or more logging devices 8.
  • the logging device may conveniently be a heat meter.
  • Heat meters normally log instant supply and return temperatures across a consumer installation, and/or flow weighted supply and return temperatures.
  • the flow weighted temperature is sometimes referred to as the integrated flow temperature product (ITFP).
  • ITFP integrated flow temperature product
  • the temperature data will simply be referred to either as supply temperature and return temperature or commonly as temperature data.
  • accumulated volume and/or instant fluid flow rate is logged by the heat meter.
  • the parameter C is a constant known from water properties.
  • Logs of T m and q can be obtained from the heat meter or another form of logging device, either directly or by suitable pre-processing, depending on the type of sensor data logged by the logging device.
  • Various embodiments of the process described herein estimate Q, or a thermal loss parameter related to Q, based on the logs of T m and q for a selected group of logging devices and, hence, for a selected group of service lines. Examples of the calculation process will be described in more detail with reference to FIG.
  • step S6 the process outputs the computed thermal losses or another loss figure of merit indicative of thermal losses of the selected service lines.
  • the process computes the thermal losses for respective observation periods, e.g. by computing the thermal loss for a new observation period as soon as new logged data is available, e.g. so as to provide estimated recent thermal losses for a sliding window.
  • thermal losses of the one or more service lines of the selected group of service lines are to be computed for further observation periods, e.g. so as to assess how the thermal losses develop over time, the process returns to step S4 where a different observation period is selected.
  • step S3 a different group of service lines is selected, e.g. so as to allow assessment and/or comparison of multiple subareas of the energy grid.
  • the process may process all selected groups for a single observation period before optionally proceeding to another observation period.
  • FIG. 4 shows a schematic flow diagram of a process of setting up and solving an optimization problem for estimating thermal losses of one or more selected service lines based on the logs of T m and q for a selected group of logging devices and for a selected observation time.
  • the process sets up suitable data structures for representing the optimization problem to be solved and populates the data structures with the fluid temperatures T m and the associated fluid flow q of the incoming fluid arriving at the logging device at respective points in time within the observation period.
  • the obtained sensor data may include time-stamped measurements from which the fluid temperature T m and the fluid flow q of the incoming fluid arriving at the logging device via the associated service line can be derived.
  • the process may receive the measured temperature T m and the measured fluid flow q from the corresponding logging device.
  • the process derives these quantities from the sensor data received from the logging device. For example, in some embodiments, the logging device may only measure the accumulated volume and not the instantaneous flow, or the flow measurements received from the logging device may be inaccurate.
  • the process may derive the fluid flow and/or the fluid temperature from the received sensor data.
  • the process may calculate the fluid flow as the gradient with respect to time of the accumulated volume using the actual time stamps of the received volume data.
  • the used gradient method is preferably a central difference method, but backward or forward difference or other methods may be used instead.
  • the received sensor data may be pre-processed in a different manner in order to derive the fluid temperature and/or the fluid flow.
  • the created data structure may be a matrix or another representation suitable for processing the optimization problem.
  • the process may arrange the data representing the fluid flow pertaining to each logging device in a corresponding time vector or other suitable representation of a time series of data points. For example, this may be done by moving the actual time stamps of the received or derived data points to the nearest time stamp which is and integer multiple of a predetermined sampling time Ts.
  • the logging rate of the logging devices may differ from logging device to logging device or even from sample to sample.
  • the selected sampling time Ts may e.g. be selected so as to match the average rate at which the logging devices log the sensor data. In one embodiment, Ts is selected to be 1 hour, but other sampling times may be used.
  • the alignment of the data to time stamps may be omitted. Accordingly, some embodiments merely load flow and temperature data pertaining to corresponding points in time within the observation period into one or more suitable data structures. Accordingly, the process is based on batches of sensor data, each batch including measurement data from one or more logging devices associated with a given measurement time. The number of batches will be denoted s. Preferably, there are a plurality of batches (s > 1), i.e. the process uses sensor data associated with multiple different points in time within an observation period.
  • the process processes the created time vector such as to remove equal time stamps for the same logging device.
  • P [P 1; P 2 , ... , P N ] T is a parameter vector where each parameter Pi relates to the thermal loss of a respective service line /.
  • the model is represented as a system of equations, linear in the thermal loss parameters Pi, which are the unknowns of the system.
  • the parameters in P may, for example, relate directly to losses Q in [W] (Watt), or to losses in [W/m] (Watt/meter), or to losses in [W/m/k] (Watt/meter/Kelvin) or to losses in [W/K] (Watt/Kelvin). Accordingly, the thermal loss parameters in P generally represent suitable thermal-loss-related figures of merit.
  • the parameter vector P may, in addition to the P ⁇ to P N service line loss related parameters, have one or more other parameters not relating to service line losses.
  • such one or more additional parameters may for example contain a temperature parameter such as T v in (eq.l).
  • step S53 the process numerically solves the model for the parameter vector P, which is returned as a result of the process.
  • the exact form of the data matrix X, the input vector y and the parameter vector P may differ for different embodiments and depend on the exact model used.
  • the input vector depends on the fluid temperature T v in the supply line, which is normally unknown.
  • the temperature T v is assumed to be equal for all service lines of the selected group. Accordingly, as described above, the group of logging devices is preferably selected such that this assumption is at least approximately fulfilled. Accordingly, an approximate solution may be found by formulating a constrained optimization problem for T v .
  • the model in (eq.2) may be solved by minimizing the squared error under the constraints that the losses are positive, that is
  • the resulting least-square solution P represents the estimated heat losses of the set of service lines during the observation time window.
  • a new data batch may be collected.
  • the solution to (eq.3) can be found by known numerical optimization algorithms.
  • the same model as above may be solved as an unconstrained least squares problem by moving the unknown temperature T v in the supply line into the parameter vector.
  • the formulation of the problem is then where and the solution to the parameter vector is calculated directly
  • the loss calculation in (eq.l) can be extended by adding some structure of the loss of the service line, e.g. according to the following modified model:
  • T g A loss proportional to the arithmetic mean temperature difference between the service line fluid and the surrounding temperature T g is assumed in (eq.5).
  • the surrounding temperature would be the ground temperature.
  • T g is now a needed input. It may be found via an outdoor temperature log and an associated model for the ground temperature, or by setting it as a predetermined constant, or by letting it be solved by the optimizer as shown below.
  • the parameter vector P for this model formulation now contains loss parameters B t with the unit [W/K],
  • T g is also a parameter to be found by the optimizer, i.e.
  • the above embodiments do not require that all logging devices of the selected group of logging devices, e.g. all heat meters in a "similar temperature area,” deliver sensor data at all time stamps. Accordingly, the process may not be able to determine the flow in the supply line feeding into the selected group of service lines.
  • the selected group of logging devices may include a logging device associated with each service line branching off a particular selected branch of the supply line, e.g. all logging devices of subarea 2b in the example of FIG. 2.
  • the process can determine the flow in the section of supply line from which the selected service lines branch off. This information can be exploited by the process to estimate the condition of that particular section of the supply line as well, e.g. by setting up a model based on the temperature decay in the supply line and the loss of the service line, e.g. as follows:
  • a is the loss parameter of the supply line and q t j is the flow in the supply line between consumers i and j.
  • This flow can be calculated exactly, if all logging devices deliver data at all time stamps, and there are no unknown by-pass flow between the supply line and the return path.
  • the model formulation is usable if these lengths are known.
  • the present embodiment estimates conditions of the supply line and the service pipes directly. solution to (eq.9) may be found in a similar manner as for the previous embodiments, e.g. as
  • the supply line loss model used in (eq.8) may be combined with the service line loss model used in (eq.5) to an optimization problem with a similar structure as (eq.9), i.e. having outer variables adhering to line loss models that are optimized and the service line losses are calculated as constraints as in (eq.3), (eq.7) or (eq.10).
  • the loss parameters calculated in (eq.3), (eq.7) or (eq.10) can easily be converted to a per-length figure of merit. This may give a more easily comparable figure when comparing service line losses between consumers as such pipes are sometimes of very different lengths within the same "similar temperature" area.
  • thermal loss of a service line in new condition is also normally given as a per-length figure of merit.
  • the service line lengths of are often not available.
  • embodiments of the process disclosed herein can assist the energy grid operator in identifying the service lines that have the highest losses, which is still very usable information for the district energy grid operator.
  • a method for estimating heat loss parameters of a district energy grid where the estimated heat loss parameters are indicative of, in particular directly related to, thermal losses of service lines.
  • Some embodiments of the method comprise the steps of: selecting a group of N logging devices, in particular logging devices within a "similar temperature area;" for example, the boundaries of the "similar temperature area” may be selected based on a street name, GIS data describing the district energy network, knowledge of valve locations in the district energy network, or based on other suitable information.
  • time-stamped sensor data logs from the N logging devices into a data structure, deriving N flow rates associated with the N logging devices from the received sensor data logs, constructing a time vector for each logging device; this may e.g. be done by moving the actual time stamps of the logging device to the nearest time stamp which is and integer multiple of a sampling time T s .
  • T s is preferably 1 hour or less, optionally processing the constructed time vector so there are no equal time stamps, for an observation time, e.g. an observation time of s sampling times, e.g.
  • the data matrix and the input vector preferably contain data for multiple time stamps of the constructed time vector.
  • Numerically solving the constructed model may comprise solving a constrained optimization problem. Solving the constrained optimization model may involve minimizing or maximizing a cost function over at least one additional parameter that is associated with a temperature and/or a loss parameter that is associated with a supply line. Alternatively or additionally, solving the constrained optimization model may involve minimizing or maximizing a cost function over at least one additional parameter that is associated with a ground temperature.
  • FIGs. 5A-F illustrate examples of measured fluid temperatures and flows of respective service lines, and examples of computed results obtained by an embodiment of the method described herein.
  • the selected group included four logging devices, labelled "Meter 1", ..., “Meter 4", respectively. It will be appreciated that other examples of selected groups of logging devices may include a different number of devices.
  • FIG. 5A shows logged temperature data for a selected group of logging devices and the estimated temperature 12 in the supply line, which has been computed by the method disclosed herein.
  • the computed supply line temperature is denoted "Estimated Main”.
  • FIG. 5B shows the logged flow data from the selected group of logging devices.
  • FIGs. 5C though 5F show resulting values of an estimated heat loss parameter for the respective service lines associated with the respective logging devices Meter 1 through Meter 4, respectively, including indications of +/- 3 standard deviations.
  • one of the logging devices (“Meter 4") has registered a relatively high flow 11 in the associated service line, while the corresponding temperature 10 registered by said logging device was relatively low. This indicates a high thermal loss of the corresponding service line which embodiments of the method disclosed herein help to identify.
  • FIG. 5F shows that the service line associated with Meter 4 has an unusually high heat loss.
  • the various embodiments described herein are applicable to thermal energy carrying pipe networks where a group of consumers are supplied from a common supply line, and where a temperature log of the flow out of the service line connecting the consumer to the supply line is available together with a log of the amount of medium flowing in the service line.
  • Examples of data logs relating to the amount of medium flowing in the service line are logs of the accumulated volume or logs of the instantaneous flow rate.
  • the thermal pipe loss or a thermal pipe loss related coefficient of the service line can be estimated by the various embodiments disclosed herein.
  • embodiments disclosed herein may be applied to energy grids of different types and sizes, including to energy grids located inside or outside of buildings. It will be appreciated that the application of the various embodiments typically is most useful in respect of energy grids where visual pipe inspection is difficult. Therefore, a district energy grid with automated data logging at the consumers and where the pipes are buried in the ground is a preferred application for embodiments of the invention disclosed herein.
  • Embodiment 1 A computer-implemented method for determining thermal loss in a district energy grid, the district energy grid comprising one or more supply lines and a plurality of service lines for transporting a fluid to a plurality of consumers, each supply line feeding the fluid into a plurality of service lines and each service line connecting a consumer with one of the one or more supply lines, the method comprising:
  • the set of equations comprising a plurality of equations relating the observed fluid temperatures and the observed fluid flows with thermal loss parameters indicative of thermal losses of respective ones of said selected one or more service lines, and
  • Embodiment 2 A method according to embodiment 1, wherein the thermal loss parameters represent unknown variables of the set of equations.
  • Embodiment 3 A method according to any one of the preceding embodiments, wherein the plurality of equations are linear in the thermal loss parameters.
  • Embodiment 4 A method according to any one of the preceding embodiments, wherein the set of equations is an overdetermined set of equations and wherein the solution of the set of equations is an approximate solution of the overdetermined set of equations.
  • Embodiment 5 A method according to embodiment 4, wherein computing the solution of the set of equations comprises computing the solution as a solution to a constrained optimization problem.
  • Embodiment 6 A method according to embodiment 4 or 5, wherein computing the solution of the set of equations comprises applying a least-square method for determining the approximate solution of the overdetermined set of equations.
  • Embodiment 7 A method according to embodiment 5 or 6, wherein computing the solution of the set of equations comprises minimizing or maximizing a cost function over at least one additional temperature parameter indicative of a temperature and/or a temperature loss associated with the supply line.
  • Embodiment 8 A method according to any one of embodiments 5 through 7 , wherein computing the solution of the set of equations comprises minimizing or maximizing a cost function over at least one ground temperature parameter that is associated with a ground temperature.
  • Embodiment 9 A method according to any one of the preceding embodiments, wherein the one or more selected service lines include a plurality of service lines, and wherein each of the plurality of equations relates a thermal loss parameter of a respective service line with an observed fluid temperature and an observed fluid flow associated with said respective service line.
  • Embodiment 10 A method according to any one of the preceding embodiments, wherein the set of equations comprises a first equation relating an observed fluid temperature and an observed fluid flow of the first service line at a first point in time with the first thermal loss parameter.
  • Embodiment 11 A method according to embodiment 10, wherein the plurality of equations further comprises another equation associated with a second point in time and relating an observed fluid temperature and an observed fluid flow of the first service line at said second point in time with said first thermal loss parameter, the second point in time being different from the first point in time.
  • Embodiment 12 A method according to embodiment 10 or 11, wherein the one or more selected service lines include a first and a second service line, different from the first service line, and wherein the plurality of equations comprises another equation associated with the second service line and relating an observed fluid temperature and an observed fluid flow of said second service line with a second thermal loss parameter indicative of a thermal loss of said second service line.
  • Embodiment 13 A method according to embodiment 12, wherein the first equation and the another equation associated with the second service line each include a temperature parameter associated with the supply line, the temperature parameter representing an unknown variable of the first equation and of the another equation associated with the second service line.
  • Embodiment 14 A method according to any one of the preceding embodiments, wherein the one or more selected service lines include a plurality of service lines all connected to a single supply line.
  • Embodiment 15 A data processing system configured to perform the steps of the method defined in any one of embodiments 1 through 14.
  • Embodiment 16 A district heating or cooling system comprising:
  • the district energy grid comprising one or more supply lines and a plurality of service lines for transporting a fluid to a plurality of consumers, each supply line feeding fluid into a plurality of service lines and each service line connecting a consumer with one of the one or more supply lines, - a plurality of sensors configured to output sensor data indicative of measured fluid temperatures and associated fluid flow of the transported fluid received at respective consumers via the respective service lines,
  • Embodiment 17 A computer program comprising computer program code configured, when executed by a data processing system, to cause the data processing system to perform the steps of the method according to any one of embodiments 1 through 14.
  • Embodiments of the method described herein may be computer-implemented.
  • embodiments of the method may be implemented by means of hardware comprising several distinct elements, and/or at least in part by means of a suitably programmed data processing system.
  • apparatus claims enumerating several means, several of these means can be embodied by one and the same element, component or item of hardware.
  • the mere fact that certain measures are recited in mutually different dependent claims or described in different embodiments does not indicate that a combination of these measures cannot be used to advantage.
  • obtaining a set of equations refers to the computer system automatically or semi-automatically, setting up the set of equations according to operational parameters of the district energy grid.
  • the computer system may take into account the size or number of users of the district energy grid, input data for one or more individual zones of the district energy grid, and/or any other relevant input parameters relating to the district energy grid, where the input data and/or other relevant input parameters may be obtained automatically and/or through user input.

Landscapes

  • Chemical & Material Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Combustion & Propulsion (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

L'invention concerne un procédé mis en œuvre par ordinateur pour la détermination d'une perte thermique dans un réseau d'énergie collectif, le réseau d'énergie collectif comprenant une ou plusieurs conduites d'alimentation et une pluralité de conduites de desserte pour le transport d'un fluide jusqu'à une pluralité de consommateurs, chaque conduite d'alimentation apportant le fluide dans une pluralité de conduites de desserte et chaque conduite de desserte reliant un consommateur à l'une desdites conduites d'alimentation, le procédé consistant à : recevoir des données de capteur mesurées en provenance d'un ou de plusieurs capteurs associés à une ou plusieurs conduites de desserte sélectionnées parmi la pluralité de conduites de desserte ; déterminer, à partir des données de capteur reçues, des températures de fluide observées et des écoulements de fluide observés respectifs du fluide transporté reçus à des instants respectifs par l'intermédiaire de conduites respectives desdites conduites de desserte sélectionnées ; obtenir un ensemble d'équations, l'ensemble d'équations comprenant une pluralité d'équations mettant en relation les températures de fluide observées et les écoulements de fluide observés avec des paramètres de perte thermique indicatifs de pertes thermiques de conduites respectives desdites conduites de desserte sélectionnées, calculer une solution de l'ensemble d'équations et calculer une valeur de résultat d'au moins un premier paramètre de perte thermique indicatif d'une perte thermique d'une première conduite de desserte desdites conduites de desserte sélectionnées à partir de la solution calculée de l'ensemble d'équations.
PCT/EP2023/062638 2022-05-12 2023-05-11 Procédé et système de détermination de perte de chaleur dans une conduite de desserte d'un système d'énergie collectif WO2023217979A1 (fr)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DKPA202270257 2022-05-12
DKPA202270257 2022-05-12

Publications (1)

Publication Number Publication Date
WO2023217979A1 true WO2023217979A1 (fr) 2023-11-16

Family

ID=86609903

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/EP2023/062638 WO2023217979A1 (fr) 2022-05-12 2023-05-11 Procédé et système de détermination de perte de chaleur dans une conduite de desserte d'un système d'énergie collectif

Country Status (1)

Country Link
WO (1) WO2023217979A1 (fr)

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3531368A1 (fr) 2018-02-21 2019-08-28 Kamstrup A/S Analyse de réseau de distribution public

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3531368A1 (fr) 2018-02-21 2019-08-28 Kamstrup A/S Analyse de réseau de distribution public

Similar Documents

Publication Publication Date Title
CN109791637B (zh) 用于布水系统的基于模型的控制的方法和装置
EP3531368B1 (fr) Analyse de réseau de distribution public
CN110168338A (zh) 用于基于模型的管道网络的泄漏检测的方法和装置
RU2428708C2 (ru) Способ и система для поверки измерительных приборов
US10508966B2 (en) Water flow analysis
JP6318053B2 (ja) 漏水分布推定装置
WO2018132138A1 (fr) Systèmes et procédés de modélisation hydraulique de sous-réseau
US20140149054A1 (en) Leak Detection Via a Stochastic Mass Balance
KR101105192B1 (ko) 관망해석에 의한 상수도 누수혐의구간 선정방법 및 그 기능을 탑재한 상수도 통합관리운영시스템
CN110274669B (zh) 一种大口径电子水表在线检定方法
Perdana et al. Water leak detection and shut-off system on water distribution pipe network using wireless sensor network
Moors et al. Automated leak localization performance without detailed demand distribution data
Farah et al. Smart water technology for leakage detection: feedback of large-scale experimentation
Li et al. Leakage localization using pressure sensors and spatial clustering in water distribution systems
CN114449373B (zh) 一种基于物联网无线传输的智慧水务系统
WO2023217979A1 (fr) Procédé et système de détermination de perte de chaleur dans une conduite de desserte d'un système d'énergie collectif
Agathokleous et al. Real-time monitoring of water distribution networks
Choi et al. Pressure measurements with valve-induced transient flow in water pipelines
CN207601643U (zh) 混凝土表面保护系统及环境指标管理系统
Arnold et al. Sensor-based cutoff strategy for border check–irrigated fields
WO2023217977A1 (fr) Procédé et système de détermination d'un paramètre de fonctionnement d'un réseau d'énergie urbain dans un système d'énergie urbain
EP4071454A1 (fr) Procédé de localisation d'une fuite dans un réseau d'alimentation en eau
CN108021157A (zh) 混凝土表面保护系统及方法
CN108763773B (zh) 直埋热水管网的热损计算方法及装置
Ayad et al. A Model-Based Approach for Leak Detection in Water Distribution Networks Based on Optimisation and GIS Applications

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 23727470

Country of ref document: EP

Kind code of ref document: A1