Search Documents

By Benoît Valley

Optimization of the mining sequence in terms of economics (maximizing net present value) often leads to multi-front mining methods generating pillars. Significant resources are tied up in these pillars, but mining them is often challenging. In order to improve our understanding of rock mass behaviour while extracting these pillars, an extensive monitoring program has been designed and implemented at Vale?s Coleman mine (Sudbury, Canada). The program focuses on existing and new technologies that have potential for monitoring deformation and rock mass property changes. It includes both active and passive methods: gravimeters, multi-point borehole extensometers, fiber optic strain meters, fixed and portable three-component seismic arrays, borehole imaging and sonic logging, and, repeated LiDAR surveys. This paper reports results from an initial project phase, when only a small amount of mining has taken place. The goal was to test and compare technologies in order to assess their sensitivity, accuracy, repeatability and suitability for underground mining conditions. Value is gained by having a broad range of monitoring devices running side by side, enabling comparisons and benchmarking.

Jun 1, 2012

By S. Makni, G. Gagnon, F. Couët, C. Rochefort

"Comminution represents a significant part of mineral processing plant operating costs. Yet, autogenous (AG) or semi-autogenous grinding mills (SAG) are difficult to control because of disturbances like feed size distribution and rock hardness. AMIRA, within project GeMIII, developed a comminution index from laboratory crusher data. The index correlated with the Bond ball mill work index and with the A×b parameter. The present work extends the use of crusher data as a grindability indicator to real-time industrial applications. The aim is to improve comminution operation using online hardness estimation. Because the operational work index is proportional to the ore work index, rock hardness variations can be estimated from the feed and product size distributions, the throughput, and the power draw of a crusher. Real-time estimation of rock hardness can help understand the dynamic behavior of AG and SAG mills. Multivariate time series were analyzed to understand the impact of hardness on a SAG operation. Results showed that the estimated hardness influences the SAG operation and its throughput. Moreover, non-linear dimension reduction was used to ease the process visualisation and search for potential control strategies. In the future, online hardness measurements should also impact the understanding of comminution circuits and improve their operation.INTRODUCTIONThe variability in hardness, and the feed size distribution, are two of the main disturbances affecting comminution circuits (Sbarbaro & del Vilar, 2010). Size distribution can be measured by imagebased granulometry systems such as WipFrag (Maerz, Palangio, & Franklin, 1996), or by 3D vision systems when the fines are difficult to discriminate by 2D segmentation (Faucher, Makni, Gagnon, Lavoie, & Roberge, 2015).There are few examples of online hardness measurements. One example was the work of Casali et al. (2001) who developed a hardness indicator. The ore grindability soft sensor used the fraction of each lithological class in the ore, as determined by image analysis, to predict the work index of the material. The model was developed using a moving window with a width of 12 months and trends were shown over many years with one data point per month."

Jan 1, 2016

By Tom Shumka

INTRODUCTION ?The Electro-Magnetic Crack Detection (EMCD) technique is an electromagnetic inspection method used for detecting and sizing surface breaking flaws. ?It works well on structural material or welds, ferrous or non-ferrous metals. ?EMCD?s rapid scanning is achieved using a hand-held probe. ?EMCD provides reliable crack detection and sizing (length and depth). Cracks can be sized accurately. ?EMCD dramatically reduces most cleaning requirements as there is no need to clean to bare metal. ?EMCD is capable of inspection through thin metallic coatings, or through non-conducting coatings up to 10mm thick. ?EMCD allows full data storage for back-up, off-line view and audit purposes.

Oct 1, 2010

By Farid Mirahadi

