ALBERT

Description: Summary Offshore oil and gas has effectively alleviated the global shortage of oil and gas resources, and drilling operations are becoming increasingly frequent. However, the cuttings discharged during surface drilling are transported and deposited to form cuttings piles, which pose a serious threat to the marine ecological environment. In this study, we consider the randomness and uncertainty of cuttings movement to divide the transport process into parabola and collision motion between the moving particles and slope particles after falling on the slope surface of cuttings piles. Through specific analysis of the stress state of a single particle in the transport process and changes in momentum distribution of the particle swarm, the evolution model of the morphological distribution of cuttings piles and the nearby flow field is established. This model can quantitatively analyze the evolution law of the morphological distribution of cuttings piles under the action of ocean current and the disturbance law of the flow field near the cuttings piles caused by the invasion of cuttings particles. Comparing the measured data at an offshore drilling field and prediction results of the model of Sun et al. (2020), the relative error of the model amounts to less than 15%, which demonstrates its rationality. The simulation results show that the morphological distribution of cuttings piles and the nearby flow field change significantly under the action of ocean current, and the intensity of evolution is related to the current velocity and cuttings size, which is of great significance for the quantitative analysis of the evolution of cuttings piles under the action of ocean currents and accurate prediction of their morphological distribution.

Description: Summary We deploy an adaptive observer recently developed for general hyperbolic partial differential equation (PDE) systems to detect and diagnose various drilling incidents. The well is modeled by a distributed PDE which, contrary to lumped models, preserves fundamental properties of well-flow dynamics enabling faster and more accurate incident detection and estimation. Wired drillpipe technology with pressure sensors is needed to locate and isolate incidents. Four realistic simulation case studies demonstrating various properties of the observer are presented. Drilling incidents treated in the simulation case studies include packoff in the annulus, formation in-flows, loss of circulation, and various combinations of these. Although simulation results show that the developed observers successfully estimate properties of the incidents that they are tailored for, they do not constitute an incident-detection system for drilling. However, they provide a part of the data on which an overall incident-detection system can rely.

Description: Summary X-ray imaging of porous media has revolutionized the interpretation of various microscale phenomena in subsurface systems. The volumetric images acquired from this technology, known as digital rocks (DR), make it a suitable candidate for machine learning and computer-vision applications. The current routine DR frameworks involving image processing and modeling are susceptible to user bias and expensive computation requirements, especially for large domains. In comparison, the inference with trained machine-learning models can be significantly cheaper and computationally faster. Here we apply two popular convolutional neural network (ConvNet) architectures [residual network (ResNet) and ResNext] to learn the geometry of the pore space in 3D porous media images in a supervised learning scheme for flow-based characterization. The virtual permeability of the images to train the models is computed through a numerical simulation solver. Multiple ResNet variants are then trained to predict the continuous permeability value (regression). Our findings demonstrate the suitability of such networks to characterize volume images without having to resort to further ad-hoc and complex model adjustments. We show that training with richer representation of pore space improves the overall performance. We also compare the performance of the models statistically based on multiple metrics to assess the accuracy of the regression. The model inference of permeability from an unseen sandstone sample is executed on a standard workstation in less than 120 ms/sample and shows a score of 0.87 using explained variance score (EVS) metric, a mean absolute error (MAE) of 0.040 darcies, and 18.9% relative error in predicting the value of permeability compared to values acquired through simulation. Similar metrics are obtained when training with carbonate rock images. The training wall time and hyperparameters setting of the model are discussed. The findings of this study demonstrate the significant potential of machine learning for accurate DR analysis and rock typing while leveraging automation and scalability.

Description: Summary Undetected gas kicks are at the root of many disastrous accidents in the oil industry. A major factor in these tragic outcomes is the suddenness with which oil and gas can be ejected from the top of the riser, which gives operators virtually no time for an adequate response. This aspect of the phenomenon is quantified by providing a physical explanation of its origin and suggesting simple, if approximate, guidelines to estimate its severity. Then a relatively straightforward and potentially practical way to detect the presence of dangerous gas volumes in a riser is described. The basic idea is to measure the pressure difference between sensors spaced along the riser. When gas occupies the space between sensors, the pressure difference undergoes large changes chiefly caused by the decreased hydrostatic head. It is shown that this difference can be enhanced by adequate signal processing and is robust in the presence of noise. The suggested detection method is supported by a set of laboratory experiments.

Description: Summary Progressing cavity pump (PCP) is the essential booster equipment in oil–gas mixing delivery. Changes in relevant parameters in PCP operations directly affect the working performance and service life of the pump. On the basis of computational fluid dynamics (CFD) in this study, we apply dynamic grid technology to establish a 3D flow field numerical calculation model for the CQ11-2.4J PCP, which is used in the field of the Hounan Operation Area in Changqing oil field, China. The effects of several operating parameters, such as oil viscosity, pump rotation speed, differential pump pressure, and void fraction of oil, on the pressure and the velocity distribution of the PCP flow field are examined. Various performance parameters in the transport of the oil–gastwo-phase mixture are used in the analysis, including volumetric flow rate, slippage, shaft power, volumetric efficiency, and system efficiency. The results show that the pressure and speed distribution in the pump chamber of the PCP is relatively homogenous under different working conditions, whereas the pressure and speed exhibited sharp changes at the stator and rotor sealing line and adjacent areas in the pump chamber. Increasing the viscosity of the oil and the speed of the rotor can effectively improve the flow characteristics of the PCP, but extremely high pump rotation speed would cause a decline in system efficiency. Increasing the differential pressure and the void fraction of oil would result in a decrease in the volumetric flow rate and efficiency of the PCP. Considering the variation law of the PCP's performance parameters, the optimal interval for each operating parameter of the PCP is as follows: Oil viscosity at 50–100 mPa·s, pump rotation speed at 200–300 rev/min, differential pressure at 0.2–0.3 MPa, and the void fraction of oil not more than 50%. This research can provide technical support for the optimization of the working conditions of the PCP on site.

Description: Summary The litmus test for downhole multiphase flowmeters is to compare the measured phase flow rates with the rates from a test separator or other surface measurement systems. In most cases, the composition of the measurand is required for flowmeters. This is typically obtained from bottomhole fluid samples. Extracting and analyzing fluid samples is an expensive process mostly done at the initial stages of field development. In some cases, the composition may be old or unavailable, leading to subpar flowmeter performance compared to surface systems. In this work, it is shown that when the data from a surface system such as a test separator are used in conjunction with the mixture sound speed measured downhole, it is possible to optimize a downhole multiphase flowmeter system without obtaining new fluid samples. The optimization process is independent of the downhole measurement device because the required flow-velocity and sound-speed measurements may be obtained from separate devices. For example, the fluid bulk velocity and mixture sound speed can be measured by a local measurement device and a distributed acoustic sensing (DAS) system, respectively. The main challenge in a flow-velocity/sound-speed measurement system is determining individual phase sound speeds so that the mixture phase fraction can be correctly determined using Wood’s mixture sound speed model. The phase fraction from the separator tests can be used as the target value to optimize the performance of the system. The system has two operation modes. In optimization mode, the individual phase sound speeds are calculated backward using the predicted phase fractions from surface measurements. Pressure and temperature variations at measurement locations, as well as pipe compliance effects, are accounted for during the process. After the adjustment of individual phase sound speeds, steady-state operation mode takes over, and a forward calculation is implemented using the same model. The final phase fraction agrees well with the actual value and can be improved further with an iterative approach. This novel method is demonstrated in a North Sea case history. A downhole optical flowmeter in a North Sea field measured mixture velocity and sound speed. Well-test results indicated that water cut from the flowmeter was underreported and phase flow rates did not match test-separator rates. Instead of halting production and going through a fluid sample analysis cycle, the test-separator water cut was used as the target value to optimize oil phase sound speed using Wood’s model in the optimization mode. The difference between the initial and optimized oil sound speeds was extrapolated to other pressure and temperature conditions, and steady-state operation mode showed that separator tests and flowmeter measurements closely matched. Subsequent flowmeter and test-separator data confirmed excellent agreement. Using surface measurements and downhole mixture sound speed to optimize phase flow rates is a novel method that has not been previously demonstrated. This method is independent of device type, is broadly applicable, and improves the understanding of multiphase flow measurement.

Description: Summary Poor hole cleaning is a major concern in coiled-tubing drilling (CTD), and it is often associated with long nonproductive time that contributes significantly to the operational cost. In this study, a transient solids transport model is developed based on transport equations of phases in the flow to predict the evolution of solids conveyed in the wellbore. The developed model is able to provide forecasts of the distribution of cuttings along the annulus, which can be important information for deciding to improve solids removal. Based on the concept of a two-layersteady-state model, a 1D time-dependent model is developed using two layers: a lower layer of solids bed and an upper layer of a solid-liquid mixture with the mechanisms of solids deposition and solids entrainment taken into account. The model is discretized by using a finite volume scheme and then solved by employing a semi-implicit numerical solution. The model’s hyperparameters, such as a deposition factor and an entrainment factor, are calibrated with experimental data conducted by the use of the large indoor flow loop (LIFL) to achieve a better match. The model is combined with a 2D cross-sectional model to handle the effect of pipe eccentricity and bed presence. Predictions from the model agree well with the experimental data acquired by using an oil-based mud for the majority of the cases.

Description: Summary Annulus pressure buildup (APB) problems in shale gas wells seriously affected on the safety and efficient exploitation of shale gas all around the world. The sealing failure of the cement sheath on interfaces caused by periodically changed fluid pressure in casing during hydraulic fracturing is treated as the main reason for APB in shale gas wells. Many methods are put forward to solve the APB problem in the field, and fortunately, the preapplied annulus backpressure (PABP) method shows an excellent utility. In this paper, an analytical model is established to explain the mechanism of the PABP method increasing the sealing ability of the cement sheath. The residual strain of the cement sheath and radial stress on interfaces are considered to analyze the factors that affect the effectiveness of the PABP method. In addition, based on the field data, an experimental device is established to test the validity of the PABP method and to certify the accuracy of the analytical model established in this paper. The analytical results show that the thickness of the casing has little effect on radial stress on interfaces. The outer diameter of the casing and the thickness of the cement sheath can temperately affect the radial stress. The elastic modulus of the cement sheath and the formation rock can significantly affect the radial stress. The higher elastic modulus of the cement sheath can dramatically increase the radial stress on interfaces. On the contrary, the higher elastic modulus of formation rock will induce smaller radial stress on the interfaces. In the field, the number of newly added shale gas wells with APB problems has dramatically decreased by using the PABP method. The work in this paper can be significantly useful for researchers and engineers to reduce the APB in shale gas wells.

Description: Summary Caprocks play a crucial role in the geological storage of carbon dioxide (CO2) by preventing its escape and thus trapping it into underlying sequestering reservoirs. An evaluation of interaction-induced alteration of the hydromechanical properties of caprocks is essential to better assess the leaking risk and injection-induced rock instability, thus ensuring a long-term viability of geological CO2 storage. We study the changes in minerals, nanopores, elastic velocities, and mechanical responses of a carbonate caprock caused by rock–water/brine–CO2 interaction (CO2 pressure: ≈12 MPa; 50°C). Before the interaction, the total and accessible porosities are 1.6 and 0.6%, respectively, as characterized by the small-angle neutron scattering (SANS) technique. SANS results show that the total porosity of the carbonate caprock increases, apparently because of rock–brine–CO2 interaction, and the increasing rate rises as brine concentration increases (2.2% for 0 M NaCl, 2.6% for 1 M NaCl, and 2.7% for 4 M NaCl). The increase in total porosity is due to the dissolution of calcite, which tends to enlarge accessible pores (by 0.8 to 1.2%) while slightly decreasing the inaccessible pores (by 0.1–0.2%). Under a CO2–acidified water environment, the compressional-wave (P-wave) and shear-wave (S-wave) velocities (5536.7 and 2699.7 m/s) of a core sample containing natural fractures decrease by 8.5 and 8.1%, respectively, whereas both P- and S-wave velocities (6074.1  and 3858.8 m/s) for an intact sample show only ≈0.5% decreases. The interaction also causes more than 50% degradation of the uniaxial compressive strength for the core sample with natural fractures. X-ray microcomputed tomography experiments on three tiny cores (diameter: 1 mm) after 5-day treatment with CO2 (12 MPa) also show that matrix erosion occurs under CO2–acidified water environment but barely occurs without a direct contact with liquid water. Our study suggests that the hydromechanical properties of carbonate caprocks could evolve over the long-term CO2–brine invasion, and it is critical to monitor the CO2–acidified brine interface for a better and long-term evaluation of the caprock integrity.

Description: Summary The objectives of this study are to perform a fundamental analysis of the mutual interactions between crude oil components, water, hydrocarbon solvents, and clays, and to determine the optimum hydrocarbon solvent in solvent steamflooding for a particular reservoir type. The water/oil emulsion formation mechanism in the obtained oil for steam and solvent steamflooding processes has been studied via intermolecular associations between asphaltenes, water, and migrated clay particles. A series of 21 steam and solvent-steamflooding experiments has been conducted, first without any clays in the oil/sand packing, and then using two different clay types in the reservoir rock: Clay 1, which is kaolinite, and Clay 2, which is a mixture of kaolinite and illite. Paraffinic (propane, n-butane,n-pentane,n-hexane,n-heptane) and aromatic (toluene) solvents are coinjected with steam. Cumulative oil recovery is found to decrease in the following order: no clay, Clay 1, Clay 2. Based on the obtained produced oil analyses, Clay 1 and Clay 2 are found to have an affinity with the water and oil phases, respectively. Moreover, the biwettable nature of Clay 2 makes it dispersed in the oil phase toward the oil/water interface, stabilizing the water/oil emulsions. Paraffinic solvent n-hexane is found to be an optimum coinjector for solvent steamflooding in bitumen recovery.

Description: Summary Modeling the dynamic fluid behavior of low-salinity waterflooding (LSWF) at the reservoir scale is a challenge that requires a coarse-grid simulation to enable prediction in a feasible time scale. However, evidence shows that using low-resolution models will result in a considerable mismatch compared with an equivalent fine-scale model with the potential of strong, numerically induced pulses and other dispersion-related effects. This work examines two new upscaling methods that have been applied to improve the accuracy of predictions in a heterogeneous reservoir where viscous crossflow takes place. We apply two approaches to upscaling to bring the flow prediction closer to being exact. In the first method, we shift the effective-salinity range for the coarse model using algorithms that we have developed to correct for numerical dispersion and associated effects. The second upscaling method uses appropriately derived pseudorelative permeability curves. The shape of these new curves is designed using a modified fractional-flow analysis of LSWF that captures the relationship between dispersion and the waterfront velocities. This second approach removes the need for explicit simulation of salinity transport to model oil displacement. We applied these approaches in layered models and for permeability distributed as a correlated random field. Upscaling by shifting the effective-salinity range of the coarse-grid model gave a good match to the fine-scale scenario, while considerable mismatch was observed for upscaling of the absolute permeability alone. For highly coarsened models, this method of upscaling reduced the appearance of numerically induced pulses. On the other hand, upscaling by using a single (pseudo)relative permeability produced more robust results with a very promising match to the fine-scale scenario. These methods of upscaling showed promising results when they were used to scale up fully communicating and noncommunicating layers as well as models with randomly correlated permeability. Unlike documented methods in the literature, these newly derived methods take into account the substantial effects of numerical dispersion and effective concentration on fluid dynamics using mathematical tools. The methods could be applied for other models where the phase mobilities change as a result of an injected solute, such as surfactant flooding and alkaline flooding. Usually these models use two sets of relative permeability and switch from one to another as a function of the concentration of the solute.

Description: Summary In this paper, an extensive series of experiments was performed to investigate the evolution of poromechanical (dry, drained, undrained, and unjacketed moduli), transport (permeability), and strength properties during reservoir depletion and injection in a high-porosity sandstone (Castlegate). An overdetermined set of eight poroelastic moduli was measured as a function of confining pressure (Pc) and pore pressure (Pp). The results showed larger effect on pore pressure at low Terzaghi’s effective stress (nonlinear trend) during depletion and injection. Moreover, the rock sample is stiffer during injection than depletion. At the same Pc and Pp, Biot’s coefficient and Skempton’s coefficient are larger in depletion than injection. Under deviatoric loading, absolute permeability decreased by 35% with increasing effective confining stress up to 20.68 MPa. Given these variations in rock properties, modeling of in-situ-stress changes using constant properties could attain erroneous predictions. Moreover, constant deviatoric stress-depletion/injection failure tests showed no changes or infinitesimal variations of strength properties with depletion and injection. It was found that failure of Castlegate sandstone is controlled by simple effective stress, as postulated by Terzaghi. Effective-stress coefficients at failure (effective-stress coefficient for strength) were found to be close to unity (actual numbers, however, were 1.03 for Samples CS-5 and CS-9 and 1.04 for Sample CS-10). Microstructural analysis of Castlegate sandstone using both scanning electron microscope (SEM) and optical microscope revealed that the changes in poroelastic and transport properties as well as the significant hysteresis between depletion and injection are attributed to the existence and distribution of compliant components such as pores, microcracks, and clay minerals.

Description: Summary In wellbore drilling, it is appreciable to devise methods to study the rheology of high-speed annulus fluid flows. In this paper, a high-speedTaylor-Couette system (TCS) was devised to explore non-Newtonian fluid flow behavior appraised by SiO2 nanoparticles toward friction reduction, power saving, and rheology modeling of nanofluids. Water-based mud (WBM) as an environmentally friendly drilling fluid is investigated by adding SiO2 nanoparticles at four low-volume concentrations of 0.05, 0.1, 0.5, and 1% at speeds from 0 to 1,600 rev/min with 200 rev/min intervals in the TCS. Five rheology models based on the Herschel-Bulkley-Extended (HBE) model and a generalized Reynolds number were optimized to fit with the experimental data. All models except the Newtonian model have fitted all nanofluids with high accuracy, especially Bingham and HBE models. Negative deviation from Darcy friction was avoided for power-law (PL) and Herschel-Bulkley (HB) models using the modification to the generalized Reynolds number. Higher energy saving and enhanced rheology is reported particularly at lower volume concentrations of SiO2 WBM nanofluids. The Darcy friction factor deviated from laminar flow at the generalized Reynolds number beyond 2,000 into turbulent, which is a good indicator for the flow condition of complex non-Newtonian nanofluids in real-lifeapplication.

