Elevator control in an elevator system may enable real-location of a call after an initial call allocation. This means that, in some cases, it would be beneficial to reassign a new elevator to an already existing call. The situation is, however, different, for example, in elevator systems using immediate call allocation. One example of the immediate call allocation is a destination control system (DCS). In the immediate call allocation, already allocated calls are not typically be reassigned. This may lead to a situation that after allocating a call to a first elevator, it may turn out that it would be more beneficial and optimal to serve the call with a second elevator. But, as already mentioned above, it may not be possible to reassign the call of the already allocated call to the first elevator since in the destination control system the serving (allocated) elevator is signaled to the passengers immediately after giving a call.
Further, an elevator may become full before it has served all the calls and passengers assigned to it. This, on the other hand, may result in reduced passenger service level especially in destination control systems.
According to a first aspect of the invention, there is provided a method for computing allocation decisions in an elevator system. The method comprises obtaining historical passenger batch journey data relating to the elevator system, wherein each passenger batch journey comprises an origin and a destination floor of the journey, the number of passengers of the journey and the time of the journey; constructing historical passenger traffic statistics based on the passenger batch journey data; modelling expected calls based on the passenger traffic statistics; and taking the modelled expected call into account in computing subsequent allocation decisions in the elevator system.
In one embodiment, the method further comprises estimating elevator load in elevators of the elevator system based on the historical passenger traffic statistics; and taking the estimated elevator load into account in computing subsequent allocation decisions in the elevator system.
In one embodiment, alternatively or in addition, the method further comprises estimating elevator load in elevators of the elevator system separately for each elevator trip based on passenger batch journey intensities and batch size distributions obtained from the historical passenger traffic statistics, and simulated service times.
In one embodiment, alternatively or in addition, the modelled expected calls comprise at least one landing and car call pair.
In one embodiment, alternatively or in addition, the modelled expected calls comprise at least one destination call.
In one embodiment, alternatively or in addition, the passenger batch journey data comprises building origin-destination matrices formed separately for each day within a predetermined day cycle.
In one embodiment, alternatively or in addition, the elevator system uses immediate call allocation.
According to a second aspect, there is provided an elevator control apparatus for computing allocation decisions in an elevator system. The apparatus comprises at least one processor and at least one memory connected to the at least one processor. The at least one memory stores program instructions that, when executed by the at least one processor, cause the apparatus to obtain historical passenger batch journey data relating to the elevator system, wherein each passenger batch journey comprises an origin and a destination floor of the journey, the number of passengers of the journey and the time of the journey; construct historical passenger traffic statistics based on the passenger batch journey data; model expected calls based on the passenger traffic statistics; and take the modelled expected call into account in computing subsequent allocation decisions in the elevator system.
In one embodiment, the at least one memory stores program instructions that, when executed by the at least one processor, cause the apparatus to estimate elevator load in elevators of the elevator system based on the historical passenger traffic statistics; and take the estimated elevator load into account in computing subsequent allocation decisions in the elevator system.
In one embodiment, the at least one memory stores program instructions that, when executed by the at least one processor, cause the apparatus to estimate elevator load in elevators of the elevator system separately for each elevator trip based on passenger batch journey intensities and batch size distributions obtained from the historical passenger traffic statistics, and simulated service times.
In one embodiment, alternatively or in addition, the modelled expected calls comprise at least one landing and car call pair.
In one embodiment, alternatively or in addition, the modelled expected calls comprise at least one destination call.
In one embodiment, alternatively or in addition, the passenger batch journey data comprises building origin-destination matrices formed separately for each day within a predetermined day cycle.
In one embodiment, alternatively or in addition, the elevator system uses immediate call allocation.
According to a third aspect, there is provided a computer program comprising program code, which when executed by at least one processing unit, causes the at least one processing unit to perform the method of the first aspect.
In one embodiment, the computer program is embodied on a computer readable medium.
According to a third aspect, there is provided an elevator system comprising a plurality of elevators and an elevator control apparatus according to the second aspect.
The accompanying drawings, which are included to provide a further understanding of the invention and constitute a part of this specification, illustrate embodiments of the invention and together with the description help to explain the principles of the invention. In the drawings:
At 100 historical passenger batch journey data relating to the elevator system is obtained. Each passenger batch journey comprises an origin and a destination floor of the journey, the number of passengers (i.e. the passenger batch size) of the journey and the time of the journey. The passenger batch journey data provides historical, realized data about the usage of elevators in the elevator system.
At 102 historical passenger traffic statistics are constructed based on the passenger batch journey data. The historical passenger traffic statistics may be based on building origin-destination (OD) matrices which in turn may be based on the passenger batch journeys discussed above. A building specific OD matrix may be formed separately for each day d, d=1, 2, . . . , D, where D is the number of days in a cycle, and fixed intervals [tk, tk+1], k=1, 2, . . . , K, where K is the number of intervals within a day. For example, if D is set to “7”, the building specific OD matrix covers each weekday.
A N×N building specific OD matrix can be denoted as Λkdb, where N is the number of floors in the building and b, b=1, 2, . . . , B, denotes a batch size (i.e. the number of passengers). The element λijkdb in the building specific OD matrix corresponds to the intensity of passenger batch journeys equal to the batch size b from an origin floor i to a destination floor j within an interval k of a day d.
In the elevator group control, each candidate solution gives the allocation of the calls for the elevators in the group. In order to calculate the cost or fitness value of a solution, the service order of the calls or passengers for each elevator has to be determined. This can be done for each elevator independently from each other, for example, as follows.
Suppose there exists a given a set of calls, some of which can be dummy calls. First, the service order of the calls is determined by the collective control principle. Then, the route of the elevator through the calls is simulated in that order. The calls can be represented with nodes. A landing call corresponds to a single origin node and a destination call to an origin and destination node. Let N=[n1, n2, . . . , nk] denote the ordered set of nodes. The service time ti of a call n1 can be recursively calculated as ti=ti−1+τi−1,i, where τij is the time to go from a node ni to a node nj, that is, the flight time from the node ni to a node nj, plus the stop time at the node ni. A node n0 is the current location of the elevator, t0=0 and τ01, is the flight time between the current location of the elevator and the first node n1.
If ni corresponds to a landing call, existing or dummy, a set of dummy car calls that can be assumed to be given when the landing call ni is served, are modelled and added to the route in right places. Table 1 lists examples of different call types and items that can be modelled.
Once the simulation is over, the service time of each call is determined. Then the fitness value of each candidate solution, that is, allocation of calls, can be calculated using an objective function. A typical objective function is the average waiting time, the average journey time or the weighted sum of these two.
At 104 expected calls (or “dummy” calls) are modelled based on the passenger traffic statistics. Let Λkd=Σb=1BΛkdb denote a matrix containing the intensity of passenger batch journeys for each pair of floors within interval k of day d. An element λijkd is the intensity of journeys from an origin floor i to a destination floor j. It is also assumed herein that the batch journeys occur according to a Poisson process. Furthermore, the batch size distributions for each pair of floors may be defined by the matrices Λkd1, Λkd2, K, ΛkdB. By leaving out the interval and day indices, the time, tij, to the next pair of a dummy landing and car call, or a dummy destination call from an origin floor i to a destination floor j, can be defined as
where λij is the intensity of the batch arrivals occurring according to a Poisson process from an origin floor i to a destination floor j in seconds and γij is the time since the previous landing or destination call from the origin floor i to the destination floor j.
Because λij is the rate parameter of a Poisson process, 1/λij is the average time between two successive arrivals, i.e. calls. The above equation then implies that even if we assume that the batch arrivals occur according to a Poisson process, the forgetfulness property of the process is assumed only if the time since the previous call is longer than the predefined time limit {circumflex over (γ)}. A suitable value for the time limit can be determined, for example, with simulation studies.
When 1/λij<γij, i.e., there is a lot of traffic from the origin floor i to the destination floor j, and the time since the previous call is longer than the average inter-arrival time, then tij<0. In this case, it is natural to assume that the time to the next future call will be very short and thus we set tij=0. If tij∈[0,{circumflex over (t)}], where {circumflex over (t)} is a predefined time horizon, e.g., elevator cycle time, a pair of a dummy landing and car call, or a dummy destination call is generated from an origin floor i to a destination floor j with the arrival time equal to tcurrent+tij, where tcurrent is, e.g., the moment of computing a new allocation decision.
There can be several destination calls per origin floor, and thus, all dummy destination calls from an origin floor i to a destination floor j such that tij∈[0,{circumflex over (t)}] are generated. Because there can be only one landing call per direction on an origin floor i at a time, among all pairs of dummy landing and car calls to the same direction such that tij∈[0,{circumflex over (t)}], the pair for which the arrival time is closest to the current time can be selected.
Let li denote this arrival time, i.e. li=arg min{tij|0≤tij≤{circumflex over (t)}}.
Hence, the arrival time of the next dummy landing call on an origin floor i to the direction defined by the dummy car calls j such that tij∈[0,{circumflex over (t)}] is tcurrent+li. Further, for example, if there is an existing landing call at the origin floor i, then only the dummy car calls to destination floors j such that tij∈[0,{circumflex over (t)}] are generated.
At 106 at least one modelled expected (or “dummy”) call is taken into account in computing subsequent allocation decisions. This improves the service level of passenger since the allocation of elevator cars becomes more optimized.
At 108 load in elevators is modelled based on the historical passenger traffic statistics. For each dummy car call discussed in more detail in
μijkd=λijkdtiE└Yijkd┘
where ti is the serving time of the landing call on a floor I, and ti becomes defined during route simulation. E[Yijkd] is the expected number of passengers related to each arrival, in other words, the expected batch size which is estimated using the batch size distribution defined by the matrices Λkd1, Λkd2, K, ΛkdB, as already illustrated earlier. λijkd is the intensity of batch arrivals from an origin floor i to a destination floor j, that is, an element in the matrix defined as Λkd=Σb=1BΛkdb.
The intensities are estimated similarly as for dummy calls, as already illustrated in the description of
The estimated intensities with their origin and destination floor numbers may be stored in a memory. For existing car and destination calls that the elevator will actually serve, the intensities may be kept in the memory as long as they are not served at the destination floor. For dummy car and destination calls, the intensities may be kept in memory as long as it would still be possible to decelerate to the destination floor defined by the intensity. The destination floors of the estimated intensities can be represented with nodes.
When the simulation of an elevator route is over, all the destination nodes of the estimated intensities are iterated through in the order defined by the route. In course of the iteration, two rules may be used to separate the route into successive elevator trips. An elevator trip ends to the current destination node if:
In one embodiment, the following steps are performed:
After the above two steps, the intensities related to each elevator trip are used to define the origin and destination floors of the elevator trip.
In one embodiment, next the load of the elevator is estimated for each elevator trip individually as follows.
Let xij denote the number of individual passenger arrivals from a node i to a node j. Assuming that the arrivals occur according to a Poisson process, xij follows a Poisson distribution with the rate parameter μij. The number of passengers in the elevator after a stop at a node k must not exceed the capacity of the elevator. Mathematically this is can written as follows:
Σi=1kΣj=k+1Pxij≤Q
where Q is the elevator capacity and P is the number of nodes on the elevator trip.
Next, since xij may not be accurately known, the uncertainty related to them is taken into account by requiring that the probability of exceeding the elevator capacity Q after a stop at the node k is smaller than some predefined small probability α. Assuming that xij are independent, this can be expressed as follows:
where μ=Σi=1kΣj=k+1Pμij. The final equation above can be considered as a penalty term since the smaller the left hand side is, the more probable is that the elevator capacity will not be exceeded during the elevator trip. The penalty term for a single elevator trip can be written as follows:
where β is a scaling factor whose value can be determined, for example, with computational experiments. It follows that the penalty term for the whole route is the sum of the above penalties for the individual elevator trips. The penalty term for the whole route may thus be used when is added to an objective function used. A typical objective function is the average waiting time, the average journey time or the weighted sum of these two.
At 110 the at least one modelled expected (or “dummy”) call and the modelled load is taken into account in computing subsequent allocation decisions. An elevator group control is then able to construct and use historical passenger traffic statistics based on passenger batch journeys to estimate load of an elevator during its route through the calls allocated to it which, when taken into account in computing the allocation decisions, help to improve passenger service level.
In addition, it is assumed that two destination calls, one 204 from the floor 8 to the floor 2 and the other one 206 from the floor 5 to the floor 1, occur about the same time and that, based on minimizing the average waiting time, the first one is allocated to the first elevator 200 and the second one is allocated to the second elevator 202. Next, it is assumed that after about 5 seconds, when the elevators 200, 202 have started to move to serve the destination calls, a new destination call 208 is given from the floor 4 to the floor 6. This has been illustrated in
Now, from the average waiting time point of view, the best allocation would be to reassign the destination call from the floor 5 to the floor 1 for the first elevator 200 and to give the new destination call 208 to the second elevator 202. Normally, however, this would not be possible since in the destination control system the serving elevator is signaled to the passengers immediately after giving a call. However, as illustrated in the embodiment of
The elevator control apparatus 300 may be an elevator control entity configured to implement only the above disclosed operating features, or it may be part of a larger elevator control entity, for example, a group controller.
The exemplary embodiments of the invention can be included within any suitable device, for example, including, servers, workstations, personal computers, laptop computers, capable of performing the processes of the exemplary embodiments. The exemplary embodiments may also store information relating to various processes described herein.
Example embodiments may be implemented in software, hardware, application logic or a combination of software, hardware and application logic. The example embodiments can store information relating to various methods described herein. This information can be stored in one or more memories, such as a hard disk, optical disk, magneto-optical disk, RAM, and the like. One or more databases can store the information used to implement the example embodiments. The databases can be organized using data structures (e.g., records, tables, arrays, fields, graphs, trees, lists, and the like) included in one or more memories or storage devices listed herein. The methods described with respect to the example embodiments can include appropriate data structures for storing data collected and/or generated by the methods of the devices and subsystems of the example embodiments in one or more databases.
All or a portion of the example embodiments can be conveniently implemented using one or more general purpose processors, microprocessors, digital signal processors, micro-controllers, and the like, programmed according to the teachings of the example embodiments, as will be appreciated by those skilled in the computer and/or software art(s). Appropriate software can be readily prepared by programmers of ordinary skill based on the teachings of the example embodiments, as will be appreciated by those skilled in the software art. In addition, the example embodiments can be implemented by the preparation of application-specific integrated circuits or by interconnecting an appropriate network of conventional component circuits, as will be appreciated by those skilled in the electrical art(s). Thus, the examples are not limited to any specific combination of hardware and/or software. Stored on any one or on a combination of computer readable media, the examples can include software for controlling the components of the example embodiments, for driving the components of the example embodiments, for enabling the components of the example embodiments to interact with a human user, and the like. Such computer readable media further can include a computer program for performing all or a portion (if processing is distributed) of the processing performed in implementing the example embodiments. Computer code devices of the examples may include any suitable interpretable or executable code mechanism, including but not limited to scripts, interpretable programs, dynamic link libraries (DLLs), Java classes and applets, complete executable programs, and the like.
As stated above, the components of the example embodiments may include computer readable medium or memories for holding instructions programmed according to the teachings and for holding data structures, tables, records, and/or other data described herein. In an example embodiment, the application logic, software or an instruction set is maintained on any one of various conventional computer-readable media. In the context of this document, a “computer-readable medium” may be any media or means that can contain, store, communicate, propagate or transport the instructions for use by or in connection with an instruction execution system, apparatus, or device, such as a computer. A computer-readable medium may include a computer-readable storage medium that may be any media or means that can contain or store the instructions for use by or in connection with an instruction execution system, apparatus, or device, such as a computer. A computer readable medium can include any suitable medium that participates in providing instructions to a processor for execution. Such a medium can take many forms, including but not limited to, non-volatile media, volatile media, transmission media, and the like.
While there have been shown and described and pointed out fundamental novel features as applied to preferred embodiments thereof, it will be understood that various omissions and substitutions and changes in the form and details of the devices and methods described may be made by those skilled in the art without departing from the spirit of the disclosure. For example, it is expressly intended that all combinations of those elements and/or method steps which perform substantially the same function in substantially the same way to achieve the same results are within the scope of the disclosure. Moreover, it should be recognized that structures and/or elements and/or method steps shown and/or described in connection with any disclosed form or embodiments may be incorporated in any other disclosed or described or suggested form or embodiment as a general matter of design choice. Furthermore, in the claims means-plus-function clauses are intended to cover the structures described herein as performing the recited function and not only structural equivalents, but also equivalent structures.
The applicant hereby discloses in isolation each individual feature described herein and any combination of two or more such features, to the extent that such features or combinations are capable of being carried out based on the present specification as a whole, in the light of the common general knowledge of a person skilled in the art, irrespective of whether such features or combinations of features solve any problems disclosed herein, and without limitation to the scope of the claims. The applicant indicates that the disclosed aspects/embodiments may consist of any such individual feature or combination of features. In view of the foregoing description it will be evident to a person skilled in the art that various modifications may be made within the scope of the disclosure.
This application is a Continuation of PCT International Application No. PCT/FI2016/050441, filed on Jun. 17, 2016, which is hereby expressly incorporated by reference into the present application.
Number | Name | Date | Kind |
---|---|---|---|
4760896 | Yamaguchi | Aug 1988 | A |
5250766 | Hikita et al. | Oct 1993 | A |
6237721 | Siikonen | May 2001 | B1 |
6328134 | Hikita | Dec 2001 | B1 |
6345697 | Siikonen | Feb 2002 | B1 |
6360849 | Hikita | Mar 2002 | B1 |
7083027 | Siikonen | Aug 2006 | B2 |
10843895 | Salmikuukka | Nov 2020 | B2 |
20060289243 | Hikita et al. | Dec 2006 | A1 |
20120090922 | Elomaa | Apr 2012 | A1 |
20170081148 | Salmikuukka | Mar 2017 | A1 |
Number | Date | Country |
---|---|---|
1056659 | Dec 1991 | CN |
1301232 | Jun 2001 | CN |
1315284 | Oct 2001 | CN |
1802303 | Jul 2006 | CN |
101830375 | Sep 2010 | CN |
102190221 | Sep 2011 | CN |
103863912 | Jun 2014 | CN |
104891284 | Sep 2015 | CN |
0 090 642 | Oct 1983 | EP |
0452225 | Oct 1991 | EP |
10-236742 | Sep 1998 | JP |
4690703 | Jun 2011 | JP |
WO 9921787 | May 1999 | WO |
Entry |
---|
International Preliminary Report on Patentability, issued in PCT/FI2016/050441, dated Oct. 16, 2018. |
International Search Report, issued in PCT/FI2016/050441, dated Dec. 30, 2016. |
Kuusinen et al, “A new concept for passenger traffic in elevators”, Informs Annual Meeting Oct. 11-14, 2009, San Diego, USA [retrieved on Dec. 27, 2016]. Retrieved from the Internet: <URL: http://sal.aalto.fi/publicattons/ppt-files/jmkuu09s.ppt>, pp. 1-20. |
Ruokokoski et al. “Passenger Allocation to Capacitated Elevator Problem, manuscript 2.1 Oct. 22, 2008” [online][retrieved on Oct. 12, 2018]. Retrieved from the internet: <URL: http://sal.aalto.fi/publications/pdf-files/mruo08.pdf>, pp. 1-26. |
Written Opinion of the International Searching Authority, issued in PCT/FI2016/050441, dated Dec. 30, 2016. |
Number | Date | Country | |
---|---|---|---|
20190106289 A1 | Apr 2019 | US |
Number | Date | Country | |
---|---|---|---|
Parent | PCT/FI2016/050441 | Jun 2016 | US |
Child | 16214565 | US |