scholarly journals Numerical Simulation on Mesoscale Mechanism of Seepage in Coal Fractures by Fluid-Sloid Coupling Method

Geofluids ◽  
2021 ◽  
Vol 2021 ◽  
pp. 1-11
Author(s):  
Kai Si ◽  
Ruidong Peng ◽  
Leilei Zhao ◽  
Yan Zhao ◽  
Yaheng Zhu ◽  
...  

Trying to reveal the mechanism of gas seepage in coal is of significance to both safe mining and methane exploitation. A series of FEM numerical models were built up and studied so as to explore the mesoscale mechanism of seepage in coal fractures. The proposed mesoscale FEM model is a cube with micron fractures along three orthogonal directions. The distribution of velocity and pressure under fluid-solid coupling was obtained, and furthermore, the seepage flow flux and an equivalent permeability of the whole model were calculated. The influences of fracture width, outlet velocity, and in situ stress level on seepage were investigated. The numerical results show that nonlinear Darcy seepage occurs during low velocity zone. The permeability is increased linearly with the increasing of facture width and outlet velocity. A certain change of lateral coefficient of in situ stress also affects seepage. The permeability is increased sharply once deviating the isotropic spherical stress state, but it is no longer changed obviously after the lateral coefficient has been increased or decreased more than 20%. The mesoscale seepage mechanism in coal fractures has been preliminarily revealed by considering fluid-solid coupling effect, and the key factors influencing fluid seepage in coal fractures were demonstrated. The proposed methods and results will be helpful to the further study of seepage behaviour in coal with more complex structures.

2011 ◽  
Vol 90-93 ◽  
pp. 2033-2036 ◽  
Author(s):  
Jin Shan Sun ◽  
Hong Jun Guo ◽  
Wen Bo Lu ◽  
Qing Hui Jiang

The factors affecting the TBM tunnel behavior in jointed rock mass is investigated. In the numerical models the concrete segment lining of TBM tunnel is concerned, which is simulated as a tube neglecting the segment joint. And the TBM tunnel construction process is simulate considering the excavation and installing of the segment linings. Some cases are analyzed with different joint orientation, joint spacing, joint strength and tunnel depth. The results show that the shape and areas of loosing zones of the tunnel are influenced by the parameters of joint sets and in-situ stress significantly, such as dip angle, spacing, strength, and the in-situ stress statement. And the stress and deformation of the tunnel lining are influenced by the parameters of joint sets and in-situ stress, too.


2012 ◽  
Vol 616-618 ◽  
pp. 538-542 ◽  
Author(s):  
Fu Xiang Zhang ◽  
Wei Feng Ge ◽  
Xiang Tong Yang ◽  
Wei Zhang ◽  
Jian Xin Peng

To alleviate the problems of casing collapse induced by the coupling effect of rock salt creep and casing wear, the effects of salt creep, attrition rate and casing abrasive position on the equivalent stress on casings in non-uniform in-situ stress field is analyzed by finite-difference model with worn casing, cement and salt formation. It indicates that, creep reduces the yield strength of worn casing to a certain extent; Equivalent stress on casings is bigger and more non-uniform when the abrasion is more serious; Wear position obviously changes the distribution of equivalent stress on casing, and when the wear located along the direction of the minimum in-situ stress, equivalent stress on casing could be the largest that leads to the casing being failed more easily. Equivalent stress on casings increases gradually with creep time increasing and will get to balance in one year or so; In addition, new conclusions are obtained which are different from before: the maximum equivalent stress on casings is in the direction of the minimum horizontal stress, only when the attrition rate of the casing is little; otherwise, it is not. This method could help to improve the wear prediction and design of casings.


2011 ◽  
Vol 90-93 ◽  
pp. 2332-2338
Author(s):  
Hai Jun Wang ◽  
Xu Hua Ren ◽  
Ji Xun Zhang

