Beyond dual-porosity modeling for the simulation of complex flow mechanisms in shale reservoirs

The state of the art of modeling fluid flow in shale reservoirs is dominated by dual-porosity models which divide the reservoirs into matrix blocks that significantly contribute to fluid storage and fracture networks which principally control flow capacity. However, recent extensive microscopic studies reveal that there exist massive micro- and nano-pore systems in shale matrices. Because of this, the actual flow mechanisms in shale reservoirs are considerably more complex than can be simulated by the conventional dual-porosity models and Darcy’s law. Therefore, a model capturing multiple pore scales and flow can provide a better understanding of the complex flow mechanisms occurring in these reservoirs. This paper presents a micro-scale multiple-porosity model for fluid flow in shale reservoirs by capturing the dynamics occurring in three porosity systems: inorganic matter, organic matter (mainly kerogen), and natural fractures. Inorganic and organic portions of shale matrix are treated as sub-blocks with different attributes, such as wettability and pore structures. In kerogen, gas desorption and diffusion are the dominant physics. Since the flow regimes are sensitive to pore size, the effects of nano-pores and micro-pores in kerogen are incorporated into the simulator. The multiple-porosity model is built upon a unique tool for simulating general multiple-porosity systems in which several porosity systems may be tied to each other through arbitrary connectivities. This new model allows us to better understand complex flow mechanisms and eventually is extended into the reservoir scale through upscaling techniques. Sensitivity studies on the contributions of the different flow mechanisms and kerogen properties give some insight as to their importance. Results also include a comparison of the conventional dual-porosity treatment and show that significant differences in fluid distributions and dynamics are obtained with the improved multiple-porosity simulation.     

A semi-implicit method for incompressible three-phase flow in porous media

In this paper, we present a semi-implicit method for the incompressible three-phase flow equations in two dimensions. In particular, a high-order discontinuous Galerkin spatial discretization is coupled with a backward Euler discretization in time. We consider a pressure-saturation formulation, decouple the pressure and saturation equations, and solve them sequentially while still keeping each equation implicit in its respective unknown. We present several numerical examples on both homogeneous and heterogeneous media, with varying permeability and porosity. Our results demonstrate the robustness of the scheme. In particular, no slope limiters are required and a relatively large time step may be taken.     

Well placement optimization subject to realistic field development constraints

This work considers the well placement problem in reservoir management and field development optimization. In particular, it emphasizes embedding realistic and practical constraints into a mathematical optimization formulation. Such constraints are a prerequisite for the wider use of mathematical optimization techniques in well placement problems, since constraints are a way to incorporate reservoir engineering knowledge into the problem formulation. There are important design limitations that are used by the field development team when treating the well placement problem, and these limitations need to be articulated and eventually formalized within the problem before conducting the search for optimal well placements. In addition, these design limitations may be explicit or implicit. In this work, various design limitations pertaining to well locations have been developed in close collaboration with a field operator on the Norwegian Continental Shelf. Moreover, this work focuses on developing constraint-handling capability to enforce these various considerations during optimization. In particular, the Particle Swarm Optimization (PSO) algorithm is applied to optimize for the well locations, and various practical well placement constraints are incorporated into the PSO algorithm using two different constraint-handling techniques: a decoder procedure and the penalty method. The decoder procedure maps the feasible search space onto a cube and has the advantage of not requiring parameter tuning. The penalty method converts the constrained optimization problem into an unconstrained one by introducing an additional term, which is called a penalty function, to the objective function. In contrast to the penalty method, only feasible solutions are evaluated in the decoder method. Through numerical simulations, a comparison between the penalty method and the decoder technique is performed for two cases. We show that the decoder technique can easily be implemented for the well placement problem, and furthermore, that it performs better than the penalty method in most of the cases.     

Reliability Analysis of Earth Slopes Considering Spatial Variability

The paper presents a computational procedure for reliability analysis of earth slopes considering spatial variability of soils under the framework of the Limit Equilibrium Method. In the reliability analysis of earth slopes, the effect of spatial variability of soil properties is generally included indirectly by assuming that the probabilistic critical slip surface is the same as that determined without considering spatial variability. In contrast to this indirect approach, in the direct approach, the effect of spatial variability is included in the process of determination of the probabilistic critical surface itself. While the indirect approach requires much less computational effort, the direct approach is definitely more rigorous. In this context this paper attempts to investigate, with the help of numerical examples, how far away are the results obtained from the indirect approach from that obtained from the direct approach. In both the approaches, it is required to use a model of discretization of random fields into finite random variables. A few such models are available in the literature for one-dimensional (1D) as well as two-dimensional (2D) spatial variability. The developed computational scheme is based on the First Order Reliability Method (FORM) coupled with the Spencer Method of Slices valid for limit equilibrium analysis of general slip surfaces. The study includes bringing out the computational advantages and disadvantages of the three commonly used discretization models. The sensitivity of the reliability index to the magnitudes of the scales of fluctuation has also been studied.     

Stabilization of Expansive Soil with Cement Kiln Dust and RBI Grade 81 at Subgrade Level

