Un sistema ultrasónico de bajo costo para la localización de equipo oceanográfico en aguas someras

Resumen

En el presente trabajo se describen las características más importantes del primer prototipo de un sistema ultrasónico capaz de detectar frecuencias de 26 kHz en profundidades ligeramente mayores de 30 m. El sistema surge de la necesidad de apoyar los trabajos de localización y recuperación de equipo oceanográfico que se fondea sin boya de señalización superficial en aguas someras. Los elementos que componen el sistema son: un circuito electrónico de transmisión y recepción de señales acústicas, un transductor acústico piezoeléctrico y una carcasa resistente a la presión. La electrónica del transceptor contiene un algoritmo que calcula la distancia entre un par de transceptores con base en el tiempo de propagación de la señal ultrasónica. La carcasa está hecha de PVC, un material de bajo costo y fácil de mecanizar. Los resultados de las pruebas operativas indican un error cuadrático medio de ± 0.4 m para la distancia entre transceptores, que es un orden de magnitud menor al error de precisión de ubicación en los sistemas de posicionamiento global. El sistema ultrasónico constituye una herramienta de apoyo idónea para disminuir el riesgo de pérdida de equipo oceanográfico por eventos hidrometeorológicos o actividades pesqueras en cuerpos de agua someros.

Abstract

This paper describes the most important features of the first prototype of an ultrasonic system capable of detecting frequencies of 26 kHz in depths slightly greater than 30 m. The system arises from the need to support the fieldwork of locating and recovering oceanographic equipment moored without surface buoy in shallow waters. The elements of the system are: an electronic circuit for transmitting and receiving acoustic signals, a piezoelectric acoustic transducer and a pressure resistant housing. The electronics in the transceiver contains an algorithm that calculates the distance between a pair of transceivers based on the propagation time of the ultrasonic signal. Housing is made of PVC, a material inexpensive and easy to machine. The results of operational tests indicate a root mean square of ± 0.4 m for the distance between transceivers, being an order of magnitude less than the precision error location in global positioning systems. The ultrasonic system is an ideal support tool to reduce the risk of loss of oceanographic equipment caused by hydrometeorological events or fishing activities in shallow water bodies.

Introducción

La ingeniería electrónica es fundamental para obtener mediciones confiables de las corrientes, las olas y las propiedades del agua en los cuerpos marinos, salobres y de agua dulce. Las aplicaciones son amplias, pues incluyen el desarrollo de tecnologías para la utilización de fuentes de energía renovable oceánicas, el diseño y construcción de instrumentos para estudios del manejo integral y desarrollo sustentable de los ecosistemas costeros y para la búsqueda de evidencia de contaminación por hidrocarburos en aguas profundas a través de vehículos autónomos submarinos, así como para la dispersión de contaminantes en aguas superficiales, entre otros. Una gran diversidad de instrumentos y carcasas para proteger los sistemas embebidos se han diseñado de acuerdo con el fenómeno o variable a medir, al método de medición, y también de acuerdo con las condiciones ambientales adversas a resistir (e.g. Chereskin y Howe, 2007; Ulloa, 2007; Vera, 2007; Joseph, 2014).

La instrumentación electrónica, como los correntómetros y las sondas multiparámetros por mencionar algunos, puede dividirse de acuerdo con el método de medición en (Chereskin y Howe, 2007): puntual (anclaje, lance, sondeo), lagrangiana (boya de deriva, dispersión de colorante) y remota (acústica, electromagnética). De particular interés son los instrumentos acústicos que se fondean en un punto fijo cerca de las ciudades costeras, bocanas de lagunas y las desembocaduras de los ríos, pues suelen presentarse dificultades durante las labores de localización/recuperación, entre las cuales se encuentra la turbiedad del agua por el efecto de las corrientes y el oleaje. Si la visibilidad es insuficiente, puede ocurrir que el buzo no localice el instrumento o que le lleve demasiado tiempo encontrarlo. Además de colocar sedimentos en suspensión que pueden cubrir el equipo de medición, las corrientes y el oleaje poseen la capacidad de mover el equipo fuera de su posición original, lo que complica el trabajo de localización/recuperación. Una solución consiste en señalizar la posición de fondeo mediante una boya superficial o subsuperficial (e.g. Trask y Weller, 2001), sin embargo, resulta inconveniente porque el instrumento queda a merced de la curiosidad de los pescadores locales, así como del arrastre de las redes que se utilizan durante la captura de camarón (e.g. INP, 2000). Otra solución es colocar en el instrumento un mecanismo de liberación acústico (Joseph, 2014), de tal modo, que al estar unido a una boya subsuperficial, pueda ascender por flotabilidad y recuperarse en la superficie del agua. No obstante, dicha solución es impráctica, pues se pierde la plataforma de trípode donde se monta el instrumento (anclaje), además de que en caso de enterramiento es probable que la boya no proporcione flotabilidad suficiente para recuperarlos. Aunque la ubicación del instrumento se determina por medio de un sistema de posicionamiento global, es mejor contar con un dispositivo auxiliar que apoye la labor de localización/recuperación del anclaje en profundidades del orden de 30 m. El objetivo de este trabajo consistió en diseñar y construir el prototipo de un sistema ultrasónico de bajo costo para localizar instrumentos fondeados en aguas someras.

