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.

Simulation of Two-Phase Flowback Phenomena in Shale Gas Wells

The gas-water two-phaseflow occurring as a result of fracturingfluidflowback phenomena is known to impact significantly the productivity of shale gas well.In this work,this two-phaseflow has been simulated in the framework of a hybrid approach partially relying on the embedded discrete fracture model(EDFM).This model assumes the region outside the stimulated reservoir volume(SRV)as a single-medium while the SRV region itself is described using a double-medium strategy which can account for thefluid exchange between the matrix and the micro-fractures.The shale gas adsorption,desorption,diffusion,gas slippage effect,fracture stress sensitivity,and capillary imbibition have been considered.The shale gas production,pore pressure distribution and water saturation distribution in the reservoir have been simulated.The influences of hydraulic fracture geometry and nonorthogonal hydraulic fractures on gas production have been determined and discussed accordingly.The simulation results show that the daily gas production has an upward and downward trend due to the presence of a large amount of fracturingfluid in the reservoir around the hydraulic fracture.The smaller the angle between the hydraulic fracture and the wellbore,the faster the daily production of shale gas wells decreases,and the lower the cumulative production.Nonplanar fractures can increase the control volume of hydraulic fractures and improve the production of shale gas wells.

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.

AQUEOUS TWO-PHASE GAS FLOATATION SPECTROPHOTOMETRIC DETERMINATION OF TRACE TETRACYCLINE IN ENVIRONMENTAL WATER SAMPLE

A green method for separating and enriching trace tetracycline （TC） in environment water by Aqueous Two-phase Gas Floatation Spectrophotometry has been proposed, the principium was discussed. In this paper, the hydrophobic complex composed of Mg（II） and TC was floated into organic phase under the optimal conditions： pH=10, the floatation equipment is home-made, n-propyl alcohol as the organic solvent, sodium chloride as the separating phase reagent. The data were obtained by spectrophotometry after floatation; The linear regression equation is A=2.33×10^5C（mol/L）＋0.2179, linear range is from 3.77×10^-7mol/L to 6.32×10^-5 mol/L, respectively, with thecorrelation coefficient （r） better than 0.9997, relative recoveries is 99.7% to 100.3%, limit of detection was 4.29×10^-5mol/L, The method can be applied to analyse the trace TC in water sample, the result is better.

Ultrasonic method for measuring water holdup of low velocity and high-water-cut oil-water two-phase flow

Oil reservoirs with low permeability and porosity that are in the middle and late exploitation periods in China＇s onshore oil fields are mostly in the high-water-cut production stage.This stage is associated with severely non-uniform local-velocity flow profiles and dispersed-phase concentration（of oil droplets） in oil-water two-phase flow,which makes it difficult to measure water holdup in oil wells.In this study,we use an ultrasonic method based on a transmission-type sensor in oil-water two-phase flow to measure water holdup in lowvelocity and high water-cut conditions.First,we optimize the excitation frequency of the ultrasonic sensor by calculating the sensitivity of the ultrasonic field using the finite element method for multiphysics coupling.Then we calculate the change trend of sound pressure level attenuation ratio with the increase in oil holdup to verify the feasibility of the employed diameter for the ultrasonic sensor.Based on the results,we then investigate the effects of oildroplet diameter and distribution on the ultrasonic field.To further understand the measurement characteristics of the ultrasonic sensor,we perform a flow loop test on vertical upward oilwater two-phase flow and measure the responses of the optimized ultrasonic sensor.The results show that the ultrasonic sensor yields poor resolution for a dispersed oil slug in water flow（D OS/W flow）,but the resolution is favorable for dispersed oil in water flow（D O/W flow） and very fine dispersed oil in water flow（VFD O/W flow）.This research demonstrates the potential application of a pulsed-transmission ultrasonic method for measuring the fraction of individual components in oil-water two-phase flow with a low mixture velocity and high water cut.

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.

Reflection and transmission of bottom simulating reflectors in gas hydrate-bearing sediments: Two-phase media models

