Researcher profile

Eric M. Dunham

Eric M. Dunham contributes to research discovery and scholarly infrastructure.

ResearcherAffiliation not importedOpen to collaborate

Trust snapshot

Quick read

Trust 15 - UnverifiedVerification L1Unclaimed author
3works
0followers
4topics
4close collaborators

Actions

Decide how to stay connected

Follow researcher0

Identity and collaboration

How to connect with this researcher

Claiming links this public author record to a researcher profile and unlocks direct collaboration workflows.

Log in to claim

Direct collaboration

Open a focused conversation when the fit is right

Claim this author entity first to unlock direct invitations.

Research graph

See the researcher in context

Open full explorer

Inspect adjacent work, topics, institutions and collaborators without jumping out to a separate graph page.

Building this graph slice

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

3 published item(s)

preprint2020arXiv

Fault Valving and Pore Pressure Evolution in Simulations of Earthquake Sequences and Aseismic Slip

Fault-zone fluids control effective normal stress and fault strength. While most earthquake models assume a fixed pore fluid pressure distribution, geologists have documented fault valving behavior, that is, cyclic changes in pressure and unsteady fluid migration along faults. Here we quantify fault valving through 2-D antiplane shear simulations of earthquake sequences on a strike-slip fault with rate-and-state friction, upward Darcy flow along a permeable fault zone, and permeability evolution. Fluid overpressure develops during the interseismic period, when healing/sealing reduces fault permeability, and is released after earthquakes enhance permeability. Coupling between fluid flow, permeability and pressure evolution, and slip produces fluid-driven aseismic slip near the base of the seismogenic zone and earthquake swarms within the seismogenic zone, as ascending fluids pressurize and weaken the fault. This model might help explain observations of late interseismic fault unlocking, slow slip and creep transients, swarm seismicity, and rapid pressure/stress transmission in induced seismicity sequences.

preprint2020arXiv

Spring-Slider and Finite Element Modeling of Microseismic Events and Fault Slip during Hydraulic Fracturing

Hydraulic fracturing increases reservoir permeability by opening fractures and triggering slip on natural fractures and faults. While seismic slip of small faults or fault patches is detectable as microseismic events, the role of aseismic slip is poorly understood. From a modeling standpoint, geomechanical analysis using the Coulomb criterion can determine if faults slip but not whether slip is seismic or aseismic. Here we propose a computational methodology to predict fault slip, and whether slip is seismic or aseismic, using rate-and-state friction. To avoid computational costs associated with resolving small faults, we use the spring-slider idealization that treats faults as points. Interaction between faults is neglected. The method is applied to study fault slip from a hydraulic fracture that grows past a fault, without intersecting it. We represent the hydraulic fracture stressing using an asymptotic expansion of stresses around the tip of a tensile crack. We investigate the effect of fault length, orientation, and distance from the hydraulic fracture. For velocity-weakening faults with stiffness smaller than a critical stiffness, slip is seismic, whereas faults with stiffness greater than the critical stiffness slip aseismically. Furthermore, we compare the spring-slider idealization with a finite element analysis that resolves spatially variable slip. The spring-slider idealization provides reasonably accurate predictions of moment and even moment-rate history, especially for faults having stiffness close to or larger than the critical stiffness. Differences appear for large faults where rupture propagation is important, though differences might still be negligible for many applications.

preprint2019arXiv

Non-stiff narrow-stencil finite difference approximations of the Laplacian on curvilinear multiblock grids

The Laplacian appears in several partial differential equations used to model wave propagation. Summation-by-parts--simultaneous approximation term (SBP-SAT) finite difference methods are often used for such equations, as they combine computational efficiency with provable stability on curvilinear multiblock grids. However, the existing SBP-SAT discretization of the Laplacian quickly becomes prohibitively stiff as grid skewness increases. The stiffness stems from the SATs that impose inter-block couplings and Dirichlet boundary conditions. We resolve this issue by deriving stable SATs whose stiffness is almost insensitive to grid skewness. The new discretization thus allows for large time steps in explicit time integrators, even on very skewed grids. It also applies to the variable-coefficient generalization of the Laplacian. We demonstrate the efficacy and versatility of the new method by applying it to acoustic wave propagation problems inspired by marine seismic exploration and infrasound monitoring of volcanoes.