Effects of Mud Properties, Hole size, Drill String Tripping Speed and Configurations on Swab and Surge Pressure Magnitude during Drilling Operations

Similar documents
Determination of Pressure Losses in Hydraulic Pipeline Systems by Considering Temperature and Pressure

u y

Lecture 13 HYDRAULIC ACTUATORS[CONTINUED]

Homogeneous and Inhomogeneous Model for Flow and Heat Transfer in Porous Materials as High Temperature Solar Air Receivers

Solved Problems. (a) (b) (c) Figure P4.1 Simple Classification Problems First we draw a line between each set of dark and light data points.

Optimisation of Pressure Loss and Flow Distribution at Pipe Bifurcation

INTRODUCING THE SHEAR-CAP MATERIAL CRITERION TO AN ICE RUBBLE LOAD MODEL

Comparison of Maximum Allowable Pump Speed in a Horizontal and Vertical Pipe of Equal Geometry at Constant Power

Chapter 1 Fundamentals

Feedback-error control

A High Precision Double Tubed Hydrostatic Leveling System for Accelerator Alignment Applications

ONE-DIMENSIONAL CALCULATIONS FOR AXIAL PULLBACK FORCE DISTRIBUTIONS IN PIPES DURING DIRECTIONAL DRILLING INSTALLATIONS

Chapter 9 Practical cycles

NUMERICAL AND THEORETICAL INVESTIGATIONS ON DETONATION- INERT CONFINEMENT INTERACTIONS

III. Flow Around Bends: Meander Evolution

Chapter 10: Flow Flow in in Conduits Conduits Dr Ali Jawarneh

The extreme case of the anisothermal calorimeter when there is no heat exchange is the adiabatic calorimeter.

ANALYTICAL MODEL FOR THE BYPASS VALVE IN A LOOP HEAT PIPE

Main Menu. Summary (1)

Wolfgang POESSNECKER and Ulrich GROSS*

A SIMPLE PLASTICITY MODEL FOR PREDICTING TRANSVERSE COMPOSITE RESPONSE AND FAILURE

Analyses of Field Measured Data With Rheology and Hydraulics Models

TRIBOLOGICAL BEHAVIOR OF THRUST BEARING USING STRUCTURED SURFACES PART 1 THEORETICAL BACKGROUND

Keywords: pile, liquefaction, lateral spreading, analysis ABSTRACT

Hydraulics and wellbore cleaning

1 atm = 1.01x10 Pa = 760 Torr = 14.7 lb / in

Development of self-adaptively loading for planetary roller traction-drive transmission

Pressure variation with direction. Pressure variation with location How can we calculate the total force on a submerged surface?

Computer arithmetic. Intensive Computation. Annalisa Massini 2017/2018

Synoptic Meteorology I: The Geostrophic Approximation. 30 September, 7 October 2014

Torque and drag modelling for Redhill South-1 in the Northern Perth Basin, Australia

CFD AS A DESIGN TOOL FOR FLUID POWER COMPONENTS

MODELING AND SIMULATION OF A SATELLITE PROPULSION SUBSYSTEM BY PHYSICAL AND SIGNAL FLOWS. Leonardo Leite Oliva. Marcelo Lopes de Oliveira e Souza

The Numerical Simulation of Gas Turbine Inlet-Volute Flow Field

International Journal of Mathematics Trends and Technology- Volume3 Issue4-2012

MODELING THE RELIABILITY OF C4ISR SYSTEMS HARDWARE/SOFTWARE COMPONENTS USING AN IMPROVED MARKOV MODEL

Hydraulics! 10 An Overview!

Approximating min-max k-clustering

Understanding DPMFoam/MPPICFoam

Factors Effect on the Saturation Parameter S and there Influences on the Gain Behavior of Ytterbium Doped Fiber Amplifier

Evaluating Circuit Reliability Under Probabilistic Gate-Level Fault Models

Application of Automated Ball Indentation for Property Measurement of Degraded Zr2.5Nb

The Noise Power Ratio - Theory and ADC Testing