The bottom simulating reflector （BSR） in gas hydrate-bearing sediments is a physical interface which is composed of solid, gas, and liquid and is influenced by temperature and pressure. Deep sea floor sediment is a porous, unconsolidated, fluid saturated media. Therefore, the reflection and transmission coefficients computed by the Zoeppritz equation based on elastic media do not match reality. In this paper, a two-phase media model is applied to study the reflection and transmission at the bottom simulating reflector in order to find an accurate wave propagation energy distribution and the relationship between reflection and transmission and fluid saturation on the BSR. The numerical experiments show that the type I compressional （fast） and shear waves are not sensitive to frequency variation and the velocities change slowly over the whole frequency range. However, type II compressional （slow） waves are more sensitive to frequency variation and the velocities change over a large range. We find that reflection and transmission coefficients change with the amount of hydrate and free gas. Frequency, pore fluid saturation, and incident angle have different impacts on the reflection and transmission coefficients. We can use these characteristics to estimate gas hydrate saturation or detect lithological variations in the gas hydrate-bearing sediments.

An analytical model for water-oil two-phase flow in inorganic nanopores in shale oil reservoirs

The existence of water phase occupies oil flow area and impacts the confined oil flow behavior at the solid substrate in inorganic nanopores of shale oil reservoirs,resulting in a completely different flow pattern when compared with the single oil phase flow.This study proposes an analytical model to describe the water-oil two-phase flow.In this model,water slippage at the solid substrate is considered while oil slip is introduced to calculate the oil movement at the solid-oil boundary in dry conditions.It is proven that the oil flow profiles of both the two-phase model and single-phase model show parabolic shapes,but the oil flow capacity drops when water takes up the flow space and the impact of water is more significant when the pore dimension is smaller than 30 nm.Also,the oil flow velocity at a pore center is found to drop linearly given a larger water saturation in wet conditions.The effects of surface wettability and oil properties on water-oil flow are also discussed.Compared with the existing singlephase models,this model describes oil flow pattern in the wet condition with the incorporation of the influence of nanopore properties,which better predicts the oil transport in actual reservoir conditions.Water-oil relative permeability curves are also obtained to improve oil yield.

Aqueous Two-Phase Interfacial Assembly of COF Membranes for Water Desalination

Aqueous two-phase system features with ultralow interfacial tension and thick interfacial region,affording unique confined space for membrane assembly.Here,for the first time,an aqueous two-phase interfacial assembly method is proposed to fabricate covalent organic framework(COF)membranes.The aqueous solution containing polyethylene glycol and dextran undergoes segregated phase separation into two water-rich phases.By respectively distributing aldehyde and amine monomers into two aqueous phases,a series of COF membranes are fabricated at water-water interface.The resultant membranes exhibit high NaCl rejection of 93.0-93.6% and water permeance reaching 1.7-3.7 L m^(−2) h^(−1) bar^(−1),superior to most water desalination membranes.Interestingly,the interfacial tension is found to have pronounced effect on membrane structures.The appropriate interfacial tension range(0.1-1.0 mN m^(−1))leads to the tight and intact COF membranes.Furthermore,the method is extended to the fabrication of other COF and metal-organic polymer membranes.This work is the first exploitation of fabricating membranes in all-aqueous system,confering a green and generic method for advanced membrane manufacturing.

Liquid holdup measurement with double helix capacitance sensor in horizontal oil-water two-phase flow pipes

This paper presents the characteristics of a double helix capacitance sensor for measurement of the liquid holdup in horizontal oil–water two-phase flow. The finite element method is used to calculate the sensitivity field of the sensor in a pipe with 20 mm inner diameter and the effect of sensor geometry on the distribution of sensitivity field is presented. Then, a horizontal oil–water two-phase flow experiment is carried out to measure the response of the double helix capacitance sensor, in which a novel method is proposed to calibrate the liquid holdup based on three pairs of parallel-wire capacitance probes. The performance of the sensor is analyzed in terms of the flow structures detected by mini-conductance array probes.

Parametric study of gas−liquid two-phase flow field in horizontal stirred tank

