AU2008201906B2 - Method for improving timekeeping and saving energy on long-haul trains - Google Patents

Description

- 1 AUSTRALIA PATENTS ACT 1990 COMPLETE SPECIFICATION FOR A STANDARD PATENT ORIGINAL Name of Applicant/s: TMG Rail Technology Pty Limited Actual Inventor/s: Philip George Howlett and Peter John Pudney Address for Service is: SHELSTON IP 60 Margaret Street Telephone No: (02) 9777 1111 SYDNEY NSW 2000 Facsimile No. (02) 9241 4666 CCN: 3710000352 Attorney Code: SW Invention Title: SYSTEM FOR IMPROVING TIMEKEEPING AND SAVING ENERGY ON LONG-HAUL TRAINS Details of Original Application No. 2003229097 dated 20 May 2003 The following statement is a full description of this invention, including the best method of performing it known to me/us: File: 34761AUP02 - la TITLE: METHOD FOR IMPROVING TIMEKEEPING AND SAVING ENERGY ON LONG-HAUL TRAINS FIELD OF THE INVENTION This invention relates to a method for the operation of trains on a rail network, 5 and has particular application in the context of long-haul rail networks. BACKGROUND OF THE INVENTION The energy costs for railways are significant. By driving efficiently, these costs can be significantly reduced. There are five main principles of efficient driving: 10 1. Aim to arrive on time. If you arrive early you have already wasted energy; if you arrive late you will waste energy making up the lost time. 2. Calculate your required average speed. On long journeys, simply dividing the distance remaining by the time remaining will give you an approximate holding speed. Recalculate during the journey to make sure you are still on target. 15 3. Aim to drive at a constant speed. Speed fluctuations waste energy. The most efficient way to drive is to aim for a constant speed. 4. Avoid braking at high speeds. Braking at high speeds is inefficient. Instead, coast to reduce your speed before declines and speed limits. 5. Anticipate hills. If the train is going to slow down on a steep incline, 20 increase your speed before the incline so that the average speed on the incline does not drop too far below the hold speed. For steep declines, coast before the decline so that the average speed does not rise too far above the hold speed. Avoid braking. A train journey can be divided into segments between "targets", that is, locations on the route where the time and speed are specified. There are many driving 25 strategies that may be used to operate a train between one target and the next. One strategy is a "speed-holding" strategy, where a constant speed is maintained, except where prevented by speed limits and steep gradients. In practice, of course, speed limits and steep gradients can disrupt a significant part of a journey. If an efficient journey for a given holding speed Vcan be determined then Vcan be adjusted to find 30 the efficient journey that satisfies the journey time constraint; if the time taken is too -2 long then V is too low. In determining an appropriate holding speed it is possible to generate points on a cost-time curve for the journey. Using this methodology a journey with holding speed V can be constructed as follows: 5 1. Ignoring speed limits and the initial and final speeds, construct a speed holding journey with holding speed V. The speed of the train will vary with steep gradients. 2. Adjust the speed-holding journey to satisfy the speed limits. 3. Construct initial and final phases to satisfy the initial and final speed 10 constraints. However, using this methodology may not result in the most energy-efficient journey. It is therefore an object of the present invention to provide a method for operating trains which overcomes or ameliorates at least one of the disadvantages of 15 the prior art, or at least provides a useful alternative. SUMMARY OF THE INVENTION A first aspect of the present invention provides a method of monitoring the progress of a train on a rail network and providing driving advice in real time to an operator of said train, said method comprising: 20 (i) estimating or determining parameters of said train; (ii) determining, by an optimal control algorithm employing an adjoint variable, an optimal journey profile for a journey from said train's current location to a target location that results in said train arriving at said target location as close as possible to a desired time and with minimum energy usage; said optimal journey 25 profile including a speed profile for the train, sequence of discrete control modes for said train, and associated switching points between said control modes; said optimal journey profile being determined by solving a system of differential equations for said speed profile of the train and for the value of said adjoint variable, said control modes being determined from the value of said adjoint 30 variable, such that said sequence of control modes is determined as said speed profile is calculated; -3 (iii) monitoring the current state of said train as it progresses to said target location; and (iv) generating said driving advice for the train operator by comparing the current state of the train to a corresponding state on said optimal journey profile and 5 displaying said advice for the train operator that will keep the train close to said optimal journey profile. A second aspect of the present invention provides a method of monitoring the progress of a train on a rail network and providing information on the progress of the train in real time to an operator of said train, said method comprising: 10 (i) estimating or determining parameters of said train; (ii) determining, by an optimal control algorithm employing an adjoint variable, an optimal journey profile for a journey from said train's current location to a target location that results in said train arriving at said target location as close as possible to a desired time and with minimum energy usage; said optimal journey 15 profile including a speed profile for the train, sequence of discrete control modes for said train, and associated switching points between said control modes; said optimal journey profile being determined by solving a system of differential equations for said speed profile of the train and for the value of said adjoint variable, said control modes being determined from the value of said adjoint 20 variable, such that said sequence of control modes is determined as said speed profile is calculated; (iii) monitoring the current state of said train as it progresses to said target location; and (iv) generating said information for the train operator by comparing the 25 current state of the train to a corresponding state on said optimal journey profile and displaying said information for the train operator to assist in keeping the train close to said optimal journey profile. A third aspect of the present invention provides a method of controlling the progress of a train on a rail network, said method comprising: 30 (i) estimating or determining parameters of said train; (ii) determining, by an optimal control algorithm employing an adjoint variable, an optimal journey profile for a journey from said train's current location to - 3a a target location that results in said train arriving at said target location as close as possible to a desired time and with minimum energy usage; said optimal journey profile including a speed profile for the train, sequence of discrete control modes for said train, and associated switching points between said control modes; 5 said optimal journey profile being determined by solving a system of differential equations for said speed profile of the train and for the value of said adjoint variable, said control modes being determined from the value of said adjoint variable, such that said sequence of control modes is determined as said speed profile is calculated; 10 (iii) monitoring the current state of said train as it progresses to said target location; and (iv) comparing the current state of the train to a corresponding state on said optimal journey profile and then controlling said train to keep the train close to said optimal journey profile. 15 Preferably, said adjoint variable evolves according to a differential equation along with the position and speed of the train. Preferably, the value of the adjoint variable is calculated directly from the speed of the train. Preferably, a numerical method is used to solve a system of differential 20 equations for said speed profile of the train and for the value of said adjoint variable. Preferably, steps (i) to (iv) are performed as required so that said driving advice automatically adjusts to compensate for any operational disturbances encountered by the train. Preferably, said discrete control modes for said train include drive, hold, coast 25 and brake modes. Preferably, said parameters include train mass and mass distribution. Preferably, said parameters further include maximum tractive efforts and maximum braking effort as functions of speed. Preferably, said parameters further include coefficient(s) of rolling resistance. 30 Preferably, said driving advice is generated and displayed by a computer located on the train.