Deformation memory effect (DME) is one of the rock memory effects. One important application of the DME is to determine the in situ stress. The in situ stress measurement methods based on the DME are commercial and permit large number of measurements. Application of DME needs the reliable and effective identification methods to detect the DME. However, the existing identification methods are insufficient and not distinct. In this paper, a new method based on fractal dimension was proposed. It takes advantage of the increase of the irregularity of the stress-strain curve after the previous maximum stress attained. Numerical models for the sandstone and granite were developed based on the contact bond model in PFC2D. Fractal dimension method was employed to detect the DME for two types of the rock. The results demonstrate that the fractal dimension method is effective and reliable in the identification of DME in rock.


Geofluids ◽  
2020 ◽  
Vol 2020 ◽  
pp. 1-13 ◽  
Author(s):  
Guannan Liu ◽  
Dayu Ye ◽  
Feng Gao ◽  
Jishan Liu

In the process of coalbed methane exploitation, the fracture and pore structure is the key problem that affects the permeability of coalbed. At present, the coupling effect of fracture and pore structure and in situ stress is seldom considered in the study of coal seam permeability. In this paper, the fractal seepage model is coupled with coal deformation, and the adsorption expansion effect is considered. A multifield coupling model considering the influence of matrix and fracture structure is established. Then, the influence of pore structure parameters of main fracture on macropermeability is analyzed, including (1) fractal dimension of fracture length, (2) maximum fracture length, (3) fractal dimension of throat diameter, and (4) fractal dimension of throat bending. At the same time, the simulation results are compared with the results of Darcy’s uniform permeability model. The results show that the permeability calculated by the proposed model is significantly different from that calculated by the traditional cubic model. Under the action of in situ stress, when the porosity and other parameters remain unchanged, the macropermeability of coal is in direct proportion to the fractal dimension of coal fracture length, the fractal dimension of throat diameter, and the maximum fracture length and in inverse proportion to the fractal dimension of coal throat curvature.


Author(s):  
Andrzej Adam Truty

Abstract Numerical modeling of deep excavations becomes a standard practice in modern geotechnical engineering. A detailed numerical model for a given case is able to reproduce major effects of soil-structure interaction by taking into account any kind of drainage conditions, strong stiffness variation due to effective stress and strain changes, creep and cracking, when reinforced concrete is used as a structural material, but also interface effects between subsoil and structure. Calibrating soil constitutive models is one of the most difficult tasks and due to several sources of uncertainty there is no one unique set of the data that should be used in numerical predictions. Lack or incompleteness of experimental data, significant mismatch between laboratory and field tests is an another source of difficulty. Contrary to several simplified methods, that are usually limited to two dimensions, numerical models allow a full 3D analysis in which many simplifications can be eliminated. This paper is devoted to the problem of in situ stress disturbance caused by diaphragm wall installation in overconsolidated quaternary sandy clays and its influence on final wall deformations.


2021 ◽  
Author(s):  
Yi Luo ◽  
Hangli Gong ◽  
Dengxing Qu ◽  
Xinping Li ◽  
Shaohua Hu ◽  
...  

Abstract The escape of toxic and harmful gases is a common disaster effect in tunnel engineering. Frequent drilling and blasting excavation disturbances under high in-situ stress environment will inevitably lead to cumulative damage effect on surrounding rock, which will increase the permeability coefficient of surrounding rock, increase the risk of toxic and harmful gas escape, and seriously endanger construction safety. In this paper, based on real-time monitoring data of harmful gases during blasting and excavation of Yuelongmen Tunnel on Chengdu-Lanzhou Railway, this study summarized laws and distribution characteristics of harmful gas escape intensified by the blasting excavation, and the effectiveness of shotcreting and grouting for water blocking to inhibit gas escape is verified. Then, taking water-containing and gas-containing voids as carriers, considering the influence of different in-situ stress, explosion load and void parameters (including void pressure, void diameter and distance between void and tunnel), to carry out research on the escape mechanism of water-soluble (H 2 S) and insoluble (CH 4 ) toxic and harmful gases under the coupling effect of stress-seepage-damage. The relationship between the amount of harmful gas escaped and the damage degree of the surrounding rock of the tunnel is analyzed, and the functional relationship between it and the in-situ stress, explosion load and cave parameters is established. The results further demonstrate that the amount of escaped harmful gases, such as methane and H 2 S is closely related to lithology of surrounding rock, occurrence conditions of the deep rock mass, development degree of structural fractures and void parameters. The damage of surrounding rock caused by dynamic disturbance during blasting excavation is the main reason of aggravating harmful gas escape. The research results can provide a theoretical reference for preventing harmful gas from escaping in the similar engineering construction.


