UNIVERSIDAD NACIONAL DE INGENIERÍA

FACULTAD DE INGENIERÍA EOLECTRICA Y ELECTRÓNICA

TESIS:

Análisis comparativo de un arreglo de antenas de GPR de doble

polarización para un modelo de subsuelo de multicapas de la
zona arqueológica de Caral

PERFIL DE TESIS PARA OPTAR EL TÍTULO DE INGENIERO DE

TELECOMUNICACIONES

AUTOR:
Tanith Alcira Guerra Huaranga

ASESOR:
Daniel Diaz Ataucuri

Lima, 30 de agosto de 2020

AGRADECIMIENTO

A mi madre Luisa, mi padre Luis, mi hermana Yelitsa y a mi hermano Cheeyo, que

en cada momento me apoyan en seguir cumpliendo mis metas. También
agradezco a las personas con las que conviví a lo largo de mi crecimiento
universitario y profesional, aportándome grandes experiencias enriquecedoras.

¡Muchas gracias!
 1 INDICE DE CONTENIDO
Contenido
INDICE DE CONTENIDO.....................................................................................................................iv
RESUMEN.........................................................................................................................................vi
INTRODUCCIÓN.................................................................................................................................7
CAPITULO I........................................................................................................................................8
1. PLANTEAMIENTO DE LA INVESTIGACIÓN..............................................................................8
1.1 ANTECEDENTES BIBLIOGRÁFICOS..................................................................................8
1.2 DESCRIPCIÓN DE LA REALIDAD PROBLEMÁTICA............................................................9
1.3 FORMULACIÓN DEL PROBLEMA..................................................................................10
1.4 JUSTIFICACIÓN E IMPORTANCIA DE LA INVESTIGACIÓN..............................................11
1.5 OBJETIVOS...................................................................................................................11
1.6 HIPÓTESIS....................................................................................................................12
1.7 VARIABLES...................................................................................................................12
1.8 INDICADORES..............................................................................................................13
1.9 METODOLOGÍA............................................................................................................13
CAPITULO II.....................................................................................................................................14
2. MARCO TEÓRICO Y CONCEPTUAL........................................................................................14
2.1 FUNDAMENTOS DEL GEO RADAR................................................................................14
2.2 ANTENAS DE GPR.........................................................................................................15
2.3 GUÍA DE ONDA CUADRADA.........................................................................................19
2.4 PROPAGACIÓN MULTICAPAS.......................................................................................29
2.5 SOFTWARE DE SIMULACIÓN........................................................................................30
CAPITULO III....................................................................................................................................32
3. DISEÑO Y SIMULACIÓN........................................................................................................32
CAPITULO IV....................................................................................................................................53
4 SISTEMA GEORADAR-SUBSUELO.........................................................................................53
4.1 ANÁLISIS PREVIO..........................................................................................................54
4.2 ANÁLISIS DE PÉRDIDAS DE RETORNO..........................................................................57
4.3 PATRÓN DE RADICACIÓN Y GANANCIA DEL ARREGLO.................................................58
4.4 COMPORTAMIENTO DEL CAMPO ELECTRICO A LO LARGO DE LA CAVIDAD................59
CAPITLO V........................................................................................................................................61
5 ANÁLISIS COMPARATIVO.....................................................................................................61
CONCLUSIONES Y RECOMENDACIONES..........................................................................................65
BIBLIOGRAFÍA..................................................................................................................................66
 2 RESUMEN

La presente tesis trata del estudio, diseño y simulación de tres tipos antenas para
Radares de Penetración Terrestre-GPR en la banda VHF / UHF, optimizando sus
parámetros con la finalidad de obtener un ancho de banda mayor o igual a los
300MHz, ancho de haz- HPBW de 60° a 70° y ganancia de 5dB a 7dB. Estos
objetivos son debido a que la señal recibida por las antenas GPR se distorsiona
por la absorción de energía electromagnética en las multicapas o múltiples capas,
la refracción y la reflexión causada por cambios en la permitividad en las múltiples
capas del subsuelo. Además, el ancho de banda reducido implica una baja
resolución en las imágenes de radar. Luego, el sistema propuesto utiliza cuatro
antenas, dos para transmisión y dos para recepción, todas con doble polarización.
Sin embargo, este sistema obtiene un acoplamiento no deseado entre las
antenas, lo que dificulta la detección en cada polarización.

En el capítulo III, se realiza un análisis comparativo entre los tres tipos de

antenas, arreglo de antena para GPR de dobles polarización y del sistema
Arreglo-Multicapas(Subsuelo), e identificar la antena más óptima para GPR de
doble polarización ortogonal utilizando los simuladores ANSYS-HFSS y COMSOL,
obteniendo un rango de penetración mayor a 12 metros con una buena resolución
a este alcance.

Esta tesis forma parte del proyecto de investigación N°140-2017-FONDECYT

titulado “Estudio de perfiles en el subsuelo mediante métodos de inversión y
modelado electromagnético utilizando un Georadar de VHF aplicado a la zona
arqueológica de Caral” es financiando por FONDECYT del CONCYTEC.
 3 INTRODUCCIÓN

La exploración del subsuelo con Radares de Penetración Terrestre está ganando

espacios muy importantes en arqueología internacional y nacional, debido que se
trata de un modo no invasivo de buscar restos arqueológicos, reduce
significativamente los costos de la investigación arqueológica y porque es mucho
más práctico y rápido de encontrar los restos arqueológicos.

Debido a la necesidad e importancia de descubrir las costumbres de la Cultura

Patrimonial como es la zona arqueológica de Caral, se busca debelar estos
misterios encontrando restos de gran importancia para la interpretación de
arqueólogos. Debido a esta necesidad se trabaja en nuevos dispositivos de GPR
con una penetración mayor de 12 metros, con una resolución apropiada para
diferenciar entierros, para ello se utilizarán antenas que estén centradas en
frecuencia de VHF/UHF y con modelos resonancia en doble polarización
mejorando la resolución y obtener más información para la interpretación de la
prospección de imágenes.
Por ello, se plantea diseñar tres tipos de sistemas de antena que cumplan o se
aproximen a los objetivos mencionados, para finalmente realizar un análisis
comparativo entre los arreglos de antenas de doble polarización que se propagan
a lo largo de múltiples capas o multicapas.

Todo el sistema es diseñado para este fin, y puesto a prueba en terrenos de Caral
para la verificación de su funcionamiento.
 4 CAPITULO I000000
1. PLANTEAMIENTO DE LA INVESTIGACIÓN
1.1 ANTECEDENTES BIBLIOGRÁFICOS

Actualmente las aplicaciones de telecomunicaciones se encuentran en diversas

áreas multidisciplinares, como medicina, educación, arqueología, militar,
aeroespacial, entre otros. En el campo de la arqueología en nuestro país aún
existen muchas zonas arqueológicas que no han sido exploradas debidos a las
dificultades que se presenta en la excavación, ya que al momento de realizar
puede dañar estructuras, restos u objetos arqueológicos, otro inconveniente que
presenta es el tiempo que demora en ubicar los restos arqueológicos, debido a
que no se conoce la ubicación exacta del objeto. Últimamente las
telecomunicaciones y la electrónica vienen desarrollando dispositivos que facilitan
las exploraciones no invasivas en las áreas arqueológicas, como los Radares de
Penetración Terrestre-GPR. En el Perú los GPR han sido aplicados en estudio de
los suelos en zonas de gran importancia histórica como Machu Picchu como se
describe en [1], también sirvió para estudios en los asentamientos de la cultura
Pachacamac y lograr encontrar restos Según se detalla en (Bray, 2015) [2]

Hoy en día, los GPR trabajan en la banda UHF con un alcance máxima de 10
metros en el subsuelo. En muchas exploraciones realizadas en zonas
arqueológicas se requieren que el alcance sea mayor a 12 metros, debido al
cambio climático en el área y a la antigüedad de la zona arqueológica. Recurrimos
a la doble polarización en el arreglo de antenas de GPR para obtener datos con
buena resolución a profundidades mayores de 12 metros o 15 metros, debido a
que la doble polarización en el arreglo de antenas permite extraer una mayor
identificación de forma y dimensiones de un sector barrido por el radar [3]
En la última década se viene desarrollando técnicas o métodos que puedan
predecir y analizar el comportamiento de la señal que emite la antena del GPR en
el subsuelo, lo cual está constituida por múltiples capas. Esto es debido a que la
propagación de las ondas electromagnéticas se ve afectada por las características
del medio en que la señal se propaga, teniendo en cuenta los fenómenos que se
presenta al momento de la propagación en múltiples trayectorias como la reflexión,
refracción y difracción de la onda propagada [4]. Autores como Hu y Chew (2001),
investigaron un algoritmo de onda plana no homogéneo para medios de múltiples
capas [5]; Mosig y Melcon (2002) presentaron una técnica eficiente para evaluar
las funciones de Green asociadas con los medios en capas [6], entre otros. La
información que proporcionan estos y otros autores permitirá interpretar el
radagrama obtenido en la toma de datos en medios no homogéneos ayudando a
reconocer el efecto en la obtención de datos ante la presencia de posibles
interferencias de naturaleza interna [7].

