Optimización multiobjetivo espacial de las prácticas de agricultura de conservación que utilizan un modelo SWAT y un algoritmo evolutivo

Summary

Este trabajo demuestra una integración de un modelo de calidad del agua con un componente de optimización utilizando algoritmos evolutivos para resolver óptimo (de menor costo) la colocación de las prácticas agrícolas para la conservación de un conjunto específico de objetivos de mejora de la calidad del agua. Las soluciones se generan utilizando un enfoque multi-objetivo, lo que permite la cuantificación explícita de compensaciones.

Abstract

Encontrar los costos eficientes (es decir, el más bajo costo) formas de orientar las inversiones en conservación de práctica para el logro de determinados objetivos de calidad del agua en todo el paisaje es de una importancia primordial en la gestión de las cuencas hidrográficas. Métodos de la economía tradicional de encontrar la solución de más bajo costo en el contexto de la cuenca (por ejemplo, ^5,12,20) asumen que el impacto fuera del sitio puede ser descrito con precisión como una proporción de la contaminación genera en el sitio. Estos planteamientos son poco probable que sea representativo del proceso de contaminación real en una cuenca, donde los impactos de las fuentes contaminantes suelen estar determinadas por complejos procesos biofísicos. El uso de los modernos basados ​​en la física, distribuidos espacialmente modelos de simulación hidrológicos permite un mayor grado de realismo en términos de representación del proceso, sino que requiere un desarrollo de un marco de simulación-optimización en el que el modelo se convierte en una parte integral de la optimización.

Evoluciónalgoritmos ary parece ser una herramienta de optimización particularmente útil, capaz de tratar con la naturaleza combinatoria de una cuenca de simulación-optimización problema y permitiendo el uso del modelo de calidad de agua completo. Los algoritmos evolutivos tratar una distribución espacial particular de prácticas de conservación en una cuenca como una solución candidata y utilizar conjuntos (poblaciones) de soluciones candidatas iterativamente aplicando operadores estocásticos de la selección, la recombinación y mutación para encontrar mejoras con respecto a los objetivos de optimización. Los objetivos de optimización en este caso son dispersas para minimizar la contaminación procedente de la cuenca, al mismo tiempo reducir al mínimo el costo de las prácticas de conservación. Un conjunto reciente y creciente de investigación está tratando de utilizar los mismos métodos e integra los modelos de calidad del agua con la definición amplia de métodos de optimización evolutivas ^{3,4,9,10,13-15,17-19,22,23,25.} En esta aplicación, se demuestra un programa que sigue Rabotyagov et al. 'S aproada e integra un moderno y de uso común SWAT modelo de calidad del agua ⁷ con un algoritmo evolutivo multiobjetivo SPEA2 ^26, y el usuario especifica un conjunto de prácticas de conservación y sus costos para buscar las fronteras compensación entre los costos completos de las prácticas de conservación y la calidad especificada por el usuario de agua objetivos. Las fronteras cuantificar las ventajas y desventajas que enfrentan los administradores de cuencas hidrográficas mediante la presentación de la totalidad de los costos relacionados con los objetivos de mejora de la calidad de agua diferentes. El programa permite una selección de configuraciones de cuencas han alcanzado determinadas metas de mejoramiento de calidad de agua y una producción de mapas de colocación optimizada de las prácticas de conservación.

Protocol

1. Preparar Modelo Cuenca y proporcionar datos de entrada para la optimización

1. Crear una base de datos i_SWAT
   1. El uso de un programa llamado "rotor", cree la base de datos desde bases de datos de entrada múltiples, incluyendo los suelos, el clima, la gestión y abono.
   2. Alternativamente, una carrera SWAT existente (posiblemente creado con ArcSWAT o AVSWAT) se pueden importar con i_SWAT.exe. En este caso, el programa "swat_rewrite" se puede utilizar para sustituir la gestión o información URH otro sobre la base de datos a nivel de campo.
   3. Calibración y validación del modelo SWAT se debe realizar en este punto. El SWAT (versión 2005) modelo incorpora dentro de este marco Cuenca del Río Raccoon modelado EA fue calibrado y validado originalmente para un total máximo de carga diaria (TMDL) estudio descrito por Jha et al. (2010). Además de calibración y validación del modelo SWAT se llevó a cabo en apoyo de la elaboración de un Plan Director de Raccoon River Watershed, como descricama en Agren, Inc. (2011), que es el modelo SWAT que se utilizó para este estudio.
   4. Usar una versión modificada de SWAT2005.exe, llamado SWAT2005GA.exe.
