0% encontró este documento útil (0 votos)
123 vistas319 páginas

Métodos de Geofísica: Gravimetría y Sismología

Este documento presenta una introducción al método gravimétrico de geofísica. Explica la ley de gravitación de Newton, cómo varía la gravedad debido a la forma de la Tierra, el concepto de potencial gravitatorio y superficies equipotenciales. También cubre cómo la gravedad varía con la latitud y altitud, y define conceptos como el esferoide normal y el geoide.
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd
0% encontró este documento útil (0 votos)
123 vistas319 páginas

Métodos de Geofísica: Gravimetría y Sismología

Este documento presenta una introducción al método gravimétrico de geofísica. Explica la ley de gravitación de Newton, cómo varía la gravedad debido a la forma de la Tierra, el concepto de potencial gravitatorio y superficies equipotenciales. También cubre cómo la gravedad varía con la latitud y altitud, y define conceptos como el esferoide normal y el geoide.
Derechos de autor
© © All Rights Reserved
Nos tomamos en serio los derechos de los contenidos. Si sospechas que se trata de tu contenido, reclámalo aquí.
Formatos disponibles
Descarga como PDF, TXT o lee en línea desde Scribd

PRIMER SEMESTRE

GEOFISICA

ASESOR: RAMON AVILES JIMENEZ


2. Geofísica

2.1. Métodos Potenciales

2.1.1. Método Gravimétrico (Gravimetría)

[Link]. Ley de Newton

[Link]. Efecto de la forma en la gravedad de la Tierra

[Link]. Potencial Gravífico

[Link]. Superficie Equipontencial

[Link]. Corrección de la deriva instrumental

2.1.2. Interpretación

[Link]. Obtención del mapa de anomalías de Bouguer

[Link]. Interpretación Geológica

[Link]. Cálculo de los mapas de anomalías regional y residual

[Link]. El método de interpretación de la segunda derivada

[Link]. Calculo de mapas de segunda derivada

[Link].1. Interpretación cuantitativa

[Link]. Efectos gravíficos de cuerpos de formas geométricas sencillas


[Link]. Determinación del regional en un perfil

2.2. Métodos Magnéticos

2.2.1. Generalidades sobre el magnetismo

2.2.2. Ley de las acciones magnéticas

2.2.3. Unidad de intensidad magnética

2.2.4. Campo magnético

2.2.5. Medición de la susceptibilidad magnética de las rocas


[Link]. Tabla de susceptibilidades de las rocas

2.2.6. Magnetometría

2.3. Sismología

2.3.1. Adquisición de Datos Sísmicos

2.3.2. Exploración Sísmica

[Link]. Equipo Geofísico

[Link].1. Principios Básicos de Operación

[Link].2. Sistema Perfilador (Batimetrico)


[Link].3. Sistema de Barrido del Fondo Marino

[Link].4. Hidrofonos

[Link].5. Selección del Equipo

[Link]. Presentación de la Información

[Link].1. Plano de posicionamiento

[Link].2. Plano Batimétrico

[Link]. Analizado

[Link].1. Plano de Isopacas

[Link].2. Plano Estructural Somero

[Link]. Tipos de Discordancias

[Link].1. Gas

[Link].2. Coral

[Link]. Métodos Directos

[Link].1. Trabajos Marítimos

[Link].2. Exploración y Muestreo

[Link].3. Pruebas de Barco


[Link]. Generalidades del método sísmico

[Link].1. Ondas Sísmicas

[Link]. Procesado de Datos Sísmicos

[Link].1. Mapas Sísmico 3D

[Link].2. Bloque del Encabezado

[Link].3. Identificación y Estructura de un Formato

[Link].4. Formato SEG-D

[Link].5. Multiplexado

[Link].6. Formatos SEG-Y

[Link].7. Construcción de Geométrica

[Link].8. Intervalo de Muestreo

[Link].9. Correcciones de Estáticas

[Link].9.1. Estáticas de elevación y de refracción

[Link].9.2. Estática de Elevación

[Link].9.3. Estática de Refracción


[Link].9.4. Reflexiones Fantasmas

[Link].9.5. Múltiples

[Link].9.6. Migración

[Link].9.7. Migración en Tiempo

[Link].9.8. Difracciones

[Link].9.9. Ruido Coherente y Aleatorio

[Link].9.10. Decovolucion

[Link].9.11. Análisis de Velocidades

[Link].9.12. Migración 2D/3D antes de apilamiento en Tiempo


(Kirhhoff)

[Link].9.13. Migración 2D/3D antes de apilamiento en


Profundidad (Kirhhoff)

[Link].9.14. Migración en tiempo post apilamiento

[Link]. Interpretación de Datos Sísmicos

[Link].1. Atributos Sísmico

[Link].1.1. Atributos Pre-Strack

[Link].1.2. Atributos Pos-Apilados


[Link].1.3. Atributos Instantáneos

[Link].1.4. Atributos de onda

[Link].1.5. Atributos Físicos

[Link].1.6. Atributos Geométricos

[Link].1.7. Atributos de reflectividad y transmisión


2. Geofísica

2.1. Métodos Potenciales

2.1.1. Método gravimétrico

El método gravimétrico de prospección se basa en la medida, en superficie, de las


pequeñas variaciones (o anomalías) de la componente vertical del campo gravífico
terrestre. Estas variaciones son debidas a una distribución irregular en profundidad
de masas de diferentes densidades, por lo que conociendo a ellas se puede llegar
a una interpretación más o menos probable de la situación de las masas en el
subsuelo, dependiendo del conocimiento del geofísico y de la distribución de
densidades en profundidad.

Tres son los aparatos empleados en el método gravimétrico: la balanza de torsión


que mide las derivadas de la gravedad, el péndulo que mide la gravedad absoluta
o la relativa, y el gravímetro que mide la gravedad relativa.

[Link]. Ley de Newton

La atracción de las masas obedece a la ley de gravitación de Newton que


establece que la atracción de dos masas m1 y m2 es directamente proporcional a
su producto e inversamente proporcional al cuadrado de su distancia.

2
F=K((m1m2)/r )

Constante de la gravitación universal.


En las mediciones de la atracción de la gravedad interviene una constante k que
es la constante de la gravitación universal, donde en la Ley de Newton, cuando
9
m1=m2=1 gramo y r=1 cm F=K. El valor de K es K=66,70.10 con dimensiones M-
1 L3 T-2.

El problema de la medición fue efectuada por Cavendish utilizando la balanza que


lleva su nombre.
[Link]. Efecto de la forma en la gravedad de la Tierra

Debido al achatamiento ( el radio de la Tierra no es constante), g varía con la


latitud.

De acuerdo con el sistema Geodésico de Referencia (G.R.S.67), el valor de g al


nivel del mar varía con la latitud , de la siguiente forma:

2 4
g  978.031,8(10,005278895 sen   0,000023462sen )mgal

En esta expresión  está medido en [rad]. Recordar que: [rad] = 180º y 1[gal] =
2
1[cm/seg ]

gEcuador= gmin = 978,049[gales] gPolos = gmax = 983,2213[gales]

En Geofísica, si se desea usar la gravedad como herramienta en determinaciones


estructurales, interesan más los valores relativos de g que los valores absolutos.
Lo que realmente interesa es poder definir contrastes de inhomogeneidades y
densidades en la parte de la Tierra que se está investigando; Esto significa
descubrir aspectos que de alguna manera se apartan de la uniformidad. En
resumen lo que interesa son las anomalías gravitatorias.

[Link]. Potencial Gravifico

El potencial gravífico terrestre consta de dos partes fundamentales; la primera y


principal la producida por la Tierra de acuerdo a la Ley de Newton. La segunda es
la debida a la rotación de la Tierra, la atracción debida a otros astros es pequeña
por eso no se considera. Como ambas partes están formadas por fuerzas
actuando sobre un punto fijo y además la magnitud de estas fuerzas depende de
la distancia a un punto las dos partes tendrán un potencial.
[Link]. Superficie Equipotencial

Estas son aquella son las que el potencial es constante; por lo tanto la fuerza de
gravedad será perpendicular a las superficies equipotenciales, ya que en caso
contrario dicha fuerza produciría u trabajo al transportar la unidad de masa a lo
largo de un superficie equipotencial o de nivel y esto no puede ser según se
deduce en la definición misma de superficie equipotencial. La superficie
equipotencial estaría representada por la superficie exterior de una masa líquida
en equilibrio sometida a la acción de la gravedad. Variación de la gravedad sobre
la superficie terrestre

a) Con la latitud.

En 1672 Richer refería, sin explicación, que la gravedad variaba de un


punto a otro de la Tierra. Poco después Newton y Huygens lo explicaban
diciendo que este fenómeno existía debido a quela Tierra no era de forma
esférica.

En efecto, como la Tierra está achatada por lo polos, la distancia a su


centro es máxima en el ecuador (valor de la gravedad mínima), y mínima en
los polos (valor de gravedad máximo). A esta variación se suma la del
efecto de rotación de la Tierra, la fuerza centrífuga que es máxima en el
ecuador y nula en los polos y siempre opuesta a la fuerza de la gravedad.
Como consecuencia de ambos efectos la asceleración de la gravedad
varía, aproximadamente de 978 cm/seg2 en el ecuador a 983 cm/seg2 en
los polos.

Esferoide normal.

Si suponemos la Tierra pastosa formada por capas homogéneas


concéntricas y sometidas a la fuerzas de gravedad que tienden a hacerla
esférica y a la fuerza centrífuga que tiende a aplastarla, la forma que
tomaría este forma ideal se llama el esferoide normal; este es ligeramente
diferente a un elipsoide de revolución.
La formula internacional de la gravedad utilizada para este esferoide data
de 1930 y es:

g0= 978,049 (1+ 0.0052884 sen2 de phi – 0.0000059 sen 2 de 2phi)


cm/seg2
Donde g0 es la gravedad a la latitud phi, y al nivel del mar; el factor 978,049
es el valor de la gravedad en el ecuador (phi=0) con esta fórmula se calcula
el valor normal o teórico de la gravedad g0 en cualquier latitud.

El geoide

Se llama geoide a la superficie de equilibrio de los mares de la Tierra si


estos se pudieran extender por los continentes a través de canales
imaginarios hechos en ellos bajo el nivel del mar. Es decir, sería la
superficie equipotencial correspondiente al nivel de los mares. Difiera del
esferoide en que los continentes y los fondos marinos están irregularmente
repartidos sobre la Tierra y, así mismo, en las discontinuidades de la
densidad de la Tierra.

Como la plomada (normal al geoide), tiende a desviarse hacia las regiones


de exceso de masa (como en los continentes) y a alejarse de las de déficit
de masa (como en los océanos), la consecuencia es que el geoide está por
encima del esferoide en los continentes y por debajo en los océanos. Fig.
del libro 11-4

b) Con la altitud

Debido a que hay mediciones que se hacen no a nivel del mar sino, sobre
la superficie terrestre, con una altura h, sobre dicho nivel, se tienen tres
correcciones debidas a la altitud: corrección de aire libre, corrección de
Bouguer, y corrección topográfica: estas correcciones permitirán reducir la
gravedad observada al nivel del mar.
1) Corrección de aire libre: Es la debida simplemente a la altura h de la

estación, es decir, a estar más alejada del centro de la Tierra.

g h  0,3086 h mgal, (esta corrección también se llama,

corrección por altura).


La estación representada por la altura h sobre el nivel del mar, está a
una distancia h más alejada del centro de la Tierra que otra estación que
estuviere al nivel del mar. Puesto que la masa de la Tierra se puede
considerar como concentrada en su centro, la ley de los cuadrados
inversos nos indica que la atracción gravitacional a una altura h será :

R2
.
g  R  h g 2
0

 
Esto significa que:

R2
g  g 0  g  1  R  h 2 g 0 
2g hR  g h2

2g h

2

 0 o 0

 

Rh
R

ya que h  R Si se reemplazan los valores de g0 y de R = 2.678 Km.


llegamos al valor indicado más arriba. Esta corrección se llama de aire
libre porque no se tiene en cuenta la atracción de ningún material
situado por encima del nivel del mar.

2) Corrección de Bouguer: Se calcula hallando el efecto gravífico del


terreno comprendido entre el nivel del mar, o sea suponiendo que la
estación y el nivel del mar, o sea, suponiendo que la estación se ha
hecho en un plano topográfico horizontal.

h mgal.
2
g B  2Gh  4,19110

3
Si :[gr/cm ] y h[m]
Esta corrección tiene en cuenta la atracción del material rocoso situado
entre el nivel del mar y la estación situada a una altura h. Se basa en la
hipótesis de que la superficie de la Tierra es horizontal en todas partes
(paralela al geoide) a una altura h por encima del nivel del mar. Las
montañas que sobresalen de esta superficie horizontal imaginaria, y los
valles que quedan por debajo falsean esta hipótesis, pero su efecto
gravitatorio se compensa por la corrección topográfica subsiguiente.
Esta corrección se resta porque en realidad estamos eliminando el
material situado entre el nivel del mar y el nivel de la estación.
3) Corrección topográfica: Para, reducir finalmente el valor observado de la
gravedad al valor que tendría al nivel del mar se debe tener en cuenta el
efecto gravífico de las masas por encima y por debajo del nivel h de la
estación, ya que en la corrección de Bouguer se suponía al terreno
horizontal. El material por encima de la estación atraerá a la masa
unidad puesta en A en el sentido opuesto a la fuerza de la gravedad,
mientras el material b (por debajo) actuará en sentido contrario. Puesto
que el efecto del material del valle b habría que restarlo a la corrección
de Bouguer y esta es negativa, la corrección total topográfica irá
siempre sumada a la gravedad observada.

c) Con el tiempo

1) Efecto de las mareas: Las atracciones gravíficas del sol y la luna causan
el efecto de las mareas que depende de la posición astronómica de
ambos y de la latitud, siendo variable con el tiempo; estas fuerzas
producen una pequeña deformación de la superficie marina. En efecto
influye también sobre los gravímetros y da lugar al efecto lunisolar o de
las mareas.

Su amplitud máxima puede llegar en total a 0.3 miligales, pero su


variación máxima es solamente de unos 0.05 mgal/hora.

El efecto de las mareas es poco importante en la prospección, ya que se


elimina parcialmente, al efectuar la corrección de deriva.
2) Deriva instrumental: El gravímetro, como todo instrumento de precisión,
da lugar a lo que se llama deriva instrumental, es decir que si
efectuamos medidas sobre una misma estación a diversos intervalos de
tiempo obtenemos valores ligeramente distintos. Esta deriva está
producida principalmente por la fatiga del sistema de muelles, o por las
variaciones de temperatura.

Observaciones de campo

Composición de un Equipo Gravimétrico

1 Jefe de equipo cuya misión es el planteamiento del trabajo a realizar, así como
su supervisión e interpretación.

1 Calculador de oficina cuya misión es el cálculo de la deriva instrumental y de


todas las correcciones para llegar a los valores del mapa de anomalías de
Bouguer.

1 Operador cuya misión es la lectura del gravímetro en los puntos de estación.

1 Ayudante de Operador, que generalmente conduce el vehículo del Operador y


ayuda a la apuesta en estación y nivelación del gravímetro. Topógrafos

Auxiliares y vehículos.

Toma de datos

Las lecturas se realizan con el gravímetro, una vez situado en estación y nivelado,
tomando siempre nota de la hora en que se realiza cada estación.
Con objeto de realizar la prospección con mayor posible se establece
primeramente una Red de Bases que cubra lo más regularmente posible, de una
manera amplia, la zona de trabajo. Una base es simplemente un punto en donde
la gravedad se determina con mucha precisión.

La red de bases se establece sobre una serie de puntos de tal manera dispuesto
que el intervalo de tiempo entre las lecturas de dos bases consecutivas sea lo
menor posible y, asimismo, que haya facilidad para observar repetidas veces las
bases a lo largo de la red.

El enlace entre las bases sucesivas B1, B2, … se observarían en el siguiente


orden: B1 -B2- B1-B3- B2-B3… Si se necesitase realizar una red de mucha
precisión, se repetirían más veces los itinerarios indicados.
Una vez establecida la red de bases, el resto de los puntos de estación podrán
referirse fácilmente a los puntos de la red, con lo que todo el trabajo quedará
perfectamente enlazado.

[Link]. Corrección de la deriva instrumental

La deriva es la diferencia entre dos lecturas realizadas en una misma estación en


un cierto intervalo de tiempo.

Estas lecturas teóricamente deberían ser iguales pero no lo son debido a la


infidelidad del gravímetro como ocurre con todos los instrumentos de precisión.
Sus diferencias en los gravímetros modernos suele ser de algunas centésimas de
miligal por día. Por lo tanto las lecturas del instrumentos deberán ser corregidas
con los valores que obtengamos de la deriva, de tal manera que los valores
corregidos se aproximen lo más posible a los que hubiéramos obtenido si todas
las estaciones las hubiéramos realizado simultáneamente y con el mismo
instrumento.

2.1.2. Interpretación

El mapa de Bouguer muestra la suma de todos los efectos debidos a todas las
masas presentes tanto en profundidad como en superficie.

Es decir, el mapa de Bouguer mostraría la suma de los efectos debido a :

 Las rocas sedimentarias próximas.



 Al basamento ígneo.

 Al magma en el interior de la Tierra.

[Link]. Obtención del mapa de anomalías de Bouguer

Al realizar observaciones en el campo medimos las variaciones relativas de la


componente vertical del campo gravífico terrestre. Los valores obtenidos, deben
ser corregidos para tener en cuenta los siguientes factores:

1) Variaciones de la gravedad a lo largo del día producidas por los cambios de


posición del sol y de la luna y por deriva instrumental.
2) Diferencias de altitud de las estaciones gravimétricas entre sí y con relación al
plano d referencia.

3) Forma de la Tierra: esferoide normal que da lugar a la corrección de latitud.

[Link]. Interpretación geológica

El mapa de Bouguer muestra la suma de todos los efectos debidos a todas las
masas presentes tanto en profundidad como en superficie. Es decir, que el mapa
de Bouguer mostraría las sumas de los efectos de las rocas sedimentarias
próximas del basamento ígneo, e incluso del magma en el interior de la tierra.

Sin embargo, los efectos de las masas anómalas muy profundas están muy
atenuados debido a su gran distancia produciendo las variaciones amplias de las
curvas isoanómalas, mientras que las estructuras sedimentarias próximas
producen variaciones relativas agudas y frecuentemente fuertes del campo
gravífico terrestre. Por lo tanto en el mapa de anomalías de Bouguer las
variaciones amplias serían debidas a contrastes de densidad profundos o
alejados, generalmente en el basamento, y las variaciones rápidas serán
producidas por contrastes de densidad más próximos.

De aquí podemos deducir el principal objetivo de la interpretación gravimétrica.

a) La separación de los efectos de las masas profundas y someras.

b) La definición de los efectos gravimétricos de las anomalías sedimentarias,


determinando la probable profundidad y dimensiones de las estructuras
interpretadas.
[Link]. Cálculo de los mapas de anomalías regional y residual

El mapa de anomalías de Bouguer podemos considerarlo como la suma de otros


dos: 1) l mapa regional, debido al efecto de las masas profundas y lejanas y 2) el
mapa residual debido a los efectos de las masas superficiales.

En prospección interesa hacer resaltar el efecto debido a las masas superficiales,


eliminando o reduciendo al mínimo el debido a las masas profundas. Por ello se
define corrientemente la anomalía regional como el efecto de todo aquello que no
interesa al prospector. La anomalía residual sería:

Anomalía residual= Anomalía de Bouguer – Anomalía regional.

A. Métodos gráficos

1) Método de suavización de curvas isoanómalas. (o métodos de suavización de


contornos).

2) Método de los perfiles. Difiere solamente del anterior en el procedimiento de


suavización utilizado.

El procedimiento consistiría simplemente en trazar áreas perfiles sobre el mapa de


Bouguer en dirección del gradiente máximo (es decir, de la variación máxima del

δg) y separados por distancias del orden de anomalías esperadas. Se considera


que en cada perfil el gradiente regional es constante (una recta) y de esta manera
se van calculando los valores de gravedad residual que permitirá generar el mapa
de anomalías residuales.

B. Métodos analíticos

1) Método de las medidas aritméticas. Supongamos que en una malla de


estaciones de observación hemos trazado las curvas isoanómalas

[Link]. El método de interpretación de la segunda derivada

El método de interpretación de la segunda derivada de la gravedad consiste en


determinar, partiendo de los valores de la gravedad observados en el mapa de
Bouguer, la segunda derivada vertical de la gravedad. Es decir, si tomamos unos
ejes coordenados con eje OZ en la vertical hacia abajo, el problema será
2
 g

determinar z 2 .

El método tiene la ventaja de que señala mejor los cambios bruscos de la


gravedad y, asimismo, que permite separar las anomalías gravimétricas
compuestas (o sea, las debidas a la superposición de los efectos de varias masa
próximas) en sus componentes separadas.

El método de la segunda derivada es un método de mucho más poder


resolutivo que los métodos anteriores ya que la doble diferenciación de g respecto
a z tiende a hacer resaltar las anomalías más pequeñas y superficiales a
expensas de las más amplias y profundas de tipo regional. Con frecuencia el
2
 g

mapa z 2 ofrece una imagen más clara y exacta de los tipos de anomalías que
interesan en prospección petrolífera.

Asimismo el método tiene diversos inconvenientes que indicaremos a


continuación:

1) No se puede aplicar si las estaciones de la red están muy espaciadas.

2) Las observaciones gravimétricas se deben efectuar con alta precisión


debiendo, por tanto, realizar con mucho cuidado las correcciones
topográficas.

3) Las segundas derivadas no están relacionadas directamente con las


estructuras que la causan, por lo que es muy difícil de entender la
relación entre las masas del subsuelo y las segundas derivadas. Por ello
la forma de una anomalía no se puede deducir del mapa de segundas
derivadas.

4) Como las derivadas segundas decrecen con las potencias de la


profundidad de las masas, las anomalías profundas quedan anuladas en
el mapa de segundas derivadas.

[Link]. Calculo de mapas de segunda derivada


Método de Rosenbach.

Sea en el plano z = 0 el punto 0, centro de una circunferencia de radio r, en el que


2g
querernos hallar el valor de

2
z

El método consiste simplemente en suponer los puntos


9 de la malla lo suficientemente próximos a 0 de modo
que podamos expresar su gravedad por un desarrollo en serie. De esta manera y
considerando los cuatro puntos de la figura adjunta podremos escribir, de acuerdo
con el desarrollo de Mc Laurin.

[Link].1. Interpretación cuantitativa

Hasta ahora los métodos de interpretación que hemos indicado trataban


solamente de separar las distintas anomalías y definirlas en términos de posibles
anticlinales, sinclinales, domos, masas de mayor o menor densidad que las rocas
circundantes… etc.

Al tratar que persigue la interpretación cuantitativa: en prospección


petrolífera determinar el tamaño y la profundidad de las posibles trampas; en
prospección minera determinar, principalmente, la masa que produce una
determinada anomalía.

La interpretación cuantitativa se basa en el primer caso en la comparación


de la anomalía obtenida con la producida por cuerpos de forma previamente
establecida. Si ambas coinciden para una determinada forma, profundidad y
distribución de densidades, ésta sería una posible solución; sin embargo, pueden
existir otras distribuciones de masas, a diferentes profundidades, que produzcan la
misma anomalía y únicamente la ulterior información geológica haría posible la
obtención de la solución más probable.
La interpretación cuantitativa en este caso es indirecta ya que trata de
asimilar las anomalías encontradas a las producidas por cuerpos de la índole de
los que puedan ser hallados en la zona objeto de la prospección. Por ello no se
podrá intentar hacer una interpretación cuantitativa sin un conocimiento geológico
previo de las condiciones aproximadas de las profundidades de las masas,
densidades en juego, forma y tamaño de las masas a esperar… etc.

