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.

Anion Content in Surface and Subsurface Water Flows in Grassland

Grassland devoted about 90% of agricultural Irish land. Anion retention capacity of most agriculture soils is less than cation retention capacity, therefore chemical analyzed for anion （NO3 （nitrate）, CI- （chloride）, SO42 （sulphate） and HPO4 （biphosphate）） in different water flow types （overland flow, interflow and vertical flow or drainage） samples. In this work, simple equipment was used to operationally distinguish between overland flow and interflow while vertical flow collected in different depths using soil water sampler equipments. Episodes of overland flow and interflow occurred even though the site is located in the lowest rainfall in Ireland and on well-drained soil. Samples of different origin showed marked differences in their anion contents, while HPO4 concentrations were almost equal to zero in all water flow types, nitrate where relatively high in overland and interflow samples and was very high in drainage samples. No sampling was carried out below the root zone but it must be assumed that the high concentrations measured in drainage samples would constitute a threat to groundwater resources. When overland flow and interflow did occur, NO3 concentrations were usually close to or in excess of the maximum admissible concentrations for drinking water and it will be constituted a threat to inland surface water bodies.

On Lagrange-Projection Schemes for Shallow Water Flows Over Movable Bottom with Suspended and Bedload Transport

In the present work we aim to simulate shallow water flows over movable bottom with suspended and bedload transport.In order to numerically approximate such a system,we proceed step by step.We start by considering shallow water equations with non-constant density of the mixture water-sediment.Then,the Exner equation is included to take into account bedload sediment transport.Finally,source terms for friction,erosion and deposition processes are considered.Indeed,observe that the sediment particle could go in suspension into the water or being deposited on the bottom.For the numerical scheme,we rely on well-balanced Lagrange-projection methods.In particular,since sediment transport is generally a slow process,we aim to develop semi-implicit schemes in order to obtain fast simulations.The Lagrange-projection splitting is well-suited for such a purpose as it entails a decomposition of the(fast)acoustic waves and the(slow)material waves of the model.Hence,in subsonic regimes,an implicit approximation of the acoustic equations allows us to neglect the corresponding CFL condition and to obtain fast numerical schemes with large time step.

Lattice Boltzmann Simulation of Free-Surface Temperature Dispersion in Shallow Water Flows

We develop a lattice Boltzmann method for modeling free-surface temperature dispersion in the shallow water flows.The governing equations are derived from the incompressible Navier-Stokes equations with assumptions of shallow water flows including bed frictions,eddy viscosity,wind shear stresses and Coriolis forces.The thermal effects are incorporated in the momentum equation by using a Boussinesq approximation.The dispersion of free-surface temperature is modelled by an advection-diffusion equation.Two distribution functions are used in the lattice Boltzmann method to recover the flow and temperature variables using the same lattice structure.Neither upwind discretization procedures nor Riemann problem solvers are needed in discretizing the shallow water equations.In addition,the source terms are straightforwardly included in the model without relying on well-balanced techniques to treat flux gradients and source terms.We validate the model for a class of problems with known analytical solutions and we also present numerical results for sea-surface temperature distribution in the Strait of Gibraltar.

A NEW NUMERICAL MODEL FOR THREE DIMENSIONAL SURFACE WATER FLOWS

A three dimensional numerical model based on the Reynolds equations is presented that can be used to predict the surface water flow in open channels.The model uses a computational mesh that conforms to the free water sur- face and the bottom of the channel so that the accuracy of boundary condition application,code complexity,and e- conomy could be enhanced.The k-ε turbulence model is used to estimate the eddy viscosity coefficient.Instead of using the“rigid-lid”approximation a 2-D equation derived from integrating the continuity equation over the total depth is adopted to determine the elevation of the free water surface.A new algorithm is presented based on the conventional SIMPLE procedure.The block correction technique is employed to enhance rate of convergence. The model presented is applied to a bottom discharge into a rectangular straight channel for three dimensional phenomena to obtain the free water surface configuration,velocities and pressure.The computed results are in good agreement with the previous experimental values.

The Water Vapor Transport Model at the Regional Boundaryduring the Meiyu Period