Description: Summary The complex pore structure and storage mechanism of organic-rich ultratight reservoirs make the hydrocarbon transport within these reservoirs complicated and significantly different from conventional oil and gas reservoirs. A substantial fraction of pore volume in the ultratight matrix consists of nanopores in which the notion of viscous flow may become irrelevant. Instead, multiple transport and storage mechanisms should be considered to model fluid transport within the shale matrix, including molecular diffusion, Knudsen diffusion, surface diffusion, and sorption. This paper presents a diffusion-based semianalytical model for a single-component gas transport within an infinite-actingorganic-rich ultratight matrix. The model treats free and sorbed gas as two phases coexisting in nanopores. The overall mass conservation equation for both phases is transformed into one governing equation solely on the basis of the concentration (density) of the free phase. As a result, the partial differential equation (PDE) governing the overall mass transport carries two newly defined nonlinear terms; namely, effective diffusion coefficient, De, and capacity factor, Φ. The De term accounts for the molecular, Knudsen, and surface diffusion coefficients, and the Φ term considers the mass exchange between free and sorbed phases under sorption equilibrium condition. Furthermore, the ratio of De/Φ is recognized as an apparent diffusion coefficient Da, which is a function of free phase concentration. The nonlinear PDE is solved by applying a piecewise-constant-coefficient technique that divides the domain under consideration into an arbitrary number of subdomains. Each subdomain is assigned with a constant Da. The diffusion-based model is validated against numerical simulation. The model is then used to investigate the impact of surface and Knudsen diffusion coefficients, porosity, and adsorption capacity on gas transport within the ultratight formation. Further, the model is used to study gas transport and production from the Barnett, Marcellus, and New Albany shales. The results show that surface diffusion significantly contributes to gas production in shales with large values of surface diffusion coefficient and adsorption capacity and small values of Knudsen diffusion coefficient and total porosity. Thus, neglecting surface diffusion in organic-rich shales may result in the underestimation of gas production.

Description: Summary Scale inhibitors have been widely used as one of the most efficient methods for sulfate-scale control. To accurately predict the required minimum inhibitor concentration (MIC), we have previously developed several crystallization and inhibition models for pure sulfate scales, including barite, celestite, and gypsum. However, disregarding the wide existence of barium-strontium-sulfate (Ba-Sr-SO4) solid solution in the oil field, no related models have been developed that would lead to large errors in MIC determination. In this study, the induction time of Ba-Sr-SO4 solid solution was measured by laser apparatus with or without different concentrations of scale inhibitor diethylenetriamine penta(methylene phosphonic acid) (DTPMP) at the conditions of barite saturation index (SI) from 1.5 to 1.8, temperature (T) from 40 to 70°C, and [Sr2+]/[Ba2+] ratios from 0 to 15 with celestite SI 

Description: Summary As the crucial step in closed-loop reservoir management, robust life-cycle production optimization is defined as maximizing/minimizing the expected value of a predefined objective (cost) function over geological uncertainties (i.e., uncertainties in the reservoir permeability, porosity, endpoint relative permeability, etc.). However, with robust optimization, there is no control over downside risk defined as the minimum net present value (NPV) among the individual NPVs of the different reservoir models. Yet, field operators generally wish to keep this minimum NPV reasonably large to try to ensure that the reservoir is commercially viable. In addition, the field operator may desire to maximize the NPV of production over a much shorter time period than the life of the reservoir under the limitation of surface facilities (e.g., field liquid and water production rates). Thus, it is important to consider multiobjective robust production optimization with nonlinear constraints and when geological uncertainties are incorporated. The three objectives considered in this paper are; to maximize the average life-cycle NPV, to maximize the average short-term NPV, and to maximize the minimum NPV of the set of realizations. Generally, these objectives are in conflict; for example, the well controls that give a global maximum for robust life-cycle production optimization do not usually correspond to the controls that maximize the short-term average NPV of production. Moreover, handling the nonlinear state constraints (e.g., field liquid production rates and field water production rates for the bottom-hole pressure controlled producers in the robust production optimization) is also a challenge because those nonlinear constraints should be satisfied at each control steps for each geological realization. To provide potential solutions to the multiobjective robust optimization problem with state constraints, we developed a modified lexicographic method with a minimizing-maximum scheme to attempt to obtain a set of Pareto optimal solutions and to satisfy all nonlinear constraints. We apply the sequential quadratic programming filter with modified stochastic gradients to solve a sequence of optimization problems, where each solution is designed to generate a single point on the Pareto front. In the modified lexicographic method, the objective is always considered to be the primary objective, and the other objectives are considered by specifying bounds on them to convert them to state constraints. The temporal damping and truncation schemes are applied to improve the quality of the stochastic gradient on nonlinear constraints, and the minimizing–maximum procedure is applied to enforce constraints on the normal state constraints. The main advantage that the modified lexicographic method has over the standard lexicographic method is that it allows for the generation of potential Pareto optimal points, which are uniformly spaced in the values of the second and/or third objective that one wishes to improve by multiobjective optimization.

Description: Summary This work aims to address a challenge posed by recent observations of tightly spaced hydraulic fractures in core samples from the hydraulic fracturing test site (HFTS) in the Middle Wolfcamp Formation. Many fractures in retrieved cores have subfoot spacing, which is at odds with conventional models in which usually one hydraulic fracture is initiated per cluster. Models assuming a single fracture at each cluster, although a common practice, often predict excessive fracture propagation that is inconsistent with microseismic observation. Here, we aim to develop a numerical approach to effectively account for densely spaced hydraulic fractures in field-scale simulations. Because it is impractical to explicitly model all aforementioned fractures, we develop a new upscaling law that enables existing simulation tools to predict reservoir response to fracture swarms. The upscaling law is derived based on an energy equivalence argument and validated through multiscale simulations using a high-fidelity code, GEOS. The swarming fractures are first modeled with a spacing that is much smaller than the cluster spacing; these fractures are then approximated by an upscaled, single fracture based on the proposed upscaling law. The upscaled fracture is shown to successfully match the energy input rate and produce the total fracture aperture and average propagation length of the explicitly simulated swarm. Afterward, the upscaling approach is further implemented in 3D field-scale simulations and validated against the HFTS microseismic data of a horizontal well. Our results show that hydraulic fracture swarming can significantly affect fracture propagation behaviors compared with the propagation of single fractures as assumed by conventional modeling approaches. Under the considered situations, the conventional treatment yields fast propagation speed that far exceeds that indicated by the microseismic data. We also illustrate that this discrepancy can be reduced readily through the implementation of the upscaling law. Our results demonstrate the importance of accounting for the fracture swarming effect in field-scale simulations and the efficacy of this approach to enable realistic predictions of reservoir responses to fracture swarms, without the need to model tightly spaced fractures individually.

Description: Summary In conventional wellbore-integrity analysis, the cement sheath’s initial state of stress and transient thermoporoelastic effects are often neglected. However, the initial state of stress is prerequisite information for accurately predicting the safe operating conditions that prevent a cemented well from being damaged. In addition, transient thermoporoelastic effects can have a profound effect on when damage will occur. In this paper, we propose a model that includes these effects to predict the safe operating pressures and temperatures that will prevent cement-sheath failure. For the initial state of stress, we proposed an empirical model using measurements. Subsequent stress changes are evaluated by a fully coupled transient thermoporoelastic model to analyze the mechanical behavior of the cement sheath. We predict the safe operating envelope (SOE) for shear, tensile, and debonding cement-sheath failures caused by pressure and temperature perturbations after the cement sets. Our model predicts that pore pressure is a key factor for cement failure, especially for rapid temperature changes. If the formation is low permeability, the transient pore pressures are amplified, causing the risk of damage to increase. Compared with conventional thermoelastic models, the thermoporoelastic model predicts a smaller SOE when heating the internal casing fluid and a larger envelope when cooling the internal casing fluid. Finally, the heating rate was considered with respect to field applications. The heating rate was also considered, and slower heating/cooling rates can prevent damage to the cement sheath. Finally, the thermoporoelastic model was applied to explain several laboratory and field experiments and achieved good matches. Correction (in press): The Fig. 4 caption has been changed from "Preliminary cyclic loading test" to "Experimental results for cement strength" to correct an error in labeling. No other content has been changed.

Description: Summary Nanoparticles have improved a surfactant's ability to create long-lasting foam. Recent studies have widely recommended the use of silica nanoparticles to enhance foam stability. This paper presents an experimental investigation of a new and highly effective alpha olefin sulfonate (AOS)–multiwalled carbon nanotube (MWCNT) system for mobility control during gas enhanced oil recovery (EOR) operations. The new AOS–MWCNT system was evaluated for its foam stability at 150°F using a high-pressure view cell. The MWCNT was obtained as solid particles of aspect ratio up to 100 and silica nanoparticles of median size of 118 nm. The foam system was optimized for its maximum half-life by varying the concentration of the AOS and the nanotube from 0.2 to 1% and 250 to 1,000 ppm, respectively. Compatibility testing with salts was done as well. Coreflood experiments with 1.5-in.-diameter, 6-in.-long Berea sandstone cores were run to calculate the mobility reduction factor at 150°F. Nitrogen foam was injected into the core at 80% foam quality in the tertiary recovery mode, and the pressure drop across the core was measured. The formation brine had a salinity of 5 wt% sodium chloride (NaCl), and the foaming solutions were prepared with 2 wt% NaCl. The optimal concentrations of the AOS solution and the nanotubes for maximum foam stability were determined to be 0.5% and 500 ppm, respectively. The optimized AOS–MWCNT system yielded 60% greater nitrogen foam half-life (32 minutes) than an optimized AOS–silica system at 150°F. The foam half-life of a stand-alone 0.5% AOS solution was 7 minutes. In the presence of crude oil, the foam half-life decreased for all the tested systems. Coreflood experiments at 150°F showed a significant increase in the mobility reduction factor when the new AOS–MWCNT system was used as the foamer instead of stand-alone AOS or AOS–silica system. The new foaming system was stable through the duration of the experiment, yielding foam in the effluent samples. There was no formation damage observed. Salt tolerance for the MWCNT nanofluid was higher than the silica nanofluid. Foam needs to be stable for long periods of time to ensure effective mobility control during gas injection for EOR. This paper investigates a new highly effective AOS-multiwalled carbon nanotube system that outperforms the AOS–silica foaming systems in terms of foam stability and mobility control at 150°F.

Description: Summary X-ray computerized tomography (CT) is a nondestructive method of providing information about the internal composition and structure of whole core reservoir samples. In this study we propose a method to classify lithology. The novelty of this method is that it uses statistical and textural information extracted from whole core CT images in a supervised learning environment. In the proposed approaches, first-order statistical features and textural grey-levelco-occurrence matrix (GLCM) features are extracted from whole core CT images. Here, two workflows are considered. In the first workflow, the extracted features are used to train a support vector machine (SVM) to classify lithofacies. In the second workflow, a principal component analysis (PCA) step is added before training with two purposes: first, to eliminate collinearity among the features and second, to investigate the amount of information needed to differentiate the analyzed images. Before extracting the statistical features, the images are preprocessed and decomposed using Haar mother wavelet decomposition schemes to enhance the texture and to acquire a set of detail images that are then used to compute the statistical features. The training data set includes lithological information obtained from core description. The approach is validated using the trained SVM and hybrid (PCA + SVM) classifiers to predict lithofacies in a set of unseen data. The obtained results show that the SVM classifier can predict some of the lithofacies with high accuracy (up to 91% recall), but it misclassifies, to some extent, similar lithofacies with similar grain size, texture, and transport properties. The SVM classifier captures the heterogeneity in the whole core CT images more accurately compared with the core description, indicating that the CT images provide additional high-resolution information not observed by manual core description. Further, the obtained prediction results add information on the similarity of the lithofacies classes. The prediction results using the hybrid classifier are worse than the SVM classifier, indicating that low-power components may contain information that is required to differentiate among various lithofacies.

Description: Summary Artificial lift systems are widely used in oil production, of which sucker rod pumps are conceptually among the simpler ones. The reciprocating movement of the plunger triggers the opening and closing of two ball valves, allowing fluid to be pumped to the surface. Their built-in ball valves are subject to long-time erosion and fail as a consequence of this damage mechanism. Understanding the principal damage mechanisms requires a thorough examination of the fluid dynamics during the opening and closing action of these valves. In this article, we present a fluid-structure interaction model that simultaneously computes the fluid flow in the traveling valve (TV), the standing valve (SV), and the chamber of sucker rod pumps during a full pump cycle. The simulations shed light on the causes of valve damage for standard and nonideal operating conditions of the pump. In particular, our simulations based on real pump operating envelopes reveal that the so-called “midcycle valve closure” is likely to occur. Such additional closing and opening events of the valves multiply situations in which the flow conditions are harmful to the individual pump components, leading to efficiency reduction and pump failure. This mechanism, hitherto unreported in the literature, is believed to constitute the primary cause of long-term valve damage. Our finite element method-based computational-fluid-dynamics model can accurately describe the opening and closing cycles of the two valves. For the first time, this approach allows an analysis of real TV speed versus position plots, usually called pump cards. The effects of stroke length, plunger speed, and fluid parameters on the velocity and pressure at any point and time inside the pump can thus be investigated. Identifying the damage-critical flow parameters can help suggest measures to avoid unfavorable operating envelopes in future pump designs. Our flow model may support field operations throughout the entire well life, ranging from improved downhole pump design to optimized pump operation or material selections. It can aid the creation of an ideal interaction between the valves, thus avoiding midcycle valve closure to drastically extend the mean time between failures of sucker rod pumps. Finally, our simulation approach will speed up new pump component development while greatly reducing the necessity for costly laboratory testing.

Description: Summary Phase behavior and physical properties including saturation pressures, swelling factors (SFs), phase volumes, dimethyl ether (DME) partition coefficients, and DME solubility for heavy-oil mixtures containing polar substances have been experimentally and theoretically determined. Experimentally, novel phase behavior experiments of DME/water/heavy-oil mixtures spanning a wide range of pressures and temperatures have been conducted. More specifically, a total of five pressure/volume/temperature (PVT) experiments consisting of two tests of DME/heavy-oil mixtures and three tests of DME/water/heavy-oil mixtures have been performed to measure saturation pressures, phase volumes, and SFs. Theoretically, the modified Peng-Robinson equation of state (EOS) (PR EOS) together with the Huron-Vidal mixing rule, as well as the Péneloux et al. (1982)volume-translation strategy, is adopted to perform phase-equilibrium calculations. The binary-interaction parameter (BIP) between the DME/heavy-oil pair, which is obtained by matching the measured saturation pressures of DME/heavy-oil mixtures, works well for DME/heavy-oil mixtures in the presence and absence of water. The new model developed in this work is capable of accurately reproducing the experimentally measured multiphase boundaries, phase volumes, and SFs for the aforementioned mixtures with the root-mean-squared relative error (RMSRE) of 3.92, 9.40, and 0.92%, respectively, while it can also be used to determine DME partition coefficients and DME solubility for DME/water/heavy-oil systems.

Description: Summary Asphaltene deposition triggers serious flow assurance issues and can significantly restrict the production capacity. Because of the complexity associated with asphaltene deposition that includes several mechanisms acting simultaneously, an accurate prediction of asphaltene blockage along the wellbore requires integration of asphaltene precipitation, aggregation, and deposition. In this work, an integrated simulation approach is proposed to predict the asphaltene deposition profile along the wellbore. The proposed approach is novel because it integrates various deposition patterns of particulate flow (which depends on hydrodynamics) with aggregation processes to investigate how the distribution of asphaltene particle size varies (governed by molecular dynamics) after being precipitated out of the oil phase (controlled by thermodynamics). To improve the predictability capability of simulations, a direct input from the wellbore flow simulator is used to update the velocity profile after the wellbore radius changes beyond a certain predefined threshold. The fraction of asphaltene precipitation is determined using the asphaltene solubility model and combined with aggregation models to feed into deposition calculations. Wellbore blockage was examined for two cases with and without the aggregation mechanism included. A sensitivity analysis was carried out to study parameters that affect the severity of blockage, such as range of pressure-temperature along the wellbore, flow velocity, and radial distribution of asphaltene particles. The simulation approach proposed in this paper provides an in-depth understanding of the wellbore flow assurance issues caused by asphaltene deposition and thus provides useful insights for improving the predictions of production performance.

Description: Summary The effect of axial flow of power-law drilling fluids on frictional pressure loss under turbulent conditions in eccentric annuli is investigated. A numerical model is developed to simulate the flow of Newtonian and power-law fluids for eccentric annular geometries. A turbulent eddy-viscosity model based on the mixing-length approach is proposed, where a damping constant as a function of flow parameters is presented to account for the near-wall effects. Numerical results including the velocity profile, eddy viscosity, and friction factors are compared with various sets of experimental data for Newtonian and power-law fluids in concentric and eccentric annular configurations with diameter ratios of 0.2 to 0.8. The simulation results are also compared with a numerical study and two approximate models in the literature. The results of extensive simulation scenarios are used to obtain a novel correlation for estimation of the frictional pressure loss in eccentric annuli under turbulent conditions. Two new correlations are also presented to estimate the maximum axial velocity in the wide and narrow sections of eccentric geometries.