Los procedimientos de interpretación serán:


1) Comparar los efectos observados con los cuerpos de formas
geométricas sencillas (esfera, cilindro horizontal… etc.)

2) Comparar los efectos observados con los calculados para cuerpos de


forma cualquiera

En prospección minera generalmente, se busca la masa que produce una


determinada anomalía. En este caso la magnitud de la masa se puede determinar
como una solución única, siempre que se conozca el contraste de densidad del
mineral que vaya buscándose.

[Link]. Efectos gravificos de cuerpos de formas geométricas sencillas

Uno de los métodos de interpretación en gravimetría es comparar las


anomalías obtenidas del mapa gravimétrico residual con el efecto producido por
cuerpos de formas geométricas sencillas enterrados en el subsuelo. Esta
comparación puede dar una orientación sobre la forma real y la profundidad de las
masas enterradas.

Por ello es conveniente estudiar el efecto producido por masas de forma


geométrica sencilla tales como la esfera, el cilindro horizontal, la capa horizontal
infinita de espesor constante y la capa horizontal semi-infinita de espesor
constante o falla.

No entraremos en el rigor matemático de las demostraciones que pueden


deducirse simplemente considerando el efecto gravífico de un elemento
infinitésimo y extendiendo su efecto al cuerpo considerando mediante la
integración correspondiente. Observemos que siempre consideramos la
componente vertical de la gravedad y que la densidad σ que adoptemos será el

contraste de densidad (o diferencia de


densidades) entre el cuerpo considerado y la
masa homogénea que le rodea.
La esfera.- una esfera de masa m y densidad σ (o contraste de densidad),
situada a la profundidad z bajo el punto o. Puede demostrarse que el efecto de
esta esfera sobre la unidad de masa colocada en un punto cualquiera de la
superficie, p (x, o) es el mismo que el que resultaría suponiendo la masa de la
esfera concentrada en su centro C.

Por tanto:
m
gk
siendo r = C P r2
La componente vertical será:

z z
g  g cos  g  km

z r r3

4 3

m  R 
Ahora bien, como 3 , resultará
3
4 R  z
g k 
z 3 (z2  x2 )3/ 2

gz
que dará en gales, cuando z, x y R estén expresados en centímetros.

gz
será máxima cuando x = 0 y llamándola g0 tendremos:
4 R 
3

g g k
z 2
0 máx . 3z

gz 1
Si representamos la curva  en función de obtenemos la
2 x
g0  3/ 2
 x 
1  y
   

  z  

 
curva. Calculando la distancia horizontal LM entre sus puntos
de inflexión

x
observamos que resulta ser igual a la profundidad de la masa  1; x  z .
y
gz
Por tanto bastara representar la curva para obtener LM, que será del orden de
g0

la profundidad del centro de la esfera.

Un domo salino puede representarse esquemáticamente por una esfera. Solo que
en este caso la anomalía sería negativa ya que, generalmente, los domos de sal
son menos densos que los terrenos que les rodean.

Cilindro horizontal.- Sea un cilindro horizontal de radio R, enterrado a la


profundidad z sobre la masa unidad colocada en el punto P(x,0). Es equivalente al
que resultaría si tuviese una masa concentrada en el eje .

2mz z
g G  2GR 2
z

r2 z2 x2

2
donde: m = R  masa por unidad de longitud.
g
z 1

2
g x
x

1  

 z 

En función de (x/z) obtenemos una curva parecida a la anterior, un poco más


achatada. Calculando la distancia entre los puntos de inflexión, esta resulta ser
1.154 veces la profundidad del eje del cilindro. Un anticlinal alargado
horizontalmente puede ser representado por esta configuración. Capa infinita
horizontal de espesor constante.

Puede demostrarse que su efecto es gz  2Gh , lo que resulta independiente de


la profundidad de la capa. (Esta formula se aplica en la corrección de Bouguer)
Capa finita horizontal de espesor h o efecto de falla.

 x 

g z  2Gh  arctg 

 2 z 
Para x = 0, gt   g 

Esto significa que la anomalía sería la mitad de la correspondiente al caso


[Link] efecto de una falla será reflejado en el mapa de anomalías por el
estrechamiento rápido de las curvas isoanómalas.

[Link]. Determinación del regional en un perfil

Se grafica la anomalía de Bouguer debida a la


sección y sobre esta se interpola la curva
regional.
NOTA: El Campo Regional coincide con la
Anomalía de Bouguer en los puntos con regional estimado
afloramiento de rocas. De modo que si conoce la
posición de la roca basal en varios puntos a lo
largo del perfil, el tratamiento del regional será roca (afloramiento)
más exacto.
Si tiene información geológica por sondajes que Anomalía Bouguer
determinan la posición de la roca basal en un

punto, se puede calcular el valor regional en dicho punto a través de la relación:

Regional (en sondaje) = Anomalía de Bouguer (en sondaje) + 2G z

donde:

 : contraste de densidad entre el material de relleno y la roca encajante en


3
gr/cm .

z : espesor de recubrimiento hasta la roca basal en metros.

G: Constante de gravitación universal.


Si conocemos el valor del regional en un sondaje y en otros puntos con
afloramiento de roca, podremos definir la pendiente de éste y con ello trazar la
curva que representará el campo regional para el perfil.

2.2. Métodos magnético

2.2.1. Generalidades sobre el magnetismo

La tierra se comporta como un gran imán y la dirección e intensidad del campo


magnético varía considerablemente de un punto a otro de la superficie de la
misma. Si se sitúa una pequeña aguja imantada, capaz de girar alrededor de un
eje vertical, en un lugar cualquiera no señalará en general hacia el polo Norte
geográfico, sino hacia otro punto llamado polo Norte magnético. La dirección hacia
donde señala la brújula, define el meridiano magnético en el punto: el ángulo
formado por éste y el meridiano geográfico se llama declinación y se representa
por la letra d. La declinación varía de un punto a otro de la superficie terrestre. A la
componente horizontal de la fuerza del campo magnético terrestre en cualquier
punto se la designa por la letra h. El valor del campo en general se expresa en
oersteds y es como veremos más adelante, la fuerza medida en dinas, ejercida
sobre la unidad de polo magnético en el punto dado.

Los dos polos de un imán presentan comportamientos opuestos. Existe una


analogía, entre los polos magnéticos y las cargas eléctricas. El polo magnético
que señala hacia el Norte terrestre (es atraido por él), se denomina polo norte o
polo positivo y el que señala hacia el Sur terrestre, polo sur o polo negativo.
Estos dos polos no se distinguen en nada al actuar sobre el hierro no magnético,
así como las cargas eléctricas no se distinguen al actuar sobre cuerpos
descargados.

Un imán, es pues un dípolo magnético a semejanza de los dipolos eléctricos, pero


existe una diferencia radical entre ellos, en un dípolo eléctrico podemos separar
las dos cargas que forman el dípolo, obteniéndose así dos cargas libres, mientras
que si un imán lo dividimos, obtenemos dos imanes cuyos polos tienen la misma
intensidad que los polos del imán primitivo. No existe por lo tanto, polos
magnéticos libres.

Designamos por p la intensidad del polo magnético. Las intensidades de los dos
polos de un imán son siempre iguales en magnitud, pero de signo contrario: polo
norte ( + p ) y polo sur ( - p ).

Si se considera un imán de gran longitud y muy fino en las proximidades de uno


de uno de los polos, se puede despreciar la influencia del otro. Tenemos,
aproximadamente, como un polo magnético aislado.

2.2.2. Ley de las acciones magnéticas

Coulomb observó que, con barras magnéticas muy largas (aproximación del
polo magnético aislado), la fuerza de repulsión entre los extremos del mismo
signo, es inversamente proporcional a los cuadrados de las distancias. Llegó así a
la conclusión de que se trataba de fuerzas newtonianas, y por tanto, se les pueden
aplicar las leyes que rigen estas clases de fuerzas:

pp´
FC en expresión escalar o

r2

pp´
FC r3r en expresión vectorial

Donde:
p y p´ = intensidad magnética de cada polo.

r = distancia entre los dos.

r
= vector de posición del polo sobre el que actúa la fuerza F , respecto del otro.

C= Constante de proporcionalidad que depende del medio en el que estén


situados los polos.
2.2.3. Unidad de intensidad magnética

Se define como unidad de intensidad magnética (unidad de polo), la


intensidad de dos polos magnéticos que a la distancia de 1 cm, se repelen con la
fuerza de 1 dina.

2.2.4. Campo magnético

Si consideramos un polo magnético aislado, de intensidad p, se define


como campo magnético H , la fuerza que experimenta un polo magnético positivo
de intensidad unidad, debido a la presencia en el espacio del polo magnético
considerado. Su expresión será por tanto:

p
H  H  p r
2
r r3
en expresión escalar ó vectorial.

El vector r , serán en este caso el vector de posición del punto del espacio
en el que estudiamos el campo magnético respecto al punto ocupado por el polo p
origen del campo. Por tanto, si se trata de un polo positivo, los vectores de campo
son divergentes y si es negativo, convergentes.

La unidad de campo magnético, en el sistema c.g.s.e.m. es el Oersteds que


se define como el campo creado por un polo magnético de intensidad unidad a 1
cm de distancia.

-5
En Geofísica se utiliza como unidad el Gamma, γ , que corresponde a 10
Oersteds, ya que el Oersteds resulta una unidad demasiado grande para
Geofísica, donde las anomalías que encontramos son corrientemente de unas
decenas de γ .

Se presentan tres fenómenos magnéticos en la material someterla a un campo


exterior:

a) Diamagnetismo. Que consiste en una variación del radio y de la velocidad de


giro de las cargas de los átomos, con lo que varía el momento magnético de éstos,
de forma que esta variación se opone al campo magnético exterior de acuerdo a la
le de Lenz. Este fenómeno lo presentan todos los átomos, pero solo se aprecia en
aquellos en los que el número de electrones es grande y están dispuestos con una
simetría tal, que el momento magnético del átomo es nulo.

Las sustancias que presentan diamagnetismo se llaman diamagnéticas y son


por ejemplo, la sal común, el bismuto, la anhidrita, etc. El resultado de este
diamagnetismo es que el campo magnético es menor en el interior de estos
cuerpos y tenemos k (susceptibilidad) <0.

Como este fenómeno es interior al átomo, es independiente de la agitación


térmica de éstos, y por lo tanto, de la temperatura del cuerpo.

b) Paramagnetismo. Se presenta en las sustancias, en las que el momento,


magnético del átomo no es nulo, dada una simetría insuficiente en la disposición
de sus orbitales. Pero los momentos magnéticos de los átomos están en todas las
direcciones, con lo que estas sustancias aparecen como no magnéticas. Pero en
presencia de un campo exterior, se ordenan de forma que refuerzan la acción de
este y presenta una susceptibilidad k>0.

El fenómeno del Paramagnetismo depende de la agitación térmica de las


moléculas y por lo tanto de la temperatura. El coeficiente de imantación x, que es

k
igual a siendo k la susceptibilidad y la d la densidad, resulta ser para las
d

x 
T-
sustancias paramagnéticas:
Donde c es una constante de proporcionalidad y θ la temperatura de Curie, a la
que deja de presentarse al paramagnetismo y es constante para cada tipo de
sustancia.

Como ejemplo de sustancias paramagnéticas se pueden citar el perclorato y


sulfato de hierro, platino, aluminio, etc.

c) Ferromagnetismo. Que se presenta solo en el estado sólido en el que las


fuerzas interatómicas son lo suficientemente grande como para producir un
paralelismo de los momentos atómicos de un conjunto de átomos próximos, que
forman los llamados recintos de Weiss. Son estos recintos los que ordenan al
someterlos a un campo exterior, en un fenómeno similar al paramagnetismo.

2.2.5. Medición de la susceptibilidad magnética de las rocas

La medición de la susceptibilidad magnética de las rocas puede hacerse en


el laboratorio con un medidor de susceptibilidades que, en esencia, consiste en un
puente de inductancia. El aparato se ajusta a cero antes de efectuar la medición;
cuando se pone la muestra pulverizada de la roca cuya susceptibilidad se quiere
medir, se introduce una inductancia cuyo valor se puede medir por el método de
cero, llevando el galvanómetro indicador a las condiciones primeras de equilibrio.

En la literatura geofísica se describen también otros sistemas para medir la


susceptibilidad de las rocas in situ.

Los valores de la susceptibilidad se miden en unidades cegesimales.

Recordando que:

J dM 1
kH  dv  H

La unidad cegesimal de susceptibilidad se define como la de un cuerpo cuya


unidad de volumen adquiere, en un campo unitario, un momento magnético
unidad.

[Link]. Tabla de susceptibilidad de las rocas

En esta tabla se indican los valores de la susceptibilidad magnética de


algunas rocas y minerales, tomados de la tabla de Heiland y Birch. La mayoría de
estas susceptibilidades se han medido en un campo magnético algo mayor que el
terrestre. Los valores se expresan en unidades cegesimales.

-6
. . .Magnetita .0,3 a 0,8 Oersted . . . . .Arenista.... 16,8 ·10 Oersted
-6
Limenita . . . . .0,135 Oersted Dolomía. . . . . . . . . 14 ·10 Oersted

Pirrotina . . . . .0,125 Oersted En el vacío . . . . . . . 0 Oersted


-6
Franklinita . . . 0.036 Oersted Calcita pura. . . . . . -1 ·10 Oersted
-6 -6
Granito. . . . . .(28 a 2.700) · 10 Sal común . . . . . . . -1,3 ·10 Oersted
-6
Oersted Grafito . . . . . . . . . . -10 · 10 Oersted
-6
Diorita. . . . . .46,8 a 10

-6
Gabro. . . . .. .(68,1 a 2,370) ·10

Oersted

Instrumentos empleados en la prospección magnética

 Balanza magnética de Schmidt



 Inductor terrestre

 Magnetómetro de Torsión

2.2.6. Magnetometría

2.3. Sismología

2.3.1. Adquisición de datos sísmicos

2.3.2. Exploración sísmica marina


Fig. 1 Esquema de un barco

1.- Brújula

2.- Puente de mando con controles vista a proa ya popa


3.- Cuarto de maquinas

4.- Oficina de trabajo

5.- Almacén de refacciones

6.- lntercomunicación entre puentes- posicionamiento- sala de graficadores

7.- Área de muestreo

8.- Cocina

.- Comedor

10.- Dormitorios (mínimo 20 )

11.- Baños (mínimo 2)

12.- Sala de descanso

13.- La vanderia

Procedimientos:

Las bases teóricas de la geofísica de exploración están relacionadas


intrínsecamente con las propiedades físicas del medio, con la propagación y
reflexión del sonido, en las diferentes capas de los medios geológicos
.

Figura 2.- xxxxxxxxxxxxxxxxxxxxxxxxx


Figura 3.- xxxxxxxxxxxxxxxxxxxxxxxxx

Cuando un material de una determinada sustancia se ve afectado por uno de


estos movimientos ondulantes compresionales y esta sustancia se encuentra en
contacto con otra diferente y por ende de distinta velocidad de propagación, las
vibraciones de la primera sustancia se comunicará a la segunda, pero parte de
estas vibraciones rebotarán en el contacto de las dos sustancias para regresar al
origen de las mismas, creándose de esta forma el fenómeno de reflexión.

Los fenómenos de reflexión siempre ocurren con cualquier cambio de velocidad de


propagación, pero las reflexiones más fuertes ocurren cuando las diferencias entre
velocidades son muy contrastantes, los rangos de velocidades normales varían
entre los 1500 y 4000 m/seg, correspondiendo la primera al agua de mar y la
segunda a rocas de alta densidad. Por lo que el conocimiento completo del suelo y
subsuelo marinos requieren del uso conjunto de tres tipos de equipos: los
destinados a investigar el tirante de agua, los del fondo marino y los que exploran
el subsuelo.

Todos los equipos son empleados simultáneamente a bordo de una embarcación


especialmente equipada, la cual de acuerdo a un recorrido programado en la zona
por estudiarse, permite llevar a cabo el levantamiento geofísico marino.
Posteriormente y mediante el auxilio de computadoras se lleva a cabo el
procesamiento de la información para que la interpretación geológica y geofísica
pueda plasmarse en planos y perfiles descriptivos del área que se estudia.

Con la conjugación de todos estos factores, se contempla la forma de los


levantamientos, la distancia entre líneas de exploración que habrá de recorrer la
embarcación, la dirección de las mismas, la distancia entre puntos coordenados
de referencia y la determinación de los parámetros y rangos de operación bajo los
cuales deberán funcionar los equipos geofísicos y los sistemas de navegación que
permitirán llevar al barco explorador en la ruta proyectada. Como ejemplos de los
tipos de levantamientos geofísicos que puedan efectuarse se pueden citar los
estudios de selección de trayectorias para tuberías submarinas, obras portuarias,
retículas de exploración para instalación de plataformas petroleras, localización de
tuberías y objetos hechos por el hombre, los estudios geofísicos se llevan a cabo
por líneas de exploración programadas en la etapa de planeación.

Estas líneas están constituidas por una alineación de puntos igualmente


separados cuyas coordenadas son determinadas de antemano, conociéndose a
estos puntos como puntos de referencia o puntos de tiro. Con todos los sistemas
operando a bordo simultáneamente la embarcación es guiada por los sistemas de
navegación a cada uno de los puntos de referencia programados, efectuándose
una marca simultánea en todos los registros de los equipos, cuando la distancia a
los puntos mencionados es mínima, dichos sistemas de navegación operan en
rangos de 200 a 3100 MHz.

[Link]. Equipo Geofísico


Todos los sistemas están diseñados en base a la teoría electromagnética de
transmisión de ondas y operan en la parte del espectro comprendida entre 20 Hz y
3100 Mhz.. Actualmente se cuenta con una gran variedad de equipos geofísicos
de alta resolución para estudios marinos, esto equipos se agrupan en cinco
categorías básicas; las primeras cuatro operan bajo el principio de reflexión
acústica (sísmico) y el último bajo un sistema magnético, como se muestra en las
tablas 1 y 2.

[Link].1. Principios Básicos de Operación

La onda compresional que se genera viaja a partir de la fuente hacia todas


direcciones y es parcialmente reflejada hacia la superficie y parcialmente
refractada o transmitida hacia bajo del fondo marino a cada interface acústica más
profunda. Tal contraste depende de la diferencia acústica de los materiales en
general (en función de la densidad y sus propiedades elásticas), según se
incremente el contraste, la cantidad de energía reflejada, también se incrementa.

El límite de penetración se alcanza cuando la energía acústica reflejada es


demasiado débil para ser detectada, el límite real varía considerablemente de área
en área y depende parcialmente del número y reflexibidad de las interfaces
acústicas implicadas.

Las ondas de las interfaces poco profundas o cercanas se reciben más rápido
que aquellas de interfaces distantes o profundas, dado que el barco se mueve a lo
largo de las líneas de investigación, a intervalos fijos se reciben series de
reflexiones los cuales son graficados para formar un perfil continuo, donde el eje
vertical representa el tiempo de viaje acústico y el horizontal representa la
distancia de la interfase, en la practica, la velocidad de la onda sísmica varia
debido a que esta en función de las propiedades acústicas del medio de
propagación, generalmente se incrementa según se aumente la densidad,
usualmente se acrecenta con la profundidad. Un sistema geofísico no acústico
usado en la investigación es el magnetómetro marino el cual detecta la intensidad
total del campo magnético de la tierra, es sensible a variaciones locales en la
intensidad del campo magnético; así como a anomalías magnéticas que lo hacen
útil para la localización de tuberías, náufragos y otros objetos metal-ferrosos en el
fondo marino.
Tabla 1.- Sistema Acústico

Tabla 2.- Sistema Magnético


El procesado de esta información en la computadora resalta y mejora las
reflexiones más profundas, suprime los múltiples y filtra el ruido. El principio físico
de las técnicas de perfilación somera y profunda es el mismo que el de la
ecosonda, o sea, el de la reflexión, de la energía emitida, por las interfaces que
separan a dos medios de diferente impedancia acústica., excepto que la señal
tiene una frecuencia, una energía y una longitud de onda que le permita una
mayor penetración.

La cantidad de energía que penetra más allá del piso marino es función de la
impedancia acústica de los medios que separan la interfase de modo que la
penetración total dependerá del número de reflectores, de la calidad acústica del
medio, de la energía de la señal emitida, de la frecuencia, de la longitud y la
amplitud del pulso, si conocemos la velocidad de propagación de la energía a
través de los diferentes medios que ha producido la reflexión, la medida del tiempo
de propagación nos permitirá conocer la profundidad y espesor de cada estrato,
de manera que si tenemos una cantidad de puntos de medición suficientemente
densa, sobre una región, se puede construir un perfil geológico o un plano de
profundidad para todos y cada uno de los horizontes reflectores.

En la practica la velocidad de la onda sísmica es compleja, depende de las


propiedades acústicas del medio y aumentan con la densidad y la profundidad, por
esta razón las profundidades calculadas a partir del sismograma para diferentes
estratos y accidentes, generalmente tiene que ser calibradas con la información
obtenida en sondeos geotécnicos en el área de estudio, en resumen, la
penetración de la señal es función inversa de la resolución, la resolución es
función directa de la frecuencia y la atenuación de la energía acústica se
incrementan en razón directa del cuadrado de la frecuencia, la máxima capacidad
de resolución es del orden de 1/4 de la longitud de onda y es inversamente
proporcional a esta ultima.

El procedimiento de operación de los equipos requieren de un barco con las


dimensiones y características ya descritas anteriormente una vez posicionado en
la red de líneas de recorrido, debe navegar en un tirante mínimo de agua que le
permita maniobrar libremente a una velocidad media de 3 nudos (1 nudo = 0.5
m/seg), misma que si se disminuye permite obtener mayor detalle. Las unidades
de transmisión y recepción de todos los equipos son remolcadas por el barco, que
no se detiene para realizar la operación, el conjunto de todos estos equipos son: la
ecosonda, el sonar de barrido lateral o mapeador del piso marino y los perfiladores
somero y profundo.

Calibración del equipo. La calibración del equipo es una de las partes más
importantes del estudio y de los resultados que se obtengan, por lo cual durante
todo el estudio se vigila constantemente la calibración de los equipos.

Para probar el equipo, se utilizan sondas graduadas con una placa metálica en el
extremo, la cual refleja las ondas emitidas por la ecosonda y mide la profundidad
de la placa, que debe coincidir con la grabación de la sonda o estar dentro del
rango de la precisión que establezca el fabricante, esta prueba se hace dos veces
en un mismo día a diferentes horas y para cada serie de mediciones, se establece
una curva de calibración y se dan las correcciones en la velocidad del sonido para
cada profundidad, otra prueba para este equipo, se hace en el momento en que se
trabaja con el oleaje máximo (olas de 1.2 a 1.5 m de altura ), donde la grabación
de la grafica no debe perder continuidad.

Una calibración general para el resto del equipo, se hace mediante las mediciones
de la ecosonda, ya que, de acuerdo a la precisión de cada equipo, en todos se
puede medir el tirante de agua y debe coincidir con lo indicado por la ecosonda.

El funcionamiento de estos equipos es similar y se basa en el principio de la


emisión de pulsos acústicos de muy alta frecuencia que al chocar con cualquier
objeto o superficie se refleja en forma de eco para ser recibido por el transductor
que lo emitió, si se mide electrónicamente el tiempo entre la emisión y recepción
(se supone constante su velocidad en el agua) se está en condiciones de conocer
el tirante de agua.

