Abstract: Method and system for determining a characteristic of a subsurface formation traversed by a wellbore. A fluid flow rate in the wellbore is established, the flow being associated with fluid movement in one or more layers in the formation. Fluid data is generated from the fluid flow. The fluid flow rate in the wellbore is selectively altered, and fluid data is re-generated pertaining to the altered flow rates. The fluid data is processed to determine desired reservoir parameters and/or characteristics of an identified stratum in the formation. A method identifies inflexion points on a single phase oil only curve plot using the fluid data to determine layer characteristics and flow stability.

Abstract: Disclosed herein are methods, systems, and devices for determining parameters of an earth formation. Pressure transient data from a formation test can be recorded and analyzed using an analytical model including one or more correction factors derived from an assumption that an induced flow within the formation is hemispherical. Regression analysis of the refined analytical model and the pressure transient data results in accurate earth formation parameters.

Abstract: Disclosed herein are methods, systems, and devices for determining parameters of an earth formation. Pressure transient data from a formation test can be recorded and analyzed using an analytical model including one or more correction factors derived from an assumption that an induced flow within the formation is hemispherical. Regression analysis of the refined analytical model and the pressure transient data results in accurate earth formation parameters.

Abstract: A method for estimating a viscosity of a fluid in a rock formation, the method including: performing a first nuclear magnetic resonance (NMR) measurement with zero magnetic field gradient on at least a portion of a sample of the rock formation to obtain a first distribution of transverse relaxation time constants; estimating a first diffusive couple factor from the first distribution; replacing the fluid of the at least a portion of the sample with another fluid; performing a second NMR measurement with zero magnetic field gradient on the at least a portion of the sample containing the another fluid to obtain a second distribution of transverse relaxation time constants; estimating a second diffusive couple factor from the second distribution; and estimating the viscosity of the fluid using the first diffusive couple factor and the second diffusive couple factor.

Abstract: Methods for analyzing the connected quality of a hydrocarbon reservoir are disclosed. A model of a portion of the reservoir is divided into cells, each cell having a volume and some attributes, and wherein a speed function is assigned to a portion of the cells. A reference cell is chosen. A connectivity between cells in the reservoir is determined by solving an Eikonal equation that describes the travel time propagation, said propagating front progressing outward from a reference cell until an ending condition is met, said Eikonal equation being solved by a fast marching method with propagation velocity as a function of spatial position being provided by the speed function. Regions of the reservoir are characterized by their connective quality to the reference cell using the connectivity.

Abstract: A method for determining fluid saturation in a formation at a plurality of radial depths near a wellbore, the method including: obtaining multi-frequency nuclear magnetic resonance (NMR) response data for the formation; and processing the data to determine simultaneously the fluid saturation at each radial depth. A computer program product is provided.

Abstract: A technique is provided to determine a flow rate of a production fluid. The technique is utilized in a well having a gas lift system. Temperatures are measured along the well to create a temperature profile. The temperature profile is used to determine the flow rate of a produced fluid.

Abstract: A method is provided for estimating the amount of supercharging in a formation penetrated by a wellbore. According to the method, pressure fluctuation measurements are obtained at a position in the formation accessible from the well bore and at an adjacent position in the wellbore, and a model is provided which relates variations in pressure with time at these positions to one or more adjustable parameters from which the amount of supercharging can be estimated. The or each parameter is then adjusted to optimize the fit between the pressure variations predicted by the model and the measured pressure fluctuations, and the amount of supercharging is estimated from the adjusted parameter(s).

Abstract: The invention relates to methods and apparatus for determining a downhole parameter in an underbalanced drilling environment which include: selectively activating a first fluid flowing from the formation through a wellbore while under balanced drilled; detecting the activated first fluid, and determining a depth at which said fluid enters the wellbore.

Abstract: A method for determining true formation pressure in formations surrounding a fluid-containing borehole having a mudcake on the surface thereof, including the following steps: at a particular location in the borehole, monitoring the electrokinetic potential across the mudcake; modifying the borehole pressure at the particular location in the borehole; and determining the true formation pressure in the formations surrounding the particular location as being the borehole pressure at which the electrokinetic potential is substantially zero.

