Blog

SUMMARY OF: ACE: An Actor Ensemble Algorithm for Continuous Control with Tree Search

• arXiv:1811.02696 [cs.LG]
• Source: [1] S. Zhang and H. Yao, “ACE: An Actor Ensemble Algorithm for Continuous Control with Tree Search,” Proc. AAAI Conf. Artif. Intell., vol. 33, no. 01, pp. 5789–5796, Jul. 2019.

Literature Review by S.Seal sayani.seal@mail.mcgill.ca

Edited by Q.Dang, D.Wu [di.wu5@mail.mcgill.ca]

Disclaimer: This website contains copyrighted material, and its use is not always specifically authorized by the copyright owner. Take all necessary steps to ensure that the information you receive from the post is correct and verified.

1.   motivation

In this paper [1], an actor ensemble learning algorithm is presented for continuous action space. In establishes improved performance over the Deep Deterministic Policy Gradient (DDPG) continuous control algorithm presented in [2]. For deterministic policies parameterized by neural networks, state-action value function is often maximized using gradient ascent method, which can get stuck at the local maxima. In [1], an actor-ensemble is used to mitigate this problem, where several actors parallelly search for the maxima in a localized manner starting from different initialization. This increase the likelihood of reaching the global maxima thereby, increasing the efficiency of the function approximator.

2.   Main Contributions

i. ACE, the actor ensemble continuous control algorithm presented in this paper is established to be a special class of the option-critic architecture [3].
ii. ACE extends the option-critic framework for deterministic intra-option policy. A deterministic intra-option policy gradient theorem and corresponding termination policy gradient theorem are proposed.
iii. The TreeQN [4] look-ahead tree search method is extended in the continuous action space and is implemented in ACE to improve the estimation of state-action value function. The multiple actors considered in the ACE are used as meta-actions to perform tree search using learned value prediction model.
iv. Compares and establishes improved performance of ACE over the baseline DDPG algorithm along with other modified implementations of this algorithm (Wide-DDPG, shared-DDPG, Ensemble-DDPG etc.) in an application in Roboschool.

3.   METHOD

A. The policy gradient method named the option-critic architecture (OCA) proposed by Bacon et al. [3] is treated as the basis of this ACE algorithm for continuous action control.

i. The OCA proposes a stochastic policy gradient which was presented using a framework of an option augmented Markov Decision Process (MDP), known as the semi-MDP. The ACE algorithm extends OCA to deterministic policy gradient framework.
ii. Derives deterministic intra-policy gradient theorem and corresponding termination gradient theorem which follows from the intra-option policy gradient theorem presented in the OCA paper [3]. The corresponding algorithm is named OCAD.
iii. By considering the corresponding termination function to 1, i.e., assuming that each policy terminates at every time step, ACE has been shown to be a special case of the OCAD architecture.
iv. Exploration noise is added for each action and experience-replay buffer is used in the ACE architecture similar to the DDPG algorithm. ACE also uses target network similar to DDPG. This leads to off-policy learning for the proposed neural network-based function approximator.

B. Model Based Enhancement
i. To improve the state-action value function used as the critic as well as the actors in the actor ensemble network, an extended TreeQN [4] for continuous control space is used to perform a look-ahead tree search method and learn using value prediction model.

4.   DATASET

Codes are available in:

5.    RESULT SUMMARIZATION

A.      Experiment Setup

i.            12 continuous control tasks from Roboschool by OpenAI.

ii.            States: Joint information of a robot presented as a vector of real number

Action: Vector, each with dimension varying between [-1,1]

B.      ACE architecture

i.            Latent stated ∈R 400

ii.            Encoding function: Single neural network layer

iii.            Reward prediction function: Single neural network layer

iv.            Transition function: Single neural network layer, residual connection used

v.            Value prediction function: Single neural network layer

vi.            Inputs: 400 latent states + m dimensional actions, 300 hidden units

vii.            Common first layer for actors with 400 latent state inputs and 300 hidden units

viii.            Activation function: tanh

ix.            Number of actors 5, panning depth 1.

C.       Baseline: ACE architecture is compared with

i.            ACE-Alt: Selected actors in ACE are updated at each time step

