Caracterización de Arenas y Gravas con Ondas Elásticas

Doctorado en Ingeniería – Facultad de Ingeniería – Universidad Nacional de Cuyo – Año 2007
Lic. Armando Luis Imhof
Se reemplazó el valor de y de la ecuación C.5 (a) en la (b), obteniéndose una ecuación polinómica
en x de orden 2:
A0x
2
donde :
+ A1x
+ A2
= 0
A0
= A + B ⋅ m + C ⋅ m
A1
= B ⋅ bb + 2 ⋅ C ⋅ m ⋅ bb + D + E ⋅ m
A2
= C ⋅bb
2
2
+ E ⋅ bb + F
(C.7)
El comando matlab roots([A0 A1 A2]) calcula las raíces (x1,x2); luego reemplazando en la
ecuación C.5 (a) se obtendrá el par (y1,y2).
Finalmente aplicando el conocido Teorema de Pitágoras se calculan las distancias d en el interior de
la anomalía para todos los rayos y el tiempo perturbado debido a los mismos. Finalmente se
representa gráficamente.
y
yc
Y
X
ϕ
x
xc
x1 x2
(a) (b)
Figura C.1. Esquema algoritmo de simulación de rayos sísmicos. (a)
Elipse rotada y trasladada. (b) Intersección con el rayo ij-ésimo.
y
d
x
299