In khaors/rves: R package for the interpretation of Vertical Electric Soundings

# The `params` object is available in the document.
library(rves)
library(knitr)
library(grid)
library(gridExtra)
library(kableExtra)

Informacion Basica del Sondeo

Las mediciones realizadas en el sondeo r params$current.ves$id[[1]] estan conformadas por r length(params$current.ves$ab2) valores tal y como se presenta en la siguiente tabla:

spacing <- params$current.ves$ab2
apprho <- params$current.ves$appres
ves.df <- data.frame("AB/2" = spacing, "Resitividad.Aparente" = apprho)
kable(ves.df, format = "html", digits = 3, align = c('c','c'), caption = "SONDEO ELECTRICO VERTICAL") %>%
  kable_styling(bootstrap_options = "striped", full_width = F)

Los valores minimo y maximo del espaciamiento de electrodos son r min(params$current.ves$ab2) and r max(params$current.ves$ab2) metros respectivamente, y algunas de las estadisticas basicas de la resistividad aparente se presentan en la siguiente tabla:

summary(params$current.ves$appres)

Interpretacion de Resultados

El grafico de las mediciones realizadas en el sondeo electrico vertical y su correspondiente interpretacion se muestran en la siguiente figura:

plot(params$current.ves, type = "ves")

rho.true <- params$current.ves$rhopar
thickness <- params$current.ves$thickpar
interpretation.ves <- cbind(rho.true, thickness)
interpretation.ves.df <- as.data.frame(interpretation.ves)
names(interpretation.ves.df) <- c("Resistividad.Real", "Espesor")
nlayers <- length(rho.true)
layers.id <- vector("character", length = nlayers)
for(i in 1:nlayers){
  layers.id[i] <- paste0("Capa", as.character(i))
}
row.names(interpretation.ves.df) <- layers.id

La interpretacion de este sondeo requiere de r nlayers capas con las siguientes caracteristicas:

kable(interpretation.ves.df, format = "html", digits = 3, align = c('r','r'), 
      caption = "MODELO DE CAPAS") %>%
  kable_styling(bootstrap_options = "striped", full_width = F)

Diagnostico del Modelo

La calidad del modelo estimado para este sondeo electrico vertical se puede evaluar de forma cualitativa por medio de la siguiente figura:

p1 <- plot(params$current.ves, type = "diagnostic")
grid.draw(p1)

khaors/rves documentation built on May 29, 2019, 3:07 a.m.