The water vapor transport model at the regional boundary in the Meiyu period is put forward through diagnostic analysis. The numerical simulation on the water vapor transport at the boundary of China in the heavy rainfall period during June–July 1998 shows that the feature of water vapor transport in June is different from that in July. The main body of the water cycle that forms the torrential rain in the Yangtze River Valley is made up of water vapor transport at the western and southern boundaries of the China region in June, whereas the water vapor flow at the western boundary in middle Tibet turns out to be the main body of water vapor sources in July. The water vapor transport at the western boundary of the Tibetan Plateau and the southern boundary of China plays an important role in the torrential rain in the Yangtze River Valley. The temporal and spatial distribution characteristics of water vapor flow at the regional boundary and their theoretical model would provide the scientific proof for the heavy rain forecasts in the Yangtze River Valley.

Forecasting water disaster for a coal mine under the Xiaolangdi reservoir

Xin’an coal mine, Henan Province, faces the risk of water inrush because 40% of the area of the coal mine is under the surface water of the Xiaolangdi reservoir. To forecast water disaster, an effective aquifuge and a limit of water infiltration were determined by rock-phase analysis and long term observations of surface water and groundwater. By field monitoring, as well as physical and numerical simulation experiments, we obtained data reflecting different heights of a water flow fractured zone (WFFZ) under different mining conditions, derived a formula to calculate this height and built a forecasting model with the aid of GIS. On the basis of these activities, the coal mine area was classified into three sub-areas with different potential of water inrush. In the end, our research results have been applied in and verified by industrial mining experiments at three working faces and we were able to present a successful example of coal mining under a large reservoir.

Three-dimensional hydrodynamic model of Xiamen waters

A semi-implicit and Eulerian - Lagrangian finite difference method for three-dimensionalshallow flow has been extended to a more complete system of equations incorporating second-moment turbulence closure model and transport equations of salinity and temperature. The simulation for flooding and drying of mudflats has been improved. The model is applied to Xiamen waters. Based on extensive survey data, water level elevation, temperature and salinity field along the eastern open boundary and at the Jiulong River inlets and runoffs are analyzed, specified and calibrated. The computed results show good agreement with the measured data, reproduce flooding, emergence of large and complex mudflat region.

Influence of underground water seepage flow on surrounding rock deformation of multi-arch tunnel

Based on a typical multi-arch tunnel in a freeway, the fast Lagrangian analysis of continua in 3 dimensions(FLAC3D) was used to calculate the surrounding rock deformation of the tunnel under which the effect of underground water seepage flow was taken into account or not. The distribution of displacement field around the multi-arch tunnel, which is influenced by the seepage field, was gained. The result indicates that the settlement values of the vault derived from coupling analysis are bigger when considering the seepage flow effect than that not considering. Through the contrast of arch subsidence quantities calculated by two kinds of computation situations, and the comparison between the calculated and measured value of tunnel vault settlement, it is found that the calculated value(5.7-6.0 mm) derived from considering the seepage effect is more close to the measured value(5.8-6.8 mm). Therefore, it is quite necessary to consider the seepage flow effect of the underground water in aquiferous stratum for multi-arch tunnel design.

Spatial matching and flow in supply and demand of water provision services: A case study in Xiangjiang River Basin

Global climate change and increased human consumption have aggravated the uneven spatiotemporal distribution of watershed water resources, affecting the water provision supply and demand state. However, this problem has often been ignored. The present study used the Xiangjiang River basin(XRB) as the study area, and the Integrated Valuation of Ecosystem Services and Trade-offs(InVEST) model, demand quantification model,supply–demand ratio, and water flow formula were applied to explore the spatial heterogeneity, flow, and equilibrium between water supply and demand. The results demonstrated significant spatial heterogeneity in the upstream, midstream, and downstream regions.The areas of water shortage were mainly located the downstream of the Changsha–Zhuzhou–Xiangtan urban agglomeration, and the Hengyang basin was the most scarcity area. Affected by terrain gradients and human needs, water flow varied from-16.33 × 10^(8) m^(3) to 13.69 × 10^(8) m^(3)from the upstream to the downstream area, which provided a possibility to reduce spatial heterogeneity. In the future, measures such as strengthening water resource system control,sponge city construction, and dynamic monitoring technology should be taken to balance the supply and demand of water in different river sections of the basin. This study can provide references for regulating water resources allocation in different reaches of the basin.

Characteristics of Underground Water Flow at Different Water Levels in Tianshengan Karst Area, Yunnan, China