Forecasting productivity of construction operations is a difficult but crucial task in planning construction projects. Over the past decades, many models have been developed to forecast productivity for different construction operations. Models made up of several functional relations and controlled by a specific number of control rules are more in line with human reasoning and logic. Neural-Network-Driven Fuzzy Reasoning (NNDFR) structure as one of these models shows a great performance for modeling datasets among which clear clusters are recognizable. Lack of the compatibility of conventional NNDFR with fuzzy clustering algorithms besides the insufficient attention paid to the optimization of number of clusters in this model, created a potential area for further research. The main contribution of the proposed model is to develop a modified NNDFR system to model construction data. To this end, Fuzzy C-Means (FCM) algorithm is substituted for K-means in NNDFR structure, and its parameters such as the number of clusters and weighting exponent are optimized through genetic algorithm. The proposed model is further verified through simulation of a construction operation in which several qualitative and quantitative factors are considered. Its implementation to the case study shows a considerable improvement of model performance with lower Mean Squared Error (MSE). The developed model assists researchers and practitioners in utilizing historical construction data to forecast productivity of construction operations with a high accuracy that could not be obtained by traditional techniques.

Aug 1, 2013

By Moatassem Abdallah

Buildings in the United States account for 72% of electricity consumption, 40% of energy consumption, 13% of water consumption, 39% of carbon footprint, and 30% of waste output. In order to minimize these negative environmental impacts, many public and private owners are requesting that their buildings be more sustainable and certified under the widely known programs such as leadership in energy and environmental design for existing buildings (LEED-EB). To accomplish this, buildings are increasingly integrating green building measures including energy efficient lighting, motion sensors, thermal pane glass, geothermal heat pumps, EnergyStar rated HVAC systems, photovoltaic systems, and wind turbines. This research paper presents the development of an automated decision support system (DSS) that is designed to optimize the selection of green building measures which can be used to upgrade existing buildings. The developed DSS incorporates two optimization models that are capable of (i) minimizing the total upgrade costs required to accomplish a specified LEED-EB certification level such as silver or gold; and (ii) maximizing the number of accredited LEED-EB points within a specified budget of upgrade costs. The DSS is designed to identify a set of optimal upgrade decisions that accomplishes these two optimization objectives. An application example is used to illustrate the capabilities of the DSS and to validate its result.

Aug 1, 2013

By Youngjib Ham

Understanding energy performance of existing buildings is vital to increasing their efficiency and reducing the overall energy consumptions. This entails facility managers to systematically monitor building energy performance and reliably identify and analyze potential problems. Currently, infrared thermography is widely used as a primary diagnostic tool for the detection of building performance problems. Nonetheless, applications of thermal images for building inspection are mainly restricted to manual and labor-intensive identification and qualitative assessment of heating or cooling loss. Automated identification of potential problems and reliable cost analysis of the associated energy loss can help homeowners to minimize financial risk of retrofits and maximize energy savings. To that end, this paper presents a new automated method for calculating the cost of energy loss for building diagnostics. In the proposed method, first, using a hand-held thermal camera, the auditors collect digital and thermal imagery from the buildings under inspection. Then, using a recently proposed method for Energy Performance Augmented Reality (EPAR) modeling, an actual 3D spatio-thermal model is generated and superimposed with a computational fluid dynamics (CFD)-based expected energy performance model. The resulting EPAR model is placed into the method proposed in this paper for cost analysis of the energy loss. Through a new 3D thermal mesh modeling using k-d tree structure and nearest neighborhood searching, performance deviations between these models are automatically calculated. Using a temperature threshold, the areas associated with potential performance problems are detected in the EPAR model and are visualized using a metaphor based on traffic light colors. Then, the actual R-values of the detected areas are measured at the level of 3D points. Based on the measured R-values and the estimated air change rate for the detected air leaks, (1) the heat loss or gain caused either by poor insulation or air infiltration/exfiltration and (2) the associated energy costs are automatically calculated. The proposed method is validated on several locations of existing residential buildings. The preliminary results show the potential of the proposed method for minimizing the inspection time as well as the risk associated with the cost analysis of retrofitting potential building performance problems.

Aug 1, 2013

By I. Clauzier