1.2 DESCRIPCIÓN DE LA REALIDAD PROBLEMÁTICA

En la fabricación del radar de penetración terrestre se tiene en cuenta el tamaño

de la antena del Georadar para ser transportada a lo largo del terreno a analizar,
el rango de frecuencia de operación de la antena y el diseño de la antena. En el
diseño de la antena podemos mencionar de ejemplo la antena de bocina, antena
Vivaldi, antena Bowtie, antena espiral, etc., donde el diseño depende del
requerimiento de la aplicación como en la detección de minas antipersonales, el
estudio de estructuras patrimoniales, la detección de restos arqueológicos,
auscultación de caminos, estudios geomorfológicos, entre otros [8]. Los
parámetros que definen la antena son estudiados y modificados a lo largo de los
años, pues estos parámetros ayudarán a obtener un buen alcance, una buena
resolución del objeto y eliminación de señales no deseadas en la antena
receptora. Estas señales no deseadas o distorsiones en la señal se deben a
fenómenos que ocurren a lo largo de la propagación de la señal en el subsuelo,
fenómenos como la absorción de la energía electromagnética [9] y la reflexión
debido a la dependencia de la permitividad en el medio que se propaga

Una posible solución en la eliminación de la señal no deseada que recibe el

Georadar es la polarización cruzada conocido como x-pol, pues la polarización de
una antena nunca es el 100% en un solo modo. Cuando realizamos la polarización
cruzada entre los canales H y V, el ruido no está correlacionado con la señal por lo
que en la correlación cruzada H & V el valor de expectativa de ruido se instalará
en cero [9] reduciendo gran cantidad de señal no deseada. También es
recomendable obtener una ratio de cross polarización mayor a los 20dB para que
las mediciones de reflectividad obtenidas en el Georadar sean precisas. Por otro
lado, es necesario optimizar el ancho de banda para incrementar la resolución en
los radagramas. Así mismo, una optimización del ancho de haz HPBW permite
también tener un mejor manejo del primer rebote y una ganancia que permita
tener una mayor penetración.

1.3 FORMULACIÓN DEL PROBLEMA

La señal que reciben las antenas de GPR se distorsiona, debido a que la señal
propagada a lo largo de las multíplices capas en el subsuelo está expuesta a
fenómenos de absorción de la energía electromagnética y al fenómeno de
refracción y reflexión producidos por los cambios de permitividad de los medios en
el que se propaga la onda. Adicionalmente un ancho de banda reducido implica
una baja resolución en las imágenes del radar, lo cual requiere trabajar con
métodos que aumente la resolución a un mayor alcance. Esto método consiste en
el diseño del GPR que consta de dos antenas en emisión y dos antenas en
recepción, este diseño genera acoplamiento entre las antenas. Por tanto, se
plantea realizar un análisis comparativo de tres tipos de arreglos de antenas para
GPR de doble polarización optimizando el ancho de haz-HPBW para la obtención
de un alcance de penetración mayor a 12 metros, una buena resolución
optimizando el ancho de banda y obtener una buena polarización cruzada
eliminando la mayor cantidad de ruido para una menor distorsión de la señal.

1.4 JUSTIFICACIÓN E IMPORTANCIA DE LA INVESTIGACIÓN

La presente investigación busca contribuir al desarrollo de los GPR realizando

un análisis comparativo de las características de radiación de diferencias
arreglos de antenas para GPR de doble polarización propagándose en
múltiples capas y así tener un registro de las optimizaciones del ancho de
banda, ancho de haz- HPBW y polarización cruzada de la antena. Con la
finalidad de eliminar considerablemente el ruido en la señal recibida, obtener
un alcance de penetración mayor a 12 metros y mejorar la resolución del objeto
detectado. La señal emitida por las antenas será analizada mediante
simulaciones, utilizando los softwares Ansys,HFSS y COMSOL. Cabe
mencionar que la presente tesis contribuirá en trabajo futuros para mejorar la
precisión en los resultados en el pre y pos procesamiento con menor inversión
de tiempo.

1.5 OBJETIVOS

1.5.1 OBJETIVOS GENERALES

Analizar y comparar un arreglo de antenas de GPR de doble polarización para ser

utilizado en el subsuelo de multicapas de la zona arqueológica de Caral.

1.5.2 OBJETIVOS ESPECÍFICOS

a.- Identificar, analizar y comparar los tipos de antenas de GPR de doble
polarización mediante los simuladores ANSYS-HFSS y COMSOL en un
medio de múltiples capas.

b.- Optimizar e implementar las antenas para obtener un mejor ancho de banda,
de haz-HPBW y una buena polarización cruzada.

1.6 HIPÓTESIS

1.6.1 HIPÓTESIS GENERALES

Analizando y comparando un arreglo de antenas de GPR de doble polarización se

podrá utilizar en el subsuelo de multicapas de la zona arqueológica de Caral.

1.6.2 HIPÓTESIS ESPECÍFICA

a.- Se podrá Identificar, analizar y comparar los tipos de antenas para GPR de
doble polarización mediante los simuladores ANSYS-HFSS y COMSOL en
un medio de múltiples capas.

b.- Se podrá optimizar e implementar los tipos de antenas para obtener un mejor
ancho de banda, de haz-HPBW y una buena polarización cruzada.

1.7 VARIABLES

1.7.1 VARIABLES INDEPENDIENTES

 La doble polarización o polarización cruzada.

 Ancho de banda.
 Ancho de haz.
 1.7.2 VARIABLES DEPENDIENTES

 Profundidad de penetración (Rango)

 Resolución de información

1.8 INDICADORES

 El ancho de banda de la frecuencia de operación de la antena

 Dimensión y forma aproximada de los objetos enterrados
 Co-polarización

1.9 METODOLOGÍA

El trabajo de tesis comprende de 3 etapas, donde se tiene una primera etapa del estudio
del estado del arte que describe el comportamiento de los Georadares, tipos de arreglo de
antenas para radares de GPR y el comportamiento del onda que propaga el Georadar a lo
largo de una cavidad, la segunda etapa describe el diseño, simulación y optimización de
los tres tipos de antenas y en la tercera etapa se discuten los resultados obtenidos
realizando un análisis comparativo de las tres antenas , dando finalmente las
conclusiones y recomendaciones para futuros trabajos. En resumen, las etapas son:
 Estudio del estado del arte
 Diseño y Simulación
 Análisis y presentación de resultados
 5 CAPITULO II000000
2. MARCO TEÓRICO Y CONCEPTUAL

2.1 FUNDAMENTOS DEL GEO RADAR

Los Georadares o GPR, son sistemas compuestos por arreglo de antenas

transmisoras y receptoras de ondas electromagnéticas, sistema de control y
almacenando datos para el procesamiento e identificación de resultados. El
aspecto a tener en cuenta para el trabajo de un Georadar es ampliado a
continuación:

2.1.1 RANGO O PROFUNDIDAD DE PENETRACIÓN

Uno de los factores más importantes en los Georadares es la limitación que

presenta en el alcance de profundidad, debido a las dispersiones que producen en
las múltiples capas que presenta el subsuelo a la onda propagada, el rango
máximo en un medio de ε r=1 (espacio libre) es de la forma [10]:

0.25
Pt G A t σ
Rmáx =
[ ( 4 π )2 S min ] (1)

Donde:
Potencia de transmisión es Pt, Ganancia de las antenas transmisoras y receptoras
es G, apertura efectiva de las antenas es A t, tamaño del target dependiendo de la
polarización de las antenas está definida por (σ ) y la
Potencia mínima captada por la antena receptora es Smin.

Para medios que presentan coeficiente de permitividad diferente a 1, la ecuación

es de la forma:
 0.25
Pt G A t σ
Rmáx =
[ ( 4 π )2 S min ] B (2)

Donde, B es una constante de corrección dependiendo del medio al que va ser

sometido, es decir depende del coeficiente de pemitividad de cada medio.

2.1.2 RESOLUCIÓN DE PROFUNDIDAD

Los Georadares presentan tres tipos de resolución, una de las cuales es la

resolución en profundidad, el cual indica cada que distancia en penetración
discrimina objetos en resolución, esta resolución está dado por [10]:

C0
∆r= (3)
2 Bw

Donde:
C 0, es la velocidad en el espacio libre
Bw, es el ancho de banda del pulso o señal a transmitir

2.2 ANTENAS DE GPR

En las últimas décadas el diseño de las antenas para radares GPR aplicados a la
arqueología, minería, agricultura, entre otros. Algunos tipos de antenas adecuadas
para los sistemas GPR es la antena de bocina, la antena Vivaldi, la antena Bowtie,
la antena espiral, la antena de bocina TEM, la antena de ranura cónica, la antena
dipolo cargada, la antena monopolo, entre otros más [11]. Estas antenas son en
su mayoría usadas para radares GPR por su alta efectividad, alta ganancia y un
buen ancho de banda.

