+

WO2014066981A1 - Procédés et systèmes pour opérations de forage améliorées utilisant des données de forage en temps réel et historiques - Google Patents

Procédés et systèmes pour opérations de forage améliorées utilisant des données de forage en temps réel et historiques Download PDF

Info

Publication number
WO2014066981A1
WO2014066981A1 PCT/CA2013/000919 CA2013000919W WO2014066981A1 WO 2014066981 A1 WO2014066981 A1 WO 2014066981A1 CA 2013000919 W CA2013000919 W CA 2013000919W WO 2014066981 A1 WO2014066981 A1 WO 2014066981A1
Authority
WO
WIPO (PCT)
Prior art keywords
drilling
data
real
time
bit
Prior art date
Application number
PCT/CA2013/000919
Other languages
English (en)
Inventor
Trent Marx
Gary William Reid
Henry Leung
Xiaoxiang Liu
Original Assignee
Resource Energy Solutions Inc.
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
Application filed by Resource Energy Solutions Inc. filed Critical Resource Energy Solutions Inc.
Priority to MX2015005568A priority Critical patent/MX360775B/es
Publication of WO2014066981A1 publication Critical patent/WO2014066981A1/fr

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

Definitions

  • TITLE METHODS AND SYSTEMS FOR IMPROVED
  • the present invention relates to drilling operations, and more particularly, to methods and systems for improved drilling operations through the use of real-time and historical drilling data.
  • drilling operations rely on "after-the-fact" analysis to determine lithology, as well as other parameters.
  • drill bit cuttings are physically analyzed at the surface by the well site geologist, or to be certain of a formations' lithology, physical parameters and geological facies, 30 feet core sections or sidewall cores are taken and analyzed at the surface.
  • samples must be physically removed from the formation and returned to the surface for evaluation.
  • in-situ tools are used.
  • the most popular is spectral analysis of the formation being drilled.
  • Spectral analysis involves the interpretation of spectra obtained from the formation drilled using logging tools including a passive gamma ray detector and a neutron induced gamma ray log.
  • a neutron source is placed alongside the formation and periodically emits bursts of high energy neutrons to excite the atoms in the formation.
  • a detector records the number of counts of returning gamma rays and segregates them according to their energies.
  • the major drawback of this method is the recorded neutron induced gamma ray spectrums which are contaminated by significant background noise due to Compton scattering.
  • Measurement-while-drilling is a type of well logging that incorporates downhole tools providing real-time information to help with steering the bit. These tools typically include sensors for measuring downhole temperature and pressure, azimuth and inclination, drilling mechanics information (e.g. torque, weight-on-bit, rotary speed, etc.) and a resistivity to determine the presence of hydrocarbons and water.
  • the present application is directed generally to methods and systems for improved drilling operations using real-time and historical drilling data. This may include using measurements obtained using measurement-while-drilling technology to make a real-time bit wear prediction, lithology prediction, pore pressure estimation, rotating friction estimation, permeability estimation, and cost estimation. These predictions may be used to optimize weight on bit and bit rotation speed in an aim to obtain a maximized drilling rate and a minimized drilling cost. An expert decision engine may also be provided to guide the operator or user to avoid safety concerns while drilling.
  • a method of drilling a formation comprising: acquiring real-time drilling data; acquiring a bit wear equation; determining a real-time bit wear prediction by using the real-time drilling data to predict a bit efficiency factor and to detect changes in the bit efficiency factor over time; determining a real-time lithology prediction by processing the real-time drilling data through a multilayer neural network; and adjusting the drilling parameters for the drilling operation in real-time based on the real-time bit wear prediction and the real-time lithology prediction.
  • a drilling system including a drilling fusion prediction engine for predicting bit wear and lithology in real-time for use in determining drilling parameters in a drilling operation, the drilling system comprising: one or more sensors for sensing data at a rig during a drilling operation; a computer system including a drilling fusion prediction engine for predicting bit wear and lithology in real-time, the drilling fusion prediction engine configured to receive the sensed data, the drilling fusion prediction engine including a software module for predicting lithology in real-time by processing the sensed data through a multilayer neural network, the drilling fusion prediction engine including a software module for determining a real-time bit wear prediction by using the real-time drilling data to predict a bit efficiency factor and to detect changes in the bit efficiency factor over time; and a communication network for transmitting the sensed data to the drilling fusion prediction engine.
  • a computer program product comprising: a storage medium configured to store computer readable instructions;the computer readable instructions including instructions for, acquiring real-time drilling data; acquiring cost data and a bit wear equation; determining a real-time bit wear prediction by using the real-time drilling data to predict a bit efficiency factor and to detect changes in the bit efficiency factor over time; determining a real-time lithology prediction by processing the real-time drilling data through a multilayer neural network; and adjusting the drilling parameters for the drilling operation in real-time based on the real-time bit wear prediction and the real-time lithology prediction.
  • FIG. la shows the system architecture of a drilling system employing the drilling fusion software in a standalone arrangement according to an embodiment of the present invention
  • FIG. lb shows the system architecture of a drilling system employing the drilling fusion software in a client-server arrangement according to an embodiment of the present invention
  • FIG. 2a is a flow diagram of the overall architecture of the drilling fusion engine according to an embodiment of the present invention
  • FIG. 2b shows the calculated values, estimated values, predicted values and raw inputs of the prediction engine shown in FIG. 2a according to an embodiment the present invention
  • FIG. 3 is a flow chart of the prediction engine according to an embodiment of the present invention.
  • FIG. 4 is a schematic representation of mode selection used in the prediction engine according to an embodiment of the present invention.
  • FIG. 5 is a block diagram of bit wear prediction according to an embodiment of the present invention.
  • FIG. 6 is a graphical diagram depicting the bit wear prediction results for bit #7 from Well #2 according to an embodiment of the present invention
  • FIG. 7 is a block diagram of a neuron according to an embodiment of the present invention.
  • FIG. 8 is a schematic diagram depicting a multilayer neural network for lithology determination according to an embodiment of the present invention.
  • FIG. 9 is a flow chart of a processing subsystem for real-time lithology prediction according to an embodiment of the present invention.
  • FIG. 10 is a screen shot of the user interface during neural network training for Well #2 according to an embodiment of the present invention.
  • FIG. 1 1a shows an intermediate Sigmalog calculations for the formation pore pressure gradient (Ibm/gal) estimation results for Well #1 within depth 2300m- 3200m according to an embodiment of the present invention
  • FIG. l i b shows an intermediate d-exponent calculations for the formation pore pressure gradient (Ibm/gal) estimation results for Well #1 within depth 2300m- 3200m according to an embodiment of the present invention
  • FIG. 1 1 c shows a formation pore pressure gradient (Ibm/gal) estimation results using d-exponent and Sigmalog approaches for Well #1 within depth 2300m- 3200m according to an embodiment of the present invention
  • FIG. 12a shows a D-exponent pore pressure estimation results for Well #1 according to an embodiment of the present invention
  • FIG. 12b shows a pore pressure results for Well #1 according to an embodiment of the present invention
  • FIG. 12c shows a pore pressure (PP) estimation by linear regression for Well #1 according to an embodiment of the present invention
  • FIG. 13a shows a D-exponent pore pressure estimation results for Well #2 according to an embodiment of the present invention
  • FIG. 13b shows a pore pressure results for Well #2 according to an embodiment of the present invention
  • FIG. 13c shows a pore pressure (PP) estimation by linear regression for Well #2 according to an embodiment of the present invention
  • FIG. 14a shows a D-exponent pore pressure estimation results for Well #3 according to an embodiment of the present invention
  • FIG. 14b shows pore pressure results for Well #3 according to an embodiment of the present invention.
  • FIG. 14c shows a pore pressure (PP) estimation by linear regression for Well #3 according to an embodiment of the present invention
  • FIG. 15a shows a trend fusion results of pore pressure estimation from Well #1 according to an embodiment of the present invention
  • FIG. 15b shows a trend fusion results of pore pressure estimation from Well #2 according to an embodiment of the present invention
  • FIG. 15c shows a trend fusion results of pore pressure estimation from Well #3 according to an embodiment of the present invention
  • FIG. 15d shows an optimized pore pressure estimation obtained from integration of the trend fusion results of pore pressure estimation from Well #1 , Well #2, and Well #3 according to an embodiment of the present invention
  • FIG. 15e shows pressure losses in the circulating system according to an embodiment of the present invention
  • FIG. 16a shows a calculated annular pressure loss (APL) for Well #1 according to an embodiment of the present invention
  • FIG. 16b shows a calculated equivalent circulatory density (ECD) for Well #1 according to an embodiment of the present invention
  • FIG. 16c shows a bottom hole pressure (BHP) for Well # 1 according to an embodiment of the present invention
  • FIG. 16d shows a calculated trip margin for Well #1 according to an embodiment of the present invention
  • FIG. 17a shows a calculated annular pressure loss (APL) for Well #2 according to an embodiment of the present invention
  • FIG. 17b shows a calculated equivalent circulatory density (ECD) for Well #2 according to an embodiment of the present invention
  • FIG. 17c shows a bottom hole pressure (BHP) for Well #2 according to an embodiment of the present invention
  • FIG. 17d shows a calculated trip margin for Well #2 according to an embodiment of the present invention.
  • FIG. 18a shows a calculated annular pressure loss (APL) for Well #3 according to an embodiment of the present invention
  • FIG. 18b shows a calculated equivalent circulatory density (ECD) for Well #3 according to an embodiment of the present invention
  • FIG. 18c shows a bottom hole pressure (BHP) for Well #3 according to an embodiment of the present invention
  • FIG. 18d shows a calculated trip margin for Well #3 according to an embodiment of the present invention.
  • FIG. 19 is a diagram depicting a force balance on a drill string element according to an embodiment of the present invention.
  • FIG. 20 is a diagram depicting measured real-time drilling torque for Well #2 according to an embodiment of the present invention.
  • FIG. 21 is a diagram depicting WOB data for Well #2 according to an embodiment of the present invention.
  • FIG. 22 is a diagram depicting estimated rotating friction factor for Well #2 according to an embodiment of the present invention.
  • FIG. 23 is a diagram depicting measured neutron porosity for an interval of Well #2 according to an embodiment of the present invention.
  • FIG. 24a is a diagram depicting estimated permeability for an interval of Well #2 according to an embodiment of the present invention.
  • FIG. 24b shows a front page summary of the tour sheet for Well #1 according to an embodiment of the present invention
  • FIG. 24c shows a portion of the tour sheet of Well # 1 relating to the drilling assembly according to an embodiment of the present invention
  • FIG. 24d shows a portion of the tour sheet of Well #1 relating to the mud record and circulation according to an embodiment of the present invention
  • FIG. 24e shows a portion of the tour sheet of Well # 1 relating to reduced pump speed according to an embodiment of the present invention
  • FIG. 24f shows calculations displayed to the user for kill sheet #1 according to an embodiment of the present invention
  • FIG. 24g shows calculations displayed to the user for kill sheet #2 according to an embodiment of the present invention
  • FIG. 24h shows calculations displayed to the user for kill sheet #3 according to an embodiment of the present invention
  • FIG. 24i shows calculations displayed to the user for kill sheet #4 according to an embodiment of the present invention
  • FIG. 24j shows a copy of the bit record of Well #1 according to an embodiment of the present invention.
  • FIG. 24k shows a copy of the bit record of Well #2 according to an embodiment of the present invention.
  • FIG 241 shows a copy of the bit record of Well #3 according to an embodiment of the present invention
  • FIG. 25 shows a decision tree of the expert decision engine relating to mud volume according to an embodiment of the present invention
  • FIG. 26 shows a decision tree of the expert decision engine relating to porosity according to an embodiment of the present invention
  • FIG. 27 shows a decision tree of the expert decision engine relating to permeability according to an embodiment of the present invention
  • FIG. 28 shows a decision tree of the expert decision engine relating to the concern of a stuck pipe according to an embodiment of the present invention
  • FIG. 29 shows a decision tree of the expert decision engine relating to bit wear according to an embodiment of the present invention
  • FIG. 30 shows a decision tree of the expert decision engine relating to environmental concerns according to an embodiment of the present invention
  • FIG. 31 shows a graphical user interface (GUI) with three panels according to an embodiment of the present invention
  • FIG. 32 shows a main display of the GUI according to an embodiment of the present invention
  • FIG. 33 shows a MWD display of the GUI in metric units according to an embodiment of the present invention
  • FIG. 34 shows a MWD display of the GUI in imperial units according to an embodiment of the present invention
  • FIG. 35 shows a hydrostatic pressure display of the GUI according to an embodiment of the present invention
  • FIG. 36 shows a dynamic pressure display of the GUI according to an embodiment of the present invention
  • FIG. 37 shows a logging-while-drilling (LWD) display of the GUI according to an embodiment of the present invention
  • FIG. 38 shows a kill sheet display of the GUI according to an embodiment of the present invention
  • FIG. 39 shows an expert decision engine prompt for a mud volume issue for Well #1 according to an embodiment of the present invention
  • FIG. 40 shows an expert decision engine prompt for potential stuck pipe due to mud density and mud volume for Well #1 according to an embodiment of the present invention
  • FIG. 41 shows an expert decision engine prompt for a mud volume issue for Well #2 according to an embodiment of the present invention
  • FIG. 42 shows an expert decision engine prompt for bit wear for Well #2 according to an embodiment of the present invention.
  • FIG. 43 shows an expert decision engine prompt for a porosity issue for Well #2 according to an embodiment of the present invention.
  • Embodiments of the present invention are generally directed to methods and systems for improved drilling operations using real-time and historical drilling data.
  • Embodiments of the invention may be used by any persons or entities which control or conduct drilling operations, such as oil and gas exploration and production companies and oil and gas service companies or mining exploration companies, mining coring companies and mining service companies.
  • the method, system and computer program product when deployed, may identify and predict geological and drilling events ahead of the drill bit, such that the user's drilling system may more precisely steer the wellbore while rotating the drill string or when using the mud motor to increase rate of penetration (ROP) and reduce drilling cost.
  • the present invention may incorporate the steps detailed below for bit wear prediction, lithology prediction, pore pressure estimation, rotating friction coefficient estimation, permeability estimation, and cost estimation using the real-time data collected by a variety of sensors, including logging tools, deployed at the well site.
  • a drilling fusion engine which may integrate drilling data in real-time from a rig site to inform adjustments to drilling parameters to improve the results of drilling operations (also referred to as "optimization").
  • This may include specific techniques that enable real-time drilling optimization: bit wear prediction, lithology prediction, pore pressure estimation, rotating friction coefficient estimation, permeability estimation, and cost estimation.
  • software may be installed and used on site as a client-server application or a standalone computer program.
  • the client-server application involves two subcomponents: 1) the client software with a graphical user interface to interact with the operator onsite; and 2) the server software running the engines described below, which may reside on a server computer onsite or offsite, such as in a remote data center or cloud infrastructure.
  • FIG. la which depicts the system architecture of a drilling system in the standalone arrangement, there is a drilling system 10 comprising a well being drilled 320.
  • Sensors 101 such as through an electronic drilling recorder, such as the Pason(tm) EDR or the NOV M/D Totco(tm), or other logging tools, acquire data and the data is transmitted over a rig site network 103 to the DFA server software 106.
  • historical database 309 may provide historical data obtained from analogous wells 302.1 to 302. n to the DFA server software 106.
  • sensors 101 acquire data, such as through an electronic drilling recorder or other logging tools, which is then communicated over a rig site network 103 to the DFA client software 108 to an operation center network 104 via a satellite 105.
  • the received data is then communicated through the operation center network 104 to the DFA server software 106.
  • the data may be communicated from the sensors 101 to the DFA client software 108 and the DFA server software 106 through various wired or wireless communication networks as appropriate according to other embodiments.
  • historical database 309 may provide historical data obtained from analogous wells 302.1 to 302. n to the DFA server software 106.
  • the drilling fusion engine 301 comprises three main engines: a drilling fusion prediction engine 20, a trend fusion engine 308 and an expert decision engine 328.
  • the prediction engine 20 (which is shown in more detail in FIG. 3) may receive real-time drilling parameters 321 as data inputs 201 from sensors 101 at the well during drilling 320. Using the data inputs 201, the prediction engine 20 may calculate the following values 324: hydrostatic pressure 3241 , ECD 2074, annular circulation pressure loss 3243, BHP 3244, volume 3245, and water saturation 3246.
  • the prediction engine 20 may estimate the following values 325: pore pressure 2053, porosity 2075, permeability 2076, and filter cake 3254.
  • the prediction engine may predict the following values 326: lithology prediction 2057 and bit wear prediction 2056.
  • certain data inputs 201 may be as used as raw inputs 327 such as: ROP (rate of penetration) 2032, torque 2036, WOB (weight on bit) 2033, RPM (revolutions per minute) 3273, pump pressure 3275, mud parameters 3276 and differential pressure 3277.
  • the calculated values 324, estimated values 325, predicted values 326 and raw inputs 327 may be used to determine adjusted real-time drilling parameters 329, which may employed in steering the drill bit during drilling or may be provided to an expert decision engine 328.
  • the expert decision engine 328 may receive the adjusted real-time drilling parameters 329 and as well as analogous well data from a historical database 309, and use this information to generate flags and control operations 330 for well control which may generate prompts to the operator through the graphical user interface.
  • Well control refers to the management of the dangerous effects caused by an unexpected influx of formation fluids (gas/oil/condensate) into the wellbore and hence potentially causing damage to personnel and surface equipment and into the nearby environment and atmosphere. Part of well control includes preventing formation fluid, referred to as kick, from entering into the wellbore during drilling.
  • the drilling fusion engine 301 comprises monitoring a well for signs of impending influx in the formation being drilled in real time, and providing predictive capability for well control by combining real-time drilling data and trends in analogous well data into an expert knowledge system.
  • the expert decision engine 328 may offer decision support towards potential safety issues during drilling and the associated operational suggestions to drilling engineers and geologists on the rig site or, according to an embodiment, be configured to interface with the drilling control system to automatically adjust drilling parameters to avoid well control hazards.
  • the trend fusion engine 308 may generate the analogous well data to populate the historical database 309.
  • the trend fusion engine 308 may optimize the pore pressure estimation results by incorporating historical well data from multiple wells using trend analysis techniques, which may reduce uncertainties due to estimation errors in individual well data.
  • Analogous wells 1 to n 302.1 to 302. n are selected by a user and the data related to the historical drilling parameters and operations 303.1 to 303. n from one or more past drilling jobs for each well are input or imported.
  • the trend fusion engine 308 may then generate expected analogous values for lithology 310, porosity 31 1, permeability 312, water saturation 313, bit wear 314, mud parameters 315, Bottom hole assembly (BHA) 331, pore pressure 316, Rate of Penetration (ROP) 332, Weight of Bit (WOB) 333, Revolutions per minute (RPM) 334, volume (mud tanks and hole volumes) 335, pump 317 and filter cake 318, and other parameters if the parameters are available in the historical data, based on the trends in the historical well data and, according to an embodiment, may also insert the current, real-time data into the historical database 309 for comparison.
  • FIG. 3 shows a flow chart of a prediction engine 20 according to an embodiment of the invention.
  • the data inputs 201 to the prediction engine 20 may be the real-time data received from the sensors 101 (e.g. MWD or logged data) as shown in FIG. 1. These data inputs are pre-processed and mapped by a data preparation and mapping component 201a.
  • the data inputs 201 from sensors 101 may be pre- processed and mapped to the real-time drilling data fields 203 and logged data fields 204 in a database 202, such as a Microsoft SQL Server(tm) database or other object- based database. Since data sources may have different data formats, the data preprocessing and mapping includes converting the data to a unified format prior to importing the data into the database 202.
  • the prediction engine 20 may be implemented by a computer program which may receive various data inputs that may be converted to a unified format using a conversion program. Additionally, a unit conversion table needs to be established. For data mining and information processing tasks dealing with a huge amount of data, data processing may be carried out to prepare data for further processing. Such data processing tasks may include operations to replace missing data with default values, and to correct incorrect data inputs, for example, through time-based interpolation and to address non-registered multiple data sources, for example, through aligning the data with the same sampling rate or representing the same context. These data processing tasks may then output time or depth aligned data.
  • this data may then also be fed into a data smoothing function to remove measurement noise in order to improve lithology prediction accuracy.
  • a data smoothing function is to use a moving average filtering technique.
  • a depth-weighted smoother with a 1 -meter window e.g. the depth point has 0.2 m sampling frequency, 1 meter window corresponds to 5 data points.
  • y (5) (x(l) + x(2) + x(3) + x(4) + x(5)) / 5
  • x(0) denotes the first data point available for smoothing operation. The starting point of averaging for the current output is found in the fifth data ahead in this case.
  • the database 202 may include pre-processed and mapped real-time drilling data 203 and logged data 204.
  • the real-time drilling data 203 may include one or more data inputs relating to depth 2031, Rate of Penetration (ROP) 2032, Weight On Bit (WOB) 2033, rotary speed 2034, flow rate 2035, and torque 2036, as shown in Table 1 , below, according to an embodiment.
  • the logged data 204 may include data inputs from a bit record 2041 , mud log 2042, lithology core database 2043, sonic log 2044, gamma ray log 2045, neutron density log 2046, resistivity log 2047, and micro log 2048.
  • the bit record 2041 may include a record of every bit run during an analogous drilling operation including bit size, IADC name and type, depth in and out, time in the hole, WOB, RPM, bit gauge, bit tooth wear, bearing wear, and total time in the hole, as shown in Tables 3 and 4, below, according to an embodiment.
  • the mud log 2042 may include mud weight, mud viscosity, filter cake thickness, chemical usage, mud properties, information found on a daily drilling report, and solids content, as shown in Table 2, below, according to an embodiment.
  • the lithology core database 2043 may include data relating to the formation geology such as depth, lithology type, permeability, porosity, water saturation, core length and core size, gross and net returns, facies, laboratory generated permeability, and fluid type associated with the formations cored.
  • the sonic log 2044 may include formation lithology as derived from calculations of sonic wave returns, sonic or interval transit times, possible fracture size and orientations.
  • the gamma ray log 2045 may include a measurement of gamma radiation: gamma ray, as shown in Table 5, below, according to an embodiment.
  • the neutron density log 2046 may include measurements of caliper, compensated neutron porosity, compensated neutron porosity dolomite/limestone/sandstone matrix, photo electric cross section, bulk density, temperature, density porosity, and total gas, as shown in Table 5, below, according to an embodiment.
  • the resistivity log 2047 may include measurements of shallow induction, medium induction, deep induction, and spontaneous potential, as shown in Table 5, below, according to an embodiment.
  • the micro log 2048 may include measurements of micro normal and micro inverse, as shown in Table 5, below, according to an embodiment.
  • embodiments of the invention may operate in the absence of some or many of these parameters, if they are not available. Additional data, including other captured or logged data parameters, however, may improve the accuracy of the predictions made by the prediction engine 20 and the recommendations of the expert system engine 328.
  • the real-time drilling data 203 and the logged data 204 may be used to calculate the drilling bits hydraulics 2051, the drag 2052 and torque 2073, the ECD 2074 & pore pressure 2053, the porosity 2075, permeability 2076 and skin factors 2054, rock compressive strength 2055, bit wear prediction 2056 and lithology prediction 2057.
  • the mud log 2042 and bit record 2041 may be used in porosity 2075 and permeability 2076 estimation, and skin factor estimation 2054.
  • the real-time drilling data 203 and logged data 204 may be used by the drilling optimization engine 205 to generate optimized drilling parameters, i.e. those adjustments made to current drilling parameters intended to improve drilling results.
  • supporting data 206 which may comprise bit dull grade 2063, cost information 2061, i.e. the cost of the bit, and a bit wear equation 2062 may be used to assist with optimizing drilling parameters with respect to bit wear replacement.
  • Manufacturers of drilling software such as Wellman(tm) may provide supporting data 206 such as bit dull grade for each bit run and cost information.
  • Bit manufacturer data may include equations and parameters to determine the bit tooth wear for the specific bit used in drilling.
  • the supporting data 206 from the manufacturer may be imported into the computer program. The supporting data may then be used for the online optimization 207 of bit wear prediction and drilling costs analysis for generating the optimized drilling parameters 2072 and cost analysis 2071.
  • the lithology prediction 2057 may use neural networks to construct a non-linear map between the mechanical drilling data and the lithology types.
  • the output of lithology types may be determined by the training data. Training data may be obtained from analogous wells.
  • the input of the training process may be a vector of features extracted from the training data, and the output may be a vector of target lithology description. For example, if we have three lithology types (sandstone, siltstone and shale), the target output vector may be represented by [1,0,0], [0,1,0], [0,0, 1] for the sandstone, siltstone and shale, respectively.
  • the prediction results will be the similar 3-element vector. With a hard decision of which value of the element is highest, we can determine the predicted lithology type. For example, if the neural network output is [0.9, 0.04, 0.01], then we may make a conclusion that the predicted lithology is sandstone.
  • online optimization may provide drilling parameters to be applied in drilling in real-time as well as associated analysis for drilling engineers.
  • the parameters and analysis may include: the estimated cost for drilling under the current drilling conditions, and adjusted drilling parameters to maximize drilling efficiency while reducing cost, estimated depth or time to pull out the drill bit.
  • sample data analysis and input data specifications are described according to an embodiment of the present invention.
  • the sample data comprises three sets of well data that were collected for three S-curve wells.
  • the sample data for analysis may be categorized into the following seven categories: Measurement While Drilling (MWD), an oilfield instrumentation and data acquisition systems Electronic Drilling Recorder (EDR), tour sheet, bit report, drilling software data, geological report, sidewall cores report, and LWD (logging while drilling) (neutron density, gamma ray, and resistivity).
  • MWD Measurement While Drilling
  • EDR Electronic Drilling Recorder
  • tour sheet bit report
  • drilling software data drilling software data
  • geological report geological report
  • sidewall cores report sidewall cores report
  • LWD logging while drilling
  • the sample data comprises data from three sources: the oilfield instrumentation and data acquisition systems EDR, geological data and drilling software data.
  • the sensors 101 such as those incorporated in the EDR, supplies well site data comprising real-time MWD data and logging data 204 such as, bit record 2041 , mud log 2042, gamma ray log 2045, neutron density log 2046, resistivity log 2047 and micro log 2048 as shown, according to an embodiment, in Tables 1 , 2, 3, and 5, below.
  • the list of parameters may be revised or expanded to match data capture capabilities.
  • the drilling data may be downloaded as functions of either drill time or well depth.
  • the geological data may comprise well path survey information, geological report data and petrophysical evaluation data.
  • the drilling software data may comprise an integrated well life cycle report comprising geology, construction, drilling, and completion. Cost information may also be provided as set out in the Table 4. Other types and organizations of data may be provided as will be appreciated by those skilled in the art.
  • the MWD data may be acquired at a data acquisition rate of up to 0.5 meters (depth-based) or 20 seconds (time-based) for the above parameters.
  • Other data acquisition rates may be used beneficially in the discretion of the operator, subject to the capabilities of the particular EDR used at the rig site.
  • the geological report may include the well summary, casing summary, daily drilling summary, bit record table, wireline logging summary, deviation/direction survey report, formation top summary, formation evaluations, and core sample descriptions. This report may be based on the information gathered from the geological descriptions and the wireline log analysis. This information may form part of the historical data from analogous wells.
  • the sidewall cores report may comprise a parameter for irreducible water saturation for each formation (average value) (in fractional (%)).
  • the neutron-density log 2046 and gamma ray log 2045 may provide data at a data acquisition rate of 0.1 meters. Other data acquisition rates may be used beneficially in the discretion of the operator, subject to capabilities of the particular logging tools employed at the drilling site.
  • the availability of the data inputs 201 to be used in the prediction engine 20 may be checked and their status marked as (1) available, (2) available but needs to be calculated indirectly, and (3) not available.
  • the units and range of data used by each product and reporting standard may be recorded.
  • the mud density information was not available from the EDR drilling data for the sample data provided herein, but it was available from the Electronic Tour Sheet (although the sampling frequency is much lower, i.e. up to 6 times per day, compared to other parameters, as low as 20 seconds per sample).
  • the different formats of reports and data sheets were imported into a Matlab(tm) workspace.
  • the Excel(tm) data was imported by a Matlab(tm) built-in method.
  • the mud record 2200 was extracted from a ETS 3.0 tour sheet and converted to text format using a custom conversion program which written using the Java(tm) programming language.
  • the .xsl style-sheet file was reprogrammed for Java(tm) and the bit record extracted and converted to .csv format.
  • the drilling software database was also linked to Matlab(tm) and the parameters imported and saved to workspace.
  • a unit conversion table was identified and a unit conversion function established in the Matlab(tm) program.
  • Tour Sheet It was observed from the Tour Sheet that, for example, in the mud record, some data entries are empty and some are inputted incorrectly.
  • Multiple data sources i.e. Tour Sheet and other sources
  • Tour Sheet and other sources have different sample rates, are either uniformly or non-uniformly sampled, are either time-stamped or depth-stamped, and these need to be aligned to form a unified data source.
  • a set of downhole mechanics measurements e.g. ROP 2032, WOB 2033, torque 2036, rotary speed 2035, were collected at a specific depth. If any of these measurements are negative or a zero value, the set of measurements is identified as invalid data by data validation, and will be removed from processing. The remaining validated data may then be fed into a data smoothing function to remove noise.
  • Unit Conversion e.g. ROP 2032, WOB 2033, torque 2036, rotary speed 2035
  • the data inputs 201 may employ metric units by default.
  • the system may also enable unit conversion between Metric and English (i.e. imperial) unit systems. Table 6 for the inputs and outputs of unit conversion between these two unit systems is provided below.
  • the unit conversion may be carried out by establishing a data structure whose sub-fields are category, field, and conversion constant and loading the data structure into the software's graphical user interface.
  • the applied bit wear prediction 2056 may be based on the bit-rock interaction model developed in T. Richard and E. Detournay, "Influence of bit rock interaction on stick slip vibrations of PDC bits" Society of Petroleum and Engineers, Texas, October 2002, which is incorporated by reference herein.
  • the model accounts for the cutting action of a single cutter and allows forces between the WOB 2033 and torque 2036 to be expressed without involving rock properties.
  • the coefficients which characterize these relationships will change.
  • the present invention may monitor these coefficients and detect their variations on a real-time basis.
  • the cutting action of each cutter is represented by two independent processes, namely, the cutting process and the friction process.
  • the forces torque, T , and WOB, w are defined by:
  • ⁇ a bit radius (m) d : depth of cut (m) ⁇ : bit shape factor (dimensionless) ⁇ : cutting force coefficient (dimensionless) ⁇ : intrinsic specific energy (Mpa) ⁇ : friction coefficient (fractional)
  • bit wear prediction 2056 may be based on the detection of a change in ⁇ and ⁇ , which leads to the change in the slope of £ as a linear function of 5 . If we define the product of ⁇ as a bit efficiency factor (BEF), the problem of bit wear prediction 2056 becomes a prediction of BEF using the measurements of the specific energy and the drilling strength.
  • BEF bit efficiency factor
  • a Rao-Blackwellised Particle Filter (RBPF) approach as described, for example, in Arnaud, E. and Memin, E., "An efficient Rao-Blackwellized particle filter for object tracking," in IEEE International Conference on Image Processing, 2005. ICIP 2005, vol. 2, 2005, II, incorporated by reference herein, may be applied to detect the change in BEF.
  • the RBPF approach is based on combining two functional modes, one of them corresponding to the levels of BEF and the other corresponding to the bit- rock interaction with respect to each level of BEF. Each of the defined functional modes is modeled by a state space system.
  • the RBPF comes from the possibility of combining the features of alman filters ( F) and those provided by the particle filters (PF).
  • KF The role of KF is to treat each of the defined state space representations and PF is used to select the one to be treated.
  • PF the underlying objective is achieved by two parallel linear regressions obtained by using RBPF and the prediction is made when the process model moves from the model which corresponds to the process when the BEF changes.
  • the RBPF algorithm may produce an estimate of all 9 models simultaneously and select one winning model and its corresponding mode level. If the winning model at time T comes from the mode level 3, then the prediction at time T for the bit wear is 3. The bit wear may have a variation trend from mode level 0 to BEF level 8, so that any change of the winning model will reflect the worn bit.
  • f x is referred to as the system model and can be non-linear and non-Gaussian. According to an embodiment, for bit wear prediction, it is linear. Then the mode selection of RBPF from S. Tafazoli, et al., "Hybrid system state tracking and fault detection using particle filters", IEEE Transaction on Control Systems Technology, vol. 14, no. 6, pp. 1078-1087, November 2006, incorporated by reference herein, may be applied as follows and its graphic illustration is shown in FIG. 4, which shows a mode selection used in the prediction engine 30:
  • x; for I,— , N .
  • ⁇ y, /1 ' J : for all ⁇ ? ,j ⁇ , ⁇ ⁇ , ⁇ ; Normalize the weights of particles in mode t .
  • Step 4 Resampling: Resample N new particles
  • x ',for 1, - -,N
  • x ',for l,- - -,N
  • Set t t + ⁇ and go to step 2.
  • the level of bit wear (i.e., the mode) may be assumed to vary from 1 to K .
  • the measurement function is not only a function of BEF, but also the mode variable, which defines the measurement model.
  • the bit wear prediction 2056 may be developed according to the block diagram 40 shown in FIG. 5.
  • the inputs are the estimate of BEF at t - 1 401 , the covariance estimate of BEF at / - 1 402, and the estimated probability density function (PDF) of the mode variable at the time instant t - 1 403.
  • PDF estimated probability density function
  • These variables may be input into a RBPF particle filter 404 to compute the estimated PDF of the mode variable at the time t 410 by the mode selection scheme proposed above.
  • computation of the mean of the mode variable at t 405 is performed.
  • This mean value along with the estimate of BEF at t - 1 401 and the covariance estimate of BEF at t - 1 402 is input into a Kalman filter (KF) 406.
  • KF Kalman filter
  • the prediction step of KF 406 yields the BEF prediction 407 and the update step of KF 406 produces the estimate of BEF at t 408 and the covariance estimate of BEF at / 409.
  • the mode variable is the value of interest.
  • the results from the bit record 2041 for bit #7 from Well # 1 were input for prediction and FIG. 6 displays the bit wear prediction results 50 for bit #7 from Well #1.
  • Complementary information may also be displayed to the user or otherwise provided from the historical bit data of the analogous wells showing what bit was used most often in other wells around the current well being drilled, the depth at which the bit was pulled, its dull coding when it was pulled, and the time the bit spent at the bottom of the well.
  • An example of complementary information for analogous Well #2 and Well #3 is shown in Table 7. This information may be part of the bit record 2041 and used in bit wear prediction 2056 to enhance the cost optimization of bit replacement. [000146] Table 7. Complementary information for analogous Well #2 and Well
  • Neural networks may be used to establish complicated non-linear mapping between inputs and outputs.
  • the relationship between formation lithology types and real-time drilling data is usually non-linear.
  • Using a neural network approach avoids an exact mathematical modelling of this relationship, and instead non-linear mapping is established using a neural network.
  • Training data may be applied to fine tune the network parameters in order to obtain one which is closest to the real world.
  • a multilayer neural network may be established to represent the mapping between the real-time drilling data and the lithology class.
  • the features may be rock compressive strength, formation drillability coefficient, dynamic drilling response, and formation shear strength.
  • Feature extraction (the calculation process) may be a factor in the prediction performance of neural networks. Poorly selected features may deteriorate the network performance and thus may result in significant errors during prediction.
  • the features should be selected to best represent the input of the non-linear mapping with less dimension, noise and redundancy.
  • One may choose measurements which may better reflect the rock properties as input features, for instance, the rock strength.
  • Rock properties may be inferred from two different sources: from elastic parameters via the seismic amplitude versus offset (AVO) data, and indirectly from mechanical drilling data by rock-bit interactive formulations.
  • AVO seismic amplitude versus offset
  • the mechanical measurements may include ROP 2032, WOB 2033, torque 2036, and rotary speed 2034 such as shown in realtime drilling data 203 in FIG. 3.
  • Other information may include mud density and bit size.
  • the raw data may be pre-processed as discussed above. After pre-processing, four features may be calculated at each point of drilling depth.
  • CCS is normalized with respect to a mud density of 9.5 ppg.
  • a, , i3 ⁇ 4 , a 3 are dimensionless correction factors determined by the laboratory experiments which may be used to estimate CCS at any mud density value.
  • neural networks for lithology prediction 2057 may be advantageous due to certain properties of neural networks, including learning by experience (neural network training), the ability to generalize (map similar inputs to similar outputs), and their robustness in the presence of noise and multivariable capabilities, which allows them to deal with the highly complex and uncertain problem of lithology determination.
  • the basic processing element of neural networks is called a neuron.
  • FIG. 7 displays a block diagram of a neuron 60.
  • Each neuron has multiple inputs p v p 2 ,' ",P k 600.1, 600.2, ... 600.& where k is limited only by desired number of features and the available computing processing power, and a single output a 620.
  • the net input n 606 is passed through the activation function F(-) 630, which may be either linear or non-linear, producing the neuron output a 620.
  • w l ,w 2 ,- - -,w k 610.1, 610.2, ... 6 ⁇ 0.k are neuron weights
  • b 604 is a neuron bias.
  • the neural network 70 may have multiple layers 640, 641 , 642, and 643. Two or more of the neurons 60 described above may be combined in a layer.
  • the neural network 70 may be supplied with inputs pi 600, p2 601, p3 602, and p4 603. The number of inputs in a layer need not equal the number of neurons in a layer.
  • Each layer has a weight matrix 610, a bias vector b 604, a function ⁇ 605, an activation function F n 630 and an output vector a 620.
  • the output vector a 620 from each intermediate layer is the input to the following layer.
  • the layers in a multilayer network 70 play different roles.
  • a layer that produces the network outputs al 621, a2 622, a3 623 is an output layer 643.
  • the other layers are hidden layers 640, 641 , 642.
  • the neural network 70 shown in FIG. 7, for example, has one output 643 layer and three hidden layers 640, 641, 642.
  • training procedures may be applied once topology and activation functions are defined.
  • the activation function may be selected from the pure linear, log-sigmoid, tangent sigmoid, etc.
  • a set of input data and correct output data i.e. lithology type data
  • the network using the set of training input, produces its own output. This output may be compared with the targeted output (the known lithology type) and the differences are used to modify the weights and biases.
  • Methods of deriving the changes that might be made in a neural network or a procedure for modifying the weights and biases of a network are known as learning rules.
  • test set i.e., a set of inputs and targeted outputs (i.e. correct output data) that were not used in training the network, is used to verify the quality of the obtained neural network.
  • the test set is used to verify how well the neural network can generalize.
  • the inputs to the neural network 600, 601, 602, 603 are the four features from feature extraction as discussed above: rock confined compressive strength (CCS), formation drillability coefficient (FDC), dynamic drilling response (DDR), and formation shear strength (FSS).
  • CCS rock confined compressive strength
  • FDC formation drillability coefficient
  • DDR dynamic drilling response
  • FSS formation shear strength
  • the learning rule used here may be back propagation learning. Backward propagation learning may be composed of the following three steps:
  • the outputs of the neurons in the last layer are considered the network
  • Step 3 weight and bias updates: finally, the weights and biases are updated using the approximate steepest descent rule:
  • the outputs are the weighting values associated with each of possible lithology types.
  • the largest value of the output weight corresponds to the predicted type of lithology.
  • the overall processing system for real-time lithology prediction 2057 may be as shown in FIG. 9.
  • the inputs to the system may include downhole mechanic measurements 801 (e.g., ROP 2032, WOB 2033, torque 2036, rotary speed 2034), mud log 2042 (e.g., mud density), and bit record 2041 (e.g., bit size, PDC cutter size).
  • the downhole mechanics measurements 801 may undergo pre-processing 810 before being fed into the feature extraction unit 820.
  • the pre-processing 810 may include data validation 81 1 and data smoothing 812.
  • the downhole mechanics measurements 801 for example, ROP 2032, WOB 2033, torque 2036, and rotary speed 2034, are collected at a specific depth. If any of these measurements has a negative or zero value, the set of measurements may be designated as invalid data by data validation 81 1 , and may be removed from processing.
  • Data smoothing techniques 812 such as moving average filtering, are then applied to the data to remove noise.
  • the pre-processed downhole mechanics measurements 801 along with mud log 2042 and bit record 2041 are inputs into feature extraction unit 820 to obtain, through further computation, the four features of interest, that is rock confined formation shear strength (FSS) 831 , dynamic drilling response (DDR) 832, formation drillability coefficient (FDC) 833, and compressive strength (CCS) 834.
  • FSS rock confined formation shear strength
  • DDR dynamic drilling response
  • FDC formation drillability coefficient
  • CCS compressive strength
  • the computation of the four features may be according to the equations discussed above.
  • These four features may be used as inputs to a multilayer neural network 70 trained as discussed above.
  • the output of the multilayer neural network is the predicted lithology type 840.
  • Pore Pressure estimation 2053 may be carried out in three stages: before drilling, while drilling and after drilling. Since the drilling fusion engine uses realtime drilling-while-processing for drilling optimization and decision support, while- drilling pore pressure estimation approaches are employed according to embodiments of the invention. According to an embodiment, two quantitative estimation techniques may be used: a D-exponent-based approach or a Sigmalog based approach.
  • d-exponent d-exponent
  • ECD Equivalent circulating density
  • the third approach is to use a linear scale for depth D but a logarithmic scale for d mod .
  • a straight-line normal pressure trend line may be obtained having intercept (d moi ) and exponent m ;
  • the d moi parameter since the d moi parameter considers only the effects of bit weight, bit diameter, rotary speed 2034, and mud log 2042, changes in other drilling variables such as bit type, bit wear, mud type, may still create problems in interpreting the obtained results. In addition, extreme changes in the variables included in the d moA calculation may create big uncertainties. Therefore any new trend may be established for the changed conditions. The utility of the d-exponent is diminished when the mud density is several pounds per gallon greater than the formation pore pressure gradient. Because of the overbalance, the penetration rates may no longer respond significantly to changes in formation pressure. Under these conditions, increases in drilling fluid density may cause an erroneous shift in the modified d-exponent plot, which may yield higher pore pressure readings.
  • Sigmalog was developed in the Po Valley in the mid-seventies as a joint venture between AGIP and Geoservices. The aim was to solve the shortcomings of the d-exponent while drilling overpressured sequences of carbonates, marls and silty shales in deep wells. Sigmalog model calculate the rock strength which is also derive from drilling parameters weight on bit 2033, rotary speed 2034, rate of penetration 2032, and bit diameter. The procedures are:
  • Step 1 Calculate total rock strength (i.e., Sigma raw)
  • AP Differential pressure of mud to formation fluid corresponding to the normal hydrostatic gradient (kg/cm 2 )
  • n Factor expressing the time required for the internal pressure of cuttings not yet cleared from the bit face to reach mud pressure
  • Step 3 Establish a normal trend line Ja r ⁇ on plotted ⁇ / ⁇ values. At any point, AP is calculated as:
  • Step 4 The formation pore pressure gradient is calculated using the following equation:
  • Sigmalog may be used in place of d-exponent, but since the method is not easy to use it may be ill-suited to unexplored basins.
  • the limitations of Sigmalog are the same as for d-exponent. According to an embodiment, its use might be entirely restricted to clays and shales.
  • Table 8 shows a list of data inputs for pore pressure estimation.
  • FIGs. 1 1a, l i b and 1 1 c illustrate the formation pore estimation results 2703 using d-exponent 2702 and Sigmalog 2701 methods. Intermediate d-exponent and Sigmalog calculations are also shown.
  • the normal pore pressure gradient used in the calculation is a constant 0.4650 psi/feet. Since the D-exponent based PP estimation equation is a function of normal pore pressure, any variation of the normal pore pressure gradient may change the real-time estimation results.
  • the software has the capability of programming the mud weight to watch if it gets to a predefined value, say, 1 or 1.25 pp, over the programmed/estimated PP. In such cases, a flag may be raised thus averting a potential safety concern.
  • the inputs to the trend fusion engine 308 are grouped into two categories: historical data and realtime data.
  • the historical data inputs as shown in Table 10 may be estimates based on mining nearby historical drilling well sites (i.e., analogous well data). This is based on the assumption that nearby analogous wells 302.1 to 302. n may have similar geological properties. Therefore, this data may provide useful information to new wells and drilling operations thereof for potential safety indications.
  • trend fusion may be performed on historical data from the historical database 309 from multiple analogous wells 302.1 to 302.3n to minimize the uncertainties and noise from individual wells in historical 303.1 to 303.3n drilling data.
  • FIGs. 12a, 12b, 12c, 13a, 13b, 13c, 14a, 14b, 14c show trend fusion performed on historical data for Well #1, Well #2, and Well #3.
  • FIGs. 12a, 13a, and 14a show D-exponent pore pressure estimation results 2801, 2901, and 3001 for Well #1, Well #2, and Well #3.
  • FIGs. 12b, 13b, and 14b show pore pressure results 2802, 2902, and 3002 for Well #1, Well #2, and Well #3.
  • FIGs. 12c, 13c, and 14c show a pore pressure (PP) estimation by linear regression 2803, 2903, and 3003 for Well #1, Well #2, and Well #3.
  • FIGs. 12c, 13c, and 14c show the D-exponent based pore estimation results (dashed curve in the Pore Pressure (PP) estimation plot) and their associated trend analysis results (overlayed by a solid line in the PP estimation plot) by linear regression (with the associated statistical parameters tabulated in Table 1 1) for Well #1, Well #2, and Well #3, respectively. Since the D-exponent based PP estimation method may use a process of line fitting to get the intercept and slope of the trend line of D-exponent, the estimated intercept and slope feeds into the PP estimation equation to obtain the PP.
  • FIG. 15a, FIG. 15b, FIG. 15c, and FIG. 15d illustrate the trend fusion engine 308 according to an embodiment in fusing pore pressure estimation results 3101, 3102 and 3103 from three different wells into optimized pore pressure estimate 3104. After performing the pore pressure estimation 316 for each well site, a centralized fusion algorithm may be carried out to integrate the results into an optimized pore pressure estimation with reduced uncertainties.
  • pore pressure estimated in Well # 1 shows a spike around 2000 meters which is caused by bit transition, where the measured mechanical sensors provide certain erroneous readings and should not affect the pore pressure estimation. After trend fusion, this effect may be effectively smoothed out. [000196] Table 11. Trend analysis parameters for three wells
  • R-square is defined as the ratio of the sum of squares of the regression (SSR) and the total sum of squares (SST)
  • real-time data inputs may be provided from ongoing drilling profiles to assist in identifying and addressing safety concerns while drilling. They may be conveyed to the trend fusion engine 308 as inputs to find any deviations from the static data inputs. As set out in Table 12, these inputs may include:
  • an expert decision engine 328 may be employed in the system to address the safety concerns, which may be based on feeding it both historical data and real-time data and determining potential safety issues according to the embedded expert decision rules.
  • drilling bit hydraulics may be another component in the drilling optimization software.
  • the well fluid column contains several sections of different fluid densities.
  • the variation of pressure with depth in this type of complex fluid column must be determined by separating the effect of each fluid segment.
  • the hydrostatic pressure 3241 at any vertical distance depth may be calculated by:
  • TVD j true vertical depth at the i-th segment, ft
  • FIG. 15e shows a schematic representation of pressure losses in the circulating system 490 having a downhole 499, a drill string 496, a bottom hole assembly 497 and a drill bit 498.
  • the factors affecting circulating pressure may include:
  • Stand pipe pressure 495 Pressure loss in drill string 492 + Pressure loss in BHA 493 + Pressure loss across the bit 494 + Pressure loss in the annulus 491
  • BHP Hydrostatic pressure + Pressure loss in the annulus [000215] Under the dynamic condition, the only effect on the bottom hole pressure may be pressure loss in the annulus.
  • Downhole hydraulics parameters including equivalent circulating density (ECD) 2074, annular velocity, annular pressure loss 3243, jet nozzle pressure loss, hydraulics horsepower, jet velocity, pore pressure gradient, and jet impact force may be calculated using the equations below:
  • API annular pressure loss
  • psi MW mud weight
  • ppg TVD rue vertical depth
  • the pressure reduction through swabbing may be relevant when tripping pipe because: the balancing pressure is the static mud hydrostatic rather than the higher equivalent circulating density, there may be repeated swabbing as each stand is pulled, and the "piston" effect may affect every permeable formation in open hole.
  • the pressure reduction may be minimized by pulling the drill string at a slower speed and keeping mud viscosity as low as possible, bearing in mind that hole cleaning and cutting lift properties may have to be maintained while drilling.
  • a safety or trip margin may be calculated to ensure that the pressure reduction does not create an underbalance. In the static condition, the trip margin may be calculated by:
  • TripMargin (psi) Hydrostatic pressure (psi) - Formation pressure (psi/ft)xTVD(ft)
  • the trip margin may be calculated by:
  • TripMargin (psi) APL(psi) + TripMargin(static)(psi)
  • FIGs. 16a, 17a, and 18a show calculated annular pressure loss (APL) 3201, 3301 and 3401 for Well # 1, Well #2, and Well #3.
  • FIGs. 16b, 17b, and 18b show calculated equivalent circulatory density (ECD) 3202, 3302 and 3402 for Well # 1, Well #2, and Well #3.
  • FIGs. 16c, 17c, and 18c show a bottom hole pressure (BHP) 3203, 3303 and 3403 for Well # 1 , Well #2, and Well #3.
  • FIGs. 16d, 17d, and 18d show calculated trip margin 3204, 3304, 3404 for Well # 1 , Well #2, and Well #3.
  • friction may also be considered in drilling operations as the drill string is tripped in or out, or rotated on or off bottom. Friction may affect the solid mechanics calculations, such as torque 2073 and drag 2052, as well as the hydraulics calculations, including surge, swab, and hook load estimation during cementing. Beyond the conventional torque and drag calculation model such as described in C.A. Johancsik et al., "Torque and drag in directional wells-prediction and measurement", June 1984, SPE 1 1380 and M.
  • the friction coefficient is defined as:
  • the drill string may be simultaneously rotated and tripped in or out and it may also be simultaneously rotated and reciprocated.
  • the drag force can thus be given as:
  • the rotary friction coefficients may be estimated from the above equation, which will lead to the calculation of F n .
  • the entire drill string may be segmented into small elements (e.g., 0.2 m per segment), and the calculation starts at the bottom of the drill string and proceeds stepwise upward.
  • Each short element of the drill string contributes small increments of axial and torsional load.
  • FIG. 19 shows the forces acting on a short, slightly-curved element.
  • the net normal force F n is the negative vector sum of normal components from the weight, w , and from the two tension forces, F t and F t +iF t .
  • the magnitude of the normal force is
  • the outside diameter of the drill collar or coiled tubing may be used to calculate the values in the above equation.
  • the outside diameter of the tool joint may be used.
  • the calculation of normal force may be done by estimation of Bouyed weight of each element together with the real-time inclination and azimuth information.
  • FIG. 19 shows a force balance 150 on a drill string element 1501.
  • Normal force 1502 results from the combination of the two tension forces, F t 1504 and F t +AF t 1503.
  • F t +AF t 1503 is comprised of component forces ⁇ + ⁇ 1505 and a + Aa 1506
  • F t 1504 is comprised of component forces ⁇ 1507 and a 1508.
  • FIGs. 20, 21 , and 22 display the measured real-time drilling data, torque 1601 (FIG. 20) and WOB 1602 (FIG. 21).
  • FIG. 22 displays the resulting estimated rotary friction factor 1603 for the sample data.
  • the porosity 2075 may be taken from the compensated neutron porosity value.
  • the permeability 2076 may be calculated based on empirical porosity-permeability correlation models as described in A. Aziz et al., "Permeability prediction: core vs log derived values", 1995 International Conference on Geology, Geotechnology and Mineral Resources of Indochina, B. Balan et al., "State-of-the-art in the permeability determination from well log data part 1 : a comparative study, model development", 1995, SPE 30978, S.
  • the model may include an irreducible water saturation ( S wj ) value, as detailed below.
  • FIG. 24 shows the calculated permeability 1702 based on the neutron porosity 1701 shown on FIG. 23 using an average S WI value of 21% for Well #1. Further, real-time S wj may be estimated if a neutron magnetic resonance log is available. That is, the water saturation may be either obtained from direct measurement data, or from the estimation using the neutron magnetic resonance log data.
  • the volumes 3245, capacities and displacements calculations may provide fundamental indications leading to potential safety issues. These calculations, together with the permeability 2076 & porosity 2075, lithology prediction 2057, pore pressure estimation 2053, may serve as the combined inputs to the expert knowledge system to address a variety of safety issues during drilling. [000237] According to an embodiment, the system may include pressure control calculations, in which the kill sheets and its related calculations may be carried out. These may involve the following calculations:
  • Table 13 lists the calculations in the kill sheet for the pressure control according to an embodiment of the present invention.
  • the calculation was performed with the three sets of sample well data and the results from Well # 1 are provided as an example.
  • the pre-recorded data for the calculation may be read from the oilfield instrumentation and data acquisition systems tour sheet and drilling software database software.
  • the snapshots of the relevant data are shown in FIGs. 24b, 24c, 24d, and 24e.
  • FIG. 24b shows the front page summary of a tour sheet for Well #1 3501
  • FIG. 24c shows the drilling assembly portion of a tour sheet for Well #1 3601
  • FIG. 24d shows the mud record and circulation portions of a tour sheet for Well #1 3701
  • FIG. 24e shows the reduced pump speed portion of a tour sheet for Well #1 3801.
  • Table 14 shows the casing summaries read from the geology report.
  • the calculation results are shown in FIGs. 24f, 24g, 24h, and 24i.
  • FIG. 24f shows a user interface display of calculations for kill sheet #1 3901
  • FIG. 24g shows a user interface display of calculations for kill sheet #2 4001
  • FIG. 24h shows a user interface display of calculations for kill sheet #3 4101
  • FIG. 24i shows a user interface display of calculations for kill sheet #4 4201.
  • Kill weight mud (KWM) [ppg] SIDPP [psi]/0.052/TVD [ft] + OMW [ppg]
  • Initial circulating pressure (ICP) [psi] SIDPP [psi] + KRP [psi]
  • one aspect of the software is to recommend drilling procedures which may result in the successful completion of the well as safely and inexpensively as possible.
  • the drilling engineer may make recommendations concerning routine rig operations such as drilling fluid treatment, pump operation, bit selection, and any problems encountered in the drilling operation.
  • routine rig operations such as drilling fluid treatment, pump operation, bit selection, and any problems encountered in the drilling operation.
  • the use of a drilling cost equation may be useful in making these recommendations.
  • the decision concerning which bit to use often based on some performance criteria, such as total rotating hours, total footage, or maximum penetration rate. Other times, the least expensive bit is chosen. This approach may be satisfactory in areas where practices and costs are constants but it may not be satisfactory where drilling costs are changing, and drilling practices and bit selection vary.
  • the usual procedure may be to break the drilling costs into (1) variable drilling costs and (2) fixed operating expenses that are independent of alternatives being evaluated. Selection of the bit best suited for a specific use may be further complicated by the variable performances and prices of the many types of bits.
  • bit selection may be based on achieving the minimum cost-per-foot. In this way, it may be possible to achieve an optimum relationship between penetration rate, bit footage, rig cost, trip time and bit cost.
  • a common application of a drilling cost formula may be to evaluate the efficiency of a bit run. A large fraction of the time required to complete a well may be spent either drilling or making a trip to replace the bit.
  • the total cost required to drill a given depth AD may be expressed as the sum of the total rotating time during the bit run, t h , the non-rotating time during the bit run, t c , and trip time, t, .
  • the cost-per-foot as related to these variables may be determined by the equation: ' AD
  • C j is drilled cost per unit depth
  • Q is the cost of bit
  • C r is the fixed operating cost of the rig per unit time independent of the alternatives being evaluated.
  • FIGs. 24j, 24k, and 241 display copies of the bit record of Well #1 4202, the bit record of Well #2 4203 and the bit record of Well #3 4204, respectively, which has the information for each bit run including the bit type, hour of running, average ROP, and depth.
  • Tables 15, 17, and 19 list the cost information read from the drilling software. The resulting cost-per- foot analysis results are shown in Tables 16, 18, and 20.
  • the above cost-per-foot calculations may also be used in generating performance proposals.
  • a target cost-per-foot may be calculated based on information from offset well performance. If the target is met or exceeded, the motor may be charged at normal list price; if it is not, then the hourly rate may be adjusted downwards (to an agreed minimum) until the target figure is reached.
  • the user should be reasonably sure that the operator will be able to achieve a minimum ROP. This ROP may be the drill rate needed to meet both the proposed cost-per-foot figure and the pricing structure. After calculating this target ROP, it may be determined whether the target cost-per-foot is attainable.
  • C f is the target cost-per-foot agreed upon with the user.
  • the decision making for expert system recommendations for safety concerns while drilling may be based on inputting both historical data from the historical database 309 and adjusted real-time drilling parameters 329 to the expert decision engine 328 which applies a set of expert decision rules in accordance with the decision trees set out below.
  • the expert decision engine 328 may include six decision trees: mud volume, porosity, permeability, stuck pipe, bit wear and environmental. Further expert decision trees may be added according to the format described below in the discretion of the operator.
  • the expert decision recommendations may appear as prompts to the user during operation of the client or standalone software through the graphical user interface.
  • the user may be warned through the graphical user interface of a danger level: green indicates normal, yellow indicates warning, and red indicates dangerous.
  • the graphical user interface display may also flash to alert the user of the issue.
  • the mud volume decision tree 4300 monitors for a change of the mud volume parameter in real time 4302 and if the mud volume parameter 4301 is greater or equal to than a defined threshold 4303, the pump is stopped 4304, a flow check is performed 4305, a visual check of the tank for input fluid, whole mud, water, oil, and for flow in the drill pipe and annulus 4306 may occur, volume gain may be established 4307 and the reason therefor such as potential blowout or need to release pressure 4308 or potential loss in circulation 4310, and a procedure in regards of potential blowout or need to release pressure 4309 or a procedure developed by the operator or contractor in regards of potential loss in circulation 431 1 may be presented to the user.
  • the defined threshold may be set according to the conditions as follows: 1 ) if losses are less than 20 barrels per hour, or if company guidelines permit losses of 20 barrels per hour while drilling with a water base mud continue drilling; if losses are greater than 20 barrels per hour, try to heal or fix the loss zone while drilling using 2 to 6 pounds per barrel loss circulation material if there are partial losses. For more extreme losses, but not full losses, stop, mix a tank of whole mud with lost circulation material of 6 plus pounds per barrel; if full losses arise, stop, and evaluate options to fix or stop loss circulation: a. set a plug. There are various types (Di-Seal, GUNK etc.), b. set cement plug across the loss zone, c. drill ahead without circulating (drill ahead blind).
  • the porosity decision tree 4400 monitors the porosity 2075 parameter for significant deviation from the historical reference value. If the deviation from the historical value is greater or equal to 60% 4402, expert decision engine 328 determines whether there has been a change in lithology 4405 by comparing the predicted lithology type with a lithology sample collected at a particular depth. Once a sample is collected, the resulting variable may be imported into the program via a user interface for comparison.
  • the expert decision engine 328 panel may flash litho-red 4412, which reflects both porosity change and lithology change (whereas red reflects porosity variations against the historical porosity records only) and prompt the user to contact the drilling engineer and geologist to re-evaluate the drilling parameters for potential hole problems 441 1. If the deviation from the historical value is not greater or equal to 60%, but greater or equal to 40% 4403, it is determined whether there has been a change in the lithology 4405. If there has been a change in the lithology, the expert decision engine 328 panel may flash litho-red and prompt the user to contact the drilling engineer and geologist 4412 to re-evaluate the drilling parameters for potential hole problems 441 1.
  • the deviation from the historical value is not greater or equal to 40%, it is determined whether it is greater or equal to 20% 4404, and if it is, then it is determined whether there has been a change in the lithology 4405. If there has been a change in lithology, the expert decision engine 328 panel may flash litho- red and yellow alternately and prompt the user to contact the drilling engineer and geologist 4414 to re-evaluate the drilling parameters for potential hole problems 441 1.
  • different thresholds than 60%, 40% and 20% may be used according to embodiments of the invention in the discretion of the operator or user.
  • the permeability decision tree 4500 monitors the permeability 2076 parameter for significant deviation from the historical reference value 312. If the deviation from the historical value is greater or equal to 60% 4502, it is determined whether there has been a change in lithology 4506 as described above, and if there has been a change in lithology then the expert decision engine 328 panel may flash litho-red and prompt the user to contact the drilling engineer and geologist 4512 to re-evaluate the drilling parameters for potential hole problems 451 1. If the deviation from the historical value is not greater or equal to 60%, but greater or equal to 40% 4503, it is determined whether there has been a change in the lithology 4506 in the same manner as described for porosity above.
  • the expert decision engine 328 panel may flash litho-red and prompt the user to contact the drilling engineer and geologist 4512 to re-evaluate the drilling parameters for potential hole problems 451 1. If the deviation from the historical value is not greater or equal to 40%, it is determined whether it is greater or equal to 20% 4504 and if yes, it is determined whether there has been a change in the lithology 4506. If there has been a change in lithology, the expert decision engine 328 panel may flash litho-red and yellow alternately and prompt the user to contact the drilling engineer and geologist 4513 to re-evaluate the drilling parameters for potential hole problems 451 1.
  • different thresholds than 60%, 40% and 20% may be used according to embodiments of the invention in the discretion of the operator or user.
  • the stuck pipe 4600 decision tree monitors the stuck pipe decision parameters 4601 such as hole parameters 4602, directions 4603, volumes 3245, filter cake/mud cake 3254, mud density 3276, and bottom hole assembly 4607 parameters for deviation from the historical reference value 4608. If there is deviation, i.e. if you get stuck on the current well where they did not get stuck on the historical well, then a determination of the differences may be made to aid in decision making; the expert system engine 328 determines whether circulation is possible 4609.
  • the stuck pipe decision parameters 4601 such as hole parameters 4602, directions 4603, volumes 3245, filter cake/mud cake 3254, mud density 3276, and bottom hole assembly 4607 parameters for deviation from the historical reference value 4608. If there is deviation, i.e. if you get stuck on the current well where they did not get stuck on the historical well, then a determination of the differences may be made to aid in decision making; the expert system engine 328 determines whether circulation is possible 4609.
  • the expert decision engine 328 reviews the difference in depth, formation, lithofacies, porosity, permeability, mud, filter cake, density, ROP, RPM, and WOB 4613, and the historical operations in the database 4615 and may analyze the probability of stuck, depth, free point, size and density of spotting pill 4617, and estimate optimized time and cost for the current operation 4619.
  • the expert decision engine 328 notes the difference in depth, formation, lithofacies, porosity, permeability, mud, filter cake, density, ROP, RPM and WOB 4613, and the historical operations 4615 in the database and may analyze the probability of stuck pipe, depth, free point, size and density of spotting pill 4617, and estimate optimized time and cost for the current operation 4619.
  • the analysis of the probability of getting stuck is made with regard to whether there is circulation or no circulation and the directionality during getting stuck: drilling ahead, tripping in the hole, tripping out of the hole, with drill string or electric logging tools.
  • drilling the free point is analyzed via an equation (provided a few paragraphs below), which conveys the depth of stuck pipe. While tripping, it is evident (number of pipe, etc.) or may be determined from sensors 101.
  • the size of the spotting pill may be done by rule of thumb. According to an embodiment, the spotting pill may be approximately 50 barrels. According to an embodiment, the density of spotting pill may be 1 to 1.5 pounds per gallon more than the drilling mud.
  • the expert decision engine 328 may therefore determine stuck pipe situations while drilling.
  • the free point constant for the stuck pipe may be determined from the depth at which the pipe is stuck and the number of feet of free pipe may be estimated by the drill pipe stretch table (see Table 21).
  • the length of free pipe may be calculated as follows: stretch (in.) x free point constant
  • FCF free point constant
  • the bit wear prediction decision tree 4700 monitors the current bit wear prediction 2056 and historical data for the same bit 314.
  • the bit manufacturer 4702, IADC (International Association of Drilling Contractors) code for bit type 4703, bit size 4704, recommended bit hours 4705, depth in 4706, time in 4707, are checked against historical data for same bit run 4709, checked against historical data for run time 4710, and checked against historical data for bit total depth 471 1 for deviation from the historical data 4713. If there is a deviation from historical data, i.e.
  • the bit is checked as to whether it is a hard formation bit 4721, if it is then is it checked to determine whether there has been a change in drilling parameters 4722 and if not, the expert system engine prompts to trip for bit 4723 to change the bit. If yes, the parameters are reviewed and changed as required to maintain the rate of penetration 4720. If there has been a rate of penetration decrease, the expert decision engine 328 also determines if RPM (revolutions per minute, which is a measure of bit rotary speed) has changed 4715. If RPM has not changed, then the expert decision engine 328 prompts to trip the bit if it is close to the bit manufacturer's recommended hours 4728.
  • RPM repetitions per minute
  • the reason is reviewed 4730 and the expert decision engine 328 determines if any changes will make a difference 4727. If changes to the drilling parameters do not make a difference then the expert decision engine 328 prompts to trip the bit 4723. If the changes make a difference, then the expert decision engine 328 continues to monitor that bit hours are under manufacturer's recommended hours 4725. If there has been a rate of penetration decrease, the expert decision engine 328 also determines if WOB has changed 4716. If WOB has not changed, then the expert decision engine 328 trips the bit if close to the bit manufacturer's recommended hours 4728. If the WOB has changed, then the reason is reviewed 4730 and the parameters may be optimized to increase ROP.
  • the system trips for bit 4729. If there has been a rate of penetration decrease, the expert decision engine 328 also determines if the drill string is vibrating 4717. If the drill string is vibrating, the geology is reviewed 4734 and the system optimizes the parameters to minimize vibration and increase the ROP, and if an increase is not possible, the system trips for bit 4731. If the expert decision engine 328 determines that the drill string is not vibrating, the system continues drilling 4733.
  • the geology review 4734 may be carried out by checking the current lithology sample to see if there is any evidence of a different geology than that indicated by the predicted lithology. For example, is there any chert or other hard or soft material at this depth that would indicated this formation is such that vibration would occur.
  • Drilling parameters may be optimized by changing pump rates (stokes per minute, SPM), RPM, and WOB to see if any minor variation causes an improvement.
  • the environmental factors decision tree 4800 may review the environmental factors 4801 to determine if there is deviation from historical values 4802.
  • the environmental factors 4801 reviewed may include what type of fluid that was spilled (e.g. water based mud, oil based mud, formation fluid, engine fluid, gas or diesel?), what was the volume of spilled fluid, how sensitive are the surrounding areas to the fluid (e.g. dry land or water, lake, pond, swamp, river, or creek, or underground water table?), is the well in an environmentally sensitive area, such as a 'wet lands', a 'national park', etc., and other environmental factors conceived of by the drilling contractor and operator.
  • the expert decision engine 328 determines how these differences differ from the historical values 4803, and then the expert decision engine 328 may compare the operations at the time of the hazard with historical wells in the area or other wells in the drilling area that have had similar incidents and the potential outcomes are analyzed 4804. The potential outcomes are analyzed by comparing them to the historical well problems in the same area, if any, or to outcomes given the contingency plans and other wells in the drilling area. The expert decision engine 328 then provides probabilities of the potential outcomes in different operations for planning 4805 and the drilling operations are then re-evaluated to make changes 4806. The probabilities of outcomes are generated by having more than one of these incidents, problems, concerns or issues occurring in the area on or by other rigs or in the historic well data. When the contingency plan is developed by the operator or contractor, they may research this data.
  • a graphical user interface (GUI) 5001 may be provided for the software.
  • the functions in the prediction engine 20, the trend fusion engine 308 and the expert decision engine 328 may be integrated into the software so that the full set of relevant drilling parameters, prediction tools, optimization tools, analytic tools and the expert decision system are available to the user through a common interface.
  • the GUI 5001 may be made up of three subpanels which govern different aspects of the software's functionality. According to the embodiment as shown in FIG. 31 , those subpanels are a Historical Reference Data Panel 5002, a While-Drilling Data Panel 5003 and Expert Decision Panel 5004.
  • the Historical Reference Data panel 5002 may display the historical drilling data obtained from the nearby analogous wells as reference information.
  • the Expert Decision panel 5004 may allow comparison between the real-time drilling parameters and operations with the historical drilling data in the Historical Reference Data panel 5002. Since the reference information in the Historical Reference Panel 5002 is time/depth dependent, there may be a slider and text box for users to provide the depth and the information may be retrieved accordingly.
  • the reference data may comprise information relating to the drill bits including bit type, depth in/out, diameter, manufacturer, manufacturer recommended hours, historical hours of running and bit wear status when pulled out; information relating to lithology including lithology type, estimated port pressure gradients, normal pore pressure for the region, porosity 2075, permeability 2076 and water saturation 3246; information relating to mud records including mud weight, mud type, viscosity, and PH; and information relating to the bottom hole assembly (BHA) including assembly components, outer diameter, inner diameter, pipe length, and weight.
  • BHA bottom hole assembly
  • the While-Drilling Data panel 5003 may show the real-time drilling data either from the real-time data acquisition system or from the data simulator.
  • the while-drilling data may be loaded from the drop-down list on the top of the While- Drilling Data panel 5003 if it is from the real-time data acquisition system.
  • the start button may initiate the data acquisition result in the plotting of the real-time drilling data as well as generation of the associated engine outputs.
  • the simulator may be enabled when no real-time data is loaded. In this case, the start button may initiate the simulator to plot the simulated data in an animated fashion.
  • the simulator may simulate various cases of potential problems in the drilling process which may not appear in the real-time data.
  • the drilling fusion prediction engine 20 and the expert decision engine 328 may be triggered immediately after the simulation or loading starts.
  • the views of the historical data panel 5002 and the expert decision panel 5004 may be enabled or disabled by clicking the check boxes on top of the while-drilling data panel 5003. [000282] Also on top of the while-drilling data panel 5003, a drop-down list may be provided for the user to switch between five different categories, which may represent five different views of current drilling parameters (Table 22). All plots for different categories may be synchronized. After clicking the start button, plotting may be executed in real-time so that switching between different categories does not interfere with the plotting processing.
  • FIGs. 32, 33, 34, 35, 36, 37, and 38 show screen shots for these five categories: GUI-main display 5100, GUI-MWD (metric units) 5200, GUI-MWD (imperial units) 5300, GUI-Hydrostatic pressure 5400, GUI-Dynamic pressure 5500, GUI-LWD 5600, GUI-Kill sheet 5700.
  • FIGs. 33 and 34 belong to the same category (i.e., MWD) but with a different unit system.
  • the unit conversion may be enabled by clicking the switch buttons on the bottom of the GUI.
  • the default unit system may be the metric unit system.
  • Real-time data may be plotted in red and the referencing historical data, if applicable, may be plotted in blue or green.
  • the category kill sheet may retrieve the current drilling assembly information to calculate the capacities and volumes along the drill string and pipe, which may be displayed in the first GUI window.
  • the GUI may also allow users to continue calculation by inputting more information in the second, third and fourth GUI windows.
  • the Expert Decision panel 5004 (right panel in FIG. 31) may list the status reports from the six expert decision engines. According to an embodiment, different drilling statuses may be indicated by using different colours on the top of the panel, for example, Normal (green), Warning (yellow), and Dangerous (red). The status may be continuously updated to indicate an ongoing problem while drilling as well as the suggested actions to take. Similar to the Historical Reference Data panel 5002, the Expert Decision panel 5004 may also enable users to examine the drilling parameters in detail at a specific depth.
  • the warning and dangerous statuses may trigger the system to be stopped immediately so that the drilling engineers or operator may examine the problem. Normal status would not interrupt the drilling as all parameters would have been determined to be within the safety requirement.
  • the current potential issues may be shown in the status window for review by the user.
  • the associated problematic parameters may also be flashed in their corresponding displays at the current drilling depth. All potential problems as determined by the expert decision engine may also be recorded in a log for subsequent viewing.
  • FIGs. 39, 40, 41, 42, and 43 show the expert decision engine prompts for mud volume issues for Well #1 5800 and Well #2 6000, a porosity issue for Well # 2 6200, a bit wear issue for Weil #2 6100 and a potential stuck pipe 5900 for Well #1 according to an embodiment. As the potential issues are determined, further suggestions towards the underlying drilling parameters may be displayed in a pop-up window.
  • optimized drilling parameters 2072 may be used by the control software of the drill string to adjust drilling operations in view of the real-time predictions or expert system engine recommendations. Drilling operations may be adjusted manually, by a human operator viewing the optimized drilling parameters 2072, according to an embodiment, on a standard computer display, and acting on them to control drilling operations; alternatively, the optimized drilling parameters 2072 may be automatically fed into the control software for operating a drilling rig, subject to a manual override by a human operator. According to an embodiment, the optimized drilling parameters 2072 may be used to update the geological model and trajectory database to ensure their contents reflect the real-time observations and measurements, or to replace planned parameter values with actual measurements in the drilling model. The prediction-optimization control may form a closed feedback loop for drilling operation, which may reduce the overall drilling cost and risk.
  • Embodiments of the invention, or aspects thereof, may be provided in a computer program comprising computer readable instructions for execution on a computer.
  • the computer program is storable on any suitable computer storage medium so as to comprise a computer program product.
  • Such a computer program may provide a graphical user interface which may allow the drilling operator to view and override adjusted drilling parameters as desired.

