13.07.2015 Views

Árbol de regresión para determinar el potencial de licuación: ARELI

Árbol de regresión para determinar el potencial de licuación: ARELI

Árbol de regresión para determinar el potencial de licuación: ARELI

SHOW MORE
SHOW LESS
  • No tags were found...

Create successful ePaper yourself

Turn your PDF publications into a flip-book with our unique Google optimized e-Paper software.

Árbol <strong>de</strong> regresión <strong>para</strong> <strong>de</strong>terminar <strong>el</strong> <strong>potencial</strong> <strong>de</strong>licuación: <strong>ARELI</strong>S. R. García, Migu<strong>el</strong> P. Romo O. & Efraín Ovando S.Instituto <strong>de</strong> Ingeniería, UNAM, MéxicoABSTRACTThis paper presents a new method for assessing the liquefaction potential of soils based on geotechnical and seismic<strong>de</strong>scriptions. Data from in situ and laboratory testing are used to <strong>de</strong>v<strong>el</strong>op a Cognitive Analysis and Computingenvironment for predicting the occurrence and non occurrence of liquefaction. No empirical equation is used toapproximate the unknown limit state curve 2D, a Regression Tree is proposed to map the soil and seismic <strong>para</strong>meters(translated and transfigured by Fuzzy Logic) to nominal and numeric conclusions (calculated by Neural Networks)about liquefaction.RÉSUMÉEn este documento se presenta un nuevo método <strong>para</strong> <strong>de</strong>terminar <strong>el</strong> <strong>potencial</strong> <strong>de</strong> licuación <strong>de</strong> los su<strong>el</strong>os basados en<strong>de</strong>scripciones geotécnicas y sísmicas. Se explotan los datos <strong>de</strong> pruebas <strong>de</strong> campo y laboratorio <strong>para</strong> construir unambiente <strong>de</strong> Cómputo y Análisis Cognitivo que permite pre<strong>de</strong>cir la ocurrencia o no ocurrencia <strong>de</strong> la licuación. Se<strong>el</strong>iminan las tradicionales ecuaciones empíricas <strong>para</strong> <strong>de</strong>finir la curva <strong>de</strong> estado límite en un espacio 2D y se incluye unÁrbol <strong>de</strong> Regresión que permite mapear las entradas <strong>de</strong>l material y la carga cíclica (traducidas y transfiguradas conLógica Difusa) hacia conclusiones nominales (calculadas con Re<strong>de</strong>s Neuronales) y numéricas <strong>de</strong> la licuación.1 INTRODUCCIÓNEl po<strong>de</strong>r <strong>de</strong>structivo <strong>de</strong>l fenómeno <strong>de</strong> licuación llamó laatención <strong>de</strong> la ingeniería mundial a partir <strong>de</strong> los dosgran<strong>de</strong>s terremotos <strong>de</strong> 1964 (Alaska, 27 <strong>de</strong> marzo,Mw= 9.2 y Niigata, 16 <strong>de</strong> junio, Mw= 7.5) don<strong>de</strong> severos<strong>de</strong>splazamientos laterales y fallas <strong>de</strong> flujo ocasionaronenormes pérdidas económicas y humanas. Lainvestigación <strong>de</strong>rivada <strong>de</strong> estos eventos se dirigió hacia<strong>el</strong> esclarecimiento <strong>de</strong> los mecanismos que controlan lalicuación, al pronóstico <strong>de</strong> su ocurrencia y la mitigación<strong>de</strong> sus efectos <strong>de</strong>vastadores.La experiencia en campo y laboratorio indica que lalicuación y la <strong>de</strong>formación <strong>de</strong>l terreno asociada a estefenómeno constituyen un problema complejo (difícil <strong>de</strong>mo<strong>de</strong>lar física y analíticamente) por lo que losprocedimientos empíricos son la práctica estándar en laevaluación <strong>de</strong> la resistencia a la licuación, predicción <strong>de</strong>la <strong>de</strong>formación <strong>de</strong>l su<strong>el</strong>o y <strong>el</strong> diseño <strong>de</strong> acciones <strong>de</strong>rehabilitación.Numerosos son los procedimientos <strong>de</strong>sarrollados<strong>para</strong> establecer <strong>el</strong> <strong>potencial</strong> <strong>de</strong> licuación <strong>de</strong> los su<strong>el</strong>osgranulares saturados sin embargo, la <strong>el</strong>ección <strong>de</strong> laecuación empírica a<strong>de</strong>cuada no es tarea fácil. La altaincertidumbre en la <strong>de</strong>finición <strong>de</strong> los ambientes sísmicosy en la caracterización <strong>de</strong> los su<strong>el</strong>os obliga a labúsqueda <strong>de</strong> mo<strong>de</strong>los simples constituidos por doscomponentes esenciales: 1) un marco analítico basadoen la experiencia (casos históricos), y 2) un índice quer<strong>el</strong>acione parámetros <strong>de</strong> pruebas in situ y lascaracterísticas <strong>de</strong> licuación <strong>de</strong>l su<strong>el</strong>o. Lasaproximaciones <strong>de</strong>rivadas <strong>de</strong> estas condiciones buscanvincular <strong>el</strong> mínimo número <strong>de</strong> parámetros <strong>de</strong> los su<strong>el</strong>osy <strong>de</strong> los sismos con una interpretación simplista <strong>de</strong>lfenómeno.Las conexiones <strong>de</strong>rivadas <strong>de</strong>l comportamiento en lanaturaleza, con un gran número <strong>de</strong> variables (conocidasy ocultas) moviéndose en rangos dinámicos(<strong>de</strong>pendientes <strong>de</strong> la variable tiempo) difíciles <strong>de</strong> <strong>de</strong>finir,hacen <strong>de</strong>l fenómeno <strong>de</strong> licuación un candidato no apto<strong>para</strong> ser estudiado con reglas lineales o <strong>de</strong> bajadimensión.El campo <strong>de</strong>l cómputo cognitivo ha mostrado ser unaalternativa po<strong>de</strong>rosa y eficiente <strong>para</strong> mo<strong>de</strong>lar problemasno-lineales multidimensionales que se <strong>de</strong>sarrollen enambientes con incertidumbre y contaminación. En estetrabajo se examina un nuevo procedimiento semiempírico<strong>para</strong> evaluar <strong>el</strong> <strong>potencial</strong> <strong>de</strong> licuación <strong>de</strong> lossu<strong>el</strong>os saturados sujetos a carga cíclica.Esta propuesta combina la lógica difusa y las re<strong>de</strong>sneuronales <strong>para</strong> alimentar e interpretar los resultados <strong>de</strong>un árbol <strong>de</strong> regresión, <strong>ARELI</strong>. Con <strong>ARELI</strong> se extraen lasr<strong>el</strong>aciones naturales <strong>de</strong>trás <strong>de</strong> las observaciones <strong>de</strong>campo y las condiciones conocidas <strong>de</strong>l su<strong>el</strong>o. Elaprendizaje neuronal y la expresión difusa se usan <strong>para</strong>la interpretación numérica y lingüística <strong>de</strong> losparámetros <strong>de</strong> entrada y las recomendaciones <strong>de</strong> salidacon las que <strong>ARELI</strong> construye <strong>el</strong> discurso “<strong>potencial</strong> <strong>de</strong>licuación”. Un aspecto esencial en este nuevo mo<strong>de</strong>lo esque, siendo extremadamente simple, generaestimaciones <strong>de</strong> gran exactitud. Adicional, pero nomenos importante, es la ventaja <strong>de</strong> realizar unase<strong>para</strong>ción natural <strong>de</strong>l efecto <strong>de</strong> las propieda<strong>de</strong>s <strong>de</strong>su<strong>el</strong>o, <strong>de</strong> los métodos <strong>de</strong> prueba, <strong>de</strong> la carga sísmica y<strong>de</strong> otros aspectos que, <strong>de</strong>bido a las fronteras y sesgos<strong>de</strong> los mo<strong>de</strong>los tradicionales, no es posible cuantificarcorrectamente.