"A new method of managing electricity to reduce costs andimprove profits is discussed. Before taking steps to better manage electricity, it is essential to understand the terms peak demand and energy and their effect on an electrical bill. Knowing how electricity is presently used in the plant is the important first step. This goes beyond knowing how much electricity is used at each time of day, the individual contribution of every plant process to the total plant load must be known. At this point, each process is ranked according to importance to production in order to identify possible savings opportunities. These improvements can then be implemented by simply improving operations or in some cases, through the use of an automatic control system. IntroductionIn today's competitive environment, controlling costs is becoming increasingly important. All costs associated with operation must be controlled in an effort to improve competitive position. Electricity is a cost component that can and should be managed.Although the term management can mean different things to different people, the intent is always the same: to make the best possible use of existing resources, whether labour, equipment or electricity.The simple approach to managing electricity presented in this paper consists of four steps.I. Understand Peak Demand and Energy. It is essential to understand the terms, peak demand and energy and their effect on your electrical bill.2. Know Your Present Electrical Use. The present use of electricity in the plant must be known. This goes beyond simply knowing how much is used at each time of day, to accurately knowing how each process is contributing to the total plant electrical load, and ultimately to the electrical bill.3. Identification of Savings Opportunities. With this information, energy savings opportunities for reducing the electrical bill and improving profits can be identified.4. Implementing Improvements. Having identified areas of potential savings, improvements can then be implemented in order of their significance. This could be a matter of simply improving operations, or in some eases may mean the use of an automatic control system.These four steps can result in significant dollar savings."

Jan 1, 1987

By P. Müller, W. J. Przybylowicz, V. M. Prozesky, D. de Bruin

"Abstract -The common style of gold mineralization in the Basal Reef and other auriferous conglomerates of the Witwatersrand Supergroup comprises clusters of gold grains occurring interstitially to quartz and pyrite pebbles as well as gold accumulations attached to the surface of and infilling hairline cracks and fractures in round, compact, detrital pyrite. The latter textural association is generally interpreted as a result of remobilization of originally detrital gold during regional metamorphism and related hydrothermal fluid circulation. There is general agreement that, although gold solubility was Iow, ore modification occurred during the metahydrothermal event, but uncertainty exists about the magnitude and scale of gold remobilization .Micro-textural and chemical investigation of compact round pyrite and associated gold accumulation has established that gold is preferentially deposited on chemically heterogeneous pyrite surfaces. Combinations of alternating high- and low-As zones in pyrite result in mixed np-type conductivity accompanied by electron exchange and the reduction-driven deposition of gold at the cathode. This mode] does not explain gold accumulations which are not immediately linked to chemical heterogeneities and it is inferred that an additional control on the deposition of gold was exerted by physical defects which lead to zones of increased charge density and/or areas of focussed conductivity. The inferred precipitation mechanisms suggest that submicrocopic physical and chemical controls are the important factors that determined the sites of gold nucleation and that precipitation mechanisms based on physicochemical changes of the hydrothermal system were not involved. The Basal Reef contains on average close to 5% of round, detrital, chemically zoned pyrite which provides a large number of electrochemically active pyrite surfaces. The extremely effective electrochemical precipitation mechanism together with the low solubility of gold in Witwatersrand fluids is the main reason for the spatially restricted remobilization of gold in the conglomerate reefs."

Jan 1, 1994

By C. R. Edwards, D. R. G. Watt

