Effect of quadratic pressure gradient term on a one-dimensional moving boundary problem based on modified Darcy's law

A relatively high formation pressure gradient can exist in seepage flow in low-permeable porous media with a threshold pressure gradient, and a significant error may then be caused in the model computation by neglecting the quadratic pressure gradient term in the governing equations. Based on these concerns, in consideration of the quadratic pressure gradient term, a basic moving boundary model is constructed for a one-dimensional seepage flow problem with a threshold pressure gradient. Owing to a strong nonlinearity and the existing moving boundary in the mathematical model, a corresponding numerical solution method is presented. First, a spatial coordinate transformation method is adopted in order to transform the system of partial differential equa- tions with moving boundary conditions into a closed system with fixed boundary conditions; then the solution can be sta- bly numerically obtained by a fully implicit finite-difference method. The validity of the numerical method is verified by a published exact analytical solution. Furthermore, to compare with Darcy＇s flow problem, the exact analytical solution for the case of Darcy＇s flow considering the quadratic pressure gradient term is also derived by an inverse Laplace transform. A comparison of these model solutions leads to the conclu- sion that such moving boundary problems must incorporate the quadratic pressure gradient term in their governing equa- tions; the sensitive effects of the quadratic pressure gradient term tend to diminish, with the dimensionless threshold pres- sure gradient increasing for the one-dimensional problem.

Numerical investigation of dual-porosity model with transient transfer function based on discrete-fracture model

Based on the characteristics of fractures in naturally fractured reservoir and a discrete-fracture model, a fracture network numerical well test model is developed. Bottom hole pressure response curves and the pressure field are obtained by solving the model equations with the finite-element method. By analyzing bottom hole pressure curves and the fluid flow in the pressure field, seven flow stages can be recognized on the curves. An upscaling method is developed to compare with the dual-porosity model （DPM）. The comparisons results show that the DPM overestimates the inter-porosity coefficient ）, and the storage factor w. The analysis results show that fracture conductivity plays a leading role in the fluid flow. Matrix permeability influences the beginning time of flow from the matrix to fractures. Fractures density is another important parameter controlling the flow. The fracture linear flow is hidden under the large fracture density. The pressure propagation is slower in the direction of larger fracture density.

Dissipative particle dynamics simulation of wettability alternation phenomena in the chemical flooding process

Wettability alternation phenomena is considered one of the most important enhanced oil recovery （EOR） mechanisms in the chemical flooding process and induced by the adsorption of surfactant on the rock surface. These phenomena are studied by a mesoscopic method named as dissipative particle dynamics （DPD）. Both the alteration phenomena of water-wet to oil-wet and that of oil-wet to water-wet are simulated based on reasonable definition of interaction parameters between beads. The wetting hysteresis phenomenon and the process of oil-drops detachment from rock surfaces with different wettability are simulated by adding long-range external forces on the fluid particles. The simulation results show that, the oil drop is liable to spread on the oil-wetting surface and move in the form of liquid film flow, whereas it is likely to move as a whole on the waterwetting surface. There are the same phenomena occuring in wettability-alternated cases. The results also show that DPD method provides a feasible approach to the problems of seepage flow with physicochemical phenomena and can be used to study the mechanism of EOR of chemical flooding.

A numerical approach for pressure transient analysis of a vertical well with complex fractures

A new well test model for a vertical fractured well is developed based on a discrete-fracture model in which the fractures are discretized as one dimensional（1-D） entities.The model overcomes the weakness of complex meshing,a large number of grids, and instability in conventional stripe-fracture models. Then, the discrete-fracture model is implemented using a hybrid element finite-element method.Triangular elements are used for matrix and line elements for the fractures. The finite element formulation is validated by comparing with the semi-analytical solution of a single vertical fractured well. The accuracy of the approach is shown through several examples with different fracture apertures,fracture conductivity, and fracture amount. Results from the discrete-fracture model agree reasonably well with the stripefracture model and the analytic solutions. The advantages of the discrete-fracture model are presented in mesh generation, computational improvement, and abilities to handle complex fractures like wedge-shaped fractures and fractures with branches. Analytical results show that the number of grids in the discrete-fracture model is 10 % less than stripefracture model, and computational efficiency increases by about 50 %. The more fractures there are, the more the computational efficiency increases.