Abstract: A method is disclosed for using a three-dimensional seismic image of a subsurface earth volume to construct a geologic model specifying the spatially-varying grain size distribution, porosity, and permeability throughout the volume. The method applies to earth volumes composed of water-lain clastic sedimentary deposits and involves, in one embodiment, (a) identifying the outline forms of geologic bodies in geologic data; (b) using the outline forms of the geologic bodies to determine the spatially-varying grain size distribution within the bodies, guided by assumptions about the nature and behavior of the paleoflow that deposited the bodies; (c) determining rock properties such as, porosity and permeability within the geologic bodies based on grain-size distribution, mineralogy and burial history information.

Abstract: Method for designing a converted mode (PS or SP) seismic survey to accomplish specified vertical and lateral resolution objectives at target depth. An equation (181) is provided for determining the minimum bandwidth required for a desired vertical resolution at a selected scattering angle, as a function of incident and reflected wave velocities, one of which is the P-wave velocity and the other is the S-wave velocity. A second equation (182) is provided for determining migration acceptance angle from the desired vertical and lateral resolutions. Source and receiver apertures may then be determined by ray tracing. Finally, a third equation (183) is provided for the maximum bin size to avoid aliasing, given the migration acceptance angle and a maximum frequency needed to achieve the bandwidth requirement. Source and receiver spacing may then be based on the maximum bin size.

Abstract: Systems and associated methods for use in the production of hydrocarbons are described. The systems include a first tubular member and a second tubular member at least partially enclosing the first tubular member. Each of the first and second tubular members have a non-permeable longitudinal section and a permeable longitudinal section. The non-permeable longitudinal section of the second tubular member is disposed adjacent to the permeable longitudinal section of the first tubular member. Similarly, the permeable longitudinal section of the second tubular member is disposed adjacent to the non-permeable longitudinal section of the first tubular member. The permeable longitudinal section of the second tubular member is separated from the permeable longitudinal section of the first tubular member by a specific longitudinal distance.

Abstract: Method for organizing computer operations on a system of parallel processors to invert electromagnetic field data (11) from a controlled-source electromagnetic survey of a subsurface region to estimate resistivity structure (12) within the subsurface region. Each data processor in a bank of processors simultaneously solves Maxwell's equations (13) for its assigned geometrical subset of the data volume (14). Other computer banks are simultaneously doing the same thing for data associated with a different source frequency, position or orientation, providing a “fourth dimension” parallelism, where the fourth dimension requires minimal data passing (15). In preferred embodiments, a time limit is set after which all processor calculations are terminated, whether or not convergence has been reached.

Abstract: Method for separating responses of multiple transmitters m a controlled source electromagnetic survey by using mutually orthogonal transmitter waveforms and transforming the combined response to the frequency domain (144) The mutual orthogonality can be based disjoint frequency spectra or on phase encoding of a common waveform element (FIG. 14).

Abstract: The method for correcting the phase of measured electric signals or magnetic signals of field data from a controlled source electromagnetic survey (CSES) by comparing the measured field data corresponding to a selected frequency to the simulated data for various signal source receiver offsets (71) and correcting the phases of the actual data based on the phase difference for a selected range of small signal offsets (76) based on a go-electric model.

Abstract: Described herein are methods of evaluating reservoirs. At least one of the methods includes providing a three dimensional reservoir framework having a plurality of cells; assigning one or more constant reservoir property values to some or all of the cells to provide a first three dimensional reservoir model; updating the first three dimensional reservoir model by populating some or all of the cells with one or more variable reservoir property values to provide a second three dimensional reservoir model; and updating the second three dimensional reservoir model by populating some or all of the cells with one or more reservoir property values derived from seismic data to provide a third three dimensional reservoir model. Other methods are also described.

Abstract: This patent delineates methods for quantifying and mitigating dip-induced azimuthal AVO effects in seismic fracture detection using Azimuthal AVO analysis by accurately accounting for the divergence correction and azimuthal dependence of the reflection angle. Solutions are provided for three cases: (1) dipping isotropic reservoirs; (2) anisotropic reservoirs with fractures aligned in arbitrary direction; and (3) anisotropic reservoirs where vertical fractures are aligned perpendicular to the dip direction.

Abstract: The invention relates to a computer system and method for simulating transport phenomena in a complex system. The computer system comprises a logic interface that enables a user of the computer system to dynamically construct logic to customize simulation of the physical system, a means for converting the constructed logic into corresponding object-oriented code, a means for integrating the object-oriented code with the main simulation system which comprises a simulation data model and simulation algorithms, resulting in an integrated simulation system, and a means for executing the integrated simulation system.

Abstract: A method and apparatus for predicting sand-grain composition and sand texture are disclosed. A first set of system variables associated with sand-grain composition and sand texture is selected (605). A second set of system variables directly or indirectly causally related to the first set of variables is also selected (610). Data for each variable in the second set is estimated or obtained (615). A network with nodes including both sets of variables is formed (625). The network has a directional links connecting interdependent nodes. The directional links honor known causality relationships. A Bayesian network algorithm is used (630) with the data to solve the network for the first set of variables and their associated uncertainties.