Desarrollo

Diseño del sistema

El sistema ultrasónico que se diseñó y construyó consta de un par de transceptores. Básicamente, un transceptor consta de un transductor acústico piezoeléctrico (e.g. Kleeman y Kuc, 2008), un circuito electrónico, las baterías y la carcasa. El circuito electrónico se compone de un microcontrolador Atmel AVR, un transmisor, un receptor y la fuente de alimentación (figura 1). En el microcontrolador se generan señales ultrasónicas con una frecuencia de (26 ± 1) kHz mediante un algoritmo (Hernández, 2015). Adicionalmente, se agregó un relevador mecánico con bobina de 5 V a fin de conmutar el sensor piezoeléctrico con el transmisor de potencia durante el tiempo que se transmite la señal ultrasónica, así como al receptor durante el tiempo de espera. En este trabajo se empleó una batería recargable sellada de ácido-plomo, de 12 V con capacidad de 1.2 A/h para la alimentación de los circuitos. La tabla 1 resume las características del circuito electrónico.

Figura 1.

Descripción ilustrada del circuito electrónico del transceptor acústico

Tabla 1.

Especificaciones del circuito electrónico

Descripción	Especificación
Tamaño del circuito	10 x 8 cm
Voltaje a la salida del regulador de voltaje	4.97 V
Corriente eléctrica consumida sin transmisión	0.002 A
Corriente eléctrica consumida con transmisión	0.117 A
Duración de la batería con transmisión	10 h
Duración de la batería sin transmisión	25 d

Uno de los transceptores, denominado el transceptor submarino (TS), se fondea junto con el anclaje donde permanece en un modo pasivo continuo (modo de espera). El otro transceptor, denominado el transceptor de rastreo (TR), se utiliza para estimar la distancia a la que se encuentra el transceptor submarino, lo cual implica ubicarse cerca de la posición del anclaje a bordo de una embarcación menor. Para la comunicación entre los transceptores TS y TR, se desarrolló el siguiente protocolo de aplicación para la transmisión-recepción. El TR emite una señal ultrasónica que se forma por una ráfaga de pulsos con una duración de 1.02 s y espera 9.18 s para recibir una señal de respuesta por parte del TS. En cuanto el TS recibe y detecta los pulsos a la frecuencia de emisión, se espera un tiempo muerto de 2 s para que lleguen todos los pulsos que transmitió el TR. Luego, el TS emite una primera señal de respuesta al TR que consiste en una ráfaga de pulsos con una duración de 1.02 s y espera 7.18 s, la diferencia entre 9.18 s y 2 s, para después emitir una segunda serie de pulsos y de ahí pasar nuevamente al modo de espera. Cuando el TR recibe la primera señal de respuesta, espera un tiempo de 2 s para que lleguen todos los pulsos del TS y queda en modo de espera de la segunda ráfaga de pulsos del TS. Este proceso se repite indefinidamente.

El circuito del TR mide el tiempo desde que se emite la señal ultrasónica hasta que recibe un par de señales de respuesta por parte del TS, el cual se denomina aquí tiempo medido (t_m). El tiempo que se almacena en la memoria del microcontrolador corresponde a la suma de un ciclo del TR (10.20 s) en tres intervalos de viaje de la señal ultrasónica entre el TR y el TS, y el (t_m). Existen dos tiempos que no se contabilizan en el microcontrolador, a saber, el tiempo de espera del TR de 2 s, así como un tiempo muerto de 0.25 s que corresponde a una medida cualitativa del tiempo de respuesta de los circuitos electrónicos inherentes en el microcontrolador (Hernández, 2015). Entonces, la distancia entre transceptores (d) se estima como

