首页 | 本学科首页   官方微博 | 高级检索  
相似文献
 共查询到20条相似文献,搜索用时 15 毫秒
1.
基于临界孔隙度模型的地震波传播   总被引:1,自引:1,他引:0       下载免费PDF全文
基于岩石物里学中临界孔隙度模型,建立一种简洁的均匀弹性流体饱和孔隙介质模型,进行地震波传播研究.首先定义了构建目标模型的基本力学模型:介绍了全孔隙度区间内基本力学模型和目标孔隙介质的含义,其中基本力学模型除了完全弹性固体模型S和完全弹性流体模型F还包括临界孔隙模型C.然后通过等效力学模型推出了目标力学模型介质本构关系的组分表达形式.文中分别通过直接求取弹性参数的表达形式和运用应力应变关系两种方法得到介质模型的本构关系,进而得到该模型波动方程的组分表达形式.最后对这种介质模型进行了地震波传播的数值模拟,结合模拟结果分析孔隙对地震波传播的影响.  相似文献   

2.
二维格子气自动机模拟孔隙介质的电传输特性   总被引:14,自引:1,他引:13       下载免费PDF全文
本文应用2-D格子气自动机模拟饱和油水两相的多孔介质的导电特性。在油水两相界面处,引入反射与透射系数来决定粒子的运动状态,通过调节反射与透射系数就可以改变油水两相的导电性差异。用模拟结果考察Archie公式的地层因素F=R0/Rw和电阻率增大系数I=Rt/R0,其中,R0为百分之百饱和水时的岩石电阻率,Rw为水的电阻率,Rt为不同流体饮和度时的岩石电阻率。结果表明F与孔隙度φ间,I与含水饱和度间都存在幂关系,并可以表示为F=aφ-m,I=bS-nw模拟结果同时证实:公式中的参数a、m的变化反映了孔隙微观结构的变化,参数b、n主要受孔隙度大小和油相分布状态的影响。  相似文献   

3.
A three-dimensional, reactive numerical flow model is developed that couples chemical reactions with density-dependent mass transport and fluid flow. The model includes equilibrium reactions for the aqueous species, kinetic reactions between the solid and aqueous phases, and full coupling of porosity and permeability changes that result from precipitation and dissolution reactions in porous media. A one-step, global implicit approach is used to solve the coupled flow, transport and reaction equations with a fully implicit upstream-weighted control volume discretization. The Newton–Raphson method is applied to the discretized non-linear equations and a block ILU-preconditioned CGSTAB method is used to solve the resulting Jacobian matrix equations. This approach permits the solution of the complete set of governing equations for both concentration and pressure simultaneously affected by chemical and physical processes. A series of chemical transport simulations are conducted to investigate coupled processes of reactive chemical transport and density-dependent flow and their subsequent impact on the development of preferential flow paths in porous media. The coupled effects of the processes driving flow and the chemical reactions occurring during solute transport is studied using a carbonate system in fully saturated porous media. Results demonstrate that instability development is sensitive to the initial perturbation caused by density differences between the solute plume and the ambient groundwater. If the initial perturbation is large, then it acts as a “trigger” in the flow system that causes instabilities to develop in a planar reaction front. When permeability changes occur due to dissolution reactions occurring in the porous media, a reactive feedback loop is created by calcite dissolution and the mixed convective transport of the system. Although the feedback loop does not have a significant impact on plume shape, complex concentration distributions develop as a result of the instabilities generated in the flow system.  相似文献   

4.
本文定义了各向异性黏弹性参数修正因子,并将其引入到黏弹性模型中以体现泥质含量对黏弹性机制的影响,同时将波传播过程中孔隙介质骨架黏弹性力学机制与两种孔隙流体流动力学机制(Biot流动和喷射流动机制)有机地统一起来处理,从而给出了描述含泥质低孔渗孔隙各向异性介质中波传播规律的黏弹性Biot/squirt (BISQ)模型.数值计算结果表明,入射波的方位角、各向异性渗透率以及泥质含量等对含流体复杂孔隙介质中波频散和衰减的影响具有显著的方位各向异性特征,在低频范围内(地震波勘探频率)黏弹性力学机制对波传播能量的衰减起主导作用.  相似文献   