Las frecuencias normales de operación de la ecosonda y el sonar lateral son de


200 y 105 KHz respectivamente.
[Link].2. Sistema Perfilador (Batimetrico)

La ecosonda es el equipo básico en el procedimiento para la medición y


exploración del tirante de agua, existen varios tipos de ellas, clasificándose,
primeramente como de navegación, de exploración, y de pesca, para los objetivos
de estudio en Campeche, los equipos se eligieron en atención a la presición fija en
:I: 0.30 m. y la profundidad especificada como máximo 100 m. requiriéndose
ecosondas de investigación.

Las cuales pueden ser de tipo digital con doble frecuencia de operación o de tipo
analógico; éstos emiten rayos de ondas acústicas, con un ángulo entre 5 y lO
grados respecto a la vertical a través de cristales piezoeléctricos y estos mismos
cristales reciben las ondas reflejadas del fondo o de cualquier otro objeto;
posteriormente, al considerar constante la velocidad del sonido en el agua (1500
m/seg), el equipo calcula el tirante de agua (las escalas horizontal y vertical de
dicha grafica generalmente son diferentes ). La profundidad real estará dada por la
siguiente relación:

H = h :I: hc -hm

Donde:

H es la profundidad real

h es la profundidad medida en la ecosonda

hc es la profundidad corregida por calibración

hm es la altura por marea

La profundidad h primero se ajusta por el movimiento vertical que tiene la


ecosonda debido al oleaje; con dicha corrección, se puede distinguir un fondo
uniforme y uno irregular. La profundidad hc, se obtiene por corrección de la
velocidad del sonido, para lo cual se utiliza un celerímetro, el monograma de
Grosso (ver figura IV.2) también se puede usar la siguiente fórmula empírica :

V = 1410 + (4.21 t- 0.037 t 2) + 1.105 s + 0.018d


Figura: 4.- Velocidad del sonido en el agua (sobre la superficie),en

términos de la temperatura y la salinidad

(De acuerdo al diagrama de grosso)

(Depto. de Ingeniería 1990 UNAM)


Donde:

s es la salinidad del agua en ppm

t es la temperatura del agua en grados centigrados

d es la profundidad en metros

La altura hm es una corrección por efecto de la marea, para lo cual se considero el


nivel de bajo del mar proporcionado por la UNAM; esta corrección influye mucho
en regiones con amplitudes de mareas grandes. Las frecuencias características de
operación de los aparatos usados son de 10 a 30 KHz, para la primera
corresponde una longitud de onda de 15 cm y para la segunda de 5 cm,
aproximadamente. La señal se radia hacia abajo a través del transductor,
simultáneamente el mecanismo medidor de tiempo comienza a funcionar y se
detiene en el instante en que la señal reflejada del piso oceánico es recibida en el
mismo transductor, determinándose. La altura o tirante de agua ya que también se
puede obtener la profundidad real de la siguiente expresión matemática :

1 Vst
Donde : h 
2

h = es la profundidad

V = la velocidad del sonido en el agua

t = tiempo de viaje

La velocidad del sonido en el agua es algo muy importante y debe ser conocida
para poder calcular la profundidad real del tirante de agua. Esta velocidad
depende de la temperatura, salinidad y presión, la velocidad del sonido en el agua
varía en un rango que va de los 1400 a los 1570 m/seg y su valor específico
depende de la región en donde se haga la medición, el movimiento del barco y la
topografía del fondo influyen en la exactitud de la medición del tirante de agua
(figura IV .3) ya que pequeñas señales del sonido conducen a ecos no solo del
piso oceánico si no directamente bajo el barco y también de otros puntos del fondo
marino.

El movimiento de cabeceo del barco alarga o acorta la profundidad, además de


éste, hay que corregir la misma por altura de marea y por la posición del sensor
con respecto a la antena del barco u otros ajustes de operación.
Figura 5.- Movimiento de cabeceo del barco
Representa un registro característico del ecosonda en la que la escala horizontal
esta dada por la velocidad del barco, la que a su vez depende del tipo de trabajo
que se éste realizando, para trabajos geotécnicos, la velocidad media es de
aproximadamente 6 nudos.

La escala vertical del registro representa la profundidad, comúnmente medida en


brazas (1 braza = 1.83 m) sin embargo también se puede medir en metros, esta
escala es generalmente exagerada con la horizontal, en el caso de fondos
conocidos como en algunas áreas de la sonda, se utilizan equipos de mayor
precisión como el caso del pinger que es una ecosonda de precisión compuesta
por un dispositivo suspendido a una profundidad cercana al piso oceánico, que
emite pulsos o pings a intervalos precisos de un segundo aproximadamente, de
estos pulsos uno va directamente y el otro después de ser reflejado en el piso
hacia la ecosonda fija en el barco, la diferencia en tiempo de arribo entre cada una
de las señales permite determinar la profundidad del pinger sobre el piso marino,
el transductor del pinger combina una frecuencia variable, un transmisor de alto
poder con ancho de banda ajustable y un receptor de ganancias de variación del
tiempo, formando así un sistema múltiple, que se puede utilizar como ecosonda,
como perfilador sísmico y como sistema de navegación. La frecuencia y potencia
de salida puede ajustarse en forma continua de; a 12 KHz y de O a 10 Kw
respectivamente.

El levantamiento batimétrico en la ecosonda se ubica fijando las líneas de


exploración con los procedimientos descritos en el capitulo de posicionamiento, en
la Sonda de Campeche los recorridos se realizan por sectores de 2 X 2
Kilometros, en líneas norte-sur a cada 150 m y oriente-poniente a cada 500 m.
(figura IV .5). Una vez posicionado el barco, define los registro por línea, los que
integrados nos proporcionan la batimetría del bloque de 2 x 2 Km, para el caso de
plataformas. Cuando se trata de levantamientos para el proyecto de ductos, los
sectores se integran en todo este recorrido de la tubería.

Como el objetivo de medir tirantes de agua se complementan con el detector de


burbujas de gas en todo el sector, en algunos casos la batimetría se complementa
con un dispositivo llamado transductor sintonizado y con el sonar de barrido
lateral. El transductor sintonizado consta de un multisensor y un registrador grafico
de 7 KHz, produce un pulso muy corto de alta energía de uno o dos ciclos de
frecuencia de señal, la señal que se refleja se registra con ganancia constante o
con ganancia de tiempo variable, lo notable de este aparato es su capacidad para
la fácil detección de burbujas en la columna de agua, algunas veces el transductor
sintonizado sustituye a la ecosonda y al perfilador somero, esto depende del
objetivo del estudio. ya que los proyectos se desarrollan en un área del piso
marino; ya que el complemento e integración de la información proporcionada por
la ecos onda es la determinación de la fisiografía y anomalías del piso marino, que
hemos establecido como segundo objetivo del estudio geofísico.

[Link].3. Sistema de Barrido del Fondo Marino

El sistema de mapeo del fondo marino de doble canal será usado para determinar
obstrucciones naturales o artificiales, variaciones del sedimento en el fondo y
fenómenos de batimetria entre líneas de rastreo. El equipo puede unirse a otro
como un sistema integrado de sonar de exploración lateral y un sub-superficial del
lecho marino el cual nos da perfiles, esta información se obtiene a través de la
interpretación de sombras e imagenes acústicas en escala corregida de la
batimétria. En su interior el equipo trae sensores capaces de almacenar la
información que después seran extraidos y llevados a una computadora. Es
posible una cobertura de 100 % usando el rango de escala recomendada de 150
m por cada canal.
Figura 6.- xxxxxxxxxxxxxxxx
Sistema Sísmico de Alta Resolución .-
1.- Perfilador de subfondo de 3.5 KHZ

2.- Perfilador Somero

3.- Perfilador Profundo análogo de 4 kj a 8 kj

1.- Perfilador de subprofundo .-El equipo cuenta con una compensación de


oleaje y potencia de salida de 3.5KHZ, es un perfilador de poca penetración, y
provee detalles de subfondo de los primeros 10 m a 30 m por debajo del fondo
marino, también proporciona información sobre las escarpas del fondo marino,
fallas cercanas a la superficie y zonas de hundimiento, y muestra los cambios de
estratos y sedimentos del fondo, imprimiendo los registros en papel con un disparo
de 0.25 segundos de velocidad. En la figura IV .8 se muestra un registro típico.

2.- Perfilador Somero .-El sistema de perfilador somero con penetración


intermedia de 300 joules de salida analógico será provisto para mostrar
condiciones geológicas, cambios estratos y sedimentos del fondo a penetraciones
tipicas de 30 m a 60 m por debajo del fondo marino, algunos tipos de perfiladores
someros se usan simultaneamente para realizar perfiles batimétricos, se
proporcionan registros con espacio de marcación de 0.25 segundos, un
miniperfilador profundo de 100 joules también estará disponible en el barco como
una alternativa al perfilador somero. El sistema de penetración intermedia del
perfilador somero (mini perfilador profundo ), el cual operará simultaneamente con
el perfilador profundo de 8 kj, turnando la ventana de grabación y disparo de cada
sistema. Los mecanismos de la fuente acústica para estos sistemas puede ser :
piezoelectricos, electromecánicos, explosión de gas, disparos de agua y el
denominado Sparker en la figura IV .9 se muestra un registro.

3.- Perfilador Profundo análogo .-El equipo cuenta con compensación de


oleaje y potencia de 4 kj a 8 kj este sistema de alta penetración incluye una fuente
de sonido, la potencia de salida de este sistema es de 400 a 8000 joules, y
proporciona los perfiles en los cuales se pueden observar las condiciones
geológicas, cambios de estratos y sedimentos del fondo marino, a profundidades
de 300 m a 600 m por debajo del fondo marino, dependiendo de las condiciones
locales, se deberá imprimir el registro con un segundo intervalo en papel termico.
Con este tipo se obtiene información menos detallada que con los perfiladores
someros. Los mecanismos de la fuente acústica para estos sistemas pueden ser
del siguiente tipo: Sparker, disparo de aire, disparo de agua y explosiones de gas.
En la figura IV .lO se muestra un registro típico Magnetometro Marino (sistema
sensible al campo magnetico) .-Continuando con el equipo utilizado para la
detección de objetos en el lecho marino mencionaremos que el magnetometro de
protones marino es un sistema que ira remolcado detrás del barco. En una
grabadora de doble canal. Se muestran los valores de campo magneticos gamma
de 1 m a lO m. Los registro se pueden valuar para determinar la presencia de
objetos ferrosos tales como tuberías, cabezales de pozos, y escombros marinos
en la figura IV .11 vemos un esquema de un sondeo
SISTEMAS

1.- Posicionamiento

2.- Medidor de Profundidad

3.- Mapeo del Fondo Marino

4.- Perfilador Somero

5.- Perfilador Profundo

6.- Transductor Sintonizado (pinger)

[Link].4. Hidrofonos

Los hidrofonos del sistema detector son dispositivos piezoeléctricos y funcionan al


recibir las señales reflejadas que presionan con diferentes intensidades,
provocando señales eléctricas que son filtradas y graficadas, y algunas veces
también grabadas en cintas magneticas, los hidrofonos generalmente se arreglan
en pares, de modo que sus salidas eléctricas se cancelan para una aceleración de
traslación, pero se suman para un pulso de presión.
Debido a que los hidrofonos tienen una impedancia eléctrica alta, frecuentemente
se usan reguladores de impedancia para cada grupo. Los hidrofonos se
encuentran montados en una especie de larga flauta (streamer) que va atada al
barco de prospección a una profundidad que varía entre los lO y 20 m, de 8 a 20
hidrofonos espaciados a intervalos de algunos centimetros se conectan en serie
de manera que se pueda sumar el voltaje generado por ellos para formar un solo
detector (figura IV.12), la información de la perfilación se graba en un papel
electrosensible. Un alambre helicoidal rodea al cilindro sobre el cual se encuentra
el papel, y una corriente eléctrica que va de los alambre helicoidales a través del
papel hasta la placa del cilindro produce la imagen sobre el papel, algunas veces
las señales acústicas se graban sobre cintas magneticas, especialmente cuando
las flautas con arreglos de hidrofonos son multicanal.

[Link].5. Selección del Equipo

La selección del equipo especifico depende de la aplicación, profundidad del agua,


costo, disponibilidad de equipo y de otras consideraciones sin embargo, los
estudios geofísicos más completos son aquellos que incluyen la aplicación de los
cincos tipos de sistemas de adquisición de datos, la eliminación de algún tipo de
equipo, con lleva a hacer suposiciones acerca de las condiciones existentes de un
sitio. En la práctica, todos los sistemas seleccionados usualmente se operan
simultáneamente para reducir el tiempo de exploración geofísica. Dependiendo del
tipo de sistema, los transductores sísmicos (fuente acústica y receptor), pueden
estar localizados en un solo sitio o desplegados separadamente, y pueden estar
fijos bajo el casco de la embarcación de reconocimiento o independientes
remolcados. También dependiendo del sistema los datos pueden ser registrados
sobre papel en tiempo real, almacenados en cinta magnética para ser procesados
posteriormente, las aplicaciones generales y descripción de algunos equipos
geofísicos de alta resolución se resumen en las tablas IV .2 y IV .3

[Link]. Presentación de la Información

La elaboración de los planos o mapas de un estudio como este, se basa en la


interpretación de los registros de los datos levantados en campo y son el
complemento del reporte definitivo. A continuación se describen los diferentes
tipos de planos que se deberán proporcionar.
[Link].1. Plano de posicionamiento

Este plano indicara en coordenadas del sistema Universal Tranversa de Mercator,


la ubicación de cada uno de los puntos de tiro. La precisión y claridad de este
plano es fundamental ya que servira de base para las siguientes etapas.
[Link].2. Plano Batimétrico

El plano batimétrico corregido por variación de mareas, debe estar referido al


plano de posicionamiento, contener las curvas batimétricas con equidistancias de
1 metro en relación al nivel del mar.
En zonas con pendientes menores de I por ciento, la equidistancia de las curvas
se reducirá a 50 cm. La precisión con respecto a la profundidad verdadera no será
mayor de 0.5 por ciento (plano IV .1 ) Plano de Fallas de la Zona .-Este plano nos
proporciona a detalle, la localización de las fallas geológicas principales, indicando
cuales son normales y cuales son inversas y si se encuentran activas o pasivas
(plano IV.2)

Plano de Riesgos Potenciales .-En este plano se muestra las posibles zonas que
puedan presentar problemas para la perforación de pozos e instalación de las
estructuras fijas; estos problemas pueden ser: burbujas saturadas de gas, fallas
geológicas de consideración, sedimentos inestables, emerción de volcanes,
topografía difícil etc. (plano IV.3). Los planos de anomalías o riesgos potenciales
se elaboran integrando la información del sonar de barrido lateral con planos de
los accidentes estructurales y cuerpos anómalos detectados con los perfiladores

.En estos planos se representan todos los eventos naturales (geológicos) y


artificiales (tendido de lineas, barcos hundidos) que representan un impedimento o
peligro para el desarrollo del proyecto y construcción de obras en el área de
estudio e interés, sin importar que estos eventos se registren sobre el piso marino
(afloramientos coralinos, zonas de alta pendiente o deslave, tuberías, etc.) o
enterrados a cierta profundidad (paleocanales, bolsas de gas, fallas, formaciones
coralinas, etc.).

Este mapa es de gran importancia para la ingeniería, ya que permite discriminar o


clasificar zonas para el trabajo de tendido de líneas nuevas, perforación de pozos
y cimentación de plataformas, es necesario agregar que todos estos planos van
acompañados de un informe escrito que describe y discute la información y
presenta, además diferentes conclusiones y recomendaciones, en concordancia
con el objetivo del estudio.

Estructuras geológicas de riesgo potencial .-Los considerados como riesgos


potenciales son los arrecifes aflorantes o encerrados por presentar una condición
de apoyo muy variable en sus propiedades y diferentes a las arcillas blandas que
constituyen la uniformidad en la sonda. Los cuerpos de arena cementados,
constituyen otro de los detalles considerados como riesgos; Las fallas geológicas
por sus desplazamiento o actividad probable son también otras de las estructuras
y se presentan con frecuencia, generalmente estos accidentes tienen asociados
algunas zonas de plegamientos; y su proyección es en la zona de campeche es
hasta de 6 Km de longitud. Otras estructuras reportadas como riesgos por su
posible inestabilidad local son los paleocanales y bolsas de gas que son
importantes. Existen también crateres que constituyen otros de los detalles de los
considerados como riesgos.

Las estructuras anteriores son las naturales, pero se deben considerar los
accidentes definidos por el hombre, como tuberías, barcos hundidos u otros
detalles que proboquen problemas al desarrollo de obras en el fondo marino
Zonificación del Perfil

[Link]. Analizado

A travez de secciones geológicas mostrar la continuidad de los estratos con sus


respectivas características como son: litología, porosidad, fracturas, pliegues y
fallas.
[Link].1. Plano de Isopacas

Las propiedades mecánicas de los suelos tales como el tamaño de grano, módulo
de poisson, densidad y porosidad, no puede ser cuantificadas mediante la
perfilación acústica, sin embargo se puede calificar estratos litológicamente
basandose en las características acústicas producidas por materiales similares
explorados y calibrados con anteriorida. En caso en que el reflector sísmico se
pueda seguir clara y fácilmente con métodos de perfilación, fuera de un área
completamente explorada y verificada, las propiedades litológicas y mecánicas de
los suelos pueden extrapolarse con poca probabilidad de error .
información obtenida con los perfiladores somero y profundo se presentan en
planos de isopacas someras y otros de isopacas profundas. En los primeros lo que
se configura son las líneas que unen a los puntos de igual espesor de la primera
unidad estratigráfica, es decir, la capa del suelo que va del piso oceánico al primer
horizonte reflector abajo de este y los cuerpos y estructuras geológicas (arrecifes y
fallas) contenidos en la unidad de suelos marinos (figura IV .4-a); en los segundos
se representa la misma información pero en este caso a mayor profundidad,
generalmente a las isopacas que van de la superficie del suelo oceánico hasta la
primera capa consolidada o el lecho rocoso (figura IV .4-b ) son las líneas que
definen las profundidades de la primera capa consolidada o lecho rocoso en el
área de estudio.

[Link].2. Plano estructural somero

Este plano presentara el horizonte reflector más significativo de las tendencias


estructurales de la zona y será elaborado mediante los datos de los perfiladores
somero y profundo Perfiles Interpretativos .-Una herramienta de, trabajo de gran
importancia por su objetividad la constituyen los perfiles interpretativos o
secciones, en las cuales se vacía la información más sobresaliente del
levantamiento geofísico. En estas secciones los diferentes ingenieros y técnicos
pueden observar sus objetivos de interés, como son: el tirante de agua, el espesor
y comportamiento de los materiales no consolidados, el comportamiento
estructural de un horizonte con capacidad de carga suficiente, fallas geológicas,
paleocanales y gas en depósitos someros, entre otros rasgo.

En general, toda la información obtenida se agrupa en un reporte técnico que


termina con las conclusiones de las condiciones del área estudiada y las
recomendaciones necesarias para trabajos subsecuentes de exploración
geotécnica, construcción de cimentaciones o de perforación petrolera.

En la sismología de reflexión es frecuente la aparición de reflexiones denominadas


multiples, que son fenómenos de repetición de eventos sísmicos provocados por
interfaces de medios de muy diferente densidad, como lo son las fronteras entre
agua -aire; agua -tierra etc. Cuando la onda se refleja una sola vez se dice que se
trata de una reflexión primaria, en cambio si después de ser reflejada bajo la
superficie donde se ha liberado la energía se refleja una o varias apareciendo con
la misma forma, se dice que es una reflexión

multiple como se muestra en la figura IV.13 .Este fenómeno puede ser muy simple
o muy complejo en si mismo o en su interpretación (figura IV.14). .La
interpretación de los multiples es muy importante para evitar la falsa identificación
y ubicación de estratificaciones, estructuras y accidentes del medio. El multiple
denominado fantasma ocurre cuando la señal es reflejada en primera instancia por
la interface que componen el techo del medio, donde se libera la energía, la que a
su vez refleja la misma señal hacia abajo, o sea hacia el horizonte inferior, del
medio donde una vez más es reflejada hacia arriba, donde finalmente es
detectada, con un tiempo de retardo respecto al tiempo utilizado por la reflexión
primaria, manifestando en el sismograma un tenue y falso horizonte a una
profundidad un poco mayor de la real, en la sismología marina los fantasmas son
importantes debido a que la superficie libre del agua es casi un reflector perfecto
de modo que la interferencia fantasma, en este tipo de operación, es más fuerte.
La probabilidad de aparición del fantasma se reduce a medida que la profundidad
del punto de liberación de energía es menor. Otro problema
común en la operación marina es la reverberación, que es un ruido indentificable
por su periodo de aparición constante, la cual es originada por la reflexión hacia
arriba y hacia abajo repetidamente en el tirante de agua, esta energía reverberada
es reforzada continuamente por la energía reflejada, en trabajos de prospección
con fines geotécnicos donde el interprete se aboca a las capas de los sedimentos
no consolidados que sobreyace a capas más consolidadas o a un lecho rocoso,
los múltiples no son un grave problema para la interpretación porque los suelos
marinos suelen ser suaves y no presentan fuertes contrastes de impedancia
acústica.
Las discordancias son reconocidas sísmicamente por la terminación lateral de los
reflectores que se traslapan, recubren, se plegan o truncan. La interpretación nos
provee de un bosquejo de sus principales rasgos que delinean las unidades
estratigraficas en cuestión y algunas veces el tipo de discordancia angular o de
contacto (figura IV .15) El primer tipo es el de más fácil de identificar ya que en él
los estratos recientes descansan sobre capas erosionadas, buzantes o plegadas,
el otro tipo tiene rasgos asociados con estructuras arrecifales y paleocanales.

[Link]. Tipos de Discordancias

Existen varios tipos de fallas (figura IV .16) y su detección en los sismogramas es


relativamente fácil, estas se manifiestan dentro del sismograma por corrimiento de
los reflectores, de una manera coherente y consecutiva, hacia arriba o hacia abajo
con respecto a el mismo(figuras IV .17) En cuanto a la determinación de la
estratigrafía, el método es tan bondadoso para este fin que basta observarse
cualquier sismograma (somero o profundo) para darse cuenta como el tipo de
señal sísmica nos va delimitando las diferentes unidades estratigraficas excepto el
cuidado que se debe tener con los multiples.

[Link].1. Gas

El gas metano comunmente encontrado en los cuerpos sedimentarios


generalmente están asociados con fallas y fracturas se encuentran disperso en la
columna sedimentaria o localizado en paquetes o lentes de arena. El efecto que la
presencia de gas tienen en los registros sísmicos es el incremento de la
atenuación acústica debido a la absorción de energía por las burbujas de gas, la
dispersión de la señal sísmica es debido al movimiento de las mismas burbujas,
en los sedimentos superficiales se atenúan tanto las propiedades acústicas del
medio que generalmente se presentan en el sismograma zonas desarticuladas o
turbias, el gas disperso afecta a las bajas frecuencias causando reflexiones
borrosas o inconsistentes, los paquetes aislados de gas se manifiestan en el
sismogramas por zonas de reflexión de gran amplitud llamadas puntos brillantes,
estas zonas están sujetas a un exceso de presión por los suelos circundantes lo
cual es causado por la expansión de las burbujas de gas dentro de los sedimentos
contenedores de gas. Uno de los afloramientos comunes a las plataformas
continentales recientes es el escape de gas hacia la superficie del mar, el cual
puede ser de origen reciente, como en el caso del gas biogenetico; o bien,
proveniente de depósitos de hidrocarburos, como parece ser común en la zona
marina petrolera de la bahía de Campeche.

