Clase N°2
Clase N°2
Clase N°2
Cuando consideramos un frente de ondas plano que se mueve en un medio 1 y choca con
una interfase, incide con un angulo θ 1. Según la ley de snell, la relación que se cumple es la
siguiente:
sin(θ1 ) sin (θ2)
=
v1 v2
Sabemos que
podemos tener las siguientes alternativas:
{
θ1=sin
V1
( )
V2
θ1 >θc donde θc =90° El rayo tiene reflexión total
−1
no hay energía por la interfaz∋transmitida
implicaθ c =90 ° se genera elrayo refractado críticamente
headwave
Hagamos análisis especial a la relación que obtenemos cuando hay un rayo refractado
criticamente:
V 1 −1
θ1=sin( ) → θc =90 °
V2
Si V2 < V1 el rayo refractado criticamente no existe, solo hay reflexiones y
transmisiones convencionales. Es decir, que la única manera de que esto se pueda
cumplir es que la velocidad de las capas aumente con la profundidad.
Si θ1 >θ cla ley de snell no puede cumplirse y por tanto no habrá rayo refractado
críticamente, este es el caso de reflexión total. Sin embargo, tanto el refractado y el
reflejado si existirán.
De toda esta gama de posibles rayos, el que nos interesa en la sísmica de refracción es el
rayo refractado criticamente, este rayo como viaja por la interfase entre los dos medios y
adquiere la velocidad del segundo, es el rayo mas rápido y llega primero. Por ende, este
rayo es el que corresponde a los primeros arribos.
Ejercicios 2
(a)
Con el ppio de Huygens podemos ver que cada punto del segmento BC que correspondería
al frente de onda del rayo criticamente refractado, se comportan como fuentes para ondas
secundarias. Estos puntos por ejemplo el P y el Q emiten ondas con frente de onda circular.
Si para un punto P cuando el rayo ya ha llegado a C se comporta como punto fuente, tendrá
un radio de curvatura V1(t-tp) siendo t el tiempo que le tomo a la onda en llegar a C, y en Q
V1(t-tq) y así sucesivamente para todo punto de ese rayo. En C obviamente el radio es nulo.
Tomando el triangulo verde dibujado tendremos la siguiente relación:
V 1(t−tp)
sin(θ c ¿ )= ¿
PC
Y a su vez:
PC=V 2(t−tp)
Entonces:
V 1( t−tp)
PC=V 2 ( t−tp )=
sin(θ c ¿ )¿
Así se tendrá que:
V1
V 2=
V1 V 1 −1
sin(θ c ¿ )→sin(θ c )= → θc =sin( ) ¿
V2 V2
b) Esta onda existe por las condiciones del medio, V2>V1 y θ1 <θ c, pero además este rayo va
a subir a la superficie por el principio de Huygens, en el que todo punto del frente de onda
de este rayo, es una nueva fuente de ondas secundarias. Este rayo ocurre teniendo también
reflexión (común) y transmisión.
Caso de capas paralelas – Una interfaz
Ejercicio 3
(a) Determinar el tiempo de viaje de la onda refractada criticamente y el rayo directo.
Lo que hay que hacer es definir los segmentos verde-rojo-amarillo en tiempos y esos
tiempos sumarlos, es decir: tiempo_total = tiempo_verde + tiempo_rojo + tiempo_amarillo
Normalmente este procedimiento es computarizado ya que a medida que tenemos mas
refracciones y mas reflexiones, será mas difícil hacer el picking.
De todas maneras para casos ideales como este, cuando lo hacemos manualmente,
aplicamos la metodología de mínimos cuadrados a cada recta de refracción. Es decir,
tomamos los primeros arribos correspondientes a la primera llegada (x, t), aplicamos
mínimos cuadrados y luego graficamos, luego buscamos los otros arribos y así
sucesivamente.
Ahora si queremos hallar una relación entre xcross y el espesor de la capa simplemente de
xcross
la expresión anterior despejamos y vemos que nos da:
Z
Explicación importante:
Si V2>V1, es evidente que V2 - V1 < V1 + V2, en tal sentido si hacemos una aproximación de
la siguiente forma:
Esta es una consideración muy importante porque resulta la distancia mínima a la que
podemos colocar un geófono de la fuente, o por así decirlo, es el offset mínimo.
¿Qué valor debe tener V2 para que xcross=z?
Al
−5
hacer las cuentas, obtenemos que V 2= V 1 esto es una incongruencia, no tiene sentido
3
porque las velocidades siempre deben ser positivas. Así que es imposible decir que xcross
será igual al espesor de la capa. Es decir, no existe valor de V2 para que esto suceda, o
simplemente no existe rayo directo en este caso.
Para el modelo de N capas, el tiempo de llegada se puede escribir como una formula
generalizada de tiempos de intercepto + tiempos de llegadas directas.
Para el modelo de tres capas: