Flow determination
11041746 · 2021-06-22
Assignee
Inventors
Cpc classification
G01F15/02
PHYSICS
International classification
G01F1/684
PHYSICS
Abstract
A method and apparatus for making a flow determination with respect to a flow through a fluid conduit is described. The apparatus comprises a first temperature sensor arranged to generate a temperature signal indicative of the temperature of the outer surface of the fluid conduit, a second temperature sensor arranged to generate a temperature signal indicative of the ambient temperature outside of the fluid conduit; and a processor arranged to make the flow determination. The flow determination is made by determining the first and second temperatures at a first time; predicting a predicted first temperature at a second time; determining the first temperature at the second time; and comparing the predicted first temperature with the determined first temperature at the second time.
Claims
1. An apparatus for making a flow determination with respect to a flow through a fluid conduit, the fluid conduit comprising a wall having an inner surface surrounding a flow space through which the fluid flows and an outer surface, the apparatus comprising: a first temperature sensor arranged to be mounted on the outer surface of the fluid conduit and arranged to generate a first temperature signal indicative of a first temperature being the temperature of the outer surface; a second temperature sensor arranged to be positioned spaced apart from the first temperature sensor and to generate a second temperature signal indicative of a second temperature being the ambient temperature outside of the fluid conduit; and a processor having inputs for the first and second temperature signals; in which the processor has an output for the flow determination and is arranged to make the flow determination by: determining the first and second temperatures at a first time; using the first and second temperatures at the first time, predicting a predicted first temperature representative of the first temperature at a second time later than the first time; determining the first temperature at the second time; comparing the predicted first temperature with the determined first temperature at the second time; and making the flow determination dependent upon the comparison.
2. An apparatus according to claim 1, wherein the processor is arranged to compare the predicted first temperature with the determined first temperature by determining whether the determined first temperature at the second time is at least one of: closer to the second temperature than the predicted first temperature, further from the second temperature than the predicted first temperature, or equal to the predicted first temperature.
3. An apparatus according to claim 1, wherein the predicted first temperature represents a zero flow state.
4. An apparatus according to claim 1, wherein the processor is arranged to determine the magnitude of the difference between the predicted first temperature and the determined first temperature, and wherein the processor is arranged to determine the level of flow based on the determined magnitude of the difference between the predicted first temperature and the determined first temperature.
5. An apparatus according to claim 1, further comprising a memory configured to hold a plurality of expected convergence curves relating to fluid conduits of different sizes and materials.
6. An apparatus according to claim 5, wherein each expected convergence curve corresponds to a model which assumes no flow in the fluid conduit.
7. An apparatus according to claim 5, wherein each expected convergence curve corresponds to a model which assumes that the difference between the first and second temperatures will decay exponentially, or a numerical approximation to an exponential decay.
8. An apparatus according to claim 5, wherein the processor is arranged to determine the expected convergence curve and to store it in the memory by recording the temperature of the fluid conduit over a period of time when there is known to be no flow, but that starts at a point where the temperature of the fluid conduit and the ambient temperature is different.
9. An apparatus according to claim 1, wherein predicting a predicted first temperature representative of the first temperature at a second time later than the first time is performed under an assumption of a zero flow state, wherein the first temperature converges with the ambient temperature following a predictable curve.
10. An apparatus according to claim 1, wherein the second time is before stable convergence of the first and second temperatures.
11. An apparatus according to claim 1, wherein making the flow determination comprises categorizing the flow as one of usage, stable convergence, stable non-convergence, predicted convergence, predicted non-convergence or undetermined.
12. A method of making a flow determination with respect to a flow through a fluid conduit, the fluid conduit comprising a wall having an inner surface surrounding a flow space through which the fluid flows and an outer surface, the method comprising: measuring, at a first time, a first temperature being the temperature of the outer surface; measuring, at the first time, a second temperature being the ambient temperature outside of the fluid conduit; using the first and second temperatures at the first time, predicting a predicted first temperature representative of the first temperature at a second time later than the first time; measuring the first temperature at the second time; comparing the predicted first temperature with the measured first temperature at the second time; and making the flow determination dependent upon the comparison.
13. A method according to claim 12, wherein comparing the predicted first temperature with the determined first temperature at the second time comprises determining whether the determined first temperature at the second time is at least one of: closer to the second temperature than the predicted first temperature, further from the second temperature than the predicted first temperature, or equal to the predicted first temperature.
14. A method according to claim 12, further comprising determining the magnitude of the difference between the predicted first temperature and the determined first temperature, wherein the magnitude of the difference between the predicted first temperature and the determined first temperature is used to determine the level of flow.
15. A sensing apparatus, comprising: a first temperature sensor; a sensor head housing the first temperature sensor, the sensor head being arranged so as to hold the first temperature sensor against a fluid conduit, with the first temperature sensor being arranged so as to produce in use a first temperature signal indicative of a first temperature being that of the fluid conduit, wherein the first temperature sensor is provided with a sensor biasing member, which tends to bias the first temperature sensor into contact with the fluid conduit in use; a second temperature sensor, arranged to produce in use a second temperature signal indicative of a second temperature being the ambient temperature; in which the sensor head comprises a engagement portion arranged to engage the fluid conduit to retain the fluid conduit, in which the engagement portion has three contact positions circumferentially around the fluid conduit and does not otherwise contact the fluid conduit.
16. A sensing apparatus according to claim 15, wherein the engagement portion comprises a biasing member arranged to bias at least one of the contact positions into contact with the fluid conduit.
17. A sensing apparatus according to claim 15, having an air gap between the engagement portion and/or the sensor head on the one hand, and the first temperature sensor on the other hand.
18. A sensing apparatus according to claim 15, wherein the first temperature sensor is mounted in the sensor head using a coupling which can be repeatedly released and engaged.
19. A kit of parts comprising a sensing apparatus according to claim 18 and a plurality of interchangeable sensor heads.
20. A fluid conduit having a circumference fitted with the apparatus of claim 15, wherein the apparatus only contacts the fluid conduit at the contact positions.
Description
(1) There now follows, by way of example only, description of embodiments of the invention, described with reference to the accompanying drawings in which:
(2)
(3)
(4)
(5)
(6)
(7)
(8)
(9) A fresh water plumbing network for a domestic dwelling is shown schematically in
(10) In order to make a flow determination—typically to determine whether there is a leak from the plumbing network—a flow determination apparatus is used. This comprises multiple housings 1, 2 at different locations on the plumbing network, and a remote processor.
(11) The size of the divergence of the temperature of the water (and hence the pipe) on the one hand and the ambient temperature may be less when far into the plumbing network. As such, for the branches 16 deeper into the network, a heating/cooling apparatus 7 can be provided which selectively provides heating 9 or cooling 8 to the pipe 16 and hence the fluid, so as to increase the divergence in temperature when water flows. Typically, the heater 9 would be provided to heat the fluid, with smaller coolers 8 to cool preferentially the pipe 16 so that heat from the heater 9 does not propagate down the pipe rather than through the water.
(12) The housings 1, 2 are identical. A main housing 1 is provided on the main supply pipe 14, whereas an auxiliary housing 2 can be provided on each branch 15, 16. The housings 1, 2 are described in more detail using the example of the main housing 1 with reference to
(13) In these Figures, it can be seen that the housing 1 houses first 110 and second 120 temperature sensors (these may correspond to first temperature sensors 10/12 and second temperature sensors 11/13 of housings 1/2 in
(14) The sensor head also comprises a pair of slightly elastic arcuate jaws 141, 142. These can be forced over a suitably-sized pipe 4, but will retain the pipe thereafter.
(15) The sensor head 140 further comprises a pair of outriggers 143, 144. The outriggers 143, 144 may prevent the sensor head 140 from levering off the pipe when the sensor head 140 is pulled.
(16) The first temperature sensor 110, when installed in the sensor head 140 will protrude slightly from bore 151 so as to define a contact face 139 for the sensor head. The jaws can surround the pipe 4 (as shown in
(17) The first temperature sensor 110 is provided with a biasing spring 149, which biases it out of the bore 147, into contact with a pipe between the jaws 141, 142.
(18) An alternative sensor head 161 is shown in
(19) It can be seen that the installer of the sensor head can be provided with multiple different sensor heads, such as different sized versions of the sensor head 140 shown in
(20) In both cases, a housing 41 depending from the sensor head 140 houses the second temperature sensor 120. Because the pipe is thermally insulated from the sensor head, there is little danger of the temperature of the pipe 4 affecting the ambient temperature which is measured by the second temperature sensor.
(21) Each housing is also provided with a transmitter—such as a Bluetooth® Low Energy transmitter—which can carry out some processing and transmits data to the remote processor. Each housing is also provided with a power source (not shown), such as a battery, to power the transmitter and the temperature sensors.
(22) The data collected by the sensors shown in
(23) The apparatus relies on the fact that, if there is no flow in the pipe 4, then the temperature of the pipe—sensed by the first temperature sensor 110 will converge with the ambient temperature—sensed by the second temperature sensor 120 following a predictable curve.
(24) When there is a substantial flow, the temperature of the pipe 4 will typically diverge substantially from the ambient temperature. This is most notable in domestic plumbing networks the closer to the point of entry of the supply pipe 14 into the premises. This is because the temperature of the fluid flowing through the pipe 4—here, water—is likely to be different to the ambient temperature. In the domestic plumbing context, this is because pipes external to the dwelling are buried in the ground. In temperate climates such as the United Kingdom, it is likely that the water flowing into a dwelling will be significantly lower than ambient temperature and this explanation will be based on that assumption, although this embodiment will function well also with water significantly above ambient (for example, in an air-conditioned home in a hot climate).
(25) This means that, in the example of a temperate climate, a substantial flow will lead to a sudden drop in temperature of the fluid flowing through the pipe 4 and so a drop in the temperature of the pipe 4 itself.
(26) Where there is a low flow, the temperature of the fluid in the pipe 4 and so the pipe 4 itself will still move towards ambient temperature. We have appreciated that the curve with which the temperature moves towards the ambient temperature with time is different from that when there is no flow, and that this can be used to determine whether there is any flow and to estimate the level of that flow.
(27) This can be demonstrated by considering
(28) In this example, at time t.sub.1, the water through the dwelling was switched off; thus it was known that there was no flow. It can be seen that the trace 31 followed a particular curve between time t.sub.1 and a later time t.sub.2. This can be used to generate an expected convergence curve (or one can be calculated depending on the size and material of the pipe using standard fluid thermodynamic techniques). At time t.sub.3 the water supply was restored to the dwelling and a small amount of divergence is seen as there is some flow to repressurise the plumbing system.
(29) Subsequent to time t.sub.3, it can be seen that there is some convergence with the ambient temperature 30. However, at time t.sub.4, a toilet which had previously been disabled was reconnected, which had a leaking cistern. This had a leak of approximately 0.06 ml/second. This caused a substantial divergence from the ambient temperature until time t.sub.5 when a tap was used, causing a substantial flow and a sudden divergence from the ambient temperature until time t.sub.6.
(30) At this point, the pipe temperature 31 begins to converge once more with the ambient temperature 30, but it can be seen that the convergence is slower than the curve between times t.sub.1 and t.sub.2. This is indicative of a small flow. In this example, the flow turned out to be another toilet that had a then undiagnosed leak.
(31) As such, in order to determine whether there is flow, the apparatus can measure the first and second temperatures at a first time t.sub.0. At that point, it can predict what the first temperature will be at a later time t.sub.1. Typically, this could be using the predictive convergence curve, or a simple exponential heating curve. At that later time t.sub.1, the first temperature can be measured again, and a comparison of prediction and actual measurements taken. Given that the prediction indicates no flow, if the actual measurement is further from the ambient temperature (in this case, colder) than the prediction, then there is some flow, and the larger the difference between prediction and actual flow can indicate the size of the flow. This can be repeated at intervals as required.
(32)
(33) As such, it can be seen that a binary determination of whether there is flow can be made based upon a determination of whether the prediction is correct. An indication of the level of flow can be made by determining the difference between the actual and expected convergence; the larger the difference, the higher the flow.
(34) In one embodiment, in order to determine whether there is any flow, a series of categorisations can be made for each batch of data. The categories determined are: Usage Stable non-convergence (low rate usage) Stable convergence (zero flow) Predicted non-convergence (low rate usage) Predicted convergence (zero flow) Undetermined
(35) The current embodiment samples the first and second temperatures every 10 seconds, but finer or coarser sampling is possible.
(36) Though the water temperature itself is not measured directly, it can be inferred by looking at the lowest point reached by the pipe wall temperature over a recent period that includes activity. Water supply temperature is very stable so in a scenario where there has been activity in the last (say) 24 hours it is a very good assumption that the lowest pipe wall temperature seen in the last 24 hours is the water supply temperature. This is not a quantity that is reported to the user but it is used as explained below.
(37) Usage
(38) There are broadly two symptoms of usage Rapid fall in pipe temperature Pipe temperature remaining at or near water supply temperature
(39) If a single pipe temperature reading is lower than the previous by an amount greater than or equal to a threshold (in this embodiment 0.3° C.) then the status for that moment is defined as “Usage”.
(40) If over a set of readings (in this embodiment, 5) a drop of at least a further threshold (in this embodiment, 0.4°) is seen, that is usage.
(41) Any reading directly after a usage reading where the value is still falling or level is regarded as a continuation of the usage.
(42) If a single reading that would otherwise not be usage is sandwiched between two which are, this can optionally be regarded as usage.
(43) Any reading which is level or falling and is within a threshold, in this embodiment, 1 degree C., of the water supply temperature may be regarded as usage.
(44) Furthermore, it is possible to determine the span between the lowest temperature recorded in the last, say, 24 hours and then look at a temperature change of a given proportion of the span as indicating a usage event. The rate of change as a proportion of the span can be used to indicate the size of the usage event, either numerically or as categories (e.g. small/medium/large).
(45) Stable Convergence
(46) Stable convergence indicates that there is no flow in the pipe.
(47) In this embodiment, it is defined as follows: A rolling average of, in this embodiment, 5 temperature readings is taken. When the difference between the rolling average and the current pipe temperature is stable and is equal or less than a threshold, in this embodiment, 0.3 degrees C. this is “stable convergence”. “Stable” means that a number of readings (in this embodiment 30) the difference does not vary more than a margin (in this embodiment, 0.2° C.).
(48) So based on these values, if the difference between pipe temperature and rolling average of ambient temperature stayed between 0.1 and 0.3° C., or between 0 and 0.2° C., for 5 minutes, that would be regarded as stable convergence.
(49) Note that the pipe and air temperature might be tracking up or down in this period, but it is the difference that is considered.
(50) Stable Non-Convergence
(51) Stable non convergence indicates that there is some flow in the pipe.
(52) It is defined as follows: A rolling average of in this embodiment 5 temperature readings is taken. When the difference between the rolling average and the current pipe temperature is stable and is greater than the threshold (0.3 degrees C. in this embodiment), this is “stable non-convergence” “Stable” means that a number of readings (in this embodiment 30) the difference does not vary more than a margin (in this embodiment, 0.2° C.).
(53) So these would count as stable non convergence: 5 minutes where the pipe temperature and the rolling average ambient temperature differed between 0.4 and 0.6 degrees. 5 minutes where the pipe temperature and the rolling average ambient temperature differed between 10 and 10.2 degrees.
(54) There are potentially scenarios that qualify both as usage and non-convergence. The test for usage should be applied first. If a situation qualifies as usage we don't bother looking for other designations.
(55) Prediction
(56) The prediction categories are needed where, for example, in some properties convergence might take 3 hours but usage occurs every 2 hours and so convergence is never given a chance to complete. In these cases, prediction comes into play.
(57) Theory
(58) In the model used, the rate of change of temperature of the pipe 4 is modelled. In this model, the change of temperature in one time interval—the time over which the algorithm is used—is given by:
ΔT=HeatGain(T.sub.ambient−T.sub.pipe)−FlowGain(T.sub.pipe−T.sub.supply)
(59) where T.sub.ambient is the ambient temperature measured by the second temperature sensor, T.sub.pipe is the temperature of the pipe wall as measured by the first temperature sensor, T.sub.supply is the temperature of the water in the pipe at the supply (which can be determined as the lowest pipe wall temperature reached, as that is the temperature that the pipe wall will reach after sustained flow).
(60) If there is no flow at all, flow gain=0. Where there is full flow usage, flow gain is very large.
(61) Heat gain, on the other hand, does not vary over time within a given installation. Once the optimum value has been found it should vary very little over time and so the previous value is likely to be a good starting point in any optimisation.
(62) If flow gain=heat gain, we might expect to see the pipe temperature stabilise exactly halfway between the water supply temperature and the ambient temperature. We might call this a gain ratio of 1.
(63) In a typical situation where ambient ˜20° C. and pipe ˜10° C., if the pipe temperature were to stabilise at 0.4 degrees below ambient, that would mean that the ratio of heat gain to flow gain (the gain ratio) was 24.
(64) Put another way, our “temperature difference” criteria for drawing the line between actual convergence and actual non-convergence discussed above equates to a gain ratio of above 24 for predicted convergence.
(65) Given a sequence of data where the temperatures are starting to converging, but will not achieve convergence before something else happens, the idea with prediction is to see what values of heat gain and flow gain, would, if used in an appropriate model, produce a trace most like that seen in the actual data. If values can be found which create a model that closely mirrors the actual data, the two gain values may be close to the reality.
(66) Practice
(67) Throughout the data set, there are patches of data that cannot be categorised as being any of the other designations. Once an unbroken run of un-categorisable data (meaning it does not satisfy the criteria for usage, convergence or non-convergence) has been captured, which is at least a predetermined length (in this embodiment, 200 values) then a prediction is attempted. If there is a run of un-categorisable data but it is less than the length then it can be designated as “undetermined”.
(68) Multiple attempts to find a model that fits are attempted (an iterative optimisation). In each case, candidate values for FlowGain and HeatGain are used. The algorithm moves through 2 dimensional (HeatGain, FlowGain) space looking for the values that give the best fit.
(69) Each cycle through the optimisation goes as follows: Values for HeatGain and FlowGain are chosen, initially by using any previously determined value. A curve for “predicted pipe temperature” is generated alongside the actual pipe temperature curve, as follows The first value of the predicted pipe temperature curve is the same as the actual pipe temperature curve. Thereafter each subsequent value of the predicted curve is worked out as follows:
T.sub.pipe(predicted)(n)=T.sub.pipe(predicted)(n−1)+increment(n−1)
(70) Where “increment”=HeatGain×(T.sub.ambient−T.sub.pipe(predicted))−FlowGain×(T.sub.pipe(predicted)−T.sub.supplyTemp).
(71) Thus for each data point in the set being analysed, a predicted pipe temperature will be generated.
(72) By the end of this exercise there will be a set “predicted pipe temperature” readings corresponding to the “actual” pipe readings. If the two data sets correspond closely, the model is a good fit, which means the values of HeatGain and FlowGain are broadly representative of what is going on in reality. The method compares the two arrays of points with a “root mean square difference” assessment, looking at the mean deviation of theoretical from actual.
(73)
(74) The optimisation process in effect hunts around in (HeatGain, FlowGain) space until it runs out of allotted attempts of finds values for which the model fits the actual data closely enough.
(75) The current algorithm optimises using binary sections on one of the gains at a time. So we might pick a midrange value for HeatGain, and hold it constant while we optimise by varying FlowGain. Then we fix the FlowGain (with whatever value the last optimisation gave us) and optimises for HeatGain. This is repeated a maximum number of times.
(76) For each of the gains in turn, the binary section method is as follows: Start with min and max possible values of gain (a and b) based on the possible range then Pick the midpoint x:x=(a+b)/2 Pick a small value e and evaluate the model deviation at x−e and at x+e If the variance at (x−e)<variance at (x+e) then b=x+e Otherwise a=x−e Repeat for the chosen number of cycles
(77) Using the Output of the Prediction
(78) For a given block of data, the output of the prediction (when the optimisation has gone as far as it can) will be A flow gain A heat gain A measure of how well the model (the best case one that was found) fits the actual data.
(79) If the measure of model quality is as shown above (i.e. minimising the value of “root mean squared deviation from the actual data”) then the measure will be in units of a degree. A fit of 0.1° C. has proved possible on a typical data set when the region of data does represent something looking like convergence.
(80) Therefore, in this embodiment, the rule for interpreting a prediction is as follows if the best model we can find has a deviation greater than a threshold (in this embodiment 0.1° C.) then we consider the model not to be a good match. Each of the points in the chunk of data should be designated as “undetermined” If on the other hand the deviation is less than the threshold then we regard this as a match. If the gain ratio (HeatGain/FlowGain) is more than a threshold, say 25, then we regard that as “predicted convergence”. Otherwise it is “predicted non convergence”.
(81) Combining all of the Algorithm Data into Output
(82) By the time a complete data set has been analysed, every data point will be designated as either Usage Stable non-convergence (low rate usage) Stable convergence (zero flow) Predicted non-convergence Predicted convergence Undetermined
(83) Events
(84) It is then possible to list all of the “usage events”. A usage event is simply a continuous unbroken sequence of points designated as “usage”. It has a start time and a duration. For example, if there are 3 consecutive data points that is regarded as 30 seconds duration.
(85) An event may qualify as an alert if it meets certain criteria. The criteria are user determined. For the engineering application we will only offer a single alert setting, “usage exceeding a set duration”. Alerts generate notifications until acknowledged).
(86) Leaks
(87) The output with respect to leaks is, in this embodiment, based on the last 24 hours. If there has been actual stable convergence in the last 24 hours, then the property is regarded as leak free. If there has not been convergence but there has been stable non-convergence, then the property is regarded as having a leak. The magnitude of the non-convergence is a measure of the severity of the leak. If neither of the above is true, then the assessment must be based on any well matched predictions in the last 24 hrs. If there is a prediction of convergence, it takes precedence over a prediction of non-convergence.
(88) Thus, looking at
(89) It will be understood that the present invention has been described above purely by way of example, and modification of detail can be made within the scope of the invention.