La importancia de los afloramientos de gas dentro de la ingeniería radica en que,


en las zonas en que aflora, altera las propiedades mecánicas de los suelos; esto
se debe a que al mezclarse con los sedimentos, ocupan los vacíos del suelo,
desplazan el agua interticial y en ocasiones la presión del afloramiento hace fallar
al suelo, el cual se conserva así sobre el fondo del mar. La identificación de este
evento se efectúa con la ayuda del sonar de barrido lateral, la ecosonda y el
perfilador somero, siendo frecuente en estos registros la aparición de montes o
volcanes de Iodos y cráteres causados por la salida de gas hacia la columna de
agua.

[Link].2. Corales

Un obstáculo frecuente para las obras marítimas en algunos casos es la presencia


de corales sobre el fondo marino, especialmente en tendidos de ductos
submarinos e instalación de plataformas petroleras.
El origen del coral se debe a la presencia de colonias orgánicas de polipos que se
desarrollan generalmente en zonas rocosas y cuya acumulación da origen a
diversas formas arrecifales que en ocasiones afloran sobre el fondo y sobre la
superficie del mar. El problema más importante de los corales es que, aún cuando
tuviese una buena capacidad de soporte, no son bien conocidas sus propiedades
mecánicos y son un elemento frágil para construcción. La identificación de estas
estructuras orgánicas es posible gracias a la ecosonda, sonar de barrido lateral y
perfilador somero, los cuales permiten dimensionarlas y detectar su forma
geométrica.

Roca .-Dentro de la ingeniería de cimentaciones es importante el conocimiento de


la profundidad del basamento rocoso para el anclaje o cimentación de estructuras,
si la roca aflora puede ser descubierta con la ayuda del sonar de barrido lateral, el
cual permite determinar las características geometricas de los afloramientos
rocosos y localización en un plano de proyección ortogonal.

Eventos arrecífales .-En algunas ocasiones la geofísica puede prestar servicios a


la localización de objetos o construcciones existentes en el fondo del mar; para
ello, tanto el magnetómetro como el sonar de barrido lateral, resultan valiosas
herramientas para estos propósitos.
Subsuelo marino .-Para la instalación de plataformas autolevantables de
perforación y cimentaciones de plataformas de producción es de gran ayuda el
conocer la existencia de un horizonte del suelo marino que ofrezca la capacidad
de carga suficiente para soportar el peso de estas estructuras. Este estrato, por lo
general, lo constituye la primera capa consolid~da, que normalmente infrayece a
los Iodos del fondo del mar, la cual aparece claramente como reflector en los
registros del perfilador somero.

Los espesores de sedimentos representados en el plano de isopacas son los


correspondientes a los materiales sin consolidarse, los cuales se determinan
restando la profundidad de la primera capa consolidada a la profundidad del fondo
del mar y en el caso que no se conozca las profundidades y se conozcan
unicamente los tiempos de reflexión de ambos eventos, se aplica la ecuación :

{2-{1

e = V- K Ecuación 1

donde: e: espesor en metros.

v: velocidad sísmica de propagación en el fondo en metros por segundos

ti : tiempo de reflexión del fondo marino en segundos

t2 : tiempo de reflexión de la primera capa consolidada en segundos

K: factor geométrico de corrección adimensional


Arrecifes .-Algunas veces los arrecifes se encuentran sepultados por sedimentos
marinos más recientes y es necesario conocer sus dimensiones y profundidades
para considerarlos en el diseño de una cimentación.

Para tal efecto se debe iniciar por localizarlos en los registros de los perfiladores y
marcar su piso y su techo, con lo cual se lograra medir el tiempo que requirieron
las señales sísmicas en recorrer ese espacio y, si se conoce la velocidad de
transmisión de las ondas, puede aplicarse la ecuación 1 modificada de la siguiente
manera :

t -t

e = V p a -K Ecuación 2
2

donde: tp : tiempo de reflexión del piso del arrecife.

ta : tiempo de reflexión del techo del arrecife

Horizontes geológicos .-Cuando es necesario localizar estratos y sus fronteras


para que puedan ser empleados como cap capas de apoyo para pilotes de punta
de cimentaciones de plataformas petroleras, mediante métodos de interpretación
de los registros del perfilador profundo se seleccionan visualmente y sus
profundidades y espesores se calculan en base a la información disponible de la
ecuación de velocidades en algun pozo petrolero cercano cuyo comportamiento se
extrapola a las zonas de estudio.

La ecuación lineal de velocidades es la siguiente :

V = Vo + Cz Ecuación 3

donde: v: velocidad en metros por segundos a la profundidad z

Vo : velocidad inicial media en metros por segundos

C : constante igual a la pendiente ~ en seg-l

Z : profundidad en metros
como Voy C de la ecuación 3 son datos del pozo observado, con la ecuación 4 se
obtiene la

velocidad Vt a un tiempo de de reflexión t dado :

Va

E.,

4t = -0- cuaclon

1--

donde: Vt : velocidad en metros por segundo al tiempo t

t : tiempo de reflexión en segundos


Para obtener la profundidad Z o espesor de los horizontes geológicos medidos se
aplica finalmente la

ecuación 5 :

ZVt .t .,

5= -Ecuación

donde Z : profundidad en metros

Para identificar las características litológicas y verificar las diferentes unidades


estratigráficas es necesario auxiliarse de la información geológica y geotécnica
disponible del área en cuestión o de las áreas vecinas y de preferencia hay que
calibrar los datos obtenidos con los perfiladores sísmicos, con los resultados de
estudios a base de exploración geotécnica directa. El inferir las características
litológicas de los suelos, únicamente de datos sísmicos, puede conducir a graves
errores sobre todo cuando no se tiene una amplia experiencia en la interpretación .

Configuración de las estructuras .-Cuando los registros obtenidos en el mar han


sido recolectados con perfiladores analógicos se hace indispensable la
identificación y marcado de los verdaderos reflejos, ya que en este tipo de
registros existen otras señales que no deben tomarse como horizontes de
reflexión, puesto que son ruidos sísmicos : ecos de reflexión, ondas refractadas,
reverberaciones y otros; los verdaderos reflejos pueden ser inidentificables con
cierta experiencia del interprete. En el caso en que los registros que sean

producto de un procesamiento digital completo, el ruido puede removerse


mediante la aplicación de programas de correcciones estáticas y dinámicas,
deconvolución predictiva, filtrado numérico, autocorrelación, punto de reflejo
común, apilamiento y edición.
Selección .-Cuando se cuenta con datos geotécnicos provenientes de sondeos
mecánicos en áreas próximas e inclusive dentro de las exploradas, se puede
efectuar una correlación entre esos datos geotécnicos y los horizontes reflectores
de los registros del perfilador profundo, con objeto de seleccionar un horizonte
cuyas características mecánicas sean adecuadas para cimentar estructuras. Este
reflector se marca en todos los registros para llevar acabo su configuración
estructural. En los casos en que la información geotécnica no exista, se elige
aquel reflector, dentro de los rangos de penetración de pilotes de cimentaciones -
50 a 200 m. bajo el fondo marino, que ofrezca continuidad y un buen contraste en
amplitud de reflejo.
Basamento o capa consolidada .-Los datos referentes al basamento rocoso o
alguna capa resistente en especial son importantes para el diseño de la
cimentación en razón de que permite a los ingenieros proyectistas escoger el tipo
de pilotes en función de la profundidad a la que se encuentre un estrato que pueda
soportar las cargas que transmitan los pilotes. Si se logra correlacionar un
horizonte de reflexión con los datos de geotecnia, ello permite conocer la
profundidad esperada de penetración de los pilotes y la longitud de los mismos,
siempre y cuando el horizonte escogido ofrezca la suficiente capacidad para
soportar el pilote de punta.

Estrato blando superficial .-El estrato blando superficial esta constituido por arcillas
calcáreas, color gris y gris olivo, de consistencia blanda a muy blanda, con
fragmentos de concha y olor acido sulfúrico. El espesor de este estrato, tiende a
aumentar hacia el noroeste, donde de este estrato y los espesores medidos en
cada sondeo, se han delimitado en cuatro zonas con espesores de 0 a 1

[Link]. Métodos Directos

[Link].1. Trabajos Marítimos

Los trabajos marítimos se dividen básicamente en dos etapas, la primera consiste


en levantamiento de los registros geofísicos y la segunda, en el muestreo de los
sedimentos no consolidados.

En la primera etapa, para las plataformas, el barco recorre las rutas del crucero
trazadas dentro de un área de 2 X 2 Km, el centro de dicha área queda localizado
en el sitio propuesto para la instalación de la plataforma y las rutas del crucero se
fijan formando una retícula; normalmente, las líneas de crucero en un sentido
quedan espaciadas a cada 150 mts. y en el otro a cada 250 mts.

En el caso de las tuberías, las rutas de crucero principales son tres, las cuales
corren a todo lo largo de la trayectoria de la tubería y equidistante a cada 150 mts.
perpendicularmente a éstas, se trazan rutas de crucero a cada 500 mts., durante
el recorrido de las rutas del crucero se operan en forma simultánea los tres
diferentes tipos de equipos electrónicos y colocados, así mismo, en forma
simultánea sobre las gráficas, se localizan los puntos de intersección,
comúnmente son llamados puntos de fijación o puntos de tiro y se enumeran en
forma secuencial.

Los receptores pueden ser remolcados por el lado de estribor y los emisores, por
el lado de babor o viceversa, esto se debe a que el ruido y la turbulencia de la
prepela del barco eJimina o disminuye la intensidad de arribo de las ondas
directas; lo anterior se incrementa cuando la velocidad de crucero es mayor de 5
nudos y produce más ruido, provocando la alteración de los registros, por lo tanto,
para evitar esto, los estudios geofísicos se realizan con olas no mayores de 1.5
mts. ni vientos mayores de 30 Km/hr, además esto dificulta el control del barco e
incrementa los errores en las gráficas.

La segunda etapa de los trabajos de campo, se refieren a la fase del muestreo del
fondo marino, el cual, consiste en la obtención de muestras de los suelos y
específicamente de los sedimentos no consolidados; la selección de los sitios de
muestreo, se hace con base al registro del perfilador electromecánico, donde es
posible identificar, por diferencia de contrastes, los sedimentos no consolidados.

[Link].2. Exploración y Muestreo

El equipo que se emplea para muestrear los suelos marinos, es el muestreador de


pistón tipo

Kullenberg ( pistón sampler ), cuya penetración máxima es de 15 mts. y recupera


muestras hasta de 11.75 cm (4 5/8") de diámetro; la profundidad de muestreo
varía de acuerdo al tipo de suelo, y para elegir dicha profundidad, se toma como
base un muestreo de hasta 3 mts. para suelos arenosos y hasta 15 mts. para
Iodos y arcillas muy blandas, a partir de esto se selecciona la profundidad de
muestreo adecuado.

El equipo de muestreo mencionado anteriormente, se presenta en la figura V .1 y


en la tabla V .1, se dan algunas características de éste y otros equipos usados con
el mismo fin. En la figura V .2, se puede ver un esquema del procedimiento de
muestreo con el muestreador de pistón; durante el viaje al fondo, el muestreador
se sumerge con una velocidad máxima aproximada de 1 m/seg, cuando la pesa
del gatillo toca el fondo, el gatillo suelta al muestreador y éste cae libremente
hasta enterrarse en el suelo.
La penetración del muestreador en el suelo es función de la energía cinética de
impacto y el tipo desuelo, por lo tanto, la velocidad de caída en el tiempo t
depende de : a) La masa y densidad media del muestreador b) El coeficiente de
razonamiento entre el agua y el muestreador (aproximadamente 0.12 V2dyn/cm2)

c) La distancia h entre el punto de disparo (en el tiempo to y velocidad V o) y la


posición del muestreador en el tiempo t.

En la gráfica de la figura V.3, se presenta la variación de la velocidad de caída del


muestreador con respecto a su masa ya la altura de caída, como se puede ver, la
energía cinética de impacto (1/2 MV2) es una función no lineal de la masa, y el
amortiguamiento de dicha energía depende de la resistencia por fricción lateral
(interna y externa) entre el suelo y el barril muestreador.

El rendimiento varía de acuerdo a las condiciones meteorológicas, al tipo de suelo


y al tirante de agua y como promedio, se obtiene de 20 a 25 muestras de 3 [Link]
largo por día en tirantes de 40 a 50 mts. de profundidad.

La calidad geotécnica de las muestras, dependen del operador del equipo, de la


forma del

muestreador y de como penetra el muestreador en el suelo, por lo anterior, se ha


observado que las mejores muestras se obtienen con muestreadores de grandes
diámetros, adaptados con buenos coeficientes de forma y constantes altas de
velocidad de penetración. Los coeficientes de forma de penetración o coeficientes
Hvorlev se define como (ver figura V.4)

Ci = ~XIOO de

Co=~XIOO dt

Ca = dw2 -de2 XIOO de2 L

Ce=-ds
donde:

de Diámetro interior de la zapata de corte


dw Diámetro exterior de la zapata de corte

dt Diámetro exterior del barril muestreador

ds Diámetro interior de la camisa muestreadora

L Longitud del barril muestreador

En la tabla V .2 se dan los valores máximos de dichos coeficientes, para obtener


una calidad geotécnica aceptables en muestras de Iodos y los valores para un
muestreador saclant CEN, considerado como uno de los mejores muestreadores
de tipo kullenberg para Iodos. La alteración que sufre la muestra, por la
inmovilización del piston en el momento de tensar el cable principal, es importante
en agua someras, ya que, en el momento de impacto del muestreador con el
suelo, la tensión repentin1l del cable produce un golpe de ondas, el cual, se
propaga a través del cable y el pistón, por lo tanto, cuando el cable es corto, las
ondas rebotan varias veces entre el pistón y el malacate, originando así una
alteración mayor de la muestra.

Realmente, el empaque de las muestras es muy sencillo, debido a que la muestra


queda almacenada en el tubo de plástico que contiene el pistón y por lo tanto, al
extraer dicho tubo de muestreador, la muestra ya queda empacada y únicamente
habrá que sellar ambos extremos del tubo para evitar perdida del material y
humedad; posteriormente, las muestras se col09an verticalmente sobre un
armazón de madera y se trasladan al laboratorio. Cuando las muestras llegan al
laboratorio, inmediatamente se les determinan algunas características y
propiedades mediante las siguientes pruebas:
1.- Descripción y clasificación según SUCS

2.- Granulometria

3.- Contenido de agua

4.- Contenido de carbonatos

5.- Peso especifico relativo de sólidos

6.- Relación de vacíos


7.- Limites de Atterberg

8.- Peso especifico sumergido

9.- Prueba de veleta miniatura

10.- Prueba de penetrometro de bolsillo

Con las pruebas mencionadas anteriormente, se puede tener una idea de la


capacidad de carga superficial, capacidad de anclaje y otras operaciones
involucradas con los suelo marinos.
[Link].3. Pruebas de Barco

Las pruebas efectuadas en el barco son: las pruebas en sitio propiamente dichas y
las realizadas en el pequeño laboratorio a bordo. Las pruebas en sitio comprenden
la de la penetración que se lleva a cabo silmultaneamente con el muestreo, como
ya anteriormente se ha descrito, y la de veleta y otras sondas que serán
presentadas a continuación :

La veleta de torsión (figura V .5) operada a control remoto, consta de una barra de
acero que lleva insertado cuatro láminas delgadas de acero en posición
longitudinal que va conectada a una serie de barras huecas hasta alcanzar la
profundidad de ensaye. En el interior de las barras van alojados los cables del
sistema electrónico de medición, una vez hincada la veleta en el fondo de la
perforación se aplica en el exterior, en la parte superior de las barras, un par de
torsión que se incrementa progresivamente hasta que el suelo falla, la magnitud
de la torsión aplicada al suelo se lee mediante el puente de Wheatstone operando
en la cubierta del barco. La veleta con todas sus barras se introduce hasta el
fondo de la perforación a través de las barras huecas del sistema de cable-guía,
esta prueba es efectuada en algunos estudios de la sonda de campeche.

El cono penetrómetro estático (figuray.6) cuya primera versión se denomina cono


holandés, ha sido empleado durantes muchos años para medir la resistencia a la
penetración de los suelos que se correlacionan en forma directa con la resistencia
al corte, una de las versiones del cono penetrómetro estático cuenta con una
superficie cilíndrica de fricción inmediatamente arriba de la punta; por medio de
este instrumento se mide en forma separada y continua la resistencia de punta y
fricción lateral obteniéndose una información adecuada para el diseño de
cimentaciones piloteadas, el cono penetrómetro ha sido utilizado en forma limitada
en la sonda de campeche.

El presiómetro de Ménard (figura V. 7) fue diseñado para determinar las


propiedades de los suelos con base en la medición de sus características
esfuerzo, deformación, empleando una celda de expansión activada en el interior
de una perforación. La interpretación de las lecturas de un ensaye presiometrico
permite estimar el valor de módulo de deformabilidad y la presiones limite que el
sub suelo puede soportar, a partir de estos parametros se calculan la capacidades
de carga de los suelos ensayados. Este aparato ha tenido, en general, una
utilización muy limitada en estudios fuera de costa y en la sonda de campeche no
se ha empleado, las pruebas ejecutadas en el laboratorio a bordo son las que se
realizan empleando la veleta miniatura, el torcómetro, el penetrómetro de bolsillo y
la de compresión simple, las tres primeras se efectúan cuando la muestra aún esta
en el interior del tubo muestreador, la veleta miniatura (figura V.8) es un
instrumento pequeño análogo en forma y operación a la yeleta convencional
accionada a control remoto.

El torcónometro (figura V .9) es una veleta muy pequeña portatil que está provista
con seis u ocho láminas de acero que se hincan en la muestra del suelo. La
torsión se aplica manualmente y su valor se leen en un circulo graduado que tiene
en la cara superior.
El penetrómetro de bolsillo .-es un instrumento manual que consta de una pieza
tubular que se aloja en su interior una barra móvil a manera de pistón que,
mediante la presión ejercida manualmente y registrada
en un manómetro, se inserta en la muestra una longitud constante. El instrumento
está calibrado para proporcionar en forma directa el valor de la resistencia al corte
de los suelos cohesivos.
La prueba de compresión simple o de compresión no confinada se efectúa con
una máquina o marco de carga (figura V .10) mediante el cual puede someterse a
compresión axialmente una probeta cilíndrica, las máquinas de compresión
pueden ser del tipo prensa hidráulica o de tomillo

Pruebas Complementarias de Laboratorio .-Las pruebas que se efectúan en los


laboratorios en territorio tiene la finalidad de complementar la información
recabada en los ensayes en el sitio y en el laboratorio del barco y así obtener los
datos para la determinación de la capacidad de carga lateral de los pilotes. A partir
de estas pruebas complementarias se evalúan las principales propiedades índice
de los suelos como: el contenido de agua, los limites de consistencia, la densidad
de sólidos, el peso volumétrico, la granulometría; las propiedades mecánicas
como son: la resistencia y deformabilidad de los suelos a partir de las pruebas
triaxiales rápidas y consolidadas -no drenadas, todas las pruebas señaladas son
las convencionales de laboratorio de mecánica de suelos.

En las pruebas de compresión triaxial se aplica al espécimen una carga lateral por
medio de un líquido a presión, además de la carga aplicada en el sentido
longitudinal. Mediante las prueba triaxiales se trata de simular, con el
confinamiento lateral, cel estado de esfuerzos inicial al que esta sometido el suelo
a diferentes profundidades. La prueba triaxial rápida representa la condición de
una falla súbIta, producidas por cargas de impacto o transitorias, en la que el suelo
es un deposito reciente, no consolidado. La triaxial consolidada- no drenada
representa la misma condición de falla repentina, pero en capas de suelo que ya
se

han consolidado bajo el peso de los depósitos suprayecentes, las arenas de


partículas calcáreas y carbonatadas en el suelo marino de la Sonda de
Campeche, son susceptible a la trituración y roturas de sus ligaduras y, en
consecuencia, a la disminución de su resistencia y de su capacidad de carga, esto
propicio que se incluyera dentro de las pruebas de laboratorio denominada II
solubilidad en ácido clorhídrico II los resultados de estas pruebas son una medida
indirecta del contenido de partículas carbonatadas, la prueba consiste en formar
una suspensión con cierta cantidad del suelo yagua con ácido clorhídrico y
después de un tiempo determinado vaciar el agua; secar y pesar la muestra; esta
operación se repite varias veces observando en cada paso la pérdida de peso de
la muestra, la prueba se termina cuando la pérdida en peso es nula o
insignificante. Se considera que a mayor contenido de carbonatos, mayor será la
solubilidad en el ácido clorhídrico
La actual situación mundial de la industria petrolera ha enfocado una mayor
tención al desarrollo tecnológico que fomente el ahorro en el uso de hidrocarburos
y al estudio de nuevas fuentes de energía.

En México, este hecho es negativo, puesto que los programas de exploración


petrolera se han visto muy reducidos.
Hoy los principales yacimientos requieren de la instalación de plataformas marinas
adicionales, por lo que las diversas actividades enfocadas a la exploración,
explotación y procesamiento de los hidrocarburos deben llevarse a cabo con
programas que impliquen menores inversiones. Esas plataformas tienen que estar
diseñadas para soportar las solicitaciones a las que estarán expuestos durante su
vida útil. Por lo tanto, es importante nontar con datos del sitio de instalación con un
alto grado de certidumbre, lo cual se puede lograr realizando estudios e integrando
datos del suelo como del viento y oleaje.

La relevancia de la integración de estos estudios se sustenta en la nueva etapa de


exploración de recursos naturales que el país tiene y concretamente en la Sonda
de Campeche (figura 1.1). En este trabajo se presenta un panorama de los
estudios que se tienen que llevar a cabo con objeto de obtener información para el
análisis y diseño de estructuras marinas en aguas profundas. Cabe aclarar que en
el ámbito petrolero en México, para propósitos prácticos, el término aguas
profundas se refiere a aquellas regiones que rebasan los 200 m de tirante. Sin
embargo, estrictamente la clasificación de aguas profundas, está regida por la
relación entre el tirante y la longitud de ola (d/L).l

A continuación se señalan algunos puntos importantes relacionados con la


obtención de la información para el diseño de plataformas, antes de ser planteada
su integración.

-Los estudios eran obtenidos por separado los cuales no siempre proporcionaban
todo lo necesario para el diseño, tal como datos Geológicos,
Localización de la Sonda de Campeche en México Oceanográfico y
Meteorológico, ya que se consideraba que los aspectos más importantes para la
cimentación y diseño eran cubiertos por los estudios Geofísico y Geotécnico.
-En otras situaciones se sobreponían las zonas exploradas o faltaba información
en otras ya que no había continuidad de un estudio con otro.

En parte la causa de ese problema era por que no había una compañía que
englobara todos los estudios ya que se especializaban en una o dos áreas
especificas.

-El alto costo de la realización de estudios que cubrieran varias ramas de la


investigación, para la obtención de datos Geológicos, Geofisicos, Geotécnicos,
Oceanográficos, Meteorológicos y Sísmicos. La realización de este último, estudio
ha sido recomendado como una norma de seguridad, debido a que México es un
país con alto riesgo sísmico.

Antecedentes

El inicio de la exploración en aguas profundas (tirante de mas de 200m ) ha