2. Prepare el archivo "cuenca presets.csv" - Este es un archivo de texto almacena muchas de las opciones específicas de la cuenca del mapache. Es leído por GeneticiSWAT y MapSWAT a continuación para configurar varios controles y campos de datos para la cuenca con un solo clic.
3. Preparar los costes de los elementos de práctica de conservación. Para este ejemplo, estos se almacenan en la tabla [Práctica Costos - Mapache por condado] en un "costos de Práctica Por subcuenca Josh.mdb" base de datos.

El costo total de una solución candidata es la suma de los costos de las prácticas de conservación aplicadas a las unidades de cuencas ("unidades de respuesta hidrológica", o HRUs). El programa de optimización considera una asignación óptima de una sola práctica de la conservación de un conjunto particular de prácticas de conservación en cada URH tierras de cultivo en la cuenca. La sets de las prácticas de conservación posibles asignados a una Dependencia de Derechos Humanos se denominan grupos de alelos.

4. Crear carpetas SWAT. Para esta carrera, se utilizaron 16 CPU, es decir, 16 copias de SWAT2005GA.exe se llevaron a cabo en 16 carpetas separadas (lo mismo se aplica para sistemas con menos CPU, por ejemplo, 4 carpetas deben ser creadas para un "quad-core" del procesador).

2. Seleccione los parámetros de optimización

1. La optimización se controla mediante el programa llamado "GeneticISWAT". Para realizar la optimización, GeneticISWAT.exe abierto.
   1. Ir a "Archivo", luego "Abrir" y seleccione la base de datos i_SWAT "Raccoon GA.mdb".
   2. Ir a "Archivo" y luego "Configuración" para asignar las rutas de acceso a archivos ejecutables del modelo (SWAT SWAT2005GA.exe).
   3. Ir a "Ejecutar" y luego "Seleccionar conjunto de alelos". Este paso determina las combinaciones de prácticas de conservación utilizados en la optimización. Para esta carrera, el alelo set # 14 fue utilizado, que cuenta con 23 combinaciones de prácticas de conservación. Avda.ailable conjuntos de alelos se almacenan en el archivo de configuración "Alleles.csv".
   4. Ir a "Ejecutar", luego seleccione "SPEA2 subconjunto de referencia Archivo Aware" para realizar la optimización multiobjetivo mediante el algoritmo evolutivo SPEA2.

Figura 1. Establecimiento de objetivos y parámetros de optimización.

Optimización de los parámetros para su selección:

Preset: Seleccione la cuenca para ser optimizado. Al hacer clic en "Aplicar" selecciona las entradas del archivo de presets "cuenca presets.csv" para llenar los valores de control en esta pantalla.

Variable de salida: Seleccione los objetivos medioambientales para la optimización. Como seleccionado (N Outlet Outlet P), esto define una función objetivo 3-dimensional: Nitrógeno (N orgánico + NO3 + NH4 + NO2) con un promedio de 5 años a la salida, Fósforo (Organic P + P Mineral) con un promedio de 5 años a la salida, y el costo total de las prácticas de conservación. Tenga en cuenta que esto creará una frontera disyuntiva 3-dimensional. Las variables de salida alternativos se pueden seleccionar, donde el programa multiobjetivo es minimizar ({} Variable de salida, Costo Total).