-3b Preferably, step (iii) involves processing data from a GPS unit and train controls to determine the location and speed of the train. Preferably, said optimal journey profile specifies the time, speed and control at each location between the current train location and the next target location on the 5 network.

-4 Preferably, said optimal journey profile is precomputed. Preferably, a plurality optimal journey profiles corresponding to different journey times are calculated and the profile that has an arrival time at the target location closest to the desired arrival time is selected. 5 Preferably, steps (i) to (iv) are performed as required so that said driving advice automatically adjusts to compensate for any operational disturbances encountered by the train. In a particularly preferred embodiment of the invention the optimal journey profile comprises driving in a hold mode (i.e. at constant speed), calculated by the 10 Pontryagin Principle of optimal control, wherever possible and where it is not possible changing as quickly as is safely possible at exactly the right location, calculated by the Pontryagin Principle of optimal control, to drive (i.e. full power), hold, coast or brake modes as necessary. Unless the context clearly requires otherwise, throughout the description and 15 the claims, the words "comprise", "comprising", and the like are to be construed in an inclusive sense as opposed to an exclusive or exhaustive sense; that is to say, in the sense of "including, but not limited to". BRIEF DESCRIPTION OF THE DRAWINGS The invention will now be described in further detail, by way of example only, 20 with reference to the accompanying drawings in which: Figure 1 shows a block diagram of the system according to a preferred embodiment of the present invention, illustrating the main data flows between various elements of the system; Figure 2 illustrates an optimal speed profile for a train over a fictitious section 25 of track; Figure 3 illustrates an optimal speed profile for a train over another fictitious section of track; Figure 4 illustrates an optimal journey for a coal train; Figure 5 shows the processing of precomputed speed profiles; and 30 Figure 6 illustrates a preferred embodiment of the system display which provides the train operator with driving advice in real time.

-5 DESCRIPTION OF PREFERRED EMBODIMENT The present invention, in one preferred form, provides a fully automatic system that monitors the progress of a train on a long-haul network, calculates efficient control profiles for the train, and displays driving advice to the train crew. In 5 a further preferred embodiment the system works in conjunction with a dynamic rescheduling tool that coordinates interactions between various trains operating on the network. The system assists the crew of a long-haul train by calculating and providing driving advice that assists to keep the train on time and reduce the energy used by the 10 train. The system performs four main tasks: (i) state estimation: monitors the progress of a journey to determine the current location and speed of the train; (ii) train parameter estimation: estimates some parameters of a train performance model; 15 (iii) journey optimisation: calculates or selects an energy-efficient driving strategy that will get the train to the next key location as close as possible to the desired time; and (iv) advice generation: generates and provides driving advice for the driver. These tasks are performed continually so that the driving advice automatically 20 adjusts to compensate for any operational disturbances encountered by the train. The system includes: e data communications between on-board units and a central control system; e automatic estimation of train performance parameters; e automatic re-optimisation of optimal journey profiles; 25 e interaction with a manual or automatic train rescheduling system; e ergonomic driver interfaces. Each of these four aspects of the methodology and system will now be discussed in further detail: State estimation 30 The station estimation task processes observations from a GPS unit and the train controls to determine the location and speed of the train and the current control setting.