Con base en los resultados <strong>de</strong> esta investigación, seconcluye que <strong>ARELI</strong> es una atractiva alternativa <strong>para</strong><strong>de</strong>scribir <strong>el</strong> fenómeno <strong>de</strong> licuación a pesar <strong>de</strong> lasobservaciones limitadas, <strong>el</strong> conocimiento parcial <strong>de</strong> lascondiciones geotécnicas y sísmicas, así como <strong>de</strong>lrazonamiento ambiguo con <strong>el</strong> que se infiere <strong>el</strong> universo<strong>de</strong> comportamientos <strong>de</strong> los su<strong>el</strong>os.2 REVISIÓN DE LAS FORMULACIONES PARAEVALUAR EL POTENCIAL DE LICUACIÓNLas variables con las que se expresa <strong>el</strong> <strong>potencial</strong> <strong>de</strong>licuación son 1) la <strong>de</strong>manda cíclica que afecta a unestrato <strong>de</strong> su<strong>el</strong>o, expresada en términos <strong>de</strong> la r<strong>el</strong>ación <strong>de</strong>esfuerzo cíclico CSR (por sus siglas en inglés cyclicstress ratio) y 2) la capacidad <strong>de</strong>l su<strong>el</strong>o <strong>para</strong> resistir lalicuación, expresada en términos <strong>de</strong> la r<strong>el</strong>ación <strong>de</strong>resistencia cíclica CRR (por sus siglas en inglés cyclicresistance ratio). El procedimiento original <strong>de</strong> Seed eIdriss (1971), <strong>para</strong> calcular <strong>el</strong> CSR, permanece como laguía básica <strong>para</strong> expresar la <strong>de</strong>manda sísmica. En loque respecta al CRR, a pesar <strong>de</strong>l gran número <strong>de</strong>refinamientos en los componentes <strong>de</strong>l marco analítico,los mayores progresos se han generado en laestandarización e interpretación <strong>de</strong> las pruebas índice insitu (prueba <strong>de</strong> penetración estándar SPT, prueba <strong>de</strong>penetración <strong>de</strong>l cono <strong>el</strong>éctrico CPT y medición <strong>de</strong>v<strong>el</strong>ocida<strong>de</strong>s <strong>de</strong> ondas <strong>de</strong> corte Vs) y <strong>el</strong> almacenamientoy manejo <strong>de</strong> las bases <strong>de</strong> datos <strong>de</strong> casos históricos.Las actualizaciones más recientes <strong>de</strong> losprocedimientos semi-empíricos, en la evaluación <strong>de</strong>l<strong>potencial</strong> <strong>de</strong> licuación, se dirigen hacia la inclusión <strong>de</strong>expresiones numéricas que minimicen la variación entr<strong>el</strong>os patrones numéricos y permitan su estudio conherramientas básicas. Dichas expresiones pre-procesanuna variable a la vez y no mejoran la interpretación físicageneral; tratan <strong>de</strong> modificar <strong>el</strong> valor numérico <strong>de</strong> unparámetro pero no están integradas con <strong>el</strong> marco <strong>de</strong>comportamiento general y por en<strong>de</strong>, no lo alteran. Entr<strong>el</strong>os ajustes más <strong>de</strong>stacados se encuentran:a) Coeficiente <strong>de</strong> reducción <strong>de</strong> esfuerzos,b) Factor <strong>de</strong> escala <strong>de</strong> magnitud,c) Factor <strong>de</strong> corrección por sobre-esfuerzo cíclico,d) Factor <strong>de</strong> corrección por sobre estimación <strong>de</strong> laresistencia a la penetración,Los cambios en la expresión <strong>de</strong> estos coeficientesestán basados en la re-evaluación <strong>de</strong> casos históricos yla re-interpretación <strong>de</strong> pruebas y corr<strong>el</strong>acionesmejoradas (SPT, CPT y Vs → propieda<strong>de</strong>s índice/clasificación <strong>de</strong> materiales) y algunas veces, en losavances computacionales disponibles.2.1 Procedimientos Semi-empíricos: aspectosbásicosEn <strong>el</strong> procedimiento <strong>de</strong> Seed e Idriss (1971) se <strong>de</strong>fine alCSR en función <strong>de</strong> los esfuerzos vertical total yefectivo (a una profundidad Z) y la ac<strong>el</strong>eraciónmáxima horizontal en la superficie <strong>de</strong>l terreno (eng’s). El coeficiente <strong>de</strong> reducción <strong>de</strong> esfuerzo se integraa esta formulación <strong>para</strong> tomar en cuenta la flexibilidad<strong>de</strong> la columna <strong>de</strong> su<strong>el</strong>o ( , correspon<strong>de</strong> alcomportamiento <strong>de</strong> cuerpo rígido). A<strong>de</strong>más, uncoeficiente (0.65) es incluido <strong>para</strong> convertir la r<strong>el</strong>ación <strong>de</strong>esfuerzo cortante cíclico máximo al niv<strong>el</strong> <strong>de</strong> esfuerzoscorrespondiente al número <strong>de</strong> ciclos equivalentes a laexcitación sísmica particular. Los valores <strong>de</strong> CSR asícalculados pertenecen al esfuerzo cortante uniformeequivalente inducido por un sismo <strong>de</strong> magnitud M. ParaSeed e Idriss (1971) <strong>de</strong>crece con la profundidad y esin<strong>de</strong>pendiente <strong>de</strong> la configuración <strong>de</strong>l perfil <strong>de</strong> su<strong>el</strong>o.Debido a que la incertidumbre en <strong>el</strong> valor <strong>de</strong> seincrementa con la profundidad, Idriss, (1999) en unaextensión <strong>de</strong>l trabajo <strong>de</strong> Golesorkhi (1989), <strong>de</strong>sarrollóvarios cientos <strong>de</strong> análisis <strong>para</strong>métricos <strong>de</strong> respuesta <strong>de</strong>sitio y concluyó que <strong>para</strong> las condiciones <strong>de</strong> mayorinterés práctico, <strong>el</strong> parámetro podía ser expresadocomo una función <strong>de</strong> Z y M, aceptando que los esfuerzos<strong>de</strong> corte inducidos en cualquier punto <strong>de</strong>l <strong>de</strong>pósito <strong>de</strong>su<strong>el</strong>o se <strong>de</strong>ben primordialmente a la propagación vertical<strong>de</strong> ondas <strong>de</strong> corte y pue<strong>de</strong>n ser obtenidos a través <strong>de</strong>lmovimiento <strong>de</strong>l terreno (intensidad y contenido <strong>de</strong>frecuencia), <strong>el</strong> perfil <strong>de</strong> v<strong>el</strong>ocida<strong>de</strong>s <strong>de</strong> ondas <strong>de</strong> corte ylas propieda<strong>de</strong>s dinámicas <strong>de</strong> los su<strong>el</strong>os.El permite ajustar <strong>el</strong> CSR inducido durante unsismo <strong>de</strong> magnitud M a uno equivalente al generadodurante un sismo <strong>de</strong> M=7.5. El es unarepresentación aproximada <strong>de</strong> los efectos <strong>de</strong> la duración<strong>de</strong> la sacudida o <strong>el</strong> número equivalente <strong>de</strong> ciclos <strong>de</strong>esfuerzo y sus valores se <strong>de</strong>rivan <strong>de</strong> 1) corr<strong>el</strong>acionesentre <strong>el</strong> número <strong>de</strong> ciclos uniformes equivalentes y lamagnitud sísmica, y 2) pruebas <strong>de</strong> laboratorio don<strong>de</strong> ser<strong>el</strong>acione <strong>el</strong> esfuerzo cíclico requerido <strong>para</strong> causarlicuación y <strong>el</strong> número <strong>de</strong> ciclos <strong>de</strong> esfuerzo uniformes. Apesar <strong>de</strong> los intentos por <strong>de</strong>purar este parámetro, laincorporación <strong>de</strong> éste en la secuencia <strong>de</strong> cálculo siguesiendo una <strong>de</strong> las etapas más complejas e inciertas <strong>de</strong>los procedimientos semi-empíricos (Idriss y Boulanger,2004).Uno <strong>de</strong> los ajustes <strong>para</strong> tomar en cuenta <strong>el</strong> efecto <strong>de</strong>la <strong>de</strong>nsidad r<strong>el</strong>ativa y los esfuerzos <strong>de</strong> confinamientoefectivos sobre <strong>el</strong> niv<strong>el</strong> <strong>de</strong> resistencia <strong>de</strong> los su<strong>el</strong>os es <strong>el</strong>factor <strong>de</strong> corrección . Su estimación impacta <strong>de</strong>manera importante <strong>el</strong> valor <strong>de</strong>l CSR y dada lacomplejidad y dimensión <strong>de</strong> su universo, existe un grannúmero <strong>de</strong> alternativas <strong>de</strong> cálculo (Boulanger, 2003b yBoulanger e Idriss, (2004, 2003). Los métodos <strong>de</strong>análisis actuales están diseñados <strong>para</strong> trabajar con losresultados <strong>de</strong> una prueba in situ. Se <strong>el</strong>ige la querepresente ventajas <strong>para</strong> los intereses <strong>de</strong>l analista (comoposibilidad <strong>de</strong> recuperación <strong>de</strong> muestras, por ejemplo) ola mejor estandarizada (Boulanger et al. 1995, 1997). Apesar <strong>de</strong> que las conclusiones prácticas <strong>de</strong> estosestudios apuntan hacia la conveniencia <strong>de</strong> acoplar losperfiles <strong>de</strong> varias pruebas in situ <strong>para</strong> i<strong>de</strong>ntificar y


caracterizar <strong>de</strong> mejor manera los estratos<strong>potencial</strong>mente licuables, sobre todo los <strong>de</strong> espesorreducido no existe un marco <strong>de</strong> análisis que permitatomar en cuenta esta recomendación (Moss, 2003). Eluso <strong>de</strong>l número <strong>de</strong> golpes N <strong>de</strong> la prueba SPT, <strong>de</strong> laresistencia a la penetración <strong>de</strong>l cono qc <strong>de</strong> la pruebaCPT y <strong>el</strong> perfil <strong>de</strong> valores <strong>de</strong> v<strong>el</strong>ocidad <strong>de</strong> ondas <strong>de</strong> corteVs como índices <strong>para</strong> caracterizar la licuación <strong>de</strong> lossu<strong>el</strong>os, ha mostrado que su efectividad se incrementaconforme las pruebas y su interpretación se mejoran.Usando las ventajas <strong>de</strong> cada prueba (recuperación <strong>de</strong>muestras o perfiles continuos <strong>de</strong> resistencia y v<strong>el</strong>ocidad,por ejemplo) se pue<strong>de</strong> minimizar <strong>el</strong> error en lai<strong>de</strong>ntificación <strong>de</strong> los contactos entre su<strong>el</strong>os con gran<strong>de</strong>sdiferencias <strong>de</strong> resistencias y los su<strong>el</strong>os interestratificadosmuy finamente (Boulanger, 1999; Andrus y Stokoe 2000;Andrus et al. 2003).3 CÓMPUTO COGNITIVO: REDES NEURONALES,LÓGICA DIFUSA Y ÁRBOLES DE DECISIÓNEn <strong>el</strong> horizonte <strong>de</strong> los <strong>de</strong>scubrimientos r<strong>el</strong>acionados conla informática surgen no sólo i<strong>de</strong>as nuevas <strong>para</strong> mejorarlas técnicas existentes, sino también <strong>para</strong>digmasradicalmente diferentes <strong>de</strong> los actuales procedimientos.Uno <strong>de</strong> <strong>el</strong>los, actualmente aplicado en casi todas lasramas <strong>de</strong> la ciencia, es <strong>el</strong> Cómputo Cognitivo CC. Enlugar <strong>de</strong> sostenerse únicamente en las habilida<strong>de</strong>s <strong>de</strong>lmo<strong>de</strong>lador, un programa <strong>de</strong> CC apren<strong>de</strong> <strong>de</strong> suexperiencia por generalización y abstracción,adquiriendo la facultad <strong>de</strong> razonamiento y aprendizaje enun ambiente incompleto, incierto, impreciso y con verdadparcial, aspectos propios <strong>de</strong> la realidad. El CC, capaz <strong>de</strong>mo<strong>de</strong>lar y controlar una gran variedad <strong>de</strong> sistemas, seconstituye como una herramienta efectiva en problemas<strong>de</strong> toma <strong>de</strong> <strong>de</strong>cisiones, control, clasificación ycompresión, reconocimiento y predicción.En <strong>el</strong> CC se combinan diferentes técnicas <strong>de</strong> laInt<strong>el</strong>igencia Artificial como la Lógica Difusa LD, lasRe<strong>de</strong>s Neuronales RNs, la Minería <strong>de</strong> Datos MD, losAlgoritmos Genéticos AGs y <strong>el</strong> RazonamientoProbabilístico; esta última compren<strong>de</strong> AlgoritmosEvolutivos, Sistemas Caóticos, Re<strong>de</strong>s <strong>de</strong> Opinión y,parcialmente, Teoría <strong>de</strong>l Aprendizaje (Za<strong>de</strong>h, 1997,Bonissone, 1997). Conviene aclarar que <strong>el</strong> CC no es unamezcla <strong>de</strong> ingredientes sino una disciplina en la quecada componente contribuye con una metodologíaespecífica <strong>para</strong> manejar problemas en su dominio <strong>de</strong>aplicación que, <strong>de</strong> otra forma, se tornarían irresolubles.En este apartado se enuncian algunos aspectos básicos<strong>de</strong> las herramientas <strong>de</strong>l CC usados en estainvestigación, <strong>para</strong> <strong>el</strong> lector interesado en profundizarsobre cada una <strong>de</strong> <strong>el</strong>las y su aplicación en la IngenieríaGeotécnica se recomienda la revisión <strong>de</strong> García (2008).La LD, en su más popular concepción, se liga a loscálculos numéricos basados en reglas difusas que tienen<strong>el</strong> propósito <strong>de</strong> mo<strong>de</strong>lar una función numérica <strong>para</strong><strong>de</strong>scribir sistemas <strong>de</strong> ingeniería (Men<strong>de</strong>l, 1995). Sinembargo, en la literatura orientada hacia lasmatemáticas mo<strong>de</strong>rnas, se le califica como una lógicamulti-valuada cuyo propósito es mo<strong>de</strong>lar valores <strong>de</strong>verdad parcial e incertidumbres (Hajek, 1995). El<strong>de</strong>sarrollo <strong>de</strong> la LD está ligado con la investigación en IA(Int<strong>el</strong>igencia Artificial) pero sobre todo con los análisis <strong>de</strong>sistemas <strong>de</strong> ingeniería, por lo que las aplicaciones aproblemas reales tratan en las primeras etapas <strong>de</strong> su<strong>de</strong>sarrollo con la representación <strong>de</strong>l conocimiento y unavez terminadas se convierten en una interfaznumérica/simbólica.Las RNs se originaron en un intento por construirmo<strong>de</strong>los matemáticos <strong>de</strong> unida<strong>de</strong>s <strong>de</strong> procesamiento<strong>el</strong>ementales anidadas en un cerebro que permitieran <strong>el</strong>flujo <strong>de</strong> señales a través <strong>de</strong> <strong>el</strong>las. Después <strong>de</strong> un periodo<strong>de</strong> estancamiento, los mo<strong>de</strong>los formales <strong>de</strong> RNsincrementaron su popularidad con <strong>el</strong> <strong>de</strong>scubrimiento <strong>de</strong>eficientes algoritmos capaces <strong>de</strong> ajustarse acomplicadas y masivas bases <strong>de</strong> datos. Las aplicacionesneuronales más exitosas se concentran en laconstrucción <strong>de</strong> arquitecturas computarizadas capaces<strong>de</strong> aproximar funciones no lineales con un gran número<strong>de</strong> variables y clasificar objetos inmersos en espaciosmultidimensionales. Básicamente, una red neuronal esuna sofisticada función no lineal, tipo caja negra, quepue<strong>de</strong> ser entrenada con datos numéricos <strong>para</strong>satisfacer tareas objetivo.La MD es un conjunto <strong>de</strong> herramientas <strong>de</strong> análisisque sirven <strong>para</strong> <strong>de</strong>scubrir patrones, validarlos yr<strong>el</strong>acionarlos <strong>de</strong>ntro <strong>de</strong> un conjunto <strong>de</strong> datos masivo(Adriaans y Zantinge, 1996). Entre estas herramientas seencuentran las <strong>de</strong>l aprendizaje automático (algoritmosque mejoran su comportamiento automáticamente através <strong>de</strong> la experiencia) entre los más exitosos seencuentran los árboles <strong>de</strong> <strong>de</strong>cisión AD. Los AD, sean <strong>de</strong>clasificación (salida nominal) o <strong>de</strong> regresión (salidanumérica), son especialmente atractivos por tresrazones, su representación intuitiva (fácil <strong>de</strong> asimilar ycompren<strong>de</strong>r por los humanos), son mo<strong>de</strong>los no<strong>para</strong>métricos (no intervienen los sesgos <strong>de</strong>l mo<strong>de</strong>lador)y son algoritmos escalables (Breimann et al. 1984;Gehrke et al., 1998; Murthy 1995; Lim et al., 1997).En esta investigación se usan las RNs <strong>para</strong> encontrarlas reglas que r<strong>el</strong>acionan entradas y salidas (altamenteno lineales), la LD <strong>para</strong> realizar la inducción y análisis <strong>de</strong>datos <strong>de</strong>s<strong>de</strong> la experiencia expresada nominalmente ylos AR (Árboles <strong>de</strong> Regresión) <strong>para</strong> establecer unambiente <strong>de</strong> estado don<strong>de</strong> existan fronteras <strong>de</strong> múltiplesdimensiones.4 COMPORTAMIENTO DE LOS REGISTROS DECAMPO: INTERPRETACIÓN CON CCLa principal ventaja <strong>de</strong> las aproximaciones semiempíricases <strong>el</strong> uso <strong>de</strong> las consi<strong>de</strong>raciones teóricas y loshallazgos experimentales <strong>para</strong> establecer <strong>el</strong> marco <strong>de</strong>análisis y sus componentes. La teoría le da sentido a lasobservaciones <strong>de</strong> campo, validando la aproximación<strong>para</strong> interpolar o extrapolar hacia zonas con datosinsuficientes y restringirla como única solución. Lasinterpretaciones empíricas <strong>de</strong> los casos históricos, sin un


marco <strong>de</strong> análisis basado en la física, llevaría a unoscurecimiento <strong>de</strong> las condiciones <strong>para</strong> las cuales estasr<strong>el</strong>aciones son aplicables. En este trabajo se presenta alCC como la herramienta capaz <strong>de</strong> cumplir con estanecesidad <strong>de</strong> abstracción y mo<strong>de</strong>lación en ambientes <strong>de</strong>trabajo inciertos y contaminados.Con base en una extensa revisión <strong>de</strong> la literatura en<strong>el</strong> tema, <strong>de</strong> las conclusiones y juicios ligados a lasobservaciones, y <strong>de</strong> las características <strong>de</strong> los estudios<strong>de</strong> campo y laboratorio que conforman las bases <strong>de</strong>datos, se <strong>el</strong>igieron los factores con más peso y mayorcertidumbre <strong>para</strong> incluirlos en <strong>el</strong> mo<strong>de</strong>lo <strong>de</strong> CC. Estosmo<strong>de</strong>los se <strong>de</strong>sarrollaron procurando mantenerloscompletos (amplia gama <strong>de</strong> comportamientos ycondiciones iniciales), simples y fundamentados eninterpretaciones físicas <strong>de</strong> las respuestas <strong>de</strong>l material.De acuerdo con <strong>el</strong> proceso empírico <strong>de</strong> creación <strong>de</strong>mo<strong>de</strong>los se plantea <strong>el</strong> siguiente esquema <strong>de</strong> solución alproblema <strong>de</strong> obtener <strong>el</strong> <strong>potencial</strong> <strong>de</strong> licuación <strong>de</strong> unsu<strong>el</strong>o:• SISTEMA: <strong>de</strong>pósito <strong>de</strong> su<strong>el</strong>o.• PROCESO: Colapso estructural <strong>de</strong>l esqu<strong>el</strong>eto <strong>de</strong>su<strong>el</strong>o <strong>de</strong>bido a esfuerzos <strong>de</strong> corte, con unacoexistente pérdida <strong>de</strong> energía por mecanismos <strong>de</strong>fricción entre partículas.• ENTRADAS: parámetros <strong>de</strong>l su<strong>el</strong>o/ parámetros <strong>de</strong> lacarga.• MODELO: Técnica Híbrida (un Sistema LD, una RNy un AR).• SALIDA: “licuación: si”/ “licuación: no” y CRR/CSR.4.1 SistemaDiversas condiciones <strong>de</strong> carga y estados iniciales en losmateriales <strong>de</strong>ben ser coinci<strong>de</strong>ntes <strong>para</strong> que la licuaciónocurra. El sistema no tiene restricciones sobre dichascondiciones. La integración <strong>de</strong> los materiales quecomponen un <strong>de</strong>pósito no requiere <strong>de</strong> caracterizacionesque dirijan la salida hacia alguna zona <strong>de</strong>l espacio <strong>de</strong>respuestas. Por ejemplo, etiquetar al su<strong>el</strong>o como“cohesivo”/“no cohesivo”, <strong>de</strong>ntro <strong>de</strong> un análisistradicional dirige <strong>el</strong> razonamiento. Si <strong>el</strong> material es “nocohesivo”entonces <strong>el</strong> término licuación se refiere alcomportamiento asociado a las presiones <strong>de</strong> poro enexceso, <strong>de</strong>formaciones significativas y la posible pérdida<strong>de</strong> resistencia, en campo se i<strong>de</strong>ntifica por <strong>de</strong>formaciones<strong>de</strong> terreno y volcanes <strong>de</strong> arena, mientras que en <strong>el</strong>laboratorio pue<strong>de</strong> ser <strong>de</strong>finida en términos <strong>de</strong> algúncriterio <strong>de</strong> falla (alcanzar <strong>el</strong> 3% <strong>de</strong> <strong>de</strong>formación porcorte, por ejemplo).Si por otro lado, la etiqueta <strong>de</strong>l material es“cohesivo”, <strong>el</strong> término licuación no es usado <strong>para</strong><strong>de</strong>scribir <strong>el</strong> fenómeno, ya que se sabe que la cargacíclica en las arcillas pue<strong>de</strong> producir exceso en lapresión <strong>de</strong> poro, pero usualmente menor al 100%.(Zergoun y Vaid 1994, Boulanger et al. 1998, Idriss1985, Idriss 1985). ¿Pero qué suce<strong>de</strong>ría en <strong>el</strong> caso <strong>de</strong>materiales que exhiben comportamiento intermedio? La<strong>para</strong>metrización rígida no permite evaluar condicionesnaturales e introduce sesgos en <strong>el</strong> patrón <strong>de</strong> respuestaestimado. En este trabajo <strong>de</strong> investigación <strong>el</strong> “sistema”es <strong>de</strong>scrito a través <strong>de</strong> parámetros índice, mecánicos ygeométricos (Figura 1), con fronteras <strong>de</strong> expresiónaproximadas. La observación ligada con <strong>el</strong> arreglonatural permite una completa exploración y explotación<strong>de</strong> la experiencia en campo y laboratorio.ENTRADASGeotécnicasvariable numéricavariable nominalGeométricasFigura 1. <strong>ARELI</strong>: entradas SUELO.Perfiles<strong>de</strong> campo(continuos)Pesovolumétrico*ClasificaciónSUCS*Descripciónmateriales*Espesor <strong>de</strong>los estratos*Z - NAF*Z - partesuperior <strong>de</strong>lestrato enestudioV s , q c , N = f (Z)5 etiquetas (membresías)3 etiquetas (membresías)Arcilla, Limo, Arena y Gravay sus combinaciones(porcentajes SUCS con7 etiquetas)Se integra a <strong>ARELI</strong> en formanumérica cerrada.4.2 Proceso


El colapso estructural <strong>de</strong>l esqu<strong>el</strong>eto <strong>de</strong> su<strong>el</strong>o, como <strong>el</strong>término general <strong>para</strong> <strong>de</strong>scribir los fenómenos asociadoscon la licuación (falla <strong>de</strong> flujo y movilidad cíclica) no semo<strong>de</strong>la en <strong>el</strong> esquema matemático expuesto en estetrabajo <strong>de</strong> investigación, sino que se vincula la expresiónnominal <strong>de</strong> una <strong>de</strong> sus manifestaciones evi<strong>de</strong>ntes (“estesu<strong>el</strong>o licuó”/“este su<strong>el</strong>o no licuó”) con algunosparámetros que caracterizan al evento generador <strong>de</strong>lProceso. La solución CC es un planteamientofenomenológico, es <strong>de</strong>cir, <strong>el</strong> mo<strong>de</strong>lo es capaz <strong>de</strong><strong>de</strong>scribir las manifestaciones <strong>de</strong>l proceso yeventualmente r<strong>el</strong>acionarlas con causas, pero no <strong>de</strong>explicar la estructura que las genera.1.1 EntradasLos parámetros <strong>de</strong> su<strong>el</strong>o que permiten construir <strong>el</strong>ambiente <strong>de</strong> análisis se <strong>de</strong>tallan en la Figura 1. Unsistema <strong>de</strong> LD funge como traductor <strong>de</strong> los parámetrosnominales a numéricos <strong>de</strong> acuerdo con la experienciaprevia <strong>de</strong> casos históricos (entrenamiento neuronal <strong>de</strong>un sistema difuso) y como transfigurador <strong>de</strong> losparámetros continuos (perfiles <strong>de</strong> resistencias yv<strong>el</strong>ocida<strong>de</strong>s) <strong>para</strong> hacer evaluaciones discretas <strong>de</strong>ntro<strong>de</strong>l ambiente “<strong>de</strong>pósito <strong>de</strong> su<strong>el</strong>o”. Una <strong>de</strong> las ventajas <strong>de</strong>usar LD en la expresión <strong>de</strong>l sistema es que involucraduplas <strong>de</strong> variables nominales y numéricas <strong>para</strong>construir la imagen <strong>de</strong>l medio. Los parámetrosseñalados como obligatorios cubren <strong>el</strong> mínimo <strong>de</strong>información geotécnica y geométrica con <strong>el</strong> que <strong>ARELI</strong>pue<strong>de</strong> comenzar un análisis, sin embargo <strong>el</strong> sistema <strong>de</strong>LD podría completar al sistema <strong>de</strong> acuerdo con laexperiencia y conocimiento <strong>de</strong>l analista.La introducción <strong>de</strong> la carga sísmica se pue<strong>de</strong> realizara través <strong>de</strong>: a) la magnitud y la ac<strong>el</strong>eración máxima o b)un ac<strong>el</strong>erograma (Figura 2). El uso combinado <strong>de</strong> <strong>el</strong>lastambién es permitido. La entrada <strong>de</strong> carga <strong>de</strong>pen<strong>de</strong>rá<strong>de</strong>l rigor necesario en <strong>el</strong> análisis y la disposición <strong>de</strong>información histórica y registros <strong>de</strong> ac<strong>el</strong>eracionespertinentes. Es motivo <strong>de</strong> investigación más profunda la<strong>de</strong>finición <strong>de</strong>l dominio en <strong>el</strong> que se expresen mejor lascaracterísticas <strong>de</strong> la carga que inci<strong>de</strong>n sobre lageneración <strong>de</strong> la licuación (series <strong>de</strong> tiempo, espectrosen <strong>el</strong> dominio <strong>de</strong> la frecuencia o espectros temporofrecuenciales).ENTRADASAc<strong>el</strong>erograma*Dominio <strong>de</strong>l tiempo¨*Dominio <strong>de</strong> la Frecuencia3 02 010001 00100- 1 0Sísmicas- 2 0- 3 0- 4 00 5 0 1 0 0 1 5 0 2 0 01010 1 10variable numéricaMagnitud MRango dinámico: 4