Description: Summary Gas kicks occur frequently in deepwater drilling because of the extremely narrow mud-weight window [minimum 0.01 specific gravity (sg)]. The traditional kick-detection method mainly relies on the driller's analysis of monitored compound comprehensive mud-logging data. However, the traditional method has significant time lag, including missed and false detection, and often leads to severe gas influxes during deepwater drilling. A novel machine-learning (ML) model is presented here using pilot-scale rig data combined with surface-riser-downhole monitoring for gas-kick early detection and risk classification. A series of pilot-scaletest-well experiments (a total of 108 tests) are performed to simulate deepwater gas kicks and produce a multisource data set through fusion of comprehensive mud-logging data from surface monitoring, acoustic data from riser-monitoring technologies, and measurement-while-drilling data [e.g., bottomhole pressure (BHP)] from downhole monitoring technologies. During these experiments, the deepwater blowout preventer (BOP) is simulated using a variable cross section of crossover (X/O; equipped with booster-flow pipes); the Coriolis flowmeter is installed in the mud-return pipe to accurately measure flow out; the acoustic wave sensors are installed outside of the riser section (X/O) to monitor gas migration; and the downhole memory pressure gauges are installed to monitor BHP. Next, data preparation and data analysis are performed including raw-data exploration, data cleaning, signal/noise-ratio (SNR) analysis, feature scaling, outlier detection, and feature engineering. Further, a novel and improved data-labeling criterion for gas-kick alarms is proposed, with six levels (displayed using different colors) instead of two-state alarms (“kick” or “no kick”). The proposed gas-kick-alarm classification is in accordance with the actual field practices. Subsequently, four ML algorithms—decision tree (DT), k-nearest neighbors (KNN), support vector machine (SVM), and long short-term memory (LSTM)—are developed through the complete workflow, beginning with the data allocation and followed by building, evaluation, and optimization of each ML model. Because the LSTM recurrent neural network (RNN) algorithm showed the best performance, it is selected and deployed to early detect gas kicks and classify the corresponding kick alarms. The recall for gas-kick levels corresponding to Risk 0, Risk 1, Risk 2, Risk 3, Risk 4, and Risk 5 are 0.92, 0.93, 0.91, 0.91, 0.92 and 0.92, respectively. Because recall for each gas-kick-alarm level is greater than 0.9, it ensures rare false negatives (FNs) during kick detection. The accuracy, precision, recall, and f1 score of the deployed LSTM model in the testing data set is 91.6%, 0.93, 0.92 and 0.92, respectively. Further, the detection time delay is approximately 2 to 7 seconds only, which provides an improved time margin to take appropriate safety measures, promptly deal with a gas kick through a well-control program, and prevent a potential blowout during deepwater drilling.

Description: Summary An explicit solution to the general 3D point-to-target problem based on the minimum curvature method has been sought for more than four decades. The general case involves the trajectory's start and target points connected by two circular arcs joined by a straight line with the position and direction defined at both ends. It is known that the solutions are multivalued, and efficient iterative schemes to find the principal root have been established. This construction is an essential component of all major trajectory construction packages. However, convergence issues have been reported in cases where the intermediate tangent section is either small or vanishes and rigorous mathematical conditions under which solutions are both possible and are guaranteed to converge have not been published. An implicit expression has now been determined that enables all the roots to be identified and permits either exact or polynomial-type solution methods to be used. Most historical attempts at solving the problem have been purely algebraic, but a geometric interpretation of related problems has been attempted, showing that a single circular arc and a tangent section can be encapsulated in the surface of a horn torus. These ideas have now been extended, revealing that the solution to the general 3D point-to-target problem can be represented as a 10th-orderself-intersecting geometric surface, characterized by the trajectory's start and end points, the radii of the two arcs, and the length of the tangent section. An outline of the solution's derivation is provided in the paper together with complete details of the general expression and its various degenerate forms so that readers can implement the algorithms for practical application. Most of the degenerate conditions reduce the order of the governing equation. Full details of the critical and degenerate conditions are also provided, and together these indicate the most convenient solution method for each case. In the presence of a tangent section, the principal root is still most easily obtained using an iterative scheme, but the mathematical constraints are now known. It is also shown that all other cases degenerate to quadratic forms that can be solved using conventional methods. It is shown how the general expression for the general point-to-target problem can be modified to give the known solutions to the 3D landing problem and how the example in the published works on this subject is much simplified by the geometric, rather than algebraic treatment.

Description: Summary Resistivity measurements are a major input into hydrocarbon reserve estimation and are usually described by Archie’s laws. In this study, we use digital rock physics to analyze the mechanisms of non-Archie and Archie behavior of formation factor (FF) and resistivity index (RI) of low-porosity Fontainebleau (FB) sandstone for ambient conditions and under high confining pressure, respectively. FB sandstone was imaged by micro-X-ray computed tomography (micro-CT) at a resolution of 1 µm. Subresolution details of the grain contact width distribution along with their length were extracted from a set of scanning electron microscope (SEM) images. The nanoscale aperture of grain contacts, which is below tomogram resolution, is accounted for in micro-CT-based numerical calculations by assigning effective porosity and conductivity to individual voxels of the extracted grain contact network. A porosity reduction of grain contacts and open pore space as a function of applied confining pressure is introduced, capturing the pressure dependence. The concept was implemented by grain contact labeling, introducing an additional phase derived from a Euclidean distance transform (EDT). Subvoxel stress-strain effects were incorporated by attributing all compressibility effects to the pore space (open pore space and grain contacts), treating the solid phase as perfectly rigid. Voxel-scale input conductivities are assigned using Archie’s law followed by solving the Laplace equation for sample-scale rock resistivity and RI directly on the segmented image using the finite element method. For the numerical modeling of the FF and RI of low-porosity FB sandstone as a function of confining pressure, which depends on subresolution features, a set of hypotheses were tested. These are based on two segmentation scenarios incorporating the measured contact aperture distribution from SEM analysis—a homogeneous aperture-based segmentation by assuming all grain contacts as an average constant value and a heterogeneous aperture-based segmentation assigning two groups of grain contact apertures. The segmentation scenarios enable homogeneous and heterogeneous morphological change of grain contacts due to confining pressure effects. Furthermore, partial saturation of grain contacts is considered. In all cases, strong water-wetness was assumed, and discretization effects were analyzed carefully. The numerical results highlight the relative contribution of each of two conductive components of FB sandstone (open pores vs. grain contacts) over the full range of partial saturations. Of importance is the connectivity of the system, with discretization effects having a significant effect on FF, but a small effect on the RI. Grain contacts and confining pressure are found to have a significant impact on RI behavior of low-porosity FB sandstone. Both the grain contact network with homogeneous aperture and the heterogeneous grain contact network are able to describe experimental observations. However, it is not sufficient to assume a homogeneous change in contact area, and an inhomogeneous deformation of grain contact zones is required to match the experiment.

Description: Summary An experimental and theoretical investigation of surfactant-stabilized oil/water emulsion characteristics was carried out under water sweep (WS) and oil sweep (OS) conditions. Both hydrophilic and hydrophobic surfactants were used, with concentrations less than and more than the critical micelle concentration (CMC). Experimental data were acquired for detection of the phase-inversion region, which was measured simultaneously by several independent methods. These include a circular differential dielectric sensor (C-DDS), a rectangular differential dielectric sensor (R-DDS) (both sensors accurately detect the phase-inversion region), a pressure transducer, and a mass flowmeter. The addition of an emulsifier surfactant to an oil/water mixture generated a stable emulsion, which resulted in a phase-inversion delay. For water-continuous to oil-continuous flow, a hydrophilic surfactant was a better emulsifier, while for oil-continuous to water-continuous flow, a hydrophobic surfactant was a better emulsifier for creating more stable emulsions. The surfactant/oil/water emulsion resulted in an increase of the dispersed-phase volume fraction required for phase inversion, as compared to the case of oil/water dispersions without surfactant. For emulsions with surfactant concentrations above CMC, the presence of micelles contributed to further delay of the phase inversion, as compared to those with surfactant concentrations below CMC. The phase-inversion region exhibits a hysteresis between the OS and WS runs, below CMC and above CMC, which was due to the difference in droplet sizes caused by different breakup and coalescence processes for oil-continuous and water-continuousflow. This research shows that the DDS is an efficient instrumentation that can be used to detect the region where the emulsion phase inversion is expected to occur. Moreover, the experimental results and the pertinent analysis and discussion provide useful insights for a more informed design of surface facilities (including emulsion separators) in oil and gas production operations.

Description: Summary North American shale drilling is a fast-paced environment where downhole drilling equipment is pushed to the limits for the maximum rate of penetration (ROP). Downhole mud motor power sections have rapidly advanced to deliver more horsepower and torque, resulting in different downhole dynamics that have not been identified in the past. High-frequency (HF) compact drilling dynamics recorders embedded in the drill bit, mud motor bit box, and motor top subassembly (top-sub) provide unique measurements to fully understand the reaction of the steerable-motor power section under load relative to the type of rock being drilled. Three-axis shock, gyro, and temperature sensors placed above and below the power section measure the dynamic response of power transfer to the bit and associated losses caused by back-drive dynamics. Detection of back-drive from surface measurements is not possible, and many measurement-while-drilling (MWD) systems do not have the measurement capability to identify the problem. Motor back-drive dynamics severity is dependent on many factors, including formation type, bit type, power section, weight on bit, and drillpipe size. The torsional energy stored and released in the drillstring can be high because of the interaction between surface rotation speed/torque output and mud motor downhole rotation speed/torque. Torsional drillstring energy wind-up and release results in variable power output at the bit, inconsistent rate of penetration, rapid fatigue on downhole equipment, and motor or drillstring backoffs and twistoffs. A new mechanism of motor back-drive dynamics caused by the use of an MWD pulser above a steerable motor has been discovered. HF continuous gyro sensors and pressure sensors were deployed to capture the mechanism in which a positive mud pulser reduces as much as one-third of the mud flow in the motor and bit rotation speed, creating a propensity for a bit to come to a complete stop in certain conditions and for the motor to rotate the drillstring backward. We have observed the backward rotation of a polycrystalline diamond compact (PDC) drill bit during severe stick-slip and back-drive events (−50 rev/min above the motor), confirming that the bit rotated backward for 9 milliseconds (ms) every 133.3 ms (at 7.5 Hz), using a 1,000-Hz continuous sampling/recording in-bit gyro. In one field test, multiple drillstring dynamics recorders were used to measure the motor back-drive severity along the drillstring. It was discovered that the back-drive dynamics are worse at the drillstring, approximately 1,110 ft behind the bit, than these measured at the motor top-sub position. These dynamics caused drillstring backoffs and twistoffs in a particular field. A motor back-drive mitigation tool was used in the field to compare the runs with and without the mitigation tool while keeping the surface drilling parameters nearly the same. The downhole drilling dynamics sensors were used to confirm that the mitigation tool significantly reduced stick-slip and eliminated the motor back-drive dynamics in the same depth interval. Detailed analysis of the HF embedded downhole sensor data provides an in-depth understanding of mud motor back-drive dynamics. The cause, severity, reduction in drilling performance and risk of incident can be identified, allowing performance and cost gains to be realized. This paper will detail the advantages to understanding and reducing motor back-drive dynamics, a topic that has not commonly been discussed in the past.

Description: Summary A continuum hydrodynamic model with immersed solid/fluid interface is developed for simulating calcite dissolution by hydrochloric acid (HCl) at the pore scale, and is most accurate for a mass-transfer-controlled dissolution regime under laminar flow conditions. The model uses averaged Navier-Stokes equations to model momentum transfer in porous media and adopts a theoretically developed mass-transfer formulation with assumptions. The model includes no fitting parameter and is validated using experimental results. The findings of previous research and existing models are briefly discussed and their shortcomings and advantages are elucidated. The present model is used in some pore-scale simulations on hypothetical but realistic cases, investigating the evolution of Darcy-scale permeability. Darcy-scale permeability exhibits totally different functionality of porosity in different dissolution regimes.

Description: Summary Interpretation of sonic data acquired by a logging-while-drilling (LWD) tool or wireline tool in cased holes is complicated by the presence of drillpipe or casing because those steel pipes can act as a strong waveguide. Traditional solutions, which rely on using a frequency bandpass filter or waveform arrival-time separation to filter out the unwanted pipe mode, often fail when formation and pipe signals coexist in the same frequency band or arrival-time range. We hence developed a physics-driven machine-learning-based method to overcome the challenge. In this method, two synthetic databases are generated from a general root-findingmode-search routine on the basis of two assumed models: One is defined as a cemented cased hole for a wireline scenario, and the other is defined as a steel pipe immersed in a fluid-filled borehole for the logging-while-drilling scenario. The synthetic databases are used to train neural network models, which are first used to perform global sensitivity analysis on all relevant model parameters so that the influence of each parameter on the dipole dispersion data can be well understood. A least-squares inversion scheme using the trained model was developed and tested on synthetic cases. The scheme showed good results, and a reasonable uncertainty estimate was made for each parameter. We then extended the application of the trained model to develop a method for automated labeling and extraction of the dipole flexural dispersion mode from other disturbances. The method combines the clustering technique with the neural-network-model-based inversion and an adaptive filter. Testing on field data demonstrates that the new method is superior to traditional methods because it introduces a mechanism from which unwanted pipe mode can be physically filtered out. This novel physics-driven machine-learning-based method improved the interpretation of sonic dipole dispersion data to cope with the challenge brought by the existence of steel pipes. Unlike data-driven machine learning methods, it can provide global service with just one-time offline training. Compared with traditional methods, the new method is more accurate and reliable because the processing is confined by physical laws. This method is less dependent on input parameters; hence, a fully automated solution could be achieved.

Description: Summary Total field strength, declination, and dip angle of the Earth's magnetic field, in conjunction with gravity, are used by magnetic-survey tools to determine a wellbore's location. Magnetic field values may be obtained from global models that, depending on the model, have a wide range of spatial resolution at the Earth's surface from large scale (3000 km) to small scale (28 km). The magnetic field varies continuously in both time and space, so no model can fully capture the complexity of all sources; hence, there are uncertainties associated with the values provided. The SPE Wellbore Positioning Technical Section/Industry Steering Committee on Wellbore Surveying Accuracy (ISCWSA) published their original measurement-while-drilling (MWD) error model in 2000. Such models and uncertainties define positional error ellipsoids along the wellbore, which assist the driller in achieving their geological target, in addition to aiding collision avoidance. With the recent update to Revision 5 of the ISCWSA error model, we have reassessed the uncertainties associated with our latest high-resolution global magnetic field model. We describe the derivation of location-specific global and random uncertainties for use with predicted geomagnetic values from high-resolution models within magnetic MWD survey-tool-error models. We propose a sophisticated approach to provide realistic values at different locations around the globe; for example, we determine separate errors for regions where the models have high spatial resolution from aeromagnetic data compared to regions where only satellite data are available. The combined uncertainties are freely available via a web service with which the user can also see how they vary with time. The use of the revised uncertainty values in the MWD-error model, in most cases, reduces the positional error ellipsoids and allows better use of the increased accuracy from recent improvements in geomagnetic modeling. This is demonstrated using the new uncertainty values in the MWD-error model for three standard ISCWSA well profiles. A fourth theoretical well offshore Brazil where the vertical magnetic field is weak shows that with drillstring interference correction relying on the more uncertain magnetic dip, the positional error ellipsoids can increase. This is clearly of concern for attaining geological targets and collision avoidance.

Description: Summary Numerical fidelity is required when using simulations to predict enhanced-oil-recovery (EOR) processes. In this paper, we investigate the conditions that lead to numerical errors when simulating low-salinity (LS) waterflooding (LSWF). We also examine how to achieve more accurate simulation results by scaling up the flow behavior in an effective manner. An implicit finite-difference numerical solver was used to simulate LSWF. The accuracy of the numerical solution has been examined as a function of changing the length of the grid cell and the timestep. Previously we have shown that numerical dispersion induces a physical retardation such that the LS front slows down while the formation water front speeds up. We also report for the first time that pulses can be generated as numerical artifacts in coarsely gridded simulations of LSWF. These effects reflect the interaction of dispersion, the effective-salinity range, and the use of upstream weighting during calculation, and can corrupt predictions of flow behavior. The effect of the size of the timestep was analyzed with respect to the Courant condition, traditionally related to explicit numerical schemes and also numerical stability conditions. We also investigated some of the nonlinear elements of the simulation model, such as the differences between the concentrations of connate water salinity and the injected brine, effective-salinity-concentration range, and the net mobility change on fluids through changing the salinity. We report that to avoid pulses it is necessary, but not sufficient, to meet the Courant condition relating timestep size to cell size. We have also developed two approaches that can be used to scale up simulations of LSWF and tackle the numerical problems. The first method is dependent on a mathematical relationship between the fractional flow, effective-salinity range, and the Péclet number and treats the effective-salinity range as a pseudofunction. The second method establishes an unconventional proxy method equivalent to pseudorelative permeabilities. A single table of pseudorelative permeability data can be used for a waterflood instead of two tables, as is usual for LSWF. This is a novel approach that removes the need for relative permeability interpolation during the simulation. Overall, by avoiding numerical errors, we help engineers to more efficiently and accurately assess the potential for improving oil recovery using LSWF and thus optimize field development. We also avoid the numerical pulses inherent in the traditional LSWF model.

Description: Summary Time-lapse-seismic-data assimilation has been drawing the reservoir-engineering community's attention over the past few years. One of the advantages of including this kind of data to improve the reservoir-flow models is that it provides complementary information compared with the wells' production data. Ensemble-based methods are some of the standard tools used to calibrate reservoir models using time-lapse seismic data. One of the drawbacks of assimilating time-lapse seismic data involves the large data sets, mainly for large reservoir models. This situation leads to high-dimensional problems that demand significant computational resources to process and store the matrices when using conventional and straightforward methods. Another known issue associated with the ensemble-based methods is the limited ensemble sizes, which cause spurious correlations between the data and the parameters and limit the degrees of freedom. In this work, we propose a data-assimilation scheme using an efficient implementation of the subspace ensemble randomized maximum likelihood (SEnRML) method with local analysis. This method reduces the computational requirements for assimilating large data sets because the number of operations scales linearly with the number of observed data points. Furthermore, by implementing it with local analysis, we reduce the memory requirements at each update step and mitigate the effects of the limited ensemble sizes. We test two local analysis approaches: one distance-based approach and one correlation-based approach. We apply these implementations to two synthetic time-lapse-seismic-data-assimilation cases, one 2D example, and one field-scale application that mimics some of the real-field challenges. We compare the results with reference solutions and with the known ensemble smoother with multiple data assimilation (ES-MDA) using Kalman gain distance-based localization. The results show that our method can efficiently assimilate time-lapse seismic data, leading to updated models that are comparable with other straightforward methods. The correlation-based local analysis approach provided results similar to the distance-based approach, with the advantage that the former can be applied to data and parameters that do not have specific spatial positions.

