Formation damage is a common problem in petroleum reservoirs and happens in different stages of reservoir development from drilling to production. The causes of formation damage include particle invasion, formation fines migration, chemical precipitation, and pore deformation or collapse. Formation damage adversely affects productivity of wells by reducing the permeability of near wellbore region. Furthermore, formation damage also affects well logging results. Therefore, understanding the mechanism of formation damage is vital to predict the extent and severity of formation damage and to control it. This thesis is focused on the study of formation damage caused by external particle invasion. A simplified numerical method based on a commercial code PFC (Particle Flow Code) is proposed to simulate the particle invasion process. The fluid-particle interaction is simplified as hydrodynamic drag forces acted on particles by fluids; the particle-grain interaction is modeled as two rigid balls on contact. Furthermore, an pore network flow model is developed in this study to estimate permeability of damaged formations, which contain two well-separated particle sizes. The effects of the particle size and the initial formation porosity on formation damage are studied in detail. Our study shows that big particles tend to occupy the formation face, while small particles invade deep into the formation. Moreover, particles which are smaller than pore throats (entrances) impair permeability more than those bigger than pore throats. Our study also indicates that a higher initial formation porosity results in more particle invasion and permeability impairment. It is suggested that, in order to reduce formation damage, mud particle size distributions should be carefully selected according to given formation properties. Although our model has some limitations, it may serve as a tool to predict formation damage according to given parameters, and to understand the mechanism of formation damage from a micro-scopic point of view. / text

The failure of brittle rocks around deep underground excavations due to the high induced stress is controlled by the crack accumulation in the rock. The study shows that the damage initiation strength, CI, corresponds to the long-term strength, and the short-term strength of the brittle rocks in-situ is the crack interaction strength, CD. Therefore the damage thresholds that are being used for the calibration and validation of numerical models are important parameters in the design of underground structures. The accurate detection of the damage thresholds is important as they define the in-situ behaviour of the brittle rocks. The two most common methods of detecting damage thresholds are the Acoustic Emission method and the strain measurement method. Apparent discrepancy that exists between the accuracy of these methods was the author’s motivation for comparing these two methods on Stanstead and Smaland granites. The author introduced two new parameters based on the measured strains for improving the strain measurement method. Based on the comparisons, the author is of the opinion that the Acoustic Emission method is a more accurate method of detecting damage thresholds. Numerical models are an important tool in the design of underground structures. The numerical methods that are able to simulate fractures explicitly have the ability to predict the brittle failure, the density and the extension of the microcracks around the opening. Itasca’s Particle Flow Code (PFC) was used in this study due to its potential to simulate fractures explicitly. Calibration of PFC models to Unconfined Compressive Strength properties of the rock does not mean that the model will behave correctly under other confining stresses or in tension. The author has tried to solve this problem by different methods and developing new procedures. Improvements in the model behaviour have been achieved but more work is required. The definition, and detection and calibrated simulation of rock damage thresholds for calibration of numerical models is helpful for a successful design of underground excavations and long term, lower bound strength, a critical design parameter for deep geological repositories for the storage of nuclear wastes, for example. / Thesis (Master, Geological Sciences & Geological Engineering) -- Queen's University, 2010-09-23 13:59:28.795