Lower bound solutions for bearing capacity of jointed rock

Actual exergy intake to perform the same task

RANDOM WALKS AND PERCOLATION: AN ANALYSIS OF CURRENT RESEARCH ON MODELING NATURAL PROCESSES

Linear diophantine equations for discrete tomography

4. Score normalization technical details We now discuss the technical details of the score normalization method.

Analysis of Pressure Transient Response for an Injector under Hydraulic Stimulation at the Salak Geothermal Field, Indonesia

Numerical simulation of bird strike in aircraft leading edge structure using a new dynamic failure model

An-Najah National University Civil Engineering Departemnt. Fluid Mechanics. Chapter [2] Fluid Statics

Shadow Computing: An Energy-Aware Fault Tolerant Computing Model

Preliminary Uncertainty Estimation of the Pressure Distortion Coefficient of a Pressure. Balance by FEM Calculations

VIBRODIAGNOSTICS OF COMPRESSOR VALVES VIA MUSIC PSEUDO-SPECTRA

PLATE GIRDERS - I 1.0 INTRODUCTION

Temperature, current and doping dependence of non-ideality factor for pnp and npn punch-through structures

Paper C Exact Volume Balance Versus Exact Mass Balance in Compositional Reservoir Simulation

MATHEMATICAL MODELLING OF THE WIRELESS COMMUNICATION NETWORK

MODELING OF CONSTANT VOLUME COMBUSTION OF PROPELLANTS FOR ARTILLERY WEAPONS

State Estimation with ARMarkov Models

Numerical and experimental investigation on shot-peening induced deformation. Application to sheet metal forming.

KEY ISSUES IN THE ANALYSIS OF PILES IN LIQUEFYING SOILS

Statics and dynamics: some elementary concepts

Unsteady Flow of a Dusty Conducting Fluid through porous medium between Parallel Porous Plates with Temperature Dependent Viscosity and Heat Source

which is a convenient way to specify the piston s position. In the simplest case, when φ

MODELING AND SIMULATION OF REFORMER AUTO- THERMAL REACTOR IN AMMONIA UNIT

An Investigation on the Numerical Ill-conditioning of Hybrid State Estimators

Numerical Methods for Particle Tracing in Vector Fields

Investigation of Drill Bit Heave Response to Drill Rig Heave Excitation

Estimation of dynamic behavior and energy efficiency of thrust hybrid bearings with active control

Parameters Optimization and Numerical Simulation for Soft Abrasive Flow Machining

Lilian Markenzon 1, Nair Maria Maia de Abreu 2* and Luciana Lee 3

Integrating Lidar and Atmospheric Boundary Layer Measurements to Determine Fluxes and Dynamics of Particulate Emissions from an Agriculture Facility

1 University of Edinburgh, 2 British Geological Survey, 3 China University of Petroleum

Compressor Surge Control Design Using Linear Matrix Inequality Approach

Research Article Comparison of HPM and PEM for the Flow of a Non-newtonian Fluid between Heated Parallel Plates

A Simple And Efficient FEM-Implementation Of The Modified Mohr-Coulomb Criterion Clausen, Johan Christian; Damkilde, Lars

Chapter 8 Internal Forced Convection

Cuttings transport with drillstring rotation

Modelling a Partly Filled Road Tanker during an Emergency Braking

A Qualitative Event-based Approach to Multiple Fault Diagnosis in Continuous Systems using Structural Model Decomposition

MODELING OF THE DYNAMIC RESPONSE OF A FRANCIS TURBINE

Meshless Methods for Scientific Computing Final Project

SELF-SIMILAR FLOW UNDER THE ACTION OF MONOCHROMATIC RADIATION BEHIND A STRONG CYLINDRICAL SHOCK WAVE IN A NON-IDEAL GAS

Location of solutions for quasi-linear elliptic equations with general gradient dependence

Can ISO MRS and ASTM HDB Rated Materials Be Harmonized

PERFORMANCE BASED DESIGN SYSTEM FOR CONCRETE MIXTURE WITH MULTI-OPTIMIZING GENETIC ALGORITHM

