13.07.2015 Views

Libro de resúmenes [revisión final, 172 páginas] - UGM

Libro de resúmenes [revisión final, 172 páginas] - UGM

Libro de resúmenes [revisión final, 172 páginas] - UGM

SHOW MORE
SHOW LESS

Create successful ePaper yourself

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

SISMOLOGÍA Geos, Vol. 27, No. 1, Octubre, 2007SIS-17ROTATING PROJECTION ALGORITHM ANDAPPLICATION FOR SEISMIC TOMOGRAPHYGrebennikov Alexan<strong>de</strong>r y Valencia Pérez TomásFacultad <strong>de</strong> Ciencias Físico Matemáticas,Benemérita Universidad Autónoma <strong>de</strong> Pueblaagrebe@fcfm.buap.mxOne of the well known approaches for recuperation of theun<strong>de</strong>rground structure by seismic methods can be formulatedmathematically as a seismic tomography problem [1], which isanalogue of X-rays tomography [2], [3]. For <strong>de</strong>tailed recognitionof investigating structures it is necessary to use a lot of measureddata and to reconstruct images by the inversion of the Radontransform. In some simplified cases, for example at preliminaryrecognition, it is possible to suppose that the investigating objecthas a “discrete” structure, i. e. the structure that consists ofsome separate elements with different characteristic insi<strong>de</strong> ofthe homogeneous (or quasi homogeneous) substance in theconsi<strong>de</strong>red domain.Un<strong>de</strong>r these suppositions we propose to use the RotatingProjection algorithm [4], [5]. It simplifies the image reconstructionfor the consi<strong>de</strong>ring case, because it does not require complicatednumerical methods as (regularized) calculation of the inverseRadon transform or solving mal conditioned systems of linearalgebraic equations [1]. The proposed algorithm was <strong>de</strong>velopedand realized as program package in MATLAB system, whichquality was <strong>de</strong>monstrated on simulated numerical examples. Thisalgorithm may be applied for another type of tomography [5],hence it has some perspectives in electromagnetic geophysicsmethods [6].Author acknowledge to SEP y CONACYT Mexico for thesupport of the part of this investigation in the frame of the ProjectCB-2006-1-57479.REFERENCES1. A. N. Tikhonov, A.V. Goncharsky at all, Seismic Tomographyin Engineer Geophysics Problems. Dokl. Acad. Nauk USSR,304(4), 840-844 (1989).2. R.A. Williams and M.S. Beck, ProcessTomography: Principles, Techniques and Applications.Butterworth-Heinemann, Oxford, 1995.3. M.S. Beck and B.H. Brown, Process tomography: aEuropean innovation and its application. Measurement Scienceand Technology, 7, 215–224 (1996).4. Alexandre Grebennikov, Linear regularization algorithmsfor computer tomography. Inverse Problems in Science andEngineering J., Vol. 14, No. 1, January, 53–64 (2006).5. A. Grebennikov, J. G. Vazquez Luna, M. A. Cruz Gama, FastLinear Algorithms for Image Processing in Electric Tomography. IIITaller Internacional sobre Física Aplicada. La Habana, Noviembre30 – Diciembre 3, 2004.6. A. N. Tikhonov. Mathematical Geophysics (collected papersedited by V.I. Dmitriev). Un. Inst. of Earth Phys., Russian Acad.of Sciences, Moscow, 1999.SIS-18DETERMINACION DE LA ESTRUCTURA RADIALDE VELOCIDADES EN LA BASE DEL MANTOUSANDO ONDAS S DIFRACTADAS DE LA RED MOMATerán Mendieta Luis Fernando 1 , Valenzuela Wong Raúl 1 ,Wysession Michael E. 2 , Fischer Karen M. 3 y Clarke Timothy J. 41 Instituto <strong>de</strong> Geofísica, UNAM2 Department of Earth and PlanetarySciences, Washington University, USA3 Department of Geological Sciences, Brown University, USA4 Corrales, New Mexico, USAluisteran@yahoo.comLa técnica <strong>de</strong> la constante <strong>de</strong> disipación ha sido usada conanterioridad para <strong>de</strong>terminar la estructura radial <strong>de</strong> velocida<strong>de</strong>sen la base <strong>de</strong>l manto a partir <strong>de</strong> ondas que rozan el núcleo o quese difractan fuera <strong>de</strong> él. En particular, Valenzuela y Wysession[1998] propusieron mo<strong>de</strong>los <strong>de</strong> velocidad para las ondas S <strong>de</strong>bajo<strong>de</strong> Siberia Oriental y para una región <strong>de</strong>bajo <strong>de</strong>l Océano Pacíficocentro-oriental. En ambos casos, los mo<strong>de</strong>los que mejor ajustanlos datos tienen un aumento súbito <strong>de</strong> la velocidad en el interfazsuperior <strong>de</strong> la región D” seguido <strong>de</strong> un <strong>de</strong>cremento gradual hastallegar al núcleo. Ahora presentamos los resultados <strong>de</strong> aplicar lamisma técnica a un volumen en D” <strong>de</strong>bajo <strong>de</strong>l Océano Pacíficocentro-oriental, localizado al noroeste <strong>de</strong> la región estudiada porValenzuela y Wysession [1998]. Empleamos datos <strong>de</strong>l sismo<strong>de</strong> Tonga <strong>de</strong>l 7 <strong>de</strong> abril <strong>de</strong> 1995 registrados en estaciones enel noreste <strong>de</strong> los Estados Unidos <strong>de</strong> Norteamérica. Escogimosun perfil que incluye 14 estaciones, cubre una distancia <strong>de</strong> 14°en D” y cuya ventana acimutal es <strong>de</strong> sólo ~1°. A diferencia <strong>de</strong>trabajos anteriores [Valenzuela y Wysession, 1998], se llegó a laconclusión <strong>de</strong> que un mo<strong>de</strong>lo <strong>de</strong> uso común, el Mo<strong>de</strong>lo Preliminar<strong>de</strong> Referencia para la Tierra (PREM, por sus siglas en Inglés)<strong>de</strong> Dziewonski y An<strong>de</strong>rson [1981], ajusta bien la constante <strong>de</strong>disipación observada. Dicho mo<strong>de</strong>lo tiene un espesor <strong>de</strong> 150 kmen D” y el gradiente <strong>de</strong> velocidad para la onda S es cercano acero. Este resultado es similar al obtenido por Schlittenhardt et al.[1985] ya que no <strong>de</strong>tectaron la existencia <strong>de</strong> una discontinuida<strong>de</strong>n una región que coinci<strong>de</strong> parcialmente con la parte suroeste <strong>de</strong>lpresente estudio. En general, nuestras conclusiones concuerdancon trabajos anteriores que <strong>de</strong>stacan la complejidad <strong>de</strong> D” <strong>de</strong>bajo<strong>de</strong>l Océano Pacífico, pues ahí coexisten regiones cercanas condiscontinuida<strong>de</strong>s, sin discontinuida<strong>de</strong>s, zonas <strong>de</strong> baja velocidad yestratos <strong>de</strong> lentitud y <strong>de</strong>lga<strong>de</strong>z extremos. Por otra parte, tambiénanalizamos el sismo <strong>de</strong>l 23 <strong>de</strong> agosto <strong>de</strong> 1995 en las IslasMarianas. Este evento nos permitió estudiar D” <strong>de</strong>bajo <strong>de</strong>l sur <strong>de</strong>Alaska y en una región al suroeste <strong>de</strong> la Península <strong>de</strong> Alaska.En este caso, las estaciones más lejanas <strong>de</strong>l perfil muestranlas mayores amplitu<strong>de</strong>s para la onda S difractada. De hecho,las mayores amplitu<strong>de</strong>s se producen a altas frecuencias, <strong>de</strong>0.025 a 0.180 Hz (periodos <strong>de</strong> 40 a 5.6 s), solamente paralas estaciones más alejadas. Valenzuela y Wysession [1998]no pudieron obtener mediciones confiables <strong>de</strong> la constante<strong>de</strong> disipación a frecuencias tan altas para otros eventos. Seemplearon 75 mo<strong>de</strong>los diferentes, incluyendo mo<strong>de</strong>los con ysin discontinuida<strong>de</strong>s así como mo<strong>de</strong>los con una zona <strong>de</strong>lgada<strong>de</strong> ultrabaja (o bien ultraalta) velocidad, pero ninguno <strong>de</strong> ellosajustó las observaciones. Las formas <strong>de</strong> onda <strong>de</strong> este eventoparecen indicar que la técnica <strong>de</strong> la constante <strong>de</strong> disipacióntiene la resolución suficiente para distinguir variaciones en laestructura <strong>de</strong> D” a escalas relativamente pequeñas. Sin embargo,es necesario realizar más trabajo al respecto.120

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

Saved successfully!

Ooh no, something went wrong!