Abstract Field observations and unmanned aerial vehicle surveys from Caldera Taburiente (La Palma, Canary Islands, Spain) show that pre‐existing dykes can capture and re‐direct younger ones to form multiple dyke composites. Chill margins suggest that the older dykes were solidified and cooled when this occurred. In one multiple dyke example, an 40 Ar/ 39 Ar age difference of 200 kyr was determined between co‐located dykes. Petrography and geomechanical measurements (ultrasonic pulse and Brazilian disc tests) show that a microscopic preferred alignment of plagioclase laths and sheet‐like structures formed by non‐randomly distributed vesicles give the solidified dykes anisotropic elastic moduli and fracture toughness. We hypothesize that this anisotropy led to the development of margin‐parallel joints within the dykes, during subsequent volcanic loading. Finite element models also suggest that the elastic contrast between solidified dykes and their host rock elevated and re‐oriented the stresses that governed subsequent dyke propagation. Thus, the margin‐parallel joints, combined with local concentration and rotation of stresses, favored the deflection of subsequent magma‐filled fractures by up to 60° to form the multiple dykes. At the edifice scale, the capture and deflection of active intrusions by older ones could change the organization of volcanic magma plumbing systems and cause unexpected propagation paths relative to the regional stress. We suggest that reactivation of older dykes by this mechanism gives the volcanic edifice a structural memory of past stress states, potentially encouraging the re‐use of older vents and deflecting intrusions along volcanic rift zones or toward shallow magma reservoirs.
Abstract Slip and slip zone development during constant‐rate fluid injection into a permeable fault in an infinite, impermeable elastic rock is studied here numerically using a hydromechanical fracture model. The interplay of slip zone growth and fault dilatancy is assumed to be affected by linear slip weakening of frictional strength and slip‐induced dilation. The discretization error is minimized by using element sizes based on a mesh sensitivity study to obtain accurate slipping and pressurized lengths. Comparisons with published results for large and zero fault permeability cases demonstrate the validity of the numerical results. Dimensional analysis is used to identify five parameters that control slip development, including the slip‐induced dilatancy factor, initial fault permeability, the product of fluid viscosity and injection rate, the background shear to normal stress ratio, and the residual coefficient of friction. Conditions leading a quasi‐static stable slip to either continue as stable aseismic slip or accelerate are discussed. Acceleration of slip occurs on less permeable faults characterized by a hybrid tensile‐shear fracture subject to higher pressure, while prolonged aseismic slip occurs along a more permeable fault subject to lower pressure. During long‐duration aseismic slip the pressure varies slightly around the value that equalizes the fault shear strength and the background shear stress. A permeability perturbation in the form of a sinusoidal fault aperture distribution does not change the tendency of the slip to be slow and stable. Like preexisting permeability, slip‐induced dilation responses play a role in limiting pressure level and extending the slow slip period.
Abstract Understanding the formation mechanisms of complex fracture networks is vitally important for hydraulic fracturing operations in shale formation. For this purpose, a hydraulic fracturing experiment under a core-plunger scale is conducted to investigate the impact of the bedding plane angle, borehole size, and injection rate on fracture initiation behaviors of laminated shale rock. The results on rock properties demonstrate that the anisotropic characteristics of shale rock are reflected not only in elastic modulus but also in tensile strength. The results of fracturing experiments show that the bedding plane dip angle and borehole size have significant effects on fracture initiation behaviors, in that fracture initiation pressure (FIP) decreases with the increase of those two factors. The impact of injection rate, by contrast, has no obvious variety regulation. The above data is further used to validate our previously proposed fully anisotropic FIP model, which shows better agreement with experimental results than those using other models under various parameter combinations. Finally, a postfracturing analysis is performed to identify the fracture growth patterns and the microstructures on the fracture surfaces. The results show that the hydraulic fractures (HFs) always grow along mechanically favorable directions, and the potential interaction between HFs and bedding planes mainly manifests as fracture arrest. Meanwhile, the roughness of fracture surfaces is physically different from each other, which in turn results in the difficulties of fluid flow and proppant migration. The findings of this study can help for a better understanding of the fracture initiation behavior of laminated shale rock and the corresponding fracture morphology.
Summary Placing multiple hydraulic fractures at intervals along horizontal wells has proved to be a highly effective method for stimulation. However, the mechanical interaction between a growing hydraulic fracture and one or more previous hydraulic fractures can affect the fracture geometry such that the final fracture array is suboptimal for stimulation. If the fracture-array geometry is idealized as a set of regular and planar fractures, history matching and production forecasting may be inaccurate. During the treatments, the fractures can curve toward or away from one another, potentially intersecting one another. A detailed parametric study of this phenomenon using a coupled 2D numerical fracturing simulator shows that the curving is associated with a combination of opening and sliding along the previously placed hydraulic fracture, as well as the previous fracture's disturbance of the local stress field because of its propped width. Dimensional analysis and scaling techniques are used to identify the key parameters that are associated with suppression of each mechanism that can lead to hydraulic-fracture curving. The analysis, which is in agreement with available data, results in a clarification of the conditions under which attractive and repulsive curving are expected, as well as the conditions under which curving is expected to be negligible or even completely suppressed. This last case of planar hydraulic-fracture growth is of practical importance and will usually be considered desirable. We present a straightforward method for determining whether planar fracture growth is expected that additionally gives insight into how design parameters can be modified to promote planar hydraulic-fracture growth.
Magmatic sulphides largely control the behaviour of chalcophile elements in evolving magmas and also play critical roles in ore-forming processes at convergent plate margins. Magmas at the back-arc spreading centre of Mariana Trough are less affected by slab-derived fluids than Mariana arc basalts. However, whether the evolution of chalcophile elements in Mariana Trough magmas resembles or differs from those in arc magmas is poorly known. Based on the detailed petrological and geochemical studies, we present high-precision analyses of whole-rock chalcophile element (S, Cu, Ag, Re, Pd, Pt, Ru, Ir) contents for a suit of basaltic lavas (n = 15) from the central Mariana Trough. The occurrence of sulphides in high-Mg olivine crystals (Fo80) indicates that at least one episode of sulphide saturation occurred during early-stage magmatic evolution. The sulphur contents of most basalts were not significantly affected by degassing processes due to low oxidation state and high hydrostatic pressures. Instead, progressive fractionation of liquid sulphide with magmatic evolution is the dominant process, as reflected by sulphide petrography, the decrease in Cu, Ag, Pd, and Pt contents, and concomitant increase in Cu/Pd ratios (22,500–708,500), as well as nearly constant Cu/Ag (mean of 3137 ± 848, 1s) with decreasing MgO. We attribute early sulphide saturation of Mariana Trough basalts to their generation and evolution under relatively reducing conditions, with the limited effect of slab components. Combined with data from other island arc and back-arc lavas, we further show that the extent of adding slab-derived components, regardless of specific arc and back-arc tectonic settings, are determinants of the sulphide evolution and chalcophile element behaviour in magmatic systems above intra-oceanic subduction zones. Early magmatic sulphide saturation at central Mariana Trough would result in barren magmatic systems that are not suitable for the large-scale Cu-Au mineralization at the seafloor.
'/%3e%3cuse xlink:href='%23a' transform='rotate(23.036 .093 25.536)'/%3e%3cuse xlink:href='%23a' transform='rotate(45.87 1.273 16.18)'/%3e%3cuse xlink:href='%23a' transform='rotate(69.945 .996 12.078)'/%3e%3cuse xlink:href='%23a' transform='rotate(20.66 -19.689 31.932)'/%3e%3c/svg%3e)
'/%3e%3c/defs%3e%3cpath fill='%23012169' d='M0 0h10080v5040H0z'/%3e%3cpath stroke='%23fff' stroke-width='.6' d='m0 0 6 3m0-3L0 3' clip-path='url(%23b)' transform='scale(840)'/%3e%3cpath stroke='%23e4002b' stroke-width='.4' d='m0 0 6 3m0-3L0 3' clip-path='url(%23c)' transform='scale(840)'/%3e%3cpath stroke='%23fff' stroke-width='840' d='M2520 0v2520M0 1260h5040'/%3e%3cpath stroke='%23e4002b' stroke-width='504' d='M2520 0v2520M0 1260h5040'/%3e%3cg fill='%23fff'%3e%3cuse xlink:href='%23d' x='2520' y='3780'/%3e%3cuse xlink:href='%23a' x='7560' y='4200'/%3e%3cuse xlink:href='%23a' x='6300' y='2205'/%3e%3cuse xlink:href='%23a' x='7560' y='840'/%3e%3cuse xlink:href='%23a' x='8680' y='1869'/%3e%3cuse xlink:href='%23e' x='8064' y='2730'/%3e%3c/g%3e%3c/svg%3e)