INSTITUTO POLITECNICO NACIONAL

Escuela Superior de Ingeniería y Arquitectura

Unidad Ticomán

Notas de Clase
Programa Académico: Ingeniería Geofísica

Alumno: López Aguirre Eric Uriel

No. De Boleta: 2021390131

Profesora: Gabriela de los Santos Cano

Materia: Pre-Procesado Sísmico

Periodo: 2024/2
 Factores de atenuación de la señal
sísmica
Afectan a la amplitud:
➢ Perdidas por transmición
➢ Perdidas por divergencia esferica
Afectan a la amplitud, frecuencia y fase:
➢ Perdidas por absorción inélastica
➢ Perdidas por dispersión

Transmición
Para una onda plana que incide normalmente en la interface de dos
capas con diferente impedancia acústica

Continuidad del esfuerzo normal

De aquí obtenemos el coeficiente de transmición y el coeficiente de
reflexión:

Coeficiente de transmición Coeficiente de reflexión

Divergencia esférica
Sí la velocidad no varía con la
profundidad, el decaimiento de
amplitud está relacionado con la
distancia de propagación.

Si la velocidad varía en la dirección vertical, el decaímiento de

amplitud tiene la forma:
 Absorción inelástica

Factor Q.
Se define como:

Suponiendo que un medio con relación lineal esfuerzo-deformación;

la amplitud de la onda es proporcional a E.

De la cual obtenemos las fluctuaciones de amplitud de acuerdo a la

𝜋
atenuación. Dada una amplitud incial 𝐴0 y 𝐴 decrementa en
𝑄
tiempos sucesivos 2𝜋/𝜔, 4𝜋/𝜔, … ,2𝑛𝜋/𝜔.

Para 𝑡 = 2𝑛𝜋/𝜔
Por lo tanto:
Pero 𝜔 = 2𝜋𝑓; (frecuencia angular) podemos obtener la perdida de
amplitud respecto a la distancia de propagación.

Asi:

Es llamado coeficiente de absorción, por lo que:

Dispersión

Se debe a que todas las componentes

de frecuencia viajan con velocidades
diferentes.

𝑣𝜔 − Distorsión de fase
Siendo:

El pulso en el dominio del tiempo estará dado por:

Sustituyendo valores:
 Ganancias (Escalamiento)
Para poder entender de mejor manera este tema debemos saber
algunos aspectos:
➢ Entre más denso sea el material la velocidad de propagación
será mayor.
➢ Puede ser mismo material, pero diferentes condiciones.
➢ No debemos alterar las amplitudes.
➢ Para las escalas siempre es recomendable multiplicarlas por un
factor constante.
La ganancia es una escala variante en el tiempo. Existen diferentes
tipos de ganancia:
➢ Balanceo
➢ PGC
➢ AGC
➢ Corrección por dispersión geométrica
Las ganancias se aplican a los datos sísmicos, para despliegue, se
usan para levantar la señal débil.
La ganancia debe usarse con cuidado, ya que puede destruir el
carácter de la señal.

Balance
Es una escala de amplitudes invariante en el tiempo. Está basado en
el criterio de amplitud
RMS. Específicamente
cada traza de un grupo
es escalada, así que
todas tienen el mismo
valor de amplitud RMS
deseado.
Valor RMS: Es la medida estadística de la magnitud de una cantidad
variable. Es la raíz cuadrada de la media aritmética de los cuadrados
de los valores, con estos cuadrados evita el problema del signo.
Amplitud RMS: Refleja el volumen percibido mejor que la amplitud
pico.

PGC
Es un tipo de ganancia simple, una función de ganancia puede estar
definida por la interpolación entre algunos valores escalares
especificados para una muestra de tiempo particular. Los valores
escalares mayores, naturalmente deben asignarse a los tiempos
mayores.
 AGC
Existen dos tipos de AGC:
➢ AGC Amplitud RMS (Normalización)
➢ AGC Amplitud instantánea (Normal)}

AGC Amplitud RMS (Normalización)