don<strong>de</strong> <strong>el</strong> analista pue<strong>de</strong> estudiar <strong>el</strong> efecto que cadaaspecto <strong>de</strong>l sistema y <strong>de</strong> la carga tienen sobre lar<strong>el</strong>ación entre la capacidad <strong>de</strong> resistencia <strong>de</strong> un su<strong>el</strong>o y<strong>el</strong> esfuerzo cíclico que le afecta (Figura 3). A<strong>de</strong>más, ensu hoja final, <strong>el</strong> árbol concluye con un valor numérico,apto <strong>para</strong> la i<strong>de</strong>a <strong>de</strong> la práctica <strong>de</strong> valorizar <strong>el</strong> riesgo conbase en este número.2 <strong>ARELI</strong>: EJEMPLO DE APLICACIÓNintegración <strong>de</strong> patrón <strong>de</strong> entrada (Figura 4) <strong>para</strong> <strong>el</strong> sitioBorah Peak (EU), durante <strong>el</strong> sismo <strong>de</strong> Idaho en 1983(Mw = 6.9) (Andrus y Stokoe, 1996). Por limitaciones <strong>de</strong>espacio se presenta <strong>el</strong> esquema <strong>de</strong> entrada más simple,pero <strong>ARELI</strong> contiene conceptos y funciones <strong>de</strong>membresía sobre otras condiciones <strong>de</strong> sitio (pruebasíndice, <strong>de</strong>scripciones geológicas, geomorfológicas ygeotécnicas) y su r<strong>el</strong>ación lingüística con un aspecto <strong>de</strong>la licuación <strong>de</strong> los materiales (<strong>de</strong>scripción <strong>de</strong>tallada <strong>de</strong>las fuentes bibliográficas en García et al., 2010).Con base en los criterios y metodología establecidos enlos capítulos anteriores, se presenta un ejemplo <strong>de</strong>RN<strong>ARELI</strong>Entradasnuméricas+Juicios <strong>de</strong>expertos+Bloques <strong>de</strong> ReglasEsquematización <strong>de</strong>la respuestaÁRBOL DE REGRESIÓNRecomendación <strong>de</strong>salidaretroalimentaciónProcesoCognitivoFigura 3. Módulo Neuronal y Árbol <strong>de</strong> Regresión: proceso <strong>de</strong> aprendizaje con retroalimentación.


