共查询到20条相似文献，搜索用时 15 毫秒
1.
The matrix–fracture transfer shape factor is one of the important parameters in the modeling of fluid flow in fractured porous media using a dual-porosity concept. Warren and Root [36] introduced the dual-porosity concept and suggested a relation for the shape factor. There is no general relationship for determining the shape factor for a single-phase flow of slightly compressible fluids. Therefore, different studies reported different values for this parameter, as an input into the flow models. Several investigations have been reported on the shape factor for slightly compressible fluids. However, the case of compressible fluids has not been investigated in the past. The focus of this study is, therefore, to find the shape factor for the single-phase flow of compressible fluids (gases) in fractured porous media. In this study, a model for the determination of the shape factor for compressible fluids is presented; and, the solution of nonlinear gas diffusivity equation is used to derive the shape factor. The integral method and the method of moments are used to solve the nonlinear governing equation by considering the pressure dependency of the viscosity and isothermal compressibility of the fluid. The approximate semi-analytical model for the shape factor presented in this study is verified using single-porosity, fine-grid, numerical simulations. The dependency of the shape factor on the gas specific gravity, pressure and temperature are also investigated. The theoretical analysis presented improves our understanding of fluid flow in fractured porous media. In addition, the developed matrix–fracture transfer shape factor can be used as an input for modeling flow of compressible fluids in dual-porosity systems, such as naturally fractured gas reservoirs, coalbed methane reservoirs and fractured tight gas reservoirs. 相似文献
2.
The occurrence of roll-waves in mud-flows is investigated based on the formulation of the marginal stability threshold of a linearized onedimensional viscoplastic (shear-thinning) flow model. Since for this kind of non-Newtonian rheological models this threshold may occur in a hypocritical flow, the downstream boundary condition may have a nonnegligible effect on the spatial growth/decay of the perturbation. The paper presents the solution of the 1D linearized flow of a Herschel and Bulkley fluid in a channel of finite length, in the neighbourhood of a hypocritical base uniform flow. Both linearly stable and unstable conditions are considered. The analytical solution is found applying the Laplace transform method and obtaining the first-order analytical expressions of the upstream and downstream channel response functions in the time domain. The effects of both the yield stress and the rheological law exponent are discussed, recovering as particular cases both power-law and Bingham fluids. The theoretical achievements may be used to extend semi-empirical criteria commonly employed for predicting roll waves occurrence in clear water even to mud-flows. 相似文献
3.
Shifang Wang Boming Yu Qian ZhengYonggang Duan Quantang Fang 《Advances in water resources》2011,34(12):1574-1580
An analytical expression is derived for the starting pressure gradient for Bingham fluids in porous media embedded with randomly distributed fractal-like tree networks based on fractal theory and technique. The proposed model relates the flow rate and the starting pressure gradient to the structural parameters of porous media and microstructural parameters of fractal-like tree networks, the yield stress and fractal dimensions of porous media and maximum mother diameter of randomly distributed fractal-like tree networks. The results show that the starting pressure gradient decreases with the increase of porosity of matrix material, fractal dimension for mother diameters, diameter ratio and permeability, and the starting pressure gradient increases with the increase of the length ratio and the yield stress. The model predictions from the present model for the starting pressure gradient are in good agreement with the available expression. 相似文献
4.
A precise value of the matrix-fracture transfer shape factor is essential for modeling fluid flow in fractured porous media by a dual-porosity approach. The slightly compressible fluid shape factor has been widely investigated in the literature. In a recent study, we have developed a transfer function for flow of a compressible fluid using a constant fracture pressure boundary condition [Ranjbar E, Hassanzadeh H, Matrix-fracture transfer shape factor for modeling flow of a compressible fluid in dual-porosity media. Adv Water Res 2011;34(5):627-39. doi:10.1016/j.advwatres.2011.02.012]. However, for a compressible fluid, the consequence of a pressure depletion boundary condition on the shape factor has not been investigated in the previous studies. The main purpose of this paper is, therefore, to investigate the effect of the fracture pressure depletion regime on the shape factor for single-phase flow of a compressible fluid. In the current study, a model for evaluation of the shape factor is derived using solutions of a nonlinear diffusivity equation subject to different pressure depletion regimes. A combination of the heat integral method, the method of moments and Duhamel’s theorem is used to solve this nonlinear equation. The developed solution is validated by fine-grid numerical simulations. The presented model can recover the shape factor of slightly compressible fluids reported in the literature. This study demonstrates that in the case of a single-phase flow of compressible fluid, the shape factor is a function of the imposed boundary condition in the fracture and its variability with time. It is shown that such dependence can be described by an exponentially declining fracture pressure with different decline exponents. These findings improve our understanding of fluid flow in fractured porous media. 相似文献
5.
《Advances in water resources》2005,28(5):493-506
Upscaling of multi-phase flow problems for a heterogeneous porous medium requires modification of constitutive functions at the grid-block scale. A particular type of heterogeneity that has important environmental consequences involves thin, continuous streaks of high permeability through lower-permeability background rocks. These streaks, which may correspond to features like abandoned wells in mature sedimentary basins, can become preferential flow paths for an invading fluid. Quantification of flow through these types of heterogeneities in deep, geological formations is necessary for estimates of migration and possible leakage of injected fluids such as hazardous liquid wastes, municipal liquid wastes, and, possibly, carbon dioxide. One of the important constitutive functions for proper estimation of flow through these flow paths is the relative permeability function. In the simple case of a single high-permeability streak in a uniform rock matrix, with both materials having identical (local) relative permeability functions, the upscaled relative permeability must be changed significantly to capture the proper leakage. Standard petroleum reservoir pseudo-functions for relative permeability capture the general features of the upscaled function, but they still produce errors of several hundred percent in the leakage estimation. Detailed three-dimensional numerical simulations and associated upscaled calculations demonstrate the proper form for the upscaled relative permeability, and provide a modified derivation of pseudo-functions to capture the leakage behavior in upscaled models. 相似文献
6.
《Advances in water resources》2007,30(4):927-936
A nonlinear model for single-phase fluid flow in slightly compressible porous media is presented and solved approximately. The model assumes state equations for density, porosity, viscosity and permeability that are exponential functions of the fluid (either gas or liquid) pressure. The governing equation is transformed into a nonlinear diffusion equation. It is solved for a semi-infinite domain for either constant pressure or constant flux boundary conditions at the surface. The solutions obtained, although approximate, are extremely accurate as demonstrated by comparisons with numerical results. Predictions for the surface pressure resulting from a constant flux into a porous medium are compared with published experimental data. 相似文献
7.
《Advances in water resources》2004,27(5):547-562
This study proposes the use of several problems of unstable steady state convection with variable fluid density in a porous layer of infinite horizontal extent as two-dimensional (2-D) test cases for density-dependent groundwater flow and solute transport simulators. Unlike existing density-dependent model benchmarks, these problems have well-defined stability criteria that are determined analytically. These analytical stability indicators can be compared with numerical model results to test the ability of a code to accurately simulate buoyancy driven flow and diffusion. The basic analytical solution is for a horizontally infinite fluid-filled porous layer in which fluid density decreases with depth. The proposed test problems include unstable convection in an infinite horizontal box, in a finite horizontal box, and in an infinite inclined box. A dimensionless Rayleigh number incorporating properties of the fluid and the porous media determines the stability of the layer in each case. Testing the ability of numerical codes to match both the critical Rayleigh number at which convection occurs and the wavelength of convection cells is an addition to the benchmark problems currently in use. The proposed test problems are modelled in 2-D using the SUTRA [SUTRA––A model for saturated–unsaturated variable-density ground-water flow with solute or energy transport. US Geological Survey Water-Resources Investigations Report, 02-4231, 2002. 250 p] density-dependent groundwater flow and solute transport code. For the case of an infinite horizontal box, SUTRA results show a distinct change from stable to unstable behaviour around the theoretical critical Rayleigh number of 4π2 and the simulated wavelength of unstable convection agrees with that predicted by the analytical solution. The effects of finite layer aspect ratio and inclination on stability indicators are also tested and numerical results are in excellent agreement with theoretical stability criteria and with numerical results previously reported in traditional fluid mechanics literature. 相似文献
8.
This paper describes the application of the Smoothed Particle Hydrodynamics（SPH） method for modeling two dimensional waves caused by dam break over a movable bed in two dimensions.The two phase SPH method is developed to solve the Navier-Stokes equations.Both fluid and sediment phases are described by particles as weakly compressible fluids and the incompressibility is achieved by the equation of state.The sediment phase is modeled as a non-Newtonian fluid using three alternative approaches of artificial viscosity and Bingham Model.In this paper,the new formulations for two-phase flows are proposed.The numerical results obtained from the developed SPH model show acceptable accuracy with comparison to experimental data. 相似文献
9.
A fundamental understanding of flow in porous media at the pore-scale is necessary to be able to upscale average displacement processes from core to reservoir scale. The study of fluid flow in porous media at the pore-scale consists of two key procedures: Imaging - reconstruction of three-dimensional (3D) pore space images; and modelling such as with single and two-phase flow simulations with Lattice-Boltzmann (LB) or Pore-Network (PN) Modelling. Here we analyse pore-scale results to predict petrophysical properties such as porosity, single-phase permeability and multi-phase properties at different length scales. The fundamental issue is to understand the image resolution dependency of transport properties, in order to up-scale the flow physics from pore to core scale. In this work, we use a high resolution micro-computed tomography (micro-CT) scanner to image and reconstruct three dimensional pore-scale images of five sandstones (Bentheimer, Berea, Clashach, Doddington and Stainton) and five complex carbonates (Ketton, Estaillades, Middle Eastern sample 3, Middle Eastern sample 5 and Indiana Limestone 1) at four different voxel resolutions (4.4 µm, 6.2 µm, 8.3 µm and 10.2 µm), scanning the same physical field of view. Implementing three phase segmentation (macro-pore phase, intermediate phase and grain phase) on pore-scale images helps to understand the importance of connected macro-porosity in the fluid flow for the samples studied. We then compute the petrophysical properties for all the samples using PN and LB simulations in order to study the influence of voxel resolution on petrophysical properties. We then introduce a numerical coarsening scheme which is used to coarsen a high voxel resolution image (4.4 µm) to lower resolutions (6.2 µm, 8.3 µm and 10.2 µm) and study the impact of coarsening data on macroscopic and multi-phase properties. Numerical coarsening of high resolution data is found to be superior to using a lower resolution scan because it avoids the problem of partial volume effects and reduces the scaling effect by preserving the pore-space properties influencing the transport properties. This is evidently compared in this study by predicting several pore network properties such as number of pores and throats, average pore and throat radius and coordination number for both scan based analysis and numerical coarsened data. 相似文献
10.
We have conducted a series of high-resolution numerical experiments using the Pair-Wise Force Smoothed Particle Hydrodynamics (PF-SPH) multiphase flow model. First, we derived analytical expressions relating parameters in the PF-SPH model to the surface tension and static contact angle. Next, we used the model to study viscous fingering, capillary fingering, and stable displacement of immiscible fluids in porous media for a wide range of capillary numbers and viscosity ratios. We demonstrated that the steady state saturation profiles and the boundaries of viscous fingering, capillary fingering, and stable displacement regions compare favorably with micromodel laboratory experimental results. For a displacing fluid with low viscosity, we observed that the displacement pattern changes from viscous fingering to stable displacement with increasing injection rate. When a high viscosity fluid is injected, transition behavior from capillary fingering to stable displacement occurred as the flow rate was increased. These observations are also in agreement with the results of the micromodel laboratory experiments. 相似文献
11.
阐述了大吨位双出杆粘滞阻尼器开发的数值与实验研究结果.基于硅油的幂律特性,通过对硅油介质的流态分析,介绍了非牛顿幂律流体模型对粘滞阻尼器耗能和阻尼出力的影响,与粘滞阻尼器的分析模型和计算公式.在商用流体软件平台上利用动网格技术对样品粘滞阻尼器阻尼出力进行了计算;综合实验研究结果,最终确定了大吨位土木用粘滞阻尼器设计的经验方法.另一方面,实验曲线和数值曲线基本吻和,验证了通过数值手段对粘滞阻尼器的出力特性进行预分析,从而提高开发精度、指导产品设计和缩短产品的开发周期的可行性. 相似文献
12.
Vittorio Di Federico Marco Pinelli Rita Ugarelli 《Stochastic Environmental Research and Risk Assessment (SERRA)》2010,24(7):1067-1076
An understanding of the interplay between non-Newtonian effects in porous media flow and field-scale domain heterogeneity
is of great importance in several engineering and geological applications. Here we present a simplified approach to the derivation
of an effective permeability for flow of a purely viscous power–law fluid with flow behavior index n in a randomly heterogeneous porous domain subject to a uniform pressure gradient. A standard form of the flow law generalizing
the Darcy’s law to non-Newtonian fluids is adopted, with the permeability coefficient being the only source of randomness.
The natural logarithm of the permeability is considered a spatially homogeneous and correlated Gaussian random field. Under
the ergodic hypothesis, an effective permeability is first derived for two limit 1-D flow geometries: flow parallel to permeability
variation (serial-type layers), and flow transverse to permeability variation (parallel-type layers). The effective permeability
of a 2-D or 3-D isotropic domain is conjectured to be a power average of 1-D results, generalizing results valid for Newtonian
fluids under the validity of Darcy’s law; the conjecture is validated comparing our results with previous literature findings.
The conjecture is then extended, allowing the exponents of the power averaging to be functions of the flow behavior index.
For Newtonian flow, novel expressions for the effective permeability reduce to those derived in the past. The effective permeability
is shown to be a function of flow dimensionality, domain heterogeneity, and flow behavior index. The impact of heterogeneity
is significant, especially for shear-thinning fluids with a low flow behavior index, which tend to exhibit channeling behavior. 相似文献
13.
基于BISQ模型的三维双相各向异性介质数值模拟 总被引：5，自引：2，他引：3
Biot-flow and squirt-flow are the two most important fluid flow mechanisms in porous media containing fluids. Based on the BISQ （Biot-Squirt） model where the two mechanisms are treated simultaneously, the elastic wave-field simulation in the porous medium is limited to two-dimensions and two-components （2D2C） or two-dimensions and three-components （2D3C）. There is no previous report on wave simulation in three- dimensions and three-components. Only through three dimensional numerical simulations can we have an overall understanding of wave field coupling relations and the spatial distribution characteristics between the solid and fluid phases in the dual-phase anisotropic medium. In this paper, based on the BISQ equation, we present elastic wave propagation in a three dimensional dual-phase anisotropic medium simulated by the staggered-grid high-order finite-difference method. We analyze the resulting wave fields and show that the results are an improvement. 相似文献
14.
Lateral deformation of liquefiable soil is a cause of much damage during earthquakes, reportedly more than other forms of liquefaction-induced ground failures. Researchers have presented studies in which the liquefied soil is considered as viscous fluid. In this manner, the liquefied soil behaves as non-Newtonian fluid, whose viscosity decreases as the shear strain rate increases. The current study incorporates computational fluid dynamics to propose a simplified dynamic analysis for the liquefaction-induced lateral deformation of earth slopes. The numerical procedure involves a quasi-linear elastic model for small to moderate strains and a Bingham fluid model for large strain states during liquefaction. An iterative procedure is considered to estimate the strain-compatible shear stiffness of soil. The post-liquefaction residual strength of soil is considered as the initial Bingham viscosity. Performance of the numerical procedure is examined by using the results of centrifuge model and shaking table tests together with some field observations of lateral ground deformation. The results demonstrate that the proposed procedure predicts the time history of lateral ground deformation with a reasonable degree of precision. 相似文献
15.
地下岩石由岩石骨架和孔隙流体组成,通常流体含黏性.地震波在地下介质中传播时受岩石骨架和黏性流体的影响会呈现出复杂的变化.本文将流、固体位移和应力连续作为边界条件,推导出含黏性流体孔隙介质分界面上反透射系数方程;通过建立上层为饱油、下层为饱盐水的砂岩孔隙介质模型,开展反透射系数特征研究,分别分析不同频率、不同黏滞系数条件下,含黏性流体孔隙介质分界面上反透射系数随入射角的变化.研究表明,孔隙介质分界面上和等效介质分界面上的反透射系数分别随入射角的变化趋势基本一致,说明方程推导和数值计算的正确性;快纵波反透射系数受频率、流体黏性的影响较小,而快横波反透射系数在一定入射角范围内受频率、流体黏性的影响比较大;由于黏性孔隙流体的作用,慢纵波和慢横波的反透射系数受入射角、频率及流体黏性的影响都很大. 相似文献
16.
《Advances in water resources》1999,23(1):83-95
We present a numerical model for two-phase porous media flow, where the phases are separated by a sharp interface. The model is based on a unified pressure equation, and an advection equation for tracking a pseudo-concentration function. The zero-level set of this function defines the interface between the fluids. The finite element method is used for spatial discretization, with local grid refinements in the vicinity of the interface. Examples on applications involving moving interface and steady-state seepage problems are investigated. 相似文献
17.
以Biot饱和多孔介质动力学理论为基础,利用Fourier Bessel级展开法,给出了饱和多孔介质半空间中半球形凹陷地形对入射平面SV波散射问题的解析解.利用这一解析解数值计算地表位移幅值,分析入射角、入射波频率对地表位移幅值的影响,并与现存的单相介质情况下三维半球形凹陷问题解析解进行对比得出如下结论：(1)本文两相介质模型与单相介质模型有很大区别;(2)入射角、入射频率对场地表面位移分量有重要影响;(3)随着入射角、入射频率的增大,地表位移分布变得更复杂,且放大效应比单相介质模拟结果更加显著. 相似文献
18.
Saturation of porous rocks with a mixture of two fluids has a substantial effect on seismic‐wave propagation. In particular, partial saturation causes significant attenuation and dispersion of the propagating waves due to the mechanism of wave‐induced fluid‐flow. Such flow arises when a passing wave induces different fluid pressures in regions of rock saturated by different fluids. Most models of attenuation and dispersion due to mesoscopic heterogeneities imply that fluid heterogeneities are distributed in a regular way. However, recent experimental studies show that mesoscopic heterogeneities have less idealized distributions and that the distribution itself affects attenuation and dispersion. Based on an approximation for the coherent wavefield in random porous media, we develop a model which assumes a continuous distribution of fluid heterogeneities. As this continuous random media approach assumes that there will be a distribution of different patch sizes, it is expected to be better suited to modelling experimental data. We also show how to relate the random functions to experimentally measurable parameters. 相似文献
19.
将时空守恒元/解元(CE/SE)方法推广到二维孔隙介质多相流问题的数值计算中,采用人工压缩法耦合速度和压力,同时结合杂交粒子水平集方法捕捉物质界面.提出一套完整的二维欧拉型孔隙介质非稳态多相不可压缩黏性流动计算方案.通过对溃坝和液滴在重力作用下的运动和变形问题的数值模拟,验证了方法的精度和有效性.在此基础上,提出了一个新的孔隙介质两相流物理模型——双层流体顶盖驱动方腔流. 相似文献
20.
基于前一篇文章中得到的关于三维层状孔隙介质中弹性波场的积分形式半解析解,本文通过离散波数法开展了数值模拟.将全空间均匀孔隙介质中单力点源和爆炸点源作用下弹性波场的解析解和我们的数值模拟结果进行对比,发现两者是完全一致的.而在一个两层半空间模型下的数值模拟,验证了固相位移Green函数的9组空间互易性情况.通过以上两种对比检验,验证了半解析解理论公式、数值模拟方法以及相应程序代码的正确性和可靠性.随后利用敏感度分析研究了不同的介质参数变化对爆炸点源在界面上会产生的反射波场的影响.通过垂直地震剖面模型的数值模拟,发现弹性波场能很好地反映孔隙介质物理性质的变化,同时也讨论了动力协调这一孔隙介质中的特殊现象.我们发展的基于半解析解的数值模拟方法可以为三维层状孔隙介质中弹性波传播特征的研究提供一种可供选择的有效工具和手段. 相似文献