Measurement of cyclone separator

Utilizing Instrumented Stand Pipe for Monitoring Drilling Fluid Dynamics for Improving Automated Drilling Operations

Churilova Maria Saint-Petersburg State Polytechnical University Department of Applied Mathematics

Pressure-sensitivity Effects on Toughness Measurements of Compact Tension Specimens for Strain-hardening Solids

Using a Computational Intelligence Hybrid Approach to Recognize the Faults of Variance Shifts for a Manufacturing Process

arxiv: v1 [physics.data-an] 26 Oct 2012

MICROPOLAR FLUID LUBRICATION OF FINITE PARTIAL ROUGH POROUS JOURNAL BEARINGS WITH SQUEEZE EFFECT

Flexible Pipes in Trenches with Stiff Clay Walls

Final report. T. Priruenrom 1), W. Sabuga 2) and T. Konczak 2) March 2013

ONE. The Earth-atmosphere system CHAPTER

Continuous, non-destructive measuring of lateral resistance of the track. Dipl.-Ing. Roel van den Bosch, Dr. Bernhard Lichtberger

EE 508 Lecture 13. Statistical Characterization of Filter Characteristics

Combining Logistic Regression with Kriging for Mapping the Risk of Occurrence of Unexploded Ordnance (UXO)

Transcription:

International Journal of Petroleum Science and Technology ISSN 0973-638 Volume 6, Number (01),. 143-15 Research India Publications htt://www.riublication.com/ijst.htm Effects of Mud Proerties, Hole size, Drill String Triing Seed and Configurations on Swab and Surge Pressure Magnitude during Drilling Oerations *Uduak Mme and Pål Skalle Deartment of Petroleum Engineering and Alied Geohysics, NTNU, Trondheim, Norway *Corresonding Author E-mail: uduak.mme@ntnu.no Abstract Factors which affect swab and surge ressures must be recisely selected in order to control kicks or blowouts and to revent loss of circulation. In this aer a sensitivity analysis is carried out to investigate the effect of drilling fluid density, lastic viscosity, yield oint, bore hole diameter and drillcollarsto-drillie length ratio on swab and surge ressures. The analysis is based on: i. The theory of hydrostatic drilling fluid mechanics, taking the effect of mud clinging constant into consideration. ii. iii. Develoing a Matlab code to comute and simlify the flow equations Transforming the outut of this code into grahical form from which effects of mud roerties, drillstring triing seed, and hole and drill string configuration; on swab and surge ressures could be redicted. The develoed code requires fairly simle inut data which can be obtained from the rig site mud engineer. Imortant conclusions and industrial imlications of the findings are 1. Mud weight, rheology and drill collars length and size are redetermined factors which have been chosen carefully and maintained at the desired level during drilling and triing oerations in order to avoid blowout or loss of circulation roblems. These can be varied only to a limited degree.. Safe triing seeds are greatly influenced by hole size. 3. The develoed comuter rogramme can be used to redict safe drill string triing.