Una parte de los Georadares está compuesta por antenas, esta es esencial y
dependiente del tipo de modulación a usar, debido a que las antenas irradian el
recorrido en su ancho de banda señaladas por el modulador y acondicionados por
el circuito de RF.

2.2.1 ANCHO DE BANDA DE LA ANTENA

El ancho de banda de la antena se define como el rango de frecuencias que

trabaja o está irradiando la antena con un funcionamiento óptimo. Uno de los
criterios más comunes para definir el ancho de banda es tomando como referencia
el parámetro de dispersión S11 de la antena menores a -10dB, es decir cuando la
antena irradia con un coeficiente de reflexión menor a 0.1 este parámetro es muy
utilizado para ver el ancho de banda, tal como se muestra en la ilustración 1. En
este caso tomando valores de S11 menores a -10dB se tendrá un ancho de banda
de 119MHz.

Ilustración 1 Parámetro S11

Cabe resaltar que el ancho de banda tiene una relación directa con la resolución
que nos puede brindar el equipo, mientras sea mayor este ancho de banda la
resolución también será mejor [10]
 2.2.2 GANANCIA DE LA ANTENA

La ganancia es el parámetro que indica la máxima densidad de potencia radiada

en cierta dirección (U (θ ,ϕ)) entre la densidad de potencia radiada de una antena
isotrópica Pin/ 4) a distancias y potencias equivalentes [12]. Para fines del
Georadar esta dirección de potencia máxima es la dirección por donde se quiere
propagar la onda. La ecuación de la ganancia está dada por:
U (θ ,φ)
G=4 π
Pin
Donde, G es la ganancia de la antena transmisora y Gr es la ganancia de la
antena receptora. La ganancia depende mucho del diseño y tamaño de la antena.

2.2.3 PATRÓN DE RADIACIÓN DE LA ANTENA

El ancho de haz-PHBW de radiación de una antena se calcula el ángulo

conformado por los puntos ubicados a -3dB del lóbulo principal del patrón de
radiación normalizado como se muestra en la ilustración 2, este haz de radiación
influye en la resolución espacial que se obtienen de las ecuaciones de GPR. Los
parámetros más importantes a mencionar son los siguientes:

 Dirección de apuntamiento: Es la dirección de máxima radiación, es la

misma que del lóbulo principal.
 Lóbulo principal: Es el margen angular en torno a la dirección de máxima
radiación. Está comprendido entre dos mínimos relativos. 
 Lóbulos secundarios: Son lóbulos menores de máximos relativos, de
valores menores al lóbulo principal. 
 Ancho de haz-PHBW: Es el margen angular de direcciones en las que el
diagrama de radiación de un haz toma el valor de la mitad del máximo. 
 Relación Font to Back (F/B): Es el cociente en dB entre el valor de máxima
radiación y el de la misma dirección y sentido opuesto.
 Ilustración 2 Patrón de Radiación [13]

La polarización de las ondas electromagnéticas radiadas influyente en la antena, ya que

una posición circular es muy usada en radares de una sola polarización como en radares
mono estáticos, pero al tratar con radares de varias antenas y polarizaciones es requerido
que las antenas posean diferente polarización y es logrado por el mismo diseño de antena
o por arreglo espacial de las antenas, donde estas polarizaciones comúnmente se desean
que sean ortogonales. La diferencias entre las polarizaciones se da debido a la rotación
de la dirección del campo eléctrico a través del tiempo como también el cambio de la
magnitud resultante, generado polarizaciones lineales, circular y elíptica, como se muestra
en la ilustración 3, donde un parámetro de medida para saber a cuál de ellos pertenecen
es el radio axial (AR), donde para la polarización lineal se obtiene que el AR tiende a ∞,
para la polarización circular se tiene que AR=1, y para la polarización elíptica se tiene el
rango de AR < 1, ∞ >.
 Ilustración 3 Tipos de polarización de ondas electromagnéticas [14]

2.2.3.1 CROSS POLARIZATIO

La polarización cruzada, es la polarización ortogonal a la polarización que se está

discutiendo. Por ejemplo, si se pretende que los campos de una antena estén
polarizados horizontalmente, la polarización cruzada en este caso es polarización
vertical. Si la polarización es polarizada circularmente a la derecha (RHCP), la
polarización cruzada está polarizada circularmente a la izquierda (LHCP).

Este término surge porque una antena nunca está 100% polarizada en un solo
modo (lineal, circular, etc.). Por lo tanto, a veces se presentan dos patrones de
radiación de una antena, el patrón de radiación co-pol (o componente de
polarización deseado) y el patrón de radiación de polarización cruzada.

La polarización cruzada se puede especificar para una antena como un nivel de

potencia en dB negativo, lo que indica cuántos decibelios por debajo del nivel de
potencia de polarización deseado es el nivel de potencia x-pol.

2.3 GUÍA DE ONDA CUADRADA

Podemos definir una guía de onda cuadrada como un tubo conductor a través
del cual se transmite la energía en forma de onda electromagnética que se
 propaga por medio de reflexiones en sus paredes internas, debido a que estas
paredes son perfectamente conductoras (PEC). La intensidad de los campos
es máxima en el centro a lo largo de la dimensión X, y debe disminuir a cero al
llegar a las paredes [15], tal como se muestra en la ilustración 4.

Ilustración 4 Guía de onda cuadrada [15].

A continuación, se explica electromagnéticamente el comportamiento de la onda a

lo largo de la guía de onda.
 2.3.1 ECUACIÓN DE HELMHOTZ
∇x⃗
E=−iωμ ⃗
H (4)
∇x⃗
H =iωε ⃗
E (5)
∇ .⃗
E =0 (6)
∇ .⃗
H =0 (7)
Sabemos que:

∇x∇x⃗ A ) −∇ 2 ⃗
A=∇ ( ∇ . ⃗ A (8)

Aplicando las ecuaciones (8), (6) y (7) a las ecuaciones (4) y (5), se tendrá:
∇x∇x⃗
E=−iωμ ∇ x ⃗
H =ω 2 με ⃗
E

ω 2 με ⃗ E ) −∇2 ⃗
E=∇ ( ∇ . ⃗ E

∇2 ⃗
E −ω 2 με ⃗
E=0 (9)
H , obteniendo así:
De igual manera para ⃗

∇2 ⃗
H + ω2 με ⃗
H =0 (10)

2.3.2 GUÍA DE ONDA

Para una guía de onda, se tendrá los siguientes parámetros:

E(x , y, z )=[ e⃗ ( x , y ) +e ( x , y , z ) e^x ] . e−iβZ

⃗ (11)

H (x , y , z) =[ ⃗h( x , y ) + hz( x , y, z ) e^z ] .e−iβZ

⃗ (12)

E(x , y, z ), corresponde al campo eléctrico propagándose al largo de la guía de onda.

⃗
Por otro lado, ⃗ H (x , y , z) , corresponde al campo magnético propagado en el eje Z, a al
largo de la guía de onda.
Recordando las ecuaciones de Maxwell, con los parámetros ⃗J =0 y ρ=0, se
tendrá:
∇x⃗
E=−iωμ ⃗
H (13)
∇x⃗
H =iωε ⃗
E (14)
El rotacional del campo eléctrico puede representarse:
 e^x e^y e^z
E= ∂
∇x⃗
∂x
Ex [ ∂
∂y
Ey
∂
∂z
Ez ]
∂ Ez ∂ E y ∂ Ez ∂ Ex ∂ E y ∂ Ex
∇x⃗
E= ( ∂y
−
∂z
e^x −
∂x) (
−
∂z
e^y +(
∂x
−
∂y z )
) e^ (15)

Comparando las ecuaciones (13) y (15)

H=
−iωμ ⃗ ( ∂∂Ey − ∂∂Ez ) e^ −( ∂∂Ex − ∂∂Ez ) e^ +( ∂∂Ex − ∂∂Ey ) e^
z y
x
z x
y
y x
z

∂E ∂ E ∂E ∂ E ∂E ∂E
−iωμ( H e^ + H e^ + H e^ )=( ) e^ −( )
z y z x y x
x x y y − z z − e^ +( x− ) e^ y z
∂y ∂z ∂x ∂z ∂x ∂y

Igualando los términos de los vectores e^x , e^y y e^z , se tendrá:

∂ Ez ∂ E y
−iωμ H x = −
∂y ∂z
∂ Ez
−iωμ H x = +iβ E y (16.a)
∂y
Análogamente para H y y H z, obtenemos:

∂ Ez
iωμ H y = +iβ E x (16.b)
∂y
∂ E y ∂ Ex
−iωμ H z = − (16.c)
∂x ∂y
De igual manera realizamos para el rotacional del campo magnético.

e^x e^y e^z

H= ∂
∇x⃗
∂x
Hx [ ∂
∂y
Hy
∂
∂z
Hz ]
∂Hz ∂H y ∂ Hz ∂ H x ∂ Hy ∂Hx
∇x⃗
H= ( ∂y
−
∂z
e^x −
∂x
−
∂z) (e^y +(
∂x
− ) e^
∂y z ) (17)

Comparando las ecuaciones (14) y (17)

 ∂ Hz ∂ H y ∂Hz ∂Hx ∂ H y ∂ Hx
iωε ( E x e^x + E y e^y + E z e^z)= ( ∂y
−
∂z ) (
e^x −
∂x
−
∂z )
e^y +(
∂x
−
∂y z
) e^

Igualando los términos de los vectores e^x , e^y y e^z , se tendrá:

∂ Hz ∂ H y
iωε E x = −
∂y ∂z
∂ Hz
iωε E x = +iβ H y (18.a)
∂y
Análogamente para E y y E z, obtenemos:
−∂ H z
iωε E y = −iβ H x (18.b)
∂x
−∂ H y ∂ H x
iωε E z = − (18.c)
∂x ∂y
Para obtener E x, E y, H x y H y , utilizamos las ecuaciones (8) y (12). Para despejar
H x, utilizamos las ecuaciones (16.a) y (18.b).

∂ Ez 1 ∂ Hz
−iωμ H x =
∂y
+iβ [
iωε ∂ x(+iβ H x ] )
2 ∂ Ez ∂Hz 2
−i ω με H x =ωε −β −i β H x
∂y ∂x

2 2 ∂ Ez ∂Hz
i H x (β −ω με)=ωε −β
∂y ∂x

2 ∂ Ez ∂Hz
i K c H x =ωε −β
∂y ∂x

i ∂ Ez ∂Hz
H x= 2
.(ωε −β ) (19.a)
Kc ∂y ∂x

Donde:

K 2c =K 2 −β2 , es el número de onda de corte.

β, es la constante de propagación.
K, es el número de onda.

Análogamente para H y , E xy E y:
 −i ∂ Ez ∂Hz
H y= 2
.(ωε +β ) (19.b)
K c ∂x ∂y

−i ∂ Ez ∂ Hz
E x= 2
(β + μω ) (19.c)
Kc ∂x ∂y

i ∂ Ez ∂Hz
E y= 2
(−β + μω ) (19.d)
Kc ∂y ∂x

En las guías de ondas existen tres maneras de propagarse, estos son: modo TEM,
Modo TE y modo TM. Para nuestro caso solo hablaremos del Modo TE.

2.3.2.1 MODO TE
EL campo eléctrico se propaga a lo largo de Z, por el cual E z=0. Esto es aplicado
a las ecuaciones (19.a), (19.b), (19.c) y (19.d), obteniendo lo siguiente:

−i ∂ H z
H x= β (20.a)
K 2c ∂ x

−i ∂ H z
H y= β (20.b)
K 2c ∂ y

−i ∂ Hz
E x= 2
μω (20.c)
Kc ∂y

−i ∂Hz
E y= 2
μω (20.d)
Kc ∂x

Para éste modo K c ≠0 , por lo tanto K 2c =K 2 −β2 .

H z( x, y , z). ∇ 2 ⃗
De la ecuación (10), calcularemos ⃗ H z +ω 2 με ⃗
H z =0

∂2 ∂2 ∂2
( 2
+
∂ x ∂ y ∂z 2
+ 2
+K 2 ⃗)
H z( x , y, z )=0 (21)

∂
Considerando ⃗ H z( x, y , z)=¿ ⃗
h z(x , y) e−iβz, aquí la derivada parcial será H z( x, y , z) se
=0 .⃗
∂z
puede reducir a una ecuación de onda bidimensional para ⃗ h z(x , y)

∂2 ∂2
( 2
∂x ∂y
+ 2 )
+ K c2 ⃗
h z( x, y)=0 (22)
Demostración:
2
∂2 ⃗
hz(x , y) e−iβz ∂ ⃗h z(x , y)e ∂2⃗
h z(x , y)e
( ∂x 2
+
∂y 2
−iβz

+
∂z 2
−iβz

+ K2⃗
h z(x , y) e−iβz =0 )
2
∂2 ⃗
hz(x , y) e−iβz ∂ ⃗h z(x , y)e
( ∂x 2
+
∂y 2
−iβz

+(−iβ )2 ⃗
hz(x , y) e−iβz + K 2 ⃗
hz(x , y) e−iβz =0 )
Factorizando ⃗
−iβz
h z(x , y) e a todos , quedaría:

∂2 ∂2
( 2
∂x ∂y
+ 2
−β 2 + K 2 ⃗ )
h z(x , y) e−iβz =0

Como K 2c =K 2 −β2 ,entonces:

∂2 ∂2
( 2
∂x ∂y
+ 2
+ K c2 ⃗ )
h z( x, y) e−iβz =0

h z❑en dos dimensiones:

Factorizando con cero e−iβz queda demostrado ⃗

∂2 ∂2
( 2
∂x ∂y
+ 2
+ K c2 ⃗ )
h z( x, y)=0 Lqqd

Por el método de separación de variables asumiremos para la ecuación (19):

h z( x , y )=X x Y y (23)

Aplicando a la ecuación (22),

∂2 X x Y y ∂2 X x Y y 2
2
+ 2
+ K c X x Y y =0
∂x ∂y

∂2 X x ∂2 Y y 2
Yy 2
+ Xx 2
+ K c X x Y y =0
∂x ∂y
Dividiendo entre X x Y y a cada miembro, obtenemos:
2 2
1 ∂ Xx 1 ∂ Y y
. + . + K c2=0 (24)
X x ∂ x2 Y y ∂ y2

Como k c2=k x2 +k y 2, la ecuación (24), se puede expresar de la siguiente manera:

d2 X 2
+k x X=0 (25)
d x2
 d2Y 2
2
+k y Y =0 (26)
dy
Para la guía de onda cuadrada en el modo TE como se muestra en la ilustración 5.

Ilustración 5. Campo electromagnético propagándose a lo largo de Z.

Con las ecuaciones (25) y (26) se puede deducir lo siguiente.

h z( x , y )=( A cos(k x x)+ B sin(k x x ) ) . ( C cos( k y y)+ D sin( k y y ) ) (27)

Condiciones de la guía de onda cuadrada:

e x ( x , y )=0 , cuando Y =0 , b

e y ( x , y )=0 , cuando X =0 , a

Aplicando la ecuación (24) a la ecuación (20.c), se tendrá:

−i ∂ Hz
E x= 2
μω
Kc ∂y

−iμω
ex= k y ( A cos( k x x )+B sin (k x x) ) . (−C sin(k y y )+ Dcos (k y y ) ) (28)
K 2c

Análogamente para la ecuación (20.d), se obtendrá:

−iμω
e y= k x ( −A sin(k x x )+ B cos (k x x ) ) . ( C cos(k y y )+ D sin(k y y ) ) (29)
K 2c

Aplicando las condiciones de contorno de la guía de onda cuadrada.

Para e x ( x , y )=0 , cuando Y =0 , b. Aplicamos esta condición a la ecuación (28)
−iμω
ex= k y ( A cos( k x x )+B sin (k x x) ) . (−C sin(k y y )+ Dcos (k y y ) )=0
K 2c
Cuando Y=0:
D cos(k y 0)=0

D=0 (30.a)
Cuando y=b:
C sin (k y b)=0

sin( k y b)=0

sin( k y b)=sin ( nπ )

k y b=nπ

nπ
k y= , n=0,1,2 … (30.b)
b

Para e y ( x , y )=0 , cuando X =0 , a. Aplicamos esta condición a la ecuación (26)

−iμω
e y= k x ( −A sin ( k x x ) + B cos ( k x x ) ) . ( C cos ( k y y ) + D sin ( k y y ) )
K 2c

Cuando X=0:
B cos ( k x x ) =0

B=0 (30.c)
Cuando X=a:
A sin ( k x x ) =0

sin( k x a)=0

sin( k x a)=sin ( mπ )

k x a=mπ

mπ
k x= , m=0,1,2 … (30.d)
a

Dado que la propagación de la onda electromagnética es constante, de las

ecuaciones (30.b) y (30.d)

k c2=k x2 +k y 2
 mπ 2 nπ 2
k c2= ( )( )
a
+
b

mπ 2 nπ 2
k c=
√( a )( )
+
b
(31)

Como β 2=K 2 −K 2c

mπ 2 nπ 2

√
β= K 2− ( ) ( )
a
−
b
(32)

Como el nuero de onda está dado por K c =ωc √ με, despejando la frecuencia de
corte se tendrá:
K c =2 π f c √ με
Kc
f c=
2 π √ με

mπ 2 nπ 2

f c=
√( a )( )
+

2 π √ με
b (33)

Para el modo más bajo TE10, se denomina modo dominante, pues con n=0 y m=0
TEmn NO EXISTE, ya que E=0 y H=0.

Para frecuencias f <f c, 10 el modo NO SE PROPAGA (modo evanescente), existe

atenuación.

π 2 2, constante de atenuación.
α 10=
√( )a
−k

Para frecuencias f >f c, 10 el modo SE PROPAGA.

π 2, constante de propagación.
√
β 10= k 2− ()
a

Con las restricciones dadas en las ecuaciones (30.a), (30.b), (30.c) y (30.d), se
logra obtener h z( x , y ) .
mπ nπ
h z( x , y )=A cos( x). C cos ( y)
a b
 mπ nπ
h z( x , y )=A mn cos( x) cos ( y) (34)
a b
H z( x, y , z)=¿ ⃗
Como el campo magnético está dado por ⃗ h z(x , y) e−iβz,

mπ nπ
H z( x, y , z)= A mn cos (
⃗ x)cos( y)e−iβz (35)
a b
De las ecuaciones (17.a), (17.b), (17.c) y (17.d) calcularemos el campo eléctrico y
magnético para el modo TEmn.

−i ∂ Hz
E x= 2
μω
Kc ∂y

mπ nπ
∂( A mn cos ( x) cos( y) e−iβz )
−i a b
E x= μω
K 2c ∂y

iμω nπ mπ nπ
E x= 2 .
⃗ A mn cos ( x) sin( y ) e−iβz (36.a)
Kc b a b

Ey,⃗
Analogamente para ⃗ Hx y ⃗
H y , se tendrá:

−iμω mπ mπ nπ
E y= 2 .
⃗ A mn sin( x) cos( y)e−iβz (36.b)
K c a a b

iβ mπ mπ nπ
H x= 2 .
⃗ A mn sin( x)cos( y)e−iβz (36.c)
Kc a a b

iβ nπ mπ nπ
H y= 2 .
⃗ A mn cos( x) sin( y)e−iβz (36.d)
K c b a b

Impedancia
Ex ωμ
ZTE = = (33)
Hy β

2.4 PROPAGACIÓN MULTICAPAS

Las ecuaciones de Maxwell en el estudio de la propgación de ondas

electromagneticas en medios homogenioos o semi-hogenios con caracteristicas
definidas como la permitividad, la permitividad, la conductividad, etc, son
necesarias para saber el comportamiendo de la señal transmitida a lo largo del
medio [15]. Si la señal elecromagnetica se propaga de un medio 1 a un medio 2,
parte de esta señal propgada se transmite y la otra parte se refleja, estas señales
reflejada y transmitida dependen del coeficiente de refleción y el coeficiente de
transmisión respecivamente, tal como se muestra en la ilustración 6.

Figura 6. Paso de la señal transmitida de un medio hacia otro medio.

Si generalizamos el caso anterior para varios medios, es decir, la propagación de

la señal incidente a lo largo de multiples capas o multiples medios continuos, se
debe tener en cuenta mucho la impedancia en los medios [15] y el análisis del
parámetro S como la Transmitancia y Reflectancia en un punto de la cavidad,
donnde la cavidad hace referencia a la profundidad de los múltiples medios en
conjunto

Ilustración 7. Medios del suelo de múltiples capas

2.5 SOFTWARE DE SIMULACIÓN

Para la simulación de los prototipos de las antenas usaremos los dos softwares,
Ansys, HFSS y COMSOL. Continuación se da una breve descripción de cada una
de ellas.

2.5.1 COMSOL Multiphysics

COMSOL Multiphysics es una herramienta de modelado y análisis para simular

virtualmente de fenómenos físicos, como: transferencia de calor, movimiento de
fluidos, electromagnetismo y mecánica estructural, entre otros. Este software
trabaja sobre un sistema de Ecuaciones Diferenciales Parciales (PDEs), una
descripción matemática de varios fenómenos físicos basados en las leyes de la
ciencia [16].

2.5.2 HFSS, ANSYS

NSYS HFSS es un software de simulación electromagnética 3D de análisis

multifísico para diseñar y simular productos electrónicos de alta frecuencia como
antenas, conjuntos de antenas, componentes de RF o microondas, entre otros.
Este programa que trabaja con elementos finitos da solución a sistemas
estáticos/cuasi-estáticos. Este método se resuelve en una sola iteración, definida a
partir de intervalos grandes de tiempo, convergiendo fácilmente y mostrando una
gran estabilidad [17].
 6 CAPITULO III
7 000000
3. DISEÑO Y SIMULACIÓN

Se trabajó tres tipos de antenas para Georadares de penetración. A continuación,

se muestra el diseño y parámetros de cada antena.

3.1 ANTENA VIVALDI

La antena Vivaldi fue desarrollada por Gibson en 1979, desarrollada para producir
un haz simétrico en la salida y lóbulos laterales pequeños. Esta antena está
formada por una placa metálica con una apertura exponencial seguida de una
ranura uniforme terminada en una cavidad, sobre un sustrato. A continuación,
describimos las principales características que presenta la antena Vivaldi.

3.1.1 APERTURA EXPONENCIAL

La curva teórica de la antena Vivaldi, desarrollada por Schaubert, está dada por la
siguiente ecuación [18]:

Z=C1 e RY +C 2 (34)

Z −Z1
2
C 1= −R Y 2 (35)
e −e− R Y 1

Z 1 e−R Y −Z 2 e−R Y
2 1

C 2= −R Y 2 − R Y1
(36)
e −e

z1 ≤ z ≤ z2 y1 ≤ y ≤ y2 (37)

Donde R se define como la tasa de apertura de la apertura exponencial. Si

consideramos R=0, la forma de apertura sería lineal constante.
 Ilustración 8. Apertura exponencial [18].

Para tener una buena radiación de la antena, la longitud (L) y la anchura(H) de la

antena deben ser mayores que λ y λ/2 respectivamente. Con una anchura de la
apertura menor a λ0/2, la onda es guiada por la ranura; por el contrario, cuando la
anchura de apertura es mayor λ0/2, la onda se propaga por e dieléctrico hasta que
es radiada al exterior, tal como se muestra en la Ilustración 8.

3.1.2 TRANSICIÓN MICROSTRIP-SLOTLINE

Una transmisión slotline-microstrip consiste en una ranura gradaba en una cara

del sustrato, atravesara por una línea microstrip en circuito abierto que está
localizada perpendicularmente en la otra cara del sustrato, tal como se muestra en
la ilustración 9.

Ilustración 9 Transmisión slotline-microstrip [18].

Es importante tener un buen diseño de una transición microstrip-slotline para
conseguir un buen ancho de banda, pues una desventaja que presenta este tipo
de transición es la reducción del ancho de banda de operación. A continuación,
mencionaremos las principales características de las líneas slotline y microstrip.

 Slotline, elemento imprescindible para acoplar la potencia a la parte

radiante. Consta de una ranura uniforme que termina en una cavidad. Parte
de la radiación, se refiere a la forma exponencial de la apertura
 Slotline uniforme, forma parte de la transición y es donde se acopla la
potencia,
 La cavidad, también forma parte de la transición, su forma dependerá del
diseño de la misma.
 Microstrip, necesaria para la conexión con el resto del circuito y para
realizar la alimentación de la antena. Consta de un stub uniforme con
terminación radial.

3.1.3 DISEÑO DE ANTENA VIVALDI

Las dimensiones de la antena Vivaldi que se diseñó es de 525mm por 750mm, el

cual es aproximadamente d 0.35 λ por 0.5 λ, donde λ es la longitud de onda a una
frecuencia de 200MHz. La antena es alimentada por una transición de
microstrip/slotline. El sustrato usado es FR4 con un espesor de 1.6mm con
constante dieléctrica de 4.4 y una tangente de pérdida de 0.02.
 Ilustración 10. Antena Vivaldi para GPR.

R1 y R2 corresponden a los factores exponenciales que están definidas por las

coordenadas (0, Ds/2+Ls, -W1/2), (0, Lr1,-W1/2) y (0, Lsub-Lr0,-Wsub/2) tal como
se muestra en la ilustración 10.

Las dimensiones de la antena Vivaldi planteada en diseño son lo siguiente:

Tabla 1. Dimensiones de la Antena Vivaldi

DIMENSIÓN VALOR (mm)

Lsub 750

Wsub 525

Lr0 50

Lr1 600

Wr1 82

Ds 88
 Rsub 69

Dfd 50

Ls 94
W1 6
Wstripline 3

Theta 110°

Ldgs 100

Wdgs 90

Gdgs 30

G_dgs 32

Dgs 236

P 13

El ángulo que presenta la sección de microstrip está definida por theta. Las
estructuras con forma de T se agregaron al diseño original de Vivaldi con la
finalidad de mejorar el rendimiento de la antena, tal como se aprecia en la
ilustración 10 y los datos se muestran en la tabla 1.

3.1.3.1 RESULTADOS ANTENA VIVALDI

La simulación fue realizada con Ansys,HFSS. Los resultados que se muestran son
e patrón de radiación, ancho de haz, ganancia y ancho de banda de la antena.

3.1.3.1.1 PATRÓN DE RADIACIÓN

En la ilustración 11a se observa el patrón de radiación que presenta la antena

Vivaldi a 240MHz presentando el lóbulo principal con una ganancia máxima de
6.85dB y dos lóbulos secundarios con un máximo de 0.4dB.
 Ilustración 11. A) Patrón de Radiación 3D y (b)patrón de radiación en el plano E y H

