Well production optimization is a complex and time-consuming task in the oilfield development.The combination of reservoir numerical simulator with optimization algorithms is usually used to optimize well production.T...Well production optimization is a complex and time-consuming task in the oilfield development.The combination of reservoir numerical simulator with optimization algorithms is usually used to optimize well production.This method spends most of computing time in objective function evaluation by reservoir numerical simulator which limits its optimization efficiency.To improve optimization efficiency,a well production optimization method using streamline features-based objective function and Bayesian adaptive direct search optimization(BADS)algorithm is established.This new objective function,which represents the water flooding potential,is extracted from streamline features.It only needs to call the streamline simulator to run one time step,instead of calling the simulator to calculate the target value at the end of development,which greatly reduces the running time of the simulator.Then the well production optimization model is established and solved by the BADS algorithm.The feasibility of the new objective function and the efficiency of this optimization method are verified by three examples.Results demonstrate that the new objective function is positively correlated with the cumulative oil production.And the BADS algorithm is superior to other common algorithms in convergence speed,solution stability and optimization accuracy.Besides,this method can significantly accelerate the speed of well production optimization process compared with the objective function calculated by other conventional methods.It can provide a more effective basis for determining the optimal well production for actual oilfield development.展开更多
The aim of this investigation is to determine the effect of fluid leak-off (suction) and fluid injection (blowing) at the horizontal base on the two-dimensional spreading under the gravity of a thin film of viscous in...The aim of this investigation is to determine the effect of fluid leak-off (suction) and fluid injection (blowing) at the horizontal base on the two-dimensional spreading under the gravity of a thin film of viscous incompressible fluid by studying the evolution of the streamlines in the thin film. It is assumed that the normal component of the fluid velocity at the base is proportional to the spatial gradient of the height of the film. Lie symmetry methods for partial differential equations are applied. The invariant solution for the surface profile is derived. It is found that the thin fluid film approximation is satisfied for weak to moderate leak-off and for the whole range of fluid injection. The streamlines are derived and plotted by solving a cubic equation numerically. For fluid injection, there is a dividing streamline originating at the stagnation point at the base which separates the flow into two regions, a lower region consisting mainly of rising fluid and an upper region consisting mainly of descending fluid. An approximate analytical solution for the dividing streamline is derived. It generates an approximate V-shaped surface along the length of the two-dimensional film with the vertex of each section the stagnation point. It is concluded that the fluid flow inside the thin film can be visualised by plotting the streamlines. Other models relating the fluid velocity at the base to the height of the thin film can be expected to contain a dividing streamline originating at a stagnation point and dividing the flow into a lower region of rising fluid and an upper region of descending fluid.展开更多
In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for c...In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for considered schemes are derived in the norm stronger than L^2-norm.展开更多
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.展开更多
The aerodynamic resistance of a train running in the open air limits the maximum speed that can be attained by the train.For this reason,evacuated tube trains(ETT)are being considered as valid alternatives to be imple...The aerodynamic resistance of a train running in the open air limits the maximum speed that can be attained by the train.For this reason,evacuated tube trains(ETT)are being considered as valid alternatives to be implemented in the future.The atmosphere in the tube,the so-called blocking ratio and the length of the streamlined nose are the key factors affecting the aerodynamic performances of these trains.In the present work,we investigate evacuated tube trains with different lengths of the streamline nose on the basis of computational fluid dynamics(CFD).The three-dimensional steady compressible Navier-Stokes equations are solved.The running speed of the ETT is 800 km/h and the blocking ratio is 0.2.Results show that with the increase of the streamlined nose length,the aerodynamic drag and lift forces of the head car decrease gradually,and the drag and lift forces of the middle car change slightly.For the tail car,the drag force decreases,whereas the absolute value of the lift force increases.At a speed of 800 km/h,a slight shock wave appears at the rear of the tail car,which affects the aerodynamic forces acting on the train.展开更多
Streamline-adjustment-assisted heterogeneous combination flooding is a new technology for enhanced oil recovery for post-polymer-flooded reservoirs.In this work,we first carried out a series of 2D visualization experi...Streamline-adjustment-assisted heterogeneous combination flooding is a new technology for enhanced oil recovery for post-polymer-flooded reservoirs.In this work,we first carried out a series of 2D visualization experiments for different chemical flooding scenarios after polymer flooding.Then,we explored the synergistic mechanisms of streamline-adjustment-assisted heterogeneous combination flooding for enhanced oil recovery and the contribution of each component.Test results show that for single heterogeneous combination flooding,the residual oil in the main streamline area after polymer flooding is ready to be driven,but it is difficult to be recovered in the non-main streamline area.Due to the effect of drainage and synergism,the streamline-adjustment-assisted heterogeneous combination flooding diverts the injected chemical agent from the main streamline area to the non-main streamline area,which consequently expands the active area of chemical flooding.Based on the results from the single-factor contribution of the quantitative analysis,the contribution of temporary plugging and profile control of branched preformed particle gels ranks in the first place and followed by the polymer profile control and the effect of streamline adjustment.On the contrary,the surfactant contributes the least to enhance the efficiency of oil displacement.展开更多
As a highly efficient production method, the technique of multi-branch horizontal well is widely used in low permeability reservoirs, heavy oil reservoirs, shallow layer reservoirs and multi-layer reservoirs, because ...As a highly efficient production method, the technique of multi-branch horizontal well is widely used in low permeability reservoirs, heavy oil reservoirs, shallow layer reservoirs and multi-layer reservoirs, because it can significantly improve the productivity of a single well, inhibit coning and enhance oil recovery. Study on sweep efficiency and parameters optimization of multi-branch horizontal well is at the leading edge of research. Therefore, the study is important for enhancing oil recovery and integral exploitation benefit of oil fields. In many applications, streamline simulation shows particular advantages over finite-difference simulation. With the advantages of streamline simulation such as its ability to display paths of fluid flow and acceleration factor in simulation, the flooding process is more visual. The communication between wells and flooding area has been represented appropriately. This method has been applied to the XS9 reservoir in Daqing Oilfield. The production history of this reservoir is about 10 years. The reservoir is maintained above bubble point so that the simulation meets the slight compressibility assumption. New horizontal wells are drilled following this rule.展开更多
The two-dimensional spreading under gravity of a thin fluid film with suction (fluid leak-off) or blowing (fluid injection) at the base is considered. The thin fluid film approximation is imposed. The height of the th...The two-dimensional spreading under gravity of a thin fluid film with suction (fluid leak-off) or blowing (fluid injection) at the base is considered. The thin fluid film approximation is imposed. The height of the thin film satisfies a nonlinear diffusion equation with a source/sink term. The Lie point symmetries of the nonlinear diffusion equation are derived and exist, which provided the fluid velocity at the base, <em>v<sub>n</sub></em> satisfies a first order linear partial differential equation. The general form has algebraic time dependence while a special case has exponential time dependence. The solution in which <em>v<sub>n</sub></em> is proportional to the height of the thin film is studied. The width of the base always increases with time even for suction while the height decreases with time for sufficiently weak blowing. The streamlines of the fluid flow inside the thin film are plotted by first solving a cubic equation. For sufficiently weak blowing there is a dividing streamline, emanating from the stagnation point on the centre line which separates the fluid flow into two regions, a lower region consisting of rising fluid and dominated by fluid injection at the base and an upper region consisting of descending fluid and dominated by spreading due to gravity. For sufficiently strong blowing the lower region expands to completely fill the whole thin film.展开更多
Allen and Liu (1995) introduced a new method for a time-dependent convec-tion dominated diffusion problem, which combines the modified method of characteristics and method of streamline diffusion. But they ignored the...Allen and Liu (1995) introduced a new method for a time-dependent convec-tion dominated diffusion problem, which combines the modified method of characteristics and method of streamline diffusion. But they ignored the fact that the accuracy of time discretization decays at half an order when the characteristic line goes out of the domain. In present paper, the author shows that, as a remedy, a simple lumped scheme yields a full accuracy approximation. Forthermore, some local error bounds independent of the small viscosity are derived for this scheme outside the boundary layers.展开更多
Vortices that develop over intakes are a hazardous hydraulic phenomenon.In this study, a 3D model was developed to study the flow field in air-core vortices.This model is based on the spiral pattern of streamlines and...Vortices that develop over intakes are a hazardous hydraulic phenomenon.In this study, a 3D model was developed to study the flow field in air-core vortices.This model is based on the spiral pattern of streamlines and the analytical solution of the momentum and continuity equations for deriving the three components of velocity.The model provides equations for free surface profiles and 3D patterns of the streamlines.Moreover, a new relationship was suggested for calculating effective viscosity and its distribution across the vortex flow field.The performance of the proposed analytical model was compared with existing experimental data and the results of previous analytical models.The outcomes indicated that the proposed model could predict characteristics of the vortex flow with good accuracy.展开更多
The peristaltic flow of a heated Jeffrey fluid inside a duct with an elliptic cross-section is studied.A thorough heat transfer mechanism is interpreted by analyzing the viscous effects in the energy equation.The gove...The peristaltic flow of a heated Jeffrey fluid inside a duct with an elliptic cross-section is studied.A thorough heat transfer mechanism is interpreted by analyzing the viscous effects in the energy equation.The governing mathematical equations give dimensionless partial differential equations after simplification.The final simplified form of the mathematical equations is evaluated with respect to the relevant boundary conditions,and the exact solution is attained.The results are further illustrated by graphs,and the distinct aspects of peristaltic flow phenomena are discussed.展开更多
This paper aims to present a new streamline diffusion method with low order rectangular Bernardi-Raugel elements to solve the generalized Oseen equations. With the help of the Bramble-Hilbert lemma, the optimal errors...This paper aims to present a new streamline diffusion method with low order rectangular Bernardi-Raugel elements to solve the generalized Oseen equations. With the help of the Bramble-Hilbert lemma, the optimal errors of the velocity and pressure are estimated, which are independent of the considered parameter ε. With an interpolation postprocessing approach, the superconvergent error of the pressure is obtained. Finally,a numerical experiment is carried out to confirm the theoretical results.展开更多
Fluid dynamics simulation is often repeated under varying conditions.This leads to a generation of large amounts of results,which are difficult to compare.To compare results under different conditions,it is effective ...Fluid dynamics simulation is often repeated under varying conditions.This leads to a generation of large amounts of results,which are difficult to compare.To compare results under different conditions,it is effective to overlap the streamlines generated from each condition in a single three-dimensional space.Streamline is a curved line,which represents a wind flow.This paper presents a technique to automatically select and visualize important streamlines that are suitable for the comparison of the simulation results.Additionally,we present an implementation to observe the flow fields in virtual reality spaces.展开更多
In this article, we use streamline diffusion method for the linear second order hyperbolic initial-boundary value problem. More specifically, we prove a posteriori error estimates for this method for the linear wave e...In this article, we use streamline diffusion method for the linear second order hyperbolic initial-boundary value problem. More specifically, we prove a posteriori error estimates for this method for the linear wave equation. We observe that this error estimates make finite element method increasingly powerful rather than other methods.展开更多
To better understand the roles natural fractures play in porous media, an embedded discrete fracture model and streamline modeling method were combined to model natural fractures and compute the flow trajectory and ti...To better understand the roles natural fractures play in porous media, an embedded discrete fracture model and streamline modeling method were combined to model natural fractures and compute the flow trajectory and time of fluid in matrix and fractures systems. The effects of fracture conductivity, number of fractures and fracture locations on fluid flow trajectory and time were examined through analyzing the differences in water breakthrough time and sweeping volume of reservoirs with different fracture networks. When other conditions are the same, compared with homogeneous reservoir without fractures, the fractured reservoir has water breakthrough time 30% sooner and swept volume 10% smaller. Although increase of single fracture can lead to faster water breakthrough and smaller swept volume, adding more fractures wouldn't necessarily reach the same effect. The effect of water flooding is also related to the strike and position of fractures. Fractures in different strikes and positions can result in 20% discrepancy in water breakthrough time and 9% gap in swept volume. The shorter the fracture, the less its effect on fluid flow trajectory and time will be. The position of fracture has a strong influence on sweeping efficiency, and the change of one fracture position could bring about 1% variation in swept volume.展开更多
基金supported partly by the National Science and Technology Major Project of China(Grant No.2016ZX05025-001006)Major Science and Technology Project of CNPC(Grant No.ZD2019-183-007)
文摘Well production optimization is a complex and time-consuming task in the oilfield development.The combination of reservoir numerical simulator with optimization algorithms is usually used to optimize well production.This method spends most of computing time in objective function evaluation by reservoir numerical simulator which limits its optimization efficiency.To improve optimization efficiency,a well production optimization method using streamline features-based objective function and Bayesian adaptive direct search optimization(BADS)algorithm is established.This new objective function,which represents the water flooding potential,is extracted from streamline features.It only needs to call the streamline simulator to run one time step,instead of calling the simulator to calculate the target value at the end of development,which greatly reduces the running time of the simulator.Then the well production optimization model is established and solved by the BADS algorithm.The feasibility of the new objective function and the efficiency of this optimization method are verified by three examples.Results demonstrate that the new objective function is positively correlated with the cumulative oil production.And the BADS algorithm is superior to other common algorithms in convergence speed,solution stability and optimization accuracy.Besides,this method can significantly accelerate the speed of well production optimization process compared with the objective function calculated by other conventional methods.It can provide a more effective basis for determining the optimal well production for actual oilfield development.
文摘The aim of this investigation is to determine the effect of fluid leak-off (suction) and fluid injection (blowing) at the horizontal base on the two-dimensional spreading under the gravity of a thin film of viscous incompressible fluid by studying the evolution of the streamlines in the thin film. It is assumed that the normal component of the fluid velocity at the base is proportional to the spatial gradient of the height of the film. Lie symmetry methods for partial differential equations are applied. The invariant solution for the surface profile is derived. It is found that the thin fluid film approximation is satisfied for weak to moderate leak-off and for the whole range of fluid injection. The streamlines are derived and plotted by solving a cubic equation numerically. For fluid injection, there is a dividing streamline originating at the stagnation point at the base which separates the flow into two regions, a lower region consisting mainly of rising fluid and an upper region consisting mainly of descending fluid. An approximate analytical solution for the dividing streamline is derived. It generates an approximate V-shaped surface along the length of the two-dimensional film with the vertex of each section the stagnation point. It is concluded that the fluid flow inside the thin film can be visualised by plotting the streamlines. Other models relating the fluid velocity at the base to the height of the thin film can be expected to contain a dividing streamline originating at a stagnation point and dividing the flow into a lower region of rising fluid and an upper region of descending fluid.
基金Project supported by National Natural Science Foundation of China and China State Key project for Basic Researchcs.
文摘In this paper, two finite difference streamline diffusion (FDSD) schemes for solving two-dimensional time-dependent convection-diffusion equations are constructed. Stability and optimal order error estimati-ions for considered schemes are derived in the norm stronger than L^2-norm.
文摘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 Sichuan Science and Technology Program(No.2019YJ0227)China Postdoctoral Science Foundation(No.2019M663550)+1 种基金National Natural Science Foundation of China(No.51605397)Science and Technolgoy program of China Railway Group Limited(No.2018-S-02).
文摘The aerodynamic resistance of a train running in the open air limits the maximum speed that can be attained by the train.For this reason,evacuated tube trains(ETT)are being considered as valid alternatives to be implemented in the future.The atmosphere in the tube,the so-called blocking ratio and the length of the streamlined nose are the key factors affecting the aerodynamic performances of these trains.In the present work,we investigate evacuated tube trains with different lengths of the streamline nose on the basis of computational fluid dynamics(CFD).The three-dimensional steady compressible Navier-Stokes equations are solved.The running speed of the ETT is 800 km/h and the blocking ratio is 0.2.Results show that with the increase of the streamlined nose length,the aerodynamic drag and lift forces of the head car decrease gradually,and the drag and lift forces of the middle car change slightly.For the tail car,the drag force decreases,whereas the absolute value of the lift force increases.At a speed of 800 km/h,a slight shock wave appears at the rear of the tail car,which affects the aerodynamic forces acting on the train.
基金financial support from the National Natural Science Foundation of China (Grant No. 51574269)the National Science Foundation for Distinguished Young Scholars of China (Grant No. 51625403)+3 种基金the Important National Science and Technology Specific Projects of China (Grant No. 2016ZX05025-003)the Fundamental Research Funds for the Central Universities (Grant No. 15CX08004A, 18CX02169A)China Postdoctoral Science Foundation (Grant No. 2017M622319)the Natural Science Foundation of Shandong Province (Grant No. ZR2018BEE004)
文摘Streamline-adjustment-assisted heterogeneous combination flooding is a new technology for enhanced oil recovery for post-polymer-flooded reservoirs.In this work,we first carried out a series of 2D visualization experiments for different chemical flooding scenarios after polymer flooding.Then,we explored the synergistic mechanisms of streamline-adjustment-assisted heterogeneous combination flooding for enhanced oil recovery and the contribution of each component.Test results show that for single heterogeneous combination flooding,the residual oil in the main streamline area after polymer flooding is ready to be driven,but it is difficult to be recovered in the non-main streamline area.Due to the effect of drainage and synergism,the streamline-adjustment-assisted heterogeneous combination flooding diverts the injected chemical agent from the main streamline area to the non-main streamline area,which consequently expands the active area of chemical flooding.Based on the results from the single-factor contribution of the quantitative analysis,the contribution of temporary plugging and profile control of branched preformed particle gels ranks in the first place and followed by the polymer profile control and the effect of streamline adjustment.On the contrary,the surfactant contributes the least to enhance the efficiency of oil displacement.
文摘As a highly efficient production method, the technique of multi-branch horizontal well is widely used in low permeability reservoirs, heavy oil reservoirs, shallow layer reservoirs and multi-layer reservoirs, because it can significantly improve the productivity of a single well, inhibit coning and enhance oil recovery. Study on sweep efficiency and parameters optimization of multi-branch horizontal well is at the leading edge of research. Therefore, the study is important for enhancing oil recovery and integral exploitation benefit of oil fields. In many applications, streamline simulation shows particular advantages over finite-difference simulation. With the advantages of streamline simulation such as its ability to display paths of fluid flow and acceleration factor in simulation, the flooding process is more visual. The communication between wells and flooding area has been represented appropriately. This method has been applied to the XS9 reservoir in Daqing Oilfield. The production history of this reservoir is about 10 years. The reservoir is maintained above bubble point so that the simulation meets the slight compressibility assumption. New horizontal wells are drilled following this rule.
文摘The two-dimensional spreading under gravity of a thin fluid film with suction (fluid leak-off) or blowing (fluid injection) at the base is considered. The thin fluid film approximation is imposed. The height of the thin film satisfies a nonlinear diffusion equation with a source/sink term. The Lie point symmetries of the nonlinear diffusion equation are derived and exist, which provided the fluid velocity at the base, <em>v<sub>n</sub></em> satisfies a first order linear partial differential equation. The general form has algebraic time dependence while a special case has exponential time dependence. The solution in which <em>v<sub>n</sub></em> is proportional to the height of the thin film is studied. The width of the base always increases with time even for suction while the height decreases with time for sufficiently weak blowing. The streamlines of the fluid flow inside the thin film are plotted by first solving a cubic equation. For sufficiently weak blowing there is a dividing streamline, emanating from the stagnation point on the centre line which separates the fluid flow into two regions, a lower region consisting of rising fluid and dominated by fluid injection at the base and an upper region consisting of descending fluid and dominated by spreading due to gravity. For sufficiently strong blowing the lower region expands to completely fill the whole thin film.
文摘Allen and Liu (1995) introduced a new method for a time-dependent convec-tion dominated diffusion problem, which combines the modified method of characteristics and method of streamline diffusion. But they ignored the fact that the accuracy of time discretization decays at half an order when the characteristic line goes out of the domain. In present paper, the author shows that, as a remedy, a simple lumped scheme yields a full accuracy approximation. Forthermore, some local error bounds independent of the small viscosity are derived for this scheme outside the boundary layers.
基金supported by the Iran National Science Foundation(INSF,Grant No.97008045)
文摘Vortices that develop over intakes are a hazardous hydraulic phenomenon.In this study, a 3D model was developed to study the flow field in air-core vortices.This model is based on the spiral pattern of streamlines and the analytical solution of the momentum and continuity equations for deriving the three components of velocity.The model provides equations for free surface profiles and 3D patterns of the streamlines.Moreover, a new relationship was suggested for calculating effective viscosity and its distribution across the vortex flow field.The performance of the proposed analytical model was compared with existing experimental data and the results of previous analytical models.The outcomes indicated that the proposed model could predict characteristics of the vortex flow with good accuracy.
文摘The peristaltic flow of a heated Jeffrey fluid inside a duct with an elliptic cross-section is studied.A thorough heat transfer mechanism is interpreted by analyzing the viscous effects in the energy equation.The governing mathematical equations give dimensionless partial differential equations after simplification.The final simplified form of the mathematical equations is evaluated with respect to the relevant boundary conditions,and the exact solution is attained.The results are further illustrated by graphs,and the distinct aspects of peristaltic flow phenomena are discussed.
基金supported by the National Natural Science Foundation of China(Nos.11271340 and11671369)
文摘This paper aims to present a new streamline diffusion method with low order rectangular Bernardi-Raugel elements to solve the generalized Oseen equations. With the help of the Bramble-Hilbert lemma, the optimal errors of the velocity and pressure are estimated, which are independent of the considered parameter ε. With an interpolation postprocessing approach, the superconvergent error of the pressure is obtained. Finally,a numerical experiment is carried out to confirm the theoretical results.
文摘Fluid dynamics simulation is often repeated under varying conditions.This leads to a generation of large amounts of results,which are difficult to compare.To compare results under different conditions,it is effective to overlap the streamlines generated from each condition in a single three-dimensional space.Streamline is a curved line,which represents a wind flow.This paper presents a technique to automatically select and visualize important streamlines that are suitable for the comparison of the simulation results.Additionally,we present an implementation to observe the flow fields in virtual reality spaces.
文摘In this article, we use streamline diffusion method for the linear second order hyperbolic initial-boundary value problem. More specifically, we prove a posteriori error estimates for this method for the linear wave equation. We observe that this error estimates make finite element method increasingly powerful rather than other methods.
文摘To better understand the roles natural fractures play in porous media, an embedded discrete fracture model and streamline modeling method were combined to model natural fractures and compute the flow trajectory and time of fluid in matrix and fractures systems. The effects of fracture conductivity, number of fractures and fracture locations on fluid flow trajectory and time were examined through analyzing the differences in water breakthrough time and sweeping volume of reservoirs with different fracture networks. When other conditions are the same, compared with homogeneous reservoir without fractures, the fractured reservoir has water breakthrough time 30% sooner and swept volume 10% smaller. Although increase of single fracture can lead to faster water breakthrough and smaller swept volume, adding more fractures wouldn't necessarily reach the same effect. The effect of water flooding is also related to the strike and position of fractures. Fractures in different strikes and positions can result in 20% discrepancy in water breakthrough time and 9% gap in swept volume. The shorter the fracture, the less its effect on fluid flow trajectory and time will be. The position of fracture has a strong influence on sweeping efficiency, and the change of one fracture position could bring about 1% variation in swept volume.