La velocidad del sonido en el océano, c, se considera constante (1500 m/s).

Circuito electrónico

La figura 2 muestra el diagrama electrónico del circuito, correspondiente a los bloques en la figura 1. El circuito electrónico contiene el microcontrolador AVR de 8 bits, ATmega32, que cuenta con las siguientes características: alto rendimiento y bajo consumo de energía, 1 kb de EEPROM, un convertidor A/D de 10 bit, así como arquitectura avanzada RISC con 131 instrucciones en lenguaje ensamblador. Además, opera en el intervalo 4.5-5.5 V, su consumo de energía es de 1.1 mA cuando está activo y 0.35 mA en modo inactivo. El microcontrolador opera con un reloj interno de 8 MHz. Los requerimientos que la electrónica debe satisfacer los impone el fabricante. Con el objeto de generar la onda acústica que se aplica al transmisor, se implementaron un conjunto de instrucciones en el entorno del programa Bascom AVR®. La ejecución de las instrucciones genera una señal cuadrada (TTL) en el pin PB0 del microcontrolador.

Figura 2.

Diagrama electrónico simplificado del circuito de la Figura 1

El receptor se compone de dos etapas: amplificador de señal y acondicionador de señal. En la primera etapa, la señal que recibe el sensor receptor se amplifica 1000 veces en dos pasos, mediante un amplificador por 100 y un amplificador por 10, tal como se observa en la figura 2. Cada paso se refiere a la relación establecida por las resistencias de retroalimentación mediante configuración inversora. Como estos dos pasos son en cascada, ambas ganancias se multiplican y de ahí una ganancia de 1000. En estas etapas de amplificación se utilizó un circuito integrado LM358; sus características de acuerdo con el fabricante son: alta ganancia en voltaje, 10⁵, frecuencia compensada en temperatura, frecuencia máxima de 1.1 MHz, operación con fuente simple y con voltaje de alimentación 3-32 V dc, bajo offset 2 mV, temperatura (0-70) °C y encapsulados DIP8, SO8 Mini SO8 y TSSOP8. Hasta este punto, la forma de onda de la señal recibida (TS) se muestra en la figura 3a.

Figura 3.

a) onda que recibe el TS, b) onda que transmite el TR

Designando la alta ganancia en voltaje del circuito LM358 como A₀ = 105, donde β es el factor de realimentación establecido por las resistencias y de acuerdo con el diagrama del circuito electrónico de la figura 2, se tiene que en el primer paso de amplificación, la ganancia finita del circuito se da por, v_o / v_i= (β+1/A_o)^–1= 100.89, donde v_o y v_i son los voltajes de salida y de entrada al circuito amplificador, respectivamente, y β = 101^–1. El error relativo debido a esta ganancia finita del amplificador operacional es 1/1000 o 0.1%. El error relativo correspondiente al segundo paso de amplificación es 0.01%.

En la segunda etapa, la señal amplificada se acondiciona en un circuito que consiste de un diodo y un comparador (figura 2). El acondicionador de señal cumple el rol de convertir toda la ráfaga de pulsos que se recibe en un pulso único. Este pulso único es el que recibe la entrada de interrupción del microcontrolador. El criterio para verificar el tren de pulsos acústicos es que los conos de transmisión/recepción se encuentren alineados lo más posible, pues así los transceptores pueden sincronizarse para calcular la distancia entre ellos. En cuerpos de agua someros, la profundidad mínima en la cual puede operar una embarcación menor es de unos 0.5 m, en tanto que la altura de la plataforma trípode donde se coloca equipo de medición acústico es aproximadamente de 0.5 m. Una profundidad de fondeo segura para evitar que la embarcación colisione con el equipo sería de unos 2 m. Es por ello que se considera que la aplicación práctica del sistema ultrasónico en el campo no requiere algoritmo anticolisión para equipo de medición fondeado en profundidades mayores de 2 m.