En el plano H, ilustración 11b, se aprecia el ancho de haz-PHBW de 70°.También

presenta F/B, es la diferencia entre la ganancia máxima y la ganancia back del
lóbulo principal, de -16dB.

3.1.3.1.2 PARÁMETRO S11

El ancho de banda que trabaja la antena Vivaldi va desde 239MHz hasta 256MHZ
para VHF.

Ilustración 12. Parámetro S11

 Figura 2. (a) Patrón de radiación de la antena Vivaldi, (b) parámetro S11 y (c) ancho de banda de la antena Vivaldi.

También presentan un rango de frecuencia de trabajo de 22MHz por debajo de

frecuencias menores a 130MHZ. En la ilustración 12 se aprecia el parámetro S11.

3.1.4 DISEÑO DEL ARREGLO DE ANTENAS VIVALDI

Para la doble polarización con Antenas Vivaldi, ya que LA Antena Vivaldi tiene
polarización lineal, se define una estructura en la cual trabajarán en pares que
tengan la misma polarización orientado en el espacio, es decir frontalmente entre
emisor y receptor, donde los dos emisores estarán rotados tal que sus
polarizaciones sean ortogonales, formando las paredes laterales de un cubo. Para
mantener fijo la estructura de antenas, se decidió usar un marco que sujetara las
antenas de material de PVC el cual no tendría un gran efecto en el
comportamiento de la radiación de las antenas Vivaldi, el cual se puede observar
en la figura 13.