"Equity Silver Mines Ltd. a member of the Placer Development group of companies operates Canada's largest silver mine near Houston in north-central British Columbia. Initially, the silver-gold-copper flotation concentrate contained excessive levels of antimony and arsenic, so a leach process was installed to remove undesirable metals and provide improved smelter terms for the leached concentrate. The leach plant which has been fully described elsewhere (1, 2, 3), was forced by deteriorating process economics to shut down in April 1984. Equipment was mothballed and alternate uses sought.The Endako Mines Division of Placer Development Limited, located some 150 kilometres east of the Equity mine site, shut down its molybdenite mining and milling facilities in the spring of 1982. Endako has continued to produce molybdic oxide (molybdenum trioxide or MTO) via intermittent roaster runs, and also a high purity, lubricant grade molybdenum disulfide in the Ultrapure Plant. With the steel industry in the doldrums, molybdenum chemicals present an alternate and growing potential for molybdenum sales. In 1982 Placer Development began testing this market. Molybdenum chemicals were produced in the United States on a toll basis using Endako MTO as process feed. The transfer of this production to Canada was being considered when Equity's leach plant shut down, and it was noted that most of the equipment and skills necessary could be provided by the leach plant. See Figure 1.A careful investigation of process requirements was completed by engineering and production staff at Equity and Endako. A provisional capital and operating cost estimate was made, and it was decided to perform a test run during July 1984. There were strong reasons to proceed with this project. A positive cash flow was expected for both Endako and Equity. Endako would have product to continue to expand in the molybdenum chemicals market with the added selling point that now the complete production process would be under direct control of Placer Development companies. Equity would further benefit since a portion of the leach plant would continue in operation thus maintaining and exercising our hydrometallurgical engineering and operating skills. Finally, Equity would have a base from which to search out other processes adaptable to the idled leach plant equipment."

Jan 1, 1986

By A. Pranowo, G. R. F. Alvear F., N. Aslin, M. S. Moats

"The toxicity of arsenic is well known and documented; however, the presence of arsenic in copper electrorefining anodes and electrolytes is critically necessary to produce high-quality cathode. Arsenic as well as antimony and bismuth concentrations vary relative to their concentration in the anode copper received from the smelter. The behaviour and benefits of arsenic in copper electrorefining and the detrimental effects of antimony and bismuth on cathode quality and tankhouse performance are discussed. This will include the minimization and mitigation of problems associated with arsenic, antimony, and bismuth—including cathode contamination, floating slimes, scaling, and anode passivation.RÉSUMÉ La toxicité de l’arsenic est bien connue et documentée; toutefois, la présence d’arsenic dans les anodes de l’électroaffinage du cuivre et dans les électrolytes est d’une nécessité critique pour produire une cathode de grande qualité. Les concentrations en arsenic, ainsi qu’en antimoine et en bismuth varient selon leur concentration dans l’anode de cuivre reçue de la fonderie. Le comportement et les avantages de l’arsenic dans l’électroaffinage du cuivre et les effets indésirables de l’antimoine et du bismuth sur la qualité de la cathode et le rendement des cuves d’électrolyse sont discutés. On traite aussi de la minimisation et de l’atténuation des problèmes associés à l’arsenic, l’antimoine et le bismuth – incluant la contamination de la cathode, les schlamms flottants, l’entartage et la passivation de l’anode."

Jan 1, 2016

By Bruce Macdonald

Operational monitoring of slope stability can be achieved by space based InSAR monitoring platforms. Natural and man-made rock slopes are monitored using a novel stack-based InSAR method that exploits both persistent point scatterers (PS) as well as coherent distributed scatterers (DS) such as areas of pavement and sparsely vegetated surfaces. The new method, Homogenous Distributed Scatterer (HDS)-InSAR identifies InSAR targets by using adaptive multilooking of statistically homogenous pixel neighbourhoods combined with thresholding the differential temporal coherence of each target candidate with respect to its surroundings. HDS-InSAR preserves PS and therefore forms a direct generalization to our proven DualScale-PSI method. The density of HDS exceeds that of PS by more than an order of magnitude often providing essentially continuous coverage over urban and industrial landscapes featuring a mix of infrastructure with mostly bare soils or rock surfaces. The improved coverage enables better definition of surface motion and underlying geology, which will lead to superior interpretation of the geophysical mechanisms causing the surface deformation, e.g. through improved integration with 3D deformation modeling. We illustrate the performance of our new HDS-InSAR method compared to DualScale-PSI using displacement results from Vancouver, Canada; Bingham mine, USA; and the Palabora mine, South Africa.

