Organic shales can split into laminae or thin layers due to their fissility and are usually characterized as transverse isotropic media. Understanding their elastic and anisotropic properties plays a significant role in geophysical modeling and imaging, reservoir geomechanics, and reservoir characterization. One can obtain their velocity anisotropy through laboratory and field measurements or rock-physics models. However, the measurements are scarce, and theoretical modeling of Thomsen’s anisotropy parameters generally requires extra inputs, which are sometimes difficult to obtain. We have compiled ultrasonic data from 159 self-sourced reservoir shale samples from the literature and our own measurements. Simple exponential models can capture the trends of the P- and S-wave anisotropy parameters ε and γ decreasing with the increasing vertical P- and S-wave velocities, with R2 being above 0.84. The error in estimated anisotropy is small for high-velocity (VP0>4 km/s) shales but relatively large for low-velocity (VP0<4 km/s) shales. We analyze the influence of multiple geologic factors on the proposed anisotropy-velocity relationships. We observe that velocity anisotropy generally increases with clay content. In addition, mature shales (0.6 < Ro (%) < 1.4) generally have stronger anisotropy strength spreading in a broader range than overmature shales (Ro (%) > 1.4). Overall, no single parameter dominates the source of velocity anisotropy, which is jointly affected by the coupled factors of clay content, organic matter, porosity, maturity level, and geologic history. The empirical model has the potential to offer a fast and straightforward method of predicting P- and S-wave anisotropy strength from vertical P- and S-wave velocities.

Understanding the elastic and attenuation signatures of shales is of considerable interest for unconventional reservoir characterization and sealing capacity evaluation for CO2 sequestration and nuclear waste disposal. We have conducted laboratory measurements on seven shale samples at seismic frequencies (2–100 Hz) to study the effects of clay-bound water (CBW) on their wave dispersion and attenuation signatures. With nuclear magnetic resonance and a helium porosimeter, the volume of CBW in the shale samples is quantified. The forced-oscillation measurement reveals that Young’s modulus exhibits a continuous dispersion trend from 2 to 100 Hz. The extensional attenuation (QE−1) shows a weak frequency and pressure dependence on effective pressure ranging from 5 to 35 MPa. The magnitude of extensional attenuation shows a positive correlation with CBW, with an R2 value of 0.89. It is found that 4% of CBW in the rock frame causes approximately a 5% modulus increase from 2 to 100 Hz. We adopt a constant Q model for assigning frequency-dependent bulk and shear moduli to the CBW in the rock-physics modeling, which can fit the experimental data of modulus dispersion and attenuation well, indicating that the bulk and shear moduli of CBW in shales might behave viscoelastically.

The static bulk modulus of unconsolidated sands is essential for predicting the in situ effective pressure to reduce drilling risks in deepwater reservoirs; however, the dynamic bulk modulus is often more broadly available from seismic or well logging data. Therefore, it is tempting to investigate the relationship between static and dynamic bulk moduli. We perform a series of ultrasonic velocity measurements on 21 deepwater reservoir sands from the Gulf of Mexico (GoM) to study the static and dynamic bulk moduli simultaneously. Both room-dry and brine-saturated ultrasonic velocities are measured under hydrostatic stress conditions to derive the dynamic bulk moduli. Under brine-saturated conditions, if the pore pressure is kept constant, the pore volume change with the confining pressure can be monitored accurately by a digital pump, which is subsequently used to estimate the static bulk modulus. The experimental results suggest that both the static and dynamic bulk moduli decrease upon pressure unloading. The pressure-dependent bulk moduli are modeled using the Hertz-Mindlin contact theory at the critical porosity and combined with the modified Hashin-Shtrikman lower bound for other porosities. The results suggest that the theoretical estimates can serve as the lower bound for the dynamic bulk modulus and the upper bound for the static bulk modulus. Under room-dry conditions, the static-to-dynamic modulus ratio decreases from a value approaching 0.8 to approximately 0.25 with decreasing differential pressure. Moreover, the effects of brine saturation on the relationship between the static and dynamic moduli are investigated using Gassmann’s equation. The brine saturation substantially reduces the difference between the static and dynamic bulk moduli, making the static-to-dynamic modulus ratio approach unity, which may be relevant to the in situ reservoir rock properties.

Diagnosing fractures under compression is of great importance in optimizing hydraulic fracturing stimulation strategies for unconventional reservoirs. However, a lot of information, such as fracture morphology and fracture complexity, is far from being fully excavated in the laboratory limited by the immature fracture identification techniques. In the current study, we propose a set of methods to analyze the fracture complexity of cylindrical cores after triaxial fracturing. Rock failure under conventional compression tests is real-time controlled by monitoring the stress–strain evolutions to ensure that the cores remain cylindrical after failure. The lateral surface of the core cylinders is scanned with a 2D optical scanner to extract the fracture parameters, surface fracture rate, and inclination dispersion, which are normalized and averaged to derive the fracture complexity. After analyzing the data for 24 shale gas reservoir cores from the Sichuan Basin, the fractal dimension of fracture images shows a good linear correlation with the surface fracture rate but has no correlation with the dip dispersion. The calculated fracture complexity has nearly no relationship with the E-v–based brittleness index but demonstrates a positive correlation with the mineral content–based brittleness index. Moreover, the fracture complexity is associated with the core mineralogical compositions. The fracture complexity is positively correlated with the content of quartz, calcite, and dolomite and negatively correlated with the content of clay minerals and has no obvious relationship with the content of feldspar. The proposed method provides an experimental basis for the evaluation of fracturability of unconventional oil and gas reservoirs.

Tensile mechanical behavior of T700 fiber bundles have been studied in the strain rate range from 0.001 l/s to 1000 l/s and complete stress strain curves were obtained. Results show that both ultimate strength and failure strain are strain rate insensitive, and T700 carbon fiber can be regarded as strain rate insensitive materials. On basis of the fiber bundles model and the statistic theory of fiber strength, a statistical constitutive model based on Weibull distribution function has been developed to describe mechanical behavior of T700 fiber. And a method for determine the statistic parameters of fibers by tensile tests of fiber bundles is established, too. The simulated stress strain curves from the model are in good agreement with the test data. Simulated results show that the strength of T700 fiber can be described by single Weibull distribution function.

Determining the crack initiation stress (Ci) for unconventional shale rocks is of critical importance in describing the entire failure process of unconventional shale reservoirs. We propose a new method to identify Ci values based on triaxial failure tests on four organic shale samples, attempting to improve the shortcomings of other methods. The new method is based on the relationship between crack development and strain energy evolution (SEE). Additionally, the proposed SEE method is compared with three widely used methods, including crack volumetric strain (CVS), moving point regression (MPR) and the lateral strain response (LSR), intending to examine the performance of different methods. The contrastive results indicate that the LSR method cannot determine Ci when the rock ruptures without volumetric dilatancy, which frequently occurs in the compression process of organic shales. Ci values obtained using the SEE method are consistent with those from the CVS and MPR methods. However, the proposed SEE method with a solid physical basis is more objective and stable than the CVS and MPR methods. The proposed method, from one aspect, compensates for the shortcomings of other methods when facing different failure modes in organic shales. From the other aspect, it provides a way to precisely determine Ci values for applications in wellbore stability evaluation and hydraulic fracturing design.