Prof. Z* H Esf.Vert Esf.Vert Ef Tipo Vs amax Licuación?Sitio NAF (m) (m) (m) (kPa) (kPa) Su<strong>el</strong>o (m/s) (g)RN ObservadaPence Ranch 1.7 1.8 1.9 57.2 46.2 Arena 107 0.36 0 01.5 1.5 2.8 52.7 40.5 gravosa 94 0.36 0 01.4 1.4 1.8 44.5 36 a 102 0.36 0 01.8 1.8 2.8 62.1 49.4 Grava 109 0.36 0 01.5 1.5 1.9 60.5 45.6 Arenosa 122 0.36 0 02 2 1.7 57.5 46.3 134 0.36 0 01.5 1.5 1.7 38.8 32.9 128 0.36 0 01.5 1.5 1.9 38.4 32.4 107 0.36 0 01.5 1.5 1.7 39.4 33.8 131 0.36 0 01.7 1.7 1.5 43.3 38.3 122 0.36 0 01.5 1.5 2.3 48.5 38.1 154 0.36 0 0Goddard 1.2 1.2 2 47.3 36 Grava arenosa 122 0.3 1 0Ranch 1.2 1.2 2 41.1 32.7 Grava arenosa 105 0.3 0 0An<strong>de</strong>rsen Bar 0.8 0.8 2.4 40.6 28.7 Grava arenosa 106 0.29 0 00.8 0.8 2.4 39 27.8 105 0.29 0 0Larter Ranch 0.8 2.2 1.3 59.9 39 Grava areno-limosa 176 0.5 0 00.8 2.2 1.3 55.4 38.4 153 0.5 0 00.8 2.2 1.3 59.9 40.5 183 0.5 0 0Whiskey 0.8 1.8 2.2 59.1 38.2 Grava limo-arenosa 181 0.5 0 0Springs 0.8 1.8 2.2 45.6 31.7 Grava limo-arenosa 210 0.5 0 0North Grav<strong>el</strong> 1 1.8 1.2 51 36 Grava arenosa 206 0.46 1 1Bar 3 3 1.3 75.2 53.5 Grava arenosa 274 0.46 1 1Mackay Dam, 2.3 2.3 2.7 66.6 57.4 Grava areno-limosa 271 0.23 0 1Z* Profundidad a la frontera superior<strong>de</strong>l estratoHEspesor <strong>de</strong>l estratoEsf.Vert Esfuerzo vertical (Total)Esf.Vert Ef Esfuerzo vertical (Efectivo)0 – Licuación SÍ1 – Licuación NOFigura 4. Ejemplo <strong>de</strong> Tabla <strong>de</strong> Entrada y Salidas calculada y observada.Con <strong>ARELI</strong> se pue<strong>de</strong> partir <strong>de</strong>l juicio <strong>de</strong> los expertos,hacerlo empatar con los conceptos almacenados en lamemoria difusa <strong>de</strong> <strong>ARELI</strong> (Youd y Perkins, 1987;correctamente contaban con información contradictoriaentre los niv<strong>el</strong>es <strong>de</strong> resistencia a la penetración y <strong>el</strong> tipo<strong>de</strong> material.Kramer, 1996, por mencionar algunos <strong>de</strong> los trabajos <strong>ARELI</strong> finalmente bosqueja <strong>el</strong> patrón <strong>de</strong>lintegrados) y <strong>de</strong>ducir <strong>el</strong> valor numérico que másconvenga al análisis. En este ejemplo la <strong>de</strong>mandasísmica se integra con a max y M.Cuando este patrón <strong>de</strong> entrada llega al móduloneuronal la salida binaria indica “licuación: sí”. Laclasificación acertada <strong>de</strong> <strong>ARELI</strong> está por encima <strong>de</strong>lcomportamiento <strong>de</strong> este caso <strong>de</strong> análisis como semuestra en la Figura 5. La com<strong>para</strong>ción <strong>de</strong> losresultados <strong>de</strong> simulación y los <strong>de</strong> <strong>ARELI</strong> rev<strong>el</strong>a la altacapacidad <strong>de</strong> predicción <strong>de</strong>l mo<strong>de</strong>lo y sus remarcablesventajas en <strong>el</strong> estudio y aprendizaje <strong>de</strong> la licuación y losprocesos r<strong>el</strong>acionados.88% <strong>de</strong> los casos usados en la prueba <strong>de</strong>l mo<strong>de</strong>loneuronal, siendo importante aclarar que un gran número<strong>de</strong> patrones que no pudieron ser evaluados