Jun 1, 2012

The original JKMRC AG/SAG mill model was developed over 30 years ago and has been upgraded a number of times since then, with the current published version being the Variable Rates model that is now over 15 years old. Since this model was developed further enhancements have been derived, some of which are applied in-house and some which have yet to be incorporated. The published and the commercially available model in JKSimMet is currently undergoing a comprehensive overhaul with a vastly extended data base for model evolution and validation. Improvements developed within the AMIRA P9 project, in-house development and a number of collaborative projects, are being incorporated into the upgraded model structure. Enhanced capabilities will include: power-based breakage rates scaling; a far wider mill filling range (from 10% to 40%) that exhibits peak in throughput as mill filling increases; higher ball loads; improved pebble discharge function; improved slurry discharge function incorporating important mill design and operational influences; and the new size dependent breakage function. The model will also be able to accommodate some influence of liner design and multi-component ores (blends of hard and soft components). This upgrade is still a work in progress but substantial changes have been made to date. This paper presents some of the new improvements and their influence on the predictive capability of the model.

Jan 1, 2012

By E. Burdukova

Long-chain primary amines are utilized as collectors in potash ore flotation. Krafft temperature of long-chain amines in brine is about 100 °C, and this makes such amines practically insoluble in brine at room temperature. When dodecylamine, melted at 70 °C and dispersed in hydrochloric acid (or acetic acid) aqueous solution, is added to a KCl-NaCl saturated brine at room temperature, the amine precipitates. The experimental results demonstrate that the precipitating amine hydrochloride particles are selectively abstracted by KCl particles, but not by NaCl particles.

Jan 1, 2009

By Andrew Bamber

A fast and objective spreadsheet-based model for the financial evaluation of innovative mining and processing strategies such as pre-concentration vs. a conventional base case is presented. The size and grade of the deposit must be chosen, from which a mining rate is selected. Expected operational data for the deposit as well as any mill testwork and results from the pre-concentration testwork are entered as input variables on the sheet, and cost and revenue impacts are calculated in the ?Output? sheet. A capital sum and operating costs for the initial mine and mill development as well as the pre-concentration facility are calculated. The net changes in capital and operating costs are calculated, and impacts on revenue through changes in metal recovery are also calculated. These outputs are used to calculate the overall impact on the cutoff grade and thus the size of the mineral reserve using a generic grade-tonnage model. The NPV and NPV profile of the two scenarios is calculated, with variations in basic parameters such as operating cost, capital requirements, metal price or grade as well as discount rate. Cost and revenue impacts are evaluated and compared to the baseline to establish a positive or negative impact overall. Different scenarios of grade, tonnage, mining rate, process performance, capital and operating cost structure can be explored using the model. Outcomes projected by the model for case studies conducted thus far have been calibrated against internal estimates by the sponsoring organizations and found to be accurate to within 30%, confirming the utility of the model in practice.

Oct 1, 2009

By Clemens Kuhnert