5.
Biofilm growth changes many physical properties of porous media such as porosity, permeability and mass transport parameters. The growth depends on various environmental conditions, and in particular, on flow rates. Modeling the evolution of such properties is difficult both at the porescale where the phase morphology can be distinguished, as well as during upscaling to the corescale effective properties. Experimental data on biofilm growth is also limited because its collection can interfere with the growth, while imaging itself presents challenges.In this paper we combine insight from imaging, experiments, and numerical simulations and visualization. The experimental dataset is based on glass beads domain inoculated by biomass which is subjected to various flow conditions promoting the growth of biomass and the appearance of a biofilm phase. The domain is imaged and the imaging data is used directly by a computational model for flow and transport. The results of the computational flow model are upscaled to produce conductivities which compare well with the experimentally obtained hydraulic properties of the medium. The flow model is also coupled to a newly developed biomass–nutrient growth model, and the model reproduces morphologies qualitatively similar to those observed in the experiment.  相似文献   

6.
《Advances in water resources》2007,30(6-7):1593-1607
An experimental and numerical investigation was conducted to study the colonization dynamics of a bioluminescent bacterium, Pseudomonas fluorescens HK44, during growth in a porous medium under steady, variably saturated flow conditions. Experiments were conducted in a thin-slab light transmission chamber filled with uniform, translucent quartz sand. Steady, variably saturated flow conditions were established using drip emitters mounted on the top of the chamber, with glucose applied through a central dripper located directly above an inoculated region of the chamber. Periodic pulses of salicylate and a dye tracer were applied to induce bioluminescence of the bacterium to monitor colony expansion and to track changes in the hydraulic and transport properties of the sand. Changes in the apparent water saturation of the sand were quantified by monitoring light transmission through the chamber with a CCD camera. The colonized region expanded laterally by about 15 cm, and upward against the flow by 7–8 cm during the 6-day experiment while apparent saturations in the colonized region decreased by 7–9% and the capillary fringe dropped by ∼5 cm. The observed data were reproduced approximately using a numerical model that accounted for the processes of water flow, solute and bacterial transport, cell growth and accumulation, glucose and oxygen consumption, and gas diffusion and exchange. The results of this study illustrate some of the complexities associated with coupled flow, reactive transport, and biological processes in variably saturated porous media, such as localized desaturation, capillary fringe lowering effects, and upstream movement of bacterial colonization, that may not readily observable using other experimental techniques.  相似文献   

7.
Unlike previous theories with velocity and/or elastic modulus averaging, we use a three-phase porous rock physics model developed by Santos for analyzing the seismic response of two immiscible fluids in saturated porous media. Considering reservoir reference pressure and coupling drag of two fluids in pores, the effects of frequency, porosity, and gas saturation on the phase velocities of the P-and S-waves are discussed in detail under field conditions. The effects of porosity and gas saturation on Vp/Vs are also provided. The data for our numerical experiments are from a sample of deep volcanic rock from Daqing. The numerical results show that the frequency dispersion effect can be ignored for deep volcanic rocks with low porosity and low permeability. It is concluded that for deep volcanic rocks the effect of gas content in pores on Vp/Vs is negligible but the effect of porosity is significant when there is a certain amount of water contained in the pores. The accurate estimate of lithology and porosity in this case is relatively more important.  相似文献   