Abstract: A system and method is provided for determining flow profiles in a deviated well. The system and method utilize temperature measurements and a modeling technique that enable the use of temperature profiles in deriving flow profiles for fluid injected into deviated wells.

Abstract: A method, system and computer program product for determining permeability of an earth formation of a reservoir are disclosed. A method may include: providing formation testing data of the earth formation; defining layering of the earth formation; determining a geological model for each layer; modeling pressure transient in the earth formation by representing pressure transient within a layer using an analytical solution and representing pressure transient between two layers using a numerical solution, the analytical solution selected based on at least one of the respective geological model and data regarding a drilled well of the reservoir, and a parameter of the analytical solution and a parameter of the numerical solution being determined based on the formation testing data; and determining the permeability of the earth formation based on the pressure transient model.

Abstract: Method and system are disclosed for determining hydrocarbon CDVs in a multiphase environment. The method and system involve calculating a total hydrocarbons volume from an apparent oil volume and an apparent gas volume. The apparent oil volume is determined using C/O data from a nuclear logging tool operating in inelastic neutron scattering mode. The apparent gas volume is determined using neutron porosity data from a nuclear logging tool operating in pulsed neutron capture (“PNC”) mode. It is also possible to determine the total hydrocarbons volume using log data from formation thermal neutron capture cross section measurements, or a cased hole formation resistivity tool. The total hydrocarbons volume is then used to back-calculate the total hydrocarbons CDV at various formation depths using the logged C/O data. The resulting independently calculated total hydrocarbon CDV may be used for proper matching of log data and simulation results for history matching purposes.

Abstract: A fully-parallelized, highly-efficient compositional implicit hydrocarbon reservoir simulator is provided. The simulator is capable of solving giant reservoir models, of the type frequently encountered in the Middle East and elsewhere in the world, with fast turnaround time. The simulator may be implemented in a variety of computer platforms ranging from shared-memory and distributed-memory supercomputers to commercial and self-made clusters of personal computers. The performance capabilities enable analysis of reservoir models in full detail, using both fine geological characterization and detailed individual definition of the hydrocarbon components present in the reservoir fluids.

Abstract: A technique for utilizing tube waves to measure sound speed of fluids and other properties in the frequency range of about 5 to 100 kHz. A drill string is equipped with a sensor tube having a cavity filled with a downhole fluid such as borehole mud or formation fluid. An acoustic transmitter and an array of acoustic receivers are mounted on the tube in direct contact with the fluid. Processing circuitry calculates a property, e.g., sound speed, of the fluid based on time-of-flight of an acoustic signal generated by the transmitter and received by the array of receivers. Alternatively change in signal phase as a function of frequency may be employed by the processing circuitry. The technique is particularly suited to measuring the sound speed of borehole mud in situ.

Abstract: In a first aspect, a method for use in a time lapse, marine seismic survey includes accessing a set of baseline seismic data; accessing a set of acquired, time lapse seismic data; and interpolating a set of time lapse seismic data from the baseline seismic data and the acquired time lapse seismic data, at least one of the baseline seismic data and the acquired time lapse seismic data being multicomponent data. In other aspects, a program storage medium is encoded with instructions that, when executed by a computing device, perform the above method and a computing apparatus programmed to perform one or more of such methods.

Abstract: Disclosed is a supervisory control and data acquisition (SCADA) system using a spread spectrum or licensed frequency data radio network and communication method therefore allowing multiple slave hosts and slave devices or remove terminal unit (RTU) the ability to communicate data connectivity in a wireless environment.

Abstract: A method for increasing production in a reservoir, comprising performing a swabbing test at a depth in a pipe, wherein the pipe is located in a wellbore and wherein a portion of the wellbore is located inside the reservoir, periodically measuring, during the swabbing test, pressure in the bottom portion of the pipe using the pressure gauge to obtain a plurality of pressure measurements, wherein the pressure gauge is affixed to an inner wall of a bottom portion of the pipe, and determining a plurality of flow rates of fluid flowing from the reservoir through perforations in the wellbore into the pipe using a flow rate equation and the plurality of pressure measurements