144 Uduak Mme and Pål Skalle 4. Triing seed is the only maniulative arameter with resect to controlling surge and swab ressure. This has been built into a grahical view for certain critical circumstances. 5. Whenever a critical situation demands that triing-out should be suorted by simultaneous uming, a grah has been develoed for this urose. Keywords: Mud roerties, Kick, Loss circulation, Mud clinging constant, Swab and Surge ressure Introduction When triing, the drillstring acts as a large iston moving through the borehole. This movement creates ressures due to friction losses between the moving string and mud. Swab ressures are created when the string moves uward and surge ressures when the string moves downward. If swab ressures are large enough, they can reduce the mud hydrostatic ressure below the formation ressure and cause a kick. Large surge ressures can fracture the formation and result in loss of circulation. Also, during the drilling oeration, it is often necessary to remove the drill string in order to change the drill bit. It may also be necessary to remove the drill string for other reasons, such as to erform desired logging oerations. After the bit is changed or after such logging oerations are comleted, the drill string must return to the borehole. The rocess of ulling out and running in the drill string is known as triing oerations. It has long been observed in the field that circulation is frequent lost after making a tri. Similarly, blowouts often occurred during the ie removal art of a tri due to formation gas (tri gas) entering the hole and lightening the mud column. Exeriment (Goins etal, 1951) and analytical (Cardwell, 1953) work showed these roblems to be largely due to down hole ressure variations caused by the iston-cylinder action of the ie and borehole. Swab and surge ressures have been studied by many researchers (Marken, 199, Mutomino, 1995; Zhong, 1995). Hussain and Sharif (009) indicated the reduction of surge ressure with the increase in eccentricity. For a artially blocked eccentric annulus with cuttings bed, the surge ressure decreases with the increase in the bed thickness. Ahmed et al. (010) resented hydraulic model to redict ressure losses while drilling and circulating drilling fluid. Hydraulic analysis of annular flow with axial motion of the inner ie have been carried out (Haige and Xisheng 1996, Fili and David 003) for different ie/borehole configurations and fluid models. Unfortunately, most of these studies have develoed models which lack in simlicity and require comlex inut arameters. In our work, effects of mud roerties, drillstring triing seed, and hole and drill string configuration; on swab and surge ressures were investigated by develoing a comuter code and transforming this code into a grahical technique. The develoed code require only simle data which could be gotten from the rig site mud engineer.

Effects of Mud Proerties, Hole size, Drill String Triing Seed 145 Mud Flow Physics The flow behaviour of drilling fluids is comlicated by the variation of aarent viscosity with rate of shear or flow. Consequently, the Newtonian fluid equations are altered for alication to tyical drilling mud systems. Swab and surge ressures can be comuted using a calculation stes based on basic mud flow Equations 1 to 10. The calculation rocedure is based on the theory of hydrodynamic viscous drag resented by Maidla and Wojtanowicz (1987) for Bingham-lastic fluids in wellbores. Hydrodynamic viscous drag is defined as the friction force between the ie string and the drilling fluid, which resist against ie movement. It deends on drilling fluid roerties, triing velocity, flow regime, ie outer diameter and wellbore diameter. The more viscous drilling fluid results the more viscous drag force.it also shows effect of clearance between string and the wellbore. The calculation rocedures includes: Ste 1 Calculation of the average effective annular velocity, V ae around the drill collars and around the drillie Average effective annular velocity is the mud velocity which roduces the viscous drag comonent of surge or swab ressure and its frame of references is the wellbore wall. The value of effective mud velocity due to moving ie wall is related to ie velocity, V by mud clinging constant C c. The value of C c deends uon the ratio of ie and hole diameters. δ + C 1 δ Vae = V c 1 Where δ reresents the ratio of ie diameter to borehole diameter δ δ lnδ 1 1 δ C c = = +, for la min ar flow (1 δ )lnδ lnδ 1 δ And 4 δ + δ δ C 1 c = + δ, for turbulent flow 3 1 δ Ste Is the mud flow laminar or turbulent? The flow critical velocity, V c around the drillcollars and around the drillie is calculated as follows: V c 1.08μ + 1.08 = μ + 9.3ρ( d ρ( d h d) h d) Y 4 Ste 3 If the flow is laminar, the ressure dro around the drillcollars and around the

146 Uduak Mme and Pål Skalle drillie is calculated as follows: L μ Vae ΔP = Y + 300( dh d) 5( dh d) 5 L is the length of drillie or drillcollars as the case may be. Ste 4 If the flow is turbulent, the ressure dro around the drillcollars and around the drillie is calculated as follows: flρvae ΔP = 6 5.8( d d) h Ste 5 The friction factor, f can be calculated using the following equations: Re 96.4 ρ V ( d d) ae h = 7 μ f = e c c1 + Re + c 3 ln(re) Where, c = 3.5378591164 c c 1 3 = 300.66099 = 0.16153971 Ste 6 Surge and swab ressure can be calculated from Equation 9 and Equation 10 resectively. Psurge = ΔPa + ΔPac + 0.05ρh 9 P swab = 0.05ρh ΔP ΔP 10 a ac Where ΔPa is ressure dro around drill ie and ΔPac is ressure dro around drillcollars. Formation fracturing ressure can be estimated using the following equation (Hubbert, 1957): ν Pf = σ ov + σ t 11 1 ν Normally, the tensile strength of reservoir rocks is neglected as a worst case could be encountered (equal zero). The overburden stress for normally stressed formations 8