requerido la ayuda de asesores con experiencia en la exploración marina en
grandes profundidades. Aunque cabe mencionar que se ha confirmado la
favorable situación que tiene México con respecto a otros países dedicados a la
explotación petrolera, ya que los mantos de hidrocarburos se encuentran a
profundidades someras,2 en México el mayor tirante en que se ha instalado una
plataforma es de aproximadamente 80 mts.

Los requerimientos de operación de un sistema para la explotación de un campo


petrolero, sí como los parámetros de diseño deben estar planeados en las bases
del usuario. La información allí presentada debe ser lo más confiable posible con
objeto de evitar cambios durante la ejecución substancial del diseño implica el
conjunto de actividades interdisciplinarias como son Geología y la conjunción
adecuada de estudios Geofísicos para el aporte de datos precisos de localización
y elección del sitio de muestreo para el estudio Geotécnico, del que a su vez se
obtendrá información de diseño para la cimentación, y de los datos
Oceanográficos y Meteorológicos, que darán información necesaria para el diseño
de la sub estructura y superestructura, se denominará estudios integrales para el
análisis y diseño de estructuras marinas. En la figura 1.2 se muestra el resultado
de un trabajo realizado con la conjunción de los estudios Geológico y Geofísico.
Objetivo de la Tesis En este trabajo se darán los lineamientos para la ejecución de
los estudios Geológicos, Geofísicos, Oceanográfícos y Meteorológicos, necesarios
para el análisis, diseño e instalación de estructuras marinas en aguas profundas.

Se tratará también la conveniencia de que estos estudios sean complemento uno


de otros.
En el lado izquierdo de la siguiente figura 1.3 se muestra un barco Geotécnico
realizando perforación y sondeo y del lado derecho un barco Geofísico obteniendo
el tirante y la configuración del lecho marino (Figura 1.3).
[Link]. Generalidades del método sísmico

[Link].1. Ondas Sísmicas

El impulso sísmico se genera ordinariamente mediante la detonación de una


carga explosiva enterrada en el suelo. Se utilizan explosivos de gran velocidad de
propagación, tales como las dinamitas gomas. La explosión puede ocurrir en un
tiempo de micro a milisegundos, dependiendo de la naturaleza, cantidad de
explosivo y material que rodea el pozo explosivo.

El impulso generado por la explosión puede considerarse como un tren de


ondas amortiguándose con el tiempo. Como incluye un espectro amplio de
frecuencias, sería la suma de todas, cada una de ellas amortiguándose con el
tiempo.

La carga explosiva se sitúa en el fondo de un sondeo perforado a poca


profundidad por medio de una sonda portátil que forma parte del equipo de campo.
El objeto de ese sondeo es situar la carga por debajo del estrato de alteración
superficial o estrato de baja velocidad weathering, con el fin de que la energía
generada se propague en mejores condiciones desde la capa consolidada situada
bajo aquél.

En algunas zonas difíciles de perforar se ha ensayado con éxito la


detonación de cargas situándolas sobre unos soportes colocados a cierta distancia
del suelo y según ciertos agrupamientos de cargas o patterns.

Otro sistema es producir las ondas sísmicas por medios mecánicos; por
ejemplo, el impacto producido por la caída de un gran peso (un bloque de acero
de unas 3 toneladas de peso) desde una altura de 2 a 3 metros. Este sistema no
es nuevo, ya que en los comienzos de la prospección sísmica se hicieron intentos
en este sentido, que no se continuaron a causa de la mala calidad de los
sismogramas. Hoy debido al avance conseguido en el instrumental geofísico y,
sobre todo, por el desarrollo que ha obtenido el registro magnético, se consiguen
hacer buenos registros sísmicos en ciertas zonas mediante este sistema. El
aparato va montado sobre un camión y es empleado con éxito por alguna de las
compañías que operan en la provincia del Sahara.

Esta conseguido ya el método de generación de ondas sísmicas en el que la


energía se produce mediante un vibrador aplicado a la superficie del suelo. Se
trata del Método Vibroseis.
Propagación de ondas sísmicas.- Las ondas sísmicas producidas por la
detonación de una carga explosiva o por cualquier otro sistema, se transmiten a
través de la tierra en forma de vibraciones. La vibración de mayor velocidad es la
de la onda longitudinal, que es la que utiliza el método sísmico.
Cuando el medio en que se propagan las ondas sísmicas no es homogéneo,
como ocurre en la tierra, se producen los fenómenos principales de refracción y
reflexión, y además, los de difracción, dispersión y scattering.

a) Difracción.- En sismología el fenómeno de difracción es análogo al de


óptica. La difracción es tanto mayor, cuanto mayor sea la longitud de la onda
considerada; por ello en la propagación de ondas sísmicas le difracción es más
frecuente que en óptica, debido a que las ondas sísmicas tiene longitudes de
ondas muchas mayores que las de las ondas luminosas y en consecuencia la
posibilidad de encontrar obstáculos de dimensiones del orden de su longitud de
onda (varias decenas de metros) es corriente. En sísmica, la difracción ocurre
principalmente en los ángulos de las fallas, y en los productos en donde los
estratos cambian de pendiente bruscamente. Las difracciones han dado lugar a
interpretaciones erróneas.

b) Dispersión.- Es la variación de velocidad de una onda con el cambio de la


frecuencia. En un medio elástico homogéneo no hay dispersión; sin embargo, en
un medio imperfectamente elástico, como la tierra, existe este fenómeno. En
sismología se ha comprobado su existencia en el caso de las ondas longitudinales

y transversales. Actualmente no hay evidencia de que exista dispersión apreciable


en prospección sísmica, excepto en las inmediaciones de la explosión, en donde
la relación entre el esfuerzo y la deformación no es lineal. La ausencia de
dispersión es prospección es muy importante, ya que permite que el impulso se
propague sin variar su longitud. En un medio que exista dispersión, las diferentes
componentes del impulso sísmico se separa y pueden aumentar la longitud del
impulso hasta tal grado que no pueda determinarse el tiempo de llegada
exactamente.

c) Scattering.- Es la formación de pequeñas ondas que propagan la energía


en todas direcciones; se produce cuando un frente de ondas choca con partículas
libres u objetos pequeños, comparados con su longitud de onda. Este fenómeno
es mayor para frecuencias altas. Parte de lo que consideramos ruido en un
sismograma puede ser debido a este fenómeno, ya que el scattering produce
energía distribuida al azar en la superficie.

La disminución de la energía sísmica con la distancia debida a la difracción,


dispersión y scattering va acompañada de pérdidas debidas a la absorción de la
energía produciendo amortiguamiento. A la distancia d la intensidad Id  Io e-d (Io

= = intensidad en el origen), donde  es el coeficiente de absorción. Este


coeficiente parece que aumenta con el cuadrado de la frecuencia; de aquí que las
frecuencias altas quedan reducidas al aumentar la distancia a la fuente y las bajas
frecuencias se conserven.

El espectro de frecuencias asociado a una explosión es muy variable,


dependiendo de las propiedades de los materiales que rodean al punto de
explosión; pero en general la energía inicial contiene una banda de frecuencias
relativamente amplia, que va desde pocos ciclos a muchos cientos de ciclos por
segundo. Cuando el impulso sísmico viaja a través de la tierra las altas
frecuencias son absorbidas más rápidamente que las bajas frecuencias, según
hemos visto anteriormente. Este es un hecho bien conocido en sismología, ya que
en los registros de los terremotos a grandes distancias solo llegan las frecuencias
muy bajas. Asimismo, cuando se trabaja por el método sísmico de refracción, en el
que se registra a distancias superficiales mucho mayores en el método de
reflexión, la frecuencia de las ondas registradas es mucho mas baja que en el
método de reflexión. Por último, si examinamos un registro sísmico de reflexión se
observa también que las frecuencias a lo largo del registro van disminuyendo con
el tiempo, es decir, con la profundidad de las capas reflejantes.

Es conocido, asimismo, que las frecuencias asociadas a la reflexión se


mantienen en una banda que generalmente va de 20 a 100 ciclos por segundo. La
ausencia de energía reflejada por debajo de los 20 ciclos se atribuye al hecho de
que las longitudes de onda correspondientes a esta clase de frecuencias son
grandes comparadas con el espesor de las capas reflejantes. Por ejemplo,
suponiendo una frecuencia de 15 c.p.s. y una velocidad media de 3000 m/seg;
=200 metros; una capa perfectamente homogénea, de un espesor de esta orden,
es muy difícil de encontrar en la realidad. Por ello la energía de baja frecuencia no
sería reflejada y solamente podría llegar a la superficie por refracción.

Por lo que a la energía reflejada se refiere, la tierra actúa como un filtro que
atenúa las componentes de las frecuencias altas, aumentando la atenuación con
la distancia y la profundidad.

Cuando la frecuencia aumenta, las longitudes de onde se hacen


comparables con las discontinuidades verticales representadas por la
estratificación, y la reflexión se produce. Si las frecuencias son muy altas la
absorción de la energía es grande, impidiendo, por tanto, su transmisión.

Respecto a la propagación de las bajas frecuencias, se admite que éstas son


sustancialmente independientes del camino recorrido.

De todo ello nace una diferencia fundamental entre los métodos de reflexión
y refracción, ya que el espectro de frecuencias reflejadas muestra la mayor
amplitud en la banda de 20 a 100 c.p.s. y el de la refractada en la de 1 a 20 c.p.s.;
por ello los geóponos y el sistema de registro han de ser diferentes en uno y otro
caso, de forma que tengan la máxima respuesta a media o baja frecuencia.
Las leyes que rigen la propagación de las ondas sísmicas son, como ya
indicamos anteriormente los principios de Huygens y de Fermat y las leyes de
reflexión y refracción.

Velocidades de propagación en las rocas.- Se incluye a continuación una


tabla con las velocidades de propagación de las ondas sísmicas en algunas rocas.
Estas velocidades se refieren a la onda longitudinal.
m/seg

Capa 300 a
meteriorizada..……………………… 900

Aluviones 350 a
modernos………………………. 1500

Arcillas……………………………… 1000 a
…….. 2000

Margas……………………………… 1800 a
…….. 3200

Areniscas…………………………… 1400 a
…….. 4500

Conglomerados…………………… 2500 a
………. 5000

Calizas……………………………… 4000 a
…….. 6000

Dolomias…………………………… 5000 a
…….. 6000

Sal…………………………………… 4500 a
…….. 6500

Yeso………………………………… 3000 a
…….. 4000

Anhidrita…………………………… 3000 a
……… 6000

Gneises……………………………… 3100 a
……. 5400
Cuarcitas…………………………… 5100 a
…….. 6100

Granitos……………………………… 4000 a
……. 6000

Gabros……………………………… 6700 a
…….. 7300

Dunitas……………………………… 7900 a
……. 8400

Diabasas…………………………… 5800 a
……… 7100

De la tabla se deduce que las rocas ígneas tienen, en general, velocidades


mayores que las rocas sedimentarias. Las rocas tienen velocidades variables
entre los límites que indicamos, dependiendo de su profundidad y constantes
elásticas.

En general, para un mismo tipo de roca, las velocidades aumenten con la


edad geológica y para una misma roca y edad, aumentan con la profundidad.

Generalidades sobre los métodos de reflexión y refracción.- En el terreno de


la geofísica aplicada a los métodos sísmicos son indirectos y consisten en producir
un pequeño sismo artificial por medio, generalmente, de una carga explosiva
enterrada en el suelo y medir los tiempos de llegada de las ondas producidas a
unos detectores o geóponos convenientemente situados en la superficie del suelo.
Conocida la ley de propagación de la velocidad en las ondas sísmicas en el
subsuelo se llega, en función de los tiempos y distancias medidos, a la posición de
los estratos en profundidad.

Las ondas producidas por la explosión se propagan en todas direcciones, y


cuando cambian las condiciones del medio, es decir, cuando la onda sísmica en
su recorrido en profundidad encuentra un medio de propagación distinto del
anterior, parte de la energía se refleja, volviendo a la superficie, y parte se refracta
siguiendo su viaje en profundidad.
Dentro de la energía refractada y una parte que experimenta la refracción
total, propagándose entonces a lo largo del contacto entre los diferentes medios y
volviendo asimismo a la superficie.

Si detectamos en superficie los tiempos de llegada de las ondas sísmicas


reflejadas en los diversos contactos del terreno, practicaremos el método sísmico
de reflexión.
Si detectamos en superficie los tiempos de llegada de las ondas refractadas
a lo largo de los diversos contactos, practicaremos el método sísmico de
refracción.

El método sísmico consiste en esencia en producir in sismo artificial y medir


los tiempos de llegada de la onda sísmica a unos receptores superficiales o
geófonos. La grafica de los tiempos de llegada de las diversas ondas registradas
desde el momento de explosión es lo que se llama sismograma.

Estudio de la grafica tiempo-distancia o dromocrona para dos capas


horizontales y velocidades constantes. Antes de entrar en el estudio de los
métodos sísmicos de reflexión y refracción, conviene examinar la curva tiempo-
distancia para el caso más sencillo, que le de dos medios horizontales y
velocidades constantes V0 y v1.

Sea P el origen de la explosión y G1 y G2… la situación de los geófonos o


detectores. Suponemos, para mayor sencillez, que disparamos y registramos en
una línea horizontal P G1. Sea h la profundidad del contacto de los dos medios de
velocidades V0 y v1.

Las ondas registradas serán las siguientes:

1) Onda directa, cuyo recorridos serán P G1, P G2…, etc.,

P G1 P G2 x
v   ... ; siendo t , t …, etc.
0 t1 t2 t 1 2

los tiempos de llegada medidos.


Su gráfica es la recta OM, cuyo coeficiente angular es el inverso de la velocidad
V0.

2) Onda reflejada.- Los recorridos serán PA1 G1, PA2 y G2…

Si el geófono Gx esta situado a la distancia x y el tiempo registrado es t,


tendremos:

2 2 2
PA G 2PA 2 x
2
 h  x  4h y como t  2 PAx resulta que
x x x 4 V

t 
2
x
2
4h
V0
3) Onda refractada.- La onda refractada que registremos será La de tiempo
mínimo, es decir, aquella que se propague a lo largo del contacto B B2, que es la
que experimenta la refracción total.

Por lo que los recorridos sucesivos serán PBB1G1, PBB2G2, etc,


verificándose la ley de Snell.

V
sen i 0

V0
sen
90  V sen i  V

1 1

El recorrido a un geófono Gx será:

PB  BBX  BXGX  2 PB  BBX , ya que PB  BXGX

h h h V1

PB   
cos i V
2
V
2
-V
2

0 1 0

2
1- V
1

V0

V1 V0

BBX  x - 2 PS  x - 2 htgi  x - 2h  x - 2h
2 2 2

v0

1- v 1
2 v 1 v 0

Con lo que el tiempo registrado será:

2 PB BB 2h V x V 2h V V 2 2
X
 x 2h V V x
1 0 1 0

t      2h    1 0

   
V V 2 2 V 2 2 2 2 VV
V V V V V V V V  V 2 2 VV V
V V
0 1 1  0 1

 1 1 0 1

0 1 0 1 1 0 1 0
1 0

2 2
2h V1  V0 x
t  
VV V
1 0 1
que es la ecuación de una recta de coeficiente angular 1/V 1, y que, además,
tangente a la hipérbola que representa la gráfica tiempo-distancia de la reflexión
en el punto T de abcisa.

2h V0

x  2h tg i   2PS PR
2 2
V V
1 0

ya que la primera onda refractada que puede llegar es aquella en la que el


recorrido BBx por el contacto es nulo, incidiendo con un ángulo i y reflejándose en
dicho contacto. De la anterior fórmula se deduce la profundidad en el punto crítico:

x V2  V1

h  2 V V
2 1

En el caso propuesto las ondas reflejadas llegarían siempre más tarde que
las ondas superficiales, ya que su grafica (O,N1) se mantiene siempre por encima
de la grafica de éstas (OM). En realidad no sucede así, como veremos a su
tiempo, ya que la velocidad de la onda directa es muy baja, debido a que su
recorrido lo hace por la capa de alteración superficial o weathering que cubre
siempre la capa consolidada, en donde las velocidades son mucho mayores; por
otra parte, los efectos de la onda directa se eliminan con filtros adecuados.

En refracción se estudian fundamentalmente las primeras llegadas; las


primeras de todas serán las debidas a la onda directa a través de la zona
consolidada, ya que la capa de baja velocidad es lenta y las distancias utilizadas
tan grandes que la onda directa por la superficie del medio de baja velocidad no se
advierte.

[Link]. Procesado de Datos Sísmicos


[Link].1. Mapas sísmicos 3D.

La información sísmica es grabada en cintas magnéticas, siguiendo un orden


determinado en dicha cinta. A esta forma de ordenar la información a medida que
se esta grabando, se le conoce como Formato de Grabación.

Los SEG-A, SEG-B, SEG-C y SEG-D se conocen como formatos básicos de


grabación de datos sísmicos, porque en la industria se utiliza una gran variedad de
formatos que, usando la estructura básica de cualquiera de los antes
mencionados, los adaptan a su conveniencia según el equipo utilizado, lográndose
de esta manera la generación de un formato SEG-A, SEG-B o SEG-C modificado,
los cuales no reciben soporte técnico de la SEG. Esto se debe al continuo avance
del diseño de equipos de grabación de datos sísmicos en el campo, los cuales
haciéndose más eficientes, requieren también de eficiencia en el formato, cuestión
que si logra con el formato SEG-D.

[Link].2. Bloque del encabezado.

Esta parte contiene información correspondiente al número de registro (FF),


calibración del instrumento y otros datos externos como área de trabajo, cliente,
etc.

El encabezado puede ser escrito de dos maneras:


-No comprimido o usando únicamente un carácter del byte.

-Comprimido o usando el byte completo.

El encabezado de un formato SEG-A se escribe no comprimido o


desempaquetado, mientras que el encabezado de los formatos SEG-B y SEG-C
se escribe comprimido o empaquetado.

[Link].3. Identificación y Estructura de un Formato.

Una vez conocidas las características generales de un formato de grabación y la


estructura básica de los formatos que existen en la industria, se puede identificar
el formato utilizado en un determinado levantamiento. La manera práctica de hacer
esto consiste en la realización de un vaciado hexadecimal (field dump) de una
porción de la cinta magnética de campo. En dicho vaciado procedemos a
identificar el bloque del encabezado y el bloque de datos.

Todo formato está constituido por dos partes bien definidas. Ellas son:

Bloque del encabezado

Bloque de datos.

[Link].4. Formato SEG-D.

El desarrollo natural de la adquisición de data sísmica, requiere cambios en las


técnicas de grabación de estos datos. Este formato incluye los mejores rasgos o
propiedades de los formatos anteriormente descritos, más características
innovadoras adicionales.

Usando este nuevo formato, los datos sísmicos pueden ser grabados
multiplexados o demultiplexados.

Los valores de los datos pueden ser grabados en uno de varios métodos, ellos
son: 8 bits ó 16 bits usando exponente cuaternario o hexadecimal, 20 bits con
exponente binario ó 32 bits con exponente hexadecimal. La diversidad de
información a ser grabada, requiere un encabezado que defina toda la data
grabada. Cada canal sísmico puede ser indudablemente identificado así como la
rata de muestreo, filtros y tipos de datos grabados.

En este formato se definen dos secciones principales:

El bloque del encabezado y


El bloque de datos de cada traza.

El bloque del encabezado contiene las siguientes partes:

Encabezado General.

Tipo de secuencia, (scan), por encabezado.

Encabezado extendido.

Encabezado externo.

Este formato viene dado por un código numérico, que define el tipo de palabra
para los datos y la designación “multiplexado” o “demultiplexado”, la cual viene
dada también por un código numérico, ubicado bajo los bytes 3 y 4 de el
encabezado general, de la siguiente manera:
[Link].5. Multiplexado.

Hay un problema asociado a los datos sísmicos de varios canales de la grabación.


Si consideramos un cierto punto a tiempo, después de que el tiro haya encendido,
todos nuestros grupos del receptor están recibiendo datos simultáneamente.

Si estos datos se registran en la cinta tan pronto como estén disponibles, los datos
sobre la cinta estarán en la orden incorrecta. Haremos que todos los canales
contemporáneamente muestreado, seguido por todos los datos en la muestra
siguiente - los datos se piden por tiempo, no como rastros sísmicos.

Este tipo de formato de grabación, conocido como grabación multiplexada.

Los formatos SEG-A, SEG-B y SEG-C fueron usados al inicio del auge de la
exploración sísmica, es decir, durante los años 70 y 80, manteniéndose en la
actualidad los paquetes de programas capaces de leer estos formatos, puesto que
aún, se realizan reprocesamientos de datos grabados en cintas en estos formatos.
A principio de los años 90 es cuando se implementa el formato SEG-D hasta los
actuales levantamientos sísmicos de adquisición.

[Link].6. Formato SEG-Y.

El formato SEG-Y es un formato estándar usado por la SEG, para contener


solamente datos demultiplexados provenientes del campo o con algún tipo de
procesamiento realizado sobre ellos.

El formato SEG-Y consta de dos encabezados principales:

Encabezado de identificación de la cinta sísmica

y Encabezado del bloque de datos.

[Link].7. Construcción de Geométrica.

Las cintas de campo contienen básicamente los registros sísmicos que se han
grabado en la ejecución de un levantamiento sísmico. El reporte del observador
debe destacar qué relación existe entre los registros secuenciales grabados en la
cinta (FILES), las estaciones de receptores fijadas en el terreno y también en que
estación se disparó; adicionalmente debe comentar sobre cualquier alteración que
haya podido surgir en la realización del levantamiento.
Reporte de observador
Distancia entre cada disparo y cada uno de sus receptores. Puede ser un valor
positivo o negativo, dependiendo de la posición de los receptores a cada lado del
disparo, (OFFSET).

Coordenadas X,Y de cada disparo, (SHT – X,Y).

Coordenadas X,Y de cada receptor, (REC – X,Y).

Punto común en profundidad o Punto medio en Profundidad, (Common Depth


Point – Common Mid Point, i.e. CDP – CMP).

Coordenadas X,Y para cada cdp en el subsuelo, (CDP – X,Y).

Distancia absoluta entre cada disparo y cada uno de sus receptores,


(AOFFSET). Número de trazas que componen un cdp, (FOLD).

[Link].8. Intervalo de Muestreo.

Es el intervalo de tiempo para tomar muestra de la señal.


Para poder reconstruir correctamente una señal sísmica el intervalo de muestreo
debe permitir tener 2 o más muestras por ciclo de la máxima frecuencia presente
en la señal.

Aliasing

Traza sísmica, muestreada aun 1 ms.

Traza sísmica, muestreada aun 2 ms.


Traza sísmica, muestreada aun 4ms.

Traza sísmica, muestreada aun 8 ms.


[Link].9. Correcciones de Estáticas.

Son correcciones en tiempo que se le aplican a los datos para compensar el


retardo que sufren los tiempos de llegada debido a variaciones de espesor y
velocidad de las capas superficiales, así como a efectos de elevación (topografía).
Es en cierto modo, simular que los datos fueron adquiridos en una superficie en la
cual no existen los efectos antes mencionados, esta superficie de procesamiento
es llamada Datum.

[Link].9.1. Estáticas de elevación y de refracción.

Muchas áreas de interés exploratorio constituyen un problema para el analista y el


intérprete, debido a las variaciones del espesor y la velocidad de las capas
someras. En la mayoría de las zonas, la superficie está cubierta por una capa
delgada de sedimentos, poco consolidados, conocida como capa meteorizada.