Abstract: A method, practiced by a Hydraulic Fracturing Simulator software, is adapted for simulating a hydraulic fracture in an Earth formation, the formation including a multilayered reservoir, a mesh adapted to overlay the fracture in the formation thereby defining a plurality of elements, the mesh and the fracture defining one or more partially active elements, the method comprising the step of: setting up an Influence Coefficient Matrix, the matrix having a plurality of numbers, the plurality of numbers of the Influence Coefficient Matrix being calculated and determined such that the method, practiced by the Hydraulic Fracturing Simulator software for simulating the hydraulic fracture, models, or take into account an existence of, the multilayered reservoir and the partially active elements.

Abstract: An apparatus and method are provided for solving a non-linear S-shaped function F=ƒ(S) which is representative of a property S in a physical system, such saturation in a reservoir simulation. A Newton iteration (T) is performed on the function ƒ(S) at Sv to determine a next iterative value Sv+1. It is then determined whether Sv+1 is located on the opposite side of the inflection point Sc from Sv. If Sv+1 is located on the opposite side of the inflection point from Sv, then Sv+1 is set to Sl, a modified new estimate. The modified new estimate, Sl, is preferably set to either the inflection point, Sc, or to an average value between Sv and Sv+1, i.e., Sl=0.5(Sv+Sv+1). The above steps are repeated until Sv+1 is within the predetermined convergence criteria. Also, solution algorithms are described for two-phase and three-phase flow with gravity and capillary pressure.

Abstract: Method and system for determining a characteristic of a subsurface formation traversed by a wellbore. A fluid flow rate in the wellbore is established, the flow being associated with fluid movement in one or more layers in the formation. Fluid data is generated from the fluid flow. The fluid flow rate in the wellbore is selectively altered, and fluid data is re-generated pertaining to the altered flow rates. The fluid data is processed to determine desired reservoir parameters and/or characteristics of an identified stratum in the formation. A method identifies inflexion points on a single phase oil only curve plot using the fluid data to determine layer characteristics and flow stability.

Abstract: The present invention provides method for detecting the possible presence of surfaces layers of diagenetic origin by establishing divergence between well signatures such as logs, particularly for use in predominantly carbonate reservoirs.

Abstract: A method intended for gradual deformation of a Boolean model allowing best simulation of the spatial configuration, in a heterogeneous underground zone, of geologic objects defined by physical quantities. The model is optimized by means of an iterative optimization process from realizations including objects whose number is a random Poisson variable of determined mean, and by minimizing an objective function. In order to impose a continuity in the evolution of the objects in size, number, positions, within the model, a combined realization obtained by combining on the one hand an initial realization comprising a number of objects corresponding to a first mean value and at least another independent realization having another number of objects corresponding to a second mean value is constructed. An application is construction of a Boolean underground reservoir model allowing simulation of the configuration of heterogeneities such as fractures, channels, etc.

Abstract: A system is provided for measuring flow in a pipeline which, in use, carries a stratified flow of a first liquid and a second, denser liquid. The system comprises: a section of the pipeline in which the stratified flow is uphill; a section of the pipeline in which the stratified flow is downhill; and a holdup meter on each of said pipeline sections for measuring the holdup in each section. The flow rates of the first and second liquids can be calculated from the measured holdups.

Abstract: Systems and processes for using the same for developing analytical device methods are provided. Also provided are computer program products for executing the subject processes.

Abstract: A method for optimizing the development of a fractured hydrocarbon reservoir wherein the network permeability is determined using a reliable compromise between numerical and analytical methods which has application to oil reservoir development. The reservoir is discretized into a set of grid cells and a geometrical description of the fracture network in each cell is elaborated. A connectivity index is then deduced within each cell for the fractures. The permeability of the fracture network of the cells whose connectivity index is above a first threshold is determined and a zero permeability value is assigned in the other cells. Other thresholds can be determined so as to choose between a numerical method and an analytical method to determine the permeability. These permeabilities are exploited in a flow simulator so as to optimize the development of the reservoir.

Abstract: The invention relates to fluid sampling in a test that is used to determine physical and chemical characteristics of the fluids in a subterranean reservoir. The method reconstructs the entire pressure history of the fluid parcel that is captured in the fluid samplers during a test. Using this reconstructed pressure history of the samples, the quality of the samples, particularly, whether there is a phase change in the samples during the test, can be accurately quantified.