El transmisor de potencia consta de un circuito integrado TDA2003 que viene en un encapsulado llamado Pentawatt. Dicho encapsulado está diseñado para funcionar con una alimentación simple desde 8 V hasta 18 V. Los componentes que se requieren para su buen funcionamiento se detallan en la hoja que proporciona el fabricante. La figura 3b ilustra la forma de la onda que se transmite (TR). El voltaje máximo de salida, que corresponde a la señal acústica que se transmite, es de 9 V pico.

Carcasa

La presión es la variable oceanográfica que determina la geometría de las carcasas de la instrumentación electrónica, por lo general cilindros o esferas (Chereskin y Howe, 2007; Joseph, 2014; MacDonald y Mullarney, 2015). Las carcasas protegen las partes internas de los transceptores (los transductores y toda la electrónica) cuando se colocan bajo el mar o en cuerpos de agua dulceacuícolas. Los materiales que componen las carcasas deben resistir o protegerse de la corrosión marina. En profundidades menores de 600 m, es frecuente utilizar materiales plásticos (Chereskin y Howe, 2007), por ejemplo, PVC (Joseph, 2014; Murphy et al., 2015).

La figura 4 muestra la carcasa cilíndrica que se diseñó para alojar los transductores. Se eligió PVC como material de construcción por su bajo costo, disponibilidad y facilidad de mecanizado en comparación con titanio o delrin®. La carcasa se compone de un tapón roscado hembra cedula 40 (6 mm de espesor) y de diámetro nominal de 100 mm, un adaptador espiga roscado macho cedula 40 (6 mm de espesor) y de diámetro nominal de 100 mm, un tramo de tubo de 100 mm de diámetro cuya longitud puede variar según se requiera, y una tapa ciega. El diseño lo completa una rueda de acrílico de 6 mm de espesor que constituye la ventana acústica, y sobre la cual se fija el transductor piezoeléctrico. La carcasa puede resistir presiones hasta de 413 686.6 Pa (60 psi o 4 atm) antes que la rueda de acrílico se fracture de acuerdo con las pruebas de estanqueidad realizadas por Hernández (2015). La figura 5 muestra una fotografía del transceptor submarino. El TS se monta sobre la plataforma de trípode que sujeta también al equipo de medición. La plataforma se baja al fondo desde una embarcación mediante un cabo y con la ayuda de buzos.

Figura 4.

Descripción ilustrada de la carcasa, 1) tapón roscado hembra, 2) adaptador espiga roscado macho, 3) tramo de tubo, 4) tapa ciega, 5) rueda de acrílico

Figura 5.

El transceptor submarino con pesos sujetos a la carcasa mediante cinchos a fin de hundirlo en el agua y realizar pruebas de campo. El transceptor de rastreo tiene la misma forma

Resultados y discusión

Se realizaron dos pruebas del funcionamiento del sistema ultrasónico en Playa Miramar, frente a la costa de Ciudad Madero, Tamaulipas (22°16’ 20.40” N, 97°44’ 4.46” O; 22°16’ 55.61” N, 97°42’ 0.07” O), a una profundidad de 20 y 33 m, respectivamente. En estos dos sitios se midió el perfil de la velocidad del sonido con un equipo denominado CastAway™-CTD (YSI, 2010). Cada prueba consistió en colocar el TR sobre la superficie del agua a bordo de una embarcación menor con motor fuera de borda anclada al fondo y orientada en dirección de la corriente promedio, para luego bajar gradualmente el TS en profundidades conocidas. En cada profundidad se efectuaron cinco mediciones y se estimó la distancia al TR. El promedio de estas se muestra en las figuras 6a y 7a con sus respectivas desviaciones estándar, y los valores correspondientes en la tabla 2.

Tabla 2.

Distancias medidas entre el TR y el TS a intervalos de profundidad de 1 m

