In f0nzie/rNodal: Oil and Gas Well Nodal Analysis for Production Petroleum Engineering

library(rNodal)

well <- "
    MD  TVD
    0       0
    600     600
    1005    1000
    4075    4000
    7700    7500
    9275    9000
"

# read string text to dataframe
well_table <- WellDeviationSurvey(well, reference = "vertical")
well_table@wds_table
getMD(well_table)
getTVD(well_table)

library(rNodal)

well <- "
    MD  TVD
    0       0
    600     600
    1005    1000
    4075    4000
    7700    7500
    9275    9000
"

# read string text to dataframe
well_table <- WellDeviationSurvey(well, reference = "horizontal")
well_table@wds_table

# in Prosper the angle is measured againt the vertical
# we obtain the agle but the reference is set as the horizontal
library(rNodal)

md_tvd_01 <- "
MD      TVD 
0        0   
600     600
1005     1000
4075     4000
7700     7500
9275     9000
"

deviation_survey <- set_deviation_survey(md_tvd_01)
compute_angle_deviation_survey(deviation_survey, reference = "horizontal")

# split the tubing in dx pieces
apply(angle_deviation_survey, 1, function(x) x["MD"] 
)

f0nzie/rNodal documentation built on May 6, 2019, 10:14 a.m.