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.

Experimental study on secondary air mixing along the bed height in a circulating fluidized bed with a multitracer-gas method

A multitracer-gas method was proposed to study the secondary air(SA)mixing along the bed height in a circulating fluidized bed(CFB)using carbon monoxide(CO),oxygen(O_(2)),and carbon dioxide(CO_(2))as tracer gases.Experiments were carried out on a cold CFB test rig with a cross-section of 0.42 m×0.73 m and a height of 5.50 m.The effects of superficial velocity,SA ratio,bed inventory,and particle diameter on the SA mixing were investigated.The results indicate that there are some differences in the measurement results obtained using different tracer gases,wherein the deviation between CO and CO_(2) ranges from 42%to 66%and that between O_(2) and CO_(2) ranges from 45%to 71%in the lower part of the fluidized bed.However,these differences became less pronounced as the bed height increased.Besides,the high solid concentration and fine particle diameter in the CFB may weaken the difference.The measurement results of different tracer gases show the same trends under the variation of operating parameters.Increasing superficial velocity and SA ratio and decreasing particle diameter result in better mixing of the SA.The effect of bed inventory on SA mixing is not monotonic.

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%.

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..

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.

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.

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.

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.

Dynamic behavior of outburst two-phase flow in a coal mine T-shaped roadway:The formation of impact airflow and its disaster-causing effect

The study of the dynamic disaster mechanism of coal and gas outburst two-phase flow is crucial for improving disaster reduction and rescue ability of coal mine outburst accidents.An outburst test in a T-shaped roadway was conducted using a self-developed large-scale outburst dynamic disaster test system.We investigated the release characteristics of main energy sources in coal seam,and obtained the dynamic characteristics of outburst two-phase flow in a roadway.Additionally,we established a formation model for outburst impact flow and a model for its flow in a bifurcated structure.The results indicate that the outburst process exhibits pulse characteristics,and the rapid destruction process of coal seam and the blocking state of gas flow are the main causes of the pulse phenomenon.The outburst energy is released in stages,and the elastic potential energy is released in the vertical direction before the horizontal direction.In a straight roadway,the impact force oscillates along the roadway.With an increase in the solid–gas ratio,the two-phase flow impact force gradually increases,and the disaster range extends from the middle of the roadway to the coal seam.In the area near the coal seam,the disaster caused by the two-phase flow impact is characterized by intermittent recovery.In a bifurcated roadway,the effect of impact airflow on impact dynamic disaster is much higher than that of two-phase flow,and the impact force tends to weaken with increasing solid-gas ratio.The impact force is asymmetrically distributed;it is higher on the left of the bifurcated roadway.With an increase in the solid-gas ratio,the static pressure rapidly decreases,and the bifurcated structure accelerates the attenuation of static pressure.Moreover,secondary acceleration is observed when the shock wave moves along the T-shaped roadway,indicating that the bifurcated structure increases the shock wave velocity.

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.

Numerical simulation of gas-liquid two-phase jet flow in air-bubble generator

Air-bubble generator is the key part of the self-inspiration type swirl flotation machines,whose flow field structure has a great effect on flotation.The multiphase volume of fluid(VOF),standard k-ε turbulent model and the SIMPLE method were chosen to simulate the present model;the first order upwind difference scheme was utilized to perform a discrete solution for momentum equation.The distributing law of the velocity,pressure,turbulent kinetic energy of every section along the flow direction of air-bubble generator was analyzed.The results indicate that the bubbles are heavily broken up in the middle cross section of throat sect and the entrance of diffuser sect along the flow direction,and the turbulent kinetic energy of diffuser sect is larger than the entrance of throat sect and mixing chamber.

A Numerical Simulation of Air Flow in the Human Respiratory System Based on Lung Model

