*******OPTIMIZACIÓN DE LA CADENA DE SUMINISTRO PARA LA OBTENCIÓN DE BIOCARBÓN A PARTIR DEL SUPER SORGO******* SET i cultivo de super sorgo /crop/ j secciones de cultivo /s1*s46/ k biorrefinería /biorref/ l almacenes /dep/ m mercado /markCARBON, markFOOD, markFUEL/ ALIAS (j,j2) ALIAS (i,i2,i3) ALIAS (k,k2) ALIAS (l,l2) ALIAS (m,m2) PARAMETERS u1(m) demanda (ton) / markCARBON 6912 markFOOD 17280 markFUEL 10368 / u2(i) rendimiento de cultivo (ton*ha-1) / crop 275.2967/ u3(i) área mínima a cultivar (ha) / crop 32/ u4(i) costo de producción primer corte ($*ha-1) / crop 22300/ u5(i) costo de producción segundo corte ($*ha-1) /crop 10700/ u6(i) cantidad de carbono fijo /crop 0.2/ u7(l) capacidad de los almacenes ton*cuatrimestre-1 /dep 85/ u8(l) costo de almacenamiento /dep 15/ u9(i) valor de la producción ($*ha-1) / crop 89376/ B(m) factor de segregación para los mercados de biocarbón / markCARBON 0.8105 markFOOD 0.0261 markFUEL 0.1634 / v(j2) cantidad de hectareas disponibles por sección / s1 6 s2 8 s3 2 s4 10 s5 3 s6 4 s7 2 s8 3 s9 14 s10 5 s11 6 s12 2 s13 3 s14 4 s15 9 s16 20 s17 14 s18 2 s19 6 s20 3 s21 5 s22 6 s23 3 s24 8 s25 4 s26 8 s27 14 s28 3 s29 14 s30 2 s31 15 s32 8 s33 6 s34 3 s35 4 s36 21 s37 8 s38 9 s39 14 s40 22 s41 2 s42 8 s43 6 s44 7 s45 4 s46 21 / table d1(j,k) distancia de cada seccion a la biorrefiniería biorref s1 85 s2 63 s3 22 s4 45 s5 44 s6 89 s7 10 s8 36 s9 85 s10 20 s11 38 s12 85 s13 41 s14 14 s15 36 s16 63 s17 98 s18 71 s19 85 s20 40 s21 54 s22 35 s23 84 s24 12 s25 36 s26 21 s27 14 s28 52 s29 25 s30 55 s31 14 s32 25 s33 96 s34 32 s35 14 s36 55 s37 21 s38 27 s39 38 s40 95 s41 35 s42 33 s43 28 s44 98 s45 74 s46 85 ; table c1(j,k) costo de transporte ($*km-1) de cada sección a la biorrefinería biorref s1 1 s2 2 s3 0.9 s4 0.8 s5 5 s6 2 s7 3 s8 1 s9 2 s10 2 s11 1 s12 3 s13 0.3 s14 0.7 s15 0.9 s16 0.8 s17 10 s18 4 s19 2 s20 3 s21 6 s22 9 s23 1 s24 0.9 s25 0.8 s26 0.6 s27 0.3 s28 0.3 s29 0.3 s30 1 s31 3 s32 2 s33 8 s34 2 s35 4 s36 6 s37 3 s38 4 s39 2 s40 5 s41 2 s42 1 s43 7 s44 2 s45 8 s46 1 ; table d2(k,l) distancia de la biorrefiniería al almacén dep biorref 50 ; table c2(k,l) costo de transporte ($*km-1) de la biorrefinería al almacén dep biorref 21 ; table d3(l,m) distancia del almacén a los mercados markCARBON markFOOD markFUEL dep 89 85 81 ; table c3(l,m) costo de transporte ($*km-1) del almacén a los mercados markCARBON markFOOD markFUEL dep 8 9 7 ; VARIABLES ls(j2) lugar de siembra z selección de secciones POSITIVE VARIABLE ls EQUATIONS AreaMax area máxima que se puede sembrar ACult area destinada al cultivo ConsGreat restricción alta ConsLow restricción baja Production1 producción de super sorgo Production2 producción de carbón fijo DepCapacity capacidad de los almacenes DepLimith1 límite de capacidad para almacenes existentes SorghumDist distribución del super sorgo a los mercados BioCarbProd producción neta de pellets de biocarbón TotalCultCt costo total de cultivo DepCost costo de almacenaje TransCost costo de transporte Sales ventas de super sorgo en los tres mercados Profit ganancia neta de la venta de super sorgo a los tres mercados ; AreaMax(i).. AreaMax.up(i)=e=u2(i)*1.5 ; ACult(j).. z=e=sum[(j2), v(j2)*ls(j2)] ; ConsGreat(i).. Consgreat.up(i)=g=sum [(i2), u3(i2)] ; ConsLow(i).. ConsLow.l(i)=l=Areamax.up(i) ; Production1(i).. Production1.up(i)=e=prod[(i2,j2),(ACult.up(j2)*u2(i2))] ; Production2(i).. Production2.up(i)=e=prod [(i2),(Production1.up(i2)*u6(i2))] ; DepCapacity(l).. DepCapacity.l(l)=g=sum [(i2), Production1.up(i2)] ; DepLimith1(l).. DepCapacity.l(l)=l=sum [(i2), u2(i2)] ; SorghumDist(m).. SorghumDist.l(m)=e=prod [(m2,i), B(m2)*u2(i)] ; BioCarbProd(k).. BioCarbProd.l(k)=e=prod [(m2,l2), SorghumDist.up(m2)*u7(l2)] ; TotalCultCt(i).. TotalCultCt.l(i)=e=sum [(i2), [prod[(i3), (Production1.up(i3)*u4(i3))]]+ [prod[(i3), (Production1.up(i3)*u5(i3))]]] ; DepCost(l).. Depcost.l(l)=e=prod [(k2,l2),(BioCarbProd.l(k2)*u8(l2))] ; TransCost(j,k,l,m).. TransCost.l(j,k,l,m)=e=prod[(j2,k2), c1(j2,k2)*d1(j2,k2)]+prod[(k2,l2), c2(k2,l2)*d2(k2,l2)]+prod[(l2,m2), c3(l2,m2)+d3(l2,m2)] ; Sales(m).. Sales.l(m)=e=prod[(m2,i2), (SorghumDist.l(m2)*u9(i2))] ; Profit(m).. Profit.l(m)=e=sum [(i2,j2,k2,l2,m2), (Sales.l(m2)-TransCost.l(j2,k2,l2,m2)-DepCost.l(l2)-TotalCultCt.l(i2))] ; model sorgo /all/ ; solve sorgo using lp maximizing z ; display z.l ;