共查询到20条相似文献，搜索用时 31 毫秒
1.
Anna Nissen Eirik Keilegavlen Tor Harald Sandve Inga Berre Jan Martin Nordbotten 《Computational Geosciences》2018,22(2):451-467
In simulation of fluid injection in fractured geothermal reservoirs, the characteristics of the physical processes are severely affected by the local occurence of connected fractures. To resolve these structurally dominated processes, there is a need to develop discretization strategies that also limit computational effort. In this paper, we present an upscaling methodology for geothermal heat transport with fractures represented explicitly in the computational grid. The heat transport is modeled by an advection-conduction equation for the temperature, and solved on a highly irregular coarse grid that preserves the fracture heterogeneity. The upscaling is based on different strategies for the advective term and the conductive term. The coarse scale advective term is constructed from sums of fine scale fluxes, whereas the coarse scale conductive term is constructed based on numerically computed basis functions. The method naturally incorporates the coupling between solution variables in the matrix and in the fractures, respectively, via the discretization. In this way, explicit transfer terms that couple fracture and matrix solution variables are avoided. Numerical results show that the upscaling methodology performs well, in particular for large upscaling ratios, and that it is applicable also to highly complex fracture networks. 相似文献
2.
Multiscale mixed/mimetic methods on corner-point grids 总被引：1，自引：0，他引：1
Multiscale simulation is a promising approach to facilitate direct simulation of large and complex grid models for highly
heterogeneous petroleum reservoirs. Unlike traditional simulation, approaches based on upscaling/downscaling, multiscale methods
seek to solve the full flow problem by incorporating subscale heterogeneities into local discrete approximation spaces. We
consider a multiscale formulation based on a hierarchical grid approach, where basis functions with subgrid resolution are
computed numerically to correctly and accurately account for subscale variations from an underlying (fine-scale) geomodel
when solving the global flow equations on a coarse grid. By using multiscale basis functions to discretise the global flow
equations on a (moderately sized) coarse grid, one can retain the efficiency of an upscaling method and, at the same time,
produce detailed and conservative velocity fields on the underlying fine grid. For pressure equations, the multiscale mixed
finite-element method (MsMFEM) has been shown to be a particularly versatile approach. In this paper, we extend the method
to corner-point grids, which is the industry standard for modelling complex reservoir geology. To implement MsMFEM, one needs
a discretisation method for solving local flow problems on the underlying fine grids. In principle, any stable and conservative
method can be used. Here, we use a mimetic discretisation, which is a generalisation of mixed finite elements that gives a
discrete inner product, allows for polyhedral elements, and can (easily) be extended to curved grid faces. The coarse grid
can, in principle, be any partition of the subgrid, where each coarse block is a connected collection of subgrid cells. However,
we argue that, when generating coarse grids, one should follow certain simple guidelines to achieve improved accuracy. We
discuss partitioning in both index space and physical space and suggest simple processing techniques. The versatility and
accuracy of the new multiscale mixed methodology is demonstrated on two corner-point models: a small Y-shaped sector model
and a complex model of a layered sedimentary bed. A variety of coarse grids, both violating and obeying the above mentioned
guidelines, are employed. The MsMFEM solutions are compared with a reference solution obtained by direct simulation on the
subgrid. 相似文献
3.
The role of shear dilation as a mechanism of enhancing fluid flow permeability in naturally fractured reservoirs was mainly recognized in the context of hot dry rock (HDR) geothermal reservoir stimulation. Simplified models based on shear slippage only were developed and their applications to evaluate HDR geothermal reservoir stimulation were reported. Research attention is recently focused to adjust this stimulation mechanism for naturally fractured oil and gas reservoirs which reserve vast resources worldwide. This paper develops the overall framework and basic formulations of this stimulation model for oil and gas reservoirs. Major computational modules include: natural fracture simulation, response analysis of stimulated fractures, average permeability estimation for the stimulated reservoir and prediction of an average flow direction. Natural fractures are simulated stochastically by implementing ‘fractal dimension’ concept. Natural fracture propagation and shear displacements are formulated by following computationally efficient approximate approaches interrelating in situ stresses, natural fracture parameters and stimulation pressure developed by fluid injection inside fractures. The average permeability of the stimulated reservoir is formulated as a function of discretized gridblock permeabilities by applying cubic law of fluid flow. The average reservoir elongation, or the flow direction, is expressed as a function of reservoir aspect ratio induced by directional permeability contributions. The natural fracture simulation module is verified by comparing its results with observed microseismic clouds in actual naturally fractured reservoirs. Permeability enhancement and reservoir growth are characterized with respect to stimulation pressure, in situ stresses and natural fracture density applying the model to two example reservoirs. Copyright © 2002 John Wiley & Sons, Ltd. 相似文献
4.
Coupled hydro-mechanical (HM) processes are significant in geological engineering such as oil and gas extraction, geothermal energy, nuclear waste disposal and for the safety assessment of dam foundations and rock slopes, where the geological media usually consist of fractured rock masses. In this study, we developed a model for the analysis of coupled hydro-mechanical processes in porous rock containing dominant fractures, by using the numerical manifold method (NMM). In the current model, the fractures are regarded as different material domains from surrounding rock, i.e., finite-thickness fracture zones as porous media. Compared with the rock matrix, these fractured porous media are characterized with nonlinear behavior of hydraulic and mechanical properties, involving not only direct (poroelastic) coupling but also indirect (property change) coupling. By combining the potential energy associated with mechanical responses, fluid flow and solid–fluid interactions, a new formulation for direct HM coupling in porous media is established. For indirect coupling associated with fracture opening/closure, we developed a new approach implicitly considering the nonlinear properties by directly assembling the corresponding strain energy. Compared with traditional methods with approximation of the nonlinear constitutive equations, this new formulation achieves a more accurate representation of the nonlinear behavior. We implemented the new model for coupled HM analysis in NMM, which has fixed mathematical grid and accurate integration, and developed a new computer code. We tested the code for direct coupling on two classical poroelastic problems with coarse mesh and compared the results with the analytical solutions, achieving excellent agreement, respectively. Finally, we tested for indirect coupling on models with a single dominant fracture and obtained reasonable results. The current poroelastic NNM model with a continuous finite-thickness fracture zone will be further developed considering thin fractures in a discontinuous approach for a comprehensive model for HM analysis in fractured porous rock masses. 相似文献
5.
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. 相似文献
6.
水力压裂是青海共和盆地干热岩地热资源开发的难点技术问题之一。本文基于升级改造的大尺寸真三轴水力压裂物理模拟实验系统模拟干热岩储层高温高压环境,利用青海共和盆地露头岩心进行水力压裂物理模拟实验,揭示干热岩储层水力裂缝的起裂和扩展规律。通过物理模拟实验发现：干热岩储层裂缝起裂可以通过文中提出的起裂模型判断起裂方式和预测起裂压力;水力裂缝在岩石基质中的扩展形态简单,仅沿最大主应力方向延伸;但是水力裂缝会受到岩石中弱面的影响,发生转向沿弱面延伸,形成较复杂的裂缝形态。因此,建议在干热岩储层实际施工中,在天然裂缝发育较丰富的层段开展水力压裂,以实现复杂裂缝网络提取地热能。 相似文献
7.
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. 相似文献
8.
岩溶化裂隙岩体是普遍发育于自然界中具有初始损伤的岩体。为了研究岩溶化裂隙岩体损伤破坏特征,本文以贵州某地赋存的溶蚀岩体为研究对象,运用损伤力学理论构建岩溶化裂隙岩体在单轴压缩条件下的损伤演化模型,并建立岩溶化裂隙岩体损伤演化方程。采用颗粒流数值软件进行单轴压缩数值试验,进一步研究岩溶化裂隙岩体试件在单轴压缩条件下的损伤演化特征,分析岩溶化裂隙岩体的微观损伤特征。结果表明：岩溶化裂隙岩体的初始损伤主要包括溶蚀损伤和裂隙损伤。岩溶化裂隙岩体的初始损伤随着溶蚀率的增加而增加,最终增加速率趋于平缓;岩溶化裂隙岩体的损伤演化曲线均呈“S”型变化,先缓慢增加,再迅速增加,最后缓慢增加至损伤值1;岩溶化裂隙岩体存在异构特征,导致破坏裂隙起源于具有初始损伤的溶蚀孔洞和裂隙处,随后裂隙经历萌发、扩张和剪切作用、数量和长度增加以及裂隙贯通4个阶段后发生宏观破坏。 相似文献
9.
We present a fully implicit formulation of coupled flow and geomechanics for fractured three-dimensional subsurface formations. The Reservoir Characterization Model (RCM) consists of a computational grid, in which the fractures are represented explicitly. The Discrete Fracture Model (DFM) has been widely used to model the flow and transport in natural geological porous formations. Here, we extend the DFM approach to model deformation. The flow equations are discretized using a finite-volume method, and the poroelasticity equations are discretized using a Galerkin finite-element approximation. The two discretizations—flow and mechanics—share the same three-dimensional unstructured grid. The mechanical behavior of the fractures is modeled as a contact problem between two computational planes. The set of fully coupled nonlinear equations is solved implicitly. The implementation is validated for two problems with analytical solutions. The methodology is then applied to a shale-gas production scenario where a synthetic reservoir with 100 natural fractures is produced using a hydraulically fractured horizontal well. 相似文献
10.
This paper presents a novel mass-conservative mixed multiscale method for solving flow equations in heterogeneous porous media. The media properties (the permeability) contain multiple scales and high contrast. The proposed method solves the flow equation in a mixed formulation on a coarse grid by constructing multiscale basis functions. The resulting velocity field is mass-conservative on the fine grid. Our main goal is to obtain first-order convergence in terms of the mesh size which is independent of local contrast. This is achieved, first, by constructing some auxiliary spaces, which contain global information that cannot be localized, in general. This is built on our previous work on the generalized multiscale finite element method (GMsFEM). In the auxiliary space, multiscale basis functions corresponding to small (contrast-dependent) eigenvalues are selected. These basis functions represent the high-conductivity channels (which connect the boundaries of a coarse block). Next, we solve local problems to construct multiscale basis functions for the velocity field. These local problems are formulated in the oversampled domain, taking into account some constraints with respect to auxiliary spaces. The latter allows fast spatial decay of local solutions and, thus, allows taking smaller oversampled regions. The number of basis functions depends on small eigenvalues of the local spectral problems. Moreover, multiscale pressure basis functions are needed in constructing the velocity space. Our multiscale spaces have a minimal dimension, which is needed to avoid contrast dependence in the convergence. The method’s convergence requires an oversampling of several layers. We present an analysis of our approach. Our numerical results confirm that the convergence rate is first order with respect to the mesh size and independent of the contrast. 相似文献
11.
K.-A. Lie O. Møyner J. R. Natvig A. Kozlova K. Bratvedt S. Watanabe Z. Li 《Computational Geosciences》2017,21(5-6):981-998
For the past 10 years or so, a number of so-called multiscale methods have been developed as an alternative approach to upscaling and to accelerate reservoir simulation. The key idea of all these methods is to construct a set of prolongation operators that map between unknowns associated with cells in a fine grid holding the petrophysical properties of the geological reservoir model and unknowns on a coarser grid used for dynamic simulation. The prolongation operators are computed numerically by solving localized flow problems, much in the same way as for flow-based upscaling methods, and can be used to construct a reduced coarse-scale system of flow equations that describe the macro-scale displacement driven by global forces. Unlike effective parameters, the multiscale basis functions have subscale resolution, which ensures that fine-scale heterogeneity is correctly accounted for in a systematic manner. Among all multiscale formulations discussed in the literature, the multiscale restriction-smoothed basis (MsRSB) method has proved to be particularly promising. This method has been implemented in a commercially available simulator and has three main advantages. First, the input grid and its coarse partition can have general polyhedral geometry and unstructured topology. Secondly, MsRSB is accurate and robust when used as an approximate solver and converges relatively fast when used as an iterative fine-scale solver. Finally, the method is formulated on top of a cell-centered, conservative, finite-volume method and is applicable to any flow model for which one can isolate a pressure equation. We discuss numerical challenges posed by contemporary geomodels and report a number of validation cases showing that the MsRSB method is an efficient, robust, and versatile method for simulating complex models of real reservoirs. 相似文献
12.
In this paper, a multiscale homogenization approach is developed for fully coupled saturated porous media to represent the idealized sugar cube model, which is generally employed in fractured porous media on the basis of dual porosity models. In this manner, an extended version of the Hill-Mandel theory that incorporates the microdynamic effects into the multiscale analysis is presented, and the concept of the deformable dual porosity model is demonstrated. Numerical simulations are performed employing the multiscale analysis and dual porosity model, and the results are compared with the direct numerical simulation through 2 numerical examples. Finally, a combined multiscale-dual porosity technique is introduced by employing a bridge between these 2 techniques as an alternative approach that reduces the computational cost of numerical simulation in modeling of heterogeneous deformable porous media. 相似文献
13.
A computational procedure for simulating heat conduction in fractured rock masses is proposed and illustrated within the context of the finite element method. The procedure makes use of simple local models for conduction in the vicinity of a single open fracture. The local models allow effective thermal properties for fractured rock to be expressed in terms of fracture characteristics such as spacing, orientation and aperture. The distributions of fractures and fracture properties within the finiteelement model are derived from a statistical representation of geological field data. The procedure is demonstrated to satisfactorily reproduce heat transfer results across single fractures and provide an efficient means of representing conduction in formations with multiple fractures. 相似文献
14.
Hydraulic fracturing is the method of choice to enhance reservoir permeability and well efficiency for extraction of shale gas. Multi‐stranded non‐planar hydraulic fractures are often observed in stimulation sites. Non‐planar fractures propagating from wellbores inclined from the direction of maximum horizontal stress have also been reported. The pressure required to propagate non‐planar fractures is in general higher than in the case of planar fractures. Current computational methods for the simulation of hydraulic fractures generally assume single, symmetric, and planar crack geometries. In order to better understand hydraulic fracturing in complex‐layered naturally fractured reservoirs, fully 3D models need to be developed. In this paper, we present simulations of 3D non‐planar fracture propagation using an adaptive generalized FEM. This method greatly facilitates the discretization of complex 3D fractures, as finite element faces are not required to fit the crack surfaces. A solution strategy for fully automatic propagation of arbitrary 3D cracks is presented. The fracture surface on which pressure is applied is also automatically updated at each step. An efficient technique to numerically integrate boundary conditions on crack surfaces is also proposed and implemented. Strongly graded localized refinement and analytical asymptotic expansions are used as enrichment functions in the neighborhood of fracture fronts to increase the computational accuracy and efficiency of the method. Stress intensity factors with pressure on crack faces are extracted using the contour integral method. Various non‐planar crack geometries are investigated to demonstrate the robustness and flexibility of the proposed simulation methodology. Copyright © 2014 John Wiley & Sons, Ltd. 相似文献
15.
In this paper, we develop a multiscale model reduction technique that describes shale gas transport in fractured media. Due to the pore-scale heterogeneities and processes, we use upscaled models to describe the matrix. We follow our previous work (Akkutlu et al. Transp. Porous Media 107(1), 235–260, 2015), where we derived an upscaled model in the form of generalized nonlinear diffusion model to describe the effects of kerogen. To model the interaction between the matrix and the fractures, we use Generalized Multiscale Finite Element Method (Efendiev et al. J. Comput. Phys. 251, 116–135, 2013, 2015). In this approach, the matrix and the fracture interaction is modeled via local multiscale basis functions. In Efendiev et al. (2015), we developed the GMsFEM and applied for linear flows with horizontal or vertical fracture orientations aligned with a Cartesian fine grid. The approach in Efendiev et al. (2015) does not allow handling arbitrary fracture distributions. In this paper, we (1) consider arbitrary fracture distributions on an unstructured grid; (2) develop GMsFEM for nonlinear flows; and (3) develop online basis function strategies to adaptively improve the convergence. The number of multiscale basis functions in each coarse region represents the degrees of freedom needed to achieve a certain error threshold. Our approach is adaptive in a sense that the multiscale basis functions can be added in the regions of interest. Numerical results for two-dimensional problem are presented to demonstrate the efficiency of proposed approach. 相似文献
16.
Liao Jianxing Gou Yang Feng Wentao Mehmood Faisal Xie Yachen Hou Zhengmeng 《Acta Geotechnica》2020,15(2):279-295
Although hydraulic fracturing has been massively studied and applied as a key technique to enhance the gas production from tight formations, some problems and uncertainties exist to accurately predict and analyze the fracture behavior in complex reservoirs, especially in the naturally fractured reservoirs like shale reservoirs. This paper presents a full 3D numerical model (FLAC3D) to study hydraulic fracturing behavior under the impact of preexisting orthogonal natural fractures. In this numerical model, the hydraulic fracture propagation direction is assumed perpendicular to the minimum principal stress and activated only by tensile failure, whereas the preexisting natural fractures can be activated by tensile or shear failure or a combination of them, and only tensile failure can open the natural fracture as well. The newly developed model was used to study the impact of preexisting orthogonal natural fractures on hydraulic fracturing behavior, based on a multistage hydraulic fracturing operation in a naturally fractured reservoir from the Barnett Shale formation, northwest of Texas in USA. In this multistage operation, two more representative stages, i.e., stage 1 with a relatively large horizontal stress anisotropy of 3.3 MPa and stage 4 with a comparatively small one of 1.3 MPa, were selected to conduct the simulation. Based on the numerical results, one can observe that the interaction between hydraulic and natural fracture is driven mainly by induced stress around fracture tip. Besides, the horizontal stress anisotropy plays a key role in opening the natural fracture. Thus, no significant opened fracture is activated on natural fracture in stage 1, while in stage 4 an opened fracture invades to about 90 m into the first natural fracture. Conversely, the hydraulic fracture length in stage 1 is much longer than in stage 4, as some fluid volume is stored in the opened natural fracture in stage 4. In this work, the shear failure on natural fractures is treated as the main factor for inducing the seismic events. And the simulated seismic events, i.e., shear failure on natural fractures, are very comparable with the measured seismic events.
相似文献17.
This paper presents a fracture mapping (FM) approach combined with the extended finite element method (XFEM) to simulate coupled deformation and fluid flow in fractured porous media. Specifically, the method accurately represents the impact of discrete fractures on flow and deformation, although the individual fractures are not part of the finite element mesh. A key feature of FM‐XFEM is its ability to model discontinuities in the domain independently of the computational mesh. The proposed FM approach is a continuum‐based approach that is used to model the flow interaction between the porous matrix and existing fractures via a transfer function. Fracture geometry is defined using the level set method. Therefore, in contrast to the discrete fracture flow model, the fracture representation is not meshed along with the computational domain. Consequently, the method is able to determine the influence of fractures on fluid flow within a fractured domain without the complexity of meshing the fractures within the domain. The XFEM component of the scheme addresses the discontinuous displacement field within elements that are intersected by existing fractures. In XFEM, enrichment functions are added to the standard finite element approximation to adequately resolve discontinuous fields within the simulation domain. Numerical tests illustrate the ability of the method to adequately describe the displacement and fluid pressure fields within a fractured domain at significantly less computational expense than explicitly resolving the fracture within the finite element mesh. Copyright © 2013 John Wiley & Sons, Ltd. 相似文献
18.
Knut-Andreas Lie Jostein R. Natvig Stein Krogstad Yahan Yang Xiao-Hui Wu 《Computational Geosciences》2014,18(3-4):357-372
A Dirichlet–Neumann representation method was recently proposed for upscaling and simulating flow in reservoirs. The DNR method expresses coarse fluxes as linear functions of multiple pressure values along the boundary and at the center of each coarse block. The number of flux and pressure values at the boundary can be adjusted to improve the accuracy of simulation results and, in particular, to resolve important fine-scale details. Improvement over existing approaches is substantial especially for reservoirs that contain high-permeability streaks or channels. As an alternative, the multiscale mixed finite-element (MsMFE) method was designed to obtain fine-scale fluxes at the cost of solving a coarsened problem, but can also be used as upscaling methods that are flexible with respect to geometry and topology of the coarsened grid. Both methods can be expressed in mixed-hybrid form, with local stiffness matrices obtained as “inner products” of numerically computed basis functions with fine-scale sub-resolution. These basis functions are determined by solving local flow problems with piecewise linear Dirichlet boundary conditions for the DNR method and piecewise constant Neumann conditions for MsMFE. Adding discrete pressure points in the DNR method corresponds to subdividing faces in the coarse grid and hence increasing the number of basis functions in the MsMFE method. The methods show similar accuracy for 2D Cartesian cases, but the MsMFE method is more straightforward to formulate in 3D and implement for general grids. 相似文献
19.
随着新时代对绿色环保资源的需求日渐旺盛,地热资源即将进入大规模开发阶段。为经济、高效获取地热资源,目前主要应用的是PDC钻头。但是,部分钻取地热井存在遭遇高温硬地层的情况,常规PDC钻头存在耐高温性能差、硬度略低、耐磨性不足等问题。本文总结了深部高温硬地层钻进时PDC切削齿受到磨粒磨损,冲击磨损,热损伤等3种磨损机理,介绍了4种提高PDC切削齿耐高温性,硬度及耐磨性的方法：涂层技术、改善钻进参数,复合钻头和钎焊切削齿。根据目前研究情况,建议单个PDC切削齿实验可以从实现独立分析与切削齿关联性相结合方面展开,通过动态建模和工作状态实时模拟优化设计方案。 相似文献
20.
Researchers have recently realized that the non-tectonic natural fractures are developed in shale formations and significant for the exploitation of shale gas. Studies have shown that the tectonic fractures in naturally fractured reservoirs have influences on the maximization of stimulated reservoir volume (SRV) during hydraulic fracturing. However, the effect of the non-tectonic randomly natural fractures on the fracturing network propagation is not well understood. Laboratory experiments are proposed to study the evolution of fracturing network in naturally fractured formations with specimens that contain non-tectonic random fractures. The influences of the dominating factors were studied and analyzed, with an emphasis on natural fracture density, stress ratio, and injection rate. The response surface methodology was employed to perform the multiple-factor analysis and optimization in the maximization of the SRV. A sensitivity study reveals a number of interesting observations resulting from these parameters on the fracturing network evaluation. It is suggested from the geometry morphology of fracturing network that high natural fracture density and injection rate tend to maximize the fracturing network. The influence of stress contrast on fracturing network is nonlinear; an optimal value exists resulting in the best hydraulic fracturing effectiveness. 相似文献