jueves, 22 de junio de 2017

International Energy Markets,Understanding Pricing, Policies, and Profit, 2nd edition, 2015 - Международные энергетические рынки, понимание ценообразования, политики и прибыли


International Energy Markets,Understanding Pricing, Policies, and Profit, 2nd edition, 2015
Международные энергетические рынки, понимание ценообразования, политики и прибыли

Work in Progress: Supplemental Material to Accompany
International Energy Markets
Understanding Pricing, Policies, and Profit, 2nd edition, 2015
written by
cadahl@mines.edu
Senior Fellow

Study Aids
















Prerequisite Tools: For some material in the text, a basic understanding of algebra, calculus, probability and statistics, microeconomics, and macroeconomics is helpful. Here are some self tests to practice these tools:












Supplementary Material by Chapter


The Journey Starts


Chapter 1. Introduction


Introduction to Markets, Modeling and Policy


Chapter 2. Energy Lessons from the Past and Modeling the Future


Chapter 3. Perfect Competition and the Coal Industry






Chapter 4. Energy Taxes, Subsidies, and Social Welfare Market Structure


Chapter 5. Natural Monopoly and Electricity Generation
















Chapter 6. Deregulation-Privatization Electricity Sector


Chapter 7. Monopoly, Dominant Firm and OPEC


Chapter 8. Transaction Costs and U.S. Natural Gas Markets




Chapter 9. Monopsony - Japan and the Asia Pacific LNG Market




Chapter 10. Game Theory in W. European Natural Gas Market






Energy and the Environment


Chapter 11. Externalities and Energy Pollution




Chapter 12. Public Goods and Global Climate Change








Chapter 13. Energy Accidents


Dynamic Models and Model Inputs and Constraints


Chapter 14. Allocating Fossil Fuel Production over Time and Oil Leasing


Chapter 15. Supply and Costs


Chapter 16. Modeling Energy Demand


Chapter 17. Refining, Transportation, and Linear Programming




Managing Risk


Chapter 18. Energy Futures Markets for Managing Risk










Chapter 19. Energy Options Markets for Managing Risk




Development and Sustainability


Chapter 20. Climbing the Energy/Development Ladder to a Sustainable Energy Future


Chapter 21. Sustainable Wealth in Fossil Fuel Rich Developing Countries


The End of the Journey


Chapter 22. Managing in the Multicultural World of Energy












Dahl, Carol A. (2007) "What Goes Down Must Come Up: A Review of the Factors behind Increasing Gasoline Prices, 1999-2006." Prepared for the American Petroleum Institute. Washington, D.C.






Last modified Februrary 6, 2016



 

Python - Cuando Nature se toma el software libre en serio: Python como alternativa real a MATLAB, Mathematica y otros programas comerciales, Когда природа свободное программное обеспечение серьезно



Cuando Nature se toma el software libre en serio: Python como alternativa real a MATLAB, Mathematica y otros programas comerciales

 http://cacheme.org/charla-nature-software-libre-python-alternativa-matlab-mathematica/

Python - Когда природа свободное программное обеспечение серьезно: Python в качестве реальной альтернативы MATLAB, Mathematica и других коммерческих программ

 

No sólo una, sino dos Tres son ya las noticias de la revista Nature donde se habla y recomienda Python para uso científico.  Si acabas de llegar y no conoces este mundillo aún, Python es un lenguaje de programación fácil de aprender, multiplaforma, libre y gratuito que está sustituyendo a software comercial (MATLAB, Mathematica, MAPLE…) en universidades norteamericanas tales como MIT, Stanford o Caltech, y que se utiliza ampliamente por empresas como Google, Dropbox o la NASA.
DiapositivasCurso de Python científico (online)

Motivación de la charla:

La famosa revista científica Nature empezó en noviembre del año pasado a promocionar una serie de herramientas libres, destacando Python entre ellas:
Nature - Notebooks interactivos: compartiendo el código
Nature: Python para para estudio del clima
¿Quieres aprender a programar? Elige Python
A nivel académico, Jupyter —la evolución de IPython Notebook—  puede mejorar la calidad de publicaciones científicas ya que permite la reproducibilidad de resultados de forma inmediata y libre. Si bien la idea de compartir el código no es nueva, Jupyter permite que cualquier persona sea capaz de reproducir tus resultados a golpe de ratón y, si lo desea, probar con otros datos, indagar en el código y/o modificarlo. ¿No te lo crees? Puedes comprobarlo online con este ejemplo de Nature 🙂
Demo interactiva de Nature y Jupyter
Demo interactiva de Nature y Jupyter (IPython Notebook), nunca había sido tan fácil detectar galaxias.
La tercera iteración de IPython, Jupyter, se encuentra en fase beta en estos momentos y soporta de forma oficial Python, R y Julia; aunque también abre la posibilidad de uso con otros lenguajes de programación. Además siendo este software libre y gratuito, Google docs ya avanzó que lo implementaría como si se tratase de un formato más dentro de sus servicios online y haciendo uso de las opciones —diseñadas en Jupyter— de trabajar de forma colaborativa con él. Es decir, que varios usuarios podrán trabajar juntos en un mismo documento interactivo y alojado en la nube. Todas estas ventajas tampoco han pasado desapercibidas en la docencia, los famosos MOOC (plataformas de cursos online abiertos) como Coursera y Edx tiene varios cursos de Python —aunque quizás el ejemplo más destacable sea el magnífico curso online de Métodos numéricos con Python de Lorena Barba y colaboradores—.

XCOS online Alternativa de Simulink - Simulacion Diagramas de Bloques, Xcos онлайн Alternative Simulink - Моделирование блок-схема

XCOS online Alternativa de Simulink - Simulacion Diagramas de Bloques

Xcos онлайн Alternative Simulink - Моделирование блок-схема

http://xcos.fossee.in/

Excelente Explicación de PID con diagramas de bloques y simulaciones con Escalon - Отлично Объяснение PID блок-схем и моделирования Эскалон


Excelente Explicación de PID con diagramas de bloques y simulaciones con Escalon

Отлично Объяснение PID блок-схем и моделирования Эскалон

http://scilab.ninja/study-modules/scilab-control-engineering-basics/module-4-pid-control/

Module 4: PID Control

This article is contained in Scilab Control Engineering Basics study module, which is used as course material for International Undergraduate Program in Electrical-Mechanical Manufacturing Engineering, Department of Mechanical Engineering, Kasetsart University.

Module Key Study Points

  • understand the basics of PID control
  • learn the effects of 3 parameters to system response
  • learn how response is degraded by integrator windup
  • How to tune PID parameters by Ziegler-Nichols Frequency Domain method
PID stands for Proportional, Integral, and Derivative. It is a standard control structure used successfully in many industrial applications. Commercial PID controllers can be bought off the shelf and installed in the system. After some setup and tuning the three PID gains, and perhaps some additional parameters, the system is up and running in no time. Of course, this means the control engineer is certain the plant can be handled by the PID controller and understand how to adjust the control parameters. Without such knowledge, it could be a frustrating experience to her.
So, in this module we discuss some basics of PID control and focus on how to analyze a feedback sysem, with Scilab, build an Xcos diagram and simulate. For a detailed study, we recommend [1]. For those interested in digital PID implementation, document [3] on this site might be helpful. The e-book [4] provides more information on independent robot joint control using PID.

Different forms of PID

The so-called “textbook” form of PID controller is described as follows
(1)   \begin{equation*}  u(t) = K\left( e(t) + \frac{1}{T_i}\int_{0}^t e(\tau)d\tau + T_d\frac{de(t)}{dt}\right) \end{equation*}
This algorithm computes the control variable u as output, given the input e, the error between the command and plant output. We see that the control variable is a function of 3 terms: P (proportional to error), I (time integral of error), and D (derivative of error), with corresponding control parameters K (proportional gain), T_i (integral time), and T_d (derivative time), respectively. Taking Laplace transform of (1), we have
(2)   \begin{equation*}  u(s) = C(s)e(s) \end{equation*}
with the transfer function of PID controller
(3)   \begin{equation*}  C(s) = K\left( 1 + \frac{1}{sT_i} + sT_d\right) \end{equation*}
Another common structure of PID algorithm is represented by
(4)   \begin{equation*}  u(t) = K_pe(t) + K_i\int_{0}^t e(\tau)d\tau + K_d\frac{de(t)}{dt} \end{equation*}
In this form, the controller gains are distributed to each of the PID terms separately, with its transfer function
(5)   \begin{equation*}  C(s) = K_p + \frac{K_i}{s} + K_ds \end{equation*}
For our discussion of PID controller, we focus on the form (4) and (5). Nevertheless, it is easy to convert between parameters of (1) and (4) by the relations
(6)   \begin{equation*}  K_p = K, K_i = \frac{K}{T_i}, K_d = KT_d \end{equation*}
Ex. 1: To get started with Xcos simulation of a PID feedback system, download file pid_feedback.zcos . This launches an Xcos diagram in Figure 1. This diagram can easily be constructed using Xcos standard palettes. It requires only a step_function, clock_c (from Sources palette), PID and CLR blocks (from Continuous time systems palette), SUMMATION (from Mathematical operations palette), MUX (from Signal Routing palette), and CSCOPE (from Sinks palette). The plant transfer function is a robot joint driven by DC motor model from module 1, with J = 10 and B = 1.



Figure 1 pid_feedback.zcos basic PID simulation diagram
Click on the Simulation Start button to see the step response. Try adjusting the PID gains and plant parameters to see their effects on the closed-loop system.
PID parameter tuning is well studied in control engineering. Here we give a brief guideline. The proportional gain K_p is a dominant quantity that normally has some nonzero value. The integral gain K_i helps eliminate steady state error but too high value could introduce overshoot and oscillation. The derivative gain K_d could help the response to reach steady state faster but could amplify high frequency noise, and could affect stability if set too high.
When the plant model is not known, adjusting these 3 gains to achieve good response could be problematic for an inexperienced user. Later we discuss a tuning procedure. Commercial PID controller products usually have auto-tuning functions for user convenience.
Ex. 2: To experiment with the tracking and disturbance attenuation performance of PID control, we setup a feedback diagram in Figure 2.



Figure 2 PID feedback diagram
which can be constructed as an Xcos model pid_dist.zcos in Figure 3. The joint inertia and friction are set to J = 10 and B = 0.1, respectively. The step disturbance at the plant input is active at t = 2 sec, with step value of 80.



Figure 3 pid_dist.zcos Xcos model for PID feedback with disturbance
Note that in this example, we have the luxury of knowing the exact plant model. So it is possible to solve for a good set of PID gains. An approach used in [2] turns off the K_i gain first and solve for K_p and K_d that yields a critical-damped system at chosen natural frequency. Then K_i is adjusted to eliminate steady-state error.
To elaborate, it is shown in [2] that when K_i is set to 0, the resulting second order closed-loop characteristic polynomial is described by
(7)   \begin{equation*}  s^2 + \frac{(B + K_d)}{J}s + \frac{K_p}{J} = s^2 + 2\zeta\omega s + \omega^2 \end{equation*}
which gives
(8)   \begin{equation*}  K_p = \omega^2J, K_d = 2\zeta\omega J - B \end{equation*}
For example, suppose we select a critical-damped closed loop system \zeta = 1, with natural frequency \omega = 8. From (8), this yields K_p = 640 and K_d = 160. Table 1 below gives calculated P and D gain values for 3 selected natural frequencies.

\omega K_p K_d
4 160 80
8 640 160
12 1440 240

Table 1 K_p and K_d gains for 3 different natural frequencies
Using these 3 sets of gain values, we simulate the step responses. The results are shown in Figure 4. As expected, the higher the natural frequency, the faster the response. We also see that the steady-state error from step disturbance cannot be eliminated with the PD controller.



Figure 4 step responses of PD control with disturbance d = 80
To fix this problem, we need to turn on the K_i gain of PID controller. With the controller transfer function described by (5), the closed-loop system is now of third order
(9)   \begin{equation*}  \theta (s) = \frac{K_ds^2 + K_ps + K_i}{\Lambda (s)}\theta_d(s) - \frac{s}{\Lambda (s)}D(s) \end{equation*}
with the characteristic polynomial
(10)   \begin{equation*}  \Lambda (s) = Js^3 + (B + K_d)s^2 + K_ps + K_i  \end{equation*}
Applying Routh-Hurwitz criteria to (10), it can be concluded that, given all positive PID gains, the feedback system is stable if
(11)   \begin{equation*}  K_i < \frac{(B + K_d)K_p}{J} \end{equation*}
The same Xcos diagram in Figure 3 can be used by putting some positive integral gain in the PID block. Let us fix K_p = 1440 and K_d = 240, corresponding to natural frequency \omega = 8 rad/s. Figure 5 shows the responses from setting some values of K_i. We see that the steady-state error to command and disturbance is eliminated by increasing K_i, though higher value could introduce more overshoot.



