setwd("/cloud/project/datos")
datos <- read.csv("Petroleo_Ontaro.csv", header=T, dec=".", sep=";")
elevacion<-datos$GROUND_ELEVATION
elevacion<-na.omit(elevacion)
# GDF
min(elevacion)
## [1] 0.3
max(elevacion)
## [1] 694
R_elevacion<-max(elevacion)-min(elevacion)
k_elevacion<-1+(3.3*log10(length(elevacion)))
k_elevacion<-floor(k_elevacion)
A_elevacion<-R_elevacion/k_elevacion
HistoSturges_elevacion<-hist(elevacion, main="Gráfica No.1:
Distribución de Elevación en el pozo",
breaks = seq(min(elevacion),max(elevacion),A_elevacion),
xlab="Elevación (m)",
ylab="Cantidad",
col="red")
HistoSturges_elevacion$counts
## [1] 37 222 88 7695 14491 1481 420 104 16 27 14 4
## [13] 0 1 1
Histograma_elevacion<-hist(elevacion, main="Gráfica No.2:
Distribución de Elevación en el pozo",
xlab="Elevación (m)",
ylab="Cantidad",
col="red")
limites_elevacion<-Histograma_elevacion$breaks
liminf_elevacion<-limites_elevacion[1:15]
limsup_elevacion<-limites_elevacion[2:16]
MC_elevacion<-(liminf_elevacion+limsup_elevacion)/2
# TABLA ELEVACIÓN
ni_elevacion<-c()
for(i in 1:k_elevacion) {
if(1==15)
ni_elevacion[i]<- length(subset(elevacion, elevacion>=liminf_elevacion[i] & elevacion<=limsup_elevacion[i]))
else
ni_elevacion[i]<- length(subset(elevacion, elevacion>=liminf_elevacion[i]& elevacion<limsup_elevacion[i]))
}
hi_elevacion<-(ni_elevacion/length(elevacion))*100
Niasc_elevacion<-cumsum(ni_elevacion)
Hiasc_elevacion<-cumsum(hi_elevacion)
Nidsc_elevacion<-rev(cumsum(rev(ni_elevacion)))
Hidsc_elevacion<-rev(cumsum(rev(hi_elevacion)))
Tabla_elevacion<-data.frame(liminf_elevacion,
limsup_elevacion,
MC_elevacion,
ni_elevacion,
round(hi_elevacion,2),
Niasc_elevacion,
round(Hiasc_elevacion,2),
Nidsc_elevacion,
round(Hidsc_elevacion,2))
Tabla_elevacion
## liminf_elevacion limsup_elevacion MC_elevacion ni_elevacion
## 1 0 50 25 49
## 2 50 100 75 239
## 3 100 150 125 70
## 4 150 200 175 14479
## 5 200 250 225 8893
## 6 250 300 275 545
## 7 300 350 325 235
## 8 350 400 375 39
## 9 400 450 425 27
## 10 450 500 475 18
## 11 500 550 525 5
## 12 550 600 575 0
## 13 600 650 625 1
## 14 650 700 675 1
## 15 700 NA NA 0
## round.hi_elevacion..2. Niasc_elevacion round.Hiasc_elevacion..2.
## 1 0.20 49 0.20
## 2 0.97 288 1.17
## 3 0.28 358 1.46
## 4 58.86 14837 60.31
## 5 36.15 23730 96.46
## 6 2.22 24275 98.67
## 7 0.96 24510 99.63
## 8 0.16 24549 99.79
## 9 0.11 24576 99.90
## 10 0.07 24594 99.97
## 11 0.02 24599 99.99
## 12 0.00 24599 99.99
## 13 0.00 24600 100.00
## 14 0.00 24601 100.00
## 15 0.00 24601 100.00
## Nidsc_elevacion round.Hidsc_elevacion..2.
## 1 24601 100.00
## 2 24552 99.80
## 3 24313 98.83
## 4 24243 98.54
## 5 9764 39.69
## 6 871 3.54
## 7 326 1.33
## 8 91 0.37
## 9 52 0.21
## 10 25 0.10
## 11 7 0.03
## 12 2 0.01
## 13 2 0.01
## 14 1 0.00
## 15 0 0.00
colnames(Tabla_elevacion)<-c("Xi-1","Xi+1","Xi (MC)","ni","hi","Ni ↓","Hi ↓","Ni ↑","Hi ↑")
# OJIVAS CONTINUAS ELEVACIÓN
plot(limsup_elevacion, Niasc_elevacion, col="blue", type="b",
main="Gráfica No.3: Distribución de Elevación en el Pozo",
ylab = "Cantidad",
xlab = "Elevación (m)")
lines(liminf_elevacion, Nidsc_elevacion, col="red", type="b")
# OJIVAS CONTINUAS PORCENTAJE ELEVACIÓN
plot(limsup_elevacion, Hiasc_elevacion, col="blue", type="b",
main="Gráfica No.4: Distribución de Elevación en el Pozo",
ylab = "Porcentaje",
xlab = "Elevación (m)")
lines(liminf_elevacion, Hidsc_elevacion, col="red", type="b")
# GDF CAJA Y BIGOTES
boxplot(elevacion,
horizontal=T,
col="darkgreen",
xlab= "Profundidad Alcanzada (m)",
main="Gráfica No.5:
Distribución de Elevación en el Pozo")
# Conclusion
conclusion <- "El comportamiento de la elevación del terreno fluctúan entre [0.3, 694] y los valores se encuentran entorno a 193.85, con una desviación estándar de 31.21711 siendo un conjunto de valores homogéneos donde los valores se agrupan en la parte izquierda de la variable y tiene un sesgo positivo con existencia de valores atípicos a partir de 250."