A transient production prediction method for tight condensate gas wells with multiphase flow

Considering the phase behaviors in condensate gas reservoirs and the oil-gas two-phase linear flow and boundary-dominated flow in the reservoir,a method for predicting the relationship between oil saturation and pressure in the full-path of tight condensate gas well is proposed,and a model for predicting the transient production from tight condensate gas wells with multiphase flow is established.The research indicates that the relationship curve between condensate oil saturation and pressure is crucial for calculating the pseudo-pressure.In the early stage of production or in areas far from the wellbore with high reservoir pressure,the condensate oil saturation can be calculated using early-stage production dynamic data through material balance models.In the late stage of production or in areas close to the wellbore with low reservoir pressure,the condensate oil saturation can be calculated using the data of constant composition expansion test.In the middle stages of production or when reservoir pressure is at an intermediate level,the data obtained from the previous two stages can be interpolated to form a complete full-path relationship curve between oil saturation and pressure.Through simulation and field application,the new method is verified to be reliable and practical.It can be applied for prediction of middle-stage and late-stage production of tight condensate gas wells and assessment of single-well recoverable reserves.

CFD modeling of gas−liquid flow phenomenon in lead smelting oxygen-enriched side-blown furnace

A validated numerical model was established to simulate gas−liquid flow behaviors in the oxygen-enriched side-blown bath furnace.This model included the slip velocity between phases and the gas thermal expansion effect.Its modeling results were verified with theoretical correlations and experiments,and the nozzle-eroded states in practice were also involved in the analysis.Through comparison,it is confirmed that the thermal expansion effect influences the flow pattern significantly,which may lead to the backward motion of airflow and create a potential risk to production safety.Consequently,the influences of air injection velocity and furnace width on airflow behavior were investigated to provide operating and design guidance.It is found that the thin layer melt,which avoids high-rate oxygen airflow eroding nozzles,shrinks as the injection velocity increases,but safety can be guaranteed when the velocity ranges from 175 to 275 m/s.Moreover,the isoline patterns and heights of thin layers change slightly when the furnace width increases from 2.2 to 2.8 m,indicating that the furnace width shows a limited influence on production safety.

Integrated numerical simulation of hydraulic fracturing and production in shale gas well considering gas-water two-phase flow

Based on the displacement discontinuity method and the discrete fracture unified pipe network model,a sequential iterative numerical method was used to build a fracturing-production integrated numerical model of shale gas well considering the two-phase flow of gas and water.The model accounts for the influence of natural fractures and matrix properties on the fracturing process and directly applies post-fracturing formation pressure and water saturation distribution to subsequent well shut-in and production simulation,allowing for a more accurate fracturing-production integrated simulation.The results show that the reservoir physical properties have great impacts on fracture propagation,and the reasonable prediction of formation pressure and reservoir fluid distribution after the fracturing is critical to accurately predict the gas and fluid production of the shale gas wells.Compared with the conventional method,the proposed model can more accurately simulate the water and gas production by considering the impact of fracturing on both matrix pressure and water saturation.The established model is applied to the integrated fracturing-production simulation of practical horizontal shale gas wells.The simulation results are in good agreement with the practical production data,thus verifying the accuracy of the model.

Effect of gas blowing nozzle angle on multiphase flow and mass transfer during RH refining process

A three-dimensional mathematical model was developed to investigate the effect of gas blowing nozzle angles on multiphase flow,circulation flow rate,and mixing time during Ruhrstahl-Heraeus(RH) refining process.Also,a water model with a geometric scale of 1:4 from an industrial RH furnace of 260 t was built up,and measurements were carried out to validate the mathematical model.The results show that,with a conventional gas blowing nozzle and the total gas flow rate of 40 L·min^(-1),the mixing time predicted by the mathematical model agrees well with the measured values.The deviations between the model predictions and the measured values are in the range of about 1.3%–7.3% at the selected three monitoring locations,where the mixing time was defined as the required time when the dimensionless concentration is within 3% deviation from the bath averaged value.In addition,the circulation flow rate was 9 kg·s^(-1).When the gas blowing nozzle was horizontally rotated by either 30° or 45°,the circulation flow rate was found to be increased by about 15% compared to a conventional nozzle,due to the rotational flow formed in the up-snorkel.Furthermore,the mixing time at the monitoring point 1,2,and 3 was shortened by around 21.3%,28.2%,and 12.3%,respectively.With the nozzle angle of 30° and 45°,the averaged residence time of 128 bubbles in liquid was increased by around 33.3%.