-6 Location is the position of the train on a given route, and is used to look up track gradient, curvature and speed limits. The state estimation task uses absolute and relative position data to determine the location of the train. Control setting is required for train parameter estimation, and for estimating 5 the energy use of the train if direct measurement of energy use is not available. Train parameter estimation The train parameter estimation task estimates parameters of a train performance model from the sequence of observed journey states. The train model used by the in-cab system has the following train parameters: 10 e train mass and mass distribution; e maximum tractive effort and maximum braking effort as functions of speed; and e coefficients of rolling resistance. Any of these parameters that are not known with sufficient accuracy before the 15 journey commences must be estimated during the journey. The unknown parameters can be estimated using a Kalman filter. If mass is to be estimated, the mass distribution is assumed to be uniform. If tractive effort is to be estimated it is assumed to take the form F(v)= ' ! -v >vo 20 where P is the maximum power of the train and vo is the speed below which maximum tractive effort is assumed to be constant. In the simplest implementation, all train model parameters are known in advance and parameter estimation is not required. Journey Optimisation 25 The optimal journey profile between a given journey state and a target journey state is found by solving a set of differential equations for the motion of the train and an additional differential equation that determines the optimal control. The optimal journey profile specifies the time, speed and control at each location of the track between the current train location and the next target. 30 Journey profiles can be precomputed or else calculated during the journey. If precomputed, several different journeys corresponding to different journey times are -7 used on the train and the journey optimisation task then simply selects the precomputed profile that has the arrival time at the target closest to the desired arrival time. If we use distance travelled, x, as the independent variable then the journey 5 trajectory is described by the state equations -= l/v (1) dx dv u--R(v)+G(x) (2) dx mv dJ (3) dx where t is elapsed time, v is the speed of the train, J is energy use, u is the controlled driving or braking force, R(v) is the resistive force on the train at speed v and G(x) is force on the train due to track gradient and curvature at 10 location x, and m is the mass of the train. We assume that R and the derivative R 'are both increasing functions. This model is based on simple physics. It does not model the complexities of traction motors, braking systems, in-train forces or wheel-rail interations. Nor does it need to; in practice, the driving advice derived from this simple model is both realistic 15 and effective. The state equations describe the motion of a point mass. In practice the length of a long-haul train can be significant. However, a long train can be treated as a point mass by transforming the track force function. Suppose the train has length L and that the density of the train at distance I from the front of the train is p(l). If we define 20 G(x) = Jp(I)G(x - )dl where G is the real track force then the motion of a point mass train on a track with track force G is equivalent to the motion of the long train on the real track. The force u is controlled by the driver, and satisfies the constraints FB(v) u FD(v) where FD(v) > 0 is the maximum drive force that can be achieved at speed v and 25 FB(v) > 0 is the maximum braking force that can be achieved at speed v.

-8 For most train journeys the speed of the train is constrained by speed limits that depend on location, and so the optimal journey must satisfy the constraint v VL(X). The optimal control is founded by forming the Hamiltonian function H = I +;r2u - R(v) +G(x) +7 5 HV =3{1 +x,u+ 17U-] my -a,[F, (v)-u] - a,[u- FI)(v)]-a,[v -V,(x)] where ;, are multipliers associated with the state equations and a, are Lagrange multipliers associated with the control and speed constraints. The complementary slackness conditions are a,[F,(v) - u] = aJu - Fl(v)] = a[v - V,(x)] = 0 10 There are three adjoint equations. The first and third adjoint equations are -l =0 and 3 -0 dx dx Ifwe let , 3 = -l and p =' my then the second adjoint equation can be written as + pR'(v)+ a, + (1- ,u)F') (v)] u = Fl(v) 15 - <["+ pR'(v) +a,] F,(v) <u < FI)(v) (4) clx + pR'(v) + a,, + (77 - p)F'H (v)] u = FB(v) This equation is found by substituting each of the three control conditions into the Hamiltonian and then differentiating. The Lagrange multiplier a is zero when the train is travelling at a speed less than the speed limit. The optimal control maximises the Hamiltonian, and so the optimal control 20 depends on the value of the adjoint variable u. An optimal strategy has five possible control modes: drive I < p -> maximum drive force u = FD(v) hold p = 1 => speed hold with 0 u : FD(v) coast 7R < P < 1 --> coast with u= 0 -9 regen p = 7R = speed hold with FB(v) < u < 0 brake P < ?7 -> brake with u = FB(v) The hold mode is singular. For this driving mode to be maintained on a non trivial interval requires dp/dx = 0. If we are not constrained by a speed limit then we have v 2 R'(v) = -ir, 5 But 7 I is a constant and the graph y = v 2 R '(v) is strictly increasing, so there is a unique hold speed V satisfying this equation. Maintaining a speed limit also requires p = 1. When a speed limit is encountered the adjoint variable p jumps to p = 1 and at the same time the Lagrange multiplier a, jumps from zero to a positive value. 10 On a track with sufficiently small gradients and no speed limits the optimal trajectory is mainly speed holding at speed V. On most tracks, however, the track gradients disrupt this simple strategy. Track intervals can be divided into four speed dependent classes: (i) steep incline: if the maximum drive force is not sufficient to maintain 15 the desired speed; (ii) not steep: if the desired speed can be maintained using a non-negative drive force; (iii) steep decline: if braking is required to maintain the desired speed; and (iv) nasty decline: if even maximum brake force is insufficient to maintain 20 the desired speed. The optimal strategy anticipates steep gradients by speeding up before a steep incline and slowing down before a steep decline. An optimal trajectory with a given hold speed V can be found by setting ir, = VR '(V) 25 and then solving the differential equations (1) and (2) while using (4) and the optimal control modes to determine the control. These differential equations are solved using a numerical method such as a Runge-Kutta method. In practice, however, the adjoint equation is unstable. To overcome this difficulty we instead search for a pair of adjacent adjoint trajectories that are lower and upper bounds for the true adjoint - 10 trajectory. The lower and upper bounds start close together, but the adjoint values eventually diverge. This does not matter while they are both indicating the same control mode, but as soon as one of the bounds indicates a control change we re search at that location to find new adjacent bounds that extend the journey. 5 The optimal journey trajectory can be constructed in this way as a sequence of trajectory segments between speed-holding phases, where speed holding can occur at the hold speed V or at a speed limit. There are two ways a non-holding optimal trajectory segment can start: I. Drive or coast with (xo, vo) known and po unknown. This occurs at the 10 beginning of the journey or at the end of a low speed limit. Calculating an initial upper bound for p is not usually possible, so instead we search for the location of the next control change. 2. Drive or coast with xo unknown but bounded, vo known and po = 1. This may occur if we are holding at the hold speed or at a speed limit. The 15 lower bound for xO is the start of the hold phase. The upper bound for xo depends on whether we are holding at the hold speed V or at a speed limit. If we are holding at the hold speed V then the upper bound for xo is the next location where either the track becomes steep or else the speed limit drops below V. If we are holding at a speed limit VL then the upper bound for xo is 20 the next location where either the track becomes steep uphill or else the speed limit drops. If a steep decline is encountered during a speed limit phase then the brakes must be partially applied to hold the train at the speed limit. There are three ways a non-holding optimal trajectory segment can finish: 1. At the end of the journey, with the correct speed. 25 2. At the hold speed with v = V, p = I and the gradient not steep. The next trajectory segment will have start type 1. 3. At a speed limit with v = Vi. The next trajectory segment will have start type 2 with control coast, or else start type I with control drive. Using these conditions, it is possible to construct a complete journey profile to 30 the next target. This journey profile will be optimal for the resulting arrival time at the target. If the resulting arrival time is beyond the desired arrival time then another journey profile, with a higher hold speed, is calculated; if the arrival time at the target - 11 is prior to the desired arrival time then another journey profile is calculated, this time with a lower hold speed. A numerical technique such as Brent's method can be used to find the hold speed that gives the desired arrival time. As a result, the optimal journey profile comprises driving in a hold mode (i.e. 5 at constant speed), calculated by the Pontryagin Principle of optimal control, wherever possible and where it is not possible changing as quickly as is safely possible at exactly the right location, again calculated by the Pontryagin Principle of optimal control, to drive (i.e. full power), coast or brake modes as necessary. Advice Generation 10 The advice generation task compares the current state of the train to the corresponding state on the optimal journey profile and then generates and displays advice for the train operator that will keep the train close to the optimal profile. Brake advice is given if braking is required to avoid exceeding a speed limit or a speed on the journey profile that has braking as the optimal control. 15 Coast advice is given if: * the speed of the train is significantly higher than the speed indicated by the optimal journey profile, or e the speed of the train is near or above the speed indicated by the optimal journey profile and the optimal control is coast. 20 Hold advice is given if the speed of the train is near or above a holding speed indicated by the optimal journey profile. The speed to be held will be either a speed limit or the journey holding speed. Power advice is given if none of the other driving modes are appropriate. These decisions can be made without considering time because the optimal 25 speed profile is automatically adjusted by the journey optimisation task to keep the train on time. For each type of trip, the optimisation software is used to calculate optimal speed profiles for six difference total journey times. Each profile is designed to minimise fuel consumption for the given journey time. As the time allowed for the 30 journey decreases the minimum possible fuel consumption increases.