"Nickelhütte Aue GmbH is specialized in the recycling of waste materials containing non-ferrous metals. We treat nickel, copper, cobalt, vanadium, molybdenum and precious metal containing materials like spent catalysts, ashes, dusts, grindings, liquids, slurries and sludge's as long as they contain one or a combination of the above mentioned metals. In our pyrometallurgical works we produce non-ferrous metal concentrates suitable for further treatment in our hydrometallurgical plant where we refine Ni-, Cu-, Co- and V-chemicals. Our production is based on synergy effects due to the combination of pyro- and hydro-metallurgy. As an example we can mention our heat recovery system and internal recycling of by-products. Furthermore, the use of alternative fuels instead of fossil fuels contributes to the economic and environ-mentally friendly treatment system.INTRODUCTION Nickelhütte Aue GmbH (NHA) can look back on nearly 400 years of tradition as a pyrometallurgical works and recycling company. The existence since its foundation in 1635 as a blue dye works shows responsible and successful action in the past. Therefore, we see highest customer service, minimization of interferences with nature by closing material cycles and maximized health and occupational safety protection by avoiding injuries and diseases, still as our future duties. These duties are evaluated and monitored by an integrated management system that meets the standards of DIN EN ISO 9001, 14001, 50001 and BS OHSAS 18001. We commit ourselves to the observance of all juridical and customized demands. The management and the department managers perceive their model function. We require responsible action from all employees. To protect the continuity of our enterprise, we will also perform in future continuous improvement of products and their quality, as well as steady progress in the areas of environment protection, health protection and industrial safety. Equally the effectiveness and propriety of the management system is checked and developed regularly. In order to obtain the objectives resulting from this policy, our company determines activities and measures by involving all employees."

Jan 1, 2017

By J. E. Gill

Discovered in 1952, the deposits have been explored by diamond drilling on and near a group of islands in lake Nipissing and from under-ground workings. The ore mineral is a uranian pyrochlore. It occurs in breccia zones in association with calcite, apatite, a?cmite, feldsiiars, biotite, magnetite, sulphides, and fluorite in varying amounts. The pyrochlore concentrations occur in a steeply-dipping gneiss series with some crystalline limestone. Trends indicate a ring structure 114 miles in diameter. The main ore zone shows evidence of repeated brecciation and a succession of replacements .giving rise to a complex of rock types. Pegmatite and lamprophyre dykes cut the ore as do calcite-fluorite veins and later faults with small displacements. Columbium and uranium occur al-most exclusively in pyrochlore, but the proportions vary. There is no tantalum and only small amounts of titanium.

Jan 1, 1957

By Renato Macciotta

The argillaceous formation known as the Opalinus Clay is host to the Mont Terri Underground Research Laboratory. The URL was developed to investigate the hydrogeological, geochemical and geotechnical properties of the formation. Underground excavations at the site are characterized by non-elastic and time-dependent deformations. Excavations trending perpendicular to the strike of the bedding presented relatively minor deformations, however, tunnels parallel to the strike of the bedding showed much larger non-elastic deformations. This paper presents an analysis of the measured response of a tunnel trending parallel to the strike of the bedding in Opalinus Clay. The measured convergence is compared to tunnel experience in weak rocks to assess the tunnel stability.

Jun 1, 2012

ABSTRACT Based on mineralogical characteristics and processing techniques required, gold ores can be classified into 12 types: gold placers, quartz vein-lode ores, oxidized ores, silver-rich ores, iron oxide copper-gold ores, porphyry copper-gold ores, iron sulphide ores, arsenic sulphide ores, antimony sulphide ores, bismuth sulphide ores, telluride ores, and carbonaceous sulphide ores. The bulk mineralogy and gold deportment of each ore type differs and may have a significant impact on ore processing and metal extraction. This paper summarizes the bulk mineralogy and gold deportment in some major types of gold ore and discusses the potential impact of gold deportment on gold metallurgy.

Jan 1, 2013

By Carricajo T., S. Nazari

