The boundary condition is a crucial factor affecting the permeability variation due to suffusion.An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundar...The boundary condition is a crucial factor affecting the permeability variation due to suffusion.An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundary effect is conducted,where the hydraulic head difference(DH)varies,and the boundary includes non-loss and soil-loss conditions.Soil samples are filled into seven soil storerooms connected in turn.After evaluation,the variation in content of fine sand(ΔR_(f))and the hydraulic conductivity of soils in each storeroom(C_(i))are analyzed.In the non-loss test,the soil sample filling area is divided into runoff,transited,and accumulated areas according to the negative or positive ΔR_(f) values.ΔR_(f) increases from negative to positive along the seepage path,and Ci decreases from runoff area to transited area and then rebounds in accumulated area.In the soil-loss test,all soil sample filling areas belong to the runoff area,where the gentle-loss,strengthened-loss,and alleviated-loss parts are further divided.ΔR_(f) decreases from the gentle-loss part to the strengthened-loss part and then rebounds in the alleviated-loss part,and C_(i) increases and then decreases along the seepage path.The relationship between ΔR_(f) and Ci is different with the boundary condition.Ci exponentially decreases with ΔR_(f) in the non-loss test and increases with ΔR_(f) generally in the soil-loss test.展开更多
Long-term permeability experiments have indicated that sorption-induced swelling can switch from internal to bulk depending on the evolutive sorption status.However,this sorption swelling switch mechanism has not been...Long-term permeability experiments have indicated that sorption-induced swelling can switch from internal to bulk depending on the evolutive sorption status.However,this sorption swelling switch mechanism has not been considered in current analytical permeability models.This study introduces a normalized sorption non-equilibrium index(SNEI)to characterize the sorption status,quantify the dynamical variations of matrix swelling accumulation and internal swelling partition,and formulate the sorption swelling switch model.The incorporation of this index into the extended total effective stress concept leads to an analytical transient coal permeability model.Model results show that the sorption swelling switch itself results in the permeability switch under stress-constrained conditions,while the confined bulk swelling suppresses the permeability recovery to the continuous reduction under displacement-constrained conditions.Model verifications show that current experimental observations correspond to the early stages of the transient process,and they could be extended to the whole process with these models.This study demonstrates the importance of the sorption swelling switch in determining permeability evolution using simple boundary conditions.It provides new insights into experimentally revealing the sorption swelling switch in the future,and underscores the requirement of a rigorous model for complex coupled processes in large-scale coal seams.展开更多
Based on Huang's accurate tri-sectional nonlin- ear kinematic equation (1997), a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeabili...Based on Huang's accurate tri-sectional nonlin- ear kinematic equation (1997), a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction: The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation (PDE) sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient (TPG) has a great effect on the external moving boundary but has little effect on the internal moving boundary.展开更多
Taking low permeability cores of Daqing oilfield for example,the flow characteristics at low velocity were studied with the self-designed micro-flux measuring instrument.Considering the throat distribution and capilla...Taking low permeability cores of Daqing oilfield for example,the flow characteristics at low velocity were studied with the self-designed micro-flux measuring instrument.Considering the throat distribution and capillary model,the thickness of fluid boundary layer under different pressure gradients was calculated,and the mechanism and influencing factors of nonlinear percolation were discussed.The results show that the percolation curve of ultra-low rocks is nonlinear,and apparent permeability is not a constant which increases with pressure gradient.The absorption boundary layer decreases with the increase of pressure gradient,and changes significantly especially in low pressure gradient,which is the essence of nonlinear percolation.The absorption boundary layer is also found to be impacted by the surface property of rocks.展开更多
The semi-analytical solutions to Fredlund and Hasan's one-dimensional (1D) consolidation for unsaturated soils with a semi-permeable drainage boundary are pre- seated. Two variables are introduced to transform the ...The semi-analytical solutions to Fredlund and Hasan's one-dimensional (1D) consolidation for unsaturated soils with a semi-permeable drainage boundary are pre- seated. Two variables are introduced to transform the two coupled governing equations of pore-water and pore-air pressures into an equivalent set of partial differential equations (PDFs), which are easily solved by the Laplace transform method. Then, the pore-water pressure, pore-air pressure, and soil settlement are obtained in the Laplace domain. The Crump method is adopted to perform the inverse Laplace transform in order to obtain the semi-analytical solutions in the time domain. It is shown that the proposed solutions are more applicable to various types of boundary conditions and agree well with the existing solutions from the literature. Several numerical examples are provided to investigate the consolidation behavior of an unsaturated single-layer soil with single, double, mixed, and semi-permeable drainage boundaries. The changes in the pore-air and pore-water pres- sures and the soil settlement with the time factor at different values of the semi-permeable drainage boundary parameters are illustrated. In addition, parametric studies are con- ducted on the pore-air and pore-water pressures at different ratios (the air permeability coefficient to the water permeability coefficient) and depths.展开更多
This study presents a novel approach for simultaneous inversion of the key reservoir parameters like horizontal permeability,vertical permeability,skin,and boundary distances for spatial distribution across the grid c...This study presents a novel approach for simultaneous inversion of the key reservoir parameters like horizontal permeability,vertical permeability,skin,and boundary distances for spatial distribution across the grid cells in a 3D single well reservoir model(SWRM).These parameters are first estimated from the standard pressure transient analysis of well test pressure and rate data,which also act as a priori for the inverse problem.A field-worthy layer cake geological model is prepared based on the prior information obtained from pressure transient analysis,followed by a sequential flow simulation of field well test operation.The simulation results provide the model pressure versus rate data as the synthetic data for this study.A cost function is defined incorporating the well test pressure data and model pressure data,which would determine the convergence.The inversion process is to optimize the spatial distribution of reservoir parameters to minimize the difference between the measured pressure transient data and the modelled one,which is obtained from the multiphase fluid flow simulator that solves the implicit black-oil fluid-flow diffusivity equations at every step.A Gauss-Newton(GN)inversion scheme is used for the inversion.The reliability of inversion results depends on the accuracy of priori reservoir parameters fed to the solver,which can be refined if required through uncertainty parameter optimization(UPO).This approach helps to obtain a faster and reliable update of reservoir parameters in a layer cake homogeneous geomodel,hereby introducing the required heterogeneity.This increases the confidence and reliability of a geomodel,which is further used for various production prediction strategies.展开更多
In the present work, the improved version of the meshless singular boundary method(ISBM) is developed for analyzing the performance of bottom standing submerged permeable breakwaters in regular normally incident waves...In the present work, the improved version of the meshless singular boundary method(ISBM) is developed for analyzing the performance of bottom standing submerged permeable breakwaters in regular normally incident waves and in the proximity of a vertical wall. Both single and dual prismatic breakwaters of rectangular and trapezoidal forms are examined. The physical problem is cast in terms of the Laplace equation governing an irrotational flow and incompressible fluid motion with appropriate mixed type boundary conditions, and solved numerically using the ISBM. To model the permeability of the breakwaters fully absorbing boundary conditions are assumed. Numerical results are presented in terms of hydrodynamic quantities of the reflection coefficients. These are firstly validated against the results of a multi-domain boundary element method(BEM) developed independently for a previous study. The agreement between the results of the two methods is excellent. The coefficients of reflection are then computed and discussed for a variety of structural conditions including the breakwaters height, width, spacing, and absorbing permeability. Effects of the proximity of the vertical plane wall are also investigated. The breakwater's width is found to have only marginal effects compared with its height. Permeability tends to decrease the minimum reflections. These coefficients show periodic variations with the spacing relative to the wavelength. Trapezoidal breakwaters are found to be more cost-effective than the rectangular breakwaters. Dual breakwater systems are confirmed to perform much better than single structures.展开更多
In this paper, the basic equations governing the flow and heat transfer of an incompressible viscous and electrically conducting fluid past a semi-infinite vertical per- meable plate in the form of partial differentia...In this paper, the basic equations governing the flow and heat transfer of an incompressible viscous and electrically conducting fluid past a semi-infinite vertical per- meable plate in the form of partial differential equations are reduced to a set of non-linear ordinary differential equations by applying a suitable similarity transformation. Approx- imate solutions of the transformed equations are obtained by employing the perturbation method for two cases, i.e., small and large values of the suction parameter. From the numerical evaluations of the solution, it can be seen that the velocity field at any point decreases as the values of the magnetic and suction parameters increase. The effect of the magnetic parameter is to increase the thermM boundary layer. It is also found that the velocity and temperature fields decrease with the increase in the sink parameter.展开更多
基金The research work described herein was funded by the National Nature Science Foundation of China(Grant No.41877213).This financial support is gratefully acknowledged.
文摘The boundary condition is a crucial factor affecting the permeability variation due to suffusion.An experimental investigation on the permeability of gap-graded soil due to horizontal suffusion considering the boundary effect is conducted,where the hydraulic head difference(DH)varies,and the boundary includes non-loss and soil-loss conditions.Soil samples are filled into seven soil storerooms connected in turn.After evaluation,the variation in content of fine sand(ΔR_(f))and the hydraulic conductivity of soils in each storeroom(C_(i))are analyzed.In the non-loss test,the soil sample filling area is divided into runoff,transited,and accumulated areas according to the negative or positive ΔR_(f) values.ΔR_(f) increases from negative to positive along the seepage path,and Ci decreases from runoff area to transited area and then rebounds in accumulated area.In the soil-loss test,all soil sample filling areas belong to the runoff area,where the gentle-loss,strengthened-loss,and alleviated-loss parts are further divided.ΔR_(f) decreases from the gentle-loss part to the strengthened-loss part and then rebounds in the alleviated-loss part,and C_(i) increases and then decreases along the seepage path.The relationship between ΔR_(f) and Ci is different with the boundary condition.Ci exponentially decreases with ΔR_(f) in the non-loss test and increases with ΔR_(f) generally in the soil-loss test.
基金supported by the Australian Research Council(Grant No.DP200101293)the National Natural Science Foundation of China(Grant No.42202286)the Zhejiang Collaborative Innovation Center for Prevention and Control of Mountain Geological Hazards(Grant No.PCMGH-2017-Z-02).
文摘Long-term permeability experiments have indicated that sorption-induced swelling can switch from internal to bulk depending on the evolutive sorption status.However,this sorption swelling switch mechanism has not been considered in current analytical permeability models.This study introduces a normalized sorption non-equilibrium index(SNEI)to characterize the sorption status,quantify the dynamical variations of matrix swelling accumulation and internal swelling partition,and formulate the sorption swelling switch model.The incorporation of this index into the extended total effective stress concept leads to an analytical transient coal permeability model.Model results show that the sorption swelling switch itself results in the permeability switch under stress-constrained conditions,while the confined bulk swelling suppresses the permeability recovery to the continuous reduction under displacement-constrained conditions.Model verifications show that current experimental observations correspond to the early stages of the transient process,and they could be extended to the whole process with these models.This study demonstrates the importance of the sorption swelling switch in determining permeability evolution using simple boundary conditions.It provides new insights into experimentally revealing the sorption swelling switch in the future,and underscores the requirement of a rigorous model for complex coupled processes in large-scale coal seams.
基金supported by the National Natural Science Foundation of China(11102237)Program for Changjiang Scholars and Innovative Research Team in University(IRT1294)+1 种基金Specialized Research Fund for the Doctoral Program of Higher Education(20110133120012)China Scholarship Council(CSC)
文摘Based on Huang's accurate tri-sectional nonlin- ear kinematic equation (1997), a dimensionless simplified mathematical model for nonlinear flow in one-dimensional semi-infinite long porous media with low permeability is presented for the case of a constant flow rate on the inner boundary. This model contains double moving boundaries, including an internal moving boundary and an external mov- ing boundary, which are different from the classical Stefan problem in heat conduction: The velocity of the external moving boundary is proportional to the second derivative of the unknown pressure function with respect to the distance parameter on this boundary. Through a similarity transfor- mation, the nonlinear partial differential equation (PDE) sys- tem is transformed into a linear PDE system. Then an ana- lytical solution is obtained for the dimensionless simplified mathematical model. This solution can be used for strictly checking the validity of numerical methods in solving such nonlinear mathematical models for flows in low-permeable porous media for petroleum engineering applications. Finally, through plotted comparison curves from the exact an- alytical solution, the sensitive effects of three characteristic parameters are discussed. It is concluded that with a decrease in the dimensionless critical pressure gradient, the sensi- tive effects of the dimensionless variable on the dimension- less pressure distribution and dimensionless pressure gradi- ent distribution become more serious; with an increase in the dimensionless pseudo threshold pressure gradient, the sensi- tive effects of the dimensionless variable become more serious; the dimensionless threshold pressure gradient (TPG) has a great effect on the external moving boundary but has little effect on the internal moving boundary.
基金Project(2008ZX05013) supported by the National Science and Technology Project of ChinaProject(10672187) supported by the National Natural Science Foundation of China
文摘Taking low permeability cores of Daqing oilfield for example,the flow characteristics at low velocity were studied with the self-designed micro-flux measuring instrument.Considering the throat distribution and capillary model,the thickness of fluid boundary layer under different pressure gradients was calculated,and the mechanism and influencing factors of nonlinear percolation were discussed.The results show that the percolation curve of ultra-low rocks is nonlinear,and apparent permeability is not a constant which increases with pressure gradient.The absorption boundary layer decreases with the increase of pressure gradient,and changes significantly especially in low pressure gradient,which is the essence of nonlinear percolation.The absorption boundary layer is also found to be impacted by the surface property of rocks.
基金Project supported by the National Natural Science Foundation of China(Nos.41630633 and11672172)
文摘The semi-analytical solutions to Fredlund and Hasan's one-dimensional (1D) consolidation for unsaturated soils with a semi-permeable drainage boundary are pre- seated. Two variables are introduced to transform the two coupled governing equations of pore-water and pore-air pressures into an equivalent set of partial differential equations (PDFs), which are easily solved by the Laplace transform method. Then, the pore-water pressure, pore-air pressure, and soil settlement are obtained in the Laplace domain. The Crump method is adopted to perform the inverse Laplace transform in order to obtain the semi-analytical solutions in the time domain. It is shown that the proposed solutions are more applicable to various types of boundary conditions and agree well with the existing solutions from the literature. Several numerical examples are provided to investigate the consolidation behavior of an unsaturated single-layer soil with single, double, mixed, and semi-permeable drainage boundaries. The changes in the pore-air and pore-water pres- sures and the soil settlement with the time factor at different values of the semi-permeable drainage boundary parameters are illustrated. In addition, parametric studies are con- ducted on the pore-air and pore-water pressures at different ratios (the air permeability coefficient to the water permeability coefficient) and depths.
文摘This study presents a novel approach for simultaneous inversion of the key reservoir parameters like horizontal permeability,vertical permeability,skin,and boundary distances for spatial distribution across the grid cells in a 3D single well reservoir model(SWRM).These parameters are first estimated from the standard pressure transient analysis of well test pressure and rate data,which also act as a priori for the inverse problem.A field-worthy layer cake geological model is prepared based on the prior information obtained from pressure transient analysis,followed by a sequential flow simulation of field well test operation.The simulation results provide the model pressure versus rate data as the synthetic data for this study.A cost function is defined incorporating the well test pressure data and model pressure data,which would determine the convergence.The inversion process is to optimize the spatial distribution of reservoir parameters to minimize the difference between the measured pressure transient data and the modelled one,which is obtained from the multiphase fluid flow simulator that solves the implicit black-oil fluid-flow diffusivity equations at every step.A Gauss-Newton(GN)inversion scheme is used for the inversion.The reliability of inversion results depends on the accuracy of priori reservoir parameters fed to the solver,which can be refined if required through uncertainty parameter optimization(UPO).This approach helps to obtain a faster and reliable update of reservoir parameters in a layer cake homogeneous geomodel,hereby introducing the required heterogeneity.This increases the confidence and reliability of a geomodel,which is further used for various production prediction strategies.
基金financially supported by the Direction Général des Enseignements et de la Formation Supérieure of Algeria(Grant CNEPRU No.G0301920140029)
文摘In the present work, the improved version of the meshless singular boundary method(ISBM) is developed for analyzing the performance of bottom standing submerged permeable breakwaters in regular normally incident waves and in the proximity of a vertical wall. Both single and dual prismatic breakwaters of rectangular and trapezoidal forms are examined. The physical problem is cast in terms of the Laplace equation governing an irrotational flow and incompressible fluid motion with appropriate mixed type boundary conditions, and solved numerically using the ISBM. To model the permeability of the breakwaters fully absorbing boundary conditions are assumed. Numerical results are presented in terms of hydrodynamic quantities of the reflection coefficients. These are firstly validated against the results of a multi-domain boundary element method(BEM) developed independently for a previous study. The agreement between the results of the two methods is excellent. The coefficients of reflection are then computed and discussed for a variety of structural conditions including the breakwaters height, width, spacing, and absorbing permeability. Effects of the proximity of the vertical plane wall are also investigated. The breakwater's width is found to have only marginal effects compared with its height. Permeability tends to decrease the minimum reflections. These coefficients show periodic variations with the spacing relative to the wavelength. Trapezoidal breakwaters are found to be more cost-effective than the rectangular breakwaters. Dual breakwater systems are confirmed to perform much better than single structures.
文摘In this paper, the basic equations governing the flow and heat transfer of an incompressible viscous and electrically conducting fluid past a semi-infinite vertical per- meable plate in the form of partial differential equations are reduced to a set of non-linear ordinary differential equations by applying a suitable similarity transformation. Approx- imate solutions of the transformed equations are obtained by employing the perturbation method for two cases, i.e., small and large values of the suction parameter. From the numerical evaluations of the solution, it can be seen that the velocity field at any point decreases as the values of the magnetic and suction parameters increase. The effect of the magnetic parameter is to increase the thermM boundary layer. It is also found that the velocity and temperature fields decrease with the increase in the sink parameter.