Ilustración 13. Arreglo de 4 antenas Vivaldi para GPR.

 3.1.4.1.1 RESULTADOS DEL ARREGLO DE ANTENAS VIVALDI

Mediante simulaciones se obtuvieron resultados favorables. En la ilustración 13 se

observa que el arreglo de cuatro antes para GRP trabaja en un rango de
frecuencias entre 219MHZ a 256MHz, con un ancho de banda 40MHz
aproximadamente.

Ilustración 13 Ancho de banda del arreglo de antena Vivaldi.

Del parámetro s11 también se puede observar el desacoplamiento que presenta

este arreglo de antes, el cual es muy favorable con los objetivos que se buscaba.
Esto quiere decir que los valores de acoplamiento están por debajo de -20dB,
como se observa en la ilustración 14.
 Ilustración 14. Acoplamiento del arreglo de Antena Vivaldi.

3.2 ANTENA LPDA

Este tipo de antena está formada por un conjunto de radiales de distinto tamaño,
ordenados de forma descendente en dirección a la máxima radiación de la antena.
Las longitudes de los elementos directores deben ser menores a los elementos
anteriores. La característica de este tipo de antena es que la impedancia de
entrada, la ganancia y el diagrama de radiación de la antena logarítmica periódica
se repite periódicamente.

Los elementos de la antena varían incrementado a medida que se reduce la

distancia al origen de la antena, suponiendo que el origen está en la parte
izquierda de la antena.

La antena presenta una relación logarítmica entre elementos, esta relación es

representada por la constante de diseño o factor de separación τ [19]:
L n+1 wn+1 Dn+1
τ= = = (38)
ln wn Dn

Otro factor que es característico de la antena LPDA es el factor de espaciamiento (

σ ) entre elementos:
 dn
σ= (39)
2l n

σ =0.25(1−τ )cot α (40)

Donde:
Ln, longitud de los elementos en el orden n.

D n, es la distancia entre los elementos n y n+1.

w n, es el ancho de los elementos en el orden n.

α , define el contorno-perfil de la antena LPDA.

Ilustración 15. Ganancia de los arreglos de dipolos log-periódicos [19].

En la ilustración 15, muestra el cuadro de la variación de la ganancia de la antena

LPDA, el cual depende de los parámetros de diseño τ y σ .

3.2.1 DISEÑO DE ANTENA LPDA

La antena LPDA fue diseñada y optimizada con el objetivo de que trabaje en el
rango de 300MHZ a 500MHz, VHF/UHF. El diseño se muestra en la ilustración 16.

Ilustración 16. Antena LPDA

Optimizando los parámetros τ =0.78 y σ =0.14 con la finalidad de mejorar el ancho

de banda y tener una ganancia por encima de los 5dB se presenta la tabla 2 con
los valores de las dimensiones de la antena LPDA.

Tabla 2. Dimensiones de antena LPDA

DIMENSIONES VALOR (mm)

Yt 750.34
Xt 525
P0 20
PP 119
D1 127
L1 525
 W1 17.8
C0 85
Xw 8.3
α 21.4977°

3.2.1.1 RESULTADOS DE LA ANTEN LPDA

Los resultaos que obtenidos de la simulación mostrados en la ilustración 17a un

buen patrón de radiación 3D, esta presenta una ganancia máxima de 5.12dB y
front to back- F/B de 11.9dB. En la ilustración 17b se visualiza el patrón de
radiación del plano E y plano H, observando el plano H se obtiene el valor del
ancho de haz que es de 62°.

Ilustración 17. A) Patrón de Radiación 3D y (b) Plano H y Plano E

Por otro lado, en la ilustración 18 se observa el ancho de banda de trabajo de la

antena LPDA es de 341MHZ, los resultados obtenidos son muy buenos debido
objetivo es diseñar una antena de tipo LPDA que trabaje en el rango de 300MHz
hasta 500MHZ.
 Ilustración 18. Parámetro S11

3.2.2 DISEÑO DEL ARREGLO DE ANTENAS LPDA

Debido a que la Antena LPDA tiene polarización lineal, se define una estructura
que estén frontalmente emisor y receptor, donde los dos emisores estarán rotados
tal que sus polarizaciones sean ortogonales, formando las paredes laterales de un
cubo. Se usa un marco que sujetara las antenas de material de PVC para
mantenerlas fijas, como se observar en la figura 19.
 Ilustración 19. Arreglo de 4 antenas LPDA para GPR.

3.2.2.1 RESULTADOS DEL ARREGLO DE ANTENAS LPDA

Mediante simulaciones se obtuvieron resultados favorables. En la ilustración 20 se

observa que el arreglo de cuatro antes para GRP trabaja en un rango de
frecuencias entre 219MHZ a 256MHz, con un ancho de banda 40MHz
aproximadamente.
 Ilustración 20. Ancho de banda del arreglo de antena Vivaldi.

En la ilustración 21 se observa el buen desacoplamiento que presenta este arreglo

de antenas LPDA, es decir que los valores de acoplamiento están por debajo de
-20dB.

Ilustración 21. Acoplamiento del arreglo de Antena LPDA.

3.3 ANTENA VIVALDI ANTIPODAL

 3.3.1 DISEÑO DE ANTENA VIVALDI ANTIPODAL

En la optimización del diseño de la antena Vivaldi Antipodal se implementó un

Algoritmo Genético Multiobjetivo-MOGA. Este algoritmo [20]. El algoritmo esta
descrita por:

a) Codificación
b) Fitness and selection probability (SP)
c) Crossover
d) Mutation
e) Elitism

Este algoritmo toma los valores de los bordes de la antena, y mediante el

algoritmo optimiza estas dimensiones. En la ilustración 22, se aprecia el diagrama
de bloques que muestra el proceso del desarrollo del algoritmo.
 Ilustración 22. Diagrama de bloques del proceso de MOGA [20]

Para el caso de la optimización de la Antena Vivaldi Antipodal, se aplicó el MOGA,

de tal manera que este algoritmo optimiza los bordes de la geometría, con la
finalidad de mejorar el ancho de banda de la antena. Adicionalmente al diseño se
agregó dos reflectores paralelos a la antena para mejora el patrón de radiación,
como se muestra en la ilustración 23.
 Ilustración 23. Antena Vivaldi Antipodal

En la tabla 3, se presenta los valores de las dimensiones de la antena LPDA que

fueron obtenidos por el algoritmo genético.

Tabla 3. Dimensiones de antena Vivaldi Antipodal

DIMENSIONES VALOR (mm)

Yt 713.3
Xt 525
P0 225.57
P1 86.23
P2 78.5
L0 259.8
h 35
A0 200
A1 357.1
 T 148.6

3.3.1.1 RESULTADOS DE ANTENA VIVALDI ANTIPODAL

Los resultaos que obtenidos de la simulación mostrados en la ilustración 24,

presenta una ganancia máxima de 6.23dB y F/B de 12dB en el plano E, del plano
H se calcula el valor del ancho de haz que es de 60°.

Ilustración 24. Plano H y Plano E del patrón de radiación.

Por otro lado, en la ilustración 25 se observa el ancho de banda de trabajo de la

antena Vivaldi antipodal es de 361MHz, los resultados obtenidos son muy buenos
debido a que el objetivo es diseñar una antena que trabaje en el rango de los
300MHz a 500MHZ.
 Ilustración 25. Parámetro S11

3.3.2 DISEÑO DEL ARREGLO DE ANTENA VIVALDI ANTIPODAL

Debido a que la Antena Vivaldi Antipodal tiene polarización lineal, se define una
estructura que estén frontalmente emisor y receptor, donde los dos emisores
estarán rotados tal que sus polarizaciones sean ortogonales, formando las
paredes laterales de un cubo. Se usa un marco que sujetara las antenas de
material de PVC para mantenerlas fijas, como se observar en la figura 26.

Ilustración 26. Arreglo de 4 antenas LPDA para GPR.