Figure 5 step responses of PID control
Applying the bound (11) to this example yields
    \begin{displaymath}  K_i < \frac{(B + K_d)K_p}{J} = \frac{(0.1 + 240)1440}{10} = 34574.4 \end{displaymath}
This can be verified by setting K_i = 34575 . The oscillating response in Figure 6 shows that the closed-loop system is at the brink of stability.



Figure 6 closed-loop response at controller gains K_p = 1440, K_i = 34575, K_d = 240

Effect from Saturation

Analysis of PID feedback control is often performed under the assumption that the system is purely linear. In real devices, nonlinear elements emerge at certain points in the feedback loop. Saturation is one of them that could cause undesirable system response. This nonlinear effect normally results from physical limits of signals and system parameters, say, maximum motor torque and current. Another common saturation is the input limit of servo amplifier. An industrial servo amp might allow input voltage command in the range \pm 10 volts, for example.
Effect from saturation is more pronounced when integral term is used in PID. The response could have excessive overshoot due to error accumulation in the integral term. This is known as integrator windup. See [1] for more detail.
Ex. 3: To see an effect from input saturation, construct an Xcos model as in Figure 7, or download pid_ilim.zcos. It represents a comparison of two closed-loop systems, which are basically the same except that the lower one has its input bounded between \pm 50 units.



Figure 7 pid_ilim.zcos Xcos diagram for Example 3
The controller gains are set at K_p = 640, K_i = 1000, K_d = 160 . Running the simulation yields the step responses in Figure 8, and the control variables (controller output) in Figure 9. We see that input saturation causes significantly larger overshoot and also worsens disturbance response.



Figure 8 effect of integrator windup on output response



Figure 9 controller output with plant input limit

Tuning the PID Parameters

As stated before, adjusting the PID gains from scratch to achieve a good response may not be trivial. As a result, most commercial PID controllers have functions to tune the 3 parameters automatically. This is normally called “autotuning” feature. One of autotuning methods suggested in the literature makes use of some relay feedback mechanisms, which is closely related to a manual tuning scheme known as Ziegler-Nichols Frequency Domain (ZNFD) method. Hence we discuss this ZNFD scheme in the last part of this module.
Note: The original ZNFD method applies the “textbook” PID equation (1). But since (1) and (4) are closely-related by (6), here we make a conversion so that standard PID block in Xcos, which assumes the form (4), can be conveniently used.
To tune a PID controller manually by ZNFD method, follow this procedure
  • turn off both the integral and derivative terms; i.e., setting K_i = K_d = 0. So now the PID is left only with the proportional gain K. We crank K up to the point that the closed-loop system starts to oscillate. At this point, the plant output will swing in a constant sinusoid motion, not growing and not dying out. Write this value down on a paper as K_u.
  • Measure the period of oscillation. Write it down as T_u.
  • Use Table 2 to tune the controller parameters
Controller Form K_p K_i K_d
P 0.5K_u
PI 0.4K_u 0.5K_u/T_u
PID 0.6K_u 1.2K_u/T_u 0.075K_uT_u
Table 2: suggested PID parameters from ZNFD method
From the above procedure, we could observe a drawback of the basic ZNFD method: the system must oscillate at some point of gain setting before it becomes unstable. This implies the Nyquist plot must cross the negative real axis at some point. For the DC motor model we are using, the gain margin is infinity. Hence it is difficult to find the oscillating point and the K_u value.
There is some trick to get around such limitation, but this is beyond the scope of this module. For our basic study, we simply modify the DC motor robot joint by adding some dynamics to it, so that the Nyquis plot of L(j\omega) crosses the negative real axis.
Ex. 4: Let us assume that the PID output is passed through a LPF with cutoff frequency 500 Hz before feeding the plant input. Using the filter design guide from module 1, we have the LPF transfer function as
(12)   \begin{equation*}  H(s) = \frac{1}{1 + 0.0003s} \end{equation*}
To simulate, construct an Xcos model in Figure 10, or download pid_znfd.zcos. With this model, we are ready to perform the ZNFD tuning procedure described above. Click on the PID block to set the I and D gains to zero. Then start increasing the P gain until the output oscillates. This is what you must do in a real application when the plant model is not known.