ii.            DDPG algorithm [2]: Instead of using ReLU with L₂ regularization to the critic tanh is used as an activation function without L₂ regularization.

 iii.            Wide-DDPG: To create equivalent number of parameters as compared to ACE, the number of hidden units were doubled in DDPG (i.e., two hidden layers with 800 and 600 units.)

 iv.            Shared-DDPG: Actor and critic shared a common bottom layer in DDPG, unlike the standard DDPG algorithm to create a comparable scenario with ACE since in ACE actor and critic shares a common representation layer, i.e. the encoder function TreeQN based Q-value estimator.

v.            Ensemble-DDPG: Tree search is removed from ACE by setting planning depth to 0, thereby generating an ensemble DDPG setup

vi.            Transition model ACE (TM-ACE): Impacts of TreeQN based value prediction model is investigated by learning a transition model instead of value prediction model.

 vii.            Different combination of actor size and planning depths are compared based on the performance.

D.      Results:

Table 1 from [1] shows the performance of ACE in comparison with the baseline.

The algorithm was trained for 1 million steps. At every 10 thousand steps 20 deterministic evaluation episodes were performed without exploration noise and the mean episode return were calculated. The performance improvement of ACE over DDPG are attributed to the actor ensemble architecture and the look-ahead tree search technique for actor and critic network update.

Figure 2 from [1] shows the trade-off for the choice of number of actors and planning depth. Since in ACE actor and critic shares the encodes latent state function approximator, increasing the number of actors dominates the training of the latent states which negatively impacts the critic network.

        Value prediction model encounters compounded error in unrolling. With increase in planning depth the error propagates through the network affecting the performance. Hence the both the planning depth and the number of actors needs to be moderated. Here, 5 actors and planning depth 1 achieves the best result as shown in the figure.

6.   SUGGESTED FUTURE WORK

Actor ensemble architecture can be implemented to train predictive models for continuous control applications.

REFERENCES:

[1]  S. Zhang and H. Yao, “ACE: An Actor Ensemble Algorithm for Continuous Control with Tree Search,” Proc. AAAI Conf. Artif. Intell., vol. 33, no. 01, pp. 5789–5796, Jul. 2019.

[2]  T. P. Lillicrap et al., “Continuous control with deep reinforcement learning,” CoRR, vol. abs/1509.02971, Feb. 2016.

[3]  P.-L. Bacon, J. Harb, and D. Precup, “The Option-Critic Architecture,” Sep. 2016.

[4]  G. Farquhar, T. Rocktäschel, M. I.- arXiv preprint arXiv …, and 2017, “TreeQN and ATreeC: Differentiable tree-structured models for deep reinforcement learning,” 2015.

[5]  V. Mnih et al., “Human-level control through deep reinforcement learning.,” Nature, vol. 518, no. 7540, pp. 529–533, Feb. 2015.

[6]  R. S. Sutton, D. A. McAllester, S. S.-A. in NIPS 2000, “Policy gradient methods for reinforcement learning with function approximation,” 1998.

————————–

<End of Review>

Written by,

Sayani Seal

sayani.seal@mail.mcgill.ca

Follow

SUMMARY OF: Deep Model Predictive Control with Online Learning for Complex Physical System

• arXiv:1905.10094v1 [cs.LG]
• Source: https://arxiv.org/abs/1905.10094 (May 2019).

Literature Review by S.Seal sayani.seal@mail.mcgill.ca

Edited by Q.Dang, D.Wu [di.wu5@mail.mcgill.ca]

Disclaimer: This website contains copyrighted material, and its use is not always specifically authorized by the copyright owner. Take all necessary steps to ensure that the information you receive from the post is correct and verified.

1.   motivation

Flow control is required in many fields of applications such as, energy, transportation, health and security. Though fluid flow has high-dimensional, multi-layer physics and nonlinear system characteristics, it can be approximated by some of the dominant low-dimensional system features. Since performance of a model predictive controller (MPC) significantly depends on the accuracy of its system prediction model, intractable complex systems pose difficulty in designing such a controller which is otherwise efficient for the particular application. The article [1] presents a DeepMPC controller where sensor-based observable low-rank system states are used to generate a recurrent neural network (RNN) based data-driven predictive system model for a real-time MPC implemented in fluid flow control.

2.   Main Contributions

 i.            DeepMPC architecture is implemented for complex fluid flow system exhibiting broadband phenomena.