Las correcciones estáticas son aquellas que se aplican a los datos sísmicos, para
compensar los efectos de la variación en las cotas o elevaciones de los puntos de
disparo (fuentes) y sus receptores, el espesor de la capa meteorizada y velocidad
de viaje de la señal sísmica en dicha capa, de tal forma que se puedan ponderar
los datos sísmicos con respecto a un nivel de referencia (Datum).
La aplicación de estas correcciones tiene dos propósitos:

Ponderar los tiempos de grabación referidos a un nivel imaginario o plano de


referencia de altura constante, llamado dátum fijo, o a un nivel imaginario o plano
de referencia de altura variable, llamado “dátum flotante”. El plano de referencia
más conocido, es el nivel medio del mar, el cual se considera que tenga una altura
cero.

Restar los tiempos de retardo ocasionados por el viaje de la señal sísmica a


través de la capa de baja velocidad, o capa meteorizada, denominada
“weathering”.

[Link].9.2. Estáticas de elevación


El cálculo de la corrección estática por elevación se realiza mediante las
siguientes expresiones matemáticas:

D  Ed
Tc Disparos

Vc
D  Er
Tc Receptores

Vc

Donde Tc: el tiempo de corrección para disparos y receptores, D: el plano de

referencia, Ed: el valor de elevación de cada disparo, Er: el valor de elevación de


cada receptor y Vc : la velocidad de corrección, donde ésta se obtiene de la
pendiente de las primeras llegadas en un registro sísmico.
Fig. Apilado bruto sin estáticas de campo
Fig. Apilado bruto con estáticas de campo.
[Link].9.3. Estáticas de Refracción

Estas condiciones están gobernadas por la Ley de Snell:

Donde

i = ángulo de incidencia

r = ángulo transmitido

V0 = velocidad de la capa superior

V1 = velocidad de la capa inferior


Forma en que viaja un pulso refractado.

(ONDAS CRITICAMENTE REFRACTADAS)


Para obtener una solución estática por el método de refracción se deben seguir los
siguientes pasos:

q Selección de las primeras llegadas (Picking).

qDefinición de refractores.

qCálculo de la velocidad del refractor.

qAnálisis del tiempo de retardo (Delay Times)


Generación del modelo geológico cercano a la superficie y cálculo de estáticas.
[Link].9.4. Reflexiones Fantasmas.
Uno de la forma más común de rayo indeseable asociada a la adquisición sísmica
marina es la reflexión del fantasma. La superficie del mar es un reflector casi
perfecto. La energía que viaja del tiro (debajo del agua) está viajando con un
medio con una velocidad de cerca de 1500 metros por segundo, y una densidad
cerca de 1.025 grams/cc. El aire sobre la superficie del mar tiene una velocidad
acústica de cerca de 350 metros por segundo, con una densidad de cerca de
0.0013 grams/cc. Tapar éstos en la ecuación de “RC” da un coeficiente de
reflexión de debajo la superficie del mar de cerca de -0.9994, o la reflexión casi
perfecta (con la revocación de la muestra).

Aquí están las trayectorias de rayo que demuestran las rutas posibles para los
reflectores del fantasma:

La energía viaja del tiro, hasta la superficie, refleja casi perfectamente (pero
muestra invertida) y después viaja una trayectoria similar al rayo directo del tiro al
reflector.

El rayo que vuelve del reflector hace exactamente la misma cosa en el receptor.

Puede haber cuatro diversos rayos producidos:


Rayo directo

Ghosted solamente en el tiro

Ghosted solamente en el receptor

Ghosted en ambos tiró y receptor

Dado que las trayectorias de rayo reales tendrán más curvatura que ésos
demostrados en el diagrama, las trayectorias de rayo del tiro a la superficie y la
parte posterior abajo serán otra vez casi verticales. Si asumimos que son
verticales, podemos calcular la diferencia del tiempo (en milisegundos) entre el
rayo directo y ghosted el rayo como aproximadamente 4/3 de las épocas la
profundidad del tiro en los metros (estamos asumiendo una velocidad del agua de
1500 m/s). El mismo cálculo se puede hacer para el receptor.
Si, por ejemplo, nuestra profundidad del tiro es 6 metros, después la diferencia del
tiempo entre el directo y ghosted rayos será cerca de 8 milisegundos. Si nos
asumimos que la energía del tiro consiste en un sistema de ondas constantes del
coseno de la frecuencia agregadas juntas (vemos más adelante que esto es una
asunción válida), podemos considerar el efecto del fantasma en los componentes
individuales de la frecuencia.
Puedes ver que las frecuencias componentes de 40 y 80 hertzios son aumentadas
realmente en el fantasma, mientras que el componente de 125 hertzios (en este
caso) se quita totalmente! El rayo directo y el fantasma son exactamente 180
grados fuera de fase, y una quita la otra.

De modo general, un fantasma causará “muescas” en el espectro de la frecuencia


en las frecuencias que son un múltiplo (áspero) de 750/(profundidad del tiro) y
750/(profundidad del receptor) hertzio.
En la práctica, las profundidades del tiro y del receptor se eligen generalmente de
modo que estas muescas de la frecuencia estén más allá de la gama de las
frecuencias requeridas para los datos sísmicos.

Típicamente, las profundidades de menos de 10 metros se utilizan, haciendo la


primera muesca > 75 hertzios.

[Link].9.5. Múltiples

Dado que la superficie del mar es un reflector casi perfecto, y que el fondo del mar
puede también ser un reflector igualmente bueno, no está sorprendiendo que
éstos causan otros problemas para el geofísico de proceso.

La capa del agua es uno de los candidatos primeros a la generación de las


reflexiones múltiples (o de los múltiplos), que ocurren cuando la energía refleja
más de una vez en cualquier reflector. Debe ser tensionado que la capa del agua
no es la única fuente de múltiplos. Cualquier capa baja, con contraste
suficientemente fuerte de la velocidad, puede convertirse en una fuente de
múltiplos en infante de marina y datos de la tierra.

Aquí están las trayectorias de rayo que demuestran las rutas posibles para los
reflectores múltiples de la capa del agua:
La energía viaja del tiro, abajo al fondo del mar, entonces hasta la superficie,
reflejando casi perfectamente cada vez (solamente la muestra invertida en la
superficie del mar).

Entonces viaja una trayectoria similar al rayo directo del tiro al reflector.
El rayo que vuelve del reflector hace exactamente la misma cosa en el receptor.
Para los múltiplos, pueden diversas combinaciones se pueden generar en el tiro y
la posición del receptor, incluyendo “despedidas repetidas” de la energía dentro de
la capa que causa en segundo lugar, tercero las generaciones del etc. de
múltiplos.

La amplitud de los múltiplos muestra-será invertida (otra vez debido a la reflexión


de la mar-superficie) y aparecerá normalmente como imagen del reflector primario,
generalmente una cierta hora constante debajo de ella. Los múltiplos
subsecuentes también exhibirán esta “periodicidad”. ¡La amplitud del múltiplo
puede ser cualquier cosa hasta una imagen invertida completa del primario!
Los múltiplos generados por la capa del agua causarán “picos” en el espectro de
la frecuencia en las frecuencias que son un múltiplo de (áspero) 750/(profundidad
del agua).

Otros múltiplos pueden producir picos en casi cualquier frecuencia, generalmente


en el extremo inferior de la gama de frecuencia.

Éstos pueden (a un mayor o a poco grado) ser atenuados en el proceso


identificando su periodicidad, o por el hecho de que la mayoría de los múltiplos
exhibirán velocidades medias más bajas que primarias - más en esto más
adelante.

Las trayectorias múltiples posibles para los rayos son sin fin, aquí son algunas
posibilidades:
[Link].9.6. Migración

En la migración todos los reflectores con una pendiente serán movidos a su


posición verdadera y todos los reflectores horizontales se conservarán
horizontales. Dentro de la migración se tienen diferentes tipos de algoritmos
basados en la solución de onda, en el siguiente cuadro sinóptico se muestra como
se puede clasificar la migración:
Migración Tiempo 2D y 3D Preapilado

Postapilado

Profundidad 2D y 3D Preapilado

Postapilado

En la migración en tiempo, no se toma en cuenta la variación de la velocidad


lateral, y se utilizan diferentes algoritmos, como la migración Kirchhoff, usando
sumas a lo largo de hipérbolas. Otro algoritmo es el de diferencias finitas (DF) y
también se tiene la migración f-k de Slot y Gazdag. Se usará migración en tiempo
si el cambio lateral de velocidades no es riguroso, este tipo de migración reduce el
tiempo de máquina, ahorrando recursos de cómputo en el centro de proceso. Este
tipo de migración sirve como un primer paso para la migración en profundidad.
Uno de los problemas de la migración en tiempo es que se debe calcular las
curvas de difracción mediante una fórmula hiperbólica, por lo que se utiliza la
migración en profundidad la cual calcula las curvas de difracción por métodos más
exactos como:

 Trazado de rayos
 Ecuación eikonal

Mientras que la migración en tiempo se usa un modelo del subsuelo simple,


generalmente basado en la velocidad de la raíz cuadrática media (V RCM),
admitiendo variaciones laterales moderadas en la velocidad, la migración en
profundidad requiere un modelo de velocidad más preciso, basado en velocidades
de intervalo (VINT).

El tiempo de viaje mediante el trazado de rayos, se fundamente en las leyes de


Snell para la reflexión y refracción, siendo entonces más exacto para modelos
complicados pero más lentos que el método de la Eikonal.

La ecuación eikonal es una ecuación diferencial parcial que describe como el


frente de onda se propaga en un medio y su solución numérica s hace por medio
del método de diferencias finitas.

Dentro de los tipos de migración, la de preapilado (como su nombre lo dice, antes


de apilar) resuelve estructuras complejas y permite obtener las velocidades de
intervalo en el proceso, pero la migración preapilado requiere una buena relación
señal- ruido y es lenta, por lo tanto cara. La migración preapilado se aplica en
tiempo para procesos especiales que requieren trazas preapilamiento. Se aplica la
migración preapilado en profundidad para objetivos estructurales debajo de
recubrimientos complejos.

[Link].9.7. Migración en Tiempo


Aunque las secciones sísmicas de reflexión generalmente proporcionan una
configuración estructural aproximada de las capas de rocas sedimentarias
depositadas debajo del subsuelo, estas pueden ser engañosas y difíciles de
interpretar geológicamente, aún cuando la calidad del registro sea excelente.

Cuando la estructura es compleja, mostrando grandes variaciones de


estratificación horizontal, la posición del subsuelo de los puntos de reflexión en la
interface no cae justamente debajo de la fuente- receptor sino que es desplazado
de un lado a otro de este punto. La dirección y magnitud de este desplazamiento
de gran pendiente depende en la dirección y magnitud de la inclinación o echado
del estrato.

Como un resultado, se convierte la sección de registro de viaje de dos caminos al


tiempo de viaje de un camino. Esta conversión puede ser completada simplemente
dividiendo nuestra sección de registro apilada por dos.

La migración puede ser realizada considerando los siguientes términos: la sección


es considerada como una condición de frontera para un campo de ondas
gobernado por la ecuación de onda. Se mide la onda viajando hacia arriba en la
superficie de la Tierra, se deja correr el tiempo del reloj hacia atrás al tiempo t=0,
para de esta manera propagar las ondas de regreso en tiempo a sus posiciones
reflectoras del subsuelo. Esta propagación de regreso al tiempo cero puede ser
vista como la progresiva detonación de los receptores dentro del subsuelo hasta
que ellos se detengan en la superficie de reflexión. De esta manera la migración
podría ser descrita tanto como la continuación hacia abajo o la depropagación.

El primer ejemplo es el caso de un reflector plano horizontal. Como una fuente


receptor coincidentes, las ondas sísmicas se propagan verticalmente al reflector
del punto fuente- receptor y entonces de regreso otra vez al mismo punto. Los
puntos de tiempo apareciendo en la sección apilada representa el tiempo de viaje
de dos caminos, esto es, el tiempo hacia abajo para el reflector más el tiempo de
regreso de viaje de dos caminos en dos. Podríamos convertir el tiempo de registro
dentro de las velocidades equivalentes a profundidades meramente multiplicando
el tiempo de un camino por la velocidad constante. Debido a que las trayectorias
de viaje son verticales, la interface que se tiene en la sección apilada es la misma
que la interface geométrica.

El objetivo de la migración en tiempo es estructural y requiere que el proceso de


datos sísmicos produzca el mejor apilado, en cuanto a velocidades, se prefiere
que se tengan las Vrcm del DMO para transformarlas el nivel de referencia final.

Técnicas de migración:

-Método de superposición del semicírculo.


-Método de suma de difracciones.

-Métodos basados en la idea de que las secciones pueden ser modeladas como
una onda de campo próxima al offset- cero generada por los reflectores
explosivos.

-Métodos de diferencias finitas.


[Link].9.8. Difracciones

Las difracciones, puede ser un problema importante.

El fondo del mar aquí es muy duro, y tiene (utilizar un término técnico) “terrones”
de una cierta clase dispersada alrededor en él. Estos terrones son, con eficacia,
los solos puntos que reflejan la parte posterior a de la energía de todas las
direcciones.

No haciendo caso para los reflectores (primarios) reales de un momento en los


datos, esto es lo que parece la sección sísmica (sintética) de a del modelo
antedicho. Todos los datos demostrados aquí son, esencialmente, ruido.

Aparece al mismo tiempo, y con las mismas velocidades


medias (de hecho, una gama entera de velocidades) como
nuestros reflectores primarios, y puede ser quitado
solamente por su inmersión evidente (la cuesta del
acontecimiento en la sección).

Cualquier punto aislado, los “terrones” demostrados arriba,


o, por ejemplo, el extremo de un interfaz que trunque
contra una avería geológica, dispersará energía de todas
las direcciones detrás hacia su fuente.

Aunque las curvas de la difracción demostradas a la


derecha son un problema, las curvas asociadas a las
roturas en reflectores son importantes y serán discutidas

otra vez cuando consideramos el proceso de la migración necesario para colocar


correctamente nuestros datos sísmicos.
Todos los datos registrados también sufren de otras formas de “ruido indeseable”.
Algo de esto puede aparecer como acontecimientos continuos en nuestra sección
(ruido coherente), o como fondo a todos nuestros datos (ruido al azar). ¡Ambos
tipos de ruido se atenúan durante el proceso, pero el nivel final del ruido puede
todavía ser suficiente enmascarar reflexiones primarias en nuestra sección
sísmica - en algunas áreas es difícil encontrar cualquier señal!

El retiro de múltiplos confía generalmente en características reconocibles as uno o


ambos.

1. Los múltiplos irán generalmente sin cesar (algo como este curso),
repitiendo con el mismo intervalo del tiempo y gradualmente disminuyendo
en amplitud (pensar en “despedidas repetidas” de la energía en una capa
del agua).
2. Los múltiplos aparecerán generalmente en nuestros frunces del CDP con
una velocidad más lenta que la velocidad primaria. Tendrán valores más
grandes del moveout (más curvatura).

Observar el uso de la palabra generalmente arriba - cualquier cosa posible con la


combinación derecha de la geometría y de la geología.

Una clase de múltiplos internos que causa el realce evidente de la señal (y


modifica el espectro)
ocurre cuando tenemos una serie de camas cíclicas con la impedancia acústica
alta y baja alterna (velocidad mide el tiempo de densidad).

Los coeficientes de reflexión en cada interfaz (demostrado por las muestras


azules) aumentan realmente la cantidad evidente de energía que va abajo a capas
más profundas. La misma cosa sucede en la trayectoria de vuelta a la superficie.
Tal lecho repetidor del tipo demostrado arriba ocurre absolutamente regularmente
dentro de la naturaleza, causando un aumento evidente en energía sísmica debajo
de la sección “cíclica”.

Todas las consideraciones antedichas nos dan algunas pistas en cuanto a cómo
podemos quitar múltiplos. El período corto (y ésos generados por el lecho cíclico)
modificarán el espectro de la frecuencia de cierta manera, múltiplos más largos del
período tendrán un período regular, y un de baja velocidad (concerniente a
primarias).

[Link].9.9. Ruido coherente y aleatorio

Registro sísmico de dos reflectores donde se aprecian las curvas hiperbolas, que
definen dos diferentes correcciones dinámicas.

Registro

sísmico de

varios

reflectores

donde se

aprecian las

curvas

hipérbolas,

que definen

dos
diferentes

correcciones

dinámicas.

A:Curva
hiperbólica

sobre un
registro
sísmico de
campo,
producida

por la
reflexión de

una capa

horizontal en el subsuelo.

B, C, D y E: Curvas hiperbólicas sobre un registro sísmico de campo, producidas


por las reflexiones de capas con buzamiento o inclinación en el subsuelo.
Cono de ruido y onda de aire
.- Registro de campo al que se le ha aplicado un filtro de corte bajo, para
minimizar el efecto del ground roll.

[Link].9.10. Deconvolución

Aquí está un ejemplo de una de las cosas que el deconvolution es útil para - quitar
múltiplos del corto-período (o reverberaciones) del rastro sísmico.
Puedes recordar de ejemplos anteriores la respuesta clásica de la capa del agua a
una fuente sísmica. La energía despide en varias ocasiones entre el fondo del mar
y la superficie del mar, produciendo la clase de respuesta del tiempo demostrada
aquí.

En este caso las reverberaciones se producen de cerca de 88 metros de agua.


Esto da un rato de dos vías alrededor de 0.117s (2 x 88 el ÷ 1500), y un espectro
de
la frecuencia (demostrado aquí en el DB) con un pico repetidor de la frecuencia
cada 8.5 hertzios (1 ÷ 0.117).

Esta secuencia múltiple será sobrepuesta en cada reflexión dentro de nuestra


sección sísmica (como cualquier otro múltiplo generado por otras capas bajas).

Si podríamos cambiar la forma de onda de reverberación demostrada arriba en un


solo punto, quitaríamos el múltiplo de cada reflexión en la sección.
El espectro del punto parece esto (recordar que un punto contiene “todas las”
frecuencias) así que realmente estamos intentando invertir el espectro de la cosa
spiky demostrada arriba en este valor constante.

¿Es probable que poder diseñar con éxito un filtro de la limitado-longitud que
podría derrumbarse una forma de onda de reverberación de la longitud con
eficacia infinita en un punto?
Aquí es la forma de onda la reverberación de la entrada (en azul) con un cociente
de la amplitud “- de R” del pico al pico.

El filtro (rojo), que recordarás es tiempo invertido, apenas consiste en las dos
muestras “1” y “R” con el período exacto de la reverberación entre ellas.

El mover de un tirón a través de los botones demuestra que algo de la


circunvolución multiplica y que agrega los productos - todo a menos que el primer
punto salga como cero! Los primeros 3 resultados de la circunvolución son:

1x1=1

1xR-Rx1=0

2
1 x R - R x R = 0 etc.
¿Recordar nuestra formulación mucho-repetida de las ecuaciones usadas en el
deconvoluciòn?

Para el deconvoluciòn estadístico utilizamos el rastro sí mismo para proporcionar


la estimación de las reverberaciones y el otro señales asociadas a cada reflexión
en la traza rastro.

De la misma forma que los tipos previamente mencionados de deconvoluciòn


nosotros ahora tienen dos opciones en cuanto a cómo especificamos a parte
dominante del lado derecho de la ecuación antedicha - la salida deseada:

1. para clavar el deconvoluciòn especificamos nuestra salida deseada como


punto, e intentamos convertir el rastro entero en un solo punto. Esto no se
puede hacer exactamente por un filtro tiempo-limitado sino tiene el efecto
de igualar las frecuencias dentro del rastro (que intenta conseguir a un
espectro “plano”).

2. para el deconvoluciòn profético utilizamos una parte más última del rastro
como la salida y el intento deseados para diseñar un filtro que haga que el
rastro parece una versión tiempo-cambiada de puesto de sí mismo. Restar
este “modelo” del rastro real es exactamente equivalente al proceso
análogo discutido en la página anterior, pero con tiempo-cambia de puesto
y escalando diseñado automáticamente.
Aquí, está de nuevo el principio del deconvoluciòn profético.

Si podemos predecir con éxito que el necesarios tiempo-cambian de puesto y los


valores de escalamiento de los datos sí mismo entonces que restan el rastro
tiempo-cambiado de puesto escalado de sí mismo quitarán todas las reflexiones
múltiples.

En la práctica, el proceso entero de ambos que clavan o el deconvoluciòn profético


es controlado por apenas cuatro parámetros:

1. La longitud del filtro. Especificado a veces en milisegundos, a veces en el

“número de muestras” esto debe ser suficientemente largo incluir por lo


menos dos “despedidas” del tiempo máximo de la reverberación que
estamos intentando quitar.

2. El boquete (otra vez en ms o muestras). Un boquete se inserta en el filtro


que evita que el filtro cambie los datos cerca de cada reflector. Un boquete
de 1 muestra o de menos implica clavar el deconvoluciòn, cualquier
boquete más alto implica el deconvoluciòn profético. Los boquetes usados
normalmente extienden a partir de 2-10 muestras de blanquear menos
espectral de los datos y de la causa (y de ruido asociado). Los altos
boquetes de la opinión el 90% del período de la reverberación se pueden
utilizar para apenas quitar ciertos múltiplos.

3. La puerta del diseño (véase abajo). Esto debe incluir solamente las
regiones de la señal en nuestros datos (que omiten las primeras roturas) y
puede tener que ser especificado en función de la compensación del rastro.
Para los resultados óptimos una puerta del diseño debe ser por lo menos
diez veces la longitud total del filtro (longitud del filtro más boquete). Esto
puede limitar la eficacia del retiro múltiple del largo-período por el
deconvoluciòn.

4. Un nivel de ruido “Blanco”.


5. Aquí está una exhibición de algún apilado sísmico de los datos antes y
después.

6. Cuando estamos utilizando el deconvoluciòn antes del apilado (DBS),


diseñado en un expediente del tiro o un frunce del CDP (u otro),
necesitamos limitar la puerta del diseño a la parte de los datos que
contienen la mejor señal (y múltiplos). Utilizamos normalmente una puerta
del diseño que dependa de la compensación del rastro, el omitir “primero-
nos rompemos” y el ruido más último.
7. Después de apilado (o, por ejemplo, en común-compensar el frunce)
podemos utilizar generalmente una ventana más grande pero podemos
necesidad inmóvil de limitar la puerta del diseño al área de la mejor señal.
Si nuestro contenido de la frecuencia cambia rápidamente con tiempo (y
nosotros puede resolver los “diez tiempos” gobiernan mencionado arriba),
podemos desear utilizar el deconvoluciòn time-variant. Diseñamos
simplemente diversas ventanas de diverso excedente de los filtros y nos
aplicamos entonces de la misma forma que otros filtros time-variant (más
en esto más adelante).

Filtros
Utilizamos los filtros de frecuencia en los puntos numerosos durante la
secuencia de proceso - esa necesitamos en cualquier momento quitar algo del
fondo “ruido”. Los filtros se aplican casi siempre, sin embargo, en el final de la
secuencia ordenada-para arriba la sección y/o la migración apiladas finales.

Hemos visto muchos ejemplos de la filtración antes de que - aquí está un


ejemplo de un filtro low-pass (que pasa solamente las frecuencias bajas)
aplicado a una fotografía.

En este caso el filtro se aplica en ambas direcciones (a través de y abajo) pero,


discutiendo la filtración espacial en un cierto detalle anterior, ahora estamos
mirando “abajo” la dirección - filtrándose en el dominio de tiempo.

¿Si se asume que tenemos cierta idea de las frecuencias que deseamos
“pasar” a través a nuestra sección final, cómo diseñamos estos filtros?