La función de ganancia con amplitud RMS con una ventana de
tiempo especifica en una traza de entrada. Esta se calcula:
1. La traza de entrada es subdividida en ventanas de tiempo
arregladas.
2. Se ajusta la amplitud de cada muestra en una ventana.
3. Se calcula en medio de este valor y se toma su raíz cuadrada
(amplitud RMS de esta ventana).
4. El radio de la amplitud RMS deseado es asignado al valor RMS
actual como función de ganancia al centro de la ventana. Por lo
tanto, la función de escalamiento para el centro de esa ventana
está dado por:

Normalmente se inicia con una

ventana que abarque la parte
superficial de la traza. La longitud
de la ventana puede ser constante
o puede incrementarse sistema-
ticamente, así como avanza la
traza.
 AGC Amplitud instantánea (Normal)
Es una de las más comunes. Esta función se calcula:
1. Se calcula el valor medio absoluto de la amplitud de la traza
dentro de una ventana de tiempo específica.
2. El radio del nivel RMS deseado es asignado como el valor de la
función de ganancia. Este valor se asigna a la n-esima muestra
de la traza.
3. Se mueve la ventana de tiempo una muestra abajo, se calcula
el valor de la función de ganancia para (n+1)-esima muestra de
tiempo, y así sucesivamente.
4. No se requiere interpolación para definir la función de
ganancia. Así la función de escalamiento g(t) de la ventana, está
dada por:

Ventanas de tiempo muy pequeñas significan pérdidas del carácter

de la señal, por otro lado, con ventanas muy grandes, se disminuye
la efectividad de la ganancia.
En la práctica, las ventanas se especifican entre los 200 y 500 ms.
 Corrección por dispersión geométrica
Un sismograma representa un campo de ondas generado por una
fuente. Un tiro es considerado una fuente puntual que tiene un
campo de ondas esférico.
La tierra tiene dos efectos en la propagación de las ondas:
1. En un medio homogéneo, la densidad de energía (E=energía x
volumen unitario) decae proporcionalmente a 1/𝑟 2 donde r es el
radio del frente de ondas.
La amplitud de la onda es 𝐴𝛼√𝐸 y decae 1/𝑟.
En la practica la velocidad incrementa con la profundidad, la cual
causa mayor divergencia del frente de onda y un decaimiento mas
rápido en la amplitud, con la distancia.
2. El contenido de frecuencias de la señal inicial cambia de un
modo variante con el tiempo.
En particular las altas frecuencias son absorbidas más
rápidamente que las bajas frecuencias, esto se debe a la
atenuación intrínseca de las rocas.
El factor que describe el decaimiento de la amplitud de la onda
como una función del radio del frente de ondas esférico es válido
para un medio homogéneo sin atenuación. Para la Tierra esférica,
el decaimiento de amplitud puede describirse por:

Velocidad: vector de cantidad que indica el intervalo de tiempo en el

cual existe un cambio de posición de una partícula.
En sísmica, la velocidad se refiere a la propagación de las ondas
sísmicas en intervalos de tiempo muy pequeños sin implicar su
dirección.
Velocidad RMS: se refiere a una trayectoria de rayo específica,
suponiendo capas horizontales.

Por lo tanto, la función de ganancia para la compensación por

dispersión geométrica está definida por:
 Ruido Sísmico
Modelo convolucional ideal:

Modelo convolucional real:

Toda aquella parte de la señal que no proporciona información útil.

El ruido generalmente se clasifica en dos categorías (coherente e
incoherente).
Las características del ruido sísmico son:
➢ Frecuencia
➢ Amplitud
➢ Fase

Ruido Coherente
Mantiene continuidad a lo largo del sismograma.
➢ Tienes de energía sistemáticos con velocidad aparente uniforme
(refracciones, difracciones, ruido de pozo, etc).
 Ruido Coherente

Se encuentra alineado entre

las diferentes trazas es
introducido por fuentes
externas
➢ Viento, actividad
humana, irregularidades en
el subsuelo, instrumental,
etc.

Filtros
Dispositivo que modifica el espectro de una señal.
Tiene como objetivo obtener una función de transferencia H(z)
realizable, aproximándola a una respuesta de frecuencias deseada;
este debe ser casual y estable.
➢ Tiempo
➢ Frecuencia
Filtro en el dominio del tiempo:
Filtro en el dominio de las frecuencias:

Dominio Transformada Dominio