Based on Fluent software,the gas−liquid two-phase flow in the horizontal stirred tank was simulated with SST k−ωturbulence model,Eulerian−Eulerian two-fluid model,and multi-reference flame method.The mixing process in the tank was calculated by tracer method.The results show that increasing the rotating speed or gas flow is conducive to a more uniform distribution of the gas phase and accelerates the mixing of the liquid phase.When the rotating speed exceeds 93 r/min,the relative power demand remains basically constant.The change in the inclination angle of the upper impeller has minimal effect on the gas phase distribution.When the inclination angle is 50°,the relative power demand reaches the maximum.An appropriate increase in the impeller distance from the bottom improves the gas holdup and gas phase distribution but increases the liquid phase mixing time.

Experimental assessment of two-phase bubble pump for solar water heating

The research goal is to develop a new solar water heater system(SWHS) that uses a solar bubble pump instead of an electric pump.The pump is powered by the steam produced from an evacuated tube collector.Therefore,heat could be transferred downward from the collector to a hot water storage tank.The designed system consists of two sets of heat-pipe evacuated tube collectors,a solar bubble pump installed at an upper level and a water storage tank with a heat exchanger at a lower level.Discharge heads of 1 and 5 m were tested.The bubble pump could operate at the collector temperature of about 90-100 ℃ and vapor gage pressure of 80-90 kPa.It is found that water circulation within the SWHS depends on the incident solar intensity and system discharge head.Experimental investigations are conducted to obtain the system thermal efficiencies from the hourly,daily and long-term performance tests.The thermal performance of the proposed system is compared with conventional solar water heaters.The results show that the proposed system achieves system characteristic efficiency of 10% higher than that of the conventional systems using electric pump if taking the consumption of electric power into account.And the former is a zero carbon system.

Gas-Water Production of a Continental Tight-Sandstone Gas Reservoir under Different Fracturing Conditions

A numerical model of hydraulic fracture propagation is introduced for a representative reservoir(Yuanba continental tight sandstone gas reservoir in Northeast Sichuan).Different parameters are considered,i.e.,the interlayer stress difference,the fracturing discharge rate and the fracturing fluid viscosity.The results show that these factors affect the gas and water production by influencing the fracture size.The interlayer stress difference can effectively control the fracture height.The greater the stress difference,the smaller the dimensionless reconstruction volume of the reservoir,while the flowback rate and gas production are lower.A large displacement fracturing construction increases the fracture-forming efficiency and expands the fracture size.The larger the displacement of fracturing construction,the larger the dimensionless reconstruction volume of the reservoir,and the higher the fracture-forming efficiency of fracturing fluid,the flowback rate,and the gas production.Low viscosity fracturing fluid is suitable for long fractures,while high viscosity fracturing fluid is suitable for wide fractures.With an increase in the fracturing fluid viscosity,the dimensionless reconstruction volume and flowback rate of the reservoir display a non-monotonic behavior,however,their changes are relatively small.

TWO-PHASE FLOW FOR A HORIZONTAL WELL PENETRATING A NATURALLY FRACTURED RESERVOIR WITH EDGE WATER INJECTION

This paper examines the two-phase flow for a horizontal well penetrating a naturally fractured reservoir with edge water injection by means of a fixed streamlinemodel,The mathematical model of the verical two-dimensional flow or oil-water for a horizontal well in a medium with double-porosity is established ,and whose accuratesolutions are obtained by using the characteristic method .The saturation distributionsin the fractured system and the matrix system as well as the formula of the time of water free production are presented .All these results provide a theoretical basis and a computing method for oil displacement by edge water from naturally fracturedreservirs.

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.

Near-source characteristics of two-phase gas-solid outbursts in roadways

Coal and gas outbursts compromise two-phase gas-solid mixtures as they propagate as shock waves and flows from their sources.Propagation is influenced by the form of the outburst,proximity to source,the structure and form of the transmitting roadways and the influence of obstacles.The following characterizes the propagation of coal and gas outbursts as two-phase gas-solid flows proximal to source where the coupled effects of pulverized coal and gas flows dominate behavior.The characteristics of shock wave propagation and attenuation were systematically examined for varied roadway geometries using experiments and numerical models.The results demonstrate that the geometry of roadway obstructions is significant and may result in partial compression and sometimes secondary overpressurization in blocked and small comer roadways leading to significant attenuation of outburst shock waves.The shock waves attenuate slowly in both straight and abruptly expanding roadways and more significantly in T-shaped roadways.The most significant attenuation appears in small angle comers and bifurcations in roadways with the largest attenuation occurring in blocked roadways.These results provide basic parameters for simplifying transport in complex roadway networks in the far-field,and guidance for the design of coal and gas outburst prevention facilities and emergency rescue.