Permeability evolution and gas flow in wet coal under non-equilibrium state:Considering both water swelling and process-based gas swelling

Accurate knowledge of gas flow within the reservoir and related controlling factors will be important for enhancing the production of coal bed methane.At present,most studies focused on the permeability evolution of dry coal under gas adsorption equilibrium,gas flow and gas diffusion within wet coal under the generally non-equilibrium state are often ignored in the process of gas recovery.In this study,an improved apparent permeability model is proposed which accommodates the water and gas adsorption,stress dependence,water film thickness and gas flow regimes.In the process of modeling,the water adsorption is only affected by water content while the gas adsorption is time and water content dependent;based on poroelastic mechanics,the effective fracture aperture and effective pore radius are derived;and then the variation in water film thickness for different pore types under the effect of water content,stress and adsorption swelling are modeled;the flow regimes are considered based on Beskok’s model.Further,after validation with experimental data,the proposed model was applied to numerical simulations to investigate the evolution of permeability-related factors under the effect of different water contents.The gas flow in wet coal under the non-equilibrium state is explicitly revealed.

Simultaneous measurement of velocity profile and liquid film thickness in horizontal gas–liquid slug flow by using ultrasonic Doppler method

Horizontal gas-liquid two-phase flows widely exist in chemical engineering,oil/gas production and other important industrial processes.Slug flow pattern is the main form of horizontal gas-liquid flows and characterized by intermittent motion of film region and slug region.This work aims to develop the ultrasonic Doppler method to realize the simultaneous measurement of the velocity profile and liquid film thickness of slug flow.A single-frequency single-channel transducer is adopted in the design of the field-programmable gate array based ultrasonic Doppler system.A multiple echo repetition technology is used to improve the temporal-spatial resolution for the velocity profile.An experiment of horizontal gas-liquid two-phase flow is implemented in an acrylic pipe with an inner diameter of 20 mm.Considering the aerated characteristics of the liquid slug,slug flow is divided into low-aerated slug flow,high-aerated slug flow and pseudo slug flow.The temporal-spatial velocity distributions of the three kinds of slug flows are reconstructed by using the ultrasonic velocity profile measurement.The evolution characteristics of the average velocity profile in slug flows are investigated.A novel method is proposed to derive the liquid film thickness based on the instantaneous velocity profile.The liquid film thickness can be effectively measured by detecting the position and the size of the bubbles nearly below the elongated gas bubble.Compared with the time of flight method,the film thickness measured by the Doppler system shows a higher accuracy as a bubble layer occurs in the film region.The effect of the gas distribution on the film thickness is uncovered in three kinds of slug flows.

Numerical simulation of gas–liquid flow in the bubble column using Wray–Agarwal turbulence model coupled with population balance model

In this paper,an improved computational fluid dynamic(CFD)model for gas-liquid flow in bubble column was developed using the one-equation Wary-Agarwal(WA)turbulence model coupled with the population balance model(PBM).Through 18 orthogonal test cases,the optimal combination of interfacial force models,including drag force,lift force,turbulent dispersion force.The modified wall lubrication force model was proposed to improve the predictive ability for hydrodynamic behavior near the wall of the bubble column.The values simulated by optimized CFD model were in agreement with experimental data,and the errors were within±20%.In addition,the axial velocity,turbulent kinetic energy,bubble size distribution,and the dynamic characteristic of bubble plume were analyzed at different superficial gas velocities.This research work could provide a theoretical basis for the extension of the CFD-PBM coupled model to other multiphase reactors..