Effects of Mud Proerties, Hole size, Drill String Triing Seed 147 can be evaluated as follows: Psi σ ov = 1 h 1 ft Formation fracture gradient for dee wells can be estimated using the following equation: ν G f = 13 1 ν Assuming that average Poisson s ratio, ν for most reservoir rocks equal to 0.5, the formation fracturing ressure gradient (comuted using Equation 13) equals to 0.67 si/ft. Formation ore ressure gradient normally constant and assumed to be 0.45 si/ft in this study. Figure 1: Flow chart for the develoed Matlab code

148 Uduak Mme and Pål Skalle Figure 1 indicates the flow chart for the Matlab code that has been develoed to do all the above mentioned calculations. Based on the above rocedure a sensitivity analysis was erformed to investigate the effect of fluctuation of the model arameters on surge and swab ressures including: i. Mud roerties (mud weight, lastic viscosity and yield oint) ii. Hole diameter iii. Drillcollars-to-drillie length ratio iv. Formation ore fluid ressure v. Formation fracturing ressure Table1: Matlab Code Inut Data Parameter Value or Range Unit Mud weight, ρ 8.5-1.5 g Mud lastic viscosity, μ 15-55 c Mud yield oint, Y 5-00 Ib/100sq.ft Well deth, h 6000-15000 ft Hole size, d h 7.875-9.875 inches Drillcollars-to-drillie ratio 0.0345-0.1538 fraction Drill ie size (ID,OD) 3.86, 4.5 inches Drillcollars size (ID,OD).813, 6.75 inches Drill string triing seed 0-18 ft/sec Formation ore fluid ressure gradient 0.45 Psi/ft Formation fracturing ressure gradient 0.67 Psi/ft Poisson s ratio 0.5 fraction By using the inut data resented in Table 1 above and the develoed comuter code, a sensitivity analysis was conducted. Results and Discussions Based on the results of the sensitivity analysis, Figures to 6 were lotted. Critical ie running seeds for both surge and swab cases were clarified in the Figures. Drilling fluid roerties including density, yield oint, and viscosity are imortant factors which greatly affect the magnitude of surge and swab ressures generated during tri in and tri out oeration resectively. Figure 3 shows the effect of mud weight on critical ie triing seeds range.

Effects of Mud Proerties, Hole size, Drill String Triing Seed 149 Figure : The effect of mud weight on surge and swab ressure magnitude The critical ie triing seed is the seed beyond which loss of circulation or blowout could haen. When the mud density increases the range of safe ie tri in seed decreases This effect is attributed to the large increase in mud ressure (in addition to the iston-cylinder action caused by the drill string) oosing the formation being drilled and the results will be serious fractures leading to loss of circulation roblem. On the other hand, the increase in mud density can easily control formation ore fluid ressure, therefore, wide range of safe ie tri-out seeds can be alied as shown in Figure. Safe ie triing seeds are affected by mud rheological roerties such as mud lastic viscosity and yield oint. Figure 3: The effect of mud lastic viscosity on surge and swab ressure magnitude

150 Uduak Mme and Pål Skalle This effect could be due to the fact that viscous mud magnify the iston-cylinder action roduced by drillstring u and down movements as shown in Figures 3 and 4. Figure 4: The effect of mud yield oint on surge and swab ressure magnitude Drillcollars has a diameter bigger than the normal drillies. Therefore it adds and extra ressure on the formation during triing in and drain the formation ore fluid into the wellbore during triing out due to the iston-cylinder action as shown in Figures 5. Figure 5: The effect of drillcollars to drillie length ratio on surge and swab ressure magnitude