[pt] As alterações no solo decorrentes de um elemento de fundação profunda e seus desempenhos sob a aplicação de carga axial são processos há muito tempo estudados na engenharia civil. Diversos fatores como, método de instalação utilizado, formato da estaca, interações solo-estrutura, mecanismos de transferências de carga, movimentação do solo e alterações na compressibilidade e tensões do solo adjacente, apresentam desafios importantes que ainda não foram totalmente compreendidos nos fenômenos de penetração e capacidade de suporte em estacas. Diversos avanços foram realizados ao longo das últimas décadas para se investigar estes comportamentos, a partir procedimentos experimentais e novas formas de instrumentação, assim como ferramentas numéricas sofisticadas com o emprego de complexos modelos constitutivos em elementos finitos. Apesar destes avanços, a modelagem numérica dos processos citados, com todas as suas complexidades, ainda encontra alguns desafios. Devido a facilidade em lidar com simulações de grandes deformações e de captar o comportamento dilatante e nãolinear de solos granulares, o Método dos Elementos Discretos apresenta uma excelente ferramenta para investigar estes processos, sem grandes complicações. O presente trabalho procurou avaliar os comportamentos obtidos a partir de diferentes processos de instalação da estaca e seus efeitos nos resultados da prova de carga estática em solos granulares. As alterações de tensão e deslocamento foram avaliadas nos diferentes modelos e discutindo sobre uma metodologia básica para obter correspondências qualitativas e quantitativas com os diferentes comportamentos de campo e laboratório. Para este estudo foram utilizados os programas PFC, na versão 2D, e o programa UDEC, da Itasca co. / [en] The disturbances experienced by the soil owing to the load applied to a deep foundation and its relative behavior consist of long time studied phenomena in civil engineering. Several factors such as the installation methods, the pile geometry, the interactions between soil and structure, the load-transfer mechanisms, the soil movements and the disturbances in the stress and compressibility fields present major challenges that have not yet been completely understood. Numerous advances have been observed throw-out the last decades, in order to investigate these behaviors starting from the different pile instrumentations, the use of calibration cameras and centrifuges and most recently the measurement of the stress and strain fields inside the soil mass in model tanks. Despite the advances the numerical modelling of those processes still faces major challenges. Due to simplified approach used by the Discrete Element Method to simulate large deformation and the dilant non-linear behavior of granular soils, it presents as an excellent tool to investigate these processes without further complications. The present work proposed to evaluate the different behaviors obtained with the variations of installation methods investigated as well as their effects in the results of the Pile Load Test. The disturbances were also evaluated in the different models considered and a basic method to achieve qualitative and quantitative comparisons was discussed. These studies were made possible with the help of the PFC2D and UDEC programs developed by Itasca co.

[pt] Diferentes mecanismos de ruptura são considerados no momento de avaliar a estabilidade de um maciço rochoso fraturado. Entre estes, os mecanismos de ruptura tipo planar, em cunha e tombamentos têm sido estudados intensivamente, existindo atualmente modelos matemáticos que permitem avaliá-los. Estes mecanismos de ruptura são restritos a taludes pequenos e com fraturas contínuas, nas quais o deslizamento ocorre ao longo destas descontinuidades. Em casos de taludes de grande altura ou quando a persistência das fraturas é pequena em relação à escala do talude, o fraturamento torna-se descontínuo. Neste caso, o mecanismo de ruptura mais provável é o tipo Step-Path, o qual, a superfície de ruptura é formada por fraturas que se propagam através da rocha intacta juntando-se entre elas. Este fenômeno de união de fraturas é chamado de coalescência. Análises de estabilidade, como os probabilísticos ou por equilíbrio limite, são usados atualmente para avaliar estes tipos de rupturas, não se tendo ainda o desenvolvimento de um modelo numérico que possa representá-lo e reforçar estas teorias. O presente trabalho avalia o uso do Método dos Elementos Discretos na modelagem do mecanismo de ruptura tipo step- path, realizando uma análise de estabilidade que permita comparar os seus resultados com o método de equilíbrio limite. Foi utilizado o programa PFC nas versões 2D e 3D, assim como o programa FracGen para a geração de fraturas tridimensionais. A análise tridimensional foi feita mediante um acoplamento PFC3D-FracGen. A pesquisa inclui a análise e modelagem dos fenômenos de coalescência em amostras, assim como a influência da anisotropia na resistência das rochas em ensaios triaxiais. / [en] Different failure mechanisms are considered when a fracturated rock mass is valued. Some of them are being subject of accurate study, like planar failure mechanism, wedges and toppling, which are currently valued by mathematical models. These failure mechanisms are restricted to small slopes and with continue fractures, where the sliding occurs along these discontinuities. To height slopes or when the fracture persistence is smaller than the slope scale, the fracturing becomes discontinuous. In this case, the most probable failure mechanism to happen is the step-path type, in which the failure surface is composed by fractures that propagate through the intact rock and that are joined together. This phenomenon of fracture union is known as coalescence. Stability analysis, like probability analysis or limit equilibrium analysis are currently utilized to evaluate this kind of failures, but its important to develop a numerical model to represent and reinforce these theories. This work aims to evaluate the use of Discrete Element Method to model step-path failure mechanism on a stability analysis and to compare the results with limit equilibrium method. The program used to simulate the slope is PFC (2D and 3D) and the program FracGen was used to generate three-dimensional fractures. Three-dimensional analysis was done by a coupling between PFC3D and FracGen. The research includes the analysis and modeling of coalescence phenomenon on rock samples, as well as the analysis of the anisotropy influence on rock strength obtained from triaxial tests.

This research consists of two main parts: development of a rock expert system (RES) as an easy-to-use and effective tool for evaluating rock properties, and modification and utilization of the three-dimensional Particle Flow Code (PFC3D) to analyze rock behavior. Because of different reasons, it is often difficult to obtain the rock property values directly. As an alternative, typical values and empirical correlations are often used to evaluate the rock property values. However, the typical values and empirical correlations come in various forms and are scattered in different sources. It is often difficult, time-consuming or even impossible for an engineer to find appropriate information to estimate the required rock properties. So in the first part of the research, the RES was developed as an easy-to-use and effective tool for evaluating rock properties by conducting detailed review and evaluation of well determined values and empirical correlations of rock properties in the published literature, and developing a central database and data application tools. The study of RES demonstrates the storage of rock property values and correlations is strongly applicable and the web based data application tool is effective to use and easy expandable. Considering its granular nature, the discrete element method (DEM) has been widely adopted to analyze the mechanical behavior of rock. The Particle Flow Code (PFC) is one of the most popular DEM softwares. The basic idea of PFC is to treat rock as an assembly of bonded particles that follow the law of motion and consider the model behavior dominated by the formation and interaction of micro cracks developed within the particle-particle cement (bond). Unlike the continuum methods, PFC can deal with the natural process from micro cracking to macro failure, without predefining a failure criterion for the rock. However, there are still issues related to the application of PFC to analyze different rock problems. For example, so far, most of the studies use PFC2D although many of the problems are three dimensional and should be better simulated with PFC3D. It is also found that the simulations using the default PFC parallel bond model extremely underestimate the ratio of unconfined compressive strength to tensile strength (UCS/T). So in the second part of the research, the important aspects related to the application of PFC3D, including model scale, particle size distribution and contact model, were studied, a new contact model was developed for addressing the limitation of the default PFC3D on obtaining unrealistically low UCS/T ratios, and finally the new contact model was used to investigate rock fracture initiation and propagation.

Doctor of Philosophy / Department of Industrial & Manufacturing Systems Engineering / Shuting Lei / This study experimentally and numerically investigates laser assisted milling (LAMill) of silicon nitride ceramics. Experiments are conducted to study the machinability of Si3N4 under LAMill. The effects of temperature on cutting forces, tool wear, surface integrity, edge chipping and material removal mechanisms are investigated. It is shown that when temperature increases, cutting force and tool wear are significantly decreased, surface integrity is improved, chip size is increased and material removal demonstrates more plastic characteristics. The mechanisms of edge chipping at elevated temperature are investigated theoretically and experimentally. When temperature is above the softening point and below the brittle/ductile transition temperature, the mechanism is mainly through softening. When temperature is above the brittle/ductile transition temperature, toughening mechanism contributes significantly to the reduced edge chipping. The coupled effect of softening and toughening mechanisms shows that temperature range between 1200 to 1400°C has the most significant effect to reduce edge chipping. Distinct element method (DEM) is applied to simulate the micro-mechanical behavior of Si3N4. First, quantitative relationships between particle level parameters and macro-properties of the bonded particle specimens are obtained, which builds a foundation for simulation of Si3N4. Then, extensive DEM simulations are conducted to model the material removal of machining Si3N4. The simulation results demonstrate that DEM can reproduce the conceptual material removal model summarized from experimental observations, including the initiation and propagation of cracks, chip formation process and material removal mechanisms. It is shown that material removal is mainly realized by propagation of lateral cracks in machining of silicon nitride. At the elevated temperature under laser assisted machining, lateral cracks are easier to propagate to form larger machined chips, there are fewer and smaller median cracks therefore less surface/subsurface damage, and crushing-type material removal is reduced. The material removal at elevated temperature demonstrates more plastic characteristics. The numerical results agree very well with experimental observations. It shows that DEM is a promising method to model the micro-mechanical process of machining Si3N4.

This work explains the effect of failure on permeability anisotropy and dilation in sands. Shear failure is widely observed in field operations. There is incomplete understanding of the influence of shear failure in sand formations. Shear plane orientations are dependent on the stress anisotropy and that view is confirmed in this research. The effect of shear failure on the permeability is confirmed and calculated. Description of permeability anisotropy due to shear failure has also been discussed. In this work, three-dimensional discrete element modeling is used to model the behavior of uncemented and weakly cemented sand samples. Mechanical deformation data from experiments conducted on sand samples is used to calibrate the properties of the spherical particles in the simulations. Orientation of the failure planes (due to mechanical deformation) is analyzed both in an axi-symmetric stress regime (cylindrical specimen) and a non-axi-symmetric stress regime (right cuboidal specimen). Pore network fluid flow simulations are conducted before and after mechanical deformation to observe the effect of failure and stress anisotropy on the permeability and dilation of the granular specimen. A rolling resistance strategy is applied in the simulations, incorporating the stiffness of the specimens due to particle angularity, aiding in the calibration of the simulated samples against experimental data to derive optimum granular scale elastic and friction properties. A flexible membrane algorithm is applied on the lateral boundary of the simulation samples to implement the effect of a rubber/latex jacket. The effect of particle size distribution, stress anisotropy, and confining pressure on failure, permeability and dilation is studied. Using the calibrated micro-properties, simulations are extended to non-cylindrical specimen geometries to simulate field-like anisotropic stress regimes. The shear failure plane alignment is observed to be parallel to the maximum horizontal stress plane. Pore network fluid flow simulations confirm the increase in permeability due to shear failure and show a significantly greater permeability increase in the maximum horizontal stress direction. Using the flow simulations, anisotropy in the permeability field is observed by plotting the permeability ellipsoid. Samples with a small value of inter-granular cohesion depict greater shear failure, larger permeability increase and a greater permeability anisotropy than samples with a larger value of inter-granular cohesion. This is estimated by the number of micro-cracks observed. / text