- 12 During the journey the system uses a GPS unit to determine the position of the train. Given the speed and position of the train and the time remaining until the train is due at the next key location, the system selects the most appropriate of the pre computed profiles. Advice is generated to keep the train as close as possible to the 5 selected profile. The crew will enter necessary information such as the arrival time at the next key location. The advice given to the driver will be one of: u Drive: drive using maximum power, subject to safety and train handling constraints; o Hold: vary the power to hold the indicated speed; or 10 o Coast: set the power to zero subject to safety and train handling constraints. Note that the driver is responsible for braking. The system is able to work with pre-computed profiles because, in practice, if the control is changed too early or too late, switching between the difference pre 15 computed profiles will automatically adjust future control changes to compensate. Energy savings can be achieved simply by demonstrating efficient control techniques to the train operator. Effective techniques can either be demonstrated on board or by using simulations. However, because of the relationship between fuel consumption and journey time some form of on-board advice system is required to 20 achieve the best possible fuel consumption, and is the reason why coasting boards by the side of the track do not work. For example, if a train is running slowly and behind schedule because of a head wind, and the driver coasts at the usual location, the train will end up even further behind schedule. Of course, drivers will take train performance into account, 25 but it is difficult for them to keep track of time and predict the effect their control decisions will have on the final arrival time. The system of the present invention achieves significant fuel savings without increasing running times because the system is an adaptive system based on optimal control theory. 30 The system can adjust the driving strategy using the actual observed train performance. All systems that rely on pre-computed profiles must take into account the current state of the train with regard to location, time and speed. Any system of - 13 non-adaptive control will give unreliable advice when the train is not in the right place at the right time doing the right speed. Whilst non-adaptive systems could possibly be used on Metropolitan railways with fixed timetables and identical trains or on tightly controlled networks with unit trains carrying consistent loads using dedicated track, 5 they are not suitable on networks where the trains and timetables vary from day to day. EXAMPLE In the following discussion of an example of the invention, the following notation is used: 10 Train m train mass (kg) FD(V) maximum drive force at speed v (N) FB(v) minimum brake force at speed v (N) R(v) resistance force at speed v (N) 15 '7R regenerative brake efficiency Route The length and mass distribution of a train can be used with a simple averaging procedure to transform the track gradients and speed limits so that the motion of a point mass train on the transformed track corresponds to the motion of the real train 20 on the real track. G(x) effective force due to gradient at distance X (N) h(x) effective elevation of the track at x (m) V(x) effective speed limit at x (ms-1) State variables 25 X distance along the route (m) t(x) time taken to reach distance x (s) v(x) speed at distance X (ms-I) J(x) energy cost at distance x (J) Control and adjoint variable 30 u applied drive force 0 s u FD(v) or brake force F(v) U < 0 (N) P an adjoint variable that determines the optimal control switching points - 14 Steep gradients and speed limits mean that travelling at a constant speed for the entire journey is usually not possible. To find the optimal control for real journeys we use Pontryagin's principle, a standard technique of optimal control theory. The method is described for trains with discrete control in the book by Howlett and 5 Pudney (1995), and for continuous control by Howlett and Khmelnitsky. The continuous control model is easier to work with, and the results from the two models are practically identical. The optimal control at any stage of the journey depends on the value of an adjoint variable " , which evolves as the journey progresses. There are five control modes in an optimal journey: drive 1<p > u=FD(v) hold p=1 =' OusF,(v) coast R= 0 regen P 1 7R => F(v) u 0 brake P<7R > u=F1(v) 10 By analysing the equations for " we can show that the control mode with l 1 corresponds to speed holding. We can also show that during any one optimal journey, speed holding must always occur at the same speed, V. W > V. The holding speed V and the regen speed W are related by the simple 15 formula qRW2R'(W) = V 2 R'(V). If regeneration is perfectly efficient then the regen speed is the same as the hold speed, and the coast mode never occurs. If the train does not have regenerative braking then the regen mode does not occur. 20 Using the same type of analysis we can show that the control mode with P =7R requires the use of regenerative braking to maintain a constant speed For a given hold speed V we can divide the track into four classes: e steep inclines, where maximum drive force is not sufficient to hold speed V; 25 e not steep, where a proportion of the maximum drive force is sufficient to hold speed V; e steep declines, where braking is required to hold speed V; and e nasty declines, where full brakes are not enough to hold speed V.