2021 ◽  
Author(s):  
Vibhas J. Pandey

Abstract Acid fracturing is a preferred method of stimulating low permeability limestone formations throughout the world. The treatment consists of pumping alternating cycles of viscous pad and acid to promote differential etching, thereby creating a conductive acid-etched fracture. Acid-type, pad and acid volumes, and the injection rates in the designed pump schedule are based on treatment objectives, rock-types and in-situ conditions such as temperatures, in-situ stress, proximity to water-bearing layers, and others. During the acid fracturing treatment, the acid-rock interaction is often marked by signature pressure responses, that are a combined outcome of acid reaction kinetics, responses to changes in fluid viscosity and densities, fluid-frictional drop in narrow hydraulic fractures, and other such parameters. This paper focuses on interpretation of bottomhole pressures during acid fracturing treatment to separate these individual effects and determine the effectiveness of the treatment. Unlike propped fracturing treatments where most fracturing treatments result in net pressure gain, acid fracturing treatments seldom result in net pressure increase at the end of the treatment because the in-situ stresses are generally relieved during the rock-dissolution and fracture width creation process that results from acid-mineral reactions. Not only is the extent of stress relief evident from the difference in the start and the end of the treatment instantaneous shut-in pressures, the loss of stresses is also apparent during the treatment itself, especially in jobs where the treatment data is constantly monitored and evaluated in real-time. The study reveals that the changes in pressure responses with the onset of acid in the formation can be successfully used to determine the effectiveness of treatment design and can aid in carrying out informed changes during the treatment. Better understanding of these responses can also lead to more effective treatment designs for future jobs. The interpretation developed in the study can be applied to most of the acid fracturing treatments that are pumped worldwide.


2021 ◽  
Author(s):  
Takuma Kaneshima ◽  
Fuqiao Bai ◽  
Nobuo Morita

Abstract Borehole stability depends on various parameters such as rock strength, rock deformations, in-situ stress, borehole trajectory, shale swelling, pore pressure change due to osmosis, overbalance mud weight and temperature. The objective of this work is to construct analytical and numerical equations to predict borehole failure including all these parameters, and to comprehensively propose a methodology to improve the borehole stability. Analytical solutions are developed for inclined wells with respect to in-situ stress, shale swelling, pore pressure change due to osmosis, overbalance mud weight and temperature. A numerical model is developed for 3D inclined wells with orthotropic formation and layered formation. Using the analytical and the numerical models, stress state around inclined wells are evaluated. The breakout angle is predicted based on Mohr-Coulomb, Mogi, Lade and Drucker-Prager failure theories. Polar diagrams of mud weights are compared to judge the effect of each parameter and the magnitude predicted by the different failure theories. Shale swelling and pore pressure change due to osmosis are the most difficult to estimate among above-mentioned parameters. The laboratory measured swelling of cores obtained from various formations showed that the magnitude to induce breakouts caused by swelling was the largest comparing with other parameters. Therefore, when shale stability problems occur, we need to estimate the magnitude of shale swelling and osmosis due to water potential difference. Then, to overcome the shale stability problem, we evaluated the sensitivity of human controllable parameters on borehole stability. The parameters which can be controlled by drilling engineers are overbalance, type of mud, borehole temperature and borehole trajectory. If the shale swelling is small, the borehole stability is improved by the mud weight. However, from the swelling tests from the cores of Nankai-Trough, we estimated unless we used a swelling inhibitor to reduce the swelling less than 0.1%, the well was not possible to drill through. Actually, the well was abandoned due to instability after trying side track several times. Unlike previous works, this paper uses all important parameters (swelling, temperature, pore pressure, orthotropic formation, layered formation) to estimate the stresses around inclined wells with the same formation conditions for quantitative analysis. Failure analysis include Mohr, Mogi, Lade and Drucker-Prager. Finally, the polar diagrams of critical mud weight are used to judge whether we can choose well trajectory, orientation with respect to bedding planes, mud weight, shale inhibitor, and temperature to stabilize the borehole.


