Planificación Ex - Ante Del Uso de La Tierra Linea Base
Planificación Ex - Ante Del Uso de La Tierra Linea Base
Tesis o trabajo de investigación presentada(o) como requisito parcial para optar al título
de:
Magíster en Ciencias Agrarias
Director (a):
[Link]. Luis Miguel Ramírez Nader
Línea de Investigación:
Producción animal tropical
Al profesor Miguel Ramírez Nader, por su permanente apoyo, por su confianza en mis
capacidades y tesón para la culminación de esta investigación.
A Héctor Fabio Espinoza Tilmans, por su colaboración para poder cuadrar los tiempos y
cumplir con mi trabajo y con mis estudios.
A Julio Castro y su familia por su colaboración con los datos de las producciones de la
Hacienda Brasilia y recibirme en su empresa, para poder desarrollar esta tesis.
A la Universidad Nacional sede Palmira, por abrirme de nuevo sus puertas para realizar
los estudios de maestría y de la cual estaré siempre orgullosa y agradecida.
Resumen
La planificación ex-ante empleando la modelación matemática es una herramienta muy útil
que permite obtener una representación lo más cercana a la realidad mediante la
evaluación de distintas variables y predecir los resultados económicos de una empresa.
En el sector agropecuario los empresarios consideran que se presentan diversidad de
objetivos y conflictos. Por esta razón, se realizó un análisis empleando la Programación
Lineal Multiobjetivo en la Hacienda Brasilia ubicada en el municipio del Cerrito, Valle del
Cauca, con el fin de realizar la distribución de un terreno para asignar a cada porción una
actividad agropecuaria de acuerdo a las metas del productor (cultivo de caña de azúcar,
cultivo de uva Isabella, cultivo de maracuyá, y producción ganadera). Para tal fin se
consideraron como funciones objetivo la maximización margen bruto total anual, la
minimización del riesgo económico predial, y la optimización de la mano de obra. Se realizó
recolección de datos con base en los registros de la hacienda desde el año 2010-2016. La
implementación del problema en computador se realizó usando el software Matlab, y se
usó el método de restricciones ϵ modificado para analizar simultáneamente todos los
objetivos. El análisis reveló que no es posible optimizar las funciones objetivo
simultáneamente puesto que estas son contradictorias entre sí. Se encontró que los
cultivos más determinantes son la caña de azúcar y la uva Isabella puesto que, el primero
de ellos tiende a abarcar grandes extensiones de tierra, mientras que el segundo es
intensivo en términos de mano de obra y capital. Posteriores análisis revelaron que, si se
quiere hacer énfasis en margen bruto, debería incrementarse las hectáreas cultivadas de
caña de azúcar y uva Isabella, trayendo consigo esto un incremento en el riesgo predial.
Por otro lado, si se quiere mejorar la utilización de la mano de obra es necesario enfatizarse
en el cultivo de caña de azúcar y la ganadería, trayendo consigo un incremento en el riesgo
económico predial. En cuanto a los puntos de operación históricos del propietario no se
pueden considerar como soluciones fronteras eficientes, sin embargo, los puntos de
operación históricos parecen estar relativamente cerca de las soluciones de frontera
eficientes. El conjunto de soluciones eficientes parece generar un exceso en el consumo
de horas de mano de obra, lo que indica al productor replantear correctamente las horas
X Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Abstract
Ex -ante planning using mathematical modeling is a very useful tool that allows obtaining a
representation as close to reality by evaluating different variables and predicting the
economic results of a company. In the agricultural sector, employers consider that there
are a diversity of objectives and conflicts. For this reason, an analysis was carried out using
Multiobjective Linear Programming at the Brasilia Farm located in the municipality of El
Cerrito, Valle del Cauca, in order to distribute a plot of land to assign an agricultural activity
to each portion according to the producer's goals (sugar cane cultivation, Isabella grape
cultivation, passion fruit cultivation, and livestock production). To this end, the maximization
of the total annual gross margin, the minimization of the real estate economic risk, and the
optimization of the labor force were considered as objective functions. Data collection was
performed based on the records of the farm from the year 2010-2016. The computer
problem was implemented using the Matlab software, and the modified ε constraint method
was used to simultaneously analyze all the targets. The analysis revealed that it is not
possible to optimize the objective functions simultaneously since these are contradictory to
each other. It was found that the most determinants crops are sugar cane and the Isabella
grape, since the first one tends to cover large tracts of land, while the second is intensive
in terms of labor and capital. Subsequent analyzes revealed that if you want to emphasize
gross margin, you should increase the cultivated hectares of sugarcane and Isabella
grapes, bringing with it an increase in property risk. On the other hand, if you want to
improve the use of labor is necessary to emphasize the cultivation of sugar cane and
livestock, bringing with it an increase in property economic risk. As for the owner's historical
operating points, efficient frontier solutions can’t be considered, however, historical
operating points seem to be relatively close to efficient frontier solutions. The set of efficient
solutions seems to generate an excess in the consumption of hours of labor, which
indicates to the producer to correctly rethink the hours of labor available. A critical
parameter for the problem is the capital available in this case study; since a decrease in
available capital can generate solutions that are not feasible by reduce the number of
efficient solutions.
Contenido
Pág.
Lista de figuras
Pág.
Lista de tablas
Pág.
Tabla 2-1. Posibles soluciones de acuerdo a los objetivos .............................................. 42
Tabla 3-1. Condiciones agroecológicas de la Hacienda Brasilia ..................................... 50
Tabla 3-2. Análisis de suelos de la Finca Brasilia............................................................ 51
Tabla 3-3. Actividades productivas de la finca Brasilia .................................................... 52
Tabla 3-4. Margen bruto por tipo de actividad [$/ha-año]. ............................................... 57
Tabla 3-5 Horas hombre por año requeridas para procesar una hectárea según tipo de
cultivo [$/ha-año]............................................................................................................. 58
Tabla 3-6 Otros datos relacionados. ............................................................................... 58
Tabla 3-7 Capital de trabajo requerido para trabajar una hectárea de cada tipo de cultivo
[$/ha-año]........................................................................................................................ 58
Tabla 3-8. Método de pagos para la selección del rango de valores 𝜖. ........................... 67
Tabla 4-1 Solución del problema de optimización para las funciones objetivo
independientes................................................................................................................ 70
Tabla 4-2 Tabla de pagos del problema multiobjetivo. .................................................... 71
Tabla 4-3. Distribución histórica del terreno .................................................................... 75
Tabla 4-4. Una muestra de soluciones eficientes considerando simultáneamente las tres
funciones objetivo. Las variables se enumeran de acuerdo al orden preestablecido de los
cultivos. ........................................................................................................................... 78
Tabla 4-5. Muestra de soluciones eficientes obtenidas al incrementar un 30 % las horas de
mano de obra disponibles al año. ................................................................................... 79
Tabla 4-6. Muestra de soluciones eficientes obtenidas al reducir el capital disponible en un
30 %, es decir, un escenario de escasez de recursos económicos. ................................ 80
Tabla 4-7. Muestra de soluciones eficientes obtenidas al reducir el terreno disponible en
un 50 %, es decir, un escenario de escasez del recurso tierra. ....................................... 81
Contenido XVII
Contenido XVIII
Abreviaturas
Abreviatura Término
ha Hectáreas
PLM Programación lineal multiobjetivo
PL Programación Lineal
SAP Sistemas agropecuarios de producción
Introducción
Los Sistemas Agropecuarios de Producción (SAP), resultan de la compleja interacción de
gran número de componentes interdependientes en los agroecosistemas mixtos, donde el
componente social está presente (Balbino, Porfírio-da-Silva, Kichel, Rosinha, & Costa,
2011; Kichel, Costa, Almeida, & Paulino, 2014).
En Colombia, al igual que sucede en otras regiones de Centroamérica y Latinoamérica los
SAP, se han establecido y funcionan, con deficientes métodos de planificación para el uso
de los recursos (URPA, 2013). Aunque este panorama ha ido cambiando favorablemente,
la falencia en la adecuada planeación persiste, es así como las decisiones relacionadas
con el empleo de las áreas de las tierras y de los tipos de actividades agropecuarias
(cultivos y/o ganadería), no obedecen a estudios científicos y técnicos previos, por el
contrario se toman como referencias, experiencias ajenas o propias, basadas en sus
conocimientos y vivencias, donde la mayoría de las decisiones se toman porque se ha
llegado a ellas por el método de “prueba y error”, a través del tiempo (CDIM-ESAP, 2017).
Para el gobierno nacional, a pesar de que el renglón agropecuario actualmente crece con
mayor rapidez si se le compara con otros sectores productivos, presenta una grave falencia
en la planificación. Se han creado planes decenales en salud y en educación, pero no se
ha hecho ninguno para planificar el agro a largo plazo, por lo tanto es una necesidad que
el sector agropecuario se convierta en uno de los motores de la economía (Fonseca, 2018).
Por esta razón, la planificación de fincas realizada a corto o mediano plazo debe ser muy
detallada, para obtener una mayor rentabilidad; debido a que los objetivos deben ser
coherentes con los del territorio al que pertenecen, con la participación de los interesados
o productores, con un enfoque de mejoramiento de las tecnologías, para adaptarse a las
circunstancias cambiantes de la naturaleza biológica, de lo económico, social, uso del
suelo, vocación agropecuaria, ambiental e institucional del entorno que lo rodea, buscando
la sostenibilidad de los agroecosistemas (Corporacion PBA & ECOPETROL, 2012).
Para la realización de una planificación es ideal que se incluyan aquellos factores que el
productor desea contemplar. Una alternativa posible para lograr lo anterior, es recurrir a
métodos como el de la Programación Lineal Multiobjetivo (PLMO), para conseguir ordenar
las soluciones factibles al problema que se plantea, maximizando y/o minimizando las
funciones objetivo y planteando las restricciones que las condicionan (Alvarado, 2011).
Las restricciones del problema, se refieren a las limitaciones de recursos, que para el
sector agropecuario serian de capital, tierra, mano de obra, agua, mercadeo entre otras.
Al reemplazar las variables que satisfacen las restricciones, se obtiene el conjunto de
soluciones factibles o alcanzables, transformándose en lo posible. Los recursos que
formarán parte de las restricciones deben estar definidos y expresados en términos
cuantitativos (Muñoz Alcaraz, 2016).
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 21
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Teniendo en cuenta lo expuesto, se puede afirmar que con los modelos para la
planificación ex-ante de las actividades empresariales como la agropecuaria, se lograrían
superar debilidades, ganando en precisión, si se analizan con la herramienta matemática
de PLMO, maximizando o minimizando las funciones que debe considerar el productor, al
momento de la toma de decisiones, para lograr alcanzar los objetivos y metas (Zunzunegui,
2017).
Con la anterior metodología, se pueden abordar por parte del sujeto decisor,
simultáneamente varios objetivos, para las decisiones acerca de la planificación de los
sistemas productivos de manera integral u holística, antes de seleccionar y desarrollar en
la realidad las actividades productivas, permitiendo de esta forma predecir los efectos de
éstos sobre la sostenibilidad de los agro ecosistemas, maximizar las utilidades o beneficio
económico, minimizar el uso de recursos como la mano de obra, del agua para riego, del
endeudamiento predial, entre otros (Pacheco & Contreras, 2008).
1. OBJETIVOS
2. MARCO DE REFERENCIA
Entre los “cuellos de botella” del sector agropecuario colombiano, tenemos: envejecimiento
de los habitantes del campo, migración de los jóvenes a centros urbanos, incumplimiento
de la legislación laboral lo cual repercute en salarios por debajo de lo legal, concentración
de la propiedad de la tierra, haciendo que 2/3 partes los productores no la posean, altos
costos de producción y bajos precios en el mercado, falta de financiación oportuna y justa,
falta el servicio de seguros a producción, falta de planificación de producción y
comercialización, falta tener una visión rural más actual acorde a los cambios
internacionales (Lopez, 2014). A pesar de lo anterior, las expectativas de crecimiento para
este sector en Colombia son positivas, y se soportan en: el ascenso en el requerimiento
universal de alimentos que se interrelacionan con las características geo climáticas que
garantizan el incremento permanente de los comodites agropecuarios; la transformación
actual de este gremio, con un crecimiento mayor que otros sectores económicos y el
aumento del aporte del presupuesto del gobierno central a este (FINAGRO, 2014).
de las cuales el 25% se encuentran en 7 países en desarrollo, entre los cuales se cuenta
Colombia (DNP, 2010). De igual forma según la FAO este país se encuentra en el puesto
22 de los 223, de los cuales se evaluó la posibilidad de ampliar la frontera agrícola sin
afectar sus reservas de bosque nativo (SAC y PROEXPORT, 2012) el Ministerio de
Agricultura y Desarrollo Sostenible (MADR), estimo el área de expansión en 16,6 millones
de hectáreas con vocación agrícola y silvo-agricola (DNP, 2011).
“Para el año 2014, Colombia contaba con 5,3 millones de hectáreas cultivadas, con un
potencial estimado de 21,8 millones de hectáreas con vocación de uso agrícola. De este
potencial total, 1,2 millones de ha se encuentran en la Altillanura y de éstas, solo 80.167
Ha (menos del 3%) se encontraban cultivadas. El total de hectáreas con potencial de
aprovechamiento agrícola, pecuario y forestal de la Altillanura se estima en 2,8 millones de
ha, según cálculos de DNP- DDRS, basados en información del IGAC, IDEAM,
INGEOMINAS, MADS, I A v H y CORPOICA” (CONPES, 2014).
haciendo según el potencial de los municipios; la uva tiene uno de los mercados más
dinámicos y debería estar dentro de las prioridades del Valle del Cauca” (El Pais, 2017).
Para el año 2017, se produjeron 2.233.831 toneladas de azúcar, esto como producto de
las 24.380.593 toneladas de caña molida, de esta producción se exportaron 702.764
toneladas de azúcar con destino a Chile, Islas del Caribe, Perú, Estados Unidos, Haití,
26 Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
México y Bolivia el 66% del total. El resto del azúcar se exportó hacia múltiples destinos
alrededor del mundo. Además, se obtuvo una producción de alcohol carburante de 336.753
litros y de melaza de 2017.947 ton, de esta última se exportaron 28.848 ton (ASOCAÑA,
2017).
Sus acciones en este aspecto se dirigen a la conservación del recurso hídrico liderando la
creación del Fondo del Agua por la Vida y la Sostenibilidad, sus actividades se extienden
a 17 cuencas hidrográficas presentes en 26 municipios del Valle del Cauca, restaurando
las cuencas más afectadas ambientalmente y fortalecimiento de las comunidades
asentadas en estos territorios en la parte alta, mediante proyectos productivos (ASOCAÑA,
2017).
toneladas por año (El Tiempo, 2016). El potencial de cultivo de la uva para el Valle esta
cimentado en las características edafoclimaticas, apoyada por el ministerio de agricultura
y el ICA que autorizaron la importación de 100.000 plántulas de uva sin semilla y de la
tecnología que permitirán el buen desempeño de esta especie, las cuales ocuparán unas
75 hectáreas sumando así 1.400 a nivel nacional, la mayoría asentadas en el Valle del
Cauca; con demanda potencial a nivel nacional e internacional, con estos cultivos se
incrementará la demanda de mano de obra y el desarrollo de otras industrias como el
hotelería en torno a estas plantaciones (Prado, 2016).
El Maracuyá es un frutal originario de Brasil, que se adapta desde los 300 metros sobre el
nivel del mar hasta los 1.300, con un rango de temperatura entre 20 – 30 grados
centígrados, es una planta que crece en forma de enredadera, se siembra por semilla en
vivero y se trasplanta cuando tiene el primer zarcillo y una altura de 40 centímetros, la
distancia es de 4 metros entre plantas y 4 entre surcos, para una densidad de 625 plantas
por hectárea, en espaldera de alambre (ICA, 2011).
En cuanto al inventario ganadero lechero según encuesta agropecuaria del Dane 2014, el
total de cabezas 20.944.801 (60.9% hembras y 39.1% machos), al realizar el análisis por
sexos y categorías se encontró que para los machos menores de 12 meses es de 9.5%,
en la categoría de 12 a 23 meses el porcentaje es de 13.7% y de 24 – 36 meses 11.5%;
en cuanto a las hembras el mayor número de cabezas se encuentra en la categoría
mayores de 36 meses con 21.8%, para la edad de 24 a 36 meses es 12.4%, en la edad de
18 a 23 meses 9.1% y en la categoría menores de 12 meses 10.1%. En cuanto al tamaño
del hato el 81% corresponden a menores de 50 cabezas. (DANE – ENA 2014).
Para el año 2014, la producción total de leche fue 17.554.680 litros, con un total de
2.801.063 vacas en ordeño, para un promedio de 6,27 litros vaca por día, procedentes de
22 departamentos, de estos se destacan el departamento de Cundinamarca con promedio
de producción por vaca por día con 15.66 litros, seguido de Antioquia con 11.35 litros,
Nariño 8.65 litros y Boyacá 6.48 litros; respecto al inventario total de cabezas ganado del
país el 6% corresponde a la lechería (DANE & ENA, 2014).
En cuanto el uso del suelo para la producción bovina corresponde al 80% (30.362.366 ha)
del total de producción agropecuaria, 20.336.865 ha en pastos y forrajes, malezas y
rastrojos 7.048.951 ha, 2.974.091 ha en sabanas, paramos y xerofítica (DANE & ENA,
2014).
En cuanto a la generación de empleo directo el sector lácteo contribuye con 717.434, de
los cuales la ganadería doble propósito aporta el 68% y la lechería especializada el 32%
(DANE & ENA, 2014).
La toma de decisiones es una situación que el hombre siempre ha tenido que enfrentar
desde los albores de la humanidad, esta práctica es tan antigua como el ser humano, con
la particularidad que las decisiones a resolver en los diferentes ámbitos de la vida, son
cada vez más complejas y plurales. Este concepto ha recibido innumerables contribuciones
durante su evolución a través de la historia, lo que ha hecho que la teoría de la decisión se
haya consolidado como una disciplina científica (Fernández & Escribano, 2014).
Para resolver un problema de programación lineal existen varios métodos tales como: el
gráfico, algebraico y el simplex con este se pueden resolver problemas con más de tres
ecuaciones simultaneas (Soler, 2013).
32 Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Este último método es un proceso matemático (algoritmo) que utiliza la suma, resta,
división y multiplicación en un singular método continuo para resolver problemas, hoy se
utilizan las computadoras, también para realizar los análisis post-óptimos o de sensibilidad
del modelo analizado (Mateo & Lahoz, 2009).
Proporcionalidad: Supuesto que indica que el cambio del valor de uno de los
símbolos o variables (recursos para producción) de la función objetivo, es
proporcional al valor de la actividad, ambas deben ser lineales. Lo anterior nos lleva
a pensar que los rendimientos son constantes (tanto los crecientes como los
decrecientes), pero en el sector agropecuario esto no es real por tratarse de
procesos biológicos. Por ejemplo: la aplicación de nitrógeno para un cultivo, el
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 33
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Los referentes del uso de la programación lineal en los problemas agrícolas, indican que
los primeros intentos presentaron pocos resultados válidos e iniciaron en la década de los
años 50 (Coronel & Araujo, 2004). El primer intento en este campo, fue realizado por
36 Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Jerome Cornfield en 1941, quien trataba de balancear raciones para los animales a costo
mínimo y llenando los requerimientos nutricionales de estos (Frank, 2001).
Posteriormente con las sucesivas mejoras a este método matemático, el desarrollo de los
sistemas de cómputo y de las rutinas prácticas del cálculo, esta herramienta se convirtió
un instrumento eficaz para planificar una empresa agropecuaria (Alvarado, 2011).
Esta herramienta proporciona métodos útiles y eficientes para la toma de decisiones sobre
problemas que incluyen diversidad de objetivos; los cuales hay que evaluar de manera
conjunta, ya que no es posible obtener la mejor u óptima alternativa, que dará respuesta a
todos los objetivos que el decisor se ha fijado, para la planeación ex – ante de la empresa
agropecuaria (Pliego, 2012).
Como una característica fundamental de las metodologías multicriterio se tiene la
pluralidad de elementos que se pueden combinar en el desarrollo de la evaluación de un
proyecto (Arancibia, Contreras, Mella, Torres, & Villablanca, 2005).
Los problemas que tienen más de un objetivo programado se les aplica la PLMO, para
resolverlos, se debe formular el modelo de tal manera que incluya las variables de decisión
que den respuesta a las restricciones que el sistema productivo genera y que el decisor,
debe lograr optimizar al encontrar los vectores o rango de funciones, que corresponda a
los objetivos del problema a resolver. Además del conjunto de soluciones eficientes, las
funciones objetivo constituyen un espacio multidimensional, llamado “espacio objetivo”,
donde se encuentra el conjunto de soluciones eficientes para resolver el problema
propuesto, de las cuales quien toma las decisiones debe elegir (Pliego, 2012).
Las técnicas de PLMO, o técnicas de optimización vectorial, deben resolver problemas que
optimicen simultáneamente varios objetivos lineales, por lo tanto, no es posible definir un
óptimo para varios objetivos en conflicto. La programación multiobjetivo, trata de encontrar
un conjunto de soluciones eficientes no dominadas llamadas óptimos de Pareto (Reyes &
Oslund, 2014).
Sujeto a: 𝑥 ∈ 𝑆.
Sujeto a 𝑥 ∈ 𝑓.
Dónde:
Optimizar 𝑓1 (𝑥),
Sujeto a 𝑓2 (𝑥) ≥ 𝐿,
𝑥 ∈ 𝑆.
Sin embargo, el método simplex multiobjetivo es más confiable porque entrega una exacta
representación del conjunto eficiente, el inconveniente es la disponibilidad del software
apropiado. “Por ello, la principal ventaja de la PLMO como método de planificación no es
tanto la conducción a un plan esquemático y sencillo, sino la consecución de un método
para analizar una variedad de decisiones alternativas” (Alvarado, 2011).
Meta: se refiere al valor al cual queremos llegar con un atributo o variable, estas
toman la forma de f(x) b, b representa esa cifra propuesta para alcanzar (Garcia,
2017).
Restricciones: son las limitaciones que debe cumplir las soluciones probables del
problema y determinan el conjunto de soluciones viables, ejemplo área máxima a
cultivar de un terreno (Martin, 2017).
Criterios: basado en los objetivos, según su jerarquía, sirve para seleccionar entre
las alternativas y toma de decisiones, pudiéndose contrastar la eficiencia (Martin,
2017).
Para los problemas de optimización se debe tener en cuenta una medida de referencia que
nos indique las bondades de las soluciones encontradas. Para un problema mono-objetivo,
la solución es entregada por la función resuelta, lográndose un conjunto de soluciones
ordenado, del cual se puede elegir fácilmente la alternativa, por ejemplo, para dos
soluciones “a” y “b”, es claro definir si a es mejor que “b” o lo contrario, otra alternativa
seria que fueran iguales.
En la figura 2-1 se observa que “la solución “a” pertenece al conjunto de soluciones
óptimas (o soluciones no dominadas) dado que no puede encontrarse una solución “b” tal
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 41
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
que mejore uno de los objetivos sin empeorar al menos uno de los otros. En este caso, la
solución “c” es dominada por “a” y por “b” “(Aranda Pinilla & Orjuela Castro, 2015).
𝑓𝑗 (𝑥1 ) − 𝑓𝑗 (𝑥2 )
𝑇𝑗𝑘 =
𝑓𝑘 (𝑥1 ) − 𝑓𝑘 (𝑥2 )
300 − 200
𝑇𝐴𝐶 = = 10.
100 − 90
El intercambio indica que cada unidad que disminuye el endeudamiento implica una
disminución de 10 u.m. en el margen bruto. En otras palabras, el costo de
oportunidad de una unidad de endeudamiento es 10 u.m. de margen bruto.
Para 1954 en la reunión anual de economistas agrarios de los Estados unidos, se presentó
un trabajo denominado la programación lineal para no matemáticos un ejemplo práctico
sobre un caso de una empresa agraria en Iowa, por el reconocido economista agrario Earl
O. Heady (1916-1987), con el cual recomendaba que al realizar estudios a gran escala se
debía recurrir a la Programación Lineal (Heady, 1954).
44 Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Swanson (1961), manifestó que para esta época el uso de la programación lineal, era muy
poco en la planificación de las empresas agropecuarias y solo se aplicaba a unas pocas
de gran tamaño, lo contrario ocurrió en el uso en el balanceo de raciones para animales
que ya se había generalizado. Además de la planificación de fincas y el balanceo de
raciones, la programación lineal en el sector agrario ha tenido otras aplicaciones como la
selección del equipo de maquinaria, la evaluación de inversiones, la mezcla de fertilizantes
y de granos, la selección de una cartera de inversiones, la planificación del riego, el trozado
de canales bovinas, la asignación de operarios a puestos de trabajo, el momento óptimo
de corte de un cultivo forestal, etc.(Frank et al., 2014).
Éstos autores publican en la revista Idia del INTA (instituto nacional de tecnología
agropecuaria de Argentina), un artículo sobre la aplicación de la programación lineal a la
empresa agraria de una manera detallada y la acompañaron de una amplia bibliografía,
esto debido a que consideraban que la falta de aplicación de esta ciencia a esta actividad
productiva, se debía a que la información de la que se disponía presentaba conceptos
puramente matemáticos, factor que limitaba su uso por parte de los técnicos y como
ejemplo práctico presentaron datos y resultados de una explotación maicera (Kohout &
Cainelli 1964).
Además de la falta de entrenamiento del técnico, se tenía la no disponibilidad de equipos
de cómputo, aun así, se desarrollaron algunos modelos de programación lineal hasta con
130 actividades agropecuarias y fue hasta la década de los 80 cuando se masifico el uso
de estos equipos razón por la que fue necesario recurrir a métodos manuales conocido
como “programme planning”, desarrollado en Europa (Rossi, 1971).
Para 1981, se incorpora un curso sobre programación lineal aplicada a las empresas
agrarias, para graduados en administración agropecuaria en la Universidad de Buenos
Aires, el cual se continuó impartiendo hasta que en 1991 se incorporó como catedra en las
carreras de administración agropecuaria, esto con la masificación del uso de los equipos
de cómputo (Frank et al., 2014).
Gutiérrez & Buitrago (1974), utilizaron la programación lineal y el computador IBlvl 370/14
con el programa MPSj360, para realizar el cálculo de raciones para cerdos en diferentes
etapas de desarrollo fisiológico, con diferentes materias primas y subproductos de uso
común en Colombia, requerimientos nutricionales (restricciones) de las tablas del NRC y a
menor costo.
lineal, este arrojo múltiple resultados, los vagones más eficientes son aquellos con
capacidad mínima de 3 toneladas, lo que conduce a la minimización de costos,
optimización de los recursos y aumento de la capacidad de molienda, este modelo además
facilita el proceso de toma de decisiones relacionadas con el número de viajes.
3. Materiales y métodos
Por su geología, la hacienda se encuentra en el cuerpo del abanico aluvial, dentro del
paisaje piedemonte y las napas de desborde y de explayamiento de la llanura aluvial de
los ríos tributarios del Cauca (IGAC, 2016). Localizado en la zona plana, influenciada por
la llanura de inundación del rio Cauca y los conos del río Amaime (Figura 3-1).
La Hacienda Brasilia fue adquirida por la familia Castro en la década de 1960, de forma
fraccionada hasta completar 70 has. La distribución inicial de la tierra fue, 56 has (80%)
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 49
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
para ganadería extensiva en producción de leche, y 14 has., (20 % del área total) para la
producción agrícola en cultivos semestrales (maíz).
PARAMETRO INDICADOR
Radiación solar 12 horas luz
Temperatura Promedia: 28-32°C (Mín.18°C-Máx. 32°C)
Precipitación Pluvial 500-600 mm /año
Altitud 994 m.s.n.m.
Relieve plano Pendiente del 2%
Vientos Moderados
Humedad relativa 70-75%
Zona de Vida según clasificación de Bosque muy seco tropical (bms-T).
Holdridge (1947)
2 Periodos de lluvia 6 Marzo a 15 Junio – 16 septiembre a 15 de
Diciembre
2 Periodos secos 16 de Junio – 15 Septiembre- 16 Diciembre a Marzo
15
Clima general Cálido durante el día 28 – 32 grados Cº, y frio en la
noche 18 – 20 ºC
Tomado de Alcaldia de Ginebra (2017)
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 51
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
3.1.4 Suelos
Los suelos que predominan en la zona de localización de la Hacienda Brasilia presentan,
de pH casi neutro a ligeramente ácido con contenidos de carbón orgánico bajo a medio,
medios a altos en calcio y potasio; no son salinos ni sódicos, y, por consiguiente, son suelos
excelentes para la agricultura; dominan los Molisoles y los Vertisoles. Las texturas
predominantes son franco–arenosos o franco–arcillosos sobre gravillosos o limosos.
(Botero, Rodríguez, & Rodríguez, 2006) (Tabla 3-2).
COMPONENTE CANTIDAD
Arenas 51.40 %
Limos 27.00 %
Arcillas 21.60 %
pH 6.4 – 6.6
Materia Orgánica 2.19%
Fosforo Bray II 17.31 ppm.
Cationes intercambiables m- ---
equiv/100: 14. 21
Calcio (Ca) mg/kg. 8. 22
Magnesio (Mg)mg/kg 0. 23
Potasio (K) 0. 12
Sodio (Na) 21.51
CIC 1.73
Ca/Mg 0.44
PSI ---
Sustancias bases
Ca cmol(+) /kg 62.61
Mg cmol(+) /kg 36.04
K cmol(+) /kg 0.87
Na cmol(+) /kg 0.47
Cationes solubles Me/L --
B mg/kg. 0.1 muy bajo
Cu mg/kg. 5.58 normal
Fe mg/kg. 54.72 normal
Mn mg/kg. 12.21 normal
Zn mg/kg. 0.77 bajo
Clasificación del suelo según
Cenicaña 3CO
Fuente: Laboratorio del ingenio providencia, 2014
52 Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Fuente: Castro, 2016 propietario hacienda Brasilia (levantamiento realizado por el topógrafo Marmolejo L).
3.2 Técnicas de recolección y procesamiento de
información para PLMO
Margen bruto por tipo de actividad: calculado a partir de los registros de producción,
ingreso bruto menos costos de producción (tabla 3-4).
Ganadería 579 772 215 658 1 405 738 1 167 407 2 567 797 2 431 180 2 795 857
Horas hombre por año necesarias para trabajar una hectárea según tipo de cultivo:
calculado según costos de la mano obra por hectárea de cada actividad productiva
para la empresa en estudio, dividido por smlv día valor del jornal de 8 horas para
cada año según histórico (Renovación Sindical, 2018) (tabla 3-5).
58 Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Tabla 3-5 Horas hombre por año requeridas para procesar una hectárea según tipo de
cultivo [h/ha-año].
Tipo de
actividad 2010 2011 2012 2013 2014 2015 2016
Caña de azúcar 312 224 220 216 375 224 431
Uva Isabella 3 238 1 609 2 934 4 009 3 691 3 387 3895
Maracuyá 2 740 2 527 2 446 2 447 2 447 2 447 2 814
Ganadería 947 923 977 897 838 963 1 107
Capital de trabajo necesario para trabajar una hectárea y otros datos relacionados:
calculado a partir de mano de obra, insumos, administración, costos de financiación
(cuotas e intereses) (tabla 3-6 y 3-7).
Disponibilidad
de horas de 35 869 28 476 30 583 30 750 44 568 33 101 45 785
mano de obra
[h/año]
Capital
disponible [$] 432 508 373 396 263 458 466 952 659 555 579 393 732 435 211 665 060 367 813 180 308
Tabla 3-7 Capital de trabajo requerido para trabajar una hectárea de cada tipo de cultivo
[$/ha-año].
Caña de
azúcar Uva Isabella Maracuyá Ganadería
El terreno total disponible para realizar las actividades agropecuarias es de 62.6 ha, las
cuales deben ser utilizadas en su totalidad.
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 59
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Variables de decisión
En primera instancia se definió maximizar el margen bruto promedio total. Dicho margen
bruto promedio total puede calcularse de la siguiente manera:
Dónde 𝑀𝐵𝑖 representa el margen bruto medio del tipo de actividad agropecuaria 𝑖.
Se tomó como referencia el modelo de Hazell (1971) y Bocco, Sayago, & Tártara (2002),
asumiendo el riesgo como la sumatoria de las desviaciones absolutas de los márgenes
brutos de cada alternativa productiva, con respecto al promedio del margen bruto de cada
una de ellas para cada año, en el periodo de 7 años de estudios (2010 a 2016). De esta
forma, este objetivo quedó formulado de la siguiente manera:
min 𝑓2 = 𝑥5 + 𝑥6 + 𝑥7 + 𝑥8 + 𝑥9 + 𝑥10 + 𝑥11 + 𝑥12 + 𝑥13 + 𝑥14 + 𝑥15 + 𝑥16 + 𝑥17 + 𝑥18.
Las variables anteriores representan las desviaciones absolutas de los márgenes brutos
(ver también Ignizio (1976)). La relación que debe existir entre dichas variables y las
variables inicialmente planteadas es la siguiente:
𝑀𝐵14 𝑥1 − 𝑀𝐵1 𝑥1 + 𝑀𝐵24 𝑥2 − 𝑀𝐵2 𝑥2 + 𝑀𝐵34 𝑥3 − 𝑀𝐵3 𝑥3 + 𝑀𝐵44 𝑥4 − 𝑀𝐵4 𝑥4 = 𝑥11 − 𝑥12 ,
𝑀𝐵15 𝑥1 − 𝑀𝐵1 𝑥1 + 𝑀𝐵25 𝑥2 − 𝑀𝐵2 𝑥2 + 𝑀𝐵35 𝑥3 − 𝑀𝐵3 𝑥3 + 𝑀𝐵45 𝑥4 − 𝑀𝐵4 𝑥4 = 𝑥13 − 𝑥14 ,
𝑀𝐵16 𝑥1 − 𝑀𝐵1 𝑥1 + 𝑀𝐵26 𝑥2 − 𝑀𝐵2 𝑥2 + 𝑀𝐵36 𝑥3 − 𝑀𝐵3 𝑥3 + 𝑀𝐵46 𝑥4 − 𝑀𝐵4 𝑥4 = 𝑥15 − 𝑥16 ,
𝑀𝐵17 𝑥1 − 𝑀𝐵1 𝑥1 + 𝑀𝐵27 𝑥2 − 𝑀𝐵2 𝑥2 + 𝑀𝐵37 𝑥3 − 𝑀𝐵3 𝑥3 + 𝑀𝐵47 𝑥4 − 𝑀𝐵4 𝑥4 = 𝑥17 − 𝑥18 ,
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 61
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Donde:
𝑀𝐵𝑖𝑗 : es el ingreso bruto promedio generado por la actividad agropecuaria i durante el
periodo j. En este caso se considerarán periodos en años para un total de seis años.
Para abordar este objetivo, fue necesario conocer las exigencias de mano de obra según
el ciclo agrícola, de cada uno de los sistemas de cultivo y el agroecosistemas de ganadería
para producción de leche, durante el periodo de estudio (2010 a 2016). En este sentido, el
tiempo de trabajo se ajustó lo más posible al tiempo disponible al año evitando tiempos
ociosos o uso de tiempo extra no programado. La función objetivo formulada para este
objetivo es la siguiente:
min 𝑓3 = 𝑥19 + 𝑥20 + 𝑥21 + 𝑥22 + 𝑥23 + 𝑥24 + 𝑥25 + 𝑥26 + 𝑥27 + 𝑥28 + 𝑥29 + 𝑥30 + 𝑥31 +
𝑥32 .
Las variables anteriores representan las desviaciones del uso de las horas de mano de
obra con respecto a su disponibilidad durante el año. La relación que debe existir entre
dichas variables y las variables inicialmente planteadas es la siguiente:
Donde:
𝐴𝑖𝑗 : son las horas-hombre requeridas para procesar una hectárea de la actividad
agropecuaria 𝑖 durante el año 𝑗.
𝐻𝑗 : son las horas disponibles en el año 𝑗.
𝑥1 + 𝑥2 + 𝑥3 + 𝑥4 = 𝑆,
Donde:
𝑆 : es el terreno total disponible.
Por otro lado, las restricciones de capital implican que el costo de operar cada actividad
agropecuaria no debe exceder el capital disponible por año. Esto es:
𝐶1 𝑥1 + 𝐶2 𝑥2 + 𝐶3 𝑥3 + 𝐶4 𝑥4 ≤ 𝐾𝑗 ,
Donde
𝐶𝑖 : es el costo de procesar una hectárea del tipo de actividad agropecuaria 𝑖
𝐾𝑗 : es el capital de trabajo disponible al inicio del año 𝑗.
Por último, se tienen restricciones acerca del número de hectáreas mínimos y máximos a
considerar en cada tipo de actividad agropecuaria.
𝑛𝑖 ≤ 𝑥𝑖 ≤ 𝑁𝑖 ,
Donde:
𝑛𝑖 : es el valor mínimo que puede tomar la variable
𝑁𝑖 : es su valor máximo.
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 63
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
min 𝑓3 = 𝑥19 + 𝑥19 + 𝑥20 + 𝑥21 + 𝑥22 + 𝑥23 + 𝑥24 + 𝑥25 + 𝑥26 + 𝑥27 + 𝑥28 + 𝑥29 +
𝑥30 + 𝑥31 + 𝑥32 ,
Sujeto a:
𝑥1 + 𝑥2 + 𝑥3 + 𝑥4 = 𝑆,
𝐶1 𝑥1 + 𝐶2 𝑥2 + 𝐶3 𝑥3 + 𝐶4 𝑥4 ≤ 𝐾1 ,
𝐶1 𝑥1 + 𝐶2 𝑥2 + 𝐶3 𝑥3 + 𝐶4 𝑥4 ≤ 𝐾2 ,
𝐶1 𝑥1 + 𝐶2 𝑥2 + 𝐶3 𝑥3 + 𝐶4 𝑥4 ≤ 𝐾3 ,
𝐶1 𝑥1 + 𝐶2 𝑥2 + 𝐶3 𝑥3 + 𝐶4 𝑥4 ≤ 𝐾4 ,
𝐶1 𝑥1 + 𝐶2 𝑥2 + 𝐶3 𝑥3 + 𝐶4 𝑥4 ≤ 𝐾5 ,
𝐶1 𝑥1 + 𝐶2 𝑥2 + 𝐶3 𝑥3 + 𝐶4 𝑥4 ≤ 𝐾6 ,
𝐶1 𝑥1 + 𝐶2 𝑥2 + 𝐶3 𝑥3 + 𝐶4 𝑥4 ≤ 𝐾7 ,
𝑛1 ≤ 𝑥1 ≤ 𝑁1 ,
𝑛2 ≤ 𝑥2 ≤ 𝑁2 ,
𝑛3 ≤ 𝑥3 ≤ 𝑁3 ,
𝑛4 ≤ 𝑥4 ≤ 𝑁4 ,
𝑥𝑖 ≥ 0.
Es de resaltar, que la función objetivo relacionada con la maximización del margen bruto
medio total se multiplicó por (−1) para obtener la función objetivo 𝑓1 que es una función
objetivo de minimización. Esto se ha realizó con el propósito de que todas las funciones
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 65
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
objetivo que se consideren sean de minimización, sin afectar el resultado final del
problema.
Uno de los métodos más conocidos y utilizados para resolver modelos de optimización
multiobjetivo es el método de restricciones 𝜖, lo que principalmente se debe a su facilidad
de implementación (Jürgen, Kalyanmoy, Kaisa, & Roman, 2008; Matthias, 2005). Es por
eso que en este trabajo se empleó dicho método para resolver el problema de optimización
planteado.
Dónde
𝐱 : representa todas las variables de decisión del problema original, esto es 𝐱 =
{𝑥𝑖 , 𝑛𝑗 , 𝑝𝑗 , 𝑟𝑗 , 𝑠𝑗 }.
𝑆 : representa la región factible generada por el conjunto de restricciones del problema P1.
El método considera que un objetivo es seleccionado para ser optimizado mientras que el
resto de objetivos son convertidos en restricciones, de esta manera se formuló el siguiente
problema denominado P2:
min 𝑓𝑙 (𝐱),
s. a. 𝑓𝑗 (𝐱) ≤ 𝜖𝑗 , ∀𝑗 = 1,2, … , 𝑘, 𝑗 ≠ 𝑙
𝐱 ∈ 𝑆,
En dónde se seleccionó el objetivo 𝑓𝑙 (𝐱), 𝑙 = 1,2, … , 𝑘, para ser optimizado, los objetivos
restantes se restringieron a una cota superior 𝜖𝑗 , 𝑗 = 1, … , 𝑘, 𝑗 ≠ 𝑙, que tomará diferentes
valores para los cuales la solución del problema anterior (P2) estará relacionada con la
frontera eficiente del problema multiobjetivo P1. Note que 𝜖𝑗 representa el peor valor que
puede tomar la función objetivo 𝑓𝑗 (𝐱), en este caso la restricción 𝑓𝑗 (𝐱) ≤ 𝜖𝑗 puede estar o
no activa en la solución óptima.
El rango de valores para cada 𝜖𝑗 está relacionado con los valores mínimo y máximo de la
función objetivo 𝑓𝑗 (𝐱). El método tradicional para obtener dichos valores se basa en una
tabla de pagos, la cual muestra el caso de dos funciones objetivo (Tabla 3.5).
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 67
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
𝑓1 𝑓2
La tabla de pagos no garantiza que los rangos anteriores produzcan puntos en la frontera
eficiente. Con el fin de garantizar soluciones de eficientes Mavrotas (2009) sugiere usar la
tabla de pagos anterior construida con base en el método lexicográfico, un método clásico
de optimización multiobjetivo. Se empleó este método de la siguiente manera:
𝑘
𝑒𝑗
min 𝑓𝑙 (𝐱) + 𝜉 ∑ ,
𝑡𝑗
𝑗=1,𝑗≠𝑙
s. a. 𝑓𝑗 (𝐱) + 𝑒𝑗 = 𝜖𝑗 , ∀𝑗 = 1,2, … , 𝑘, 𝑗 ≠ 𝑙,
𝐱 ∈ 𝑆,
𝑒𝑗 ≥ 0,
Dónde
𝜉 : es un número pequeño seleccionado adecuadamente
𝜉 ∈ [10−3 , 10−6 ], 𝑒𝑗 corresponde a la variable de holgura de la restricción impuesta a la
función objetivo 𝑗, y 𝑡𝑗 es el rango de la 𝑗-ésima función objetivo obtenido de la tabla de
pagos.
4. Resultados
La solución de cada objetivo individual generó los resultados mostrados en la Tabla 4-1.
Para resolver el problema de optimización, considerando el análisis de las funciones
objetivas de forma individual, se seleccionó la de interés y se eliminaron las demás. Para
el caso de la maximización del margen bruto, el cultivo predominante es la caña de azúcar,
seguido del de uva Isabella y el maracuyá. Esto debido a que este cultivo tiene unos costos
de producción bajo, si se compara con las demás producciones en estudio.
Por su parte, para la minimización del riesgo económico predial, la actividad predominante
debe ser la ganadería, seguido del cultivo de la caña de azúcar. Esto se debe a que el
riesgo económico predial es una medida de la variabilidad del margen bruto a lo largo de
los años, la ganadería es la actividad agropecuaria que menos variabilidad tiene. Es decir,
esta actividad genera un rendimiento bruto a lo largo de los años más estable que las otras
actividades, lo que se puede evidenciar de los datos presentados en la metodología. A su
vez la ganadería y el cultivo de caña de azúcar son los que tienen consumos de mano de
obra por hectárea más estables, lo que explica la predominancia de dichas actividades en
la optimización de la mano de obra requerida.
El consumo del capital de trabajo usado comparado con el capital de trabajo disponible
resultó ser mayor cuando se pretende maximizar el margen bruto (ver tabla 4-1). Lo
anterior no se cumple para el resto de funciones objetivo en donde apenas se utiliza
apenas el 50 % en promedio del capital disponible. La razón para esto es que para la
70 Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
maximización del margen bruto el cultivo de la uva Isabella alcanza 10 hectáreas, el mayor
valor obtenido, siendo que este cultivo es intensivo en capital.
Tabla 4-1 Solución del problema de optimización para las funciones objetivo
independientes.
Margen bruto [$/año] 481 093 796 168 994 896 304 878 322
Riesgo económico predial 5 012 441 558 344 803 088 4 950 123 389
[$/año]
* Porcentaje de capital usado promedio calculado con base al capital disponible cada año.
** Porcentaje del tiempo de trabajo usado calculado con base en el tiempo disponible cada
año.
El porcentaje de tiempo de trabajo usado, calculado con base en las horas de trabajo anual
disponibles, en la mayoría de los casos supera el programado. Cuando se pretende
maximizar el margen bruto, las horas de trabajo utilizadas son cerca de un 140% de las
horas disponibles programadas, lo que se debe principalmente al cultivo de uva Isabella
que es intensivo en requerimientos de tiempo de mano de obra, y es uno de los cultivos
predominantes en este caso.
Un resultado similar es obtenido para la minimización del riesgo económico predial, donde
la ganaderia es la actividad de mayor demanda de mano de obra, requiriendo más horas
hombre de las disponibles. Por su parte, para la función optimización de la mano de obra,
las horas hombre trabajadas son similares a las horas hombre disponibles.
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 71
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Para tal fin, se construyó la tabla de pagos cuyo resultado se muestra en la tabla 4-2. Las
funciones objetivo se enumeran siguiendo el orden descrito en la metodología. Es de
mencionar que, dado que se pretende maximizar el margen bruto, esto es equivalente a
minimizar el negativo de dicho margen bruto.
La figura 4-1 muestra la frontera eficiente para el margen bruto y el riesgo económico
predial, los números entre paréntesis representan las hectáreas operadas según cada tipo
72 Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Figura 4-1. Frontera de Pareto para las combinaciones de dos objetivos: margen bruto y
el riesgo económico predial
8
x 10
9
3
-5 -4.5 -4 -3.5 -3 -2.5 -2 -1.5 -1 -0.5
8
f1 x 10
Figura 4-2. Frontera de Pareto para las combinaciones de dos objetivos: margen bruto y
uso de mano de obra
4
x 10
12
(46, 2, 1, 13), f2=5.8109
10
(56, 6, 1, 0), f2=1.01010
8
(55, 7, 1, 0), f2=5.2109
6
f3
Figura 4-3. Frontera de Pareto para las combinaciones de dos objetivos: riesgo predial y
uso de mano de obra
5
x 10
2
(0, 0, 0, 63), f1=3.0108
1.8
1.6
(7, 0, 0, 56), f1=1.9108
1.4
1.2
(14, 0, 0, 49), f1=1.7108
1
f3
0.6
0.2
(46, 2, 1, 13), f1=1.0108
0
3 3.5 4 4.5 5 5.5 6 6.5 7
f2 8
x 10
El propietario del terreno durante los años 2010 a 2016 ha considerado la distribución de
tierras como se muestra en la tabla 4-3. Claramente el propietario ha tendido a favorecer
el cultivo de uva Isabella en los últimos años, mientras ha reducido el cultivo de la caña de
azúcar y ha puesto menor interés a la ganadería. La tendencia del propietario muestra que
la estrategia actual es más intensiva en capital y recursos de mano de obra, y sugiere una
distribución deseada del terreno así: 46 has de caña de azúcar, 10 ha cultivo de uva
Isabela, 1 ha cultivo de maracuyá y 6 ha para actividades ganaderas. Esta estrategia se
analizó a la luz de las soluciones eficientes del problema multiobjetivo.
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 75
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
𝒙𝟒
𝒇𝟐 𝒇𝟑 𝒙𝟏 𝒙𝟐 𝒙𝟑 %
𝒇𝟏 (ha) % horas
Año Riesgo Horas/ha/ (ha) (ha) (ha) capital
Margen bruto Gana usadas
económico año Caña Uva Maracuyá usado
dería
2010 288 954 664 422 646 182 8 284 54 1 1 13 52% 97%
2011 288 954 664 422 646 182 8 284 54 1 1 13 52% 97%
2012 288 954 664 422 646 182 8 284 54 1 1 13 52% 97%
2013 288 954 664 422 646 182 8 284 54 1 1 13 52% 97%
2014 330 200 982 495 076 928 21 565 54 4 1 10 59% 111%
2015 322 512 525 493 531 954 11 888 52 4 1 6 56% 101%
2016 416 014 723 652 078 831 104 390 46 10 1 6 76% 152%
Media 317 792 412 475 896 063 24 425 53 3 1 10 57% 107%
1.8
1.6
1.4
1.2
f3
0.8
0.6
0.4
0.2
5
x 10
1.5
1
f3
0.5
0
-5
-4
-3 9
8
x 10 8
-2 7
6 8
-1 5 x 10
4
f1 0 3 f2
Por su parte, la tabla 4-4 muestra un conjunto de soluciones eficientes para el caso del
análisis simultáneo de las tres funciones objetivo. Esta información es de gran utilidad para
los actores interesados permitiendo hacer una mejor distribución del terreno al tener en
cuenta los tres criterios de decisión.
Por esta razón, es necesario que el propietario tenga en cuenta este aspecto para planear
una mayor disponibilidad de horas hombre al año. Finalmente, una lista más detallada de
soluciones eficientes se encuentra disponible en el Anexo C.
Tabla 4-4. Una muestra de soluciones eficientes considerando simultáneamente las tres
funciones objetivo. Las variables se enumeran de acuerdo al orden preestablecido de los
cultivos.
Solucio 𝒇𝟏 𝒇𝟐 𝒇𝟑 𝒙𝟏 𝒙𝟐 𝒙𝟑 𝒙𝟒 % %
nes ($) ($) (horas/ (ha) (ha) (ha) (ha) capital horas
eficient Margen Riesgo hombre) caña Uva Marac Gana usado usadas
es bruto económico uyá dería
Con el fin de evaluar el efecto que tienen ciertos parámetros en el conjunto de soluciones
eficientes se han considerado las siguientes variaciones:
Este análisis se realizó variando un parámetro a la vez, dejando los demás constantes.
Un incremento en las horas de mano de obra disponibles por año, genera que el porcentaje
de uso de estas sea más razonable, como lo muestra la tabla 4-5.
El impacto que tiene esto es una reducción de las hectáreas cultivadas de caña de azúcar,
y un aumento en el uso del terreno para la ganadería, como se observa en los valores
promedio de las Tablas 4-5 y 4-6, respectivamente. Otro cambio que se observa es una
reducción en el margen bruto y el riesgo económico predial.
Una reducción en el capital de trabajo disponible ocasiona que gran parte de las
formulaciones generadas por el método de restricciones 𝜖 no sean factibles, indicando que
este parámetro es un recurso crítico para la solución del problema. Un efecto de esta
variación es el incremento del porcentaje de uso del mismo, como lo muestra la tabla 4-6.
Esto también ocasiona que el porcentaje de horas de mano de obra usadas se reduzca, lo
que se justifica en el sentido de que actividades menos intensivas en recursos tienden a
ser predominantes, como se observa en una reducción de las hectáreas cultivadas de uva
Isabella al comparar los resultados de las Tablas 4-5 y 4-7.
% %
Soluciones
𝒇𝟏 𝒇𝟐 𝒇𝟑 𝒙𝟏 𝒙𝟐 𝒙𝟑 𝒙𝟒 capital horas
eficientes ($) ($) (horas/ (ha) (ha) (ha) (ha) usado usadas
Margen Riesgo hombre) Caña Uva Mara Gana
bruto medio económico cuyá dería
A 417 864 520 819 474 379 21 241 50 7 1 5 63% 98%
B 374 567 718 572 250 368 67 266 43 6 1 12 61% 101%
C 481 093 796 819 474 379 67 266 52 10 1 0 71% 108%
D 99 832 790 325 026 356 113 290 0 0 0 63 49% 132%
E 374 567 718 572 250 368 113 290 43 6 1 13 61% 101%
F 481 093 796 819 474 379 113 290 52 10 1 0 71% 108%
Media 343 339 706 634 056 370 78 772 36 6 1 20 61% 108%
5. Conclusiones
Se logró formular y validar el modelo de Programación Lineal Multiobjetivo
propuesto. Con él se puede planificar de forma ex ante las áreas de terreno y la
asignación de las producciones, mediante una modelación en la que se pueden
analizar de forma simultanea múltiples objetivos y las restricciones propias del
sistema, la simulación de las combinaciones más eficientes según el tipo de análisis
que se proyecte, ayudando al agente decisor a tomar las decisiones previas a la
ejecución de las actividades productivas.
Para el caso de estudio las soluciones eficientes mostraron las combinaciones para
las tres funciones de forma simultánea (maximización del margen bruto promedio,
minimización del riesgo económico predial y la optimización de la mano de obra),
independientemente de las condiciones edafoclimaticas de la finca, estas le
servirán al propietario como guía, puesto que se lograron mediante el análisis de
los datos históricos de sus producciones durante los 7 años estudiados.
El análisis demostró que las actividades más determinantes para lograr un mejor
margen bruto son la caña de azúcar y la uva Isabella, siendo el primero el que
ocupa mayor área y el segundo el de mayores costos de producción. Para lograr
mayor optimización de la mano de obra las actividades deben ser caña de azúcar
y ganadería.
Para realizar este tipo de estudio los mayores limitantes que se tienen para aplicar
los modelos de programación lineal multiobjetivo, son la falta de información en la
empresa agropecuaria y de técnicos agropecuarios con los conocimientos
matemáticos y de Tics, al igual que de equipos interdisciplinarios que permitan
abordar simultáneamente mayor información, para la planeación ex – ante de forma
holística y sostenible.
84 Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
6. Recomendaciones
La aplicación de este modelo multiobjetivo, permite realizar una planificación del
uso de la tierra, de los recursos disponibles y sus restricciones para que los
empresarios agropecuarios que tengan la información requerida para su aplicación,
puedan tomar las decisiones más acertadas sin que se requiera haber ejecutado
sus proyectos con los riesgos y consecuencias que se corren en el sector
agropecuario, logrando así mejores rendimientos y una adecuada planificación de
sus producciones, teniendo como base sus históricos de producción.
Otra aplicación debería ser para la planificación de los entes gubernamentales que
se encargan de realizar la planificación del sector agropecuario de nuestro país.
Planificación ex - ante del uso de la tierra con programación multiobjetivo, en un 85
agroecosistema de cultivos y ganadería bovina en el Valle del Cauca, Colombia
Anexo A: Formulación compacta del
problema de optimización
Conjuntos
Parámetros
𝑚𝑖𝑗 Margen bruto [$/ha] generado por la actividad agropecuaria 𝑖 obtenido en el año 𝑗,
𝑖 ∈ 𝐶, 𝑗 = 1,2, … , 𝑇.
𝑚̅𝑖 Margen bruto medio [$/ha] generado por la actividad agropecuaria 𝑖 en el año 𝑗,
̅ 𝑖 = (1/𝑇) ∑𝑇𝑗=1 𝑚𝑖𝑗 .
𝑖 ∈ 𝐶, 𝑗 = 1,2, … , 𝑇. Note que 𝑚
𝑎𝑖𝑗 Horas hombre por hectárea por año [h/ha-año] requeridas para operar la actividad
agropecuaria 𝑖 en el año 𝑗, 𝑖 ∈ 𝐶, 𝑗 = 1,2, … , 𝑇.
Variables de decisión
max ∑𝑚
̅ 𝑖 𝑥𝑖 .
𝑖∈𝐶
Se tomará como referencia el modelo de Hazell (1971) y Bocco, Sayago, y Tártaro (2002),
asumiendo el riesgo como la sumatoria de las desviaciones absolutas de los márgenes
brutos de cada alternativa productiva, con respecto al promedio del margen bruto de cada
una de ellas para cada año, en el periodo de 6 años de estudios (2010 a 2016). De esta
forma, este objetivo queda formulado de la siguiente manera:
𝑇
min ∑ |∑ 𝑚𝑖𝑗 𝑥𝑖 − 𝑚
̅ 𝑖 𝑥𝑖 | .
𝑗=1 𝑖∈𝐶
Observe que la función objetivo anterior es una función no lineal; sin embargo, se conoce
una formulación lineal equivalente para dicha función. Para esto sean 𝑛𝑗 y 𝑝𝑗 , 𝑗 = 1,2, … , 𝑇,
dos nuevas variables de decisión (ver también Ignizio (1976)). La formulación equivalente
requiere también la adición de un nuevo conjunto de restricciones como sigue:
𝑇
min ∑ 𝑛𝑗 + 𝑝𝑗 ,
𝑗=1
s. a. ∑ 𝑚𝑖𝑗 𝑥𝑖 − 𝑚
̅ 𝑖 𝑥𝑖 = 𝑛𝑗 − 𝑝𝑗 , ∀𝑗 = 1,2, … , 𝑇.
𝑖∈𝐶
Bibliografía 89
Para abordar este objetivo, es necesario conocer las exigencias de mano de obra según
ciclo agrícola, de cada uno de los sistemas de cultivo y el agroecosistema de ganadería
para producción de leche, durante el periodo de estudio (2010 a 2016). En este sentido,
se pretende que el tiempo de trabajo se ajuste lo más posible al tiempo disponible al año
evitando tiempos ociosos o uso de tiempo extra no programado. La función objetivo
quedaría formulada como sigue
𝑇
min ∑ |∑ 𝑎𝑖𝑗 𝑥𝑖 − ℎ𝑗 | .
𝑗=1 𝑖∈𝐶
Al igual que la función objetivo anterior, esta función es no lineal, su transformación a una
función lineal equivalente se realiza añadiendo las variables adicionales 𝑟𝑗 y 𝑠𝑗 , 𝑗 = 1,2, … , 𝑇,
y un conjunto de restricciones adicionales. En este sentido el problema modificado es:
𝑇
min ∑ 𝑟𝑗 + 𝑠𝑗 ,
𝑗=1
s. a. ∑ 𝑎𝑖𝑗 𝑥𝑖 − ℎ𝑗 = 𝑟𝑗 − 𝑠𝑗 , ∀𝑗 = 1,2, … , 𝑇.
𝑖∈𝐶
∑ 𝑥𝑖 = 𝑆.
𝑖∈𝐶
Por otro lado, las restricciones de capital implican que el costo de operar cada actividad
agropecuaria no debe exceder el capital disponible por año. Esto es:
∑ 𝑐𝑖 𝑥𝑖 ≤ 𝐾𝑗 , ∀𝑗 = 1,2, … , 𝑇.
𝑖∈𝐶
90 Planificación ex-ante con análisis multiobjetivo, para optimizar el uso de la tierra
en un agroecosistema agropecuario en Cerrito, Valle del Cauca, Colombia
Por último, se tienen restricciones acerca del número de hectáreas mínimos a considerar
en cada tipo de actividad agropecuaria.
𝑥𝑖 ≥ 𝑁𝑖 , ∀𝑖 ∈ 𝐶.
min 𝑓1 = − ∑ 𝑚
̅ 𝑖 𝑥𝑖 ,
𝑖∈𝐶
𝑇
min 𝑓2 = ∑ 𝑛𝑗 + 𝑝𝑗 ,
𝑗=1
𝑇
min 𝑓3 = ∑ 𝑟𝑗 + 𝑠𝑗 ,
𝑗=1
s. a. ∑ 𝑚𝑖𝑗 𝑥𝑖 − 𝑚
̅ 𝑖 𝑥𝑖 = 𝑛𝑗 − 𝑝𝑗 , ∀𝑗 = 1,2, … , 𝑇,
𝑖∈𝐶
∑ 𝑎𝑖𝑗 𝑥𝑖 − ℎ𝑗 = 𝑟𝑗 − 𝑠𝑗 , ∀𝑗 = 1,2, … , 𝑇,
𝑖∈𝐶
∑ 𝑐𝑖 𝑥𝑖 ≤ 𝐾𝑗 , ∀𝑗 = 1,2, … , 𝑇,
𝑖∈𝐶
∑ 𝑥𝑖 = 𝑆,
𝑖∈𝐶
𝑥𝑖 ≥ 𝑁𝑖 , ∀𝑖 ∈ 𝐶,
𝑥𝑖 , 𝑛𝑗 , 𝑝𝑗 , 𝑟𝑗 , 𝑠𝑗 ≥ 0, ∀𝑖 ∈ 𝐶, 𝑗 = 1,2, … , 𝑇.
Bibliografía 91
clear all
close all
clc
load Parameters
LB = zeros(size(F,1),1);
k = size(F,2); % Number of objectives
n = 3; % Number of Pareto points
options = optimoptions('linprog', 'display','off'); % Ajuste del
algoritmo de programación lineal
A = cat( 2, A, zeros(size(A,1),1) );
LB = [LB; 0];
for i = 1:n
for j = 1:n
AeqNew = cat(2, Aeq, zeros(size(Aeq,1),1));
AeqNew = cat(1, AeqNew, [F(1:end-1,2:end)', ones(k-1,1)]);
beqNew = cat(1,beq,[eps(i,1);eps(j,2)]);
[x, ~, exitflag] = linprog(F(:,1), A, b, AeqNew, beqNew, LB, [],
[],options);
end
end
figure;
mesh(fval1, fval2, fval3)
xlabel('$f_1$','Interpreter','latex')
ylabel('$f_2$','Interpreter','latex')
zlabel('$f_3$','Interpreter','latex')
Bibliografía 93
Bibliografía
Alcaldia de Ginebra, V. del C. (2017). Sitio Oficial del Municipio de Ginebra, Valle del
Cauca, Colombia. Recuperado 25 de enero de 2017, a partir de [Link]
[Link]/informacion_general.shtml
Anghinoni, i., moraes, a., carvalho, p. C. F., souza, e. D., conte, o., & lang, c. R. (2012).
Beneficios da integração lavoura-pecuária sobre a fertilidade do solo em sistema
plantio direto. Recuperado 23 de mayo de 2018, a partir de
[Link]
integracao-lavoura-pecuaria-sobre-a-fertilidade-do-solo-em-sistema-plantio-direto
Arancibia, S., Contreras, E., Mella, S., Torres, P., & Villablanca, I. (2005). Evaluación
Multicriterio: aplicación para la formulación de proyectos de infraestructura deportiva.
Santiago de Chile. Recuperado a partir de
[Link]
Aranda Pinilla, J. A., & Orjuela Castro, J. A. (2015). Optimización multiobjetivo en la gestión
de cadenas de suministro de biocombustibles TT - Una revisión de la literatura
Multiobjective optimization in biofuel supply chain management. A review of the
literature. Ingeniería, 20(1), 21-47. Recuperado a partir de
[Link]
750X2015000100004&lang=pt
Bibliografía 97
Balbino, L. C., Porfírio-da-Silva, V., Kichel, A. N., Rosinha, R. O., & Costa, J. A. (2011).
Documentos 303: Manual Orientador para Implantação de Unidades de Referência
Tecnológica de Integração Lavoura-Pecuária- Floresta – URT iLPF. Embrapa
Cerrados. Documentos 303. Planaltina-DF.
[Link]
Banco Agrario de Colombia. (2015, noviembre). Colombia: Perspectivas macro del sector
agropecuario en 2015 – 2016. Investigaciones Economicas, 4. Recuperado a partir
de
[Link]
[Link]
Bocco, M., Sayago, S., & Tártara, E. (2002). Modelos multicriterio: una aplicación a la
selección de alternativas productivas. Agricultura Técnica, 62(3), 450-462.
[Link]
Botero, P. J., Rodríguez, J., & Rodríguez, C. (2006). Paisajes y territorio ancestral de
Palmira Valle del Cauca,. En Territorio ancestral, rituales funerarios y chamanismo en
Palmira prehispánica. Departamento de Antropología. Recuperado a partir de
[Link]/1552/2/01P
Coronel, M., & Araujo, P. (2004). La Programación Lineal Aplicada Al Manejo Forestal,
Economía y Administración forestal, Facultad de ciencias Forestales, Universidad
Nacional de Santiago del Estero, Argentina. Recuperado a partir de
[Link]
[Link]
DANE, & ENA. (2014). Boletin Tecnico: Encuesta nacional agropecuaria. Bogota,
Colombia: DANE. Recuperado a partir de
[Link]
na_2014.pdf
DNP. (2011). Bases del Plan nacional de desarrollo 2010-2014. Bogota, Colombia:
Departamento Nacional de Planeacion. Recuperado a partir de
[Link] PND 2010-2014 Versión 5 14-04-
2011 [Link]
Eiselt, H. A., & Sandblom, C. L. (2007). Linear programming and its applications. Springer
Science & Business Media. Recuperado a partir de
[Link]
ming_and_its_Applications.pdf
Eisgruber, L. M., & Reisch, E. (1961). A note on the application of linear programming by
agricultural economics departments of land grant colleges. Journal of Farm
Economics, 43(2), 303-307.
El Pais. (2017, mayo 7). Estos diez cultivos podrían transformar al Valle del Cauca.
Periodico El Pais noticias de Economia, p. 1. Recuperado a partir de
[Link]
[Link]
Estupiñan, F. (2017, enero 7). Más sombras que luces para el agro en el 2017. El
Espectador, p. 1. Recuperado a partir de
[Link]
2017-articulo-673662
Fernández, B. G., & Escribano, C. M. (2014). Decision Theory: state of the art after three
hundred years. Boletín de Estadística e Investigación de Operaciones, 30(3), 292-
312. Recuperado a partir de [Link]
/files/BEIOVol30Num3_HistoriayEnsenanza.pdf
Frank, R, G., George, B., & Dantzig. (2014). La programación lineal y el sector agrario,
Reconocimiento de la Sociedad Argentina a los Referentes del Cambio realizada a
través del programa “Argentina, Sociedad Anónima”. Recuperado a partir de
[Link]
Dantzig, la programación [Link]?sequence=1%0A
Getán, J., & Boj, E. (2014). Programación lineal. Barcelona: Universitat de Barcelona.
Recuperado a partir de
[Link] [Link]
Gobernación de Antioquia. (2014). Manual técnico del cultivo de maracuyá bajo buenas
practicas agricolas (1.a ed.). Medellín, Colombia: Francisco Velez. Recuperado a partir
de [Link] BPA_0.pdf
Gutiérrez, H., & Buitrago, J. (1974). Cálculo de raciones de mínimo costo para cerdos en
zonas tropicales, CIAT Cali Colombia. Recuperado a partir de
10 Planificación ex-ante con análisis multiobjetivo, para optimizar el uso de la tierra
2 en un agroecosistema agropecuario en Cerrito, Valle del Cauca, Colombia
[Link]
Guzman, J. (2016, abril 12). Colombia duplicaria su produccion de uva y el valle del cauca
es el líder. Agronegocios, p. 1. Recuperado a partir de
[Link]
y-el-valle-del-cauca-es-el-lider-2621888
Hausmann, R., Rodrik, D., & Velasco, A. (2005). Growth Diagnostics. Estados Unidos: John
F. Jennedy School of Government Harvard University. Recuperado a partir de
[Link]
pdf
ICA. (2012). Manejo fitosanitario del cultivo de la vid (Vitis vinifera y [Link]): Medidas
para la temporada invernal. Bogota, Colombia: Produmedios. Recuperado a partir de
[Link]
57b7984c5afb/[Link]
Jürgen, B., Kalyanmoy, D., Kaisa, M., & Roman, S. (2008). Multiobjective optimization.
Berlin Heidelberg:Springer.
Kichel, A. N., Costa, J. A. A., Almeida, R. G., & Paulino, V. T. (2014). Sistema de integração
lavoura-pecuária-floresta (ILPP) - Experiencia no Brasil. Boletim de Indústria Animal.
[Link]
Leibovich, J., & Estrada, L. (2009). Competitividad del sector agropecuario colombiano
Capitulo IV. Colombia. Recuperado a partir de [Link]
content/uploads/2017/05/[Link]
Lopez, C. (2014). Los 10 cuellos de botella del sector agropecuario | Finanzas | Economía
| Portafolio. Recuperado 24 de mayo de 2018, a partir de
[Link]
60292
MADR. (2017). Implementacion politica para mejorar la competitividad del sector lacteo
nacional. Bogota, Colombia: Ministerio de agricultura y desarrollo rural, Ministerio de
Comercio. Recuperado a partir de
[Link]
Manuel/INFORME UE FASE 2 -TERCER TRAMO VARIABLE %28JUNIO 7 2017%29
%281%[Link]
Maino, M., Pittet, J., & Kobrich, C. (1993). Programación Multicriterio: Un Instrumento Para
El Diseño De Sistemas De Producción. Red Internaciona de Metodologia de
investigacion de sistemas de produccion. Santiago de Chile. Recuperado a partir de
[Link]
[Link]/bitstream/handle/10625/14853/[Link]?sequence=1&isAllo
wed=y
investigacion-operativa-en-ingenieria/Contenido/Materiales/t8.1toma-de-decisiones
Martínez, A. L. (2009). NRC vs. INRA to formulate rations based on dry roughages and
pelleted concentrates for leisure horses. Archivos de zootecnia, 58(223), 333-344.
Recuperado a partir de [Link]
67650489390&partnerID=40&md5=7fe9ae530b4b73c5b8984b9c199a4871
Mateo, P. M., & Lahoz, D. (2009). Programación Lineal III. Análisis Post-optimal. Zaragoza:
Programación Lineal III. Recuperado a partir de
[Link]
operativa/ficheros/[Link]
MMA, & A. (2014). Guia para la evaluacion ex ante de proyectos de riego (1.a ed.). Bolivia:
Ministerio de Medio Ambiente y Agua. Recuperado a partir de
[Link]
f?sequence=1&isAllowed=y
Muñoz, R. ., Ramón Nonato, B., & Albelo, M. (2005). Modelo de programación lineal para
optimizar alternativas agroproductivas en fincas del piedemonte amazónico,
Universidad de la Amazonia, Colombia. Revista Momentos de Ciencia, 3(1), 2006.
Recuperado a partir de [Link]
Núñez, M., & Aspitia, M. (2013). Manual para Desarrollar Capacidades Institucionales en
la Gestión del Riesgo Agroempresarial (IICA- Instituto Interamericano de Cooperación
para la Agricultura). San Jose, Costa Rica: AGRIS. Recuperado a partir de
[Link]
[Link]
[Link]?sequence=1
Praguer, M. M., Restrepo, J. M., Angel, D. I., Malagon, R. M., & Zamorano, A. M. (2002).
Agroecologia: Una disciplina para el estudio y desarrollo de sistemas sostenibles de
producción agropecuaria (1.a ed.). Palmira: DIPAL. Recuperado a partir de
[Link]
Ramos, A., Sánchez, P., Ferrer, J. M., Barquín, J., & Linares, P. (2010). Modelos
matemáticos de optimización (Publicación ténica). Universidad Pontificia Comillas.
[Link]
Reyes, O., & Oslund, F. (2014). Teoría del bienestar y el optimo de pareto como problemas
microeconómicos. Revista Electrónica de Investigación en Ciencias Económicas,
2(3), 217-234. [Link]
Reyes, O., & Oslund, F. (2014). Teoría del bienestar y el óptimo de pareto como problemas
microeconómicos. Revista Electrónica de Investigación en Ciencias Económicas,
2(3), 18. Recuperado a partir de [Link]
[Link]
Seibane, C., & Riachi, J. (2014). Rol profesional de los ingenieros agrónomos y forestales.
En Introducción A Las Ciencias Agrarias Y Forestales, Una Primera Aproximación A
La Realidad Teoría Y Metodología Para Una Mirada Problematizadora Y Crítica De
Los Distintos Componentes De La Realidad Agropecuaria Y Forestal (1.a ed., p. 354).
Recuperado a partir de
[Link]
quence=3
URPA. (2013). Proyecto: estudios y planificación del uso del suelo rural, con criterios
sociales, económicos y ambientales, Unidad De Planificación De Tierras Rurales,
Adecuación De Tierras Y Usos Agropecuarios. Recuperado a partir de
[Link]/..uso.../39499890-aa61-47b5-b2c2-98c1dccd3978
[Link]?sequence=1