The 2016 Kumamoto earthquake sequence started with a M J (Japan Meteorological Agency magnitude) 6.5 event on April 14, and culminated in a M J 7.3 event on April 16. Associated with the sequence, approximately 34-km-long surface ruptures appeared along the eastern part of the Futagawa fault zone and the northernmost part of the Hinagu fault zone. We carried out an urgent field investigation soon after the earthquake to map the extent and displacement of surface ruptures with the following results. (1) The rupture zone generally consisted of a series of left-stepping en echelon arrays of discontinuous fault traces of various lengths. (2) Slip exceeding 100 cm occurred on previously unrecognized fault traces in the alluvial lowland of the Kiyama plain and on the western rim of the Aso volcano caldera. (3) Large slip with maximum dextral slip of 220 cm was measured throughout the central section of the rupture zone along the Futagawa segment, and the slip gradually decreased bilaterally on the adjoining northeastern and southwestern sections. (4) The surface rupture mostly occurred along fault traces mapped in previous active fault investigations. (5) Most of the surface ruptures were produced by the mainshock, and significant postseismic slip occurred after the mainshock.

Abstract:This paper presents a modeling approach for studying hydro-mechanical coupled processes, including fracture development, within geological formations. This is accomplished through the novel linking of two codes: TOUGH2, which is a widely used simulator of subsurface multiphase flow based on the finite volume method; and an implementation of the Rigid-BodySpring Network (RBSN) method, which provides a discrete (lattice) representation of material elasticity and fracture development. The modeling approach is facilitated by a Voronoi-based discretization technique, capable of representing discrete fracture networks. The TOUGH-RBSN simulator is intended to predict fracture evolution, as well as mass transport through permeable media, under dynamically changing hydrologic and mechanical conditions. Numerical results are compared with those of two independent studies involving hydro-mechanical coupling: (1) numerical modeling of swelling stress development in bentonite; and (2) experimental study of desiccation cracking in a mining waste. The comparisons show good agreement with respect to moisture content, stress development with changes in pore pressure, and time to crack initiation.The observed relationship between material thickness and crack patterns (e.g., mean spacing of cracks) is captured by the proposed modeling approach.

No abstract

No abstract

Lattice models provide discontinuous approximations of the displacement field over the computational domain, which facilitates the modeling of fracture and other discontinuous phenomena. By discretizing the domain with two-node elements, however, ordinary lattice models cannot simulate the Poisson effect in a local (intra-element) sense, which is problematic for some types of analyses.Furthermore, such methods are limited in the range of Poisson ratio values that can be simulated. We present a new approach to remedy such known, yet underappreciated, shortcomings of lattice models. In this approach, the Poisson effect is modeled through the introduction of fictitious stresses into a regular lattice. Capabilities of the new approach are demonstrated through compressive test simulations of homogeneous and heterogeneous materials. The simulation results are compared with theory and those of continuum finite element models. The comparisons show good agreement for arbitrary Poisson ratios (including ν ⩾ 1/3) with respect to nodal displacement, intra-element stress, and nodal stress. This form of discrete method, supplemented by the proposed fictitious measures of stress, retains the simplicity of collections of two-node elements.

No abstract

This study involves the development of the auxiliary stress approach for producing elasticallyhomogeneous lattice models of damage in geomaterials. The lattice models are based on random, three-dimensional assemblages of rigid-body-spring elements. Unlike conventional lattice or particle models, the elastic constants of a material (e.g., Young's modulus and Poisson's ratio) are represented properly in both global and local senses, without any need for calibration. The proposed approach is demonstrated and validated through analyses of homogeneous and heterogeneous systems under uniand tri-axial loading conditions. Comparisons are made with analytical solutions and finite element results. Thereafter, the model is used to simulate a series of standard laboratory tests: (a) splitcylinder tests, and (b) uniaxial compressive tests of sedimentary rocks at the Horonobe Underground Research Laboratory in Hokkaido, Japan. Model inputs are based on physical quantities measured in the experiments. The simulation results agree well with the experimental results in terms of pre-peak stress-strain/displacement responses, strength measurements, and failure patterns.

[1] The problem of solute transport through a water-saturated single fracture in a permeable rock matrix is examined using an analytical modeling approach. A closed-form analytical solution is obtained that accounts for transverse and longitudinal advective transport in the fracture and matrix and transverse diffusion in the matrix. The solution also accounts for both diffusive and advective solute exchange between the fracture and matrix and a general solute source position in either the fracture or matrix. The novel features are the incorporation of advective transport in the matrix and a general source position into a closed-form solution for the solute-transport problem. Examples of the solution behavior are presented, which demonstrate the effects of matrix advection in combination with advection along the fracture, transverse diffusion in the matrix for solute release in the fracture and matrix. A semianalytical solution in the form of a superposition integral is also derived that includes these transport features, plus independent levels of longitudinal diffusion and dispersion in the matrix and fracture, respectively. Examples are presented that include advective transport in the fracture and matrix, longitudinal and transverse diffusion in the matrix, longitudinal dispersion in the fracture, as well as solute release from the fracture and matrix. An approximate criterion is proposed to evaluate the significance of longitudinal diffusion and dispersion relative to longitudinal spreading caused by fracturematrix interaction.