A New Model for Predicting Liquid Loading in Shale Gas Horizontal Wells

Current critical flow rate models fail to accurately predict the liquid loading statuses of shale gas horizontal wells. Therefore, a new critical flow rate model for the whole wellbore of shale gas horizontal wells is established. The results of the new model are compared to those of current models through the field case analysis. The new model is based on the dynamic analysis and energy analysis of the deformed liquid-droplet, which takes into account the liquid flow rate, the liquid-droplet deformation and the energy loss caused by the change of buildup rate. The major axis of the maximum stable deformed liquid-droplet is determined based on the energy balance relation. Meanwhile, the suitable drag coefficient equation and surface tension equation applied to shale gas horizontal wells are chosen. Finally, the critical flow rate equation is established and the maximum critical flow rate of the whole wellbore is chosen as the criterion for liquid loading prediction. The precision of liquid loading prediction of the new model is compared to those of the four current models, including Belfroid's model, modified Li's model, liquid film model and modified Wang's model. Field parameters of 29 shale gas horizontal wells are used for the comparison, including parameters of 18 unloaded wells, 2 near loaded-up wells and 9 loaded-up wells. Field case analysis shows that the total precision of liquid loading prediction of the new model is 93.1%, which is higher compared to those of the current four models. The new model can accurately predict the liquid loading statuses of loaded-up wells and near loaded-up wells, while the prediction precision for unloaded wells is high enough for the field application, which is 88.9%. The new model can be used to effectively estimate the field liquid loading statuses of shale gas horizontal wells and choose drainage gas recovery technologies, which considers both the complex wellbore structure and the variation of flowback liquid flow rate in shale gas horizontal wells. The results of the new model fill the gap in existing studies and have a guiding significance for liquid loading prediction in shale gas horizontal wells.

Part Load Instability and Rotating Stall in a Multistage Low Specific Speed Pump

A new diffuser design is developed for a low specific speed, multistage pump. In this design the diffuser and the de-swirl vanes are integrated into single vanes. This creates diffuser channels that extend from behind the impeller exit through the cross-over, up to the eye of the next stage impeller. Experiments show the occurrence of a saddle type instability in the head curve. At a critical flow rate of close to 50% of the flow rate at Best Efficiency Point (BEP), the head drops by 7% of the head at BEP. In this study Computational Fluid Dynamics (CFD) are used in an effort to understand the underlying flow phenomena. The head curve that is obtained with the transient CFD simulations contains a saddle type instability at a flow rate that is approximately the same as in the experiments, but with a lower magnitude. At flow rates higher than the critical flow rate, the predicted head and power are in very good agreement with the experimental data. At flow rates lower than the critical flow rate, the head and power are slightly over-predicted. An analysis of the pressure distribution in the pump reveals that the head loss at different flow rates in the diffuser shows a discontinuity at the critical flow rate. Since both the impeller head and the head loss in the vaneless gap increase continuously for decreasing flow rate, this is an indication that the cause of the head instability lies in the diffuser. Moreover, a strong increase in the variability of head and power at flow rates below the critical flow suggests that the phenomenon is unsteady. Flow patterns in the impeller and in the diffuser, as calculated by CFD, show a high degree of periodicity and are very similar for flow rates down to the critical flow rate. However, for lower flow rates the flow pattern changes completely. A single rotating stall cell is observed that causes two or three neighboring diffuser channels to stall, leading to a significantly lower flow rate or even a reversed flow. This stall pattern rotates in the direction of impeller rotation at a very low frequency of approximately 3.3% of the impeller rotation frequency.

Role of Shearing Dispersion and Stripping in Wax Deposition in Crude Oil Pipelines

