Pressure transient characteristics of non-uniform conductivity fractured wells in viscoelasticity polymer flooding based on oil-water two-phase flow

Polymer flooding in fractured wells has been extensively applied in oilfields to enhance oil recovery.In contrast to water,polymer solution exhibits non-Newtonian and nonlinear behavior such as effects of shear thinning and shear thickening,polymer convection,diffusion,adsorption retention,inaccessible pore volume and reduced effective permeability.Meanwhile,the flux density and fracture conductivity along the hydraulic fracture are generally non-uniform due to the effects of pressure distribution,formation damage,and proppant breakage.In this paper,we present an oil-water two-phase flow model that captures these complex non-Newtonian and nonlinear behavior,and non-uniform fracture characteristics in fractured polymer flooding.The hydraulic fracture is firstly divided into two parts:high-conductivity fracture near the wellbore and low-conductivity fracture in the far-wellbore section.A hybrid grid system,including perpendicular bisection(PEBI)and Cartesian grid,is applied to discrete the partial differential flow equations,and the local grid refinement method is applied in the near-wellbore region to accurately calculate the pressure distribution and shear rate of polymer solution.The combination of polymer behavior characterizations and numerical flow simulations are applied,resulting in the calculation for the distribution of water saturation,polymer concentration and reservoir pressure.Compared with the polymer flooding well with uniform fracture conductivity,this non-uniform fracture conductivity model exhibits the larger pressure difference,and the shorter bilinear flow period due to the decrease of fracture flow ability in the far-wellbore section.The field case of the fall-off test demonstrates that the proposed method characterizes fracture characteristics more accurately,and yields fracture half-lengths that better match engineering reality,enabling a quantitative segmented characterization of the near-wellbore section with high fracture conductivity and the far-wellbore section with low fracture conductivity.The novelty of this paper is the analysis of pressure performances caused by the fracture dynamics and polymer rheology,as well as an analysis method that derives formation and fracture parameters based on the pressure and its derivative curves.

Transient pressure analysis of polymer flooding fractured wells with oil-water two-phase flow

The oil-water two-phase flow pressure-transient analysis model for polymer flooding fractured well is established by considering the comprehensive effects of polymer shear thinning,shear thickening,convection,diffusion,adsorption retention,inaccessible pore volume and effective permeability reduction.The finite volume difference and Newton iteration methods are applied to solve the model,and the effects of fracture conductivity coefficient,injected polymer mass concentration,initial polymer mass concentration and water saturation on the well-test type curves of polymer flooding fractured wells are discussed.The results show that with the increase of fracture conductivity coefficient,the pressure conduction becomes faster and the pressure drop becomes smaller,so the pressure curve of transitional flow goes downward,the duration of bilinear flow becomes shorter,and the linear flow appears earlier and lasts longer.As the injected polymer mass concentration increases,the effective water phase viscosity increases,and the pressure loss increases,so the pressure and pressure derivative curves go upward,and the bilinear flow segment becomes shorter.As the initial polymer mass concentration increases,the effective water phase viscosity increases,so the pressure curve after the wellbore storage segment moves upward as a whole.As the water saturation increases,the relative permeability of water increases,the relative permeability of oil decreases,the total oil-water two-phase mobility becomes larger,and the pressure loss is reduced,so the pressure curve after the wellbore storage segment moves downward as a whole.The reliability and practicability of this new model are verified by the comparison of the results from simplified model and commercial well test software,and the actual well test data.

Characterization of Flame Retardancy and Oil-Water Separation Capacity of Superhydrophobic Silylated Melamine Sponges

A silylated melamine sponge(SMS)was prepared by two simple steps,namely,immersion and dehydration of a melamine sponge coated with methyltrichlorosilane.The silylated structure of SMS was characterized by FT-IR(Fourier-transform infrared)spectroscopy,SEM(Scanning electron microscopy)and in terms of water contact angles.Its oil-water absorption and separation capacities were measured by FT-IR and UV-visible spectrophoto-metry.The experimental results have shown that oligomeric silanol covalently bonds by Si-N onto the surface of melamine sponge skeletons.SMS has shown superhydrophobicity with a water contact angle exceeding 150°±1°,a better separation efficiency with regard to diesel oil(by 99.31%(wt/wt%)in oil-water mixture and even up to 99.99%(wt/wt%)for diesel oil in its saturated aqueous solution.Moreover,SMS inherited the intrinsicflame retardancy of the melamine sponge.In general,SMS has shown superhydrophobicity,high porosity,excellent selectivity,remarkable recyclability,and better absorption capacity for various oils and organic solvents,and a high separation efficiency for oil in saturated aqueous solutions.

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.

Arbitrary High-Order Fully-Decoupled Numerical Schemes for Phase-Field Models of Two-Phase Incompressible Flows

Due to the coupling between the hydrodynamic equation and the phase-field equation in two-phase incompressible flows,it is desirable to develop efficient and high-order accurate numerical schemes that can decouple these two equations.One popular and efficient strategy is to add an explicit stabilizing term to the convective velocity in the phase-field equation to decouple them.The resulting schemes are only first-order accurate in time,and it seems extremely difficult to generalize the idea of stabilization to the second-order or higher version.In this paper,we employ the spectral deferred correction method to improve the temporal accuracy,based on the first-order decoupled and energy-stable scheme constructed by the stabilization idea.The novelty lies in how the decoupling and linear implicit properties are maintained to improve the efficiency.Within the framework of the spatially discretized local discontinuous Galerkin method,the resulting numerical schemes are fully decoupled,efficient,and high-order accurate in both time and space.Numerical experiments are performed to validate the high-order accuracy and efficiency of the methods for solving phase-field models of two-phase incompressible flows.

Identifying the enhancement mechanism of Al/MoO_(3) reactive multilayered films on the ignition ability of semiconductor bridge using a one-dimensional gas-solid two-phase flow model

Energetic Semiconductor bridge(ESCB)based on reactive multilayered films(RMFs)has a promising application in the miniature and intelligence of initiator and pyrotechnics device.Understanding the ignition enhancement mechanism of RMFs on semiconductor bridge(SCB)during the ignition process is crucial for the engineering and practical application of advanced initiator and pyrotechnics devices.In this study,a one-dimensional(1D)gas-solid two-phase flow ignition model was established to study the ignition process of ESCB to charge particles based on the reactivity of Al/MoO_(3) RMFs.In order to fully consider the coupled exothermic between the RMFs and the SCB plasma during the ignition process,the heat release of chemical reaction in RMFs was used as an internal heat source in this model.It is found that the exothermal reaction in RMFs improved the ignition performance of SCB.In the process of plasma rapid condensation with heat release,the product of RMFs enhanced the heat transfer process between the gas phase and the solid charge particle,which accelerated the expansion of hot plasma,and heated the solid charge particle as well as gas phase region with low temperature.In addition,it made up for pressure loss in the gas phase.During the plasma dissipation process,the exothermal chemical reaction in RMFs acted as the main heating source to heat the charge particle,making the surface temperature of the charge particle,gas pressure,and gas temperature rise continuously.This result may yield significant advantages in providing a universal ignition model for miniaturized ignition devices.

A Note on an Order Level Inventory Model with Varying Two-Phased Demand and Time-Proportional Deterioration

The main purpose of this paper is to generalize the effect of two-phased demand and variable deterioration within the EOQ (Economic Order Quantity) framework. The rate of deterioration is a linear function of time. The two-phased demand function states the constant function for a certain period and the quadratic function of time for the rest part of the cycle time. No shortages as well as partial backlogging are allowed to occur. The mathematical expressions are derived for determining the optimal cycle time, order quantity and total cost function. An easy-to-use working procedure is provided to calculate the above quantities. A couple of numerical examples are cited to explain the theoretical results and sensitivity analysis of some selected examples is carried out.

Complex network analysis in inclined oil-water two-phase flow

Complex networks have established themselves in recent years as being particularly suitable and flexible for representing and modelling many complex natural and artificial systems. Oil-water two-phase flow is one of the most complex systems. In this paper, we use complex networks to study the inclined oil water two-phase flow. Two different complex network construction methods are proposed to build two types of networks, i.e. the flow pattern complex network （FPCN） and fluid dynamic complex network （FDCN）. Through detecting the community structure of FPCN by the community-detection algorithm based on K-means clustering, useful and interesting results are found which can be used for identifying three inclined oil-water flow patterns. To investigate the dynamic characteristics of the inclined oil-water two-phase flow, we construct 48 FDCNs under different flow conditions, and find that the power-law exponent and the network information entropy, which are sensitive to the flow pattern transition, can both characterize the nonlinear dynamics of the inclined oil-water two-phase flow. In this paper, from a new perspective, we not only introduce a complex network theory into the study of the oil-water two-phase flow but also indicate that the complex network may be a powerful tool for exploring nonlinear time series in practice.

Relative permeability estimation of oil-water two-phase flow in shale reservoir

Oil-water two-phase flow is ubiquitous in shale strata due to the existence of connate water and the injection of fracturing fluid.In this work,we propose a relative permeability model based on a modified Hagen-Poiseuille(HP)equation and shale reconstruction algorithm.The proposed model can consider the nanoconfined effects(slip length and spatially varying viscosity),oil-water distribution,pore size distribution(PSD),total organic matter content(TOC),and micro-fracture.The results show that the increasing contact angles of organic matters(OM)and inorganic minerals(iOM)increase the relative permeability of both oil and water.As the viscosity ratio increases,the relative permeability of oil phase increases while that of water phase decreases,due to the different water-oil distribution.The effective permeability of both oil and water decreases with the increasing TOC.However,the relative permeability of water phase increases while that of oil phase decreases.The increasing number and decreasing deviation angle of micro-fracture increase the effective permeability of oil and water.However,microfracture has a minor effect on relative permeability.Our model can help understand oil-water twophase flow in shale reservoirs and provide parameter characterization for reservoir numerical simulation.

Volume fraction measurement of oil-water two-phase flow using a coaxial conductivity sensor

The accurate measurement of volume fraction of oil-water two-phase flow on line is important in the oil field.This paper presents a new coaxial conductivity sensor for measuring the volume fraction of oil-water two-phase flow.This structure may get the more uniform sensitivity field and the vertical installation may get the more axial symmetry of the flow field,which improve the measurement accuracy.In order to minimize the influence of the edge effect,guard electrodes were designed.An anti-edge effect degree Ae was defined to optimize the length of the guard electrode.Different models of effective conductivity of two materials were used in calculating the oil volume fraction of oil-water two-phase flow.The experimental results indicate that Maxwell model is the best model under the condition of oil volume fraction less than 50%and the mean value of the calculation results using Maxwell model and Bruggeman model possesses higher accuracy in the range of oil volume fraction(50%-70%).The experimental results show that the sensor obtains similar measurement performance in both vertical upward and downward flow conditions.The accuracy of the sensor system is 2%when the oil volume fraction less than 50%,and the accuracy is about 5%when the oil volume fraction between 50%and 70%.

Numerical modeling and parametric sensitivity analysis of heat transfer and two-phase oil and water flow characteristics in horizontal and inclined flowlines using OpenFOAM

Estimating the oil-water temperatures in flowlines is challenging especially in deepwater and ultra-deepwater offshore applications where issues of flow assurance and dramatic heat transfer are likely to occur due to the temperature difference between the fluids and the surroundings. Heat transfer analysis is very important for the prediction and prevention of deposits in oil and water flowlines, which could impede the flow and give rise to huge financial losses. Therefore, a 3D mathematical model of oil-water Newtonian flow under non-isothermal conditions is established to explore the complex mechanisms of the two-phase oil-water transportation and heat transfer in different flowline inclinations. In this work, a non-isothermal two-phase flow model is first modified and then implemented in the InterFoam solver by introducing the energy equation using OpenFOAM® code. The Low Reynolds Number (LRN) k-ε turbulence model is utilized to resolve the turbulence phenomena within the oil and water mixtures. The flow patterns and the local heat transfer coefficients (HTC) for two-phase oil-water flow at different flowlines inclinations (0°, +4°, +7°) are validated by the experimental literature results and the relative errors are also compared. Global sensitivity analysis is then conducted to determine the effect of the different parameters on the performance of the produced two-phase hydrocarbon systems for effective subsea fluid transportation. Thereafter, HTC and flow patterns for oil-water flows at downward inclinations of 4°, and 7° can be predicted by the models. The velocity distribution, pressure gradient, liquid holdup, and temperature variation at the flowline cross-sections are simulated and analyzed in detail. Consequently, the numerical model can be generally applied to compute the global properties of the fluid and other operating parameters that are beneficial in the management of two-phase oil-water transportation.

A two-phase type-curve method with multiscale fluid transport mechanisms in hydraulically fractured shale reservoirs

The quantitative understanding of hydraulic fracture(HF)properties guides accurate production forecasts and reserve estimation.Type curve is a powerful technique to characterize HF and reservoir properties from flowback and long-term production data.However,two-phase flow of water and hydrocarbon after an HF stimulation together with the complex transport mechanisms in shale nanopores exacerbate the nonlinearity of the transport equation,causing errors in type-curve analysis.Accordingly,we propose a new two-phase type-curve method to estimate HF properties,such as HF volume and permeability of fracture,through the analysis of flowback data of multi-fractured shale wells.The proposed type curve is based on a semianalytical solution that couples the two-phase flow from the matrix with the flow in HF by incorporating matrix influx,slippage effect,stress dependence,and the spatial variation of fluid properties in inorganic and organic pores.For the first time,multiple fluid transport mechanisms are considered into two-phase type-curve analysis for shale reservoirs.We analyze the flowback data from a multi-fractured horizontal well in a shale gas reservoir to verify the field application of the proposed method.The results show that the fracture properties calculated by the type-curve method are in good agreement with the long-time production data.

Pore-Scale Investigation of Coupled Two-Phase and Reactive Transport in the Cathode Electrode of Proton Exchange Membrane Fuel Cells

A three-dimensional multicomponent multiphase lattice Boltzmann model(LBM)is established to model the coupled two-phase and reactive transport phenomena in the cathode electrode of proton exchange membrane fuel cells.The gas diff usion layer(GDL)and microporous layer(MPL)are stochastically reconstructed with the inside dynamic distribution of oxygen and liquid water resolved,and the catalyst layer is simplifi ed as a superthin layer to address the electrochemical reaction,which provides a clear description of the fl ooding eff ect on mass transport and performance.Diff erent kinds of electrodes are reconstructed to determine the optimum porosity and structure design of the GDL and MPL by comparing the transport resistance and per-formance under the fl ooding condition.The simulation results show that gradient porosity GDL helps to increase the reactive area and average concentration under fl ooding.The presence of the MPL ensures the oxygen transport space and reaction area because liquid water cannot transport through micropores.Moreover,the MPL helps in the uniform distribution of oxygen for an effi cient in-plane transport capacity.Crack and perforation structures can accelerate the water transport in the assembly.The systematic perforation design yields the best performance under fl ooding by separating the transport of liquid water and oxygen.

Preparation of Silica Nanosphere with Vertical Pore and Its Application in Oil-water Separation

Uniform monodispersed mesoporous silica nanospheres with vertical pores were successfully synthesized using chiral amphiphilic small molecule L-16Ala5PyClO4and solvents as dual templates via solgel transcription.The morphologies and pore sizes of silicas are adjustable by changing the type and amount of solvents in the reaction systems.With the increase of the organic solvent content,the morphologies of the obtained silica changed from nanospheres with vertical pore structures to nanosheets structures.When 1 mL of benzene,cyclohexane or toluene were used as solvents,only silica nanospheres were obtained,the BET surface areas of silica nanospheres reached 600.7,669.5,and 560.8 m^(2)/g,respectively.The pore sizes were 3.51,3.54,and 3.46 nm,respectively.Significantly,these ordered silica nanospheres/poly(vinyl alcohol-co-ethylene)(PVAco-PE)nanofiber membranes have high separation efficiencies(>99%)for n-hexane/water mixtures.

Wax deposition modeling in oil-water stratified pipe flow

Wax deposition in oil-water stratified flow is commonly encountered onshore and offshore oil production pipe systems,and typically reduces transportation capacity of oil.The accurate predicted model of wax deposition has becomes an indispensable approach to design effective remediation strategies.However,a reliable mechanistic model for wax deposition prediction in oil-water two-phase stratified pipe flow is lacking to validate the deposition process.In this work,a three-dimensional(axial,radial,and angular)robust wax deposit model for oil-water stratified circular pipe flow was developed.The model of formation of a gel deposit based on the first principles of rheology was developed,associated with the results obtained from hydrodynamics and heat/mass transfer simulations.The predictions for wax deposition are found to compare satisfactorily with experimental data with two different oils for single phase and four different water cuts for oil-water stratified pipe flow.It can be seen from the wax gelation mechanism that an increase in water cut can help to reduce the wall/oil-deposit interface shear stress,thereby leading to an increase in the degree of gelation as well as the deposit rate.Furthermore,a local deposit analysis in the circumferential direction was conducted,for water cut 75%and total flow rate 5 m3/h,which provided insights to understand that the thickness on pipe wall was roughly uniformly distributed locates near the top of the pipe and the nearer the position gets close to two points,where the oil-water interface contacts the inner wall,the deposition thickness quickly dropped to 0.It was attributed to the fact that a roughly uniformly thickness far away from the oil-water interface contact the inner wall resulted in the slowly changes temperature along the circumferential pipe wall wetted by oil.

Study on two-phase cloud dispersion from liquefied CO_(2) release

The knowledge of two-phase cloud dispersion mechanism from HLG(hazardous liquefied gas) release is the prerequisite for accurate assessment and precise rescue of such accidents. In this paper, an experiment of two-phase cloud dispersion from liquefied CO_(2) hole release is performed. The source terms, such as vapour mass fraction, release velocity and mean droplet diameter, are calculated based on thermodynamic theory. Taking phase transition of CO_(2) droplets to gas into account, CFD(computational fluid dynamics) model for two-phase cloud dispersion is established. The predicted cloud temperatures at the downstream agree well with the experimental data, with the maximum relative error of 5.8% and average relative error of 2.3%. The consequence distances in the downstream direction and in the crosswise direction calculated through two-phase model are larger than those through single-phase model,with the relative differences of 57.8% and 53.6% respectively. CO_(2) concentration calculated by twophase model is smaller in the vicinity of release hole, and larger beyond 0.135 m downstream. A smaller leakage rate results in a lower CO_(2) concentration and a higher cloud temperature.

Preparation of Hydrophobic Fly Ash by Surface Modification and Oil-water Separation Devices

A simple and effective superhydrophobic mesh was designed and made to separate oil-water mixture. Alkali-activated fly ash reacted with 1-bromooctadecane to prepare superhydrophobic modified fly ash (MFA) with low surface energy through Williamson ether synthesis. The MFA powder was then coated uniformly on a stainless steel mesh (SSM) along with the epoxy resin E44 and curing agent T31 to give the superhydrophobic MFA-modified stainless steel mesh (MFA-SSM). The MFA-SSM has a high static water contact angle (CA) of 150.1°and can separate various oil or organic solvent from water with>95%separation efficiency. The oil-water separation efficiency remained high after 30 runs of petroleum ether/water separation.The developed superhydrophobic stainless steel mesh is expected to have wider use in oil-water separation.

Experimental and Numerical Analysis of Oil-Water Flow with Drag Reducing Polymers in Horizontal Pipes

The well-known frictional effect related to liquid-liquid two-phaseflow in pipelines can be reduced using drag-reducing additives.In this study,such an effect has been investigated experimentally using a mixture of oil and water.Moreover,numerical simulations have been carried out using the COMSOL simulation software.The mea-surements were taken in a horizontal pipe with the length and diameter equal to 3 and 0.125 m,respectively.Moreover,Polyethylene oxide with 150 ppm was exploited to reduce the drag effect while considering different water-to-oil fractions(0.3,0.4,0.5,and 0.7)and a constant totalflow velocity of 2.3 m/s.As made evident by the results,a significant reduction can be obtained in terms of pressure drop,which becomes even more significant as the water to oil fraction is increased.The maximum achieved drag reduction is 70%with a water fraction of 0.7.The results also show that the addition of polymer additives can also have an impact on theflow pattern.Com-parison of experimental and numerically determined pressure drop indicates that the error is smaller than 7%.

Liquid-liquid two-phase flow in a wire-embedded concentric microchannel: Flow pattern and mass transfer performance

In this work, flow pattern and mass transfer of liquid-liquid two-phase flow in a wire-embedded concentric microchannel are studied using toluene-water system. Droplet flow, slug flow, oval flow and annular flow are observed in the wire-embedded concentric microchannel. The effects of embedded wires and physical properties on flow patterns are investigated. The embedded wire insert is conducive to the formation of annular flow. The flow pattern distribution regions are distinguished by the Caaq(capillary number)±We_(org)(Weber number) flow pattern map. When Weorg<0.001, slug flow is the main flow pattern, and when Weorg>0.1, annular flow is the main flow pattern. Oval flow and droplet flow are between We_(org)= 0.001-0.1, and oval flow is transformed into droplet flow with the increase of Caaq. The effect of flow rate, phase ratio, initial acetic acid concentration, insert shape and flow patterns on mass transfers are studied. Mass transfer process is enhanced under annular flow conditions, the volumetric mass transfer coefficient is up to 0.36 s^(-1) because of the high interfacial area and interface renewal rate of annular flow.

Simulation of gas–liquid two-phase flow in a flow-focusing microchannel with the lattice Boltzmann method

A lattice Boltzmann method for gas–liquid two-phase flow involving non-Newtonian fluids is developed. Bubble formation in a flow-focusing microchannel is simulated by the method. The influences of flow rate ratio, surface tension,wetting properties, and rheological characteristics of the fluid on the two-phase flow are analyzed. The results indicate that the flow pattern transfers from slug flow to dry-plug flow with a sufficiently small capillary number. Due to the presence of three-phase contact lines, the contact angle has a more significant effect on the dry-plug flow pattern than on the slug flow pattern. The deformation of the front and rear meniscus of a bubble in the shear-thinning fluid can be explained by the variation of the capillary number. The reduced viscosity and increased contact angle are beneficial for the drag reduction in a microchannel. It also demonstrates the effectiveness of the current method to simulate the gas–liquid two-phase flow in a microchannel.