- 15 We will assume that there are no nasty declines, nor any inclines so steep that the train can not get up them even at low speed. The key to handling steep grades is to anticipate the grade. For steep inclines, the speed of the train should be increased before the start of the incline; for seep declines, speed should be reduced before the 5 start of the decline. Figure 2 shows an optimal journey segment on a fictitious section of track. The holding speed is 70km/h. The steep sections are each 1% grades. The optimal journey has the train coasting 2km before the start of the decline, and driving 500m before the start of the incline. The grey curve shows the adjoint variable used to determine the optimal control; it has been scaled and shifted to make it easier to see. 10 For both the drive and the coast phases the adjoint variable starts and finishes at pt = 1. Where steep grades are close together the correct switching sequence and switching points are more difficult to find, but they can be calculated using the adjoint equation. In Figure 3 the steep sections are once again 1% grades. The control is switched from power to coast as the adjoint variable R passes through g = 1, before the 15 top of the hill. The same principle can be used to find an optimal speed profile for more complex journeys. Figure 4 shows an optimal journey for a coal train. The hold speed is 70km/h. The elevation profile has been smoothed to compensate for the length and mass distribution of the train. 20 This is a particularly difficult journey; there is only one short period of speed holding, indicated by the dark shading at 220km. The lighter shading indicates periods of coasting. The dark shading at the end of the journey indicates braking. On long journeys the adjoint variable can be difficult to calculate. The light curves show lower and upper bounds for the adjoint variable. We have to search for a 25 more accurate value whenever the bounds become too far apart, or whenever one bound indicates a control change but the other does not. The method used to calculate an optimal journey is easily extended to handle speed limits (Pudney & Howlett, 1994; Howlett & Pudney, 1995; Cheng et al, 1999; Khmelntisky). Whenever the speed profile meets a speed limit there is no choice but 30 to apply partial braking to hold the speed of the train at the speed limit. At the point where the speed limit is encountered the value of the adjoint variable jumps by an amount that can be calculated. The optimal journey can be found as before, using the - 16 adjoint variable to determine the control and calculating the adjoint jump each time a speed limit is encountered. To find the optimal strategy for a given journey time we need to find the appropriate hold speed. Simply dividing the journey time by the journey distance 5 gives an initial guess. In most cases this guess will be an underestimate of the holding speed required; speed limits, gradients and the initial and final phases of a journey tend to reduce the actual average speed. The time taken for an optimal journey with hold speed V decreases as V increases. We simply use a numerical search technique to find the hold speed that 10 gives the correct journey time. As a by-product we generate a sequence of points (T, J) that describe the energy cost J of an optimal journey that takes time T. These points describe a cost-time curve that can be used for calculating timetables that take into account energy costs. It may appear that the speed-holding strategy for long-haul trains is different to 15 the drive-coast-brake strategy for suburban trains, but this is not so. On suburban journeys, the hold speed required to achieve the timetable on short journey sections is usually greater than the maximum speed that can be achieved before coasting and braking are required. The suburban drive-coast-brake strategy is simply a subset of the speed holding strategy used on longer journeys. 20 The invention is designed to work on a train with optimisation working as a background task continually updating the optimal speed profile from the current state of the journey to the next target. Advice is provided from the result of comparing the current state to the optimal journey and generating appropriate control advice. 25 Figure 5 shows the processing of precomputed speed profiles, and Figure 6 shows a typical advice task provided to a driver. Figure 6 shows a preferred embodiment of the driver display 10 for providing real time driving advice to the driver of the train. The target location 12 (in this example, "Crystal Brook") is selected by the 30 driver and represents the next destination the train must reach by a certain time. The estimated time of arrival (ETA) 14 is calculated by the system and represents the predicted time the train will reach the target location 12 based on the current location - 17 of the train, the distance to the target location, and the selected journey profile (i.e. driving strategy). If the calculated ETA does not satisfy the driver's requirements (i.e. by being too early or too late) the driver can select a "faster" or "slower" journey profile from a series of journey profiles. These profiles may be selected by the driver 5 from a graduated scale 16. In the preferred embodiment depicted the driver has a choice of seven (7) different journey profiles. As may be appreciated, the slower the journey profile the less fuel is used, whilst the faster the journey profile the less fuel efficient the journey. Line 18 on the display illustrates the vertical profile of the section of track on 10 the display, whilst line 20 depicts the track curvature, or horizontal profile, of the section of track on display. The line 22 represents the train, with the vertical line 24 denoting the location of the front of the train and vertical line 26 denoting the location of the rear of the train. The train 22 progresses from left to right on the display. 15 Line 28 and associated numbers 30 indicate the speed limits (in km/hr) in various zones of the section of track on display. In this example, the speed limit over the first zone is 75km/hr, then the speed limit reduces to 50km/hr, then increases to 60km/hr, then reduces to 50km/hr, and finally increases to 55km/hr. The coloured line 32 indicates the recommended driving profile for the train 20 over the various zones of track. The colour of line 32 at any point denotes the control mode the driver is required to use at that point on the track (i.e. brakes, coast, or power). In the display depicted, red represents "brake" mode, white represents "coast" mode, and green represents "power" mode. The shade of the colour varies to indicate the degree of braking or power required. The darker the shade of colour, the greater 25 the degree of braking or power required at the particular point on the track. This is particularly useful when the control mode is 'hold' which is, by nature, somewhere between full power and coast modes. Indicator 34 provides a visual indication to the driver as to how the train is progressing against the recommended speed profile. The indicator comprises a pair of 30 spaced apart arrows which move horizontally across the display as the train progresses and vertically to indicate how the train is progressing against the recommended speed profile. Ideally, the pair of arrows will span the recommended speed profile. If the - 18 train is travelling too slowly the arrows will fall below the line 32, whilst if the train is travelling too quickly the arrows will lie above the line 32. In the example shown, the arrows lie slightly below the line, indicating that the train is travelling slightly slower then recommended. 5 Advantageously, by using colour to advise the driver of the control mode, and colour shades to advise the driver as to the intensity of power or braking required, it is possible to move beyond a prescriptive and somewhat inefficient "power-hold-coast" form of driving advice and provide more intuitive driving advice to the driver. Advantageously, the present invention at least in the preferred form provides 10 one or more of the following benefits: e efficient driving strategies which can reduce energy costs whilsting improving time keeping and network performance. e improved on-time running and shorter waits at crossing loops; e reduced air braking, lower brake wear, reduced wear on traction motors, 15 extended service life, lower maintenance costs; e improved consistency between drivers; e accelerated driver training. Although the invention has been described with reference to specific examples, it will be appreciated by those skilled in the art that the invention may be 20 embodied in many other forms.

