A new approach to upscaling fracture network models while preserving geostatistical and geomechanical characteristics

A new approach to upscaling two-dimensional fracture network models is proposed for preserving geostatistical and geomechanical characteristics of a smaller-scale “source” fracture pattern. First, the scaling properties of an outcrop system are examined in terms of spatial organization, lengths, connectivity, and normal/shear displacements using fractal geometry and power law relations. The fracture pattern is observed to be nonfractal with the fractal dimension D ≈ 2, while its length distribution tends to follow a power law with the exponent 2 < a < 3. To introduce a realistic distribution of fracture aperture and shear displacement, a geomechanical model using the combined finite-discrete element method captures the response of a fractured rock sample with a domain size L = 2 m under in situ stresses. Next, a novel scheme accommodating discrete-time random walks in recursive self-referencing lattices is developed to nucleate and propagate fractures together with their stress- and scale-dependent attributes into larger domains of up to 54 m × 54 m. The advantages of this approach include preserving the nonplanarity of natural cracks, capturing the existence of long fractures, retaining the realism of variable apertures, and respecting the stress dependency of displacement-length correlations. Hydraulic behavior of multiscale growth realizations is modeled by single-phase flow simulation, where distinct permeability scaling trends are observed for different geomechanical scenarios. A transition zone is identified where flow structure shifts from extremely channeled to distributed as the network scale increases. The results of this paper have implications for upscaling network characteristics for reservoir simulation

Role of natural fractures in damage evolution around tunnel excavation in fractured rocks

This paper studies the role of pre-existing fractures in the damage evolution around tunnel excavation in fractured rocks. The length distribution of natural fractures can be described by a power law model, whose exponent a defines the relative proportion of large and small fractures in the system. The larger a is, the higher proportion of small fractures is. A series of two-dimensional discrete fracture networks (DFNs) associated with different length exponent a and fracture intensity P21 is generated to represent various scenarios of distributed pre-existing fractures in the rock. The geomechanical behaviour of the fractured rock embedded with DFN geometry in response to isotropic/anisotropic in-situ stress conditions and excavation-induced perturbations is simulated using the hybrid finite-discrete element method (FEMDEM), which can capture the deformation of intact rocks, the interaction of matrix blocks, the displacement of natural fractures, and the propagation of new cracks. An excavation damaged zone (EDZ) develops around the man-made opening as a result of reactivation of pre-existing fractures and propagation of wing cracks. The simulation results show that when a is small, the system which is dominated by large fractures can remain stable after excavation given that P21 is not very high; however, intensive structurally-governed kinematic instability can occur if P21 is sufficiently high and the fracture spacing is much smaller than the tunnel size. With the increase of a, the system becomes more dominated by small fractures, and the EDZ is mainly created by the coalescence of small fractures near the tunnel boundary. The results of this study have important implications for designing stable underground openings for radioactive waste repositories as well as other engineering facilities that are intended to generate minimal damage in the host rock mass

Ground state fidelity in bond-alternative Ising chains with Dzyaloshinskii-Moriya interactions

A systematic analysis is performed for quantum phase transitions in a bond-alternative one-dimensional Ising model with a Dzyaloshinskii-Moriya (DM) interaction by using the fidelity of ground state wave functions based on the infinite matrix product states algorithm. For an antiferromagnetic phase, the fidelity per lattice site exhibits a bifurcation, which shows spontaneous symmetry breaking in the system. A critical DM interaction is inversely proportional to an alternating exchange coupling strength for a quantum phase transition. Further, a finite-entanglement scaling of von Neumann entropy with respect to truncation dimensions gives a central charge c = 0.5 at the critical point.Comment: 6 pages, 4 figure

Thermodynamic analysis of BN, AlN AND TiN Precipitation in boron-bearing steel

In this paper, the precipitation behavior of BN, AlN and TiN particles in boron-bearing steel was studied based on thermodynamic calculation. During solidification process, precipitation amount of BN has a positive relationship with boron content, while has negative relationship with temperature. The binding capacity of boron and nitrogen is greater than that of aluminum and nitrogen, BN is preferentially precipitated as boron added to steel. BN particle reduces the free nitrogen content in steel and then prevents the formation of AlN particle. Combination of titanium and nitrogen element is more precedence than of boron and nitrogen element. Formation of TiN particle precedes BN particle, and the precipitation amount of BN is significantly reduced by adding titanium element to boronbearing

Thermodynamic analysis of BN, AlN AND TiN Precipitation in boron-bearing steel

In this paper, the precipitation behavior of BN, AlN and TiN particles in boron-bearing steel was studied based on thermodynamic calculation. During solidification process, precipitation amount of BN has a positive relationship with boron content, while has negative relationship with temperature. The binding capacity of boron and nitrogen is greater than that of aluminum and nitrogen, BN is preferentially precipitated as boron added to steel. BN particle reduces the free nitrogen content in steel and then prevents the formation of AlN particle. Combination of titanium and nitrogen element is more precedence than of boron and nitrogen element. Formation of TiN particle precedes BN particle, and the precipitation amount of BN is significantly reduced by adding titanium element to boronbearing

Ground-State Fidelity and Kosterlitz-Thouless Phase Transition for Spin 1/2 Heisenberg Chain with Next-to-the-Nearest-Neighbor Interaction

The Kosterlitz-Thouless transition for the spin 1/2 Heisenberg chain with the next-to-the-nearest-neighbor interaction is investigated in the context of an infinite matrix product state algorithm, which is a generalization of the infinite time-evolving block decimation algorithm [G. Vidal, Phys. Rev. Lett. \textbf{98}, 070201 (2007)] to accommodate both the next-to-the-nearest-neighbor interaction and spontaneous dimerization. It is found that, in the critical regime, the algorithm automatically leads to infinite degenerate ground-state wave functions, due to the finiteness of the truncation dimension. This results in \textit{pseudo} symmetry spontaneous breakdown, as reflected in a bifurcation in the ground-state fidelity per lattice site. In addition, this allows to introduce a pseudo-order parameter to characterize the Kosterlitz-Thouless transition.Comment: 4 pages, 4 figure