Inter-layer interference for multi-layered tight gas reservoir in the absence and presence of movable water

Due to the dissimilarity among different producing layers,the influences of inter-layer interference on the production performance of a multi-layer gas reservoir are possible.However,systematic studies of inter-layer interference for tight gas reservoirs are really limited,especially for those reservoirs in the presence of water.In this work,five types of possible inter-layer interferences,including both absence and presence of water,are identified for commingled production of tight gas reservoirs.Subsequently,a series of reservoir-scale and pore-scale numerical simulations are conducted to quantify the degree of influence of each type of interference.Consistent field evidence from the Yan'an tight gas reservoir(Ordos Basin,China)is found to support the simulation results.Additionally,suggestions are proposed to mitigate the potential inter-layer interferences.The results indicate that,in the absence of water,commingled production is favorable in two situations:when there is a difference in physical properties and when there is a difference in the pressure system of each layer.For reservoirs with a multi-pressure system,the backflow phenomenon,which significantly influences the production performance,only occurs under extreme conditions(such as very low production rates or well shut-in periods).When water is introduced into the multi-layer system,inter-layer interference becomes nearly inevitable.Perforating both the gas-rich layer and water-rich layer for commingled production is not desirable,as it can trigger water invasion from the water-rich layer into the gas-rich layer.The gas-rich layer might also be interfered with by water from the neighboring unperforated water-rich layer,where the water might break the barrier(eg weak joint surface,cement in fractures)between the two layers and migrate into the gas-rich layer.Additionally,the gas-rich layer could possibly be interfered with by water that accumulates at the bottom of the wellbore due to gravitational differentiation during shut-in operations.

Numerical Simulation of Oil and Gas Two-Phase Flow in Deep Condensate Gas Reservoirs in Bohai Buried Hills

The BZ19-6 gas field is characterized by high temperature and high pressure (HTHP), high condensate content, little difference between the formation pressure and dew point pressure, and large amount of reverse condensate liquid. During the early stage of depletion development, the production gas-oil ratio (GOR) and production capacity remain relatively stable, which is inconsistent with the conventional reverse condensate seepage law. In view of the static and dynamic conflict in development and production, indoor high-temperature and high-pressure PVT experiment was carried out to reveal the mist-like condensation phenomenon of fluids in the BZ19-6 formation. And the seepage characteristics of condensate gas reservoirs with various degrees of depletion under the condition of HTHP were analyzed based on production performance. The change rule of fluid phase state was analyzed in response to the characterization difficulties of the seepage mechanism. The fluid state was described using the miscible mechanism. And the interphase permeability interpolation coefficient was introduced based on interfacial tension. By doing so, the accurate characterization of the “single-phase flow of condensate gas-near-miscible mist-like quasi single-phase flow-oil-gas two-phase flow” during the development process was achieved. Then the accurate fitting of key indicators for oilfield development was completed, and the distribution law of formation pressure and the law of condensate oil precipitation under different reservoir conditions are obtained. Based on research results, the regulation strategy of variable flow rate production was developed. Currently, the work system has been optimized for 11 wells, achieving a “zero increase” in the GOS of the gas field and an annual oil increase of 22,000 cubic meters.

Simulation of Gas-Water Two-Phase Flow in Tight Gas Reservoirs Considering the Gas Slip Effect

A mathematical model for the gas-water two-phase flow in tight gas reservoirs is elaborated.The model can account for the gas slip effect,stress sensitivity,and high-speed non-Darcy factors.The related equations are solved in the framework of a finite element method.The results are validated against those obtained by using the commercial software CMG(Computer Modeling Group software for advanced recovery process simulation).It is shown that the proposed method is reliable.It can capture the fracture rejection characteristics of tight gas reservoirs better than the CMG.A sensitivity analysis of various control factors(initial water saturation,reservoir parameters,and fracturing parameters)affecting the production in tight gas wells is conducted accordingly.Finally,a series of theoretical arguments are provided for a rational and effective development/exploitation of tight sandstone gas reservoirs.