Wax deposition during crude oil transmission can cause a series of negative effects and lead to problems associated with pipeline safety. A considerable number of previous works have investigated the wax deposition mechanism, inhibition technology, and remediation methods. However, studies on the shearing mechanism of wax deposition have focused largely on the characterization of this phenomena. The role of the shearing mechanism on wax deposition has not been completely clarified. This mechanism can be divided into the shearing dispersion effect caused by radial migration of wax particles and the shearing stripping effect caused by hydrodynamic scouring. From the perspective of energy analysis, a novel wax deposition model was proposed that considered the flow parameters of waxy crude oil in pipelines instead of its rheological parameters. Considering the two effects of shearing dispersion and shearing stripping coexist, with either one of them being the dominant mechanism, a shearing dispersion flux model and a shearing stripping model were established. Furthermore, a quantitative method to distinguish between the roles of shearing dispersion and shearing stripping in wax deposition was developed. The results indicated that the shearing mechanism can contribute an average of approximately 10% and a maximum of nearly 30% to the wax deposition process. With an increase in the oil flow rate, the effect of the shearing mechanism on wax deposition is enhanced, and its contribution was demonstrated to be negative; shear stripping was observed to be the dominant mechanism. A critical flow rate was observed when the dominant effect changes. When the oil flow rate is lower than the critical flow rate, the shearing dispersion effect is the dominant effect; its contribution rate increases with an increase in the oil flow temperature. When the oil flow rate is higher than the critical flow rate, the shearing stripping effect is the dominant effect; its contribution rate increases with an increase in the oil flow temperature. This understanding can be used to design operational parameters of the actual crude oil pipelines and address the potential flow assurance problems. The results of this study are of great significance for understanding the wax deposition theory of crude oil and accelerating the development of petroleum industry pipelines.

Numerical Simulation of Vibration of Composite Pipelines Conveying Pulsating Fluid

Vibration problems of pipelines made of composite materials conveying pulsating flow of gas and fluid are investigated in the paper. A dynamic model of motion of pipelines conveying pulsating fluid flow supported by a Hetenyi’s base is developed taking into account the viscosity properties of the structure material, axial forces, internal pressure and Winkler’s viscoelastic base. To describe the processes of viscoelastic material strain, the Boltzmann–Volterra integral model with weakly singular hereditary kernels is used. Using the Bubnov–Galerkin method, the problem is reduced to the study of a system of ordinary integro-differential equations (IDE). A computational algorithm is developed based on the elimination of the features of IDE with weakly singular kernels, followed by the use of quadrature formulas. The effect of rheological parameters of the pipeline material, flow rate and base parameters on the vibration of a viscoelastic pipeline conveying pulsating fluid is analyzed. The convergence analysis of the approximate solution of the Bubnov–Galerkin method is carried out. It was revealed that the viscosity parameters of the material and the pipeline base lead to a significant change in the critical flow rate. It was stated that an increase in excitation coefficient of pulsating flow and the parameter of internal pressure leads to a decrease in the critical flow rate. It is shown that an increase in the singularity parameter, the Winkler base parameter, the rigidity parameter of the continuous base layer and the Reynolds number increases the critical flow rate.

Study of Stable Combustion Ranges of Diffusion Flame Microjet of Hydrogen Effluent from around Micronozzle at Addition of Inert and Reactive Gases into Hydrogen or Air

The effect addition of inert (He, N2, Ar, CO2) and reacting (СН4, O2, CF3Br, (CH3O)PO) gases in a hydrogen or in coflowing air stream on lift-off diffusion flame conditions of hydrogen micro jet effluxed from the round micro nozzle was experimentally studied. Using the schlieren technique a critical hydrogen flow rate was established at which the flame of the hydrogen microjet detaches from the nozzle when introducing additives of the studied gases into both air and hydrogen. It has been established that the addition of all the studied gases to the hydrogen leads to a decrease in the velocity range of its microjet at which flame stabilization is possible, regardless of whether the gases introduced into hydrogen are inert or reactive. It is shown that in the case of the addition of various gases to the hydrogen, the main factor determining the critical flow rate at flame lift-off from the micro nozzle is the average molecular weight of the H2 gas mixture with additives. At addition of the studied gases into the coflowing air, the critical flow rate of H2 is determined by their affect on the chemical reactions of hydrogen oxidation (inhibition effectiveness), as well as by a decrease in oxygen concentration due to dilution of air by additives. The data obtained are of interest to hydrogen energy in terms of determining the limits of sustainable combustion of the hydrogen microjet, as well as determining the minimum phlegmatizing concentrations of additives of inhibitors and fire suppressant in the air, preventing the ignition and explosion of hydrogen in case of emergency at its leak.

Effect of Sand Bed Deposits on the Characteristics of Turbulent Flow of Water in Horizontal Annuli