ii.            Instead of using assumptions of full system states, the “surrogate” predictive system model uses only observable system states for future prediction. Thus, the method achieves a trade-off between accuracy and efficiency in capturing the essential physical system mechanisms.

iii.            The proposed learning approach for the RNN utilizes limited past information from the sensors.

Figure 1: DeepMPC with surrogate RNN prediction model presented in [1].

3.   METHOD

A.      DeepMPC

i.            Finite Open loop control problem with quadratic cost. Penalties assigned on deviation from reference trajectory, control input and any variation in the control input. The last component among the three restricts sudden change in the control input.

ii.            Surrogate system state prediction model, based on deep RNN architecture, is generated using control relevant observable sensor-based system states. For this flow control model, the states are lift and drag.

1.       RNN based predictive model design:

a.       Decoder:

i.        Performs actual prediction task

ii.        N-cell for N time steps in the prediction horizon

b.       Encoder:

i.        Predicts latent states and thereby accounts for long-term dynamics.

2.       RNN based MPC problem is solved using gradient based optimization method.

3.       The gradient information with respect to the control inputs is calculated using backpropagation-through-time.

 iii.            Training RNN:

1.       Offline three-stage training [2] with time-series data of observable system states.

2.       Training data, i.e., a time-series data of the lift and the drag, is generated using random but continuously variable control sequence of rotation force on cylinder(s).

4.   Result Summarization

A.      Setup:

A detailed simulation model of the full system is used instead of a real physical system. It is solved by OpenFOAM solver using finite volume discretization.

B.      Experiments:

Objective: The objective is to control the cylinder(s) such that

Four flow (laminar regime) control models with different complexity levels are considered:

i.            One cylinder: Flow around a single cylinder

1.       RNN prediction evaluated on exemplary control input sequence which showed accurate prediction for both lift and drag except for a very small duration at the start of the experiment.

2.       Successful showcase of tracking control of maintaining a schedules lift sequence for 20 sec with bounded rotation control input.

3.       Reynolds number () is assumed to be 100.

4.       Training dataset:

a.       Random rotation between -2 to +2 chosen at every 0.5 sec. Thus, high input frequencies are avoided

b.       Intermediate control inputs are computed using spline interpolation for every 0.1sec.

c.       A time-series with 110 000 datapoints are used for RNN training corresponding to 11 000sec.

ii.            Fluidic Pinball: Control the flow around three cylinders, two of which can be rotated the third one is fixed, as shown in Figure 1.

Figure 2: System is controlled by rotating cylinders 1 and 2 with respective angular velocities  and  [1].

1.       Objective is to follow three given lift trajectories for each cylinder by rotating cylinders 1 and 2.

2.        considered as the base case, other two chaotic cases with  and  are analyzed.

3.       Training dataset:

a.       Random rotation between -2 to +2 chosen for each cylinder at every 0.5 sec.

b.       Intermediate control inputs are computed using spline interpolation for every 0.005sec.

c.       Time series with 150 000, 200 000 and 800 000 are used for  and  respectively.  

4.       In order to improve performance for more chaotic systems with  and , knowledge regarding physical system characteristic is used by incorporating symmetric input and corresponding lift data along the horizontal axis. This reduces the tracking error by 50%.

5.       Robustness of the system is tested by performing five identical experiments with , using 10%, 15% and 100% of symmetrized training data points. No trend is observed with respect to the amount of training data.

Figure 3: DeepMPC lift tracking performance for laminar flow around rotating cylinders [1].

Figure 4: Re = 100 with online update [1].

6.       Finally, online data is collected from the feedback loop at each time step and new data collected over 25sec for each update. These 500 datapoints within each interval is used to further train the RNN surrogate model. This has significantly improved the performance of the DeepMPC as compared to (a) [1] in Figure 3. Online update of the RNN system reduces both tracking error and control cost.

5.   SUGGESTED FUTURE WORK

The surrogate RNN prediction model proposed for the DeepMPC in this article can be very usefully implemented for many practical engineering problems where the complete system description is too complicated and poses significant difficulty in solving related control problems. This method can be used for system modelling with targeted observable states which predominantly define respective system behaviour. This improve real-time implementation of MPC for complex nonlinear systems.

REFERENCES:

