+

WO2018017724A1 - Procédés de commande de trajectoire de trépan par prédiction de marche de trépan et forage en spirale - Google Patents

Procédés de commande de trajectoire de trépan par prédiction de marche de trépan et forage en spirale Download PDF

Info

Publication number
WO2018017724A1
WO2018017724A1 PCT/US2017/042876 US2017042876W WO2018017724A1 WO 2018017724 A1 WO2018017724 A1 WO 2018017724A1 US 2017042876 W US2017042876 W US 2017042876W WO 2018017724 A1 WO2018017724 A1 WO 2018017724A1
Authority
WO
WIPO (PCT)
Prior art keywords
drill bit
bha
bit
model
drilling
Prior art date
Application number
PCT/US2017/042876
Other languages
English (en)
Inventor
Reed W. SPENCER
Timothy P. Uno
Jonathan M. HANSON
Original Assignee
Baker Hughes Incorporated
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Priority claimed from US15/348,003 external-priority patent/US10296678B2/en
Application filed by Baker Hughes Incorporated filed Critical Baker Hughes Incorporated
Priority to CA3031465A priority Critical patent/CA3031465C/fr
Priority to GB1901944.7A priority patent/GB2567376B/en
Priority to CN201780054409.3A priority patent/CN109661502B/zh
Priority to BR112019001148-0A priority patent/BR112019001148B1/pt
Publication of WO2018017724A1 publication Critical patent/WO2018017724A1/fr
Priority to SA519400930A priority patent/SA519400930B1/ar
Priority to NO20190197A priority patent/NO20190197A1/no

Links

Classifications

    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B44/00—Automatic control systems specially adapted for drilling operations, i.e. self-operating systems which function to carry out or modify a drilling operation without intervention of a human operator, e.g. computer-controlled drilling systems; Systems specially adapted for monitoring a plurality of drilling variables or conditions
    • E21B44/02—Automatic control of the tool feed
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B10/00—Drill bits
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B44/00—Automatic control systems specially adapted for drilling operations, i.e. self-operating systems which function to carry out or modify a drilling operation without intervention of a human operator, e.g. computer-controlled drilling systems; Systems specially adapted for monitoring a plurality of drilling variables or conditions
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B45/00—Measuring the drilling time or rate of penetration
    • E—FIXED CONSTRUCTIONS
    • E21—EARTH OR ROCK DRILLING; MINING
    • E21B—EARTH OR ROCK DRILLING; OBTAINING OIL, GAS, WATER, SOLUBLE OR MELTABLE MATERIALS OR A SLURRY OF MINERALS FROM WELLS
    • E21B7/00—Special methods or apparatus for drilling
    • E21B7/04—Directional drilling
    • E21B7/06—Deflecting the direction of boreholes
    • G—PHYSICS
    • G01—MEASURING; TESTING
    • G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
    • G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G06F30/23—Design optimisation, verification or simulation using finite element methods [FEM] or finite difference methods [FDM]