This paper discusses a unique dynamic solution for optimizing comminution area in real time. The proposed technology is a combination of model predictive controller, optimization methods and real-time digital twin for mineral processing. For grinding area, the ultimate goal of an effective control strategy is to maximize the mineral processing, reduce energy consumption and stabilized size classification. To have a continuous and effective comminution operation, a reliable control system that minimizes the need for human intervention is necessary. The implemented advanced control strategy for comminution area consists of two layers in hierarchical structure. The lower layer objective is to stabilize the process stabilization using model based predictive controller (MPC) technique. The unique MPC used in this layer is called BrainWave™ chosen for its robustness, convergence and delay handling; furthermore, the developed platform offers stable and repeatable solution and has short implementation time. The upper layer is a supervisory control, which sends the preferred process operation points to the lower layer to drive the process towards optimal performance. The goal of this layer is to optimize the process targets, as well as making decision. To increase process performance, automatic ore characteristic adaptations are implemented in this layer. A digital twin of comminution area is built in the second layer that runs amongst the process in parallel. By communicating, the process data in real time, a virtual densitometer is built that is being used in our advance process control strategy. Our control strategy has been developed, tested and implemented in various areas of a copper concentration plant. The targets for implementing advance process control are usually increased throughput, decrease energy consumption and variability reduction in particle size distribution. Overall, the results showed in different mine sites are increase throughput between 2 and 5%, specific energy consumption reduction about 15% and process variability reduction up to 60%. In this paper we discuss the results of the recent implementation of this technology on Minera Los Pelambres, Minera Collahuasi and CAP Minería. Their concentrator plants have installed control systems with the most updates optimization features, for SAG mill, ball mills and thickeners.

Jan 1, 2019

By A. Osouli, B. Moradi Bajestani

"Increase in permeability of rock mass is vital for both petroleum and geothermal applications. Hydraulic fracturing is a common method to create and increase the permeability for production. Due to the complexity of fracture propagation in a naturally fractured medium, understanding the interaction between hydraulic fracture (HF) and natural fracture (NF) is of great importance to optimize the fracturing process. The possibility of interaction between HF and NF is partly a function of orientation. Various scenarios are predictable for hydraulic fracture regarding the intersection with NF such as being arrested by NF, crossing NF or crossing the fissure while injecting water into it. Different parameters such as geomechanical and geometrical properties of surrounding rock mass, the operational parameters such as injection rate and the angle between created hydraulic fracture and natural fracture (approach angle) can affect the result of hydraulic fracture and natural fracture intersection. In this paper, using a series of flow-coupled Distinct Fracture Network (DFN) simulations, hydraulic fractures were analyzed to investigate the influence of the aforementioned parameters on the hydraulic fracture propagating. Results show that the most significant parameters affecting the interaction of natural and hydraulic fracture are orientation of the intersection as well as injection rate. The interaction and fracture propagation in NF and HF intersections are also studied and discussed in detail.INTRODUCTIONThe hydraulic fracture technology has substantially changed the oil, gas and geothermal industries. With the help of this technology, extraction of fossil resources from tight shale reservoirs will be economical, feasible and sustainable. Due to presence of massive oil and gas trapped in these reservoirs, this technology has been the center of attention for the past 40 years (King, 2012). Hydraulic fracturing is the process in which high pressure water is pumped into shale unconventional reservoirs with low to ultralow permeabilities. With this process paths are created in these mediums for delivering oil and gas to horizontal wellbore and increase the surrounding rock medium permeability. Economical and sustainable stimulation of these tight reservoirs is highly dependable on our understanding in regard to fracture propagation in the highly complex rock medium. Hydraulic fracturing process in rock mediums with preexisting natural fractures is more complex because of the high possibility of Hydraulic Fracture (HF) being arrested or deviated by Natural Fracture (NF). Also, HF and NF intersection can influence the required pumping pressure, injection rate, the outcome of horizontal wellbore and the affected geological underground zones. Since the 1950s, numerous studies have been conducted to predict and plan fracture propagation with the help of modeling and simulation methods. Having a thorough and accurate prediction of fracturing process will not only prevent stimulation of unwanted geological zones or pollution of aquifers, but also will improve the efficiency of the fracturing process. Previously, several studies have been conducted in order to investigate and analyze the interaction of hydraulic fracture and natural fracture (Cottrell, Hosseinpour, & Dershowitz, 2013; Doe, Lacazette, Dershowitz, & Knitter, 2013; Shakib, 2013; Zhang & Jeffrey, 2006)."

Jan 1, 2015