EARLY DETECTION OF KICKS

[0001] This application claims the priority of U.S. Utility Application Serial No. 15/045,362, filed February 17, 2016 and U.S. Provisional Application Serial No. 62 117,061, filed February 17, 2015, and hereby incorporates that application by reference for all purposes as if set forth verbatim herein.

[0002] This section of this document introduces various information that may be related to or provide context for some aspects of the technique described herein and/or claimed below. It provides background information to facilitate a belter understanding of that which is disclosed herein. This is therefore a discussion of "related" art. That such art is related in no way implies that it is ako ''prior" art. The discussion in this section is to be read in this light, and not as admissions of prior ait.

[0003] The efforts of the oil and gas industry to discover and bring into production new or additional hydrocarbon deposits has led to ever more sophisticated and demanding technical environments. This sophistication and demand is reflected in the costs of the endeavor. ( >nc part of this evolution in the industry responsive to these concerns is unproved techniques for monitoring and managing phenomena such as "kicks". Kicks arc unplanned subsurface fluid or gas flow influxes from the geological reservoir into the wellbore during oil and gas drilling, tripping, and completion or intervention operations. Drilling mud, completion fluids, and drilling cement serve as barriers against pressurized hydrocarbons in the reservoir and keep them sealed in the reservoir until production commences. In the event that wellbore fluid pressures become less than that of an exposed subsurface formation, a kick may occur. Drilling operations and unanticipated high pressure gas pockets in porous rock formations can lead to pressure imbalances between wellbore fluids and reservoir fluids, causing gas influx into the wellbore or loss of drilling mud into the reservoir.

[0004] One issue in kick detection is that the conditions indicating that a kick has occurred arc typically not readily detectable by the human eye. A fair portion of this fact is that many of the conditions used to detect or predict a kick are downhole, and so are not readily discernible directly to the human eye. Some factors may be deduced at the surface but the delay caused by the change in conditions propagating to the surface works against the need for a quick detection. Accordingly, the industry typically instruments a string downhole as well as at the surface to monitor condition which might indicate that a kick has occurred.

[000SJ However, even with automated monitoring systems, many techniques for detecting and managing kicks suffer from a number of drawbacks. It is not uncommon for them to rely on lagging rather than leading indicators, which can delay an otherwise timely response. They arc also subject unpredictable human error. For example, many of the measured parameters may be correlated to the unplanned influx of formation fluids into the wellbore without being indicative, or the operator may miss the significance of a piece or stream of information.

[0006] The presently disclosed technique is directed to resolving, or at least reducing, one or all of the problems mentioned above. As set forth above, several techniques for monitoring well conditions and detecting kicks arc known to the art and arc competent for their intended purposes. The art, however, is always receptive to improvements or alternative means, methods, and configurations. Therefore the art will consequently well receive the technique described herein.

[0007] The presently disclosed technique presents to the art a well monitoring system particularly useful in detecting kicks m the well. ΊΊκ well monitoring system comprises a well, a well system, and a computing apparatus. The well defines a wellbore and the well system includes at least one sensor measuring at least one well condition. The computing apparatus includes a processor, storage, a bus system over which the processor communicates with the storage, a data structure residing in the storage, and a well monitoring software component residing in the storage. The data structure stores real-time data acquired by the sensor.

[0008] The well monitoring software component when executed by the processor over the bus system, performs a method to detect a kick in a well. The method comprises: storing a set of real-time data from a measurement of a well condition by the sensor, the measurements being correlative to an unplanned fluid influx into the well; modeling the operation of the well with a physics-based, state space model of a well system of the well to obtain an estimate of the well condition; accessing the stored real-time data set and applying the accessed real-time data set and the estimate to a probabilistic estimator to yield a probability of an occurrence of a kick and a confidence measure for the probability.

[0009] Other aspects of the presently disclosed technique include a computer-implemented method to detect a kick, a non-transitory program storage medium encoded with instructions that, when executed, perform such a computer-implemented method, and a computing apparatus programmed to perform such a method.

[0010] The above presents a simplified summary of the invention in order to provide a basic understanding of some aspects of the subject matter disclosed herein and claimed below. This summary is not an exhaustive overview of that which is claimed. It is not intended to identify key or critical elements of the claimed subject matter or to delineate its scope. The sole purpose of this summary is to present some concepts in a simplified form as a prelude to the more detailed description that is discussed later.

[0011] The claimed subject matter may be understood by reference to the following description taken in conjunction with the accompanying drawings, in which like reference numerals identify like elements, and in which: [0012] Figure 1 depicts a drilling operation in which one particular embodiment of the presently disclosed technique is practiced in a partially sectioned, plan view.

[0013] Figure 2 presents one particular embodiment of a method practiced in accordance with the technique disclosed herein.

[0014] Figure 3 conceptually illustrates selected portions of the hardware and software architecture of a computing apparatus such as may be employed in some aspects of the present invention.

[0015] Figure 4 graphically illustrates the performance of the method of the disclosed technique in one particular embodiment.

[0016] Figure 5-Figurc 6 convey how combining multiple models/predictions of the same quantity gives significantly reduced uncertainty in the estimated value.

[0017] Figure 7 depicts selected portions of a well system for purposes of illustrating a particular model thereof.

[0018] Figure 8 is a bond graph model from which process and measurement equations may be obtained for the wcUborc and well reservoir hydraulics of the well system of Figure 7.

[0019] Figure 9 illustrates the efficacy of the presently disclosed technique.

[0020] While the invention is susceptible to various modifications and alternative forms, the drawings illustrate specific embodiments herein described in detail by way of example. It should be understood, however, that the description herein of specific embodiments is not intended to limit the invention to the particular forms disclosed, but on the contrary, the intention is to cover all modifications, equivalents, and alternatives falling within the spirit and scope of the invention as defined by the appended claims.

[0021] Illustrative embodiments of the subject matter claimed below will now be disclosed. In the interest of clarity, not all features of an actual implementation are described in this specification. It will be appreciated that in the development of any such actual embodiment, numerous implementation-specific decisions must be made to achieve the developers ^* specific goals, such as compliance with system-related and business-related constraints, which will vary from one implementation to another. Moreover, it will be appreciated that such a development effort, even if complex and time-consuming, would be a routine undertaking for those of ordinary skill in the art having the benefit of this disclosure.

[0022] The technique disclosed herein and claimed below employs a cyber-physical approach to the detection, monitoring, and managing of kick in wells. For present purposes, a "cyber-physical" technique is one in which a model of the well system for the well is coupled to the well system in operation. The model and well system arc coupled in that the model incorporates system knowledge and physical knowledge of the well system developed during the well system's design and implementation. The model then resides and operates in a virtual environment to model the well system's operation in real time while the well system is operating based on information acquired by interacting with the well system through the coupling. In this sense, the model "mirrors" the operation of the well system and can continuously track and provide information regarding the well system's operation that is not always amenable to direct observation. This informatton can then be analyzed to determine whether a kick is actually occurring or even is imminent before it happens.

[0023] Unlike conventional practice, the cybcr-physical approach combines multiple measurements by linking the measurements of the operation with the physics of the operation. This provides for natural scaling of the measurements relative to each to other for making predictions of output variables. It ako provides for natural filtering or smoothing of the estimate. Conventional practice, on the other hand, relies on ad hoc smoothing or averaging of the measured data. The presently disclosed technique furthermore does not just trigger on a pattern in the data but provides a quantifiable estimate of a kick with quantifiable uncertainty.

[0024] This technique uses multiple real-time measurements of conditions in the well environment that can be linked, or correlated, to kick. In a drilling context, commonly available variables include mud pit volume, return flow, input flow, standpipc pressure, drilled depth, hook load, gas content, and others. These measurements arc combined with physics-based, state space models of the operation. It is applicable in a wide variety of wells including both on-shore and off-shore in which there arc a variety of types and accuracies of measurements and physical configurations.

[0025] One principle of the technique is that combining multiple measurements of even very noisy and uncertain measurements reduces the uncertainty in estimated values provided by the models. In some embodiments, these measurements arc then combined with estimates made by a physics-based state space model to produce even more accurate estimated values representing a probability. For example, a typical output estimated value of interest in early kick detection is amount (mass or moles) of hydrocarbon influx. This combination uses measurements that are numerically quantified by the states of the model. In order to combine measurements and model estimates this approach also quantifies the uncertainties in the measurements and the modcL Model uncertainty includes uncertainty in both model inputs and in model parameters.

[0026Ί Once this has been done a real-time probabilistic estimator is then used to estimate the states of the model, which give probabilistic estimates or, a probability of outputs such as hydrocarbon influx. The estimator gives not only a most likely value but also the uncertainty of the value. These procedures allow estimation of values that cannot be easily measured.

[0027] The physics of the model allow construction of a relationship between measured quantities and kick. In one embodiment, a simple incompressible hydraulic model allows us to link the pump pressure to the bottom hole pressure and with a model of the formation permeability. This allows a prediction of influx rate. [0028] Higher fidelity models, which predict variables with more accuracy, can also be used. There is a trade-oil between higher fidelity and simulation time. Some embodiments may seek prediction in realtime. If the model runs slower than real-time there arc at least two remedies. One is to develop a lower order model that captures the important physics of the high fidelity model. The second is to use modern computer architecture and hardware that can run parallel processes. These systems are becoming available at very low cost. Λ graphics processing unit is an example of some this new computer hardware.

[0029] The presently disclosed technique will now be described with reference to the attached figures. Various structures, systems and devices are schematically depicted in the drawings for purposes of explanation only and so as to not obscure the present invention with details that arc well known to those skilled in the art. Nevertheless, the attached drawings are included to describe and explain illustrative examples of the present invention.

[0030] Turning now to Figure 1. a drilling operation 100 includes a hydrocarbon well 103 drilled through the earth's surface 106 and into and through a subterranean formation 109 surrounding the hydrocarbon well 103. The hydrocarbon well 103 includes a string 112 shown run into the wellbore 115. The wellbore 115 is also filled with drilling fluids 118 in a manner known to the art for purposes well known to the art The drilling fluids 118 may be any kind of drilling fluid known to the art and suitable for the purpose for which it is introduced. For example, the drilling fluids 118 may be a drilling "mud" introduced to maintain the hydrostatic pressure of the well 103 at a desired level. The wellbore 1 IS passes through a portion of the formation 109 containing deposits of formation fluids 121, such as water or brine, or a hydrocarbon such as natural gas or petroleum. The identity of the formation fluids 121 is not material to the practice of the technique disclosed and claimed herein although it may be significant in a given embodiment

[0031] Those skilled in the art having the benefit of this disclosure will appreciate that the illustration in Figure 1 is highly idealized. For example, the subterranean formation 109 is illustrated in a manner from which one might infer it is of a homogeneous composition. Those in the art will understand that this is unlikely to be the case and that the subterranean formation 109 will contain many strata (not shown) of varying geophysical characteristics. Similarly, mere may be many deposits of formation fluids 121 in the subterranean formation 109 or, in some circumstances, none. These and other such variations which have been suppressed for the sake of clarity will be readily recognized by those skilled in the art.

[0032] The welbore 1 IS is "cased", as is evident from the casing 116. Most wells will be cased as shown. However, the presently disclosed technique is not limited to cased wclb. It may also be applied to what are known as "open holes", or those wells whose wellbores remain uncased or from which previously installed casing has been removed. It may also be applied to cased wells that are open at the bottom. [0433] The drill string 112 includes, for example, a bottom hole assembly 124 comprised of a bit 127, data and crossover sub 130, and sensor apparatus 133. The drill string 112 also includes other conventional string components that arc not indicated such as tools, jars, stabilizers, drill collars, and drill pipe. The constitution, assembly, and deployment of the drill string 112 may accord with conventional practice using principles and techniques well known to those in the art.

[0034] Those in the art might infer from the presence of the bottom hole assembly 124 that the operation depicted in Figure 1 is a drilling operation. However, the presently disclosed technique is not necessarily limited to use in drilling operations. The presently disclosed technique may be used in practically any phase of well operations in which kick is of interest

[0035] It is well known to instrument the drill string 1 12 with a variety of sensors 136 (only one indicated) to monitor conditions throughout the wellbore 115. For example, the data and crossover sub 130 may house an acceleromeler (not otherwise shown) useful for gathering real-time data from the bottom of the wellbore 115. For example, the accelerometer can give a quantitative measure of bit vibration. Many types of data sources may and typically will be included. Exemplary measurements that may be of interest include hole temperature; the pressure, salinity and pH of the drilling mud; the magnetic declination and horizontal declination of the bottom-hole assembly; seismic look-ahead information about the surrounding formation; electrical resistivity of the formation; pore pressure of the formation; gamma ray characterization of the formation, and so forth.

[0036] Any given embodiment will typically be more interested in some quantities than in others. In particular, as is described further below, the inputs to the models should be correlated in some way to kick. Thus, quantities such as mud pit volume, return flow, input flow, slandpipe pressure, drilled depth, hook load, gas content, etc. will be of particular interest

[0037] To this end, a variety of instrumented tools 139 (only one indicated) for gathering information regarding downhole drilling conditions will be included in the drill string 112. However, not all sensors 136 will necessarily be disposed on or in an instrumented tool 139. The sensors 136 may be disposed anywhere throughout the drill string 112 in any manner suitable to those skilled in the art that is known to the art.

[0038] Note that the embodiments illustrated herein are intended for use with quantities that are already sensed and whose measurements are already available through well monitoring software. The technique is therefore suitable for retrofit onto existing wells. However, there is no need to limit other embodiments to those quantities that arc already sensed and whose measurements arc available. Some embodiments may contemplate the use of quantities not typically sensed such that additional sensors may be added to the string over and above those that are conventionally used [0409] Information sensed by the sensors 136 is communicated back to the surface 106 where it is collected. In the illustrated embodiment the information is communicated electronically over a line 142 to a computing apparatus 14S. The sensed information is converted into digital data at the sensor 136 and electronically transmitted over the line 142. In some embodiments, the data transmission is interleaved on the line 142. Some embodiments may employ more than one line 142 to avoid or alleviate operational constraints imposed by using a single line 142. Some embodiments may even transmit some or all of the information wirelessly. There are still other techniques known to the art by which the sensed information may be communicated to the eurface. .Any such technique known to the art suitable for the purpose may be employed in alternative embodiments.

[0040J It is abo known to instrument surface operations. For example, in Figure 1 there is conceptually shown a mud pit 141 from which the mud 118 is pumped into the wellbore 115 and to which mud 118 is returned from the wellbore 115. Sensors 137 measure various aspects of the well 103*8 operation with respect to the mud pit 141 such as mud volume in the mud pit 141 and the rate of flow out of the mud pit 141. The measurements arc then also communicated to the computing apparatus 14S over a line not shown in Figure 1. Those in the art will appreciate that many aspects of surface operations arc monitored in this fashion and that the mud pit operations arc merely illustrative of surface operations in general.

[0041] Figure 2 illustrates a method 200 in accordance with one aspect of the presently disclosed technique. The method 200 is, in this particular embodiment, performed at least in part by the computing apparatus 14S. Λ brief description of those portions of the computing apparatus 145 pertinent to that performance shall therefore now be discussed before returning to Figure 2.

[0042] Figure 3 shows selected portions of the hardware and software architecture of one particular embodiment of the computing apparatus 145. The computing apparatus 145 includes in this embodiment a processor 300 communicating with storage 303 over a bus system 306. The storage 303 may include a hard disk and/or random access memory ("RAM") and/or removable storage such as a floppy magnetic disk 309 and an optical disk 312.

[0043] The processor 300 may be any suitable processor known to the art. Those in the art will appreciate that some types of processors will be preferred in various embodiments depending on familiar implementation specific details. For example, some processors arc more powerful and process faster so that they may be more preferred where large amounts of data are to be processed in a short period of time. On the other hand, some processors consume more power and available power may be severely limited in some embodiments, how power consumption processors may therefore be preferred in those embodiments.

[0044] These kinds of factors are commonly encountered in the design process and will be highly implementation specific. Because of their ubiquity in the art, such factors will be easily reconciled by those skilled in the art having the benefit of this disclosure. Those in the art having the benefit of this disclosure will therefore appreciate that the processor 300 may be a micro-controller, a controller, a microprocessor, a processor set. or an appropriately programmed application specific integrated circuit ("ASIC") or field programmable gate array ("FFGA"). Some embodiments may even use some combination of these processor types.

[00451 As with the processor 300, implementation specific design constraints may influence the design of the storage 303 in any particular embodiment For example, it is well known that certain types of types of memory (eg., cache) have much faster access times than other types (e.g., disk memory). Some types of memory will also consume more power than others. Some embodiments may wish to only temporarily buffer acquired data whereas others may wish to store H for a more prolonged period. As with the processor 300, these kinds of factors are commonplace in the design process and those skilled in the art having the benefit of this disclosure will be able to readily balance them in light of their implementation specific design constraints.

[0046] The storage 303 is encoded with a data structure 3 IS in which the data 318 received from the one or more sensors 136 over the line 142 may be buffered or otherwise stored. As is apparent from the discussion above, the data 318 comprises information regarding the (trilling conditions in the wellbore 1 IS, the drilling fluids 118, the wellbore 115, and the surrounding formation 109. The data 318 therefore represents tangible, real world object namely, the wellbore 115, drilling fluids 118, and the formation 109. The data structure 315 may be any suitable data structure known to the art, such as a buffer, a string, a linked list a database, etc. The data 318 may be buffered or it may be stored more long term even archived— depending on the embodiment The data structure 315 may even be a composite of constituent data structures (not shown) if, for example, it is desired to have a separate data structure for each set of data generated by different sensors 136. The disclosed technique admits wide variation in the implementation of the data structure 315.

[0047] A well monitoring software component 321 that performs the software-implemented method described below is also encoded on the storage 303. The well monitoring software component 321 may be coded in any suitable manner known to the art. The well monitoring software component 321 is, in mis particular embodiment, an application. Note, however, that there is no requirement that this functionality be implemented in an application. For example, the well monitoring software component 321 may be implemented in some other kind of software component such as a daemon or utility. 'Ilie functionality of the well monitoring software component 321 also need not be contained in a single software component and may be separated into two or more components. The functionality may be aggregated into a single component or distributed across more than two components.

[0048] The storage 303 is also encoded with one or more physics-based state space model(s) 324 of the well system and a probabilistic estimator 327. The model(s) 324 and probabilistic estimator 327 are used by the well monitoring software component 321 as described below to implement the software implemented aspects of the presently disclosed technique. The modeKs) 324 and the probabilistic estimator 327 arc also described in more detail below. Just as the well monitoring software component 321 may be implemented in wide variation across embodiments, so may the modcl(s) 324 and the probabilistic estimator 327. For example, rather than being stand-alone components called by the well monitoring software component 321, either one or both of the modeKs) 324 or the probabilistic estimator 327 may be incorporated into the well monitoring software component 321. ()r, they may be separate from the well monitoring software component 321 but combined with each other into another component

[0049] In particular, the modeKs) 324 model the well system of the well 100 that arc pertinent to a kick. For example, in the detection of kick, the pertinent parts of the well system that should be modeled include the hydraulics, the mechanics of the system, and the formation. They hydraulics would include information such as the physical characteristics (e.g.. weight, temperature, pi L, gas content), the volume, and the rate of circulation of the drilling fluids as well as return flow and input flow. The mechanics of the system includes such things as the mud pit volume, the drilled depth of the wellborc, the cased diameter of the wellborc, the rate of penetration, standpipc pressure, the hook load, and other information pertaining to the physical characteristics of the wellborc. The formation would include geophysical characteristics such as those listed in Table 3 below. The various part of the well system may be separately modeled and then interfaced or all integrated into a single model. Thus, the model(s) 324 may be a single model or a plurality of models.

[0050] The storage 303 is also encoded with an operating system 330 and user interface software 333. The user interface software 333, in conjunction with a display 336, implements a user interface 339. The user interface 339 may include peripheral 1 0 devices such as a keypad or keyboard 342, a mouse 34S, or a joystick 348. The processor 300 runs under the control of the operating system 330, which may be practically any operating system known to the art. The well monitoring software component 321 is invoked by the operating system 330 upon power up, reset, or both, depending on the implementation of the operating system 330. The application 465, when invoked, performs the method of the present invention. The user may also invoke the monitoring software component 321 in conventional fashion through the user interface 339 in some embodiments.

[0051] One aspect of the presently disclosed technique that separates it from many computing applications is the computationally intensive nature of the tasks to which it is assigned. The software processes voluminous real-time data through a model of the well system and quick resolution and reporting arc typical objectives. It is unlikely that a general purpose computing apparatus will meet these performance considerations. The process 300 should be implemented as a processor set that will include some degree of parallel processing. The storage 303 should be designed for rapid read'write operations, which favors RAM and cache of removable storage. The model(s) 327 should be designed or selected with a suitable balance of resolution and speed. These and other design considerations mitigate for a computing environment that is much more computationally robust than in a general purpose computing environment.

[0052] As is evident from the discussion above, some portions of the detailed descriptions herein are presented in terms of a software implemented process involving symbolic representations of operations on data bits within a memory in a computing system or a computing device. These descriptions and representations are the means used by those in the art to most effectively convey the substance of their work to others skilled in the art. The process and operation require physical manipulations of physical quantities that will physically transform the particular machine or system on which the manipulations arc performed or on which the results are stored. Usually, though not necessarily, these quantities take the form of electrical, magnetic, or optical signals capable of being stored, transferred, combined, compared, and otherwise manipulated. It has proven convenient at times, principally for reasons of common usage, to refer to these signals as bits, values, elements, symbols, characters, terms, numbers, or the like.

[0053] It should be borne in mind, however, that all of these and similar terms arc to be associated with the appropriate physical quantities and arc merely convenient labels applied to these quantities. Unless specifically stated or otherwise as may be apparent, throughout the present disclosure, these descriptions refer to the action and processes of an electronic device, that manipulates and transforms data represented as physical (electronic magnetic, or optical) quantities within some electronic device's storage into other data similarly represented as physical quantities within the storage, or in transmission or display devices. Exemplary of the terms denoting such a description arc, without limitation, the terms "processing." "computing," "calculating," "determining," "displaying," and the like.

[0054] Furthermore, the execution of the software's functionality transforms the computing apparatus on which it is performed. For example, acquisition of data will physically alter the content of the storage, as win subsequent processing of that data. The physical alteration is a "physical transformation" in that it changes the physical state of the storage for the computing apparatus.

[0055] Note also that the software implemented aspects of the invention are typically encoded on some form of non-transitory program storage medium or implemented over some type of transmission medium. The program storage medium may be magnetic {e.g., a floppy disk or a hard drive) or optical {e.g.. a compact disk read only memory, or "O) ROM"), and may be read only or random access. Similarly, the transmission medium may be twisted wire pairs, coaxial cable, optical fiber, or some other suitable transmission medium known to the art. 'fhc invention is not limited by these aspects of any given implementation.

[0056] Another thing that will typically separate the computing aspects of the technique from general purpose computing is the environment of the well system. The computing apparatus 145 nominally appears as a work station in Figure 1. Those in the ait having the benefit of this disclosure will appreciate that many, if not most, rigs are equipped with computers of some kind. These computers are hardened against vibration, dust, and other environmental conditions encountered in a drilling environment but not in more sedate office and residential environments. Some of these computers may be rack mounted rather than a stand-alone workstation. The computing apparatus 14S may be, in some embodiments, a computer already on a rig retrofitted to implement the technique disclosed herein. Alternatively, rigs may be equipped with new computers not only programmed to implement the present technique but also finished out in accordance with practices well known to the art to adapt them to the drilling environment,

[0057] There also is no theoretical or operational requirement that the computing apparatus 145 he implemented in a single, unitary, integrated package. For example, some embodiments might choose to store the data 318 locally while hosting the well monitoring software component 321 offstte at another location. In these embodiments, the data 318 can be accessed by the well monitoring software component 321 for analysis remote .from the location at which it is collected. Information output by the well monitoring software component 321 can then be utilized at that remote location, or locally at the location where it is collected, or at yet a third location.

[0058] Referring now to both Figure 2 and Figure 3, the method 200 is performed by the well monitoring software component 321 when invoked by the processor 300 over the bus system 306. The method 200 assumes that well monitoring through, for example, the sensors 136 and 137 is ongoing in a manner known to the ait and mat the sensed measurements arc being stored in the data structure 315 as data ("DATA"). The data is therefore real-time data. Note that some embodiments may also employ near realtime or even archived data in addition to real-time data.

[0059] The method 200 begins, in this particular embodiment, with the well monitoring software component 321 storing (at 210) a set of real-time data from a measurement of a well condition acquired during the operation of the well, the measurements being correlative to an unplanned fluid influx into the well 103. The measured well condition may be a downholc condition or a surface condition. Typically, a plurality of measured conditions is used and that plurality will include both downholc and surface conditions. The conditions themselves, as well as their measurements, may be independent of one another or they may be related. Again, most embodiments will typically include both independent and related measurements.

[0060] The well monitoring software component 321 also models (at 220) the operation of the well 103 with a physics-based, state space model 324 of well system of the well 103 to obtain an estimate of the well condition, the model being cyber-physically coupled to the well system. It also accesses (at 225) the stored real-time data set The accessing (at 225) and the modeling (at 220) may be performed sequentially or simultaneously and, if sequentially, the order in which they are performed is not material. The method 200 then applies (at 230) the accessed real-time data set and the estimate to a probabilistic estimator to yield a probability of an occurrence of a kick and a confidence measure for the probability.

[0061] Once the probability and its confidence measure arc obtained (at 230), it may be used in a variety of ways. In one embodiment, it is communicated to a drilling engineer or some other operator who then decides whether corrective action is warranted and. if so, what that action might be. Or, the process may be automated so that when the probability breaches a specified threshold within a specified confidence measure, certain corrective actions are automatically taken. What these corrective actions might be will be implementation specific and will depend on the circumstances of (he kick within the context of the well. The probability and its confidence measure may also be archived for review at a later date.

[0062] Turning now to Figure 3 and Figure 4, the process flow 400 encompasses the computer- implemented method 200 of Figure 2. In operation, the sensors 136 and 137 sense their respective quantities and communicate those values as described above. The well system model 324 is previously constructed using a priori knowledge 403 of the well, such as the well geometry, the formation structure, etc. and is a physics-based, state space model of the well system. (Examples of two suitable models are given below.) Inputs to the well system model 324 can be defined as prescribed boundary conditions of the model. For example, these can be pressures, flow rates, temperatures, geometry, and mole fractions. These arc generally values that one can set in the operation of the well and can be static (i.e.. constant) or dynamic (time-varying).

[0063] Both the well system model 324 and the real-time information 410 will have uncertainties associated with them. More particularly, the welt system model 324 includes model parameter uncertainties 42S and the real-time information 410 includes measurement uncertainties 420. Model parameter uncertainties 42S will typically arise from variability in mud and formation properties. Measurement uncertainties 420 will typically arise from margins for error in the sensors used to take the measurements.

[0064] The data 315 comprises measurements of conditions in the wellbore 115 of the well 103 and at the surface as described above. The real-time information 410 is selected from the data 3 IS because it is correlative to a kick. Thus, the identity of the real-time information 410 will depend not only on what data 315 is available, but on its relationship to the presence or absence of kick. The real-time information 410 is "real-time" in the sense that it is input to the well system model as soon as it is available. Different sensors will sample at different rates, and thus some of the real-time information 410 will be fresher than will some other information. But the real-time information 410 constitutes the freshest information available at the time given the rates at which the data 315 is sampled.

[0065] The physics-based, slate space well system model 324 generates an estimate of a modeled condition correlative to a kick. A kick can generally be represented by a downhole or surface condition that is quantifiable but not amenable to direct measurement. For example, a kick may be indicated by an influx of formation fluids that cannot be directly measured, but that will affect the values of quantities that can be measured, such as those discussed below. The well system model 324 estimates a value for just such a quantifiable, not directly measured, condition.

[0066] In the illustrated embodiments. Ac real-time simulation 415 also yields an uncertainty measure, which is a measure in the confidence of the estimated value. The uncertainty measure is a function of the model parameter uncertainties 42S. This information will be known from the implementation of both the well 103 and the well system model 324 and the formulation of the model. For example, certain assumptions may underlie the design of the model and introduce uncertainties into the results. One such set of assumptions is discussed further below in connection with a particular model.

[0067] The estimate from the real-time simulation 415 obtained from the well system model 324 and its model parameter uncertainties 425 are then applied along with the real-time information 410 and its measurement uncertainty 420 to a probabilistic estimator 327. The probabilistic estimator 327 then yields a probability of an occurrence of a kick and a confidence measure for the probability 352. In the illustrated embodiments, the uncertainties arc represented by Gaussian distributions but other types of distributions may be used as well. Furthermore, the probabilistic estimator 327 is a Baycsian estimator although alternative embodiments may employ different probability theories.

[0068] The probability and the confidence measure 352 are then communicated to a drilling engineer in this particular embodiment. The manner in which the communication is performed and to whom the communications is made will be implementation specific. For example, the probability and the confidence measure 352 may be communicated by rendering it into a graphic display in human perceptible form for viewing by an operator of the well. Alternatively, the probability and the confidence measure 352 may be communicated to an alarm that automatically sounds if the value of the probability and the confidence measure 352 exceeds some predetermined threshold.

[0069] A more detailed disclosure of one particular implementation of the techniques described herein shall now be provided to further an understanding of the subject matter claimed below. The technique detects a "kick", which as described above is an unwanted penetration of fluids from the formation into the wellborn. The embodiment now being described is concerned with kicks arising from the influx of gas from the formation. When the gas enters the wellbore, it can rise up the annulus either as free gas or dissolved gas in drilling mud. As it encounters lower pressure regions at the top of the annulus. it expands, and dissolved gas comes out of solution. In detecting kicks early, well control personnel can isolate the influx and circulate it out while re-balancing the well for continued operation

[0070] In the context of Figure 1, a kick may occur when the formation fluids 121 penetrate into the wellbore 115. Such a condition may be caused in a number of ways. For example, the volume or density of the drilling fluids 118 might drop so that the hydrostatic pressure exerted by the drilling fluids 118 is less than the pressure to which the formation fluids 121 are subject Or, motion of the drill string 112 in the wellbore 115 might cause the hydrostatic pressure to effectively decrease, thereby creating a pressure differential leading to a kick. Those in the art may appreciate other ways in which such a pressure differential may be created and, thus, other ways in which a lack may be initiated.

[0071] There arc known downhole conditions that may be considered 'kick indicators". For present purposes, indicators can be cither primary or secondary. Primary indicators are those changes that are attributable to kicks alone, while secondary indicators may be caused by other drilling anomalies or well maneuvers. Primary kick indicators may include an increase in outflow rate, mud pit gain, incorrect fluid fill while (ripping, positive flow while pumps arc off, etc. Secondary kick indicators may include a decrease in stand pipe pressure and pump pressure, an increase in gas content in outflow mud, increase in rate of penetration, etc. Still other indicators may be known to those in the art having the benefit of this disclosure.

[0072] These quantities are considered "indicators" because they are correlated to a kick. For example, an increase in outflow rate may be an indicator because sustained deviation between known inflow rate and measured outflow rate could be caused by a kick For a mud pit gain, the closed mud loop serves to circulate mud around the well with the mud pit serving as a storage tank. An increase in the volume of fluid in the mud pit could be an indication of influx from the reservoir. On the other hand, for incorrect fluid fill while tripping, if pit volume does not reduce by an amount equal to the volume of steel being removed while tripping out, a kick may be occurring.

[0073] Similar correlations can be made for secondary kick considerations. Λ decrease in stand pipe pressure and pump pressure can be caused by gas influx into the armulus, which causes a decrease in the density of annulus fluid, and consequently a decrease in the hydrostatic pressure that creates a pressure deferential between drill pipe and annulus. This forces fluid from drill pipe to annulus, effectively reducing standpipe pressure. Note that the pump pressure should increase initially through exposure to the influx fluid, then decrease with continuous influx. For an increase in gas content in outflow mud, the percentage of gas in mud increases with kick, although this may also mean that a gas-bearing formation has been drilled through. An increase in rate of penetration (drilling break) occurs when more porous rock formations are encountered, which comes with increased risk of gas kick.

[0074] The well system model 324 in this embodiment incorporates a model of a kick in the context of the well system 103. Inherent in a model-based approach is the assumption that all computational parameters and variables, whether surface or downholc, can be transferred in real-time to calculation servers and that results from the computer models arc immediately available for application. A first, detailed model approach for the well flow system and formation in a discretized distributed flow model will now be discussed An alternative will be discussed afterward. [0075] The first model is expressed mathematically in a series of equations using a number of variables. As those in the art will appreciate, mathematical expressions are simply stand-ins for verbal descriptions. For example, one might verbally refer to "gravity" while using the symbol "g" to represent it mathematically. Both expressions represent the same thing. The variables and the quantities they represent used in the equations below likewise represent physical, real world quantities in the downholc environment, both measured and calculated. They therefore are not abstractions and the equations representing them are not abstractions, but rather descriptions of tangible, physical objects and conditions. Each variable will be defined as it appear* in the course of the discussion. However, for convenience, they are also collected in Tabic 3 toward the end of this detailed description.

[0076] The equations discussed below model the transient hydraulics and well-formation interactions in single and multiphase flow. The drill string and annuhis will be spatially discretized and balance relations and closure equations are defined for each discrete space. The physical effects estimated in the model are the frictional pressure loss, both for single and two phase flows; pressure loss in bit; viscosity variations with pressure, temperature and composition of the mud; density variations with pressure, temperature and gas content of the mud; dynamics of gas dissolution in mud (non -equilibrium): rise in gas velocity as it expands up the annul us; and simple reservoir dynamics including permeability and porosity of reservoir (when a reservoir model is included).

[0077] The model assumes that all variables are dependent on only one spatial coordinate - length along flow line. Effects from cross-sectional, non -uniform velocity and mass distribution profiles arc neglected It is also assumed that temperature at each point along the flow line is known. ("ITiis is an input to the model based on estimates or measurements made elsewhere.) Additional assumptions include that gas in the flow line can exist either as free gas or dissolved gas; gas and mud pressures at the same point are assumed to he equal; and gas is insoluble in water-based mud, hence single phase flow. The system is treated as a black oil system, one that is able to predict compressibility and mass transfer effects between phases in a reservoir as it is depleted.

[0078] The conservation of mass and momentum for a compressible fluid form the fundamental governing equations for this engineering problem. Though flow is assumed to be onc-dimcnsional inside a cylindrical pipe or annular region, cross-sectional geometric flexibility is accommodated through a variable diameter formulation. This discussion begins first by presenting the relevant equations for a single-phase flow and then expands on this formulation for multi-phase flow.

[0079] In one dimension, the governing equations for a single-phase (mud) flow in conservative form arc as follows: and

where is the mud density, A is the local hydraulic diameter, t is time, is the mud flow velocity in the x direction, x is a spatial coordinate, is a frictions 1 force term discussed further below and included to model viscous effects, P is the pressure of the fluid, and g is a gravity acceleration.

[0080] For single-phase flow, the density of mud, is derived from correlations for slightly

compressible fluids as follows:

where p c are the density, temperature and pressure of mud at standard conditions, respectively,

and T is temperature, P is pressure, and E is a volume modulus The parameter c, is the mud compressibility constant . These parameters are considered model input such that pm can be calculated directly given the local pressure and temperature, (P, T).

[0081] Additional pressure losses incurred from fluid friction arc calculated using empirical correlations. For single phase flow, the friction loss term, fa ^' given by:

where d _h is the local hydraulic diameter and/ is a friction factor that is determined separately depending

on whether the local flow is laminar or turbulent For laminar flow (Re < Ret, where Re is the Reynolds number and is the highest Reynolds number limit for laminar flow).

and for turbulent flow (Re >Re _T , where Re is the Reynolds number and Re _T is the lower Reynolds number for turbulent flow), where

where is die yield point and μ _ρ is the fluid plastic viscosity of the fluid in the wellbore for a Bingham plastic model of the fluid.

[0063] The local Reynolds number (Re) and associated transition points arc computed as follows:

where

and

for the drill string, and

for the annuhis.

[0083] To model flows which may include some distribution of kick hydrocarbon in the early flow dynamics using oil-based mud, a multi-phase flow solver is desired. The model used here is based on tracking three constituents: the free gas in the system, the gas dissolved in the drilling mud, and the drilling mud itself. The drilling mud is made up of water. oiL weighting solids, and dissolved gas. The governing principles are conservation of mass and conservation of momentum. Three conservation of mass equations arc used: one each for the free gas, dissolved gas, and drilling mud. Conservation of momentum is expressed via a single partial differential equation governing the momentum of the entire mixture.

[0084] The governing equations take the following form. The various quantities used in Eq. (16) through Eq. ( 19) arc collected in Table 1 for convenience. For mass conservation for the mud, the equation is:

And, for the momentum conservation for the entire mixture:

[0085] The models use a variety of variables that can be categorized as follows. The model employs two independent variables, time f (sec), and position x (ft). There are four "state variables": pressure p mud velocity volume fraction of fixe gas a, and mass fraction of

dissolved gas in mud ψ. There are also six derived quantities: mud density rate of free gas

, velocity of the free gas free gas

injection term q (Ibm/ft secX and force due to frictional effects To close the system, the models define all the "derived quantities" in terms of the state variables (or other derived quantities that can be computed explicitly from the state variables) as well as the "given" quantities cross-sectional area A (ft ² ), temperature 7'(°R), acceleration due to gravity g (ft/sec ² ), and the wellbore angle from the vertical Θ.

[1086] These models are generally referred to as "submodels" or "closure models" and there are six of them: the mud density pm-, the free gas dissolution rate nig, the free gas density pg, the free gas velocity ug, the injection source q, and the frictional force Ff. There is a lot of literature concerning various options for these submodels. The choices used in this work will be discussed in detail below.

[0087] To write the partial differential equations above in a compact form, l

denote the primitive variables and

denote the conserved variables. Given appropriate submodels, the conserved variables can be computed from the primitive variables and the given quantities. Letting the given quantities be denoted by W= [A. T \νκ have

(21)

[0088] Further, denote the fluxes and sources, respectively, in Eq. (16)

through Eq. (19). Specifically,

Then, Eq. ( 16) and Eq. ( 19) can be written in the following compact form:

Procedures for disardi/ing liq. (2()) are described below.

Table 1: Variables appearing in Eq. (16) through Eq. (19).

0089] The multi-phase governing equations given above benefit from closure relationships for a number of quantities. These quantities include: the density of the free gas, the density of the mud, pm \

the velocity of the free gas, ug; the friction or viscous force, F, the rate of gas dissolution, rog; and the gas influx rate, q. As promised above, the model* used in this work will now be discussed in detail.

[0091] Except for r, Eq. (24) is the ideal gas law. Thus, the compressibility factor is a dimensionless number that accounts for the departure of the gas from ideal gas behavior. It is computed from the Ha 11- Yarborough correlation:

where /v is the pseudo reduced pressure ratio, is the reciprocal pseudo reduced temperature ratio, and>>

is the reduced density. The pseudo reduced pressure ration and inverse pseudo reduced temperature ratio are given by

[0092] Correlations are then used to compute the critical pressure and temperature:

Finally, y is given by solving the following nonlinear equation:

[0093] Turning now to mud density, for clarity and simplicity, wc first consider the case whac there is no gas dissolved in the oil. Wc then generalize to include the effect of the dissolved gas.

[0094] When there is no gas dissolved in the oil that makes up part of the mud, the mud density is given by

where are the mass fractions of the water, oil, and solids (weighting materials) within the mud and arc the respective densities. The weighting materials arc incompressible and thus

is a constant, but the densities of the water and oil depend on p and T. For example, the following correlations have been proposed:

where

Note that these correlations use for density, (°F) for temperature, and (psia) for pressure. Thus, appropriate unit conversions are performed to use these results.

[009SJ In the case where gas is dissolved in the oil, Eq. (34) through Kq. (36) are modified to accordingly. It is assumed here that the gas is insoluble in the water component of the mud. Thus , all the dissolved gas in the mud is dissolved into the oil. In this situation, recalling that φ denotes the mass fraction of the dissolved gas in the mud, we have

[0096] It is common to express the density of the oil with dissolved gas in terms of two ratios: the

correlations for the gas/oil ratio in terms of p and Γ, these are unnecessary here because the mass fraction of dissolved gas is known.

[0098] To complete the model, a correlation for j is desired. Many such correlations are available in the literature. In the illustrated embodiment, the correlation for B ₀ depends on the bubble point pressure /¾ and the formation volume factor at the bubble point pressure, Specifically.

(43)

[0099] Correlations for , and C arc given below. Note that the correlation used for C depends

on whether the local pressure is greater or less than the bubble point pressure. The bubble point pressure is computed from

where the parameters are

and

Finally, note that T in Eq. (44) is in Ψ.

[00100] The formation volume factor at bubble point pressure is computed from

where

[00101] The free gas moves relative to the drilling mud. Thus, a free gas velocity model is used to close both the free gas mass conservation and the momentum equations. The model used here expresses the free gas velocity as

(49)

where is the free gas velocity, is the gas/mud mixture velocity, the slip velocity. The mixture

velocity is given by

(50)

Substituting Eq. (SO) into Eq. (49) and solving for gives

(51)

[00102] To complete the model, a correlation for the slip velocity is introduced. Here, a very simple slip velocity is used. Experiments by others in the art indicate that the gas rise velocity in oil-based mud is essentially independent of the volume fraction of free gas, perhaps because the large slug-type bubbles form at very low free gas volume fractions. Based on these experiments, \vc adopt the following simple model:

where d is the outer diameter of the annulus and r is the ratio of the inner diameter to the outer diameter

[00103] This model is, for vertical wells and can represent deviated wells through an angle correction. The angular correction has not been implemented here but those in the art having the benefit of this disclosure will be able to add it if it is found necessary or desirable. Similarly, there are other models known to the art that may be suitable, These may be ueed in alternative embodiments. Indeed, any suitable model known to the art may be used.

[00104] Turning now to the friction factor, the force due to friction on the right-hand side of the momentum equation is modeled as

where is the hydraulic diameter and /is the friction factor. The friction factor is determined as

(54)

where is the "no-slip" friction factor. For the no-slip friction factor.

(55) where

(56)

and are the viscosities of the mud and free gas, respectively. Note that this friction factor is just

that given by the "smooth wall" curve on the Moody diagram. Further, this friction factor is based on data for pipe flow of Newtonian fluids. Thus, some embodiments may choose to use a correction to account for the non-Newtonian nature of the drilling mud.

[00105] The quantity 5 in the correction that accounts for multi-phase flow is given by

[00106] Equation (63) has a singularity near y 1 (slightly greater than 1), and so, when 1<y<1.2, one can replace Eq. (63) with

At y-1, this switch is continuous. At y=1.2 it is not (but the discontinuity appears fairly small). The derivatives with respect to y are not continuous on either side.

[00107] The rate at which free gas dissolves into the mud is dependent on many factors, including the solubility of the gas in oil (as measured, e.g., by the gas/oil ratio at saturation), the "distance" from the saturated state (as measured by the difference between the actual gas/oil ratio and the saturation gas/oil ratio), and many other factors. Unfortunately, the literature on gas kick simulation docs not fully specify an appropriate model for this effect. The illustrated embodiments employ a non-equilibrium model primarily based on dimensional analysis and some assumptions. This will allow us to begin simulations and investigate the sensitivity of the results to features of this model. [00108] As a starting point for the model, we consider the following model form for dissolution:

( ⁶⁵⁾

where C is the concentration (moles of solute per unit volume of solution), C is the concentration at saturation (i.e., the solubility), and k is a rate constant Multiplying Eq. (65) by the molar mass of solute M (gas in our case) gives

Since

We have

(68) The mass fraction of dissolved gas at saturation <p, can be computed from a correlation for the gas oil ratio at saturation and hence will depend on p and 7'.

[0 0109) When the mud is entirely made up of oil (such that φ is (he mass fraction of dissolved gas in the oil), we have the following:

Thus,

(70)

where the gas 'Oil ratio at saturation can be computed via a correlation:

For hydrocarbon gas in base oil

where ft is the specific gravity of the gas (e.g., -0.6409 for natural gas) and the temperature Γ is given in F. Note that this correlation gives in scf/bbl. To convert this divide the result by 5.61458.

[00110] To complete the model, specify k, the rate constant The illustrated embodiments uses the following definition:

where is a constant.

[00111] The model for gas influx is specified at known location and time. Recall from Eq. ( 18):

Gas influx rate, q, is specified using a simple linear model:

where is a reservoir constant specified to give a desired flow rate for pressure, P{y), at varying reservoir depths, y. Alternative embodiments may employ alternative models.

[00112] The partial differential equations documented here can be discretized using a large variety of different methods. This section describes some of those methods will now be discussed. For the purposes of compactly describing the different methods, the notation introduced in Eq. (23) is used throughout this discussion.

[00113] denote a partition of the domain Ω, and let

Further, let

where and U, is the cell-average of u on the Ah cell. Integrating Eq. (23) over Ωι ^' gives

where F, <s and F, arc the flux at.s »; and x. respectively, and £ is the ccll-avcragc source term. Tlicsc fluxes and sources cannot be computed exactly given only the ccll-avcragcd quantities. Instead, generic ccll- centered finite volume methods arc derived by developing approximations for these terms.

[00114] The Lax-Friedrichs method is written as follows:

where is the f lux evaluated using the state in cell j at time n . The method can be shown to be first- order in both space and time and is monotone. Further, it is very easy to implement and very robust. Thus, it represents a good scheme to start with, allowing development and testing of the physical models described earlier. However, it is well-known to be very diffusive, even compared to other first -order methods.

[00115] The Roc scheme can be written as follows: where is a Roc-flux function. The Roc flux can be mitten as follows:

where

and is any matrix such that

is diagonatizablc with real eigenvalues; and

[00116] For some systems - e.g., the Euler equations - it is straightforward to analytically construct such a matrix A . However, for the general multi-phase flow equations with complex submodels used here. this is not at all trivial. Instead, we use a '"numerical" Roc matrix. The matrix A is computed by the following procedure.

[00117] First, compute an average conserved state vector U _m by computing the conserved state corresponding to the average of the primitive states on the left and right of the interface. For this is given by

[00118] Next, evaluate the flux Jacobian at this average state:

(83)

[00119] Then modify' to satisfy the first criterion set forth above above. Specifically, find a diagonal matrix such that

where R are the eigenvectors of A _m . Letting Λ be the diagonal matrix of eigenvalues of /¾_, . we find Λ by letting

[00120] The above may be referred to as a "Distributed Hydraulics ModcP, or "DHM" and may be employed in some embodiments. However, alternative embodiments may use other types of models such as the "Lumped Parameter Model". The Lumped Parameter Model, or "LPM", provides a real-time tool for monitoring well processes as well as detection of reservoir influx at the bottom hole. It models well hydraulics and combines it with well measurements in an optimal way that accounts for uncertainties in each as shown in Figure 5 and Figure 6. It also incorporates a Confidence Interval on the Expected Value which establishes a bound on the estimated variables including any influx. This serves to help eliminate false positives. The LPM selectively combines several subsidiary techniques including flow measurement and well monitoring systems, flow models for predictive systems, and probabilistic models.

[00121] Flow measurement and well monitoring systems include flow meters, mud pit volume sensors and stand pipe pressure gages. Typically, a kick threshold for any or all of these parameters is set and the system generates an alarm if the set maximum is exceeded. Many different types of flow meters are in use today. In practice, the kick threshold for outflow rate is set at a specific value of outflow minus inflow, known as delta flow. This precludes the need for continual resetting of alarm levels when drilling conditions demand a change in the inflow rate.

[00122] Flow models for predictive systems include process models, which have found increasing use in kick prediction with the availability of high speed computers. Real-time, advanced mathematical models incorporating multi-phase flow, torque and drag models as well as several sub-models, compute flow out and other well parameters as the drilling process progresses using inputs from installed sensors along the flow line. This is men compared to real-time well data and any discrepancy is used as a predictor of kick or other drilling anomalies.

[00123] Probablistic models use a model matching framework based on Bayesian probability. Kicks of different types and rates arc modeled and compared to real-time data using Baycs rule. Other rig activities are also modeled to reduce incidences of false positives. The system outputs the kick probability at each data point and when it exceeds a set threshold (90%), an alarm is raised. It uses flow out/flow in comparison as the primary kick indicator. It is claimed to have high, adaptable sensitivity with low false alarm rate. It is also rig independent, requires little or no calibration and can use crude flow meters like the paddle meter.

[00124] For the LPM model, rather than the complex multi-phase flow models described above, which involve solving partial differential equations of mass, momentum and energy conservation, wellbore hydraulics is simplified into time-only dependent mathematical models with the wellbore lumped as a single block as illustrated in Figure 7. Process and measurement equations are obtained from a bond graph model such as that in Figure 8 of wellbore and well reservoir hydraulics. The equations arc then linearized and transition matrices obtained. These matrices form a basic component of the Linearized Kalman filter used.

[00125] Several assumptions arc made in simplifying the wellbore hydraulics. These include, for example, influx enters the wellbore at the same density as the drilling mud and remains at this density for the early stage of kick detection. Hence, only a single liquid phase is considered. The fluid is incompressible. This proceeds from the assumption of a single liquid phase. Nonlinear, square law pressure drop assumed for drill pipe and bottom hole assembly, RA* and aimulus, Reservoir pressure is modeled as a non-zero mean random walk where the bias and diffusion strength are known from experimental data. Mud inflow rate is known hence there is no need to include the drill string fluid momentum subsystem.

[00126] The result is two state functions. The first one, the fluid momentum, /'. describes the wellbore- reservoir hydraulics, and the second, the mud pit volume change, , as a result of well influx.

[00127] The proposed model is uncertain due to the simplifications assumed in the construction of the bond graph and the inherent measurement uncertainties in the data supplied from the wells. The dynamic system is augmented to include formation pressure, as a shaping filter for the random walk process.

where formation pressure, is modeled as a random walk Gaussian process with zero mean and variance = , resulting in the evolution equation of the form with / as a

deterministic mapping of the state vector and _* as the additive noise associated with the

process.

[001281 Three measurements are used for the process estimation as a basis for comparison with the state vector. These are the pump pressure, , return flow, and the mud pit volume They make up the observation vector, which has the form

The measurement noise vector is also modeled as an additive Gaussian process noise with zero mean and variance given by

[00129] The process estimation process consists of the calculation of the probability distribution of that is, the states, given all available measurements and the nonlinear models.

which has an initial, in the form of a random vector of mean and

A random vector, captures the uncertainlies in the model while

another random vector, captures the noise in the measurements. Both of them are described by:

One solution for linear Gaussian models is the Kalman filter. For nonlinear and/or non-Gaussian models, sequential Monte Carlo methods are used to construct approximate solutions.

[00130] The Kalman filter is based on a linear Gaussian model. For nonlinear, Gaussian systems, the Linearized Kalman filter and the Extended Kalman Alter may be used to approximate the solution. These methods are based on linearization of the state and measurement functions about a steady state value, resulting in the following state and measurement matrices:

For a stcady-etatc linearization about the inflow rate, we get a constant matrix for the Linearized Kalman case, given by

The measurement matrices become:

which represent a continuous system of the form

where A are the 3x 3 continuous matrices above. These are converted to discrete time system

using Zero Order Hold transformation to obtain

where are discrete matrices. This linearized discrete system is used in the Linearized Kalman

Filter in Tabk 2.

Table 2: Linearized Kalman Filter

[00131] The submodels collect such information as well geometry, formation characteristics, mud properties, and information on current drilling maneuvers to calculate parameters used in process estimation and to make decisions on whether changes in the kick indicators are attributable to influx or to current well operations. The sub-models arc described below:

[00132] The annular pressure drop is given by where the annular pressure loss coefficient,

is a constant obtained at the steady state inflow rate which is a known input into the system. The

rheological model used to develop the friction pressure loss sub-model is the non-Newtonian, Power Law model. A preliminary annular pressure loss is calculated in field units as

where the friction factor, f . depends on whether the flow is laminar, turbulent or in transition as determined by the value of the dimcnsionlcss Reynold's number, Re/ is found for laminar and turbulent flows as

For transition flow, /is interpolated between die two values above and given as

where,

The annular pressure loss coefficient ίκ then calculated as

[00133] For practical use in drilling operations, the model has to accommodate changing wellbore geometry for each bit run. Wellbore length or depth is calculated at each new time step by monitoring the rate of penetration ("ROP"), such that

Alongside the depth, wellbore area is also continuously monitored at each time step. Sections of uniform area have the same fluid inertia given by

The different sections with differen t areas are aggregated to get the total fluid inertia:

[00134] The rate of penetration is determined using the following model:

Where D is the true vertical depth, a, to a, are constant coefficients to be determined and to J¾ are drilling parameters. Eq. ( 126) can be written as

The function f _t models the effect of parameters such as formation strength, mud type, bit typo and solid content. This is given by, The functions model the effect of compaction thusly.

The functions model the effects of overbalance, weight on bit (WOB). and rotary speed respectively. Thus,

Lastly, the functions model bit tooth wear and bit hydraulics

lite estimator adopts a simplified form of Eq. (127) based on Eq. (131). the overbalance function. This is shown in Hq. ( 136) below:

Where the effects represented by functions barring have been concentrated in a nominal ROP,

is a nominal pressure variation function, and is one of the state variables

obtained from Eq. (106) at every time step.

[00135] In general, the LPM is advantageous relative to the DHM in that it uses existing rig process measurement data and continually updates this at every new data point as drilling progresses. No additional measurement parameter or equipment is needed. The system works within the uncertainties of sensors in current use, including the inaccurate flapper used for flow measurements. Set uncertainties for important variables increase noise tolerance and help keep fake alarm rates at a minimum, if not totally dim mated. Rig and process specific data collection is minimal. It works on a broad range of rigs, from land rigs to deepwater well drilling. It uses mud pit volume increase as the primary kick indicator.

[00136] The volume of influx that trips the alarm can be set to any level acceptable to the drilling crew thereby accommodating differences in rig types and peculiarities. Even for deepwater wells, the procedure ensures that there is no time delay between an occurrence at the bottomholc and observation at the wellhead. Kicks or losses bottomholc cause immediate changes in the pump pressure which is used as the primary driver of the prediction process. Hence it ends up being a faster means of kick detection than outflow rate. The volume of influx taken in is known in real-time, with a confidence interval on the accuracy of results. Advantages of using pressure as the primary driver are harvested, ^' fhese include: sensors do not fail due to gas flow; high accuracy of measurements; can predict flow rate as well; are a normal part of the rig Bystem; fast reaction time to downhole changes.

[00137] On the other hand, the assumption of incompressible flow in the wellboare annul us may lead to over predicting the rate of influx into the well bore for slower kicks when some gas phase may be present. Increased friction pressure loss associated with mis assumption may dampen this effect. Incompressible flow assumptions also give rise to immediate topside response to well bore influx, which may not be realistic when well breathing effects (elasticity in the mud/formation interaction arc significant, or when gas phase material is present), or when significant topside mud fill and drainage occurs (within piping between the outflow meter and the mud pits).

[00138] The current LPM includes a model of the resistance to flow between the well bore and the formation which is linearized and therefore independent of the direction of flow. A non-linear resistance, which is dependent on flow direction can be added to the LPM Estimation of the resulting non-linear model can be obtained by non-linear estimation methods such as statistical linearization and Unscented Kalman Fitter methods. Mud is intended to providing sealing effect with the formation and increase the resistance to outflow or mud loss, which is non-linear. The LPM does not resolve effects along the length of the annular region. It therefore is insensitive to where in the open hole an influx may occur, and assumes that it occurs at the bottom hole region.

[00139] Turning now to Figure $ and Figure 6, as mentioned above, this particular embodiment includes an update/correction feature. Figure 5-Figurc 6 convey how combining multiple models/predictions of the same quantity gives significantly reduced uncertainty in the estimated value. More particularly, this embodiment employs a technique by which even noisy or poor estimates and measurement can be combined arrive at predictions that arc less noisy and better than either of the those that were combined. In this context, "noise" is "uncertainty" in either the estimates or the measurements as discussed above.

[00140] Figure 5 includes three curves 500. 503, 306, each representing an uncertainty distribution. The distributions are Gaussian but for illustrative purposes only as any kind of distribution that is suitable to the data may be used. The curve 500 represents the uncertainty distribution for a first measurement and the curve 503 represents the uncertainty distribution for a first estimate. The curve 506 represents the combined measurement and estimation uncertainty distribution. Notice how reduced the uncertainty in the combination is despite relatively large uncertainties in both the measurement curve 500 and the estimate curve S03. Figure 6 illustrates how the principle can be extended through a second iteration. Thus, embodiments employing this technique for updating estimates can combine a first estimate with a first uncertainty and a measurement with a second uncertainty to obtain a second estimate with a third uncertainty, the third uncertainty being lees than the first uncertainty and the second uncertainty.

[00141] The presently disclosed technique docs not just trigger on a pattern in the data but provides a quantifiable estimate of a kick with quantifiable uncertainty. Since it is based on physics prediction as compared to empirical models and methods, it should be more adaptable to new configurations and changing environments. It combines multiple measurements of drilling operations by linking the measurements with the physics of the operation. This provides for natural scaling of the measurements relative to each to other for making predictions of output variables. It also provides for natural filtering or smoothing of the estimate, sometimes called "physical filtering", instead of ad hoc smoothing or averaging of the measured data as found in conventional practice. Note that not all these characteristics will necessarily be found in all embodiments and, where found together, may not all be manifested to the same extent

[00142] The efficacy of the presently disclosed technique is illustrated in Figure 9. The trace 900 represents die performance of the presently disclosed technique. The trace 905 represents the performance of a conventional measured mudpit technique. Note that the kick is detected at time 910 for the disclosed technique (i.e., when the trace 900 crosses the alarm threshold 915) sooner than does the conventional technique, which detects the kick at time 920 (i.e., when the trace 905 crosses the alarm threshold 915). This earlier detection of the kick will typically be advantageous in responding to its occurrence.

[00143] In the embodiments set forth above, the sensors 136, 137 and the computing apparatus 145 (including well monitoring software component 32 land data 318) comprise a well monitoring system. The technique can also be integrated into well management and monitoring techniques such as arc known to the art, primarily by retrofitting the software architecture with the functionality of the well monitoring software component 321 described above. The embodiments disclosed above are presented in isolation from other wells and/or operations that might be happening nearby.

[00144] For example, wells arc typically drilled in a field containing other wells. Well management and monitoring techniques are sometimes implemented across multiple weds, for example a number of wells within a field. Thus, well monitoring and management techniques such as those disclosed in U.S. Application Serial No. 14/196307, U.S. Application Serial No. 13/312,646, and U.S. Letters Patent 8,121,971, may be modified to implement the techniques disclosed herein. The manner in which such techniques known to the art may be modified to implement this technique will be readily apparent to those skilled in the art having the benefit of this disclosure.

Table 3. Summation of Values Employed Above

00145) The following patents referenced above are identified more completely:

[00146] U.S. Application Serial No. 14/196,307. entitled, "System and Console for Monitoring and

Managing Well Site Operations," filed March 4, 2014, in the name of the inventors Fereidoun Abbaesian el al., and published September 4, 2014, as U.S. Patent Publication 2014/0246238.

[00147] U.S. Application Serial No. 13/312.646, entitled, "Geological Monitoring Console," filed

December 6, 2011, in the name of the inventor Paul J. Johnston and published June 6, 2013, as U.S. Patent

Publication 2013 0144531.

[00148] U.S. Letters Patent 8.121,971. entitled, 'Intelligent Drilling Advisor", and issued February 21, 2012, to DP Corporation North America Inc., as assignee of the inventors Michael L. Edwards et al.

[00149] This concludes the detailed description. The particular embodiments disclosed above are illustrative only, as the invention may be modified and practiced in different but equivalent manners apparent to those skilled in the art having the benefit of the teachings herein. Furthermore, no limitations arc intended to the details of construction or design herein shown, other than as described in the claims below. It is therefore evident that the particular embodiments disclosed above may be altered or modified and all such variations are considered within the scope and spirit of the invention. Accordingly, die protection sought herein is as set forth in the claims below.