Tamaño de la población: Fijar el tamaño inicial de la población. Esto determina el número inicial de soluciones candidatas. Cuando "Semilla con cada alelo" está seleccionada, las soluciones de candidatos que representan a una aplicación uniforme de cada práctica de conservación especificado en el alelo establece para todo el HRUs tierras de cultivo en la cuenca se creó por primera vez. Las soluciones candidatas restantes son creados por una asignación aleatoria de las prácticas de conservación del alelo establecido en tierras de cultivo HRUs. Cuando se selecciona la "Semilla con cada alelo" opción, uno debe asegurarse de que el tamaño de la población inicial es al menos tan grande como el número de alelos de un conjunto de alelo (23 en esta demostración).

Número de generaciones: Ajuste el número deseado de generaciones (iteraciones) para el ciclo de optimización (tenga en cuenta que la carrera se puede reiniciar).

Probabilidad de cruce: Cuando dos soluciones candidatas son seleccionadas para la creación de nuevas soluciones candidatas, probabilidad de cruce indica la probabilidad de que distintas soluciones se crean nuevas (puesto a 1 para esta demostración).

Tamaño de la población temporal: Determina el número de soluciones candidatas nuevos creados. Los recursos del procesador se utilizan más eficientemente cuando este valor es un múltiplo entero del número de hilos de procesador (16 en esta demostración).

Probabilidad de mutación: Especificar la probabilidad de un cambio en la asignación al azar URH a otra práctica de la conservación del conjunto alelo. (Se establece en 0,03 para esta demostración).

Número deTemas: Seleccione el número de procesadores o hilos utilizados. 16 se utiliza en esta demostración.

No. Curve factor de calibración: Se proporciona desde la calibración del modelo SWAT.

Población en Guardar archivo de texto: Es importante seleccionar si se desea reiniciar la optimización correr tras el número especificado de iteraciones se ha completado. Al seleccionar esta opción genera un archivo de texto con los valores de los alelos de cada Dependencia de Derechos Humanos en todas las soluciones candidato sobreviviente (individual). Esto se puede leer de nuevo para reiniciar y continuar la carrera.

Los parámetros secundarios de optimización

Primer Año: Debe establecerse un año después del inicio de la información del tiempo histórico, ya más de 7 años antes del final de estos datos.

Precio del maíz: Se utiliza con la ecuación de pérdida de rendimiento para estimar el costo de la reducción de fertilizantes.

Método de calificación: SPEA2 Archivo. Anotaciones determina la probabilidad de un individuo que sobrevive es para ser seleccionado para cruce.

Método de purga: Dominada. Las personas que son peores en las 3 dimensiones son dominados y purgado.

Fuente HUC: Set de "ubicación especificada", es decir, el valor "7100006" en el siguiente campo "Cuenca del HUC" se utiliza para buscar una fila en la tabla HUC Zone. El valor "07100006" es el código de ocho dígitos HUC para la cuenca del mapache.

Fuente Costo: Set de "Condado (Código URH Vivo)" para indicar que los costos no sean la PCR será determinado por los códigos FIPS de condado en la tabla de costos de la práctica anterior.

Costo Fuente CRP: Ajustado a "1 Lugar" para indicar que el costo CRP será determinado por los códigos FIPS de condado en la tabla de costos de la práctica anterior.

3. Los resultados representativos

GeneticiSWAT.exe produce un archivo de registro que muestra la configuración y los resultados para todas las soluciones de candidatos (personas), así como un "Guardar" de archivos que codifica los resultados de la iteración algoritmo final y que puede ser usado para reiniciar el ciclo de optimización.

En este punto, se puede visualizar todo el conjunto de Pareto-eficientes soluciones (la frontera compensación) siguiendo los siguientes pasos:

1. Ejecutar GeneticiSWAT;
   1. Ir a "Archivo", luego "Abrir" para abrir la base de datos i_SWAT "Raccoon GA.mdb".
   2. Ir a "Ejecutar", luego "Exportar lista Dependencia de Derechos Humanos". Guarde el archivo como "HRU.txt Alelo Raccoon".