Landscapes

  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Geology (AREA)
  • Mining & Mineral Resources (AREA)
  • Physics & Mathematics (AREA)
  • Environmental & Geological Engineering (AREA)
  • Fluid Mechanics (AREA)
  • General Life Sciences & Earth Sciences (AREA)
  • Geochemistry & Mineralogy (AREA)
  • Earth Drilling (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

La présente invention concerne des procédés et des systèmes pour des opérations de forage améliorées grâce à l'utilisation de données de forage en temps réel pour prévoir l'usure de trépan, la lithologie, la pression des pores, un coefficient de frottement en rotation, une perméabilité, et un coût en temps réel et pour régler des paramètres de forage en temps réel en fonction des prévisions. La prévision de la lithologie en temps réel est réalisée en traitant les données de forage en temps réel au moyen d'un réseau neuronal à couches multiples. La prévision de l'usure de trépan en temps réel est réalisée en utilisant les données de forage en temps réel pour prévoir un facteur de rendement de trépan et pour détecter des changements du facteur de rendement de trépan au fil du temps. Ces prévisions peuvent être utilisées pour régler des paramètres de forage dans l'opération de forage en temps réel, en fonction d'une annulation par l'opérateur. Les procédés et les systèmes peuvent également comprendre la détermination de divers paramètres hydrauliques en fond de trou et d'un facteur de frottement en rotation. Des données historiques peuvent être utilisées en association avec des données en temps réel pour fournir une assistance de système spécialisée et pour identifier des problèmes de sécurité.
PCT/CA2013/000919 2012-10-31 2013-10-25 Procédés et systèmes pour opérations de forage améliorées utilisant des données de forage en temps réel et historiques WO2014066981A1 (fr)

Priority Applications (1)

Application Number Priority Date Filing Date Title
MX2015005568A MX360775B (es) 2012-10-31 2013-10-25 Métodos y sistemas para operaciones de perforación mejoradas usando datos de perforación históricos y en tiempo real.

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
CA2794094 2012-10-31
CA2794094A CA2794094C (fr) 2012-10-31 2012-10-31 Procedes et systemes pour operations de forage ameliorees utilisant des donnees de forage historiques et en temps reel

Publications (1)

Publication Number Publication Date
WO2014066981A1 true WO2014066981A1 (fr) 2014-05-08

Family

ID=50605935

Family Applications (1)

Application Number Title Priority Date Filing Date
PCT/CA2013/000919 WO2014066981A1 (fr) 2012-10-31 2013-10-25 Procédés et systèmes pour opérations de forage améliorées utilisant des données de forage en temps réel et historiques

Country Status (3)

Country Link
CA (2) CA2794094C (fr)
MX (1) MX360775B (fr)
WO (1) WO2014066981A1 (fr)

Cited By (50)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2598785C2 (ru) * 2014-11-05 2016-09-27 Общество с ограниченной ответственностью "ТатАСУ" Система статистического и нейросетевого анализа данных телеметрии нефтепромысловых объектов
RU2598786C2 (ru) * 2014-08-27 2016-09-27 Общество с ограниченной ответственностью "ТатАСУ" Способ нейросетевого анализа состояния телемеханизированных нефтепромысловых объектов
WO2016209271A1 (fr) * 2015-06-26 2016-12-29 Landmark Graphics Corporation Visualisation de données quantitatives d'opérations de forage concernant un événement de tube coincé
CN106545327A (zh) * 2016-12-09 2017-03-29 北京四利通控制技术股份有限公司 智能司钻钻机控制系统
WO2017070069A1 (fr) * 2015-10-23 2017-04-27 Schlumberger Technology Corporation Utilisation de données de forage pour conditionnements techniques
EP3177806A4 (fr) * 2014-08-04 2018-02-21 Landmark Graphics Corporation Modélisation de l'usure d'un tubage/d'une colonne montante et d'un coefficient de frottement en utilisant des techniques d'inversion discrètes
US9920624B2 (en) 2015-04-09 2018-03-20 Joy Mm Delaware, Inc. System and method of detecting dull and worn cutter bits
WO2018106346A1 (fr) * 2016-12-07 2018-06-14 Safekick Americas Llc Forage basé sur un modèle automatisé
CN108612518A (zh) * 2018-04-20 2018-10-02 重庆地质矿产研究院 一种煤层气井径向微井眼钻孔、水力压裂参数的确定方法
CN109062248A (zh) * 2018-08-02 2018-12-21 西北工业大学 基于自组织网络的空间非合作目标姿轨一体化参数估计方法
CN110261906A (zh) * 2019-06-27 2019-09-20 中国石油化工股份有限公司 多参数岩性信息融合储层预测方法
CN110941866A (zh) * 2019-12-06 2020-03-31 中国石油集团川庆钻探工程有限公司 一种基于固井大数据的环空水泥浆界面设计方法
CN111479982A (zh) * 2017-11-15 2020-07-31 吉奥奎斯特系统公司 具有滤波器的现场操作系统
US10773352B2 (en) 2017-06-05 2020-09-15 Joy Global Underground Mining Llc System and method for determining efficiency of an industrial machine
CN112004989A (zh) * 2018-03-07 2020-11-27 贝克休斯控股有限责任公司 用于显示钻地工具的可靠性的钻地工具监测系统及相关方法
US10883346B2 (en) 2015-12-18 2021-01-05 Schlumberger Technology Corporation Method of performing a perforation using selective stress logging
CN112412390A (zh) * 2019-08-22 2021-02-26 中国石油化工股份有限公司 基于深度学习模型评价固井第二界面的方法及装置
CN112418489A (zh) * 2020-11-05 2021-02-26 中国石油天然气股份有限公司 一种天然气井修井暂堵剂用量神经网络预测方法
NO345487B1 (en) * 2015-05-12 2021-03-01 Halliburton Energy Services Inc Enhancing oilfield operations with cognitive computing
CN112714820A (zh) * 2018-09-14 2021-04-27 埃尼股份公司 通过钻井设备估算待钻的地质构造中孔隙压力值的方法
CN112836406A (zh) * 2021-01-13 2021-05-25 中国地质大学（武汉） 一种考虑不确定阻尼系数的钻进过程钻压建模方法
CN113065664A (zh) * 2021-02-20 2021-07-02 北京三一智造科技有限公司 一种地质层识别模型构建方法和地质层识别方法
CN113073959A (zh) * 2021-02-24 2021-07-06 西南石油大学 一种基于模糊决策树的钻井漏失机理预测及辅助决策方法
US11066917B2 (en) * 2018-05-10 2021-07-20 Baker Hughes Holdings Llc Earth-boring tool rate of penetration and wear prediction system and related methods
US11078785B1 (en) 2020-06-17 2021-08-03 Saudi Arabian Oil Company Real-time well drilling evaluation systems and methods
CN113310829A (zh) * 2021-04-14 2021-08-27 西南石油大学 一种冲击方式下岩石的可钻性测试装置及实验方法
CN113378344A (zh) * 2020-02-25 2021-09-10 中国石油化工股份有限公司 一种用于优化钻井起下钻速度阈值的方法及系统
US11131184B1 (en) 2020-04-29 2021-09-28 Saudi Arabian Oil Company Method and system for determining a drilling hazard condition using well logs
CN113586028A (zh) * 2021-07-21 2021-11-02 宁夏天地奔牛实业集团有限公司 基于数字孪生的反井钻机扩孔刀盘的智能监控系统
CN113792936A (zh) * 2021-09-28 2021-12-14 中海石油(中国)有限公司 一种随钻岩性智能识别方法、系统、设备和存储介质
US11255180B2 (en) * 2017-12-22 2022-02-22 Halliburton Energy Services, Inc. Robust early kick detection using real time drilling
CN114116213A (zh) * 2021-11-24 2022-03-01 上海派拉软件股份有限公司 分布式系统中超额负载估计方法、装置、设备及存储介质
CN114169656A (zh) * 2020-09-11 2022-03-11 中国石油化工股份有限公司 一种基于邻井历史数据的钻井卡钻风险预警方法及系统
CN114357887A (zh) * 2022-01-06 2022-04-15 中国地质大学（北京） 一种基于bp神经网络的复杂井况钻前泥浆漏失预测方法
CN114427430A (zh) * 2020-09-21 2022-05-03 中国石油化工股份有限公司 一种多井实时协同的钻井参数优化方法与系统
CN114526052A (zh) * 2021-12-31 2022-05-24 中国石油天然气集团有限公司 一种钻完井工程风险预测方法及装置
CN114607354A (zh) * 2021-12-14 2022-06-10 中国石油天然气集团有限公司 一种预报钻井复杂情况的方法
US11384612B2 (en) 2017-07-11 2022-07-12 Equinor Energy As Method and system for monitoring influx and loss events in a wellbore
CN114912301A (zh) * 2022-07-18 2022-08-16 中国空气动力研究与发展中心低速空气动力研究所 一种低速风洞全机模型测力试验数据处理与修正系统
CN114997485A (zh) * 2022-05-26 2022-09-02 中国石油天然气集团有限公司 溢流工况预测模型训练方法、装置及溢流工况预测方法
WO2022186893A1 (fr) * 2021-03-03 2022-09-09 Landmark Graphics Corporation Analyse et corrélation d'usure et de comportement de trépan de forage
CN115059448A (zh) * 2022-06-01 2022-09-16 中国石油大学(华东) 一种基于深度学习算法的地层压力监测方法
US11459882B2 (en) 2020-03-20 2022-10-04 Saudi Arabian Oil Company Systems and methods for the determination of lithology porosity from surface drilling parameters
WO2022225537A1 (fr) * 2021-04-23 2022-10-27 Landmark Graphics Corporation Logiciel d'exploitation minière permettant de générer un flux de traitement pour former un puits de forage
US11542788B2 (en) * 2019-03-01 2023-01-03 Brigade Energy Services Drilling out frac plugs
CN116181261A (zh) * 2023-03-29 2023-05-30 安徽理工大学 一种煤矿立井钻井法凿井钻头吸渣口布置的确定方法
CN117034665A (zh) * 2023-10-10 2023-11-10 中国石油大学(华东) 一种基于钻井施工预演的井身结构设计方法
CN118309368A (zh) * 2024-06-07 2024-07-09 陕西延长石油矿业有限责任公司 一种煤矿厚基岩地层一钻完井智能建井方法
WO2024155870A1 (fr) * 2023-01-19 2024-07-25 Baker Hughes Oilfield Operations Llc Détection probabiliste de scénarios de forage
US12291955B2 (en) 2022-03-18 2025-05-06 Saudi Arabian Oil Company Real-time model of rig and bit hydraulics efficiency

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN109519116B (zh) * 2017-09-20 2020-12-25 贵州远东兄弟钻探有限公司 矿产资源勘查中硬地层快速钻进工艺方法
CN110795888A (zh) * 2019-09-12 2020-02-14 四川宏华电气有限责任公司 一种石油钻井风险预测方法
CN111537663B (zh) * 2020-04-20 2022-10-04 中国石油天然气集团有限公司 一种岩性识别剂携带装置及基于其的岩性识别系统和方法
US11946366B2 (en) 2021-02-10 2024-04-02 Saudi Arabian Oil Company System and method for formation properties prediction in near-real time
CN113236223A (zh) * 2021-06-10 2021-08-10 安徽理工大学 一种用于煤矿井下瓦斯防治钻孔的智能设计系统及方法
CN113810025B (zh) * 2021-07-23 2024-05-10 西安交通大学 一种摩擦系数的自适应趋势滤波处理方法
CN115726759B (zh) * 2021-08-26 2024-08-23 中国石油化工股份有限公司 一种用于控制钻井的智能分析决策系统
WO2023028645A1 (fr) * 2021-09-06 2023-03-09 V2 Innovations Pty Ltd Fonctionnements de machines et systèmes de gestion
US12264572B2 (en) * 2022-03-29 2025-04-01 Saudi Arabian Oil Company System and method to predict and optimize drilling activities
CN116446842A (zh) * 2023-03-06 2023-07-18 中建八局第二建设有限公司 一个基于物联网、云架构技术的智能钻机平台检测系统
CN119475817B (zh) * 2025-01-13 2025-03-25 湖南科技大学 一种油区储层参数定量表征方法及系统

Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6424919B1 (en) * 2000-06-26 2002-07-23 Smith International, Inc. Method for determining preferred drill bit design parameters and drilling parameters using a trained artificial neural network, and methods for training the artificial neural network
US20020120401A1 (en) * 2000-09-29 2002-08-29 Macdonald Robert P. Method and apparatus for prediction control in drilling dynamics using neural networks
US20070284147A1 (en) * 2005-02-01 2007-12-13 Smith International, Inc. System for optimizing drilling in real time
US20110174541A1 (en) * 2008-10-03 2011-07-21 Halliburton Energy Services, Inc. Method and System for Predicting Performance of a Drilling System

Patent Citations (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6424919B1 (en) * 2000-06-26 2002-07-23 Smith International, Inc. Method for determining preferred drill bit design parameters and drilling parameters using a trained artificial neural network, and methods for training the artificial neural network
US20020120401A1 (en) * 2000-09-29 2002-08-29 Macdonald Robert P. Method and apparatus for prediction control in drilling dynamics using neural networks
US20070284147A1 (en) * 2005-02-01 2007-12-13 Smith International, Inc. System for optimizing drilling in real time
US20110174541A1 (en) * 2008-10-03 2011-07-21 Halliburton Energy Services, Inc. Method and System for Predicting Performance of a Drilling System

Cited By (84)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP3177806A4 (fr) * 2014-08-04 2018-02-21 Landmark Graphics Corporation Modélisation de l'usure d'un tubage/d'une colonne montante et d'un coefficient de frottement en utilisant des techniques d'inversion discrètes
US10678966B2 (en) 2014-08-04 2020-06-09 Landmark Graphics Corporation Modeling casing/riser wear and friction factor using discrete inversion techniques
RU2598786C2 (ru) * 2014-08-27 2016-09-27 Общество с ограниченной ответственностью "ТатАСУ" Способ нейросетевого анализа состояния телемеханизированных нефтепромысловых объектов
RU2598785C2 (ru) * 2014-11-05 2016-09-27 Общество с ограниченной ответственностью "ТатАСУ" Система статистического и нейросетевого анализа данных телеметрии нефтепромысловых объектов
US10227869B2 (en) 2015-04-09 2019-03-12 Joy Global Underground Mining Llc System and method of detecting dull and worn cutter bits
US9920624B2 (en) 2015-04-09 2018-03-20 Joy Mm Delaware, Inc. System and method of detecting dull and worn cutter bits
NO345487B1 (en) * 2015-05-12 2021-03-01 Halliburton Energy Services Inc Enhancing oilfield operations with cognitive computing
US11168558B2 (en) 2015-05-12 2021-11-09 Halliburton Energy Services, Inc. Enhancing oilfield operations with cognitive computing
US10597995B2 (en) 2015-06-26 2020-03-24 Landmark Graphics Corporation Visualization of quantitative drilling operations data related to a stuck pipe event
GB2554304A (en) * 2015-06-26 2018-03-28 Landmark Graphics Corp Visualization of quantitative drilling operations data related to a stuck pipe event
GB2554304B (en) * 2015-06-26 2021-05-26 Landmark Graphics Corp Visualization of quantitative drilling operations data related to a stuck pipe event
WO2016209271A1 (fr) * 2015-06-26 2016-12-29 Landmark Graphics Corporation Visualisation de données quantitatives d'opérations de forage concernant un événement de tube coincé
WO2017070069A1 (fr) * 2015-10-23 2017-04-27 Schlumberger Technology Corporation Utilisation de données de forage pour conditionnements techniques
US10883346B2 (en) 2015-12-18 2021-01-05 Schlumberger Technology Corporation Method of performing a perforation using selective stress logging
EA038033B1 (ru) * 2016-12-07 2021-06-25 Сейфкик Америкас Ллс Автоматизированное бурение на основе модели
US10781681B2 (en) 2016-12-07 2020-09-22 Safekick Americas Llc Automated model based drilling
WO2018106346A1 (fr) * 2016-12-07 2018-06-14 Safekick Americas Llc Forage basé sur un modèle automatisé
CN106545327A (zh) * 2016-12-09 2017-03-29 北京四利通控制技术股份有限公司 智能司钻钻机控制系统
CN106545327B (zh) * 2016-12-09 2017-11-28 北京四利通控制技术股份有限公司 智能司钻钻机控制系统
US10773352B2 (en) 2017-06-05 2020-09-15 Joy Global Underground Mining Llc System and method for determining efficiency of an industrial machine
US11384612B2 (en) 2017-07-11 2022-07-12 Equinor Energy As Method and system for monitoring influx and loss events in a wellbore
CN111479982A (zh) * 2017-11-15 2020-07-31 吉奥奎斯特系统公司 具有滤波器的现场操作系统
CN111479982B (zh) * 2017-11-15 2023-01-17 吉奥奎斯特系统公司 具有滤波器的现场操作系统
US11255180B2 (en) * 2017-12-22 2022-02-22 Halliburton Energy Services, Inc. Robust early kick detection using real time drilling
CN112004989A (zh) * 2018-03-07 2020-11-27 贝克休斯控股有限责任公司 用于显示钻地工具的可靠性的钻地工具监测系统及相关方法
CN112004989B (zh) * 2018-03-07 2024-05-10 贝克休斯控股有限责任公司 用于显示钻地工具的可靠性的钻地工具监测系统及相关方法
CN108612518A (zh) * 2018-04-20 2018-10-02 重庆地质矿产研究院 一种煤层气井径向微井眼钻孔、水力压裂参数的确定方法
US11066917B2 (en) * 2018-05-10 2021-07-20 Baker Hughes Holdings Llc Earth-boring tool rate of penetration and wear prediction system and related methods
CN109062248A (zh) * 2018-08-02 2018-12-21 西北工业大学 基于自组织网络的空间非合作目标姿轨一体化参数估计方法
CN112714820A (zh) * 2018-09-14 2021-04-27 埃尼股份公司 通过钻井设备估算待钻的地质构造中孔隙压力值的方法
US11852010B2 (en) 2018-09-14 2023-12-26 Eni S.P.A. Method for estimating a pore pressure value in geological formations to be drilled by a drilling apparatus
CN112714820B (zh) * 2018-09-14 2024-06-07 埃尼股份公司 通过钻井设备估算待钻的地质构造中孔隙压力值的方法
US11542788B2 (en) * 2019-03-01 2023-01-03 Brigade Energy Services Drilling out frac plugs
CN110261906A (zh) * 2019-06-27 2019-09-20 中国石油化工股份有限公司 多参数岩性信息融合储层预测方法
CN110261906B (zh) * 2019-06-27 2021-06-22 中国石油化工股份有限公司 多参数岩性信息融合储层预测方法
CN112412390A (zh) * 2019-08-22 2021-02-26 中国石油化工股份有限公司 基于深度学习模型评价固井第二界面的方法及装置
CN112412390B (zh) * 2019-08-22 2022-09-02 中国石油化工股份有限公司 基于深度学习模型评价固井第二界面的方法及装置
CN110941866B (zh) * 2019-12-06 2023-04-07 中国石油集团川庆钻探工程有限公司 一种基于固井大数据的环空水泥浆界面设计方法
CN110941866A (zh) * 2019-12-06 2020-03-31 中国石油集团川庆钻探工程有限公司 一种基于固井大数据的环空水泥浆界面设计方法
CN113378344A (zh) * 2020-02-25 2021-09-10 中国石油化工股份有限公司 一种用于优化钻井起下钻速度阈值的方法及系统
US11459882B2 (en) 2020-03-20 2022-10-04 Saudi Arabian Oil Company Systems and methods for the determination of lithology porosity from surface drilling parameters
US11131184B1 (en) 2020-04-29 2021-09-28 Saudi Arabian Oil Company Method and system for determining a drilling hazard condition using well logs
WO2021221700A1 (fr) * 2020-04-29 2021-11-04 Saudi Arabian Oil Company Procédé et système pour déterminer un état de danger de forage à l'aide de diagraphies de puits
US11078785B1 (en) 2020-06-17 2021-08-03 Saudi Arabian Oil Company Real-time well drilling evaluation systems and methods
CN114169656A (zh) * 2020-09-11 2022-03-11 中国石油化工股份有限公司 一种基于邻井历史数据的钻井卡钻风险预警方法及系统
CN114427430B (zh) * 2020-09-21 2024-05-07 中国石油化工股份有限公司 一种多井实时协同的钻井参数优化方法与系统
CN114427430A (zh) * 2020-09-21 2022-05-03 中国石油化工股份有限公司 一种多井实时协同的钻井参数优化方法与系统
CN112418489A (zh) * 2020-11-05 2021-02-26 中国石油天然气股份有限公司 一种天然气井修井暂堵剂用量神经网络预测方法
CN112836406B (zh) * 2021-01-13 2022-06-14 中国地质大学（武汉） 一种考虑不确定阻尼系数的钻进过程钻压建模方法
CN112836406A (zh) * 2021-01-13 2021-05-25 中国地质大学（武汉） 一种考虑不确定阻尼系数的钻进过程钻压建模方法
CN113065664A (zh) * 2021-02-20 2021-07-02 北京三一智造科技有限公司 一种地质层识别模型构建方法和地质层识别方法
CN113065664B (zh) * 2021-02-20 2024-05-14 北京三一智造科技有限公司 一种地质层识别模型构建方法和地质层识别方法
CN113073959A (zh) * 2021-02-24 2021-07-06 西南石油大学 一种基于模糊决策树的钻井漏失机理预测及辅助决策方法
GB2617039B (en) * 2021-03-03 2024-12-04 Landmark Graphics Corp Drill bit wear and behaviour analysis and correlation
GB2617039A (en) * 2021-03-03 2023-09-27 Landmark Graphics Corp Drill bit wear and behaviour analysis and correlation
WO2022186893A1 (fr) * 2021-03-03 2022-09-09 Landmark Graphics Corporation Analyse et corrélation d'usure et de comportement de trépan de forage
US11578583B2 (en) 2021-03-03 2023-02-14 Landmark Graphics Corporation Drill bit wear and behavior analysis and correlation
CN113310829A (zh) * 2021-04-14 2021-08-27 西南石油大学 一种冲击方式下岩石的可钻性测试装置及实验方法
CN113310829B (zh) * 2021-04-14 2022-10-25 西南石油大学 一种冲击方式下岩石的可钻性测试装置及实验方法
US12098631B2 (en) 2021-04-23 2024-09-24 Landmark Graphics Corporation Process-mining software for generating a process flow for forming a wellbore
WO2022225537A1 (fr) * 2021-04-23 2022-10-27 Landmark Graphics Corporation Logiciel d'exploitation minière permettant de générer un flux de traitement pour former un puits de forage
GB2618971B (en) * 2021-04-23 2024-12-25 Landmark Graphics Corp Process-mining software for generating a process flow for forming a wellbore
GB2618971A (en) * 2021-04-23 2023-11-22 Landmark Graphics Corp Process-mining software for generating a process flow for forming a wellbore
CN113586028A (zh) * 2021-07-21 2021-11-02 宁夏天地奔牛实业集团有限公司 基于数字孪生的反井钻机扩孔刀盘的智能监控系统
CN113586028B (zh) * 2021-07-21 2024-03-29 太原理工大学 基于数字孪生的反井钻机扩孔刀盘的智能监控系统
CN113792936A (zh) * 2021-09-28 2021-12-14 中海石油(中国)有限公司 一种随钻岩性智能识别方法、系统、设备和存储介质
CN114116213A (zh) * 2021-11-24 2022-03-01 上海派拉软件股份有限公司 分布式系统中超额负载估计方法、装置、设备及存储介质
CN114607354A (zh) * 2021-12-14 2022-06-10 中国石油天然气集团有限公司 一种预报钻井复杂情况的方法
CN114526052B (zh) * 2021-12-31 2023-09-19 中国石油天然气集团有限公司 一种钻完井工程风险预测方法及装置
CN114526052A (zh) * 2021-12-31 2022-05-24 中国石油天然气集团有限公司 一种钻完井工程风险预测方法及装置
CN114357887A (zh) * 2022-01-06 2022-04-15 中国地质大学（北京） 一种基于bp神经网络的复杂井况钻前泥浆漏失预测方法
CN114357887B (zh) * 2022-01-06 2025-02-28 中国地质大学（北京） 一种基于bp神经网络的复杂井况钻前泥浆漏失预测方法
US12291955B2 (en) 2022-03-18 2025-05-06 Saudi Arabian Oil Company Real-time model of rig and bit hydraulics efficiency
CN114997485A (zh) * 2022-05-26 2022-09-02 中国石油天然气集团有限公司 溢流工况预测模型训练方法、装置及溢流工况预测方法
CN114997485B (zh) * 2022-05-26 2023-10-03 中国石油天然气集团有限公司 溢流工况预测模型训练方法、装置及溢流工况预测方法
CN115059448A (zh) * 2022-06-01 2022-09-16 中国石油大学(华东) 一种基于深度学习算法的地层压力监测方法
CN114912301A (zh) * 2022-07-18 2022-08-16 中国空气动力研究与发展中心低速空气动力研究所 一种低速风洞全机模型测力试验数据处理与修正系统
WO2024155870A1 (fr) * 2023-01-19 2024-07-25 Baker Hughes Oilfield Operations Llc Détection probabiliste de scénarios de forage
CN116181261A (zh) * 2023-03-29 2023-05-30 安徽理工大学 一种煤矿立井钻井法凿井钻头吸渣口布置的确定方法
CN116181261B (zh) * 2023-03-29 2024-03-19 安徽理工大学 一种煤矿立井钻井法凿井钻头吸渣口布置的确定方法
CN117034665B (zh) * 2023-10-10 2024-01-09 中国石油大学(华东) 一种基于钻井施工预演的井身结构设计方法
CN117034665A (zh) * 2023-10-10 2023-11-10 中国石油大学(华东) 一种基于钻井施工预演的井身结构设计方法
CN118309368A (zh) * 2024-06-07 2024-07-09 陕西延长石油矿业有限责任公司 一种煤矿厚基岩地层一钻完井智能建井方法
CN118309368B (zh) * 2024-06-07 2024-10-18 陕西延长石油矿业有限责任公司 一种煤矿厚基岩地层一钻完井智能建井方法

Also Published As

Publication number Publication date
CA3064241A1 (fr) 2014-04-30
MX2015005568A (es) 2016-02-03
MX360775B (es) 2018-11-16
CA3064241C (fr) 2022-12-13
CA2794094A1 (fr) 2014-04-30
CA2794094C (fr) 2020-02-18

Similar Documents

Publication Publication Date Title
US9970266B2 (en) Methods and systems for improved drilling operations using real-time and historical drilling data
CA3064241C (fr) Procedes et systemes pour operations de forage ameliorees utilisant des donnees de forage historiques et en temps reel
US12050981B2 (en) Petroleum reservoir behavior prediction using a proxy flow model
US20180171774A1 (en) Drillstring sticking management framework
US7878268B2 (en) Oilfield well planning and operation
US20130341093A1 (en) Drilling risk avoidance
US8229880B2 (en) Evaluation of acid fracturing treatments in an oilfield
US20120203525A1 (en) Three-dimensional modeling of parameters for oilfield drilling
US10920561B2 (en) Drilling assessment system
EP3619399B1 (fr) Structure d'ensemble de train de tiges de forage
EP2788579A2 (fr) Console de contrôle géologique
US20170002630A1 (en) Method of performing additional oilfield operations on existing wells
GB2458356A (en) Oilfield well planning and operation
CN116615600A (zh) 马达效率和退化解释系统
Zhdaneev et al. Predictive systems for the well drilling operations
CA3107639A1 (fr) Diagraphie synthetique en temps reel permettant l'optimisation du forage, de la direction et de la stimulation
EP3250922B1 (fr) Procédé et système informatiques pour traiter et transformer des données collectées pour améliorer la productivité de forage
US9482088B2 (en) Mean regression function for permeability
US20240426178A1 (en) Orchestration framework to compute dynamic risk and impact maps
US11640012B2 (en) Virtual high-density well survey
US12252974B2 (en) Equipment health monitor

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: 13850354

Country of ref document: EP

Kind code of ref document: A1

NENP Non-entry into the national phase

Ref country code: DE

WWE Wipo information: entry into national phase

Ref document number: MX/A/2015/005568

Country of ref document: MX

122 Ep: pct application non-entry in european phase

Ref document number: 13850354

Country of ref document: EP

Kind code of ref document: A1

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