Effect of gas flow on the nanoparticles transport in dusty acetylene plasmas

This article presents simulation results on the effects of neutral gas flow for nanoparticle transport in atmospheric-pressure,radio-frequency,capacitively-coupled,and acetylene discharge.The acetylene gas is set to flow into the chamber from the upper showerhead electrode.The internal energy of the gas medium therein is transferred into kinetic energy so the gas advection can be triggered.This is represented by the pressure volume work term of the gas energy converse equation.The gas advection leads to the gas temperature sink at the gas inlet,hence a large gas temperature gradient is formed.The thermophoresis relies on the gas temperature gradient,and causes the profile of nanoparticle density to vary from a double-peak structure to a single-peak one.The gas advection influences the properties of electron density and temperature as well and causes the drift-ambipolar mode profile of electron density asymmetric.In the bulk region,i.e.away from the inlet,the gas advection is more like one isovolumetric compression,which slightly increases the temperature of the gas medium at consuming its kinetic energy.

Gas liquid cylindrical cyclone flow regime identification using machine learning combined with experimental mechanism explanation

The flow regimes of GLCC with horizon inlet and a vertical pipe are investigated in experiments,and the velocities and pressure drops data labeled by the corresponding flow regimes are collected.Combined with the flow regimes data of other GLCC positions from other literatures in existence,the gas and liquid superficial velocities and pressure drops are used as the input of the machine learning algorithms respectively which are applied to identify the flow regimes.The choosing of input data types takes the availability of data for practical industry fields into consideration,and the twelve machine learning algorithms are chosen from the classical and popular algorithms in the area of classification,including the typical ensemble models,SVM,KNN,Bayesian Model and MLP.The results of flow regimes identification show that gas and liquid superficial velocities are the ideal type of input data for the flow regimes identification by machine learning.Most of the ensemble models can identify the flow regimes of GLCC by gas and liquid velocities with the accuracy of 0.99 and more.For the pressure drops as the input of each algorithm,it is not the suitable as gas and liquid velocities,and only XGBoost and Bagging Tree can identify the GLCC flow regimes accurately.The success and confusion of each algorithm are analyzed and explained based on the experimental phenomena of flow regimes evolution processes,the flow regimes map,and the principles of algorithms.The applicability and feasibility of each algorithm according to different types of data for GLCC flow regimes identification are proposed.

Simulation of dilatancy-controlled gas migration processes in saturated bentonite using a coupled multiphase flow and elastoplastic H2 M model

Dilatancy-controlled gas flow in preferential pathways plays a key role in the safety analysis of radioactive waste repositories.This is particularly the case for bentonite,an often-preferred barrier material.Gas flow in preferential pathways is characterized by localization and spontaneous behavior,which is challenging to simulate in numerical models due to strong hydro-mechanical coupling.To analyze a laboratory experiment in the framework of the DECOVALEX-2023 project,this study introduced a new approach of combining continuous modelling methods with spatial material properties derived from material heterogeneities and experimental observations.The proposed model utilized hydro-mechanical spatial distributions,namely Young’s modulus and gas entry pressure,and elastoplasticity combined with a linear swelling model.A conceptual strain-dependent permeability approach simulated dilatancycontrolled gas flow based on hydro-mechanical coupling.To test the effectiveness of the presented approach,a gas injection test in a compacted,saturated bentonite sample was simulated using the opensource code OpenGeoSys 5.8 and compared with experimental observations.The presented methodology is capable of simulating localized gas flow in preferential pathways.The spatial distributions of Young’s modulus and gas entry pressure affect the swelling pressure,relative permeability and,in combination with the strain-dependent permeability model,also the intrinsic permeability.

Quantitative research of the liquid film characteristics in upward vertical gas, oil and water flows

