共查询到20条相似文献，搜索用时 593 毫秒
1.
We consider the impact of using time-lapse seismic data in addition to production data for permeability estimation in a porous
medium with multiphase fluid flows, such as a petroleum reservoir under water-assisted production. Since modeling seismic
wave propagation in addition to modeling fluid flows in the reservoir is quite involved, it is assumed that the time-lapse
seismic data have already been inverted into fluid saturation differences (pseudoseismic data). Because an inversion process
often leads to considerable error growth, we will consider pseudoseismic data with large uncertainties. The impact of pseudoseismic
data is assessed through permeability estimation with and without such data and through application of some uncertainty measures
for the estimated parameters. A multiscale algorithm is used for the parameter estimations, so that potential differences
in attainable permeability resolution will be easily revealed. The numerical examples clearly indicate that the permeability
estimation problem is stabilized at a higher level of resolution when pseudoseismic data are applied in addition to production
data, even if the pseudoseismic data have large associated uncertainties. Use of the parameter uncertainty measures confirm
these results. 相似文献
2.
We describe a new approach for simulation of multiphase flows through heterogeneous porous media, such as oil reservoirs.
The method, which is based on the wavelet transformation of the spatial distribution of the single-phase permeabilities, incorporates
in the upscaled computational grid all the relevant data on the permeability, porosity, and other important properties of
a porous medium at all the length scales. The upscaling method generates a nonuniform computational grid which preserves the resolved structure
of the geological model in the near-well zones as well as in the high-permeability sectors and upscales the rest of the geological
model. As such, the method is a multiscale one that preserves all the important information across all the relevant length
scales. Using a robust front-detection method which eliminates the numerical dispersion by a high-order total variation diminishing
method (suitable for the type of nonuniform upscaled grid that we generate), we obtain highly accurate results with a greatly
reduced computational cost. The speed-up in the computations is up to over three orders of magnitude, depending on the degree
of heterogeneity of the model. To demonstrate the accuracy and efficiency of our methods, five distinct models (including
one with fractures) of heterogeneous porous media are considered, and two-phase flows in the models are studied, with and
without the capillary pressure. 相似文献
3.
We consider the inverse problem of permeability estimation for two-phase porous-media flow. The novel approach is based on
regularization by zonation, where the geometry and size of the regions are chosen adaptively during the optimization procedure.
To achieve this, we have utilized level-set functions to represent the permeability. The available data are sparsely distributed
in space; hence, it is reasonable to confine the estimation to coarse-scale structures. The level-set approach is able to
alter the boundaries between regions of different permeability without strict restrictions on their shape; however, when the
data are sparse, a reasonable initial guess for the permeability is required. For this task, we use adaptive multiscale permeability
estimation, which has the potential of identifying main permeability variations. These are described by a piecewise constant
function, where the constant values are attained on rectangular zones. In the current work, we develop a level-set corrector
strategy, assuming adaptive multiscale permeability estimation as a predictor. 相似文献
4.
In this paper, we introduce a novel stochastic model for the permeability tensor associated with stationary random porous media. In the light of recent works on mesoscale modeling of permeability, we first discuss the physical interpretation of the permeability tensor randomness. Subsequently, we propose a nonparametric prior probabilistic model for non‐Gaussian permeability tensor random fields, making use of the information theory and a maximum entropy procedure, and provide a physical interpretation of the model parameters. Finally, we demonstrate the capability of the considered class of random fields to generate higher levels of statistical fluctuations for selected stochastic principal permeabilities. This unique flexibility offered by the parameterization of the model opens up many new possibilities for both forward simulations (e.g. for uncertainty propagation in predictive simulations) and stochastic inverse problem solving. Copyright © 2011 John Wiley & Sons, Ltd. 相似文献
5.
The problem of multiphase phase flow in heterogeneous subsurface porous media is one involving many uncertainties. In particular, the permeability of the medium is an important aspect of the model that is inherently uncertain. Properly quantifying these uncertainties is essential in order to make reliable probabilistic-based predictions and future decisions. In this work, a measure-theoretic framework is employed to quantify uncertainties in a two-phase subsurface flow model in high-contrast media. Given uncertain saturation data from observation wells, the stochastic inverse problem is solved numerically in order to obtain a probability measure on the space of unknown permeability parameters characterizing the two-phase flow. As solving the stochastic inverse problem requires a number of forward model solves, we also incorporate the use of a conservative version of the generalized multiscale finite element method for added efficiency. The parameter-space probability measure is used in order to make predictions of saturation values where measurements are not available, and to validate the effectiveness of the proposed approach in the context of fine and coarse model solves. A number of numerical examples are offered to illustrate the measure-theoretic methodology for solving the stochastic inverse problem using both fine and coarse solution schemes. 相似文献
6.
In this paper, we propose multilevel Monte Carlo (MLMC) methods that use ensemble level mixed multiscale methods in the simulations of multiphase flow and transport. The contribution of this paper is twofold: (1) a design of ensemble level mixed multiscale finite element methods and (2) a novel use of mixed multiscale finite element methods within multilevel Monte Carlo techniques to speed up the computations. The main idea of ensemble level multiscale methods is to construct local multiscale basis functions that can be used for any member of the ensemble. In this paper, we consider two ensemble level mixed multiscale finite element methods: (1) the no-local-solve-online ensemble level method (NLSO); and (2) the local-solve-online ensemble level method (LSO). The first approach was proposed in Aarnes and Efendiev (SIAM J. Sci. Comput. 30(5):2319-2339, 2008) while the second approach is new. Both mixed multiscale methods use a number of snapshots of the permeability media in generating multiscale basis functions. As a result, in the off-line stage, we construct multiple basis functions for each coarse region where basis functions correspond to different realizations. In the no-local-solve-online ensemble level method, one uses the whole set of precomputed basis functions to approximate the solution for an arbitrary realization. In the local-solve-online ensemble level method, one uses the precomputed functions to construct a multiscale basis for a particular realization. With this basis, the solution corresponding to this particular realization is approximated in LSO mixed multiscale finite element method (MsFEM). In both approaches, the accuracy of the method is related to the number of snapshots computed based on different realizations that one uses to precompute a multiscale basis. In this paper, ensemble level multiscale methods are used in multilevel Monte Carlo methods (Giles 2008a, Oper.Res. 56(3):607-617, b). In multilevel Monte Carlo methods, more accurate (and expensive) forward simulations are run with fewer samples, while less accurate (and inexpensive) forward simulations are run with a larger number of samples. Selecting the number of expensive and inexpensive simulations based on the number of coarse degrees of freedom, one can show that MLMC methods can provide better accuracy at the same cost as Monte Carlo (MC) methods. The main objective of the paper is twofold. First, we would like to compare NLSO and LSO mixed MsFEMs. Further, we use both approaches in the context of MLMC to speedup MC calculations. 相似文献
7.
Micro-CT and FIB–SEM imaging and pore structure characterization of dolomite rock at multiple scales
Rajakumar S. Devarapalli Amina Islam Titly Farhana Faisal Mohamed Sassi Mustapha Jouiad 《Arabian Journal of Geosciences》2017,10(16):361
With improvements of imaging techniques and computational power, Digital Rock Physics (DRP) has been increasingly used to determine transport and elastic properties of reservoir core plugs. Since numerical computations highly rely on accurate 3D representations of the porous microstructure of the rocks, the imaging technique and the scale at which the imaging is performed is a critical parameter. In this paper, we introduce a multiscale imaging workflow that uses both micro-X-ray tomography (micro-XCT) and focused ion beam combined with scanning electron microscope (FIB–SEM) to characterize a dolomite rock from the microscale to the nanoscale. This allows for the accurate capture of the different heterogeneities that exist in the carbonate (texture, mineralogy, pore size). The reconstructed microporous structures were then used to successfully predict elastic and permeability properties of selected carbonate. 相似文献
8.
黄土作为一种典型的多孔介质,多孔性是其重要特性之一。多孔性不仅影响着黄土的物理力学及化学特征,也严重影响着黄土的水理性质,特别是黄土的渗透特性。然而,作为多孔介质渗流理论的经典模型——毛细管模型,已经广泛应用于油气田开采、煤成气开采以及岩土工程等领域,但将多孔介质渗流理论模型引入表征黄土的多孔构造,这方面的资料尚欠缺。因此,本文以黄土的渗透性作为研究示例,在总结分析目前相对较为成功的多孔介质毛细管渗流模型的基础上,选取适于定量表征黄土渗透性的渗流模型,结合计算结果和渗透试验结果的对比,探讨将多孔介质毛细管渗流模型引入描述黄土渗透性这一方法的适用性。并提出孔隙的“香肠”构造（或称为“莲藕”构造）,为黄土的渗透性研究提供可借鉴的理论依据。 相似文献
9.
Victor Ginting Felipe Pereira Michael Presho Shaochang Wo 《Computational Geosciences》2011,15(4):691-707
In this paper, we develop a procedure for subsurface characterization of a fractured porous medium. The characterization involves
sampling from a representation of a fracture’s permeability that has been suitably adjusted to the dynamic tracer cut measurement
data. We propose to use a type of dual-porosity, dual-permeability model for tracer flow. This model is built into the Markov
chain Monte Carlo (MCMC) method in which the permeability is sampled. The Bayesian statistical framework is used to set the
acceptance criteria of these samples and is enforced through sampling from the posterior distribution of the permeability
fields conditioned to dynamic tracer cut data. In order to get a sample from the distribution, we must solve a series of problems
which requires a fine-scale solution of the dual model. As direct MCMC is a costly method with the possibility of a low acceptance
rate, we introduce a two-stage MCMC alternative which requires a suitable coarse-scale solution method of the dual model.
With this filtering process, we are able to decrease our computational time as well as increase the proposal acceptance rate.
A number of numerical examples are presented to illustrate the performance of the method. 相似文献
10.
Behnam V. Damirchi Marcelo R. Carvalho Luís A. G. Bitencourt Osvaldo L. Manzoli Daniel Dias‐da‐Costa 《国际地质力学数值与分析法杂志》2021,45(1):83-107
A new discrete fracture model is introduced to simulate the steady‐state fluid flow in discontinuous porous media. The formulation uses a multi‐layered approach to capture the effect of both longitudinal and transverse permeability of the discontinuities in the pressure distribution. The formulation allows the independent discretisation of mesh and discontinuities, which do not need to conform. Given that the formulation is developed at the element level, no additional degrees of freedom or special integration procedures are required for coupling the non‐conforming meshes. The proposed model is shown to be reliable regardless of the permeability of the discontinuity being higher or lower than the surrounding domain. Four numerical examples of increasing complexity are solved to demonstrate the efficiency and accuracy of the new technique when compared with results available in the literature. Results show that the proposed method can simulate the fluid pressure distribution in fractured porous media. Furthermore, a sensitivity analysis demonstrated the stability regarding the condition number for wide range values of the coupling parameter. 相似文献
11.
The prediction of fluid flows within hydrocarbon reservoirs requires the characterization of petrophysical properties. Such characterization is performed on the basis of geostatistics and history-matching; in short, a reservoir model is first randomly drawn, and then sequentially adjusted until it reproduces the available dynamic data. Two main concerns typical of the problem under consideration are the heterogeneity of rocks occurring at all scales and the use of data of distinct resolution levels. Therefore, referring to sequential Gaussian simulation, this paper proposes a new stochastic simulation method able to handle several scales for both continuous or discrete random fields. This method adds flexibility to history-matching as it boils down to the multiscale parameterization of reservoir models. In other words, reservoir models can be updated at either coarse or fine scales, or both. Parameterization adapts to the available data; the coarser the scale targeted, the smaller the number of unknown parameters, and the more efficient the history-matching process. This paper focuses on the use of variational optimization techniques driven by the gradual deformation method to vary reservoir models. Other data assimilation methods and perturbation processes could have been envisioned as well. Last, a numerical application case is presented in order to highlight the advantages of the proposed method for conditioning permeability models to dynamic data. For simplicity, we focus on two-scale processes. The coarse scale describes the variations in the trend while the fine scale characterizes local variations around the trend. The relationships between data resolution and parameterization are investigated. 相似文献
12.
This paper is aimed towards investigating the filtration law of an incompressible viscous Newtonian fluid through a rigid non-inertial porous medium (e.g. a porous medium placed in a centrifuge basket). The filtration law is obtained by upscaling the flow equations at the pore scale. The upscaling technique is the homogenization method of multiple scale expansions which rigorously gives the macroscopic behaviour and the effective properties without any prerequisite on the form of the macroscopic equations. The derived filtration law is similar to Darcy's law, but the tensor of permeability presents the following remarkable properties: it depends upon the angular velocity of the porous matrix, it verifies Hall–Onsager's relationship and it is a non-symmetric tensor. We thus deduce that, under rotation, an isotropic porous medium leads to a non-isotropic effective permeability. In this paper, we present the results of numerical simulations of the flow through rotating porous media. This allows us to highlight the deviations of the flow due to Coriolis effects at both the microscopic scale (i.e. the pore scale), and the macroscopic scale (i.e. the sample scale). The above results confirm that for an isotropic medium, phenomenological laws already proposed in the literature fails at reproducing three-dimensional Coriolis effects in all types of pores geometry. We show that Coriolis effects may lead to significant variations of the permeability measured during centrifuge tests when the inverse Ekman number Ek−1 is 𝒪(1). These variations are estimated to be less than 5% if Ek−1<0.2, which is the case of classical geotechnical centrifuge tests. We finally conclude by showing that available experimental data from tests carried out in centrifuges are not sufficient to determining the effective tensor of permeability of rotating porous media. Copyright © 2004 John Wiley & Sons, Ltd. 相似文献
13.
Multiscale estimation of the Freundlich adsorption isotherm 总被引：1，自引：1，他引：0
M. N. Nounou Ph.D. H. N. Nounou Ph.D. 《International Journal of Environmental Science and Technology》2010,7(3):509-518
Adsorption plays an important role in water and wastewater treatment. The analysis and design of processes that involve adsorption rely on the availability of isotherms that describe these adsorption processes. Adsorption isotherms are usually estimated empirically from measurements of the adsorption process variables. Unfortunately, these measurements are usually contaminated with errors that degrade the accuracy of estimated isotherms. Therefore, these errors need to be filtered for improved isotherm estimation accuracy. Multiscale wavelet based filtering has been shown to be a powerful filtering tool. In this work, multiscale filtering is utilized to improve the estimation accuracy of the Freundlich adsorption isotherm in the presence of measurement noise in the data by developing a multiscale algorithm for the estimation of Freundlich isotherm parameters. The idea behind the algorithm is to use multiscale filtering to filter the data at different scales, use the filtered data from all scales to construct multiple isotherms and then select among all scales the isotherm that best represents the data based on a cross validation mean squares error criterion. The developed multiscale isotherm estimation algorithm is shown to outperform the conventional time-domain estimation method through a simulated example. 相似文献
14.
The empirical Darcy's law of water transport in porous media, Fick's law of chemical diffusion, and Fourier's law of thermal transport have been widely used in geophysics/geochemistry for over 150 years. However, the strong couplings between water, temperature, and chemicals in a membrane porous medium have made these laws inapplicable and present a significant hurdle to the understanding of multiphase flow in such a material. Extensive experiments over the past century have observed chemical osmosis and thermal osmosis, but a model for understanding their underlying physicochemical basis has remained unavailable, because of the highly cross‐disciplinary and multiscale‐multiphase nature of the coupling. Based on the fundamental principles of nonequilibrium thermodynamics and mixture coupling theory, a rigorously theoretical and mathematical framework is proposed and a general model accounting for all of the coupled influences is developed. This leads to a simple and robust mathematical matrix for studying multiphase couplings in a membrane porous medium when all chemical components are electrically neutral. 相似文献
15.
裂隙岩体的渗透特性受控于裂隙的发育特征、连通特性和充填情况,并与岩体的地应力水平具有显著的相关性。基于裂隙岩体渗透性的影响因素,并考虑现有渗透系数估算模型的不足,利用钻孔压水试验和钻孔电视图像资料,建立考虑埋深(Z)、岩石质量指标(RQD)以及充填物指标(FSD)等3个指标的渗透系数估算ZRF模型,并应用到牙根二级水电站及其他工程区的渗透系数估算中。结果表明,与已有的渗透系数估算模型相比,ZRF模型较好地反映了岩体渗透性的影响因素,且模型参数物理意义明确,便于获取,对分析裂隙岩体渗透性具有一定的工程参考价值。 相似文献
16.
测井多尺度分析方法用于层序地层划分研究 总被引：1，自引：0，他引：1
测井数据包含了丰富的地质信息,是研究地层多尺度沉积旋回的主要资料.本文阐述了小波变换及多尺度分析方法,探讨了测井多尺度分析方法在层序地层划分中的应用.以东营凹陷某井为例,选取Morlet小波基函数对GR测井曲线进行连续小波变换,将测井信号与深度的关系转换为与深度和尺度域的变化关系.通过研究多种伸缩尺度下小波系数曲线表现出的周期性振荡特征,并结合不同测井曲线多尺度分解后的高频信号特征,划分出各级层序界面,与传统方法所划分的界面基本一致. 相似文献
17.
非常规油气资源的孔隙结构及其连通性非常复杂,其孔隙尺度从毫米到纳米跨越多个量级.多孔介质中气体的输运过程不仅依赖于介质的多尺度微观结构特征,还依赖于气体的相关属性.气体在多尺度多孔介质中的输运过程包括无滑流、滑脱流和过渡流,涉及分子扩散和努森扩散等多种机制,因此很难用唯一的连续介质理论来描述气体的输运特征.大量的数据表明真实多孔介质中的内部孔隙具有分形标度特征,因此采用分形几何表征多尺度多孔介质的孔隙结构,引入孔隙分形维数和迂曲度分形维数定量表征多孔介质的微结构和弯曲流道,建立多尺度多孔介质气体输运过程的细观模型;推导了多尺度多孔介质中气体的有效渗透率和有效扩散系数,并讨论了多尺度多孔介质微结构参数和气体属性对于气体等效输运特性的定量影响.该研究不仅可以丰富渗流理论,且有利于深入理解非常规油气藏的产出机制. 相似文献
18.
In the analysis of petroleum reservoirs, one of the most challenging problems is to use inverse theory in the search for an optimal parameterization of the reservoir. Generally, scientists approach this problem by computing a sensitivity matrix and then perform a singular value decomposition in order to determine the number of degrees of freedom i.e. the number of independent parameters necessary to specify the configuration of the system. Here we propose a complementary approach: it uses the concept of refinement indicators to select those degrees which have the greatest sensitivity to an objective function quantifying the mismatch between measured and simulated data. We apply this approach to the problem of data integration for petrophysical reservoir charaterization where geoscientists are currently working with multimillion cell geological models. Data integration may be performed by gradually deforming (by a linear combination) a set of these multimillion grid geostatistical realizations during the optimization process. The inversion parameters are then reduced to the number of coefficients of this linear combination. However, there is an infinity of geostatistical realizations to choose from which may not be efficient regarding operational constraints. Following our new approach, we are able through a single objective function evaluation to compute refinement indicators that indicate which realizations might improve the iterative geological model in a significant way. This computation is extremely fast as it implies a single gradient computation through the adjoint state approach and dot products. Using only the most sensitive realizations from a given set, we are able to resolve quicker the optimization problem case. We applied this methodology to the integration of interference test data into 3D geostatistical models. 相似文献
19.
We review and perform comparison studies for three recent multiscale methods for solving elliptic problems in porous media
flow; the multiscale mixed finite-element method, the numerical subgrid upscaling method, and the multiscale finite-volume
method. These methods are based on a hierarchical strategy, where the global flow equations are solved on a coarsened mesh
only. However, for each method, the discrete formulation of the partial differential equations on the coarse mesh is designed
in a particular fashion to account for the impact of heterogeneous subgrid structures of the porous medium. The three multiscale
methods produce solutions that are mass conservative on the underlying fine mesh. The methods may therefore be viewed as efficient,
approximate fine-scale solvers, i.e., as an inexpensive alternative to solving the elliptic problem on the fine mesh. In addition,
the methods may be utilized as an alternative to upscaling, as they generate mass-conservative solutions on the coarse mesh.
We therefore choose to also compare the multiscale methods with a state-of-the-art upscaling method – the adaptive local–global
upscaling method, which may be viewed as a multiscale method when coupled with a mass-conservative downscaling procedure.
We investigate the properties of all four methods through a series of numerical experiments designed to reveal differences
with regard to accuracy and robustness. The numerical experiments reveal particular problems with some of the methods, and
these will be discussed in detail along with possible solutions. Next, we comment on implementational aspects and perform
a simple analysis and comparison of the computational costs associated with each of the methods. Finally, we apply the three
multiscale methods to a dynamic two-phase flow case and demonstrate that high efficiency and accurate results can be obtained
when the subgrid computations are made part of a preprocessing step and not updated, or updated infrequently, throughout the
simulation.
The research is funded by the Research Council of Norway under grant nos. 152732 and 158908. 相似文献