Podríamos decidir simplemente (probando) sobre las frecuencias que


deseamos, y después aplicamos simplemente un filtro en el dominio de la
frecuencia que los pasos justos esos filtros. Por ejemplo puede ser que
apliquemos un filtro de 20-50 hertzio que pasa solamente las frecuencias a
partir del 20 a 50 hertzios.

El problema con esto es los atajos agudos en 20 y 50 hertzios. El filtro


resultante (en el dominio de tiempo) es áspero equivalente a la diferencia en
dos de las funciones de la sinc. que mencionamos anterior - una ondìcula sucia
bonita.

Para evitar los “borde-efectos” asociados a atajos agudos, utilizamos


generalmente los filtros que se inclinan por intervalos en el dominio de la
frecuencia.
El primer de éstos, conocido a veces como Ormsby o filtros trapezoidales,
aplicaciones cuatro frecuencias de definir las cuestas.

Todas las frecuencias debajo de F0 y sobre F3 se quitan, mientras que ésos de


F1 al F2 se pasan completamente. Las amplitudes se afilan encendido de F0 a
F1 y apagado del F2 a F3. Cuando utilizamos cuestas lineares (como se
muestra aquí) se conocen a veces como ahusamientos de Bartlett.

Aunque estas cuestas reducen grandemente cualquier “bamboleo indeseable”


en el dominio de tiempo, las esquinas agudas pueden causar problemas
(especialmente cuando el passband total es estrecho). Una solución a esto es
utilizar el “liso-hacia fuera anterior discutida los ahusamientos” de Hanning (a
veces llamado Lechuga-Ajustado) las cuestas.

Esto tiende para mejorar la respuesta del tiempo pero todavía tiene un
problema importante - las frecuencias fuera de F0 a F3 se quitan totalmente,
ellas se van por siempre! De hecho podemos alcanzar los atajos agudos
mencionados arriba (y los otros efectos) haciendo F0 igual a F1, y el F2 igual a
F3.
Puede ser que deseemos tan utilizar los filtros de Butterworth, sobras de

electrónico filtrando eso proporcionamos un diverso acercamiento al diseño (otra

vez, con cuatro

parámetros).
Para estos filtros especificamos los puntos del atajo (FA y FB), que serán
reducidos por 3dB (el 70%) por el filtro, y las cuestas bajos y altos (en dB/Octave)
en cualquier extremo. Las cuestas típicas varían a partir 6 a 72 dB/Octave.

Apilamiento de trazas

Los datos de diversos tiros y grupos del receptor (y, para datos 3D, de diversas
líneas) se recopilan juntos en los grupos de los rastros que comparten un punto
mediano común - el frunce del CMP (generalmente llamado el frunce común del
punto de la profundidad que no es verdad para los acontecimientos que
sumergen).
Repasemos todos los procesos geométricos aplicados a nuestros datos hasta
ahora.

Los datos de diversos tiros y grupos del receptor (y, para datos 3D, de diversas
líneas) se recopilan juntos en los grupos de los rastros que comparten un punto
mediano común - el frunce del CMP (generalmente llamado el frunce común del
punto de la profundidad que no es verdad para los acontecimientos que
sumergen).

Hemos aplicado correcciones estáticas (y el estática residual para los datos de la


tierra), para mover con eficacia nuestros tiros y geofonos (o flotando) a un dato
fijo.

Ganancia

Apenas como hay muchos ejemplos en nuestras vidas donde estaría útil un cierto
nivel de la igualación, en muchos casos de sísmico procesando un cierto nivel de
la igualación de las amplitudes de las trazas puede ser muy beneficioso.

Si, al igual que a menudo el caso, no estamos intentando conservar las amplitudes
relativas absolutas dentro de nuestros datos, la igualación de amplitudes a través
de los datos podemos mejorar el cociente de la señal/ruido.

La igualación hace apenas qué dice que lo hace - iguala las amplitudes de
remontar-a-remonta y traga escala de tiempo.

Se aplica normalmente antes del apilado para mejorar la respuesta y para eliminar
los grandes ruido de la amplitud, y/o después de apilado, al “balance” la sección
final. Utilizaremos una traza, como ejemplo para demostrar cómo trabaja.

Aquí tenemos 3 segundos de la traza sísmica (cualquiera antes o después del


apilado). Notarás que las altas amplitudes al inicio (la reflexión inferior del agua?),
la porción más débil en el centro, y la alta parte otra vez en el extremo.
Recuperación de amplitud
Imaginar una esfera de la energía el ampliarse del tiro sísmico en todas las
direcciones

Toda la energía contenida inicialmente en el tiro se separa hacia fuera sobre un


área más grande y más grande con el tiempo. Esto causa una de las pérdidas
posibles de energía en un expediente del campo, y se refiere generalmente como
divergencia esférica.

El área superficial de una esfera (el excedente del área que la energía se separa)
es tiempos de 4/3 pi que el radio ajustó (tuve que ir mirar el para arriba, que las
demostraciones justas cuánto tiempo hace mí estaba en escuela!).
Podemos calcular el radio mientras que la distancia el frente de onda viaja en una
época dada, que es igual al tiempo (t) multiplicado por el excedente medio de la
velocidad ese intervalo del vez (v). Esperaríamos, por lo tanto, que la pérdida de
energía debido a la divergencia esférica fuera proporcional (VT) a 2.

Otra causa de la pérdida de energía se conoce como atenuación inelástica. Esto


es simplemente el perdido energía debido a las partículas de la tierra a través de
las cuales los recorridos de la onda que no son perfectamente elásticos - algo de
la energía se absorbe y altera permanentemente la posición de las partículas.

La escala izquierda del diagrama superior está en el DB, y demostraciones de este


diagrama la función que fue aplicada a los datos. Intentamos un cierto
escalamiento logarítmico directo de DB 3, 6 y 9 por segundo, de dos funciones
basadas en tiempo y velocidad, y de dos funciones exóticas basadas en una
aproximación de las pérdidas inelásticas (la función de e0.002TV), y de una
función compleja que aplica una diversa corrección a cada rastro en el expediente
dependiendo de su compensación (distancia) del tiro.

El ejemplo con 6dB/second es razonable, al igual que la curva simple del T2. Los
ejemplos que usan “V” (particularmente el pasado) son buenos, pero, como I dicho
arriba, nosotros todavía no saben que las velocidades para esto tiraron (engañé
en este caso!). La función exponencial (e0.002TV) es obviamente mucho
demasiado fuerte.

Fuera de todos los éstos, una de las funciones más simples (y razonables) es la
curva del T2. Esto es simple aplicarse (multiplicar cada muestra por tu tiempo-
ajustado) y tiene la ventaja que podemos quitarla más adelante en la secuencia
(cuando tenemos información de la velocidad) y aplicar una aproximación mejor.
El retiro de la función requiere simplemente la división por tiempo-ajustado.

Aunque un tipo función de “dB/sec” causa menos distorsión dentro de los datos,
una curva del T2 es en la práctica generalmente suficiente balancear los datos
sobre la gama del tiempo. Si hay ruido excesivo en el fondo del expediente (un
tiempo más último), podemos necesitar parar el aplicar de la función en algún

momento (manteniéndola constante de entonces encendido). Esto es determinada


lo más mejor posible, de nuevo, probando.

Debemos esperar que la función del aumento sea bastante constante sobre un
área dada (proporcionada las condiciones de la fuente y de la grabación siguen
siendo similares).
Disparo crudo de campo (izquierda) y con corrección por divergencia esférica
aplicada (derecha).
Balanceo espectral

Asumiendo un sismograma de entrada con amplitudes que decaen en el tiempo.


El componente de baja frecuencia de la traza FL presenta una rapidez de
decaimiento menor que el conjunto de altas frecuencias presentes en la traza FM.
Igualmente, el conjunto de altas frecuencias intermedias presentes en la traza FM,
presenta también una rapidez de decaimiento menor que el conjunto de altas
frecuencias presentes en la traza FH; entonces una serie de funciones de
atenuación o de decaimiento, tales como G1, G2, y G3, pueden ser calculadas
para describir las rapidez del decaimiento o de atenuación por cada banda de
frecuencia.

Esto se realiza mediante el cómputo de múltiples filtros pasa – banda. El inverso


de estas funciones de atenuación, denominado atenuación o ganancia, es
aplicado a cada banda de frecuencia y luego se suman los resultados.

La principal ventaja entre el balanceo espectral y el proceso convencional de


deconvolución, es que el balanceo espectral logra normalizar mejor el espectro de
amplitudes, compensando la mayor atenuación de las altas frecuencias presentes,
en relación con la menor atenuación que sufren las bajas frecuencias. El objetivo
es obtener una traza de salida cuyo espectro de frecuencia asociado, sea lo más
plano posible.

Aplicación

El programa que realiza estos cálculos, trabaja en el dominio de frecuencias de la


transformada de FOURIER y su operación consiste en analizar bandas de
frecuencia muy estrechas, medir la amplitud media de cada banda de frecuencia y
finalmente, calcular una función que tienda a normalizar las amplitudes del
espectro. Una vez logrado el espectro normalizado, sólo falta retornar al dominio
del TIEMPO, para contar con la traza sísmica ya espectralmente balanceada.

El algoritmo divide el proceso en tres partes:

Descompone las trazas sísmicas, en series de filtros pasa - banda, de bandas


estrechas.
Aplica un escalar (i.e., multiplica por un número) a cada traza sísmica filtrada,
para normalizar o balancear los niveles de sus amplitudes relativas.

Suma las trazas ya normalizadas, a la traza sísmica de entrada, para formar la


traza de salida final.

La aplicación del balanceo espectral como procedimiento, depende de una


apropiada descomposición de las trazas de entrada y de la normalización que se
pueda aplicar a éstas.
Mediante la utilización de un análisis espectral, se puede determinar cual es la
mínima y máxima frecuencia presentes en los datos.
De acuerdo a ello, se diseña la serie de filtros pasa banda, que descompongan la
traza de entrada, desde la mínima frecuencia deseada hasta la máxima frecuencia
deseada.

AMPLITUD

ESPECTRO PLANO
0 DB

-0.5 DB 1 2 3 4 n n+1

10 15 20 25 30 35 40 45 50 55 60 70 75 80 85

FRECUENCIA (Hz)

Resolución vertical
Una de las preguntas más comunes hechas por la gente a mí sobre la sismología
de la reflexión está en el nivel del detalle que puedo ver. La respuesta a esta
pregunta aparentemente fácilmente contestada primero requiere que tenemos una
definición clara de la “resolución”.

La resolución sísmica es la capacidad de distinguir características separadas; la


distancia mínima entre 2 características para poder definir los dos por separado
más bien que como una. Pensamos normalmente en la resolución en el sentido
vertical, pero hay también un límite a la anchura horizontal de un objeto que poder
interpretar de datos sísmicos.

El criterio de Rayleigh la “medida sísmica” es una longitud de onda. Para que dos
interfaces reflexivos próximos que se distinguirán bien, ellos tengan que ser cerca
de 1/4 longitud de onda en el grueso (criterio de Rayleigh). Éste es también el
grueso donde los criterios de la interpretación cambian (AAPG Explorere, esquina
geofísica). Para gruesos más pequeños que 1/4 wavlength confiamos en la
amplitud para juzgar el grueso de la cama. Para los gruesos más en gran parte de
1/4 longitud de onda podemos utilizar la forma de la onda para juzgar el grueso de
la cama.

V 

Donde:  es la frecuencia y  es la longitud de onda

Resolución horizontal

La resolución horizontal se refiere a cómo los dos puntos de reflejo cercanos se


pueden situar horizontalmente, pero se reconozca como dos puntos separados
más bien que uno.

Si pensamos solamente en rayos entonces nunca tenemos cualquier problema


con la resolución del grado lateral de características porque un rayo es
infinitamente fino, tenemos frecuencias infinitas, y podemos detectar todos los
cambios. Nos ocuparemos del concepto de rayos contra ondas más adelante en el
semestre.
Sin embargo, el efecto que vemos en datos normales verdaderos de la incidencia
es explicado mejor por conceptos de la onda:

Comenzar, una reflexión no es energía de apenas un punto debajo de nosotros.


Una reflexión es la energía que despide detrás en nosotros de una región. Pues
las formas de onda son realmente no-planar, las reflexiones de una superficie se
vuelven sobre de una región y de un excedente un intervalo del tiempo. Señalar
que viene adentro en el tiempo casi igual no se puede separar en componentes
individuales temporal cortos. Así pues, predecimos que las reflexiones que se
pueden considerar como casi pero no absolutamente coincidente a tiempo
interfieren con uno a.

El área que produce la reflexión se conoce como la primera zona de Fresnel: la


zona de reflejo en la superficie inferior insonified por el primer trimestre de una
longitud de onda. Si la longitud de onda ha grande entonces el terminado de la
zona de el cual las vueltas reflejadas vienen es más grande y la resolución es más
baja.

La resolución horizontal depende de la frecuencia y de la velocidad.

Estudio de caso: El estudio de caso del sheriff en el explorador de AAPG


demuestra que las reflexiones en modems normales de la incidencia recogen
energía sobre un área finita que tamaño dependa de la primera zona de Fresnel.
La energía reflejada dentro de esta huella contribuye constructivo para construir un
meneo en el modem. Energía reflejada fuera de las cancelaciones de la huella
hacia fuera en el modem. Por esta razón el borde físico de los objetos que tienen
terminaciones laterales agudas no coincidimos directamente con la bajada de la
amplitud en los datos.

Factor Q

Las ondas sísmicas que viajan a través de la tierra experimentan la absorción,


significando la atenuación y la dispersión, debido la inelasticidad y a la
heterogeneidad del medio (Ricker, 1953; Futterman, 1962; White, 1983; Kneib y
Shapiro, 1995). Entender, estimar, y la compensación de la absorción de ondas
sísmicas son importantes en la búsqueda mejorar la resolución de imágenes
sísmicas, mejor interpretar los efectos de AVO, e invertirlos para las
características materiales.

La atenuación es medido por una cantidad llamada el factor de calidad Q. Este


factor puede cambiar en profundidad y en la dirección lateral. Si tenemos una
forma analítica para una función de atenuación, entonces podría fácilmente
compensar este efecto.
Los métodos para estimar Q de los datos sísmicos superficiales no están bien
desarrollados (Dasgupta y Clark, 1998). Una cierta investigación tiene, sin
embargo, publicado referente a la valoración de Q de datos bien sísmicos
verticales del perfil y de la cruz (Tonn 1991). Casi todos estos métodos utilizan la
amplitud de señales recibidas, pero ésta la información de la amplitud es a
menudo inexacto debido al ruido, a los efectos geométricos el separarse, de la
dispersión, y otro.
Cuadro 1: (a) Un acontecimiento en un frunce del CMP es generado por un
reflector en un medio absorbente, en donde la amplitud se ha normalizado
respecto a la amplitud máxima. (b) Espectros de la amplitud del acontecimiento de
la firma de la fuente de la entrada, rastro 1, 11 y 21 respectivamente de la derecha
a la izquierda.

(a) Un frunce sintético del CMP de dos reflexiones con se han normalizado la
absorción y los ruidos al azar del 10%, rastros basaron en su amplitud máxima. (b)
Valores de los factores de calidad de la entrada. (c) Factores de calidad
computados.
a) Un WTVS de un rastro sísmico verdadero. El pico escogido los puntos de la
frecuencia son marcados por los círculos y conectados por las líneas rectas.

(b) los factores de calidad estimados, el valor de la capa pasada son infinitos.

[Link].9.11. Análisis de Velocidades

Uno de los métodos mas comunes de comprobación espacial está examinando


contornos de la ISO-Velocidad a lo largo de la línea.
Las hachas aquí son iguales que nuestra sección sísmica (espacio y tiempo), y los
colores representan el valor interpolado de nuestras funciones de la velocidad en
los intervalos discretos a través de la sección.

En este caso, el rojo representa velocidades y púrpura bajas los altos valores - las
velocidades aquí son muy lisas y no esperábamos ver mucho de la manera de
variaciones estructurales en la sección.
Podríamos también computar velocidades del intervalo hacia fuera de las
funciones (que usan los intervalos pequeños del excedente de la ecuación del
DIX), y trazamos éstos en función de tiempo y de espacio.

De nuevo, este demostración un sistema liso de funciones, con una rotura obvia
de lento a ayunar velocidades del intervalo sobre llanura intermedia la sección.

Puesto que nuestros coeficientes de reflexión son dependientes en cambios en


velocidad del intervalo (y densidad), esperábamos ver una reflexión importante en
el medio de nuestra sección.

En algunos casos necesitamos estar enterados de exactamente cómo la


computadora interpola las velocidades entre nuestras funciones escogidas. A
menos que tengamos muy cuidados y sigamos cada acontecimiento del análisis al
análisis usando un apilado o un NTG inicial (cerca de frunce del rastro), somos
poco probables escoger exactamente los mismos acontecimientos en cada
análisis. El software de proceso tiene que tomar algunas decisiones en cuanto a
cómo interpolar los valores entre el CDPs analizado.
Este diagrama demuestra uno de los problemas comunes en datos sísmicos del
Mar del Norte. El reflector terciario bajo es muy fuerte y causa una “rotura
repentina” en nuestra función de la velocidad. Si, como se muestra aquí,
escogimos velocidades en CDP 500 y 700, y las inmersiones terciarias bajas
fuertemente entre estas localizaciones, un error de la interpolación pueden ocurrir.

La función generada interpolando nuestros intervalos constantes del tiempo del


excedente de las funciones (conocidos a veces como interpolación de IsoTime) da
la función negra en CDP 600. La función verdadera debe seguir la línea azul.

Dependiendo de la magnitud del error, la solución más simple es generalmente


agregar simplemente más funciones - para analizar CDPs adicional.

Un error grueso demostrará generalmente para arriba bastante fácilmente en una


exhibición de la ISO-Velocidad.

¡En este caso las velocidades fueron entradas de una sección de papel
previamente procesada, y han tratado al mecanógrafo de!

Antes de acabar nuestra discusión larga sobre velocidades, miremos más formas
exóticas de análisis y de QC.
Muchos otros métodos (y están siguiendo siendo) se han utilizado a la selección y
a las velocidades sísmicas de QC. Con la cosecha y el análisis cuidadosos de
apilar velocidades y de modelar, es posible establecer absolutamente un buen
campo de los datos sísmicos sí mismo de la velocidad del intervalo.
Es común examinar un grupo cercano de la traza, o el apilado inicial del CDP
antes de escoger el sistema final de velocidad funciona.

Si hacemos una cierta interpretación geológica de nuestra sección sísmica,


identificando las reflexiones principales a través del examen (o aún apenas a lo
largo de una línea) podemos procurar escoger el mismo sistema de horizontes en
cada análisis - conocido a veces como velocidades “rodilla-constantes” (hay el
mismo número de momentos cruciales o “rodillas” en cada función).

Podemos incluso hacer uso estos horizontes escogidos para hacer los análisis
detallados de las velocidades a lo largo de cada horizonte. En vez de computar la
apariencia de la velocidad traza para CDP seleccionado sobre el conjunto - medir
el tiempo de la gama, nosotros producen las exhibiciones para cada CDP sobre un
intervalo limitado del tiempo alrededor de cada horizonte.

Aunque es desperdiciador de tiempo (y costoso) en computadora y hombre medir


el tiempo esta técnica puede proporcionar la información detallada referente al
intervalo entre dos reflectores.

Tomografías
El tomografía de la reflexión determina velocidades e interfaces del intervalo
simultáneamente, permitiendo proyección de imagen sísmica en áreas de la
geología compleja. Puede potencialmente mejorar análisis de la velocidad de la
migración para superar tres problemas: la propagación de error de bajo a capas
más profundas, no-llanura de CRP recolecta en los puntos del análisis pobre de la
velocidad, campo altamente oscilatorio de la velocidad debido a la determinación
local del análisis de la velocidad de la migración.

El método propuesto tiene varios pasos:

1) La definición y el parametrizaciòn de verdad y se refieren al modelo de la


profundidad del velocidad-interfaz,

2) Irradiar remontar usando algoritmo del Shortest-Path junto con el montón


binario para el almacenaje de datos,
3) Construcción de los núcleos para las perturbaciones de la lentitud y del
interfaz,

4) Inversión usando algoritmo de LSQR.

5) El modelo se reparte en bloques básicos con velocidad constante. Los


interfaces entre los bloques son planar sabio del pedazo, describiendo
irregularidad y permitiendo la convexidad del trazo del rayo. Cada bloque
modelo tiene cuatro nodos de la esquina (en [Link]) y cada interfaz se
subdivide en segmentos usando el acoplamiento semi regular.
Trazado de rayos

Utilizamos el método del Shortest-Path propuesto por Moser (1991). El trazo del
rayo se basa en la solución del problema que remonta de tres puntos. Se requiere
que el rayo de la fuente a dispersar y de la dispersión al receptor sigue el rayo del
Shortest-Path.

El montón binario se utiliza para almacenar el tiempo de viaje para cada nodo de
la rejilla a lo largo de cada rayo. Este método es eficiente y se estima la hora de
cómputo de ser N*log (N), donde está número N de nodos.
[Link].9.12. Migración 2D/3D antes de Apilamiento en
Tiempo (Kirhhoff)
El Kirchhoff-tipo tiempo o la migración de la profundidad se ha utilizado como
workhorse por industria de petróleo desde el trabajo pionero de Hagedoorn (1954),
adentro cuál las superficies de la convexidad máxima fueron relacionadas más
adelante con la ecuación de onda escalar y llegó a ser familiar en la literatura de la
geofísica como “migración de Kirchhoff” (Schneider, 1978). En los tres pasados
décadas, sin embargo, migración de Kirchhoff en el significado indicado arriba
desarrollado de un solo operador de la proyección de imagen solamente a un
operador de la inversión. Esto ha conducido al desarrollo de un misceláneo de
otras técnicas de la proyección de imagen (Tygel y otros., 1993), donde están
capaces los operadores que resultan de manejar datos de AVO/AVA, esos
importantes en la caracterización de reservas conteniendo aceite.

[Link].9.13. Migración 2D/3D antes de Apilamiento en


profundidad (Kirhhoff)

La teoría de la migración/de la inversión implica que la remuneración del ángulo de


la iluminación se debe explicar en Kirchhoff migración de la profundidad del
prestack, aunque se descuida generalmente. Las comparaciones de los resultados
de la migración de Kirchhoff del subsalt con y sin la remuneración del ángulo de la
iluminación demuestran que puede la inclusión de esta amplitud que carga término
mejorar dramáticamente la imagen emigrada, justificando el aumento en esfuerzo
y coste de la puesta en práctica.

Migración de la profundidad del prestack de Kirchhoff de un cuerpo profundo de la


sal sin la remuneración del ángulo de la iluminación. Observar los artefactos de la
migración del domo de sal donde indicado.

Migración de la profundidad del prestack de Kirchhoff de un cuerpo profundo de la


sal con la remuneración del ángulo de la iluminación. Los artefactos de la
migración del domo de sal presentes en en el cuadro anterior se atenúan
grandemente.
[Link].9.14. Migración en Tiempo Post Apilamiento

La migración es el proceso que mueve los datos sobre nuestra sección sísmica
apilada en su posición correcta en espacio y tiempo.
Incluso después de nuestras correcciones de NMO y de DMO, las reflexiones de
acontecimientos que sumergen se trazan en nuestra sección del apilado en el
lugar incorrecto. Necesitan ser movidas “para arriba-sumergen” a lo largo de una
curva hiperbólica para ponerlas en el lugar derecho, la forma de esta hipérbola
dependiendo del campo de la velocidad.