t Fourier 1D F
t-x Fourier 2D F-K x
t-x τ−p τ − px

Transformada de Fourier

Transformada inversa de Fourier

Muestreando siempre para preservar todas las frecuencias hasta la
frecuencia de Nyquist.

La transformada de Fourier de una serie de tiempo de longitud 2N

se puede hacer a través de dos transformadas de Fourier de longitud
N.
Tenemos como requisito:
La serie de tiempo tiene que ser potencia de 2 (2𝑛 ).

Espectro de amplitud del filtro:

TF
Espectro de frecuencias de la traza sísmica. Diseño del filtro
Estos filtros no tienen una buena respuesta de esta forma, es
necesario relajar las características de bajada del filtro en las bandas
de paso y rechazo de banda (es decir, suavizarlos).
Truncar la respuesta al impulso de un filtro para hacerlo realizable,
crea oscilaciones en el dominio de las frecuencias, estas oscilaciones
son conocidas como el fenómeno de Gibbs.
Cuando una función tiene una discontinuidad de salto en un punto,
su serie de Fourier tiene un comportamiento especial en dicho punto.
En procesamiento de imágenes se les llama artefactos de anillo.
En audio se perciben en instrumentos de persecución (sonido
producido por un golpe de tambor, de batería, platillos, etc).
Entre mayor sea el número de coeficientes (M) que usemos para
nuestro filtro, mayor será el número de ondulaciones, sin embargo, el
ancho de las ondulaciones disminuye. En cambio, si el número de
coeficientes del filtro incrementa la amplitud de las oscilaciones
disminuirá.
El lóbulo principal se hace más angosto, conforme se incrementan los
coeficientes de filtro.
El fenómeno de Gibbs es el resultado
de la operación de ventaneo
(multiplicar la respuesta al impulso
del filtro ideal con una función de
ventana regular).
 COMPLEMENTO