The study of liquid film characteristics in multiphase flow is a very important research topic, however,the characteristics of the liquid film around Taylor bubble structure in gas, oil and water three-phase flow are not clear. In the present study, a novel liquid film sensor is applied to measure the distributed signals of the liquid film in three-phase flow. Based on the liquid film signals, the liquid film characteristics including the structural characteristics and the nonlinear dynamics characteristics in three-phase flows are investigated for the first time. The structural characteristics including the proportion, the appearance frequency and the thickness of the liquid film are obtained and the influences of the liquid and gas superficial velocities and the oil content on them are investigated. To investigate the nonlinear dynamics characteristics of the liquid film with the changing flow conditions, the entropy analysis is introduced to successfully uncover and quantify the dynamic complexity of the liquid film behavior.

Gas kinetic flux solver based finite volume weighted essentially non-oscillatory scheme for inviscid compressible flows

A high-order gas kinetic flux solver(GKFS)is presented for simulating inviscid compressible flows.The weighted essentially non-oscillatory(WENO)scheme on a uniform mesh in the finite volume formulation is combined with the circular function-based GKFS(C-GKFS)to capture more details of the flow fields with fewer grids.Different from most of the current GKFSs,which are constructed based on the Maxwellian distribution function or its equivalent form,the C-GKFS simplifies the Maxwellian distribution function into the circular function,which ensures that the Euler or Navier-Stokes equations can be recovered correctly.This improves the efficiency of the GKFS and reduces its complexity to facilitate the practical application of engineering.Several benchmark cases are simulated,and good agreement can be obtained in comparison with the references,which demonstrates that the high-order C-GKFS can achieve the desired accuracy.

Depositional process of hyperpycnal flow deposits:A case study on Lower Cretaceous Sangyuan outcrop in the Luanping Basin,Northeast China

Sedimentary process research is of great significance for understanding the distribution and characteristics of sediments.Through the detailed observation and measurement of the Sangyuan outcrop in Luanping Basin,this paper studies the depositional process of the hyperpycnal flow deposits,and divides their depositional process into three phases,namely,acceleration,erosion and deceleration.In the acceleration phase,hyperpycnal flow begins to enter the basin nearby,and then speeds up gradually.Deposits developed in the acceleration phase are reverse.In addition,the original deposits become unstable and are taken away by hyperpycnal flows under the eroding force.As a result,there are a lot of mixture of red mud pebbles outside the basin and gray mud pebbles within the basin.In the erosion phase,the reverse deposits are eroded and become thinner or even disappear.Therefore,no reverse grading characteristic is found in the proximal major channel that is closer to the source,but it is still preserved in the middle branch channel that is far from the source.After entering the deceleration phase,normally grading deposits appear and cover previous deposits.The final deposits in the basin are special.Some are reverse,and others are normal.They are superimposed with each other under the action of hyperpycnal flow.The analysis of the Sangyuan outcrop demonstrates the sedimentary process and distribution of hyperpycnites,and reasonably explain the sedimentary characteristics of hyperpycnites.It is helpful to the prediction of oil and gas exploration targets in gravity flow deposits.

Insights into the hydrogen evolution reaction in vanadium redox flow batteries:A synchrotron radiation based X-ray imaging study

The parasitic hydrogen evolution reaction(HER)in the negative half-cell of vanadium redox flow batteries(VRFBs)causes severe efficiency losses.Thus,a deeper understanding of this process and the accompanying bubble formation is crucial.This benchmarking study locally analyzes the bubble distribution in thick,porous electrodes for the first time using deep learning-based image segmentation of synchrotron X-ray micro-tomograms.Each large three-dimensional data set was processed precisely in less than one minute while minimizing human errors and pointing out areas of increased HER activity in VRFBs.The study systematically varies the electrode potential and material,concluding that more negative electrode potentials of-200 m V vs.reversible hydrogen electrode(RHE)and lower cause more substantial bubble formation,resulting in bubble fractions of around 15%–20%in carbon felt electrodes.Contrarily,the bubble fractions stay only around 2%in an electrode combining carbon felt and carbon paper.The detected areas with high HER activity,such as the border subregion with more than 30%bubble fraction in carbon felt electrodes,the cutting edges,and preferential spots in the electrode bulk,are potential-independent and suggest that larger electrodes with a higher bulk-to-border ratio might reduce HER-related performance losses.The described combination of electrochemical measurements,local X-ray microtomography,AI-based segmentation,and 3D morphometric analysis is a powerful and novel approach for local bubble analysis in three-dimensional porous electrodes,providing an essential toolkit for a broad community working on bubble-generating electrochemical systems.