Claims (44)

1. A method of monitoring the progress of a train on a rail network and providing driving advice in real time to an operator of said train, said method comprising: (i) estimating or determining parameters of said train; 5 (ii) determining, by an optimal control algorithm employing an adjoint variable, an optimal journey profile for a journey from said train's current location to a target location that results in said train arriving at said target location as close as possible to a desired time and with minimum energy usage; said optimal journey profile including a speed profile for the train, sequence of discrete control modes for 10 said train, and associated switching points between said control modes; said optimal journey profile being determined by solving a system of differential equations for said speed profile of the train and for the value of said adjoint variable, said control modes being determined from the value of said adjoint variable, such that said sequence of control modes is determined as said speed profile 15 is calculated; (iii) monitoring the current state of said train as it progresses to said target location; and (iv) generating said driving advice for the train operator by comparing the current state of the train to a corresponding state on said optimal journey profile and 20 displaying said advice for the train operator that will keep the train close to said optimal journey profile.

2. The method as claimed in claim I wherein said discrete control modes for said train include drive, hold, coast and brake modes.

3. The method of monitoring the progress of a train on a rail network as claimed 25 in claims I or 2 wherein said adjoint variable evolves according to a differential equation along with the position and speed of the train.

4. The method of monitoring the progress of a train on a rail network as claimed in claims I or 2 wherein the value of the adjoint variable is calculated directly from speed of the train. 30

5. The method of monitoring the progress of a train on a rail network as claimed in any one of claims I to 4 wherein a numerical method is used to solve said system of -20 differential equations for said speed profile of the train and for the value of said adjoint variable.

6. The method of monitoring the progress of a train on a rail network as claimed in any one of claims 1 to 5 wherein steps (i) to (iv) are performed as required so that 5 said driving advice automatically adjusts to compensate for any operational disturbances encountered by the train.

7. The method of monitoring the progress of a train on a rail network as claimed in any one of claims I to 6 wherein said parameters include train mass and mass distribution. 10

8. The method of monitoring the progress of a train on a rail network as claimed in claim 7 wherein said parameters further include maximum tractive efforts and maximum braking effort as functions of speed.

9. The method of monitoring the progress of a train on a rail network as claimed in claims 7 and 8 wherein said parameters further include coefficient(s) of rolling 15 resistance.

10. The method of monitoring the progress of a train on a rail network as claimed in any one of claims I to 9 wherein said driving advice is generated and displayed by a computer located on the train.

I . The method of monitoring the progress of a train on a rail network as claimed 20 in any one of claims I to 10 wherein step (iii) involves processing data from a GPS unit and train controls to determine the location and speed of the train.

12. The method of monitoring the progress of a train on a rail network as claimed in any one of claims I to 11 wherein said optimal journey profile specifies the time, speed and control at each location between the current train location and the next 25 target location on the network.

13. The method of monitoring the progress of a train on a rail network as claimed in any one of claims 1 to 12 wherein said optimal journey profile is precomputed.

14. The method of monitoring the progress of a train on a rail network as claimed in any one of claims I to 12 wherein a plurality optimal journey profiles 30 corresponding to different journey times are calculated and the profile that has an arrival time at the target location closest to the desired arrival time is selected. - 21

15. A method of monitoring the progress of a train on a rail network and providing information on the progress of the train in real time to an operator of said train, said method comprising: (i) estimating or determining parameters of said train; 5 (ii) determining, by an optimal control algorithm employing an adjoint variable, an optimal journey profile for a journey from said train's current location to a target location that results in said train arriving at said target location as close as possible to a desired time and with minimum energy usage; said optimal journey profile including a speed profile for the train, sequence of discrete control modes for 10 said train, and associated switching points between said control modes; said optimal journey profile being determined by solving a system of differential equations for said speed profile of the train and for the value of said adjoint variable, said control modes being determined from the value of said adjoint variable, such that said sequence of control modes is determined as said speed profile 15 is calculated; (iii) monitoring the current state of said train as it progresses to said target location; and (iv) generating said information for the train operator by comparing the current state of the train to a corresponding state on said optimal journey profile and 20 displaying said information for the train operator to assist in keeping the train close to said optimal journey profile.

16. The method as claimed in claim 15 wherein said discrete control modes for said train include drive, hold, coast and brake modes.

17. The method of monitoring the progress of a train on a rail network as claimed 25 in claims 15 or 16 wherein said adjoint variable evolves according to a differential equation along with the position and speed of the train.

18. The method of monitoring the progress of a train on a rail network as claimed in claims 15 or 16 wherein the value of the adjoint variable is calculated directly from speed of the train. 30

19. The method of monitoring the progress of a train on a rail network as claimed in any one of claims 15 to 18 wherein a numerical method is used to solve said system - 22 of differential equations for said speed profile of the train and for the value of said adjoint variable.