An experimental program was conducted to investigate turbulent flow of water over the stationary sand bed deposited in horizontal annuli. A large-scale horizontal flow loop equipped with the state-of-the-art particle image velocimetry (PIV) system has been used for the experiments. Experiments were conducted to measure the instantaneous local velocity profiles during turbulent flow and examine the impact of the presence of a stationary sand bed deposits on the local velocity profiles, Reynolds shear stresses and turbulence intensities. Results have shown that the existence of a stationary sand bed causes the volumetric flow to be diverted away from the lower annular gap. Increasing the sand bed height causes further reduction of the volumetric flow rate in the lower annulus. Velocity profiles near the surface of the bed deposits showed a downward shift from the universal law in wall units indicating that the flow is hydraulically rough near the sand bed. The equivalent roughness height varied with flow rates. At flow rates less than the critical flow rate, the Reynolds stress profile near the bed interface had slightly higher peak values than that of the case with no sand bed. At the critical flow rate, however, the peak Reynolds stress values for the flow over the sand bed was lower than that of the case with no bed. This behavior is attributed to the bed load transport of sand particles at the critical flow rate.

Determination of Optimum Concentration of Nanofluid for Process Intensification of Heat Transfer Using Corrugated Plate Type Heat Exchanger

In this study cooling of hot water is taken up using a compact heat exchanger such as corrugated plate type heat exchanger and with utility fluid as a nanofluid prepared from mixing Al2O3 in water. In general a monotonic increase of 30 % to 70 % in overall heat transfer coefficient is observed for increase in nanofluid concentration as well as its flow rate. An optimum concentration of nanofluid is hence not possible to be found as heat transfer coefficient exhibited a monotonic trend. But there is a penalty for using nanofluid of higher concentrations in heat exchangers in the form of additional hydraulic power required to supply the nanofluid due its higher viscosity. Hence, as a novel approach, a target temperature drop of 15 °C for hot fluid (with constant flow rate) is assumed and the minimum critical flow rate of cold nanofluid of various concentrations required for achieving this is determined using simulation. For such a critical flow rate at various nanofluid concentrations, the determined hydraulic power (product of pressure drop and flow rate) exhibited a global minimum around 0.75 % volume concentration of Al2O3 in water. Thus this article presents the process intensification procedure for the heat exchangers using nanofluids as heat transfer enhancement option.

Effect of condensate flow rate on retention angle on horizontal low-finned tubes

The paper reports experimental results using simulated condensation on eight hor?izontal integral finned tubes with different fin spacing but same root diameter. Condensation was simulated with low approaching zero vapor velocity of condensate using three liquids (water, ethylene glycol, and R141b) supplied to the tube via small holes between the fins along the top of the tubes. Controlling parameters of the investigation were fin spacing of condensation tubes, flow rate of condensate and surface tension to density ratio of the condensate. The results indicate that the retention angle (measured from the top of the tube to the position where the inter-fin space is completely filled with liquid) increases with the increase in fin spacing. Also, retention angle increases as the density of the condensate increases but retention angle decreases with increase in surface tension. Interesting finding is seen as retention angle remains constant with increase in condensate flow rate, starting from very low (nearly zero) flow rate to the flow rate at which the tube gets fully flooded. The critical flow rate for eight tubes of defined fin density against three working fluids is measured. Results obtained from simulated condensation for almost zero condensate velocity are in good agreement with earlier data and theoretical model for retention angle on such tubes

Investigation of Diameter Effect on Critical Flow

Discrepancy has long existed about nozzle diameter effect on critical flow rate in two-phase critical flow analysis. The issue is of great importance because it involves whether the scaling test can accurately simulate the actual process during loss of coolant accident. A series of transient critical flow experiments has been performed in China Institute of Atomic Energy (CIAE), which aims to study the nozzle diameter effect. The diameter of adopted test nozzle is 5mm, 10mm and 15mm, respectively. The experiment result shows that the discharge mass flow rate decreases as nozzle diameter increases. This is contradictory with the results obtained from steady state critical flow experiment. Comprehensive analysis shows that the diameter effect observed in transient critical flow experiment is transient effect caused by system pressure change, and the mechanism is the dynamic unbalance in the transient period. It is presumed that for the same test section with a certain diameter, if the pressure change velocity is different, the measured critical flow rate will be different. The conclusion is validated in the pressurization and depressurization blowdown test. It| is arrived that nozzle diameter has no direct effect on critical flow rate, and thus the discrepancy on diameter effect is clarified.