Profundidad (m)	Distancia (m)	Desviación estándar (m)	Distancia (m)	Desviación estándar (m)
0.5	0.6	0.4	0.7	0.0
1.5	2.0	0.5	1.5	0.6
2.5	2.6	0.2	2.8	1.2
3.5	3.4	1.3	3.8	0.7
4.5	4.2	0.5	4.4	0.3
5.5	4.7	0.3	5.9	0.8
6.5	5.8	1.1	6.7	0.7
7.5	7.7	1.0	7.6	0.8
8.5	8.5	1.3	8.9	0.1
9.5	9.0	0.3	9.6	0.3
10.5	10.2	1.6	10.9	0.2
11.5	11.3	0.7	12.1	1.0
12.5	12.0	0.8	12.2	0.4
13.5	13.6	1.5	13.7	0.5
14.5	14.1	1.4	14.8	0.5
15.5	15.8	0.2	15.7	0.4
16.5	16.4	1.1	17.2	2.0
17.5	17.8	0.6	18.2	1.9
18.5	18.3	0.5	18.6	0.1
19.5	18.6	0.9	19.4	0.6
20.5			21.6	2.7
21.5 22.5			21.7 22.5	0.8 0.6
23.5			23.5	0.0
24.5			24.7	0.8
25.5 26.5 27.5 28.5 29.5 30.5 31.5 32.5			26.0 26.3 27.9 28.7 29.5 29.9 31.8 32.6	0.6 0.7 0.1 0.1 0.9 1.9 0.6 0.1

Las variaciones en las desviaciones estándar (figuras 6a y 7a) sugieren movimientos verticales a causa del oleaje y el efecto advectivo del flujo superficial promedio, así como variaciones en el alineamiento de los conos de transmisión entre el TS y el TR dentro del área de emisión ultrasónica. En ausencia de olas, dicha área es del orden de 300 m² a una profundidad de 30 m. Las ubicaciones de los puntos estuvieron razonablemente lejos de la influencia de la corriente de descarga del río Pánuco.

Figura 6.

a) Distancias entre transceptores en función de la profundidad, b) perfil de la velocidad del sonido. Mediciones realizadas en la posición 22°16’ 20.40” de latitud norte y 97°44’ 4.46” de longitud oeste

Figura 7.

a) Distancias entre transceptores en función de la profundidad, b) perfil de la velocidad del sonido. Mediciones realizadas en la posición 22°16’ 55.61” de latitud norte y 97°42’ 0.07” de longitud oeste

En los primeros 3 m de profundidad, el perfil de la velocidad del sonido muestra la influencia dominante de la variación vertical de la salinidad sobre la variación correspondiente a la temperatura (figuras 6b y 7b). En el resto de la columna de agua, los cambios que se observan en la velocidad del sonido se deben principalmente a las variaciones verticales de la temperatura. La velocidad del sonido en la columna de agua (figuras 6b y 7b) fue mayor al valor nominal del océano y, por ello, se utilizaron valores promedio locales para el cálculo de la distancia entre transceptores. Las velocidades del sonido promedio a 20 m y 33 m son en orden respectivo, 1530.5 m/s y 1530.0 m/s, con una variación medida en el campo de ± 4.7 m/s y ± 4.1 m/s, respectivamente.

El valor del error cuadrático medio de las distancias entre transceptores, evaluado como: fue de ± 0.40 m en una columna de agua con un tirante de aproximadamente 20 m (figura 6a) y de ± 0.37 m en una columna de agua de unos 33 m de profundidad (figura 7a). Estos errores son inherentes a la estimación de las profundidades, donde se situó el TS, pero son un orden de magnitud menor al error de precisión de ubicación en los sistemas de posicionamiento global con tecnología del Sistema de Aumentación de Área Amplia que es del orden de ± (1−3)m (Madry, 2015).

Aunque en el cálculo de la distancia entre transceptores se consideró la velocidad del sonido constante, en realidad es una función de la temperatura, la salinidad, la presión estática y sus perfiles dependen de la localidad, la estación del año, del clima y los procesos de mezcla que se producen en la capa superficial del mar por la acción del viento y el oleaje (Waite, 2005). La propagación del sonido depende a su vez de las fronteras (superficie, fondo) y los dispersores acústicos presentes en la columna de agua, por ejemplo, burbujas de aire, sedimentos en suspensión y organismos (Medwin, 2005). Entonces, como los transceptores requieren tener previamente asignado un valor promedio constante de la velocidad del sonido, la localización de equipo oceanográfico en cualquier época del año dependerá de contar con una buena estimación local de la velocidad del sonido. Lo anterior se puede lograr con la caracterización anual de la velocidad del sonido en la zona de estudio donde se instalen anclajes. Por otra parte, la inhomogeneidad de la columna de agua produce una pérdida exponencial de la señal ultrasónica durante las etapas de transmisión y recepción, ocasionando a su vez que la amplitud de la señal sea altamente variable. Por ejemplo, para la señal acústica de 9 V pico que transmite el TR, el TS únicamente recibe voltajes del orden de cientos de mV. Dicha pérdida no puede evitarse y es la razón por la cual se necesita, en el diseño de los transceptores, un circuito amplificador de señales.