Antes de que comencemos a mirar algunos de los métodos de migración, miremos


el aspecto de algunas estructuras simples en nuestra sección apilada.
En cada uno de las exhibiciones abajo, presionar los botones para cambiar entre

1) la manera que la estructura debe mirar (después de la migración) y 2) la


manera mira en nuestra sección apilada. Las dos cajas en la elasticidad derecha
una descripción del aspecto de la estructura antes de la migración y, en la caja
más baja (en rojo), de los errores podemos contar con si nuestro campo de la
velocidad es incorrecto. Recordar, él es la sección que ves cuando se presiona el
botón del número 2 que representa nuestra sección apilada - los intentos del
proceso de la migración para corregir esto de nuevo a su posición verdadera
(botón 1).
Un colmo local, o, en términos geológicos, un anticlinar aparece más de par en par
en la sección del apilado que debe ser.

La cantidad de “estiramiento” depende de las inmersiones en cada flanco - más


altas son las inmersiones más ancha la estructura aparecerá en nuestra sección
apilada.

Si nuestras velocidades son incorrectas entonces la estructura emigrada final


pueden ser más estrechas o más de par en par que la estructura verdadera - esto
podría conducir al incorrecto estima de cualquier reserva del aceite o del gas por
debajo de este “colmo”.
Un boquete dentro de un acontecimiento de otra manera continuo también genera
difracciones en cada rotura.

Nota, de nuevo, la revocación de polaridad y la interferencia entre los dos


difracciones en el centro del boquete.

Si nuestras velocidades son incorrectas en este caso, el boquete será mal


reflejado. En casos extremos “sonríe” de excedente emigró los datos, o de
explosiones de ruido al azar en los datos puede ocultar totalmente el boquete.
Finalmente una estructura muy extraña. Si tuviéramos un acontecimiento
perfectamente parabólico en la tierra, con su centro en nuestro CDP, veríamos
solamente un punto en nuestra sección apilada - el acontecimiento sería como un
espejo parabólico, concentrando toda la energía reflejada en un punto.

¡Nuestro apilado sintético no demuestra absolutamente que esto - los errores


causados por las inmersiones muy escarpadas en el programa de la “inverso-
migración” que utilicé soy igual a que ésos nosotros encontrarán en los programas
delanteros de la migración - son todos los (hasta cierto punto) una aproximación!

Agregar un error de la velocidad a los errores introducidos ya por el proceso de la


migración y esto será muy difícil a la imagen. Si tenemos un solo “punto” de la
energía en nuestra sección apilada que no corresponda a un acontecimiento
“verdadero” (por ejemplo un punto), ésta emigrará en una “sonrisa parabólica” en
nuestra sección emigrada.

Las ecuaciones que gobiernan la cantidad de corrección requerida para la


migración son relativamente simples. La cambio espacial (delta-x) y el nuevo
tiempo (emigrado) pueden ser calculados escogiendo inmersiones en una sección
apilada y usando las velocidades que apilan. Las ecuaciones son:

Donde está Tiempo en segundos, sumergir en segundos por el metro y la


velocidad en metros por segundo. Se utiliza nuestra convención normal de la
inmersión (el tiempo de aumento con el aumento de la distancia es positivo).
Debes poder ver del antedicho que cualquier acontecimiento se moverá “para
arriba-sumerge”. La ecuación del tiempo es idéntica a ésa usada para NMO, salvo
que utilizamos mitad de la velocidad (o dos veces la distancia) en el cálculo.

Un modelo sísmico apilado es un poco como la vista del mundo considerado por
alguien los cristales gruesos que usan con las lentes agrietadas! Cada cambio de
la velocidad en las capas bajas causa una curva en los rayos que introducen (y
que se van) la tierra, torciendo todo debajo de esa capa.

Los métodos usados para corregir para esto, y para emigrar los datos a su
posición correcta, confían en algo del trabajo hecho por los pioneros tempranos en
la investigación óptica que explican efectos similares sobre rayos ligeros.
Si se asume que un acontecimiento en tiempo de dos vías de 2 segundos con una
velocidad de 2000 m/s en nuestra sección del apilado, he trazado estas
ecuaciones en la derecha para todas las inmersiones posibles.

¡Observar que, incluso solamente el permitir sumerge hasta 45 grados, el


acontecimiento se podría colocar mal realmente por hasta 2000 metros y 2
segundos!

Si combinamos las ecuaciones antedichas, quitando el término de la inmersión,


podemos trazar todas las posiciones posibles para el acontecimiento emigrado.

Esta curva nos demuestra que que podría este solo “punto” en 2 segundos
(dependiendo de su inmersión) emigrar a cualquier posición respecto a esta curva

- la forma de la curva que es controlada solamente por la geometría y las


velocidades.
La migración y las variaciones de Kirchhoff en las otras técnicas de la ecuación de
onda se pueden modificar para considerar la flexión de los raypaths (por la
refracción) dentro de una sección compleja. Estas técnicas modificadas son
generalmente migraciones referidas de la profundidad - así como proyección de
imagen los datos en la posición correcta también convierten escala de tiempo
vertical de nuestra sección a la profundidad. Por consistencia convertimos a
menudo la sección final de la profundidad nuevamente dentro del tiempo desde
entonces que es la escala generalmente para una sección sísmica.

Las migraciones del tiempo y de la profundidad son solamente aproximaciones


cuando están aplicadas a los [Link] datos sísmicos. A menos que la 2.a línea
funcione exactamente a través de las inmersiones máximas en la estructura
(normalmente llamada una línea de la “inmersión”), algunos de los datos se deben
emigrar del plano de la sección en tres dimensiones. La única solución a esto es,
por supuesto, el proceso de 3D (y la migración 3D) que discutiremos abajo.
Dejar nos comenzar mirando algunos ejemplos verdaderos de la migración, y
algunos de los problemas comunes a todos los programas de la migración.

Esta sección demuestra los dos segundos inferiores de una línea sísmica que fue
procesada a cuatro segundos con la migración de la ecuación de onda.

Los botones demuestran la época (en segundos) de la “continuación hacia abajo”


computada por la aproximación a la ecuación de onda sobre ese punto. 0 es no
migrado (el apilado) y 4 (representando 4 segundos) se emigra completamente.

Nota cómo la incompatibilidad compleja en el fondo de esta sección llega a ser


gradualmente reflejada mientras que el frente de onda baja con eficacia la sección.
Todos los algoritmos de la migración sufren cuando el muestreo espacial de los
datos no es bastante correctamente a la imagen ninguna inmersiones escarpada
en la sección.

Aquí está un pedazo pequeño de la línea emigró en su intervalo del CDP del metro
del ½ de la original 12, y entonces (simplemente procesando cada otro rastro) en
un intervalo de 25 metros. Puedes ver claramente la migración “ruido” introducido
por el debajo-muestreo (o el aliasing espacial).

Antes de apilar la migración de la profundidad de una estructura compleja (y muy


escarpada) puede producir una migración muy exacta.

El campo inicial de la velocidad para estos datos fue modificado (con algoritmos
genéticos) durante el proceso para proporcionar la “mejor” migración.

Los colores en la exhibición final representan el campo computado final de la


velocidad.

[Link]. Interpretación de datos sísmicos

[Link].1. Atributos sísmicos

Son obtenidos de los datos Sísmicos, por medición directa.

El estudio y la interpretación de los atributos sísmicos nos proveen información


cualitativa de la geometría y los parámetros físicos del subsuelo. Es sabido que el
contenido de amplitud de los datos sísmicos es el factor principal para la
determinación de los datos sísmicos, como la impedancia acústica, los
coeficientes de reflexión, las velocidades, la absorción, etc. Los componentes de
la fase son el factor principal en la determinación de la forma de los reflectores, así
como de las configuraciones geométricas, etc. Los atributos sísmicos pueden ser
de velocidad, de amplitud, de frecuencia, y el índice de cambio de cualquiera de
estos con respecto al tiempo o al espacio tiene como principal objetivo proveer
información detallada al intérprete en parámetros estructurales, estratigráficos y
litológicos.

La clasificación de atributos
Los atributos pueden ser generados sobre sísmica pre apilada o pos apilada, y
antes o después de la migración en tiempo. Los atributos pueden ser clasificados
de muchas formas diferentes. Aquí damos una clasificación basada en las
características del dominio de los atributos:

[Link].1.1. Atributos Pre-Strack

Los datos de entrada son de CDP o trazas de gathers de imagen. Estos tendrán
información relacionada al azimut y al offset. Estos procesos generan grandes
cantidades de datos; sin embargo no son prácticos para estudios iníciales sin
embargo contienen grandes cantidades de información que puede ser
directamente relacionada al contenido de fluidos y a la orientación de las fracturas.
AVO, velocidades y variación de azimut son algunos de los incluidos en esta
clasificación

[Link].1.2. Atributos Pos-Apilados

El Apilado es un promedio con el cual se elimina la información relacionada al


offset y al Azimut. Los datos de entrada podrían ser apilados CDP o migrados.
Uno notaria que los datos migrados en tiempo mantendrán su relación en tiempo
así como su frecuencia y las dimensiones físicas. Para las secciones migradas en
profundidad, la frecuencia es remplazada por el numero de onda, la cual e una
función de la propagación de la velocidad y la frecuencia. Los atributos pos-
apilado son más manejables ya que se obtienen cantidades de datos en
investigaciones de reconocimiento inicial. Para estudios detallados, los atributos
pre-apilados deben ser incorporados.

Clasificación de los atributos por sus características computacionales:

[Link].1.3. Atributos instantáneos

Estos atributos son procesados muestra por muestra y representan variaciones


instantáneas de varios parámetros. Los valores instantáneos de atributos tales
como la envolvente de la traza, sus derivadas, la frecuencia y la fase deben ser
determinados a partir de trazas complejas.

[Link].1.4. Atributos de onda

Los atributos instantáneos son calculados en el pico de la envolvente de la traza y


tiene una relación directa con la transformada de Fourier en la vecindad de la
ondícula de la envolvente del pico. Por ejemplo la frecuencia instantánea en el
pico de la envolvente es igual a la frecuencia media del espectro de amplitud de la
ondícula. La fase instantánea corresponde a la intersección de la fase. Estos
atributos son llamados atributos respuesta. Todos estos atributos representan la
parametrización de Fourier de la ondícula.
[Link].1.5. Atributos físicos

Los atributos físicos responden a propiedades físicas ejemplos.:

1. La magnitud de la traza de la envolvente es proporcional a los constrantes


de impedancias acústicas

2. La frecuencia se relaciona con las capas delgadas, difracción de la ondícula


y absorción
3. Los atributos Instantáneos y la velocidad promedio está directamente
relacionado con las propiedades de la roca.

Consecuentemente, estos atributos son sobre todo usados para la clasificación de


litología y la caracterización de yacimientos Ejemplos: envolvente, frecuencia, Q,
absorción, etc.

[Link].1.6. Atributos Geométricos

Los atributos geométricos son una manifestación de como los atributos varían en
espacio y tiempo

Las medidas de continuidad lateral por semblanza es un buen indicador de


similaridad de capas así como las discontinuidades, echados y curvaturas que
contienen información de yacimientos.

Los atributos geométricos son una excelente herramienta para la interpretación


estratigráfica definen características y sus relaciones espaciales que están
directamente relacionada con los patrones de depósito y litología.

Ejemplo. Echados, bordes, terminaciones estructurales, etc.

[Link].1.7. Atributos de reflectividad y transmisión

La mayoría de los atributos, instantáneo o de ondicula; para estudiarlos asumen


las características sísmicas reflejadas de la ondicula. Se consideran las interfaces
entre dos capas. Sin embargo, la velocidad y la absorción se miden como
características que ocurren entre dos interfaces, o dentro de una capa. Por lo
tanto, podemos dividir las cualidades en dos categorías básicas basadas en su
origen. Reflectividad o Transmisión.
Fig. Datos de entrada originales
Fig. Filtro pasa bandas Butterworth

Fig. Parte Real de la traza compleja


[Link] Compleja
Fig. Envolvente

Clase de atributos

• Parte Real de la Traza Compleja


• Parte Imaginaria de la Traza Compleja

• Traza Envolvente (Energía Total)

• Derivada en el tiempo de la envolvente (Índice de cambio de energía en el


evento)

• Segunda derivada en el tiempo (aceleración) de la envolvente.

• Fase Instantánea

Usos de la Traza envolvente

• Representa contrastes de I.A. por lo tanto reflectividad

• Puntos brillantes (Bright spot)

• Posible acumulación de gas

• Limites de secuencias

• Incompatibilidades de la capas delgadas


• Inconformidades

• Cambios importantes de litología

• Cambios importantes en ambientes de deposito

• Cambios laterales indicando fallamiento

• Correlación espacial para porosidad y otras variaciones litológicas.

Segunda derivada de la envolvente

• La segunda derivada del envolvente da una medida de agudeza del pico ,


cuál puede ser más útil como una exhibición principal de la cualidad. En
exhibiciones negras y blancas, demuestra todos los picos del sobre cuál
corresponde a cada interfaz de reflejo dentro de la anchura de banda
sísmica. Puesto que demuestra todos los interfaces, el segundo derivado
del sobre tiene una anchura de banda mucho más amplia, una resolución
por lo tanto mejor.

• Demuestra todos los interfaces de reflejo visibles dentro de anchurade


banda sísmica

Cont. De la segunda derivada de la envolvente

• Demuestra la agudeza de acontecimientos

• Indica cambios agudos del litología

• Cambios grandes del ambiente deposicional, incluso cuando la amplitud


correspondiente de la envolvente puede ser baja.
• Presentación muy buena de la imagen de la superficie inferior dentro de la
anchura de banda sísmica.

• Todos los picos positivos se pueden utilizar como interfaces de reflejo


iniciales en un esquema escaso de la inversión del punto.
Fig. Segunda derivada de la envolvente
Cont. Atributos

• Frecuencia instantánea (que corresponde a la frecuencia mala del espectro


de la amplitud)

• Ancho de envolvente (alisado)

• el indicador capas finas de la frecuencia instantánea (puntos aislados de la


freq. Inst.)

• 2da. Derivada del tiempo (aceleración) de la fase.

• frecuencia dominante instantánea (que corresponde a la frecuencia mala


del espectro de energía

Frecuencia Instantánea
• La frecuencia instantánea se relaciona con la propagación de la onda y el
ambiente deposicional, por lo tanto son cualidades físicas y pueden ser
utilizados como discriminadores eficaces:

• Corresponde a la frecuencia media (del centro) del espectro de energía de


la ondícula sísmica.

• Correlacionador de carácter sísmico en la dirección lateral

• Indica los bordes de las capas finas de la impedancia baja

• Indicador del hidrocarburo por anomalía de la frecuencia baja. Esto afecta


es algunas veces acentuadas cerca de arenas sin consolidar debido al
contenido del aceite en los poros.

Cont.. Frecuencia Instantánea

• Indicador de la zona de la fractura. Las fracturas pueden aparecer como


zonas de una frecuencia más baja

• Indicador caótico de la zona de la reflexión debido a la dispersión excesiva

• Indicador del grueso de la capa Frecuencias más altas indican interfaces


agudos o lecho fino de la pizarra; frecuencias más bajas indican lecho
arena-propenso.

• Relación arenas/arcilla indicador de ambientes clásticos.


Fig. Frecuencia Instantánea
Fig. Ancho de la envolvente frecuencia promedio
Indicador de capas delgadas

• La una sola pieza de la información que podemos extraer es las


localizaciones adonde las frecuencias instantáneas saltan o viajan en la
dirección contraria. Estas inflexiones negativas son indicativas de ondículas
reflejadas. Por lo tanto, la derivada en el tiempo de la función de la fase
contendrá los indicadores para las capas finas, en la forma de variaciones
grandes de la frecuencia instantánea.

• El indicador fino de la capa es, por lo tanto, computado como la diferencia


entre las frecuencias instantáneas y el tiempo-hechas un promedio: Esta
cualidad demuestra las zonas de interferencia en fase.

• Es una cualidad física puesto que se relaciona con los acontecimientos de


espaciados cercanos.

Apliación de capas delgadas

• Puede ser utilizado en estudios detallados computado de puntos grandes


de la frecuencia instantánea

• indica el traslapado de los acontecimientos indica capas finas, cuando


lateralmente es continuo

• indica zona de no-reflejo, cuando aparece lateralmente al azar, produciendo


un.
Fig. Indicador de capas delgadas
Fig. Frecuencia dominante instantánea

Cont. Atributos

• Ancho de banda instantáneo


• Factor Q

• Amplitud Normalizada (coseno de la fase instantánea)

• Envolvente (fase instantánea modulada)

• Impedancia Acústica relativa (integración de la parte real de la traza


sísmica)
Fig. Ancho de banda en frecuencia instantánea
Fig. Estimador de Q instantáneos
Fig. Amplitud normalizada (solo picos)
Fig. Amplitud Normalizada
[Link] Acústica Relativa
Atributos Geométricos

• Continuidad de eventos

• Indicador de arcilla

• Reflexiones caóticas.

• Indicador de capas paralelas

• Zonas de discordancias

• Echados de semblanza máxima

• Variación de echados

• Semblanza
[Link] de eventos.
Fig. Semblanza máxima de echados
Fig. Semblanza lateral máxima de echados
[Link] de echados
Fig. Variación de Semblanza
Fig. Indicador de capas paralelas
Fig. Indicador de capas caóticas
[Link] de zonas caóticas
[Link] de arcillas (azules indica arcilla)
Fig. Indicador Discontinuidades
Atributos Ondícula (wavelet)

• Envolvente de traza (energía total)

• Primera Derivada en tiempo de la envolvente (cambios de energía en


eventos)

• Segunda Deriva en el tiempo de la envolvente (aceleración)

• Fase de ondícula

• Frecuencia de ondícula
[Link]ícula envolvente
Fig. 2da. Derivada en el tiempo de la ondícula
Fig. Fase de ondícula
Fig. Frecuencia de ondícula
Fig. Frecuencia promedio de ondícula
Atributos de ondícula cont…

• Aceleración de fase de la ondícula

• Frecuencia dominante de la ondícula (RMS)

• Ancho de banda de la ondícula

• Factor Q de la ondícula relativa

• Polaridad aparente de la ondícula.


Fig. Frecuencia (RMS) de la ondícula dominante
Fig. Ancho de banda de la ondícula
Fig. Q relativa de ondícula
Fig. Polaridad aparente de la ondícula
Papel de la prospección sísmica en el conjunto de un programa de
exploración.

Los métodos sísmicos, sobre todo el método sísmico de reflexión, son


métodos bastante exactos, pero son muy caros. Como por otra parte son métodos
de detalle, su empleo debe ser cuidadosamente estudiado y planeado.

El planteamiento de una prospección sísmica es inicialmente, como el de


toda prospección geofísica, un planteamiento geológico. Asimismo es muy
importante también, desde el punto de vista económico, considerar el objetivo que
se persigue y el posible valor económico de la información que pueda
determinarse.

Antes de iniciar una prospección sísmica convendrá, pues, conocer con el


mayor detalle posible la geología de la zona estudiar. Las posibilidades de
utilización de los métodos sísmicos varían según los problemas concretos que se
quieran investigar. En este sentido es oportuno tratar separadamente las
prospecciones petrolíferas, las prospecciones mineras y los trabajos de ingeniería
civil, las tres grandes aplicaciones generales de los métodos geofísicos de
prospección.
a) Papel de los métodos sísmicos en la exploración petrolífera.- En este el
campo mas amplio de utilización de los métodos sísmicos y al que deben el
extraordinario desarrollo a que han llegado en la actualidad.

Una vez escogida una cuenca sedimentaria o una parte de ésta, en donde
se supone o se conoce la existencia de formaciones porosas aptas para
almacenar hidrocarburos, se debe seguir como norma general, el siguiente
escalonamiento en la prospección petrolífera de la misma:

1) Reconocimiento geológico general para estudiar las posibilidades


petrolíferas de la cuenca en lo que se refiere a la presencia de posible
formaciones de rocas madre y de rocas almacén, formaciones impermeables que
puedan cubrir estas últimas, espesores de los sedimentos, etc. y estudios
tectónicos regionales que establezcan la posibilidad de formación de estructuras
anticlinales, fallas, etc., que puedan servir de trampas para la acumulación de
hidrocarburos.

2) Reconocimiento geológico y fotogeológico de detalle sobre la parte del área


que los estudios anteriores hayan determinado ser de mayor interés. Estos
estudios incluirán la preparación de un mapa geológico de detalle del área elegida,
así como los consiguientes estudios estratigráficos (para determinar la serie
estratigráfica de la zona y su posible variación de espesores y facies a través de la
misma), tectónicos, paleontológicos y paleogeográficos.

3) Prospección geofísica de la zona por algún método geofísico de


reconocimiento, previo estudio general de las propiedades físicas (densidad y
susceptibilidad magnética principalmente) de las rocas de los afloramientos mas
importantes de la zona. Los métodos de más utilidad son el aeromagnético y el
gravimétrico. Incluso en algunos casos concretos el telúrico y el eléctrico de
resistividades. Cabe aclarar que en prospección petrolífera los métodos citados
suelen emplearse como métodos de reconocimiento y para determinar la
profundidad del basamento.

El objeto de intercalar alguno de estos métodos de reconocimiento (que son


relativamente baratos) antes de los métodos sísmicos es doble; por una parte,
avanzar en el conocimiento estructural de la zona y por otra, reducir la zona a
investigar por los métodos sísmicos.
4) Prospección sísmica, Con los datos de los estudios anteriores se eligen las
áreas que presentan mayor interés; hecho esto, se puede proceder a su estudio
por los métodos sísmicos.
La prospección sísmica se lleva a cabo de una manera sistemática,
trazando primeramente unas líneas de reconocimiento general que luego se
van completando con las líneas de detalle correspondientes.

Como norma general se suele emplear el método sísmico de reflexión,


ya que es mucho más preciso, cómodo y fácil de interpretar que el sísmico de
refracción. No obstante, en zonas particulares (en donde, por ejemplo, el
sísmico de reflexión no proporciona buenas reflexiones, o en zonas dudosas)
se utiliza también el sísmico de refracción. A este respecto modernamente se
viene empleando cada vez más el método sísmico de refracción, si bien el
grueso de la prospección sísmica se realiza por reflexión.

5) El paso final de la prospección petrolífera, antes de iniciar la


campaña de sondeos exploratorios, es la interpretación de todos los datos
obtenidos anteriormente de manera que encajen lo mejor posible uno con
otros. Por ello la interpretación o interpretaciones finales debe realizarlas un
geofísico geológico.

b) Papel de los métodos sísmicos en la prospección minera.- El empleo


de los métodos sísmicos en la prospección minera está normalmente limitado
a casos particulares, ya que las condiciones del yacimiento generalmente
imposibilitan su localización por dichos métodos. La explicación es sencilla,
ya que la utilización con éxito de los métodos sísmicos supones 2
condiciones fundamentales:

1) Regiones de tectónica tranquila, con superficies de contacto planas o


casi

planas.

2) Existencia de contrastes de velocidad entre los sucesivos medio del


subsuelo.

De estas dos condiciones, la primera no suele cumplir casi nunca


prospección minera.
Por ello los métodos sísmicos serán útiles generalmente en la
prospección de yacimientos minerales sedimentarios y en zonas poco
plegadas. En este sentido pueden ser muy útiles en la prospección de aguas
subterráneas, en la localización de potasa y a veces en la prospección de
carbones. Indirectamente son utilizados algunas veces para determinar el
espesor de la capa meteorizada.

También podría gustarte