Atenuacion:
La atenuación reúne los efectos de la absorción anelástica y del
scattering y expresa el decaimiento de la amplitud de las ondas
cuando se propagan.
Un sismograma local puede considerarse en un principio, compuesto
de ondas directas que han sufrido absorción y ondas secundarias que
además de absorción han sufrido desviaciones de sus trayectorias con
la consecuente redistribución de energía.
Si consideramos un medio homogéneo, la
amplitud de una onda que se propaga
disminuye, en primera instancia, por la
geometría de su frente de onda. Esto es
debido a la divergencia esférica, en el
caso de ondas internas, o divergencia
cilíndrica, en el caso de ondas
superficiales.
La atenuación inelástica de las ondas sísmicas es el proceso por el
cual la energía de las ondas que viajan en las rocas se convierte en
calor. Depende de varias características como la composición de la
matriz, la saturación de fluidos, la porosidad, la densidad de
fracturas y otros. La atenuación de ondas sísmicas generalmente se
representa como el factor de calidad Q ó el coeficiente de atenuación.
Ganancia:
La ganancia es una función de tiempo escalada, en la cual, el
escalamiento de la función está basado en un criterio deseado. La
señal al propagarse por el interior de la tierra sufre una atenuación
de las altas frecuencias.
Además, el frente de onda aumenta proporcionalmente a la distancia
epicentral, de manera que la energía se tiene que distribuir en este
frente de ondas. Esto hace que la amplitud disminuya con el tiempo.
Para compensar esta pérdida de energía los datos se deben de
corregir por este efecto. La corrección es proporcional al tiempo
(distancia recorrida). Esta corrección se denomina divergencia
esférica.
Simplificando la fuente como un punto generador de ondas esféricas,
en un medio homogéneo la densidad de energía decae de la forma
1/𝑟 2 , donde r es el radio del frente de ondas. La amplitud de la onda
es proporcional la energía, la cual decae 1/r.
La corrección que se aplica para compensar el decaimiento de la
amplitud con el tiempo es la “geometric spreading correction”.
Esta corrección realza la señal de los tiempos más tardíos. Esta
corrección no restaura las amplitudes de las altas frecuencias ya que
estas están sometidas a una fuerte atenuación.
Para remover el efecto de atenuación se procede a modificar el
espectro de amplitud de la señal, con etapas de procesado conocidas
como deconvolución y, compensación espectral de tiempo variable.
Existen diferentes tipos de ganancia que se pueden aplicar. Una
función de ganancia g (t) consiste en una serie de escalares que se
multiplican con las amplitudes y se pueden derivar de los datos y son
Para ello existen tres tipos de ganancia: control de ganancia
programada (programed gain control, PGC), el control de ganancia
automática (AGC, automatic gain control) con amplitud cuadrática
media (root mean square RMS). Esta última es el tipo de ganancia
más comúnmente usado.
Ejemplo de ganancia de tiempo
escalado basada en una función g (t).
En la parte superior se ilustra antes
de la aplicación de la ganancia y en la
parte inferior la respuesta luego de
multiplicar la función g(t) por la
amplitud de la traza de entrada.
Las AGC, PGC, etc., utilizan una medida estadística de la amplitud
dentro de una ventana. Esta medida estadística puede ser, el RMS,
la media, la mediana, etc.
El cálculo de la medida estadística se realiza dentro de la ventana, la
ventana recorre toda la traza a intervalos de una muestra, generando
una serie numérica de la misma longitud que la traza.
Esta serie se multiplica por la serie de amplitudes de la traza. Si el
tamaño de las ventanas es muy pequeño puede causar una pérdida
significante del carácter de la señal. En la práctica utilizan ventanas
entre 256 y 1024 ms.
En algunos casos, se tiene, que las características de la amplitud
varían de traza en traza. Esto puede ser debido a diferencias de
amplificación de los sensores, diferencias de acoplamiento entre el
sensor (geófono) y el suelo, etc. Para compensar, corregir estas
diferencias se aplica un balanceado entre trazas. Esto significa que
todas las trazas deben tener un nivel de amplitud RMS similar.
El espectro sísmico, especialmente para frecuencias altas, está sujeto
a la absorción de la propagación debido a la atenuación. Las altas
frecuencias se atenúan con la profundidad.
En consecuencia, el contenido en frecuencia de la señal varía con el
eje de tiempos, de manera que los filtros deben variar también con el
tiempo. En la práctica, los filtros son aplicados entre ventanas
adyacentes para obtener una transición suavizada de las regiones del
filtro pasobanda.
Dispersión Geométrica:
Una buena analogía de estos
mecanismos la obtenemos de la Óptica. A
medida que nos alejamos de una lámpara
en la calle la intensidad o brillo de la luz
disminuye debido a algunas razones
como: la luz se propaga desde la lámpara
en forma de frentes de ondas esféricas.
Por conservación de la energía, la energía por unidad de área del
frente de ondas a medida que se aleja de la fuente decrece como 𝑟 2 ,
donde r es el radio de la esfera o la distancia a la lámpara.
La energía del frente de ondas esféricas, en el campo cercano, que se
originan en una fuente puntual se distribuye sobre una superficie
esférica de radio creciente.
Ruido sísmico:
En el contexto de la prospección sísmica; Dobrin y Savit (1988)
definen ruido como señales sísmicas no deseables del subsuelo no
asociadas con reflexiones.
El ruido sísmico es cualquier cosa, menos la señal deseada. El ruido
incluye las perturbaciones de los datos sísmicos causadas por una
energía sísmica indeseada, tales como la onda superficial generada
por los tiros, las ondas de superficie, las múltiples, los efectos del
intemperismo y de la actividad humana, o las ocurrencias aleatorias
de la Tierra.
El ruido puede ser minimizado mediante el empleo de arreglos de
fuentes y receptores, la generación de un nivel mínimo de ruido
durante la adquisición, y el filtrado y apilamiento de los datos
durante el procesamiento.
El término señal se usa para denotar cualquier evento en el registro
sísmico del cual se desea obtener información. Todo lo demás es
ruido, incluyendo eventos coherentes que interfieren con la
observación y medición de las señales. La relación entre señal y ruido,
abreviada como S/R, es la relación de la energía de la señal en una
porción específica del registro y la energía total del ruido en la misma
porción.
Cuando la relación señal/ruido es pequeña, se obtienen registros muy
pobres; el nivel a que ésta se considera pequeña es, en gran medida,
un juicio subjetivo. Sin embargo, cuando S/R es menor que la unidad,
usualmente la calidad del registro es marginal y se deteriora
rápidamente a medida que decrece la relación.
El ruido sísmico puede ser: coherente, o bien, incoherente. El ruido
coherente se puede seguir cuando menos en unas cuantas trazas; el
ruido incoherente es disímil en todas las trazas y no se puede predecir
cuál traza será parecida, aun conociendo las trazas vecinas.
A veces la diferencia entre ruido coherente e incoherente es
simplemente cuestión de escala y si se tuvieran hidrófonos
espaciados más estrechamente, el ruido incoherente se pudiera
observar cómo ruido coherente. No obstante, el ruido incoherente se
define con respecto a los registros que se estén usando, sin importar
qué tan cerrado debería ser el espaciamiento para revelarlo.
El ruido incoherente se denomina a menudo ruido aleatorio
(espacialmente aleatorio), lo cual implica no sólo impredictibilidad
sino también ciertas propiedades estadísticas; con mayor frecuencia
el ruido no es verdaderamente aleatorio. Debe observarse que la
aleatoriedad espacial y temporal pueden ser independientes; la traza
sísmica usual tiende a ser aleatoria en tiempo, ya que no se sabe
cuándo ocurrirá la reflexión con base en lo que ha mostrado
previamente la traza, con la excepción de los múltiples.
Algunas veces el ruido coherente se subdivide en: energía que viaja
sobre todo horizontalmente, y energía que llega al tendido
verticalmente. Otra importante distinción es entre: ruido que es
repetible, y ruido que no lo es; es decir, si se observa el mismo ruido
al mismo tiempo sobre la misma traza cuando se repite un tiro.
Las tres propiedades coherencia, dirección de viaje y repetibilidad
forman las bases de la mayoría de los métodos para mejorar la
calidad del registro.
El ruido coherente incluye ondas superficiales, reflexiones o
refracciones reflejadas de estructuras cercanas a la superficie, como
planos de falla o cauces de corrientes subterráneas, refracciones
conducidas por intercalaciones delgadas de alta velocidad, ruido
causado por tráfico vehicular o tractores agrícolas, multiples, etc.
Todos los anteriores, excepto los múltiples, viajan esencialmente
horizontales y todos excepto el ruido vehicular es repetibles en tiros
sucesivos.
El ruido incoherente que es espacialmente aleatorio y también
repetible se debe a dispersión proveniente de irregularidades
cercanas a la superficie e inhomogeneidades como cantos rodados,
fallas de pequeña escala, etc.
Debido a que hay muchos tipos de ruido, es necesario emplear varios
métodos de atenuación de ruido. Todos están basados en las
diferencias entre las propiedades del ruido y de la señal. Puesto que
la naturaleza de la “señal” es algo subjetiva y las propiedades tanto
de la señal como del ruido no son completamente conocidas, la
atenuación de ruido no puede ser completamente objetiva.
La relación señal a ruido (S/R) generalmente depende de la
frecuencia. Si el ruido presenta energía fuera de la principal gama de
frecuencias de la señal, el filtrado de frecuencia se usa
ventajosamente. Las componentes de frecuencia muy baja (como las
ondas superficiales de baja energía ricas en frecuencias bajas) se
pueden filtrar durante el registro inicial ya que las bajas frecuencias
se separan suficientemente de las frecuencias de reflexión.
Sin embargo, si el espectro del ruido traslapa con frecuencia al
espectro de la señal, entonces el filtrado de frecuencia tiene un valor
limitado para mejorar la calidad del registro.
 Referencias
-Divergencia esférica - SEG Wiki. (s. f.).

https://wiki.seg.org/wiki/Dictionary:Spherical_divergence/es

-GrilledSeismic. (s. f.). Sismología de exploración: Procesamiento de datos sísmicos, parte

1. https://geocafemx.blogspot.com/2016/02/sismologia-de-exploracion-

procesamiento.html

-Ruido sísmico, https://glossary.slb.com/es/terms/n/noise

-Sheriff, R.E. & Geldart, L.P. (1991), Exploración sismológica: Historia, Teoría y

obtención de datos, (1° edición, Vols. 1-2).