The present research work deals with an expansive high plastic clayey soil with cement kiln dust (CKD) and stabilizer (RBI Grade 81). The physical and engineering properties of soil are plasticity, compaction, unconfined compressive strength (UCS), consolidation and California bearing ratio (CBR) of the clayey soil and clay treated with CKD and stabilizer were determined. Soil chemistry was examined before and after treatment using scanning electron microscope (SEM) and elemental dispersive spectrometer. The clay mixed with CKD, CKD and RBI Grade 81 was found that optimum contents are 10 % (CKD), 15 % CKD with 4 % RBI Grade 81, respectively. The result indicates that CKD alone will decrease maximum dry density and increase optimum moisture content. CKD with RBI Grade 81 slightly increases maximum dry density and decreases optimum moisture content. UCS increased with CKD alone and CKD with RBI Grade 81 from 88.3 to 976 kN/m², respectively. CBR values were increased by the addition of CKD, CKD with RBI Grade 81 from 1.65 to 21.7 %. With the curing time of 3, 14 and 28 days, UCS and CBR values were increased due to pozzolanic reaction from cementations material. The treated soil has considerable reduction in compression index. SEM images clearly indicate the formation of CSH and CAH gel.     

The Blasting Test and Blasting Vibration Monitoring of Vertical Crater Retreat Mining Method in the Luohe Iron Mine

Because the land requisition and demolishing became difficult more and more, the mining scheme of Luohe iron mine was changed from caving method to filling method. In order to ensure the safety of the residence and the underground tunnel cavern within the mobile belt of the underground mining, the Luohe iron mine did the blasting test on the vertical crater retreat mining method and blasting vibration monitoring. The blasting experiments use common emulsified oil explosives and non-electric initiation system. The way of caved ore adopts the cutting groove and bench side. The NUBOX-6016 intelligent vibration monitor was chosen in the blasting vibration monitoring. Twice experiments on the blasting vibration monitoring were done on the surface or in the underground refuge cavern. The first test select the three monitoring points on the ground and the second select two monitoring points on the ground and a monitoring point in the underground refuge cavern. The blasting vibration monitoring data were conducted by the regression analysis in the Sodev’s empirical formula. The vibration attenuation formula about the underground blasting vibration transmitting in three directions is derived. The blasting test on the vertical crater retreat mining method and the blasting vibration were analyzed. It is estimated if the vibration damage possibly the surface buildings and related facilities of mine.     

Study on Strength Influence Mechanism of Fiber-Reinforced Expansive Soil Using Jute

This study was undertaken to research the effects of jute fiber content, fiber length, water content and dry density of reinforced and unreinforced soil on the strength influence mechanism by implementing a series of laboratory tests and analysis. The most efficient fiber reinforcement effects was achieved by means of adding jute fiber with content of 0.6 % and length of 6 mm into expansive soil specimen prepared at maximum dry density and optimum moisture content. The cohesion of reinforced specimens increased first with increasing fiber content and fiber length and then decreased with further increase in fiber content and fiber length. The internal friction angle of reinforced specimens were not affected significantly by fiber content and fiber length. Higher water content reduces the fiber reinforcement effects by means of acting as lubricant in the interface of fiber and soil particles. Fiber reinforcement effects is more prominent for specimens prepared at higher dry density by increasing the effective contact area of fiber/soil. The application prospect of soil reinforcement using natural fiber is impeded by the hydrophilic nature and biodegradability of natural fiber, thus, studies on using chemical additive to do surface treatment for natural fiber are needed to improve the interfacial interaction of fiber/soil so as to widen the application of natural fiber.     

Two Methods for Predicting Reliability Index and Critical Probabilistic Slip Surface of Soil Slopes

Two methods of reliability analysis of soil slopes are studied, and the representative flow charts of both methods are illustrated. Method 1 can predict the reliability index and the critical probabilistic slip surface directly and it is computational efficient, but it needs the development of new codes for integrating the reliability analysis code and the slope stability code. Method 2 makes the reliability analysis code call the slope stability analysis code directly, and each code can be considered as an intact part. The main result of Method 2 is the reliability index of soil slope. Combined with the proposed method for locating the critical slip surface, Method 2 can also predict the probabilistic slip surface. Although Method 2 needs much more callings of the subprogram of slope stability analysis code, it needs not the developing of new computer program. Thus, Method 2 is easy to use and can be applied to different reliability analysis methods and slope stability analysis methods.     

Evaluation of Subsurface Lithologic Units for Groundwater Exploration in the Barind Tract of Bangladesh

Barind Tract, located north western part of Bangladesh, is one of the most diversified physiographic units of the country. The surface water supply in this part is particularly limited, so the irrigation is almost entirely depends on groundwater. However, over exploitation indicates falling groundwater heads in this area. The objective of present study is to examine the nature of the aquifer system Barind Tract of Bangladesh in order to assess the sustainability of groundwater yield. Borehole lithology data were collected, processed and analyzed for this purpose. Representative panel diagram, 3-D stratigraphic and cross-sectional views were also prepared for necessary assessment of the variation of individual subsurface stratum in different locations. The study identified three subsurface geologic formations namely, a top clay layer, sand layer of different grain size and at the bottom an impermeable clay zone. Maps of formation thickness and index revealed that aquifer thickness is low in the north-western corner and in some places of middle of south-western corner. The thickness of sand formation in other places is recorded above 20 m. It is expect that the finding of the study will help groundwater resources development, planning and management in the area.