kitagawa (version 3.1.2)

owrsp-methods: Generic methods for objects with class 'owrsp'.

Description

An object with class 'owrsp' is a list containing the response information, and the mechanical, hydraulic, and material properties used to generate the response for an open well.

Usage

# S3 method for owrsp
as.data.frame(x, ...)
data.frame.owrsp(x, ...)
# S3 method for owrsp
print(x, n = 3, ...)
# S3 method for owrsp
summary(object, ...)
# S3 method for summary.owrsp
print(x, ...)
# S3 method for owrsp
lines(x, series = c("amp", "phs"), ...)
# S3 method for owrsp
points(x, series = c("amp", "phs"), pch = "+", ...)
# S3 method for owrsp
plot(
  x,
  xlims = c(-3, 1),
  ylims = list(amp = NULL, phs = 185 * c(-1, 1)),
  logamp = TRUE,
  ...
)

Arguments

x

'owrsp' object

...

optional arguments

n

numeric; the number of head and tail to print

object

'owrsp' object

series

character; the series to plot (amplitude or phase)

pch

point character, as in par

xlims

limits for x-axis (applies to both amp and phs frames)

ylims

optional list of limits for y-axis (i.e., list(amp=c(..),phs=c(...)))

logamp

logical; should the amplitude be in log10 units

Author

A. J. Barbour <[e-mail]>

Details

The response information is a matrix with frequency, complex response [\(\omega\), \(Z_\alpha (\omega)\)] where the units of \(\omega\) will be as they were input. The amplitude of \(Z\) is in meters per strain, and the phase is in radians.

See Also

open_well_response

kitagawa-package

Examples

S. <- 1e-5  	# Storativity [nondimensional]
T. <- 1e-4		# Transmissivity [m**2 / s]
frq <- 1/(1:200)
# Defaults to the Rojstaczer formulation
W <- open_well_response(frq, T. = T., S. = S., Rs. = 0.12, freq.units="Hz")
# (warning message about missing 'z')
W <- open_well_response(frq, T. = T., S. = S., Rs. = 0.12, freq.units="Hz", z=1)
str(W)
print(W)
print(summary(W))
plot(rnorm(10), xlim=c(-1,11), ylim=c(-2,2))
lines(W)
lines(W, "phs", col="red")
points(W)
points(W, "phs")
#
Wdf <- as.data.frame(W)
plot(Mod(wellresp) ~ omega, Wdf) # amplitude
plot(Arg(wellresp) ~ omega, Wdf) # phase
plot(W)
# change limits:
plot(W, xlims=c(-4,0), ylims=list(amp=c(-7,-3), phs=185*c(-1,1)))