Effects of Mud Proerties, Hole size, Drill String Triing Seed 151 Hole size have a large effect on safe drillstring triing seeds. This is because during triing out, large hole diameter allows the drilling fluid to raidly fill in the lace which was occuied by the drillstring, therefore the formation ore fluid ressure is easily controlled and avoid blowouts. During ie tri-in, large hole diameter rovides bigger assage area for the mud, therefore, the iston-cylinder action acting oosite to the formation is minimized and the formation fracturing can be avoided as shown in Figure 6 (Hole size). Figure 6: The effect of hole size on surge and swab ressure magnitude Thus, by lotting data obtained from the develoed comuter rogramme, effects of mud roerties, hole size, drill string triing seed and configurations on swab and surge ressure magnitude can be redicted during drilling oerations Conclusions Based on the erformed analysis, the following conclusions can be drawn: Mud weight, rheology and drill collars length and size are redetermined factors which have been chosen carefully and maintained at the desired level during drilling and triing oerations in order to avoid blowout or loss of circulation roblems. These can be varied only to a limited degree. Safe triing seeds are greatly influenced by hole size. The develoed comuter rogramme can be used to redict safe drill string triing. Triing seed is the only maniulative arameter with resect to controlling surge and swab ressure. This has been built into a grahical view for certain critical circumstances. Whenever a critical situation demands that triing out should be suorted by simultaneous uming, a grah has been develoed for this urose.

15 Uduak Mme and Pål Skalle References [1] Ahmed, R., Enfins, Miftah-El-Kheir, H., Laget, M., and Saasen, A., (010): The Effect of Drillstring Rotation on Equivalent Circulation Density: Modeling and Analysis of Field Measurements, aer SPE 135587 resented at the SPE Annual technical Conference and Exhibition held in Florence, Italy, 19- Setember. [] Cardwell, W. T., Jr.,(1953): Pressure Changes Caused by Pie Movement, API Drilling and Production Practices,.97. [3] Fili, P. and David, J. (003): Axial Couette-Poiseuille flow of Power Law viscolastic Fluids in concentric Annuli, Journal of Petroleum Science and Engineering, Vol. 40,. 111-119. [4] Goins, W. C., Weichert, J. P., Burba, J. L., Dawson, D.D., and Telitz, A.J., (1951): Down the Hole Pressure Surges and Their Effect on Loss of Circulation, API Drilling and Production Practices,.15. [5] Haige, W. and Xisheng, L. (1996): study on surge Pressure for Yield- Pseudolastic Fluid in a Concentric Annulus, Alied mathematics and mechanics, Vol.17, No. 1, 15-3. [6] Hubbert, M and Willis, D. (1957): Mechanics of Hydraulic Fracturing. [7] Hussain, Q. E. and Sharif M. A. R. (009): Viscolastic fluid flow in irregular eccentric annuli due to axial motion of the inner ie, The Canadian journal of Chemical Engineering, vol. 75, issue 6,. 1038-1045. [8] Maidla, E.E., Wojtanowicz, A.K. (1987): Field Method of Assessing Borehole Friction for Directional Well Casing, Society of Petroleum Engineers, SPE Middle East Oil Show, Manama, Bahrain, March 7-10, 15696. [9] Marken, C. D., Xiaojun, H. and Arild, S, (199): The Influence of Drilling Conditions on Annular Pressure Losses, SPE aer no. 4598 resented at the 67 th Annual Technical Conference and Exhibition of the Society of Petroleum Engineers held in Washington, Dc., October 4-7. [10] Mutomino, M., Rubertone, M. and Cesare, P. (1995): An Innovative Models for drilling Fluid Hydraulics, SPE aer no. 959, Asia and Pacific Oil & Gas Conference held in Kuala Lumur, Malaysia, 0- March. [11] Zhong, B., Zhou, K. and Yuan, Q. (1995): Equations hels Calculate Surge and Swab ressures in Inclined Wells, Oil & Gas Journal, Setember 18, 74-77.