Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the ef...Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the effects of complex pore structures and wettability.To address this issue,based on the digital rock of low permeability sandstone,a direct numerical simulation is performed considering the interphase drag and boundary slip to clarify the microscopic water-oil displacement process.In addition,a dual-porosity pore network model(PNM)is constructed to obtain the water-oil relative permeability of the sample.The displacement efficiency as a recovery process is assessed under different wetting and pore structure properties.Results show that microscopic displacement mechanisms explain the corresponding macroscopic relative permeability.The injected water breaks through the outlet earlier with a large mass flow,while thick oil films exist in rough hydrophobic surfaces and poorly connected pores.The variation of water-oil relative permeability is significant,and residual oil saturation is high in the oil-wet system.The flooding is extensive,and the residual oil is trapped in complex pore networks for hydrophilic pore surfaces;thus,water relative permeability is lower in the water-wet system.While the displacement efficiency is the worst in mixed-wetting systems for poor water connectivity.Microporosity negatively correlates with invading oil volume fraction due to strong capillary resistance,and a large microporosity corresponds to low residual oil saturation.This work provides insights into the water-oil flow from different modeling perspectives and helps to optimize the development plan for enhanced recovery.展开更多
Laboratory scale model of DMMBF （dual mixed media biofilter） were designed and installed in AI-Mustansiriya University Environmental Hydraulic Lab. Experiments were conducted using two mixed layers through PVR colum...Laboratory scale model of DMMBF （dual mixed media biofilter） were designed and installed in AI-Mustansiriya University Environmental Hydraulic Lab. Experiments were conducted using two mixed layers through PVR column--2.2 m height and 300 mm diameter. The first mixed media filter of depth 640mm mixed of sand, rice husk and granular activated carbon. The percentage volume mix is 1：1：1. While the other mixed media of depth 740 mm, consisting of coal, crash porcelinaite, rock and granite with equally percentage volume. Fifty samples were collected during the experiments, which was spread over a period of forty two weeks. The obtained results indicate that when the flow loading raised from 0.15 L/min to 2.7 L/rain, the removal efficiency of BOD decreased 8%-11%, and the removal efficiency of COD deceased 3%-4%, while the removal efficiency of turbidity increased with the decreasing of hydraulic loading. The results showed that the removal efficiency of turbidity is more than 95% at the lower discharge （0.15 L/min）. Therefore, infiltration should be conservatively designed using low loading rates.展开更多
The tight oil formation develops with microfractures and matrix pores,it is important to study the influence of matrix physical properties on flow characteristics.At first,the representative fracture and matrix sample...The tight oil formation develops with microfractures and matrix pores,it is important to study the influence of matrix physical properties on flow characteristics.At first,the representative fracture and matrix samples are selected respectively in the dual media,the fracture and matrix digital rocks are constructed with micro-CT scanning at different resolutions,and the corresponding fracture and matrix pore networks are extracted,respectively.Then,the modified integration method is proposed to build the dual network model containing both fracture and matrix pore-throat elements,while the geometric-topological structure equivalent matrix pores are generated to fill in the skeleton domain of fracture network,the constructed dual network could describe the geometric-topological structure characteristics of fracture and matrix pore-throat simultaneously.At last,by adjusting the matrix pore density and the matrix filling domain factor,a series of dual network models are obtained to analyze the influence of matrix physical properties on flow characteristics in dual-media.It can be seen that the matrix system contributes more to the porosity of the dual media and less to the permeability.With the decrease in matrix pore density,the porosity/permeability contributions of matrix system to dual media keep decreasing,but the decrease is not significant,the oil-water co-flow zone decreases and the irreducible water saturation increases,and the saturation interval dominated by the fluid flow in the fracture keeps increasing.With the decrease in matrix filling domain factor,the porosity/permeability contributions of matrix system to dual media decreases,the oil-water co-flow zone increases and the irreducible water saturation decreases,and the saturation interval dominated by the fluid flow in the fracture keeps increasing.The results can be used to explain the dual-media flow pattern under different matrix types and different fracture control volumes during tight oil production.展开更多
Considering the influence of quadratic gradient term and medium deformation on the seepage equation, a well testing interpretation model for low permeability and deformation dual medium reservoirs was derived and esta...Considering the influence of quadratic gradient term and medium deformation on the seepage equation, a well testing interpretation model for low permeability and deformation dual medium reservoirs was derived and established. The difference method was used to solve the problem, and pressure and pressure derivative double logarithmic curves were drawn to analyze the seepage law. The research results indicate that the influence of starting pressure gradient and medium deformation on the pressure characteristic curve is mainly manifested in the middle and late stages. The larger the value, the more obvious the upward warping of the pressure and pressure derivative curve;the parameter characterizing the dual medium is the crossflow coefficient. The channeling coefficient determines the time and location of the appearance of the “concave”. The smaller the value, the later the appearance of the “concave”, and the more to the right of the “concave”.展开更多
Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consum...Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consuming. In streamline method, transport equations are solved on one-dimensional streamlines to reduce the computation time with less memory for simulation. First, pressure equation is solved on an Eulerian grid and streamlines are traced. Defining the "time of flight", saturation equations are mapped and solved on streamlines. Finally, the results are mapped back on Eulerian grid and the process is repeated until the simulation end time. The waterflooding process is considered in a fractured reservoir using the dual porosity model. Afterwards, a computational code is developed to solve the same problem by the IMPES method and the results of streamline simulation are compared to those of the IMPES and a commercial software. Finally, the accuracy and efficiency of streamline simulator for simulation of two-phase flow in fractured reservoirs has been proved.展开更多
基金supported by National Natural Science Foundation of China(Grant No.42172159)Science Foundation of China University of Petroleum,Beijing(Grant No.2462023XKBH002).
文摘Multiphase flow in low permeability porous media is involved in numerous energy and environmental applications.However,a complete description of this process is challenging due to the limited modeling scale and the effects of complex pore structures and wettability.To address this issue,based on the digital rock of low permeability sandstone,a direct numerical simulation is performed considering the interphase drag and boundary slip to clarify the microscopic water-oil displacement process.In addition,a dual-porosity pore network model(PNM)is constructed to obtain the water-oil relative permeability of the sample.The displacement efficiency as a recovery process is assessed under different wetting and pore structure properties.Results show that microscopic displacement mechanisms explain the corresponding macroscopic relative permeability.The injected water breaks through the outlet earlier with a large mass flow,while thick oil films exist in rough hydrophobic surfaces and poorly connected pores.The variation of water-oil relative permeability is significant,and residual oil saturation is high in the oil-wet system.The flooding is extensive,and the residual oil is trapped in complex pore networks for hydrophilic pore surfaces;thus,water relative permeability is lower in the water-wet system.While the displacement efficiency is the worst in mixed-wetting systems for poor water connectivity.Microporosity negatively correlates with invading oil volume fraction due to strong capillary resistance,and a large microporosity corresponds to low residual oil saturation.This work provides insights into the water-oil flow from different modeling perspectives and helps to optimize the development plan for enhanced recovery.
文摘Laboratory scale model of DMMBF （dual mixed media biofilter） were designed and installed in AI-Mustansiriya University Environmental Hydraulic Lab. Experiments were conducted using two mixed layers through PVR column--2.2 m height and 300 mm diameter. The first mixed media filter of depth 640mm mixed of sand, rice husk and granular activated carbon. The percentage volume mix is 1：1：1. While the other mixed media of depth 740 mm, consisting of coal, crash porcelinaite, rock and granite with equally percentage volume. Fifty samples were collected during the experiments, which was spread over a period of forty two weeks. The obtained results indicate that when the flow loading raised from 0.15 L/min to 2.7 L/rain, the removal efficiency of BOD decreased 8%-11%, and the removal efficiency of COD deceased 3%-4%, while the removal efficiency of turbidity increased with the decreasing of hydraulic loading. The results showed that the removal efficiency of turbidity is more than 95% at the lower discharge （0.15 L/min）. Therefore, infiltration should be conservatively designed using low loading rates.
基金This work was supported by National Natural Science Foundation of China(No.51704033,No.51804038)PetroChina Innovation Foundation(No.2018D-5007-0210).
文摘The tight oil formation develops with microfractures and matrix pores,it is important to study the influence of matrix physical properties on flow characteristics.At first,the representative fracture and matrix samples are selected respectively in the dual media,the fracture and matrix digital rocks are constructed with micro-CT scanning at different resolutions,and the corresponding fracture and matrix pore networks are extracted,respectively.Then,the modified integration method is proposed to build the dual network model containing both fracture and matrix pore-throat elements,while the geometric-topological structure equivalent matrix pores are generated to fill in the skeleton domain of fracture network,the constructed dual network could describe the geometric-topological structure characteristics of fracture and matrix pore-throat simultaneously.At last,by adjusting the matrix pore density and the matrix filling domain factor,a series of dual network models are obtained to analyze the influence of matrix physical properties on flow characteristics in dual-media.It can be seen that the matrix system contributes more to the porosity of the dual media and less to the permeability.With the decrease in matrix pore density,the porosity/permeability contributions of matrix system to dual media keep decreasing,but the decrease is not significant,the oil-water co-flow zone decreases and the irreducible water saturation increases,and the saturation interval dominated by the fluid flow in the fracture keeps increasing.With the decrease in matrix filling domain factor,the porosity/permeability contributions of matrix system to dual media decreases,the oil-water co-flow zone increases and the irreducible water saturation decreases,and the saturation interval dominated by the fluid flow in the fracture keeps increasing.The results can be used to explain the dual-media flow pattern under different matrix types and different fracture control volumes during tight oil production.
文摘Considering the influence of quadratic gradient term and medium deformation on the seepage equation, a well testing interpretation model for low permeability and deformation dual medium reservoirs was derived and established. The difference method was used to solve the problem, and pressure and pressure derivative double logarithmic curves were drawn to analyze the seepage law. The research results indicate that the influence of starting pressure gradient and medium deformation on the pressure characteristic curve is mainly manifested in the middle and late stages. The larger the value, the more obvious the upward warping of the pressure and pressure derivative curve;the parameter characterizing the dual medium is the crossflow coefficient. The channeling coefficient determines the time and location of the appearance of the “concave”. The smaller the value, the later the appearance of the “concave”, and the more to the right of the “concave”.
文摘Streamline simulation is developed to simulate waterflooding in fractured reservoirs. Conventional reservoir simulation methods for fluid flow simulation in large and complex reservoirs are very costly and time consuming. In streamline method, transport equations are solved on one-dimensional streamlines to reduce the computation time with less memory for simulation. First, pressure equation is solved on an Eulerian grid and streamlines are traced. Defining the "time of flight", saturation equations are mapped and solved on streamlines. Finally, the results are mapped back on Eulerian grid and the process is repeated until the simulation end time. The waterflooding process is considered in a fractured reservoir using the dual porosity model. Afterwards, a computational code is developed to solve the same problem by the IMPES method and the results of streamline simulation are compared to those of the IMPES and a commercial software. Finally, the accuracy and efficiency of streamline simulator for simulation of two-phase flow in fractured reservoirs has been proved.