Figure 10 pid_znfd.zcos Xcos diagram for ZNFD tuning method
Well, in this example we do know the plant model, don’t we? So, to save us some time, the K_u value that causes output oscillation can be found, for example, using the method from Module 2. First from the loop transfer function containing the plant, LPF, and PID controller with K_p = 1.
-->kp = 1;
-->ki = 0;
-->kd = 0;
 
-->s=poly(0,'s');
-->P=syslin('c',1/(10*s^2+0.1*s));  // plant
-->H = syslin('c',1/(1+0.0003*s));  // filter
-->C=syslin('c',kp + ki/s + kd*s);  // controller
-->L = C*H*P; // loop t.f
Then find the gain where the closed-loop pole touches the j\omega axis.
->[kmax,s]=kpure(L)
 s  = 
    5.7735027i  
 kmax  = 
    333.33433
resulting in K_u = 333.
Verify this by setting the P gain to 333 in PID block. This yields the desired oscillatory response as shown in Figure 11. The oscillation period can be roughly measured from the plot to yield T_u = 1 sec.



Figure 11 oscillatory response with K_p = 333, K_i = K_d = 0
So from the last row of Table 2, ZNFD method suggests the three PID parameters as
    \begin{displaymath} K_p = 0.6 \times 333 = 200, K_i = 2 \times 200 = 400, K_d = 333/8 = 42. \end{displaymath}
Setting these gains gives the response in Figure 12.



Figure 12 response from PID gains suggested by ZNFD method
Note that the overshoot is quite excessive (70%). In a sense, ZNFD just gives us some good values to start with. We may want to fine-tune the PID gains to improve the response further. For example, decreasing the Ki gain would bring the overshoot down. Figure 13 shows the response from original PID gains (green), compared with the case when Ki is reduced to 300 (blue), and 200 (magenta).



Figure 13 reduce overshoot in the response by decreasing K_i
The ZNFD method could be explained using a Nyquist diagram in Figure 14. The diagram shows how a point x on the curve is moved related to the P, I , and D terms. Using the P term alone, x could be moved in radial direction only. The I and D terms help provide more freedom to move perpendicular to the radius. It can be shown that by using ZNFD method, the critical point (-1/K_u, 0) is moved to the point -0.6 - 0.28i. The distance of this point to the critical point is 0.5. So the sensitivity peak is at least 2. This explains the high overshoot in the step response.



Figure 14 How a point on Nyquist curve is moved with PID control

Summary

PID is a simple control structure that is still used widely in various industrial applications. It is a close relative to the lead-lag compensator explained in module 2, except that its functionality may be more user-friendly. With some knowledge and practice, an engineer or technician would be able to tune and operate a plant equipped with PID control.
In this module we discuss the basics of PID feedback systems, with emphasis on Xcos simulations to show how the responses are related to three control parameters, as well as effect from input saturation that could worsen the response. Without some good starting values, tuning the PID gains can be cumbersome for a novice. So at the end, we mention a manual tuning procedure known as the Ziegler-Nichols frequency domain method. Some auto-tuning scheme of a commercial PID controller, such as the relay feedback method, is based on the ZNFD manual tuning.

References

  1. K.J. Astrom and T.Hagglund. PID Controllers, 2nd ed., Instrument Society of America, 1995.
  2. M.W.Spong, S. Hutchinson and M. Vidyasagar, Robot Modeling and Control. John Wiley & Sons. 2006.
  3. V. Toochinda. Digital PID Controllers, a previous online document available at controlsystemslab.com. 2009.
  4. V.Toochinda. Robot Analysis and Control with Scilab and RTSX. Mushin Dynamics, 2014.

Comments

Module 4: PID Control — 2 Comments