Three tracing tests from the same injection point executed at low, medium, and high water levels in the karst aquifer near Tianshengan village, Lunan Stone Forest, Yunnan Province, China, have revealed the basic properties of underground water flow. They showed the general directions of water flows; tracer concentrations were observed at six successive points allowing for the calculation of apparent dominant flow velocities at these sections towards the Dalongtan karst spring. For the high water level, the discharge between single sections was between two and 10 times greater than that at low water level. For the medium water level, the flow velocity at different sections was between 1.4 and 3.7 times faster than that at low water level; and for high water level, it was between 1.3 and 2.7 times faster than that at medium water level. The fastest water flow appeared at the first section （23 cm/s at medium water level）; and the slowest （0.6 cm/s at low water level） appeared where water flow must cross the Tianshengan fault （north-south direction）, and later, a layer of 20-30 m thickness of quartz sandstone and shale clay-stones. It was also possible to calculate the recovery of the tracer for point 4, Dakenyan, where discharge was measured. At the medium water level, 50% of the injected tracer was detected a half-day after its first appearance and at low water level after more than 3 days. The previously published research illustrates the transport velocities of possible contaminants and their solubilities in water at different hydrological conditions.

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.

The effects of water flow and temperature on thermal regime around a culvert built on permafrost

Temperature and water flow through a culvert beneath the Alaska Highway near Beaver Creek,Yukon,were measured at hourly intervals between June and October 2013.These data were used to simulate the effect of the culvert on the thermal regime of the road embankment and subjacent permafrost.A 2-D thermal model of the embankment and permafrost was developed with TEMP/W and calibrated using field observations.Empirical relations were obtained between water temperatures at the entrance to the culvert,flow into the culvert,and water temperatures inside the structure.Water temperatures at the entrance and inside the culvert had a linear relation,while water temperatures inside the culvert and water flow were associated by a logarithmic relation.A multiple linear regression was used to summarize these relations.From this relationship,changes in the flow rate and water temperatures at the entrance of the culvert were simulated to obtain predicted water temperatures in the culvert.The temperatures in the culvert were used in the thermal model to determine their effects on the ground thermal regime near the culvert.Variation of ±10% in water flow rate had no impact on the thermal regime underneath the culvert.Variation of water temperature at the entrance of the culvert had a noticeable influence on the thermal regime.A final simulation was conducted without insulation beneath the culvert.The thaw depth was 30 cm with insulation,and 120 cm without insulation,illustrating the importance of insulation to the ground thermal regime.

Comparative experimental investigation of chemical grouting into a fracture with flowing and static water

We present a series of experimental tests on chemical grouting into a fracture with flowing and static water,using a transparent fracture grouting experimental device.Variations of seepage pressure and grout propagation were compared in our investigation.The results show that flowing water results in drops of seepage pressure,development of penetration radii in the upstream side and drops of propagation area during the same period,compared with grouting in static water.The propagation area in static water is always round before grouts reach the joint boundaries.However,the propagation shape changes from round to an elliptic shape for grouting into a fracture with flowing water.A theoretical model for the grout penetration radius in a fracture considering flowing velocity was developed and validated by our experimental results.These results are helpful in improving understanding of fracture grouting mechanism and in guiding engineering practices.

A New Wetting and Drying Method for Moving Boundary in Shallow Water Flow Models

To deal with the moving boundary hydrodynamic problems of the tidal flats in shallow water flow models, a new wetting and drying （WD） method is proposed. In the new method, a ＂predicted water depth＂ is evaluated explicitly based on the simplified shallow water equations and used to determine the status （wet or dry） together with the direction of flow. Compared with previous WD method, besides the water elevation, more factors, such as the flow velocity and the surface shear stress, are taken into account in the new method to determine the moving boundary. In addition, a formula is deduced to determine the threshold, as critical water depth, which needs to be preset before simulations. The new WD method is tested with five cases including three 1D ones and two 2D ones. The results show that the new WD method can simulate the wetting and drying process, in beth typical and practical cases, with smooth manner and achieves effective estimation of the retention volume at shallow water body.

Simulation of Wetting and Drying Processes in A Depth Integrated Shallow Water Flow Model by Slot Method

A particular porosity method named ＂slot method＂ is implemented in a depth-integrated shallow water flow model （DIVAST） to simulate wetting and drying processes. Discussed is the relationship between the shape factors of the ＂slot＂ and the preset depth used in ＂wetting-drying＂ algorithm. Two typical tests are conducted to examine the performance of the method with the effect of the shape factors of the ＂slot＂ being checked in detail in the first test. Numerical results demonstrate that： 1 ） no additional effort to improve the finite difference scheme is needed to implement ＂slot method＂ in DIVAST, and 2） ＂slot method＂ will simulate wetting and diying processes correctly if the shape factors of the ＂slot＂ being selected properly.