20. The method of monitoring the progress of a train on a rail network as claimed in any one of claims 15 to 19 wherein steps (i) to (iv) are performed as required so that 5 said driving advice automatically adjusts to compensate for any operational disturbances encountered by the train.

21. The method of monitoring the progress of a train on a rail network as claimed in any one of claims 15 to 20 wherein said parameters include train mass and mass distribution. 10

22. The method of monitoring the progress of a train on a rail network as claimed in claim 21 wherein said parameters further include maximum tractive efforts and maximum braking effort as functions of speed.

23. The method of monitoring the progress of a train on a rail network as claimed in claims 21 or 22 wherein said parameters further include coefficient(s) of rolling 15 resistance.

24. The method of monitoring the progress of a train on a rail network as claimed in any one of claims 15 to 23 wherein said information is generated and displayed by a computer located on the train.

25. The method of monitoring the progress of a train on a rail network as claimed 20 in any one of claims 15 to 24 wherein step (iii) involves processing data from a GPS unit and train controls to determine the location and speed of the train.

26. The method of monitoring the progress of a train on a rail network as claimed in any one of claims 15 to 25 wherein said optimal journey profile specifies the time, speed and control at each location between the current train location and the next 25 target location on the network.

27. The method of monitoring the progress of a train on a rail network as claimed in any one of claims 15 to 26 wherein said optimal journey profile is precomputed.

28. The method of monitoring the progress of a train on a rail network as claimed in any one of claims 15 to 26 wherein a plurality optimal journey profiles 30 corresponding to different journey times are calculated and the profile that has an arrival time at the target location closest to the desired arrival time is selected. --23

29. A method of controlling the progress of a train on a rail network, said method comprising: (i) estimating or determining parameters of said train; (ii) determining, by an optimal control algorithm employing an adjoint 5 variable, an optimal journey profile for a journey from said train's current location to a target location that results in said train arriving at said target location as close as possible to a desired time and with minimum energy usage; said optimal journey profile including a speed profile for the train, sequence of discrete control modes for said train, and associated switching points between said control modes; 10 said optimal journey profile being determined by solving a system of differential equations for said speed profile of the train and for the value of said adjoint variable, said control modes being determined from the value of said adjoint variable, such that said sequence of control modes is determined as said speed profile is calculated; 15 (iii) monitoring the current state of said train as it progresses to said target location; and (iv) comparing the current state of the train to a corresponding state on said optimal journey profile and then controlling said train to keep the train close to said optimal journey profile. 20

30. The method as claimed in claim 29 wherein said discrete control modes for said train include drive, hold, coast and brake modes.

31. The method of controlling the progress of a train on a rail network as claimed in claims 29 or 30 wherein said adjoint variable evolves according to a differential equation along with the position and speed of the train. 25

32. The method of controlling the progress of a train on a rail network as claimed in claims 29 or 30 wherein the value of the adjoint variable is calculated directly from speed of the train.

33. The method of controlling the progress of a train on a rail network as claimed in any one of claims 29 to 32 wherein a numerical method is used to solve said system 30 of differential equations for said speed profile of the train and for the value of said adjoint variable. - 24

34. The method of controlling the progress of a train on a rail network as claimed in any one of claims 29 to 33 wherein steps (i) to (iv) are performed as required so as to automatically adjust to compensate for any operational disturbances encountered by the train. 5

35. The method of controlling the progress of a train on a rail network as claimed in any one of claims 29 to 34 wherein said parameters include train mass and mass distribution.

36. The method of controlling the progress of a train on a rail network as claimed in claim 35 wherein said parameters further include maximum tractive efforts and 10 maximum braking effort as functions of speed.

37. The method of controlling the progress of a train on a rail network as claimed in claims 35 or 36 wherein said parameters further include coefficient(s) of rolling resistance.

38. The method of controlling the progress of a train on a rail network as claimed 15 in any one of claims 29 to 37 wherein step (iii) involves processing data from a GPS unit and train controls to determine the location and speed of the train.

39. The method of controlling the progress of a train on a rail network as claimed in any one of claims 29 to 38 wherein said optimal journey profile specifies the time, speed and control at each location between the current train location and the next 20 target location on the network.

40. The method of controlling the progress of a train on a rail network as claimed in any one of claims 29 to 39 wherein said optimal journey profile is precomputed.

41. The method of controlling the progress of a train on a rail network as claimed in any one of claims 29 to 39 wherein a plurality optimal journey profiles 25 corresponding to different journey times are calculated and the profile that has an arrival time at the target location closest to the desired arrival time is selected.

42. A method of monitoring the progress of a train on a rail network and providing driving advice in real time to an operator of said train substantially as herein described with reference to any one of the embodiments of the invention illustrated in the 30 accompanying drawings and/or examples.

43. A method of monitoring the progress of a train on a rail network and providing information on the progress of the train in real time to an operator of said train - 25 substantially as herein described with reference to any one of the embodiments of the invention illustrated in the accompanying drawings and/or examples.

44. A method of controlling the progress of a train on a rail network substantially as herein described with reference to any one of the embodiments of the invention 5 illustrated in the accompanying drawings and/or examples.