Description: Summary Horizontal drilling and hydraulic fracturing are recognized as the most efficient techniques to enhance recovery in shale-gas reservoirs. Because of the exploitation difficulties and complex flow mechanism in shale gas, it is imperative to focus on the optimization of fracturing parameters. However, most of the current heuristic algorithms follow the principle that the variable dimension is constant during iteration, which leads to poor performance when dealing with dimension-varying problems. The optimization of fracturing parameters can be regarded as a typical dimension-varying problem when considering the difference among fracture properties such as half-length and conductivity. Thus an improved algorithm named modified variable-lengthparticle-swarm optimization (PSO) (VPSO) (MVPSO) was proposed to automatically select the optimal fracturing parameters: the number of fractures as well as the corresponding fracture properties. Then, MVPSO was verified and compared with VPSO by several benchmarks. In addition, a gas/water two-phase model considering gas-adsorption and Knudsen-diffusion effects was used to describe the shale-gas flow in matrix and fracture domains. An embedded discrete-fracture model (EDFM) was applied to model the hydraulic-fracture geometries and fractal methods were adopted to generate the fracture networks. The results indicated that MVPSO showed better performance in both convergence speed and accuracy than that of VPSO, which also provided a new perspective for the optimization of fracturing parameters. Besides, the multispindle-shapedfracture-distribution pattern reached a higher net-present-value (NPV) contrast to that of homogeneous fracture distribution. The decrease of gas price leads to smaller and more nonuniform half-lengthdistribution.

Description: Summary It is common to inject acidic stimulation fluids into oil-bearing carbonate formations to enhance well productivity. This process of matrix acidizing is designed to maximize the propagation of wormholes into the formation by optimizing the injection parameters, including acid-injection rate and volume. Previous studies have suggested that saturation conditions, permeability, heterogeneity, temperature, and pressure can significantly affect the design of matrix-acidizing treatments. However, laboratory studies’ results are inconsistent in their conclusions and are mostly limited to water-saturated cores. In this work, we designed a systematic experimental study to evaluate the impact of multiphase flow on the acidizing process when injecting 15 wt% hydrochloric acid (HCl) into crude-oil-saturated Indiana Limestone cores. The results reveal the following: Contrary to published literature for water-saturated cores, acidizing in partially oil-saturatedhigh-permeability cores at high pressure requires less acid volume than in low-permeability cores; lower-pressure acid injection results in more efficient wormhole propagation in low-permeability cores compared to high-pressure acid injection; acidizing in low- and high-permeability cores at low pressure leads to similar efficiency; and wormholing is more effective in partially oil-saturated cores, resulting in multiple parallel branches as compared to inefficient leakoff in water-saturatedcores.

Description: Summary In this paper we present a methodology to superimpose the American Petroleum Institute (API) uniaxial and triaxial limits on tubular design limits plots (API TR 5C3 2018). Complications caused by a recent change of axis are resolved, producing a practical design limits plot that avoids the horizontal shift of the API vertical limits, which is currently the industry standard. The commonly used slanted ellipse is compared against an adaptation of the circle of plasticity in the form of a horizontal ellipse, showing the convenience of this last one with examples. After the current official collapse formulation was made part of the main body of standard API TR 5C3 (2018), the horizontal axis on the standard industry well tubular design limits plot changed. The present study evaluates this redefinition of the horizontal axis. One consequence of this modification is a difficulty plotting the API tension and compression limits. The API horizontal limits (uniaxial burst and collapse) are found to be independent of load situation, whereas the API vertical design limits (uniaxial tension and compression) are dependent on inside and outside tubular pressures. The approaches used by commercial software and industry publications to solve this challenge are reviewed. A new design methodology is developed to link API uniaxial limits to the triaxial theory. One main objective of the study is to establish a mathematical relationship between API tubular design limits and the von Mises triaxial theory (API TR 5C3 2018). A methodology that allows plotting the API uniaxial force limits on the design limits plot is developed. The study also shows that the results obtained from the industry standard slanted ellipse are identical to those obtained from the horizontal ellipse and circle. One important difference is that the slanted ellipse is based on the zero axial stress datum, whereas the horizontal ellipse/circle uses the neutral axial stress datum. The horizontal ellipse/circle is well suited for calculations involving buckling, compatible with the information used in field operations, and its formulations are less complicated than the tilted ellipse. Therefore, attention is called to the use of the horizontal ellipse/circle in well tubular design.

Description: Summary Numerical simulation of coupled multiphase multicomponent flow and transport in porous media is a crucial tool for understanding and forecasting of complex industrial applications related to the subsurface. The discretized governing equations are highly nonlinear and usually need to be solved with Newton’s method, which corresponds with high computational cost and complexity. With the presence of capillary and gravity forces, the nonlinearity of the problem is amplified even further, which usually leads to a higher numerical cost. A recently proposed operator-based linearization (OBL) approach effectively improves the performance of complex physical modeling by transforming the discretized nonlinear conservation equations into a quasilinear form according to state-dependent operators. These operators are approximated by means of a discrete representation on a uniform mesh in physical parameter space. Continuous representation is achieved through the multilinear interpolation. This approach provides a unique framework for the multifidelity representation of physics in general-purpose simulation. The applicability of the OBL approach was demonstrated for various energy subsurface applications with multiphase flow of mass and heat in the presence of buoyancy and diffusive forces. In this work, the OBL approach is extended for multiphase multicomponent systems with capillarity. Through the comparisons with a legacy commercial simulator using a set of benchmark tests, we demonstrate that the extended OBL scheme significantly improves the computational efficiency with the controlled accuracy of approximation and converges to the results of the conventional continuous approach with an increased resolution of parametrization.

Description: Summary We studied the applicability of a gradient-boostingmachine-learning (ML) algorithm for forecasting of oil and total liquid production after hydraulic fracturing (HF). A thorough raw data study with data preprocessing algorithms was provided. The data set included 10 oil fields with more than 2,000 HF events. Each event has been characterized by well coordinates, geology, transport and storage properties, depths, and oil/liquid rates before fracturing for target and neighboring wells. Each ML model has been trained to predict monthly production rates right after fracturing and when the flows are stabilized. The gradient-boosting method justified its choice with R2 being approximately 0.7 to 0.8 on the test set for oil/total liquid production after HF. The developed ML prediction model does not require preliminary numerical simulations of a future HF design. The applied algorithm could be used as a new approach for HF candidate selection based on the real-time state of the field.

Description: Summary Fly ash, which is a pozzolan generated as a byproduct from coal-powered plants, is the most used extender in the design of lightweight cement. However, the coal-powered plants are phasing out due to global-warming concerns. There is the need to investigate other materials as substitutes to fly ash. Bentonite is a natural pozzolanic material that is abundant in nature. This pozzolanic property is enhanced upon heat treatment; however, this material has never been explored in oil-well cementing in such form. This study compares the performance of 13-ppg heated (dehydroxylated) sodium bentonite and fly-ash cement systems. The raw (commercial) sodium bentonite was dehydroxylated at 1,526°F for 3 hours. Cement slurries were prepared at 13 ppg using the heated sodium bentonite as partial replacements of cement in concentrations of 10 to 50% by weight of blend. Various tests were done at a bottomhole static temperature of 120°F, bottomhole circulating temperature of 110°F, and pressure of 1,000 psi or atmospheric pressure. All the dehydroxylated sodium bentonite systems exhibited high stability, thickening times in the range of 3 to 5 hours, and a minimum 24-hour compressive strength of 600 psi. At a concentration of 40 and 50%, the 24-hour compressive strength was approximately 800 and 787 psi, respectively. This was higher than a 13-ppg fly-ash-based cement designed at 40% cement replacement (580 psi).

Description: Summary As wells in modern operations are getting longer and more complex, assessing the effect of casing wear becomes ever more crucial. Degradation of the tubulars through mechanical wear reduces the pressure capacity significantly. In this paper, we use the finite element method (FEM) to analyze the stress distribution in degraded geometries and to assess reduction in collapse strength. A model for the collapse strength of the casing with a crescent-shaped wear groove is developed and its performance evaluated in relation to experimental data. The model was created by using the Buckingham Pi theorem to make generalized empirical expressions for yield and elastic collapse of tubulars. Finite element analysis (FEA) of 135 geometries was used in the development of the model. The results show that the generalized expressions capture the trends observed in the FEA accurately and match the experimental data from six tubular collapse tests with an average relative difference in collapse pressure of 5.2%.

Description: Summary Polymer flooding has been widely used to improve oil recovery. However, its effectiveness would be diminished when channels (e.g., fractures, fracture-like channels, void-space conduits) are present in a reservoir. In this study, we designed a series of particular sandwich-like channel models and tested the effectiveness and applicable conditions of micrometer-sized preformed particle gels (PPGs, or microgels) in improving the polymer-flooding efficiency. We studied the selective penetration and placement of the microgel particles, and their abilities for fluid diversion and oil-recovery improvement. The results suggest that polymer flooding alone would be inefficient to achieve a satisfactory oil recovery as the heterogeneity of the reservoir becomes more serious (e.g., permeability contrast kc/km 〉 50). The polymer solution would vainly flow through the channels and leave the majority of oil in the matrices behind. Additional conformance-treatment efforts are required. We tried to inject microgels in an attempt to shut off the channels. After the microgel treatment, impressive improvement of the polymer-flooding performance was observed in some of our experiments. The water cut could be reduced significantly by as high as nearly 40%, and the sweep efficiency and overall oil recovery of the polymer flood were improved. The conditions under which the microgel-treatment strategy was effective were further explored. We observed that the microgels form an external impermeable cake at the very beginning of microgel injection and prevent the gel particles from entering the matrices. Instead, the microgel particles could selectively penetrate and shut off the superpermeable channels under proper conditions. Our results suggest that the 260-µm microgel particles tested in this study are effective to attack the excessive-water-production problem and improve the oil recovery when the channel has a high permeability (〉50 darcies). The gels are unlikely to be effective for channels that are less than 30 darcies because of the penetration/transport difficulties. After the gels effectively penetrate and shut off the superpermeable channel, the subsequent polymer solution is diverted to the matrices (i.e., the unswept oil zones) to displace the bypassed oil. Overall, this study provides important insights to help achieve successful polymer-flooding applications in reservoirs with superpermeable channels.

Description: Summary A new classification of gas-hydrate deposits is proposed that takes into account their location (marine vs. permafrost), porosity type (matrix vs. fracture), and gas origin (biogenic, thermogenic, or mixed). Furthermore, by incorporating currently used Classes 1 through 4, which describe the nature of adjacent strata, a total of 16 classes of hydrate deposits have been identified. This new classification provides detailed information on the properties of the hydrate-bearing layer and adjacent strata that can be used for both scientific research and ranking of field-development potential. Using this new classification system, a qualitative ranking of field-development potential for different classes of hydrate deposits according to likely productivity, capital, and operating costs can be conducted. Finally, we demonstrate the usefulness of this new classification by applying it to 11 well-knowngas-hydrate deposits worldwide.

Description: SummaryMultiphase flowmetering is a requirement across a range of process industries, particularly those that pertain to oil and gas. Generally, both the composition and individual phase velocities are required; this results in a complex measurement task made more acute by the prevalence of turbulent flow and a variety of flow regimes. In the current review, the main technical options to meet this metrology are outlined and used to provide context for the main focus on the use of nuclear magnetic resonance (NMR) technology for multiphase flowmetering. Relevant fundamentals of NMR are detailed as is their exploitation to quantify flow composition and individual phase velocities for multiphase flow. The review then proceeds to detail three NMR multiphase flowmeter (MPFM) apparatus and concludes with a consideration of future challenges and prospects for the technology.

Description: Summarybp’s (“the company’s”) wells organization manages its operational risks through what is known as the “three lines of defense” model. This is a three-tiered approach; the first line of defense is self-verification, which wells assets apply to prevent or mitigate operational risks. The second line of defense is conducted by the safety and operational risk function using deep technical expertise. The third line of defense is provided by group audit. In this paper, we discuss the wells self-verification program evolution from its first implementation and share case studies, results, impact, lessons learned, and further steps planned as part of the continuous improvement cycle.The company’s wells organization identified nine major accident risks that have the potential to result in significant health, safety, and environment (HSE) impacts. Examples include loss of well control (LoWC), offshore vessel collision, and dropped objects. The central risk team developed bowties for these risks, with prevention barriers on cause legs and mitigation barriers on consequence legs. Detailed risk bowties are fundamental to wells self-verification, adding technical depth to allow more focused verification to be performed when compared with the original bowties, because verification is now conducted using checklists targeting barriers at their component level, defined as critical tasks and equipment. Barriers are underpinned by barrier enablers (underlying supporting systems and processes) such as control of work, safe operating limits, inspection and maintenance, etc. Checklists are standardized and are available through a single, global digital application. This permits the verifiers, typically wellsite leaders, to conduct meaningful verification conversations, record the resulting actions, track them to closure within the application, and gain a better understanding of any cumulative impacts, ineffective barriers, and areas to focus on.Self-verification results are reviewed at rig, region, wells, and upstream levels. Rigs and regions analyze barrier effectiveness and gaps and implement corrective actions with contractors at the rig or region level. Global insights are collated monthly and presented centrally to wells leadership. Common themes and valuable learnings are then addressed at the functional level, shared across the organization, or escalated by the leadership.The self-verification program at the barrier component level proved to be an effective risk management tool for the company’s wells organization. It helps to continuously identify risks, address gaps, and learn from them. Recorded assessments not only provide the wells organization with barrier performance data but also highlight opportunities to improve. Leadership uses the results from barrier verification to gain a holistic view of how major accident risks are managed. Program evolution has also eliminated duplicate reviews, improved clarity of barrier components, and improved sustainability through applying a systematic approach, standardization, digitization, and procedural discipline.

Description: Summary In this paper, we propose a hydromechanical model to simulate hydraulic fracture propagation in deep shale formations. The Drucker-Prager plasticity theory, Darcy’s law, Reynolds lubrication theory, and Kirchoff’s laws are adopted to describe the plastic deformation, matrix flow, fracture flow, and wellbore flow, respectively. A global embedded cohesive zone model is constructed to achieve the free evolution of hydraulic fractures and the characterization of natural fractures. The finite element method (FEM) and finite volume method (FVM) are used for the spatial discretization of the stress field and pressure field. On the basis of Newton-Raphson iteration, fixed-stress iteration, and Picard iteration, a mixed numerical scheme is built up to solve the strong nonlinear coupling problem. The proposed model is verified against several reference cases and experimental results. Finally, some numerical cases are carried out to investigate the influences of rock properties, natural fracture distribution, and fracturing fluid properties on the complex hydraulic fracture development. The results show that rock plasticity leads to a decrease in stimulated fracture area, an increase in average fracture width, and an increase in propagation pressure. As the cluster number increases, the adverse effect of rock plasticity on multiple hydraulic fracturing in deep shale formations increases significantly. In addition, appropriate optimization of cluster spacing could weaken the adverse effect of rock plasticity on fracturing treatment to a certain extent by using the stress interference effect.

Description: Summary This paper gives a critical review of miscibility-measurement techniques published in the open literature along with recommendations and lessons learned. Many of these published methods violate the inherent assumptions for multicontact miscibility (MCM). The confusion often arises from a failure to distinguish between first-contact miscibility (FCM), in which two fluids can be mixed in all proportions without forming two phases, and MCM, in which fluid compositions that arise during the flow of two phases in a porous medium approach a specific critical point within the constraints of the MCM definition. There are many analytical, numerical, correlational, and experimental methods available to estimate the minimum miscibility pressure (MMP) for MCM flow. The numerous available methods, some of which are quite inexpensive, have caused significant misunderstandings in the literature and in practice regarding their ability to estimate MMP. Our experience has shown that the best methods are those that honor the multicontact process (MCM), in which flow interacts with phase behavior in a prescribed way. Good methods that achieve this are slimtube experiments, detailed slimtube simulations, multiple-mixing-cell calculation methods, and the method of characteristics (MOC). Techniques such as the rising-bubble-apparatus (RBA) and vanishing-interfacial-tension (IFT) (VIT) experiments are subject to significant uncertainties, although they can still provide useful information. Numerous MMP correlations have been developed. They should be used with caution for systems similar to those used to develop the correlation. Use for other fluid systems can lead to significant errors. We discuss the advantages and disadvantages of most current methods and show that various combinations of methods can reduce uncertainty.

Description: Summary Formation damage mechanisms in general lower the quality of the near wellbore, often manifested in the form of permeability reduction, and thus reducing the productivity of production wells and injectivity of injection wells. Asphaltene deposition, as one of the important causes, can trigger serious formation damage issues and significantly restrict the production capacity of oil wells. Several mechanisms acting simultaneously contribute to the complexity associated with prediction of permeability impairment owing to asphaltene deposition; thus, integration of modeling efforts for asphaltene aggregation and deposition mechanisms seems inevitable for improved predictability. In this work, an integrated simulation approach is proposed to predict permeability impairment in porous medium. The proposed approach is novel because it integrates various mathematical models to study permeability impairment considering porosity reduction, particle aggregation, and pore connectivity loss caused by asphaltene deposition. To improve the accuracy of simulation results, porous media is considered as a bundle (different size) of capillary tubes with dynamic interconnectivity. The total volume change of interconnected tubes will directly represent permeability reduction realized in porous media. The prediction of asphaltene deposition in porous media is improved in this paper via integration of the particle aggregation model into calculation. The simulation results were verified by comparing with existing experimental data sets. After that, a sensitivity analysis was performed to study parameters that affect permeability impairment. The simulation results show that our permeability impairment model—considering asphaltene deposition, aggregation, and pore connectivity loss—can accurately reproduce the experimental results with fewer fitting or empirical parameters needed. The sensitivity analysis shows that longer aggregation time, higher flow velocity, and bigger precipitation concentration will lead to a faster permeability reduction. The findings of this study can help provide better understanding of the permeability impairment caused by asphaltene deposition and pore blockage, which provides useful insights for prediction of production performance of oil wells.