3.3.2.1 RESULTADOS DEL ARREGLO DE ANTENAS LPDA

Los resultados de obtenidos del arreglo de antena Vivaldi Antipodal no fueron lo
esperado, pues el rango de frecuencia que irradia es de 31MHz que va desde
420MHz a 457MHz, tomando en cuento valores menores a los -10dB con la
finalidad de tener poca perdida de retorno.

Ilustración 27. Ancho de banda del arreglo de antena Vivaldi.

 8 CAPITULO IV

4 SISTEMA GEORADAR-SUBSUELO

El sistema Georadar-Subsuelo esta compuesto por un arreglo de antes de GPR,

se habla en el capítulo III, y el subsuelo. El arreglo de antena está envuelto por el
aire seguido del subsuelo y rodeado por una capa de PML, tal como se muestra
en la ilustración 28.

Ilustración 28. A) Vista frontal del sistema arreglo de antena-subsuelo y (b) el sistema visto desde la parte superior.

El subsuelo consta de once tipos de medios con diferente coeficiente de

permitividad, en la tabla 3 se observa los diferentes tipos de medios y sus
características electromagnéticas. Los tipos de materiales fueron obtenidos del
“ESTUDIO DE SUELOS, CANTERAS Y FUENTES DE AGUA” [21].

Tabla 4. Tipos de Materiales

Permitividad Permeabilidad
N° Tipo de Medio Conductividad
Relativa Relativa
Aire Aire 1 1 0
M1 Arena arcillosa semidenso 6 1 0.001
M2 Grava arenosa 15 1 0.0001
M3 Grava arenosa 16 1 0.0002
M4 Grava arenosa 17 1 0.0003
M5 Suelo limoso seco semidenso 3 1 0.001
M6 Grava arenosa 18 1 0.00045
M7 Grava arenosa 19 1 0.0005
 M8 Arena limoso saturado 13 1 0.0012
M9 Grava arenosa 20 1 0.00065
M10 Arena arcillosa densa 9.5 1 0.001
M11 Grava arenosa 22 1 0.00075

4.1 ANÁLISIS PREVIO

En primera instancia se tenía la confusión en la aplicación de las condiciones de

frontera, la opción de aplicar open región o PML en la frontera con la finalidad de
obtener buenos resultados evitando la reflexión en los bordes de las capas que
simulara regiones que se extienden al infinito. A continuación, mencionaremos una
breve comparación teórica de estos:

Una open región o región abierta es un volumen que rodea un modelo abierto en
el que se calculan los campos cercanos en el vacío [21]. Las caras externas de la
región tienen límites de radiación PML aplicados, y estos límites absorben las
ondas electromagnéticas salientes, tal como se muestra en la Ilustración 29a.

Ilustración 29. A) Open region y (b) PML

La capa perfectamente combinada (PML) son materiales ficticios que absorben los
campos electromagnéticos [21], PML es utilizado para crear un modelo abierto en
las caras extremas del espacio de solución. PML son más apropiadas que los
límites de radiación-open region para las simulaciones de antena, debido a que
pueden reducir el volumen de la solución, ya que se pueden colocar a una
distancia 1/10 veces la longitud de onda cerca de cualquier estructura radiante,
como muestra la ilustración 29b.

En primera instancia simulamos el sistema Radar-subsuelo, donde el subsuelo

está compuesto por 4 diferentes capas de materiales, como se en la ilustración 30.

Ilustración 30. A) sistema con open region y (b) sistema con PML

Al realizar la simulación electromagnética de la geometría descrita en la ilustración

31, se obtuvo los siguientes resultados:

Ilustración 31. A) S11 del sistema con open region y (b) S11 del sistema con PML

En la ilustración 31, se observa el parámetro de reflexión que radia en desde los

213 MHz hasta 255MHz. Los resultados varían ligeramente con una profundidad
mínima en el sistema de región abierta de -18dB y en el sistema con PML muestra
un mínimo de -15dB.
 Ilustración 32. A) acoplamiento lateral del sistema con open region y (b) acoplamiento lateral del sistema con PML

El sistema con open region presenta acoplamiento en el rango de frecuencia de

trabajo de 220MHz-221MHz como se muestra en la ilustración 32a, teniendo en
cuenta que para tener un buen desacoplamiento tomamos valores por debajo de
-20dB. Por otro lado, el sistema con PML presenta un buen desacoplamiento en
las presencias de trabajo que va de 213MHz hasta 255MHz, tal como se muestra
en la ilustración 32b.

Ilustración 33. A) acoplamiento dual del sistema con open region y (b) acoplamiento dual del sistema con PML

En la ilustración 33a, se observa el acoplamiento que presenta el sistema –open

region en el rango de frecuencia objetivo que debe trabajar el sistema. Por otro
lado, el sistema –PML presenta buen desacoplamiento por debajo de los -20dB,
como se observa en la ilustración 33b.

Los resultados presentados anteriormente se complementan con el análisis del

campo eléctrico a lo largo de los medio Aire, M1, M2 y M3, tal como se muestra en
la ilustración 34.
 Ilustración 34. Magnitud del campo eléctrico a lo largo de la cavidad.

La curva roja es la manigtud del campo electrico del sistema con open region y la
curva azul es la magnitud del campo electrico del sistema con PML, se observa que
ambas presentan atenuación al pasar de un medio a otro debido a las caracteriscas
de los medios.

La magnitud del campo electrico de la curva roja en el medio M3 aumenta debido a

la reflexipon de la señal en el borde, el cual no asemeja a superficies infinitas. Esto
es debido a que se trabaja con materiales con caracteristicas diferentes al vacio.

Por otro lado, la magnitud del campo eléctrico de la curva azul no presenta un
incremento significativo debido a que en los bordes no existe señal que regresa, es
decir, no esxites reflexión de la señal con lo cual simula superficies infinitas, que es
lo que se busca en el análisis.

4.2 ANÁLISIS DE PÉRDIDAS DE RETORNO

La ilustración 35 muestra el parámetro S11, el rango de frecuencia de trabajo del

sistema va desde los 216MHz a 254MHz.
 Ilustración 35. Parámetro S11 del sistema.

Para evitar el acoplamiento en un arreglo de antenas debemos tener pérdidas de

retorno bajas, es decir que los valores del parámetro S11 deben ser menores a los
-20dB.

Ilustración 36. A) S11 de la antena frontal y (b) S11 de antenas laterales.

La ilustración 36, muestra valores menores a -20dB para el rango de frecuencia

que trabaja el sistema, por tanto, hay desacoplamiento en el sistema simulado.

4.3 PATRÓN DE RADICACIÓN Y GANANCIA DEL ARREGLO

 Ilustración 37 a) Plano E, b) Plano H y (c) Patrón de radiación

El patrón de radiación del arreglo de antenas en el sistema planteado, como se

observa en la ilustración 37, en el Plano H presenta una reducción del 20%, esto
es debido a las ondas electromagnéticas que se reflejan de las multicapas
generando un aumento en el back de 3dB.

4.4 COMPORTAMIENTO DEL CAMPO ELECTRICO A LO LARGO DE LA

CAVIDAD
 Ilustración 38. Variación del campo eléctrico a lo largo de la cavidad- multicapas.

La variación de la magnitud del campo eléctrico a lo largo de la cavidad o múltiples

capas se muestra en la ilustración 38.

En la cavidad se presenta atenuación debido a los diferentes valores de

conductividad de cada medio, como se observa en la tabla 3. También se observa
un ligero incremento en la magnitud del campo eléctrico en medios cuyo
coeficiente de reflexión es menor, se visualiza en medios M5 y M10 con
coeficiente de permitividad relativa de 3 y 9.5 respectivamente.

En los primeros cuatro medios los campos transmitidos presentan ángulos

incidentes pequeños. El campo, al pasar del medio 5 al medio 6 presenta ángulo
de reflexión pequeña y más directivo al ángulo incidente generando un incremento
en la magnitud del campo total en el medio 5, de igual manera sucede en el medio
10.

El arreglo de cuatro antenas se encuentra en el aire, debido a esto las ondas

reflejadas por las multicapas varía el patrón de radiación y ganancia del arreglo de
cuatro antenas para GPR. Por otro lado, la magnitud campo disminuye a lo largo
de la cavidad, llegando a 10 metros con valores de 1,3 V/M.
 9 CAPITLO V

5 ANÁLISIS COMPARATIVO

En este capítulo se realiza un análisis comparativo de los resultados obtenidos en

el capítulo tres, para esto lo dividiremos en tres partes:

 Analisis comparativo de antenas

 Analisis comparativo de arrgelos de antena
 Analisis comparativo de sistemas Georadar-subsuelo.


5.1 ANÁLISIS COMPARATVO DE LOS TRES TIPOS DE ANTESN PARA GPR

5.1.1 ANCHO DE BANDA DE LAS ANTENAS

Ilustración 39. Ancho de banda de las tres antenas.