Definitions

  • Embodiments of the present disclosure relate to methods of controlling drill bit trajectory in a subterranean formation to improve drilling operations and plans for drilling by accounting for and predicting bit walk and wellbore spiraling.
  • Geologic formations are used for many applications such as hydrocarbon production, geothermal production, and carbon dioxide sequestration.
  • boreholes are drilled into the formations to provide access to them.
  • a predetermined drill path extending into the earth formation may be provided for an operator to follow. The drill path may curve, turn, or otherwise be nonlinear, requiring the operator to control a direction in which an earth- boring tool proceeds into the earth formation.
  • Components for inclusion in a drillstring such as a bottomhole assembly (BHA) and one or more earth-boring tools (e.g., an earth- boring drill bit, a reamer, or another tool configured to remove earth material when forming or enlarging a borehole), may be selected for their ability to perform within, and to create, a nonlinear borehole.
  • BHA bottomhole assembly
  • earth-boring tools e.g., an earth- boring drill bit, a reamer, or another tool configured to remove earth material when forming or enlarging a borehole
  • the drill bit When a drill bit walks in such a way that the direction angle increases, the drill bit may be said to walk to the right or to exhibit "right walk.” Similarly, when a drill bit walks in such a way that the direction angle decreases, the drill bit may be said to walk to the left or to exhibit “left walk.” When a drill bit does not walk or drift away from the straight path that is parallel to the longitudinal axis of the well borehole at the bottom thereof, the bit may be referred to as an "anti-walk” drill bit and may be said to exhibit "neutral walk.” In a similar manner, when a drill bit drifts in a direction such that the inclination angle increases, the drill bit is said to exhibit a tendency to "build,” and when a drill bit drifts in a direction such that the inclination angle decreases, the drill bit is said to exhibit a tendency to "drop.”
  • WOB weight-on-bit
  • a wellbore design may be created using three-dimensional modeling software by inputting design variables associated with drilling parameters and lithology data and using the computational software to estimate by mathematical calculations the reactive forces and torques applied to the drill bit and BHA by the surrounding subterranean formation during drilling to estimate the trajectory of the drill bit and BHA through the subterranean formation.
  • BHA bottomhole assembly
  • the method includes: constructing a BHA model of the BHA, the BHA model having (a) dimensions, geometry, mass distribution, material density, and material stiffness of the BHA and (b) dimensions and geometry of the borehole as predicted to be drilled, the BHA model being configured to calculate (c) one or more points of contact between the BHA and a wall of the borehole and (d) forces imposed on the BHA at the one or more contact points, the forces comprising a side force vector imposed on the drill bit, and (e) drill bit tilt; calculating confined compressive strength of the rock formation using an axial motion drill bit model that receives drilling parameters for the drill rig, drill bit design information, and lithology information comprising formation rock strength; calculating lateral motion of the drill bit using a lateral motion drill bit model that receives (i) the calculated confined compressive strength, (ii) the lithology information, (iii) the drill bit design information, and (iv) the drill bit side force vector and drill bit tilt from the BHA model; calculating
  • a non-transitory computer-readable medium having computer-executable instructions for predicting a path of a borehole that will be drilled in an earth formation by a bottomhole assembly (BHA) having a drill bit coupled to a drill tubular, the BHA being operated by a drill rig, by implementing the following steps.
  • the steps include: constructing a BHA model of the BHA, the BHA model having
  • the BHA model being configured to calculate (c) one or more points of contact between the BHA and a wall of the borehole and (d) forces imposed on the BHA at the one or more contact points, the forces comprising a side force vector imposed on the drill bit, and (e) drill bit tilt; calculating confined compressive strength of the rock formation using an axial motion drill bit model that receives drilling parameters for the drill rig, drill bit design information, and lithology information including formation rock strength; calculating lateral motion of the drill bit using a lateral motion drill bit model that receives (i) the calculated confined compressive strength, (ii) the lithology information, (iii) the drill bit design information, and (iv) the drill bit side force vector and drill bit tilt from the BHA model; calculating a ratio of lateral motion to axial motion; calculating an inclination angle and an azimut
  • a method of controlling drill bit trajectory in a subterranean formation may include receiving drilling parameters for operating a specific bottomhole assembly (BHA), constructing, with a computer processor, a directional drill-ahead simulator that may include a computer model of the BHA and the subterranean formation, calculating, with the computer processor, axial motion and lateral motion of a drill bit connected to a bottom end of the BHA using at least one formation parameter and at least one drilling parameter, predicting, with the computer processor, bit walk of the drill bit by accounting for and calculating contact forces and frictional forces between the BHA and a wall of a borehole in the subterranean formation using the computer model of the BHA, and determining, with the computer processor, an adjusted drill bit traj ectory to account for the predicted bit walk.
  • the method may include determining adjusted drilling parameters for operating the BHA to substantially follow the adjusted drill bit trajectory and operating the BHA according to the adjusted drilling parameters.
  • a method of planning and drilling a wellbore in a subterranean formation may include defining a target in a designated subterranean formation.
  • the method may include predicting wellbore spiraling and bit walk of a drill bit connected to a specific bottomhole assembly (BHA) that may include using a computer processor programmed to execute a directional drill-ahead simulator that may include a computer model of the BHA and the designated subterranean formation, receiving with the computer processor lithology data and drilling parameters for operating the BHA in the designated subterranean formation, calculating with the computer processor a ratio of lateral motion to axial motion using a lateral motion drill bit computer model and an axial motion drill bit computer model, predicting drill bit trajectory with the computer processor by accounting for and calculating lateral contact forces, angular displacement, and frictional forces using the computer model of the BHA, and adjusting the drill bit trajectory with the computer processor based at least in part on predictions from the computer model of the BHA.
  • the method may include adjusting the drilling
  • subterranean formation may include receiving lithology data for a specific subterranean formation, and receiving one or more drilling parameters for operating a bottomhole assembly (BHA), the one or more drilling parameters that may include weight-on-bit, torque, rotary speed, rate of penetration, drilling fluid flow rate, or lateral aggressiveness of a drill bit.
  • BHA bottomhole assembly
  • the method may include predicting wellbore spiraling and bit walk of a drill bit of the BHA that may include constructing, with a computer processor, a directional drill- ahead simulator that may include a dynamic computer model of the BHA and the subterranean formation, accounting for the wellbore spiraling and the bit walk of the drill bit by rotating a direction of normal contact forces on the drill bit by a constant angle using the computer processor, calculating combined forces on the BHA by adding torsional frictional forces to the normal contact forces in the dynamic computer model of the BHA in each iteration of a finite element analysis using the computer processor, predicting the drill bit trajectory, with the computer processor, based at least in part on calculating the combined forces on the BHA, and adjusting the drill bit trajectory, with the computer processor, based at least in part on predictions for the drill bit trajectory.
  • the method may include adjusting the one or more drilling parameters based at least in part on predictions for the wellbore spiraling and the bit walk of the drill bit.
  • FIG. 1 illustrates a cross-sectional view of an exemplary embodiment of a drillstring and drill bit disposed in a borehole penetrating the earth;
  • FIG. 2 depicts aspects of output from a bottomhole assembly (BHA) model
  • FIG. 3 illustrates inclination measurement convention
  • FIG. 4 illustrates a process flowchart for a method for predicting inclination of a borehole being drilled
  • FIGS. 5A and 5B collectively known as FIG. 5, illustrates a more detailed process flowchart for the method for predicting inclination of a borehole being drilled
  • FIG. 6 depicts aspects of a front steering angle model
  • FIG. 7 depicts aspects of the front steering angle model being applied to a BHA
  • FIG. 8 depicts aspects of a rear steering angle model being applied to the BHA
  • FIG. 9 is a flowchart for a method for predicting a path of a borehole that will be drilled in an earth formation by a BHA having a drill bit coupled to a drill tubular where the BHA is operated by a drill rig;
  • FIG. 10 depicts aspects of a frame of reference of the drill bit
  • FIG. 11 depicts aspects of a first cutter cutting formation rock
  • FIG. 12 depicts aspects of a second cutter cutting the formation rock
  • FIG. 13 depicts aspects of a third cutter cutting the formation rock
  • FIG. 14 depicts aspects of the formation rock after being cut by three cutters
  • FIG. 15 depicts aspects of pickets used as a virtual representation of the drill bit and the formation rock
  • FIG. 16 depicts aspects of a cutter face and cutter chamfer represented by a picket
  • FIG. 17 depicts aspects of rock cut by the cutter face and the cutter chamfer
  • FIG. 18 depicts aspects of a projected area of cut rock on a R-Z plane by a series of cutters on one cutter blade on the drill bit;
  • FIG. 19 depicts aspects of a chamfer cutter coming out of the page with two superimposed pickets;
  • FIG. 20 depicts aspects of a cutter represented by a picket at a rock interface
  • FIG. 21 depicts aspects of three orthogonal forces imposed on a cutter by a local rock surface interfacing with the cutter
  • FIG. 22 depicts aspects of circumferential (CIR) and side (SID) forces imposed on a cutter by a local rock surface interfacing with the cutter;
  • FIG. 23 depicts aspects of a relationship between a cutter and a local rock surface interfacing with the cutter
  • FIG. 24 depicts aspects of intrinsic specific energy required to cut rock as a function of cutter back rake angle
  • FIG. 25 depicts aspects of forces acting on a cutter
  • FIGS. 26 A through 26C depict aspects of various back rake angle of a cutter
  • FIG. 27 depicts aspects of force model parameter ⁇ as a function of back rake parameter ⁇ for Vosges sandstone
  • FIG. 28 depicts aspects of a comparison of predicted versus measured weight-on-bit in a first laboratory drilling simulator test
  • FIG. 29 depicts aspects of a comparison of predicted versus measured weight-on-bit in a second laboratory drilling simulator test
  • FIG. 30 depicts aspects of a comparison of predicted versus measured weight-on-bit in a third laboratory drilling simulator test
  • FIG. 31 depicts aspects of a comparison of predicted versus measured rate of penetration at an actual field location
  • FIG. 32 is a flowchart for a method for predicting an amount of axial motion of a drill bit having one or more cutters for drilling a borehole into a formation rock;
  • FIG. 33 is a flowchart for depicting aspects of a lateral motion drill bit model
  • FIG. 34 depicts aspects of a drill bit having gage pads configured to wear or crush formation rock and gage trimmers configured to cut formation rock;
  • FIG. 35 depicts aspects of a relationship of a gage pad with gage trimmers
  • FIG. 36 depicts aspects of a simulation of laterally removing the formation rock with the gage pad and the gage trimmers
  • FIG. 37 depicts aspects of lateral displacement as a function of axial displacement for the simulation
  • FIG. 38 depicts aspects of a comparison of predictions using the lateral motion drill bit model with measured laboratory side load drilling experiments in Bedford limestone;
  • FIG. 39 is a flowchart for a method for predicting a change in lateral displacement with a change in axial displacement of a drill bit drilling in a rock formation
  • FIG. 40 is a process flowchart for a directional drill-ahead simulator
  • FIG. 41 is a simplified drawing of a downhole view of the borehole shown in
  • FIG. 1 is a diagrammatic representation of FIG. 1 ;
  • FIGS. 42 A through 42K are a series of graphs of results of drilling simulations accounting for bit walk while building;
  • FIGS. 43 A through 43K are a series of graphs of results of drilling simulations accounting for right-hand bit walk while dropping;
  • FIG. 44A is a simplified drawing of a top view of a wellbore spiral
  • FIG. 44B is a simplified drawing of a side view of a wellbore spiral
  • FIG. 45 is a simplified transverse cross-sectional view of a BHA in a wellbore, and illustrates force vectors acting on the BHA.
  • FIG. 46 is a series of graphs of results of drilling simulations accounting for left-hand bit walk in build direction while turning.
  • a drill ahead method for predicting or simulating a geometry of a borehole that will be drilled or is being drilled based on operating parameters that may be applied to a drillstring by a drill rig.
  • the operating parameters can be selected so that the actual borehole has a desired geometry such as trajectory or path, bend radius, and final end location.
  • the computer processing system 12 may be programmed with a software program developed and employed commercially by Baker Hughes Incorporated of Houston, TX is known as the "Directional Drill Ahead Simulator (DDAS)."
  • DDAS Directional Drill Ahead Simulator
  • the drill ahead method uses a drillstring steering model (i.e., steering algorithm) to predict an inclination of the drill bit and thus the borehole at the point it is being drilled.
  • the DDAS may include additional models to predict axial motion and lateral motion and may further include models to predict drilling performance for a specific bit design in a realistic drilling environment.
  • dynamics of the drill bit are modeled using an axial motion drill bit model, referred to as DRILLBIT, and a lateral motion drill bit model referred to as SIDECUT, which are disclosed in U.
  • drilling parameters may be received that are specific to the BHA and specific to a defined target in a designated subterranean formation.
  • FIG. 1 illustrates a
  • FIG. 1 the drill tubular 5 is a drillstring made of drill pipes 6 serially coupled together.
  • a drill bit 7 is disposed at a distal end of the drill tubular 5.
  • the drillstring together with the drill bit 7 may be referred to as a bottomhole assembly (BHA).
  • a drill rig 8 is configured to conduct drilling operations such as rotating the drill tubular 5 and thus the drill bit 7 in order to drill the borehole 2.
  • the drill rig 8 is configured to pump drilling fluid (sometimes referred to as drilling mud) through the drill tubular 5 in order to lubricate the drill bit 7 and flush cuttings from the borehole 2.
  • a mud motor (not shown) may be included in the BHA. The mud motor is configured to convert energy from the drilling fluid into additional rotational energy for rotating the drill bit 7.
  • a steering device 9 is coupled to the drill tubular 5 and is configured to steer the drilling of the borehole 2 in a desired or intended direction using extendable pads 13 for example. Other steering configurations may also be used such as one that bends the BHA.
  • the intended direction may include an inclination direction (i.e., up or down with respect to the earth's surface) and/or an azimuthal direction (i.e., with respect to a reference azimuth such as true or grid north).
  • the steering device 9 is disposed close to the drill bit 7, within three feet (0.914 meters), for example, so that the force or portion of the force being applied to the drill tubular 5 by the steering device 9 is also being applied to the drill bit 7.
  • the steering device 9 in one or more embodiments may be considered as part of the BHA.
  • rotary steerable drilling systems may be used to steer the BHA through the borehole 2.
  • the RSS may include AutoTrakTM eXact high-build rotary steerable system (RSS) ("AutoTrakTM eXact RSS”) commercially available from Baker Hughes
  • a controller 11 is configured to control the steering device 9 to steer the drilling of the borehole 2 in the desired direction.
  • the controller 11, which may include downhole electronics, may also act as an interface with telemetry to communicate data or commands between downhole components and a computer processing system 12 disposed at the surface of the earth 3.
  • Non-limiting embodiments of the telemetry include pulsed-mud and wired drill pipe.
  • System operation, control and/or data processing operations may be performed by the controller 11 , the computer processing system 12, or a combination thereof.
  • the BHA may include a sensor 10 configured to sense various downhole parameters, which may be transmitted uphole according to programmed instructions to the computer processing system 12 for data recording, processing or display to be utilized by an operator or analyst to control the drilling operations.
  • the computer processing system 12 may be configured to accept inputs (e.g., via the sensors 10 or via a user input device) including simulating drilling operations, to improve aspects of an active drilling operation through corrective measures comprising alteration of operating parameters, to provide recommendations for equipment selection and operation for a subsequent, planned drilling operation, or both.
  • the sensor signals may be provided at selected time intervals, at depth intervals along the drill path, at reduced intervals during drilling of nonlinear portions of the borehole, or a combination thereof.
  • the controller 1 1 may receive the signals from the downhole sensors 10 and any other sensors used in the drilling assembly and process the signals according to programmed instructions.
  • the controller 1 1 may send the results of the processed signals (e.g., current downhole conditions, current position, position relative to the predetermined drill path, current operating parameters, recommended operating parameters, current equipment deployed, and recommended equipment for deployment) to an electronic display of the computer processing system 12, which may be utilized by an operator to control the drilling operations.
  • measured characteristics of the earth formation 4 may be utilized by the computer processing system 12.
  • the measured characteristics may be extrapolated using the computer processing system 12 by accessing a database of characteristics of geographically closest earth formations and accepting estimated characteristics of the earth formation 4 (e.g., using linear, polynomial, or other known extrapolation techniques) at the computer processing system 12.
  • estimated characteristics include BHA inclination, BHA acceleration, and formation parameters that may be logged such as mineralogy.
  • the teachings disclosed herein may be implemented in real time by the computer processing system that receives sensor data or the teachings may be implemented by another computer processing system that does not receive sensor data in real time.
  • various analysis components may be used, including a digital and/or an analog system.
  • the steering device 9, the downhole sensor 10, the controller 1 1, or the computer processing system 12 may include digital and/or analog systems.
  • the computer processing system 12 may have components such as a processor, storage media, memory, input, output, communications link (wired, wireless, pulsed mud, optical or other), user interfaces, software programs, signal processors (digital or analog) and other such components (such as resistors, capacitors, inductors and others) to provide for operation and analyses of the apparatus and methods disclosed herein in any of several manners well appreciated in the art.
  • the signal receiving device may be a display monitor or printer for presenting the result to a user.
  • the signal receiving device may be memory or a storage medium. It can be appreciated that storing the result in memory or the storage medium will transform the memory or storage medium into a new state (containing the result) from a prior state (not containing the result). Further, an alert signal may be transmitted from the processor to a user interface if the result exceeds a threshold value.
  • drilling a subterranean formation is often a three-dimensional process, as the drill bit not only penetrates the formation linearly along a vertical axis, but is either purposefully or unintentionally drilled along a curved path or at an angle relative to a theoretical vertical axis extending into the subterranean formation in a direction substantially parallel to the gravitational field of the earth.
  • drilling operation means and includes any operation performed during the formation or enlargement of a borehole in a subterranean formation.
  • drilling operations include drilling, reaming, and other formation removal processes.
  • earth-boring tool means and includes any type of tool used for earth removal during the formation or enlargement of a borehole in a subterranean formation and includes, for example, fixed-cutter (i.e., "drag") bits, roller cone bits, percussion bits, core bits, eccentric bits, bicenter bits, reamers, mills, hybrid bits, and other drilling bits and tools known in the art.
  • fixed-cutter i.e., "drag” bits
  • roller cone bits i.e., percussion bits
  • core bits i.e., eccentric bits
  • bicenter bits reamers, mills, hybrid bits, and other drilling bits and tools known in the art.
  • BHA model relates to a finite element model or beam model that models the dynamics of the BHA (i.e., drillstring and drill bit) in a borehole.
  • the user can create a hole geometry by specifying the length of a hole section and the inclinations at each end of the hole section.
  • Each section may be represented by a series of circles where each circle represents a cross-section of the borehole.
  • the model automatically represents sections of the BHA between the two ends of the hole section that the user has specified. Once a hole is created, a drillstring or BHA can be created and put in the hole.
  • the finite element model or beam model computes how the drillstring lays in the hole, where it contacts the hole, how the drillstring is curved or flexed and magnitude and direction of the contact forces where it touches the borehole wall based on the operating parameters applied to the BHA by the drill rig.
  • the BHA model can also calculate the bit tilt of the drill bit, which is the angle between the longitudinal axis of the drill bit and the longitudinal axis of the borehole or the difference between the inclination of the drill bit and the inclination of the borehole.
  • One example of a BHA model developed and employed commercially by Baker Hughes Incorporated of Houston, TX is known as the BHASYS PRO.
  • FIG. 2 illustrates an example of output provided by the BHA model, which illustrates the direction and magnitude of the contact forces and the geometry of the BHA. In that BHA models are known to persons skilled in the art, further details of these models are not discussed in further detail.
  • Inclination ⁇ is measured from the vertical depth axis as illustrated in FIG. 3.
  • BUR is measured as change in inclination with respect to measured depth. In other words, it is the "hole curvature” projected onto a vertical plane.
  • Dog Leg Severity (DLS) is the change in hole angle with respect to measured depth. It is simply the "curvature of the hole,” regardless of how the hole is oriented. It could be curving to the side, or curving upward, etc.
  • the inclination will be denoted with the symbol ⁇ .
  • the subscript "new" used with Gnew refers to the inclination at the very end (bottom of the hole) of the survey section that is being created by the predictive BUR algorithm.
  • Gold refers to the survey section at the top of the newly created survey section. It is called “old” because in BHASYS PRO, the inclinations where survey sections join are always equal to one another.
  • the inclination at the bottom of the hole for survey section N will be equal to the inclination at the top of the survey section for survey section N + 1.
  • a directional drill ahead simulation method 40 for adding a new section of borehole to the BHA model is discussed. Aspects of the method 40 are depicted in FIG. 4.
  • the BHA model gives its forces, bit tilt, moments, and curvature information to the axial and lateral motion drill bit models. It also gives hole curvature and BHA alignment within the hole to the steering model.
  • the axial motion drill bit model runs using input from formation lithology logs and drill rig operating parameters. It calculates a rock strength (e.g., confined compressive strength) and a rate of penetration (ROP). Then, the ROP and rock strength are given to the lateral motion drill bit model.
  • a rock strength e.g., confined compressive strength
  • ROP rate of penetration
  • the lateral motion drill bit model uses this data and the data such as bit side force and bit tilt, which was supplied by the BHA model, to calculate the ratio of lateral ROP to axial ROP or lateral displacement to axial displacement for the drill bit for a certain time period. This is called dL/dZ.
  • a front steering angle, 5_front, of the BHA is derived from dL/dZ and is defined as:
  • the steering model is next invoked (explained in detail below).
  • the steering model calculates inclination and azimuth for the end of a new hole section to be created by simulated drilling using 5front.
  • a new hole section is created based on the calculated inclination and azimuth.
  • the BHA and bit are moved down to the bottom of the new hole section and the three models (BHA, axial motion drill bit model and lateral motion drill bit model) are invoked again.
  • the drill ahead simulation method 40 runs the BHA model first, then feeds its results to the axial and lateral drill bit models, which are incorporated into the drill ahead simulation method 40.
  • the drill ahead simulation method 40 then predicts the location and geometry of the next hole section.
  • There are alternate ways or sequences of running the models but the main point is that the steering model is fed information by the other models, and a new simulated hole section is placed in the direction that the bit/BHA system wants to drill as calculated by the steering model.
  • FIGS. 5A and 5B are a flowchart depicting further details of the drill ahead simulation method 40.
  • the general directional behavior of the drill bit will be governed or influenced by the following factors: WOB/ROP relationship; lateral bit force/lateral ROP relationship; drill bit tilt, drill rotary speed (e.g., RPM);
  • CCS confined compressive strength
  • the drillstring steering model includes a front steering angle model. Steering accuracy may be improved by adding a rear steering angle model to the front steering angle model.
  • the front steering angle model is based on the kinematic motion of a two-wheeled vehicle, such as a bicycle.
  • FIG. 6 illustrates a simple diagram of such a vehicle in a top view and aspects related to a front steering angle.
  • the following equation describes the angular rate of change with respect to time of the vehicle as a function of the velocity V of the vehicle, the length L of the wheel base, and a front steering angle ⁇ .
  • the steering angle is the angle between the vehicle's longitudinal axis and the direction that the front wheel is pointing.
  • the front steering angle model makes an analogy between this two-wheeled vehicle and the portion of the BHA from the bit to the first point of contact with the borehole wall.
  • One difference between a bicycle path and directional drilling is that the back wheel of a turning bicycle does not follow the front wheel's path.
  • the second point of contact of the BHA does follow the path of the bit, which makes the first point of contact with the formation. For simplicity, this difference is neglected at this point.
  • FIG. 7 depicts aspects of the front steering angle as applied to the BHA. The following parameters are used to describe the BHA drilling the borehole:
  • the steering angle ⁇ The steering angle of the bit. (Note that this is not bit tilt. No physical part on the bit has necessarily been tilted.
  • the steering angle ⁇ can be conceptually thought of as the angle between: the axis of the tool (from bit to first point of contact after the bit) and the instantaneous trajectory of the bit.)
  • the steering angle of the BHA is influenced by certain drilling and lithology parameters. All of these parameters do change in time during a normal drilling run. For now 5(t) will be denoted as ⁇ , with the assumption that it can change throughout the run as time or depth change because the parameters mentioned above will change with time or depth.
  • L denote the distance from the drill bit to the first point of contact of the BHA after the drill bit
  • ⁇ /ds is the change in bit axis inclination (and hole inclination) with respect to depth.
  • the curvature of a circular section is mathematically defined as: ds ⁇ ⁇ R) ⁇ k where R is the radius of the circular survey section and k is the curvature of the circular survey section (i.e., radius of curvature of the borehole at the survey section).
  • This equation can be manipulated to find the new inclination angle 0 new for the new survey section that the drill ahead simulation method predicts in an iteration.
  • 0 o n is the inclination at the hole bottom prior to the creation of the new survey.
  • 0 mw may then be written using the front angle steering model as:
  • FIG. 8 depicts aspects of the rear steering angle model.
  • the BHA has an imaginary steering front wheel and an imaginary steering rear wheel.
  • the rear wheel must follow the path of the borehole already drilled as if the rear wheel was guiding the BHA at the steering axis of the rear wheel.
  • the rear steering angle model is very similar to the front steering angle model. The only two differences are that now a rear steering angle appears in the formula immediately above, and the cosine of the front steering angle modifies the prediction as well in the following equation: [tan ⁇ from- tan 5 rear ]
  • the rear steering angle S rear is the angle between the borehole inclination #at the first point of contact of the BHA after the drill bit contact and the tool axis, or it can be defined as the angle between the borehole inclination #at the first point of contact after the drill bit and the borehole inclination ⁇ where the bit contacts the wall. Using the latter definition, the rear steering angle may be written as:
  • FIG. 9 is a flowchart for one example of a method 90 for predicting a path of a borehole that will be drilled in an earth formation by a bottomhole assembly (BHA) having a drill bit coupled to a drill tubular, the BHA being operated by a drill rig.
  • Block 91 calls for constructing a BHA model of the BHA.
  • the BHA model includes (a) dimensions, geometry, mass distribution, material density, and material stiffness of the BHA and (b) dimensions and geometry of the borehole as predicted to be drilled.
  • the BHA model is configured to calculate (c) one or more points of contact between the BHA and a wall of the borehole and (d) forces imposed on the BHA at the one or more contact points, the forces comprising a side force vector imposed on the drill bit, and (e) drill bit tilt.
  • Block 92 calls for calculating confined compressive strength of the rock formation using an axial motion drill bit model that receives drilling parameters for the drill rig, drill bit design information, and lithology information comprising formation rock strength.
  • Block 93 calls for calculating lateral motion of the drill bit using a lateral motion drill bit model that receives (i) the axial ROP from the axial motion drill bit model, (ii) the lithology information, and (iii) the drill bit design information, and (iv) the drill bit side force vector and drill bit tilt from the BHA model.
  • Block 94 calls for calculating a ratio of lateral motion to axial motion.
  • Block 95 calls for calculating an inclination angle and an azimuthal direction of the BHA using a BHA steering model that receives the ratio. The inclination angle together with the azimuthal direction (e.g., direction or angle with respect to true north) provides a three-dimensional direction.
  • Block 96 calls for iterating the above steps by updating the BHA model to include extending the borehole an incremental distance in the direction of the inclination angle and the azimuthal direction and displacing the BHA the incremental distance in the extended borehole.
  • Each of the blocks in the method 90 may be implemented by a processor such as in a computer processing system.
  • data used by the method 90 as input to the various above discussed models may be updated in real time as an actual borehole is being drilled according to drilling parameters that were used to obtain the predicted borehole path.
  • the updated data may be obtained from one or more sensors disposed on a BHA that drills the actual borehole.
  • Sensors may include a borehole caliper sensor and/or a formation sensor that is configured to sense data from which the formation lithology may be derived. Examples of formation sensors include a natural gamma-ray sensor and neutron tools that emit neutrons and sense neutrons or gamma-rays resulting from neutron interactions with the formation.
  • the axial motion drill bit model (i.e., DRILLBIT) is discussed.
  • the axial bit model is part of the Directional Drill-Ahead Simulator. It allows for predicting drilling performance for a specific bit design in a realistic drilling environment and, in particular, uses a PDC cutter force model to calculate the forces on the drill bit. Given operating parameters, formation characteristics and a specific bit design, the axial model predicts the rate of penetration (ROP) if the weight-on-bit (WOB) is prescribed or predicts the WOB if ROP is prescribed.
  • ROP rate of penetration
  • WOB weight-on-bit
  • the overall approach is as follows: (a) from bit design and operating parameters, compute the areas of cut (projected onto a vertical plane) across the face and chamfer for each of the cutters representing the bit; (b) from details of the rock formation, drilling depth and mud weight, estimate the confined compressive strength of the rock formation; (c) given these estimated areas of cut, the detailed geometry of how the edge of the cutter is engaging with the rock and the CCS of the rock formation, calculate forces on the face and chamfer of the cutters using a "force model" (discussed below); (d) and sum the forces on all of the cutters to find the net force on the drill bit.
  • ROP is prescribed, these forces are easily translated to a WOB (weight-on-bit) and TOB (torque-on-bit). If WOB is prescribed, an iterative procedure is used by which ROP is adjusted until the predicted WOB matches the prescribed WOB.
  • the Directional Drill-Ahead Simulator (DDAS) is drilling with a prescribed ROP
  • the CCS and some information on the bit design are passed back and used by the DDAS.
  • SIDECUT the lateral motion drill bit model
  • the axial motion drill bit model assumes that the bit is drilling "on-center" and the bit axis of revolution is always coincident with a fixed Z axis in the rock frame of reference as illustrated in FIG. 10.
  • the Z axis lies on the bit axis of rotation and positive Z is in the direction of nose to shank.
  • the rock surface is defined in this plane as a simply-connected line defined in this vertical plane.
  • FIG. 10 illustrates this frame of reference with the vertical plane superimposed.
  • the drill bit in FIG. 10 is a poly crystalline diamond compact (PDC) bit.
  • drill bit rotation and rock surface update are now discussed.
  • the drill bit and rock surface are modeled using nodes and the node locations are updated based on the drill bit-rock interaction.
  • the cutters on the bit are first sorted in terms of increasing angle around (angular position on the bit). Consider a bit with 3 cutters and parameters listed Table 1.
  • the rock surface is updated during one revolution of this bit as follows: (1) the intersection of cutter 1 with the rock is calculated and the rock surface is updated accordingly; (2) the bit is rotated through an angle (45-0) and, based on RPM and ROP, the bit vertical position is modified accordingly; (3) the intersection of cutter 2 with the rock is calculated and the rock surface is updated accordingly; (4) the bit is rotated through an angle (270-45) and, based on RPM and ROP, the bit vertical position is modified accordingly; (5) the intersection of cutter 3 with the rock is calculated and the rock surface is updated accordingly; and (6) the bit is rotated through an angle (360-270) and, based on RPM and ROP, the bit vertical position is modified accordingly. This sequence is repeated until run termination criteria (e.g., a number of iterations) are met.
  • FIGS. 11 through 14 illustrate the progression of the rock surface update.
  • the axial model is computationally fast because in one or more embodiments it is required only to keep track of the rock surface on one 2D plane and calculate cutter interaction with the rock surface on this plane. This is possible for the situation where the bit is drilling on-center. As a cutter passes the X-Z plane, it clips off the tops of the pickets in order to update the represented rock surface as illustrated in FIG. 15.
  • the PDC cutters' edges and chamfers are represented initially as ellipses (circles with non-zero back rake projected onto a vertical plane) as they pass through the rock "plane.” These are exact (analytic) representations for an unworn cutter with the design back rake and side rake.
  • the rock surface is also broken up into an identical set of pickets that match the location (along the X-axis) of the cutter pickets as illustrate in FIG. 15.
  • FIG. 16 illustrates this configuration.
  • FIG. 16 illustrates this configuration.
  • FIG. 16 shows the cutter face (with chamfer) with a non-zero back rake coming out of the page at the viewer. If the cutter is allowed to wear, that edge of the cutter in contact with the rock will deviate from an ideal ellipse. This deviation, applied on a picket-by -picket basis, is governed by a wear model discussed in the disclosure further below.
  • the areas of cut are delineated in FIG. 16.
  • the upper patch is the (projected) area of cut on the cutter face and the lower patch is the (projected) area of cut on the chamfer. It is noted that these areas are dependent upon the location of the rock surface with respect to the cutter at a particular picket.
  • the rock surface may only intersect part of a chamfer and none of the face at a picket location.
  • the areas of cut on the chamfers and face are kept track of independently as each will behave differently regarding cutter forces. This is also discussed later. It is also possible for a cutter picket to be completely worn through and no longer exist.
  • FIG. 17 shows the updated rock surface for a new cutter. If the cutter is worn, the updated rock surface will simply be the bottom-most end of a picket when the picket takes a cut. If the picket does not take a cut, the rock surface at that location does not change.
  • FIG. 18 illustrates an example. This example represents the areas of cut for all cutters on a single blade of a PDC bit. Areas of cut on a cutter, and the forces associated with these areas of cut, depend both on the cutter layout on the bit and operating parameters (e.g., RPM, ROP). After a few revolutions of the bit, the magnitude and shape of the areas of cut stabilize to constant values for a fixed ROP and RPM. Note that the image in FIG.
  • the effective back rake of a cutter is now discussed.
  • the effective back rake is an angle that plays a role in cutting forces.
  • Effective back rake is not the design back rake found in a bit design file. The two are associated but that is all. Effective back rake is the angle between the (normal to the) local cutting surface and the local rock surface.
  • FIG. 19 illustrates a chamfer cutter coming out of the page. Two pickets are superimposed on the image. The two large arrows represent the surface normals of the cutter chamfer at the locations of the pickets. The small arrow is the surface normal at the cutter face at the location of one of the pickets. The dots represent vectors coming out of the page and are parallel with the local rock surface that the cutter will see as it passes through the rock. The effective back rake (EBR) is that angle between this vector and the arrows. As is evident, the EBR can depend upon location on the cutter, ... both along the edge of the cutter and whether or not the cut is on the chamfer or on the face.
  • FIG. 20 is a close-up of a bottom of a picket edge (at an edge of a cutter).
  • the edge of the cutter is delineated by the contact between the picket and the triangle.
  • the vector arrow into the cutter is the normal to the local cutter surface (either on a chamfer or on the face). Note that this normal vector is in three dimensions and has three components, each associated with an axis in a right-handed coordinate system. There is a component of the normal vector coming out of the page.
  • the angle ⁇ is referred to as the "cut" angle.
  • a simple two-dimensional (2D) rotation operator is constructed that will rotate the bottom edge of the picket into the horizontal plane. This operator is applied to the three-dimensional (3D) normal to the local cutter surface (red vector depicted in FIG. 20). If the resulting rotated normal vector (normalized) has components (Nx, Ny, Nz), where X and Z have been defined previously Frame") and Y is "out of the screen, then the effective back rake (EBR) and effective side rake (ESR) are given by:
  • ESR SIN ⁇ l (N x ) .
  • the forces on a cutter are computed by summing the forces on individual pickets that span the cutter. Forces on an individual picket will depend upon the area of cut on the chamfer, the area of cut on the face, the effective back rake, the rock formation, confined compressive strength (CCS), drilling depth and mud weight. Details of the force model are given in a different section. Once the normal (Fn) and circumferential (Ft) forces are computed on a picket, they are transformed back into the bit frame of reference. The net force on the bit is a summation of the individual cutter forces which are themselves the summations of individual picket forces. Consider FIG.
  • Fn is always normal to the local rock surface (regardless of the orientation of the local rock surface with respect to the earth).
  • Ft is in the direction normal to the X-Z plane and Fr is perpendicular to the plane containing Fn and Ft.
  • TRQ (torque on picket) R P i C CIR (see below for definition of R P i C )
  • helix angle is computed for each picket and is given by the following equation: where ROP is penetration rate (ft/hr), RPM is bit rotary speed, and R P i C is radial position of picket (distance from bit axis to picket in horizontal plane).
  • FIG. 22 These force components on a picket are illustrated graphically in FIG. 22.
  • AA is the angle around of the cutter and CIR and SID are defined as above. Note that the view is from a perspective looking down onto the horizontal plane from shank to nose.
  • the forces on a picket are given by the following equations.
  • the total force on a cutter is simply a summation of these forces for all (N) pickets associated with that cutter, as given by:
  • N is the number of pickets associated with a particular cutter.
  • the net force and moment on the bit is thus a summation of forces and moments on all of the cutters on the bit.
  • the forces on a PDC cutter are given by two orthogonal components: a "tangential" component that is parallel with the local rock surface and a "normal” component that is perpendicular to the local rock surface. Note that when evaluating such things as weight-on-bit and torque-on-bit, Fn and Ft must be transformed to the bit frame of reference.
  • the model is defined by the following equations:
  • ⁇ ratio of Fn to Ft (dependencies to be defined).
  • FIGS. 23 and 18 define the PDC cutter in relationship to the rock.
  • FIG. 23 is a side view of a cutter moving through the rock, with Fn and Ft indicated.
  • FIG. 18 shows one blade of a PDC bit with the cutters coming out of the page.
  • FIG. 18 illustrates what in meant (in the following) by the projected areas of cut (blackened areas) associated with the rock removed.
  • the cutter force model is extended to a practical form by identifying the nature of "intrinsic specific energy” and developing a means of characterizing its value, on a site-specific basis, from well logs, drilling depth and mud weight.
  • the intrinsic specific energy term ⁇ is the effective cutting strength (or resistance to be cut) of the rock. This effective strength will depend upon the following: wellbore pressure; hole depth; mud weight; lithology; cutter orientation on the bit; cutter geometry (face/chamfer); and cutter interfacial friction coefficient.
  • the cutter orientation on the bit and the orientation of the formation rock with respect to the bit provide an orientation of the cutter with respect to the formation rock to be cut, which is accounted for in the intrinsic specific energy as disclosed herein.
  • the confined compressive strength is taken to be the unconfined compressive strength (UCS) modified by a term associated with the confining pressure and is defined as follows:
  • CCS confined compressive strength
  • UCS unconfined compressive strength
  • FA internal angle of friction
  • CP confining pressure
  • the intrinsic specific energy cannot be equated to the CCS because ⁇ can be backed out as a function of cutter back rake ⁇ as illustrated in FIG. 24 for Vosges sandstone under atmospheric conditions. From FIG. 24, it is clear that the strength of the rock increases with increasing cutter back rake. This is associated with the additional confinement of the rock by the cutter.
  • a 2 is a multiplicative parameter included in order to allow adjustment during model fitting to lab or field data.
  • FIG. 25 illustrates the forces acting on a cutter.
  • the cutter back rake is given by ⁇ and ⁇ is the inclination due to friction of the net cutting force from the perpendicular to the cutting face. From FIG. 25, the ratio may be expressed as:
  • ⁇ ( ⁇ , ⁇ ) ⁇ ( ⁇ + ⁇ ) .
  • the effective back rake ⁇ is a well-defined parameter. ⁇ on the other hand is not so straight-forward. Depending upon the degree of back rake, the rock can "flow" upward along the cutter face or downward beneath the cutter as illustrated in FIGS. 26A through 26C. Based on laboratory experiments on Vosges sandstone, the relationship of ⁇ with back rake ⁇ is illustrated in FIG. 27.
  • the confining pressure (CP) model is one of the key elements of a successful "drill ahead" simulation.
  • CP ⁇ [( ⁇ D eff -0.05)/0.15] + BHP[ (0.2- ⁇ ⁇ 1 ⁇ )/0.15].
  • ⁇ is the measured porosity of the formation or porosity inferred from well log data. Note that the last entry above is an interpolation between the maximum allowed CP (namely BHP) and the minimum allowed CP (namely ⁇ ) and that this interpolation is based on effective porosity.
  • a mud weight correction factor is implemented in DRILLBIT.
  • the correct interpretation of this effect is via the "chip hold-down effect" that has to do with the dislodgement of a chip created in the cutting process. If fluid cannot migrate between the produced chip and the formation from whence it came, the chip is subject to the entire load of the fluid column and the formation is tougher to drill. However, if fluid can migrate into the fractures and equilibrate the stresses on the chip, the chip can be removed easily. Since a fluid migrates through generated cracks in this mechanism, the migration will depend upon mud viscosity (related to mud weight) and will be independent of whether the formation is initially permeable or impermeable.
  • the mud correction factor is:
  • the force model is applied independently to the cutter chamfer and the cutter face. Depending upon the operating parameters all or part or none of the chamfer pickets and face pickets will be engaging the rock.
  • the associated "bookkeeping" for this is an integral part of the axial bit model. Many well logs do not give a specific lithology (sandstone, limestone, etc.) at a particular depth but rather a mix of lithologies in terms of fractional content. In general, the sum of the fractions at a given depth is unity.
  • FIGS. 28 through 30 represent a comparison of predicted and measured weight on a few laboratory drilling simulator tests. During each drilling test, the prescribed WOB was changed and the measured ROP was recorded. The upper and lower bounds in these plots were computed by propagating uncertainties in porosity, friction angle and UCS through DRILLBIT.
  • the drilling fluid was water and the rock was Carthage limestone.
  • the drilling fluid was 11 lb. (4.990 Kg) water-based drilling mud and the rock was Carthage limestone.
  • the drilling fluid was 16 lb. (7.257 Kg) oil-based drilling mud and the rock was Mancos shale.
  • FIG. 31 Another example is a comparison of predicted with measured ROP at an actual field location as illustrated in FIG. 31.
  • the legend for this plot is: SHA: fraction of shale; SSA: fraction of sandstone; LSO: fraction of limestone; DOL: fraction of dolomite; POR: formation porosity; CP: predicted confining pressure; and UCS: unconfined compressive strength.
  • FIG. 32 is a flowchart for a method 320 for predicting an amount of axial motion of a drill bit having one or more cutters for drilling a borehole into a formation rock.
  • Block 321 calls for receiving lithology data for the formation rock.
  • Block 322 calls for receiving drilling parameters that include mud weight and depth of drill bit for a drill rig that will operate a drill tubular coupled to the drill bit and drill bit design information.
  • Block 323 calls for calculating confined compressive strength of the formation rock using the received lithology data, drilling parameters, and drill bit design information.
  • Block 324 calls for calculating an area of cut into the formation rock by each of the one or more cutters.
  • Block 325 calls for calculating an effective back rake of each of the one or more cutters.
  • Block 326 calls for calculating a force applied to each of the one or more cutters using the calculated confined compressive strength of the formation rock, the area of cut and the effective back rake of the one or more cutters, and by accounting for an orientation of each of the one or more cutters with respect to a surface of the rock to be cut.
  • Block 327 calls for summing the calculated forces applied to each of the one or more cutters to calculate the WOB and TOB (torque-on-bit).
  • the above method may be implemented by a processor. If ROP is prescribed, then the method 320 may be used to calculate the WOB and TOB that provides the prescribed ROP. Alternatively, if WOB is prescribed, then the ROP and TOB may be calculated to provide the prescribed ROP.
  • the method 320 may call for representing the drill bit with a first set of virtual pickets and the formation rock with a second set of virtual pickets.
  • the first set of pickets lines up with the second set of virtual pickets at the interface of the drill bit with the formation rock. That is, if the first set of pickets was extended downward and the second set of pickets was extended downward, the two sets of pickets would overlap.
  • the area of cut, the effective back rake, and the force applied to each cutter are calculated on a picket-by-picket basis, and the WOB and TOB are calculated by summing the forces on each of the pickets.
  • SIDECUT is a program/procedure for calculating the lateral migration dL (with respect to the borehole) of a drill bit drilling a rock formation as a function of drilling depth (dZ) along the borehole.
  • FIG. 33 is a flowchart depicting aspects of the SIDECUT procedure.
  • the drill bit can have a tilt (with respect to the borehole) and an applied side load that forces the drill bit against the borehole wall.
  • SIDECUT may be used for PDC bits, roller cone bits, and hybrid bits combining aspects of PDC bits and roller cone bits.
  • FIG. 34 depicts aspects of a PDC bit in three-dimensional view.
  • the procedure allows for arbitrary values of: bit rate of penetration (ROP) through the rock formation; bit rotary speed (RPM); bit tilt with respect to the borehole side load applied to the drill bit; bit geometry (gage pad length, gage pad recess, fraction of bit circumference occupied by gage pads); and confined compressive strength (CCS) of rock formation.
  • the CCS is calculated in DRILLBIT, which provides the CCS to SIDECUT.
  • gage pad recess defined as the radial distance that gage pads are recessed inward from the position of the outermost edge of the ground gage (PDC) trimmers, can vary with respect to vertical position on the gage pad. This allows for tapered, stepped, and other gage pad geometries. FIG. 35 illustrates these aspects in a two-dimensional representation.
  • SIDECUT uses a two-dimensional (2D) description of the drill bit and rock formation. The reason for this is computation speed. However, a full three-dimensional (3D) model can be implemented using more computational time. There is excellent agreement between the simplified 2D model and the more sophisticated 3D model.
  • the 3D characterization of the PDC bit illustrates the PDC cutters (small cylinders) and the gage pads (rectangular structures).
  • the collection of PDC cutters and associated gage pads is referred to as a "blade.”
  • FIG. 35 the 2D characterization illustrates a single blade and shows the gage pad and PDC cutters ("ground gage trimmers") below the gage pad. The blade is being forced to the right against the rock with an applied side load.
  • the dotted line in FIG. 35 is the initial rock surface in SIDECUT. The rock is simply "shrink wrapped" to the bit in order to initialize the rock surface and is indicated by the filled region to the right. Note that the outer surface of the gage pad in FIG.
  • gage pad recess is significant in governing the directional drilling properties of a PDC drill bit.
  • the vertical (Z) extent of the gage pad is called the “gage pad length” and the fraction of the drill bit circumference that is occupied by gage pads is simply referred to as “fraction” in the following discussion.
  • SIDECUT accounts for both bit geometry and rock wear in computing the lateral migration of a bit in a bore hole.
  • Rock wear is estimated using a sliding wear model. The amount of rock worn off at a particular position on the bore hole wall is proportional to the total sliding distance of the gage pads across that position as well as the contact stress applied to that location and the confined compressive strength of the rock formation.
  • Sliding distance is governed by the bit RPM, ROP, gage pad length and "fraction.”
  • Contact stress is governed by the applied side load and instantaneous rock-pad contact area.
  • Drill Bit Construction Construct the 2D bit by constructing a series of equally spaced (vertically) points or nodes from the top of the gage pad to the bit nose. Differentiate those points associated with the gage pad from those point associated with the PDC cutters. Let the vertical spacing be a prescribed value DZ. The totality of PDC Cutters below the gage pads are represented by a single vertical construction with a prescribed back rake. This series of equally spaced nodes can have an arbitrary shape in R-Z space (see FIG. 35). This generality can therefore encompass tapered and stepped gage pads. The methodology also can allow for a 2D bit profile although the implementation described herein simply assumes a continuous vertical blade as representative of the PDC cutting structure.
  • the initial rock surface is a series of equally spaced (vertically) points or nodes that is a clone of the initial drill bit representation.
  • the rock node locations are initially identical to the bit node locations.
  • K stiffness
  • the rock stiffness can depend upon rock properties or other parameters as provided by a lithology log, for example.
  • the associated equation for pad-rock contact stress is:
  • the processing loop includes three procedures discussed in the following: force balancing, rock removal and hole growth. At each vertical position of the drill bit with respect to the rock, all three of these procedures are applied in the order given below for this example. The bit is then moved by an amount DZ (see FIG. 15) and the process is repeated until some prescribed drilling distance is achieved. This distance is generally a multiple of the gage pad length. A. Balance restoring force of rock on bit to equal the applied side load on bit.
  • the contact force on a single pad node is:
  • AOC A dZ
  • MSE mechanical specific energy
  • CCS the confined compressive strength
  • AOC area of cut
  • EBR effective back rake
  • ⁇ depth of cut (penetration depth)
  • A2 and A3 force model calibration coefficients.
  • Rock removal is done on a node-by-node basis for all of the nodes representing the bit (pads and cutters).
  • Rock removal by pads On a node-by node basis on the set of rock nodes that are adjacent to the pads, check to see if the contact stress exceeds the CCS. If so (and this is optional by the user), the rock node can be adjusted to simply reflect displacement. That is, the new R location of the rock node is at the R position of the respective pad node. If the contact stress is less than the CCS, the new rock node R position is governed by the sliding wear model. The amount ⁇ a rock node is moved in the +R direction due to wear is given by:
  • the final piece of the computation is to grow the borehole by the amount dZ. This is done only at the end of going through all of the rock nodes and removing rock according to the previous section. Recall that dZ is the prescribed vertical node spacing for both the bit and the rock. So once the rock is removed according to the previous section, a single rock node is added. If at this point there are N rock nodes, then:
  • R (N + 1) R (N) + dZ tan(TILT)
  • the indices of the rock are simply adjusted.
  • the penetration depth ⁇ of the pad at rock node k is:
  • FIG. 35 illustrates rock and drill bit geometry half-way through a SIDECUT simulation of drilling a borehole.
  • FIG. 36 illustrates a graph of axial depth versus lateral side cut for the simulation illustrated in FIG. 35.
  • FIG. 37 illustrates results of force balancing at each depth step during a SIDECUT simulation with a prescribed side load of 800 lbs.
  • FIG. 37 illustrates the sensitivity of gage pad contact area to bit RPM and gage pad length.
  • FIG. 38 illustrates aspects of a comparison of SIDECUT predictions with measured laboratory side load drilling experiments in Bedford limestone.
  • FIG. 39 is a flowchart for a method 390 for predicting a change in lateral displacement with a change in axial displacement of a drill bit drilling in a formation rock.
  • Block 391 calls for constructing a virtual representation of the drill bit and the formation rock, the drill bit comprising a gage pad configured to remove rock by wearing or crushing the rock during moving contact and a cutter configured to cut into rock during moving contact.
  • Block 392 calls for adjusting lateral penetration depth until rock reactive force equals a side force applied to the drill bit to provide a lateral displacement of the drill bit to provide an adjusted lateral penetration depth.
  • Block 393 calls for removing formation rock to where the pad and the cutters contact the formation rock at the adjusted lateral penetration depth and to a selected axial displacement of the drill bit.
  • Block 394 calls for moving the drill bit in a drilling direction to the end of the currently drilled borehole.
  • Block 395 calls for iterating the adjusting, the removing, and the moving to predict the change in lateral displacement with the change in axial displacement of the drill bit.
  • the method 390 may also include using a first set of nodes as representation of the drill bit and a second set of nodes as a representation of the rock formation.
  • the method 390 may further include setting the selected axial displacement as the axial distance between two nodes in the axial direction.
  • the method 390 may further include initializing the nodes such that the nodes at the boundary of the drill bit with the rock formation overlap.
  • DDAS Downlink Bit Walk Model
  • Constant Bit Walk Model supplements the drillstring steering model of the DDAS by changing the contact force (i.e., the applied side load) by rotating the contact force by a constant angle.
  • a “Torsional Friction Bit Walk Model” supplements the drillstring steering model of the DDAS by accounting for torsional friction between components of the drillstring and the walls of the wellbore in calculating the net contact force (i.e., the applied side load) in each calculation of the processing loop.
  • the Constant Bit Walk Model historical laboratory drilling data is used to determine the actual walk behavior of a particular bottomhole assembly. For example, historical data may show that a particular BHA will walk 12° to the left.
  • the contact force i.e., the applied side load
  • the contact force in the drillstring steering model of the DDAS may be rotated by the same or substantially the same degree (e.g., 12°) so as to account for bit walk in the simulation calculations.
  • the direction of the contact forces on the drill bit may be rotated, for example, in a
  • the constant angle may be based on properties of a designated subterranean formation and/or on testing for bit walk characteristics in a specific drill bit. In some embodiments, the constant angle may be obtained by calculating the constant angle with an empirical formula. In other embodiments, the constant angle may be obtained from a table of known values, for example, formation properties and/or bit walk characteristics.
  • the constant angle may be obtained by using a value equal to an angle between the contact forces and the frictional forces at the drill bit.
  • the finite element analysis (FEA) computations used in the "force balancing" component of the BHA finite element program/procedure may include calculating the torsional friction between components of the drillstring (e.g., gage pads) and the walls of the wellbore, which modifies the direction and the magnitude of the net contact force (i.e., the applied side load) in each calculation of the processing loop.
  • the torsional frictional forces may be calculated using the formula:
  • F f F n *FricCoefficient, where Ff represents the torsional frictional forces, F n represents the normal contact forces and
  • FricCoefficient represents a frictional factor
  • an individual frictional factor may be assigned to one or more individual components of the BHA. Further, the accounting of the torsional friction in the calculations may also affect displacements such as angular displacement (i.e., tilt angle) of the drill bit, which may further affect bit walk, as described in further detail below with respect to FIG. 45.
  • FIG. 40 is a process flowchart for a DDAS method 400, including the additional models (i.e., the Constant Bit Walk Model and the Torsional Friction Bit Walk Model) to account for and predict bit walk and/or wellbore spiraling.
  • the BHA model may provide forces, drill bit tilt, moments, and curvature information to the axial motion drill bit model and lateral motion drill bit model at 402.
  • the BHA model may also provide hole curvature and BHA alignment within the hole to the steering computer model (i.e., drillstring steering model or steering algorithm). Further, the BHA model may operate in a drill-ahead depth domain within the DDAS.
  • the BHA model may calculate BHA contact forces (F n ) and/or drill bit angular displacements, including tilt. Both the tilt of the drill bit (with respect to the borehole) and an applied side load that forces the drill bit against the borehole wall are accounted for in the SIDECUT procedure.
  • BHA frictional forces (F f ) may be calculated by applying frictional force values to one or more components of the BHA (e.g., gage pads, ribs, steering devices, etc.) either individually or collectively.
  • the BHA model may then apply the calculated frictional forces to the BHA in each iteration of the DDAS using the FEA with error-reducing iterations and may calculate final values for angular displacement, contact forces, and/or frictional forces. In this manner, the torsional frictional forces between the borehole wall and the drill bit and/or the BHA may be calculated and applied to predict BUR and walk rate.
  • the axial motion drill bit model may use inputs from formation lithology data (shown at 412) and drill rig operating parameters (shown at 410). As mentioned above, the axial motion drill bit model may calculate rock strength and ROP, which calculations may then be given to the lateral motion drill bit model.
  • the lateral motion drill bit model may use the rock strength and ROP data from the axial motion drill bit model along with the data supplied by the BHA model, such as bit side force and bit tilt, for example, to calculate the ratio of lateral ROP to axial ROP, or in the alternative, to calculate the ratio of lateral displacement to axial displacement for the drill bit for a certain depth or certain time period, as described in greater detail above with reference to FIG. 4.
  • the steering computer model may include the front steering angle model, or optionally, the rear steering angle model, as explained in greater detail above with reference to FIGS. 5A through 8.
  • the steering computer model may, in some embodiments, rotate the contact force at the drill bit by a constant walk angle.
  • the steering computer model may account for frictional forces by applying the lateral motion drill bit model side-cutting results in the direction of the rotated contact force at the drill bit.
  • the steering computer model may then use the updated models to calculate inclination and azimuth for the end of a new hole section to be created.
  • a new hole section may be created based on the calculated inclination and azimuth.
  • the drillstring i.e., BHA and drill bit
  • the three models BHA model, axial motion drill bit model and lateral motion drill bit model
  • the cycle continues until the simulation is stopped.
  • the drill- ahead simulation method runs the BHA model first, then feeds its results to the axial motion and lateral motion drill bit models, which models may be incorporated into the DDAS method.
  • the DDAS method then predicts the location and geometry of the next hole section.
  • the steering computer model may be supplied information by the other models, including the Constant Bit Walk Model and the Torsional Friction Bit Walk Model, and a new simulated hole section may be placed in the direction that the bit BHA system wants to drill as calculated by the steering computer model. It is to be appreciated that initial conditions and past behavior along with forces on the drill bit and/or BHA affect the future behavior of the drill bit trajectory.
  • any combination of model or models for accounting for and predicting bit walk and/or wellbore spiraling may be used by the steering computer model.
  • the Constant Bit Walk Model and the Torsional Friction Bit Walk Model may be applied separately or in combination.
  • the original DDAS method i.e., not accounting for bit walk and/or wellbore spiraling
  • the torsional frictional forces (F f ) may be applied to the BHA model in order to supplement the original DDAS method.
  • a user may input a constant angle to rotate the contact force at the drill bit by the predicted walk angle.
  • a computer program may calculate a walk angle with, for example, an empirical formula or a table of known values for formation properties and/or drill bit configurations.
  • the walk angle may be updated automatically within the computer program according to a defined function of certain parameters of the BHA. For example, gage pad length, lateral force on the drill bit, friction factor, ROP, RMP, etc., may be factored into an equation, either separately or in combination, for calculating the bit walk angle.
  • the walk angle may be calculated by assuming the walk angle is equal to the angle between the contact force and frictional force at the drill bit. It is to be appreciated that any type of user inputted or calculated value may be used to adjust the angle of the contact force.
  • FIG. 41 is a simplified drawing of a downhole view 500 of the borehole 2 shown in
  • FIG. 1 An end view of the drill tubular 5 followed by the drill bit 7 is depicted from looking downhole. As illustrated, the drill bit 7 is tilted in an arbitrary direction due to a lateral contact force 504 in an arbitrary direction.
  • the lateral contact force 504 calculated and employed in the BHA finite element program/procedure of the DDAS may be rotated by a constant angle ⁇ to a rotated contact force 506 in each iteration of the processing loop.
  • the value of the constant angle ⁇ may be selected empirically based on laboratory or field data of actual bit walk angles, such that the results of the simulation align with the historical data for the same or substantially similar BHA designs.
  • FIG. 1 An end view of the drill tubular 5 followed by the drill bit 7 is depicted from looking downhole. As illustrated, the drill bit 7 is tilted in an arbitrary direction due to a lateral contact force 504 in an arbitrary direction.
  • the rotated contact force 506 is shown to be rotated counterclockwise of the lateral contact force 504.
  • the change in the rotation of the lateral contact force 504 may be the only programming change made to the DDAS.
  • the constant angle ⁇ may be between about 10° and about 20°, and more particularly, between about 12° and about 15°. As one non-limiting example embodiment, the constant angle ⁇ may be about 12°.
  • FIGS. 42A through 42K contain a series of graphs illustrating results for the system walk for the entire BHA system using the Constant Bit Walk Model in the DDAS. Results are shown for left-hand bit walk (negative turn rate) while building with the AUTOTRAKTM eXact Rotary Steerable System (RSS), which is commercially available from Baker Hughes Incorporated of Houston, TX.
  • the DDAS may be programmed with a constant 12° bit walk angle.
  • Various results versus measured depth (MD) are illustrated.
  • FIG. 42A illustrates inclination versus MD
  • FIG. 42B illustrates azimuth versus MD
  • FIG. 42C illustrates build-up-rate (BUR) versus MD
  • FIG. 42D illustrates rum rate versus MD
  • FIG. 42A illustrates inclination versus MD
  • FIG. 42B illustrates azimuth versus MD
  • FIG. 42C illustrates build-up-rate (BUR) versus MD
  • FIG. 42D illustrates rum rate versus MD
  • FIG. 42E illustrates dogleg severity (DLS) versus MD
  • FIG. 42F illustrates first contact force versus MD.
  • Various parameters versus MD are illustrated.
  • FIG. 42G illustrates rate of penetration (ROP) versus MD
  • FIG. 42H illustrates drill bit rotary speed (RPM) versus MD
  • FIG. 421 illustrates weight-on-bit (WOB) versus MD
  • FIG. 42J illustrates rib force versus MD
  • FIG. 42K illustrates rib orientation versus MD.
  • FIGS. 43 A through 43K contain a series of graphs illustrating results for the system walk rate for the entire BHA system using the Constant Bit Walk Model in the DDAS.
  • Results are shown for right-hand bit walk (positive rum rate) while dropping with the AUTOTRAKTM eXact RSS.
  • the DDAS may also be programmed with a constant 12° bit walk.
  • FIG. 43 A illustrates inclination versus MD
  • FIG. 43B illustrates azimuth versus MD
  • FIG. 43C illustrates BUR versus MD
  • FIG. 43D illustrates turn rate versus MD
  • FIG. 43E illustrates DLS versus MD
  • FIG. 43F illustrates first contact force versus MD.
  • FIG. 43G illustrates ROP versus MD
  • FIG. 43H illustrates RPM versus MD
  • FIG. 431 illustrates WOB versus MD
  • FIG. 44A is a simplified drawing of a top view of a wellbore spiral 600
  • FIG. 44B is a simplified drawing of a side view of the wellbore spiral 600 of FIG. 44A.
  • a spiraled (i.e., corkscrew) hole is a well-known phenomenon in the industry and may be measured by lateral distance 602 (side movement) per measured depth 604, as shown in FIG. 44B.
  • a right- handed spiral 606 is illustrated in FIG. 44B.
  • Wellbore spiraling can usually be detected from wellbore logs.
  • attempts have been made to predict and correct for wellbore spiraling which is more pronounced in horizontal legs and in laterally aggressive systems.
  • the Constant Bit Walk Model may be utilized to predict wellbore spiraling.
  • a 12° correction factor may be implemented in the Constant Bit Walk Model.
  • modeling may be used to predict wellbore spiraling behavior in addition to bit walk behavior. In the lab, spiraling behavior was observed in lateral simulations with the AUTOTRAKTM eXact RSS when a rib force was set to 0 percent.
  • a right-handed thread spiral wellbore was primarily observed. Once spiraling behavior is predicted through modeling, the results may be used to account for the predicted behavior and to take corrective action. For example, longer-gauge bits or bits that are less aggressive laterally may tend to help reduce wellbore spiraling.
  • the computations used in the "force balancing" component of the BHA model may include calculation of the torsional friction between components of the drillstring and the walls of the wellbore, which modifies the direction and the magnitude of the net contact force (i.e., the applied side load) in each calculation of the processing loop.
  • dynamic modeling software including BHASYS, BHASYS PRO AND BHASYS TD software of Baker Hughes Incorporated
  • FIG. 45 is a simplified transverse cross-sectional view 700 of a BHA in a wellbore, and illustrates force vectors acting on the BHA.
  • the drill tubular 5 and drill bit 7 are disposed in a horizontal leg of the borehole 2.
  • Longitudinal axis 702 represents the vertical leg of the borehole 2, prior to building to the horizontal leg.
  • the normal force vector 704 represents the direction and magnitude of normal contact forces on the drill bit 7 as calculated in the previous DDAS software.
  • a torsional friction force vector 706 is calculated and summed with the normal force vector 704, resulting in a resultant vector 708 having a direction and magnitude different from the normal force vector 704.
  • the torsional friction force vector 706 represents the direction and magnitude of frictional forces that are tangential to the borehole 2 and perpendicular to the normal contact forces and represents frictional forces that occur between components of the BHA and the wall of the borehole 2.
  • the drill bit 7 may roll to the left, which is illustrated as the normal force vector 704 being displaced from vertical in FIG. 45.
  • the torsional friction forces due to bit rotation may also cause further displacement of the drill bit 7 from alignment with the longitudinal axis 702.
  • the normal force vector 704 and the torsional friction force vector 706 may be added as is commonly done in vector addition.
  • the resultant vector 708 represents the combined forces and may be used in the Torsional Friction Bit Walk Model of the DDAS software, which has been found to accurately predict and account for bit walk and/or wellbore spiraling.
  • the tangential friction forces are calculated using a fiction factor.
  • the drill bit 7 may also exhibit tilt as a result of the forces, which, in rum, may add additional displacement forces. In drilling, and more particularly in direction drilling, it is well known in the art that as more force is applied to a lateral side of the drill bit 7, the faster the drill bit 7 will drill laterally.
  • FIG. 46 contains a series of graphs illustrating results for bit walk using the Torsional Friction Bit Walk Model in the DDAS. Results are shown for left-hand bit walk in the build direction while executing a 90° turn in a lateral direction with an RSS. Various results versus MD are illustrated in FIG. 46, including hole inclination versus MD, azimuth versus MD, BUR versus MD, and turn rate versus MD. Modeling may be provided comparing a friction coefficient to walk rates in building and dropping scenarios. For example, models may include results of walk rates while turning 90° and 270° rib orientation. Alternatively, models may include results of walk rates while building and dropping 0° and 180° rib orientation.
  • results may show that a slight build is natural for a certain BHA when turning left or right.
  • the results of this scenario may show a 1.5° build for a left or a right rum, which may be due to gravity impacting a lateral side of the drill bit 7 and/or due to tilt in the build direction.
  • a 90° right-turning BHA may walk left where it is being pushed, for example, resulting in a greater build.
  • a 270° left-turning BHA may walk "less right” resulting in less build and may eventually exhibit a "true” walk left at a friction coefficient of about 0.25, for example.
  • adjustments and/or modifications may be made to the wellbore plan, the BHA configuration, tool components, etc., so as to ensure that the results of the drilling process match the wellbore plan.
  • a drilling services company may employ the DDAS as disclosed herein to make recommendations to drilling operators to modify the design of the BHA, drill bit, drilling parameters or the wellbore plan to improve the directional drilling process.
  • selecting or adjusting parameters may include adjusting weight-on-bit, torque applied to a drillstring of the BHA, rotary speed of the drillstring, rate of penetration, or drilling fluid flow rate.
  • adjusting lateral aggressiveness of the drill bit may include, by way of non-limiting example, adjusting gauge of the drill bit, blade configuration, gage pad length, gage pad location, or cutting element placement on the drill bit.
  • the drill bit trajectory may be adjusted during drilling with a rotary steering drilling system connected to the BHA, which may be configured, for example, to adjust the drill bit trajectory based at least in part on calculations from the BHA model.
  • receiving parameters for operating the BHA may include receiving drilling parameters specific to the BHA and receiving lithology data may include receiving lithology data specific to a defined target in the subterranean formation.
  • the drilling parameters and lithology data may be received either prior to drilling and/or in real time during drilling operations.
  • adjustments and/or modifications to the wellbore plan and the BHA system may be made prior to drilling and/or in real time during drilling operations.
  • Embodiment 1 A method of controlling drill bit trajectory in a subterranean formation, the method comprising: receiving drilling parameters for operating a specific bottomhole assembly (BHA); constructing, with a computer processor, a directional drill- ahead simulator comprising a computer model of the BHA and the subterranean formation; calculating, with the computer processor, axial motion and lateral motion of a drill bit connected to a bottom end of the BHA using at least one formation parameter and at least one drilling parameter; predicting, with the computer processor, bit walk of the drill bit by accounting for and calculating contact forces and frictional forces between the BHA and a wall of a borehole in the subterranean formation using the computer model of the BHA; determining, with the computer processor, an adjusted drill bit trajectory to account for the predicted bit walk; determining adjusted drilling parameters for operating the BHA to substantially follow the adjusted drill bit trajectory; and operating the BHA according to the adjusted drilling parameters.
  • BHA specific bottomhole assembly
  • Embodiment 2 The method of Embodiment 1, further comprising predicting wellbore spiraling by accounting for and calculating the contact forces and the frictional forces using the computer model of the BHA.
  • Embodiment 3 The method of Embodiment 1, wherein constructing the directional drill-ahead simulator comprising the BHA computer model comprises: constructing a three- dimensional model of the BHA and the subterranean formation; configuring the BHA computer model to be a dynamic model in a depth domain; configuring the BHA computer model to predict a build-up-rate based at least in part on the predicted bit walk of the drill bit using error-reducing iterations in each step of a finite element analysis; and configuring the BHA computer model to operate in real time during drilling operations.
  • Embodiment 4 The method of Embodiment 1, wherein: calculating the axial motion and the lateral motion of the drill bit comprises using a lateral motion drill bit model based at least in part on the axial motion of the drill bit calculated using an axial motion drill bit model separate from the lateral motion drill bit model and the computer model of the BHA; and determining the adjusted drill bit trajectory comprises using a steering computer model separate from the lateral motion drill bit model and the axial motion drill bit model, wherein the steering computer model is configured to apply results from the lateral motion drill bit model adjusted by the predictions of the computer model of the BHA.
  • Embodiment 5 The method of Embodiment 1, wherein determining the adjusted drill bit trajectory further comprises: calculating a new inclination angle and a new azimuth angle of the drill bit; adding a distance to the borehole in the subterranean formation; moving the drill bit in a drilling direction to an end of the borehole; and iterating the calculating, the adding, and the moving to update the computer model of the BHA.
  • Embodiment 6 The method of Embodiment 1 , wherein predicting the bit walk of the drill bit by accounting for and calculating the contact forces and the frictional forces between the BHA and the wall of the borehole further comprises accounting for angular displacement of the drill bit.
  • Embodiment 7 The method of Embodiment 1, wherein predicting the bit walk of the drill bit by accounting for and calculating the contact forces and the frictional forces between the BHA and the wall of the borehole further comprises calculating torsional frictional forces between at least one component of the BHA and the wall of the borehole.
  • Embodiment 8 The method of Embodiment 1, wherein determining the adjusted drill bit trajectory comprises applying calculations of the lateral motion of the drill bit in a direction of the contact forces and rotating the direction of the contact forces on the drill bit by a constant angle.
  • Embodiment 9 The method of Embodiment 8, wherein rotating the direction of the contact forces on the drill bit by a constant angle further comprises obtaining the constant angle by at least one of calculating the constant angle with an empirical formula, obtaining the constant angle from a table of known values, or using a value equal to an angle between the contact forces and the frictional forces at the drill bit.
  • Embodiment 10 The method of Embodiment 8, wherein rotating the direction of the contact forces on the drill bit by a constant angle comprises rotating the direction of the contact forces on the drill bit by about 12° in a counterclockwise direction from a perspective looking down the borehole.
  • Embodiment 11 The method of Embodiment 1, wherein determining the adjusted drilling parameters for operating the BHA comprises: adjusting lateral aggressiveness of the drill bit by adjusting at least one of gauge of the drill bit, blade configuration, gage pad length, gage pad location, or cutting element placement on the drill bit; adjusting at least one of weight-on-bit, torque applied to a drillstring of the BHA, rotary speed of the drillstring, rate of penetration, or drilling fluid flow rate; and adjusting the drill bit trajectory during drilling with a rotary steering drilling system operably connected to the BHA.
  • Embodiment 12 A method of planning and drilling a wellbore in a subterranean formation, comprising: defining a target in a designated subterranean formation; predicting wellbore spiraling and bit walk of a drill bit connected to a specific bottomhole assembly (BHA), comprising: using a computer processor programmed to execute a directional drill- ahead simulator comprising a computer model of the BHA and the designated subterranean formation; receiving with the computer processor lithology data and drilling parameters for operating the BHA in the designated subterranean formation; calculating with the computer processor a ratio of lateral motion to axial motion using a lateral motion drill bit computer model and an axial motion drill bit computer model; predicting drill bit trajectory with the computer processor by accounting for and calculating lateral contact forces, angular displacement, and frictional forces using the computer model of the BHA; and adjusting the drill bit trajectory with the computer processor based at least in part on predictions from the computer model of the BHA; adjusting the drilling parameters for operating the BHA to substantially
  • Embodiment 13 The method of Embodiment 12, further comprising updating the drilling parameters during drilling using information received from at least one sensor connected to the BHA.
  • Embodiment 14 The method of Embodiment 12, wherein predicting the drill bit trajectory comprises calculating the frictional forces between at least one component of the BHA and a borehole wall and adding the frictional forces to the lateral contact forces.
  • Embodiment 15 The method of Embodiment 12, wherein adjusting the drill bit traj ectory comprises rotating a direction of a lateral contact force on the drill bit by a constant angle, wherein the constant angle is based on properties of the subterranean formation and on testing for bit walk characteristics in a specific drill bit.
  • Embodiment 16 The method of Embodiment 12, wherein adjusting the drilling parameters for operating the BHA further comprises: selecting components of the BHA based at least in part on the predictions for the wellbore spiraling and the bit walk of the drill bit; adjusting at least one of weight-on-bit, torque applied to a drillstring of the BHA, rotary speed of the drillstring, rate of penetration, or drilling fluid flow rate; and adjusting lateral aggressiveness of the drill bit by adjusting at least one of gauge of the drill bit, blade configuration, gage pad length, gage pad location, or cutting element placement on the drill bit.
  • Embodiment 17 The method of Embodiment 12, wherein adjusting the drilling parameters for operating the BHA comprises adjusting the drilling parameters using a rotary steering drilling system.
  • Embodiment 18 A method of controlling drill bit trajectory in a subterranean formation, the method comprising: receiving lithology data for a specific subterranean formation; receiving one or more drilling parameters for operating a bottomhole assembly (BHA), the one or more drilling parameters comprising at least one of weight-on-bit, torque, rotary speed, rate of penetration, drilling fluid flow rate, or lateral aggressiveness of a drill bit; predicting wellbore spiraling and bit walk of a drill bit of the BHA, comprising: constructing, with a computer processor, a directional drill-ahead simulator comprising a dynamic computer model of the BHA and the subterranean formation; accounting for the wellbore spiraling and the bit walk of the drill bit by rotating a direction of normal contact forces on the drill bit by a constant angle using the computer processor; calculating combined forces on the BHA
  • predicting the drill bit trajectory with the computer processor, based at least in part on calculating the combined forces on the BHA; and adjusting the drill bit trajectory, with the computer processor, based at least in part on predictions for the drill bit trajectory; and adjusting the one or more drilling parameters based at least in part on predictions for the wellbore spiraling and the bit walk of the drill bit.
  • Embodiment 19 The method of Embodiment 18, wherein calculating the combined forces on the BHA comprises calculating the torsional frictional forces using:
  • Ff F n *FricCoefficient, wherein Ff represents the torsional frictional forces, F n represents the normal contact forces and FricCoefficient represents a frictional factor; and wherein an individual frictional factor may be assigned to one or more individual components of the BHA.
  • Embodiment 20 The method of Embodiment 18, wherein adjusting the drill bit trajectory comprises adjusting the drill bit trajectory with a rotary steering drilling system operably connected to the BHA, the rotary steering drilling system configured to adjust the drill bit trajectory based at least in part on calculations from the dynamic computer model of the BHA.

Landscapes

  • Engineering & Computer Science (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Mining & Mineral Resources (AREA)
  • Geology (AREA)
  • Physics & Mathematics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Environmental & Geological Engineering (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Fluid Mechanics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • General Engineering & Computer Science (AREA)
  • Geometry (AREA)
  • Evolutionary Computation (AREA)
  • Computer Hardware Design (AREA)
  • Geophysics (AREA)
  • Mechanical Engineering (AREA)
  • Earth Drilling (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)
  • Geophysics And Detection Of Objects (AREA)

Abstract

La présente invention concerne un procédé de commande de trajectoire de trépan dans une formation souterraine. Ledit procédé comprend la réception de paramètres de forage pour actionner un ensemble en profondeur de forage (BHA) spécifique, la construction, avec un processeur informatique, d'un simulateur de forage directionnel qui comprend un modèle informatique du BHA et de la formation souterraine, le calcul du mouvement axial et du mouvement latéral d'un trépan relié à une extrémité inférieure du BHA en utilisant des paramètres de formation et des paramètres de forage, la prédiction de marche de trépan du trépan en tenant compte et en calculant des forces de contact et des forces de frottement entre le BHA et une paroi d'un trou de forage dans la formation souterraine en utilisant le modèle informatique du BHA, et la détermination d'une trajectoire de trépan ajustée pour tenir compte de la marche de trépan prédite. Le procédé comprend la détermination de paramètres de forage ajustés pour actionner le BHA pour suivre sensiblement la trajectoire de trépan ajustée et actionner le BHA selon les paramètres de forage ajustés.
PCT/US2017/042876 2016-07-20 2017-07-19 Procédés de commande de trajectoire de trépan par prédiction de marche de trépan et forage en spirale WO2018017724A1 (fr)

Priority Applications (6)

Application Number Priority Date Filing Date Title
CA3031465A CA3031465C (fr) 2016-07-20 2017-07-19 Procedes de commande de trajectoire de trepan par prediction de marche de trepan et forage en spirale
GB1901944.7A GB2567376B (en) 2016-07-20 2017-07-19 Methods of controlling drill bit trajectory by predicting bit walk and wellbore spiraling
CN201780054409.3A CN109661502B (zh) 2016-07-20 2017-07-19 通过预测钻头游动和井筒螺旋来控制钻头轨迹的方法
BR112019001148-0A BR112019001148B1 (pt) 2016-07-20 2017-07-19 Método para controlar trajetória de broca de perfuração em uma formação subterrânea
SA519400930A SA519400930B1 (ar) 2016-07-20 2019-01-19 طرق للتحكم في مسار لقمة حفر من خلال التنبؤ بمسار لقمة الحفر وحلزونية حفرة البئر
NO20190197A NO20190197A1 (en) 2016-07-20 2019-02-12 Methods of controlling drill bit trajectory by predicting bit walk and wellbore spiraling

Applications Claiming Priority (4)

Application Number Priority Date Filing Date Title
US201662364833P 2016-07-20 2016-07-20
US62/364,833 2016-07-20
US15/348,003 US10296678B2 (en) 2013-10-18 2016-11-10 Methods of controlling drill bit trajectory by predicting bit walk and wellbore spiraling
US15/348,003 2016-11-10

Publications (1)

Publication Number Publication Date
WO2018017724A1 true WO2018017724A1 (fr) 2018-01-25

Family

ID=60992854

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/US2017/042876 WO2018017724A1 (fr) 2016-07-20 2017-07-19 Procédés de commande de trajectoire de trépan par prédiction de marche de trépan et forage en spirale

Country Status (7)

Country Link
CN (1) CN109661502B (fr)
BR (1) BR112019001148B1 (fr)
CA (1) CA3031465C (fr)
GB (1) GB2567376B (fr)
NO (1) NO20190197A1 (fr)
SA (1) SA519400930B1 (fr)
WO (1) WO2018017724A1 (fr)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109472058A (zh) * 2018-10-16 2019-03-15 成都泰隆游乐实业有限公司 一种水滑梯轨迹分析方法
CN112052586A (zh) * 2020-09-02 2020-12-08 四川省建筑科学研究院有限公司 一种顶力预测方法及顶力预测装置
US20230203935A1 (en) * 2021-12-29 2023-06-29 Halliburton Energy Services, Inc. Method for real-time pad force estimation in rotary steerable system
EP4034749A4 (fr) * 2019-09-26 2023-11-15 Services Pétroliers Schlumberger Dispositif de commande d'exploitation de forage par vitesse de pénétration

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112412435A (zh) * 2020-10-23 2021-02-26 中国石油集团工程咨询有限责任公司 一种用于螺旋井眼的密度曲线失真校正方法
CN113944456B (zh) * 2021-09-26 2024-04-16 浙江省工程勘察设计院集团有限公司 一种钻进深度测量方法、系统、钻机及存储介质
CN117514128B (zh) * 2024-01-08 2024-03-19 山西郎腾信息科技有限公司 一种煤矿井下水平定向智能化钻探控制方法
CN118350137B (zh) * 2024-03-11 2025-01-24 西南石油大学 一种螺旋刀翼pdc钻头切削结构的设计方法及系统

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4804051A (en) * 1987-09-25 1989-02-14 Nl Industries, Inc. Method of predicting and controlling the drilling trajectory in directional wells
US6438495B1 (en) * 2000-05-26 2002-08-20 Schlumberger Technology Corporation Method for predicting the directional tendency of a drilling assembly in real-time
US6443242B1 (en) * 2000-09-29 2002-09-03 Ctes, L.C. Method for wellbore operations using calculated wellbore parameters in real time
US20100259415A1 (en) * 2007-11-30 2010-10-14 Michael Strachan Method and System for Predicting Performance of a Drilling System Having Multiple Cutting Structures
US20130043077A1 (en) * 2005-08-08 2013-02-21 Shilin Chen Methods and systems for designing and/or selecting drilling equipment using predictions of rotary drill bit walk

Family Cites Families (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN1282818C (zh) * 2001-08-16 2006-11-01 中海油田服务股份有限公司 水平井钻头前进方向的预测方法、控制方法及其控制系统
EP2118441B1 (fr) * 2007-01-08 2016-08-10 Baker Hughes Incorporated Composants de forage et systèmes pour contrôler de manière dynamique des dysfonctionnements en termes de forage et procédés de forage d'un puits avec ceux-ci
US7957946B2 (en) * 2007-06-29 2011-06-07 Schlumberger Technology Corporation Method of automatically controlling the trajectory of a drilled well
EP2179134B1 (fr) * 2007-08-06 2018-01-10 GeoMechanics International Système et procédé de guidage de puits de forage à base de champ de champ de contrainte
WO2009079371A1 (fr) * 2007-12-14 2009-06-25 Halliburton Energy Services, Inc. Procédés et systèmes de prédiction du pas d'un trépan rotatif et de conception de trépans rotatifs et d'autres outils de fond de trou
US20100185395A1 (en) * 2009-01-22 2010-07-22 Pirovolou Dimitiros K Selecting optimal wellbore trajectory while drilling
US20120132468A1 (en) * 2010-11-30 2012-05-31 Baker Hughes Incorporated Cutter with diamond sensors for acquiring information relating to an earth-boring drilling tool
US10132119B2 (en) * 2013-10-18 2018-11-20 Baker Hughes, A Ge Company, Llc Directional drill ahead simulator: directional wellbore prediction using BHA and bit models
CN105003245B (zh) * 2015-07-13 2017-10-03 中国海洋石油总公司 一种井下定向动力钻具工具面的动态控制系统及方法

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4804051A (en) * 1987-09-25 1989-02-14 Nl Industries, Inc. Method of predicting and controlling the drilling trajectory in directional wells
US6438495B1 (en) * 2000-05-26 2002-08-20 Schlumberger Technology Corporation Method for predicting the directional tendency of a drilling assembly in real-time
US6443242B1 (en) * 2000-09-29 2002-09-03 Ctes, L.C. Method for wellbore operations using calculated wellbore parameters in real time
US20130043077A1 (en) * 2005-08-08 2013-02-21 Shilin Chen Methods and systems for designing and/or selecting drilling equipment using predictions of rotary drill bit walk
US20100259415A1 (en) * 2007-11-30 2010-10-14 Michael Strachan Method and System for Predicting Performance of a Drilling System Having Multiple Cutting Structures

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109472058A (zh) * 2018-10-16 2019-03-15 成都泰隆游乐实业有限公司 一种水滑梯轨迹分析方法
EP4034749A4 (fr) * 2019-09-26 2023-11-15 Services Pétroliers Schlumberger Dispositif de commande d'exploitation de forage par vitesse de pénétration
CN112052586A (zh) * 2020-09-02 2020-12-08 四川省建筑科学研究院有限公司 一种顶力预测方法及顶力预测装置
US20230203935A1 (en) * 2021-12-29 2023-06-29 Halliburton Energy Services, Inc. Method for real-time pad force estimation in rotary steerable system
US11788400B2 (en) * 2021-12-29 2023-10-17 Halliburton Energy Service, Inc. Method for real-time pad force estimation in rotary steerable system

Also Published As

Publication number Publication date
SA519400930B1 (ar) 2023-01-02
BR112019001148A2 (pt) 2019-04-30
GB201901944D0 (en) 2019-04-03
BR112019001148B1 (pt) 2023-04-11
NO20190197A1 (en) 2019-02-12
CN109661502B (zh) 2022-09-23
GB2567376A (en) 2019-04-10
GB2567376B (en) 2021-12-01
CA3031465A1 (fr) 2018-01-25
CA3031465C (fr) 2021-03-02
CN109661502A (zh) 2019-04-19
BR112019001148A8 (pt) 2023-02-14

Similar Documents

Publication Publication Date Title
US10296678B2 (en) Methods of controlling drill bit trajectory by predicting bit walk and wellbore spiraling
CA3031465C (fr) Procedes de commande de trajectoire de trepan par prediction de marche de trepan et forage en spirale
US10132119B2 (en) Directional drill ahead simulator: directional wellbore prediction using BHA and bit models
US9951560B2 (en) Axial motion drill bit model
US8775145B2 (en) System and apparatus for modeling the behavior of a drilling assembly
US10282497B2 (en) Model for estimating drilling tool wear
CA2964228C (fr) Procedes et systemes pour la modelisation d'un ensemble de fond de trou tridimensionnel de pointe
US10012025B2 (en) Lateral motion drill bit model
GB2431753A (en) A method of designing a drill bit for specific formation segments
US20150152689A1 (en) Rotary drill bits with back-up cutting elements to optimize bit life
EP2521830A1 (fr) Système et procédé de modèle d'interaction alésoir et foret
US20200011138A1 (en) Identification of weak zones in rotary drill bits during off-center rotation
US10781642B2 (en) Rotary drill bit including multi-layer cutting elements
Tortrakul et al. Eccentric String Reamer for ECD Reduction in the Depleted Field of Gulf of Thailand GoT
Jardine et al. Lessons Learned from Drilling Horizontal Wells on Didon Field; Offshore Tunisia
CA2566640A1 (fr) Optimisation de la priorisation des formations

Legal Events

Date Code Title Description
121 Ep: the epo has been informed by wipo that ep was designated in this application

Ref document number: 17831797

Country of ref document: EP

Kind code of ref document: A1

ENP Entry into the national phase

Ref document number: 3031465

Country of ref document: CA

NENP Non-entry into the national phase

Ref country code: DE

REG Reference to national code

Ref country code: BR

Ref legal event code: B01A

Ref document number: 112019001148

Country of ref document: BR

ENP Entry into the national phase

Ref document number: 201901944

Country of ref document: GB

Kind code of ref document: A

Free format text: PCT FILING DATE = 20170719

ENP Entry into the national phase

Ref document number: 112019001148

Country of ref document: BR

Kind code of ref document: A2

Effective date: 20190121

122 Ep: pct application non-entry in european phase

Ref document number: 17831797

Country of ref document: EP

Kind code of ref document: A1

点击 这是indexloc提供的php浏览器服务，不要输入任何密码和下载