Two new three-dimensional rock mass strength criteria are developed in this dissertation by extending an existing rock mass strength criterion. These criteria incorporate the effects of the intermediate principal stress, minimum principal stress and the anisotropy resulting from these stresses acting on the fracture system. In addition, these criteria have the capability of capturing the anisotropic and scale dependent behavior of the jointed rock mass strength by incorporating the effect of fracture geometry through the fracture tensor components. Another significant feature of the new rock mass strength criterion which has the exponential functions (equation 6.7) is having only four empirical coefficients compared to the existing strength criterion which has five empirical coefficients; if the joint sets have the same isotropic mechanical behavior, the number of the empirical coefficients reduces to two in this new strength criterion (equation 6.10). The new criteria were proposed after analyzing 452 numerical modeling results of the triaxial, polyaxial and biaxial compression tests conducted on the jointed rock blocks having one or two joint sets by the PFC3D software version 5. In this research to have several samples with the same properties a synthetic rock material that is made out of a mixture of gypsum, sand and water was used. In total, 20 joint systems were chosen and joint sets have different dip angles varying from 15 to 60 at an interval of 15 with dip directions of 30 and 75 for the two joint sets. Each joint set also has 3 persistent joints with the joint spacing of 42 mm in a cubic sample of size 160 mm and the joints have the same isotropic mechanical behavior. The confining stress combination values were chosen based on the uniaxial compressive strength (UCS) value of the modeled intact synthetic rock. The minimum principal stress values were chosen as 0, 20, 40 and 60 percent of the UCS. For each minimum principal stress value, the intermediate principal stress value varies starting at the minimum principal stress value and increasing at an interval of 20 percent of the UCS until it is lower than the strength of the sample under the biaxial loading condition with the same minimum principal stress value. The new rock mass failure criteria were developed from the PFC3D modeling data. However, since the joint sets having the dip angle of 60 intersect the top and bottom boundaries of the sample simultaneously, the joint systems with at least one of the joint sets having the dip angle of 60 were removed from the database. Thus, 284 data points from 12 joint systems were used to find the best values of the empirical coefficients for the new rock mass strength criteria. λ, p and q were found to be 0.675, 3.16 and 0.6, respectively, through a conducted grid analysis with a high R2 (coefficient of determination) value of 0.94 for the new criterion given by equation 6.9 and a and b were found to be 0.404 and 0.972, respectively, through a conducted grid analysis with a high R2 value of 0.92 for the new criterion given by equation 6.10. The research results clearly illustrate how increase of the minimum and intermediate principal stresses and decrease of the joint dip angle, increase the jointed rock block strength. This dissertation also illustrates how different confining stress combinations and joint set dip angles result in different jointed rock mass failure modes such as sliding on the joints, failure through the intact rock and a combination of the intact rock and joint failures. To express the new rock mass strength failure criteria, it was necessary to determine the intact rock strengths under the same confining stress combinations mentioned earlier. Therefore, the intact rock was also modeled for all three compression tests and the intact rock strengths were found for 33 different confining stress combinations. Suitability of six major intact rock failure criteria: Mohr-Coulomb, Hoek-Brown, Modified Lade, Modified Wiebols and Cook, Mogi and Drucker-Prager in representing the intact rock strength was examined through fitting them using the aforementioned 33 PFC3D data points. Among these criteria, Modified Lade, Modified Mogi with power function and Modified Wiebols and Cook were found to be the best failure criteria producing lower Root Mean Square Error (RMSE) values of 0.272, 0.301 and 0.307, respectively. Thus, these three failure criteria are recommended for the prediction of the intact rock strength under the polyaxial stress condition. In PFC unlike the other methods, macro mechanical parameters are not directly used in the model and micro mechanical parameter values applicable between the particles should be calibrated using the macro mechanical properties. Accurate calibration is a difficult or challenging task. This dissertation emphasized the importance of studying the effects of all micro parameter values on the macro mechanical properties before one goes through calibration of the micro parameters in PFC modeling. Important effects of two micro parameters, which have received very little attention, the particle size distribution and the cov of the normal and shear strengths, on the macro properties are clearly illustrated before conducting the said calibration. The intact rock macro mechanical parameter values for the Young’s modulus, uniaxial compression strength (UCS), internal friction angle, cohesion and Poisson's ratio were found by performing 3 uniaxial tests, 3 triaxial tests and 5 Brazilian tests on a synthetic material made out of a mixture of gypsum, sand and water and the joint macro mechanical parameter values were found by conducting 4 uniaxial compression tests and 4 direct shear tests on jointed synthetic rocks with a horizontal joint. Then the micro mechanical properties of the Linear Parallel Bond Model (LPMB) and Modified Smooth Joint Contact Model (MSJCM) were calibrated to represent the intact rock and joints respectively, through the specific procedures explained in this research. The similar results obtained between the 2 polyaxial experiments tests of the intact rock and 11 polyaxial experimental tests of the jointed rock blocks having one joint set and the numerical modeling verified the calibrated micro mechanical properties and further modification of these properties was not necessary. This dissertation also proposes a modification to the Smooth Joint Contact Model (SJCM) to overcome the shortcoming of the SJCM to capture the non-linear behavior of the joint closure varying with the joint normal stress. Modified Smooth Joint Contact Model (MSJCM) uses a linear relation between the joint normal stiffness and the normal contact stress to model the non-linear relation between the joint normal deformation and the joint normal stress observed in the compression joint normal stiffness test. A good agreement obtained between the results from the experimental tests and the numerical modeling of the compression joint normal test shows the accuracy of this new model. Moreover, another shortcoming associated with the SJCM application known as the interlocking problem was solved through this research by proposing a new joint contact implementation algorithm called joint sides checking (JSC) approach. The interlocking problem occurs due to a shortcoming of the updating procedure in the PFC software related to the contact conditions of the particles that lie around the intended joint plane during high shear displacements. This problem increases the joint strength and dilation angle and creates unwanted fractures around the intended joint plane.