Figura 5. Ejemplo <strong>de</strong>l patrón <strong>de</strong> comportamiento y su interpretación con <strong>el</strong> esquema <strong>de</strong> análisis <strong>ARELI</strong>.3 CONCLUSIONESExiste un gran número <strong>de</strong> métodos que <strong>el</strong> ingenieropue<strong>de</strong> s<strong>el</strong>eccionar cuando analiza, clasifica o resu<strong>el</strong>veproblemas <strong>de</strong> regresión. Las técnicas cognitivas cuando“trabajan” y “producen” predicciones numéricas oclasificaciones acertadas basadas en un reducidonúmero <strong>de</strong> condiciones lógicas, representan una granventaja sobre otras alternativas escasas y “oscuras”, <strong>de</strong>difícil acceso y con poca utilidad fuera <strong>de</strong> las bases <strong>de</strong>datos con las que fueron construidas.Los resultados finales <strong>de</strong> <strong>ARELI</strong> se pue<strong>de</strong>n resumircomo una serie <strong>de</strong> condiciones lógicas si-entonces(nodos <strong>de</strong>l árbol), siendo la RN la que ratifica yretroalimenta dichas condiciones y la LD la que traduce<strong>el</strong> conocimiento y la experiencia <strong>para</strong> completar <strong>el</strong>análisis. En <strong>ARELI</strong> no existen hipótesis implícitas sobr<strong>el</strong>as r<strong>el</strong>aciones entre las variables <strong>de</strong> entrada y las<strong>de</strong>pendientes, por lo que esta herramienta es i<strong>de</strong>al <strong>para</strong><strong>el</strong> estudio <strong>de</strong> fenómenos r<strong>el</strong>acionados con los sismos ylos su<strong>el</strong>os, dón<strong>de</strong> existe poco conocimiento a priori sobr<strong>el</strong>as ligas <strong>para</strong>métricas o <strong>el</strong> conjunto <strong>de</strong> teorías parece enalgunos sentidos incongruente (importancia <strong>de</strong> variablesy cómo están r<strong>el</strong>acionadas).La interpretación <strong>de</strong> los resultados <strong>de</strong> <strong>ARELI</strong> es muysimple. Esta simplicidad es útil no sólo cuando serequiere una rápida clasificación (o predicción) <strong>de</strong>nuevas observaciones, sino porque tener un mo<strong>de</strong>lo mássimple <strong>para</strong> explicar las observaciones y los resultados<strong>de</strong> clasificación (o predicción) permite <strong>el</strong>aborarconceptos, conclusiones y conocimientos “naturales”.Se <strong>de</strong>muestra con los resultados presentados en estetrabajo que, en <strong>el</strong> estudio <strong>de</strong>l <strong>potencial</strong> <strong>de</strong> licuación, losmétodos <strong>de</strong>l CC pue<strong>de</strong>n ayudar a rev<strong>el</strong>ar r<strong>el</strong>acionessimples, explotando las variables <strong>de</strong> fácil y económicaobtención en campo y laboratorio.4 REFERENCIASAdriaans P y Zantinge D, (1996). Data Mining , NewYork: Addison Wesley.Aghda S, Teshnehlab A, Suzuki T., Akiyoshi T yKitazono, (1998), Liquefaction Potential AssessmentUsing Multilayer Artificial Neural Network, J. Sci. I.R.Iran, Vol.9, No.3.Andrus, R.D., y Stokoe, K.H.,(2000) II. "Liquefactionresistance of soils from shear-wave v<strong>el</strong>ocity," Journalof Geotechnical and Geoenvironmental Engrg.,ASCE, 126(11), 1015-1025.Andrus, R.D., Stokoe, K.H., II, Chung, R.M., Juang, C.H.(2003), "Gui<strong>de</strong>lines for evaluating liquefaction resistanceusing shear wave v<strong>el</strong>ocity measurements andsimplified procedures." NIST GCR 03-854, NationalInstitute of Standards and Technology, Gaithersburg,MD.Andrus, R.D.y Stokoe, K.H., (1996), “LiquefactionResistance Based on Shear Wave V<strong>el</strong>ocities”, Proc.NCEER Workshop on Eval. Liquefaction Resistanceof Soils, Eds. Youd and Idriss, NCEER-97-0022.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!