8.
Diffusive mass exchange into immobile water regions within heterogeneous porous aquifers influences the fate of solutes. The percentage of immobile water is often unidentified in natural aquifers though. Hence, the mathematical prediction of solute transport in such heterogeneous aquifers remains challenging. The objective of this study was to find a simple analytical model approach that allows quantifying properties of mobile and immobile water regions and the portion of immobile water in a porous system. Therefore, the Single Fissure Dispersion Model (SFDM), which takes into account diffusive mass exchange between mobile and immobile water zones, was applied to model transport in well‐defined saturated dual‐porosity column experiments. Direct and indirect model validation was performed by running experiments at different flow velocities and using conservative tracer with different molecular diffusion coefficients. In another column setup, immobile water regions were randomly distributed to test the model applicability and to determine the portion of immobile water. In all setups, the tracer concentration curves showed differences in normalized maximum peak concentration, tailing and mass recovery according to their diffusion coefficients. These findings were more pronounced at lower flow rates (larger flow times) indicating the dependency of diffusive mass exchange into immobile water regions on tracers' molecular diffusion coefficients. The SFDM simulated all data with high model efficiency. Successful model validation supported the physical meaning of fitted model parameters. This study showed that the SFDM, developed for fissured aquifers, is applicable in porous media and can be used to determine porosity and volume of regions with immobile water. Copyright © 2015 John Wiley & Sons, Ltd.  相似文献   

9.
Determination of aqueous phase diffusion coefficients of solutes through porous media is essential for understanding and modeling contaminant transport. Prediction of diffusion coefficients in both saturated and unsaturated zones requires knowledge of tortuosity and constrictivity factors. No methods are available for the direct measurement of these factors, which are empirical in their definition. In this paper, a new definition for the tortuosity factor is proposed, as the real to ideal interfacial area ratio. We define the tortuosity factor for saturated porous media (tau5) as the ratio S/S(o) (specific surface of real porous medium to that of an idealized capillary bundle). For unsaturated media, tortuosity factor (tau(a)) is defined as a(aw)/a(aw),o (ratio of the specific air-water interfacial area of real and the corresponding idealized porous medium). This tortuosity factor is suitably measured using sorptive tracers (e.g., nitrogen adsorption method) for saturated media and interfacial tracers for unsaturated media. A model based on this new definition of tortuosity factors, termed the interfacial area ratio (IAR) model, is presented for the prediction of diffusion coefficients as a function of the degree of water saturation. Diffusion coefficients and diffusive resistances measured in a number of saturated and unsaturated granular porous media, for solutes in dilute aqueous solutions, agree well with the predictions of the IAR model. A comparison of permeability of saturated sands estimated based on tau(s) and the same based on the Kozeny-Carman equation confirm the usefulness of the tau(s) parameter as a measure of tortuosity.  相似文献   

10.
基于Biot理论,考虑液相的黏弹性变形和固液相接触面上的相对扭转,提出了含黏滞流体VTI孔隙介质模型.从理论上推导出,在该模型中除存在快P波、慢P波、SV波、SH波以外,还将存在两种新横波-慢SV波和慢SH波.数值模拟分析了6种弹性波的相速度、衰减、液固相振幅比随孔隙度、频率的变化规律以及快P波、快SV波的衰减随流体性质、渗透率、入射角的变化规律.结果表明慢SV波和慢SH波主要在液相中传播,高频高孔隙度时,速度较高;大角度入射时,快P波衰减表现出明显的各向异性,而快SV波的衰减则基本不变;储层纵向和横向渗透率存在差异时,快SV波衰减大的方向渗透率高.  相似文献   

11.
A common way to simulate fluid flow in porous media is to use Lattice Boltzmann (LB) methods. Permeability predictions from such flow simulations are controlled by parameters whose settings must be calibrated in order to produce realistic modelling results. Herein we focus on the simplest and most commonly used implementation of the LB method: the single-relaxation-time BGK model. A key parameter in the BGK model is the relaxation time τ which controls flow velocity and has a substantial influence on the permeability calculation. Currently there is no rigorous scheme to calibrate its value for models of real media. We show that the standard method of calibration, by matching the flow profile of the analytic Hagen-Poiseuille pipe-flow model, results in a BGK-LB model that is unable to accurately predict permeability even in simple realistic porous media (herein, Fontainebleau sandstone). In order to reconcile the differences between predicted permeability and experimental data, we propose a method to calibrate τ using an enhanced Transitional Markov Chain Monte Carlo method, which is suitable for parallel computer architectures. We also propose a porosity-dependent τ calibration that provides an excellent fit to experimental data and which creates an empirical model that can be used to choose τ for new samples of known porosity. Our Bayesian framework thus provides robust predictions of permeability of realistic porous media, herein demonstrated on the BGK-LB model, and should therefore replace the standard pipe-flow based methods of calibration for more complex media. The calibration methodology can also be extended to more advanced LB methods.  相似文献   