La ilustración 39, muestra el rango que irradian las antenas. La curva roja
representa la antena Vivaldi trabaja en el rango de frecuencia de 239MHz a
257MHz. La antena LPDA presenta un ancho de banda mayor a los 200MHZ que
va desde 250MHZ a 530MHz, por último, la antena Vivaldi Antipodal presenta
ancho de banda mayos a los 200 MHZ que va desde 250MHZ hasta 510MHZ,
teniendo en cuenta valores de S11 menores a -10dB para una baja perdida de
retorno.

9.1.1.1.1
5.1.2 PATRÓN DE RADICACIÓN Y GANANCIA DE LAS ANTENAS

Ilustración 40. a) Plano E y b) Plano H

El ancho de haz- PHBW se visualiza en la ilustración 40a. La antena Vivaldi

presenta un PHBW de 70°, la antena LPDA con 67° y la antena Vivaldi Antipodal
con 60°.

Tabla 5. Cuadro comparativo de las tres antenas

A. Vivaldi
A. Vivaldi B. LPDA
Antipodal
Ancho de banda (dB) 40MHz 277MHz 260MHz
HPBW 70° 65° 60°

Ganancia (dB) 6.84 5.56 6.23

 F/B (dB) -15 -11 -

El cuadro comparativo de la tabla 4 presenta un resuenen de los parámetros de

cada antena, la antena que más se adecua a los objetivos es la antena LPDA,
pues presenta un ancho de banda mayor a los 200MHz, F/B de 11bB y un buen
ancho de haz.

5.2 ANÁLISIS COMPARATIVO DE LOS TRES ARRGELOS DE ANTENAS

PARA GPR

5.2.1 ANCHO DE BANDA DE ARREGLOS DE ANTENA

Ilustración 41. A) S11 de la antena Vivaldi y b) S11 de la antena LPDA y antena Vivaldi Antipodal.

El ancho de banda del arreglo de antena Vivaldi es de 40MHz que va de los

219MHz a 256MHz, tal como se muestra en la ilustración 41a.

En la ilustración 41b muestra el ancho de banda de las antenas LPDA y Vivaldi

Antipodal, los valores para la antena LPDA va de los 250MHz a 550MHZ, y de la
antena Vivaldi Antipodal va de 420MHZ a 458MHZ, teniendo en cuenta valores de
S11 menores a -10dB para una baja perdida de retorno.
 5.2.2 ACOPLAMIENTO EN LOS ARREGLOS DE ANTENA

Ilustración 42. A) Lat de la antena Vivaldi, b) last. De la antena LPDA y c)

Para evitar el acoplamiento en un arreglo de antenas debemos tener pérdidas de

retorno bajas, es decir que los valores del parámetro S11 deben ser menores a los
-20dB. La ilustración 42a, 42b y 42c muestra que, en los tres arreglos de antenas,
esta presenta un buen acoplamiento en ellas.

10
 11 CONCLUSIONES Y RECOMENDACIONES

En el capítulo II se realizó un análisis comparativo de las tres antenas radiada

en el vacío, de acuerdo a los resultados obtenidos y al objetivo que se planteó,
la antena LPDA es la que mejores resultados presentes, pues el ancho de
banda en que trabaja la antena es mayor a los 200MHz que va de los 300MHz
a los 500MHz, ancho de haz de 65° y el lóbulo principal presenta una ganancia
mayor a los 5dB. Por otro lado, la antena Vivaldi presenta el mejor F/B y mejor
ganancia de 6.8dB, el ancho de banda es solo de 40MHz.

Para el caso de análisis comparativo del arreglo de antena para GPR, el

arreglo más óptimo que cumple con los objetivos planteados es el arreglo de
antenas LPDA. Este arreglo presenta resultados de ancho de banda mayor a
los 200MHZ y un buen desacoplamiento en el rango de frecuencia que va de
los 300MHz a 500MHZ. El segundo arreglo es de la Vivaldi, esta trabaja en el
rango de frecuencia de 40MHZ. El ultimo arreglo es de Vivaldi Antipodal, no se
logró los objetivos en el ancho de banda y también presenta acoplamiento
entre las frecuencias que va de 300MHZ a 500MHz.

Para el caso del sistema-subsuelo, se puede concluir que la presencia de las

multicapas afecta 0.02% al ancho de banda del arreglo cuatro antenas para
GPR en el caso de la antena Vivaldi comparada con el arreglo de antenas
Vivaldi radiada en el vacío. También presenta un buen desacoplamiento
tomando valores menores a -20dB del parámetro S11. El patrón de radiación
se ve afectado presentando un incremento en el back de +3dB, generando un
mayor valor en el F/B de -0.5dB.
 12 BIBLIOGRAFÍA

[1] Mamani R.M. et al. (2005) “The Application of Ground Penetrating Radar (GPR)
at Machu Picchu, Peru (C101-1)”. In: Sassa K., Fukuoka H., Wang F., Wang G.
Landslides. Springer, Berlin, Heidelberg. https://doi.org/10.1007/3-540-28680-
2_5

[2] Pachacamac-old wak’a or Inca syncretic deity? Imperial transformation of

sacred landscape in the Lower Ychsma (Lurín) Valley, T. Bray, “The
Archaeology of wak’as: Explorations of the sacred in the pre-Columbian Andes”,
127-166, University of Florida Press, Gainesville.

[3] Jan Igel and Kazanori Takahashim, “Influence of antenna polarization on GPR
signals of linear targets” DGG Tangung 2014

[4] Churio, Julio; Ochoa, Ángel; Morantes, Zoraida "SIMULACIÓN DE LA

PROPAGACIÓN DE ONDAS ELECTROMAGNÉTICAS EN AMBIENTES
CERRADOS UTILIZANDO FRONTERAS ABSORBENTES TIPO PML EN 3D"
Télématique, vol. 9, núm. 1, 2010, pp. 76-97

[5] W.C. Chew; B. Hu; Y.C. Pan; J.S. Zhao “Fast algorithm for complex structures”
APMC 2001. 2001 Asia-Pacific Microwave Conference (Cat. No.01TH8577), doi:
10.1109/APMC.2001.985592, 2001.

[6] Y. Brand; A. Alvarez-Melcon; J.R. Mosig; R.C. Hall “Large distance behavior of
stratified media spatial Green's functions”, IEEE Antennas and Propagation
Society International Symposium 1997. Digest, doi: 10.1109/APS.1997.625437,
1997.

[7] Xianyang Gao, Frank J. W. Podd, Wouter van Verre, David J. Daniels, and
Anthony J. Peyton "Investigating the Performance of Bi-Static GPR Antennas for
Near-Surface Object Detection",IEEE,doi: 10.3390/s19010170, 2019.

[8] Carlos Jara Bustamante “Exploración Geofísica mediante la técnica del radar de
penetración terrestre. Un estado de conocimiento”, tesis de grado, 2015.

[9] Almeida I. R., Borges W. R., Lanniruberto M., Bittencourt A. C. and Cunha L. S.,
et. al. “Modelagem GPR de Canais Fluviais Atuais”. Twelfth International
Congress of the SBGf, Sociedade Brasileira de Geofísica, Brasil, 2011.
[10] Daniels David J. (2004), “Ground Penetrating Radar” Institution of Engineering
and Technology, London, second edition.

[11] Karthikeyan.R, Chandramouli A and G. Srivatsun “Ground Penetrating Radar

(GPR) Antenna Design: A Comparative Study”, vol. 8, pp. 2249-8958, 2018

[12] Constantine A. Balanis “advanced engineering electromagnetics”, 2nd edición.

[13] Antenas, 3Cu Electrónica. Retrieved from

https://sites.google.com/site/3cuelectronica/home/radio-enlaces-1/antenas

[14] Steven Rojas, Polarización de Ondas Electromagnéticas. Retrieved from

https://medium.com/@steven98.sr/polarizaci%C3%B3n-de-ondas-electromagn
%C3%A9ticas-67bce51a203d

[15] David M. Pozar “MICROWAVE ENGINEERING” 4ta edición.

[16] ANSYS. (2020). ANSYS HFSS: High Frequency Electromagnetic Field

Simulation Software. Retrieved from
https://www.ansys.com/products/electronics/ansys-hfss

[17] COMSOL Multiphysics. Retrieved from https://www.comsol.com

[18] Lizette Lorraine Tovar “DISEÑO DE UNA ANTENA VIVALDI ANTIPODAL

MODIFICADA PARA APLICACIONES GPR EN LA BANDA DE 700 MHz - 9
GHz.”,2014.

[19] Marek Bugaj, Marian Wnuk, Rafal Przesmycki “Microstrip broadband LPDA
antenna on VHF and UHF bands”, 2nd URSI AT-RASC, Military University of
Technology, Warsaw, Poland, 2018.

[20] Marco Pumallica, Jorge Arizaca, Mark Clemente “A Multiobjective Genetic

Algorithm for Analysis, Design and Optimization of Antipodal Vivaldi Antennas”,
Published in: 2019 IEEE-APS Topical Conference on Antennas and
Propagation in Wireless Communications (APWC) IEEE IDO:
10.1109/APWC.2019.8870420, 2019.

[21] Estudio de suelos, canteras y fuentes de agua, MTC-Provincia descentralizado

municipal distrital de Super, Lima-Perú.