Laboratory investigation into the oil diffusion from submarine pipeline under water flow

A physical model test has been conducted to study the oil diffusion from the submarine pipeline under water flow.The crude oil in the flume is spilled from a leakage point of the pipeline and diffused from the seabed to the surface. By the non-contact optical measuring technology, an image acquisition and data analysis system is designed to explore the spilled mechanism and characteristic. The oil trajectory, velocity and the rising time to the surface are obtained through this system. The influence of the water flow and the spilled discharge on the behavior of the spilled oil are analyzed from both qualitative and quantitative perspectives. The sensitivity study of the characteristic physical quantities to various factors are presented afterward. The spilled oil under water is mainly distributed in the form of the scattered particles with different sizes. The rising process of the oil can be divided into three stages: full, dispersion and aggregation period. The spilled discharge is the primary factor affecting the rising time of the oil particles. In the rising process of the oil particles, the vertical velocity of the oil is mainly affected by the spilled discharge, and the transverse velocity is more dependent on the water velocity. The deviation of the transverse oil velocity is much larger than that of the rising time and the vertical oil velocity. The study can provide a theoretical reference for the prediction system of oil spill emergency.

Wall Sticking of High Water-Cut, Highly Viscous and High Gel-Point Crude Oil Transported at Low Temperatures

Some crude oils with high water cut have the capability to flow below the oil gel point, while the oil particles adhere to the pipe wall in the form of paste, also called "wall sticking". Wall sticking is a serious problem during the pipeline transportation, leading to partial or total blockage of the pipeline and energy wastage. In this paper, a series of laboratory flow loop experiments were conducted to observe the wall sticking characteristics of crude oil with high water cut, high viscosity and high gel point at low transportation temperatures. The effects of shear stress and water cut on the wall sticking rate and occurrence temperature were investigated. Experimental results indicated that the wall sticking rate and occurrence temperature were lower under stronger shear stress and higher water cut conditions. A criterion of wall sticking occurrence temperature(WSOT) and a regression model of wall sticking rate were then established. Finally, the software was developed to calculate the pressure drop along the pipelines of crude oils with high water-cut. It was able to predict the wall sticking thickness of gelled oil and then calculate the pressure drop along the pipelines. A typical case study indicated that the prediction results obtained from the software were in agreement with actual measured values.

Prediction of curved oil–water interface in horizontal pipes using modified model with dynamic contact angle

In this study,interface shapes of horizontal oil–water two-phase flow are predicted by using Young-Laplace equation model and minimum energy model.Meanwhile,the interface shapes of horizontal oil–water twophase flow in a 20 mm inner diameter pipe are measured by a novel conductance parallel-wire array probe(CPAP).It is found that,for flow conditions with low water holdup,there is a large deviation between the model-predicted interface shape and the experimentally measured one.Since the variation of pipe wetting characteristics in the process of fluid flow can lead to the changes of the contact angle between the fluid and the pipe wall,the models mentioned above are modified by considering dynamic contact angle.The results indicate that the interface shapes predicted by the modified models present a good consistence with the ones measured by CPAP.

Water invasion and remaining gas distribution in carbonate gas reservoirs using core displacement and NMR

Water invasion is a common phenomenon in gas reservoirs with active edge-and-bottom aquifers.Due to high reservoir heterogeneity and production parameters,carbonate gas reservoirs feature exploitation obstacles and low recovery factors.In this study,combined core displacement and nuclear magnetic resonance(NMR)experiments explored the reservoir gas−water two-phase flow and remaining microscopic gas distribution during water invasion and gas injection.Consequently,for fracture core,the water-phase relative permeability is higher and the co-seepage interval is narrower than that of three pore cores during water invasion,whereas the water-drive recovery efficiency at different invasion rates is the lowest among all cores.Gas injection is beneficial for reducing water saturation and partially restoring the gas-phase relative permeability,especially for fracture core.The remaining gas distribution and the content are related to the core properties.Compared with pore cores,the water invasion rate strongly influences the residual gas distribution in fracture core.The results enhance the understanding of the water invasion mechanism,gas injection to resume production and the remaining gas distribution,so as to improve the recovery factors of carbonate gas reservoirs.