Description: Summary Shale oil formations have a very low primary yield despite advances in multistage hydraulic fracturing and horizontal drilling. In so doing, gas huff ‘n’ puff (HnP), among other improved oil recovery methods, is implemented to recover more liquid hydrocarbons. Gas HnP has proven to be an effective recovery process in shales taking into account the fracture properties, fluid-fluid interactions, and gas diffusion controlled by matrix properties. However, a laboratory-scale understanding of the gas HnP mechanism proves challenging. At this scale, measuring saturation before and after HnP tests in a nonintrusive and nondestructive manner, and understanding rock properties that affect diffusion is essential. In addition to ascertaining how the multiscale pore systems and varying mineral composition of shales affect its evaluation. The nuclear magnetic resonance (NMR) is considered a suitable tool for estimating fluid content in shales and understanding rock/fluid interactions. Generally, synthetic oil samples are used on either outcrop core plugs or crushed reservoir samples for NMR measurements, which may not be representative of rock/fluid interactions in bulk shales. This study is focused on carrying out NMR tests with dead oil on reservoir core plugs at relatively different depths to determine an effective means of saturation and understand oil production due to gas HnP. Gas HnP experiments were performed at reservoir conditions (3,500 psi and 125°C) on representative rock types from the Lower Eagle Ford (LEF) interval. Low field NMR measurements were subsequently carried out on the LEF core plugs at different states: as-received, saturated, and after gas HnP. The results show that oil recovery due to gas HnP occurred mainly in the organic pores (OPs) and inorganic pores (IPs) and ranged from 48 to 56% of the oil-in-place with indications of adsorbed/trapped methane (CH4) and remaining heavier components. This plays a vital role in evaluating the HnP process to know the extent of invasion and remaining oil components. In saturating the core plugs, the optimum saturation period was found to be 2 weeks for the LEF shale at current conditions. This presents an idea of how long to saturate a shale oil core effectively before it is tested for gas HnP. On the basis of the impact of varying mineral composition on the recovery mechanism, we observed the LEF core plug with the highest clay content to have the least recovery. This is in line with a high T1/T2 ratio alluding to reduced mobility of fluids in the presence of clay minerals with relatively small sizes of clay porosity and adsorptive surfaces.

Description: Summary The transport behaviors of both single-phase gas and single-phase water at nanoscale deviate from the predictions of continuum flow theory. The deviation is greater and more complex when both gas and liquid flow simultaneously in a pore or network of pores. We developed a pseudopotential-based lattice Boltzmann (LB) method (LBM) to simulate gas/water two-phase flow at pore scale. A key element of this LBM is the incorporation of fluid/fluid and fluid/solid interactions that successfully capture the microscopic interactions among phases. To calibrate the model, we simulated a series of simple and static nanoscale two-phase systems, including phase separation, a Laplace bubble, contact angle, and a static nanoconfined bubble. In this work, we demonstrate the use of our proposed LBM to model gas/water two-phase flow in systems like a single nanopore, two parallel nanopores, and nanoporous media. Our LBM simulations of static water-film and gas-film scenarios in nanopores agree well with the theory of disjoining pressure and serve as critical steps toward validating this approach. This work highlights the importance of interfacial forces in determining static and dynamic fluid behaviors at the nanoscale. In the Applications section, we determine the water-film thickness and disjoining pressure in a hydrophilic nanopore under the drainage process. Next, we model water imbibition into gas-filled parallel nanopores with different wettability, and simulate gas/water two-phase flow in dual-wettability nanoporous media. The results showed that isolated patches of organic matters (OMs) impede water flow, and the water relative permeability curve cuts off at water saturation [= 1–volumetric total organic carbon (TOC)]. The residual gas saturation is also controlled by the volumetric TOC, ascribed to the isolation of organic patches by the saturating water; therefore, the gas relative permeability curve cuts off at water saturation (= 1–volumetric TOC).

Description: Summary Miniaturized transponder systems are under development for monitoring unconventional reservoirs, mapping hydraulic fractures, and determining other wellbore parameters. These gadgets are an extension of radio-frequency identification (RFID) and are known as fracture robot (FracBot) nodes to recognize wireless underground sensor networks (WUSNs) for characterization and mapping of hydraulic breakages in unconventional reservoirs. 3D constellation maps of proppant-bed placement are generated by autonomous localization algorithms as FracBots are injected during hydraulic-fracturing operations. To investigate this model, a FracBot platform was established to explore this concept, and three basic functions have been explained. First, we have developed an innovative cross-layer communication model for magnetic-induction (MI) networks in altering underground environments, coupled with selections of coding, modulation, and power control and a geographic forwarding structure. Second, we have developed an innovative MI-based localization framework to capture the locations of the randomly deployed FracBot nodes by exploiting the exceptional properties of the MI field. Third, we have proposed an energy model for a linear FracBot network scheme that provides reasonable data rates while preserving collected energy limitations. Finally, to examine the functionalities of FracBot nodes in air, sand, and stone media, a physical MI-based WUSN test bed was implemented. Experiments indicated that the constructed FracBots can form a communication link and transfer data over amplitude-shift keying (ASK) modulation with 1.6 kbit/sec as a data rate and a minimum receiver sensitivity of −70 dBm. The performance of near-field-communication (NFC) antennas was affected by sand and stone media, which ultimately affect MI signal propagation and decrease the energy transfer. In sand or stone media, augmented mismatch between transmitter and receiver antennas was detected, leading to the decision that an advanced matching circuit design or an adaptive-frequency feature should be integrated into the FracBot design. This permits an optimal energy transmission and consistent communication link through sand and stone media.

Description: Summary The sudden influx of reservoir fluids (i.e., reservoir kick) into the drilling annulus is one of the common abnormal events encountered in drilling operations. A kick can lead to a blowout, causing loss of lives, assets, and damage to the environment. This study presents a framework for real-time kick monitoring and management in managed-pressure-drilling (MPD) operation. The proposed framework consists of three distinct steps: the unscented Kalman filter (UKF) is used to detect and estimate the kick's severity; the estimated kick size and optimal control theory are used to calculate the time to mitigate the kick in the best-case scenario; and on the basis of the total predicted influx and pressure rise, the monitoring system generates a warning and activates the mitigation strategy. Thus, the proposed method can estimate, monitor, and manage kick in real time, enhancing the safety and efficiency of the MPD operation. The developed method was validated and demonstrated using a simulated MPD system, a pilot-scale experimental setup, and field data collected from an MPD operation in western Canada.

Description: Summary Wells are sometimes deformed due to geomechanical shear slip, which occurs on a localized slip surface, such as a bedding plane, fault, or natural fracture. This can occur in the overburden above a conventional reservoir (during production) or within an unconventional reservoir (during completion operations). Shear slip will usually deform the casing into a recognizable shape, with lateral offset and two opposite-trending bends, and ovalized cross sections. Multifinger casing caliper tools have a recognizable response to this shape and are especially useful for diagnosing well shear. Certain other tools can also provide evidence for shear deformation. Shear deformations above a depleting, compacting reservoir are usually due to slip on bedding planes. They usually occur at multiple depths and are driven by overburden bending in response to reservoir differential compaction. Shear deformations in unconventional reservoirs, for the examples studied, have been found to be caused by slip on bedding planes and natural fractures. In both cases, models, field data, and physical reasoning suggest that slip occurs primarily due to fluid pressurization of the interface. In the case of bedding plane slip, fracturing pressure greater than the vertical stress (in regions where the vertical stress is the intermediate stress) could lead to propagation of a horizontal fracture, which then slips in shear.

Description: Summary Reaction kinetics between calcite and acid systems have been studied using the rotating disk apparatus (RDA). However, simplifying assumptions have been made to develop the current equations used to interpret RDA experiments to enable solving them analytically in contrast to using numerical methods. Previous work has revealed inadequacies of some of these assumptions, which necessitates the use of a computational fluid dynamics (CFD) model to investigate their impact on the RDA results. The objectives of the current work are to develop a calibrated CFD and proxy model to simulate the reaction in the RDA and use this model to estimate the diffusion coefficient and the reaction rate coefficient of the reaction in the RDA. The present work developed the first calibrated CFD model to determine the diffusion coefficient and the reaction rate coefficient in the RDA with minimum assumptions in the hydrochloric acid (HCl) carbonate reaction. More specifically, the model relaxes the constant fluid properties, infinite acting reactor boundaries, and constant reaction surface area assumptions. The proxy model obtained results in reduced computational time with minimal compromise on accuracy. Finally, the proposed model showed an improvement of 63% in predicting the reaction kinetics between calcite and HCl compared to traditional methods.

Description: Summary A newly formulated chemical additive from a group of amines has been tested and applied to in-situheavy oil thermal recovery. Switchable-hydrophilicity chemical additives were successfully synthesized from N,N-dimethylcyclohexylamine in the form of homogeneous and hydrophilic solution. Fundamentally, tertiary amines comprise functional groups of hydrophilic and hydrophobic components. These unique features enable this chemical additive to wet both water and heavy oil, yielding potential interfacial tension (IFT) improvement. Furthermore, the reversible chemical reaction of this chemical additive yields both positive and negative ions. An ion pair formed due to the adsorption of cations—[C8H17NH+]—on the surface of heavy oil, whereas the anions—[HCO3−]—promoted solid-phase surface charge modification, therefore, resulting in the repulsive forces between heavy oil and the rock surface—substantially improving water-wetness and restoring an irreversible wettability alteration due to the phase change phenomenon during steam injection. In this research, two types of heavy oil acquired from a field in western Alberta encompassing the viscosity of 5,616  and 46,140 cp at 25°C was utilized in each experiment. All experiments were performed and measured at high-pressure, high-temperature (HPHT) steam conditions up to 200 psi and 200°C. We perceived that favorable IFT reduction was achieved, and irreversible wettability could be restored after combining switchable-hydrophilicity tertiary amines (SHTA) with steam as a result of the solid-phase surface charge modification to be more negatively charged. Phase distribution/residual oil in the porous media developed after steam injection was able to be favorably recovered, indicating that capillary forces could be reduced. Consequently, more than 80% of the residual oil could be recuperated post-SHTA injection, presenting favorable oil recovery performance. In addition to this promising evidence, SHTA could be potentially recovered by switching its reversible chemical reaction to be in hydrophobic form, hence, promoting this chemical additive to be both reusable and more economically effective. Comprehensive studies and analyses on interfacial properties, phase distribution in porous media, and recovery performance exhibit essential points of view in further evaluating the potential of SHTA for tertiary recovery improvement. Valuable substantiations and findings provided by our research present useful information and recommendations for fields with steam injection applications.

Description: Summary Current relative permeability models rely on labeling a phase as “oil” and “gas” and cannot therefore capture accurately the effect of compositional variations on relative permeabilities and capillary pressures in enhanced oil recovery processes. Discontinuities in flux calculations caused by phase labeling problems not only cause serious convergence and stability problems but also affect the estimated recovery factor owing to incorrect phase mobilities. We developed a fully compositional simulation model using an equation of state (EoS) for relative permeabilities (kr) to eliminate the unphysical discontinuities in flux functions caused by phase labeling issues. The model can capture complex compositional and hysteresis effects for three-phase relative permeability. Each phase is modeled separately based on physical inputs that, in part, are proxies to composition. Phase flux calculations from one gridblock to another are also updated without phase labels. The tuned kr-EoS model and updated compositional simulator are demonstrated for simple ternary cases, multicycle three-phase water-alternating-gas (WAG) injection, and three-hydrocarbon-phase displacement with complex heterogeneity. The approach improves the initial estimates and convergence of flash calculations and stability analyses, as well as the convergence in the pressure solvers. The new compositional simulator allows for high-resolution simulation that gives improved accuracy in recovery estimates at significantly reduced computational time.

Description: Summary In hydrocarbon production and processing, choke and control valves mix and emulsify petroleum phases. The consequence is often that the efficiency of separation processes is affected and finally that the quality of oil and water phases is degraded. Over the last few years, low-shear valves targeting petroleum processes have emerged on the market. This paper presents four separate live-fluid experiences from low-shear valve installations, each surveyed and documented by an independent third party. Three of the installations refer to choke valves, whereas the fourth installation refers to a control valve. For each installation, standard choke and control valves were used as reference valves. In terms of downstream separation efficiency, the low-shear choke valves reduced oil-in-water concentrations respectively by 70, 45, and 60%, by total average. In the control valve application, the low-shear valve, which was located between the hydrocyclones and a compact flotation unit, reduced the oil-in-water concentration by 23%. In sum, the field installations have demonstrated that low-shear valves significantly and consistently reduce oil-in-water concentrations and thus improve the produced water quality. The results signify that low-shear valves may be used in debottlenecking separation and produced water treatment processes, reducing the environmental influence from produced water discharges. Because the low-shear technology enables processing of petroleum phases with less effort, energy, and chemicals, it also reduces emissions to air.

Description: Summary Matrix acidizing improves productivity in oil and gas wells. Hydrochloric acid (HCl), because of its many advantages such as its effectiveness, availability, and low cost, has been a typical first-choice fluid for acidizing operations. However, HCl in high-pressure/high-temperature (HP/HT) wells can be problematic because of its high reactivity, resulting in face dissolution, high corrosion rates, and high corrosion inhibition costs. Several alternatives to HCl have been tested; among them, emulsified acid is a favorable choice because of its inherent low corrosion rate, deeper penetration into the reservoir, fewer asphaltene/sludge problems, and better acid distribution due to its higher viscosity. The success of the new system is dependent upon the stability of the emulsion, especially at high temperatures. The emulsified acid must be stable until it is properly placed, and it must also be compatible with other additives in an acidizing package. This study develops a stable, emulsified acid system at 300°F using aliphatic nonionic surfactants. This paper introduces a new nonaromatic, nonionic surfactant to form an emulsified acid for HP/HT wells. The type and quality of the emulsified acid were assessed through conductivity measurements and drop tests. The thermal stability of the system was monitored as a function of time through the use of pressure tubes and a preheated oil bath at 300°F. A LUMisizer® (LUM GmbH, Berlin, Germany) and Turbiscan® (Formulaction, S. A., L’Union, France) were used to determine the stability and the average droplet size of the emulsion, respectively. The viscosity of the emulsified acid was measured at different temperatures up to 300°F as a function of shear rate (1 to 1,000 s−1). The microscopy study was used to examine the shape and the distribution of acid droplets in diesel. Coreflood studies at low and high flow rates were conducted to determine the performance of the newly developed stable emulsified acid in creating wormholes in carbonate rocks. Inductively coupled plasma and computed tomography (CT) scans were used to determine the dissolved cations and wormhole propagation, respectively. Superior stimulation results with a low pore volume of acid to breakthrough (PVBT) were achieved at 300°F with the newly developed emulsified acid system. The wormhole propagation was narrow and dominant compared to branched wormholes resulting from conventional emulsified acid systems. Results indicate that a nonionic surfactant with optimal chemistry, such as a suitable hydrophobe chain length and structure, can form a stable emulsified acid. In this study we introduce a new and effective aliphatic nonionic surfactant to create a stable emulsified acid system for matrix acidizing at HP/HT conditions, leading to a deeper penetration of acid with low pore volume to breakthrough. The successful core flood studies in the laboratory using carbonate cores suggest that the new emulsified acid system may efficiently stimulate HP/HT carbonate reservoirs.

Description: We enter the new year together as an industry and as SPE members, but also as individuals who have traveled on varying roads and experienced personal detours. While this also can be said of years past, 2020 brought with it unanticipated upheavals within the oil and gas industry, global societies, and our personal and professional lives (and livelihoods). Outlooks for the coming year depend in large part on your own worldview: Are you generally an optimist or a pessimist? An optimist believes that problems are temporary and will get better. A pessimist is convinced that the problem is here to stay and can only get worse. Optimists go into new situations with high expectations, while pessimists hang onto low expectations to prepare for negative outcomes. Do you see the glass as half full or half empty? The objective truth is that the water in the glass is at the halfway mark. The rest is up to our interpretation of that truth. The truth itself doesn’t change, but how we interpret it can have a huge effect on our actions. In his column this month, SPE President Tom Blasingame, a self-described optimist, wrote, “It’s time to look at the horizon” and “open our sails.” The metaphor describes taking action after the worst of a storm has passed or is passing and to make adjustments to get back on course. None of us are continuously optimistic or pessimistic. Life happens and moves the needle in either direction, but opening our sails may help us recover our optimism when it falters. Tapping into our resources is critical to our worldview. Am I consistently a cheerleader with a rosy view, no matter what happens? Certainly not, but I tap into my resources and mightily try to move the needle back toward optimism. (Warning: Success in doing so may not be immediate.) Your personal resources vary, and this is a reminder that SPE is one of those resources. This list is intended to serve as an “SPE Guide to Optimism.” These offerings can help to move the needle for you. Best wishes for 2021 from the JPT editorial team.

Description: Summary As drilling sections become deeper and longer, transferring more weight downhole to improve rate of penetration is the primary concern for the operator. Drillstring dynamics and buckling are some primary limiters for drilling efficiency. Aggressive drilling parameters may lead to severe downhole dynamics, which leads to cutter breakage and tool damage. When axial compression exceeds a certain threshold, the drillstring buckles sinusoidally inside the wellbore first, followed by helical buckling. Buckling leads to accelerated joint wear, tool fatigue failures, and lower drilling efficiency. To better manage drillstring dynamics and buckling, we propose a method of simulating drillstring dynamics motion and postbuckling state using an advanced transient dynamics model. An analysis methodology was developed on the basis of the finite element transient dynamics model. The model captures the enriched physics of drillstring dynamics and loading: the large deformation of buckled drillstring, the strong nonlinearity of contact and friction forces, and the dynamically triggered instability caused by drilling rotation. Transient dynamics simulations are conducted for drillstring with the actual well trajectory and rotation speed. The weight on bit (WOB) is ramped up gradually, and the drillstring deformation is monitored to detect the onset of buckling or dynamics instability. To conduct the model validation, the buckling inception loads predicted by the model are compared against the analytical equation of critical buckling loads. A field extended reach drilling (ERD) job was simulated by the model. The downhole weight and torque data from the measurement-while-drilling (MWD) tool was used to validate the weight transfer prediction by the model. Most existing buckling theories use the analytical equations of critical buckling load, which were normally derived on the basis of the idealized assumptions, such as perfect wellbore shape and uniform tubular geometry. The proposed method simulates the drillstring behaviors in the field drilling conditions and aims to capture effects of wellbore friction and string rotation. The transient dynamics model is capable of simulating drillstring dynamics movement (whirling and snaking) and weight lockup under severe helical buckling. An automatic method is proposed to interpret the drillstring behaviors from the simulation results. Using the transient dynamics model, the procedure presented in this article can simulate the dynamics and buckling behaviors of drillstring and help mitigate associated risks in well-planning and execution phases.