The lung is an important organ that takes part in the gas exchange process. In the study of gas transport and exchange in the human respiratory system, the complicated process of advection and diffusion (AD) in airways of human lungs is considered. The basis of a lumped parameter model or a transport equation is modeled during the inspiration process, when oxygen enters into the human lung channel. The quantitative measurements of oxygen are detached and the model equation is solved numerically by explicit finite difference schemes. Numerical simulations were made for natural breathing conditions or normal breathing conditions. The respiratory flow results for the resting conditions are found strongly dependent on the AD effect with some contribution of the unsteadiness effect. The contour of the flow rate region is labeled and AD effects are compared with the variation of small intervals of time for a constant velocity when breathing is interrupted for a negligible moment.

Air pressure law of a reservoir constructed in karst sinkholes

Karst sinkholes with natural negative landform provide favorable conditions for the pumped storage reservoir construction for less excavation work.However,the construction of the reservoir would plug the natural karst channels for water and air,which would cause remarkable air pressure in karst channels when the groundwater level fluctuates.A large laboratory simulation test was carried out to study the air pressure variation of a reservoir built on the karst sinkhole.The air pressure in the karst channel and inside the model was monitored during the groundwater rising and falling process.Result showed that the variation of air pressure in the karst channel and the surrounding rock exhibited a high degree of similarity.The air pressure increased rapidly at the initial stage of water level rising,followed by a slight decrease,then the air pressure increased sharply when the water level approached the top of the karst cave.The initial peak of air pressure and the final peak of air pressure were defined,and both air pressure peaks were linearly increasing with the water level rising rate.The negative air pressure was also analyzed during the drainage process,which was linearly correlated with the water level falling rate.The causes of air pressure variation in karst channels of a pumped storage reservoir built on the karst sinkhole were discussed.The initial rapid increase,then slight decrease and final sudden increase of air pressure were controlled by the combined effects of air compression in karst channel and air seepage into the surrounding rock.For the drainage process,the instant negative air pressure and gradual recovering of air pressure were controlled by the combined effects of negative air pressure induced by water level falling and air supply from surrounding rock.This work could provide valuable reference for the reservoir construction in karst area.

Analysis of the Erosion-Corrosion Mechanism of the Air Cooler in a Hydrocracking Unit:A Numerical and Experimental Study

Corrosion leakages often occur in the air cooler of a hydrocracking unit,with the failure sites mainly located in the entrance area of the tubes.An analysis of the macroscopic morphology and corrosion products confirmed that the damage was caused by erosion-corrosion(E-C).Numerical and experimental methods were applied to investigate the E-C mechanism in the air cooler.Computational fluid dynamics(CFD)was used to calculate the hydrodynamic parameters of the air cooler.The results showed that there was a biased flow in the air cooler,which led to a significant increase in velocity,turbulent kinetic energy and wall shear within 0.2 m of the tube entrance.A visualization experiment was then performed to determine the principles of migration and transformation of multiphase flow in the air cooler tubes.Various flow patterns(pure droplet flow,mist flow,and annular flow)and their evolutionary processes were clearly depicted experimentally.The initiation mechanism and processes leading to the development of E-C in the air cooler were also determined.This study provided a comprehensive explanation for the E-C failures that occur in air coolers during operation.

Development and technology status of energy storage in depleted gas reservoirs

Utilizing energy storage in depleted oil and gas reservoirs can improve productivity while reducing power costs and is one of the best ways to achieve synergistic development of"Carbon Peak–Carbon Neutral"and"Underground Resource Utiliza-tion".Starting from the development of Compressed Air Energy Storage(CAES)technology,the site selection of CAES in depleted gas and oil reservoirs,the evolution mechanism of reservoir dynamic sealing,and the high-flow CAES and injection technology are summarized.It focuses on analyzing the characteristics,key equipment,reservoir construction,application scenarios and cost analysis of CAES projects,and sorting out the technical key points and existing difficulties.The devel-opment trend of CAES technology is proposed,and the future development path is scrutinized to provide reference for the research of CAES projects in depleted oil and gas reservoirs.

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.

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.