2. Producir una animación mediante la ejecución de Mapswat.exe, seleccionando la opción "Ejecutar" y luego en "Animación 3D".

Figura 2. Screenshot para crear "instantáneas" para la visualización frontera 3-dimensional.

La salida es una serie de archivos que se pueden prestar de una vez en archivos de imagen mediante el uso de POV-RAY programa y seleccionando la opción "Render", luego "cola de archivos". Las imágenes se pueden utilizar por sí mismos o combinados en una película que muestra la progresión algoritmo.

Figura 3. Visualización estática de la frontera compensación.

Si se desea, una película que muestra la progresión algoritmo puede ser creado mediante la ejecución de "Framescanner.exe" y siguiendo estos pasos:

1. Ir a "Archivo", luego "Nuevo" y luego "Archivo", luego "Importar" y luego "Archivos PNG". Seleccione las imágenes estáticas.
2. Para crear una película, ir a "Archivo", luego "Exportar", entonces "AVI".
3. Seleccione codec "DIB" para crear archivos AVI a partir de lotes de imagenarchivos.

Cada punto de la frontera representa una configuración de cuencas (una asignación específica de prácticas de conservación de un paisaje). Los mapas de estas configuraciones se pueden ver por toda la frontera, siguiendo estos pasos:

1. Ejecutar Mapswat.exe, seleccione "Ejecutar", luego "Animación Mapa".
2. Seleccione "mapache" en la lista predefinida y haga clic en "Aplicar".
3. Seleccione "Diseño 7 (mapache)" de la lista Map Layout, y luego haga clic en "OK".

Figura 4. Captura de pantalla de la creación de un mapa de cada individuo en la última frontera.

Exportar configuraciones de cuencas específicas (personas) de interés.

A menudo, una cuestión de interés es la selección de configuraciones específicas de las cuencas hidrográficas (individuos) han alcanzado determinados objetivos de calidad del agua. Por ejemplo, se puede desear encontrar un indiviDual en la frontera, lo que reduce en un 30% de nitrógeno y fósforo en un 20% con respecto a las cargas iniciales. MapSWAT permite buscar la frontera para la persona con la menor distancia euclídea para el objetivo especificado. Esto se puede hacer de la siguiente manera:

1. Abra MapSWAT.exe. Seleccione "Execute" | "Buscar".

Figura 5. Capturas de pantalla de búsqueda para un individuo específico en la frontera basado en objetivos de calidad del agua.

2. Introduce mínima y máxima reducción objetivos Tmin y Tmax, así como un tinte intervalo. También se introduce un porcentaje de reducción específica en nitrógeno (nspec) desde el inicio de la "reducción%" espacio al lado de línea de base N, y por ciento de reducción en fósforo (Pspec) en el "% de reducción" junto a la línea base de fósforo. El programa genera una salida en una pantalla emergente:

Figura 6. Captura de pantalla de resultados de la búsqueda

3. Haga clic en "Copiar texto" y pegar en una hoja de cálculo. Tres mesas se producen. En el primero están las personas más cercanas a N y P objetivos del mismo porcentaje de reducción, que va desde Tmin hasta Tmax por Matiz. Justo debajo de este el más cercano solo individuo a la diana (nspec, Pspec) aparece. En segundo lugar, un cuadro en el que los individuos más cercano donde los intervalos de referencia P desde Tmin hasta Tmax mientras que N se mantiene constante cerca nspec. En tercer lugar, un cuadro en el que los individuos más cercanos N objetivos que van desde Tmin hasta Tmax mientras que P se mantiene constante cerca Pspec. En este caso, el más cercano individuo a una reducción del 30% de N fue ID 8423 con un valor de N de 14,637,279.60. Aquí está el mapa que muestra la distribución espacial de las prácticas de conservación y la ubicación de esta configuración hito en la frontera de compensación:

files/ftp_upload/4009/4009fig7.jpg "alt =" Figura 7 "/>
Figura 7. Captura de pantalla de un mapa de ejemplo que describe la persona seleccionada en la frontera. Haga clic aquí para ampliar la cifra .