Description: Summary High-temperature (HT) deep carbonate reservoirs are typically drilled using barite (BaSO4) as a weighting material. Primary production in these tight reservoirs comes from the network of natural fractures, which are damaged by the invasion of mud filtrate during drilling operations. For this study, weighting material and drilling fluid were sampled at the same drillsite. X-ray diffraction (XRD) and X-ray fluorescence analyses confirmed the complex composition of the weighting material: 43.2 ± 3.8 wt% of BaSO4 and 47.8 ± 3.3 wt% of calcite (CaCO3); quartz and illite comprised the rest. The drilling fluid was used to form the filter cake in a high-pressure/high-temperature (HP/HT) filter-press apparatus at a temperature of 300°F and differential pressure of 500 psig. Compared with the weighting material, the filter cake contained less CaCO3, but more nondissolvable minerals, including quartz, illite, and kaolinite. This difference in mineral composition makes the filter cake more difficult to remove. Dissolution of laboratory-grade BaSO4, the field sample of the weighting material, and drilling-fluid filter cake were studied at 300°F and 1,000 to 1,050 psig using an autoclave equipped with a magnetic stirrer drive. Two independent techniques were used to investigate the dissolution process: analysis of the withdrawn-fluid samples using inductively coupled plasma-optical emission spectroscopy, and XRD analysis of the solid material left after the tests. The dissolution efficiency of commercial K5-diethylenetriaminepentaacetic acid (DTPA), two K4-ethylenediaminetetraacetic acid (EDTA), Na4-EDTA solutions, and two “barite dissolvers” of unknown composition was compared. K5-DTPA and K4-EDTA have similar efficiency in dissolving BaSO4 as a laboratory-grade chemical and a component of the calcite-containing weighting material. No pronounced dissolution-selectivity effect (i.e., preferential dissolution of CaCO3) was noted during the 6-hour dissolution tests with both solutions. Reported for the first time is the precipitation of barium carbonate (BaCO3) when a mixture of BaSO4 and CaCO3 is dissolved in DTPA or EDTA solutions. BaCO3 composes up to 30 wt% of the solid phase at the end of the 6-hour reaction, and can be dissolved during the field operations by 5 wt% hydrochloric acid. Being cheaper, K4-EDTA is the preferable stimulation fluid. Dilution of this chelate increases its dissolution efficiency. Compared with commonly recommended solutions of 0.5 to 0.6 M, a more dilute solution is suggested here for field application. The polymer breaker and K4-EDTA solution are incompatible; therefore, the damage should be removed in two stages if the polymer breaker is used.

Description: Summary In a recent paper, we published a machine learning method to quantitatively predict reservoir fluid gas/oil ratio (GOR) from advanced mud gas (AMG) data. The significant increase of the model accuracy compared to traditional modeling approaches makes it possible to estimate reservoir fluid GOR based on AMG data while drilling, before the wireline operation. This approach has clear advantages because of early access, low cost, and a continuous reservoir fluid GOR for all reservoir zones. This paper releases further study results to predict other reservoir fluid properties in addition to GOR, which is essential for geo-operations, field development plans, and production optimization. Two approaches were selected to predict other reservoir fluid properties. As illustrated by the reservoir fluid density example, we developed machine learning models for individual reservoir fluid properties for the first approach, similar to the GOR prediction approach in the previous paper. As for the second approach, instead of developing many machine learning models for individual reservoir fluid property, we investigated the essential properties for equation of state (EOS) fluid characterization: C6 and C7+ composition and the molecular weight and density of the C7+ fraction. Once these properties are in place, the entire spectrum of reservoir fluid properties can be calculated with the EOS model. The results of reservoir fluid property prediction are satisfactory with both approaches. The reservoir oil density prediction has a mean average error (MAE) of 0.039 g/cm3. The accuracy is similar to the typical density derived from the pressure gradient from wireline logging data. For the essential fluid properties required for EOS model prediction, the overall accuracy is less than the laboratory measurements but acceptable as the early phase estimations. The reservoir fluid properties predicted from the EOS model are similar to the predictions from individual machine learning models. We applied the field measured AMG data into the reservoir fluid property models and achieved good results, as illustrated by the reservoir fluid density example. The previous paper completed the methodology to predict all reservoir fluid properties based on AMG data. This work paves the way to generate a complete reservoir fluid log for all relevant reservoir fluid properties while drilling. The method has a significant business impact, providing full coverage of reservoir fluid properties along the well path in the early drilling phase. The advantage of providing reservoir fluid properties in all reservoir zones while drilling far outweighs the limitation of somewhat reduced reservoir fluid property accuracy.

Description: Summary Field-development optimization and optimization at the pattern scale are crucial to maximize the value of thermal enhanced-oil-recovery (EOR) projects. Application of a field net-present-value (NPV)-based pattern optimization algorithm honoring field-scale surface and subsurface constraints for in-situ-upgrading (IUP) projects has been described in the recent past. In this paper, we describe the development and application of a novel field-development-optimization capability, including the optimization of the ramp-up phase to accelerate the production to achieve a faster cash flow and high surface-facility utilization. We integrate this new capability into a robust field NPV optimization platform. A two-stagefield-development optimization algorithm is developed in this work. First, the steady-state pattern is optimized using the field-scale pattern optimization algorithm while honoring field-scale constraints and using a combined surface and subsurface performance-indicator-driven objective function. Ramp-up pattern designs are optimized separately using a solely pattern-scaleperformance-driven objective function in this stage. A preliminary pattern-delay time optimization follows next to precondition the problem for the subsequent field-scale optimization stage. The ramp-up pattern and pattern-delay times are optimized using a constant steady-state pattern in the second step of the algorithm. An appropriately penalized field-NPV-based objective function is used in this step to enforce field-scale surface and subsurface constraints. Optimization results on a realistic example application indicate that the time to oil-rate plateau could be significantly reduced on the order of multiple years while honoring the surface production constraints. This requires the use of an optimized ramp-up pattern in conjunction with the optimal steady-state pattern. The ramp-up pattern is approximately two patterns wide and features an increased heater density to deliver production acceleration. It is also notably more robust against the effects of subsurface uncertainties.

Description: Summary Most of the existing slug (SL) to churn (CH) or SL to pseudo-slug (PS) transition models (empirical and mechanistic) account for the effect of the SL liquid holdup (HLS). For simplicity, some of these models assume a constant value of HLS in SL/CH and SL/PS flow transitions, leading to a straightforward solution. Other models correlate HLS with different flow variables, resulting in an iterative solution for predicting these transitions. Using an experimental database collected from the open literature, two empirical correlations for prediction HLS at the SL/PS and SL/CH transitions (HLST) are proposed in this study. This database is composed of 1,029 data points collected in vertical, inclined, and horizontal configurations. The first correlation is developed for medium to high liquid viscosity two-phase flow (μL 〉 0.01 Pa·s), whereas the second one is developed for low liquid viscosity flow (μL ≤ 0.01 Pa·s). Both correlations are shown to be a function of superficial liquid velocity (VSL), liquid viscosity (μL), and pipe inclination angle (θ). The proposed correlations in a combination with the HLS model of Abdul-Majeed and Al-Mashat (2019) have been used to predict SL/PS and SL/CH transitions, and very satisfactory results were obtained. Furthermore, the SL/CH model of Brauner and Barnea (1986) is modified by using the proposed HLST correlations, instead of using a constant value. The modification results in a significant improvement in the prediction of SL/CH and SL/PS transitions and fixes the incorrect decrease of superficial gas velocity (VSG) with increasing VSL. The modified model follows the expected increase of VSG for high VSL, shown by the published observations. The proposed combinations are compared with the existing transition models and show superior performance among all models when tested against 357 measured data from independent studies.

Description: Summary As part of studying miscible gas injection (GI) in a major field within the Green Canyon protraction area in the Gulf of Mexico (GOM), asphaltene-formation risk was identified as a key factor affecting a potential GI project. The industry has not conducted many experiments to quantify the effect of asphaltenes on reservoir and well performance under GI conditions. In this paper we discuss a novel laboratory test for evaluating the asphaltene effect on permeability. The goals of the study were to define the asphaltene-precipitation envelope using blends of reservoir fluid and injection gas, and measure permeability reduction caused by asphaltene precipitation in a core under GI. To properly analyze the effect of GI, a suite of fluid-characterization studies was conducted, including restored-oil samples, compositional analysis, constant composition expansion (CCE), and differential vaporization. Miscibility conditions were defined through slimtube-displacement tests. Gas solubility was determined through swelling tests complemented by asphaltene-onset-pressure (AOP) testing. The unique procedure was developed to estimate the effect of asphaltene deposition on core permeability. The 1-ft-long core was saturated with the live-oil and GI mixture at a pressure greater than the AOP, and then pressure was depleted to a pressure slightly greater than the bubblepoint. Several cycles of charging and depletion were conducted to mimic continuous flow of oil along the path of injected gas and thereby to observe the accumulation of asphaltene on the rock surface. The test results indicated that during this cyclic asphaltene-deposition process, the core permeability to the live mixture decreased in the first few cycles but appeared to stabilize after Cycle 5. The deposited asphaltenes were analyzed further through environmental scanning electron microscopy (ESEM), and their deposition was confirmed by mass balance before and after the tests. Finally, a relationship was established between permeability reduction and asphaltene precipitation. The results from the asphaltene-deposition experiment show that for the sample, fluids, and conditions used, permeability is impaired as asphaltene flocculates and begins to coat the grain surfaces. This impairment reaches a plateau at approximately 40% of the initial permeability. Distribution of asphaltene along the core was measured at the end by segmenting the core and conducting solvent extraction on each segment. Our recommendation is numerical modeling of these test results and using this model to forecast the magnitude of the permeability impairment in a reservoir setting during miscible GI.

Description: Summary We are interested in the development of surrogate models for the prediction of field saturations using a fully convolutional encoder/decoder network based on the dense convolutional network (DenseNet; Huang et al. 2017), similar to the approaches used for image/image-regression tasks in deep learning. In the surrogate model, the encoder network automatically extracts the multiscale features from the raw input data, and the decoder network then uses these data to recover the input image resolution at the output of the model. The input of multiple influencing factors is considered to make our surrogate model more consistent with the physical laws, which has achieved good results in the prediction of output fields in our experiments. Various reservoir parameters including the static reservoir properties (i.e., permeability field) and dynamic reservoir properties (i.e., well placement) are used as input features, and the water-saturation distributions in different periods are taken as the output. Compared with traditional numerical reservoir simulation, which has a high computational cost and is time consuming, not only does it present the same precision, but it costs less time. At the same time, it can also be used for production optimization and history matching.

Description: Summary Imbibition is one of the most common physical phenomena in nature, and it plays an important role in enhanced oil recovery, hydrology, and environmental engineering. The imbibition in a capillary is one of the fluid transports in porous media, and the effect of a dynamic contact angle that changes with the imbibition rate on liquid-liquid imbibition is not clear. In this paper, the molecular kinetic theory (MKT) is used to study the effect of a dynamic contact angle on spontaneous capillary-liquid-liquid imbibition at a micrometer scale. The results show that: Using a scaling time, the effects of various forces in different imbibition systems can be compared, the influence of a dynamic contact angle on imbibition can be characterized by a frictional effect of the three-phase contact line, and the proposed model considering the effect of a dynamic contact angle is better than the model neglecting the effect of a dynamic contact angle. As the displacing phase viscosity increases, the influence of a dynamic contact angle on imbibition strengthens, which is attributed to a decrease in the viscous effect and an increase in the frictional effect during the imbibition process; as the displaced phase viscosity increases, the influence of a dynamic contact angle on imbibition weakens, which is attributed to an increase in the viscous effect and a decrease in the frictional effect during the imbibition process. As the interfacial tension increases, the frictional effect increases, with the result that the effect of a dynamic contact angle on imbibition increases. As the capillary becomes more hydrophilic, the effect of a dynamic contact angle on imbibition becomes stronger because of a decreasing viscous effect and an increasing frictional effect. As the capillary length increases, the viscous effect increases, whereas the frictional effect decreases, leading to a decrease in the dynamic contact angle effect. As the capillary radius increases, the frictional force decreases, whereas its proportion in total resistance or the frictional effect increases, resulting in an increase in the effect of a dynamic contact angle. This work sheds light on the effect of a dynamic contact angle on capillary-liquid-liquid imbibition, including displacing phase viscosity, displaced phase viscosity, interfacial tension, capillary wettability, length, and radius. It will provide new insights into manipulating a capillary imbibition process and provide a fundamental theory for enhanced oil recovery by imbibition in conventional or unconventional reservoirs. Supplementary materials are available in support of this paper and have been published online under Supplementary Data at https://doi.org/10.2118/205490-PA. SPE is not responsible for the content or functionality of supplementary materials supplied by the authors.

Description: Summary Pulsating hydraulic fracturing (PHF) is a promising fracturing technology for unconventional reservoirs because it could improve the hydraulic fracturing efficiency through inducing the fatigue failure of reservoir rocks. Understanding of the pressure wave propagation behavior in wellbores and fractures plays an important role in PHF optimization. In this paper, a transient flow model (TFM) was used to describe the physical process of pressure wave propagation induced by PHF, and this model was solved by the method of characteristics (MOC). Combination of the TFM and MOC was validated with experimental data. The impacts of controlling factors on the pressure wave propagation behavior were fully discussed, and these factors include the frequency of input loading, an injection mode, an injection position, and friction. More than 10,000 sets of pressure wave propagation behaviors in different scenarios were simulated, and their differences were illustrated. In addition, the generation mechanisms of different pressure wave propagation behaviors were explained by the Fourier transform theory and the vibration theory. The important finding is that there is resonance phenomenon in the propagation of the pressure wave, and the resonance frequencies are almost equal to the natural frequencies of a fluid column. As a consequence of resonance phenomenon, the amplitudes of bottomhole pressure (BHP) and fracture tip pressure will increase sharply when the input loading frequency is close to the resonance frequency and less than 5 Hz; otherwise, the resonance phenomenon will disappear. Furthermore, an injection mode can alter the resonance frequency and the amplitude and frequency of the induced pressure wave. In addition, a friction effect can significantly decrease both the resonance frequency and the resonance amplitude. These findings indicate that the optimized input loading frequency should be close to the natural frequency of a fracturing fluid in a wellbore to enhance its BHP.

Description: Summary Hole cleaning is a concern in directional and horizontal well drilling operations where drill cuttings tend to settle in the lower annulus section. Laboratory-scale experiments were performed with different non-Newtonian fluids in a 6.16-m-long, 114.3- × 63.5-mm transparent annulus test section to investigate cuttings transport behavior. This experimental study focused on understanding the cuttings transport mechanism in the annulus section with high-speed imaging technology. The movement of cuttings in the inclined annular section was captured with a high-speed camera at 2,000 frames/sec. Also, cuttings bed movement patterns at different fluid velocities and inner pipe rotations were captured with a digital single-lens reflex video camera. The electrical resistance tomography (ERT) system was used to quantify the cuttings volume fraction in the annulus. Different solid bed heights and cuttings movements were observed based on fluid rheology, fluid velocity, and inner pipe rotation. The mechanistic three-layer cuttings transport model was visualized with the experimental procedure. This study showed that solid bed height is significantly reduced with an increase in the inner pipe rotation. This study also identified that cuttings bed thickness largely depends on fluid rheology and wellbore inclination. The image from the high-speed camera identified a downward trend of some rolling particles in the annulus caused by gravitational force at a low mud velocity. Visual observation from a high-speed camera identified a helical motion of solid particles when the drillpipe is in contact with solid particles and rotating at a higher rev/min. Different cuttings movement patterns such as: rolling, sliding, suspension, helical movement, and downward movement were identified from the visualization of a high-speedcamera.

Description: Summary Temperature-dependent irreducible water saturation has great implications for heavy-oil production. Especially in processes using thermal methods, the irreducible water saturation varies significantly when temperature rises from the initial reservoir condition to the temperature of injected hot fluids. In this work, the irreducible water saturation retained in a heavy-oil/oil-sands reservoir has been theoretically analyzed as a function of temperature in the view of thermodynamics. This analysis involves oil/water interactions, thermodynamic stability, pendular rings between particles, and a dense random-packing theory. The temperature-dependent irreducible water saturation in two heavy-oil reservoir samples (Coalinga and Huntington Beach) and two oil-sands reservoir samples (Cat Canyon and Peace River) have been analyzed using an oil/water/silica system. The computed results have been compared with published experimental data. The good agreements of the comparison demonstrate the feasibility of the proposed analysis to describe the irreducible water saturation in a heavy-oil/oil-sands reservoir up to 300°C. Through these analyses, the theoretical understandings of temperature-dependent irreducible water in a heavy-oil/oil-sands reservoir have been provided. As temperature increases, the mutual water/oil solubilities are increased by enhanced molecular interactions, as well as the surface energy at an oil/water connecting interface. As a result, the oil/water interfacial tension (IFT) decreases, which diminishes the contact angle and enlarges a water-filled pendular ring between particles at elevated temperatures. Thus, the irreducible water saturation is increased by the enlarged pendular rings in a dense packing porous medium. In addition, this study demonstrates the possibilities to alter the irreducible water saturation appropriately in a heavy-oil/oil-sands reservoir to enhance oil recovery, decrease water cut, save costs of surface oil/water separation, and reduce heat consumption.

Description: Summary In this paper, we introduce the entropy weight method (EWM) to establish a comprehensive evaluation model able to quantify the brittleness of reservoir rocks. Based on the evaluation model and using the adaptive finite element-discrete element (FE-DE) method, a 3D model is established to simulate and compare the propagation behavior of hydraulic fractures in different brittle and ductile reservoirs. A failure criterion combining the Mohr-Coulomb strength criterion and the Rankine tensile criterion is used to characterize the softening and yielding behavior of the fracture tip and the shear plastic failure behavior away from the crack tip during the propagation of a fracture. To understand the effects of rock brittleness and ductility on hydraulic fracture propagation more intuitively, two groups of ideal cases with a single failure mode are designed, and the fracture propagation characteristics are compared and analyzed. By combining natural rock core scenarios with single failure mode cases, a comprehensive evaluation index BIf for reservoir brittleness and ductility is constructed. The simulation experiment results indicate that fractures in brittle reservoirs tended to form a complex network. With enhanced ductility, the yielding and softening of reservoirs hamper fracture propagation, leading to the formation of a simple network, smaller fracture area (FA), larger fracture volume, and the need for higher initiation pressure. The comprehensive index BIf can be used to define brittleness or ductility as the dominant factor of fracturing behavior. That is, 0 