Abstract: A method of measuring a dry gas flow from hydrocarbon wells. The hydrocarbon wells produce a stream which contains liquids and natural gas. The method comprises providing an differential flow measurement device (orifice) located at a well site containing the hydrocarbon wells. The method further includes testing the stream for gas composition and physical characteristics. Additionally, the method comprises measuring a static temperature (T) of the stream, measuring a static pressure (Ps), and measuring a pressure differential (dP) in the differential flow measurement device (orifice). The method includes computing the liquid-gas ratio (LGR) of the stream, computing the wet gas flow rate (Qw) of the stream, computing a first dry gas flow rate (Qd1) and a second dry gas flow rate and utilizing a iterative process to correct and obtain the dry gas flow rate. The data may be transmitted to another location.

Abstract: A technique for fast performance prediction of multivariable model predictive controllers (MPC) for large-scale spatially-distributed dynamic systems is provided. When operating an MPC for a two-dimensional sheetmaking process, the controller's performance can be separated into spatial performance and dynamic performance. The steady-state spatial performance for different sheet property modes is first predicted, thereafter, the dynamic performance is predicted on the basis of the steady-state actuator profiles. During the steady-state spatial performance prediction, the original MPC's cost function is approximated to be a steady-state cost function. Then the steady-state measurement profiles and the actuator profiles are calculated through the steady-state cost function. The actuator profiles in time series can be efficiently solved by minimizing the difference between the actuator profiles and the steady-state actuator profiles in the face of constraints.

Abstract: A system and method using inverse-vector processing to iterate through a loop of three steps: set a guide direction, invert opposite vectors, and average vectors to update the guide direction, for smoothing seismic amplitude data. The inverse-vector method can overcome instabilities where the traditional structure-tensor approach fails. The inverse-vector smoothing is simple to implement and more computational efficient. The resultant dips and azimuths are spatially consistent and thus more interpretable and suitable for calculation of curvature and other dip based attributes.