Analysis of Maximum Liquid Carrying Capacity Based on Conventional Tubing Plunger Gas Lift

China’s unconventional gas fields have a large number of low-productivity and low-efficiency wells, many of whichare located in remote and environmentally harsh mountainous areas. To address the long-term stable productionof these gas wells, plunger-lift technology plays an important role. In order to fully understand and accurately graspthe drainage and gas production mechanisms of plunger-lift, a mechanical model of plunger-liquid column uplift inthe plunger-lift process was established, focusing on conventional plunger-lift systems and representative wellboreconfigurations in the Linxing region. The operating casing pressure of the plunger-lift process and the calculationmethod for the maximum daily fluid production rate based on the work regime with the highest fluid recovery ratewere determined. For the first time, the critical flow rate method was proposed as a constraint for the maximumliquid-carrying capacity of the plunger-lift, and liquid-carrying capacity charts for conventional plunger-lift withdifferent casing sizes were developed. The results showed that for 23/8 casing plunger-lift, with a well depth ofshallower than 808 m, the maximum drainage rate was 33 m3/d;for 27/8 casing plunger-lift, with a well depth ofshallower than 742 m, the maximum drainage rate was 50.15 m3/d;for 31/2 casing plunger-lift, with a well depthof shallower than 560 m, the maximum drainage rate was 75.14 m3/d. This research provides a foundation for thescientific selection of plunger-lift technology and serves as a decision-making reference for developing reasonableplunger-lift work regimes.

Rate transient analysis methods for water-producing gas wells in tight reservoirs with mobile water

Tight gas reservoirs with mobile water exhibit multi-phase flow and high stress sensitivity.Accurately analyzing the reservoir and well parameters using conventional single-phase rate transient analysis methods proves challenging.This study introduces novel rate transient analysis methods incorporating evaluation processes based on the conventional flowing material balance method and the Blasingame type-curve method to examine fractured gas wells producing water.By positing a gas-water two-phase equivalent homogenous phase that considers characteristics of mobile water,gas,and high stress sensitivity,the conventional single-phase rate transient analysis methods can be applied by integrating the phase's characteristics and defining the phase's normalized parameters and material balance pseudotime.The rate transient analysis methods based on the equivalent homogenous phase can be used to quantitatively assess the parameters of wells and gas reservoirs,such as original gas-in-place,fracture half-length,reservoir permeability,and well drainage radius.This facilitates the analysis of production dynamics of fractured wells and well-controlled areas,subsequently aiding in locating residual gas and guiding the configuration of well patterns.The specific evaluation processes are detailed.Additionally,a numerical simulation mechanism model was constructed to verify the reliability of the developed methods.The methods introduced have been successfully implemented in field water-producing gas wells within tight gas reservoirs containing mobile water.

PIV MEASUREMENT FOR SWIRLER FLOW FIELD IN GAS TURBINE COMBUSTOR

The characteristics of swirler flow field, including cold flow field and combustion flow field, in gas tur- bine combustor with two-stage swirler are studied by using particle image velocimetry （PIV）. Velocity compo- nents, fluctuation velocity, Reynolds stress and recirculation zone length are obtained, respectively. Influences of geometric parameter of primary hole, arrangement of primary hole, inlet air temperature, first-stage swirler an- gle and fuel/air ratio on flow field are investigated, respectively. The experimental results reveal that the primary recirculation zone lengths of combustion flow field are shorter than those of cold flow field, and the primary reeir- culation zone lengths decrease with the increase of inlet air temperature and fuel/air ratio. The change of the geo- metric parameter of primary hole casts an important influence on the swirler flow field in two-stage swirler com- bustor.