Description: Summary The highly exothermic reaction between ammonium chloride (NH4Cl) and sodium nitrite (NaNO2) has an important application in the area of flow assurance. Because of the high heat generation, this reaction has been used as a heat source for the fluidization of low-melting-point deposits formed during oil and gas production. Because this reaction is strongly pH dependent, the incorrect choice of pH can result in an uncontrollable temperature increase caused by the system’s inability to dissipate the large amount of heat generated in a short time, causing accidents such as structural damage and explosions. Thus, the aim of this work was to study a method that involved adjusting the pH over time to ensure controlled heat generation, with high calorimetric conversion, and avoid the development of a thermal-runaway reaction (pH-based control of the kinetics and process safety). The kinetics and thermodynamics of this reaction were studied using heat-flow reaction calorimetry and attenuated total reflection (ATR)-Fourier-transform infrared (FTIR) (ATR-FTIR) spectroscopy. Following a semiempirical approach, calorimetric and spectroscopic data were fitted to a kinetic equation using nitrite, ammonium (NH4+), and hydronium concentrations. The molar enthalpy calculated was –322.92 kJ/mol, and the Arrhenius parameters were determined as the frequency factor [ln(A)] = 22.21 and the apparent activation energy (Ea) = 63.40  kJ/mol. The kinetic model constructed made it possible to properly evaluate the pH profile that should be maintained to control the kinetics (heat-generation rate) and process safety [time to maximum rate under adiabatic conditions (TMRAD)] of the reaction. The strategy of adjusting the pH over time ensured controlled heat generation and high calorimetric conversion, which cannot be achieved by simply adding catalyst at the beginning of the reaction, and minimized the risk of developing a runaway reaction. However, in real applications, the pH control must be made using the balance between the thermal risk (TMRAD) and the performance of the method (qr), because although it is possible to decrease the thermal risk (increase the value of TMRAD) by increasing the pH, this increase is accompanied by a decrease in the heat-generation rate. Thus, from the proper balance of these factors (qr and TMRAD), pH control can ensure adequate levels of heat production within an acceptable thermal risk. Supplementary materials are available in support of this paper and have been published online under Supplementary Data at https://doi.org/10.2118/205389-PA. SPE is not responsible for the content or functionality of supplementary materials supplied by the authors.

Description: Summary The ability of geochemistry techniques in reservoir-continuity studies has already been proved. Most of the traditional methods mainly involve analyzing nonpolar components of crude oil and overlooking polar components. Despite valuable information obtained from nonpolar components, these compounds are sometimes affected by various alterations or likely provide only a piece of the reservoir-compartmentalization puzzle. In this paper, an integrated geochemical approach that uses nonpolar (i.e., saturates and aromatics) and polar (i.e., asphaltenes) components of crude oil was performed to evaluate reservoir continuity efficiently. The Shadegan Oil Field in the Dezful Embayment in southwest Iran was investigated for reservoir-continuity studies to show the efficiency of this proposed technique. The selected interparaffin peak ratios and light hydrocarbons [the C7 oil correlation star diagram (C7CSD)] from whole-oil gas chromatography (GC) (WOGC) chromatograms were used to obtain oil fingerprints from the nonpolar fraction of crude oils. The Fourier-transform infrared (FTIR) spectroscopy of asphaltenes was applied to obtain oil fingerprints from the polar fraction of crude oils. The pairwise comparison of studied wells by each technique was summarized in a similarity matrix with green, yellow, and red colors to show connectivity, limited connectivity, and disconnectivity according to oil fingerprints. Finally, a compartmentalization model was prepared from the integrated results of different techniques considering the worst-case scenarios regarding the occurrence or absence of reservoir continuity when relying on individual methods for the studied field. Results show that the Shadegan Oil Field comprises three zones in the Asmari Reservoir and two zones in the Bangestan Reservoir. Reservoir-engineering data, including pressure data and pressure/volume/temperature (PVT), completely corroborated the obtained results from the geochemical approach. The consistency of results suggested FTIR oil fingerprinting of asphaltene as a novel and straightforward technique, which is a complementary or even alternative method with respect to previous geochemical methods.

Description: Summary Two-phase flow is a common occurrence in pipes of oil and gas developments. Current predictive tools are based on the mechanistic two-fluid model, which requires the use of closure relations to predict integral flow parameters such as liquid holdup (or void fraction) and pressure gradient. However, these closure relations carry the highest uncertainties in the model. In particular, significant discrepancies have been found between experimental data and closure relations for the Taylor bubble velocity in slug flow, which has been determined to strongly affect the mechanistic model predictions (Lizarraga-García 2016). In this work, we study the behavior of Taylor bubbles in vertical and inclined pipes with upward and downward flow using a validated 3D computational fluid dynamics (CFD) approach with level set method implemented in a commercial code. A total of 56 cases are simulated, covering a wide range of fluid properties, pipe diameters, and inclination angles: Eo ∈ [10, 700]; Mo ∈ [1×10–6, 5×103]; ReSL ∈ [–40, 10]; θ ∈ [5°, 90°]. For bubbles in vertical upward flows, the simulated distribution parameter, C0, is successfully compared with an existing model. However, the C0 values of downward and inclined slug flows where the bubble becomes asymmetric are shown to be significantly different from their respective vertical upward flow values, and no current model exists for the fluids simulated here. The main contributions of this work are (1) the relatively large 3D numerical database generated for this type of flow, (2) the study of the asymmetric nature of inclined and some vertical downward slug flows, and (3) the analysis of its impact on the distribution parameter, C0.

Description: Summary Accurately predicting wax deposits in a crude pipeline through empirical formulas or numerical modeling is unreliable because of the incomplete mechanism and the time-dependent unsteady actual operating conditions. With the help of the data collected by the supervisory control and data acquisition system of pipelines, wax deposit prediction is made possible by developing the time-dependent data mining method. In this article, the data from a typical long-distance crude pipeline in China operating over a 4-year time period was investigated. The inlet temperature prediction was first conducted by developing the long short-term memory (LSTM)-recurrent neural networks (RNNs) model, during which the feature sequencing, overfitting problems, and optimal hyperparameters were fully considered. Because of the time sequence cell, the accuracy of the LSTM-RNN model, as well as the time consumption, is much better than the RNN model when dealing with a great deal of data over a long period of time. Taking the inlet temperature prediction results as input features, the prediction model of average wax deposit thickness was established based on the backpropagation (BP) neural network and optimized by the particle swarm optimization (PSO), chaos particle swarm optimization (CPSO), and adaptive chaos particle swarm optimization (ACPSO) algorithms. The conclusions and associated algorithm from this article help to determine the reasonable pigging circle of long-distance pipelines practically. It could also be applied to guide the wax deposit prediction in the wellbore or oil-gatheringpipes.

Description: Summary The historical challenges and high failure rate of using standalone screen in cased and perforated wellbores pushed several operators to consider cased-hole gravel packing or frac packing as the preferred completion. Despite the reliability of these options, they are more expensive than a standalone screen completion. In this paper, we employ a combined physical laboratory testing and computational fluid dynamics (CFD) for laboratory scale and field scale to assess the potential use of the standalone screen in completing the cased and perforated wells. The aim is to design a fit-to-purpose sand control method in cased and perforated wells and provide guidelines in perforation strategy and investigate screen and perforation characteristics. More specifically, the simultaneous effect of screen and perforation parameters, near wellbore conditions on pressure distribution and pressure drop are investigated in detail. A common mistake in completion operation is to separately focus on the design of the screen based on the reservoir sand print and design of the perforation. If sand control is deemed to be required, the perforation strategy and design must go hand in hand with sand control design. Several experiments and simulation models were designed to better understand the effect of perforation density, the fill-up of the annular gap between the casing and screen, perforation collapse, and formation and perforation damage on pressure drop. The experiments consisted of a series of step-rate tests to investigate the role of fluid rate on pressure drop and sand production. There is a critical rate at which the sand filling up the annular gap will fluidize. Both test results and CFD simulation scenarios are comparatively capable to establish the relation between wellbore pressure drop and perforation parameters and determine the optimized design. The results of this study highlight the workflow to optimize the standalone screen design for the application in cased and perforated completions. The proper design of standalone screen and perforation parameters allows maintaining cost-effective well productivity. Results of this work could be used for choosing the proper sand control and perforation strategy.

Description: Summary Decline curve analysis (DCA) has been the mainstay in unconventional reservoir evaluation. Because of the extremely low matrix permeability, each well is evaluated economically for ultimate recovery as if it were its own reservoir. Classification and normalization of well potential is difficult because of ever-changing stimulation total contact area and a hyperbolic curve fit parameter that is disconnected from any traditional reservoir characterization descriptor. A new discrete fracture model approach allows direct modeling of inflow performance in terms of fracture geometry, drainage volume shape, and matrix permeability. Running such a model with variable geometrical input to match the data in lieu of standard regression techniques allows extraction of a meaningful parameter set for reservoir characterization, an expected outcome from all conventional well testing. Because the entirety of unconventional well operation is in transient mode, the discrete fractured well solution to the diffusivity equation is used to model temporal well performance. The analytical solution to the diffusivity equation for a line source or a 2D fracture operating under constrained bottomhole pressure consists of a sum of terms, each with exponential damping with time. Each of these terms has a relationship with the constant rate, semisteady-state solution for inflow, although the well is not operated with constant rate, nor will this flow regime ever be realized. The new model is compared with known literature models, and sensitivity analyses are presented for variable geometry to illustrate the depiction of different time regimes naturally falling out of the unified diffusivity equation solution for discrete fractures. We demonstrate that apparent hyperbolic character transitioning to exponential decline can be modeled directly with this new methodology without the need to define any crossover point. The mathematical solution to the physical problem captures the rate transient functionality and any and all transitions. Each exponential term in the model is related to the various possible interferences that may develop, each occurring at a different time, thus yielding geometrical information about the drainage pattern or development of fracture interference within the context of ultralow matrix permeability. Previous results analyzed by traditional DCA can be reinterpreted with this model to yield an alternate set of descriptors. The approach can be used to characterize the efficacy of evolving stimulation practices in terms of geometry within the same field and thus contribute to the current type curve analyses subject to binning. It enables the possibility of intermixing of vertical and horizontal well performance information as simply gathering systems of different geometry operating in the same reservoir. The new method will assist in reservoir characterization and evaluation of evolving stimulation technologies in the same field and allow classification of new type curves.

Description: Summary Low-frequency distributed-acoustic-sensing (LF-DAS) strain data are direct measurements of in-situ rock deformation during hydraulic-fracturing treatments. In addition to monitoring fracture propagation and identifying fracture hits, quantitative strain measurements of LF-DAS provide opportunities to quantify fracture geometries. Recently, we proposed a Green’s function–based algorithm for the inversion of LF-DAS strain data (Liu et al. 2020b) that shows an accurate estimation of fracture width near the monitor well with single-cluster completions. However, multicluster completions with tighter cluster spacings are more commonly adopted in recent completion designs. One main challenge in the inversion of LF-DAS strain data under such circumstances is that strain measurements at fracture-hit locations by LF-DAS are not reliable, which makes the individual contribution of each fracture to the measured strain data indistinguishable. In this study, we first extended the inversion algorithm to handle multiple fractures, investigated the uncertainties of the inversion results, and proposed possible mitigation to the challenges raised by completion designs and field data acquisition through a synthetic case study. Ideally, there are available data on both sides of each fracture so that the inverted width of each fracture can be obtained with a negligible error. In reality, the strain data are usually limited, providing less constraint on the width of individual fracture. Nevertheless, the inversion results provide an accurate estimation of the width summation of all fractures. To evaluate the individual fracture width, a time-dependent constraint is added to the inversion algorithm. We assume that the width at the current timestep is dependent on the width at the previous step and the width variation between the two timesteps. The width variation can be roughly estimated from LF-DASstrain-rate data at the fracture-hit location. This extra constraint helps to improve the inversion performance. Finally, a field example is presented. We show the width summation of all fractures and the width of each individual fracture as a function of treatment time. The time-dependent width profiles show consistent trends with the LF-DASstrain-rate data. The calculated strains from the inverted model match well with the LF-DAS measured strain data. The findings demonstrate the potential of LF-DAS data for quantitative hydraulic-fracture characterization and provide insights on better use of LF-DAS data. The direct information on fracture width helps to calibrate fracturing models and optimize the completion designs.

Description: Summary The surface dynamometer card is composed of ground load and ground displacement, which is of great significance to reflect the operation of rod pumping and the exploitation of crude oil. However, the current method of obtaining the surface dynamometer by sensors is a huge financial investment on the sensor installations and maintenance. In this paper, we propose an innovative method based on deep learning to reproduce the surface dynamometer card directly from electrical parameters. In our method, the convolution neural network is used as the basic layer to automatically extract the spatial characteristics of input data. A long short-term memory (LSTM) network as the core component is used for the output layer to consider the time dependence of the dynamometer card. Finally, the experimental shows that the proposed method achieves the mean relative error (MRE) of 4.00% on the real oil well data in A-oilfield, and the dynamometer card calculated by our model is basically consistent with the field data. In addition, the method has been tested in new wells with a rod pumping system, and the results show that the accuracy of the model is close to 90%, which has already greatly outperformed the previous methods.

Description: Summary Two of the most important parameters to monitor during a primary cementing job are the flow rate in and return flow rate measurements. To achieve optimum job results of a primary cementing job, measuring annular return rates and comparing them with simulated data in real time will provide a better understanding of job signatures and result in the best possible top of cement (TOC) estimation prior to running any cement evaluation log or making a decision to continue drilling the next section of the well. The return rate job signature along with the wellhead pressure is essential to understanding the behavior and discrepancies between simulated and acquired surface data. Therefore, to assess the risk of job issues, such as unsuspected washout and lost circulation among others, accurate measurements of the return rate are critical. Historically, the cement job evaluation has been limited by the fact that most drilling rigs do not have an accurate flowmeter installed on the annulus return line, and a simple verification of mud tanks volume vs. pumped volume, as reported by drillers or mud loggers, more often than not results in an unreliable assessment of the volume lost downhole, due to the unfamiliarity with the U-tubing effect and lack of data consolidation from the cement unit (flow rate in) and the rig (flow rate in and flow rate out). In this paper, we will review a solution developed to mitigate the lack of a direct flow-rate measurement by computing and displaying the return rate using either a paddle meter measurement or the derivative over time of the volume observed in the rig tanks.

Description: Summary An approach to the analysis of production data from waterflooded oil fields is proposed in this paper. The method builds on the established techniques of rate-transient analysis (RTA) and extends the analysis period to include the transient- and steady-state effects caused by a water-injection well. This includes the initial rate transient during primary production, the depletion period of boundary-dominated flow (BDF), a transient period after injection starts and diffuses across the reservoir, and the steady-state production that follows. RTA will be applied to immiscible displacement using a graph that can be used to ascertain reservoir properties and evaluate performance aspects of the waterflood. The developed solutions can also be used for accurate and rapid forecasting of all production transience and boundary-dominated behavior at all stages of field life. Rigorous solutions are derived for the transient unit mobility displacement of a reservoir fluid, and for both constant-rate-injection and constant-pressure-injection after a period of reservoir depletion. A simple treatment of two-phase flow is given to extend this to the water/oil-displacement problem. The solutions are analytical and are validated using reservoir simulation and applied to field cases. Individual wells or total fields can be studied with this technique; several examples of both will be given. Practical cases are given for use of the new theory. The equations can be applied to production-data interpretation, production forecasting, injection-water allocation, and for the diagnosis of waterflood-performanceproblems. Correction Note: The y-axis of Fig. 8d was corrected to "Dimensionless Decline Rate Integral, qDdi". No other content was changed.

Description: Summary In this study, an investigation of sand transport in heavy-oil/water multiphase flow is performed. The study is conducted in three multiphase-flowpipeline-test facilities with internal diameters (IDs) of 1, 1, and 3 in. The pipeline orientations relative to the horizontal in the facilities are 0, +30, and 0°, respectively. Oil viscosity of 3.5 and 10.0 Pa·s with sand volume fractions from 0.010 to 0.100 vol% were used in the study. The effects of oil viscosity, upward inclination, sand volume fraction, pipe ID, and water cut on the sand-transport mechanism in pipelines are investigated. In the horizontal test section, flow patterns—namely, dispersed flow (DF), plug flow (PF), plug flow with moving sand bed (PFM), and plug flow with stationary sand bed (PFS)—were identified through flow visualization. In addition to the aforementioned, two flow patterns—stratified wavy flow with moving sand bed (SWM) and stratified wavy flow with dunes (SWD)—were observed in the inclined pipeline orientation. The pressure gradient measured decreased with a decrease in water cut until a minimum value was reached. Beyond the minimum pressure gradient, further reduction in water cut led to an increase in pressure gradient. The sand minimum transport condition (MTC) in the oil/water/sand test was largely the same for the 1-in. 30° upward inclined and the 1-in. horizontal test section. In contrast, that of the 3-in. horizontal test section was considerably higher. An improved MTC predictive correlation is proposed for multiphase heavy-oil/water/sand flow. The proposed correlation outperforms the existing models when tested on the heavy-oil/water/sand data set.

Description: Summary As an economical and efficient artificial lift method, plunger lift can be used to unload the accumulated liquids from the bottom of gas wells, which helps lower the bottomhole pressure, resulting in higher gas production rate. However, the transient flow behavior of the plunger-lift-aided production system is still not well understood due to the lack of a reliable and accurate prediction model. In this study, a transient mechanistic model is developed to simulate the comprehensive dynamic process of a plunger-lift system that is cyclically paced by a surface control valve. Starting from the Gasbarri and Wiggins (2001) dynamic plunger-lift model, four stages in the cyclic movement of a plunger can be identified and calculated using a set of specific governing equations. Considering the gas flows with a plunger in the tubing, the model can calculate the instant velocities of the plunger during its rising and falling movement. The classical inflow performance relationship (IPR) is employed as the reservoir model to obtain the fluid flow rates from the reservoir to the wellbore. The proposed new model can capture the essential parameters of plunger-lift cycles, including plunger velocity/acceleration, tubing/casing pressure, production rates, etc. Compared to previous models, the predicted rising and falling speeds of the plunger are improved. The hydrocarbon mixture properties in the gas well are computed by a compositional model in this study, which provides more accurate and reasonable predictions of tubing and casing pressure. Several parametric studies are presented in the paper. These studies will help to understand the influence of different parameters on the process of plunger lift. An appropriate combination of casing and tubing pressure should be taken into consideration. A reservoir coefficient term is introduced and defined. A larger reservoir coefficient will improve the ultimate profitability of the well by increasing the production rate at the beginning and accelerate the depletion of gas wells. If the gas/liquid ratio (GLR) is too low, liquid loading may be triggered. The parametric study shows that an adequate GLR is necessary for reliable plunger-liftperformance.