12.
The semi-empirical Kozeny–Carman (KC) equation is the most famous permeability–porosity relation, which is widely used in the field of flow in porous media and is the starting point for many other permeability models. However, this relation has many limitations from its inception, and the KC constant is an empirical parameter which was proved to be not a constant. In this paper, we briefly reviewed the KC equation, its modifications and various models for the KC constant. We then derived an analytical expression for the permeability in homogeneous porous media based on the fractal characters of porous media and capillary model. The proposed model is expressed as a function of fractal dimensions, porosity and maximum pore size. The analytical KC constant with no empirical constant is obtained from the assumption of square geometrical model. Furthermore, a distinct linear scaling law between the dimensionless permeability and porosity is found. It is also shown that our analytical permeability is more closely related to the microstructures (fractal dimensions, porosity and maximum pore size), compared to those obtained from conventional methods and models.  相似文献   

13.
Critical path analysis (CPA) is a method for estimating macroscopic transport coefficients of heterogeneous materials that are highly disordered at the micro-scale. Developed originally to model conduction in semiconductors, numerous researchers have noted that CPA might also have relevance to flow and transport processes in porous media. However, the results of several numerical investigations of critical path analysis on pore network models raise questions about the applicability of CPA to porous media. Among other things, these studies found that (i) in well-connected 3D networks, CPA predictions were inaccurate and became worse when heterogeneity was increased; and (ii) CPA could not fully explain the transport properties of 2D networks. To better understand the applicability of CPA to porous media, we made numerical computations of permeability and electrical conductivity on 2D and 3D networks with differing pore-size distributions and geometries. A new CPA model for the relationship between the permeability and electrical conductivity was found to be in good agreement with numerical data, and to be a significant improvement over a classical CPA model. In sufficiently disordered 3D networks, the new CPA prediction was within ±20% of the true value, and was nearly optimal in terms of minimizing the squared prediction errors across differing network configurations. The agreement of CPA predictions with 2D network computations was similarly good, although 2D networks are in general not well-suited for evaluating CPA. Numerical transport coefficients derived for regular 3D networks of slit-shaped pores were found to be in better agreement with experimental data from rock samples than were coefficients derived for networks of cylindrical pores.  相似文献   

14.
基于非饱和多孔隙介质BISQ模型的储层参数反演   总被引:22,自引:10,他引:12       下载免费PDF全文
Biot流动和喷射流动是含流体多孔隙介质中流体流动的两种重要力学机制.文中基于同时包含这两种力学机制的非饱和多孔隙BISQ模型,利用小生境遗传算法实现了储层参数(孔隙度、渗透率、含流体饱和度等)的反演.结果表明,本文方法在储层参数的反演过程中目标函数收敛快、且具有较强的抗“噪声”干扰性能,当观测噪声ε≤5%时,储层参数的反演精度很高.最后,通过实测数据的反演应用验证了该方法处理储层参数反演问题的有效性和非饱和BISQ模型的准确性.  相似文献   

15.
Streaming potential is the result of coupling between a fluid flow and an electric current in porous rocks. The modified Helmholtz–Smoluchowski equation derived for capillary tubes is mostly used to determine the streaming potential coefficient of porous media. However, to the best of our knowledge, the fractal geometry theory is not yet applied to analyse the streaming potential in porous media. In this article, a fractal model for the streaming potential coefficient in porous media is developed based on the fractal theory of porous media and on the streaming potential in a capillary. The proposed model is expressed in terms of the zeta potential at the solid?liquid interface, the minimum and maximum pore/capillary radii, the fractal dimension, and the porosity of porous media. The model is also examined by using another capillary size distribution available in published articles. The results obtained from the model using two different capillary size distributions are in good agreement with each other. The model predictions are then compared with experimental data in the literature and those based on the modified Helmholtz–Smoluchowski equation. It is shown that the predictions from the proposed fractal model are in good agreement with experimental data. In addition, the proposed model is able to reproduce the same result as the Helmholtz–Smoluchowski equation, particularly for high fluid conductivity or large grain diameters. Other factors influencing the streaming potential coefficient in porous media are also analysed.  相似文献   