En áreas costeras donde la actividad pesquera es activa, el equipo oceanográfico generalmente se fondea sin boyas superficiales que indiquen su ubicación. El prototipo del sistema ultrasónico constituye una buena herramienta de apoyo para la recuperación del equipo, pues permite posicionar la embarcación prácticamente sobre el anclaje de manera que el buzo lo encuentre con relativa rapidez, minimizando así el riesgo de pérdida cuando se haya movido fuera de su posición original a causa de la dinámica costera o de las actividades propias de la pesca de arrastre. El costo de los transceptores es aproximadamente 40% menor al costo de los sistemas comerciales, los cuales no son de manufactura nacional. El sistema actual tiene la ventaja de que el circuito receptor puede modificarse fácilmente para medir el nivel de ruido submarino, tema de interés en la evaluación de propuestas para la construcción de infraestructura en habitat marinos (Williams et al., 2005; Farcas et al., 2016). En términos de mejoras al prototipo, se requiere incrementar el rendimiento de la batería con el objeto que el sistema permanezca montado en anclajes oceanográficos someros por un periodo mayor de 30 días, así como disminuir el error cuadrático medio de las distancias entre transceptores. La carcasa, aunque funcional, no es estética por lo que una mejora al diseño y menores dimensiones es parte de un trabajo futuro.

Conclusiones

La realización de mediciones en el medio marino requiere un equipo con un alto grado de sofisticación tecnológica y por ende su adquisición resulta onerosa. Los gastos de importación representan en general, un porcentaje importante del costo del equipo. El riesgo de pérdida, cuando el equipo opera en aguas costeras durante la evolución de eventos hidrometeorológicos o en aquellas regiones donde las actividades pesqueras son importantes, es alta. En consecuencia, el mantenimiento y cuidado del equipo de medición es una actividad primordial en toda actividad de investigación limnológica y oceanográfica.

El presente trabajo contribuye con el diseño y construcción de una herramienta auxiliar que facilita las labores de localización y recuperación de equipo oceanográfico fondeado en aguas someras, con que se pretende disminuir el riesgo de pérdida. La herramienta auxiliar es el primer prototipo de un par de transceptores capaces de detectar frecuencias de 26 kHz hasta profundidades ligeramente mayores de 30 m a través de un protocolo de aplicación de diseño propio para la transmisión-recepción de señales ultrasónicas. El costo de los transceptores es 40% menor al costo de equipo comercial de manufactura extranjera. Las pruebas operativas efectuadas al sistema ultrasónico demostraron que es posible calcular la distancia entre ambos transceptores con un error del orden de ± 0.4 m, un buen logro respecto al error de precisión de ubicación de los GPS comerciales. Este logro se debe en gran medida a no utilizar la velocidad del sonido promedio del océano. La velocidad del sonido en las zonas estuarinas y áreas costeras adyacentes es variable, caracterizarla, requiere al menos un año de mediciones con una frecuencia mensual en toda el área de interés. La obtención de mediciones confiables es por tanto cara, peligrosa dependiendo de las condiciones de la mar y lenta en términos del tiempo total, pero relativamente rápida en cuanto a la duración de un muestreo mensual. A pesar del largo tiempo que precisa obtener resultados en el medio marino, el trabajo de campo y el procesamiento de datos son de naturaleza práctica, es la manera correcta de trabajar en un sustrato que está en movimiento continuo.

En México, no existe una industria dedicada al diseño y fabricación de equipo oceanográfico aún cuando se tiene la necesidad de generar tecnología propia para la explotación de hidrocarburos en mar profundo, energía eléctrica mediante fuentes de energía oceánica, así como para atender las necesidades de la industria en el monitoreo de vertimientos en áreas fluviales y costeras. Por ello es importante impulsar los desarrollos tecnológicos. Este trabajo avanzó un paso pequeño para el logro de tal fin.

Agradecimientos

El apoyo financiero para la realización del presente trabajo fue gracias a los proyectos SIP-20120616 y CONACyT-FOMIX Tamaulipas 175692. Los comentarios de dos revisores anónimos ayudaron a enriquecer el manuscrito original.