Description: Summary With the flourishing shale gas exploitation producing more oil-based mud (OBM) cuttings, the hard-to-treat hazardous wastes heavily burden the local environment. However, the problems of treating OBM cuttings, such as huge energy consumption, tremendous treatment costs, and high risk of secondary contamination, still remain unsolved with the current treatment technologies, such as thermal desorption, incineration, and chemical extraction. In this study, we introduce a new method and equipment based on cyclone desorption to recover oil from OBM cuttings. The technological process includes viscosity reduction in heated gas, cyclone deoiling, condensation and recycling of the exhaust, and separation of oil and water in the coalescer. Based on the analysis of the physicochemical properties and the oil distribution inside the OBM cuttings samples collected from the Chongqing shale gas field, we designed this cyclone oil desorption technology and built the pilot-scale equipment to conduct the deoiling experiments. The results showed that the deoiling efficiency of OBM cuttings improved as the processing time increased. To be precise, after 2.7 seconds of treatment, the oil content of the cuttings samples fell sharply from 17.9 to 0.16%, which is about one-half of the maximum allowable oil content in pollutants of 0.3%, specified in the national standard (GB 4284-84 1985) promulgated by the People’s Republic of China. The foundation of the technology is that the particles have a high-speedself-rotation (more than 30,000 rad/s) coupled with a revolution in the cyclone in which a generated centrifugal force removes the oil from the pores of the particles. This process is purely physical and involves no phase change of the oil, so it is free of chemical addition and high heating temperature. The application of this newly developed cyclone oil desorption technology is expected to lower the treatment costs, enhance the processing efficiency, contribute to the energy development, and eventually benefit the local environment where the shale gas exploitations take place.

Description: Summary The accumulation of rock cuttings, proppant, and other solid debris in the wellbore caused by inadequate cleanout remarkably impedes field operations. The cuttings removal process becomes a more challenging task as the coiled-tubing techniques are used during drilling and fracturing operations. This article presents a new hole cleaning model, which calculates the critical transport velocity (CTV) in conventional and fibrous water-based fluids. The study is aimed to establish an accurate mechanistic model for optimizing wellbore cleanout in horizontal and inclined wells. The new CTV model is established to predict the initiation of bed particle movement during cleanout operations. The model is formulated considering the impact of fiber using a special drag coefficient (i.e., fiber drag coefficient), which represents the mechanical and hydrodynamic actions of suspended fiber particles and their network. The dominant forces acting on a single bed particle are considered to develop the model. Furthermore, to enhance the precision of the model, recently developed hydraulic correlations are used to compute the average bed shear stress, which is required to determine the CTV. In horizontal and highly deviated wells, the wellbore geometry is often eccentric, resulting in the formation of flow stagnant zones that are difficult to clean. The bed shear stress in these zones is sensitive to the bed thickness. The existing wellbore cleanout models do not account for the variation in bed shear stress. Thus, their accuracy is limited when stagnant zones are formed. The new model addresses this problem by incorporating hydraulic correlations to account for bed shear stress variation with bed height. The accuracy of the new model is validated with published measurements and compared with the precision of an existing model. The use of fiber drag and bed shear stress correlations has improved model accuracy and aided in capturing the contribution of fiber in improving wellbore cleanout. As a result, for fibrous and conventional water-based fluids, the predictions of the new model have demonstrated good agreement with experimental measurements and provided better predictions than the existing model. Model predictions show a noticeable reduction in fluid circulation rate caused by the addition of a small quantity of fiber (0.04% w/w) in the fluid. In addition, results show that the existing model overpredicts the cleaning performance of both conventional and fibrous water-basedmuds.

Description: Summary To determine which salt-based cement system (potassium chloride or sodium chloride) was suitable for cementing across halite and anhydrite salt sections in West Africa, eight slurry recipes were tested to assess how formation salt contamination would affect slurry properties. The formation salt used for testing was sampled from a deepwater, presalt well in Angola. The recommendations developed from the laboratory study were implemented in 10 projects across West Africa over 5 years with 100% operational and well integrity success. A candidate deepwater well was selected in which the surface and intermediate strings penetrated salt formations. Four slurry designs (a lead and tail slurry used on each casing string) were programmed. Each slurry was designed and tested as two distinct systems using potassium chloride and sodium chloride salt, respectively, yielding a total of eight slurry designs. Using the methodology and data presented by Martins et al. (2002), the mass of dissolved formation salt that each slurry may receive during placement was estimated and duly incorporated into each slurry design. Subsequently, the salt-contaminated slurries were tested and compared with the properties of the initial uncontaminated slurries. On the basis of these results, conclusions were then made on which salt slurry system (potassium chloride or sodium chloride) exhibited better liquid and set properties after contamination with formation salt. Subsequently, this knowledge was applied to 10 projects across three countries in West Africa. This study showed that when the contact time of liquid cement slurry to salt formation was low—typically when the salt-formation interval across which the cement slurry flowed was less than 100 m thick—the level of formation salt dissolution entering the slurry during placement was limited. In this case, a potassium chloride salt-based slurry delivered improved liquid and set properties as compared with a sodium chloride salt-based slurry. In the field, this knowledge was applied in all oilfield projects cemented by an oilfield service company between 2015 and 2020. This included deepwater, shallow offshore, and onshore wells. All related salt-zone cement jobs, including sidetrack plugs, placed across the salt formations were successful on the first attempt. In an absence of industry consensus around salt-formation cement slurry design, this paper validates a guideline for West Africa, based on results from laboratory testing and 5 years of field application. In contrast to current literature that recommends only sodium chloride salt-based slurry designs across halite or anhydrite salt intervals, this work demonstrates that potassium chloride salt-based slurry systems can effectively be used to achieve well integrity where a halite or anhydrite salt interval is less than 100 m (328.1 ft) thick.

Description: Summary The subsurface safety valve (SSV) is an essential device of a subsea production system. The hydraulic SSV is widely used, but it has some drawbacks when used in deep water: It needs a long pipeline from platform to wellhead, and brings additional pressure loss. The system needs higher pressure and manufacturing costs, but the response is slower. To solve the problem, a new all-electricsurface-controlled SSV (E-SCSSV) system consisting of an E-SCSSV body and a control system is developed. The innovative structural designs include electric-drive mechanisms and a magnetic coupler that can transfer linear motion. The mechanical property of E-SCSSV body is analyzed to determine the ability to resist well pressure. The coupling rule of the magnetic coupler is studied through experiment and finite-element analysis. The dynamics of the failure-safety mechanism is investigated to obtain the optimal performance and combination of structural parameters. Using sensors on the E-SCSSV body, the operation status of the E-SCSSV can be monitored. An E-SCSSV system prototype has been developed to validate the function of the E-SCSSV. Reliability requirements are discussed. Compared with existing SSVs, the biggest advantage of the E-SCSSV is its quick response.

Description: Summary Hydraulic fracture (HF) modeling is a multiscale and multiphysics problem. It should capture various effects, including those of in-situ stresses, poroelasticity, and reservoir heterogeneities at different length scales. A peridynamics (PD)-based hydraulic fracturing simulator has been demonstrated to reproduce this physics accurately. However, accounting for such details leads to a reduction in computational speed. In this paper, we present a novel coupling of the PD-based simulator with numerically efficient finite element methods (FEMs) and finite volume methods (FVMs) to achieve a significant improvement in computational performance. Unlike classical methods, such as FEM and FVM that solve differential equations, PD uses an integral formulation to circumvent the undefined spatial derivatives at crack tips. We implemented four novel coupling schemes of our PD-based simulator with FEM and FVM: static PD region scheme, dynamic PD region scheme, adaptive mesh refinement scheme, and dynamic mesh coarsening scheme. PD equations are solved using a refined mesh close to the fracture, whereas FE/FV equations are solved using a progressively coarser mesh away from the fracture. As the fracture grows, a dynamic conversion of FE/FV cells to PD nodes and adaptive mesh refinement are incorporated. To improve the performance further, the dynamic mesh coarsening scheme additionally converts the fine PD nodes back to coarse FE/FV cells as the HF grows in length. The coupling schemes are verified against the Kristianovich-Geertsma-de Klerk (KGD) fracture propagation problem. No spurious behavior is observed near the transition between PD and FE/FV regions. In the first three coupling schemes, the computational runtime for single fracture propagation is reduced by up to 10, 20, and 50 times, respectively, compared to a pure PD model. Laboratory experiments on the interaction of an HF with a natural fracture (NF) are revisited. The model captures complex fracture behavior, such as turning in the case of low stress contrast and low angle of interaction, kinking for higher stress contrast or higher angle of interaction, and fracture crossing for near-orthogonal NFs. Moreover, several previously reported phenomena, including fracture propagation at an angle to the principal stress directions, competing fracture growth from multiple closely spaced clusters, and interaction with layers of varying mechanical properties are successfully modeled. Thus, the coupling of PD with FEM and FVM offers an innovative and fundamentally comprehensive solution to alleviate the high computational costs typically associated with the pure PD-based hydraulic fracturing simulations. At the same time, these coupling schemes retain the versatility of the nonlocal PD formulation at modeling the evolution of arbitrary material damage, commonly observed during HF propagation in complex heterogeneous reservoirs.

Description: Summary The extensive depletion of the development target triggers the demand for infill drilling in the upside target of multilayer unconventional reservoirs. However, such an infill scheme in the field practice still heavily relies on empirical knowledge or pressure responses, and the geomechanics consequences have not been fully understood. Backed by the data set from the Permian Basin, in this work we present a novel integrated reservoir-geomechanics-fracture model to simulate the spatiotemporal stress evolution and locate the optimal development strategy in the upside target of the Bone Spring Formation. An embedded discrete fracture model (EDFM) is deployed in our fluid-flow simulation to characterize complex fractures, and the stress-dependent matrix permeability and fracture conductivity are included through the compaction/dilation option. After calibrating reservoir and fracture properties by history matching of an actual well in the development target (i.e., third Bone Spring), we run the finite element method (FEM)-based geomechanics simulation to model the 3D stress state evolution. Then a displacement discontinuity method (DDM) hydraulic fracture model is applied to simulate the multicluster fracture propagation under an updated heterogeneous stress field in the upside target (i.e., second Bone Spring). Numerical results indicate that stress field redistribution associated with parent-well production indeed vertically propagates to the upside target. The extent of stress reorientation at the infill location mainly depends on the parent-child horizontal offset, whereas the stress depletion is under the combined impact of horizontal offset, vertical offset, and infill time. A smaller parent-child horizontal offset aggravates the overlap of the stimulated reservoir volume (SRV), resulting in more substantial interwell interference and less desirable oil and gas production. The same trend is observed by varying the parent-child vertical offset. Moreover, the efficacy of an infill operation at an earlier time is less affected by parent-well depletion because of the less-disturbed stress state. The candidate infill-well locations at various infill timings are suggested based on the parent-well and child-well production cosimulation. Being able to incorporate both pressure and stress responses, the reservoir-geomechanics-fracture model delivers a more comprehensive understanding and a more integral solution of infill-well design in multilayer unconventional reservoirs. The conclusions provide practical guidelines for the subsequent development in the Permian Basin.

Description: Summary The characteristics of hydraulic fractures in the near-wellbore region contain critical information related to the production performance of unconventional wells. We demonstrate a novel application of a fiber-optic-based distributed strain sensing (DSS) technology to measure and characterize near-wellbore fractures and perforation cluster efficiency during production. Distributed fiber-optic-based strain measurements are made based on the frequency shift of the Rayleigh scatter spectrum, which is linearly dependent on strain and temperature changes of the sensing fiber. Strain changes along the wellbore are continuously measured during the shut-in and reopening operations of a well. After removing temperature effects, extensional strain changes can be observed at locations around the perforation cluster during a shut-in period. We interpret that the observed strain changes are caused by near-wellbore fracture aperture changes caused by pressure increases within the near-wellbore fracture network. The depth locations of the measured strain changes correlate well with distributed acoustic sensing (DAS) acoustic intensity measurements that were measured during the stimulation of the well. The shape and magnitude of the strain changes differ significantly between two completion designs in the same well. Different dependencies between strain and borehole pressure can be observed at most of the perforation clusters between the shut-in and reopening periods. We assess that this new type of distributed fiber-optic measurement method can significantly improve understanding of near-wellbore hydraulic fracture characteristics and the relationships between stimulation and production from unconventional oil and gas wells.

Description: Summary Pulse hydraulic fracturing technology can greatly improve the effect of fracture propagation in rock and form complex fracture networks in reservoirs. The interaction mechanism between hydraulic fractures and pre-existing fractures under pulse hydraulic pressure is unclear. The induced laws of pre-existing fractures on the propagation direction of hydraulic fractures under different pulse frequencies and pulse hydraulic pressures are revealed in this work. We have carried out traditional hydraulic fracturing (THF) tests and pulse hydraulic fracturing tests with rock-like specimens. We compared the interaction between hydraulic fractures and pre-existing fractures in the two hydraulic fracturing tests. Acoustic emission (AE) characteristics of the interaction between hydraulic fractures and pre-existing fractures during pulse hydraulic fracturing are analyzed. The results show that pre-existing fractures in the rock-like specimen can induce the direction of propagation of hydraulic fractures. The influence of pre-existing fracture tips on hydraulic fracture propagation is greater with low pulse frequencies than with traditional hydraulic pressures and high pulse frequencies. When the pulse frequency is 1 Hz, hydraulic fractures are easily induced by pre-existing fracture tips. With increasing pulse frequency, the hydraulic fracture propagation direction gradually moves away from the pre-existing fracture tips and extends perpendicularly to the direction of the minimum principal stress. Under pulse hydraulic loading, more hydraulic fractures are generated around the wellbore than under THF and extend to the pre-existing fracture, and more hydraulic fractures around the wellbore are created with low-frequency pulse loading than with high-frequency pulse loading. Compared with traditional hydraulic pressures, hydraulic fracture propagation with low pulse frequencies (1 and 3 Hz) is more complex than hydraulic fracture propagation with traditional hydraulic pressures and high pulse frequencies (5 Hz). Under high pulse hydraulic pressure and pulse frequency, hydraulic fractures easily extend along the direction perpendicular to the direction of the minimum principal stress like propagation under traditional hydraulic pressure. The study of the interaction mechanism between hydraulic fractures and natural fractures under pulsating hydraulic pressure can provide a method for the formation of fracture network systems in large-scale fracturing and may improve the fracturing efficiency.

Description: Summary In this study, we used two-photon polymerization 3D printing technology to successfully print the first true pore-scale rock proxy of Berea sandstone with a submicrometer resolution. Scanning electron microscope (SEM) and computed tomography (CT) images of the 3D-printed sample were compared with the digital file used for printing to verify the rock’s internal structures. Petrophysical properties were estimated with a digital rock physics (DRP) model based on the 3D-printed sample's initial pore network. The results show that our 3D-printing workflow was able to reproduce true-scale 3D porous media such as Berea sandstone with a submicrometer resolution. With a variety of materials and geometric scaling options, 3D printing of nearly identical rock proxies provides a method to conduct repeatable laboratory experiments without destroying natural rock samples. Rock proxy experiments can potentially validate numerical simulations and complement existing laboratory measurements.

Description: Summary Enhanced oil recovery (EOR) in fractured carbonate reservoirs is challenging because of the heterogeneous and oil-wet nature. In this work, a new application of using polymer nanospheres (PNSs) and diluted microemulsion (DME) is presented to plug fractures and enhance water imbibition to recover oil from the tight, naturally fractured carbonate reservoirs. DME with different electric charges is compared through contact-angle and core-imbibition tests to evaluate their performances on EOR. The cationic DME is chosen because it has the fastest wettability-alteration rate and thus the highest oil recovery rate. Migration and plugging efficiency tests are conducted to identify the screened particle sizes of PNSs for the target reservoir cores. PNSs with a particle size of 300 nm are demonstrated to have the best performance of in-depth propagation before swelling and plugging after swelling within the naturally fractured cores are used in this study. Then coreflooding experiments are conducted to evaluate the EOR performance when PNSs and DME are used together, and results indicate that the oil recovery rate is increased by 24.3 and 44.1% compared to using PNSs or DME alone. In the end, a microfluidic experiment is carried out to reveal how DME works with PNSs.

Description: Summary Achieving effective results using conventional thermal recovery technology is challenging in the deep undisturbed reservoir with extra-heavy oil in the LKQ oil field. Therefore, in this study, a novel approach based on in-situ combustion huff-and-puff technology is proposed. Through physical and numerical simulations of the reservoir, the oil recovery mechanism and key injection and production parameters of early-stage ultraheavy oil were investigated, and a series of key engineering supporting technologies were developed that were confirmed to be feasible via a pilot test. The results revealed that the ultraheavy oil in the LKQ oil field could achieve oxidation combustion under a high ignition temperature of greater than 450°C, where in-situ cracking and upgrading could occur, leading to greatly decreased viscosity of ultraheavy oil and significantly improved mobility. Moreover, it could achieve higher extra-heavy-oil production combined with the energy supplement of flue gas injection. The reasonable cycles of in-situ combustion huff and puff were five cycles, with the first cycle of gas injection of 300 000 m3 and the gas injection volume per cycle increasing in turn. It was predicted that the incremental oil production of a single well would be 500 t in one cycle. In addition, the supporting technologies were developed, such as a coiled-tubing electric ignition system, an integrated temperature and pressure monitoring system in coiled tubing, anticorrosion cementing and completion technology with high-temperature and high-pressure thermal recovery, and anticorrosion injection-production integrated lifting technology. The proposed method was applied to a pilot test in the YS3 well in the LKQ oil field. The high-pressure ignition was achieved in the 2200-m-deep well using the coiled-tubing electric igniter. The maximum temperature tolerance of the integrated monitoring system in coiled tubing reached up to 1200°C, which provided the functions of distributed temperature and multipoint pressure measurement in the entire wellbore. The combination of 13Cr-P110 casing and titanium alloy tubing effectively reduced the high-temperature and high-pressure oxygen corrosion of the wellbore. The successful field test of the comprehensive supporting engineering technologies presents a new approach for effective production in deep extra-heavy-oil reservoirs.