16.
 A stochastic simulation is performed to study multiphase flow and contaminant transport in fractal porous media with evolving scales of heterogeneity. Numerical simulations of residual NAPL mass transfer and subsequent transport of dissolved and/or volatilized NAPL mass in variably saturated media are carried out in conjunction with Monte Carlo techniques. The impact of fractal dimension, plume scale and anisotropy (stratification) of fractal media on relative dispersivities is investigated and discussed. The results indicate the significance of evolving scale of porous media heterogeneity to the NAPL transport in the subsurface. In general, the fractal porous media enhance the dispersivities of NAPL mass plume transport in both the water phase and the gas phase while the influence on the water phase is more significant. The porous media with larger fractal dimension have larger relative dispersivities. The aqueous horizontal dispersivity exhibits a most significant increase against the plume scale.  相似文献   

17.
本文综合考虑了在波传播过程中孔隙介质的三种重要力学机制——"Biot流动机制一squirt流动机制-固体骨架黏弹性机制",借鉴等效介质思想,将含水饱和度引入波动力学控制方程,并考虑了不同波频率下孔隙流体分布模式对其等效体积模量的影响,给出了能处理含粘滞性非饱和流体孔隙介质中波传播问题的黏弹性Biot/squirt(BISQ)模型。推导了时间-空间域的波动力学方程组,由一组平面谐波解假设,给出频率-波数域黏弹性BISQ模型的相速度和衰减系数表达式。基于数值算例分析了含水饱和度、渗透率与频率对纵波速度和衰减的影响,并结合致密砂岩和碳酸盐岩的实测数据,对非饱和情况下的储层纵波速度进行了外推,碳酸盐岩储层中纵波速度对含气饱和度的敏感性明显低于砂岩储层。  相似文献   

18.
19.
An equivalent medium model for wave simulation in fractured porous rocks   总被引:3,自引:0,他引:3  
Seismic wave propagation in reservoir rocks is often strongly affected by fractures and micropores. Elastic properties of fractured reservoirs are studied using a fractured porous rock model, in which fractures are considered to be embedded in a homogeneous porous background. The paper presents an equivalent media model for fractured porous rocks. Fractures are described in a stress‐strain relationship in terms of fracture‐induced anisotropy. The equations of poroelasticity are used to describe the background porous matrix and the contents of the fractures are inserted into a matrix. Based on the fractured equivalent‐medium theory and Biot's equations of poroelasticity, two sets of porosity are considered in a constitutive equation. The porous matrix permeability and fracture permeability are analysed by using the continuum media seepage theory in equations of motion. We then design a fractured porous equivalent medium and derive the modified effective constants for low‐frequency elastic constants due to the presence of fractures. The expressions of elastic constants are concise and are directly related to the properties of the main porous matrix, the inserted fractures and the pore fluid. The phase velocity and attenuation of the fractured porous equivalent media are investigated based on this model. Numerical simulations are performed. We show that the fractures and pores strongly influence wave propagation, induce anisotropy and cause poroelastic behaviour in the wavefields. We observe that the presence of fractures gives rise to changes in phase velocity and attenuation, especially for the slow P‐wave in the direction parallel to the fracture plane.  相似文献   

20.
针对饱和多孔介质中热弹性波的传播特性问题,基于多孔介质理论和广义的热弹性模型,研究平面S波在饱和多孔热弹性介质边界上的反射问题.以考虑流-固耦合的饱和多孔介质波动方程和热-弹耦合的广义热弹性基本方程出发,建立饱和多孔介质的热-流-固耦合弹性波动模型.通过引入势函数并考虑自由透水和绝热的边界条件,经过理论推导最终给出在饱...  相似文献   

设为首页 | 免责声明 | 关于勤云 | 加入收藏

Copyright©北京勤云科技发展有限公司  京ICP备09084417号