Anlisis de Sistemas de Potencia II UNI-FIEE
Mtodo del Gradiente Aplicado al Despacho Econmico
Mtodos de optimizacin El despacho econmico de carga es un elemento importante dentro de la produccin, planeacin y control de un sistema elctrico de potencia, por lo tanto se hace necesario buscar un mtodo matemtico que proporcione una solucin ptima al problema del despacho. A continuacin se presenta el mtodo del gradiente.
MTODO DEL GRADIENTE PARA EL DESPACHO ECONMICO EN GENERACION TERMICA
La funcin a minimizar es la funcin de costo total, teniendo como restriccin solamente el equilibrio de potencia, ver ecuaciones:
( )
Para resolver el problema del despacho econmico, aplicamos la tcnica del Gradiente a la funcin de Lagrange que se define a continuacin: ( ) Y la gradiente de la funcin de Lagrange ser: ( )
El procedimiento del despacho econmico requiere valores inciales de , P1, P2, P3. Pi para encontrar el gradiente de L luego se encuentra los nuevos valores de , P1, P2, P3. Pi los cuales se encuentran de la siguiente manera:
( )
Hernando Cruz Astete - UNI - ING. ELECTRICA
Anlisis de Sistemas de Potencia II UNI-FIEE
Donde x es el vector de valores inciales y es un escalar que nos garantiza que el proceso De convergencia. El vector x tiene la siguiente estructura:
[ ] El proceso converge cuando el L es igual a cero, lo que indica que se ha encontrado el valor ptimo del despacho. ( ) ( )
INICIO
LEER P C , Fi
O , 0.1
(k )
P ( k ) 1 2ci
Pi
( k ) bi
2ci
( k 1) ( k 1) ( k )
p ( k ) PC Pi ( k )
NO
p ( k )
SI
Pi , Fi , , FT
FIN
Hernando Cruz Astete - UNI - ING. ELECTRICA
Anlisis de Sistemas de Potencia II UNI-FIEE
DIAGRAMA DE FLUJO METODO DEL GRADIENTE
PROGRAMACIN EN MATLAB METODO DEL GRADIENTE
display('ANALISIS DE SISTEMAS DE POTENCIAII UNI-FIEE 2014'); display('METODO DEL GRADIENTE PARA DESPACHO ECONOMICO ') display('FUNCION DE COSTO DEL GENERADOR 1') F1=input('funcion costo gen1 [a1,b1,c1]='); display('FUNCION DE COSTO DEL GENERADOR 2') F2=input('funcion costo gen2 [a2,b2,c2]='); display('FUNCION DE COSTO DEL GENERADOR 3') F3=input('funcion costo gen3 [a3,b3,c3]='); % poner aqu los datos CARGA=input('ingrese el valor de la carga en MW '); display('LIMITES DE GENERACION DE POTENCIA') Lmin1=input('ingrese el LIMITE MINIMO DE P1 EN MW '); Lmax1=input('ingrese el LIMITE MAXIMO DE P1 EN MW '); Lmin2=input('ingrese el LIMITE MINIMO DE P2 EN MW '); Lmax2=input('ingrese el LIMITE MAXIMO DE P2 EN MW '); Lmin3=input('ingrese el LIMITE MINIMO DE P3 EN MW '); Lmax3=input('ingrese el LIMITE MAXIMO DE P3 EN MW '); display('VALORES INCIALES') P1=input('ingrese el valor de P1 EN MW= '); P2=input('ingrese el valor de P2 en MW= '); P3=input('ingrese el valor de P3 en MW= '); a=input('ingrese el valor de alfa= '); n=input('ingrese el nmero de iteraciones del programa '); %DERIVADAS DE LAS FUNCIONES y1=polyder(F1); y2=polyder(F2); y3=polyder(F3); for i=1:n; alfa=a; %valores de lambda (i); lambda1=polyval(y1,P1); lambda2=polyval(y2,P2); lambda3=polyval(y3,P3); lambda_t=1/3*(lambda1+lambda2+lambda3); %Primera iteracin; P1(i)=P1-(lambda1-lambda_t)*alfa; P2(i)=P2-(lambda2-lambda_t)*alfa; P3(i)=P3-(lambda3-lambda_t)*alfa; lambda_t(i)=lambda_t-(CARGA-(P1(i)+P2(i)+P3(i)))*alfa; if P1(i)<Lmin1; P1(i)=Lmin1; end if P1(i)>Lmax1; P1(i)=Lmax1; end if P2(i)<Lmin2; P2(i)=Lmin2; end if P2(i)>Lmax2; P2(i)=Lmax2; end if P3(i)<Lmin3; P3(i)=Lmin3; end if P3(i)>Lmax3;
Hernando Cruz Astete - UNI - ING. ELECTRICA
Anlisis de Sistemas de Potencia II UNI-FIEE
P3(i)=Lmax3; end P1=P1(i) P2=P2(i) P3=P3(i) lambda_t=lambda_t(i) end
Analizando el programa: Mediante el comando polyder se extraen las derivadas de F1 con respecto a P1, De F2 con respecto a P2 y la derivada de F3 con respecto a P3. %DERIVADAS DE LAS FUNCIONES y1=polyder(F1); y2=polyder(F2); y3=polyder(F3); Cada elemento del vector gradiente es igual a cero, de modo que se establecen los valores de lambda1, 2 y 3, El valor de lambda1 por ejemplo es igual a la derivada de F1 con respecto a P1 y evaluada en P1, lo mismo se aplica para obtener lambda 2 y lambda 3. De los tres valores de lambda extraigo el promedio que se convierte en mi lambda inicial
Ejemplo 1:
VALORES TOMADOS DEL LIBRO POWER GENERATION OPERATION AND CONTROL (ALLEN J. WOOD) FUNCIONES DE COSTO F= [ c b a ]
F1= [0.001562 ,7.92, 561] F2= [0.00194 ,7.85 ,310] F3= [0.00482 ,7.97 ,78] Carga= 800 MW LIMITES DE GENERACIN
Valores inciales
Hernando Cruz Astete - UNI - ING. ELECTRICA
Anlisis de Sistemas de Potencia II UNI-FIEE
Una mirada al matlab , insertando los datos del ejemplo 1
Hernando Cruz Astete - UNI - ING. ELECTRICA
Anlisis de Sistemas de Potencia II UNI-FIEE
Los valores de convergencia son : P1=370.0021 MW P2=315.3876 MW P3=114.6103MW Lambda = 9.0748
Hernando Cruz Astete - UNI - ING. ELECTRICA
Anlisis de Sistemas de Potencia II UNI-FIEE
Ejemplo 2:
Problema #4 tomado de la 1ra practica
FUNCIONES DE COSTO
F1=[0.001962,7.92,561] F2=[0.00234,7.85,310] F3=[0.00522,7.97,78] Carga= 900 MW
LIMITES DE GENERACIN
Valores inciales
Hernando Cruz Astete - UNI - ING. ELECTRICA
Anlisis de Sistemas de Potencia II UNI-FIEE
Los valores de convergencia son : P1=401.5881 MW P2=352.1415 MW P3=146.2704 MW Lambda = 9.4970
Hernando Cruz Astete - UNI - ING. ELECTRICA
