A Semi-Analytical Geomechanical Approach for Forecasting Production Performance in Multifractured Composite Systems

Stress-dependence of reservoir matrix and fractures can strongly affect the performance of multifractured horizontal wells (MFHWs) completed in unconventional hydrocarbon reservoirs. In order to model fluid flow in unconventional reservoirs exhibiting this stress-dependence, most traditional reservoir flow simulators, and many simulators described in published work, use conventional reservoir fluid flow model formulations. These formulations typically neglect the influence of the rate of change of volumetric strain of the reservoir matrix and fractures, even though reservoir stress and pressure change significantly during the course of production. As a result, the effect of matrix and fracture deformation on production is neglected, which can lead to errors in predicting production performance in most stress-sensitive reservoirs. To address this problem, some studies have proposed the use of porosity and transmissibility multipliers to model stress-sensitive reservoirs. However, in order to apply this approach, multipliers must be estimated from laboratory experiments, or used as a history-match parameter, possibly resulting in large errors in well performance predictions. Alternatively, fully-coupled, fully numerical geomechanical simulation can be performed, but these methods are computationally costly, and models are difficult to setup. This paper presents a new fully-coupled, two-way analytical modeling approach that can be used to simulate fluid flow in stress-sensitive unconventional reservoirs produced through MFHWs. The model couples poroelastic geomechanics theory with fluid flow formulations. The two-way coupled fluid flow-geomechanical analytical model is applied simultaneously to both the matrix and fracture regions. In the proposed algorithm, a porosity-compressibility coupling parameter for the two physical models is setup to update the stress- and pressure-dependent fracture/matrix properties iteratively, which are later used as input data for the fracture-matrix reservoir fluid flow model at each iteration step. The analytical approach developed for the fully-coupled, two-way analytical model, using the enhanced fracture region conceptual model, is validated by comparing the results with numerical simulation. Predictions using the fully-coupled enhanced fracture region model are then compared with the same enhanced fracture region model but with the conventional pressure-dependent modeling approach implemented. A sensitivity study performed by comparing the new fully-coupled model predictions with and without geomechanics effects accounted for reveals that, without geomechanics effects, production performance in stress-sensitive reservoirs might be overestimated. The study also demonstrates that use of the conventional stress-dependent modeling approach may cause production performance to be underestimated. Therefore, the proposed fully-coupled, two-way analytical model can be useful for practical engineering purposes.

Damage Evaluation of Concrete under Uniaxial Compression Based on the Stress Dependence of AE Elastic Wave Velocity Combined with DIC Technology

This study presented evaluation of a concrete damage process by the acoustic emission (AE) technique under uniaxial multi-step compressive loading procedure combined with digital image correlation (DIC). The results showed that AE elastic wave velocity had good stress dependence in the damage process of concrete specimens with different sizes (cube, prism) and coarse aggregate characteristics (volume fraction, maximum size), and the effects of specimen sizes and coarse aggregate characteristics on the stress dependence can be nearly neglected. The standard deviation of 32 AE elastic wave velocities was used as the criterion to evaluate the relative stress ratio of concrete under different damage states, and the damage process of concrete was divided into three damage stages according to this criterion. When the standard deviation is below 70, in the range of 70 to 1700, and greater than 1700, the concrete damage process is defined as steady damage process, accelerated damage process and buckling damage process, respectively. The accuracy of the presented evaluation methodology was demonstrated by comparative results with digital image correlation. The results indicate that the standard deviation of AE elastic wave velocities can potentially serve as a reliable, convenient, and non-destructive evaluation criterion of concrete damage state under uniaxial compressive loading.

Flow and Yield Characteristics of Yield Stress Fluids Using Hysteresis Loop Test Below Slip Yield Point

The flow characteristics of angel O/W emulsion, which is a yield stress fluid, was investigated. The hysteresis loop test was conducted for the strain below the slip yield point, and the single relaxation Maxwell model was used to fit the experimental data. Using these methods, the shear-rate dependence, stress dependence, and time dependence of the viscoelastic properties of the sample were evaluated in the region below the slip yield point. The shear-rate dependence induced by the stress-ramp rate and the stress dependence from the maximum applied stress influence the viscoelastic characteristics below the slip yield point in terms of the flow history. However, the time dependence of the viscoelastic characteristics could not be confirmed for any creep time. The yield stress measured in the stress-ramp test increases with the stress-ramp rate owing to the contribution of the viscous strain from the flow history.