Abstract: A “near wellbore modeling” software will, when executed by a processor of a computer, model a localized area of a reservoir field which surrounds and is located near a specific wellbore in the reservoir field by performing the following functions: (1) receive input data representative of a reservoir field containing a plurality of wellbores, (2) establish a boundary around one specific wellbore in the reservoir field which will be individually modeled and simulated, (3) impose an “fine scale” unstructured grid inside the boundary consisting of a plurality of tetrahedrally shaped grid cells and further impose a fine scale structured grid about the perforated sections of the specific wellbore, (4) determine a plurality of fluxes/pressure values at the boundary, the fluxes/pressure values representing characteristics of the reservoir field located outside the boundary, (5) establish one or more properties for each tetrahedral cell of the unstructured grid and each cylindrical grid cell of the structured grid, (6

Abstract: A system and method for managing a production from a wellbore that includes taking measurements relating to one or more selected parameters for each of the production zones over a time period and determining the occurrence of the cross flow in the wellbore using a trend of the one or more of the measurements. The system includes a processor that receives information relating to the measurements made over time relating to the selected parameters, wherein the processor determines or predicts the occurrence of the cross flow from a trend of at least one of the measurements.

Abstract: The invention relates to a method of estimating the flow restart conditions of a paraffinic hydrocarbon in a pipe, wherein the following stages are carried out: selecting a theological model, determining the temperature and shear stress fields in the steady hydrocarbon flow, evaluating the conductive thermal exchanges during circulation standstill phases, from the steady flow temperatures and the standstill time taken into account, selecting various sections of the pipe and recording the thermomechanical scenarios undergone by the hydrocarbon in these sections, experimentally measuring the rheological parameters of the model on a sample of the hydrocarbon, according to each one of said scenarios, and taking into account these experimental measurements to evaluate the value of the pressure that restarts circulation of the paraffinic hydrocarbon in the pipe.

Abstract: A method including: a dividing step (A) of dividing external data into a plurality of cells (13) having boundaries orthogonal to each other, the external data including boundary data of an object which contacts incompressible viscous fluid; a cell classifying step (B) of classifying the divided cells into an internal cell (13a) positioned inside or outside the object and a boundary cell (13b) including the boundary data; a cut point determining step (C) of determining cut points in ridges of the boundary cell on the basis of the boundary data; a boundary face determining step (D) of determining a polygon connecting the cut points to be cell internal data for the boundary face; and a analyzing step (E) of applying a cut cell finite volume method combined with a VOF method to a boundary of a flow field to analyze the flow field.

Abstract: Physical models of wells are used to determine rate and phase composition for fluid produced from or injected into the wells on a near real time and continuous basis. The method can be used to alert operators or others of conditions in a well or field to permit more immediate response. The method can be more tolerant of sensor faults, sensor drift, anomalous data or other occurrences which could otherwise lead to incorrect output. More than one model can be used to permit a hierarchy of calculation such that the rate and phase determination is more tolerant of anomalous data. The calculated rate and phase values from one or more wells can be reconciled against facility data.

Abstract: The present invention recited a method and apparatus for measuring a parameter of a fluid within a fluid channel using a tracer element injection by an injection element and detected by a detection element at a location separate from the injection element location. The injection element and detection element of the present invention are in communication with a computational element, wherein the computation element is capable of measuring a parameter of a fluid in a fluid channel based upon time of flight data derived from the injection element and the detection element.

Abstract: Method for predicting lithology and porosity of subsurface rocks from seismic reflection data. The seismic data is inverted to yield elastic properties of the rocks such as the compressional and shear impedances. A rock physics model is built to relate porosity, the shale volume fraction, the fluid content of the rock and the elastic properties of the rock. The model is run backward in a second inversion process to solve for porosity and lithologic properties such as the shale volume fraction.

Abstract: Methods and apparatus for improving the determination of at least one multiphase flow parameter of an earth formation comprising propagating within the formation a first salinity front, determining a first value of the at least one multiphase flow parameter, propagating in the formation a second salinity front and improving the determination of the at least one multiphase flow parameter from the first value for the at least one multiphase flow parameter and a saturation profile associated with the first and second salinity fronts within the formation.

Abstract: This invention relates to computing numerical solutions of linear systems of equations, specifically to implementing preconditioning of the coefficient matrix of such a system. The preconditioning applies to any coefficient matrix, dense or sparse, based on the solutions of a physical problem of unknown functions, commonly referred to as basis or interpolation functions, where the basis function spans more then one mesh element. Examples of such linear systems can result from, as examples, an electromagnetic analysis of printed circuit boards or field scattering in radar applications, fluid mechanics and acoustics. A method and system to compute a preconditioner for a coefficient matrix A that is compatible with the linear system of equations that provides basis function support over at least two mesh elements. Coupling of the preconditioner between partitions of a portioned mesh representation is only through basis functions at the partition boundaries.

Abstract: To determine fluid flow information associated with a well, data from plural discrete temperature sensors in zones without radial fluid flow in the well are received. The fluid flow information associated with the well is derived by combining the received data from the plural discrete temperature sensors with reservoir data that is independent of the fluid flow. That information can be further enhanced by monitoring the changes of temperature data over time.

Abstract: A method and computer program product of either preventing or minimizing pore pressure increase near the wellbore wall within argillaceous formations through which a borehole has been drilled. By interpreting relevant drilling experience data, the type, extent and time-dependency of wellbore instability mechanisms are determined. The impact of drilling mud designs on the time-dependent wellbore instability and hole enlargement is determined by back-analyzing observed drilling events. At least one field-based criterion relationship between net mud weight reduction percentage ratio and hole enlargement is determined. A maximum allowable percentage ratio(s) of net mud weight reduction and either breakout mud weight or mud weight used for the adopted maximum hole enlargement that the wellbore may experience during drilling is determined. Drilling mud salinity value and salt type to satisfy maximum allowable percentage ratio(s) is then determined.

Abstract: Dynamic gain setting is provided for a multiphase flow measuring device that measures in-phase and quadrature signals. The average current at a plurality of measure electrodes is determined and subtracted from the individual measure currents to give a measurement with improved dynamic range.

Abstract: An apparatus for automated fluid rheology measurement includes a probe, a torque tube engaging an inner cylinder of the probe, a motor that drives the inner cylinder, and a circuit board programmed to drive the motor at variable shear rates in a time-dependent manner. A method of continual automated fluid rheology measurement uses an automated viscometer for measuring fluid rheology at differing shear rates. The method involves (1) obtaining and storing data at the differing shear rates, (2) calculating fluid rheology values, and repeating these steps until a user-determined endpoint.

Abstract: Techniques for evaluating physical aspects of a formation fluid from within a wellbore include changing at least one of a pressure on and a temperature of a sample of the formation fluid and transmitting at least one acoustic pulse through the fluid sample and analyzing acoustic information collected. Apparatus and methods for the evaluating involve using at least one acoustic transducer. Analyzing typically involves use of formulae that relate equation (s) of state and other properties for the fluid to a change in the sound speed in the fluid.

Abstract: A method for measuring the wettability of rocks by low-field nuclear magnetic resonance is described which has application for hydrocarbon reservoir engineering or development, or in civil engineering. The method essentially determines the water wet pore surface and the oil wet pore surface when the sample is saturated with water and oil, by measurements of relaxation times (T1, T2) of the sample placed in a nuclear magnetic resonance device, previously brought to various water or oil saturation states, and calculating the wettability index by combination of values the relaxation times obtained for said surfaces.

Abstract: The present invention includes a method and system for automated choke control on a hydrocarbon producing well. The method includes sensing parameters such as sand production, well head tubing pressure, flow rate through the well system, and pressure differential between downhole tubing and an adjacent producing formation. The opening size of the choke is automatically changed, without need for human intervention, in response to changes in the sensed parameters to optimize operation of the well system. Ideally, a programmable logic controller is used to determine when the opening size on the choke should be changed. An advisory communication is preferably sent to a remotely located monitoring station advising that the opening size of the choke has been automatically adjusted. A decision is then made whether an operator should go to the well site and manually well test and shake out the well before placing the system on automatic operation again.

Abstract: The invention allows construction of the architecture of sedimentary layers deposited on a topographic surface as a result of the passage of sedimentary currents, using cellular automata. The topographic surface is discretized into a set of cellular automata. Then, for each current and each automaton, a stationary physical state resulting from an iterative process wherein intermediate physical states are constructed by material and energy exchanges between the automata is determined. These automata are thus constructed so as to be representative of the action of the sedimentary currents on the construction of the sedimentary architecture. A set of geologic variables allowing the architecture of the sedimentary layers to be reconstructed is thus deduced from these stationary physical states. Preferred fields of application of the method are oil exploration and production and reservoir characterization.

Abstract: A method, practiced by a Hydraulic Fracturing Simulator software, is adapted for simulating a hydraulic fracture in an Earth formation, the formation including a plurality of layers, a slip zone, which represents a debonding between adjacent layers of the Earth formation, possibly existing between the adjacent layers, the method comprising the step of: setting up an Influence Coefficient Matrix, the matrix having a plurality of numbers, the plurality of numbers of the Influence Coefficient Matrix being calculated and determined such that the method, practiced by the Hydraulic Fracturing Simulator software for simulating the hydraulic fracture, will model, or take into account an existence of, the slip zone and/or the debonding between the adjacent layers of the Earth formation.

Abstract: The present invention is a method of generating a model of a random field which has directionally varying continuity. First, a tentative model for the random field is specified. Second, connected strings of nodes within the model are identified. Third, a spectral simulation on each of the strings of nodes is performed to determine updated values for the random field. Finally, the tentative model is updated with the data values from the spectral simulations.

Abstract: One aspect of the invention relates to a method (10) for characterising a well (700) using distributed temperature sensor (DTS) data to optimise a well model (12). The method comprises providing a well model of flow and thermal properties (12a, 12b) of the well (700), where the well model (12) has a plurality of adjustable physical parameters, providing a data set made up of a plurality of DTS temperature profiles of the well (700) taken at different times during operation of the well (700), and running the well model (12) with different combinations of the plurality of adjustable physical parameters to match to the plurality of DTS temperature profiles. The DTS temperature profiles may also be pre-processed to make them consistent with one another.