Optimizing Two-Phase Flow Heat Transfer:DCS Hybrid Modeling and Automation in Coal-Fired Power Plant Boilers

In response to escalating challenges in energy conservation and emission reduction,this study delves into the complexities of heat transfer in two-phase flows and adjustments to combustion processes within coal-fired boilers.Utilizing a fusion of hybrid modeling and automation technologies,we develop soft measurement models for key combustion parameters,such as the net calorific value of coal,flue gas oxygen content,and fly ash carbon content,within theDistributedControl System(DCS).Validated with performance test data,thesemodels exhibit controlled root mean square error(RMSE)and maximum absolute error(MAXE)values,both within the range of 0.203.Integrated into their respective automatic control systems,thesemodels optimize two-phase flow heat transfer,finetune combustion conditions,and mitigate incomplete combustion.Furthermore,this paper conducts an in-depth exploration of the generationmechanismof nitrogen oxides(NO_(x))and low oxygen emission reduction technology in coal-fired boilers,demonstrating a substantial reduction in furnace exit NO_(x) generation by 30%to 40%and the power supply coal consumption decreased by 1.62 g/(kW h).The research outcomes highlight the model’s rapid responsiveness,enabling prompt reflection of transient variations in various economic indicator parameters.This provides a more effective means for real-time monitoring of crucial variables in coal-fired boilers and facilitates timely combustion adjustments,underscoring notable achievements in boiler combustion.The research not only provides valuable and practical insights into the intricacies of two-phase flow heat transfer and heat exchange but also establishes a pioneering methodology for tackling industry challenges.

Comparison of Isobaric and Non-isobaric Operations on a Two-Phase Dynamic Model for PEFC Cathode Gas Diffusion Layer

Key words,： Two 1-D dynamical and isothermal models of cathode gas diffusion layer（GDL） with isobaric and non-isobaric operations for polymer electrolyte fuel cells（PEFCs） were developed and implemented in COMSOL Multiphysics v3.5.The artificial diffusion coefficient was introduced as well to make the numerical computation be stable.In the non-isobaric model,the pressure of gas mixture was obtained by summing up the governing equations of gaseous components,instead of Navier-Stoks equation.Comparison of the two models were carried out with the steady-states and dynamical simulations under given conditions.The corresponding analysis based on the simulated results was also given simultaneously.This paper is contributed to finding the differences between the isobaric and non-isobaric operation in the two-phase model of cathode GDL.

NUMERICAL MODELING OF TURBULENTEVAPORATING GAS-DROPLET TWO-PHASE FLOWS IN AN AFTERBURNER DIFFUSOR OF TURBO-FAN JET ENGINES

The two-dimensional turbulent evaporating gas-droplet two-phase flows in an afterburner diffusor of turbo-fan jet engines are simulated by the k-ε turbulence model and the particle trajectory model. Comparison of predicted gas velocity and temperature distributions with experimental results for the cases without liquid spray shows pretty good agreement. Gas-droplet two-phase flow predictions give plausible droplet trajectories, fuel-vapor concentration distribution, gas-phase velocity and temperature field in presence of liquid droplets. One run of computation with this method is made for a particular afterburner. The results indicate that the location of the atomizers is not favorable to flame stabilization and combustion efficiency. The proposed numerical modeling can also be adopted for optimization design and performance evaluation of afterburner combustors of turbo-fan jet engines.

SIMULATION OF GAS AND LIQUID TWO-PHASE FLOW THROUGH THE BLAST FURNACE DROPPING ZONE

A three-dimensional mathematical model,based on differential balances of mass and momentum,hasbeen developed to describe the two-phase flow of gas and liquid through the dropping zone of the blast fur-nace.Agreement between observed and calculated values verifies the validity of this model.On the basis of this model,various parameters for the surrounding of the dry zone of Blast FurnaceNo.I-BF of the Beijing Iron and Steel Company have been computed,from which a diagram for demar-cation of fluidization of coke and flooding of slag has been proposed.