Exportación de datos de mapas para su posterior análisis es posible siguiendo estos pasos:

1. Ejecutar Mapswat.exe, seleccione "Ejecutar", luego en "Map Data Export".
2. Seleccione "mapache" en la lista predefinida y haga clic en "Aplicar".
3. Introduzca el ID de una configuración de cuenca (individual) (8423 muestra), marque "Mostrar Propiedades de alelos" y "Costos" Mostrar y haga clic en "OK". Estos datos se pueden utilizar para crear mapas personalizados de la configuración cuenca seleccionada utilizando un programa de GIS.

Nombre del Programa	Fuente	Descripción
Rotator	TARJETA
Swat2005GA.exe	USDA Prado, Soil & Water Research Laboratory	Modelo de simulación de cuencas
i_SWAT.exe	TARJETA	SWAT interfaz de base de datos
GeneticISWAT.exe	TARJETA	Algoritmo evolutivo controlador SWAT. Incorpora Galib de http://lancet.mit.edu/ga/ .
MapSWAT.exe	TARJETA	Lee las bases de datos i_SWAT y shapefiles, produce imágenes de las generaciones y las personas.
POV-Ray	Povray.org	Trazador de rayos Persistence of Vision.
Framescanner.exe	Todd Campbell	Imagen PNG a AVI convertidor
Windows Live Movie Maker	Microsoft	Se utiliza para comprimir AVI WMV

Crea y rellena una base de datos con los datos i_SWAT suelo, el clima y la gestión de una cuenca hidrográfica.

Tabla 1. Tabla de programas requerido.

Nombre del archivo	Tipo	Descripción
Raccoon GA.mdb	Acceso a base de datos	Estructura y gestión de las descripciones de las cuencas hidrográficas Raccoon. Leído por GeneticiSWAT y MapSWAT.
cuenca presets.csv	Texto	Presets de ajuste para GeneticiSWAT.exe y MapSWAT
Alleles.csv	Texto	Lista de alelo establece para el algoritmo evolutivo.
HRU.txt Alelo Raccoon	Texto	Archivo creado por GeneticISWAT lista de los alelos determinados como tierras de cultivo. Leído por MapSWAT.
Costos de Práctica Por subcuenca Josh.mdb	Acceso a la base de datos	Costos por prácticas de gestión y del condado.
Terraza Zones.mdb	Acceso a la base de datos	Tabla [HUC datos] tiene la terraza y los números de rendimiento de la zona de la cuenca.
NRI Budgets.mdb	Acceso a la base de datos	Leído por GeneticISWAT.exe para mesas de cultivo y máquina que no se utilizan en esta carrera.
phucrp 2008-12-15.dat	Texto	Calor Planta Unidad de búsqueda en la tabla, no se utiliza en esta carrera.
Management.mdb	Acceso a la base de datos	Rotación tabla de consulta, no se utiliza en esta carrera.
Mapache 28/09/2011 1313.log GA, GA Raccoon 09/29/2011 0732.log, Mapache GA 07/10/2011 0644.log	Texto	Los archivos de registro de carrera GeneticISWAT.
Raccoon GA.wmv	Animación	Visualización 3D de los individuos por generación
Subbasin.shp	ESRI Shapefile	Los esquemas de subcuencas en la cuenca.
Raccoon Map.wmv	Animación	Visualización de los alelos dominantes para cada subcuenca para cada individuo en la frontera.

Tabla 2. Tabla de archivos de muestra requerido.

Discussion

Construimos un sistema integrado de simulación-optimización marco para buscar el sentido de Pareto eficientes conjuntos de configuraciones de las cuencas de menor costo que implican mezcla y ubicación de las prácticas agrícolas de conservación para lograr una serie de objetivos de nivel de cuencas de reducción de nutrientes. Un diagrama conceptual del sistema de simulación-optimización se presentan en la Figura 8. Simulación de cuencas, incluyendo la simulación de los impactos en la calidad del agua de ...