[1]    K. Bieker, S. Peitz, S. L. Brunton, J. K.- arXiv preprint arXiv, and 2019, “Deep model predictive control with online learning for complex physical systems,” 2012.

[2]    I. Lenz, R. Knepper, and A. Saxena, “DeepMPC: Learning Deep Latent Features for Model Predictive Control,” in Robotics: Science and Systems XI, 2015.

————————–

<End of Review>

Written by,

Sayani Seal

sayani.seal@mail.mcgill.ca

Follow

REVIEW ON: Markov Chain Monte Carlo Simulation of Electric Vehicle Use for Network Integration Studies

• Source: [1] Y. Wang, D. Infield, Markov Chain Monte Carlo simulation of electric vehicle use for network integration studies, International Journal of Electrical Power & Energy Systems, Vol.99, 2018, Pages 85-94

• https://doi.org/10.1016/j.ijepes.2018.01.008.

Literature Review by Q.Dang [qiyun.dang@mail.mcgill.ca]

Edited by D.Wu [di.wu5@mail.mcgill.ca]

Disclaimer: This website contains copyrighted material, and its use is not always specifically authorized by the copyright owner. Take all necessary steps to ensure that the information you receive from the post is correct and verified.

1. Paper Motivation

As the penetration of electric vehicles (EVs) increases, their patterns of use need to be well understood for future system planning and operating purposes. Using high resolution data by 10 minutes, accurate driving patterns were generated by a Markov Chain Monte Carlo (MCMC) simulation.
However, previous MCMC simulation works was not complete in the sense that model results were not subject to verification and uncertainty analysis for practical network assessment was not undertaken. The present paper includes both these important elements.

2. Methods

Method Name: Time-inhomogeneous Markov Chain Monte Carlo (MCMC) simulation

Description: The EV movement was simulated using a discrete-state, discrete-time Markov chain to define the states of all the EV at each time step of T minutes. It was assumed that, at every unit of time, one and only one event from a set of a finite number of events can occur to a given EV.

Four events were considered: {D, H, W, C}, correspond to ‘driving’, ‘parking at home’, ‘parking at workplace’, and ‘parking at commercial areas’ respectively

Proposed Markov Chain Diagram:

Fig. 1. Markov Chain diagram of possible vehicle state transitions at time t

From time step t-1 to t, the associated transition probability is given for each possible transition at this specific time stamp. For instance, P^t_H->D indicates the probability of the vehicle being ‘D’ at t given being ‘H’ at time t-1.

3. Paper structure

1) Review Previous Markov Chain Simulation of Electric Vehicle →

2) Introduce the survey data, the 2000 UK Time of Use Survey (TUS) data →

3) A matrix representation of the transition diagram at time t, Tt, is shown by Eq. (1) →

An example of the state transition matrix at 8:40 am (t = 29, t0=4am, 4am+29*10min=8:50am) is shown in Eq. (2),

→ Verification of proposed MCMC method by convergence analysis. →

4) Distribution grid case study by OpenDSS software (Case 1 commercial, Case 2 residential).

Case 1: A University building at Strathclyde, accommodates up to 300 workers, and has a nominal parking availability for approximately 100 cars. This building is supplied by a dedicated 1000 kVA transformer.

Case 2: low-voltage single-phase domestic network that consists of 17 households.

Fig. 2. Case 2 Single phase distribution network layout.

3.Paper Results

Results Description: 24 hour Load (KVA) profile in grids, before and after EV connected.

Fig. 3. (Upper)Aggregate demand of workplace EV charging. And (Lower) averaged voltage profile for Household 17 with 99% CI under full EV penetrations.

Notes:

Case 1 : An office building, approximately 100 cars, 100% EV penetration level, that is, 100 out of 100 cars are EV. This building is supplied by a 1000 kVA transformer.

For Case 1, a 1000 kVA transformer would easily survive the extra EV load for both standard and fast charging cases. A more typical transformer for this building with rating of nearer 500 kVA would, however, fail to supply the EV related load in the fast charging scenario.

Case 2: low-voltage single-phase community consists of 17 households.

For Case2, EV penetration in this case causes a severe voltage violation of the network (with specified tolerance of [−0.06 +1.10] p.u.,

4. Summarization

1) Markov Chain Monte Carlo simulation, as a numerical approach, can be used to generate different electricity load profiles according to various EV charging schemes.

