共查询到20条相似文献，搜索用时 78 毫秒
1.
J.D. Jansen 《Computational Geosciences》2002,6(2):195-206
This paper concerns the computation of near-well flow in numerical reservoir simulation with unstructured grids. In particular, it uses spherical trigonometry to derive analytical expressions for the flow towards a well modeled as either a number of point sources or a constant-flux line source. The expression for the point source representation is based on projections of the grid block boundaries on spheres with unit radius around the sources. The expression for the line source is based on projection on a prolate spheroid. The computation of the surface area is done through transformation to prolate spheroidal coordinates and subsequent projection on a sphere at infinity. The point source expression for a single source is exact for grid block boundaries with straight edges; the line source expression is an approximation. Both representations are fully volume conserving, such that the sum of the fluxes through the grid block boundaries surrounding a source adds up exactly to the total source flow rate. Both representations can be used to accurately model complicated wells in the form of segments. The point source representation is simpler to implement and not necessarily less accurate than the line source representation. 相似文献
2.
This paper presents a finite-volume method for hexahedral multiblock grids to calculate multiphase flow in geologically complex reservoirs. Accommodating complex geologic and geometric features in a reservoir model (e.g., faults) entails non-orthogonal and/or unstructured grids in place of conventional (globally structured) Cartesian grids. To obtain flexibility in gridding as well as efficient flow computation, we use hexahedral multiblock grids. These grids are locally structured, but globally unstructured. One major advantage of these grids over fully unstructured tetrahedral grids is that most numerical methods developed for structured grids can be directly used for dealing with the local problems. We present several challenging examples, generated via a commercially available tool, that demonstrate the capabilities of hexahedral multiblock gridding. Grid quality is discussed in terms of uniformity and orthogonality. The presence of non-orthogonal grid and full permeability tensors requires the use of multi-point discretization methods. A flux-continuous finite-difference (FCFD) scheme, previously developed for stratigraphic hexahedral grid with full-tensor permeability, is employed for numerical flow computation. We extend the FCFD scheme to handle exceptional configurations (i.e. three- or five-cell connections as opposed to the regular four), which result from employing multiblock gridding of certain complex objects. In order to perform flow simulation efficiently, we employ a two-level preconditioner for solving the linear equations that results from the wide stencil of the FCFD scheme. The individual block, composed of cells that form a structured grid, serves as the local level; the higher level operates on the global block configuration (i.e. unstructured component). The implementation uses an efficient data structure where each block is wrapped with a layer of neighboring cells. We also examine splitting techniques [14] for the linear systems associated with the wide stencils of our FCFD operator. We present three numerical examples that demonstrate the method: (1) a pinchout, (2) a faulted reservoir model with internal surfaces and (3) a real reservoir model with multiple faults and internal surfaces. 相似文献
3.
Michael G. Edwards 《Computational Geosciences》2002,6(3-4):433-452
Locally conservative flux-continuous, full-tensor, discretization schemes are presented for general unstructured grids. The schemes are control-volume distributed, where flow variables and rock properties are assigned to the polygonal control-volumes derived from the primal grid. A relationship between these finite volume schemes and the mixed finite element method is established. An extension for unstructured grids is described that leads to a general symmetric positive definite discretization matrix for both quadrilateral and triangular grids. A novel flow based gridding approach for unstructured mesh generation is also proposed for heterogeneous reservoir domains. Results computed with the flux continuous schemes on unstructured flow-based grids demonstrate the advantages of the methods. 相似文献
4.
?ystein Pettersen 《Computational Geosciences》2012,16(2):211-230
The layering in reservoir simulation grids is often based on the geology, e.g., structure tops. In this paper we investigate
the alternative of using horizontal layers, where the link to the geology model is by the representation of the petrophysics
alone. The obvious drawback is the failure to honor the structure in the grid geometry. On the other hand, a horizontal grid
will honor the initial fluid contacts perfectly, and horizontal wells can also be accurately represented. Both these issues
are vital in thin oil-zone problems, where horizontal grids may hence be a viable alternative. To investigate this question,
a number of equivalent simulation models were built for a segment of the Troll Field, both geology-based and horizontal, and
various combinations of these. In the paper, it is demonstrated that the horizontal grid was able to capture the essentials
of fluid flow with the same degree of accuracy as the geology-based grid, and near-well flow was considerably more accurate.
For grids of comparable resolution, more reliable results were obtained by a horizontal grid than a geo-grid. A geo-grid with
local grid refinement and a horizontal grid produced almost identical results, but the ratio of computing times was almost
20 in favor of the horizontal grid. In the one-phase regions of the reservoir, relatively coarse cells can be used without
significant loss of accuracy. 相似文献
5.
We simulate direct current (DC) borehole resistivity measurements acquired in steel-cased deviated wells for the assessment
of rock formation properties. The assumed data acquisition configuration considers one current (emitter) and three voltage
(collector) electrodes that are utilized to measure the second difference of the electric potential along the well trajectory.
We assume a homogeneous, 1.27-cm-thick steel casing with resistivity equal to 10 − 5 Ω· m. Simulations are performed with two different numerical methodologies. The first one is based on transferring two-dimensional
(2D) axisymmetric optimal grids to a three-dimensional (3D) simulation software. The second one automatically produces optimal
3D grids yielded by a 3D self-adaptive goal-oriented algorithm. Both methodologies utilize high-order finite elements (FE)
that are specially well-suited for problems with high-contrast coefficients and rapid spatial variations of the electric field,
as it occurs in simulations that involve steel-cased wells. The method based on transferring 2D-optimal grids is efficient
in terms of CPU time (few seconds per logging position). Unfortunately, it may produce inaccurate 3D simulations in deviated
wells, even though the error remains below 1% for the axisymmetric (vertical) well. The method based on optimal 3D grids,
although less efficient in terms of CPU time (few hours per logging position), produces more accurate results that are validated
by a built-in a posteriori error estimator. This paper provides the first existing simulations of through-casing resistivity
measurements in deviated wells. Simulated resistivity measurements indicate that, for a 30° deviated well, measurements in
conductive layers 0.01 Ω· m) are similar to those obtained in vertical wells. However, in resistive layers (10,000 Ω· m),
we observe 100% larger readings in the 30° deviated well. This difference becomes 3,000% for the case of a 60° deviated well.
For this highly-deviated well, readings corresponding to the conductive formation layer are about 30% smaller in magnitude
than those in a vertical well. Shoulder effects significantly vary in deviated wells. 相似文献
6.
The problem of calculating equivalent grid block permeability tensors for heterogeneous porous media is addressed. The homogenization method used involves solving Darcy's equation subject to linear boundary conditions with flux conservation in subregions of the reservoir and can be readily applied to unstructured grids. The resulting equivalent permeability tensor is stable as defined relative to G-convergence. It is proposed to use both conforming and mixed finite elements to solve the local problems and compute approximations from above and below of the equivalent permeability, respectively. Comparisons with results obtained using periodic, pressure and no-flux boundary conditions and the renormalization method are presented. A series of numerical examples demonstrates the effectiveness of the methodology for two-phase flow in heterogeneous reservoirs. 相似文献
7.
The problem of calculating equivalent grid block permeability tensors for heterogeneous porous media is addressed. The homogenization method used involves solving Darcy's equation subject to linear boundary conditions with flux conservation in subregions of the reservoir and can be readily applied to unstructured grids. The resulting equivalent permeability tensor is stable as defined relative to G-convergence. It is proposed to use both conforming and mixed finite elements to solve the local problems and compute approximations from above and below of the equivalent permeability, respectively. Comparisons with results obtained using periodic, pressure and no-flux boundary conditions and the renormalization method are presented. A series of numerical examples demonstrates the effectiveness of the methodology for two-phase flow in heterogeneous reservoirs. 相似文献
8.
为了改进计算区域离散化问题,本文利用自适应非结构化网格有限单元法求解二维地电结构下大地电磁场满足的加权余量表达式。在有限元求解电磁场的过程中,网格剖分越精细、计算精度越高,计算量也会越大。此外,结构化网格难以适应任意地形以及复杂地质构造。而自适应非结构化网格在电性变化剧烈的区域会自动加密,在电性缓变的区域则生成粗疏的网格,从而优化网格质量与数量。因此,文中引入COMSOL Multiphysics软件,以实现若干地电模型的构建及非结构化自由四边形单元网格化。将网格数据信息导入本文算法,计算大地电磁场响应,并与解析解及数值解对比。结果表明,基于非结构化网格的正演模拟精度高、适应性强,为计算区域网格化提供了新的方法。 相似文献
9.
The increasing use of unstructured grids for reservoir modeling motivates the development of geostatistical techniques to
populate them with properties such as facies proportions, porosity and permeability. Unstructured grids are often populated
by upscaling high-resolution regular grid models, but the size of the regular grid becomes unreasonably large to ensure that
there is sufficient resolution for small unstructured grid elements. The properties could be modeled directly on the unstructured
grid, which leads to an irregular configuration of points in the three-dimensional reservoir volume. Current implementations
of Gaussian simulation for geostatistics are for regular grids. This paper addresses important implementation details involved
in adapting sequential Gaussian simulation to populate irregular point configurations including general storage and computation
issues, generating random paths for improved long range variogram reproduction, and search strategies including the superblock
search and the k-dimensional tree. An efficient algorithm for computing the variogram of very large irregular point sets is developed for
model checking. 相似文献
10.
11.
直井开发煤层气钻井和压裂成本高,控制面积小,单井产气量低;煤层内水平井钻进难度大,风险高,薄煤层中井眼轨迹控制难度大,钻井液有害固相对储层伤害严重,采收率低。基于此,分析贵州织金区块煤系地质构造,提出在煤系地层内稳定的非储层内布水平井,通过压裂造缝沟通水平井上下煤层同时开发多层煤层的新思路。与常规开发方式相比,非储层内水平井具有钻井风险小、储层伤害小、单井产量高的优点,同时还可以开发煤系致密气和页岩气,提高非常规天然气利用率。研究非储层内水平井开发贵州织金煤层气技术,为解决贵州煤系地层煤层多而薄、层间距小等特性煤层气开发难题以及综合利用煤系气提供新的方式。 相似文献
12.
水平井的开发技术关键之一是建立油气藏的精细三雏地质模型,据此可以确定明确的地质目标.常规的储层地质建模忽视了水平井与地层在空间上的独特配置关系,无法体现和充分利用其在表征储层属性空间展布特征上的优势.这里着重讲述了针对水平井开发的储层地质建模,引入了水平井储层概念模型,探讨了平面单元网格划分模式,并以塔河1区三叠系下油... 相似文献
13.
赵晓博 《物探化探计算技术》2011,33(5):517-521,463
利用Delaunay三角化这种网格非结构化方法。通过编程实现了二维模型的非结构化三角形网格剖分,并编写了中心回线法瞬变电磁2．5维有限元正演程序。与前人计算结果对比,在取得相同计算精度的情况下,与结构化网格相比,非结构化网格所需网格和节点数量大大减少,计算效率更高。通过将非结构化网格法引入到瞬变电磁2．5维正演模拟中,实现了对复杂二维地电模型的有限元数值模拟,提高了现有有限元算法的应用范围。 相似文献
14.
井筒与油藏耦合的压裂水平井非稳态产能计算模型 总被引：1，自引：0，他引：1
为了准确预测压裂水平井的产量,综合考虑流体在储层、井壁、压裂裂缝、射孔孔眼以及水平井筒中的流动过程,根据势叠加原理,建立了考虑射孔孔眼、压裂裂缝生产干扰作用下的井筒、油藏耦合压裂水平井非稳态产能计算模型。结果表明:只考虑裂缝生产与裂缝、孔眼同时生产时的水平井产量在生产初期差别较大,生产后期差别变小;由于裂缝、孔眼的干扰,水平井筒两端裂缝的产量高于中间裂缝的产量;由于井筒摩阻作用,关于水平井筒中心对称的裂缝流量不相等;随着裂缝长度增加,裂缝、射孔孔眼、水平井产量增大;随着裂缝导流能力增加,裂缝产量增大,射孔孔眼产量降低,水平井产量增大;射孔段越长,孔眼与裂缝的干扰愈强,单个孔眼的产量降低,但由于射孔数增加,水平井产量增大;射孔深度增加,压裂水平井产量增加。 相似文献
15.
水平井分段多簇压裂是非常规油气藏开发的关键技术,在合理利用压裂诱导应力增大储层改造体积的同时,避免井间干扰所导致的裂缝砂堵和压裂窜扰,是压裂工艺优化中的关键科学问题。文章针对超临界CO2分段多簇压裂的缝间干扰和井间干扰问题,采用流固耦合的扩展有限元方法研究单井及多井裂缝诱导应力演化特征,充分考虑裂缝内超临界CO2流动和滤失,从非常规油气储层岩性特征、地应力场分布及施工工艺等多方面对压裂扰动应力进行系统研究,揭示单井分段多簇压裂缝扩展机制及应力扰动特征,在此基础上研究多井井间压裂缝干扰规律。结果表明：高水平应力差、高弹性模量的储层中压裂干扰界限较大,低水平应力差、低弹性模量地层需适度增大簇间距,减小簇间干扰;老井压裂后,其邻井压裂缝非对称系数随井间距呈先增大、后减小的趋势;当井间距等于压裂干扰界限时,非对称系数λ达到最大,且井周改造范围最大,但裂缝两翼的非对称性可能导致储层动用不充分。本研究为水平井细分切割压裂和立体式井网设计优化提供理论基础,在“双碳”战略背景下对非常规油气资源高效开发具有重要意义。 相似文献
16.
Knut–Andreas Lie Stein Krogstad Ingeborg Skjelkv?le Ligaarden Jostein Roald Natvig Halvor M?ll Nilsen B?rd Skaflestad 《Computational Geosciences》2012,16(2):297-322
Accurate geological modelling of features such as faults, fractures or erosion requires grids that are flexible with respect
to geometry. Such grids generally contain polyhedral cells and complex grid-cell connectivities. The grid representation for
polyhedral grids in turn affects the efficient implementation of numerical methods for subsurface flow simulations. It is
well known that conventional two-point flux-approximation methods are only consistent for K-orthogonal grids and will, therefore, not converge in the general case. In recent years, there has been significant research
into consistent and convergent methods, including mixed, multipoint and mimetic discretisation methods. Likewise, the so-called
multiscale methods based upon hierarchically coarsened grids have received a lot of attention. The paper does not propose
novel mathematical methods but instead presents an open-source Matlab? toolkit that can be used as an efficient test platform for (new) discretisation and solution methods in reservoir simulation.
The aim of the toolkit is to support reproducible research and simplify the development, verification and validation and testing
and comparison of new discretisation and solution methods on general unstructured grids, including in particular corner point
and 2.5D PEBI grids. The toolkit consists of a set of data structures and routines for creating, manipulating and visualising
petrophysical data, fluid models and (unstructured) grids, including support for industry standard input formats, as well
as routines for computing single and multiphase (incompressible) flow. We review key features of the toolkit and discuss a
generic mimetic formulation that includes many known discretisation methods, including both the standard two-point method
as well as consistent and convergent multipoint and mimetic methods. Apart from the core routines and data structures, the
toolkit contains add-on modules that implement more advanced solvers and functionality. Herein, we show examples of multiscale
methods and adjoint methods for use in optimisation of rates and placement of wells. 相似文献
17.
由于存在地层缺失、井眼穿层不断变化、重点层段测井曲线较少、井眼环境复杂、测井响应值与储层参数关系较为复杂等特点,页岩气储层水平井测井解释难度较大。针对该问题,本文采用“水平井钻遇地层模型—水平井间及水平井与直井响应差异分析—水平井曲线校正—基于直井的参数评价模型”的评价技术流程,进行页岩气水平井TOC测井评价。在测井评价中,针对响应规律复杂的特点,采用直方图法进行测井响应规律分析,同时考虑到水平井测井曲线较少,为了综合利用测井曲线信息,结合曲线重叠方法提出双重信息融合评价模型。结果表明,本方法在充分分析测井响应差异的基础上对水平井测井曲线进行校正,且提出的双重信息融合模型具有能充分利用测井信息、评价精度高以及操作简单等优点。 相似文献
18.
在分析苏里格气田辫状河体系及储层结构特征的基础上,提出了剖面储量集中度的概念,建立了厚层块状型、垂向叠置泛连通型、分散局部连通型、分散孤立型4种砂体分布模式,探讨了不同储层结构下的水平井采出程度,提出了水平井提高采收率技术对策。研究结果表明,辫状河沉积体系复合有效砂体由于“阻流带”的存在,直井动用不完善,水平井能克服“阻流带”的影响,提高层内储量动用程度;但由于砂体多层状分散分布,水平井开发会导致纵向含气层系储量动用不充分,影响层间采出程度。对于剖面储量集中度高的厚层块状型、垂向叠置泛连通型储层,采用水平井整体开发,Ⅰ+Ⅱ类井比例达70%以上,可显著提高储量动用程度和采收率。对于剖面储量集中度低的分散局部连通型储层,采用直井井网开发后进行甜点式优选水平井井位加密部署,可提高采收率10%以上。 相似文献
19.
水平井和大斜度井中阵列侧向测井响应数值模拟 总被引：1，自引：0，他引：1
针对水平井和大斜度井中阵列侧向电极系的工作原理,利用多电场叠加方式进行电场合成,采用三维有限元方法模拟仿真各个分电场的场分布,进而利用电场线性叠加原理得到阵列侧向测井响应。在基于计算机仿真的基础上,得到阵列侧向五条测井曲线的径向探测深度,阵列侧向径向探测深度要小于深侧向探测深度。考察了三维地层模型下井斜和侵入深度变化对阵列侧向测井响应的影响,分析了水平井和大斜度井中阵列侧向测井响应特征。模拟结果表明,在井斜小于15°时,阵列侧向测井响应受井斜影响小,可以不进行井斜校正;井斜超过60°的大斜度井以及水平井中,阵列侧向测井响应视地层厚度逐渐增大,测井响应值与直井条件下响应值差别较大,必须进行井斜校正。 相似文献
20.
A key ingredient in simulation of flow in porous media is accurate determination of the velocities that drive the flow. Large‐scale irregularities of the geology (faults, fractures, and layers) suggest the use of irregular grids in simulation. This paper presents a control‐volume mixed finite element method that provides a simple, systematic, easily implemented procedure for obtaining accurate velocity approximations on irregular (i.e., distorted logically rectangular) block‐centered quadrilateral grids. The control‐volume formulation of Darcy’s law can be viewed as a discretization into element‐sized “tanks” with imposed pressures at the ends, giving a local discrete Darcy law analogous to the block‐by‐block conservation in the usual mixed discretization of the mass‐conservation equation. Numerical results in two dimensions show second‐order convergence in the velocity, even with discontinuous anisotropic permeability on an irregular grid. The method extends readily to three dimensions. 相似文献