1983 ◽  
Vol 20 (2) ◽  
pp. 276-287 ◽  
Author(s):  
A. M. Crawford ◽  
J. W. Bray

Rock wedges in the roofs of underground excavations may be wholly or partly self-supporting due to the mobilization of shear resistance on discontinuities bounding such wedges. The extent of the mobilization, which occurs as the wedge deforms, is markedly influenced by the magnitude of the stress field tangential to the opening and the relative stiffness of the intact rock and the shear and normal stiffnesses of the discontinuities. Analytical and numerical models are described for determining the upper and lower bounds of the failure loads for two-dimensional asymmetric and symmetric rock wedges. Experiments conducted on physical models indicate that the stability of rock wedges may be substantially reduced by loosening and that the analytical solutions may overestimate the load necessary to cause failure. Keywords: rock wedge stability, in situ stresses, joint stiffness.


SPE Journal ◽  
2015 ◽  
Vol 20 (06) ◽  
pp. 1276-1286 ◽  
Author(s):  
Mojtaba P. Shahri ◽  
Trevor T. Oar ◽  
Reza Safari ◽  
Moji Karimi ◽  
Uno Mutlu

Summary Drilling depleted reservoirs often encounters a host of problems leading to increases in cost and nonproductive time. One of these problems faced by drillers is lost circulation of drilling fluids, which can lead to greater issues such as differential sticking and well-control events. Field applications show that wellbore strengthening effectively helps reduce mud-loss volume by increasing the safe mud-weight window. Wellbore-strengthening applications are usually designed on the basis of induced-fracture characteristics (i.e., fracture length, fracture width, and stress-intensity factor). In general, these fracture characteristics depend on several parameters, including in-situ stress magnitude, in-situ stress anisotropy, mechanical properties, rock texture, wellbore geometry, mud weight, wellbore trajectory, pore pressure, natural fractures, and formation anisotropy. Analytical models available in the literature oversimplify the fracture-initiation and fracture-propagation process with assumptions such as isotropic stress field, no near-wellbore stress-perturbation effects, vertical or horizontal wells only (no deviation/inclination), constant fracture length, and constant pressure within the fracture. For more-accurate predictions, different numerical methods, such as finite element and boundary element, have been used to determine fracture-width distribution. However, these calculations can be computationally costly or difficult to implement in near-real time. The aim of this study is to provide a fast-running, semianalytical work flow to accurately predict fracture-width distribution and fracture-reinitiation pressure (FRIP). The algorithm and work flow can account for near-wellbore-stress perturbations, far-field-stress anisotropy, and wellbore inclination/deviation. The semianalytical algorithm is modeled after singular integral formulation of stress field and solved by use of Gauss-Chebyshev polynomials. The proposed model is computationally efficient and accurate. The model also provides a comprehensive perspective on formation-strengthening scenarios; a tool for improved lost-circulation-materials design; and an explanation of how they are applicable during drilling operation (in particular, through depleted zones). Sensitivity analysis included in this paper quantifies the effect of different rock properties, in-situ-stress field/anisotropy, and wellbore geometry/deviation on the fracture-width distribution and FRIP. In addition, the case study presented in this paper demonstrates the applicability of the proposed work flow in the field.


Sign in / Sign up

Export Citation Format

Share Document