Abstract: A method and system for drilling a wellbore is described. The system includes a wellbore with a variable density drilling mud, drilling pipe, a bottom hole assembly disposed in the wellbore and a drilling mud processing unit in fluid communication with the wellbore. The variable density drilling mud has compressible particles and drilling fluid. The bottom hole assembly is coupled to the drilling pipe while the drilling mud processing unit is configured to separate the compressible particles from the variable density drilling mud. The compressible particles in this embodiment may include compressible hollow objects filled with pressurized gas an configured to maintain the mud weight between the fracture pressure gradient and the pore pressure gradient. In addition, the system and method may also manage the use of compressible particles having different characteristics, such as size, during the drilling operations.

Abstract: Method for controlling the phase spectrum of seismic data to match assumptions inherent in subsequent processing steps. The source signature, after processing with the same initial processing steps used on the data, is used to design a phase control filter that shapes the seismic data to have a minimum phase spectrum or whatever other phase spectrum the subsequent processing algorithms may assume. The processed data is then filtered with a second phase-control filter, also designed using the parallel-processed signature, to shape the data to zero phase or whatever other phase may be desired for interpretation of the data.

Abstract: Method for survey design including configuring electrodes to reduce near-surface noise in the seismic response from an electroseismic survey of a subterranean formation. Different embodiments of the invention include (1) selective measurement of the surface noise to remove it from the data; (2) suppressing surface noise generation by reducing electric fields in the vicinity of some of the electrodes; (3) creating source signature differences between the near-surface seismic response and the deep response enabling the near surface response to be removed in data processing; (4) applying an external near-surface magnetic field to modulate the near-surface seismic response, enabling it to be removed in processing; and (5) constructing a partial Faraday cage to shield a near-surface region from fields generated by the electrodes.

Abstract: The present invention is a method of processing seismic data in which one or more seismic vibrators are activated with one or more pilot signals and vibrator motions are recorded along with seismic data. Vibrator signatures are computed from measured vibrator motions, such as the ground force signal. A desired impulse response is specified from either a measured vibrator motion or from test data or field data from a location near the location from which the seismic data was acquired. A deconvolution filter is computed from the impulse response and the vibrator signature. Alternatively, a single separation and deconvolution filter is derived from the impulse response and from vibrator signatures from multiple vibrators and sweeps. The deconvolution or deconvolution and separation filter is used to process the seismic data. The vibrators are then moved to a new location, and the activation is repeated.

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: Method for enhancing resistive anomalies in electromagnetic geophysical survey data. Scaled values of a measured electromagnetic field parameter are plotted on a depth section at locations related to corresponding source/receiver locations. Scaling is performed relative to a reference signal selected to represent a baseline free of unknown resistive bodies. Scaled values are represented by a color scale in the display, and the color scale may be adjusted to enhance perceived anomalies. The method may be employed in either the frequency domain or the time domain.

Abstract: The present invention is a method for constructing a geologic model of a subsurface reservoir which is a composite sedimentary body (30) composed of many smaller fundamental sedimentary bodies (32). In one embodiment, the fundamental sedimentary bodies (32) within the composite body (30) are specified by the properties of the flow which built them, including the flow properties that existed at the local inlet (36) of each fundamental body (32). The statistical distribution of these local inlet properties is characterized throughout the composite sedimentary body (30) using either outline forms of some of the fundamental sedimentary bodies (32) or a well penetration which samples the composite sedimentary body (30).

Abstract: A switch core is set forth that comprises a plurality of duplex switches that are interconnected with a interconnection fabric to implement, for example, strictly non-blocking operation of the switch core for reciprocal traffic. In one embodiment, an N-way reciprocal switch is implemented. The N-way reciprocal switch comprises a plurality of duplex switches numbering N of at least a 1×(N?1) switch type (e.g., the duplex switches have at least N?1 ports available for connection to implement the interconnection fabric). The interconnection fabric interconnects the plurality of duplex switches so that each duplex switch is connected to every other duplex switch used in the interconnection fabric by a single connection. A similar architecture using switches numbering N of at least a 1×N switch type are also set forth.

Abstract: An echo canceller circuit for use in an echo canceller system is set forth that provides sensitive double-talk detection. The echo canceller circuit comprises a second digital filter having adaptive tap coefficients to simulate an echo response occurring during the call. The adaptive tap coefficients of the second digital filter are updated over the duration of the call using a Least Mean Squares process having an adaptive gain a. A channel condition detector is used to detect channel conditions during the call. The channel condition detector is responsive to detected channel conditions for changing the adaptive gain a during the call. For example, the channel condition detector may detect the presence of a double-talk condition and set the adaptive gain a to zero. Similarly, the channel condition detector may detect the occurrence of a high background noise condition and set the adaptive gain a to a level less than 1 that is dependent on the detected level of the background noise.

Abstract: A method to convert seismic traces to petrophysical properties is disclosed. One embodiment of the method comprises (a) deriving a combined log seismic response from at least one petrophysical log from a well and at least one seismic trace substantially near the well, (b) convolving the combined log seismic response filter with each seismic trace to convert the seismic traces to logs of petrophysical properties, (c) outputting the petrophysical properties. Equations to derive the combined log seismic response and to convolve the combined seismic response filter are disclosed. The method may be used to output the petrophysical data as a two-dimensional cross-section or as a three-dimensional volume cube.