2) The impact of the additional EV charging loads on the local distribution network can be assessed by identifying the expected value and associated uncertainty, as measured by the standard deviation, for various grid operational metrics, such as thermal loading, voltage profiles, transformer loss of life, energy losses, and harmonic distortion levels.

3) The uncertainty identification of these different metrics requires large number of trials from MCMC simulation to achieve convergence. These uncertainties could not be generated directly by sampling from the original TUS dataset due to its size limitation.

4) Also, the same steps of MCMC approach, as described in this work, can be applied to new data sets for extracting their own inherent statistical characteristics.

5. SUGGESTED FUTURE WORK

The EV movement was simulated using a discrete-state, discrete-time Markov chain for four events {D, H, W, C}, correspond to ‘driving’, ‘parking at home’, ‘parking at workplace’, and ‘parking at commercial areas’ respectively

The model can be extending to EV Charging States, including V2G and G2V, and further implemented in reinforcement leaning problems.

APPENDIX

• Useful Datasets download link (1&2):

1. National Household Travel Survey : http://nhts.ornl.gov/index.shtml

2. The United Kingdom 2000 Time Use Survey. National Statistics Technical Report; 2003.https://beta.ukdataservice.ac.uk/datacatalogue/studies/study?id=4504

Screenshot:

• Review of previous Markov Chain × EV works by author

Table 1. Summary of relevant literature works.

• A: Fine data resolution (less or equal to 10 min per step).
• B: Vehicle status definition.
• C: Vehicle movement simulation.
• D: Vehicle use pattern verification.
• E: Detailed network impact analyses considering charging location.
• F: Uncertainty analysis of detailed network impact.
• ✓: model feature is included in a suitable manner.
• ✗: model feature not included.
• —: not relevant.

References

[1] T.-K. Lee, Z. Bareket, T. Gordon, Z.S. FilipiStochastic modeling for studies of real-world PHEV usage: driving schedule and daily temporal distributions IEEE Trans Veh Technol, 61 (4) (May 2012), pp. 1493-1502

[2] F.J. Soares, J.P. Lopes, P.R. Almeida, C.L. Moreira, L. SecaA stochastic model to simulate electric vehicles motion and quantify the energy required from the grid PSCC, Stockholm, Sweden (2011)

[3]  Iversen EB, Møller JK, Morales JM, Madsen H. Inhomogeneous Markov models for describing driving patterns. IEEE Trans Power Syst.

[4]  A. Lojowska, D. Kurowicka, G. Papaefthymiou, L. van der Sluis Stochastic modeling of power demand due to EVs using copula IEEE Trans Power Syst, 27 (4) (2012), pp. 1960-1968

[5]  A. Ashtari, E. Bibeau, S. Shahidinejad, T. MolinskiPEV charging profile prediction and analysis based on vehicle usage data IEEE Trans Smart Grid, 3 (1) (2012), pp. 341-350

[6]  A.D. Hilshey, P.D. Hines, P. Rezaei, J.R. DowdsEstimating the impact of electric vehicle smart charging on distribution transformer aging IEEE Trans Smart Grid, 4 (2) (2013), pp. 905-913

[7]       F. Rassaei, W.S. Soh, K.C. ChuaDemand response for residential electric vehicles with random usage patterns in smart grids IEEE Trans Sustain Energy, 6 (4) (2015), pp. 1367-1376

[8]       Fluhr J, Ahlert KH, Weinhardt C. A stochastic model for simulating the availability of electric vehicles for services to the power grid. In: System Sciences (HICSS), 43rd Hawaii International Conference on. IEEE; 2010. p. 1–10.

[9]       S. Shafiee, M. Fotuhi-Firuzabad, M. RastegarInvestigating the impacts of plug-in hybrid electric vehicles on power distribution systems IEEE Trans Smart Grid, 4 (3) (2013), pp. 1351-1360

[10]     Wang Y, Huang S, Infield D. Investigation of the potential for electric vehicles to support the domestic peak load. In: Electric Vehicle Conference (IEVC), IEEE. Dec. 2014. p. 1–8.

————————–

<End of Review>

Written by,

Qiyun(Kevin) Dang

qiyun.dang@mail.mcgill.ca

Follow