rc76_edge_disloc_poro: Rice and Cleary solutions

In abarbour/deform: Crustal Deformation Modeling Tools in R

Description

Rice and Cleary solutions

Usage

rc76_edge_disloc_poro(
  Times,
  slip,
  radial_dist,
  flt_angle = 30,
  diffusiv = 1,
  nuu = 0.33,
  nu = 0.25,
  B = 0.6,
  mu = 3e+10,
  ...
)

Arguments

Times	numeric; times to calculate at depth [s]
slip	numeric; relative slip along the edge dislocation [m]
radial_dist	numeric; distance from xxx to xxx [m]
flt_angle	numeric; angle counter clockwise from fault strike
diffusiv	numeric; in plane hydraulic diffusivity [m^2/s]
nuu	numeric; undrained Poisson's ratio [-]
nu	numeric; drained Poisson's ratio [-]
B	numeric; Skempton's coefficient [-]
mu	numeric; elastic shear modulus [Pa]
...	additional arguments

Details

rc76_edge_disloc_poro is the 2D solution to an edge dislocation (e.g., fault) embedded in a poroelastic medium

Value

tibble

Examples

# Pore pressure transient associated with near-field fault slip

ti <- c(-100, -1, 0, 10**seq(-4, log10(6598), length.out=301))
Mw <- 4.33
M0 <- Mw2M0(Mw) #3.935501e+15 N*m for Mw 4.33
bx <- 10^(-3.22 + 0.69 * Mw) #0.5
rd <- 250
mu <- 2e8
Diffu <- 60
rc <- rc76_edge_disloc_poro(ti, bx, rd, diffusiv=Diffu, mu=mu)
nrc <- rc76_edge_disloc_poro(ti, bx, -rd, diffusiv=Diffu, mu=mu)
all.equal(rc$pp, -nrc$pp)

## Not run: 
plot(pp / 1e3 + 5 ~ I(Times / 3600), rc, type='l', lwd=2, ylab='kPa', 
     ylim=c(-8.8,36.9), yaxs='i', xaxs='i')
axis(4)
abline(v=1*60/3600, lty=3)
rect(0,0,100,100)
SrcArea <- M0 / (mu * bx)
SrcRad <- sqrt(SrcArea / pi)
print(c(Area=SrcArea, Rad=SrcRad))

## End(Not run)

abarbour/deform documentation built on Feb. 15, 2022, 6:24 p.m.